--- trunk/OOPSE/libmdtools/SimSetup.cpp	2004/04/13 16:26:03	1104
+++ trunk/OOPSE/libmdtools/SimSetup.cpp	2004/04/14 15:37:41	1108
@@ -190,6 +190,8 @@ void SimSetup::makeMolecules(void){
   set<int> skipList;
 
   double phi, theta, psi;
+  char* molName;
+  char rbName[100];
 
   //init the forceField paramters
 
@@ -206,6 +208,7 @@ void SimSetup::makeMolecules(void){
 
     for (i = 0; i < info[k].n_mol; i++){
       stampID = info[k].molecules[i].getStampID();
+      molName = comp_stamps[stampID]->getID();
 
       molInfo.nAtoms = comp_stamps[stampID]->getNAtoms();
       molInfo.nBonds = comp_stamps[stampID]->getNBonds();
@@ -265,7 +268,7 @@ void SimSetup::makeMolecules(void){
 
 #ifdef IS_MPI
 
-        molInfo.myAtoms[j]->setGlobalIndex(globalIndex[j + atomOffset]);
+        molInfo.myAtoms[j]->setGlobalIndex(globalAtomIndex[j + atomOffset]);
 
 #endif // is_mpi
       } 
@@ -414,6 +417,9 @@ void SimSetup::makeMolecules(void){
         // Create the Rigid Body:
 
         myRB = new RigidBody();
+
+        sprintf(rbName,"%s_RB_%s", molName, j);
+        myRB->setType(rbName);
         
         for (rb1 = 0; rb1 < nMembers; rb1++) {
 
@@ -1346,10 +1352,14 @@ void SimSetup::mpiMolDivide(void){
   int localMol, allMol;
   int local_atoms, local_bonds, local_bends, local_torsions, local_SRI;
   int local_rigid;
+  vector<int> globalAtomIndex;
+  vector<int> globalMolIndex;
 
   mpiSim = new mpiSimulation(info);
 
-  globalIndex = mpiSim->divideLabor();
+  mpiSim->divideLabor();
+  globalAtomIndex = mpiSim->getGlobalAtomIndex();
+  globalMolIndex = mpiSim->getGlobalMolIndex();
 
   // set up the local variables 
 
@@ -1363,7 +1373,7 @@ void SimSetup::mpiMolDivide(void){
   local_bends = 0;
   local_torsions = 0;
   local_rigid = 0;
-  globalAtomIndex = 0;
+  globalAtomCounter = 0;
 
   for (i = 0; i < n_components; i++){
     for (j = 0; j < components_nmol[i]; j++){
@@ -1376,8 +1386,8 @@ void SimSetup::mpiMolDivide(void){
         localMol++;
       }      
       for (k = 0; k < comp_stamps[i]->getNAtoms(); k++){
-        info[0].molMembershipArray[globalAtomIndex] = allMol;
-        globalAtomIndex++;
+        info[0].molMembershipArray[globalAtomCounter] = allMol;
+        globalAtomCounter++;
       }
 
       allMol++;
@@ -1445,15 +1455,15 @@ void SimSetup::makeSysArrays(void){
 #else // is_mpi
 
     molIndex = 0;
-    globalAtomIndex = 0;
+    globalAtomCounter = 0;
     for (i = 0; i < n_components; i++){
       for (j = 0; j < components_nmol[i]; j++){
         the_molecules[molIndex].setStampID(i);
         the_molecules[molIndex].setMyIndex(molIndex);
         the_molecules[molIndex].setGlobalIndex(molIndex);
         for (k = 0; k < comp_stamps[i]->getNAtoms(); k++){
-          info[l].molMembershipArray[globalAtomIndex] = molIndex;
-          globalAtomIndex++;
+          info[l].molMembershipArray[globalAtomCounter] = molIndex;
+          globalAtomCounter++;
         }
         molIndex++;
       }