From 0373ffdfaaaa3845306ca71243d535fdffd941d4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Mon, 16 Feb 2015 12:40:45 -0500
Subject: Initial commit

---
 synthe/airtovac.for | 30 ++++++++++++++++++++++++++++++
 1 file changed, 30 insertions(+)
 create mode 100644 synthe/airtovac.for

(limited to 'synthe/airtovac.for')

diff --git a/synthe/airtovac.for b/synthe/airtovac.for
new file mode 100644
index 0000000..67f78a2
--- /dev/null
+++ b/synthe/airtovac.for
@@ -0,0 +1,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
-- 
cgit