diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-realistic-same-intrinsic.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-realistic-same-intrinsic.py
index cb6372a665985099e1e0b064bbbc38be26984370..667b6828a954f03b60b34741e68805a8b5f02ce3 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-realistic-same-intrinsic.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-realistic-same-intrinsic.py
@@ -30,35 +30,35 @@ date = datetime.datetime.now()
 datestr = date.strftime("%Y-%m-%d")
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-R-layered-soil-with-inner-patch-realistic-same-intrinsic"
+use_case = "TP-R-layered-soil-with-inner-patch-realistic-same-intrinsic-write-iterations"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-layered_soil_with_inner_patch-realistic-same-intrinsic.py"
 
 # GENERAL SOLVER CONFIG  ######################################################
 # maximal iteration per timestep
-max_iter_num = 1000
+max_iter_num = 2000
 FEM_Lagrange_degree = 1
 
 # GRID AND MESH STUDY SPECIFICATIONS  #########################################
-mesh_study = True
+mesh_study = False
 resolutions = {
-                1: 2e-6,  # h=2
-                2: 2e-6,  # h=1.1180
-                4: 2e-6,  # h=0.5590
-                8: 2e-6,  # h=0.2814
-                16: 1e-6, # h=0.1412
-                32: 1e-6,
-                64: 5e-7,
-                128: 5e-7
+                # 1: 2e-6,  # h=2
+                # 2: 2e-6,  # h=1.1180
+                # 4: 2e-6,  # h=0.5590
+                # 8: 2e-6,  # h=0.2814
+                # 16: 1e-6, # h=0.1412
+                32: 7e-7,
+                # 64: 5e-7,
+                # 128: 5e-7
                 }
 
 # starttimes gives a list of starttimes to run the simulation from.
 # The list is looped over and a simulation is run with t_0 as initial time
 #  for each element t_0 in starttimes.
-starttimes = {0: 0.0}
+starttimes = {450: 0.45}
 # starttimes = {0: 0.0, 1:0.3, 2:0.6, 3:0.9}
 timestep_size = 0.001
-number_of_timesteps = 1000
+number_of_timesteps = 50
 
 # LDD scheme parameters  ######################################################
 
@@ -115,10 +115,10 @@ analyse_condition = False
 # when number_of_timesteps is high, it might take a long time to write all
 # timesteps to disk. Therefore, you can choose to only write data of every
 # plot_timestep_every timestep to disk.
-plot_timestep_every = 3
+plot_timestep_every = 1
 # Decide how many timesteps you want analysed. Analysed means, that
 # subsequent errors of the L-iteration within the timestep are written out.
-number_of_timesteps_to_analyse = 5
+number_of_timesteps_to_analyse = 2
 
 # fine grained control over data to be written to disk in the mesh study case
 # as well as for a regular simuation for a fixed grid.
@@ -149,8 +149,8 @@ else:
     write_to_file = {
         'space_errornorms': True,
         'meshes_and_markers': True,
-        'L_iterations_per_timestep': False,
-        'solutions': True,
+        'L_iterations_per_timestep': True,
+        'solutions': False,
         'absolute_differences': True,
         'condition_numbers': analyse_condition,
         'subsequent_errors': True