Merge pull request #665 from RemiLehe/remove_fortran_push_position

Remove fortran push position in electrostatic code

1 files changed, 0 insertions, 110 deletions

diff --git a/Source/Particles/push_particles_ES.F90 b/Source/Particles/push_particles_ES.F90
index 60833d456..a22ee5a62 100644
--- a/Source/Particles/push_particles_ES.F90
+++ b/Source/Particles/push_particles_ES.F90
@@ -151,114 +151,4 @@ contains
 
   end subroutine warpx_push_leapfrog_2d
 
-
-!
-! This routine advances the particle positions using the current
-! velocity. This is needed to desynchronize the particle positions
-! from the velocities after particle initialization.
-!
-! Arguments:
-!     particles : a pointer to the particle array-of-structs
-!     ns        : the stride length of particle struct (the size of the struct in number of reals)
-!     np        : the number of particles
-!     xx_p      : the electric field in the x-direction at the particle positions
-!     vy_p      : the electric field in the y-direction at the particle positions
-!     vz_p      : the electric field in the z-direction at the particle positions
-!     dt        : the time step
-!     prob_lo   : the left-hand corner of the problem domain
-!     prob_hi   : the right-hand corner of the problem domain
-!
-  subroutine warpx_push_leapfrog_positions_3d(particles, ns, np,      &
-                                             vx_p, vy_p, vz_p, dt,   &
-                                             prob_lo, prob_hi)       &
-       bind(c,name='warpx_push_leapfrog_positions_3d')
-    integer, value,   intent(in)    :: ns, np
-    real(amrex_particle_real), intent(inout) :: particles(ns,np)
-    real(amrex_particle_real), intent(inout) :: vx_p(np), vy_p(np), vz_p(np)
-    real(amrex_real), intent(in)    :: dt
-    real(amrex_real), intent(in)    :: prob_lo(3), prob_hi(3)
-
-    integer n
-
-    do n = 1, np
-
-       particles(1, n) = particles(1, n) + dt * vx_p(n)
-       particles(2, n) = particles(2, n) + dt * vy_p(n)
-       particles(3, n) = particles(3, n) + dt * vz_p(n)
-
-!      bounce off the walls in the x...
-       do while (particles(1, n) .lt. prob_lo(1) .or. particles(1, n) .gt. prob_hi(1))
-          if (particles(1, n) .lt. prob_lo(1)) then
-             particles(1, n) = 2.d0*prob_lo(1) - particles(1, n)
-          else
-             particles(1, n) = 2.d0*prob_hi(1) - particles(1, n)
-          end if
-          vx_p(n) = -vx_p(n)
-       end do
-
-!      ... y...
-       do while (particles(2, n) .lt. prob_lo(2) .or. particles(2, n) .gt. prob_hi(2))
-          if (particles(2, n) .lt. prob_lo(2)) then
-             particles(2, n) = 2.d0*prob_lo(2) - particles(2, n)
-          else
-             particles(2, n) = 2.d0*prob_hi(2) - particles(2, n)
-          end if
-          vy_p(n) = -vy_p(n)
-       end do
-
-!      ... and z directions
-       do while (particles(3, n) .lt. prob_lo(3) .or. particles(3, n) .gt. prob_hi(3))
-          if (particles(3, n) .lt. prob_lo(3)) then
-             particles(3, n) = 2.d0*prob_lo(3) - particles(3, n)
-          else
-             particles(3, n) = 2.d0*prob_hi(3) - particles(3, n)
-          end if
-          vz_p(n) = -vz_p(n)
-       end do
-
-    end do
-
-  end subroutine warpx_push_leapfrog_positions_3d
-
-  subroutine warpx_push_leapfrog_positions_2d(particles, ns, np,      &
-                                              vx_p, vy_p, dt,         &
-                                              prob_lo, prob_hi)       &
-       bind(c,name='warpx_push_leapfrog_positions_2d')
-    integer, value,   intent(in)    :: ns, np
-    real(amrex_particle_real), intent(inout) :: particles(ns,np)
-    real(amrex_particle_real), intent(inout) :: vx_p(np), vy_p(np)
-    real(amrex_real), intent(in)    :: dt
-    real(amrex_real), intent(in)    :: prob_lo(2), prob_hi(2)
-
-    integer n
-
-    do n = 1, np
-
-       particles(1, n) = particles(1, n) + dt * vx_p(n)
-       particles(2, n) = particles(2, n) + dt * vy_p(n)
-
-!      bounce off the walls in the x...
-       do while (particles(1, n) .lt. prob_lo(1) .or. particles(1, n) .gt. prob_hi(1))
-          if (particles(1, n) .lt. prob_lo(1)) then
-             particles(1, n) = 2.d0*prob_lo(1) - particles(1, n)
-          else
-             particles(1, n) = 2.d0*prob_hi(1) - particles(1, n)
-          end if
-          vx_p(n) = -vx_p(n)
-       end do
-
-!      ... y...
-       do while (particles(2, n) .lt. prob_lo(2) .or. particles(2, n) .gt. prob_hi(2))
-          if (particles(2, n) .lt. prob_lo(2)) then
-             particles(2, n) = 2.d0*prob_lo(2) - particles(2, n)
-          else
-             particles(2, n) = 2.d0*prob_hi(2) - particles(2, n)
-          end if
-          vy_p(n) = -vy_p(n)
-       end do
-
-    end do
-
-  end subroutine warpx_push_leapfrog_positions_2d
-
 end module warpx_ES_push_particles