anchor.diff:      real                                                        ::  dh, dh1, dh2, flop, vcs, EX, EY, EZ, rei
       real                                                        ::  u_ref, v_ref, w_ref, rix, rjx, rkx, m
       real                                                        ::  u_bc, v_bc, w_bc, u_bc_ref, v_bc_ref, w_bc_ref
  +    real                                                        ::  u_bc_ref2, v_bc_ref2, w_bc_ref2
       real, dimension(1-g:sz(1)+g, 1-g:sz(2)+g, 1-g:sz(3)+g, 3)   ::  v0, wv
       real, dimension(0:3)                                        ::  v00
       real, dimension(3)                                          ::  vec
  
positive.code:       end do
        end do
  
        flop = flop + m*15.0
        ! flop = flop + rix*15.0
        
      case (Y_minus)
        m = 0.0
  
negative.code:           Vp0 = v0(i  ,j  ,k  ,2)
            Wp0 = v0(i  ,j  ,k  ,3)
  
            Us1  = u_bc_ref2 - Up0
            Vs1  = v_bc_ref2 - Vp0
            Ws1  = w_bc_ref2 - Wp0
          
            EX = Us1 - Up0
            EY = Vs1 - Vp0
  