Commit 641436c1 authored by Norbert Grunwald's avatar Norbert Grunwald Committed by Dmitry Yu. Naumov
Browse files

<Tests/TH2M> add two thermo-hydro-mechanic tests

parent a0027be9
......@@ -27,6 +27,44 @@ AddTest(
result_TH2M_M_ts_2_t_2.000000.vtu result_TH2M_M_ts_2_t_2.000000.vtu displacement displacement 1e-8 1e-8
)
# TH2M THM point_heatsource benchmark
AddTest(
NAME TH2M_THM_point_heatsource
PATH TH2M/THM/sphere
RUNTIME 40
EXECUTABLE ogs
EXECUTABLE_ARGS point_heatsource.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
result_point_heatsource_ts_5_t_100000.000000.vtu result_point_heatsource_ts_5_t_100000.000000.vtu gas_pressure_interpolated gas_pressure_interpolated 4e-7 1e-8
result_point_heatsource_ts_5_t_100000.000000.vtu result_point_heatsource_ts_5_t_100000.000000.vtu temperature_interpolated temperature_interpolated 1e-8 1e-8
result_point_heatsource_ts_5_t_100000.000000.vtu result_point_heatsource_ts_5_t_100000.000000.vtu displacement displacement 1e-8 1e-8
)
# TH2M Thermohydromechanics in a slab
AddTest(
NAME TH2M_THM_THM_1d_dirichlet
PATH TH2M/THM/slab
RUNTIME 15
EXECUTABLE ogs
EXECUTABLE_ARGS THM_1d_dirichlet.prj
WRAPPER time
TESTER vtkdiff
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
result_1d_dirichlet_slab_ts_5_t_100000.000000.vtu result_1d_dirichlet_slab_ts_5_t_100000.000000.vtu gas_pressure_interpolated gas_pressure_interpolated 1e-8 2e-8
result_1d_dirichlet_slab_ts_5_t_100000.000000.vtu result_1d_dirichlet_slab_ts_5_t_100000.000000.vtu temperature_interpolated temperature_interpolated 1e-8 1e-8
result_1d_dirichlet_slab_ts_5_t_100000.000000.vtu result_1d_dirichlet_slab_ts_5_t_100000.000000.vtu displacement displacement 1e-8 1e-8
)
# TH2M Heatpipe w/ static gas phase in radial domain
AddTest(
NAME TH2M_TH_unsaturated_heatpipe_radial
......
<?xml version="1.0" encoding="ISO-8859-1"?>
<OpenGeoSysProject>
<meshes>
<mesh axially_symmetric="false">domain.vtu</mesh>
<mesh axially_symmetric="false">boundary_left.vtu</mesh>
<mesh axially_symmetric="false">boundary_right.vtu</mesh>
<mesh axially_symmetric="false">boundary_bottom.vtu</mesh>
</meshes>
<processes>
<process>
<name>TH2M</name>
<type>TH2M</type>
<integration_order>3</integration_order>
<constitutive_relation>
<type>LinearElasticIsotropic</type>
<youngs_modulus>E</youngs_modulus>
<poissons_ratio>nu</poissons_ratio>
</constitutive_relation>
<jacobian_assembler>
<type>CentralDifferences</type>
<component_magnitudes>
1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0
</component_magnitudes>
<relative_epsilons>
1e-03 1e-03 1e-03 1e-03
1e-03 1e-03 1e-03 1e-03
1e-06 1e-06 1e-06 1e-06
1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08
1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08 1e-08</relative_epsilons>
</jacobian_assembler>
<reference_temperature>T0</reference_temperature>
<process_variables>
<gas_pressure>gas_pressure</gas_pressure>
<capillary_pressure>capillary_pressure</capillary_pressure>
<temperature>temperature</temperature>
<displacement>displacement</displacement>
</process_variables>
<secondary_variables>
<secondary_variable internal_name="velocity_gas" output_name="velocity_gas"/>
<secondary_variable internal_name="velocity_liquid" output_name="velocity_liquid"/>
<secondary_variable internal_name="sigma" output_name="sigma"/>
<secondary_variable internal_name="epsilon" output_name="epsilon"/>
<secondary_variable internal_name="liquid_pressure" output_name="liquid_pressure"/>
<secondary_variable internal_name="liquid_density" output_name="liquid_density"/>
<secondary_variable internal_name="gas_density" output_name="gas_density"/>
<secondary_variable internal_name="porosity" output_name="porosity"/>
<secondary_variable internal_name="saturation" output_name="saturation"/>
</secondary_variables>
<specific_body_force>0 0</specific_body_force>
</process>
</processes>
<media>
<medium id="0">
<phases>
<phase>
<type>AqueousLiquid</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>0.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>density</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>thermal_expansivity</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1</value>
</property>
</properties>
</phase>
<phase>
<type>Gas</type>
<properties>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>4280.0</value>
</property>
<property>
<name>molar_mass</name>
<type>Constant</type>
<value>.028949</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>0.6</value>
</property>
<property>
<name>density</name>
<type>Linear</type>
<reference_value>999.1</reference_value>
<independent_variable>
<variable_name>temperature</variable_name>
<reference_condition>273.15</reference_condition>
<slope>-4.e-4</slope>
</independent_variable>
</property>
<property>
<name>viscosity</name>
<type>Constant</type>
<value>1.e-3</value>
</property>
</properties>
</phase>
<phase>
<type>Solid</type>
<properties>
<property>
<name>density</name>
<type>Constant</type>
<value>2290.0</value>
</property>
<property>
<name>thermal_conductivity</name>
<type>Constant</type>
<value>1.838</value>
</property>
<property>
<name>specific_heat_capacity</name>
<type>Constant</type>
<value>917.654</value>
</property>
<property>
<name>thermal_expansivity</name>
<type>Constant</type>
<value>1.5e-5</value>
</property>
</properties>
</phase>
</phases>
<properties>
<property>
<name>permeability</name>
<type>Constant</type>
<value>2.e-20</value>
</property>
<property>
<name>biot_coefficient</name>
<type>Constant</type>
<value>1.0</value>
</property>
<property>
<name>saturation</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>relative_permeability</name>
<type>Constant</type>
<value>0</value>
</property>
<property>
<name>relative_permeability_nonwetting_phase</name>
<type>Constant</type>
<value>1</value>
</property>
<property>
<name>porosity</name>
<type>Constant</type>
<value>0.16</value>
</property>
</properties>
</medium>
</media>
<time_loop>
<processes>
<process ref="TH2M">
<nonlinear_solver>basic_newton</nonlinear_solver>
<compensate_non_equilibrium_initial_residuum>false</compensate_non_equilibrium_initial_residuum>
<convergence_criterion>
<type>PerComponentDeltaX</type>
<norm_type>NORM2</norm_type>
<abstols>1e-3 1e-3 1e-5 1e-6 1.e-6</abstols>
</convergence_criterion>
<time_discretization>
<type>BackwardEuler</type>
</time_discretization>
<time_stepping>
<type>FixedTimeStepping</type>
<t_initial>0</t_initial>
<t_end>100000</t_end>
<timesteps>
<pair>
<repeat>5</repeat>
<delta_t>20000</delta_t>
</pair>
</timesteps>
</time_stepping>
</process>
</processes>
<output>
<type>VTK</type>
<prefix>result_1d_dirichlet_slab</prefix>
<timesteps>
<pair>
<repeat>1</repeat>
<each_steps>5</each_steps>
</pair>
</timesteps>
<variables>
<variable>gas_pressure</variable>
<variable>gas_pressure_interpolated</variable>
<variable>capillary_pressure</variable>
<variable>capillary_pressure_interpolated</variable>
<variable>temperature</variable>
<variable>temperature_interpolated</variable>
<variable>displacement</variable>
<variable>sigma</variable>
<variable>epsilon</variable>
<variable>velocity_gas</variable>
<variable>velocity_liquid</variable>
<variable>liquid_pressure</variable>
<variable>liquid_density</variable>
<variable>gas_density</variable>
<variable>porosity</variable>
<variable>saturation</variable>
</variables>
</output>
</time_loop>
<parameters>
<parameter>
<name>E</name>
<type>Constant</type>
<value>5.e9</value>
</parameter>
<parameter>
<name>nu</name>
<type>Constant</type>
<value>0.3</value>
</parameter>
<parameter>
<name>T0</name>
<type>Constant</type>
<value>273.15</value>
</parameter>
<parameter>
<name>displacement0</name>
<type>Constant</type>
<values>0 0</values>
</parameter>
<parameter>
<name>pressure_ic</name>
<type>Constant</type>
<values>0</values>
</parameter>
<parameter>
<name>dirichlet0</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>temperature_ic</name>
<type>Constant</type>
<value>273.15</value>
</parameter>
<parameter>
<name>pressure_bc_left</name>
<type>Constant</type>
<value>0</value>
</parameter>
<parameter>
<name>temperature_bc_left</name>
<type>Constant</type>
<value>320</value>
</parameter>
<parameter>
<name>temperature_bc_right</name>
<type>Constant</type>
<value>273.15</value>
</parameter>
</parameters>
<!-- ## Displacement equation #########################################-->
<process_variables>
<process_variable>
<name>displacement</name>
<components>2</components>
<order>2</order>
<initial_condition>displacement0</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>boundary_bottom</mesh>
<type>Dirichlet</type>
<component>1</component>
<parameter>dirichlet0</parameter>
</boundary_condition>
<boundary_condition>
<mesh>boundary_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>dirichlet0</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<!-- ## Gas pressure equation #############################################-->
<process_variable>
<name>gas_pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>pressure_ic</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>boundary_right</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>pressure_bc_left</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<!-- ## Capillary pressure equation #######################################-->
<process_variable>
<name>capillary_pressure</name>
<components>1</components>
<order>1</order>
<initial_condition>dirichlet0</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>domain</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>dirichlet0</parameter>
</boundary_condition>
</boundary_conditions>
</process_variable>
<!-- ## Temperature equation ##############################################-->
<process_variable>
<name>temperature</name>
<components>1</components>
<order>1</order>
<initial_condition>temperature_ic</initial_condition>
<boundary_conditions>
<boundary_condition>
<mesh>boundary_left</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>temperature_bc_left</parameter>
</boundary_condition>
<boundary_condition>
<mesh>boundary_right</mesh>
<type>Dirichlet</type>
<component>0</component>
<parameter>temperature_bc_right</parameter>
</boundary_condition>
</boundary_conditions>
<source_terms>
</source_terms>
</process_variable>
</process_variables>
<nonlinear_solvers>
<nonlinear_solver>
<name>basic_newton</name>
<type>Newton</type>
<max_iter>20</max_iter>
<linear_solver>general_linear_solver</linear_solver>
</nonlinear_solver>
</nonlinear_solvers>
<linear_solvers>
<linear_solver>
<name>general_linear_solver</name>
<lis>-i bicgstab -p ilu -tol 1e-16 -maxiter 10000</lis>
<!-- <eigen>
<solver_type>BiCGSTAB</solver_type>
<precon_type>DIAGONAL</precon_type>
<max_iteration_step>10000</max_iteration_step>
<error_tolerance>1e-25</error_tolerance>
</eigen>-->
<eigen>
<solver_type>SparseLU</solver_type>
<scaling>true</scaling>
<!--
<solver_type>PardisoLU</solver_type>
-->
</eigen>
<petsc>
<prefix>sd</prefix>
<parameters>-sd_ksp_type cg -sd_pc_type bjacobi -sd_ksp_rtol 1e-16 -sd_ksp_max_it 10000</parameters>
</petsc>
</linear_solver>
</linear_solvers>
</OpenGeoSysProject>
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
<UnstructuredGrid>
<Piece NumberOfPoints="3" NumberOfCells="1" >
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="402" offset="0" />
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="0" offset="40" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.2" offset="56" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="160" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="200" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="216" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_GAAAAJIBAAAAAAAAAwAAAAAAAAAAAAAAAAAAAA==CAAAAAAAAAAAAAAASAAAAAAAAAAAAAAAmpmZmZmZuT8AAAAAAAAAAAAAAAAAAAAAmpmZmZmZyT8AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA==GAAAAAEAAAAAAAAAAgAAAAAAAAAAAAAAAAAAAA==CAAAAAMAAAAAAAAAAQAAABU=
</AppendedData>
</VTKFile>
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32">
<UnstructuredGrid>
<Piece NumberOfPoints="3" NumberOfCells="1" >
<PointData>
<DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="400" RangeMax="1002" offset="0" />
</PointData>
<CellData>
<DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="199" RangeMax="199" offset="40" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="2.4" RangeMax="2.4083189158" offset="56" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="160" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="200" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="216" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_GAAAAOoDAAAAAAAAkAEAAAAAAACRAQAAAAAAAA==CAAAAMcAAAAAAAAASAAAADMzMzMzMwNAmpmZmZmZuT8AAAAAAAAAADMzMzMzMwNAAAAAAAAAAAAAAAAAAAAAADMzMzMzMwNAmpmZmZmZyT8AAAAAAAAAAA==GAAAAAEAAAAAAAAAAgAAAAAAAAAAAAAAAAAAAA==CAAAAAMAAAAAAAAAAQAAABU=
</AppendedData>
</VTKFile>
This diff is collapsed.
<?xml version="1.0"?>
<VTKFile type="Collection" version="0.1" byte_order="LittleEndian" compressor="vtkZLibDataCompressor">
<Collection>
<DataSet timestep="0" group="" part="0" file="result_1d_dirichlet_slab_ts_0_t_0.000000.vtu"/>
<DataSet timestep="100000" group="" part="0" file="result_1d_dirichlet_slab_ts_5_t_100000.000000.vtu"/>
</Collection>
</VTKFile>
<?xml version="1.0"?>
<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
<UnstructuredGrid>
<FieldData>
<DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="20" format="appended" RangeMin="45" RangeMax="103" offset="0" />
</FieldData>
<Piece NumberOfPoints="1003" NumberOfCells="200" >
<PointData>
<DataArray type="Float64" Name="HydraulicFlow" format="appended" RangeMin="0" RangeMax="0" offset="84" />
<DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="172" />
<DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="1002" offset="268" />
<DataArray type="Float64" Name="capillary_pressure" format="appended" RangeMin="0" RangeMax="0" offset="2384" />
<DataArray type="Float64" Name="capillary_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="2472" />
<DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="2560" />
<DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="2656" />
<DataArray type="Float64" Name="gas_density" format="appended" RangeMin="999.1" RangeMax="999.1" offset="2772" />
<DataArray type="Float64" Name="gas_pressure" format="appended" RangeMin="0" RangeMax="0" offset="2924" />
<DataArray type="Float64" Name="gas_pressure_interpolated" format="appended" RangeMin="0" RangeMax="0" offset="3012" />
<DataArray type="Float64" Name="liquid_density" format="appended" RangeMin="0" RangeMax="0" offset="3100" />
<DataArray type="Float64" Name="liquid_pressure" format="appended" RangeMin="0" RangeMax="0" offset="3188" />
<DataArray type="Float64" Name="porosity" format="appended" RangeMin="0.16" RangeMax="0.16" offset="3276" />
<DataArray type="Float64" Name="saturation" format="appended" RangeMin="0" RangeMax="0" offset="3428" />
<DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="0" RangeMax="0" offset="3516" />
<DataArray type="Float64" Name="temperature" format="appended" RangeMin="0" RangeMax="273.15" offset="3632" />
<DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="273.15" RangeMax="273.15" offset="3736" />
<DataArray type="Float64" Name="velocity_gas" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="3840" />
<DataArray type="Float64" Name="velocity_liquid" NumberOfComponents="2" format="appended" RangeMin="0" RangeMax="0" offset="3936" />
</PointData>
<CellData>
<DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="4032" />
<DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="199" offset="4100" />
<DataArray type="Float64" Name="saturation_avg" format="appended" RangeMin="0" RangeMax="0" offset="4612" />
</CellData>
<Points>
<DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="2.4083189158" offset="4684" />
</Points>
<Cells>
<DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="9792" />
<DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="13156" />
<DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="13584" />
</Cells>
</Piece>
</UnstructuredGrid>
<AppendedData encoding="base64">
_AQAAAAAAAAAAgAAAAAAAABQAAAAAAAAAHAAAAAAAAAA=eF4z0zPRM9A1tDTWTU80MTc1t7SwMAAALQIEZg==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAALA+AAAAAAAAJgAAAAAAAAA=eF7twTEBAAAAwqD1T20Gf6AAAAAAAAAAAAAAAAAAAADgNz6wAAE=AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAEgYAAAAAAAA=eF412MMWIAgCA8Bp27Zt27Zt27Zt27Zt27Ztu2cPW5NLfUJe8s8//08ABmQgBmYQBmUwBmcIhmQohmYYhmU4hmcERmQkRmYURmU0RmcMxmQsxmYcxmU8xmcCJmQiJmYSJmUyJmcKpmQqpmYapmU6pmcGZmQmZmYWZmU2ZmcO5mQu5mYe5mU+5mcBFmQhFmYRFmUxFmcJlmQplmYZlmU5lmcFVmQlVmYVVmU1VmcN1mQt1mYd1mU91mcDNmQjNmYTNmUzNmcLtmQrtmYbtmU7tmcHdmQndmYXdmU3dmcP9mQv9mYf9mU/9ucADuQgDuYQDuUwDucIjuQojuYYjuU4jucETuQkTuYUTuU0TucMzuQszuYczuU8zucCLuQiLuYSLuUyLucKruQqruYaruU6rucGbuQmbuYWbuU2bucO7uQu7uYe7uU+7ucBHuQhHuYRHuUxHucJnuQpnuYZnuU5nucFXuQlXuYVXuU1XucN3uQt3uYd3uU93ucDPuQjPuYTPuUzPucLvuQrvuYbvuU7vucHfuQnfuYXfuU3fucP/uQv/uYf/uW//K/4AzAgAzEwgzAogzE4QzAkQzE0wzAswzE8IzAiIzEyozAqozE6YzAmYzE24zAu4zE+EzAhEzExkzApkzE5UzAlUzE10zAt0zE9MzAjMzEzszArszE7czAnczE38zAv8zE/C7AgC7Ewi7Aoi7E4S7AkS7E0y7Asy7E8K7AiK7Eyq7Aqq7E6a7Ama7E267Au67E+G7AhG7Exm7Apm7E5W7AlW7E127At27E9O7AjO7Ezu7Aru7E7e7Ane7E3+7Av+7E/B3AgB3Ewh3Aoh3E4R3AkR3E0x3Asx3E8J3AiJ3Eyp3Aqp3E6Z3AmZ3E253Au53E+F3AhF3Exl3Apl3E5V3AlV3E113At13E9N3AjN3Ezt3Art3E7d3And3E393Av93E/D/AgD/Ewj/Aoj/E4T/AkT/E0z/Asz/E8L/AiL/Eyr/Aqr/E6b/Amb/E27/Au7/E+H/AhH/Exn/Apn/E5X/AlX/E13/At3/E9P/AjP/Ezv/Arv/E7f/Anf/E3//Av/+V/gz8AAzIQAzMIgzIYgzMEQzIUQzMMwzIcwzMCIzISIzMKozIaozMGYzIWYzMO4zIe4zMBEzIREzMJkzIZkzMFUzIVUzMN0zId0zMDMzITMzMLszIbszMHczIXczMP8zIf87MAC7IQC7MIi7IYi7MES7IUS7MMy7Icy7MCK7ISK7MKq7Iaq7MGa7IWa7MO67Ie67MBG7IRG7MJm7IZm7MFW7IVW7MN27Id27MDO7ITO7MLu7Ibu7MHe7IXe7MP+7If+3MAB3IQB3MIh3IYh3MER3IUR3MMx3Icx3MCJ3ISJ3MKp3Iap3MGZ3IWZ3MO53Ie53MBF3IRF3MJl3IZl3MFV3IVV3MN13Id13MDN3ITN3MLt3Ibt3MHd3IXd3MP93If9/MAD/IQD/MIj/IYj/MET/IUT/MMz/Icz/MCL/ISL/MKr/Iar/MGb/IWb/MO7/Ie7/MBH/IRH/MJn/IZn/MFX/IVX/MN3/Id3/MDP/ITP/MLv/Ibv/MHf/IXf/MP//Jf/nf0B2BABmJgBmFQBmNwhmBIhmJohmFYhmN4RmBERmJkRmFURmN0xmBMxmJsxmFcxmN8JmBCJmJiJmFSJmNypmBKpmJqpmFapmN6ZmBGZmJmZmFWZmN25mBO5mJu5mFe5mN+FmBBFmJhFmFRFmNxlmBJlmJplmFZlmN5VmBFVmJlVmFVVmN11mBN1mJt1mFd1mN9NmBDNmJjNmFTNmNztmBLtmJrtmFbtmN7dmBHdmJndmFXdmN39mBP9mJv9mFf9mN/DuBADuJgDuFQDuNwjuBIjuJojuFYjuN4TuBETuJkTuFUTuN0zuBMzuJszuFczuN8LuBCLuJiLuFSLuNyruBKruJqruFaruN6buBGbuJmbuFWbuN27uBO7uJu7uFe7uN+HuBBHuJhHuFRHuNxnuBJnuJpnuFZnuN5XuBFXuJlXuFVXuN13uBN3uJt3uFd3uN9PuBDPuJjPuFTPuNzvuBL/g/7ne+4AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAALA+AAAAAAAAJgAAAAAAAAA=eF7twTEBAAAAwqD1T20Gf6AAAAAAAAAAAAAAAAAAAADgNz6wAAE=AQAAAAAAAAAAgAAAAAAAAGB9AAAAAAAANgAAAAAAAAA=eF7twTEBAAAAwqD1T20LL6AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPgZ9YAABAQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAATwAAAAAAAAA=eF7tyLEJgDAARNG1LLOhG1kENEQCrqTF2Qg2VhbvmuO/ZVyb5lLzW359ca211lprrbX+2ve3/PH4PT8455xzzjnnnHPOOf+B9/wJF4qYFw==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAATwAAAAAAAAA=eF7tyDERgDAQAEET8ZMyghiGEAZLUYERnKT56ksqiq1ubrcy2zueukd7anbvvffee++9//pH+jt1RC/OOeecc84555xzzn/gZ3QB2SsxDw==AQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAHwAAAAAAAAA=eF7twYEAAAAAw6D5U5/gBlUBAAAAAAAAAM8AH1gAAQ==AQAAAAAAAAAAgAAAAAAAAGB9AAAAAAAANgAAAAAAAAA=eF7twTEBAAAAwqD1T20LL6AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPgZ9YAABAQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAALQAAAAAAAAA=eF7txTENADAMA7BiGKKiyV/+z4ZiUiX7cfKc6di2bdv25wsAAAAAFrgB3FNPAQAAAAAAAAAAgAAAAAAAAFgfAAAAAAAAKwAAAAAAAAA=eF7txSEBAAAIA7AOJCLNPf0NeCpsZsmp6di2bdu2bdu2bdu27fcCUtnKwg==AQAAAAAAAAAAgAAAAAAAALA+AAAAAAAAJgAAAAAAAAA=eF7twTEBAAAAwqD1T20Gf6AAAAAAAAAAAAAAAAAAAADgNz6wAAE=AQAAAAAAAAAAgAAAAAAAALA+AAAAAAAAJgAAAAAAAAA=eF7twTEBAAAAwqD1T20Gf6AAAAAAAAAAAAAAAAAAAADgNz6wAAE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAEAAAAAAAAAA=eF5jYBgFo2AU4AIAAyAAAQ==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAXgEAAAAAAAA=eF4txddCCAAAAMAQGUWylVFIyajMrBQNo5CQ7MxkZ2Qle2QVIitkZpdQsslveXD3cgEB/7VxW7dzoNu7g4Pc0Z3c2V0c7BB3dTeHurvD3MM93cu93cd93c/9He4ID/BAD/JgRzrKQzzUwxzt4Y5xrEc4ziM9yqM9xvFOcKLHepzHe4InepKTPNlTPNXTPN3JnuEUp3qmZznN6c5wpmd7jud6nrOc7fle4IXO8SLnerGXeKnzvMz5Xu4VXulVXu01XusCr/N6b/BGb/JmF3qLi7zV27zdO7zTu7zbxd7jvd7n/S7xAR/0IR/2EZf6qMt8zMd9wid9yqd9xmd9zudd7gu+6Eu+7CuucKWv+pqvu8o3fNPVvuXbvuO7vuca3/cDP3StH/mxn/ipn/m56/zCL/3Kr/3Gb/3O9W7wezf6gz/6k5vc7M9u8Rd/9Td/9w//9C//9h+3+q//AdFeTb0=AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAFQAAAAAAAAA=eF5jYBgFo2AUjIJRMApIBwAGQAABAQAAAAAAAAAAgAAAAAAAAAheAAAAAAAA1A4AAAAAAAA=eF51m33MlmUZh4lYc84xx5yRlZlZ2ZxjrTHnyIHZl5GZ0VxrrLnmmmtlZWZEVisztTLLvijTCjOxHGNkhlYMUfRVFJWPgMzIiAy/8FsQsZ77Ovfr3Pk7ruf9z4Ofx973ee77uq/rvM9z0qT+z/Nf/vifzv7+hbPH8Z/+ZPi5Hf/uPPKLxngWwTNr28kX//jMxciL9/1LxviXIL/i0PMfPmXFUuTF+7/P8jG/z3LkXznj/iNe9tXfIy/e//1XjPn9VyC/dd1z87Zf9Ufkxft/78oxf+9K5N97zEXHn7ZjFfLi/c9n9ZjPZzXyR11w8OhXvRl58f7nuWbM57kG+f03X/np4/a/DXnx/uc/Mebzn0B+5xvevN+1r7oDefH+97V2zPe1Fvnbz1n580Nm3IW8eP/7XTfm+12H/Ae3X3XZLYfdg7x4/3pgXtzzb//WJx7ce829yIv3r5/1Y66f9ci/pf1sQF68f70xL+751923Z/TNb0RevH99bhpzfW5Cftp5q+745Tv+grx4/3rePOZ63oz85PaHMS/ev/63jLn+tyD/5PqTPnrgh7YiL96/X/6KvLjnH1h40NJ3PcC8eP/+ug95cc/fe/jwRf4NefH+/Xg/8uKev2n0LV7/FPPi/fv378iLe37ZWWd879FztyEv3r/f/4G8uOffP7rrpp3HvHh/fWBe3PPvji8SefH+esK8uOdnz/zIjSce8E/kxfvrD/Pino9lj3nx/nq1HXlxz7dl44fMi/fXN+bFPf/G4bb78L+QF++vh8yLe/41w2V76A7kxfvrJ/Pinn95u7GZF++vt/9GXtzzU49efvmuq5kX76/PDyIv7vnRl3jC6NJFXry/njMv7vkXNy4YLa3/QV68v/4zL+75Z4+cf8nqZ5gX93w8F3YiL+75x84d7lTmxT0fz5GHkBf3/I57hoWeeXHPx7LEvLjn2216wsPIi3s+nlPMi3t+0/A17vcI8uKej+ca8+Kev+vO4Yd5cc/Hc/BR5MU9v+aw9qRCXtzz8dxkXtzzo4t2+pRTH0Ne3PPxnGVe3PPXTQx/2C7kxT0fz2XmxT0/WoRPn7WNeXHPx3P8ceTFPf+B9ocxL+75eO4zL+75943u6mOvY17c87FPYF7c8+85YFj5nkBe3POxr2Be3PPvvP6pfRsXMC/u+diHMC/u+bedNqw0zIt7PvYtTyIv7vnj2oOHeXHPxz6HeXHPH9sudObFPR/7IubFPd/+rMlPIS/u+dhHMS/u+bbNOZ15cc/Hvot5cc+3bcKtzIt7PvZpzIt7/k3DY/bIp5EX93zs65gX9/zr20aEeXHPxz6QeXHPv7YtfMyLez72jc8gL+75Vw/L5MnMi3s+9pnMi3v+FVNOHW0FmRf3fOxLmRf3/LCrWHjQs8iLez72scyLe35YhSfOYV7c87HvZV7c820Z3sq8uOdjn/wc8uKeb7fpLObFPR/7aubFPd++xsuZF/d87MOZF/d8/OxGXtzzsW9nXtzzL1zzxGhnzby452Ofz7y453fPGzYWzIt7Ps4FzIt7/ul9wx+2B3lxz8c5gnlxzz9+9a1zj7mIeXHPx7mDeXHPP3LKsJFiXtzzcU55Hnlxzw+bivUnMS/u+TjXMC/u+bZtXMa8uOfjHMS8uOfbtmLaXuTFPR/nJubFPd8eO2czL+75OGcxL+75KLsyL+75OJe9gLy45zcMl+2xzIt7Ps5xzIt7/u7hz7qMeXHPx7mPeXHPrx2OZfuYF/d8nBP3IS/u+dvmDjtH5sU9H+dK5sU9f/Mzw8LEvLjn4xzKvLjnW9X4iBeRF/d8nFuZF/f8sEs+6gLmxT0f51zmxT3fjkE7mRf3fJyL/4u8uOfbNnku8+Kej3M08+Keb9uopcyLez7O3cyL9+/3SXM8L+75OKczL95fH5gX93yc65kX768nzIt7PuoAzIv31x/mxT0fdQPmxfvrFfPino86A/Pi/fWNeXHPR13iJciL99dD5sU9H3UM5sX76yfz4p6Pugfz4v31lnlxz0edhHnx/vrMvLjno67CvHh/PWde3PNRh5mMvHh//Wde3PNRt2Fe3PPxXGBe3PNR52Fe3PPxHGFe3PNRF2Je3PPx3GFe3PNRR2Je3PPxnGJe3PPxNb4UeXHPx3ONeXHPR52KeXHP6yf4Df/ncV5bgJy4e2ofV3re2r73S+ER73sWwRPf1xXwiPevw8XwHNLqMFfBI97/fZbAE/v538Ij3r8Ol8ITz4Vl8Ij3/67l8ER54HfwiLun9o+lJ+okf4BHvP/5rIDnoXbevBEecffUfrP0tMt8xp/hEe9/zivhif0P+9DE3VP709ITz5Gb4BHvf1+r4Yn7kX1r4u6p/WzpibroLfCI97/3NfBEfe9WeMTdU/vf0hN1KvbLifevnwl4ot7CdUrcPbVfLj1RN2B/nXj/OlwLT+xX74RH3D21vy49q9pzn/144v3reR08cR65Gx5x99R+vPTEvpr9eOL9+4Ke2B+yT0/cPbV/Lz1z2oOQ/Xvi/ftrPTzxa9Ij7p647zbAE+9l2O8n3r9P6Yn3I+wDFHdP7Q9MT7ynYH+geP9+3wTPwe1Cp0fcPbWfMD1Rt2f/oXh/3dgMT9TP2Wco7p7af5ieqGOz/1C8v/5sgSfqyfSIu6f2K6Yn6rpb4RF3T+1jTE/UV9nHKO6e2t+Ynqhzsr9R3D217zE9UX9g36O4e2o/ZHriHE2PuHtqn2R64jzIPklx99T+yfTEuYb9k+LuqX2V6Yn9Ofsqxd1T+y3TE3UhesTdU/sw0xP1GfZVirun9memJ+ok9Ii7p/ZtpifqFezDFHdP7edMT9QN6BF3T+3zTE+c3+kRd0/t/0xPnKPZzynuntoXmp44z9Ij7p7aL5qeOFeyn1PcPbWPND3Rz0CPuHtqf2l6om+BHnH31L7T9ER/AvtIxd1T+1HTE30I9Ii7p/appif6DegRd0/tX01P9BWwH1XcPbWvNT3RP0CPuHtqv2t6ok+A/avi7ql9sOmJfgB6xN1T+2PTE+/96RF3T+2bTU+832cfrLh7aj9teuI9Pj3i7ql9tjfY/8c+WHH31P7b9MR7eXrE3VP7ctMT79/pEXdP7ddNT7xn3wmPuHtqH2964n06PeLuqf296Yn35vSIu6f2/aYn3o+zj1fcPbUfOD3xHpwecffUPuH0xPtu9v2Ku6f2D6cn3mvTI+6e2lecntjm0yPuntpvnJ54D8X+YXH31D7k9MT7IHrE3VP7k9MT72XoEXdP7VtOT7wfYR+yuHtqP3N64j0FPeLuqX3O6Yn3BexbFndP7X9OT9Tt6RF3T+2LTk/Uz+kRd0/tl05PvOdl/7O4e2ofdXrifS494u6p/dXpieuc/cni7ql9135f0CPuntqPnZ64X+gRd0/t0/b7ix5x99T+7fTEfUePuHtqX7ffp/SIu6f2e6cn7l/2V4u7p/aB+/1Oj7h7an94emIdoEfcPbVv3NcNesTdU/vJ0xN1S3rE3VP7zH39YV+3uHtq/3l6Yl2iR9w9tS/d1zF6xN1T+9XTE+sbPeLuqX3svh7SI+6e2t+enlgn2R8u7p7a9+7rKj3i7qn98OmJ9ZYecffUPnlfn+kRd0/tn09PrNv0iLun9tWn5zfDdvtA9qWLu6f226dnybC9PZ4ecffUPvz0/HrYTp5Fj7h7an9+en41bN+upEfcPbVvPz2Lh2V4Iz3i7qn9/On5RXtxy354cffUPv/0XDFctsfQI+6e2v+fnp8NH+cZ9Ii7p84FpCf+mx5x99R5gfQsOnM4kdMj7p46R5CeHw2n3730iLunzhekZ1iFpx7Nfn5x99S5g/SMFofRo4QecffUeYT0DHfd/EvoEXdPnVNIT9sWraJH3D11fiE98bimR9w9da4hPe0xcjjnCMTdU+cd0tOWt3n0iLunzkGkp91259Mj7p46H5GedjlcT4+4e+rcRHqinEOPuHviZzc8XxnKJ9M5jyDunjpnkZ4vtY0RPeLuqfMX6fliu/HoEXdPnctIzxda4YMecffUeY30fL5t4OgRd0+d40jP59oCsQcecffU+Y70fHbYVsymR9w9de4jPZ9pG1B6xN1T50HS86m2wNEj7p46J5KeT7bCEz3i7qnzI+lpl+0Uzl+Iu6fOlaSnfZwz6RF3T503Sc/H2g894u6pcyjpOb1t9OkRd0+dT0lPO7ZO0CPunjq3kp52nNpDj7h76jxLeto2/yjOj4i7p865pKdtP+fTI+6eOv+SnrYtupgecffUuZj0tMf1SnrE3VPnZdLTHiO76BF3T52jSU/UYTiHIu6eOl+TnqjD0CPunjp3k56ow9Aj7p46j5OeqMPQI+6eOqeTnqjD0CPunjq/k56ow3D+Rdw9da4nPVGHoUfcPXXeJz1Rh6FH3D11Dig9UYehR9w9dT4oPVGHoUfcPXVuKD3RT8W5G3H31Hmi9EQdhh5x99Q5o/REHYYecffU+aP0RB2GHnH31Lmk9EQdhh5x99R5pfREHYbzPuLuqXNM6Yk6DD3i7qnzTemJOgw94u6pc0/piToMPeLuqfNQ6Yk6DD3i7qlzUumJbT494u6p81PpiXMr53HE3VPnqtLzg3aepUfcPXXeKj1x/qVH3D11Dis9l7ZzMT3i7qnzWemJczQ94u6pc1vp+W47X9Mj7p46z5WeOI/TI+6eOueVnu+0czo94u6p81/piXM9PeLuqXNh6fl2O+/TI+6eOi+WnqgPcP5I3D11jiw932x1A3rE3VPny9ITdQZ6xN1T587Sc2GrP9Aj7p46j5aeqFfQI+6eOqeWnm+0OgY94u6p82vpiboHPeLuqXNt6fl6q4fQI+6eOu+Wnqif0CPunjoHl56vtboKPeLuqfNx6YmfyXMmdX84b1Xn5tLTyh7TJ89xg7h76jxdeqKeQ4+4e+qcXXrasf5EesTdU+fv0hN1IXrE3VPn8tLTjq0L6RF3T53XS0/Ul+gRd0+d40tPe1xfS4+4e+p8X3qiTkWPuHvq3F96og2ZHnH31HnA9ES9ix5x98wsc4LpaZfDVM79ibunzg+mJ+pm9Ii7p84Vpud/9ahbwA==AQAAAAAAAAAAgAAAAAAAAAAyAAAAAAAAuAkAAAAAAAA=eF5d2UOQrIkWReG+tm3btm3btm3btm3btm37dl9bPXhrv4jeNfkGdfaqzPwnFZF//fW/nwAYEAPhJH4xBSfjxAD/3QXGIKjOVA6m4zTraRcUg6E6M7ibhTOtp11wDIHqzOZuLs6xnnYhMRSqM4+7BTjfetqFxjCozkLuFuMi62kXFsOhOku4W4ZLraddeIyA6iznbiWusJ52ETESqrOKuzW42nraRcYoqM5a7tbjOutpFxWjoTobuNuEG62nXXSMgeps5m4rbrGedjExFqqzjbsduN162sXGOKjOTu524y7raRcX46E6e7jbh3utp118TIDq7OfuIB6wnnYJMRGqc4i7I3jYetolxiSozlHujuMx62mXFJOhOie4O4UnraddckyB6pzm7iyesZ52KTEVqnOOuwt43nrapcY0qM5F7i7jJetplxbToTpXuLuGV62nXXrMgOpc5+4m3rCedhkxE6pzi7s7eNt62mXGLKjOXe7u4z3raZcVs6E6D7h7hA+tp112zIHqPObuKT6xnnY5MReq84y7F/jcetrlxjyozkvuXuMr62mXF/OhOn9z9wb/sZ52+bEAqvOWu/f4znraFcRCqM4H7j7hR+tpVxiLoDqfufuKX6ynXVEshup84+4HfreedsWxBKrzk7vf+Mt62pXEUqjOH+4C6B8aVE+70lgG1QnIfWAMhOppVxbLoTpBuAuGQa2nXXmsgOoE5y4khrCedhWxEqoTirswGNp62lXGKqhOWO7CYzjraVcVq6E6EbiLhBGtp111rIHqROYuKkaxnnY1sRaqE427GBjdetrVxjqoTkzuYmMs62lXF+uhOnG4i4dxraddfWyA6sTnLiEmsJ52DbERqpOIuySY2HraNcYmqE5S7pJjMutp1xSboTopuEuFKa2nXXNsgeqk5i4tprGedi2xFaqTjrsMmN562rXGNqhORu4yYybradcW26E6WbjLhlmtp1177IDqZOcuJ+awnnYdsROqk4u7PJjbetp1xi6oTl7u8mM+62nXFbuhOgW4K4QFraddd+yB6hTmrigWsZ52PbEXqlOMuxJY3Hra9cY+qE5J7kpjKetp1xf7oTpluCuHZa2nXX8cgOqU564iVrCedgNxEKpTibsqWNl62g3GIahOVe6qYzXraTcUh6E6NbirhTWtp91wHIHq1OauLtaxnnYjcRSqU4+7BljfetqNxjGoTkPuGmMj62k3FsehOk24a4ZNrafdeJyA6jTnriW2sJ52E3ESqtOKuzbY2nraTcYpqE5b7tpjO+tpNxWnoToduOuEHa2n3XScgep05q4rdrGedjNxFqrTjbse2N162s3GOahOT+56Yy/raTcX56E6fbjrh32tp918XIDq9OduIA6wnnYLcRGqM4i7ITjYetotxiWozlDuhuMw62m3FJehOiO4G4UjrafdclyB6ozmbiyOsZ52K3EVqjOOuwk43nrarcY1qM5E7ibjJOtptxbXoTpTuJuGU62n3XrcgOpM524mzrCedhtxE6ozi7s5ONt62m3GLajOXO7m4zzrabcVt6E6C7hbhAutp9123IHqLOZuKS6xnnY7cReqs4y7Fbjcetrtxj2ozko9Z1xlPe324j5UZw136/S8rafdfjyA6qzXc8IN1tPuIB5CdTZxt0XPy3raHcYjqM5Wfc64zXraHcVjqM4O7nbp87aedsfxBKqzW58T7rGedifxFKqzj7sD+rysp91pPIPqHNT7xEPW0+4snkN1jnB3TO/XetqdxwuoznG9TjxhPe0u4iVU5xR3Z/R6rafdZbyC6pzV38Fz1tPuKl5DdS5wd0l/z3raXccbqM5l3eEV62l3E2/h/zvc3dC99bS7jXdQnZvc3cZb1tPuLt5Dde5wdw/vWk+7+/gA1bnP3UN8YD3tHuIjVOcRd0/wsfW0e4xPtMOn3D3HZ97THT5DdV5w9wpfWk+757pHdV5z9w/+bT3tXuIrVOcNd+/wrfW0e41/ozrvufuIH6yn3T/4BtX5xN0X/Gw97d7iO1TnK3ff8Zv1tHuPH1CdH9z9wp/W0+4jfkJ1fnOnL5T/WE+7z/gF1QnALhAG1BfTgf67+4rfUJ3A3AXFINbT7jv+QHWCcRcCg1tPu5/4C9UJyV1oDGU97X7jH1QnDHfhMKz1tNMX8QFQnfDcRcQI1tMuIAayTiTuomBk62kXGINYJyp30TGa9bQLisGsE4O7WBjTetoFxxDWic1dXIxjPe1CYijrxOMuAca3nnahMYx1EnKXGBNZT7uwGM46SbhLhkmtp114jGCd5NylxBTW0y6inrd1UnGXBlNbT7vIGMU6ablLj+msp11UjGadDNxlwozW0y46xrBOZu6yYhbraRcTY1knG3c5MLv1tIuNcayTk7vcmMt62sXFeNbJw10+zGs97eJjAuvk564gFrCedgkxkXUKcVcEC1tPu8SYxDpFuSuOxaynXVJMZp0S3JXCktbTLjmmsE5p7spiGetplxJTWaccdxWwvPW0S41prFORu8pYyXrapcV01qnCXTWsaj3t0mMG61TnribWsJ52GTGTdWpxVwdrW0+7zJjFOnW5q4/1rKddVsxmnQbcNcKG1tMuO+awTmPummIT62mXE3NZpxl3LbC59bTLjXms05K71tjKetrlxXzWacNdO2xrPe3yYwHrtOeuI3awnnYFsZB1OnHXBTtbT7vCWMQ6Xbnrjt2sp11RLGadHtz1wp7W0644lrBOb+76Yh/raVcSS1mnH3cDsL/1tCuNZawzkLvBOMh62pXFctYZwt0wHGo97cpjBesM524kjrCedhWxknVGcTcGR1tPu8pYxTpjuRuP46ynXVWsZp0J3E3CidbTrjrWsM5k7qbiFOtpVxNrWWcadzNwuvW0q411rDOTu9k4y3ra1cV61pnD3Tycaz3t6mMD68znbiEusJ52DbGRdRZxtwQXW0+7xtjEOku5W47LrKddU2xmnRXcrcKV1tOuObawzmru1uIa62nXEltZZx13G3C99bRrjW2ss5G7zbjJetq1xXbW2cLdNtxqPe3aYwfrbOduJ+6wnnYdsZN1dnG3B3dbT7vO2MU6e7nbj/usp11X7GadA9wdwoPW06479rDOYe6O4hHradcTe1nnGHcn8Lj1tOuNfaxzkrvTeMp62vXFftY5w905PGs97frjAOuc5+4iXrCedgNxkHUucXcFL1tPu8E4xDpXubuO16yn3VAcZp0b3N3Cm9bTbjiOsM5t7u7iHetpNxJHWecedw/wvvW0G41jrPOQu8f4yHrajcVx1nnC3TN8aj3txuME6zzn7iW+sN6/0B4DhQ==AQAAAAAAAAAAgAAAAAAAAEAGAAAAAAAAHgEAAAAAAAA=eF4txRu0rAAAAMDe94RhGIZhGIZhGC4uhmG4uBiGYRiGYRiGYXihGZkweESOnTh15tyFS1eu3fjlt1t37v3x14NHT569ePXm3YdPX74d/HoKHTl24tSZcxcuXbl245ffbt2598dfDx49efbi1Zt3Hz59+Xbw+yl05NiJU2fOXbh05dqNX367defeH389ePTk2YtXb959+PTl28Gfp9CRYydOnTl34dKVazd++e3WnXt//PXg0ZNnL169effh05dvB3+fQkeOnTh15tyFS1eu3fjlt1t37v3x14NHT569ePXm3YdPX74d/HsKHTl24tSZcxcuXbl245ffbt2598dfDx49efbi1Zt3Hz59+Xbw/yl05NiJU2fOXfgH2lRgNw==AQAAAAAAAAAAgAAAAAAAAMgAAAAAAAAADAAAAAAAAAA=eF4TFx8eAAAPDRH5
</AppendedData>
</VTKFile>