anchor.diff:       &                                          qcoarse,
        &                                          fpcoarse0,fmcoarse1,
        &                                          gpcoarse2,gmcoarse3,
  -     &                                          fm0, fp1,gm0, gp1,
  +     &                                          fmfine0, fpfine1,
  +     &                                          gmfine2, gpfine3,
        &                                          efc0, efc1, efc2, efc3,
        &                                          eff0, eff1, eff2, eff3,
        &                                          maskfine, qfine_dummy,
  
positive.code:       double precision gpcoarse2(mx,meqn)
        double precision gmcoarse3(mx,meqn)
  
        double precision fm0(my,meqn)
        double precision fp1(my,meqn)
        double precision gm0(mx,meqn)
        double precision gp1(mx,meqn)
  
        double precision eff0(my,meqn)
        double precision eff1(my,meqn)
        double precision eff2(mx,meqn)
        double precision eff3(mx,meqn)
  
        double precision efc0(my,meqn)
        double precision efc1(my,meqn)
        double precision efc2(mx,meqn)
        double precision efc3(mx,meqn)
  
        integer*8 transform_cptr
  
  
negative.code:                if (.not. skip_this_grid) then
                    fineval = qfine_dummy(i2(1),j2(1),mq)
                    if (iface_coarse .eq. 0) then
                       deltac = fineval + fpcoarse0(jc,mq) + efc0(jc,mq,0)
                       areac = area0(jc)
                    else
                       deltac = fineval + fmcoarse1(jc,mq) + efc1(jc,mq,0)
                       areac = area1(jc)
                    endif
                    qcoarse(ic,jc,mq) = qcoarse(ic,jc,mq) + deltac/areac
  