RES - Reservoir Package

  • Scroll to top of this topic Scroll to Top of Page

  • Print the current topic Print Topic

  • Expand/collapse expanding elements on page Show/Hide Expanders

Abbreviation in Name file

RES

Purpose

The Reservoir Boundary package simulates leakage between a reservoir and an underlying groundwater system as the reservoir area expands and contracts in response to changes in reservoir stage.  The user can specify beginning and ending stages for the reservoir in each stress period.

Unlike the Lake package, the head in reservoirs is unaffected by the interaction between the reservoir and the groundwater.

Documentation

Fenske, J.P., Leake, S.A., and Prudic, D.E., 1996, Documentation of a computer program (RES1) to simulate leakage from reservoirs using the modular finite-difference ground-water flow model (MODFLOW): U.S. Geological Survey Open-File Report 96-364, 51 p.

Related Packages

•Head-Dependent Flux Boundary packages

•River package

•Lake package

Supported in

•MODFLOW-2000

•MODFLOW-2005

•MODFLOW-LGR

•MODFLOW-CFP

•MODFLOW-NWT

•MODFLOW-OWHM

Input Instructions

Input for RES1 is read from the file that has file type "RES" in the name file.

FOR EACH SIMULATION

Data Set 1

NRES IRESCB NRESOP IRESPT NPTS

Format: I10 I10 I10 I10 I10

NRES is the number of reservoirs to be simulated. The minimum number of reservoirs that can be simulated is one and the maximum number is the number of active cells in the model grid.

IRESCB is a flag and unit number.

•If IRESCB > 0, it is the unit number on which cell-by-cell flow terms will be recorded whenever ICBCFL is set or "Save Budget" is specified in Output Control (see McDonald and Harbaugh, 1988, chap. 4, p. 14–15).

•If IRESCB ≤ 0, cell-by-cell flow terms will not be recorded.

NRESOP is a flag to determine the layer to which reservoirs are connected.

•If NRESOP = 1 all reservoirs are connected to layer 1.

•If NRESOP = 2 all reservoirs are connected to layer specified in array IRESL.

•If NRESOP = 3 all reservoirs are connected to the highest active cell in the vertical column of cells at each reservoir location.

IRESPT is a flag.

•If IRESPT > 0 reservoir stage, area, and volume will be printed each time step.

•If IRESPT ≤ 0 reservoir stage, area, and volume will not be printed each time step.

NPTS is the number of values in printed table of stage, volume, and area for each reservoir. First and last stage value are minimum and maximum elevations within area of potential inundation. A value of 15 or greater is recommended for detailed representation of stage-volume and stage-area relations.

Data Set 2

IRES(NCOL,NROW)

Module: U2DINT

IRES is an array that specifies the model cells that are included in each reservoir. Array values of “0" or less indicate the absence of reservoirs at corresponding cell locations. Values greater than “0” indicate that the cell is part of the corresponding reservoir. For example, values of “1” indicate the location of the first of NRES reservoirs and values of “2" indicate the location of the second of NRES reservoirs.

Data Set 3

IRESL(NCOL,NROW)

Module: U2DINT

IF NRESOP DOES NOT EQUAL 2, SKIP DATA ITEM 3

IRESL is an array that specifies the model layer to which each reservoir cell is connected.

Data Set 4

BRES(NCOL,NROW)

Module: U2DREL

BRES is an array that specifies land-surface elevation for all reservoirs. Values in the array are assigned to the reservoir denoted by the corresponding element in the IRES array. Values are not used for cells where the corresponding element in the IRES array is less than or equal to zero.

Data Set 5

HCres(NCOL,NROW)

Module: U2DREL

HCres is an array that specifies vertical hydraulic conductivity of the reservoir bed for all reservoirs. Values in the array are assigned to the reservoir denoted by the corresponding element in the IRES array. Values are not used for cells where the corresponding element in the IRES array is less than or equal to zero.

Data Set 6

Rbthck(NCOL,NROW)

Module: U2DREL

Rbthck is an array that specifies thickness of reservoir-bed sediments for all reservoirs. Values in the array are assigned to the reservoir denoted by the corresponding element in the IRES array. Values are not used for cells where the corresponding element in the IRES array is less than or equal to zero.

FOR EACH STRESS PERIOD

FOR EACH OF NRES RESERVOIRS

Data Set 7

Ststage Endstage

Format: F10.0 F10.0

Ststage is the value of reservoir stage at the start of the stress period. Because the package assigns stage values on the basis of linear interpolation to the end of each time step, the assigned reservoir stage will never equal Ststage unless Ststage and Endstage are equal. A total of NRES records with values of Ststage and Endstage are read at the start of each stress period. Stage is measured relative to the same datum used for the groundwater heads.

Endstage is the value of reservoir stage at the end of the stress period. Endstage is the value that will be assigned as stage for the corresponding reservoir for the last time step in the stress period. Stage is measured relative to the same datum used for the groundwater heads.