anchor.diff:     V(1,:) = 6.0_ReKi*( Slope(2,:) - Slope(1,:) )
   
      DO I=2,NumPts-1
  +      if ( equalRealNos( U(I-1), 0.0_ReKi ) ) then
  +         CALL ExitThisRoutine ( ErrID_Fatal, RoutineName//':XAry must be monotonic.' )
  +         RETURN
  +      ENDIF
  +            
         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
  
positive.code: 
     DO I=1,NumPts-1
        DelX (I  ) =   XAry(I+1  ) - XAry(I  )
        Slope(I,:) = ( YAry(I+1,:) - YAry(I,:) )/DelX(I)
     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
  