anchor.diff:        if (ErrStat >= AbortErrLev) return
      
      if (UAMod == 2) then
  -       denom = ( alpha-alpha0 )*(alpha)    !NOTE: Added back (alpha) because idling cases with alpha 90º show problems with tan(alpha), the code should match steady state if the formulation in the calculation of Cc is in agreement with this formulation
  +      denom = ( alpha-alpha0 )*(alpha)    !NOTE: Added back (alpha) because idling cases with alpha 90º show problems with tan(alpha), the code should match steady state if the formulation in the calculation of Cc is in agreement with this formulation
      else
  -   denom = ( alpha-alpha0 )*tan(alpha)  !NOTE,NOTE: On 8/27/15 GJH Added back tan(alpha) because results for Fy did not match steady state without it. ! NOTE: We removed the tan(alpha) term from the equation and repace with another (alpha-alpha0) term, per Rick's suggestion 8/13/2015.    *tan(alpha)
  +      denom = ( alpha-alpha0 )*tan(alpha)  !NOTE,NOTE: On 8/27/15 GJH Added back tan(alpha) because results for Fy did not match steady state without it. ! NOTE: We removed the tan(alpha) term from the equation and repace with another (alpha-alpha0) term, per Rick's suggestion 8/13/2015.    *tan(alpha)
      endif
      
      !denom = ( alpha-alpha0 )**2  !testing again On 9/16/15
  
positive.code: 
        ! Determine if any vortex we are tracking is continuing to separate or is reattaching to the airfoil
     if (p%UAMod == 2) then   !Added specific logic for UAMod=2
        TESF = ((ABS(Cn_prime)>ABS(xd%Cn_prime_minus1(i,j))))
     else
        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   
     end if
  
negative.code: !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++  
  ! This code is taken from ADv14 but doesn't reflect the original intent of the UA theory document
  #ifdef TEST_THEORY
  IF ( p%UAMod /= 2 ) THEN
     IF ( alpha_filt_cur * Cn_prime_diff < 0. ) THEN
        T_f   = T_f0*1.5
     ELSE
        T_f   = T_f0
     ENDIF
  ENDIF
  #endif   
  !++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++  
        ! Compute alpha_f using Eqn 1.34
  