make_kpts.f90

Go to the documentation of this file.
00001 subroutine makekpts(Ndiv,Nblk,NSK_BAND_DISP,SKI,SK)  
00002 implicit none 
00003 integer::Nblk,Ndiv,NSK_BAND_DISP 
00004 real(8)::SKI(3,Nblk) 
00005 real(8)::SK(3,NSK_BAND_DISP)
00006 integer::i,k  
00007 do k=1,Nblk-1 
00008  do i=1,Ndiv 
00009   SK(:,i+Ndiv*(k-1))=SKI(:,k)+(SKI(:,k+1)-SKI(:,k))/dble(Ndiv)*dble(i-1)  
00010  enddo!i 
00011 enddo!k   
00012 SK(:,Ndiv*(Nblk-1)+1)=SKI(:,Nblk) 
00013 !write(6,*) NSK_BAND_DISP
00014 !do i=1,Ndiv*(Nblk-1)+1
00015 !write(6,*) SK(:,i)
00016 !enddo 
00017 return 
00018 end