output the densities after we converge (#69)

AMReX-Astro · May 20, 2024 · 15618e9 · 15618e9
1 parent 69554e1
commit 15618e9
Showing 1 changed file with 5 additions and 1 deletion.
diff --git 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;