FUNCTION get_random_N() USE omp_lib IMPLICIT NONE REAL(rp) :: get_random_N get_random_N = random_get_number_N(state) END FUNCTION get_random_N