Merge branch 'prerelease-3.2.1' into czarek

[unres.git] / source / unres / src_CSA_DiL / randgens.f

diff --git a/source/unres/src_CSA_DiL/randgens.f b/source/unres/src_CSA_DiL/randgens.f
deleted file mode 100644 (file)
index 0daeb35..0000000
--- a/source/unres/src_CSA_DiL/randgens.f
+++ /dev/null
@@ -1,99 +0,0 @@
-C $Date: 1994/10/04 16:19:52 $
-C $Revision: 2.1 $
-C
-C
-C  See help for RANDOMV on the PSFSHARE disk to understand these
-C  subroutines.  This is the VS Fortran version of this code.
-C
-C
-      SUBROUTINE VRND(VEC,N)
-      INTEGER A(250)
-      COMMON /VRANDD/ A, I, I147
-      INTEGER LOOP,I,I147,VEC(N)
-      DO 23000 LOOP=1,N
-      I=I+1
-      IF(.NOT.(I.GE.251))GOTO 23002
-      I=1
-23002 CONTINUE
-      I147=I147+1
-      IF(.NOT.(I147.GE.251))GOTO 23004
-      I147=1
-23004 CONTINUE
-      A(I)=IEOR(A(I147),A(I))
-      VEC(LOOP)=A(I)
-23000 CONTINUE
-      RETURN
-      END
-C
-C
-      DOUBLE PRECISION FUNCTION RNDV(IDUM)
-      DOUBLE PRECISION RM1,RM2,R(99)
-      INTEGER IA1,IC1,M1, IA2,IC2,M2, IA3,IC3,M3, IDUM
-      SAVE
-      DATA IA1,IC1,M1/1279,351762,1664557/
-      DATA IA2,IC2,M2/2011,221592,1048583/
-      DATA IA3,IC3,M3/15551,6150,29101/
-      IF(.NOT.(IDUM.LT.0))GOTO 23006
-      IX1 = MOD(-IDUM,M1)
-      IX1 = MOD(IA1*IX1+IC1,M1)
-      IX2 = MOD(IX1,M2)
-      IX1 = MOD(IA1*IX1+IC1,M1)
-      IX3 = MOD(IX1,M3)
-      RM1 = 1./DBLE(M1)
-      RM2 = 1./DBLE(M2)
-      DO 23008 J = 1,99
-      IX1 = MOD(IA1*IX1+IC1,M1)
-      IX2 = MOD(IA2*IX2+IC2,M2)
-      R(J) = (DBLE(IX1)+DBLE(IX2)*RM2)*RM1
-23008 CONTINUE
-23006 CONTINUE
-      IX1 = MOD(IA1*IX1+IC1,M1)
-      IX2 = MOD(IA2*IX2+IC2,M2)
-      IX3 = MOD(IA3*IX3+IC3,M3)
-      J = 1+(99*IX3)/M3
-      RNDV = R(J)
-      R(J) = (DBLE(IX1)+DBLE(IX2)*RM2)*RM1
-      IDUM = IX1
-      RETURN
-      END
-C
-C
-      SUBROUTINE VRNDST(SEED)
-      INTEGER A(250),LOOP,IDUM,SEED
-      DOUBLE PRECISION RNDV
-      COMMON /VRANDD/ A, I, I147
-      I=0
-      I147=103
-      IDUM=SEED
-      DO 23010 LOOP=1,250
-      A(LOOP)=INT(RNDV(IDUM)*2147483647)
-23010 CONTINUE
-      RETURN
-      END
-C
-C
-      SUBROUTINE VRNDIN(IODEV)
-      INTEGER IODEV, A(250)
-      COMMON/VRANDD/ A, I, I147
-      READ(IODEV) A, I, I147
-      RETURN
-      END
-C
-C
-      SUBROUTINE VRNDOU(IODEV)
-C       This corresponds to VRNDOUT in the APFTN64 version
-      INTEGER IODEV, A(250)
-      COMMON/VRANDD/ A, I, I147
-      WRITE(IODEV) A, I, I147
-      RETURN
-      END
-      FUNCTION RNUNF(N)
-      INTEGER IRAN1(2000)
-      DATA FCTOR /2147483647.0D0/
-C     We get only one random number, here!    DR  9/1/92
-      CALL VRND(IRAN1,1)
-      RNUNF= DBLE( IRAN1(1) ) / FCTOR
-C******************************
-C     write(6,*) 'rnunf  in rnunf = ',rnunf
-      RETURN
-      END