From d6fb7600b97c7c22c7c65f1c31545ad845b954cb Mon Sep 17 00:00:00 2001 From: danieldouglas92 Date: Tue, 4 Jun 2024 13:10:27 -0600 Subject: [PATCH 1/3] Add cookbook for water cyclying in a fixed 2D slab --- .../fixed_tian_slab/coupled-fixed-slab.prm | 262 ++++++++++++++++++ .../fixed_tian_slab/uncoupled-fixed-slab.prm | 257 +++++++++++++++++ 2 files changed, 519 insertions(+) create mode 100644 cookbooks/fixed_tian_slab/coupled-fixed-slab.prm create mode 100644 cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm diff --git a/cookbooks/fixed_tian_slab/coupled-fixed-slab.prm b/cookbooks/fixed_tian_slab/coupled-fixed-slab.prm new file mode 100644 index 00000000000..1602f2420c3 --- /dev/null +++ b/cookbooks/fixed_tian_slab/coupled-fixed-slab.prm @@ -0,0 +1,262 @@ +# This is a cookbook which simulates a simple 2D subduction system with +# a hydrated sediment layer, MORB layer, gabbro layer, and peridotite +# layer. The layers dehydrate as they are subducted, with the solid-fluid +# reactions governed by an approximation published by Tian et al., 2018. +set World builder file = /home/900358141/2024_aspect_hack/fixed_slab_cookbook/2D_tian_slab.wb +set Adiabatic surface temperature = 1600 +set Nonlinear solver scheme = iterated Advection and Stokes +set Output directory = tian_fixed_slab +set Max nonlinear iterations = 50 +set Nonlinear solver tolerance = 1e-5 +set Dimension = 2 +set End time = 1e6 +set Surface pressure = 0 +# If this is larger than the melt generation then that means that the fluid will instantly dehydrate +# Perhaps this is unstable, especially if the solver tolerance is 1? Make this 10x less than the +# reaction time step in the volatile material model. +set Maximum time step = 1e3 +set Use operator splitting = true +set Resume computation = auto + +subsection Solver parameters + subsection Stokes solver parameters + # set Stokes solver type = block GMG + set GMRES solver restart length = 150 + set Number of cheap Stokes solver steps = 2000 + set Use full A block as preconditioner = true + set Linear solver tolerance = 1e-8 + set Maximum number of expensive Stokes solver steps = 100 + end + subsection Operator splitting parameters + set Reaction time step = 1e2 + set Reaction time steps per advection step = 10 + end +end + +subsection Discretization + # We choose relatively large values for the stabilization parameters: + # # However, note that in an application model with a higher resolution, + # # we would choose much smaller values for the stabilization parameters. + subsection Stabilization parameters + set beta = 0.5 + set cR = 1 + end +end + +subsection Checkpointing + set Steps between checkpoint = 100 + # set Steps between checkpoint = 0 +end + +subsection Compositional fields + set Number of fields = 7 + set Names of fields = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone + set Compositional field methods = field, field, static, static, static, static, static +end + +subsection Temperature field + # set Temperature method = field + set Temperature method = static +end + +subsection Initial composition model + set List of model names = world builder + subsection World builder + set List of relevant compositions = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone + end +end + +subsection Boundary composition model + set Fixed composition boundary indicators = left, top + set List of model names = initial composition +end + +subsection Geometry model + set Model name = box + subsection Box + set X extent = 1000e3 + set Y extent = 500e3 + set Y repetitions = 1 + set X repetitions = 2 + end +end + + +subsection Gravity model + set Model name = function + subsection Function + set Variable names = x,y + set Function expression = 0; -9.81 + end +end + + +subsection Initial temperature model + set Model name = world builder +end + +subsection Boundary temperature model + set Fixed temperature boundary indicators = top, bottom, left, right + set List of model names = initial temperature +end + + +subsection Material model + + set Model name = reactive fluid transport + # set Model name = visco plastic + set Material averaging = harmonic average only viscosity + + subsection Reactive Fluid Transport Model + set Base model = visco plastic + set Reference fluid density = 1000 + set Shear to bulk viscosity ratio = 0.1 + set Reference fluid viscosity = 1 + set Reference permeability = 1e-6 + set Exponential fluid weakening factor = 27 + set Fluid compressibility = 0 + set Fluid reaction time scale for operator splitting = 5e4 + set Fluid-solid reaction scheme = tian approximation + set Maximum weight percent water in peridotite = 2 + set Maximum weight percent water in gabbro = 1 + set Maximum weight percent water in MORB = 2 + set Maximum weight percent water in sediment = 3 + end + + subsection Visco Plastic + set Viscosity averaging scheme = harmonic + set Viscous flow law = composite + set Prefactors for diffusion creep = background:4.5e-15, \ + porosity:4.5e-15, \ + bound_fluid:4.5e-15, \ + peridotite:4.5e-15, \ + gabbro:4.5e-15, \ + MORB:4.5e-15, \ + sediment:4.5e-15, \ + weak_zone:5e-19 + + set Stress exponents for diffusion creep = 1.0 + + set Activation volumes for diffusion creep = background:8.2e-6, \ + porosity:8.2e-6, \ + bound_fluid:8.2e-6, \ + peridotite:8.2e-6, \ + gabbro:8.2e-6, \ + MORB:8.2e-6, \ + sediment:8.2e-6, \ + weak_zone:8.2e-6 + + set Activation energies for diffusion creep = background:375e3, \ + porosity:375e3, \ + bound_fluid:375e3, \ + peridotite:375e3, \ + gabbro:375e3, \ + MORB:375e3, \ + sediment:375e3, \ + weak_zone:375e3 + + set Prefactors for dislocation creep = background:7.4e-15, \ + porosity:7.4e-15, \ + bound_fluid:7.4e-15, \ + peridotite:7.4e-15, \ + gabbro:7.4e-15, \ + MORB:7.4e-15, \ + sediment:7.4e-15, \ + weak_zone:1e-50 + + set Stress exponents for dislocation creep = 3.5 + + set Activation volumes for dislocation creep = background:14e-6, \ + porosity:14e-6, \ + bound_fluid:14e-6, \ + peridotite:14e-6, \ + gabbro:14e-6, \ + MORB:14e-6, \ + sediment:14e-6, \ + weak_zone:14e-6 + + set Activation energies for dislocation creep = background:530e3, \ + porosity:530e3, \ + bound_fluid:530e3, \ + peridotite:530e3, \ + gabbro:530e3, \ + MORB:530e3, \ + sediment:530e3, \ + weak_zone:530e3 + + set Angles of internal friction = 0 + + set Cohesions = background:1e50, \ + porosity:1e50, \ + bound_fluid:1e50, \ + peridotite:1e50, \ + gabbro:1e50, \ + MORB:1e50, \ + sediment:1e50, \ + weak_zone:1e6 + + set Minimum viscosity = 5e19 + set Maximum viscosity = 5e23 + # set Use plastic damper = true + # set Plastic damper viscosity = 1e21 + + # set Strain weakening mechanism = plastic weakening with plastic strain only + # set Friction strain weakening factors = 0.5 + # set Cohesion strain weakening factors = 0.5 + # set Start plasticity strain weakening intervals = 0 + # set End plasticity strain weakening intervals = 1 + # set Thermal expansivities = 0 + end +end + +subsection Mesh refinement + set Coarsening fraction = 0.05 + set Refinement fraction = 0.8 + set Initial adaptive refinement = 2 + set Initial global refinement = 3 + set Strategy = isosurfaces, composition threshold, minimum refinement function + set Time steps between mesh refinement = 2 + subsection Isosurfaces + set Isosurfaces = max, max, bound_fluid: 0.005|0.04 + end + + # minimum of 5 global refinements + subsection Minimum refinement function + set Function expression = 3 + end + + # refine where the porosity is bigger than 1e-6 + subsection Composition threshold + set Compositional field thresholds = 1e-6, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0 + end +end + +subsection Boundary velocity model + set Tangential velocity boundary indicators = top, bottom, right + # set Zero velocity boundary indicators = bottom +end + +subsection Boundary traction model + set Prescribed traction boundary indicators = left: initial lithostatic pressure + subsection Initial lithostatic pressure + set Representative point = 0, 0 + end +end + +subsection Melt settings + set Include melt transport = true +end + +subsection Postprocess + set List of postprocessors = visualization, composition statistics, velocity statistics + subsection Visualization + set List of output variables = density, viscosity, thermal expansivity, named additional outputs, melt material properties, melt fraction + subsection Melt material properties + set List of properties = fluid density, permeability, fluid viscosity, compaction viscosity + end + set Output format = vtu + set Time between graphical output = 10e3 + set Interpolate output = true + set Number of grouped files = 0 + end +end diff --git a/cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm b/cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm new file mode 100644 index 00000000000..d277a1ef2c6 --- /dev/null +++ b/cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm @@ -0,0 +1,257 @@ +# This is a cookbook which simulates a simple 2D subduction system with +# a hydrated sediment layer, MORB layer, gabbro layer, and peridotite +# layer. The layers dehydrate as they are subducted, with the solid-fluid +# reactions governed by an approximation published by Tian et al., 2018. +set World builder file = /home/900358141/2024_aspect_hack/fixed_slab_cookbook/2D_tian_slab.wb +set Adiabatic surface temperature = 1600 +set Nonlinear solver scheme = iterated Advection and Stokes +set Output directory = tian-slab-uncoupled +set Max nonlinear iterations = 50 +set Nonlinear solver tolerance = 1e-5 +set Dimension = 2 +set End time = 1e6 +set Surface pressure = 0 +# If this is larger than the melt generation then that means that the fluid will instantly dehydrate +# Perhaps this is unstable, especially if the solver tolerance is 1? Make this 10x less than the +# reaction time step in the volatile material model. +set Maximum time step = 1e3 +set Use operator splitting = true +set Resume computation = auto + +subsection Solver parameters + subsection Stokes solver parameters + # set Stokes solver type = block GMG + set GMRES solver restart length = 150 + set Number of cheap Stokes solver steps = 2000 + set Use full A block as preconditioner = true + set Linear solver tolerance = 1e-8 + set Maximum number of expensive Stokes solver steps = 100 + end + subsection Operator splitting parameters + set Reaction time step = 1e2 + set Reaction time steps per advection step = 10 + end +end + +subsection Discretization + # We choose relatively large values for the stabilization parameters: + # # However, note that in an application model with a higher resolution, + # # we would choose much smaller values for the stabilization parameters. + subsection Stabilization parameters + set beta = 0.5 + set cR = 1 + end +end + +subsection Checkpointing + set Steps between checkpoint = 100 + # set Steps between checkpoint = 0 +end + +subsection Compositional fields + set Number of fields = 7 + set Names of fields = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone + set Compositional field methods = darcy field, static, static, static, static, static, static + # set Compositional field methods = darcy field, field, field, field, field, field, field +end + +subsection Temperature field + set Temperature method = static +end + +subsection Initial composition model + set List of model names = world builder + subsection World builder + set List of relevant compositions = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone + end +end + +subsection Boundary composition model + set Fixed composition boundary indicators = left, top + set List of model names = initial composition +end + +subsection Geometry model + set Model name = box + subsection Box + set X extent = 1000e3 + set Y extent = 500e3 + set Y repetitions = 1 + set X repetitions = 2 + end +end + + +subsection Gravity model + set Model name = function + subsection Function + set Variable names = x,y + set Function expression = 0; -9.81 + end +end + + +subsection Initial temperature model + set Model name = world builder +end + +subsection Boundary temperature model + set Fixed temperature boundary indicators = top, bottom, left, right + set List of model names = initial temperature +end + + +subsection Material model + + set Model name = reactive fluid transport + set Material averaging = harmonic average only viscosity + + subsection Reactive Fluid Transport Model + set Base model = visco plastic + set Reference fluid density = 1000 + set Shear to bulk viscosity ratio = 0.1 + set Reference fluid viscosity = 1 + set Reference permeability = 1e-6 + set Exponential fluid weakening factor = 27 + set Fluid compressibility = 0 + set Fluid reaction time scale for operator splitting = 5e4 + set Fluid-solid reaction scheme = tian approximation + set Maximum weight percent water in peridotite = 2 + set Maximum weight percent water in gabbro = 1 + set Maximum weight percent water in MORB = 2 + set Maximum weight percent water in sediment = 3 + end + + subsection Visco Plastic + set Viscosity averaging scheme = harmonic + set Viscous flow law = composite + set Prefactors for diffusion creep = background:4.5e-15, \ + porosity:4.5e-15, \ + bound_fluid:4.5e-15, \ + peridotite:4.5e-15, \ + gabbro:4.5e-15, \ + MORB:4.5e-15, \ + sediment:4.5e-15, \ + weak_zone:4.5e-15 + + set Stress exponents for diffusion creep = 1.0 + + set Activation volumes for diffusion creep = background:8.2e-6, \ + porosity:8.2e-6, \ + bound_fluid:8.2e-6, \ + peridotite:8.2e-6, \ + gabbro:8.2e-6, \ + MORB:8.2e-6, \ + sediment:8.2e-6, \ + weak_zone:8.2e-6 + + set Activation energies for diffusion creep = background:375e3, \ + porosity:375e3, \ + bound_fluid:375e3, \ + peridotite:375e3, \ + gabbro:375e3, \ + MORB:375e3, \ + sediment:375e3, \ + weak_zone:375e3 + + set Prefactors for dislocation creep = background:7.4e-15, \ + porosity:7.4e-15, \ + bound_fluid:7.4e-15, \ + peridotite:7.4e-15, \ + gabbro:7.4e-15, \ + MORB:7.4e-15, \ + sediment:7.4e-15, \ + weak_zone:1e-50 + + set Stress exponents for dislocation creep = 3.5 + + set Activation volumes for dislocation creep = background:14e-6, \ + porosity:14e-6, \ + bound_fluid:14e-6, \ + peridotite:14e-6, \ + gabbro:14e-6, \ + MORB:14e-6, \ + sediment:14e-6, \ + weak_zone:14e-6 + + set Activation energies for dislocation creep = background:530e3, \ + porosity:530e3, \ + bound_fluid:530e3, \ + peridotite:530e3, \ + gabbro:530e3, \ + MORB:530e3, \ + sediment:530e3, \ + weak_zone:530e3 + + set Angles of internal friction = 0 + + set Cohesions = background:1e50, \ + porosity:1e50, \ + bound_fluid:1e50, \ + peridotite:1e50, \ + gabbro:1e50, \ + MORB:1e50, \ + sediment:1e50, \ + weak_zone:1e6 + + set Minimum viscosity = 5e19 + set Maximum viscosity = 5e23 + # set Use plastic damper = true + # set Plastic damper viscosity = 1e21 + + # set Strain weakening mechanism = plastic weakening with plastic strain only + # set Friction strain weakening factors = 0.5 + # set Cohesion strain weakening factors = 0.5 + # set Start plasticity strain weakening intervals = 0 + # set End plasticity strain weakening intervals = 1 + # set Thermal expansivities = 0 + end +end + +subsection Mesh refinement + set Coarsening fraction = 0.05 + set Refinement fraction = 0.8 + set Initial adaptive refinement = 2 + set Initial global refinement = 4 + set Strategy = isosurfaces, composition threshold, minimum refinement function + set Time steps between mesh refinement = 2 + subsection Isosurfaces + set Isosurfaces = max, max, bound_fluid: 0.005|0.04 + end + + # minimum of 4 global refinements + subsection Minimum refinement function + set Function expression = 4 + end + + # refine where the porosity is bigger than 1e-6 + subsection Composition threshold + set Compositional field thresholds = 1e-6, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0 + end +end + +subsection Boundary velocity model + set Tangential velocity boundary indicators = top, bottom, right, left +end + +# subsection Boundary traction model +# set Prescribed traction boundary indicators = left: initial lithostatic pressure +# subsection Initial lithostatic pressure +# set Representative point = 0, 0 +# end +# end + +subsection Melt settings + set Include melt transport = false +end + +subsection Postprocess + set List of postprocessors = visualization, composition statistics, velocity statistics + subsection Visualization + set List of output variables = density, viscosity, thermal expansivity, named additional outputs + set Output format = vtu + set Time between graphical output = 0 + set Interpolate output = true + set Number of grouped files = 0 + end +end From deeea8f3ff7c4111a94ee80b14c54c49b215bcaf Mon Sep 17 00:00:00 2001 From: danieldouglas92 Date: Tue, 4 Jun 2024 13:54:48 -0600 Subject: [PATCH 2/3] fix stuff --- .../fixed_tian_slab/uncoupled-fixed-slab.prm | 47 ++++++++++--------- 1 file changed, 25 insertions(+), 22 deletions(-) diff --git a/cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm b/cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm index d277a1ef2c6..748393f6339 100644 --- a/cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm +++ b/cookbooks/fixed_tian_slab/uncoupled-fixed-slab.prm @@ -49,9 +49,9 @@ subsection Checkpointing end subsection Compositional fields - set Number of fields = 7 - set Names of fields = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone - set Compositional field methods = darcy field, static, static, static, static, static, static + set Number of fields = 8 + set Names of fields = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone1, weak_zone2 + set Compositional field methods = darcy field, static, static, static, static, static, static, static # set Compositional field methods = darcy field, field, field, field, field, field, field end @@ -60,10 +60,14 @@ subsection Temperature field end subsection Initial composition model - set List of model names = world builder + set List of model names = world builder, function subsection World builder - set List of relevant compositions = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone + set List of relevant compositions = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone1 end + subsection Function + set Coordinate system = cartesian + set Variable names = x, y + set Function expression = 0;0;0;0;0;0;0;if(x<=50e3 & y>=450e3, 1, 0) end subsection Boundary composition model @@ -76,8 +80,8 @@ subsection Geometry model subsection Box set X extent = 1000e3 set Y extent = 500e3 - set Y repetitions = 1 - set X repetitions = 2 + set Y repetitions = 10 + set X repetitions = 5 end end @@ -132,7 +136,8 @@ subsection Material model gabbro:4.5e-15, \ MORB:4.5e-15, \ sediment:4.5e-15, \ - weak_zone:4.5e-15 + weak_zone1:4.5e-15, \ + weak_zone2:4.5e-15 set Stress exponents for diffusion creep = 1.0 @@ -143,7 +148,8 @@ subsection Material model gabbro:8.2e-6, \ MORB:8.2e-6, \ sediment:8.2e-6, \ - weak_zone:8.2e-6 + weak_zone1:8.2e-6, \ + weak_zone2:8.2e-6 set Activation energies for diffusion creep = background:375e3, \ porosity:375e3, \ @@ -152,7 +158,7 @@ subsection Material model gabbro:375e3, \ MORB:375e3, \ sediment:375e3, \ - weak_zone:375e3 + weak_zone1:375e3 set Prefactors for dislocation creep = background:7.4e-15, \ porosity:7.4e-15, \ @@ -161,7 +167,8 @@ subsection Material model gabbro:7.4e-15, \ MORB:7.4e-15, \ sediment:7.4e-15, \ - weak_zone:1e-50 + weak_zone1:1e-50, \ + weak_zone2:1e-50 set Stress exponents for dislocation creep = 3.5 @@ -172,7 +179,8 @@ subsection Material model gabbro:14e-6, \ MORB:14e-6, \ sediment:14e-6, \ - weak_zone:14e-6 + weak_zone1:14e-6, \ + weak_zone2:14e-6 set Activation energies for dislocation creep = background:530e3, \ porosity:530e3, \ @@ -181,7 +189,8 @@ subsection Material model gabbro:530e3, \ MORB:530e3, \ sediment:530e3, \ - weak_zone:530e3 + weak_zone1:530e3, \ + weak_zone2:530e3 set Angles of internal friction = 0 @@ -192,7 +201,8 @@ subsection Material model gabbro:1e50, \ MORB:1e50, \ sediment:1e50, \ - weak_zone:1e6 + weak_zone1:1e6, \ + weak_zone2:1e6 set Minimum viscosity = 5e19 set Maximum viscosity = 5e23 @@ -226,7 +236,7 @@ subsection Mesh refinement # refine where the porosity is bigger than 1e-6 subsection Composition threshold - set Compositional field thresholds = 1e-6, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0 + set Compositional field thresholds = 1e-6, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0 end end @@ -234,13 +244,6 @@ subsection Boundary velocity model set Tangential velocity boundary indicators = top, bottom, right, left end -# subsection Boundary traction model -# set Prescribed traction boundary indicators = left: initial lithostatic pressure -# subsection Initial lithostatic pressure -# set Representative point = 0, 0 -# end -# end - subsection Melt settings set Include melt transport = false end From 1f22ea92f89f4a1ac083cd42b2b35bd03657f13a Mon Sep 17 00:00:00 2001 From: danieldouglas92 Date: Tue, 4 Jun 2024 22:53:32 -0600 Subject: [PATCH 3/3] Add dynamic 2D two-phase fluid subduction cookbook --- .../2D_tian_slab.wb | 80 ++++++ .../doc/dynamic_subduction_with_two_phases.md | 12 + .../doc/fixed_slab_temperature_fluid.png | Bin 0 -> 121188 bytes .../uncoupled-fixed-slab.prm | 252 ++++++++++++++++++ .../user/cookbooks/geophysical-setups.md | 1 + 5 files changed, 345 insertions(+) create mode 100644 cookbooks/dynamic_subduction_with_two_phases/2D_tian_slab.wb create mode 100644 cookbooks/dynamic_subduction_with_two_phases/doc/dynamic_subduction_with_two_phases.md create mode 100644 cookbooks/dynamic_subduction_with_two_phases/doc/fixed_slab_temperature_fluid.png create mode 100644 cookbooks/dynamic_subduction_with_two_phases/uncoupled-fixed-slab.prm diff --git a/cookbooks/dynamic_subduction_with_two_phases/2D_tian_slab.wb b/cookbooks/dynamic_subduction_with_two_phases/2D_tian_slab.wb new file mode 100644 index 00000000000..44c8c08ede5 --- /dev/null +++ b/cookbooks/dynamic_subduction_with_two_phases/2D_tian_slab.wb @@ -0,0 +1,80 @@ +// Composition 1 is bound fluid +// Composition 2 is peridotite +// Composition 3 is gabbro +// Composition 4 is MORB +// Composition 5 is sediment +// Composition 6 is weak zone +{ + "version": "0.6", + "gravity model":{"model":"uniform", "magnitude":10}, + "cross section":[[0,0],[50e3,0]], + "surface temperature":273, "potential mantle temperature":1573, + "thermal expansion coefficient":3.1e-5, "specific heat":1000, "thermal diffusivity":1.0e-6, + "features": + [ + {"model": "mantle layer", "name": "peridotite mantle", "coordinates": [[-500e3, 100e3], [-500e3, -100e3], [2000e3, -100e3], [2000e3, 100e3]], + "min depth": 0, "max depth":10000e3, + "composition models": + [{"model": "uniform", "min depth":0.0, "max depth":10000e3, "compositions": [2]}]}, + + {"model": "oceanic plate", "name": "Overriding Plate", + "coordinates": [[300e3, -100e3], [300e3, 100e3], [2000e3, 100e3], [2000e3, -100e3]], + "min depth": 0, + "max depth": 100e3, + "composition models": [{"model": "uniform", "compositions": [2], "min depth":0, "max depth": 100e3}], + "temperature models": [{"model": "half space model", "ridge coordinates": [[[1500e3,-100e3],[1500e3,100e3]]], + "spreading velocity": 0.05, "min depth": 0, "max depth":100e3}] + }, + + {"model": "oceanic plate", "name": "Subducting Plate", + "coordinates": [[300e3, -100e3], [300e3, 100e3], [-500e3, 100e3], [-500e3, -100e3]], + "min depth": 0, + "max depth": 150e3, + "composition models": + [ + {"model": "uniform", "compositions": [1], "min depth":0.0, "max depth": 10e3, "fractions": [0.02]}, + {"model": "uniform", "compositions": [1], "min depth":10e3, "max depth": 20e3, "fractions": [0.01]}, + {"model": "uniform", "compositions": [1], "min depth":20e3, "max depth": 30e3, "fractions": [0.005]}, + {"model": "uniform", "compositions": [5], "min depth":0.0, "max depth": 10e3, "operation":"add"}, + {"model": "uniform", "compositions": [4], "min depth":10e3, "max depth": 20e3, "operation":"add"}, + {"model": "uniform", "compositions": [3], "min depth":20e3, "max depth": 30e3, "operation":"add"}], + + "temperature models": [{"model": "half space model", "ridge coordinates": [[[-500e3,-100e3],[-500e3,100e3]]], + "spreading velocity": 0.05, "min depth": 0, "max depth": 150e3, "bottom temperature": -1}] + }, + {"model":"subducting plate", "name":"Interface", + + "coordinates":[[310e3,-100e3],[310e3,100e3]], + "dip point":[1e7,0],"max depth":50e3, + "segments":[{"length":150e3,"thickness":[150e3],"top truncation":[-50e3],"angle":[0,45]}], + + "composition models":[ + {"model":"uniform", "compositions":[6], "min distance slab top":0, "max distance slab top":20e3}] + }, + + {"model":"subducting plate", "name":"Slab", + + "coordinates":[[300e3,-100e3],[300e3,100e3]], + "dip point":[1e7,0],"max depth":1000e3, + "segments":[{"length":150e3,"thickness":[150e3],"top truncation":[-50e3],"angle":[0,45]}], + + "composition models":[ + {"model":"uniform", "compositions":[1], "min distance slab top":0, "max distance slab top":10e3, "fractions": [0.02]}, + {"model":"uniform", "compositions":[1], "min distance slab top":10e3, "max distance slab top":20e3, "fractions": [0.01]}, + {"model":"uniform", "compositions":[1], "min distance slab top":20e3, "max distance slab top":30e3, "fractions": [0.005]}, + {"model":"uniform", "compositions":[5], "min distance slab top":0, "max distance slab top":10e3, "operation":"add"}, + {"model":"uniform", "compositions":[4], "min distance slab top":10e3, "max distance slab top":20e3, "operation":"add"}, + {"model":"uniform", "compositions":[3], "min distance slab top":20e3, "max distance slab top":30e3, "operation":"add"}], + + "temperature models":[{"model":"mass conserving", + "reference model name": "half space model", + "density":3300, "thermal conductivity":3.3, "adiabatic heating":true, + "plate velocity":0.05, + "ridge coordinates":[[[-500e3,-100e3],[-500e3,100e3]]], + "coupling depth":80e3, + "forearc cooling factor":5.0, + "taper distance":0, + "min distance slab top":-200e3, "max distance slab top":300e3}] + } + ] +} diff --git a/cookbooks/dynamic_subduction_with_two_phases/doc/dynamic_subduction_with_two_phases.md b/cookbooks/dynamic_subduction_with_two_phases/doc/dynamic_subduction_with_two_phases.md new file mode 100644 index 00000000000..f7d55495609 --- /dev/null +++ b/cookbooks/dynamic_subduction_with_two_phases/doc/dynamic_subduction_with_two_phases.md @@ -0,0 +1,12 @@ +(sec:cookbooks:dynamic_subduction_with_two_phases)= +# Dynamic 2D Subduction Model with Reactive Two Phase Fluid Flow + +*This section was contributed by Daniel Douglas.* + +In this cookbook we extend the simplified kinematic subduction model built in INSERT_LINK_HERE to a more dynamic model of subduction in a 2D Cartesian box. The model begins with a slab that has just started to subduct with a leading slab dip of 45 degrees, and a weak plate interface to facilitate subduction. The initial conditions for this model were constructed using the Geodynamic World Builder, see this [tutorial](https://gwb.readthedocs.io/en/latest/user_manual/cookbooks/simple_subduction_2d_cartesian/doc/README.html) for an in depth explanation on how to setup your own World Builder file, and the file used to generate this model can be found [here](https://www.github.com/geodynamics/aspect/blob/main/cookbooks/dynamic_subduction_with_two_phases/WORLDBUILDER_FILE.WB). While we still kinematically impose a convergence velocity of 5cm/yr, this is only prescribed on the upper left boundary. We use the 'box with lithosphere boundary indicators' geometry model to only impose this velocity on the shallowest 100 km of the left boundary, and for depths larger than 100 km the boundary in left open to allow inflow/outflow. On the right boundary, we impose a free slip boundary condition on the shallowest 100 km, and an open boundary for depths greater than 100 km. For the top and bottom boundaries, a free slip boundary condition is imposed. With the exception of the entire left boundary, the composition is not held fixed. We fix the temperature on the boundaries for the top, bottom, and the shallowest 100 km on both the left and right boundaires, and allow it to dynamically evolve for depths greater than 100 km on the left and right boundaries. The subducting plate (slab and unsubducted oceanic lithosphere) has a layered composition comprised of (from shallowest to deepest): 10 km of sediment containing 2 wt% water, 10 km of mid-ocean ridge basalt (MORB) hydrated to 1 wt% water, and 10 km of gabbro hydrated to 0.5 wt% water. The rest of the model is composed of peridotite, and there is no free fluid at the beginning of the model. + +```{figure-md} fig:initial-temperature-hydration + + +Schematic diagram showing the model region of a subduction zone, and the actual schematic diagram of the model. +``` diff --git a/cookbooks/dynamic_subduction_with_two_phases/doc/fixed_slab_temperature_fluid.png b/cookbooks/dynamic_subduction_with_two_phases/doc/fixed_slab_temperature_fluid.png new file mode 100644 index 0000000000000000000000000000000000000000..670bc29bf0e25ff9e115469d462d5a25437e8df8 GIT binary patch literal 121188 zcmeFZgx(;s-JO&VR1u!;Xmlb4kS&>h_<@mRJA(&>+zNcQ-_B1~}-c zv{bw7>p#tio+3(_%RfE`{tM>8B0>f#FpP76KyOI_n|fE~pcWl`xo9SQy5~7-p6S2r zNtxv_F$zBuF4Px2FK|E>OpgSA4>HHjc_2rw$VF&FRm3ra82mee#VN)V-*)#`4c#6z zH~kb4XgJA?i5P|jA6K69wQ_MqT2mE0)qS1+@fG&_3#@(Ig6(wMzJJmN_u4#V&)2JG z|BiAY;dwf32N?i?-o;p-O!qqld!te}^NE;VltBMTcw`dUW+Be@ui-&HR0>*DpI_77 zQzm0|CsT5N)U8a*uAI@~UfD_q0&U?#@^!U2ixS6FxJs*&!Q57f0==x$-ovNr*QScL|2;_HE|7f?5v3 zR~0($!;SUz9qx!w-d%QG9vJN{iz ziF%Z3FrN!f%fd-TitO&WrYyo~C}hNM)sNMthH@r10?g2FvgIofZuMqVz7Ea2X9DT!-YhAia7+ik;n^W8>##2g-8hpN6&HC3Rbh#|X6buLI+a888L_Sf2k^H^CFERW*{9+*F z^rrgjLDVXX%Nk2(#m%Dh+p}8N01Ky-JMsl2QJpWVA2EjdE+llS8|0h5druz9?(fe{ ztZx1fICtR@wBUz`LOoxBZKpa~e4h_9>^nN^5ggE7L5RPzcSz8%dkM4)>=Z51MO6yb z2RlI$Csa6-ywmA)%XHCxW&Ru|=4eDFdbNKO5`)#wou{yYaO)Z&IflAL41ah-w+wz= z`ljTUA{BEq1^9kx*o55ShMU>j!N;#-e*4ixZ-pRYQH3F~2@r#|(AT2oX@8IZWoM+@ z_j1!?M5Y6?an)QebOOws2Nmd|nt!kxYo$7{bP{qq>SShIrrWpf=ujagT=mL`Lyv3R z1Hu)GNEgcP2D^diz{r1|T}Rw;3bhJ@wn5vx6edXdo1R6Jjku1PuZ`_(X~JL7XBr4{ z)#~7^j-d*oFJ471Erlt0nuMnYfzr?2@(T*FBWencn0rXt3_4plIXdHF<$5D zdk%^Y_5~(IT<7^(F)3!5BhHFfBzJac-WLIS1{U9i?;frN%^^N?#J4@!)Gt@Sa*1}oFX0N6L$fNMeXUV z9H&;(U}A7fmLeuZJ8kWnnDKYZ(QBIi$SPi4 zex`{&gl7jwNmPClA=vZ!IR7jG5|T{M-qXvUz@ogu5)sd0b~W=emDJ$8!vd{~1$vj$ zj=|NZU#_+{dA@1Q~CA4jq1@CqBmRc!+J-DLPFW-tl$L#(Jac>EVs|@MRjLe=3+J;r z9+JIk-$ceF9o9=5Bi4Se8Jy(pn#MbxMJ`q>Ho|VF9_rCnVG9E*c!zp{vE?bF6zL=B zRF;V}XSRd^|7kVFxY(}G8nXPZHO*=BCKTvJyG-fi4pObB^kJm0pW5v;ig?b@Pi#nP zTpV7|*~i5wET_0#FYb32GGrKfC6uoz$%+g51$g&$yX|b(kD3SNWu~5n@rly?CA9B? z|2ohVU3onk))!Zs-P!EZMLmR^$rDg?B?e+SIdw8Qt)T8%z;=fX5FJNTAsa9&SkI{bGxO4 zs0Pd79dl=2(HstUKOWGpnsyX1&f4_PN&d$s5$^{W8wwY1cp zaE-o3psJteWPR&{pX_@$^0^RW49r{UD)=i1^nFO&Pt?VlXzipmVbrp~GWfC?`>-f|fNi;ue`y^^8jUfH0h1IrC*b?T3gpB+L76oSHH!!R#~V%eC9whU)jqC)HJL zG<>rbM_sZv?KARmV~$Suwoa(+`eoO9gf^w~0yn6WMm0`P>vimEH1&6KF~gOJC^T_8 zsRMV=9zY^C4cQ-n=1)9jy%r2m$DbW^K=1Cbq@YN&dlTS^a>|?#U2P-Q@rcWkG%5Gb z($VAq2Yu3MtDi{P9l%&Vi{{T3%XNjxV@y5zw_SpXJ|3f!mliO2t42LZ|I3Fq@jT|2 zh~MzkD1H;)!KxinQn~Qo74Jvd+#f++klinxys33j(*7MOCn8?wKYDaT$@jSVZ1Q|h zF5@L@f`3AD)ZJ&BC{%WEqR8EnxDp(TsfrcAzxKo`efgVC-b=y})u*Ec`hUQ68Y|#s zv9`y{Qymp;UZI;^l-b#2Iv1piT+-438#c|3aRW(AwbMJnp12?nAW;ODkaZJ1CNEfif-iNVbuYdjCX!$=~05 zZ)DM)t}LJE?iPXeqY?mX-eI!5cX#B84_T4+s?~z6+2#pq6eG5^it0l)A1uPaLKA=97L24r^9_Oj-wH5B5J>*R6FBKw9xu|6o0<~`ja2V3Ou4*H)~Tvj4Kdaj z0K)B*V1ie)UA}fVR&63!`9&M>rhqW-;DaOwA;223o=I+u)&7nfDQj1>=BALtBmZ$I zMcllW+6HwnaDKr^bJNzzq_I+%k(J@~M7)kam6p6S<4yDz?W^ZA=OBHdtn;qQ)Skux z@fqQn;a*3Z#IeXQAYZW1{Az{XMysZax0WAo1ISM0%tN@C_QWjso78&$(LjpJ(^;=d zZOZXaJd@T?#vOGdr(Q?XWXHj7*=~BycZny6Sl0bg)x;j~s@M56Cje^iq0l3uk`7Wv zAb33|TFl?Famba3rF2XUJQYzI6Mr zebxS{SE$5??OK#S7EeA0ES=^_Al4L#6#req2B&QQCHr>tL7-Pby;6u|LR+nK{X1Cf z8|RV1l`Ik)1e3ovYPz`(xK z0L<#0Y1b&r?KXqiafb2`KBH|~tIuX}h$!?3x~)C)LOzykAvOP7Tp28oeWo@wOZw<% z?RcE-{lzI>?iy+%t{GMFfc3BO)y}}O=br(m{N4O)$P<*1b!*NaQlnR_IxtzF3XJwE zCFOm&9R6Z{M_-s*lN?#4!$;$Sf=&GLn^MGYxtS-wY*f=D?TZrQ!MywbRAy8`;|l}7 z+yjB0p@Efyo;piB*h`7o}3-rPlc=4cXTJ;W5o ze{ALFCU1@3+FhpJ<>4@D%vNZd`F}t*;ZMn9Dh{krT}95~ZA>y95J{#SF9HCb_w@c{ zpc9V^>FQGr$MrH{xY)wUBJl}xJwtU{an$EI{$ODj@33jQ*o}7xMd4q_+ z%S);0cO_#lpp1ROMlRw0IumJ?;FQtR0|u1W41@~AKIxJd9S`Oz&DXg8ejhj>I17ZV z^7a`XC&a{mXyo;hHWzaW1Ah}5kxwl58xqq02~7b4{h`#YulXiYwfi#)I#m5=bFK%` zKRUjuG}3j+#C}+IvqQhL6st1mik?c}5DK9*yWUAemZ{(yY5iuyW!G<747}9{+BJ1D z`v7NxB!X2+XJU+v)aP+P76xy*mF_mRV#pj$3p5*stX5(mYJWBNq5t;fPHT zi%npOT^%SN)>q+zBqnNemL`s0?6A0e%alnv`G=;Qtcu^+2%LaJ6jd0 z^_3&j#{$P^o-)UuxHadR!LnyLL-keTPq?`YX>#(Xx&~JRx)-*WhKEde@Wq!UBCNu! z!h`{Eo+K$F2N#;o3p@I|q9Jaho&@eZVPxYQ5%a-+tP1n@FPIK4W=}SILBA%4Y42F5 zVzdbhB6D_%F%q=?iC10iE&#N3UxtboYUv@gcEUXg}3M*GuzZH10~2e5ZS z{|?B1JBt1hgzlETtsLMvSa|ny9A7lVXx!1^<1Ez936T}KbQ1CGpQa@&1NLithbrT| zg3Ze2BX4ox6p?UNw%u>h5cgGGT?fEe+5Zk`SL8SVJEfmKu~~|GlK7mX{Z(jAG!u0v zQ>+4Q9|whdhIb+FDNWAQS!#Cs>pIE0K&CHo(D(|=qY9^lVZL%q3>xfq`wID4GGPNG zZtBza5ny@!$CWe0DEc5<8}bCNj0~$vjhQDn6Fsc{eFF`nO{IVbH^gF6g^N#O6k0ra zFFPv2nWz<+2xl41?!7<5ns%N(T;SxEgo`!Ved7Pn7xs!rk%^+mA+J@koD9i#c-cfE z`(^x7$v+CEjQ=m5-QKJL^ZmtH47GBe7=hu&Cg#!=K_pN$uN#xQDV+ekAiRU2u`Q2` zO4ntO;znjOyW-xT?65vBM=b%zYwiI*GWm~<@8V4=n<&*u?^huw(Pi-cgZ1*R|f3HMnWh8Q|p%1nIA2{1f;J;2fnzEes!Ql%`FIgosStZ&WmHMcv!Z|*&F~jXMJ)!TQ z<`$Wk@q+EIQ(tfZsqmBt=Ah(%T%Maa ztkXZfM&Xv9D}BljD)!Gqr1JM(*$y?(R+7oAkpvDCEc)4u>~A`potYdbq4}h=RcPjH zc#aSFI3Z)&oSaFzV6MDm@I_=KFF7M7a!+FHe$Q+E=UP}0_F?tjezT9BXG!1n%RLIc zQgRuzjo!TTwVA^=3B=KhH-6?i_fnVmg9TD~@&{aF{4M#edu#iO>Z@B#TOJ+8DzhfV zYpWLJYpYO0lIc~8Sc|KG$zNf#Tu{Dh%&JQ9oZ(=PaVoU7X)~t1>Hr2lJ!n*_38Y3a z;n@)PJ|%n2k@r%uBAl6|TQKs)%yL;UGi_|%tIy0I`mWPFCNwYNa2Q}B+ZW|};=;Ru zHxx!c0sb_Rmb;LMp;QV`M(*?KVxN1ejXIpX5N#EH-f4Z6p7rVhOBWz!wf?21ro>?t z0}cD#t@a05vhZF8hIp^VS@Led*cVh0fcLSqzYKj1%ZYmOCHM)q@aBYD;XLG!@a>%x zp&a`Txp2g^CVt4Ou*2#+3pXJB{q|7TrvdH`EdboTXe02()6|^3too#@is#r&*#Q!* z3w9`;O@+0aMIPg3&d398ru?-658cS|P3{PmR{;TV06$yt_s8k&zRimZKK=hLPHBCJ-t!!I*9HWiWm zd+U_g37yrF`#gGOndgO~FwgOqeZAWaj<6%EMP4krRs&ou@5Z`Jg{Um~`V$oVE}7 zVX{Rf;?L)W#@+qdM8_dS(~wAl|F(k&xMC&X&BGqrpLc^R|2wO^9v2lU;7*$GHgG)r zy1t-VF+18NuxEIekh(s1=TYDt5ZxDk*p9q|P}8e&kr71z{mg3Rg{9_XPnNnD9z}a2 zzqFSMMxDzuAOoAG!Y3AYte;5&mVcE0C9W*QRC^0mPAx4h>9<`z{$O$hM^tdGuPI`_ zrDf`>avg%ZpGyR4xj0irba@{E3YG=`Y%tNJqF-U${IT8S74!VzH{)E#nBnGn)??9g z_Uk8#N6TDGsYjw^YLU~`6}&p5p3~;O0o`_{T&0z`QMb+6=gIbY)p;?hdtEIv{&#H| zf#hWtFk{~hKv(*s?}YgM-qN=D@_Z520G-{1>y?_-Ayt3YK(U(@jusL;-`iV%koC8! zEWdUtULSKl7+D0ulfy?Rsx)fGl}=Q z{huuAKX>^y?;;fv$xd2v%ay&Cpm#kSCvJb+t&f#RNtDO({mwgeFu)VwNzX7GsuS<7 zlW-x1#Sj9)?}`dQiK_=gOd#nhlXju7Xks}o>a|9gtgp+{g zDRX~$&hg>d)qCffJ)gtGq|#D`Mc$8kEjFx=&JFRt?$O*w-;zsKEM9W;VVh_SG%z6h z(*kk`&pg~#9Yr?{hw_Om1AdGHOu`M57kr}CWU)|9gX;UvY~AegJvcSX@WOJk9O}wa zrzMqSt6-~;5ZJG6oeF-g=SUjj7gaQtl>N`mcVPjLp8GHSI$|cywN=$Y$}Dcc`%2!F z05YmX1SV#}8c25{>W5WBbm)x1I1Bj66iN^65h7(5Cw!`nqHIaJDj3Y%?_L8a$U7p8 z%(z9X7qle|&M|@Jy`14=iqKNEncL~Wcwm+8kkwFxhil?czge)>Q=>VAjA?#xOci>q zs==f@rnu`pq(+Qm_QOMy|2&)g|yMFSw|17C{_0 zo#**x*o9q-@TLGdV*EFbP5xfTaw+b%E#v&gBXgD`dRZAX8@4kYa)oS|HF>bu$FT%s z<=n{I&;$p9GToDK89STu{g%V`lW%Z%jW(~5uQZxJ0_hR>W7pNf+ zgE#WBm08rWRoswx74FickzTPb!_pXA(MOov=Wkux>3f_{Ck~wnS1NVS`U6|HfXnex z7lolx?Q!t(sk#2RI`q&ch`trK;a%T#x}5dLhk&|4Hex%;T4?K;f~v_#a|s#u6$$u$ z3iC=C5p3o!JOe4cc7axOHP~``)3C0C&~JdH#vWE0R#eekzY~8bI#Nw8$^Nvn% zZ_FU~yF@o0?tHJ^pZ&kE)?$rwGC90u(YCQwX`kUI6mO6(RW)#XQmtBX%akKC`Ozsl zRR>Bk+=aI5!+yfGVwMw#p~--chD|9xJQZaE6ZVqM#+ATq_DKH1a**}$fFb^G?R%u1 z45;n#rn$p_suja7K%*y70>U=geOdcj50ZN$4|O57Znxeurrxa}@;63N3SaOT@GG@4v*93ItqeT zW=-KJN*qz-hG;ZaN&gD$JfZva0K+TJ$kYSX$A{-`oXlu9n#6V$vxrI~q|?RNywxn# zExjyf+_?5NnnaF!_+fEU!_HTNYH5Yj)yHNu-|zB;0z41FO3A|Pr}di^nXzSgFWX3L z_0yuLCo7`@S>KbLPOgo;Z=2Rv2PqwbyI}s!8quf9=>XB&=J>TGda`&k!aiav63@bo zk|l{CR=$4C^yn)~MEI(h`1=}d+!&8ua9CR_eeozSnW^soM1NNv zU(VlY&<^PhkX`N9{dujnt~PtpUd`>@Z9#rSieJx?t%@%$5Fcmm@mzcF#^#Q#Eo#pE zicZ@6XQ-=No6qqxxNJ&*@uY23me;jtFBV8bF#!9QGi;eAAC0Vrh`hdu>+9S7V6ZN*5NwTn)-CSc`V2Awrz zWu`h%?BdDuT*~UhPf;Te8MaT@gJB4uFuZ3IQ1!`;xV3cC zdWZ@+$Y4JDcw6gQ|6K7{$o$Ls6^;lfEvMz3R-`_sjh$d^i|j-Iwy7|HCVs zoU@>9T)`tqv-kGJZq8QKLzi1)V4i{flG=;@Yas2pJJ`#9GRE1-3NoyDjCp zGv|I}DYf)i9`DAa!pXt|!G8da^g*cG*_irdMuNylkb5I^2+}uRAaHd7+idnR7j~g| z=5rJX@*e(G@Me`S(`~ui9l8CP-uIA3BR*#@tLpDPEF)rivsGAwbnb91MJ%NE7^@T|CaUm; z6eW5{b@$Z`C>9$zUGbXZE!V|*@wPtx7xvD^xU)#tg}VNST6+;+NNjb_ktwFuy%RX^ zx3_#vAB0$TJb2+?0HOjHG0H@)JdIr4*JC^@e7(Q)xjZ;kHY={mb!~d8s=|pm&#Iec zf<>Rb;C5nN_JgSAIL);z%1_dOL{1?V+~4hBaHcCu@rE&!c$k}ygfvf?kNISpC52pL z((#bmEx>$ZkO9gddUv)}3AS)}EZxy#>;$RS2RoF;n1UTk^{rjOKr-G^gs}{h^Yn2eS9yDEcxJT*vb{{K`jIQH$uTkHal-P zqciJ~1mdNwNoMv|AbT>e7MQ%t^J7*cB)*}Ew!?mx3gQLxacXiJdI_YLg#WaOep=i$ zPo5MEUTRgMP+3=jfXjGXESMn9xJ8P`tTYRtDbEE&dDdH$P}gDbud{iE!NpsP6rZ5^ z^jP*|4G|B<%**tgy`|y@U2}#QN+R#J5OLu``O~f$ed<+SPYX}cYLWaxb~1cku8L3R zdU!yaE_s*|GMioOMF9^ZMn;sl8=wpKBp~BksNBatB<2>FXC=BM^wS}eGQDM4o|sw` zFiG%^}M%oO zp1br1{4$8)-H#OyqIPiJLS&9p%n6-?KBl4~max7wCpaV0-Cjo0l{MA~o!TuP!3OdM z24w(hz1Ah`$NvRIB!t|gtohJ8=Hj)<7*bOk-H5V9dh;E??5edF=R zVd>339HG_E@qH!4X~GrR0H;rgcyK8-EJSc!9#SsR@I|DqW>`L~`s?7i%elyeVYgVB z{04@tVNkkR4XN)h`b&+w2s)uTIeEtc$j-}807+dYFN zqVBpi48{PqrbhfR?! zD#YjQYNFIH=U!zq7Q#i@FzJd+=A>G?wr7wT~zxb1Ke{#aAEN9i8_rBhPq9$lkpOR1G$XSBEO z*34len(||&!{UF>uC}j79$GJ{hRn-OT<5>L)I=piUiF9N(z6%FVIKzNalEbj-+*<` z*8H+cUdiM_`FsK_Ubea}#lKoG05GWUllgTSntFFECOJq{skuowftLd{C&^=KN9xDa zRMcv?n#c~PX5lgh20BiEvpAXZ*i~Ni@Kc==2LAQtQ^Vs41_xR9_Kr^xF^@>u*eKXk z-txT7_+hxkqE#`qt}8&R9dc|q?GN7I)6^e57{e~KJ?79zUi^4Y_!idytn`{J)csfa zQA_r%*X@B!M1_kpDPDyyhA!K@L4F8GxxO!& zzAkf+J(BWyoBxeQvI=@==ULV-p*5jMO;gY2p}vX#<$o<>p-N~*j(g2acdH^pgT{HHFg}hyUE{Gv2hGwnSeZB39VTiwHGIYX@5SMsdCEy2~dFqQv4qb1`GN?94~kh++n<3vX{Tygv{4^}xh zdz=v{u;~w*m;)-U8dM$mQ3F7h+pL=1W!?CN7fSlkrKU?Zmor!@} z8_IY)yXF9GV?Wn>BVKf%nQw}9P^??bS;ys4=!uNdq)kVN)q3^u0%N$F#T#bdJn~)y z;?e4=W=^M(Xm10s4o}y;=_$|L44QSJsuXiFHHyh5gF>}>)R-*7d2Ny+Yo}vygq^|p zeWge@Tm7E)J^qaqRkJT-l2TZsgZ;9rDXx`Or{vpZ_=}dfB2wEyhzmj2nlG$=SRkAS z%9H`r0N?*LBguz!EI1j^p}7>aytLq~8<>>2yr}aiDQ~;<*T*^Lh+P$nNU#HpneN#Qm5VY1nCr(o}L4{1`#|xR?58_yIafxT+Ra(J+oFLHDfEr0;Suxl z?Y&FudD#a8Mk!5SS%5B6t&GqBuh)o{7e_h-b?z>x{$V%tuyCn{>MoA%zlxCltZ(#P zczj?!E7-dxeKB7?KK50(JAhZ34F%5r)0K^@C3T0lXS!c%kD^e_+6XlO_*g z98e|9^XlDjG7p|e#r4R;wkuh}W2E2;E?Qy8bd$PLfK1@szS5<64JFaveb3;08S!2g zJwTvD5KFwW@LART_JY*s#-0zuXuspRniN+js^It?#>&F|*)OD~TmA*Qr><83m_s+_ zuIMTPhEn!iQOVtRwPHhwd)^N(nHqRC0{T?t@W^SDBfMm@Z#pb&6m>=5V340K_79t0 z3IP*A6Oda`0UQ(`BLq<4@KYmagEoGj7$jKFOU=uO2obIV+;TPOw6t92dhg?IR;3}& z+S&0f_05wq;SXTupCTU|nN3O%pX2nnp})x4lKc&eu#>J7jqRc-G9KQ{U@I0U9X75i zd197hY65~gzXIaC$3>%1@C}@M6el3(&1p*iuu(Id_PY;z&A#g$N;ww0Y!aYYTV13f zwQb5UFn{}OMmKW!_QmSA>Y`HWblP>=zN8syFhkt10?>2()bMX%jmnU`@ubLnHl-Xn z@ZHQ6`gow4JCvA7Bc21F;MnS?*I9Jcpv-l~rmb|k=z8>~UsifgU6*rBNon-zxc{r7 z&%igA(9`vak8!;8Sq`+BZggq=GV5nI%?)pIrM<6OBa(ckauUM$3h7!9bdl{Mb|XKkdjx^svQXX_`dK!ksg~c ztI@J4GD&m)LG?pV3d%6z(y6Oq+mMvfM20@UdOoF{k9Kq~X{dFzfy*(-1iOmQj47td zk#;jV*~sex(3=vM3g~tBO6E=tf`3HulTO!_ffTwXgO4G=**#%7Kpi^QQyf+*V<{PASJ?8XFsZg~3jeAC%nVS$Zym#V!W| zY*9!!Qm4Bi7fs-641lm0gS84PPQ@DHrSSMjVcL*xiPUieAW5>rh`li*-adTei zRR;qLyL00^&C)1)riO`+I4>t0O#Kf;SHg<+Ll1==(Du@U*A>d~Yt@Jcj4vMXcOI|b z&Fb+r^JWq9R%G0b_KP|1^`%lQoC}8b1m&sp94<&+C%Mrhk9B!n8a=T{T4&i`m$736 z%8J^1R!AaS35(5QO=RBCc%!XkAph+A-33adF&_INK4ZK|jr+M%ya}62pv%GB&FR81 z+yXvr1n8B9ld)=LF!^EqbF>V5*?5+fRQJqk#q8ViYfO=p;qT%u;rmc#WCFO|^imd8 z3Zu_`_b@d!z<(w4y0oiYLD}bmm5L`8-%J3>O4FZDs2GZc;)WXFYZKMF7}rxls!w{vhXA5`CkNdg-HKX;Q#_-h(~~ zORl}p?Bk<2lId9#IqZqoLHm@eOGqb;rSfM6j9w>94`$C$r}7xqN-}I+y~$4S-X3g- z-K47h47&yBJ|Qk<>RDd{H|9mn3dU!48xI=L$CAN_-@pNx?p=QhqZ^4{iQXk+&+_kE zcG`WNqOlAOH1R8DCn6ar|H&O+p8(bmz?c?3M zTH|S#w=lg@J{lrDx8I67I8Gf0raViqz34ejM=(}U5I^DjHzbO|oG#j{KB%1O<1>-8p9^1NKE2{X*P5H63QYSK zb1#8&T^oKiX0^;2v@jiG9-=;joeoJWUB^368=^%N^E*E^$@ndqx-ZjIi#}7HpS`WD z%A)JD!=FWE0sdwoO_#-j1v0cHUF_=Zh=?W4#7n|UdYSYx6^}FjN9l94i}awdS#$` zuGRfx!FG55DVXjtzjf1rv-$(T_|Jb&nx1a})>pj4w>OFllGY)l7*VWIb35yTp`E-L!#+EfN6*LHumd5A+==cy8VhgO1{<^{Vmw?k zJb|h-po7hEWT!UWO0*>(seH{K>iy^l+w2?YaLtm&+)Xl^fopr3*rX(5_^+WpeEw(-x(@s zZVO=QjSMqBT{R27t^3%rJ6jQ?3IgNz3y?(UuWnqxsve^H!1Z%1F*>cfmt|A8v_S6w@a2aDXiN=N6h-FqRo z)U7%MzJV5Ns5m$))&A4L8SVXk+3P^zMMI&$+I@W#g*5+?Q}tnBnfo}PxMD|2)*kKo zG!A7ul~ThyFJ@}J_TZxpZ(eTT{8#{LC51VJaC_!JE2YLuy5@XMncl;)YliFAhxEQ{ z!_4cw6lOSZR#U>uXJ|=b;L#TWYZ0UN7-4g_pr<6e1zihoplHv}@LWe7pyg8Jk7O>P z{+(LtbaAn(`Adw1{wrwL(dw*u{!!E!!?`Q@(fq(-M6XC$hDDY~b+Xi%%^4k##WEK2&rQ4H_s{`_Z6XX%m{;aGtQfih^>?wV&R>y}s;hNs;6cs?ROmAFC3*Vlf($Oeg?t1LD=C}AQ#8i6k zO+ut=MQ;f5x1E%qWqD2+%^aN;PVPIWVxUP=0|Eii76(TR?P-fF` z8}UcyYWoF~Y^!vrGknH9lkFbSpg)626LfF(e_OLN5Ix8cC>o0mb@(X?I=k=4|7(&M z{&|tx?4oV>YvZn1tW}Bd1eH^?fWNp|?a&`p*Z>9a6yPsI^V9TXCtuobWu^n;pkYp_(t+XN>_m9(a+dYNAjZCByulVxi)EUS_2lpc)wszCa60Lq7G#ugdDO%Z4Y7(aUsx2qP^;@ zbvIFw+^lep2==I$+9r`6uYI3gk2|$hkRBk`5`gokwVK}3GoRcA z=`#gBY|aiG_B^dGSv7Quflg30RWBKHmqN8woW{>5y<9m12^f2m!1da<$IMD_^exP4 z{1z(GlJ%ms+w4@9hCzjf+PBOR)xWm5F);ExtyIx-y|(65w6CR>nK1WhWMO&8i`}zh zi(KN^&5scpHG@-+^orxqjvp^;ZGNLU3LQ!noDa0^FoQ0WhfiBknW5rYo5w{ zuuOS=g9=ZG=i+&e0Cl|z9%Gwzm23idH)6D_dSNe+>!+wUSs0M*>V_UJ7q0NRF-*Iz z?^_24rn^~1*Gzz&XuCH&rlZW?qWn`*xW_%}!C)Q+Nf84uPf8yU4OI4Ue%(2dZg5Rt zw=3KJo2qzc3S`&Lo&xX3bRN(e2P|XYin9Zu_Rm?Ak1u%p*p$(NS`LHq+R*RTA~YIn>Nu z&2=3TvwGcaGNhcY;xn=Ka)Jicy!CT6$w;|~tP^-#>Ws2Y)Vl+xuRmfQ^a%$$Fx+LX zHA%k~he4x*!f{69C>qT_bOqM-=W#h=!K$WIF<*Cc_qg+rE$Shl1!Yd%?M(bz+HCS2 zc_3mbijRohI6S+^91>nY?6A~^ozE1Dc@hlycQbq@5&$2K6&RTL*t92<*bPYp=U%^K zV!Z=AIQ8$}1sKtGWKDtVj?lYJ18pP_p!F_?1ql?a&V(@&ELu)qAhIz0+I;yzj}B-X zBoA4c7)grxamqDGu8+gda9ho;>HmD3jp&aVN044i%4x04Kd6rnn>ner1NBb)b(^$q z`YCOiymv&2bWj>@i_Z&=Cni}s2SZY(ks~^;SDwb@aWOaM7Y9-rBR^%%&_>y%H3wM- z({@y!^oT>l>mG$ZspNfO7?9MEsl8pk%t*iU`dYn zl`d^f3rL8teQ$nr9x8Qxn&BB7A|)*>J+F7Jt7kAdavL;S9(?{bo6Gj6|Ls}Pp#W?e zJ{hy&X^UGE!~ZCnq*F*M-9qjD1fKF;AE-~E>EDOyT0lUn3B`*?a7N%+pl5ei0HspJ ztA7W2uN>|&q^B_#oV1IVAhx+y`z>FQ#bx1pzub25-B=>mlKz`L&&jLyOf-nSCSlft zzrty;bbe^@2PHP`VcaW=`s|=`lVKEZKYLU(_yhHP6!h1sb!IPa1<)yQ2#(t&O(9)UD`$rG}8)Vy3>$yW4+IwIUeH#Q3?sn3jZ1T7)=J0Q9ETA@$d@bVgF8w9!%i51i` zJ>@tRReFQxum***q+YzvsPP`&{uJxEF{6|vjc6c!QCNG)r5+{X)+Wc%r|?*bYIW=7 zUcqC?WAd-gZ;0rWgcXI^+6K{%3ZWkyol(0LyzT-^va|F}NS6HC`b;^nB2wLsnEzM3 z#^*{ZtH+@ZlydGo9$_L*0fS>F`n|Rno$YvoZ$z9Ot&(Ld9gSkBqa739s8F@u10o~v zkxV*u3@{`RkI}Euz^}wATg#<+tH+$O?!isRSE%P-It(AN;3eO?ABp{3pk1I-pq-@s zt8#`xs%&xDB@;=QyX-@`?PdBmuJ(yS(OQYaKqI_U-typ(RO&}DK5yxjh7RxO?rfH_ z^hWqKO$(*N*`w!Y54 zpU!SiN{dej=2et*sVP#w|Am<+=qVOPWdr>xc}{Cycil)cnoJn)X>e06t?B4qL%QfC|BAqil;(mI?j8na@v6(ledWb@!bP@1d5NL18UBZu zsemwoFOd6@aI9tuZiey zoL_bE7#+6TMqVhUP@6tBVN;iKqW3?f$*TyQ5gvcx;xHKPL?N!z#(e@7=7=SIfJW;*~Ovu zw(YVPQ)?*R*$|~6;{0WTH;AEBz~$6stOWsPJUnI<`d;==ZZo^u4nHwP7Jk;QhM zNV+x#Gw!Er&UObM^C?}SdK&QPjuuj8|VyLr=*>GsJ!EbkbkM*Gaebi34y4O6suGH;X69}728r2ohpS9sKIyx1`kO&jvMY8?IPB~JRU zV>}Z|)z-QU@CylUt$aA%rM(r2IBAR~?P9w`6AlZ9v;5|bTv zR6n}zmFZ)KlCN8nT?fOt@(u)XTbs{0ZeO)*Zn~-cst>l_T?}O5sg|uyt1~-@4G7n{ z>>FL$I82@q{YiCj*sR36lsu-oE`;7-*vR%0uk;R8w~EG+kLW|&*WQ|z-PZlJ zyg#)f*yZ`W4%DG;e<<+IEG)F-I9P$Q+5J4=Z&rwnQmsB775{&zddsk=+b>!e69nmQ z5h-awLIk9{Q@Xnbq!~u(MnJkG1tp|=ke2R_p@$y2W`OsH&;Ok3J(pibKFr*?W39FK zJ?+k(&TZ`5UEI5>C383Nx4_7e8U&FbgB%8U_}y7 zny=24jGp6F?r|EFO3&S)jYe+pw7#zc*Kr?O#KqJH?E^Utu)l25j7B~$`~n%UAKt{f zv^A7dzP4~z8t3g6T#EvAljTkNE3xF{o5$WHt1BfhJjIHq$dWv?u5GD4D<(bPqxH<7 zu_i)GqA6N&`xa1SLLS_1BWMJIwZXRL1_%=bdMH&~sko@Psnn?4sNPVy0-p*E8zER; z!zcXWe?1=-^3dn)k|ZSECJY;svX_7<^|N2rN*Itf(SYnA?E39jl!)-ZUu4TDT}9jy z``Arjo=kHWW>lbO^FZo@_L=X)3!@?6vXAh;7@t1sg+cio7nwlJf0}o!h{#dQxi*|hpr3gGa`ti%$5D-Pz~9c^LaYKZ;il+Wz7K|4!A&B6Il~xKn!Kg)7 z32J=-qob>d7uGiR2q1q(a?8-K-HaeUeCDE}IO$LLfnX5aYP&nm%C7-(B1miirt6DLLeF@=mE($ZGQ-=Y^mt6U!aGOMl-lEW-G z!s1r0+Q(+~C%=^Wm5me!O|J)UfUo)8AvK?t+T*TwF_=WyR<}-8yZqv^mOP`o!`~XB z`m8j&?dSt`3yybNQCckqVimZDhf5XQyVd0f?&3OvsdsF6h-Z3Qq;? zgKv@}x(crVSyhG1y3Tk~wc?MJ$seNMrm35@RAQ4iIX^!%Nh|=x9w_#voqEPy!QPet`m6EdAJVMWbO)FWUP-E=6`N)1SfvHuTDEuNK%q(^ zcSM03)HP3wt#L-ao!_ni%p=cNN?~>q)9-}qaPbuF`Qxp2KyLc*Uga1QzzseA+0RGa zWZE;Nt_$I8kd9jSw7xX+Ri9SZ(@ zlHr;*%XNIh3Fcy}hVbL}3@%-D#}zXi@+{y*ZBpqJ zXeV!X$ETkF9__!tQ{y1cX%ityiQnAvu#8pQQ#!)BP~;-3cK+gTpXs8_&ShT$e^eoR zvAws`MaU3P=Y;k-IPkub46QmCRI;1GpZNaIvfeRuv)-SE?jB5uB-4W$83>nmi#S%= z09UqL##%I!AP!j!SqyWGAPzZ>U$IcM**NMq2{n74^IGQfrz}rrNmyZdfxWw{1 z8t_gN85Y{O89bQt)iF);*E^2bZIie1?H<4d@&%Mv1$x?^Z2}f8zf=6_Xyp=`vY-MD zgg{R5yKL*(X!LN;*|LtGo#M6i6Ro$)OM3%pA;c7Yx<@ zX#s{+8&;WRUsoBYB{m9e&1apqyf$C|lvMk3{`-BM9Gc0C`-)^Vw0F^XLC+Y;@B=V{ z$3>EVGElWCh@FO^?B;$tEekYSj6ag_A*Wl{4j4(_&f$8=V%}pE^mX~Ku^4gYmEvR#YBBme8t>ivV|}YPiL9K7v*(^8`nBNgR`9s1 zLy&6Fr{Yb%gY!9ldM~Z+t*atJ=0#vAXk?Q2b`I^ck^MHF^84L0j9hQpkMY#aN8TxJ zZ0fB1x?Yp^G(+*n6ZQLAn=xx!zaB3ALPg^p2aqYSxYOWTujV+s&j^=xBh+aXq_96x zl+yE1xm?)y>#Nln8LMcffD+zhwCZESit_f6Ks%qtZ2^(qAzWC|A~4oBf`*CQWZO${ zzyHRf`w(!y1&EGTe_}Wa9xHq_7NgVWVn&B%!jdZJ>t=+<^w=~2qx({pTcwF(M6RDN zfKwegWOKX0B+S;2z8Q-6c}xCqsROkyBxAX1!`8Z%6e~^d2<1O!L9P!knC{-94l}p( zZlz0!4*A+W#rs`SD_$e)tT*Mf3e_7&0Cq*zCl01mOY#aX-n^2kqOypto0e}o6hj=n_s&9Wj|foPS#XL``C1i^#R)Ff)6bDMDYaA zzDP!`Vo2p+J|e#4Q8E3VizkIoz#adPB|I^$fz9G5cvMptj49Vs$aN9XW^+t5F=x~= zGy^qm*3Q9!3+8Nbz-O@Ma6%}-lZ=n;qbMlyf?fpWYZiIdcQQ00ymc@B z?(^@4g}3QsxQWfBJu$svsOfrkN-6(LI0^_jfezWI?=w@b`WY>d4v&0Y?Umoem7LGV!GTD27*wRt6va1;cm+Qk#ZT;!7``47~hQcZ>=b*Ff&N1(M7?Ax;Abweo)3j<^=)MjTK_<0P#~$ z2ucb?r%Ha}EoGx^wS%S8FOK9xMg#L8y`*)=Ls54e~NaCb{i732Bx-E60~^yNR;v`5S8u>Ju20Et1bHy z*ivWTy&o!jULfl^2H>C!B+$}jqKKf)Y|$7-n7IL#cu_yPx@8tANz69^H}ZcixWK*r z2L3%O6I#iKNPqv-@R?FDl2Wo8WI_=OxqYKD_S>*p{2N;yYp6JT-!LA^(*rgU8=k0 zf)^L(A|h#TY+#Qr5E$YFQm;rPF(|Ou9>>MYhR41y9=h8DuErzlO$9WmUN?EIM`B+g^&{rHQ?4+H6ZQc+ScYYx__p>B!_q3Yuo<{+E)lb0Sl{H?* z3kSCkI3@dWq!S&g+iLK2ut$C*d5fJ(sW}~9hCjp&4+0I4o0i-_7Zfi!boN|l{jY6_ z`{CV8jMP0XyF+UMFKz;WuS>H?I}+lYF}1@N$D~`|enOTo7WE}&`NEsOkbV1UOrApQ zP1|oZX$&bb%~TInNYsaIqgxS&i_j&C67rakH>0YIuXV>9+@j8?I?tw1`y_l?Vi-~+ zN8K@sUI^t|r%Lb-9fTkkPMM#Ow0vq|wW1|OL>*s7K25*JlJ{V~qtw!>w8%aI zkLx@A$tjmZpZ~7qcf_k64VHqh;_7q^aq9B&`LE*D<-P=|ss#n9s=r`~iwFLMQ~ur$ zlus1o%r!&Sk%8Y0Oxmr)UUE9zRJG%Kk`!BHbhIrka!{#?GoJUrtoXf%h~U5A_3H4l z3*jJ$lQVPj&X0x^aAiTok-O#ZEw)zZnzrN^3nY69>B8_7na!SSL|_Cn8jS0jW#SSJF;-^P{iCxRz@xio9v*_1A-=x>Zj+ zaI`&CPy_oj9{smh5!qb_Qw3mI)kywQ8~84 z0ohl6XI`dx8qGF&VlRDM5)Oh;Fi?_mdHdn~I|hH^#~O?>&LI?Jc+fUL#K!4_lc%%E zdk_R4Jp#F1u%-(8&Y^MuA@K(V%#TN3JgF-4Kt0a&<0R&i%!OaucPSPymK6U5=;JrP zhW#{W6RJ8qCyy7BfVXnD9 z^N-zgsQ(Pw(mcws#Q11kM}(lZ=|@V^Vj{$|F5C80Wy|VTdG4b9muIdOQu~V*^!81v zeHB>0_9h0>lsLoF3hGUpyB__fbU}%gG}XCu*!}OI>EeZTnS20(ECtwXdA4p@!Qonrvh@uTn>e5$a?c4+Qa36DQcq z!o9GUNj}p2(^fP$12wm*>2>cxEvWsA0nuI87$5|^n!{N=ndhG+UHg-9oyWT2k`6P~ zt39dJNAR-p)Xnr#N>i2_Cvv-92MpB(-SfM`Xa!`Sjqa-!fdfjz6}^=|xR!x6#E=&c z7`j|df|t+KWU^fJVGD4XzvBFoIL%^F%aae4$qjiL8`*SVki$Un0$hjB~ZeX=Q4h zSF+X-I~d=9(8tWTK^f=pAt}E%|}gvfG-w6|>HU+9xVe8KzdJTqIdTb+=u}4i7kZ_)8v! zuLwwZ&pYhv9>lNpx{#v|k@my5R_)KPeM*&Q>Twu!-0_@X6sB!WYJTax=R{P#R_c9;J=Kcf2eP*n zL30VRjp#3DvvL!cw{2Ps0||5MK0%2fx~my7JLOJ?*8By`wDH897+mCKqJJG+Jv{_Y zDS2UJN1_I>exCTc0H7N*u@?57?Yk#QXuRrtr_Q^vRGCc=)Ra+ zZhgJm>vuzMCf%nc=_ot{3hsEWfMUd}9uvObU%{)4xPT$<_-Lm>eJoE{kY?}G^q@KY zN@(Rt#hC7Ui!O;-F9E{nO~x{RH{K+bCIiQk^{i9ty1WSCnp9+-3^-&T;M^}$5VM=v zULOm1l_-0kb$kt>W9#=mV`-O9edIKI-gLV`%Py+hWW|u{jyquSxV(>y@y8j-#N(Bw z*bhSYX&DLaeg2L1;>n;Mn?3&cpOTtRpot`*)Nyz^W~lt5DJZs~k@qN;e^jazw`CZ| zO_RSI$zgCcqU4|zJ{c?Tn`?Bbe0q*HZ~%9P304}uGx|7lTSprmveP}=eOc``V7$1Y zdykE5iQgW4f4WaeXJKJ=+{jfbk!e^I+ z$0M3njSpgM7F+MoHV};S{FQ@8yBq|wb3BSafdc}hNpSq7_wkn1qf9{FTU7vtJ;H}3 z;=RVy#4Gv}YVyy()s9|dyQ|LuxIM6(L3r`vd!LYca)v60|u7yMA=bts3j1 zY>}>(?;GLO!bClvhM_sw_K_KKL3rwz&%WfUK7nFBdS#YRB<&?9?V=Rit0DW2XWToE zEQ(P@EE9X0!kVdMrfyi2E`DFJ8ZQE>Ia@6T8*&Mwc=^A4gt@@NQ-NyQ>zF8*dk3?S zyO5hDzqn$J2I2vmW$FQPk=OF@rQ-A&YR&zh#V?(K4LoweIdCBu6FA2ZjdSbA~=Db?7Ih33+PHYF)AyOX&uUG}- zbMnNEG2}j$3)uD0mH*w%k|*gQ49V+w(rkZ~{>Eza6XWWnCBgH64uOo)RPJ=3tnYLVEPfue7*1Of;Pg!wF7il zL)&^8H97zxGmYvkA;pl(A{_d~xwig9o{KJUez3)B(^|JvP9P zVAkp-G>_-*YT#uD&5;W$C}+OG(@_0MuUhgk-Hz4;s&4JwPZslOZyRu#$#ex8r9E~W z`qw8aT-FT3 zq$I~<7P00fb+C-ishoM?sI`aHT2t_%tD++sr_mt0PWQ%erRvOXBK*_7*BuS<#U{=D zg#H}gmp=e73m%!dI2S{hiQ-?ioQ(K&r4#~5W+$Ir;|nFCUPs(CV%#Q?^M`w(0%0zA zIb7{=A3#y`U(CMve#0bPydqUDGuYR;WA^>@(o6=BOR2D}m-Ig!B5!~ccv#?gQ4J~i zs<8zY@uQn~&_kGjJ{Kz|sG5#ZSw5;hNmf@|QX8?A7{t*c5>H%~E+~$vvARIGIyPfH z9B)SxfM|D(^X!3Z`tOYZ@^AEk7~^Tp>*A&#&ypC}=Dp{Wd}Enj4qJ)TvGRIqc#GCK zZ_k;ziZ^=9i38I@`;7A+QphF-ObDZph^eWIg)YNfnwx>yo568~tGbOCizt#_PB_Zy zUDN;L4`vFCiOE8)4{JvI)U-YJk6amId{xDT$^bR`4ta`=sZvJ5R$l@Q)!@4v+F1tJ zukkc1GCs{vU>A9IIZGx7Qr|CM*-XZ!K>gIC{X+jOz3-E zBPmn|)koQ`UQubT|KBV>Jt`w{k>|>jTub2|U@j!8vW)8@B$RUP7z{3}K^SO*S9e&X z*QmvsdD}fFzmoz^!8rka8hWlqsKa=kv=fZzd;?H2ZWBh^_o!6&F6kKbN&5LoGQIiYbrF_6ClPxLAPwD)uvSH9g{YndLOKA5XJP^epY^p44N z>OET34w(cOC8*!-ObVDe~1q+Pjy(M z*h-GIg0rxe9=mgzvaCe9&)QMy-4gkB#j4Cf-mU0GWg93N==vgSqf5@&vg){$FYv2@ z4c~SyMklkgY;t2=|gov$&S)_n4D&^x0<5$*S>(*A*5Ojn6+H^fJw&>jg@qJL+5A^;C?>n+OHf{HOA$=rf=EZ)oi>3-^*Ctz1i9!Qdi{qQQjNbSvY z$Z_UNjkUeL*^-Ex?mItK3A=q+kB*$)H2bLy#e^TtxSa3(deASz&GuV?NdWP|@JhDG zabWe^%+WEtj?rfPrlE*!a}+LC@n|3QqKU z@yD42kZ34W=jr_s#9;xNU=D#?0**8*RSaY|WI#mj$TEkottTIEkA4{|aor9Da*9@g ze=b#OD#G=olf~a9-44KEv@z9pA|D)>Jf~aE&cS(K>|x?7xt$MF3u#D8fOs;|3Q&pe zRn;vHVHjoef1E4tCz*lfwKU}$D(eGC2M`fZf&uP84fJ}Zs^5R5;ds1NQ6KC2>9_k` zb!_`)G8En&KD$OFqV}pz*hc&P-t`RDZ`#!t?P5Xg82=?95l>>83>h_2lIpj{B($wI zx7v$ZRZzcnIv!@#JA>Lig(;3WbP@rl$Q(6SIp{X{s}EvALyg!d=X3$TWdt-H~4nCpj1u?~0ucyd*TGFDwuB)h^Y7-mmv`S83ICP5Hf$3ih;G72 z6K9-x*B3l4ds(JILoq0LUd=CTBn{?4txt>^CPQy@_)GLtQ5At#^H^qV4DpzBn(WJ? zy8f$ab&}>+)e7`>zCC{1$~PZ+4clJt&H~_*?jC9}tGmb0E_k5Ra=#AI48NVa@x&D( z4%pry&I_`^CgIEnhYHoUeY*F3-#f0MjJRJpsGY*`2z2}{ye5C<6R!eDbWge z*)-~wY3{jpmH8)MpFay*G#U8-be#U|w0%M(0+WneaB<5?-xtc3p!IA0-Ere2@D;Xw zZv7xTyOp#Utgq(+FkBle6PE|{ZVBdefwGc@kT9{g5X!W7d{tlvhdu0a-$Ao=XV|mN zW<~%@t3SH7n&anu$K-#GL)$RUWCOH;%EdkL^g%*V)MIv4U~4~w`NfjpPLq0nDeF3^ z?WWw}d(iSCibzHCsBr6?KPWa0E6^D z+DmERfi;8YpiZ5DcPr>p`q2#3y@R54$JMq*5Fxbd+j`T|RF(5GsC-yIDt|%5vD*2c z#;%_t-(utf7+@>E$6UtrhrOpJ%VHkrXs5KfATygzUgVEJL^flqubXLQMf)^qb=B6MoCno|R}!_gg1Rf$N$aXUdtftEsF3(QSXgpPanyP=Q^>{r+m} zhaZFoYDU!|x}oW{csXDCV6RlIx)$z)k|X`nkx;j+3-k*=nU~!d%8s!8Am@~O5g_$r zYVyy@yJD9iOW0^zWXuS(&%8GyZ4zHRl50evfR<;U?uqhBTDxk*Gup{b(?X9op+MKV3$V5;O>b9g>BO(cmt<1_Jhn*li6J(0BQQnfYB#} z7Fftc{hp{OpxTW5?r2au`2w6>?VaJ))?fj&A}7=OCNKVqE&>Nz?I(s<(DSiaw;PQ= ze4#7f1Be?zj#~A!?}WDH|td@gS%VUie6>YhJ}n( zaaK^9?`m{_3gZMnGSvY3Oo7v{znIDKWl;0!@e-heI^8=F+BLi~+X_Vu3Mh%(y>s1= zIA;=1G_oRP?u9+|P+ZNh5$5(KiCB|ae1B6p2OTxrbS~ruI>Ph{Sc{x&;PUezd~sIm zXU@WQ3Ff)|uQUg5mT!_;t?$(od#taB#F^)sp{U!)U1JWzG3fNh^MooNdp&H}Gh&(1 zGHR|p%FiSpB!Mj``vU;A}!VP-?Wx;KbJ5Xw>!lcL|TZMJJMbF?gVzHq~_x=^!>5I zfLiIZvk_?9=Pv<(+V{K2Kd>O_tJSrTbr6A~?V0ynVt^$kkGxW#Ze$aA^Ot+eUJO|F zgunlHW7XPgSzE)%!QP$DaP14HoGKfzh{IwKyH1pM9uymJBl;YL??>RFWVI;eBrUap;!{zIGpZir_Hv65_;Cf6guOv6$E{ z#CgIAnN&u-f3QO@`hiJI1Tc}mi2%IzW_KcX@k4zrB8OI4q~Wb2I~iFdxX!8m(ClLQ zi&UfQZ_UB&<$;lg;)Jf*eXOTCqCdeVX+r%Ft0gDm^;dxI-*(OocWq*Y)lBa|78mfcck8I)f zoD17z(3Y}ukL}xDnEO1xxU_WkoxJJCP>u%Kv~CptVA5F&tPwG^!llb!86S@`{65?R zVX)c(O)u{7fAOA{PEsi_<=y6Uo+Q6X)e2QmNv>{)>Q^B2{_F&t&+F-3(K{VIQ)X}P zI^;ZINP5T3>8fe{iC6XmEJ?DB!Xb9TmjWrAv(J@z z!1!%Fk?SGpn1<<;S^YT|#?-_&kKU#$60}X18@jh2rN_)SgU-|AF_vn+ECaH9N;pDh zgZ<}F_q^2>uOIRn-Nh1t;%DzgQAcmyAuw=zL0EGKpgE8$Y@gw*A$p^ejAIRs1zw;y2kpLzaNM zHfgw@;WurvcgW$^u#kpcxVbspJ%FSS`xvQqg=L_=R zEIzrVdS zhdVrS_J48t1fF$s%!;fEJ{aP6bl=PI2Eb^xzx}x@Ih|_y+du@Fgq6_ss~gjeqb$EOh^^b-^fw+> zgFdi|#ed7b8TsWLA1lP|c0-v0KRo^tB~C8&5+eu^`R7NY7HZXMbL@o&;;n98l9=^SP}_Q`Co;DQ}6yA z3OrWWo_(`vk1wuKu*_qOlnOl0)Q{^gt2nsp_=WOABB9~?>$M9Cd*6xNiIW9+=P1WC*y+?mm{op#C}{$mA+%Ep zKso#kJK4zTc`Q{~G{q=yyP#s@geX(GXK<_rg+()d(Np`lGq{2@DD7RQ<>lx4gom8p#Z^eygU3 z9E!1UPsoWn@3=dH{tj=}S4r5XJd?R~Dg=+S$T23SFeel8G;eY~TEB3a@d*%(3hCb49=FDBqyuK*R&`K7;?F zm9qW>V=Zp>FL>=_Y7A`s8KQW87R;=P@3=$&{n;n>(9fJizHmhZYwNf+0z*o@q7x4R>K1iQSwh$l^Z z4g(Bm?)%KzpGxf>!TaL(IRZfP1DcRRNko0p24-5fpml~@H?w(^n*9IOOCtja*`)?t zyEYThJ60#a`o<}TtC5V2*B|=ZCU1Oe3!2>3I{2H1KZ^l5SPGu_v%*hSQf*Mv z2MkP1S!z`C4j)4-(~>M*v)+*BiILYa;mG$ef?gOfc~A%8RlS+$<6`Z1eRBGt0$ccq*!H}@v3Yq7 z`#=X%2F)CUHDu8-5dtZ9FH$t1^?h(IanA~>lC8*o7wC@rQ1M6jxlADoM~lPjgv?aj zgDdW|w|J=agjmtr8Eg!?0!=)cNpf}#{!fE%qLV%L!G|k7fDBqql=yI#@(*pCvE%Pb zM?e}B%zcj7lFVUYifXjAKSIklTX>oudqwDxm23qvuF3ZD3ijsIN}3@>IgH_kyI zBU&$Zx$~_1_CUy?GJs2*tl<|WIA$jlwTJ>Gw^uT9LREG5af2GpzBI*?@;z)2QoWzzn{2*ctTvOKIFI)mDhstX61gQ?Z^YD#G}GAT8pE)Z3u zKuj>68d0cP5Z;&oJfp^Q<51*L(tAa&93K6yXL*`yq3M?cC)KRFL>D3EM!ba60@j7L zDIrd2-F+FLW1EoFsmJ}Pg{G+?wuPcx5D)V5cy(n(+8XOBinc_&p)5lUmD>^_%0CF~ zo(0*RM(w}W{85`KhWvJLY5PK#pb-~s`ftx*;lfNs-!ggwUwwM{W6Z?-5|nJ1-MqX9 z)0(>~@w9=oBt~~V>{^w@#?4EWJO~AVm)*69P-owf3}%m8a5wLyjOZ(Do~$CprlkdG zX)V9iXFXY-2b6Y3QzLKtClXLM`v{cDmjo4po77f&Jrlxd+smbK+M#$?Q0^-A(8$f# zRu$p1Z`07t{Mqtyr<*Lc>+GA>Sbis%n#`BbD zjGPuE&cUqji|xb>K>Hxi$IIA7!ncEyWOXoAix{^mqhepaexUns^OMWZBmWsLk)|k0 zV6_`_c0+PvH=+`Bbd{2h>iwxc)>#D;J!H|EgF0Ss_be_N&Y@+cB($sO3$>p-8K&ri z+nr9m3hrlq-^`8r?9$Kdf-n9)A+NDvr{n}9^O+ka9byD{EJdX{vz0QPIiAvqgH0OA zu|67FM!aET=6|W&ykBaa?XtlA78Zm9cb-1EyD38;7DbJmKPo(BfN=f>Tl3yk64GjK z+$@bgXbindH#v<{(57j%-9TA}lx9D-?>Jb z`wRv-^Xl%;SK}D^#5t@^N1UOMnRm=p-eQ%51MlO?sNJ{R9DZMXPdWQ{?9o1J|BKQ^ zcb%m}TmlAMcrd^mGGHt3Qp+dy4gA+{N+`r zQgD|5*EsKFK{D;jYLChGm*N=IO|~~7XZFWe4vglOJV}L<%;52^V~JR&ea%~Gu*RR= z4QIte`XhGCnmeJ_(IL+I-n3YD3uYBx%V>MGm6M$2_s+q`BQb8U<2j|ClI+$f_4s}( zYGc1!M0N^C-4qmG97@<}A5sVV1m^uFrypt4MSb(5NuqRAY}W0R55fd5b~?O(Qbq7}L1tG>ZR2GBNJY8x%lP_C&?&4Jd%s%5>W(8lT4f z;~Mfz4pU$wYA3^0nViA)QB#xEMPZdPWQ~%N4Wok@c#7zGAQy!c)RlGYqVrVpvM&;w zV^V3*;dyLs{iY><2_@l651fP2J-?~!J(t@hekBNE6SbaO?F?#@LQkmTfXlXDP^EQO zk+UXL9AE9t^o_=?w5y6l;-zJ9wA8O9-YuPvp(3kGolp1~XP7gtWRECR@{<>QcnLC+ zTRUF79n073JIix&x><4%(YRSFG@r_LK^0y-$-sM!e>XOr?UV;Lm+#$apW-y2BpT}h zJam08`{sJ#R0aEyf z?fhHDa-Bgg4V`Q|i;YD(IA}%TF4E)$85(&#Z~b;^ND2%5g=Vfgo08~);#buf5 zdX?|p)|)u~-BOiFKO|0;@d&>;*ZMshj7I&X)n(yb5MDb5eH6Di;3`~*+wO^B#uQM* zGsgA5ySV8MMDk_#1Pxlh{)PQzY%?m$JI*4k5C64xg;y9bfp-VpXrWRyNG&*MTyf}f zBf}PuH7l+8-9G=CIiF<+-z)NFC_Of(!qa106K!T`BI8kcE|KB7`oLb7pWQU2ofip; zwH75hd<8KN(1|Yjn>JXFb>5ccZO$B{u)^ydDp@&L!O zDKlS*{?(qs_}%Y~E0S7&!g%92({~c@kSo?pq#NUP?+NF*Yel3SLhJ8~Y6n^@ceXoF znbpp3SZqI#O$XxHKxX;RhRM4hJ$M1m8zYVnP1$e40 zd;|*j`I`g*0NIy51W&JP)%`t9j#-#T@s#^`t4>JAJK{I?aDm1s?Hpy_*fl_d|lnu-+WmEA%-J4adA4&E@;@jc=e_C|*pppD003SDB1K@=VOtrN&HY6K%8TJU>nhYK~L4$>%uU_*X-v zoE*}L3QFqQ$#7qi1q69NKJoi`qOZFHpH1*g6 zTd-wPAF<}T4z?w>m2KykO$+PC@8K^S_Qk4JO{aT|KHmU^LLTq;9?zVSZU`}r6CvvT zUZ)dcTbVUKd-(dn!eMVXCFn_fRX{Ku;rF&&&O4`C4(Nv4{RlY5ge|LF} zO*xE8H1-CmxX`=!+58_Z0IZ|Qp2MDvYEby1Yqx#WV6MF}Y`64!&hs*-(qe`-xBcf) z7VH;laUGPb0nd!YCnzj_qwP>#g+xvv5g#xdpI^KP!Y>13k1iz_7#3d@cK@#yStH5p z@_>SQ?|vC%#4w$g<&3il9H=8JHT85Y!`5p;>PAYsDoMDp2S`{U(`(l=&wj6&^8fJ< zYe8!gDf7ClMIQ_|LtfMJuK*i zzLE0s6N@WMXU7-q?9JfuN?(-m0-bPEwOHhF4Xh5?45kmlqf08F4>Ie(3&E!h`CYP2 zs1h|;81fmiTLuv)E4SY8S)rtRWhk`ND$`ZENc5jjMYm7$ct)z=n>Z#QYu1!!Vb~aO zqmyF$_EpBn&1-Prx?_0kRs}=slIFSU3~-y-biTIFhuq7ysz?S4ao-wg{-+)Q)MI`$ zh8NHxBIQ#pKt1MtYb1*t+PQt(dd%$l>N{DYj7@o~^I>!pHKl>uYQ2R~lAyw?(GRq0Egv?bUZ1adE3%fR~@Mu9U2#il$B7m_+c_6-(Et<4Zy@^5@BFD7vZ7~L2(ZF|n448@c z=f6NBt^TK<4jD&Fx?EEah+GO(r7FM-*Vg>{30N~1jufMEJxVSlZXD!^C0{3bdmAha zKuf?~jDq&TqQ>9m{bF?;yudMF(oCntc?yEOg(S<_!51G4!W7}b6(uzRFl6z?lZHY@ z-w8(-+IeX!+ml)kztd)c`PB2p?n38_E<;-dV05PsHLbW*f!p;P z<9ij*qyqLCC+jEH69bm@PBJ=PL%QrnhE|~D4jZWTg-3ArrmPyIB?;_D0k3Ham8%@{ zuT5A)5ximNLv`Lz2hTNc{Bwiq<={MyO@wldbD3set00u}eXW-Tkvv+dZEz z6d8K#FZAQAJHI&*LK6pVf5V;^@pNC?h{m7qg?V#&AgjT7@2;ARklmrUanVX{!&ie@ z?pM(WNV(Xh)7iKF!{wa%g^7J2*q}xWC;!ZF-vRYI``5V#XbSBsI8H7aOv3LjD{$G3 z9;l#N5qiG3quZFNWy`kdUTv%X++0dWL#%~v4XG(QpmspFBQbp9dE74E&MzLewnd!| z;d1zfFRq4{_`~i~?*riV?*D%+1x024=rNx6s*nAx2No?S7PW7lA{h%}^1jq@52@T) z*l@96HGffhcjXu~(rtpvsJPNFy`YXHDWlyJPN~I2&VpQ=O$0UJu4H_eU`GbyX7Pvm3Jk z7~HNcVJ$SyV#}19_&kGKKQ#QS?AN)>Nw2#K7 zfY%29{xOo+b;52Zzv)TZeIoN(@skgZ7hSB#@dMkQP#>~?^?>QNFu{)2ohTWt4b?!P zSojoWXo+iC3)!FtYET^M%%rUM2Vo#Le;&ZN4Wbq^&u>bb?h1@9rAq;#@)Y};q7#6b z$sYgvWlk{{^zY6e9+&4PwP7}NEdf;z?^N?!-8B3w#{rIV(o#o}WwnrfFXVOr7Ls%G z7YS3ZNOKmhj0z^DJ&_PJ^!$F}5XgJkgd>jeQOgU(S&ntGy^$=iG*Wk* zpsqR>jOBCyU2bL^ZfrjL!cN(7(z7*aY>h zQ!Ihgj$Ze3YR=<_fBVqT4)3Zw;VYo5vz!aW#;oUvRlKn{tA&5ChGy0Nkz=WTrOaYjc5VQ4uM_u3_@K=+^9rvr&MP{3CW*oO)!_Nw zup8_G+4HyZeOjuuhiyI#KGtO^H>?|D)jb2dy7CMhdUCIi{KJOUM*Zpg8ooh-X7FoE ziLr&L+=n!n85?6ZNvLw)sKQ||ksEjfmt+liv-e|e8=Q`=L39WCHK92C1RNiJs z`9*|>$cDP!-Lvd^#lS%`$C%f0KqY2#YAQ0|m-5Sjb`ya=+$tf2{B&E*#QN8GuLeQU10UI2KLROHuf0=~sHFrw!^Q)B; zpQX0iNNPa0S}W`zICK1azjj^Fck*pO!4Q;iaat7zYke~^T*TJp*sx8d z^WdEA*2xLU@tpwFv?g@WyoiXLd8+ZCqGbJs$s-3t4<(xwe6NgGCEk8aM`U-stThRk z^Ohgk0uA8o%UJBExUM{+H;lw#UmHS0=N6#g{$~!!{~uRh9TjEwg^LM@w3L(}r8Edg zhlF&Oq>@sjghLIZG)gE)cXxNqfFNB;3_WzikV6b{U-0|gyViBN*8HdEocEj^&-3iP z@4Zb2t3!ZppFl7z_)m|j^5m{^1=BMfO)JG59Y)tnUF9L;gHNpPQ{k;uCJr9s&Sl#` z1RWmaNU=(L-nM$AlT-74S+QY2ayMOtv71OHoZ`gCh}5Wcs~Z_s1&Y55_t>ax_M~I5 zToz;Z`Wxz^zbP0` z67c{iRK!tEkYa@hJ0;E!{jeiwNP1+gdNpn^qN_^uvLulgQg|9c4{QLd3j3xCA&um@ z=5u1joUaf3=F=Ds_J)w3uwijas*A`fE53< zylz~0K(LUIRPSC4+cZy$xd_k<&VfsR8fR)_K8l3Y;A&zH#BQ;O5AB$PSmPHl~T`c>A9O(oP z{@f(6RJTZPRQxXHzR#Eed9vc$*EFOmT_KqQ5w8LwWG}qUVR8C+ty49(v(0%_ugqPX zi)Jf%-AtxkUA}fPJ3C_fM-IVn4-*URq4PfBAtGkEn)%R+irf#E?J(53$j1ZRbXn7p zgKwge^Q;94o?Id0YWrV0%S_Nky}mKwKKhhr@|rHl(6yjaDEDxI$mafRF1^V-zWP-i zj^Er2WRzR1{1`^7cU0+1pf?_n(T>lo2bJA?q$+u*l1a2Z$fYNCCnU^f(|afHx!2XG=tmB& z{K^HZkoX>ClF|&Go}}Q(l%ByVw-n3es0pF%c#*fEmtg9Pa!4mR)55nfo6A*&X*>W1~i0J0jKw+`({pA!`T(^?8Q!tdT~c!i1?T+)EcUZS#*f|s0#nad&Kr{+zbHMxWeN)=$+&mYUB2KF^Sh_^jqFD` znZ#H*LJ-uAE;jtc&4_m*Tx*(8eCg@sa;ZETP2ZQTzrQq!C#qY@IgPE5{?^F%M_f81 zzTl1VR}UyuX#6zqe{OXpy`^A@6a?H+{VCJK@xyA?fMnK}vK!wK7wA7*@x7;2&1}Qs zbt@x>&$$4-pOJlYX^#;7PW9G3&Q z+YnjzG2!?;qpJ3ibyQt?^VQd|Qj|p1>YPSE;s~d;%ENQ(wlAyBaQfht&<*968s(WJ z3e0Mu)`LT4DQAYmdp1{VwT}^&hKwY*x*iyP5MILu^dTwm-Kzo&beU%+_|Ir_A9uUZ zJb|-b?gSOicxVm$#s+&IwL%NLRpSIMEG(i+^Eo~kyI9hEN(459@M=VhLK<ZAnSYhSyrSytEiNDNuaMEG-x@WKhQ_oI~I1?^YesYr0*Y1JTugs$@NxW?@4j>B@VtNvyk(;ES);FNLv^Kde-{+@{&pP-$TlBL z1~9Yo>|A@PK5U+wj)9+}aCPFy3D1YBAycl!Nl}C2HS6^?t6sNkGI)lSR;H&v=CVD$ z(|LC|rnozWm=R~H?AEl#Z(k;}*}*YxLLtY)T3G(1leop)+u-9ZW*c}A%R5(&)Kz|7 z#f`Ih7?|K89jzkMyGSPLg6W2JlwyLcK!_(R=&9kCS{ibp=te~=-&YC zu*)$aoyHO}wv-5;iyDJUk;6KHrvUWz$oc1h;qbxnMmC=Ix^a$GJjnCXR5A5sA1H2? zy^z@XJ))esmsH`0jand3o6|uWWIC|l&q~H-I@IWRe3HXoN503joDf;RrsH@l&7yzA z{@w(`MFGp@7Ly*HtgUy5|Goho@VV;i9TacMl}u$nm9i^?9Ps$x)WP;)-?%6&4S1BD z`>a4^O|*;gRM4YYAD8;inAUczI;LtaD0M;It^<03S_!fV7pHGN6_jx6PJ2RDZhEZwAfyswTYCrgtx$J4{u!hz|k)zLUqbW_?7^mgcKoY|-l`L)ra#lL9^iIt#;(@vSS8~;Kd&z%~b zY95eO`=j|h$Jh*G{gD3UbfxDXO=5r9UvaXieA^%Sr=oNb%u#!FB{d>Y^>80jnV3(D z)#BKMahrEcr*1{uQyBtp>-I)NNYJ^^8vI1{H!sdO>v$e69}0bfd(BR2giU%d_Fg_K z{CR=>aWtexZo6?;Z`m*w5P9(V&FY6Qjd3q(M#E{0{?wF+f`?}j?_fszw^eX{Fk6qE zyVrD9`U;#o!&|p;{QqQ zd~O|Fvxtf4pl=87>)0nZbP*eY~Z=n;PR zpYfWgcKnQ^b9>vF7RchpY&#v?z>`4KwMtr!>4)s6Ub? zQW@b)P>wqjy9vGp(*xm6ERnMWGI7UB^G$sYoMNegWPJp zy@BsESiUD^usfh^=841JVTq%NW~`2KSbl&C%|=hJ8Af)9_FID-@4JcdPr;5)zbvRr znC$Q&3qGkeu_zU}79rl*{ZSy1X_Jtlxo1)imi$0A{FOV;6e>}93Im~6l_ftmrRx;G zLdL&pE4wM3j3fMqubT@-K2p?Vl8a&vP^yr9ExuTrmLzP@vKTlLofb&$_O=>!koS>E zuG2Ydm)y+*zog#B!F;d%#E!DY$wUEnyAOHKkbArXPV1)2{9vzy)}4`I6%Q473$JP6 zf|X&hgudnlRHz(6L8R1_auO~II>u~3L$WdO6b%9bH5j#P@?qqBXf5ce#Ga^1-un9R zSWp9yu5Ws#jnTAX`uij==QbfLvKaD(c4$6MyqEiH|DrxsWgdl&X79cGfS4Mv-Y*%5 zWMpvs{fTaeC>rB1yrl32qIF-q0d?>5y2VLOdI3igb2Frz+?m@=T5A)Wfk~h~tYPjXK`?tRG+iLi#8|!xOE>m6@wrKz7;h+x? zFZSX^4fJCvyCaNZkOR9WX^QT93t#XC#4D2l{bgsN!@`3jN$<0~E@9{CnPSNN z3;az={kjve6sc;c|K?C@q6HV(?$|QQS0~Wh$3_3?wrIRmYW*t%#^_$J;;kyH_!qN5 z#S`D*g$~HTg_%z~p8~8RXuE&nA@eu^`|uZbg(iSXCvO`K7}%UoTzmmVq-Oklh|IgU z-DrFwuUiH3*8rbqt>%=zZ{)7aNZ6AHTEUzI@H8ko9oYLcP$p>+WbRo7h}J0hFS+Q*80Xq`O2eOP#Y2aw3Syni3%5a zZ+ZV1uXJnD&D?e)dIonfN7V*2apvi#@J{eA%x$9abcDj)x=iMTYpw1NDh{pkw{$u`&TaQH_-f~pq0J{KXV8_W7E*~Jcrf@kG;l>9nCecgGu9+cmZi!j;&F# z<>iji_v!kpwK21;ruo18v0OWuKVK9=GM}9qMdnXdf1Enu)T3? z&rcqXyH59-<^Bl$)qLl)=PE43ij~eg$$;2KP_Pz4@sne2_-5a^H6S6zObV~P zHALRoL45ylv@KGAQ%sRwtf#OM7t=ogtLNYFW}Vcz(+1Q!!ve z#arCoW)w{I2c}th06qd>ko|2VMcj+HGKzMeqd&@bI6Ze7_t1{vp&tG2lt9@0O$vNJ z4QO^FE&ikkxHED5=9Eo(E}5T&e6Xxrh)Do~CFZHwc#?x95=?AjlQ7_#OvPI>ADT2f z`^f3FF#usQ2oHK-tSq+Xc|FQFA}Beymm( z@^OvpyJCJ*#D_gG_pZ9lph_VqDPtjDG#kT9m!1qE*!hn<>Ys8d?0Z8oUG>22L5JCO zBbldalj3psfp{OYp?teq_64_A2}HksD=~<}%kpj@0M9c^yWTgymYrKa#n}W7-&*!9 zrSY(C@Z4zJ9iRb63LNjhc>h(0Gbu>3HZR7h?D(o9ju;OQ<>`~iXu<#mQ4WEc&|fpu zHI+QQA{%nHAm9^B1PIi#CX`Tdn$=Uu)vFl{o{3Cli58OyzL##FJbp60(a(xsjZ(Zz zf4Zc*ld+Q~O77a}>*S};dnra%z6HC&u&J&R%lMj?%snwZyEQS5Hi4YTlYvv7_}|u@ za+_y{kzcfL4chF~qc$&2UFS`{Wc?y+%Iu6IgMg(l_su`~u06!EB?ox96yF-a0HyH- z&*>&y?}#0$E@0Tgn&>{l+4I+JFJ67$Jv7_-{qSBU7lLsmy3thg?e1b8u@63Xm zv+s24IQS0V`cNfmw_@%t8YR&HllF9lxYvr8|^ zj|5aJT6L}N?ATcn)|Grql;jDtE&OVy7_ck;ezeERUfo2X`m?rDWJuiNy{i+y+=~9f zX4Ut|ge7F%zF~@l&VFU`%QNwAqw(qklmwNloS4#bX!RC7T@8toZgxTv9f~=j7?T}ju)bf+wW6Xf{{!?{Z?HO_0`R1KLfpQ?vqVca0 zxoXk8K=M*y>w}TNc#2*{T9LEjju+)@mpx`a)5`KyAUlznV~HM*>VuQ8T9A&#Ww-4{Z_~EPJ4M>*55pPBP%;>vvRxN`S$9QndB8zR2b}Ia5_+y+P%!0zq zj=}P`|LnLG{El40Ul?A2*{R7}WUXMAu#joy$Ds_nNDJbN5sFPo~X;*~GZP((%5E9CzI2*S2POWqS5c{@tN6H<@m zF3i+e$D&OlbjCe=MeGiZgfE{k-SzsoeN;NLqLA#Scw?bAR)JMc3fZ~u0IfiCJ$Ns+ zRb%hk9fESRC&6oUbS3$agYVgdbYTqz-IR+k(i2lGw`+vP_Q6e8Vg;$U*Z`ZsyUUo| zc$L+9Zt}h7yI1NTZ7RN)M0X%_sCG@BIZRvxL}5e=#_m4&*`*T21IQ?!t!XkhMF8sp z8DHwuS~L3J<}P!|Mwo)f*1%tovC|%_+dgI)`&R9r&aRG#-*ZF|Ej%EIjC)2**pyR+ z%ISkpoD|exekpr%o=qT9)j=@8^jdo^Ko}6b$H@}zRpL`>VytDVo#5MD)4*b-T&|xP zmImS47S;%MczaTT3wBFlBs@2PV=hWn>^c={I02nbLUo8?h&-*~)5y>qPF(AqH*oV)yBK$LMxsDCQ=vOl6Ky{KMlpyUT*8(7-! zv@vHljB%n+sgG5c+Oa}0QU@>>>y~zE(^jMwH6*xnI!o)W!GeA;M@d^ z2`}7Hn5{I4>wNBy@(riHKw~0BFL$D()!PkY_v_pyW9iG7YtMGW9@`_yehC%4TcZ{j zf?*C8kqyUJ95pThyJ2Zt^w%z!Y7ZAPR5L>24J zYaHUUzVasT4g&8$+(^r38$cUgJQ;3KPd`L9BS+#qRlSuS632hR>O5U_C*T+ zUNQQXJXw#gh{9aH47<1B!?>)Z!)e@e4T)2sTL{AHXNV!5Em}HY?AFg=B}smnF1@2w zA>+f@nNrG2S+oCPS-)x7>vs1A9OOza*YcfaB{LJL{b+ zB!1NfVF+G%?|!a%m(AJx=SykG)RDIvU zoXgY8(qZ`>`BA9LBgf5#1EFmQ=SpKhgt<&`LD}=Kg7RJcdLl&C1SzDJo%^SuoCQE_ zGw{cj=!OkEh?NRp@jnV*)wD#vkXqhQ3uua?k($LZCz={(-LN4mSutl9*Jjy^4eo z#9-ZjRXWnu1JYUi7`r^s@PYt$;i--E)+nI~x!#_Z_Wy7j{s3^HF*?}?9ISA(_quH0 z+Th1M`$xk~CcU2L(_9LHJkp0AH>i?Z2qIWNjB!YdfjwB*qT>P8Q;vHr3BXUT{4+L`hV0;lc8UL0f!Bo7u zB70BJyObpm!_boVo;|7ydY){VU1$%0Az{~!PEF`d;;*SO7%y|qV4T$J-PC|%Srbzk z2B`h`-;nA7cayf?_FM&Gdll?l$d7FgWTJCTB^nOQYxaL@3&g?K_acm znDphkQzlj)&dt-7dWF+^d%jY_t#Ed#SFcn{*J?pkQ>S_bx^DKfD|b-_=veRHKEgl! zfY$f}mp`TVIt6wrgy6wRQZAJpM}DEFP0_9QjE0D+ojMZHXa4=xVM{2*tNl9OmrBDP zkpqiWiW{are90x1iIU&LJE>&4(644ILRM;)0RUtSk4R?f1c1+rWII|h`rkVV=F`;(ya7piFU%Wx%P@;|P^+}H>2lNr!s;f`YA?g5bt-gHB40OW{LdNEKQ$$5`Eu#be)_YV^7w*dlb=5k`eF6g`AY!JS|A9o zTs(GA94oI%w0G?nx_SQQIsbM3^DF3?!<0Nbo3p&d-Sq5hQP2K@vAjTN1G}(=gWZut z$@gEsBjxb!uKIp^H3hvKP)|{7KC=r+z1>Gvz||Rfc20wZD!O1*?M6n?5ay}^$KaY0 zy3z$1+{jlsw2p9gH94`lN@z-EmLW!5d&;Fp$wV;Gb#Kq_xtw}zLp}xoY*N)K^8O9+djbDW_Wp6z7 zf>yp1{c5R-6F*TR$@j)-6_?Nu-Gvh|8LUNxce(=-8K~yl8NB0UiGTzl8K(a!WIJQFQJx*gwy;qUGhbL8upX)7VOQhC; zZH^56{35{-Eiql{1PE59XH3GjZmF{ju{z-o&%Q1MdV*O1t(dHu(^}EA+S7ziC-GhM z><^H6U;m)ox-SrBknRfR3tM>wxF2k33VC8VL|f<`FXmKog_Nw&DIMDl+i<+``L&R% zt&D!X)*B+6uLke zk1>MT!Xn)&Feh5r`()K?)}$J>9cf0~$9z!v8LV;HL85#AEJ5H*{Qb*UuptA?*o%)B zAR6sayGjSY$ryO0X<71BCfnY-x79iuW{oGTBXrMCC*CN);d{-sQ`0b9=KYKmv*}Y@ zN#gP+uX`rHU9A&W_T~e5Uc)3%SQfTnkCgS<`>Ohdr1?$(qUUnB8H`?LH66^GhTZAx z$InsnDexS+nGS&dM42~o>TGv)|m$-f(wT&&f(5WG-; z@6X?|xnp|eBF!-C)#&pANVjft`M<{i6rgeJ^i?^bNo??Y+rY7u_8SR7?$wEf)89OA z=JMpFi->ofCEkX@0?VkbT@#I;aKIbWq{$f`M589(_d@Bvubon zAiPg0jgx(_6tScDqT(Q{%kq{wQ(M@obp}_d_U3+V3tM{rgY5MWbFA@$ff<*%Qp*7J zWajm^<_}09-9Au-`IlwGQ7H&7V+HnY?nz6yvIDj#0VC68`*3tA=oLnOcpsb&S;Sgb zcqzk{@*?GOclis=g)edhJrkG?l%=iSwP-8$e10uTF@JvZVz6HR7Zpf7Nd2}liUl~Y zH(=usxzM_S`#`RiUwl*zUb;F>wmh}p1YfJp?x^BGnXB_4-=d($p1)v)eF3_ELfotn zPnSQXq>)GXAds-Vc02|hMklaCXL9N;1=Dk6hHrbOx{`u2NG!Dgc5d?aN^xlo3^s1l zPQg(G%=6y0c}_~0gf%~H>)TT@9s8-xMf^{l6}I;KjU)HTwXrgHBE|X0^-IHD0&HR} z+V?ST_29{p@6#_t9XPKZpbXx&&*sE^yQ=s5{tJO~0qgTC( zNc9@sF>7xY6R@ERdRAIYbiLF=z2bAOd#cERe(>`qqAC2%3Dnb1n|;< zdeAeQ(yRT8@hX00CVR!~BIC_z@G>d)(;(QNZW|Lu{{|h8DLLJK$nl$r#;XH%<^9at z-2(t}%|9_SH!T$U&sp8Q(br(*O)TlupXwtl!^n{#-uy}pb!KI@x4{83-oLR#9p&oYXI`CC7i^Cn96o&2OZd#Wt@jU0&czk2mS!f}R(dZ&H8N z*2kd54m*r&j~-LZ2Ts~sIobnhgje6FW(?fPLZZ{&`C=#!t3A2aNyKmMumD&0YE3rkB_ zD45(3-#P#Kaevl}Z{TX;A=~lbzzL8FCVF@x0VBG}m164|u(;}r9(`w@h)_Em%Go026 z6~8hEM%n|=8=YF4*UAH)pI zN;ihA0Ud3EJz!>roUZ+Db@fxX=hAN+q=vD+#qG%+|EuslQoBtE_s8pdHZXDLxkF|S zyWFT^TO&kr<9oo-t(53Hui{2$5AAsvYzxu$cr z7|6lk&*%NScYxwV; zQHBgRASYyXd#2Z|SPbSs8&IuhE^d^eA*8{b6IGM(3(X%Vcwp zKtn8D{GgTnws6D#@hX2J3S$kGXmdbN(>(J~ovW*F)#+AO6%tI^OL~0Sw)upWM2Rn- ziPV`f|C3vB(>aO%P>BK^J>+OKFYy`)aUuQpIROPDtuc$vub=O3dudu@_iI$gWdNdw z-}1f$%s1#P;c4TOXoh2n%v?T84xv4Aw{KBv_Z!iH77ZWY?Noas9Cf@Vg=E+^Rqyry z;9PWRsH?}Z`q^^puZVlAN&LO|GJ2nOJ^ug5ALNbN9;-&(eokB~>T}pefq{gi3GE!W z`^`?mM%ybIeNHHO)vG{U`e&+*Iq0ySRn*4RaaLB21ALcHJDraDr>d3 zXRs5==R$E0XuxNV-|A{4B_aWIR1lB5fvwz`S@#iRudHr9FsUfL;49I)6jmJkGxFz| zPiNQiQib}#gmjqS+!;Pfz?qoWu$_Bo;Q%^<+qhrulE^$^@Fxm(X%}AgzL5reX>WJL ztNkZ1`m83!lk-R$mz^YQcS1&sN$$sCUCzdz#VFd`r~9>hLqXXq@Y}pMT+;HUTQX$s z)r4#$gTENJ>*h?8W;S{WWnKNBJ{-`xRbM#tJo3%l0Fp#>`X&i4K2SD#>ZQJg^w^aZ z+*V@~GOIZ$cYGE9Ub<-qEk0v;>lj1wiXPUEQr3Z`7!3Svq>Jo0z+uRZ^=u}N=u4Oc zYp=lY8gwRfN6nHU(gNF!gMAK{mk1JTa=4c)52O3u6>9gpRYd4+*_NcIIm2M+${X=b zaoU(Hp@C-S#X|1jg4lB4Jg~@K6w22A_wHS@NS_(Ad{r@5K6$9+COk{#cRu4|kkKP$ zRw)lxRsskdA+THxpfJC>;D0UjYyPlqT}L2gI(Wbi*{}t1So{>79YOn;+aC~nQT(d9 z5CI}`{PhUQk8jkqCS2Iyr*Nzq7&q3tl!iPIJgi6z9WRpTSi6B8001*ut(@v(@08KWo_t1CG-&~HR zKxw4D(bpnI+f6uTWOa#d*AN(u*m!E$0LUSBc7P^cb0&(1`BU}lxsyaNyvcYtwG;9( ztT1J9F^IJdDBHq&i>N~{!BThNd%seal|At&9AHpX%B$1JT1QWe^@R!ia)s5;vYqaO zuw!wr*V&5y$**aW&0FmYB&>H!rvz|~7_O8MQF`%g|Z@I3Sj`-2}$jGDo7 zjsmRMEg`atcbY3d9a2|2?A6jnwh5mqAGdc@5)VE#u|Z7pWKdh_)78pQ@7ET+_H_Pj z@d|Dca0KNHGtg(DAL|YweS*h`NGZ}ozNvOz1`q^G!+|Aa*LAX*-K-q>27b0gA-N&j z4Nk1^9GT3+gcfYVm@CkE%8a9ZdPlP)A2MObpqxc;iPz+d)ARy>hzBVoHC3EG`~x z;<1g}x~48!Ua;mj3e;m|ovX5KzoA4t{QB^;TN3s=A_2pa6YcJf$7^qkyvK~)vgw4( zTNogH78DbqXg%$|Q_HnmwU*3!aoV@y)mwu)AIbwIpXkwxf0qF~nyA$0`v>IzU3~4K z)~WY9w4a5055O)$sqPL7xO^#3&8Cw=r;q<%9T92BZjsBI>Z^3OQzJd5coG)b0ld_5 zIGLLVvq*W#Yl#d}m$TgGrQ)aKIkVcatluOIY>{LH4IFW^u*QX+pi;w&{nCi$W^R*x z7CD@3JC*21Ux^+FoUWSHbFPKUP2ae*{P#Q5y0xcwN8P8b>d~CbMAdQ;EQyL1m{E%R z&W9k405t7j=r=!wuMz_2f%$Iw#u_0Bt){X2SUiK$ko##D|LzpllNZ^3ZO*%6@m?79 z72Qj&*EdEdo3FYWm0tpLa+qfX@0*zV^wkp-+at9v#sN)Tq>|2p@yzieW%A2dzFoD4 z_OV7@6%XM#)aAY0!Zv#uml^70li&b9F1)Xk$0#4$ey=gEwRJUEHR#-oZk#uz`0X?7 z^!CjVRNcWRI5+XfY%B299B-tg=4m|LUROS!cP2|U%J-LHEL>#uHsF~mfP00xJmPk4 zThyAB#d+)R-vV1+3X;bIs$C)ovJJ=kum!Qw4AUg;@Ra@o=djoW>uG2huZ3n8td)`g z=&=X(+cp|*&qbMkJIiib47L25f;JImb#|Sbihsa7v*kr(k_7S=)uHWOys8~ zz|fXdK92jH-^h>PtD#3gRJv zlDoE8Qq=$LDKCqtRW~|s_)(GZ|Fi%~y3O}>K;ihG z2UpH#4u2!@YfzTAA!zVxF)s8&f*MBe%$h2Y4hvL&0){_<=`x(MgnGN~wub~qeH}|t zaYW|nT$WTrHoqE7K~BaZrhf(tjltlTig2L?BTlQLVfdlc+M7jt%hK9O{PiE(jr$8I zy;>9<|G`|X50Z|qpSIuAgOH(JCWWfsL(R1U+dSs;;X|yr1E6cxnaLFbT5)-Hyff7R zLB!-9JI>I#39mgAf*zDIsX`8dB7Q*t9oTUG16^IpJ0Mv73oMDmnDI(!KA|q>XN#!*m;&-&>e5E% ziAW#+?}J%TsJnhNz)6P$%z%6r2_)a^z*N+X*kPyNtzM1>fgX**?14d;J?P@R(j>1R zrLw7R@nr|c7kdiw$26q3?>9?}r-=fRqe6gEJ;Ee&A|sBH-Osmxy`YA>ti%cqUc64m`%{NA4h9W!+Rz+?m0EY#&|Pqiu)q5%KY z%M&B$x@?G{R*_XOy-_oi&XZhgpV;Ale%_WSvz;sAA5&XFSr{IIB50F*AU_tjsgone z;qH4G$0lxX`-V0RPK(#<92A{5uP*Ci*_?Xh^41T1j(>cAL|XK(BGRo)C-H&Ohy^j!HsXn*l14Y&98_9M6#y$^=(cPs#8+cVf;~lGm5TbE{z4WI2 zWv|dZhvW$F8zeD1Vaa8!3Rw+k(l)aTt;f^qJ_#LUi{lhxL6Z zo;}(Pxbtm5#AOD6%86to-L02eG;|>og0%4!(Iyrjf5VsCgvQGEjkKame4f^5Ozi(Q zn!bkTCPJ|3S@qVh{=TjbMEQ0lvknKrMn5unLpD_bSM=B8hr8#w3w<4f$D!>=K-%3k zK!t+!ra<^cZm7SWfyC$_BIVxb%b=7By8UOPWpRWjl45&BCz2BTjV#vxfFs;J_{)-q zKE|&GlyrKhp57w7O&93POYSfBkXS}GGLM7oGPZ|LmrPNYNjKU}b)E{<}|hw~3oVR{AcZ3=`qh)81M*NP-`H!sn}rVGXM=d_NP&Mev^m zbl7C%wdx|BU3>L50J}YeZCrfU1E+< zAMpBRPx}X2b~v+c2-x1!c%k5YShD%Ix_lqxx)KPe06VHeodE=TGlaWp>7=adm~s%Z zd%N^&JxbdLwhi8?W8Nq*WZ-*&k$8ne73fd254rybn298%u=#m+#?FPM2zj5h8iQv4P?m`3p`#!WO$%wLsgh9rjbyF9-V%XWrbN&IN!c z=HNM1UK@03-^#fOa)p$k4C?*h^q&~~(AGsbuyq!32B_^i(46e7bmfDM;Ar^yUEor& z(EoWo0Z-I1@X}i9gS~AAA8{Y+2ponzHag|Z|5~IagUc|kzZ?rkhop(>owEG%pcXrg zKNh(Fi`_%Hp~vj6VW){?k4^H+DDqrig?nKXGD{MO;#1ysTu*$u*0yo-KIOgm3qL0< zEHu)zG>ifMGCfy}lzo=%ucdU|4PI-NNlh^c;p&c8<5J)9>6iB$cR`ov-GunnDO*3?l``pX z>81k7p_*bqa|NV2p}(OYXMeNx;9sFu2QL3yd@s0#MP@l5ym~Is8`HGUSWZNq0Y1HpCln#ijhnq3 zUHN1h8TT#wd7IEzJETAWNT(cvk#|4@P+&VB7DVAoo|JA%VlX#)V63zmKn(nz^7xPU zh@6!<2M)3S6B!mCsr)oq4MmL}=I+EF@1O2S3R2N*&q8U;FmGlccLU+C|KB(F-gI}{Qu4UQt)cp(BTpCiuMf9UQ|I>J9 zgbo~?rKhd$13k|a8vQuUEYpQj-Z^g+Jx>Wnl_iw}nnY}4r(AW-0z%}PJBG_{npx5y z>kPeZKe>a)wO75iD&xJoR*aY-qW8bSdOjOs{}Kw;b4Dk&IEm41fpuVg3*7fsYybWJloOf(P$)_y? zZzOCfoHnG{t9t%zKH4G}*o*8tl{z}S+u^v8wwXsZDdr4d4(`1R%B=}s1-eK(+g`EX z!NWv_P%@uaxWo|9Rcfsl1ipC@ifFzZ`0))9WS33SFTf1FeB(BFi{bG_Z92ys4#VRk zkS$a8i6-sI{7BpD_!3n8#)d|ik0w4a!VyS*dRyCt^5{jqXwAlfkunxk-mX6%hRM%l z;M$o$2<#$-F%CA}B#xG{)Q`0~I^Vtl%I2-_KKzVGA3jd)15_OSv=;=zJx@KxbztY5 zZer9Q53`-h>`SQsSZo8CSTWrNS09Acm5g;@I-%jW^`C1g-?dzlg!j`>Ml@18qH$|; zv+**acY(J!SM+oYwg>evS6BjyScfAo0?Va5&%f2yg^;?CF%$QyiJiQB@6eLb6N2(} zc0s0>5IRjh+93AZz`5DU5H1jvp5i?QTAO{IQ27y3qCJjfN~#;%h>6tM7(qO=o?nS# zak}902eD`fMqF*}0Q~?4NdNydwA-=vbs^AzJec(F%#7p?xA*g-bs$^Sal<$rO-706 z`4sYUj%6!7pD#G?OGQ4eKk$#!TbSo;FM>l5KopygE_P!kF!RE{8MbUb0LZ|lNTpr@ zBK#o{9$-1of9$|4SHpdt`dgZOBRq^`={FRT)0T(73TJ8DGAJp4Amx}8Z65%k3V0a<=wHr*8GepJFm> z4`mR~n?u9#-Zn~iB|zEMM$nK&H+?RnzY#?H@rNivl&6=O$i^5B!_qM`c-EuwpI##g zcIED$0|0Og^E`p^TH8-8`MfD~;^BDPiBuiqVF*?6jU1YrnY~ijUjTa)C#RJ1^bL|a zi3Lb2Tplj|I;gV>ddr5>`cuB>b;_#9#hP-f-29L-qxR$oDa&T(ytoGoJ{JXPp&@|5I zjWWsWx@cH~EKVMnreY}fBP2&tj>TWcC}5>BT;j}`oa@r5XC>;o5~&28%o-Jmdl)ZT z_glU)2GI^Y2)#i_CQ$pVI{@j_J${vY_0TYX+7`ASU2N*5)!Bo(4dF;~Tj^|!ok7$#wyYl2WWB#F0l9smAcpxcmiDBj>|O&d*Y zK)MvAVpX+AUS|1s(VR`MTiI2MIdBikftq!oWk6q5aKQ_SGuq*HGMjrwNDZ3Ol-Lhq z^v!$ZZ<;$uPqG9kRMTfFMRoI`s zKOpX``X?T)SL~Z=CYjGkZj2qQ;;De}@l2GCQ?t`(V_S^ZK2zP+yzi%W<6bv*x-T2E_Xe4kbq8yOlR=|P%|W_V zqop^I`qW5}W#D1b+fJcWjaQ&uMCS9g^U4Xd=Ms4FN)mn&rJfGq=Sl$Tr>(XX|F5AoNQV#tR{A{!!{J0kk*Ey> zRdH)E+bn{7#yEGW3UN#avo$fOhuAw`HffJGV~FwCjmLsRHfSR#B^Ca8wD5a?I){e> zSh&dSlbkzPA_RW7O#{>hk+Ap51!6#$RNrF2OpzR9>X|$?^?$hf5_qV-_y1{EB4o*0 zwveq*_Q)C%*>}lKj3vf8RLCy-mLXZPFWI*#h3tgHVC-XQjCBTM%>NF3zOUcw-}-vh zz4x4Z?m5r$e%{aXv~|xxVwU(_j9onz1s%Nb2E3SiNh5g5zoOzH9P{Mvk50*lOP|?s zb)`p7W7+)BE1Q3+^UpP@SG2wwo1`;;L}lhV_afTevL1kxjvSWr=y(%2obOC*c9x$* zl~+u)g>v|>|K`P0&gUE~9Gy<_07RleB+(mpItRW@paeN&7F}R@+AdJ!3G``~##PTL z3O9ab?{p6}cdvt^DH;c`tnWGBMcbCWL<;sHZ#Y)xb-&U)+GA~u;()+H6eQQO74GhY zP+wtW=Sq&74XLp~ciZ5f;?<~aaDQi7H%7-kd`Sy4r1)94(D8yvVWmdj96ovff!ACF zP(A1Qa)ds@o26-mXpb5Av>-JoJqX9O;>+mdt(NX=&~hhC;-*FT_^Vl^w{K4jbFJnj z8bw~`NyO$&yZ_yY70t@myu&n@xA5d?wi<*+SVuL1b%LfA$ni8bQ`X^l^d9sG6q>Ky zTh59U^p4(iF%mVq|2`B)|M)Z4TT7XMjj%=t?5C(=Ee(G)iKF4|#};XcM7 zWBS~ETr-v*R=OCC2_g8+jmH`mxEIKXbI8G;CIcnQ1C4`&k3Tjkc!l;BiPb-q{_K_- z7+OhI-#pdygt*9g%|^5_3N!mymE$Q{`eV}|EKZf@7XzFi{P^sWuL4>A2Sx+R<5g>e zlEK0xJI~>ptzVZs@@P)JYQ`m!95~wt>1&A=e^PKxjBVUHY<03u)%Quz;H`!cWtsyR7*#PSV4C*l?ZYbj3^-3*Va0_XmaMtE5DX{8p;r4XuXli8RA>I zdE3M69Zs0=*9sn0;F;#&WAZzya0=`B`|v)x+b3E*@(1;r)BwlShR)?$S0m_zIfrKx z^%cJ2J}q-Bx)HYCP`)5lr|6-wZpy>sikx-Yeli$Hoyg^kF8+fkOtw<&@oDL2abXj# zEINISE2P(We($*iYtGZh9QyheDQ}V>47&1#@O7u%%3)c^Bv~|({qWN z)frC?7_{_l;n2(Yz|Bgj!~pYn2=N0kDp6&fslZm!*YJ)*f1MU$k7-Pa-@V@OvZ4B7)AWt0tXGN?>{-*VGQgS z54sg=?@^}U8Nf;A4L``d643BPUynaAlcBt& zmd8`0U?ox7I==>Elq9p&km;pT!qBF+hrsoA28kITL|G~%Y>_)7ZHm6V*cWm$(-QyN zQ!cwYFGm$IGnU?9x6o8R@S_IjUDMSyvP?X{bw!eE6;xdxwhlVxmflFr>)86+BzN#p zh|x9#OtN8uZ-?mffIw8w+cHfndoMhHQ10(eysR4W6`%T}fab>%&z|>4Z~Nr+I~q&V zTOKOxq0?VF{Py0L;>+akX?pDQ)_WvBlCa-MB!3oD%cCn89}Dwzwa@Ddl65iMfBtep zXHNjLs(Z7d6D?%A0%GUN92wl*<5BHulA_7=Dbhr};v8nO_~q_6*DaRoD2#+St!%`2 zBG-TE=0_mL!(y|lk`4&g%6kAC%i$@B$~0s~_3cc7hpb(O#Zw5&J9m~V`?irJh+Z94 z|FX$Iri*V8{7iSb)}i^1%(9&a%xXKsq{#BgV=KtZI!igOgL0xbOF6ElB_LENlJ08D zh+1E<#WB4j)U#{TIItDn@Bg*&tL;(ak-kRW2nS?>h|~{Qf@FB*8ykM0ysl6!NtpT? zn&KnC(pIxAUC68<^{McZh;`y}a;j*jnDUg*CgGTQT+*{1tlVb-S2>=&YLD<_$-HtR z9vs<$6eal`h4YtqJ>2CN4?ipz!#iD)#U5D zBJ-R+Az7Bn)gLpP4#3jPnuD=|S4~R%@+{(;OKA;Tfcm+&jmNpxBZsXJ-d8>$Gm7bV z>aEfYD;#Bap*@vT{1tTiYRL^=7&c5RlrZIqd_r3`#-Bg)a9|sZJh8p{uaelI!SB=S z3T}ledMccO)9l=MkS*|7SryLO$tTk2UqblJJ`6!cxQx8^WLiqM+Kg4a%&s3fJLUk) zjO)==HSa>7H5=??H3TN{gsH}0@OPVeoW3^_axci#S&rWVq1&Z2WmIQo?(qZf|A^ZJO{`|M;IA?VcmaQ4sgH7_pCzz(L{FOL;>bP$C$P zYo*oao|;hS_sj0yb~`KVb=$|~5oA^%mY0WT74tWon+KKkR)0)i&rsi7aM5)t$)J|; z;`?6R3520ut*F7`70TkIY;a*c5l{DAhUl;wOs3*$ll$8ip^}2 zMqG&t!W-ThQhlv=^D* zPa6#~Foh3Xe)pMcyXC87EDEw`cBf@K)7($Dd*#xAB9x)8!|Io$=l99-s^10{I$Y!K zvbcG09aDWXR|wUkF$y7Dm1=LB`J<0J3<+~?yfnXmS47{G&4Sf!hWJ{{^?R8K+1G76 zyyuR;Xqag4@2Ws{B1`u-MtZ{8^JQmjp03PCo$OE0iD3x3L;iMd=qm3LSNxqT!V$K0 zy`8s8UdT`{g?F5Ks8E*zl9X|Igp;rk>oo!XZ7on@N|%lQ{RQaolUqIU>4S@l3x0)Ap9$Lx z3M9VOvn$8kH2qJE_66=Z2BO8=v^@RlV(e0)&YxJ(4?`kIQcJzb<8DB#cwe-l{5=(x zRK@~l1ly7%EF0P4@||odP?(uCRlSWVNc@5u&kX8;3s*s=>Ff{4Gj3FyeJ)wcss&cV ze4uq`kPU_oh+o($RPj##o-;lJ8iltk4MzTh8$gU%gA0BkUYq8jMGWI|W zi}c~hbQc7@(+rSIa<=bkP#r2qu^O$%$auk|_HlE08*fwXgZnBbgwQW12 zt@i;BW3b3khxEdTM`?`DP!IX?_I}@NaZ;N8TZccL9cT`L6vq;;FA3wB4)@CgFJ%YO zT>xqjtewx07`<#fuV3qn%RYZEXuNy*x`l4u!1?-0%kG2A#}J507yW44bll6E<+H6q zl^FPNtT99n$$t(;>%;x8QNKpwU`wx2(Hyr zAn}8h8ENI%Zbg0PaNa(%`Fsg{XLD}gw-h76AjJ#lZn;i5o}JQmTvPWRNc7wwAi`>& z*MT3=HeEMotXSt!w0xCLLF{an|K zzb_`j@ZLorY5LI{0m?C57EtQ^mh(;G(U3F-hWs*)$-PKjdSuS?2}y3|pDWE`J*5-| z1z-mSZ`k=Nk_gu{rMjR`<YMM{;0(U(nIw(vX*)OA7GO8{2)+P_&^_`;EO$Z| z8vKcWhhe5;{ zyTmtEq@5s7@?Yk}HKGjP^Upta+8$bMYsFJS-HO4pN7~1KR!6^T`#Fe7OT}yD$2=M6 zVSLclSibHHl%IDGS?Fo`?tg~=U|nXq(!-a3bgaCt)9*o9vY6N3zF?#>>8&2)@`}Jr zwMXTkw6Qby>uN;5X?O^MuS8%quYsML{j48mb(=TLLL4jqsn|CXnCRx8X!qhmQPl01PHUNM@B@nx&X+?p=z!Zb{dEEE3QdD@$|+X62~+63$x zz;zj(sr$lofZr(JJDV-H$liR7%Dokt%N4=NX&o21vGB#*)!6%XX7BRu?vKTx6(6ZA zO?!{kCkBCz71pbT&i0nWh`R`ArGYPE`i;eD` z1;o(t$-zJ^_s{F~r9@hG7bvPg{?tIDxY34sw?LL3VsT{IsjI*@=gYAfQm90`A0MIG zM$pgyQdV~V;R6_xTCDjE^@5R!mrn<6JX5?@d({tyF3oJ-06MsWf|uV?QgY`eh*byT z6JP_GsowPDK8dCBa~CpJTrVr(Z|~LY`M-v=e0V0JoIORYky3d{hKbwqs3o&5jGwtp zEwBzeJlW>Kz>I;51YB2;>3sFwdrBu&KIL;Xw$<`b@t4ZaATI8@Q`@&-YUN}lwf3PR8&L{I{hNjVNb11@|XuFsNS#Fcb@p;{O zeK==otNN)5oWZB7cj9W{imvpF{Hb%g`V(DjYTLQ15Vx7)$+qsKBf;MJk^}f=IC1Sp z{3Na`DChm|%S$ZdBNWZmpR!y!9{$G#K<>C8@%eIGAxVGQP9$}P@uqix@Pn>BgI9_3 zX%SY#bNBqVOE(AHrP3BfBWXx0pz$qa}Jf3 z_WGW?l@eFC0I0W=}VaHWHrY#DR|jPRE0$y)W5HPQRe&N{n|J1 ztD1J=U@xQ$*Wrn+jESx$*~-bKtsztq#!+9kO{w<1To^9)d)1+Ks$@E|#VbrW_9CZj zMHiHl44_N(yb(Lb2bvcTq%X5`Y5lZNaOG6W+~YbG{fhh8E*kfNy9Wz_?Zdi9-CuUw zc~i#-LV#=!2I^f~shB+|$&yivV;v*}43d9u+?`n;Jqj=w4bX^se@b3QB}_!-r2vA^Fhuu7X7ROw+`)ztJWb`G8fTP_Rq^ ze>=*fS>w^HsiIlj zklcC>U}92(y{oa_+j)?c#~*{Qncqmc@NkA{9GF!N*w4olF~+>J`sDl7;7hK%h@Q1K z0^7O7^5{F;Mi&s`fKSYI#~=5y7Ae8C3W&?P)vuB+#&;3*is3HH>G!MEk5(Q1n~`*0 z(`~{Vl6-e+Fv%g&RNLkpe|sR(mqbcZ%Am&`s?dYWr_S0tR#I_PZH7ul7NPlrKM5%# zreK~B0n1;$3tpLiDHdjDVd0?{XJ#O4mO&3(E;pn+p4uL(pi-DQjv^Mv#abLBdGQ;W zo?0BW(}(6U$;fWnKo5I?^li*)@$$olK&3IEj)^Nm&E;!7n4_e#x#vuHcpyE@zy-8-m1-T_aY_3J-erClSR7x zs`xw{B^*|0Kls5J*EjT+U)z=N9f@f>ZgMk3RX1;}1T6@fUh6YxdL9GykJZJOA@z+8 zSatrhm*?O+i)Pt%R-c$-bRI1-T3alI!AhF8n;S3y4g7=u$81>{oUcOo_a9%Dn_@IN zrKNn-i;A@1ttgJcz*L>q%j_3j%tg0zk$k}XPU``Hv@I*iFOac^RDIY08WaUAMq`hw z;sfy(4e!(4FGQjn;t?1T0gYMT`%K1uD0HIz;hWm)GU29h=jMaF%I(Y^nzEq^D>gLZ zV{73UqMW?r`*t3zAGB^lvwfhi5qp4J?DbDsHG21UM^7>7CscSxEGYvDMUd&D5;ZUHkrdo)!p zua3rJKnv~GV6X)(k&dDxS*5r-#=PR5=&7 z1|<`*l`j5J``#)EXD#(wJfpF0CWA5+>h+Y;>kC}45ATi4=Tq=yPhvv z>ey14qfm9kHSS9%E%vG*<70&rl z^FMQfrWJ%Tq+N)n3zB8dM)_)QE|X58YAhPxI~xN{oLxLM5+`T267+zRxz%u7v;R)I z-8%7@e_X$x@UzZJ-;+R@QZ$J+)wN;w>tFwr-7QHCUr;(R6pO{gyV4`;G2^9$N3m2V zrdK8CK6pew^yCZKGCA$9kh&2z#E!P57RSnKcg+Y_+r4_Y`R9UH=>5e zV9E>ibwloESug|MRehQc1TZ)CW9J#I6Q(YF^9G|@x+p~)m9Dpe*%>?&8m3bo^jd6NQ#pOKa!+_zj=*C&N6D z&ICThVDL}^#_!}enG0f6Vm1w*OXM`9)ypwdZ@c!WSeFPGk%U&0nCx$ zNX0bQ6UgWUGh)7aJ!GIkp3L3%TZ)Zof}nA~vz2MW$kuswrIFE%x_t^+7r}8e&*MgW zyk)>JhHu@3a4mSed!@d*)c+p9RIaoz62uulYoah_UNAex^C$oH|&9&r*D9IQoHk{`S58Jy`gQU4P0bWZNg| z?M%43{+O&gw32Bqb5s{ksn%Ke0A0}O;3RP$kuEYR4ya9eR)yrS*^}K5Y9M2kS>cK= zz9Xzb?vnNTE4kS4w9X(JKoDO#+PffdX8}eQm+o#aV2J$1$Ey6a&i1-pI(H`&T+NI2 zugBG3T7@pNM?MdHwK~|a`<2ChdF;f{=n1|SwH@RqZ>e3qyo+; z6nPmi3WGb}a#3h|p5`Dvr22NQSPiJsiH@iUt&puU-!L;doE#YX!-ayRe9pkfyL>NR zK4dj}zg7!wj}oiiN?x?zqe+NS_k<{C>pQw`SjVofhRGsUNUI_s!sQrz`!t=Z`exk13zey!!H>&=fyQ}?t)6eWw`0B3%(%ivr{=19|69XPfl&!x-KX-pZ+KeVZ(D$Tpie(Yd5(XU})>FhWR_>A3IRW`*( zgO$aN*JNc1LO--?3Hm>xZec<<=eSk$*0Wm$gVPR27k6 zwCC(ppU!{+PEI074|K&RSi47=^TjzJb@O2PdJR}1Jh$4Ckjos~(7mfyu*Okb@3#PJf|J_j34K~fyRTMkdG`6H17

)3P@VOLcHvHp#*%fb-6t~teEIskngNl=;V0>BvQ{UafpM-F;Hn-xLV9l_i_Md#!_70*|g#$>_o9m>AD zP%g{(VHRF3Dpi+*MMG`yZuuKtYQ^Jil21`7Ek!~3Nop?fCb+J)3bzqr)VnovMTdSs zn}<-fIj=@(1qGtfriHTM0%S2fzh(}{7~mQQy~mGz1IF7v$pw~RvT!2ZTN6M{VTb3y ztA2jHBd!8RRldi%G?)dM9DU0A7z*%ZJuyJ-@-bihazn~{1%=)qvCZD51R#WObfmlyZ%aQH)BBPa0KjVfDc)6)oGn?k9Ev_dV=} zel7~gf>D}?yR0KGVi%0+tS8DB8tN7%V%Hm5zUn@kResHI{t?XuHCh`r9h+y#+?w%< zPM^zc*BRLOpeyCEDf4mZsC%7HTfArg40>~tE9lGak1k|6r>IJsd}4J8dB6M9$j>I= z<@5daOQ?d@0oGm*L4>E-xVIT#K$6<8B; zB<4);a>!W=Fq1@3FUhi56Ac>o20)2qWLWP~_}%omGsApe3V1o5 zomV$!6$}+fW7OQ+kn>yN0q<@=A+b9k+cb=G8+Wd@dOEU-BbG{q%-XYclANz)(@757 zgTTn(u2`&rw44014|N7>oeeF&*!}%G`_;ct?vJO#?rx(E7>@MXBfHSDQeHv)mV(A% z+p2JzIHK+BR{)XGAO3UD>?)o;3cOTwGLmxPirkH?Jj)4cO`+02xyrIw2+Zxn=A83w zYpLcwR?nNlaZb0heb(tG)U^CrxS0T3u_t&lF#TQIp0t60fLbo^72wxdVE`T2k~V3@dm-YU-5+ z&_Os9aY2>fkHwqo(Iq_oL>;d)fO-CP^Q`8p`cZp)!P%Rplg1_3*7vnM-L!;vsbe;g zvH8Mn2j=&(iW_g@L+=T`+FHLkV&ZVV!VM35@v8>FmA0hs%VoG*;WsQtquBvzpnV(H z8DnZFU|Q|&^zjE`qCm_%_U~L$2 zDVn8U`Z1%nk+$zqz>=IDB!6>@EP$O;J%Lw4hr?b4Vfq~{IEOV*(~TM*;J>i!(eyB_ zFQ!oF!;0(>q6;3_hq(f}xgA#Gv7T~}m7&ike4Jcdw_V_@)?I;}Uzd@Vs&Ym)`(&L9 z7DizO1*pkGYX-1&6!8x8leI2oTSb@NAiiZv1{Go&IIF$)Wb4h5UrGx%A9K(bkAWme zw*7;d6^JaDKkD}yJA2F7!Ckei2$y=0@weU%fZWb$I}J^kf(WbDA(%Xd_TcV^Z#5*{c#Hg#?( zyWjSZ2|ZJ7wRTs6ujFlzd)Li~$lbLd{3A@;chBt7u!4SMX^Hbuim8d|i<=WoRe{rM z?QhY7@t&6`+jelohW&%)XR7dZZ?sy!A=I!qkb(O|{pn0X#MGe& zm(^DfGDk!wKG3x>(vy0o5(D@jM z5KEw3;)gW9=N(>#MBruUm-yK8Ajr9!0fAkMn>UPKX=cJNv+E^H2?NP;wY{lk^^Kd4 z%K!ss*MDowZ3!|_-~0ja#7JPiNjsaHzU{T37fIWI;`05hga2)gxxX%*Keg`h0EokR zDNOn#dluhB+&+^Ne=A@-fw0wFC>KITM$vP;Ql=CiDVNWn()jIJ%fOlE@w_y9e@3k~ z?>jlvQIf=GAnMOtcFROFe!3PiEXlpOw~-=x2abF7sx*(J==`}cfSwy^U$*0K<)AFd8Q8x+hf+^hhQBYsywIQo$F66kX-e`;>u?{C)Z z+>-HpM@I4akB#*8Cp=q4FU7WD>uV(0 zu6EY)crczF51{!EEEIq;Ex2Qs?>E0dhqaezF+o+JESj&qi5uxHMCx%JJ zqeU?KV1Q)|jh#7nY4 zAR;?xm%&b&&m*t3b!imw4npWBokEW~>#^T`2h9>7=%CB^xD7FaOGooMlW@?^iEuQ< z-=?VAM%XeCbHA=|5wci^Z5+|~*Dia_J9%#8A!oaEp}E9WmUo4RKM|HIOb9dt;KG4& z0*>q)50PsmfW^V&gRSWJf)>mx4z@i&8R49YDDS1_=uLIi)1IYH=w^JO@pAfHeOc>j zMg)2H6TuDtB7dDZ%}<*3cVl|(Y7=?7hc=si@vktA2EF0lkhKCqi1*B9QJfdJ^>!+Y z#x+WF-?@G`c7O5f1^z!Mnbi=08;n0W7_egr&E78~vup)~a?HIsP&^SzA{?{VXR}p0z6zPh&~* zML>!G^wYMwDe}q{nK>OhSNSAz#qq^*xHZinXA(~zr(NtKkD6U{UYKFN<+|F234$m` zfC0b+==q~l_x5g7bNQEr)MXCNcx^G7*sk1&x)e7fa>opM9mR}{k)4BYHizv$4+yHU z@M=gb%9KS@`kcl5abmu>D`9aYD1}uGN*srZODB9F94r@Z6~Nbhmmt=q0IKrO^zI?r zek%Su4xLEU=SgQ9WR67wFoN3#VrX@PfirA_$Q)01$V-%kld52X=uY~Uf$oYz7SwZp zD$%C#(h^gt6!ur$RSBzzr-bQ`O_JpMQ+TZ#);K7dUr@uJjGYdQ%%k3c1ePD&t_4_p ziCEUC1Yzj6{y&__SCDxM=;qzEGWSOGS`<_2jx-rCi1wLNPb#Lgdj~h|oGRSt>N)of zn-Y%leXD5whN~l&U723=`zPU`^+EnN3NsodTst zK;-`kEmnTT#++0O^CFi%sUE^)?CA;;E+sBgvN8F1VoDMZzL9f2)12e>eZ{drJ>TH^ zRQ$;ZWoV7_!vpQX-_DmPSxokOFdQTX5ueDQd(2QbWbIL0^=T{EuoApWf$AF`G=j;jVqSELe>%wd3-`E5H#*PBC z7XVD}VX;#Xh^ak(w2RUDW-PR1y&h>8ibYS>@u<-P%*>_J17YwMo7L?b*6qA@R;~u^ zZWS}tI_9_j{xiN7j3`&f!wio}@B}FM@ZCwNk0G-R(V*tsS5JP0g(W6%buJ7SVLZ&l z!Og|RrGc9pGv?q4zm7?9#|A)R#j3)cN7^o1)qUBl<)PsgL3Pk#U}eA8-`P&{76l5I zIN)z{$3agpZ24QK*`Z>e?0S)-UfaE?n#=tot*9Aa+{z!HNa}*0)x42FgAt>dYV&p% zpG{aiPHUa>6Le7gT1-8VUjdBtr;VXu8{l*;2(^XMdJodC*Z!J>5RF|3QamiV2|usL zHdtKXXaJv}Eui3>(;Px$(#HUSGCPV#(dGRc48R-NmiM!Y&mR^LAR}W#wV& zANjbM+OTF*%inF0A-TalRbZC5;6pa82jP9Stk8WTG1B;>ElIxTRG-3`zB65 z$}RbZE=B$$QZ^V!UiUGXF`Sd)&E+>Gyydw~8BXQK)<3V15FKn?N*SwQmT86T&J=mq z@mzNL`ho1o5y=P6pFlIJ4_rxwv#0mnu70oT_6|jVI_9Qj4;KxAzEBHU0YL7$+ zlmmmK|IXs;wDEN*(rx%m+%<()ZR17$CwOPy{kHMm zF~IdJzL?pa_HEzWoU?kEQoxYyiXb`4V#QBM10GM~jymrnZ*?WU7M z-jVq`XRZGo%^uewTN5MvvINK0lX+vky^ZM;JrV=(I;*eLg3VEu2jc@Wq-YZ8{PmOS zT%bPYZAiy`z%=|`%L7{&wL^Ee1G$rrt#CDmL?fkgK0q;1OLHG=N)9oVx4#ILG~zOGBKm<2%DDK|I9d-+<60F+BJ4S(eIOR z&d37g*tV`Rvt*PnWFt}c)EfCNzaO?-O@P2hYh9OW5gRQS@rQ`ojXSX(2i>LlC_?!; zGPZ!sOCYDS2hG(?)>{zPa@nH(9&lV^`Um^8mI9&0g*!Bjj!}4n%V>mq!DWkkS}d!m25No@8z*1m8OBJPx|kHW`PR zglqy9bm+OS^BSZhseMN+Z#dkD$8%&Uo&arsl?E>f;6M)NR8yS^VbocOQK&#s}@}RFyxF%^KDO-uA0&otL3X2gBlrk z-qEV*pEpjRHP*9L(_MT|?-T_-IO3enN^NbIRDbsMOwR51aQ;voAtq|>VYt_CL5rM- z4f+}k0QYn)00pEKZj0>Vyj4^+oRQaFG8@?cqv;L#P_=y&k~!QRA4$^HSJtP~)YsGB zYJM?e!f|D;p(<%`e1!yPoQ3{PuD>RO-%gjkgnD>QJ4`kExas3#5;2TYEnx+6BqB(i zyg(oZ^vXPtx+g12C)I%R4XMI;qot4L*Lf_I&s+*t4n1o@{pcoP`2v%clC2k%c+E)m z(MzQAe{n+Mr2N>_rTVzf8M*->fept1FDvKi%fmmiX_~xL65n0AB)??LwehKBSFPu{> z#Ii^Q?c2BtEJ3{9Ubo?cqYYd+JOd3;k%G*N@4dIjT2*muHs2RC>R$e7Aq;$%gDxD+ zB|^2%fi%z0F|WbCUVnDWcc;3GwtIi2Y=n}2M|==KZ!W+X&U3$d=lmZR;1NKMqE&j9 z#N}>d^S5d+@I#+a>WHV`L*Ne{{D8!7HxHTAk)`avTYoh^-2;trT(Rllcj;u1WLKgM zm9$u_<|+egy#kwQsso)eQfw_HwxQk*yY0g`IiEcj5;Z9K^--%E#ZcS23p(QS9e_!t z%to+SBo)a=$~@BR-(Y6jk_lU*&>qa4MzhArwWv@O4ZHCUb6f{E)Dv5?(dcQy#-YD! z)9Nsqg?$jZliv3;4JmqGz=F0f^TBk@i-|+LYb^E^jl!>+AAAY^ax602ZOD8e)(LE5 z>SyZbVPrDbecE@FW-n_+@NRNmj|;$;sN0*GItRRzxu3wo%%pBt?N4>=1S#SN@(+;E+*VbcbDO5alIXv99uAf1qOwXFR7}mA&hwFv780Y~)gP{0pvjGAPVa{CKd~+$-L=lKFUo<+1!J5Zo zWyQY9)rv|PDBulfy$Laim}8j)RUkl+!VuZR};teVh1YBc9O_|pGtb(7{F}(UhSH6;HY$BbzeO7(?DLTIIcyE$X(#6 zwJur-yj{tRvd&SyeM6ehXfC)vQ!TUqA!BFviU%B#A%eWxO4a9m;HUtx7>n#&W~S0KJu zoVDH{+SJU;WiLD6p#ZTylRJ>|l&1**L8L{HZ&kz#DFVJ5z{0P=T&TlHKaj4tnXnB# zrAvdDuAbUK=Rp9UYY0@|tc17l!P9elEYPE3vZD%Q*HKg_Knc^ee-~by#0AiSA!@!t zIS0)%K-b6yEc$M&??5a6B@2IHS1D1MbEkF-*XdUzqlt#5e%mu*57MBam|sq8YMQp4 zovcUsz1%!|uS_~uHP$>_rExTCat2>BZu&EPK*1`bRm+$0->5hrhfpLb10TTU5~-@8 zs?Vj=8y%FAe)x8cNPQ#4Fx$y%iF+qhfxPxy%8#CBD+Jb;mvre-eH2Zk`y{61HYB~B zWKOfKawPvPX)t>b>4n_Dei2yn7x*@#R!!+K(#5Wbjg_=}Yaw1*8F-*{5GeS+PyG7b zk*wwGwpH8KTJR%!-)EcBrjpz>n2BQCOB87~|1GW<*&^oP%7&&eIodQy>CuO-_I0P{ zY=HTYSg)kHmDOIJCW)xGN2nlu$E$Vig-KC!ko?GLcsDo<=vu@(o6rC9;WL|uZvQqC zUD5I^r8v&C+PHI>3G@(3|7IoygnacKkSXT)a2erHtjn?s z{H=cdPc*`F>hGom`;Q^bUhG-iHGr85%9;Q7V6~5dgZ(>3PQTPOFTUJGc(qi09Q~tR z3@mcg@3yF$KZi6G!macOoB$~eG%ELRAPZn1+X)LQRd|1k`88)P`j3HH1ggeYFe;jc zAS2*P(AYm&&`){b@A{3E;JS~M-Rtge0p>Jl{q3-eRHR><`^b)O0|L1Jzz%*8&yM z=zG0(1-|#H*I{C)>p(lZx%0=JmdANCmrW(72biH|L51!JDbo;1Ugq+M-;`Nb&*X$0(gC8Bl>vxp}7P zn5Vt{c!qR&YPh&Z6YJB@zSL1Tl=+ zp2UsfjetxodE3j80s(8swi;5`6A$h~%x*@#yBY_$8_dg|NWbM1YYzC9ai{wU6p#}* zcsm^3%=R1BS5NX=zUETv<)g;s0{aWRWT&%2RYh*@wjL0V$qtx7XWtPT1|B<#VSZ z6U03%d-)Smw0#OPA=6cLvd5%kzxl^zHr`FSsU;dmeWHnjeT%n(PtLP^o`2a`yO zkRrauumV?+bZkzR_!YVPI}c7o#n+a94UwNi?;gs2YM-jVr46fCy1Lgg+C3${ljc0> zA*Q##T9ZF>oii~#$s+Akfr5CcPgRn|r0~yR4M^7?0LC-BtDR`iqvy}%n8DGm} z`TojKb3dc)bO6-tg&`Nf~5HaqLxh!PY!6T?9~_)p+q68ek@bNsxZ7> z3aefU-*?QvEsU6aR!<+agsX{jm^c;QUlmTZ7to07#Zq#?-jrWm8rpxYYg^_4TKc>X zJ4p>Co%uJW$$No{S4-1O>%Zq~z@lKSvY7?G-fH1MrIrxVoO-^?vYE`CN+M8C>WOrq z1f)aD1M&AC{z-?;Y(@`O;mlvGU0&vcUS3rvzD zeTdNBildZMwIM(8?_kiM+s27e{BG-S$zsw>m1T=P0g<$gkFCi`Qvg-Et)z{J^$lDf z-|7t0u~y_e-E@bitR7J1V{Oqw$6&0vPfJ62-qjP>?SE3X(nlrhZuPR>YgeQ$M#K2x zAOJm|Bw2#=&Uv4hO`y?>bx?+9Cy3w0=R)sDI2?C~4Ugv%g+%;aLr?Y-C?^t_#3bak zt^WdF$$g;0Z-fkE8_Nq3J zanh>%n-=iWjp=Qk>Cma^=Dl}X+D>sb^8J3#1tA9DAxjjZse}0TRGfw|pOnUq{!a-d^Z~ryRF9nw^!R(+As$bI+I!6tCI}#*p1e9|d_;)QrcVfu< zY2rfY2$(R1@$qcjiLm@D?GrFNdl~*#`{5!+z6LXzXFD-7WlvAj993RAlXqoKdWs&7 zAD7$dL_+06MK&O$M$tLt?l)sxH>JO5(-EDn(j$3@+s=||+WN-&w94`N`t%_l^grsA z^kNG|+oF75zIqWRH~kM0v917och3qy!{K_YfhW=7X0v4)kL%X2$hk&^#=`4}StyG~ zXLv$3BX?KfMH0U0dbl=#uUECKNNW^WcZ)gQ3L>W-E?3SEK?N-!eu?CSZQ=9R7H#tM!Il(b*b|pEHTE4dTX}Wv7{c z=*>fVxjnK-mcYLsW3j(ie%hdlZL9Ll+4Vp4Wkg!X?PkgCwmYM)aoV=`0b7}F+~w~~ zQ+Sq*ZtU-t}b$(R#$ke>b$$##tA-}z6Nwyr$jeHCn28-(@4jy_Y;a_upi9vA%p9D6i_V<(!f^DEe_ z&c@`7-0p(-#W?7DwFibg^R`7x8gT%)^2TlenecewUN=T>?7sJ&5YSi@Wb`D!NOYt5 zHC1(9p{FOuaO3eE(Jbk$qw-I3jWtt@Bu*1WzDY>2tDFzHuF$bnn5cm3%LK|eQYMbP zEWeOKy1`e^>#H*aRZ7BQzoYeR)j-`PD8B0LAZBfPJXrb0fxo+$CqA}|OKrE&<%#}% z#nW+tn3N)b9&3MIJ1{NhmJ!Bgp}MG>&wy?{w-W#;9?HDT8?C+wjBLEoZ)LL8X~^O` zq7}gF`14nbQST=*D#sB4abbXoY5n&!dEjZC^<-bx0CAw~IdwJBZE9#JE$8vp#YVrN7A3y z2YP?n>9#Q138fN8zh4tDpINwP2ler^Q;t_ixN*hVo|1(yd=vBX=P=kq`S^ColMfPv zNgOoz$edae(RP`g6lHQBR)Gs5H0=M{mCRglAZ3S8X#)Y4J#ygTIrjj<;9q9(x*wn5 zyHSY3){^(|{#V$b!1Z)deB6MK?>V5|5W=8BKvI@8EIw$m&G`n|v~` zsumxBpi0*-z@US2*W>y(abyAbMMV-u6j#GjtX zfP&YaG{PExhR}608@i4%oUT6D!W>aAIr{nh1qbcL0k%HS zH`@>0{NoY$*U(rq2u(7I&XVd47hkYq1wCiNFF>v#EUqtHr9(CY=)mPlTj`17|TOYfWbXs6^G3D88T81y=I9hUSx*&v`G ztb}oK6tp-e6KEOKD8n9ta?HJ7{PP!(g$D%ok)HwwHAL&dI>X7xe9QwrhlKyeI71ra zu4srCDA);T9MivU2yS)v+`1dyzUE7w4ag0nJI3~ae)?&53z`Fr=+hSzW$6XWO+%ElZE07&I zNuHQ_cZ5|toZz8p3%C-%&bc=^x~bHcrVAX6oa|)Zbyzh_>E8E;_+$<{>AadGDmu zK?a$+l>FWBv*aZK19$7EiyfVE%a)eCoD+z0{y5f~UV#$Oqd$k+0fpGKsE7X?4hsD2DQ z2MVrK23+YMLQ)dJpH=Y>=MTRQJh=SP52tP$bG?_|%)7hsVrwkvW!f)^*PjP2biH)t z)|%Ia(JJ4r|3}$-M>Q2Ral=6rB8niWARQGDkS5Zr2qGXLARry-z4sPW5R|4=>77WI zUP1^+lPa>WDWV+qBLyrnV7RW1Cdx#B zHi&%W8g}%2@p?WUYkp?H6G*-}wYED6ejyB%1v;5x1d}j}Zb_y?n5G(%{o6|db#3kh z?Vg^+xkrXC4!0f)SLBz~aXk*|-WstUan1+offI|ZC+FfzKv0|D4KAWc)>};Ft;SfO zxLJz*Q;G&j_ZQvFI!=eHJ=>TXmBM>6)D*=m>HZ`AZvkU!E+ZDtlqJei%Z#DfmrBXI zV%!?!<6xKO=Mw%3Nw~1)_9RsdzH(O&)?zbisnXw_lBR9&3sL$7_)^OUd}>dM*ZC%$ zC6!3v9Qh|bMAS9Qr*8Ao#jI;5vR%e0F9u+cp*tMGOqM zOS>1Gy{Ap!B=J9 z+;Ko2jbNYk{u#QkT@IxR@1AoL(a!|qQnkA?*U_wcmP)?RVt7UFOo|mBn{LmS(bk^D zA#ObuuPP8V$yH}28(B{*t9O*Rxvu(d>{^&gT53xCrmZdHqLGchS1?~e`ibZPR7iKe zSM1`bbB^mo-IM?{2u3~ZdM3r=e)e;ZwNR^cWVOns_+@&0hF0M-n=0L6^+HYc8;2#Z z03yO>3%wnqpE=NXgdYOIM4k@gdrS4o*y`J zLP&#=!7n2l6^M$VE3TIZ*SHcF0r6jGo7Z6JT({dzQu6CvfkX8lFyzV4^Y;>I$U z9bD61VAa2y%{H_grX8fHq^WwDt;TUzzH|)(&``Z9ABC089vB!k zz^pH{pY?)WH+g6cLZ++)*Lzt%O;{;bue(iVYm3EGybUeel60e@^;|y4CPD~g)><@O2AMyYM<}s36LQ3%6+1@H%>rZDU$BD5FqLUDOpz!i2VcRPW$9ayN6DZ=rnCJlDgN~ z(Y9}0A~=j)-@Z;RC`)>BFNQp;xVR#BWOUljwO(H(VQS5-4ycyfy>)Oxx8dn)Fn3h^ z>gE?y?}26~qdw3l@hQy zSh7ri%DqxKRr>C9VIv9&YFzz10z;k`WoPE5kDNLqqaYch8@KFvRug~0Mi_(-r`=Td_`;U_Z7tDB&Q0nV772ZuJ+b_>`Cpx*?^%^#d72Y|;x$|QO`lQv3)^e&Ae zoi$xfbDOL}C3_>IsY*`PnkF%acf*3o(-CM7Bf_qI4BdFWglBWg#7YTRbCclpeQH-l z_2pK5X4zicVD)KxT8G=8yZZ;Yo$7og)ta0~?(S}*@pRLNH9`o#bE9Ts-e;bNm7PkZ z2*1I#8&C4t8K4!7m3C7wMnClL!n6`*+@p%JqWv!;V-gKZ=fhkCHI+ug^0aT0W#b>@Ls!rl<4H-(Q@sz!v0j> z<_UI@?A}g=TY7pebNe5o1nmNriqW6*(kLB$=Dq}@5E(0Re?PFvgq9n?Ab{;3DX#fA z4r);5)7-XVu^2JdkvK_{)OCQd^S7NrdrGYa2wnLewo9*h$EUcaWTdL}tKS#XfhwrJ z4H%U2O4f<#wX793&xbJdDck{MfLcjOF(z=;A3(g^59+{e0#zOx@hKURF=c4ryc{0LlScV5Wi0laA(v* z$$n48$kCrhZX=jr@(Z6u%))YS!!}uzj(Q9%9S}Q-8qqFPJt8}4Q^@lp{tN%J2?948 zC7b;yq>nT3kY}_hDxLywi}afLrcruT1CN8-3iP4EK|4R^ORQpPu7IRvMdxMjJSQQ3 zsnuHLSRfSmc+{1!wvEDbPNX|XW?*ve4S2#_6{TJ6-(Me+zGY_@r@Kr$vsGI$G1|`y zD^><)3Q!XXnj4p#&AiBlDtk=OBQ9{H@+pii3+WK(ppB&8a$D-f$!McQIsOJx^k zS45#a0vr7nQk4QuG~gl*xV9BP^olmDG`(T@j5ZyjU-H(iy0Y!Ax)>?t zE#sGVPIeah876fy-f(pErAcBwR#>LP#66og(;G08{i+fj3uoIBt$}mShY7wfF3(PU zK!I)kaQoi#%k6(JNH9=rzIRE3)gtOXU>x31oBCVO-&MOwCIvYfC5|fd?~X&|WR^rm z27LB5a9|FS$)_HVB@jOE5`Y?&>4Nv5v^yo75`L`}HmT@NCCBNq?d-qFTcpZZ$%}&7 zmLp69wrg;M$H=WM$(Z*T*q`{Rk|ET!T?gii=gg`_6$P0|NvYk8FT6M}gtmU$ub-Ym zdP4>@QRyDbS=+l34#=L3^z#V){PF!L$i-0~wtmRN_k{)^a_BcCX~cXdqmTjCfGl?i z*u$t6)!|knrnJwuGXte|(Vr0hI~d4%V~2(Biy<6OBkr)jH05Gps=p2!Uq_V!&udtZ zz(QDpz(iLvSaSK z-1v^@ByRIr{rG+86r+THDb_Sr!oMpnO6qZt*Up!fY6P}gRi%TSlF&>{&)dj!8ygwn7*#c6QPz zWr`ZkPs1Dbo#N(I!?R|1KzeDZ4uElOGg@`IlP+viBxKDo^03%SY8 zXP!rp>?6ntwgBePW5C;y6mL>!yLvJVVZ1QC@LlnD=wH@G)-NVVr2_3=vx=c3&?9pG z$hnG=!(aUingyEHbs_>-B`nQ7`k(D;;rLvMEV^rS)VgyMc;G=!-v?i*1_U2g90}!y z)m7B>n|4jD)B3`Z1H%+Pv%jh^qD@Gf8$?&0*6A0Iphb;&0i7%py4{nSk_uE;sb8$r z-dsBGS!Xx5L*Jf*9Tpn+lWuFNy?YcCs-(!v`VQ+2r3X-IT`BZqv$3APsfpHS?!Mu% zAK#zHwHxCDoYF8Qhh*`OXdf~YB;CKjbz`bunggEUTs%=hxC&NsQ4TTzg+ z+03nL=ixUGog2I`0V~E#e>XtPvb^N6r2PYki6Pc7hUDIn(}l=t`Jm;%{v8iR$mpd1 zfvSegSy_l<)W)A7xbM2n4*7!@xoQ}Jz(*0g@_4#M*=lJ2@Y-zC9!UsVd7mh-*|pm9 zGsDlyq=3)&dz@v!@_N{&Ji4$TwY8dt!<^ULR5Ia7(QE?Zn}6`MR*OE zS=W}9NO0-h^x+CBkz$@~=P@0>>-xOy*yyIz@m>Xwf-=cNs8L;eEIH%u3CyuU$M6?s zxQc2iRTQzJutY&$=gWCo<30J*-)C-464l8J?MYz{TD7tlGjk;Z3d+@dUXd3AMAZX2 z&sW5eauxl;ze#--A@ZP5AFdI1nNUDgo!*a;@c&AlwKcT>CFYe=S5)y>O7o9O@$Zzv zp-O|vZ6WTcO7pNBR?@`n{&OGP{+~NCVA$)CjS2~f_mI5*4iyH>EC@e%ZXNc}ZEq;5 z5fH8U3+?Y#`Zg|HyY&LjD-!ml>K8hRBLi$ksuB(%IpONL`aYc}GNGlC<4&8N)mKd} zN@|#okF8y*db`wl zFaC4Eq~&(NPppr>+)+EE5kLiJC$ivZ4LD4ZF=k>3~tn*l>nU@jg)sd-#TY26_I+c--m*nZ{p}-^z^BFonQF-guRba$IK(U z(e)DjQoClD!}|S#m2LmUQ4|v4btH1%R30O+#6xzry^1Q~xr-(9kV1PmkQsVlH8w83 zi;5mll}fnwi)ee#qy;QTQdenZk8I3S6lJ@(wyUQ5onGU7+7WeICYCBe((Pnxgt$A3 z8l|sy(Ue;wd4^7B?N2{B)84?HcW9&SCIYPb1A^5PC9KbMw2RXa`!^F`$0Mva_@&mD zQYe{Q`8}oFrW|(=d*6LkB^*znjzwaXu-fV(4%5w5XEtSh=cTp|CqZ(Iw{eW3-7&%m4N-9{TmzZr(=u`Rt5Rk zxgj#uBq7tmhWqm38*XlJ^oAwMXJ1YQLlpTEFDt2IA2ycfej@yYIt9kfyZE3)HmboL z60lWBCin!tcy_``#bTM4T(sO*9a&%A1#H$`oDpiXvWl!Tc0s*tXZE4?#qI~8G=uWx z=>kIY(hIam{{0qiYH)1|wYq@lp%qB+75fxw2@M-S93ziCUgG<@ngVGT9-< zM@%N`VbgE5KA9-570>4bWE(a3fw0L$AmFW|4W)0j%HZkD9yFT3CX$4 z9Dcm4PSpYs&<T6B>qy=aTRzYyo)?WlS!$cn}LA;r~TtE}gjOfaAW5Qvc9=Wx}Sm-;3V z|I?>A?N?Um@osQ34VM{G1T-HV5% zFaH6$NmrwBCD8C}w{I#>g(f{IDJ44kNzhM6D?2MIj?yAz4WhEHuBsUSf)~T#RU@ob zgPyVUa4CYQ_bcHFGX}Tg(G6?v((w% zHN%*uS&55|4tLQR38@n~e+ft?h-vh8BS8k0DsH=Mxbj){fbmnBVH?d-(8Zjk-<2rHX;WnC3GkqpgS<{VQLB~f*Z^Gj zoYhd~SPs16yqwdKCi3XRkNs&zh?M{2aHPkS>dplV<01i5`{ht|i_vE{xNk{l6>73( zH(NdX5q^qB-z#jhD|`3ko~y1+feLiX^81_nI|2OUV8JR64@Bbkc%~rm)hp?clOvQ7 z1RJ9&!~nknGX7)W-+L6Mk}AOrlqt%ptzCs;BwpJD<*)%~q9Ug!sz>8FOX_m-^L1>F z_v%iQXu0c-i)w`_w`i%wQ&l{NlxvD~>s}RW-KCn`Prxtnz#`&V|%4Va@TdtK@_TcU&- z4)-_2JntTg0*9Qwy-Cy(AvDuqvc~T0oX0J{lAC~%wvoMX6P@X;nXY!3$Vd0}(K^;u zl$9yiv^(=A_fi>FenyktjH|7!y}_-|6;!)lesCme;H);lm=5`JMhfB!1+dDC<<+Mp z>g-28sfN&@jjB$Cy$XxhTaar}9S6-JK-$wD6w62iqJC5%vi&=3RNt@p>#>)TZ}OI< z)C$`?GjlAIN7(X9|BCNnSx3#zIE%o7@ZZJ{%tAK9Kiibzc_jV&grw$Q;D`dp0ak1v zGZW?0y4S2!<=ez0i9IZ5E{Xc!_fHckC@4$er(Y|>-+p@T_;=$b5xiw!Y1KuP&LDLQ z9bx346k1bRS$&;ZbaG0_5AFQ%@5#pi1o-q50RnvdsTpwZe{9Y;MiJ&q-2t^etNYd& ze7tO0{*B%2Y4A7ppx5j%`9b=_^X`yC`)>H-oK zdT)#m=j$IHR%!t6vU)l$HdR`?zHDkl7pil9(Q8BH5Yy~+A~5zOMUchSFMjnqA|~!7 z&V@jFIEE{%I|d+!Q(uXXt2p)P($HUWSUa3b^#=lu1h00*mh4==WhtssK>KR#jyOJb zQXkf)8~9shbaM$@XNxupy%(%vGwh^BgcbDzFPs_ z_MTTIB!Ted?6L4t0e4*&`C^L_Y*?i^c;Md4g=0;!U&umYWXA^ZI zu3BkNt^0D~>n-$Ngbmw)q;nql5jKJM1gl})G}Gp!$WuWGG-i<0*i3klgZKnyL=YwP z#_gxs`RLP#2cExc*)3Bg1NvgWrnWhJG(EN!&m0kJ-P!Y4u`cf9UFkV_@dFy-?#$Fq z{~a;pS^?!r1~wIoYPMvb4wDwACUawPUZAK7>~ZV!)Z`8D=olY&bo}CGiiN8V5xS!U{Ml$d8|vR zb>E6na}WIrJpcS!UVCpq7kKP|8|{wm6iLiHGbqm;NovtE@DRE?a9H1d%rz|J~XH0!2DXD|hXt=e94y z=-|t=iNJUA*lo>=jK_a9#TI})NbTm` zK){>~v%{tJ!mI#s>U^oLk1aNYuAN!Ubfo35u~AbN#ZdK^ioJ1oRc?5^8dYxbT9F|^ z^Kgzz-FCKC?#Ql$ibLVkf6+_dDJewuIzxT~}Qf_xB3E1~J1l!nbZ~ z!(aS&7VB*F{EIyj^+%CQ;_nQr3Ns!zn^d#D`$W=q?5#oO2k%(O@h6mSy`(H>P~F*0 zlj%r%Y)2nK)LWY}U|6ck8Ze6R9wQhv!OldbJz~tre}CBU*vvpx_yjphLn8G z&g6@cysO>jh}!qlMPK6NLhMPo)Pa;ykgU`tjh^U%;FRaDJT+=Jk2+y}t&umVpF-=( zFF(Bw<%arBr?3+v#@E*7r@t&v`U6*jc>CyAPi8O*_hJnq|yOoEGXhRl{#6 z5p}^@4*VRNB{vK~8n>_~Y7UMz04!yS@#nOYSlE`nMYAU3}ei z-St1_I!Pp7x^EuhY2wDeOq@3`H6Q6cCK{ib>sjWb^wczC?823LXI4D{|AzfIc;9pLo)aQ{eI{x=A1MP?b%jCe zUSw($k96&SU0=4)b8#`T+NkN?Ql^hn*YT^hC&!k<1z}w#;Gc0zcLjQTizTRPT03$9 zv%ZNE7;yKcu%W(JmVy{4{84~p#&IW9ln+czXZ|fTc<@G3jemX0RCSucNL_b+o0QOT z=y^J#W4^iIYm%DI8b_RBXh!!It=~I<7$0hNxr1(IQ=>qU4t($6F!vS5Izt0pgBYtS zYGop3SzWF5*{MTMs?u{0=ZV=C;6TU!-qxdU2ve=_=`QtXMEBBTNAURO(QjBy?K3t` zekQZ)?0dd~+uxEnidW|4+Hy_Hp8{B2SIcV1vk*|Gpuoqncw4XO&wE0z*Ua9Gy|$i{ z7)>8ry}~OcDMhc9bk*$@CrE}5ZvE-(IuBYnrC=7x8f}}wvi-@_w1T{X3^oNzauRb}dkst1o!Os$mad+|aL-E+}24EN3$(*2FhJGEiw*_0K78aC(zrC#A% znQa(Hb|3}mf^=RZOFLF;KL$=D zIBwxNMh8=;Ioq|@4!HKEXxCcH z=Ttf}%L1r`TUfYq1g2-TxBRu99oqMkt2mnG*6dYV)i9;u`$GrquZ=_@gDDDJ?!7Mz zzlB6`jacn=S3p&N_^tO`=Vz;@>tlX0d_-KhUe(-_)9SsFT4uM^wb^MZ!QD?@LdK-q zc}W$*+k-D(m7~}NuD)NF)pTSIYR^&Cm-O9QuZUpH?%(!WAC;e?u~&})F>O8;Sf8#^ zQcaW;06-&oR%XuTr-V^@ZgD079{rzHNuAoIDMHpVgFlC3&1~nE4btJbxi{x0aWqx# zZAE+&_xU>_=p&B^WEn3rKbXRu_7#FeKaA*s=P+Ipn#v*5kqvtjx-oq7Xh0*qYTCiW zGJOPGYp$lCYOQYKdlhzS=U{AS+~`+;FI<$V#;c(i=|$%(fi1 zkst?^uR=c}?!wUaUIz=zXNZhZ8&(7OeolB#?V~V}0o^L78&+NqZ9g+3VLWm~6;S7p zU7T5!E(en3!{1T+An*Qvxd0$0y0RlN^&o|#dz)AyaswF{BCPT8<8>>0E4$9u%83wU zR)Y>@>I)w~^sjX!2$_>$ZGcPBiHSsJ*nbr&t;2Q|oSSoZI73>*Dgn-f{k8Of84E4w zJ1+is9%0k12~*a{wMzG#sNRM{q-|cRBy9POrjmWkun8!O{NGKF0yo|F8|vd}iV_O^ zd1n$7Fc}^b{W3QvuN6pOhb-u&$J+5kCsDnWp>2X%Sg=NQ{Yu*Mj;@Z-CaAj4k)df6|j%?b`1HruGU8N|41@~he1&T0fi0*<0eevS)$ zD;|3Iu<|)({SsM!Z66KK+dx}aM?}OkUyBkh{<^Sw4h%W|LLZuy-8%i;K=m& z9>meQPoXYX3`6=~Snh_MfJ-Z~X^Gkl(ig zz;;TZCwPli`p!log3HF<#{u0M0~L{wjOr*HEwD@%K^V7qU}h)xz~JLE+|9zypJQv4 zHn|*+Wi1mls4P=pcF%FA_6!MnHXm?*dVFq!wVny>AYl-r&dkZvmpXb0lV&|hPqi4p zqYMS{x*+Liv*CEHEK99&E?M>);p4f`;0RLVaYJK7EV;HDA5h6`A4#}kO}E}da;m&N zO2pfq1pVD`SGYdP=NuotHzC7nb4wF`ai+IDssh+q8UK=MWf>+hUnriBhymhhU2OXbC&ngHYnds^{avAw`B`bR5N$0JTK4=X5`(~tZU6b7 zfbIgxG63k6`9X=#jXrS#a>DS1L9&`^0LxAdPUgM&!xb-_sL$~Bao&?83GcqePuE!H zAe&#vWv+fiosF${UIIZ|rNYy60?_%Ta%FlTvdFQI3AKK^O6L+c{F7G>*l9ER5%!H* zFfh?%5^%%=t;h0Cak1^4VvpkI(5JMnc@aD(KBfIzak@mHEK+=s4a8LLyxHj#)*6A7 zU+AF(mAxn#x|-q;rzbSAarq|NA5DK4w{&rMuf=LuZcbEc*x*HUkl4S5IzX-drYM79082fb1mf%} zoj<-Oq@~7x`>wC!Ha`0I0bnoyv}2mf^g_N5XZ}x2a(qI9CTsYSy#D|kc+SPw9Y9!F zxw{ox>J>oeTQkp7_YWjv5uLEstd;SD$mJJ&BmCv6d=eq*@jX&M?|7vCwL&|K|0Q-3|)H6$4+vxjUZ6Z;a%QRDID38!&V7_1n8dKCwwa zgd*_Qzhvv6P$a&s>B^3aNEhuDL0|apy}Jp^X6Ek(qXB5){%j&A511;Z2J%K``!>xL z8m_$4eck*Iz>M7HN!l0c^fhkqfC7M|$)eOrj5VCzUb*rPM2$%RpA$GISN()r)#m1d zmU7#^nmpFk-n_%`M;bCI119A)?vwVZ=xP;TJ!o;KmTPTZMODiHa;@6aeNuF@S|y(a zv3C}U2w}KmH(fTz$hdnUCE5I|LO3ET^n z3))#N+nl^XF&GoBrJ?T=ZCH9^CcIHwD@LW{+YRPT`_xp4F0YXfU8ASVk;~PQ%Yw{7 z0xZZRx>AAF638L3OP%;LA#!l+=40gXGTSv;&%-Jpb75dEoKq{+xl<7UR!$jbL)s`d zNJri)JJj2qOcrg9Cky;~7|mCtbQzos=gqZi-jc9|d_BfoVAlbT-Ni2}MZzu76pRU0 z`I8#Dznm>k(gyqDJ@UP!i;kYu(e!O_fUGaws!^c%#j^jdj6NKFC09`a z7pe1o1PG1aW?KMg$Oit^17?5$4kD7y9eG{va?+Y}=Tlp(xG2}Gn`VO6GNiiLODS{s zEmPO+M$?b4i*Ncgia^U@;TvwcaFxO#>lb9NRS+!ffDvb+Qf_m2(LFDjj=;zTlgs2h zO2XK4o-Q7dDK<_H_45J|*91u}W`}?@z}tuLvLrNm)tV9d2bir_IecrB3H^EaQ7Gh& z1V9Ji`bc39onZ-FTx0xoXcl6l_q!0$z9c9+-hu8B*Zw%a6qLWv~wqK+YI zY-H(ep=NsH;dTjG&7l-Euq*F&q4u;6tocD!IH0uj zLVy4i>ZTYW>bW&2SCI!@HhZM2jIE(M#{_!{Up(>%0%~FV$P)l5gT}%^w`4@DhabB} z3vf3_r`XxsS(UXX8_-8;Lr$(h#~#$*Ze(sQcNra^A|%}^X%|Gy+;^(|)an3|-dsyW z1e24pm@_t1MxN0>{9Wv?0na7^lB0P2WUl$O4Qkv~1Re{wtc~JXqrDuD)NaGfD6;@< z=&P?~NwWETq1G;+wB8ST(@su@Y#0#Rq_)MCd~)Be*Mi7v|6RKKA#fNF6SA3sC{(Ob z(#6D>I|)LrSh>G9tLq~UewyABr%ZV9O`)BId{5szPM9S0+Z%Ws2LVKLfm-}_M(I&btZW9eufeY&2@$`lJq3lib50_;V zZ{7gB=+J@*_gbU@a@GFV2}d$uWE-6LZmiRA`FzqUEu8iRd~>n|%+>%kCg?BTFH0o- zPp1Gg9Jmm_$v7MFy3U^^GVfAG7{X5NvO8BR#M0`I2C85Fv#H11BJr@(_}xta@xvMh zr~vZ~@XS$N0v2BV$TI&BAE8d#bNqG((fz|Qyy2&gD7(FO=e0ur4i??M`n9adK@LwU zdSHR!tIh%a@lN~L1Eadx1Ngbn(RC355mYGOkstXBc`2gVt6^7w?FtI~S83s4h9!`( zVfqa>*4SZQNP4O^#0Jh)hiY0MP7bGp|L(6=;sQD(Kfj#LrdQz9dFB1~Z-CS3!{_Q2 zZr{EiYBFGEGwvE~sHTnvB}dq;DH`2J%fM)5#p-})DDaXVcb1c@#F77ALIZKd+N(-)J}USuBEE?KwC&ItS+w7NvpT|kyQtvj{kS-sJ@y|Ik-UT= zKwIhI%k>BVnSIZ?l>Y*rg8CBdIAu*oLTr4kZ_BOAT(7Kqmsk+=;0NHVDY8-5y%`sw z45SaS@0F3X>*M?pA@-`=TgMDujj3>Ry%we>Z zPc)33yE%)W<{3>KY6XM!1CnC4tJ$^MC8?W4<+tN?WpBpYY0otDkFnp8lOLv#i7 zAaHEEa;I<_k^4)NElR_7tG61jsTDvf4M01}JHG+cD?ko%dTrJ{zn{|{lt7zaFQ06z z5?Vv&JYzH1dDq1=0X8%`Ffln`y#W~D`ND4|kY*cp)}EE_l|nT~^~Qe~5x8WTOmdk` zlwlng(m$k{Jo}Iwk?{H}Tgykq1t;G=7A47IpWDn@`5R z!8wS(!F3kyB>Z9~0}(wcZ}F76%>5v7tCNd*;=sE~N|ZW~A++kAZ{x3pq}KQ-@(j6h zfT_?`Z^(B#C{Z=S{_<0`mTCc%&lSZ8pUvE`l}LKAfYC?tuo+`PtDo)@$qvx~W>lTKOY35ig4e zehSzCNk<+q35TeTEp+V+LsZ95IilP})p}s%M7^63L{2)a+nZCUo%YOTxc0vHjC_U@ zSQv0&g278W0H6SzvrG`eCQz%@@wNiJHe~5~yww-jLVtLAns@^w z9llx^2sxc3wkcjlv+aEEG-DA9yFkw#RDagWZv%vjQ&D6eibv`gTK!UT;YuKmyb(Ye z$uTs$$A6{)ksJ2V;nRJ@#MdR5b6=7U}cDL}xToN-5~ z0&W@+*XO;N1@XWpj=LPxl>mrFfHe+kG{%0ez7(@mIo4$9G`DtIryRB8df7Vaf;yEz`z|Gr9)XzdZHe>+aFTfO1Kk- zp8^*U)&k;|c;AxXkJQA0R0Jh|jHB0Pcwv=InNBN490Ibln;sv`rafbN>4gKkNQI4<2+jCe0&ZV^jBmA z9VgCal0EC{R!9PWwEmx7Em9;5^Rzi5>NCinGWY4bxp=6=gNY=#YgD*>Ppdm>f!5Uj zYrULrfAscQouVz=UGTkUoTvL{7fMeXm2-65d!PMF_D~lI+3_?(zB!1g;6Hi~#3Y(u zxxSacJt8pSiLiVlm5v!WbO%RX)-AA8*oAs3H z>7n0p6HYFq&G`)MOVC64d@U~6@>w9t(bnoa%?*$mK53ze+)A16$H%Xm)|4(5D*6Gv zC#I}R0JV$+kOiRM&(E@~Z>XeJ%%1V@TPRlwR7WES{<(KvWEAH3?5w~@?3E*Wi#bAk zH;WWRh7RJ5e4BqLh%`z!O|0z-&s;_bK`mK@MU75?2QuOhLmb~lgPeB-t0;Y69*2iD zDAc0d+}T^*6j-Z!YSscVCAxfk?B2sce9sgAu0He0&3c9~o08Z2T60IS7r*{K2diiI z8Nd!J@3XN)C@Ls_L7jW!+G`})Ck!76kbxeM0eU>}ea@iGAdj$!2>W{%AD_O;>kKt) zES=?mwZ=1MpA`tq#%$0(G561IViEf)<#!^r2zVAhJ4X2my+rT7*gPKW(Z-Z1Ac+Ld|P@=AFMWIb>sB>`lzsW z$ij$G=$XpFM9Zd+qE)N%Z<>dX~)O_oD*OY(Kf=8m=6in7h;s zZIOb@7Uy$C*p?W^H8)*qb4KyckBEB~BeayT3_-i6{w}5JS}aTaKLiOqyr0|-1UAAC zLM43r3f@0m)sq|-SE7|4ME;yvlwJisA?$hbu`7DN`emd90E(M-qkB;Ci^pLQH}9p! zrwJ|WRn(xMYxwoY_O5A#%zF3hFw}tK(kdjKteq!q!sYZ>NL9}bM^nS zM!IL(QpX8`c>cS~e5gDJe-Ktx9UUp#Q~~~PN|~#{yVpJ)B&6&AdM{R!ML=Ht-+j88 z3d_|f(eeHzvY`nv9i38zZ0CR^)tt ziM2x4t(=>NI)@Kz$;Df~MD@=tHm#l-to~t*Z>)El!4H<=fAa1dzA((cS76ZpHF6M{Z$Zj$ zhlBZ+QoiX?hs~d+5(j!&DwoY*TVK4GIdAkXod$RW6sR;cCZJ1b~8 zHldvcsjJD?a}WXPN033z`F&EC*;FNjBm`M%JiO~xT**rHjj9W5ju7IwsH`X;8)t$9zh=BR3S{zM^e=B$B~{XHevyj`4M-PR^INQ zJT5F5df3^2FrF3R*H7`GOKUiBXie1zc5VF5$YMhd+ltL6e zCI?57i7w?i{>sx$%6jg?)f`BAE=Vx=&ynI=efc zs3>-F@_fCwTN7|T8GP-e$rJ;}4Sr|0e1g--9@r?~o`_C}DSN=)79DSaForKQ$?ljoOf~?;u`{6~0zwJ_= z3;mrwpS)-&7}}6L=%HeLsn6qgG(ml#`_$*ApS@w40GC4H*nK+-pnP9w&*5TUJ({X< z-vC#ne+55l?rxGFOFMgY${gax7BqjpC1HLbV|C|iWo1rZy~n=ild$ET)R7JTFxO!6 z>{NfmNq>Tg6iW<(K%WdKGyNC1o3R&_fh%T~!zJ7md6OeVspK*&vKB90R`P+%oX(72 zR0!)fjcLa&Et+3vmX-Z^N9LGN>a456pX*kgwLT>SP@c0kGUts~Ua{k`nf%n=Py}UP zVq)F=ydPjigocuDdzUs_Adea)+BUS@$2@msi>Mz3x0JM}eSX0fqRIS1x(hfVDrrlz zP%u_GO8I4afVD#P%*8|qu^855JwcTuxPz81B1AQELsnDS{Af#W%ubh|pO2qUoLh(l zw8KsZCL-$SKVIopBsseG|FQPoQB8GG`ylG8bU~#fASxnVYNRU)hzKZ1k*@UKi?qN? zRRjbCq>D7^y%U)bxPu=_3&)#Pr!M@u1 zeirG6wka;e1X@p3wRWA!Mb~Q+^d&`YUoR~iTe?hyP&)1Dr}xxbAeQ_J8|ZKvc~m4S zqCve^mgVKzZ%01U`uJmb?HeUz^M^q8c%O)~Y(~f=bMgjhH+#Fy%79fQY&zLj-yJROtMqKV1yG|Wj| zy~b_Y6Nm8Owv*k;f>&KlrB9&yP>|0Z*it6BEVi6omS34a_^%GqnWq$FH@QWi@bvNP zI+w*`bujq?=m<2!P%r{8g=UT;TLyKJ&)O$5#;v1TqO=P%a@DWJSRV8+KIhUrNZKMh zJaP1#eDD_fGEJjL`8U~%7tYSf7Uq_A=9bb^e|&6BwToUK=?XeB&Tu-65EB@#YK>Cm zk`D5EwHgc2gDO$|a^!d{qDm#RFoO4zVCfEQQ zk@H8~{u=CPAvV_~URKm^yoLq3S+ml^&f3d2*ZJ9l2Wyx34(hUB@C8 zvt(SbaD!KC&f1AR@5;0}clY{+t5{BVxJdin(c@_;-@LR@D@Cd5^}%bOmm$oEVb?t%F908r&^LDddBTZg{*`fV@H{5$Zg3usF;Tp(t91EJIzbJ z-i{2BrW2q$P_MI02p>K&LgtJ4ewD@OKYB=+WQDYo^6JW=UR>nLa~&f?lPQa`FDwDg zPDAL1?_!ho>dI%)rfO@#G9gn{*1X(ZH)Nireb)HswcMF)_HH-ak$^uG(?CmedMI(y zq<1CyZN~gYDAq# zeK@Be4X^b#(4zD=;jxy|+aJ&r;#b`3f%i)5DITwbzc-P#h~-|g1T||n>7AJbn;4+( z7e$`g(^h~ygTh8ndNGar<^vN^zvm{~p0F>9fM=?8|2``OZtB5QnH2GL^LBuSAnn$Z zU6ng*KF>^BAh<{??io!VtIHIKFqA9AbDY|kWy!^1+CkVEH_bWp;?6RMAKXjSpjoWT zGQE9mc8`$Uqd*v-i#3j#)nBw%|6yBpF6z_M{lenix=u<*u5Fc@`~=^&lKXf71+KzR z(B7dV=+&OzHFqMP?Z@QAAAPT2lU;vK%dC6P%PA$%N<{at>(`e$cN>cx4=Ld1G;}x{ zD=_f!UEdTZDx^9sOSF5I4TM%IWUUk@4D1#-!gJu9NrHk3og!tO?pWFMSDN*O<(ux8 zSXclzw0RCDi#uMhiPj3ioBBh#h2tYratD3{Qc6MlQ>=ULcg?ok&~tVC^=)&Vg~3h9 z)VTh1XXp5$0!QVixtqKaU=@yqpWEx9+%pNpwWjKF1u9RkVEPu9&E(;Q(v-hDdg#oN zpewsJBRjmF6gK_J@A|e!Z1K%*?B4UekUh->ft& zhQ)${)S@>{%M`%2oFL%0Nu(I-1S8VI+~{V={eVAD$Q;w&)}B=O>9JICJp@|3<^s+xfL%_l!wk3`?5J?njDE2_!|snUJAS0yAV zeVw?x(k+f^xhR(x=}{J#tqYcY-dTJ+f$uu+brSgqA2eoJEle&@ThhVSt5=u@QuXPd zDSfx)Z$gsLQE|+H+iRC#@3ZK*$$|XMWiCU-nsu4Uye+a z`!sRf>Lom}NqtiQ%Pw)lXW}u;~cTmRgG+Lh#_%Kfu4W$fGrkwY7 zC`+YIPq5Fvl@FYNIOqoAnOE$DfXhcV2z%LZY*JKoiuUKsv31jIGrss1+8-mS$D_*? z9-Bpt=)GPw8`=ijDg8+O&lMtq;==5bY~oXf{JqM3axGCjt%WxR13A3>dQ-kcdiAkO z+t0gx2M}wQf!>)#1w9#Pk-~_05{~HTT+w_CH3a} z+IxjA9k4cOmCvgJTj_C;Nx=uB(oyP3(NQs|96N92oY&8zhSkTgT~2&@ zr7x|3lSy1e7SgA>u7;>9GpK%#&215vMm6l$e+q=Olcx<>32sd!MW(1GB!Po<4HZgy z`f8Ol$~fs(q^8~s1Ml3UqWblyv>!8y0j;*Oab-}}1DAGYmg_W!uk+%L;<@?D!rX;y*8Unb37 z`L`V5p5tIf@i0oG8`*#hEv!Xd#!{2Iq0C7%bJVYK2pqTf?D2Uzal`hm)&H^XzW-Z+ zuOlHoX<)7VUf6C6Z;iB<&LdU0Af)kGA3sYUi}f)_PGQ{SB7AVOfd+0P zJH(<3vU}*Zs2DD;?tlGe%r#lisRRzdrBy#4#o=Z0+)`o062E<-*4*1vTeZD&)W92Q zW|qhlH-S{qb=(^)-=~f-kqK$46X|iJ) zQN=uGlf6AZFUXKQ`b|CYcU$2x3W=L)R1!UU473aZh|G517Z>iFcr%9b=fe>~QfC~0 zxJ-P3=Q))Mip0b@;)p57h!9OC9Ei#WXM=)HOA&^Bx%6rC2SltA`bF*Yib|D`n?CCu zDo`X#Vk~?~B}rx0gaXOM(W|1MR5{M=k-YWYgjXJChuvFjh)q$n{g=)G3-XGnXO44$ z-6q&a503A_l&$)6B?Y&Zu_)ZS)hwJ2RIKUx2jF`OJGnF9iEuw!p=!)I+YHAWT6q6SVqC5OV~eewd)tkal(_l} zoh%`T|3*pQhw$`iz54cCBo$%s+e6Y!eldu$E()bd%KAd?iL0xevTIO}sc@cf<^Do^ zNvD&Zcet8O0npDV{In^4_+7rSY9zm}vu7?Xu9}2)3K|ueGVHvu&*F3PQ6u|wxJpV| zkpHwN!Y3X@{)++VTyClP1h1H$3;}Se2uYYW(X3B}=n82^4^T-2)x4m%KuWOIQ&yhv zi5`dLQA!dtqSyn%7gc79IQ#0A0d(d7Aai_yO1=I3`Uv zPf;eVSSDxx^1o3phqvGo9q6Ru&W zE4ZEzPj^paj^x?wwsZTV-h`M~W1an&5X!V?#%K(hs))xSX>0qv7by~;wZCjnlfZ3Y zXoRQP_F3<*pbH4y3bYz?^9rXR$S~gOqnSd#@G?@ne@&E0f-x7jf+`C?bgqVktl+cB zK_JMX1TtYvrBJ;G=geOB=7pK#a|bek?;@MWp2Gf$o4#@^8&8fYBDREe`S6?5NY+R1 zi0RD@mXwm^;nr6td#03Zy z1z>FPSNH4vv^clvQnu{721IX1t%BRL9ZOZSlcM{)cC92%ug9sRX>c4&6vOidtQ46i zjOg6P8ogfI5;eK}#|d^P$QALN_D_TS&eWXh`_#a0`Xzp56X1(apd#r*T}|DovCGw> z_BQF=ROoy0rYgO5_1V->UE>7;ho>A>XB8p1)vRC{BTY#s0r9}bINd4C__UxRwRc%Z zd+?$^h~WJ7;S5(26i==?sR`$dPM7)jmJmmG8uEj@Wlpv#T4{gEvnvnW{SougaNrAk zat;{?3~NZ0e+S-Mredhjd?pFI)Ens#JxNF+~xml?vdBjPc2h1LmT_1Aoj{gQn=%wTA3=` z!0JmP$Ts^=#pF5l>ATn=xUfrBb1&fGoq*AAWgGf@5D?)+{Sh`QP$PWaXb@UTp85BF z@Si#2bDur<`oGF97JrtzcVFm2kL7rFIj^!8dcpY3#-#fW^?oX9& zP4`Jk{F(cK+5CX#_-;IV^tYA9=lYtX6ZOT|({CV3|8uhbXJK(k&v>R2HNssFEn^9c&g!psxe>GX64rrTe-`iJ>Ar zZGFKGVGd1Ku%+4oijVtpn7+<7R!22zQ&ii?QrY9VtH&)mvf!Ur(pm5HNJ)ybMBTpP zubiifC!ioYAS*4%bXL*5i9bsm$m{OwJTvR0BpUe!)Z3Q@1~+@!8TFw@D}k3!{-eK^ z7cOFM=~FiNF+NHk%bG6fJcBNJRl_8Kb1hZ5)E~=*Jli`f4Wuw-D>?ROr4n#|RHWH{ z!iQ}$6!F5*Do52 zi!!#8aNOOFFd%+Wv*>z07n`W7fB0r+X|=PP22jS+cj{$N7#yy-Ww&Q7pdM`i5d+a0 z935#f)ir!-Tq#R*MZ>Y=9zomJME!q!In2e~9O#BuU9`gJ2Xt#^1Vn z-(PZE>|G6!zJK&2^3gMrfA08wp2VNB0WKPl!S4M`cWVDQk4U0y_D~Ac`@7k8y8$M} z)6&)N&98g^apk^?y|pkI6MiCAXR(4F(~>tlvXvQgU8P3SK8pvbSbbXt{{NOL0l~BT z&uh*n3wXt^zX#WXuebg`@a_M%>+*k{%KzUkI+wckb&YV7OBH(dQta?ug{UR7Nk=QZ z?gN)0E3hP*R6Bqr;fc`*lZ%dt6p~tzfKS6czME8=8RkQu+EP{XR{CsL0(;vc$5q~F ztP_#4;tr-gm=m{_ez+t^V=(OpmLq^}kOttp3XIt;ycNse%r7`TJR~kAH7f~86Ij`^ zkj?GK)&2u@?}Xi?2e#dTZp6VNP@aBH4^vxk|^6zbMZoQ zB8B^@;}EDQKkNFmjmaXHqXv!jgc{trycdFu#><{(l;<*{DQa+~Vo7futUn+oE~c}= zZi&NgS-E;>qz@RI-VLJ=!|=K=pUkV^0-LSf7i;3gR4RCiMeQ?%Jf_yuzZ#NOdp7MSWm;H}p|%AgYE2LGx9*U>5j z>I25c)W7werDxfqoAjZqre57-TT={>KyT323Mz`c3eS6Zc>3v_#bV~H6;)2tqP=;- zYd*y9aZyOP7H6kLBb$K;(6WDvxmqWQJKI3eeQB1_VX5V+er#mT_tF>mBo=*qFS)2J zEKuRmHjQ2=Gf)J2nz-6yfG7n|hl>is9$-%Rx94$txy{yICyA-*_o;yBLf{ZMk;Mf8 zCBP2&!`3I!V>L)1$?k%v$T|I^(S!A^8a6J-kg<^N{yojUCN4aNoN(%VJTn1^Z}>EM z%{cu!DXB(2AL-#qPN&A$LXt+2u(@zU%@Ll~_Q8VAn6$@-Ar?<=+bjN(9Gh7YRjMwj zNr(O9UrN}2Uq4-95lHC~%}Pg-7|qQYcH9DQ@{_$A$YuvNs(5P-r9~H=qSH9a#MmL? zx1ZjpLIH3^S%>$BW)1E*8F|@&m=4=-_P6NhZc5L)#l=vq@??d1S|)~%VwPZzPT=)B zAAvpkPus_fE~RfwL;|MtW5#FBkt6wqP?drqukG&C*Kf#ZLdW0n!}fwgzOK-p42-3? zsPr_Tg@E!opq50Rb`T;*84o1Zfl+o>1dm1Ut)oa{7B_X=K$hB8>r$=Xfi#`~a+k2+ zjw@2JKC=#B#MXOaBOSM$C}=hE_m%p5mOt3x8Tv;u=U}U8_xq?t7jwhr)~Wa%DQSF= z=iR$o`0bDYhRKpB3Z7xR#h{O%IXdB~s2FQ`7Pk^gu#q_f#5Jto ztq=d6mA%L{(a;wJqh4f7L{7ZYP@f*^Ura=I;eqxqp}*vp)iT`ejl!$v%580u6~GT> z0hp~@g@2ihn$Sw5pbz*xa{Apwb#~>dd9vNKnNkJvq#v7de6x9b^-6}uBX?TK!jZX`|_0UM>+`! z$i2EDvU`dij;60%;R8cw3**0AZrxcy)rwycN0=u)xNoC_|g97UrPl1MlAauzy{OEy@d(!VZJGw#r-4^8mSh=^6Gv@@&CYr z1jV{HPv(NFro9y((4l&J2=- zN>S{cJo@uk{)yRgPNiZ4HRSmq6C*kw#sRyfW5yr!^R55Dx9y`VVi;|NHN7j@r-O(N z8~kX3>(<-M$v0(V4Drwv1mAuytd!Sjj`T#u0}(mCyD7@>btlX>}CVE>wnKhG3!d|#s#*PD2Yr>J`uuH9Uc+h3vi@!8Uv`lAMf#Y@R2iEj zPmRPK!JTkFE|hjRjlRQM5(7=$FE9Z2RE45;XVhj)MOi1?qX`& zYV_I+@KB`>aeg83f1It(!9$}4PvikUM1b?|Emoczd<2vT-uGnnTOc4nS*U2LA#Dj| z1DqrxEYhb?><&h)_mC)PcNTNIHO7P*z4<_=N6vrRd>Cx9WDRIvKk)F%QK&>nVWi0$ zC0VNo-?biHZOEGOTJ?E`{cY2;yGW@kDUm4)317iTUMXmU zKyWRZ4UDvLP+f<^(d-GKJSTP2pKJ&*G^Wk|fpaH>)zl*DHyoE*C+H{`CI@T&21&u0 z28}ISi^Dx@%Qg=!YJgBUyY7Bj;otU0Z|ogix_vu9=6+zu3ze;WD=6y{@?-Vfb^+J3 zh=9yzx9%QYx@PclmvHeOC1E~%Om{bteoHgazoDtl0DuMj^?7?(tjFKry{^Te4!3D& zD@a>OW8b;?ekW(y=%_G+>naeTXP2)@3m&Gv$8_YkzAOEzL=#%mkz!o>W(|6({S*0a zf6nSZKnuX=f6fH(B&Eto?X7J8L>L{tuiv^F130rH9G0)>myyC#}H#al?P9|?faUva3PaG zv1?hf9Ii)bbGK%2?JeM7E??QV-Op7AXT1%V@ZkjoyZ>-DDcxK0MGl{nwrXk3fL?Bg~GFyLrK3b5Q*dKt}Bdh;iH!`?Bflt$b!lachMV z#NL8GHAJ)WW=V@ z16kcf=OAk?h$?molN^3DK3RX0HDHVu_+tTN1jdVFIw-3(`zf`twq-qGF3(R>u&`OwU4K|W{8c$FO&DdegFQw?4wwKC$td6| zRc?@}svNd!a09GzPgN`72<*>j9h6-%ox_d^(Pu+Yny*%5z0=fEBE+^~d9*aQG^a#C z3i(f=?53o+-O-FqUT&9XVy*X{{+U0v`_*o^i%H%`r>>7<$U;4sgYXgKv%6qCALVcF zhl?;t31>x6EiHkNq39EG-O%2u1U)&IgQc#pIxAq-;dj~+)F&=)n0H4WUsL~6IARt@ z6V>SV6d&vPi<09evz@(AgGUs2^bQAi&%>Zv$+8%A>N6pLzJAJ(u{xr5C%0JY)}rQ@ z{Bix(IMzoO@ixl4t4EAlRI`sBJb0*#Uz(q(!M&RmZI$DlLt&dbm1Fq`Jv9yWrDXJ3 zB;ti%w*jWOB^TM~XLT+Iyi^}C&`iG_KY>l&V` z1kC(Q)R8zH*%W4#!nWQfZ{9ihj9J%e05XD{k#p2K&l;K@k23~?Yb@JHX7O-;G=$Td zYHgqTeVo_vA~s@roq=I~dCgY~PoDn@BL!Z^9!5^vLRGr`D|MNs$fYHtY~uG#7x8W! zPv}r&?Ln1o?xfy416?x1Vlgo1ed|g0s`T>uczHf8A&kB^(OVw7%`L)Q zN#TB}*m;U9fQavbs&9)_DC->eZ`rKdb@6o4%KUk8JX^w%IBr;`b1yj>x4&$VG8fHt z%M>5nK^b>qYggpiwVE)zR<6*A3ei?q^&I*&vr-$a99?26?04S^7@a^p)~Ymgy))fYf{`Wmc>K z^<8otq%n4wpPtd(2P` zA5$ni9xo$+2qOr9x(_bXUG4MmVrY%`OFIr4osHx_$XBrvE|{v|V1sxU&s%q}e1TtS za!HHT0T1{FEGF*d8Ovm{pTyM}u~|I*(6k$)`$}jI{VA+gKyXNip@Q;u-K|i3)DVH@YJ}T^ z-U2p{Cjk@ z)du$KAd573<@M2hj#w+kK$SPb4`f9{eRJ3xB{p#RX26p@v=;>J#AL1 z=;gMB*h@{NPH-batF4Ivghv~42f+PXM#xkIx+cQeScB)Ai|1LZVr8DslN$ea}to*G>Q)?@!2DM{{W|{ z?^7U?eHFnIp^ILa!uvWo!XOC|m}1rvGvtL&KV!>{=Rau&PV~@S>*;1!1T_{$ zCSi9ZrJtCIOyd@UN2gUs+z+V+2%8+AfLcXDr{RwI(;#>;Pf=v65TA3yml1#dKzHx( z%G2(F<62GqGL5*x#u_oI56@(VH{gR_`6(-x3~NAAoX}_g@t`K^vIkoUS$AA~z@7y# z(`IUcmIF-Y)j+-EQ`tnW4Og9~LTJPc;xl@r8%PM_DB7rsm*x3D%*pJ zWX5((P`2FrQcytk1uziBUhyP?p4sf{jM~BSZDh?MlI$=q__2D7iJ8t#XBL8Pl1)af zrvs?X?VbV-oj0Ci;k#=i;X%rU_>Dk{_!BO=;}F90PGp#yKcmCbdX2Jm%yiS0HT9LV z#;sRuTsxF5J(9t<Se6lc>UGoa8E^cb&Gcwn%ZhI1a*mU3R#-*{_ z9$T1C2)$R8s#XH_TL27lvz-wI+u?hgp53sjSoN5`PbzT1D66f>X{TB>TdJnkbG-<{ zT{mLhch=A|y4qTc`Xu0idF(lTd3a+0rB~wp#Padz=p{%TqF$;f-~L;HTUX}=`x(o_ zWf)KUTC(=5hFatKt%P-iSOc1SD%?!F22(?xaefA95Uz?oR#%@;M@a=I? zN!A}c`w1bO0r5ye^K@&d-yp5eYh)q4BU#3x$uydl!h~j7=wz0s2B-;NuJhnJ-nTxc zyjCtH1C}&0hbPUy>-YUo(-!OzOP2~&_q&j`eH*ndd611M)t7p3h=1I3x|?>tKzGF9dwhCMle!}OwbVum6HDsvXDVUXhk@VRkWHvsSa z;x{MxgwOKfpC4R;cx!=6Xsu=Ud)9w}v)t#MSo)7d!xkj*bQ~oVx~OZG3JhH&$HtE5 zakAfs@o&$aDN+;_D-(IL^8d~E6FCbQ4fO{Lg@8bpmxH^8!Sly-)Ax`iI|~n=rzGW- zF0|AOY-V?>m=`%&fjL5SHe-N+nNkMRjo$)Psuo5D335gve$>0ub*7 zmt$Kn4TmeEp04{ME_lz9l~HdqpN~6yYUj8YPeL{c<45UTtKok<4=gv=+QT~PB#hDa zUU#jH*(UIl@U$zp>q({bl-j~ysll=<7u`=_d*ERL`!LXMU#(eSpoR~Jty2fqh*r2E zb=B3E9!_1`K3-e`fvpc7PH)B+UHu2| zzX5BW2X6%|Jxj9}f)YALyD|N#!7h>Of<|;KqrgU`WX+QA*4-2F)yK+Yc!w21Q)%X3 z3D!hz)b~R4utvUAb16dX_$q9B9Z#we&(ywL7c!whXy}#7q*F?J9}^DkL^rye$u)VF zTKV0}4ST;Jyvvd}RAg$bh6Z>2k{RL%+vvTT6F6r0D#O(5Q{mg}$=!c; zh86qNxBi@@LUDY2*!f1-SOT9Ux`&An4bBx}cY|(Q>x_Bgw^%mqW_$tItSL-LQf%k6 z5d`0z%%M49*fiW}r!@3fLGh!KEs>TNE3|(n`bneod$IJ4nH{qRLlu?S#1Lw%XvBJz zOCN4!`8f;~jk~Ub6fO3}q@JvV`(y4zb_}s#xN-qvBbWxk9RwjRX^)CGJsUkxdpMw) z(LSizdSR~JvHXXJaZpo5#i+z9G(GUkM}esv`?U{k5SHoKdE8-~Mh%i#vSi&49?jK@ z8*$%IV*N$YfZ8NAI6*fawY4M?jP4g#Ys(QH^d`d5P;)gXF=|k4QmyQcw|`jsw8x?}&gx`xrEV;H z*#uX7@cC?p0hfw@D$aeR(57=WXRplZk{+zld$y6{bYXGtjnm@9?j;s)`%{z>$D38D z2&AqGVkZuY?lv>n_w=USvQVmoE+D120^Ct{sje!02{xllNO+Q-jfwQedD^!MT>uk0s=S(L!Do8|0&kw{`Lc{{c7V+$+th-Bs8xnslVJ- z%Hep#$p%B0dP9&b-Hef6-gI1vT11~tN=L;ne;S7ce}4Z^R83FET>#w_o1ECFZJ?(; zu&Zb}8aeiT#M}(w@-xI@zKc{bwSKKv&`6-hYxvWLJ7zCJp3|)?BGZS~a-VlaY)Dq0 z!j*FB(4?$1MY>WJ%}GjmV-2VMWf4C6yN7vZDWa46A&RU|H3kPT?Yb^Y(DLM9eIJJf zipJ017D-vKPNzjJ8+-@42tU#%_u`~cJ9*95-d4nns2F%hXy-&Rfa&9_NZSe__=d?* z7TbEcu#-p0yjAG$f!>mQ;MnD@Y3G5lQn&FRe%bVHjWc(Q7ju8Djl;<9tL@P6mwFpq z@te4~MVi!STWfvef-VTN3KY{Ybm4O_!khI0cg9gjc4I_#qe|mj{~|rF<)Dbse2+C@ zn8w)JpOt{>V1(y%XJ+uKni~tB8%xoJ^y_lUsfTj~p%s=rlb&|9n>F&*R3cfv$^QOU z?{uBtP3DUiL3eY-orZ8n$4p7dS<}v$l<#^Jk!^R|B=>vwQm|4zW|AlM;A0dFVvIWU z3>uG{YatpHwffa}idy;Zt_(!zqmUzNM6V3gJ(qGq9+xOiLrXNBQ~_L+c%Z|r@<_eOkNNP!W>EiLN!|w0^R~ z12$>Qm(XAxd&D&OkxtlYaHG!d6TOK0>Zic2V-2j!8fai?MU??=df5|l`pRB{z1r?i ze`j1@9dhZ~P&;C$WpMwr(~TGUf~dr0DoBzM+9%;W%-&Ii>tH2~onEKCdAJ+YXB(C} zp`{5O6~4i)=Y#v0y(qG~v0w!CI#qLhu07B>Inx09N%&NyX11&UYXfPj&wVe${oTOO zzKQ&L|Dwf)vDAI!NZIz&t-@3hB;v$uZ|5!Wg)Z+wHa0d7#7_BAy`7Bv{6L8!B(exU zeup@URIk&dNWpvm$OptwQwdYNKPYTBsKXGm2Dp7owFpCzm(0W8jC&VO9IJ{SUUg0x zaLUrHJ5*-<)g_;*GUfe<6EAlaneek1K9p0@E_Ec!dD)pn3AVd%oczYS#Y^U57uUL% zm&6iJ&D(>)s?0V$udLDO9a9Xw<|h~89rUg~cbP|JvZk`cF4JMX=ZCr47d+NO`i4dbMjYW`(e-#oJcmgu%;&jrd8EK$x{P-DtO$ z%7@+F`2;UacKlhNm6a8d>akf@rI~n_0C*cJ%MYx=57S)5EMtYud zq&IJmmqHr7Ju3zbAI^6Xiujxy{40-aSeg}=?@LjQQakJ%DiV%ZJG3GW4Vx(U*d&Z- z%uu&Vpa^Q)_x3ticSWJtBTVZ(mgm7Ap;XQesfyD&C7lNCM)$?t&mG0!v!tdZm%h3N zGU=|Cd~cD_K2;URUUiNt{_}$k*@!=)YPB|1`VAK~-Y@E2%?4CJ`$PGEHpRYXyv_}v zT#RvnVbJ8EAo^>V&3r9}e!WsdeTeQRvX zr+$R3DLz;g=ed|{nQ8fFxWuR^VOqd`2A$kf;8r#~g@F_u=@Y%RvUA&x(8K>LsoYxa zg4lYYMuRD1U{5Z)fgZXoH&&@cKM!+qR5b>xc~0c8cuM-@dy7&#&#tyIQ@5AdPN`s? zG3zfqKFO__1uZ@866dxSK^>%NJ67uCR}%bnSXXSem`XR_$Zg~!Gy8*G_oa2WO(v%B z4bke|$mPsA&|F_iDm>6)dL83cUy{Gno(kpT;JqR9dPJmG+qodGs&2TNxB0#Ci^p55 zuqdwFwI>)fF(cgo_Osvti)cn))7N|HI|Z^T^!qlva=p*Rj)#J+0w{-UYibwA3)V#{ zB|z`)oUS4k{$Mw=G%nRaPIv?C_jb28>Fq~zzUX0I;94ua4tBdHKizO-Jc>c{2r-p( zSiNc74J|_{uL|Ds=wJL9NNLCIQCgVWGqTTCb2$%-vD7gWe{$9;-`4T_)+r|^NUE3RY|teX^;;PVd7}QEGuL18*l_dzo!6#OjSK z;7A$^cqQz4j5g6Ly}X#qko z>otPE(0)n6ek(o*)uc+kTJ7;Z$;EQ&<0y;n_>1fT+apUD#?5tyCIu9hrOOas5P36k2F z@k2xNjWb8}w7BYZNc2_s&q8L#O1xtSEr0)biX$*U zO0hPdw(LbDPkpTHk`~e;x+ZPAur%SEvtjXg7RDJ0(q5^m?-r>fj=X*yFZH9S&_?*B zraB0y>soO1ln*Xb&fHYAbm@p!&d4BUo%F8Pg-KX)<#O+8YfE_vJQUC+ibSE`vdU%; zESVYyx5_D8=1E_Svy60WDsB5+cUvf{i1POOc4SENDOlW{H*!|n)dn-#fWjp>n`ktd znpR~}h$Q6bopSRx(bdqtqN_Y>IdQrOfYmDTg*79gdJ)VuC4-t2557^W2@C zneA`nC$tl$>aPXwAjBP$I?b(1Ep`|`U4cC=P_3fD#hq;293r&UXRe@DPu|zom7EHs zJP>wiA(Lg_Fe}E4DF1RWGD8)OwWcO%HjZXKLFHGK~ffBCeTHdRZtvZ&jP#x)QGJPMV+$ zE}xi?aGgD@VQZi+GYAl@Gph3Vq*T7pft(69>5a(e#Na5?Mq4`p22T`m*e8~LpWPHq z?4HHo65+If8;1j-Qscn{YK)4@n)*gzGS3UXv>XNJ2`hfdA%Xk#V~bANm8nUU+m`RU zy>%wCEY0&dil}7vKa$yT32v9B+1R)Ys_>$m7KB22>*tvM1sC37d>J-S^ZeDqSjh6rfn1 z06M2dYMaZZQ}4ZXvbEPTz~X4o!fwrOdYU*M8L>33G2Xn7*EUU=q4=%N1k1?zZ2cnM zJTfA3ebgqQ-&H*qFoC45WsH>q*)>e_0Aqg3qASO-1#bdpi%-J z&Vs(frEc!S|05=4of#{#g2+1ltPimfsge*+S}RVVHXZEUqi zt@P2NP9@4!^$w`C{fF=8U72^E^F6B8lhJbRn4Ixs7QtxiZR#+WF`$2;p`qFDLo^?B zX&tSv)z~#K&av9pIITCiDsCj?>DDO_*n;X^j3_v?=_%JK{J2r z3QtV+YLa`qw|s0M{ob@SK3#0JBI#k&&Rl__GlNahk8KD}%d?LQw+fM%?tu@r{^YUj z@D&GWe+g-T-nel*mQ^4QJyBYFhl_G}a8RVA6G%zXh&f4Tfk_y?{x`8IJF#l{TX!%I zX6f`jdMYW#x%~T&*J}c}75s3mapuKMZTsryh0;-ZBJ*%DXTbM*=I68h51R_opG(ar zPvoWCIY51}PpW#LE0uzMBJ^R5(r+-$^B@W0ab=;Q1I7ogudj8oE#B8o#zZCP-2cW* zJ**1{+r?Jdc@OT6e*(e^g{?oB?IZ<;4l~1$+i3Dw)uf+scftXZ_l+xJU8O>-o{0~f zattsoy|B|9uztcDJ(l9LKdLP97#wbQGfLT5S?A$LDj&>Bq;%y8mRdyxwwhZG610ux z$DJEHR{=~UAP~z1_|*SKky4V7$9_^y$xlMhV+;d%Q0d8B29qbm8>OV*Mw~+m&Q9-v zPRBm9IT;b;1Gx?E*T>~$dwRRrH18yx>R>19O}%Yo)pah3uT~hdF~O{U<aGzJ=tB+VEy&)3jcC>>sgo^*UvL83VtUb$dkc~qqD1+ z^GqF`eg8lI_T7W5X3@XBbSFc#Nm{plvhYpq1JhrRa&FmUty?x>nF8SUM-apGVuu}j z>6EBSXm_f1g&ba4a~rkrlm++E(^9jj>uIR#X>(8${@Q+H6x_?qD4LHuksZ5oK%T=< zWV+NeXB@4Yo`Ek4U)p-T9J*@HweF~PoJ2>~LZg?(%11OlicOd(PE4E-<3pN#b%Xl4 z$OF~vkr?Bp4eVQ;V}l&kc>HF^^-v$cfDX>M8_aV`uwWtT|4oYTX&s?MFPpg!*5C%V!vDcxp3`iy`d|dX0jKzK|I=0_Y8de#b-?h)GbtZyx82Buk5Pjpp6T58#kg%-idR_q!w$k0 zk1|zh33~C?7UY>hMXQze6YwvX~zQUr}!KE>uxzd^~?Wx^gQUBbRq(1K?{U7TMCw>)8CO zd!FzZeRUUQ-zo7Gb>ze{@`qE-HFirFlD+wyN9DLjQTZ3tMT%`Q~Q?K z{gef4W=BG&y|R*kBr-LtXNn+R?EGkBv`F`-m=1Sh}Pf)yLivP`5`ue!@r$fjiEDnq$G(=VT+A8P-Z9mWL^ zqjIv1^599lwa@j36sEme%}V!c)22(VbvNSv?3VzFNvE|+mls!Ytfi)==3>z*EpH-F z1Ka!%&A!eRO6K8T9LY!?^yk8*i~o9%vL<|L3!3l{C-n4`Cu9v6F!ws$V#jJ>)b>{1 zqj~DR$FIlh^2P^LI5II=&CjhYEYQ)-z`5=|=*Kw|+giUY^gP&`ywje^yG`NL@XAST zTD^;iioD3yHvu}y)t{oHB7 zB844PVevK*q6?rzgQ8%gQok;xg$BdF zLsID_J-c?kl*`U9G45Y+;1xW(-6uyv+gYHhz_C9m<46!cdJdbvZER%*rfdq8?D7to zA8_uBD5b_N(oRu|WAoLNFX-F;OgFmILfG2wadh0;%0|R$na5SCdKKo1aS}FB=Vi9# zJ1qxj`YT{O>g0ZRDRk4w022;1UcnVe3AmGzli1(`?PSXDBX{!?Ti_$Th$mZyX6wtq z4t8|h-$8}!9-vpLi;-{y8dod7cCp|AE|~I$Eoa@sn%6U`La5wTnGiQxKEcx_)_^~| zBNT+VKh!jg)PKaQMP27_pu24r>?F5`DUi^&3CS`TaT9hb&UUo|L*&%GK>N zy}Yq216rKEc({b0!BvWo(TE*L&9iF6+Tb%l7Ib zIj2vb?%jK@wR-R55?~y=4Hor;LEkrJ?r^#9sdN<62ps$rfyOvbTKqGV{?!K|=sw7z zy#{h&2PkR?n)j}pPdF|0Ox{?hqKu1+d;V7uAKSn|+oz$=v#&C&^2iy7R;TvFhW3^H z`vjtf$?<;v-kuK}HJ-@urv;A2ca|(BlpQ_pBJ;|<)=r~p2Z+e1ANn!~d^0~TF+n1p zR|w4Nne7TMC&2Qv>DXp{0nnUU9DdXmW%BV>Tgl9ca^*jdrJ01q>J56~_yBL@bRN}y zjNN~Gmp%`D<`wCipYfu{Tl){_jvl4oZXyn(wdo>dyQ8;qAz$0Wmg?$j>K*hgbK|p%qloPR7+Lt zn71X=-4ro+4W>&wQ;4fahje1Bh0DD{7BI4ff}q8ZWhe#TkNLLG6rJXVXFQT)1^EIU zR5&g1pC#HZO&fX^ULJfiA*|PZTTs-xpwGypv~Pncj2DO?j;nT zY{3k%&*nutd!_hJ${or{Kckl-r4ba~v8Um}9@BNx`e z|Lx*rXp+VYZq;8fE^Ql(=U7@+y`x7+s1_pi)YlrdZ;v{-E6Q-;rIZpC7jB6Kk)F;d zT2*?3#cJAEK*sGn(b}6O1fU*eW|e}8UnMYDn|*t;1gq~`?$V>Yp^~16kdTn0w$iJ0 zSbCdFL}WwKR*$?6Cl^6%6`Gk*;anbo3JWE>%Wt-j77g`Et15iOzN`QA;p~k+cvO3+ zg7ILv1x^|c^tyxY7AN~S4lE@GmeqMai9X@8S+T^m{J&A3N6U;TrAR~>EJ#Fg9JCDvhYb7c z;MAe&aHQBmNDWK=MT%sjEdPV?v9Y-CW(f8q)LMsoZvCq*UCIlUFqpOBOpwk%=@y zJ2f2z6NM{Yz`xi$Kho~(6;br@wt`g^&Yd_3d8}`=4Cf!V4uCb7+-y~+%4oZr#xn%G z#nn)(b=@mVtX0>}0Mz|GJSwqj7YbwPb%2qRqs~en0nh^tK)8S)S(?PLG=A8M{X8@qhaD^7nM5hBJ;1+eL(sl(j2{(_*CYhb_w7}|9m z2&8u&-N5j8o70+SYSop86;YW<6$e7JuNwYzXe3>Q^8g!gD+p=_R+inQ5HIKBlP@t` zHAg8jbmq)nLyW)!RZ~*eQUx^5<0AQ8!ij|*lilz-KmE-m>!z+2J%*G$$f5`rRg;%l5uTFQ_~6a3R->Y{;whxnae1n0Fl|AHPaMwo$L zht(|0kUPE|!I3CozNm*hJrvqC$vBd%n<!v7~HwNJ?>mH@>@AjGoCaO=7Iq zEl3H=bcM8stJ39PV9apoRqu4lG5AUM?PAZ8>*#z~R! z&jA&No}=FRb!Zev}u4CaPbFpU0Z$dG4V4 z3h^AWSpI`tvrV-e^wPCt7o*dD)u8!I8OkK9_HmJ~t z{~SFp7pH*{u~;Omy-9m2<`d8Tf>%lT6r#zm86pL;Wy(``+!ZZ%aV?3%TcWZshi3MJ zcZ2pS>ST>j=$9cYF}smm3I@JR^>%x)T!pKPOyxXzh*V-JnqpfL#v*ERfCMMYn*}e|@rLT<#w)>0zun ztsQ)f(H*+S+u>0r$;U{3(ZE`-Njk7nF-Y6u8rqJEL4WU%bv7~U-W4^ZmJ-Z-+Z{~0%cAPMC}Z0)M-*LmFG-;L-546W@{4dnI;RJNI01KEl-s}y?(0X zY`U$txuDoAqfq?s6Kc%#|e6M8G zH9Dx8gEWbB&ocyl95=72@qc~(iY4&;i$|Li{9mvQ>_sOuy>cyeQrAf-OXw`)qiYMY ze{83ZN5>K5yY7(C{)u7E6s)vLPxJwIZRmD9Iu1Jt3nUfIJ?>#%jJwXV`VgV~_g)gs z=;L7qSmzUyt&!Vd1)h_G{n~?BoqXEz>`Tt8!Zc$lzz^evF+g;Cqc&r=yr^7GTLoluP=+1tzK>?weoD;qlmV=j8j>? zGeK3(KYQDoY})6uH}cqopLwCzO};j^$k)KHKt!?2GL7h5Mx#q9j>cPjn~z;7{2Bn_3{(f&CnLF z6Jv**P^sxYcN_nORqP%~0qK8Lp)9_*DVj40oYq#7$-JaH0!Dzx8&h+B*>WSV_w zYi*lgI~zH>AWw)(_*-FQ&Nu-bNY%Q-^XF$oUM%kNKKa-^yE4D6a(n;}TA{ch18Tmk z2V4dMa-1Lv2=h2y>8mApS~1i7qT^d}jIaT^8|xo7mGMNQWAmlJ1CurJD$pDGi8G*L z@W#5!1w{|uec``|HU&CEru3an-O^gymel5Z!eziyvmXu)1Jv2r*yzpJ-qk&<(M(lV z4_c$e8UWZVm@iW7H&w&Zo(c&62SnNMe?Tz-Txm3&+v@A7$OtJ<`8PtjWW4aS+GQa3 z8NE``KnN3z1p>Kx z>tLF$ObL;@q`X8OuH-*z!IO>Q*=x;z9E|k33dxOO>g3OcLXkQj8|lK{AZ)(YzS0Z5 zdSiNiX^mp&g{Qlyh+~q`>T>hj9*A1mS5J=_I=o=XKp0DY;oAO{cg#U(ii$a{i>*g- z2%nr%8ifzxp&vU9nvFY#kgBP971Yn(!F{UzzJZA|ZSjfaq}osR@sIEQ>OqY<_Jf%D;X*u zeubsJzjJIv(>=9$yQC+JJA&Sb9*WepV#9ud_a&Zv;M0CtpB{sJd&W7owcO>Uh$lfJ z()4?5*B*qu&95|^zpH_Pd2R}o6*U~JAQqDkDJ-&*yn~w zUSi+x>K4ZK+N}B%8BSzVXzjv&0O6`l82ZtQF?5q=b3$DtWuUwzotHv)aJ}r z;o3==gU8gvm8m_6*RAR3?GuIAzvl;OSX@URMBk#-suoNGf`Zdvxjexd$8maHPo$xV z+paX>#L2>=jhZAyAFz?2HwVa6MoW(SxQ9GY(qN$I8g9nxx=z)UQGb?MqURi7aSh>G z57$){-f14&7%%v{)GZ5{n^OG?O^7tWNEJv>Oz! zanp@^hsLX|Z&#=%+F#9is~RfI z2%onc>VH#?LLb0qI$BsPsfOz2hJmW+80h$CFE`8i)I6-zV9v@K-Hp@ zJ6)YW76$|?gYrhjv894XYhJ?lj2Dl@{ShYzsaTMGDY&rE-{%Di!kt?w!~BS_eST%<*AD>z!Y$0{^cVuiLQg_(X*q52{V^k&CyDJ; zHeC@yVg^w6^CYaE2}VvFz}c3-_xK+Q!PP91kEEfFneTk>%N3k^&RF(mg6n1Iz(|VR z5Uzh7PL}AT3Lc&1|NNL`xrdaZu28GAs4rqfp4f96xk$qXBxq-(5_M{5npT*9GV$e zCtGWY+*>_Jn7G=1JC#VIVl9cAW|?F>^XwM{&VTtt5pW$|9Okg+#p#ca_q77n_{83nKp<)dKC`u9n?x0u`eAg=Lk;v6WU{^CH0ms}MNd zKaFs+Z>@<)9WUKf$sQ2Ji4C7!yHOBMejQRbDl?wselX|uyrR?IcY|Uf=ZIenr_XxH z8S(iPS+=1A+Mtr0&&IFl^E3l{R259tw>E2Y7iDKVBmQS5{HzRFxBQ6~O*YDP{txgL>`a zxDm}6d|p0XikPfV8qEt{{On|9LT))te$4%T`W^J``AiN?&@E{{qIneuI~t9tg)V=3 zi1)1!5+rGN6kJ5{c1g>l2%1Q~KQEYGGTA1{eR=nQ#d8vkR7SXL4dw}^nAb?)9;ZKp8h zxxatECeVTJq!i)pX7i*4ThT%GSbm~xzC<%&l>LkJ{kV%LqMu?EoVtH;aVgRs4(VlE z<|)q!kSrfQ%Mio8d%7&{)Px>dn2)y#W}U$tza#m-`MTyMkZeec13aQI1NpHOn7MwA^DBs?13kNlmtz=&SvF&}M>`DPFEQcA~oj&hR}-7mT? zdjfUBp~)bY-%R^zLT|EhVF?9PsmXu-CQ&CfXQI47G4L*}1=ZR$way(T})<+W*63d#rih>q}$v`ldDjt#T%6roys9{`RP z;wAs@*YMIQUCCTgXu^8oxwWtlK8AEhw>M?nTnC%g_p7wwzaC&o^}p}9|5-Qx|9AE! zlU`1=M!7z_5%Z0RVzRMOmNvShWQ>!|J;83S94L_+8fG7|)mmlC7kU}V^__RQ?S>|1 zF&0Dtd(7i?z$%9b516RZa}3ug-ub9z-kzN%ff-+qM|{WKq>Zd;$@ctr3^Jb&I&?Gn zd~Kr?oc#G3e}~|pO`0Dzyp$$!?VgghV(i`t@(KNp8~ztE9#)%fs$Gdra?-oEzsEZ& zrYcK^08{tEM>X?bE^j=$&C)k*$rXu6^vVSIVK;`FN+MA3wwty6lzG23^cswFp4Y|C z(W+Z#W5e%$W}5%idyeq)Zbm%mrObgrK`Gbv3Y-M<{Pq+a>p)~1+x14_O%1Xbjk|Skc?M%_1 z&tm`j)wVUrZS9@k)YUSZUl!fIe~cn6m{N~MH_VzE<1Y37{t&T}%1_uNPm##3|6%&G zyBai1zqmPpCxMm|Ez5#flcQ(L%V@B_@P;f^EW3m9Na4B1alm4FVm4!51UW{{LG%f2 zc~Z>(G7KKYR_He#-=gKv@URz^IZiSRQG9xNr3Z70=vk&k*+jn+EH@V~bg+ktD?U6v zz-(=T5UXO8p$gPzqPI9gjM7QY&xQD^pFQ})h8mINLKhT(KqN6XNjoCDL_<>{v-q|; zmY+VYAinH#s8`Rx`i*pqps_AJ6w{j3=M0?f35o^|2B8`CEq&??7TH1VtITGXvd6k7 zA(*pr`Y|Gl7YB4X*Csxly0^=bto(L`iJ)(1z5HF)HspIdI84d0ws#1vw&?vRTbcgp z@{gD>0euL0M^d4TIRO87Z(!U-%zyuH%<%sc`_M!;fzV_E%382jKL6LP{K8+~)S{K8 zus-(Rbu10HinZh*%a)GidilD`-MaoKchvRQ{s0+P1453bo)7?#Om`55mcKVNvcdAw zvtf*J9#4qZEw6Qwi9vi&9^^JtA&1?bB)c@5B_?x`#LY;KN3L&V92NcH3QtIGysAt1 z*4T_!ZLx6}$l_*2i;#G^^y$kjcyvH z=c-9elE`@1-24@5+p)ktu~|C2kf1O{&g|Tx>(_mW+H*gFSwM4WE)1 zhVN1le`B}#pzlcC@#f>MF$Y$|pKKV~I-p?TcGdoE@G(CR5cb6Qk_B7z20JnSB9=l&wnJA|BDhq2B+Sj(v z*K3qG=pC)1Jj!L(IBU=_^N=3NmsxRzhMn)clyoClD)7F5C+I32@!|7G91I!IG8;xbNJy4GRb<%uv$->{k8}e*lkD}d3+MrMnq}qm1JO_A4vU zwpCr5@r0mH&(B{q9+U0+{}qePr13XR?GQ?$xP3dLpJLC|M+0n~+Dr1(xk(&An%KG- zS2=w-it2bjY}DK!_;R0}i`LSqNCE~{zJ%)awkeyNK&iDwNdK}$a;Tp*$x7{M4XyB_ z#Z&n5i-OG*PX9j#>)d=oOhQtFSxh}_AyQwe9haZ_(BsJ>N-s|H3GF>H0`!NB7qIqk zmD!XO7yl*^jah-F63rIG5pg|uxN02-|7Cvm>W?Ztx}UAPW!#dqyaqWOs&^KhhLhfO zM*!5pcQMw&>Jho2W4`w$!J$PHduQBeW- zXYmBn%0%#m{P$u)j@o6 zqHy;0=cR(^X02?FbYN5XvExjy{C%X~zLjQqgKwv$>jOK7;T(l|B#_PVR`~t+QIDOf zKo@Op?@IqV;P7+f_cEHPVQ;-ALoVA?jZD+GDp6|lkrC-8NpIw2<$DAabQEA^->?pk zoDW^^87)d`FI?5X91sb@V0C1UMo^NEzgAM}nPs(+n>eik`*@yll!j%LCeiC$y=S6v zCU_QPlObeYYxPRev~*$Tq$SA zEXF`^ET5`OGwPBxAor?ePd0{q4cIs;b?HFqpbP~E!^Izpl^LyB*8Rk6<{de1&K2@H zJjD-~qs3ES-nftNf2`D0FZNSg%N2(f@U8`~fv1_Zs|Dc*UcbAX9Tr_zK3yKu^5@GW_N-fUX%$rW}p{^Cl%n0(Y(mtJT-& z+jiF=Ea~BJ0ki#=g2%+Qs@Gj)R2&R^&spEs&2PiJsNbJMsOI=moKzJQBqb|CLMG}F zh@qiUNu1wu!tYmFPje^62rfcHeWJ~EjKo*^I`=w>=hGu?DJlXo-NaJcyGLyf2fo4h zkCgDVWG=&xUIpwtxp>{A^jr!K9>;+%eS$Y$?=RDAAo&>p{;Mcv1AFnoVc*onuVwef zm-~u3D)hb|J}93Yx?uaF1WM+ius+CGDl9je(PCwl+*VpDjIyh0U&ztAmxn`gB_AlXMb6pHsIzYk}Cf zrw;dT%3`Pf&`IZ0`Pt2?XhNz&Dd=F)0+_|~yVIUmgb;rBc|kzl%lMBv@|}N3;`-jj z#NSG2mJGP9tZais=9zO_e3WYtd2W19?lf^r2;SDN>G8!Q+YsZy8Is>5oa^|&b;$c? zl)v*gvDkbcQj}ZLfJBa~>NERp9qj8h634XFkFP6(RBw_Mc7&PKf(b^V9E^9xn(3*> z69>NKfodh?9G7}eZ=)=!oZpcSd=%dzb(dzBH)dvXm5ZlucBZ$p4XX^UCa%jbkY-SK zZSC(YP5<=){<8P3*rnR^p21tQh>G;)xaJwU^}dqaEgivqiD{~-NcdR_=uxsr^7+pX zjf_A3NQ1x4=;&axR~p&hI>5u-jD`C!0xX9YvyT7*xrKJ?p6`TkhVZUn)6kpf8z`wL z>-5NVvB!&~DjDBXPEUA)Ow=+)4m1v(&((9J)-#~J)P&~7rx;g*n2~a=p@4cZoj8Dd zO_`e^edZ^M*u0#`tJ}s zh5vsF5M0Vs(az<+<5+K%UpLE8QD}!0?byHnEBhkOH;Xkd0{1E@j@=z^qK3ha_0VV< zGCl`GF-tKQPndF5zk++JG#TISP9W_Wo~wtpcnCl)_4L)AYV*7wTMB<&GmrwM=_a1L9VXZ)Y%2^JsSF_No-o@SFZGJ_Q45|u^xseRLUM|D zhz&4xiCfe$*7QBOYrH_dcP!+SdhMZ`N`&$f+HC~s7#Ta;sjBtNW_njHC1tU#Q8s^%?fLt=rnA`CX=>CtEo8pQO(p$2Iqs zX8|wBggtwO;y3AmrxI5irxepQEgoEcO;pEwAIgLxg7h4zY&78AVZFN@gKF(^^*KLL z6Ee@uj|Mvhyn{Q3>RD&}CG9fh*|P_9Ys-%frJ@FKQtWuKoGS>z?ufEfS%o&JD#S(Z zw#4bi#MV8L)`zZgnOgU}tR<`UPP$zkXV}nPk5H%Wk?WqQJt{y;wn+;iK~Q^fTTQ>k zrSnv5RP-{Qdca@Y7q({YIqBF}AbhX%p%f0Y5fAuVCNEAR(W5=T2hr4KZhM*s_IUm- zWA5raGCQK@G-9D{5SKW^PkEwwg`N+;EBwnhmc z28%(55ObY(kwA35B4s%<9ksN=f#05LB+l#iBy^0f?18QuzqqGys7_8a2v&J$Fsu&h zm148eu$*K9&)U+4VNx^}gcyIi7lyEHhs(F%J|Diz~#Hk zjJHur@Sh@lcHt|RZr^JsYAa7_WA9pBOCPCLLsziDgjgKCw5+S6=?P599>YLW=6@d4tRLe2HJ%*4EZmlvm1a z#HwA7nO4bHyN#|HqJwyEu8Je)iS6b`a#n&Xme^U@SI;x&}Xs&admu)_= zwV+1-{RP++O4_ZSQ#+YJh1-!WsFk3l<)$2|nqO{kf3@bLbhrnbuHS3I*hEp|oVVrt z+%!(h(`z#`lv9iA#$QP@6UQp`UylzTlAtNUbOu!oVjFFj^ub5!k1v7Pl>G$U}sjh~eMhKk1{0Nn)Z<$$!0 zz0?6ZHOo()=+laT?aU%Su2i(hsYY+HKZHQwd>As*X|G4jD0Ihn+J~5q;Av~S!~Q|} z%D4DC8KD9;JYspu(%RRSUS)C%kGe1z@99zaiE?l*<4c6dEzZpm3G08@3H2e*VDJw8 zaoBwE+sR2-Hqi#s;>PfSitTAs^d|}u2lkm)?`60QqMoYHd_Utq5k9zasK$bDF>yyy zz#bG^bU)L-+afI@ean(?606ML;C+JasXJ@3K<(mj!e-Mhq)x=F2jrj}zr`PUP6+?r zqG;0XrB%f}D|{V7;CH538$2Gbq-4V8X?N6_-lm*!GB?sXr$%aGb^6m|cL*yHR>ouPrg=M1@C%1u`?CO0Ymgd9cD>N8tjB?4%1r;S86<(4k*h= z^@n&niN1plkPfRV)&CiRxu(Xxd90&scFkw>M?cM(>zU2~N0s_|1tPR2iQY>vk+sv( z!h$Z`GS=RS8OL!Yqyy;*vSAgVb2aL;>AcYL+B)o8Z|HSCN70z1mR^Hw1{ESkbl%;k z>~S)Jqtw8aX+jCuI)xDL)cKw7A+n4sigsIU539y-QCuCGx6&h*8?<11(Cj{`WHPc? zvc!vJeY>*T4tRCa&J!TFBjkPl=I&>`2ZA$DvC3}_M=Dfaj2a`DKkJk9FN7~6r0(|xY-`>m}7cbJ$AMBPL z;J?(WKocV1`h$2*(Ze~r($l+Uj|R65x07!aH`Uizm_U05qDiSI_iZ^LHkf>H%iTI<;`Z~t(@VF*Ou448~>?* y^+0v(=VdvhI{nLhP!J+!Od}zXU;l$IF0Pn9D;>9UOgTZj8j|92ukuCpKK(EAQta3O literal 0 HcmV?d00001 diff --git a/cookbooks/dynamic_subduction_with_two_phases/uncoupled-fixed-slab.prm b/cookbooks/dynamic_subduction_with_two_phases/uncoupled-fixed-slab.prm new file mode 100644 index 00000000000..173ef6d45c9 --- /dev/null +++ b/cookbooks/dynamic_subduction_with_two_phases/uncoupled-fixed-slab.prm @@ -0,0 +1,252 @@ +# This is a cookbook which simulates a simple 2D subduction system with +# a hydrated sediment layer, MORB layer, gabbro layer, and peridotite +# layer. The layers dehydrate as they are subducted, with the solid-fluid +# reactions governed by an approximation published by Tian et al., 2018. +set World builder file = 2D_tian_slab.wb +set Adiabatic surface temperature = 1600 +set Nonlinear solver scheme = iterated Advection and Stokes +set Output directory = tian-slab-uncoupled-no-smoothin +set Max nonlinear iterations = 50 +set Nonlinear solver tolerance = 1e-5 +set Dimension = 2 +set End time = 1e6 +set Surface pressure = 0 +# If this is larger than the melt generation then that means that the fluid will instantly dehydrate +# Perhaps this is unstable, especially if the solver tolerance is 1? Make this 10x less than the +# reaction time step in the volatile material model. +set Maximum time step = 1e3 +set Use operator splitting = true +set Resume computation = auto + +subsection Solver parameters + subsection Stokes solver parameters + set GMRES solver restart length = 150 + set Number of cheap Stokes solver steps = 2000 + set Use full A block as preconditioner = true + set Linear solver tolerance = 1e-8 + set Maximum number of expensive Stokes solver steps = 100 + end + subsection Operator splitting parameters + set Reaction time step = 1e2 + set Reaction time steps per advection step = 10 + end +end + +subsection Discretization + subsection Stabilization parameters + set beta = 0.052 # 0.5 + set cR = 0.11 # 1.0 + end +end + +subsection Checkpointing + set Steps between checkpoint = 100 +end + +subsection Compositional fields + set Number of fields = 7 + set Names of fields = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone1 + set Compositional field methods = darcy field, field, field, field, field, field, field +end + +subsection Temperature field + set Temperature method = field +end + +subsection Initial composition model + set List of model names = world builder + subsection World builder + set List of relevant compositions = porosity, bound_fluid, peridotite, gabbro, MORB, sediment, weak_zone1 + end +end + +subsection Boundary composition model + set Fixed composition boundary indicators = left, left lithosphere + set List of model names = initial composition +end + +subsection Geometry model + set Model name = box with lithosphere boundary indicators + subsection Box with lithosphere boundary indicators + set Lithospheric thickness = 100e3 + set X extent = 800e3 + set Y extent = 400e3 + + set X repetitions = 8 + set Y repetitions = 3 + set Y repetitions lithosphere = 1 + + set Box origin X coordinate = 100e3 + set Box origin Y coordinate = 0 + end +end + + +subsection Gravity model + set Model name = function + subsection Function + set Variable names = x,y + set Function expression = 0; -9.81 + end +end + + +subsection Initial temperature model + set Model name = world builder +end + +subsection Boundary temperature model + set Fixed temperature boundary indicators = top, bottom, left lithosphere, right lithosphere + set List of model names = initial temperature +end + + +subsection Material model + + set Model name = reactive fluid transport + set Material averaging = harmonic average only viscosity + + subsection Reactive Fluid Transport Model + set Base model = visco plastic + set Reference fluid density = 1000 + set Shear to bulk viscosity ratio = 0.1 + set Reference fluid viscosity = 1 + set Reference permeability = 1e-6 + set Exponential fluid weakening factor = 30 + set Fluid compressibility = 0 + set Fluid reaction time scale for operator splitting = 5e4 + set Fluid-solid reaction scheme = tian approximation + set Maximum weight percent water in peridotite = 2 + set Maximum weight percent water in gabbro = 1 + set Maximum weight percent water in MORB = 2 + set Maximum weight percent water in sediment = 3 + end + + subsection Visco Plastic + set Viscosity averaging scheme = harmonic + set Viscous flow law = composite + set Prefactors for diffusion creep = background:4.5e-15, \ + porosity:4.5e-15, \ + bound_fluid:4.5e-15, \ + peridotite:4.5e-15, \ + gabbro:4.5e-15, \ + MORB:4.5e-15, \ + sediment:4.5e-15, \ + weak_zone1:4.5e-15 + + set Stress exponents for diffusion creep = 1.0 + + set Activation volumes for diffusion creep = background:8.2e-6, \ + porosity:8.2e-6, \ + bound_fluid:8.2e-6, \ + peridotite:8.2e-6, \ + gabbro:8.2e-6, \ + MORB:8.2e-6, \ + sediment:8.2e-6, \ + weak_zone1:8.2e-6 + + set Activation energies for diffusion creep = background:375e3, \ + porosity:375e3, \ + bound_fluid:375e3, \ + peridotite:375e3, \ + gabbro:375e3, \ + MORB:375e3, \ + sediment:375e3, \ + weak_zone1:375e3 + + set Prefactors for dislocation creep = background:7.4e-15, \ + porosity:7.4e-15, \ + bound_fluid:7.4e-15, \ + peridotite:7.4e-15, \ + gabbro:7.4e-15, \ + MORB:7.4e-15, \ + sediment:7.4e-15, \ + weak_zone1:1e-50 + + set Stress exponents for dislocation creep = 3.5 + + set Activation volumes for dislocation creep = background:14e-6, \ + porosity:14e-6, \ + bound_fluid:14e-6, \ + peridotite:14e-6, \ + gabbro:14e-6, \ + MORB:14e-6, \ + sediment:14e-6, \ + weak_zone1:14e-6 + + set Activation energies for dislocation creep = background:530e3, \ + porosity:530e3, \ + bound_fluid:530e3, \ + peridotite:530e3, \ + gabbro:530e3, \ + MORB:530e3, \ + sediment:530e3, \ + weak_zone1:530e3 + + set Angles of internal friction = 0 + + set Cohesions = background:1e50, \ + porosity:1e50, \ + bound_fluid:1e50, \ + peridotite:1e50, \ + gabbro:1e50, \ + MORB:1e50, \ + sediment:1e50, \ + weak_zone1:1e6 + + set Minimum viscosity = 5e19 + set Maximum viscosity = 5e23 + end +end + +subsection Mesh refinement + set Coarsening fraction = 0.05 + set Refinement fraction = 0.8 + set Initial adaptive refinement = 2 + set Initial global refinement = 3 + set Strategy = isosurfaces, composition threshold, minimum refinement function + set Time steps between mesh refinement = 2 + subsection Isosurfaces + set Isosurfaces = max, max, bound_fluid: 0.005|0.04 + end + + # minimum of 4 global refinements + subsection Minimum refinement function + set Function expression = 3 + end + + # refine where the porosity is bigger than 1e-6 + subsection Composition threshold + set Compositional field thresholds = 1e-6, 10.0, 10.0, 10.0, 10.0, 10.0, 10.0 + end +end + +subsection Boundary velocity model + set Tangential velocity boundary indicators = top, bottom, right lithosphere + set Prescribed velocity boundary indicators = left lithosphere: function + subsection Function + set Function expression = 0.05; 0.0 + end +end + +subsection Boundary traction model + set Prescribed traction boundary indicators = left:initial lithostatic pressure, right:initial lithostatic pressure + subsection Initial lithostatic pressure + set Representative point = 100e3, 0 + end +end + +subsection Melt settings + set Include melt transport = false +end + +subsection Postprocess + set List of postprocessors = visualization, composition statistics, velocity statistics + subsection Visualization + set List of output variables = density, viscosity, thermal expansivity, named additional outputs + set Output format = vtu + set Time between graphical output = 1e4 + set Interpolate output = true + set Number of grouped files = 0 + end +end diff --git a/doc/sphinx/user/cookbooks/geophysical-setups.md b/doc/sphinx/user/cookbooks/geophysical-setups.md index a7b5c2dfd77..5dee572994d 100644 --- a/doc/sphinx/user/cookbooks/geophysical-setups.md +++ b/doc/sphinx/user/cookbooks/geophysical-setups.md @@ -111,6 +111,7 @@ cookbooks/kinematically_driven_subduction_2d/doc/kinematically_driven_subduction cookbooks/allken_et_al_2012_rift_interaction/doc/allken.md cookbooks/inclusions/doc/inclusions.md cookbooks/subduction_initiation/doc/subduction_initiation.md +cookbooks/dynamic_subduction_with_two_phases/doc/dynamic_subduction_with_two_phases.md cookbooks/vankeken_subduction/doc/vankeken_subduction.md cookbooks/tomography_based_plate_motions/doc/tomography_based_plate_motions.md cookbooks/future/README.md