HIJING/hijing1_36/hirnd.F

   1 * $Id$
   2 C
   3 C
   4 C
   5         FUNCTION HIRND(I)
   6 #include "hijhb.inc"
   7         SAVE
   8         RX=RLU_HIJING(0)
   9         JL=0
  10         JU=202
  11 10      IF(JU-JL.GT.1) THEN
  12            JM=(JU+JL)/2
  13            IF((RR(I,201).GT.RR(I,1)).EQV.(RX.GT.RR(I,JM))) THEN
  14               JL=JM
  15            ELSE
  16               JU=JM
  17            ENDIF
  18         GO TO 10
  19         ENDIF
  20         J=JL
  21         IF(J.LT.1) J=1
  22         IF(J.GE.201) J=200
  23         HIRND=(XX(I,J)+XX(I,J+1))/2.0
  24         RETURN
  25         END