Close lowdin file if wrote to it

This commit is contained in:
Sasha Fonari 2021-09-01 14:19:56 -04:00
parent db8b54e746
commit cfe899f0d7
1 changed files with 5 additions and 1 deletions

View File

@ -1394,12 +1394,16 @@ SUBROUTINE projwave( filproj, filowdin, lsym, lwrite_ovp )
!
IF (TRIM(filowdin) /= ' ') THEN
lowdin_unit = find_free_unit()
OPEN(unit=lowdin_unit, file=trim(filowdin), status='unknown', form='formatted')
OPEN(unit=lowdin_unit, file=trim(filowdin), status='unknown', &
form='formatted')
ELSE
lowdin_unit = stdout
END IF
!
CALL print_proj( lmax_wfc, proj, lowdin_unit )
!
IF (TRIM(filowdin) /= ' ') CLOSE( unit=lowdin_unit )
!
CALL write_proj_file ( filproj, proj )
!
END IF