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
new file mode 100644
index 0000000000000000000000000000000000000000..c4f9167a586c5adbfdf377dd87b40b0831657247
--- /dev/null
+++ b/Two-phase-Richards/multi-patch/layered_soil_with_inner_patch/TP-R-layered_soil_with_inner_patch-all-params-one-finer-dt.py
@@ -0,0 +1,3025 @@
+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'>
+ type of dS is Piecewise
+ type of dtpc is <class 'sympy.core.add.Add'>
+ type of dS is Piecewise
+ type of dtpc is <class 'sympy.core.add.Add'>
+ type of dS is Piecewise
+ type of dtpc is <class 'sympy.core.add.Add'>
+ type of dS is Piecewise
+ type of dtpc is <class 'sympy.core.add.Add'>
+ type of dS is Piecewise
+#!/usr/bin/python3
+"""Layered soil simulation with inner patch.
+
+This program sets up an LDD simulation
+"""
+
+import dolfin as df
+import sympy as sym
+import functools as ft
+import LDDsimulation as ldd
+import helpers as hlp
+import datetime
+import os
+import pandas as pd
+
+# check if output directory exists
+if not os.path.exists('./output'):
+ os.mkdir('./output')
+ print("Directory ", './output', " created ")
+else:
+ print("Directory ", './output', " already exists. Will use as output \
+ directory")
+
+
+date = datetime.datetime.now()
+datestr = date.strftime("%Y-%m-%d")
+
+# init sympy session
+sym.init_printing()
+# solver_tol = 6E-7
+use_case = "TP-R-layered-soil-with-inner-patch-all-params-one"
+# name of this very file. Needed for log output.
+thisfile = "TP-R-layered_soil_with_inner_patch-all-params-one-finer-dt.py"
+
+max_iter_num = 150
+FEM_Lagrange_degree = 1
+mesh_study = False
+resolutions = {
+ # 1: 2e-6, # h=2
+ # 2: 2e-6, # h=1.1180
+ # 4: 2e-6, # h=0.5590
+ # 8: 2e-6, # h=0.2814
+ # 16: 7e-6, # h=0.1412
+ # 32: 5e-6,
+ 64: 2e-6,
+ # 128: 2e-6
+ }
+
+# GRID #######################
+# mesh_resolution = 20
+timestep_size = 0.00025
+number_of_timesteps = 4000
+plot_timestep_every = 4
+# decide how many timesteps you want analysed. Analysed means, that we write
+# out subsequent errors of the L-iteration within the timestep.
+number_of_timesteps_to_analyse = 5
+starttimes = [0.0]
+
+Lw1 = 0.25 # /timestep_size
+Lnw1 = Lw1
+
+Lw2 = 0.25 # /timestep_size
+Lnw2 = Lw2
+
+Lw3 = 0.25 # /timestep_size
+Lnw3 = Lw3
+
+Lw4 = 0.25 # /timestep_size
+Lnw4 = Lw4
+
+Lw5 = 0.25 # /timestep_size
+Lnw5 = Lw5
+
+Lw6 = 0.25 # /timestep_size
+Lnw6 = Lw6
+
+lambda12_w = 4
+lambda12_nw = 4
+
+lambda23_w = 4
+lambda23_nw = 4
+
+lambda24_w = 4
+lambda24_nw= 4
+
+lambda25_w= 4
+lambda25_nw= 4
+
+lambda34_w = 4
+lambda34_nw = 4
+
+lambda36_w = 4
+lambda36_nw = 4
+
+lambda45_w = 4
+lambda45_nw = 4
+
+lambda46_w = 4
+lambda46_nw = 4
+
+lambda56_w = 4
+lambda56_nw = 4
+
+include_gravity = False
+debugflag = False
+analyse_condition = False
+
+# if mesh_study:
+output_string = "./output/{}-{}_timesteps{}_P{}".format(
+ datestr, use_case, number_of_timesteps, FEM_Lagrange_degree
+ )
+# else:
+# # for tol in resolutions.values():
+# # solver_tol = tol
+# output_string = "./output/{}-{}_timesteps{}_P{}_solver_tol{}".format(
+# datestr, use_case, number_of_timesteps, FEM_Lagrange_degree, solver_tol
+# )
+
+
+# toggle what should be written to files
+if mesh_study:
+ write_to_file = {
+ 'space_errornorms': True,
+ 'meshes_and_markers': True,
+ 'L_iterations_per_timestep': True,
+ 'solutions': True,
+ 'absolute_differences': True,
+ 'condition_numbers': analyse_condition,
+ 'subsequent_errors': True
+ }
+else:
+ write_to_file = {
+ 'space_errornorms': True,
+ 'meshes_and_markers': True,
+ 'L_iterations_per_timestep': True,
+ 'solutions': True,
+ 'absolute_differences': True,
+ 'condition_numbers': analyse_condition,
+ 'subsequent_errors': True
+ }
+
+
+# global domain
+subdomain0_vertices = [df.Point(-1.0,-1.0), #
+ df.Point(1.0,-1.0),#
+ df.Point(1.0,1.0),#
+ df.Point(-1.0,1.0)]
+
+interface12_vertices = [df.Point(-1.0, 0.8),
+ df.Point(0.3, 0.8),
+ df.Point(0.5, 0.9),
+ df.Point(0.8, 0.7),
+ df.Point(1.0, 0.65)]
+
+
+ # interface23
+interface23_vertices = [df.Point(-1.0, 0.0),
+ df.Point(-0.35, 0.0),
+ # df.Point(6.5, 4.5),
+ df.Point(0.0, 0.0)]
+
+interface24_vertices = [interface23_vertices[2],
+ df.Point(0.6, 0.0),
+ ]
+
+interface25_vertices = [interface24_vertices[1],
+ df.Point(1.0, 0.0)
+ ]
+
+
+interface32_vertices = [interface23_vertices[2],
+ interface23_vertices[1],
+ interface23_vertices[0]]
+
+
+interface36_vertices = [df.Point(-1.0, -0.6),
+ df.Point(-0.6, -0.45)]
+
+
+interface46_vertices = [interface36_vertices[1],
+ df.Point(0.3, -0.25)]
+
+interface56_vertices = [interface46_vertices[1],
+ df.Point(0.65, -0.6),
+ df.Point(1.0, -0.7)]
+
+
+
+
+interface34_vertices = [interface36_vertices[1],
+ interface23_vertices[2]]
+# interface36
+
+interface45_vertices = [interface56_vertices[0],
+ df.Point(0.7, -0.2),
+ interface25_vertices[0]
+ ]
+
+
+# interface_vertices introduces a global numbering of interfaces.
+interface_def_points = [interface12_vertices,
+ interface23_vertices,
+ interface24_vertices,
+ interface25_vertices,
+ interface34_vertices,
+ interface36_vertices,
+ interface45_vertices,
+ interface46_vertices,
+ interface56_vertices,
+ ]
+
+adjacent_subdomains = [[1,2],
+ [2,3],
+ [2,4],
+ [2,5],
+ [3,4],
+ [3,6],
+ [4,5],
+ [4,6],
+ [5,6]
+ ]
+
+# subdomain1.
+subdomain1_vertices = [interface12_vertices[0],
+ interface12_vertices[1],
+ interface12_vertices[2],
+ interface12_vertices[3],
+ interface12_vertices[4], # southern boundary, 12 interface
+ subdomain0_vertices[2], # eastern boundary, outer boundary
+ subdomain0_vertices[3]] # northern boundary, outer on_boundary
+
+# vertex coordinates of the outer boundaries. If it can not be specified as a
+# polygon, use an entry per boundary polygon.
+# This information is used for defining
+# the Dirichlet boundary conditions. If a domain is completely internal, the
+# dictionary entry should be 0: None
+subdomain1_outer_boundary_verts = {
+ 0: [subdomain1_vertices[4],
+ subdomain1_vertices[5], # eastern boundary, outer boundary
+ subdomain1_vertices[6],
+ subdomain1_vertices[0]]
+}
+
+#subdomain1
+subdomain2_vertices = [interface23_vertices[0],
+ interface23_vertices[1],
+ interface23_vertices[2],
+ interface24_vertices[1],
+ interface25_vertices[1], # southern boundary, 23 interface
+ subdomain1_vertices[4], # eastern boundary, outer boundary
+ subdomain1_vertices[3],
+ subdomain1_vertices[2],
+ subdomain1_vertices[1],
+ subdomain1_vertices[0] ] # northern boundary, 12 interface
+
+subdomain2_outer_boundary_verts = {
+ 0: [subdomain2_vertices[9],
+ subdomain2_vertices[0]],
+ 1: [subdomain2_vertices[4],
+ subdomain2_vertices[5]]
+}
+
+
+subdomain3_vertices = [interface36_vertices[0],
+ interface36_vertices[1],
+ # interface34_vertices[0],
+ interface34_vertices[1],
+ # interface32_vertices[0],
+ interface32_vertices[1],
+ interface32_vertices[2]
+ ]
+
+subdomain3_outer_boundary_verts = {
+ 0: [subdomain3_vertices[4],
+ subdomain3_vertices[0]]
+}
+
+
+# subdomain3
+subdomain4_vertices = [interface46_vertices[0],
+ interface46_vertices[1],
+ interface45_vertices[1],
+ interface24_vertices[1],
+ interface24_vertices[0],
+ interface34_vertices[1]
+ ]
+
+subdomain4_outer_boundary_verts = None
+
+subdomain5_vertices = [interface56_vertices[0],
+ interface56_vertices[1],
+ interface56_vertices[2],
+ interface25_vertices[1],
+ interface25_vertices[0],
+ interface45_vertices[1],
+ interface45_vertices[0]
+]
+
+subdomain5_outer_boundary_verts = {
+ 0: [subdomain5_vertices[2],
+ subdomain5_vertices[3]]
+}
+
+
+
+subdomain6_vertices = [subdomain0_vertices[0],
+ subdomain0_vertices[1], # southern boundary, outer boundary
+ interface56_vertices[2],
+ interface56_vertices[1],
+ interface56_vertices[0],
+ interface36_vertices[1],
+ interface36_vertices[0]
+ ]
+
+subdomain6_outer_boundary_verts = {
+ 0: [subdomain6_vertices[6],
+ subdomain6_vertices[0],
+ subdomain6_vertices[1],
+ subdomain6_vertices[2]]
+}
+
+
+subdomain_def_points = [subdomain0_vertices,#
+ subdomain1_vertices,#
+ subdomain2_vertices,#
+ subdomain3_vertices,#
+ subdomain4_vertices,
+ subdomain5_vertices,
+ subdomain6_vertices
+ ]
+
+
+# if a subdomain has no outer boundary write None instead, i.e.
+# i: None
+# if i is the index of the inner subdomain.
+outer_boundary_def_points = {
+ # subdomain number
+ 1: subdomain1_outer_boundary_verts,
+ 2: subdomain2_outer_boundary_verts,
+ 3: subdomain3_outer_boundary_verts,
+ 4: subdomain4_outer_boundary_verts,
+ 5: subdomain5_outer_boundary_verts,
+ 6: subdomain6_outer_boundary_verts
+}
+
+
+isRichards = {
+ 1: True,
+ 2: True,
+ 3: False,
+ 4: False,
+ 5: False,
+ 6: False
+ }
+
+# isRichards = {
+# 1: True,
+# 2: True,
+# 3: True,
+# 4: True,
+# 5: True,
+# 6: True
+# }
+
+# Dict of the form: { subdom_num : viscosity }
+viscosity = {
+ 1: {'wetting' :1,
+ 'nonwetting': 1},
+ 2: {'wetting' :1,
+ 'nonwetting': 1},
+ 3: {'wetting' :1,
+ 'nonwetting': 1},
+ 4: {'wetting' :1,
+ 'nonwetting': 1},
+ 5: {'wetting' :1,
+ 'nonwetting': 1},
+ 6: {'wetting' :1,
+ 'nonwetting': 1},
+}
+
+# Dict of the form: { subdom_num : density }
+densities = {
+ 1: {'wetting': 1, #997
+ 'nonwetting': 1}, #1}, #1.225},
+ 2: {'wetting': 1, #997
+ 'nonwetting': 1}, #1.225},
+ 3: {'wetting': 1, #997
+ 'nonwetting': 1}, #1.225},
+ 4: {'wetting': 1, #997
+ 'nonwetting': 1}, #1.225}
+ 5: {'wetting': 1, #997
+ 'nonwetting': 1}, #1.225},
+ 6: {'wetting': 1, #997
+ 'nonwetting': 1} #1.225}
+}
+
+gravity_acceleration = 9.81
+# porosities taken from
+# https://www.geotechdata.info/parameter/soil-porosity.html
+# Dict of the form: { subdom_num : porosity }
+porosity = {
+ 1: 1, #0.2, # Clayey gravels, clayey sandy gravels
+ 2: 1, #0.22, # Silty gravels, silty sandy gravels
+ 3: 1, #0.37, # Clayey sands
+ 4: 1, #0.2 # Silty or sandy clay
+ 5: 1, #
+ 6: 1, #
+}
+
+# subdom_num : subdomain L for L-scheme
+L = {
+ 1: {'wetting' :Lw1,
+ 'nonwetting': Lnw1},
+ 2: {'wetting' :Lw2,
+ 'nonwetting': Lnw2},
+ 3: {'wetting' :Lw3,
+ 'nonwetting': Lnw3},
+ 4: {'wetting' :Lw4,
+ 'nonwetting': Lnw4},
+ 5: {'wetting' :Lw5,
+ 'nonwetting': Lnw5},
+ 6: {'wetting' :Lw6,
+ 'nonwetting': Lnw6}
+}
+
+
+# interface_num : lambda parameter for the L-scheme on that interface.
+# Note that interfaces are numbered starting from 0, because
+# adjacent_subdomains is a list and not a dict. Historic fuckup, I know
+# We have defined above as interfaces
+# # interface_vertices introduces a global numbering of interfaces.
+# interface_def_points = [interface12_vertices,
+# interface23_vertices,
+# interface24_vertices,
+# interface25_vertices,
+# interface34_vertices,
+# interface36_vertices,
+# interface45_vertices,
+# interface46_vertices,
+# interface56_vertices,
+# ]
+lambda_param = {
+ 0: {'wetting': lambda12_w,
+ 'nonwetting': lambda12_nw},#
+ 1: {'wetting': lambda23_w,
+ 'nonwetting': lambda23_nw},#
+ 2: {'wetting': lambda24_w,
+ 'nonwetting': lambda24_nw},#
+ 3: {'wetting': lambda25_w,
+ 'nonwetting': lambda25_nw},#
+ 4: {'wetting': lambda34_w,
+ 'nonwetting': lambda34_nw},#
+ 5: {'wetting': lambda36_w,
+ 'nonwetting': lambda36_nw},#
+ 6: {'wetting': lambda45_w,
+ 'nonwetting': lambda45_nw},#
+ 7: {'wetting': lambda46_w,
+ 'nonwetting': lambda46_nw},#
+ 8: {'wetting': lambda56_w,
+ 'nonwetting': lambda56_nw},#
+}
+
+
+# after Lewis, see pdf file
+intrinsic_permeability = {
+ 1: 1, # sand
+ 2: 1, # sand, there is a range
+ 3: 1, #10e-2, # clay has a range
+ 4: 1, #10e-3
+ 5: 1, #10e-2, # clay has a range
+ 6: 1, #10e-3
+}
+
+
+# relative permeabilty functions on subdomain 1
+def rel_perm1w(s):
+ # relative permeabilty wetting on subdomain1
+ return intrinsic_permeability[1]*s**2
+
+
+def rel_perm1nw(s):
+ # relative permeabilty nonwetting on subdomain1
+ return intrinsic_permeability[1]*(1-s)**2
+
+
+# relative permeabilty functions on subdomain 2
+def rel_perm2w(s):
+ # relative permeabilty wetting on subdomain2
+ return intrinsic_permeability[2]*s**2
+
+
+def rel_perm2nw(s):
+ # relative permeabilty nonwetting on subdomain2
+ return intrinsic_permeability[2]*(1-s)**2
+
+
+# relative permeabilty functions on subdomain 3
+def rel_perm3w(s):
+ # relative permeabilty wetting on subdomain3
+ return intrinsic_permeability[3]*s**3
+
+
+def rel_perm3nw(s):
+ # relative permeabilty nonwetting on subdomain3
+ return intrinsic_permeability[3]*(1-s)**3
+
+
+# relative permeabilty functions on subdomain 4
+def rel_perm4w(s):
+ # relative permeabilty wetting on subdomain4
+ return intrinsic_permeability[4]*s**3
+
+
+def rel_perm4nw(s):
+ # relative permeabilty nonwetting on subdomain4
+ return intrinsic_permeability[4]*(1-s)**3
+
+
+# relative permeabilty functions on subdomain 5
+def rel_perm5w(s):
+ # relative permeabilty wetting on subdomain5
+ return intrinsic_permeability[5]*s**3
+
+
+def rel_perm5nw(s):
+ # relative permeabilty nonwetting on subdomain5
+ return intrinsic_permeability[5]*(1-s)**3
+
+
+# relative permeabilty functions on subdomain 6
+def rel_perm6w(s):
+ # relative permeabilty wetting on subdomain6
+ return intrinsic_permeability[6]*s**3
+
+
+def rel_perm6nw(s):
+ # relative permeabilty nonwetting on subdomain6
+ return intrinsic_permeability[6]*(1-s)**3
+
+
+_rel_perm1w = ft.partial(rel_perm1w)
+_rel_perm1nw = ft.partial(rel_perm1nw)
+
+_rel_perm2w = ft.partial(rel_perm2w)
+_rel_perm2nw = ft.partial(rel_perm2nw)
+
+_rel_perm3w = ft.partial(rel_perm3w)
+_rel_perm3nw = ft.partial(rel_perm3nw)
+
+_rel_perm4w = ft.partial(rel_perm4w)
+_rel_perm4nw = ft.partial(rel_perm4nw)
+
+_rel_perm5w = ft.partial(rel_perm5w)
+_rel_perm5nw = ft.partial(rel_perm5nw)
+
+_rel_perm6w = ft.partial(rel_perm6w)
+_rel_perm6nw = ft.partial(rel_perm6nw)
+
+subdomain1_rel_perm = {
+ 'wetting': _rel_perm1w,
+ 'nonwetting': _rel_perm1nw
+}
+
+subdomain2_rel_perm = {
+ 'wetting': _rel_perm2w,
+ 'nonwetting': _rel_perm2nw
+}
+
+subdomain3_rel_perm = {
+ 'wetting': _rel_perm3w,
+ 'nonwetting': _rel_perm3nw
+}
+
+subdomain4_rel_perm = {
+ 'wetting': _rel_perm4w,
+ 'nonwetting': _rel_perm4nw
+}
+
+subdomain5_rel_perm = {
+ 'wetting': _rel_perm5w,
+ 'nonwetting': _rel_perm5nw
+}
+
+subdomain6_rel_perm = {
+ 'wetting': _rel_perm6w,
+ 'nonwetting': _rel_perm6nw
+}
+
+# dictionary of relative permeabilties on all domains.
+relative_permeability = {
+ 1: subdomain1_rel_perm,
+ 2: subdomain2_rel_perm,
+ 3: subdomain3_rel_perm,
+ 4: subdomain4_rel_perm,
+ 5: subdomain5_rel_perm,
+ 6: subdomain6_rel_perm
+}
+
+
+# definition of the derivatives of the relative permeabilities
+# relative permeabilty functions on subdomain 1
+def rel_perm1w_prime(s):
+ # relative permeabilty on subdomain1
+ return intrinsic_permeability[1]*2*s
+
+
+def rel_perm1nw_prime(s):
+ # relative permeabilty on subdomain1
+ return -1*intrinsic_permeability[1]*2*(1-s)
+
+
+def rel_perm2w_prime(s):
+ # relative permeabilty on subdomain2
+ return intrinsic_permeability[2]*2*s
+
+
+def rel_perm2nw_prime(s):
+ # relative permeabilty on subdomain2
+ return -1*intrinsic_permeability[2]*2*(1-s)
+
+
+# definition of the derivatives of the relative permeabilities
+# relative permeabilty functions on subdomain 3
+def rel_perm3w_prime(s):
+ # relative permeabilty on subdomain3
+ return intrinsic_permeability[3]*3*s**2
+
+
+def rel_perm3nw_prime(s):
+ # relative permeabilty on subdomain3
+ return -1*intrinsic_permeability[3]*3*(1-s)**2
+
+
+# definition of the derivatives of the relative permeabilities
+# relative permeabilty functions on subdomain 4
+def rel_perm4w_prime(s):
+ # relative permeabilty on subdomain4
+ return intrinsic_permeability[4]*3*s**2
+
+
+def rel_perm4nw_prime(s):
+ # relative permeabilty on subdomain4
+ return -1*intrinsic_permeability[4]*3*(1-s)**2
+
+
+# definition of the derivatives of the relative permeabilities
+# relative permeabilty functions on subdomain 5
+def rel_perm5w_prime(s):
+ # relative permeabilty on subdomain5
+ return intrinsic_permeability[5]*3*s**2
+
+
+def rel_perm5nw_prime(s):
+ # relative permeabilty on subdomain5
+ return -1*intrinsic_permeability[5]*3*(1-s)**2
+
+
+# definition of the derivatives of the relative permeabilities
+# relative permeabilty functions on subdomain 6
+def rel_perm6w_prime(s):
+ # relative permeabilty on subdomain6
+ return intrinsic_permeability[6]*3*s**2
+
+
+def rel_perm6nw_prime(s):
+ # relative permeabilty on subdomain6
+ return -1*intrinsic_permeability[6]*3*(1-s)**2
+
+
+_rel_perm1w_prime = ft.partial(rel_perm1w_prime)
+_rel_perm1nw_prime = ft.partial(rel_perm1nw_prime)
+_rel_perm2w_prime = ft.partial(rel_perm2w_prime)
+_rel_perm2nw_prime = ft.partial(rel_perm2nw_prime)
+_rel_perm3w_prime = ft.partial(rel_perm3w_prime)
+_rel_perm3nw_prime = ft.partial(rel_perm3nw_prime)
+_rel_perm4w_prime = ft.partial(rel_perm4w_prime)
+_rel_perm4nw_prime = ft.partial(rel_perm4nw_prime)
+_rel_perm5w_prime = ft.partial(rel_perm5w_prime)
+_rel_perm5nw_prime = ft.partial(rel_perm5nw_prime)
+_rel_perm6w_prime = ft.partial(rel_perm6w_prime)
+_rel_perm6nw_prime = ft.partial(rel_perm6nw_prime)
+
+subdomain1_rel_perm_prime = {
+ 'wetting': _rel_perm1w_prime,
+ 'nonwetting': _rel_perm1nw_prime
+}
+
+
+subdomain2_rel_perm_prime = {
+ 'wetting': _rel_perm2w_prime,
+ 'nonwetting': _rel_perm2nw_prime
+}
+
+subdomain3_rel_perm_prime = {
+ 'wetting': _rel_perm3w_prime,
+ 'nonwetting': _rel_perm3nw_prime
+}
+
+
+subdomain4_rel_perm_prime = {
+ 'wetting': _rel_perm4w_prime,
+ 'nonwetting': _rel_perm4nw_prime
+}
+
+subdomain5_rel_perm_prime = {
+ 'wetting': _rel_perm5w_prime,
+ 'nonwetting': _rel_perm5nw_prime
+}
+
+subdomain6_rel_perm_prime = {
+ 'wetting': _rel_perm6w_prime,
+ 'nonwetting': _rel_perm6nw_prime
+}
+
+
+# dictionary of relative permeabilties on all domains.
+ka_prime = {
+ 1: subdomain1_rel_perm_prime,
+ 2: subdomain2_rel_perm_prime,
+ 3: subdomain3_rel_perm_prime,
+ 4: subdomain4_rel_perm_prime,
+ 5: subdomain5_rel_perm_prime,
+ 6: subdomain6_rel_perm_prime,
+}
+
+
+
+# S-pc-relation ship. We use the van Genuchten approach, i.e.
+# pc = 1/alpha*(S^{-1/m} -1)^1/n, where we set alpha = 0, assume
+# m = 1-1/n (see Helmig) and assume that residual saturation is Sw
+# this function needs to be monotonically decreasing in the capillary pressure
+# pc.
+# Since in the richards case pc=-pw, this becomes as a function of pw a mono
+# tonically INCREASING function like in our Richards-Richards paper. However
+# since we unify the treatment in the code for Richards and two-phase, we need
+# the same requierment
+# for both cases, two-phase and Richards.
+# def saturation(pc, n_index, alpha):
+# # inverse capillary pressure-saturation-relationship
+# return df.conditional(pc > 0, 1/((1 + (alpha*pc)**n_index)**((n_index - 1)/n_index)), 1)
+#
+# # S-pc-relation ship. We use the van Genuchten approach, i.e. pc = 1/alpha*(S^{-1/m} -1)^1/n, where
+# # we set alpha = 0, assume m = 1-1/n (see Helmig) and assume that residual saturation is Sw
+# def saturation_sym(pc, n_index, alpha):
+# # inverse capillary pressure-saturation-relationship
+# #df.conditional(pc > 0,
+# return 1/((1 + (alpha*pc)**n_index)**((n_index - 1)/n_index))
+#
+#
+# # derivative of S-pc relationship with respect to pc. This is needed for the
+# # construction of a analytic solution.
+# def saturation_sym_prime(pc, n_index, alpha):
+# # inverse capillary pressure-saturation-relationship
+# return -(alpha*(n_index - 1)*(alpha*pc)**(n_index - 1)) / ( (1 + (alpha*pc)**n_index)**((2*n_index - 1)/n_index) )
+##
+# # note that the conditional definition of S-pc in the nonsymbolic part will be
+# # incorporated in the construction of the exact solution below.
+# S_pc_sym = {
+# 1: ft.partial(saturation_sym, n_index=3, alpha=0.001),
+# 2: ft.partial(saturation_sym, n_index=3, alpha=0.001),
+# 3: ft.partial(saturation_sym, n_index=3, alpha=0.001),
+# 4: ft.partial(saturation_sym, n_index=3, alpha=0.001),
+# 5: ft.partial(saturation_sym, n_index=3, alpha=0.001),
+# 6: ft.partial(saturation_sym, n_index=3, alpha=0.001)
+# }
+#
+# S_pc_sym_prime = {
+# 1: ft.partial(saturation_sym_prime, n_index=3, alpha=0.001),
+# 2: ft.partial(saturation_sym_prime, n_index=3, alpha=0.001),
+# 3: ft.partial(saturation_sym_prime, n_index=3, alpha=0.001),
+# 4: ft.partial(saturation_sym_prime, n_index=3, alpha=0.001),
+# 5: ft.partial(saturation_sym_prime, n_index=3, alpha=0.001),
+# 6: ft.partial(saturation_sym_prime, n_index=3, alpha=0.001)
+# }
+#
+# sat_pressure_relationship = {
+# 1: ft.partial(saturation, n_index=3, alpha=0.001),
+# 2: ft.partial(saturation, n_index=3, alpha=0.001),
+# 3: ft.partial(saturation, n_index=3, alpha=0.001),
+# 4: ft.partial(saturation, n_index=3, alpha=0.001),
+# 5: ft.partial(saturation, n_index=3, alpha=0.001),
+# 6: ft.partial(saturation, n_index=3, alpha=0.001)
+# }
+
+def saturation(pc, n_index):
+ # inverse capillary pressure-saturation-relationship
+ return df.conditional(pc > 0, 1/((1 + pc)**(1/(n_index + 1))), 1)
+
+
+def saturation_sym(pc, n_index):
+ # inverse capillary pressure-saturation-relationship
+ return 1/((1 + pc)**(1/(n_index + 1)))
+
+
+def saturation_sym_prime(pc, n_index):
+ # inverse capillary pressure-saturation-relationship
+ return -1/((n_index+1)*(1 + pc)**((n_index+2)/(n_index+1)))
+
+
+S_pc_sym = {
+ 1: ft.partial(saturation_sym, n_index=1),
+ 2: ft.partial(saturation_sym, n_index=1),
+ 3: ft.partial(saturation_sym, n_index=2),
+ 4: ft.partial(saturation_sym, n_index=2),
+ 5: ft.partial(saturation_sym, n_index=2),
+ 6: ft.partial(saturation_sym, n_index=2)
+}
+
+S_pc_sym_prime = {
+ 1: ft.partial(saturation_sym_prime, n_index=1),
+ 2: ft.partial(saturation_sym_prime, n_index=1),
+ 3: ft.partial(saturation_sym_prime, n_index=2),
+ 4: ft.partial(saturation_sym_prime, n_index=2),
+ 5: ft.partial(saturation_sym_prime, n_index=2),
+ 6: ft.partial(saturation_sym_prime, n_index=2)
+}
+
+sat_pressure_relationship = {
+ 1: ft.partial(saturation, n_index=1),
+ 2: ft.partial(saturation, n_index=1),
+ 3: ft.partial(saturation, n_index=2),
+ 4: ft.partial(saturation, n_index=2),
+ 5: ft.partial(saturation, n_index=2),
+ 6: ft.partial(saturation, n_index=2)
+}
+
+#############################################
+# Manufacture source expressions with sympy #
+#############################################
+x, y = sym.symbols('x[0], x[1]') # needed by UFL
+t = sym.symbols('t', positive=True)
+
+
+p_e_sym = {
+ 1: {'wetting': (-7.0 - (1.0 + t*t)*(1.0 + x*x + y*y)),
+ 'nonwetting': 0.0*t },
+ 2: {'wetting': (-7.0 - (1.0 + t*t)*(1.0 + x*x + y*y)),
+ 'nonwetting': 0.0*t },
+ 3: {'wetting': (-7.0 - (1.0 + t*t)*(1.0 + x*x)),
+ 'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2)*y**2 },
+ 4: {'wetting': (-7.0 - (1.0 + t*t)*(1.0 + x*x)),
+ 'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2)*y**2 },
+ 5: {'wetting': (-7.0 - (1.0 + t*t)*(1.0 + x*x)),
+ 'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2)*y**2 },
+ 6: {'wetting': (-7.0 - (1.0 + t*t)*(1.0 + x*x)),
+ 'nonwetting': (-1 -t*(1 + x**2) - sym.sqrt(2+t**2)*(1+y)*y**2)*y**2 },
+}
+
+pc_e_sym = dict()
+for subdomain, isR in isRichards.items():
+ if isR:
+ pc_e_sym.update({subdomain: -p_e_sym[subdomain]['wetting']})
+ else:
+ pc_e_sym.update(
+ {subdomain: p_e_sym[subdomain]['nonwetting']
+ - p_e_sym[subdomain]['wetting']}
+ )
+
+
+symbols = {"x": x,
+ "y": y,
+ "t": t}
+# turn above symbolic code into exact solution for dolphin and
+# construct the rhs that matches the above exact solution.
+exact_solution_example = hlp.generate_exact_solution_expressions(
+ symbols=symbols,
+ isRichards=isRichards,
+ symbolic_pressure=p_e_sym,
+ symbolic_capillary_pressure=pc_e_sym,
+ saturation_pressure_relationship=S_pc_sym,
+ saturation_pressure_relationship_prime=S_pc_sym_prime,
+ viscosity=viscosity,
+ porosity=porosity,
+ relative_permeability=relative_permeability,
+ relative_permeability_prime=ka_prime,
+ densities=densities,
+ gravity_acceleration=gravity_acceleration,
+ include_gravity=include_gravity,
+ )
+source_expression = exact_solution_example['source']
+exact_solution = exact_solution_example['exact_solution']
+initial_condition = exact_solution_example['initial_condition']
+
+# Dictionary of dirichlet boundary conditions.
+dirichletBC = dict()
+# similarly to the outer boundary dictionary, if a patch has no outer boundary
+# None should be written instead of an expression.
+# This is a bit of a brainfuck:
+# dirichletBC[ind] gives a dictionary of the outer boundaries of subdomain ind.
+# Since a domain patch can have several disjoint outer boundary parts, the
+# expressions need to get an enumaration index which starts at 0.
+# So dirichletBC[ind][j] is the dictionary of outer dirichlet conditions of
+# subdomain ind and boundary part j.
+# Finally, dirichletBC[ind][j]['wetting'] and dirichletBC[ind][j]['nonwetting']
+# return the actual expression needed for the dirichlet condition for both
+# phases if present.
+
+# subdomain index: {outer boudary part index: {phase: expression}}
+for subdomain in isRichards.keys():
+ # if subdomain has no outer boundary, outer_boundary_def_points[subdomain]
+ # is None
+ if outer_boundary_def_points[subdomain] is None:
+ dirichletBC.update({subdomain: None})
+ else:
+ dirichletBC.update({subdomain: dict()})
+ # set the dirichlet conditions to be the same code as exact solution on
+ # the subdomain.
+ for outer_boundary_ind in outer_boundary_def_points[subdomain].keys():
+ dirichletBC[subdomain].update(
+ {outer_boundary_ind: exact_solution[subdomain]}
+ )
+
+
+# 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()
+
+
+for starttime in starttimes:
+ for mesh_resolution, solver_tol in resolutions.items():
+ # initialise LDD simulation class
+ simulation = ldd.LDDsimulation(
+ tol=1E-14,
+ LDDsolver_tol=solver_tol,
+ debug=debugflag,
+ max_iter_num=max_iter_num,
+ FEM_Lagrange_degree=FEM_Lagrange_degree,
+ 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.initialise()
+ output_dir = simulation.output_dir
+ # simulation.write_exact_solution_to_xdmf()
+ output = simulation.run(analyse_condition=analyse_condition)
+ for subdomain_index, subdomain_output in output.items():
+ mesh_h = subdomain_output['mesh_size']
+ 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 norm_type, errornorm in error_dict.items():
+ data_dict.update(
+ {norm_type: errornorm}
+ )
+ errors = pd.DataFrame(data_dict, index=[mesh_resolution])
+ # check if file exists
+ 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
+ )
+ else:
+ errors.to_csv(
+ filename,
+ sep='\t',
+ 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:
+[33minitialise meshes and marker functions ...
+[0m
+[33minitialise interfaces ...
+[0m
+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.
+[33minitialise subdomains ...
+[0m
+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.
+[33mcreating xdmf files for each subdomain ...
+[0m
+[33mcreate initial values ...[0m
+[33mcreate source Expressions ...[0m
+[33mwriting exact solutions and/or source terms for all time steps to xdmf files ...[0m
+[32mStart time stepping[0m
+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