SUBROUTINE initialize_random_N(seed) USE omp_lib IMPLICIT NONE INTEGER, INTENT(IN) :: seed state = random_construct_N(seed) END SUBROUTINE initialize_random_N