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:                         do jj = 1,refratio
                             ifine = (i-1)*refratio + ii
                             jfine = (j-1)*refratio + jj
                             kf = areafine(ifine,jfine)
                             volf = kf*dz
                             qf = qfine(ifine,jfine,k,m)
                             sum = sum + volf*qf
                          enddo
                      enddo
  
                      kc = areacoarse(i+ic_add,j+jc_add)
                      volc = kc*dz
                      qc = qcoarse(i+ic_add, j+jc_add,k,m)
                      cons_diff = (qc*volc - sum)/r2
  
                      do ii = 1,refratio
                          do jj = 1,refratio
                             ifine  = (i-1)*refratio + ii
                             jfine  = (j-1)*refratio + jj
                             kf = areafine(ifine,jfine)
                             volf = kf*dz
                             qfine(ifine,jfine,k,m) = qfine(ifine,jfine,k,m) + cons_diff/volf
                         end do
                      end do
                  end do !! j loop
  