From 18119b0b2a986120e53ef616a34ec3be0b81198a Mon Sep 17 00:00:00 2001
From: David Seus <david.seus@ians.uni-stuttgart.de>
Date: Fri, 16 Aug 2019 17:36:00 +0200
Subject: [PATCH] adjust examples
---
.../TP-TP-layered_soil_with_inner_patch.py | 4 ++--
TP-TP-layered-soil-case/TP-TP-layered_soil.py | 12 ++++++------
2 files changed, 8 insertions(+), 8 deletions(-)
diff --git a/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py b/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
index 1e5e0a6..b46fbab 100755
--- a/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
+++ b/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
@@ -25,7 +25,7 @@ use_case = "TP-TP-layered-soil-with-inner-patch"
solver_tol = 1E-6
############ GRID #######################ü
-mesh_resolution = 50
+mesh_resolution = 30
timestep_size = 0.001
number_of_timesteps = 1500
# decide how many timesteps you want analysed. Analysed means, that we write out
@@ -43,7 +43,7 @@ include_gravity = True
debugflag = False
analyse_condition = False
-output_string = "./output/nondirichlet_number_of_timesteps{}_".format(number_of_timesteps)
+output_string = "./output/test-after-bugfix-nondirichlet_number_of_timesteps{}_".format(number_of_timesteps)
# global domain
subdomain0_vertices = [df.Point(-1.0,-1.0), #
diff --git a/TP-TP-layered-soil-case/TP-TP-layered_soil.py b/TP-TP-layered-soil-case/TP-TP-layered_soil.py
index a56a37a..4eb4c59 100755
--- a/TP-TP-layered-soil-case/TP-TP-layered_soil.py
+++ b/TP-TP-layered-soil-case/TP-TP-layered_soil.py
@@ -22,7 +22,7 @@ import helpers as hlp
sym.init_printing()
use_case="TP-TP-layered-soil"
-solver_tol = 5E-7
+solver_tol = 1E-6
############ GRID #######################ü
mesh_resolution = 30
@@ -43,7 +43,7 @@ include_gravity = True
debugflag = True
analyse_condition = False
-output_string = "./output/number_of_timesteps{}_".format(number_of_timesteps)
+output_string = "./output/test-after-bugfix-number_of_timesteps{}_".format(number_of_timesteps)
# global domain
subdomain0_vertices = [df.Point(-1.0,-1.0), #
@@ -406,10 +406,10 @@ x, y = sym.symbols('x[0], x[1]') # needed by UFL
t = sym.symbols('t', positive=True)
p_e_sym_2patch = {
- 1: {'wetting': -3 - (1+t*t)*(1 + x*x + y*y),
- 'nonwetting': -1-t*(1-y - x**2)**2 - sym.sqrt(2+t**2)*(1-y)**2},
- 2: {'wetting': -3.0 - (1.0 + t*t)*(1.0 + x*x),
- 'nonwetting': -1-t*(1- x**2)**2 - sym.sqrt(2+t**2)*(1-y)**2},
+ 1: {'wetting': -7 - (1+t*t)*(1 + x*x + y*y),
+ 'nonwetting': -1-t*(1.1 + y + x**2)**2},
+ 2: {'wetting': -7.0 - (1.0 + t*t)*(1.0 + x*x),
+ 'nonwetting': -1-t*(1.1 + x**2)**2 - sym.sqrt(2+t**2)*y**2},
}
p_e_sym = {
--
GitLab