SUBROUTINE AIR2 ( FSET, DISP, SPACING, evapor )
C
C   Calculate the appropriate dispersions for the two color method.
C   Modified for use with the e-option in cfi; by C. Hummel
C
      IMPLICIT UNDEFINED (A-Z)
      INTEGER *4 I, J, K, FSET
      REAL *8 WAVE(2), WREF(2), T, P, E, INDEX(2), S1(2), S2(2)
      real*8 evapor
      REAL *8 LAMBDA(4), FR1, FR2
      REAL *4 SPACING(9,2), DISP(4)
      DATA WAVE / .550, .500 /
      DATA WREF / .700, .800 /
      DATA T, P, E / 273., 1013., 0. /
      write(6,'(a,5f8.3)') ' Wavelengths are: ', wave, wref(fset)
      e=evapor
      CALL AINDEX ( P, T, E, WREF(FSET), INDEX(2) )
      DO I = 1, 2
          CALL AINDEX ( P, T, E, WAVE(I), INDEX(1) )
          DISP(I) =  INDEX(1) / ( INDEX(1) - INDEX(2) )
          SPACING(1,I) = 0
          DO J = 1, 2
              FR1 = FLOAT(J) * WREF(FSET)
              DO WHILE ( FR1 .GT. 0. )
                  FR1 = FR1 - WAVE(I)
              END DO
              SPACING(4*J-2,I) =  FR1
              SPACING(4*J-1,I) =  FR1 + WAVE(I)
              SPACING(4*J  ,I) = -FR1
              SPACING(4*J+1,I) = -FR1 - WAVE(I)
          END DO
          DO J = 1, 9
              SPACING(J,I) = SPACING(J,I) * DISP(I)
          END DO
      END DO
      if(fset.eq.2)then
       call aindex ( p, t, e, wref(1), index(1) )
       disp(4) = index(1) / ( index(1) - index(2) )
      endif
      CALL AINDEX ( P, T, E, WAVE(1), INDEX(2) )
      CALL AINDEX ( P, T, E, WAVE(2), INDEX(1) )
      DISP(3) =  INDEX(1) / ( INDEX(1) - INDEX(2) )
      RETURN
      END