The behavior of Ensemble-Stat is controlled by the contents of the configuration file passed to it on the command line. The default Ensemble-Stat configuration file may be found in the $MET_BASE/config/EnsembleStatConfig_default file. The configuration used by the test script may be found in the met-8.0/scripts/config/EnsembleStatConfig file. Prior to modifying the configuration file, users are advised to make a copy of the default:
The configurable items for Ensemble-Stat are broken out into two sections. The first section specifies how the ensemble should be processed to derive summary fields, such as the ensemble mean and spread. The second section specifies how the ensemble should be verified directly, such as the computation of rank histograms and spread/skill. The configurable items include specifications for the following:
- Section 1: Ensemble Processing (ens dictionary)
- The ensemble fields to be summarized at the specified vertical level or accumulation interval.
- The threshold values to be applied in computing ensemble relative frequencies (e.g. the percent of ensemble members exceeding some threshold at each point).
- Thresholds to specify how many of the ensemble members must actually be present with valid data.
- Section 2: Verification (fcst and obs dictionaries)
- The forecast and observation fields to be verified at the specified vertical level or accumulation interval.
- The matching time window for point observations.
- The type of point observations to be matched to the forecasts.
- The areas over which to aggregate statistics - as predefined grids or configurable lat/lon polylines.
- The interpolation or smoothing methods to be used.
You may find a complete description of the configurable items in the MET Users Guide or in the $MET_BASE/config/READMEfile. Please take some time to review them.
For this tutorial, we'll configure Ensemble-Stat to summarize and verify 24-hour accumulated precipitation. While we'll run Ensemble-Stat on a single field, please note that it may be configured to operate on multiple fields. The ensemble we're verifying consists of 6 members defined over the west coast of the United States.
Open up the $MET_TUTORIAL_DATA/config/EnsembleStatConfig_APCP_24 file and edit it as follows:
In the ens dictionary, set
{
name = "APCP";
level = [ "A24" ];
cat_thresh = [ >0, >=5.0, >=10.0 ];
}
];
To read 24-hour accumulated precipitation from the input GRIB files and compute ensemble relative frequencies for the thresolds listed.
In the fcst dictionary, set
{
name = "APCP";
level = [ "A24" ];
}
];
To also verify the 24-hour accumulated precipitation fields.
- In the fcst dictionary, set message_type = [ "ADPSFC" ];
To verify against surface observations. - In the mask dictionary, set grid = [ "FULL" ]; To accumulate statistics over the full model domain.
In the mask dictionary, set
"${MET_BASE}/poly/SWC.poly" ];
To also verify over the northwest coast (NWC) and southwest coast (SWC) subregions.
In the output_flag dictionary, set
ecnt = BOTH;
rhist = BOTH;
phist = BOTH;
orank = BOTH;
ssvar = BOTH;
relp = BOTH;
}
To write to the ".stat" output file as well the optional "_type.txt" file, a more readable ASCII file sorted by line type.
Save and close this file.
