HIJING/hijing1_36/hirnd2.F

   1 * $Id$
   2 C
   3 C
   4 C
   5 C
   6 C       This generate random number between XMIN and XMAX
   7         FUNCTION HIRND2(I,XMIN,XMAX)
   8 #define BLANKET_SAVE
   9 #include "hijhb.inc"
  10         SAVE
  11         IF(XMIN.LT.XX(I,1)) XMIN=XX(I,1)
  12         IF(XMAX.GT.XX(I,201)) XMAX=XX(I,201)
  13         JMIN=1+200*(XMIN-XX(I,1))/(XX(I,201)-XX(I,1))
  14         JMAX=1+200*(XMAX-XX(I,1))/(XX(I,201)-XX(I,1))
  15         RX=RR(I,JMIN)+(RR(I,JMAX)-RR(I,JMIN))*RLU_HIJING(0)
  16         JL=0
  17         JU=202
  18 10      IF(JU-JL.GT.1) THEN
  19            JM=(JU+JL)/2
  20            IF((RR(I,201).GT.RR(I,1)).EQV.(RX.GT.RR(I,JM))) THEN
  21               JL=JM
  22            ELSE
  23               JU=JM
  24            ENDIF
  25         GO TO 10
  26         ENDIF
  27         J=JL
  28         IF(J.LT.1) J=1
  29         IF(J.GE.201) J=200
  30         HIRND2=(XX(I,J)+XX(I,J+1))/2.0
  31         RETURN
  32         END