add start time for RR

RevathiJambunathan · Dec 20, 2023 · 5da05fb · 5da05fb
1 parent 8ff3c78
commit 5da05fb
Show file tree

Hide file tree

Showing 3 changed files with 7 additions and 2 deletions.
diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp
@@ -2854,7 +2854,8 @@ PhysicalParticleContainer::PushP (int lev, Real dt,
     amrex::Real BC_width = Pulsar::m_BC_width;
     amrex::Real gather_buffer_boxmin = Pulsar::m_gatherbuffer_min;
     amrex::Real gather_buffer_boxmax = Pulsar::m_gatherbuffer_max;
-    amrex::Real re_scaledratio = Pulsar::m_re_scaledratio;
+    amrex::Real re_scaledratio = 1.;
+    if (cur_time > Pulsar::m_RR_start_time) { re_scaledratio = Pulsar::m_re_scaledratio;}
     int do_zero_uperpB_driftframe = Pulsar::m_do_zero_uperpB_driftframe;
     amrex::Real rmax_zero_uperp_driftframe = Pulsar::m_rmax_zero_uperpB_driftframe;
 #endif
@@ -3169,7 +3170,8 @@ PhysicalParticleContainer::PushPX (WarpXParIter& pti,
     amrex::Real BC_width = Pulsar::m_BC_width;
     amrex::Real gather_buffer_boxmin = Pulsar::m_gatherbuffer_min;
     amrex::Real gather_buffer_boxmax = Pulsar::m_gatherbuffer_max;
-    amrex::Real re_scaledratio = Pulsar::m_re_scaledratio;
+    amrex::Real re_scaledratio = 1.;
+    if (cur_time > Pulsar::m_RR_start_time) { re_scaledratio = Pulsar::m_re_scaledratio;}
     int do_zero_uperpB_driftframe = Pulsar::m_do_zero_uperpB_driftframe;
     amrex::Real rmax_zero_uperp_driftframe = Pulsar::m_rmax_zero_uperpB_driftframe;
 #endif

diff --git a/Source/Particles/PulsarParameters.H b/Source/Particles/PulsarParameters.H
@@ -715,6 +715,7 @@ public:
     static amrex::Real m_re_scaledratio;
     static int m_do_zero_uperpB_driftframe;
     static amrex::Real m_rmax_zero_uperpB_driftframe;
+    static amrex::Real m_RR_start_time;
 private:
 };
 

diff --git a/Source/Particles/PulsarParameters.cpp b/Source/Particles/PulsarParameters.cpp
@@ -155,6 +155,7 @@ amrex::Real Pulsar::m_beta_rec_RR;
 amrex::Real Pulsar::m_re_scaledratio = 1.;
 int Pulsar::m_do_zero_uperpB_driftframe = 0;
 amrex::Real Pulsar::m_rmax_zero_uperpB_driftframe;
+amrex::Real Pulsar::m_RR_start_time;
 
 Pulsar::Pulsar ()
 {
@@ -408,6 +409,7 @@ Pulsar::ReadParameters () {
     pp.get("do_zero_uperpB_driftframe",m_do_zero_uperpB_driftframe);
     amrex::Print() << " do zero uperp B " << m_do_zero_uperpB_driftframe << "\n";
     pp.get("rmax_zero_uperpB_driftframe",m_rmax_zero_uperpB_driftframe);
+    pp.get("RR_start_time",m_RR_start_time);
 }