diff --git a/docs/Users_Guide/config_options.rst b/docs/Users_Guide/config_options.rst index ae5ca26a40..ceb19d6b4c 100644 --- a/docs/Users_Guide/config_options.rst +++ b/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/src/libcode/vx_statistics/ens_stats.cc b/src/libcode/vx_statistics/ens_stats.cc index d0bcfde0bd..d7c0e2c088 100644 --- a/src/libcode/vx_statistics/ens_stats.cc +++ b/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;