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

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

156 Processing <strong>Modflow</strong><br />

P 10 = 0.5 x P1 & 0.5 x P2 stat 0.2 flag 2 plot 4<br />

P 1.03E-3 = 1.0 x P1 & 1.45 x P2 & -3.9 x P3 & 2.123 x P4 stat 2.0 flag 2 plot 3<br />

P 1.344 = 2.12 x P3 & 3.2 x P6 stat 2.2 flag 2 plot 1<br />

< Control Data<br />

The control data are used to control the sensitivity and regression calculations and define the<br />

inversion algorithm. The items <strong>of</strong> the control data (Fig. 3.56) are described below.<br />

- Convergence criterion (TOL): When parameter values change less than this fractional<br />

amount between regression iterations, parameter estimation converges (0.01 is<br />

recommended).<br />

- Convergence criterion (SOSR): When the sum-<strong>of</strong>-squared, weighted residuals change<br />

less than this fractional amount over three regression iterations, parameter estimation<br />

converges. Ideally, for the final results convergence is achieved by satisfying the TOL<br />

criterion so that SOSR can be equal to 0.0 (in which case SOSR is not used as a<br />

convergence criteria). Values <strong>of</strong> SOSR <strong>of</strong> 0.01 and even 0.1 can be useful, however, in<br />

the early stages <strong>of</strong> model calibration because it stops the regression when it is not<br />

progressing.<br />

- Maximum number <strong>of</strong> regression iterations (MAX-ITER) is self-explanatary. Starting<br />

with twice the number <strong>of</strong> parameters is recommended.<br />

- Maximum fractional parameter change (MAX-CHANGE) is the maximum fractional<br />

change <strong>of</strong> a parameter value allowed in one regression iteration. For example, if MAX-<br />

CHANGE = 2.0, a parameter value <strong>of</strong> 1.0 will not be allowed to change by more than 2.0<br />

(MAX-CHANGE times the parameter value). Consequently, the new value will be<br />

between -1.0 and 3.0. A parameter value <strong>of</strong> 2.0 will not be allowed to change more than<br />

a value <strong>of</strong> 4.0 (again, MAX-CHANGE times the parameter value), and the new value will<br />

be between -2.0 and 6.0). This maximum change is applied to the physical parameter<br />

value, not its log transformed value. Exceptions are discussed in Hill (1998, Appendix<br />

B). MAX-CHANGE = 2.0 is common, but smaller values may help an oscillating<br />

regression to converge.<br />

- Differencing method controls the method used to calculate sensitivities during the<br />

parameter-estimation iterations. Starting with the forward differencing method is<br />

recommended.<br />

- Apply quasi-Newton update when sum-<strong>of</strong>-squared, weighted residuals changes less<br />

than 0.01 over three regression iteration. According to Hill (1998), applying the quasi-<br />

Newton update may facilitate convergence <strong>of</strong> highly nonlinear problems.<br />

3.6.6 UCODE (Inverse Modeling)

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

Saved successfully!

Ooh no, something went wrong!