diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py
index 269cfbcfc274b11511b3431823630b3f04893d4f..0612f280dfd66836b1775ddd543444094d270bcd 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py
@@ -63,13 +63,13 @@ number_of_timesteps = 400
 
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.025  # /timestep_size
+Lw1 = 0.05  # /timestep_size
 Lnw1 = Lw1
-Lw2 = 0.025  # /timestep_size
+Lw2 = 0.05  # /timestep_size
 Lnw2 = Lw2
-Lw3 = 0.025  # /timestep_size
+Lw3 = 0.05  # /timestep_size
 Lnw3 = Lw3
-Lw4 = 0.025  # /timestep_size
+Lw4 = 0.05  # /timestep_size
 Lnw4 = Lw4
 
 lambda12_w = 4
@@ -274,10 +274,9 @@ t = sym.symbols('t', positive=True)
 
 p_e_sym_2patch = {
     1: {'wetting': -6 - (1+t*t)*(1 + x*x + y*y),
-        'nonwetting': 0.0*t},  # -1-t*(1.1 + y + x**2)**2},
+        'nonwetting': 0.0*t},
     2: {'wetting': -6.0 - (1.0 + t*t)*(1.0 + x*x),
         'nonwetting': (-1-t*(1.1+y + x**2))*y**2},
-        # 'nonwetting': (-1-t*(1.1 + x**2)**2 - sym.sqrt(5+t**2))*y**2},
 }
 
 p_e_sym = {