diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py
index fd212760a7106e1dcd88feab6f7ddbb7e765409a..bfb754e15fa1bdcce6fba3fd0597a7290c906054 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/TP-TP-layered_soil-different-intrinsic.py
@@ -47,8 +47,8 @@ resolutions = {
                 # 2: 1e-6,
                 # 4: 1e-6,
                 # 8: 1e-5,
-                # 16: 5e-6,
-                32: 2e-6,
+                16: 5e-6,
+                # 32: 2e-6,
                 # 64: 2e-6,
                 # 128: 1e-6,
                 # 256: 1e-6,
@@ -65,12 +65,12 @@ number_of_timesteps = 1
 
 # LDD scheme parameters  ######################################################
 Lw1 = 0.007 # /timestep_size
-Lnw1 = 0.005
+Lnw1 = 0.001
 Lw2 = 0.007  # /timestep_size
-Lnw2 = 0.005
-Lw3 = 0.0007  # /timestep_size
+Lnw2 = 0.001
+Lw3 = 0.007  # /timestep_size
 Lnw3 = 0.0005
-Lw4 = 0.0007  # /timestep_size
+Lw4 = 0.007  # /timestep_size
 Lnw4 = 0.0005
 
 lambda12_w = 0.5
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 a176766ddf9f28e91b533a00ff721a7aa50e3a49..64ccc3c5582119efe111e6cf1e88b3afec765285 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
@@ -37,7 +37,7 @@ thisfile = "TP-TP-layered_soil.py"
 
 # GENERAL SOLVER CONFIG  ######################################################
 # maximal iteration per timestep
-max_iter_num = 10
+max_iter_num = 1000
 FEM_Lagrange_degree = 1
 
 # GRID AND MESH STUDY SPECIFICATIONS  #########################################
@@ -47,8 +47,8 @@ resolutions = {
                 # 2: 1e-6,
                 # 4: 1e-6,
                 # 8: 1e-5,
-                # 16: 5e-6,
-                32: 5e-6,
+                16: 5e-6,
+                # 32: 2e-6,
                 # 64: 2e-6,
                 # 128: 1e-6,
                 # 256: 1e-6,
@@ -57,30 +57,31 @@ 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}
+starttimes = {0: 0.0, 1: 0.3, 2: 0.6, 3: 0.9}
+# starttimes = {0: 0.0}
 timestep_size = 0.001
-number_of_timesteps = 5
+number_of_timesteps = 1
 
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.025  # /timestep_size
-Lnw1 = Lw1
-Lw2 = 0.025  # /timestep_size
-Lnw2 = Lw2
-Lw3 = 0.025  # /timestep_size
-Lnw3 = Lw3
-Lw4 = 0.025  # /timestep_size
-Lnw4 = Lw4
-
-lambda12_w = 40
-lambda12_nw = 40
-lambda23_w = 40
-lambda23_nw = 40
-lambda34_w = 40
-lambda34_nw = 40
+Lw1 = 0.007
+Lnw1 = 0.005
+Lw2 = 0.007
+Lnw2 = 0.005
+Lw3 = 0.007
+Lnw3 = 0.005
+Lw4 = 0.007
+Lnw4 = 0.005
+
+lambda12_w = 0.05
+lambda12_nw = 0.05
+lambda23_w = 0.05
+lambda23_nw = 0.05
+lambda34_w = 0.05
+lambda34_nw = 0.05
 
 include_gravity = False
-debugflag = True
+debugflag = False
 analyse_condition = False
 
 # I/O CONFIG  #################################################################
@@ -385,7 +386,7 @@ if __name__ == '__main__':
             #     parameter=simulation_parameter
             #     )
 
-        LDDsim.join()
+        # LDDsim.join()
         if mesh_study:
             simulation_output_dir = processQueue.get()
             hlp.merge_spacetime_errornorms(isRichards=isRichards,