anchor.diff:              DO I = 1,N  ! Loop through all nodes where the line position and tension are to be computed
   
                  IF ( ( s(I) <  0.0_DbKi ) .OR. ( s(I) >  L ) )  THEN
  -                  ErrStat = ErrID_Warn
  +                 ErrStat = ErrID_Warn
                    ErrMsg = ' All line nodes must be located between the anchor ' &
                                    //'and fairlead (inclusive) in routine Catenary().'
                    RETURN
  
positive.code: 
              ! Compute the determinant of the Jacobian matrix and the incremental
              !   tensions predicted by Newton-Raphson:
  
              DET = dXFdHF*dZFdVF - dXFdVF*dZFdHF
  
              dHF = ( -dZFdVF*EXF + dXFdVF*EZF )/DET    ! This is the incremental change in horizontal tension at the fairlead as predicted by Newton-Raphson
              dVF = (  dZFdHF*EXF - dXFdHF*EZF )/DET    ! This is the incremental change in vertical   tension at the fairlead as predicted by Newton-Raphson
  
  
negative.code:       
        ! check for NaNs - is this a good place/way to do it?
        DO J = 1, Nx
           IF (Is_NaN(REAL(x%states(J),DbKi))) THEN
              ErrStat = ErrID_Fatal
              ErrMsg = ' NaN state detected.'
           END IF
  