*     F16JSF 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 ..
      COMPLEX *16      X(1+(NMAX-1)*ABS(INCMAX))
*     .. External Subroutines ..
      EXTERNAL         BLAS_ZAMAX_VAL
*     .. Intrinsic Functions ..
      INTRINSIC        ABS
*     .. Executable Statements ..
      CONTINUE
*
      WRITE (NOUT,*) 'F16JSF/BLAS_ZAMAX_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 = ARGMAX(ABS(Re(X))+ABS(Im(X))) and
*             R = MAX(ABS(Re(X))+ABS(Im(X))).
*
         CALL BLAS_ZAMAX_VAL(N,X,INCX,K,R)
*
         WRITE (NOUT,*)
         WRITE (NOUT,99999) K
         WRITE (NOUT,99998) R
      END IF
*
99999 FORMAT (1X,'Index of absolutely largest component of X is',I3)
99998 FORMAT (1X,'Absolutely largest value is',F12.5)
      END