Latest revision as of 18:51, 20 May 2008

17 - GSSHA File Formats

GSSHA User's Manual

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
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	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

Difference between revisions of "Template:GUM17"