diff --git a/Rechnungen_auszuwerten.txt b/Rechnungen_auszuwerten.txt
index b4fed71710b63af688b1ae7e950cd39e80003661..486461d463497af067a516049379a5f1f5d5ba53 100644
--- a/Rechnungen_auszuwerten.txt
+++ b/Rechnungen_auszuwerten.txt
@@ -2,6 +2,8 @@ Rechungen, die noch nicht ausgewertet wurden
 
 - TP-one-patch (purely positive pc)
 - TP-TP-layered-soil-case-with-inner-patch
+- TP-R-two-patch-test-case
+
 
 Rechnungen, die neu gestartet werden müssen:
     -TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting
diff --git a/TP-R-two-patch-test-case/TP-R-2-patch-test.py b/TP-R-two-patch-test-case/TP-R-2-patch-test.py
index f0710f0d532f03c1d7cdda0f1505c04ae6b66e79..168870d1569f1ce255675f364d4aa152b23a2592 100755
--- a/TP-R-two-patch-test-case/TP-R-2-patch-test.py
+++ b/TP-R-two-patch-test-case/TP-R-2-patch-test.py
@@ -18,24 +18,24 @@ solver_tol = 5E-7
 
 ############ GRID #######################ü
 mesh_resolution = 40
-timestep_size = 0.00001
+timestep_size = 0.000001
 number_of_timesteps = 15
 # 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 = 10
 starttime = 0
 
-Lw = 2.5  #/timestep_size
+Lw = 5 #/timestep_size
 Lnw=Lw
 
-l_param_w = 40
-l_param_nw = 40
+l_param_w = 100
+l_param_nw = 100
 
 include_gravity = True
 debugflag = True
 analyse_condition = False
 
-output_string = "./output/before_reimplementing_gravity_term_number_of_timesteps{}_".format(number_of_timesteps)
+output_string = "./output/after_reimplementing_gravity_term_number_of_timesteps{}_".format(number_of_timesteps)
 
 ##### Domain and Interface ####
 # global simulation domain domain
@@ -325,7 +325,7 @@ t = sym.symbols('t', positive=True)
 p_e_sym = {
     1: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x + y*y))}, #*(1-x)**2*(1+x)**2*(1-y)**2},
     2: {'wetting': (-5.0 - (1.0 + t*t)*(1.0 + x*x)), #*(1-x)**2*(1+x)**2*(1+y)**2,
-        'nonwetting': (-1-t*(1.1+y + x**2) - sym.sqrt(2+t**4))*y**2}, #*(1-x)**2*(1+x)**2*(1+y)**2},
+        'nonwetting': (-1-t*(1.1+y + x**2))*y**2}, #*(1-x)**2*(1+x)**2*(1+y)**2},
 } #-y*y*(sym.sin(-2*t+2*x)*sym.sin(1/2*y-1.2*t)) - t*t*x*(0.5-y)*y*(1-x)