anchor.diff:              deltam = fmfine0(jj1,mq) + fmfine0(jj2,mq) +
        &               eff0(jj1,mq,1) + eff0(jj2,mq,1)
        
  -            deltap = fpfine1(jj1,mq) + fpfine1(jj2,mq) +
  +            deltap = (fpfine1(jj1,mq) + fpfine1(jj2,mq)) -
        &               (eff1(jj1,mq,1) + eff1(jj2,mq,1))
   
   c           # Put the same value in four ghost cells;  grab the first
  
positive.code:             do jj = jj1,jj2
                 do ii = 1,mbc
                    qfine_dummy(1-ii,jj,mq) = -deltam
                    qfine_dummy(mx+ii,jj,mq) = -deltap
                 enddo
              enddo
           enddo
  
negative.code:             ii2 = 2*(i-1) + 2
              deltam = gmfine2(ii1,mq) + gmfine2(ii2,mq) +
       &               eff2(ii1,mq,1) + eff2(ii2,mq,1)
  
              deltap = gpfine3(ii1,mq) + gpfine3(ii2,mq) -
       &               (eff3(ii1,mq,1) + eff3(ii2,mq,1))
  
              do ii = ii1,ii2
                 do jj = 1,mbc
                    qfine_dummy(ii,1-jj,mq) = -deltam
  