diff --git a/Exec/science/wdmerger/wdmerger_util.cpp b/Exec/science/wdmerger/wdmerger_util.cpp index 1e2b56d505..65b09f89e5 100644 --- a/Exec/science/wdmerger/wdmerger_util.cpp +++ b/Exec/science/wdmerger/wdmerger_util.cpp @@ -97,8 +97,8 @@ void set_wd_composition (Real mass, Real& envelope_mass, Real core_comp[NumSpec] } core_comp[iHe4] = 1.0_rt; - - amrex::Print() << "Created a pure He " << star_type << "." << std::endl; + + amrex::Print() << "Created a pure He " << star_type << "." << std::endl; for (int n = 0; n < NumSpec; ++n) { envelope_comp[n] = core_comp[n]; @@ -117,11 +117,11 @@ void set_wd_composition (Real mass, Real& envelope_mass, Real core_comp[NumSpec] core_comp[iC12] = problem::hybrid_wd_c_frac; core_comp[iO16] = problem::hybrid_wd_o_frac; - envelope_mass = problem::hybrid_wd_he_shell_mass; + envelope_mass = problem::hybrid_wd_he_shell_mass; + + amrex::Print()<< "Creating " << star_type << "with CO core with mass fractions C = "<< problem::hybrid_wd_c_frac << "and O = " + << problem::hybrid_wd_o_frac <<" and a He shell of solar mass =" << problem::hybrid_wd_he_shell_mass << "." << std::endl; - amrex::Print()<< "Creating " << star_type << "with CO core with mass fractions C = "<< problem::hybrid_wd_c_frac << "and O = " - << problem::hybrid_wd_o_frac <<" and a He shell of solar mass =" << problem::hybrid_wd_he_shell_mass << "." << std::endl; - if (envelope_mass > 0.0_rt) { if (iHe4 < 0) { @@ -180,7 +180,7 @@ void set_wd_composition (Real mass, Real& envelope_mass, Real core_comp[NumSpec] core_comp[iO16] = problem::onemg_wd_o_frac; core_comp[iNe20] = problem::onemg_wd_ne_frac; core_comp[iMg24] = problem::onemg_wd_mg_frac; - + amrex::Print()<<"Creating an ONeMg " << star_type << "." < 0.0_rt) {