Skip to content

Commit

Permalink
Per #1908, ci-run-unit more changes to Ensemble-Stat. Rename ensemble…
Browse files Browse the repository at this point in the history
…_flag to nc_orank_flag and update its contents. Also make it configurable separetly for each vx task. Also update the User's Guide with this info.
  • Loading branch information
JohnHalleyGotway committed Sep 19, 2022
1 parent 660effe commit 7254cdc
Show file tree
Hide file tree
Showing 17 changed files with 207 additions and 133 deletions.
10 changes: 7 additions & 3 deletions data/config/EnsembleStatConfig_default
Original file line number Diff line number Diff line change
Expand Up @@ -215,6 +215,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = NONE;
Expand All @@ -235,13 +236,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
28 changes: 17 additions & 11 deletions docs/Users_Guide/ensemble-stat.rst
Original file line number Diff line number Diff line change
Expand Up @@ -345,27 +345,33 @@ The **output_flag** array controls the type of output that is generated. Each fl
_____________________

.. code-block:: none
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}
The **ensemble_flag** specifies which gridded verification output types should be written to the Observation Rank (**_orank.nc**) NetCDF file. Setting the flag to TRUE produces output of the specified field, while FALSE produces no output for that field type. The flags correspond to the following output line types:
The **nc_orank_flag** specifies which gridded verification output types should be written to the Observation Rank (**_orank.nc**) NetCDF file. This output file is only created when gridded observations have been provided with the -grid_obs command line option. Setting the flag to TRUE produces output of the specified field, while FALSE produces no output for that field type. The flags correspond to the following output line types:

1. Grid Latitude and Longitude Fields

2. The grid area weights.
2. Ensemble mean field

3. Raw observation values

4. Observation ranks

3. Observation values, ranks, and probability-integral transform values.
5. Observation probability-integral transform values

4. Ensemble Valid Data Count
6. Ensemble valid data count

5. Ensemble Mean Field
7. Grid area weight values

__________________

Expand Down Expand Up @@ -410,13 +416,13 @@ ensemble_stat_PREFIX_YYYYMMDD_HHMMSSV.stat where PREFIX indicates the user-defin
The output ASCII files are named similarly:


ensemble_stat_PREFIX_YYYYMMDD_HHMMSSV_TYPE.txt where TYPE is one of ecnt, rps, rhist, phist, relp, orank, and ssvar to indicate the line type it contains.
ensemble_stat_PREFIX_YYYYMMDD_HHMMSSV_TYPE.txt where TYPE is one of elements of the **output_flag** configuration option to indicate the line type it contains.


When verification against gridded analyses is performed, Ensemble-Stat can produce output NetCDF files using the following naming convention:


ensemble_stat_PREFIX_YYYYMMDD_HHMMSSV_orank.nc contains gridded fields of observation ranks when the -grid_obs command line option is used. Its contents are specified by the **ensemble_flag** configuration option.
ensemble_stat_PREFIX_YYYYMMDD_HHMMSSV_orank.nc contains gridded fields of observation ranks when the -grid_obs command line option is used. Its contents are specified by the **nc_orank_flag** configuration option.

The Ensemble-Stat tool can compute the following statistics for the fields specified in the fcst and obs dictionaries of the configuration file:

Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig
Original file line number Diff line number Diff line change
Expand Up @@ -181,6 +181,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = BOTH;
Expand All @@ -201,13 +202,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig_MASK_SID
Original file line number Diff line number Diff line change
Expand Up @@ -173,6 +173,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = STAT;
Expand All @@ -193,13 +194,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = FALSE;
rank = FALSE;
pit = FALSE;
vld_count = FALSE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig_climo
Original file line number Diff line number Diff line change
Expand Up @@ -202,6 +202,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = BOTH;
Expand All @@ -222,13 +223,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = FALSE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = FALSE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig_grid_weight
Original file line number Diff line number Diff line change
Expand Up @@ -190,6 +190,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = STAT;
Expand All @@ -210,13 +211,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = TRUE;
mean = TRUE;
raw = FALSE;
rank = FALSE;
pit = FALSE;
vld_count = FALSE;
mean = TRUE;
weight = TRUE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
6 changes: 4 additions & 2 deletions internal/test_unit/config/EnsembleStatConfig_one_cdf_bin
Original file line number Diff line number Diff line change
Expand Up @@ -195,6 +195,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = BOTH;
Expand All @@ -214,9 +215,10 @@ output_flag = {
////////////////////////////////////////////////////////////////////////////////

//
// Ensemble product output types
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = false;
nc_orank_flag = false;

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

Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig_python
Original file line number Diff line number Diff line change
Expand Up @@ -177,6 +177,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = STAT;
Expand All @@ -197,13 +198,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = FALSE;
rank = FALSE;
pit = FALSE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig_qty_inc_exc
Original file line number Diff line number Diff line change
Expand Up @@ -190,6 +190,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = BOTH;
Expand All @@ -210,13 +211,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig_single_file_grib
Original file line number Diff line number Diff line change
Expand Up @@ -226,6 +226,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = STAT;
Expand All @@ -246,13 +247,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
10 changes: 7 additions & 3 deletions internal/test_unit/config/EnsembleStatConfig_single_file_nc
Original file line number Diff line number Diff line change
Expand Up @@ -232,6 +232,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = STAT;
Expand All @@ -252,13 +253,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
10 changes: 7 additions & 3 deletions scripts/config/EnsembleStatConfig
Original file line number Diff line number Diff line change
Expand Up @@ -227,6 +227,7 @@ interp = {

//
// Statistical output types
// May be set separately in each "obs.field" entry
//
output_flag = {
ecnt = BOTH;
Expand All @@ -247,13 +248,16 @@ output_flag = {

//
// Gridded verification output types
// May be set separately in each "obs.field" entry
//
ensemble_flag = {
nc_orank_flag = {
latlon = TRUE;
weight = FALSE;
mean = TRUE;
raw = TRUE;
rank = TRUE;
pit = TRUE;
vld_count = TRUE;
mean = TRUE;
weight = FALSE;
}

////////////////////////////////////////////////////////////////////////////////
Expand Down
4 changes: 3 additions & 1 deletion src/basic/vx_config/config_constants.h
Original file line number Diff line number Diff line change
Expand Up @@ -739,8 +739,9 @@ static const char conf_key_nmep_smooth[] = "nmep_smooth";
static const char conf_key_skip_const[] = "skip_const";
static const char conf_key_rng_type[] = "rng.type";
static const char conf_key_rng_seed[] = "rng.seed";
static const char conf_key_ensemble_flag[] = "ensemble_flag";
static const char conf_key_nc_orank_flag[] = "nc_orank_flag";
static const char conf_key_rank_flag[] = "rank";
static const char conf_key_pit_flag[] = "pit";
static const char conf_key_ssvar_bin[] = "ens_ssvar_bin_size";
static const char conf_key_phist_bin[] = "ens_phist_bin_size";
static const char conf_key_prob_cat_thresh[] = "prob_cat_thresh";
Expand All @@ -760,6 +761,7 @@ static const char conf_key_ens_field[] = "ens.field";
static const char conf_key_ens_ens_thresh[] = "ens.ens_thresh";
static const char conf_key_ens_vld_thresh[] = "ens.vld_thresh";
static const char conf_key_normalize[] = "normalize";
static const char conf_key_ensemble_flag[] = "ensemble_flag";
static const char conf_key_mean_flag[] = "mean";
static const char conf_key_stdev_flag[] = "stdev";
static const char conf_key_minus_flag[] = "minus";
Expand Down
Loading

0 comments on commit 7254cdc

Please sign in to comment.