Download we include it here.

Transcript
Advamced Regularisation and SVD-Assist
8-22
original model batch file which are transferred directly to the new model batch file
svdabatch.bat.)
As was mentioned above, in the original base PEST control file the command to run the
model must be the name of a batch file – not the name of an executable program;
SVDAPREP needs to modify this batch file to produce svdabatch.bat. Consequently, if the
suffix “.bat” does not comprise the last four letters of the model command, SVDAPREP will
cease execution with an appropriate error message.
8.5.4.3 PARCALC
The user never needs to run the PARCALC utility program. This is part of the model run by
PEST when implementing SVD-assisted parameter estimation. Nevertheless its functionality
will be described.
PARCALC reads the current values of super parameters being estimated by PEST. It also
reads base parameter eigenvalues as calculated by PEST on the basis of the JCO file cited in
the “svd assist” section of the super PEST control file; recall that super parameter definition
is based on these eigenvectors. Using this information PARCALC calculates current base
parameter values. These are then provided to the model through the template files what were
originally used by the base PEST control file.
In calculating current base parameter values, PARCALC does the following:1. It respects upper and lower base parameter bounds as provided in its input file
parcalc.in, these having been extracted from the base PEST control file.
2. The values of fixed base parameters are unaltered. (These parameters are not used in
the calculation of super parameters.)
3. Parameter value ratios for tied parameters are maintained. (Tied parameters are also
not used in the calculation of super parameters.)
PARCALC is run by typing the command “parcalc” at the screen prompt. It reads an input
file named parcalc.in. This file is actually written by PEST on the basis of the template file
parcalc.tpl, also written by PEST. Use of the template file mechanism for the provision of
PARCALC input data allows parallel PEST to transfer PARCALC input data to multiple
computers involved in a parallelisation process.
8.5.4.4 PICALC
PICALC requires no user input. It reads an input file named picalc.in and writes an output
file named picalc.out. The first of these files is written by PARCALC on the basis of the
SVDAPREP-prepared template file picalc.tpl. This contains current values of base
parameters, as well as parameter coefficients used in prior information equations. PICALC
evaluates these equations on the basis of current base parameter values and records the
outcomes of these calculations to a file named picalc.out.