aboutsummaryrefslogtreecommitdiff
path: root/synthe/airtovac.for
diff options
context:
space:
mode:
authorJoseph Hunkeler <jhunkeler@gmail.com>2015-02-16 12:40:45 -0500
committerJoseph Hunkeler <jhunkeler@gmail.com>2015-02-16 12:40:45 -0500
commit0373ffdfaaaa3845306ca71243d535fdffd941d4 (patch)
tree194c3c278d7e352e39d555d31aae93c0be2dfc03 /synthe/airtovac.for
parent01b51f73bd06b2d6eabb776ba6cc69e4abfaa0b3 (diff)
downloadkasym-0373ffdfaaaa3845306ca71243d535fdffd941d4.tar.gz
Initial commit
Diffstat (limited to 'synthe/airtovac.for')
-rw-r--r--synthe/airtovac.for30
1 files changed, 30 insertions, 0 deletions
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