diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py
index ec01e7288e1d79924de47f196a6f29c4e8bc5396..dd7584a7848da27eb0a17d6568483c13362712dd 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py
@@ -28,7 +28,7 @@ datestr = date.strftime("%Y-%m-%d")
 sym.init_printing()
 # solver_tol = 6E-7
 use_case = "TP-TP-layered-soil-realistic"
-max_iter_num = 50
+max_iter_num = 1000
 FEM_Lagrange_degree = 1
 mesh_study = False
 resolutions = {
@@ -45,11 +45,11 @@ resolutions = {
 ############ GRID #######################
 # mesh_resolution = 20
 timestep_size = 0.005
-number_of_timesteps = 15
+number_of_timesteps = 240
 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 = 8
 starttime = 0.0
 
 Lw = 0.025 #/timestep_size
@@ -59,8 +59,8 @@ lambda_w = 40
 lambda_nw = 40
 
 include_gravity = False
-debugflag = True
-analyse_condition = False
+debugflag = False
+analyse_condition = True
 
 if mesh_study:
     output_string = "./output/{}-{}_timesteps{}_P{}".format(datestr, use_case, number_of_timesteps, FEM_Lagrange_degree)