setup layered soil case mesh study

Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case/mesh_study/TP-TP-layered_soil_pure_dd_mesh_study.py

@@ -32,35 +32,35 @@ max_iter_num = 1000
 FEM_Lagrange_degree = 1
 mesh_study = True
 resolutions = {
-                # 1: 1e-7,  # h=2
-                2: 1e-7,  # h=1.1180
-                # 4: 1e-7,  # h=0.5590
-                # 8: 1e-7,  # h=0.2814
-                # 16: 1e-7, # h=0.1412
-                # 32: 1e-7,   # h=0.0706
-                # 64: 1e-7,
-                # 128: 1e-7,
+                1: 1e-6,  # h=2
+                2: 1e-6,  # h=1.1180
+                4: 1e-6,  # h=0.5590
+                8: 1e-6,  # h=0.2814
+                16: 1e-6, # h=0.1412
+                32: 1e-6,   # h=0.0706
+                64: 1e-6,
+                128: 1e-6,
                 # 256: 1e-7
                 }
 
 ############ GRID #######################
 # mesh_resolution = 20
-timestep_size = 0.001
-number_of_timesteps = 300
+timestep_size = 0.005
+number_of_timesteps = 100
 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 = 0
-starttimes = [0.0, 0.25, 0.5]
+starttimes = [0.0]
 
 Lw = 0.025 #/timestep_size
 Lnw=Lw
 
-lambda_w = 80
-lambda_nw = 80
+lambda_w = 8
+lambda_nw = 8
 
-include_gravity = True
-debugflag = True
+include_gravity = False
+debugflag = False
 analyse_condition = True
 
 if mesh_study:
@@ -76,10 +76,10 @@ if mesh_study:
         'space_errornorms': True,
         'meshes_and_markers': True,
         'L_iterations_per_timestep': False,
-        'solutions': False,
-        'absolute_differences': False,
+        'solutions': True,
+        'absolute_differences': True,
         'condition_numbers': analyse_condition,
-        'subsequent_errors': False
+        'subsequent_errors': True
     }
 else:
     write_to_file = {
@@ -262,7 +262,7 @@ densities = {
          'nonwetting':1},  #1.225}}
 }
 
-gravity_acceleration = 9.81
+gravity_acceleration = 1
 # porosities taken from
 # https://www.geotechdata.info/parameter/soil-porosity.html
 # Dict of the form: { subdom_num : porosity }