diff --git a/inputs/precipitator_lowres.in b/inputs/precipitator_lowres.in index 0d82a203..482d12e2 100644 --- 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 index 2f248391..cda3b19e 100644 --- a/src/pgen/precipitator.cpp +++ b/src/pgen/precipitator.cpp @@ -331,8 +331,8 @@ void TurbSrcTerm(MeshData *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;