From 99513cfe7458ea75f179228f376f7d45c46a17bd Mon Sep 17 00:00:00 2001 From: kamischi <karl-michael.schindler@web.de> Date: Sat, 30 Dec 2023 21:35:04 +0100 Subject: [PATCH] More comments --- source/f90/calltree.txt | 4 +++- source/f90/doboson.f90 | 5 ----- source/f90/eels.f90 | 6 +++--- source/f90/o1.f90 | 4 ++++ source/f90/o2.f90 | 4 ++++ 5 files changed, 14 insertions(+), 9 deletions(-) diff --git a/source/f90/calltree.txt b/source/f90/calltree.txt index a248480..90bb8c9 100644 --- a/source/f90/calltree.txt +++ b/source/f90/calltree.txt @@ -22,4 +22,6 @@ BOSON sintr sintr rcffi - respon \ No newline at end of file + respon + o1 + o2 diff --git a/source/f90/doboson.f90 b/source/f90/doboson.f90 index 8ae4b5f..f2ee66d 100644 --- a/source/f90/doboson.f90 +++ b/source/f90/doboson.f90 @@ -57,11 +57,6 @@ subroutine doboson(t, width, gauss, asym, emin, emax, wmin, wmax, np, p, xout, y ! respon called when the input parameter gauss is < 0 or > 1. ! with 0 <= gauss <= 1, r1 and r2 are not used. -! *** rational approximations for ei(u) * exp(-u) + e1(u) * exp(+u) -! *** in the intervals (0, 1.3) and (1.3, infinity) <accuracy : 4.e-04> -! *** used for fourier transforming half-lorentzian functions - - debug = .false. data fm1 / 0.0d0 /, fp1 / 0.0d0 / pi = 4 * datan(1.0d0) diff --git a/source/f90/eels.f90 b/source/f90/eels.f90 index d6ac13f..ae45580 100644 --- a/source/f90/eels.f90 +++ b/source/f90/eels.f90 @@ -30,11 +30,11 @@ program eels call change_working_dir() open(unit = 11, file = 'eelsin') -! impact energy (ev) +! impact energy (eV) read(11, *) e0 -! incidence angle (%) +! incidence angle (°) read(11, *) theta -! angular apertures of the elliptic detector (%) +! angular apertures of the elliptic detector (°) read(11, *) phia read(11, *) phib ! energy-loss interval and step size (cm**-1) diff --git a/source/f90/o1.f90 b/source/f90/o1.f90 index a346fc0..16c9bf0 100644 --- a/source/f90/o1.f90 +++ b/source/f90/o1.f90 @@ -1,5 +1,9 @@ double precision function o1(u) +! *** rational approximations for ei(u) * exp(-u) + e1(u) * exp(+u) +! *** in the intervals (0, 1.3) and (1.3, infinity) <accuracy : 4.e-04> +! *** used for fourier transforming half-lorentzian functions + implicit none double precision, intent(in) :: u diff --git a/source/f90/o2.f90 b/source/f90/o2.f90 index 1d4fc3b..76c562b 100644 --- a/source/f90/o2.f90 +++ b/source/f90/o2.f90 @@ -1,5 +1,9 @@ double precision function o2(u) +! *** rational approximations for ei(u) * exp(-u) + e1(u) * exp(+u) +! *** in the intervals (0, 1.3) and (1.3, infinity) <accuracy : 4.e-04> +! *** used for fourier transforming half-lorentzian functions + implicit none double precision, intent(in) :: u -- GitLab