diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py
index 9b07c00cf0bebf0d55e6606d15009a998ab3f9ad..7220ca87ba8c41018658e902bad090d7a1f94ab2 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py
@@ -28,29 +28,29 @@ datestr = date.strftime("%Y-%m-%d")
 sym.init_printing()
 # solver_tol = 6E-7
 use_case = "TP-TP-layered-soil-pure-dd"
-max_iter_num = 1000
+max_iter_num = 200
 FEM_Lagrange_degree = 1
 mesh_study = True
 resolutions = {
-                1: 1e-6,  # h=2
-                2: 1e-6,  # h=1.1180
-                4: 1e-6,  # h=0.5590
+                1: 5e-4,  # h=2
+                2: 5e-4,  # h=1.1180
+                4: 3e-4,  # h=0.5590
                 8: 1e-6,  # h=0.2814
                 16: 1e-6, # h=0.1412
                 32: 1e-6,   # h=0.0706
-                64: 1e-6,
-                128: 1e-6,
-                # 256: 1e-7
+                64: 1e-6,   # 0.03535
+                128: 1e-6,  # 0.01768
+                # 256: 1e-6  # 0.00884
                 }
 
 ############ GRID #######################
 # mesh_resolution = 20
 timestep_size = 0.005
-number_of_timesteps = 100
-plot_timestep_every = 2
+number_of_timesteps = 140
+plot_timestep_every = 1
 # decide how many timesteps you want analysed. Analysed means, that we write out
 # subsequent errors of the L-iteration within the timestep.
-number_of_timesteps_to_analyse = 0
+number_of_timesteps_to_analyse = 6
 starttimes = [0.0]
 
 Lw = 0.025 #/timestep_size