fix param bug

parthenon-hpc-lab · Aug 30, 2024 · 494e24c · 494e24c
1 parent 7d29aac
commit 494e24c
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/src/pgen/precipitator.cpp b/src/pgen/precipitator.cpp
@@ -287,14 +287,15 @@ void TurbSrcTerm(MeshData<Real> *md, const parthenon::SimTime /*time*/, const Re
   IndexRange kb = pmb->cellbounds.GetBoundsK(IndexDomain::interior);
 
   const auto sigma_v = hydro_pkg->Param<Real>("sigma_v");
-  const auto vertical_driving_only = hydro_pkg->Param<bool>("vertical_driving_only");
-
   const Real h_smooth = hydro_pkg->Param<Real>("h_smooth_heatcool");
 
   // the maximum height at which to drive turbulence
   const Real max_height_driving = hydro_pkg->Param<Real>("max_height_driving");
 
   if (sigma_v > 0) {
+    // NOTE: this param only exists if sigma_v > 0
+    const auto vertical_driving_only = hydro_pkg->Param<bool>("vertical_driving_only");
+
     // generate perturbations
     auto *few_modes_ft = hydro_pkg->MutableParam<FewModesFT>("precipitator/few_modes_ft_v");
     few_modes_ft->Generate(md, dt, "tmp_perturb");