Data Set 0

[#Text]
Item 0 is optional  “#” must be in column 1. Item 0 can be repeated multiple times.
Text  is a character variable (up to 199 characters) that starts in column 2. Any characters can be included in Text. The “#” character must be in column 1. Text is printed when the file is read.


Data Set 1

NPTH NLOC IOUTADV KTFLG KTREV ADVSTP FSNK [IMPATHOUT IMPATHOFF]
NPTH  is the total number of advectivetransport paths that will be used (that is, the total number of repetitions of ITEM 3).

NLOC  is the total number of advectivetransport observation locations that will be used (that is, the total number of repetitions ITEM 4).

IOUTADV – is a flag and a unit number.
•  If IOUTADV > 0, it is the unit number to which the particletracking information will be written. 
•  If IOUTADV = 0, particletracking information will not be written. 

KTFLG – is a flag indicating how the particletracking time step is calculated..
•  If KTFLG = 1, particles are displaced from one cell face to the next, and the timestep length varies. 
•  If KTFLG = 2, particles are displaced using the time steps specified by ADVSTP. 
•  If KTFLG = 3, both of the above are included. 

KTREV – is a flag indicating the direction of particle displacement.
•  If KTREV = 1, particles are displaced in a forward direction. 
•  If KTREV = 1, particles are displaced in a backward direction. 

ADVSTP – is the particletracking timestep length to be used when KTFLG equals 2 or 3.

FSNK – is a flag and a fraction indicating how weak sinks are treated. Performs identically to Option 3 for terminating pathlines in MODPATH (Pollock, 1989, p. 34).
•  If FSNK < 0, particles will be discharged at cells with any amount of discharge to boundary conditions. 
•  If FSNK > 0, particles will be discharged at cells where discharge to sinks is larger than the specified fraction of the total inflow to the cell. 

IMPATHOUT – is a flag and a unit number.
•  If IMPATHOUT > 0, it is the unit number to which particle tracking information will be written using the conventions of the MODPATH pathline file. 
•  If IMPATHOUT = 0, no file is produced. 

IMPATHOFF – is a flag indicating the convention used to interpret the offsets in items 3 and 4 (see Anderman and Hill, 2001, p. 49). The offsets for particle starting locations are defined in input item 3 by variables SLOFF, SROFF, and SCOFF. The offsets for observation locations are defined in input item 4 by variables LOFF, ROFF, and COFF. If IMPATHOFF is omitted or set to zero, the offsets in the ADV2 input file need to be defined using the ADV2 convention. If IMPATHOFF equals 1, the offsets need to be defined using the MODPATH convention.



If porosity is to be defined using parameters, Item 2 must be entered as shown below and followed by Items 2a. and 2b.
The capability to define some porosity values using parameters and others without using parameters is not provided. Thus, if parameters are used to define the porosity, then they must define porosity for all model cells and Quasi3D confining beds. Conversely, if porosity is defined using layer arrays, then no porosity parameters can be defined.

Data Set 2 (with parameters)

PARAMETER NPADV IPFLG
IPFLG – is the print code for all ADV parameters, which controls the printing of the porosity arrays created when the parameters are substituted into the arrays.



Repeat Items 2a2b for each parameter to be defined (that is, NPADV times).

Data Set 2a

[PARNAM PARTYP Parval NCLU]
PARNAM – is the name of a parameter to be defined. This name can consist of up to 10 characters and is not case sensitive.
All parameter names must be unique.

PARTYP – is the type of parameter to be defined. For the ADV2 Package, the supported parameter types are:
•  PRST – defines the effective porosity of the material represented by a model layer. 
•  PRCB – defines the effective porosity of the material represented by a quasi3D confining unit. 

NCLU – is the number of clusters required to define the parameter. Each repetition of item 2b is called a cluster, with variables Layer, Mltarr, Zonarr, and IZ.



Data Set 2b

[Layer Mltarr Zonarr IZ]
Each repetition of Item 2b is called a parameter cluster. Repeat Item 2b NCLU times.
Layer – is the modellayer number to which a cluster definition applies. Even if the HydrogeologicUnit Flow package (Anderman and Hill, 2001) is used, the effective porosity is defined by model layer, not hydrogeologic unit.

Mltarr – is the name of the multiplier array to be used to define array values that are associated with a parameter. The name “NONE” means that there is no multiplier array, and the array values will be set equal to Parval.

Zonarr – is the name of the zone array to be used to define array elements that are associated with a parameter. The name “ALL” means that there is no zone array and that all elements in the model layer are part of the parameter.

IZ – is up to 10 zone numbers (separated by spaces) that define the array elements that are associated with a parameter. These values are not used if Zonarr is specified as “ALL”. Values can be positive or negative, but 0 is not allowed. The first zero or nonnumeric value terminates the list.



If porosity is not to be defined using parameters, Item 2 must be entered as shown

Data Set 2 (without parameters)

PRST(NCOL,NROW)
Item 2 is read using the arrayreading utility module U2DREL (Harbaugh and others, 2000, p. 8688. Item 2 is repeated for each model layer and quasi3D confining unit in the grid. Thus, the number of PRST arrays must equal the number of model layers plus the number of quasi3D confining beds. These layer variables are read in sequence going down from the top of the system.
PRST – is the porosity of the porous medium represented by a model layer or quasi3D confining unit.


Items 34 are repeated as a group for each advectivetransport path to be defined (that is, NPTH times).

Data Set 3

NPNT SLAY SROW SCOL SLOFF SROFF SCOFF
NPNT – is the number of observation points along the path being defined.

SLAY – is the layer of the initial location.

SROW – is the row of the initial location.

SCOL – is the column of the initial location.

SLOFF – is the layer offset used to locate the initial location within the finitedifference cell (must range between –0.5 at the bottom of the layer and 0.5 at the top of the layer).

SROFF – is the row offset used to locate the initial location within the finitedifference cell (must range between –0.5 and 0.5; see fig. 2 of Hill and others, 2000).

SCOFF – is the column offset used to locate the initial location within the finitedifference cell (must range between –0.5 and 0.5; see fig. 2 of Hill and others, 2000).



Data Set 4

OBSNAM LAY ROW COL LOFF ROFF COFF XSTAT IXSTAT YSTAT IYSTAT ZSTAT IZSTAT TADV PLOTSYMBOL
Item 4 is repeated for each advectivetransport observation point to be defined along the path (that is, NPNT times).
OBSNAM – 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. The current version of MODFLOW2000 (v 1.17.02) accepts duplicate observation names, but if duplicate names are found, a warning concerning the duplication is written to the Global file.

LAY – is the layer of the advectivetransport observation.

ROW – is the row of the advectivetransport observation.

COL – is the column of the advectivetransport observation.

LOFF – is the layer offset used to locate the advectivetransport observation within the finitedifference cell (must range between –0.5 at the bottom of the layer and 0.5 at the top of the layer).

ROFF – is the row offset used to locate the advectivetransport observation within the finitedifference cell (must range between –0.5 and 0.5; see fig. 2 of Hill and others, 2000).

COFF – is the column offset used to locate the advectivetransport observation within the finitedifference cell (must range between –0.5 and 0.5; see fig. 2 of Hill and others, 2000).

XSTAT, YSTAT, ZSTAT – are values from which the weights of the observed locations are calculated depending on how IXSTAT, IYSTAT, and IZSTAT are set.

IXSTAT, IYSTAT, IZSTAT – are flags to indicate what XSTAT, YSTAT, and ZSTAT are and how the observation weights are calculated. If the variable STATISTIC represents XSTAT, YSTAT, or ZSTAT and the variable STATFLAG represents IXSTAT, IYSTAT, or IZSTAT, respectively, the following definitions apply:
•  STATFLAG = 0, STATISTIC is a scaled variance [L2], weight = 1/STATISTIC, 
•  STATFLAG = 1, STATISTIC is a scaled standard deviation [L], weight = 1/STATISTIC2, and 
•  STATFLAG = 2, STATISTIC is a scaled coefficient of variation [dimensionless], weight = 1/(STATISTIC × Observation value)2. The observation value is calculated as the distance in the x, y, or zcoordinate directions between particle starting position, as defined in input item 3, and the observation positions, as defined in item 4. 

XSTAT, YSTAT, ZSTAT – are values from which the weights of the observed locations are calculated depending on how IXSTAT, IYSTAT, and IZSTAT are set.

IXSTAT, IYSTAT, IZSTAT – are flags to indicate what XSTAT, YSTAT, and ZSTAT are and how the observation weights are calculated. If the variable STATISTIC represents XSTAT, YSTAT, or ZSTAT and the variable STATFLAG represents IXSTAT, IYSTAT, or IZSTAT, respectively, the following definitions apply:
•  STATFLAG = 0, STATISTIC is a scaled variance [L2], weight = 1/STATISTIC, 
•  STATFLAG = 1, STATISTIC is a scaled standard deviation [L], weight = 1/STATISTIC2, and 
•  STATFLAG = 2, STATISTIC is a scaled coefficient of variation [dimensionless], weight = 1/(STATISTIC × Observation value)2. The observation value is calculated as the distance in the x, y, or zcoordinate directions between particle starting position, as defined in input item 3, and the observation positions, as defined in item 4. 

XSTAT, YSTAT, ZSTAT – are values from which the weights of the observed locations are calculated depending on how IXSTAT, IYSTAT, and IZSTAT are set.

XSTAT, YSTAT, ZSTAT – are values from which the weights of the observed locations are calculated depending on how IXSTAT, IYSTAT, and IZSTAT are set.

TADV – is the time of this advectivefront observation, relative to the beginning of the steadystate simulation, in consistent units.

PLOTSYMBOL – is an integer that is written to output files intended for graphical analysis to allow control of the symbols used to plot data.



Data Set 5

IOWTQAD – is a flag indicating whether a full weight matrix will be specified.
•  If IOWTQAD ≤0, a full weight matrix will not be specified and the STATISTIC values in ITEM 4 will be used to calculate weights. 
•  If IOWTQAD > 0, a full weight matrix will be specified. 


Items 67 are required if IOWTQAD>0.

Data Set 6

FMTIN IPRN
FMTIN—is the Fortran format to be used in reading each line of the variancecovariance matrix used to calculate the weighting. The format needs to be enclosed in parentheses and needs to accommodate real numbers.

IPRN—is a flag identifying the format in which the variancecovariance matrix is printed. If IPRN is less than zero, the matrix is not printed. Permissible values of IPRN and corresponding formats are:






1

10G12.3

6

5G12.3

2

10G12.4

7

5G12.4

3

9G12.5

8

5G12.5

4

8G13.6

9

4G13.6

5

8G14.7

10

4G14.7




Data Set 7

WTQ(NTT2,NTT2) (format FMTIN)
WTQ—is an array containing the variancecovariance matrix on advectivetransport observations [L 2 ]. The array is square with dimensions equal to the total number of advectivetransport observations (for example, NLOC times two if the tracking is two dimensional or three if the tracking is three dimensional) and for each observation location is entered in the order X, Y, and Z, if applicable. For elements WTQ(I,J), if I ≠J, WTQ(I,J) is the covariance between observations I and J; if I = J, WTQ(I,J) is the variance of observation I. Note that the variancecovariance matrix is symmetric, but the entire matrix (upper and lower parts) must be entered.

