diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py
index 04f7fc22c74157316698b5f191b06ee330a8aa3d..ed44ca41e08e0c7f2a571ef87d313b5187bcf174 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil.py
@@ -57,27 +57,27 @@ resolutions = {
 # 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, 1: 0.3, 2: 0.6, 3: 0.9}
-# starttimes = {0: 0.0}
+# starttimes = {0: 0.0, 1: 0.3, 2: 0.6, 3: 0.9}
+starttimes = {0: 0.0}
 timestep_size = 0.001
-number_of_timesteps = 1
+number_of_timesteps = 1500
 
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.007
+Lw1 = 0.0025
 Lnw1 = 0.005
-Lw2 = 0.007
+Lw2 = 0.0025
 Lnw2 = 0.005
-Lw3 = 0.007
+Lw3 = 0.0025
 Lnw3 = 0.005
-Lw4 = 0.007
+Lw4 = 0.0025
 Lnw4 = 0.005
 
-lambda12_w = 4
+lambda12_w = 1
 lambda12_nw = 2
-lambda23_w = 4
+lambda23_w = 1
 lambda23_nw = 2
-lambda34_w = 4
+lambda34_w = 1
 lambda34_nw = 2
 
 include_gravity = False
@@ -88,10 +88,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 = 1
+plot_timestep_every = 5
 # 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 = 1
+number_of_timesteps_to_analyse = 10
 
 # 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.