diff --git a/simulations/scripts/run_non_adaptive_cyclod_5ns.py b/simulations/scripts/run_non_adaptive_cyclod_5ns.py
index fdb1c2a..2a67105 100644
--- a/simulations/scripts/run_non_adaptive_cyclod_5ns.py
+++ b/simulations/scripts/run_non_adaptive_cyclod_5ns.py
@@ -60,8 +60,8 @@ def main() -> None:
     calc.wait()
 
     # Analyse
-    calc.set_attr_values("_equilibrated", True)
-    calc.set_attr_values("_equil_time", 0)
+    calc.recursively_set_attr("_equilibrated", True)
+    calc.recursively_set_attr("_equil_time", 0)
     calc.analyse()
     calc.analyse_convergence()
     results_df = calc.get_results_df()