[cig-commits] commit 2446 by heister to /var/svn/dealii/aspect
dealii.demon at gmail.com
dealii.demon at gmail.com
Wed Apr 9 07:50:19 PDT 2014
Revision 2446
merge from mainline
_U branches/freesurface/
A branches/freesurface/include/aspect/boundary_temperature/constant.h
A branches/freesurface/source/boundary_temperature/constant.cc
_U branches/freesurface/tests/
A branches/freesurface/tests/constant_temperature_boundary_conditions/
D branches/freesurface/tests/constant_temperature_boundary_conditions/screen-output
A branches/freesurface/tests/constant_temperature_boundary_conditions/screen-output
D branches/freesurface/tests/constant_temperature_boundary_conditions/solution-00000.0000.vtu
A branches/freesurface/tests/constant_temperature_boundary_conditions/solution-00000.0000.vtu
A branches/freesurface/tests/constant_temperature_boundary_conditions.prm
http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2446&peg=2446
Diff:
Index: branches/freesurface
===================================================================
--- branches/freesurface 2014-04-09 14:49:18 UTC (rev 2445)
+++ branches/freesurface 2014-04-09 14:50:17 UTC (rev 2446)
Property changes on: branches/freesurface
___________________________________________________________________
Modified: svn:mergeinfo
## -2,4 +2,4 ##
/branches/compositional:1141-1251
/branches/fully-nonlinear:542-728
/branches/j-dannberg:1679,1875-1882
-/trunk/aspect:2386-2438
+/trunk/aspect:2386-2445
\ No newline at end of property
Copied: branches/freesurface/include/aspect/boundary_temperature/constant.h (from rev 2445, trunk/aspect/include/aspect/boundary_temperature/constant.h)
===================================================================
--- branches/freesurface/include/aspect/boundary_temperature/constant.h (rev 0)
+++ branches/freesurface/include/aspect/boundary_temperature/constant.h 2014-04-09 14:50:17 UTC (rev 2446)
@@ -0,0 +1,109 @@
+/*
+ Copyright (C) 2014 by the authors of the ASPECT code.
+
+ This file is part of ASPECT.
+
+ ASPECT is free software; you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation; either version 2, or (at your option)
+ any later version.
+
+ ASPECT is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with ASPECT; see the file doc/COPYING. If not see
+ <http://www.gnu.org/licenses/>.
+*/
+
+
+#ifndef __aspect__boundary_temperature_constant_h
+#define __aspect__boundary_temperature_constant_h
+
+#include <aspect/boundary_temperature/interface.h>
+#include <map>
+
+
+namespace aspect
+{
+ namespace BoundaryTemperature
+ {
+ /**
+ * A class that implements a temperature boundary condition for an arbitrary
+ * geometry, where it returns a constant value for a given boundary indicator.
+ * Takes an ordered list of temperatures, where the position in the list
+ * corresponds to the indicator for that temperature.
+ *
+ * @ingroup BoundaryTemperatures
+ */
+ template <int dim>
+ class Constant : public Interface<dim>
+ {
+ public:
+ /**
+ * Return the temperature that is to hold at a particular location on the
+ * boundary of the domain. This function returns the constant temperatures
+ * read from the parameter file.
+ *
+ * @param geometry_model The geometry model that describes the domain. This may
+ * be used to determine whether the boundary temperature model is implemented
+ * for this geometry.
+ * @param boundary_indicator The boundary indicator of the part of the boundary
+ * of the domain on which the point is located at which we are requesting the
+ * temperature.
+ * @param location The location of the point at which we ask for the temperature.
+ **/
+ virtual
+ double temperature (const GeometryModel::Interface<dim> &geometry_model,
+ const unsigned int boundary_indicator,
+ const Point<dim> &location) const;
+
+ /**
+ * Return the minimal the temperature on that part of the boundary
+ * on which Dirichlet conditions are posed.
+ *
+ * This value is used in computing dimensionless numbers such as the
+ * Nusselt number indicating heat flux.
+ */
+ virtual
+ double minimal_temperature (const std::set<types::boundary_id> &fixed_boundary_ids) const;
+
+ /**
+ * Return the maximal the temperature on that part of the boundary
+ * on which Dirichlet conditions are posed.
+ *
+ * This value is used in computing dimensionless numbers such as the
+ * Nusselt number indicating heat flux.
+ */
+ virtual
+ double maximal_temperature (const std::set<types::boundary_id> &fixed_boundary_ids) const;
+
+ /**
+ * Declare the parameters this class takes through input files.
+ * This class declares the inner and outer boundary temperatures.
+ */
+ static
+ void
+ declare_parameters (ParameterHandler &prm);
+
+ /**
+ * Read the parameters this class declares from the parameter
+ * file.
+ */
+ virtual
+ void
+ parse_parameters (ParameterHandler &prm);
+
+ private:
+ /**
+ * Temperatures at the inner and outer boundaries.
+ */
+ std::map< types::boundary_id, double> boundary_temperatures;
+ };
+ }
+}
+
+
+#endif
Copied: branches/freesurface/source/boundary_temperature/constant.cc (from rev 2445, trunk/aspect/source/boundary_temperature/constant.cc)
===================================================================
--- branches/freesurface/source/boundary_temperature/constant.cc (rev 0)
+++ branches/freesurface/source/boundary_temperature/constant.cc 2014-04-09 14:50:17 UTC (rev 2446)
@@ -0,0 +1,178 @@
+/*
+ Copyright (C) 2014 by the authors of the ASPECT code.
+
+ This file is part of ASPECT.
+
+ ASPECT is free software; you can redistribute it and/or modify
+ it under the terms of the GNU General Public License as published by
+ the Free Software Foundation; either version 2, or (at your option)
+ any later version.
+
+ ASPECT is distributed in the hope that it will be useful,
+ but WITHOUT ANY WARRANTY; without even the implied warranty of
+ MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
+ GNU General Public License for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with ASPECT; see the file doc/COPYING. If not see
+ <http://www.gnu.org/licenses/>.
+*/
+
+
+#include <aspect/boundary_temperature/constant.h>
+#include <deal.II/base/utilities.h>
+#include <limits>
+
+
+namespace aspect
+{
+ namespace BoundaryTemperature
+ {
+// ------------------------------ Constant -------------------
+
+ template <int dim>
+ double
+ Constant<dim>::
+ temperature (const GeometryModel::Interface<dim> &geometry_model,
+ const unsigned int boundary_indicator,
+ const Point<dim> &location) const
+ {
+ types::boundary_id bid = boundary_indicator;
+ std::map< types::boundary_id, double>::const_iterator it = boundary_temperatures.find(bid);
+ if (it != boundary_temperatures.end())
+ return it->second;
+ else
+ {
+ Assert (false, ExcMessage ("Unknown boundary indicator."));
+ return std::numeric_limits<double>::quiet_NaN();
+ }
+ }
+
+
+ template <int dim>
+ double
+ Constant<dim>::
+ minimal_temperature (const std::set<types::boundary_id> &fixed_boundary_ids) const
+ {
+ std::map<types::boundary_id, double>::const_iterator it = boundary_temperatures.begin();
+ double min = it->second;
+ ++it;
+
+ for( ; it != boundary_temperatures.end(); ++it)
+ if ( it->second < min )
+ min = it->second;
+
+ return min;
+ }
+
+
+
+ template <int dim>
+ double
+ Constant<dim>::
+ maximal_temperature (const std::set<types::boundary_id> &fixed_boundary_ids) const
+ {
+ std::map<types::boundary_id, double>::const_iterator it = boundary_temperatures.begin();
+ double max = it->second;
+ ++it;
+
+ for( ; it != boundary_temperatures.end(); ++it)
+ if ( it->second > max )
+ max = it->second;
+
+ return max;
+ }
+
+
+
+ template <int dim>
+ void
+ Constant<dim>::declare_parameters (ParameterHandler &prm)
+ {
+ prm.enter_subsection("Boundary temperature model");
+ {
+ prm.enter_subsection("Constant");
+ {
+ prm.declare_entry ("Boundary indicator to temperature mappings", "",
+ Patterns::Map (Patterns::Integer(0, std::numeric_limits<types::boundary_id>::max()),
+ Patterns::Double()),
+ "A comma separated list of mappings between boundary "
+ "indicators and the temperature associated with the "
+ "boundary indicators. The format for this list is "
+ "``indicator1 : value1, indicator2 : value2, ...'', "
+ "where each indicator is a valid boundary indicator "
+ "and each value is the temperature of that boundary." );
+ }
+ prm.leave_subsection ();
+ }
+ prm.leave_subsection ();
+ }
+
+
+ template <int dim>
+ void
+ Constant<dim>::parse_parameters (ParameterHandler &prm)
+ {
+ prm.enter_subsection("Boundary temperature model");
+ {
+ prm.enter_subsection("Constant");
+ {
+ //get the list of mappings
+ const std::vector<std::string> x_boundary_temperatures
+ = Utilities::split_string_list(prm.get ("Boundary indicator to temperature mappings"));
+
+
+ for (std::vector<std::string>::const_iterator it = x_boundary_temperatures.begin();
+ it != x_boundary_temperatures.end(); ++it)
+ {
+ // each entry has the format (white space is optional):
+ // <id> [x][y][z] : <value (might have spaces)>
+ std::string comp = "";
+ std::string value = "";
+
+ std::stringstream ss(*it);
+ int b_id;
+ ss >> b_id; // need to read as int, not char
+ types::boundary_id boundary_id = b_id;
+
+ char c;
+ while (ss.peek()==' ') ss.get(c); // eat spaces
+
+ if (ss.peek() != ':')
+ {
+ Assert(false, ExcMessage("Cannot parse boundary temperature list, format"
+ "``boundary_id : value'' appears to be missing"));
+ }
+ else
+ ss.get(c); // read the ':'
+
+ while (ss.peek()==' ') ss.get(c); // eat spaces
+ std::getline(ss,value); // read until the end of the string
+
+ AssertThrow (boundary_temperatures.find(boundary_id) == boundary_temperatures.end(),
+ ExcMessage ("Boundary indicator <" + Utilities::int_to_string(boundary_id) +
+ "> appears more than once in the list of indicators "
+ "for constant temperature boundary conditions."));
+
+ boundary_temperatures[boundary_id] = Utilities::string_to_double(value);
+ }
+ }
+ prm.leave_subsection ();
+ }
+ prm.leave_subsection ();
+ }
+ }
+}
+
+// explicit instantiations
+namespace aspect
+{
+ namespace BoundaryTemperature
+ {
+ ASPECT_REGISTER_BOUNDARY_TEMPERATURE_MODEL(Constant,
+ "constant",
+ "A model in which the temperature is chosen constant on "
+ "a given boundary indicator. Parameters are read from the "
+ "subsection 'Constant'.")
+ }
+}
Index: branches/freesurface/tests
===================================================================
--- branches/freesurface/tests 2014-04-09 14:49:18 UTC (rev 2445)
+++ branches/freesurface/tests 2014-04-09 14:50:17 UTC (rev 2446)
Property changes on: branches/freesurface/tests
___________________________________________________________________
Added: svn:ignore
## -0,0 +1 ##
+
Copied: branches/freesurface/tests/constant_temperature_boundary_conditions/screen-output (from rev 2445, trunk/aspect/tests/constant_temperature_boundary_conditions/screen-output)
===================================================================
--- branches/freesurface/tests/constant_temperature_boundary_conditions/screen-output (rev 0)
+++ branches/freesurface/tests/constant_temperature_boundary_conditions/screen-output 2014-04-09 14:50:17 UTC (rev 2446)
@@ -0,0 +1,39 @@
+-----------------------------------------------------------------------------
+-- This is ASPECT, the Advanced Simulator for Problems in Earth's ConvecTion.
+-- . running in DEBUG mode
+-- . running with 1 MPI process
+-- . using Trilinos
+-----------------------------------------------------------------------------
+
+Number of active cells: 1,024 (on 6 levels)
+Number of degrees of freedom: 13,764 (8,450+1,089+4,225)
+
+*** Timestep 0: t=0 seconds
+ Solving temperature system... 12 iterations.
+ Rebuilding Stokes preconditioner...
+ Solving Stokes system... 7 iterations.
+
+ Postprocessing:
+ Writing graphical output: output-constant_temperature_boundary_conditions/solution-00000
+ Temperature min/avg/max: -2.3e-16 K, 0.03098 K, 3 K
+ Heat fluxes through boundary parts: 2.083e-06 W, -7.792e-05 W, -0.000228 W, -0.000342 W
+
+Termination requested by criterion: end time
+
+
++---------------------------------------------+------------+------------+
+| Total wallclock time elapsed since start | 0.603s | |
+| | | |
+| Section | no. calls | wall time | % of total |
++---------------------------------+-----------+------------+------------+
+| Assemble Stokes system | 1 | 0.119s | 20% |
+| Assemble temperature system | 1 | 0.0574s | 9.5% |
+| Build Stokes preconditioner | 1 | 0.131s | 22% |
+| Build temperature preconditioner| 1 | 0.00454s | 0.75% |
+| Solve Stokes system | 1 | 0.0185s | 3.1% |
+| Solve temperature system | 1 | 0.00474s | 0.79% |
+| Initialization | 2 | 0.0184s | 3% |
+| Postprocessing | 1 | 0.127s | 21% |
+| Setup dof systems | 1 | 0.0866s | 14% |
++---------------------------------+-----------+------------+------------+
+
Copied: branches/freesurface/tests/constant_temperature_boundary_conditions/solution-00000.0000.vtu (from rev 2445, trunk/aspect/tests/constant_temperature_boundary_conditions/solution-00000.0000.vtu)
===================================================================
--- branches/freesurface/tests/constant_temperature_boundary_conditions/solution-00000.0000.vtu (rev 0)
+++ branches/freesurface/tests/constant_temperature_boundary_conditions/solution-00000.0000.vtu 2014-04-09 14:50:17 UTC (rev 2446)
@@ -0,0 +1,39 @@
+<?xml version="1.0" ?>
+<!--
+# vtk DataFile Version 3.0
+#This file was generated
+-->
+<VTKFile type="UnstructuredGrid" version="0.1" compressor="vtkZLibDataCompressor" byte_order="LittleEndian">
+<UnstructuredGrid>
+<FieldData>
+<DataArray type="Float32" Name="TIME" NumberOfTuples="1" format="ascii">0</DataArray>
+</FieldData>
+<Piece NumberOfPoints="4096" NumberOfCells="1024" >
+ <Points>
+ <DataArray type="Float64" NumberOfComponents="3" format="binary">
+AQAAAACAAQAAgAEAMxsAAA==eNqlncuOK8eVRf1Lsj6gfsb2N1kaGfAPFMAJBxoUIAImLLpBUDRNEAQv3+9HWp7bjRbRfbGSvRhZmzkMLezYOrkjeDIryPub3/z/17e/Xn9++83T60//999/C/7ruOncx1sP46bzW/B/ekv98HoX/v0t9fNcp/UwbnXo/Xr9IPV5nJe8jdfr/2Nh3vr/14+N6/Bc57EOP4F/f0v9mA759zgPvNrCt+M8PNdpPYx/1g+vD+E/Yj/PdVqF8feG9/Fx3r+Cb8d5MB3y7diP1YE+P2I/pkP+4620D9yvruw/P8q+2i2M1/eZXmHe+n7Sa7wPmI7d9//8ev1F1ntzP6ZD3sbr9ekX6lb30y+M1/f/nwvz1uv2c+M6PPdjdegXxpv66ct6/7nxPmD30dYdc9KO82A65NuxH9sH/g2fH7Ef0/kJ/MeLeeiJn35hvJwH0yHfjv08r0O/MN7UT1/2/5/l/jbvB3h1hO+82A+YTuth/LN+7HOJfDf281ynVRj/aPg5/jgv+4RO3A+YDvlO7MfqQJ/d2I/pkO/Geah9WgrfezEPptN6GP+sH1594fuxn+c6rcJ4t+F9fJyXOenFeTAd8r3Yj9WBPvuxH9Mh34/7AVt35DtxP2A6fwXfif3YfSffjf2YDvlu3A/YvOQ7cT9gOuQ7sR+rA/lu7Md0yHfjPNi6Y5/fi/NgOuR7sR/bB+izH/sxHea8H+fB5iXfi/NgOuR7sR+rA/l+7Md0yPffSu8Bvv/uf66/y/ufx+du8jZev8aFeevvc8aN3wOYjj33fQf+/S31Yzrkbfzr9Ydfr0mhbl/nJW/j9fc808K89fcq08Z1MB2rA/n3t9SP6ZB/j/Ngz93k23EeTId8O/Zjn/t/BP8R+zGd++fIP+R9QvM82Lzk23EeTId8O/ZjdSD/EfsxHfI2Xt/nvxT2jfq++qUwXteZF+atr/d5433guY7tA18K42U/pkPex
uvXslC3+r66LIzX3/+vCvPW/79WjetgOlaH+7pYyHpv7sd0yL/HebB19w34dpwH0yHfjv3YPkCfH7Ef0yH/EefB5iXfjvNgOuTbsR+rA/mP2I/pkP+I+wF77ibfifsB0+HnbCf2Y8997HO6sR/TId+N+wGbl3wn7gdMh3wn9mN1IN+N/ZgO+W6cB3vuZk56cR5MhznpxX7sPQD5fuzHdP4Nvh/nweYl34vzYDrke7Efq8PvwPdjP6bzE/h+3A/YumOf0In7AdMh34n9PN8HvjyMf9aP6ZDvxv2AzUu+E/cDpkO+E/uxOpDvxn5Mh3w3zoOtO+akF+fBdMj3Yj+2D9BnP/ZjOuT7cR5sXvK9OA+mQ74X+7E6kO/HfkyHfD8+D8BrIPzgxfMAptN6GP+sH15D4Yexn+c6rcJ4v+Hf8R/n5TmBQXwewHTID2I/Vgf6HMZ+TIf8MM4Dr5HwoxfzYDqth/HP+rG/O5Efx36e67QK48OG9/FxXuZkFOfBdMiPYj9WB/ocx35Mh/w4Pg9g6478ID4PYDr8O/sg9mP3nfww9mM65IfxeQCbl/wgPg9gOuQHsR+rA/lh7Md0yA/jPNi64zmZUZwH0yE/iv3YPkCf49iP6TDn4zgPNi/5UZwH0yE/iv1YHciPYz+mQ34c9wO8JsJPXuwHTKf1MP5ZP7ymwk9jP891WoXxccPP8cd52SdM4n7AdMhPYj9WB/qcxn5Mh/w0zgOvmfCzF/NgOq2H8c/64TUXfh77ea7TKoxPG97Hx3mZk1mcB9MhP4v9WB3ocx77MR3y87gfsHVHfhL3A6bDz9lJ7MfuO/lp7Md0yE/jfsDmJT+J+wHTIT+J/VgdyE9jP6ZDfhrnwdYd++RZnAfTIT+L/dg+QJ/z2I/pMOfzOA82L/lZnAfTIT+L/VgdyM9jP6ZDfh6fB7DnbvKD+DyA6ZAfxH7suY/8MPZjOuSH8XkAm5f8ID4PYDp/BT+I/VgdyA9jP6ZDfh
jnwZ67yY/iPJgO+VHsxz73yY9jP6ZDfhznweYlP4rzYDrkR7EfqwP5cezHdMiP4/MAtu54TmAQnwcwHfKD2I/dd/LD2I/pkB/G5wFsXvKD+DyA6ZAfxH6sDuSHsR/TIT+M82DrjvwozoPpkB/FfmwfID+O/ZgO+XGcB5uX/CjOg+mQH8V+rA7kx7Ef0yE/jvsBe+4mP4n7AdMhP4n92HMf+Wnsx3R4Tnga9wM2L/lJ3A+YDvlJ7MfqQH4a+zEd8tM4D/bcTX4W58F0WOdZ7MfeA5Cfx35Mh/w8zoPNS34W58F0yM9iP1YH8vPYj+mwz5/H/YCtO/KTuB8wHfKT2I/dd/LT2I/pkJ/G/YDNS34S9wOmQ34S+7E6kJ/GfkyH/DTOg6078rM4D6ZDfhb7sX2Az+nz2I/pkJ/HebB5yc/iPJgO+Vnsx+pAfh77MR3y89rv/6zkd2DW8vsMK/ndlXVh/Ot1/12dTWHe+u/5bBr/DoDp2Pe+yb+/pX5Mh7yNf73+9mtOtoW6fZ2XvI3Xr31h3vrvOewb18F0rA7k399SP6ZD/j3Og33vnnw7zoPpkG/Hfux7f+Q/Yj+m8x34jzgPNu9/gW/HeTAd8u3Yj9WhB/4j9mM65G3863X/XDgU9o36v/NyKIzX9/ljYd76vnpsvA+Yjt138u9vqR/TIW/jX697Tk6FutV/d+5UGP963XNyLsxb/92ec+M6mI7Vgfz7W+rHdMi/x3mwdUe+HefBdMi3Yz+2D5D/iP2YDvmPOA827x/Bt+M8mA75duzH6kCfH7Ef0yH/EfcD9r178p24HzAd8p3Yj33vm3w39mM65LtxP2Dzku/E/YDpkO/EfqwO5LuxH9Mh343zYN+7J9+L82A6fwTfi/3Y3/1/D74f+zEd8v04DzYv+V6cB9Mh34v9WB3I92M/pnN/L7ST33No3g/YuiPfifsB0yHfif3YfSffjf2YDvlu3A/YvOwTOnE/YDrkO7EfqwN9dmM
/pkO+G+fB1t1vwffiPJgO+V7sx/YB8v3Yj+mQ78d5sHmZk16cB9Mh34v9WB3osx/7MR3yNv71uufkUnhvUP/d3UthvP53hFth3vp7+1vj9wCmY8995N/fUj+mQ97G6/WpCnWr+6kK4/X7WBXmrdetalyH536sDlVhvKmfSp73q8bvAew+2nP3v8G34zyYDvl27Mc+9+nzI/ZjOuQ/XszDTfxUhfFyHkyHfPtFP1aHqjDe1E8l/V5pvP6e+Z+FfaP+XvefhfH69Uth3vp7+18a7wOmY/ed/Ptb6sd0yNt4/e+8vxTqVv+7amm8/neWfxXmrf9d41+N62A6Vgfy72+pH9Mh/x7nwdYd+XacB9Mh34792D5A/iP2YzrkP+I82Lzk23EeTId8O/ZjdSD/EfsxHfIfcT9gz93sEzpxP2A65DuxH3vuI9+N/ZjON+C7L/YDN/FTFcbL/YDpkO/Efp7XoSqMN/VTyfN+1fg9gN1He+4m34vzYDqscy/2Y+8BuE77sR/TId9/MQ838VMVxst5MB3yvRf93KQ+VWG8qZ9K3vOUxsv9gK078p24HzAd8p3Yj9138t3Yj+mQ78b9gM1LvhP3A6ZDvhP7sTqQ78Z+TId8N86DrTvyvTgPpkO+F/uxfYB8P/ZjOuT7cR5sXvK9OA+mQ74X+7E6kO/HfkyHfD8+D2Dfu+c5gUF8HsB0yA9iP/a9b/ocxn5Mh/wwPg9g85IfxOcBTId/ZxzEfqwO5IexH9MhP4zzYN+7Z05GcR5Mh/wo9mPf+yM/jv2YDs/5jOM82LzkR3EeTIf8KPZjdSA/jv2YDvlxfB7A1h35QXwewHTID2I/dt95TmMY+zEd8sP4PIDNy3MCg/g8gOmQH8R+rA70OYz9mA75YZwHW3fkR3EeTIf8KPZj+wD5cezHdMiP4zzYvMzJKM6D6ZAfxX6sDvQ5jv2YDvlx3A/Y9+7JT+J+wHTIT2I/9r1v8tPYj+mQn8b9gM1LfhL3
A6ZDfhL7sTqQn8Z+TIf8NM6Dfe+e/CzOg+mwz5zFfux3APicMo/9mA75eZwHm5f8LM6D6ZCfxX6sDuTnsR/TYZ8/j/sBW3fkJ3E/YDrkJ7Efu+/kp7Ef0yE/jfsBm5d9wiTuB0yH/CT2Y3Wgz2nsx3TIT+M82Lrjc+IszoPpkJ/FfmwfID+P/ZgO+XmcB5uXOZnFeTAd8rPYj9WBPuexH9MhP4/PA9hzN88JDOLzAKZDfhD7sec+8sPYj+nwnMPwxfMAN/FTFcbL5wFMh/wg9vO8DlVhvKmfSp73q8bvAew+2nM3+VGcB9NhnUexH/vc5zodx35Mh/z4xTzcxE9VGC/nwXTIj170c5P6VIXxpn4q6fdK4+XzALbuyA/i8wCmQ34Q+7H7Tn4Y+zEd8sP4PIDNS34QnwcwHfKD2I/Vgfww9mM65IdxHmzdkR/FeTAd8qPYj+0D5MexH9MhP47zYPOSH8V5MB3yo9iP1YH8OPZjOuTHcT9gz93sEyZxP2A65CexH3vuIz+N/ZgO+5zpi/3ATfxUhfFyP2A65Cexn+d1qArjTf1U8rxfNX4PYPfRnrvJz+I8mA7rPIv92HsArtN57Md0yM9fzMNN/FSF8XIeTIf87EU/N6lPVRhv6qeS9zyl8XI/YOuO/CTuB0yH/CT2Y/ed/DT2Yzrkp3E/YPOSn8T9gOmQn8R+rA7kp7Ef0yE/jfNg6478LM6D6ZCfxX5sHyA/j/2YDvl5nAebl/wszoPpkJ/FfqwO5OexH9MhP5ffFah/7vy5wC+EXxR17uOtt+fzLoRfxH54LYVfxn6e67QK4/V/r+eHQt0WtX8H5IeH8ec6P8q/B/SD/H/92LgOz3Ue60Cfy9iP6ZBfxnngtRJ+9WIeTKf1MP5ZP7zWwq9jP891WoXxZcP7+Dgvc7KK82A65FexH6sDfa5jP6ZDfl3cB+7PfX8p7BsL4RdFnbvPXmHe+n7Sa7wPmI7dd/LL2I/pkF8W63Af7
sv+/Dgv+UVRx/4dQNMhv4j9WB3IL2M/pkN+GefB1t1/wK/iPJgO+VXsx/YB+lzHfkyHOV/HebB5ya/iPJgO+VXsx+pAfh37MR3y67gf4LURfvNiP2A6rYfxz/rhtRV+G/t5rtMqjK8bfo4/zss+YRP3A6ZDfhP7sTrQ5zb2Yzrkt3EeeO2E372YB9NpPYx/1g+vvfD72M9znVZhfNvwPj7Oy5zs4jyYDvld7MfqQJ/72I/pkN/H/YCtO/KbuB8wHX7ObmI/dt/Jb2M/pkN+G/cDNi/5TdwPmA75TezH6kB+G/sxHfLbOA+27tgn7+I8mA75XezH9gH63Md+TIc538d5sHnJ7+I8mA75XezH6kB+H/sxHfL74nuA+/cO/l54b7AQflHUuefkH4V5F8IvYj/23Ed+GfsxHfLLYh1+9+vf0SaFui2EXxR17N8BNx3+O+OL2I/Vgfwy9mM65JdxHuy5m/wqzoPpkF/Ffuxzn/w69mM65NdxHmxe8qs4D6ZDfhX7sTqQX8d+TIf8urgP3HPyRd4/P667b8Avijr3z6m5nAf7Iu915433AdOx+05+GfsxHfLLYh3ufdeiULeF8Iuijv0OgOmQX8R+rA7kl7Ef0yG/jPNg6478Ks6D6ZBfxX5sHyC/jv2YDvl1nAebl/wqzoPpkF/FfqwO5NexH9Mhv477AXvuJr+J+wHTIb+J/dhzH/lt7Md0/gN+G/cDNi/5TdwPmA75TezH6kB+G/sxHfLbOA/23E1+F+fBdFjnXezH3gOQ38d+TIf8Ps6DzUt+F+fBdMjvYj9WB/L72I/psM/fx/2ArTvym7gfMB3ym9iP3Xfy29iP6ZDfxv2AzUt+E/cDpkN+E/uxOpDfxn5Mh/w2zoOtO/K7OA+mQ34X+7F9gM/p+9iP6ZDfx3mwecnv4jyYDvld7MfqQH4f+zEd8vv4PACvg/CHF88DmE7rYfyzfngdhT/Gfp7rtArj+4Z/x3+cl+cEDvF5AN
Mhf4j9WB3o8xj7MR3yxzgPvE7Cn17Mg+m0HsY/64fXWfhz7Oe5Tqswfmx4Hx/nZU5OcR5Mh/wp9mN1oM9z7Md0yJ/j8wC27sgf4vMApsO/sx9iP3bfyR9jP6ZD/hifB7B5yR/i8wCmQ/4Q+7E6kD/GfkyH/DHOg607npM5xXkwHfKn2I/tA/R5jv2YDnN+jvNg85I/xXkwHfKn2I/Vgfw59mM65M9xP8DrIvzlxX7AdFoP45/1w+sq/DX281ynVRg/N/wcf5yXfcIl7gdMh/wl9mN1oM9r7Md0yF/jPPC6CX97MQ+m03oY/6wfXlXt9wH+l69iP6bzLfgqzoPN2wN/i/NgOuRvsR+rA/kq9mM634Kv4n7A1t396j6Mf7YfMJ0e+Evsx+47+5xr7Md0yF9f7Ad64qdfGC/3A6ZD/vKin57Up18Yb+qnL+v958b7gN1HW3fMyS3Og+mQv8V+bB8gX8V+TId89WIeeuKnXxgv58F0yN9e9GN16BfGm/rpy/7/s+z/zc8D2HM3+UN8HsB0eE7gEPux577vwB9jP6Zzv8bynqH5eQCb9/fgD/F5ANPpgT/EfqwOfwB/jP2YDvljnAd77iZ/ivNgOuRPsR/73GfOz7Ef02FOznEebF7ypzgPpkP+FPuxOpA/x35Mh/w5Pg9g647nBA7xeQDTIX+I/TzfB748jH/Wj+mQP8bnAWxe8of4PIDpkD/EfqwO5I+xH9Mhf4zzYOuOOTnFeTAd8qfYj+0D9HmO/ZgO+XOcB5uX/CnOg+mQP8V+rA7kz7Ef0yF/jvsBe+4mf4n7AdPh5+wl9mPPfexzrrEf0yF/jfsBm5f8Je4HTIf8JfZjdSB/jf2YDvlrnAd77v4e/C3Og+mQv8V+7D0A+Sr2YzrkqzgPNi/5W5wH0yF/i/1YHchXsR/TIV/F/YCtO/YJl7gfMB3yl9jP833gy8P4Z/2YDvlr3A/YvOQvcT9gOuQvsR+rA/lr7Md0yF/jPNi
6I3+L82A65G+xH9sHyFexH9MhX8V5sHnv11L+/tI8D6ZD/hb7sTqQr2I/pvMt+Kr4OwD3967rwu8GLGqfO+uHcdO599Ub+fcm1vK7OpvGvwNgOva9b/pcxn5Mh/yyWId7TraFui2EXxR17N8BNx2+J1rEfqwO5JexH9Mhv4zzYN+7Z05WcR5Mh/wq9mPf+yO/jv2Yznfg13EebF7yqzgPpkN+FfuxOpBfx35Mh/y6uA/cc3KQf/9lL/vqoTBe3+ePhXkXwi9iP3bfvwW/jP2YDvllsQ73nJwKdavvq6eHcdOxfwfQdMgvYj9WB/pcxn5Mh/wyzoOtO/KrOA+mQ34V+7F9gPw69mM65NdxHmxe5mQV58F0yK9iP1YH+lzHfkyH/DruB+x79+Q3cT9gOuQ3sR/73jf5bezHdMhv437A5iW/ifsB0yG/if1YHchvYz+mQ34b58G+d09+F+fBdNhn7mI/9nd/PqfsYz+mQ34f58HmJb+L82A65HexH6sD+X3sx3TY5+/jfsDWHflN3A+YDvlN7MfuO/lt7Md0yG/jfsDmZZ+wifsB0yG/if1YHehzG/sxHfLbOA+27vicuIvzYDrkd7Ef2wfI72M/pkN+H+fB5mVOdnEeTIf8LvZjdaDPfezHdMjvi+8B7jm5FN4b1H/X/SK/x3uR9zw3+X37i7zPuTV+D2A69txHfhn7MZ1vwC+LdbiPV4W6LYRfFHXs3wE3HfKL2M/zOlSF8aZ+Knnerxq/B7D7aM/d5FdxHkyHdV7Ffuxzn+t0HfsxHfLrF/NwEz9VYbycB9Mhv3rRz03qUxXGm/qppN8rjdffM/+zsG8shF8Ude7XL4V5F8IvYj9238kvYz+mQ35ZrMP9HN0vhbothF8Ude596b8K8y6EX8R+rA7kl7Ef0yG/jPNg6478Ks6D6ZBfxX5sHyC/jv2YDvl1nAebl/wqzoPpkF/FfqwO5NexH9Mhv477AXvuZp+wifsB0yG/if3Y
cx/5bezHdNjnbF/sB27ipyqMl/sB0yG/if08r0NVGG/qp5Ln/arxewC7j/bcTX4X58F0WOdd7MfeA3Cd7mM/pkN+/2IebuKnKoyX82A65Hcv+rlJfarCeFM/lbznKY2X+wFbd+Q3cT9gOuQ3sR+77+S3sR/TIb+N+wGbl/wm7gdMh/wm9mN1IL+N/ZgO+W2cB1t35HdxHkyH/C72Y/sA+X3sx3TI7+M82Lzkd3EeTIf8LvZjdSC/j/2YDvl9fB7AvnfPcwKH+DyA6ZA/xH7se9/0eYz9mA75Y3wewOYlf4jPA5gOuUPsx+pA/hj7MR3yxzgP9r175uQU58F0yJ9iP/a9P/Ln2I/p8JzPOc6DzUv+FOfBdMifYj9WB/Ln2I/pkD/H5wFs3ZE/xOcBTIf8IfZj953nNI6xH9Mhf4zPA9i8PCdwiM8DmA75Q+zH6kCfx9iP6ZA/xnmwdUf+FOfBdMifYj+2D5A/x35Mh/w5zoPNy5yc4jyYDvlT7MfqQJ/n2I/pkD/H/YB97578Je4HTOd78JfYj33vm/w19mM65K9xP2Dzkr/E/YDpkL/EfqwOfwN/jf2YDvlrnAf73j35W5wH0yF/i/3Y7wCQr2I/pkO+ivNg85K/xXkwHfK32I/VgXwV+zEd8lXcD9i6I3+J+wHTIX+J/dh9J3+N/ZgO+WvcD9i87BMucT9gOuQvsR+rA31eYz+mQ/4a58HWHflbnAfTIX+L/dg+QL6K/ZgO+SrOg81L/hbnwXTI32I/VgfyVezHdMhX8XkAe+7mOYFDfB7AdMgfYj/23Ef+GPsxHZ5zOL54HuAmfqrCePk8gOmQP8R+ntehKow39VPJ837V+D2A3Ud77iZ/ivNgOqzzKfZjn/tcp+fYj+mQP7+Yh5v4qQrj5TyYDvnTi35uUp+qMN7UTyX9Xmm8fB7A1h35Q3wewHTIH2I/dt/JH2M/pkP+GJ8HsHnJH+LzAKZD/hD7sTqQP8Z+TIf8Mc6Dr
TvypzgPpkP+FPuxfYD8OfZjOuTPcR5sXvKnOA+mQ/4U+7E6kD/HfkyH/DnuB+y5m33CJe4HTIf8JfZjz330eY39mA7564v9wE38VIXxcj9gOuQvL/qxOlSF8aZ+Knnerxq/B7D7aM/d5G9xHkyH/C32Y+8ByFexH9MhX72Yh5v4qQrj5TyYDvnbi36sDlVhvKmfSt7zlMbL/YCtO/KXuB8wHfKX2I/dd/LX2I/pkL/G/YDNS/4S9wOmQ/4S+7E6kL/GfkyH/DXOg6078rc4D6ZD/hb7sX2AfBX7MR3yVZwHm5f8Lc6D6ZC/xX6sDuSr2I/pkP86/t/yc/Zs
+ </DataArray>
+ </Points>
+
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" format="binary">
+AQAAAABAAAAAQAAAVxYAAA==eNot29MSIMq2BMAztm3btm3btm3btm3btm3b9szNiH0f8hM6unutqv/973//C0AgAhKYIAQjKMEJQShCEpowhCMs4YlAJCISmShEIyrRiUEsYhKbOMQjLvFJQCISkpgkJCMpyUlBKlKSmjSkIy3pyUAmMpKZLGQjK9nJQS5ykps85CMv+SlAIQpSmCIUoyjFKUEpSlKaMpSjLOWpQCUqUpkqVKMq1alBLWpSmzrUoy71aUAjGtKYJjSjKc1pQSta0po2tKMt7elAJzrSmS50oyvd6UEvetKbPvSjL/0ZwCAGMpghDGMowxnBKEYymjGMYyzjmcAkJjKZKUxjKtOZwSxmMps5zGMu81nAIhaymCUsYynLWcEqVrKaNaxjLevZwCY2spktbGMr29nBLnaymz3sYy/7OcAhDnKYIxzjKMc5wSlOcpoznOMs57nAJS5ymStc4yrXucEtbnKbO9zjLvd5wCMe8pgnPOMpz3nBK17ymje84y3v+cAnPvKZL3zjK9/5wS9+8ps//OPv//47/AEIREACE4RgBCU4IQhFSEIThnCEJTwRiEREIhOFaEQlOjGIRUxiE4d4xCU+CUhEQhKThGQkJTkpSEVKUpOGdKQlPRnIREYyk4VsZCU7OchFTnKTh3zkJT8FKERBClOEYhSlOCUoRUlKU4ZylKU8FahERSpThWpUpTo1qEVNalOHetSlPg1oREMa04RmNKU5LWhFS1rThna0pT0d6ERHOtOFbnSlOz3oRU9604d+9KU/AxjEQAYzhGEMZTgjGMVIRjOGcYxlPBOYxEQmM4VpTGU6M5jFTGYzh3nMZT4LWMRCFrOEZSxlOStYxUpWs4Z1rGU9G9jERjazhW1sZTs72MVOdrOHfexlPwc4xEEOc4RjHOU4JzjFSU5zhnOc5TwXuMRFLnOFa1zlOje4xU1uc4d73OU+D3jEQx7zhGc85TkveMVLXvOGd7zlP
R/4xEc+84VvfOU7P/jFT37zh3/8DfDfpR+AQAQkMEEIRlCCE4JQhCQ0YQhHWMITgUhEJDJRiEZUohODWMQkNnGIR1zik4BEJCQxSUhGUpKTglSkJDVpSEda0pOBTGQkM1nIRlayk4Nc5CQ3echHXvJTgEIUpDBFKEZRilOCUpSkNGUoR1nKU4FKVKQyVahGVapTg1rUpDZ1qEdd6tOARjSkMU1oRlOa04JWtKQ1bWhHW9rTgU50pDNd6EZXutODXvSkN33oR1/6M4BBDGQwQxjGUIYzglGMZDRjGMdYxjOBSUxkMlOYxlSmM4NZzGQ2c5jHXOazgEUsZDFLWMZSlrOCVaxkNWtYx1rWs4FNbGQzW9jGVrazg13sZDd72Mde9nOAQxzkMEc4xlGOc4JTnOQ0ZzjHWc5zgUtc5DJXuMZVrnODW9zkNne4x13u84BHPOQxT3jGU57zgle85DVveMdb3vOBT3zkM1/4xle+84Nf/OQ3f/jH34D/PfwDEIiABCYIwQhKcEIQipCEJgzhCEt4IhCJiEQmCtGISnRiEIuYxCYO8YhLfBKQiIQkJgnJSEpyUpCKlKQmDelIS3oykImMZCYL2chKdnKQi5zkJg/5yEt+ClCIghSmCMUoSnFKUIqSlKYM5ShLeSpQiYpUpgrVqEp1alCLmtSmDvWoS30a0IiGNKYJzWhKc1rQipa0pg3taEt7OtCJjnSmC93oSnd60Iue9KYP/ehLfwYwiIEMZgjDGMpwRjCKkYxmDOMYy3gmMImJTGYK05jKdGYwi5nMZg7zmMt8FrCIhSxmCctYynJWsIqVrGYN61jLejawiY1sZgvb2Mp2drCLnexmD/vYy34OcIiDHOYIxzjKcU5wipOc5gznOMt5LnCJi1zmCte4ynVucIub3OYO97jLfR7wiIc85gnPeMpzXvCKl7zmDe94y3s+8ImPfOYL3/jKd37wi5/85g//+Bvovw9/AAIRkMAEIRhBCU4IQh
GS0IQhHGEJTwQiEZHIRCEaUYlODGIRk9jEIR5xiU8CEpGQxCQhGUlJTgpSkZLUpCEdaUlPBjKRkcxkIRtZyU4OcpGT3OQhH3nJTwEKUZDCFKEYRSlOCUpRktKUoRxlKU8FKlGRylShGlWpTg1qUZPa1KEedalPAxrRkMY0oRlNaU4LWtGS1rShHW1pTwc60ZHOdKEbXelOD3rRk970oR996c8ABjGQwQxhGEMZzghGMZLRjGEcYxnPBCYxkclMYRpTmc4MZjGT2cxhHnOZzwIWsZDFLGEZS1nOClaxktWsYR1rWc8GNrGRzWxhG1vZzg52sZPd7GEfe9nPAQ5xkMMc4RhHOc4JTnGS05zhHGc5zwUucZHLXOEaV7nODW5xk9vc4R53uc8DHvGQxzzhGU95zgte8ZLXvOEdb3nPBz7xkc984Rtf+c4PfvGT3/zhH3//f9gXgEAEJDBBCEZQghOCUIQkNGEIR1jCE4FIRCQyUYhGVKITg1jEJDZxiEdc4pOARCQkMUlIRlKSk4JUpCQ1aUhHWtKTgUxkJDNZyEZWspODXOQkN3nIR17yU4BCFKQwRShGUYpTglKUpDRlKEdZylOBSlSkMlWoRlWqU4Na1KQ2dahHXerTgEY0pDFNaEZTmtOCVrSkNW1oR1va04FOdKQzXehGV7rTg170pDd96Edf+jOAQQxkMEMYxlCGM4JRjGQ0YxjHWMYzgUlMZDJTmMZUpjODWcxkNnOYx1zms4BFLGQxS1jGUpazglWsZDVrWMda1rOBTWxkM1vYxla2s4Nd7GQ3e9jHXvZzgEMc5DBHOMZRjnOCU5zkNGc4x1nOc4FLXOQyV7jGVa5zg1vc5DZ3uMdd7vOARzzkMU94xlOe84JXvOQ1b3jHW97zgU985DNf+MZXvvODX/zkN3/4x98g/w35AxCIgAQmCMEISnBCEIqQhCYM4QhLeCIQiYhEJgrRiEp0YhCLmMQmDvGIS3wSkIiEJCYJyUh
KclKQipSkJg3pSEt6MpCJjGQmC9nISnZykIuc5CYP+chLfgpQiIIUpgjFKEpxSlCKkpSmDOUoS3kqUImKVKYK1ahKdWpQi5rUpg71qEt9GtCIhjSmCc1oSnNa0IqWtKYN7WhLezrQiY50pgvd6Ep3etCLnvSmD/3oS38GMIiBDGYIwxjKcEYwipGMZgzjGMt4JjCJiUxmCtOYynRmMIuZzGYO85jLfBawiIUsZgnLWMpyVrCKlaxmDetYy3o2sImNbGYL29jKdnawi53sZg/72Mt+DnCIgxzmCMc4ynFOcIqTnOYM5zjLeS5wiYtc5grXuMp1bnCLm9zmDve4y30e8IiHPOYJz3jKc17wipe85g3veMt7PvCJj3zmC9/4ynd+8Iuf/OYP//gb9L9FXwACEZDABCEYQQlOCEIRktCEIRxhCU8EIhGRyEQhGlGJTgxiEZPYxCEecYlPAhKRkMQkIRlJSU4KUpGS1KQhHWlJTwYykZHMZCEbWclODnKRk9zkIR95yU8BClGQwhShGEUpTglKUZLSlKEcZSlPBSpRkcpUoRpVqU4NalGT2tShHnWpTwMa0ZDGNKEZTWlOC1rRkta0oR1taU8HOtGRznShG13pTg960ZPe9KEffenPAAYxkMEMYRhDGc4IRjGS0YxhHGMZzwQmMZHJTGEaU5nODGYxk9nMYR5zmc8CFrGQxSxhGUtZzgpWsZLVrGEda1nPBjaxkc1sYRtb2c4OdrGT3exhH3vZzwEOcZDDHOEYRznOCU5xktOc4RxnOc8FLnGRy1zhGle5zg1ucZPb3OEed7nPAx7xkMc84RlPec4LXvGS17zhHW95zwc+8ZHPfOEbX/nOD37xk9/84R9/g/234A9AIAISmCAEIyjBCUEoQhKaMIQjLOGJQCQiEpkoRCMq0YlBLGISmzjEIy7xSUAiEpKYJCQjKclJQSpSkpo0pCMt6clAJjKSmSxkIyvZyUEucpKbPOQjL/kpQCEK
UpgiFKMoxSlBKUpSmjKUoyzlqUAlKlKZKlSjKtWpQS1qUps61KMu9WlAIxrSmCY0oynNaUErWtKaNrSjLe3pQCc60pkudKMr3elBL3rSmz70oy/9GcAgBjKYIQxjKMMZwShGMpoxjGMs45nAJCYymSlMYyrTmcEsZjKbOcxjLvNZwCIWspglLGMpy1nBKlaymjWsYy3r2cAmNrKZLWxjK9vZwS52sps97GMv+znAIQ5ymCMc4yjHOcEpTnKaM5zjLOe5wCUucpkrXOMq17nBLW5ymzvc4y73ecAjHvKYJzzjKc95wSte8po3vOMt7/nAJz7ymS984yvf+cEvfvKbP/zj7/+HewIQiIAEJgjBCEpwQhCKkIQmDOEIS3giEImIRCYK0YhKdGIQi5jEJg7xiEt8EpCIhCQmCclISnJSkIqUpCYN6UhLejKQiYxkJgvZyEp2cpCLnOQmD/nIS34KUIiCFKYIxShKcUpQipKUpgzlKEt5KlCJilSmCtWoSnVqUIua1KYO9ahLfRrQiIY0pgnNaEpzWtCKlrSmDe1oS3s60ImOdKYL3ehKd3rQi570pg/96Et/BjCIgQxmCMMYynBGMIqRjGYM4xjLeCYwiYlMZgrTmMp0ZjCLmcxmDvOYy3wWsIiFLGYJy1jKclawipWsZg3rWMt6NrCJjWxmC9vYynZ2sIud7GYP+9jLfg5wiIMc5gjHOMpxTnCKk5zmDOc4y3kucImLXOYK17jKdW5wi5vc5g73uMt9HvCIhzzmCc94ynNe8IqXvOYN73jLez7wiY985gvf+Mp3fvCLn/zmD//4G+K/UF8AAhGQwAQhGEEJTghCEZLQhCEcYQlPBCIRkchEIRpRiU4MYhGT2MQhHnGJTwISkZDEJCEZSUlOClKRktSkIR1pSU8GMpGRzGQhG1nJTg5ykZPc5CEfeclPAQpRkMIUoRhFKU4JSlGS0pShHGUpTwUqUZHKVKEaValODWpRk9rUoR51q
U8DGtGQxjShGU1pTgta0ZLWtKEdbWlPBzrRkc50oRtd6U4PetGT3vShH33pzwAGMZDBDGEYQxnOCEYxktGMYRxjGc8EJjGRyUxhGlOZzgxmMZPZzGEec5nPAhaxkMUsYRlLWc4KVrGS1axhHWtZzwY2sZHNbGEbW9nODnaxk93sYR972c8BDnGQwxzhGEc5zglOcZLTnOEcZznPBS5xkctc4RpXuc4NbnGT29zhHne5zwMe8ZDHPOEZT3nOC17xkte84R1vec8HPvGRz3zhG1/5zg9+8ZPf/OEff0P+F+wNQCACEpggBCMowQlBKEISmjCEIyzhiUAkIhKZKEQjKtGJQSxiEps4xCMu8UlAIhKSmCQkIynJSUEqUpKaNKQjLenJQCYykpksZCMr2clBLnKSmzzkIy/5KUAhClKYIhSjKMUpQSlKUpoylKMs5alAJSpSmSpUoyrVqUEtalKbOtSjLvVpQCMa0pgmNKMpzWlBK1rSmja0oy3t6UAnOtKZLnSjK93pQS960ps+9KMv/RnAIAYymCEMYyjDGcEoRjKaMYxjLOOZwCQmMpkpTGMq05nBLGYymznMYy7zWcAiFrKYJSxjKctZwSpWspo1rGMt69nAJjaymS1sYyvb2cEudrKbPexjL/s5wCEOcpgjHOMoxznBKU5ymjOc4yznucAlLnKZK1zjKte5wS1ucps73OMu93nAIx7ymCc84ynPecErXvKaN7zjLe/5wCc+8pkvfOMr3/nBL37ymz/842+o/wL9AQhEQAIThGAEJTghCEVIQhOGcIQlPBGIREQiE4VoRCU6MYhFTGITh3jEJT4JSERCEpOEZCQlOSlIRUpSk4Z0pCU9GchERjKThWxkJTs5yEVOcpOHfOQlPwUoREEKU4RiFKU4JShFSUpThnKUpTwVqERFKlOFalSlOjWoRU1qU4d61KU+DWhEQxrThGY0pTktaEVLWtOGdrSlPR3oREc604VudKU7PehFT3
rTh370pT8DGMRABjOEYQxlOCMYxUhGM4ZxjGU8E5jERCYzhWlMZTozmMVMZjOHecxlPgtYxEIWs4RlLGU5K1jFSlazhnWsZT0b2MRGNrOFbWxlOzvYxU52s4d97GU/BzjEQQ5zhGMc5TgnOMVJTnOGc5zlPBe4xEUuc4VrXOU6N7jFTW5zh3vc5T4PeMRDHvOEZzzlOS94xUte84Z3vOU9H/jERz7zhW985Ts/+MVPfvOHf/z9/zJPAAIRkMAEIRhBCU4IQhGS0IQhHGEJTwQiEZHIRCEaUYlODGIRk9jEIR5xiU8CEpGQxCQhGUlJTgpSkZLUpCEdaUlPBjKRkcxkIRtZyU4OcpGT3OQhH3nJTwEKUZDCFKEYRSlOCUpRktKUoRxlKU8FKlGRylShGlWpTg1qUZPa1KEedalPAxrRkMY0oRlNaU4LWtGS1rShHW1pTwc60ZHOdKEbXelOD3rRk970oR996c8ABjGQwQxhGEMZzghGMZLRjGEcYxnPBCYxkclMYRpTmc4MZjGT2cxhHnOZzwIWsZDFLGEZS1nOClaxktWsYR1rWc8GNrGRzWxhG1vZzg52sZPd7GEfe9nPAQ5xkMMc4RhHOc4JTnGS05zhHGc5zwUucZHLXOEaV7nODW5xk9vc4R53uc8DHvGQxzzhGU95zgte8ZLXvOEdb3nPBz7xkc984Rtf+c4PfvGT3/zhH3/D/FfiC0AgAhKYIAQjKMEJQShCEpowhCMs4YlAJCISmShEIyrRiUEsYhKbOMQjLvFJQCISkpgkJCMpyUlBKlKSmjSkIy3pyUAmMpKZLGQjK9nJQS5ykps85CMv+SlAIQpSmCIUoyjFKUEpSlKaMpSjLOWpQCUqUpkqVKMq1alBLWpSmzrUoy71aUAjGtKYJjSjKc1pQSta0po2tKMt7elAJzrSmS50oyvd6UEvetKbPvSjL/0ZwCAGMpghDGMowxnBKEYymjGMYyzjmcAkJjKZKUxjKtO
ZwSxmMps5zGMu81nAIhaymCUsYynLWcEqVrKaNaxjLevZwCY2spktbGMr29nBLnaymz3sYy/7OcAhDnKYIxzjKMc5wSlOcpoznOMs57nAJS5ymStc4yrXucEtbnKbO9zjLvd5wCMe8pgnPOMpz3nBK17ymje84y3v+cAnPvKZL3zjK9/5wS9+8ps//ONv2P+KvAEIREACE4RgBCU4IQhFSEIThnCEJTwRiEREIhOFaEQlOjGIRUxiE4d4xCU+CUhEQhKThGQkJTkpSEVKUpOGdKQlPRnIREYyk4VsZCU7OchFTnKTh3zkJT8FKERBClOEYhSlOCUoRUlKU4ZylKU8FahERSpThWpUpTo1qEVNalOHetSlPg1oREMa04RmNKU5LWhFS1rThna0pT0d6ERHOtOFbnSlOz3oRU9604d+9KU/AxjEQAYzhGEMZTgjGMVIRjOGcYxlPBOYxEQmM4VpTGU6M5jFTGYzh3nMZT4LWMRCFrOEZSxlOStYxUpWs4Z1rGU9G9jERjazhW1sZTs72MVOdrOHfexlPwc4xEEOc4RjHOU4JzjFSU5zhnOc5TwXuMRFLnOFa1zlOje4xU1uc4d73OU+D3jEQx7zhGc85TkveMVLXvOGd7zlPR/4xEc+84VvfOU7P/jFT37zh3/85f8AnCZweQ==
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" format="binary">
+AQAAAAAQAAAAEAAA9AQAAA==eNol1yHXOQgWwGG7+98du2N3BEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQBEEQhDcIgiAIgiAIwj5z5pz7fIBfuvf+CgQCQUKEiRAlRpwESVKkyZAlR54CRUqUqVClRp0GTVq06dClR58BQ0aMmTBlxpwFS1as2bBlx54DR06c+eHClRt3Hjx58ebDl8DfAoFfBAkRJkKUGHESJEmRJkOWHHkKFClRpkKVGnUaNGnRpkOXHn0GDBkxZsKUGXMWLFmxZsOWHXsOHDlx5ocLV27cefDkxZsPXwJ/10+QEGEiRIkRJ0GSFGkyZMmRp0CREmUqVKlRp0GTFm06dOnRZ8CQEWMmTJkxZ8GSFWs2bNmx58CRE2d+uHDlxp0HT168+fAl8A/9BAkRJkKUGHESJEmRJkOWHHkKFClRpkKVGnUaNGnRpkOXHn0GDBkxZsKUGXMWLFmxZsOWHXsOHDlx5ocLV27cefDkxZsPXwK/DEFChIkQJUacBElSpMmQJUeeAkVKlKlQpUadBk1atOnQpUefAUNGjJkwZcacBUtWrNmwZceeA0dOnPnhwpUbdx48efHmw5fAP/UTJESYCFFixEmQJEWaDFly5ClQpESZClVq1GnQpEWbDl169BkwZMSYCVNmzFmwZMWaDVt27Dlw5MSZHy5cuXHnwZMXbz58CfxLP0FChIkQJUacBElSpMmQJUeeAkVKlKlQpUadBk1atOnQpUefAUNGjJkwZcacBUtWrNmwZceeA0dOnPnhwpUbdx48efHmw5fAb/oJEiJMhCgx4iRIkiJNhiw58hQoUqJMhSo16jRo0qJNhy49+gwYMmLMhCkz5ixYsmLNhi079hw4cuLMDxeu3Ljz4MmLNx++v/21+H8RJESYCFFixEmQJEWaDFly5ClQpESZClVq1GnQpEWbDl169BkwZMSYCVNmzFmwZMWaDVt27Dlw5MSZHy5cu
XHnwZMXbz58/zx6/q2fICHCRIgSI06CJCnSZMiSI0+BIiXKVKhSo06DJi3adOjSo8+AISPGTJgyY86CJSvWbNiyY8+BIyfO/HDhyo07D568ePPhS+A/+gkSIkyEKDHiJEiSIk2GLDnyFChSokyFKjXqNGjSok2HLj36DBgyYsyEKTPmLFiyYs2GLTv2HDhy4swPF67cuPPgyYs3H74EftdPkBBhIkSJESdBkhRpMmTJkadAkRJlKlSpUadBkxZtOnTp0WfAkBFjJkyZMWfBkhVrNmzZsefAkRNnfrhw5cadB09evPnw/f2vo/8XQUKEiRAlRpwESVKkyZAlR54CRUqUqVClRp0GTVq06dClR58BQ0aMmTBlxpwFS1as2bBlx54DR06c+eHClRt3Hjx58ebD98+H57/6CRIiTIQoMeIkSJIiTYYsOfIUKFKiTIUqNeo0aNKiTYcuPfoMGDJizIQpM+YsWLJizYYtO/YcOHLizA8Xrty48+DJizcfvgT+p58gIcJEiBIjToIkKdJkyJIjT4EiJcpUqFKjToMmLdp06NKjz4AhI8ZMmDJjzoIlK9Zs2LJjz4EjJ878cOHKjTsPnrx48+FL4A/9BAkRJkKUGHESJEmRJkOWHHkKFClRpkKVGnUaNGnRpkOXHn0GDBkxZsKUGXMWLFmxZsOWHXsOHDlx5ocLV27cefDkxZsP3z/+evj/DxdNFi8=
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="binary">
+AQAAAAAEAAAABAAAEQAAAA==eNrj5BwFo2AUjFQAABo4JAE=
+ </DataArray>
+ </Cells>
+ <PointData Scalars="scalars">
+ <DataArray type="Float64" Name="velocity" NumberOfComponents="3" format="binary">
+AQAAAACAAQAAgAEAEjYAAA==eNqlnXd4FcX79hcCKfSEQAoQJr2fnFBCKAkbQqR3pEhbFFCaKFWRskiRKh2psoD0rhQFclhEOipdUcoqHQSV3nlfnvl+z/U7v+vc73XO++Q/h/Ez93n2ntlnZmdnFQX/fd3ledb/45/FH7ue5v+P/7bGP3vyur41OOlVvgcc8+G/RbLdlCOOGfw71TfuxLryvdXz3/r9nxV2/M//Hh33lHS+t8yl3Fs9iGOsqEv1xdJtLuUoDsqGzsHuymG7z7OpvjH9R9ffCzjWykbl3JUjjvZpZaqv9VMUB0MPioM652Fp0nPBz9fB0AM528xw4lQr5lLurR+cvzusYqBL/b/6kU+UZTmBDoYfIKf1wYDX9dUhHQM5epz1u3Yr63Id448+fl1ffTCmLEcP4qixO/1f17eubizr4PgBtCtadaP64sDkQI4fICeoHNUX4QsCWXpAHHR1X9Drcj3qalmOHsSxTsWHUf24nBCHB+OAHhAV6a4c9ruyRai+2rucSzniKLkbYtyVI445YCbVV6fWdyn3Vg+67ua1E4L0vJ3l62DoQRx9ZI78vbe/cx1/QBzM+fUS3Jajdp8lUH3toxau5Yjz8Fay23LAMW4fpPpi7weexQFwUByMOT/GU7/IfejrYOhBHKXxu1RfzUryY/kB9DuzQA+qbzY5HMjxA+JYdw9HUxwGRgRx9MBx4My2ihSf+iNCOHoQR+sdJuu/cSSE4wcYh8OVqb668/0gjh8QR3/fTCKOz9ogjh4UBzWmfxxdl99qh3L0II515HgixXnlhlBOPuCM6/4n5Vz03Jrw6HV9M7xUeU4+ADmz8/0oDhMjy3P0ONtfYQgXn6yq95Dy/HeXCI4exNHSG/jSdf/zW8HJB1C7otqVQBrHWgwoz8kHEMf6ICiUdP72d3mOHhQH89L1UlT/s18FSw/gaFWXlqX6p/wjOX5wtj8pONpF57xv7pPP36sYzfED5Mw7V
oiu48Ya0Rw9zt/bp22sS/m6u3df1ze+ahPL0gM4yuznBSn+DYbHcvyA2tXOdylJ41v1btEcPyCO+klwGeov1qhojh4UB2XIpmJ0n737TSxHD+IIpYqcR5y8EsvJB1C/s3pUjqDyVfUrcPIBxBHf5EURZ8H4Chw98LpnfVKe4rA6J5KjB3FERCrNI6yWn0Vy8gHUrhrWO5Z8HlUpgpMPII61eDnl4Wb4jxEcPSgOSlQ65YHi7tFIjh7EMbpWpzzKtEVFcfyA+p35TnfKw/UuR6M5fkAckXCB+pEZGB/D0YPGAf1WJPUjbUdqHEcP4ug1y1cgvx0fH8fxA2rX7LWbfGhceD+G4wfEER0WkQ/1Kd/HcPSgOKin4ygORpHf4zh6EEdkHaRxRpmREe/wYB1An70kzV05mndrt4dQfXPtMtf5PuIkbkt3u74EOOaURVTfSrnk0ToA4qB5n5nsS/WtP7b5ORh6IGdiZ6pvOFL9HR7EQV+/pLLb9RnUbvmpVF+vV7GAwxPObqOq23LAEVvHyvo3+xRwMPSgOJjd3qD65o7v/B0cPYAjbBlUX2zMCOD4Ac27zVv1bKQztUhpjh8QR6/+j510po4qzdGD7vv623mpxOkUFMbSAzjqgG5U3+qxIIzjB9hu1MlKdH0bO0pz/IA4atTWKlTud7s0Rw+Kg7E3heqrpR6GcfRAzv5XFActp2K4w4NxQLloVHP7fBD0O+3cSFn/0jiXcsRR7Qurux3/AUf/chTVV3dN8GgcQBx03Y3rCVTf2NYlwMHQgzjGywiqbxwZ5VKO4mB0XVnT7X0BtGud06m+OL7ZVQ/gmHFrstw+FwYc80Nd1n+w3qM4QD0gDnrVsJpyHJ4X4ODoARxFFJP1M0IDOH5A/U5tNiOD+nWtmGCWHwBHLByRSePwhHrBHD1oHBBDf6I4GF++G87RgziG7QsZh7Kzwjl+gHHo3q0GlSfnBHP8gDhKnRa1KM4lTpXm6EFx0M8MoD
iIL0aGc/QgjjmuLsXBvJQXzskH0LzbSP2H8iizeqOKnHwAcbSSzSmPMppdqMjRg+Z9+sftKY8ydo6M4uhBHHP3eIqD2uZUFCcfgO1W70R5lPo8WXDyAcTRq5aRebh/Y8HRg+IgbnaiOOi3w6I5ehBHWfNIzlMmto7m+AGuA2wdJX34U41Yjh8QRx37IfnQ+uSLWI4etA6gjOtF64TqByviOXoQRxsTROOMdqpYAscPqF3tl+40TxF9L8Ry/IA4+gfrJWeoiOPoQXFQp+8kjta4awJLD+CoZfra5Di8OoGTD8B+l3qA8hBt5ieCkw8gjrVjJOUhWp9xgqMHXXexYppcJ2nyWTRHD+LozeIoDpp9djQnH0Dtmv6ZlIcYRYcJTj4AOYXuE0fz6y1YekAcxLD9xDHyPonm6EEcdVMn4pitm0Rz/ID6nSifIO+z99vEcfwAOc9XEEcs7BnH0YPGAaPKbMkpdTmBowdx1Oy/Kc7i2cMEjh9Qu1rZSDmPG984juMHxDEqzyGOCLwTy9GD4qBcqi856vkEjh7EUV/OpnHMXNkngbMfwPm77/R3ub8bQ57/Q/sEavaJ5+wHQBx9QHAB8uHHo+NZev77d7tFokv9e0F3aJ/P+taJHD2II7qXffm6XBTvlsjZD4Da1R9/Q/uQzbL74jn7ARDHHJVJ+3DUkFIJHD0oDsbAxbQPxDq0IJGjB3KeVylB9TOvJnL84OyXz60kl/7V5+CN1/W1b28mcfyAOGLEN7QPX1tSPJmjxxm/wbYUF52OslfI5z/XTOHogZzv79M+ImNyhxSOH1C72oZKPjTOx2Qmc/yAOIryrfThgoHJHD0oDtZHnV8RZ+DEFI4exLEmX6J9SmbuTyms/QCg3ynjltI+BHVSuwTOfgDE0UpWoX0IVjsjgaUHXHcrT6E4mG0ykjh6EMco8KysHIenJ3H2A6B2lfYnaJ+SXq1gImc/AOSs+oE42vcDEjl6UByUN8vQexzq6id
JHD2IY0XY5D6Nd3sls/yA+t3ZiqWkD08nc/wAOZU+k+8x9ayTwtIDxgFrlHwPS69WJZWjB3GMDsdoH6PqtyqV4wfUrnrUJ4R8GLgyheMHxNHbjySfa51KpHL0oDiIvdton7bydSkbRw/iiB0jaBxT/hpi4+QDzn8vU8P1Oo5+dZHu4xWyUjn5AOIoq/0pDxe9W6Zy9Djj0fGRS7l+VD9N+7SjX6Ry9CCOfm8KxcF6s7iNkw+gdrV6f9N7CtboKamcfAByWt8gjthxPZWjB8ahbCXiaAl5No4exFH2zH9APvxhpY3jB6cvjr3rykn6ZR/p+b6vjeMHxLH+LXmGfteXY2wcPc64Dp7hUm76Xh/0ur4+eB5LD+JYDbM/p3Jzk43jB9SufnuERfx1O2wcPyCOFmu/SuXDg9I4elAcrPD8eTRPLHWVpwdwFDV1Pl3HCzXSWPkA6nftChFHzGlm4+QDiGMWakXvSYm2220cPei6i4rpz6jcNziNowdxjJ29FcrDfxmdxskHULt63kS5zrMmPo2TDyCOuqYd5YHm/blpHD0wDtWzaZ6rTH2SxtGDOGahGLkOUDPDzvED6nfmolW3aBybPT6N4wfEMZ6uleuQDwvYOXrQOKBGLzSI/3gzSw/iqOq+r0inLczO8QNqV4/8/h7lUQdr2jl+QBxrxV16j08d9qGdowfFQRs5eD3VV1qx9CCOeW7tZrqOdSbaOfsB0LxbbVWB3rMw0v5I5OwHgBy1n+RMaJjE0QPXAbotkudUdL+QzNGDONY9g+JsHmmYwtkPAOOQdUSekzB1YxJnPwDiWJ1rSs7je0kcPSgOSsg3tI9Fvb0phaUHcMTcY3KfT7sHKRw/oHm3lnOG1vGsrp+nsvwAOFqDr+gcBrWbj42jB933RcsatE4rPn9h4+hBHOVioDxn5mjfNI4fULta8b7yPanTb9s4fkAcc5Y9Rr5nsd7G0gPioDZ/IdeB4/ancfQgjnnlHdKppxW3
c/YDoH6n+q+lfVx6hD2Zsx8Acays4sRR6qnJHD3ouisHb8t9sD/aUzl6EEfs9U0hn4zKSeXsB0Dtikt15Hs0R0OTOfsBICe1j9wHYq5I4uhBcbCiksmH5nH/VJYewDH8Q0mnODUsheMH1O+UlzflOHPyso3jB8SxVm2R59hUvmHj6IHjwLkDdJ/SF0faOXoQR8tLpnHe6Bpl5/gBtavv+pjibLzvOo/z1g+I838TdhpnjG8L2jh6UBzM+bPkOVqNXqVx9CCOUu4OXUfRsVkaJx9A826rx0BaZzAvu/rK63wAcZ6upHMghU+enaMHzfvUwSXoOZoSPsDO0gM4xuhzRan+H5vtnHwAttvxuDzHYGV/OysfAByt/026juqpL+wcPSgO5h/H6DmsMu8cSw/iGNuOkE/0JCWd5Qcw77Z2r5P7QJK2sPyAONpHw+kcXb3mPZYetA6gtn2whepPOczSAznbK2+n8lz/dI4fULvGC4PWgVURkc7xA+IYa4sXkOuQGekcPSgO+tmSVK4UqsTSgzi6cmYNrQv9kZfOyQdQvzOyH8tzQoZ9befkA4ijvVuAxjE9bx1LD7ru2qpSQXIfS0g6Sw/gqMP+c45r5+B0Tj6A2lX2FaV1DGPEBDsnH0Acbcw08olIttk5elAc1E/3yXOqW1+3c/QgjjU+tKy8jhN5fgD9zuxblPbxaocap3P8gDjiSG9axzYfNGbpQeOA0lP/hJ4LfNaWpQdxrEMzJ70+v9tSOqZz/IDa1QLz6Tmdcj42neMHxDEX3y5M9VMOsPSgOBh75gx7XW4szEjn6EEcfdzbC1+Xa13Pu44P4BwAremC2m7PDQDv3WsnP6D6+m/1XcsBR7u9R3VbDjjq22upvrZ0vuJg6IHnAOx7mk315zX2dzD0QM7pHlTfGN7tf50v5D4Oypz1OW7PzQDtGkGzZP1/xrme/w84ytDP6rg9dwJwzEO9qb441PtFPkMPioNavBHVV6f8UdjB0
QM45rAIWf/lNB+OH9B799bCOKqv9W9ZmuMHxNH73JL9KGNzEEcPPAegty/FQR0/K4yjB3FUdStxxNwJoQ6GH1C7aoefiSM2LA50MPwAOQsWEke8dbckSw+KQ6VUOQ4khoVw9CCOcv4alVsf33Y9RwKMA1Zi41y34z/od/q5FKpvOT5y/Q4L4Ciht92WI47x9CcqN3+LfJTP0IOuux70ax25zm+6nl/npR7EUSv3l5yamS4cGIcP+9R1Oz6jOIQ1lfXHZt3N94AjAve4LUccMecLKtfD513LZ+hBcdCPXcuV68BZ/+Yz9CCOUjSK6lv1HCc4foD9rkwDWf/VzaIOhh8QRyl2msq1ERN9OXrQOCCuTZLj5MWgIAdDD+KIoiVlnBcYxRwMP6B29ZBc6cOx2kuOHxBHXHtf9tPLM65y9KA4WEPeoXJz2rqCDoYexBH7RsjyE02v5DPyAfTevf50HnG0swUFJx+AnKPViGP16hvB0QPf+356S55jFn00iqMHcdRf3pHnUJ39J5KTD6B2rW1bJKfp2vKcfAByKvhSfbVV4XCOHhQH8+oeymP16d8Jjh7EUUa/yJLfkRERHD+g9+7N0U/pHA/1XqdYjh8gZ3ML4ui9+8dw9KDn/uoPx+Q5DOqOeI4exNGXF5DnKA5aFMfxA253Jp2TZnR5FsXyA+AYDdeRTmvLM5YeFAfl29hM+Z7snRiOHsRRJoSSTuvB6ihWPoD63ewFVK7vyS3LyQcQx9pQgjj6NC2Qowddd2VSEZrHmc33hHH0II51/ACVa3vOlXYw8gHY7tGWVK7FFvR1MPIByKnXRpZ3fe9yPkcPiINqlaE4m9GF/B0MPYhjJEqfGEFbf89n+AH1O3XRNBpnjEV/lef4AXHMa2kyj/roURmWHjAOiIFnZR5VqkNFjh7EsT5oT+OY+VWTEI4fULtmxgdULh6WCnAw/IA4+sEqVG6MKbY7n6EHxmHxECpXz9cNcDD0II6SHS3vI5
tW/ZzlwTqA8ctR9+sGYN6thn5ZV74n9ebJfA84+vVBbssRRysp1wHM4uN/yWLoQfM+I6GiLA+/djWLoQdy6jjkPDSphpLtQRyEst/tuhCMQ+hnsv7t23ezPODoeaPcrvshjpZTmeqriwYVyGboQXFQdl2W/BvVimcz9CCOsbI33S+Mk++GZDP8gObdRkozOU+vuvwWxw+IYx08K+/Xlab6ZjP0oPu+frWR9EP/sw+zGHoQR8nYKON/vn2pbIYfYLuPLsl8YN+gMhw/II5qbybz5F6TI1l6UBz+XCPn6S0yIzh6EMco5i+f44xLdPkeLBoH1I7T3T+fBf3OLF9Z1r/sGmfEMd7KdX8uPeCIKxnyvtk81eV6easHXXft2TiKj/plfmQ2Rw/iZJyV876P+iVlexAHcbp+JbfloF3j1lD5HYqOZ6KyPeCgP8Rx/vuno6pkM/SgOOi9J8j9zDcWVMlm6EEcp95ea3I4fkD9Tinzt1zHWP9zEscPiKMXkudAKoseVOboQeOAkZkuv0MxSK3E0YM45s/D5fsadfbW4PgBtWu+VY324etr0rI5fkAc53h/Oq8eRw+Kg2h+Re4DbDMyl6MHcZx+9I9qnM3IB9C8W8svKstTJr/KYuQDiGO9945chxHvudwHvdWD5n2a1oHmodb3pwpz9CCOPvGmnK+dq16Okw/AOFyZRP4wnkyI5uQDiGPes8txoMw6G0cPjIO6VH7n5UIwSw/iaNMayO907F7qMi556wc071bO9Kd1EqXgpWIsPwCO0s5X5glNFgmOHrgO0KeHXMdudDGQowdxxP2H8jlCoa9YfkDtWrX3y+c4c8ancPyAONrNiTTO64ePVuXoQXEwg0/Lc56bvm/n6EEc44dz9Du1YiszOfkA6nei0cfyO0H/HM7g5AOQM7myzAPL9Fc5etB114ZdqCz3mVerydGDOOqoHnT9zL2ZuZx8ALVrtQ2Q52A0q8fKBxDH+e8fhzXn6EFxMGZ
vke+pZU1pyNIDOM763U+04PgB9Tv9x17yOwKHO9Tm+AFxlLw98jsso7blcfSgcUCczpX9Z9aEHI4eyMnIovdNzEW3WX6A7fbdQt8x17a2aMLxA+I4//1ztRVHD4qDOH2d8ic9Znszjh7Ecfql/9DW2Yz9AOi9e723Ls8ZeONYImc/AOJoa6dUke9BD03g6EHvfVtZRaXP35mdzNGDOFqZG1Su30tM4uwHQO2KRmNknC8viePsB0Aco0hb+Z0XvyCWHhQHs9JOeR1/CEvg6IGc+v3kd44m5MZy/IDeu1cSDqTI78iXSeX4AXG0JXXld6Bm8vTAcwBS5HfQrC2NbBw9iGM9+oDGSTFwQwrHD6hds8ki+g6dMF2/4+CtHxBHLOwg8+kJsXEcPSgOyowxdM6DOTcniaMHcZTubWl9TUu6GMfZDwD7Xei/VG5mLROc/QCIYz7PlfOmiWooSw+47mrcVZmnpQ6M5OhBHHGgIo2TunI+lLUfALRrzW8hvyPm/4+/Jxy0DoA4WvPnMg5+uVeyGHpQHMwGFeR3bPdf9nNw9ACO8ccBeR/vt/LfLIYfUL8zD38jy0f4RXH8gDiizwUq1yZVC+PoQeOA2uhPKtdP943i6EEctXii/I5qOS2M4wcYh87H5X38wEBfB8cPgGNtHSy/05oz5XEWQw+KgzL9ifxO6JJ+hR0MPYhjPawvr8u5HJd1bG/zAfTevZU5k+bp5otqaZx8AHG0AqPp+Z+W/nEqRw9679taMp6e/+nv307j6EEc0dJB51Qo60JtnHwAtvtbQ3mOYoX/dR6ml/kA4uibr0bLfdRd4jl6UByMb9fTPFff/GEyRw/kBH1C53iYo7+O5/gBvXdvWY39qHxIkp3lB8R5eyWds21eHWzj6EHnAKhBG1a9LlfnZts5ehBHy3uw9nW56DDJxvIDaFcc/YzOWdI2rE/m+AFx9OvPqdzKPBfP0QPjMGb5ciof+GMyRw/iaO13fkfPh+YX
TuDkA6jfWc0byHNmHiyL4uQDiKMX6krl5tGuYRw96Lob/or8TpDfsSiOHsQRhYtQuRbTMIyTD8B2L9aSOpdXLuRg5AOIIzbtls8RrBoFsxl6YByi/5TX5dHdgg6GHsRRsmrKc5JrNSyUzfAD6nfGsrZF5Drhv1EcP0COtbyIPGc1NYyjB40DZu2JDhq3fR9HcfQgjrL58rbX5WaCa32v/QDaVd8ZSuXWnwkFHQw/QE7LKXIcXtbK5Tm+t3pQHJTuq7bTuT1v9i/oYOhBHNHv1rfU7tvbCnP2A8B5d3Q1+Zy0XNcynP0AiKP8asl53KoP4jh60LxPzyou5+nXvw7l6EEcpfePtM6p/n4/gbMfALY766p8rtekbyXOfgDE0f9qLP//ro1qcvSgOOi1S8l5dEq1Khw9iKN9uorus9a/67I4fkDzbqvXArlOK74ux/ED4hh/D6ZzyNVNeckcPfC+/0U1OU9MXxzB0YM42thBtA6gFWmayvEDatdsM0XuH216oirHD4hjDf+e/G3W31WbowfFQfweSeswmplejaMHck4+pnUM9eXXKms/AOh31uxI0qmuiczl7AdAHG1zFVpnMHcebMDRg667lV+S4qkPmlyXowdxzPMbqT9YN7s04uwHQO2K9aPonEmjw9XmnP0AiOP897sz3uToQXFQb0yl76ApHx9uwdGDOM76q0e2YfkB9DtlRAn5fsLkvXkcPyCO2ao7rWOoQ83GHD1oHNCM9uVJz1jzDY4exLHst0rL7yidbcLxA2x31GDKk9Vin7Zk+QFwnLyx1dty9KA4KL4D5HcGwwNbcfQgjrN+3QNtOfkAmnebN86TT8SuroKTDyCO+sMs6duRD2wcPXAd4EBBuc/qfkIkRw/iWBvby/eLxvnZOfkAalf/c4w8J/mSmsnJByBn+Cp5DnbAmzkcPSgO1qQLwfI63svk6EEcK2swzT/Fj/tyOH5A8271ulwHFve/jOT4AXHMKum0Dq9md
LZz9KB1AKuYH63bW+fPRHL0II7V/9EyarfzcDvHD6hdpXwmXT91V7vqLD8AjjI0g94vVbcUrcPSA+Kg9vxwM5WHTqvO0YM4Wue/6LmAUjylDicfQP1Ob/tSnrO9YlQ9Tj6AOFZAtnz/dv/XTTl60HXX666k71AYx/3rc/RAzslldH8yIho24+QDqF3zjTUv6Lse80e34uQDkPPf3zGzejuOHhQH/cBgWn8xtp9pxdGDOE69ETPbcfwAx4HwCvSdXz24fX2OHxDHvBdPv0sduqoZRw8aB5Q3M3fQc8OWn9Rn6UEcq2E+PRdI+LkZyw+gXTG5En0nXf2zdGuOHxDH2X76mXYcPTAOvZvOp/IfU1pz9CCO82/8nXZu1wEGz3DdV+B7fRBd98HzXPcD//ffHw9y3RcU+Hjv63b1UUM84lgNsz+ncnOTZ5z4kFPk2zETWXqc496Owq77Dd6ZSnxtRwBLD+KYzWbS+8nmryEexcEKz5/3utwsddWjds0Y/QJd94wfPOKg74AjjlZ5AekXLyqmcfSgOKjVU/+mvLFBCxtHD+IYh/68T+0O/IblB2c/K9HM9buQzQpfJD13mqZy/IA4ZnoQ/S5RumMqR49zXLnk+r1dMWO81DOlSQpHD+Ioh/vL+J/uksLyA2rX3oTOH9f/nJPK8QPiGHnV5Pcnaj9K5eiBcfC/Tuc760Vnp7D0AI46cqPcP7H1Voon44AavdCg+87jzWke9bsBP92g31VrhkccVd33FY17tjC7JxyzyN47NA7UKGzn6EHX3eoXSOeP629Fp3H0II7WsSVdL7PTFx7FQRs5eD3lLUoruyftiszydD6+eaC2ZxzwHUDEUfyKPqB9BbOG8vSAOCj1xtJ+BXO0L0sP4ugZTWn/k1Yz187xA+p31tYp9P0A/eW7No4fEMdcPIX2r5urfrNx9KBxQCv/N70forXpnMrRgzjii3D6PoF55GQqxw+oXaXFRPo+gVq1SRrHD4ij3v
qG9KslzDSOHhQH60htOh9Zv55j4+hBHDPsFZ3vLDI32zj5gPPvvI/L/lK1d00aJ5UDAcmcfAByajZ4QuPMswrJHD3OOPl/7Po+Xek/bhP/q48SOXoQx+h4XX6f5tPPEln5AGo3LpD6qbq0aTInH0Actd5X1K+tMouSOXpQHCzbDPl+ywAzkaMHcbTadUvI8rJJHD844xq6zvW9+Kxw+r6IEbY6nuUHwDHMdPn9uWP/q10v9Tjj0Wa263uvjWPu0TzliGu5t3oQx9gfLr+vM2hpLMcPqF1923F534m8Gs/xA+KIbk3pfHPrYd0Ejh4UB3HxOJ3/rlX/KZajB3GUyk3ofG29WUgcJx9A/U5vUbQUlTcqnsLJBxBHzFNJv7ZkbApLD7ju2skn8r7ZalgSRw/kNDhB3w/QK1xN4uQDqF2l2+fy+wSb/kph5QOAYx6W31cw3nVdZ/BWD4qD0XIX8UUp13HeWz2IoxXZRef76122JXP8APvd3J10/rtuLU1g+QFw9J+qyu9bvFASOXrQOKBtr0Pv36pfDI3j6EEc89+cCsR58Hscxw+oXWPgBrq+YnXPRI4fEEcd8iW9v202PpnI0YPioCbZia/Vzojn6EEcNXkPfd/IyJ0T78k6gNr2wRb6Tv2Uwx7Nu9WxfeQ6VYzDI461vTJ9p17L9U/3hKPdPE15pv7qPksPmvdph57QPE48Gm3n6EEccW59EbpeZXZ5FAf9bEkqVwpVSvekXXVYKcqj9NshHnHQd8ARR3zUjr5fqDniWHpQHMwPlxUnf17+1c7RAzmLutF9XE+9xvMDmHcrs17QOKyXjGL5AXGMYWF0v1YPtWfpQfd946tjdP6IObxqGkcP4lgn59E4ae1ancbxA2pXZNlknjBmsJ3jB8SxZneT51OMHcXSg+Kgzf2Y7gvqLSuNowdxrFbJdH6NGng3zZNxQOmpf0LPBT5r69k4MKgUrbOJltkecaxDMyfRcxmlo0ccY8pBmje
JynksPXAcyPGlPM16+cLO0YM41tWZlHep/gEexcHYM2cY3dcWZnjUrvV2o8IyT4vyiIPOAUAcsXMpzaPVNa73cW/1oDgo3YtSnqntvcbSgzhm26nUj0TGZDvHD7DfKTupX6v3v7Bz/IA4SruX9P0tcdZg6UHjgDV3KeWNeoMSdpYewBFb4ij/1AeG2zl+QO0KezU6X0y9O47lB8SxhjgobtaWdJYeFAdl+qoY6ZMCPD2I4xNJ8w0rskMaJx9A8279jYF03xG7z6Ry8gHEURd9IH9vgaY2jh64DhDbk3xrPElJ4ehBHHPFEDn/erkwhZMPoHaVxEx6314/q9s4+QDiqN39yVdmiZE2jh4Yh/6TaL1SL3YihaMHcsYbNO+1fC6kcPyA5t3iw3/od1kb1CSOHxBHr9+G+q94Oj+JowetA1jtWxHfWPoynqMHcdTgv2n9yDrSMYHjB9hu0DfEFYfOJnH8ADl/lJX+GHgviaMHxUFV1xBf2zQ9gaMHcdT7XaifmE9XJnDyAdjvth6iOGufb7Ox8gHAMacsiJfrhKdsHD3wuhfYTf3HHPkshaMHcUztEl0XMbhCKicfQO2KiqOpXevEXBsnH0Aca1Aj6av0J6kcPSgORmhRqqc3fZDC0YM4uu0/8f3rnRSWH0C/s74aQTqsK8WSOX5AHPPyfbpOesVSyRw9aBzQsmR+qXf6IYGjB3H0oPOU31sDjydw/IDaNdtUJr72VEnm+AFyDvWsRPprz0ni6EFxUKOyK1McGm5P4OhBHH3WjCrE+VVLYO0H+G+/XNMq2oW/8DztY9EPvBXN2Q+AOGrMrUJyHbVfNEeP88/nost5ztpnnR7SPg3dtdxbPYijHRwg9yFM/1tw9gOgdtXUvrQOrB+bGs3aDwA4xpGKNI82HzqiOXpQHLSh9+Q6zPZykRw9iGN8spHm+8bjppEcPzh9UXG2y/fQjRtb6DmRdWVmeZYfAEffe5aeI+jlp5Xn6HH6wu9W
WZf63U8/pn04+8+V5ehBHLXxCTq/W315ryzHD6hdc+y/9JxClLxRnuMHxBG1Sst1Nq1nBY4eFAdR/hCtE4rDqSEcPYijzvrPeXPPBoVw9gOgfidSPpPPwUuHxXD2AyCO/rMP+UMfOyqGowddd3HzM4qvfndKJEcP4iiDash1gH1nIzn7AVC76oSTtE4i3jwWw9kPADm1v6P7h7koKJajB8XBCG8o53EbYqM4ehDH3KbR71QX9I9i+QH0O0PLoHUkK+hoBY4fEEfZ2lDqLp4QwdGDxgFd+5rWmUW3P0M4ehBHtIsivwk9L5TjBxiHjLHkA6vergiOHxBH+ehnypv1JR0rcvSgOJi1Jsr5+OFVoRw9iGNmvZDrg/X/17mjXuYDTh81+SDQReeE2bQPxNqxNJCTDyCO2sCU+zyn5wdy9Dh1Df7F5TxktU4B+Z7OkKmFOXoQxyqdRfsclDzXcm/zARiHze/J52jzjwSy8gHAMaaG0H3QnH81kKMHxcG8+YjqW8Ncr4u3ehDHurld7gNseq0wxw/O+m+Uf+HyveBiM6RPWjwqwPED4pjNt8l9MrvuF+DoQX/mxT3P892tC3ipB3HE4wtF3L534aUfYLu/VKbnmMaSjIIcPyCOoaSUo/Idgwpy9KA4aPXqBbt9r8NLPYhj6I3KuR03vc0HQL9TwrpRHqjdzgli5QOAY8Sdpvrmls1BHD3wukcep/pGxTxfjh7EUUu2ozxEpP/uy8kHULviYluZxw99EcTJBxBH7/oPrTOqGXVKc/SgOOi77lF9M7mpH0cP4mjF5iTLfZJb/Dh+QP1OHAim/NXaNKEgyw+Ao1RZSfX1GcsLcvTAceD3dLfrYt7qQRyrrRnjtl0v/YDaVbIbUz6sLnc9599bPyCOtaeQfE7wRgsfjh4UB2PVe27X5b3VAznjE9y+x+vtfgA079Yuj5XPqavqsZz9AIhjFG1P10vfejWWowfN+9TTg2k80S7uj2LpARz11
nh6vqXGRURz9gOgds2FjeVzvRHF4jj7ARBHeW+jfF620jeOpQfEwarwFvGtErnRHD2Io4Q8pedo6keNozl+QPNuUaYc9S/xxbWKHD8gjijTn54ziv/s1/j/1YPu+0b9kfI579hlYRw9iKPmf0rXRV1fJpzjB9SuWbaX3D96sJVg+QFw1Fvl6bmq0aW74OiBcXjcgJ7nGuPVcI4exBHbg+j5r6rXC+fsB4D9rnU8taufrxLH2Q+AOFaXL+Tvvd45jqMHXXd1zDz53PxU92iOHsRRhiVUJc6s6dGc/QCoXWt6KeJrh+rFcfYDII6+fkIGtbv9YixHD4qD9vK7DHmOyqxojh7E0S41zaQ4j24QzfED6nf682NV5D7wTwXLD4CjFBhDv8sqPU1w9KBxQA/5hfSIuZ+Hc/QgjlbcqEY63zoZzvEDalesT6frK8YuFRw/II6x7WJ14gSMEiw9IA7arr7V5TrPznCOHsQxIuw16b/f7BnOyQfQvFsbKM9/tH7aUJqTDyCONrA89VNjQlAwRw+a9xm2RnQfN1tm+XP0II4ydbwcZ5Y5/Dn5AGx3x2O5H2x+ejAnH4CcwF/k/aJTgWCOHhiHoQNpHLP2Hvfn6EEca8VbclwasMSf4wc071av5qfJfQ4TfFh+ABx1/n6ZBz5f5MPRA9cB9t9Mc1vupR7IOXI23e26hJd+QO1qp7dQfbPhZh+WHxCn3Uqqb5b9zoelB8RBbDlY2e1zAS/1II4126zqtl0v8wHU75Q+q+k+qBbICmblA4BjrBlB90FjkRHM0YOuu/p7TdJjTlvuz9IDOGKAoPuyuqFCACcfQO1aJ7vWIP7gScGcfABxjPNqLfpdf/oFc/SgOKjNS5IedXVCAEcP4uj2V1TfmjrYn+MH1O/00M+pvri91YfjB8Qxn3xC9YXfBB+OHjQOGD+sreZ2X6mXehDHSppb3e3446UfYBxaD6T6RoThw/ED4hjr3s+icebvQz4cPSgOxt
4varq9b3qpB3GMV3Oy3F5fcA6AGrSBvuOjzs22e/LevXq9LsVR/JrsEUfLe7CWzm3uMMnmCUf/dQVdP1HkfRtHD3rvW3vxgp5Pi8X30jh6EEfbV5fmP2JTnM2TOIgxy+k7HWLgj8metKv2ctC+AmPiN8mecNB3wCGnzDw6h9ZIux/P0YPioJ7vLs8Z6PNeMkcP4hhf16Xnnpq2IZ7jB/TevfntOfmeY+VqaRw/II5+azc9P9DqTE7l6EHv/al/TpJxuaWy9CCOYR6Vz9FKu77n660fULv64PHyPcr0+0kcPyCONaWmfM77qkM8Rw+MQ6fTtA5j9KqcxNGDOEr+d/TfIvpSnCfjgFl7ooPGbd/HUR6NA6mV6D0jy/F3lCccZfPlbXTucYLruhPi6CmzqFw8Twnj6EHXXRlpUrl18Pcojh7IUZ6Xk9e9U5gncVC6r9pO5/a82b+gJ+1adZvL97xGTi3oCQd+BxBw1PwUirPo29ulvrd6UByM+j/K61K5qo+DoQdxjFWKLO+X4/KdR2/9gPqdNVu+Z62sXhPF8QPkLN9J5erod8JYesA4oN5pJ9/DPdMriqMHcpImyvJnuWEcP6B2lf1fyjiUqlDIwfAD5HwZJd+jD41z+c6st3pQHKxRch+FMn5hYQdDD+JojRx03fU3V77MYuQD6L17Jewp3c+sKYVY+QDiKGI+/f/G9KHJHD3ovW+jVJhch4wanszRgzhqpg+tY6tKMCsfQO3qiS/ke+j9qyZy8gHEURV5Pqz+pS2OowfFwSz6p3zOvrNgAkcP4ihhM+U66uaUWI4f0Hv32u6pxFfbbUnk+AFxtNgVtM6p1+6ewNGDnvvr13+ldSF13OZ4jh7EsX4tS+tCWuEpcRw/oHbNj5fIdbCNs1h+QByR+AmtAysnAmI4elAczKNVaB3VPHkzhqMHcfQrSXKd6sBiVj6A+p26/rrsR5EBUZx8AHFE29LUj4Sva57prR543R/503U0Ow2
N5OiBnK51qVy5fzOUkw+gdvVjBdLl7zrh62DkA5AzbDnFQctv8jiLoQfFQQ2tJcelSw/8HAw9iGMdviDHsTul/s1i+AH1O/F7qHyO82KW4PgBcbRiPWT5FlsoRw8aB8xLt+V+ko8TK3L0II62sDM9x1F2RIRw/ADbHfu+HFdHRQQ4GH5AHOVUkCzvU+NSFkMPioNZZwSVW/vfD3Aw9EDOvmpUbvb4eE+WJ+sAxfxo3d46fyYy25N1gEWvqD2z0JZITzhW/0fLqN3Ow+2ecMxDheQ6/NNOdo4eNO8TX+2m8csIjonk6EEcUfoHut5GUBG7J3FQe364mcpDp1X3pF1FD6HxXS3Yu7onHPgdcMARtq60/qJvC6nD0QPjcM5fnuMR/iSTowdx9MS79H6gevNgDscPaN6t1v7POX69ugiOHyDn3nZ5ToLvXRtHD7rvq/HTKb8Xv62J4OhBHK2OnA+YLd5L5fgBtWteWUf3G6tWRibHD5Cz8zg9R9DaN8/h6EFxMAbVpvxD1FercfRAztpH8pzYdbtUj8aBNzN30HPDlp/U96TfaR1n0HcGzaHd63vCMa2G+fRcIOHnZp5wRJcl8hy28N3NOHrgdT98sZg8B6ZsfY4exFFs++R3VBf0bOZJHPTeTedT+Y8prT1pV38Wdonusw/Lt/aE4/wbf6edJxxnPDqcb8fRA+Nw7SqdE6ivvtmKowdxnONT7vJ2HD+gfqe1iJTjTOep9Th+QBxteHP6jolS5nRTlh4wDojFyXROkdr8lzc4eiDn0x/k/uk2r5pw/IDaNUKu0vkXxph5rTh+QBynX863bsfRg+IgGvWk9X3j/dhWHD2I4+QN/L0tJx9A827z9i25TisiynPyAcQxv9xO6+/6iqXJHD1o3qdNDJPl95LCOHoQx9p4Qq5jtHwvkZUPgHbFitXyPdMeBTM4+QDiCPtuec55o/u1OXpQHMTaQPkeRK1xVTh6EMfo3k8+Jxrsm83xA5p3
K1G5tA6jdZhbhuMHxFGCr8j3sNLOxnH0oHUAY1tvWrex2s8N5OhBHP30PdJpXd0fzfEDalf0k+v8avTOShw/II5yXe5fNJqursnRg+Kg19kv34PbvdrO0YM4Rq1z8jlI4/uZnHwA9Ts964I8Z7LxzTxOPoA46oeXaH3NiLremKMHXXdj4A80TqrnHXU5eiCnVVN5jmLk+EacfAC1a9ZKoX2S2h/TW3LyAcRxth+d15ajB8VBr1iDznm2lp5pwdGDOE7ekIltOH5A/c6M8pX3qe+b5nL8gDiiQViUXIf/twFLDxgHRKW35HP8pt/kcPQgjlibK7/DUr9MfY4fULtaWCf5nbWFT5pz/IA4Tt7TpW9y9KA4GFcv0DxLjPqlGUcP4jj/feYE1/UBL/cDoPfu9fi/5Xsidd6K5ewHQBzNqE95hX5jaAxHD3rvW5v0O/G1vH1RHD2II0Ibyfcy7haN4uwHQO2qKZPlc6VP/aI5+wEQR09cRvp11T+SowfFQftwpXyf6+s7gqUHcKxlN7LkOaLjIjh+QO/dG1vHkh4zs4Tg+AFyrpaV7/sU2BjB0YPe+9NSL1HcxF83wjh6EEdUnJZNnAP3Qzl+gHH4drLkry9QgeMHxFHG3KxN5YMHh3P0oDhYhUJU0uPTK4SjB3KmH6VyY+i7ZTj7AVC/E59Pojgbe46U5+wHQBx9SSzF2RryXRmOHnTd1cj7FGerTaFwjh7EUYdsltfr2D+lOfsBULvmkO5UbjSMC/BoXwFYB0AcLbBKlnzPa+ZP+Qw9KA7GPwHSnyNa+zsYehBHneIjy/1u/5HP8APsd08mE996d39Zjh8QR7/yjMqV05sCOXrgOHC4dw7FocE3QRw9iKO2/pfK9YDrxTl+QO0abVTiWwHbfB0MPyCOvq8RlZs7/7qRz9CD4mD0aFSH6sft9HEw9CCOcvQ9WX57/c18Tj4A3rsX1wqTz9XvfyjNygcQ55EpfeLbpjRHD3rvW+9hy
fr6JT+OHsSxXrUkn5h+Xf04+QBqV6zeKfnDooM4+QDiqFG69O2t90qx9IA4aHvTpZ/PXSnM0oM4cwqT/43ouT4cP6D37tVWb1EcrLMtfTh+QBxtxVxZf8zwghw96BwAkTCyttvzB7zUgzii2ybV7bkQXvoBtjt2GNXXWhcryPED5AS3ovpaarsCHD0oDvqzuTlu4+alHsQxVvWp47bcy3wA9jt7ktRTO7M4Jx9AHDVkZx05Dyroz9GDrrsas1vGJyWyAEsP4Ci32+ZS/WfjnuUz8gHUrtEjSfLNgYqDkw8AjrmtPZXrrcfdymfogXHY8ZPkH9n1L0sP4GiDAuvS/x8w+GQ+xw+g35lflZY6n0x+lc/wA+JYji1ULi786lLfWz1oHNDybbluvz/irR7AUUadcV/upR9Qu2Z8Zbq+hvL2I44fEEd7NLKufM+lAEsPioOqtKnrNv5e6kEc9cVX7vle7gdA827lcn/Ko9TL14ty9gMgjrnrkVzH6LZAcPSgeZ9Z+22Zx47LL8TRgzhm7EuZd6VfDufsB0DtGvN3yHNi7cdTOPsBEEdkfi7PDWvztCpHD4zDP9vlOVedKiVw9CCOfqq95LzxUyWOH9C8W/uwhJynp33/PIvhB8g51FPO45YfDeboQfd9bUUdeR8xx9zLYuhBHP2tr2UeWLNmSY4fULvW14tpnmVMHhPN8QPiKN/WqinXYV7YWHpAHJTE9XJe0Lh4BEcP4hizykhOQK0Ezn4A1O+MZn3lOfZ9ltbm7AdAHGvzd/I7OPV+y+PoQdddGX1ZntP+4v2aHD2Iow7qLc+RSGucy9kPgNrVqu2l9TXrZK8mnP0AiOPUe6ZRK44eFActZRY93xK15jbk6EEc51/ery04fkD9zvpwhjwHfuPjDI4fEEf/siHt8zGLjFI5etA4IFJz5XP8Q60qcfRAzoixcpxxXK3B8QNs94qv/E5ZvRb1OH5AHKfeI+HNOXpQHLTffpb7D05+ns
vRgzjOv3diXPbpeZsPwHn3+rpUrraMupLFyAcQx6p/gjjibrLLe75e6wHzPmtdSSrXHmxxqe+tHsQxbdty5XtAIa+yOPkAaFdff0PmaW1PBnPyAcjp3UXmOd2qR3L0oDiIJbfkOtLhp8U4ehBH2zBS5smTc0KyGX5A826z6CQqVz/1vZTP8APiqMVzZHnG4qkcPWgdQF26zf36gJd6ICeln/tyL/2A2rWCx8l1mAErb2Yx/IA41qlMOY7VFa7vSXmpB8VB1M13uy7krR7E0fYNd7/+72U+AMeBxv5yn8za3UmcfABxtNryO0HqqpJVOHrgdV+8SMah9i+RHD2Io529I/eRfudI4uQDqF19TSvaXyt+7pXNyQcQx+nfmz3qsfSAOGjvLJDfkaxnVWHpARyn30f/lMPyA+h3+oUa8j5b71URjh8QR4xuSetsRlJ4WY4eNA7o6iz3zyW91AM5RZu6PZfeWz+gdsXLafI9r649ojh+QBxnvI9+ytKD4qBWy6nk9r7gpR7Ecfv3fwCkoPpk </DataArray>
+ <DataArray type="Float64" Name="p" format="binary">
+AQAAAACAAAAAgAAApQcAAA==eNrtWk2olkUUnl2Ui/y5Xv+9409a0D5a+DngqqVQoCU1lamrIoOyRTQUblu1KASZdkGLLq1cvQwiJOEiIsIiYvq1HytxUdtutznvne/5Zjzf39tc/Wb14Hkez3M4Zzjv9VyFEG5O3DgoVlDfL/6MkeMRxX7xR4wcX9p/oJ77xO8xcjyi3iuux8jxpf3r/Gd7/ohyj/gtRo5H9LvErzFyfGl/RCPFLzFyPKJaED/HyPGl/e/U+buMv1ul83eZ/K5j/7r/6/e/zr9+///v/es62v9uRP9h96/rav925C8z/hL40v0v7T/t/ssR+1/af+DnkR3iWowcj2i3i59i5PjS/gP1bBM/xsjxiGKr+CFGji/tX+c/2/MfqGez+D5GjkfUm8R3MXJ8af+BeubFtzFyPOLSz6k+Ro4v7V/nP9vzXy37X2byS4av+38yf1m4/6X8zW02f1O//3X/1+9/nf+U/ev99/a8/9oR/e0ddv9XI/or4P2+9Pwt9IvjLdNfjh/Wn+qdtn+9/872/TfXX5HpvwAe9VivYHguP8eLCf1X6/x1xl8Dj/qB+TI87lcF+5LjEWlPEXL8sP62I39E3Jccn9uvKrN/FfDj+puMv8r4mxyf6a8MiLwEHvXtXAOi3gGfyy/AXwekuN/fn1dk8qO+9cV3GNBAvxzDG6a/HO+G9Lcd+eN7xX2pGH7S/as68rcZfws8t185Xk+6f6Ee1FuGn9R/3P773en+20z/SY/8wD7c3a/neETMz/G2cP9L+0/7/mrgXsnxpf0Nc181Hd9fS/vX+XfTf78l3X+KE07bH/NzfL3/1vt//B6wX/hept1/Ln/X/rM+/3H3rxtx/7oh978DPcc7Jj/Hl96/pf1L97/Ofzx/OaK/zPCl92/d/2X9S/VfdTR/1dH8VVf//6787//738M+8U2MHI+o9oqvY+T40v6Idpf4KkaORxQL4mqMHF/av85/tuePq
LeLL2LkeES3RXweI8eX9keUm8RnMXI8opkTn8bI8aX96/xne/51/9fvf51//f7X/V+//3X+szN/C7q2nj3iSowcP6DfKS7HyPHF/QHVFnEpRo5HFHPCxcjxXfnLjL8Evs5/uv2Xo/Y/k19n/HXGV2fya8Yf36PdIS7EmOMN6AjVBrG47BsQ9Rp4fK+YF+vFfKjXI/pjvWatsMs/L977H6Ke4oSod5+cO938W++aL59rEnqKy/OXTjZD5Of4gfwhr7zn0ZNJ/xDXV86daFLvccj+2xH7n+O5/BxvJ/Sfdv+pr3bzkydS/W/jj39w4lbzF2vT/hQX8O7aP4e87syDy4h6iptX/jre3GK/+sz+9Zm962G/2oCot8CP6+8z/jbjn+OxXtynqNfAD7t/czyXn+P1hP5cfzneT9r/If1tR/6l+1/KX4Gu1P7N5XcZfwc86tvveEDUG+C5/crxiLRX9eknjqf2L8XVxy8dbxL1kk5SPuAl8KinvO4992yT0FNcv3Pz6WaI/nI8198cbzPzJ50Hf3qvHnjUW8jvGT7XXwX9J70CHvXUV/PIoWVEPcXta143iXraPAGR18AP6ENe9caZp5qUPsTdgatHm3r/rff/Ov96/5/B+6/v6P5r4edVjs/lV+BP9SrgUY/5FcNPu/8e6uX4XH6T8TfAox7zG4b3TH85Xk3Y/9L+df+X2b+yo/0vIT/H1/mP1383Yv9dhh/ox4j+bkJ/br/ie0XeTLh/JZNfduxfav+6Cfe/zvg7yI965Ev3f7XNX0O/OF4z/eX40v6Wua/aju6vjnSZ/GJv2p/ihKh3kB/1yOfuqwbur6Q3wKNewb3KMLxj8nO8GdNf072V6S/Hc/3N8T4zf9LJjL8EHvUe8kuGz/XXQv/beoFHvcb+Mrxn8ud4BzrO3wFv4J3gfVXsDL/PB57ihKinO6WjezboHfBcfnyvFJegw/we/Envgcd62zsq3cuBt8AP6MNd1R549UiT0F
Ncb7v5WJOot81Dv88Fvo2jL9UT8vq7jx1uEnqK+92nDjWJ9yKhXzlegQ77K+bCPR30bTwg5iedBn96rxQ3kJfiEvODnuJqDt4J/Znu+uvC71cgH8XVuv6+Upz6Kj/68OEmoae4PXvX2iZRL+kM6aFeijvIS3HK615466FeQk9xf3TN4d4Q+xV/XqW4Ah3u16V/N1xO6SlOOPDzKug5XjH5c7xc6N+rbT2wLwfqpb0MOkID7xX1ZuW761L5PegH6qX41n4doYX8qKe4W8Gx+q8z/ZdQL+rbeGb+GvQcj6ggP8djvdQnQqy3jYOOkPrqab6gb+Mbw77DekDP8YiUl5Dj8b3SnrS0N+E9UdzBXrUr32va64spfRtfwb78pPOUF94rxZfmciGjXwz6xaSe9u3G/r1KddGelOvD9xTqpbhZ379XKU571Vx880gvoae4/fvYM71EvaTzlB/qpbje0J+X4pTXH7x2qpfQU1y/+/yLvcR7pT7pgPheKW4X+vtK8bavAVFPcRsQ3yPpfEDk27jsz0vxNu986BfoKS7n4R1sbP+fB333bbLeEF+qx6b0bV/FAy/3EnqKa3vx9V6i3jbvfPh7UG8b3wx1hDjlNdffP9tL6Cnuzt94u/cPfuWwZA== </DataArray>
+ <DataArray type="Float64" Name="T" format="binary">
+AQAAAACAAAAAgAAA5h0AAA==eNp1nXd4FFXbxgdFRIp0pAUGpIMUkeLOCEMLvQkBVMpEQEBEaiCENiAdAelN5dCRYgFE2BnJUJUqgYQmwlBFKQJSRKJ817U79/n8Jve3/7zv/vZ3Pd48PJycnD3ZKMr/eRj/+f/q05L9nQA3/sOtNRsmO/+PH3ne5d+GDvHB3a/qmqw+PPNSznE6eR3cnLdjjk581PWm9umkEx9cWbClr07y4iG6DyrK8oJbpVvWYr7s5dC0ezbxwY1rS150SF54nshdjeUFtz5aWZ/5qGvdzJGF+uC/1yxE++v3SayNE+x1cGvxM5/ozPdzqhv6DdOJD65seNyL1Zfemb870vnzuakvmKoz369rbZ5W1WG+z5Xus+JZXvTJq7S/NcsLbmyuUJf56Ku3vXgZ5oOLuC35dZIXnli0+jmH5AVXGxepoxNf1m217Reb+JIvmZZbJ/OKh3X6zAGbzCu4mSXTVebL5/u+X24TH9ztm+rYZF7hmWqj8xqZV3Bl3md3NOKjrjpi31aN+OCix9H9GsmLh/ngrUSWF9yMOTCX+XI9abutLfPBrebl+9ksr++ZG8d/zPKCmxtbLWM+6lqZZ/dgPrjx3cgEjcwr+uRW/iddI/MK7u7ve1MjPnJ6Ba+e0IgPbtweu0cj8wrPfbnnVpvMK7jZees1jfioa+T7Z7pNfHD1VoVklhd9Ur5uvYHllfzkniXMR1+NNd+PZ77kXboN1kheeKJDzR42y+tzr0qP+RrxUVcd2DVkE1/y0u/019j676+TbuwAuv6Du4dSl+jM99cp9eTMYTrxwZWnfQfT+r4nSsRMoF9/fC7uVJqtM9+vay5e+o5OfHDv4OzeLC/WSXXQw9YsL7j3Vmwz5mNOrTP1yjAfXPm4c0mWF5734b1XWV5w905lg/moayzLlI354O6PNQrS/vp9Ul8rtYq9Dm59mD5PJz5yijkrRurM97lb9OgH9Ou/7xmvblzKX
ge3Ks+coxMfdb3j/wzWiQ+uPJPYl+VFn5Q6j96keX3uHtjbkPnoq1cta0XmS/52qaIsLzxRP09zlhfcEEn1mY+64vDlksyX/NA/L7F5levk6or/aGRewZXnfnqkEV+uUy0yp2rEB3cfdDymkXmF5+Y5c1Ej8yr58BG3NeKjrtl/7HaN+OBm/t/3srxYJ42JIzawvOBuQv811PfXVZH95/HMl3z5j2NoXt9Tnj/2CcsLbtzY+inzUddaka8388FVa8sQjcwr+uTN/SUTm1fJ6yy7pxEfOc3PJp7WiA9u9pl1QCPzCs8qN/MvjcwruJdreGTfEfRR1yx19phGfHB34uL9GsmLPqmxqV+zvJJ/NW0589FXMazwZOaDWyfzJbK88JTxJdeyvOBqseeWMR91vd3JY5kv+Y1mCRrZr8pHpztl2X4V3Gi5sgHz5X640ainYeKDi6fzCttkvwrPa/5ziO5Xfa6kjGrFfNRVEz7My3xw8bVeWmN5/Ye7cWlamOQFtz6aeZf58vupFes2hYkPbuUt/GOY5IXnjW75W4jkBRd/7H8SIj7qKke77AoRH9yqnvtEiO1X0afKO95i+1VwNb5WM+Yjp/HW2urMB3cXtCqpkf2q9PIWymuT/Sq4mFemkUZ81DUvpF0PEx/cqrW2KMuLPnlVtawsL7iZnvwoRHz01a2x3AsRH1ysanksRPLCUztn3xUmecGV3n/fCxEfddXavRaGiQ+uTM5xMETmVX4/XPKdCWxewc27UwXz5fO9WbowH1y96QwLk3mVdY20hWxewY3EiuuZL+vWjR/IfHB1w5/jQiQvHl6WxjVZXnDlYFJb5sv1tEn2qswHt2ZtrxdmeX3PmDetCcsreb6unanv1xXter7GfHClwpjSITKv6JPX6sQ2Nq/gXtnua5iPnOJc21nMBxcvfjUmROZV1g3lGhgm8wruvu8tDxEfdd0vS8SGiQ+uLpiUyPLKPu2t2Ivm9bm3fXIH5qOvyo
LUSswHN/d46a+TvLLuyWtqmOQFF30Xtg0RX9Y9taVwmPjgap6WL4XIfhXrpDHx4ltsvwru9VHjmI91Ssm2pzrzwd3RTyuz/ar0tm+sx/ar4OaybC2Yj7rqXacA88HdP2aWYnnlulotPivLC25suvYs9TGnjW0vRHxwcbrtzyGSF57ocPJWiOQFV3KOeBwiPuoa+ZP2hpjvc6vDyJQQ2a/KPm3f2p3tVyWvNqkt85HTHLC7JvV9bjX7sizbr8Jznz0Yx/ar4EZSnVbMR12l+fVXmA+uHs5ZmuVFn5RLv+RkecG9T0s+DREffXVbjL0SIr7kQ7qmhUheeCJhUmaWF9zcOSm67wj6fl2vaOVzIeb7XD14+XiIzCvWSVHq0bYQmVdwkb5zM/OxTnmHJ81iPrhx5NupdF7hjSvyKZtXcGvvv2up79f1KuhDmA8uYkZYNC/WyXG7e9G8PncHVolnPtZVo2ZafeaDK1YpheWFp3pPm9C8PjcGeR2Zj7rW0O9rhIkveerNmiEyr+iT2baUw+YV3MvbfhPzkVOcbz2X+eBeasyEEJlXeMaEuC1sXsENcesL6vt13e7TpzFf8pXXx7G86JOSbez7NK/P1UsF3mE++qpcz96A+eDqqSGFWF54bvr9eJYXXP3wcmfmo655q0pZ5oMbF/tXDbHzV5yTVs4zjZ2/gnu7R6/Qme+fU5o7L8frxAdXRsaPoPV9zxhceR49//W5cSp/M535fl2x/Ot+OvHB3dm1SrK8OCdV8s+sQ/OCV+zXjvk4p1IKnM3NfHB17vVyLK/0zisNaV6fKz1E9Pw14KOuOqBAEeaDi7u1jrH3f9Ano+jubrT/PhePfxzjEB85jT+K1tSZ73OxSo1j7+/DE6+5ifT9CZ+LV5pP0YmPuqaV3lYnvuRTNpgsL/rkdnmcQyd5wa2DTys7xEdf3QGXLmvEB1fuXHvWIXnhmYfblmN5wY0rF2sxH3XN5xWF+eDq2vdysXnFOak
1oOSvGplXcDG2j8J8nFOKL5KTNeKDiykTT2pkXuG566/dZfMq+ar6azTio65XeM6PzAe3es0cw/LKc9V7B+ZrJC+4uem3L5mPfaqodb0/88GtmNoTWV54yi/nBMsLrr6e3oH5su6MMsOZD+4dfi6yrwzOq+xT3xa2RuYVXHzd9meb+DJntZ1zNOKDG4/Cm20yr7Ju/cdpGplXcO+P+Vc14qOu8WHn6N9PwAf31lbcyfKiT96NvX01khfcG7R6qk189NX6S2mgER/c+KtUvE3ySu/3ahNpXp+LMXvnU9+vq26vF898cO/88WFs/cc66XV7eTb9+uxzddVra3XiY52y2k/qrRNf8ktnx9L3l33PXXZhOb0f4HNxbWFPnfioaw6dkagTH1ytX7guyyvX1QJlDZZX8nkX4pgv19XyZQoyH1z95uErLC889ezvbVlecOPnG0WYj7ri4oWyzAc32iTvpn//fp/Eu9em0a/PPrdyvzfSIT5yWk3UeIf44EqZmxXo/U/fM3u+3i/w+p16/5dHnmfw/brq6E5DdeKDy+eBvOiTsXNXZZYX3BSjm+nER1+NP75Sqe9zs/a412le3zOLXYjeRwnkBZfPAz7qqqWKRnIHfXD5PDCvcl1tuSTyvm5wXsGtD37JzHy5TsXN2a8RH9wrsdTTyLzKupWHZmLzCq7uPHdfIz7qGmVa/8Z8cKtjzed0ltdfJ93khSs0khfcaJKWzHysq2bS+dnMB/eO19imkbzwjDw/p7G84GJh4t/MR11rXulzGvHBhZcts07mFX0yyxcsxOYVXLwxpaFOfJlzXJcqOvHBXaN5Rzav0tu8LXIvJjiv4Hge9FFXecmwdOKDy/UjkFf2qX31l1lecK9alpbMR1+V8tVeYT64KNM4juWFpzr9huokr+T+86CPumr7f8foxAeX/QjsV3FOagyf1IjtV8G9c8W7Ml+eqzavXIz5kr+dXoPtV+GZ3/7Qhu1Xwd06BzMzH3WN
UillmA8uRib/EiJ5cU5qHen4Z4jl9bmwB2VnPs6prOTJB0PEB3cHHLsUInml16XgvyGSF1wZnuKGiI+6VnKfVOaDixtnoueLgf0q+mR9ahTVyH4V3NiXUsMmPnKKO8Pvh4gPLmrmymKT/So8s7mIZftVcK/2jN7MR12vWJZyzAf3zhRowfKiT8aGtV6I5fW52mEs9/2+mjvfXRkmPriZ7ysvTPLCM5T3yrO84MbSOu2Yj7rKgmOZbOKDWzHbK9pkXuW56oVaK0JkXsFFu9M7mC/PVdVhScwHtz4+P5fNq/Sea7CRzSu4d/V8b+r7db3FxRcwX/LH+maWV56Tfu12ZHnB3erxQ5mPc1Uxo+qbzAdXpi5OZHnhmWnVprO84Kox7Bvmo66xcrhgPrh37N4PbF7RJ7PJXx+xeQX3VvZKCBEfOd0NY86GiA9uddOeZfMKTx2xY1+IzavPzSaf3AoRX9Ytljf6vmrAB7d2Do6ehwTyok/q3P5nWF5w0f+jZzTio69K6qC7IeKDi9vvF9ZIXnhW9vSKGskLbrnr2jMfdc1rExowH1xJf/I+269inXRHvzCD7VfBvYZjtzIf65Ry6sYg5oMb8V8sZvtV6T0p/Qvbr4KrieWi+/eAj7rqd4N3a8QHdxef+J3lletk6MdEllfy6plXMB9zqjz7ZXub+OBiSniUzfL6njWn5kGN5AX3Cr70J/NR1135YJVNfHDj3vSjNtmvok+i9vEqbL8Krq4f2pH5yGlWnp6P+eBqmySd7ldR9/Uulk72q+DyedBH3R7JPXXig+N5MC/6ZPx1rCjLC64816gR89FXL7b+Q5v5Pvd214p+vxXIC0+xJvfXSV5w+Tzgo65SRDRxmA+O738C84p1Urn0Z342r+Bu13frasTHOuW+cmi0Rnxwd9fW1Wxe4RlDXnyPzSu4OSbXJxrzUdcWP2nEBzevZnvE8mKdVKqvHsXygiufnV/FfKyr7sl1c5kPbr272mF5p
Zfc5yjLK/mScg+Zj7rqnvNXmA9uPe7zIptX2dej732nkXkFd1fnuKgRHznVLw6oOvHBrbLDmrJ5hefuPpVVJ/MqOZ4HfNQVbY4P0okPLs/PAnllX/OUjjwP5gUX5yY0YT76quwrV4v54F7V0ibLK72JpQbpJC+4/PMEfNQVL0ycrDMfHPUC91Xl8/WP2tH7qj4XD4t0Z778efgeT0YwH9z4tuOkMLmvCs+99tOb7L4quNH0af0w8VHXS72TxHxwY2jTD1he+bhycwXLC25tbvQ18+V94GlzDzEf3O259VSY5JXe2c9WhkhecFVZOy9MfNQVeecfDhEf3Pqq0M4wua+KPimHhsWFyX1VcOP08cjfa9CXfU2bYoWJD26eOJPE7qtK76/CdcPkviq491nfesxHXeX+/l5h4oMb7Zr3YXnRJ+N8yXUsL7iRf/5K5qOv4o35KcwHN7NOPczywjPmhGfRvD43BoZnMx91rVzFv2M+uLIkbQedV/x7jn3pPptXcKti40w28eW/pzaxMTbxJS9RsYJN5hWemnI8eu4SmFdwY8Ctq2Hio646rVNxjfjg3pA/c9kkr/zzPNOhMcsr+YdxbzJf/jzQsXwDmC952rYkm+SFJ3IVjZ5/BfKCm1Wb1LaZ79c1bhUZqBEf3FjwwLTZvKKvFz5/zOZV8qoH7zMfOa3x60vZxAd3s7YobpN5hackn/LovIIPP36J+bJu8rSsNvHBvfz9s9skL/rkHX21hU3ygquP8sQyH301fzg2hPo+F63yD7RJXnhu3MXqNskL7tZ7UoP5qGtszPk288FFYkxXm9xXxTpprnrYkN1XBRfb4rqHme+vU0Ib1Zf54Gq/TyaHyX1VeIZ3pGeY3FcFd3/I1ChEfNQ19i2YGSY+uHWjcT+WV66Tr4+eEyJ5wa2nVb4JEx9zau75IxwiPrj5TLYzLC889/Pm28IkL7i3PHVuiPioK5a+4YWJD26euuKw+6rok/l5kzzsvi
q4OTjuzTDxkdNr2750mPjg3q7LI8Pkvio89dz3kX1D8L4quBsa0CZMfNQVCZusMPHBzfSziWGSF30yBrVsy/KCi1mdVjEffVXOHx7OfHD10MojLK+s2779OpYX3EopsIL5qOte/y6F+eBmu5YHw2xe/XVSOb3scojMq+Q56j5jE1+uq9qiHBrxJd87L3LuEJxXWbfA+udtMq/g5uMcV0PER11v5RvVbeJL3nzjixrJi3VSVEp7TSN5wdWVXdvbxMe6av77djeN+ODu0dqjWF54okyJt2ySF9w62rmWRnzUNQa0Gm8TH9z4rYmpkXlFn9R+nZazeQU3YhY9YL7Meb/RAeZL/uLwEjaZV+m1SnnM5hVc7C7/J/NRVwn9VMomPri16feiNsmLPln5+t5jeSWfdjPWJj76qlgnijAf3J2+JfJ1OkNe33Nb92lhk7yS7+/YiPmo67qPhjAfXBxa0Z/tV+XPw+qrFrD9Krg5dNhK5sv98PkxLvPBvUfvHLbJfhWeOX34Qo3sV8HV8tum2MSXdVMK7dKIL/ldZTPLKz8PoNIL11lecDP2i/vMl/34bXweh/jg7gtFizssr+8Zo3+7rpG84O527axNfNS1OpXNqxMfXNlTLHLvNrhfRZ/ceRc+pftVn3tX+i+kvp/TrPDCPuaDGwfy7GL7VXhurSfj2X4VXBi1JlLfr+v2XLeB+eBGwvUvaV70ddea2ywvuKh74Trz0Vf3hw4FHeb73PWK53VIXnjeoLKpLC+4t3rCSer7dUXvLv9Q3+fibG/FIfMqHwUmRO5lBecVXP3ndmPmy38Pfy14l/ng5soHAxwyr/CM2oNDOptXn7uxZyo5xJd1hzXqoRMfXMnVvAPNi0eHgtNZXnCxu+lC5v/nkdXgflbjP/8bzAvPVPN9rJO84EbLT0Y7xM9QP+AHXw/Oq+zT340Mh8wruNqpd4j5yOmlbO3NfHBjmt7DIfMKTy3Ss4xD5hXcXfFrOYf4qKs
0qtSG+eDGyPfbsbyyrwmNZ7O84O685z9mfrC/QT/4ejCvrPvmt8NYXnB3pDmC+sH6AT9DvuB+1V8njfkxkzSyXwVXbs9cZRNfrlN7ikU+9yjog4s194+w/So895UN622yXwU3X9oxWSM+6noP+5+wiQ9uzGn6DcuLdVJ9PemURvKCW8+deGATH3Pq7i0SvXcc8CW/fLiEQ/LKuovapNskL7iVt+oZjfio6zUQpR3ig5svvPWMTvarsq9vT2/I9qvg1q24RcxHTuWDpA+YD+6mhXex/So8r1uhz9h+Fdx7d+F86vt11amv7mM+uPJvQjLLiz4Zw5Lnsbzg6uy035iPvqrdtJ3Ml7zj8eh6GsgLz812+TbLC27uuXGN+bLuFiu6bgd8cPH35twOmVesk2bczfI6mVdwNW+lJg7xsU55fctHPvcw6IOrRS4OdMi8whO/5mztkHkFNws4FXXio67odzXBIT642ySpPcuLdVLVE5J0ltfnVutFixziB9fXoJ/h9UBeeOrmX5c5JC+4Mj1+lE78DPUDfob9QWBe0SdvWpZrbF7B3dQPNIf4yOn+nT038yXP8STydS/DvPqe+Lun4bB5BS8cV4f5qGsOvtyb+eCWFo53SF70yd0TX5vlBReVY2YwP9jfDH7g9Qx5fc/r33s2yyv516lTmZ+hfsDPsD8I3FeV56Q1Gr7M7quCe2c+ac98nFOaKd0bh4kvedsR0feXA/dVpbet0HB2X1XyhA9WMR913Yc5poWJD260qCrCJK88Vx1do1+Y5AU3fug3hPr+OZXRuOFs5oMbW/dPC5O88LzN+TaFSV7J9ZrfMh91jTpjfqI++KMhU9l9VdnX6auPhMh9VXDryuIHzEdOY3ko8rluGXyfi0krIp9bmOG+qu+pn4wrwe6rSl7hfizzZd2RE6Kfqxjwwb2uCyOfGxnMK/uk9PTCJC+4kpIaOVcO+uiru6fuDur7XJ1R9WKY5IWn1ioSfT8vkBdc
OXKjM/NR10hMi56HBnxwpdiF6PudwXn1z0m9pTE72LyCm9mabmG+PKfM2fQS88HFpqs/s3mFZ/U+95DNK7h4o1P0fYmAj7rq+M9Vm/jgypi510IkrzxXHbUmm03yghs5e2W2iY99qlGobw2b+OCqteYVm+SFZx14o6nN8vrcTZoUObfO4Pt13eJrBtnEBzeU+NoamVfZp9T3I+/LBOcVXBTbrdjER053Q8UnzAf3OoVKsnmFZ72VKTubV3CR73i0vwFf1l0xOvrzRwEfXIxNKKKRvOiT+UfdciwvuDGmVGub+OircfXj1swHt+5PGcLyyrq9FzVjecHNf+8PtomPum7WV960iS/5F82i7z8F7qtinVTq7xjI7quCq826LWI+1imRo0D084IDvuRnwpHfWxG8rwrPu+nuYvdVwY1B539jPuq6S5Mjv1cj6IOrv9y9YrO8WFfblRlvk7zgytMB65mPOTWc8c2ZD26sqjLEJnnhieMNTtgkL7ixI3c681FXvHx3KfPBlV/f3muT+6rok3XkdF56X9XnrjFMY77MeXRC9Pe6BHxwJYce/b1Fgfuq8NTfM/XUyX1VcHn/N+Cjrvvg1cjvVcrg+1zWC+RFn4yhgyPfdwbzgpvlD7ViPvrq7V53yyY+uLt1aQGH5IUnDm5NcEhecHl/Oej7dd3YUfUc4oPL+7+BecU6aSmZ3mbzCm62mTWe+Vin1Bu939aID27kHhX9eeLAvMJzf8y8gc0ruLJh0wnmy7rxQ77SiA9u5PrzpEbyynV140fzaV7w2G3f28THumqeThxmMx+86NoEjeSFZ2VZdJXlBRd1BuZymI+6FbetYj64MqBN9PO/AvOKPinVNqbbbF597q66HJ33gI+cyuPtkc/tyeD73GvcLXq+FJhXeN6LA1s7ZF7B5TwHfNRVlF2Rz0UL+uDy8wQCedEnkfvXWiwvuNnrq+jvgwv6fl9FsXvR33cR9H1uPTOupE7ywrNS2kxxS
N7/5f6//4CPuqLg5Ojvmwn44PJ5YL+Kc1JjQtMubL8K7pVaEsd8ea4as3QC88HV4klj2X4Vnvned4ttsl8FNzY8H/258IAv68Y7u23ig7uTjk7RSF6ck4o+Wb5kecGVg7fXMB/nVOaDC2nU97lZ5dtjLC889XTLGzbJC269k3ezRnzU9XZezucQH9ytEHdWI/tV9El5odEwtl8F9x6P+JT5yOmu0wTzJY/bs5ftV+GZefYksf0quOX8tIr5qCsqX1zJfHBv4cWjLK/sq7buB5YXXB3S4Db1/b561fS7zAf3ZmWPnGsG88q6JQsfYXnBrUqLHjIfddWUqg+o73PjmyGqQ+YV56Sec+Epm1dw0abzX8zHOaVytE45h/jgxhynpEPmFZ6o/UR3yLyCm02rPqsTH3XdFjG9HOKDuw3uV9JJXpyTWkVvtHVIXnCv0Olm1A+crwb9DOevgbzwvJf2znRIXnDxXdcOOvEznP8G/Azn04F5RZ/UfDGFHTKv4GJ9LsMhPnKKWy0aMB/cKt07cq4dnFd4ylivuEPmFdz9/XETh/ioKx6PiaW+z90rwwexvOiTmb7ufZYX3DqfYzbzM7y/EvAzzEcgLzx1T4+BLC+4WnDoYuZneP8n4Gf4+w/sV+U6WW7GcrZfBRdbjkQ+FyDoY51Slqcctokv+Q+717L9Kjz1WoW9GtmvgquFE+8wH3WtEq9tYj64uaH7HxrJi3XSvTXjgU3ygrvX259nPubU7Ncy8r5m0AdXhqZmdkheeN6iFjdsltfn5tmYP23io661/dXo+1YBH1xZOvCxTfar6JPRr0H0c7sC+1XJ2++OfO5x0EdOc/27hZgPri5r0EAn+1V41q/X++pkvwou968BH3VN6/D7OvHB5fNAXvTJS78XeV8ymBdcOGdqOsRHX72Cm3IyH1zcW9SC5ZV15yzp6rC8Ppd//oCPuuaRLsN14oPL73+C8+qvk+6oE7EOmVfJS1yrwnysUyJt+U
DmgyvzZ3V0yLzC8za/+rpD5lXymrNfZT7qqpb1rkN8cLe83YXlxTppJOqLWF5w5dchFvMzfH0P+Bnu3wXywhPxOaezvOCizYiJzM9wvyvoB//7gXmVfd3X3GDzCm42K/AS82XOGPGBQ3xwtU5sHzqvvmccyp+gs3n1uZzngI+6Rtydj3Tig8vzg0Be9MnKP2CBQ/KCm/9cWsL8DF/fA36Gr8+BvPCM765Ffm4/mBdc/vkDfob5C/j/ef1Ovf8BiUc4KQ== </DataArray>
+ </PointData>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
Copied: branches/freesurface/tests/constant_temperature_boundary_conditions.prm (from rev 2445, trunk/aspect/tests/constant_temperature_boundary_conditions.prm)
===================================================================
--- branches/freesurface/tests/constant_temperature_boundary_conditions.prm (rev 0)
+++ branches/freesurface/tests/constant_temperature_boundary_conditions.prm 2014-04-09 14:50:17 UTC (rev 2446)
@@ -0,0 +1,73 @@
+set Dimension = 2
+set CFL number = 1.0
+set End time = 0
+set Start time = 0
+set Adiabatic surface temperature = 0
+set Surface pressure = 0
+set Use years in output instead of seconds = false # default: true
+set Nonlinear solver scheme = IMPES
+
+
+subsection Boundary temperature model
+ set Model name = constant
+ subsection Constant
+ set Boundary indicator to temperature mappings = 0 : 0.0, 1 : 1.0, 2:2.0, 3:3.0
+ end
+end
+
+
+
+subsection Gravity model
+ set Model name = vertical
+end
+
+
+subsection Geometry model
+ set Model name = box
+
+ subsection Box
+ set X extent = 1.2 # default: 1
+ set Y extent = 1
+ set Z extent = 1
+ end
+end
+
+
+subsection Initial conditions
+ set Model name = function
+ subsection Function
+ set Function expression = 0
+ end
+end
+
+subsection Material model
+ set Model name = simple
+
+ subsection Simple model
+ set Reference density = 1 # default: 3300
+ set Reference specific heat = 1250
+ set Reference temperature = 1 # default: 293
+ set Thermal conductivity = 1e-6 # default: 4.7
+ set Thermal expansion coefficient = 2e-5
+ set Viscosity = 1 # default: 5e24
+ end
+end
+
+
+subsection Mesh refinement
+ set Initial adaptive refinement = 0
+ set Initial global refinement = 5
+end
+
+
+subsection Model settings
+ set Fixed temperature boundary indicators = 0, 1, 2, 3
+ set Prescribed velocity boundary indicators =
+ set Tangential velocity boundary indicators = 0, 1, 2, 3
+ set Zero velocity boundary indicators =
+end
+
+subsection Postprocess
+ set List of postprocessors = visualization, temperature statistics, heat flux statistics
+end
+
More information about the CIG-COMMITS
mailing list