Mapping Table:Constituent Mapping Tables

From Gsshawiki
Revision as of 21:45, 9 January 2017 by Cdowner (talk | contribs) (Contaminant Transport Table Format)

(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to: navigation, search

Simple constituents that can be uptaken and can decay, but not converted into other constituents being simulated, are created using the contaminant transport mapping table. This table is a part of the grid mapping table file (.cmt)

Mapping Table Parameters

The contaminant transport table is a bit larger than most other tables but all parameters for each index entry are on one line.

Table Name # Values Parameter Units
CONTAMINANT_TRANSPORT 8 Overland Dispersion Coefficient m2s-1
Soil and Overland Decay Coefficient d-1
Overland Uptake Coefficient m d-1
Overland Loading Kg/cell or (mg/Kg) soil
groundwater concentration mg/L
initial concentration mg/L
soil water partition in the soil L Kg-1
maximum concentration (solubility) mg/L

Contaminant Transport Table Format

The contaminant transport table is somewhat different from the other tables in that it is more like a group of tables than a single table. In addition to the NUM_IDS card, there is a NUM_CONTAM card that identifies how many constituents are being modeled. Each constituent has its own 1) short description, 2) index map, and 3) base output file name. Additionally, the concentration of each constituent in the rainfall must be specified (constant value). This is done through the PRECIP_CONC card. The PRECIP_CONC is followed by the soil water partition coefficient (L Kg-1) for the aqueous phase (overland and channel domains) with the PARTITION card. These two values apply for the contaminant over the entire domain. Other parameters, listed in the table above, are distributed over the domain using the standard mapping table format.

The table header information consists of the table card and the NUM_CONTAM card

CONTAMINANT_TRANSPORT
NUM_CONTAM  ##

After the table header information, each constituent has its own table. Each table starts with a line that has the contaminant description, the contaminant index map, and the base file name for output files for that contaminant. Following that line, the PRECIP_CONC for the constituent, the PARTITION card, and then the NUM_IDS card, followed by a text line, followed by standard index entries in the mapping table. For the parameter "load" the value is the total loading in the cell (Kg). If the project file card SOIL_CONTAM is included, the value of load is the concentration in the soil (mg/Kg). This concentration applies to the first layer, that is, the MIXING_LAYER_DEPTH. Each of these lines are repeated for each constituent you wish to model.

Note: The load for the first constituent can be specified with a map by specifying the map in the project file with the CONTAM_MAP card, where a load is specified for each grid cell. As described above, the load is total mass in the cell (Kg), unless the SOIL_CONTAM card is specified, then the load is soil concentration (mg/Kg). All other values for constituent one come from the mapping table. A value for load must still be supplied in the table, even though this value is not used in GSSHA.

"Contam 1 description" "index_map_1" "contam1.out"
PRECIP_CONC  ##.##
PARTITION ##.##
NUM_IDS  ##
Text line...
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
...
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
"Contam 2 description" "index_map_2" "contam2.out"
PRECIP_CONC  ##.## 
PARTITION ##.##
NUM_IDS  ##
Text line...
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
...
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
...
"Contam N description" "index_map_N" "contamN.out"
PRECIP_CONC  ##.## 
PARTITION ##.##
NUM_IDS  ##
Text line...
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]
...
[ID]    [DESC...]           ... [Disp]  [decay] [uptake]  [load]  [gw_conc]  [init_conc]  [sw_part]  [max_conc]


An example of two contaminants is shown below.

CONTAMINANT_TRANSPORT
NUM_CONTAM  2
"Industrial constituent" "mixed_land_use" "industrial1.out"
PRECIP_CONC  0.1
PARTITION 91.4
NUM_IDS 3
ID   Desc...                     ...   disp   decay  uptake load  gw_conc  init_conc  sw_part max_conc
1    forest areas, natural             0.001  0.01   0.1    2.0    3.0       0.0      91.4     10000.0
2    residential areas                 0.001  0.001  0.1    10.0   8.0       2.0      91.4     10000.0
3    industrial areas                  0.001  0.001  0.1    50.0   23.0      8.0      91.4     10000.0
"Agricultural constituent" "rangeland" "ag2.out"
PRECIP_CONC  2.0
PARTITION 1.0
NUM_IDS 2
ID   Desc...                     ...   disp   decay  uptake load  gw_conc  init_conc  sw_part max_conc
1    grazed area                       0.001  0.01   0.1    50.0   6.0       2.0      10.0     10000.0
2    other areas                       0.001  0.01   0.1    3.0    2.0       0.0      10.0     10000.0


GSSHA User's Manual

12 Mapping Table
12.1     File Description
12.2     Index Maps
12.3     Grid-based Mapping Tables
12.4     ID Line Format
12.5     Example Mapping Table File
12.6     Stream Mapping Tables
12.7     Sediment Erosion Mapping Tables
12.7     Constituent Mapping Tables