anchor.diff:     
      
         ! Compute Cc_pot using eqn 1.28
  -   Cc_pot          = Cn_alpha_q_circ*tan(alpha_e+alpha0)
  -   
  +   !Cc_pot          = Cn_alpha_q_circ*tan(alpha_e+alpha0)
  +   Cc_pot          = Cn_alpha_q_circ*alpha_e !testing 9/16/15 without tangent
      
      if (OtherState%FirstPass(i,j)) then
         Cn_pot_minus1 = Cn_pot
  
positive.code:          OtherState%sigma3(i,j) =  4.0_ReKi
        end if
        
        
        if ((.not. OtherState%TESF(i,j)) .and. (Kq*dalpha0 < 0.0_ReKi)) then
           OtherState%sigma3(i,j) = 1.0_ReKi
  
negative.code:          if  ( Kafactor < 0.0_ReKi ) then ! .AND.  .NOT. OtherState%VRTX(OtherState%iBladeNode, OtherState%iBlade )  )then
              y%Cc = eta_e*Cc_pot*(sqrt(fprimeprime_c) - f_c_offset)
           else         
              !y%Cc = eta_e*Cc_pot*(sqrt(fprimeprime_c) - f_c_offset) + Cn_v*tan(alpha_e)*(1-xd%tau_v(OtherState%iBladeNode, OtherState%iBlade)/(2.0*T_VL))
              y%Cc = eta_e*Cc_pot*(sqrt(fprimeprime_c) - f_c_offset) + Cn_v*alpha_e*(1-xd%tau_v(OtherState%iBladeNode, OtherState%iBlade)/(2.0*T_VL)) !testing without tan 9/16/15
           end if
           
        elseif ( p%UAMod == 2 ) then
  