Line data    Source code

       1              : !--------------------------------------------------------------------------------
       2              : ! Copyright (c) 2017 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_bs_comfort
       8              :          use m_juDFT
       9              :       CONTAINS
      10              : 
      11            0 :       SUBROUTINE bs_comfort(eig_id,input,noco,nkpt,param)
      12              : 
      13              :       USE m_types
      14              :       USE m_eig66_io, ONLY : read_eig
      15              : 
      16              :       IMPLICIT NONE
      17              : 
      18              :       
      19              :       TYPE(t_input),     INTENT(IN)   :: input
      20              :       TYPE(t_noco),      INTENT(IN)   :: noco
      21              : 
      22              : C     .. Scalar Arguments ..
      23              :       INTEGER, INTENT (IN) :: eig_id
      24              :       INTEGER, INTENT (IN) :: nkpt
      25              :       INTEGER, INTENT (IN) :: param
      26              : 
      27              : C     .. Local Scalars ..
      28              :       INTEGER              :: i, jsp, k, nspins
      29              : 
      30              : C     .. Local Arrays ..
      31            0 :       REAL,    ALLOCATABLE :: eig(:,:,:)
      32            0 :       INTEGER, ALLOCATABLE :: ne(:,:)
      33              : 
      34            0 :       call timestart("bs_comfort")
      35              : c---> pk non-collinear
      36            0 :       IF (noco%l_noco) THEN
      37              :          nspins = 1
      38              :       ELSE
      39            0 :          nspins = input%jspins
      40              :       ENDIF
      41              : c---> pk non-collinear
      42              : 
      43            0 :       ALLOCATE (eig(input%neig,nkpt,nspins))
      44            0 :       ALLOCATE (ne(nkpt,nspins))
      45              : 
      46            0 :       DO jsp = 1, nspins
      47            0 :          DO k = 1, nkpt
      48            0 :             CALL read_eig(eig_id,k,jsp,neig=ne(k,jsp),eig=eig(:,k,jsp))
      49              :          END DO ! k = 1,nkpt
      50              : 
      51            0 :          DO i = 1, input%neig
      52            0 :             DO k = 1, nkpt
      53            0 :                WRITE(776+jsp,*) param,k,eig(i,k,jsp)
      54              :             END DO
      55            0 :             WRITE(776+jsp,*)
      56              :          END DO
      57              : 
      58            0 :          DO i = 1, input%neig
      59            0 :             DO k = 1, nkpt
      60            0 :                IF (k.eq.param) WRITE(778+jsp,*) param,k,eig(i,k,jsp)
      61              :             END DO
      62            0 :             WRITE(778+jsp,*)
      63              :          END DO
      64              :       END DO ! jsp = 1,nspins
      65              : 
      66            0 :       DEALLOCATE (ne,eig)
      67              : 
      68            0 :       call timestop("bs_comfort")
      69              : 
      70            0 :       END SUBROUTINE bs_comfort
      71              :       END MODULE m_bs_comfort