anchor.diff:                          jfine = (jbc-1)*refratio + jj
                           if (iface_coarse .eq. 2) then
                              qf = qfine(ifine,jfine,mq)
  -                           kf = auxfine(ifine,jfine,1)
  +                           kf = areafine(ifine,jfine)
                           else
                              qf = qfine(ifine,my+1-jfine,mq)
  -                           kf = auxfine(ifine,my+1-jfine,1)
  +                           kf = areafine(ifine,my+1-jfine)
                           endif
                           sum = sum + kf*qf
                        enddo
                     enddo
                     if (iface_coarse .eq. 2) then
  -                     kc = r2*auxcoarse(i+ic_add,1-jbc,1)
  +                     kc = r2*areacoarse(i+ic_add,1-jbc)
                        qcoarse(i+ic_add,1-jbc,mq) = sum/kc
                     else
  -                     kc = r2*auxcoarse(i+ic_add,my+jbc,1)
  +                     kc = r2*areacoarse(i+ic_add,my+jbc)
                        qcoarse(i+ic_add,my+jbc,mq) = sum/kc
                     endif
                  enddo
  
positive.code:                         ifine = (ibc-1)*refratio + ii
                          jfine = (jbc-1)*refratio + jj
                          qf = qfine(ifine,jfine,mq)
                          kf = auxfine(ifine,jfine,1)
                          sum = sum + kf*qf
                       enddo
                    enddo
                    kc = r2*auxcoarse(1-ibc,1-ibc,1)
                    qcoarse(1-ibc,1-jbc,mq) = sum/r2
                 enddo
              enddo
  
negative.code:                         jfine = (j-1)*refratio + jj
                          if (iface_coarse .eq. 0) then
                             qf = qfine(ifine,jfine,mq)
                             kf = areafine(ifine,jfine)
                          else
                             qf = qfine(mx+1-ifine,jfine,mq)
                             kf = areafine(mx+1-ifine,jfine)
                          endif
                          sum = sum + qf*kf
                       enddo
                    enddo
                    if (iface_coarse .eq. 0) then
                       kc = r2*areacoarse(1-ibc,j+jc_add)
                       qcoarse(1-ibc,j+jc_add,mq) = sum/kc
                    else
                       kc = r2*areacoarse(mx+ibc,j+jc_add)
                       qcoarse(mx+ibc,j+jc_add,mq) = sum/kc
                    endif
                 enddo
  