--- trunk/OOPSE/libmdtools/calc_dipole_dipole.F90 2003/07/15 22:22:41 619 +++ trunk/OOPSE/libmdtools/calc_dipole_dipole.F90 2003/07/16 02:11:02 621 @@ -31,7 +31,6 @@ contains return end subroutine initialize_dipole - subroutine do_dipole_pair(atom1, atom2, d, rij, r2, pot, u_l, f, t, & do_pot, do_stress) @@ -84,11 +83,8 @@ contains call getElementProperty(atypes, me1, "dipole_moment", mu1) call getElementProperty(atypes, me2, "dipole_moment", mu2) - if (rij.le.rrf) then - - if (rij.lt.rt) then taper = 1.0d0 dtdr = 0.0d0