Skip to content
Snippets Groups Projects
Commit e6259ed9 authored by David's avatar David
Browse files

setup another layered soil TP-R simulation

parent 4cab014a
No related branches found
No related tags found
No related merge requests found
...@@ -63,13 +63,13 @@ number_of_timesteps = 400 ...@@ -63,13 +63,13 @@ number_of_timesteps = 400
# LDD scheme parameters ###################################################### # LDD scheme parameters ######################################################
Lw1 = 0.025 # /timestep_size Lw1 = 0.05 # /timestep_size
Lnw1 = Lw1 Lnw1 = Lw1
Lw2 = 0.025 # /timestep_size Lw2 = 0.05 # /timestep_size
Lnw2 = Lw2 Lnw2 = Lw2
Lw3 = 0.025 # /timestep_size Lw3 = 0.05 # /timestep_size
Lnw3 = Lw3 Lnw3 = Lw3
Lw4 = 0.025 # /timestep_size Lw4 = 0.05 # /timestep_size
Lnw4 = Lw4 Lnw4 = Lw4
lambda12_w = 4 lambda12_w = 4
...@@ -274,10 +274,9 @@ t = sym.symbols('t', positive=True) ...@@ -274,10 +274,9 @@ t = sym.symbols('t', positive=True)
p_e_sym_2patch = { p_e_sym_2patch = {
1: {'wetting': -6 - (1+t*t)*(1 + x*x + y*y), 1: {'wetting': -6 - (1+t*t)*(1 + x*x + y*y),
'nonwetting': 0.0*t}, # -1-t*(1.1 + y + x**2)**2}, 'nonwetting': 0.0*t},
2: {'wetting': -6.0 - (1.0 + t*t)*(1.0 + x*x), 2: {'wetting': -6.0 - (1.0 + t*t)*(1.0 + x*x),
'nonwetting': (-1-t*(1.1+y + x**2))*y**2}, 'nonwetting': (-1-t*(1.1+y + x**2))*y**2},
# 'nonwetting': (-1-t*(1.1 + x**2)**2 - sym.sqrt(5+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