{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", " \"Open\n", "" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "

\n", "\n", "# Tutorial: basics of FARGOpy" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Not running in Colab, skipping installation\n" ] } ], "source": [ "try:\n", " from google.colab import drive\n", " %pip install -Uq git+https://github.com/seap-udea/fargopy\n", "except ImportError:\n", " print(\"Not running in Colab, skipping installation\")\n", " %load_ext autoreload\n", " %autoreload 2\n", "!mkdir -p ./gallery/\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If you are in `Google Colab`, install the latest version of the package:" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For this tutorial you will need the following external modules and tools:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Running FARGOpy version 1.1.0. A major refactor has been done in versions 1.1.X. Please check the documentation for more information.\n" ] } ], "source": [ "import fargopy as fp\n", "import numpy as np \n", "import matplotlib.pyplot as plt\n", "from celluloid import Camera\n", "from IPython.display import HTML\n", "from tqdm import tqdm" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Let's `FARGOpy`" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "First we need the data. Let's download a precomputed simulation. Check the list:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "fargo:\n", "\tDescription: Protoplanetary disk with a Jovian planet [2D]\n", "\tSize: 55 MB\n", "p3diso:\n", "\tDescription: Protoplanetary disk with a Super earth planet [3D]\n", "\tSize: 220 MB\n", "p3disoj:\n", "\tDescription: Protoplanetary disk with a Jovian planet [3D]\n", "\tSize: 84 MB\n", "fargo_multifluid:\n", "\tDescription: Protoplanetary disk with several fluids (dust) and a Jovian planet in 2D\n", "\tSize: 100 MB\n", "binary:\n", "\tDescription: Disk around a binary with the properties of Kepler-38 in 2D\n", "\tSize: 140 MB\n", "pds70iso:\n", "\tDescription: PDS-70c - isothermal protoplanetary disk and circumplanetary disk in [3D]\n", "\tSize: 2940 MB\n" ] } ], "source": [ "fp.Simulation.list_precomputed()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Download, for instance, the 2D simulation of a disk with a Jovian planet:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Precomputed output directory '/tmp/p3disoj' already exist\n" ] } ], "source": [ "simpath = fp.Simulation.download_precomputed('p3disoj')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Once download it, we need to connect a `Simulation` with the directory with the simulation results:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Your simulation is now connected with '/Users/jzuluaga/fargo3d/'\n", "Now you are connected with output directory '/tmp/p3disoj'\n", "Found a variables.par file in '/tmp/p3disoj', loading properties\n", "Loading variables\n", "85 variables loaded\n", "Simulation in 3 dimensions\n", "Loading domain in spherical coordinates:\n", "\tVariable phi: 128 [[0, np.float64(-3.117048960983623)], [-1, np.float64(3.117048960983623)]]\n", "\tVariable r: 64 [[0, np.float64(0.5078125)], [-1, np.float64(1.4921875)]]\n", "\tVariable theta: 32 [[0, np.float64(1.4231400767948967)], [-1, np.float64(1.5684525767948965)]]\n", "Number of snapshots in output directory: 11\n", "Planets found in summary.dat:\n", " Name: Jupiter, Initial pos: [1.0, 0.001, 0.0], Mass: 0.001\n" ] } ], "source": [ "sim = fp.Simulation(output_dir=simpath)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's read the field we want to animate, eg. the gas density. It is important to notice that we will animate the field in a slice corresponding to $z=0$:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "fields_spherical = sim.load_field(\n", " fields='gasdens',\n", " snapshot=[1, sim.nsnaps-1],\n", " coords='spherical',\n", " slice='theta=1.56'\n", ")\n", "\n", "fields_cartesian = sim.load_field(\n", " fields='gasdens',\n", " snapshot=[1, sim.nsnaps-1],\n", " coords='cartesian',\n", " slice='theta=1.56'\n", ")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And the plot can be done using:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "((64, 128), (64, 128), (64, 128), (64, 128), (64, 128))" ] }, "execution_count": 86, "metadata": {}, "output_type": "execute_result" } ], "source": [ "snapshot = 5\n", "\n", "phi = fields_spherical.var1_mesh[snapshot]\n", "r = fields_spherical.var2_mesh[snapshot]\n", "\n", "x = fields_cartesian.var1_mesh[snapshot]\n", "y = fields_cartesian.var2_mesh[snapshot]\n", "\n", "gasdens_plane = fields_cartesian.gasdens_mesh[snapshot]\n", "\n", "x.shape, y.shape, gasdens_plane.shape, r.shape, phi.shape" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAABBwAAAIRCAYAAADp1a7hAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQABAABJREFUeJzsvQucJFV59/9U9/Rcdy573+UOYrgodwRBY1BRY4gJ0fgSTQRR9K+vGAXjBV/EKAkkJgoaUVRUYhRF83qL+BIRRERQFMULCgpiFvY+uzv3nZme7v5/nlN9qk+dPqfqVHV1T3fP7/v5nN3uupw6XV3d0+dXz/N7vEqlUiEAAAAAAAAAAACADMll2RkAAAAAAAAAAAAABAcAAAAAAAAAAAA0BUQ4AAAAAAAAAAAAIHMgOAAAAAAAAAAAACBzIDgAAAAAAAAAAAAgcyA4AAAAAAAAAAAAIHMgOAAAAAAAAAAAACBzIDgAAAAAAAAAAAAgc3poBVMul2nbtm00PDxMnuct93AAAACscCqVCk1PT9MBBxxAuRzuCWQB/tYDAEB7/m2an5+nxcVFaga9vb3U39/flL5BMla04MBiw8EHH7zcwwAAAABCPP7443TQQQfhrGQA/tYDAED7/W1isWHdwADNUnPYtGkTPfbYYxAd2oAVLThwZANz0hEXUj7XS52Kly9QN1MpFZd7CF1JpbywLMf1cn3LclzQed89K/GzXyov0k9/9+ng71M3cfXVV9OXv/xleuihh2hgYIDOPPNM+ud//mc66qijIvf70pe+RO9617vo97//PT35yU8W+/zJn/yJ83HlueQfyiMjIw2/DgAAWGlMTU2Jm7RZ/m3iyAYWGy7O5yjrX4b8C/fDO3aIYyDKYflZ0YKDTKNgsaEn38mCQ+eO3YUKId2lOee1QstBt1+vK4lmv5cr+bPfjWl+3/3ud+kNb3gDPe1pT6OlpSV65zvfSc9//vPpV7/6FQ0NDRn3ueeee+hlL3uZECv+9E//lG666SY699xz6Sc/+Qk99alPTXQuWWyA4AAAABn9bZqZS7bzqkHj4iGvQn0Z/83rqSzPb1xgZkULDgAAAABoDbfeemvo+Y033kgbNmyg+++/n571rGcZ9/ngBz9If/zHf0xvfetbxfMrr7ySbrvtNvrwhz9M119/fUvGDQAAwMD4RCaCQ84jymesscMBqb2A4AAAAACAljM5OSn+X7NmjXWbe++9ly699NLQshe84AX01a9+1brPwsKCaGooMAAAgIw5cEP9suIS0a694XWLRaLd+3D6VzAQgAAAAADQUrhyxJvf/GZ6xjOeEZkasWPHDtq4cWNoGT/n5TY4/WJ0dDRoMIcGAIAmUOipb/nq1LInX1vGjyPo8ZrTQPsAwQEAAAAALYW9HH75y1/SF77whcz7vuyyy0T0hGxsFgkAAACA5QEpFR1CJxntteNYK6Xm1PjtZLw8qkWA9v4cNzoGfO7bk4svvpi+8Y1v0F133RVbXo3Lmu3cuTO0jJ/zcht9fX2iAQAAaH/yTfBwiI6pAK0GEQ4AAAAAaDqVSkWIDV/5ylfojjvuoMMPPzx2nzPOOINuv/320DI2jeTlAAAAlhFTJYhcjqh/+W9YgPYCEQ4AAAAAaEkaBZe1/NrXviZquUsfBvZZGBgYEI/PP/98OvDAA4UPA/OmN72J/uiP/oje//730znnnCNSMH784x/Txz/+cbxjAACwnDyxk2hwgGjVAFFfVWRgz4ZN68Lb8bKNa63dIMKh+0GEAwAAAACazkc/+lHhqXDWWWfR5s2bg3bzzTcH22zZsoW2b98ePD/zzDOFSMECwwknnED/+Z//KSpURBlNAgAAaAG9BaLpWaLt40TbdhFNzhCVSuaohwGkua1kEOEAAAAAgJakVMRx55131i176UtfKhoAAIA2gqMWlpZ8oWF6jmhximjflC8urBokGuwn8uLNGfJepQkeDvF/b0DrgOAAAAAAAAAAACDhTLKHaGjAFxw2rCHaP080u59o/wJRzqumXAxG+jrkm1DGEqaR7QUEBwAAAAAAAAAA6WFRgaMa1owSzbHwMOe3mTmiQp7owI04uysUCA6g7cvm5QrD1nWV8kLisUaVykMZPbBSaWUZzGYeS/0MpzkOvgMAAACABuA0ChYe2LuBMxs42qFo8Haoks81oSwmMiraCggOAAAAAAAAAAAao1gkmtnvRzYslYlYSGDxgdMqwIoFggMAAAAAAAAAgOTIaIIde4gWi/7jvgLRyLBfMpMjHSJAWczuB4IDAAAAAAAAAAB3OFViZtb3a2DKJaLRVX40QwFTTFADVwMAAAAAAAAAAHd27vErUXCViqFBvxxmChDh0P1AcAAAAAAAAAAA4M66sZo5JAARQHAAAAAAAAAAAOBORkaQiHDofiA4VEsrVgLHk+bg5dOFGbmUamtlObssy89F4eX6Em/nOoZKLqIsZn7B+Po6uVTeclwfnXy+Vir6e5bldePal+vnvt3I4nqvlPzvnkoZnx0AAAAdwPi+ZNuvW21cjLKY3Q8EBwAAAAAAAAAA7szXbtIJ+N5tqUzUkzMvtwgOoPuB4AAAAAAAAAAAwJ2DNoWfLywSbR8nOnAjkefVL7eQq6ZVZEmuuYHrICFw+QAAAAAAAAAAAEDmIMIBAAAAAAAAAEB6lkr+/yKtIl9bzs+9FptGZtwfaAxEOAAAAAAAAAAASEdxiWjflP94cpqoUs1pKJeJpmaICgWc2RUMIhwAAAAAAAAAACSjVCKanCGanuVaU0TrVxONTxDNzftRDixElCtEa8esXeRyfssS3FFvLyA4AAAAAAAAAABwhyMapmb9aIaBPqI1o0SFHl9oYBGCxYa+XqJVg0RDAzizKxgIDtXa71wnXtZBbwaN9u3l+1LVgOfXlSVZ92ejUl6IfL8aIdc7ZD9uSXl9SvRXuThtPedR5z/NNRD1XqdBHV+r3r80x3E9j1nQzM+6K1m/z1mzHNdNFC27drXvl6jvoqREf1f7x/WiEl0BAACAdoFFhcE+otFhX1iQ8OMNa5y7gYdD9wPBAQAAAAAAAACAOwesJ+pt3JsBgkP3A8EBAAAAAAAAAIA7utjAVSrY04HJ58OVKsCKBoIDAAAAAAAAAIDkTM/5lSjYs0GF/RxGV/keDhHkc37LEkgd7QUEBwAAAAAAAAAAydgz4QsOHO3A4oKMauBoh/0LfsWKhSLR2lGc2RUMBAcAAAAAAAAAAO5w6UsWG9aMEI2sql+/mvzIh71TfhWLwX5jNznPb1mSdX+gMVCmFAAAAAAAAACAOzOzvohgEhskvI634W3BigURDinK1C1HSb20JRVtZdjSlpjLuj+XvrPAdXxe3uy2m8+Hy/tUSkXjuNXymWlfT5rry/3azfb9a/R9z6L0YprzbDtfrfxsZ32sdi+zmYY010SW16RrGd4sS2l6VEm9LwAAANAyOFVijUOqxNAA0d5J6+pczqN8xiEJ/h11/D1tFxDhAAAAAAAAAADAnXI57PbIzxcMN4F4G14HViyIcAAAAAAAAAAA4I7nEVWUKAKuUrF9nOjQzf46SbkSfq4BD4fuB4IDAAAAAAAAAAB3uDLFYpFowGwGGTC/QFQwpywzKIvZ/UBwAAAAAAAAAADgDgsNXIWCIxsYU9pEqUQ0M+eXzAQrFggOAAAAAAAAAADcGR70xQSOYAhmlpo94Pyibxo5PGTtBikV3U9bmEbedddd9KIXvYgOOOAA8jyPvvrVrzrv+/3vf596enroxBNPbOoYAQAAAAAAAACI8hJEB24gOmhTuKl+DSw2rB3zt21j/v7v/17MQdV29NFHB+vn5+fpDW94A61du5ZWrVpFL3nJS2jnzp2hPrZs2ULnnHMODQ4O0oYNG+itb30rLS1Voz+q3HnnnXTyySdTX18fHXnkkXTjjTfWjeW6666jww47jPr7++n000+n++67L7TeZSztRltEOMzOztIJJ5xAr3rVq+jFL36x834TExN0/vnn03Of+9yGTnSud4hy+T5reb26MmkdWj5THXdUKcE0ZeWy7i/tsZLur49NLYsZNe7wuqHQtSQpL4ZrDttKZmZxnaQtm2oaT/0+vU0rs+m6Xauu13b+bKcdUxblMtOUL3XezqHspN+fPf/TBbWUbfRxsv2+Slo+E2UxAQAArCTaxcPhKU95Cn37298OnvMNbckll1xCt9xyC33pS1+i0dFRuvjii8WclW98M6VSSYgNmzZtonvuuYe2b98u5qiFQoGuuuoqsc1jjz0mtnnd615Hn/vc5+j222+niy66iDZv3kwveMELxDY333wzXXrppXT99dcLseHaa68V6x5++GEhYriMpR1pC8HhhS98oWhJ4Tfs5S9/OeXz+URREQAAAAAAAAAAGkT6NHD6xFLJX9aTJ+rvJVo1lL2a0CRYYGDBQGdycpI++clP0k033UTPec5zxLJPf/rTdMwxx9APfvADevrTn07f+ta36Fe/+pUQLDZu3Cgi76+88kp6+9vfLqInent7hYhw+OGH0/vf/37RB+9/99130zXXXBMIDh/4wAfoNa95DV144YXiOe/D4sKnPvUpesc73uE0lnakM64AA3xyf/e739G73/1u530WFhZoamoq1AAAAAAAAAAAJGT/PNHWXUT7pv2KFdKQgR/zsq07ifZHR4XmPI9yuYxbRBlOG7/97W9Fev8RRxxBf/3Xfy1SJJj777+fisUinX322cG2nG5xyCGH0L333iue8//HHXecEBskLCLwXPPBBx8MtlH7kNvIPhYXF8Wx1G1yuZx4LrdxGUs70hYRDmkuCFZ5vve974XCXeK4+uqr6T3veU9TxwYAAAAAAAAAXU2xSLRrn18ec80IUZ+WkriwSLR3imjXXqID1hMVWj/t1G8us3cCNx1OX2A/haOOOkqkQ/B88Q//8A/pl7/8Je3YsUNEKIyNjYX2YXGB1zH8vyo2yPVyXdQ2PMb9+/fTvn37RGqGaZuHHnoo6CNuLO1Ix0U48BvBaRR8IfzBH/xBon0vu+wyEYoi2+OPP960cQIAAAAAAABAVzI544sIm9bWiw0ML+N1vM1kzcdMRwZFZN2Ygw8+WPgcyMY3n01wav9LX/pSOv7440XUwTe/+U3hFfjFL34xs9O1kum4CIfp6Wn68Y9/TD/96U+FSQZTLpepUqmIaAfOoZE5LTo2VQsAAAAAAAAAgCOcKrF6JFyVQofXjQwR7ZtqrWlkxf+fby6PjIwEy13ngRxBwDe2H3nkEXre854n0h1YgFAjC7hggfR84P/1ahKyoMEmZRu9yAE/5/ENDAwIT0Jupm3UPuLG0o50XIQDvym/+MUv6IEHHggam0dyCAw/5pAYAAAAAAAAAABNolz2zSHj4G1422WaN6rNVXCYmZmhRx99VFSQOOWUU0S1Ca4qIeGqEezxcMYZZ4jn/D/PT3ft2hVsc9ttt4ljHnvsscE2ah9yG9kHp0rwsdRt+KY6P5fbuIylHWmLCAd+U1lBknDZEBYP1qxZI0wwOBVi69at9JnPfEaYZzz1qU8N7c9lQrhWqb7cFS/XK8qVqSXL1HJlepm0xssyLk8ZxKxL+aUpaeheajJbQu+n+j5HvJdqGT29JJ86VtvjXO+q0D7lRXPJTLVcpj6mrN+zdNdN8pKU9aVkG3tvm3ltuBK+jrMoNdm80ppZjC9L3Etf2svURm2XJa7lM6NIUzpUXl9epXpbBgAAAGhncrlaVYooeBve1oJXNXrMEi/hn9K/+7u/oxe96EV06KGH0rZt20RRAo42eNnLXiZSMV796leLcpU8N2UR4Y1vfKOY4MuqEM9//vOFsPCKV7yC3ve+9wk/hcsvv5ze8IY3BCIH3yD/8Ic/TG9729voVa96Fd1xxx0iZYOrUEj4GBdccAGdeuqpdNppp4mymLOzs0HVCpextCNtIThwisSzn/3s4DmfRIZPOBt4sHmHdAoFAAAAAAAAALCMDPQRTc0SDQ3Y0ypYROdteNs25oknnhDiwp49e2j9+vX0zGc+U5SZ5McMl67km94veclLRNVD9nn4yEc+EuzP4sQ3vvENev3rXy8m/0NDQ2Ie+973vjfYhktisrhwySWX0Ac/+EE66KCD6IYbbghKYjLnnXce7d69m6644gohWnB5zVtvvTVkJBk3lnbEq7D5wQqFXUFZKTr9qe+gnny/9Y64jmuEg2275YpwcMV1nzR3GdvhTnXU3Vbb+FwjHKIoL84se4TDcl0D7fC+Nxql1GhkU7dEOKR5z9shwsH9e7vxCIc0yPEtlebph7/8J2FsrOadgsb/1uOcAgBAht+jXKVi2zhRbw/RmlFLlYpJosUlY5UK2eedR+dpVT7bCIeZUoXOeqiE7/02oS0iHAAAAAAAAAAAdAiFAtGGNUS79xJtH/edH6WnA6dRlMr+Mt5mGUpigvYB7z4AAAAAAAAAgGRwqsSBG4lm5ojmF2qeDr0Fov5eolVDsSUo2N4hwuIhFbkVG7/fnkBwAAAAAAAAAACQHBYURlf5DQADEBwAAAAAAAAAALScXBOqVCDCob2A4CDKFw5RLt+vlUGsL1cWZ4CmG03a+nAtsxllKuda7jL7UpjmcompDCQdzRuzNuuLwvaaXI3kvN5a6Uud/FDNYTY3UDtOef++0HaqoWRpfs+yG0g2E1czwWYeN8ogNutrMuvSmu2GzSjS3Wiy4NZ3VH9qH5bv9Kj3MOsyyK7I154rOdQ072Duuusu+pd/+Re6//77RQWqr3zlK3Tuuedat7/zzjtDVawkvO+mTZuaPFoAAAAqhT2TRHMJ/i5uWhe7CVIqup+MM2YAAAAAAMxwPfETTjiBrrvuukSn6OGHHxYig2wbNmzAKQYAgOWmuEQ0v+h7N3DhQ/6fn/Ny+CiAKohwAAAAAEBLeOELXyhaUlhgGBsba8qYAAAAuFFcO0oDsizm3Lxf9nLzunBJTC6HuXufs6cDIhy6H0Q4AAAAAKCtOfHEE2nz5s30vOc9j77//e8v93AAAADsmyIaGw6LDQw/5+W8HgBEOAAAAACgXWGR4frrr6dTTz2VFhYW6IYbbqCzzjqLfvjDH9LJJ59s3Ie34yaZmsKPXgAAyJylJXs9S17O6x2AaWT3g5QKAAAAALQlRx11lGiSM888kx599FG65ppr6D/+4z+M+1x99dX0nve8p4WjBCudM47+W2oH7n3oQ8s9BLCS6Okhmp4lGugj8pQqE+zlMDPrrwcAggMAAAAAOonTTjuN7r77buv6yy67jC699NJQhMPBBx/cotGBbqRdBIVGxwlBAmQKp02wV8O23USD/UT5HFGp7Hs7sGnk+tVO3XBFTFugRFpy5Wz7A40B6QkAAAAAHcMDDzwgUi1s9PX1iQZAtwkKjQJBAmTK0IAvMkxME03N+FUpvKqHw9q1RP34HgY+EBxEvfheUXtdrb+u1l7Xa8NXlNruWRM+bu2DWinV8lF1XLdbLtTz6uX6nLZTyfUOWfdxfS/U85oF1uMqh/Eixp0rDBkfM+XibG2d0kdpfm94u8WZzN539RqK3q430+1cybo/W9+u14l+HVfKy/+5y/oadznn7tdDoeE+rGh98/d5gPo4ar+I7xGX8WV57r1Kd9+WmZmZoUceeSR4/thjjwkBYc2aNXTIIYeI6IStW7fSZz7zGbH+2muvpcMPP5ye8pSn0Pz8vPBwuOOOO+hb3/rWMr4KsNIEBv4eyPo7thPOCyIiQCwsKmyq/iYqlYjy+cQnrSkeDiiL0FZAcAAAAABAS/jxj39Mz372s4PnMvXhggsuoBtvvJG2b99OW7ZsCdYvLi7SW97yFiFCDA4O0vHHH0/f/va3Q30A0AyBAUCAAAlRxQb2ceDUCo6CiPu85bIXCLhP0D5AcAAAAABAS+AKExX+IWqBRQeVt73tbaIBsNwCQ9w2iIAAK4HCnkmi/gGiXkvk4lLJN4ycnvP9HBwEB9D9QHAAAAAAAAAdCyIY2uO8IwWj+/HYDJJFBV1w2L/gV6zYP+9XrGChYXjIrc+cJ1qm40SEQ1sBwQEAAAAAAHQUZz7l71qaIrESIhiyFCAgPnQnlUIP0Z4JotKwX5lidr8fzcBCRF+BaM0Y0VA/TBRACAgOAAAAAACga0QGeDAsPxAfupPi2DANFMtEnFrBjenvJVq3zq9OkQL2b8i8LCYiHNoKCA4AAAAAAKBjhQYIDJ0hPiDqoQvgMpirx/yIhpk5v80vEu2dIhoe9FMpOKUCAAUIDilK5akl3tTyiPofvFCJS6WMnl5Cr9E/lGpJRNfyhlFE9WEba9RriCqF6UJU6cuocnuNhkhGrQtfH/ElMgV5dbtaXpvXOxjerHd97fFQ7XFudndou/L+fcHjpdkd1nGnuT6yLn+Z5hpI87lIcz1EfYazJup7wIZzqc4Wuaq7ftYjt2s3B/gGv0ey+GmFcG0A0qdMZAU+h80FUQ9dBKdWrB4hGhv2q1Gwf8P4BNHeSaJVg774UHD724qymN0PBAcAAAAAALAiRYZOQb15YCLuhkLWN6caBeJDlyANIrkVi76fA0c9TM36qRab1i33CEEbAMEBAAAAAAAsGxAami9IRG273AIEUi66BI5oWDPqRz4IM8lZp91QpaL7geAAAAAAAABaCkSGZESJAknEhjT7t0qQQNRDF0U9cFoFNwAgOAAAAAAAgJUqNMC3oT0jIhD1sHJAlYruBxEOAAAAAABgRQgN3SgwJPFv6DQBAsJD94OUiu4HggMAAAAAAOhqoWElowoCWYsPap8QHgAAJiA4NFEBt5XWjCoRaCubF11GbpXxOK5jS0KryhtmUTKzmeUzbec59D6HymDyc2Wdsl1OH+vAaO1xYSB43DN2SHhD5XluYnXweGnyidBm5eI0NUK7lb50L3FZaEkpzbYs89hEIr+/lPPQ6GcuivDnLIPjOH6PAADcgdDQnuiiQJYCBIQHkIac54nSmFmScXegQSA4AAAAAACAFSE0sDCa5KaSLihH7btc4nOSMenbNiMdoxUlOJFqAUDnAMEBAAAAAAB0tdDgSpxo0IiokHbfOIEkSb9JBJQs0jGaHfUA4aHz8XJ+y7pP0D5AcAAAAAAAACtGaGhWJMJy9dtIGnBU33HREEkECAgPXcjMHFF/L1FPzHRyfIJosN9vYEUCwQEAAAAAALREaEia0tBuNE1UiPDHUT29GhUNsvAiq63vSy08qPtnCSIeWggLCetXxwsOC4tElYpVcMjlPdGyJFfOtDvQIBAcAAAAAABA04UG9XGc70AzRQmX46c1SI4TB9KQxpTZhay8LNIID+r2EB66nL6CLzqAFQsEBwAAAAAAkInYoE/mk0YENBpBIPevTwdI782QZMKfKwxTK4n1YLBUP3OJYEh6ntP6PdiiHnh5o2IERzzc+9CHGuoDRDAxTTQ1G32KSiWipZJ1tZfzRMsSeDi0FxAcGgxPa2bpPaKhFPuAJGRRUs+lRGZUycy67Yq1L+7c0PrgcX5obfgAhVpoWu+BJ9S2G94c2mxpzyO1x7M7aDnusnRb2chmlnysP1b6H4RJ7r5lUdbUdl6a+v5rJS0r6nFz7X3ddXJYOVhZuEY1yM96ms98/feJOQrBPNGtP14zBYZ2+puWyINBe11R5dfTlFnPMupBf96I8IA0iyYTpxX05P1m2z3viZbpkJBS0VZAcAAAAAAAAA0LDWlw3dckSCQ+VgJhNe1raqUgzVQ00bUZwkiU8KOuNy9LJzzUjwHCQ1syNkw0NLDcowBtDgQHAAAAAADQdkJDbJRCri/SMyFp5Ja7+NFaUSHtWKLEiLgIOjU1pEzTzv3IvvS+Gy2vqe+LiIcuogkpFYSymG0FBAcAAAAAAJDIpyELbIJBlO+Avj5qXUPmkG0kKrRCjIg6F7ovhXpOje+PlvZi67tR4QH+DsvMujHfEBKAGCA4AAAAAACscFoR1RAnGsR6JzjuEzl57h2KmHSvnMmT+lpNkRBR61VBwcUXIs7AUxceWEhwESMQ7bDMrBrMpBuUxex+IDgAAAAAAKxgWh3VoN89j4pCMK23GxX2RgoMdWOJEBiSmDF2OnFCi018MBp4GqJTkrxftXXu3g9ZCQ+oZgFAc4DgAAAAAACwAllOoSFJFYm4iIUkkQppX08SkSHOOyFL0kZluJ6HerGnXnywCgqG6Ie4KiP1x2ud8IBqFgmZmcsk0gFVKrofWGoAAAAAAKww4sQGngRmITbwpFP3VQiLCeEJs7qe15km1HK5XKc+t/WX5vXw5Fe2doUn/aYWRZLzEHX+TOc96lyr14J1G8vYWERwFRIarYghhQcQw/iE3wCIAREOwRdmgi/fiPC8jiED9T3NH+B2/qMdRfTdEnOoYdRrDRkuadupfZQXZ2rLF6tKcpX8yKbaGEbW15avPUTbbkPwOLf918Hjpakn7GNQczAdXb7bqS75cuIahrtcucKNvk/NuHPYLCrlxr9vmvmdZfyxXSk17XgAJBEbGkUXGaK+S/QohXrRwPy94+oR0Gpa9Vsn+u9Nut95SSJPbNEPpsgHdXt5bagRD6H1lvc1ScRDFtEOSLGIYdNaygLPo8yrVHheJdP+QGNAcAAAAAAAWAGod27td5J7nctOpkmdsIsH7iKs3meud5Wzz0McrtummdAnPZdxwn/aChxR4o6rx4JpfU0wCAsP+va68BCsrz4vF6eNVUbU560QHpBiEUN/es+MpqdUZNwfaAwIDgAAAAAAXU5cmHidMCBD32NEB9tk2FVocKW+v+jnrkRHI6aLEkgrLCTtLyoC0VWMiPK/SBo5oos8icwmFTFBluFk4cH2PuvCg79soWnCA6IdAEgPBAcAAAAAgBUkNkRN1p3T6JogNOh3tOvGlkBUcE1rTCsyyLvwrUi1sL3uNEKQPo76KIXsox9U8cEkPOjPdeFBfZ1RaRcu5TR5PUSHDCmXffPIuXmi4pL/XLx5OaJCD9Fgv28qyc8tIMKh+4HgAAAAAACwAsUGW1SDjhrlkFZoMIkJtn3jxtlKkcEmLLTSkyrLY9VVjmih+GAXHszeD2rqhV5i0yY8uKRaQHTIiKUS0c5x//++Xl9ckMICCw+LS0R7p4imZok2rSPqyWd1ZNBhQHAAAAAAAFjBYoNLVEPUNrIv3/wxOkIhbqKapdAQNmguZuDbkO12rTMp7s003cMmPqQRHkz9moQH1WhSP5bfX2+iaAeIDhmwd9J3fDxwA1GPZUrJUQ+79/rbblhj3MTL50TLEg/aRlsBwQEAAAAAoMvIUmxQ97NNZFlssB3PZbnrepX6Kk/2qk9x+6bZbrkrb7lGjLhGI0RVqIr2hCikEh5s/ereD1HCg+xPig5Joh0aER2YFe/rML9AtG61XWxgOK1ibIRofF/i8wy6BwgODH+xcWnMnO2PY0T+YatK2znmFLrmHjoXi4noL60DtAtpjZoapXFTq0LEHYBi4h9M4ce1MpZMeb725Z2fnwoe96w9NDymodHgce+RTw8e57bWSmQWx38T7tvydi53qbEkRBlhNUra89BJ5y/Na2jV64sqK+vqS90OdyBDyOuzstTU8YCVQVQ1ikbEBtNzV6Ehaf+tFhmyTJ2I8xHICttkOa1IwtU+kooP/mQ/nB6RNt3CvL6QKOIhabQDzCSXH3g4dD8QHAAAAAAAVpDYEOXDYL7jHI5kkBNAXWxoltAQJzIkEfWzFBnSCAtJ76i7VF5o9BiS8uJMavGhNuFPJz5ECQ9xqRau0Q5IsWhCWcyJKaLeHnuUw9KSv01GJTRBZwLBAQAAAACgi8UGl6iGQIyISJuQE79WCA1mUSFeaEgjMMTt1wphIU0/cXftGxmTi/gg+pSRBoZ0ippAEE650LcLj7F+ea3vWrSD3DZJtANEh4xZM0K0Yw/R1l2+aSSnT6imkezfML/om0XythYQ4dD9QHAAAAAAAOhyscHF9NG2PJwu2NtysSEubSKJyGA3mGxVGkSy8xftndCXiRARBR/DJj7o55fLWJpEBV148JfVohDix2CvcpE02sFFdJCvO+nnb8V5OnBUwwHra2UxuellMVloiCuLmfOE6JAl3CdoHyA4AAAAAACsQLHBNuFT0yayFBuyFBoiIxLKC5FeL1mJDFlFMUQfIztvnDSvV91HFx/0sdnEB5vw4GoyaRMeXKMd1ONI0UF/babXDdHBARYSRlb5DQALEBwAAAAAALq0EoUqNuiTLxNyMqgb7bZSbEgjNOih/lmKDGmEBT0SIGuiIx+SRlEsNLydTXxwFR6SjLkmHLhFO5iu+2ZFO6CCRXJQFrP7geAAAAAAALACxAZ9Ox1bNZ9GxIZmCg26yMCod+BbITI0W1hIel5dhQi9ikOac6VOyG3+DK7CQ7OjHUzeJC6ig3ydiHawGEIyUWUxAYDg4MPlMPWSmKE/uo7l9NKW3Wu0BGRkCShbfUOdFO7OabYL7+Na6nOx7Wpap/mBZv9Dby+fGfVDrLS/9ry8OFt7PD8R2q6w4ajgcW79wcHjniNOro2hMBDaZ3H7z4zj6WRcXofrZzjLO0tJaGYp2qhjNbJNq0nzneV6jaeZjEWapLWqrDJYMWKDzRzSdI3WT9Lqr8dmRjXIz0acEWQrhYaoSWVSYSFJydE0mM5BEiEiS/FBjwQwiQcuwoPL60gS7RAcpxrtUF+5AqJDQzyxy///sAMa6gamkd0PJCkAAAAAgC4RG/SJrk1siHpu288VVxGuWREN8eNLbqCZlbCQlWmk7bhphAi30qQLTRceTOJDEvFbj3bQRYdg+1xfMA6b6GDyeUCkg8bYcOx7AwADwQEAAAAAYAWKDa2MashSaFgOkaHZwkLSfpohRERFQYSPHZ+CsBzCQ23bWrSDLZWCjyvHoZtJhvsKCxEQHZogOHCFioyrVGTeH2gICA4AAAAAAF0kNsQZQur9JM2fzyp9IonQ4CoymO5SJxmvnIi6kjgiImEqlS3lqxlCRNbig/5eRAkP4jqsjkV9D3RvBvux6l+3jHYwpbjK7flYUnSw91MvOsjlrqzIspkAVIHgAAAAAADQJWUvk4oNUfs0K6oha6HBNOk1361OLzIkOa9ZklX/LkJEUvFBjwpwrfhgjmioLVOjH1yFBz2SobY87O1g2l4eR/o8mK45k7lk0miHrhMdKhX/f0+JJiguEU3PEi2V/PV9vUTDQ0T5nLUbeDh0PxAcAAAAAAC6VGxoRuWJtFENyyE0pBEZshAXmm2qaze5rh+j6S6/qZ+k4kO9CWO970GWwoP/OvyKFvo+/Dw6SiFskiqPob9+ft1sEArRwYFde4kKPURrRv3ns/uJxvfxmSQq5IlYj5if9gWIzeus1SxQFrP7sctNAAAAAAAZctddd9GLXvQiOuCAA8jzPPrqV78au8+dd95JJ598MvX19dGRRx5JN95444p5TzpRbOAJXJTYoK4P1pUXggkur4sTG3jSqk9qeQIrJ7xy8mkaI09e5QRW3c62D59D2aLORVQfzSDJMaNeg60fvrZks22vLzO9Fy7vofGaUJZJ4cF0HanXWlxqkOkc1I29+nptVUlMrytpNRQ9PapjWSwS9VfPHUcz7J0k6u8jOmgj0QEbiA6stlyOaO/Uco8WLCOIcBB/6BapUs7bz1KdImr+o5N1+UDnspFRpS9TlLuMXpf8NWZdvs9mfJTFGKw/oLRjhu4CWMre1SvpqlFX1A8DtY+CtT/1uOof4/K+2mPxfD+rzT6F2b3B457DTqg9flLtsU7WJTJbWeZUxc2FO93rU9/b5Xp9WRNZztGwzXIRlROeto+k23l5t78RdZjOcZeUobUxOztLJ5xwAr3qVa+iF7/4xbHbP/bYY3TOOefQ6173Ovrc5z5Ht99+O1100UW0efNmesELXkDdTJZiQ1Y+DY16NZiEBn3fKEwig8u41GgGt6oH2Qg1zRIekpbwtpXidvGHENEDjpEPzYp4kNEOaqSD/lqyjHZApIMj5TKfLP/x0hJRqexHO6jpExzVwOaSe8Il20PkmmAayX2CtgGCAwAAAABawgtf+ELRXLn++uvp8MMPp/e///3i+THHHEN33303XXPNNV0tOMSJDSouYoO8Y9uoSJjUq8EmNiQRGqJEhqhxxZ23ThAWkh4vqRBhw+aXoIoPdTdiLKKWq89DlPBgEh3kOG2ig9w/ztsh3JdZdFD7VMUUk6eDXLciyOeJikWiAX69Xr2fg4SXVe0ewMoEKRUAAAAAaEvuvfdeOvvss0PLWGjg5TYWFhZoamoq1LpNbJCTPlexoVFzyKiQfVNIu/+8GCk2RIXRy3B7vcqByZ/BFNqvLjN5NcSlF8Sdh1anT7iSdHy29S7nx5RyYU5JMadb2CbleqpFbXktvUJPsbBFaujCQxTy9eqfM/2zpr+WRuj41IqhAaLJGd8gkr0cegtEk9M1M0kZBTE1Q9QXlZLkNaWB9gGCAwAAAADakh07dtDGjRtDy/g5iwj79+837nP11VfT6Oho0A4++GDqFJopNrTSq8Hk16D7NNj7snszqGNSX59pQql6NUT5GNi9G5okLPDxk7QMiHst9uXR5yx4rokPLsJD3ITd5O+g+zqYvB3SpAeqr9e0XVLRIYmnQ0eLDpwqwULDtl1+ygT7OUzPEW3dRbR7n28qyY/Z62H1yHKPFiwjEBwAAAAA0DVcdtllNDk5GbTHH3+cOgHbxCNrscF18uwS1aA+1u80m4QG077h5/aIBte792FTyPQiQyKaKSA0qe8sxAfjdor4YOrLFPFgQ70ebAKXKjzIa9B2LdrGrL9O4zghOhhOlke0cS3R6Cqi/QtEU7P+co544IoVvIxFCK5QweUxbXA0Avs+ZNoQ4dBOwMMBAAAAAG3Jpk2baOfOnaFl/HxkZIQGBgaM+3A1C26dhC421OfA9zqJDbaoiKRE7WssZ2mZ3Ol5/tFeD2GxwTYOm8iQ1pchlbjQrpjG5mBEq19v8cvDxpNyO/09lT4Ipr74eZQfQm3bBc0DQt/X71eKDurnRPV4MJla2s9HIdLXQS5T+7N5OrimXfB3wL0PfYg6UnQYHfYbCw2lkr+cK1P05M2eDmDFAcEBAAAAAG3JGWecQd/85jdDy2677TaxvFuIEht04iIb4vaXRE24kogNtueNCg2mcZhSJurHbvYdsPXR0eJCk0QI0zmKmqyrFR9MAoXJZNI2Yfcf1wsPuhmjKjao/cjjCGEgECbCokMS5GuziQ56nytWdJCwwMAtIc3wXICHQ3sBwYG/DBbnqMKmJgoVmrV++VbaoGpZo2XcnEtuuh4nRanKtMdqZt/u5TMtTtDl+vA7U99qGb2oH4tRP5TUkplqZVT9vVia3aFsNxu67iWFw08J7aOWyawUa3nSizsftL6mRsniGlKJOv/AcL4ifli3QylT1++yNGNt9Nqrc2sv9TmVzwyvq+Yll+apm5mZmaFHHnkkVPbygQceoDVr1tAhhxwi0iG2bt1Kn/nMZ8R6Lof54Q9/mN72treJUpp33HEHffGLX6RbbrmFugGX/G2X8G/jdrzO8Lc+jdhgjmowiwg2scHch9nVP04k0MWGONElkdDQDSJDkteYIAKibnIfTPxrEQHqPqHfPZrwYIpSsEUKJI12MIkO6jizEh1Mr3XFiw5pECkVGUdCIKWirYCHAwAAAABawo9//GM66aSTRGMuvfRS8fiKK64Qz7dv305btmwJtueSmCwucFTDCSecIMpj3nDDDV1bEjPJRDmNQWQjYoNNYJBCtsmrwSY2qJUIkogNNiNIdXv1XDj7MWRs0thRJHjt+jmt92Nw83hQ1+v9+/+7VbGINC81VENxvdnmaibZrOoVAHQbiHAAAAAAQEs466yzqKKWTNO48cYbjfv89Kc/pW4jLpXCZFjnLDbI9UqUQ9LoH2sEg0VsaGZUgzQgdInsQDRDA8jz2mDUgx7xoKcfqNEOdesSpFnYoh30SIfw+nTRDnqkg+5PYfKlUMe6olIrEuLlc6Jl26f97wxoPYhwAAAAAABoo1SKKHd803ZGsUF57iI2JBEY1MeuUQ1xZS7j7ojrVRLiqh9YWcnRDE2MejA9N71noX21Upq2fpNEO6j/i8fWCinJox3012GK1jB/Tu3CSVeWywRAA4IDAAAAAECHiA11k54osUHp09VM0i48FAOhwSY2mPqIQg/PN01Gc71DotleS2KhASQjA+EhLD5EiUrmMp0uaRZR12GU6JBUeIDo0ARyVQ+HLBv3CdoGCA4AAAAAAJ0oNvBEXE2fMKEZ+bliEhvUCZxNbHAh0gcgS6Ghk6IZ+By7tOXC4Vya3n9X4UGPdrBHTbiJDvo1aRMd/OcQHUCYf/qnfyLP8+jNb35zsGx+fp7e8IY30Nq1a2nVqlX0kpe8pK5sM3sQnXPOOTQ4OEgbNmygt771rbS0tBTa5s4776STTz5ZlG8+8sgjjamE1113HR122GHU399Pp59+Ot13332h9S5jaScgOAAAAAAALKPYEHXHN1JskOsKtcchUk5Q48QG07ZJxAb5v+l160KDaR9noaEbaQcRwnJ+o01O9cgFN+EhSbRDfUlKiA4dQdbRDQ1WvfjRj35EH/vYx+j4448PLb/kkkvov/7rv+hLX/oSffe736Vt27bRi1/84mB9qVQSYsPi4iLdc8899O///u9CTLiiaoosKzPxNs9+9rNFhSYWNC666CL67//+72Cbm2++WRgqv/vd76af/OQnwjCZjZJ37drlPJZ2A4IDAAAAAMAyYU0pcBQbgm1yvZmKDWq4ue0ucRKxwTV9wnWfTolmqCzOurXq+YxrbSVCGM55lCBkTpmweHKkjHaI83VQo3OyjHSwiyN2TwfQnqWb//qv/5o+8YlP0OrVq4Plk5OT9MlPfpI+8IEP0HOe8xw65ZRT6NOf/rQQFn7wgx+Ibb71rW/Rr371K/rsZz9LJ554Ir3whS+kK6+8UkQrsAjBXH/99aL6EldcOuaYY+jiiy+mv/zLv6RrrrkmOBYf4zWveQ1deOGFdOyxx4p9OGLiU5/6lPNY2g1UqeDv6fl95OX7I/6Y1UyS2p00NejT1qFPe6ys9k92rIWW9Bf+IzKjrVP+YJaVxyXdPKuY2Ilb/Z1ZXrSf59L8nuDx/PYf1/Yphq/xviedGTzuOeLk2v5T20LbFadq5euWC9t1pC9P5FzepM+P7jQPEpxjyw/Buh+MEWGzy41+najXgxxruTTf8nGB5YtuME1KbBUpbGKDKjTw4wr/EXD8vNjW60KDeJwyhSJSSFFeq54+Ydu/jjYSGFg4aPoxLGkLdZiugWaeK63ChV59QqW+MkWtqkWoSoVWCSK8j6wMYTZkVKtC1FeQ6A0qWNRX2eAxFBJXr2C4ZGu5OB16nfqxbWNc0ZUrGoxIsPZJRFNTU6HFnMbAzQanKXAEwtlnn03/8A//ECy///77qVgsiuWSo48+mg455BC699576elPf7r4/7jjjqONGzcG23Bkwutf/3p68MEHRQlo3kbtQ24jUzdYmOBjXXbZZcH6XC4n9uF9XcfSbkBwAAAAAABYJrEhtCwulDxCbIij/m5ufRnC5RAb/LvDhWRiQzuIDIYJfdYip+sddy9ff9xc7yonQSTzc6zup5SRNKGv04UHXXTQ9wk/rhceTKKD+jgr0UHFVXTgcaFcpjxpTTB5rPZ38MEHhxZzmsLf//3fG3f5whe+IFIYOKVCZ8eOHdTb20tjY2Oh5Swu8Dq5jSo2yPVM3DYsjOzfv5/27dsnUjNM2zz00EPOY2k3IDgAAAAAALQQ06Q6yrdBXx4rNvB+ymQ1aiLMFSfUiV4zxAZbVEOUIGGlWWKDPF+m/g0T//JiOIqxUdKUaYza1zQ+0/nVhSen8+FK9TpMEu2gRg2oooPYRol2MIkOJpKIDsbzkDDKIUp0qN8PokOzefzxx2lkZCR4botu4O3e9KY30W233SaMGkG2wMMBAAAAAGAZq1K4mESGwtUtYoNIp5AEIoIplDsczRAnNqTxazC9RpvYsGxGkCbPA4sfAk/g1dbuSA8O3YvD5glh9Ylo1BNCed+S+DuY/BGirp/a42hfB1PZV7Ools7PQRUd4iKX1PFmnfrbUeRzzWlEQmxQm01w4DQFNmXk6hE9PT2isRnjhz70IfGYowc43WFiYiK0H1eG2LRpk3jM/+uVIuTzTTHb8NgGBgZo3bp1lM/njduofcSNpd1oC8Hhrrvuohe96EV0wAEHiBIkX/3qVyO3//KXv0zPe97zaP369eINOuOMM0LungAAAAAA7UiqihShCAEHsSEC04SKoxzUdX5/9gmaDd3k0TRRrHs9FhEh2KaZQoMDWQgMnN4Q1/IDq+uarNihtrSowkN4eb3I0ArhwUZchI+tfGac6MAsh+ggxyoFiCxEhzghEyTnuc99Lv3iF78QlSNkO/XUU4WBpHxcKBTo9ttvD/Z5+OGHRRlMnoeK9+WMM0QfajUJjpgYGRkR5o9yG7UPuY3sg1Ml2ARS3aZcLovnchteHzeWdqMtUipmZ2dFyY9XvepVTiU9WKBgweGqq64S+SvszMmCxQ9/+ENhyAEAAAAAsJyYJgX2u7KOYoOt/GVCj4FwZEPVsLQ4XRe6nlRs0McbfVfaPPkMPc9SaEgweUyajqCSG1jd0LHj0EUHKRa5or736gTZlFIQTnPorX8tad4fxzQL1VBSjjm0XLlWTT4JMl3BlmIRlV4RWtaAiaQ+ZplqYUqxSFq5oqtMJJtoGunK8PAwPfWpTw0tGxoaorVr1wbLX/3qV4tylWvWrBEiwhvf+EYxwZcmjc9//vOFsPCKV7yC3ve+9wk/hcsvv1wYUcrIite97nX04Q9/mN72treJee8dd9xBX/ziF+mWW24JjsvHuOCCC4TIcdppp9G1114r5spctYIZHR2NHUu70RaCA5cN4eYKn3gVFh6+9rWviXqkEBwAAAAA0M7Y7vzr60PLIsSGuugGY06/79VgEhuCbap55/qEKEuxwckgsklCg24OKJfpuObxG8UFA/mxQ8ILyktO++l32EtTW+vHYIh6MIkQUak14f7qx6F7JtQGVEwvOlT3twkPYTHALDqI51VTybSig/6abKJDbVzxr1f3RDELDNHeE6D94NKVXDHiJS95CS0sLIjqEh/5yEeC9ZwK8Y1vfENUpeDJPwsWLBy8973vDbbhkpgsLlxyySX0wQ9+kA466CC64YYbRF+S8847j3bv3k1XXHGFEC24xOatt94aMpKMG0u74VUqlQq1EZxS8ZWvfIXOPfdc53041OSwww4TahHXM7XBbwg3CTuCsnvpacdeSj1tXhO3Xb6UshzHcuWrLddxVfU6SslW//CFcgAttav1dfoPDdW4SL1ro54HfTx9G08IHvceckptRS6sUc797P8ajwPSk8awqtm0czlPvdRkO+RWN/IeLpUW6EcPf1jU2VaNrkB6+G893xFq5Tl1iW6wV6Qw5K7HRDaEBAdT9IKSMmGa4OkRDbbtTJjTPtogqsEirLi+Lp1c/+o6Yaew9klOYfd126QUHIwYtinN7qpftn8fJcEkZERWEGnkfTP4h6iExS/z+xoVlRNXyjzJNew/d3+ttmgi2+tLSquiHJrxPSr7nPj48TQymM+kz6DvuRKNvfbn+FvaJrRFhEOj/Ou//ivNzMzQ//pf/ytyu6uvvpre8573tGxcAAAAAFhZJEmliMJVbAgRITaoxIkNpu1s46t77CCk2PrJhAST1ziBwQYLDTbyGw4jKoUFhUqpFH9AkwhR10/9NqU9v68fw9CGOuGBPSGSCBDyulGFh8hokAyqWtirOYRTLPyxhEtuqpEOcqy1qAe39IqQ0KZEOtTv555mYYrO4JtJ6ueukWiHrkqtAF1LxwsON910kxAROKViwwb/C9bGZZddJvJd9AgHAAAAAIBWiQ1xk/LgeYZig353VU7QGhEbXAww1deVRnhxQrvL38gd46RCgxAY0hxnveYmXy7XbVPesSW2n/zaw6wChEl40AWINMKDXM7L6qJHkgoP2nsXlWIRVTozztehEdHBnIrTmOigv64VnWLRBh4OoLl0tODwhS98gS666CL60pe+RGeffXbs9mzYYSuHAgAAAACwXL4NKly1wBU9zD9KbFApze8Rx0lS7jJKKIiqtuEsNrhM4HTfgAixwTRRdBUZelbX35DKrde8GJhyffSCt/GA8HODmEAOUQ+59QdGRkqUdz9Rt09+9UGx/aoihBr5ECU+qNeTPKdSdDCeb1M6iOn9rZpI6th8D+RxbKKDqQ9VdPD3X3AWHYyvrUHRQRf8GqXjoxxynt+y7hO0DR0rOHz+858X7p4sOpxzzjnLPRwAAAAArGDiohtc88Kl2OAS3RCIDdqEzVa5QE5wpNcJ/5/kzqpJNJBeQ/qELHFUg6vYYHmexOSyIYFBw9u4iahkEBQUKqtrfkwBZc1CTe/DJFLs3hMe39rN9d3ufrxuWX4kHE1RmtoRRD9IlqaecK6AoZot6hEQkSKPrcJFhOigHk9/rIsOwfGVyXxd1IIl2sEmOphMJNOiCyRRhpcAdBNtITiw/8IjjzwSPH/sscdEzVMu9XHIIYeIVIitW7fSZz7zmSCNgl0/2d3z9NNPFw6ezMDAgDAfAQAAAABYTrFBpam+BRFVJ9RlutggKS/KyhR9qcWGJKU9O0Jo2HR49TiG6IX162OHWxlNkAqTAG8s/Bu3sndv3Ta5NeHoivLuLVYBgoWHYNmQ74Bfmt1Z66sqIiQVHowVLeIiVZTKFTomX4a4lAWT6BDexy46qK8vTnRIWjYzjrSiQydHOVTynmhZ9wnah7YQHH784x/Ts5/97OC59FlgUeHGG2+k7du305YttS/Mj3/847S0tCTqmnKTyO0BAAAAAJaTqOiG+m0biG5wSKWwiQ1Rk60kYoP+OjITGyKqNOiTMlv4fSqhQfaxZi1RPtewwGCa/NQv0Y5jOm4xbBrpjY7VH2uvHgURTstgynu2hoSHpX2PNyw8qGkWtf2q44yKeIhIj9GJigRwFR3CBpQ2Xwc9IiI70SEuykEdKwDdQlsIDmeddRZFVefURYQ777wz0+Pzl2M5XwvHSks7fDksV8nHbhxr1uOLKkOpoqvrteXh6zN8Z8leI9orqSG8i8bxlBenQvss7PyZcpzaD6nC4aeEfVEOOCl4vP9/7jKOGySj0TzqlTIml7G5fobDn7nl/V5yKoMHOi66IUkqhQt6KkWcSaR/B7pebIi7w9sJYkOU0OAZSjs6CQ2GiXxwvMH+uvzwVt1NrRS0n+25cNqFN7efvOFwycLKZL0vQ27MFxXKEzszi3gwRTVYjSVVEn7nNVt00FMrbKJD/etv3M+hUaGhY6Mc4OHQ9bSF4AAAAAAA0C1ig9VYMWK5nNQ4V6aIKP+oL3MRG+LGm5nYkCaFIoHQ4CwysD9DPh8fMcACg74sRmDIjYdF/GZR6dPOcb82GRYChJaGMbGnTngo7XmizutBig9SeFDFB9XrQRcfdOEh0lgyJVmLDjpRokPU60iaXqFX2FBNJFdiagXoXiA4AAAAAABkhD5Jd0mlCHCYrKipFOqERE78QsuqJTDrJ1J2sUFuEyU2yHHrE6yGxYYGhAaT2KCWsTQZQBoFhuHB8IJcfH55bu+M2eQxgsreca0TTfwYs5fnDLbRUizq0jBMldkGw4JRZXKcciNhX4ry1G6L0aQ56iEu4iGRsWRGooN/nPrUBdO+LiaSwfIYE0lX0cFWYWMlplbAw6H7geAAAAAAALBMqRR1E5aqoODJBPgEVSlMYkN9tMNUrEGk3E+8BovYUP+4uWJD3MRLFRt6xg6ppaqYhAb2ZoiLZIgoqycEhqSCQkIq01PRgsSwofJFTLUMsd9QddI/618z3rAvbFSma6kXuSF/WXlWXeafs/LsnjrhwSY6mIQHVXQIjSuhAGGbpPvrwtdSo6KDeswsK1cgtQKsFCA4AAAAAABkQJpUCtvdeRYXTKKDLZUiTmzQ89T949cLD+qyKLHB9DqWQ2zQhYaQWaIpaiCXi02Z0MUGb3q+YYGhMjsZXpDX/RjCzz0t3UOPTKC5ufqDDGrRGTpKFEYgPFSFjUB4mNhdJzyUpnbVCQ+l6e2UG6gKE/v3ORtL6ikW+npXoiIDorZPKjroY2tG5YpmpFZ0FPyRzDWhT9A2QHAAAAAAAGhidEP9es3c1zAB06Md4lIpXMSGuDu4JgGiqWJDihSKYFzVya44xPDm0PjyGw6LTU+o9FfHG1WFYsEyvulpouJiMoGhUfTj5Qfqt1lYiE6pkGKLKjwMatdetRRoZbpWcjM/sqFOeMgPrafSrC9OuAoPqugg9tNFNpdymg6igx45YNveRXRotolks1Ir4OUA2gkIDgAAAAAALYpusKVSRJEklcJVbKgfZ3jipFbNaKbYYDLj019D3bhUcWF4c/iQitggKjZokQKB0KBXfFApV0IlKL2ZWaosaud4LpkxpNc/FBPhoEU06BEOVSEgOL72XEREqMsKBaKiNunmZarwoK8XkQ/hailG4WFiW53oILodWC1EB/HYIjzYUiz0bZolOkRFEdgm9mlNJJc7taJTqOSzr/TCfYL2AYJD1QCmQvaynGlZ7vJq7TKGbhxro+OOKpHpWj7TVmZKdxAPraPaOluoLVPaX/sBMb/9x7W++sPltvJrlfDVbbW+lxb9HyMgW2wlU9uF5SifWT/Rqg8bbzWNHFe9iw26M7pBbhMX3aBTWZwNTUBcTSLV/5OKDaaKFFmLDeHxOogN6iEUsSG/9qBgEu+NrrZXdODJthbVUCnkyCtqJSYnw+cy2HZ6gpwoGFI1GqGgC1cNzqg4ykHvQxEgAuGhXKqL2MgN+8JDeXqXEB0YNdpBig5xppIm0SEcbZC96GBKrdD7yMpEcrlTKzomysHzIj1TUvcJ2gYIDgAAAAAALSiDWZdKkVB8k+Ut5STKxSQySmwwplBoYoNpwhQpNjimUNjSJlzEBlVo6NlweG1cUUKDIaqBhQYdb9rgjcDMTddFFUTh9faFoxb0yb0e0aDv36igoCOjG/RzoVa6kNuowoNWVrM8tSdWeFBTLFxFh3B1CVv1CYORquV6SSM6hPfLxkSyVqKz0HBqRRo6RnQAXQ0EBwAAAACABlArOiSJbnDuP6LsoItJpHkcfalMIrMWG6KW1QkNiikkR9ipk3JvTbi0Y2VwoCGhoTI7Q1TUzmFZK0MZlzYRtW1vzHWQQNwwoqdL6AKGjPJIKDzkSv52ZTXiYWgNlWf91ItAeJjaFrx3LJS5pFioEQHmkpPpjCVdRYcsTSRlf0miHWypFfqYu40Kl57NOMIh6/5AY8DDEwAAAACgwVQKxqUMpm291UBRS53QJ05JTCJrx+6zCwnV19BOYkN+5IA6sSE0LlVs6O+riQ08mVbEhkp/j2hin1KFvLnFoLERZGXvnqCJ7YsLoabjFfpCzZgGIVsuT16hN2ix8KTf1tLA4oFsBv8KrtahNpF6wqaTmvGkNzwm/s8NjYomRRgWHVRyVeFBPO5dJZq8g2/CdH3HXStxk+9wpEDE9WwQCWvXujnt1TQOPVqiPtrI/Bm3Hdv2vBnfXQA0E0Q4AAAAAAC0gLqJgz4JkhMSy+TdlEqhP3fxbZCh4bUoB/fIBiMJxQYX931VbAgdShEbvLG1NQPE/j57RMNgwTeCrAoNdeze424IqZWvDBEnJKhiQV2KhXYPUPWZmNfeS110MJg/hvvSjjW33/9fiQLxiktUUcWZ4cFa5IcUHaqlOKXowH4WuWG/TGZ5ek8gOpSmdvovaWg9lRVTyfzAairt3xdhziijC2ppCLZIhzToEQQhQcJg1JjURFL2Y3pNUajRDysxyqGS81vWfYL2AW8HAAAAAECDRpHR6xNGN0hhwVKFolGTSNNYXAwijWOPIoPIBkl+9QGB2JBbu0G0QGhQxIaQZ4MUG6QxpC42sNCgiA0sNITEhmr6QCRsVCmbFsXAaRNqS418jdXGERxqi42A4PQMtanCgxQfqqIDt1BUyPCgLz6wADIcNu2UwgPDwoMUH/IjG0UTy4fWh6IdWHQQ+8ZEOsQtc8F0jKi0hlp0jx5hUB8BYYp0MH0m6/dTPxNFRURURIZ8ITLKIU3EA6IcwHKCCAcAAAAAgAaJmrBHRjfEpFJwzrt61zOoUpHCJNI2JnXsepnCRic6+hjjJo+BaaU6SZVCA0c0qHfrVZ8GTWyQQoN4rJXc83b5fgNSfqhMjtcPxCI21Jk59oXHEEVleMgtokGsr43Zm9nv3K83PRsfAcGig/o6tIgHb27e75dTK+QxBvv95VJ0mJ4JRAe1cgeLDhztIF5SVXTgiAc12iFJpEOUp4ML+jHCE3uDMaW1ZGa8iaSEP5vqZ8pmJKlj83no9jKZiHDofiA4AAAAAAA0KbpB3ybNpKnm4m/fN5nQUJskRYkN4X0SjtuQ/tGQ2LBmY3iyrIgNzkLD5CzRglI6lyfKuhGkSWhQ0yiqkQxWBgeVbbMJJK6s0oQV5XWx/0Ro3ahWGntmLiwuzPtiQhDpkEB4CAQIjoJY40c3eHsnQikWuuggDjG0lkqze9pGdIhapxtIhreNN5G0kdzwsvaaxXNDRY2koGIFWC4gOCgsV832bhhru42nk15vXJ6tRFfRS/treZF6mSabw7eaW6guLy/a81VL+2t3foq7fqX1Xfth1XfQ04LHSw9/zfo6bJjuEoDO+py5voftdocm6nNm2y4NuMZXTnRDbJnMuOgGGcWgOfrznVP+HtejG5KkUsjrMNcbP2anCZIpRD2h2CArfNSJDWMHhMUGiUFsUIWGUr8/0c8V/RiG/M7a37jKxJ5QecvK4jx5LCKEXlNPfYlL29j7a1EAJqEgJHqkdc6v+k8YjzGohd7PhwWTyqiWBqEbRiopFcZlgwP1EQ9KZYuKJjyoooMY+sROJ9HB5k2Sheigm0faKlbo6FUrTKKDyc/BFOXgPtawn0P9+OvH1OnwZ0QXB7PoE7QP8HAAAAAAAGhSdIM1asDRsT5pjngS0ct1QmR9nRmJDWIs/WMhsaFn7aH++tF1tY05nL8a0i+qKWhiw9JgLhAbRJ+T84HYUNm5TTSVyryfFmAdm6X6hDe0Kmg6ldGhoCUWFWRrhEK+1gwRFmqahHi+ShunnurBhpVV00oWHqT4EKoCUuihykb/fWLRQUQ8VH0tcmsPDDZn0YGvD36vg2UDqyN9HVyrV9hIOinXTR8jK7PEwJ9VNSrBdSzS26FOGIkpudupXg4ypSLrBtoHRDgAAAAAACTENhFxjm6wYItusJXcc41uCKVRJIhuaIXYoBOIDWs2+Qs0s0I5aS4P9wZGkCw2hPoY9+/SixKXWklLXWiolJbCUQ65nrBXQ7WspdU/wVBmMlhfsM98jBUzhHljObxAESH0/mx9GCMqqv0EosP8Ykh08GRkgy46iG2r57C/T4gOunAh+pGiw95J8kZXU2Vyn99dVXQo79lK+ZFNVJraIUSH8nzN+yF4PYY797ZIB3WfurFEXnfuUQ6mda5RDsF69l9oQCyoP373RTmA7gaCAwAAAABAE4jyRHAhEBWqYoM57NucVx6bUhGUwJTCQqGxu7raBC5uIqT2nR/yK07kx2p3w6XY4I2NmsWGwdpPWFVs6JleIlqojmXb9tC+5dlJ+4AK2gS6X/Fi0OHqEIrIUFpl/zntKdpBbl6pEMEUledlZcOSJjj091o9G0LbKtv5B4wIK+d1g8okuFiq93+Q0Qyq6WS1ioXJWFJSWTNaJzqIQ64/hMq7t4REB/ESZneJz4rqVSIjeeREXS+Z6W+XPs0iSnRQ/RJME/0o6g0iq+Ouig5pfSiCMWTg5dBuwDSy+0HACQAAAABAitDj1NENCb0bbNulTaWIEhScxAY9uiFhuLtRbFh7WG392Po6sUGkUMSIDYXJqn/Bzl0hsYGFBpvY4PX2ixY87x8UrW67sdGgiT5H+4Omw94RsnG0gmxpkakMos0vBk0XJtTthDDAwotseqoGix2q4MHpF1oKhiiHychymzLqg/vmpooSBtFBrB9dLZqMEGHRQRxuZFPdNaAKdFIQ01MSZPOfm687l2vQVikiuP4tEQm1z0dfSz2WbN810v8kiYjRjmkVoLtBhAMAAAAAQALcKlMUUt3NlGKDLbqhkRBqdZIUrpxRSG8SaZj0uaZRmMSG3NrNRrFBFxpKfbW79z1zyoR+y9bQ8coTu6oPlkLVJkQKhZJGIVIoFL8Gr38gVF6yPDoYGTGgRjJY4Ym/IhJ4alRDCrzqpD9AEQyEIKFQ0aMmAuPHUnhfm+hQ7c9bMETZyGgHpVoI9xOIDtNz5K1Z66e3yEiHvdsov/ogKu17wt98aENdpAOLDtIHIYh00CIFbGkWacSG4PWETCKTpy20KsqhW1IqEOHQ/SDCAQAAAAAgIXF3QtNGN4SWOYROu3o31KpTaFULDGJDIt8GB7FBTvCMYsPqg+rFhjVrrGJDueDViQ0BithQmd4bFhtCr1mrQlEVGlhkMJlBlteN1L0mHodtouTNFZW2GLRMkJEIppKbLCrIpo9JEwr4uWwinUMXL3Sq6RosKghhweRbofahjEOKFiw6SHLrD657/+MiHQLxrU6EczNgdUGNoJCfbXt0QfIoh0bSIeJEikZEDACaCSIcOpQ05dWaGe6VRbm3di/55/J6m/kaXPtWy2Xq+6lGYao7eX5AcQHnENWZrcb9i1NbQtv1TCk/FJQ7VD1Dfimz2n6PUbNAqcFsyOLa7faykeHys53zfQWSExVyHJs2YSLGKFJiS6WIixywXY+62OBSxjPydfDETn0tDr4NLmJDbsx/Xic2rBoQqQCq2CC2KfsTfVVs8HbtJZk0UJnYbRQa/HGYzSFFZQV97KOrQiaNagUMndxczIRdoggFFZNoIA5kj3wIRUWo+0fsI/0hAtFBEws4GqLCgoIUDGwmmNIjolQW51oVMSrSaLLq8RAaF0c7DA/WIh0mq6Uz1x9M5d2PB9dBaWon5UcOpNLUVmOkg0oj0QIm/waXyJwkUUa2KId2+o6796EPUTtQ9ojKuez7BO0DBAcAAAAAgATE53lbBAbH6AZxl9VSlSK8LF74ClenMEc31J47TIik6GDNny/WTepCd6z7V4fFhjU1gdobHQuLDVVDRhYYpNgQ9FOseRLkto77+0zvIyr5k+ZKUfMWkOkULDgoqRVsDulxFQp5XDZOVA0hB83VKURXc5oJpMQ1VUKfmJseawQT+yTig1p1Qo9mqK6TKRjCDFO91uQY1f655CabTMoUCjkOeVxduFBEB5GasW4d0fh4KNKBRSKuTrK0538Siw7isZJa0YjoEGXSaEphsFWsiELtLym6saWtggUA7QRSKgAAAAAAmhDdoE/wbaHbLtENtecLQZPPzWPqM44lKrrBGQffhkixoWoYGBIbquH27JUQiA3Vyb4qNnAZSFkKMrdrMhAbyru31Y6tig25msCgmkOK50Phc1FZE45yWBot1Ikc+fly0ALjRW48qVbMFEP0FWqNJ+Gy2VC3iUihYPFBthDqPvp++raaOFKXeiFNIWWKhPSxYNFBfX3qcU2vT+4vl7PowI3NJHP5oCqJLInKooMtvaLuPFhSK+K9GtwEN11YjPvM1D6b9UJiO9Eu5pEyNSnrBtoHRDgAAAAAADQruiGhd4NrdINLWcy66hRKCHqq6AbL2G0TPpPYEDw3iA0qUmyo5KtpFFWRQdKz1Q/LL+/ZFUqdqBMbTEIDV6GQFReYYU6d6DEKDULgUObk3kJVXDAhJ94RqMaXLuQmlddjScGoM4/U0aMfpOggxQab6CCjH6QpJKdULJRrr9MkOkjUChg6fK4NY2bRobJ3hxAdShM1ASkq0sFkIqlii3xI4vugfh5dohyS9IkoBz+dIvOUCggObQXeDgAAAAC0jOuuu44OO+ww6u/vp9NPP53uu+8+67Y33ngjeZ4XarxfO5EmukEnSXSDeQx9TtENqapSxOBi3KeKDT2rDw6PRREbymtWGcUGnUBs2BGuSFGZn4kXGwZ9A0hveET4RUjPiNLq/qDpkRSRPg086ZZNodyfN7akqOPSxxi87kJP0GJRIx9ySqs7MEduFM3RDmqpTa3cZvCeqefDcH6C6I31YbFJRDpUfTZ6VivVS3qHRIsyV7V7nLhHPrh+JuqNJPs6MsoBgFYAwQEAAAAALeHmm2+mSy+9lN797nfTT37yEzrhhBPoBS94Ae3aVa0mYGBkZIS2b98etP/5n/9pq3SKNNEN9olRsuiG8LH7tLKX9dENqYwiLZh8G2zRDbrYwGURdaTYUO5TPAqqE3514p8fn6sTG1hoCIkNVXNIVWzw+odqYsNYONqitK4myJhEjtx8KSQ2VPp7ghZ6DQ7CwsJosmZjad1A0FKLD6EXqadlVF8Diw7cpPHk3Lzfiks1AaIqPMhzFxIddOFBokZVKKKDVxUU82sPqRMdoqKL4qo/8DUqW2IxIWFKhSui8kZKAUI/F2nH1A5pFUip6H4gOAAAAACgJXzgAx+g17zmNXThhRfSscceS9dffz0NDg7Spz71Kes+HNWwadOmoG3cuHHZ3624UpitjG6oP/ZI5Fh0ISCLVIr6MdSnUpjEhsC3wSA2BNsYxIbKXt+3QTyem6jfp6CnUAzVohoUsaG8eigQG3iSrIoNuYWyMIWUxpBcrUI2FY7GkE0lqYBgY3HY3JzEh2oJy8DcMTRwQ1qIHu3AooMUHjjiQWV+ISRAcH/evCLK8Pk0nC8hOkifB5X1awOxIbfxoDrRQV6nPas2Rn4GG40g0L1V2jnKIeqz204VMQCAh8MKot1Kzuk/ntptfFG4hrU20ldWfZQXp2L31/8w5QfWG8ts6n0V99XKXXq9fp1tpu+Ak0LbLc1uM47V9XylOa/4YxtP2Nm6cz5/rSxj2Y3nZblYXFyk+++/ny677LJgWS6Xo7PPPpvuvfde634zMzN06KGHUrlcppNPPpmuuuoqespTnmLdfmFhQTTJ1FT8d6Ar9RMMh1KYCaIbXJYluW6l2BAnkDTbtyEQG9YeSFQuVU0Cw2IDT9j1FAZJfnKxlvO/3f97IiIaTGUvFbHBK/QF4fneaDiqobzBn/2X+rzAo0EYQSopAraUDl0YKdb+/FHFEtxQzplfmySXsI5fcbC+v555j4obaqJDYVypnjConBcWCOKqWoQqYOR9wcEkOvRX77DPzFFl1WAgOojoDz6XOV94UM9nEC2ieTlUhofImw4LcCw6lPZto/zYIVSa2GKtINFoqcxYrxTlWI0KBvrY2q1k5nJEOGTtuQDTyPYCEQ4AAAAAaDrj4+NUKpXqIhT4+Y4dO4z7HHXUUSL64Wtf+xp99rOfFaLDmWeeSU888YT1OFdffTWNjo4G7eCDw54BLUunWIboBpNApooNUgwIwrFTTnJcJly62BAak0FssCHEBklVbChP7bEctKc+qmFwVUhsqKweDsSGpcHaz2AhNkTAIoNsUmSQLbRdrmJsJu6+IBe0uy70gmbDK9U3lVKhEjQxxnV9QVNRxQf/xeeSRTvookNVNGDRIRjrXPV9UwWcnME4U4t0YNFBHH7jQX4aTP8Q9Wx+sj+EsUNCUQ5ZRSq4bmsS7tJGObimgrjQDWkVoLtBhAMAAAAA2pIzzjhDNAmLDccccwx97GMfoyuvvNK4D0dQsE+EGuHgKjpk8eM+rhRm1PIsJh+2cUgH/zTCiMRW+tKWrsHlL+UxRXQDj2PdhnRiw/i4m9hQ6CevOjH2hsNlLivrxuqEhjqxgcs+KpNkNZpBFxckcoKvw0JCGqJEB+aPblDSTIq1bVVhQx0Tb1Nc7b/vhX3FkOgQ+DBI0cE12kGFn8tzzpEOQd+LVBnsFedTnkcWHTiihEUHr1j1hdBsFUSkw/yiuFYq476/C4sOS9t/60c6TG+PjHIIrYuJHnAWJWIqVTRKM6IcmjHOZlDJVUTLuk/QPkBwAAAAAEDTWbduHeXzedq5c2doOT9nbwYXCoUCnXTSSfTII49Yt+nr6xMtCS7RDYkmAzGlMGsTjILBQV9PWUgf3WDeLl10Q5xvg0pdCUwHsUFPp5DeCYLtfgRMZXZSmEFWSku1dIpQVENVZBkcDkQHsd/ocGBYWBwN//TNFesnJqZx6WJD1iJDEr57UU1kePbHzWkZqvjAaR4yGkIKD0xhcokqq2opGN7Mfl94MIkOolOZe2LwdVAec1qEjFSQogN7Yqiig0ivKOTF+8zrg4gIWeazOg5ddCjv2Up8pJ7Rg2hp8gmj6BC8bkuKRZJJuIsIqE/s9TKZ/NhfJkt3auOpjjMN9cdOJzLwd+C9D30o1RgAiAMpFQAAAABoOr29vXTKKafQ7bffHizjFAl+rkYxRMEpGb/4xS9o8+bNtBzYohtk1EASs0hJI3cgZWUKWxlM+3gbFxts0Q1SbODohtAx16wLiQ2hdQ5iQwiT2KCVvwzGuc4ve7k03BMSG/guO4sNuuBQ6q/3Z1DFBjVlQUWmRrSa77zWbzoc1cBNihAsOuj+Err4EogPphQLxlQ+U8JihGxV0UE0xaOBRQedoAQqR0LUvYhcTajqG/AbCyVrnyT+Z9FBYjWRVHwXkqRZqPuaiPpMpfUMaui7wPL6O8Ebgv0bmtFA+4AIBwAAAAC0BE51uOCCC+jUU0+l0047ja699lqanZ0VVSuY888/nw488EDhw8C8973vpac//el05JFH0sTEBP3Lv/yLKIt50UUXte87ltAsMqpEXtKJS9zdTVOpyqzKYNrEhvyGWlnDyuiQEBf47racaAalL8v+nW9VbPDG9xIVeqlSXKTKhG9gXClWUwAUQuUvq2kUlTW1dIqlNVXhRZvvqs91oWGpz5yeoLIcAoMNVXR4zkdrj3WfBx0pOnC0QyA6zC2ExAVVNAhFO6jpFWrqhZqCwasm56g8OhiIDnoaTWAY2tdP+X3zoqRncExN/MgdcCSVtz0iRIfinkeF6FCa2RUpEMjlLlEEWaY1hfrVohxMx83K3FX23wnpFEwzBAIIDu0FBAcAAAAAtITzzjuPdu/eTVdccYUwijzxxBPp1ltvDYwkt2zZIipXSPbt2yfKaPK2q1evFhES99xzjyipmRWNplM0Et3gahrZqHdDI8T5NuhChklskNENEhlar0c28GSUxQhe7m2vVUiqTO9zExv6q6EI631TSqY0avHaqEY2qEKDKjKIdZZIhnbnjtdbxIcIX0yOAGGE4DPcT9507VyHBACJFBhsooNkbj/R4IAQHaS3Q89kmZZGC6HzXqj6TZZWG0QHHvvYKFUmJkNds+iwtO/x2pAKw1QuTltfY9SkXu6X5PPSiEeCVXhI4eUQpGlYUktcQVoFaBYQHKqKZyeEHHUbrvmmWZSNbLdye2lKQ0b14bqdrURmfmBt+Hl/+Llt/9J8zbyrNFv7cdhTdZI29b8041Yi03Ze0n5W8RlPf0465S5JFkR9HtXzlTYfOPb4lKw0Xidy8cUXi2bizjvvDD2/5pprRGsnsiiFWffclH8e8x1f/53ZG/nbRq9QkRV6KoUkv/agkNjA1SFC41PEhrrIg/GJ4HF57043sWFwpDYxpdrkVTUqFONUUijKBc9ZbOgEoSFKfHjudW7bc7SB8FUYrho/VoUHIQBIMUF6PCQRHQoF8iZnqDLqC2I9k8WQ6MBpKyHRYWYpVDJTFR04yqGyt5Zq0zNyEC1NPZF68h4lUpiwfefX/BnCHg6uUQ5piRtPuxNVyaWRPkH7AMEBAAAAAEDDNIHPEj26Ic6UUU5W0uCbUyaf5NQbWKbzbVDRw+lZABBRDdWKBSpSbDAhDSJF+ctcNT1j/Xp/v9H+WglGRVSQYoO6zH+uvOZ8bbLSqSKDidvfUHscJz7I90gKD4HokMv5ooNa0UIVHcSykll0KBbrRAdJz4InBB/VKyM/41fSCKpoKKIDCxhyOlnYfAwVt/86EB30KIcog0az0Gf/nDS7MoXN5LJROkV4AN1L93yTAgAAAABkXH/e5Yd/mlKYwfoUIdC26Ibg/wzEEdcJih7ZoEY3BH1VoxtsJTBVsSG300+fKE/sCUwhK/Mz4e0L/TWxQUxCV4smjrFuVSAqqMKCjKiwiQ2qqeJyGUC2UnzgJlJXtKYSVJWoRjuIx7mcaAIpPPDzYFm+fpkqOvB7MTkjohd6xveHRAfV1HJ+Q/W9r6Zh+PvXoh5ym8KRlKKPkYNCZV9tXgbCLyXB5y6p2aR/7GwFynb6Tswa9jXN3DSy+4MFO4ru/TYFAAAAAEiJaYJhq06RtN847wY9HNuWXhHtLZHNHVLbREs3oKyLblgffm4SG/QKEVJsCB1/cd4sNhTqJ3Tl1WZTTFNkAwsNUmxQ0yfuunDlzFS+/bf1r1UXH4TokPOoMlornxnCVNFCKUlaB4sOVeGBUUUHRq+kEYgOMpqCKRRCogNHOeQGRkUTh+9fU9+HUqVCv6b1z0rUtnHEfe7kZzk+oildRIIuNiKdFLQDEBwAAAAAACxIM7bY6AaLf0MUasm+tPDxeVJhi26QokiSiUczohuC8eheCmX/fzYVpH5/7OXdOyIjG4zlNjf4k02m1OcFx5LH47QNk9CwUsWGSNGhHG7y3LHoIISHfE54OgSokQ6OpTS9yemQ6MDXQ99kWHTgKIdKX140W38sOnhDtfe+sO4PwsdJkR6R5vOYVqDIJuIp/rPdzsKD9HDIuoH2AR4OAAAAAABN+MGeZPIR5d3gH79P+79+PFn7TDQa3SBz9et8Gxak6lA/2a1M7osVG4JUitHVkWJD0KciNgTbVYWGlSgyRIkOz7s2PFGTZUxl9RCmsqqPvJmF+vKVuqdDlJGkFB36qtfsaFhMYtGhVChT32T1vVs3TLnx6bDooIgPPQceQ0tbfx085yiH0vxe53PAn6fyYviaaxSbeaRO1p4N9eOAhwNYXhDhAAAAAIAVR6O5yhw5kPSHvMu2WVRmqnk6hKMbXCY1WUY3MOX+fDD5ryvLqEY3qGNQxYayVqlAig1r1gcRFNIDYmkwZy3/CLHBjdveXO97YfJ3YNFB/K9HOnCTy6JSKxgpSPCmO/0qWBzloN6hnt2oVBVZF46WCcZZFZ7UKAeXMrV+etNM0PR1zUJPq8iaRoXHVvs4ZO7fUG2gfUCEg/KjAbQn7V66z2YMFPWjMU1ZTX27NMe1lcjU+8pXw1qjSmTqlPfXcm9LufB7VhipmTuV9tdKaS7XtZL1XcBux/V8NVL/Oyts5SpNObpJ+7NdX1HXXeL838B7HXQcjqkU0r/B5fOifmerVSqSCAhpf9/Ul/IsRkY3BGUwteiGYBxSdNCqUXhzi8Gd8srOWqlDVWyoM4kcHvOPsc7/n4HYkK3owJEOalSDEIe0iBRjpIOERQcZ8cCeDWqUgyI00Nwc0eBgIDqU1/B1U41SKXFkRUWIDn3TvngRXM0cGbGwUBflUJ7wK5uUpraJ3x7FqS1GYTDtb0ZbNEKWUQRJq1V0egUKYRqZcaARTCPbC+g/AAAAAAAa8oc+l9lLg+skIIu8bxdhMJi4pDC6jItusIkN5cHwfa3cfKnuTrmkstcPfy9PaaJ0tVKFFBuC7deMBhESxVH7/TNENjQW6SDSUarnWfyvNhnpwJUr+iImxlWTR5ufgxAdFDPJwV21a4RFh2J/TcxYPKSWPuOPqSz28wZXEfXVjC0Lm59ad8Mj7WezEX8G12oVrRANcHMVLBcQHAAAAAAA2uCHukvEW1x0g55OYYTXGdZHTXq4vzjvBlNkg7EvGao/p911np0kL99Tn0pR6KdK9a64jG6QRIkNKvBsaEB0KFimC4rwILbTRQe1goUUHWJSLWSlEhYdcmVPNGZmveGz0ddHlcXwNZTbeFjdZiw6pBEObZ/HZGlUC8uSVtFJYoOfApG1aeRyvyqggrcDAAAAACuKrPwbksJpCSL0OmE6ReRYcn3W6Ias0imSRDcEffTlQ9EMHN0QGtv8Uu2O985d4YNIsaEa3RAw6E8aKxv9yhRx6GUvYRDZmOigezkEJpL9EaKPXjZTvue66KCUymTRQRWj+qd7Q6IDRzlU+ntFE8eXYpTWJ0c55PrHRHNPg1sI2nLSaAWbTvNxAN0NBAcAAAAAgJQT9UbuHCYJ03YJzY71GYmIekg6sXGNbrClUATr9074+07stqRShPusrBkjqvoFFFeHX49qFgmxIXvRgQ1AxXlWhAdddGgotYJFB06RqHo9DGyZr9uERYe5NYpfyvp6nymOcvAGx0QLLY+IcrCJDK5VJrKIVEgtLrRx9IILlSYYRnKfoH3A2wEAAAAAEPFDPulEIM3EwWVik+sddppgNGKGbRs792eNbqhWjAj66NPKYGrRDVQMP6/MTou0CfG4WJtker39NeNILZViaRRiw7KJDloqRazokCK1QpTMVOAoB7V6xeTR4agKjnKozO8nKrLRYq3fngOfGvPZThfN4OrrECUULkdahdoAaBUQHAAAAAAAHIn0RsiIJOkUsdvYJhaOVTWcoxtWDfpVCeLgiWp1AuqN7w33UTKnUnhDw9a72XXjqP6yRRpFc0QHFhcCQSlCdEidWqFUsGDRYfA3U1SY80SzYopyGDNfK2m9HPTWCXRKRa6yV2lKA+0DBAcAAAAArBiyyE1OUpZSL4dZLk4nvqvJd0m5RZs6NpbakVl0AxsMqhNRrQymfge7PGGuShFsN1w97prVgaCxtG4gvE05LDbIu+DwbMieb73FP8lCdOD3Q4lWUVMtjKkVNtFBh0UH2RRGt/tRMBKOciiP9osmKFSjFjjSQYGjHAaO+CPKD6wWLWmqUqMCg2ulikZ8HJoRsQAfB5AVbta+XQ5Ci7oD1y/bTqpXrP6RcwnLM20X55DMlBenrOeyZ6h2J0t3KJc/oP3+anfLyvP7wv0p+8mQYH//cMhkmvfWmrPsqOwjrLDzPoPOx8nZj1PJN/gD0vX6ctxO/rD0CHdllptGvxOalU5RJ3ZYjSIdIjAyjm4wUq7UxIZc9e60IkQwlekJEf7uh8PPhNZxOoU3OOJPOpW74KXV1dSLnFlsQGRDa0SH57+/TJXBAnlzXM6S36Nc7b0u5MUyFh28Be3zILYrhEwixfuriQsBc3NUeIKoeNBoIDpMb1yghQG/32GyV0XhKIfKnH9dJflmrf9N1Gf8zLr+LnP93hCGstX/W4E83nJTyvkt6z5B+4C3AwAAAACgA8XIuHxsl9egTjhM25uWWaMbbOH0mtDg7Z0M76eJDXVjGPX9G0rrBp3EBtA6WHTwT345/F6z6JDPUWUwHJVQhy21QqPwxCT1TC+JprLrydq1w1EOxQUhZnFT6d18vL/JyCGh5apw4HKTxkaalItOSc9oJlkbRsoG2ge8HQAAAAAATikU6fwb/LuWBf//amSCq6CRJBy72f4StuiGyqhDXrwmOkQRRDeou68ZDFWhCI5t+CX73Ysi8v1BpqkVIVh0MPh4sOhQUdMnOMpBNY+0pVZI5uaIFsIT87794c9P+cgDyRsaEi08Jl+gyK09sGmiYrN8HdQopuUqjwlAFkBwAAAAAMCKoF1zkl1FhVzvqmSpOqZJSoJ0ish0Ni26ge9oyxz+qFKY3vScmGxWpv2w9fLsZKTY4K0Jm/+VC76YoIsPMroBYsMy+DlwlINSzjLwdZCpNFUiRQdLlINIt6mmW3i79taJDkuFJdp27ARNb65FPnjr7Wk/HOWQH9pA/ZtPtRpIun4esxYa4rxUutXHAaaR3Q8EBwAAAACAlGaQWdO0EOsGoh94gtaz+mBjBIUtusGbD4e+m6jMz5KXt9uJeWOra5NYfV9ENrRvpIMUHQxRLUJ0sERCuKRWMKsfyYsymdxUJk62+DkoUQ5ev3mbOD8GbuztkKXIUCuLuXxpFZ2QOgY6H5hGAgAAAAAY4EiC5fxBLidBLndcszCLdLmzyv4NddEN1uNVJ5WyDCZHN/Bx9mqVKZToBhPlNf4ksdSn3zGvHgbeDctOZVUfeTMLYYFIrUrR30s0b7i+ClUDSRYiDOJSYCq6yGlJefK27yYaXkW00b8m1m4doz0HTtTvt34TlXdsJSrO+8+HfNNJnXz/WirN11+P/JlTzSPTVpow0S4CQzsYRjIwjWw9X//61xPv87znPY8GBsIVglyB4AAAAAAA4DihbzQKIkmpu0bSKZzG4lBmk6MbrPsPD/mTRD0HXymVGDC/SJVCD3lFc+SDiHQoLfkTw3KJKBe+062LDcFyiA3tJTpM7q8vfxlHQUvJsIgOKoOPztDck+qjFTjKYWRL9fgsOBjo2XAkLe16xFp+Vq9SkQVxVS2M+5QXGvq8y+8aAHTOPfdcSoLnefTb3/6WjjjiCEoDBAf+I15YRbl8f9sofaC5hEx4uvDLWFfOXdyX9T+u6nZqjmPfASeFt1ucU/ap3Tkr7vtd+AD79xn7Uz9zWSj+UX+Y00wOmm3AtlKxnVf1Gmq0L32d3neWd3myuAMu+8jlcM11Mo1eS2nKY/qPC9HbOHyXuYw9N7bRvIIrEsjweS1vX598VqanRSWByuR49bgWAWLY93EorwubR9qAd8Pyl8kMRbbw+87iUlCFIhdEOViFJxntYCmTyaIDiw86HOUwtW6OFgb3i+cjW1YrffYHUQ46veuOonJxlnpGD6LivkfFsvLiTKRgECUcuPzeStOvK506j+Erp5yxz6slYQco7NixgzZssAvKKsPDDsbAEcDDAQAAAAAgI/QylY0LWsknIUnLYZqiKbiPqOgGK6boBusYon0eKquHIw0oYRTZhgaSo9WQa+nToAsLLDpEoRtJmpieIZpfoMEH91D/eFk0V0QETb6HejYfba26kuuNF7lUQUFGRiTZZyUJCqD9uOCCCxKlR/zN3/wNjYy4ib8mEOEAAAAAANDkcnLlxdnU6RRpSBLdEAdPzKzRDX0R/VtC62V0g0QaR3r9fklDb1S5Q81z10G7kSAiG9o40kEiRIee0PVQF+WgezgkjHJg+uYGgiiHEIV+Kk/spMqivy434uhBUnfsmmEkixJZ+jrUH4s9K3obTqvoBMrsI5rLvk9g59Of/jQl4aMf/Sg1At4OAAAAAIAEZGUkWSuDl215PSNKWlEjgoc0jBT+Dc7HLhPxBHNmNr7/4dHQxHNpGPfGOo0gykFHq0zBooNTfwbhQaReWHj85H00/uQFGv+LIyg3tlY0G4W1TzIuZ0EhP1ATJpr5GVWrVazEKAZOp2hGA2b2799PW7fW+5s8+OCD1CwgOAAAAACg63nuKW+jTsI1TDsLUk9yeMLI3g2cRmErc6gzN01eoc83iFRKYsrohmBMo8ORYgPMIttfdPBMJpB8nbimVTiUyJRseKifRncO0obH1lKpx+IL0msPIS+sfhIVRg6hvnXHipQiadIaRRpzyeWsUJEUKbI8/cn/33IPBWTEf/7nf9KTn/xkOuecc+j444+nH/7wh8G6V7ziFdQsIDgAAAAAACgT72akUzSCCK3mEOvquJpZqtNlomWFQ+Qjqg2YUEUHG7L8pQTeDe3v5ZDq2tCrncRRLFJu+x7KTc6JlhRvaCyIcsgPpUuzsPad76Nc77Bo8nm70WjFnawoec1poJ5/+Id/oPvvv58eeOABkVbx6le/mm666SaxrlKx++U0CuLUAAAAAAAif5CHS2JmSdo7ns2uopPrH6v5N6w90F8YY/Jo9W/gxxEREIF3w/CYn7cfcWcbkQ2dQ3n1EOX21cS2isHTI6pUqo0oHwcb+QOeRKVtj4ooh9KeLUHlovyQPd1CrB9YT6X9u50/w/kBvz8WB8vF6URjBKDZFItF2rjR9+M55ZRT6K677qK/+Iu/oEceeUSUvmwWEByqf7Rly6pUGwDNRndIdtnOZbm+bnHfQ8Hj/kPPDG1XOLb2vLK39ge5NLsztN3ivt/W/TFuF9L8aF/uuwGAWlY+s9XXl0fuDv9gebB9/tsl91qvkiG8GzISJwL/hsGB5PvOzFFlelJUCIgTLjidojTaKypUVPLmH8Ewi+xA80gJCwwiHae6TdTXbj5Pns3HIYKZ1VOUK+fo9+fm6eCfVaMXtvmlLxuFfyOxCMG/vTitgv0eXCpbNBsZCRU8z/V1RPl3NnjM2uQRppFmuAzmz3/+c5FOwaxZs4Zuu+02UbWClzcLpFQAAAAAAGSAOtGXj1nEWnaRMiPDyLrwd1n2MEEKhRhDccEXHly2xS/VrkFEMshymZHXVSl5egURDfx+hvp3LdHh319LI3v8CJ1G4BsknBJRGD2iLQQFABrlP/7jP4TooNLb20uf//zn6bvf/S41C3yNAwAAAAA0Gb0kJv/fDgZymQkQ6mRSihEuFPr9u9eDw7F3sZFO0ZlpFYIIoaESZyJp2kdeK/z/3H7/+kuYmqGS6x+h/MgBlBtaT32bTxEGktxajVqpIvi/A6IUOt3Dgcs+8l3/kZER0c444wz6f//v/wXr5+fn6Q1veAOtXbuWVq1aRS95yUto585wNO+WLVuEGePg4KCY1L/1rW+lpaXwNXnnnXfSySefTH19fXTkkUfSjTfeWDeW6667jg477DDq7++n008/ne67777Qepex2DjooINo06ZNxnXPeMYzqFkgpQIAAADoUr7+9a8n3ud5z3seDQwkD5kHjaGnxi2HyZw0jORJlyCjFCNvbt594+FVVB7tJ69sj25AOkVnpFU879qKe1pFDOzZkDStIrq/ApVmd1N5fiJYlrVxZPTx+0KCIz9uR2PJlQJPxP/pn/5JVHBg88R///d/pz//8z+nn/70p/SUpzyFLrnkErrlllvoS1/6Eo2OjtLFF19ML37xi+n73/++2L9UKgmxgSfz99xzD23fvp3OP/98KhQKdNVVV4ltHnvsMbHN6173Ovrc5z5Ht99+O1100UW0efNmesELXiC2ufnmm+nSSy+l66+/XogN1157rVj38MMPB5EJcWNJA4sYnFKxa9cuKmufyT/7sz+jRoHgAAAAAHQp5557bqLt2TTqt7/9LR1xxBFNGxNIhkzHyBWGFY+G5hpGSvKrD6LcmG8w5o2ubqivyvQUUaGXKtP73PfJe5HCA+ggqsahnFbBRpFOcKQMT35YbFhMH4nz+Am7qbDQR3T0cXTgd8rEdpOLD9xBzYSjEvgz2zN6UEi8K05t8Z/u39PU43cSZc9vWfeZhBe96EWh5//4j/8ooh5+8IMfCDHik5/8pKjm8JznPEes5woPxxxzjFj/9Kc/nb71rW/Rr371K/r2t78tTBlPPPFEuvLKK+ntb387/f3f/71IW2AR4fDDD6f3v//9og/e/+6776ZrrrkmEBw+8IEP0Gte8xq68MILxXPeh8WFT33qU/SOd7yDJicnY8eSlFtvvVWII+Pj48bfBCymtExwwF0SAAAAoPPYsWNHXc6mjeFhf1ILOgOTN0Sz/SKEYWS+OhGUIfKG6gOxsHFkoV889IaGicqctx9fnQLRDZ3DbW/2RJRDad0qyo/PRG7LaRXefFVU4AoUKSc53uSMECrW/Yhozbp1tPW4SWoWHJVQLs0HJpLs88D/c2ULpmfoAOpk1CiM55z0pqYdp8ST2owrJMj+pqamQss5lYFb5L6lkogemJ2dFakVXEaSqzucffbZwTZHH300HXLIIXTvvfeKST7/f9xxxwUVIBgWEV7/+tfTgw8+SCeddJLYRu1DbvPmN79ZPF5cXBTHuuyyy4L1uVxO7MP7Mi5jScob3/hGeulLX0pXXHFFaPxZ4iw44C4JAAAA0Fmw83SS9Ii/+Zu/EfmrKxVZ4UHenew2mlZJQ83Nl2JECtM/jwUvcRe8TJUCwhpWBDKKIQWh8pjFYu2aS3Ht2ehdfaRS5eVAWtz7azERLy02T8hYjsoWeoWLbuHggw8OPX/3u98tIg5M/OIXvxACA6cXsDfCV77yFTr22GPpgQceEBEKY2NhI1KenLOgz/D/+mRdPt8Rsw2LIvv376d9+/YJscO0zUMPPRT0ETeWpLD/A6dxNEtsSJxS0a13SbycVkIKdDVq2btued9djcds+YH6/urz4szjyvJwPu85//o3weP/+M+vBY8HHw2H3pa2mWtYu+Yrhso8Ob5nUdvZwpFbeT10y7XXKGne22aijyHLH2Aur8+rpPvhbYNDLZPAIaTA/H1RV25SuT7a4dpNgzr+zCcbMuphfpEqPAHk6AWOYnCkXPDvUCKdoovga6K4RN7CovBviHB4qMPr7W0orcIF9nSoLM5Srn81eYUh8noHxXLV56Gd0T/Dqtmkvq6dvrOamVLx+OOPh0T0qOiGo446SogLnLbwn//5n0Kwb2blhnbhL//yL4WZ5ZOe9KTlFxxwlwQAAAAAwEzwA74qUPDERSLuIOb6qFycbtvTpxtG5kc2Bv4NWVGZnSbK9RAV55OPr+zRd16b6XBAC9MqEiPTKhqIgDCx9dk5WrV3iEb2Hi+eL/zmLvF/ef++tpuIM0H0QXlBfIeAZMiqEy5w5ABXjmBOOeUU+tGPfkQf/OAH6bzzzhPpDhMTE6HIAo4MkBUf+H+9moSsHLFJ2UavJsHPeXwciZjP50UzbaP2ETeWpHz4wx8WKRXf+973RFoIG12q/O3f/i21THDAXRIAAACgc3nve98buZ7zN0Fr6LTQZW941J/0DQ6QVy4nujMd23f/oJhYemOrxTEqq4fJK2V5BNCucKSDaiAZ8nFoBE6vyOcpNz5NB38nR+V1ftT14yfupXZh6MnniP8Xtv1U/L80u80oNKwEyk3wcOA+G+6jXKaFhQUhPvAknKtKcAlKhqtGcBlMTsFg+H82muQqDzIb4LbbbhNiwrHHHhts881vfjN0DN5G9sGCBx+LjyOtDHgM/JwrUTAuY0nK5z//eWF6yWU4OdKBjSIl/LilgoMKfrQAAAAAnQXno6qw8RSX6erp6RGhlBAcOqNChb5uuWEBQqXCYfOuVQhikOkUFbuXJOgQAuNIvjbmF2o+C1kUXOFICI6I4GtRuzu7HLBx5MABZ1J5cZaWJp+g/ICfZsqpGr1rj6H8yCaqFPcv9zCBAhs1vvCFLxTmi9PT06IKBE++//u//1uUnnz1q18tfA7WrFkjRAQ2WuQJvjRpfP7zny+EhVe84hX0vve9T9gQXH755fSGN7whSOPgcpgcTfC2t72NXvWqV9Edd9xBX/ziF0UVCgkfg7MKTj31VDrttNNEWUw2r5RVK1zGkpT/83/+D73nPe8RVTDYpLIZpPqLgB8tAAAAQGfB9cR12Kzqla98Jf3FX/zFsowJZE9b3RllM8lCD3n7Jvn2nV/esLhAlO/xq1QAEEUD1SqypnfdMZQf3kwDhz8nEAsqi3NUmvXD3wsb/LvYpentoRQNEE/J81uWJO2PIxO4NOT27dvFpP74448XYsPznvc8sZ5LV/JknKMKOOqBq0t85CMfCfbnVIhvfOMboioFT/6HhoaEcKDepD/88MOFuHDJJZeIVA0ut3nDDTcEJTEZTt/YvXu3uAHAogWX1+SylaqhY9xYksIpGnzcZokNqQUH/GgBAAAAOh++O8J3NrgGOd+ZASuHOsPLfIFy/WHn86bCokNZER1KJaqs4dSNClEu49kHaFs4hYJTKUJk7NsQRd8fPIsqiwu0+Ojd4nlp/z4RmSDJr9pA+aqvCehePvnJT0au53SD6667TjQbhx56aF3KhM5ZZ51lnEercPqETKFIO5YksDBy88030zvf+U5qFtnEvOFHCwAAANCRsCM3NwAk+bED60/GoF9eVVScMKRSZElpMC98HCqK8ADDyM42jjz7Q5V6n4Vmpz9Ur9HcrkkhaB0y56cm7Txxgfb/2QG0dsuQeJ6/75fNHQdwqFKRtYcDTrorXIqT00A4ooMjO3TTyA984APUNoIDgx8tAAAAQHvyoQ99KPS8UqmI8NH/+I//ELmrKxEuASzL85pKXuZ6h0LrOs3ssdUIMUKkUVQX5HPkTXLefsGfYALQAKbSmBU2HJWpF/w/qP9cGr631AgnWRZX/Q6sre9zLr0OOpNf/OIXdNJJJ4nHv/xlc8S3VIJDt/1o4RrbbZPvqNCNP2yiznOrXq+sqd4MKqXaD6qsr6ksaiirfzTkj2x9eeT+81Oh508+5Lbg8a1Pr5Xjeenw50LblRenjMfN9Y4kHmvWuJ7HdvyO6Aa64bx2ync1532qcL7m+vXrRTglG2aBbMphCkpFsby8OCP+5qjXiHy80n7EV2ZZeOgjYg+HKl7/EFG5VKuCAVY8olJFcSl8PbQwxWK5GDrqz8X/C0/8SPxvK6Eb5dEihYNOo9SEKhVZ99fNfOc732n6MVIJDvjRAgAAAHQWXJECdAa6yFopLyiVKpbfhb9pcHRElXLBnzCgQkX3sLRhFfXsmiHq7yOa218TEeT7zhVOOoCe5/smu3/3v/2qAN8d+Yn4/y1/MUx//Mod4vHqb/5O/L/42I/F/0tTT9TExrL/P3tD9B1wEuWGN1B5tn1KdraaEnmiZd0ncOPqq68WppRcOUPlU5/6lDCwfPvb306Nkkv7o0Vtjz76KP3gBz+gq666ioaH/T+ISbjrrruEYdUBBxwg6n1+9atfjd2HS5WcfPLJotTIkUceSTfeeGOalwIAAACsKH71q18J1+uvf/3rodYq2OjqsMMOE8ZXp59+Ot13332R23/pS1+io48+Wmx/3HHHxZpyATNRdz69IHWkEBY0cnk/+kBOCktl4d3QTP8GAJoCX7dcLWXnbtr03f2iLfVVRPOOeBL1Pe8C6n3SM32TyFUbgt1KM7toceeDVJraRqWJLXhzQNfxsY99TPyN1XnKU55C119/fSbHyNTDIS1cX/SEE04QysqLX/zi2O1Z5DjnnHNEPdPPfe5zdPvtt9NFF11EmzdvDpUWAQAAAIDP7373O1H+kvM1WdzndEiGH0vjqGbDTthcP5x/xLDYwDXG+e/2ww8/TBs21H7kS+655x562cteJu7A/Omf/qmojX7uuefST37yE3rqU5+Kt7aNIhIYeU+xInwcqvn0hepPTZlfn9OqUwBggq8bTq2wUS4JQawOaUYpS/xBHGt72DAye9NIRDi4wuU3eQ6twymXbJmw7IID3yXZsmWLqN+p8md/9meJ+mHfhyTeD/xDhWuZvv/97xfPjznmGLr77rtFqgcEBwAAAKCeN73pTeJvJ4v0/D9HFuzZs4fe8pa30L/+67+25JSx2/VrXvMauvDCC4O/51yXnEM33/GOd9Rtz7XK//iP/5je+ta3iudXXnkl3XbbbfThD384szsvoA2QYle5TOXRfvIQQNF1VPKGCWArKlUAACI5+OCD6fvf/774XaDCyzj7YNkEh+W+S3LvvffS2WefHVrGQsOb3/zmph4XAAAA6FT4b+cdd9xB69atE4aR3J75zGeK6IG//du/ja0N3ih8c+L+++8PGVTyGPjvOY/NNmaOiND/3rukXgIAAGh/SpQTLes+gRt8E4Dn0MVikZ7znOeIZXxj4m1ve5u4IbFsgsNy3yXh0A82t1Dh51NTU7R//34aGPBrRessLCyIJuHtAQAAgJUA3wyQPkssOmzbto2OOuooOvTQQ0VKQ7MZHx8XYzD9/X7ooYcS/b3n5Tbwtx4AAABwgyMIeR7/v//3/w6yFtgzic0is6pg1dOJd0nSwuN7z3vesyxl29KUqemGcnEr4fWq762re3gzy2emQS9PqZaxpNJ88NDrD28nU3SZ93+lVsKpeNKvwv2pfah9R7iyZ10KM42ze6r3ppsd5JuN8rloF9rh85kV7Hnws5/9TNwoYP+E973vfdTb20sf//jH6YgjjqBuIe3ferUufSN4uV7y4z7VvmabVrZOHW9LS+Jp/g2hZSJv3v8DUal6OHgFWT40TxXp4VDy/69UI2O5HKJXCs5ek18AWFZkNHQuF5iQVlpRpaJcovJuX7Dsq/5Ph/qh5N6hh1Nv9dpc2u6LoKXZXf7YZIWJ4n7yHtsmHl/yvf8W/3/ooHXi//t+MEWv+KJf3vKWH3zW73NLIdi/NL/HP96sf9z8jN93YXGO8iN+afHB4/z9S3t8g8qFbT8VVWri0D/7Lvu0C+y3kHUZS3g4uMMZCv/8z/9M73rXu+jXv/61uHH/5Cc/WRRmyIpcVndJmFbdJdm0aRPt3LkztIyfj4yMWKMbGFZpJicng/b44483fawAAABAO3D55ZdTufrD/r3vfa8wYP7DP/xDUfXhQx/6UNOPz78X8vm88e83/11P8vfetj2Dv/UAAABANFdccYVIc5SsWrWKnva0p4mbE1mKDakjHJb7LskZZ5xRVxaLTaR4eRR88rI+gQAAAEAnoJoqczlpTmPYu3cvrV69OvBgaib8O+GUU04R6ZhcaYJhAYSfX3zxxcZ9+O86r1c9muL+3uNvPQAAdA5lyomWdZ8gmieeeEIUbeC/zS960YtE0YfnPve54nnW5NrhLsnMzAw98MADojHcHz/mChjybsX5558fbM/lMNm4ks0s+AfTRz7yEfriF79Il1xySZqXAwAAAKxI1qxZ0xKxQcIGkJ/4xCfo3//930Xo5utf/3pRGltWreC/9WrOKHtG3XrrraIqFf+9//u//3v68Y9/bBUoQJvAqRVz1ZBuTqngxmUxReshr3/Iul9ufJqoXBGVKlCtonuo5GTLiSaIqlARVxITdA2cTtGMBqLh6lDsh/T5z39eZC6wsM+RiC95yUvoM5/5jLghsawRDlnfJeEfD89+9rOD59KR+oILLqAbb7xR1ACV4gPDkRVcRosFBi6ZddBBB9ENN9yAkpgAAACAws9//nMRlcheSy48+OCDwkiyp6ehqtlWzjvvPNq9e7cI5eQfOieeeKIQFKQxJP+tV8d65pln0k033SRudLzzne8UeaVcoYJfE8iOyuIseb1mEaAyO0PecNi/JwTfgHK8vgDIlKrwUPNMydeLGPLaZKFjzvdSoapPQ2V2Qvxf+vHvKDe8gXJHHYc3CKwocrmcCBrgxhkLfCPgv/7rv+hjH/sYvfa1r6XTTjtNRD687GUvowMPPDD1cXqy+tHCd0nS/mg566yzgtKaJlh0MO3TruaUAAAAQDtw0kkniYn9+vXrnbbnVAWOMGxmeiRHJ9giFO688866ZS996UtFA/FEGV42w7CyFeSKRGV48XYFhe3VCX+bUn74F5Qb2yAe59ceFiyXxpHlxVkqj/+GSlNbxfMND66m3NBauuop/0WVed8Q2/tTX5wrHL0ML6BDQVnM9uGYY44RjbMIdu3aJcSHr3/962Ld3/3d3zVfcGjHHy0AAAAAsMNiPjtPDw4OOp0mWRILNEcIYOf4XO8qMfGXAkClVO8mz8v0Sj3l4jR5uT6xT84SjdAR5KrljXI9QXUKmp4gb2i4Wt2CRKWKnuklWhpuTqQNWF68yWq1qur7zak0oppJqRxUrOh0ln77U7+iRe8gFcd/Iz6zskqFrARW2j8uqoQtzW6nwsghYllhw7HUs/nJ1HPgMVRZXKDB9UfQ4qN3+33OhA10AWgmGzZsoFe/+tWiNYrzN3k3/2jhEoUVL9/U8meN9t2JdyVWCq7vrWv5TNeSmWrJo4bLuDmWoCzPVks6Vfnbj40Hj6/5w1qq1T1/4odHm1B/XKvlN/XSnGngH+PLAkphts95bGZpTXV8GZS2bcX3+rOe9axE1aP4ZkFUtaeVQFAeM9cn7mj6wkDROsmX4gGXwwyVxiwvijQFTlfoerRJYlDmkP/J+am2lQ1ryNtXK58MQNLrq2KaW0gBawXg/Huz+vtQ/Ruj/71pp5LPiHBYXqSVgQ7bJPT394tURk6r0LMZmiI44EcLAAAA0FmYUhSAXQR1FV9XJMWinx/PhpDLwLM/TvSd1y7LoUEG7x01WRfwopz1+doFzsJ3WKioj8AC3cVPf/pT+slPfkKlUknYITC/+c1vRBnro48+WhRnYFHi7rvvpmOPPba5ggN+tAAAAAAAdBcyQoMpze6m/NB6Kk/vEiZ6LaVUosrEPvLGVotIiMrq4dYeHzSfvOIDxwaPDItYrU6jyPVQeWJ7KHKzuNWvlNc39vzWjgXQkpcXLUuWPLs3IAjz53/+5yJ64dOf/jSNjPjRxpOTk3TRRRfRM5/5THrNa15DL3/5y0Wxhv/+7/+mNCA5DgAAAABgBdHqNM1K3k+r8Ao95K1fT5W9e4VnQ2UWKRZAu1ZyOfIWqten9PhoIQsPfotyQ+uD9Ff2YPDme6k45VfLW5rZ5g9t/24aPvqvWj4+ALLmX/7lX+i2224LxAZmdHRUlKF+/vOfL8pTc2UpfpwWCA4AAAAAAB0A52ZLrxrVeLKd8rEFxSVRhtCbnKXK8CB55QpVCglLZ5Yr1DNZpKXRAipVdAG9Oxf8a2J8b+o+uOxlpdg+nmbTD32BSouTdcsLqw4W//euPrqzDV5bRJk84eOQdZ/ADY5m4IoUeroEl7CemvK91sbGxhryZ0ThZAAAAGAF8Pjjjy/3EICBKN8IGYnQ1sbRlnLpSanMz/n/T+xbljvboMXwe8zpFAxXqMjCv0H2l9E1mSVsQJvvXxuYZLM3AjeOlGAD7YXxX9LS7A5aiSxRrikNuKdUvOpVr6KvfOUr9MQTT4jGj7k6xbnnniu2ue++++gP/uAPKC2J3w02lPjqV79K09MIgwMAAAA6BTZ/4rDIuTl/YgeaC1eqkJUr5F3OXMH3JWi7iAQD7ONQntglHlempzIx3+PSh6L84cYN5BV6xR1rbonMB0FH8Uc3uOXS15XDzEp0ql6zlTm/Wkxleh9ReYnalf2P3UHFx+6nyvysGKc3PEZ9J/4p5UcOoN6NT1nu4YEu5GMf+xg997nPpb/6q7+iQw89VDR+zMuuv/764PfDDTfckPoYiSVFdqx82cteRg8++CAND3eHoY/I0crij79rSbcUZdta+eOkVXdSOuEHV5bnyPZ660sVFVpeIlPfJ9db+2yX9tfuvpXn94W2O+jrtXJv//qLH9T66/1KaLt872jsGKKckOEcDxr+Hs66XGYG3/emz6pXad6dXc7RZNOnT37yk/SP//iP9MpXvrJpx1pR5TOraQ3icWEo9D2ufoe3O1y+k8t4SuPIhqlOIL3JGaqsGiTqy1NlVR95M46u99VqGD17F6g02ktUQIh0p1LOVYRhpLdzPMXO5ZqpaIbpFPwZLU35EQWl2V1BFEK5OE19WVz/TYLHLb9X9N+P6m/Bto6K0ihRXrRs+1yeajqdyKpVq+gTn/gEXXPNNfS73/1OLDviiCPEcsmJJ57Y0DFSxZs87WlPo8cee6yhAwMAAACgdZx55pn0wx/+kK6++mp617veRaeccgp973vfw1ug4Ic5+z/UeeLBP+D5f/F8cbajfsRHiVxNF/z57nSaO9TlUrWVydtXnxsPQIhcc2ptLmy/nxZ3/pIWtv2UiuO/ofLijIhO6sTPP4+5/sZW597wA9nw85//nMpaVBELDMcff7xoqtjAcKDB0tJSawWHN77xjfTOd74T+aAAAABAh3H++efTww8/TOeccw698IUvpL/8y79csTcRuiV6Sk4qWjYh4vKGsqXNmed0irG1RPke8nqrRpiL8dEPSKvoHJ716YqDj0NPUJ2CCSpUpPWDaLJ/Q35gnRAm2YfBNYKzk1C/Q9TI0zt++sGmHXOJ8k1pwM5JJ51Ee/b4JWFdOOOMM2jLFr9SSxpSubScd9554v+nPOUp9Gd/9md01llniYEfd9xx1CuNWwAAAADQtnCJK3ag/rd/+ze65ZZbxM0E9njQ72yA9kJEXSzKVI5ikL7R7LuW7OPgDY9Uc+ILNbEhYWnMymg/5eaWqLJuNXnj4TS9Oqp34PL75olW91OpH2kVnUYlTzSws5QsnUJGx+i+Dg36NzSLXL4/JGKyCSQLEhwdhSoVoB2pVCoi0nFwcNBp+0YqVKQWHPhOyM9+9jN64IEHxP8cnvn73/+eenp66KijjhJhGgAAAABoH9j86Uc/+pFov/71rymXy9FTn/pUet3rXkcnnHACfeELXxBlsb785S/TqaeeutzDXfHw3UWevCxH+Uvp48DGkUxubEND/Xkzc5Qrl6k8WpuYRVGZniRvcBV5eyeosmasoWOD9qJSFRO8frdrIfBvyMBEUhhGso4xtZvKs+53dxuBK09Is1iGK1WU5ltz7E6hGVUlUKUimmc961ki0jFJhMPAwAC1VHCQDpYc3SDhqhUsQEBsAAAAANoPNoo8/fTTRUrF05/+dOHhoP6AeO1rX0tXXXWVMJP85S9/uaxjXQlI8cAXElhYcDRTbAKBqFE1jmwIdWLIj+OqUFTz8L3+QVEak9MqZIqFTn6+IqIcOK3iO69tbJigNekU/eMRUQr91VQa1zQIvlbY40M+Fp+jxkLnpWFkaf++kOliK5Gf/aWZrZQfWE/lYnMjMgC48847W3oSGit8q8AVK/7wD/9QNAAAAAC0F48//njsNlx3m8MsQYMEaQ4Fcd+uUpyl/ECvmNTw3U5pRLmc+dlqdY2G4NQKrijBk0bHEPhKwZ9gVjavpdx4gSoT9tQKEeUwOkL58TkqrXML/wXthbd9t/g/LkpB+DekiWRg/wbVx6FBVD8U+VnlyIRWsrjlfup70plEc9PC66Rw9NNp6dGfUTdSbkKVijKqVLQVmQkOoAll11rIcjjvRh2z3R101fE1eu701xr+QWgukWnaL9guQqF3Hbf6h1Xdh13aVYq//WFtu9/X7pR6hXDYVe+aY4LHrODXxtD43QTXayVcYrSQ6sd54s++a+lE0Bl08He8Kxs2bKA77rhjuYfRUagpD0yuf8y/Vrr588+58TbDTb7jbJo0lsqU2ztH5TWDQnDwitHiRGVuRqRV2ECUQ/vyzH8vU47sfhuei9+bUg6zEf+GrJBpEOVF9mUYET4NywX/hpG/x/RSmO3++9nEUiVPeTb7yLRPlMVsJ5pr5QoAAACAjsHzPPqjP/qj5R5G2yFKZJYWg9KYWYr0jU4QpJgqx5S1mM8+DuWJXYFxpPuOZb/xxK9UIm9uv31bFim4asXQsB8iXzKXX/PKRD1zMZUPQFtQzlWof28lJiqh7JROEURGyHSKBEjDSNW/IStkpQpOiZJtuVG/nzqxjCfoTiA4AAAAAAAokVfyR3t5cSaT8+IVGvRFqBJMJMr+GLMkeM379wXGkalh0SBmIslRDtxKm8fIW7PGMqaqWeDefeSVy5Qf998PiA6dQ+6JqlhVjL5eveKSezRDg/4NpakdQqiT/g1ZI8WH0v7xzD6nWaZTtJsQ4ZtGZl0WE1PcdgLvBgAAAABWLO1wV9KGTDtrtwlCHbJEpsPdaq9UCVpouS46cJTDXLTXBadVgPZLp7BSLhnTKYTYENouQTi89G1w8G+wRTeIyKUUhpGuaRWcjiFEwmpr17QH9XtmOU1sQfcBwQEAAAAAXc/t97+P2hWegKQRPtIIEepkp6GJz/xC4vQKb2aWcvtmo+9cD9ZKCAo00QFRDp3BwD6PclvHQ8u8/qq/U19fWKhKUQ4zMrqhAf8GNR0h3efL7XPcs+rAUCpGu6RkmPjBbz/W1P6zj27wG4jmOc95Dk1MTFjXj4+P0xFHHEFZAMEBAAAAAKCFcKUKpl3vdJoQPg7lElUm7RUl6rBEPOiiQ6WvVzTauN5foIsOvM3EpLEvRDm0D8/4bCkQGwLyeXM6RVVsqItukOhCg1YOM5F/w4Q9PUg3w67royo8+ClWrakus7CzO6tRgPYrjbm4WPtsLiyEReRSqUT/8z//k8mxIDgAAAAAADQJr3fQL4/ZP0b5gdWJ908T2txwCgYbUXKO+9Q2P9d9olZdyJm4/PpyRVSqKBc8WloTc2e3GuUgRQeZklGY9kPvITq0j9jQNxOeWlTmfbNQrxAjrsmoBId0ijTeDaV91Wt5akfDn59WpRssPHoPUd4vKNjzpBMy67fdhM5SE6Ibsi6z2e388pe/pCOPPJIuv/xyqlSyN+aF4AAAAAAA0ACqyWS5OJu6jGratAq9UkWj2Mwyy1N7zGPIsQFkj3/Xmv+PinLgO9rFJcqNm/PfOcdftLG1wWSrrrzmeGvuNIN0sIiU37LbHLXAXgtxqRQ29OgGB/8GW3RDeX5fbHRDnMjAVSrUVIhWpET0bvaFh1xvNka0AHzve9+jZz3rWfTsZz+bPvGJT9DZZ59Nu3Zla6hq+CYHzc6ZbAfa3YCq3ce3HPAdMhX5AzPJ9eV6XtVjFQYOse5fmtpqHEN+5IDQdr3rjjEeRw1PXC6DIvU1RZ3HVKXmUkw46tDe964ji3MEQAsJf/cWGkqr4HKbyY+/KL5zhPlcrvEJjugr4vutPLubckPrtYXVu9Bc6nJwQFSQCJU3ZNFBDZWX4fF8Z5rvZCsTxPxciUqDeaKCP5msbFhD3q69tX1ZdGADyWpIfWVygrzRMSE6lNfVUi84yuE7r015EkAm0Q2SwqO7QtENAf394eclLZpBjW7gVAwZ7m0qh2kTGVL4N4jPVK7PyTjS9bcKG0q20pOhk0thLlV6KFfJdkq6VElgPLqC+cpXvkKXXnopXXnlleL/rVu30nnnnUcnnngiXXPNNZkdBxEOAAAAAAARd/tF9EJVHGj0B33atIr4MG/3KAcXkTqqJGhl2uynEBIdoiaH2qSwkiMqjvaY72hzK/QRzc/Vrx6fDtIqwPLTN52nwmNho8ggnSKqioRBJAjEBonNu8HSb2XaXkGCoxviTVzdbxRGiRCl/budPr8m+v7gWYn30celtnYFppHLx5ve9Cb6+Mc/LsQG5sADD6Tvfve79LKXvYxe/vKXZ3YcCA4AAAAAAAqN/khvdVqFM5aIDJfJlfBxKC9ReWJneAVHOZgQKRZ5uxFgqUS5nft8L4ZJPxoiEB14ElkokDemiTNV0YGjHEKHgpfDsnHGTUuB2BDAKTYm7waRUtGTzLtBi24Q/g2O0Q0ynaI0tSs+MsAQ3SA/F1lN1mUkBUc3cd+ckpGUxe2+oWSSdJB2jrgGy8cFF1xA//Vf/0V//dd/HVqez+fp/e9/P/3f//t/6fzzz8/kWBAcAAAAALCiWc5ydHIC4DIRSBI2HbdNkigHTqtIHeVgEh24WULfy/15Km0eq41Tig6cWsGtWvEgEB3KFdEKk/7EFAaSrRcbeuZzxuoQdVEJ/W6fs7roBptZZMLohvLsntjohshxVYUHFwFCfqfEiQqqD0Ta76Es0inUfVsdDbFUyTelgWg+/elP0/Oe9zzr+nPPPVdskwUQHAAAAACwIrj3oQ81tH8az4WotApZHjMrnNIqGvCdKE1p0Q2Suf3CxyEN3ng4WqHU55dUVEUHIVqYqhyU2UCyNrn0Sn57zkdTDQWkEBskpUKFen9fP5mvSKPIjPwWrFj6skU3uEQILIcfAk/2RfRDoT1MIX/w248t9xBAFwDBAQAAAAAg64kHT+zl5L76P5fGzBIZBp7FxMgW8ZCJlwMLBlGlDIslKowv1IkOlb5e0WjdOn+FJjrIO9lSdOiZMZgLgqZSWPDf14HfToSjG6p40ihSFx3mlbvoUWJV1cejGdENST43UXf9XaMTpMDoYk7JFB+7X1zznJ7CUR+VoiV9qc2juOJYolwTfBwwxW0n8G4AAAAAAGREeb8/qSnPh+/cd2pahSnKoc7HQfanTBxFqUwWC1T0SWP1rrQ36UeOqKIDl1Zc2mi5y1tNq4gSHRDl0DyefnMxLDY8apjkm0we9aoUJrFB9fkwEWU+qVCZ3psousG14otrukHW20kK6/5ARDLJSjUAdAIoi1kNQax4zc316cQyNaA1NPvaSFM2Tf1jzH/cJJXFsEt4abb2x7xU/ZFtCtlV7+r1rj4yeFyc2mK9i+b+x9q11GevsaxdGvQf6U19Dx3Ldi5Lmc0WlrRcKd+hjV6bgBq6A5j0x7/p/XItlclpFaHvTadQ6774cpliMmIZA49XRl4Yx24vkWksj1mNcvCGR2sL8jnzBJMnihEh9Cw6VEaHQ5UrvDIJ0aFn52zNfJDv+M7N+KJDdRmLDt7wiBAdyutqOfMsOtzxeushQQNig2TwN1WxgQUm6atRFRu8Xl1wytVHNziQRXRDM0rGunxfcHnMJOaQud5h6t986orxb2BKlR5RGjPbPhHt1E4gwgEAAAAAK4YsfBxSp1g0Ma2imeiCsDHKoVqtwiuGc/tdoxz8A5WpsMsQNl7oocq6NcFTb3CVP6k13A1Xq16A5pEr+2kvEqvYwBEMtqiErKIb9OoUhugGmU5h825oyecooQdM4Q/O8B/k8rT0u5847ZNlGUz4N4CsgOAAAAAAABDrTL/Y1LQK7l8eIy7f2uWupnUbGdmQItqprlpFaUm0ykT4DnIgOsg72uzpwK2vzyw68MSzXK6lVlRFh3LBb4sb/f1U0aE2Bq1s4r5qH1XRAakV2Uc35Iv++zrw2Fx8GkXwnmtTDotvgzCZjIpusJB1dIM/FjXCsC+RF0LUdhx5KlqKqMXFnQ8mFi/a2b+BQZWK7geCAwAAAADAMqXFSCO5XO8qh76j71rWRAb3MfCkRzb53NU8sjy9xxjlUDeuQk/9BDQYcylUBlGKDpxOkVNehhQdpAGlN7TKOuGE6NA8scEr+ZENq34z7593jj6R0Q3y/NtSKRgtAibAFK2iig0JohtsuFSmyNosktMpOE3Cvn44s7KYnUqJeprSQPsAwQEAAAAAIOGdwEZ+5Mu0Clkes1HiwsGznJCIKIdSkUr7nggfIybKISQ6iOX50N3skOhQXKLCdl94UEWHSn8PlTevrW3HooMyKRWiA/dZKgXlNtVIB0Q7pBMa6sSGR5XraXzcLjboIoG8JiypFHp0g7Efg/CgRzfIdIok0Q1xfldxwoMLLmVweZvF8YfJ6x2kpKRJo1hu/wawMoDgAAAAAIAVRVofB/njPGkutiTJJKLRtAphiG0TGixpFQ1HOfCEcdpQRlMPp9eiHEL7K31I0aGcq4i2sNrvR4gO3EdfH3mjYU8MtSyjEB3KFSrsq6kWEB3SG0QyI7+V107OTWzQ33uLcahTKoUuNiwsOEU36OUw/eOZPxtqtI95nNlMyl2+Q5Z+/zP/wULry2G20r8BKRXdDwQHAAAAAIAI0kwykvg4NCOtwrRc9YlwwRct/KYiX1NdlMN09bVWBYOQgWQ+R5X+3sjQ+9CkU/bBJpBzJBrDokO5L0dLm2vO/0J0YP+AqodASHTYOyn+h+iQgdjw+0rtyc6wKaPXP1AvDMioFr4OEqRShAj6cqtMoUc3SLEh7rrP2jQy60pbPQceQ+X9k9QzclB6wTOLSlcApAAJLl2aD9XqL5duOV/dQFTJRuc608ofXvWHZn7tYdZ9youPBY9LM9oPkVxtTGpZNfXng1oi0++PGqKZjtNZlC50LZ0XPu5i45/1FGVE2+07oR2+b/DDrXtJUx5TJ7I0pZJWwRP3pOUxzcdbTF4G2aFEpt+3JjYszoSFkdKSqFiRG9sY7oAFg+EIAYUjFLQ71P7xSrW726Uy9WyfCsQFFh3mRytUGiXqn/TEcl4v0jXWrCbau88XHcqlQHTwBoeE6FBZMypEh+Jq/31ByUw3saGcL1OulKPRR6sLdtenKQQlS21igy26wSI2OJXB1K6d8sQu8vI9kZEN/Jl0+futm0XGi3zR68uL05QfqKUC2cj3+9v0HXBS7LZynFmVw1xOlio58ir5zPsE7QPeDQAAAACAFKKmLJHp+sNdCBEpcrPTVs6wpVUkES+lcGKbqJWmdoj/ZZnMIMpB7q+XyVSjHFRxwcScH9bAogKTK1ZosBrFz8JDnTcEiw4agfCASIdUYgMz8kT1Pdq20y42qOKBi9ig98NGoC5Cg0JlVivXOrHNvJ3heo9PnQjvk0WVCv7OSCJW9xx2AlHR9bulUZEU/g2geUBwAAAAAMCKo1EfBxeRoVlpFa5jtD13pebpUKiLcjC9Jik6BJj8HFQGffElNNHk6AR1klUVHVhsCHarig6lfo8WDxoSrTbmfPiOey5Plfn5OtGBG1fCeO51flvpqOaQS4WSaMzY//jvvbcz7NdQB4sDXPpUYkuhiItu0CtcOEY31J4sUWl6u8W3wR7JEy/uLSQ3iewdsa8rcHWKXvF5dxEhFh/5gf//nl87f6Y7wb+B4YoSSxk3VKloLyA4AAAAAACkvANoS99KknqVplqFPh4RaWEYS+Q4eF2CSVhclEOwXZIoh6rowEKDaHIbVXQoFim3fa9oKsW+WrUDITr091Jl4zr/mCw6VD0dRH+a6CDSNfYuiMasZNFBCg2LfUXRJCPb/Pcpt2VnWBiqvk+BsKOKAqrokDCVIqg4wv/bohuqYoMe3VCeDqdySuKjgLJPK9An+jKdIqpKRZo0vSRjRxogWE4gOAAAAAAApKgQEbfMhkyrkOUxgx9l2oTE5Q5lXTRDVXSIi3JIOtFqOMqBJ6Ky6aKD7gMgx1jkdJVSIBaIbrbvJW+hRENbS4HoINvsgVXTSBYdCgXy+vtF88efDya+gehQRRUdVpLwcMZNS6IxqtDAjD3hX3u5rbXIhsr8/iDE3yg2yMelsnMqhREpWsRUpaiLbuBDz+6mymLNNDQL76V6cS99+kHU55BTtHK9Q1Qc/w3lhtYk6KOxcpjLTamSb0oD7QMEBwAAAACAlFEFria1cnLuMvFxTatw9XLI1og4ZZTDgjJGXXSQqRUG0UHeTReiQ7ksmjfuh8uz6JAre6JJWHSo9OWpvHmNaKLfquggkKLD5HTQVNFhJQgPqtDAzA/6783orkHRhsf989Xz+91hsaFajSJSbFDf4yRig4xuMIkN/L4ryOgGKTbERTfYrtmoz28jEQEmAUAXE52MraviTmU2LJB1GyiL2f1AcAAAAADAiiTKxyFt/rPNxyEuykGmVUSFXcdhm0DJCVeSKAfThMsU5VCa3eUc5RASHeQYkogOi0r1gKroIJGiQyVfodmNfhO7sujAk1eZusET26pwEfRlEB26UXh4xmdLoqnsH/Zf86p9tbKWfb8Zp54te4xig4DPZ5TYYIlGcBYbVOT7FOXbEBHdkDSyweTRwt8Ftu+DqOiC8uJU5PeI7fPXd+jT65YVNh9j9KWwkfS48jW22r8BrAxQFlP8YeaQvdZpL1mU1LOh/xjI+q6HTfFN4vQLsnovCk7Xl8t7Vu/GXFtXnq3d4SgcFi7V5PXW7hqVi7PWspjqH8megdHgcX7kgNpG2utZmqzVd3etOR3+gdDbsjKZaXD9XLjeZbF9r6Qpv+luSOU6tuX5DkjzXet6vpr1mvjvEWhveILhl8pbjC5DnO+tK4/JUQ56GkUccWX54sajlrdM+n3gvwa/D/W1qX3xchYd8kMbRJRDfmRTdXmJKhN7KDe2NlQmk0WHSh/7K+SCySSLDt68cu4KvWEPB/U1LS7WjAWLSzSwxb9Tv/+QfiE6lPI1c0kWHfqm80SHjFLvlklfdFAnrnz86iSXRYfK6HBIdFha40/aWHS4/Q3U0ahCgxQZJCPjgyGxQSXwSjBFNZiex4kNNv8GVWwIvCKqYoMSKVGZngx5c8joBhYb0vwNS+aDkK5kLper5eNwukTU8XqG/M+Oih7dIPfnUpvuFTMaK3PdbEqVHvIqPZn3CdoHRDgAAAAAADRIlKipRzlIMczm5ZBFtQpXL4c4XMtk6pEOlWl/4lqeqN0pD/qQkQ76HW0ZhWCKdFBMJUWkA5fU3DshGjOwZZ7yRY9652qT0VKhTHNrqoaIh4wKYYNFBdn046spFoxuKtmJ0Q7P/PeyaMzs6LxokoHpPtEkfY/612llespvutigl63Un7PQkFRssPWppVH445q0+jaI9dXoBvcytYsxUQ7pvBukEJAfWF+3jkWHKAFARjfkxjaGli88do9hnAsoZwk6AggOAAAAAGg6e/fupb/+67+mkZERGhsbo1e/+tU0MxNdNvGss84iz/NC7XWve13bpFXwxITTCmzrTESVyLSlVbjeyRSVKiImW6aonzTRVrrooPahejnUiQ5KmUxddNBTK/zOSrWmCQ9qBIQQHopL1LfFFwtYdGCxQcKiw1J/RURBcBP0FaiyzhxtoooOqvDATQoP7Sw+POvTlaDNrNkvmio0MKrQwAw+VPUYmZ6qLWShoSo2BFElUggwiQ1RRIkNUvjR+1SiGwKxYaomYJm8G+ojNgupI/2yIN9frVDRWyvdKqMeZGlMGzbvBo5uSEvU5z1tqeBGgYdD94N4EwAAAAA0HRYbtm/fTrfddhsVi0W68MIL6bWvfS3ddNNNkfu95jWvofe+973B80F1QrpMqHc6VTGAowp4IhHe1k8/4CiHXET5S45y0MUIjnKwCRr6eELjkGkQ1fGY0i3C+0evN6VW8GsyRTnUpVZMj5M3vE6IDiK9oppaEX6hfnpFkFrB7/HcnPAMkP4BdRPWfF6IDiL9gtdNTJI3NipEh+LmYRoa9yeZs+uqEQ6DfupF71yPEB0GdlajJdaNkTepVTVgY8qZOX/9qvD1JiMeSqO9dLYyP/v239aMK5cDFhdUZlZXz5tGT7E2qR/9VXibyoQSiaMLDcFyQwnMOHNI9b1jUYGfpxAbVGypFDb06zWJ4au/vYwYctveFN3g7+9ehlaPbihObXE6NgDtBgQHAAAAADSVX//613TrrbfSj370Izr11FPFsn/7t3+jP/mTP6F//dd/pQMOUPxcNFhg2LSpPre5nWBRwJQCETWRl0IEp1VUFv3JrYxyKBmMJtPmj5uOmyavO87PQa5TRYfckH93tzKxnXJjm+tEh8DPwVV0UKlOWoXYIMfCosPQEBW2+OeveMhqITxI0UEKD17Jo+Ih/vORLURLm0fE457tU1QphH8ae3O1yIDKYM23KD+pRHUI8SE84W+mCPHsj9cel3PV4+btIoPK2l/IcSqT7717wxvFiQ1qKooUGxRPDitSVOB+FP8MV7FBmJHmeoxiA6dTxEXrmESyLLwd4qKR9OgGiS52sCBZ3PUw9R54grX/NNENcYLiclOu5DIvY8l9gvYBggMAAAAAmsq9994r0iik2MCcffbZlMvl6Ic//CH9xV/8hXXfz33uc/TZz35WiA4vetGL6F3veldklMPCwoJokqmpqfjxPfQhOuPovzWuSzrRN0U56LiYR3K4NRvmJo1yUEUB3Rw4yXOXiZicwNlEh/LsnkB0CF67aiRZTa0IiQ7Vib2Y6MeJDqLDquCgGAlKWHgobxilkSf8dVMHlYTYoDJ1CNFQ1SeRhYf8eE38IfZ8YAFkeo4onwtV2QjGbBAfJCYRIqkQIVM3QvMn5aXOrjFfm8X+2sQ6X/R/7m98kK/L8ITbJDZ4ugCgCgymShLigEW7eGAre2kzkJRj08UGPr8TW8W1llRsENvFiA2mChXR2y8YxQb+nytU5AfWGY/vNNbF9OKivaJG+6VTgJUBBAcAAAAANJUdO3bQhg0bwj9AenpozZo1Yp2Nl7/85XTooYeKCIif//zn9Pa3v50efvhh+vKXv2zd5+qrr6b3vOc91CrkRD8qykE61OtpFXqUg0yryDLKQRVAamJE8ooVcXeA9T510aFcjXLwt62lQNSJDrI/g+gglqvCg5ywilKXJSE6VGb99AiOdOBJcG5r1UfiwHVCeJjeHL4LX+ot0dQBfj9De/uodFDtbnTveHWCODwYjnTgKAj1bn5wpz5H+ZmlWt+rzD+zbUJE0H8+WpDYP2qeOC8M1q6NXDksDPhiQ5WdfoRA3SjyebvYoFejUNMokooN8rltey1FIxAbpqr/ZyA2xEUxmMQEl8+ePtm3Rjhp0Q2F1U8S/4eiGwr9VNqzhXrWHkkLT/wo9tj2MbVvdENQUQJVKroaCA7CuKU39YcxTchVmj/yUapsM81wsviSavcvumaVvVuu153uuOZwPx3dGTp31HHB495c7etk8dG7w/spd+fUHwqFtU+q9dXvh7QGzwtDRvdz/Ud4luUu9b7i7lK2gqhSn277Fxv+vmj0Wm5l2dxGyw5ncb4aOY9eJSYkuc14xzveQf/8z/8cm06RFvZ4kBx33HG0efNmeu5zn0uPPvooPelJte8Plcsuu4wuvfTSUITDwQcf3NIoBxeSlMjMysuhfvv0UQ61/2th6jLaQYaQ2yIdKtMT5A2P+SkQVdHB+Ho00UEc3xTtID0BlGgHFh7UlAAWHirrVtPI7/3nU4d5lC/mRFscKorUhOl1vqgwPN5Pi4MlWjykNmkf2ln7m+TNKNeCloJRO6BH+blSnYBQ7ksW6l0crMkCqhGmyuyof23kSuEJ/OYHFfPRJ7bFHitUGUQVGGxigyoM2MSDuMgIR98GKTaoyKoUaYirUFG/fCGR0KBGN8jPRON/nxa7MrpBpFNknFKRdYoGaAwIDgAAAABIxVve8hZ65StfGbnNEUccIdIhdu0KC5ZLS0uickUSf4bTTz9d/P/II49YBYe+vj7RlgM1JJsn5Ookv7w46xTloFes4LSKpOKHKhwkFRXSpFboyNeqig5iuRrloIoOw8NCWFA9EhKJDrJkppwwy/KZ1W1EdASX0dzpRztUNq6jkd9XaPZAf/Pe2QLND9deDwsPfXO1SXaxf4kmDq0dbninH23B9EwqXgiFnFFYyC3UhIJcsSYglAs1ISIkSlh0zokNtfSgct6cjrDx0dq1lduyvdY/pRAbdKFBjepIKzZEbWvzbVCHUDVWTVoCU8VkGOkWJaFuvxAvOlQrVEh80aG+kkxddEPB/xxwdAOjRjc0Q/wEoNlAcAAAAABAKtavXy9aHGeccQZNTEzQ/fffT6eccopYdscdd1C5XA5EBBceeOAB8T9HOrSapBP9tFEOpmoVYnlKL4dguRbloPsupHkdpigHk+jA62TlCqPoMD3duOjAcKlMnjhrVRF4O3VCzcJDZXSEhh5VNjtkMOSLsDBYez25ck0MWOop0b4Da8LQqsGa+FDsqx1X+kKI5aO1SXaxrzZp7zdXPaS9B9ZEpnLeHNUwtrN23RQWav33/ra+VGTt4Ivh9AlVZNCFAJvYMB8ur5mZ2MD9auv1VApJI2JD0IdFbDCnC4U/93HfA7p3g8k0UlSPSRChWZqvlQNNQiMGmK0CEQ7dDyw8AQAAANBUjjnmGPrjP/5jUeLyvvvuo+9///t08cUX01/91V8FFSq2bt1KRx99tFjPcNrElVdeKUSK3//+9/T1r3+dzj//fHrWs55Fxx9/fFPGmTa02DQBEXcytQkFT8LF/wZ/BoajHFQ4yoHNI13Dpm3jqZ9ALSY2sasfQ69z+lF53n+9LDqI5xPKnffp6h3raX+SrXolGGHRoWoaKn0dhGGkNI1UJ9XioNVoh+KiaOIxL9u7z29VD4mB306Ixow90SdaMP5cJWih15Wv0NS6uaCpTBy0EDQbMxuXgja5eT5oNgZmBoIWOiUPTQRNh8UWIcyo54WFBnm+eJKvTvQ5OkgVG1hoSCM2sNCQRGxgSqUgusEkNqhCXJJrNiqVwfbZqN8uWVSBHt0g+5bfARLT59sW3RA3jrjvhXZMpwArA0Q4AAAAAKDpcLUJFhnYg4GrU7zkJS+hD32o9mO3WCwKQ8i56h3s3t5e+va3v03XXnstzc7OCh8G3ufyyy9v63crbZSDOmFvNMrBNhaTgaR124RRDv7jcKSDGuXAokOuf3XqSAf5PBAkTGaSPImW6RXihOXrzCWF6KBWtGDRgQ0mq7DosHSgHz3AooMqGPQUeb/avvODStWKXJnmRmsRF4WF2vu5b3Ntklkq1Ewlyzlz9IJq+CgrTOiM/speArMyG3GNREU16KlIUmjQTBwT+zW4iA087rmZWLGhEd8GV6JSJ6Im/LwuP7A+cUnM0swuKqz7A1ra8z/Us+FJmUU3dArlSp68zMtiwsOhnYDgAAAAAICmwxUpbrrpJuv6ww47jCqV2h1kFhi++93vtvydiTOPZMwRDWHDRn9ZdSJuMGzUvRyCY2Tk5eBqIKlWl6gJB2b/hzgDSRMNiw5y0quXzGRMKRZSTJDCgyou6OaS8hxM17wRuJ+e39fMjdcs1O5UTx0SFghW7a1FGkyvq02EvXKOlgr+MXIlT7RgCAl9aEPVJRYi7tRPmKNmBGo6SVzqhE4big2NRDeon0k9AsmWSpHUMNIUteAS3RCIDUVzlEua6IZ2NosEKwcIDgAAAAAACYif7Nsn6uoEvNVRDq5RDEmiNKKqVkS9ZmfRoV8aQZbdRQfGIdrBeL40o8ncEzVPhNG5sN/G7GG1PtZsHTZGNOgMzNQiN9isUlJYqE3Ue+bU1A2zOaQY6+6aMBJCe211IoNJaNCFggXt+lZ9MfT9NFHINg6T2CCFhmaIDVHY/BOS+jWoqNENshSuy7E5uqGur5ENtLj1Z7QSKFVyTahSAdeAdgLvBgAAAABAk7wcxPLyQtCWw8uh7s5t3Z3dYmTFjTikOFH7P97PwdnTYX5RtJCfQLkshIhAjKh6OoR8HXiyy3f2tcoVtYPbJ/Li+JohpUjF2LkraDymVQ9NB01l/aODtfarQqgl5okdtbZ7T3D8yt694e2kF0N1ks+GkLJlJjboAoK6n+7XoG/Lfch+EooNkkbFhii/hjijSBtxlSlsY7ZGN5T8lBspNnBkky/mwbuhmVx99dX0tKc9jYaHh2nDhg107rnnihQ/lfn5eXrDG95Aa9eupVWrVokUv507w4amW7ZsoXPOOYcGBwdFP29961tFRSaVO++8k04++WRRTenII4+kG2+8sW481113nYj66+/vF8bK0tsoyVjaCUQ4VP+wez1hAx5XPFLUywbr60Z9MSXJB3Uhk/4arE/v5RofQ6XcPOddr5TtOc/y3GV9LqN+HKqUp8Pu1+pPFO+IWt6hPprF3//Q2F9ptpaXqOct9vSvCh7n5zfV/UA1/SBR7/rZwiib8VlwPX9ufS3TZ93xNaS51qI+p2m+Q6N+dGZ9/lxIfUzDOfdca9WDZce1PJ1rtIBLlANPUngCkiTKQa9aYYtyUFMr6teFUy2iMEU6yGoYcr1MrQi9/ohIB7Hf/GJNYFDH2N9LXnGJaLj6N2N6xo900E0S+bkuOsjnemWL4PUrooM2ga5M1t4Xb3CI+n6jlqVQSmVuVCoWzC/S6p8r50oNnVdTF9T99UFpFThCqOkjcSKDKZVCHYN+PmxigxQXZLqKCbWvCLFBJSvPBlPEjS2VIi66IerzbvNuUD8zoiyu9npM0Q0rDeHhUF5eDwdO3+MJPIsOLBC8853vpOc///n0q1/9ioaqHi+XXHIJ3XLLLfSlL32JRkdHhSfRi1/8YmGCzJRKJSE2cKnne+65h7Zv3y6MjguFAl111VVim8cee0xs87rXvU74Gt1+++100UUXicpLL3jBC8Q2N998M1166aV0/fXXC7GBfYx4HQsgLGK4jKXd8CpqwuQKY2pqSrxJzzjz36gnpeAQoomCQ9ZAcGjN+9ltgoN+163n1OfWnuRrE6XK7x51EhzUcnB1RklVlV/0N1/7YQLBobsEhxAZCA7LQZaCw9LSfvr+PW+kyclJGhkZaXxwIPhbn/Sc2nwcVOIEh1zvsLmiQ9VDQU5KpJeDXiJTfibk5KtUjYhg0SEssLrnmPNYgkgExcvBFp1QX43Cfr2b7hirooPehxQdOLVCIkUHsd2wfz6k6BD0bRAehOggmasJBSFfB2ZBSblg9JQLV9TJtTY5D6IsZEUMiSZ4eKNjZpFB8ZSIxBS9IPvu7TULBDbPBnn8OM8G18oTJtGi2reIYsn31IkNMrrBJDbECc326hLF+muSo420CB6TUWSc4KB+rlhwUKMb+LOtpjKJ16W8HrmOBQdbdMPi+K+r+4UjaJajOkXa71GXPg/73Y2UGw7/vmyU8vQc/f6IV6Ye7+7du8XknoUIrozE/XAJaPYh+su//EuxzUMPPSQqMN1777309Kc/nf7f//t/9Kd/+qe0bds22rhxo9iGRYO3v/3toj82QubHLBT88pe/DI7F1Zq4bPStt94qnrPIwMLHhz/8Yf+1lMvC0+iNb3wjveMd73AaS7uBWykAAAAAACl+kLv+0I8yVEwDRzpwlIMrtjKZpjz2uOiwpBFGNUPK+omfXi4zLr0i6FOmWKjjKihBu4MDftMm/4K+Ab9JZElN6feg+j5EIVMETJER1VKUelpG3XZ799Ta9FTQIpFpIgaxgb0aZKvfrxAWG3jyb2pJohrSig1VjGLD7O7EYoP8X7/+9Gu5fv2ie1pUhNhQZxqpiA02bB4PSX0bUAozXtRQ24KeLmSBJ/XS7Jjh8sxcSenss88OtuEyzocccoiY5DP8/3HHHReIDQxHJkxNTdGDDz4YbKP2IbeRfSwuLopjqdtwVSd+LrdxGUu7AcEBAAAAAKAJ2O5M2gzr6kwiq5N1GfkgK1Y0Mumou7NbHUvUpC6N6KBOuNKIDlJ4UEUHVXgQokOpXGu66MCw6JDLkdffL1oIddKu+z0kjXaIglM5ZFNFClPKgktLIjJwmUvZNA+M2vgsQoPaHwsNqtgQh/ralP7le1me3tOQ2GC6xnSifEnql7tVo4j67NnGYopuUD/LuX7z3XfX6IZuoFzJNaUxHBnAURSysVdD7HjKZXrzm99Mz3jGM+ipT32qWLZjxw4RoTA2Fo5CY3GB18ltVLFBrmfitmFRYv/+/TQ+Pi5SM0zbqH3EjaXdgIcDAAAAAMAyeTnIiYipRKaKyc/B338mcZnMOh+HiFKZtvFHeVOYymWaUvfkNiw65Ic2ULnq65MbWhvydZATVenr4D8eDtIohNAgRYfqxNiTk+r+6uuaXwhEh4pSjlGIDGpEg8XPwZmoCXnafpMIIKqwogoLJrJMnzAJDQqq2BBsrhpESrEhwsgxTmxwEc3SVKewRTfI//MD64JUJ06r0McWF8lUntpNXv8wVeanmx7dsNJ4/PHHQykVbNQYB3s5cMrD3Xff3eTRrRwQ4QAAAAAA0KQ69a5RDrJihS3KIXiqCRNZpFaYnpvuENf3lyzSIaqP0uwu0RgpPNhSLPzH02b/BrleN2Bl4aEqPsiIB9E4OkCt6MCTezVtQU9diIk4iCRB5EIigUFtMhrBZkDLYgs3KQjoY9D3dx1fSrGBhYakYoMJfXtV4LJF87j479hEPNNkX69QYUN+hgsbjvH76q95lPSMHWKtXtMu32lZU6nkm9IYFhvUFic4sPniN77xDfrOd75DBx10ULCcjSA53YG9FlS4MgSvk9volSLk800x2/DYBgYGaN26dZTP543bqH3EjaXdgOAAAAAAANAASbwcZBPPtTKZceipFTLKQYoOLuOIK5UZJTqYc96zEx2YVKJDuUzewqJooX5zubDwwFEQqpcBT6TlRF2Oubc3ZLgoBIn+gaDVoQsHEd4OaZFjCpoimIQEApPIEOXRoKLvn0RoaEBsCLZrQGwwIa/dqNShLFIsZHSDFBt0Xwb9M2RLi5LRDYvbfQ+HnqFNMZ8tRDdkCddQYLHhK1/5Ct1xxx10+OGHh9afcsopotoEV5WQcNUILoN5xhlniOf8/y9+8QvatatW0e22224TYsKxxx4bbKP2IbeRfXCqBB9L3YZTPPi53MZlLO0GUir4wz+4jvKFxt1RXe4GZFfNotiacn0p9s+yRGCrcX4PW4Tzucz1ZNp3+M5A7ZzUhdj+z2PB44U/fErwuGedn/Mm6eut/WEsbvm5sT8OKVTJbzistt3QaO3x8LrwGOYUA6rZWl3y8v5Je6UEQx5xK6o9uFR4SP35sewX2Z9y3bheD1kQzuOO+MyVl8zbRezTqs+w8zlxrf4hc9yLc40MCzTpjqBLxYoo1O+ZqO8SvrPJqRWyRGYAXx8x13aaUplybIGbfjW1IipdwrXMZ1R6RVy/LDrIFAs9vUJsq6VYeNWydWLZwiJV+nprKRVqmkU17ULAooOaciBFh+qEORAdIipQiL5VY0hpQCmJMp/U0iTqPBiiJvu2ShOMoylebRwGkcI14sJiDFlnDtkmYkN4m3B1Cj0lKUps0Cf7cZENpmvfFN2gUpzaQishuoEpc0nMrMtiJuyP0yi46sPXvvY1Gh4eDrwQ2PeBIw/4/1e/+tWiXCUbSbKIwFUjeIIvq0JwGU0WFl7xilfQ+973PtHH5ZdfLvqWkRVcDpOrT7ztbW+jV73qVULc+OIXvygqV0j4GBdccAGdeuqpdNppp4mymLOzs3ThhRcGY4obS7sBwQEAAAAAoEGiPBRME/zg/+okn6Mc+O6oFB3qqIoO0suB75RymUyOcuDccZdx1I8rPBbzusb8HLISHXgMlenxkPDMk9rc2Fo/RUARDUSkg2IgGYgPsnyzFB7kZFsVHtTJPE+o9Qm5ui0bUg7WxI6QNwQf15SOIbddjBG7HXLNA6IiFxiTQaWJJKkdlqgGm9DQDmKDbpiaJqoh1ztifGyKbpCRS7YUDik26NEN+phA8/noRz8q/j/rrLNCyz/96U/TK1/5SvH4mmuuERUjXvKSl4hqF1xd4iMf+UiwLadCcDrG61//ejH5HxoaEsLBe9/73mAbjpxgceGSSy6hD37wgyJt44YbbhB9Sc477zxRRvOKK64QosWJJ54oSmaqRpJxY2k3vArHkKxQZP3XZ539eepBhIMZRDgsK+0e4aBGG4QiHBa80HY9D/w2NsLBKwxYIxxUKotaDjIiHKonDBEO+vXaiREOS8U5uuvbL8u01vlKJ6v68S5RDnETBN3RPvi/atgoJyxScAhFOTDV61v1eWDRQRUcZJSDy2QlbjymShPh/bUJVILUiahlaj8sOoj/R8K5yXq0mxAeTMaJiqGk35GWYqFiM1mUk2t9Qh43kdcEiLrKD1H72jwYTONUx8HrbOkcJkGBj5Ox0NBqsaEmJITFhrhUCl1kcBEdpMgg/1dNIl0Fh56Rg8T1bBMcZHRDkiiLZkc3ZPU9aurzgIe/SLnhxiPNVcrTc7TtqP+Fv6VtAjwcAAAAAABaiK00pauBpC5ESC+HRv0ckppIJimXaZtA2qoNqGaS4v+pHaHSmTLaQVKe2FNvjMjohpJaKU2nig42PwPpB6H6QqgMDoabWqZS3zfKi0EtZ2kap+rRYBIbbGUuo0QNE9W+OZJDjeaIEhtYaNDFBhbF1PSfLMWG0DoH34ZGxQbZpy7IRYkN/vPw9ohuAN0MBAcAAAAAgAzuEMZN8sM54vWTLNU80upUb6lakVZ0MI3HFIaelehgIs4vgoUXKb6owgOLDqrwEBIdVOHBNlGPQjV/dDWBjDNwbARVVJCvSxVWVEwVMGzjSmgMqaeN6CkUUmzQhQZVbDAJS1mKDXq6j/mxm1+Dji42mCIbTMczUmThJlzFhqMbXCOlOt27QVIu55rSQPsADwcAAAAAgIyI81BQ/Rz06ALhx1D1cpDUGUhWJ1heYYgqRV+gMPk5uJpImsajekuoy+x9uHs6BAaVdWUMzcvlaxCvRzkXUnTg0HQpOuTGNoYMJQM0j4fUYoAqOvBEXRUxdEFCj3iIEzz0/ZMKJDpxrzFhqUtdaKC5aapEmEJKTEJDFFlENkRtHxXN4JrCYBIbbCKHKbqhZ/XBobEvTWRvENlJVCq5oIxlln2C9gHvBgAAAADAMtwp1O/wql4MUakVXq+f78yig44a6ZA2yiHacC/5BM9EnCChIyes+rlQ0yzKEztFE8efnqgL9Q+EBz1FIc3kPq7cJFeLUFujxEVa6BENUett2xgiOvT0CRYaRNMiGpZDbHBdHyU2JPM6CX+WTJENtvQkKTboSLGBfVtK83tXXHQDWBkgwoE/5MPrKWf4o90wSkm3SErx22VthOZueJbiEklhYNg2uL5nrcLxXHpp3ifHfTx1DBHnZ+Ch2o+OHX8UNoAcWF2rZzz6s5rRV+l3D1r7q8zWylp6/UNmYzDeTrmLlVvYaNxfPJ+v/eApL87Fl14UA7R87tKWkHQoQxm5v/qeRVwboetBf58t+9WXZFOPlcLRPCr8WC0Vp3z/Veru8C1ZtluyX5OW79O2/A41vBe56h1r0NllMtNEOaiTLTXKwVa1gkWHivJdJqMcZH/h6hDu47FV0UhaucK2LHyezFUrbMt54iojHRg12kE1lZSigzHiQY92sCG/31xSKRhddIgzkGyUuAiFtEaQERENNkNIXWxQjSFdxIYklSjqowcKWqnlqFKYdrHB9vmQnwuZSuHiQ6KmRuno0Q2S4r5HqbS//rw2QseIDRyNkHVEAiIc2gpEOAAAAAAAtBjbBEdGOeiTljoDySpqlIP0c1BDvmsVKJL7SySNdHCtQhEaf77XOokzLVcNB23eDoJ8T+gOfF3Eg83/wGXi7hIRkdT/oVHS+EcYxpYkoiHOq6HZYoO/zC42mMpRNiI25HqHg8+oXtnFNj5/G12cS5dKkTS6AYB2AYIDAAAAAECLDSR1sk6tUEUHfUw8wZKtVaJDs1IsYoUHbaKcWniIwyVloV0EBhehYWG/3xIKDc1IobDhIjboj5MaRkqxIe6zJY+hC4UyUqk0uzP4vOpiw4qObhDvWa4pDbQPHRz7DgAAAADQuUSlVvDEqTS/QPn+NalTK2wpCmqKRf0YmpNeEbXcNk7XFIvguTSVzBeCyXB+pJZqJyfNueG17qkWpugE0wRfj3rIWnRwERX0YzpEVtSlTlRFBjU1JdSlxacheJ5xVINdNEgmNiSJbDD7NQw3mEpxuDWVwgVEN4BOBoIDAAAAAMAyeDlEISfnbCQnRQdb1QpxLEPVCr4TW5rfY/RziBIdTONIIzrI45r6sx3LhksVi+D8DK0P1rsKD8FxaCy62kRaQaDRqhM6JkGjQYEhidAgtmtRVENUadYkYkN5ccpxjGGTSFVskNENtlKYKnIbKTbk+v1oiUoxfM7johu63iiyzBVfMhbosu4PNAQEBwAAAACAJtGIgWSdmWQ1yiEQHaqTL2mcFyc66GOyiQ5yWVSkg1if64sUHUyvQS5LO/mMinZghKlkdSKcRHjwV/ZQZc7vxxusRk6YKlC4ihAqaUtxRhEnMJjGHiU0TEWbQTYiNLRSbDBVpUgrNqhEiQ362HSxIdiuKjbkh9ZTcfw3TmMCoNNBggsAAAAAQAqyupOoeyroE2qOctD9HPQqOnKybzKRNE34oiZW6rgivRnqJnrFRJEL9RUHzEaROlHbSL8Ibrq3gO49IDF5FLDwIMWHxCUxW4FJbHAZl+LNoAoNuthgPE/V8ykNIZNENWQpNqieIFFiQzDuBsQGl1KYMpVCf42q2CCjGyRSbIg7h10f3aBWqci6gbYBEQ78JTC6zpiblRi11FuDZTD1EnEetY668ngupCnL2O44vk8tO0e20oSJxtC8EDM1THPjA+HP05YzayGrM8+u3RnatPkZtaE9ppl8TZhDDCvz4R9KXr9SgnO09uPaGx4N77dY+xGSU35sVeZnre95pTjvdj3Y3lut7KG1fKmy3Cv0uZWn1N9LZV3oM1y3Xc68LuruW96yruQYIqyHElvuynkR24W+D0NlNR3Lbzb7891wOWH/vcgths0CQXeQJLXCFvFQl1pRnK0r6W2KdGBMqRXquOT/ptSKuBKeMsUiTaSDbbktikHfRu3DRGDkp4gOMupBnUzLqAc9lUCU1dREhyDyQUWf3KeJgEiDg9hhEk1MkQxi+fSu+mWKUWmSaIZGTSGjxIbgcYzYEBdZZB5zX51hJPdrMorUfRtskTu2VIri1JaGU6+6AqRUdD1dOEsEAAAAAGgfL4c0OKVWxIgO5fl94o4sT4psokNtv3jRQU+tSCM6yOPHLXcRHdRtZT/mUoi1Za7pFvrEPDey1jqBrxMhXKIe4kSJlJETtqgMm8jgIjRUIswQm1N9whwxk0ZsME3kXSb3anUKKTbo13SU2KCamepig55K4SJ+dG10A1gRtFW8yXXXXUeHHXYY9ff30+mnn0733Xdf5PbXXnstHXXUUTQwMEAHH3wwXXLJJTSfRWkjAAAAAIAMSTJ5SJJawaKDjVz/6lAYuHqXVk6MTOHjUZM06wTPMb3CtK9puavYYAvdN5XiDJblC0EpzXBJTT+NwJhOYEg70NMvItMw4tIfUqRp6Me1RTPo42aBQW3itc/uDjXR/+JsIrEhKnVCrosTI5yunyaJDTWTyHqxQX6GbGVf/f17Q6lMLG7pIoUUG2R0w4qObJBUqqaRWTbuE7QNbRPhcPPNN9Oll15K119/vRAbWEx4wQteQA8//DBt2LChbvubbrqJ3vGOd9CnPvUpOvPMM+k3v/kNvfKVryTP8+gDH/jAsrwGAAAAAKw8mhHloEYV+M8XxV1TNbXCFukgoxyk6GCLdKhtXwvpbkWkg0uKha0qhQuxk1qZc6/k40vRQVYA4XHrd/5zwxsiIwWC7ag1uIzFFsEgkeKCiklkSJLiErXcTWzKTmxwLYOpigxRJpH6dWkTG8TjkQMiUymSlOWMA9ENoJ1pmwgHFgle85rX0IUXXkjHHnusEB4GBweFoGDinnvuoWc84xn08pe/XERFPP/5z6eXvexlsVERAAAAAABZ4/KDP+kkQjWSVHPojSaSdccqxEY6mEwk48ZpinSIDnUvRppJ2u4YqxO5RkP0bZju4KsRDzp6ZIBtIi+jCqKiItLg0q/rGNVIhqjzYXuPJLb3R18e1Y+LqGQudxmuRpGF2CAfJxEbJCaxQUU931mKDR1P2SMq5zJurXS/Ax0R4bC4uEj3338/XXbZZcGyXC5HZ599Nt17773GfTiq4bOf/awQGE477TT63e9+R9/85jfpFa94hfU4CwsLokmmptxcawEAAAAAsiCNQZwp2iCJiaQa6SCRkQ4mE0n1uHK5aTyqKBDynFAiHWrrayU8k0Q7yP2TeDokxXRH32Q5y6UMdUxREHXbZCg6xB2/kSiGYJ3DeY4SgVyjGlyOZRKz6pbHCAxJxQbdf8EkNqj+DTpRYkPUfmlBdANod9pCcBgfH6dSqUQbN4bNevj5Qw89ZNyHIxt4v2c+85lUqVRoaWmJXve619E73/lO63Guvvpqes973pP5+AEAAAAAXFMrGhEdOMpBToik6OBiIimRqRVxooPq52BLr6htX59iISeHJuEhiegg17UaU9qFPnF3ESBsIkSzxAWbwJCF0NBoyVIXTBUpbEKD/9we1WB67io28DFzhWFnsUFe01Fiw9LkE9axrNjoBlF5Ki9a1n2C9qFtUiqScuedd9JVV11FH/nIR+gnP/kJffnLX6ZbbrmFrrzySus+HEExOTkZtMcff7ylYwYAAABAd+N6tzHNxEJOnkzlCaNMJIXoUF4MJs9qmLiaXqGncKjjtE3cXM0kXdIs4sL3l4Mo40TdaNE20TelOaRpjYwhKl0ibdqEvj4qSsU4pjrxwL8m1PE0S2wwlb9UIxuSiA3q5ylObKhViunL5DsB0Q2gE2iLCId169ZRPp+nnTvDjsD8fNOmTcZ93vWud4n0iYsuukg8P+6442h2dpZe+9rX0v/5P/9HpGTo9PX1iWYqZ+T1rYquQ58Gvaa8DVtN+ohd1Jr0WeDlV44SmMW5y/R8ZdFX1tduo9f45HTo6cEP1MJ+tx43GTx+/ITaD6PBw0ZD+6z7be2OkDdeS3+qzM5a38/Q+6KVHPMGB2pPSsO15YobeEWrcuOVlWtFvW7U5WLdEjmRV75yc8pY1XGry8XdtV7z+6yXVMvX1lUKPfZrQ9mOcilyHMsV+7pS2XxtqMs1PMfvSetII/pO9X1c13+G37UOY/B6Is4v6BhaEekgiTORDFIZcr2B6MATT1OkQ653mMqL09ZxmiId1DGpkQ7+fuEUC39bt2gHdf9mok5io45nEh3U6AdJeX/tb5yN3ED4750NvS9+/1xIY/qYdTRD0jQKVWyIT6FwT6cwj7smokWJDVGVKYJtVY8HLbooKrIhi+oU3SI2eJW8aFn3CdqHthAcent76ZRTTqHbb7+dzj33XLGsXC6L5xdffLFxn7m5uTpRgUULhlMsAAAAAADanbgIAhum1AqT6KCmVESJDhIX0UEdd5yvQ/BYER387dNXsmgWUUKHcdJpWFQhwx3wqmlnElHCOkYHD4BGokTizrfL+5EmqkFfZ46Mqa0vLyo3IlKIDbbKFFFigx6VE4oUGtoYvLde72BwTatig4xuMI1xJaZSBEijx6z7BG1DWwgODJfEvOCCC+jUU08VJpBcFpMjFrhqBXP++efTgQceKHwYmBe96EWissVJJ50kymg+8sgjIuqBl0vhAQAAAACgE0plJol2iPJzUEWHWt/RokOJt69OtMrFaaPoYApXd/F1CIkLiq9DaHmMoaS/rjnCg8ng0HXCbtrONE7VrDNrGk1BcT2vrRIaTJVObGJDkhQKexrFcGZigzhG76C//eJcaPusxYZuiW4AK4O2ERzOO+882r17N11xxRW0Y8cOOvHEE+nWW28NjCS3bNkSimi4/PLLyfM88f/WrVtp/fr1Qmz4x3/8x2V8FQAAAAAA0aKDXi6wGaIDo6ZWRIkOspyfFB500SHKv8ElxUJ9nbrwEN5m+YSHuPfGFX2frMebhcdFkjFlLTTY0ieaLTao6GKDFBpcxQYVk9hA1e31krUQG8zANLL7aRvBgeH0CVsKBZtEqvT09NC73/1u0QAAAAAAOgF9AqeXezSJDnFpDKroIDH5OegVK8SyqujAsPBgEh2i0j6iTCZ10cFfHl0+09+mPdIsdLEkDe1kgtlKoSFtVEOrIxtUscEsdNWP2XRudLFhaeqJQMBYmt0ROy5XENkAOhEkuAAAAAAANAl9guAyAdUnIi4Tflm5gqMcoipXyEkVRzkEy9Tw8Gq0Qy28fDiYnLlMkEz59Hqagp6rb76rba5kYeqv2ahVGOIqNrQLacbsum3c+a9/v2vvZd067f03XScsNDQqNvD/8nGayAazKNYbKTaIx1WxwSbMJaFbxQavlG9KA+0DBAcAAAAAgGUUHdwmg9GiA7fS/j1OokPQp4PooAoP6jiSiA4m4cElb7/dhId2FSEaHUuSfZIIDf4ys9AgllmEBrlMFxoaERtkdEPw2CA2xBmGquvqSmFGiA0mE9akdKvYAFYGbZVSsWz09XLNzHT7qiXmWoiXpixcM0rGmci4ZGfWeHo5wVbRqJnpcpW+zOC6yY3X/tge+ItaObC9h9d+YEytnQjts//k2tfT6Hhtn+HtY+G+J5VSlvOOPzzVEpL9yo+GYf0HRO31eUWl9KVSStP4POk1oL63+jahdcq4o8pdqmUx645V+/ao5FOUxVT28Up6RSDzWClnv04qab5fU7yG+rGmKPXpeu1b9vGc9mm8VBroDOJSA0ypFXHlKRkWHfIDa62VK1SRQfVzUNMreGxqBQtTikXcZE8fo/o8rnymaYLXLmkWUWQxaU/TXxKS9JlEzInyPHBNn1Cf26pQRC2LExskcWKDWrklSmwIqo80WWzodnLlnGiZgioVbUUbzmAAAAAAALqLuDuUxrxwg6jgUkYzKtJBigsukQ7qZMwU6RAX7aDflQ4/Dkc7BI8jQ+zjox3ayTMhSURCMyIm0vSZ5Dzq2+rvT9L0Cb/ZrxkX1OvSJDbk+9dmKjaYRLCsxQZEN4BOB4IDAAAAAECH+Dmoy+TkzDQpSyU68GSxVBQGlNKE0iY6JE2xqKVV6GHx5jQLm7dDnPCg97lSaFSwSFoOVBcZooQGsSyh0BBXHcXdr6FebNDPT5TYoBIlNugpS1FiQxL/hpUgNnjlXFCpIruGKW47gZQKAAAAAIBlQk0H0CtWSNKkMEjRgYUBece1Z2hTkF4hymIqkQ1ChJATqOrkkUUHvwJGfXqFv949xUIdo6nqRlSaRf35iK9mIfvU9+l0snwdaVMm/Of2SJMkqRP+4/B7rKZTpBEbJDaxobZP7XGU2KBeY1FiA3+uliZraRVR4wTKOankMhcIuE/QPuDdAAAAAABo49SKqMoVKnGT/iDU2xDpwOJDIECoE6yISAd/uXukgynawSXNwhaO77f6O+ztHvXgmlLRLGPKLKMZ1G1Eq75PcRENpvdfNYmsH4c5kiet2KBWpbBVUFG3Cx/PHNkgxQbbuU0qNqyE6AawMoDgAAAAAADQYVUrXEUHjkBQw7ujRAfRZ4zoEISVW0SHtMJDeF10NQv9znnchNjUd9biQztVrZCEhZnsvBmM20S8L+oym3iQNqohSmxQU4NsYoNpfPp2teOF31OT2FCa32MdaxJWktiQK+eb0kD7AMEBAAAAAKADRIcoP4dwX/XeDnGigx7tYBId9GgHXXRIKzyoYw6vM0c7+GO2ixRJxYdm0cwoBZ00woKtjzhvBlvkibqdvd96n4YorwbTOl1osIkN4v/CcKzYwOkUcWKD7M9FbEhqdLmSxQawMoCHgyy9ppdfiyo16VgKs65sXcL9MymDaCvDFvH6QqXbos5DlqU026HsZDvSjqUwXcanXxvKc7VE5tqFweDx0IaNoV0mN84Fj6dXzwSP968Kf20NzPQHjwf3DgSPe+bCY/AWlHKtaonLKJTykpV+xx+JSt+hUpr658n2+Yl6z9UynfpncyH559G5KKYypopaclP/LlPPV0Et05mzlrWMwlbyMio107ZPVFlMrxyxnZJXGl6Xt5fVzJnf56gSoN4yl1puBf/4j/9It9xyCz3wwAPU29tLExPhUrgmKpUKvfvd76ZPfOITYvtnPOMZ9NGPfpSe/OQnUzfAE4szjv5bo58DIz0Nap4Fi5HlMuO8HVh0kKIAiw6qp4PYVvF1CMplsuigTDp1XwfV00H1dZDjkmOIQvV00P0dVE+HKH8HLqdpOof+Mn/8UWU1W1lCMuvymUmx9eviyxBaHuHRoC+LMoSMWqdeQ3FpRSaxIbzeLDbo2MQGrzBUd55cxAaYRLqYRmbs4QDTyLaie3/ZAAAAAKBtWFxcpJe+9KX0+te/3nmf973vffShD32Irr/+evrhD39IQ0ND9IIXvIDm5+dpJUY61Mwl3SIdan26RTqIbW2RDhG+Djyxs0U7yLE1q5pFnMdDuA83vwf9mFmnYdhSHVp1rIb6i/FoUI/r/2+vPhG1Tr9m9Me25yaxgYUmKTbFiQ26GKSLDcF2ud7gcwOxAYBoIDgAAAAAoOm85z3voUsuuYSOO+44p+05uuHaa6+lyy+/nP78z/+cjj/+ePrMZz5D27Zto69+9au0UstlRokOUTQiOthSLPwqEbU0AZPwEB579mkWNo+H+Amxu/DQKTRLtHAXdBad0yeiol6irhO70FBLndDFBnV8LmKDvC5sYgNTmt1Z99qSvA6dlZxG4ZXylMu4cZ+gfYDgAAAAAIC247HHHqMdO3bQ2WefHSwbHR2l008/ne69917rfgsLCzQ1NRVqnUCSCYdNdIiLJNBFByk82ESHKF8HnrypJpKhPHlFeDBFO2RZzcJVfHCJekgaAbFcZBER4fKa9aiSqGgGPX1CT5mIEpFUbNdHtG/DcM2voZpaY0uhsb1W3ech+IxFiA38uVHTiVxeh4mVLDaAlQE8HAAAAADQdrDYwGzcGPZ34edynYmrr75aRFN0sqeDS45/zd8hPLGRue5Rng7+/mFfB5OnQ7SvQ22CKf0cVK8JnuSrz+WkkI8XNb60/g617RWRRF2uGRuatgkftxh7lzy+j/YowylxEVJCUQCWybS+XRqfBldsaRWMKmTpfg02scF8jHpTSc9QDlMXG/RoIds4o4DYwLZHOdGypAIPh7YCEQ4AAAAASMU73vEO8jwvsj300EMtPbuXXXYZTU5OBu3xxx+nTiJNpEN4WX2ZQNcUCzXSIc7XgSdk6qTMFu2gP09TycJ2l7z23O3Oe1x6QNr0i2Z6MGRBXASDPmYWGkxigz1KJOq9iY5oaFRsqE/jGYoVG0yiVOj6Va5r9dp3ERuSALEhbBqZdQPtAyIcAAAAAJCKt7zlLfTKV74ycpsjjjgiVd+bNm0S/+/cuZM2b94cLOfnJ554onW/vr4+0TqZex78VzrzKX/ntK1avaG2zI8GiIskMFWwkJMpWwULdQImjtU7RBUlH16tXiH6rU5c1XHq0Q5yLC7URzmoERCW6AaHqAc1FD+pWGCrftFqklaZCG1jiQCxV7RYSPQ8CVHVKKKiGvTqEnFig2m/RiIbTGMHAEBw8OHyclHlH3Uct7UVfosqjRYqieZadjKqjF7O0l9EGVC1nKfn+tpLSsnBVpbPVF9f1BjauWRm1Lij1mXJ/9/enUC5Td2LH7/2jGcmC9nIDiQEKBD2pQca/pQtQAI8SoBHgfLKcmgoBVqWFEigJIQtbdlaeOwPSHpaoKUHKI+dJlAeEHhseSwlKWELoQTKkoQsM+PY+p/fDbavZEmWZdmWPd/POU7GlnSla8mW9PO991fL9ydtNFNdUUh32WuNPf7Z/lUhZWbnYOPvDewXC+tShfdo1ZDCsZbqtJfXajxPGetKGukkE+ms9/EeNFWuLTVka/npM53rMd4v22fJcWxY5vNsxvsYMqcFlTSOj5QxYn6b49dV8ybPTCPabr8YzxppMrOpRKB0l2bqStvrjjSWfukvg5QXmJkK07nPAqYWtqUgblBDhgzRj2oYM2aMDjrMnTs3H2CQ8RgkW0U5mS4aVZRBh1I3hmbQQTi7WOjXeg0sCjbk15VLN/hN4CF3A+cMPDi3M0wKTWc93JaLMvhgX2flrRbKDU6UM55E2eM4lBlkCNJNopJAg1+wIUygwXz/SgUa9LrcyjFah0QZbKB1g7EfMkk90GOUrAwtHOKEvQEAAKpuyZIlasGCBfr/TCaj/5bHqlWFwOPWW2+t7r//fv23dMc466yz1GWXXaYefPBB9cYbb6jjjz9ejRw5Uk2aNKlH7DEJOkTRvaLUa86bRLOLhW1ASbnxyj3ctsFxw2be0OVuEs2sFpWk0DS33SsbgrPLRamsCq6DTTqyXrhvQ/DuFF6DVJY7eGVZ6wxQn2z3qvyjeF2F99P53jrnC8s/1WW4YENhu8IHG8S6VZ/ms1oQbADKR5cKAABQddOnT1dz5szJP9955531/0899ZTaZ5999N+LFi3S4y7knHfeeWr16tXqlFNOUcuXL1d77rmneuyxx1RHR0eP2WNRtHRwawHgfM1rMEmRb+3wzc2oThcoN2POX+u/SZO5vrzukq0d7PMVWjt4baMft/q5ve420KS5HW7T8vO4BB3KbQ0RZEDQIOV4zu8TGHHjFmAoZ8DHSls0+I3VUDRugy3dZen30S31ZVEXCrdyjECPBBtK71NaNlSCQSObHwEHAABQdbNnz9YPP5bl6CKTSKhLLrlEP3qycoMO7jfRxeM5uL3mNq6DlJnp/FK/1tIxSN+k5oMO+sWU/f9vAg/27hPrx3bQf6c2sI3tkNtWc2wHcxtz2xWUM3jifL3coI25nUXr8sp+UcVuGX7rrSSwEOT1UtPK4RdcCBJs8OueUk6wwda6wSPYkOn8wrZdXseYH7pRoKci4AAAANBEQQevm2evoINfawe5+W/ptWF+mgQeckEHz9YOOm2mvbXD+hSe6+eRrgJug0pWK/BgLl9qWpCxHyoJRFSbV2ChXkEGvxY2zr/dnuvjKz/NGEPII9jgFWhYX1ZxFpV8sMHRfcUt2OBEsCHCFg4Rj7lAWsx4YQwHAACAJhvTIac4RaX7L7J+Yztk1n6hHzm51g62G1zn2A5yQ/jNTWHxNqTyN4zOPvnmGA/O8R1y2xkmE0CpsQfcx4IIMiZDfVJilrte7zqW/74E4bafzNe8ppuBBrdgg3nsBA02mMcUwQag9mjhAAAA0MNaOpjsqSW7PKdJ0CHX2iHX0kEEae2Qu3k0f4GXG0e5Scx1P7BlivjmBlHmdxuToVSLB6+WHM5pYbpdBGrhUIOgQ1StF8LM58dvsNJSA5maQYbirhTugQa/wTW9ul+4tWxwtg7xatlQDrpRlJbIJvUjSlGXh8oQcDBFnaLMo7yEX2rH4NmPKuNMi+mlnHShgcrLVK+8KFI7VlpGmPoFPe6qmSLTTL1YyzSZZt0d25DsLFz09F5RuBjp1beXbb51GxS+xtK9C/2/s0l7X3DLqFK6kGVTtRp1Tabtn4tkp/Gep82/19k/00ZqzYr5HQ/GMWBLg6m3ybi4TZtNo51pMc0UnPZ62LR4nB6Srd7HSao9UPrMFiN9ZtInfablkT4zKDP1ZVG6TCOtpW8qTb/0l/l5fNKp2rbHb99mq/Odi6YURdDBPq3dN+hgju2Qa+kggQcz6CCCjO3gDDyYfepzWTfNwINzfmfgIEggwfkrfdDgg7P+USh30EivrAhRbme1ggylpheyUNiDDEEDDX7BBr/uF6WCDbluPpUi2BBMMpPQjyhZEZeHyhBwAAAAaDD1DDqIXGsHczDJQK0d9MLuLR4KGS1KBx5y25nbLvO5+Zrb6+W3YvAaANK7rGRbP49l3N//crepXoEFU5CuLbUKNHiNvWGWYQv2eAQbcmlgK0WwASgg4AAAANCDgg7CPRODPeiwfj7vLBahWjsYf1vSncJleyoNPJivBX29nPmCyHavLCsQEefAQk45Y2eUG2gws3r4BRr8gg1uY4UU/jbGcHAJNjgDDZV0iSHYUB7SYjY/Ag4AAAA9KOjgF3hwtm4olTqzrNYOLqkzc9vg3B65qXTeQOYCD+b85jJeQZIogw+V8gpExCWoEHWQoZxAg9m9ppJWDeufFw9Y6pzHLdhQ6dgbBBuAYgQcAAAAemDQwSvw4NZVwW8wycCtHfSMad+uHoXtKcxXFHgoCpLYgw9+gYQox2/w62rRyKIIMpQTaKgk2OA3JoZXsME8fnLBhigG+STYEE4ym9CPKFkRl4fKEHAAAADowUEHr5v6csd1CNzawWf9uRu/QreKQqrDfLpMZU+VaQ7y53cDGjT4ULxt7vM1emAhTJAhSIAhp5qBBr99XTRIpEuwIbP2K9cywyLYAHgj4AAAANAEKg06uLU6MFs7+AUd1i/Trh9mcMEt8OAVdDDXm+MXeMgHJ4yBBp0ZBsrNCuFVt6jGdqi3cgIGYctze8+dgYZygw1BgwJewQZz+VywgbSX8ZDIJFUyE20ay2zE5aEyBBwAAACaRNRBB2drB69xEnKvmTegZnDBDDzYM1O4/Aru8ot0kMCDjPHg1XTeq9xKulp4qXcgIorWCmGWCxJoKBVkcKo00ODVqoFgA1A7BBxMSSMa5pcv3eTMcR9ENcsOq6XF/X3wmy/I635lF+Wur7C+tXy/Kt23IbbbqmX9Aq4r4bffo9yGtNGXd81a22ypFYWLmdYO4wKjt6N5cKrQny9r/G2bx3HoW6nCCwm//Ww+T69zr4PbchXsC+d7bymjvsa0RNE2GM8z67yPr6xRD4/XLec8afeLbCvZ6r3tKeNi1FGnZEevwt+pwsWk1W7UNeU4jaUKZVgt9OFEzw06iFq3dnAbcNLZlSIXeHAGH0p1tVi/3u6iG1YJQOQDD0Y2C7dlg08r3Yqh3kGFKFskVLqc+V6arU3CBxmKx/nwWl89gg10n4hWIpPQj6jLRHwQcAAAAGhCtWzt4HZDana1cOtK4dbqoVR3C6+WD278xnbwGwvAbQBNr7q5b2f5A1DWI7DgnCdoyw7n++UWZAjTlaXcQEOQ6WaLGEGwIX6SVvSDRkqZiA8CDgAAAE0qyqCDX2sHP7kb2azRkMAMLJQbeMhthz2NpmyT0RIqkw48tkM5rR8qDSpEmWIz6kBEpUEGt/nCBhqCDA5Zar0EG4B4IOAAAADQxKoRdHBLQRm0G4JbVwvp65+7QZRWC87pXjegboMCrn89XPDBvr3dPtNqE1QIq5zgQZDAi1+QwatMZ1eHqAMNXuVHHWgQdKOoHrpUND8CDgAAAE0u6qCDvbWD/abfb2DJ3Ou5eXKBhZYOc7DJ4sCDW6sDe+DDr1uE/cbUTKvpDD446+ddZrTdJqo7OGTwVgelxmMIUmalgYZyyzfnNYMNZqChktSXBBuAyhBwAAAA6AGiHkyyVGuDUhktcvPIcuaNf+5G1ww8uN2Q+t1E+neZKNywJjI+gYMyBqGMXwuHtrKnyetu6SuDlBkm0OD2HoYNNJjHi7NFQ9hgA4GG2pDxGyIfwyHi8lAZAg4AAAA9SDVaO5jKDTwU3Th+E3xwBh7MQSLLa3rvHShwy3qRX87n5jvpsUq/m9uwN75hxkUIuoxfHUuVUyrAELQ1Q5B1+c1nlkerBiB+CDhEmTLQL+VdpakPzVR2Icv2ZUtd6Z3q0JbOzivdpTNVojnNL+Wm1/sXNLWg3/sQh5SZtdoev2MlYpaq/3tnpn1MdBkjbH/tONaMtIrJDve0imYaTCczxWKi3X6RZXkc4wkzRaYjvaftOI7gePD6bFrJTKBjJVT6TGdaTGM+k9XdaX/usTmJFvspyTKfG6k1E23GBXLKkR7NeJ4w9rkymmvrss10mrbvL7+0wMlgryezrvvZcs6XKUxLBFk/0EBBB7dWCF6BB6/ggTPw4Gz1UGp5r21yWy7HLNc5FkCQG/WgrS6iEioYUeUgQ7mBhnJ47T8z3Wql66JlQ20lMusfUZeJ+CDgAAAA0ANVo4uFk/vYC94ZLpzjM4hS3S1KLR9k28xl8+vySLnpF4jwG/fAb7wIP2EDFaUCC1EGGMrpNuG9rPsgnV77OOpgA4EGoDoIOAAAAPRgtWjt4NfVwo1X4MCvu4Vbl4uwAyaGbangt5xfMMKPcyyJMEGFagYZ3ETRmqHcQENugNEwCDbUD2M4ND8CDgAAAD1cLVo7lEpl6cYre4Fbdwu38R4aoYuD8N3OmAYVogo0FKXfLBFocAYbCDQ0NrpUND8CDgAAAKhpa4dyWhSYzPmSbX2LM0kYN+derR+81hl80MJUWWMXVFMU21yJSuprz3BSvH1ux4S5vkoCDYJWDUBtEHAAAABAXYIOYVo95OZ3GyzSDEA4Wwa4ZbuIflDDcMGI6rS2qF2QoZz3rlSgwa085zoJNjSPRHb9I+oyER8EHAAAAFC3LhZBbrxL3dC6jdng1/rBeQNbrZvzKNRi2/zeD79ASfhUn6lAZZrrNlNehkWrBqD2CDjIl1l7m7IkTZ6ZosyWsi5kmMwsw0xZF3VKxKBlZ0NMc6TItLLmulrcU/I501Z6pcx0psj0SpkZNJWmmQIvbCrNMPvGmQa00mU8tsH2HjvYUhX6pDWNPJVmwOOmUn51tzHfB+f7aKSkNNNnmsdNwkyVWE7qxFSL6zT9vWLbPuM4NFJm2r5vnKk0K/y+8H3vjLKL0md67NtE0LSoHikyixipNS1nmk2Pa1yrs7WyVJrClj7TTJPquAg2ntuOB/Nv5/FgPvdLc2msyvrm2LCs6H/lBBqttUMlLSD8Bou0tX4oylCRrvoNftyCGv6BhHTVx2ooJbcNzkADWSiaSyJrqUTGirxMxAcBBwAAANSstUOpdJTlLOvFb3wGaelQjV/m48i/G0fw7h+VdDcpJ9Dg16KBQAPQmAg4AAAAoGatHZxjL5Q7fkPF6yyR9cEvWFHql/96BCSCDtxY7rZVuj/cAg3ug0MaA0F+k30kqu2gC0X8MYZD8yPgAAAAgLq0djC7QNQi8OA2toPJLRjhNjhlXLJUVBpkqMX2eg8M6R9oqORYINAAxIdP51IAAADAPfCQCz5EQW4sncGHamRuMNdlrjPssmHLifqG3ushN/V+D7+6uQVqzEfQ7fIjgQZnsEEyUMgjzHsqgQaCDY1Fxm+oxqNczzzzjDr00EPVyJEjVSKRUA888IBtumVZavr06WrEiBGqV69eav/991fvvPOObZ4vv/xSHXfccapfv35qwIAB6uSTT1arVtnTt77++uvqu9/9ruro6FCbbLKJ+vWvf120Lffee6/aeuut9Tzbb7+9euSRR8reljgh4AAAAIBYBR7cWj6Yz6MMSsQhaFCJoAGFMHUvJ8AgqUbNhxdJTyoP5zpzgYYwCDSgUqtXr1Y77rijuuGGG1ynS2DguuuuUzfffLN68cUXVZ8+fdSECRNUZ2dnfh4JNrz11lvqySefVA899JAOYpxyyin56StXrlQHHnigGj16tHrllVfUlVdeqS6++GJ166235ud5/vnn1bHHHquDFa+99pqaNGmSfrz55ptlbUucJCwJkfRQstP79++v9pu6TLV29AuepSJoxoOg83llG7CVFXCk+rBZKrw4sw3Yskx4ZKlwjoofNEtFGFXcF1XNUuEnxDZ5HjfVzlLhpV5ZKoIyyzOzEoTNUpE0cjf4jYxchywVnse+o+yiY8g8JtLdrtknipdZ5zGfT8YKZ2aKIJLeWSpUHLJUhPHNsbGuc6Wad8VQtWLFCv0LCaI71/Oe1kY1MlqUN+hg9OMO2Kb7jP9QrRYZlXSLKLsFR4nAgrP+foEFJwkyuK0vbJBB0Jqhcb9Hc2Xue8mX6+/DIiTn0qemDwq9vdLC4f7779c3+kJul6Xlw5QpU9TPf77+O07KHjZsmJo9e7Y65phj1Ntvv6222WYb9dJLL6lvf/vbep7HHntMHXzwwWrp0qV6+ZtuukldeOGFatmyZaqtbf33xdSpU3VrioULF+rnRx99tA5+SMAi5zvf+Y7aaaeddIAhyLbEDWM4iPbU+ofBdsvgvIEwbxp8AhOeNxE+N9q2GyvzYt5xwxXqBjNMak7njaj53LiptMxznd/NoVeQwrmcV5DCOS3KVJrCM+AT7KbNV9Cb5hA310FTH4ZZZ8XBjHoFFYIy0mXa/tYXVh7HYdjAhPk8ZeSDDxqkCPp947W8swxjW4vef1vKzBbX4EOixZlW1vhOsAUFnEEUIxjhvtX+wYigqTSNbbA67Re0ZqDCSnV4vw/GfrIFKb65UHAPUhhlm8eD47ixHw8tRe8h0IjM1g7lBh/sgzN2B0p16VeGs5xKlRpssqKyKxhsspI6lhtgCBNUcQs0eI3ZEARBhuYStgtEqTJzQQ1Te3u7fpTr/fff10EC6bqQI8GS3XffXc2fP1/f5Mv/0o0iF2wQMn8ymdStEA4//HA9z1577ZUPNghpmfCrX/1KffXVV2rgwIF6nnPOOce2fpkn18UjyLbEDQEHAAAANETwoVzVanXgptR2VrItbjftQQXpDuGHIAMalYyRYJoxY4buwlAuucEX0orAJM9z0+T/oUOH2qa3traqQYMG2eYZM2ZMURm5aRJwkP9LrafUtsQNAQcAAADUPaWmVzaIWgYNomwNYJtWog71DChU0pJDunPkWmd41aHSVJe0aGhy0srTr6Vn2DKVUh999JGtS0WY1g2oHAEHAAAAxCqlZjndKOLC78Y/6qBAlCrtKuIMNPh1lyDQgFqSYEMUY04MHz5c///pp5/qzBA58lzGVsjN89lnn9mWW7dunc5ckVte/pdlTLnnpeYxp5falrghSwUAAABimdmikbNH1Ipfqk55SEDB7+EUJgtILr2lV7ChnP1IxokeRsaaqsYjQtINQm70586dm39NxoeQsRnGjRunn8v/y5cv19kncubNm6ey2aweXyE3j2SuSBtjhklGi6222kp3p8jNY64nN09uPUG2JW5o4QAAAICGaPGQE7eWD/UOilTyfoRdNqqWDDl0nUA9rVq1Si1evDj/XAZnXLBggR6DYdSoUeqss85Sl112mfrWt76lb/ovuuginS0il8li7NixauLEiWry5Mk6m4QEFc444ww9iKPMJ37wgx+omTNn6pSX559/vk51+dvf/lZde+21+fWeeeaZau+991ZXX321OuSQQ9Q999yjXn755XzqTMmgUWpb4oaAAwAAABoi8NAoAYioVVq/qN4fvxYMYRFo6OEke5ZfNriwZZZJbur33Xff/PNcpogTTjhBp5s877zzdLrKU045Rbdk2HPPPXXay46OQqarP/zhDzrIMH78eJ2d4sgjj1TXXXedLZvEE088oU4//XS16667qsGDB6vp06frMnP22GMPddddd6lf/OIX6oILLtBBBclQsd122+XnCbItcZKwJJlnD5XL/7rfjC/8879GnRbTL82grbxM+akKnWksvYRNi2kyU+WpOqXFDCNMistapsWMWoVpLZsyLWbQdUWdFtNWdqLwd73SYpr71rm8mRbTnM9Ii1n0/WCmu7SV7ZMW0/g7cFrMoIy0mE5mWkwVh7SY3xwPkjt83swNI8113tNVI388qq+c4EMcgg3VbuFQqo6lxmGo9D1yBhqy3fZUw+UiyNBYqvE9mr8Pm7rM/z4sBH0u/eVwvvdjghYOoiutVCJtvxg0bgaslOOmod3INd9i3DQ42HLKGjcURblm0+4X/QlbkMJxc2H0/fELUlR8s+i8+A5Snt885k1M1lF2OkRgwuQXpPCar5xpJvM9N282ggaTaqnCG/maBgLixvMm3r6fbZ9V4xiyBSKKbjA9AhNmIEI/N79vfL6yvXJYBwyY2urgmM/2fWPcaFvdjgtsIwCRsAUp7IOlWUYgINHiE1QIEpjwC0r4TLPMaeZ6/AITtiCFcXHv3C9GYML2XWYGKZzfHd8cK4mutZ7bDPQklaTUrIdSN/SVBiQqHdixku4RUQQYBEEGuB9c0Y+5ULdrbrgi4AAAAIDYcg4w6QxANEL3irhtU6kAQxQtNggwIAhpoenZSjOkqMtDZQg4AAAAoGE0YgCi2i0USiHAAKBeCDgAAACgaQedjGMAImpWtivyFgvlpLEEQqtCGsvIy0NFCDgAAACgx4z7UI8ARKmAQKUtIKLoAlFOGQQZAARFwAEAAAA9dtDJUjfacWgRUapLRC2CEgQZ0MxpMVE9BBwAAEDVXX755erhhx9WCxYsUG1tbTp3eCknnniimjNnju21CRMm6HzjQNgxH5wqDUhUI+1kuRjkEUBcEXCQk8RXK1SiPWtPH2ikK0s4UtuZqe7Mac70mebzbIeRKi8ZbDckjOCcM5VmMm2k2ewy08/Z08DZUt3Z/k57p6s009kZae408z0Kk3IzaIpFY72Wcxs80mcGTuVoTvNLg+k3n9dyQdNqhonCVpritNEEPVbi8L6Y22D8XTRKskcqTMsrRaYwv2+Ml50peW3Pnak1bdzfVyvbFiitr/k9UpRK0+s7JtPLMZ9x8e73fWOmwrTNFyyVpskzrWYJtuW8/namxexcXfjbSAFadEzbUmuun2Z1V/4rZlx1d3ero446So0bN07dfvvtgZebOHGiuvPOO/PP29vrO/geet4glEFUu5VEFAEFJ1osID5jOER8LccYDrFCwAEAAFTdzJkz9f+zZ88uazkJMAwfPrxKWwVUJwBRbsCgGgEFJwIMAOqBgAMAAIitp59+Wg0dOlQNHDhQ7bfffuqyyy5TG264oef8XV1d+pGzcuXKGm0pmlWlXTJqhYACGhJjODQ9Ag4AACCWpDvFEUccocaMGaPeffdddcEFF6iDDjpIzZ8/X7V4dLuaNWtWvjUFEIeABAAf0p0i8i4VMehui7wQHc4BAACUmjp1qkokEr6PhQsXhn6rjjnmGPW9731Pbb/99mrSpEnqoYceUi+99JJu9eBl2rRpasWKFfnHRx99xK4CAKBOaOEAAABCmTJlis4k4WezzTaL7N2VsgYPHqwWL16sxo8f7znmAwNLAkBjkMHhbQPWR1Qm4oOAAwAACGXIkCH6UStLly5VX3zxhRoxYkTN1gkAAMIj4CBRsFUrlNWdsac2M9MtOvqJJlJtrtMSRipNZzq7ZLuxTMpeXraj8DybKqSzyxrFZR1ZwDK2VRWWT2btu7QlXSgvaWSpa+m0p+uzpdbsLIyUnOhyjJpsS3VnpK4009wFTBkYSXpEM4Jp7DO/SKktYaDf9vilq/RKf+m3rV7LONKuevKbL2j6n6ApOOshTErRSparx3tirsv4O+FXHyPVpJk+M+GTPtNMi+lMn+k5zTGf8zsnv17jO8aZrtdMn5lIZ8tP1+v8PHp839jSavql0rS97pjPL7VmkPKCzmNMsxzfzyrdVZQ+0+pepZrVkiVL1Jdffqn/z2QyasGCBfr1LbbYQvXt21f/vfXWW+sxGA4//HC1atUqPRbDkUceqbNUyBgO5513np5/woQJda4NACASch6PukUCLRxihYADAACouunTp6s5c+bkn++88876/6eeekrts88++u9FixbpcReEDAr5+uuv62WWL1+uRo4cqQ488EB16aWX0mUCAIAGQcABAABU3ezZs/XDj2UVWsn06tVLPf744+wZAGhmZKloemSpAAAAAAAAkaOFAwAAAACg9mRMq8jHcIjxmGU9EAEHAAAAAEDtyeDKQQZoLrdMxAZdKgAAAAAAQHO3cLjhhhvUlVdeqZYtW6Z23HFHdf3116vddtvNc34ZtfrCCy9U9913n061NXr0aPWb3/xGHXzwwWWt10p3KSuRck1RVjL1oVcqTed8xrREW5t9NiOdpj19ZqFsq8O+Pdn2Qqwo01FIZ5dutzchyvQuPM+kvJsXJTKFMlrSvYy/e9vmS60xqtRlpPVbY6TSNNJqFqW2M/4uSl1pNqcKk0rTnM8nPaXlkUrTKRFmXWYTLmd6Q696ZEOk0tTrNaY5UyQazFSKYSSiaJYWNG1npXzeh0DL1Go7y2FsU8LcPuNjVVb6THOakaLXmT4z4ZE+0zKKy6b83m/jOy9j//5KGNUwU2vaUmnqOmXc04h6pdUsSp/pkVZTlxfwO8Y2n/svJn5peO1lrSs9LWi6TgAAmoCcQwOfR8soE/ERm4DDH//4R3XOOeeom2++We2+++46cCB5tiVF1tChQ4vm7+7uVgcccICe9uc//1lttNFG6sMPP1QDBgyoy/YDAAAAAIAYBhyuueYaNXnyZHXSSSfp5xJ4ePjhh9Udd9yhpk6dWjS/vC6tGp5//nmV+qaFwKabblrz7QYAAAAAhCAt+6Ju3UdrwViJxRgO0lrhlVdeUfvvv3/+tWQyqZ/Pnz/fdZkHH3xQjRs3Tp1++ulq2LBharvttlNXXHGFytCEBgAAAACAuotFC4fPP/9cBwokcGCS5wsXLnRd5r333lPz5s1Txx13nHrkkUfU4sWL1WmnnabS6bSaMWOG6zJdXV36kbNy5cqIawIAAAAACNwaIeqsErRwiJVYtHAII5vN6vEbbr31VrXrrruqo48+Wg8gKV0xvMyaNUv1798//9hkk01qus0AAAAAAPQUsQg4DB48WLW0tKhPP/3U9ro8Hz58uOsyI0aMUFtuuaVeLmfs2LE6w4V00XAzbdo0tWLFivzjo48+irgmAAAAAIDgWSrWRfwgS0WcxKJLRVtbm26lMHfuXDVp0qR8CwZ5fsYZZ7gu8//+3/9Td911l55PxnsQ//jHP3QgQspz097erh9OiVS7fpj8DlTbNJ/5Eh6pGa0ux9tupmY0lkmk2jzLShr1aDHSZ7Z12OtutReWW9e7kH5zXUchFZ1Itxfq0d2nkEoum7TPpwa5Vkm1pgvb0Nppfy9TXYW4VuuaQnnJNY7mU0Y6TXvaO0f+P6/UmkHT3Jn89rNvusqWkmV47f8iQdNOOtNbpr1SsNrns6VSDJg20kyraEurGSbtpCgceuGESFcZKp1n0PpFkT4z7HtZbt193gcrHTB9pjnNTJfpSIvplT7TmXLTMg5Xq92c5jh2s63u6TM90moWpdbMWj7pMz3SbDrfL9t3jDHN7zwQJuVmbluc37kAADQzBo1serFo4SAkJeZtt92m5syZo95++231k5/8RK1evTqfteL444/XLRRyZLpkqTjzzDN1oEEyWsigkTKIJAAAAAAAqK9YtHAQMgbDv/71LzV9+nTdLWKnnXZSjz32WH4gySVLluRbMggZf+Hxxx9XZ599ttphhx3URhttpIMP559/fh1rAQAAAAAIRAaMTEY8aGTUg1CiOQIOQrpPeHWhePrpp4tek7SYL7zwQg22DAAAAAAANGzAAQAAAADQM1iZtLKS6cjLRHzEZgwHAAAAAADQPGjhAAAAAACoOZ3GMvIWDozhECe0cAAAAAAAAJGjhYPo3Veptg1sb0yigjzqrsxImzPqZj438s5b6e7C60kjcb3oXOs6LdFiny+RSuX/bjP+TnW02ebr1V6YluldmJbpSNjm6+5dqG+6o7Dd3e2FyGRnb2O7pR5JI999thDjSmbsZbeuK9S9bXWvwrZ22eNirWsKeeqTncb731XYhkTa8R6bz9Npz31rmc/99q25bzxYzmCtY98UynI8d+5rD7Z9bdbJaz26bI8Yo/O48VrG2H9OllfZunyPaX7LmFLlf1UVjpJvZArHoaeszzzm8j7bnfAro0K+73G5730pZj3Mv41tSKTtn5GEua5k4SiyWuyfdXNaNuUzn7k5xnx2ztfd65vIpuzPM4UjJGFUL5HOes6nzPoa70nR903G471znktcpiVaio5cAACaV3bd+kfUZSI2CDgAAAAAAGrOyqYjH+RRykR80KUCAAAAAABEjhYOAAAAAIDak9YNiYhbJJAWM1Zo4QAAAAAAACJHCwcAAAAAQH3SYkbcwoG0mPFCCwcAAAAAABA5WjhIGrL+/VSivZ93KjMnW+pEn5RnXukWw6bWDLC85SzLWG/CSJ2Y6Oqyz2ekums10me2ONIRthnpNC0zlWZHe/7vdG97Wrd0u3sqzXUp+7au7ejM/7263xr79pmbaqRmbEkXtq81XUjtmOospNXUz7sK01o7C2n0Wtc4UuCZKfGCptkMeDwETrkZ8Hgo2tf5bQuWZtMzraZfak2/dJCh0nEmg6V8NNNxBk2/aaReXP+8pXT6Recyntvjk7rQTKMYdLmg3z1BhUpRGrDuBr80lkEl04X3xPKpupUMtl5zPr/XLa/6tjtSxNq2qbVkWk3n94gtraZz/7vsd6uT3wEAAD2Hle1WVrYl8jIRH1zZAAAAAACAyNHCAQAAAABQpywVEd+SkqUiVgg4AAAAAABqzsqklZWIuEsFAYdYoUsFAAAAAACIHC0cAAAAAAD1GTQyE+1v4AwaGS+0cAAAAAAAAJGjhYNkIxvQR2U7+thSltlSmTmZqc380tkZ0xK+aTY9pgVNj+eXcjMMM32mcxuMdJCJNYV4VcJIn9nqSKXZkSrMl20vpM/MdNjT0qXbC+ta15HxTJ+Zbi9sXyZV2J50eyEFzur+3u9JMtPimmIzbJrNlnShHi2dlmu6P5Hoynik1XSm5lwXLO2q1/ERxfHgVYZvetds+WkajTRIiYBpMUOnzzTSQXolVSz61BupE81UjFbKXoJXmkavFI3lcKZcDPId5bVMFMLUKYr0mUGFqbtfnbzTbCa815lqCb3Pssb3GAAAPWIMB8UYDs2MFg4AAAAAACBytHAAAAAAANSclelWVsS/gUuZiA9aOAAAAAAAgMjRwgEAAAAAUJ8sFYlox3oiS0W8EHAAAAAAANSpS0XEAQe6VMQKXSoAAAAAAEDkaOEgaRd7JZXqlQyVTs03faY5XzZEGWb6zQiWCZSKMwrONI/G3y3Gdie77NHMlJFmLmuk1rRaUrb5sqmOwqpShfIyqazr387Umtmkd92zLZbrMs7UnGs3KPydzCRc02wms/b6JTNGWsZMynO+lrQ5n/l6wvOYMlNwFh0bXulevdK7FqXj9H6/fNO9VolznbY0mUGP60ywdJnmO+l49z3nNFNAOj/3ZopFy8gAlU16f27N+eLAPCaDC/a95Pws1Gdbq5cq1Eyl6SUT8XsAAED802JG3cIhHWl5qAwtHAAAAAAAQORo4QAAAAAAqDkr26WsgC0hg5dJWsw4oYUDAAAAAACIHC0cAAAAAAC1l+lSlhVtCwdFC4dYoYUDAACoqg8++ECdfPLJasyYMapXr15q8803VzNmzCGt4mQAABdxSURBVFDd3f7NXjs7O9Xpp5+uNtxwQ9W3b1915JFHqk8//ZS9BQBAg6CFAwAAqKqFCxeqbDarbrnlFrXFFluoN998U02ePFmtXr1aXXXVVZ7LnX322erhhx9W9957r+rfv78644wz1BFHHKGee+459hgANAEr0x15CwcrS5aKOCHgAAAAqmrixIn6kbPZZpupRYsWqZtuuskz4LBixQp1++23q7vuukvtt99++rU777xTjR07Vr3wwgvqO9/5DnsNAICYI+DwTZ773KNsqSjyxiZqnk8+kY24vEzEfa98yk4a71cyXfi7NVnYgVaLfWdmk4VD3W8/Z5OFdVktluvrTuZ8QZnLZBzLe60rmfU+ThIZ72nJbLLkMeVfdvnHkd/xEObYq+bxZbJaKv88m9vqLM+r7i3G/is6Po1l/I/DYNvnV0Y9OuKZx57zs1CpRJ06Dfp9nkoJdR5qUBJQGDRokOf0V155RaXTabX//vvnX9t6663VqFGj1Pz58z0DDl1dXfqRs3Llyoi3HAAQFUuP4RDtjQktHOKFgAMAAKipxYsXq+uvv963O8WyZctUW1ubGjBggO31YcOG6WleZs2apWbOnBnp9gIAqoOAQ/Nj0EgAABDK1KlTVSKR8H3I+A2mjz/+WHevOOqoo/Q4DlGbNm2abj2Re3z00UeRrwMAAARDCwcAABDKlClT1Iknnug7j4zXkPPPf/5T7bvvvmqPPfZQt956q+9yw4cP11ksli9fbmvlIFkqZJqX9vZ2/QAANMqgkVF3qVgXaXmoDAEHAAAQypAhQ/QjCGnZIMGGXXfdVQ/+mEz6N7KU+VKplJo7d65OhylkoMklS5aocePGsccAAGgABBwAAEBVSbBhn332UaNHj9bjNvzrX//KT8u1VpB5xo8fr373u9+p3XbbTafBPPnkk9U555yjB5fs16+f+ulPf6qDDWSoAIDmkM12qoQV7S1p1qKFQ5z06IBDLufrus74j2Ad+ywV2dpkERBW0n0EeCvpPdK7PftEiCwViWizVPjxyjiRtMJmqUjUPktFtmdlqTCPPfvrifKXCXh8lrNc0DIaLaNDrb83a1Gn3Pko6pzk9fbkk0/qgSLlsfHGG9um5eoqGSmkBcOaNWvy06699lrdEkJaOEjmiQkTJqgbb7yxrHXnyidbBQCEk/v+rMa5KWNlGqJMhJewmu2qpgxLly5Vm2yySb03AwAAGxno0HljjnDee+89tfnmm/P2AUCMzk2dnZ1qzJgxvlmHKiGt595//33V0dFRlfIRXI8OOGSzWT2A1QYbbKBH0q535FCCH/JBlmajzaqn1FNQ1+bUU/ZrT6ln3Ooqp+Svv/5ajRw5suQYBwhGBp0cOHCgHvtBumk0mzgdv9XS7HWkfo2v2fdhtc5NEnSQwYGrQdIqE2yIhx7dpUI+MHH7BUm+pJrxi6qn1lNQ1+bUU/ZrT6lnnOrajDfF9ZS7OJb3NQ77t9mP32pq9jpSv8bXzPuwGucmCQgQFGh+/HwCAAAAAAAiR8ABAAAAAABEjoBDTLS3t6sZM2bo/5tZT6mnoK7Nqafs155Sz55W156o2fdvs9evJ9SR+jW+Zt+HQCV69KCRAAAAAACgOmjhAAAAAAAAIkfAAQAAAAAARI6AAwAAAAAAiBwBBwAAAAAAEDkCDjV2+eWXqz322EP17t1bDRgwwHWeRCJR9Ljnnnts8zz99NNql1120aPhbrHFFmr27NmqEeu6ZMkSdcghh+h5hg4dqs4991y1bt26hqur06abblq0D3/5y1/a5nn99dfVd7/7XdXR0aE22WQT9etf/1o1ohtuuEHXV+qx++67q//93/9Vje7iiy8u2n9bb711fnpnZ6c6/fTT1YYbbqj69u2rjjzySPXpp5+qRvDMM8+oQw89VI0cOVLX64EHHrBNl3GEp0+frkaMGKF69eql9t9/f/XOO+/Y5vnyyy/Vcccdp/r166c/2yeffLJatWqVaqR6nnjiiUX7eOLEiQ1XT4Q//zgFOSYauX5BPttxEeazt88++xTtv1NPPVU16rny3nvv1ecdmX/77bdXjzzyiIqzcuon13HOfSXLxVWp84mbRrx2BaqFgEONdXd3q6OOOkr95Cc/8Z3vzjvvVJ988kn+MWnSpPy0999/X9+k77vvvmrBggXqrLPOUj/60Y/U448/rhqprplMRtdD5nv++efVnDlz9BeyXBA1Wl3dXHLJJbZ9+NOf/jQ/beXKlerAAw9Uo0ePVq+88oq68sor9U3urbfeqhrJH//4R3XOOefoVFCvvvqq2nHHHdWECRPUZ599phrdtttua9t/zz77bH7a2Wefrf77v/9bXxD+7W9/U//85z/VEUccoRrB6tWr9X6Si0M3Evi67rrr1M0336xefPFF1adPH71PJciSIzcCb731lnryySfVQw89pC/GTjnlFNVI9RRyM2nu47vvvts2vRHqicrOteUeE41cvyCf7bgI+9mbPHmybf/FJZBf7rlSromOPfZYHWh57bXX9DWgPN58803VLNcCEkwy99WHH36o4irI+cTUyNeuQFVIWkzU3p133mn179/fdZrslvvvv99z2fPOO8/adtttba8dffTR1oQJE6xGqusjjzxiJZNJa9myZfnXbrrpJqtfv35WV1dXQ9Y1Z/To0da1117rOf3GG2+0Bg4cmK+nOP/8862tttrKaiS77babdfrpp+efZzIZa+TIkdasWbOsRjZjxgxrxx13dJ22fPlyK5VKWffee2/+tbffflt/bufPn281Eud3TTabtYYPH25deeWVtvq2t7dbd999t37+97//XS/30ksv5ed59NFHrUQiYX388cdWHLl9p55wwgnWYYcd5rlMI9YT5Z1rnUodE41cvyCf7bgI+9nbe++9rTPPPNNqhnPl97//feuQQw6xvbb77rtbP/7xj61mqF85n8u4KXWN3sjXrkC10MIhpqS59uDBg9Vuu+2m7rjjDt0UMmf+/Pm6KaRJIsnyeiOR7ZVmgsOGDbPVQ379l182Gr2u0oVCmtzvvPPOugWD2VVEtn+vvfZSbW1ttnotWrRIffXVV6oRyC9s0jrD3D/JZFI/b4T9U4o0NZbmk5tttpn+tU26/wipczqdttVbmr2OGjWq4estv8osW7bMVrf+/fvr5rG5usn/0sT529/+dn4emV/2vfxq2kikyat05dpqq630L8VffPFFfloz1RPRHBPN/tmOi0o+e3/4wx/0tdN2222npk2bptasWaMa8VzZSNc+Ya8FpIuMtPKULqWHHXZY/rqvGTTS/gNqobUma0HZTfH3228/3TfziSeeUKeddpr+Yv7Zz36mp8tFg3mTLuS53KivXbtW981sBF71yE1r5LrKvpK+e4MGDdJNI+XCR5oMXnPNNfl6jRkzxrPuAwcOVHH3+eef624xbvtn4cKFqpHJRbh075GbDtlvM2fO1ONtSHNW2T8SKHL2m5Z6547bRpXbfrd9an4m5YbM1Nraqo/1Rqq/NJ2XbjDyOXz33XfVBRdcoA466CB9QdjS0tI09UR0x0Szf7bjIuxn7wc/+IG+gZVAsYyRdP755+sg/n333aca7Vzpde0Tt30Vtn5ybpUf03bYYQe1YsUKddVVV+kxSSTosPHGG6tG16jXrkC1EHCIwNSpU9WvfvUr33nefvtt26Bzfi666KL83/LruPQdk1/IcwGHZqprIymn7tKXMUdOqHKD+uMf/1jNmjVLDyCEeJObDHP/SQBCLmT/9Kc/caHQJI455pj839LSSvbz5ptvrn/hHj9+fF23DfU5/9T7mGj282vQ+oVljvEg+08Gx5T9JsEj2Y+Ij3HjxulHjgQbxo4dq2655RZ16aWX1nXbAESPgEMEpkyZoke39iPNssOSmx35Au7q6tI3q8OHDy8aEV+eywA81Y6aRllXqYdzFONcvWRa7v961TXKuss+lC4VH3zwgY7se9XLrHvcSbNV+dXPrR6NUoegpDXDlltuqRYvXqwOOOAA3YR0+fLltlYOzVDv3PZLXeRiPUee77TTTvl5nAOBybEto8o3cv3lsyvHtOxjuUlp1no2smqfa0sdE41cvyCf7bjUL6rPnpx3hey/egYcwpwrva4R4vjdE8W1QCqV0j+wyb5qBnG6dgXigIBDBIYMGaIf1SIj3EoT+9wv4xIVdqZHkpGczWhxI9RVtldSe8mFRa75pNRDvpC32Wabutc1yrrLPpQ+jbl6yvZfeOGFeiwAOdHm6iXBiEboTiGk1cauu+6q5s6dm8+iks1m9fMzzjhDNRPp0iS/kv3whz/UdZZ9JvWUdJhCmu3KGA/1OC6jJE3J5UJJ6pa7CZEmoNJvOjcavtRRgi3SZ1feCzFv3jy973MX+I1o6dKlur9+7masWevZyKp9ri11TDRy/YJ8tuNSv6g+e3LeFbXaf1GeK+U9kOmS3aDe1z61uBaQLhlvvPGGOvjgg1UziNO1KxALVRuOEq4+/PBD67XXXrNmzpxp9e3bV/8tj6+//lpPf/DBB63bbrvNeuONN6x33nlHZzPo3bu3NX369HwZ7733nn7t3HPP1aPj33DDDVZLS4v12GOPNVRd161bZ2233XbWgQceaC1YsEBv/5AhQ6xp06Y1XF1Nzz//vM5QIXV69913rd///ve6Xscff7xtdPBhw4ZZP/zhD60333zTuueee3Q9b7nlFquRyHbLKOezZ8/WI4ufcsop1oABA2yZRxrRlClTrKefftp6//33reeee87af//9rcGDB1ufffaZnn7qqadao0aNsubNm2e9/PLL1rhx4/SjEcjnL/dZlFPANddco/+Wz6v45S9/qffhX/7yF+v111/Xo/aPGTPGWrt2bb6MiRMnWjvvvLP14osvWs8++6z1rW99yzr22GOtRqmnTPv5z3+us4rIPv7rX/9q7bLLLroenZ2dDVVPhDv/CMkKdN999+m/gx4TjVq/oJ/tuCj12Vu6dKmun0wXixcvti655BL9fSz7T+q42WabWXvttZfVCOdKuRaYOnVqfn4577S2tlpXXXWVvvaRzEmSHUmuDeOo3PrJcfv444/ra6RXXnnFOuaYY6yOjg7rrbfeshrxvCl1kzo28rUrUE0EHGpM0m7Jl5Xz8dRTT+VTP+200076AqJPnz46Nd/NN9+sUwyZZH6Zr62tTZ9UJcVQo9VVfPDBB9ZBBx1k9erVS9/QyY1eOp1uuLqa5OQp6ask5ZOcQMeOHWtdccUVRRet//d//2ftueee+iS90UYb6YvBRnT99dfrm2/ZP5Ia64UXXrAanaSvGjFihK6T7Bt5Lhe0OXKBftppp+nUpnJRcfjhh1uffPKJ1Qjk8+T2uZTPay593kUXXaQDYnJsjh8/3lq0aJGtjC+++EJf/Mv3lKSxPemkk2w3OnGv55o1a3SgUwKBchEvaWwnT55cFChrhHoi/PlHnufOJ0GPiUatX9DPdlyU+uxJUMGs75IlS3RwYdCgQbpuW2yxhb7ZW7FihdUI50pJ6Zn7Ds7505/+ZG255ZZ6fkmx+PDDD1txVk79zjrrrPy8cjwefPDB1quvvmo16nlT/pc6NvK1K1BNCfmn3q0sAAAAAABAc0nWewMAAAAAAEDzIeAAAAAAAAAiR8ABAAAAAABEjoADAAAAAACIHAEHAAAAAAAQOQIOAAAAAAAgcgQcAAAAAABA5Ag4AAAAAACAyBFwAAAAAAAAkSPgAPRQY8eOVf/1X//lO88+++yjEomEfixYsMB3vrPOOivybTzxxBPz63/ggQciLx8AgFoJek4Ng/MlgLgi4AD0QGvXrlXvvPOO2nHHHUvOO3nyZPXJJ5+o7bbbTtXab3/7W71uAACaQbXOqZwvAcRVa703AEDtvfnmm8qyrEAXPL1791bDhw9X9dC/f3/9AACgGVTrnMr5EkBc0cIB6EGkCed+++2n9txzT5XNZtWoUaPUb37zm7LKWL16tTr++ONV37591YgRI9TVV19dNI+UPWvWLDVmzBjVq1cv3ZLiz3/+s22er7/+Wh133HGqT58+upxrr722al0zAACI2t13363PcWZLvJNOOkntsMMOasWKFYHLeeyxx/R5ecCAAWrDDTdU//Zv/6befffd/PRNN9206Fy90047qYsvvjiimgBA9RBwAHoIuXjZe++9dcDhe9/7njriiCPUlClT1Nlnn11WX9Jzzz1X/e1vf1N/+ctf1BNPPKGefvpp9eqrr9rmkWDD7373O3XzzTert956S6/jP/7jP/RyOeecc4567rnn1IMPPqiefPJJ9T//8z9F5QAAEFfHHHOM2nLLLdUVV1yhn8+YMUP99a9/VY8++mhZrfMkkC/nxJdfflnNnTtXJZNJdfjhh+vgPQA0OrpUAD3EqaeeqoMMv/jFL9Tuu++ujj76aN2aQIIDcrMvv5aUsmrVKnX77ber3//+92r8+PH6tTlz5qiNN944P09XV5e++JKLrnHjxunXNttsM/Xss8+qW265RQc9pHWDLHfXXXfly7nzzjvVyJEjq1Z/AACiJIM/Xn755erf//3fdTeJ66+/Xp9PN9poo7LKOfLII23P77jjDjVkyBD197//vS7jJwFAlAg4AD3AsmXL1Lx589Tzzz+vMpmMeuONN3SgQX5FaWlpUW1tbYFbSXR3d+uARc6gQYPUVlttlX++ePFitWbNGnXAAQfYlpXldt55Z/33e++9p9LptNptt93y0+XXILMcAADiTro/bLPNNuqSSy7Rrf623XbbssuQQZynT5+uXnzxRfX555/nWzYsWbKEgAOAhkfAAegBXnjhBX0BI60YFi1apLNUyN8ffPCB+uqrr9Qee+wR2bqkFYR4+OGHi37laW9vj2w9AADUm4y/sHDhQh3MHzZsWKgyDj30UDV69Gh122236ZZ+cr6Wlg0SqBfy44AM9GySoD0ANALGcAB6gNxFS2dnp3rttdf0hY20TJAxFuSiZvvttw9Uzuabb65SqZT+FSZHAhb/+Mc/8s/llx4JLMgvM1tssYXtsckmm+S7WEg5L730Un45GWDLLAcAgDiTcYe+//3v666G0j3woosuKruML774Qv8QIN0dpYyxY8fq86pJuleYA1OuXLlSvf/++5HUAQCqjRYOQA8gYym0trbqJp/SAkFu+P/zP/9T9zd95plnApcjmSlOPvlkPXCkjKQ9dOhQdeGFF+pfX3I22GAD9fOf/1wPFCm/0sjI2xJMkAEi+/Xrp0444QQ9j/wv5UjgQ8qRwbakHOkTCwBAnEkLwUMOOURdcMEF6thjj9XnVTnXShBil112CVzOwIED9fn01ltv1RmbJFg/depU2zwy2PPs2bN1SwjJZCHdL6Q7JAA0AgIOQA8gLQtkEKrzzz9f/0oiwQcZZ0Gagu66665llXXllVfqoIVc+EjgQDJdONN/XXrppfoXGRknQsZrkAskuQCTC7Oca665Rg9kKf1fJRBx3nnnqY8++kh1dHREVm8AAKL25ZdfqokTJ6rDDjssHxyQsY0OOuggfZ6Tc2tQEmi/55571M9+9jPd4lDGMrruuut0muicadOm6RYNcr6U8Y7kHEsLBwCNImE5O4UBaGrSokB+KZHUmKXIBY+M9eDM/10NkhZMxny4+uqrdSsKk7R6uP/++9WkSZOqvh0AAFRDLc6pnC8BxA1jOAA9yNKlS3Xf0HLSbN144426K4VktoiSjCVx991368wX0gT1uOOO06/LL0Y50gJC1g0AQDOo1jmV8yWAuKKFA9CDPProo+qoo45SX3/9daCxEj7++GOd0UKMGjUqcPrMoAGHH/3oR3qwLClXunZINwtzAMvPPvtMD44lpG9rnz59Ils/AAC1VM1zKudLAHFFwAEAAAAAAESOLhUAAAAAACByBBwAAAAAAEDkCDgAAAAAAIDIEXAAAAAAAACRI+AAAAAAAAAiR8ABAAAAAABEjoADAAAAAACIHAEHAAAAAAAQOQIOAAAAAAAgcgQcAAAAAABA5Ag4AAAAAAAAFbX/D7DDKw69JutxAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.close('all')\n", "fig,axs = plt.subplots(1,2,figsize=(12,6))\n", "cmap = 'turbo'\n", "\n", "ax = axs[0]\n", "ax.pcolormesh(phi*fp.RAD,r*sim.UL/fp.AU,gasdens_plane*sim.USIGMA,cmap=cmap)\n", "ax.set_xlabel('$\\phi$ [deg]')\n", "ax.set_ylabel('$r$ [au]')\n", "\n", "ax = axs[1]\n", "c = ax.pcolormesh(x*sim.UL/fp.AU,y*sim.UL/fp.AU,gasdens_plane*sim.USIGMA,cmap=cmap)\n", "ax.set_xlabel('$x$ [au]')\n", "ax.set_ylabel('$y$ [au]')\n", "ax.axis('equal')\n", "fp.Plot.fargopy_mark(ax, frac=1/4)\n", "\n", "axc = fig.colorbar(c)\n", "axc.set_label(\"$\\Sigma$ [g/cm$^2$]\")\n", "plt.savefig('gallery/fargopy-tutorial-animations_0.png') # Save figure\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Notice that we have used the units of the simulation to convert all quantities from simulation units to physical units." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now, to create the animation we repeat the previous procedure at each snapshot:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "100%|██████████| 10/10 [00:00<00:00, 288.47it/s]\n" ] } ], "source": [ "plt.ioff()\n", "fig,axs = plt.subplots(1,2,figsize=(12,6))\n", "cmap = 'turbo'\n", "\n", "# Create a Celluloid camera\n", "camera = Camera(fig)\n", "\n", "for snapshot in tqdm(range(sim.nsnaps-1)):\n", " # For each snapshot get the desnity and meshslice it\n", " phi = fields_spherical.var1_mesh[snapshot]\n", " r = fields_spherical.var2_mesh[snapshot]\n", " x = fields_cartesian.var1_mesh[snapshot]\n", " y = fields_cartesian.var2_mesh[snapshot]\n", " gasdens_plane = fields_cartesian.gasdens_mesh[snapshot]\n", "\n", " ax = axs[0]\n", " ax.pcolormesh(phi*fp.RAD,r*sim.UL/fp.AU,gasdens_plane*sim.USIGMA,cmap=cmap)\n", " ax.set_xlabel('$\\phi$ [deg]')\n", " ax.set_ylabel('$r$ [au]')\n", "\n", " ax = axs[1]\n", " c = ax.pcolormesh(x*sim.UL/fp.AU,y*sim.UL/fp.AU,gasdens_plane*sim.USIGMA,cmap=cmap)\n", " ax.set_xlabel('$x$ [au]')\n", " ax.set_ylabel('$y$ [au]')\n", " ax.axis('equal')\n", " fp.Plot.fargopy_mark(ax)\n", "\n", " # Take a snapshot of the figure\n", " camera.snap()\n", "\n", "# Create the animation object\n", "animation = camera.animate()\n", "plt.ion()\n", "plt.savefig('gallery/fargopy-tutorial-animations_1.png') # Save figure\n", "animation.save('gallery/fargopy-tutorial-animations_2.gif', writer='pillow') # Save animation\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Once you have snapped all the frames we can create the animation, save it and/or visualize it:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 89, "metadata": {}, "output_type": "execute_result" } ], "source": [ "HTML(animation.to_html5_video())" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Save it into a gif:" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "animation.save('gallery/fargo.gif', fps=5)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "And that's it!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Animations preview" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Since the animations in this notebook will not be displayed on the web we include here the animation figures which are available in the `GitHub` repository:" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "

\"Animation\"\"/

" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "---\n", "*Powered by fargopy*. For more examples see [fargopy GitHub repo](https://github.com/seap-udea/fargopy/tree/main/examples). \n", "\n", "Jorge I. Zuluaga, Alejandro Murillo-González and Matías Montesinos © 2023-present\n" ] } ], "metadata": { "kernelspec": { "display_name": "fargopy-env", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.12" } }, "nbformat": 4, "nbformat_minor": 2 }