anchor.diff:        do  i=mm,1,-1
   
            best = z( j1(i1,1), j2(i1,1), j3(i1,1) )
  -c         call mpi_allreduce(best,temp,1,dp_type,
  -c     >        mpi_max,mpi_comm_world,ierr)
            call shmem_real8_max_to_all(temp,best,1,0,0,nprocs,pwrk1,
        >                               psync1)
            best = temp
  
positive.code:       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 )
  
negative.code: 
        integer i3, i2, i1,temp
  
  
        buff_id = 3 + dir
        indx = 0
  