anchor.diff:        else
            denom = eta_e*c_nalpha_circ*( alpha-alpha0 )*tan(alpha)
         endif
  -      Get_f_c_from_Lookup =  min(fc_limit, (  Cc / denom  + Gonzales_factor ) **2 )
  +      Get_f_c_from_Lookup =  min(fc_limit, (  Cc / denom  + Gonzalez_factor ) **2 )
                     
      end if
  -      ! Apply an offset of 0.2 to fix cases where f_c should be negative, but we are using **2 so can only return positive values
  +      ! Apply an offset of Gonzalez_factor = 0.2 to fix cases where f_c should be negative, but we are using **2 so can only return positive values
         ! Note: because we include this offset, it must be accounted for in the final value of Cc, eqn 1.40.  This will be applied
         ! For both UA_Mod = 1,2, and 3 when using Flookup = T
      
  
positive.code:    integer(intki), parameter :: UA_Gonzalez      = 2   ! UAMod = 2 [Gonzalez's variant (changes in Cn,Cc,Cm)]
     integer(intki), parameter :: UA_MinemmaPierce = 3   ! UAMod = 3 [Minemma/Pierce variant (changes in Cc and Cm)]
     
     real(ReKi),     parameter :: Gonzales_factor = 0.2_ReKi   ! this factor, proposed by Gonzales (for "all" models) is used to modify Cc to account for negative values seen at f=0 (see Eqn 1.40)
     
     contains
     
  
negative.code:    character(*),     intent(  out) :: ErrMsg        ! Error message if ErrStat /= ErrID_None
     
     
     real(ReKi), parameter           :: fc_limit = (1.0_ReKi + Gonzalez_factor)**2    ! normally, fc is limited by 1, but we're limiting (sqrt(fc)-Gonzalez_factor) to 1, so fc is limited to 1.44 instead (when Gonzalez_factor is 0.2)
     real(ReKi)                      :: Cc, Cl, Cd, Cm, Cd0, Cm0, denom
     ErrStat = ErrID_None
     ErrMsg  = ''
        ! NOTE:  This subroutine call cannot live in Blade Element because BE module calls UnsteadyAero module.
     
    
        ! in cases where denom is zero, Get_f_c_from_Lookup = min(fc_limit, inf)
     if (EqualRealNos(real(alpha,SiKi), 0.0_SiKi)) then
        
        Get_f_c_from_Lookup = fc_limit
        
     elseif (EqualRealNos(real(alpha,SiKi), real(alpha0,SiKi))) then
        
        Get_f_c_from_Lookup = fc_limit
        
     else if (EqualRealNos( real(c_nalpha_circ,SiKi), 0.0_SiKi )) then
        
        Get_f_c_from_Lookup = fc_limit
  
     else
           
  