{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Tutorial 3\n", "\n", "## Characterization of Stellar Activity Components" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This tutorial demonstrates how `ARVE` can be used to characterize stellar activity components (such as oscillations and granulation phenomena) by fitting the components with analytical functions in the velocity power spectral density (VPSD)." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We begin by importing the `ARVE` package and other useful packages." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "# Import packages\n", "import arve\n", "import matplotlib.pyplot as plt\n", "import pandas as pd" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For our analysis, we make use of some synthetic RVs computed in [Tutorial 0](https://arve.readthedocs.io/en/latest/tutorials/tutorial_0.html) and downloadable [here](https://github.com/almoulla/arve/tree/main/docs/tutorials/example_data)." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# Read RVs\n", "df = pd.read_csv(\"example_data/tutorial_3/vrad.csv\")\n", "time_val = df[\"time_val\"].values\n", "vrad_val = df[\"vrad_val\"].values\n", "vrad_err = df[\"vrad_err\"].values" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The RVs are fed to an `ARVE` object and plotted for visualization." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA90AAAC+CAYAAAAoRmzvAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAZyJJREFUeJztnQmcFMXZ/5/ZRVFBUTxAsoAoENGoUUTU10RE/x4xGo3x1iAxiIpEg6gourxeoKwHiqASwSMxEUw0vsQrRja+iseioJEXQzCIWTSs8UIWFRX6/3nK7ramp6q7uqd7unvm9/18CnZnZ3qqq+t6zipYlmURAAAAAAAAAAAAYqcu/ksCAAAAAAAAAACAgdANAAAAAAAAAAAkBIRuAAAAAAAAAAAgISB0AwAAAAAAAAAACQGhGwAAAAAAAAAASAgI3QAAAAAAAAAAQEJA6AYAAAAAAAAAABICQjcAAAAAAAAAAJAQHZK6cK2yYcMGevfdd2nzzTenQqGQdnUAAAAAAAAAACSAZVm0Zs0a6tGjB9XV6e3ZELpjhgXunj17pl0NAAAAAAAAAAAVoLW1lRoaGrR/h9AdM2zhdhp+iy22oCzyySefCMVAluvIDBo0iBYsWEBZJev1y0Md89AXs96GDOpY/f0w622YlzpmvX7oi7VRvzzUMQ99MettyKCO1d8PP7Hr6MiAOiB0x4zjUs4dI6udwyHrdayvr0f9aqCOWe+LeWhD1LH6+2Fe2jDrdcx6/RzQF6u7fnmpY9b7Yh7aEHWs/n7oEBRWjERqILOMGjWKskzW65eXOmadPLQh6lgb5KENs17HrNcvL2S9HbNev7zUMevkoQ1RR+BQsDj6G8TqYtClSxfq37+/0BxxR85aZ3bquHr16sxrjUB1g74IsgD6IcgK6IsgK6AvgizwSYb74bRp00RZv349/eMf/wisI9zLE4JjI7LWORw6duxIEyZMEP8DkCboiyALoB+CrIC+CLIC+iLIAh0z3A8dw6qjGAgClu4a0sgAAAAAAAAAfFi7lqhz569/bm8n6tQp7RqBKpD9ENMNAAAAAAAAAF5Y+GYhHIAygdANAAAAAAAAAAAkBIRuAAAAAAAAAAAgISB0AwCqC3YD47MSucAlDAAAAAAApAyEbgAAAKBagNIJAADKZiURNdv/RwJzMfAAoTshBg0aRLvssos4vw2AVMHEDwAAAABgxMx776XeRDSUiHoRUdOUKWlXCWQQlvFY1mOZzwQcGRYzODIMZI5aO/qi1u4XABn0fwAAiMzKlSupd+/etGHDhqLXm5qaaOzYseYXwlxcM3yCI8MAACXwAlADFm/XLeydd9KuCgDwNgEAgJywbNmyEoGbueSSS4RADkBUIHQDAKrWLaz3gAE0c+bMtKsEah1Z0IbQDQAAmaVfv35UYAWpBxbE33zzzVTqBKoDCN0AgKqBtdBnjR5NG6RFcuTIkdBOg5oCnh4AABCNhoYGuv6KK0per6+vp759+4a6FuZiIAOhG4AaoOwsnDl2C1u/fj2006BmgKdHDZJ2+ELa3w9AzFw0ahQ1SUISC9x33nmnEMhNwVwMvCCRWswgkRrIGjOnT6ezRo0S1l9eQGYQ0ZlVmtRDlQCFF8sVK1aEWiwBiA05mQ7T1ka03XaJfBX6f42SdsKmtL8fgLh57z2ibt2EoYJV9n2XLqWG/v2NP465uLb4BInU0gVHhoFMulsT0VlEtOCVV6ga4cVsxtSpVO+nnYZVpjap1HPXfI/rbfLuu4nVB54eAAAQH7xzGML/f+tboT6Hubg2mBbyyDAI3QmxYMECWrJkCY0aNSrtqoAaRjXx82/7DhlSta5OZw4bRitsAWfFkiV05plnpl0lUOPMtM96FWe+7rVXYmOPEwDV1dWVHYcIAABAgj05WDlqqCzFXBzA2uowfrCMx7Iey3wmQOgG+aBKBmilUU38zAbLquoEY8baaV2/qsb+Vo33lIN25QQ6I4jIiePiiK4RI0YkklinYautaMaGDf6eHgAAABLNhWPkdQdqDgjdAFQxzsSvGuhwdcoJEJbzc169Yw2ReP6ZZ1yBmyTB+4WXXqIkYL+OmvX0qOGxknaW5LS/H4DI84D3/WvX0sxu3b5JgmZ7K4UFXncg90I3+8/vsMMOtMkmm9DgwYOppaXF9/0PPvgg7bzzzuL9u+22Gz322GMlm5/GxkbafvvtadNNN6VDDjlEuOTK8PfxmX1yue666xK5P1BDk36FNog88b+oGOy17OpUEl8LgqlhgSYR1q375ueY2zNqHGJNUUX9Oe0syWl/PwBxsnLZMpH7Rs6FM5JfZ6WqnBTTAMzFILdC9+zZs2nMmDE0YcIEWrhwIe2xxx502GGH0XucZVDB888/TyeffLLQLi1atIiOOeYYURYvXuy+Z/LkyXTrrbfSHXfcQS+99BJ16tRJXPPzzz8vutZVV11F//73v90yevToxO8XgLgYZGcth6uTZ4M4cGAkDTbIJ5U6Ok+2+r2tCOHgkI/99t77m/d165Z7wQ9kJFnmhg0VDR2K/P0pJzYENY5Pv1i2ZInbnx3W21nMS64BQLUK3TfddJOIhRs+fLjIFseC8mabbUazZs1Svv+WW26hww8/nC666CIaMGAAXX311bTXXnvRbbfd5lq5p0yZQpdffjn96Ec/ot13353uu+8+evfdd+mPf/xj0bU233xz6t69u1tYOAcgT2jdTmtoU6LcIJoIYTXURpkggfYuUrZEdBc0+p7773e/p9fOO9MlV19d8p7rr7+enmxuLq7PvfcmVCNQzSz7299SzZKcyyzNtTaf19r9lkm/bt1KhKMCEXUK6ynnzRMDqNZDUnIjdH/xxRf0yiuvCPdv2VrAv7/wwgvKz/Dr8vsZtmI773/rrbdo1apVRe/hc9bYbd17TXYn33rrrWnPPfekpqYm+uqrrwLPbJPLOtmVEESiFgdo3LiuTltumRn390qi3CByGIqiX1Vjf6vGezJJiLeyUCg5Om9kjO3gtCvnLz1r7Fj3e1ixy8ULn98qv0/U5xe/iMU6WWTNr5Jxm6m4/YzNj/169CgNHeJzhQ866OuzhhOuI7I0g2qjYaONhGcgC9oOPIsPJqIzbCWpUJbuuSfNnD7d7KLszRTHWMzIvFPrISnr1q0rkfOqSuh+//33hfa0G3dcCf6dBWcV/Lrf+53/g675i1/8gh544AFqbm4WblMTJ06kiy++2Le+PXv2FAK8UyZNmhTyjkFmBmjYSY43OiHjflJPIqL6uYayuY/x9KuS/iZbIZPYzAdlwI5BWE5sDOVgE8BZOlTKljiSmcntuo/ie0ghEFmffRaPddDT9iXW/N//nmqNSoUQZIWGHj1KQ4dsBWvirF1LDT17ImN+1qmhNT7qmul9/+7SiRMO/DvvBIyUpdzOHtkCZCMkJg5YppNlPJb5qkroThOOIx8yZIhwPz/77LPpxhtvpKlTp/par1tbW2n16tVuufTSS6kmiWFDntcBajLpx2J5DGpjDoVga1t7e/TviLM+aWdzv+EGd4Mo4/QrPm+xpL/xwur9QIXuLQ5hOa64y7xayvvZnlFeTjrjjLKUD952NeE0Itp/333Ltw7ys5GUeytfe43OGjXK3Hqe4XFaqRCCvPbnItauLQ4devJJ8XslCZ0x39N3k6bkOcv9PWmPiCSowrGb5po5c9asovefcfXVQoFqQomy1Hk2aRpeku4/sgeNJqdW1YWkeGCZTpbxWOarKqF7m222EZuStra2otf5d46xVsGv+73f+T/MNRl2P2f38hUreJlRs8UWWxSVjh07GtwlqJYBylN80KQfi+XRu3lRbWZUi0ANa2DPPPVUsUG8SfE37lfPPfecur9R5dEKyyEXTO0Yev31SLHKeXMHY5sbK1u8C15kBZ49ppb17Olr2WbBmk+7kPmN/f/pxxxT9PqxRx8dzjroTf6zYkVp8h9+xqyBr4GNeclYCQghyIR7YxzC06efFocOeULqKkWULM2VUHoYP+coR0tFfU85lKuwqPKkjWEUzPzKHA4JGj++6P33/uUvxt8nQjl69Iip9jXA2rXUb+jQ3J+mwzKdV86rKqF74403poEDB9LTTz/tvsaDg3/fb7/9lJ/h1+X3M0899ZT7/j59+gjhWn4P++VzFnPdNZlXX31VbKa22267GO6sBgm5KOUtZoyTaxQdN8GT/s9//rWgZGsFU7PeexRMuT9CS9OXgjZzvEE8XnOM2gEHHKDsb51ScFvVCsshr6McQ7xZ2Gkno8/zPRfFICfdX2WruucIx7A41zlszBj6neLv5Sjw2ILuFaodbiai3919d0lMt+PW/mtPss7fP/ww3XDDDZHqIeqyww7q2F6qDXT5Gt785z9L3ms0/xqsU1mxlGfKpd7Q8qVVTCvOTI4qxGqf87vvJtdmcbpyq9oirMBtuzlnqo8kxdq1SkVo0Rxvt+nMzp1F/zvRICQoyHOpSNHkeebKds+C8kZHBb63oZZP07FyxAMPPGB17NjRuueee6wlS5ZYZ511lrXllltaq1atEn8//fTTrXHjxrnvnz9/vtWhQwfrhhtusN544w1rwoQJ1kYbbWS9/vrr7nuuu+46cY1HHnnE+tvf/mb96Ec/svr06WN99tln4u/PP/+8dfPNN1uvvvqq9c9//tP6zW9+Y2277bbWT3/6U2UdV69ezTss8T+wLKu9nbecXxf+WfdaAHdNm2bVfx1SY9XX11t33XWXVTFC1nfeH/4g6uktzfz5trav3zNvnvo9zc3R6+ZcX/7dW5Yvd39uJbLmEVlN48ZZdfb319XVfdO27e3ue1qXLk2svcq6pvf+iay7brpJfT8ObW3f3Dv3J6df8WenTSvtb0TWsBNO+Oaa/L4K3V9ra6u4B7mPcH1aI3y/957ce7D7pF+9uK20/TWB5x34DD31k/uz3Ff5nuXnNtn+v6g96+tFO4dC+t7JjY0l7eI8I66L9/nx98/WtGed/GyD2tQ71hcvFs+03nO/2msEPbckxnGCqMYKt0HLM8+UvNdo/g24/6K+5ddH/Yihje8aP94qSPcwWR4LixZFm8NNUcy/Yo0JgOulHYeqa0ZsI91zHjtiROl87jcPqu7XZFyWe03v31RtE9Qm7e3i/hJdv4JIei7xrAMlfcuZV+02VL0naqnzji1pjyW3O4/Rkfbcz/0/8Dn67HnKHtNh+lwcfVxz/VZ7X5zI3FRhTGW/XAndzNSpU61evXpZG2+8sbXPPvtYL774ovu3Aw880Bo2bFjR++fMmWP1799fvH/XXXe1Hn300aK/b9iwwbriiiusbt26CYH+4IMPtpZKHeCVV16xBg8ebHXp0sXaZJNNrAEDBlgTJ060Pv/889oSunUDMWjwqT4XZcDaE00qAzTkJKfcUDiTvn2/YjMe4+Zf3mQV/e6d2FtaShYDncBgLPhkSOgWi6lXSPW0a9F92ZtU0a88iiGnv7Uorik283EKKz59rEThFHUTKo8hr6LG5zPadnXqGscmQEKpaNCNDbt+RZtKu6/qFBZFipaoCjz5Gbe1iWvWqRQavOG9+OIioYjLZfYmTKucCyF0u20/d674eaR0bXeDrbpGlLk748hjhbxzl7zZ8xP6DO4/VB8N6DvltDF/X6FQKOlDzvPnvxXKVQwE3EvJWhMkdLe3myvxyhS6dYoY5Vw2f35JPUu+N0NCt+mcq9xrePYjiRPXXGK4D5CVj/Xc7z1tqOt/UcvY88//pi6LF2uFf3d8FgrF63gI4TcuZZ+2/yQtdC//RikR+RoZpGqF7qwDoTtgcOsWVANtbVobwLCTHFtFtROs00a29jmOzb9X4Bim02rbRQiRAYvIHAPh1W9jUDKhl/v8DBdb382czyasyHJseM2CZBk3qmvUPiYLy14NeUihW7uZ9QoA0mtarwBHqIxjEyDVL+gZej+jtGzU11uzZ8/WCrVlK/AUgpNSodHebk0+5xxlPY70UXgZPePly4uEfZ0QX2TlkZGerehXfvfoN+68a4BuLJn01XLnibY2a7yqDXju4met25g786+hQFyWp1KMQreuj+tKJMWuD15PEiOLsc+YjWzp9uk3rGSV562xujVv6tTia6r2NpUWur3X8lEy6pj36KPaebBahe6i+fjBB5Vztbf/8fzZSGSdc8opoYVuFuzlvsvXPyHoM6q5Xr4vRT8qS9knUbJmh52z29rEPvJG3k8+8US4Z7gcQncgW221VajStWtXa8WKFVYtAqG7+P0liwMLKXkQuv0sfD6TXEtLi3KCbdEs3mVt/g3dpOrt73cEJ93mXC7TQwo+8rPRTuhhn7nPdyhf17mWSc9L627oI3Sr+oC7UBcKxX3BK7ga3GPgQurZpLfaLmqum1qENlR6QiiEbl+vAI37dOSNfdCY0zx/nZDOXk6+rvlR5xKNd4mqcPuoLJFOe46WxmORhdygjpPHjTPeFM7RCCO+G/dKC91+3xf0fvu9/Cy0HgSPPabfmLOV0/v9PgJxnJbucrxEwgrd2jk8AtrQl6A5yVYUy8+pjhXTUliT6jn5KnEDBFdZIcZF5x1QpERVzeUmXj1ByqyoQrffGqdTrPGe5Jln9HuSKhe6S9Y4W+Brtdf8eunZm+yLAseVom/7fiak0K1VoIQY074ha4ZzttewNCzoucr30tJS+mwgdBfDE9Qtt9wiYqmDyt13321tuummIv65FnEanl3a2RX9tttus2pG6JY1WPbCo1scSgZd3BuymO43lNXNsqwbb7xR+f5ferWCcdxPe7tWa6/a4IdZQLSWbtlapBJQVUJYkEbX0+ah/qZwb9R6ENjvVS2ITpyt63LnWbj9lBUs3JUjdAdazaRr3jVpUlE9hBeFn7Ck2Jho4/s8ddeNX0eBI9omhk1ASZ3t+9Qq6xTPX6dokd2NhdKgsbE866vUx0wUS0FujLKV2o2/9hbFhoTvLeymkIVs9/Mm4S1+fVluE3nel3/228CpXo9B6Nbl0+Bn5LV0F23MNUK3n4AVOc+I9D1axYdhH9UJVJWwdOvmrduvusr3c6q+J56PpGSU53SlVTeM0O2xrOmEbmc+deuh6P8mXj2hw7IMFAby+sbhI8o5V9NXdHP02BT2U2V/l8E+1BGo3RALeY1bvrwktO7UmOK7m5qalH3bb+539xwaodt99n/5i35NDjmmtfsNaSz5zXtaw5Iid4Yyzv2CC6LlU8goLOOxrMcyX6xCd1uIRuncuXPNC921ZOku2rgYbDh5QjQedIoNr24jmrqlu71dCCRBk63IOxCHVYetaDEsFrrFoKmx0d/9UiF0+yaRCyt0G1iE5E1Qwe5brV4Llkcg8FNU1DnCiUI4PSUhoTvQaia5xBX8XNu8bagQaHRudWLh92xOdeO3yPJ9zTXlWbo149vXQ0YhdI/1UbQ0BXnbhFT4BSkjZKVEa8hxN90eey1ebwbP2Ihi4ZQFm0Blideabyp02/GMRfN0VKE76Nko5gud50/TNdeU9O8gLyyvl0dJKEnUPCMma4uhkDL73nt9n3nBm2tAFQ4TEdW85a5xp56qvfcgRbY8pzv3EEr56yN0t9rKb782c+dzrxLSwKtHFWNfkmgrTD4PTx/8iabOfjHaPEfr5psSD7ykiCvvh+469v5AFra9xZlX40qgRqqQpoDxKPfpgmLfIvffEuWNNG+pwrxisXS3tQWGLugMSzf7zVdynHuhdHz4Cuw5ATHdKVHNQrfKxcvP4qPaAIkkJqrJXxcXotvcVUjoFpPcpEl6i4bCzXKYwaTLk0zRhtYgDs67qdBtnMOW8bKblSfxjutObOLW5yQo0kzogd4Nmvi1kg2wT6Ispz/Jmm3vRj7IJZ+fsWpxrtMs6GJD7yC5r2k3GIpn6Ws1s9vFz2rqFZZ0Ao3uGocp2sokdIHbn7Pfh7b46ZRpftp8Rb9RKl2kNvdNIhRWsDNURsi/XxbD+BTeDB6BUJf9PKi4uQ107aKZ24UlXu7nnnnZa5FJTOjW9Rt+3f6uIiWL7EEguTW2BPQxI2WrqfVOMa/5CX9hTgQI2uT/3JtrIGaLJvdL3XezRUzVFpM172cFqcmcUxIqYCB0OwKZybVlobs1pFePThkmwjs07a+znvspNbRWa0+IX1B7ulbgMPuPCHkIIuf98DxX3XV4f2BigEjCSCGX26dMUb7uzEPcD27XPBPZABV0aklJ/hDdc9CMC11yVpMEk4EhlD5C9zxNuymVmjkDQndKVJXQLQs9Cu2XytJalJRKoZXTWRlFTJdqI5WG0O218tibATeRladtVGUukXUukXWm5n5PPO64YiuK5PpZcn3NpiKMK5PfIuQIxMKl3DvhepUCmsVXjmErmtBtwV3p3aDbnHuEbt0G2M/a52rxVUk7ArK3czkqRBsWCSse9zXdkWUl/ViVWVxhpQ5t6fY8N92Gl8vc224rEkq8ydN0GxY+fqfIsyDEGFMJ3FphxNv//WJKHQH5/vv9rxU0lyjqaSoYxFnEM/bUIcp15H7iVQ4OM5hbXAWTJ8Ge3N8ny8o1WSCXFVCqzaA3G3sYodv+rhJFgTzmx48PfG7O+PMViP3uIaj/sCXJY8ktx9Id5F5e5FqbwLrJz1T33XzUasn758/XziOsvDPJKq2zdOv6WZgYW9e93PM5py8FCSRhhW4/LyedG7CqzFXFkBtm6S4YxuFr1+oAwTuu/Ac6TwMWAk2VE0kWrsv0m29W/o0FbROFqdi3PPyw3htEN/eFFLqL9hvSsYKzZ8wwChdTxnTr5hhp/zK7gkkeq0bo5rjtP/3pT+7vF110kThOa7/99qvZ5GlVK3T7ZHr0m2CcBUbewDubGT/Lj9IClYZ7ucKC724kDYRuk0U+VDyl12VTEsZME6P5Fed5+U703jpIi6/KTdOZ0JsuvFAbQ91iKHSHSZQlF7842aaYF1xhebE3laGerUqZoRgDsoDj68rpvY5zfI9h0j2VEu0nAZ8Tbsvy90URRoKOsbO9D+RNtTY27aGHvtmoKTb5XsWgb101z8UraMbZl7R9zDAu1a9wvUVf8LN0+2zWHRd12fNE2y+4bx59tHnMsicbu86l26sAcV35VWNPUiqaPCdObBcoEKvqYtjPVZ5A8vXd9jHMgK0L51G1QRhByRQ/byvl+ejTp/vOI0Huv3JMtyxka63FIRXTQgBtbdWGbpWc4ODpnyr3cjd8R9Hufvk8dMfB6caaaoyZKgjnzJiRmNAdV6Z/3Zx00003hZ4H4y6O4eSUE09U/n2iLXCONLwW71sKPgYS5XGwqmcR4PnmzZ8iezoGJXVtsV3KXYWPbr6yXdYLFUzyWFVCNweLP/300+Ln559/3tpss82sO++80zrqqKOsY4891qp1qlHo1k12jZrsuY0aF7/pJoMuQOiWLapJ4bdQKxO5eDfDAfd4qMmEoxC6vZp8XvDjcJe6RKq71tItCd3yM/BzBQ7a6BZZvFVWLrleGlerIIu1K2hI1zJZ+MIU94zXKM82SOj2LIrO4qt8RpKLbVGeBfv6cZ9Nqhy3QUK3z7jRCju2Jt0kvMArYPDZ1cabF537s0eBKHttFGxFYl2FLd1Rn6UTX8rnymqfo22d1iYOLON+ZUuu16I9ecwY3/f7KWn4/7E//7l2PYrqjk9+sZMRhG6dkHyzJkY1aK3jv5nUf06c2YJl66PCHdU5Dk819oPWx2afOV2O/SwKK2FBwSOcutbiCCFYY8eO1brtny15hc3WJdiTFFK6eF0TK3AYodtvzJiskacdd1x2LN0+3nSqvQb3iUp7HikV/I2NsV1PF8r2wyCPwZBCt2pvxd9rFE7p6VtuX/fMV37eLSTfs+yxl0MSE7o5M/nbb78tfr744out008/Xfy8ePFia5tttrFqnWoUuv0muzhjZJQWKI8bo+vuJbujB1mowrjeGAgnYsHQHBOkUyxMtDdWLaYxqx6hu8TyPnlyrBa2ve26KZN0+Gx0tUKmoWAgC0mqrNrus5eTunnqFZS8ruAkhksgpsuty+LF2nqUZE6WvQa8QqpKGDV4zm7m0bY2beIwkyR/kcZtoSC03SI7f9hNQMCZps6Y0XkPeM/g9fbXebNmKa+rDXfwjj2DTSt/r+ps6DjL8cccU2zVjWjp5jLnvvu0n22S+qfuHPFyBVhWQHnPd+Z4cV2dSmKcNX2C+3pQ/GqU+k5QZRlWxaEbeEr45byQP2+U1dy2IJncQ6GcbMF+oU4az66S47ek6/h5GcnJCJVjlAV3w1hnx7su7HP38+Krs8e6an6V5w9vDHmJgOR51qr42nKUpMLbR3pOJmO2KDeJoq8VKYFCCN3iHuWQs6DkXz4Cfsn+hNu+rc0oj07SJel4cV1/LJqbvONTc+xdkPK20VYslSiVZIWbTjnmCambfeutRvfSxCeLmMw/tSZ0b7vtttbChQvFz9/97net++67T/z85ptvWp06dbJqnVwK3bpOLR9V5NV+2QtNnIJMkRuWwYbXfb9qIPok5fLdMBm6ZBW5sPq4/pbEdCnasyjZj1NfWcP7l78oN5hJTOK8eJXEFgdMssq4uBBCXrPm2nLmXaeNVPXqn8JiJx834ngC6OJATzzxxOK+tnhxsSXaR+g23Xg51nxVzDaP1yQUDk7Z+zvfKelDJkK3szENui/fbMe2VdbtF96MwwpLt9O3grwFQvX5hPsZjymvRWHobrtFutb0668PVHxO9nke5cS062Iw/eYzIQxI/SbJvqx6rlrljN/a6envLZJSSidkOZ/3tQzK48fQglTyDMJ6innvUXG0oPJ7vJYrT9iW3zji9wz1vO4IaiaxznJYQSGiYipUe8rtY89JJWu2LtmUPIdJiQijKtbqJM8r2UMnMLRIPi7N89y14XYmgpBkVTVK/uU5htbrkVR0Hd6vKcJK0ipx7oXD3JOslNIqS1Tro8/3iGtOmqQV2A/1GxOSJ8P0q64KrP+xfn2p1oXuU045xdprr72sM888U7iWv//+++L1Rx55xNp1112tWic353TrtGEeodsZsGzJusJOCOEkyDFxWwpbbvJou4OuPyeK0K3R2MttIVsLSyYUe1HzarP9irNZlS1VI1UWCEnodpNb3HCDelKMaOkKKtM9C5ybvdV0wbCFbtMkLs4xSVFiE7lfUsrFiR3XKRmES6S8gbnyyuIzRO0jyrxKHNdF37BtWnwW/aQ2JXvusov6nr3jyiOkmMZ41dn9Qxsn792UebLHt9xyi/G9yN4CTr2DlAJ5LHOmTvXtD/K8VPJZA2Woto/aytqwFjzZU8QvEWAlSolVOij2ur3dOt5PsSlbDe1r+mbJlq1UGi8Ok9L829+ab2QDhG6/uVsOqyn3iEsn4ZdKKcFrjkqREWX8iv62aFEo75Wi8JrFi/2zNCuEblkhK7ftZSNGRGsnhefYDww+qzr+0m8NKjnVIKj/6NYFzXvl0Cqv95szDnXJMtMqTj8st4wMOceq9kfeUAdVjhvZAERBHpiSEjFwTNhj/5Qf/ci4396lymCecaE7kXO6ZT766CNr1KhR1tFHH209/vjj7uuNjY3WNUGDrwbIjaVbNxF63J5Ui6SzaUhau2iyQBdZkDVW+lBCtyJOhRNJeV2jZMHJ9F7ka/pOmuw2KE+WmuQWfoqBcotIguRZtMPErkZ17TO5rvyszkm4D4Ypuk2afM6yynrhKClUfW9yCAHnpgy0geySq4yXDpHcx1sKijizkk3ZpEmRlYFeRYHOgieXrFhYTIuwtj38sK9bsm8iK2f+C3nvImZZsryE9dzhRIV+SqVKFlc5o1pnPGusLnFo0fnInjGiEk7lOcLtnw8+GPkeRgcJP7p9Qkihuyizt+EZxiYJv0qOPJKSdzqKjKjrD3tWmCQKlMdEGK+wotMuGM0JG0Ft69tOEeemoizrBkdVum3W1BSr0N1q37sq3rhovufnHtHjI6lyWorz1ByD0DQ5x41TWnz2Ut6TPvi6pxt6TXAJU/96VVx3xoXuxCzdM2fOtP7zn//EUbeqJvdCt92xg5J4mCRFS7qULHYKly0TF96gxZIn0GudGJeEFQ7OpsGv7eVze5OK0y2nyInUYr2uwgKvc+nOUmmeNcvtf7pN1ByffsXCvInVZW4Klm5dUcViOvNOOTHBsmVFZfEKiusNqrMcd+lXT8fFOyuCYNj7LMdiXIhwz14h0xtec0rQWfCNjZlpZ1NLt+9xWjpBXSN0F50Q4MyBd91V1n0cX4bQLbstzzaME752woRYzkF26uUK2XbGfe9aHjUmenLI8C2VEONVoCpPuwjI5TEv4j7L8R6L3LdDKBACk2B5cxr4Ja+03x9GoSeMD5qEvmmVNOcpeQ3z87aQ91F+3preBJ6s1A6TFC3MsXeky2Beq0L3QQcdZHXs2FEcDXbddddZb7zxRhz1rDqqRej2OwPZsWRVaiOvG+TjfaxpJUm5PC68ukXAb6F24sqSzAJ9Pgv4F14Y+D7hspxwXaIWXgTjrpubfGX5cvHcs7IBN6m3vMBdoXnf7T7JxExLIeX4dt82kDLux6XUSaLvOxsSPyWHHJtpmsgqa6XS4+coRUyhqaeO6nzkNMtI2VLpk2Xc7zgvv4SDWvdyz5o1e+zYsu9FGWOs2CfI98heXsZWYEkYO/h734tV6A5Snsd9SoWqHKgTJn1yvDjeJrp9SBwebFEUYyVnunvy1QS5TTer+pImLtx1dfb0O102fN8+llCYXV5LvaQM8wsTGhvQR7yJScPmGGDhmfN3qP62g2leAbsPueNDcVpAVcd0f/jhh9avf/1r6/jjj7c233xzq2/fvtaYMWOsZ555xlq/fn25da4KciF0azIIy3EyQUJ3SwUT2hxrak2zrQbKM2jlc34Vi4Bzjcllxs5WqpzIR3xkIMZRVebccINoI52AGVVz32q7+lOOimNNMk0glNdymmm4QUwbJF3yvThKkOLivBA5J+Isee8jshUmzLObnuF7KQpv8GTt9T1Oy0fo1q1f3jjdOMrNBkJ3kbCkOJ4rqDhnTsdRX05wptvDyMJdJRPtcTLNohwdBp5J8lhIsr6hM7dr9kdcjrGzWo8+4QT/2F+PAcTxyKszyJietfjsvJZy54hfSlnRXSVfyGuwwB02pKSgELqLTrrwzLE1kUjNYd26dSKm+5xzzrEaGhqsrbfeWhwf9uCDD1rtGXYBSJpcCN2es3ybFFZhP/fyvT3HZmShCO3xE0/4WhhkLasqPtk0Y2RjBu6Xn00WkoihBJe5GRsrld7cOZbhuCzTsstbmkfFFCr8Xd/PwHMu+9lJmbyz2tampc4vwZ+9zuoEqrNUrrySoCKvycLidOWViSiZ3KSZmj2bSgEQqo3KcDNVFeecbp3QHVcbNVSg/2RFgR9Wocl7j594+6ec5Vrhdai7x6JcBW1t2nC5rLVRnkrYTOjez8oKpVD9qLnZmj1jRujvHCv1odBnvFez0O1lwYIF1hVXXGHtvvvu1lVXXWXVKrkQuj3nPpOhVr13BiaQIDejuzIa54yCUsvFiSMud/MkNnjjx3/j7paBe0MJ3xcq+blKF+FmayfHClqLdNmES9xwFy9OzF26KAO1R5gtV1F21IEHxjpWXUuXSuhevjyT4VZZV+B7y2wDq6ajMHCP7vJkXZeTAtaZClltbVV5WkQWysQKf5+TEC1Kkrt6OTO9RllXEvddi0K3zBdffGHVKnkQulXnPpd06gq7acU22MvI+JlFC2lzxjJ0o+SvFDI0NssJSenjsQ7mbYONoi+8Hu2egXqEKQW/OGlbsR02EZfKwsjjZvzw4Ynfj6iHJzP55Bg2+3FdS2npkoXv5ctzpXCvtxO0Uk7LHMX52d6jG03maPdM9blzM7NWVVOpdJs6WeUF9tHCUU+IaIWl27I2bNggzvJjt/LjjjvOOvbYY93y4x//2Kp18nBOt8n5nllzHw9Tbs9AHcotR2bAfRYl3pLWeNqXyNovA/fvLKa6BENhBJTWnG9YUb4pQ3K61pwQ03GJBSmhWZqKpDrnpI8YXbXl2OU4leFu0i4pkR0n8MyL0Mbtwi7aeamvt3iP7vIeL+q4m08O8TyhRM1/ud0brmKfIlSIaOlmuC8VHQ+YsZjuxM/p/sUvfiGymB9++OHWsGHDrDPOOKOo1Dq5sHTPnRt8BEQGBnCUUl+hrKVJFxYu8qS1r2Tpm4E6RC0/TOl7987AvRdiiG8vlCm0o6DEVa7VWKWjZNZ3hNMkj6M0Pt4wZgWAcyKA37FwYYqctCvMedoo8ZUgIYqtk6bx6k7/iBqCsG8G2gPl62dXdJKQ7YVyYJhrFL4OEbUUgnuzoxTMuezXgULy61//mh566CH6wQ9+EPajICM8+eKLtCHgPUF/zyrHEtEMzd8K9v88MrLOJUQ0Ou1KZJQ3Kb/8KaXvfZnSh8fdUTFc4yKfv+9CREvK/A4ATBgvrSvcL+uJaBIRLSOitSGvxevtWUT0NhFdR0QXU+WpI6K+vXrRSvseOkv3Vg4b7HnvtzHV8c7166mBSNTzrDFjcrtXyTNBfWLDhg10g2HfudTub09H7GsvRvgMiB9+di8Q0TZE1K9zZzFGeY58Jsw1LIt2J6JmvsY771BD//7i9Qa70PbbU97hOSwUXbp0oR133DGZ2oDEWblyJY2YOJGqlT/4TNz8+rWUD3gjcUvalQAgZ0DgBpWG15WzbWF5HBENJaKjI87542zhPQ2uZ4V8czP1tu9hsF3i4H8M3uMoxf04kYjOtH9mQQ0Cd3b5vcF7drMVVUNztDcD/uOTnyXPIU1EFEXSGOxcY8AAmjlzpnhtpS2Ir/z3vyn3hDWh33PPPdZJJ51kffrpp+VY4quWrLuXB52/Xe3l1AzUAQUFBQWlukqe3Zz3tWOu076H7wT83TkHO87jwVCSK70yUAeUdEohplCSpgsvrKpzugv8Txgh/bPPPqNjjz2W5s+fTzvssANttNFGRX9fuHAh1TKffPKJ8AZYvXo1bbHFFpQ1Jk6cSOPHO45xAAAAAAD5cc9kl+J2InrF9g5Yn3alQOLkKTwQxEddoUAbJDG1vr6eVqxYQQ0NwuE8d7Jf6JjuYcOG0SuvvEKnnXYadevWjQoFE6cgkBU+//zztKsAAAAAABAadinfx/65YLvFD7LzVvjlewD5ZQ8ieh3hBLkibD6IgvinIOK6Herq6kR+AJn169fTm2++mTmh25TQlu5OnTrRk08+SQcccEBytcoxWbd0L1iwgPbZx1myAAAAAADyS4steO9FRIvSrgwANU7Bjs0Ok+TuNCIaQkQjbc8VtmhfcsklwjvXS0tLCw0axCM+f7Jf6ERqPXv2zKQwmTW4Q+yyyy40bdo0yhLbV0H2PwAAAAAAhs0Ix0PgBiATWBGyyt9PRIcR0Qo7adqKJUvokEMOUb537dqw50MkB8t4LOuZKgFCW7offfRRmjp1Kt1xxx0iphvky9Ld3NxMQ4dybkAAAAAAAAAASJeJ9hFygvZ2WvnRR9S7d+8iF/Oai+nmWO5PP/2UdtppJ9pss81KEql9+OGH0WoMKkK/rbYS7g2IjQEAAAAAAACkzWVEtJSI7uFf7LO+Tz/kELr3L38pkkGzJnCHIbR7+c0330wzZsygWbNm0W233SZ+l0slTPlsYd9kk01o8ODBwrffjwcffJB23nln8f7ddtuNHnvssaK/s6G/sbFRuF1vuummwp1h2TI+AbJYkXDqqacK7cWWW25JZ555JrW3c+7M/NHQpQudnnYlAAAAAAAAAMDmXs49Zf/M53Pf9/TTJPPrX/+aVq7kv+SUOM8pS/rs7gceeMDaeOONrVmzZln/93//Z40YMcLacsstrba2NuX758+fL855mzx5srVkyRLr8ssvtzbaaCPr9ddfd99z3XXXWV26dLH++Mc/Wq+99pp19NFHW3369LE+++wz9z2HH364tccee1gvvvii9eyzz1p9+/a1Tj755Fye0906f37q5/ehoKCgoKCgoKCgoKCQVG5m0ZTImq35+5w5c6ysYSr7hRa6R48erXy9vb3dGjJkiJUk++yzjzVq1Cj39/Xr11s9evSwJk2apHz/CSecYB155JFFrw0ePNgaOXKk+HnDhg1W9+7draamJvfvH3/8sdWxY0frd7/7nfidhXVuyAULFrjvefzxx61CoWC98847uRO6591/f+oDCgUFBQUFBQUFBQUFhaQyvoqF7tDu5ZxIbcKECSWZ5A4//HD66quvKCm++OILcT64nM2Oz3Dj31944QXlZ/h1b/a7ww47zH3/W2+9RatWrSp6DwfCs9u68x7+n13K9957b/c9/H7+7pdeeonyRj8kvwMAAAAAAABkjEm2a/n+zvndEoVCgfbbbz/KK6ETqf35z3+m733ve7TVVlvRBRdcQGvWrBGCbIcOHejxxx9PppZE9P7774tD0bt161b0Ov/+97//XfkZFqhV7+fXnb87r/m9Z7vttiv6O99r165d3ffoMtnJdOzYUZS0adh+ezqSlSdpVwQAAAAAAAAAbDjR8zVEdCIRDSSil50/EAmBOwuJ1NatWyeKTubTEdrSzVnLn3jiCbr66qvp1ltvpUMPPZQ23nhjIXB36tQp7OWqFj7PnK3mTpk0iXU36XP8uedC4AYAAAAAAABkjjuJaKhH4Gaef/55WrDASbWWHizTyTIey3yJCN3M7rvvTn/605/osssuE8eGVULg3mabbcT5bG1tbUWv8+/du3dXfoZf93u/83/Qe957772iv7MbPWc0130v09raKs5rc8qll7qnz6UGd9TfP/FE2tUAAAAAAAAAgFDMnz+f0oZlOlnGY5kvNvfyPffcU/jRe2F36XfffZf+67/+y31t4cKFlARsTR84cCA9/fTTdMwxx4jX+MB0/v28885TfobdEPjv7Abv8NRTT7nxAH369BGCM7/nu9/9rusiwLHa55xzjnuNjz/+WMST8/cz8+bNE9/Nsd86+HgxvwPS02Du3LlpVwEAAAAAAAAAQtO3b19Km6ghw0ZCtyPkps2YMWNo2LBhIqnZPvvsQ1OmTBFJ3IYPHy7+/tOf/pS+9a1vua7c559/Ph144IF044030pFHHkkPPPAAvfzyy+KccYYVCSyQX3PNNdSvXz8hhF9xxRXUo0cP954HDBggksSNGDGC7rjjDvryyy+FkH/SSSeJ9+UJPoscAAAAAAAAAPJG586dKbdYOWPq1KlWr169xHndfIQYn53tcOCBB1rDhg0rej+nlu/fv794/6677mo9+uijRX/nY8OuuOIKq1u3buKosIMPPthaunRp0Xs++OADcS53586drS222MIaPny4tWbNGmX9snxkWGtra+pHAaCgoKCgoKCgoKCgoFCIUl9fL2SZrGEq+xX4n7QF/2qC3dM5qJ59/LPmXs7MvPxy+vm116ZdDQAAAAAAAAAwoqmpicaOHUt5lf2MEqnx8Vh8ZJcpvXr1orffftv4/aBydOvaNe0qAAAAAAAAAEBtuJabxnRzIjHOUM5SvAkffPCBOFO7lhk0aJDItj5q1ChRssK1s2alXQUAAAAAAAByyZFEOH43BVatWkVZYtq0aaKYyrxG7uV1deFPFnvzzTdpxx13pFojy+7lK1euND5LDgAAAAAAgFpl/KhRNHHaNBFQzPA5TtcT0RAi2iflutUiLS0twqhZ1e7lfDxW2FKLAnfW4UPlAQAAAAAAAHr2J6JrLryQ/kVEc+zCP19ERM+mXbkaZdq0aZRnwpuwAQAAAAAAAKBKeam+nlZ++SU1ENHxduGfme+lXLda5d5776UFCxZQXoHQXUPsvz/r7SrH/6votwEQjQlpVwDUDIMTvDYWcwAAiA+O033zn/9U/o0dnIeVOV8fW8bna5lHH81vND3W6RqioaGBjh06tGLf9xcimkxEZ1fsG2ubgWlXIIfUE9HGIT9TWdUVqCYW2DGBSbAhoetWKzcT0ci0K1GjJDUGAIgTTobcd6edtH+/h4jmRhS8XiSi0WXVrnbp3r075RUI3TXG/9t774p9FyeeuJSIfoZFtiK8QkS9065EzgTuO4no85Cfu8B2M6tmqnVhGJ7yXMSC8YWeOvDPmB+TG+M/8Xn9zBTqBL72+OBnAKLj9OGwc8ceVLv8MkR7if3BrbdSw7e+VfT6SiJqtv9nfkhEYyOsA2uJqF8Vr7VJ8sMfcqvnE+Pn/eGHHyZbE1ARtjY89i0uOIn+3bYADpLnbSK6OO1K5ABeeF+wN91HRfg8J1RpsV3MqoU6abMxI8JGIg/sdtFFqc5F3Lbn28l45MQ8v0pQCLmEahPuvyuI6EEiapL7d12dULZxbGY7VRfHHxVlNqs8bOX7adqVyDnsmvxQyL0VH3N1GdUuUwzaq06aO84cNoyoUyciPuTJsmjmtGnCsMH+ovz/TPszJ4SsB8/1fe05iNdaKF3NGTZsmPDazS2WIR07drROPPFE689//rPpR2qS1atX85i2+vfvbw0YMMC67bbbrCzROneuqB9K9ZZfEll1itebGhut4zNQv6yU5t/+1l5Kydp/t92MP9d0zTXiM60ZuIe4CveX6UTWHPu++P5aMlCvuMvcG2+0Cgl/x+12O44nKvquukLBusvdvkll+XK3PzVz/4q5PkdloN3TKPVSX3b68038/+9+575WTWOYS/NDD4l7ut1eB263x3Xa9ULJRmmuwj4fZznBbh8u8/j/pUu/2Tu3tlp1dXUlcwzP12HWFP6Mdx2YnYF7z0upr68XzyIrsIzHsh7LfFw/lgH9MBa677vvPmvo0KGi0/Xu3duaMGGC9dZbb8VR56oUuoMaPi3uGj8+9UGDkmy53Z7U6yWBavLo0VbrokVKYbzaCt/3sB//2L1/VamTFlSewAuFQqjrt+Z8M7uz1D/43gtSuzgbgnmazyYttMrPKO5rzrnhBuP694xw/bGezVSrLYDPcfpbe7tW6HbeH/d9Fyr4zLIoZFhtbaJPO+3Ke5jJzqbafmZUBUXMaYsWlfQvbOiDC8+FO2WgHkmXJrtPTM5AXbJWCvZ84J0r7rrrLrFPmDdvnvZzpt9xtqwI5LWAWb48kXm/mktzc7OVV9nPWOh2WL58udXY2GjtsMMOQuNwyCGHWA888IC1bt26cupbNWRZ6A4rXKDks7R4Nvy86WIPh3l/+EPqdatE4Q2Fs5DN0Vn9pUVPt5gGbeav8Pn7tzPQDkGKGbd9PHOCo0nWKWkO3H332Ovj3bg4FoQ4NyJ8zdk33JCs0KOyZHs3WV7Bu63N/Vmn6EirDCeyRmagHlGft7Ba+SgbhUKySja87MmkUupk1bJZl6F6/MTn77JCsn8G6huX90fc82veC7eHSvh110OFpTtskZX9Lvb8f1cIAZ69WC6I6b7zppAtFAqZsnQnLnTLPPXUU9Ypp5xibbbZZlbXrl2t0aNHW7VOloXu2bNnpz5gUCpk3VFobCc3NpYsKDzh+lmEK1m4LkPj2ljMn+9uOuXNBX/HWc7GwxaCwiqjxMI5f36u3a9bA4Q8oUlua1NaRLxCuly2iVCXQzyeGfVsXZBcguPYFPA17po2TYyBJNpT5TJoJHR7XsuaxSNLdQlbb9djI0DZmLdNp6pwKIPoQypPClYsjBmTeh295dBDD63o9xU0r50X0AfOtC2U1dBPuDQ/9ljR2ph2fTKjqA9aD9lTdNq0svuBaH8Ze8xODqtMjKlPji7js3txyBYr4iN+fmiEzwwZMsTKIhURuh1+//vfC6GbN/a1Ti0I3XEtPtWyiGWpuFrZ+fOVGlu2hrjCTX29WERa7VjHStZzlL3gu3VxNsktLULQ2q/M68+ZMaNE8VBktZDcxpi7rrzSuD+6G1zLsgbvsUfqzzyqG50r5Hnj1JyYqfZ27SYkTpdcXrTlmObWv/wldssv33NLS0vZlgpdzPRcyV05UOhWIQlL3Bfl0JCy7jtAmZS1OXjHiM/2MiKrkcia6LjyL1rkxmW2PPNM6OeetXYJKhP9hG5+ra2tRMBK+x6TGIt+ZWzCSiQesz/IQF/wraMzt7e1ZU7Bl1ZxPZTa29Vx21IMMVupCwbXGxvG0s3fG2I8OiECzv6m3PtvzFn40k9+8hOrJoXuFStWiLhu2c38d5ygpMbJstDNm85KDFBnAHK/mJyCQBd14qCcF5GoyYk/uv9+30Q7zU78n4+FjTfsSVnBHSGlRUruIhaSlpZYNLhO3GbBNCGH5I7uJMKa7edJYAtQs2+9taLPWPZMcNwiC2V4Q4iF+6abihUxkjJi8rhxynaLywItJ7Bp0cQ4x/UdHBYV97PwtoGrzNEJPybYm7DmiK7ddVIfZoFTJ3jXZWzTfSqPuRkzYrFwc5+WvXyGnXpq6LmskFElg7Ku0tyvxBPX7lh4037me++6a8W+i9e7uOPbhTLPHqvjzz03lmSWQd8n3ltXFzpUsKiPtLVlLpQlzSKv6WyI0K2H8x59VHuNgUTWJZIVOiisTSZMiJu8dperONHNcVlaF8hTJk6caNWM0P35559b999/v3XwwQeLztirVy+xkUFCtXwI3VFiV+VSb1v5gt53uyNI2S6+KmtaJUvB4O+quD6vgMOblAlE1qAMTDzKdr/8cvdZ6yzdQssqCwGyhc1reb7ySlcQjbOeu8pu71633JaWWK2bxgk5NLGQqj4xW07CNn9+xRQ2jjeAqzSR6sljLkw9xo8frxTygrK1Otnbk9qwDUtQ6ObFOpbr+OQK0I6zMEiKMNONtco13/lulUtknR1mEffzkwXdKN4XLBiVO5748yprVRjhqGDHs1OVZPVVxbWnnd8liuBYTmEvrzitu3JISVzz1BzTvdWiRdZdkyaFUiTdPmVKyZxfDcaG0HODN4eJ1/qsWQ9NLd37K/ZTIncEn3qiQbfWqnKdyInY5s2aFXsbOXtdymhpaWmxakLoPuecc6ytttpKHB12wgknWE8++aS1YcOGOOpaVWRZ6C4nEQQP9iDLoTNgWxUJgsIuEPL1Jlcwa2VB4bYojpkpIxusbF1IdDKSPE3Y0hNogVO4IxYtNvbzi1vAKlEGyH1m8eJIG6NCuZtUjWWS+63q2m57Ll9eYkEamcCzvtlpo8WLRV+80emTvABpFvnvf+c7ZslIVDHGPko6jklL0jWx5Ykn4ncvLxRiCa9x5rd5IeImowrd80Ju1mVFp3eMe704wh5zE9bq1/Lww6EEKhFzf9NNoq5+872Yk08/PVLfcyxEOi+gas7qW8kkmqb9auTIkRVtG56znXW+HA8ux31YDiWJa57yHjPou7eyvbPEeAto94ImgVcls5inLeA7SpIizy5n3lHMv8q5m/cEBvdSEjJlsA7I9RL73nHjSgV3TpYowaFYKsOAqk4H2J4Ut0dQ0mepzJkzx6oJoXu33XazpkyZYr3//vva93z66adWrZP1c7qjCr+jDRaWkmRCktAdlE1aO8CkY12SnCQcd+KgLLfz7AnV9Lr19oIY9r7LSdChSwzmxOn6Cd2qv8UpYBWCkotImTxN+6kuA7GfJcXrNqYTulUeA/I1ZG+OpGIGXaVEe7s4Ck3+G/8uK0zCCPzu4qXZZPjGt9mbjyTCD26+6qpv7mnRorL7mvCm4NwFZWaf9Vq2gizdqrNejQg57sT3qeZcH8VaEp5HRULB8uXK0ITA+1i61HedmWMrnoL6nnfcy20Utl5xxf0f73nNL2N27JZuRd+Pqw8I66FUB342LYbzj0kddo2xjZo9wpCc48REIC4KhfKMqTis9n4Cd8neyjPW/QToQRoB0qrg6QQ/tBPSVeK7vG1apCSx49m1wnCA0O08b56LdtF856gRI8KHF8kWdsnjS04EW5KLxmMUcPZCSkOB9H+ej0mckzGhO7FzuoPczm+44QarW7duVq2TZUu3YPnySEdFFGxh029RkK3BKqHb+dn0+73HBvF3H5zQpBynUCnHtAedeRxXkV2kdBY9MVnpFhR7wneTQUmb9TBHWZiUQPcuqT6OFp81tN/TXK9F5R5vC9Xe2M7JGrexEqHbbpuxP/95sGXJts7HlmBQiqWX+xDH5+ruP8ozChK6S+LbbOFV/ozzjMIoooIK36fTF3lTHOUaE885pzgxm3M/N91k1E6uECH3Gc/8phL83H43bZp2o2SMR7HheBvxhqnEUhNC6E5yLmr2KKzC9knHi8Iv2y2PN6fv6XKMHHjAASUJI5NUOJgqX6fbIUrOnFXu5rdorQnoYyUWPs7mX6YCwvHAKfKQMsjG7CgIvNa9kjCIQiG2UyJ4LJcIzLYAxgJUwS/8bNw4tXJang80XlFxFLZMloxrz1j3G9e6BF6uwqBC40B4YAXkdfB7FmHHhnM0pnZfqhKE/YRuT7vpvnvugw8W76cMhW5V3wxzhJmj+PR636nCi+I4/uz2MpRNdWwUifC5pqYmqyYs3XwO97hx46yBAwda++23n/Xwww+L12fNmmVtv/32VkNDg3XddddZtU7WhW4/y53JIPP7e/Nvf+sr1Pllk1aVkRVK6BDHBF9na3HliU7E8FTgKCDhmim1t6/QraFISHCEK1vwVbkvRQ5T4Azq11yjTVTi11/q7GQlJTHA0iJaEosla7Vl11sDl9+go8Rcy9Ly5bEl6BFtP2mSss43amKSJ0RRohm4lzt/01oFPJuEYUcfXfb9Dxs2rKgvRi0td9/9zaZHErqd/nBYwHhyx/CYMfqNttS33ISAnDU7IAOuMR7Fhrx5VCVDNNpEhhQ8eVO1h+nY9mwwgxRW2vHU1uYrjI4fNSpw48sCVotCEEwzeVRJHou2NusyTXy5/F4jK56JN4UqVnX58hILvPHzUgk0HiWPLgfM8ccco5xfSvKKTJoUy9xapPiSx4O0doz12YfIRgOdgovrGvSsoggoBWcO8/NSM8hGrjuqqpLWbpM+FZRMzrTN7gpzZKNuvvT5u24u6d+3b+l+KorQzTHbmsRtrOzXhn951gjd+Bl77rll7X2bIs6n4+W8BCFOjam5c7ovvvhiq0uXLtZxxx0nhOwOHTpYI0aMEG7nnLX8q6++iqPeuSfrQndSsV3u2cg6dAJcoaBcKPi1uDONqsopMV2HXQWDEpd5NxT7hpik/P4+/LjjAoVFEc8WwvXQ2QDrJvaxI0aEdi12LaU+iUqC+ouy37HgoVtgNd4WJkK3X+LBImVBe3tsfVVkmJUtl1KddZbu6dOnh/sOfg6mlleNB4Dyb4sX+8YkBmYvf+KJsl3A3f552mlqS7PUH1goOzHAIlBkHePP+rWN/R7tZkgTb2vc9qpiEi8uCaiO5cVr9VSFZzj3bmJp9CZw43Ed9H5WmCm9KGwBwm+MyPPq2DPOULe3oo38+pfjdsn9MI6x7Hv/AZnlb7dziZgIIEpPoaD+ZLdJVCV80fP2zqvS9+g25e565O3jUpyyGHdtbdb0K6+M1MZzbavpHEWiLJXgrEvi5+5tfITuIEXWBPs9k0ePDnUPrsDGbez3PKWwLNXzVCr9JO+2OATduEohbaE7xHyqamse08rwliBBUbHWcr8NY+mWFZdFfVM1/zzxRFnKbWeMmq7XTpy6JfdZO0dNZC+9aha6+/TpYz3yyCPi59dff110quHDhyOZWs6EbtUgjqM0qTalMj5n/s7RuEzGYR2+O8Jn6uw6jQ/x/XWG2bFlC+Zcw6zq/Dk/N0CxOfZM6CpXQh1+QoJyYrfrL2+O2O1Tt3mUrYDGi55Pf9Fq8FXXLUPoVt17QbOJi+KiV7CPGCkExehKDDvhhBLLcBhBdYL8HEwwiG9zNwiaDfz+pq6TPgJr6DGs2/R43a8NLALNBkK3K9AmYOkuW+j2erKwwMQx7h4LOgvBRX3r6KPd7xkWZGn1eJH4PcdRnpMtSrwo7Pu+LGjMO22vSCZUNIYCQiZKwgf4eK0EXYWd0qiZ06OEO5VYMYP6k90mYZXw7vw3d65+XvWMi4Lf2qgQur3jKeqxjEVx0DploTwf+HhYCC8+7zwgCa2zDT0CwmabFjkMgvZWDva9OPlFgrzJ5DkhylnLWU/EVeKJEUXADhhHOiWHtr8bXleeC1Vhc/Iz1CkurYBjQXVW9AkGz1aeX+UQvsBEsMuXlwjdTi6aUOOiFoTujTbayFq5cqX7+yabbGL97W9/K6+WVUjWhW7vIDYtQZt6+ZxDHdrNqGfz62i4nUktarImFlCEJSSCBt8bE/OTo44KnFh0Gkmn3VWbC69AMlDWzsuCwfz54SZ02a06QMgKEhJKJnbv5GkisITVNIexdDsCsGbRUm6wQggrJfeu2sgrYoVNNiW6M+xVVjrnubbYi5jI8q2oZ50TmqFSlphu4CII3VGyQjsnBzjPJ8gSWc5GT7nJ9wjfYYQ3ZQiEbVX3O+s1TNuX5FkIeiYR5l2je7eVZzd7LYjys5MSx/F3+mb25T5o4EVxpMmRbJ7EanIMt+65qVz2i94vJf5MUrhQWbrDfl85lu4oVv1m7xwSEFKmiu2WEzJ6BXbveNK54Bu3TcA6I6/Jyro6yjRPYkRvGIzqWXqTy4bxLPCGa4R6vt4QFEOFsqrvyUenlt1X6+uLPFxCjxc7qa/piTAm+9KyMPRyi6R0lb5D6xmoC//yKJNU11Ba0X32cm7+CDuUQr6+o+hRell41glL6tNwL/eBB+h7773n/t65c2drucrlpcbJvNDNtLeHSk7CnXzYqadqJ8owE0qJECO70srF4x4zxzAe6nZHIOHjsyIkCpkzY4Z2A9oSMLGwtV9rXVYkyXB+nmtbf5xjJlTFT6BRTkLeSTfMc1Fk9pZjVx0LhKzt1dWtKLleBKHBW68il1S/DOQxupIpN+gBLqsmfdVJrhRK0NM9V0U95ZMKijZ/cQndHqKcsSw8ZDz18j5z2RIpC1dh4yPF+Z4B1mM/q4Kub+jCSgJDKALQJmML2aeDYv+EwkRj8SxSmvm1mydhIdd17Pnn66/pjUnVCGxc5JAF5TORrKPeGO5AZZGP0O38HtYlONRRaeyKGnLMiDY27aNBlm6Ntct3HEUQApsM66wSBMoNNzEVvFRnmYs5qqmpZCxOnjzZ91xl77FirsA+f75RIll3vg4zV3ufsc8coZsTHK/DOk89orgAK62d9tgs58hC05xAZQm6IdpbFzZZttJV+g6jucyvj5ta0TUKWCehYZFRTHFtpZeFbl+/dGmkHAfjx4+3skjsQjc3zg9+8APr2GOPFYVjug899FD3d6fUOnkRusMmQHCEzuYg4dLgu7WbI8/g9GZ/DHL5K9Ko8abJ4wYUeI+2MKlrGzlBhdJ9y0AjWWK1sjPxBhU/gSYOoTuo7k4bFm2sDRLEKS3dYfDWK4QFPxZ0/VO6F90Gxq9dZCE4lKDnI3SX1NFzjqvb98Ju5EyJcO5rUf/wuG26z9g7FqS/DVOMhbFjxwZ7gyjayxmfRQm4AvA7y7ycfh+bi7qJpTsgdlBWXhZl45XaTXV9jllWWVLc+c97XUXcrfO3ubZXiOzd4aLK2O6HzvtFI3TzHB1GSBh/wQVF6+MpJ57o+/4owsecMhU6DqyICvvdoa2H3tAqn3k7ihKgLEu3/N0axZPqaDOdIsB1v/co8upCKFlcV9wEhe5QXocxJCCUn4MuxNErgJUj5Jct6IZsb5WnTblK10oI3SXjUrUn4d+vvNJIQapao8Pum8inHDlkiJVVYhe6zzjjDKNS62T9nO5yYlDl2EalS4vhdwcOXNbAqbI/2kKELPDKk3qRRs0jdDsTwVE+9+dkwjZxtSyaWPzcm/2yg3uFbrbEaQQBnrwmh3QvDyt0+03sumQcRrFs5bh3KdwPyxLiy/l+jXu6buPKFrp6b5ImjVur8QIdUuguUZaEVZJFaK+ic0XtTOG+bmd+Gwlvmyv+1mLHod3u5BoII7CqFEqm2Wb9Nq+y+3OEfhpbMjadJ4vC/VoV9+d+ftKk0vYJ2JCz4kHevDlnyMrPUuWa736nIlZQefyan8Vchckc4hG6LcNkU/Ia4oxzXQJE537CZrwO7XbsQ9iNbyGMpVvV3gGfjeo1FIfgpVM86U4D0caqy/M6XzOEpT62OGSDPl4yJ+i8DgOOfwsbW69TrKiE7rqI+XjKFnQjtLOxp01WhO4ABWTRHt9rHAoQuoP2TIxfCJKqDHf2NBkjlXO6Qf4s3c5G01TwLknuFHVhMJgk/JJ3eSeEIiHGm6TBISCTY9F3SPFaRm7w8n1o7s1XCJDr7JlQStxLNZnMdZlJQz2fgIldmwRvzhztpoiTfJVNloRuzff5ue9qY0d9FE5lC93ONQLO+UwESXsuuykWbOG7ZEz5XMNU6C5Z5AMESOP5xnDjpgzNKLOfxmnpFpicGKDxItG2j31vQYqHIiWMJDT73aNfbH8s851Be3n72WzT/ATe85M1QoZwU77mGl/Bgt/TNG6cOkwkhnuNclpA6HPnQypFivJT8Fm+PkK3LJiVKDbDCl5ed1tbQanrp+Mvuki9hzB49kXtKc8dcT1bk/Gg8CJT7W9Mcqv43Zs3xEwZQuDzfH3HnNeNu5J7A287J7VPiUPoDntt73357fn9hO6AOo3VhCBpkwpmUOhOzNINqk/otgxjdUosIwkL3VohRs4gqio6odszKfhlnJQzYRu5wRsI3b5WK43Q7ef+FRjjGIUolm52IW1p0W6KYhHuciB0G1k74xK6w9TTR1kS1WJqWg+dZ4Q2LiyK0M0/e+ci55ohwhB8XcQN7zk2V0KJWJKxhdyEh56//OrqF/Pucxyh3xm02v6b5Lwg9elC2GR9GqWEEEaeeaZE0NOte0bKu4jIuR9MS6i5PYIngnO/s++9V18Hey2PrW3s750tu4g7XiCa/l0Sq+4VWAMs3XLYXsmxZEmvb35rqyRw32jQHxwvlrEqzxbZSip5uJQcWRghdr+psVHtylwpodukLStdlwSFbmUokO4zBvfeGuI0pdMgdFeWDz74wDrllFOszTffXJwX/rOf/cxas2aN72c+++wz69xzz7W6du1qderUyfrxj39srVq1qug9b7/9tohV33TTTa1tt91WxAN++eWX7t95EVV1gH//+9/5FbodfDblzoSpFDrLEboNKFuI0dyr/P6WMNnGVdfWCDmqz0SxWmnPx/bUmzO7VgrVxiP05rhKFzOlZdVUA5zUPQUIPYlg30tQXgTf+wyj2dfNRSHCK8q1dCf9/BKPCwx4j9H8pamrLgGk9jjCqJbuJJHaxU9hq83IzfPD+PGl1mpp3Wi1wyNU1xVx7HJfjru/2QqqOfapB0VW5rjm9jB1lt6rElpdF2JJMFTlGgjbNt4jtCYrnpEqnCxI6JPXTVZQ1+myQHtOBKm40O153eQ4rB8qYti94UVKjyZV4k9pDTU+hkvOnRHhdJLYyMk+xRjVHKZL7GnwWb/vuStA6ehVsmSVqhO6Dz/8cGuPPfawXnzxRevZZ5+1+vbta5188sm+nzn77LOtnj17Wk8//bT18ssvW/vuu6+1//77u3//6quvrO985zvWIYccYi1atMh67LHHrG222ca69NJLS4TupUuXCkHbKevXr68KoVuZqVCXdVr6XNKLQ6DLZplCt7OohMqErbuewf2HtVoprSN1ddrsqhVBdexEa6vW0s31zeLRDomgsqyaLsQJL9CxWkzLtHTHFYvqohO6Ex6flRK6Y7tuGUK3UfvoBHa/BG0B15X/VvEESV6k++Okbqr57siDDtJ/3j6Ptsgi6+mrOvf1s3/6U2MFUiT4XHJPKIhz9rpO8KqE4s5pI+2xjbpcDBHGjE7J466tQUr2gHElrw2+sb9ywtWk45J9xrtJ2AG3tzfk0FjJq9rLSWuoboyV7C9kQ0yWhO68Y+CtE6cnY4uhoiWrx4VVndC9ZMkScTMLFixwX3v88cfFA3jnnXeUn/n444/F2eIPPvig+9obb7whrvPCCy+I31nI5o4kW79vv/12a4sttrDWrVtXJHR/9NFHRnXNm9BdInw6C5ec1CsGoTNK3YKyaUcS/D2f1bqgxn1/EaxW3g2pLv6lYsKtpk0mNzaqtZKVUgZkAZVlNSsLcULuz5HmlbStOEm0UQ0I3YHt46NQClJuBp34UJS7o1IJknS0tWndbX03hKpTKgyFbjG/JzhudOf1ijVROoO30oo7ryDqdSGPU7GnPYHCWVvLFLq1gqFCwWBkTYwDn/Gua4+is5o9AjP/bxzOFNCeJrHwYryVmbASpJPY08V+5mGy47M3chapKqF75syZ1pZbbln0GruA8+T/0EMPKT/D1m2VsNyrVy/rJtvd5YorrhDWcxk+e5w/t3DhwiKhu3fv3lb37t2FVfy5556rKqFbFj7do4Xk+KKUhO7Erm8itCcgdEdREMgbUj+rcmIxuib30NZW4lLGgnhNETZbfCVJS/j3y4uQQaE7ct3yIHQnXZeoAnvA5yIrV5OirU2skboNoXYe9sZEKoSNVr/rJnj/OgFHuO4q9ggVV3yo+kEYAc+AwLW1AkK3ypW+4qEUASE3utM3yrJ0K9rTNOt7uUczAjOStnS3hk3qnEFrt6nsV0c5YNWqVbTddtsVvdahQwfq2rWr+JvuMxtvvDFtueWWRa9369bN/Qz/z797/+78jdl+++3pjjvuoD/84Q+i9OzZk4YMGUILFy70rfMnn3xSVNatW0eZo1Onr6er9nZqIKJ/EtG+RDSUiHp/73s0M+DjK4momf9/5x2qOpy24cI/pwQ/lyH8/7e+RQ0NDXT91VeXvKe+vp769u1LaTKWiN62+wP/f9HFF1NNIfeRFPtL1pD7b0XHVLV+V1YJaIOifpBnOnWiQUT0E8WffOfhTp3EetrbWV/r6mjmAw+UtNHkMWPU121tTax/9evbl7wbwXoi6rvTTtl7hm1t7o+diUrrHXEtjLS28rNobw/9Xbpxs+ydd2jDhg1Ff16/fj29+eabVGm4PWZMnSr6gdMOdxKJvi/6gOdZOPfR0N5OM+z+437uzjvF9ULNF9/6Fs244Qb3OnV1dVRQPZvddsPcm1Z/UD3XsHQKP4bSGhNeWKbzynkmpCp0jxs3jgqFgm/5+9//nmYV6dvf/jaNHDmSBg4cSPvvvz/NmjVL/H/zzTf7fo6F8y5durhl0qRJlGVYgD6LiJwpnyf/kfX1tJIXe8VkNvPee7/ZQAwYQDNnBono2VQ25I2LLriAmqSBG9vkZ4JuobR/djdlVINA8AJxgz6VWR4kovFEriAQNA+zYrpkfR05klauXFn0nC8677yKz+8N/fvTjGnTvtlQEwkBK5MKErutZk6bJgwEsohablv5rq1+a1+YMerz/n79+gnhMisK9TOHDaMVtiJ9xcsv05km99KpE53Z1lb8uTPPjPb9p57qXuftZ5+lX5kI86Ay/WHJksjPVcUy24xtQhaMTAzLdLKMxzKfCR0oRS688EI644wzfN+z4447Uvfu3em9994rev2rr76iDz/8UPxNBb/+xRdf0Mcff1xk7W5ra3M/w/+3tLQUfY7/7vxNxz777EPPPfecb71bW1tpiy22cH/v2LEjZRnu9Bs0GiXvxMYbhbNGjy7ZQBx22GGYBCtkVT6JiFjX13fJErFpAgCUgbOBBNlpH+c7164l6sx2zQzg1Om99+iabt3obGcefvllavjud7UfE1ZMw/U1jfmdN9SHjRpFL9ib3/2XLv36XrntM4Z3/0G28uOFefNo0Pe/X9a1i9o+4JnGPXYca+LIUaNofUYES/5m8e39+n1tpHDGIf/so2RwP9ejh/8XyGObrym309q131xn++2F0H9Y3M8GhMJ9HnEq5Dp1on4833z720Zv/+Uvf5kJOePSSy+lMZJnElu6TQTvVIXubbfdVpQg9ttvPyE8v/LKK8LizMybN08Ie4MHD1Z+ht+30UYb0dNPP03HHXeceG3p0qX0r3/9S1zPue61114rBHrHff2pp54SwvIuu+yirc+rr74q3M794GvIQnfW6WdreDcYaJSWLVumdYOKbTBUYpNXyY1kzN+VyOQXFwELMqgwEChBnvtmhvtvkVDig2PFlNdN5foqeQ1Ven5/UvJ2qxswgGbMmPG1NSttxYfn+S9raSnZf/Bf5zz0UNlCt3y9NHCUH5lUqJuOw7VrheckG3L68Z7QExoaFfeaPDYCxhuoTuoKBTr//PMpC7AhNYoxNRcx3QMGDKDDDz+cRowYISzT8+fPp/POO49OOukk6mFr0t555x3aeeedXcs1m/t5wWBNRHNzsxDYhw8fLgTtffdlxySiQw89VAjXp59+Or322mv05JNP0uWXX06jRo1yG3PKlCn0yCOPCIFy8eLFdMEFFwiBn99TNXTqRA2WVexi5qNlzZobVNWSF/fSOGLbQPWRl/4LkgX9ILmYyJjwdX/P2DPk/QeHHnq5eerUb+obkaK4+4EDUwmby0TsfBnMnD37mzbk3EBR21DqZzP/+MdvrmmHN4LqYZlBjDZLHDNuuy0zc2ZkrJzwwQcfiHO5O3fuLI70Gj58uLVmzRr372+99VZJ5srPPvvMOvfcc62tttrK2myzzaxjjz1WnLEts2LFCuuII46wNt10U3FG94UXXigyoztcf/311k477WRtsskmVteuXa0hQ4aI9PlVkb28jGNzKnr2L/gG04ypAAAAkkF3Pny562tK87s2g3klTsWIgO74zHLqG3SmfEXI04kFcWS5Nji9QGREz0hW95omwb7Zqhh7XEYTWXPsUnTiQwYxlf0K/E/agn81wX79bGVfvXp1rtzLlfE1fhptdiHq3PlrN6ilS7PlBlXthHlOAAAA0p+HTd+fwvzOFuLevXuXuL+vWLEik5allf/4B/X+9rdLwuHKqS97RA4dOlT5Op9YUxFyvraHbkOD+83EcwHJsnYtzezcmUZyqKpt1b6Okxt635fhMWEq++XCvRxkk7y7QeWWDLn6AQBATZLUPJzC/J5193flkVIxZ7NG2Fz5RGnDoKNn8Vxqg8OI6Lecm8E5djZHAncYIHQDAAAAANQwSR4JlARcuzjrmwnFQ84V6qINpfO16+vqfNvQ5OjZTDwXkCgz7X5won16ACd1rFYgdCfEoEGDRJK2adOmpV0VAAAAAICq8l6Lu755UzxkFTchn0/0qu7oWVUyPDyX6mWltx8QCTfz8lIiVg6W8VjWY5nPBMR0x0yuY7oBAAAAUHvkKZ7Ye4RZXPXNUxtkkDC5AULFauO5VC3Nmn4wZ8YMOn7ECMoLiOkGAAAAAADVBY6pzCTLli0rOUN9/fr14shdL4jVBrp+wJx09tmpHNmXNBC6AQAAAABqmZzHE4P0CSNII1YbyP3AK4z6hRvkGQjdAAAAAAAgP0BJkP1EagGCNGK1gdMPfkel6Lwk8gxiumMGMd0AAAAAAKDmWLuWVnbuTCwq9V26lBr69/d9L2K1Adl9hjOYy8EJunwAWQQx3QAAAAAAAICKkbcs+CAbfWYGC9pVHm4AoTshcGQYAAAAAAAAAPjDwQV5CzfAkWEpA/dyAAAAAABQc4RxGYd7OVAdAZjD/mAq+3WoaK0AAAAAAAAAtY2TDA/UNp3sfqASvqsMCN0AAAAAAACA8oAgDYAWxHQDAAAAAAAAAAAJAUs3AAAAAAAAAIB06FT9XhKwdNcg69ato//+7/8W/wOQJuiLIAugH4KsgL4IsgL6IsgC66qoHyJ7eQ1mL89DHUFtgL4IsgD6IcgK6IsgK6AvgizwSQ76oWkdYelOCJzTXT5Zb7us1y8vdcw6eWhD1LE2yEMbZr2OWa9fXsh6O2a9fnmpY9bJQxuijtXLNJzTnS7VpJFJG+7IS5YsoayS9frloY556ItZb0MGdaz+fpj1NsxLHbNeP/TF2qhfHuqYh76Y9TZkUMfq74ef4JzudHB0GPwAsopTtyzXkVm/fn2m65j1+uWhjnnoi1lvQwZ1rP5+mPU2zEsds14/9MXaqF8e6piHvpj1NmRQx+rvh5/YdQuyY8PSHTMrV66knj17pl0NAAAAAAAAAAAVoLW1lRoaGrR/h9AdMxs2bKB3332XNt98cyoUCmlXBwAAAAAAAABAArAovWbNGurRowfV1enTpUHoBgAAAAAAAAAAEgLZywEAAAAAAAAAgISA0A0AAAAAAAAAACQEhO4aPVduhx12oE022YQGDx5MLS0taVcJ1BCTJk0SZxpy3oPtttuOjjnmGFq6dGna1QKArrvuOpGL44ILLki7KqAGeeedd+i0006jrbfemjbddFPabbfd6OWXX067WqCG4CzWV1xxBfXp00f0wZ122omuvvrqwKzMAJTL//7v/9JRRx0l4qJ5Hf7jH/9Y9Hfug42NjbT99tuLvnnIIYfQsmXLKE9A6K4xZs+eTWPGjKEJEybQwoULaY899qDDDjuM3nvvvbSrBmqEZ555hkaNGkUvvvgiPfXUU/Tll1/SoYceSmvXrk27aqCGWbBgAd155520++67p10VUIN89NFH9F//9V+00UYb0eOPPy7OzL3xxhtpq622SrtqoIa4/vrr6fbbb6fbbruN3njjDfH75MmTaerUqWlXDVQ5a9euFTIJGwZVcD+89dZb6Y477qCXXnqJOnXqJOSXzz//nPICEqnVGGzZZisjT6hOtnU+4mz06NE0bty4tKsHapD//Oc/wuLNwvj3v//9tKsDapD29nbaa6+9aPr06XTNNdfQd7/7XZoyZUra1QI1BK+/8+fPp2effTbtqoAa5oc//CF169aNZs6c6b523HHHCcvib37zm1TrBmqHQqFADz/8sPCEZFhUZQv4hRdeSGPHjhWvrV69WvTVe+65h0466STKA7B01xBffPEFvfLKK8Ilw4FT2/PvL7zwQqp1A7ULT5xM165d064KqFHY8+LII48smhsBqCT/8z//Q3vvvTcdf/zxQgm555570q9+9au0qwVqjP3335+efvpp+sc//iF+f+211+i5556jI444Iu2qgRrmrbfeolWrVhWt0V26dBGGxDzJLx3SrgCoHO+//76I12HNkAz//ve//z21eoHahT0tOH6W3Sq/853vpF0dUIM88MADItSG3csBSIvly5cLt14O/7rssstEf/zFL35BG2+8MQ0bNizt6oEa8rj45JNPaOedd6b6+nqxZ7z22mvp1FNPTbtqoIZZtWqV+F8lvzh/ywMQugEAqVoYFy9eLDTpAFSa1tZWOv/880VuAU4sCUCaCki2dE+cOFH8zpZunhs5fhFCN6gUc+bMofvvv59++9vf0q677kqvvvqqUIyzay/6IQDlAffyGmKbbbYRmsu2trai1/n37t27p1YvUJucd9559Kc//Ymam5upoaEh7eqAGoTDbTiJJMdzd+jQQRTOLcDJWvhntvIAUAk4I+8uu+xS9NqAAQPoX//6V2p1ArXHRRddJKzdHCPL2fNPP/10+uUvfylOHQEgLbrbMkre5RcI3TUEu6kNHDhQxOvI2nX+fb/99ku1bqB24IQYLHBzkox58+aJo0kASIODDz6YXn/9dWHNcQpbG9mVkn9mJSUAlYBDbLxHJ3Jcbe/evVOrE6g9Pv30U5HrR4bnQd4rApAWffr0EcK1LL9wGARnMc+T/AL38hqD48XYRYg3lvvss4/I0Mtp+ocPH5521UANuZSz69ojjzwizup24nE4KQZnSAWgUnD/8+YS4GNI+Jxk5BgAlYStiZzEit3LTzjhBGppaaEZM2aIAkCl4HOSOYa7V69ewr180aJFdNNNN9HPfvaztKsGauAUkTfffLMoeRorvznJLvdHDnPg00X69esnhHA+T57DHpwM53kAR4bVIHxcWFNTkxB2+GgcdqXkDIAAVOooCBV33303nXHGGRWvDwAyQ4YMwZFhIBU43ObSSy+lZcuWiU0lK8lHjBiRdrVADbFmzRohzLAnGofesFBz8sknU2Njo/CWBCAp/vrXv9JBBx1U8jobCvlYMBZXJ0yYIBSRH3/8MR1wwAHimM/+/ftTXoDQDQAAAAAAAAAAJARiugEAAAAAAAAAgISA0A0AAAAAAAAAACQEhG4AAAAAAAAAACAhIHQDAAAAAAAAAAAJAaEbAAAAAAAAAABICAjdAAAAAAAAAABAQkDoBgAAAAAAAAAAEgJCNwAAAAAAAAAAkBAQugEAAIAa5YwzzqBjjjmm4t97zz33UKFQEOWCCy7wfe8OO+xAU6ZMKfrd+ezHH39cgdoCAAAA5dGhzM8DAAAAIIOwUOrHhAkT6JZbbiHLsigNtthiC1q6dCl16tQp1OcWLFhAzz77LB133HGJ1Q0AAACIEwjdAAAAQBXy73//2/159uzZ1NjYKIRch86dO4uSplKge/fuoT+37bbbUteuXROpEwAAAJAEcC8HAAAAqhAWaJ3SpUsXV8h1CgvcXvfyIUOG0OjRo4XL91ZbbUXdunWjX/3qV7R27VoaPnw4bb755tS3b196/PHHi75r8eLFdMQRR4hr8mdOP/10ev/990PX+b333qOjjjqKNt10U+rTpw/df//9sbQFAAAAkCYQugEAAADgcu+999I222xDLS0tQgA/55xz6Pjjj6f999+fFi5cSIceeqgQqj/99FPxfo6rHjp0KO2555708ssv0xNPPEFtbW10wgknhP5uVgK0trZSc3Mz/f73v6fp06cLQRwAAADIMxC6AQAAAOCyxx570OWXX079+vWjSy+9lDbZZBMhhI8YMUK8xm7qH3zwAf3tb38T77/tttuEwD1x4kTaeeedxc+zZs0SgvM//vEP4+/l97IFnS3r++67Lw0cOJBmzpxJn332WYJ3CwAAACQPYroBAAAA4LL77ru7P9fX19PWW29Nu+22m/sau48zjgX6tddeEwK2Kj78n//8J/Xv39/oe9944w3q0KGDELYdWIjfcssty7ofAAAAIG0gdAMAAADAZaONNir6nWPB5decrOgbNmwQ/7e3t4s47Ouvv77kWttvv33i9QUAAACyDoRuAAAAAERmr732oj/84Q/i/Gy2VEeFrdpfffUVvfLKKzRo0CDxGmdbx1ncAAAA8g5iugEAAAAQmVGjRtGHH35IJ598sjhDm13Kn3zySZHtfP369cbX+fa3v02HH344jRw5kl566SUhfP/85z8XmcwBAACAPAOhGwAAAACR6dGjB82fP18I2JzZnOO/+cgxjsWuqwu3zbj77rvF9Q488ED68Y9/TGeddRZtt912idUdAAAAqAQFy7KsinwTAAAAAAAR3XPPPUIwj+o6/te//pUOOugg+uijj5BoDQAAQOaBpRsAAAAAFWf16tUi4/kll1wS6nO77rorHXHEEYnVCwAAAIgbWLoBAAAAUFHWrFlDbW1t4me2VPM54Ka8/fbb9OWXX4qfd9xxx9Au7AAAAEClgdANAAAAAAAAAAAkBNTDAAAAAAAAAABAQkDoBgAAAAAAAAAAEgJCNwAAAAAAAAAAkBAQugEAAAAAAAAAgISA0A0AAAAAAAAAACQEhG4AAAAAAAAAACAhIHQDAAAAAAAAAAAJAaEbAAAAAAAAAABICAjdAAAAAAAAAAAAJcP/B1ihcTC97hQOAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Initiate ARVE object and plot RVs\n", "example = arve.ARVE()\n", "example.star.target = \"Sun\"\n", "example.star.get_stellar_parameters()\n", "example.data.add_data(time_val=time_val,\n", " vrad_val=vrad_val,\n", " vrad_err=vrad_err)\n", "fig = example.data.plot_vrad(figsize=(10,2))\n", "plt.show(fig)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We perform the VPSD analysis in 3 steps:\n", "1. Compute the VPSD from the input data.\n", "2. Add the expected VPSD components based on the span and cadence of the input.\n", "3. Fit the coefficients of the analytical VPSD components to the computed VPSD.\n", "\n", "The final results can thereafter be plotted." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAA0I9JREFUeJzsnQV4U/f7xU9Td6NQpbi7u47BYGxjTJnCfGz7Tdj+c+bCjDF3xtx9DDbc3V2LVahQ17T/53zDTdOQpGmbtkn7fvrcJ3Zzc5MmuSevnNetrKysDIIgCIIgCILLo6vvHRAEQRAEQRAcgwg7QRAEQRCEBoIIO0EQBEEQhAaCCDtBEARBEIQGggg7QRAEQRCEBoIIO0EQBEEQhAaCCDtBEARBEIQGggg7QRAEQRCEBoJHfe+AK1NaWopTp04hMDAQbm5u9b07giAIgiA0QDhLIjs7G9HR0dDpbMfkRNjVAIq6uLi4+t4NQRAEQRAaAcePH0dsbKzNdUTY1QBG6kjr1q3h7u6OW265Bbfeeiuckb59+2LDhg1O/xjV3UZV7mfPupWtU53bs7Ky1A8BfjCDgoLgzMj7pWrryvtF3i9VWVfeL/J+sXfdDz/8EI8++iiaN2+OQ4cOGXWHLUTY1QAt/bp582an/yBReNb2PjriMaq7jarcz551K1unJrfzenm/yPvF3tvl/eK4x5D3i3Mg7xf7150xYwY+/fRTrF27FsHBwXaVfUnzRCNh+vTpLvEY1d1GVe5nz7qVrVPT250deb9UbV15v8j7pSrryvtF3i+OfL+Y41bGijyhWjD0TQXdrl07pbr54rv6B06o/fdLZmam0/+iFuofeb8IVUHeLw2Td955Ry16vR779++36/8rqVgHwPy4fJCEyvD29sbMmTPVqSBUhrxfhKog75eGyfSzASNNuNuDROxqgPxCEgRBEATBmfSG1NgJgiAIgiA0EETYOQC2K3fq1EnlwQVBEARBEBwBdQX1BXWGvUgqtgZIKlYQBEEQhNpGUrGCIAiCIAiNEBF2giAIgiAIDQQRdg5AauwEQRAEQXA0UmNXx0iNnSAIgiAIzqQ3xKC4gUBX6hUrViAxMRFRUVEYOnSomobhbLjKfgqCIAiCKyLCrgHw888/44EHHsDRo0eN17Vo0QKvvfYaLr30UqcRaPWxn9XdV0EQBEFwRaTGzsWhWLrsssvQtWtXrFmzBtnZ2eqUl3k9b6+Nx2zTpg1GjhyJKVOmqFNetvVY9bGf1d1XU0G4dOlSfPPNN+qUlwVBEATBqWGNnVA9MjMzWZ9Y1q5du7KOHTuWvf3223X6+CUlJWUtWrQomzhxYpler69wGy/z+pYtW6r1HMVPP/1U5ubmpra9Zs2asuzsbHXKy7yetzvDflZlX7kPp0+fLtu5c2fZokWLyr7++uuy+++/v6x58+bq/6stfA6Wnp8gCIIg1AbUFdQX1Bk8DlF3VIY0T7hw8wSjSIxAMfI1YMCACreVlJTgn3/+wcSJEzFr1iy0bdsWbm5uaiHaeUuLtdtLS0tx3XXXoVWrVnjuuedUOlO7jW+jRx55BIcPH8ZXX30Fnc4QDOb1W7Zswf333485c+agY8eO6jrTZffu3SpF++KLL6rn0bt3bwQGBtbotWF0jZE5RgR//fVX4/4UFhbi5MmTuOGGG3DgwAE89dRTSEtLQ3FxsfG+e/bswffff48OHTqoodoXXnghdu7ciRdeeAF//vknfvzxx1pNHQuCIAiCKdI8UcdMmDABXl5eRpFDEWF6au18TW8/dOiQevwPPvgA7777LvLz841LQUGBMXX48ccfo2nTpup+1Vm0xzx9+jSSkpLQp08fo3gzXaKjo7Fq1Sp8+OGHSvxR+PH6HTt2qP1gjRtFlDkUW2TDhg1qv1euXKnE3ZAhQ6ot8FhTx1o+plG5D0VFRVi9erXaP4q4du3aqcdZu3atqvPz9PRUr1FERITa/x49euCiiy7C3r171XUU0BSIl1xyCWbMmIGLL75Y6vQEQRAEp0Midg5Q0IJ1KH64UFjxtfLz81MiikJYWyjsGCXr2bOnup3Ci9f7+PigdevW6NatmxJXAQEBxoWCTzvP+2gROQ0KOtbU8X/EKOKiRYvUecL94HLTTTepaOa0adMQGhqqtqFFQSn6WAdI4UcY/Zs8eTK2bt2KQYMGYcmSJRgxYkS9vKaCIAhC4yJLInZ1y/PNIuGt06EUZYZirDJDUZZP714InjxZpRv1+QVIfPbZswVbZcZ1Ss9e9u7YEUGXXGJYt6QEya+8UmFbhnUNGjwrMBBHmjdX665btw7XBQTAvax8W6b3SykrQ5sbrjdsV69H6g8/oqSwAPqyMpSUAXqUQX/2tCwgEF49uqs0Lpec7dtRUlSEEqY2y4DislIUFhdD5+GBUk9PKiTjukVZWSgsKUGhnlsqh4+pRQ75huRiDaZszfnvv//s+h/wDR8WFmZcuE9k7NixSgT6+voiPDwcw4cPR79+/XD8+HF1O00feb0Go4qke/fuSjQyCvnHH3/g4MGDKpLH6KzpeoIgCILgTEjEzgEKeteNUxHg5QmUUk1RYRkklv/gIQi/aZpatzQvD8duufWs4jKsQwlWWqJHcVER8tu1xbFhw5CcnIzTyckY9edfcKNEKkOFU52bDnmtW6Fw2jQ0a9YMGzduRMxzzyPYzw+BAQHwcHdHSXEx8nJzUVRYCM9uXdHPpAN0/+Ah0GdknN3Hivj26IEW335jvHxg+AiUJCdbfO5e7dqh9e+/GS8fHHcBis/amJSUlaGwrAzFPC0thT4qCjsvvwz33nsvBg8ejKkFhShJSkJRWZlxKSwrNZz38YHPddepaBnTvmnLlqEoNxc5paXILtEjt1SP/NIyFFBkskavBv8/ppjj4+MRGRlpXJgKnjt3rorknX/++ep2Xvfdd98hIyNDCToKvNqK2Ik1iyAIglCTiJ0IOwe80KzX4sF3+vTparEEI0ipqalKuJkuOTk5FtfXar4o3rSFl5l2tMcfrmXLlnj11VetFvmXUdgxusalRA/oDREud5PUctHx4ygrKkJZSYlauP7SRYvw5OOPo++gwbjumafRpUsX1Vjw0+OPY9/Gjbh3+nT07NIFZYVFKCssRFlRIXQBgQifNtW4n1fk5qE507E6NwT7+qFN8+YI8fdX67uHhqLFN18b9+HIpZNRsHu3xeeQ7eaG/nv3qBpCCsZd996H0zt3IlOvR2apHmf0emTpS5FRqke2zh35bVorschFi+jZAz9EzZs3h4eHB44cOaJSxxSpV1xxhfrf+/v7o7pQNC5YsEBFBo8dO6bq9+ra508QBEFw3pFiXPijf//+/SLs6kNB8+Xk9ZpwS0lJUacUdewqtQRTh5pw00ScVvPljJGe6gjJ6u5nUWoqRvTth25t22DG9Luwc8MGpJ06CU92sbq54Z19+9Tryy7W0y/PQkBKCrwLC+Gn18OzoECJUZJeUoIhhw4a93Vem7YoOHQQp0v0SKXoLinBaX0J0sqA5LIyLMtIV3V+rA2sjNjYWCXwTBcKXopBrcvYGkw1s56PNYY//PADxo8fj8cff9womKUTVxAEQciSiF3dvtAsuM/NzTWKOUZhLMFmANMInCbmKCBcjboSkpYsXRISEtT1jJ6xVu7TTz9V9iWMbjVp0kSlUGnvQkpz81CcnoYNK1bghE5n3Nczn89Dwd49KDl9GiUpp9Vp6dnmCveIJth5331G8fp8ZCSaeHjgVHExTpWU4JTOHUcLC3AsPx/ZJjYp5vC9QbsVLmwA0RbW7hFG/t544w3k5eXhrbfeUu+F22+/HVdeeaWK4Kn9Ly1VnbgUebRnqcprzO1v375dPaa3t3eN/g+CIAhC/SHCro5f6IcffrjCgZORNgoMcxHHIv7KIjgCLHa3suZOE0QaFHj06qMYuuqqq1R6lFYs1RWYpQUFBoGXkwOfjh2N4tXzvfcQdCoRHqxNNEvhJnp44Kfhw5S9S3p6Orb/9jsOJifhYEYGSixMquD/n15+bOBg0wY/pBT2tKv53//+h5CQELX/48aNU8+F61PUVqcT999//1X2LrSOoZ+hIAiC4JpIV2wdw4M6o0VaBI6ijvVYQs1hhI0wYmVuwszGBkajyNSpU9G/f/8aPZbOxwdecXHGyxRYSkidFVOsRyxJSUHxiRMoOnYcOfv24fi+fcqKhRHbSZMm4cbVawxdwmHhOFJYiP1FhTjo6YkDAPZmZSHpzBllyMzF9HEI78ft8IfBX3/9pZ4z31OstyS8D2sJtcuVRVAZ4dNeOwpFS/cTBEEQGhYSsXPhyRONAWsTJGqapnQUFGO///67SnmSkb6+6BgQgKK9+1Cwdy+Kz9qqEN9eveDz+mvKiJl1c+t//Q2HMs8g+6xBsyn8gcAaPf5ooCCjITTTzZyGQTHLKBxFn7WaR4peCl1GBwnn8bJuTxAEQXA9JGInNBgo1tgVSmFCEcexZVpjAUeQaY0F9WUJwsgsI3WM2tEEeUl+PhKiohA58UKkDRqI6JAQ9AoLQ/G+fXAPCUFwZKQaUZZ05Aj+t2w5EBKCE8XFWJOfhxXFJUjw9MDBpCTVdMOFdjba4zD1zNQsa/LWr1+PuLg49SFnJy23ybS19tqwPIANJezcpbjbtm2bCDtBEIRGgETsaoBE7OqOmnTi1hXsbH322Wcxf/58nDlzxqplCUfBfffxx+hw6DB6+vuhaMdOVdenwQ7enZ064T9vL1VDaG6Jw1pNRvLYIMI6Okbw2ExiGs385ZdflLij9x5rECkMOa/XvE5REARBaFh6w34/DcEqnF7QqVMn5TUj1A4URZz+wAaCr7/+Wp0y/eosok4Tdt9++62qceO4stmzZ6v9ZBqZEUeKU8KZtUXe3vCdeiNafvYZ2q1bixY//ojkC8ZhvV4PLzc3LF6+XEUiGQn89KWXMKt7d0zu2ROhQUGqkYQRON7OD/nJkyfx0UcfGadh8LcaBTBn7XK+L2cHM23NSJ4gCILgOlBXUF9QZ9iLROxqgETsBEu1gJxM8cUXX6hGCNbBXXvttbjmmmuUsGKDA0UYo2v33HOP6oI1386KpUuRdOoUIuPiVBNE5ldfIfmFF41TPba6u2NFYAB+378fySaRQcK6OqZlGanjGLXnn39eRRFZC0jBedttt9Xp6yIIgiDUHLE7qSNE2AnW/PYYKfv8889VGpWNEGx6GD16NJ5++mkVPWO9G6N49qA/cwY5q1Yhd/ly5CxdBv3Zebv86O4oyMcHQUE4WVSkXMk1WIvHFDC9/n766Sfs2rVLicY77rjD2HShti0jzARBEJweScUKQh2jpUG1BgWmUG+88UZVD8cmiC1btqjrd+zYoU5pW2IvquliwgREv/wy2q5cgeZzP0PoNdfAMyoS7Xx8kVZUhOuuu06lX1+f8SA6tGypRB9FHWHEkLV6TF1v3rxZCcvDhw/jySefRExMjBKk9ArkKaOOWspYEARBcD2kK1YQaslvj36GFHeM3DFiRjhXlk0f2vpVxc3TE/4DBqil2eOP4c8PP8SOO+5AwRdfqNFmF+7Zg3Fe3ljboiXeSU5CYlAQTiUnq7o+Lr/99puq1eDM4b///luNP2PqllHFtLQ0nDhxQkUSZYSZIAiCayKp2BogqVjBHr89zgmm0TGbHO6++25cf/31al1HwTo6dgcXZGbixagojPAPgPfZx+eH+3CnjvjLzQ1fLliAnNxc4/34nuW0Dk7BYGMK6/I4ju2ll16qljcgR+lxikZVZhwLgiAIlSM1dnWECDvBFKYwGe1iBMyS397NN9+sBN7VV1/t0NFyrOOjsGMqlmPKcpKTcUmzZuiUchp5GzYY1ysdMwY/NglXVijsqtVg9PCCCy5Q72FG8Lp3765SxVUZYcbHZmSSNXqjRo1y2HMTBEEQIAbFglAfMHXJFCb99jjbVYOp19pMbdKbjs0ZGqzrCx03DvEDBqAoIQEZ336HMz//jLgrrsDMoUOUeLv7uutw7TXX4Jt//lH1gZ9++qmKtLGtfubMmRXqBu1h69atqq5v7969IuwEQRDqERF2guBAKN4uvvjiOu807dWrlxJVGjQwJl7x8Wj2fw8h4p674ebtra7jPt0QGoab163HncNHYE2HDpi7Yrnq7GWE8fLLL1frMRXLFHNl+05Bp82lZdq5uLhY5tIKgiDUE1IMIwgOhkJIS7nytC7sQ1izp4XnGcFjV64pOl9fuJ2tfaPQbBIcjBI3N+i3b0e/77/Hx/7++PXxx9G9a1fjfRi5Y3TvrbfeUibH1mDtID37CDtu2QUsCIIg1A8i7AShAcA0KqN2msizVcNHodl+1ssYffAAloaGoszHB0UHD6HdDz9iTlExJgYGYcyYMQgLC1O2KDRSZjr5lVdeUVMvzNGidRpJSUm18AwFQRAEexBhJwgNBEbiJk2apDpb7UkZv//dd5iVlIgBO7bj9dMpSC0pAeN8l3btqjp3jx8/jnfffVcZHScnJ+Ohhx5S5znJwnQWrmaMrEUMRdgJgiDUH9IV64AuFaarGAWZPn26WgTBVagweSI8HK0SjuHbffuQHxigrFn8UlKQuWYt/tKX4KVZs4zROYq4iy66CH369FHmy82bN8f48eOVETL99NgBLAiCINR8ViwXflfzR7TYndQyYnciNES++uorJeCY2u3593zkLF0Kj6gohN52KxYWF+Phxx5T0TxTQkND8frrryvbEzZO0O7FHj87fnZolizNFoIgCNaRkWKCIFSbYcOGqdNtW7fCvX8/eDRtipLERJx+6ml0fv8DtM/IQJ/evVWkWoNfNlOnTlUTNtgVm56eXunjMEo4e/ZsNQ1DEARBcAwi7ARBqEBcXJxqltCXlmJrs2ZovXABmj3yMNzDw+CVloY5MbH4oV07rJk3T82hjYyMVN2w5KeffsKqVatU5K4yOO1C876jGBQEQRBqjgg7QRCsRu02b96M3OJihN1wAxKfeALvpqaizMsL+Zu34NDff6N169bKFuW7775TdXYUav/995+aZPHhhx+eI9hYJ0K/vG+++UatR0FYUlJilxAUBEEQKkeEnSAI58DuVwo1iq7Vq1er605lZODttFRE//oLQqdOxcqAAHX9wIEDcXHfvti8dq267Ovrq4yKb7vtNjXJglE8Cj6OXKMVy8iRIzFlyhQ888wzyiNvz549OHTokOqmnTdvHk6cOFGvz10QBMGVEWEnCMI50AdPi9pt3LgRGRkZamIF2ZuSgvSJFyIzNxf+/v7o1K4dTky/C0mXXY7RAQGYfOmlqmOWJslMt3J+LgXe5MmT0bVrV6xZs0aZGN90001o2rQpvv/+ezVy7ffff1e+ef/++289P3tBEATXRYSdIAgWYZqVkTumU3/55RcMHjxYXX7hhReUOCO0OylNTEJpfj50aWl4KyYWU0+cxNA2bbBo0SI8+eSTKoKnjTtjByxr8nJycpQtyrRp01QTBqN6WqQuISEBaWlp9frcBUEQXBURdoIgWI3aXXLJJfD29saxY8ewdu1avPbaa/jzzz8xa9YsNUqsffv22JychIf9/fBBWhpK3d0RlZSEcfP/wclZr+DJhx7Cxx9/bNwea/F4n8cffxwFBQVKKDKSR8NjPoaHh2F8Nb3xBEEQhKojwk4QBKuEhISoRgiyZMkSJcYuv/xylUqlYIuJicGgQYOwZfduDPvgfbT960/4DBoE99JSxG/YgD1jzseaP/9U97/lllvQu3dvFBUV4euvv8acOXNU5G/AgAHqdl4/ceJEdX7r1q1qfJnWbSsIgiDYhxgU1wAxKBYaA/yK+OGHH7B7925j5K179+4qrcqUKWvvOM6M01e09Te8/Q50n32G7MBAfN4iHp/Pm6dq6igEmYalCTLHlBE2aTBa9/nnn+Oaa67BG2+8odYhnFd73XXXKQNkQRCExkpWFfSGCLsaIMJOaCzk5eWpNCrTsqNHj0azZs1srs9u2nffeAO5aWnwj45WKdwOzZvj6nbtcTQ2BqVlZcpKZd26daqDlrDhgqKOjRq0QmF0kF9PrMmjKJTpFIIgNFZE2NURIuwEwTpM19LKpG/fvmqG7PY7p+PKkBAcDQ3Fpr59kFxWprpmeRvHjzHt6uXlhRkzZqiRZLz8wQcfIDc3V6VwtTStIAhCYyNLRopVn3379qFHjx7GhR19v/76a33vliC4HLQyGT58uErZTpo0CSOvvAJFZWVokZGBC/5ZgIwvv8L+fftUR+yOHTtw3nnnqTo7dt2ywYJNGmzeIBxVJr9BBUEQKkcidjZgnQ+79mi/QL8ucyRiJwhVI//gQey/9z54HTyoLvsNGIDoF1+AZ1SUEm6cG3v//ffjyJEj6nY2ZtD7junYe++9VzVVcMaseV2fIAhCQyZLUrGOgZ17PNCwtsgSIuwEoeqUlZYi45tvkPLKqygrKIAuIACxb78F/7PdsaytY00eI3es7WOzBmfX5ufnK1GnwR9dXO/SSy+tx2cjCIJQ+zToVOzy5ctVrU10dLT6wreUJn3nnXfUl76Pjw/69++P9evXV+ux6Ih/5ZVXOmCvBUHQcNPpEHbNNWj16y/w7d4dbl5e8G7b1ng7P7ePPfaYKou4+uqrVSSPEylYs3ffffep5grapDCSx6kWHFUmCIIguGjEbv78+Vi1apUqpuYvdTria3U4hNG166+/Hu+//74SdbNnz1ZWDTxIsOaHsHaOXXvmLFy4UAlGTR1zriVtGHigsYRE7AShZpSVlKDo2DF4t2plvC5/xw74du2qzuv1evW5ZQNFYWGhuq5Lly5qxixHnvGzv3PnTtWkQdi1a/WxysqUqTLTuJLCFQTBlWg0qVhG7MyFHcUcu/DefvttdZmddXFxcbj77rvx8MMP273tL774AgsWLMCXX35Z6Qt9/PjxCi80Dy62DjCCIFgma/58nLzvfgRffDGaPfE4VmzciJEjR2Lq1KlK5LGhIj09Xa3LaDojevz8P/HEE+oH2PTp09Vn0hKsz2OEv2fPnrj44ouN1/Pzy/FnY8eONc7DFQRBqE/4Q1b7MavpDWqZBpmKtQU76jZt2qS66zRoo8DL2mzL2kjD8sXmwURbXnzxxSrvuyAIQFHCMX5okfnbbzhyySRkbtigrmfzBCN1+/fvxx133KE+14zOT5kyRd3OqB0//0ePHrW67aVLl1ocV7Z69Wp1P3bmCoIgOAPUEaa6gjrDXhqUsKPRKX/Vm5un8nJSUpLd26EiZl0ef8HbA3/x8z7aQg8uQRCqTpPbb0P8l1/AMzoaxSdOoO1nc3FzWBhSkpNVOpZRuXfffVf9gBs8eLBqriA0NKbo06ZZWIJTLDRYp6dF9LUOXDZnCIIgOAPUEaa6gjqjUQo7R0F1zAMEzVLtgWFR00XSsIJQffx69ULLX39B0PgL4FZaivsjmuL81avhnZdnnFLRqVMnPPTQQ2ocGaN3nCv7zTffqOsOHDhQYXv8scdoHRuvGJmjmGMzBmGXLbtwnV3Y0XqJtcUUt4IgNHy8vb3P0Rb24oEGRJMmTVRRtPmvdl5mKqe2YE0fH5f1PVwEQagZ7kFBiH7tNfgPHoyTTz2NrsXF+OW339C0fXucf/75qimK9XZsbqL5MT/ja9euxbZt29C5c2dce+21mDNnjmqIeuCBByqkaENCQlSHLb8TNIFnLuwo/ljDy8UZoGhl4wcnddxwww31vTuCINQRdPngwh+ojVLYMcLGblkWQmsNFfyC5uW77rqr1h53w4YN0hUrCA6Goipk8mT49uyFnx9/DP/Mn49vb77ZeHtoaCiuuOIKdOzYUV1mUwTHkx06dAifffaZaqw6c+YMLrzwQtUM9ffffytBt2LFCjWqjMXII0aMMG5PS+syOvbee++hdevWTuORR1FHtLSxIAiNg+lnA0Zas2aDTMXyS5fdbVy0Lzqe5y93Qtf6jz76CJ9//rmyQGChNdMX7KoTBMH18G7VEpd+9hlef/11Fa2aeuFE/DNoMPatWKGmTxCmY1u1aoVrrrlGNT0xKkdRR5im1X740d+S4o9RPVonnThx4pyIHetr+Z2xffv2OnuOjDjyO0uro+G+V+UXuiAIgssKu40bN6pf5lw0IcfzTz75pLrML/VXX31VXaZfHUUfv8jNGyocnYplzQ/DpYIg1E69yXXXXac+37e5Ac3T0pA+dRqGeBrqYGNjY9VnnFG+Dh06YMKECcb7LVu2TM2sZcTOw8NDTbGYNm2aEk/8Qejp6WkUdnR/0hortIh/XcCuXv5IpZhk2pip5nnz5p2zHgWsIAiNh3feeUfpC+qMBpuKZeqkMus9pl1rM/VqjqRiBaH2oSi74IILUNy9u/K6y9+2Dd5vvIHLLpyAyAkTsGL1auO6WgMTo2A//vijWvg5ZRQ/PDxcmRtr0bx27dph165dKkJGyxQt0qelZwMCAmr9uRUXF6tTPj73k3BGtTlirCwIjYvpjSEVKwhC44ZWKPFfzEPYjTeqy+5//oWc/92LQe3aISIiQtXF8RcuYUSOk2feeOMN1VzFUo5bb71VNVeQwMBAtG/fXolGTciZWiNxfXvgj03W8dH/sibCjqe2UrAi7ATB8Hlj2YL2uREqIsJOEASXg/Nlmz38f4h5aw50gYHI37oVmTfdjJtGj0a3bt2UKTnr7L766iuVTmVU7vbbb1cLo+scMUjYaMH1fH19jWkP04OFvfYiTN9yW7t377Z5sKFooyeVOYzU2SPsNAEqCI0ZdsB/8skn+Prrr+t7V5wSEXYOQGrsBKF+CBozBi1//gk+nTrBIzISXi1aqOs5Goy2KKxZY4c805ucD92nTx8MGDDAeP+VK1eia9euajINBZX5DGl7I3amtXi2hB3tVxg9NJ+QYRqxs1XXJxE7QYAa90lYl8oou9Yg1RB5pzHU2DkjUmMnCPWHV1wc4r/5GvozZ6Dz8VHXhYWGokf79sDll6tf93/88YdxfTZP/PTTT4iOjlYzpNmQxctM444bNw5jxoxREy44YozWKIzCUSiyA9easDIVhIy++fn5WVxv3bp1RoHHlLClGjtbxugi7ITGjmZLpDFr1izjpIaGOBxgejVq7ETYCYLg8ui8vaEz6XzP+OgjjF2wED4DB2DqBx/g22+/VQO1b7nlFtUhqwkkCq1PP/1UddefPn1a1cnReoSp2y+//LJCIwWtUl577TWL3nZaKtVWxM7UANlcoJlG7GylWx0l7CoTkILgrGiTYszhZ7U23S9cCUnFCoLQoCjNz0fmL7/CNzMT5/37H05+970SZZdddhlGjRpVQRzRPuTmm29WUyr69eun7FIYTXv77bfV+ZkzZ6o0D1O1TNlyGz///PM5j2kq5kxFnikcX6bBX9+W7l8XNXaMRL7wwgvG6KEguBJpaWkWr3eWKTHOgAg7ByA1doLgPOh8fdHi++9Q1LkTPPR6DFqzBj22bEG71q2t3odiilYqjOgxDas1RNAUmVMsONHm119/VVMsZsyYodZnBI5lGDw1FXOWhB3r5kyHeLOBghFES6nY2q6xY9qZ0KBZEFwNWhc1JmH3TjVq7ETYOQB+ubMOR+bECoJz4BEaitDXXsPus+PG2u/bD/9XX0WJlV/7rK8jFFtM9fzf//2fskdhtO6ee+5Bly5dlLBjVI8F26y948iyv/76S9Xv2UrFssbvxRdfVPexllKyFbEz9e2UGjuhsWP6g6gxCLvp06crfaH5W9qD1NgJgtAgad6iBdzffgtZCxbA7a23UbhxE45OmYLWf/yh7FJMGT16tKo5i4yMVKbGt912m7q8efNmLFmyBPv378fkyZMRFxdnTKvyOsIvXTZkWIvYcfKNJUwFnC1hZ3pZhJ0gWEamspQjwk4QhAZLTEwMYqZNQ+GwYThx190Iv/mmc0Qd8ff3Vx2xS5cuVZdpUkyxx4MFa+sYNaNNiZZO/fjjj9X19MAzj9IxNcvtUPyxm5apVe2go3XLsrPPkrCzFJEwvU187AQBjSpiVx3kW0IQhAaPd5s2aPnrL0Y7FFKcmAiPiAi4mYglWpqw0YKpU6Ze6XNH/ve//6mo3EMPPaQaHxYvXqzEG2vvOJ5Mi9JxZNmHH35YoVGC4o+eeh07dlRGyNq6tEhJT09Xj2Eq3sy7/lzJXZ8CWA6wQn0g77tyJHYpCEKjwFTUlWRkIOHa63D81tugN5kEwVQnLU3+/PNPZWyckpKiImjLly/HBx98oEQdfe5atWqlInGse3nzzTfx/vvvY+vWrWqkGAUgu2i1btqmTZuq6yn62JihpVN5Hbtxme4l3B6Ni7ds2aJOeZn3MZ1+Udmc7PqEZs587TTzWHtw5ucjOB+23i+asMvNzUVqaioaMxKxcwDsVuGXtWYkKAiCc1O4bz9K0tNRfPIkjl5xJWLfexferVqp2+hTx867Bx54oIKxcVhYGK644goVeRs0aJBqomDk7sSJE/juu+/UOuHh4arxQptuwdOrrrpK+ejRRmX8+PHGiBy7bjUo4Hi7qW8eI30cfWbaDefMQmj16tVK3FHMjh07ttL1T506hblz56qUd//+/etkHwXXxtI4Pg2ttOGVV15Rp/fdd5/dhr7O3hXLxZYNkjkSsXMA0hUrCK6F/4D+aPHN1/CIjkJRQgKOXnkVcky6VinuDh48iEcffVSdv+GGG9Tnm6JOg5G5adOm4eqrr1aCTvPY4rqMxGm1cowkDBkyRIm2w4cPn9MAQVHH6B0jezfddJNy0OepFukz7YazZYVCGCWsL/FX2b6Zw1Q309JiuyI44j3Gz5xpo9KpU6fQWLtiRdgJgtAo8enQAS1/+AG+fXqjNDsbx2+7HWmfzTUKIwowzpZlkwTr7ix13VG0cUoF6+xIaGioEnBaTR6NjllLR5GmpYlMhR0PVIzUcRuM7MXGxqpuXJ7yMq/n7doBzZZoW79+vUqFag0gzi7sBKGqmEa0LUFrIUGEnSAIjRiP8HDEf/opQi6/jMoEKS+/jPTPPzfebm3mqznarGjW5dHEODAwUDVQcBYtIwma2OLII1Nhd+zYMXWwYtOGefG3aaSP61Umnv7++291umzZMtQVjFAyXc20sgg7oTah7dC8efPqezdcAqmxEwShUUP7k8hnnoF3u/bI+OYbhFxyiV3CjmKM6RHWzDVv3lzVxK1atUpF2rp3766aIlh3xrogRhIo1FiXR9FnmjolWkTPHO16bT1bETtGFM3FFaOFFF2sNTp06JBq+nAkrJHjvrHZQxBqk7r8weLqiLATBKHRQ9EVdt21CLnyCuhMfO4C8vOt3ocpU9bfff3110pU0dKENXFslGCkjeKOwow1ZOyupSh77733VHSPEyzYFKGJPN7O9Ks5vJ5o69mKilkSdoymbdu2zXiZqV1HoglOa/M769qegj6DfL05Hk4zk64qnBBCn0F7GkAEwRmRVKwDkFmxgtAwMBV1Gd9+C58ZDyLm+Al1+eabb1Y1dBqenp5KKN1xxx3qMhsr2DVLMfbpp5/ipZdeUlMs2CDAqRUTJ05Uhsm0TPnvv/8we/ZslWKl0KOYMI/G8TI97hgJZESQULjx/tq69LhjVE7tu4UaQFNRR7RpGc5AbTR5fPLJJ6ponqfVZdGiRaqztzpiVXAclc1Nbiy8U41ZsRKxcwBMx2g1NoIguD4UHTnLV/DogiGrVmFb927K7oRiTkM7r3XEauKuffv2SrAxmsVIG0WZJrpmzpypplYwskdzYtYNseaOguuLL77AyJEjVZSP4pCijtdTLGr3p2h5/fXX0bNnT1XL9/LLL6vI4YMPPigjlRyMKxlDNzQYMZ01axaio6Nx6623ojEz/ayNGn/Q2WvfIsJOEATBQpowds6bOP7008j94Ud037YdZ55/Hl4m6T0KKm3MF9OuWnSMAotdtJdffjl++OGHCtulGBw1ahQiIiKUhcHGjRuRkJCgbmP9HRcNRuo03zxzaGI8ePBgFa3jwjo+04kVp0+friA4rT3HmsD6QWeCY+AE14R1oD///LN6T3fo0MEYWeYPGb7Xf/vtt2r9OCtrpJNQ5CeeIAiCBThqrPmzz8Kb/pQ6HTJ/+hldfv4ZXmf96Uyjd5MmTTKmZDUo+NgFaz6TlhE6ir8uXbrg1VdfxXXXXYfWrVtXWI8zZmlmzIOcNUzTVEzrmqdvOAnDFlqErzrpLjr704alJjj6gMtIqCNxZjPohgaFG+sjWZ9KTM14qyPqCKPiTMk3xv+jCDtBEAQbtLr7LsS99y50/v4IPnYcoxctgk6vryDsiPllCjtG7VhXx6gbp1CwK9XU7oRCj9dde+21uO2225RnHgUPrVLYlMExZuyutZQW5G222LFjh83b+ThMeVFc0iy4KpiOOXMWtFpDwfXIN2lS4nvdUe+vEydOVNh2Y0FSsYIgCJUQMHw44r/+Gvuuvx5HWsSj1N3dopAzv9ykSRPccsstFa43FXamdiqRkZGqy5apWtqjUNAlJyerzlY2W/Tu3VsVUGv1vJWNGKJorCxix8eguGN0jx589mJvlI91QRSQphYvGrYiKUwl83l6e3vbvU9Cw+Czzz5rMFMj6gsRdoIgCHbg074dDt99F/YeOGCssSstKIDOx0ddNhd65qPDLF3v6+t7jl0Ja+vGjRuH4cOHq/oiTpRgDR2bKeiTxw45zlalPYqWzuT9zBs2du7cafP58D4UddXBnvQWhScbPcjjjz9+jvC1BlNyTKFR2N1///3V2j/BdakNUVdaWmosPeCoP0aoO3furMohGiIi7ARBEOzEndGys2LKvahIzZgNGDECEf+7p9IIXmURO0anTNNGFH2DBg1SKdx9+/YpgUcj4F27dqmFdXgUeNweLTpMxy1RHNJXz1LjhWnqsroNEPZE7GhXoUHrEPN6Q9MaO9MidxbRa9G+uqYx1mM5A7XZ4JCZmYk33ngDvXr1Uv6G/IHE+cxcRNgJVmF6hF+uWluyIAgNE1Pxpl+zBoX79qml6OhRRL34gjpAaeLAHmFnGrFjBNBSPRAjDRRoXNj5SYHH+jnW4Wm1cfTXYwctmzBolUJfPBaPW+uqtYeTJ08iICDAosWCNuLMXig6zYWdKazzY10V7S3YIekMiMirXRiZ/eeff1R0ujZZtWqVqttbt26d+pzQnsyVYCMUl8pKL0wRYecAxMdOEBoHpl5xwRddBB9vbyQ+8SSyFyxA8alT8G/fDjlnBZ09ws5UKJo3VVgqIGcd3kUXXYTzzjtPWaVwBi0FCMUQhVybNm3Qp08fJeh4mZ2r9NWzx+PONF3FGrePPvpInX/qqafOWXf58uVVEkaVRfi05yq1VY0HzUSap3xf1wX//PMPXAV+ZihIq+NjJ12xgiAIdmIqVhhh41zZ+E8/gXtwMAp27MDIfxYgMDPLbmGneeER085X88YH1gOZwhQu6+i4P6NHjzbapRw8eFBZRsyZM0dF2xgpsze6ZtpVymhKTZGIl2AO34/yvoBd0TcK3hdffNE4tq8qiLATBEGwE9ODkhbd8uvbFy2++xae8c3hl5OD8/77D+GpqVWO2Jl2gJo3XpjW4pH77rvP+IXPubO0S7n77rsxcOBAld7lr3t2vBJ21NL4uLIDqibs2JFrKyLnCNHnCjR0AUJD67p8juy8pt+iJWud+qinrC8YNXzuueeUF2RlpRBk7969VX4MEXaCIAg1LPj2atECLb75BhlNm6LEwwN5fn52dcWaCjvWwrHAm7Yj5qlTn7Odtxq8XbMQYU0d4cgzNkywk5SGyRxNph0g5s2bp+p0KNqs+XoxYsj0Dw88po0YpvB2e+ewmooGnmcRO2ewmk7IqIzt27cr42HzAz8bMz788EMsW7YMDRF2bnKpLZhq5yxjeiXWFSwbsDS/mFS3O9sedlbSHV7X8DPIz0NtvndF2AmCINiJrQiHR1gY1l84AUtGjUS+n59d9h5MxWrrtW3bVtXP9ejRo4JoZNOBefcehR1Tsex+ZaOE6X5xezQ65m0Uf/S/4+OwM3XBggXKgoRRE0beTO/H2yurheMYNHsxF3affvqpenwKisOHDyuPvspghyxNZs1rozZt2qTq8Thrt7pQaJr/P50hSsf/AVNwXKozFcQetAaCA2ete4TqwR8p/IHBBg1nQponBEEQHESZpydyzkbSKM6y/lmA3DVrEPnE42pEmXl9DSN2TKump6cjzmQOranAmDZtGnJycio8DrdNcccIHZskWFc3ZMgQFaVjBI+WDpy3qXXFjhkzRnXSsuGCgopREy5cv2XLluqUQu/tt9+2+fyq0g1rKWKnwQhiVaAAMbVEMa0H5PNntNIWjPCxm5ivBWfoUqDyfhTMl112WYXGkfrGNKLK85UZTQv1x5o1a9QPDC6cc+ssiLATBEGwk8oiOqaRtpL0dJx65BGU5eerjtmYN96Ae4B/hSgM07I8cJsfvE3XYQTOXHTwMv3tCMUbu18ZEdNgtM7U6oT1e+yWZfSOqVkKPKaoKAK1VC75+++/MWzYMKsWKREREVYNkU338d1338XQoUPhKJgmZo1Wz549z/k/WIsi0utP+59QxGn1hhTKjHISvgacDkIhzAkh2vOzBlPUfC1Zx0iB/MMPP6hJITSNFizD/xUjrNbS+86MXq9X3pF8f5v7VBJLo/6cARF2giAIdkJxxDQW06aVCTumZmNeexUn738AuStWIOG66xD3/nt2pddM16FgMjdw5XVjx45VUTmKMFqa2BJapvvHiRVcjxE7ChlGs7RoGn3yKIL4PBnlM+3a1YQojV0pJG0ZIlMs/vTTTw5NcbIZRBN29kDRZQkKYFN7Da32699//8WUKVOM1y9evFjVNjKqx9eMsPhfs4Dh9lkAz9fLkiWMMxv31tZ2tf8zm3UYbWQ3N821//zzzwrrMULtCixcuFD533Xo0AFXXXUVXAURdoIgCHbCereHHnqogrGwrQNm4KhRiP9iHo7ffgcK9+xRkyrcbr6p0sehHx0jQiNGjDCmXc0fhwujTTxI8vYWLVrY9RwoGnnAateunfFgRVHICBfFE6MUjOgxddutWzeVqtRq/DjhgkKG9508ebIx9VuZIbIjTIcdKUasiWvaxWgcOnRInfJ1YZSOz9faZA1bsHbxt99+UxFM7QcBIz2WIkDOUudXXdiMweerCbeYmJgKEWEN2vG4AuvWrat2Z2p94hxFBYIgCC4CrUeqIjJ8u3ZVdiherVqhJCkJYa++hmZJSTbvQwPiBx98UAk7YiliR6x13hLOmjVFq+GjiGO0jUJDE4jx8fEYP368qvdjtI6wK5PRSTZiMO3LujiOZqKoY4SKUSxG9HhKgcjrKRgtiSZGvxxJTcWPpftTgFjrEuVrUd3HZAMIX/OvvvpKXaZgfv75542ioS5h92ltdqCyFtI0GmfJZFuofQEvws5BI8VYY0E7AUEQGi/WBJ9XbCxafPO18rzTFRQgMinJZlE8t2N6u6WIXWXCzrwrV0urav53mh2KKXxMikrCYnBG3/jYFCE33HCDGmPG7Zh3tHJ/2LxRFUPk+sSS+KzMS838QGyPuGf001zcaCnq+fPnoy5hw8msWbPU4shuW0ZjtUivOUzba6luoSL8HNI6iJY+tmDKn/qCOsNeJBXrAJ757hn4BRoMRP849Ifx+rjAOPRo2kOdLy4txsKjC61uI8o/Cr2a9TJ+gfxz1ProkwjfCPSJNPyqJtxuaZnJB9Xk+ybcJxx9I8vfEIuPLVb7Ur5q+crB3sHoH9XfeHn5ieUo1BdaXNff0x8DowcaL68+tRr5xWb+WGdX93X3xaCYQcarNyRtQHaRZTdtL3cvDIkZYry8JWULMgvLu+lMcXdzx9DY8gLtHad3IKPQespnWOww4/ndabuRlp9mdd1B0YPgrjMcNPdn7MfpvNNW1+0X2Q+e7oa0yuEzh5GUZz0a06tpL/h4GDzJErISkJibaHXdbk26wc/T8L46kX0Cp3Ksj1vqFN4JAV4B6nxiTiJO5Jyo8H/TuRnqtPjXJqSNcd30gnSk5KUY11Hr8Y/pPzedeq9p+5BXnIesoqxz1tHu6+fhZ3wdSkpL1GL6uPx/1eawb2fA1vPjdIq4Tz7Gme9/QOjoUbggNLTa27VH2FkTeqb+d1rtmCla6ozpYKZ3eQDasmWLEneM9rDhgJE93sboHuuPuB+aUKyOU749OPK9U51oSVXv88cffyjBY+/9KIIoAvl61oYth6mxdGVpZNbHWSs3MOfNN980PoY5bEpxBfR6vXpf8z1t79guW80T2v/c1ns2ISFBnfL/wpIHUwHOukQNfkbZIFSVkWIi7BzA02ufhrvvuV+wl7S5pFzY6Yvx8IqHrW5jbIux5cIOZXho+UM2BYqpsHt05aMVBJi56DAVdjNXz8SZQsvdSV2bdMXXE8pTEc+ufRZJuZZFCsXBLxf/Yrz80vqXcCTziMV1YwJi8M/kcqH66sZXlbCyRJhPGJZdWW7cOHvTbGxOMTjom+Pr4Yv116w3Xn5n2ztYddKynxCFxfYbyn8ZfbzjY/yb8C+sseGaDUZhN3fnXPxxuFywm7P8yuUIdTccpL/a8xW+32/owLPEwskLERUQpc7/sO8HfL77c6vr/nbxb2gV0kqd//Xgr/hg+wdW1/12wrfo3MQwdurvI39j9mZDkbclPh37qfE98c+Rf/Di+hetrvvu6HeN4nnB0QV4cvWTVtd9bfhrOL/F+er8fwn/4cHlD56zjoebh3pdZw6ciYmtJ6rr1iWuwxOrnlDCz0PnoU65Dk89dZ64scuNGBM/Rq27L30f3trylnEdru/t7q0W/igYGTfS+Nwo3BcdW6Su99J5GdfRTvljqpm/YTA9fxjpy/Tq8WpTfOi8vBB27TXGy6VFRUj//HOE33AD3MwaFWoq7Mzvowk7U/87plBN1+MBiQdj3s71NCHITlkeYL744guV0qW3HLsFuTDKR++96Oho4/ravlVlcHldUlv+cESzZWEnaFX47rvv1GllnbnV4a233qoQObT1PmWEjQtNrrt3737O7RQibLxhI4tpR7ImVOoSvr/4v7RWr2jv/4q2JeyYJmysufHGG88xBaev4O23345Qkx9kmh+gti12nNNO58svv1TimbWpnAxjr0gm9GY09cWrTqOJCDsH0D+yP7z8vSyKHw1GLgZEGewJLNE2pO0527RGu9B2FS73bta7QhTO1ro9Inogu9jyL+pWwQYRodElvAui/Q1f1pbEmikdwjog2Mvyr4kIv4hz9sn0AEohqxHkFVRh3dYhrVXkx9K6PDib0iKoBTILMi2uaxptJLEBsSrKZU+ahUKsY5hl+wftf6tBodA+tL3VdSlETF+XtqGWuyuJFv3SBK/p+8kc09cixDukwv+SrwOfH8ULFx93nwrimFE5Xm9cj3+MAJdV3F++Jvy/VVjHhAoCweS1N6WkrAQl+vL/J8ktzrUZucwoKI/CMsK47IR1x/Zmfs2Mwu549nH148Qad/a4E3d0v0Od54+SS367RD0/RigZfVTL2fMTWk3ApLaT1Lo5RTn4dt+36r0a5B2kTvne52mhW6F67ubvN2skPTkTmb/+irw1axEzZ46yQ7GENY81Fqrbiybs7PW/M39MRjPoF8eDFMtOeF9GpOixpx2I+BiM6tW2L1xN65OqY71hz2Nq/ngTJkywe7tsrNDGRxHTGjhHRSlt1bqxg5VWMhRrFPTaa8NOVkvCTps8wgab+oL/C9ZtUpDxPEXmuHHj7DIFN//8sAFJa5TROsO53ZEjR55TZ8nopLUOaH5uvvnmm3NEGssW+HmyF82mpyaIsHMAb456E0FBFQWJOUy/fXT+R3Ztj0Lh47Ef2/34H4yxHskx563Rb9m97hsj37B73VnDZtm97rODrR9szXlyoPUIkTkP97MeETXn/j73273u3T3vVos93NrtVrXYww2db1CLPUzpOEUt9jC53WS12APFiiZYKoMRaC6m8EuVfxR5pgL3vPjzsHbKWqNg5IGeETEuFOoUnxqMPjPiqERfaQn0pfoK59uFlf84oWB9etDTxlQvf9BwKdIXqag1o84agV6BGN18tLqet2vrqPOlRQj1Lv/lzTQz4baY+jdP//dsVm6zwdT1m5sN6adz8Adah7VGj3RDpJ6p6ydWPoFw33DD4hOOJr5NjOf9x4yE24IFyF29GgnXX4fmH3wADwsRG2sHeNOIHQWXqbmttYgdYe2cPf53plCoTZw4EZ9//jk+++wz1VTRv39/FZnimCQKEqaRWBPEqB3Ho3Gp7LvRXuo7jW+PsNP88v766y+7tsmJGkxz1yf8f1oTvLzMjlD+H82tb+qzm5evGw2nNfge5HuP4/iqAveZ71eKW1NMDbDtxZqfIqPa9sBUNusuHeH3J8JOEIRqo9XPmYo6wsiXvWlNRrq0NHJlMCJ6adtL7VqX0d7ZI62npE1h9HblVSuRX5KvRF5eSV6FU27L9EcaBS7FH4WbWgoNp7y/V6mXigIQ1mYuPm69I/Sq9lfh/nmf4/htt6Nw9x5snXQB9j5+OYLbdlKpYkaLm/o1tRr9YpSCZsCElhzW5sAS82hGVfzvNJhaovBjqs40DUVByC5c+uEx+sPtcRbm8uXLVc0YBSC3Xd/izJFU57mY17ZREDuz3QnfW9xnpgMDAgIcEk2qDqaG2Kw100TdxRdfrH7QMI3NiCN989pa8Zi0BBtZGKFzBJZm4JpGYCnYaABuDUYgrW2DkUV2U9uLCDtBEBo9rNdj8xCXyogOiLYadS4oKsC+A/vQpoUhbR7qE4rH+z+OtII0VfOXmp9qPM9TRu+UHcq33+DQtBvhfyIR7R6ei5cud8eB2LN1dG7u6jFDg0LRNstw0GK0kqKxWYShRpBUNjTeUpqqKv53muigIGRtHqMmNOg1F4RMYTHKw1Qd667YFMCFtUsUeBSHVU2ZaZEPHtwYKawP8WP6mBSw5ulNe/bJUSKiLmDUytRw+PRp6w1ktYklQ2xCAcfaTkI7HqZPaTLdpk0bu0W3tf8HBaS1sXfVKTNgVJyvpWkDizm2InWskazs822KCDtBEAQH4ePlg+6du1eojbyyw5VW12eqmXg1b47gue/i2G23I/hQMh75TYfnHojB8aJklR5mvWCwW7noZEc1awLZCOIT7YPgomCElIQgCEHqvLfe+5yDW3ULzC01HVCY/e9//zOmHk3hbRRvXFi3x25aWjrwIMp6sgULFigxSYHGaF5VDpKMsGhmyXWNaXqOfn6m0AbGVjSmulE+0/WZ6mONI+sheZC3VpBP0cn0d02jo84QOaSoMzXEprUKfQH5nqFnHm/nDw0297Dek+Jz//79KhJdE2xZkHCe8j333FOl7dW1BZAIO0EQhHpC67wm0bEdEPnD3zj18CNofs01+K1/P1WjyJo+CrkfPikfkZWcl6waW1gzWOhdiEzvTBxD+cGjS3oXjHMbp86XuJWgSFekIhw19RTThB0PrPYc+ClCWJc3evRoFeGjpQRriRjR40KxyYMyGzjshamt+hAdv/76q9XbPvjA/jrnyrBkIMxIjyaiKYoZvaTBNNfl+auvvlql/1lrRpsV2tBceOGFDtun+ojWmU9IoVDljwIycOBA9PI5jrsj/0FE2Sqs8umLlX36qCaeNWvW1FjY2aI2xqEx8u1IRNgJgiA4CTo/P8TOKW/MYO1iWHoxmsX2w/yS+RV8FmnJQ8/CZ99/Fplemcj0zFSnOZ45CCg2+BSS0z6nsTpyNTas2gDPpp4IKwwzLh5lVTsEaPYlVY0GMSXLNCrFJe1SGFWhZQp9wDimjHV5nLLB6Etl3nyvvPIKXBV7XjfzuarEdCyXVoxPuxqN33//HdOmTVOvJWEa/NSpUyp6d8EFF1RrP6qKo61ttAkpjNRxf3mZ7xm+P67uH4k7AnZAp55GNiZgMRL7DgabsxMSEpT4YnTTVWAUkALWUYiwEwRBcFIKDxzA0WuuVTNn3bw8UWaStmTELi4oDqOaj1LRr9atWyvbBkboNANpkueRp+r0WN/Hrt1Ef4O1jFuZmxJ33dK7qdOqROyqIgwsRV4YbaLw0NK0LA5nmpUpRNbh9e7dWzWDCIY5pZVFWpkmppgxFVd8fcmPP/5ocyauo3B0FNV8QgojcaRHty64PnC1Gpv12ZYihMS0Rcum/sgO7oA2bQ6o58ZOY0aJawtantCvzpHwx46jkJFigiAITkrB3n0ozc1VXndDl6+AhwWne1o8MN2pDalnFM69rDzq1Tq7NdZMWYPPx32OrmldEZsTC98SX5S5lSHNJw0epeW/7xN9E7E7ZDfSvdItehGapmLtxXw2rQYNjWmCS7FHaAjLjkcWwLOGjWKwsjFfjQH6BVY2zYMibs6cOXYX2NfGdBBHCzvTCSkU/hS45K6B/ghELpKL/XH7XwX4La8XfsVYZCJI2bIQRoBr04S6qh6SdY1E7ARBEJyU4IkXQhcYgJP33Y+opCSMXLwEK4aVj9HTBBEjXOaYiigaUXOyTbssQ7qHoi3XI1cJu8BiwwGUHA84rpY9oXtUA0azvGaIzI9Es/xmysbF0rbN98V8tJSt2bRE68odO3asOhgzMsOaI54yosdmCdZU8f5VsWZxRhqS3Ys55l5wNcV0QgqtVgijvhdHGGpJX1znAb/A8gkp2u2+vr7KNJsp61atKpruNxZE2FmAvxY//vhj9QuEA7HpNt2QP5CCIDgvgSNGIP7zudh73fUIy8jA6P8WofDIEXi3bGnzfvT0ouEpU7TmMFUbUBKAgJzyWjwSnReNUrdSJPsmo9C9EMcCj6mFadsmBU0wOHmwigby+9BShIbRQ0aYKuy/nbNpeRCnyOMUAXY8Utjx4MxULRd225p2pnJ9TtCwZKbcUGjMxx3TCSna68CmkCdO6BCRtARvLz6OyZeXT0iJRhIGuG/B6o4xWLz5oGrUaazCzrV+7tQB7P5hISO7i1jsq7mqC4Ig1Be+3bph0XmjkR0QgIDcXCRcPQWFJmOQLMEox2OPPYZrr732nNs439USsbmxGJAyABMTJmJY4jC0O9MOQUVBKm1boisxpnh5oN14ZiPy3A0dnJyHOXXqVItCxDTyYi4GLc2m5TYYebnhhhtwyy23qGYLook6zlJlJy1PedCn5YUtGAXUPPB4WtspOls0ZqFWHSjaO3XqZHzffP3113jnky8xe1WuEnWmor4jDqIb9uLaLobXmO+L6kyQaAhIxM4CfDNo6QR2bVlLIQiCINQVOYGBStwNW74csW3awvOs4LGFuXcdvedoncED5S+//GL1fjroEFEQoZauGV1V2paWKRoFKMDsw7Ohb65XkbxubbshPCocp46eOndb1ZxNS2hqzNRry5Yt1X00rzIujATS4oM2KrS3sHR/S+a2jSHS11Dg1AstxcvJJmxYsJaG3422GIoNuDo+Ff8LCEB2To5qJqpN6xNnxeUidhxRw0JhFt7y148lbyEOqGZIn/Ue7LAynSlXGfwVOGPGDPXG4WMwFWsplSEIglDXFPr4YMnIkYh9523ozs7urErRemhoKLp166YWfrfZi3+JP0KLymfr5pTmoI2fYbpGqk8qZq6eiZHfj8TsA7NVGte88UKbTUsxRxf+l156SZ3ysrXZtKaNF6NGjVLj0+677z7V7choJAUfa/HYYMGOWlp7WDK3pSC86aab8Mgjj6hTXrYn0lcb1Nf0BleFTRAck8f37YPDAjG76y5c0KLYoohPRFOkIxg+Oj2Gdo5R1zEd2xhxuYgdx7h0795defZceum5MyM5M+7+++/H+++/r0Td7NmzVVHuvn37jJE3jiGxFKLlLzsWXtJHiCF7nqcHEMUkQ/+CIAj1TYmnJ9zP1q2RlJdnQefniyZ33233NvijmJGz//77r1r7EFoWiqc6PIV5v8xTzRa5LXNxJPMIVqSuACKBHqk9VDduTWfTmjde8DuZ+81mCqZWV69ercQSpzJQqLHRYtCgQWp9c4sVwho/XmbkkLdbi/TVFraipM46/aG+0BpptJFh3XX7EYsk1RFrGTfsRRsMwiZc1dUTf6+DOu4z6ud19kdQY8HlhB2FliXDRY3XX39d1WWw3oNQ4P3111/q1+HDDz9s/BVgjR9++EHNmtPMDSdMmKBq7ETYCYLgbORt2YL0uXPV+ZLTpw1NDQ4QKhREtny1tDmifno/tM9sj5kXz8TO1J34YssXWHxiMWLyDBETkuqdCr2bHk0LmlZ5Nq21xgua1PIHOr+nP/vsM5Wy5cgyij0uUVFRKtLHH//mdW2aqOUxgSKzKvtT11ibV9oY2LVrl/of+vn5YXjP1ojFKhUH3g1DpNgSB9BCCbuLo0+rKB9HkFHccXxdY8LlUrG2oDJns4NpioFfJLysKf/KYKEufwWyxo5mjzSGrCxHz1SA6VKVYb2CIAjVLbz369kTkU8/zS86nPnhRwxetQruDigYr6zI33zKANfvGtEVt7e5HROOTYCP3sd4267QXVgZtRILYxbiQNABFLud68VnjcoaLzhCirfzx/ytt96qInbcF85uJRyvtW3btnP2V4sA1oafmyPRTIYbG1pTDWHmrYenYdpGAmKQg4qd3KYcQzSK4IEgt3wM6GLoiKXQd0X444lawnRpsBE7W7Degh/gZs2aVbielzVzw8pgyHf8+PGq5Z6ikPUcF110kc37aF1bGjNnzlTO1IIgCI6C30eWxjaFXnkFPMLDcPKBGYg5eQrDly7DyqH2zV5lWQuFD9OXzFRoESJ7hB2jYubwfvzTKEUpgouCccb7DHK8crA9fDv2hOxBy+yWaJPVBr5630qfs72NF9wfmjTzO5tpVqZmmaZlHTZHbWkTLVh7rVmsaBFBwbng9Aj+j5hC7du3LzrDUEu/G7bHbunhgSOIQxgyMaxbPOav2KS2pXXIupL/Id/fy5Ytq9Z9G5SwcxTPP/+8WuyFXWYchaMho3AEQagrYUcCzzsPzT/5GAemTkNEaipGLVqM4qQkeEZG2tzmxRdfrMpMmNKsig0I94MdqWxG0MxjLQlCdtf2SO+BzhmdVS0eI3YUePtD9uNA8AF0yuiEDpkdbD6W1nhBscb0qQYjdZYaL3j9ZZddpvxHKQyYfWFkjvWErJemuGPa1tRiRXCuaB3/T4T/q0jfYlVbV1kaVuMHXIgSSpsmQFjYUjVqjT8MXK0rWqsl1WDEjh67jU7Y8YuGtRfJyckVrudl1mDUFvyFyMedPn26WgRBEBwNv2Nov2QNvz59sGvKFLT//nsEZWUhf8eOSoUdxaI289I0imFvKtY8W2Htfp5lnmiV3UpF6kpblmJFwQo19YLdthrspDWN9tWk8YLXs4uWB3RGIvn9z6yNNtGCxMfHq6iQtWMDha6rT7pwRRiFPXHihDKkprDphP12pWE1Ss7KGkbpKOo0bcB0Pf/fTOvzfWGrG9sZ4PPnQlePDRs2VOmHV4MSdvx1RoXPsDsd0AlfDF6+6667au1x+aKbRuwEQRAcjT2iYuJd07GlfTvEFBcjaMyYWjPPtRY5rPQx4IZ2unZwT3RHunc6QgpDjLcdCjqkrFIYxTO1VtGoauOFaaTPdOi9NsEiISEBH3zwgbKzooDglALtNRD/u/qBx+vFixcby6IoqJPRBHvRGvvRskrbWfLfAnRp1xK7DyYoQU87nPruiq4ONP/mwogdbYJqTdj9/vvvVb7PmDFjVKt6TeEMONMPKc0L2eXKVAJ/UdHqhI7lHD3CF4N2J7RI0bpkBUEQXDViVxk8EA4zs4FiSrZw/34EOLCz35qws0ccagfSsEKD84AWrdsfvB/5HvlI8ktCZF6kVYFXFaxF+piKZZMcbVJoYsuFkTsKPL7OP/74o+oMZs2eVtPnKpEeV4aj4/hasw5y8ODB6rrDiFdLVehy5j8cua0UiwuaYsZfnioKyG0zu+ZKXdHVpVrCTouG2QtfSM7+c8Tcto0bN2LkyJHGyxRyhGJu7ty5uPLKK1XB7JNPPqk+vGyJpzO5eUOFIAiCK1GdyII+KwvHb7kFhYePIOqZpxEyebJd96tMoFkb1VTZ/aw9B0byhiYOxd6QvTgWcEyJOy5RuVHoktEFQcXVz4hYivTRfJ51eLTDoJ3Vli1b1PGCPnN8DjSqnzRpkhIYxBUjPa4GayE1X0UKL+21rw6peaXwD3NDV48zyoxbE3bUDvy/uUpXdHWp9juTHwKGO+1Z6EPjKEaMGKGKK80XijoNpl0ZZmfoct26daobqjZh1w7n2XHihSAIQn1F7MzR+frCp0tXhtiQ+NjjSH3/A4eY3lY3FcuDqrXHDywJRN/Uvjj/xPlont2cYTwk+ifi35h/cTCoPEvjSOh1Rl9UTrTgdAuKCe4fgwPM9lDEUfwRLdLD9CwjPYJjYQqW2TXWwzENS3pgF0JRng63lyMw1H7GuiWhc/tWKltIKzIGmIgrdUWzxo7a4sMPP6zdiB2jY1VJq3IIdUOuQZMaO0EQapvqRDDcPD0R9cLz8IiIQNqHH+L07NkoSU1Fs0cfgZuNiFNlkTeWuVQHe6JcmsCj8TE98E75nVLzaGsTHs+GDh2qDvS//fabaihJS0tTjRZc2rZtq56zZpJsK9IjTRdVh80SPI4S2o3xR0wQsnEJFqpu2FdwG/Jgf4DIN7oDTmQvRGxgGVq6Jyv7Mqbe+RhMsdNKxFW6ouusxo5O31Xhvffeq87DCIIgCGdhapAjEzkMvSpQpDW9/z54NAlH8gsvIuPLL1GSlorol182zputapMa01vVgfvCBgTTOmlrMP06MGUgsj2yldjT2B2yG956b9VhSzsVR8KDvWYDwxmlFALcV0Z6uGg/4K1ZWknTRdWhYPn555/Veb6vWrY0NEl0giG6dgwxVRJ1RKdzxyF9NGJxCnm7FyI62lCvx1pKejVyZKjmf9gQaVBdsfUFU7FidyIIQm3Cuq+adPeHXX893MPDcerhR5A9/x8kh4QgauZMq+uz0JzNBUxNsvRGg6my6qJEZtOmuPnmm/Hxxx/bdR9TUZfrkavq8MrcynAk8Ai6p3dHREEEHIU26YIRHdbUMbrDyB0FHpv0mM4jbK7gmCpGUrT6bYo6NldI0wWqFM1kDTzT3cHBwRXGhXY+a3OyqxJTYmtkhHTh7A70Ds3FnR/9aLyeU0lc6X9Rr3YndVHL5qxIKlYQBFcgeMIEeISGIvnFl9Dk9tttCjCmJrmY1i87iqpYq5jiW+KL7mndVYo20zsTy6OWIzYnFl3Tu6q5tTXF2qSLzp07K4HH6B0FCK0zNm/erJaYmBjVpEcxSFFHQag9P2m6sB3NZAqckVFi2qwSjCzEIdFuU2JLHIUhbd47Sodbb7wGO/cfUelY/lBxpU7YOkvFWuLyyy+XglJBEAQnx3/QILT89Re4mTRj6HNy4R5QbhZsSk2aLZg2Nh2LxM7HmmyTqdfW2a0RmxuL3aG7cTjwME4EnECiX6KqyWuX2Q7uZVVvMqnKpIsOHTqoYx0jKTQ9PnnypFoIRSAnEdG4WRN3jcFeozIsRTMZAZ0/f766ndFPGkZraGnYBMTaZUpsiTMIxjZ0wGm3cLSMj0Vk89YqFcuBBfzfVbWkwZWokrDjm9oS/JBqDs+CIAiCc2Mq6jL/+ktF8OJMaqGr04FrSbDRXmLTpk3Kf9SReJd6o2daT1Vnty18G1J9UrEveB/is+MdErmrbNIFRQgXPi/O2qVlCs9rtXhMVzOKR8HCbI499hoNtemCz4si2TSaSfHLwQGEz5XPm+tpz7emaViNX1Ce2qXOpsD+6aef1P+LnbcNdfxnlYQdPWa++OKLCrMBzWe7CYIgCK5BWWkp0ufNgz41FQk33IALbrsNq/NylfWH1ftUMdpmKe3qCMsVElIUgmGJw3DC/wT0bvoKoq7ErQQeZdVPStkz6YLHQhrp0hePRfnsnmVhPicd8HjJhc0ATNfastdoyE0XFG18XozU8b3ADtivvvoKRUVF6rVhup+vnRbN9EU+miINpXCrdhrWGrQlW7p0qUqrM2rHx26IeFTVQ045m1twMK9ul1RDQJonBEFwRWh50vyTT3Di7ruRt2Ytgt5+G9NefEHVkWlUxevL3to5Rwk79ZhwQ1xuxZm1KT4pWNd0naq9i8+JtzqD1lEwekchxudP0/xdu3apSB4jU5yOxIUwUkSvNkavPD09G0XThRal5PNiMw5NoGlwzdfs6quvNjYFaOvlw1fZm8QgCbmwXB5QFQKQg5Y4gX1ohSKdl9Iv3AfWRNIGxTxQ1RCaJ6oU52VLsiVRR/799180Vvii8w0rok4QBFfDPSAAcR98gMALxgHFxTg140GkmTRMjBs3rsL69nbF0jKkuhE788esKoeDDqPIvQibIjZheeRyZHvW7oQBremCEw4oGihirrnmGlV7bvp67du3T3XUzpo1S1nXcOIFu0K1NCWbLWgnozVd8HpG8qpyUHc2tB8GrKf74YcflKhjZHPKlClK3FoyCy6GJ46eNRmuKVPxAyZjPuJxQl1mepwRVkYMGb1zdtg4QW1x66232n2fGiXwTVvgBUEQBNeEfnYxr72G0OuuU5dTXnoZya+8ogQYIxoUGBqmlhS2BBtrzGrSAVsT+qX0Q5f0LnAvdUeqbyr+i/5P1eBxJm1toTVdUKiwUYIdjJqQ4fV33HGHqvFiZI/XsfGCc9dpocJoFZ0ltCkXVZ10QeHHFPCOHTvUqTMJQU6eYkaLzRJahouilQKW7xdTs2AdHL/fR892xzJqZ+qlSNjVfOrUKTQ0atQVyxeH89cEQRAE10/LciIFp1Scfv11dQDURJlpatbeEZG2BJ0jU7HWumfZJcvu2S3hW5Dsl4ydYTvVFIveqb1rNHu2Jk0X9AZk/SI7MynsKHZonUJvNS4LFixQY85at26tlqioqEqbLmpSn1ebDRu052DtPdPPmtCMjIxUEbPi4mIlgCnqGOXUzIIHYx26Yh9WoB92oIND9uMo4tAbO9ECx43XMQ1MCxumzDlphNGw6jYMNThhV9sfTkEQBKHuoBhrcust8O3RHX59+xqvpxhh6qo6tdTVScU66tjiX+KPwcmDkRCQoLpn033SkemVWWvCzp6mC74eFDhcuN7nn3+uPGCZAWNNHqN2GzduVIv22rHTlh2ctFExHedZk/q82mrYYA0ho4+sDaO4I9w/pl9XrVpl0ULG8Hhl6IHdCMcZuKN6s4htReyikAIfFKAAPsbIM2sf+XpRgLKDu6FQI2FXHyF2Z0SaJwRBaEj4m8yCLS0sRNbs2bjojjvgfnbkVk2py6AAGyda5LRA0/ymOBZwTEXxNPTQwx31F6nRJl1QzHEGO8UzU6n0Wzt8+LDq3iRMsXLRJpBQ4LHTdsmSJUowVdUUuaYNG+aRPtasUSSxYYR1hFqEjvWFY8aMMabye/XqZTVCGItEJeqK4IHdaOuw1zgbAUhFKJogA/E4iX1ora739/dX4o72J3zejOK1atUKzka9Tp5ozMjkCUEQGipJTz+DzJ9/Rs6qVWj+8cfwjIysdWuT2hB+tELpkFme3ivUFWJx9GK0yWqjltrunLV30gWjeEx3M0VLYdenTx9Vk8doHi+fPn1aLawPIxSBH3zwgYoAcrxZeHg4wsLCMHDgQBUNNDdFtuQrV5UpGRSFTBlz/6xB0UkbGJo5m74nbEUzGa1T20dbFKHqM4wri9o1QYZKx2rCjjAdy9ePKXEKPKZkTcsO0NgnTwiCIAgNj/CpNyJ31SoUHTyEo1dPQfNPPoa3lcgGD+iM3rBOzBaar1t9cjTwKPI887A9fDuS/JLQ53Qf+OrL05x1RWWTLkwjZ0xzUuBRrHG8GQUehRpr9riYogmqv//+W9XqsQmGUTKO8GL6lR5ujNKxiUGzXqGgpgUIO3YpdigQuT4fl1FFPi596KxBsci0vTY/1148UGI0Jd6KTnA0RxCHPthhbKAwfY3Gjx+v0uBc6K83derUCuluV6RGwq4hFRsKgiAI5+Ldti1afP0Vjt18C4qOHEHClGsQ9+EH8LVQb8eIzC233FIhSmMp+uYMjv8cP+ZR6oEdYTuQ4puC/2L+Q6/UXojJq3vRWVnThQbTh4yCcaGIYkSO4o+vN4UJhR6jepwExQYFokX4zPnjjz9s7lNltzOqSMHICBwbPZgWplBkqriqtMch+KIQmQh0mM2JKYfRHD/hAmO9nSkUtVdeeaUS1XydKO6uu+46db2rTgKpkbCjB48gCILQsPGMiUH811/h+K23oWDHDiTcOBWxc+YgYMjgStOxFCGsE3K2GjymXjl3NqIgAhsiNuCM9xmsbbYWLbJboHta9xpNraitSReW6vMYWWMKla+zBqN4FCgUWoygMeLGkWdcOBWD0T2KFb1er+r6mOolWic07+/j46OifDzVLEuYhj3vvPNUCpMpS9P/NSOA1Z2Hq6Vht6EjymohJZ4PX5tdtnwdr732WsydO1fN/WVqm68NrWhccRKIpGIFQRCESvEIDUX83M9w4u57kLt6NU49/DDaLFwAXSX2Jyyer6qwq0vYITvy1EjsCt2F/cH7VYrWs9QT3dKde5qSpfo8rQmCNiJsvrDUBEHR9tZbb6l1tRo7CmntlNviNu6+++4KESo2b1DYsVmQ6Vtz7JmHa41V6IMieNZKGtZeuP8UdxxvxrQzI3aTJk1SgtnVJoFUS9hxZIq9vP7662joSFesIAiNAZ2/P+Lefw+JM59CyOWXVyrqiFa/5cxWWvS965rRFZH5kdgduhsdM5z7wF2d+jx7BaGpr5wp2mQIrsdGC3MsTZCwF6ZfayMFa4ovCtALOxCKTPyJ82AJNqBQtDJax1Q209E8ZTevPY0lLt0Va56CZXcOQ7l8soRvDAqd3r17ozEgXbGCIDQW3Ly8EP3iCxWuKzp+HJ6xsXVugTV8+HAsW7bMYdtjWnZ44nDjZU6q4MQKpmd9Sg3+Z86GvfV5NRWEWuqXkSvTblpiPkHCGXFDGcZgpTq/BIOQi3N/lGivIcfBUVDRP/DPP/9URsZssqAIrm662em7YlkkaRqR4xuJRZxaJxTDmOwsYc5dEARBaLgU7NmDhOuuR9D48Yic+STc6qiprlOnTspU1pHCzpzDgYexK2wXDgUdQt/TfdG0wJBudDaqWp9XHUFY3UhfZU0TtCDZiO5Ig+1O6pqSB18koQkikaoecxcMgShTtDQyXwc2hKxZs0bpHXr0vffee+jevXuF9Rpsjd1rr72mVL9pezvPP/fcc+pN8MADD9T0IQRBEAQnpWDvPpTm5uLM999Dn5GB6Fdfga6GXa+WUrE04mUEpS5hBC+oKAhZXllYEbkCHc90VEt9eN45gyCsTqTPFgOxGS1wAoXwwlIMQm1zFHFK2LW0IuzM082DBg1Sz2n+/PnqvadlK2nETJ1D2x7zKHVtjmmrM2HHrhFLrdS8ztlVrSAIglAzQiZdomrtTs2Ygex//1Wds7HvvA33gIAKtUu046gJ7MSsa2GnNVZwHBmbKvaE7sFpn9Pod7pfvXjeuWrq1xLNcFqJulK4YTO6oi44hHgMwBa0RoJKtDNBW1m6mQJuypQpKmr3ww8/KF8/mhpzoV8fI8dcOGXDWeb21ljYsWuEaVdG7pgHJpwT9+CDD+LSSy+t6eYFQRAEJydo7PlwD/4IJ6ZPR966dUi4/no0//BDeDRpUu1t0iSWB1FTYUf/O5rn1mWDBW1Peqf2RkR+BLY02YJU31TledcvpR+aFVTNiLehUJ3UrzkDsUmdcnxYFqrecFEdjiIWJXBHKLIQjgykIczudPPatWvV+5H2MbSMYd2dZgzNdC3FGEUZvf04qoyvT33N7a2xsHv//fcxY8YMpWg1Q0QPDw/cdNNNeOWVV2q6eUEQBMEF8B/QH83nfY7jt9yKwt17cPSaaxA/7wt4NqteXRqjfIySmB50q3KQ4+xV1n47iua5zRFaGIp1Tdep1Cw7aYXqEYRsdMU+dX416q7JshieOIZotMJxtEHCOcKuKunmsWPHYu/evUqU0VpGy1AmJibim2++UUKPIo9RPd7GdC5n0Zqbc9d0bm+tCDsaF7777rtKxDE0SVh0SIdsQRAEofHg27mzYUrFTTfDKzYOHqEhdt+XBzAeyOylsg7cli1bwtEElgRiZOJIpHqnqvo70+7ZhlJ3Vxf0xxa4oxRHEItTsD57uDY4iBaIQ6KyP6lJupnahzYoXNg4QjHHqLJmAM37ceFtGuxq5bY4qo1CkTqJdXsUgJyty8uMRLN2z9xepV4MirlD3SyMmGkMiI+dIAiCAa8WLRD/zddw9/dX1ij2EBcXV2FygiXhZn65rr3uNNzL3CukYLM8s1QUj7NmQ4tqt7OzIeCNQvTGjjqP1mlsRDesRw+UVCJ/qpJuphUJueiii5QPHi9T3DF6x9pS7TzRBF9CAuv8DDDF+9lnn1V4bE78YPaTfQwvvvhilertZPKEA/jwhR8R4G+oEdj0z1Hj9eExAWjR1VBjoi8pxdb/jlndRmikP1r1iDB+YW1eUP5PNyc4wg9tepenN7YsPGbVvDAwzAft+pX/Itq26DhKivUW1/UP9kaHgVHGyzuWnkBRgWHUjDm+AV7oNCTaeHnXipMoyDWk4s3x9vVAl+HlhpZ7ViciL8vwQTDHw8sd3UeVG1XuW5eEnAzLv6x07jr0HFPumXRgYzKyUstrcszpPa78Q3p4y2lkJOdaXbfHmOZwdzd8kI5uT0XaqRyr63YbFQdPL4PFw7FdaTh93HrTEF8Hvh7kxN50JB8tH1ljDl9fvs7k1IEzSDxUXn9hDv9v/P+RpMOZOLk/w+q6fD/wfUFSErJwfE+61XX5PuP7jaSeyEHCzlSr6/L9y/cxSU/MxZFt5zZVafBzwc8HOZOSh0ObDeamlmjeKRwRzQ2fr6y0fBzYUHHYuSkx7UMR2TJYnc89U4i9aw1fppaIahOC6DaGiFJ+dhF2rzpldd1mLYMR295w0C7MK8bO5Setrst95T6T4kI9ti85bnXdhv0dUYJOQwwRiFYrVgB+XZAZ3NT4HembU/4516cEVRBu6Yf0KDwVWGEdPgd1v9xo5PuX/6+885tCpz+3C5eP45sbg3z/8v+Vd34T6Gw0PeQHlP+vvArC4V5i3YA53/8E4FaGnaE7EZIbh4wyd/gW+6t0rXlRvtoHN8Pr71kYCo/i8sYScwr8TqFMZ3j9PQtD4FFsvf6swC8RZTrD6+9ZFASPIsN73xKFvskodS9S5z2KAuFZZD2aWuibglJ3w3c095X7bHVdn9Mo9TB8R7sX+8Or8Nz0pkaRTyrgUajEVZOSYpwsGAZfK5HOIu806D3z1HldiQ+8TSKk5hR7Z6DE0/AdrdN7q/eE1XW9zqDEy/AdrdN7wTvfep1ksVcmSrwM39FupZ7wybMcXWwbEoC20YeN3bQ+nn5oHzYY7flSdAYyM89g48aN6Ni5I0rKCpGSfQwZeYk4evQoDh8+grjIFsjJy0JWzhmU6EvUZzUvz/DciTbyzV5E2DmADX8dha/XuannDoOiKnxpr/31sNVt8EtY+9Jms46tdeO7hlf40l73x2Hoiy1/ace0D6nwpb3x76NWBVjTFkEVvrR54MjJsCzAwqL9Kwg7HgwyksrfiKYENfGpIOx4MDh9zLL48Q30rCDsKBgTD2ZaXNfD272CsNu7JhHHdlkRKW4Vhd3+9Uk4tMW68FD7cNaO6+CmFCUwrcHXQRN2h7elYpeNgz7/F5qwS9iZhq3/WT/ot+weYRR2FF/831kjrmOYUdhRBNp6/0S1DjYKu6TDWTbXbRIbaBR2FIG21uV6mrBLO5ljc13/EO9yYZeUZ3Ndvl5GYXc63+a6gy5tYxR22RkFNtftN7GlUdjlZRXZXLfn+c2Nwq4gt8Tmul2Hx1QQdrbWbQzfEdn//IO4TZtxsu8EBLhHGvc7AOWp0pKz9dkaaftLUZQejACUCxXtfn665krYaULQJy8KXsXnChqu7+/WooKw88mPhFeRZeHBdKqpsPPJbwrvQutiwrDdMuVvp89vhrgzXWwLsLPCzrsgHL755d+dlkRVGQzCzqswFH55cTaFkl4TdoWh8M+Ntyl+NGHnWRSMgBzrqWqKGaOwKwqyua7eMwdFZ4WdZ3GAzXUzPfJR6JGP/zAU3sVNEJTTyuq6We6FRmHnUeJnc7vZuhKjsHMv8bW5bk7gYaOw89K721w3N+CoUdjp9J5W1w3wbImebRKwYsVClULVlVZcN8AdiOnf03i5ZbOTyAk6pERgyql0PHTRR8bbikoKkVeYhbzCHOQVZGPz4WXQh5xWtXmLFy+GPYiwcwDt+kfC3zfA4gFUQ+fupr7ErdEsvuKvMlvrNjl7QNRoPyASpXrLaYnQZhV/cbbt10wdbCwRFF7RWZ0HhoI8y78UAkIq/kLmASc3y/ClYY6vf8WRQi27N0F4rOVfrF4+Fc1N47uEI7ip5V/N7h66c8SN31lxY46bhciO51mBZXF9Xfk9otuGwM3dev2M6X5EtgpSB2hraALQeJC08X/WBCCJiAu0ua6PyWvMA6qtdX0Dy9NjoVF+NtelANMIaWp7XQp44/lwX5vrBjcpj5oEhHnbXDfkrFgk/P/aWjcspnxdimKbnyOT96C3n4fNdZvGB1V4j9r8LLcq/9y7e+psrtsYviMCx4xBYqdOiEjdiqKsYPh06Qzvdu2wZetW47qBwb4VInaBMToUexhSVho9exoOjJu3bqiQii3yTofe49xIPdffunVzheuKvDOgPytuKqPI6wxKz0bOLGN4fM8yTwSUeiI5aB8yvA2Rcq9SL9VF61GmfS7LX/tiJRSsp9XK3Mofk2Il3zfJIeuWupX/n/SeubbX1ZW/RnqPPJvr6iusm297u+7l2Re9e6Ht7VZp3fKgQmkl65Z45CISKZiM+cjSNcUvvtaFc4lHeaaGkVFb220SF4D9y/erurgRQ0bDxz8G+QUFKg2blZmpaj6DQww/JIs9y42Q/QJ8sStxNVqxJvTsZ8A7EPAu80HGkUR4eOkwfNQo1Zxqr7BzK6tmocKTTz6Jiy++uNGMDbMEc9/BwcHIzMyUkWKCIAhWeP+999B04UJ03LNXXQ678Ua8l59nPJDxADdt2jQ89dRT6jItJQxpqvKopHabdsoC9yuvvNJ42Rxez4MqJwbUFSk+KVjfdD0K3QvhqfdU48mCLUQTG2tt3SQswBr0QgLOnTVbl/ghHw/iffWD/zXcgmxYT41XherYlph2xVqb5sH7aiPF7NEb1Y7YnThxQnm1sFBw4sSJqmhw9OjR6rIgCIIgaLjpdNjevTsKvb3RY+s2pM+di6tGjMB3TSNQptNVqxHCnrm0PEjWJRw5NurkKNVMoSvTIdBGfVxjoz+2ogMOKf+4d3G96iOuL/Lgi1Nohhgkow2OYgusp9GdcW5vZVRb2HEHWOC3atUq/PHHH7j33ntV18eYMWNUJO/CCy9ULb2CIAhC40YTbvs6dMAFV12FxMceR9ny5WgyYjhO2yG+rrvuOqvbrKn4czR+ej8MSxwGvU5v9Lor5Z9bqTI7boz4Iw+DsVGdX4b+9SrqNPajlRJ27XDYYcKurub2VkaN3mV80KFDh6pl1qxZKqRIkffBBx/g1ltvVZMoGMm7+uqrlS+LIAiC0LgJueQSuIeEQJ+aCo+0NDq6GoerW4PeqK6EO/9Ky+tpd4TtUKPIBqYMhH9J4/N4HYE18EYRTqGpxRmt9cE+tMJIrFHjxTxQUqn9iStM89Bw6DOh6uTy0EMPqVmxv//+u1oIp1MIgiAIjQ/zyFngiBHqdGpRkfL5inR3R7HJLFlG4yqLttVHNK46FOoKcTzguKq7Wxy9GANSBlQwN27oNEGa0bduAYY7RbSOJCECmarvOgctcQwHYL1D19WotZkoERERaqzYb7/9JqJOEARBOAfWZEcHBOD4TTfj6NVTEJhl3dfRnPoyKK4q3qXequ4upDAERe5FWBG5AocCDyl7lcbA+VgBHcqwV8XG6rdpoiJu2AdDJLg9rFsHuSIy7E4QBEGoN8oKDLYWJYmJGP3fIoQxPdvAYN0dO2TjcuJQ5laGrU22YnOTzdCf9atrqMThJNrhCFht+C+GwtnYhXbYgk7YgzZoSIiwEwRBEOoNz+hoxH/9FXy6doV3URFGLFkKr927K72fq6RiNdg4QTPjrmldla3d0cCjWB25ukFH7o4jGr/gfCzDAKTB+catMYL4G8biEBxT2+YsNM4WHQcjs2IFQRCqj0dYGOLnfoZl4ycgMjkZwW+/g9CJFwI+FQ2RXTEVa4ob3NAuqx2CioOwIWIDWmW1Utc1XNywjTO1hGqzfv16bNiwwepIQEuIsHMAfNHFoFgQBKH66Pz9sWLYUPRftw7Njx1Hq19/Q8rAATgWb31MlqsSmR+JccfHqakVGsVuxRUuuzLByEKBsiS2PAnImXBDGaKQjFgkYj3Kx345C3QX4aIZFNdKKjY/Px8nT547C3PXrl1V3ZQgCIIgGCl1d8eagQORN2IEikJCkGLD487VUrHmmIq4XI9cLIxdiL3Be10+NatDKS7HX5iOz1WNnbPjiwLcjG8xHksRAstzyV2NKgm7H3/8EW3btsWECRPQrVs3rFu3zqaBpCAIgiBUCTc3ZF95BfbcNA0Fvr5WU6+umIq1xkm/kyjwKMCusF0qRas3mQHragzGBsQiCV4oQSacP5OVB19jt24nHECjE3bPPfccNm3ahK1bt+Kzzz5TdiZff/11g/uQCYIgCI7Dx0atnEXc3KA3EXVnfv0VJ++9D6VF5QPnnYVevXrVeBusu+uZ2hNuZW7K82555HIU6Azdwq5EJFIwAmvV+b8xAllwjZFqu9C2QQm7KtXYFRcXo1mzZup87969sXz5ckyaNAkHDx50+bC4IAiCUDtwAhEzPoMGDaryfb0LCpD0zLMoy8vD8TNn4BEbgxJPT6c55kRFRTlkO62yWyGgOABrm65Fuk86FscsxqDkQQgpCoErwOkNl+IfuKMUu9EG21H1Gaf1xV60wQQsVpFG1ge6QqTRYRE7DlTevn278TJnwf77779qlJjp9YIgCIKgERoailtuuQWdO1e9Q7LQxwdx774DnZ8f8tatw8jFS5TYq48skbe3NwYPHlzhOkeOy2xa0BQjT41EQFEA8j3ysTRqKdK8XcPXbyyWoSnSkAM//InRqi3BVciBf4NKx1ZJ2H3xxRdK3Jk7h3/zzTdYtmyZo/dNEARBEOA/YACaz5sH97AwhGVkYNSixfA8c6bO94Pi1NwBwdGRw8CSQIxMHImI/AgVwQsuCoaz0wEH0RfbVdvHrxiLPPjB1djVgNKxVRJ2sbGxiIyMtJiibd68Ofbt24f09HRH7p8gCIIgwLdLZ8R/9SVy/fwQlJ2NNp98isID9X8Q9vNzvIjxKvXCkKQhGJo0VBkbE3bLlsJ+L7O65AjilDBajv446KJmv3vQVgnTMJxRRi2NcvJEdnY23nvvPQwfPlz9gmnRogU6duyoZsTGx8erXzb0dxMEQRCEqsA0q6VImHfLllh03mhkBgfBMysL2f/9h/omODgYF198scO3q4NOzZnV2Be8T82ZLdI5XwMJ/ep+wAQsxcBaf6yQkJBaS8d+hKvxGm51Cf89hwu7119/XQk5dsaed955+PXXX1Wn7P79+7FmzRrMnDkTJSUlOP/88zFu3DgccIJfVYIgCILrk+/nh8WjRuHU+ecj/Pbba7w9rSHQHqylXXv2rF1jW3bI7gvZh1TfVCyJWoIcjxzUN24oRVfsVbFE7ZqyOqir46Sn2uIUIlHaACatVmvyBCNx7Ii1VghLl+Rp06bh/fffV+JvxYoVyv9OEARBEGpKkbc3Unv0MAot95ISRJw+jSQHdahWB19fX2XgXxv4lPpgxKkRWBW5CjleOVgSvUR1zIYXhqO+OB/LMRBb0ALH8QfGoCHhhjLV5VsMz8Yj7NgsYW8H0e0O+EUlCIIgCJZwKy3FwNVrEH3qFDb36oWD7aoWRHAW25TKCC4OVh2zq5utxhnvM8rrrm9qX8TmGro565I+2KZEHTmEhjXyrQd2YSRWYzO6YFkdpJZrA9ePOdYCr776qopGdunSBV9++WV9744gCEKjQ/O869Chg831ytzckOfvp5KAvTdvRucdO1mkd856TZo0sXz/Ktim1LcI9NX7YnjicETlRqFUV4p1Tdep2ru6pBP2YzyWqPOLMAi70a5OH79r1661uv1S6BCMHHSrkGZuBBE7Uxh65gdD6wxKSEjAL7/8gk6dOqkaO1djx44dapoGJ2zweY0cORIXXnhhrRVsCoIgCOfSsmVLzJgxA/7+/rZXdHNTkboCbx903bkTXXbtgk9hAcr0eri5u+ORRx5Rzg3fffedQ/bLmrirK189dskOTBmI7WHbcTD4ILz0XqgrWiMBkzEfOpRhE7pgBfqhrjG3m3E0e9EaRfBAOM4ow+ITqL/0fr1F7NgNNG/ePHX+zJkz6N+/P1577TV1PbtmXQ2aLQ8cOFCNwGHNRPfu3fHPP//U924JgiA0OgICAuyLkrm5YXeXztjYu7eKsbQ5eAgnZ8xQI8hYEsTtuHoq1hQ3uKF7endVd9cyp2WdPGYsTuFK/K4mS9DaxNVMiO2lCF5qEgXphj1wRWos7DZv3oyhQ4eq8xwZww4jRu0o9ubMmQNHw6aNiRMnIjo6Wn0g2ZFrzjvvvKO6dinOKDTXr19v9/aZfl26dKkSqRkZGer8yZMnHfwsBEEQBEcyevRo6MechzWDBkKv0yF7/j9IfPgR4+0NcZ65afNEoa4Qq5uuRq5Hbq08lh/ylag7iHj8jAtQ1oArubadHYfWBfuggx6NLhWbl5eHwEDDoN+FCxfi0ksvhU6nw4ABA5TAczS5ubkqisauWz6WOQy333///aojl6Ju9uzZGDt2rDJP1qZm9OjRQ9mxmMP9Zwr5nnvuwahRo5Q/EZ+Hu7u7w5+HIAiCYJnqiDAt8na8eXMUeXlh9K5dCLvxBoful63oXn0Lxy1NtiDRPxHp3umqYzasKMyh29+P1vgSk1RqUo+GfUw8guZqNFoA8tAGCdiPVnAlaiy527Rpo6Jmx48fx4IFC4x1dSkpKbWSC7/gggvw3HPPYdKkSVY99miOPHXqVCXSKPBY//fpp58a16Hn3s6dO89ZGAUkt912m4pELlmyBJ6enpVatWRlZVVYCgtd27VaEATBlUmOjESbf/+Fb7du5VeWOmZqA6csOSPd0rohuDAYhR6FWB61HCf9ap5paoljCMWZCoLHVS1AqtpAsQPt1fnu2I36gMEnagnTpc6E3ZNPPqkKXJn6pH8d69O06FdtmzaaU1RUpJoeaJqsweghL9M42V4oSgmjfEzjMuJni7i4OBXd05YXX3yxBs9CEARBqCk6Hx/j+fwdO9Bz7ucIzMyq8XY5VpMlO6R169bG64cNG4b6xE/vpzpmm+U1g16nx9qma3EgqPrDAWg+fA1+xQ34CYGof0PkumYrOmMDumElas8Q2RYrV67ESy+9ZFzeeOON2k/FUtCxQeKyyy7DkCFDkJiYqFKkpvUO1qJqtUVqair0ev05TuK8vHcvW5ftg88rMzNTdWPRYNnDw/bLxGilaXSSxbqCIAhC/cMUafLLL8MvPR2jFy3C8uHDkB4eXqNU7OTJk1UQw/RYE1WP5sganmWeKg27LXwbDgcdxvbw7cjzyEO39G6q4cI+yjAEG3AeVqlLJ9EMefBFYyMZEfhLNYjUD9RVWqCMMGJnr7irtrA7ceKESot6eXmpZoaLLrpIeb/xMmH0zlWpSnSPUNTVdgu2IAiCUD0xFvvWW9g6+TL4nzqFEUuWIuzll/BxDWaZc5sxMTFOVWNnOmO2R1oP+JX4YWfYTiT6JaLjmY7wKq3cFoWNAvSo64Md6vJq9Ma/GFono8KEijCgVFlQyeHCjjVrpaWlWLVqFf744w/ce++9Kmo3ZswYFfGi91tYmGOLNyuDBpRsdEhOTq5wPS8zfF6bs+v4uNOnT1eLIAiCUH0cLZI8QkNx6Prr0PSzuYhMTkb2w48grm8f1WhhDwxYUMix1MYVYHSufWZ7+Jf4q7o7e0SdP3JxOf5CC5xUljHzMQLrUbflVM5IDJLQGzvUJIr68LRjORjHuFJv1UmNHevXaHUya9YsVY+2bt061Yn6wQcfqEYE1hxwikNd2YXww9e7d28sWrTIeB1fDF42DWk6Gr7ou3fvFlEnCIJQj/CYZI1Sb2+sGDYUx+LigOJiNYas9YGDdm13ypQpuOGGG2x2xdp67JowePBg5dJQHThuLLDE4FpBTvidsGqHMgYrlKgrhBe+xUUi6kzGp/XCTiXu6gNmP6ktbr31Vrvv49B3YseOHfHQQw+pKB7rzvhBWLFihd2zZe0hJydHdbVyIUeOHFHnjx07pi7T6uSjjz7C559/rsyG77jjDmWRwi5ZQRAEoWFgSWRNmDBBlcWMHz/e4n1K3d2xduAAhFx1pUoucr6spfFjpuOrGKSwpxM2Pj6+1jJRjmjMSPZJxvqm67EkagkyvDLOuX0BhuMAWuAjXI19KG8KaexshqFRpjP2wRuFjcPHzhoRERG46aab1OJINm7cqMZ8aVDIEYrIuXPn4sorr8Tp06dVc0dSUpLyrOPkCPOGCkciqVhBEIT6T9eGh4cbjwlW76fTIXLmTPx76BAS4uMRYSMKxyYJPo490ylqK2LnKIKKgxBUFIRM70wsi1qGESk9cUF+Ktagl0re5sMXX6FuGx5dgeOIxmmEIQLpyrB4E0wsdJw0FesQYVdQUIDt27crmxDzB2dThSMZMWJEpfUXd911l1rqCr7o0jwhCILgGlCoHTaxKmHUrsWRozgW31xF9czXtRdGDOfPn1+lg3Bd4av3VXYoa5utRYpvCpY224ChqenonuOLbehU37vnxLipqN1YLFdNJZvQtU5HqTEVy4VdsbQ9qRNhx2jY9ddfr6xGzOEHgvYjgiAIguBIHDnjtfPOXeiyaxeaHz+GVYMHQ1/NbkRmb1jn/cwzz1T5vgwO0OC+Ngks0+P5pAL80SQHvwcG4KmIcPTxyEXzM+x7lc5Xa2xFJ4zCKkQhBbFIxAkYhhk4KzWOHd999924/PLLVUcsf6WYLiLqBEEQBGfvnE0LD0eJuzuiEpOUHYpXDaYHVTclyxr12qQjDmA65mEAduK51HRcmGEwcN4YekxZogjWYZp6Jzqo8/2wDc5OjSN2tBJhTUNt1rA5O1JjJwiC4LpMnPkkkpYuQ+lLL6FJWhpGLV6MZcOHo6EwFksxEFvU+TSE4HeMgfeZWPTSH1GNFFF59W+u7OysR3e0QgKSENHwa+w4eWLp0qUVRqs0NqTGThAEoW6jcY5MxdKjLuzSSXh/w3oMW7YcwZlZGP3fIhQeOQLvli3h6rDLlZGmVeiD5eiPkrOH/pbZLcE/jRK3EhTrilU9nlCRRDTDbNyEMseaiThnjd3bb7+tUrG0NWFruKdnxQHB99xzT00fQhAEQRBqfdJDZkgIFo0ejRHLliIwOwfHrr8BrRf8A52fH+rr+VRVwLqhDF2wF94owkYYxnweRZwSJdkIsHq/UpRiXdN1yPTKxJCkIaqLVqhIXYu66lJjYUePuoULF8LHx0dF7kzfhDwvwk4QBEGobRwVwcsL8FfijpG7zvfdVyeizjGUoR2OYCTWqCL/IniqSJ0m5myJOlLkXoQcjxzke+RjadRSDEwZiIiCuk07ugI66NEBh9T8XArmBinsHnvsMTz99NN4+OGHnd7Hp7aQGjtBEIT6FXKOjOBFtGkDz4svRojJxKLSoiLozs5Crw2qv/9laIVjGIXViEWSuqYQnliBfsiHoUHCHnz0PhiROAJrmq1Bmk8aVkauRN+UvojNi63mfjVMBmETzsMqJCAan+HKhlljV1RUpEyBG6uoI1JjJwiC4PqpWA0e0wIDy0dxlaSmIuG66xE2bSpCL78czkJTpGICFiEep9TlYnhgHXqoWjp2clYV71JvDE0aivUR63HK/5RKzRakF6BNVpta2HvXZBs6qagoX3POkT2J2ptDX90auxqrMU58+O6772q6GUEQBEGotmhzZDOFOWd+/AlFR44g6YknkfrhR7UqKqtCAbxVlK4E7liLnngT0/AfhlZL1Gm4l7ljQMoAtMpqpXx4t4Vvw97gvQ7db1cmGwHYcdb6ZBA2whmpccSOXnWzZs3CggUL0K1bt3OaJ15//fWaPoQgCIIg1Bvht92K0txcpH30EU6//jr06elo+tCDcHNgpqoysRiAXDWMPgRZyq6EZCEQP2McjiGm0hq6qkCz4h5pPeBb4ot9IfvQNL+pw7bdEFiN3uiB3eiIgwjFGWQgBA1K2O3YsQM9e/ZU53fu3OmIfRIEQRAEp4HRwKYP3A/30FCkzJqF9Llzoc/IQNRzz8LNLJjhWMoQmLoZ+P513Ic/4A5DndUa9MRpNFHnd6F9rTwyxV2HzA6Iz4mvYH9SBplSkYImOIAWaIujGIjN+Buj0KCE3ZIlS9DYkeYJQRCEhk/4tKlK3CU+/jgyf/sN+sxMxMx+Azof+5sU7CEAOeiOPeiJXWiyMkNdxwm2xxClLEzSEYq6wlTUpXulY2fYTvRP6a/q8RozjNpR2PF/tAQDa5T+drrmCdqdXH311RZve/DBB/HKK6+goSPNE4IgCI2DkEmXwD04GCfvu0/V3ZXm5TlE2JmmYlviOMZgpTqv9/CFe48peG9jIZLreOpBhf1DGTZGbES2VzaWRC9RXncBJY5L/7oaRxCHRDRFMdzhj7xaE3b10jxxxx13YP78+edcf9999+HLL7+s6eYFQRAEwakIHDUSzT/9BHGffAKPsLBqb8cLReiE/bgMf6FT1jLj9XvRBgcRj99xHnaP+wW48PV6FXWE6Vd62/kV+yHXMxdLo5eqcWSNFzd8jsn4FFchFeFwJmos7L766isVsVu50vDrgtx99934/vvvJU0rCIIgNEj8eveGV2yM8XLmX3+h6NixSu/niwJ0xy5chd/wEN7HFfgLXbAfLdIp7AxRu2J44ktcis3oilJP5zFIDiwOVF53wYXBKHQvxPKo5Uj2Sa6zx6e2MOXCCy9ESEj9NS4UVMEn0KWE3YQJE/Duu+/ioosuwqZNm3DnnXfi559/VqKuQwdDS7AgCIIgNFRyli3DqQcfwtEp16Bgr3VrkMvxJx7E+5iEheiAw/CAHukIxkr0QfGEN+EKsOZueOJwRORHoERXglWRq3DMv3JB6wjCwytGxgICAmrV5sZefFCAwdgAd+jhDNS4xo5MmTIFZ86cweDBgxEREYFly5ahTRsxNBQEQRCqjrP4xNmLT6dO8G7XDoV79yoj45i+3RDSNBfNcQrzMUKl7UgpdNChDElogj1ogz1oixSVxnND7yi6Syx02D5R8NTW6+hZ5onBSYNVzd2JgBM44X8CcblxjbRbtgw341s0QYaa9MEoq0sKu/vvv9/i9RR1vXr1UhE8DfGxEwRBEBoyHr5liH/0chx/6h3kH07H0GUrETMoHYExhdiErsoegyzDACzGoDrxPYuPj8fRo0drbfvucEe/0/0QXhiOFtktGqmoI27qfzwWyzEEG7AVnZWAdzlht2XLFovXM0qXlZVlvN0ZQqR1gdidCIIgNBKU7UQZoKMBCYBVbwL/PqnsSJr3Ak4WhCHnlA9OrAxDTt8QFLUq97lLRdUbLZz5OEoxZzpujJ2zRwOOonlOcyX8qppmTUtLU+nVnJwcl4rubkQ3JerCkImu2KvGjrmc3Yk0RVRE7E4EQRBcG9N55xVmn+ekACc3ASc2Gk5PbQYmfQC0v8Bwe9OzB/FmXaCLH4zYK/rj34feR9zREwhYnwn3oGKcDdg1ePaE7MGe0D04HnAcA5MHqpStvVx77bVYsWKFKul666234EoUw1P52tGiZijWYzs6oMxBUbvq2J1US9gdO3YMzZs3t3v9kydPIiamvHtIEARBEJwBLSLm7e2tRAUjQP45CcBvTwNJO4GsE+feiQJPE3YthwEPHQH8DNE4bm11/+3o7r0N7qV6pJkV/DdkmJb1KPXAad/TWBa1DIOTB1cwOLZFaGioasJ0VTagOwZjo6q164wD2FlLE0HsQVfd1ONtt92mIlXWyMzMxEcffYQuXbrgp59+qsk+CoIgCIJDcGMa9cxxtEIC+mMLhp35AfhwJLD2fYwZMwbnn38+4O4J7P/nrKhzAyI6Aj2vBS58A7htOTD8/8o36OFtFHXlD+KGbT17YHOvXuo8cS8pgVsV0mmuSLP8ZhiWOAzeem9kemdiadRSZHtkozFQBC+shWG86jCshdvZ8W/1QbUidrt378bzzz+vPgQ+Pj7o3bs3oqOj1fmMjAx1+65du1QjxaxZszB+/HjH77kgCIIgWKMoFyjOB/wNeVDf4gzchq8QjnRg9mxcr62Xf3YJb11+37BWwAWvAJFdgMiugHdg9fbhrKjT6fUYvHIV9O7uWDNoIErd3V2qlq4qhBaFYsSpEVgZudJoZMzIXVih/fWFEydOxB9//FHhuvr0q7OX9eiB/tiK44iGJ0qU2HMZYcciR3a7Utz99ddfypw4ISEB+fn5aNKkCa655hqMHTtWResEQRAEoSrYXRCvL0ZLHEMoMhF/8BSQ+AmQkQCcSQByTxuibBe/o1YtdA9AM5xWdiPQeeB0aRDSEIrcgBbofcF1QLQh2qJgY0T/Wx32fELOnEHTlBS4l5Zi2LLlWDl0CEo87a8/czU4aozibnXkamR4Z2BF5ApccPwCeJXaJ3RatWp1znWeLvB6FcAHb2JavQk6h/jY+fr64rLLLlOLIAiCIDgCtzI9kHkCyEoEsk8BWeXLECRiJfoZVizV4wacLfU5YmFD2UnGs6U6T3yFSTiDINz92Kt459nn1fWRAZHo3XlSrT6f9PBwLBs+HENXrECzlBSMWrxYXS50wIxZZ8Wn1EelZdc2Xas87uwVda5OUT2LOocZFAuCIAiCTRiFKzgD5KYaomnsNuUpLwfHAL0MyVEd9Dhv3bXAWssu/q0QVy7sPH1wFDGqKzEgtjOiOg4AQuOBkHjDqW9ohfseQrzhDGvo6pjTzZpi8ehRGL50GUIzzmD0f4uwbMRw5AYEVHub/v7+yM3NhbPiUeah0rCmHndFuqJqiTzTVDWbNxMTE9GyZUu7o7tRUVHqPnVFBFIxEFvwN0aipI6lVrUebfv27SrNWqEl3Aast2vfvj08PBqmjhQfO0EQGh08oKYfNoi1/Awg3+SU14W2APrdYliXTQMvRAElBZa31WKoUdiVwh3FHgHw1ucCgVFAUHT5aVA01izcWOGuc3GFOh3TcQyiBg+u8tOoy9q2M6GhWHTeaCXuAnNyjOIu00b9mLX9Y0kULTCcyd+N6dLi4uIK15mKugL3AtVQEZMbgy4ZXayaGle271OnTlW+bjzu2suoUaPUbPu6gI0TU/AbQpGF0wjDGvR2fh+7nj17IikpSU2asIeBAwdi69atFvPmDQHxsRMEwWXQlwBF2UBhDlCYbejq1BoHePBY977herVONlCQVS7eYnoDE14r39bbfQGmTa2JNU3YMQjABgQKO+9gQ0ODfwQQEGE41bzgzrKq++sYNf5Sw/3MOLDwqWo9bWdpTsgJDDSKO//cXOiq2SnLOvbff/8dzkRlgizZN1k1VOwP2Y9C90L0Su2lhqxV539ZFVFX19DDbjn642L8q3ztNqMLCuHt3D52/Oc98cQT8PPzs2v9oqKi6jyMIAhC4xReFEBcivMMnZ089Q4qF2AlhcDWr8/edrb7U1uPp2wEGHBH+brvDjwr1nIM65jS4ULgqrORDIqff58ASkss75uXSdqQ6zKFSmHClKdvyNmF50OB8LYV73v7KsAnWKVPK6PII8CiqKuJcHOGCQUaBb6+Ki0bmJ2NjLCqT6MgzJg5i1i1l/iceDWdYnOTzUgITFDirn9Kf3Rp3/AaLbehk9HXbgA2YxkG1tljV0vYDRs2DPv27bN7fUbs2GghCIJQ5/CAri8GSovPnuoN5z39AJ+zkXaKodQDZ9cpMZxS3GjnQ1sCTTsY1mWqcfev5bfrC4GSorOnBUDzgUCHCYZ1WT/2x/8M4qrCemdPWbQ/6vHydV9tC5RZieB0nwJMes9wns/hz3utP2eKN03YuXsBGUfPjay5ewPeAYbXQYNCoccUwM3dcBvFJCNtFGo+IYZ0qCn37rD//xDYzP51GwHFXl6qqUKjcPt2ND+agGMtztYB1hBnFX0tcloon7t1TdchyS9Jdcy2d2/vEvteFTgvdgkG4nL8jcHYhE3ohhz4w2mF3dKlSx2/J64Mv+g9Tb40tV+G/ELll6N2XV66yZ3Mfj1yXe0gw3VZVGy6LdP7cF1TQ0x2jmm3mf8qZZrlrI+TIvOkyYHDwj4ERpqse8Jw8LK2D8Gx5VefOWY4aJq/BmpdD0O9jQbtCHigM318bX3aDDQx+aXPA1KRaYTB5D5uOqBpx/LL6UcMEQlL+0Ciupmsy9qgTBvr9iiPGHC7tv53Ud3Li7G5LovCuQ63qV7rs+d5GtuvPGJBIUFbhjKT203vQ0d7zT8rZQ9weq/J7ai4btvzy98TdMvn2CPzx9fuxwhNUJThfOJ24OgKK+uWGYQHPb3UutuAffPLb1cLRZLecL771QbfL3JiE7D587O3m653dt1+twEtztZDHV8PLH2xfDtcTNcdfC/Q6aLydX+bbrKuyfZ5ecTDQJ9phnWPrQM+u8B6qpCCatiD5f+LD4bCKkPuB86baTifl2YQa9agaNOEHT8Te/+0vm62STE3P1Pmos7DB/D0rShCtevbTwC8/MpvNz01jZbxIDl1vmFdvp+8Ag3fS/xusMRFrjXOqSHgl5uLpLvuxsCsLHgXFuJA+3Y13mZdRChZlmVtdrwtovKjMDRpKFY1W4V0n3R8lv8ZxuSOQaS/yfGnAbAL7TAQmxGLJIzEavyBMXXyuA2zm6Gumd0F8LbwC6PHtcAlBg8lJThesVFjyAPo5XMN5/mB5C93a7QbB0z5rvzynB62i5JvNDmwvD/YUCtjCdbP3LK4/PInYy2P0yGsiblzTfnlLycDqfstr8sOtXu3l1/+/jqDSLAE620ePFh++Zc7gGOrraeFHj1ZfvnvGcDB/yyvyyLdp86UX/z3SWBPRQPMCjyWBOjORpmXvgRs/9b6uibjhLB6DrDxU+vr3rsTCIkznN/4GbD27PvDEtPXAxFnf8nu/AlY/or1dW9dWr4PBxYCi562vi7/d5qwO7YGWPCo9XUju5kIu+0GAWaNuH7lwi7jiEHYWaMdxzENLhdKh0zed+bkUijDJLJm5X1GWDemwR8J1kQdfxSYHvgolAIiDQJd51F+qp03jVRRHHH/+YNF52kQSBRl2ikjdhpMTU54/ext3oZT0/VNf0jx/fzAvvLbPHytpyN5/dVfw26a94crwMwO/VDZbOdo6jMKdMEFF2D+/PlWb8/z80PAxAuR9dXX6LVlixJ3O7t2MRocW6O+I1sBNejo5fixEYkGI+NSz1J48X3vhGnzmuGGBRiGm/A9emKXmkxxug4GB4uwc1rMPrAVPsBmt2kHIPPbeB/j9SYHMC6WtsUDjykqAuBveR94myk8KLEo2tJTYDrHFF42syEwrmx+PWty/JpYfg0YgTCF9+XBucJm3Sy/Zr5hQFCM5dfM/DpGPEPMZyNbug9/eocbopMUDlyHt5meN/1/UCxQOKn7W1jXNKJCcRw/+OztOHddRmE0KMQo/ivsg8ljmEZ7w9sAXS+3sr9n91GjSTtDNEy7ncKJ57nwPLel0awLMPJxgwBhWs+47tnz/BGhwdfgkvfP3Z62rmlhPaOjN/5ltk2T9U2FErd7/14zsaYJNjPBFNEOmGFneUlAU2CKDaFvCj8nfW+yb13uk+n+N0LuueceNb2Ik4wcTX2KBQpVW8KOn7vwBx/Eyu070G3HDnTevVuJu829exlXmTJlCvbs2VOtCJmzElQcpIyM23dojzCf6tUaOrsIPI4Y1TzB7th01M30DBF2joARG2NXrJUDPoXPkxmWbzP/1cUveNPoUmWYRq0q44G99q97d0VbAZvcusT+dU0jiJVh7wGUTP7Y/nUvmmP/umOfNyz2wPSeVjNVGYPuMiz20Os6w2IPnS8xLPbQZrRhsTfqY2/kh/VoWk1aZbAAv8fV9q3LCFiLIfat6+FVHpkUXCZi11jrsRl929O5E4q8vdB74ya0OXQI3kWFwCWGz3K7du3U4szCztvbGyUlVhpvrOCn90OwW3lQ4JcDv8C7xNvpopM14fc6SsFqVL3PWDgXRgKMi0f5wiiChorW6MoXLYriwm9WQRAEwbEcatMGqwcNgl6nQ9zxE9D99pvVdetb7Jg//kUXna2FrSabkjfhqTVP4dENjyIhIAENER1Klc9d7T6GIAiCIDRw6lsEVYUTzeOwYthQnIqOQlkNxVJdEhoaarcNmiW6RXTD+JbjoS/TY2PERuwPslFP64K0RgLuwBfogd3OK+zohPzpp5/iwgsvVJMounbtqhT7vHnznD7vLQiCIDQe6vOYZI+oNF8nOTISK4axM97buP8l6eXd+TTo7d27+hMNqgNTwbWJp84Tzw95Hle0MkwT2RG+AztCdyjvu4ZABNIQgXSMxip4w/rUkHoTdnyTUcTdfPPNOHnypBJ1nTt3RkJCAm688UZMmlS7Q5WdbaRYp06d8M47NjocBUEQBMFEmN199912r3/69ddxZNKlGN6qlWrGaNGiBQabjFBr2rQpapu6EJI6Nx3u7HQnuqQbOuw5pWJTk03Qw0qHuwuxAd2RilAEIA/DsdbukWLUFh9++GHtN0/MnTsXy5cvx6JFizBy5MgKty1evBiXXHKJitxdf71h/l9DRkaKCYIg1C+BgSad4S6QiuX+2DtvvTQ3F9mLl6AkORkx77yLAR9+oO5v+pyYOeOQe5vdt7Ug9KgDHA2fV/vM9srIWJtScbjwMFwdPdzxD0bgWvyC/tiKLapbttyk2lEjxaodsfvmm2/w6KOPniPqtGG7Dz/8cJ0N3BUEQRAaJ1dffbWKXLEcqKGmhnX+/oj/8gv4dOsGfWYmEqZOQ87KVRWEHTtS+/evW7/C4OBgPPbYY7U6pWJg8kC0PdMWbUtteLu6EAfRAnvRGu4oxQWgm4Tj08zVFnbbt2/HuHHjbBoybttmxYRWEARBEBwA05Jjxoxxuoico/EIDUX8Z5/Cf/BglOXl4fgddyDr7/l1WlNo6TX28Khd1zROqeiW0Q1uZ63E8orzkJybDFdmAYahBO5ohePoCBND/voWdunp6WjWzPrsP95Go0lBEARBEGoOI3dx772LoPEXAMXFODVjBlofcLwwcFaK9cW4b+l9uHb+tTh8xnVTsxkIwUr0Uefb4bDzCDu9Xm9TqbMwtKpGhYIgCIIgWI+SuXl5IfqVVxBy9VVqLF6Jp+E43BicKDKLMpGYm4ik3CRc/8/12HbadbOCq9AX3+Ai/IbzHb7tasdQ+SZi9yvz+pZgoZ8gCIIgNHb8/U1GMzoAN3d3RD75JEIuvhjf/fqr1fVatmyJI0eOOPSxbe5XLafDm/g2wbxx8zB90XRsT92OmxfcjCGeQxAMs3GWLkAxPLEPrWtl29WO2N1www2qvZrFk5YW3tYYOmIFQRAEwRJt27ZVjYS1UYdGEeXbo4fxsm9BAZJfehk6vb7W7Enqs45Re+wQnxB8dP5HGBIzBAX6AiwKXOTyUyp8UYCB2OSwRopqv9s+++wzh+yAIAiCINQ1zDYxs9SqVatafQwfHx+71582bZoy/a8Kt956KwoLCpD+wAwU7NqFIZGRWDWk3N/OFalMQPp5+mHOqDl4ctWT+PPwn2pKRd9BfZGyMAWuhjtKcBu+RAiykQM/7EDHGm9TRooJgiAIjY4777xTeb+NGDGiTh6PViQdO9o+aDdv3rzSdcyJjo5Gy1atEHHvvXDz9UVUUhJGLFkKj/x8uCohISEqjdyhQ4dKp1Rc3+l6+Hv6Y2T7c63XzLnyyivhbOjhgU3ops6PwzIVvas3YbdmzRr8+eefFa6jITH/GUzDql8RUmcnCIIgOCEsGerTpw88PT3r5PFoAVabwiJg6BBlh1Lo5YUmaWno/M23QFoaXDVix3Kvq666qtIpFQ/2fRC/XPQLOoZXLog72BCK9clq9EYKwuCPfJyHFfUn7J555hns2rXLeHnHjh246aabcN5556magj/++AMvvvhijXdQEARBEBoi1ZkhawvW3C0ePQp5vr7wS02F7vEnEJCVjYZOVECU8XyadxrWNl2LErcSl5pI8SfOU+d7Yyea40T9CLutW7di9OjRxsvffvutCjV/9NFHuP/++zFnzhx8//33Ndo5QRAEQWhM1NS2JCs4GIvOG4380FC4nT6Nvhs2KFsUV8cegct5suuarsNJ/5NYGbkSRboiq+tec801aNOmDZyFY4jBRnRV5ydikaq9q3NhR/NhU4PiZcuWqVCzRt++fXH8+HE4K5MmTUJoaCguu+yyc25jiplu5uxo+vjjj+tl/wRBEITGhyM6T/P8/bFjytUo69cXawcM4EYdNnO3KrRr1w4333yzw7ZXGe5wR7+UfvCBD9J80rA8ejny3fOt1vE5G/9hiGqgiEC6miVb58KOok7zxykqKsLmzZsxgG+gs2RnZ9dZ7UJ1+N///qdqAs2hqTIjjosXL8aWLVvwyiuvIM1F6xQEQRAE14Lj0Rzhe1fi54fShx5Cvr+f8Tq/3NxajRiybtGUKVOmIDY2FnVJk8ImuNbjWkT4RiDTMxNLo5Yi28M10tEF8ME/GIEt6IzN6FL3wm78+PGqlm7FihV45JFH4Ofnh6FDh1aYJdu6de2Y7zkCdkJZ+mWyfv16dO7cGTExMQgICFBRyIULF9bLPgqCIAiNA01QhYWFYcaMGQ7ffu/CQoz/62+0PHTYodFEXvfQQw/hgQcegJeXF5yBprqm+GL8FwhzC0OeZx6WRS/DoZxDcAV2or2aRkGRV+fC7tlnn1Wmi8OHD1d1dR9++GGFfyq9eM4/v3qjMpYvX46JEyeqNm6+aX614Kz9zjvvoEWLFsojiLV9FGSO4NSpU0rUafD8yZMnHbJtQRAEQaiKgGLJkCMIPXkK7qWl6LdhAzrs2ePQujsGdhyZwnUEMQExmOY7DSGFISh0L8SfiX86jdmy/ZQhAql1Z1DcpEkTJcAyMzNVZIuzYU354Ycf1PXVITc3F927d1dmjZdeeuk5t3/33XcqXfr+++8rUTd79myMHTsW+/btU1YrpEePHhZn1TL6RsEoCIIgCM4KXSZ4fI2KKu/4rAmJF05AcnY2Ou3Zg+7btsO7oBDbenSvUv1dfcyj7dKlC3bu3IkhQ4ZU+b4BugAMSxyGvSF7cceAO5SY69atG/Lz8xEeHg5nxgtFuBx/oSWO431ci5Pwr31hV1paqurPfv/9d1Vjxw7ZmTNnwtfX1xhOri5Mf5o2Ypjz+uuv45ZbbsHUqVPVZQq8v/76S0UJmR7WunarA0WfaYSO5/v162fzPllZWee4jVuboSsIgiAIlREXF6eWmmAqXtx0Ouzo3g1F3l7osXUbOuzbB++iQmzo2xdlOvuSd/UR5Zo8ebLK/gUFBVXr/p5lnuia0RVeOi9j4+TG5I1OH7HLK3FDmU4PD50eE0v/wZzC8bWfin3++efx6KOPqqgc05Vvvvkmpk+fjtqGInLTpk3KL09Dp9OpyzRNrikUcfx1QEGXk5OD+fPnq2igLfjhM52TK/59giAIQmWYli+ZZ70cAaN9V1xxhRoYoEXb9nXogHX9+6HUzQ0tjxzF4FWr4FZaCmeFAqy6os4Sb2x+A9MWTMNH2z+qlwikvaxcuQqj39yNrMIyxOuSkb/i3dqP2LGj9N1338Vtt92mLv/333+YMGGCsgeh0KotUlNTodfrK1itEF7eu3ev3duhENy2bZtK+7Jrh6njgQMHqrrB1157DSNHjlRRSRaFVhaypa2L6RtPonWCIAiCPbVpl1xyiRJ1PPbUBp06dTrnuqMtW6LIywuDVq1GTkAAypw8euUoysrK4OFmeJ3nbJmD3m69EY94uMH5nj9Tz9Qk/+l24VIswzOjA/Dexgy77lvtd9KxY8dUZ6ypUKKyZvNBXbc3VwcKUWtcdNFFarEXijpH/qIQBEEQGgesB68LzFOPp2JisHDs+cjiscvkNmeOYjniNbin1z0I8wnDyxtexqayTUiJSEGf032gq34Cs1ag0OeyHT3RA4fRxO2Y3fet9jNhYwI7Uk2hb11xcTFqEzZt8NdNcnJyhet5OTIyEvUBzZj5q4iduoIgCILgCnBKhSbqdHo9+q5fDz+zmvGGyLWdrsWLQ1+Eu5s7jgccx+pmq512BNn69Rtw2bwkfLjNfrlW7YgdVf2NN95YIe1YUFCA22+/vYK54s8//wxH1yT07t0bixYtUiFswpQpL991112oDzZs2CARO0EQBMFl6b5tG1odPoLoxCQsHT4MmU44mcGRXNjqQpTlleGJDU8g2S8Zq5qtwrCkYU6XllXNm/36Ia2wEJj/Uu0KuxtuuOGc66699lo4AjYtHDx40HiZEy7Y5cpO2+bNmyurEz5+nz591JOm3Qlr5bQuWUEQBEFozFQ1pbqnY0c0TU5BSGYmRi1ajBXDhiI1IgINmX4R/TA0aaiK2LXKbuV0oq66VFvY0WWa/jK1wcaNG1XzggaFHKGYmzt3Lq688kqcPn0aTz75JJKSklSNwj///HNOQ4UgCIIgCJVT4OuLxaNHYcTKVQhLScHwpcuwevAgJJr4vjq7RUh1CC8Mx9gTY+FVWt6hXIYylxZ51a6xo8kfZ8Ny6gTnwjp63Bd/bZgvFHUaTLsmJCSgsLAQ69atU0bF9YXU2AmCIAiuTrGXF9ZeMA6noqLgoddjyIqViD961OK6119/fZ3vX00It+FuYSrq8t3zsSRqCTK87OtArW04VYvagtO9al3YLVu2TIkZRu7olcNoGufGNkZYY7d79+468fETBEEQBEfw4IMPqpImGv5r6N3dsXLoEByNj4eurAw9N2+BF+u7zNCmPDkjbhYii/YOTdgZuhMZPhlYHrUcKT4pqG9YbkZtQS/CWhd2Q4cOVZMeEhMT8dZbb+Ho0aNqbmy7du3w8ssvqxSpIAiCIAjn4gy2IvTRu/DCCyvMRyecRLFuQH/s6dABK4cMQZGDvVlNGyydjR5pPRCRH4ESXQlWRa7CCb8TcDV0jvgHsWmBEbz9+/fj8ssvV2FDNjlUxQvOlZFUrCAIgnNTkzGXjRI3N2zv0R2pTcsbKMpOnHDqKRWOwLPME4OTBiMmNwalbqVY13QdDgcedqlUrEOtrtu0aaPGjMXHx+ORRx5R81sbA2J3IgiC4Nwwo0Sf1c6dO9f3rrgkIRkZKL7/AQwIC1PRPGeIONYW7nBH/5T+2BK+BUeCjmBLky0odC9EhzMd6rypgqlYLuwneOkl++xOHGa1vHz5cuVrR5Ng5u0vvfRSrFq1ylGbFwRBEIRqQ89Vjr1s0aJFfe+KS9ImKBgoLETz48cxdPkKlObloSHjBjf0TOuJDhkd1OXj/sehd9PDFahRxI7jw9ipyoW+c4MGDcKcOXPU0GFnzqELgiAIQn3irNYhliJxHAZAW7H0USNx4u57EJmcjNN33An/Tz6GR2goGrK463ymM/xK/NAsvxk8ympnnq+jqXbE7oILLlApVzZOTJo0CXv27MHKlStVvZ2IOkEQBEGwjqukMhnl1ObZ+g4YgKUjR6DQywvFu3cj4ZprUZyYiIZOy5yW8NP7GS+f8j3ltCPISLXlJ+fC/vjjj6qjhrNbGzNsnuBrwJZksTwRBEEQGqpwTA8Px6LRozFx0yYUHT6Mo1OuQfwXX8ArtmJnbUONmp7wP4F1EesQVhiGQcmD4F3q2I5hS80TrOPn6NRaF3a///57de/a4JDmCUEQBMGchlrPlx0chIhPPkb6Pf+DR2QzeJh0zloLBLFxpSHgW+ILz1JPpPukY1nUMgxJGlIhmtegmicEQRAEQQDuvfdeTJ48Gb169YIzR6dqUufnERmJ+K+/Quxbb0PnVT65wRI333xzg+lGDi8Mx4jEEUrgZXtlY1n0MmR7Onb6Vk0RYScIgiAIDiQkJARdu3aFTtewD7FsnHAP8DemflNeew1RR84dQcY57vS4bSgEFQdheOJwBBQFIM8jD0ujliLdKx3OQsN+1wmCIAhCIyI6Olqdag0PdUX2v/8i7aOP0XvJErQ6eAgNHf8SfxW5Cy0MRZF7kRpBluuRC2fANXp3nRxpnhAEQRDsgYJr69atGDZsWK1sn36yycnJiI2Ndfi2baVuA0ePRsiVV+LMd9+h78aN8C4sVFE8Z7V1cQRsnBiaOBRrm61FcFEwmvk2Q052jus2TwjlSPOEIAhC3RIRYbtg31nhqM0BAwao9GRt4OXlhbi4OId109orEN3c3RH51Eys37sX7bZtQ7cdO5Dy0kto+n//B7cGnJL2PDuCjJ53boEGEVuKUugclBCV5glBEAShQcNC/KFDh6rFFWHdHSc0uUIka9q0aYiKirJ7fT6nfb17YXPPnupy+ufzcOrhh1FWTx2xzWpJPHfr1q3CZYo4bdQYp1OsjFyJPSF7UIb68SqUiJ0gCILgMjCCVBtpRuFcmjdvXi2fvAPt26HI2wsDNmxE1h9/ImTyZfDv3w91zYABA6DX69Uce0diS5Qn+iXitO9ptRTqCtE9vXudz5cVYScIgiAIgkNJaNECl914I0qSU+pF1BEPDw8MHz4cdUlsbiwKUguwLXwbDgUfQqF7Ifqe7lshNcsoaGItTuyQVKwgCIIgCA6fWBE4YgRCr7zCeNk3Lw8++flo6LTJboN+p/vBrcwNJwJOYFWzVRVGkDGSWJtIxM4BSFesIAiCIFinJCMDw5cug7tej/XjL0BuLY4idbOzfjE4ONh4Pj4+Hunp6cjOdozZcFxunJpQsSFqA1L8UrA8cjkGJw+u8ggy6YqtJ6QrVhAEQRCsU5qbB11pKQJyczHoz7+wbNhQnAkNVbeFhYXVS9TR29sb//vf/1RgJjAwUN3vmWeesbjuAw88gNdee61KjxuZH4k7m9yJd5PeRY5nDgrcC6os7KQrVhAEQRCEWqEmnbxesTFYPHoUMkJCVDp25OIliEhJUUGR+uxwDg0NVfvA58aO5VtvvdXiehR/1XkdWni3UFMqGK0LLjZECMPDw1GbSMSuDmBXTkMZgCw0HDiY2/TLShAEwVa0q6b1dwW+vlgyaiSGr16N8KRkDFu2HB79+6vvImeb3FFVLr74YixatMjqCDJTkj2T0X1cd2z7ZxtqAxF2tQg/BElJSThz5kx974ogWJ1p6SqeWoIguD7FXl5YP24cui38FzGnTgEvz0JWVBSCzj8froKb2fclL/fs2dOqsDMlyzMLdy+6GyVlJejt0xvNChzvtSfCrhbRRF3Tpk3h5+cnB0/BqX505OXlISUlRV2uigmpIAgNl7o4TpV6eGDVkMHou34DWublwbd73c61rQ1bFWv4+PhUuOxX4ocuTbpgTeIarIpcpbpnaZHi0P1x6NaECulXTdTVdj5dEKqDr6+vOqW44/tU0rKC0HCxV7DZSrc6UvSV6XRY378fWo4bB89mTeFMuLm5Ocz2hY0PR48exYEDB9RljzIPvD36bTy68lEsOLoA6yLWKSPj1tmt4SikeaKW0GrqGKkTBGdFe39KDaggCHUOhWJIiPFi5h9/IvHJmSjT6+FKuNkQvKwfvOaaaypc5+XuhZeHvoxWWa3AoRRbm2zF7pDdDhtBJhG7Wvaxk/Sr4MzI+1MQGgeONB6uDYqTU5D42GMoKyqCPiMD0a++Ap131axB6opwB2Th3HXu6JHWA956b+wJ3aMWpmlb5LSosY+dROwcAF/03bt3izmxIAiC0ChwtFBkOpZizs3TE9n//ovjt94GfU5OvfxgLavkuXXv3h2OgDNkO53phB6pPdA0rynicuIspnKpLazZsFhChJ1Q57DegB+8rVu3qstLly5Vl7Xu4blz56puzZriqO0IgiC4Oq4QnWdnbNxHH0Ln54e8detw7IYbUZKWBmdDp9M5tOGM9XVDkofAHYY6Z6ZkTUeQVXn/HLZnQoPi+PHjmDZtmvL08fLyUuNW6NCd5oAPWVxcnBqA3KVLFziKFi1aYPbs2RWuu/LKK7F//36HPYYgCIJQu/gPGIDm8+bBPTQUBbt2IeGaa1F88iScGS8vL3U6aNAgdVodXz5G7zRRtz1sO1ZErlBNFdVBhJ1wDocPH0afPn1UF88333yDgwcP4v3331cePQMHDlTz9GoC6xHpnWarRdxRXZ/s9hQEQRDqJsXK7/aaPoZvl86I//oreERHoejoUdVU4cxMmTJFnfL4eNttt+Gyyy6r9rY4duxYwDGk+6RjWdQy5LnnVXkbIuyEc2A+n79AFi5ciOHDh6N58+a44IIL8N9//+HkyZN47LHH1Hrvvvsu2rZtq3x6mjVrVuHNzELPWbNmoU2bNmoeH7fx/PPPW0zFVsahQ4eUqzcfIyAgQDWrcF80RowYgYSEBNx3331qu1rKwVIq9r333kPr1q3V82vfvj2++OKLCrfzvh9//DEmTZqkOkb5/H7//Xfj7RkZGarDKSIiQglH3v7ZZ59V63UWBEFoaPC7k4GBO+64o0bb8W7ZEi2++QZN7r4L4bfZX19WV/AYYD6tgscPpmhrYh3lq/dVI8h8S3yR7ZWNtS3XKlPjqiDCrq4pyrW+FBdUYd18+9atIozGLViwAHfeeafR58z0lxhFzXfffacaRu655x41MHnfvn34559/MGzYMOO6jzzyiBpY/MQTT6jGkq+//loJs+qQk5OD8ePHq4jhli1bMG7cOEycOBHHjh1Tt//888+IjY1V+8IULxdL/PLLLyqdzGHOO3fuVL+spk6diiVLllRY7+mnn8YVV1yB7du3q8flc9ailNrzmT9/Pvbs2aOEYpMmTar1vARBEBpanV5gYCAuvPBCm9/39tb7eTZrhojp043rlxYWIm/zFjgD7rXo+8kRZCNOjUBAUQAy9BkqcpfhmWH3/cXupK55wcYcurbnA9f8UH75lTZAsZUwbPwQYOpf5ZdndwXyLNS/PZVZpd1j+pVh8o4dO1q8ndczanXkyBH4+/urDzA/yKzB40gVkp2djTfffBNvv/02brjhBnUdo2RDhgxBdWAHkmkX0rPPPqtEGiNpd911F8LCwtSHjPthKw3w6quv4sYbb1Sildx///1Yu3atun7kyJHG9bjO1Vdfrc6/8MILmDNnjmo5p6CkmOTz5C9SrbZPEARBqF3KSkpw8v4HkLN8OaJfehHBEyagIeOn98OIxBE40usIdqXtUlMq7EUidkK1aiD69++vxFyrVq1w3XXX4auvvlIjqggjWYWFhRg9erRD9oURuxkzZihRydQq07F8DC1iZy+8z+DBgytcx8u83pRu3boZz1O8BgUFGUdvMb3w7bffokePHnjooYewevXqGj03QRAEZ8Jp/e5KS6Hz9qKbOk7NeBDpX3+Nho53qTc+GfsJIvIjqmReLBG7uubRU9ZvczML7T540Ma6Zpr83h1wBKyJY9ibYoe1Eubw+tDQUFUzt3nzZmVVwlq8J598Ek899ZRK0ZqncGsKRd2///6rImvcP26f9XxFRUWoDcw7mvh6aOaQrDVkPd/ff/+t9onilTWJ3DdBEITGQl0LQDcvL0S/8grcQ0KQ8fU3SH7mWWVk3OTOO53OysXDgY2B/p7+GJw0GKfLTmMXdtl1H4nY1TVe/tYXT58qrGsmnqytVw1H7TFjxqjGiPz8inV8SUlJKjJHGxF+kPjmPe+881STBOvR2BSxePFi1VBA8cWaOEewatUqlR6l0OzatatKt/KxTGEzBOfz2oIRP27LfNudOnWqctEsU8xffvmlslj58MMPq3R/QRAEoeq4ubuj2RNPoMnZYQCpb72N5OeeR1kVpjI4Cpb+WIPZrA4dOlS7/Mgc+tuFFoXavb5E7Gp5pJgrwto4+vGMHTsWzz33HFq2bIldu3bhwQcfRExMjOpu/fPPP5UtChsmGMFjBItRLXaaskv2//7v/1SqkoKL6c7Tp0+rbdx0001V3h8KRTZIsGGCgpINDObjVVjrtnz5clx11VWqC9dSQwP3n00RrJGjIP3jjz/Udk07bCuDkcnevXujc+fOKt3M18FaPaIgCEJDpb6iZHzciLvvUj53yc8/j4yvvkJZqR5RM2c67DGGDh2KFStW2IxQUrSx3pzHAkv7yGORI3jnnXfUUpWRYiLsHADTj6zDaihQSG3cuBEzZ85UQogdoYySXXLJJeo6Niuw1o2iiOnXgoICdR963mlvcoovRvQohE6dOqVawG+//fZq7c/rr7+uzJIpNinYKBqzsiq2f7Mjll2ubNKg4LKUJuD+s6mDaVN2x1Kw0qqEdin2QqHKjl9GDBmV5BcAa+4EQRCEuiPs2mtUWjbpmWcQcvHFDt32yJEjVcTto48+snksmDx5MmobBosYGOFxjU4T9iDCTrAaSqYPnDX4a4X1dbZGrtDvTvO8M4+umQovCivTy0y7cjFdnyleU8wjowMGDMC2bdsqXGe+Ha35wZa/kiVBqI06I48//rhaBEEQGjOOqrFj+U91Cb5wAgKGDYW7aWClrOb7xeMXs1Ouigg7QRAEQRAUvXr1UvZObJCrTZhhoTVWTacDmYq60PR09Nq0GSmxsWjMiLATBEEQhEYEU4isRbYEy25Yj+xodwNzWJ7DxZERxL4bNiA04wyCPvkEheedp6ZXNEakK1YQBEEQGhF0F2jXrp3V2+nfyXSkK8GGhZWDhyA7MABeZzKRcM21yN9lnz1IQ8O1/nOCIAiCIFQZZ/N6qw3yAvyxaPRo5EdFQp+ejmPX34DcdevR2BBhJwiCIAhCgxCHhT4+OHTjjfDr2xelubk4fsstyK6CpVVDQISdIAiCIAgNhlIfH8R9/BECzhuNsqIinPnhR+cdlVYLSPOEIAiCIAgNKrKo8/ZG7OzZSPv0M+V55wrRRkchETtBEARBEKqEK0TA3Dw80OTWW6Dz8zPuM9OyrrDvNUGEnSAIgiAIDZ7Ud97FibvuRuJjj6OspAQNFRF2gkvCkV4MrW/dutUptiMIgiA4N570zdPpkPnzzzjxv3tRWlBgc/1OnTrBFWnUwm7SpElqgP1ll11WpdsaA0lJSWqeaps2beDj44NmzZph8ODBeO+995CXlwdXhOPFOC/WlLi4OCQmJqJLly71tl+CIAhC7RMy+VLEvjUHbl5eyFm0CMdvuRX67Gybfn6uSKMWdhQu8+bNq/JtDZ3Dhw+jZ8+eWLhwIV544QVs2bIFa9aswUMPPYQ///wT/1lpHS8uLoar4e7urpzWPTykj0gQBMHT09Pqbabfk4GBgXW0R44lcPRo1TGr8/dH3oYNSLj+BpSkplpct3379uqUwQ1XolELOw6ft/bmtHVbQ+fOO+9UH+CNGzfiiiuuQMeOHdGqVStcfPHF+OuvvzBx4kS1HlOYjOBddNFF6pfN888/D71ej5tuugktW7ZUI2n4wXjzzTctRs5effVVNVKGQ6CnT59eQRhy27/++muF+4WEhGDu3LkW97myx33qqafw+eef47ffflPb5rJ06VKLqdhly5ahX79+auQO9+/hhx9GiUk9Bt8b99xzjxK6YWFhShhy+4IgCK7K+eefjw4dOthMP/K78tFHH8Ujjzzi0j+G/fv1Q/wX8+AeHo7CPXuQcN31KC0qOme91q1b4+abb1aBHlfCaYXd8uXLlYCIjo62eJAn77zzDlq0aKHUdP/+/dXgYmcnrzjP6lKoL7R73YKSirUB1tarKmlpaSpSR6FlLQxt2jZOQcO09Y4dOzBt2jSUlpYiNjYWP/zwA3bv3o0nn3xSfRF8//33FbaxZMkSHDp0SJ1ScFGwWRNt9lDZ486YMUOJ1HHjxqnUK5dBgwads52TJ09i/Pjx6Nu3L7Zt26aE6yeffILnnnuuwnrcZ74+69atw6xZs/DMM8/g33//rfb+C4Ig1CY8ltqC34dXXXWVymLYwsvLy+qcWVfCp1MntPjqS3jGxiL8pmnQeXlZPNbxuFLbc3MdjdNK7tzcXHTv3l2JhUsvvfSc27/77jvcf//9eP//27sT8Jju9Q/gbxJJSGKJXWQrRdUWS6naiYsq5dZaf0LboNfWWqr+LZ5y23DVGje1t2jvn7rKRVVVKnob1L1cpZcKHooKsaSJkIgs/+f7a890ZpJJZiazZfL9PM885pw5c85vxiFv3ve3rFqlgrply5ZJr1695Ny5c1KzZk11TEREhEGmRYPApbib3F7a/a2dydc61e0kcZFxuu2un3aVzJzMQo9tU6uNfNj7Q9127+29JfVhaoHjTkedtqh9Fy5cUEPBtRS0pnr16pL1W0dTBH0LFy5Uz1988UUZM2aMwbHvvPOO7jkyaCjjIsBCYKVB/8WVK1eq/0TwW2Lfvn0lPj5eoqOjxdryQVHXDQgIUP84Hz58qDJspsTFxal+d2gb/lGjbdevX5eZM2eqYFFbP7F58+Yyd+5c9bxBgwbqeLS/Z8+eVrWfiMie8MsqIBniLKhw3L17t8DPF2fxCQ+Xert3iade4Jaflycedlgnd/jw4apiVFg8U2YCuz59+qiHKUuWLFFBgBZUIMBDmXDDhg2qdAYc6Wg7yIYiKzZixAgVHGnatGlTaCYVfw9XrlyRzMxMyc7OVkG2viZNmhj8ZoiSJ7J+JWHOdYtz9uxZad++vUFWEoNGMjIy5Nq1axIaGqoL7PSh/SkpKSVqPxGRveCXUiRBnAmJmqSkJJcarOapF9Tl3LkjV16JlpqvvyYBnTsbHKffz86aPncIZmfMmOGQiZJdNrArCn5gHz9+XNX59W/ayMhIlaVxtPT0dIPtotLU3734ncnXvDwNU+AJQxJMHuvpYfgbxb4X9oktYBQsbjxkPvWhjx0Yp6SNy7VbtmxRZc/FixerAAn9FBctWqRKlkV10MU1ETjqbxtPIlnU4Axzr2srxbWfiIgMoXLSqlUrcVV31q1Xfe6u/mmCBMW8J5V/608OSESgX7X23BqWBHWoNuJnIOILJFP0EypuGdjdvn1bdZbHFBz6sP3jjz+afR4EguhHhbKv1j8LQUFxrxlD6U4fSnRa1tCYn/evM2Cbw17HFgUDGVBORGlx0qRJFg/3TkxMVH01MABDg750lqpRo4bqB6c5f/58kdOsmHNd9A3BfVMUDBTZvn27+gel/SPEuREo4j4gIiL3VHPq62qEbPru3XJ9xhuSm/qLVB01Uve6328rWDjCt99+qwbyYWYKS5XKwM5WTE3bUdxrxq5evSqVKlUyyNiV5iVL0M8M5UeUWTE4AmVHZET/9a9/qcC5devWJt+L/maYJubLL79U/dw2b96s3ofnlujevbsKLhFMIxhDH7eihuGbc130LcHryEYigK1cuXKB8yAwRH9NBLUTJ05UxyJQR39OrX8dERG5Hw9vbwlauEC8qlSR1M2b5eZ770nuL6lSfdIku5ZQ8bPFuOLTsWNH2bVrl3oeExOjMnZLly4173xSCqEjP1KhN2/eNNiP7aI6xtsLgjr9R2kfMYQh3pi7DllLlLsxiAVBXmxsrCp3zp8/3+R7x40bpzqHDh06VPXnwChb/SyauVBSRSa0U6dOaoAGrlvUb0vmXBd9MtHPAZ8FGUFk4ozVrVtX9u7dq/oU4nOPHz9eTaPy9ttvW/wZiIiodPHw9JRa/ztLakyZrLZvx30gN955R/KLqfZYA9Wxp59+utDBnJhORj+msCSu8MgvBaklRMo7duwwWDUAP7wx1xiCDUC0i47tyLKYKoPaGmrfyPo0bNhQBZoYLYoHYATppUuXVMaotE1uSGUH71MichfafKJPPvmkwSwM1kr9v/+TG/Pmi3dIiIRv3SLlAgPFHtatW6cG5xX2eTAoEH8ixsGI4rS0NIMKYakqxWIUIqbe0OCHD0a5Yrg0AjiUxqKiolT2BQEeymfoD2c89YYjoORX3BdNREREpUfg8OFqEuPyjRvbLagrDpJFt27dUqXYBQsWmPUelw3ssOpBt27ddNsI5ADBHCayRckNHxZzi2FdU0xrsW/fvgIDKoiIiIisUekPfzDYvvf1QanQMsJpgV6pDuywbFNxVWKUXfFwhYkfjUuxRERE5DionKWnpxe5LFpJ3Dt4UK5NmqQmNg5dt1a869QRe0MpFoMB3X66E1fDUiwREZFzvfrqq6qSZzwFma34hIZKuRo1JPviRbn84ggJXb9OfH+b49VetISR1qffbUfFEhEREenDBPqhoaF2m5rEt359Cf/bJ+Lz2GOSk5wsP704QjJLuGKSPTCwIyIiIjKDd1CQhH3ysZRv2lRyf/lFrkSNlvuHD0tJYGoTW2JgZ6M+dqjpoxZORERE7qscZuf46CPxa/+05D14IFfHjZespCSrz9dPb+kyY4grEF8gzjC7fVa3hHTYx46IiKjs8Arwl5DVq+X69BniFRgovg0aWH0urIRkyz52DOyIiIiILOTp4yN1ly5Rz7V+ffmPHqG2atclyIptl9OuTOSmRo8ebbBKirPPQ0RE9uHh5aUeWlB3deJESVmwUPKN1n51JAZ2NuBufewwXBzDxjG6COvTYf3dXr16Fbq2KpXc5cuX1W93WFlF3/Lly9Vk3ERE5PruHz0q9w99I3c3bpTkWbN+zd6VEPvYOYm79bF74YUXJDs7WzZu3Cj16tWTmzdvSnx8vNy5c8ep7UKbfHx8HPY+ZzO3PwURETlfQKdOUmdBjCS/9bak/WOX5KalS91lS8WzBOtwcx47KrFffvlF/vnPf8rChQvVkm5hYWFqLd5Zs2ZJ//79C80u4T3Yl5CQoLbxJ7Y///xzad68uVpc/umnn5YffvjB4FrffvutdOrUSc09hAklJ0+erNb71YSHh8v8+fNl1KhRKnAeO3as2r927Vp1vJ+fnwwcOFCWLFkiVapU0b0PCyZjiTksrKy/uD2WnOvYsaM6Fp1Vn3vuObl48aLufdpn++yzz9Rnx/lbtGghR44cKXBufVinGG01pbjroo3QsmVLdX2sulJYKRYzj+M7qlmzpvpMOCd+qdBo3zuCcKyhjPY/88wzcu7cuWL+1omIyBaqDBggwStjxcPXVzISEuTKK69Ibnq6OBIDOwfD0GiTD6MlQ4o8NivLrGMtFRAQoB47d+60aAmTwsyYMUMWL16sgo8aNWqoId2PfktNI7Dp3bu3yg6eOnVKtm7dqgI94yXi3n//fRVc/ec//5HZs2ercvD48eNlypQpKrjs2bOnvPvuuwWufeHCBdm+fbsK0rQgFEEj1hzGOsQIfjw9PVVgmGfUF+Ktt96S6dOnq/c1bNhQhg8fLjk5OVZ/D8Vd99ixY+rPAwcOSHJysmpzYd544w31mZBJPXHihDz++OOqRH737t0C7cf3juthfqSXXnrJ6rYTEZFlKnbrJqEb1otnxYqS+e/j8tOoKMm5dUschaVYBzvXqrXJ1/y7dJbQ1at120kdOkp+Zmahx/o99ZSEbd6k277QI1JyU1MLHNf4x7MWtQ+BAPp1RUdHy6pVq6RVq1bSpUsXGTZsmMq+WWLu3Lkq8AIEI8HBwbJjxw4ZMmSIxMTEyIgRI+S1115Trzdo0EBWrFihrvXBBx/osmzdu3eXadOmGQQtffr0UYEXIPA6fPiw7Nmzp0D5ddOmTSqg1CCI1Ldhwwb1+pkzZ6Rp06a6/Th337591fN33nlHmjRpogLFJ554wqLPb+51tTYim4f+jKaCQ3wv+LvB59cyl1999ZWsX79eBdEaBLr4HuHNN99UnyUrK0v3nRIRkX35tW4tYR9vVhm77CtX5NHNFLUcmSMwY2cD7jZ4AoHI9evXZdeuXSqrhhIfAjxLO/K3b99e97xq1arSqFEjOXv210Dz+++/V+fTMoR4IPuELNalS5d070NJUR/KiigN6zPeBpSQ9YM6OH/+vMq+od8gSrta+fTKlSsGx+kHsHV+W+Q5JSXFos9uzXWLggwnsp0dOnTQ7fP29lafXftO7dV+IiKyXPlGjST8b3+TkFUfSIWmTaw4AwdPlIrBE41OHDf94m9DpjUNE781faynYUz+ePwBsSVkd5BtwwMl0FdeeUVl4ND/DvLz83XHauVVS2RkZMi4ceNUnzFjGI2r8ff3t6r9hb0PpWAEfMh0BQUFqSASGTNk9/QhYNJocxFpZVOUUfU/uzmf39zr2kpR7SciIsfxCQlRDw3Wls1/+FD8jJIWpnCC4lLA08/P6cdaA78xoN+dlgVDXzB09gfjaTo0R48e1QVpqampkpSUJI0bN1bbyACiFIl+YpZA1k9/wAAYbxcGI3qR7UNwhQEbgD59lsLnv3HjhgrutKDJ1Oc397raiN3c3FyT56lfv746Dn0MESRqASU+u1bOJiIi15X9009yNXqs5GVmqtGy6ItnDwzsqEAgMnjwYNXhHiW9ihUrqk74f/nLX+T5559XI1gxwnXBggVqNCdKfG+//Xah55o3b57qN1arVi3VN6569eq6UZ4zZ85U58FgCWQDkWFDoIc+YytXrjTZvkmTJknnzp3VSFhkwr7++mv54osvip3lOzAwULVlzZo1qjyJMij6n1kKI1Yxzx++j0GDBqkRr7i+qYytOdfFKFd8rzgX+iEiW2r8mxm+H8wtiL50KGsjYEYbHjx4IC+//LLFn4OIiByrHP6vj4hQo2WvTZwkdd79sxpFa2vsY0cG0NetXbt2snTpUhVAoWSIUiwGU2gBFzr/Y5Ro69atVbboz3/+c6HnQvCH0as4Dlmu3bt367JTCBoPHTqksnjIZCH7N2fOHFWqLAr6mGFQBwI7jJZFMPT6668XOzAAJdQtW7bI8ePH1WfCexYtWmTx94OMY1xcnOr3gOtjRKs2kMPa62LACgaOrF69Wn1+BNCmvk/0fxw5cqTKeGJAx5dffqmCRyIicm2eFSpIcOwKqYz/43NzJfnNWXLHDpPQe+Qbdxgis2k177S0tAIZG4xCxCAA/XnUygoMtsA8cCi/6s8vZy8IOn/88Udd/z8yX1m+T4mInAHLjaX8ZZHc/S2oqzZunMT9kopO0QbHYd5Uc+INY8zY2YC7jYp1dZjbDqNqkbGKjY1VU6lERUU5u1lERETF8vD0lJoz35AaU6eq7TurV0v9C79PWq+Po2KdxN2WFHN1KH+if9m9e/fUFCIoY6KfHhERUWmAfuHVx0aLV5XKkr57j1wO+302CH0cFUsuAQMM7Fnh//TTT+12biIiIkcJHDJEqgwaJLnz5v26Iz9fvHJzJbec9eEZS7FERERETizNapqe/kG6f/21+BotG2oJBnZEREREToZgrv7Fi1L1bqp0j/9aHl2/btV5GNgREREROdnD8uXl6x7d5b6fn1S6d08uD39RHl64YPF5GNgRERERuYB7lSpJfGQPSatUSXJu3pSfRvyPZH7/vUXnYGBHRERE5CIy/fxU5q58i+aSm5YmP415STKOHjX7/QzsbIDz2BEREZGtZPv6StiGDXI3JETyHzyQ2RMnmv1eTndiA5zHrnDh4eFqyTEuUk9ERGQZT39/af/5HrmfmChLW7eWDWbOY8eMHRUwevRoNXki1ibVt3PnTrXfkoB37NixdmghERGR+/P08ZGK3bpZ9h67tYZKNawbunDhQrXeq7Vq1Kghfn5+Nm0XERERmcbAjgoVGRkptWvXlpiYGJPHbN++XZo0aSK+vr6q7Lp48WKD17Fv2bJl6jlWosCCxqGhoer4oKAgmTx5su7Yhw8fyvTp06Vu3bri7+8v7dq1k4SEBDt+QiIiIvfDwM7BHj3MNfnIeZRr/rHZ5h1rLS8vL3nvvfckNjZWrl27VuD148ePy5AhQ2TYsGFy+vRpFbTNnj1bPvroI5NB4NKlS2X16tVy/vx5VdZt1qyZ7vWJEyfKkSNHZMuWLXLq1CkZPHiw9O7dWx1LRETkzrpZWG4tCgdPONiaKYdMvhbWtJo8N7GFbnvDjH9KTnZeoccGNagiA6e10m1veuuwZGU8KnDchFXdrW7rwIEDJSIiQubOnSvr1683eG3JkiXSo0cPFcxBw4YN5cyZM7Jo0SLVR8/YlStXVAYQmUBvb2+VuWvbtq3utQ8//FD9iUweIHu3b98+tR8BJhERkbvq0qWLHDx40CbnYsaOioR+dhs3bpSzZ88a7Md2hw4dDPZhGxm23NyCmUJk4DIzM6VevXoSHR0tO3bskJycHPUaMn54D4LDgIAA3ePQoUNy8eJFO39CIiIi98GMnYONXd7F5GseRmH2S4s6mT7WaHDqqHefEXvo3Lmz9OrVS2bNmlVoJs5cISEhcu7cOTlw4IB89dVX8qc//Ull9xC8ZWRkqNIvyrv4Ux8CPCIiIjIPAzsH8/b1cvqxlsK0JyjJNmrUSLevcePGkpiYaHActpF1Mw7ONBUqVJB+/fqpx4QJE+SJJ55Q2bqWLVuqjF1KSop06mQ6mCUiIqKiMbCjYmGQw4gRI2TFihW6fdOmTVMrbsyfP1+GDh2qBj6sXLlS4uLiCj0HBlUgeMNoV0yB8vHHH6tALywsTKpVq6bOP2rUKDWyFoHerVu3JD4+Xpo3by59+/Z14KclIiIqvdjHjswyb948ycv7fSBHq1at5NNPP1WjWJs2bSpz5sxRx5gq11apUkXWrl2r+uEhWENJdvfu3SqoAwySQGCHgBGZwQEDBqgJjjHIgoiIiMzjkY8Jxsgq6enpUrlyZV35EeVFPCArK0suXbokjz32mJrsl8gV8T4lInINmDbMeBtr0OOBildSUpKkpaUVu4QpS7E2wLViiYiIyNa0hJGWSDIHS7FEREREboKBHREREZGbYGBHRERE5CYY2BERERG5CQZ2RERERG6CgR0RERGRm2BgR0REROQmGNgRERERuZCBAwda/V4GduR0Hh4esnPnTmc3g4iIyOmwElCLFi2sfj8DOzIIsIp6GC93ou/y5cvqmJMnTzq0zURERPQ7LilGOsnJybrnW7dulTlz5si5c+d0+wICApzUMiIiIjKHZ1mvYQcGBsqgQYMM9l+9elW6du0qTz75pDRv3ly2bdsmZUHt2rV1D6xJhwyctl2zZk1ZsmSJBAcHi6+vr0RERMi+ffsMUsfQsmVL9T58f9o6uj179pTq1aurc3bp0kVOnDjhtM9IRETkzsp0xm7KlCny0ksvycaNGw32lytXTpYtW6aClxs3bkjr1q3l2WefFX9/f6uvlZ+fLw8ePBBn8PPzU8FWSSxfvlwWL14sq1evVsHbhg0bpH///vLf//5XGjRoIMeOHZO2bdvKgQMHpEmTJuLj46Ped+/ePYmKipLY2Fj1HeAc+C7Pnz8vFStWtNEnJCIiIinrgR2ySgkJCQX216lTRz0A2Spkm+7evVuiwA5BnbNKmRkZGSVqO7z//vsyc+ZMGTZsmNpeuHChHDx4UAXAf/3rX6VGjRpqf7Vq1dR3punevbvBedasWSNVqlSRQ4cOyXPPPVeiNhEREVEpKcV+88030q9fPwkKCjI5ahIBRXh4uJQvX17atWunska2dvz4ccnNzZWQkBApq9LT0+X69evSoUMHg/3YPnv2bJHvvXnzpkRHR6usHkqxlSpVUoHmlStX7NxqIiKissdlM3b3799Xw31RKv3jH/9Y4HV07p86daqsWrVKBXXIHPXq1Ut19kd/MEApNScnp8B79+/frwLG4iBLN2rUKFm7dq1NyqEIaJwB13YWlGHv3LmjSrlhYWGqf1779u0lOzvbaW0iIiJyVVpXJrcL7Pr06aMepqAjPzJBY8aMUdsI8D7//HPV9+vNN99U+0oy9cbDhw9lwIAB6lzPPPNMsRktfQhejCHrWNJyqLMgy4ZAODExUQ1+0GAb/er0b0RkN/XhmLi4ONWvThuYcvv2bYe2n4iIyNUhiXX06FH18xIxCB6m4oxSWYotCrI9KJFGRkbq9nl6eqrtI0eOlPj86OQ/evRo1T9s5MiRxR6PMi3KjNojJiZG3M2MGTNUvzpkSpEVRcCLwBkDUABZ0goVKqiRsii/pqWlqf0owW7evFmVbL/77jsZMWKEOo6IiIh+h1k4xo4dq4sj9OMKS7qDlcrADhkfZIZq1aplsB/bGMVqLgSCgwcPlr1796ppPLSgEFkmBDDo14dyLh6nT582eR5koRDIaI9Zs2aJu5k8ebIqfU+bNk2aNWumArhdu3apwE0bSbxixQo1ahbZveeff17tX79+vaSmpkqrVq1UkIzzaKVyIiIiKghxhH5cgTjDXB75SE+5OJQxd+zYoUqjgI78devWlcOHD6v+Wpo33nhDjbZEZsgRkBpFJI0vHeVKfVlZWXLp0iU1vxsGdxC5It6nRESur6h4wy0ydph+xMvLS5X89GFbf6oNR3nqqafUZMYYpUtERERkC4grEF8gzjBXqQzs0FEfkwbHx8fr9uXl5alt/Qyeo2B1hTNnzsiECRMcfm0iIiJyTxMmTFDxBeKMUj8qFlODXLhwQbeNchE661etWlVCQ0NVfy9MpdGmTRs1MhPTnWCKFG2ULBEREVFZ47KB3b///W/p1q2bbhuBHCCY++ijj2To0KFy69YttVA9Bkxoa5caD6hwBKRIURpGZM2sHREREdmqFIuH8VRipX7whKsyZ/AEVsbg9B7kqjIzM+Xy5cscPEFE5MLcfvBEaeDt7a1bI5bIVWn3p3a/EhFR6eaypdjSDqVZLHafkpKiW9YL07YQuQIk6hHU4f7EfYr7lYiISj8GdnbsY6dNvaIFd0SuBkGdM6YIIiKi4rGPnYvWvPEX8ujRI4e2jag4KL8yU0dE5F597JixcwD88OQPUMKCzlj/D0vF+Pr6Ors55OJ4v5AleL+Qhhk7B0XQRLxfyBK8X8gSvF/cWzpHxTpWaVhSzBFts8U1rD2HJe8z59jijinp666O94tlx/J+4f1iybG8X3i/WPJ3XadOHYuWFMPoOLJSWloasp3qT1fXuHHjUnENa89hyfvMOba4Y6x5nfeL7a/B+8U18H6x7FjeL7xfLL1fLPn7ZR+7EtCq2EiRujoM4LB3O21xDWvPYcn7zDm2uGOseV3b5v1iu2vwfnENvF8sO5b3C+8Xa+8Xc3rPsY9dCVy7dk1CQkKc3QwiIiIqA65evSrBwcFFHsPArgTy8vLk+vXrUrFiRU4+TERERHaBUO3evXsSFBQknp5FD49gYEdERETkJjgqloiIiMhNMLAjIiIichMM7IiIiIjcBAM7IiIiIjfBwI7IBezZs0caNWokDRo0kHXr1jm7OeTiBg4cKIGBgTJo0CBnN4VKwfQYXbt2VasjNW/eXLZt2+bsJpGdcVQskZPl5OSo/3QPHjyo1gJs3bq1HD58WKpVq+bsppGLSkhIUFMfbNy4Uf7+9787uznkwpKTk+XmzZsSEREhN27cUP+/JCUlib+/v7ObRnbCjB2Rkx07dkyaNGkidevWlYCAAOnTp4/s37/f2c0iF4YMDObPJCoO1hlFUAe1a9eW6tWry927d53dLLIjBnZEJfTNN99Iv3791MSRmKh6586dhS7kHB4eLuXLl5d27dqpYE6DSa4R1Gnw/Oeff3ZY+6l03S9Uttjyfjl+/LhanoorJrk3BnZEJXT//n1p0aKF+s+1MFu3bpWpU6fK3Llz5cSJE+rYXr16SUpKisPbSs7H+4Wccb8gSzdq1ChZs2aNg1pOToM+dkRkG/gntWPHDoN9bdu2zZ8wYYJuOzc3Nz8oKCg/JiZGbScmJuYPGDBA9/qUKVPyP/nkEwe2mkrT/aI5ePBg/gsvvOCwtlLpvV+ysrLyO3XqlL9p0yaHtpecgxk7IjvKzs5W5Y/IyEjdPqzzh+0jR46o7bZt28oPP/ygyq8ZGRnyxRdfqN+4qewx534hsuR+QTw4evRo6d69u4wcOdKJrSVHYWBHZEe3b99WfVpq1aplsB/bGKEG5cqVk8WLF0u3bt1UJ+dp06ZxRGwZZc79AvjBPXjwYNm7d68EBwcz6CujzLlfEhMTVbkWffPw/wsep0+fdlKLyRHKOeQqRFSk/v37qweROQ4cOODsJlAp0bFjR8nLy3N2M8iBmLEjsiNMLeDl5aXmkdKHbUw9QKSP9wtZgvcLFYaBHZEd+fj4qAlB4+Pjdfvw2zO227dv79S2kevh/UKW4P1ChWEplqiEMODhwoULuu1Lly7JyZMnpWrVqhIaGqqmIoiKipI2bdqogRLLli1TUxiMGTPGqe0m5+D9Qpbg/UIWc9JoXCK3gWkn8E/J+BEVFaU7JjY2Nj80NDTfx8dHTU9w9OhRp7aZnIf3C1mC9wtZimvFEhEREbkJ9rEjIiIichMM7IiIiIjcBAM7IiIiIjfBwI6IiIjITTCwIyIiInITDOyIiIiI3AQDOyIiIiI3wcCOiIiIyE0wsCMiIiJyEwzsiIiIiNwEAzsiIhc2cOBACQwMlEGDBjm7KURUCjCwIyJyYVOmTJFNmzY5uxlEVEowsCMiMlNsbKyEhYVJuXLlZPr06Qavde3aVTw8PNTj5MmTJs+B41577TWzr4njK1asWOhro0eP1l1z586dFnwSInJXDOyIiMzw/fffy9SpU+WDDz6Qq1evyrx58wocEx0dLcnJydK0aVOHtGn58uXqekREmnK6Z0REZNKePXukbdu28uyzz5o8xs/PT2rXrm3ReSMiIiQnJ6fA/v3790tQUFCR761cubJ6EBFpGNgRERXj8ccfl4sXL6rnKHuOHDnSrH5v9+/fl1dffVU+++wzVU41Lt9CUWVbIiJLsRRLRFSMw4cPS7169WTRokWq9BkXF2fW+2bMmCGHDh2Sf/zjHyoDl5CQICdOnLB7e4mo7GLGjoioGAEBAXL58mXp2LGj2aXWjIwMWb9+vXz88cfSo0cPtW/jxo0SHBxs0bUjIyNV/z5k//Debdu2Sfv27a36HETk/hjYEREV49SpU+rPZs2amf0elG6zs7OlXbt2un1Vq1aVRo0aWXTtAwcOWHQ8EZVtLMUSERUD/eDQz87f39/ZTSEiKhIDOyIiMwK7Fi1aWPSe+vXri7e3t3z33Xe6fampqZKUlGSHFhIR/YqlWCIiMwK7/v37W9wv7+WXX1YDKKpVqyY1a9aUt956Szw9+fs0EdkPAzsioiLk5eXJ6dOnZfbs2Ra/F6NoMYiiX79+arqTadOmSVpaml3aSUQEDOyIiIqADBtGpFoDWbvNmzerhwYZPCIie2FNgIjIRjC/HYI5ZPgcYfz48ep6REQaj/z8/HzdFhERWeXnn3+WzMxM9Tw0NFR8fHzsfs2UlBRJT09Xz+vUqcNRu0TEwI6IiIjIXbAUS0REROQmGNgRERERuQkGdkRERERugoEdERERkZtgYEdERETkJhjYEREREbkJBnZEREREboKBHREREZGbYGBHRERE5CYY2BERERG5CQZ2REREROIe/h9E2jHOOI++PQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Compute VPSD, add VPSD components, fit VPSD coefficients and plot everything\n", "example.star.compute_vpsd()\n", "example.star.add_vpsd_components()\n", "example.star.fit_vpsd_coefficients()\n", "fig = example.star.plot_vpsd_components()\n", "plt.show(fig)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "For further analysis, the best-fit VPSD coefficients can be accessed in the `ARVE` object. For example, the granulation timescale is obtained like this:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Granulation coefficient names:\n", "['amplitude', 'timescale', 'decay rate']\n", "Granulation timescale [hr]:\n", "0.962004134983532\n" ] } ], "source": [ "# Print granulation timescale\n", "print(\"Granulation coefficient names:\")\n", "print(example.star.vpsd_components[\"granulation\"][\"coef_name\"])\n", "print(\"Granulation timescale [hr]:\")\n", "print(example.star.vpsd_components[\"granulation\"][\"coef_val\"][1]*24)" ] } ], "metadata": { "kernelspec": { "display_name": "arve", "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.16" } }, "nbformat": 4, "nbformat_minor": 2 }