--- branches/development/src/perturbations/ElectricField.cpp 2012/08/29 18:13:11 1787 +++ branches/development/src/perturbations/ElectricField.cpp 2013/01/25 15:21:04 1839 @@ -78,7 +78,7 @@ namespace OpenMD { Vector3d EFfrc; Vector3d pos; RealType chrg; - RealType pot, fieldPot, moment; + RealType pot, fieldPot; RealType chrgToKcal = 23.0609; RealType debyeToKcal = 4.8018969509; bool isCharge; @@ -92,6 +92,8 @@ namespace OpenMD { isCharge = false; chrg = 0.0; + atom->addElectricField(EF * chrgToKcal); + FixedChargeAdapter fca = FixedChargeAdapter(atom->getAtomType()); if ( fca.isFixedCharge() ) { isCharge = true; @@ -143,5 +145,4 @@ namespace OpenMD { snap->setLongRangePotential(longRangePotential); } } - }