Download Getting Started with PEST

Transcript
model input files using the input template files. Then the model (batch file) is called. If
the model runs successfully, the model will generate a set of output files. The modelgenerated outputs in the model output files, which will be compared against the
corresponding observations, are then read by PEST through using model output
instruction file. At this stage, the objective function and Jacobian matrix are calculated,
based on which the PEST will make its decision for next iteration until one of its stopping
criteria is met. The stop criteria are specified in PEST control file as you may suspect.
PEST Control File
Input Template
File
Model Input
(input files)
PEST
Optimization
Algorithm
Model (model
batch file)
Output
Instruction File
Model Output
(output files)
Figure 1. Schematic diagram of PEST optimization process
1.1 Model batch file
Model batch file is simple. It can be as short as one line – the model DOS command. But
if the model output files are not text format files; or you plan to use TSPROC as your
model post-processor, the model batch file will get slightly longer, like in Example 1.
The purpose of each command in Example 1 is briefly explained as follows. The first line
of Example 1 (@echo off) is to tell the computer system not to display the commands
while executing them. The second line (del basins.rch > nul) is to delete the
SWAT output file basins.rch that has been generated by previous modeling run. If it
had not been deleted, even though the following SWAT run is not successful, PEST
would still have read model outputs from the existing basins.rch file, which
obviously should be avoided. The suffix of the second line (> nul) is to suppress the
deleting command to be displayed in screen so that the messages issued by system
commands (or model) will not interfered with the messages issued by PEST. The third
line is a SWAT model command line. The fourth line is a short Fortran program I coded
to convert partial information contained in the SWAT basins.rch file to a Site
Sample File (SSF), a file format that the TSPROC can read. The information of what
2