在Fortran中，有没有一种用随机数生成器进行并发计算的简单方法？

program random_concurrent
    use iso_fortran_env, only: real64
    implicit none

    integer :: i, nthreads
    real(real64) :: random_number

    ! 设置线程数
    nthreads = omp_get_max_threads()
    print *, "Number of threads: ", nthreads

    !$omp parallel do private(i, random_number)
    do i = 1, 100
        call random_number(random_number)
        print *, "Thread ", omp_get_thread_num(), " - Random number: ", random_number
    end do
    !$omp end parallel do

end program random_concurrent