diff --git a/src/libcode/vx_statistics/ens_stats.cc b/src/libcode/vx_statistics/ens_stats.cc index d7c0e2c088..e9811b5fd3 100644 --- a/src/libcode/vx_statistics/ens_stats.cc +++ b/src/libcode/vx_statistics/ens_stats.cc @@ -177,7 +177,7 @@ void ECNTInfo::clear() { othresh.clear(); n_ens = n_pair = 0; - crps_emp = crpscl_emp = crpss_emp = bad_data_double; + crps_emp = crpscl_emp = crpss_emp = crps_emp_fair = bad_data_double; crps_gaus = crpscl_gaus = crpss_gaus = bad_data_double; ign = bad_data_double; me = rmse = spread = bad_data_double; @@ -199,7 +199,8 @@ void ECNTInfo::assign(const ECNTInfo &c) { crps_emp = c.crps_emp; crpscl_emp = c.crpscl_emp; crpss_emp = c.crpss_emp; - + crps_emp_fair = c.crps_emp_fair; + crps_gaus = c.crps_gaus; crpscl_gaus = c.crpscl_gaus; crpss_gaus = c.crpss_gaus; @@ -230,6 +231,8 @@ void ECNTInfo::set(const PairDataEnsemble &pd) { // Compute empirical CRPS scores crps_emp = pd.crps_emp_na.wmean(pd.wgt_na); + // Stub in for now? + crps_emp_fair = pd.crps_emp_fair_na.wmean(pd.wgt_na); crpscl_emp = pd.crpscl_emp_na.wmean(pd.wgt_na); crpss_emp = (is_bad_data(crps_emp) || is_bad_data(crpscl_emp) ||