diff --git a/geomagio/TimeseriesFactory.py b/geomagio/TimeseriesFactory.py
index 0096adeed16623550d449fc51d74aa9f153c5938..50adb2fd9159e45902f54d57f56c8c103e471e41 100644
--- a/geomagio/TimeseriesFactory.py
+++ b/geomagio/TimeseriesFactory.py
@@ -110,6 +110,8 @@ class TimeseriesFactory(object):
             try:
                 data = Util.read_url(url)
             except IOError as e:
+                print("Error reading url: %s, continuing" % str(e),
+                        file=sys.stderr)
                 continue
             try:
                 timeseries += self.parse_string(data,
diff --git a/geomagio/WebService.py b/geomagio/WebService.py
index e0e16610577b2087699b711ba2cd5dc5eb4e57b7..e304ade0831390d5077198feacd48dd1d9c62893 100644
--- a/geomagio/WebService.py
+++ b/geomagio/WebService.py
@@ -6,6 +6,7 @@ from cgi import escape, parse_qs
 from collections import OrderedDict
 from datetime import datetime
 from json import dumps
+import sys
 
 from geomagio.edge import EdgeFactory
 from geomagio.iaga2002 import IAGA2002Writer
@@ -109,6 +110,8 @@ class WebService(object):
             if isinstance(timeseries_string, str):
                 timeseries_string = timeseries_string.encode('utf8')
         except Exception as e:
+            print("Error processing request: %s" % str(e),
+                    file=sys.stderr)
             message = "Server error."
             error_body = self.error(500, message, environ, start_response)
             return [error_body]
diff --git a/geomagio/algorithm/SqDistAlgorithm.py b/geomagio/algorithm/SqDistAlgorithm.py
index 69db1d9138ecd1373d2ea0fa467c5d5fa31c977e..3822a73c7ac9a123746a8f155e309c1709462e81 100644
--- a/geomagio/algorithm/SqDistAlgorithm.py
+++ b/geomagio/algorithm/SqDistAlgorithm.py
@@ -11,7 +11,7 @@
     Use of fmin_l_bfgs_b to estimate parameters inspired by Andre Queiroz:
         https://gist.github.com/andrequeiroz/5888967
 """
-from __future__ import absolute_import
+from __future__ import absolute_import, print_function
 
 from .. import StreamConverter
 from .Algorithm import Algorithm
@@ -82,7 +82,7 @@ class SqDistAlgorithm(Algorithm):
             # state is up to date, only need new data
             return (start, end)
         # state not up to date, need to prime
-        print('WARNING: missing or incompatible state...' + \
+        print('WARNING: missing or incompatible state...' +
               'reinitializing.', file=sys.stderr)
         return (start - 3 * 30 * 24 * 60 * 60, end)