diff --git a/src/postProcessing/functionObjects/forces/forces/forces.C b/src/postProcessing/functionObjects/forces/forces/forces.C
index 539e14b..0c4fc15 100644
--- a/src/postProcessing/functionObjects/forces/forces/forces.C
+++ b/src/postProcessing/functionObjects/forces/forces/forces.C
@@ -378,12 +378,16 @@ void Foam::forces::applyBins
         forAll(dd, i)
         {
             label binI = floor(dd[i]/binDx_);
-            force_[0][binI] += fN[i];
-            force_[1][binI] += fT[i];
-            force_[2][binI] += fP[i];
-            moment_[0][binI] += Md[i]^fN[i];
-            moment_[1][binI] += Md[i]^fT[i];
-            moment_[2][binI] += Md[i]^fP[i];
+            
+            if(0<=binI && binI<force_.size())
+            {
+              force_[0][binI] += fN[i];
+              force_[1][binI] += fT[i];
+              force_[2][binI] += fP[i];
+              moment_[0][binI] += Md[i]^fN[i];
+              moment_[1][binI] += Md[i]^fT[i];
+              moment_[2][binI] += Md[i]^fP[i];
+            }
         }
     }
 }
