--- trunk/OOPSE/libmdtools/calc_sticky_pair.F90 2003/03/24 21:55:34 394 +++ trunk/OOPSE/libmdtools/calc_sticky_pair.F90 2003/04/02 22:19:03 443 @@ -9,7 +9,7 @@ !! @author Matthew Meineke !! @author Christopher Fennel !! @author J. Daniel Gezelter -!! @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 $ +!! @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 $ module sticky_pair @@ -240,6 +240,10 @@ contains #endif ! Now, on to the forces: + write(*,'(a4,3es12.3)') 't1= ', t(1:3,atom1) + write(*,'(a4,3es12.3)') 't2= ', t(1:3,atom2) + write(*,*) + ! first rotate the i terms back into the lab frame: #ifdef IS_MPI