anchor.diff:    do while (time .lt. tfin)
       iter = iter + 1
   
  -    ! Simulation the fluid flow
  -    !$OMP PARALLEL
  +    ! Simulate the fluid flow
       call step
       time = time + dt
   
  
positive.code:   integer,dimension(:,:),allocatable :: nbs
    integer,dimension(:,:,:),allocatable :: part_count
    integer,dimension(:,:,:,:),allocatable :: binpart
  
    ! Expected max particles in single bin
    integer :: nbinx, nbiny, nbinz
  
negative.code:     gwmag = sxGW*sxGW + syGW*syGW + szGW+szGW
      gwmag = sqrt(gwmag)
  
      fx(i) = fx(i) - STC*mass*LW*sxGW/gwmag
      fy(i) = fy(i) - STC*mass*LW*syGW/gwmag
      fz(i) = fz(i) - STC*mass*LW*szGW/gwmag
  
    end do
  
  