FUNCTION get_random_U() USE omp_lib IMPLICIT NONE REAL(rp) :: get_random_U get_random_U = random_get_number_U(state) END FUNCTION get_random_U