GWF-LAK
Structure of Blocks
FOR EACH SIMULATION
BEGIN OPTIONS
[AUXILIARY <auxiliary(naux)>]
[BOUNDNAMES]
[PRINT_INPUT]
[PRINT_STAGE]
[PRINT_FLOWS]
[SAVE_FLOWS]
[STAGE FILEOUT <stagefile>]
[BUDGET FILEOUT <budgetfile>]
[BUDGETCSV FILEOUT <budgetcsvfile>]
[PACKAGE_CONVERGENCE FILEOUT <package_convergence_filename>]
[TS6 FILEIN <ts6_filename>]
[OBS6 FILEIN <obs6_filename>]
[MOVER]
[SURFDEP <surfdep>]
[MAXIMUM_ITERATIONS <maximum_iterations>]
[MAXIMUM_STAGE_CHANGE <maximum_stage_change>]
[TIME_CONVERSION <time_conversion>]
[LENGTH_CONVERSION <length_conversion>]
END OPTIONS
BEGIN DIMENSIONS
NLAKES <nlakes>
NOUTLETS <noutlets>
NTABLES <ntables>
END DIMENSIONS
BEGIN PACKAGEDATA
<ifno> <strt> <nlakeconn> [<aux(naux)>] [<boundname>]
<ifno> <strt> <nlakeconn> [<aux(naux)>] [<boundname>]
...
END PACKAGEDATA
BEGIN CONNECTIONDATA
<ifno> <iconn> <cellid(ncelldim)> <claktype> <bedleak> <belev> <telev> <connlen> <connwidth>
<ifno> <iconn> <cellid(ncelldim)> <claktype> <bedleak> <belev> <telev> <connlen> <connwidth>
...
END CONNECTIONDATA
BEGIN TABLES
<ifno> TAB6 FILEIN <tab6_filename>
<ifno> TAB6 FILEIN <tab6_filename>
...
END TABLES
BEGIN OUTLETS
<outletno> <lakein> <lakeout> <couttype> <invert> <width> <rough> <slope>
<outletno> <lakein> <lakeout> <couttype> <invert> <width> <rough> <slope>
...
END OUTLETS
FOR ANY STRESS PERIOD
BEGIN PERIOD <iper>
<number> <laksetting>
<number> <laksetting>
...
END PERIOD
Explanation of Variables
Block: OPTIONS
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.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_STAGEkeyword to indicate that the list of lake stages will be printed to the listing file for every stress period in which “HEAD PRINT” is specified in Output Control. If there is no Output Control option and PRINT_STAGE is specified, then stages 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.STAGEkeyword to specify that record corresponds to stage.stagefilename of the binary output file to write stage 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.PACKAGE_CONVERGENCEkeyword to specify that record corresponds to the package convergence comma spaced values file.package_convergence_filenamename of the comma spaced values output file to write package convergence information.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 LAK 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 LAK package.MOVERkeyword to indicate that this instance of the LAK Package can be used with the Water Mover (MVR) Package. When the MOVER option is specified, additional memory is allocated within the package to store the available, provided, and received water.surfdepreal value that defines the surface depression depth for VERTICAL lake-GWF connections. If specified, SURFDEP must be greater than or equal to zero. If SURFDEP is not specified, a default value of zero is used for all vertical lake-GWF connections.maximum_iterationsinteger value that defines the maximum number of Newton-Raphson iterations allowed for a lake. By default, MAXIMUM_ITERATIONS is equal to 100. MAXIMUM_ITERATIONS would only need to be increased from the default value if one or more lakes in a simulation has a large water budget error.maximum_stage_changereal value that defines the lake stage closure tolerance. By default, MAXIMUM_STAGE_CHANGE is equal to 1 x 10-5. The MAXIMUM_STAGE_CHANGE would only need to be increased or decreased from the default value if the water budget error for one or more lakes is too small or too large, respectively.time_conversionreal value that is used to convert user-specified Manning’s roughness coefficients or gravitational acceleration used to calculate outlet flows from seconds to model time units. TIME_CONVERSION should be set to 1.0, 60.0, 3,600.0, 86,400.0, and 31,557,600.0 when using time units (TIME_UNITS) of seconds, minutes, hours, days, or years in the simulation, respectively. CONVTIME does not need to be specified if no lake outlets are specified or TIME_UNITS are seconds.length_conversionreal value that is used to convert outlet user-specified Manning’s roughness coefficients or gravitational acceleration used to calculate outlet flows from meters to model length units. LENGTH_CONVERSION should be set to 3.28081, 1.0, and 100.0 when using length units (LENGTH_UNITS) of feet, meters, or centimeters in the simulation, respectively. LENGTH_CONVERSION does not need to be specified if no lake outlets are specified or LENGTH_UNITS are meters.
Block: DIMENSIONS
nlakesvalue specifying the number of lakes that will be simulated for all stress periods.noutletsvalue specifying the number of outlets that will be simulated for all stress periods. If NOUTLETS is not specified, a default value of zero is used.ntablesvalue specifying the number of lakes tables that will be used to define the lake stage, volume relation, and surface area. If NTABLES is not specified, a default value of zero is used.
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 stage for the lake.nlakeconninteger value that defines the number of GWF cells connected to this (IFNO) lake. There can only be one vertical lake connection to each GWF cell. NLAKECONN must be greater than zero.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: CONNECTIONDATA
ifnointeger value that defines the feature (lake) number associated with the specified CONNECTIONDATA data on the line. IFNO must be greater than zero and less than or equal to NLAKES. Lake connection information must be specified for every lake connection to the GWF model (NLAKECONN) or the program will terminate with an error. The program will also terminate with an error if connection information for a lake connection to the GWF model is specified more than once.iconninteger value that defines the GWF connection number for this lake connection entry. ICONN must be greater than zero and less than or equal to NLAKECONN for lake IFNO.cellidis the cell identifier, and depends on the type of grid that is used for the simulation. For a structured grid that uses the DIS input file, CELLID is the layer, row, and column. For a grid that uses the DISV input file, CELLID is the layer and CELL2D number. If the model uses the unstructured discretization (DISU) input file, CELLID is the node number for the cell.claktypecharacter string that defines the lake-GWF connection type for the lake connection. Possible lake-GWF connection type strings include: VERTICAL–character keyword to indicate the lake-GWF connection is vertical and connection conductance calculations use the hydraulic conductivity corresponding to the K_{33} tensor component defined for CELLID in the NPF package. HORIZONTAL–character keyword to indicate the lake-GWF connection is horizontal and connection conductance calculations use the hydraulic conductivity corresponding to the K_{11} tensor component defined for CELLID in the NPF package. EMBEDDEDH–character keyword to indicate the lake-GWF connection is embedded in a single cell and connection conductance calculations use the hydraulic conductivity corresponding to the K_{11} tensor component defined for CELLID in the NPF package. EMBEDDEDV–character keyword to indicate the lake-GWF connection is embedded in a single cell and connection conductance calculations use the hydraulic conductivity corresponding to the K_{33} tensor component defined for CELLID in the NPF package. Embedded lakes can only be connected to a single cell (NLAKECONN = 1) and there must be a lake table associated with each embedded lake.bedleakreal value or character string that defines the bed leakance for the lake-GWF connection. BEDLEAK must be greater than or equal to zero, equal to the DNODATA value (3.0E+30), or specified to be NONE. If DNODATA or NONE is specified for BEDLEAK, the lake-GWF connection conductance is solely a function of aquifer properties in the connected GWF cell and lakebed sediments are assumed to be absent. Warning messages will be issued if NONE is specified. Eventually the ability to specify NONE will be deprecated and cause MODFLOW 6 to terminate with an error.belevreal value that defines the bottom elevation for a HORIZONTAL lake-GWF connection. Any value can be specified if CLAKTYPE is VERTICAL, EMBEDDEDH, or EMBEDDEDV. If CLAKTYPE is HORIZONTAL and BELEV is not equal to TELEV, BELEV must be greater than or equal to the bottom of the GWF cell CELLID. If BELEV is equal to TELEV, BELEV is reset to the bottom of the GWF cell CELLID.televreal value that defines the top elevation for a HORIZONTAL lake-GWF connection. Any value can be specified if CLAKTYPE is VERTICAL, EMBEDDEDH, or EMBEDDEDV. If CLAKTYPE is HORIZONTAL and TELEV is not equal to BELEV, TELEV must be less than or equal to the top of the GWF cell CELLID. If TELEV is equal to BELEV, TELEV is reset to the top of the GWF cell CELLID.connlenreal value that defines the distance between the connected GWF CELLID node and the lake for a HORIZONTAL, EMBEDDEDH, or EMBEDDEDV lake-GWF connection. CONLENN must be greater than zero for a HORIZONTAL, EMBEDDEDH, or EMBEDDEDV lake-GWF connection. Any value can be specified if CLAKTYPE is VERTICAL.connwidthreal value that defines the connection face width for a HORIZONTAL lake-GWF connection. CONNWIDTH must be greater than zero for a HORIZONTAL lake-GWF connection. Any value can be specified if CLAKTYPE is VERTICAL, EMBEDDEDH, or EMBEDDEDV.
Block: TABLES
ifnointeger value that defines the feature (lake) number associated with the specified TABLES data on the line. IFNO must be greater than zero and less than or equal to NLAKES. The program will terminate with an error if table information for a lake is specified more than once or the number of specified tables is less than NTABLES.TAB6keyword to specify that record corresponds to a table file.FILEINkeyword to specify that an input filename is expected next.tab6_filenamecharacter string that defines the path and filename for the file containing lake table data for the lake connection. The TAB6_FILENAME file includes the number of entries in the file and the relation between stage, volume, and surface area for each entry in the file. Lake table files for EMBEDDEDH and EMBEDDEDV lake-GWF connections also include lake-GWF exchange area data for each entry in the file. Instructions for creating the TAB6_FILENAME input file are provided in Lake Table Input File section.
Block: OUTLETS
outletnointeger value that defines the outlet number associated with the specified OUTLETS data on the line. OUTLETNO must be greater than zero and less than or equal to NOUTLETS. Outlet information must be specified for every outlet or the program will terminate with an error. The program will also terminate with an error if information for a outlet is specified more than once.lakeininteger value that defines the lake number that outlet is connected to. LAKEIN must be greater than zero and less than or equal to NLAKES.lakeoutinteger value that defines the lake number that outlet discharge from lake outlet OUTLETNO is routed to. LAKEOUT must be greater than or equal to zero and less than or equal to NLAKES. If LAKEOUT is zero, outlet discharge from lake outlet OUTLETNO is discharged to an external boundary.couttypecharacter string that defines the outlet type for the outlet OUTLETNO. Possible COUTTYPE strings include: SPECIFIED–character keyword to indicate the outlet is defined as a specified flow. MANNING–character keyword to indicate the outlet is defined using Manning’s equation. WEIR–character keyword to indicate the outlet is defined using a sharp weir equation.invertreal value that defines the invert elevation for the lake outlet. Any value can be specified if COUTTYPE is SPECIFIED. 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.widthreal value that defines the width of the lake outlet. Any value can be specified if COUTTYPE is SPECIFIED. 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.roughreal value that defines the roughness coefficient for the lake outlet. Any value can be specified if COUTTYPE is not MANNING. 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.slopereal value that defines the bed slope for the lake outlet. Any value can be specified if COUTTYPE is not MANNING. 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.
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.numberinteger value that defines the lake or outlet number associated with the specified PERIOD data on the line. NUMBER must be greater than zero and less than or equal to NLAKES for a lake number and less than or equal to NOUTLETS for an outlet number.laksettingline of information that is parsed into a keyword and values. Keyword values that can be used to start the LAKSETTING string include both keywords for lake settings and keywords for outlet settings. Keywords for lake settings include: STATUS, STAGE, RAINFALL, EVAPORATION, RUNOFF, INFLOW, WITHDRAWAL, and AUXILIARY. Keywords for outlet settings include RATE, INVERT, WIDTH, SLOPE, and ROUGH.STATUS <status> STAGE <stage> RAINFALL <rainfall> EVAPORATION <evaporation> RUNOFF <runoff> INFLOW <inflow> WITHDRAWAL <withdrawal> RATE <rate> INVERT <invert> WIDTH <width> SLOPE <slope> ROUGH <rough> AUXILIARY <auxname> <auxval>
statuskeyword option to define lake status. STATUS can be ACTIVE, INACTIVE, or CONSTANT. By default, STATUS is ACTIVE.stagereal or character value that defines the stage for the lake. The specified STAGE is only applied if the lake is a constant stage 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 rate (LT-1) 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.evaporationreal or character value that defines the maximum evaporation rate (LT-1) 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.runoffreal or character value that defines the runoff rate (L3 T-1) 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.inflowreal or character value that defines the volumetric inflow rate (L3 T-1) 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. By default, inflow rates are zero for each lake.withdrawalreal or character value that defines the maximum withdrawal rate (L3 T-1) 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.ratereal or character value that defines the extraction rate for the lake outflow. A positive value indicates inflow and a negative value indicates outflow from the lake. RATE only applies to outlets associated with active lakes (STATUS is ACTIVE). A specified RATE is only applied if COUTTYPE for the OUTLETNO is SPECIFIED. 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. By default, the RATE for each SPECIFIED lake outlet is zero.invertreal or character value that defines the invert elevation for the lake outlet. A specified INVERT value is only used for active lakes if COUTTYPE for lake outlet OUTLETNO is not SPECIFIED. 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.roughreal value that defines the roughness coefficient for the lake outlet. Any value can be specified if COUTTYPE is not MANNING. 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.widthreal or character value that defines the width of the lake outlet. A specified WIDTH value is only used for active lakes if COUTTYPE for lake outlet OUTLETNO is not SPECIFIED. 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.slopereal or character value that defines the bed slope for the lake outlet. A specified SLOPE value is only used for active lakes if COUTTYPE for lake outlet OUTLETNO is MANNING. 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
PRINT_INPUT
BOUNDNAMES
PRINT_STAGE
PRINT_FLOWS
STAGE FILEOUT lak-1.stage.bin
BUDGET FILEOUT lak-1.cbc
END OPTIONS
BEGIN DIMENSIONS
NLAKES 1
NOUTLETS 1
END DIMENSIONS
BEGIN PACKAGEDATA
# ifno strt lakeconn boundname
1 110.00 57 LAKE_1
END PACKAGEDATA
BEGIN CONNECTIONDATA
# ifno iconn layer row column ctype bedleak belev telev dx width
1 1 1 7 6 HORIZONTAL 0.1 0 0 500 500
1 2 1 8 6 HORIZONTAL 0.1 0 0 500 500
1 3 1 9 6 HORIZONTAL 0.1 0 0 500 500
1 4 1 10 6 HORIZONTAL 0.1 0 0 500 500
1 5 1 11 6 HORIZONTAL 0.1 0 0 500 500
1 6 1 6 7 HORIZONTAL 0.1 0 0 500 500
1 7 2 7 7 VERTICAL 0.1 0 0 0 0
1 8 2 8 7 VERTICAL 0.1 0 0 0 0
1 9 2 8 7 HORIZONTAL 0.1 0 0 250 500
1 10 2 9 7 VERTICAL 0.1 0 0 0 0
1 11 2 9 7 HORIZONTAL 0.1 0 0 250 500
1 12 2 10 7 VERTICAL 0.1 0 0 0 0
1 13 2 10 7 HORIZONTAL 0.1 0 0 250 500
1 14 2 11 7 VERTICAL 0.1 0 0 0 0
1 15 1 12 7 HORIZONTAL 0.1 0 0 500 500
1 16 1 6 8 HORIZONTAL 0.1 0 0 500 500
1 17 2 7 8 VERTICAL 0.1 0 0 0 0
1 18 2 7 8 HORIZONTAL 0.1 0 0 250 500
1 19 3 8 8 VERTICAL 0.1 0 0 0 0
1 20 3 9 8 VERTICAL 0.1 0 0 0 0
1 21 3 10 8 VERTICAL 0.1 0 0 0 0
1 22 2 11 8 VERTICAL 0.1 0 0 0 0
1 23 2 11 8 HORIZONTAL 0.1 0 0 250 500
1 24 1 12 8 HORIZONTAL 0.1 0 0 500 500
1 25 1 6 9 HORIZONTAL 0.1 0 0 500 500
1 26 2 7 9 VERTICAL 0.1 0 0 0 0
1 27 2 7 9 HORIZONTAL 0.1 0 0 250 500
1 28 3 8 9 VERTICAL 0.1 0 0 0 0
1 29 3 9 9 VERTICAL 0.1 0 0 0 0
1 30 3 10 9 VERTICAL 0.1 0 0 0 0
1 31 2 11 9 VERTICAL 0.1 0 0 0 0
1 32 2 11 9 HORIZONTAL 0.1 0 0 250 500
1 33 1 12 9 HORIZONTAL 0.1 0 0 500 500
1 34 1 6 10 HORIZONTAL 0.1 0 0 500 500
1 35 2 7 10 VERTICAL 0.1 0 0 0 0
1 36 2 7 10 HORIZONTAL 0.1 0 0 250 500
1 37 3 8 10 VERTICAL 0.1 0 0 0 0
1 38 3 9 10 VERTICAL 0.1 0 0 0 0
1 39 3 10 10 VERTICAL 0.1 0 0 0 0
1 40 2 11 10 VERTICAL 0.1 0 0 0 0
1 41 2 11 10 HORIZONTAL 0.1 0 0 250 500
1 42 1 12 10 HORIZONTAL 0.1 0 0 500 500
1 43 1 6 11 HORIZONTAL 0.1 0 0 500 500
1 44 2 7 11 VERTICAL 0.1 0 0 0 0
1 45 2 8 11 VERTICAL 0.1 0 0 0 0
1 46 2 8 11 HORIZONTAL 0.1 0 0 250 500
1 47 2 9 11 VERTICAL 0.1 0 0 0 0
1 48 2 9 11 HORIZONTAL 0.1 0 0 250 500
1 49 2 10 11 VERTICAL 0.1 0 0 0 0
1 50 2 10 11 HORIZONTAL 0.1 0 0 250 500
1 51 2 11 11 VERTICAL 0.1 0 0 0 0
1 52 1 12 11 HORIZONTAL 0.1 0 0 500 500
1 53 1 7 12 HORIZONTAL 0.1 0 0 500 500
1 54 1 8 12 HORIZONTAL 0.1 0 0 500 500
1 55 1 9 12 HORIZONTAL 0.1 0 0 500 500
1 56 1 10 12 HORIZONTAL 0.1 0 0 500 500
1 57 1 11 12 HORIZONTAL 0.1 0 0 500 500
END CONNECTIONDATA
BEGIN OUTLETS
# outletno lakein lakeout couttype invert width rough slope
1 1 0 SPECIFIED 0 0 0 0
END OUTLETS
BEGIN PERIOD 1
1 RAINFALL 0.0116
1 EVAPORATION 0.0103
END PERIOD
BEGIN PERIOD 100
1 STATUS CONSTANT
1 STAGE 110.
END PERIOD
Available Observation Types
| Stress Package | Observation Type | ID1 | ID2 | Description |
|---|---|---|---|---|
| LAK | stage | ifno or boundname | -- | Surface-water stage in a lake. If boundname is specified, boundname must be unique for each lake. |
| LAK | ext-inflow | ifno or boundname | -- | Specified inflow into a lake or group of lakes. |
| LAK | outlet-inflow | ifno or boundname | -- | Simulated inflow from upstream lake outlets into a lake or group of lakes. |
| LAK | inflow | ifno or boundname | -- | Sum of specified inflow and simulated inflow from upstream lake outlets into a lake or group of lakes. |
| LAK | from-mvr | ifno or boundname | -- | Inflow into a lake or group of lakes from the MVR package. |
| LAK | rainfall | ifno or boundname | -- | Rainfall rate applied to a lake or group of lakes. |
| LAK | runoff | ifno or boundname | -- | Runoff rate applied to a lake or group of lakes. |
| LAK | lak | ifno or boundname | iconn or -- |
Simulated 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. |
| LAK | withdrawal | ifno or boundname | -- | Specified withdrawal rate from a lake or group of lakes. |
| LAK | evaporation | ifno or boundname | -- | Simulated evaporation rate from a lake or group of lakes. |
| LAK | ext-outflow | outletno or boundname | -- | External outflow from a lake outlet, a lake, or a group of lakes to an external boundary. If boundname is not specified for ID, then the external outflow from a specific lake outlet is observed. In this case, ID is the outlet number outletno. |
| LAK | to-mvr | outletno or boundname | -- | 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 outletno. |
| LAK | storage | ifno or boundname | -- | Simulated storage flow rate for a lake or group of lakes. |
| LAK | constant | ifno or boundname | -- | Simulated constant-flow rate for a lake or group of lakes. |
| LAK | outlet | outletno or boundname | -- | Simulated outlet flow rate from a lake outlet, a lake, or a group of lakes. If boundname is not specified for ID, then the flow from a specific lake outlet is observed. In this case, ID is the outlet number outletno. |
| LAK | volume | ifno or boundname | -- | Simulated lake volume or group of lakes. |
| LAK | surface-area | ifno or boundname | -- | Simulated surface area for a lake or group of lakes. |
| LAK | wetted-area | ifno or boundname | iconn or -- |
Simulated wetted-area for a lake or group of lakes and its aquifer connection(s). If boundname is not specified for ID, then the wetted area of a specific lake connection is observed. In this case, ID2 must be specified and is the connection number iconn. |
| LAK | conductance | ifno or boundname | iconn or -- |
Calculated conductance for a lake or group of lakes and its aquifer connection(s). If boundname is not specified for ID, then the calculated conductance of a specific lake connection is observed. In this case, ID2 must be specified and is the connection number iconn. |
Example Observation Input File
BEGIN OPTIONS
DIGITS 10
PRINT_INPUT
END OPTIONS
BEGIN CONTINUOUS FILEOUT my_model.lak.csv
l1stage stage 1
l1vol volume 1
vflow lak 1 1
hflow1 lak 1 2
hflow2 lak 1 3
hflow3 lak 1 4
hflow4 lak 1 5
lakflow lak lake_1
END CONTINUOUS