diff --git a/src/units.hpp b/src/units.hpp
index 3290e87e..7579354d 100644
--- 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