anchor.diff:        !! NEW CODE 2/19/2015
         VRTX = (xd%tau_V(i,j) <= 2.0_ReKi*T_VL) .and. (xd%tau_V(i,j) > 0.0_ReKi)
         
  -      
  -      if (( xd%tau_V(i,j) >= (T_VL + T_sh) ) .and. TESF)  then !.and. (TESF) RRD added 
  -         xd%tau_V(i,j)     = 0.0_ReKi
  -        ! LESF=.FALSE. !also added this
  -      end if
  +
         
   !bjj: update sigma1 to value at t + dt      
         OtherState%sigma1(i,j) = 1.0_ReKi
  
positive.code:             end if
              
           end if
        !!!if (.not. LESF ) then  !RRD: trying to emulate the old AD14 SEPAR.f90 with SHIFT=NOT(LESF), go back to original commented!!! out above when done
        !!!        OtherState%sigma1(i,j) = 0.667_ReKi
        !!!end if
        
  !bjj: update sigma3 to value at t + dt      
           
  
negative.code:          xd%tau_V(i,j)          = xd%tau_V(i,j) + 2.0_ReKi*p%dt*u%U / p%c(i,j)  
        end if
     
        
           ! If we a have been tracking a vortex and 1) it is now past the chord [T_VL] and 2) we have gone beyond the next shedding period [T_sh], and 
           !   3) we are continuing the flow serparation, we will shed the existing vortex so that we can create a new one at the leading edge
        if (( xd%tau_V(i,j) >= (T_VL + T_sh) ) .and. TESF)  then !.and. (TESF) RRD added
           xd%tau_V(i,j)          = 0.0_ReKi
        end if
        
  #ifdef UA_OUTS
     m%TESF(i,j) = TESF  
  