For Each Simulation:
|
1.
|
NLAKES ILKCB
If data are being read using the fixed format mode, then each field should be entered using I10 format.
NLAKES Number of separate lakes.
Sublakes of multiple-lake systems are considered separate lakes for input purposes. The variable NLAKES is used, with certain internal assumptions and approximations, to dimension arrays for the simulation.
|
ILKCB Whether or not to write cell-by-cell flows (yes if ILKCB> 0, no otherwise). If ILKCB< 0 and "Save Budget" is specified in the Output Control or ICBCFL is not equal to 0, the cell-by-cell flows will be printed in the standard output file.
ICBCFL is specified in the input to the Output Control Option of MODFLOW.
|
|
|
2.
|
THETA [NSSITR SSCNCR]
If data are being read using the fixed format mode, then the data should be entered using format (F10.4,I10,F10.4).
THETA Explicit (THETA = 0.0), semi-implicit (0.0 < THETA < 1.0), or implicit (THETA = 1.0) solution for lake stages.
| • | THETA is automatically set to a value of 1.0 for all steady-state stress periods. THETA for transient simulations has been revised and is now limited to range from 0.5 to 1.0. A value of 0.5 represents the stage midway between the previous time step and the end of the current time step. A value of 1.0 (fully implicit) represents the lake stage at the end of the current time step. A THETA of less than 0.5 does not perform well and a zero value is undefined in the Newton iteration method. Slight errors in the solution of lake stage and seepage may result when THETA is greater than 0.5 (Fread, 1993). Values greater than 0.5 have been recommended for damping oscillations in streamflow-routing equations (Fread, 1993). A value of 0.5 represents a semi-implicit method that is often called Crank-Nicolson (Wang and Anderson, 1982, p. 81). Wang and Anderson present results for drawdown for a simple confined aquifer with pumping in which the Crank-Nicolson method produced the best results when compared with the Theis solution for different times. A value of 0.5 is generally recommended. |
| • | In MODFLOW-2000, ISS is not part of the input. Instead NSSITR or SSCNCR should be included if one or more stress periods is a steady state stress period as defined in Ss/tr in the Discretization file. |
| • | SSCNCR and NSSITR can be read for a transient only simulation by placing a negative sign immeditately in front of THETA. A negative THETA sets a flag which assumes input values for NSSITR and SSCNCR will follow THETA in the format as described by Merritt and Konikow (p. 52). A negative THETA is automatically reset to a positive value after values of NSSITR and SSCNCR are read. |
|
NSSITR Maximum number of iterations for Newton’s method solution for equilibrium lake stages in each MODFLOW iteration for steady-state aquifer head solution. Only read if ISS (option flag input to BCF Package of MODFLOW indicating steady-state solution) is not zero.
| • | NSSITR and SSCNCR may be omitted for transient solutions (ISS = 0). |
| • | In MODFLOW-2000, ISS is not part of the input. Instead NSSITR or SSCNCR should be included if one or more stress periods is a steady state stress period as defined in Ss/tr in the Discretization file. |
| • | SSCNCR and NSSITR can be read for a transient only simulation by placing a negative sign immeditately in front of THETA. A negative THETA sets a flag which assumes input values for NSSITR and SSCNCR will follow THETA in the format as described by Merritt and Konikow (p. 52). A negative THETA is automatically reset to a positive value after values of NSSITR and SSCNCR are read. |
| • | If NSSITR = 0, a value of 100 will be used instead. |
|
SSCNCR Convergence criterion for equilibrium lake stage solution by Newton’s method. Only read if ISS is not zero. In MODFLOW-2000, ISS is not part of the input. Instead SSCNCR should be included if one or more stress periods is a steady state stress period as defined in Ss/tr in the Discretization file.
| • | NSSITR and SSCNCR may be omitted for transient solutions (ISS = 0). |
| • | In MODFLOW-2000, ISS is not part of the input. Instead NSSITR or SSCNCR should be included if one or more stress periods is a steady state stress period as defined in Ss/tr in the Discretization file. |
| • | SSCNCR and NSSITR can be read for a transient only simulation by placing a negative sign immeditately in front of THETA. A negative THETA sets a flag which assumes input values for NSSITR and SSCNCR will follow THETA in the format as described by Merritt and Konikow (p. 52). A negative THETA is automatically reset to a positive value after values of NSSITR and SSCNCR are read. |
| • | If SSCNCR = 0, a value of 0.0001 will be used instead. |
|
|
|
For the First Stress Period Only:
|
3.
|
STAGES [SSMN SSMX] [CLAKE(1)..........CLAKE(NSOL)]
This data set should consist of one line for each lake, where line 1 includes data for lake 1, and line n includes data for lake n. There must be exactly NLAKES lines of data.
If data are being read using the fixed format mode, then each field should be entered using F10.4 format.
STAGES The initial stage of each lake at the beginning of the run.
|
SSMN Minimum stage allowed for each lake in steady-state solution.
| • | SSMN and SSMX are not needed for a transient run and must be omitted when the solution is transient. |
| • | When the first stress period is a steady-state stress period, SSMN is defined in record 3. For subsequent steady-state stress periods, SSMN is defined in record 9a. |
|
SSMX Maximum stage allowed for each lake in steady-state solution.
| • | SSMN and SSMX are not needed for a transient run and must be omitted when the solution is transient. |
| • | When the first stress period is a steady-state stress period, SSMN is defined in record 3. For subsequent steady-state stress periods, SSMN is defined in record 9a. |
|
CLAKE The initial concentrations in each lake at the beginning of the model run. Values are entered for NSOL constituents. The value of NSOL is passed from MOC3D or GWT. CLAKE values are ignored if entered in MODFLOW runs.
|
|
|
For Each Stress Period:
|
4.
|
ITMP ITMP1 LWRT
If data are being read using the fixed format mode, then each field should be entered using I10 format.
Lake definition data are restricted to cells for which IBOUND and WETDRY values have been set to zero.
ITMP > 0, read lake definition data (records 5-7, and, optionally, records 8 and 9);
ITMP = 0, no lake calculations this stress period;
ITMP < 0, use lake definition data from last stress period.
|
ITMP1 ≥ 0, read new recharge, evaporation, runoff, and withdrawal data for each lake, and associated concentrations if needed for MOC3D runs;
ITMP1 < 0, use recharge, evaporation, runoff, and withdrawal data, and concentrations, if needed, from last stress period.
|
LWRT > 0, suppresses printout from the lake package.
ICBCFL ≤ 0 or not specifying "Save Budget" also suppresses printout from the lake package. ICBCFL or "Save Budget" are specified in the input to the Output Control option of MODFLOW.
|
|
|
If ITMP > 0 read lake definition data (records 5-7, and, optionally, records 8 and 9).
|
5.
|
LKARR(NCOL,NROW)
A NCOL by NROW array is read for each layer in the grid by MODFLOW module U2DINT.
LKARR A value is read in for every grid cell.
If LKARR(I,J,K) = 0, the grid cell is not a lake volume cell.
If LKARR(I,J,K) > 0, its value is the identification number of the lake occupying the grid cell. LKARR(I,J,K) must not exceed the value NLAKES. If it does, or if LKARR(I,J,K) < 0, LKARR(I,J,K) is set to zero.
Lake cells cannot be overlain by non-lake cells in a higher layer.
The Lake package can be used when all or some of the model layers containing the lake are confined. The authors recommend using the Layer-Property Flow Package (LPF) for this case, although the BCF and HUF Packages will work too. However, when using the BCF6 package to define aquifer properties, lake/aquifer conductances in the lateral direction are based solely on the lakebed leakance (and not on the lateral transmissivity of the aquifer layer). As before, when the BCF6 package is used, vertical lake/aquifer conductances are based on lakebed conductance and on the vertical hydraulic conductivity of the aquifer layer underlying the lake when the wet/dry option is implemented, and only on the lakebed leakance when the wet/dry option is not implemented.
|
|
6.
|
BDLKNC(NCOL,NROW)
A NCOL by NROW array is read for each layer in the grid by MODFLOW module U2DREL.
BDLKNC A value is read in for every grid cell. The value is the lakebed leakance that will be assigned to lake/aquifer interfaces that occur in the corresponding grid cell.
If the wet-dry option flag (IWDFLG) is not active (cells cannot rewet if they become dry), then the BDLKNC values are assumed to represent the combined leakances of the lakebed material and the aquifer material between the lake and the centers of the underlying grid cells, i. e., the vertical conductance values (CV) will not be used in the computation of conductances across lake/aquifer boundary faces in the vertical direction.
IBOUND and WETDRY should be set to zero for every cell for which LKARR is not equal to zero. IBOUND is defined in the input to the Basic Package of MODFLOW). WETDRY is defined in the input to the BCF or other flow package of MODFLOW if the IWDFLG option is active.
|
|
7.
|
NSLMS
If data are being read using the fixed format mode, then NSLMS should be entered using format I5.
NSLMS The number of sublake systems if coalescing/dividing lakes are to be simulated (only in transient runs). Enter 0 if no sublake systems are to be simulated.
|
|
If ITMP > 0 and NSLMS > 0:
A pair of records (records 8a and 8b) is read for each multiple-lake system, i.e., NSLMS pairs of records. However, IC = 0 will terminate the input.
|
8a.
|
IC ISUB(1) ISUB(2) ............ ISUB(IC)
If data are being read using the fixed format mode, then each field of Record 8a should be entered using I5 format
IC The number of sublakes, including the center lake, in the sublake system being described in this record.
|
ISUB The identification numbers of the sublakes in the sublake system being described in this record. The center lake number is listed first.
|
|
|
8b.
|
SILLVT(2) ............. SILLVT(IC)
If data are being read using the fixed format mode, then each field of Record 8b should be entered using F10.4 format
SILLVT Sill elevation that determines whether the center lake is connected with a given sublake. One value is entered in this record for each sublake in the order the sublakes are listed in the previous record.
|
|
If ITMP1 ≥ 0
At least one of the above records will be read for each lake; i.e., NLAKES records, or sets of records, will be read. If MODFLOW is being run, only the first record is read. If MOC3D or GWT is being run, a set of two or more records will be read for each lake (see next note).
If record 9b is included because solute transport is being simulated, then 9b should consist of one record (line) for each solute; each record must contain two or three values; and there must be as many records as the number of solutes being simulated (NSOL). The order of records must be that all necessary lines for 9b are listed for a given lake before line 9a for the next lake. For example, if the Lake Package is representing three lakes and the solute transport package is representing two solutes, then the order of data for Record 9 would be 9a, 9b, 9b, 9a, 9b, 9b, 9a, 9b, and 9b.
|
9a.
|
PRCPLK EVAPLK RNF WTHDRW [SSMN] [SSMX]
PRCPLK The rate of precipitation per unit area at the surface of a lake (L/T).
|
EVAPLK The rate of evaporation per unit area from the surface of a lake (L/T).
|
RNF Overland runoff from an adjacent watershed entering the lake. If RNF > 0, it is specified directly as a volumetric rate, or flux (L3 /T). If RNF < 0, its absolute value is used as a dimensionless multiplier applied to the product of the lake precipitation rate per unit area (PRCPLK) and the surface area of the lake at its full stage (occupying all layer 1 lake cells).
When RNF is entered as a dimensionless multiplier (RNF < 0), it is considered to be the product of two proportionality factors. The first is the ratio of the area of the basin contributing runoff to the surface area of the lake when it is at full stage. The second is the fraction of the current rainfall rate that becomes runoff to the lake. This procedure provides a means for the automated computation of runoff rate from a watershed to a lake as a function of varying rainfall rate. For example, if the basin area is 10 times greater than the surface area of the lake, and 20 percent of the precipitation on the basin becomes overland runoff directly into the lake, then set RNF = -2.0.
|
WTHDRW The volumetric rate, or flux (L3 /T), of water removal from a lake by means other than rainfall, evaporation, surface outflow, or ground-water seepage. A negative value indicates augmentation. Normally, this would be used to specify the rate of artificial withdrawal from a lake for human water use, or if negative, artificial augmentation of a lake volume for esthetic or recreational purposes.
|
SSMN Minimum stage allowed for each lake in steady-state solution.
| • | SSMN and SSMX are not needed for a transient run and must be omitted when the solution is transient. |
| • | When the first stress period is a steady-state stress period, SSMN is defined in record 3. For subsequent steady-state stress periods, SSMN is defined in record 9a. |
|
SSMX Maximum stage allowed for each lake in steady-state solution.
| • | SSMN and SSMX are not needed for a transient run and must be omitted when the solution is transient. |
| • | When the first stress period is a steady-state stress period, SSMN is defined in record 3. For subsequent steady-state stress periods, SSMN is defined in record 9a. |
|
|
|
9b.
|
If solute transport is also being simulated (Ftype “CONC” exists), then for each solute the following data are read:
CPPT(NSOL) CRNF(NSOL) [CAUG(NSOL)]
If data are being read using the fixed format mode, then each field should be entered using F10.4 format.
CPPT The concentration of solute in precipitation onto the lake surface.
|
CRNF The concentration of solute in overland runoff directly into the lake.
|
CAUG The concentration of solute in water used to augment the lake volume.
It is implicitly assumed that no solute is present in water evaporated from the lake surface.
|
|
|
