Per #1706, add bugfix to the develop branch. Also add a new job to un…

…it_stat_analysis.xml to test out the aggregation of the ECNT line type. This will add new unit test output and cause the NB to fail. (#1708)
dtcenter · Mar 10, 2021 · d80aafa · d80aafa
1 parent 592c937
commit d80aafa
Show file tree

Hide file tree

Showing 2 changed files with 17 additions and 0 deletions.
diff --git a/met/src/tools/core/stat_analysis/aggr_stat_line.cc b/met/src/tools/core/stat_analysis/aggr_stat_line.cc
@@ -2597,6 +2597,7 @@ void aggr_ecnt_lines(LineDataFile &f, STATAnalysisJob &job,
          m[key].ens_pd.var_oerr_na.add(square(cur.spread_oerr));
          m[key].ens_pd.var_plus_oerr_na.add(square(cur.spread_plus_oerr));
          m[key].ens_pd.wgt_na.add(cur.total);
+         m[key].ens_pd.skip_ba.add(false);
 
          //
          // Store the summary statistics

diff --git a/test/xml/unit_stat_analysis.xml b/test/xml/unit_stat_analysis.xml
@@ -74,6 +74,22 @@
     </output>
   </test>
 
+  <test name="stat_analysis_AGG_ECNT">
+    <exec>&MET_BIN;/stat_analysis</exec>
+    <param> \
+      -lookin &OUTPUT_DIR;/ensemble_stat/ensemble_stat_SKIP_CONST_20120410_120000V.stat \
+      -job aggregate -line_type ECNT -by FCST_VAR -obs_thresh NA -vx_mask NWC,GRB \
+      -dump_row &OUTPUT_DIR;/stat_analysis/AGG_ECNT_dump.stat \
+      -out_stat &OUTPUT_DIR;/stat_analysis/AGG_ECNT_out.stat \
+      -set_hdr VX_MASK NWC_AND_GRB \
+      -v 1
+    </param>
+    <output>
+      <stat>&OUTPUT_DIR;/stat_analysis/AGG_ECNT_dump.stat</stat>
+      <stat>&OUTPUT_DIR;/stat_analysis/AGG_ECNT_out.stat</stat>
+    </output>
+  </test>
+
   <test name="stat_analysis_AGG_STAT_ORANK_RHIST">
     <exec>&MET_BIN;/stat_analysis</exec>
     <param> \