--- trunk/src/applications/sequentialProps/ContactAngle2.cpp 2015/03/07 21:41:51 2071 +++ trunk/src/applications/sequentialProps/ContactAngle2.cpp 2015/03/07 23:52:07 2073 @@ -122,11 +122,12 @@ namespace OpenMD { // z is possibly symmetric around 0 z = pos.z(); - std::size_t whichRBin = int(r / dr); - std::size_t whichZBin = int( (zLen/2.0 + z) / dz); - - if ((whichRBin < nRBins_) && (whichZBin >= 0) && (whichZBin < nZBins_)) + int whichRBin = int(r / dr); + int whichZBin = int( (zLen/2.0 + z) / dz); + + if ((whichRBin < int(nRBins_)) && (whichZBin >= 0) && (whichZBin < int(nZBins_))) { histo[whichRBin][whichZBin] += sd->getMass(); + } }