14.07.2013 Views

User's guide of Proceessing Modflow 5.0

User's guide of Proceessing Modflow 5.0

User's guide of Proceessing Modflow 5.0

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

150 Processing <strong>Modflow</strong><br />

- NPHINORED is the first termination criterion. A value <strong>of</strong> 3 is <strong>of</strong>ten suitable. If PEST<br />

has failed to lower the objective function over NPHINORED successive iterations, the<br />

program stops.<br />

- RELPARSTP and NRELPAR represent the second termination criterion. If the<br />

< Options:<br />

magnitude <strong>of</strong> the maximum relative parameter change between optimization iterations<br />

is less than RELPARSTP over NRELPAR successive iterations, the program stops. The<br />

relative parameter change between optimization iterations for any parameter is calculated<br />

using equation 3.55. For many cases, a value <strong>of</strong> 0.01 for RELPARSTP and a value <strong>of</strong> 3<br />

for NRELPAR is adequate.<br />

When the optimization process is complete, one <strong>of</strong> the termination criteria having been met<br />

(or perhaps another termination criterion such as zero objective function or zero objective<br />

function gradient, for which no user-supplied settings are required), PEST writes some<br />

information concerning the optimized parameter set to its run record file PESTCTL.REC.<br />

This file is saved in the data directory <strong>of</strong> your model. It tabulates the optimal values and the<br />

95% confidence intervals pertaining to all adjustable parameters. It also tabulates the<br />

model-calculated values based on these parameters, together with the residuals, i.e., the<br />

differences between measured and model-calculated values. If you wish, PEST will write<br />

the parameter covariance matrix, the parameter correlation coefficient matrix and the<br />

matrix <strong>of</strong> normalised eigenvectors <strong>of</strong> the covariance matrix to the run record file<br />

PESTCTL.REC.<br />

If the option Save data for a possible restart is checked, PEST will dump the contents<br />

<strong>of</strong> many <strong>of</strong> its data arrays to a binary file at the beginning <strong>of</strong> each optimization iteration; this<br />

allows PEST to be restarted later if execution is prematurely terminated. If subsequent PEST<br />

execution is initiated using the "/r" command line switch (see the PEST manual for details),<br />

it will recommence execution at the beginning <strong>of</strong> the iteration during which it was<br />

interrupted.<br />

If the option Include decimal point even if redundant is not checked, PEST will omit<br />

the decimal point from parameter values on model input files if the decimal point is<br />

redundant, thus making room for the use <strong>of</strong> one extra significant figure. If this option is<br />

checked, PEST will ensure that the decimal point is always present.<br />

3.6.5 PEST (Inverse Modeling)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!