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

tmaxr deleted in hermitian_2

Browse Source
This commit is contained in:
Daniele Micheletti 2015-05-22 12:36:31 +00:00
parent fd86e1d821
commit 50ca63a2b9
1 changed files with 2 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
src/dispersion.f90
View File

@ -595,7 +595,7 @@ subroutine diel_tens_fr(xg,yg,mu,npl,e330,cr,ci,epsl,lrm,fast)
integer :: n,m,ih,i,k,n1,nn,llm,neval,ier,last,ihmin integer :: n,m,ih,i,k,n1,nn,llm,neval,ier,last,ihmin
integer, dimension(liw) :: iw integer, dimension(liw) :: iw
real(dp) :: mu2,mu4,mu6,npl2,dt,bth,bth2,bth4,bth6 real(dp) :: mu2,mu4,mu6,npl2,dt,bth,bth2,bth4,bth6
real(dp) :: t,x,s,sy1,sy2,sy3,resfh,epp,tmaxr real(dp) :: t,x,s,sy1,sy2,sy3,resfh,epp
real(dp), dimension(lw) :: w real(dp), dimension(lw) :: w
real(dp), dimension(npar) :: apar real(dp), dimension(npar) :: apar
! !
@ -618,7 +618,6 @@ subroutine diel_tens_fr(xg,yg,mu,npl,e330,cr,ci,epsl,lrm,fast)
n1=1 n1=1
if(fast.gt.10) n1=-llm if(fast.gt.10) n1=-llm
! !
tmaxr = real(tmax,dp)
apar(1) = yg apar(1) = yg
apar(2) = mu apar(2) = mu
apar(3) = npl apar(3) = npl
@ -634,7 +633,7 @@ subroutine diel_tens_fr(xg,yg,mu,npl,e330,cr,ci,epsl,lrm,fast)
do ih=ihmin,2 do ih=ihmin,2
apar(7) = real(ih,dp) apar(7) = real(ih,dp)
! call dqagi(fhermit,bound,2,epsa,epsr,resfh, ! call dqagi(fhermit,bound,2,epsa,epsr,resfh,
call dqagsmv(fhermit,-tmaxr,tmaxr,apar,npar,epsa,epsr,resfh, & call dqagsmv(fhermit,-tmax,tmax,apar,npar,epsa,epsr,resfh, &
epp,neval,ier,liw,lw,last,iw,w) epp,neval,ier,liw,lw,last,iw,w)
rr(n,ih,m) = resfh rr(n,ih,m) = resfh
end do end do