anchor.diff:        
         if ( p%flookup ) then 
         !if ( p%UAMod == UA_Gonzalez ) then
  -         Cc_FS = BL_p%eta_e*KC%Cc_pot *(sqrt(KC%fprimeprime_c) - Gonzales_factor)                                                                  ! Eqn 1.40  
  +         Cc_FS = BL_p%eta_e*KC%Cc_pot *(sqrt(KC%fprimeprime_c) - Gonzalez_factor)                                                                  ! Eqn 1.40  
         else
            Cc_FS = BL_p%eta_e*KC%Cc_pot * sqrt(KC%fprimeprime_c)                                                                                     ! Eqn 1.40 without Gonzalez's modification for negative Cc
         end if
  
positive.code:       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 )
                    
     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
        ! 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
     
  
negative.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 :: Gonzalez_factor = 0.2_ReKi   ! this factor, proposed by Gonzalez (for "all" models) is used to modify Cc to account for negative values seen at f=0 (see Eqn 1.40)
     
     contains
     
  