diff --git a/docs/source/pyplots/guides/output_file.py b/docs/source/pyplots/guides/output_file.py
index b022adb9..8c75e832 100644
--- a/docs/source/pyplots/guides/output_file.py
+++ b/docs/source/pyplots/guides/output_file.py
@@ -40,7 +40,7 @@
     delta.uw, cmap='plasma')
 ax[1].set_xticks([])
 ax[1].set_yticks([])
-ax[0].set_title('velocity')
+ax[1].set_title('velocity')
 
 plt.tight_layout()
 plt.show()
diff --git a/pyDeltaRCM/_version.py b/pyDeltaRCM/_version.py
index ccd1061d..50c8f67d 100644
--- a/pyDeltaRCM/_version.py
+++ b/pyDeltaRCM/_version.py
@@ -4,4 +4,4 @@ def __version__():
     Private version declaration, gets assigned to pyDeltaRCM.__version__
     during import
     """
-    return '2.1.0a'
+    return '2.1.1'
diff --git a/pyDeltaRCM/init_tools.py b/pyDeltaRCM/init_tools.py
index f8e1d2a5..77070c72 100644
--- a/pyDeltaRCM/init_tools.py
+++ b/pyDeltaRCM/init_tools.py
@@ -54,8 +54,8 @@ def init_logger(self):
         self.logger.setLevel(logging.INFO)
 
         # create the logging file handler
-        fh = logging.FileHandler(
-            self.prefix_abspath + '/pyDeltaRCM_' + timestamp + '.log')
+        fh = logging.FileHandler(os.path.join(
+            self.prefix_abspath, 'pyDeltaRCM_' + timestamp + '.log'))
         formatter = logging.Formatter(
             '%(asctime)s - %(levelname)s - %(message)s')
         fh.setFormatter(formatter)
diff --git a/pyDeltaRCM/model.py b/pyDeltaRCM/model.py
index 2abf4146..62c674d6 100644
--- a/pyDeltaRCM/model.py
+++ b/pyDeltaRCM/model.py
@@ -1353,6 +1353,7 @@ def time_step(self, new_time_step):
             self.sigma = (self.sigma / self._dt) * new_time_step
 
         self._dt = new_time_step
+        self.init_sediment_routers()
 
     @property
     def time_iter(self):