From 888b0c8fb0ba3baaf7ad6c4cce64c12958a0cb1c Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Fri, 16 Oct 2020 14:45:22 +0200
Subject: [PATCH] setup analysing solver at times for diff perm type
---
.../layered_soil/TP-R-layered_soil.py | 25 ++++++++++++-------
1 file changed, 16 insertions(+), 9 deletions(-)
diff --git a/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py b/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
index c58def0..9aa98b7 100755
--- a/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
+++ b/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,
--
GitLab