[3] | 1 | FUNCTION REFLIQ(W) |
---|
| 2 | DIMENSION WAVENO(55),XIMG(55) |
---|
| 3 | DATA WAVENO/0., 10., 40., 60., 80., 100., 120., 140., |
---|
| 4 | &160., 180., 200., 220., 240., 260., 280., 300., 325., |
---|
| 5 | &340., 350., 360., 380., 400., 450., 480., 500., 520., |
---|
| 6 | &540., 560., 600., 640., 684., 720., 760., 780., 800., |
---|
| 7 | &840., 875., 920., 960., 1000., 1040., 1080., 1120., 1160., |
---|
| 8 | &1200., 1220., 1240., 1260., 1280., 1300., 1320., 1340., |
---|
| 9 | &1400., 1800., 2000./ |
---|
| 10 | DATA XIMG/1.0E-5, 5.0E-5, 1.1E-3, 1.5E-3, 1.9E-3, 2.3E-3, |
---|
| 11 | &2.4E-3, 2.5E-3, 2.4E-3, 2.2E-3, 1.8E-3, 1.3E-3, 1.1E-3, |
---|
| 12 | &8.6E-4, 6.8E-4, 5.0E-4, 4.0E-4, 2.9E-4, 2.0E-4, 1.8E-4, |
---|
| 13 | &1.2E-4, 5.0E-5, 2.0E-5, 2.1E-5, 2.3E-5, 2.7E-5, 3.1E-5, |
---|
| 14 | &3.5E-5, 5.0E-5, 6.2E-5, 9.0E-5, 1.2E-4, 1.7E-4, 2.0E-4, |
---|
| 15 | &2.4E-4, 3.3E-4, 4.5E-4, 6.2E-4, 9.0E-4, 1.3E-3, 2.0E-3, |
---|
| 16 | &2.9E-3, 4.4E-3, 6.6E-3, 1.0E-2, 1.3E-2, 1.9E-2, 3.0E-2, |
---|
| 17 | &7.0E-2, 1.6E-1, 7.0E-2, 6.0E-3, 6.0E-3, 6.0E-3, 6.0E-3/ |
---|
| 18 | DO 100 I=2,55 |
---|
| 19 | IF (W .GT. WAVENO(I)) GO TO 100 |
---|
| 20 | FACTOR= (WAVENO(I) - W )/(WAVENO(I) - WAVENO(I-1)) |
---|
| 21 | REFLIQ=XIMG(I) + FACTOR*(XIMG(I-1) - XIMG(I)) |
---|
| 22 | RETURN |
---|
| 23 | 100 CONTINUE |
---|
| 24 | REFLIQ=XIMG(55) |
---|
| 25 | RETURN |
---|
| 26 | END |
---|