--- trunk/OOPSE/libmdtools/calc_dipole_dipole.F90 2003/03/21 17:42:12 378 +++ trunk/OOPSE/libmdtools/calc_dipole_dipole.F90 2003/03/21 22:11:50 388 @@ -23,7 +23,7 @@ contains rt = this_rt rrfsq = rrf * rrf ! pre converts from mu in units of debye to kcal/mol - pre = 14.38362d0 + pre = 14.38362_dp dipole_initialized = .true. @@ -81,7 +81,7 @@ 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