diff --git a/geomagio/algorithm/SqDistAlgorithm.py b/geomagio/algorithm/SqDistAlgorithm.py
index c3686fb8c0542ab69c1ceb77f117b789a77e1bc5..1141a1fdc68218fdd42f1b38ad797728089bb0e7 100644
--- a/geomagio/algorithm/SqDistAlgorithm.py
+++ b/geomagio/algorithm/SqDistAlgorithm.py
@@ -436,6 +436,14 @@ class SqDistAlgorithm(Algorithm):
                 # performance.
                 r[i + 1] = gamma * (1 - alpha) * et / m
 
+                ## this properly smooths SQ locally; figure out if we need an
+                ## input parameter to control this, or if we just fix it to a
+                ## reasonable, but static value; also, figure out how/if r
+                ## needs to be changed  -EJR 8/2016
+                #s[i + m] = s[i] + gamma * (1 - alpha) * et / 61.
+                #s[i+m-30:i+m] = s[i+m-30:i+m] + gamma * (1-alpha) * et / 61.
+                #s[i+1:i+30+1] = s[i+1:i+30+1] + gamma * (1-alpha) * et / 61.
+
                 # update and append to s using equation-error formulation
                 s[i + m] = s[i] + gamma * (1 - alpha) * et