program main

  use omp_lib

  integer, parameter :: n = 10000000

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

  wtime = omp_get_wtime ( )

!$omp parallel

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

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

!$omp end parallel

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