--- branches/development/src/primitives/DirectionalAtom.cpp	2012/08/29 18:13:11	1787
+++ branches/development/src/primitives/DirectionalAtom.cpp	2013/01/22 16:20:11	1838
@@ -106,10 +106,10 @@ namespace OpenMD {
   
   void DirectionalAtom::setA(const RotMat3x3d& a) {
     ((snapshotMan_->getCurrentSnapshot())->*storage_).aMat[localIndex_] = a;
- 
+
     if (atomType_->isMultipole()) {
       RotMat3x3d atrans = a.transpose();
-      
+
       if (atomType_->isDipole()) {
         ((snapshotMan_->getCurrentSnapshot())->*storage_).dipole[localIndex_] = atrans * dipole_;
       }
@@ -147,7 +147,8 @@ namespace OpenMD {
     Vector3d force;
     Vector3d torque;
     Vector3d myEuler;
-    RealType phi, theta, psi;
+    RealType phi, theta;
+    // RealType psi;
     RealType cphi, sphi, ctheta, stheta;
     Vector3d ephi;
     Vector3d etheta;
@@ -159,7 +160,7 @@ namespace OpenMD {
     
     phi = myEuler[0];
     theta = myEuler[1];
-    psi = myEuler[2];
+    // psi = myEuler[2];
     
     cphi = cos(phi);
     sphi = sin(phi);