Now reading

Author: n8xm

Source/Driver/Nyx.H

@@ -996,6 +996,11 @@ protected:
   // Seed mass of SMBH
   //
   static amrex::Real mass_seed;
+
+  static int min_halo_n_cells;
+  static bool halo_negate;
+  static amrex::Real halo_component_threshold;
+
 #endif
 
   // Previous maximum number of steps for sundials

Source/Driver/Nyx.cpp

@@ -204,6 +204,9 @@ Real Nyx::he_species        = 0.24;
 #ifdef REEBER
 Real Nyx::mass_halo_min     = 1.e10;
 Real Nyx::mass_seed         = 1.e5;
+Real Nyx::halo_component_threshold = 81.66;
+int Nyx::min_halo_n_cells = 10;
+bool Nyx::halo_negate = true;
 #endif
 
 #ifdef _OPENMP
@@ -445,6 +448,11 @@ Nyx::read_params ()
 #ifdef REEBER
     pp_nyx.query("mass_halo_min", mass_halo_min);
     pp_nyx.query("mass_seed", mass_seed);
+
+    ParmParse pp_reeber("reeber");
+    pp_reeber.query("min_halo_n_cells", min_halo_n_cells);
+    pp_reeber.query("halo_component_threshold", halo_component_threshold);
+    pp_reeber.query("negate", halo_negate);
 #endif
 }
 

Util/reeber/src/reeber.cpp

@@ -598,10 +598,7 @@ void Nyx::runReeberAnalysis(Vector<MultiFab*>& new_state,
 
     diy::DiscreteBounds diy_domain(3);
 
-    // TODO: take rho, min_halo_n_cells as parameters
-    Real min_halo_n_cells = 10;
-    Real rho = 81.66;
-
+    Real rho = halo_component_threshold;
     Real absolute_rho = (Nyx::average_dm_density + Nyx::average_gas_density) * rho;
     bool negate = true;  // sweep superlevel sets, highest density = root
     int finest_level = parent->finestLevel();