anchor.diff:  
      NumPts  = SIZE( XAry )
      NumCrvs = SIZE( YAry, 2 )
  -   NumCrvs = SIZE( YAry, 2 )
  +   !NumCrvs = SIZE( YAry, 2 )
   
      ALLOCATE ( Res( NumCrvs ) , STAT=ErrStatLcl )
      IF ( ErrStatLcl /= 0 )  THEN
  
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:    ZHi(:) = 0.0_ReKi
  
     DO I=NumPts-1,1,-1
        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) )                     ! bjj: already checked for DelX(I) == 0
        ZHi(:)      = ZLo(:)
     END DO ! I
  
  