--- trunk/OOPSE/libmdtools/ZConstraint.cpp	2003/08/20 14:34:04	701
+++ trunk/OOPSE/libmdtools/ZConstraint.cpp	2003/08/20 14:50:32	702
@@ -525,6 +525,7 @@ template<typename T> void ZConstraint<T>::calcForce(in
   double zsys;
   double COM[3];
   double force[3];
+  double zSysCOMVel;
 
   T::calcForce(calcPot, calcStress);
 
@@ -532,14 +533,20 @@ template<typename T> void ZConstraint<T>::calcForce(in
     zeroOutVel();
    forcePolicy->update();
   }  
+
   zsys = calcZSys();
-  cout << "---------------------------------------------------------------------" <<endl;
-  cout << "current time: " << info->getTime() << endl;
-  cout << "center of mass at z: " << zsys << endl;     
-  //cout << "before calcForce, the COMVel of moving molecules is " << calcMovingMolsCOMVel() <<endl;
-  cout << "before calcForce, the COMVel of system is " << calcSysCOMVel() <<endl;
+  zSysCOMVel = calcSysCOMVel();
+#ifdef IS_MPI
+  if(worldRank == 0){
+#endif
+    cout << "---------------------------------------------------------------------" <<endl;
+    cout << "current time: " << info->getTime() << endl;
+    cout << "center of mass at z: " << zsys << endl;     
+    cout << "before calcForce, the COMVel of system is " << zSysCOMVel <<endl;
 
-  //cout <<  "before doZConstraintForce, totalForce is " << calcTotalForce() << endl;
+#ifdef IS_MPI
+  }
+#endif
 
   //do zconstraint force; 
   if (haveFixedZMols())
@@ -549,8 +556,6 @@ template<typename T> void ZConstraint<T>::calcForce(in
   if (haveMovingZMols())
     this->doHarmonic();
 
-  //cout <<  "after doHarmonic, totalForce is " << calcTotalForce() << endl;
-
   //write out forces and current positions of z-constraint molecules
   if(info->getTime() >= curZconsTime){    
    for(int i = 0; i < zconsMols.size(); i++){
@@ -571,9 +576,16 @@ template<typename T> void ZConstraint<T>::calcForce(in
     fzOut->writeFZ(info->getTime(), zconsMols.size(), indexOfZConsMols, fz, curZPos);
    curZconsTime += zconsTime;
   }
-  
-  //cout << "after calcForce, the COMVel of moving molecules is " << calcMovingMolsCOMVel() <<endl; 
-  cout << "after calcForce, the COMVel of system is " << calcSysCOMVel() <<endl; 
+
+  zSysCOMVel = calcSysCOMVel();  
+#ifdef IS_MPI
+  if(worldRank == 0){
+#endif
+    cout << "after calcForce, the COMVel of system is " << zSysCOMVel <<endl;
+#ifdef IS_MPI
+  }
+#endif
+
 }
 
 
@@ -786,7 +798,7 @@ template<typename T> void ZConstraint<T>::doZconstrain
       } 
       totalFZ_local += fz[i];
 
-      cout << "Fixed Molecule --\tindex: " << indexOfZConsMols[i] 
+      cout << "Fixed Molecule\tindex: " << indexOfZConsMols[i] 
              <<"\tcurrent zpos: " << COM[whichDirection] 
              << "\tcurrent fz: " <<fz[i] << endl;
 
@@ -903,7 +915,7 @@ template<typename T> void ZConstraint<T>::doHarmonic()
 
     if (states[i] == zcsMoving){
       zconsMols[i]->getCOM(COM);
-      cout << "Moving Molecule --\tindex: " << indexOfZConsMols[i] <<"\tcurrent zpos: " << COM[whichDirection] << endl;
+      cout << "Moving Molecule\tindex: " << indexOfZConsMols[i] <<"\tcurrent zpos: " << COM[whichDirection] << endl;
     
     diff = COM[whichDirection] -zPos[i];