diff --git a/src/pgen/cluster/agn_feedback.cpp b/src/pgen/cluster/agn_feedback.cpp index bcf63c24..7ef1f848 100644 --- a/src/pgen/cluster/agn_feedback.cpp +++ b/src/pgen/cluster/agn_feedback.cpp @@ -361,8 +361,8 @@ void AGNFeedback::FeedbackSrcTerm(parthenon::MeshData *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 *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