--- trunk/OOPSE/libmdtools/NPT.cpp	2003/10/29 00:19:10	837
+++ trunk/OOPSE/libmdtools/NPT.cpp	2003/11/07 17:09:48	857
@@ -1,4 +1,5 @@
 #include <math.h>
+
 #include "Atom.hpp"
 #include "SRI.hpp"
 #include "AbstractClasses.hpp"
@@ -97,6 +98,9 @@ template<typename T> void NPT<T>::moveA() {
   tStats->getCOM(COM);
 
   //evolve velocity half step
+
+  calcVelScale();
+  
   for( i=0; i<nAtoms; i++ ){
 
     atoms[i]->getVel( vel );
@@ -223,6 +227,7 @@ template<typename T> void NPT<T>::moveB( void ){
 
     this->evolveChiB();
     this->evolveEtaB();
+    this->calcVelScale();
 
     for( i=0; i<nAtoms; i++ ){