setup analysing solver at times for diff perm type

Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.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-different-intrinsic-vanG-Mualem"
+use_case = "TP-R-layered-soil-different-intrinsic-vanG-Mualem-analyise-timesteps"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-layered_soil.py"
 
@@ -75,9 +75,16 @@ resolutions = {
 # starttimes gives a list of starttimes to run the simulation from.
 # 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}
+starttimes = {
+    190: 0.19,
+    440: 0.44,
+    690: 0.69,
+    940: 0.94,
+    1190: 1.19,
+    1440: 1.440
+    }
 timestep_size = 0.001
-number_of_timesteps = 1500
+number_of_timesteps = 60
 
 
 # LDD scheme parameters  ######################################################
@@ -108,7 +115,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.
@@ -116,17 +123,17 @@ 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,
+        'absolute_differences': False,
         # analyise condition numbers for timesteps determined by
         # number_of_timesteps_to_analyse and save them over time to csv.
         'condition_numbers': analyse_condition,
@@ -432,7 +439,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,