Please enable JavaScript to view this site.

Online Guide to MODFLOW

Navigation: Parameter Estimation Process

PES - Parameter-Estimation Process Input File

Scroll Prev Top Next More

Abbreviation in Name file

PES

Purpose

The Parameter-Estimation Process is used to calibrate a model by modifying parameter values until the fit between the simulated values and observations is optimal.

Documentation

Hill, M.C., Banta, E.R., Harbaugh, A.W., and Anderman, E.R., 2000, MODFLOW-2000, the U.S. Geological Survey modular ground-water model -- User guide to the Observation, Sensitivity, and Parameter-Estimation Processes and three post-processing programs: U.S. Geological Survey Open-File Report 00-184, 210 p.

Supported in

MODFLOW-2000

Other Notes

The Parameter-Estimation Process requires that the flow, observation, and sensitivity processes all be used.

See Package Incompatibilities

The Parameter Estimation Process is not included in MODFLOW-2005 or MODFLOW-LGR.  External programs such as UCODE must be used instead.

Input Instructions

Input for the Parameter-Estimation Process is read from a file that is specified with file type "PES" in the name file.

Data Set 0

[#Text]

Item 0 is optional and can include as many lines as desired. Each line needs to begin with

the “#” character in the first column.

Data Set 1

MAX-ITER MAX-CHANGE TOL SOSC (free format)

Item 1 includes variables that are most often changed by the user.

Data Set 2

IBEFLG IYCFLG IOSTAR NOPT NFIT SOSR  RMAR RMARM IAP (free format)

Item 2 includes variables that control the modified Gauss-Newton calculations and are sometimes changed by the user.

Data Set 3

IPRCOV IPRINT LPRINT (free format)

Item 3 includes variables that control printing.

Data Set 4

CSA FCONV LASTX (free format)

Item 4 includes variables that control the modified Gauss-Newton calculations and are rarely changed by the user.

Data Set 5

NPNG IPR MPR (free format)

Item 5 includes variables that indicate whether additional items need to be read. The additional possible items are as follows.

Data Set 6

If NPNG is greater than zero, read item 6 once.

PARNEG(1), PARNEG(2), . . . , PARNEG(NPNG) (free format)

Data Set 7

If IPR is greater than zero, read item IPR repetitions of item 7. Parameters that appear in item 7 may not appear in item 10.

NIPRNAM BPRI PLOT-SYMBOL (free format)

If IPR is greater than zero, read items 8 and 9 once. The size of item 9 depends on IPR

Data Set 8

IWTP (free format)

Data Set 9

WTP(1,1), WTP(1,2), ..., WTP(1,IPR) (free format)

WTP(2,1), WTP(2,2), ..., WTP(2,IPR) (new line; free format)

...

WTP(IPR,1), WTP(IPR,2), ..., WTP(IPR,IPR) (new line; free format)

Data Set 10

If MPR is greater than zero, read MPR repetitions of item 10. Examples are shown after the next section. Parameters that appear in item 7 may not appear in item 10.

EQNAM PRM "=" [SIGN] [COEF "*"] PNAM [SIGN [COEF "*"] PNAM [SIGN…]] "STAT" STATP STAT-FLAG PLOT-SYMBOL (free format: maximum of 200 characters; one or more spaces must separate all words, numbers, and symbols)

Example Input File

# Example PES file

#

15 2.0 1E-5 1E-4 ITEM 1: MAX-ITER MAX-CHANGE TOL SOSC

0 0 0 0 0 0. .001 1.5 0 ITEM 2: IBEFLG IYCFLG IOSTAR NOPT NFIT SOSR RMAR RMARM IAP

2 0 0 ITEM 3: IPRCOV IPRINT LPRINT

0.08 0.0 0 ITEM 4: CSA FCONV LASTX

0 0 6 ITEM 5: NPNG IPR MPR

EQ-Q_3 -97000. = WQ_3 STAT 1940 1 9 ITEMS 10: 6 (MPR) PRIOR-INFO EQUATIONS

EQ-Q_4 -51000. = WQ_4 STAT 1020 1 9

EQ-HK_2 420. = HK_2 STAT 84. 1 9

EQ-RCH_1 0.0004 = RCH_ZONE_1 stat 1.2E-4 1 9

EQ-RCH_3 1.7E-4 = RCH_ZONE_3 stat 5.1E-5 1 9

EQ-KRB_1 0.08 = RIVERS stat 0.008 1 9

 

In this example, WQ_3, WQ_4, HK_2, RCH_ZONE_1, RCH_ZONE_3, and RIVERS are

parameter names, which need to be defined in input files for Groundwater Flow Process

packages and need to be listed with ISENS > 0 in the Sensitivity Process input file. The program

matches parameter names among the various input files in a case-insensitive manner. EQ-Q_3,

EQ-Q_4, EQ-HK_2, EQ-RCH_1, EQ-RCH_3, and EQ-KRB_1 are names assigned to the 6 (MPR)

prior-information equations.

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