Per #1735, since we are writing multiple output line types to the sam…

…e .stat file, adding stat_row counter to the job class. This way all functions that write to it can increment that row counter when needed.

met/src/libcode/vx_analysis_util/stat_job.cc

@@ -186,6 +186,8 @@ void STATAnalysisJob::clear() {
    if(dump_row)  { delete [] dump_row;  dump_row  = (char *) 0; }
    if(stat_file) { delete [] stat_file; stat_file = (char *) 0; }
 
+   stat_row = 0;
+
    out_line_type.clear();
 
    out_fcst_thresh.clear();
@@ -307,8 +309,8 @@ void STATAnalysisJob::assign(const STATAnalysisJob & aj) {
    wmo_fisher_stats     = aj.wmo_fisher_stats;
 
    column_thresh_map    = aj.column_thresh_map;
-   column_str_inc_map       = aj.column_str_inc_map;
-   column_str_exc_map       = aj.column_str_exc_map;
+   column_str_inc_map   = aj.column_str_inc_map;
+   column_str_exc_map   = aj.column_str_exc_map;
 
    by_column            = aj.by_column;
 
@@ -354,6 +356,8 @@ void STATAnalysisJob::assign(const STATAnalysisJob & aj) {
    mask_poly            = aj.mask_poly;
    mask_sid             = aj.mask_sid;
 
+   stat_row             = aj.stat_row;
+
    set_dump_row (aj.dump_row);
    set_stat_file(aj.stat_file);
 
@@ -1880,6 +1884,8 @@ void STATAnalysisJob::open_stat_file() {
 
    close_stat_file();
 
+   stat_row = 0;
+
    if(!stat_file) return;
 
    stat_out = new ofstream;
@@ -1904,9 +1910,10 @@ void STATAnalysisJob::setup_stat_file(int n_row, int n) {
    int i, c, n_col;
 
    //
-   // Nothing to do unless output STAT file stream is defined
+   // Nothing to do if no output STAT file stream is defined
+   // or if the output table has already been setup
    //
-   if(!stat_out) return;
+   if(!stat_out || stat_at.nrows() > 0) return;
 
    //
    // Check for a single output line type
@@ -1916,6 +1923,11 @@ void STATAnalysisJob::setup_stat_file(int n_row, int n) {
    out_lt = (out_sa.n() == 1 ?
              string_to_statlinetype(out_sa[0].c_str()) : no_stat_line_type);
 
+   //
+   // Multiply number of rows by the number of output line types
+   //
+   n_row *= out_sa.n();
+
    //
    // Loop through the output line types and determine the number of
    // output columns
@@ -2029,6 +2041,11 @@ void STATAnalysisJob::setup_stat_file(int n_row, int n) {
          break;
    }
 
+   //
+   // Initialize the stat row counter after writing the header line
+   //
+   stat_row = 1;
+
    return;
 }
 

met/src/libcode/vx_analysis_util/stat_job.h

@@ -249,6 +249,7 @@ class STATAnalysisJob {
       char        *stat_file; // dump output statistics to a STAT file
       ofstream    *stat_out;  // output file stream for -out_stat
       AsciiTable  stat_at;    // AsciiTable for buffering output STAT data
+      int         stat_row;   // Counter for the current stat row
 
       StringArray  out_line_type;        // output line types
       ThreshArray  out_fcst_thresh;      // output forecast threshold(s)