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