fix driving normalization

parthenon-hpc-lab · Mar 4, 2024 · 96ff9be · 96ff9be
1 parent 1ff3999
commit 96ff9be
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 2 deletions.
diff --git a/inputs/precipitator_lowres.in b/inputs/precipitator_lowres.in
@@ -107,6 +107,8 @@ perturb_exponent = 1    # P(k) ~ k^{-2}
 #sigma_v = 0
 sigma_v = 0.0001   # kpc Myr^{-1}
 k_peak = 2
+#num_modes = 40
+#vertical_driving_only = false
 num_modes = 12
 vertical_driving_only = true
 sol_weight = 1.0

diff --git a/src/pgen/precipitator.cpp b/src/pgen/precipitator.cpp
@@ -331,8 +331,8 @@ void TurbSrcTerm(MeshData<Real> *md, const parthenon::SimTime /*time*/, const Re
           Real dv_x = 0;
           Real dv_y = 0;
           if (!vertical_driving_only) {
-            dv_x = perturb_pack(b, 0, k, j, i);
-            dv_y = perturb_pack(b, 1, k, j, i);
+            dv_x = perturb_pack(b, 0, k, j, i) / v_norm;
+            dv_y = perturb_pack(b, 1, k, j, i) / v_norm;
           }
           const Real dv_z = perturb_pack(b, 2, k, j, i) / v_norm;