--- trunk/OOPSE/libmdtools/SimInfo.cpp	2004/06/01 17:15:43	1212
+++ trunk/OOPSE/libmdtools/SimInfo.cpp	2004/06/01 18:42:58	1214
@@ -453,13 +453,13 @@ void SimInfo::refreshSim(){
   
   isError = 0;
   
-  getFortranGroupArray(this, mfact, ngroup, groupList, groupStart);
+  getFortranGroupArrays(this, FglobalGroupMembership, mfact);
   //it may not be a good idea to pass the address of first element in vector
   //since c++ standard does not require vector to be stored continuously in meomory
   //Most of the compilers will organize the memory of vector continuously
   setFsimulation( &fInfo, &n_global, &n_atoms, identArray, &n_exclude, excl, 
                   &nGlobalExcludes, globalExcludes, molMembershipArray, 
-                  &mfact[0], &ngroup, &groupList[0], &groupStart[0], &isError);
+                  &mfact[0], &ngroup, &FglobalGroupMembership[0], &isError);
   
   if( isError ){
     
@@ -567,12 +567,13 @@ GenericData* SimInfo::getProperty(const string& propNa
 }
 
 
-void getFortranGroupArray(SimInfo* info, vector<double>& mfact, int& ngroup, 
-                          vector<int>& groupList, vector<int>& groupStart){
+void SimInfo::getFortranGroupArrays(SimInfo* info, 
+                                    vector<int>& FglobalGroupMembership,
+                                    vector<double>& mfact){
+  
   Molecule* myMols;
   Atom** myAtoms;
   int numAtom;
-  int curIndex;
   double mtot;
   int numMol;
   int numCutoffGroups;
@@ -584,39 +585,34 @@ void getFortranGroupArray(SimInfo* info, vector<double
   double totalMass;
   
   mfact.clear();
-  groupList.clear();
-  groupStart.clear();
+  FglobalGroupMembership.clear();
   
-  //Be careful, fortran array begin at 1
-  curIndex = 1;
 
+  // Fix the silly fortran indexing problem
+#ifdef IS_MPI
+  numAtom = mpiSim->getNAtomsGlobal();
+#else
+  numAtom = n_atoms;
+#endif
+  for (int i = 0; i < numAtom; i++) 
+    FglobalGroupMembership.push_back(globalGroupMembership[i] + 1);
+
   myMols = info->molecules;
   numMol = info->n_mol;
   for(int i  = 0; i < numMol; i++){
     numCutoffGroups = myMols[i].getNCutoffGroups();
-    for(myCutoffGroup =myMols[i].beginCutoffGroup(iterCutoff); myCutoffGroup != NULL; 
-                                                  myCutoffGroup =myMols[i].nextCutoffGroup(iterCutoff)){
+    for(myCutoffGroup =myMols[i].beginCutoffGroup(iterCutoff); 
+        myCutoffGroup != NULL; 
+        myCutoffGroup =myMols[i].nextCutoffGroup(iterCutoff)){
 
       totalMass = myCutoffGroup->getMass();
       
-      for(cutoffAtom = myCutoffGroup->beginAtom(iterAtom); cutoffAtom != NULL; 
-                                           cutoffAtom = myCutoffGroup->nextAtom(iterAtom)){
+      for(cutoffAtom = myCutoffGroup->beginAtom(iterAtom); 
+          cutoffAtom != NULL; 
+          cutoffAtom = myCutoffGroup->nextAtom(iterAtom)){
         mfact.push_back(cutoffAtom->getMass()/totalMass);
-#ifdef IS_MPI        
-        groupList.push_back(cutoffAtom->getGlobalIndex() + 1);
-#else
-        groupList.push_back(cutoffAtom->getIndex() + 1);
-#endif
       }  
-                              
-      groupStart.push_back(curIndex);
-      curIndex += myCutoffGroup->getNumAtom();
+    }
+  }
 
-    }//end for(myCutoffGroup =myMols[i].beginCutoffGroup(iterCutoff))
-
-  }//end for(int i  = 0; i < numMol; i++)
-
-
-  //The last cutoff group need more element to indicate the end of the cutoff 
-  ngroup = groupStart.size();
 }