diff --git a/Bioindustrial-Park b/Bioindustrial-Park
index 8f86dfec..b727e367 160000
--- a/Bioindustrial-Park
+++ b/Bioindustrial-Park
@@ -1 +1 @@
-Subproject commit 8f86dfecc3ac029fe4340fab647a69ae7de595ab
+Subproject commit b727e367911e8718afe986bfa14499bbc8bdcdf5
diff --git a/How2STEAM b/How2STEAM
index e9c0e3d0..3774b831 160000
--- a/How2STEAM
+++ b/How2STEAM
@@ -1 +1 @@
-Subproject commit e9c0e3d0d5ddf25e8d5c6b65a1f4eadfaa382511
+Subproject commit 3774b83195020387707aa7a7ddf14e09189a91c3
diff --git a/biosteam/units/aerated_bioreactor.py b/biosteam/units/aerated_bioreactor.py
index 51ecedcf..fd0c4d7e 100644
--- a/biosteam/units/aerated_bioreactor.py
+++ b/biosteam/units/aerated_bioreactor.py
@@ -20,7 +20,7 @@
 import numpy as np
 from scipy.constants import g
 import flexsolve as flx
-from warnings import catch_warnings
+from warnings import filterwarnings, catch_warnings
 from scipy.optimize import minimize_scalar, minimize, least_squares
 from biosteam.units.design_tools import aeration
 
@@ -582,7 +582,8 @@ def total_power_at_substrate_flow(F_substrates):
                 return total_power
             
             f = total_power_at_substrate_flow
-            with catch_warnings(action='ignore'):
+            with catch_warnings():
+                filterwarnings('ignore')
                 results = minimize(f, 1.2 * SURs, bounds=bounds, tol=SURs.max() * 1e-6)
                 load_flow_rates(results.x / x_substrates)
         else:
@@ -598,9 +599,11 @@ def gas_flow_rate_objective(F_substrates):
             
             f = gas_flow_rate_objective
             bounds = bounds.T
-            with catch_warnings(action='ignore'):
+            with catch_warnings():
+                filterwarnings('ignore')
                 results = least_squares(f, 1.2 * SURs, bounds=bounds, ftol=SURs.min() * 1e-6)
-                load_flow_rates(results.x / x_substrates)
+            load_flow_rates(results.x / x_substrates)
+            
         
         # self.show()
         # breakpoint()