anchor.diff:                          do m = 0,r2-1
                              sum = sum + qfine(i2(m),j2(m),mq)
                           enddo
  -                        qcoarse(ic,jc,mq) = sum/r2
  +                        qcoarse(ic,jc,mq) = sum/dble(r2)
                        endif
                     endif
                  enddo
  
positive.code:                         do m = 0,r2-1
                             sum = sum + qfine(i2(m),j2(m),mq)
                          enddo
                          qcoarse(ic,jc,mq) = sum/r2
                       endif
                    endif
                 enddo
  
negative.code:                   do m = 0,r2-1
                       sum = sum + qfine(i2(m),j2(m),mq)
                    enddo
                    qcoarse(i1,j1,mq) = sum/dble(r2)
                 endif
  
  
  