pythia6416.f.patch

--- pythia6416.f.orig	2008-03-14 09:37:15.000000000 +0100
+++ pythia6416.f	2011-02-03 20:03:16.622831025 +0100
@@ -70354,8 +70354,24 @@
 C...PYR
 C...Generates random numbers uniformly distributed between
 C...0 and 1, excluding the endpoints.
+C...Calls back to an externally defined function EXT_PYR
  
       FUNCTION PYR(IDUMMY)
+      DOUBLE PRECISION PYR, EXT_PYR
+      INTEGER IDUMMY
+  100 PYR=EXT_PYR(IDUMMY)
+      IF(PYR.LE.0D0.OR.PYR.GE.1D0) GOTO 100
+      RETURN
+      END 
+ 
+C*********************************************************************
+ 
+C...OLD_PYR
+C...Generates random numbers uniformly distributed between
+C...0 and 1, excluding the endpoints.
+C...This is the original Pythia RNG
+ 
+      FUNCTION OLD_PYR(IDUMMY)
  
 C...Double precision and integer declarations.
       IMPLICIT DOUBLE PRECISION(A-H, O-Z)
@@ -70424,7 +70440,7 @@
         MRPY2=MRPY2+1
         MRPY3=0
       ENDIF
-      PYR=RUNI
+      OLD_PYR=RUNI
  
       RETURN
       END