--- trunk/OOPSE/libmdtools/ZConsWriter.cpp	2003/09/25 19:27:15	787
+++ trunk/OOPSE/libmdtools/ZConsWriter.cpp	2004/03/16 19:22:56	1091
@@ -54,14 +54,14 @@ ZConsWriter::~ZConsWriter()
 /**
  *
  */
-void ZConsWriter::writeFZ(double time, int num, int* index, double* fz, double* curZPos){
+void ZConsWriter::writeFZ(double time, int num, int* index, double* fz, double* curZPos, double* zpos){
 
 #ifndef IS_MPI
   output << time << endl;
   output << num << endl;
   
   for(int i = 0; i < num; i++)
-    output << index[i] <<"\t" << fz[i] << "\t" << curZPos[i] << endl;
+    output << index[i] <<"\t" << fz[i] << "\t" << curZPos[i] << "\t" << zpos[i] <<endl;
 
 #else
   int totalNum;
@@ -76,6 +76,7 @@ void ZConsWriter::writeFZ(double time, int num, int* i
   enum CommType { RequesPosAndForce, EndOfRequest} status;
   double pos;
   double force;
+  double zconsPos;
   int localIndex;
   MPI_Status ierr;
   int tag = 0;
@@ -107,11 +108,12 @@ void ZConsWriter::writeFZ(double time, int num, int* i
         status = RequesPosAndForce;
         MPI_Send(&status, 1, MPI_INT, whichNode, tag, MPI_COMM_WORLD);
         MPI_Send(&globalIndexOfCurMol, 1, MPI_INT, whichNode, tag, MPI_COMM_WORLD);
-        MPI_Recv(&force, 1, MPI_DOUBLE_PRECISION, whichNode, tag, MPI_COMM_WORLD, &ierr);
-        MPI_Recv(&pos, 1, MPI_DOUBLE_PRECISION, whichNode, tag, MPI_COMM_WORLD, &ierr);
+        MPI_Recv(&force, 1, MPI_DOUBLE, whichNode, tag, MPI_COMM_WORLD, &ierr);
+        MPI_Recv(&pos, 1, MPI_DOUBLE, whichNode, tag, MPI_COMM_WORLD, &ierr);
+        MPI_Recv(&zconsPos, 1, MPI_DOUBLE, whichNode, tag, MPI_COMM_WORLD, &ierr);
       }
 
-     output << globalIndexOfCurMol << "\t" << force << "\t" << pos << endl;
+     output << globalIndexOfCurMol << "\t" << force << "\t" << pos << "\t"<<  zconsPos << endl;
               
     } //End of Request Loop
     
@@ -142,8 +144,9 @@ void ZConsWriter::writeFZ(double time, int num, int* i
              break;
            }
     
-           MPI_Send(&fz[localIndex], 1, MPI_DOUBLE_PRECISION, 0, tag, MPI_COMM_WORLD);    
-           MPI_Send(&curZPos[localIndex], 1, MPI_DOUBLE_PRECISION, 0, tag, MPI_COMM_WORLD);       
+           MPI_Send(&fz[localIndex], 1, MPI_DOUBLE, 0, tag, MPI_COMM_WORLD);    
+           MPI_Send(&curZPos[localIndex], 1, MPI_DOUBLE, 0, tag, MPI_COMM_WORLD);    
+           MPI_Send(&zpos[localIndex], 1, MPI_DOUBLE, 0, tag, MPI_COMM_WORLD);     
            break;
        
         case EndOfRequest :