hijing/src/hirnd.f

0001 C
0002 C
0003 C
0004         FUNCTION HIRND(I)
0005         COMMON/HIJHB/RR(10,201),XX(10,201)
0006         SAVE  /HIJHB/
0007         COMMON/RANSEED/NSEED
0008         SAVE  /RANSEED/
0009         RX=ATL_RAN(NSEED)
0010         JL=0
0011         JU=202
0012 10      IF(JU-JL.GT.1) THEN
0013            JM=(JU+JL)/2
0014            IF((RR(I,201).GT.RR(I,1)).EQV.(RX.GT.RR(I,JM))) THEN
0015               JL=JM
0016            ELSE
0017               JU=JM
0018            ENDIF
0019         GO TO 10
0020         ENDIF
0021         J=JL
0022         IF(J.LT.1) J=1
0023         IF(J.GE.201) J=200
0024         HIRND=(XX(I,J)+XX(I,J+1))/2.0
0025         RETURN
0026         END