anchor.diff:        double precision qfine(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
         double precision qcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
   
  -      integer mq,r2
  +      integer mq
         integer i, i1, i2, ibc, ii, ifine
         integer j, j1, j2, jbc, jj, jfine
         integer ic_add, jc_add, ic, jc, mth
  
positive.code:       subroutine mb_interpolate_block_corner_ghost(mx,my,mbc,meqn,
       &      refratio, qcoarse, qfine, icorner_coarse, iblock)
        implicit none
        integer mx, my, mbc, meqn, icorner_coarse, iblock, refratio,maux
        double precision qcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
        double precision qfine(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
  
        integer mq, ibc, jbc, ii, jj, r2, ifine, jfine
        integer ic, jc, mth
        double precision gradx, grady, shiftx, shifty
        double precision sl, sr, qc, value, compute_slopes
  
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
  