fixed e330 calculation in dispersion module (trunk)

2016-06-09 07:41:23 +00:00 · 2016-06-09 07:41:23 +00:00 · feaa736c0e
commit feaa736c0e
parent bede98b0ae
1 changed files with 2 additions and 2 deletions
--- a/src/dispersion.f90
+++ b/src/dispersion.f90
@ -129,6 +129,8 @@ subroutine warmdisp(xg,yg,mu,npl,nprf,sox,lrm,err,nprr,npri,fast,imx,ex,ey,ez)
  else
    call diel_tens_fr(yg,mu,npl,a330,epsl,lrm,fast)
  end if
+  a330= xg*a330
+  e330= one + a330
 !
  do
    do i=1,imxx
@ -150,14 +152,12 @@ subroutine warmdisp(xg,yg,mu,npl,nprf,sox,lrm,err,nprr,npri,fast,imx,ex,ey,ez)
      a33 = xg*sepsl(3,3)
      a13 = xg*sepsl(1,3)
      a23 = xg*sepsl(2,3)
-      a330= xg*a330
 !      a31 = a13
 !      a32 =-a23
 !
      e11 = one + a11
      e22 = one + a22
      e12 = a12
-      e330= one + a330
 !      e33 = e330 + npra2*a33
      e13 = npra*a13
      e23 = npra*a23