*     F16JRF Example Program Text
*     Mark 22 Release. NAG Copyright 2007.
*     .. Parameters ..
      INTEGER          NIN, NOUT
      PARAMETER        (NIN=5,NOUT=6)
      INTEGER          NMAX, INCMAX
      PARAMETER        (NMAX=10,INCMAX=NMAX)
*     .. Local Scalars ..
      DOUBLE PRECISION R
      INTEGER          I, INCX, K, N
*     .. Local Arrays ..
      DOUBLE PRECISION X(1+(NMAX-1)*ABS(INCMAX))
*     .. External Subroutines ..
      EXTERNAL         BLAS_DAMIN_VAL
*     .. Intrinsic Functions ..
      INTRINSIC        ABS
*     .. Executable Statements ..
      CONTINUE
*
      WRITE (NOUT,*) 'F16JRF/BLAS_DAMIN_VAL Example Program Results'
*
*     Skip heading in data file
*
      READ (NIN,*)
*
*     Read N and INCX from data file
*
      READ (NIN,*) N, INCX
*
      IF (N.LE.NMAX .AND. ABS(INCX).LE.INCMAX) THEN
*
*        Read X from data file
*
         READ (NIN,*) (X(I),I=1,1+(N-1)*ABS(INCX),INCX)
*
*        Find K = ARGMIN(ABS(X)) and R = MIN(ABS(X)).
*
         CALL BLAS_DAMIN_VAL(N,X,INCX,K,R)
*
         WRITE (NOUT,*)
         WRITE (NOUT,99999) K
         WRITE (NOUT,99998) R
      END IF
*
99999 FORMAT (1X,'Index of absolutely smallest component of X is',I3)
99998 FORMAT (1X,'Absolutely smallest value is',F12.5)
      END