Commit b7c8069e authored by Daniel Wortmann's avatar Daniel Wortmann

Put read of old inp into inpgen2

parent 3631692a
......@@ -79,8 +79,8 @@ PROGRAM inpgen
!CALL read_old_input()
l_fullinput=.TRUE.
ELSEIF (judft_was_argument("-inp.xml")) THEN
STOP "not yet"
!CALL r_inpXML()
call read_old_inp(input,DIMENSION,atoms,sphhar,cell,stars,sym,noco,vacuum,forcetheo,&
sliceplot,banddos,enpara,xcpot,kpts,hybrid, oneD)
l_fullinput=.TRUE.
ELSEIF(judft_was_argument("-f")) THEN
!read the input
......
......@@ -34,7 +34,6 @@
USE m_inpnoco
USE m_constants
USE m_types
USE m_inv3
USE m_setlomap
IMPLICIT NONE
! ..
......@@ -169,10 +168,11 @@
cell%amat(:,1) = a1(:)
cell%amat(:,2) = a2(:)
cell%amat(:,3) = a3(:)
CALL inv3(cell%amat,cell%bmat,cell%omtil)
cell%bmat(:,:) = tpi_const*cell%bmat(:,:)
cell%bbmat=MATMUL(cell%bmat,TRANSPOSE(cell%bmat))
cell%omtil = ABS(cell%omtil)
call cell%init(-1.)
!CALL inv3(cell%amat,cell%bmat,cell%omtil)
!cell%bmat(:,:) = tpi_const*cell%bmat(:,:)
!cell%bbmat=MATMUL(cell%bmat,TRANSPOSE(cell%bmat))
!cell%omtil = ABS(cell%omtil)
IF (input%film .AND. .NOT.oneD%odd%d1) THEN
cell%vol = cell%omtil/cell%amat(3,3)*vacuum%dvac
......
TYPE(t_input) :: input
TYPE(t_atoms) :: atoms
TYPE(t_cell) :: cell
TYPE(t_sym) :: sym
TYPE(t_noco) :: noco
TYPE(t_vacuum) :: vacuum
TYPE(t_banddos) :: banddos
TYPE(t_hybrid) :: hybrid
TYPE(t_xcpot_inbuild_nf)::xcpot
TYPE(t_oned) :: oned
TYPE(t_sliceplot):: sliceplot
TYPE(t_stars) :: stars
MODULE m_read_old_inp
IMPLICIT NONE
CONTAINS
SUBROUTINE read_old_inp(input,DIMENSION,atoms,sphhar,cell,stars,sym,noco,vacuum,forcetheo,&
sliceplot,banddos,enpara,xcpot,kpts,hybrid,&
oneD)
USE m_types_input
USE m_types_atoms
USE m_types_cell
USE m_types_sym
USE m_types_noco
USE m_types_vacuum
USE m_types_banddos
USE m_types_hybrid
USE m_types_xcpot_inbuild_nofunction
USE m_types_forcetheo
USE m_types_kpts
USE m_types_enpara
USE m_types_oneD
USE m_types_sliceplot
USE m_types_stars
USE m_fleur_init_old
TYPE(t_enpara) :: enpara
TYPE(t_input),INTENT(OUT):: input
TYPE(t_atoms),INTENT(OUT):: atoms
TYPE(t_cell) ,INTENT(OUT):: cell
TYPE(t_sym) ,INTENT(OUT):: sym
TYPE(t_noco) ,INTENT(OUT):: noco
TYPE(t_vacuum) ,INTENT(OUT):: vacuum
TYPE(t_banddos),INTENT(OUT):: banddos
TYPE(t_hybrid) ,INTENT(OUT):: hybrid
TYPE(t_xcpot_inbuild_nf),INTENT(OUT)::xcpot
TYPE(t_oned) ,INTENT(OUT):: oned
TYPE(t_sliceplot),INTENT(OUT):: sliceplot
TYPE(t_stars),INTENT(OUT):: stars
TYPE(t_enpara) ,INTENT(OUT):: enpara
TYPE(t_forcetheo):: forcetheo
TYPE(t_kpts) :: kpts
TYPE(t_kpts) ,INTENT(OUT):: kpts
TYPE(t_coreSpecInput)::coreSpecInput
call fleur_init_old(&
CALL fleur_init_old(&
input,DIMENSION,atoms,sphhar,cell,stars,sym,noco,vacuum,forcetheo,&
sliceplot,banddos,enpara,xcpot,kpts,hybrid,&
oneD,coreSpecInput)
END SUBROUTINE read_old_inp
END MODULE m_read_old_inp
......@@ -33,7 +33,6 @@ CONTAINS
USE m_calculator
USE m_constants
USE m_inpeig
USE m_inpnoco
USE m_sort
USE m_types_xcpot_inbuild
#ifdef CPP_LIBXC
......@@ -2060,12 +2059,6 @@ input%preconditioning_param = evaluateFirstOnly(xmlGetAttributeValue('/fleurInpu
CALL enpara%READ(atoms,input%jspins,input%film,.FALSE.)
END IF
! Read in nocoinp file iff available
l_nocoinp = .FALSE.
INQUIRE (file ='nocoinp',exist= l_nocoinp)
IF (l_nocoinp) THEN
CALL inpnoco(atoms,input,vacuum,noco)
END IF
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment