Skip to content

Commit

Permalink
Per #1575, add tests for Grid-Stat and Point-Stat to exercise the mpr…
Browse files Browse the repository at this point in the history
…_column and mpr_thresh config file options.
  • Loading branch information
JohnHalleyGotway committed Apr 1, 2021
1 parent 8acc282 commit af73fb0
Show file tree
Hide file tree
Showing 4 changed files with 549 additions and 0 deletions.
274 changes: 274 additions & 0 deletions test/config/GridStatConfig_mpr_thresh
Original file line number Diff line number Diff line change
@@ -0,0 +1,274 @@
////////////////////////////////////////////////////////////////////////////////
//
// Grid-Stat configuration file.
//
// For additional information, please see the MET User's Guide.
//
////////////////////////////////////////////////////////////////////////////////

//
// Output model name to be written
//
model = "GFS";

//
// Output description to be written
// May be set separately in each "obs.field" entry
//
desc = "NA";

//
// Output observation type to be written
//
obtype = "GFSANL";

////////////////////////////////////////////////////////////////////////////////

//
// Verification grid
//
regrid = {
to_grid = NONE;
method = NEAREST;
width = 1;
vld_thresh = 0.5;
}

////////////////////////////////////////////////////////////////////////////////

//
// May be set separately in each "field" entry
//
censor_thresh = [];
censor_val = [];
mpr_column = [];
mpr_thresh = [];
cat_thresh = [];
cnt_thresh = [ NA ];
cnt_logic = UNION;
wind_thresh = [ NA ];
wind_logic = UNION;
eclv_points = 0.05;
nc_pairs_var_name = "";
nc_pairs_var_suffix = "";
rank_corr_flag = FALSE;

//
// Forecast and observation fields to be verified
//
fcst = {

name = "TMP";
level = "Z2";

field = [
{
desc = "NO_MPR_THRESH";
nc_pairs_var_suffix = desc;
},
{
mpr_column = [ "OBS-FCST" ];
mpr_thresh = [ >=-5&&<=5 ];
desc = "OBS_FCST_DIFF";
nc_pairs_var_suffix = desc;
},
{
mpr_column = [ "ABS(OBS-FCST)" ];
mpr_thresh = [ <=5 ];
desc = "ABS_OBS_FCST_DIFF";
nc_pairs_var_suffix = desc;
},
{
mpr_column = [ "ABS(OBS-CLIMO_MEAN)" ];
mpr_thresh = [ <=5 ];
desc = "ABS_OBS_CLIMO_MEAN_DIFF";
nc_pairs_var_suffix = desc;
},
{
mpr_column = [ "CLIMO_CDF" ];
mpr_thresh = [ >=0.25&&<=0.75 ];
desc = "CLIMO_CDF_IQR";
nc_pairs_var_suffix = desc;
}
];
}
obs = fcst;

////////////////////////////////////////////////////////////////////////////////

//
// Climatology mean data
//
climo_mean = fcst;
climo_mean = {

file_name = [ ${CLIMO_MEAN_FILE_LIST} ];

regrid = {
method = BILIN;
width = 2;
vld_thresh = 0.5;
}

time_interp_method = DW_MEAN;
day_interval = ${DAY_INTERVAL};
hour_interval = ${HOUR_INTERVAL};
}

climo_stdev = climo_mean;
climo_stdev = {
file_name = [ ${CLIMO_STDEV_FILE_LIST} ];
}

//
// May be set separately in each "obs.field" entry
//
climo_cdf = {
cdf_bins = 1;
center_bins = FALSE;
write_bins = TRUE;
}

////////////////////////////////////////////////////////////////////////////////

//
// Verification masking regions
//
mask = {
grid = [ "FULL" ];
poly = [];
}

////////////////////////////////////////////////////////////////////////////////

//
// Confidence interval settings
//
ci_alpha = [ 0.05 ];

boot = {
interval = PCTILE;
rep_prop = 1.0;
n_rep = 0;
rng = "mt19937";
seed = "";
}

////////////////////////////////////////////////////////////////////////////////

//
// Data smoothing methods
//
interp = {
field = BOTH;
vld_thresh = 1.0;
shape = SQUARE;

type = [
{
method = NEAREST;
width = 1;
}
];
}

////////////////////////////////////////////////////////////////////////////////

//
// Neighborhood methods
//
nbrhd = {
width = [ 1 ];
cov_thresh = [ >=0.5 ];
vld_thresh = 1.0;
shape = SQUARE;
}

////////////////////////////////////////////////////////////////////////////////

//
// Fourier decomposition
//
fourier = {
wave_1d_beg = [];
wave_1d_end = [];
}

////////////////////////////////////////////////////////////////////////////////

//
// Gradient statistics
// May be set separately in each "obs.field" entry
//
gradient = {
dx = [ 1 ];
dy = [ 1 ];
}

////////////////////////////////////////////////////////////////////////////////

//
// Distance Map statistics
// May be set separately in each "obs.field" entry
//
distance_map = {
baddeley_p = 2;
baddeley_max_dist = NA;
fom_alpha = 0.1;
zhu_weight = 0.5;
}

////////////////////////////////////////////////////////////////////////////////

//
// Statistical output types
//
output_flag = {
fho = NONE;
ctc = NONE;
cts = NONE;
mctc = NONE;
mcts = NONE;
cnt = NONE;
sl1l2 = STAT;
sal1l2 = NONE;
vl1l2 = NONE;
val1l2 = NONE;
vcnt = NONE;
pct = NONE;
pstd = NONE;
pjc = NONE;
prc = NONE;
eclv = NONE;
nbrctc = NONE;
nbrcts = NONE;
nbrcnt = NONE;
grad = NONE;
dmap = NONE;
}

//
// NetCDF matched pairs output file
//
nc_pairs_flag = {
latlon = FALSE;
raw = FALSE;
diff = TRUE;
climo = FALSE;
climo_cdp = FALSE;
weight = FALSE;
nbrhd = FALSE;
fourier = FALSE;
gradient = FALSE;
distance_map = FALSE;
apply_mask = FALSE;
}

////////////////////////////////////////////////////////////////////////////////

grid_weight_flag = COS_LAT;
tmp_dir = "/tmp";
output_prefix = "${OUTPUT_PREFIX}";
version = "V10.0.0";

////////////////////////////////////////////////////////////////////////////////
Loading

0 comments on commit af73fb0

Please sign in to comment.