diff --git a/src/eps.f90 b/src/eps.f90 index 4edba0b..6348a3a 100644 --- a/src/eps.f90 +++ b/src/eps.f90 @@ -59,7 +59,7 @@ function epsAu(wavelength) result(eps) end function epsAu ! - function AuRaschke(wavelength) result(eps) + function epsAuRaschke(wavelength) result(eps) !-------------------------------------------------------- ! Single-crystal (SC) dielectric function dataset for Au from Raschke et al PHYSICAL REVIEW B 86, 235147 (2012) !---------------------------------------------------------- @@ -70,11 +70,95 @@ function AuRaschke(wavelength) result(eps) real(8), dimension(164) :: lambda integer :: i - 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 = (/ -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.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 /) + 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 = (/ & + & -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.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)) @@ -107,7 +191,7 @@ function epsAg(wavelength) result(eps) end function epsAg ! - function AgRaschke(wavelength) result(eps) + function epsAgRaschke(wavelength) result(eps) !-------------------------------------------------------- ! Template-stripped dielectric function dataset (Cc) for Ag from Raschke et al Phys. Rev. B 91, 235137 (2015) !---------------------------------------------------------- @@ -118,11 +202,95 @@ function AgRaschke(wavelength) result(eps) real(8), dimension(167) :: lambda integer :: i - 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 = (/ 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.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 /) + 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 = (/ & + & 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.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)) diff --git a/src/termsProgram.f90 b/src/termsProgram.f90 index 05ebca6..7439d69 100644 --- a/src/termsProgram.f90 +++ b/src/termsProgram.f90 @@ -2513,14 +2513,12 @@ subroutine calcEpsilon() do j = 1, size(wavelen) escat(i, k, j) = epsAu(wavelen(j)) end do - elseif (trim(labels(i, k)) == 'AuRaschke') then - do j = 1, size(wavelen) - escat(i, k, j) = epsAuRaschke(wavelen(j)) - end do - elseif (trim(labels(i, k)) == 'AgRaschke') then - do j = 1, size(wavelen) - escat(i, k, j) = epsAgRaschke(wavelen(j)) - end do + elseif (trim(labels(i, k)) == 'Au2') then + escat(i, k, 1:size(wavelen)) = epsAuRaschke(wavelen) + + elseif (trim(labels(i, k)) == 'Ag2') then + escat(i, k, 1:size(wavelen)) = epsAgRaschke(wavelen) + elseif (trim(labels(i, k)) == 'Ag') then do j = 1, size(wavelen) escat(i, k, j) = epsAg(wavelen(j)) @@ -2544,6 +2542,7 @@ subroutine calcEpsilon() elseif (trim(labels(i, k)) == 'Cr') then escat(i, k, 1:size(wavelen)) = epsCr(wavelen) + elseif (trim(labels(i, k)) == 'Water') then escat(i, k, 1:size(wavelen)) = epsWater(wavelen) diff --git a/test/epsilon/input3 b/test/epsilon/input3 new file mode 100644 index 0000000..249210f --- /dev/null +++ b/test/epsilon/input3 @@ -0,0 +1,13 @@ +ModeAndScheme 2 3 +Wavelength 400 800 20 + +Medium 1.7689 + +OutputFormat HDF5 xsec + +Scatterers 2 +Ag2 0 0 0 10 +Au2 0 30 0 10 + + +