From 006f9a1e6fe60ee1bf9e60659e40ef9be4a30166 Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Fri, 16 Oct 2020 15:09:21 +0200
Subject: [PATCH] setup TPR analysing timesteps
---
.../TP-R-layered_soil-g-but-same-perm.py | 25 ++++++++++++-------
1 file changed, 16 insertions(+), 9 deletions(-)
diff --git a/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm.py b/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm.py
index c842b6f..e0eee48 100755
--- a/Usecases/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-g-but-same-perm.py
+++ b/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,
--
GitLab