C Illustrates passing arrays as arguments.

C Created:  Joe Zachary, November 1, 1992
C Modified:


C This subroutine computes the dot product of two REAL arrays
C (X and Y) of length 5.  The answer is communicated back
C via Z.

      SUBROUTINE DOTPROD (X, Y, Z)

      IMPLICIT NONE
      INTEGER I, SIZE
      PARAMETER (SIZE=5)
      REAL X(SIZE), Y(SIZE), Z(SIZE)

      DO I = 1, SIZE
         Z(I) = X(I) * Y(I)
      ENDDO

      RETURN
      END


C The main program exercises DOTPROC

      PROGRAM TESTER

      IMPLICIT NONE
      INTEGER I, SIZE
      PARAMETER (SIZE=5)
      REAL A(SIZE), B(SIZE), C(SIZE)

      DO I = 1, SIZE
         A(I) = I
         B(I) = SIZE
      ENDDO

      PRINT *, 'A   = ', A
      PRINT *, 'B   = ', B

      CALL DOTPROD(A, B, C)

      PRINT *, 'A*B = ', C

      STOP
      END