FOR EACH SIMULATION

Data Set 1

IBCFCB HDRY IWDFLG WETFCT IWETIT IHDWET
These 6 variables are free format if the option “FREE” is specified in the Basic Package input file; otherwise, the variables all have 10character fields.
IBCFCB—is a flag and a unit number.
•If IBCFCB > 0, it is the unit number to which cellbycell flow terms will be written when "SAVE BUDGET" or a nonzero value for ICBCFL is specified in Output Control. The terms that are saved are storage, constanthead flow, and flow between adjacent cells.
•If IBCFCB = 0, cellbycell flow terms will not be written.
•If IBCFCB < 0, cellbycell flow for constanthead cells will be written in the listing file when "SAVE BUDGET" or a nonzero value for ICBCFL is specified in Output Control. Cellbycell flow to storage and between adjacent cells will not be written to any file.
The flow terms that will be saved are the flows through the right, front, and lower cell face. Positive values represent flows toward higher column, row, or layer numbers.

HDRY—is the head that is assigned to cells that are converted to dry during a simulation. Although this value plays no role in the model calculations, it is useful as an indicator when looking at the resulting heads that are output from the model. HDRY is thus similar to HNOFLO in the Basic Package, which is the value

IWDFLG—is a flag that determines if the wetting capability is active.
•If IWDFLG = 0, the wetting capability is inactive.
•If IWDFLG is not 0, the wetting capability is active. 
WETFCT—is a factor that is included in the calculation of the head that is initially established at a cell when it is converted from dry to wet. (See IHDWET.)

IWETIT—is the iteration interval for attempting to wet cells. Wetting is attempted every IWETIT iteration. If using the PCG solver (Hill, 1990), this applies to outer iterations, not inner iterations. If IWETIT is less than or equal to 0, it is changed to 1.

IHDWET—is a flag that determines which equation is used to define the initial head at cells that become wet:
•If IHDWET = 0, equation 33A (equation 3a in McDonald and others, 1992) is used:
h = BOT + WETFCT (hn  BOT)
•If IHDWET is not 0, equation 33B (equation 3b in McDonald and others, 1992) is used:
h = BOT + WETFCT (WETDRY)



Data Set 2

Ltype(NLAY)
Read one value for each layer. These values are free format if the word FREE is specified in Item 1 of the Basic Package input file; otherwise, they have fixedformat fields. The fixed format fields are each 2 characters wide with 40 values per line. Use only as many lines as required for the number of model layers.
In MODFLOW2005 v 1.5 and earlier versions of MODFLOW2005, Ltype must always be in free format.
Ltype—contains a combined code for each layer that specifies both the layer type (LAYCON) and the method of computing interblock conductance. Use as many records as needed to enter a value for each layer. Values are twodigit numbers:
The left digit defines the method of calculating interblock transmissivity. The methods are described by Goode and Appel (1992).
•0 or blank—harmonic mean (the method used in MODFLOW88). (This is most appropriate for confined and unconfined aquifers with abrupt boundaries in transmissivity at the cell boundaries or for confined aquifers with uniform hydraulic conductivity.)
•1—arithmetic mean (This is most appropriate for unconfined aquifers with uniform hydraulic conductivities.)
•2—logarithmic mean (This is most appropriate for confined aquifers with gradually varying transmissivities.)
•3—arithmetic mean of saturated thickness and logarithmicmean hydraulic conductivity. (This is most appropriate for unconfined aquifers with gradually varying transmissivities.) If the transmissivity method is 3 and LAYCON (see below) is 0 or 2, the transmissivity method will be changed to 2 (logarithmic).
The right digit defines the layer type (LAYCON), which is the same as in MODFLOW88:
•0—confined—Transmissivity and storage coefficient of the layer are constant for the entire simulation.
•1—unconfined—Transmissivity of the layer varies. It is calculated from the saturated thickness and hydraulic conductivity. The storage coefficient is constant. This type code is valid only for layer 1.
•2—confined/unconfined—Transmissivity of the layer is constant. The storage coefficient may alternate between confined and unconfined values. Vertical flow from above is limited if the layer desaturates.
•3—confined/unconfined—Transmissivity of the layer varies. It is calculated from the saturated thickness and hydraulic conductivity. The storage coefficient may alternate between confined and unconfined values. Vertical flow from above is limited if the aquifer desaturates. 

Data Set 3

TRPY(NLAY)  U1DREL
TRPY—is a onedimensional variable containing a horizontal anisotropy factor for each layer. It is the ratio of transmissivity or hydraulic conductivity (whichever is being used) along a column to transmissivity or hydraulic conductivity along a row. Set to 1.0 for isotropic conditions. This is a single variable with one value per layer. Do not read a variable for each layer—that is, include only one array control record for the entire variable. Note that the LPF Package can be used to represent horizontal anisotropy that varies in magnitude within a model layer.


A subset of the following twodimensional variables are used to describe each layer. The variables needed for each layer depend on the layertype code (LAYCON, which is defined as part of the Item2 Ltype), whether the simulation has any transient stress periods (at least one stress period defined in the Discretization File specifies Ss/Tr as “TR”), and if the wetting capability is active (IWDFLG not 0). If a variable is not needed, it must be omitted. In no situation will all variables be required. The required variables (Items 49) for layer 1 are read first; then the variables for layer 2 and so forth.

Data Set 4

[Sf1(NCOL,NROW)]U2DREL
If there is at least one transient stress period.
Sf1—is the primary storage coefficient. Read only if there are one or more transient stress periods specified in the Discretization File. For LAYCON equal to 1, Sf1 will always be specific yield, whereas for LAYCON equal to 2 or 3, Sf1 will always be confined storage coefficient. For LAYCON equal to 0, Sf1 would normally be confined storage coefficient; however, a LAYCON value of 0 can also be used to simulate watertable conditions where drawdowns are expected to remain everywhere a small fraction of the saturated thickness, and where there is no layer above, or flow from above is negligible. In this case, specific yield values would be entered for Sf1. (The confined storage coefficient is specific storage times the layer thickness.)


Data Set 5

[Tran(NCOL,NROW)]  U2DREL
If LAYCON is 0 or 2 (see Ltype ), then read item 5.
Tran—is the transmissivity along rows. Tran is multiplied by TRPY to obtain transmissivity along columns. Read only for layers where LAYCON is 0 or 2.


Data Set 6

[HY(NCOL,NROW)]  U2DREL
Otherwise, if LAYCON is 1 or 3 (see Ltype), read item 6.
HY—is the hydraulic conductivity along rows. HY is multiplied by TRPY to obtain hydraulic conductivity along columns. Read only for layers where LAYCON is 1 or 3.


Data Set 7

[Vcont(NCOL,NROW)]  U2DREL
If not the bottom layer.
Vcont—is the vertical hydraulic conductivity divided by the thickness from a layer to the layer below. The value for a cell is the hydraulic conductivity divided by thickness for the material between the node in that cell and the node in the cell below. Because there is not a layer beneath the bottom layer, Vcont cannot be specified for the bottom layer. See p 516 to 517 of the MODFLOW2005 documentation for more details.


Data Set 8

[Sf2(NCOL,NROW)]  U2DREL
If there is at least one transient stress period and LAYCON (see Ltype) is 2 or 3.
Sf2—is the secondary storage coefficient. Read only for layers where LAYCON is 2 or 3 and only if there are one or more transient stress periods specified in the Discretization File. The secondary storage coefficient is always specific yield.


Data Set 9

[WETDRY(NCOL,NROW)]  U2DREL
If IWDFLG is not 0 and LAYCON is 1 or 3 (see Ltype).
WETDRY—is a combination of the wetting threshold and a flag to indicate which neighboring cells can cause a cell to become wet.
•If WETDRY < 0, only the cell below a dry cell can cause the cell to become wet.
•If WETDRY > 0, the cell below a dry cell and the four horizontally adjacent cells can cause a cell to become wet.
•If WETDRY is 0, the cell cannot be wetted.
The absolute value of WETDRY is the wetting threshold. When the sum of BOT and the absolute value of WETDRY at a dry cell is equaled or exceeded by the head at an adjacent cell, the cell is wetted. Read only if LAYCON is 1 or 3 and IWDFLG is not 0.
The Lake package requires that WETDRY be set to 0 in Lakes.

