program hello

      use omp_lib
      implicit none
      integer :: i, n
      
      call omp_set_num_threads(2)
      
!$omp parallel private(i,n)
      n = omp_get_num_threads()
      i = omp_get_thread_num()
      print *, 'Hello World from thread ', i, ' of ', n
!$omp end parallel

      end program hello