From 8b789409ffe37d60164515e0300c9b2afa5b3709 Mon Sep 17 00:00:00 2001 From: Forrest Glines Date: Tue, 19 Jul 2022 19:09:19 -0400 Subject: [PATCH] Removed DEBUGGING statements --- src/hydro/srcterms/tabular_cooling.cpp | 19 ------------------- src/pgen/cluster/agn_triggering.cpp | 6 ------ 2 files changed, 25 deletions(-) diff --git a/src/hydro/srcterms/tabular_cooling.cpp b/src/hydro/srcterms/tabular_cooling.cpp index 91eaa081..72696266 100644 --- a/src/hydro/srcterms/tabular_cooling.cpp +++ b/src/hydro/srcterms/tabular_cooling.cpp @@ -269,15 +269,6 @@ void TabularCooling::SubcyclingFixedIntSrcTerm(MeshData *md, const Real dt // TODO(pgrete) with potentially more EOS, a separate get_pressure (or similar) // function could be useful. - // DEBUGGING(forrestglines) please remove - Real total_e = cons(IEN, k, j, i); - Real kinetic_e = 0.5 * (SQR(cons(IM1, k, j, i)) + SQR(cons(IM2, k, j, i)) + - SQR(cons(IM3, k, j, i)) / rho); - Real magnetic_e = 0.5 * (SQR(cons(IB1, k, j, i)) + SQR(cons(IB2, k, j, i)) + - SQR(cons(IB3, k, j, i))); - Real pres = prim(IPR, k, j, i); - // END DEBUGGING - Real internal_e = cons(IEN, k, j, i) - 0.5 * (SQR(cons(IM1, k, j, i)) + SQR(cons(IM2, k, j, i)) + @@ -310,16 +301,6 @@ void TabularCooling::SubcyclingFixedIntSrcTerm(MeshData *md, const Real dt // temperature is already below floor. const Real dedt_initial = DeDt_wrapper(0.0, internal_e_initial, dedt_valid); - // DEBUGGING(forrestglines) please remove - // Testing if current cooling rate would drop the internal energy below the floor - // if ( internal_e_initial + dt*dedt_initial < 0 ) { - // printf("Fast Cooling\n"); - // } - // if ( dedt_initial == 0 ) { - // printf("Zero initial cooling\n"); - // } - // END DEBUGGING - if (dedt_initial == 0.0 || internal_e_initial < internal_e_floor) { return; } diff --git a/src/pgen/cluster/agn_triggering.cpp b/src/pgen/cluster/agn_triggering.cpp index cbf810e7..be82c4e4 100644 --- a/src/pgen/cluster/agn_triggering.cpp +++ b/src/pgen/cluster/agn_triggering.cpp @@ -585,12 +585,6 @@ AGNTriggeringFinalizeTriggering(parthenon::MeshData *md, } } - // DEBUGGING (forrestglines) - if (agn_triggering.GetAccretionRate(hydro_pkg.get()) > 0) { - printf("Accreting cold gas!\n"); - } - // END DEBUGGING - return TaskStatus::complete; }