diff --git a/src/balance.cpp b/src/balance.cpp
index 79d7ec97b7..dcbbf665f8 100644
--- a/src/balance.cpp
+++ b/src/balance.cpp
@@ -488,13 +488,15 @@ int *Balance::bisection(int sortflag)
   rcb->invert(sortflag);
   
   // NOTE: this logic is specific to orthogonal boxes, not triclinic
-  // insure comm->rcbcut will be exactly equal to sub-box boundaries it touches
+  // NOTE: check this: insure comm->rcbcut will be exactly equal to sub-box boundaries it touches
 
   comm->rcbnew = 1;
 
   int idim = rcb->cutdim;
-  double split = (rcb->cut - boxlo[idim]) / prd[idim];
-  comm->rcbcut = boxlo[idim] + prd[idim]*split;
+  //double split = (rcb->cut - boxlo[idim]) / prd[idim];
+  //comm->rcbcut = boxlo[idim] + prd[idim]*split;
+
+  comm->rcbcutfrac = (rcb->cut - boxlo[idim]) / prd[idim];
   comm->rcbcutdim = idim;
 
   double (*mysplit)[2] = comm->mysplit;