anchor.diff:        double precision efneighbor3(mx,meqn,0:1)
   
   
  -      double precision delta, deltap, deltam, areathis
  +      double precision delta, deltap, deltam, area
   
         double precision neighborval
   
         integer*8 transform_cptr
   
   
  -      integer i,j,ibc,jbc,mq
  +      integer i,j,mq
         integer i1,j1, i2, j2
   
         idir = this_iface/2
  
positive.code:          idir = 1
           do k = 0,1
              do m = 1,maux
                 auxvec(m) = aux(i,1-k,m)
              enddo
              do m = 1,meqn
                 qvec(m) = q(i,1-k,m)
  
negative.code: 
  c     # High side of 'qthis' exchanges with low side of
  c     # 'qneighbor'
        if (idir .eq. 0) then
            do mq = 1,meqn
                do j = 1,my
  c                 # Exchange at low side of 'this' grid in
  c                 # x-direction (idir == 0) 
                    j1 = j
  