*********************************************************
*

      real function rangauss(xd)
*
*  return normalized gaussian random number Gauss(0,1)
*
*  last modified on 2/3/98 by:  R. Holzmann 
*********************************************************
      dimension xrand(2)
      call grndm(xrand,2)

      theta = 6.28318531*xrand(1)
      rangauss = cos(theta)*sqrt(-2.0*alog(xrand(2)))

      return
      end