anchor.diff:        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 auxcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,1)
  -      double precision auxfine(1-mbc:mx+mbc,1-mbc:my+mbc,1)
  +      double precision areacoarse(-mbc:mx+mbc+1,-mbc:my+mbc+1)
  +      double precision areafine(-mbc:mx+mbc+1,-mbc:my+mbc+1)
   
         double precision sum, kf, kc, qf
   
  
positive.code:                         jfine = (j-1)*refratio + jj
                          if (iface_coarse .eq. 0) then
                             qf = qfine(ifine,jfine,mq)
                             kf = auxfine(ifine,jfine,1)
                          else
                             qf = qfine(mx+1-ifine,jfine,mq)
                             kf = auxfine(mx+1-ifine,jfine,1)
                          endif
                          sum = sum + qf*kf
                       enddo
                    enddo
                    if (iface_coarse .eq. 0) then
                       kc = r2*auxcoarse(1-ibc,j+jc_add,1)
                       qcoarse(1-ibc,j+jc_add,mq) = sum/kc
                    else
                       kc = r2*auxcoarse(mx+ibc,j+jc_add,1)
                       qcoarse(mx+ibc,j+jc_add,mq) = sum/kc
                    endif
                 enddo
  
negative.code:                         ifine = (ibc-1)*refratio + ii
                          jfine = (jbc-1)*refratio + jj
                          qf = qfine(mx+1-ifine,jfine,mq)
                          kf = areafine(mx+1-ifine,jfine)
                          sum = sum + kf*qf
                       enddo
                    enddo
                    kc = r2*areacoarse(mx+ibc,1-jbc)
                    qcoarse(mx+ibc,1-jbc,mq) = sum/r2
                 enddo
              enddo
  