HIJING/hijing1_36/hirnd.F

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