hijing/src/hirnd2.f

0001 C
0002 C
0003 C
0004 C
0005 C       This generate random number between XMIN and XMAX
0006         FUNCTION HIRND2(I,XMIN,XMAX)
0007         COMMON/HIJHB/RR(10,201),XX(10,201)
0008         SAVE  /HIJHB/
0009         COMMON/RANSEED/NSEED
0010         SAVE  /RANSEED/
0011         IF(XMIN.LT.XX(I,1)) XMIN=XX(I,1)
0012         IF(XMAX.GT.XX(I,201)) XMAX=XX(I,201)
0013         JMIN=1+200*(XMIN-XX(I,1))/(XX(I,201)-XX(I,1))
0014         JMAX=1+200*(XMAX-XX(I,1))/(XX(I,201)-XX(I,1))
0015         RX=RR(I,JMIN)+(RR(I,JMAX)-RR(I,JMIN))*ATL_RAN(NSEED)
0016         JL=0
0017         JU=202
0018 10      IF(JU-JL.GT.1) THEN
0019            JM=(JU+JL)/2
0020            IF((RR(I,201).GT.RR(I,1)).EQV.(RX.GT.RR(I,JM))) THEN
0021               JL=JM
0022            ELSE
0023               JU=JM
0024            ENDIF
0025         GO TO 10
0026         ENDIF
0027         J=JL
0028         IF(J.LT.1) J=1
0029         IF(J.GE.201) J=200
0030         HIRND2=(XX(I,J)+XX(I,J+1))/2.0
0031         RETURN
0032         END