anchor.diff:  c     # to 'qcoarse' at face 'iside'  in direction 'idir' of 'qcoarse'
   c     # Assume this is mapped.
         subroutine mb_average_face_ghost(mx,my,mbc,meqn,qcoarse,
  -     &      qfine,auxcoarse, auxfine, maux,
  +     &      qfine,auxcoarse, auxfine,
        &      idir,iface_coarse,p4est_refineFactor,refratio,igrid)
         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 auxcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,maux)
  -      double precision auxfine(1-mbc:mx+mbc,1-mbc:my+mbc,maux)
  +      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 sum, kf, kc, qf
   
  
positive.code: 
  c     Average fine grid to coarse grid or copy neighboring coarse grid
        subroutine mb_average_corner_ghost(mx,my,mbc,meqn,
       &      refratio,qcoarse,qfine,auxcoarse, auxfine, maux,
       &      icorner, is_block_bdry)
        implicit none
  
        integer mx,my,mbc,meqn,refratio,icorner, maux
        integer is_block_bdry(0:3)
        double precision auxcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,maux)
        double precision auxfine(1-mbc:mx+mbc,1-mbc:my+mbc,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)
        double precision sum
  
negative.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
        double precision qcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
        double precision qfine(1-mbc:mx+mbc,1-mbc:my+mbc,meqn)
  
  