anchor.diff:      call verify_array_shape_3D(grid(lev)%recvW,nz,ny,nh)
       recvW => grid(lev)%recvW
   
  -    call verify_array_shape_3D(grid(lev)%sendSW,nz,nh,nh)
  -    sendSW => grid(lev)%sendSW
  -    call verify_array_shape_3D(grid(lev)%sendSE,nz,nh,nh)
  -    sendSE => grid(lev)%sendSE
  -    call verify_array_shape_3D(grid(lev)%sendNW,nz,nh,nh)
  -    sendNW => grid(lev)%sendNW
  -    call verify_array_shape_3D(grid(lev)%sendNE,nz,nh,nh)
  -    sendNE => grid(lev)%sendNE
  -
  -    call verify_array_shape_3D(grid(lev)%recvSW,nz,nh,nh)
  -    recvSW => grid(lev)%recvSW
  -    call verify_array_shape_3D(grid(lev)%recvSE,nz,nh,nh)
  -    recvSE => grid(lev)%recvSE
  -    call verify_array_shape_3D(grid(lev)%recvNW,nz,nh,nh)
  -    recvNW => grid(lev)%recvNW
  -    call verify_array_shape_3D(grid(lev)%recvNE,nz,nh,nh)
  -    recvNE => grid(lev)%recvNE
  +!!$    call verify_array_shape_3D(grid(lev)%sendSW,nz,nh,nh)
  +!!$    sendSW => grid(lev)%sendSW
  +!!$    call verify_array_shape_3D(grid(lev)%sendSE,nz,nh,nh)
  +!!$    sendSE => grid(lev)%sendSE
  +!!$    call verify_array_shape_3D(grid(lev)%sendNW,nz,nh,nh)
  +!!$    sendNW => grid(lev)%sendNW
  +!!$    call verify_array_shape_3D(grid(lev)%sendNE,nz,nh,nh)
  +!!$    sendNE => grid(lev)%sendNE
  +!!$
  +!!$    call verify_array_shape_3D(grid(lev)%recvSW,nz,nh,nh)
  +!!$    recvSW => grid(lev)%recvSW
  +!!$    call verify_array_shape_3D(grid(lev)%recvSE,nz,nh,nh)
  +!!$    recvSE => grid(lev)%recvSE
  +!!$    call verify_array_shape_3D(grid(lev)%recvNW,nz,nh,nh)
  +!!$    recvNW => grid(lev)%recvNW
  +!!$    call verify_array_shape_3D(grid(lev)%recvNE,nz,nh,nh)
  +!!$    recvNE => grid(lev)%recvNE
   
       comm(:) = 0
       req(:)  = MPI_REQUEST_NULL
  
positive.code:     integer(kind=ip) :: nx, ny, nz
      integer(kind=ip) :: nh
      integer(kind=ip) :: south, east, north, west
      integer(kind=ip) :: southwest, southeast, northeast, northwest
  
      integer(kind=ip) :: sntag, ewtag, nstag, wetag
      integer(kind=ip) :: swnetag, senwtag, nwsetag, neswtag
  
      integer(kind=ip) :: i, j
      integer(kind=ip) :: icount
  
negative.code:     east      = grid(lev)%neighb(2)
      north     = grid(lev)%neighb(3)
      west      = grid(lev)%neighb(4)
  !!$    southwest = grid(lev)%neighb(5)
  !!$    southeast = grid(lev)%neighb(6)
  !!$    northeast = grid(lev)%neighb(7)
  !!$    northwest = grid(lev)%neighb(8)
  
      call verify_array_shape_3D(grid(lev)%sendS,nz,nh,nx)
      sendS => grid(lev)%sendS
  