Skip to content
Snippets Groups Projects
Commit 399f5f98 authored by Karl-Michael Schindler's avatar Karl-Michael Schindler
Browse files

Correct sign in Kurosawa formula

parent 8a851176
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Igor Experiment/Test new seteps vs original.h5xp
+ 0
− 0
View file @ 399f5f98
No preview for this file type
source/f90/seteps.f90
+ 3
− 3
View file @ 399f5f98
......@@ -52,8 +52,8 @@ subroutine seteps(neps, nos, oscType, osc, epsinf, wn, eps, layers)
yTO = osc(2,k)
wLO = osc(3,k)
yLO = osc(4,k)
nomi = nomi * dcmplx(wLO**2 - wn2, yLO * wn)
deno = deno * dcmplx(wTO**2 - wn2, yTO * wn)
nomi = nomi * dcmplx(wLO**2 - wn2, -yLO * wn)
deno = deno * dcmplx(wTO**2 - wn2, -yTO * wn)
case (Drude)! Lambin 1990, equ. 8
wp = osc(1,k)
y = osc(2,k)
......@@ -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, (yp - y0) * wn) / dcmplx(wn2, wn * y0)
addDrude = addDrude - dcmplx(wp**2, (y0 - yp) * wn) / dcmplx(wn2, wn * y0)
case default
write (*,'(A)') '*** Error in oscillator type! ***'
end select
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment