From f327af31a2105b6051eeeb6c97b0317c91937e9f Mon Sep 17 00:00:00 2001 From: Michael Zingale Date: Mon, 20 May 2024 12:53:12 -0400 Subject: [PATCH] output the densities after we converge --- sub_chandra/init_1d.H | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/sub_chandra/init_1d.H b/sub_chandra/init_1d.H index f876d84..635b049 100644 --- a/sub_chandra/init_1d.H +++ b/sub_chandra/init_1d.H @@ -462,7 +462,6 @@ AMREX_INLINE void init_1d() { } // end loop over zones - mass_he = 0.0; mass_wd = 0.0; @@ -568,6 +567,11 @@ AMREX_INLINE void init_1d() { amrex::Error("ERROR: mass did not converge"); } + std::cout << "converged" << std::endl; + std::cout << "central density of WD: " << rho_c << std::endl; + std::cout << "density at base of He layer: " << rho_he << std::endl; + + std::cout << "final masses: " << std::endl; std::cout << " mass WD: " << mass_wd / C::M_solar << std::endl; std::cout << " mass He: " << mass_he / C::M_solar << std::endl;