X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?p=unres.git;a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fprng_32.F;fp=source%2Funres%2Fsrc_MD-M%2Fprng_32.F;h=9448f312872b217a542976f4f5e8653ac7db7e2e;hp=21cac760c31ffa52a400e3fed907c028d2fc6f7c;hb=2a226bfc86eabc6e4eae0c3ad1cbc3cb5417a05a;hpb=a0e685f844163003749ba91dfbf4644bcc8cfa30 diff --git a/source/unres/src_MD-M/prng_32.F b/source/unres/src_MD-M/prng_32.F index 21cac76..9448f31 100644 --- a/source/unres/src_MD-M/prng_32.F +++ b/source/unres/src_MD-M/prng_32.F @@ -1,7 +1,7 @@ #if defined(AIX) || defined(AMD64) - real*8 function prng_next(me) + real*8 function prng_next(mel) implicit none - integer me + integer me,mel c c Calling sequence: c = prng_next ( ) @@ -47,8 +47,11 @@ c parameter ( ishift = 9007199254740991 ) ! 2**53-1 crc g77 doesn't support integer*8 constants m = dint(34522712143931.0d0) ishift = dint(9007199254740991.0d0) - if(me.gt.nmax) me=mod(me,nmax) - + if(mel.gt.nmax) then + me=mod(mel,nmax) + else + me=mel + endif c RS6K porting note: ishift now takes 64-bit integers , with a warning if ( 0.le.me .and. me.le.nmax ) then next = iparam(1,me)*m + iparam(2,me) @@ -125,9 +128,9 @@ c c me is the particular generator being restarted c iseed is a 8-byte integer containing the "l"-values c - logical function prng_restart (me, iseed) + logical function prng_restart (mel, iseed) implicit none - integer me + integer me,mel integer*8 iseed integer nmax @@ -135,7 +138,11 @@ c parameter(nmax=1021) common/ksrprng/iparam(2,0:nmax) - if(me.gt.nmax) me=mod(me,nmax) + if(mel.gt.nmax) then + me=mod(mel,nmax) + else + me=mel + endif if (me .lt. 0 .or. me .gt. nmax) then prng_restart=.false. return