Difference between revisions of "Mapping Table:Constituent Mapping Tables"

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

Revision as of 20:02, 4 August 2009

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 Dispersion Coefficient m2s-1
Decay Coefficient d-1
Uptake Coefficient m d-1
Loading Kg/cell or (mg/Kg) soil
groundwater concentration mg/L
hot start concentration mg/L
soil water partition L Kg-1
maximum concentration 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. There is only one PRECIP_CONC card, and all concentrations are listed in order on that line after the PRECIP_CONC card. The PRECIP_CONC is followed by the soil water partition coefficient (L Kg-1) for the aqueous phase with the PARTITION card.

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