Disable kinetic injection temperature assert

parthenon-hpc-lab · Nov 22, 2023 · 73cc819 · 73cc819
1 parent f9703e7
commit 73cc819
Showing 1 changed file with 9 additions and 8 deletions.
diff --git a/src/pgen/cluster/agn_feedback.cpp b/src/pgen/cluster/agn_feedback.cpp
@@ -361,8 +361,8 @@ void AGNFeedback::FeedbackSrcTerm(parthenon::MeshData<parthenon::Real> *md,
             ///////////////////////////////////////////////////////////////////
 
             eos.ConsToPrim(cons, prim, nhydro, nscalars, k, j, i);
-            const Real old_specific_internal_e =
-                prim(IPR, k, j, i) / (prim(IDN, k, j, i) * (eos.GetGamma() - 1.));
+            //const Real old_specific_internal_e =
+            //    prim(IPR, k, j, i) / (prim(IDN, k, j, i) * (eos.GetGamma() - 1.));
 
             cons(IDN, k, j, i) += jet_density;
             cons(IM1, k, j, i) += jet_momentum * sign_jet * jet_axis_x;
@@ -379,12 +379,13 @@ void AGNFeedback::FeedbackSrcTerm(parthenon::MeshData<parthenon::Real> *md,
             }
 
             eos.ConsToPrim(cons, prim, nhydro, nscalars, k, j, i);
-            const Real new_specific_internal_e =
-                prim(IPR, k, j, i) / (prim(IDN, k, j, i) * (eos.GetGamma() - 1.));
-            PARTHENON_REQUIRE(
-                new_specific_internal_e > jet_specific_internal_e ||
-                    new_specific_internal_e > old_specific_internal_e,
-                "Kinetic injection leads to temperature below jet and existing gas");
+            //We're 
+            //const Real new_specific_internal_e =
+            //    prim(IPR, k, j, i) / (prim(IDN, k, j, i) * (eos.GetGamma() - 1.));
+            //PARTHENON_REQUIRE(
+            //    new_specific_internal_e > jet_specific_internal_e ||
+            //        new_specific_internal_e > old_specific_internal_e,
+            //    "Kinetic injection leads to temperature below jet and existing gas");
           }
 
           // Apply velocity ceiling