Skip to content

Commit

Permalink
aligned 5.0.2 patch to the newer (mass correction)
Browse files Browse the repository at this point in the history
  • Loading branch information
Iximiel committed May 28, 2024
1 parent 99e9255 commit 0b587d4
Showing 1 changed file with 9 additions and 4 deletions.
13 changes: 9 additions & 4 deletions patches/qespresso-5.0.2.diff/PW/src/plugin_forces.f90
Original file line number Diff line number Diff line change
Expand Up @@ -19,19 +19,21 @@ SUBROUTINE plugin_forces()
USE plugin_flags
!
USE cell_base, ONLY : alat, at
USE ions_base, ONLY : tau, nat,amass
USE ions_base, ONLY : tau, nat, amass, ityp
USE force_mod, ONLY : force,sigma
USE control_flags, ONLY : istep
USE ener, ONLY : etot
!
IMPLICIT NONE
!
INTEGER:: i,j
INTEGER:: i,j,ia
REAL(DP) :: at_plumed(3,3)
REAL(DP) :: virial(3,3)
REAL(DP) :: volume
REAL(DP), ALLOCATABLE :: tau_plumed(:,:)
REAL(DP) :: masses_plumed(nat)
!
masses_plumed = 0.0_DP
IF(use_plumed) then
IF(ionode)THEN
at_plumed=alat*at; ! the cell, rescaled properly
Expand All @@ -44,9 +46,12 @@ SUBROUTINE plugin_forces()
-at_plumed(1,2)*at_plumed(3,3)*at_plumed(2,1) &
-at_plumed(1,3)*at_plumed(3,1)*at_plumed(2,2)
virial=-sigma*volume

! the masses in QE are stored per type, see q-e//Modules/ions_base.f90
do ia=1,nat
masses_plumed(ia)=amass(ityp(ia))
end do
CALL plumed_f_gcmd("setStep"//char(0),istep)
CALL plumed_f_gcmd("setMasses"//char(0),amass)
CALL plumed_f_gcmd("setMasses"//char(0),masses_plumed)
CALL plumed_f_gcmd("setForces"//char(0),force)
CALL plumed_f_gcmd("setPositions"//char(0),tau_plumed)
CALL plumed_f_gcmd("setBox"//char(0),at_plumed)
Expand Down

0 comments on commit 0b587d4

Please sign in to comment.