651 |
|
|
652 |
|
|
653 |
|
!! Does EAM pairwise Force calculation. |
654 |
< |
subroutine do_eam_pair(atom1,atom2,d,rij,r2,pot,f,do_pot,do_stress) |
654 |
> |
subroutine do_eam_pair(atom1, atom2, d, rij, r2, sw, vpair, pot, f, & |
655 |
> |
do_pot, do_stress) |
656 |
|
!Arguments |
657 |
|
integer, intent(in) :: atom1, atom2 |
658 |
|
real( kind = dp ), intent(in) :: rij, r2 |
659 |
< |
real( kind = dp ) :: pot |
659 |
> |
real( kind = dp ) :: pot, sw, vpair |
660 |
|
real( kind = dp ), dimension(3,nLocal) :: f |
661 |
|
real( kind = dp ), intent(in), dimension(3) :: d |
662 |
|
logical, intent(in) :: do_pot, do_stress |
789 |
|
pot_Row(atom1) = pot_Row(atom1) + phab*0.5 |
790 |
|
pot_Col(atom2) = pot_Col(atom2) + phab*0.5 |
791 |
|
end if |
792 |
+ |
vpair = vpair + phab |
793 |
|
|
794 |
|
f_Row(1,atom1) = f_Row(1,atom1) + fx |
795 |
|
f_Row(2,atom1) = f_Row(2,atom1) + fy |
803 |
|
if(do_pot) then |
804 |
|
pot = pot + phab |
805 |
|
end if |
806 |
+ |
vpair = vpair + phab |
807 |
|
|
808 |
|
f(1,atom1) = f(1,atom1) + fx |
809 |
|
f(2,atom1) = f(2,atom1) + fy |