BASEMENT Reference Manual: Description of 1D Geometry File Structure and Parameters


Please read the INTRODUCTION first

Input block structure

CROSS_SECTION

INTRODUCTION
In this doc_BASECHAIN_GEOMETRY.html file the structure of the geometry file for BASEchain is described. The topography of a 1D simulation is defined by the repeated definition of cross sections. The topography file itself is the root block called BASEMENT_GEOMETRY. Thus the block BASEMENT_GEOMETRY must not be specified in the topography file.

List of Blocks


BASECHAIN_GEOMETRY
BASECHAIN_GEOMETRY compulsory repeatable
Description:
Inner Blocks: CROSS_SECTION


CROSS_SECTION
CROSS_SECTION compulsory repeatable
Super Blocks: BASECHAIN_GEOMETRY
Description: Every element in the computational grid receives a cross section, which defines the topographical information used to describe the hydraulic and sediment properties within the element. The cross sections are separated by the difference of their 'distance_coord' which is the coordinate along the river.
Each cross section is defined within its own cross section block. The definition of the cross section requires at least the cross section nodes, the name of the cross section as a unique identifier and the distance coordinate.
Node coordinates are defined locally, given by two coordinates: first, the distance to the most left point of the cross section with respect to the flow direction and second, the elevation. Other properties such as specific friction values for specific parts of the cross section, the range of the main channel and other properties can also be defined. These properties are defined by giving the ranges over which they span over the cross section. These ranges can either be defined by giving the actual node coordinates along the cross section (which have to correspond exactly to a given node coordinate) or by referencing the slice directly using slice indexes. A slice refers to the segment between two adjacent nodes.
Concerning Friction: For 1D, only Strickler Values are a valid input. In case you are used to Manning Values, please make the conversion k_str = 1/n_manning on your own. Note that you can define the friction coefficient on several stages. Your first option is to define a default_friction in the BMC file within the FRICTION block. Your second option is to define a default friction in the cross section itself which overwrites the coefficient from the BMC file. Further, it is possible to assign different frictions to different ranges.
Inner Blocks: SOIL_DEF

Variable type comp. condition units default values
name string yes - - - -
distance_coord real yes - [km] - -
sternberg real no - - 0.0 -
reference_height real no - [m] 0.0 -
left_point_global_coords comma separated list of 3 reals no - [m] - -
orientation_angle real no - [deg] 0.0 0 <= x <= 360
node_coords comma separated list of sublists of 2 reals yes - [m] - -
interpolation_fixpoints comma separated list of integers no - - - Node index (starting from 1).
water_flow_range comma separated list of 2 reals no exclusive: 'water_flow_slice_indexes' [m] - exact coordinates of a defined node
water_flow_slice_indexes comma separated list of 2 integers no exclusive: 'water_flow_range' - - Two slice indexes (starting from 1). The indexes supplied are included in the range.
main_channel_range comma separated list of 2 reals no exclusive: 'main_channel_slice_indexes' [m] - exact length coordinates of a defined node
bottom_range comma separated list of 2 reals no exclusive: 'bottom_slice_indexes' [m] - exact length coordinates of a defined node
main_channel_slice_indexes comma separated list of 2 integers no exclusive: 'main_channel_range' - - Two slice indexes (starting from 1). The indexes supplied are included in the range.
bottom_slice_indexes comma separated list of 2 integers no exclusive: 'bottom_range' - - Two slice indexes (starting from 1). The indexes supplied are included in the range.
active_range comma separated list of 2 reals no exclusive: 'active_slice_indexes' [m] - exact length coordinates of a defined node
active_slice_indexes comma separated list of 2 integers no exclusive: 'active_range' - - Two slice indexes (starting from 1). The indexes supplied are included in the range.
bed_form_factor real no - - 1.0 0 < x <= 1
theta_critic real no - - -1.0 -1 <= x <= 1
friction_calibration_factor real no - - 1.0 -
default_friction real no - - - x >= 0
friction_coefficients comma separated list of reals no - m^(1/3)/s - -
friction_ranges comma separated list of sublists of 2 reals no exclusive: 'friction_slice_indexes' [m] - exact length coordinates of a defined node
friction_slice_indexes comma separated list of sublists of 2 integers no exclusive: 'friction_ranges' - - Two slice indexes (starting from 1). The indexes supplied are included in the range.
table_fixed_points comma separated list of reals no - [m] - -
table_fixed_lower_limit real no - [m] - -
interpolated boolean no - - no {no, yes}

Description of Variables

name - The unique name of the cross section.
distance_coord - Distance coordinate of cross section, measured from up to downstream. Note that the units of the coordinate is in [km]!
sternberg - Sternberg-factor for the reduction of volume by abrasion V = V0e**(s*Dx)
reference_height - The value of the reference height is added during input to the node elevation.
left_point_global_coords - Global coordinates of the left cross section point, which is used as the geographical reference point. If this value is provided together with the 'orientation_angle' for all cross sections, the modeled river is displayed in the region view of the graphical user interface. Further, this information is used by the cross section interpolator to estimate curvatures etc.
orientation_angle - In-plane angle of the cross section line with left point as origin. The angle is measured in a mathematical way, counter-clockwise: The angle 0 is given by the (1,0) vector.
node_coords - Coordinates of the points defining the cross section. The first coordinate is the distance to the left point, the second the absolute elevation.
interpolation_fixpoints - The interpolation fixpoints are used to create new cross sections in between existing ones. Fix points are in interpolations connected to each other. This allows to cleanly transfer slopes etc. into the interpolated cross section. The first and last nodes of the cross sections are for obvious reasons always included. Counting starts at 1.
water_flow_range - Defines which part of the cross section is considered to be flowing. Only this part will be used for solving of the momentum equation, rest of the cross section is considered to be only storage area in the continuity equation. The default is that the whole cross section area is flowing. Try to not define storage areas in the main channel, as in this case a stable behaviour is not assured, i.e. when only the non flowing part of the channel is wetted! Also the non flowing areas are not taken in accont for sediment deposition and erosion!
water_flow_slice_indexes - Same as tag 'water_flow_range', except the fact that not coordinates are supplied but slice indices.
main_channel_range - Range defining the span of the main channel. Everything left and right of the channel is treated as floodplain. If not defined, the whole cross section is assumed to belong to the main channel.
bottom_range - Range defining the span of the bottom/bedload activ zone.
main_channel_slice_indexes - Same as 'main_channel_range', except the fact that not coordinates are supplied but slice indices.
bottom_slice_indexes - Same as 'bottom_range', except the fact that not coordinates are supplied but slice indices.
active_range - This range defines the active part of the cross section and must span from the left to the right dike. If not defined, the whole cross section is assumed to be active and the elevation of the first and last point will be assumed to correspond to the dike elevation.
active_slice_indexes - Same as 'active_slice_indexes', except the fact that not coordinates are supplied but slice indices.
bed_form_factor - Factor considering the influence of bed forms on bottom shear stress (ripple factor). This factor reduces the bottom shear stress due to bedforms with a constant value. If this factor is set to 1.0 (default), there is no influence and the standard shear stress is used.
theta_critic - User defined dimensionless critical bottom shear stress (Shields paramter) for this cross section. If specified by the user (values between 0.0 and 1.0), this local theta critic is used for this cross section.However, for negative values of this local theta critical, the standard theta critical is used defined in the BEDLOAD-PARAMETER block (default behaviour)
friction_calibration_factor - Calibration factor between 0.0 and 10.0 to modify the friction value locally for this cross section.
default_friction - Default Strickler friction value which will be used for the cross section where no individual friction coefficient is defined.
friction_coefficients - Strickler Friction values which should be assigned to the 'friction_ranges' (or 'friction_slice_indices').
friction_ranges - Ranges to which the 'friction_coefficients' are applied.
friction_slice_indexes - Same as 'friction_ranges', except the fact that not coordinates but slice indices are supplied.
table_fixed_points - Some fixed points which should be included in the table calculation. Note that these values are only used if you use lookup tables for hydraulic properties (see SECTION_COMPUTATION within the command file).
table_fixed_lower_limit - This variable allows to fix the lowest water surface elevation used in the 'lookup-table' for A(z) to a different value than the lowest elevation in the cross section. This option is useful if you include a deep lake in your model containing a large water volume not participating in the dynamics of your river system.
interpolated - This is an internal tag used to indicate that the cross section is the result of an interpolation of two primary cross sections.

SOIL_DEF
SOIL_DEF compulsory repeatable
Super Blocks: BASECHAIN_GEOMETRY > CROSS_SECTION
Description: This block is used to assign a soil to a part of the cross section. Soil assignments must not overlap and must reference a valid soil definition, defined in the input bmc file.

Variable type comp. condition units default values
index integer yes - - - x >= 0
range comma separated list of 2 reals yes exclusive: 'slice_indexes' [m] - -
slice_indexes comma separated list of 2 integers yes exclusive: 'range' - - -

Description of Variables

index - Index of the soil, as defined within the SOIL_ASSIGNMENT block in the bmc input file.
range - Range in cross section to which the soil with the corresponding 'index' should be assigned.
slice_indexes - Same as 'range', but refers to slice index rather than coordiantes.