39 INTEGER a1,a1vw,a1w,a2,a2vw,a2w,m1,m2
42 INTEGER cg1(numg),cg2(numg),ig1(numg),ig2(numg),lg1(numg),
52 EXTERNAL qrgnin,mltmod
58 COMMON /globe/m1,m2,a1,a2,a1w,a2w,a1vw,a2vw,ig1,ig2,lg1,lg2,cg1,
66 IF (qrgnin()) go to 10
67 WRITE (*,*)
' INITGN called before random number generator ',
68 +
' initialized -- abort!'
70 + (
' INITGN called before random number generator initialized')
73 IF ((-1).NE. (isdtyp)) go to 20
78 20
IF ((0).NE. (isdtyp)) go to 30
82 30
IF ((1).NE. (isdtyp)) go to 40
83 lg1(g) = mltmod(a1w,lg1(g),m1)
84 lg2(g) = mltmod(a2w,lg2(g),m2)
87 40 CALL xstopx(
'ISDTYP NOT IN RANGE')
subroutine initgn(isdtyp)