File | Description |
---|---|
02.ibound | IBOUND array used to define active model cells |
11.etsurf | Ground surface array used to calculate depth-to-water for ET |
12.bot | Aquifer bottom array used to calculate aquifer saturated thickness |
12.top | Top of saturation array used to calculate aquifer saturated thickness |
annual.lpf | Layer Property Flow package input file |
12p.k | Aquifer hydraulic conductivity array |
12.ss | Aquifer specific storage array |
12.sy | Aquifer specific yield array (not used) |
12p.pcg | Preconditioned Conjugate Gradient (PCG2) solver control parameters |
annual.drn | Drain package input file |
12.hyd | HYDMOD file for V12p7 streams |
12s.hyd | HYDMOD file for V12s streams |
annual.oc | Output Control |
It should never be necessary to modify any of these files.
The extent of the active aquifer is defined by 02.ibound. This file remains unchanged from model version 2.
The aquifer is simulated as having a constant saturated thickness, calculated as the difference between the top and bottom of the saturated aquifer specified in the files 12.top and 12.bot, respectively. The aquifer hydraulic conductivity is specified in the file 12p.k in units of feet per day, which the LPF package then uses to calculate an aquifer transmissivity.
The aquifer is in reality a phreatic aquifer. However, in order to simulate the aquifer as having a constant transmissivity, the LPF package treats the aquifer as being confined. The specific yield of the aquifer is contained in the file 12.sy. However, the LPF package requires a specific storage, defined as specific yield divided by aquifer thickness. These values are contained in the file 12.ss and it is this file that is read by the LPF package.
Note that the top of saturation 12.top and ET ground surface 11.surf are quite different. The ET ground surface array represents the actual ground surface, while the top of saturation array represents an average water table elevation.
Except for changes to the LPF package file necessitated by changes to the directory structure, all the aquifer property files remain unchanged from V12p7.
The files 12.top and 12.bot are referenced from the MODFLOW-2000 discretization package. However, since this package also contains the time step information which differs year to year because of leap years, the discretization package file is generated for each year.
The V12p7 model consisted of a steady state pre-development period, and a transient period from 1918 to 2000. The drain package contained a different set of data for the pre-development steady state and the transient periods. The file annual.drn contains the V12p7 drain parameters as well as a repeat flag for 1199 stress periods so that this file can be used in simulations for up to 100 years of monthly stress periods.
The 12.hyd file contains the set of locations corresponding the the V12p7 stream package. An error in the V12p7 stream package was discovered after the model was officially adopted by the RRCA. A corrected model, known as V12s, was adopted for use starting with model year 2001. The correction in the streams required a different set of stream locations, which are contained in 12s.hyd.
For model simulations from 2001 onwards, the 12s.hyd should be used. For the year 2000, the 12.hyd was used to verify that model simulations exactly match those obtained during model calibration, while 12s.hyd was used to evaluate the differences introduced by the V12s stream correction.
In order to provide for a stable solution algorithm, a damping factor of 0.8 is used. This allows most numerical oscillations to be overcome without excessively slowing down the convergence.
A maximum of 100 outer iterations are permitted. Occasionally, this maximum number of iterations will be reached without satisfying the head and residual criteria, but with a very good mass balance, usually better than 0.01 percent. The version of MODFLOW-2000 used has the Mass Balance Override feature enabled to continue the simulation when this occurs.