Please enable JavaScript to view this site.

Online Guide to MODFLOW-NWT

Navigation: Groundwater Flow Process > Subsidence

SUB - Subsidence and Aquifer-System Compaction Package

Scroll Prev Top Next More

Abbreviation in Name file

SUB

Purpose

The  Subsidence and Aquifer-System Compaction package is used for simulating the drainage, changes in groundwater storage, and compaction of aquifers, interbeds and confining units that constitute an aquifer system. Delays in the release of groundwater from interbed storage, and thus delays in aquifer-system compaction, can be simulated. Delayed drainage and compaction in confining units can also be simulated.

Documentation

Höffmann, J., Leake, S.A., Galloway, D.L., and Wilson, A.M., 2003, MODFLOW-2000 Ground-Water Model--User Guide to the Subsidence and Aquifer-System Compaction (SUB) Package: U.S. Geological Survey Open-File Report 03-233, 44 p.

release.txt

Schmid, Wolfgang, and Hanson R.T., 2009, The Farm Process Version 2 (FMP2) for MODFLOW-2005 - Modifications and Upgrades to FMP1: U.S. Geological Survey Techniques in Water Resources Investigations, Book 6, Chapter A32, 102p. (http://pubs.er.usgs.gov/publication/tm6A32)

Hanson, R.T., Boyce, S.E., Schmid, Wolfgang, Hughes, J.D., Mehl, S.M., Leake, S.A., Maddock, Thomas, III, and Niswonger, R.G., 2014, One-Water Hydrologic Flow Model (MODFLOW-OWHM): U.S. Geological Survey Techniques and Methods 6–A51, 120 p., http://dx.doi.org/10.3133/tm6A51.

Related Packages

Interbed Storage package

SWT Subsidence and Aquifer-System Compaction Package

Supported in

MODFLOW-2000

MODFLOW-2005

MODFLOW-LGR

MODFLOW-CFP

MODFLOW-NWT

MODFLOW-OWHM

Other Notes

The Interbed Storage Package can not be used in the same model with the Subsidence and Aquifer-System Compaction Package.

Subsidence cannot be used in simulations in which stress periods other than the first are steady-state.

Comment lines can be inserted into the SUB input file anywhere except where arrays are to be read by one of the array-reader subroutines (U2DREL or U2DINT).   Each comment line must have the "#" character as its first character and may be up to 200 characters in length.

The UZF package does not take into account changes in the land surface due to subsidence.

 

A property from Data Set 7 through 12 or 14 through 18 can either be defined there or by a parameter, but not both. So if parameter PTYP is used, then its corresponding data set item is not read.

See the Frequently Asked Questions for information on how to read data from binary files generated by MODFLOW.

Input Instructions

Input for the SUB Package is read from the file that has the type “SUB” in the name file. Optional variables are shown in brackets. All single-valued variables in data items 1, 15, and 16, layer assignments for systems of interbeds in data items 2 and 3, and material properties in data item 9 are read in free format. Data items 1, 2, 3, and 15 consist of at most one record. Two-dimensional arrays in data items 4-8 and 10-14 are read with MODFLOW-2000 utility array readers U2DREL and U2DINT. For instructions on use of array readers, refer to Harbaugh and others (2000).

FOR EACH SIMULATION

Data Set 0

[#Text]

Item 0 is optional—“#” must be in column 1. Item 0 can be repeated multiple times.

 


Data Set 1

ISUBCB ISUBOC NNDB NDB NMZ NN AC1 AC2 ITMIN IDSAVE IDREST

(Enter integers for variables other than AC1 and AC2, which are floating-point variables.)

Data Set 2

[LN(NNDB)] if NNDB > 0

(Enter NNDB integers separated by one or more spaces or by commas.)

Data Set 3

[LDN(NDB)] if NDB > 0

(Enter NDB integers separated by one or more spaces or by commas.)

 


 


Data Set 4

[RNB(NCOL,NROW)] U2DREL if NDB > 0

(One array for each of the NDB systems of interbeds)

The following four arrays are needed to describe the initial conditions and properties of each of the NNDB systems of no-delay interbeds. All of the arrays (items 5–8) for system 1 are read first; then all of the arrays for the remaining systems.

Data Set 5

[HC(NCOL,NROW)] U2DREL if NNDB > 0

Data Set 6

[Sfe(NCOL,NROW)] U2DREL if NNDB > 0

Data Set 7

[Sfv(NCOL,NROW)] U2DREL if NNDB > 0

Data Set 8


[Com(NCOL,NROW)] U2DREL if NNDB > 0

 


 

 


Data Set 9

[DP(NMZ,3)] if NDB > 0

(Use one record for each material zone. Data item includes NMZ records, each with a value of vertical hydraulic conductivity, elastic specific storage, and inelastic specific storage.)

The following five arrays are needed to describe the initial conditions and properties of each of the NDB systems of delay interbeds. All of the arrays (items 10-14) for system 1 are read first; then all of the arrays for the remaining systems.

Data Set 10

[Dstart(NCOL,NROW)] U2DREL if NDB > 0 and IDREST 0

Data Set 11

[DHC(NCOL,NROW)] U2DREL if NDB > 0 and IDREST 0

Data Set 12


[DCOM(NCOL,NROW)] U2DREL if NDB > 0

 


 

 


Data Set 13

[DZ(NCOL,NROW)] U2DREL if NDB > 0

Data Set 14

[NZ(NCOL,NROW)] U2DINT if NDB > 0

Data Set 15

[Ifm1 Iun1 Ifm2 Iun2 Ifm3 Iun3 Ifm4 Iun4 Ifm5 Iun5 Ifm6 Iun6] if ISUBOC > 0

(Data item 15 consists of one record with 12 integers separated by one or more spaces or by commas)

Data Set 16

[ISP1 ISP2 ITS1 ITS2 Ifl1 Ifl2 Ifl3 Ifl4 Ifl5 Ifl6 Ifl7 Ifl8 Ifl9 Ifl10 Ifl11 Ifl12 Ifl13] if ISUBOC > 0.

Data item 16 consists of ISUBOC records with 17 integers separated by one or more spaces or by commas. Please see the section entitled “Package Output” for a detailed explanation of the use of data item 16.

The variables ISP1, ISP2, ITS1, ITS2, and Ifl1 through Ifl13 are used to control printing and saving of information generated by the SUB Package during program execution. The use of some of these variables is explained in more detail in the section entitled Package Output. The default condition for flags Ifl1 through Ifl13 is to not print or save the indicated information, except for printing budgets for no-delay interbeds for the last time step of each stress period.

Created with Help & Manual 8 and styled with Premium Pack Version 4 © by EC Software