diff --git a/src/eps.f90 b/src/eps.f90 index efe4495..4edba0b 100644 --- a/src/eps.f90 +++ b/src/eps.f90 @@ -70,19 +70,18 @@ function AuRaschke(wavelength) result(eps) real(8), dimension(164) :: lambda integer :: i - lambda = (/ 300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 1300, 1310, 1320, 1330, 1340, 1350, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 2000 /) + lambda = (/ 3.00000000E+02, 3.10000000E+02, 3.20000000E+02, 3.30000000E+02, 3.40000000E+02, 3.50000000E+02, 3.60000000E+02, 3.70000000E+02, 3.80000000E+02, 3.90000000E+02, 4.00000000E+02, 4.10000000E+02, 4.20000000E+02, 4.30000000E+02, 4.40000000E+02, 4.50000000E+02, 4.60000000E+02, 4.70000000E+02, 4.80000000E+02, 4.90000000E+02, 5.00000000E+02, 5.10000000E+02, 5.20000000E+02, 5.30000000E+02, 5.40000000E+02, 5.50000000E+02, 5.60000000E+02, 5.70000000E+02, 5.80000000E+02, 5.90000000E+02, 6.00000000E+02, 6.10000000E+02, 6.20000000E+02, 6.30000000E+02, 6.40000000E+02, 6.50000000E+02, 6.60000000E+02, 6.70000000E+02, 6.80000000E+02, 6.90000000E+02, 7.00000000E+02, 7.10000000E+02, 7.20000000E+02, 7.30000000E+02, 7.40000000E+02, 7.50000000E+02, 7.60000000E+02, 7.70000000E+02, 7.80000000E+02, 7.90000000E+02, 8.00000000E+02, 8.10000000E+02, 8.20000000E+02, 8.30000000E+02, 8.40000000E+02, 8.50000000E+02, 8.60000000E+02, 8.70000000E+02, 8.80000000E+02, 8.90000000E+02, 9.00000000E+02, 9.10000000E+02, 9.20000000E+02, 9.30000000E+02, 9.40000000E+02, 9.50000000E+02, 9.60000000E+02, 9.70000000E+02, 9.80000000E+02, 9.90000000E+02, 1.00000000E+03, 1.01000000E+03, 1.02000000E+03, 1.03000000E+03, 1.04000000E+03, 1.05000000E+03, 1.06000000E+03, 1.07000000E+03, 1.08000000E+03, 1.09000000E+03, 1.10000000E+03, 1.11000000E+03, 1.12000000E+03, 1.13000000E+03, 1.14000000E+03, 1.15000000E+03, 1.16000000E+03, 1.17000000E+03, 1.18000000E+03, 1.19000000E+03, 1.20000000E+03, 1.21000000E+03, 1.22000000E+03, 1.23000000E+03, 1.24000000E+03, 1.25000000E+03, 1.26000000E+03, 1.27000000E+03, 1.28000000E+03, 1.29000000E+03, 1.30000000E+03, 1.31000000E+03, 1.32000000E+03, 1.33000000E+03, 1.34000000E+03, 1.35000000E+03, 1.43000000E+03, 1.44000000E+03, 1.45000000E+03, 1.46000000E+03, 1.47000000E+03, 1.48000000E+03, 1.49000000E+03, 1.50000000E+03, 1.51000000E+03, 1.52000000E+03, 1.53000000E+03, 1.54000000E+03, 1.55000000E+03, 1.56000000E+03, 1.57000000E+03, 1.58000000E+03, 1.59000000E+03, 1.60000000E+03, 1.61000000E+03, 1.62000000E+03, 1.63000000E+03, 1.64000000E+03, 1.65000000E+03, 1.66000000E+03, 1.67000000E+03, 1.68000000E+03, 1.69000000E+03, 1.70000000E+03, 1.71000000E+03, 1.72000000E+03, 1.73000000E+03, 1.74000000E+03, 1.75000000E+03, 1.76000000E+03, 1.77000000E+03, 1.78000000E+03, 1.79000000E+03, 1.80000000E+03, 1.81000000E+03, 1.82000000E+03, 1.83000000E+03, 1.84000000E+03, 1.85000000E+03, 1.86000000E+03, 1.87000000E+03, 1.88000000E+03, 1.89000000E+03, 1.90000000E+03, 1.91000000E+03, 1.92000000E+03, 1.93000000E+03, 1.94000000E+03, 1.95000000E+03, 1.96000000E+03, 1.97000000E+03, 1.98000000E+03, 1.99000000E+03, 2.00000000E+03 /) - re_eps = (/ -0.896, -0.781, -0.633, -0.513, -0.459, -0.494, -0.615, -0.788, -0.952, -1.05, -1.08, -1.06, -1.07, -1.13, -1.19, -1.21, -1.19, -1.23, -1.42, -1.82, -2.38, -3.04, -3.74, -4.44, -5.13, -5.81, -6.47, -7.13, -7.79, -8.45, -9.11, -9.77, -10.4, -11.1, -11.8, -12.5, -13.2, -13.9, -14.6, -15.3, -16, -16.7, -17.5, -18.2, -18.9, -19.7, -20.5, -21.2, -22, -22.8, -23.6, -24.4, -25.2, -26, -26.9, -27.7, -28.6, -29.4, -30.3, -31.1, -32, -32.9, -33.8, -34.7, -35.6, -36.6, -37.5, -38.4, -39.4, -40.3, -41.3, -42.3, -43.3, -44.3, -45.3, -46.3, -47.3, -48.3, -49.4, -50.4, -51.5, -52.5, -53.6, -54.7, -55.8, -56.9, -58, -59.1, -60.2, -61.4, -62.5, -63.7, -64.8, -66, -67.2, -68.4, -69.6, -70.8, -72, -73.2, -74.4, -75.7, -76.9, -78.2, -79.5, -80.7, -91.3, -92.6, -94, -95.2, -96.5, -97.7, -99, -100, -102, -103, -104, -105, -107, -108, -109, -110, -112, -113, -114, -116, -117, -118, -120, -121, -122, -124, -125, -126, -128, -129, -130, -132, -133, -134, -136, -137, -138, -140, -141, -142, -144, -145, -147, -148, -149, -151, -152, -154, -155, -156, -158, -159, -161, -162, -163, -165, -166, -168 /) + re_eps = (/ -8.96000000E-01, -7.81000000E-01, -6.33000000E-01, -5.13000000E-01, -4.59000000E-01, -4.94000000E-01, -6.15000000E-01, -7.88000000E-01, -9.52000000E-01, -1.05000000E+00, -1.08000000E+00, -1.06000000E+00, -1.07000000E+00, -1.13000000E+00, -1.19000000E+00, -1.21000000E+00, -1.19000000E+00, -1.23000000E+00, -1.42000000E+00, -1.82000000E+00, -2.38000000E+00, -3.04000000E+00, -3.74000000E+00, -4.44000000E+00, -5.13000000E+00, -5.81000000E+00, -6.47000000E+00, -7.13000000E+00, -7.79000000E+00, -8.45000000E+00, -9.11000000E+00, -9.77000000E+00, -1.04000000E+01, -1.11000000E+01, -1.18000000E+01, -1.25000000E+01, -1.32000000E+01, -1.39000000E+01, -1.46000000E+01, -1.53000000E+01, -1.60000000E+01, -1.67000000E+01, -1.75000000E+01, -1.82000000E+01, -1.89000000E+01, -1.97000000E+01, -2.05000000E+01, -2.12000000E+01, -2.20000000E+01, -2.28000000E+01, -2.36000000E+01, -2.44000000E+01, -2.52000000E+01, -2.60000000E+01, -2.69000000E+01, -2.77000000E+01, -2.86000000E+01, -2.94000000E+01, -3.03000000E+01, -3.11000000E+01, -3.20000000E+01, -3.29000000E+01, -3.38000000E+01, -3.47000000E+01, -3.56000000E+01, -3.66000000E+01, -3.75000000E+01, -3.84000000E+01, -3.94000000E+01, -4.03000000E+01, -4.13000000E+01, -4.23000000E+01, -4.33000000E+01, -4.43000000E+01, -4.53000000E+01, -4.63000000E+01, -4.73000000E+01, -4.83000000E+01, -4.94000000E+01, -5.04000000E+01, -5.15000000E+01, -5.25000000E+01, -5.36000000E+01, -5.47000000E+01, -5.58000000E+01, -5.69000000E+01, -5.80000000E+01, -5.91000000E+01, -6.02000000E+01, -6.14000000E+01, -6.25000000E+01, -6.37000000E+01, -6.48000000E+01, -6.60000000E+01, -6.72000000E+01, -6.84000000E+01, -6.96000000E+01, -7.08000000E+01, -7.20000000E+01, -7.32000000E+01, -7.44000000E+01, -7.57000000E+01, -7.69000000E+01, -7.82000000E+01, -7.95000000E+01, -8.07000000E+01, -9.13000000E+01, -9.26000000E+01, -9.40000000E+01, -9.52000000E+01, -9.65000000E+01, -9.77000000E+01, -9.90000000E+01, -1.00000000E+02, -1.02000000E+02, -1.03000000E+02, -1.04000000E+02, -1.05000000E+02, -1.07000000E+02, -1.08000000E+02, -1.09000000E+02, -1.10000000E+02, -1.12000000E+02, -1.13000000E+02, -1.14000000E+02, -1.16000000E+02, -1.17000000E+02, -1.18000000E+02, -1.20000000E+02, -1.21000000E+02, -1.22000000E+02, -1.24000000E+02, -1.25000000E+02, -1.26000000E+02, -1.28000000E+02, -1.29000000E+02, -1.30000000E+02, -1.32000000E+02, -1.33000000E+02, -1.34000000E+02, -1.36000000E+02, -1.37000000E+02, -1.38000000E+02, -1.40000000E+02, -1.41000000E+02, -1.42000000E+02, -1.44000000E+02, -1.45000000E+02, -1.47000000E+02, -1.48000000E+02, -1.49000000E+02, -1.51000000E+02, -1.52000000E+02, -1.54000000E+02, -1.55000000E+02, -1.56000000E+02, -1.58000000E+02, -1.59000000E+02, -1.61000000E+02, -1.62000000E+02, -1.63000000E+02, -1.65000000E+02, -1.66000000E+02, -1.68000000E+02 /) - im_eps = (/ 5.71, 5.91, 5.99, 5.97, 5.89, 5.77, 5.67, 5.64, 5.68, 5.76, 5.81, 5.78, 5.67, 5.54, 5.4, 5.24, 4.96, 4.52, 3.95, 3.37, 2.86, 2.45, 2.16, 1.95, 1.79, 1.67, 1.58, 1.5, 1.43, 1.37, 1.32, 1.28, 1.24, 1.21, 1.18, 1.16, 1.14, 1.13, 1.12, 1.11, 1.1, 1.1, 1.11, 1.11, 1.12, 1.13, 1.14, 1.15, 1.17, 1.18, 1.2, 1.22, 1.25, 1.27, 1.3, 1.32, 1.35, 1.38, 1.42, 1.45, 1.48, 1.52, 1.56, 1.59, 1.63, 1.67, 1.71, 1.76, 1.8, 1.85, 1.89, 1.94, 1.99, 2.04, 2.09, 2.14, 2.19, 2.25, 2.3, 2.36, 2.42, 2.47, 2.53, 2.59, 2.65, 2.72, 2.78, 2.84, 2.91, 2.98, 3.04, 3.11, 3.18, 3.25, 3.32, 3.4, 3.47, 3.54, 3.62, 3.7, 3.77, 3.85, 3.93, 4.01, 4.1, 4.18, 4.88, 4.97, 5.06, 5.16, 5.26, 5.37, 5.47, 5.57, 5.68, 5.79, 5.9, 6.01, 6.13, 6.24, 6.36, 6.48, 6.6, 6.72, 6.85, 6.98, 7.11, 7.24, 7.37, 7.51, 7.65, 7.79, 7.93, 8.07, 8.22, 8.37, 8.52, 8.68, 8.84, 9, 9.16, 9.33, 9.5, 9.67, 9.84, 10, 10.2, 10.4, 10.6, 10.8, 11, 11.1, 11.3, 11.5, 11.8, 12, 12.2, 12.4, 12.6, 12.8, 13.1, 13.3, 13.5, 13.7 /) + im_eps = (/ 5.71000000E+00, 5.91000000E+00, 5.99000000E+00, 5.97000000E+00, 5.89000000E+00, 5.77000000E+00, 5.67000000E+00, 5.64000000E+00, 5.68000000E+00, 5.76000000E+00, 5.81000000E+00, 5.78000000E+00, 5.67000000E+00, 5.54000000E+00, 5.40000000E+00, 5.24000000E+00, 4.96000000E+00, 4.52000000E+00, 3.95000000E+00, 3.37000000E+00, 2.86000000E+00, 2.45000000E+00, 2.16000000E+00, 1.95000000E+00, 1.79000000E+00, 1.67000000E+00, 1.58000000E+00, 1.50000000E+00, 1.43000000E+00, 1.37000000E+00, 1.32000000E+00, 1.28000000E+00, 1.24000000E+00, 1.21000000E+00, 1.18000000E+00, 1.16000000E+00, 1.14000000E+00, 1.13000000E+00, 1.12000000E+00, 1.11000000E+00, 1.10000000E+00, 1.10000000E+00, 1.11000000E+00, 1.11000000E+00, 1.12000000E+00, 1.13000000E+00, 1.14000000E+00, 1.15000000E+00, 1.17000000E+00, 1.18000000E+00, 1.20000000E+00, 1.22000000E+00, 1.25000000E+00, 1.27000000E+00, 1.30000000E+00, 1.32000000E+00, 1.35000000E+00, 1.38000000E+00, 1.42000000E+00, 1.45000000E+00, 1.48000000E+00, 1.52000000E+00, 1.56000000E+00, 1.59000000E+00, 1.63000000E+00, 1.67000000E+00, 1.71000000E+00, 1.76000000E+00, 1.80000000E+00, 1.85000000E+00, 1.89000000E+00, 1.94000000E+00, 1.99000000E+00, 2.04000000E+00, 2.09000000E+00, 2.14000000E+00, 2.19000000E+00, 2.25000000E+00, 2.30000000E+00, 2.36000000E+00, 2.42000000E+00, 2.47000000E+00, 2.53000000E+00, 2.59000000E+00, 2.65000000E+00, 2.72000000E+00, 2.78000000E+00, 2.84000000E+00, 2.91000000E+00, 2.98000000E+00, 3.04000000E+00, 3.11000000E+00, 3.18000000E+00, 3.25000000E+00, 3.32000000E+00, 3.40000000E+00, 3.47000000E+00, 3.54000000E+00, 3.62000000E+00, 3.70000000E+00, 3.77000000E+00, 3.85000000E+00, 3.93000000E+00, 4.01000000E+00, 4.10000000E+00, 4.18000000E+00, 4.88000000E+00, 4.97000000E+00, 5.06000000E+00, 5.16000000E+00, 5.26000000E+00, 5.37000000E+00, 5.47000000E+00, 5.57000000E+00, 5.68000000E+00, 5.79000000E+00, 5.90000000E+00, 6.01000000E+00, 6.13000000E+00, 6.24000000E+00, 6.36000000E+00, 6.48000000E+00, 6.60000000E+00, 6.72000000E+00, 6.85000000E+00, 6.98000000E+00, 7.11000000E+00, 7.24000000E+00, 7.37000000E+00, 7.51000000E+00, 7.65000000E+00, 7.79000000E+00, 7.93000000E+00, 8.07000000E+00, 8.22000000E+00, 8.37000000E+00, 8.52000000E+00, 8.68000000E+00, 8.84000000E+00, 9.00000000E+00, 9.16000000E+00, 9.33000000E+00, 9.50000000E+00, 9.67000000E+00, 9.84000000E+00, 1.00000000E+01, 1.02000000E+01, 1.04000000E+01, 1.06000000E+01, 1.08000000E+01, 1.10000000E+01, 1.11000000E+01, 1.13000000E+01, 1.15000000E+01, 1.18000000E+01, 1.20000000E+01, 1.22000000E+01, 1.24000000E+01, 1.26000000E+01, 1.28000000E+01, 1.31000000E+01, 1.33000000E+01, 1.35000000E+01, 1.37000000E+01 /) do i = 1, 164 eps_(i) = complex(re_eps(i), im_eps(i)) end do call interp1(lambda, eps_, wavelength, eps) -end function AuRaschke - +end function epsAuRaschke function epsAg(wavelength) result(eps) ! @@ -119,18 +118,18 @@ function AgRaschke(wavelength) result(eps) real(8), dimension(167) :: lambda integer :: i - lambda = (/ 270, 280, 290, 300, 310, 320, 330, 340, 350, 360, 370, 380, 390, 400, 410, 420, 430, 440, 450, 460, 470, 480, 490, 500, 510, 520, 530, 540, 550, 560, 570, 580, 590, 600, 610, 620, 630, 640, 650, 660, 670, 680, 690, 700, 710, 720, 730, 740, 750, 760, 770, 780, 790, 800, 810, 820, 830, 840, 850, 860, 870, 880, 890, 900, 910, 920, 930, 940, 950, 960, 970, 980, 990, 1000, 1010, 1020, 1030, 1040, 1050, 1060, 1070, 1080, 1090, 1100, 1110, 1120, 1130, 1140, 1150, 1160, 1170, 1180, 1190, 1200, 1210, 1220, 1230, 1240, 1250, 1260, 1270, 1280, 1290, 1300, 1310, 1320, 1330, 1340, 1350, 1430, 1440, 1450, 1460, 1470, 1480, 1490, 1500, 1510, 1520, 1530, 1540, 1550, 1560, 1570, 1580, 1590, 1600, 1610, 1620, 1630, 1640, 1650, 1660, 1670, 1680, 1690, 1700, 1710, 1720, 1730, 1740, 1750, 1760, 1770, 1780, 1790, 1800, 1810, 1820, 1830, 1840, 1850, 1860, 1870, 1880, 1890, 1900, 1910, 1920, 1930, 1940, 1950, 1960, 1970, 1980, 1990, 2000 /) + lambda = (/ 2.70000000E+02, 2.80000000E+02, 2.90000000E+02, 3.00000000E+02, 3.10000000E+02, 3.20000000E+02, 3.30000000E+02, 3.40000000E+02, 3.50000000E+02, 3.60000000E+02, 3.70000000E+02, 3.80000000E+02, 3.90000000E+02, 4.00000000E+02, 4.10000000E+02, 4.20000000E+02, 4.30000000E+02, 4.40000000E+02, 4.50000000E+02, 4.60000000E+02, 4.70000000E+02, 4.80000000E+02, 4.90000000E+02, 5.00000000E+02, 5.10000000E+02, 5.20000000E+02, 5.30000000E+02, 5.40000000E+02, 5.50000000E+02, 5.60000000E+02, 5.70000000E+02, 5.80000000E+02, 5.90000000E+02, 6.00000000E+02, 6.10000000E+02, 6.20000000E+02, 6.30000000E+02, 6.40000000E+02, 6.50000000E+02, 6.60000000E+02, 6.70000000E+02, 6.80000000E+02, 6.90000000E+02, 7.00000000E+02, 7.10000000E+02, 7.20000000E+02, 7.30000000E+02, 7.40000000E+02, 7.50000000E+02, 7.60000000E+02, 7.70000000E+02, 7.80000000E+02, 7.90000000E+02, 8.00000000E+02, 8.10000000E+02, 8.20000000E+02, 8.30000000E+02, 8.40000000E+02, 8.50000000E+02, 8.60000000E+02, 8.70000000E+02, 8.80000000E+02, 8.90000000E+02, 9.00000000E+02, 9.10000000E+02, 9.20000000E+02, 9.30000000E+02, 9.40000000E+02, 9.50000000E+02, 9.60000000E+02, 9.70000000E+02, 9.80000000E+02, 9.90000000E+02, 1.00000000E+03, 1.01000000E+03, 1.02000000E+03, 1.03000000E+03, 1.04000000E+03, 1.05000000E+03, 1.06000000E+03, 1.07000000E+03, 1.08000000E+03, 1.09000000E+03, 1.10000000E+03, 1.11000000E+03, 1.12000000E+03, 1.13000000E+03, 1.14000000E+03, 1.15000000E+03, 1.16000000E+03, 1.17000000E+03, 1.18000000E+03, 1.19000000E+03, 1.20000000E+03, 1.21000000E+03, 1.22000000E+03, 1.23000000E+03, 1.24000000E+03, 1.25000000E+03, 1.26000000E+03, 1.27000000E+03, 1.28000000E+03, 1.29000000E+03, 1.30000000E+03, 1.31000000E+03, 1.32000000E+03, 1.33000000E+03, 1.34000000E+03, 1.35000000E+03, 1.43000000E+03, 1.44000000E+03, 1.45000000E+03, 1.46000000E+03, 1.47000000E+03, 1.48000000E+03, 1.49000000E+03, 1.50000000E+03, 1.51000000E+03, 1.52000000E+03, 1.53000000E+03, 1.54000000E+03, 1.55000000E+03, 1.56000000E+03, 1.57000000E+03, 1.58000000E+03, 1.59000000E+03, 1.60000000E+03, 1.61000000E+03, 1.62000000E+03, 1.63000000E+03, 1.64000000E+03, 1.65000000E+03, 1.66000000E+03, 1.67000000E+03, 1.68000000E+03, 1.69000000E+03, 1.70000000E+03, 1.71000000E+03, 1.72000000E+03, 1.73000000E+03, 1.74000000E+03, 1.75000000E+03, 1.76000000E+03, 1.77000000E+03, 1.78000000E+03, 1.79000000E+03, 1.80000000E+03, 1.81000000E+03, 1.82000000E+03, 1.83000000E+03, 1.84000000E+03, 1.85000000E+03, 1.86000000E+03, 1.87000000E+03, 1.88000000E+03, 1.89000000E+03, 1.90000000E+03, 1.91000000E+03, 1.92000000E+03, 1.93000000E+03, 1.94000000E+03, 1.95000000E+03, 1.96000000E+03, 1.97000000E+03, 1.98000000E+03, 1.99000000E+03, 2.00000000E+03 /) - re_eps = (/ 0.165, 0.514, 0.985, 1.7, 1.62, 0.561, -0.392, -1.13, -1.7, -2.25, -2.77, -3.28, -3.78, -4.27, -4.77, -5.26, -5.77, -6.27, -6.78, -7.31, -7.84, -8.38, -8.92, -9.47, -10, -10.6, -11.2, -11.8, -12.4, -13, -13.6, -14.2, -14.9, -15.5, -16.1, -16.8, -17.5, -18.2, -18.9, -19.6, -20.3, -21, -21.7, -22.5, -23.2, -24, -24.8, -25.6, -26.4, -27.2, -28, -28.8, -29.6, -30.5, -31.4, -32.3, -33.1, -34, -35.2, -36.1, -37.1, -37.9, -39.1, -40, -41.1, -41.7, -42.8, -43.8, -44.8, -45.9, -46.9, -47.7, -49, -50, -51.2, -52.1, -53.3, -54.4, -55.5, -56.6, -57.8, -59, -60, -61.3, -62.5, -63.6, -64.9, -66.1, -67.3, -68.6, -70, -71, -72.2, -73.6, -74.8, -76.2, -77.4, -78.8, -80.1, -81.7, -82.7, -84.1, -85.6, -87, -88.4, -89.9, -91.2, -92.7, -94.1, -106, -108, -109, -111, -113, -114, -115, -117, -119, -120, -122, -124, -125, -127, -129, -130, -132, -134, -135, -137, -139, -141, -142, -144, -146, -148, -149, -151, -153, -155, -157, -159, -160, -162, -164, -166, -168, -170, -172, -175, -176, -179, -180, -182, -185, -187, -188, -190, -192, -195, -197, -198, -201, -203, -205, -208, -209, -211 /) + re_eps = (/ 1.65000000E-01, 5.14000000E-01, 9.85000000E-01, 1.70000000E+00, 1.62000000E+00, 5.61000000E-01, -3.92000000E-01, -1.13000000E+00, -1.70000000E+00, -2.25000000E+00, -2.77000000E+00, -3.28000000E+00, -3.78000000E+00, -4.27000000E+00, -4.77000000E+00, -5.26000000E+00, -5.77000000E+00, -6.27000000E+00, -6.78000000E+00, -7.31000000E+00, -7.84000000E+00, -8.38000000E+00, -8.92000000E+00, -9.47000000E+00, -1.00000000E+01, -1.06000000E+01, -1.12000000E+01, -1.18000000E+01, -1.24000000E+01, -1.30000000E+01, -1.36000000E+01, -1.42000000E+01, -1.49000000E+01, -1.55000000E+01, -1.61000000E+01, -1.68000000E+01, -1.75000000E+01, -1.82000000E+01, -1.89000000E+01, -1.96000000E+01, -2.03000000E+01, -2.10000000E+01, -2.17000000E+01, -2.25000000E+01, -2.32000000E+01, -2.40000000E+01, -2.48000000E+01, -2.56000000E+01, -2.64000000E+01, -2.72000000E+01, -2.80000000E+01, -2.88000000E+01, -2.96000000E+01, -3.05000000E+01, -3.14000000E+01, -3.23000000E+01, -3.31000000E+01, -3.40000000E+01, -3.52000000E+01, -3.61000000E+01, -3.71000000E+01, -3.79000000E+01, -3.91000000E+01, -4.00000000E+01, -4.11000000E+01, -4.17000000E+01, -4.28000000E+01, -4.38000000E+01, -4.48000000E+01, -4.59000000E+01, -4.69000000E+01, -4.77000000E+01, -4.90000000E+01, -5.00000000E+01, -5.12000000E+01, -5.21000000E+01, -5.33000000E+01, -5.44000000E+01, -5.55000000E+01, -5.66000000E+01, -5.78000000E+01, -5.90000000E+01, -6.00000000E+01, -6.13000000E+01, -6.25000000E+01, -6.36000000E+01, -6.49000000E+01, -6.61000000E+01, -6.73000000E+01, -6.86000000E+01, -7.00000000E+01, -7.10000000E+01, -7.22000000E+01, -7.36000000E+01, -7.48000000E+01, -7.62000000E+01, -7.74000000E+01, -7.88000000E+01, -8.01000000E+01, -8.17000000E+01, -8.27000000E+01, -8.41000000E+01, -8.56000000E+01, -8.70000000E+01, -8.84000000E+01, -8.99000000E+01, -9.12000000E+01, -9.27000000E+01, -9.41000000E+01, -1.06000000E+02, -1.08000000E+02, -1.09000000E+02, -1.11000000E+02, -1.13000000E+02, -1.14000000E+02, -1.15000000E+02, -1.17000000E+02, -1.19000000E+02, -1.20000000E+02, -1.22000000E+02, -1.24000000E+02, -1.25000000E+02, -1.27000000E+02, -1.29000000E+02, -1.30000000E+02, -1.32000000E+02, -1.34000000E+02, -1.35000000E+02, -1.37000000E+02, -1.39000000E+02, -1.41000000E+02, -1.42000000E+02, -1.44000000E+02, -1.46000000E+02, -1.48000000E+02, -1.49000000E+02, -1.51000000E+02, -1.53000000E+02, -1.55000000E+02, -1.57000000E+02, -1.59000000E+02, -1.60000000E+02, -1.62000000E+02, -1.64000000E+02, -1.66000000E+02, -1.68000000E+02, -1.70000000E+02, -1.72000000E+02, -1.75000000E+02, -1.76000000E+02, -1.79000000E+02, -1.80000000E+02, -1.82000000E+02, -1.85000000E+02, -1.87000000E+02, -1.88000000E+02, -1.90000000E+02, -1.92000000E+02, -1.95000000E+02, -1.97000000E+02, -1.98000000E+02, -2.01000000E+02, -2.03000000E+02, -2.05000000E+02, -2.08000000E+02, -2.09000000E+02, -2.11000000E+02 /) - im_eps = (/ 3.76, 3.72, 3.55, 2.91, 1.39, 0.526, 0.271, 0.277, 0.274, 0.266, 0.259, 0.26, 0.262, 0.265, 0.27, 0.274, 0.278, 0.284, 0.296, 0.304, 0.315, 0.326, 0.339, 0.352, 0.368, 0.378, 0.398, 0.413, 0.429, 0.447, 0.462, 0.485, 0.507, 0.522, 0.545, 0.564, 0.585, 0.608, 0.633, 0.649, 0.674, 0.698, 0.717, 0.737, 0.773, 0.797, 0.814, 0.84, 0.866, 0.902, 0.937, 0.971, 0.991, 1.03, 1.08, 1.06, 1.09, 1.13, 1.14, 1.18, 1.22, 1.24, 1.29, 1.32, 1.36, 1.39, 1.43, 1.47, 1.51, 1.56, 1.61, 1.64, 1.7, 1.74, 1.79, 1.84, 1.89, 1.94, 2, 2.05, 2.09, 2.16, 2.22, 2.3, 2.33, 2.39, 2.46, 2.54, 2.61, 2.68, 2.74, 2.78, 2.84, 2.92, 3.03, 3.07, 3.15, 3.22, 3.28, 3.36, 3.4, 3.58, 3.64, 3.76, 3.78, 3.86, 3.93, 4, 4.2, 4.89, 5.04, 5.28, 5.3, 5.21, 5.28, 5.65, 5.9, 5.79, 6.05, 6.07, 6.07, 6.26, 6.47, 6.58, 6.71, 6.78, 6.9, 6.88, 7.07, 7.33, 7.43, 7.54, 7.63, 7.89, 7.8, 7.99, 8.19, 8.4, 8.69, 8.69, 8.96, 9.1, 8.96, 9.46, 9.77, 9.81, 10.3, 10.2, 10.2, 10.5, 10.6, 11.1, 10.9, 11.5, 11.4, 11.6, 11.9, 12.1, 12.1, 12.3, 12, 12.9, 13, 13.5, 13.4, 13.6, 13.9 /) + im_eps = (/ 3.76000000E+00, 3.72000000E+00, 3.55000000E+00, 2.91000000E+00, 1.39000000E+00, 5.26000000E-01, 2.71000000E-01, 2.77000000E-01, 2.74000000E-01, 2.66000000E-01, 2.59000000E-01, 2.60000000E-01, 2.62000000E-01, 2.65000000E-01, 2.70000000E-01, 2.74000000E-01, 2.78000000E-01, 2.84000000E-01, 2.96000000E-01, 3.04000000E-01, 3.15000000E-01, 3.26000000E-01, 3.39000000E-01, 3.52000000E-01, 3.68000000E-01, 3.78000000E-01, 3.98000000E-01, 4.13000000E-01, 4.29000000E-01, 4.47000000E-01, 4.62000000E-01, 4.85000000E-01, 5.07000000E-01, 5.22000000E-01, 5.45000000E-01, 5.64000000E-01, 5.85000000E-01, 6.08000000E-01, 6.33000000E-01, 6.49000000E-01, 6.74000000E-01, 6.98000000E-01, 7.17000000E-01, 7.37000000E-01, 7.73000000E-01, 7.97000000E-01, 8.14000000E-01, 8.40000000E-01, 8.66000000E-01, 9.02000000E-01, 9.37000000E-01, 9.71000000E-01, 9.91000000E-01, 1.03000000E+00, 1.08000000E+00, 1.06000000E+00, 1.09000000E+00, 1.13000000E+00, 1.14000000E+00, 1.18000000E+00, 1.22000000E+00, 1.24000000E+00, 1.29000000E+00, 1.32000000E+00, 1.36000000E+00, 1.39000000E+00, 1.43000000E+00, 1.47000000E+00, 1.51000000E+00, 1.56000000E+00, 1.61000000E+00, 1.64000000E+00, 1.70000000E+00, 1.74000000E+00, 1.79000000E+00, 1.84000000E+00, 1.89000000E+00, 1.94000000E+00, 2.00000000E+00, 2.05000000E+00, 2.09000000E+00, 2.16000000E+00, 2.22000000E+00, 2.30000000E+00, 2.33000000E+00, 2.39000000E+00, 2.46000000E+00, 2.54000000E+00, 2.61000000E+00, 2.68000000E+00, 2.74000000E+00, 2.78000000E+00, 2.84000000E+00, 2.92000000E+00, 3.03000000E+00, 3.07000000E+00, 3.15000000E+00, 3.22000000E+00, 3.28000000E+00, 3.36000000E+00, 3.40000000E+00, 3.58000000E+00, 3.64000000E+00, 3.76000000E+00, 3.78000000E+00, 3.86000000E+00, 3.93000000E+00, 4.00000000E+00, 4.20000000E+00, 4.89000000E+00, 5.04000000E+00, 5.28000000E+00, 5.30000000E+00, 5.21000000E+00, 5.28000000E+00, 5.65000000E+00, 5.90000000E+00, 5.79000000E+00, 6.05000000E+00, 6.07000000E+00, 6.07000000E+00, 6.26000000E+00, 6.47000000E+00, 6.58000000E+00, 6.71000000E+00, 6.78000000E+00, 6.90000000E+00, 6.88000000E+00, 7.07000000E+00, 7.33000000E+00, 7.43000000E+00, 7.54000000E+00, 7.63000000E+00, 7.89000000E+00, 7.80000000E+00, 7.99000000E+00, 8.19000000E+00, 8.40000000E+00, 8.69000000E+00, 8.69000000E+00, 8.96000000E+00, 9.10000000E+00, 8.96000000E+00, 9.46000000E+00, 9.77000000E+00, 9.81000000E+00, 1.03000000E+01, 1.02000000E+01, 1.02000000E+01, 1.05000000E+01, 1.06000000E+01, 1.11000000E+01, 1.09000000E+01, 1.15000000E+01, 1.14000000E+01, 1.16000000E+01, 1.19000000E+01, 1.21000000E+01, 1.21000000E+01, 1.23000000E+01, 1.20000000E+01, 1.29000000E+01, 1.30000000E+01, 1.35000000E+01, 1.34000000E+01, 1.36000000E+01, 1.39000000E+01 /) do i = 1, 167 eps_(i) = complex(re_eps(i), im_eps(i)) end do call interp1(lambda, eps_, wavelength, eps) -end function AgRaschke +end function epsAgRaschke function epsPt(wavelength) result(eps) !