program main

      include 'omp_lib.h'

      integer n
      parameter ( n = 10000000 )

      integer i
      double precision wtime, x(n), y(n), s

      wtime = omp_get_wtime ( )

      s = 123.456

c$omp parallel

c$omp do
      do i = 1, n
        x(i) = mod ( i, 17 )
        y(i) = mod ( i, 31 )
      end do
c$omp end do

c$omp do
      do i = 1, n
        y(i) = y(i) + s * x(i)
      end do
c$omp end do

c$omp end parallel

      wtime = omp_get_wtime ( ) - wtime
      write ( *, '(g14.6,a)' ) wtime, ' seconds.'

      stop
      end