New report and code for final submission.HEADmaster

1 files changed, 0 insertions, 33 deletions

diff --git a/class/putParticleInBox.f90 b/class/putParticleInBox.f90
deleted file mode 100644
index d98a116..0000000
--- a/class/putParticleInBox.f90
+++ /dev/null
@@ -1,33 +0,0 @@
-   ! intakes coordinates and outputs index
-   integer function indexfxn(px, py) result(idx)
-     double precision, intent(in) :: px, py
-     double precision :: boxSize,yCoord
-      integer :: ix, iy
-
-     !> Check if particle is in box
-     if(abs(px).GT.L/2 .OR. abs(py).GT.L/2) then
-         !> Particle is outside, don't track it
-         idx=b*b ! This is the "recycle bin" for untracked particles.
-     else
-      ! the sector size      
-      ! M = floor(L / rc)
-      yCoord = - py
-      ! coordinate + half length of boundary [-L/2,L/2] -> [0,L] 
-      ! divided by effective particle size to give coordinate
-      boxSize=L/real(b,8)
-      ix = floor((px + L/2.0d0)/boxSize)
-      iy = floor((yCoord + L/2.0d0)/boxSize)
-
-      ! for particles on the boundary
-      ix = min(max(ix, 0), b-1)
-      iy = min(max(iy, 0), b-1)
-      if (ix < 0 .or. ix >= b .or. iy < 0 .or. iy >= b) then
-         print *, "ERROR: particle outside sector grid"
-         print *, "px, py =", px, py
-         stop
-      end if
-
-      ! call rey's function
-      idx = TWOtoONE(ix,iy)
-     end if
-   end function indexfxn