Removed DEBUGGING statements

src/hydro/srcterms/tabular_cooling.cpp

@@ -269,15 +269,6 @@ void TabularCooling::SubcyclingFixedIntSrcTerm(MeshData<Real> *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<Real> *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;
         }

src/pgen/cluster/agn_triggering.cpp

@@ -585,12 +585,6 @@ AGNTriggeringFinalizeTriggering(parthenon::MeshData<parthenon::Real> *md,
     }
   }
 
-  // DEBUGGING (forrestglines)
-  if (agn_triggering.GetAccretionRate(hydro_pkg.get()) > 0) {
-    printf("Accreting cold gas!\n");
-  }
-  // END DEBUGGING
-
   return TaskStatus::complete;
 }