1 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   2 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   3 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   4 | c$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
   5 | c.......................................................................
   6 | c
   7 | c     tirage de nombres aleatoires en double precision
   8 | c
   9 | c.......................................................................
  10 |       subroutine unra(xunra)
  11 |       implicit double precision (a-h,o-z)
  12 |       real urand
  13 | c.......................................................................
  14 |       common/comsee/iseed
  15 | c.......................................................................
  16 |       if (iseed.le.0) iseed=14396
  17 | 1000  continue
  18 |       xunra=dble(urand(iseed))
  19 |       if ((xunra.eq.0.d0).or.(xunra.eq.1.d0)) goto 1000
  20 | c.......................................................................
  21 |       return
  22 |       end