diff --git a/source/f90/doeels.f90 b/source/f90/doeels.f90 index 559087080740e3ccb76578ec574b4959d0a4272c..2ca42d30e50e5692722b18167e066beb410eff13 100644 --- a/source/f90/doeels.f90 +++ b/source/f90/doeels.f90 @@ -8,6 +8,28 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, comment, comment_size, ! * statified medium made from isotropic materials in specular * ! * geometry using the dielectric theory of eels. * ! * * +! * e0: impact energy (eV) * +! * theta: incidence angle (°) * +! * phia, phib: angular apertures of the elliptic detector (°) * +! * wmin, wmax, dw: energy-loss interval and step size (cm**-1) * +! * comment: comment lines * +! * comment_size: number of comment lines * +! * layers: number of layers * +! * neps: number of epsilon values * +! * nper: number of periodic layers * +! * name: layer names * +! * name_size: number of layers * +! * thick: layer thickness * +! * epsinf: epsilon at infinite frequency * +! * nos: number of oscillators * +! * osc: oscillator parameters, wTO, Q, lambda * +! * osc_size: number of oscillators * +! * contrl: 'image' for image-charge screening * +! * mode: 'kurosawa' for kurosawa model * +! * wn_array: frequencies * +! * f_array: spectrum * +! * wn_array_size: number of frequencies * +! * * ! ****************************************************************** use quanc8_mod diff --git a/source/f90/fint1.f90 b/source/f90/fint1.f90 index 5fde9fa0ad21131650c64358d481e1bb36e08d9f..01b9dde14b199f827064aadde47debe1bb162af9 100644 --- a/source/f90/fint1.f90 +++ b/source/f90/fint1.f90 @@ -6,6 +6,12 @@ double precision function fint1(u, eps, thick, layers, nper, eps_size) ! * * ! * integration over the azimutal angle from 0.0 to pi * ! * * +! * eps: epsilon * +! * thick: thickness * +! * layers: number of layers * +! * nper: number of periodic layers * +! * eps_size: number of layers * +! * * ! ****************************************************************** use surlos_mod diff --git a/source/f90/fint2.f90 b/source/f90/fint2.f90 index 846c7d908fa72e46afd2d304564ef072df84d6ae..c86fd6997e7ae6c83d7090c05b3f93efaa29073b 100644 --- a/source/f90/fint2.f90 +++ b/source/f90/fint2.f90 @@ -6,6 +6,12 @@ double precision function fint2(u, eps, thick, layers, nper, eps_size) ! * * ! * integration over the azimutal angle from 0.0 to phi < pi * ! * * +! * eps: epsilon * +! * thick: thickness * +! * layers: number of layers * +! * nper: number of periodic layers * +! * eps_size: number of layers * +! * * ! ****************************************************************** use surlos_mod diff --git a/source/f90/fint3.f90 b/source/f90/fint3.f90 index 71dc07ee499f4e8a20c9b968ef43b0e7c6846971..796a54e94620450d31830f14bc164ce149cba598 100644 --- a/source/f90/fint3.f90 +++ b/source/f90/fint3.f90 @@ -6,6 +6,12 @@ double precision function fint3(u, eps, thick, layers, nper, eps_size) ! * * ! * integration over the azimutal angle from phi1 > 0 to phi2 < pi * ! * * +! * eps: epsilon * +! * thick: thickness * +! * layers: number of layers * +! * nper: number of periodic layers * +! * eps_size: number of layers * +! * * ! ****************************************************************** use surlos_mod diff --git a/source/f90/quanc8.f90 b/source/f90/quanc8.f90 index b4cfbf8c92e1e9761d386fc5aad6592a9e026db8..c34529b61d6b8baf00f81ab76cbd146c8d64b7ab 100644 --- a/source/f90/quanc8.f90 +++ b/source/f90/quanc8.f90 @@ -26,6 +26,13 @@ subroutine quanc8(fun, a, b, abserr, relerr, result, errest, nofun, flag, eps, t ! xxx.yyy , then xxx = the number of intervals which have ! not converged and 0.yyy = the fraction of the interval ! left to do when the limit on nofun was approached. +! +! additional input .. +! +! eps epsilon +! thick thickness +! layers number of layers +! nper number of periodic layers implicit none diff --git a/source/f90/queels.f90 b/source/f90/queels.f90 index 22ff1114e32620b1353fbc35b9e9ae8dfd565652..650000a838de7934d8246ee93f1c729eae65f5dc 100644 --- a/source/f90/queels.f90 +++ b/source/f90/queels.f90 @@ -12,6 +12,11 @@ subroutine queels(x, f, aerr, rerr, facru, eps, thick, layers, nper) ! * facru*x is the units of wavevectors omega/v_perpendicular * ! * f is the q-integral multiplied by (2/pi)**2 * ! * * +! * eps: epsilon * +! * thick: thickness * +! * layers: number of layers * +! * nper: number of periodic layers * +! * * ! ****************************************************************** use quanc8_mod diff --git a/source/f90/seteps.f90 b/source/f90/seteps.f90 index fdda8ae5a948452e519727e69cea7102682a8beb..e9604f4e343bb7583e23e4c7f5573e5d4e41cb95 100644 --- a/source/f90/seteps.f90 +++ b/source/f90/seteps.f90 @@ -7,6 +7,16 @@ subroutine seteps(neps, nos, osc, epsinf, wn, name, eps, layers, mode) ! * set up long-wavelength dielectric functions of the layers for * ! * the present frequency wn (in cm**-1) * ! * * +! * neps: number of epsilon values * +! * nos: number of oscillators * +! * osc: oscillator parameters, wTO, Q, lambda * +! * epsinf: epsilon at infinite frequency * +! * wn: frequencies * +! * name: layer names * +! * eps: epsilon * +! * layers: number of layers * +! * mode: 'kurosawa' for kurosawa model * +! * * ! ****************************************************************** implicit none diff --git a/source/f90/surlos.f90 b/source/f90/surlos.f90 index 59ec547b55219747eef12950ba8e7b76bc4c7ca9..71bd2e9e76dd43c03ada019c79311b35e3cb0c82 100644 --- a/source/f90/surlos.f90 +++ b/source/f90/surlos.f90 @@ -6,6 +6,12 @@ double precision function surlos(dk, eps, thick, layers, nper) ! * * ! * eels surface loss function for an arbitrary multilayered target* ! * * +! * dk: * +! * eps: epsilon * +! * thick: thickness * +! * layers: number of layers * +! * nper: number of periodic layers * +! * * ! ****************************************************************** implicit none