typo

lanl · Mar 21, 2024 · 4658fa0 · 4658fa0
1 parent 7f3bc17
commit 4658fa0
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 1 deletion.
diff --git a/src/analysis/history.hpp b/src/analysis/history.hpp
@@ -123,6 +123,10 @@ Real ReduceInGain(MeshData<Real> *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
@@ -90,7 +90,7 @@ std::shared_ptr<StateDescriptor> 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();