{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import clustertools as cts\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Setup" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "As discussed in the documentation, a wrapper has been written around the LIMEPY code (Gieles, M. & Zocchi, A. 2015, MNRAS, 454, 576 & Claydon, I., Gieles, M., Varri, A.L., Heggie, D.C., Zocchi, A. 2019, MNRAS, 487, 147) to set up a StarCluster that represents a Galactic Globular Cluster, where the structural information is taken from either de Boer, T. J. L., Gieles, M., Balbinot, E., Hénault-Brunet, V., Sollima, A., Watkins, L. L., Claydon, I. 2019, MNRAS, 485, 4906 (default) or Harris, W.E. 1996 (2010 Edition), AJ, 112, 1487. Orbital information is taken from Vasiliev E., 2019, MNRAS, 484,2832 via ``galpy``'s ``Orbit.from_name()`` feature Bovy J., 2015, ApJS, 216, 29. \n", "\n", "To setup NGC6101, for example: \n" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "LOAD_CLUSTER MADE USE OF:\n", "Gieles, M. & Zocchi, A. 2015, MNRAS, 454, 576\n", "Vasiliev E., 2019, MNRAS, 484,2832 \n", "de Boer, T. J. L., Gieles, M., Balbinot, E., Hénault-Brunet, V., Sollima, A., Watkins, L. L., Claydon, I. 2019, MNRAS, 485, 4906\n", "1936 11.8971760895 1936.332\n", "-2005.21009948 -9348.18148437 -3945.46817625 -296.181213344 82.7743019402 -190.8475368\n" ] } ], "source": [ "cluster=cts.setup_cluster(ctype='limepy',gcname='NGC6101')\n", "print(cluster.ntot,cluster.rm,cluster.mtot)\n", "print(cluster.xgc,cluster.ygc,cluster.zgc,cluster.vxgc,cluster.vygc,cluster.vzgc)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Unless otherwise specified, the cluster is setup to be in ``pckms`` units in clustercentric coordinates. The number of stars is set using ``mbar`` variable which has a default of 1.0 solar masses." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Clusters can easily be viewed as they would be in the sky using the ``skyplot`` command:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABjsElEQVR4nO39eZhs21XYCf7WmWPKjLyZecf37nt672kWg9AtjI3VNiCpXbTVwlANkgc11WVkGyi7CgqXbNyiPky5Ufnjo8plt0FFlcGULGMbhA3YlhE2dpVscD/pSUhC45PeeKecIjPGM+7+45wT92TciMyIzIiMiMz9+774Ms+89xn22mvttdcSpRQajUaj0SwaxrwLoNFoNBrNMLSA0mg0Gs1CogWURqPRaBYSLaA0Go1Gs5BoAaXRaDSahUQLKI1Go9EsJFpAaTSAiHxWRP7ovMuh0WgeoAWU5kIgIq3CLxGRbmH5TymlXq+U+u15l3MUkvJ+EdnJfu8XETli/z8pIs+LSFtEflVELhW2XRKRD2fbnheRP7kIx2o0g2gBpbkQKKWq+Q94AXh7Yd0H512+MXgP8B3A1wFfC7wd+HPDdhSR1wM/C/wZ4ArQAf6/hV3+DhBk2/4U8HezY+Z2rEYzFKWU/unfhfoBzwFvGbUO+O+Afwz870AT+DTwKuCvAPeBF4G3FY5dBf5X4A7wMvATgDnlMv974D2F5f8C+J0R+/4N4B8Ulp8kFQw1oJL9/6rC9l8EfnKex+qf/g37aQ1KoxnO20kb0DXgGeAjpBaHG8CPk2oKOT8PRMBTwBuBtwF/dthJMxNY44jfzRHleT3wqcLyp7J1x+6rlHqWTDhkv0gp9cUR55rXsRrNQ2gBpdEM5/9QSn1EKRWRalObpL39EPiHwOMiUheRK8C3A/+VUqqtlLoP/DTwzmEnVUr9A6VU/YjfCyPKUwX2C8v7QHXEONTgvvn+tWzbwYht8zxWo3kIa94F0GgWlHuF/7vAtlIqLixD2uheB2zgTkFWGKRmwGnSAlYKyytASyk1LNrz4L75/k0gOWLbPI/VaB5Ca1Aazel4EfCBjYIWtKKUGmq6EpE/NeBROPgbZeL7LKmDRM7XZeuO3VdEngBc4IvZzxKRV44417yO1WgeQgsojeYUKKXuAP8K+CkRWRERQ0SeFJE/MmL/D6qCR+GQ3ygT398HfkhEbojIdeCHSce+hvFB4O0i8mYRqZCOmf2KUqqplGoDvwL8uIhUROSbgXeQjrfN81iN5iG0gNJoTs+7AQf4fWAP+CfAtSlf42eBXyP1KPwM8BsUHDUy7evNAEqpzwJ/nlRg3Ccd5/n+wrm+Hyhl2z4E/IXsmLkdq9EMQ4absDUajUajmS9ag9JoNBrNQqIFlEaj0WgWEi2gNBqNRrOQaAGl0Wg0moXkQk3U3djYUI8//vi8i6HRaDSajI9//OPbSqnNYdsulIB6/PHHefrpp+ddDI1Go9FkiMjzo7ZpE59Go9FoFhItoDQajUazkGgBpdFoNJqFRAsojUaj0SwkWkBpNBqNZiHRAkqj0Wg0C8mFcjPXaDTzJ4gSwjjBNg0cS/eRNaPRAkqj0ZwZQZTw8l4HBQhwY62shZRmJPrN0Gg0Z0YYJyig4lqobFmjGYXWoDQazZlhmwYCtP0IyZbPAm1WXE60gNJoNGeGYxncWCufqbDQZsXlRT8ljUZzpjiWQcW1zkxIaLPi8qIFlEajOdfMy6yoOT3axKfRaM418zAraqaDFlAajebc41haMC0j+olpNBqNZiHRAkqj0Wg0C4kWUBqNRqNZSLSA0mg0Gs1CogWURqPRaBYSLaA0Go1Gs5BoAaXRaDSahUQLKI1Go9EsJFpAaTQajWYh0ZEkNJopo1M7aDTTYS5fj4j8koh8Mvs9JyKfzNa/VUQ+LiKfzv5+64jjL4nIb4rIl7K/a2daAY1mBHlqh3sHPV7e6xBEOnK2RnNS5iKglFLfo5T6eqXU1wO/DPxKtmkbeLtS6muA/yfwiyNO8V7gt5RSrwR+K1vWaOaOTu2g0UyPudofRESA7wY+BKCUekYpdTvb/FmgJCLukEPfAfxC9v8vAN8x46JqLihBlND2o7E1oUVJ7TBpuTWaRWTeY1BvBu4ppb40ZNt3AZ9QSvlDtl1RSt3J/r8LXBl1ARF5D/AegJs3b56yuJqLxEkysS5CagedQVZzXpjZWysiHxWRzwz5vaOw27vItKeBY18PvB/4c8ddRymlAHXE9g8opW4ppW5tbm6eoCaai8pJzXVnnTF2EG1m1JwXZqZBKaXectR2EbGA7wTeNLD+EeDDwLuVUs+OOPyeiFxTSt0RkWvA/WmUWaMpsijmuklZ1nJrNIPM08T3FuDzSqmX8hUiUgd+A3ivUupjRxz7z0idKH4y+/tPZ1hOzQVlEcx1J2FZy63RDDLPN/edPGze+0HgKeB9BTf0ywAi8nMicivb7yeBt4rIl0gF3U+eVaE1F4t5m+tOyrKWe5nRjinTR9IhnIvBrVu31NNPPz3vYmg0mnOGdkw5OSLycaXUrWHb9B3UaDSaU6IdU2bDvN3MNZqpoMMLacZhVu+JdkyZDVpAaWbCWQqMi2Re0YL45MzyPdGOKbNBCyjN1DlrgVE0r7T9iDBOzmUDMY37OqmAO08CcdbviWMt/z1aNLSA0kydsxYYF8W8ctr7OqmAO6uOxlkJwYvynpwntIDSTJ2zbgguinnltPd1UgF3Fh2Ns9S2L8p7cp7QAkozdebRECyaeWUWWsFp76ttGkRxwnbTxzblWAF3Fh2Ns9a2F+090RyNFlCamXCRG4JZD8af1X09i46GNrtpjkILKI1myozSCk6rVbV6EZ0gouxYVL3JP90wTrBMg9Xy+NrKrAWiNrtpjkILKI1mygzTCnKtKogT4lhxc71C1bPGFlqtXsQzz++SoDAQ3vjYpYmF1KJqKxdZ29YcjRZQGs2UGaYVtP2IIE446IZ0ggiAm+sVtpq9sUyBnSAiQbFe9dhp9egE0cQCSmsrmmVDCyiNZgYMagW2aRDHKjXRuTamKXSCaGwHgbJjYSDstHoYCGXnZJ+u1lY0y4QWUBrNGeBYBjfXKwCYpuCYBmUnFUzjmNyqnsUbH7t0qjEojWbZ0G+5RnNGVD2Lp67UDpnYHGt8k1vVuziCadyxufMU6WLanId7czHeds2FY1E/zkET27xMbot6f2B8N/2LFINxUs7LvVm+Ems0x5B/nPcOeqnn3DEJ5C5aorlJ789ZM27qCp3iYjTn5d5oAaVZWkYJlkk+zkVrrM9CWC564zWuO/yius0vAufl3mgTn2YpOcqEMcnHOYtQOyc1n50kmOtJrrPojde47vDj7rfI5sxZcV6mFGgBpZkbwxqOcRuTowTLJB/ntBvrVi/ihZ1231NvEtv/JMLyNGMMy9B4jTs2d9x+52Us5iSchykFWkBp5sKwhgMYuzE5TrBM0sBN2liPEqJBlPDCTputVo+ya7PiWcdqZMVznaXmdxaN1ySay6y0nIuSK+y8ogWUZi7kDYdtGhx0Q9p+hGMZDzUm+b6DDdc0tYBJGuujeuRhnGCa6STajh9Ssc0jhcywc81L85s2k2guJ9VyxhFqi36fNEejBdSYXEQ79izJUz+8tNsGBNsUrtfLhxoTpY7WqOZhwjiqR26bBo5psFKyqTgWN9crY5spG52ARiegXnaouMd/lotupptEczmJljOuUFv0+zRrlr3d0gJqDC6yHXtWOJbBZs0jjBWrZZsgShDhUGOyiOaZo3rkkzaG+bkanYD7+z3Izjvu+zWOgJ5XAzWJ5nISLWeSd+M8jMWchPPQbmkBNQaL2FCeByquRdkxU+EEhegKD+7tWZlnxm3IjxNCwxrDUefOz9XoBADUy85U3695NlCTCut62QHSd2KcMmrT3fGch3ZLC6gx0B/DbBinsT8L88ykDfm0xqzyc+WC6bj3qyjoYPjYXJF5N1DjanjF+zOOeTM/90U23Y3DeWi3tIAaA/0xzI7jGrGzMM/MsiEf59zjvF/FhjzKnEesrAEaJVCn3UDNwlx4mnt/UU1343Ie2i0toMZEfwznl1n2NMc993HvV7Eh3276AMdmxp1mAzVKEzyt0DoPvfxFZtnbLS2glpxl99KZFZPcl1n2NKd17mJDbpsCjNeoj9tAHXe/hmk6MP68taPKt+y9fM3s0AJqiTkPXjqz4CT3ZZY9zdOcuyg4ig05cOj/VHCd7Drj3K9hms60TKP5/cnjEJ6mHlrQnS/m8hRF5JdE5JPZ7zkR+WS2/q0i8nER+XT291tHHP//EJHPikgiIrfOtPALxKIH/ZwX5+W+DAayhQdebo5l9B0KThvsdpz7lWs6V1a8vgCbpnnutEF7Fy3or2Y6zEWDUkp9T/6/iPwUsJ8tbgNvV0rdFpE3AB8Bbgw5xWeA7wR+dtZlXWS0/X44w+7LMvSuB8s4joZy1D7j1vmk42TjmufGKcdptbF5eyxqZsNcTXwiIsB3A98KoJR6prD5s0BJRFyllF88Tin1uez4syrqQnJR7PeTCpfB+wKHx0o2ax4iTP2enUYIDjOzjSM4lIJeEKcNchbPb9T5jprIetL36Djz5bjlOG1nS3fWzifzHoN6M3BPKfWlIdu+C/jEoHCaFBF5D/AegJs3b57mVAvJsnvpHMdJx9mK96XtR4dCCr2w08ZzzKmO240Kfjtuoz9MA6i41pGCI4gStpo9TFOIY8Vm3QNgrx3QCSKCOKFedg6FUZo02sJpNc9xNZvTdrbm0VlbBq181sz6HsxMQInIR4GrQzb9qFLqn2b/vwv40JBjXw+8H3jbacuhlPoA8AGAW7duqdOeT3O2TMN0U9Qy4lhhmjJ1U9BgOdt+RKMTjC1YR2kAR3VAwjghyDSngLShuN3ocHe/R5ik0TnCOGFrv0c3iGl0Ah7fqA4VdKOis4+ben1UI5XHXNxu+timTMXjcNzjZ9l4agels7kHMxNQSqm3HLVdRCzScaQ3Dax/BPgw8G6l1LOzKp9mOTit6aaoZXSDmM2aeyhqg1KTecCNavQGywlMJFhPogEoBff3eyQoDIS1skMYK0qOiYeJbRi4lkGoEg56ATttxWbNw7GcQ/UpTgDerHl9R4xxOgeL2lCftlwncbtfhHqfJWdxD+Zp4nsL8Hml1Ev5ChGpA78BvFcp9bF5FUyzOJzWdJN/RGXH4t5+lzTvUtqY26bBVrM3UQbbUY3esHGvRicYKViHNYCTzlkK44TLqx6OZRTmJim6QQwo6qslyrZFox3i2DFBmBAOeLfl98exDF7a7RDGirJjjj0GdlwjFcYJlmkcO6l4UmYpQE7qdn/ROIt7ME8B9U4eNu/9IPAU8D4ReV+27m1Kqfsi8nPAzyilnhaRPwH8z8Am8Bsi8kml1P/1zEquOVPGabiP02wOuiEg/cjp+T5547zfyXNSOUPPD8c3evn/w+YtTct0NqjxSLZurxVkwtfgNddWsE2jX871mottCGGisEc0tPudEFCslOy+8DtuDKx4/KhGanC7UukYGYwfGHbYfZmlACk+51HjdxfFQekozuIezE1AKaW+d8i6nwB+YsT+f7bw/4dJzYAazViaTR6BoRg5HchyUnUAhZ2NTY360I4bTxnliZeaFJmK6Wzw2NSsl2pFedDZimsVgq5aPLpWJowf1K/I4P0J48P357jOwXGNVHG7UvTHyEBxdbU0dEzsOKYR3/C4cbNx0qCcdwelcZj1PZi3F59Gc2rG0WwcK00EONgo5TmpiprDST+4wZ73VrNHox2y0+4BwnrV4Xp9fNNZkJUlb0whbVijOOkfmwucUdHQHcvg8Y3qkb3c4v1p+9HE9R5HiD0wQ6ZjZOmEYHWi+33a+IbjRJifZRqUs+I8eBlqAaVZek7TYOU5qQY1h2EcN54y2PPulmzuH/hUvbRBTqNBCGXHPHYu1qADxOWqx/2DdLwMYC3LvJsfO2mOqiJ5Q6YUfc/DRic4tcPDYAOZj//1x8hK9onGLaY1LnmcBjZuGpRFZFGdVyZFCyjNwjKtJIJHMcmxx5n4BnveFdfi/kGPnYMA0wLPMg9lDz4q91EYJ9RKNuVM2+iEEZ0gZrVs9681qbnpuDGtXhBjmjIVjWFUA/n4RpXNmte/Pyc9/2lMS5N0aJZ1nOm8eBlqAaVZSGaRRPCozLbwIAbdaT7kwZ63bQhra2mDPJg9eFQZ237E7b0uB92Ag25qGsxNhTstn6ur3onc7Z/bbvXHovKxn2JDls8TO6kLfpFRDWR6Hqt/r+fBJIJnWceZzouXoRZQmoVk2j3AowTeOMIwiJK++Wsjm0t1XFSEwTGMtcwTbFSjmJejE8TstHs8slahG8bUyzZ+lPCKco2DbpjNZZrsXrT9iLv7PUqOyUEvpOJabNa8Qw2ZYxps1lPTo1JM5IKflz9v9Ec1kPMwPZ3GpX9ZWWbtr4gWUJqZsGiJ7I4SeEdtyzWavLHOvboc0zhSwxg2hnGcSavtp2a8smuy0xK6YUzZMVktOWw1e4RxQtkxx06L/jCpU8K9/R4lO63r4xvVoQ1Z248ectIYJrRzp4phc8py78AiZ216Oi9jMSdhFkL4rB0vjn3TReQy8M3AdaBLGkn8aaWUjmevGco0GoWi+/NpynFcjx5GC8NDGk2rxys2a1xe9VgrO5Qd61gNw7EMNmse+93g0PUGG/YwTifQbrVSE95OK/f4KxXSaxwWIkXHhmHOFsMcFMq2zf1Wj2YvpNkN6YYR9ZJDxbMeOn7QSeNGvfzQvX1uu8Xd/S4grJZsPMc8NH5lm8ZDThdHPYdFSymvOcw8hP1IASUi3wK8F7gEPAPcBzzgO4AnReSfAD+llDqYaQk1S8c0G4ViAzdJFPJhH9Mok8coc0hej9WyzU7L56AbUs4a4XHnMqXzftJG/Oqqx/V6+VC8PD+M0hIqwXUMXrFRpRvEXK+XWKs4h8o4aJIM4oT7+700mkQ2MTgXXsXxpvyaux2fr241udvooRJFxbUp2yYrJefQuBSkQi8/bxCnTh3FejU6AZ0gxnOsdLwKDo1f5YJ33AC4s2r8zstYzCIwD2F/lAb17cD3KaVeGNyQxdH748BbgV+eUdk0C8g4vdxpNQqD84ryKOSDMeOOO7bYOI7zQR3STrLlq6te/5r5Pt0gohvE2KYMNfel86oUlmkQJQmdIKYTPJgLFHYVL+/1WCnbWAIKm24QH2vGy+sG0OxF1EoxOA+cPLaaPV7c61DzbHpBlDlAKDzbZK3sEidQLlm4psFeJ8AyjSxOYRqnL4iSvjOHgqFpPII4Ybftg1JYpkm5arJaskFgtfQg6sK4AXBP2vgd9z6el7GYRWAewn7kV6CU+pEjtkXAr86iQJrpMcqt+KQf67i93MFGAU7mDVb8IPIo5LZp8NJum2LMuGHnnORjGgwfBGBlxw9qba1exAs7bRIU9w96XKq4gHC70ekfU8w4i4KXdtNzx1XFY5cq/blAvTDGtQwsI3Vfv7LqHTLrHXVfun7Ms/ebPLfTZq8bsFl1qWaNeyeI2Wn6eLYJpNEhUrNeBAL1is1G2WWt6tAJ4kzYqf69yLUvUA91BHJBUi+n2l01E6Qv7LZ5fqeNa5s8ulbua2ObNY9OEFF2jq/TpI3fJO+jFkynZx7CfpwxqCvA3wBuKKX+mIi8DviDSqn/deal05yYUfmJTmNGmaSXmzcKpzHdFD+I3KtsMKbecZ5043xMxXptN9P0Y/lk3OJ8pSBKeGGnzVarh4GglGK1bNPxY8JYPTSB17EMrq+VCJOkP5/JsR/MBcrHncJYIcCTl2tUPYtWL03VUXYsqt6DT7So2UVJQpjE1CsO9bJNFCvuNHp0gohH18tsVF1c08Ry4CtbLZRSiKH4Q09uUHatvvB8Ybed1r2SusI3OkHf2y81NR72dixGsnBMg9WSw5fvN/nS3RZBFLNWcbhUdg9pc4rMS9Aa7h15XOzCYeRmxjznlR5fOhvOWtiP4w7088DfA340W/4i8EuAFlALzDBhApOlgBjkJL3c09qtix+EY42OqXcaivWyzXSwZVgdwzjBNIWyY7HfC1FJ2lAOHlM091Vci9WS3RfQeeObB6WtV5xDjXKrF/HM87t954Q3PnaJqmc9NKG27Uc0uzE7HZ8gTLi5UeJS1aGzG9LxYx65VO5rL/vdkPVqiZ1Wj7Jr9Rv0rVYPLxNEYRJzu9GlE0REcYzChL4h8XCHB6DqWNj5PCoFFdegG8W0g7CvcR737IfNzRrHQ3FwDA4OmyHPgll5s52H8ETTZBwBtaGU+kci8lcgNe+JSDzjcmlOyShhchob8klU/GnarfOGfVhMvUFOk/IchmfCtU0DxzQo2RZxori66rFacvqNalHTy6+7WfP65rBi3LwHgurwNTpBRIJiveqx00o1oqpnHWrs99shn365QaMXYQpcX/Oouen41XrV5Xq9lN53AVRa1tt7XaxM5nz5XpMgTjjoBjyyVmE39NltB1yueRz0QlZLDqZh4JTSlCR5Tq1GJ2SjlmpH2y2/Px7omAb1ikvJsbhc83jqcu3IMaic4tys4hjYcQyaGdfKzpHZgqfNJO/WJALnIrvEj2IcAdUWkXWy7pSIfBOwP9NSaU7NKGEyroA5KurCpBrQNOzWg+U57jyTam6D5xxlNtyseRx0m3T8dCzq6mpCxa32j0/TrachidrZPnl6+WFzhQavYxsGQZhwb7+DY5qUnfQTLQp6P46xLIPVksXdvR4tP8IyBc+2qLjmoetEcULNs+iFCs8W7ux3aWRu770w4cv3W0RJzH4noOrYxDFsVF06Wdij240OnSCNbNHoBFypl1ivuJQyR462H7FZc7i+VgIOz/Ua79mrQ2Ng4zA4ufgshROM/25NKnC0S/zDjCOgfgj4Z6Su5R8jzcH0n820VJqpMKwhH6dxn3ZP7rR265OUp5jmfZrmH5HUgaKWme2KEbmDKOH2XpeXG13uHnTZqLh4hYa8E0SHGqC0ETL6c5mUgkY34Fq9TC+KeXKj2h+DKjb2rmnwla0WYaTAgHrJBmVQdk0s0zh0ndt7XbYOfKolm61mkEWPMGn1IsqOiWublB2PIE641+xRL1vsd8N+uV/c7XBnr4tpCI6dCszr9dLYE5CPevYV1+LqaokwTgPHjjsB+TSdnmmY0Ma1Cow7FWGcuXoXlWPfCKXUJ0TkjwCvJm0fvqCUCmdeMs3cWKSe3EkGw4tp3uNYsVmfPDTQKNJGQ7F1kOY0qmxUD3kq7rR7VD2TVi9mo+oSxA9SY5SdB4IpipND0Skur3rEsSJRpFqQlXrf3T/oYRsGjm30x7SUgqpr8tx2kyCM+dSLDdZrLr0g5vqlEmslpx9VvdH1ub3fox7FdHsxJcsijCNCBUoZNFo+B6bQaIWsVVyqnoNpCr0g5t5+j1YvLW/ZNYmThLKdzgPL54KNa7oapY0flwpkFCfp9Eyr4zWugDxO4EwyV++0LOvY1jhefGVSLeoxpdT3icgrReTVSqlfn33xNPNgWj25034UJx0ML45R5J54k5bvuLKnBikD1IOTh3FCL0qoWgaebVB2LTYL42XwYNwkjBOe224TK0U7TO9zO4jYbqaODH6Y8IJqIZbBXjPgqatVXNPk5nqFME6olRyu1suEkeKg51N1TBq9HmovwTaEjarLvYMee02fg26QGtKUEKmYph9T8yyavQglikdqJaIEyrZJGMUY2Fxe9bBMg0fWynx1p9W/7zcvVfrlH2cKwaggtTln6RU2zY7XOOU+TpCdZq7eJCzz2NY4OvXfAz4O/MFs+WXgHwNaQJ1TpjFuNI2PYpzB8FYvNZ0VtYxxBGxR+MWx4uZ6pW9OO6rs+eTbmmdhmQbIgwmyjU6AJdDqRtxYK4YpOuxu3/XTYLCffKFBIglRCO1ehGUYiCg2qh4H3YDdtk/NstgPQl7c6VL2TII49Ro86AQ88/wevTDCNAzKjsVGzaMbJTx7v8VXtlsYCDudkF6UELVDlEq4UnP46naLSxWH1bJLzclMjKSzkm3DYKPqpq7o3YBuGLNedhEB0xSe3WoCQq2URpAAHpr/VeSkjhD5c5hmr38eJrSjBNlZhX1aJIvIpIwjoJ5USn2PiLwLQCnVERnVJ9WcF07bsz3qoxj34ztuMDx3yfaThL1mwGtvrFBxrLFMJW0/db9u9kK6QZq99nXXV3Gso92jlYK9ls/9lo8Aq9edvku1ZRq88uoqB92Q6/XSoes22gE7LZ+aZ/P8Tpvdto8pisurJXb2A9pBgqiYRBTPbzepeQ57rYCdVsizd5scrARcrrhUbBO75HB51Uu1NMclihVRlLqeN9oRNc+k5YeISnh2q0sUxazVHDzLodELMcWg5trESYJt21Rck0Y35G6zRxylUS7iWNEL476jhWtZREnMczttUAavvFLBs1MBPCqBY/9dSBII0r9hFodwHCedUR2IkzKtjte0BMeo8kxb41nmsa1xnnggIiUeePE9CfgzLZVm6Rn1UZzG/Xtwv9wle7Vks9PMJs/CsaaSfIzq9l6br251uLLmYbZTIVLxrH6Io0YnII4VquBgJgLXL5W5sV6m48dsrqRu10GYcNANszlBh0MVtXoRn325wf2Wn0ZuEEXJNNnphJiWRdlLx7V2Oj2qrkM3TKhXDNYqHoah+PJ9uLfvc3/fx7FNHlmD23tdJBEiUbzc6FCyTZQyEJXw1e0e+50AxxYcEyzLwo8UcRJSsSt0HJM7Bx3a3Zit/R7P7bR5fruF6zh0egErJYuDbsR2J6BkG4RRQjmb91R2LWzDoNmLsUwDI07d6tOJy6Pfg16YAIqtVo9WYB377MM4yVzhQzpBGlj3qSu1qQiF/PzF5XGYhalsWEdw2hrPtDxp58E4AurHgH8JPCoiHySNbP69syyUZvkZ9VGc1v272IMtOxYGwn43xMAgQY3VQ8y1nZsbVXZaUd+EeGe/Sy2wEaBecriz38U0ha1mjzwKQj4XSgFlB7YOfF6Ou9xpdIhihWkJj69XDl2vE0QYpvCKjSpfvtfEsoRH1sqYBlxeKaGUYqcdsNOCjh8TJQlPbtb4cqdHqxdzbz+g5Jgo0sgTe20/NQtGIURQdW1Kno0YijgxCKOISxUP1zRoBz4iJlfrHnGcsNcNMS3Y2w9xHIudTsh+LyKOIY4j/Chmv5uO2xkCnTDBEoOrqx5JBK5r4FkGKyWbkm2hgDuNDtfqpUP3KUfkQUO+30mjgIzz7G3TII5VGibJtTFNmYpp6jRC5qxMZbPQeM5yrG+aHCmgRMQA1oDvBL6J9Jn+JaXU9hmUTTMjpmGmGOccwz6K03x8g43LZs3j1ddWCOOE8qNWfwxqEhf0jZrNquvgWHLIJTxM0qgReTTvYviiPA3IdivVPmxT+MLdJpdWbEqRxX43DVWUmyTLTjpZ9kv39wnidMym0QlphzG9KGbroIdC0eqFmEaCGSg++eIe9/a72JbBfjegFxkYYrB74BMLuNl9u98JMETxmRf2uLLi8uSVGiXHIk4SlBLWKh4QE8cJB72QbpiGIzroxbhxjCgDTFCiMDG5eSmN6dcLE1QMcZRgO9Dpxbz2+irX10r9Z7bXCZDsOVdcq+9xWTTF9gW6SjMKkz37KE6yKByjQ1XdzAR9/hym0VCfRsiM++6e9vtaZo1n2hwpoJRSiYj8ZaXUPwJ+44zKtHQskwvnac0UxQR+Rw2Oj+I0H1+xcSlGN0+F1fiRyreaPRIF202fq/USrmlybbVEoxs8aHwMg/v7PfwkQSWKy1WPivvgPFvNHl++1+Sr2y1WPIcgitlqKBIC9johJceg7UfcWCtT9SxefXWFRiek5gpfvHtArRzS7kU8dqnS1+jaQYIhCbYJu22f7bbPSsmiGyUoERwjoexafO7eQTZmpegEIY5lYxiKjZrHk5drXK553D/ocP8g4Pn7Tdp+TNm2uFYvk5DQDRSrrtD0E8IwpGMKK2WXJzfLbK541LxUO3p5vwMKLlUcFPDYRqWfAqTo/h8niu2mz17bp7uSalLX6+W+mXVY4OCtZo+9TnBkGpWqZ/HUldrCOEqM8+5O05V90duSs2AcE99HReS/IY2/185XKqV2Z1aqJWLZXDhP04M8nMDP54nL1SMDto7ipB9fsXHJo5sPxho8rjF7UH+z3/MPo3TQqdj4hHHCWtVhtxMShhF3Mm1GJO+QKOoVh6thiSRW1Mo2ohTtMMY2ySbCWv17Y2cRyxudmLsHPTZXXF7qBnzqxT3iWFFyLNZKDo1uQC9UWEbMfividqNLsxNQqZfwbAM/UnimieUKdxtd4kjheiAq4X6zywt7LVwr9Zj73J0mu90eCsGzIizDpFoyqZdsNlbLSLPLvijWKi69KOZu08ez09QckHrw7XcDDFNwjcMaTFGT7AYRe+2Qe/s+rmVy7yANgLtasvtBiovHhbGBlc3pGuxoDH4/026oT6udHFeeZfaYW0TGEVDfk/39gcI6BTwx/eIsH8v2Qp6mB5nXdaVks9Pqsd9JE/idlVdQsXHJY94VA7SO01HIzXuJgjhRvLTTwTIF2xQq7uFgpQaCShJWSg4J6lA+KlBEseJS2ca2TDZWXPY7IbudgJpr4ccJvSBOJ8t2AgDWay6GGeBaJkGsMBFKtsXlNZe7Bz0cE1zbyILLKkxT4XdiOkHIdstgRxRV28SPYlY8B+hRcg16fsT9dpd7zYBnt9o8vl6mXna53/ZpdXyUCFdqHmVPuFL1sI1Uu+vF0OiGhDGIgqod0Q6itG5JQr3i0u6FbGShjQZ9dx3LoO3DQTdCodjt+pTbJmGicCyDTlb/th8dei7HdTRm/f3MUjtZZo+5RWScSBKvOIuCLCvL9kKepgeZ1zWME66ulo5NGjgLio1LMQ36cR2F3CR1p5E6PhgivPrKCo1uyErJzuY3JYVzHx4DyaM8hHFCx4+5US9xvZ5e//Zely/eO6AXxZgGPLZexbENgjjmY1/aRpGwVnGouha7seLKist+J8B1LPxQ8fl7TSqegW0bGL5gGELXjyk7Np4TUbItqrZJO4zAEFpBzHa7SaiSNHZfFCNiYIgQBDHP73TYbgWIilkp2wRhTJyk9T/oRNQrNpeqDiuuxT2EbhhiYvDCToc7+z1c2+CpyzU2ah6WaaLUwxOkc805zIS1ZQhxomj7MUkc8+Jum5KdCvPBVPDFrLqDHY1F+H5OY7LX40fTZZxIEt85ZPU+8Gml1P3pF2m5WMYX8qQ9yEWr62A9jpr0+Nx2K03u1/Z59FKV9apD2bX6DhDFQXugL6zyMZAgTHjmxV22mj0EA8cSXnllJU174QeYhrBastmseTy+mQq2L9xtYphQsuw0Rbof0fJjqmWbg3ZAzTU58EN6QUwUKWJJIAHHMLjTCWn1ArZa6VygZjYh1w9ikkTR7gR0I/CjCMcU4jhmuxcRhlBy4UrdI04STMukXvNY81xKtsVeNyJohkRxwvV6lcfXa7R6AXudHt0w4Yl6DVOEMEkQEV51tcpmzaVs52GaDge9PeiEiAhRnIZ2euRSiVY3ZMWzeeRShU6QakmDLvuDHY22H03tvTiNgMkTUuaOGScx2evxo+kxjonvvyCNIvFvsuU/ShpZ4hUi8uNKqV+cUdmWhov0Qi5qXY8Snnn0h6pn0eyGbB/4VDKPvVzjygftt5ppSKXBjLoIXKq4WKaBZaQrGp2AL9w94P/80hb3DwJWyxZ/4LFL2FdX2Wr16AQht3e7RCohiRVV1yJSikglNHsxB52IRCW0snBMQaIIgphuFNHoRERRTKMTE8cxvYMI1xI+H0Q0g4AwAD+GOAGD1CUcI10+6EKc9Cg7FmXT5ErNo9OLuZ2Fi1or2dRcm8urNl+80+IL9w9wTSFOErYOfLpRmh/qOc/itVfrtPyIL9w9AAVRHCOGQdWzWHEdPv1yg5WSTaIUGxWXqmtnHntWXzhtVF22W/5DLvtFGp0Alf09zTjuOGPCowRYMSFl2bH6mvUivu8XhXEElAW8Vil1D/oZdv8+8AeAfwdceAG1KCyTN+EsGCU80/BHaZr1th9RctK4d+1emgI9j/5tmwa7LR/LNLhWdw4N4EdZiCGUIooBFC9sd/jsy/s0OhGuKdQrNpWyQ5hkUSUurxDHCfvdkIpj0osUW80ufqB4ZK3EgR/QaMckhDTaAd1eTDeOub3fwTENXtrrsNVOUgEEhEk62bbtQ5DAIZ1DAYUsbYEPCRHVkkXFsrh2pcTtRif1uGsFVD2LL95psd3qEYQJRmJQr9isli3WTRfXMnh5p4Nrmbzqygq7nYCOH/LyXpc4FmwrzQ6cJHCp6mIZwkrJpuJY1Es2m1WvP49su+WTKCgNuOznTHMcdxxT71FhrExTKLs2HT+k4lgLYXK8yIwjoB7NhVPG/WzdroicKKq5iPwSaXR0gDrQUEp9vYi8FfhJwAEC4EeUUv96yPF/E3h7ts+zwH+ulGqcpCznhWXzJjwNkwZ5daw0arYI7DZ9kiThq9sttps+lSx8ThwnmFljZCj64X6KA/ibNY/NqpeG7lGw22qmHmydHvvdGEyDVi/ANgy6QcTn7+zzwm439ejzTFY8OztXKhS6oSJOFCXL4qvtDge9EM+EbphGUfCzrMGQPtNumDp2BJnWdBwGqdDthTFxS+GHCSIGfhhjkMbV6/mp0N5PUg3voBtiGTEJcLnmstfy+aLa5+PP7bHT9tn3Q65Xy2yuehgChii2mz0u11xuXipjG0Yq9MMI0xTqZYetZo87+10cy8BAuFE/7Nk3zflFx53rKAGWz9la8SwqdhqY97x9Q8vWiR1HQP22iPw6aYBYSHNB/baIVIDGSS6qlMo9AxGRn+JBAsRt4O1Kqdsi8gbgI8CNIaf4TeCvZNl93w/8FeC/PUlZzgvL5k04KfmHNZixtiiIjxPSbT+mGyr8OKTi2IRKYRsGJcckStL/r9Y92n7EWtmh7FgPDeDnc39e2GrR6Ib4seLmpSoiwmuv1bi2WgGBZi/i0y/tc3u/hykGr7xS5qAb8tx2k8/faeHZBp5pUHUsDEvhhyEiigM/SdOmJxDGYBug0qGpVGOK079WVsdRaf4csvxSHZ8kgQSFZ5jYjsKPEr663SZKEtZKDn6UEKk0CaIlBldXXLY6AakcNvCjhJJtslF1QKAdBjht4dGNCmXPwjRSIb7XCdja75FkUjUXrgbCtXqpP5l3mDfgtOYXOZbRT3Vfdh524MkF2LAwVos2xjptlrETO46A+gHSSBJ/OFv+BeCXlVIK+JbTXDwLOvvdwLcCKKWeKWz+LFASEVcpdSj2n1LqXxUWfwedQHHpvAknofhh9bJMr8NyQ40S0rkHn2UYPL5e5rmdTj/zbKuXRo24uuJhm+m+xcC0wzwFRWC3G2bR1QM82yBR4FkWTnbfu37IbidkpxUQJwkl22C1YtH1ExKlaHYj7kcRrp3O/0EJlqmwTDABwxY6gcI0U+Gksl9usjjOpUABjQBUEhOE3bQhNgViwTAUFceg6tmUHYNLFYdSlFAv2YghrFU9rl4qU7FNbjd6NP3U9KUiYa3kcKNeYbVicansUivZtLpphIrddkA7illxLTwnjdlXFPS5N6BS6XuaJ2rMhcGoxnKSnGD5ROx8PGvQ09SxDOolh69stXAd46ExsUUdY50Gy9iJHcfNXInI08C+UuqjWX6oKtCcwvXfDNxTSn1pyLbvAj4xKJyG8P8inUQ8FBF5D/AegJs3b560nAvPee79FT+sME6jWw8TxMOEdDEq9l7b51LVZaVsUy/bbDV9UIJjSeoabhl9b7LBTKe59iak8fKEbM6Sa3NltcR+p0e1ZBEm6eTWRi/EEEXZNjAMAyHBEqHphwRRRDuIEVOwEsWKa1MvuTSDdFLwfhLS8VUqmOLUjj3xPcv+7kdAlKoJFUtRcw08x6JWdlnxLJ7aKNMKEzqZ1lhyLWqZd+NBN8ZzTGxTuHxtBdc2aPohlyoepoJqKY0w0QtDbMuk64ckccJ20ydMety8VME20snNubNJrgHnOb4ur3pHessNywkm2fphYZLyd8WxDF7a7RDGirJj9s8fRAl39ru0gpAEq/9sz9P3Mopl7MSO42b+faQN/CXgSVKT288A33bMcR8Frg7Z9KNKqX+a/f8u4ENDjn098H7gbcdc40dJO5MfHLWPUuoDwAcAbt26Ncoici5Yxt7fJOMKuVnm2mqpH3cPDifMK86vCTPX8WJOKdcyAJPtpt83A3mO1Tc75T31vPEczHm0WfOol510nCmMuFz3uLListVKG91nnm8QpF4UlG2L6iWHspfOUXpsvcx+N6bsWLy018IwDDp+TJIkKBTtrt+PiJGo1OchmeK9jiLoSoJtKcI4JoqEdhhjYrBWdVhxHTzH4FLFA0MRxQntbkSgYlzL5Bseu0TLj/qCPIrTffbbEYaRakb1ss2nX9pnq+Wz1/Rp9kIe26j0hVBRgCSofuDdkak64sPPr5r1/vMwSYOCLX9X0sC06qE5bhfZEWIZO7Hjmvi+EfhdAKXUl0Tk8nEHKaXectR2EbFITYdvGlj/CPBh4N1KqWePOP57gT8OfFtmbtQsGZOOK+TzUxrdgBteOtA+6vh8fZQ1+HmA0pf3OkQxbLe6PLKWztPJG6li4+knCVGsUpOf9SDnUW6Sqlds6spmpWTTCVKNyg9iXtjpYFnpHJrrdY+Sa7LmuTy/26YbJIRJwmuvrhBEikQlXK0JJcei40d0w5hWkBBGqWAyeGDamwYRYJlgkab7COKE2/tdNislMExqJYuKY3O/3SOKU1Ndww/ZKNusr7h9J4IgSvCDBNc2WS3ZNOpptIluGLHXDthth/TChOd3O6xULV55tdbvMOSRPPw4zuI6xn1z6zAGc4KVHQs/SkaaqYohmOwsAnpRWziNI8SyORgMY9k6seMIKF8pFeQ5CjPBMo1v5i3A55VSL+UrRKROGpT2vUqpj406UET+GPCXgT+ilOpMoSyaOTCJTVyEw9HGM+2oE8Sslu1DMQEHz7uWjSdtHfSIYqhXbPY6PraZhhUqNlJ5Y7jXTA1rBnBlxeubRYIw4TMvN0hQGAhvuFHveyt88f4BnTCi24qIkoTXXV+h6lqUXIPVik3Vtam6Jq5r4tlCMxDWax5tP2S7ExImiiRJww65FiQJGEaqTXXiobdlIkxADGjHCX4cEaEQw6BecRBDZeGeTMxIiMxUa9zvBDQ6Pi0/xjFMbl4q0235GCbc2e+wUfEwgN1WDyXCQcenG4Z0wtSzo9FOvQQrjtU37/lhwouNLldX0snEFdcdWebBXj/QD580ykyVNsJO3yQ86NF5Ei1iGR0MzgPjCKh/KyJ/ldRh4a3A9wO/NoVrv5OHzXs/CDwFvE9E3pete5tS6r6I/BzwM0qpp4G/DbjAb2aC83eUUn9+CmXSnCGT2MSL+0ZxQrsXsdXqsdPy2Wn5XF31DvWSi+fN4+u1gzQNxl6nx1rF4cnL1UOpIfLGq9EJ4Ea63PHj/kB77sWXoFiveuy0eoRJwuUVD4ADP+TRSz322hFf3W6y1fL58v0WX3NzDUHYqDo8tlHBDyIe36jS9NP5Uy/utfGDtKyGZGa9GFRuXxzSHRTA5uHxqWIydVXYLzHTc7uWRTeIOGiHdGKFqsDn7+7zDY9e4sZqie12SKPdY7cdUi3ZrJYc4iihFylu73XoRjFV1yRKIEkS2kHEE5tV2n7aUbhUdokUNNohq2Wb115b5XJmFg3j1H2+5Yd0g5CWb9L0U0GYR34fpUEX148rYEZpCyfRIpbRweA8MI6Aei9pNIlPA38O+OfAz532wkqp7x2y7ieAnxix/58t/P/Uaa+vmT+T9GaLpputZo+tls9OK+CRSxW6QSpEjuol77UDQPi6m3Wa3Yib6+W+YBm8Ttmx6AYxz2+nk3TLWSy5fJuBsNNKs8iWnfQTqrgWZdvEswxMAy6VbO7s+9xudNnrRLz6SoVvemKdR9fK/M5XtvBj6AUJjV6Iaxh4jkmrq7ANsD3o+ul8LNMEzwInhk6YCiQh1YaGjU8FpEKq5KZjTnEMyky1MseCMI4xRRDDgDhipeTQ6ES83EgNESXXJAaiBFq9kG6Y0PFDTMMkDGMSFPf2EkquRZQoTCumkQUNDqKEetnmba+7zgt7HTzLpOKa/fGjIEroBjF+kuA5Nt0gwjatobEQj3tvBiffztr0towOBueBcbz4EuB/yX4azVQZ1ZsdNeE2T9VQcgx2Wj26QUzZOZxivXjeIErYawfcbnTYaQWA4upqqd9oDrvu7UaHrYNUO3t8o4KCfhw6gDfcqBMmaRr0qvfgupZhUHFtVqNUu9jrtHFMQYyEZjY+lWpMXfww5vZBmzBOJ9C2/JA48/YT9cBJwogBA4I49cyzAMeEtYpBEqeTetvBA6+9tBxQc216VkQcpRNwgzhho+rRi2JClbDf6dGNFC/ttlmvlljxbHbaPgc7MZFK6PRCHluvIsAjlyq0eiGKhCABxxL2OgH1ksuqZ1NyzH4HIX9e9YpzyCz33HYrS+Oh2Ky4XF1JA/BahvHQONFR78Cw53XS0EaTMKozdR7GpRaZkQJKRD7NEWNNSqmvnUmJNBeeoxqdSSKq54E/U5fpkEfXy+x3woeEWX7NfFwrz/UUJIpeqOgG8UMC7vGN6sAk4IiDXsjmqkc7iHndtRp3DrpELdjvxJTtgGfvNolF+NzLB5iGwZ1GD8dO01JEkWKj5HDnoEccp/UuOeCaUPFcgtBnp51O3DVtqJUcur6iG/nYZjqpt0gvinBEMDwDpQTLhChKaHZ8DDHwY4XnGjiGYNmKIFZUSha2ZRBGCTXbYmPV4aAdogDXtogUqDjEVyYlyyRMYkp2Ovdr8BkUOx577YC7+z1KTpqn6jXXVvomUzicw2ucCdlFThPaaFKGaW56XGq2HKVB/fHsb54HKo+596eZnmORRvMQRzU645oFi4E/bdMkimO2mz6NdoBlCkGU9JPpFTMEd4OIbhDTC2JWXJP1qsu1eonbjS6gsEyDMJuHFcaHU8yHScJ+O8AURcm1eeMja8TxLlGSYIiw7wckiUGsFK4FppVGsVACFcei5Bnsdn2csonvR7iWhWUZBFGEH4FlgWuZlG3YqLqoCsQ7ChXHWH6czjMS8CxBxCBRMbaycBwhiE1AcCwb00zodtPxKNOyuFby0utbBpKkguzqqsOVssuVqospFmVb6EYJnSDiXqPHI2seZcdirZKOU+Ua5uj0KyprNBS2+SBFfD7vrDhn7bgJ2YfOmnkFhnFyaBJw/lxmOXakx6Vmz0gBpZR6HkBE3qqUemNh038rIp8gHZvSnHPmYcIYNeG2WI6jBFOuCZmm9KNqr5QcokQhprDXDmj2HsznOeiG7PdCntisstMKqHk2l6ou1+upKTCIEvZaPrcPeiSJ4rXXVh9KeV9xLVZLDs9vt3Ht1BHBsw0eXa+QxGm8um6UsOJavO7aKn4UcXuvQ6xMLBGevFwljBTBWowhJjv42GaaGypKElzXxFGCAI5lcqns4tomUSI0eyEYAVUnHQNrxwlxkLDXi/GjiBXDwTAMDEvhKANTLGpezCP1MrtNn4Mwxm76XMHhxlqJ/V7EetXh0Y0qj12qsN3yCeKEnbbP5opLJ4ipV1zsLMvu7UaHu/s9RmmXFdfi6mop1UxLdl84PbfdOnTcZs0ba0J28VlvNXv9fF31qvOQ1jXLsaNlHpdaFtPkOE4SIiLfnLt9i8gfYrxYlZolZ14mjGGuxeOUo1jeKBvXyCNsb9Rc9joBYRTzqRcbrHgWW/s+ppmmW7+3380cHhQbNfdQeKMwTthc9VAC3SD1IhzsOVdci1XPohfEuJbghwZXV8u0Q4UfxTiWwcZq6pL91GaN/a5PoxuxVrK51/b55lduYonBRz9/l51WwEbVoRtFxDFYkmpdhim4poltCwe9iMu2hWdC4pgY4vCK9Spf8+gqv/dCI02hERvEcYJlmqyX04gXQRz1BYFjGnwuUVRci8RQ7PZCSm0TwxRWvXTeU8WzqFfSyO67bZ+DbpoiRBSslFLhv9NKtVLLNAljNXRu0vV6+VB8vPS+qVSDBDpBTCe7t/mcp826dygU0iDFSbxtPw1ZNey5zGpy6jJOfIXlMk2Omw/qfxOR1Wy5QRpeSHPOmacJo6gl5SnDjytHXl7bNOgGMfWyjZ1N7swbRdM0WSnZvObaCnf2utzZ99msApkbeD6hNIqTvpYUxUnfjLRW9ah6Vr93nyc6DMKAL99vEiSK/U7E9Usl6hWHb6g47LcDOkFC1TPpRjFPbFb53O2IThDRi2LiJM1Ia5vCay5XaK5V8MOY57aa3Gv6VCwTI0nYqLnsdkJCFdMOYuoli1h5vLTXo2QLiUp4pO5x0K6y2wmRRFFxHVY9myeuVlkpWTimyUbFZq8XYQIH3ZAoEfwgZrVmc6Ve5qVGm3sHPWwT1OVa//56tkGj20MpodEJafohL2x32e8GmKZwdcWlXrKHOjvkmk3bj3Cs8qEUKFEcg0iaVh5YKztjZWq2s2ez3Uy1zbKTJ1U8rNHMcnLqsk18heUyTY7jxfdx4OtyAaWU2j/mEM05YRYmjJOYFsYtR95gvbTbJooV9w66fdf0G2vlvimuF0Z0gwjHEq6velQ9m2rJYqPq9c1LQZSw1wmouFY6LwpBSMeoNipOP917nujw/kEPP1S8Yr3MdifAtUy8zLvttik0t9p0w7QRBvATxRuurdDoBniWyU4n4Pn7LZp+jGUKN9Y81qoOQZjg2MJmzaPmWST32yRJQr3iYpkmhgiGkRCE0Aoi7jUDImJuXvIwTbha8/BckxXP4lLV5bVXVnluu8Vm2U3DG3k2fpxQsS1c2+JSNQ139Nh6pa9BbjV7NLoB//EruwSxQiUJl2sO95sBq55D2TO5XPN4YrPKZi113S+OA7X9iE4QH3Inr7gWj29U+5HHW37U14RO2ujPU6NZFpMZLJdp8igvvj8N/IPMzfwhwSQiTwLXlFL/52yLqJkX0/7gT2paGLcceUikMFbYlnCn0T0U663iWlxe8fhGZ4NOEGEbBlutXrp/lvcpbxxtM+lHLIhjRa1kcWXV46AbsllLtai2H2FlwWQPuqlGoZSglKLimmxlwU23mz4l16Ts2lwq2yBQskwu18u0s5BBdxpdbh/0KNkWniO8arNGq5vQ8dN5RiKCwmC9atMOYh67VOHyqoMfKvygzYEfIOLw5XtNwiSh5th4VohnG/SCmJ0k5O6+z347ZKWcZtSNgUfX0zxZVc/qT2ou2SaC9OMQKlLnjFrJwhBJkxeGMa5pYJhCL4yxTIPVUuq6X3zGmzWPrWY+obrH1dXSgGbjYJsG+52QRifAMY2xG8wwThND5mGo8md81gJi3Pd6UYTYMpkmj9Kg1oFnROTjpCnetwCPNNLDHyHN3aQdJc450zRhnMa0UJzXVOydD1JxrXTSaJxgIASZd1ex0XMsAxErG2OpHgouW7xWMfDsVrNHGCeH5lzlHmQHcYhnm/yBJ9d5ebfLQS+kZFscxCGuZbBR9Wj6aaRyA4fVksOVVY9GN+By3cMyDJ7batHx07EqEgvLMil7FjfXqyAdbqyVaHRjbm5U2W+HvPpKhWrJ5X7DZ9WzCaIY1zGJkgTHNIlRXF31KJsWZUfRDtJJtvt+yI01D8cSSBQi4JomrV6EKUIYK7YOumyulrBNs9/b9qMYP1BUPZP1isOtx9f7oaa2D3rUvDSlRr3sHHrGnSAV4k9crrLfCQ9NqIaHHR02697Y78SiaALjvNeLNu6zLKbJo7z4/icR+dukuZq+GfhaoAt8DvgzSqkXzqaImvPCaRuUcT7yomC5US8/NMg+7By2aRxyrqiXnb4rdC6MinmhckGZN6zEsF51sAyDmmvTzcarbEMo2xYv7nbo+AmGkUZir3oW1+tlgkixWQ3YbUckGwmGIbiWyVrJZHPFpbprpibFqkM3jGh2YzzboOJYVEsOrV5AyTVxrTS303rZ5lVXV7ANk2YvZLVqEwTpJOBmEGMZQqsb8eJuB9sw2aha3D/wcZ00UeFBN+DapTK73ZDHLlexslQZuWn0Dz4luFYahfx6vdwP/bRSsqmXU0eKoqODQH9cKIiSoROq++OGhkGz69Pxo0OTn49iUTSBwTBcw1KBLNO4zyJx5JuglIpJs9f+5tkUR3OeOW2DMu5HflTvcNg5gL5zxXNbTZ7baeNZBldXS1wvCLli41r0ILPNNCBtGCfstnwsC1q9mBv1EraVmhjLrtl3tCC7fsU12e2knmyvqaxQsQ0SJcRKYRoGr7u2wgtOmyRRPL/TZaOSuryXHDN1cReTsmtx7VKJzchhs1ri6mqJG/US3SB19LAsYaPi8uX7Tb5yv02oEsJIgR3TCk1sSWj3YKNmYhoGkrmyd/yY1ZLR1yzLjkW9bD+I75dHjMjGjbaaPe7sdbm2VsIxjQFHh6Mz3Hb9mM/ebqCAZjek7G48JKRGmcf6c9AK2u9Zk7/X+X0YlgpkUbS9ZWO8ropGMyWKwmNSm/w4H/lx5xx1DiH3aINqycIy0ggPL+yk8fiGRbQoepDZpsGdRpd2EFH1HMq2yfW1Uv+6zW5IL4q5s9+lXnEOjZdZhvCFuweUHYeSaxJG6cRe2zS5uV7DMYVurNjvBOzth7Q6IV+8f4BKhJINaxWXimNyqeLw+uur2KbBV+6nwWo9y6Rkpxrby/tdJIQgiggD2AkDVkom11ZLXKuX8GyTtbLD1RWP62vpWFHufTeoWRYnTtdLDi/tdUiU0A0Pz1XLQ0el43zBQ3Ok0rBINqtlm8srJQ66qRZWFFB5RBDTlIeSGy6K6SydFJyG4TrNBHPNYbSA0syFkzQsx33kk5oAi+fIe8Cg2GkFRES4JQczc54YTCHf9qNUE8k8CfqJ8LKJwY5n9zWujYpLlM2l8hyzf5403I/w3HaLRifENiVNVd726QUh9w581souK2WLyxWHmm1TL7u4psFaxaHpBwSRouLAKzYqrJVd9jshz++0+P07B+lYm4J7zR5397t0/IgoFtYrDquVtHxVx+JNj11KQzuFyaEYg7l7v20avLSbxg0cNNPlGWr9OMGPIva7ULEf5Hdq+1E/zNFBLw0zNTgOtVpyKNsWB93gUADe/Px5RJCy8yCwbFFzWhTT2XEdqEnGfRbFoWLeaAGlmQsnbVgmNd+NawLMPcoqrsX1+oPwO1vN3kMRLZ7bbnHQjdjvBbzqysoD85Jp9CcGF3NMPXm5hm0afQ2g6MW2WfPY74Y0/ZhuGNMLY1Y8i91uwBe3Wjy6ltCLXf7wUxt0gph7zS6NdkgUx8RROqZ1Z6+HJQYtP0Ql0PJjgjDi5b0gHR9b9ej4IY9uVAiihCtVB9O0uLGWanBhkgrcRifop77I5yrlmiXIQ3m3gijppx/JPfjqJWdIEsB08u5O02er4tPOXMpzbazqWbzxsUvsd4OHGuS+4B+RAXdRpkLA9LSkRdEKF4FxUr7/DeB/UEo1suU14IeVUn9txmXTnGNm0bBM45y5oHqw/HDajrv7XQBu7/ZY9WzWq27foSIXjPlcKscyqHoWT12pDW24Kq5FxbGolyw2Vzwsge1mgGlEuKaBa1n4QUSYKJ68XOPmeoUwSuiEEc9utah5Np1ezO1GO/WkawWEUULZcVirpOkz6p7Lp1/o8txWG6WgaluUnYhGK8CyhGe3WulkYRFedW2l30DnURjy7LRBlPSdAJSKuN3o0PZj9jo+axWHNc/h0bV0zC6vex7m6KAbsl5zWSnZfHW7xXPbbTzbOBQaKYgS/Ew7zRtl2xyeAbcoRBZhKsTh9+d0ZVgkrXDejKNB/adKqb+aLyil9kTk2wEtoDSHmKTnOQub/KzOOXieMFG0ejGGqeiEMTcLJq+tZo+7+2nK9PWqx7XVEmGSlqeYHLF4/pvrFQBMM3VQKNsWkHCpbNNop152L+20sAxhs+qBpPusZ6bDiIQYhW1brJUSOpHiUs1lr2Ox5qVRI157Y4Wya2MAj21UcG2TJElTwJtZyot2L2Q/y+1UDLq6VnH6jWXuBNDshmy3fWqeTRQremFCzbX5wt0DLq96h8aKHt+o9o/tBmnkjKpn9QPv5sJ8WHbkYc90mBAZFqH+JCyCcNAOFQ8Y56maIuIqpXwAESmRZrPVaPqcdEzpNB//qJxRo845Dbt+xbWyMaUul8plgjhhq/nAbBXGCssUdjsxB702z7y4m4b1MQxef2OVV15ZAQ6nmMg1rH5k77rB9bUStZLD3b0uhiVUSw6NTsjzO20Mgd1WyFNXakSJYrPqcL+p2Gv7lG2L11ytUK/YtP0I1zLYbqWJE7thQhjFRFGMZRisVx1Qwk67Bwg31spcrx92jig+y6ITwEE3JE5UOgQnULLT5IQ7bYVjPfD+y59Hbj49NM4XR9RLdn+e2U6WIbmYHXnwmeZmxSBz2pi2EFkE4aAdKh4wjoD6IPBbIvL3suX/HPiF2RVJs4ycdc9zUoF4UtPNsCjqT2bx6dpBjMRRv9cPYJtCqxfhh6k3X5II1XLaCHeDB2M9w8pRXF9xLeolB1Pgzn6PRtun5NrESlEve2y3fESg44dEKk1QuFKyeWKjSjuIUCpN4bFZ8wgjRbscc9k00nlWFZuyY1HOxnM2ay5hkrBacg45Rww+y2Jqi7JjcqNeAoSykzpF9CdHRw9PjoaHx/mKAXkt0+AVm7V+pI5hzyZ/hkGccD+L0jFJ5ImjnmuxjIsgHJZlIu2sGScW3/tF5FPAW7JVf10p9ZHZFkuzbJx1z3NSgXgSATpKqDlWmscpUdALo3707XQcqkq97HC70UUlsN8JaPkhtpFmAQaGlqNYvkYn4KDRpdkLCRPFY+sVTEOwxODL95uI6mKQurkrEW6sps4Oj61XWMvmShUb2NyEmKDYawX4kcndRovLq14/nJGVaU1F54jis8xdvRPSVL/5WFAx4nwYJ1yuen1PwFH3N11vHZocnTNsMu/gM8yzIa+VnaFm0+MoCro4Vtxcr4w9Ofioc85bqJ1Hxn0qnwMipdRHRaQsIjWlVHOWBdMsF7PueQ42AJMKxJMI0FFCLe/xX19zKLvmQw3l5RUvM/clPL5RoROmprtcYxlWjuK8qjBOJ+Ou11K38XrZxo+SfsPtOgYbldTK/tmXG2y3fAxkIMbdg/pVvdSrcLuVah0V12Kn5eOY6Vwv4KF4dsVnGUQJn797wE7Lp+yYrFUcRIb38u8fHI5aPq7H5Vp2/456d4r3zhkxpjcOYZz0syx3gvTZPHWlNnJ866QdGc3pGceL7/uA9wCXgCeBG8DPAN8226Jplo1RZonT9i5HNQCTCMSTCNBBYZI7Dqhs+lPeUJYd66FIBvm9qLhQiVJtwY+CfgDVo/Ic2ZmzRB4eaLWUJuJrZNpY2Uvj+FVci81Vr58AUOShU/Xv31azRxAn7LUCgH6cQttMDxqsY+4Cns9D2mn5vLjTwcsiYjy+Xn3oOpNoqYP3dpwgr9PqBNmmQRyrNLqFa2OaMlSTHVfTXgTHivPKOBrUDwDfCPwugFLqSyJyeaal0pwbptG7HNUATGqnP8lEyVyY5AP5xSjdw9YPq99g+UXoa0PFDMDFyNyDGoVjlbO0H+mxX7mfupi3emHfFXuUVjjMNPZ4FsW8aJ4bVpcHE5BNPNfkcs3lWr00VBgO0w5PMtZzVIdmGmMzg56TxXt3Ek17ERwrzivjCChfKRVI9kaKSJp2VKMZg2n0Ls+6ARjMzJvnOBomZEY5FLR6UT+lB0I/gGoxmChw6DphnNANYoqpP3IcK9XUtps+202fPPPvMBPjIOOYxhzLGFqXXPitVRx6Qcy1eumhybLFc0ySCXmYsDkrc9mouWnjammDQnQRHCvOI+MIqH8rIn8VKInIW4HvB35ttsXSnBemNXl2kkZDKY40oY06Lt8/H6MAeH6rw0E3ouSkvexGJyCOFUqNrl+rF/HM87v4ScJeM+C1N1ZwTINq5gCRBxMtpqZodALyvFRhpDJ37AdCqhg9XWJYr3ipMBxjLGbc+5fXpVjH4rFFresoB4h827iZkIucpblslDZ2nJZ2lPOMZrqMI6DeS5r2/dPAnwP+OfBzsyyU5vwwrd7luI1G7oI8OFn0uOOKjY1ScH+/x0Ev4rntFl/jrNIJ4MnNGvvdENOUvsfbsPo9CP9js9PsYSBYmWAuBhOFB2NZcZxGNS87Fl/daqaTcR2zMAfpgZmubY7nVHDU/Rs1h2yz5vUDsxbrOKlp1DYPO7J0g8y9XnGkx9wymMv0mNPZMY6beSIivwr8qlJqa/ZF0iwS03CfnbR3eZJr5o2GY6X5ioqZdI9zPw/iB7Hlwjh1Nri86rFSiml0/VSKKAEBL3ODPmosrOxYGAj73RCDtDzF3EhFx4A8xXw+/jMs7l1+T4p5loY5FYx7344yo8mIOgJ9s2UeTHac895YS8fOnttu9b0N3/jYpZFCapiZ8KgElfNgGYToeeGolO8C/Bjwg4CRrYuB/1kp9eNnUzzNPJnW5NazuGbeaATR6Ey6w8i1pQSFgXCjnmkMpgEOXF8pseI5h7zpjmuY8uCn+RiUYx92dhimuaR/y5lQiLnb6FFy0smx+f0ADuVZKt5nOHq8p7jvURrAqMY3N1sGcUySwK3H11mrPohZCMM1i1wARzFcqjhst322Ww80s2HkAn9R3bf1mNPZcZQG9V+TZtL9T5RSXwUQkSeAvysi/7VS6qfPooCa+XHaya25g8E4LsSnuSYcn0l3FLm25GSREHINKvfSGzzXMAEzjKo3XMsoNr5FzaCvLSnYbQdEMVi9tBzF+zGq8R5MtR4WtMI8Dp6VCZ/Nmnekt92wxrcTRARxTKxgr+Pz7HaLr/Xqh+7BKC++Rieg0fHZbvZA4OqKx8t7nWMFziKb0vSY09lwlID6M8BblVLb+Qql1FdE5E8D/wrQAuqcc5rJrbZ5OIfQpJrQScwnJ2k0ci+1IE7oBTG3Gx1KjnVkPqnTNkyDwmWz5vXdu/daqbPEparDfjdMBSYP34/BxhsYKhxe3uvQCWJ2Wj1esVnrC+DjvO0GozmUHYskSYVTybbwLPMhgTFMs2j7ESXH4tYr1nl5r4ubjXONI3AuqilNR6V4wFECyi4Kpxyl1JaI2Ke5qIj8EvDqbLEONJRSX595Cf4k4AAB8CNKqX895Pi/DrwDSID7wPcqpW6fpkyahzmJKSNvVEblEJrFNU9D0TEgUbDTCnjisnNsmU/TiAwKl07wwNttu+nz0m6HrZaPaxiUH7XYrFkPXWvYRNfcnFYUDnm6+ZcbaZSKetkuaIPje9tVPYtbj6/z7HYLz0oDw45yNR+mVSngyorbv944AucimtIW1aw5L44SUMEJtx2LUup78v9F5KeA/WxxG3i7Uuq2iLwB+Ahp5IpB/qZS6v+dHf8XgfcBf/40ZdIMZ1KNIW9UijmEzkITOg25Y0CtZNPsBexn2W3z+UrDnBGGNSLjCq1B4ZI7TzQ6AfudgKuZyXGj5uHYwyckFxvvPHJ4Hv2heJ0oTri738MSiGM1NM7dOJpKECU4tsFrrq5M5MI/zOlhEoFz0Uxpi2zWnAdHCaivE5GDIesF8KZx8cwR47uBbwVQSj1T2PxZ0rlX/VQfOUqpYrkq6InDC0XaqDgP9egXlbyBDuOEq6sl6mWHRifgfrNH3Hg4mOiwRgSOdlIoMkwzyCNFhGsJ3TCh44e41mgnj+Kcr1GRLHLtMIwVJdvkpb12PzXI4H5HaSrTSOKX37dBIao5zEU1a45i5JuilDLP4PpvBu4ppb40ZNt3AZ8YFE45IvLfA+8m1b6+ZdQFROQ9pLEEuXnz5qkLrBmfZen9DjbQYZyw1wk46EV0/BB4EEwUhjcik/Z8B++NYxn9/Ea2meAYwkZteNq1osDoBTGmKSNzI1Vci7JjZkFhR5tcj3pWp+3Va7PV+FxEs+ZRzKz2IvJREfnMkN87Cru9C/jQkGNfD7yfdGLwUJRSP6qUepQ0X9UPHrHfB5RSt5RStzY3N09eIc25xrGMvrdhP5ioH1J2rH4w0eK+N9bKXFnxDqUmH9bzzb3o8nxRx5XhxlqZyzWPkmPS8qN08vHAsUWBYZpCHKuHnCPya+bnvF4vcXXVO5S2fZwywel79cXy5nPTNKMpvosXnZnp2kqptxy1PYvp953AmwbWPwJ8GHi3UurZMS71QdLoFj92wqJqNIdwrNHBRIv7HOXBBrDXDvrmt3HzDjnW4cy1wzSWwdh6m3XvUODX4WF4Hk7b3ugEY2c+Ps3kWW220pyUeRqD3wJ8Xin1Ur5CROrAbwDvVUp9bNSBIvLKglnwHcDnZ1lQzcVjVDDRoxico9QJYu7tdym5FkGU5l0qmgpHcVyDfpQZ6CiPvHGE37h1m8Rcp81WmpMyzzflnTxs3vtB4CngfSLyyex3GUBEfk5EbmX7/WRmLvw94G3AXzqzUmsuDCc1teQmrZWSTZQkNLvRUFPhUdcdNCGOW7bjhNtpTZHjmOuGnesszVaTmFU1i40odXEc4G7duqWefvrpeRdDc84pahndIEq96BxzrOC107r+UdrK4PZJtKLj9p23Q8S8r6+ZHBH5uFLq1rBt2t9To5kyp537M43rTzLBeBIvvePMdfOexzPv6y8bix61QgsojWYGDHOiOCtGRR0fpV1M6sRwlACct0PEvK+/TCyDtqkFlEZzjsijjufR2YupLYopSfY7YRZ81pmqE8O8HSLmff1lYhm0TS2gNDNnUjPCac0Oi262mJRJ6tMJIhIU61WPuwfdQ6kt8tBHL+12gDR7b+64MM1J1fOeoD3v6y8Ly6BtagG1hCxTAzypGeG0ZoezMluc1TOYtD55ssS7B132mgGtlehQaos89NFKySaMhwfEXab3a1LOc90mZRm0TS2gloxlsBsXmdSMcFqzw1mYLc7yGUxanzxZ4narR2sleii1RR76KE/jMdhrXrb3axLOS92mKWQXXdvUAmrJWAa7cZFJzQinNTuchdniLJ/BSepT9Swcq8zLe52HjhvVa84bvSB6kM/roPtgnOo8sGzfzjDOi5AdFy2gloxlsBsXmdSMcFqzw1mYLU7zDCbt/U5an3ySKjzICjx43GCveTALchgnvNTyATk0TrXsLNu3M4zzIGQnQQuoJWMZ7MaDTGpGOK3ZYdZmi5M+g5P2fsetTxAlPLfd4u5+FxCurno8vlE99tjBRi9NhzF5sslFZxm/nUHOg5CdBC2glpBFtxtfBE7yDGbd+02dHhRelrI+jNVY1yg2elGcUHWtEyebXHSW/ds5D0J2ErSA0swU7TWVkqe3iOJkZr1f2zSwTaEXRICwWrLHukbe6OWRzluZiXCt7MzUvKffjZOx7EJ2ErSA0syMizagO4rifYDZNfyOZfD4RpXNWprwepJrONbDkc5n2RDqd0MzDvqN0MwMnagupXgfrH6K99l8eo5lsFZxWKs4E1/jLMc35v1u6Ijn4zHv+6Q1KM3MuGgDuqNQKk3NHsbJ0OSHi8I0xzeOM9/N893Q2tt4LMJ90gJKMzMu2oDuMIIoYavZ66dm36x7C30fTqPd5UJJKfqZhEc1bPN8Ny6aq/ZJWYT7pAWUZqYs44DuNAfv84+8XnZSbUGmU8bjOGsHhGJvuxfEmKb06zyqYZvXu6E1+/FYhPukBZRGU2DaZo15fOTzMM0Ue9thnBDHaqI6n6VA1Zr9eCzCfdICSqMpMG2zxlEf+awa5XmYZoqC2DENNuvDo1gMYx4CdRk1+3kw7/ukBZRmqiz73JZZaDzDPvJZNcpnMd9qGKfpbS/CWIdmMdECSjM1FsHr57SclVljFo3yWc23GsVJe9uLMNahWUy0gNJMjfPSE56GJ9txwm0WjfLg/Z+3eWZcFmGsQ7OYaAGlmRonbXSX3SyYM4kGOYtGeZk1kWURppqzRQsozdQ4SaN7HsyCOZNqkNNulLUmojlvaAGlmSqTNrrnxSwIZ6fBHKVxak1Ec57QAkozV5bZLDXIWWgw50nj1GiOQwsozVw5a7PUrMe7Zq3BnCeNU6M5Di2gNHPnrMxSg9rHqJToi8x50jg1muPQAmoJOC9ebvOmqH00OgEv7LTxHHOpTGXaEUJzkdACasHRYw7To6h9xLHCNGUpTWXT0Dh1p0ezDMzlzRSRXxKRT2a/50Tkk9n6t4rIx0Xk09nfbz3mPD8sIkpENs6k4HNg3ondzhO59nFlxePmegXHNC6kqSzv9Nw76PHyXkcn7dMsLHPRoJRS35P/LyI/Bexni9vA25VSt0XkDcBHgBvDziEijwJvA16YcXHnih5zmC5F7cOxLqapTDtaaJaFuZr4RESA7wa+FUAp9Uxh82eBkoi4Sil/yOE/Dfxl4J/OvKBzRI85zI6LOmdId3o0y8K8x6DeDNxTSn1pyLbvAj4xTDiJyDuAl5VSn5KzygA3Ry5qQ6qZHoNjTrrTo1kGZiagROSjwNUhm35UKZVrPe8CPjTk2NcD7yc14Q1uKwN/ddi2EeV4D/AegJs3b45Vdo3mPDHK0UYLJs2iMzMBpZR6y1HbRcQCvhN408D6R4APA+9WSj075NAngVcAufb0CPAJEflGpdTdIeX4APABgFu3bqnB7RrNeUePOWmWlXma+N4CfF4p9VK+QkTqwG8A71VKfWzYQUqpTwOXC8c8B9xSSm3PtLQazZKix5w0y8o839R38rB57weBp4D3FdzQLwOIyM+JyK2zLqRGs+wU3ev1PDrNMiFKXRyr161bt9TTTz8972JoNBqNJkNEPq6UGqp86K6URqPRaBYSLaA0Go1Gs5DMex6URqM5I3T8veXkIj83LaA0mguADjq8nFz053ZxaqrRXGB00OHl5KI/Ny2gNJoLgJ4LtZxc9OemTXwazQVAx99bTi76c9MCSjN3LvIg8Fmi4+8tJxf5uWkBpZkrF30QeJHQHQXNoqEFlGau6ECmi4HuKGgWEf0GaubKRR8EXhQuureYZjHRGpRmrlz0QeBFQXcUNIuIFlCauXORB4EXBd1R0CwiWkBpNBpAdxQ0i4d+GzUajUazkGgBpdFoNJqFRAsojUaj0SwkWkBpNBqNZiHRThIajUZTQEfUWBy0gNJoNJoMHVFjsdB3XqPRnFuCKKHtRwTReJExdESNxUJrUBqN5lxyEm1IR9RYLLSA0mg055KTBCLWETUWCy2gNBrNueSk2pCOqLE4aAGl0WjOJVobWn60gNJoNOcWrQ0tN/rJaTQajWYh0QJKo9FoNAuJFlAajUajWUi0gNJoNBrNQqIFlEaj0WgWEi2gNBqNRrOQiFJq3mU4M0RkC3h+3uUYgw1ge96FmALnoR7noQ6g67FInIc6wPTq8ZhSanPYhgsloJYFEXlaKXVr3uU4LeehHuehDqDrsUichzrA2dRDm/g0Go1Gs5BoAaXRaDSahUQLqMXkA/MuwJQ4D/U4D3UAXY9F4jzUAc6gHnoMSqPRaDQLidagNBqNRrOQaAGl0Wg0moVEC6gZICKPisi/EZHfF5HPishfGtj+wyKiRGQjW/6jIrIvIp/Mfu8bcd7/o7DPbRH51UmOn2UdCuX4ZLb/vx1x3leIyO+KyJdF5JdExMnWu9nyl7Ptj5+2DjOuxwdF5Asi8hkR+d9ExC4cuyzP4udF5KuFsn59tl5E5G9lz+L3ROQbTluHGdfjzL6Lk9RDRH6kUIbPiEgsIpeGnPfMvo0Z1mG634VSSv+m/AOuAd+Q/V8Dvgi8Llt+FPgI6YThjWzdHwV+fcJr/DLw7pMeP4M61IHfB25my5dHnPcfAe/M/v8Z4C9k/38/8DPZ/+8EfmnB6/HtgGS/DxXqsUzP4ueB/2xE3f5FVrdvAn53ketxlt/FSeoxcOzbgX89729jhnWY6nehNagZoJS6o5T6RPZ/E/gccCPb/NPAXwZO7J0iIivAtwK/erqSjuYEdfiTwK8opV7Ijrk/pNySlfufZKt+AfiO7P93ZMtk278t23/h6pGt/+cqA/iPwCOnLesoZlWHI3gH8Pez6v0OUBeRa6epQ1aOmdbjLL6LrByn+b7fRdpwH+Ksv41Z1CE711S/Cy2gZkymjr8R+F0ReQfwslLqU0N2/YMi8ikR+Rci8vpjTvsdwG8ppQ5OePxEjFmHVwFrIvLbIvJxEXn3kFOtAw2lVJQtv8SDj+IG8CJAtn0/238R61E8pw38GeBfFlYvw7PI+e8zM95Pi4ibres/i4zic5oKs3gWnPF3ARN934hIGfhjpFreIHP7NqZYh+J+0/kuTqsq6t+RanQV+DjwnUAZ+F1gNdv2HA9MGStAVT1Qkb90zHn/BfBdheWJjp9RHf428DtAhTRG15eAVw2cawP4cmH5UeAz2f+fAR4pbHuWIeaFRajHwHn/F+B/XLZnke13jdQU45L20N+Xrf914A8X9vst4Nai1mMe38Uk9Sjs/z3Ar40411y+jWnWYRbfhdagZkTWg/hl4INKqV8BngReAXxKRJ4jVX0/ISJXlVIHSqkWpCoyYEthoHjgvBvANwK/ka+b5PhZ1YG0x/cRpVRbKbUN/Dvg6wZOuUNqLrKy5UeAl7P/Xyb9KMm2r2b7n5oZ1CM/748Bm8AP5euW6FmgUjOPUkr5wN8jfa+g8Cwyis9p4eqRnffMvosT1CPnnYwwjTGHb2MGdcjPO73vYpo9Cv3r9xYE+PsUehBD9nmOBz3FqzyYNP2NwAv58pDj/jzwCwPrxj5+hnV4LWlP2yLtiX0GeMOQY/4xhweCvz/7/wc4PBD8j+b0LMatx58F/j1QWuJnca1w/v8R+Mls+f/GYSeJ/7jIz+Isv4uT1CNbXgV2gcoRx5zZtzHDOkz1uzj1S6d/Qx/SHyYdYPw94JPZ79tHPXzgB4HPAp8iNWn8ocJ+/xy4Xlj+beCPDZxr5PFnVYds+UdIva4+A/xXw+oAPEE6ePrl7IN0s/VetvzlbPsT83gWE9QjIjW15OfMzWPL9Cz+NfDpbJ//nQcmGAH+Tla/TzMl896s6nGW38Up6vG9wD8ccq65fBszrMNUvwsd6kij0Wg0C4keg9JoNBrNQqIFlEaj0WgWEi2gNBqNRrOQaAGl0Wg0moVECyiNRqPRLCRaQGk0UyCL7pxHev41EalPePw1Efn1Edt+W0RuTaWgEyIiXyMiPz+Pa2s0WkBpNNOhq5T6eqXUG0gnM/7AhMf/EGl4mJlRiFIwNkqpTwOPiMjNGRRJozkSLaA0munzH8gCfYrIN4rIfxCRZ0Tk34vIq0cc811kgTVFpCQi/1BEPiciHwZK+U4i8rbsfJ8QkX8sItVs/beLyOezoKp/K9fGROS/E5FfFJGPAb8oIpsi8ssi8v/Lft+c7VeRNH/Pf8zK+o5C2X6NNIKBRnOmaAGl0UwRETGBbwP+Wbbq88CblVJvBN4H/I0hx7wC2FNpTDyAvwB0lFKvBX4MeFO23wbw14C3KKW+AXga+CER8YCfBf5TpdSbSOOgFXlddsy7gP8J+Gml1H9CKhR/LtvnR0lz/Hwj8C3A3xSRSrbtaeDNJ70nGs1JmVjl12g0QymJyCdJNafPAb+ZrV8FfkFEXkkaWsYecuw1YKuw/H8B/haAUur3ROT3svXfRCpsPiZpOiCHVFt7DfAVpdRXs/0+BLyncL5/ppTqZv+/BXidPEgntJJpYW8D/u8i8t9k6z3gZlaX+8D18W6DRjM9tIDSaKZDVyn19Vm+nI+QjkH9LeCvA/9GKfUnsrw7vz3sWFKBcBwC/GamCT1YmaVqP4J24X8D+CalVG/gHEKaquILQ473sjJqNGeKNvFpNFNEKdUB/iLww4XUCHnahO8dcdgXgccLy/+ONJssIvIG4Guz9b8DfLOIPJVtq4jIq4AvAE9kAhDSnD2j+FfAf5kvFITbR4D/MhNUiMgbC8e8ijRYq0ZzpmgBpdFMGaXUM6RRot8F/A/A/0dEnmGExUIp1QaezQUP8HeBqoh8Dvhx0oRyKKW2SIXchzKz338AXpOZ774f+Jci8nGgSZp1dRh/EbglaRbd3ydNUwGppmcDvycin82Wc76FQp4ljeas0NHMNZoFQET+BPAmpdRfO+HxVaVUK9OA/g5pxtKfnkK5XODfkmbYjY7bX6OZJlqD0mgWAKXUh0nz75yU78ucND5Lalb82SkUC1JHifdq4aSZB1qD0mg0Gs1CojUojUaj0SwkWkBpNBqNZiHRAkqj0Wg0C4kWUBqNRqNZSLSA0mg0Gs1C8v8H1hY6EwQGHDcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "cts.skyplot(cluster)" ] } ], "metadata": { "celltoolbar": "Edit Metadata", "kernelspec": { "display_name": "Python 3", "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.9.1" }, "toc": { "base_numbering": 1, "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": true, "toc_window_display": false } }, "nbformat": 4, "nbformat_minor": 4 }