makekpts.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))& 
00010 =SKI(:,k)+(SKI(:,k+1)-SKI(:,k))/dble(Ndiv)*dble(i-1)  
00011 enddo 
00012 enddo  
00013 SK(:,Ndiv*(Nblk-1)+1)=SKI(:,Nblk) 
00014 !write(6,*) NSK_BAND_DISP
00015 !do i=1,Ndiv*(Nblk-1)+1
00016 !write(6,*) SK(:,i)
00017 !enddo 
00018 return 
00019 end