diff --git a/gillespy2/solvers/numpy/tau_hybrid_solver.py b/gillespy2/solvers/numpy/tau_hybrid_solver.py
index 1d1a270d3..b9018e6e2 100644
--- a/gillespy2/solvers/numpy/tau_hybrid_solver.py
+++ b/gillespy2/solvers/numpy/tau_hybrid_solver.py
@@ -430,7 +430,7 @@ def __integrate(self, integrator, integrator_options, curr_state, y0, model, cur
                     model.listOfAssignmentRules]
         rhs = lambda t, y: TauHybridSolver.__f(t, y, *int_args)
         if 'min_step' in integrator_options:
-            tau_step = max(min_step, tau_step)
+            tau_step = max(integrator_options['min_step'], tau_step)
         else:
             tau_step = max(1e-6, tau_step)
         if pure_ode: