anchor.diff:        end do
         end do
   
  -      flop = flop + m*15.0
  -      ! flop = flop + rix*15.0
  +      flop = flop + m*12.0
         
       case (Y_minus)
         m = 0.0
  
positive.code:           Vp0 = v0(i  ,j  ,k  ,2)
            Wp0 = v0(i  ,j  ,k  ,3)
            
            Uw1  = 2.0*u_bc_ref - Up0
            Vw1  = 2.0*v_bc_ref - Vp0
            Ww1  = 2.0*w_bc_ref - Wp0
            
            EX = Uw1 - Up0
            EY = Vw1 - Vp0
  
negative.code:       end do
        end do
        
        flop = flop + m*12.0
        
      case (X_plus)
        m = 0.0
  