diff --git a/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py b/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py
index 84614ce80a79dbea9bfadca85fb709e9aed828c4..4178c0f50ccd8ffd3aa26fa94734003e42640509 100755
--- a/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py
+++ b/Two-phase-Two-phase/two-patch/TP-TP-2-patch-test-case/TP-TP-2-patch-same-intrinsic-perm.py
@@ -36,7 +36,7 @@ thisfile = "TP-TP-2-patch-same-intrinsic-perm.py"
 
 # GENERAL SOLVER CONFIG  ######################################################
 # maximal iteration per timestep
-max_iter_num = 300
+max_iter_num = 400
 FEM_Lagrange_degree = 1
 
 # GRID AND MESH STUDY SPECIFICATIONS  #########################################
@@ -47,7 +47,7 @@ resolutions = {
                 # 4: 1e-6,
                 # 8: 1e-6,
                 # 16: 5e-6,
-                32: 3e-6,
+                32: 1e-5,
                 # 64: 2e-6,
                 # 128: 1e-6,
                 # 256: 1e-6,
@@ -61,16 +61,16 @@ timestep_size = 0.001
 number_of_timesteps = 1000
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.25 #/timestep_size
+Lw1 = 0.5 #/timestep_size
 Lnw1= 0.25
 
-Lw2 = 0.25 #/timestep_size
+Lw2 = 0.5 #/timestep_size
 Lnw2= 0.25
 
 lambda_w = 4
 lambda_nw = 4
 
-include_gravity = False
+include_gravity = True
 debugflag = False
 analyse_condition = False
 
@@ -120,16 +120,9 @@ else:
     }
 
 # OUTPUT FILE STRING  #########################################################
-if mesh_study:
-    output_string = "./output/{}-{}_timesteps{}_P{}".format(
-        datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
-        )
-else:
-    for tol in resolutions.values():
-        solver_tol = tol
-    output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
-        datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
-        )
+output_string = "./output/{}-{}_timesteps{}_P{}".format(
+    datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
+    )
 
 
 # DOMAIN AND INTERFACE  #######################################################
@@ -243,8 +236,8 @@ lambda_param = {#
 }
 
 intrinsic_permeability = {
-    1: 0.1,
-    2: 0.1,
+    1: 0.01,
+    2: 0.01,
 }