check added

lanl · Mar 13, 2024 · 9146f83 · 9146f83
1 parent 1a32819
commit 9146f83
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/src/analysis/history.hpp b/src/analysis/history.hpp
@@ -98,6 +98,8 @@ Real ReduceInGain(MeshData<Real> *md, int idx = 0) {
   static auto desc =
       MakePackDescriptor<Varname, iv::GcovHeat, iv::GcovCool>(resolved_pkgs.get());
   auto v = desc.GetPack(md);
+  PARTHENON_REQUIRE_THROWS(v.ContainsHost(0, iv::GcovHeat(), iv::GcovCool()),
+                           "Must be doing SN simulation");
   const int nblocks = v.GetNBlocks();
 
   Real result = 0.0;

diff --git a/src/fluid/fluid.cpp b/src/fluid/fluid.cpp
@@ -208,8 +208,8 @@ std::shared_ptr<StateDescriptor> Initialize(ParameterInput *pin) {
   physics->AddField(p::entropy::name(), mprim_scalar);
   physics->AddField(p::cs::name(), mprim_scalar);
   physics->AddField(diag::ratio_divv_cs::name(), mprim_scalar);
-  // physics->AddField(diag::localization_function::name(), mprim_scalar);  // (MG)
-  // currently not in use, turn on when needed.
+  // physics->AddField(diag::localization_function::name(), mprim_scalar);
+  // (MG) currently not in use, turn on when needed.
   physics->AddField(diag::entropy_z_0::name(), mprim_scalar);
   physics->AddField(p::gamma1::name(), mprim_scalar);
   if (ye) {