FOR EACH SIMULATION

Data Set 1

NSRF ISTRAT NOBS ISWIZT ISWIBD ISWIOBS [OPTIONS]
Number of active surfaces (interfaces). This equals the number of zones minus one.

Flag indicating density distribution.
•  0 – density varies linear between surfaces. 
•  1 – density is constant between surfaces. 

Number of OBS observation locations.

Flag and a unit number for ZETA output.
•  If ISWIZT > 0, unit number for ZETA output 
•  If ISWIZT ≤ 0, ZETA will not be recorded. 
If ISWIZT > 0, ZETA output is written with the same frequency as HEAD output specified in the Output Control (OC) file.

Flag and a unit number for BUDGET output. When this option is selected, corrections to the cell by cell flows computed by MODFLOW will be written to the same or different file (depending on the unit number). Corrections are called SWIADDTOFLF, SWIADDTOFRF, and SWIADDTOFFF, for the lower face (LF), right face (RF) and front face (FF), respectively
•  If ISWIBD > 0, unit number for BUDGET 
•  If ISWIBD ≤ 0, BUDGET will not be recorded. 

Flag and a unit number for OBS output
•  If ISWIOBS > 0, unit number for ASCII OBS 
•  If ISWIOBS = 0, OBS will not be recorded. 
•  If ISWIOBS < 0, ISWIOBS unit number for binary OBS. 

Is an optional list of character values
•  “ADAPTIVE” —defines that adaptive SWI2 time step variables will be specified in dataset 3b. 
•  "FSSSOPT" DEBUG OPTION  ADDITION DATA DURING TIP/TOE 



Data Set 2a

NSOLVER IPRSOL MUTSOL
Flag indicating solver used to solve for ZETA surfaces
•  If NSOLVER = 1, the MODFLOW DE4 solver will be used. 
•  If NSOLVER = 2, the MODFLOW PCG solver will be used. The PCG solver should be used for large problems when the time required for the DE4 solver to obtain a solution is excessive. 

is the printout interval for printing convergence information. If IPRSOL is equal to zero, it is changed to 999. The maximum ZETA change (positive or negative) and residual change are printed for each iteration of a time step whenever the time step is an even multiple of IPRSOL.

is a flag that controls printing of convergence information from the solver.
•  If MUTSOL = 0, tables of maximum head change and residual will be printed each iteration. 
•  If MUTSOL = 1, only the total number of iterations will be printed. 
•  If MUTSOL = 2, no information will be printed. 
•  If MUTSOL = 3, information will only be printed if convergence fails. 



IF NSOLVER = 2

Data Set 2b

MXITER ITER1 NPCOND ZCLOSE RCLOSE RELAX NBPOL DAMP [DAMPT]
Maximum number of outer iterations—that is, calls to the solution routine.

Maximum number of inner iterations—that is, iterations within the solution routine.

flag used to select the matrix conditioning method for the MODFLOW PCG solver (NSOLVER = 2).
•  If NPCOND = 1, is for Modified Incomplete Cholesky (for use on scalar computers). 
•  If NPCOND = 2, is for Polynomial (for use on vector computers or to conserve computer memory). 

is the ZETA change criterion for convergence, in units of length. When the maximum absolute value of ZETA change from all nodes during an iteration is less than or equal to ZCLOSE, and the criterion for RCLOSE is also satisfied (see below), iteration stops.

is the residual criterion for convergence, in units of cubic length per time. The units for length and time are the same as established for all model data. (See LENUNI and ITMUNI input variables in the Discretization File.) When the maximum absolute value of the residual at all nodes during an iteration is less than or equal to RCLOSE, and the criterion for ZCLOSE is also satisfied (see above), iteration stops.

is the relaxation parameter used with NPCOND = 1. Usually, RELAX = 1.0, but for some problems a value of 0.99, 0.98, or 0.97 will reduce the number of iterations required for convergence. RELAX is only used if NSOLVER is 2 and NPCOND is 1.

is only specified when NSOLVER = 2 and used when NPCOND = 2 to indicate whether the estimate of the upper bound on the maximum eigenvalue is 2.0, or whether the estimate will be calculated. NBPOL = 2 is used to specify the value is 2.0; for any other value of NBPOL, the estimate is calculated. Convergence is generally insensitive to this parameter. NBPOL is not used if NPCOND is not 2.

is the damping factor and is only specified when NSOLVER = 2. It is typically set equal to one, which indicates no damping. A value less than 1.0 and greater than 0.0 causes damping.
•  If DAMP > 0, applies to both steadystate and transient stress periods. 
•  If DAMP < 0, applies to steadystate periods. The absolute value if used as the dampening factor. 

is the damping factor for transient stress periods and is only specified when NSOLVER = 2. DAMPT is enclosed in brackets indicating that it is an optional variable that only is read when DAMP is specified as a negative value. If DAMPT is not read, then the single damping factor, DAMP, is used for both transient and steadystate stress periods.



Data Set 3a

TOESLOPE TIPSLOPE [ALPHA] [BETA]
Maximum slope of toe cells.

Maximum slope of tip cells.

Fraction of threshold used to move the tip and toe to adjacent empty cells when the slope exceeds userspecified TOESLOPE and TIPSLOPE values. If ALPHA is not specified, a default value of 0.1 is used. ALPHA must be greater than 0.0 and less than or equal to 1.0.

Fraction of threshold used to move the toe to adjacent nonempty cells when the surface is below a minimum value defined by the userspecified TOESLOPE value. A default BETA value of 0.1 is used if ALPHA is not specified. BETA must be greater than 0.0 and less than or equal to 1.0.



IF OPTIONS = ADAPTIVE

Data Set 3b

NADPTMX NADPTMN ADPTFCT
Maximum number of SWI2 time steps per MODFLOW time step. If NADPTMX is less than 1, it is changed to 1.

Minimum number of SWI2 time steps per MODFLOW time step. If NADPTMN is less than 1, it is changed to 1. NADPTMN must be less than or equal to NADPTMX.

is the factor used to evaluate tip and toe thicknesses and control the number of SWI2 time steps per MODFLOW time step. When the maximum tip or toe thickness exceeds the product of TOESLOPE or TIPSLOPE the cell size and ADPTFCT, the number of SWI2 time steps are increased to a value less than or equal to NADPT. When the maximum tip or toe thickness is less than the product of TOESLOPE or TIPSLOPE the cell size and ADPTFCT, the number of SWI2 time steps is decreased in the next MODFLOW time step to a value greater than or equal to 1. ADPTFCT must be greater than 0.0 and is reset to 1.0 if NADPTMX is equal to NADPTMN.



Data Set 4

NU(ISTRAT=0: NSRF+2, ISTRAT=1: NSRF+1)
Values of the dimensionless density
•  ISTRAT = 1 – Density of each zone (NSRF+1 values). 
•  ISTRAT = 0 – Density along top of layer, each surface, and bottom of layer (NSRF+2 values) 


FOR EACH SURFACE (FROM 1 to NSRF)

FOR EACH LAYER

Data Set 5

ZETA(NCOL,NROW) U2DREL
Initial elevations of the active surfaces.


FOR EACH LAYER

Data Set 6

SSZ(NCOL,NROW) U2DREL

FOR EACH LAYER

Data Set 7

ISOURCE(NCOL,NROW) U2DINT
Source type of any external sources or sinks, specified with any outside package (i.e. WEL Package, RCH Package, GHB Package). There are three options:
•  If ISOURCE > 0 – Sources and sinks have the same fluid density as the zone ISOURCE. If such a zone is not present in the cell, sources and sinks have the same fluid density as the active zone at the top of the aquifer. 
•  If ISOURCE = 0 – Sources and sinks have the same fluid density as the active zone at the top of the aquifer. 
•  If ISOURCE < 0 – Sources have the same fluid density as the zone with a number equal to the absolute value of ISOURCE. Sinks have the same fluid density as the active zone at the top of the aquifer. This option is useful for the modeling of the ocean bottom where infiltrating water is salt, yet exfiltrating water is of the same type as the water at the top of the aquifer. 


IF NOBS > 0

FOR EACH NOBS

Data Set 8

OBSNAM LAYER ROW COLUMN
is a string of 1 to 12 nonblank characters used to identify the observation. The identifier need not be unique; however, identification of observations in the output files is facilitated if each observation is given a unique OBSNAM.

is the layer index of the cell in which the ZETA observation is located.

is the row index of the cell in which the ZETA observation is located.

is the column index of the cell in which the ZETA observation is located.


