diff --git a/src/pgen/cluster/agn_feedback.cpp b/src/pgen/cluster/agn_feedback.cpp index 40bc0123..7523f663 100644 --- a/src/pgen/cluster/agn_feedback.cpp +++ b/src/pgen/cluster/agn_feedback.cpp @@ -361,8 +361,6 @@ 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.)); cons(IDN, k, j, i) += jet_density; cons(IM1, k, j, i) += jet_momentum * sign_jet * jet_axis_x; @@ -379,13 +377,6 @@ void AGNFeedback::FeedbackSrcTerm(parthenon::MeshData *md, } eos.ConsToPrim(cons, prim, nhydro, nscalars, k, j, i); - // 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