MET Online Tutorial for METv3.0

Ensemble-Stat Tool: Configure

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 METv3.0/data/config/EnsembleStatConfig_default file. The configuration used by the test script may be found in the METv3.0/scripts/config/EnsembleStatConfig file. Prior to modifying the configuration file, users are advised to make a copy of the default:

    cp data/config/EnsembleStatConfig_default tutorial/config/EnsembleStatConfig_APCP_24

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 to compute rank histograms. The configurable items include specifications for the following:

  • Section 1: Ensemble Processing
    • 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: Rank Histograms
    • The forecast 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 comments of the configuration file itself. 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 tutorial/config/EnsembleStatConfig_APCP_24 file and edit it as follows:

  • Set ens_field[] = [ "APCP/A24" ];
    To read 24-hour accumulated precipitation from the input GRIB files.
  • Set ens_thresh[] = [ "gt0.0 ge5.0 ge10.0" ];
    To compute ensemble relative frequencies for these thresholds.
  • Set fcst_field[] = [ "APCP/A24" ];
    To also verify the 24-hour accumulated precipitation fields.
  • Set message_type[] = [ "ADPSFC" ];
    To verify against surface observations.
  • Set mask_grid[] = [ "FULL" ];
    To verify over the full domain.
  • Set mask_poly[] = [ "MET_BASE/data/poly/NWC.poly",
    "MET_BASE/data/poly/SWC.poly" ];

    To also verify over the northwest coast (NWC) and southwest coast (SWC) subregions.

Save and close this file.