anchor.diff:      end if
   
     end do
  +  !$OMP END DO
   
     return
   end subroutine enforce_BCs
  
positive.code:   use state
    use neighbors
    use posvel
    implicit none
    integer :: i, j, q, p, nb, l
    integer :: binx, biny, binz
  
negative.code:   implicit none
    integer :: i
  
    !$OMP DO
    do i = 1, n
  
      if(px(i) .lt. 0.0_WP) then
  