* $Id$ C C C C C This generate random number between XMIN and XMAX FUNCTION HIRND2(I,XMIN,XMAX) #include "hijhb.inc" SAVE IF(XMIN.LT.XX(I,1)) XMIN=XX(I,1) IF(XMAX.GT.XX(I,201)) XMAX=XX(I,201) JMIN=1+200*(XMIN-XX(I,1))/(XX(I,201)-XX(I,1)) JMAX=1+200*(XMAX-XX(I,1))/(XX(I,201)-XX(I,1)) RX=RR(I,JMIN)+(RR(I,JMAX)-RR(I,JMIN))*RLU_HIJING(0) JL=0 JU=202 10 IF(JU-JL.GT.1) THEN JM=(JU+JL)/2 IF((RR(I,201).GT.RR(I,1)).EQV.(RX.GT.RR(I,JM))) THEN JL=JM ELSE JU=JM ENDIF GO TO 10 ENDIF J=JL IF(J.LT.1) J=1 IF(J.GE.201) J=200 HIRND2=(XX(I,J)+XX(I,J+1))/2.0 RETURN END