hijing/src/romg.f

0001 C
0002 C
0003 C
0004         FUNCTION ROMG(X)
0005 C               ********This gives the eikonal function from a table
0006 C                       calculated in the first call
0007         DIMENSION FR(0:1000)
0008         DATA I0/0/
0009
0010         SAVE I0, FR                             !Uzhi
0011
0012         IF(I0.NE.0) GO TO 100
0013         DO 50 I=1,1001
0014         XR=(I-1)*0.01
0015         FR(I-1)=OMG0(XR)
0016 50      CONTINUE
0017 100     I0=1
0018         IF(X.GE.10.0) THEN
0019                 ROMG=0.0
0020                 RETURN
0021         ENDIF
0022         IX=INT(X*100)
0023         ROMG=(FR(IX)*((IX+1)*0.01-X)+FR(IX+1)*(X-IX*0.01))/0.01
0024         RETURN
0025         END