openacc_main.f90

program main
  use saxpy_mod
  integer, parameter :: N = 2**20
  real, dimension(N) :: X, Y

  X(:) = 1.0
  Y(:) = 0.0

  !$acc data copy(y) copyin(x)
  call saxpy(N, 2.0, x, y)
  !$acc end data

  print *, y(1)
end program