diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-gravity-but-same-intrinsic-perm.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-gravity-but-same-intrinsic-perm.py
index ba5052cdfc0fb8590d3753f3f1733f92ca455ae9..d83a480e7de26030621fad994f53452e6600827e 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-gravity-but-same-intrinsic-perm.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-gravity-but-same-intrinsic-perm.py
@@ -30,7 +30,7 @@ date = datetime.datetime.now()
 datestr = date.strftime("%Y-%m-%d")
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-R-2-patch-realistic-same-intrinsic-perm-maybe-complement-08-07"
+use_case = "TP-R-2-patch-realistic-same-intrinsic-perm-control"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-2-patch-realistic-gravity-but-same-intrinsic-perm.py"
 
@@ -47,8 +47,9 @@ resolutions = {
                 # 4: 1e-5,
                 # 8: 1e-5,
                 # 16: 5e-6,
-                # 32: 2e-6,
-                64: 2e-6,
+                32: 5e-6,
+                32: 2e-6,
+                # 64: 2e-6,
                 # 128: 2e-6,
                 # 256: 1e-6,
                 }
@@ -57,19 +58,19 @@ resolutions = {
 # 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 = {1100: 1.1}
+starttimes = {1: 0.0}
 timestep_size = 0.001
-number_of_timesteps = 400
+number_of_timesteps = 1500
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.6
-Lnw1 = 0.5
+Lw1 = 0.25
+Lnw1 = 0.25
 
-Lw2 = 0.6
-Lnw2 = 0.3
+Lw2 = 0.5
+Lnw2 = 0.25
 
-lambda_w = 1
-lambda_nw = 1
+lambda_w = 40
+lambda_nw = 40
 
 include_gravity = True
 debugflag = False
@@ -335,7 +336,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,