hdiff output

r33037/watermethane.f90 2017-07-19 16:30:08.285831077 +0100 r33036/watermethane.f90 2017-07-19 16:30:08.505833998 +0100
153: data deltaang8(5,:)/4.415080d0,4.415080d0,4.415080d0,4.415080d0,3.962102d0,3.741763d0,3.741763d0,3.741763d0,3.741763d0/153: data deltaang8(5,:)/4.415080d0,4.415080d0,4.415080d0,4.415080d0,3.962102d0,3.741763d0,3.741763d0,3.741763d0,3.741763d0/
154: !T1 and T2:154: !T1 and T2:
155: data deltaang8(6,:)/methanedof*0.0d0/155: data deltaang8(6,:)/methanedof*0.0d0/
156: data deltaang8(7,:)/methanedof*0.0d0/156: data deltaang8(7,:)/methanedof*0.0d0/
157: 157: 
158: contains158: contains
159: !-------------------------------------------------------------------------------------------159: !-------------------------------------------------------------------------------------------
160: subroutine clath_pot(x, grad, ereal, gradt, nwater)160: subroutine clath_pot(x, grad, ereal, gradt, nwater)
161: implicit none161: implicit none
162: 162: 
163: integer::              i,j,k, nwater 
164: double precision::     x(nwater*24+3*methanedof), grad(nwater*24+3*methanedof), ereal, ewater, ewatermeth163: double precision::     x(nwater*24+3*methanedof), grad(nwater*24+3*methanedof), ereal, ewater, ewatermeth
165: double precision::     xwatermeth(51), tempgrad(51), xwater(nwater*9), gradwater(nwater*9)164: double precision::     xwatermeth(51), tempgrad(51), xwater(nwater*9), gradwater(nwater*9)
 165: integer::              i,j,k, nwater
166: logical::             gradt166: logical::             gradt
167: 167: 
168: ereal=0.0d0168: ereal=0.0d0
169: grad(:)=0.0d0169: grad(:)=0.0d0
170: do i=1, nwater170: do i=1, nwater
171:    xwatermeth(:)=0.0d0171:    xwatermeth(:)=0.0d0
172:    do j=1, waterdof+1172:    do j=1, waterdof+1
173:       do k=1, 3173:       do k=1, 3
174:          xwatermeth(3*(j-1) + k)= x(3*((waterdof+1)*(i-1) + j -1) +k)/autoang174:          xwatermeth(3*(j-1) + k)= x(3*((waterdof+1)*(i-1) + j -1) +k)/autoang
175:       end do175:       end do
342: subroutine gser(gamser,a,x,gln)342: subroutine gser(gamser,a,x,gln)
343:   integer:: ITMAX343:   integer:: ITMAX
344:   double precision:: a,gamser,gln,x,EPS344:   double precision:: a,gamser,gln,x,EPS
345:   PARAMETER (ITMAX=100,EPS=3.e-7)345:   PARAMETER (ITMAX=100,EPS=3.e-7)
346:   ! Returns the incomplete gamma function P (a, x) evaluated by its series representation as346:   ! Returns the incomplete gamma function P (a, x) evaluated by its series representation as
347:   ! gamser. Also returns ln Γ(a) as gln.347:   ! gamser. Also returns ln Γ(a) as gln.
348:   integer:: n348:   integer:: n
349:   double precision:: ap,del,sum349:   double precision:: ap,del,sum
350:   gln=gammln(a)350:   gln=gammln(a)
351:   if(x.le.0.0d0) then351:   if(x.le.0.0d0) then
352:      if(x.lt.0.0d0) stop 'gser'352:      if(x.lt.0.0d0) stop('gser')
353:      gamser=0.0d0353:      gamser=0.0d0
354:      return354:      return
355:   endif355:   endif
356:   ap=a356:   ap=a
357:   sum=1.0d0/a357:   sum=1.0d0/a
358:   del=sum358:   del=sum
359:   do n=1,ITMAX359:   do n=1,ITMAX
360:      ap=ap+1.360:      ap=ap+1.
361:      del=del*x/ap361:      del=del*x/ap
362:      sum=sum+del362:      sum=sum+del
390:      b=b+2.390:      b=b+2.
391:      d=an*d+b391:      d=an*d+b
392:      if(abs(d).lt.fpmin) d=fpmin392:      if(abs(d).lt.fpmin) d=fpmin
393:      c=b+an/c393:      c=b+an/c
394:      if(abs(c).lt.fpmin) c=fpmin394:      if(abs(c).lt.fpmin) c=fpmin
395:      d=1.0d0/d395:      d=1.0d0/d
396:      del=d*c396:      del=d*c
397:      h=h*del397:      h=h*del
398:      if(abs(del-1.d0).lt.EPS) goto 1398:      if(abs(del-1.d0).lt.EPS) goto 1
399:   enddo399:   enddo
400:   stop 'gcf'400:   stop('gcf')
401: 1 gammcf=exp(-x+a*log(x)-gln)*h401: 1 gammcf=exp(-x+a*log(x)-gln)*h
402:   ! Put factors in front.402:   ! Put factors in front.
403:   return403:   return
404: end subroutine gcf404: end subroutine gcf
405: 405: 
406: function gammln(xx)406: function gammln(xx)
407:   double precision:: gammln,xx407:   double precision:: gammln,xx
408:   ! Returns the value ln[Gamma(xx)] for xx > 0.408:   ! Returns the value ln[Gamma(xx)] for xx > 0.
409:   integer j409:   integer j
410:   double precision:: ser,stp,tmp,x,y,cof(6)410:   double precision:: ser,stp,tmp,x,y,cof(6)


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0