anchor.diff:              end do
           end do
           !$omp end do
  +        deallocate(u1d, v1d)
           !$omp end parallel
       end subroutine spatial_winds
   
  
positive.code:         integer :: step, dpos, npos, spos, nexts, nextd, nextn
          integer :: north, south, east, west, top, bottom, n
          real :: u, v
          real :: dweight, nweight, sweight, curspd, curdir, curnsq, wind_first, wind_second
          real :: blocked
  
  
negative.code:         !$omp barrier
          !$omp do
          do k=1, nyv
  
              uk = min(k,ny)
              do i=1,nxu
                  vi = min(i,nx)
                  u1d(i)  = sum(domain%u(i,:,uk)) / nz
                  v1d(i)  = sum(domain%v(vi,:,k)) / nz
              enddo
  
  
              do j=1, nz
                  do i=1, nxu
  
  