diff --git a/LDDsimulation/domainPatch.py b/LDDsimulation/domainPatch.py
index d2e1bb66188f4bf0cc492b0567a71bdad21c398d..3c50a01076c6957c17d29712f514e3f15f7955ab 100644
--- a/LDDsimulation/domainPatch.py
+++ b/LDDsimulation/domainPatch.py
@@ -1,4 +1,4 @@
-[interf_ind]"""Class definitions for subdomains representing soil patches for LDD-TPR simulation
+"""Class definitions for subdomains representing soil patches for LDD-TPR simulation
 
 This file contains class definitions for the class DomainPatch which defines the
 subdomains used for the LDD simulation. Each subdomain of this class holds
@@ -581,7 +581,7 @@ class DomainPatch(df.SubDomain):
                     phi_a  = self.testfunction[phase]
                     gli_pressure_part = df.assemble(-2*dt*Lambda[phase]*pa_prev*phi_a*ds)
                     gli_p_prev_part = gli_pressure_part.get_local()
-                    gli = gli_p_prev_part - gli_prev_neighbour
+                    gli = gli_p_prev_part - gli_prev_of_neighbour
                     # check if there are shared dofs with another interface.
                     if dofs_in_common_with_other_interfaces[phase].size != 0:
                         for common_dof in dofs_in_common_with_other_interfaces[phase]:
diff --git a/RR-2-patch-test-case/RR-2-patch-test.py b/RR-2-patch-test-case/RR-2-patch-test.py
index 76d6a07035837ca394733cfffcb96673ce85df74..728b86c9fe7242070aca7201bf39962677bda6b6 100755
--- a/RR-2-patch-test-case/RR-2-patch-test.py
+++ b/RR-2-patch-test-case/RR-2-patch-test.py
@@ -85,10 +85,10 @@ isRichards = {
     }
 
 ##################### MESH #####################################################
-mesh_resolution = 40
+mesh_resolution = 6
 ##################### TIME #####################################################
 timestep_size = 4*0.0001
-number_of_timesteps = 5000
+number_of_timesteps = 200
 # 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 = 11
@@ -115,8 +115,8 @@ L = {#
 
 lambda_param = {#
 # subdom_num : lambda parameter for the L-scheme
-    1 : {'wetting' :100},#
-    2 : {'wetting' :100}
+    1 : {'wetting' :4},#
+    2 : {'wetting' :4}
 }
 
 ## relative permeabilty functions on subdomain 1