Fixed missing magnetic energy from temperature ceil

parthenon-hpc-lab · Jul 17, 2023 · 828288d · 828288d
1 parent e695577
commit 828288d
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/eos/adiabatic_glmmhd.hpp b/src/eos/adiabatic_glmmhd.hpp
@@ -156,7 +156,7 @@ class AdiabaticGLMMHDEOS : public EquationOfState {
     const Real eff_pressure_ceiling = gm1 * u_d * e_ceiling_;
     if (w_p > eff_pressure_ceiling) {
       // apply temperature ceiling, correct total energy
-      u_e = (u_d * e_ceiling_) + e_k;
+      u_e = (u_d * e_ceiling_) + e_k + e_B;
       w_p = eff_pressure_ceiling;
     }