From c94f31a4acaffd1e830b2559f3067bf068ed1f14 Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Fri, 3 May 2019 20:25:57 +0200
Subject: [PATCH] debug calc_gli_term
---
LDDsimulation/domainPatch.py | 4 ++--
RR-2-patch-test-case/RR-2-patch-test.py | 8 ++++----
2 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/LDDsimulation/domainPatch.py b/LDDsimulation/domainPatch.py
index d2e1bb6..3c50a01 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 76d6a07..728b86c 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
--
GitLab