fix plasma beta

parthenon-hpc-lab · Oct 10, 2024 · 358ae04 · 358ae04
1 parent ce1eb09
commit 358ae04
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/inputs/precipitator_mhd_lowres.in b/inputs/precipitator_mhd_lowres.in
@@ -83,7 +83,7 @@ cfl = 1.0
 
 <precipitator>
 hse_profile_filename = inputs/profile_tc_tff_3.0.txt
-bfield_gauss = 1.0e-8           # 10 nG
+bfield_gauss = 1.0e-6           # 1.0 uG
 uniform_init = 0
 enable_heating = magic          # options: [magic, none]
 thermostat_temperature = 1.0e7  # K

diff --git a/src/pgen/precipitator.cpp b/src/pgen/precipitator.cpp
@@ -1279,7 +1279,7 @@ void UserMeshWorkBeforeOutput(Mesh *mesh, ParameterInput *pin,
           const Real by = prim(IB2, k, j, i);
           const Real bz = prim(IB3, k, j, i);
           const Real P_mag = 0.5 * (SQR(bx) + SQR(by) + SQR(bz));
-          const Real beta = P_mag / P;
+          const Real beta = P / P_mag;
 
           const Real v1 = prim(IV1, k, j, i);
           const Real v2 = prim(IV2, k, j, i);