From 28d6f0893d730a0e5ee14572c152102caf465bbf Mon Sep 17 00:00:00 2001
From: David <forenkram@gmx.de>
Date: Fri, 12 Jun 2020 14:48:34 +0200
Subject: [PATCH] correct missing gravity

---
 LDDsimulation/helpers.py                      |   12 +-
 ...ame-wetting-phase-as-RR-zero-nonwetting.py |    1 +
 .../TP-R-layered_soil-all-params-one.py       |    4 +
 .../layered_soil/TP-R-layered_soil.py         |    3 +
 .../TP-R-layered_soil-mesh-study.py           |    3 +
 ...ith_inner_patch-all-params-one-finer-dt.py | 2014 +----------------
 ...ed_soil_with_inner_patch-all-params-one.py |    3 +
 ...layered_soil_with_inner_patch-realistic.py |    3 +
 ...h_inner_patch-all-params-one-mesh-study.py |    3 +
 ...oil_with_inner_patch-pure-dd-mesh-study.py |    3 +
 .../TP-R-2-patch-pure-dd-realistic.py         |  104 +-
 .../TP-R-2-patch-pure-dd.py                   |  104 +-
 .../TP-R-2-patch-realistic.py                 |  103 +-
 .../TP-R-2-patch-test.py                      |    6 +
 ...2-realistic-parameters-densities-scaled.py |    1 +
 ...2-patch-mesh-study-all-params-one-but-g.py |    1 +
 .../TP-R-2-patch-mesh-study-all-params-one.py |    8 +-
 .../mesh_studies/TP-R-2-patch-mesh-study.py   |    1 +
 18 files changed, 238 insertions(+), 2139 deletions(-)

diff --git a/LDDsimulation/helpers.py b/LDDsimulation/helpers.py
index c73f81b..efcaf68 100644
--- a/LDDsimulation/helpers.py
+++ b/LDDsimulation/helpers.py
@@ -111,20 +111,20 @@ def generate_exact_solution_expressions(
 
             if phase == "nonwetting":
                 # x part of div(flux) for nonwetting
-                div_flux_x = -1*(-1/mu*dka(1-S)*dS*dxpc*dxpa + 1/mu*dxdxpa*ka(1-S))
+                div_flux_x = -1/mu*(-1*dka(1-S)*dS*dxpc*dxpa + dxdxpa*ka(1-S))
                 # y part of div(flux) for nonwetting
                 if include_gravity:
-                    div_flux_y = -1*(-1/mu*dka(1-S)*dS*dypc*(dypa + rho*g) + 1/mu*dydypa*ka(1-S))
+                    div_flux_y = -1/mu*(-1*dka(1-S)*dS*dypc*(dypa + rho*g) + dydypa*ka(1-S))
                 else:
-                    div_flux_y = -1*(-1/mu*dka(1-S)*dS*dypc*dypa + 1/mu*dydypa*ka(1-S))
+                    div_flux_y = -1/mu*(-1*dka(1-S)*dS*dypc*dypa + dydypa*ka(1-S))
             else:
                 # x part of div(flux) for wetting
-                div_flux_x = -1*(1/mu*dka(S)*dS*dxpc*dxpa + 1/mu*dxdxpa*ka(S))
+                div_flux_x = -1/mu*(dka(S)*dS*dxpc*dxpa + dxdxpa*ka(S))
                 # y part of div(flux) for wetting
                 if include_gravity:
-                    div_flux_y = -1*(1/mu*dka(S)*dS*dypc*(dypa + rho*g) + 1/mu*dydypa*ka(S))
+                    div_flux_y = -1/mu*(dka(S)*dS*dypc*(dypa + rho*g) + dydypa*ka(S))
                 else:
-                    div_flux_y = -1*(1/mu*dka(S)*dS*dypc*(dypa) + 1/mu*dydypa*ka(S))
+                    div_flux_y = -1/mu*(dka(S)*dS*dypc*(dypa) + dydypa*ka(S))
 
             div_flux[subdomain].update({phase: div_flux_x + div_flux_y})
             contructed_rhs = dtS[subdomain][phase] + div_flux[subdomain][phase]
diff --git a/Two-phase-Richards/multi-patch/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting.py b/Two-phase-Richards/multi-patch/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting.py
index ccc8614..6855e6b 100755
--- a/Two-phase-Richards/multi-patch/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting.py
+++ b/Two-phase-Richards/multi-patch/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting/TP-R-multi-patch-same-wetting-phase-as-RR-zero-nonwetting.py
@@ -477,6 +477,7 @@ simulation.set_parameters(use_case=use_case,
                           exact_solution=exact_solution,
                           densities=densities,
                           include_gravity=include_gravity,
+                          gravity_acceleration=gravity_acceleration,
                           write2file=write_to_file,
                           )
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py
index 60f94ca..468ebb8 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil-all-params-one.py
@@ -612,12 +612,15 @@ for subdomain in isRichards.keys():
                 {outer_boundary_ind: exact_solution[subdomain]}
                 )
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
 print(f.read())
 f.close()
 
+
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -656,6 +659,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
index 8f3f087..4392f91 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/TP-R-layered_soil.py
@@ -720,6 +720,7 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
@@ -727,6 +728,7 @@ print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -765,6 +767,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil-mesh-study.py b/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil-mesh-study.py
index 0e9c966..bbb024b 100755
--- a/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil-mesh-study.py
+++ b/Two-phase-Richards/multi-patch/layered_soil/mesh_study/TP-R-layered_soil-mesh-study.py
@@ -720,6 +720,7 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
@@ -727,6 +728,7 @@ print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -765,6 +767,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one-finer-dt.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one-finer-dt.py
index c4f9167..ddd85b3 100644
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one-finer-dt.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one-finer-dt.py
@@ -1,4 +1,4 @@
-Directory  ./output  created 
+Directory  ./output  created
  type of dtpc is <class 'sympy.core.mul.Mul'>
  type of dS is Piecewise
  type of dtpc is <class 'sympy.core.mul.Mul'>
@@ -924,6 +924,7 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
@@ -931,6 +932,7 @@ print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -969,6 +971,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
@@ -1014,2012 +1017,3 @@ for starttime in starttimes:
                         encoding='utf-8',
                         index=False
                         )
-
-
- ###############################################
- ############# Begin Simulation ################
- ###############################################
-
-Init LDD simulation...
-Set internal fenics parameter...
-Set default parameter...
-
-The following timesteps will be used for plotting:
-initialise meshes and marker functions ...
-
-initialise interfaces ...
-
-Iterating over subset, found 98 entities out of 24359.
-Iterating over subset, found 48 entities out of 24359.
-Iterating over subset, found 32 entities out of 24359.
-Iterating over subset, found 16 entities out of 24359.
-Iterating over subset, found 32 entities out of 24359.
-Iterating over subset, found 16 entities out of 24359.
-Iterating over subset, found 24 entities out of 24359.
-Iterating over subset, found 39 entities out of 24359.
-Iterating over subset, found 40 entities out of 24359.
-initialise subdomains ...
-
-Iterating over subset, found 98 entities out of 2636.
-Iterating over subset, found 98 entities out of 2636.
-Iterating over subset, found 98 entities out of 2636.
-Iterating over subset, found 98 entities out of 2636.
-Iterating over subset, found 98 entities out of 2636.
-Iterating over subset, found 98 entities out of 2636.
-Iterating over subset, found 98 entities out of 2636.
-Iterating over subset, found 98 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 98 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 98 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 98 entities out of 9726.
-Iterating over subset, found 98 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 98 entities out of 9726.
-Iterating over subset, found 98 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 48 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 32 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 16 entities out of 9726.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 48 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 32 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 16 entities out of 2149.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 32 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 24 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 39 entities out of 1637.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 16 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 24 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 40 entities out of 1883.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 16 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 39 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-Iterating over subset, found 40 entities out of 6673.
-creating xdmf files for each subdomain ...
-
-create initial values ...
-create source Expressions ...
-writing exact solutions and/or source terms for all time steps to xdmf files ...
-Start time stepping
-The simulation interval is [0.0, 1.0]
-
-the following timesteps will be analysed: [   1 1000 2000 3000 4000]
-
-
-will save output of this simulation to path ./output/2020-06-11-TP-R-layered-soil-with-inner-patch-all-params-one_timesteps4000_P1mesh_res64_dt0.00025solver_tol2e-06_t0-0.0/
-on subdomain1 I have boundary part 0
-on subdomain2 I have boundary part 0
-on subdomain2 I have boundary part 1
-on subdomain3 I have boundary part 0
-on subdomain5 I have boundary part 0
-on subdomain6 I have boundary part 0
-initialising exact solutions for subdomain1 and phase wetting
-initialising exact solutions for subdomain2 and phase wetting
-initialising exact solutions for subdomain3 and phase wetting
-initialising exact solutions for subdomain3 and phase nonwetting
-initialising exact solutions for subdomain4 and phase wetting
-initialising exact solutions for subdomain4 and phase nonwetting
-initialising exact solutions for subdomain5 and phase wetting
-initialising exact solutions for subdomain5 and phase nonwetting
-initialising exact solutions for subdomain6 and phase wetting
-initialising exact solutions for subdomain6 and phase nonwetting
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  1 at time t = 0.000000
-analyising timestep ...
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-Calling FFC just-in-time (JIT) compiler, this may take some time.
-all phases on all subdomains reached a subsequent error = 1.918211520822758e-06 < tol = 2e-06 after 17 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  2 at time t = 0.000250
-all phases on all subdomains reached a subsequent error = 1.979818869329363e-06 < tol = 2e-06 after 21 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  3 at time t = 0.000500
-all phases on all subdomains reached a subsequent error = 1.9948422135833415e-06 < tol = 2e-06 after 23 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  4 at time t = 0.000750
-all phases on all subdomains reached a subsequent error = 1.9862731302447824e-06 < tol = 2e-06 after 27 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  5 at time t = 0.001000
-all phases on all subdomains reached a subsequent error = 1.9865421015146577e-06 < tol = 2e-06 after 29 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  6 at time t = 0.001250
-all phases on all subdomains reached a subsequent error = 1.988467041204908e-06 < tol = 2e-06 after 29 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  7 at time t = 0.001500
-all phases on all subdomains reached a subsequent error = 1.9922923313248306e-06 < tol = 2e-06 after 31 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  8 at time t = 0.001750
-all phases on all subdomains reached a subsequent error = 1.9958250875920524e-06 < tol = 2e-06 after 33 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  9 at time t = 0.002000
-all phases on all subdomains reached a subsequent error = 1.99870089546709e-06 < tol = 2e-06 after 33 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  10 at time t = 0.002250
-all phases on all subdomains reached a subsequent error = 1.9998257300672894e-06 < tol = 2e-06 after 47 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  11 at time t = 0.002500
-all phases on all subdomains reached a subsequent error = 1.999981304639206e-06 < tol = 2e-06 after 47 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  12 at time t = 0.002750
-all phases on all subdomains reached a subsequent error = 1.9990873260943718e-06 < tol = 2e-06 after 29 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  13 at time t = 0.003000
-all phases on all subdomains reached a subsequent error = 1.999915365433545e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  14 at time t = 0.003250
-all phases on all subdomains reached a subsequent error = 1.9973080053977883e-06 < tol = 2e-06 after 27 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  15 at time t = 0.003500
-all phases on all subdomains reached a subsequent error = 1.999752558056508e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  16 at time t = 0.003750
-all phases on all subdomains reached a subsequent error = 1.997442541709269e-06 < tol = 2e-06 after 27 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  17 at time t = 0.004000
-all phases on all subdomains reached a subsequent error = 1.9995915883209195e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  18 at time t = 0.004250
-all phases on all subdomains reached a subsequent error = 1.9975845775074303e-06 < tol = 2e-06 after 27 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  19 at time t = 0.004500
-all phases on all subdomains reached a subsequent error = 1.9999018447902792e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  20 at time t = 0.004750
-all phases on all subdomains reached a subsequent error = 1.999621591031126e-06 < tol = 2e-06 after 29 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  21 at time t = 0.005000
-all phases on all subdomains reached a subsequent error = 1.9999444465398397e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  22 at time t = 0.005250
-all phases on all subdomains reached a subsequent error = 1.9998073438200985e-06 < tol = 2e-06 after 29 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  23 at time t = 0.005500
-all phases on all subdomains reached a subsequent error = 1.999634615238944e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  24 at time t = 0.005750
-all phases on all subdomains reached a subsequent error = 1.9997415432443886e-06 < tol = 2e-06 after 29 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  25 at time t = 0.006000
-all phases on all subdomains reached a subsequent error = 1.9999790280514205e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  26 at time t = 0.006250
-all phases on all subdomains reached a subsequent error = 1.9997382527972418e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  27 at time t = 0.006500
-all phases on all subdomains reached a subsequent error = 1.9998550430676235e-06 < tol = 2e-06 after 29 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  28 at time t = 0.006750
-all phases on all subdomains reached a subsequent error = 1.9997351087279943e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  29 at time t = 0.007000
-all phases on all subdomains reached a subsequent error = 1.9999830591665673e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  30 at time t = 0.007250
-all phases on all subdomains reached a subsequent error = 1.9998629552902427e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  31 at time t = 0.007500
-all phases on all subdomains reached a subsequent error = 1.9996829369415055e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  32 at time t = 0.007750
-all phases on all subdomains reached a subsequent error = 1.9999358659121504e-06 < tol = 2e-06 after 47 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  33 at time t = 0.008000
-all phases on all subdomains reached a subsequent error = 1.999717421720876e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  34 at time t = 0.008250
-all phases on all subdomains reached a subsequent error = 1.9995981136695903e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  35 at time t = 0.008500
-all phases on all subdomains reached a subsequent error = 1.999936513598357e-06 < tol = 2e-06 after 47 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  36 at time t = 0.008750
-all phases on all subdomains reached a subsequent error = 1.9997811208357047e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  37 at time t = 0.009000
-all phases on all subdomains reached a subsequent error = 1.9997088287184935e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  38 at time t = 0.009250
-all phases on all subdomains reached a subsequent error = 1.9996781471329576e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  39 at time t = 0.009500
-all phases on all subdomains reached a subsequent error = 1.999661860678882e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  40 at time t = 0.009750
-all phases on all subdomains reached a subsequent error = 1.9996561835448854e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  41 at time t = 0.010000
-all phases on all subdomains reached a subsequent error = 1.9996608773638082e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  42 at time t = 0.010250
-all phases on all subdomains reached a subsequent error = 1.99967556308115e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  43 at time t = 0.010500
-all phases on all subdomains reached a subsequent error = 1.999699453049972e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  44 at time t = 0.010750
-all phases on all subdomains reached a subsequent error = 1.9997315267865606e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  45 at time t = 0.011000
-all phases on all subdomains reached a subsequent error = 1.9997707098243382e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  46 at time t = 0.011250
-all phases on all subdomains reached a subsequent error = 1.9998159850578517e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  47 at time t = 0.011500
-all phases on all subdomains reached a subsequent error = 1.99986644753625e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  48 at time t = 0.011750
-all phases on all subdomains reached a subsequent error = 1.9999213210334746e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  49 at time t = 0.012000
-all phases on all subdomains reached a subsequent error = 1.99997995564464e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  50 at time t = 0.012250
-all phases on all subdomains reached a subsequent error = 1.999618170399563e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  51 at time t = 0.012500
-all phases on all subdomains reached a subsequent error = 1.9998157320666863e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  52 at time t = 0.012750
-all phases on all subdomains reached a subsequent error = 1.9999306578395437e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  53 at time t = 0.013000
-all phases on all subdomains reached a subsequent error = 1.999596620722835e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  54 at time t = 0.013250
-all phases on all subdomains reached a subsequent error = 1.9998242486928033e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  55 at time t = 0.013500
-all phases on all subdomains reached a subsequent error = 1.999967751396954e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  56 at time t = 0.013750
-all phases on all subdomains reached a subsequent error = 1.9996569161546213e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  57 at time t = 0.014000
-all phases on all subdomains reached a subsequent error = 1.9999115573093767e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  58 at time t = 0.014250
-all phases on all subdomains reached a subsequent error = 1.9996459644607683e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  59 at time t = 0.014500
-all phases on all subdomains reached a subsequent error = 1.999922963903211e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  60 at time t = 0.014750
-all phases on all subdomains reached a subsequent error = 1.999672609950191e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  61 at time t = 0.015000
-all phases on all subdomains reached a subsequent error = 1.999967978199895e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  62 at time t = 0.015250
-all phases on all subdomains reached a subsequent error = 1.999731053210406e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  63 at time t = 0.015500
-all phases on all subdomains reached a subsequent error = 1.999609507359447e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  64 at time t = 0.015750
-all phases on all subdomains reached a subsequent error = 1.999961230752958e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  65 at time t = 0.016000
-all phases on all subdomains reached a subsequent error = 1.9997528589456998e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  66 at time t = 0.016250
-all phases on all subdomains reached a subsequent error = 1.9996573015902194e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  67 at time t = 0.016500
-all phases on all subdomains reached a subsequent error = 1.9996010410906844e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  68 at time t = 0.016750
-all phases on all subdomains reached a subsequent error = 1.999998725608488e-06 < tol = 2e-06 after 49 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  69 at time t = 0.017000
-all phases on all subdomains reached a subsequent error = 1.9998238375265143e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  70 at time t = 0.017250
-all phases on all subdomains reached a subsequent error = 1.9997618784840593e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  71 at time t = 0.017500
-all phases on all subdomains reached a subsequent error = 1.9997368501276124e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  72 at time t = 0.017750
-all phases on all subdomains reached a subsequent error = 1.99972500607984e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  73 at time t = 0.018000
-all phases on all subdomains reached a subsequent error = 1.9997233337240725e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  74 at time t = 0.018250
-all phases on all subdomains reached a subsequent error = 1.9997318989630227e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  75 at time t = 0.018500
-all phases on all subdomains reached a subsequent error = 1.9997505134974254e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  76 at time t = 0.018750
-all phases on all subdomains reached a subsequent error = 1.9997785401397994e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  77 at time t = 0.019000
-all phases on all subdomains reached a subsequent error = 1.9998150913284465e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  78 at time t = 0.019250
-all phases on all subdomains reached a subsequent error = 1.999859210872033e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  79 at time t = 0.019500
-all phases on all subdomains reached a subsequent error = 1.999909983203223e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  80 at time t = 0.019750
-all phases on all subdomains reached a subsequent error = 1.999966585673252e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  81 at time t = 0.020000
-all phases on all subdomains reached a subsequent error = 1.9995781986492383e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  82 at time t = 0.020250
-all phases on all subdomains reached a subsequent error = 1.9998034691196397e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  83 at time t = 0.020500
-all phases on all subdomains reached a subsequent error = 1.9999214680981326e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  84 at time t = 0.020750
-all phases on all subdomains reached a subsequent error = 1.9995667031003187e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  85 at time t = 0.021000
-all phases on all subdomains reached a subsequent error = 1.9998280164667358e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  86 at time t = 0.021250
-all phases on all subdomains reached a subsequent error = 1.999979406015239e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  87 at time t = 0.021500
-all phases on all subdomains reached a subsequent error = 1.9996516370086485e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  88 at time t = 0.021750
-all phases on all subdomains reached a subsequent error = 1.9999436622700264e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  89 at time t = 0.022000
-all phases on all subdomains reached a subsequent error = 1.9996621445958534e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  90 at time t = 0.022250
-all phases on all subdomains reached a subsequent error = 1.999978355994722e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  91 at time t = 0.022500
-all phases on all subdomains reached a subsequent error = 1.999713786985629e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  92 at time t = 0.022750
-all phases on all subdomains reached a subsequent error = 1.9995928841999657e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  93 at time t = 0.023000
-all phases on all subdomains reached a subsequent error = 1.999970861257241e-06 < tol = 2e-06 after 51 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  94 at time t = 0.023250
-all phases on all subdomains reached a subsequent error = 1.999741058938296e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  95 at time t = 0.023500
-all phases on all subdomains reached a subsequent error = 1.9996526185826705e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  96 at time t = 0.023750
-all phases on all subdomains reached a subsequent error = 1.9996045791577748e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  97 at time t = 0.024000
-all phases on all subdomains reached a subsequent error = 1.9995740971779432e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  98 at time t = 0.024250
-all phases on all subdomains reached a subsequent error = 1.9995584146170313e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  99 at time t = 0.024500
-all phases on all subdomains reached a subsequent error = 1.999557173392789e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  100 at time t = 0.024750
-all phases on all subdomains reached a subsequent error = 1.9995695574324034e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  101 at time t = 0.025000
-all phases on all subdomains reached a subsequent error = 1.9995943042138673e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  102 at time t = 0.025250
-all phases on all subdomains reached a subsequent error = 1.999629978988537e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  103 at time t = 0.025500
-all phases on all subdomains reached a subsequent error = 1.9996751732235406e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  104 at time t = 0.025750
-all phases on all subdomains reached a subsequent error = 1.9997286075530424e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  105 at time t = 0.026000
-all phases on all subdomains reached a subsequent error = 1.999789168592374e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  106 at time t = 0.026250
-all phases on all subdomains reached a subsequent error = 1.999855910618377e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  107 at time t = 0.026500
-all phases on all subdomains reached a subsequent error = 1.9999280408119495e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  108 at time t = 0.026750
-all phases on all subdomains reached a subsequent error = 1.999540512719676e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  109 at time t = 0.027000
-all phases on all subdomains reached a subsequent error = 1.9997999252415238e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  110 at time t = 0.027250
-all phases on all subdomains reached a subsequent error = 1.9999316363266304e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  111 at time t = 0.027500
-all phases on all subdomains reached a subsequent error = 1.99957692311099e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  112 at time t = 0.027750
-all phases on all subdomains reached a subsequent error = 1.9998722832072025e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  113 at time t = 0.028000
-all phases on all subdomains reached a subsequent error = 1.9995687465105934e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  114 at time t = 0.028250
-all phases on all subdomains reached a subsequent error = 1.9998935571561708e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  115 at time t = 0.028500
-all phases on all subdomains reached a subsequent error = 1.9996126116117475e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  116 at time t = 0.028750
-all phases on all subdomains reached a subsequent error = 1.9999622130530517e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  117 at time t = 0.029000
-all phases on all subdomains reached a subsequent error = 1.9997000325781756e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  118 at time t = 0.029250
-all phases on all subdomains reached a subsequent error = 1.999600612615121e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  119 at time t = 0.029500
-all phases on all subdomains reached a subsequent error = 1.999541998587699e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  120 at time t = 0.029750
-all phases on all subdomains reached a subsequent error = 1.999972333922848e-06 < tol = 2e-06 after 53 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  121 at time t = 0.030000
-all phases on all subdomains reached a subsequent error = 1.9997598173166918e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  122 at time t = 0.030250
-all phases on all subdomains reached a subsequent error = 1.9997079725122286e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  123 at time t = 0.030500
-all phases on all subdomains reached a subsequent error = 1.999694082124124e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  124 at time t = 0.030750
-all phases on all subdomains reached a subsequent error = 1.999695040249009e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  125 at time t = 0.031000
-all phases on all subdomains reached a subsequent error = 1.9997085016849328e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  126 at time t = 0.031250
-all phases on all subdomains reached a subsequent error = 1.9997342918929782e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  127 at time t = 0.031500
-all phases on all subdomains reached a subsequent error = 1.9997717530533028e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  128 at time t = 0.031750
-all phases on all subdomains reached a subsequent error = 1.9998198027053527e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  129 at time t = 0.032000
-all phases on all subdomains reached a subsequent error = 1.9998772013250687e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  130 at time t = 0.032250
-all phases on all subdomains reached a subsequent error = 1.999942735550735e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  131 at time t = 0.032500
-all phases on all subdomains reached a subsequent error = 1.99954689861795e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  132 at time t = 0.032750
-all phases on all subdomains reached a subsequent error = 1.9998156529644257e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  133 at time t = 0.033000
-all phases on all subdomains reached a subsequent error = 1.9999448823875604e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  134 at time t = 0.033250
-all phases on all subdomains reached a subsequent error = 1.9995866273923327e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  135 at time t = 0.033500
-all phases on all subdomains reached a subsequent error = 1.9998953875992114e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  136 at time t = 0.033750
-all phases on all subdomains reached a subsequent error = 1.9995889765344648e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  137 at time t = 0.034000
-all phases on all subdomains reached a subsequent error = 1.999927988893059e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  138 at time t = 0.034250
-all phases on all subdomains reached a subsequent error = 1.999645305954545e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  139 at time t = 0.034500
-all phases on all subdomains reached a subsequent error = 1.999539227895406e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  140 at time t = 0.034750
-all phases on all subdomains reached a subsequent error = 1.9999452225470063e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  141 at time t = 0.035000
-all phases on all subdomains reached a subsequent error = 1.999702778716144e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  142 at time t = 0.035250
-all phases on all subdomains reached a subsequent error = 1.9996343193479073e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  143 at time t = 0.035500
-all phases on all subdomains reached a subsequent error = 1.9996048996286966e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  144 at time t = 0.035750
-all phases on all subdomains reached a subsequent error = 1.9995926849598777e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  145 at time t = 0.036000
-all phases on all subdomains reached a subsequent error = 1.9995958368374966e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  146 at time t = 0.036250
-all phases on all subdomains reached a subsequent error = 1.9996139352906315e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  147 at time t = 0.036500
-all phases on all subdomains reached a subsequent error = 1.9996458818665915e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  148 at time t = 0.036750
-all phases on all subdomains reached a subsequent error = 1.9996901596211916e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  149 at time t = 0.037000
-all phases on all subdomains reached a subsequent error = 1.999745161620822e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  150 at time t = 0.037250
-all phases on all subdomains reached a subsequent error = 1.9998093827960617e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  151 at time t = 0.037500
-all phases on all subdomains reached a subsequent error = 1.9998815001561013e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  152 at time t = 0.037750
-all phases on all subdomains reached a subsequent error = 1.9999603895661895e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  153 at time t = 0.038000
-all phases on all subdomains reached a subsequent error = 1.999578043220159e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  154 at time t = 0.038250
-all phases on all subdomains reached a subsequent error = 1.999865205249265e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  155 at time t = 0.038500
-all phases on all subdomains reached a subsequent error = 1.9995355235963594e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  156 at time t = 0.038750
-all phases on all subdomains reached a subsequent error = 1.999855478502922e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  157 at time t = 0.039000
-all phases on all subdomains reached a subsequent error = 1.9995556006634248e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  158 at time t = 0.039250
-all phases on all subdomains reached a subsequent error = 1.999907155465218e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  159 at time t = 0.039500
-all phases on all subdomains reached a subsequent error = 1.999632944147031e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  160 at time t = 0.039750
-all phases on all subdomains reached a subsequent error = 1.99954271896798e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  161 at time t = 0.040000
-all phases on all subdomains reached a subsequent error = 1.999959670107604e-06 < tol = 2e-06 after 57 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  162 at time t = 0.040250
-all phases on all subdomains reached a subsequent error = 1.999723893401076e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  163 at time t = 0.040500
-all phases on all subdomains reached a subsequent error = 1.9996694832534492e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  164 at time t = 0.040750
-all phases on all subdomains reached a subsequent error = 1.9996519664555705e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  165 at time t = 0.041000
-all phases on all subdomains reached a subsequent error = 1.999650501054114e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  166 at time t = 0.041250
-all phases on all subdomains reached a subsequent error = 1.999663885164175e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  167 at time t = 0.041500
-all phases on all subdomains reached a subsequent error = 1.999691788202695e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  168 at time t = 0.041750
-all phases on all subdomains reached a subsequent error = 1.999733084551157e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  169 at time t = 0.042000
-all phases on all subdomains reached a subsequent error = 1.999786232969025e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  170 at time t = 0.042250
-all phases on all subdomains reached a subsequent error = 1.9998496285205893e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  171 at time t = 0.042500
-all phases on all subdomains reached a subsequent error = 1.999921789521836e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  172 at time t = 0.042750
-all phases on all subdomains reached a subsequent error = 1.9995419951359836e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  173 at time t = 0.043000
-all phases on all subdomains reached a subsequent error = 1.9998274399461375e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  174 at time t = 0.043250
-all phases on all subdomains reached a subsequent error = 1.9999600636701487e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  175 at time t = 0.043500
-all phases on all subdomains reached a subsequent error = 1.9996179119131356e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  176 at time t = 0.043750
-all phases on all subdomains reached a subsequent error = 1.999943741120217e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  177 at time t = 0.044000
-all phases on all subdomains reached a subsequent error = 1.999650164391626e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  178 at time t = 0.044250
-all phases on all subdomains reached a subsequent error = 1.9995451882488753e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  179 at time t = 0.044500
-all phases on all subdomains reached a subsequent error = 1.999939776972972e-06 < tol = 2e-06 after 59 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  180 at time t = 0.044750
-all phases on all subdomains reached a subsequent error = 1.999691694351409e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  181 at time t = 0.045000
-all phases on all subdomains reached a subsequent error = 1.9996309384648153e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  182 at time t = 0.045250
-all phases on all subdomains reached a subsequent error = 1.9996075917466506e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  183 at time t = 0.045500
-all phases on all subdomains reached a subsequent error = 1.9996017086641324e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  184 at time t = 0.045750
-all phases on all subdomains reached a subsequent error = 1.9996122668166585e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  185 at time t = 0.046000
-all phases on all subdomains reached a subsequent error = 1.9996386208826395e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  186 at time t = 0.046250
-all phases on all subdomains reached a subsequent error = 1.999679298285692e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  187 at time t = 0.046500
-all phases on all subdomains reached a subsequent error = 1.9997324829684775e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  188 at time t = 0.046750
-all phases on all subdomains reached a subsequent error = 1.999796373144754e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  189 at time t = 0.047000
-all phases on all subdomains reached a subsequent error = 1.999869352243608e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  190 at time t = 0.047250
-all phases on all subdomains reached a subsequent error = 1.9999500412128177e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  191 at time t = 0.047500
-all phases on all subdomains reached a subsequent error = 1.999587193415156e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  192 at time t = 0.047750
-all phases on all subdomains reached a subsequent error = 1.9998795021928555e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  193 at time t = 0.048000
-all phases on all subdomains reached a subsequent error = 1.9995659168519177e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  194 at time t = 0.048250
-all phases on all subdomains reached a subsequent error = 1.999890801758034e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  195 at time t = 0.048500
-all phases on all subdomains reached a subsequent error = 1.9996078119730168e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  196 at time t = 0.048750
-all phases on all subdomains reached a subsequent error = 1.999964249873366e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  197 at time t = 0.049000
-all phases on all subdomains reached a subsequent error = 1.999706718672856e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  198 at time t = 0.049250
-all phases on all subdomains reached a subsequent error = 1.9996375036737568e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  199 at time t = 0.049500
-all phases on all subdomains reached a subsequent error = 1.9996043318277947e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  200 at time t = 0.049750
-all phases on all subdomains reached a subsequent error = 1.9995896125053296e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  201 at time t = 0.050000
-all phases on all subdomains reached a subsequent error = 1.9995931961027965e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  202 at time t = 0.050250
-all phases on all subdomains reached a subsequent error = 1.999614372055212e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  203 at time t = 0.050500
-all phases on all subdomains reached a subsequent error = 1.999651356757948e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  204 at time t = 0.050750
-all phases on all subdomains reached a subsequent error = 1.999702004987821e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  205 at time t = 0.051000
-all phases on all subdomains reached a subsequent error = 1.9997642360075426e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  206 at time t = 0.051250
-all phases on all subdomains reached a subsequent error = 1.9998362170244817e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  207 at time t = 0.051500
-all phases on all subdomains reached a subsequent error = 1.99991640651597e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  208 at time t = 0.051750
-all phases on all subdomains reached a subsequent error = 1.9995661486446614e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  209 at time t = 0.052000
-all phases on all subdomains reached a subsequent error = 1.9998556604208633e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  210 at time t = 0.052250
-all phases on all subdomains reached a subsequent error = 1.9999903745577267e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  211 at time t = 0.052500
-all phases on all subdomains reached a subsequent error = 1.9996758460368895e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  212 at time t = 0.052750
-all phases on all subdomains reached a subsequent error = 1.9995682129884643e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  213 at time t = 0.053000
-all phases on all subdomains reached a subsequent error = 1.999934072903313e-06 < tol = 2e-06 after 63 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  214 at time t = 0.053250
-all phases on all subdomains reached a subsequent error = 1.9996861912560116e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  215 at time t = 0.053500
-all phases on all subdomains reached a subsequent error = 1.999627128718713e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  216 at time t = 0.053750
-all phases on all subdomains reached a subsequent error = 1.9996038152579714e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  217 at time t = 0.054000
-all phases on all subdomains reached a subsequent error = 1.999598808750617e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  218 at time t = 0.054250
-all phases on all subdomains reached a subsequent error = 1.999611733930811e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  219 at time t = 0.054500
-all phases on all subdomains reached a subsequent error = 1.9996415742987686e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  220 at time t = 0.054750
-all phases on all subdomains reached a subsequent error = 1.999686410394373e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  221 at time t = 0.055000
-all phases on all subdomains reached a subsequent error = 1.999744090180209e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  222 at time t = 0.055250
-all phases on all subdomains reached a subsequent error = 1.9998125957211685e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  223 at time t = 0.055500
-all phases on all subdomains reached a subsequent error = 1.9998901828792654e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  224 at time t = 0.055750
-all phases on all subdomains reached a subsequent error = 1.9999754020611916e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  225 at time t = 0.056000
-all phases on all subdomains reached a subsequent error = 1.99964276644852e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  226 at time t = 0.056250
-all phases on all subdomains reached a subsequent error = 1.999932460904878e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  227 at time t = 0.056500
-all phases on all subdomains reached a subsequent error = 1.9996444349323642e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  228 at time t = 0.056750
-all phases on all subdomains reached a subsequent error = 1.999965558630062e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  229 at time t = 0.057000
-all phases on all subdomains reached a subsequent error = 1.9997078703610373e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  230 at time t = 0.057250
-all phases on all subdomains reached a subsequent error = 1.9996365950241453e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  231 at time t = 0.057500
-all phases on all subdomains reached a subsequent error = 1.9996006668725633e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  232 at time t = 0.057750
-all phases on all subdomains reached a subsequent error = 1.99958487043382e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  233 at time t = 0.058000
-all phases on all subdomains reached a subsequent error = 1.9995893243776807e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  234 at time t = 0.058250
-all phases on all subdomains reached a subsequent error = 1.9996127141295597e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  235 at time t = 0.058500
-all phases on all subdomains reached a subsequent error = 1.9996526902296884e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  236 at time t = 0.058750
-all phases on all subdomains reached a subsequent error = 1.9997067156851116e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  237 at time t = 0.059000
-all phases on all subdomains reached a subsequent error = 1.9997724708467888e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  238 at time t = 0.059250
-all phases on all subdomains reached a subsequent error = 1.999847985917955e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  239 at time t = 0.059500
-all phases on all subdomains reached a subsequent error = 1.9999316458175833e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  240 at time t = 0.059750
-all phases on all subdomains reached a subsequent error = 1.999613310413941e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  241 at time t = 0.060000
-all phases on all subdomains reached a subsequent error = 1.9998961595037612e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  242 at time t = 0.060250
-all phases on all subdomains reached a subsequent error = 1.999620823217655e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  243 at time t = 0.060500
-all phases on all subdomains reached a subsequent error = 1.999935092434966e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  244 at time t = 0.060750
-all phases on all subdomains reached a subsequent error = 1.9996903531788075e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  245 at time t = 0.061000
-all phases on all subdomains reached a subsequent error = 1.999627458631204e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  246 at time t = 0.061250
-all phases on all subdomains reached a subsequent error = 1.999597973907292e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  247 at time t = 0.061500
-all phases on all subdomains reached a subsequent error = 1.999995230814462e-06 < tol = 2e-06 after 67 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  248 at time t = 0.061750
-all phases on all subdomains reached a subsequent error = 1.9998137613438176e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  249 at time t = 0.062000
-all phases on all subdomains reached a subsequent error = 1.9998107744339653e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  250 at time t = 0.062250
-all phases on all subdomains reached a subsequent error = 1.9998340234381767e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  251 at time t = 0.062500
-all phases on all subdomains reached a subsequent error = 1.9998687998412288e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  252 at time t = 0.062750
-all phases on all subdomains reached a subsequent error = 1.9999162181233994e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  253 at time t = 0.063000
-all phases on all subdomains reached a subsequent error = 1.9999759957193107e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  254 at time t = 0.063250
-all phases on all subdomains reached a subsequent error = 1.999653550049691e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  255 at time t = 0.063500
-all phases on all subdomains reached a subsequent error = 1.9999136210398168e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  256 at time t = 0.063750
-all phases on all subdomains reached a subsequent error = 1.9996392129170667e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  257 at time t = 0.064000
-all phases on all subdomains reached a subsequent error = 1.999935831286618e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  258 at time t = 0.064250
-all phases on all subdomains reached a subsequent error = 1.999696478105679e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  259 at time t = 0.064500
-all phases on all subdomains reached a subsequent error = 1.999634886586723e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  260 at time t = 0.064750
-all phases on all subdomains reached a subsequent error = 1.9999971015640522e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  261 at time t = 0.065000
-all phases on all subdomains reached a subsequent error = 1.9998042119719585e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  262 at time t = 0.065250
-all phases on all subdomains reached a subsequent error = 1.9997860472633068e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  263 at time t = 0.065500
-all phases on all subdomains reached a subsequent error = 1.9997950412877168e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  264 at time t = 0.065750
-all phases on all subdomains reached a subsequent error = 1.9998180500774463e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  265 at time t = 0.066000
-all phases on all subdomains reached a subsequent error = 1.999856172878552e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  266 at time t = 0.066250
-all phases on all subdomains reached a subsequent error = 1.99990862819898e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  267 at time t = 0.066500
-all phases on all subdomains reached a subsequent error = 1.999973643085192e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  268 at time t = 0.066750
-all phases on all subdomains reached a subsequent error = 1.9996719228577e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  269 at time t = 0.067000
-all phases on all subdomains reached a subsequent error = 1.999929016302054e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  270 at time t = 0.067250
-all phases on all subdomains reached a subsequent error = 1.999672082011331e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  271 at time t = 0.067500
-all phases on all subdomains reached a subsequent error = 1.9999641742471928e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  272 at time t = 0.067750
-all phases on all subdomains reached a subsequent error = 1.9997410831112406e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  273 at time t = 0.068000
-all phases on all subdomains reached a subsequent error = 1.999689383568917e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  274 at time t = 0.068250
-all phases on all subdomains reached a subsequent error = 1.999668192176524e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  275 at time t = 0.068500
-all phases on all subdomains reached a subsequent error = 1.9996661872366067e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  276 at time t = 0.068750
-all phases on all subdomains reached a subsequent error = 1.9996841269483617e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  277 at time t = 0.069000
-all phases on all subdomains reached a subsequent error = 1.9997203109700195e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  278 at time t = 0.069250
-all phases on all subdomains reached a subsequent error = 1.9997720969224895e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  279 at time t = 0.069500
-all phases on all subdomains reached a subsequent error = 1.9998368413762715e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  280 at time t = 0.069750
-all phases on all subdomains reached a subsequent error = 1.9999122443417694e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  281 at time t = 0.070000
-all phases on all subdomains reached a subsequent error = 1.9999964210238816e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  282 at time t = 0.070250
-all phases on all subdomains reached a subsequent error = 1.9997262208290415e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  283 at time t = 0.070500
-all phases on all subdomains reached a subsequent error = 1.9999891414972806e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  284 at time t = 0.070750
-all phases on all subdomains reached a subsequent error = 1.9997567838468343e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  285 at time t = 0.071000
-all phases on all subdomains reached a subsequent error = 1.9996908454275507e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  286 at time t = 0.071250
-all phases on all subdomains reached a subsequent error = 1.999656991300086e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  287 at time t = 0.071500
-all phases on all subdomains reached a subsequent error = 1.9996452422680723e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  288 at time t = 0.071750
-all phases on all subdomains reached a subsequent error = 1.9996560423684903e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  289 at time t = 0.072000
-all phases on all subdomains reached a subsequent error = 1.9996870240585326e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  290 at time t = 0.072250
-all phases on all subdomains reached a subsequent error = 1.9997349864309286e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  291 at time t = 0.072500
-all phases on all subdomains reached a subsequent error = 1.999796879241701e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  292 at time t = 0.072750
-all phases on all subdomains reached a subsequent error = 1.9998701204634652e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  293 at time t = 0.073000
-all phases on all subdomains reached a subsequent error = 1.999952630717778e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  294 at time t = 0.073250
-all phases on all subdomains reached a subsequent error = 1.9996975992018606e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  295 at time t = 0.073500
-all phases on all subdomains reached a subsequent error = 1.9999513711099027e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  296 at time t = 0.073750
-all phases on all subdomains reached a subsequent error = 1.9997329645261487e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  297 at time t = 0.074000
-all phases on all subdomains reached a subsequent error = 1.9996742584806638e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  298 at time t = 0.074250
-all phases on all subdomains reached a subsequent error = 1.9999887015854958e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  299 at time t = 0.074500
-all phases on all subdomains reached a subsequent error = 1.999822281893168e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  300 at time t = 0.074750
-all phases on all subdomains reached a subsequent error = 1.9998135834455234e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  301 at time t = 0.075000
-all phases on all subdomains reached a subsequent error = 1.9998286058594926e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  302 at time t = 0.075250
-all phases on all subdomains reached a subsequent error = 1.999857847695671e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  303 at time t = 0.075500
-all phases on all subdomains reached a subsequent error = 1.9999027259090126e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  304 at time t = 0.075750
-all phases on all subdomains reached a subsequent error = 1.9999619601749817e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  305 at time t = 0.076000
-all phases on all subdomains reached a subsequent error = 1.9997042109328517e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  306 at time t = 0.076250
-all phases on all subdomains reached a subsequent error = 1.999934814648315e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  307 at time t = 0.076500
-all phases on all subdomains reached a subsequent error = 1.999719136264237e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  308 at time t = 0.076750
-all phases on all subdomains reached a subsequent error = 1.999985372609699e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  309 at time t = 0.077000
-all phases on all subdomains reached a subsequent error = 1.9998040263274404e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  310 at time t = 0.077250
-all phases on all subdomains reached a subsequent error = 1.99977369563707e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  311 at time t = 0.077500
-all phases on all subdomains reached a subsequent error = 1.9997689759414882e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  312 at time t = 0.077750
-all phases on all subdomains reached a subsequent error = 1.99978216004034e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  313 at time t = 0.078000
-all phases on all subdomains reached a subsequent error = 1.99981436745637e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  314 at time t = 0.078250
-all phases on all subdomains reached a subsequent error = 1.999863586457804e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  315 at time t = 0.078500
-all phases on all subdomains reached a subsequent error = 1.9999270324051345e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  316 at time t = 0.078750
-all phases on all subdomains reached a subsequent error = 1.999688719109432e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  317 at time t = 0.079000
-all phases on all subdomains reached a subsequent error = 1.9999141305108624e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  318 at time t = 0.079250
-all phases on all subdomains reached a subsequent error = 1.9997150881805906e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  319 at time t = 0.079500
-all phases on all subdomains reached a subsequent error = 1.9999747271706312e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  320 at time t = 0.079750
-all phases on all subdomains reached a subsequent error = 1.999808803602481e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  321 at time t = 0.080000
-all phases on all subdomains reached a subsequent error = 1.9997863052474913e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  322 at time t = 0.080250
-all phases on all subdomains reached a subsequent error = 1.9997876432446134e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  323 at time t = 0.080500
-all phases on all subdomains reached a subsequent error = 1.999806271973168e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  324 at time t = 0.080750
-all phases on all subdomains reached a subsequent error = 1.999843417830195e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  325 at time t = 0.081000
-all phases on all subdomains reached a subsequent error = 1.9998970125267586e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  326 at time t = 0.081250
-all phases on all subdomains reached a subsequent error = 1.9999642736702697e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  327 at time t = 0.081500
-all phases on all subdomains reached a subsequent error = 1.999744216406675e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  328 at time t = 0.081750
-all phases on all subdomains reached a subsequent error = 1.9999647040095988e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  329 at time t = 0.082000
-all phases on all subdomains reached a subsequent error = 1.9997817656757766e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  330 at time t = 0.082250
-all phases on all subdomains reached a subsequent error = 1.999738981800342e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  331 at time t = 0.082500
-all phases on all subdomains reached a subsequent error = 1.9997234960806416e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  332 at time t = 0.082750
-all phases on all subdomains reached a subsequent error = 1.9997292429049815e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  333 at time t = 0.083000
-all phases on all subdomains reached a subsequent error = 1.9997566087329033e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  334 at time t = 0.083250
-all phases on all subdomains reached a subsequent error = 1.9998026584683656e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  335 at time t = 0.083500
-all phases on all subdomains reached a subsequent error = 1.999863971526888e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  336 at time t = 0.083750
-all phases on all subdomains reached a subsequent error = 1.999937531954931e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  337 at time t = 0.084000
-all phases on all subdomains reached a subsequent error = 1.999737483587208e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  338 at time t = 0.084250
-all phases on all subdomains reached a subsequent error = 1.9999541046762252e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  339 at time t = 0.084500
-all phases on all subdomains reached a subsequent error = 1.999788042278747e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  340 at time t = 0.084750
-all phases on all subdomains reached a subsequent error = 1.999754465155324e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  341 at time t = 0.085000
-all phases on all subdomains reached a subsequent error = 1.999746373995622e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  342 at time t = 0.085250
-all phases on all subdomains reached a subsequent error = 1.999758702588046e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  343 at time t = 0.085500
-all phases on all subdomains reached a subsequent error = 1.9997919023890435e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  344 at time t = 0.085750
-all phases on all subdomains reached a subsequent error = 1.9998430219528443e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  345 at time t = 0.086000
-all phases on all subdomains reached a subsequent error = 1.999908690685874e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  346 at time t = 0.086250
-all phases on all subdomains reached a subsequent error = 1.999985977137155e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  347 at time t = 0.086500
-all phases on all subdomains reached a subsequent error = 1.9998031746330366e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  348 at time t = 0.086750
-all phases on all subdomains reached a subsequent error = 1.9997473793962925e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  349 at time t = 0.087000
-all phases on all subdomains reached a subsequent error = 1.9999828134183467e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  350 at time t = 0.087250
-all phases on all subdomains reached a subsequent error = 1.9998564440843108e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  351 at time t = 0.087500
-all phases on all subdomains reached a subsequent error = 1.9998564013723356e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  352 at time t = 0.087750
-all phases on all subdomains reached a subsequent error = 1.9998766388034826e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  353 at time t = 0.088000
-all phases on all subdomains reached a subsequent error = 1.99991241020748e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  354 at time t = 0.088250
-all phases on all subdomains reached a subsequent error = 1.9999647725285033e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  355 at time t = 0.088500
-all phases on all subdomains reached a subsequent error = 1.9997760469386455e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  356 at time t = 0.088750
-all phases on all subdomains reached a subsequent error = 1.9999646793873242e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  357 at time t = 0.089000
-all phases on all subdomains reached a subsequent error = 1.9998147922511992e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  358 at time t = 0.089250
-all phases on all subdomains reached a subsequent error = 1.999785844366449e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  359 at time t = 0.089500
-all phases on all subdomains reached a subsequent error = 1.9997816418015575e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  360 at time t = 0.089750
-all phases on all subdomains reached a subsequent error = 1.9997983028247067e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  361 at time t = 0.090000
-all phases on all subdomains reached a subsequent error = 1.99983586581303e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  362 at time t = 0.090250
-all phases on all subdomains reached a subsequent error = 1.9998910136753246e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  363 at time t = 0.090500
-all phases on all subdomains reached a subsequent error = 1.9999602368144265e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  364 at time t = 0.090750
-all phases on all subdomains reached a subsequent error = 1.9997985340360856e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  365 at time t = 0.091000
-all phases on all subdomains reached a subsequent error = 1.999990778783559e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  366 at time t = 0.091250
-all phases on all subdomains reached a subsequent error = 1.9998616082672365e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  367 at time t = 0.091500
-all phases on all subdomains reached a subsequent error = 1.9998450401859057e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  368 at time t = 0.091750
-all phases on all subdomains reached a subsequent error = 1.999850778606466e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  369 at time t = 0.092000
-all phases on all subdomains reached a subsequent error = 1.999875957983952e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  370 at time t = 0.092250
-all phases on all subdomains reached a subsequent error = 1.9999207975217796e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  371 at time t = 0.092500
-all phases on all subdomains reached a subsequent error = 1.9999821057221016e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  372 at time t = 0.092750
-all phases on all subdomains reached a subsequent error = 1.9998271149441662e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  373 at time t = 0.093000
-all phases on all subdomains reached a subsequent error = 1.9997800763807596e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  374 at time t = 0.093250
-all phases on all subdomains reached a subsequent error = 1.9999839847374263e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  375 at time t = 0.093500
-all phases on all subdomains reached a subsequent error = 1.999888537198252e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  376 at time t = 0.093750
-all phases on all subdomains reached a subsequent error = 1.9998999561465595e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  377 at time t = 0.094000
-all phases on all subdomains reached a subsequent error = 1.9999288681614146e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  378 at time t = 0.094250
-all phases on all subdomains reached a subsequent error = 1.9999728237277873e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  379 at time t = 0.094500
-all phases on all subdomains reached a subsequent error = 1.9998154264321716e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  380 at time t = 0.094750
-all phases on all subdomains reached a subsequent error = 1.999978146917168e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  381 at time t = 0.095000
-all phases on all subdomains reached a subsequent error = 1.9998602901163085e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  382 at time t = 0.095250
-all phases on all subdomains reached a subsequent error = 1.9998450798500277e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  383 at time t = 0.095500
-all phases on all subdomains reached a subsequent error = 1.9998522706277382e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  384 at time t = 0.095750
-all phases on all subdomains reached a subsequent error = 1.999879697585836e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  385 at time t = 0.096000
-all phases on all subdomains reached a subsequent error = 1.999927001831275e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  386 at time t = 0.096250
-all phases on all subdomains reached a subsequent error = 1.999990618908554e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  387 at time t = 0.096500
-all phases on all subdomains reached a subsequent error = 1.999861580281795e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  388 at time t = 0.096750
-all phases on all subdomains reached a subsequent error = 1.9998270138646097e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  389 at time t = 0.097000
-all phases on all subdomains reached a subsequent error = 1.9998145075992755e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  390 at time t = 0.097250
-all phases on all subdomains reached a subsequent error = 1.9998263362928846e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  391 at time t = 0.097500
-all phases on all subdomains reached a subsequent error = 1.9998618901056036e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  392 at time t = 0.097750
-all phases on all subdomains reached a subsequent error = 1.9999166798193944e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  393 at time t = 0.098000
-all phases on all subdomains reached a subsequent error = 1.999986393091332e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  394 at time t = 0.098250
-all phases on all subdomains reached a subsequent error = 1.9998733769090655e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  395 at time t = 0.098500
-all phases on all subdomains reached a subsequent error = 1.9998476123228405e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  396 at time t = 0.098750
-all phases on all subdomains reached a subsequent error = 1.9998424592133683e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  397 at time t = 0.099000
-all phases on all subdomains reached a subsequent error = 1.999860728232292e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  398 at time t = 0.099250
-all phases on all subdomains reached a subsequent error = 1.9999017594619593e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  399 at time t = 0.099500
-all phases on all subdomains reached a subsequent error = 1.9999611193354577e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  400 at time t = 0.099750
-all phases on all subdomains reached a subsequent error = 1.9998510534190394e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  401 at time t = 0.100000
-all phases on all subdomains reached a subsequent error = 1.9998242849015364e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  402 at time t = 0.100250
-all phases on all subdomains reached a subsequent error = 1.999818960250468e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  403 at time t = 0.100500
-all phases on all subdomains reached a subsequent error = 1.9998377676652446e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  404 at time t = 0.100750
-all phases on all subdomains reached a subsequent error = 1.999879556966103e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  405 at time t = 0.101000
-all phases on all subdomains reached a subsequent error = 1.9999396548777105e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  406 at time t = 0.101250
-all phases on all subdomains reached a subsequent error = 1.999840435329867e-06 < tol = 2e-06 after 103 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  407 at time t = 0.101500
-all phases on all subdomains reached a subsequent error = 1.9999910372069267e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  408 at time t = 0.101750
-all phases on all subdomains reached a subsequent error = 1.999923569257241e-06 < tol = 2e-06 after 103 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  409 at time t = 0.102000
-all phases on all subdomains reached a subsequent error = 1.999934120814012e-06 < tol = 2e-06 after 103 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  410 at time t = 0.102250
-all phases on all subdomains reached a subsequent error = 1.999962340057767e-06 < tol = 2e-06 after 103 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  411 at time t = 0.102500
-all phases on all subdomains reached a subsequent error = 1.999844062975549e-06 < tol = 2e-06 after 105 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  412 at time t = 0.102750
-all phases on all subdomains reached a subsequent error = 1.9999677116656813e-06 < tol = 2e-06 after 103 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  413 at time t = 0.103000
-all phases on all subdomains reached a subsequent error = 1.999893040222926e-06 < tol = 2e-06 after 105 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  414 at time t = 0.103250
-all phases on all subdomains reached a subsequent error = 1.9998953876542195e-06 < tol = 2e-06 after 105 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  415 at time t = 0.103500
-all phases on all subdomains reached a subsequent error = 1.9999176883158945e-06 < tol = 2e-06 after 105 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  416 at time t = 0.103750
-all phases on all subdomains reached a subsequent error = 1.9999594748495e-06 < tol = 2e-06 after 105 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  417 at time t = 0.104000
-all phases on all subdomains reached a subsequent error = 1.9998647060560403e-06 < tol = 2e-06 after 107 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  418 at time t = 0.104250
-all phases on all subdomains reached a subsequent error = 1.999994830794208e-06 < tol = 2e-06 after 105 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  419 at time t = 0.104500
-all phases on all subdomains reached a subsequent error = 1.999937525775106e-06 < tol = 2e-06 after 107 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  420 at time t = 0.104750
-all phases on all subdomains reached a subsequent error = 1.999950290930565e-06 < tol = 2e-06 after 107 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  421 at time t = 0.105000
-all phases on all subdomains reached a subsequent error = 1.999980813251091e-06 < tol = 2e-06 after 107 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  422 at time t = 0.105250
-all phases on all subdomains reached a subsequent error = 1.999882956956837e-06 < tol = 2e-06 after 109 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  423 at time t = 0.105500
-all phases on all subdomains reached a subsequent error = 1.999999569194975e-06 < tol = 2e-06 after 107 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  424 at time t = 0.105750
-all phases on all subdomains reached a subsequent error = 1.999943943757847e-06 < tol = 2e-06 after 109 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  425 at time t = 0.106000
-all phases on all subdomains reached a subsequent error = 1.999955579752383e-06 < tol = 2e-06 after 109 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  426 at time t = 0.106250
-all phases on all subdomains reached a subsequent error = 1.999985650379983e-06 < tol = 2e-06 after 109 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  427 at time t = 0.106500
-all phases on all subdomains reached a subsequent error = 1.9998961164726202e-06 < tol = 2e-06 after 111 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  428 at time t = 0.106750
-all phases on all subdomains reached a subsequent error = 1.9998715983169138e-06 < tol = 2e-06 after 111 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  429 at time t = 0.107000
-all phases on all subdomains reached a subsequent error = 1.999871539892865e-06 < tol = 2e-06 after 111 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  430 at time t = 0.107250
-all phases on all subdomains reached a subsequent error = 1.9998977533511366e-06 < tol = 2e-06 after 111 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  431 at time t = 0.107500
-all phases on all subdomains reached a subsequent error = 1.9999469915545982e-06 < tol = 2e-06 after 111 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  432 at time t = 0.107750
-all phases on all subdomains reached a subsequent error = 1.9998837101828165e-06 < tol = 2e-06 after 113 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  433 at time t = 0.108000
-all phases on all subdomains reached a subsequent error = 1.999876903352534e-06 < tol = 2e-06 after 113 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  434 at time t = 0.108250
-all phases on all subdomains reached a subsequent error = 1.9998904827744e-06 < tol = 2e-06 after 113 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  435 at time t = 0.108500
-all phases on all subdomains reached a subsequent error = 1.999927277888849e-06 < tol = 2e-06 after 113 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  436 at time t = 0.108750
-all phases on all subdomains reached a subsequent error = 1.9999846936424323e-06 < tol = 2e-06 after 113 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  437 at time t = 0.109000
-all phases on all subdomains reached a subsequent error = 1.9999351442507176e-06 < tol = 2e-06 after 115 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  438 at time t = 0.109250
-all phases on all subdomains reached a subsequent error = 1.999936810655722e-06 < tol = 2e-06 after 115 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  439 at time t = 0.109500
-all phases on all subdomains reached a subsequent error = 1.9999575286495832e-06 < tol = 2e-06 after 115 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  440 at time t = 0.109750
-all phases on all subdomains reached a subsequent error = 1.9998843089920348e-06 < tol = 2e-06 after 117 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  441 at time t = 0.110000
-all phases on all subdomains reached a subsequent error = 1.999981759451093e-06 < tol = 2e-06 after 115 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  442 at time t = 0.110250
-all phases on all subdomains reached a subsequent error = 1.999953512270096e-06 < tol = 2e-06 after 117 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  443 at time t = 0.110500
-all phases on all subdomains reached a subsequent error = 1.9999771475176597e-06 < tol = 2e-06 after 117 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  444 at time t = 0.110750
-all phases on all subdomains reached a subsequent error = 1.999908069077264e-06 < tol = 2e-06 after 119 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  445 at time t = 0.111000
-all phases on all subdomains reached a subsequent error = 1.999997912410784e-06 < tol = 2e-06 after 117 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  446 at time t = 0.111250
-all phases on all subdomains reached a subsequent error = 1.9999726618894347e-06 < tol = 2e-06 after 119 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  447 at time t = 0.111500
-all phases on all subdomains reached a subsequent error = 1.999996832621749e-06 < tol = 2e-06 after 119 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  448 at time t = 0.111750
-all phases on all subdomains reached a subsequent error = 1.9999348512697148e-06 < tol = 2e-06 after 121 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  449 at time t = 0.112000
-all phases on all subdomains reached a subsequent error = 1.9999202333604004e-06 < tol = 2e-06 after 121 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  450 at time t = 0.112250
-all phases on all subdomains reached a subsequent error = 1.9999310214263565e-06 < tol = 2e-06 after 121 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  451 at time t = 0.112500
-all phases on all subdomains reached a subsequent error = 1.9999681909765497e-06 < tol = 2e-06 after 121 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  452 at time t = 0.112750
-all phases on all subdomains reached a subsequent error = 1.9999300270520026e-06 < tol = 2e-06 after 123 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  453 at time t = 0.113000
-all phases on all subdomains reached a subsequent error = 1.9999340810342608e-06 < tol = 2e-06 after 123 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  454 at time t = 0.113250
-all phases on all subdomains reached a subsequent error = 1.999959091872133e-06 < tol = 2e-06 after 123 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  455 at time t = 0.113500
-all phases on all subdomains reached a subsequent error = 1.999915209628613e-06 < tol = 2e-06 after 125 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  456 at time t = 0.113750
-all phases on all subdomains reached a subsequent error = 1.9999139406710575e-06 < tol = 2e-06 after 125 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  457 at time t = 0.114000
-all phases on all subdomains reached a subsequent error = 1.9999358896228633e-06 < tol = 2e-06 after 125 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  458 at time t = 0.114250
-all phases on all subdomains reached a subsequent error = 1.9999819563977935e-06 < tol = 2e-06 after 125 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  459 at time t = 0.114500
-all phases on all subdomains reached a subsequent error = 1.999961564140234e-06 < tol = 2e-06 after 127 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  460 at time t = 0.114750
-all phases on all subdomains reached a subsequent error = 1.9999762881613478e-06 < tol = 2e-06 after 127 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  461 at time t = 0.115000
-all phases on all subdomains reached a subsequent error = 1.9999286106939567e-06 < tol = 2e-06 after 129 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  462 at time t = 0.115250
-all phases on all subdomains reached a subsequent error = 1.9999208921035833e-06 < tol = 2e-06 after 129 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  463 at time t = 0.115500
-all phases on all subdomains reached a subsequent error = 1.9999396806411574e-06 < tol = 2e-06 after 129 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  464 at time t = 0.115750
-all phases on all subdomains reached a subsequent error = 1.9999845580807553e-06 < tol = 2e-06 after 129 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  465 at time t = 0.116000
-all phases on all subdomains reached a subsequent error = 1.999973496947399e-06 < tol = 2e-06 after 131 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  466 at time t = 0.116250
-all phases on all subdomains reached a subsequent error = 1.9999931609695902e-06 < tol = 2e-06 after 131 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  467 at time t = 0.116500
-all phases on all subdomains reached a subsequent error = 1.999959403378963e-06 < tol = 2e-06 after 133 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  468 at time t = 0.116750
-all phases on all subdomains reached a subsequent error = 1.999960576932561e-06 < tol = 2e-06 after 133 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  469 at time t = 0.117000
-all phases on all subdomains reached a subsequent error = 1.9999870905908227e-06 < tol = 2e-06 after 133 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  470 at time t = 0.117250
-all phases on all subdomains reached a subsequent error = 1.999970279925692e-06 < tol = 2e-06 after 135 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  471 at time t = 0.117500
-all phases on all subdomains reached a subsequent error = 1.9999847629257167e-06 < tol = 2e-06 after 135 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  472 at time t = 0.117750
-all phases on all subdomains reached a subsequent error = 1.999956788098896e-06 < tol = 2e-06 after 137 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  473 at time t = 0.118000
-all phases on all subdomains reached a subsequent error = 1.9999613248094363e-06 < tol = 2e-06 after 137 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  474 at time t = 0.118250
-all phases on all subdomains reached a subsequent error = 1.9999914201284696e-06 < tol = 2e-06 after 137 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  475 at time t = 0.118500
-all phases on all subdomains reached a subsequent error = 1.9999855222859186e-06 < tol = 2e-06 after 139 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  476 at time t = 0.118750
-all phases on all subdomains reached a subsequent error = 1.9999490702449077e-06 < tol = 2e-06 after 141 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  477 at time t = 0.119000
-all phases on all subdomains reached a subsequent error = 1.9999977420608695e-06 < tol = 2e-06 after 139 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  478 at time t = 0.119250
-all phases on all subdomains reached a subsequent error = 1.999957891832747e-06 < tol = 2e-06 after 143 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  479 at time t = 0.119500
-all phases on all subdomains reached a subsequent error = 1.9999507836287367e-06 < tol = 2e-06 after 143 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  480 at time t = 0.119750
-all phases on all subdomains reached a subsequent error = 1.9999687257376737e-06 < tol = 2e-06 after 143 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  481 at time t = 0.120000
-all phases on all subdomains reached a subsequent error = 1.9999640462757792e-06 < tol = 2e-06 after 145 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  482 at time t = 0.120250
-all phases on all subdomains reached a subsequent error = 1.9999853192084456e-06 < tol = 2e-06 after 145 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  483 at time t = 0.120500
-all phases on all subdomains reached a subsequent error = 1.999981642055573e-06 < tol = 2e-06 after 147 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  484 at time t = 0.120750
-all phases on all subdomains reached a subsequent error = 1.999956948743824e-06 < tol = 2e-06 after 149 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  485 at time t = 0.121000
-all phases on all subdomains reached a subsequent error = 1.9999579836496074e-06 < tol = 2e-06 after 149 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  486 at time t = 0.121250
-all phases on all subdomains reached a subsequent error = 1.9999883576112147e-06 < tol = 2e-06 after 149 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  487 at time t = 0.121500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  488 at time t = 0.121750
-all phases on all subdomains reached a subsequent error = 1.9997683682400495e-06 < tol = 2e-06 after 45 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  489 at time t = 0.122000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  490 at time t = 0.122250
-all phases on all subdomains reached a subsequent error = 1.999493346745467e-06 < tol = 2e-06 after 55 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  491 at time t = 0.122500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  492 at time t = 0.122750
-all phases on all subdomains reached a subsequent error = 1.999781275476984e-06 < tol = 2e-06 after 61 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  493 at time t = 0.123000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  494 at time t = 0.123250
-all phases on all subdomains reached a subsequent error = 1.99984524396632e-06 < tol = 2e-06 after 65 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  495 at time t = 0.123500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  496 at time t = 0.123750
-all phases on all subdomains reached a subsequent error = 1.9995440725520006e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  497 at time t = 0.124000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  498 at time t = 0.124250
-all phases on all subdomains reached a subsequent error = 1.999912490217293e-06 < tol = 2e-06 after 69 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  499 at time t = 0.124500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  500 at time t = 0.124750
-all phases on all subdomains reached a subsequent error = 1.999866810941661e-06 < tol = 2e-06 after 71 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  501 at time t = 0.125000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  502 at time t = 0.125250
-all phases on all subdomains reached a subsequent error = 1.9997528397636112e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  503 at time t = 0.125500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  504 at time t = 0.125750
-all phases on all subdomains reached a subsequent error = 1.999979544984062e-06 < tol = 2e-06 after 73 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  505 at time t = 0.126000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  506 at time t = 0.126250
-all phases on all subdomains reached a subsequent error = 1.999869031936433e-06 < tol = 2e-06 after 75 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  507 at time t = 0.126500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  508 at time t = 0.126750
-all phases on all subdomains reached a subsequent error = 1.99971130651244e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  509 at time t = 0.127000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  510 at time t = 0.127250
-all phases on all subdomains reached a subsequent error = 1.9998710244786948e-06 < tol = 2e-06 after 77 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  511 at time t = 0.127500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  512 at time t = 0.127750
-all phases on all subdomains reached a subsequent error = 1.9997518182405254e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  513 at time t = 0.128000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  514 at time t = 0.128250
-all phases on all subdomains reached a subsequent error = 1.9999217330100747e-06 < tol = 2e-06 after 79 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  515 at time t = 0.128500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  516 at time t = 0.128750
-all phases on all subdomains reached a subsequent error = 1.9998231170669717e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  517 at time t = 0.129000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  518 at time t = 0.129250
-all phases on all subdomains reached a subsequent error = 1.9999929479804307e-06 < tol = 2e-06 after 81 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  519 at time t = 0.129500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  520 at time t = 0.129750
-all phases on all subdomains reached a subsequent error = 1.999910803537671e-06 < tol = 2e-06 after 83 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  521 at time t = 0.130000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  522 at time t = 0.130250
-all phases on all subdomains reached a subsequent error = 1.9997915977493266e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  523 at time t = 0.130500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  524 at time t = 0.130750
-all phases on all subdomains reached a subsequent error = 1.999928570780129e-06 < tol = 2e-06 after 85 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  525 at time t = 0.131000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  526 at time t = 0.131250
-all phases on all subdomains reached a subsequent error = 1.9998609186416554e-06 < tol = 2e-06 after 87 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  527 at time t = 0.131500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  528 at time t = 0.131750
-all phases on all subdomains reached a subsequent error = 1.9997645079690154e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  529 at time t = 0.132000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  530 at time t = 0.132250
-all phases on all subdomains reached a subsequent error = 1.99989925113314e-06 < tol = 2e-06 after 89 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  531 at time t = 0.132500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  532 at time t = 0.132750
-all phases on all subdomains reached a subsequent error = 1.9998592401151995e-06 < tol = 2e-06 after 91 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  533 at time t = 0.133000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  534 at time t = 0.133250
-all phases on all subdomains reached a subsequent error = 1.9997920376425966e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  535 at time t = 0.133500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  536 at time t = 0.133750
-all phases on all subdomains reached a subsequent error = 1.9999310087306925e-06 < tol = 2e-06 after 93 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  537 at time t = 0.134000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  538 at time t = 0.134250
-all phases on all subdomains reached a subsequent error = 1.999919416962591e-06 < tol = 2e-06 after 95 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  539 at time t = 0.134500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  540 at time t = 0.134750
-all phases on all subdomains reached a subsequent error = 1.999881006098991e-06 < tol = 2e-06 after 97 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  541 at time t = 0.135000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  542 at time t = 0.135250
-all phases on all subdomains reached a subsequent error = 1.9998317214649547e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  543 at time t = 0.135500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  544 at time t = 0.135750
-all phases on all subdomains reached a subsequent error = 1.9999669998179473e-06 < tol = 2e-06 after 99 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  545 at time t = 0.136000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  546 at time t = 0.136250
-all phases on all subdomains reached a subsequent error = 1.999987175394182e-06 < tol = 2e-06 after 101 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  547 at time t = 0.136500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  548 at time t = 0.136750
-all phases on all subdomains reached a subsequent error = 1.9999835376725228e-06 < tol = 2e-06 after 103 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  549 at time t = 0.137000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  550 at time t = 0.137250
-all phases on all subdomains reached a subsequent error = 1.9999705596576647e-06 < tol = 2e-06 after 105 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  551 at time t = 0.137500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  552 at time t = 0.137750
-all phases on all subdomains reached a subsequent error = 1.999958424766745e-06 < tol = 2e-06 after 107 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  553 at time t = 0.138000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  554 at time t = 0.138250
-all phases on all subdomains reached a subsequent error = 1.999951734262714e-06 < tol = 2e-06 after 109 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  555 at time t = 0.138500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  556 at time t = 0.138750
-all phases on all subdomains reached a subsequent error = 1.999952598827091e-06 < tol = 2e-06 after 111 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  557 at time t = 0.139000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  558 at time t = 0.139250
-all phases on all subdomains reached a subsequent error = 1.9999619909021433e-06 < tol = 2e-06 after 113 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  559 at time t = 0.139500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  560 at time t = 0.139750
-all phases on all subdomains reached a subsequent error = 1.999980317006799e-06 < tol = 2e-06 after 115 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  561 at time t = 0.140000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  562 at time t = 0.140250
-all phases on all subdomains reached a subsequent error = 1.9998926580991274e-06 < tol = 2e-06 after 119 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  563 at time t = 0.140500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  564 at time t = 0.140750
-all phases on all subdomains reached a subsequent error = 1.9999835272659975e-06 < tol = 2e-06 after 119 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  565 at time t = 0.141000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  566 at time t = 0.141250
-all phases on all subdomains reached a subsequent error = 1.9999525818204977e-06 < tol = 2e-06 after 123 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  567 at time t = 0.141500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  568 at time t = 0.141750
-all phases on all subdomains reached a subsequent error = 1.9999688171019093e-06 < tol = 2e-06 after 125 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  569 at time t = 0.142000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  570 at time t = 0.142250
-all phases on all subdomains reached a subsequent error = 1.9999234672016076e-06 < tol = 2e-06 after 129 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  571 at time t = 0.142500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  572 at time t = 0.142750
-all phases on all subdomains reached a subsequent error = 1.99993774843289e-06 < tol = 2e-06 after 131 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  573 at time t = 0.143000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  574 at time t = 0.143250
-all phases on all subdomains reached a subsequent error = 1.9999916462146503e-06 < tol = 2e-06 after 133 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  575 at time t = 0.143500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  576 at time t = 0.143750
-all phases on all subdomains reached a subsequent error = 1.999994842172298e-06 < tol = 2e-06 after 137 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  577 at time t = 0.144000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  578 at time t = 0.144250
-all phases on all subdomains reached a subsequent error = 1.999973281010627e-06 < tol = 2e-06 after 141 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  579 at time t = 0.144500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  580 at time t = 0.144750
-all phases on all subdomains reached a subsequent error = 1.9999492008937564e-06 < tol = 2e-06 after 145 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  581 at time t = 0.145000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  582 at time t = 0.145250
-all phases on all subdomains reached a subsequent error = 1.999990673255318e-06 < tol = 2e-06 after 147 iterations.
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  583 at time t = 0.145500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  584 at time t = 0.145750
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  585 at time t = 0.146000
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  586 at time t = 0.146250
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  587 at time t = 0.146500
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
-entering timestep  588 at time t = 0.146750
-LDD simulation use case: TP-R-layered-soil-with-inner-patch-all-params-one
diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py
index 618d2e1..5b066bf 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one.py
@@ -911,6 +911,7 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
@@ -918,6 +919,7 @@ print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -956,6 +958,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
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 0fb9ba7..ec5f1d5 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
@@ -914,6 +914,7 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
@@ -921,6 +922,7 @@ print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -959,6 +961,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py
index 622ac56..89e3366 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-all-params-one-mesh-study.py
@@ -912,6 +912,7 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
@@ -919,6 +920,7 @@ print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -957,6 +959,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-pure-dd-mesh-study.py b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-pure-dd-mesh-study.py
index fbd4a8a..9b11c1a 100755
--- a/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-pure-dd-mesh-study.py
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/mesh_study/TP-R-layered_soil_with_inner_patch-pure-dd-mesh-study.py
@@ -911,6 +911,7 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
 # read this file and print it to std out. This way the simulation can produce a
 # log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
@@ -918,6 +919,7 @@ print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -956,6 +958,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd-realistic.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd-realistic.py
index 511a99b..9ae8b6b 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd-realistic.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd-realistic.py
@@ -449,6 +449,15 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
+# read this file and print it to std out. This way the simulation can produce a
+# log file with ./TP-R-layered_soil.py | tee simulation.log
+f = open(thisfile, 'r')
+print(f.read())
+f.close()
+
+
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -461,33 +470,35 @@ for starttime in starttimes:
             mesh_study=mesh_study
             )
 
-        simulation.set_parameters(use_case=use_case,
-                                  output_dir=output_string,
-                                  subdomain_def_points=subdomain_def_points,
-                                  isRichards=isRichards,
-                                  interface_def_points=interface_def_points,
-                                  outer_boundary_def_points=outer_boundary_def_points,
-                                  adjacent_subdomains=adjacent_subdomains,
-                                  mesh_resolution=mesh_resolution,
-                                  viscosity=viscosity,
-                                  porosity=porosity,
-                                  L=L,
-                                  lambda_param=lambda_param,
-                                  relative_permeability=relative_permeability,
-                                  saturation=sat_pressure_relationship,
-                                  starttime=starttime,
-                                  number_of_timesteps=number_of_timesteps,
-                                  number_of_timesteps_to_analyse=number_of_timesteps_to_analyse,
-                                  plot_timestep_every=plot_timestep_every,
-                                  timestep_size=timestep_size,
-                                  sources=source_expression,
-                                  initial_conditions=initial_condition,
-                                  dirichletBC_expression_strings=dirichletBC,
-                                  exact_solution=exact_solution,
-                                  densities=densities,
-                                  include_gravity=include_gravity,
-                                  write2file=write_to_file,
-                                  )
+        simulation.set_parameters(
+            use_case=use_case,
+            output_dir=output_string,
+            subdomain_def_points=subdomain_def_points,
+            isRichards=isRichards,
+            interface_def_points=interface_def_points,
+            outer_boundary_def_points=outer_boundary_def_points,
+            adjacent_subdomains=adjacent_subdomains,
+            mesh_resolution=mesh_resolution,
+            viscosity=viscosity,
+            porosity=porosity,
+            L=L,
+            lambda_param=lambda_param,
+            relative_permeability=relative_permeability,
+            saturation=sat_pressure_relationship,
+            starttime=starttime,
+            number_of_timesteps=number_of_timesteps,
+            number_of_timesteps_to_analyse=number_of_timesteps_to_analyse,
+            plot_timestep_every=plot_timestep_every,
+            timestep_size=timestep_size,
+            sources=source_expression,
+            initial_conditions=initial_condition,
+            dirichletBC_expression_strings=dirichletBC,
+            exact_solution=exact_solution,
+            densities=densities,
+            include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
+            write2file=write_to_file,
+            )
 
         simulation.initialise()
         output_dir = simulation.output_dir
@@ -495,26 +506,39 @@ for starttime in starttimes:
         output = simulation.run(analyse_condition=analyse_condition)
         for subdomain_index, subdomain_output in output.items():
             mesh_h = subdomain_output['mesh_size']
-            for phase, different_errornorms in subdomain_output['errornorm'].items():
-                filename = output_dir + "subdomain{}-space-time-errornorm-{}-phase.csv".format(subdomain_index, phase)
-                # for errortype, errornorm in different_errornorms.items():
-
-                    # eocfile = open("eoc_filename", "a")
-                    # eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
-                    # eocfile.close()
-                    # if subdomain.isRichards:mesh_h
+            for phase, error_dict in subdomain_output['errornorm'].items():
+                filename = output_dir \
+                    + "subdomain{}".format(subdomain_index)\
+                    + "-space-time-errornorm-{}-phase.csv".format(phase)
+                # for errortype, errornorm in error_dict.items():
+
+                # eocfile = open("eoc_filename", "a")
+                # eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
+                # eocfile.close()
+                # if subdomain.isRichards:mesh_h
                 data_dict = {
                     'mesh_parameter': mesh_resolution,
                     'mesh_h': mesh_h,
                 }
-                for error_type, errornorms in different_errornorms.items():
+                for norm_type, errornorm in error_dict.items():
                     data_dict.update(
-                        {error_type: errornorms}
+                        {norm_type: errornorm}
                     )
                 errors = pd.DataFrame(data_dict, index=[mesh_resolution])
                 # check if file exists
-                if os.path.isfile(filename) == True:
+                if os.path.isfile(filename) is True:
                     with open(filename, 'a') as f:
-                        errors.to_csv(f, header=False, sep='\t', encoding='utf-8', index=False)
+                        errors.to_csv(
+                            f,
+                            header=False,
+                            sep='\t',
+                            encoding='utf-8',
+                            index=False
+                            )
                 else:
-                    errors.to_csv(filename, sep='\t', encoding='utf-8', index=False)
+                    errors.to_csv(
+                        filename,
+                        sep='\t',
+                        encoding='utf-8',
+                        index=False
+                        )
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd.py
index 481c7f3..2a42078 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-pure-dd.py
@@ -442,6 +442,15 @@ for subdomain in isRichards.keys():
                 )
 
 
+# LOG FILE OUTPUT #############################################################
+# read this file and print it to std out. This way the simulation can produce a
+# log file with ./TP-R-layered_soil.py | tee simulation.log
+f = open(thisfile, 'r')
+print(f.read())
+f.close()
+
+
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -454,33 +463,35 @@ for starttime in starttimes:
             mesh_study=mesh_study
             )
 
-        simulation.set_parameters(use_case=use_case,
-                                  output_dir=output_string,
-                                  subdomain_def_points=subdomain_def_points,
-                                  isRichards=isRichards,
-                                  interface_def_points=interface_def_points,
-                                  outer_boundary_def_points=outer_boundary_def_points,
-                                  adjacent_subdomains=adjacent_subdomains,
-                                  mesh_resolution=mesh_resolution,
-                                  viscosity=viscosity,
-                                  porosity=porosity,
-                                  L=L,
-                                  lambda_param=lambda_param,
-                                  relative_permeability=relative_permeability,
-                                  saturation=sat_pressure_relationship,
-                                  starttime=starttime,
-                                  number_of_timesteps=number_of_timesteps,
-                                  number_of_timesteps_to_analyse=number_of_timesteps_to_analyse,
-                                  plot_timestep_every=plot_timestep_every,
-                                  timestep_size=timestep_size,
-                                  sources=source_expression,
-                                  initial_conditions=initial_condition,
-                                  dirichletBC_expression_strings=dirichletBC,
-                                  exact_solution=exact_solution,
-                                  densities=densities,
-                                  include_gravity=include_gravity,
-                                  write2file=write_to_file,
-                                  )
+        simulation.set_parameters(
+            use_case=use_case,
+            output_dir=output_string,
+            subdomain_def_points=subdomain_def_points,
+            isRichards=isRichards,
+            interface_def_points=interface_def_points,
+            outer_boundary_def_points=outer_boundary_def_points,
+            adjacent_subdomains=adjacent_subdomains,
+            mesh_resolution=mesh_resolution,
+            viscosity=viscosity,
+            porosity=porosity,
+            L=L,
+            lambda_param=lambda_param,
+            relative_permeability=relative_permeability,
+            saturation=sat_pressure_relationship,
+            starttime=starttime,
+            number_of_timesteps=number_of_timesteps,
+            number_of_timesteps_to_analyse=number_of_timesteps_to_analyse,
+            plot_timestep_every=plot_timestep_every,
+            timestep_size=timestep_size,
+            sources=source_expression,
+            initial_conditions=initial_condition,
+            dirichletBC_expression_strings=dirichletBC,
+            exact_solution=exact_solution,
+            densities=densities,
+            include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
+            write2file=write_to_file,
+            )
 
         simulation.initialise()
         output_dir = simulation.output_dir
@@ -488,26 +499,39 @@ for starttime in starttimes:
         output = simulation.run(analyse_condition=analyse_condition)
         for subdomain_index, subdomain_output in output.items():
             mesh_h = subdomain_output['mesh_size']
-            for phase, different_errornorms in subdomain_output['errornorm'].items():
-                filename = output_dir + "subdomain{}-space-time-errornorm-{}-phase.csv".format(subdomain_index, phase)
-                # for errortype, errornorm in different_errornorms.items():
-
-                    # eocfile = open("eoc_filename", "a")
-                    # eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
-                    # eocfile.close()
-                    # if subdomain.isRichards:mesh_h
+            for phase, error_dict in subdomain_output['errornorm'].items():
+                filename = output_dir \
+                    + "subdomain{}".format(subdomain_index)\
+                    + "-space-time-errornorm-{}-phase.csv".format(phase)
+                # for errortype, errornorm in error_dict.items():
+
+                # eocfile = open("eoc_filename", "a")
+                # eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
+                # eocfile.close()
+                # if subdomain.isRichards:mesh_h
                 data_dict = {
                     'mesh_parameter': mesh_resolution,
                     'mesh_h': mesh_h,
                 }
-                for error_type, errornorms in different_errornorms.items():
+                for norm_type, errornorm in error_dict.items():
                     data_dict.update(
-                        {error_type: errornorms}
+                        {norm_type: errornorm}
                     )
                 errors = pd.DataFrame(data_dict, index=[mesh_resolution])
                 # check if file exists
-                if os.path.isfile(filename) == True:
+                if os.path.isfile(filename) is True:
                     with open(filename, 'a') as f:
-                        errors.to_csv(f, header=False, sep='\t', encoding='utf-8', index=False)
+                        errors.to_csv(
+                            f,
+                            header=False,
+                            sep='\t',
+                            encoding='utf-8',
+                            index=False
+                            )
                 else:
-                    errors.to_csv(filename, sep='\t', encoding='utf-8', index=False)
+                    errors.to_csv(
+                        filename,
+                        sep='\t',
+                        encoding='utf-8',
+                        index=False
+                        )
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic.py
index 3853c34..bfc3030 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-realistic.py
@@ -435,7 +435,15 @@ for subdomain in isRichards.keys():
                 {outer_boundary_ind: exact_solution[subdomain]}
                 )
 
+# LOG FILE OUTPUT #############################################################
+# read this file and print it to std out. This way the simulation can produce a
+# log file with ./TP-R-layered_soil.py | tee simulation.log
+f = open(thisfile, 'r')
+print(f.read())
+f.close()
 
+
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -448,33 +456,35 @@ for starttime in starttimes:
             mesh_study=mesh_study
             )
 
-        simulation.set_parameters(use_case=use_case,
-                                  output_dir=output_string,
-                                  subdomain_def_points=subdomain_def_points,
-                                  isRichards=isRichards,
-                                  interface_def_points=interface_def_points,
-                                  outer_boundary_def_points=outer_boundary_def_points,
-                                  adjacent_subdomains=adjacent_subdomains,
-                                  mesh_resolution=mesh_resolution,
-                                  viscosity=viscosity,
-                                  porosity=porosity,
-                                  L=L,
-                                  lambda_param=lambda_param,
-                                  relative_permeability=relative_permeability,
-                                  saturation=sat_pressure_relationship,
-                                  starttime=starttime,
-                                  number_of_timesteps=number_of_timesteps,
-                                  number_of_timesteps_to_analyse=number_of_timesteps_to_analyse,
-                                  plot_timestep_every=plot_timestep_every,
-                                  timestep_size=timestep_size,
-                                  sources=source_expression,
-                                  initial_conditions=initial_condition,
-                                  dirichletBC_expression_strings=dirichletBC,
-                                  exact_solution=exact_solution,
-                                  densities=densities,
-                                  include_gravity=include_gravity,
-                                  write2file=write_to_file,
-                                  )
+        simulation.set_parameters(
+            use_case=use_case,
+            output_dir=output_string,
+            subdomain_def_points=subdomain_def_points,
+            isRichards=isRichards,
+            interface_def_points=interface_def_points,
+            outer_boundary_def_points=outer_boundary_def_points,
+            adjacent_subdomains=adjacent_subdomains,
+            mesh_resolution=mesh_resolution,
+            viscosity=viscosity,
+            porosity=porosity,
+            L=L,
+            lambda_param=lambda_param,
+            relative_permeability=relative_permeability,
+            saturation=sat_pressure_relationship,
+            starttime=starttime,
+            number_of_timesteps=number_of_timesteps,
+            number_of_timesteps_to_analyse=number_of_timesteps_to_analyse,
+            plot_timestep_every=plot_timestep_every,
+            timestep_size=timestep_size,
+            sources=source_expression,
+            initial_conditions=initial_condition,
+            dirichletBC_expression_strings=dirichletBC,
+            exact_solution=exact_solution,
+            densities=densities,
+            include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
+            write2file=write_to_file,
+            )
 
         simulation.initialise()
         output_dir = simulation.output_dir
@@ -482,26 +492,39 @@ for starttime in starttimes:
         output = simulation.run(analyse_condition=analyse_condition)
         for subdomain_index, subdomain_output in output.items():
             mesh_h = subdomain_output['mesh_size']
-            for phase, different_errornorms in subdomain_output['errornorm'].items():
-                filename = output_dir + "subdomain{}-space-time-errornorm-{}-phase.csv".format(subdomain_index, phase)
-                # for errortype, errornorm in different_errornorms.items():
-
-                    # eocfile = open("eoc_filename", "a")
-                    # eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
-                    # eocfile.close()
-                    # if subdomain.isRichards:mesh_h
+            for phase, error_dict in subdomain_output['errornorm'].items():
+                filename = output_dir \
+                    + "subdomain{}".format(subdomain_index)\
+                    + "-space-time-errornorm-{}-phase.csv".format(phase)
+                # for errortype, errornorm in error_dict.items():
+
+                # eocfile = open("eoc_filename", "a")
+                # eocfile.write( str(mesh_h) + " " + str(errornorm) + "\n" )
+                # eocfile.close()
+                # if subdomain.isRichards:mesh_h
                 data_dict = {
                     'mesh_parameter': mesh_resolution,
                     'mesh_h': mesh_h,
                 }
-                for error_type, errornorms in different_errornorms.items():
+                for norm_type, errornorm in error_dict.items():
                     data_dict.update(
-                        {error_type: errornorms}
+                        {norm_type: errornorm}
                     )
                 errors = pd.DataFrame(data_dict, index=[mesh_resolution])
                 # check if file exists
-                if os.path.isfile(filename) == True:
+                if os.path.isfile(filename) is True:
                     with open(filename, 'a') as f:
-                        errors.to_csv(f, header=False, sep='\t', encoding='utf-8', index=False)
+                        errors.to_csv(
+                            f,
+                            header=False,
+                            sep='\t',
+                            encoding='utf-8',
+                            index=False
+                            )
                 else:
-                    errors.to_csv(filename, sep='\t', encoding='utf-8', index=False)
+                    errors.to_csv(
+                        filename,
+                        sep='\t',
+                        encoding='utf-8',
+                        index=False
+                        )
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-test.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-test.py
index b85e062..116b1e1 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-test.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-patch-test.py
@@ -444,11 +444,16 @@ for subdomain in isRichards.keys():
                 {outer_boundary_ind: exact_solution[subdomain]}
                 )
 
+
+# LOG FILE OUTPUT #############################################################
+# read this file and print it to std out. This way the simulation can produce a
+# log file with ./TP-R-layered_soil.py | tee simulation.log
 f = open(thisfile, 'r')
 print(f.read())
 f.close()
 
 
+# RUN #########################################################################
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -487,6 +492,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-realistic-parameters-densities-scaled.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-realistic-parameters-densities-scaled.py
index 8bef020..df79a3d 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-realistic-parameters-densities-scaled.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/TP-R-2-realistic-parameters-densities-scaled.py
@@ -500,6 +500,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one-but-g.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one-but-g.py
index 20e49f3..8a6c2eb 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one-but-g.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one-but-g.py
@@ -542,6 +542,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one.py
index 64512a6..a6139d8 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study-all-params-one.py
@@ -48,8 +48,8 @@ resolutions = {
                 8: 5e-6,
                 16: 5e-6,
                 32: 5e-6,
-                64: 5e-6,
-                128: 5e-6,
+                64: 2e-6,
+                128: 2e-6,
                 # 256: 1e-6,
                 }
 
@@ -58,7 +58,7 @@ resolutions = {
 #  for each element t_0 in starttimes.
 starttimes = [0.0]
 timestep_size = 0.001
-number_of_timesteps = 800
+number_of_timesteps = 20
 
 # LDD scheme parameters  ######################################################
 Lw1 = 0.25 #/timestep_size
@@ -504,6 +504,7 @@ f.close()
 
 
 # RUN #########################################################################
+
 for starttime in starttimes:
     for mesh_resolution, solver_tol in resolutions.items():
         # initialise LDD simulation class
@@ -542,6 +543,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
diff --git a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study.py b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study.py
index fd1b6c0..07eb650 100755
--- a/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study.py
+++ b/Two-phase-Richards/two-patch/TP-R-two-patch-test-case/mesh_studies/TP-R-2-patch-mesh-study.py
@@ -542,6 +542,7 @@ for starttime in starttimes:
             exact_solution=exact_solution,
             densities=densities,
             include_gravity=include_gravity,
+            gravity_acceleration=gravity_acceleration,
             write2file=write_to_file,
             )
 
-- 
GitLab