From 7ad04ae410804b1d5306251795f95b7f03a936f5 Mon Sep 17 00:00:00 2001 From: johnhg Date: Thu, 9 Jun 2022 22:08:41 -0600 Subject: [PATCH] Bugfix #2170 main_v10.1 agg_ecnt (#2176) --- met/docs/Users_Guide/config_options.rst | 7 +++++++ met/src/libcode/vx_statistics/ens_stats.cc | 12 ++++++++++-- 2 files changed, 17 insertions(+), 2 deletions(-) diff --git a/met/docs/Users_Guide/config_options.rst b/met/docs/Users_Guide/config_options.rst index ae5ca26a40..ceb19d6b4c 100644 --- a/met/docs/Users_Guide/config_options.rst +++ b/met/docs/Users_Guide/config_options.rst @@ -58,6 +58,13 @@ The configuration file language supports the following data types: and 10 and "==1||==2" defines numbers exactly equal to 1 or 2. * Percentile Thresholds: + + * A threshold type (<, <=, ==, !=, >=, or >), followed by a percentile + type description (SFP, SOP, SCP, USP, CDP, or FBIAS), followed by a + numeric value, typically between 0 and 100. + + * Note that the two letter threshold type abbreviations (lt, le, eq, ne, + ge, gt) are not supported for percentile thresholds. * Thresholds may be defined as percentiles of the data being processed in several places: diff --git a/met/src/libcode/vx_statistics/ens_stats.cc b/met/src/libcode/vx_statistics/ens_stats.cc index d0bcfde0bd..d7c0e2c088 100644 --- a/met/src/libcode/vx_statistics/ens_stats.cc +++ b/met/src/libcode/vx_statistics/ens_stats.cc @@ -257,9 +257,17 @@ void ECNTInfo::set(const PairDataEnsemble &pd) { } } - // Compute ensemble mean based statistics, if possible + // Check if the ensemble mean based statistics were + // already computed by Stat-Analysis + if(!is_bad_data(pd.me)) { + me = pd.me; + rmse = pd.rmse; + me_oerr = pd.me_oerr; + rmse_oerr = pd.rmse_oerr; + } + // If not, compute them from the pairs, if possible // HiRA stores the ensemble mean as bad data - if(!pd.mn_na.is_const(bad_data_double)) { + else if(!pd.mn_na.is_const(bad_data_double)) { // Compute ME and RMSE values fbar = obar = ffbar = oobar = fobar = 0.0;