Correct ran function.

diff --git a/HIJING/hijing1_36/hijing.F b/HIJING/hijing1_36/hijing.F
index 790d8e5e5a7e550e0e689f30016ae5fcf03971a2..9df2e9215acce7fd7722fba4a6c5f7e9ddf33a27 100644 (file)
--- a/HIJING/hijing1_36/hijing.F
+++ b/HIJING/hijing1_36/hijing.F
@@ -189,9 +189,9 @@ C
 5      R=HIRND(1)
 c
         if(IHNT2(1).EQ.2) then
-           rnd1=max(RAN0(NSEED),1.0e-20)
-           rnd2=max(RAN0(NSEED),1.0e-20)
-           rnd3=max(RAN0(NSEED),1.0e-20)
+           rnd1=max(RLU_HIJING(NSEED),1.0e-20)
+           rnd2=max(RLU_HIJING(NSEED),1.0e-20)
+           rnd3=max(RLU_HIJING(NSEED),1.0e-20)
            R=-0.5*(log(rnd1)*4.38/2.0+log(rnd2)*0.85/2.0
      &          +4.38*0.85*log(rnd3)/(4.38+0.85))
         endif
@@ -246,9 +246,9 @@ C******************************
 15     R=HIRND(2)
 c
          if(IHNT2(3).EQ.2) then
-            rnd1=max(RAN0(NSEED),1.0e-20)
-            rnd2=max(RAN0(NSEED),1.0e-20)
-            rnd3=max(RAN0(NSEED),1.0e-20)
+            rnd1=max(RLU_HIJING(NSEED),1.0e-20)
+            rnd2=max(RLU_HIJING(NSEED),1.0e-20)
+            rnd3=max(RLU_HIJING(NSEED),1.0e-20)
             R=-0.5*(log(rnd1)*4.38/2.0+log(rnd2)*0.85/2.0
      &           +4.38*0.85*log(rnd3)/(4.38+0.85))
          endif
@@ -645,7 +645,7 @@ C       ****** identify the mother particle
                    V3MIN1=RARB/GAMMA0
                    V3MIN2=1.0/MAX(1.0,5.08*ABS(PATT(I,3)))
                    VATT_MIN=MAX(V3MIN1,V3MIN2)
-                   VATT03=VATT03+(0.5-RLU(0))*VATT_MIN
+                   VATT03=VATT03+(0.5-RLU_HIJING(0))*VATT_MIN
                    amt2=P(I,1)**2+P(I,2)**2+P(I,5)**2
                    IF(amt2.GT.0.0) THEN
                       tauf=0.2*2.0*P(I,3)/amt2
@@ -733,7 +733,7 @@ C             of the produced hadrons
                    V3MIN1=RARB/GAMMA0
                    V3MIN2=1.0/MAX(1.0,5.08*ABS(PATT(I,3)))
                    VATT_MIN=MAX(V3MIN1,V3MIN2)
-                   VATT03=VATT03+(0.5-RLU(0))*VATT_MIN
+                   VATT03=VATT03+(0.5-RLU_HIJING(0))*VATT_MIN
                    amt2=P(I,1)**2+P(I,2)**2+P(I,5)**2
                    IF(amt2.GT.0.0) THEN
                       tauf=0.2*2.0*P(I,3)/amt2