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             : 
       9             :       CONTAINS
      10             : 
      11           0 :       SUBROUTINE bs_comfort(eig_id,DIMENSION,input,noco,nkpt,param)
      12             : 
      13             :       USE m_types
      14             :       USE m_eig66_io, ONLY : read_eig
      15             : 
      16             :       IMPLICIT NONE
      17             : 
      18             :       TYPE(t_dimension), INTENT(IN)   :: DIMENSION
      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             : c---> pk non-collinear
      35           0 :       IF (noco%l_noco) THEN
      36             :          nspins = 1
      37             :       ELSE
      38           0 :          nspins = input%jspins
      39             :       ENDIF
      40             : c---> pk non-collinear
      41             : 
      42           0 :       ALLOCATE (eig(DIMENSION%neigd,nkpt,nspins))
      43           0 :       ALLOCATE (ne(nkpt,nspins))
      44             : 
      45           0 :       DO jsp = 1, nspins
      46           0 :          DO k = 1, nkpt
      47           0 :             CALL read_eig(eig_id,k,jsp,neig=ne(k,jsp),eig=eig(:,k,jsp))
      48             :          END DO ! k = 1,nkpt
      49             : 
      50           0 :          DO i = 1, DIMENSION%neigd
      51           0 :             DO k = 1, nkpt
      52           0 :                WRITE(776+jsp,*) param,k,eig(i,k,jsp)
      53             :             END DO
      54           0 :             WRITE(776+jsp,*)
      55             :          END DO
      56             : 
      57           0 :          DO i = 1, DIMENSION%neigd
      58           0 :             DO k = 1, nkpt
      59           0 :                IF (k.eq.param) WRITE(778+jsp,*) param,k,eig(i,k,jsp)
      60             :             END DO
      61           0 :             WRITE(778+jsp,*)
      62             :          END DO
      63             :       END DO ! jsp = 1,nspins
      64             : 
      65           0 :       DEALLOCATE (ne,eig)
      66             : 
      67           0 :       END SUBROUTINE bs_comfort
      68             :       END MODULE m_bs_comfort