Bugfix #2170 main_v10.1 agg_ecnt

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:

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;