anchor.diff:           return
         end if
         
  -
  -   
  +      
      ! Compute coefficients for cubic spline fit, clamped at two ends
      CALL AllocAry(SP_Coef,InputFileData%kp_total-1,4,4,'Spline coefficient matrix',ErrStat2,ErrMsg2)
         CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName )
         if (ErrStat >= AbortErrLev) then
            call cleanup()
            return
  -      end if
  -      
  -   SP_Coef(:,:,:) = 0.0D0
  +      end if   
  +      
  +   SP_Coef(:,:,:) = 0.0_BDKi
      temp_id = 0
      temp_id2 = 0
      DO i=1,InputFileData%member_total
  
positive.code:                hhx,hpx,TZw,Jacobian,&
                 ErrStat,ErrMsg)
  
     REAL(ReKi),         INTENT(INOUT):: GL(:)     ! GL(Gauss) point locations
     REAL(ReKi),         INTENT(IN   ):: GLL(:)     ! GLL point locations
     REAL(ReKi),         INTENT(IN   ):: uuN0(:,:) ! Initial position vector
     INTEGER(IntKi),     INTENT(IN   ):: quadrature ! Quadrature method
     INTEGER(IntKi),     INTENT(IN   ):: node_elem  ! Nodes per element
     INTEGER(IntKi),     INTENT(IN   ):: elem_total ! Total number of elements
  
negative.code: ! It also calcuates the linearized matrices Mi, Gi, and Ki for N-R algorithm
  !---------------------------------------------------------------------------
  
     REAL(BDKi),    INTENT(IN   ):: m00
     REAL(BDKi),    INTENT(IN   ):: mEta(:)
     REAL(BDKi),    INTENT(IN   ):: rho(:,:)
     REAL(BDKi),    INTENT(IN   ):: vvv(:)
     REAL(BDKi),    INTENT(IN   ):: aaa(:)
     LOGICAL,       INTENT(IN   ):: fact
     REAL(BDKi),    INTENT(  OUT):: Fi(6)
     REAL(BDKi),    INTENT(  OUT):: Mi(6,6)
     REAL(BDKi),    INTENT(  OUT):: Gi(6,6)
     REAL(BDKi),    INTENT(  OUT):: Ki(6,6)
     INTEGER(IntKi),INTENT(  OUT):: ErrStat       ! Error status of the operation
     CHARACTER(*),  INTENT(  OUT):: ErrMsg        ! Error message if ErrStat /= ErrID_None
  
     REAL(BDKi)                  :: beta(3)
     REAL(BDKi)                  :: gama(3)
     REAL(BDKi)                  :: nu(3)
     REAL(BDKi)                  :: epsi(3,3)
     REAL(BDKi)                  :: mu(3,3)
     REAL(BDKi)                  :: ome(3)
     REAL(BDKi)                  :: omd(3)
     REAL(BDKi)                  :: tempV(3)
     REAL(BDKi)                  :: tempA(3)
     INTEGER(IntKi)              :: i
  
     ErrStat = ErrID_None
     ErrMsg  = ""
     Fi(:)   = 0.0_BDKi
     Mi(:,:) = 0.0_BDKi
     Gi(:,:) = 0.0_BDKi
     Ki(:,:) = 0.0_BDKi
  
     ome(:) = vvv(4:6)
     omd(:) = aaa(4:6)
  