anchor.diff:        implicit none
   
         integer mx,my,mbc,meqn,refratio,igrid,idir,iface_coarse
  -      integer p4est_refineFactor, maux
  +      integer p4est_refineFactor
         double precision qfine(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
         double precision qcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
         double precision areacoarse(-mbc:mx+mbc+1,-mbc:my+mbc+1)
  -      double precision   areafine(-mbc:mx+mbc+1,-mbc:my+mbc+1)
  -
  -      double precision sumfine, areaf, areac, qf
  -
  -      integer mq,r2
  +      double precision areafine(-mbc:mx+mbc+1,-mbc:my+mbc+1)
  +
  +      double precision sum, kf, kc, qf
  +
  +      integer mq
         integer i, ic_add, ibc, ii, ifine
         integer j, jc_add, jbc, jj, jfine
   
   c     # 'iface' is relative to the coarse grid
  -
  -      r2 = refratio*refratio
   
   c     # Average fine grid onto coarse grid
         do mq = 1,meqn
  
positive.code:          return
        endif
  
  
        do mq = 1,meqn
           if (iface .eq. 0) then
              do j = 1,my
  
negative.code:             do i = 1,mx
                 do jbc = 1,mbc
                    qthis(i,1-jbc,mq) = qneighbor(i,jbc,mq)
                     qneighbor(i,1-jbc,mq) = qthis(i,jbc,mq)
                 enddo
              enddo
           elseif (iface .eq. 3) then
  