RIV - River Package

Top  Previous  Next  Print

Abbreviation in Name file

RIV

Purpose

The River package is used to simulate head-dependent flux boundaries.  In the River package if the head in the cell falls below a certain threshold, the flux from the river to the model cell is set to a specified lower bound.

Documentation

Related Packages

Supported in

MODFLOW-2000
MODFLOW-2005

Other Notes

Special care must be taken when constructing the RIV package input file if the observations are used.  See discussion of "Observations at Cells Having More Than One Head-Dependent Boundary Feature Represented by the Same Package" in the RVOBinput file.

Input Instructions

Input to the River (RIV) Package is read from the file that has file type "RIV" in the name file. Optional variables are shown in brackets. All variables are free format if the option “FREE” is specified in the Basic Package input file; otherwise, the non-optional variables have 10-character fields and the optional variables are free format.

FOR EACH SIMULATION

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.

1.

[ PARAMETER NPRIV MXL]

This optional item must start with the word “PARAMETER”.

MXL—is the maximum number of river reaches that will be defined using parameters. MXL must equal or exceed the sum of NLST×N for all parameters, where N is the greater of 1 and NUMINST.

2.

MXACTR IRIVCB [Option]

MXACTR—is the maximum number of river reaches in use during any stress period. MXACTR includes reaches that are defined using parameters as well as reaches that are defined without using parameters.

IRIVCB—is a flag and a unit number.

If IRIVCB > 0, it is the unit number to which cell-by-cell flow terms will be written when "SAVE BUDGET" or a non-zero value for ICBCFL is specified in Output Control.
If IRIVCB = 0, cell-by-cell flow terms will not be written.
If IRIVCB < 0, river leakage for each reach will be written to the listing file when "SAVE BUDGET" or a non-zero value for ICBCFL is specified in Output Control.

Option—is an optional list of character values.

AUXILIARY abc” or “AUX abc”—defines an auxiliary variable, named "abc", which will be read for each river reach as part of Items 4 and 6. Up to five variables can be specified, each of which must be preceded by "AUXILIARY" or "AUX." These variables will not be used by the Ground-Water Flow Process Package, but they will be available for use by other processes. The auxiliary variable values will be read after the Rbot variable.

Auxiliary parameters are used with MODPATH and the Ground-Water Transport process.  In MODPATH, the variable IFACE is specified with an auxiliary parameter.  In the Ground-Water Transport process, concentration is specified with an auxiliary parameter.  The option COMPACT BUDGET AUXILIARY must be used in the Output Control when IFACE is used for MODPATH.

CBCALLOCATE” or “CBC”—indicates that memory should be allocated to store cell-by- cell flow for each river reach in order to make these flows available for use in other packages. (This option is not available in MODFLOW-2005.)
NOPRINT—suppresses the printing of cell lists.

Repeat Items 3 and 4 for each NPRIV parameter.

3.

[PARNAM PARTYP Parval NLST [INSTANCES NUMINST]]

If PARNAM is to be a time-varying parameter, the keyword “INSTANCES” and a value for NUMINST must be entered.

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 RIV Package, the only allowed parameter type is RIV, which defines values of riverbed hydraulic conductance.

Parval—is the parameter value. This parameter value may be overridden by a value in the Sensitivity Process input file or by a value generated by the Parameter-Estimation Process.

NLST—is the number of river reaches that are included in a non-time-varying parameter or in each instance of a time-varying parameter.

INSTANCES – is an optional keyword that designates a parameter as time varying. The keyword is case-insensitive; that is, it may be entered in any combination of upper- and lower-case letters. If INSTANCES is present, it must be followed by a value for NUMINST. If INSTANCES is absent, PARNAM is non-time-varying and NUMINST should not be present.

NUMINST – is the number of instances that are included in the definition of a time-varying parameter, where each instance is a list of river reaches and associated properties. If the keyword INSTANCES is present, NUMINST must be present and must be at least 1. If the keyword INSTANCES is absent, NUMINST should not be present.

4a.

[INSTNAM]

Item 4a is read only if PARNAM is time-varying. NUMINST repetitions of Item 4 (parts a and b) are read. After each repetition of Item 4a, NLST repetitions of Item 4b are read.

INSTNAM – is the name of an instance associated with the parameter PARNAM specified in the corresponding Item 3. The 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. Names entered for INSTNAM must be unique for any given parameter, but names may be reused for instances associated with different parameters.

4b.

Layer Row Column Stage  Condfact Rbot [xyz]

NLST repetitions of Item 4b are required; they are read by module ULSTRD. (SFAC of the ULSTRD utility module applies to Condfact). The NLST repetitions of Item 4b follow each repetition of Item 4a when PARNAM is time-varying.

Layer—is the layer number of the cell containing the river reach.

Row—is the row number of the cell containing the river reach.

Column—is the column number of the cell containing the river reach.

Stage—is the head in the river.

Condfact—is the factor used to calculate riverbed hydraulic conductance from the parameter value. The conductance is the product of Condfact and the parameter value.

Rbot—is the elevation of the bottom of the riverbed.

[xyz]—represents any auxiliary variables for a river reach that have been defined in Item 2. The auxiliary variables must be present in each repetition of Items 4 and 6 if they are defined in Item 2.

FOR EACH STRESS PERIOD

5.

ITMP NP

ITMP—is a flag and a counter.

If ITMP < 0, non-parameter river data from the last stress period will be reused.
If ITMP 0, ITMP will be the number of non-parameter reaches read for the current stress period.

# NP—is the number of parameters in use in the current stress period.

6.

Layer Row Column Stage  Cond Rbot [xyz]

ITMP repetitions of Item 6 are read by module ULSTRD if ITMP > 0. (SFAC of the ULSTRD utility module applies to Cond.) Item 6 is not read if ITMP is negative or 0.

Layer—is the layer number of the cell containing the river reach.

Row—is the row number of the cell containing the river reach.

Column—is the column number of the cell containing the river reach.

Stage—is the head in the river.

Cond—is the riverbed hydraulic conductance.

Rbot—is the elevation of the bottom of the riverbed.

[xyz]—represents any auxiliary variables for a river reach that have been defined in Item 2. The auxiliary variables must be present in each repetition of Items 4 and 6 if they are defined in Item 2.

7.

[Pname [Iname]]

Item 7 is repeated NP times. It is not read if NP is negative or 0. Iname is read if Pname is a time-varying parameter.

Pname—is the name of a parameter that is being used in the current stress period. NP parameter names will be read.

Iname – is an instance name that is read only if Pname is a time-varying parameter. Multiple instances of the same time-varying parameter are not allowed in a stress period.