From 88705d823120fae1d7932ffc9b9c53e78d325141 Mon Sep 17 00:00:00 2001
From: Lorenzo Figini <lorenzo.figini@istp.cnr.it>
Date: Thu, 9 Jun 2016 08:32:14 +0000
Subject: [PATCH] warmdisp bug fixed in trunk propagated to branches/refl

---
 src/dispersion.f90 | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/src/dispersion.f90 b/src/dispersion.f90
index 9fbfec8..cc8deaa 100644
--- a/src/dispersion.f90
+++ b/src/dispersion.f90
@@ -134,6 +134,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
@@ -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)
       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