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:       integer ifine, jfine
        double precision kf, qf, kc
  
        integer get_block_idx, get_patch_idx
        logical debug_is_on
  
        include "debug.i"
  
        block_idx = get_block_idx()
        patch_idx = get_patch_idx()
        if (block_idx .eq. 0 .and. patch_idx .eq. 11) then
           call set_debug_on()
        endif
  
        if (debug_is_on()) then
           write(6,'(A,I2)') 'Before corner averaging (icorner = ',
       &         icorner,')'
           do i = 1-mbc,2
              do j = 1-mbc,2
                 write(6,100) i,j,qcoarse(i,j,1)
              enddo
              write(6,*) ' '
           enddo
           write(6,*) ' '
        endif
    100 format(2I5,F24.16)
  
        r2 = refratio*refratio
        do mq = 1,meqn
           if (icorner .eq. 0) then
              if (debug_is_on()) then
                 write(6,*) 'branch : icorner == 0'
              endif
              do ibc = 1,mbc
                 do jbc = 1,mbc
  c                 # Average fine grid corners onto coarse grid ghost corners
  