anchor.diff:            Vp0 = v0(i  ,j  ,k  ,2)
             Wp0 = v0(i  ,j  ,k  ,3)
   
  -          Ub1  = 2.0*u_bc_ref - Up0
  -          Vb1  = 2.0*v_bc_ref - Vp0
  -          Wb1  = 2.0*w_bc_ref - Wp0
  +          Ub1  = u_bc_ref2 - Up0
  +          Vb1  = v_bc_ref2 - Vp0
  +          Wb1  = w_bc_ref2 - Wp0
             
             EX = Ub1 - Up0
             EY = Vb1 - Vp0
  
positive.code:           Vp0 = v0(i  ,j  ,k  ,2)
            Wp0 = v0(i  ,j  ,k  ,3)
  
            Ut1  = 2.0*u_bc_ref - Up0
            Vt1  = 2.0*v_bc_ref - Vp0
            Wt1  = 2.0*w_bc_ref - Wp0
            
            EX = Ut1 - Up0
            EY = Vt1 - Vp0
  
negative.code:       end do
        end do
  
        flop = flop + m*12.0
        
      case default
      end select FACES
  