From 92b42f972bdcae3067630c5f5bd9471e2bb21820 Mon Sep 17 00:00:00 2001
From: Michael Zingale <michael.zingale@stonybrook.edu>
Date: Wed, 17 Jul 2024 15:42:02 -0400
Subject: [PATCH] NSE now requires SDC (#72)

---
 massive_star/GNUmakefile |  2 ++
 massive_star/init_1d.H   | 11 +----------
 2 files changed, 3 insertions(+), 10 deletions(-)

diff --git a/massive_star/GNUmakefile b/massive_star/GNUmakefile
index e6690f6..eb75261 100644
--- a/massive_star/GNUmakefile
+++ b/massive_star/GNUmakefile
@@ -15,6 +15,8 @@ USE_REACT = TRUE
 USE_CXX_EOS = TRUE
 USE_CXX_REACTIONS = TRUE
 
+USE_SIMPLIFIED_SDC = TRUE
+
 USE_FORT_MICROPHYSICS = FALSE
 
 MAX_ZONES := 65536
diff --git a/massive_star/init_1d.H b/massive_star/init_1d.H
index 5dbc716..a84bf6b 100644
--- a/massive_star/init_1d.H
+++ b/massive_star/init_1d.H
@@ -63,14 +63,7 @@ constexpr Real temp_fluff = 1.e5_rt;
 AMREX_INLINE void
 set_aux(eos_t& eos_state) {
 
-    burn_t burn_state;
-    burn_state.rho = eos_state.rho;
-    burn_state.T = eos_state.T;
-    for (int n = 0; n < NumSpec; ++n) {
-        burn_state.xn[n] = eos_state.xn[n];
-    }
-
-    bool nse_check = in_nse(burn_state);
+    bool nse_check = in_nse(eos_state);
 
     if (nse_check) {
 
@@ -212,8 +205,6 @@ AMREX_INLINE void init_1d() {
 
     for (int i = 0; i < problem_rp::nx; ++i) {
 
-        // we need to fill a burn_t in order to check for NSE
-
         eos_state.rho = model_mesa_hse(i, model::idens);
         eos_state.T = model_mesa_hse(i, model::itemp);