Lakes and streams can be designated as having gaging stations located on them. For each such designated lake or stream, the time and a variety of other information will be written to a separate output file after each time step (and each transport time increment) to facilitate graphical postprocessing of the calculated data.
Cells of the model grid can be designated as “stream gaging station” locations. At each designated cell (or stream reach), the time, stage, streamflow out of that reach, the streambed seepage, unsaturated storage, change in unsaturated storage, and groundwater recharge will be written to a separate output file to facilitate model output evaluation and graphical post processing of the calculated data. Several options are available to print additional information. The input file for specifying gaging station locations is read if the file type (Ftype) “GAGE” is included in the MODFLOW name file. The output file will contain two header lines that provide relevant information (the text will be contained within quotes).
Data Set 2 must include exactly NUMGAGE lines (or records) of data. If NUMGAGE > 1, it is permissible to interleaf in Data Set 2 records for stream gaging stations (according to the format specified in the documentation for the Stream Package) with records for gages on lakes. Data lines (records) within Data Set 2 can be listed in any arbitrary order.
LAKE Negative value of the lake number of the lake where the gage is located.
•
The Gage Package can also be used in conjunction with the Stream Flow Routing (SFR) Package to specify the location of a gaging station on a stream. Therefore, to guarantee that the code can distinguish between input for lakes and that for stream locations, lake numbers are specified as their negative value.
UNIT Unit number for output file. Specify negative value of unit number to invoke reading of optional parameter OUTTYPE for listing additional information for lakes.
•
The user should specify a unique unit number for each gaging station and match those unit numbers to DATA file types and file names in the MODFLOW name file (Harbaugh and McDonald, 1996b).
OUTTYPE Code for type of expanded listing desired in output file:
0.
Use standard default listing of time, stage, volume, and concentration.
1.
•
Default values plus all inflows to and outflows from lake (as volumes during time increment), and total lake conductance.
2.
•
Default values plus changes in lake stage, lake volume, and solute concentrations.
•
3.
All of the above except time-step budget error.
Total lake conductance (OUTTYPE options 1 and 3 for a lake gaging station) is the sum of the conductances of each seepage interface for each lake. Changes in lake stage, volume, and solute concentrations (OUTTYPE options 2 or 3) are listed as incremental changes from previous time increment and as cumulative change since start of simulation.
Note 1: A unique unit number must be specified for each stream gaging station and the MODFLOW name file (see Harbaugh and others, 2000, p. 42-44) must specify the file for each unit. The file type in the name file must be “DATA”.
Note 2: If the LAK3 Package is also active, a gaging station may be placed on a lake. In this case, only two required variables and one optional variable are read. The first variable should be the negative value of the lake number and the second is the unit number for the output. The original description of the Gage Package (Merritt and Konikow, 2000, p. 57) did not include additional print options listed for OUTTYPE. For a gaging station on a lake, a minus sign on the unit number is a flag indicating that OUTTYPE for a lake will be read.
Note 3: Data Set 2 must include exactly NUMGAGE lines (records) of data. If NUMGAGE > 1, it is permissible to interleaf the Item 2 lines for stream gaging stations with lines for lake gages. Data lines (records) within Item 2 can be listed in any arbitrary order.
Note 5: Unsaturated flow is not simulated beneath lakes.
OUTTYPE An integer value that is a flag for type of expanded listing desired in output file:
•
0 Use standard default listing of time, stage, outflow, and solute concentration.
•
1 Default values plus depth, width, and flow at midpoint.
•
2 Default values plus streambed conductance for the reach, head difference across streambed, and hydraulic gradient across streambed.
•
3 Default values plus solute load in stream (if GWT is active).
•
4 All of the above.
•
5 Use for diversions to provide a listing of time, stage, flow diverted, maximum assigned diversion rate, flow at end of upstream segment prior to diversion, solute concentration, and solute load.
•
6 Used for unsaturated flow routing to provide a listing of time, stage, groundwater head, streambed seepage, change in unsaturated storage, and recharge.
•
7 Used for unsaturated flow routing to provide a listing of time and the unsaturated water content profile beneath the stream. Two profiles are printed. The first is the volume average water content of ISUZN unsaturated cells beneath a stream reach. The second profile is the volume averaged water content beneath the low flow channel only. If ISUZN =1 then both water content profiles will be the same.
Note 1: A unique unit number must be specified for each stream gaging station and the MODFLOW name file (see Harbaugh and others, 2000, p. 42-44) must specify the file for each unit. The file type in the name file must be “DATA”.
Note 3: Data Set 2 must include exactly NUMGAGE lines (records) of data. If NUMGAGE > 1, it is permissible to interleaf the Item 2 lines for stream gaging stations with lines for lake gages. Data lines (records) within Item 2 can be listed in any arbitrary order.
Note 4: Solute is not routed through the unsaturated zone. For this reason, solute routing between streams and ground water cannot be simulated when using the unsaturated flow option in SFR2.
