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_get_kpts
       8              :       use m_juDFT
       9              :       USE m_types
      10              :       contains
      11           10 :       subroutine wann_get_kpts(input,kpts,
      12              :      >               l_bzsym,film,l_readkpts,
      13              :      <               nkpts,kpoints)
      14              : c********************************************************
      15              : c     Read in the k-points from kpts/w90kpts file.
      16              : c
      17              : c     Frank Freimuth
      18              : c********************************************************
      19              : 
      20              :       USE m_constants
      21              : 
      22              :       implicit none
      23              : 
      24              :       TYPE(t_input), INTENT(IN) :: input
      25              :       TYPE(t_kpts), INTENT(IN)  :: kpts
      26              :       logical,intent(in)  :: l_bzsym,film
      27              :       logical,intent(in)  :: l_readkpts
      28              :       integer,intent(out) :: nkpts
      29              :       real,intent(inout),allocatable  :: kpoints(:,:)
      30              : 
      31              :       real             :: scale
      32              :       integer          :: at,j
      33              :       integer          :: iter,len,num_wann,num_bands,nn,i
      34              :       logical          :: l_file
      35              : 
      36           10 :       call timestart("wann_get_kpts")
      37           10 :       if(l_bzsym)then
      38            0 :          inquire(file='w90kpts',exist=l_file)
      39            0 :          IF(.NOT.l_file) CALL juDFT_error("where is w90kpts?",calledby
      40            0 :      +        ="wann_get_kpts")
      41            0 :          open(987,file='w90kpts',status='old',form='formatted')
      42            0 :          read(987,*)nkpts, scale
      43            0 :          write(oUnit,*)"wann_get_kpts: nkpts=",nkpts
      44            0 :          if(l_readkpts)then
      45            0 :             IF(SIZE(kpoints,1)/=3) CALL juDFT_error("wann_get_kpts: 1"
      46            0 :      +           ,calledby ="wann_get_kpts")
      47            0 :             IF(SIZE(kpoints,2)/=nkpts)CALL juDFT_error("wann_get_kpts:2"
      48            0 :      +           ,calledby ="wann_get_kpts")
      49            0 :             do iter=1,nkpts
      50            0 :                read(987,*)kpoints(:,iter)
      51              :             enddo
      52              :          endif
      53            0 :          close(987)
      54              :       else
      55           10 :              nkpts = kpts%nkpt
      56           10 :             write(oUnit,*)"wann_get_kpts: nkpts=",nkpts
      57           10 :             if(l_readkpts)then
      58           45 :                do iter=1,nkpts
      59          165 :                   kpoints(:,iter) = kpts%bk(:,iter)
      60              :                enddo
      61              :             endif
      62              :       endif
      63              : 
      64           10 :       IF (l_readkpts) THEN
      65           45 :          do iter=1,nkpts
      66           45 :             write(oUnit,*)kpoints(:,iter)
      67              :          enddo
      68              :       END IF
      69              : 
      70           10 :       call timestop("wann_get_kpts")
      71           10 :       end subroutine wann_get_kpts
      72              :       end module m_wann_get_kpts