diff --git a/inputs/precipitator_mhd_lowres.in b/inputs/precipitator_mhd_lowres.in
index 355d7710..0e35d9ef 100644
--- 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
index 14bd5418..4453d299 100644
--- 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);