diff --git a/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py b/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py
index 441c93be4843e5a2defe82e6a8133bceee54d258..a4d21306fad10c41d1b33f7cc7a9cc9e346e39ea 100755
--- a/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py
+++ b/TP-TP-2-patch-pure-dd-avoid-interface-at-origin/mesh_study_convergence/TP-TP-2-patch-pure-dd-convergence-study.py
@@ -562,11 +562,12 @@ for mesh_resolution in resolutions:
 
     simulation.initialise()
     # simulation.write_exact_solution_to_xdmf()
-    errornorms = simulation.run(analyse_condition=analyse_condition)
-    for subdomain_index in errornorms.keys():
-        for phase, different_errornorm in errornorms[subdomain_index].items():
-            for errortype, errornorm in errornorms[subdomain_index][phase].items():
-                eoc_filename = "{}_error".format(output_string)
-                eocfile = open("eoc", "a")
-                eocfile.write( str(dx) + " " + str(err) + "\n" )
+    output = simulation.run(analyse_condition=analyse_condition)
+    for subdomain_index in output.keys():
+        mesh_h = output['mesh_size']
+        for phase, different_errornorms in output[subdomain_index]['errornorm'].items():
+            for errortype, errornorm in different_errornorms.items():
+                eoc_filename = "{}_{}-errornorm-for-{}".format(output_string, errortype, phase)
+                eocfile = open("eoc_filename", "a")
+                eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
                 eocfile.close()