diff --git a/LDDsimulation/LDDsimulation.py b/LDDsimulation/LDDsimulation.py
index 69769d35a4218a9acd03c525c0d3530e171b79b6..421643d638f656f2b5fd2b37d0969fcba54ef550 100644
--- a/LDDsimulation/LDDsimulation.py
+++ b/LDDsimulation/LDDsimulation.py
@@ -409,7 +409,6 @@ class LDDsimulation(object):
         # print(timing_table.__get_value__())
         # with open(self.output_dir+"timings.txt", "w") as timefile:
         #     timefile.write(df.timings)
-
         df.dump_timings_to_xml(self.output_dir+"timings.xml", df.TimingClear.keep)
         # df.info(colored("start post processing calculations ...\n", "yellow"))
         # self.post_processing()
diff --git a/Two-phase-Two-phase/two-patch/mesh_studies/TP-TP-2-patch-mesh-study.py b/Two-phase-Two-phase/two-patch/mesh_studies/TP-TP-2-patch-mesh-study.py
index 90f26667d7055c72cce43c559a0b488d353135aa..889d50fe374ea121c4f99a147150d832f86a9960 100755
--- a/Two-phase-Two-phase/two-patch/mesh_studies/TP-TP-2-patch-mesh-study.py
+++ b/Two-phase-Two-phase/two-patch/mesh_studies/TP-TP-2-patch-mesh-study.py
@@ -458,6 +458,10 @@ for subdomain in isRichards.keys():
 #
 # sa
 
+f = open('TP-TP-2-patch-mesh-study.py', 'r')
+print(f.read())
+f.close()
+
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class