Commit d83e3039 authored by giannozz's avatar giannozz

My change to the newly introduced "force theorem" computation of magnetic

anisotropy broke projwfc.x in that case. Here a hopefully better fix. It 
should be sufficient not to read 'config' from file in that specific case.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/[email protected] c92efa57-630b-4861-b058-cf58834340f0
parent 0f2a35bc
......@@ -1373,16 +1373,16 @@ SUBROUTINE iosys()
wfc_dir = tmp_dir
ENDIF
!
! IF ( lmovecell ) THEN
at_old = at
omega_old = omega
! ENDIF
!
! ... Read atomic positions and unit cell from data file, if needed,
! ... overwriting what has just been read before from input
!
ierr = 1
IF ( startingconfig == 'file' ) ierr = read_config_from_file(nat, at_old,omega_old, lmovecell, at, bg, omega, tau)
IF ( startingconfig == 'file' .AND. .NOT. lforcet ) &
ierr = read_config_from_file(nat, at_old, omega_old, lmovecell, &
at, bg, omega, tau)
!
! ... read_config_from_file returns 0 if structure successfully read
! ... Atomic positions (tau) must be converted to internal units
......
......@@ -971,8 +971,6 @@ MODULE pw_restart
USE scf, ONLY : rho
USE lsda_mod, ONLY : nspin
USE mp_bands, ONLY : intra_bgrp_comm
USE spin_orb, ONLY : lforcet
USE control_flags, ONLY : lscf, lbands
USE mp, ONLY : mp_sum
!
IMPLICIT NONE
......@@ -1139,21 +1137,6 @@ MODULE pw_restart
!
END SELECT
!
!-- To do Force Theorem calculation (AlexS)
!
IF ( lforcet.and. .NOT.lscf .AND. .NOT.lbands ) THEN
IF ( what.eq.'config' ) THEN
ierr = 1
ELSE
IF (lrho) CALL read_rho( rho, nspin )
ENDIF
RETURN
ENDIF
!--
IF ( .NOT. lheader .AND. .NOT. qexml_version_init) &
CALL errore( 'pw_readfile', 'qexml version not set', 71 )
!
......
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