Commit ea8fd320 authored by Frank Freimuth's avatar Frank Freimuth

Update wann_uHu.F

parent b7a59e7d
......@@ -1030,36 +1030,36 @@ c if (mod(i_rec-1,isize).eq.irank) then
& nbands,eigg,zzMat)
do i = 1,nbands
if((eigg(i).ge.e1s.and.nslibd.lt.numbands
& .and.wann%l_bynumber)
&.or.(eigg(i).ge.e1s.and.eigg(i).le.e2s.and.wann%l_byenergy)
&.or.(i.ge.wann%band_min(jspin)
& .and.
& (i.le.wann%band_max(jspin))
& .and.
& wann%l_byindex))then
nslibd = nslibd + 1
eig(nslibd) = eigg(i)
we(nslibd) = we(i)
if(l_noco)then
! funbas = nv(1) + nlotot
! funbas = funbas+nv(2) + nlotot
else
! funbas = nv(jspin) + nlotot
endif
IF (zzMat%l_real) THEN
! do j = 1,funbas
zMat%data_r(:,nslibd) = zzMat%data_r(:,i)
! enddo
ELSE
! do j = 1,funbas
zMat%data_c(:,nslibd) = zzMat%data_c(:,i)
! enddo
END IF
endif
enddo
c$$$ nslibd=0
! do i = 1,nbands
! if((eigg(i).ge.e1s.and.nslibd.lt.numbands
! & .and.wann%l_bynumber)
! &.or.(eigg(i).ge.e1s.and.eigg(i).le.e2s.and.wann%l_byenergy)
! &.or.(i.ge.wann%band_min(jspin)
! & .and.
! & (i.le.wann%band_max(jspin))
! & .and.
! & wann%l_byindex))then
! nslibd = nslibd + 1
! eig(nslibd) = eigg(i)
! we(nslibd) = we(i)
! if(l_noco)then
! ! funbas = nv(1) + nlotot
! ! funbas = funbas+nv(2) + nlotot
! else
! ! funbas = nv(jspin) + nlotot
! endif
! IF (zzMat%l_real) THEN
! ! do j = 1,funbas
! zMat%data_r(:,nslibd) = zzMat%data_r(:,i)
! ! enddo
! ELSE
! ! do j = 1,funbas
! zMat%data_c(:,nslibd) = zzMat%data_c(:,i)
! ! enddo
! END IF
! endif
! enddo
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment