add a comment

AMReX-Astro · Dec 29, 2023 · 2ff803f · 2ff803f
1 parent ca00655
commit 2ff803f
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/massive_star/init_1d.H b/massive_star/init_1d.H
@@ -540,6 +540,11 @@ AMREX_INLINE void init_1d() {
 
                eos_state.aux[AuxZero::iye] = ye;
 
+	       // note: if we are in NSE, then as we adjust density in
+	       // this Newton loop, the composition will change, so
+	       // convergence should mean that the density and
+	       // composition have converged.
+
                set_aux(eos_state);
 
                eos(eos_input_rt, eos_state);
@@ -577,6 +582,7 @@ AMREX_INLINE void init_1d() {
                amrex::Error("Error: HSE non-convergence");
            }
 
+
            if (temp_zone < temp_fluff) {
                temp_zone = temp_fluff;
            }