Per issue #2206, in set() function, added checks for crps_emp and crp…

…s_emp_fair, for values of 0.0. SL
dtcenter · Sep 12, 2022 · f9e17db · f9e17db
1 parent 8db8ce8
commit f9e17db
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/src/libcode/vx_statistics/ens_stats.cc b/src/libcode/vx_statistics/ens_stats.cc
@@ -231,8 +231,11 @@ void ECNTInfo::set(const PairDataEnsemble &pd) {
 
    // Compute empirical CRPS scores
    crps_emp   = pd.crps_emp_na.wmean(pd.wgt_na);
-   // Stub in for now?
+   if(is_eq(crps_emp, 0.0)) crps_emp = 0.0;
+
    crps_emp_fair = pd.crps_emp_fair_na.wmean(pd.wgt_na);
+   if(is_eq(crps_emp_fair, 0.0)) crps_emp_fair = 0.0;
+
    crpscl_emp = pd.crpscl_emp_na.wmean(pd.wgt_na);
    crpss_emp  = (is_bad_data(crps_emp)   ||
                  is_bad_data(crpscl_emp) ||