The behavior of Point-Stat is controlled by the contents of the configuration file passed to it on the command line. The default Point-Stat configuration file may be found in the $MET_BASE/config/PointStatConfig_default file. The configuration used by the test script may be found in the $MET_BASE/PointStatConfig file. Prior to modifying the configuration file, users are advised to make a copy of the default:
Open up the $MET_TUTORIAL_DATA/config/PointStatConfig_tutorial file for editing with your preferred text editor.
The configurable items for Point-Stat are used to specify how the verification is to be performed. The configurable items include specifications for the following:
- The verification grid.
- The forecast fields to be verified at the specified vertical levels.
- The threshold values to be applied.
- The economic cost-loss value ratios to be evaluated.
- The reference climatological dataset.
- 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, configurable lat/lon polylines, or individual stations.
- The confidence interval methods to be used.
- The interpolation methods to be used.
- The types of verification 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/README file. Please take some time to review them.
For this tutorial, we'll configure Point-Stat to verify the model temperature at two vertical levels and winds at the surface. However, Point-Stat may be configured to verify as many or as few model variables as you desire. The sample input forecast file is not on the NCEP Grid 212 domain. However, we'll use the NCEP Grid 212 domain to define a masking region for our data. Edit the $MET_TUTORIAL_DATA/config/PointStatConfig_tutorial file as follows:
- Set wind_thresh = [ >0.0, >=1.0, >=5.0, >=8.0 ];
To indicate that we'd like VL1L2 lines computed using these thresholds on the wind speeds. - In the fcst dictionary, set
field = [
{
name = "TMP";
level = "Z2";
cat_thresh = [ >278, >283, >288 ];
},
{
name = "TMP";
level = "P750-850";
cat_thresh = [ >278 ];
},
{
name = "UGRD";
level = "Z10";
cat_thresh = [ >=5.0 ];
},
{
name = "VGRD";
level = "Z10";
cat_thresh = [ >=5.0 ];
}
];To verify 2-meter temperature, 10-meter winds, and temperature fields between 750hPa and 850hPa and apply the categorical thresholds listed. TMP is in Kelvin and the U and V components of wind are in m/s.
- Retain the default settings of obs = fcst;
To use the settings from the fcst dictionary above. - Set message_type = [ "ADPUPA", "ADPSFC" ];
To verify using these 2 observation types. - In the mask dictionary, set grid = [ "G212" ];
To accumulate statistics over NCEP Grid 212 domain. - In the mask dictionary, set poly = [ "${MET_BASE}/poly/EAST.poly", "${MET_BASE}/poly/WEST.poly" ];
To accumulate statistics over the regions defined by the EAST and WEST polyline files. - In the interp dictionary, set
type = [
{
method = NEAREST;
width = 1;
},
{
method = UW_MEAN;
width = 5;
}
];To indicate that the forecast values should be interpolated to the observation locations using the nearest neighbor method and by averaging the forecast values over the 5 by 5 box surrounding the observation location.
- Set
output_flag = {
fho = NONE;
ctc = BOTH;
cts = BOTH;
mctc = BOTH;
mcts = BOTH;
cnt = BOTH;
sl1l2 = BOTH;
sal1l2 = NONE;
vl1l2 = BOTH;
val1l2 = NONE;
vcnt = NONE;
pct = NONE;
pstd = NONE;
pjc = NONE;
prc = NONE;
ecnt = NONE; // Only for HiRA.
eclv = NONE;
mpr = NONE;
}To indicate that the contingency table counts (CTC), contingency table statistics (CTS), multi-category contingency table counts (MCTC), multi-category contingency table statistics (MCTS), continuous statistics (CNT), scalar partial sum (SL1L2), and vector partial sum (VL1L2) line types should be output.
Save and close this file.
