Line data    Source code

       1             : c*************************c
       2             : c  routine to set up the  c
       3             : c  composite matrix anglmom c
       4             : c*************************c
       5             :       module m_wann_gwf_anglmom
       6             :       use m_juDFT
       7             :       implicit none
       8             :       contains
       9             : 
      10           0 :       subroutine wann_gwf_anglmom(nkpts,nqpts,l_unformatted)
      11             :       use m_wann_gwf_tools
      12             : 
      13             :       implicit none
      14             :       
      15             :       ! input parameters
      16             :       logical,intent(in) :: l_unformatted
      17             :       integer,intent(in) :: nkpts,nqpts
      18             : 
      19             :       integer :: nwfs,nbnd,nkqpts
      20             :       integer :: d,i,j,k,q,kq,t1,t2,t3,t4
      21             :       real :: tempr,tempi 
      22             : 
      23             :       complex,allocatable :: anglmom(:,:,:,:)
      24             :       character(len=12) :: fending
      25             : 
      26           0 :       call timestart("wann_gwf_anglmom")
      27           0 :       nkqpts = nkpts*nqpts
      28             : 
      29             :       ! get number of bands and wfs from proj
      30           0 :       open(405,file='proj',status='old')
      31           0 :       read(405,*)nwfs,nbnd
      32           0 :       close(405)
      33           0 :       write(*,*)'nbnd=',nbnd
      34           0 :       write(*,*)'nwfs=',nwfs
      35             : 
      36           0 :       allocate(anglmom(3,nbnd,nbnd,nkqpts))
      37             : 
      38             :       ! read in separate files
      39           0 :       do q=1,nqpts
      40           0 :        WRITE(fending,'("_",i4.4)')q
      41           0 :        open(405,file='WF1'//trim(fending)//'.anglmom')
      42           0 :        read(405,*) !header
      43           0 :        read(405,*) !nbnd and nkpts
      44           0 :        do k=1,nkpts
      45           0 :         kq=get_index_kq(k,q,nkpts)
      46           0 :         do i=1,nbnd
      47           0 :          do j=1,nbnd
      48           0 :           do d=1,3
      49           0 :            read(405,*)t1,t2,t3,t4,tempr,tempi
      50           0 :            anglmom(d,j,i,kq) = cmplx(tempr,tempi)
      51             :           enddo
      52             :          enddo
      53             :         enddo
      54             :        enddo
      55           0 :        close(405,status='delete')    
      56             :       enddo 
      57             : 
      58             :       ! write file either unformatted or formatted
      59           0 :       if(l_unformatted) then
      60           0 :        open(405,file='WF1_gwf.anglmom',form='unformatted')
      61           0 :        write(405)anglmom
      62           0 :        close(405)
      63             :       else
      64           0 :        open(405,file='WF1_gwf.anglmom')
      65           0 :        write(405,*)'Matrix elements of angular momentum'
      66           0 :        write(405,'(3i5)')nbnd,nbnd,nkqpts
      67           0 :        do kq=1,nkqpts
      68           0 :         do i=1,nbnd
      69           0 :          do j=1,nbnd
      70           0 :           do d=1,3
      71           0 :            write(405,'(4i5,3x,2f18.12)')d,j,i,kq,
      72           0 :      >                                  anglmom(d,j,i,kq)
      73             :           enddo
      74             :          enddo
      75             :         enddo
      76             :        enddo
      77           0 :        close(405)
      78             :       endif
      79             : 
      80           0 :       deallocate(anglmom)
      81             : 
      82           0 :       call timestop("wann_gwf_anglmom")
      83           0 :       end subroutine wann_gwf_anglmom
      84             :       end module m_wann_gwf_anglmom