1 files changed, 0 insertions, 109 deletions
diff --git a/Exec/test_particle_injection/ParticleProb.cpp b/Exec/test_particle_injection/ParticleProb.cpp
deleted file mode 100644
index 7cdee1d46..000000000
--- a/Exec/test_particle_injection/ParticleProb.cpp
+++ /dev/null
@@ -1,109 +0,0 @@
-
-//
-// Each problem must have its own version of PhysicalParticleContainer::InitData()
-// to initialize the particle data.  It must also initialize charge and mass.
-//
-
-#include <cmath>
-
-#include <AMReX_BLProfiler.H>
-
-#include <ParticleContainer.H>
-#include <WarpXConst.H>
-
-using namespace amrex;
-
-void
-PhysicalParticleContainer::InitData()
-{
-    BL_PROFILE("PPC::InitData()");
-
-    charge = -PhysConst::q_e;
-    mass = PhysConst::m_e;
-
-    const int lev = 0;
-
-    const Geometry& geom = Geom(lev);
-    const Real* dx  = geom.CellSize();
-
-    Real weight, ux, uy, uz;
-    Real particle_xmin, particle_xmax, particle_ymin, particle_ymax, particle_zmin, particle_zmax;
-    int n_part_per_cell;
-    {
-      ParmParse pp("langmuirwave");
-      n_part_per_cell = 1;
-      pp.query("num_particles_per_cell", n_part_per_cell);
-      weight = 1.e25;
-      pp.query("n_e", weight);
-      #if BL_SPACEDIM==3
-      weight *= dx[0]*dx[1]*dx[2]/n_part_per_cell;
-      #elif BL_SPACEDIM==2
-      weight *= dx[0]*dx[1]/n_part_per_cell;
-      #endif
-
-      particle_xmin = particle_ymin = particle_zmin = -2.e-5;
-      particle_xmax = particle_ymax = particle_zmax =  2.e-5;
-      pp.query("particle_xmin", particle_xmin);
-      pp.query("particle_xmax", particle_xmax);
-      pp.query("particle_ymin", particle_ymin);
-      pp.query("particle_ymax", particle_ymax);
-      pp.query("particle_zmin", particle_zmin);
-      pp.query("particle_zmax", particle_zmax);   
- 
-      ux = 0.;
-      uy = 0.;
-      uz = 0.;
-      pp.query("ux", ux);
-      pp.query("uy", uy);
-      pp.query("uz", uz);
-
-      ux *= PhysConst::c;
-      uy *= PhysConst::c;      
-      uz *= PhysConst::c;
-    }
-
-    std::array<Real,PIdx::nattribs> attribs;
-    attribs.fill(0.0);
-    attribs[PIdx::w ] = weight;
-    attribs[PIdx::ux] = ux;
-    attribs[PIdx::uy] = uy;
-    attribs[PIdx::uz] = uz;
-    
-    for (MFIter mfi = MakeMFIter(lev); mfi.isValid(); ++mfi)
-    {
-        const Box& tile_box = mfi.tilebox();
-        RealBox tile_real_box { tile_box, dx, geom.ProbLo() };
-
-        const int grid_id = mfi.index();
-        const int tile_id = mfi.LocalTileIndex();
-
-        const IntVect& boxlo = tile_box.smallEnd();
-        for (IntVect iv = boxlo; iv <= tile_box.bigEnd(); tile_box.next(iv))
-        {
-            for (int i_part=0; i_part<n_part_per_cell;i_part++)
-            {
-		Real particle_shift = (0.5+i_part)/n_part_per_cell;
-#if (BL_SPACEDIM == 3)
-                Real x = tile_real_box.lo(0) + (iv[0]-boxlo[0] + particle_shift)*dx[0];
-                Real y = tile_real_box.lo(1) + (iv[1]-boxlo[1] + particle_shift)*dx[1];
-                Real z = tile_real_box.lo(2) + (iv[2]-boxlo[2] + particle_shift)*dx[2];
-#elif (BL_SPACEDIM == 2)
-                Real x = tile_real_box.lo(0) + (iv[0]-boxlo[0] + particle_shift)*dx[0];
-                Real y = 0.0;
-                Real z = tile_real_box.lo(1) + (iv[1]-boxlo[1] + particle_shift)*dx[1];
-#endif   
-
-		if (x >= particle_xmax || x < particle_xmin ||
-		    y >= particle_ymax || y < particle_ymin ||
-		    z >= particle_zmax || z < particle_zmin )
-                {
-                    continue;
-                }
-                else
-                {
-                    AddOneParticle(lev, grid_id, tile_id, x, y, z, attribs);
-                }
-	    }
-        } 
-    }
-}