anchor.diff:                          sum = sum + kf*qf
                        enddo
                     enddo
  -                  kc = r2*areacoarse(1-ibc,1-jbc)
  +                  kc = areacoarse(1-ibc,1-jbc)
                     qcoarse(1-ibc,1-jbc,mq) = sum/kc
                  enddo
               enddo
  
positive.code:                      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
  
negative.code:                         sum = sum + kf*qf
                       enddo
                    enddo
                    kc = areacoarse(mx+ibc,my+jbc)
                    qcoarse(mx+ibc,my+jbc,mq) = sum/kc
                 enddo
              enddo
           endif
        enddo
  
        end
  
  
  