cleanup hydro_driver

src/hydro/hydro_driver.cpp

@@ -8,22 +8,19 @@
 #include <memory>
 #include <string>
 #include <sys/types.h>
-#include <utility>
 #include <vector>
 
 // Parthenon headers
 #include "amr_criteria/refinement_package.hpp"
 #include "bvals/comms/bvals_in_one.hpp"
-#include "prolong_restrict/prolong_restrict.hpp"
 #include <parthenon/parthenon.hpp>
 // AthenaPK headers
 #include "../eos/adiabatic_hydro.hpp"
-//#include "../pgen/cluster/agn_triggering.hpp"
-//#include "../pgen/cluster/magnetic_tower.hpp"
+#include "../pgen/cluster/agn_triggering.hpp"
+#include "../pgen/cluster/magnetic_tower.hpp"
 #include "glmmhd/glmmhd.hpp"
 #include "hydro.hpp"
 #include "hydro_driver.hpp"
-#include "srcterms/tabular_cooling.hpp"
 #include "utils/error_checking.hpp"
 
 using namespace parthenon::driver::prelude;
@@ -92,7 +89,6 @@ TaskCollection HydroDriver::MakeTaskCollection(BlockList_t &blocks, int stage) {
   auto num_task_lists_executed_independently = blocks.size();
 
   const int num_partitions = pmesh->DefaultNumPartitions();
-#if 0
   // calculate agn triggering accretion rate
   if ((stage == 1) &&
       hydro_pkg->AllParams().hasKey("agn_triggering_reduce_accretion_rate") &&
@@ -128,7 +124,6 @@ TaskCollection HydroDriver::MakeTaskCollection(BlockList_t &blocks, int stage) {
       prev_task = new_remove_accreted_gas;
     }
   }
-#endif
 
   for (int i = 0; i < blocks.size(); i++) {
     auto &pmb = blocks[i];
@@ -191,7 +186,7 @@ TaskCollection HydroDriver::MakeTaskCollection(BlockList_t &blocks, int stage) {
         },
         hydro_pkg.get());
   }
-#if 0
+
   // calculate magnetic tower scaling
   if ((stage == 1) && hydro_pkg->AllParams().hasKey("magnetic_tower_power_scaling") &&
       hydro_pkg->Param<bool>("magnetic_tower_power_scaling")) {
@@ -233,7 +228,6 @@ TaskCollection HydroDriver::MakeTaskCollection(BlockList_t &blocks, int stage) {
         hydro_pkg.get());
 #endif
   }
-#endif
 
   // First add split sources before the main time integration
   if (stage == 1) {