Merged in plasma_injector (pull request #14)

Plasma injector

Approved-by: Weiqun Zhang <WeiqunZhang@lbl.gov>
Approved-by: Remi Lehe <remi.lehe@normalesup.org>

1 files changed, 3 insertions, 12 deletions

diff --git a/Source/WarpXEvolve.cpp b/Source/WarpXEvolve.cpp
index 16b957dec..29336e66b 100644
--- a/Source/WarpXEvolve.cpp
+++ b/Source/WarpXEvolve.cpp
@@ -284,20 +284,11 @@ WarpX::InjectPlasma (int num_shift, int dir)
 
         const Real* dx = geom[lev].CellSize();
 
-        for (int i = 0; i < num_injected_species; ++i)
-        {
-            int ppc = injected_plasma_ppc[i];
-            Real density = injected_plasma_density[i];
-#if BL_SPACEDIM==3
-            Real weight = density * dx[0]*dx[1]*dx[2]/ppc;
-#elif BL_SPACEDIM==2
-            Real weight = density * dx[0]*dx[1]/ppc;
-#endif
-
+        for (int i = 0; i < num_injected_species; ++i) {
             int ispecies = injected_plasma_species[i];
             WarpXParticleContainer& pc = mypc->GetParticleContainer(ispecies);
-
-            pc.AddParticles(lev, particleBox, weight, ppc);
+            auto& ppc = dynamic_cast<PhysicalParticleContainer&>(pc);
+            ppc.AddParticles(lev, particleBox);
         }
     }
 }