From 871bd72fc19c8396b8c4b9f660cb76d786b45ede Mon Sep 17 00:00:00 2001
From: Karl-Michael Schindler <karl-michael.schindler@web.de>
Date: Wed, 11 Dec 2024 15:41:42 +0100
Subject: [PATCH] Update doeels.f90
revert accidental change
---
source/f90/doeels.f90 | 13 ++++++++++++-
1 file changed, 12 insertions(+), 1 deletion(-)
diff --git a/source/f90/doeels.f90 b/source/f90/doeels.f90
index 829415c..677c81b 100644
--- a/source/f90/doeels.f90
+++ b/source/f90/doeels.f90
@@ -204,7 +204,18 @@ subroutine doeels (e0, theta, phia, phib, wmin, wmax, dw, comment, comment_size,
wn = wmin + (iw - 1) * dw
! if (debug) write (*,*) 'wn: ', wn
if (wn >= 0.0d0) then
-q
+ if (wn /= 0.0d0) then
+ if (.not. user) call seteps(neps, nos, oscType, osc, epsinf, wn, eps, layers)
+
+ x = wn / (2 * ener * psia)
+ if (rational) then
+ f = qrat(x, alpha, beta, c1, c2) * dimag(-2 / (1.0d0 + eps(1)))
+ else
+ call queels(x, f, aerr, rerr, facru, eps, thick, layers, nper)
+ endif
+ f = prefac * f / wn
+ endif ! wn /= 0.0d0
+
wn_array(iw) = wn
f_array(iw) = f
--
GitLab