Difference between revisions of "Module DrainageNetwork"

Revision as of 16:36, 15 February 2010

Overview

Main Processes

Other Features

Creatig a drainage network file

• Obtain a drainage network with MOHID GIS in Delineate Basins

• Alternatively can also obtain a drainage network with program BasinDelimiter (from SourceSafe) with a basin.dat with the keywords, adapted to the study watershed:

TOPOGRAPHIC_FILE              : ..\..\GeneralData\DTM\MDT200mSD.dat
TRESHOLD_AREA                 : 100000
DELINEATE_BASIN               : 1
OUTLET_I                      : 1
OUTLET_J                      : 44
WRITE_REACHES                 : 1
REACHES_FILE                  : ..\..\GeneralData\DrainageNetwork\DrainageNetwork.dnt

• Define the cross sections for the drainage network file in MOHID GIS.

DO NOT FORGET TO SELECT SAVE ALL IN MOHID GIS TO UPDATE THE FILE!!!

(Currently in this version of the cross sections program keyword TERRAIN_LEVEL is not written, instead is written BOTTOM_LEVEL (old version). TERRAIN LEVEL = BOTTOM_LEVEL + HEIGHT. Change it for example with excel, new program should be written)

• Define the network file just created in the drainage network data file with:

NETWORK_FILE : ..\..\GeneralData\DrainageNetwork\DrainageNetwork.dnt

Outputs

References

Data File

DrainageNetwork_X.dat is the input file to this module, where X is the simulation number. The following tables describe the keywords that can be used, their data type, and the default values (in case of omission)

General Keywords

Keyword                   : Data Type         Default     !Comment
NETWORK_FILE              : char              -           !Path to drainage network file
CHECK_NODES               : 0/1               [1]         !Ckeck nodes consistency
CHECK_REACHES             : 0/1               [1]         !Check reaches consistency
GLOBAL_MANNING            : real              -           !Rugosity in Channels
GEO_CONVERSATION_FACTOR   : real              [1.]        !Lat to Meters rough estimation
OUTPUT_TIME               : int int...        [-]

Time Series Keywords

The file specified in TIME_SERIE_LOCATION has to have the numbers of the nodes to output time series between blocks <BeginNodeTimeSerie> / <EndNodeTimeSerie>.

TIME_SERIE_LOCATION       : char              -           !Path to time serie especification nodes
MAX_BUFFER_SIZE           : 1000
COMPUTE_RESIDUAL          : 1
DT_OUTPUT_TIME            : 1200
TIME_SERIE_BY_NODES       : 0/1               [0]         !Keyword to see if the user wants the time series to be written by nodes, i.e.,
                                                          !One file per node, with all variables in the headers list
                                                          !if FALSE, its one file per variable with nodes in the headers.

Stabilization Keyword

STABILIZE                 : 0/1               [0]         !Restart time iteration if high volume gradients
STABILIZE_FACTOR          : real              [0.1]       !max gradient in time steps as fraction of old volume
MAX_ITERATIONS            : int               [100]       !Max iterations for stabilized check
DT_FACTOR                 : real              [0.8]       !Factor for DT Prediction
MAX_DT_FLOOD              : real              [10.0]      !Max DT if channel water level exceeds full bank

Hydrodynamic Keywords

HYDRODYNAMIC_APROX        : int               [1]         !1 - KinematicWave, 2 - DiffusionWave, 3 - DynamicWave
NUMERICAL_SCHEME          : int               [0]         !0 - ExplicitScheme, 1 - ImplicitScheme
MASS_ERR                  : real(8)           [0.001]     !Max error in mass conservation
MIN_WATER_DEPTH           : real              [0.001]     !Min water depth in nodes (For h < MIN_WATER_DEPTH water stops flowing)
INITIAL_WATER_DEPTH       : real              [0.0]       !Initial water depth
MINIMUM_SLOPE             : real              [0.0]       !Minimum Slope for Kinematic Wave

Processes Keywords

MIN_WATER_DEPTH_PROCESS   : real              [0.01]      !Water Quality Process / Surface Fluxes shutdown
DISCHARGES                : 0/1               [0]         !Use module discharges (WWTP, etc)
TRANSMISSION_LOSSES       : 0/1               [0]         !If user wants to use transmission losses
HYDRAULIC_CONDUCTIVITY    : real              -           !Hydraulic Conductivity to calculate transmission losses
REMOVE_OVERTOP            : 0/1               [0]         !Removes Water if channels are overtoped
AERATION_METHOD           : int               [-]         !1 - PoolAndRifle, 2 - ChannelControled_
T90_DECAY_MODEL           : 0                 [1]         !0 - Constant, 1 - Canteras, 2 - Chapra
T90                       : real              [7200.]     !if T90_DECAY_MODEL = Constant
SHADING_FACTOR            : real              [1.]        !0-1 fraction of riparian shading
FRACTION_SEDIMENT         : 0/1               [0]
GLOBAL_TOXICITY           : char              ['SUM']     !Global Toxicity Computation Method : SUM,MAX,RISKRATIO

Downstream Boundary Keywords

DOWNSTREAM_BOUNDARY       : int               [1]         !0 - Dam, 1 - ZDG, 2 - CD, 3 - ImposedWaterDepth, 3 - ImposedVelocity
DEFAULTVALUE              : real              -           !Default value for water depth at downstream boundary
FILE_IN_TIME              : char              [NONE]      !If DOWNSTREAM_BOUNDARY = ImposedWaterDepth, this can be NONE or TIMESERIE
FILENAME                  : char              -           !If FILE_IN_TIME = TIMESERIE, this is the name of timeserie file for the downstream boundary
DATA_COLUMN               : int               -           !Number of column with data in FILE_IN_TIME

<beginproperty>
  NAME                    : cohesive sediment
  UNITS                   : mg/L
  DESCRIPTION             : cohesive sediment
  INITIAL_VALUE           : 100.00
  MIN_VALUE               : 0.0
  ADVECTION_DIFUSION      : 1
      ADVECTION_SCHEME    : 1       !Upwind
      DIFFUSION_SCHEME    : 5       !CentralDif
      DIFFUSIVITY         : 1E-8    !m2/s
      VIRTUAL_COEF        : 0.01
  WATER_QUALITY           : 0
  BENTHOS                 : 0
  Decay                   : 0       !uses T90 decay model for fecal coliform
  TIME_SERIE              : 1
<endproperty>

Network file:

<BeginNode>
  ID                      : int               -           !Node ID number
  COORDINATES             : real real         -           !Node coordinates
  GRID_I                  : int               -           !I position of node, if grid
  GRID_J                  : int               -           !J position of node, if grid
  TERRAIN_LEVEL            : real              -           !Bottom level of cross section
  MANNING_CHANNEL         : real         GLOBAL_MANNING   !Node rugosity
  WATER_DEPTH             : real      INITIAL_WATER_DEPTH !Node initial water depth
  CROSS_SECTION_TYPE      : int               [1]         !1 - Trapezoidal, 2 - TrapezoidalFlood, 3 - Tabular
  1 - Trapezoidal, 2 - TrapezoidalFlood
      BOTTOM_WIDTH        : real              -           !Bottom width of cross section
      TOP_WIDTH           : real              -           !Top width of cross section
      HEIGHT              : real              -           !Max height of cross section
  2 - TrapezoidalFlood
      MIDDLE_WIDTH        : real              -           !Middle width of cross section
      MIDDLE_HEIGHT       : real              -           !Middle height of cross section
  3 - Tabular
      N_STATIONS          : integer           -           !number os stations that define the cross section
      STATION             : real real ...     -           !station values
      ELEVATION/LEVEL     : real real ...     -           !elevation values
<EndNode>
<BeginReach>
  ID                      : int               -           !Reach ID Number
  DOWNSTREAM_NODE         : int               -           !Downstream node ID
  UPSTREAM_NODE           : int               -           !Upstream node ID
<EndReach>

Sample

NETWORK_FILE              : ..\..\GeneralData\DrainageNetwork\DrainageNetworkAlcabrichel_100ha_terrain.dnt
CHECK_NODES               : 1
CHECK_REACHES             : 1
HYDRODYNAMIC_APROX        : 1
MASS_ERR                  : 0.001
GLOBAL_MANNING            : 0.06
MIN_WATER_DEPTH           : 0.01
INITIAL_WATER_DEPTH       : 0.0
XS_CALC                   : 1
DOWNSTREAM_BOUNDARY       : 1
TIME_SERIE_LOCATION       : ..\..\GeneralData\Timeseries\TimeSerieLocation_200_new.dat

OUTPUT_TIME               : 0 864000

DISCHARGES                : 1 (1-if have a discharge file; 0 if not discharging)
CONTINUOUS                : 0

MAX_DT_FLOOD              : 100
STABILIZE                 : 1
DT_FACTOR                 : 0.75
STABILIZE_FACTOR          : 0.05

MINIMUM_SLOPE             : 0.0001

If want to model Properties

Than additionaly to the above options need to define properties to be modelled. One block for each property (see a list

<beginproperty>
NAME                    : salinity
UNITS                   : psu
DESCRIPTION             : salinity
DEFAULT_VALUE           : 0.3
MIN_VALUE               : 0.0
ADVECTION_DIFUSION      : 1       !1-want to transport property; 0 - no transport
DISCHARGES              : 0       !1-property is discharged;0 - no discharges
TIME_SERIE              : 1
WATER_QUALITY           : 0       !1-use water quality model for property transformation;0-no transformation
<endproperty>