set up realistic example

Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py

@@ -23,16 +23,16 @@ datestr = date.strftime("%Y-%m-%d")
 # init sympy session
 sym.init_printing()
 # solver_tol = 6E-7
-use_case = "TP-R-layered-soil-realistic"
-max_iter_num = 500
+use_case = "TP-R-layered-soil-realistic-new-lambda"
+max_iter_num = 600
 FEM_Lagrange_degree = 1
 mesh_study = False
 resolutions = {
                 # 1: 1e-7,  # h=2
                 # 2: 2e-5,  # h=1.1180
                 # 4: 1e-6,  # h=0.5590
-                8: 1e-6,  # h=0.2814
-                # 16: 1e-6, # h=0.1412
+                # 8: 1e-6,  # h=0.2814
+                16: 1e-6, # h=0.1412
                 # 32: 1e-6,
                 # 64: 5e-7,
                 # 128: 5e-7
@@ -40,23 +40,29 @@ resolutions = {
 
 # GRID #######################
 # mesh_resolution = 20
-timestep_size = 0.001
-number_of_timesteps = 5
-plot_timestep_every = 1
+timestep_size = 0.005
+number_of_timesteps = 500
+plot_timestep_every = 2
 # 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 = 2
+number_of_timesteps_to_analyse = 5
 starttimes = [0.0]
 
-Lw = 0.025  # /timestep_size
-Lnw = Lw
+Lw1 = 0.25  # /timestep_size
+Lnw1 = Lw1
+Lw2 = 0.025  # /timestep_size
+Lnw2 = Lw2
+Lw3 = 0.025  # /timestep_size
+Lnw3 = Lw3
+Lw4 = 0.025  # /timestep_size
+Lnw4 = Lw4
 
 lambda12_w = 4
 lambda12_nw = 4
-lambda23_w = 4
-lambda23_nw = 4
-lambda34_w = 4
-lambda34_nw = 4
+lambda23_w = 40
+lambda23_nw = 40
+lambda34_w = 40
+lambda34_nw = 40
 
 include_gravity = False
 debugflag = True
@@ -298,14 +304,14 @@ porosity = {
 
 # subdom_num : subdomain L for L-scheme
 L = {
-    1: {'wetting': Lw,
-        'nonwetting': Lnw},
-    2: {'wetting': Lw,
-        'nonwetting': Lnw},
-    3: {'wetting': Lw,
-        'nonwetting': Lnw},
-    4: {'wetting': Lw,
-        'nonwetting': Lnw}
+    1: {'wetting': Lw1,
+        'nonwetting': Lnw1},
+    2: {'wetting': Lw2,
+        'nonwetting': Lnw2},
+    3: {'wetting': Lw3,
+        'nonwetting': Lnw3},
+    4: {'wetting': Lw4,
+        'nonwetting': Lnw4}
 }
 
 # interface_num : lambda parameter for the L-scheme on that interface.