synthe/airtovac.for


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30

FUNCTION AIRTOVAC(W)
      IMPLICIT REAL*8 (A-H,O-Z)
C     W IS AIR WAVELENGTH IN NM
C     WAVEN IS AIR WAVENUMBER WHICH IS USUALLY GOOD ENOUGH 
c
c     Edlen, B. 1966. Metrologia 2, 71.
c     quoted by 
c     Peck, E.R. and Reeder, K. 1972.  JOSA 62, 958-962.
c
C     MUST ITERATE FOR EXACT SOLUTION
      WAVEN=1.D7/W
      WNEW=W*(1.0000834213D0+
     1 2406030.D0/(1.30D10-WAVEN**2.)+15997.D0/(3.89D9-WAVEN**2.))
C    1(1.000064328+2949810./(1.46E10-WAVEN**2)+25540./(4.1E9-WAVEN**2))         
      WAVEN=1.D7/WNEW
      WNEW=W*(1.0000834213D0+
     1 2406030.D0/(1.30D10-WAVEN**2.)+15997.D0/(3.89D9-WAVEN**2.))
      WAVEN=1.D7/WNEW
      AIRTOVAC=W*(1.0000834213D0+
     1 2406030.D0/(1.30D10-WAVEN**2.)+15997.D0/(3.89D9-WAVEN**2.))
      RETURN
      END
      FUNCTION VACTOAIR(W)
      IMPLICIT REAL*8 (A-H,O-Z)
C     W IS VACUUM WAVELENGTH IN NM
      WAVEN=1.D7/W
      VACTOAIR=W/(1.0000834213D0+
     1 2406030.D0/(1.30D10-WAVEN**2)+15997.D0/(3.89D9-WAVEN**2))
      RETURN
      END