program main

  use omp_lib

  integer, parameter :: n = 10000000

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

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

!$omp parallel

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

!$omp end parallel

  stop
end