I don’t know how RANDOM_NUMBER() actually works behind the scenes, but if someone has some insight of that I’ll be interested since I do stochastic simulations as well.
However, if having it inside the parallel is a problem, what about “precomputing” it before entering the parallel region and storing it into an array?