From 4658fa09ebfe72b7600a132ebebaa42f577a8c80 Mon Sep 17 00:00:00 2001 From: mari2895 Date: Wed, 20 Mar 2024 20:56:28 -0400 Subject: [PATCH] typo --- src/analysis/history.hpp | 4 ++++ src/progenitor/progenitordata.cpp | 2 +- 2 files changed, 5 insertions(+), 1 deletion(-) diff --git a/src/analysis/history.hpp b/src/analysis/history.hpp index 55dfdae3..44491e59 100644 --- a/src/analysis/history.hpp +++ b/src/analysis/history.hpp @@ -123,6 +123,10 @@ Real ReduceInGain(MeshData *md, bool is_conserved, int idx = 0) { net_heat_threshold); // checks that in the gain region bool is_outside_pns = (v(b, fluid_prim::entropy(), k, j, i) > outside_pns_threshold); // checks that outside PNS + std::cout << "is_outside_pns=" << is_outside_pns + << "outside_pns_threshold=" << outside_pns_threshold << std::endl; + std::cout << "is_netheat=" << is_netheat + << "net_heat_threshold=" << net_heat_threshold << std::endl; const auto &coords = v.GetCoordinates(b); const Real vol = coords.CellVolume(k, j, i); if (is_conserved) { diff --git a/src/progenitor/progenitordata.cpp b/src/progenitor/progenitordata.cpp index 535a1424..39de2da3 100644 --- a/src/progenitor/progenitordata.cpp +++ b/src/progenitor/progenitordata.cpp @@ -90,7 +90,7 @@ std::shared_ptr Initialize(ParameterInput *pin) { 2.42e-5); // corresponds to entropy > 3 kb/baryon Real inside_pns_threshold = pin->GetOrAddReal("progenitor", "inside_pns_threshold", 0.008); // corresponds to r < 80 km - Real net_heat_threshold = pin->GetOrAddReal("progenitor", "inside_pns_threshold", + Real net_heat_threshold = pin->GetOrAddReal("progenitor", "net_heat_threshold", 1e-8); // corresponds to r < 80 km UnitConversions units(pin); Real LengthCGSToCode = units.GetLengthCGSToCode();