anchor.diff:                          sum = sum + kf*qf
                        enddo
                     enddo
  -                  kc = r2*areacoarse(mx+ibc,1-jbc)
  +                  kc = areacoarse(mx+ibc,1-jbc)
                     qcoarse(mx+ibc,1-jbc,mq) = sum/kc
                  enddo
               enddo
  
positive.code:       integer mq, ibc, jbc, ii, jj, r2, ifine, jfine
        double precision sum,qf,kf,kc
  
        r2 = refratio**2
        do mq = 1,meqn
           if (icorner .eq. 0) then
  
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
  
  
  