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