diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py
index dd7584a7848da27eb0a17d6568483c13362712dd..0b8104efdb6d235f21d8be3319a0b19684a69054 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch-realistic.py
@@ -38,18 +38,18 @@ resolutions = {
                 # 8: 1e-6,  # h=0.2814
                 # 16: 5e-7, # h=0.1412
                 # 32: 4e-7,   # h=0.0706
-                64: 7e-7,
+                64: 1e-6,
                 # 128: 5e-7
                 }
 
 ############ GRID #######################
 # mesh_resolution = 20
-timestep_size = 0.005
-number_of_timesteps = 240
+timestep_size = 0.001
+number_of_timesteps = 1100
 plot_timestep_every = 1
 # decide how many timesteps you want analysed. Analysed means, that we write out
 # subsequent errors of the L-iteration within the timestep.
-number_of_timesteps_to_analyse = 8
+number_of_timesteps_to_analyse = 10
 starttime = 0.0
 
 Lw = 0.025 #/timestep_size
@@ -655,13 +655,13 @@ p_e_sym = {
     2: {'wetting': -5.0 - (1.0 + t*t)*(1.0 + x*x + y*y),
         'nonwetting': (-1 -t*(1.1 + y + x**2)) },
     3: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
     4: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
     5: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
     6: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)),
-        'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2) },
+        'nonwetting': (-1 -t*(1 + x**2) - sym.sin(2+t**2)**2*(1+y)*y**2) },
     # 2: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)),
     #     'nonwetting': - 2 - t*(1 + (y-5.0) + x**2)**2 -sym.sqrt(2+t**2)*(1 + (y-5.0))},
     # 3: {'wetting': 1.0 - (1.0 + t*t)*(10.0 + x*x + (y-5.0)*(y-5.0)*3*sym.sin(-2*t+2*x)*sym.sin(1/2*y-1.2*t)),