anchor.diff:  
      DO I=1,NumPts-1
         DelX (I  ) =   XAry(I+1  ) - XAry(I  )
  +      
  +      if ( equalRealNos( DelX(I), 0.0_ReKi ) ) then
  +         CALL ExitThisRoutine ( ErrID_Fatal, RoutineName//':XAry must have unique values.' )
  +         RETURN
  +      ENDIF
  +      
         Slope(I,:) = ( YAry(I+1,:) - YAry(I,:) )/DelX(I)
      END DO ! I
   
  
positive.code:    V(1,:) = 6.0_ReKi*( Slope(2,:) - Slope(1,:) )
  
     DO I=2,NumPts-1
        U(I)   = 2.0_ReKi*( DelX (I-1) + DelX (I)     ) - DelX(I-1)*DelX(I-1  )/U(I-1)
        V(I,:) = 6.0_ReKi*( Slope(I,:) - Slope(I-1,:) ) - DelX(I-1)*   V(I-1,:)/U(I-1)
     END DO ! I
  
negative.code: 
     NumPts  = SIZE( XAry )
     NumCrvs = SIZE( YAry, 2 )
     !NumCrvs = SIZE( YAry, 2 )
  
     ALLOCATE ( Res( NumCrvs ) , STAT=ErrStatLcl )
     IF ( ErrStatLcl /= 0 )  THEN
  