diff options
Diffstat (limited to 'Source/Particles/RigidInjectedParticleContainer.cpp')
| -rw-r--r-- | Source/Particles/RigidInjectedParticleContainer.cpp | 11 |
1 files changed, 3 insertions, 8 deletions
diff --git a/Source/Particles/RigidInjectedParticleContainer.cpp b/Source/Particles/RigidInjectedParticleContainer.cpp index 00e67b5c2..71bddabdd 100644 --- a/Source/Particles/RigidInjectedParticleContainer.cpp +++ b/Source/Particles/RigidInjectedParticleContainer.cpp @@ -487,27 +487,22 @@ RigidInjectedParticleContainer::PushP (int lev, Real dt, getExternalE(ip, Exp, Eyp, Ezp); getExternalB(ip, Bxp, Byp, Bzp); + amrex::Real qp = q; + if (ion_lev) { qp *= ion_lev[ip]; } + if (do_crr) { - amrex::Real qp = q; - if (ion_lev) { qp *= ion_lev[ip]; } UpdateMomentumBorisWithRadiationReaction(uxpp[ip], uypp[ip], uzpp[ip], Exp, Eyp, Ezp, Bxp, Byp, Bzp, qp, m, dt); } else if (pusher_algo == ParticlePusherAlgo::Boris) { - amrex::Real qp = q; - if (ion_lev) { qp *= ion_lev[ip]; } UpdateMomentumBoris( uxpp[ip], uypp[ip], uzpp[ip], Exp, Eyp, Ezp, Bxp, Byp, Bzp, qp, m, dt); } else if (pusher_algo == ParticlePusherAlgo::Vay) { - amrex::Real qp = q; - if (ion_lev){ qp *= ion_lev[ip]; } UpdateMomentumVay( uxpp[ip], uypp[ip], uzpp[ip], Exp, Eyp, Ezp, Bxp, Byp, Bzp, qp, m, dt); } else if (pusher_algo == ParticlePusherAlgo::HigueraCary) { - amrex::Real qp = q; - if (ion_lev){ qp *= ion_lev[ip]; } UpdateMomentumHigueraCary( uxpp[ip], uypp[ip], uzpp[ip], Exp, Eyp, Ezp, Bxp, Byp, Bzp, qp, m, dt); |