38 INTEGER a1,a1vw,a1w,a2,a2vw,a2w,m1,m2
41 INTEGER cg1(numg),cg2(numg),ig1(numg),ig2(numg),lg1(numg),
52 EXTERNAL mltmod,qrgnin
58 COMMON /globe/m1,m2,a1,a2,a1w,a2w,a1vw,a2vw,ig1,ig2,lg1,lg2,cg1,
77 IF (.NOT. (qrgnin())) CALL
inrgcm()
82 ig1(g) = mltmod(a1vw,ig1(g-1),m1)
83 ig2(g) = mltmod(a2vw,ig2(g-1),m2)
subroutine setall(iseed1, iseed2)
subroutine initgn(isdtyp)