--- trunk/OOPSE/libmdtools/SimInfo.cpp	2004/05/12 15:02:03	1164
+++ trunk/OOPSE/libmdtools/SimInfo.cpp	2004/05/27 00:48:12	1198
@@ -62,7 +62,8 @@ SimInfo::SimInfo(){
   useReactionField = 0;
   useGB = 0;
   useEAM = 0;
-  
+  useThermInt = 0;
+
   haveCutoffGroups = false;
 
   excludes = Exclude::Instance();
@@ -591,29 +592,6 @@ void getFortranGroupArray(SimInfo* info, vector<double
   myMols = info->molecules;
   numMol = info->n_mol;
   for(int i  = 0; i < numMol; i++){
-    numAtom = myMols[i].getNAtoms();
-    myAtoms = myMols[i].getMyAtoms();
-
-    
-    for(int j = 0; j < numAtom; j++){
-
-    
-#ifdef IS_MPI      
-      atomIndex = myAtoms[j]->getGlobalIndex();
-#else
-      atomIndex = myAtoms[j]->getIndex();
-#endif
-
-      if(myMols[i].belongToCutoffGroup(atomIndex))
-        continue;
-      else{
-        // this is a fraction of the cutoff group's mass, not the mass itself!
-        mfact.push_back(1.0);
-        groupList.push_back(myAtoms[j]->getIndex() + 1);
-        groupStart.push_back(curIndex++);   
-      }
-    }
-      
     numCutoffGroups = myMols[i].getNCutoffGroups();
     for(myCutoffGroup =myMols[i].beginCutoffGroup(iterCutoff); myCutoffGroup != NULL; 
                                                   myCutoffGroup =myMols[i].nextCutoffGroup(iterCutoff)){
@@ -632,6 +610,9 @@ void getFortranGroupArray(SimInfo* info, vector<double
     }//end for(myCutoffGroup =myMols[i].beginCutoffGroup(iterCutoff))
 
   }//end for(int i  = 0; i < numMol; i++)
-  
-  ngroup = groupStart.size();
+
+
+  //The last cutoff group need more element to indicate the end of the cutoff 
+  groupStart.push_back(curIndex);
+  ngroup = groupStart.size() - 1;
 }