NSE now requires SDC (#72)

AMReX-Astro · Jul 17, 2024 · 92b42f9 · 92b42f9
1 parent f4d3200
commit 92b42f9
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 10 deletions.
diff --git 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
@@ -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);