← Back to branch summary

~siesta-ts/siesta/trunk_ts_soc

« back to all changes in this revision

Viewing changes to Util/Contour/grid1d.f

  • Committer: Nils Wittemeier
  • Date: 2019-02-14 07:45:07 UTC
  • mfrom: (746.1.15 trunk)
  • Revision ID: nils@4wittemeier.de-20190214074507-1mvzbmj9kw19gllr
Merged trunk 761

Show diffs side-by-side

added added

removed removed

Lines of Context:
Util/Contour/grid1d.f
1
1
! 
2
 
! Copyright (C) 1996-2016       The SIESTA group
 
2
! Copyright (C) 1996-2016       The SIESTA group
3
3
!  This file is distributed under the terms of the
4
4
!  GNU General Public License: see COPYING in the top directory
5
5
!  or http://www.gnu.org/copyleft/gpl.txt.
207
207
         do i = 1, 3
208
208
            celmod(n) = celmod(n) + cell(i,n)*cell(i,n)
209
209
         enddo
210
 
         sid(n) = celmod(n)/dfloat( mesh(n) )
 
210
         sid(n) = celmod(n)/real( mesh(n), kind=kind(1.d0) )
211
211
      enddo
212
212
 
213
213
 
214
214
c scan the segment in fractional coordinates of av --------------------------
215
215
 
216
216
      xincr = 0.d0
217
 
      if (nx .gt. 1) xincr = 1.d0 / dfloat(nx -1)
 
217
      if (nx .gt. 1) xincr = 1.d0 / real((nx -1),kind=kind(1.d0))
218
218
 
219
219
c loop in 1D grid -----------------------------------------------------------
220
220
 
252
252
 
253
253
c fractional coordinates of the point within the small 'cube' ---------------
254
254
 
255
 
            alp(n) =  rn(n)*mesh(n) - dfloat(m(n)-1) 
 
255
            alp(n) =  rn(n)*mesh(n) - real((m(n)-1),kind=kind(1.d0))
256
256
 
257
257
c bring the points to the unit cell (using the traslational symmetry) -------
258
258