From 8c629e620359fb99f5a3a87a6af73c0c1e5e5a7a Mon Sep 17 00:00:00 2001
From: John Halley Gotway <johnhg@ucar.edu>
Date: Tue, 3 Sep 2024 14:10:45 -0600
Subject: [PATCH] Per #2887, update ctc_by_row() logic to create reproducible
 results with the develop branch.

---
 src/libcode/vx_statistics/contable_nx2.cc | 11 +++++++----
 src/libcode/vx_statistics/met_stats.cc    |  2 +-
 2 files changed, 8 insertions(+), 5 deletions(-)

diff --git a/src/libcode/vx_statistics/contable_nx2.cc b/src/libcode/vx_statistics/contable_nx2.cc
index 42f4742c9b..55a9187259 100644
--- a/src/libcode/vx_statistics/contable_nx2.cc
+++ b/src/libcode/vx_statistics/contable_nx2.cc
@@ -543,10 +543,10 @@ TTContingencyTable Nx2ContingencyTable::ctc_by_row(int row) const {
       exit(1);
    }
 
-   double sy;
-   double sn;
+   double sy = 0.0;
+   double sn = 0.0;
 
-   for(int j=(row + 1), sy=sn=0.0; j<Nrows; ++j) {
+   for(int j=row+1; j<Nrows; ++j) {
       sy +=    event_count_by_row(j);
       sn += nonevent_count_by_row(j);
    }
@@ -554,7 +554,10 @@ TTContingencyTable Nx2ContingencyTable::ctc_by_row(int row) const {
    tt.set_fy_oy(sy);
    tt.set_fy_on(sn);
 
-   for(int j=0, sy=sn=0.0; j<=row; ++j) {
+   sy = 0.0;
+   sn = 0.0;
+
+   for(int j=0; j<=row; ++j) {
       sy +=    event_count_by_row(j);
       sn += nonevent_count_by_row(j);
    }
diff --git a/src/libcode/vx_statistics/met_stats.cc b/src/libcode/vx_statistics/met_stats.cc
index dc6b4e9419..bcdd37e23e 100644
--- a/src/libcode/vx_statistics/met_stats.cc
+++ b/src/libcode/vx_statistics/met_stats.cc
@@ -3470,7 +3470,7 @@ double PCTInfo::get_stat_pjc(const string &stat_name,
 
 double PCTInfo::get_stat_prc(const string &stat_name,
                              ConcatString &col_name) const {
-   int i = 0;
+   int i = bad_data_int;
    double v = bad_data_double;
    col_name = stat_name;
    TTContingencyTable ct;