anchor.diff:              qfine(iff,jff,mq) = value
            enddo
   
  +      enddo
  +      enddo 
         enddo
   
         end
  
positive.code:          if (idir .eq. 0) then
  c           # this ensures that we get 'hanging' corners
  
              if (iface_coarse .eq. 0) then
                 ic = 1
              elseif (iface_coarse .eq. 1) then
                 ic = mx
              endif
              do jc = 1,mx
                 i1 = ic
  
negative.code:                   enddo
                 endif
              enddo
              enddo
           else
              do jbc = 1,mbc/2
              if (iface_coarse .eq. 2) then
                 jc = jbc
              elseif (iface_coarse .eq. 3) then
                 jc = my - jbc + 1
              endif
              do ic = 1,mx
      1          i1 = ic
  