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 |