GWT-LKT
Structure of Blocks
FOR EACH SIMULATION
BEGIN OPTIONS
[FLOW_PACKAGE_NAME <flow_package_name>]
[AUXILIARY <auxiliary(naux)>]
[FLOW_PACKAGE_AUXILIARY_NAME <flow_package_auxiliary_name>]
[BOUNDNAMES]
[PRINT_INPUT]
[PRINT_CONCENTRATION]
[PRINT_FLOWS]
[SAVE_FLOWS]
[CONCENTRATION FILEOUT <concfile>]
[BUDGET FILEOUT <budgetfile>]
[BUDGETCSV FILEOUT <budgetcsvfile>]
[TS6 FILEIN <ts6_filename>]
[OBS6 FILEIN <obs6_filename>]
END OPTIONS
BEGIN PACKAGEDATA
<ifno> <strt> [<aux(naux)>] [<boundname>]
<ifno> <strt> [<aux(naux)>] [<boundname>]
...
END PACKAGEDATA
FOR ANY STRESS PERIOD
BEGIN PERIOD <iper>
<ifno> <laksetting>
<ifno> <laksetting>
...
END PERIOD
Explanation of Variables
Block: OPTIONS
flow_package_namekeyword to specify the name of the corresponding flow package. If not specified, then the corresponding flow package must have the same name as this advanced transport package (the name associated with this package in the GWT name file).auxiliarydefines an array of one or more auxiliary variable names. There is no limit on the number of auxiliary variables that can be provided on this line; however, lists of information provided in subsequent blocks must have a column of data for each auxiliary variable name defined here. The number of auxiliary variables detected on this line determines the value for naux. Comments cannot be provided anywhere on this line as they will be interpreted as auxiliary variable names. Auxiliary variables may not be used by the package, but they will be available for use by other parts of the program. The program will terminate with an error if auxiliary variables are specified on more than one line in the options block.flow_package_auxiliary_namekeyword to specify the name of an auxiliary variable in the corresponding flow package. If specified, then the simulated concentrations from this advanced transport package will be copied into the auxiliary variable specified with this name. Note that the flow package must have an auxiliary variable with this name or the program will terminate with an error. If the flows for this advanced transport package are read from a file, then this option will have no effect.BOUNDNAMESkeyword to indicate that boundary names may be provided with the list of lake cells.PRINT_INPUTkeyword to indicate that the list of lake information will be written to the listing file immediately after it is read.PRINT_CONCENTRATIONkeyword to indicate that the list of lake concentration will be printed to the listing file for every stress period in which “CONCENTRATION PRINT” is specified in Output Control. If there is no Output Control option and PRINT_CONCENTRATION is specified, then concentration are printed for the last time step of each stress period.PRINT_FLOWSkeyword to indicate that the list of lake flow rates will be printed to the listing file for every stress period time step in which “BUDGET PRINT” is specified in Output Control. If there is no Output Control option and “PRINT_FLOWS” is specified, then flow rates are printed for the last time step of each stress period.SAVE_FLOWSkeyword to indicate that lake flow terms will be written to the file specified with “BUDGET FILEOUT” in Output Control.CONCENTRATIONkeyword to specify that record corresponds to concentration.concfilename of the binary output file to write concentration information.BUDGETkeyword to specify that record corresponds to the budget.FILEOUTkeyword to specify that an output filename is expected next.budgetfilename of the binary output file to write budget information.BUDGETCSVkeyword to specify that record corresponds to the budget CSV.budgetcsvfilename of the comma-separated value (CSV) output file to write budget summary information. A budget summary record will be written to this file for each time step of the simulation.TS6keyword to specify that record corresponds to a time-series file.FILEINkeyword to specify that an input filename is expected next.ts6_filenamedefines a time-series file defining time series that can be used to assign time-varying values. See the “Time-Variable Input” section for instructions on using the time-series capability.OBS6keyword to specify that record corresponds to an observations file.obs6_filenamename of input file to define observations for the LKT package. See the “Observation utility” section for instructions for preparing observation input files. Tables ref{table:gwf-obstypetable} and ref{table:gwt-obstypetable} lists observation type(s) supported by the LKT package.
Block: PACKAGEDATA
ifnointeger value that defines the feature (lake) number associated with the specified PACKAGEDATA data on the line. IFNO must be greater than zero and less than or equal to NLAKES. Lake information must be specified for every lake or the program will terminate with an error. The program will also terminate with an error if information for a lake is specified more than once.strtreal value that defines the starting concentration for the lake.auxrepresents the values of the auxiliary variables for each lake. The values of auxiliary variables must be present for each lake. The values must be specified in the order of the auxiliary variables specified in the OPTIONS block. If the package supports time series and the Options block includes a TIMESERIESFILE entry (see the “Time-Variable Input” section), values can be obtained from a time series by entering the time-series name in place of a numeric value.boundnamename of the lake cell. BOUNDNAME is an ASCII character variable that can contain as many as 40 characters. If BOUNDNAME contains spaces in it, then the entire name must be enclosed within single quotes.
Block: PERIOD
iperinteger value specifying the starting stress period number for which the data specified in the PERIOD block apply. IPER must be less than or equal to NPER in the TDIS Package and greater than zero. The IPER value assigned to a stress period block must be greater than the IPER value assigned for the previous PERIOD block. The information specified in the PERIOD block will continue to apply for all subsequent stress periods, unless the program encounters another PERIOD block.ifnointeger value that defines the feature (lake) number associated with the specified PERIOD data on the line. IFNO must be greater than zero and less than or equal to NLAKES.laksettingline of information that is parsed into a keyword and values. Keyword values that can be used to start the LAKSETTING string include: STATUS, CONCENTRATION, RAINFALL, EVAPORATION, RUNOFF, EXT-INFLOW, and AUXILIARY. These settings are used to assign the concentration of associated with the corresponding flow terms. Concentrations cannot be specified for all flow terms. For example, the Lake Package supports a “WITHDRAWAL” flow term. If this withdrawal term is active, then water will be withdrawn from the lake at the calculated concentration of the lake.STATUS <status> CONCENTRATION <concentration> RAINFALL <rainfall> EVAPORATION <evaporation> RUNOFF <runoff> EXT-INFLOW <ext-inflow> AUXILIARY <auxname> <auxval>
statuskeyword option to define lake status. STATUS can be ACTIVE, INACTIVE, or CONSTANT. By default, STATUS is ACTIVE, which means that concentration will be calculated for the lake. If a lake is inactive, then there will be no solute mass fluxes into or out of the lake and the inactive value will be written for the lake concentration. If a lake is constant, then the concentration for the lake will be fixed at the user specified value.concentrationreal or character value that defines the concentration for the lake. The specified CONCENTRATION is only applied if the lake is a constant concentration lake. If the Options block includes a TIMESERIESFILE entry (see the “Time-Variable Input” section), values can be obtained from a time series by entering the time-series name in place of a numeric value.rainfallreal or character value that defines the rainfall solute concentration (ML-3) for the lake. If the Options block includes a TIMESERIESFILE entry (see the “Time-Variable Input” section), values can be obtained from a time series by entering the time-series name in place of a numeric value.evaporationreal or character value that defines the concentration of evaporated water (ML-3) for the lake. If this concentration value is larger than the simulated concentration in the lake, then the evaporated water will be removed at the same concentration as the lake. If the Options block includes a TIMESERIESFILE entry (see the “Time-Variable Input” section), values can be obtained from a time series by entering the time-series name in place of a numeric value.runoffreal or character value that defines the concentration of runoff (ML-3) for the lake. Value must be greater than or equal to zero. If the Options block includes a TIMESERIESFILE entry (see the “Time-Variable Input” section), values can be obtained from a time series by entering the time-series name in place of a numeric value.ext-inflowreal or character value that defines the concentration of external inflow (ML-3) for the lake. Value must be greater than or equal to zero. If the Options block includes a TIMESERIESFILE entry (see the “Time-Variable Input” section), values can be obtained from a time series by entering the time-series name in place of a numeric value.AUXILIARYkeyword for specifying auxiliary variable.auxnamename for the auxiliary variable to be assigned AUXVAL. AUXNAME must match one of the auxiliary variable names defined in the OPTIONS block. If AUXNAME does not match one of the auxiliary variable names defined in the OPTIONS block the data are ignored.auxvalvalue for the auxiliary variable. If the Options block includes a TIMESERIESFILE entry (see the “Time-Variable Input” section), values can be obtained from a time series by entering the time-series name in place of a numeric value.
Example Input File
BEGIN OPTIONS
AUXILIARY aux1 aux2
BOUNDNAMES
PRINT_INPUT
PRINT_CONCENTRATION
PRINT_FLOWS
SAVE_FLOWS
CONCENTRATION FILEOUT gwt_lkt_02.lkt.bin
BUDGET FILEOUT gwt_lkt_02.lkt.bud
OBS6 FILEIN gwt_lkt_02.lkt.obs
END OPTIONS
BEGIN PACKAGEDATA
# ifno STRT aux1 aux2 bname
1 0.00000000 99.00000000 999.00000000 MYLAKE1
2 0.00000000 99.00000000 999.00000000 MYLAKE2
3 0.00000000 99.00000000 999.00000000 MYLAKE3
END PACKAGEDATA
BEGIN PERIOD 1
1 STATUS ACTIVE
2 STATUS ACTIVE
3 STATUS ACTIVE
END PERIOD 1
Available Observation Types
| Stress Package | Observation Type | ID1 | ID2 | Description |
|---|---|---|---|---|
| LKT | concentration | ifno or boundname | -- | Lake concentration. If boundname is specified, boundname must be unique for each lake. |
| LKT | flow-ja-face | ifno or boundname | ifno or -- | Mass flow between two lakes connected by an outlet. If more than one outlet is used to connect the same two lakes, then the mass flow for only the first outlet can be observed. If a boundname is specified for ID1, then the result is the total mass flow for all outlets for a lake. If a boundname is specified for ID1 then ID2 is not used. |
| LKT | storage | ifno or boundname | -- | Simulated mass storage flow rate for a lake or group of lakes. |
| LKT | constant | ifno or boundname | -- | Simulated mass constant-flow rate for a lake or group of lakes. |
| LKT | from-mvr | ifno or boundname | -- | Simulated mass inflow into a lake or group of lakes from the MVT package. Mass inflow is calculated as the product of provider concentration and the mover flow rate. |
| LKT | to-mvr | outletno or boundname | -- | Mass outflow from a lake outlet, a lake, or a group of lakes that is available for the MVR package. If boundname is not specified for ID, then the outflow available for the MVR package from a specific lake outlet is observed. In this case, ID is the outlet number, which must be between 1 and NOUTLETS. |
| LKT | lkt | ifno or boundname | iconn or -- |
Mass flow rate for a lake or group of lakes and its aquifer connection(s). If boundname is not specified for ID, then the simulated lake-aquifer flow rate at a specific lake connection is observed. In this case, ID2 must be specified and is the connection number iconn for lake ifno. |
| LKT | rainfall | ifno or boundname | -- | Rainfall rate applied to a lake or group of lakes multiplied by the rainfall concentration. |
| LKT | evaporation | ifno or boundname | -- | Simulated evaporation rate from a lake or group of lakes multiplied by the evaporation concentration. |
| LKT | runoff | ifno or boundname | -- | Runoff rate applied to a lake or group of lakes multiplied by the runoff concentration. |
| LKT | ext-inflow | ifno or boundname | -- | Mass inflow into a lake or group of lakes calculated as the external inflow rate multiplied by the inflow concentration. |
| LKT | withdrawal | ifno or boundname | -- | Specified withdrawal rate from a lake or group of lakes multiplied by the simulated lake concentration. |
| LKT | ext-outflow | ifno or boundname | -- | External outflow from a lake or a group of lakes, through their outlets, to an external boundary. If the water mover is active, the reported ext-outflow value plus the rate to mover is equal to the total outlet outflow. |
Example Observation Input File
BEGIN options
DIGITS 7
PRINT_INPUT
END options
BEGIN continuous FILEOUT gwt_lkt02.lkt.obs.csv
# obs_name obs_type ID ID2
lkt-1-conc CONCENTRATION 1
lkt-1-extinflow EXT-INFLOW 1
lkt-1-rain RAINFALL 1
lkt-1-roff RUNOFF 1
lkt-1-evap EVAPORATION 1
lkt-1-wdrl WITHDRAWAL 1
lkt-1-stor STORAGE 1
lkt-1-const CONSTANT 1
lkt-1-gwt1 LKT 1 1
lkt-1-gwt2 LKT 1 2
lkt-2-gwt1 LKT 2 1
lkt-1-mylake1 LKT MYLAKE1
lkt-1-fjf FLOW-JA-FACE 1 2
lkt-2-fjf FLOW-JA-FACE 2 1
lkt-3-fjf FLOW-JA-FACE 2 3
lkt-4-fjf FLOW-JA-FACE 3 2
lkt-5-fjf FLOW-JA-FACE MYLAKE1
lkt-6-fjf FLOW-JA-FACE MYLAKE2
lkt-7-fjf FLOW-JA-FACE MYLAKE3
END continuous