--- trunk/OOPSE/libmdtools/NPTi.cpp	2003/11/06 22:01:37	855
+++ trunk/OOPSE/libmdtools/NPTi.cpp	2003/11/07 17:09:48	857
@@ -70,16 +70,20 @@ template<typename T> void NPTi<T>::getVelScaleA(double
 		 (p_convert*NkBT*tb2));
 }
 
+template<typename T> void NPTi<T>::calcVelScale(void) {
+  vScale = chi + eta;
+}
+
 template<typename T> void NPTi<T>::getVelScaleA(double sc[3], double vel[3]) {
   int i;
 
-  for(i=0; i<3; i++) sc[i] = vel[i] * ( chi + eta );
+  for(i=0; i<3; i++) sc[i] = vel[i] * vScale;
 }
 
 template<typename T> void NPTi<T>::getVelScaleB(double sc[3], int index ){
   int i;
 
-  for(i=0; i<3; i++) sc[i] = oldVel[index*3 + i] * ( chi + eta );
+  for(i=0; i<3; i++) sc[i] = oldVel[index*3 + i] * vScale;
 }