| 9 |
|
!! @author Matthew Meineke |
| 10 |
|
!! @author Christopher Fennel |
| 11 |
|
!! @author J. Daniel Gezelter |
| 12 |
< |
!! @version $Id: calc_sticky_pair.F90,v 1.2 2003-03-24 21:55:34 gezelter Exp $, $Date: 2003-03-24 21:55:34 $, $Name: not supported by cvs2svn $, $Revision: 1.2 $ |
| 12 |
> |
!! @version $Id: calc_sticky_pair.F90,v 1.3 2003-04-02 22:19:03 mmeineke Exp $, $Date: 2003-04-02 22:19:03 $, $Name: not supported by cvs2svn $, $Revision: 1.3 $ |
| 13 |
|
|
| 14 |
|
module sticky_pair |
| 15 |
|
|
| 240 |
|
#endif |
| 241 |
|
! Now, on to the forces: |
| 242 |
|
|
| 243 |
+ |
write(*,'(a4,3es12.3)') 't1= ', t(1:3,atom1) |
| 244 |
+ |
write(*,'(a4,3es12.3)') 't2= ', t(1:3,atom2) |
| 245 |
+ |
write(*,*) |
| 246 |
+ |
|
| 247 |
|
! first rotate the i terms back into the lab frame: |
| 248 |
|
|
| 249 |
|
#ifdef IS_MPI |
