quantum-espresso/GWW/head/close_gwl.f90

78 lines
2.4 KiB
Fortran

!
! Copyright (C) 2001-2013 Quantum ESPRESSO group
! This file is distributed under the terms of the
! GNU General Public License. See the file `License'
! in the root directory of the present distribution,
! or http://www.gnu.org/copyleft/gpl.txt .
!
!
!
!----------------------------------------------------------------------------
SUBROUTINE close_phq( flag )
!----------------------------------------------------------------------------
!
! ... Close all files.
! ... Called at the end of the run with flag=.TRUE. (removes 'recover')
! ... or during execution with flag=.FALSE. (does not remove 'recover')
!
USE io_global, ONLY : ionode, stdout
USE uspp, ONLY : okvan
USE units_ph, ONLY : iubar, iudrhous, iuebar, iudrho, &
iudvscf, iucom, iudvkb3
USE units_lr, ONLY : iuwfc, iudwf
USE control_ph, ONLY : zue, epsil
USE recover_mod, ONLY : clean_recover
USE output, ONLY : fildrho, fildvscf
USE ramanm, ONLY : lraman, elop, iuchf, iud2w, iuba2
USE wannier_gw, ONLY : l_head
!
IMPLICIT NONE
!
LOGICAL :: flag
LOGICAL :: exst
!
!
CLOSE( UNIT = iuwfc, STATUS = 'DELETE' )
!
IF (flag) THEN
CLOSE( UNIT = iudwf, STATUS = 'DELETE' )
CLOSE( UNIT = iubar, STATUS = 'DELETE' )
!
IF ( okvan ) CLOSE( UNIT = iudrhous, STATUS = 'DELETE' )
!
IF ( epsil .OR. zue .OR. l_head) THEN
CLOSE( UNIT = iuebar, STATUS = 'DELETE' )
IF (okvan) CLOSE( UNIT = iucom, STATUS = 'DELETE' )
IF (okvan) CLOSE( UNIT = iudvkb3, STATUS = 'DELETE' )
ENDIF
ELSE
CLOSE( UNIT = iudwf, STATUS = 'KEEP' )
CLOSE( UNIT = iubar, STATUS = 'KEEP' )
!
IF ( okvan ) CLOSE( UNIT = iudrhous, STATUS = 'KEEP' )
!
IF ( epsil .OR. zue .OR. l_head) THEN
CLOSE( UNIT = iuebar, STATUS = 'KEEP' )
IF (okvan) CLOSE( UNIT = iucom, STATUS = 'KEEP' )
IF (okvan) CLOSE( UNIT = iudvkb3, STATUS = 'KEEP' )
ENDIF
ENDIF
!
IF ( ionode .AND. &
fildrho /= ' ') CLOSE( UNIT = iudrho, STATUS = 'KEEP' )
!
IF ( flag ) CALL clean_recover()
!
IF ( fildvscf /= ' ' ) CLOSE( UNIT = iudvscf, STATUS = 'KEEP' )
!
IF (lraman .OR.elop) THEN
CLOSE ( UNIT=iuchf, STATUS = 'keep' )
CLOSE ( UNIT=iud2w, STATUS = 'keep' )
CLOSE ( UNIT=iuba2, STATUS = 'keep' )
ENDIF
!
RETURN
!
END SUBROUTINE close_phq