add sqrt(4pi) to cgs B-field units

parthenon-hpc-lab · Nov 4, 2024 · b7394c0 · b7394c0
1 parent c7273f9
commit b7394c0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/units.hpp b/src/units.hpp
@@ -90,7 +90,7 @@ class Units {
     return code_energy_cgs() / kev_cgs * code_length_cgs() * code_length_cgs();
   }
   parthenon::Real code_magnetic_cgs() const {
-    return sqrt(code_mass_cgs()) / sqrt(code_length_cgs()) / code_time_cgs();
+    return std::sqrt(4.0 * M_PI) * sqrt(code_mass_cgs()) / sqrt(code_length_cgs()) / code_time_cgs();
   }
 
   // Physical Constants in code units