ran1.F

Go to the documentation of this file.

!
! $Id: ran1.F 1403 2010-07-01 09:02:53Z fairhead $
!
      FUNCTION ran1(IDUM)
      dimension r(97)
      save r
      save iff,ix1,ix2,ix3
      parameter(m1=259200,ia1=7141,ic1=54773,rm1=3.8580247e-6)
      parameter(m2=134456,ia2=8121,ic2=28411,rm2=7.4373773e-6)
      parameter(m3=243000,ia3=4561,ic3=51349)
      DATA iff /0/
      IF (idum.LT.0.OR.iff.EQ.0) THEN
        iff=1
        ix1=mod(ic1-idum,m1)
        ix1=mod(ia1*ix1+ic1,m1)
        ix2=mod(ix1,m2)
        ix1=mod(ia1*ix1+ic1,m1)
        ix3=mod(ix1,m3)
        DO 11 j=1,97
          ix1=mod(ia1*ix1+ic1,m1)
          ix2=mod(ia2*ix2+ic2,m2)
          r(j)=(REAL(ix1)+REAL(ix2)*rm2)*rm1
11      CONTINUE
        idum=1
      ENDIF
      ix1=mod(ia1*ix1+ic1,m1)
      ix2=mod(ia2*ix2+ic2,m2)
      ix3=mod(ia3*ix3+ic3,m3)
      j=1+(97*ix3)/m3
      IF(j.GT.97.OR.j.LT.1)pause
      ran1=r(j)
      r(j)=(REAL(ix1)+REAL(ix2)*rm2)*rm1
      RETURN
      END