anchor.diff:  
       integer(kind=is) :: etag, wtag, ntag, stag
       integer(kind=is) :: ierr,status1 
  -    integer(kind=is) :: ilev
  +    !!integer(kind=is) :: ilev
   
       real(kind=rl),dimension(:,:,:),allocatable :: sendN,recvN,sendS,recvS
       real(kind=rl),dimension(:,:,:),allocatable :: sendE,recvE,sendW,recvW
  
positive.code:     endif
  
      if (southeast.ne.MPI_PROC_NULL) then
         p(:,1-nh:0,nx+1:nx+nh) = recvE
      else !!Homogenous Neumann  
         p(:,1-nh:0,nx+1:nx+nh) = p(:,nh:1:-1,nx:nx-nh:-1)
      end if
      !
      deallocate(sendSE)
  
negative.code:     if (east.ne.MPI_PROC_NULL) then
         p(:,1:ny,nx+1:nx+nh) = recvE
      else !!Homogenous Neumann  
         p(:,1:ny,nx+1:nx+nh) = p(:,1:ny,nx:nx-nh+1:-1)
      end if
      !
      deallocate(sendE)
  