From 2ff803f52bc3afc42ccec5b0ac85ddd15c91e597 Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Fri, 29 Dec 2023 14:50:54 -0500 Subject: [PATCH] add a comment --- massive_star/init_1d.H | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/massive_star/init_1d.H b/massive_star/init_1d.H index 0158abc..5dbc716 100644 --- 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; }