From 6acf85b9b00d68b07cae4ef685fababd0091e6b0 Mon Sep 17 00:00:00 2001
From: Michael Zingale <michael.zingale@stonybrook.edu>
Date: Sun, 22 Oct 2023 15:07:53 -0400
Subject: [PATCH] fix compilation

---
 integration/nse_update_simplified_sdc.H | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/integration/nse_update_simplified_sdc.H b/integration/nse_update_simplified_sdc.H
index a6b1c285ad..0a424d5dba 100644
--- a/integration/nse_update_simplified_sdc.H
+++ b/integration/nse_update_simplified_sdc.H
@@ -279,6 +279,8 @@ void sdc_nse_burn(BurnT& state, const Real dt) {
         rhoX_new[n] = state.y[SFS+n] + dt * state.ydot_a[SFS+n] + dt * rhoX_source[n];
     }
 
+    burn_t nse_state;
+
     for (int iter = 0; iter < 3; iter++) {
 
         // update (rho e)^{n+1} based on the new energy generation rate
@@ -316,18 +318,19 @@ void sdc_nse_burn(BurnT& state, const Real dt) {
             rhoX_tilde[n] = state.y[SRHO] * nse_state.xn[n] - dt * state.ydot_a[SFS+n];
         }
 
-        dyedt = 0.0_rt;   // we can update this in the future by calling actual_rhs()
+        Real dyedt = 0.0_rt;   // we can update this in the future by calling actual_rhs()
 
         // we want to compute (rho eps) = - N_A c^2 sum{m_i (rhoX_tilde - rhoX_old) / A_i}
         rho_enucdot = 0.0;
         for (int n = 0; n < NumSpec; ++n) {
             rho_enucdot += (rhoX_tilde[n] - rho_old * X_old[n]) *
-                network::mion(n+1) * aion_inv(n+1);
+                network::mion(n+1) * aion_inv[n];
         }
         rho_enucdot *= C::Legacy::enuc_conv2;
 
         // now get the updated neutrino term
 
+
         // update the new state for the next pass -- this should
         // already implicitly have the advective portion included,
         // since it was there when we called the NSE state
@@ -343,6 +346,7 @@ void sdc_nse_burn(BurnT& state, const Real dt) {
     // the new mass fractions are just those that come from the table
     // make sure they are normalized
     Real sum_X{0.0_rt};
+    Real X[NumSpec] = {0.0_rt};
     for (int n = 0; n < NumSpec; ++n) {
         X[n] = amrex::max(small_x, amrex::min(1.0_rt, nse_state.xn[n]));
         sum_X += X[n];