anchor.diff:           OtherState%sigma3(i,j) =  4.0_ReKi
         end if
         
  +      !!!if ( (.not. OtherState%VRTX(i,j)) .OR. (.not. OtherState%TESF(i,j)) ) then !RRD: trying to emulate the old AD14 SEPAR.f90 with SHIFT=NOT(TESF), go back to original commented!!! out above when done
  +      !!!       OtherState%sigma3(i,j) =  2_ReKi
  +      !!!   else
  +      !!!      OtherState%sigma3(i,j) = 1.0_ReKi
  +      !!!endif
         
         if ((.not. OtherState%TESF(i,j)) .and. (Kq*dalpha0 < 0.0_ReKi)) then
            OtherState%sigma3(i,j) = 1.0_ReKi
  
positive.code:    call GetSteadyOutputs(AFInfo, alpha, Cl, Cd, Cm, Cd0, ErrStat, ErrMsg)
        if (ErrStat > ErrID_None) return
     
     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)
     if (abs(denom) < .015 ) then
        Get_f_c_from_Lookup = 1.0_ReKi
        return
  
negative.code:    
     
        ! Compute Cc_pot using eqn 1.28
     !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
  