take sqrt

RevathiJambunathan · Oct 23, 2023 · ff4182b · ff4182b
1 parent ab64254
commit ff4182b
Showing 1 changed file with 1 addition and 0 deletions.
diff --git 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";