anchor.diff:     ZHi(:) = 0.0_ReKi
   
      DO I=NumPts-1,1,-1
  -      ZLo(:)      = ( V(I,:) - DelX(I)*ZHi(:) )/U(I)
  +      ZLo(:)      = ( V(I,:) - DelX(I)*ZHi(:) )/U(I)                             ! bjj: already checked for u(I) == 0
         Coef(I,:,1) = Slope(I,:) - DelX(I)*( ZHi(:)/6.0_ReKi + ZLo(:)/3.0_ReKi )
         Coef(I,:,2) = 0.5_ReKi*ZLo(:)
  -      Coef(I,:,3) = ( ZHi(:) - ZLo(:) )/( 6.0_ReKi*DelX(I) )
  +      Coef(I,:,3) = ( ZHi(:) - ZLo(:) )/( 6.0_ReKi*DelX(I) )                     ! bjj: already checked for DelX(I) == 0
         ZHi(:)      = ZLo(:)
      END DO ! I
   
  
positive.code: 
     NumPts  = SIZE( XAry )
     NumCrvs = SIZE( YAry, 2 )
     NumCrvs = SIZE( YAry, 2 )
  
     ALLOCATE ( Res( NumCrvs ) , STAT=ErrStatLcl )
     IF ( ErrStatLcl /= 0 )  THEN
  
negative.code: 
     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
  
  