anchor.diff:        include 'mpinpb.h'
         include 'globals.h'
         include 'mpp/shmem.fh'
  -c      include 'new_global.h'
   
         integer axis, dir, n1, n2, n3
         double precision u( n1, n2, n3 )
  
positive.code: c---------------------------------------------------------------------
  c     exchange boundary points
  c---------------------------------------------------------------------
  c      print *,"comm3 called from psinv"
        call comm3(u,n1,n2,n3,k)
  
        if( debug_vec(0) .ge. 1 )then
  
negative.code: 
              call shmem_double_put(buff(1,buff_id1),buff(1,buff_id),
       >           buff_len, nbr( axis, dir, k ))
  
  
           endif
  