Line data    Source code

       1              : !--------------------------------------------------------------------------------
       2              : ! Copyright (c) 2016 Peter Grünberg Institut, Forschungszentrum Jülich, Germany
       3              : ! This file is part of FLEUR and available as free software under the conditions
       4              : ! of the MIT license as expressed in the LICENSE file in more detail.
       5              : !--------------------------------------------------------------------------------
       6              : 
       7              :       MODULE m_wann_plot_symm
       8              :          use m_juDFT
       9              :       CONTAINS
      10            0 :       SUBROUTINE wann_plot_symm(jspin,mrot,ikpt,kplot,l_conjugate)
      11              : 
      12              :       implicit none
      13              :       integer,intent(in)::jspin,ikpt,kplot
      14              :       logical,intent(in)::l_conjugate
      15              :       integer,intent(in)::mrot(3,3)
      16              :       character(len=10)::vandernameold
      17              :       character(len=10)::vandernamenew
      18              :       integer::grid(3),nslibd,ikpt_copy,nbnd
      19              :       integer::ix,iy,iz,ixx,iyy,izz
      20              :       complex::xdnout
      21            0 :       real,allocatable::xdnreal(:,:,:,:),xdnimag(:,:,:,:)
      22              :       real::point(3)
      23              :       real::rotpoint(3)
      24              :       integer::opoint1,opoint2,opoint3
      25              : 
      26            0 :       call timestart("wann_plot_symm")
      27            0 :       WRITE (vandernameold,201) ikpt,jspin
      28            0 :       write (vandernamenew,201) kplot,jspin
      29              :  201  FORMAT ('UNK',i5.5,'.',i1)      
      30            0 :       open(55,file=vandernameold)
      31            0 :       read(55,*)grid(1),grid(2),grid(3),ikpt_copy,nslibd
      32            0 :       allocate(xdnreal(0:grid(1)-1,0:grid(2)-1,0:grid(3)-1,nslibd))
      33            0 :       allocate(xdnimag(0:grid(1)-1,0:grid(2)-1,0:grid(3)-1,nslibd))
      34            0 :       do nbnd=1,nslibd
      35            0 :        do iz=0,grid(3)-1
      36            0 :         do iy=0,grid(2)-1
      37            0 :          do ix=0,grid(1)-1
      38            0 :             read(55,*)xdnreal(ix,iy,iz,nbnd),xdnimag(ix,iy,iz,nbnd)
      39              :          enddo
      40              :         enddo
      41              :        enddo
      42              :       enddo
      43            0 :       close(55)
      44            0 :       open(666,file=vandernamenew,form='formatted')
      45            0 :       write(666,'(5i4)')grid(1),grid(2),grid(3),kplot,nslibd
      46            0 :       do nbnd=1,nslibd
      47            0 :        do iz=0,grid(3)-1
      48            0 :         do iy=0,grid(2)-1
      49            0 :          do ix=0,grid(1)-1
      50            0 :           point(1)=real(ix)/grid(1)
      51            0 :           point(2)=real(iy)/grid(2)
      52            0 :           point(3)=real(iz)/grid(3)
      53            0 :           rotpoint(:)=matmul(point(1:3),1.0*mrot(1:3,1:3))
      54            0 :           where(rotpoint.lt.0)
      55              :              rotpoint=rotpoint+1.0
      56              :           endwhere
      57            0 :           opoint1=int(grid(1)*rotpoint(1)+0.01)
      58            0 :           opoint2=int(grid(2)*rotpoint(2)+0.01)
      59            0 :           opoint3=int(grid(3)*rotpoint(3)+0.01)
      60            0 :           if(l_conjugate)then
      61            0 :             write(666,*)xdnreal(opoint1,opoint2,opoint3,nbnd),
      62            0 :      &               -xdnimag(opoint1,opoint2,opoint3,nbnd)
      63              :           else
      64            0 :             write(666,*)xdnreal(opoint1,opoint2,opoint3,nbnd),
      65            0 :      &               xdnimag(opoint1,opoint2,opoint3,nbnd)
      66              :           endif
      67              :          enddo
      68              :         enddo
      69              :        enddo
      70              :       enddo
      71            0 :       close(666)
      72              : 
      73            0 :       call timestop("wann_plot_symm")
      74            0 :       END SUBROUTINE wann_plot_symm
      75              :       END MODULE m_wann_plot_symm