--- trunk/OOPSE/libmdtools/notifyCutoffs.F90	2003/07/16 21:30:56	626
+++ trunk/OOPSE/libmdtools/notifyCutoffs.F90	2004/05/11 21:31:15	1160
@@ -2,36 +2,51 @@ module notifyCutoffs
   
   use definitions
   use do_Forces, only:      setRlistDF
-  use dipole_dipole, only:  setCutoffsDipole
   use reaction_field, only: setCutoffsRF
   use lj, only:             setCutoffLJ
-
+  use eam, only:            setCutoffEAM
+  use switcheroo, only:     set_switch
   implicit none
 
   PRIVATE
+
+#define __FORTRAN90
+#include "fSwitchingFunction.h"
   
   public::cutoffNotify
 
   contains
     
-    subroutine cutoffNotify( this_rcut, this_rlist, this_ecr, this_est )
+    subroutine cutoffNotify( this_rcut, this_rsw, this_rlist )
       
-      real(kind=dp), intent(in) :: this_rcut, this_rlist, this_ecr, this_est
+      real(kind=dp), intent(in) :: this_rcut, this_rsw, this_rlist
 
-      real(kind=dp) :: rtaper, rcut, rlist, ecr
+      real(kind=dp) :: rsw, rcut, rlist
       integer :: localError
-      
-
+      logical :: do_shift
+     
       rcut   = this_rcut
+      rsw    = this_rsw
       rlist  = this_rlist
-      ecr    = this_ecr
-      rtaper = this_ecr - this_est
 
+      if (rcut .lt. rsw) then
+         write(*,*) 'warning, rcut = ', rcut, ' but rsw = ', rsw
+      endif
 
+      if (rlist .lt. rcut) then
+         write(*,*) 'warning, rlist = ', rlist, ' but rcut = ', rcut
+      endif
+
+      do_shift = .false.
+      if (abs(rcut-rsw) .lt. 0.0001) then
+         do_shift = .true.
+      endif
+      
       call setRlistDF( rlist )
-      call setCutoffsDipole( ecr, rtaper )
-      call setCutoffsRF( ecr, rtaper )
-      call setCutoffLJ( rcut, localError )
+      call setCutoffsRF( rcut, rsw )
+      call setCutoffLJ( rcut, do_shift, localError )
+      call setCutoffEAM( rcut, localError)
+      call set_switch(GROUP_SWITCH, rsw, rcut)
 
     end subroutine cutoffNotify