From e6259ed9e73ca679b57a752f132200d75abcfc7a Mon Sep 17 00:00:00 2001 From: David <forenkram@gmx.de> Date: Tue, 30 Jun 2020 11:08:33 +0200 Subject: [PATCH] setup another layered soil TP-R simulation --- ...layered_soil-g-but-same-perm-coarse-dt-longterm.py | 11 +++++------ 1 file changed, 5 insertions(+), 6 deletions(-) diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py index 269cfbc..0612f28 100755 --- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py +++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm-coarse-dt-longterm.py @@ -63,13 +63,13 @@ number_of_timesteps = 400 # LDD scheme parameters ###################################################### -Lw1 = 0.025 # /timestep_size +Lw1 = 0.05 # /timestep_size Lnw1 = Lw1 -Lw2 = 0.025 # /timestep_size +Lw2 = 0.05 # /timestep_size Lnw2 = Lw2 -Lw3 = 0.025 # /timestep_size +Lw3 = 0.05 # /timestep_size Lnw3 = Lw3 -Lw4 = 0.025 # /timestep_size +Lw4 = 0.05 # /timestep_size Lnw4 = Lw4 lambda12_w = 4 @@ -274,10 +274,9 @@ t = sym.symbols('t', positive=True) p_e_sym_2patch = { 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), '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 = { -- GitLab