Difference between revisions of "Template:GUM17"

From Gsshawiki
Jump to: navigation, search
 
(No difference)

Latest revision as of 18:51, 20 May 2008


17.1 Project File Format

The project file is the heart of a GSSHA simulation. It controls what processes run, when they run, and what output is returned. The project file references many other files. These files often differ in name by not much more than the extension. The extensions shown here are the defaults used in WMS, but they are not fixed. Filenames or paths, if they contain spaces, should be enclosed in double quotes.

Note: This section has not been kept up to date with developments. Please refer to Chapter 3 of the Manual, which describes the project file, and projects cards in detail.

Comments in the project file begin the line with a #.

Header

A GSSHA project file always begins with the same header

Card Input Req'd Comments
GSSHAPROJECT none Y Identifies the file as a GSSHA project file.


WMS info

Next comes info that WMS needs when reading the project

Card Input Req'd Comments
WMS version Y This simply conveys the version of WMS used to create the project.
FLINE filename.map Y The WMS map file associated with this project. WMS reads this file in to create the GIS data.


General Simulation Data

The general simulation data provides the framework for how the simulation is to be set up and run. The mask file controls what cells are on and off, and all other grid files must have the exact same dimensions as the mask file.

Card Input Req'd Comments
WATERSHED_MASK filename.msk Y The watershed mask controls the simulation domain.
PROJECT_PATH path N Path to all project files. This path is prepended to all files that GSSHA tries to open.
METRIC none N output hydrographs are in cms (default)
ENGLISH none N output hydrographs are in cfs,either specify METRIC or ENGLISH, not both
GRIDSIZE value Y value in meters of the grid dimensions
ROWS value Y number of rows in the simulation, should match the value in the watershed mask
COLS value Y number of colums in the simulation, should match the value in the watershed mask
TOT_TIME value Y simulation length (non-long-term) in minutes
TIMESTEP value Y overland and stream time step in seconds
OUTROW value Y outlet row number, for overland flow only
OUTCOL value Y output column number, for overland flow only
OUTSLOPE value Y outlet slope, for overland only
TIME_SERIES_FILE filename.ts N read file containing time series; this is optional and it may be repeated to include several files
EVENT_MIN_Q value Y the minimum flow to signal the end of event
END_TIME year month day hour min N if this time is reached in the simulation then the simulation automatically ends.
NUM_THREADS ## N If using GSSHA compiled with OpenMP, then this will change the default thread level (which is to use all logical processors).

The simulation will end at the end of hmet data when the flow has fallen below event_min_q (for long-term runs) or when the tot_time has been reached (uniform or gaged precip, non-long-term runs) or if the current simulation time is after the end time.

Replacement cards

GSSHA has the ability to replace any part of any line of any input file with a value from another file. These files must be set up for batch mode processing.

Card Input Req'd for Replacement Comments
REPLACE_PARAMS filename.ext Y Replacement parameter file.
REPLACE_VALS filename.ext Y Replacement value file.
REPLACE_FOLDER path N Puts all output files in the specified folder.
REPLACE_LINE ## N Use a certain line from the values file.

Output cards

These cards control the output of GSSHA. There are several output formats-- gridded data, link/node data, time series data, and summary files.

Card Input Req'd Comments
SUMMARY filename.sum Y Summary output file. Most important output file of them all.
OUTLET_HYDRO filename.otl Y outlet hydrograph time series
MAP_FREQ value Y gridded output frequency, in minutes
HYD_FREQ value Y hydrograph, link/node dataset output frequency, in minutes
MAP_TYPE value Y gridded output format, 0=ascii grid sequence, 1=WMS gridded dataset
DEPTH filename.dep or path N overland flow depth maps, if map_type==0 then put a path else put a filename.
SURF_MOIST filename.smo or path N soil moisture output, see DEPTH
GW_OUTPUT filename.gwh or path N groundwater heads, see DEPTH
DIS_RAIN filename.drn or path N rainfall rate as it is being applied, see DEPTH
INF_DEPTH filename.idp or path N infiltration depth, see DEPTH
DISCHARGE filename.dsc or path N overland discharge, see DEPTH
CHAN_DEPTH filename.cdp N channel flow depth, link/node format
CHAN_DISCHARGE filename.cdq N channel discharge, link/node format
CHAN_STAGE filename.cst N channel stage, link/node format
CHAN_VELOCITY filename.cvl N channel velocity, link/node format
STRICT_JULIAN_DATE none N force hydrograph output to be in strict julian date
QUIET none N do not output to screen as model is running
SUPER_QUIET none N do not output to screen
IN_HYD_LOCATION filename.ihl N describes hydrograph output locations
OUT_HYD_LOCATION filename.ohl Y if .ihl output hydrographs from .ihl file
IN_THETA_LOCATION filename.ith N describes soil moisture output locations
OUT_THETA_LOCATION filename.oth Y if .ith output soils moisture time series from .ith
FLOOD_GRID filename.fgd N outputs the maximum water depth on the overland grid
FLOOD_STREAM filename.fst N outputs the maximum water depth on the stream network
OPTIMIZE filename.opt N outputs the values of event peak flow and event volume for use by external optimization methods

Overland Flow

There are several cards to set up overland flow. The elevation file must be specified, an overland solver selected, and roughness values set up.

Card Input Req'd Comments
ELEVATION filename.ele Y surface elevation grid file
OVERTYPE EXPLICIT, ADE, or ADE-PC Y overland flow solver (EXPLICIT is being deprecated in v8.0 and will not be an option in future versions)
OVERBANK_FLOW none N allow streams to flow back to overland
MAX_COURANT_NUMBER value N used to change the maximum courant number allowed from the default (0.06)
INITIAL_DEPTH filename.idp N read in a hot start file
READ_OV_HOTSTART filename.idp N read in a hot start file, same as INITIAL_DEPTH
WRITE_OV_HOTSTART filename.odp N write the overland depths at the end of the simulation
ROUGHNESS filename.ovn N GSSHA reads a gridded file of n values; if not present and MANNING_N not present, then GSSHA assumes to use the mapping table. If the mapping table is not properly set up you will see an error.
MANNING_N value N use uniform manning's n value for overland flow, see ROUGHNESS
OV_BOUNDARY none N do overland boundaries according to the OVERLAND_BOUNDARY table in the mapping table file.
EMBANKMENT filename.dik N set up no-flow embankments between overland grid cells.
LOWSPOT_FILE filename.lsp N set up low spots for overtopping on the embankments

Stream Flow

There are several input files that set up the stream network. The most important of these are the channel input file (.cif) and the grid/stream file (.gst.) The stream network is not required to run a GSSHA simulation, so the required column lists those cards that are required to run the stream network.

Card Input Req'd for Streams Comments
NON_ORTHO_CHANNELS none Y Run the channel model
DIFFUSIVE_FLOW none Y Use the diffusive flow numerical model, currently the only supported numerical model for streams.
STREAM_CELL filename.gst Y The grid-stream file tells GSSHA how the streams link to the grids.
CHANNEL_INPUT filename.cif Y The channel input file describes the entire stream network, e.g. cross-sections, thalweg profile, roughness values, etc.
HEAD_BOUND none N Use a head boundary at the channel outlet
BOUND_DEPTH value if HEAD_BOUND then either BOUND_DEPTH or BOUND_TS Use a constant value (m) for the outlet channel depth (not stage.)
BOUND_TS time series name if HEAD_BOUND then either BOUND_DEPTH or BOUND_TS Use a time series to specify the outlet boundary depth (m) (not stage.)
WRITE_CHAN_HOTSTART filename.ext N Create two channel hot-start files at the end of the simulation. Depth file has extension .dht, flow has extention .qht.
READ_CHAN_HOTSTART filename.ext.ext N Read in a hot-start file, created by WRITE_CHAN_HOTSTART.
HOTSTART_IS_DEPTH none N Treat hot-start file as depth of flow when reading it in.
STREAM_LOSS none N Do stream losses even if groundwater is not turned on
M_RIVER value N Uniform value (m) for stream sediment thickness for stream losses. Only use if not specified in the channel input file.
K_RIVER value N Uniform hydraulic conductivity (cm/hr) of streambed for stream losses. Only use if not specfied in the channel input file.
EXPLIC_BACKWATER N
FIRST_ORDER_CHAN_Q N


Infiltration

These are the four different infiltration options. One of these are required to run infiltration.

Card Input Req'd for Groundwater Comments
GREEN_AMPT none Run single or two layer Green & Ampt. (Single event model)
INF_REDIST none Use G&A w/ soil moisture redistribution (Long-term model, may have 1 or 2 layers)
INF_LAYERED_SOIL none Use G&A layered model. (Single event model)
INF_RICHARDS none Use Richard's Equation model. (Long-term model)


These parameters are required for the Richard's Equation infiltration model.

Card Input Comments
RICHARDS_WEIGHT ## Usually 0.8 to 1.2
RICHARDS_K_OPTION geometric or arithmetic
RICHARDS_C_OPTION havercamp or brooks
RICHARDS_UPPER_OPTION green_ampt or average
RICHARDS_ITER_MAX
RICHARDS_DTHETA_MAX
AQUIFER_DELTA_Z ##.##

These are used for the G&A and the G&A w/ Redistribution to create the soil moisture layers. If top_layer_depth is not included, then only one layer will be used.

Card Input Comments
SOIL_MOIST_DEPTH ##.## Full depth for soil moisture computations.
TOP_LAYER_DEPTH ##.## Upper layer depth (included in soil_moist_depth) for just an upper layer soil moisture layer.

There are two cards for hot-starting.

Card Input Comments
WRITE_SM_HOTSTART <filename> Write out the soil moisture at the end of the simulation as a hot start file.
READ_SM_HOTSTART <filename> Read initial soil moistures.

Groundwater

Card Input Req'd for Groundwater Comments
GW_SIMULATION none Y Run groundwater simulation.
WATER_TABLE filename.wte Y Specifies the initial water table elevation.
AUIFER_BOTTOM filename.aqe Y Specifies the lower extent of the aquifer.
GW_ASSIGN_THETA
GW_LSOR_DIR
GW_UNIF_POROSITY
GW_UNIF_HYCOND
GW_TIMESTEP
GW_LSOR_CON
SINGLE_UNSAT_SAT
GW_LEAKAGE_RATE
GW_RELAX_COEFF
GW_FLUX_BOUNDTABLE
GW_BOUNDFILE
GW_OUTPUT
GW_WELL_LEVEL
OUT_WELL_LOCATION


Precipitation

Card Input Req'd Comments
PRECIP_FILE filename.gag either PRECIP_FILE or PRECIP_UNIF input gaged or radar precip data
PRECIP_UNIF none either PRECIP_FILE or PRECIP_UNIF use uniform precipitation
RAIN_THIESSEN none either RAIN_THIESSEN or RAIN_INV_DISTANCE use theissen polygons to distribute rainfall
RAIN_INV_DISTANCE none either RAIN_THIESSEN or RAIN_INV_DISTANCE use inverse distance weighting to distribute rainfall
RAIN_INTENSITY value Y if PRECIP_UNIF Set the precipitation intensity for uniform precip. mm/hr.
RAIN_DURATION value Y if PRECIP_UNIF Set the precipitation duration for uniform precip. Minutes.
START_DATE year month day Y if PRECIP_UNIF starting simulation time for uniform precip
START_TIME hour minute Y if PRECIP_UNIF starting time for uniform precip

Constituent Transport

This section controls all constituents, simple and NSM. To include NSM controls, see the Nutrients section.

Card Input Req'd for Transport Comments
OV_CON_TRANS none N Turns on routing in the overland. Can route in overland w/out streams, soils.
CHAN_CON_TRANS none N Turns on routing in the streams. Can route in streams w/out routing in overland.
SOIL_CON_TRANS none N Turns on routing in the soils. Must have overland routing turned on.
SEDIMENT_CONTAM none N Turns on sorption/desorption of constituents to sediments.
SOIL_STATIC_CONC none N Keeps the soil concentrations constant.
OUT_MASS_LOCATION filename.ext N Filename for output mass fluxes, uses .ihl locations
OUT_CON_LOCATION filename.ext N Filename for output concentrations, uses .ihl locations

Nutrients

Setting up nutrients requires several cards. The transport options need to be specified for whichever flow options are turned on. The aquatic environment file, with its stream index map and overland index map, needs to be set up with the kinetic rates.

General Setup

Card Input Req'd for Nutrients Comments
NUTRIENTS none Y Tells GSSHA to do NSM nutrients kinetics and transport
MIXING_LAYER_DEPTH value Y if SOIL_CONTAM The depth of the surface soil layer containing the constituents
OV_CON_TRANS none Y Tells GSSHA to do overland constituent transport
CHAN_CON_TRANS none Y if streams Tells GSSHA to do channel constituent transport
SOIL_CONTAM none Y if infiltration Tells GSSHA to do vadose zone nutrient transport
AQ_ENV filename.aqe Y contains the controlling kinetic rates
CHAN_PH value Y if computing Kjeldahl N The pH of the water (constant)
WATER_TEMP value Y The temperature of the water (constant)
SOURCEFILE filename.src N describes point and non-point loadings
OV_DISPERSION none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values.


Soils Initialization

Card Input Req'd for Nutrients Comments
SOILS_NSM_P none N Initialize cell loadings from mapping table values
SOILS_NSM_N none N Initialize cell loadings from mapping table values
SOILS_NSM_C none N Initialize cell loadings from mapping table values


Overland Initialization

Values of constituents in water produced by hotstart or init_cell_elev cards.

Card Input Req'd for Nutrients Comments
OV_INIT_NO2 none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_NO3 none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_NH4 none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_ON none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_OP none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_DP none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_ALG none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_CBOD none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_INIT_DO none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_DOC_INIT none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.
OV_DOC_FOC none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in mg/l.

Stream Initialization

Card Input Req'd for Nutrients Comments
ST_MAPPING_TABLE filename.smt N Mapping table for streams.
ST_INIT_NO2 none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_NO3 none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_NH4 none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_ON none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_OP none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_DP none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_PO4 none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_ALG none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_CBOD none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.
ST_INIT_DO none, uniform value, or filename.ext N If blank, use mapping table. If filename, apply grid of values. Value is in kg per cell.


Overland Nutrient Output

Card Input Req'd for Nutrients Comments
total_ov_n filename.ton or path N gridded output of total N mass and concentration (will add .mass and .conc)
total_ov_p filename.top or path N gridded output of total P mass and concentration (will add .mass and .conc)
out_ov_no2 filename.on2 or path N gridded output of overland NO2 mass and concentration (will add .mass and .conc)
out_ov_no3 filename.on3 or path N gridded output of overland NO3 mass and concentration (will add .mass and .conc)
out_ov_nh4 filename.on4 or path N gridded output of overland NH4 mass and concentration (will add .mass and .conc)
out_ov_on filename.oon or path N gridded output of overland ON mass and concentration (will add .mass and .conc)
out_ov_op filename.oop or path N gridded output of overland OP mass and concentration (will add .mass and .conc)
out_ov_dp filename.odp or path N gridded output of overland DP mass and concentration (will add .mass and .conc)
out_ov_po4 filename.op4 or path N gridded output of overland PO4 mass and concentration (will add .mass and .conc)
out_ov_alg filename.oal or path N gridded output of overland Algae mass and concentration (will add .mass and .conc)
out_ov_cbod filename.ocb or path N gridded output of overland CBOD mass and concentration (will add .mass and .conc)
out_ov_do filename.odo or path N gridded output of overland DO mass and concentration (will add .mass and .conc)


Stream Nutrient Output

Card Input Req'd for Nutrients Comments
total_chan_n
total_chan_kn
total_chan_p
out_st_no2
out_st_no3
out_st_nh4
out_st_on
out_st_op
out_st_dp
out_st_po4
out_st_alg
out_st_cbod
out_st_do


Nutrient Loadographs

Card Input Req'd for Nutrients Comments
out_hyd_no2
out_hyd_no3
out_hyd_nh4
out_hyd_on
out_hyd_op
out_hyd_dp
out_hyd_po4
out_hyd_alg
out_hyd_alg
out_hyd_cbod
out_hyd_do


Sediment Transport

Card Input Req'd for Sediments Comments
SOIL_EROSION Y
SOIL_ERODABILITY N Uniform erodability factor
CHAN_SED_FLUX N
IN_SED_LOC N
VOL_SED_SUSP N
MAX_SED_FLUX N
NET_SED_VOLUME N
OUTLET_SED_FLUX
OUT_SED_LOC
SED_POROSITY Y
OUT_SED_LOC

Analyzers

Analyzers look into the GSSHA code as it is running and print out specified data values.

Card Input Req'd for Analyzers Comments
ANALYZERS filename.ext Y Turn on output analyzers.


17.2 Mapping Tables

Mapping tables are the heart of how parameters are set up in GSSHA. Mapping tables work in conjuction with index maps to create a spatially-defined parameter set. The mapping tables refer to an index map, that has an ID value for each active cell in the grid, and a set of paramters for each id for a particular process. The mapping tables have a defined format, but can have as many IDs as needed or desired. The mapping table is divided into two parts, the header and the tables.


Mapping Table Header

The mapping table file must always begin with the line

GSSHA_INDEX_MAP_TABLES

After that line comes a set of index map indentifiers of the following format:

INDEX_MAP "filename.idx" "index map name"
INDEX_MAP "filename.idx" "index map name"
INDEX_MAP "filename.idx" "index map name"
...


These lines associate a particular file with a name. This name is used in the tables that follow. .idx is the usual extension for the index maps, but it can be anything.

General Table Format

The basic format of a mapping table is as follows:

TABLE_NAME "index map name"
NUM_IDS [###]
ID Description ... Parameter descriptions... (text line)
[##] text goes here... [###.###] [###.###] [###.###] ...
[##] text goes here... [###.###] [###.###] [###.###] ...
[##] text goes here... [###.###] [###.###] [###.###] ...
...


TABLE_NAME must be one of the pre-defined tables. "index map name" must refer to an index map (see above.) NUM_IDS must be the count of the number of table entries. The line "ID Description..." must be present but is merely a text line to assist in manually reading and editing the file. The ID can be up to 6 characters, and the description must start at the 7th character. There are two description fields, each of 40 characters. After that begins the parameter value listing.

Standard Format Tables

17.3 Stream Network Format

File Formats:Stream Network Format

17.4 GRASS files

GRASS ASCII grid files are much like Arc/Info ASCII grid files. They are text files that contain a few header lines describing the size and location of the data followed by the data for the grid (row-ordered data with rows on each line, separated by spaces.) The data can either be integer values or floating-point (single or double) values. For example, index maps use integer values while the elevation map uses floating point values.

The first four header lines define the a bounding box for the grid cells, in UTM coordinates. The next two lines define the number of rows and columns there are in the grid. To obtain the cell dimensions in meters, take the difference between the appropriate bounding box coordinates and divide by the number of cells in that direction. For example, the following grid has a y-dimension grid cells size of (4969474.648439-4969114.648439)/12=30.0m. The header lines must have the identifying card and colon.

north: 4969474.648439
south: 4969114.648439
east: 552376.687907
west: 551986.687907
rows: 12
cols: 13
0.000000 372.779999 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 
0.000000 372.380005 372.420013 372.410004 0.000000 0.000000 0.000000 0.000000 371.950012 371.839996 0.000000 0.000000 0.000000 
0.000000 372.359985 372.380005 372.299988 372.320007 372.250000 372.109985 371.940002 371.790009 371.670013 371.589996 0.000000 0.000000 
0.000000 372.209991 372.209991 372.070007 372.119995 372.049988 371.940002 371.730011 371.519989 371.380005 371.350006 0.000000 0.000000 
0.000000 372.040009 372.000000 371.959991 371.880005 371.839996 371.630005 371.359985 371.140015 371.100006 371.029999 370.880005 0.000000 
0.000000 0.000000 371.989990 371.920013 371.850006 371.779999 371.279999 370.850006 370.630005 370.799988 370.709991 370.429993 0.000000 
0.000000 0.000000 372.269989 371.950012 371.959991 371.670013 371.059998 370.380005 370.019989 370.119995 370.079987 369.929993 0.000000 
0.000000 0.000000 372.519989 371.989990 371.809998 371.390015 370.959991 369.859985 369.410004 369.100006 368.959991 368.980011 0.000000 
0.000000 372.649994 372.670013 372.079987 371.809998 371.179993 370.769989 369.269989 368.480011 367.910004 367.920013 368.040009 0.000000 
372.209991 372.230011 372.209991 372.160004 371.239990 370.950012 369.970001 368.559998 367.350006 367.089996 367.679993 367.730011 0.000000 
0.000000 372.040009 372.089996 371.980011 370.790009 369.529999 368.519989 367.809998 367.109985 367.029999 367.109985 367.250000 0.000000 
0.000000 0.000000 371.519989 371.510010 371.230011 368.339996 367.100006 367.089996 366.950012 366.899994 367.320007 367.350006 0.000000 

17.5 Nutrient Files

There are a few files that specify the nutrient models setup outside of the mapping table files (.cmt, .smt). These are primarily the aquifer environment file and the point source file.

These file formats are for GSSHA v3.0b.


The Aquifer Environment File

The aquifer environment file is set up to look much like a mapping table. The file has the header

AQUATIC_ENV_NSM10

followed by two index map specifier cards. These index maps do not need to be referenced in the mapping table files to be valid here.

OV_MAP "ov_map.idx" "map name here"

ST_MAP "st_map.idx" "map name here"

These index map specifiers are followed by the more typical mapping table entries:

NUM_IDS ##

ID ##

The main difference here being that the ID lines do not have a list of parameter values, but rather the parameter values follow on the next line or lines.

There are two options for specifying the aquifer environment paramter values. The first is to simply put:

DEFAULT

on the line following the ID. This assumes that all default parameter values are applicable. The other option is to specify all of the parameter values. This comprises 11 lines of values.

...


The Point Source Files

There are two point source formats, one that contains constant point source inputs and another that contains time variable point and non-point sources.

The time varying point and non-point source file can actually specify constant values as well, and is somewhat easier to work with when specifying several points, and must be used for any non-point source loadings. Please refer to the Point/Non-Point Constituent Source File section of this manual for more information on setting up this file.

The constant point source file is a little more difficult to work with in that you need to know the exact order of all the constituents as they are in GSSHA. The order is as follows: the generic constituents (in the same order as in the mapping table file) NO2, NO3, NH4, ON, OP, DP, ALG, CBOD, DO. There are 9 NSM constituents. For example, if you have two generic constituents and the NSM constituents, you should have 11 concentrations specified.

[# of point source locations]
[node] [link] [Q] [C1] [C2] [C3] ... [Cn]
[node] [link] [Q] [C1] [C2] [C3] ... [Cn]
...
[node] [link] [Q] [C1] [C2] [C3] ... [Cn]

The order of the concentrations is: first the simple constituents in their listed order in the mapping table file, then the nutrient concentrations in the following order: NO2, NO3, NH4, Organic Nitrogen, Organic Phosphorous, Dissolved Phosphorous, PO4, Algae, CBOD, Dissolved Oxygen.

To add a point source file to the project file, include the card:

CHAN_CON_INPUT filename.ext

The units of Q are cms, and the units of concentration are mg/l (ppm).


17.6 Time and Elevation Series Files

Time Series

Time and elevation series can be used in many places in GSSHA to specify time-varying model inputs. For example, stream boundary conditions, overland flow boundary conditions, and well pumping rates. Time series are set up in separate files from the rest of the project. These files may have one or more time series in them, and there may be more than one file. To include a time series in the project, create it in a text file and include it in the project file using this card:


TIME_SERIES_FILE "filename.ts"


This is the only card that may be repeated in the project file so that you can include as many time series files as needed.


Time series are simply listings of date/time and value but there are two types of time series that may be specified -- relative and absolute. These differ in how the date/time is interpreted. For relative time series, the date and time specified is added to the simulation starting time; for absolute the values are simply converted to julien times internally to be compared against the current simulation time. Relative time series work best for short project where you are initializing some feature or testing against a hypothetical set of data. Absolute time series are better for observed data and long run times.

The time series begins with the following card:


GSSHA_TS


next comes the name of the time series. The time series will be referenced by this name throughout the project.


time series name


Then either the card ABSOLUTE or RELATIVE follows to indicate how to treat the date/time formats. (See above.) After the ABSOLUTE or RELATIVE card, the date/time and value lines appear, as many as needed.

year month day hour minute value
year month day hour minute value
year month day hour minute value
...
year month day hour minute value


Finally the card END_TS indicates the end of the date/time and value lines. There may be as many time series in one file as desired; simply place empty lines between the time series.

Here is an example time series


GSSHA_TS
Storm1
RELATIVE
0000 0 0 0 0 0.0
0000 0 0 0 30 1.0
0000 0 0 1 0 2.0
0000 0 0 1 30 3.0
0000 0 0 2 0 2.0
0000 0 0 2 30 1.0
0000 0 0 3 0 0.0
END_TS


Elevation Series