--- trunk/OOPSE/libmdtools/InitializeFromFile.cpp	2004/04/22 21:33:55	1131
+++ trunk/OOPSE/libmdtools/InitializeFromFile.cpp	2004/06/11 17:16:21	1268
@@ -229,7 +229,7 @@ void InitializeFromFile :: readInit( SimInfo* the_simn
       simError();
     }
 
-    for (i=0 ; i < mpiSim->getTotNmol(); i++) {
+    for (i=0 ; i < mpiSim->getNMolGlobal(); i++) {
       which_node = MolToProcMap[i];
       if(which_node == 0){
        //molecules belong to master node
@@ -308,7 +308,7 @@ void InitializeFromFile :: readInit( SimInfo* the_simn
       simError();
     }
   
-    for (i=0 ; i < mpiSim->getTotNmol(); i++) {
+    for (i=0 ; i < mpiSim->getNMolGlobal(); i++) {
       which_node = MolToProcMap[i];
       
       if(which_node == worldRank){
@@ -518,6 +518,12 @@ char* InitializeFromFile::parseDumpLine(char* readLine
 
     qSqr = (q[0] * q[0]) + (q[1] * q[1]) + (q[2] * q[2]) + (q[3] * q[3]);
 
+    if (fabs(qSqr) < 1e-6) {
+      sprintf(painCave.errMsg,
+          "initial quaternion error (q0^2 + q1^2 + q2^2 + q3^2 ~ 0).\n");
+       return strdup(painCave.errMsg);
+    }
+    
     qLength = sqrt( qSqr );
     q[0] = q[0] / qLength;
     q[1] = q[1] / qLength;
@@ -649,7 +655,7 @@ void initFile::nodeZeroError( void ){
   int j, myStatus;
 
   myStatus = 0;
-  for (j = 0; j < mpiSim->getNumberProcessors(); j++) {
+  for (j = 0; j < mpiSim->getNProcessors(); j++) {
     MPI_Send( &myStatus, 1, MPI_INT, j,
 	      TAKE_THIS_TAG_INT, MPI_COMM_WORLD);
   }