setup TPR analysing timesteps

Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm.py

@@ -49,7 +49,7 @@ datestr = date.strftime("%Y-%m-%d")
 
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-R-layered-soil-realistic-g-same-perm-vanG-Mualem"
+use_case = "TP-R-layered-soil-realistic-same-perm-vanG-Mualem-analyise-timesteps"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-layered_soil-g-but-same-perm.py"
 
@@ -78,10 +78,17 @@ resolutions = {
 # that this initial timestep should be counted as.
 # Format:
 # timestep_num: initial_time
-starttimes = {0: 0.0}
+# starttimes = {0: 0.0}
+starttimes = {
+    190: 0.19,
+    440: 0.44,
+    690: 0.69,
+    940: 0.94,
+    1140: 1.14,
+    # 1440: 1.440
+    }
 timestep_size = 0.001
-number_of_timesteps = 1200
-
+number_of_timesteps = 60
 
 # LDD scheme parameters  ######################################################
 Lw1 = 0.025  # /timestep_size
@@ -111,7 +118,7 @@ analyse_condition = False
 plot_timestep_every = 1
 # 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 = 1
+number_of_timesteps_to_analyse = 2
 
 # 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.
@@ -119,14 +126,14 @@ if mesh_study:
     write_to_file = {
         # output the relative errornorm (integration in space) w.r.t. an exact
         # solution for each timestep into a csv file.
-        'space_errornorms': True,
+        'space_errornorms': False,
         # save the mesh and marker functions to disk
-        'meshes_and_markers': True,
+        'meshes_and_markers': False,
         # save xdmf/h5 data for each LDD iteration for timesteps determined by
         # number_of_timesteps_to_analyse. I/O intensive!
         'L_iterations_per_timestep': False,
         # save solution to xdmf/h5.
-        'solutions': True,
+        'solutions': False,
         # save absolute differences w.r.t an exact solution to xdmf/h5 file
         # to monitor where on the domains errors happen
         'absolute_differences': True,
@@ -428,7 +435,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,