anchor.diff:     Ax = kappa*A
      Ay = kappa*A
      
  -   K(1:12,1:12) = 0
  +   K(1:12,1:12) = 0.0_ReKi
      
      IF (Shear) THEN
  -      Kx = 12.0*E*Iyy / (G*Ax*L*L)
  -      Ky = 12.0*E*Ixx / (G*Ay*L*L)
  +      Kx = 12.0_ReKi*E*Iyy / (G*Ax*L*L)
  +      Ky = 12.0_ReKi*E*Ixx / (G*Ay*L*L)
      ELSE
  -      Kx = 0.0
  -      Ky = 0.0
  +      Kx = 0.0_ReKi
  +      Ky = 0.0_ReKi
      ENDIF
         
      K( 9,  9) = E*A/L
  -   K( 7,  7) = 12.0*E*Iyy/( L*L*L*(1.0 + Kx) )
  -   K( 8,  8) = 12.0*E*Ixx/( L*L*L*(1.0 + Ky) )
  +   K( 7,  7) = 12.0_ReKi*E*Iyy/( L*L*L*(1.0_ReKi + Kx) )
  +   K( 8,  8) = 12.0_ReKi*E*Ixx/( L*L*L*(1.0_ReKi + Ky) )
      K(12, 12) = G*Jzz/L
  -   K(10, 10) = (4.0 + Ky)*E*Ixx / ( L*(1.0+Ky) )  
  -   K(11, 11) = (4.0 + Kx)*E*Iyy / ( L*(1.0+Kx) )
  -   K( 2,  4) = -6.*E*Ixx / ( L*L*(1.0+Ky) )
  -   K( 1,  5) =  6.*E*Iyy / ( L*L*(1.0+Kx) )
  -   K( 4, 10) = (2.0-Ky)*E*Ixx / ( L*(1.0+Ky) )
  -   K( 5, 11) = (2.0-Kx)*E*Iyy / ( L*(1.0+Kx) )
  +   K(10, 10) = (4.0_ReKi + Ky)*E*Ixx / ( L*(1.0_ReKi+Ky) )  
  +   K(11, 11) = (4.0_ReKi + Kx)*E*Iyy / ( L*(1.0_ReKi+Kx) )
  +   K( 2,  4) = -6._ReKi*E*Ixx / ( L*L*(1.0_ReKi+Ky) )
  +   K( 1,  5) =  6._ReKi*E*Iyy / ( L*L*(1.0_ReKi+Kx) )
  +   K( 4, 10) = (2.0_ReKi-Ky)*E*Ixx / ( L*(1.0_ReKi+Ky) )
  +   K( 5, 11) = (2.0_ReKi-Kx)*E*Iyy / ( L*(1.0_ReKi+Kx) )
      
      K( 3,  3)  = K(9,9)
      K( 1,  1)  = K(7,7)
  
positive.code:    ! Local variables
     REAL(ReKi) :: DC(12, 12)
  
     F(1:12) = 0  ! init 
     F(3) = +T0  
     F(9) = -T0 
  
     DC = 0
     DC( 1: 3,  1: 3) = DirCos
     DC( 4: 6,  4: 6) = DirCos
     DC( 7: 9,  7: 9) = DirCos
  
negative.code: 
     t = rho*A*L;
  
     M(1:12,1:12) = 0.0_ReKi
  
     M( 1,  1) = 13._ReKi/35._ReKi * t
     M( 2,  2) = 13._ReKi/35._ReKi * t
  