Random Gaussian Number in Fortran

I am looking for a routine which returns a sequence of normally distributed pseudorandom numbers in Fortran 90. I've already heard about the routine NORMAL (https://people.sc.fsu.edu/~jburkardt/f_src/normal/normal.html), but how can I call it? And how can I set my Input?

Or do you know another routine?

Thanks for your time