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