--- trunk/OOPSE/libmdtools/ZConstraint.cpp	2003/09/02 14:30:12	738
+++ trunk/OOPSE/libmdtools/ZConstraint.cpp	2003/09/15 16:52:02	763
@@ -279,8 +279,8 @@ template<typename T> ZConstraint<T>::ZConstraint(SimIn
       zPos.push_back((*parameters)[searchResult].zPos);
 //       cout << "index: "<< (*parameters)[searchResult].zconsIndex 
 //              <<"\tzPos = " << (*parameters)[searchResult].zPos << endl; 
-      kz.push_back((*parameters)[searchResult]. kRatio * zForceConst);
       
+      kz.push_back((*parameters)[searchResult]. kRatio * zForceConst);
       molecules[i].getCOM(COM);
     }
     else
@@ -394,7 +394,6 @@ template<typename T> void ZConstraint<T>::update()
       zconsMols.push_back(&molecules[i]);      
       zPos.push_back((*parameters)[index].zPos);
       kz.push_back((*parameters)[index].kRatio * zForceConst);
-      
       massOfZConsMols.push_back(molecules[i].getTotalMass());  
        
       molecules[i].getCOM(COM);
@@ -932,8 +931,8 @@ template<typename T> void ZConstraint<T>::doHarmonic()
   MPI_Allreduce(&totalFZ_local, &totalFZ, 1, MPI_DOUBLE,MPI_SUM, MPI_COMM_WORLD);  
 #endif
 
-  cout << "before substracting harmonic force from moving molecuels "
-        << "total force is " << calcTotalForce()  << endl;
+  //cout << "before substracting harmonic force from moving molecuels "
+  //      << "total force is " << calcTotalForce()  << endl;
 
   force[0]= 0;
   force[1]= 0;
@@ -951,8 +950,8 @@ template<typename T> void ZConstraint<T>::doHarmonic()
      }
   }   
 
-  cout << "after substracting harmonic force from moving molecuels "
-        << "total force is " << calcTotalForce()  << endl;
+  //cout << "after substracting harmonic force from moving molecuels "
+  //      << "total force is " << calcTotalForce()  << endl;
 
 }