From cde441a29f2034746dfa453e67775fad78eb2189 Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Tue, 23 Jun 2020 16:43:55 +0200
Subject: [PATCH] setup new inner patch TPR simulation

---
 ...layered_soil_with_inner_patch-realistic.py | 16 ++++-----
 .../TP-TP-layered_soil_with_inner_patch.py    | 36 +++++++++----------
 2 files changed, 26 insertions(+), 26 deletions(-)

diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py
index 8d330bf..d02bb1f 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-realistic.py
@@ -31,13 +31,13 @@ date = datetime.datetime.now()
 datestr = date.strftime("%Y-%m-%d")
 
 # Name of the usecase that will be printed during simulation.
-use_case = "TP-R-layered-soil-with-inner-patch-realistic"
+use_case = "TP-R-layered-soil-with-inner-patch-realistic-split-interface45"
 # The name of this very file. Needed for creating log output.
 thisfile = "TP-R-layered_soil_with_inner_patch-realistic.py"
 
 # GENERAL SOLVER CONFIG  ######################################################
 # maximal iteration per timestep
-max_iter_num = 700
+max_iter_num = 1000
 FEM_Lagrange_degree = 1
 
 # GRID AND MESH STUDY SPECIFICATIONS  #########################################
@@ -47,8 +47,8 @@ resolutions = {
                 # 2: 2e-6,  # h=1.1180
                 # 4: 2e-6,  # h=0.5590
                 # 8: 2e-6,  # h=0.2814
-                16: 8e-6, # h=0.1412
-                # 32: 2e-6,
+                # 16: 8e-6, # h=0.1412
+                32: 5e-6,
                 # 64: 2e-6,
                 # 128: 2e-6
                 }
@@ -58,7 +58,7 @@ resolutions = {
 #  for each element t_0 in starttimes.
 starttimes = [0.0]
 timestep_size = 0.01
-number_of_timesteps = 5
+number_of_timesteps = 400
 
 # LDD scheme parameters  ######################################################
 
@@ -107,15 +107,15 @@ lambda46_nw = 4
 lambda56_w = 4
 lambda56_nw = 4
 
-include_gravity = True
-debugflag = True
+include_gravity = False
+debugflag = False
 analyse_condition = False
 
 # I/O CONFIG  #################################################################
 # when number_of_timesteps is high, it might take a long time to write all
 # timesteps to disk. Therefore, you can choose to only write data of every
 # plot_timestep_every timestep to disk.
-plot_timestep_every = 4
+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 = 5
diff --git a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
index 5a4362c..7be6b7a 100755
--- a/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
+++ b/Two-phase-Two-phase/multi-patch/TP-TP-layered-soil-case-with-inner-patch/TP-TP-layered_soil_with_inner_patch.py
@@ -79,32 +79,32 @@ Lnw5 = 0.25
 Lw6 = 0.25  # /timestep_size
 Lnw6 = 0.25
 
-lambda12_w = 4
-lambda12_nw = 4
+lambda12_w = 40
+lambda12_nw = 40
 
-lambda23_w = 4
-lambda23_nw = 4
+lambda23_w = 40
+lambda23_nw = 40
 
-lambda24_w = 4
-lambda24_nw= 4
+lambda24_w = 40
+lambda24_nw= 40
 
-lambda25_w= 4
-lambda25_nw= 4
+lambda25_w= 40
+lambda25_nw= 40
 
-lambda34_w = 4
-lambda34_nw = 4
+lambda34_w = 40
+lambda34_nw = 40
 
-lambda36_w = 4
-lambda36_nw = 4
+lambda36_w = 40
+lambda36_nw = 40
 
-lambda45_w = 4
-lambda45_nw = 4
+lambda45_w = 40
+lambda45_nw = 40
 
-lambda46_w = 4
-lambda46_nw = 4
+lambda46_w = 40
+lambda46_nw = 40
 
-lambda56_w = 4
-lambda56_nw = 4
+lambda56_w = 40
+lambda56_nw = 40
 
 include_gravity = False
 debugflag = True
-- 
GitLab