--- trunk/src/brains/Thermo.cpp	2008/09/11 19:40:59	1291
+++ trunk/src/brains/Thermo.cpp	2008/10/22 20:01:49	1313
@@ -249,16 +249,21 @@ namespace oopse {
         Vector3d pos1, pos2, rab;
 
 #ifdef IS_MPI        
+        std::cerr << "tap = " << tap.first << "  " << tap.second << std::endl;
 
-        mol1 = info_.globalMolMembership_[tap.first];
-        mol2 = info_.globalMolMembership_[tap.second];
-        
+	int mol1 = info_->getGlobalMolMembership(tap.first);
+	int mol2 = info_->getGlobalMolMembership(tap.second);
+        std::cerr << "mols = " << mol1 << " " << mol2 << std::endl;
+
         int proc1 = info_->getMolToProc(mol1);
         int proc2 = info_->getMolToProc(mol2);
 
-        if (proc1 == worldRank) {
-          RealType data[3];
+        std::cerr << " procs = " << proc1 << " " <<proc2 <<std::endl;
+
+	RealType data[3];
+        if (proc1 == worldRank) {
           StuntDouble* sd1 = info_->getIOIndexToIntegrableObject(tap.first);
+          std::cerr << " on proc " << proc1 << ", sd1 has global index= " << sd1->getGlobalIndex() << std::endl;
           pos1 = sd1->getPos();
           data[0] = pos1.x();
           data[1] = pos1.y();
@@ -268,10 +273,11 @@ namespace oopse {
           MPI_Bcast(data, 3, MPI_REALTYPE, proc1, MPI_COMM_WORLD);
           pos1 = Vector3d(data);
         }
-          
+
+
         if (proc2 == worldRank) {
-          RealType data[3];
           StuntDouble* sd2 = info_->getIOIndexToIntegrableObject(tap.second);
+          std::cerr << " on proc " << proc2 << ", sd2 has global index= " << sd2->getGlobalIndex() << std::endl;
           pos2 = sd2->getPos();
           data[0] = pos2.x();
           data[1] = pos2.y();