istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

warmdisp bug fixed in trunk propagated to branches/refl

Browse Source
This commit is contained in:
Lorenzo Figini 2016-06-09 08:32:14 +00:00
parent 0aa1bb3537
commit 88705d8231
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/dispersion.f90
View File

@ -134,6 +134,8 @@ subroutine warmdisp(xg,yg,mu,npl,nprf,sox,lrm,err,nprr,npri,fast,imx,ex,ey,ez)
else else
call diel_tens_fr(yg,mu,npl,a330,epsl,lrm,fast) call diel_tens_fr(yg,mu,npl,a330,epsl,lrm,fast)
end if end if
a330= xg*a330
e330= one + a330
! !
do do
do i=1,imxx do i=1,imxx
@ -155,14 +157,12 @@ subroutine warmdisp(xg,yg,mu,npl,nprf,sox,lrm,err,nprr,npri,fast,imx,ex,ey,ez)
a33 = xg*sepsl(3,3) a33 = xg*sepsl(3,3)
a13 = xg*sepsl(1,3) a13 = xg*sepsl(1,3)
a23 = xg*sepsl(2,3) a23 = xg*sepsl(2,3)
a330= xg*a330
! a31 = a13 ! a31 = a13
! a32 =-a23 ! a32 =-a23
! !
e11 = one + a11 e11 = one + a11
e22 = one + a22 e22 = one + a22
e12 = a12 e12 = a12
e330= one + a330
! e33 = e330 + npra2*a33 ! e33 = e330 + npra2*a33
e13 = npra*a13 e13 = npra*a13
e23 = npra*a23 e23 = npra*a23