quack/src/utils/read_basis_pyscf.f90

37 lines
786 B
Fortran

subroutine read_basis_pyscf(nBas,nO,nV)
! Read basis set information from PySCF
implicit none
include 'parameters.h'
! Input variables
integer,intent(in) :: nO(nspin)
! Local variables
! Output variables
integer,intent(out) :: nV(nspin)
integer,intent(out) :: nBas
!------------------------------------------------------------------------
! Primary basis set information
!------------------------------------------------------------------------
open(unit=3,file='int/nBas.dat')
read(3,*) nBas
close(unit=3)
write(*,'(A28)') '------------------'
write(*,'(A28,1X,I16)') 'Number of basis functions',nBas
write(*,'(A28)') '------------------'
write(*,*)
! Number of virtual orbitals
nV(:) = nBas - nO(:)
end subroutine