diff --git a/src/tespy/networks/network.py b/src/tespy/networks/network.py
index 45737ad6a..c0a7c3cf0 100644
--- a/src/tespy/networks/network.py
+++ b/src/tespy/networks/network.py
@@ -1121,7 +1121,7 @@ def presolve_fluid_topology(self):
                 }
                 self.num_conn_vars += 1
 
-    def _reset_topology_reduction_specifications(self):
+    def reset_topology_reduction_specifications(self):
         for c in self.conns["object"]:
             if hasattr(c, "_m_tmp"):
                 value = c.m.val_SI
@@ -1963,7 +1963,7 @@ def solve(self, mode, init_path=None, design_path=None,
         self.initialise()
 
         if init_only:
-            self._reset_topology_reduction_specifications()
+            self.reset_topology_reduction_specifications()
             return
 
         msg = 'Starting solver.'
@@ -1974,7 +1974,7 @@ def solve(self, mode, init_path=None, design_path=None,
         self.solve_loop(print_results=print_results)
 
         if not prepare_fast_lane:
-            self._reset_topology_reduction_specifications()
+            self.reset_topology_reduction_specifications()
 
         if self.lin_dep:
             msg = (
diff --git a/tutorial/advanced/starting_values.py b/tutorial/advanced/starting_values.py
index 357b0403e..ad6798fe5 100644
--- a/tutorial/advanced/starting_values.py
+++ b/tutorial/advanced/starting_values.py
@@ -107,7 +107,8 @@
     nw.solve("design")
 except ValueError as e:
     print(e)
-    nw._reset_topology_reduction_specifications()
+    nw.reset_topology_reduction_specifications()
+
 # %%[sec_4]
 import CoolProp.CoolProp as CP