anchor.diff:      double precision :: qcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,1-mbc:mz+mbc,meqn)
       double precision ::   qfine(1-mbc:mx+mbc,1-mbc:my+mbc,1-mbc:mz+mbc,meqn)
   
  -    double precision :: areacoarse(-mbc:mx+mbc+1,-mbc:my+mbc+1,-mbc:mz+mbc+1)
  -    double precision ::   areafine(-mbc:mx+mbc+1,-mbc:my+mbc+1,-mbc:mz+mbc+1)
  +    double precision :: areacoarse(-mbc:mx+mbc+1,-mbc:my+mbc+1)
  +    double precision ::   areafine(-mbc:mx+mbc+1,-mbc:my+mbc+1)
   
       integer :: ii, jj, i,j, i1, i2, j1, j2, ig, jg, mq, mth, k
       integer :: ic,jc,ic_add, jc_add, ifine, jfine
  
positive.code:     end do mq_loop
  
      if (manifold .ne. 0) then
          call fclaw2d_clawpatch46_fort3_fixcapaq2(mx,my,mz,mbc,meqn, & 
                          qcoarse,qfine, areacoarse,areafine,igrid)
      endif
  
negative.code: 
      double precision ::  qcoarse(1-mbc:mx+mbc,1-mbc:my+mbc,1-mbc:mz+mbc,meqn)
      double precision ::    qfine(1-mbc:mx+mbc,1-mbc:my+mbc,1-mbc:mz+mbc,meqn)
      double precision ::  areacoarse(-mbc:mx+mbc+1,-mbc:my+mbc+1)
      double precision ::    areafine(-mbc:mx+mbc+1,-mbc:my+mbc+1)
  
      integer :: i,j,k,ii, jj, ifine, jfine, m, ig, jg, ic_add, jc_add
      double precision :: kf, kc, r2, sum, cons_diff, qf, qc, volf, dz, volc
  
      p4est_refineFactor = 2
      refratio = 2
  
      !! This is a bogus value, since we shouldn't end up here (yet). 
      dz = 1
  
      !! # Get (ig,jg) for grid from linear (igrid) coordinates
      ig = mod(igrid,refratio)
  