]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | FUNCTION ROMG(X) | |
6 | C ********This gives the eikonal function from a table | |
7 | C calculated in the first call | |
8 | DIMENSION FR(0:1000) | |
9 | DATA I0/0/ | |
10 | SAVE | |
11 | IF(I0.NE.0) GO TO 100 | |
12 | DO 50 I=1,1001 | |
13 | XR=(I-1)*0.01 | |
14 | FR(I-1)=OMG0(XR) | |
15 | 50 CONTINUE | |
16 | 100 I0=1 | |
17 | IF(X.GE.10.0) THEN | |
18 | ROMG=0.0 | |
19 | RETURN | |
20 | ENDIF | |
21 | IX=INT(X*100) | |
22 | ROMG=(FR(IX)*((IX+1)*0.01-X)+FR(IX+1)*(X-IX*0.01))/0.01 | |
23 | RETURN | |
24 | END |