--- trunk/OOPSE/libmdtools/do_Forces.F90	2003/07/16 21:30:56	626
+++ trunk/OOPSE/libmdtools/do_Forces.F90	2003/07/17 19:25:51	631
@@ -4,7 +4,7 @@
 
 !! @author Charles F. Vardeman II
 !! @author Matthew Meineke
-!! @version $Id: do_Forces.F90,v 1.21 2003-07-16 21:30:55 mmeineke Exp $, $Date: 2003-07-16 21:30:55 $, $Name: not supported by cvs2svn $, $Revision: 1.21 $
+!! @version $Id: do_Forces.F90,v 1.22 2003-07-17 19:25:51 chuckv Exp $, $Date: 2003-07-17 19:25:51 $, $Name: not supported by cvs2svn $, $Revision: 1.22 $
 
 module do_Forces
   use force_globals
@@ -603,6 +603,51 @@ contains
   end subroutine do_pair
 
 
+
+  subroutine do_preforce(i, j, rijsq, d, do_pot, do_stress, u_l, A, f, t, pot)
+   real( kind = dp ) :: pot
+   real( kind = dp ), dimension(3,getNlocal()) :: u_l
+   real (kind=dp), dimension(9,getNlocal()) :: A
+   real (kind=dp), dimension(3,getNlocal()) :: f
+   real (kind=dp), dimension(3,getNlocal()) :: t
+   
+   logical, intent(inout) :: do_pot, do_stress
+   integer, intent(in) :: i, j
+   real ( kind = dp ), intent(inout)    :: rijsq
+   real ( kind = dp )                :: r
+   real ( kind = dp ), intent(inout) :: d(3)
+   
+   logical :: is_EAM_i, is_EAM_j
+   
+   integer :: me_i, me_j
+   
+   r = sqrt(rijsq)
+   
+#ifdef IS_MPI
+   if (tagRow(i) .eq. tagColumn(j)) then
+      write(0,*) 'do_pair is doing', i , j, tagRow(i), tagColumn(j)
+   endif
+   
+   me_i = atid_row(i)
+   me_j = atid_col(j)
+   
+#else
+   
+   me_i = atid(i)
+   me_j = atid(j)
+   
+#endif
+   
+   if (FF_uses_EAM .and. SimUsesEAM()) then
+      call getElementProperty(atypes, me_i, "is_EAM", is_EAM_i)
+      call getElementProperty(atypes, me_j, "is_EAM", is_EAM_j)
+      
+      if ( is_EAM_i .and. is_EAM_j ) &
+           call calc_EAM_prepair(i, j, d, r, rijsq )
+   endif
+  end subroutine do_preforce
+  
+  
   subroutine get_interatomic_vector(q_i, q_j, d, r_sq)
     
     real (kind = dp), dimension(3) :: q_i