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

Extended Drude works 

Change the sequence of parameters for consistency.
parent 25e5c7aa
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 @ adaa8a9a
No preview for this file type
source/f90/seteps.f90
+ 2
− 2
View file @ adaa8a9a
......@@ -60,8 +60,8 @@ subroutine seteps(neps, nos, oscType, osc, epsinf, wn, eps, layers)
addDrude = addDrude - wp**2 / dcmplx(wn2, y * wn)
case (extendedDrude) ! Gervais 2002, equ. 2.14
wp = osc(1,k)
yp = osc(2,k)
y0 = osc(3,k)
y0 = osc(2,k)
yp = osc(3,k)
addDrude = addDrude - dcmplx(wp**2, (y0 - yp) * wn) / dcmplx(wn2, wn * y0)
case default
write (*,'(A)') '*** Error in oscillator type! ***'
......
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