3
parameter (m1=259200,ia1=7141,ic1=54773,rm1=3.8580247e-6)
4
parameter (m2=134456,ia2=8121,ic2=28411,rm2=7.4373773e-6)
5
parameter (m3=243000,ia3=4561,ic3=51349)
8
if (idum.lt.0.or.iff.eq.0) then
11
ix1=mod(ia1*ix1+ic1,m1)
13
ix1=mod(ia1*ix1+ic1,m1)
16
ix1=mod(ia1*ix1+ic1,m1)
17
ix2=mod(ia2*ix2+ic2,m2)
18
r(j)=(float(ix1)+float(ix2)*rm2)*rm1
22
ix1=mod(ia1*ix1+ic1,m1)
23
ix2=mod(ia2*ix2+ic2,m2)
24
ix3=mod(ia3*ix3+ic3,m3)
26
if(j.gt.97.or.j.lt.1)then
27
write(*,*) 'j is bad in ran1.f',j, 97d0*ix3/m3
31
r(j)=(float(ix1)+float(ix2)*rm2)*rm1