diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py
index 9f3899cc074378794e662fb2cca1510c20c2c1d1..618d2e15f9caf0374d6e346db0b9ff265ebfa121 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py
@@ -102,19 +102,19 @@ lambda56_w = 4
 lambda56_nw = 4
 
 include_gravity = False
-debugflag = True
+debugflag = False
 analyse_condition = False
 
-if mesh_study:
-    output_string = "./output/{}-{}_timesteps{}_P{}".format(
-        datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
-        )
-else:
-    for tol in resolutions.values():
-        solver_tol = tol
-    output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
-        datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
-        )
+# if mesh_study:
+output_string = "./output/{}-{}_timesteps{}_P{}".format(
+    datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
+    )
+# else:
+#     # for tol in resolutions.values():
+#     #     solver_tol = tol
+#     output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
+#         datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
+#         )
 
 
 # toggle what should be written to files
diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py
index 38789d7997eacc3fe1ca81a22ac2f1b7aab0d71d..c4652a277e4f717d18fc99b3ef3fd6cf9dfb9af6 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py
@@ -30,26 +30,27 @@ sym.init_printing()
 # solver_tol = 6E-7
 use_case = "TP-R-layered-soil-with-inner-patch-all-params-one"
 # name of this very file. Needed for log output.
-thisfile = "TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py"
+thisfile = "TP-R-layered_soil_with_inner_patch-all-params-one.py"
 
-max_iter_num = 300
+max_iter_num = 150
 FEM_Lagrange_degree = 1
-mesh_study = True
+mesh_study = False
 resolutions = {
-                1: 5e-6,  # h=2
-                2: 5e-6,  # h=1.1180
-                4: 5e-6,  # h=0.5590
-                8: 5e-6,  # h=0.2814
-                16: 5e-6, # h=0.1412
+                # 1: 2e-6,  # h=2
+                # 2: 2e-6,  # h=1.1180
+                # 4: 2e-6,  # h=0.5590
+                # 8: 2e-6,  # h=0.2814
+                # 16: 7e-6, # h=0.1412
                 32: 5e-6,
-                64: 5e-6,
-                128: 5e-6
+                64: 2e-6,
+                128: 2e-6,
+                256: 2e-6
                 }
 
 # GRID #######################
 # mesh_resolution = 20
-timestep_size = 0.001
-number_of_timesteps = 1000
+timestep_size = 0.00025
+number_of_timesteps = 4000
 plot_timestep_every = 4
 # decide how many timesteps you want analysed. Analysed means, that we write
 # out subsequent errors of the L-iteration within the timestep.
@@ -102,19 +103,19 @@ lambda56_w = 4
 lambda56_nw = 4
 
 include_gravity = False
-debugflag = True
-analyse_condition = True
-
-if mesh_study:
-    output_string = "./output/{}-{}_timesteps{}_P{}".format(
-        datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
-        )
-else:
-    for tol in resolutions.values():
-        solver_tol = tol
-    output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
-        datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
-        )
+debugflag = False
+analyse_condition = False
+
+# if mesh_study:
+output_string = "./output/{}-{}_timesteps{}_P{}".format(
+    datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
+    )
+# else:
+#     # for tol in resolutions.values():
+#     #     solver_tol = tol
+#     output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
+#         datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
+#         )
 
 
 # toggle what should be written to files