FOR EACH SIMULATION
|
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 (199 characters) that starts in column 2. Any characters can be included in Text. The “#” character must be in column 1. Except for the name file, lines beginning with # are restricted to these first lines of the file. Text is printed when the file is read.
|
|
Data Set 1
|
NPHFB MXFB NHFBNP [MXHFBNP] Option
MXFB—is the maximum number of HFB barriers that will be defined using parameters.
|
NHFBNP—is the number of HFB barriers not defined by parameters for current stress period. Item 4 will be read NHFBNP times as specified in item 1. If NHFBNP is set to -1 in item 7, then NACTHFB and items 8 and 9 are not read, and the barriers defined in the previous stress period are used. Item 8 will be read NHFBNP times as specified in item 7.
|
the maximum number of HFB barriers not defined by parameters used during any stress period. If MXHFBNP is not specified, it is set equal to NHFBNP in item 1. MXHFBNP is only included in MODFLOW-OWHM.
|
Option—is an optional list of character values.
•NOPRINT—suppresses the printing of cell lists. |
|
|
Repeat Items 2 and 3 for each NPHFB parameter. Items 2 and 3 are not read if NPHFB is negative or zero.
|
Data Set 2
|
PARNAM PARTYP Parval NLST
PARNAM—is the name of a parameter. This name can consist of 1 to 10 characters and is not case sensitive. That is, any combination of the same characters with different case will be equivalent. All parameter names must be unique.
|
PARTYP—is the type of parameter. For the HFB Package, the only allowed parameter type is HFB, which defines values of the hydraulic characteristic of the barrier.
|
NLST—is the number of horizontal-flow barrier cells included in the parameter.
|
|
|
Data Set 3
|
Layer IROW1 ICOL1 IROW2 ICOL2 Factor [Layer2]
NLST repetitions of Item 3 are required; they are read by module ULSTRD. (SFAC of the ULSTRD utility module applies to Factor).
Layer—is the number of the model layer in which the horizontal-flow barrier is located.
|
IROW1—is the row number of the cell on one side of the horizontal-flow barrier.
|
ICOL1—is the column number of the cell on one side of the horizontal-flow barrier.
|
IROW2—is the row number of the cell on the other side of the horizontal-flow barrier.
|
ICOL2—is the column number of the cell on the other side of the horizontal-flow barrier.
|
Factor—is the factor used to calculate hydraulic characteristic from the parameter value. The hydraulic characteristic is the product of Factor and the parameter value. If final hydraulic characteristic is negative, then it acts as a multiplier to the conductance between the two model cells specified as containing a barrier. For example, if the final result was –1.5, it would multiply by 1.5 the conductance between the two cells. Negative Factor values are only supported in MODFLOW-OWHM.
|
the model layer of the second cell which contains a barrier. This cell is indexed as [IROW2, ICOL2, Layer2]. This variable is optional and if not present is set equal to Layer. If it is equal to Layer, then no flow routing occurs, and the barrier functions identically to previous versions of HFB. To use the flow routing capability, the NWT solver must be used. Layer2 is only included in MODFLOW-OWHM.
If Layer2 is not equal to Layer, either IROW1 must equal IROW2 or ICOL1 must equal ICOL2. Typically either the absolute value of IROW1 - IROW2 equals 1 or the absolute value of ICOL1 - ICOL2 equals 1
For vertical barriers (IROW1 equals IROW2 and ICOL1 equals ICOL2), Layer and Layer2 must specify adjacent layers.
|
|
If IROW1=IROW2 and ICOL1=ICOL2 and ABS(Layer2- Layer)=1, then there is no flow routing, but instead the barrier is changed to a vertical barrier between layers and will work with any solver package.
|
FOR FIRST STRESS PERIOD OR ALL STRESS PERIODS
|
Data Set 4
|
Layer IROW1 ICOL1 IROW2 ICOL2 Hydchr [Layer2]
NHFBNP repetitions of Item 4 are read. Item 4 is not read if NHFBNP is negative or zero.
Layer—is the number of the model layer in which the horizontal-flow barrier is located.
|
IROW1—is the row number of the cell on one side of the horizontal-flow barrier.
|
ICOL1—is the column number of the cell on one side of the horizontal-flow barrier.
|
IROW2—is the row number of the cell on the other side of the horizontal-flow barrier.
|
ICOL2—is the column number of the cell on the other side of the horizontal-flow barrier.
|
Hydchr—is the hydraulic characteristic of the horizontal-flow barrier. The hydraulic characteristic is the barrier hydraulic conductivity divided by the width of the horizontal-flow barrier. If hydraulic characteristic is negative, then it acts as a multiplier to the conductance between the two model cells specified as containing a barrier. For example, if the final result was –1.5, it would multiply by 1.5 the conductance between the two cells. Negative Hydchr values are only supported in MODFLOW-OWHM.
|
the model layer of the second cell which contains a barrier. This cell is indexed as [IROW2, ICOL2, Layer2]. This variable is optional and if not present is set equal to Layer. If it is equal to Layer, then no flow routing occurs, and the barrier functions identically to previous versions of HFB. To use the flow routing capability, the NWT solver must be used. Layer2 is only included in MODFLOW-OWHM.
If Layer2 is not equal to Layer, either IROW1 must equal IROW2 or ICOL1 must equal ICOL2. Typically either the absolute value of IROW1 - IROW2 equals 1 or the absolute value of ICOL1 - ICOL2 equals 1
For vertical barriers (IROW1 equals IROW2 and ICOL1 equals ICOL2), Layer and Layer2 must specify adjacent layers.
|
|
If IROW1=IROW2 and ICOL1=ICOL2 and ABS(Layer2- Layer)=1, then there is no flow routing, but instead the barrier is changed to a vertical barrier between layers and will work with any solver package.
|
Data Set 5
|
NACTHFB—is the number of active HFB parameters.
|
|
Data Set 6
|
Pname
NACTHFB repetitions of Item 6 are read. Item 6 is not read if NACTHFB is negative or zero.
Pname—is the name of a parameter to be used in the simulation. NACTHFB parameter names will be read.
|
|
OPTIONAL: FOR EACH SUBSEQUENT STRESS PERIOD Data Sets 7-9 are only included in MODFLOW-OWHM
|
Data Set 7 (only in MODFLOW-OWHM)
|
NACTHFB NHFBNP
If NACTHFB is set to –1, then NHFBNP and items 8 and 9 are not read. Instead the barriers defined in the previous stress period are used. If MXHFBP is 0, then NHFBNP is not read in since it cannot be larger than 0.
the number of active HFB parameters for the current stress period.
|
the number of HFB barriers not defined by parameters for current stress period. Item 4 will be read NHFBNP times as specified in item 1. If NHFBNP is set to -1 in item 7, then NACTHFB and items 8 and 9 are not read, and the barriers defined in the previous stress period are used. Item 8 will be read NHFBNP times as specified in item 7.
|
|
|
Data Set 8 (only in MODFLOW-OWHM)
|
Pname
NACTHFB, from item 7, repetitions of Item 8 are read. Item 8 is not read if NACTHFB is negative or zero.
the name of a parameter to be used in the current stress period. NACTHFB parameter names will be read.
|
|
Data Set 9 (only in MODFLOW-OWHM)
|
Layer IROW1 ICOL1 IROW2 ICOL2 Hydchr [Layer2]
NHFBNP, from item 7, repetitions of Item 9 are read. Item 9 is not read if NHFBNP or MXHFBP is negative or zero.
Layer—is the number of the model layer in which the horizontal-flow barrier is located.
|
IROW1—is the row number of the cell on one side of the horizontal-flow barrier.
|
ICOL1—is the column number of the cell on one side of the horizontal-flow barrier.
|
IROW2—is the row number of the cell on the other side of the horizontal-flow barrier.
|
ICOL2—is the column number of the cell on the other side of the horizontal-flow barrier.
|
Hydchr—is the hydraulic characteristic of the horizontal-flow barrier. The hydraulic characteristic is the barrier hydraulic conductivity divided by the width of the horizontal-flow barrier. If hydraulic characteristic is negative, then it acts as a mulitplier to the conductance between the two model cells specified as containing a barrier. For example, if the final result was –1.5, it would multiply by 1.5 the conductance between the two cells.
|
the model layer of the second cell which contains a barrier. This cell is indexed as [IROW2, ICOL2, Layer2]. This variable is optional and if not present is set equal to Layer. If it is equal to Layer, then no flow routing occurs, and the barrier functions identically to previous versions of HFB. To use the flow routing capability, the NWT solver must be used. Layer2 is only included in MODFLOW-OWHM.
If Layer2 is not equal to Layer, either IROW1 must equal IROW2 or ICOL1 must equal ICOL2. Typically either the absolute value of IROW1 - IROW2 equals 1 or the absolute value of ICOL1 - ICOL2 equals 1
For vertical barriers (IROW1 equals IROW2 and ICOL1 equals ICOL2), Layer and Layer2 must specify adjacent layers.
|
If IROW1=IROW2 and ICOL1=ICOL2 and ABS(Layer2- Layer)=1, then there is no flow routing, but instead the barrier is changed to a vertical barrier between layers and will work with any solver package.
|
|
