From 8a851176100dba627f9c387c3c3bf37e9ba04d60 Mon Sep 17 00:00:00 2001
From: Karl-Michael Schindler <karl-michael.schindler@web.de>
Date: Tue, 19 Nov 2024 17:16:46 +0100
Subject: [PATCH] Update seteps.f90
correct sign
---
source/f90/seteps.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/source/f90/seteps.f90 b/source/f90/seteps.f90
index 1d797da..464c95a 100644
--- a/source/f90/seteps.f90
+++ b/source/f90/seteps.f90
@@ -62,7 +62,7 @@ subroutine seteps(neps, nos, oscType, osc, epsinf, wn, eps, layers)
wp = osc(1,k)
y0 = osc(2,k)
yp = osc(3,k)
- addDrude = addDrude - dcmplx(wp**2, (y0 - yp) * wn) / dcmplx(wn2, wn * y0)
+ addDrude = addDrude - dcmplx(wp**2, (yp - y0) * wn) / dcmplx(wn2, wn * y0)
case default
write (*,'(A)') '*** Error in oscillator type! ***'
end select
--
GitLab