diff --git a/Source/Particles/PulsarParameters.cpp b/Source/Particles/PulsarParameters.cpp
index 1c940cee984..255c6ef8b0a 100644
--- a/Source/Particles/PulsarParameters.cpp
+++ b/Source/Particles/PulsarParameters.cpp
@@ -398,6 +398,7 @@ Pulsar::ReadParameters () {
         m_BLC = m_B_star/(S*S*S); // Bfield at LC
         m_re_scaled = m_beta_rec_RR * (3./2.) * PhysConst::q_e
                     / ( 4 * MathConst::pi * PhysConst::ep0 * m_gammarad_RR * m_gammarad_RR * PhysConst::c * m_BLC);
+        m_re_scaled = std::sqrt( m_re_scaled);
         amrex::Print() << " re_scaled at gammarad : " << m_gammarad_RR << " is : " << m_re_scaled << "\n";
         m_re_scaledratio = m_re_scaled / PhysConst::r_e;
         amrex::Print() << " ratio : re_scaled/re : " << m_re_scaledratio << "\n";