What subroutine should I use to sample from a Gaussian distribution?
You can use one of the following:
CALL FLNRRN (RGAUSS)
returns one number normally distributed
CALL FLNRR2 (RGAUS1, RGAUS2)
returns two such numbers, independent of each other