| 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 |
|---|
