anchor.diff:        xd%Cn_v_minus1(i,j)         = Cn_v
         xd%C_V_minus1(i,j)          = C_V
         OtherState%FirstPass(i,j)   = .false.
  -   
  -      if ( xd%tau_V(i,j) > 0.0 .or. LESF ) then                !! TODO:  Verify this condition 2/20/2015 GJH
  -         xd%tau_V(i,j)          = xd%tau_V(i,j) + 2.0_ReKi*p%dt*u%U / p%c(i,j)  
  + 
  +         ! If we are currently tracking a vortex, or we are in the stall region, increment tau_V
  +      if ( xd%tau_V(i,j) > 0.0 .or. LESF ) then
  +         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
  
positive.code:       VRTX = (xd%tau_V(i,j) <= 2.0_ReKi*T_VL) .and. (xd%tau_V(i,j) > 0.0_ReKi)
        
  
        
  !bjj: update sigma1 to value at t + dt      
        OtherState%sigma1(i,j) = 1.0_ReKi
        Kafactor      = Kalpha_f*dalpha0
       ! if ( fprimeprime < 0.7 .AND.  fprimeprime > 0.3 ) then 
           if ( TESF ) then  ! Separating flow
              if (Kafactor < 0.0_ReKi) then
                 OtherState%sigma1(i,j) = 2.0_ReKi  ! This must be the first check
              else if (.not. LESF ) then
                 OtherState%sigma1(i,j) = 1.0_ReKi !.4_ReKi    ! Leading edge separation has not occurred
              else if (xd%fprimeprime_minus1(i,j) <= 0.7_ReKi) then ! For this else, LESF = True
                 OtherState%sigma1(i,j) = 2.0_ReKi !1.0_ReKi 
              else
                 OtherState%sigma1(i,j) = 1.75_ReKi
              end if
           else ! Reattaching flow
          
              if (.not. LESF ) then
                 OtherState%sigma1(i,j) = 0.5_ReKi
              end if
  
              if ( VRTX .and. (xd%tau_V(i,j) <= T_VL) ) then  ! Still shedding a vortex?
                 OtherState%sigma1(i,j) = 0.25_ReKi
              end if
              if (Kafactor > 0.0_ReKi) then
                 OtherState%sigma1(i,j) = 0.75_ReKi
              end if
              
           end if
        
  !bjj: update sigma3 to value at t + dt      
           
        OtherState%sigma3(i,j) = 1.0_ReKi
        
        if ( (xd%tau_V(i,j) <= 2.0_ReKi*T_VL) .and. (xd%tau_V(i,j) >= T_VL) ) then
           OtherState%sigma3(i,j) = 3.0_ReKi
           if (.not. TESF) then
              OtherState%sigma3(i,j) = 4.0_ReKi
              if ( VRTX .and. (xd%tau_V(i,j) <= T_VL) ) then
                 if (Kafactor < 0.0_ReKi) then
                    OtherState%sigma3(i,j) = 2.0_ReKi
                 else
                    OtherState%sigma3(i,j) = 1.0_ReKi
                 end if
              end if
           end if
        else if (Kafactor < 0 ) then 
           OtherState%sigma3(i,j) = 4.0_ReKi
        end if
        
        if ((.not. TESF) .and. (Kq_f*dalpha0 < 0.0_ReKi)) then
           OtherState%sigma3(i,j) = 1.0_ReKi
        end if
     
        
       
        !---------------------------------------------------------
        ! Update the Discrete States, xd
        !---------------------------------------------------------
  
negative.code:       call SetErrStat(ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName)
        if (ErrStat>= AbortErrLev) return
              
  
     T_sh = 2.0_ReKi*(1.0_ReKi-fprimeprime) / St_sh  ! The shedding frequency for a vortices which have extended past the chord (it takes T_VL amount of time to travel along the chord)
  
        ! Are we in the stall region?  If yes, then LESF is set to true.  This means we may have conditions for tracking a new vortex at the leading edge.
     LESF = (Cn_prime > Cn1) .or. ( Cn_prime < Cn2 ) ! LE separation can occur when this is .true.; assumption is that Cn2 <  0.0 and Cn1 > 0
  
        ! Determine if any vortex we are tracking is continuing to separate or is reattaching to the airfoil
     TESF = fprimeprime < xd%fprimeprime_minus1(i,j) ! Separation point is moving towards the Leading Edge when .true.; otherwise separation point is moving toward trailing edge   
              
        ! Process VRTX-related quantities
        !!!!!!!!!!!!!!!!!!!!!
  