Skip to content
Snippets Groups Projects
Commit 18119b0b authored by David Seus's avatar David Seus
Browse files

adjust examples

parent 1ea0cfe0
No related branches found
No related tags found
No related merge requests found
...@@ -25,7 +25,7 @@ use_case = "TP-TP-layered-soil-with-inner-patch" ...@@ -25,7 +25,7 @@ use_case = "TP-TP-layered-soil-with-inner-patch"
solver_tol = 1E-6 solver_tol = 1E-6
############ GRID #######################ü ############ GRID #######################ü
mesh_resolution = 50 mesh_resolution = 30
timestep_size = 0.001 timestep_size = 0.001
number_of_timesteps = 1500 number_of_timesteps = 1500
# decide how many timesteps you want analysed. Analysed means, that we write out # decide how many timesteps you want analysed. Analysed means, that we write out
...@@ -43,7 +43,7 @@ include_gravity = True ...@@ -43,7 +43,7 @@ include_gravity = True
debugflag = False debugflag = False
analyse_condition = 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 # global domain
subdomain0_vertices = [df.Point(-1.0,-1.0), # subdomain0_vertices = [df.Point(-1.0,-1.0), #
......
...@@ -22,7 +22,7 @@ import helpers as hlp ...@@ -22,7 +22,7 @@ import helpers as hlp
sym.init_printing() sym.init_printing()
use_case="TP-TP-layered-soil" use_case="TP-TP-layered-soil"
solver_tol = 5E-7 solver_tol = 1E-6
############ GRID #######################ü ############ GRID #######################ü
mesh_resolution = 30 mesh_resolution = 30
...@@ -43,7 +43,7 @@ include_gravity = True ...@@ -43,7 +43,7 @@ include_gravity = True
debugflag = True debugflag = True
analyse_condition = False 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 # global domain
subdomain0_vertices = [df.Point(-1.0,-1.0), # subdomain0_vertices = [df.Point(-1.0,-1.0), #
...@@ -406,10 +406,10 @@ x, y = sym.symbols('x[0], x[1]') # needed by UFL ...@@ -406,10 +406,10 @@ x, y = sym.symbols('x[0], x[1]') # needed by UFL
t = sym.symbols('t', positive=True) t = sym.symbols('t', positive=True)
p_e_sym_2patch = { p_e_sym_2patch = {
1: {'wetting': -3 - (1+t*t)*(1 + x*x + y*y), 1: {'wetting': -7 - (1+t*t)*(1 + x*x + y*y),
'nonwetting': -1-t*(1-y - x**2)**2 - sym.sqrt(2+t**2)*(1-y)**2}, 'nonwetting': -1-t*(1.1 + y + x**2)**2},
2: {'wetting': -3.0 - (1.0 + t*t)*(1.0 + x*x), 2: {'wetting': -7.0 - (1.0 + t*t)*(1.0 + x*x),
'nonwetting': -1-t*(1- x**2)**2 - sym.sqrt(2+t**2)*(1-y)**2}, 'nonwetting': -1-t*(1.1 + x**2)**2 - sym.sqrt(2+t**2)*y**2},
} }
p_e_sym = { p_e_sym = {
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment