restart same intrinsic permeability with flattening parameters

Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic-same-intrinsic-perm.py

@@ -30,7 +30,7 @@ date = datetime.datetime.now()
 datestr = date.strftime("%Y-%m-%d")
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-R-2-patch-realistic-same-intrinsic-perm-different-L"
+use_case = "TP-R-2-patch-realistic-same-intrinsic-perm-flattening-L"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-2-patch-realistic-same-intrinsic-perm.py"
 
@@ -47,7 +47,7 @@ resolutions = {
                 # 4: 1e-5,
                 # 8: 1e-5,
                 # 16: 5e-6,
-                32: 5e-6,
+                32: 3e-6,
                 # 64: 2e-6,
                 # 128: 2e-6,
                 # 256: 1e-6,
@@ -57,21 +57,21 @@ resolutions = {
 # The list is looped over and a simulation is run with t_0 as initial time
 #  for each element t_0 in starttimes.
 starttimes = {0: 0.0}
-timestep_size = 0.01
-number_of_timesteps = 8
+timestep_size = 0.005
+number_of_timesteps = 1000
 
 # LDD scheme parameters  ######################################################
-Lw1 = 0.025
-Lnw1 = 0.025
+Lw1 = 0.01
+Lnw1 = 0.01
 
-Lw2 = 0.025
-Lnw2 = 0.025
+Lw2 = 0.005
+Lnw2 = 0.005
 
-lambda_w = 4
-lambda_nw = 4
+lambda_w = 0.5
+lambda_nw = 0.5
 
 include_gravity = False
-debugflag = True
+debugflag = False
 analyse_condition = False
 
 # I/O CONFIG  #################################################################
@@ -81,7 +81,7 @@ analyse_condition = False
 plot_timestep_every = 4
 # Decide how many timesteps you want analysed. Analysed means, that
 # subsequent errors of the L-iteration within the timestep are written out.
-number_of_timesteps_to_analyse = 5
+number_of_timesteps_to_analyse = 10
 
 # fine grained control over data to be written to disk in the mesh study case
 # as well as for a regular simuation for a fixed grid.
@@ -334,7 +334,7 @@ if __name__ == '__main__':
             #     parameter=simulation_parameter
             #     )
 
-        LDDsim.join()
+        # LDDsim.join()
         if mesh_study:
             simulation_output_dir = processQueue.get()
             hlp.merge_spacetime_errornorms(isRichards=isRichards,