From f1b1164bfd4230f0ffa32cef2a326794163f8f55 Mon Sep 17 00:00:00 2001 From: RevathiJambunathan Date: Wed, 25 Oct 2023 11:49:07 -0700 Subject: [PATCH] add flag to follow Bfield --- Source/Particles/PulsarParameters.H | 1 + Source/Particles/PulsarParameters.cpp | 3 +++ 2 files changed, 4 insertions(+) diff --git a/Source/Particles/PulsarParameters.H b/Source/Particles/PulsarParameters.H index cda38107a66..86a03f62d9c 100644 --- a/Source/Particles/PulsarParameters.H +++ b/Source/Particles/PulsarParameters.H @@ -713,6 +713,7 @@ public: static amrex::Real m_BLC; static amrex::Real m_beta_rec_RR; static amrex::Real m_re_scaledratio; + static int m_do_zero_uperpB_driftframe; private: }; diff --git a/Source/Particles/PulsarParameters.cpp b/Source/Particles/PulsarParameters.cpp index 4dc0945f454..1ea378ee3cc 100644 --- a/Source/Particles/PulsarParameters.cpp +++ b/Source/Particles/PulsarParameters.cpp @@ -153,6 +153,7 @@ amrex::Real Pulsar::m_gammarad_RR; amrex::Real Pulsar::m_BLC; amrex::Real Pulsar::m_beta_rec_RR; amrex::Real Pulsar::m_re_scaledratio = 1.; +int Pulsar::m_do_zero_uperpB_driftframe; Pulsar::Pulsar () { @@ -403,6 +404,8 @@ Pulsar::ReadParameters () { m_re_scaledratio = m_re_scaled / PhysConst::r_e; amrex::Print() << " ratio : re_scaled/re : " << m_re_scaledratio << "\n"; } + pp.get("do_zero_uperpB_driftframe",m_do_zero_uperpB_driftframe); + amrex::Print() << " do zero uperp B " << m_do_zero_uperpB_driftframe << "\n"; }