diff --git a/pymule/errortools.py b/pymule/errortools.py
index 2cbd25ad04a2ec2c0d202625b9a5f747c6da7ebe..acb36c1223ab315db55346d1cdea4bc5ecb8f0b4 100644
--- a/pymule/errortools.py
+++ b/pymule/errortools.py
@@ -91,3 +91,30 @@ def combineNplots(func, plots):
     for plot in plots[1:]:
         accum = func(accum, plot)
     return accum
+
+
+def mergebins(p, n):
+    if p[0,0] == -np.inf:
+        return np.concatenate([
+            [p[0,:]],
+            mergebins(p[1:,:], n)
+        ])
+    if p[-1,0] == +np.inf:
+        return np.concatenate([
+            mergebins(p[:-1,:], n),
+            [p[-1,:]]
+        ])
+
+    # Make sure the length fits
+    if len(p) % n:
+        short = p[:-(len(p) % n)]
+    else:
+        short = p[:]
+    part = np.split(short, len(short)/n)
+    return np.array([
+        (
+            sum(i[:,0])/n,
+            sum(i[:,1])/n,
+            np.sqrt(sum(i[:,2]**2))/n
+        ) for i in part
+    ])