removing particles if inside R_Star

RevathiJambunathan · Jan 6, 2021 · 955a3e6 · 955a3e6
1 parent 93b0dfd
commit 955a3e6
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 5 deletions.
diff --git a/Source/Initialization/InjectorPosition.H b/Source/Initialization/InjectorPosition.H
@@ -139,8 +139,8 @@ struct InjectorPosition
     bool
     insidePulsarBounds (amrex::Real r, amrex::Real R_star, amrex::Real dR_star) const noexcept
     {
-    //    return (r>=R_star and r<=(R_star+dR_star));
-        return (r<=R_star and r>=(R_star-dR_star));
+        //return (r<=R_star and r>=(R_star-dR_star));
+        return (r<=(R_star+dR_star) and r>=(R_star-dR_star));
     }
 #endif
 

diff --git a/Source/Particles/PhysicalParticleContainer.cpp b/Source/Particles/PhysicalParticleContainer.cpp
@@ -934,8 +934,8 @@ PhysicalParticleContainer::AddPlasma (int lev, RealBox part_realbox)
                        }
                        //if (sigma_inj < 0 and q_pm >0) {p.id()=-1; continue;}
                        //if (sigma_inj > 0 and q_pm <0) {p.id()=-1; continue;}
-                       //if (sigma_inj < 0 and q_pm >0) {p.id()=-1; continue;}
-                       //if (sigma_inj > 0 and q_pm <0) {p.id()=-1; continue;}
+                       if (sigma_inj < 0 and q_pm >0) {p.id()=-1; continue;}
+                       if (sigma_inj > 0 and q_pm <0) {p.id()=-1; continue;}
                        // if rho is too smal -- we dont inject particles
                        if (std::abs(rho_GJ) < 1.0E-20) {
                           p.id() = -1;
@@ -2254,7 +2254,8 @@ void PhysicalParticleContainer::PulsarParticleRemoval() {
                       Real r = std::sqrt((x-xc)*(x-xc)
                                        + (y-yc)*(y-yc)
                                        + (z-zc)*(z-zc));
-                      if (r<=(PulsarParm::R_star-PulsarParm::dR_star)) {
+                      //if (r<=(PulsarParm::R_star-PulsarParm::dR_star)) {
+                      if (r<=(PulsarParm::R_star)) {
                           pp[i].id() = -1;
                       }
             });