{ "cells": [ { "cell_type": "markdown", "id": "e800e54d", "metadata": {}, "source": [ "# Searching by distance from Sun\n", "\n", "Import the CDF library and open the orbit file, which is available for download from https://s2e2.cosmos.esa.int/confluence/display/SOSP/Trajectory+Overview+-+10+February+2020+Launch" ] }, { "cell_type": "markdown", "id": "93aed497-d081-40dd-a828-cf22519aa16a", "metadata": {}, "source": [ "## Setup" ] }, { "cell_type": "code", "execution_count": 6, "id": "d612df0b", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Today is 2023-10-12 13:22:24.671310\n" ] } ], "source": [ "import cdflib\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import numpy as np\n", "import datetime as dt # date\n", "\n", "AU = 149597870.700 #km\n", "\n", "today = dt.datetime.now()\n", "print(f'Today is {today}')\n", "\n", "# Distances requested in AU\n", "min_dist = 0.6\n", "max_dist = 0.7\n", "# Could be made the min distance after data has been read in below" ] }, { "cell_type": "markdown", "id": "69e958f8-e656-4356-bf2f-3106095182e1", "metadata": {}, "source": [ "## Read the file and print info" ] }, { "cell_type": "code", "execution_count": 7, "id": "21685582", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "All zVariables: ['EPOCH', 'HCI_POS', 'HCI_VEL', 'HEE_POS', 'HCENTRIC_DIST', 'HGRAPH_LAT', 'ECLIP_LAT', 'HCI_LON', 'SHURST_LON', 'LABL_HCI_POS', 'REP_HCI_POS', 'LABL_HCI_VEL', 'REP_HCI_VEL', 'LABL_HEE_POS', 'REP_HEE_POS'] \n", "\n", "zVar, FIELDNAM, CATDESC: \n", "\n", "EPOCH, Time, Interval Centred Timetags\n", "HCI_POS, HCI Position, Spacecraft Position in Heliocentric Inertial Coordinates\n", "HCI_VEL, HCI Velocity, Spacecraft Position in Heliocentric Inertial Coordinates\n", "HEE_POS, HEE Position, Spacecraft Position in Heliocentric Earth Ecliptic Coordinates\n", "HCENTRIC_DIST, Heliocentric Distance, Spacecraft Heliocentric Distance\n", "HGRAPH_LAT, Heliographic Latitude, Spacecraft Heliographic Latitude\n", "ECLIP_LAT, Ecliptic Latitude, Spacecraft Ecliptic Latitude\n", "HCI_LON, HCI Longitude, Spacecraft Heliocentric Inertial Longitude\n", "SHURST_LON, Stonyhurst Longitude, Spacecraft Stonyhurst Longitude\n", "LABL_HCI_POS, HCI POS Axis Label, Axis Label for Spacecraft HCI Position\n", "REP_HCI_POS, HCI POS Representation, Vector Representation for Spacecraft HCI Position\n", "LABL_HCI_VEL, HCI VEL Axis Label, Axis Label for Spacecraft HCI Velocity\n", "REP_HCI_VEL, HCI VEL Representation, Vector Representation for Spacecraft HCI Velocity\n", "LABL_HEE_POS, HEE POS Axis Label, Axis Label for Spacecraft HEE Position\n", "REP_HEE_POS, HEE POS Representation, Vector Representation for Spacecraft HEE Position\n", "\n", " Time resolution of orbit file is 1 hour: \n", "['2020-02-10T04:55:50.000000']\n", "['2020-02-10T05:55:50.000000']\n" ] } ], "source": [ "fileAndPath = ('/Users/hmiddleton/Library/CloudStorage/OneDrive-ESA/SolarOrbiter/Data/'\n", " 'solo_ANC_soc-orbit_20200210-20301120_L011_V1_00200_V01.cdf')\n", "cdf_file = cdflib.CDF(fileAndPath)\n", "\n", "# Read in all the metadata:\n", "info_dict = cdf_file.cdf_info()\n", "zvars = info_dict.zVariables\n", "print('All zVariables: ', zvars, '\\n')\n", "print('zVar, FIELDNAM, CATDESC: \\n')\n", "\n", "for z in zvars:\n", " y = cdf_file.varattsget(z)\n", " print(f\"{z}, {y['FIELDNAM']}, {y['CATDESC']}\")\n", "\n", "# EPOCH is in CDF_TIME_TT2000 (TT2000 for short) which is nanoseconds since J2000 with leap seconds.\n", "t = cdf_file.varget(\"EPOCH\")\n", "d = cdf_file.varget(\"HCENTRIC_DIST\")\n", " \n", "print('\\n Time resolution of orbit file is 1 hour: ')\n", "print(cdflib.cdfepoch.to_datetime(t[0]))\n", "print(cdflib.cdfepoch.to_datetime(t[1]))" ] }, { "cell_type": "markdown", "id": "6f76b20b", "metadata": {}, "source": [ "*Note that the CATDESC for HCI_VEL is position and not velocity - this will be corrected.*\n", "\n", "## Put into DataFrame and plot" ] }, { "cell_type": "code", "execution_count": 10, "id": "96f01e1e", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAB7BklEQVR4nO29ebhk11Uf+ls1V935dt/W1LIk27KMbIOx2gYTIAaeE5shjoGAnS8xMYOfAOdB3iMvJg+S8EzyMAl8CWAQ4sUYEoIfhMEmCBtiwAbjQW1blixbQ0uypNbUc9+h5qr9/thnn7PPOXtY+/Ste7tvnd/36VPfumefoW7VWnut31q/RUIIlChRokSJ+UVlv2+gRIkSJUrsL0pHUKJEiRJzjtIRlChRosSco3QEJUqUKDHnKB1BiRIlSsw5avt9A6E4fPiwuPHGG/f7NkqUKFHiisKnP/3pM0KIDdPvrjhHcOONN+L48eP7fRslSpQocUWBiB63/a5MDZUoUaLEnKN0BCVKlCgx5ygdQYkSJUrMOUpHUKJEiRJzjtIRlChRosSco3QEJUqUKDHnKB1BiRIlSsw5SkdwhWI8meK/fOJxPHuxX2j95568gL944NQu31WJEiWuRJSO4ArFb3/qCfzkH34eP/a7nwte2xtO8B2/8jd463vvxkPPbc3g7kqUKHElYWaOgIjeQ0SniOjzlt8TEf0CEZ0gonuJ6BWzupeDiA9Hu/lPPHoWveEkaO1nnzyP8VQOJCqjghIlSswyIngvgNc5fv96ADdH/70NwK/M8F4OHD735AWsduoYT0Xwrv6zT1wAACy1arj35MUZ3F2JEiWuJMzMEQghPgrgnOOQNwD4TSHxCQCrRHTNrO7nIOFid4Tz3RG++WXy7Xrw2TBH8MipbVy70sJX3XQIj5zensUtlihR4grCfnIE1wF4Uvv5ZPRaDkT0NiI6TkTHT58+vSc3dznjsbM7AICvfeFhVAg4eb4btP7k+R6OrnfwvPUOnjzXRTm3ukSJ+cZ+OgIyvGa0SEKIO4UQx4QQxzY2jCqqc4XHI0dw85FFXLXcwlMXwiqHTp7v4uhaG89bb2NnOMG5neEsbrNEiRJXCPbTEZwEcL3281EAT8/ygp989CzuefLCJZ/n3M4Qf3LfMxiOp5d+UwVwemsAADiy3MJ1q208dYEfEYwmUzy72cfR1TauW+sAkBFCiRIl5hf76Qg+AOAtUfXQVwO4KIR4ZlYXu/fkBXz3nZ/A33/3x4JTKVn8b7/9Wfzgb30Gv/Dhh3fp7sKw2RuBCFhq1nDtahtPB0QEz17sYyqA69baOLLUBACc2R7M6lZLlChxBWCW5aO/DeDjAG4hopNE9H1EdDsR3R4dcheARwGcAPBrAH5oVvcCyJ3wTYcXAAB/9Lni/ubpCz389YkzAIDf/fST+5Jfv9gbYblVR6VCOLzYxNkAQ342SgMdXmzicOQIVIRRokSJ+cQsq4beLIS4RghRF0IcFUL8ZyHEHUKIO6LfCyHEDwshXiCEeJkQYqZjx267YR1/8WOvwS1XLeFvHjlT+Dwff+QsAOAtr74Bz20O8PjZS4suiuBib4Tlthwud2ixgZ3hBP0Rr5fg3I40+usLDRxaaAAoI4ISJeYdc9dZ/PLrV3H/05uFd/KfffI8llo1fNcxSW/sBucQiou9EVbadQCIjflZJuF7dnsYrWuiVa9iuVUrI4ISJeYcc+cIXnzNEs7tDAsbv0dO7eCFRxbx4quXUK8SHtwHiQbdEaxHjuDcNs8RqAqhtQW5/vBSE2eYa0uUKHEwMXeO4JarlgAADz1XrJHqkdPbeMHGImrVCq5f7+Cx0zu7eXssbPbHWG5FEcFiRPju8BzbuZ0hGtUKFpsytXR4oVmmhkqUmHPMnSO4fl2VTIbn9jf7I5zaGuAFG4sAgOcfXsCXzu6DI4jIYiBJDYVEBOsLDRDJNo7ldh0Xe6PZ3GiJEiWuCMydI7hmpYVqhfBkAUfw+Bm5RlUf3XhoAY+d2cF0ureVQ93hBJ1mFQCwvhg5AiZHcG5niLXIeQDASruOzdIRlCgx15g7R1CrVnDtagtPngtvonrmolxz7WoLAHB0rY3BeIrz3b3LsQsh0B2OsdCQqZ3FRg1EwFafZ8zP7gzjKAKQjiA0IhhPpoXI9lHBdYNxmLqqQtGGP24FVokSBwVz5wgA4Ohqp1BE8NymbNy6elk6gqtX5P+fKTgc5qkLPdzxkUdwIcCRDMZTTAXQbsiIoFIhLDVr2OyPWes3eyOsdOrxzyvtOnaGE4wmPKN5ZnuAV//Mn+P//O/3su8ZAB49vY2X/9Sf4uf+9KGgdZ9+/Bxe9q//FL/+sceC1v3ZF57DS/71B/E/7g1rVv+vn3gcL/s3H7qkEuMSJa40zKUjuHa1jecKGO9nN/uoVSgmaK+KHIJyEKH4sd/5HH7mTx7Auz74IHtNN5o9sBA5AkDm+bnpHUk01+KfV6J+BO76P/vCczi9NcDvfvokdgY85wPIJr6d4QS/+tFHMAlIpf3eZ57CMJrGFoL3feoJjCYC7/vUk/6D9XV3y3W/e/xk0LoSJa5kzKUj2Fhq4vT2IDhN8czFPo4sNVGtSKL1mpU2AOkgQnF2e4BPPCab0z74+WfYxlEZ304zMeZLrTo7ItgejOKKIQBxdMBNDx3/0vn43/c9xZ9lcPxxqUg+moggov7ux+S6R0/vYDvA8Xw26u+458kL7L9zdzjGF5+R5cCfO3mBfa0SJa50zK0jGE1EcG78uc0+rorSQQBweLGBCqHQ3OBPPHoOQsgO5fPdEU6c4pWzqoigo0cErRo2GRzBaDJFfzTFUiudGgLAdiRfOruD61alA3w0oHT20dM7Mcn+MLN0dzoVePxcFzcckpVe3FLdC90hzu0Mcf16G9uDMbtn5PGzXUymAs/fWMDjZ7uFuYkSJa40zK0jAIBTgU1lz17s4xrNEdSqFRxZahVyBA8+t4UKAf/gNtmh/MVnNlnrukNpsBVZDPBTQ9uRsU9FBO2wiOBLZ3bwNS84hEatwi6dHU2meOZiD19/82EAwGNneOue3exjOJ7iG245AgB49AzPgTwanf8bo3VfYsqAKBXWr795A5OpKFRQUKLElYi5dARHCoqtndoc4MhSK/Xa4aVGoYasE6e28Lz1Dl58zRIa1UqAI5C71HYqIqhji7GjV8cstYo5gq3+CGd3hnj+xiKOrrXZKZ5nLkjF05dcu4J2vcpOpSkdp69+/iF5HqbDffKcXPf1L9qIzsNzPOp5XnnjOoDi3E+JElca5tIRbBRwBIPxBFuDMQ4vNlKvry80Cw12OXFqGy88soR6tYLnHeqwd9cJWaxzBLzU0NZgFB+vsNiUjmCb4UhUBHXNSgtHlprs908Z2KPrbVyzwo+gVLnui65axGKzxjbM6r5eet2K/JnpqE+e76Fdr+LWa5cBFEv5lShxJWKuHcGpLf4X/UJXGlG9GQuQnb3nAvsIplOBL53p4vkbMmcud9e8NIRKDbUzVUPbg7G3sS2JCBKOYCFqTONUAMWCdYsNHFlqsVNrz20lZbdXLbfYEYG63uGlJo4sN3Fqk3e901sDNKoVHFlqYrHJF9V7brOPq1dacXlwkSKAEiWuRMylI1hq1tCsVYLE1tSuf72TdgRrnQZb3iE+V3eI4WSKayO+4fq1ToAjiCKCZposFgLYHrqN+bYhNdSJIgtORY5KgR1ebOLIkjTMnIqcczvSiR5aaOLqgIjgbKSLtNSs4chSkx8RbA+wsdQEEckKMaYjUPIb7YZUZS0jghLzgrl0BESEtU4D5wNSOudj1c5MRBA4DwBIUg5XR+WnR9fauNgbsdI7cfloJjUEwMsTqNSQThZXK4R2vcqMCKRBPbTYwJHlJnqjCcuBnN8ZolohLLVqWF9osDuxz24PYl2kkAjkzPYwTuFtLIY7AkCK+YVGeiVKXKmYS0cAAKudOi4ElI8qo7C+kOUIwrR+gIT0VBVI161Jh/D0BX9UYCofVU6hx44I6qnXF5o17Az9juz09hBEMio6HDXVnWVEQ2d3hljryIlqa506usMJqzTz3M4QhyKDvr7QYHdgn9kaxPenekY4ON8dxhHfaqce1PFdosSVjPl2BAFf9Dgi6Fy6I3g2IkGVI1CVSJyda3c4QaNaQb2a/OmSPL/buG4aUkMAsNjkRwRrnQZq1UpQtdF5bae9Gr1/inNx4Yy2bqUtm+Y4jXentxNHIP/O/msJIVKCfDJiLMX4SswHZuoIiOh1RPQgEZ0goncYfr9GRH9ARPcS0aeI6KWzvB8da50GzjMMhILKc6920rvpohGBLlWh0hg8RzBOEcVAEhHseCKCrf4Y9SqhWUv/2ReaNTZZfGgh2TEDYEVV57rD2IGq/3PSQ+d20gYd8EthCCFwoZsYdCWq5+MytgdjjCYi9XxlRFBiXjDL4fVVAO8G8HoAtwJ4MxHdmjnsXwK4Rwjx5QDeAuA/zep+sljtNFg7RYXz3SGWWrXUThwoGBFspqUqVBUTpx+hO5ykdIaApJS064kIlLyEmkUQr2/WWLn+C71hbJBDIgI9974Wrefsts9uJw6E63j6oylGExHPdF5p1zGZCm/qS92PciDrgRuFEiWuZMwyIngVgBNCiEeFEEMA7wPwhswxtwL4MAAIIR4AcCMRXTXDe4qhdnxcHZrz3WGOHwCSKqIQR3Bmexgbf0CSt616pXhEoFJDnoigO5hgoVnLvb7QqHrXAjKiUPzCSls+90XGrtmcGnKvG02m6A4nscNZbfPWKTnu5VaYwzobTXhbj0Z4ri000BvxiwCOf+kcPvLQadaxCp978gL+4oFTQWvuPXkBf/lg2JoHnt3En33huaA1JeYLs3QE1wHQpR9PRq/p+ByAbwcAInoVgBsAHM2eiIjeRkTHiej46dNhXzYb1jp1jKeCLWR2bmeY4wcAfsWODt0wAggqc+wO88Y8jgg8u97ucIJ2vZp7faFZ80YTgHIEyU4b8BtYIQQuaDOW1axk325bpari6zEjAlV5tdzOOALP9dRzKAcXRyCMqKA/muA77/g4vuc9n2J3mQ/HU7zxlz+Gt773bnaZ6mQq8A9/7ZP4J79+d9w97YMQAm+68xP4gd88jkdOFxvPWuLgY5aOgAyvZbffPwNgjYjuAfBPAXwWQM6iCiHuFEIcE0Ic29jY2JWbCyEtAXtEUIvm/4YI2GWnhAFRmSMnNTSYpCqGAC0i8Di17ii/FpARCcchbvWTEZmNWgWdRtX7/vVHU0ymIo4k1M7e935l5TCUE/YZ9CwhznVY6vmXM+s4Jb26IqtSS/XhvqcuQPHen3j0LGvNQ89txff51yd48xIePrUd/40+xlxTYv4wS0dwEsD12s9HAaSmhAghNoUQbxVCvBySI9gAEDaBpCDWAh3BhW6yq82Cq/6poJcpKhxm1rvvDMepHgIA6ES7fF9E0B9OcmklgEcWCyFSEQHAm24W9y5E61r1CqoV8l4v6wjUe+9LDSkyWTmsZa4jUIJ80fVUrwXHQT58aiv+9+ef5klzP6QpsN7PXHOvJo390HNb9gP1e9Ov8xRPz6rE/GGWjuBuADcT0U1E1ADwJgAf0A8gotXodwDw/QA+KoTYk0+rCv25zU1bmYEuOkIGwPdHE3SHk1xEsMYkr3vD/K6+Vq2gUav4OYJR3okAiiOYOCUqeqMJxtrOHpDG2ZeqiXsXIsNKRKwIROX61fW4hlk5kNzOnhkRqOsoB8TRYDpxahvLrRquX2/jCaZi6cPPbaNdr+KFRxZjcT0fTp7voVohvPjqJbZsudKwevHVS3iCmU4qMX+YmSMQQowBvB3AhwB8EcDvCCHuJ6Lbiej26LAvA3A/ET0AWV30I7O6nyzWAsof5W54lGvEUgiZEHbe0pi2uiCNqo+8lhGBYVffqHrz/C6OAJCpIxtMyqWLjEgia2DVv71d0JnrNWoVNKoVbHt7JTIcQYeX4lHXi2dBKzE+RkTwxLkubjq8gOvXOniKqch68nwX16+3cXStjacv8pzHU+d7uHq5hRdsLOIpRvMhIOdAXL3cwi1XLxUaz1piPmDe4u4ShBB3Abgr89od2r8/DuDmWd6DDcqoc4a+7wwnmArEJYlZLLfqbEnmc5bGtNV2A8PxFL3RxLhrV+gOzb/vNGreiMCWGlLTzrqDccpg60h26MnvF5o1b6omm3IBFCcRtkNX5/CtSyIC+fddYGopbUfPXolKelWTHiciOL01wPXrHay26+zKoVNbUtL8utU2PhdNU/PhqQs9XLvawpHlJj7yEF9I75rVFq5f6+B/3PsMplMRP2OJEgpz21kcUu2TTVNksdLmzQMAknr1XETAqFIRQkSOwLSrZ0QEFrJYcQz9kX2A/WbGwAI8knnLatDDUkOAfEZv93RvhFqF0KrLj7bSUvLxJ9v9cUrIbymKCLYYEcHpLSlyd81KC6e3B7zu560Bjiw1ce1qG+e7I1aZqnQEbVy13ML2YMxqAjwTdVlvLDUxmQp2KrTEfGFuHUG7XkW1QqyIYLNnlmZQWG7zq4YSzaK0U1ljOILhRFbgmHoBOBFB1xIRqNd6gakhjmE2KZ5KB+JetzUwXK/BcSBjLLfrqaa5hWaVHRHoa/T7t2E8meJcd4iNxSbWFxoQwk9oCyGk81huBjUkKudx1TJ/wp4U4GvGmk1nC8zOKHHwMbeOgEiqYe5GRLDckvMAxhP7jlrBplmk6tcv9OxfVLXjt0YEjl3vZCowHE+NHIF6ze0ITDv0ghxBq4ZtRs4+K4ex2Kx5DfNmf5Rz2Nz7XNSerVatoF2velNR53aGEELOTFiP5DB8Rv1ib4ThZIojSy22I+gNJxiMp1jtNGJtKp8s92QqIpmOBg4t8LvXS8wf5tYRAAhwBOlKlCxUZQqHWFRf+GwpKic1pMhcY3qn4TZ2PcfaZpRG6TkciYqKdJ5kIYpCXAR37Ai0926JWTW01Erv7Bdb/qhnZzBOTW8D1Hvjk98Yx5VN+vV896l25RuLzVinyLfrjtcs8deolM5ap8Eufb7QHWIq5PAkpWcVMoOjxPxgvh1Bs85LDfkiggDdnYs9uWOtZTSLOF/urmEWgUKnUXXu6JPJZvm17Zgj8K9faKbJ4qlwRxLbgzEa1QqatcQBcXb225meBXU9n2E2cSgcddXtfp4o56Sw1A57Y6mRpF88xlZFheudhhYRuHfqiSOox5uGi47oUd6bmijXjAUOz5YRQQkD5toRLLZqMQnqwqYnIlCvq12zC7IfIe9QElE1+5d7xzCLQKFVqzoNudrtG1NDDI4gnoNQ1w26X/46S8ICamc/cZKq24ad/WLD70BMPIict8BJDRkcgWejkEhT1NlGXX2e0mt8u/tEFC9prvPoJ2mDhFbbdRAhaBhTifnBXDuC5UCOYNnSWRyixLllyGEDQKteRbNWcUcEQ3tE0KpXnFU/poE2CjFH4EgNqTkItdQchEj+2rHbthlY9TvX9YwOxLOzNzXcLXjSZoD8u5giAm6/w3Krrkls85zHcruG5VYd1QoFRAQNdBpV1Kvk7YFR11ltN1CpyEY+zsanxPxhrh3BUouZGuqZdfwVlKHzEYuAmcxU8Gngu8jiVt2dGlK/M1YNMcjinkH1dIFh0Lf6+Z29WudzPNk0lpqk5uqANnVP+6qbhJDig/lUlL/sVG9gq1craNUr3s+ULoNRiUZ4+qJJ5VzWOpI3WWn7O9GzlV7LrXqQFEqJ+cGcOwJ+RJAlLlPnUTXnTOLZxjUsRdVHNiiyOLtTBqQjGI6nViPZM6R24rUNDkeQ70pWBt612+4O87l3juPpjyZo19Mfz0WG3LZRlM9TWjsYTzEVeSfpc66A3CRIDkTeq9xc+Cub5LGJnIWP+7gQpXRWtfkMPo5gMyPJLTvgy4igRB5z7wi2B+6qFyAtv2zCAlP9U53LxjX4HJOLLG5FxnUwNqeH1M7WFRE4HYGhGW2BYZh7o3zOXv3cdRl0Qwd1LIXhiSTyZLFMDdn+zjYn2WlUnVELkKT61CaBs7nY7EnnqNJsnKKFC70ROo0qGpHDWW37R3Cq+1iMI4IwccQS84M5dwRyepUv/Nfll01QXzTOAHiXZpEvJ+3K86tOWtsONuEX8mvr1QpqFXLv0A0kbJLrd5PUrYyB5VUp2dfZjPN0KiLHkykfbVYxFfbOaVvarNOoOZ0VIIlfnTvipF82+6PUZmCJUbSQraJaYTqChUY1noQXoolVYr4w546Ap0Oz6YkImjVJ3vl2giYpZx3LHs7CRRb7jGs/Nnbma7frVfSGbrI5PwdBnqvvcIAyxZPfaatzutZlr+dLKfXHZkepHJYtclHnyzqeVr3qJOABme/PGnV/RDBKOQ9OOml7OM6V7vorodKbjmXGdUrMJ+bcEfCE52yVPjo43av90TQn5azDp92zM5ygVqE4PaCj5XEEpvLP1HpfH4Jhp80imQ2OoOXZ2UtNpXF+nafM1RYxKcdp02KyldZ2GlUMJ1Nnx/hm32RsORGB7ghq3jU7OQkMTiVUetMhSekyIiiRx5w7gqj+n1EiaDPeChxJ5ixJaLof147NVBqpoFJDth2siyMApBF09yGMc07El45S92zjCGzrhhMzeRtHPRYHYjPo8Xsz9kVLYREIoLSNQiOC/BpfVJrtmF5k6CdlHcFyu46twZglildivjDXjkCF9N6w3CHPrLDYrHmVKn39CIutGrqORqudgXmwDAA0PUZrEL1uK4GVqaEwErZV8/cf9EdTQ+7dvc5m0H2GOYkILJGL7Xojy/U89wmo1FB2d8/RQ8qvcUt1THKpof7IHa1k+ahlZip0qz/CD/3Wp/E3j/BGW/7FA6fw9v/2GVaxxHQq8K/f/3n8t08+wTp3ib3BXDuCuArF8QEWQjh34vq5/BGBW8VUfWlt3bPd0SSeT5yFMmIDmyMYT9GsVawlsL7UkGlnX6nI3gpbJDGeTDGc5IXu1M82jsCmi+SLJGyEuM+BKENvJacd70u2L2SpVUdvNMHIYaCzGwtVtOC6jkwNpWU6AHeBwlY/3cyXVF25P6cf+NzTuOu+Z/ETf/h553EK/9cf3If/ce8z+ND9z3qP/cwT5/EbH38c//IP7mOJNJbYG8y3I2AMLRlOZF7fJP2sg5Ma8onXKdGzLUtjWtcgu6AQcwSW9IdyBDa063aDDiAamGOWp7Ct8+60PTt7a6rGF0lk1jU9/Im9asjtsCZTgf5omvpsKKfg7q1Iv5ec2Rg7gzxZ7LvOZqZUWV3TJ8B330k5Q/mJs12nQwPke/70RamCyhmwc492zJeYIzpLzB7z7Qg4demObl4dIakhe0OZ2yDY5gkAfo5gMM6XY+pwcQRCCCPpq9ZZd9qqGidzz41qBRUqkBryNL7ZtJh8FVU9yzo/OZ2PQHzOYzyZYjieptJXHMmN/LwEjrxHOjWkNhG+iODB57bkvU4FnvTMOVYzkfV1LjxyOjn+YcbxJfYGM3UERPQ6InqQiE4Q0TsMv18hoj8ios8R0f1E9NZZ3k8W6kvrrtRJz7K1YYGhcGka7qJj0ZPD7Q4nWHCQvYDduA5G01hu2rjekRrqj6YQwlx6Krtvzc6nH5WjZg06EaHTqHkdiDXXH5ga8jlJW+TSYUcgyX22PcbWJCXu+9uNJlMMxunIQ6WJbJ+VOFpp5CMCX9/ME2e7uOWqJQBpQ2/Cl87I3z9/YwFPX3DPRwCAx85s4+YjiwCAk+d5c5dLzB4zcwREVAXwbsih9LcCeDMR3Zo57IcBfEEI8RUAXgPg54iogT1Cs1ZBtULeDlcA1ty8wmKzzhqjCMDanOYrZ5WD64ulhvrjSUoKOre+ZncErma0loNkthlYtc7KEcQGNv3xTMpOLQbdYJjlzx6OwNJH4COLTSW5HR//Yfg8xeWtljVqg5F2BPXod3yepcPgCEbRxLWXHV0BAJzadIvhPXVBGvNX3rCOZzf73i79Zy/28eJrltGqV7yDdUrsHWYZEbwKwAkhxKNCiCGA9wF4Q+YYAWCJJIO5COAcgD3reCEiLDTcgmTxl9ATEahyPpcg2lZ/jGqFrGkmFfrbUkPO8tGaSn9YUkMjN0fQatgbytzyFBUMLM4nyb3nryvlG+yRjzx3+j2vRj0U3qqhwI7k/nAConxFVbyDDuhb8O26TWt8khvJlLdkzYInIkjmT2hrGBzB2W05ce3Wa5YByPGYLpzZHqJeJdx81SKG4ylDGnuIw4sNXL3cwrOlI7hsMEtHcB2AJ7WfT0av6fglAF8G4GkA9wH4ESFEzhoR0duI6DgRHT99+vSu3qSv2scl65A9D2A3GoD80nYaVWvljq+8bydab0IzTn8UI4ubtQqGHoNuI4t9qRMTN+HiFmx1/fE6T4ewSTwOsEdLiv/I/l3iddaIJ29s1a7b9p6oz5ru5HzltMpwpyMCN0dg4lk45bDK8B9da2OlXffORT63M8D6QgPXrLQBAM9t2Y17fzTB1mCMw4tNHFlulRHBZYRZOgKTtctul/8ugHsAXAvg5QB+iYiWc4uEuFMIcUwIcWxjY2NXb3KhWXPmTE1huQmLrGoRe9WPfg5rRDCaxIYmC1ka6nIEbrK4UatYBetsZCrANOiWYTghu2bO9QYj885e/Ww36GYivOPL9w/zBjquzPE4q4UmP4rYNnwGfRLgJucdq8U6UkOnIkO+sdTEkaVm/LMN53aGWF9oYm2hHv/sOhaQozM3FpveOc0l9g6zdAQnAVyv/XwUcuev460Afl9InADwGIAXz/CeclhouDs0XbthHb60DhCVDTq4hnZdCoSZOILheIrRRFjJYiJyTinzRwRVDMZTY46369jZu6SafTMQ7HIYY+v1JKltIactvRJEJAf3WB3d1Bq1yOfwpMz0XbeHI3CmhmzVUHFqSHMEHl4h/psZU1D+iGBjqYn1hQbOeyannYlSPWooz0VHauisNjpzpVNnDXIqsTeYpSO4G8DNRHRTRAC/CcAHMsc8AeCbAICIrgJwC4BHZ3hPOfgUJnfi8lFP1RBLm9/dmObiLGxEqA7XlLL+yE0WKycxNNSNKw7A5gh8kg+dev6eXRFB3+F8XeT0wPGMrs7pvkEqW15LRlnWVJSDI7Cny/LCgYkWkjvfr18njnJ8JbH19BpfcYSaerYejcT0KanKiKARj1p1TWc7kxmdeaE78pLLJfYGM3MEQogxgLcD+BCALwL4HSHE/UR0OxHdHh32TgBfQ0T3AfgwgH8hhOD1te8SFjwDyuMvoa9qiJka8jqUptkxJWWsjsofV9pk7C4fVYbFlB5Sr5kiina9at9px30EJpE8RypqJMX16lXT9eyNb4PxNC4VNd6nIwIxpYaIZOe0b8aD0ah78v2m8lG7Y8yX4fq6unuGe5Nlu+7iiM3eCLUKoV2vYqXt37UrR6AiAtfMbZUKWu9IxzGeCpZ0e4nZw22VLhFCiLsA3JV57Q7t308D+DuzvAcf5DhCRkTgyK8DWmrIExEcWnBXx3YaVeOXw1QFkoXL2PmqhlT37dDpCMLIYhdH0KpVrdVG/ZE5VRNfz8FJ2CICXwrLdI+AfGZfCkv/m6gowtdHoK+pRkY9NMXmei+S6+QroVyfdyVUR0RY9jiC0WSK7cEYa50Ga+a2rrW12o4cR3fo1fHqDSeoVMyfPxMudkdY6bhFIhXGkyn646n3Hg465rqzGIjGGHoigmYtPbTdBI6Oi0s0Tj+PKUUQE5OO9Y2aPTXkI4udEYFDsE4ZWFOI76oaatYrGDjutYg4nosHaTmcZG80zXU/x/fJigiStUSEjrNHwlyO7DLQ6r5bGUPoei96sZPKTHlruIsjdEG8lXYd3aFdN2k70yDpm7m9rXEdylD7yk3Hkym+5Rf/Cm/4pY85S7MV/vLBU3j5O/+ULWr3I++7B1/3rj93chvzgLl3BIuejuCdzEAQGzhdmxzxOltEYEopZNGsV405foBXPgqYReviiMCQdonF7gzGsjuaoF41p3gUOW3CsKBBd6W/nPzJMD8fOb5PRwqrOzQ/X9thbNXfMT/7wL7GGhG4ohzbjIWm3UkBaVntlUgl1zbDYDtDYq+2G06OYKsvN1WNWgWr0bl9qafPP72JR0/v4IFnt/DwqW3nsQDwR597BkIAv/+Zk95jL3ZH+OP7nsH57ggfeXh3y9KvNMy9I1BSBzbpZ9MwdBPiigxXc9pw4nUqCxbyWtWs28pHAaBZNfcCCCHYZLGLIzBX1thJS9c15U7bZdDN61xlroPxJLdrju/TkUbpjezRkiuF1bNwC65mOZWGqlTSlU3uFJuZo3E5RZOUBSCJe9fGZ6s/wlIziQgAu7FWjkBFBIstz7kHyXyEZaYjuOeJ8/G/7z15wXksAHz2SXn8F57Z9BLRD53aiv/9xWc2vec+yJh7R7DoSenseGr/FXzSAoBZyjl3nmbN6EzipiLH+katYszxj6cCU2GfRQAk+VezI3CnhgDzc48mU+M0NXUuW7mqi89wRRJ9h56Sq7R2NJmiYUn9NR2RhKwCy382bFGdXGNuCuw4qqgGo0nEPRich40vsXRLNx1ltEB6aI76v88RKLkLX3OmPuCJI7QHAE+e78UihT4BvOlU4OS5Xvxent52N8M9+OxWfC/zLoDndQREdJSIfoyI3k9EdxPRR4nol4noW4joinckqhrIVY/tqxgCgFpVhrw2hzIcS21+lyEHpKE3NfzYJJZ12PLZrtSOvhYwp4aUITQZS3VOkwMajYXDwFYhhL1c1e4I7B3QA4eekmvegtth2SOCrkWau+PY3XcH5s2AM99viVjajpnK3aG5W7pdr1pnVgDpYTbLsfaVRfoiel1VzC15Rq1u90exA+DMAgGk8b/xcAfXrrbxhMcRnNkeYDiZ4mtecAgA8LhH5vrJ8100qhW88sY1lmDeQYbTkBPRrwN4D4AhgHcBeDOAHwLwPwG8DsBfE9HXz/omZwlf/b9U/ORVFCw4dnWcPgAgyhWbIgKGCqotIlBffCdZXHelhibWoTauSGI0maJeM8tp+FJRVsPsyNm7IolWzd7vMBxPjTyGuk8bqd2zbBKcKRsLT9RpVNEd2clim3Cf1eFYnJTr3oCII4gcgK8UdivDEUgFXve51bFJB7a7fPTk+R6OrnVwdK3tVStVAni33bAOQArcuXB6a4CNpSauXmnPvdyFz8L9nBDCNKbo8wB+P2oUe97u39beIdF1t3d1+ko+FVyEn/qSeyOCpowIhBApw8tRQW3WKpYdtr0PIFnrSA05DKza8Rsd0MRtYNW50crfr7Wcs5qklLKOSfYR2LWYbET60BERtOpVx3yIsbFZrlWvWhuxBmO7nEVvaN7B9izltC1HT0V/aI4iWnV7mepkKlJ5/IU4WnZHBMnx7tTQ9mCM5613AMi/f61C3tTQMxd7+MrnreJib4T7n3bn8dWu/isi5VSfYJ5yBNestHB2Z+iMKA86nBGBxQnovx9G0hBXLGJFRoeeDIcsBtwlgHHVj4cs7jRqmIq8Qe4OpXKpLdUCRESqYfeqjIXrQ64Moa2PwEbeJpGEgSMYO3LvseMxVSk5UkOq38Fg1CU5bXdYtkhiNBGoV+2Ri6vs1JRuazrKePujqfHv4Ip0+o7UkEuJ1abVZLu3LPmrIgLbrn07mqSnNlOLzRq2h/bZy/roTCLyOo7JVOBCb4RDi00cXmzijMewK12kF129hFqF4k5mG+KIYFnuRHyS2wcZTqsUdfzqf1UB4AyAvwDwH4QQV3w81Wn6UkNjr/GOz+UUUouqfjyNaWoXtjMYp778O1H1kk25FIhSQ4UjArtBd+fs7Y1ozty7KxXlIH31lFLWoLpKZG2GdjIVmEwFGtXwPoLBaILWUjP3eqtu5xUGY3PlmGxcczkCQ+mup6HMNkjI1yCnp3oA+3djuz8GUfKZXmjWIESUTjU841Z/lJrFsejhFC72RhACWOvU0agStgZjq1MEEnmLtU4DhxYb3ojg1NYAt92whvUo4j+3M8T1UcQyb/BZuG81vLYO4HsA/CKAH9j1O9pjKI13265nZ2CfCpZF2xERcAfc6HnZQ9rrPQZX0ayZiUAeWeyqGnKkhhy5/qEzNRRdz2D8TEY+WeeKXOxGQpbW5lNKqlnKzmXYjbqt38Gt+TTF+oIhInCU0/ZH5rJYp36SpTeiWZdVV9OpyJWwZudOtGpV2SVtMdZbgzEWG7X4PPr4zKwjEEIYxm1WneXWsSTFQiNOp53dGeK61bbx+IvdIZZbNVQrhI2lpjMiGE2mOLczxJGlVqyTdGGORfCclkUI8bjh5ccBfJaIPjubW9pbuAS/plMRDW3nksU167CNRDjM10dgTlXtWMoOdVgjAktnqg5nQ5klnZFaZ0wNOVIujpQSJwLJOh7ZK+GKCJKUkv4s6v2ypbBcmkg2TsLmkOWaidEht+pVKyndH02xsZSXTGg5jHpvNMHhxTy31dIisWz1UrYJrVKRXdLW1JCW6gFk1RAgU0xHMsf2RhNMRV5K2yWJrbqU1zoNdBryvTm7PbA6gvPdEdai3f1ax93cpprk1hbqWO0kchfziksp/7ziS0cBt667STveBVdTUCIR4e8jAPLkNaeMtVGtYDQRuVZ8VkTgqxqyrPVFBA2vAykYgWQMrTLo1kY0C6k9in62Ri7OTmazw3LV6tueTRH9JhmFniM1pM6ZxdBSeaWMvOmZTB3MHYsIIiC/M6YZCS71XH0z40sNqYhgrZOom7oa0M53h3HH8rJHOXUzIrqXW/UkImDITOwMxvj5P30wrlDy4VOPncN7/voxlsrqdCrwy395AvedvMg6927CxxG8wvDyGoB/BOCjM7mjPYZLGmKHuYtXWGjYdzhdJlm8EN9P+jwc5dK4pn8yRauSfOFYZHHVk7N31PXb1slGLXvKxXk9qwic+Xo+HkR3dEupe5RfUE7jW5afsRK/taoxDaXW2OS81f1ld+oushiIupUza0YWJ+ya1mbSTnKNcpWl1dqxjvGZpvnVC42as2zzvIoIFurYGci/z2bPzSkoFdSVdt0qjSHPowTwarHzOM+ICH7j41/CL/z5CTx8ahu/8o9ucx4rhMD3vfdubA3GeOl1K3jVTevO4z/8wCn87AcfxG8uP45P/Mtv8t7LbsJbPpr5WQA4C+AvAdw5ixvaa9SjRjCTAe8OwiMCq8YMkyyOKzUMEcG6p4xVN+a64eCQxbWq1Kq3pWrWLNd2kcW++nwgv7MXQniqhiyOYFQsIhj6IoJaBUJIh9GoZctVbWke898hXmOJCNTviziCLAbjqTEtp+7NFLkmqaG0rLYtIsg2uiXDcvLHx8J5KcfhFnzUyd9qRf7btcs/3x3ipsMLAOROf7OXL8NWUOdZbtVRq1aw1KqxIoKPP3IWAPDJx85Zz63w2JmduNfiYyfOeB3BJx6V5352s48z2wMcXswXIswKPo7gG2y/iwbJPLfrd7QPWGiYSavgiCAS9DJ9QLhksa12e2cwxvVr7ooGm5S0S0Y6td7WkFaQLB45yGJlkLIpFJ8chs3xJFFPmANRKSUbl5Hs1CepqGEyFRhNhJksriXpl7wjMEcR8U7dWP5r7qtoxkY9b3glF+JKDeWvY5rBvOgw1v1RenOi1hnPPczPVFA9Mzac7w5Rr8oZCiqx4trlX+imI4LhZCrfO0M6VkUWSvNotVNnRQSql+HczhDndoY45DDWX9D0ix7WdI1suO+pJCX08HPbe+oIgvL8RLRCRN9LRP8TwGdmdE97jo4lpcORfs6eZzIVRsK2Oxyj5ukDUOcA8lVMHOXSZjXZVepIJoy5r+2SqChEFk+EU7oBMOT6o+u7UjWm6/miHlUemuMIJp6UUjwJzBxJ2CS29XtSmE4FhpbhOa730cYRuJyHzQk3A1NDHce8jmw1kytCMaWGXFVPgCSjl1r1eHJfhexyF+PJFFv9cSyUp3SSbBGEHhEAwKEF/wzlncEY53aG+KpoZ3/Co4aqOqFfddM6Hj294zwWAJ4638Orn38oOrffcewmOFpDbSL6biJ6P2RH8c8D+Gmk5xFf0bCVsan6aZ9QnEL8RTAOlpEhvyuUVPcC5KuYdhiOwNYUFqdNvBGBuXJl4GjUcpdzMlJDgdFLQhZn17l5EJuhjctHAxvfnEJ8WkSgIya0AyKC0WSKyVQYI4L4723YeAwtzXzxvblSQ5mB97aqoSw3oRygm4hO7qndSKqeTJBlqPL8alCOzbArXkLt8H3KqTpHAMgZzb6+A2XYv/aFh1M/2/DkuS5WO3W89NoVr07SaDLFMxd7uO2GNVy70sJfn9jTQY1eraHfAvAQ5BSxXwJwI4DzQoi/FELYJQyvMHgjAiZHsODoSegOeJpFqnY7ew5OY5vPuNp22fH6urmW3aXzT0TWrt0iZLHLwLrXqR26O4UVzBFYdvd9h3P18hgejiB9HbtOVNPiFAGl82SvNDJFBKYZBq161dpH0BumORBnRZJhSFHbEZ0AwHbmO7PUqllTQzuZqjy107cdv9lPRnICwOFFjiOQxvy2G9YAyIY0F54838P1ax1ctdxEdzhxVkg9e7GPqQCuX2/j61+0gb955CxrEM9uwRcRvBTAeciZww8IISZIdxo7QUSvI6IHiegEEb3D8Pt/TkT3RP99nogmRORmVGaABctwGvUaNzWkOjlNOVtOHwCQ1G7rX77heIrRRHhLT33GzusICqSG1LrincWZnXZM+oalhnyVUTYnOeRGBNYIhB8R9Md2o65XDenoOR1B0huhYxrxF8aIIN61G9KXhkFC7Ya9FDYrhtdyRMTqvUjPd7ZX7AHy+7ek9Skst+px2afpWCApYfVHBMlITkBGBOe6Q4wtelQA8HRUMvrCI4tYarorngDg2Ys9XLPSwlWxhIX9eFWOeu1qG19+dBVb/TG7RHU34LQMQoivAPBdAJYB/E8i+isAS0R0te/ERFQF8G4ArwdwK4A3E9GtmfP/eyHEy4UQLwfw4wA+IoQ4V+hJLgG28X29+MPLjAgcH2ybUqUJnWY6HOcql9pSBaPJFBWSs3F9613qo+51Jo6AKTqXupY7NWSTvfb1SticpLd8NCa1LZxEQO7eGRFY0ipqjTsiyKS7pnbHb3NSgPycZa/jmuOQlb6oV6WQHJcjcDkOIN+nsNyqx3OPTccCyMlc29Jam/1RnEYCgEMLDQjh7i4+t6Oa0Bo4styMtY1cxx9abOJIJEPynEPLSPETG0tN3HL1IgDgoT2ckeDlCIQQDwgh/pUQ4hYA/wzAbwL4FBH9jWfpqwCcEEI8KoQYAngfgDc4jn8zgN9m3veuwtbhGA+DYWoNqbDbVGWxw+gDiO8nI1XBVS617V5dBjm7PmvQk+oYd0SQvabaldquS0RGx1M4NaQMpjUiMOf642gp1GGxUkMWXsFFFgdEHtYox/FM6jNqNNYGHkppE2UbokaTKcYG7sImamdKO7nuBchLVSy2anYl2MwoVxvXprDZG6WiDV1vyIbzkYRFvVrBkaWW07BPpwLnu0OsL9RxREUEDscRy2l0Gnjhhux04RDMu4WgqiEhxHEhxP8B4AbIHbwL1wF4Uvv5ZPRaDkTUgZxv8HuW37+NiI4T0fHTp3d/tmjHUj7aHY5RMUx5sp8nSg0ZdOXDVEzTJXtc5dIkIsiTlL60EGA26C7jFa8zzEp27UpT17NW/3jI4gAjm15nrhpyaQ25rmes5qmZIwIXr2Br9HLxO7ZSWhXlGPsILPcGqBkG6c9Yq17BVCTnVLDNUbYNAIpTXBpZrL4ProhgUbsf17yP7UxqyKecupPhH5TU/NltuyM4uzOMHcb6YsNZbrrVH2MyFVhfaMZSH04nE/1utdPAcruGxWYNT1+8TFJDRPQTppy9kPgIEX0jEZmE6QDA9M2y8QvfBuBjtrSQEOJOIcQxIcSxjY0N1y0XgisiWGjUvJU+8XkcEUHQgJsMZ8FVLrVV8LhGMWbXZw26K52h0KjmHUiccnFc1zR2ks0RWHbo3p29xRG4tIZM61wOyxcROMtHA95Hm3NLOCH7vbFTQxYH1TeQv4CKCMwcQYXSz+EacwpEjkDbtbcdzW07mdRQzD/YVIVHac5ujRMR7Azj41bbdWcD2rmuEsyrRyWwbgmLc90hlpo1NGpyANS1q62Yk9gL+CzTfQD+iIj6kH0DpyHHiNwM4OWQk8r+nWXtSaRLTI8CeNpy7JuwT2khQH5o+iNZpqfn0WWlDm8XD2ihrrFqaMwuQ11s1nBW+0Bym9FchoGTGmoYtPc5zWimaqOkGsfuRE1a/77UUK1CqJCd9PU5kPDOYnMqytXAtpsRQeyoAiqNfGuIbCWe+YKGplYJpEtI9y3cRatesZdPZ0Znql276V6mU5GTs3ZFBMpBqOMTxQDb8ZNUhK0ignOOXf65nSGuXZVpntVOHRe6Q6PgnzxWpo1kVzRhuVV36iRd6I6wupC8v9eutvd0fKaPLH6/EOJvAbgdwP0AqgA2AfxXAK8SQvwzIYQtV3M3gJuJ6KZoktmbAHwgexARrQD42wDeX/wxLg221vidgF286zyADLu5ctYLzVo8/Uk/n+9ebHpBo4mwpj5S62vVXArAZ5jVdU0ENQBjCWNyvQrGueu5HQ8RSS2fiWUXHGjQhw6jKddZSFwnWWxbw8j3WxrXzM1hbuLc5ISJyEoAm9KXLct9mchf9bMtNZTdCLka0FSEvqhtfqQA3sRYVrltkIPJcm2p+xlOUhF2HBE4UkPnu8O4c3mt08BUANuW8yti+dCCJIqV47Dh3M4Q652kS/ualfaeVg2xrJwQ4mEAD4ecWAgxJqK3A/gQpAN5jxDifiK6Pfr9HdGhbwTwp0KIvWNGMtAH2C9pu56QXTygkcWWPgLugJulVlqVcSdDhNngklHgpIbqVSpUetqs51NDPsNsu54ylv7rhTWGWauGCvYRuBrY/L0O/IjAJYFhd/zudJ5tXkJvOMnJGrS0iCB1rKFBTB1vTA0Z0k5OwcfoM7/YrOeO74/zXMbOQPJ5umPKcm3Z4/XvU71awWKzhgs9s7EWQkhjvZBIWADAxW562I5CHBEsJA1urookSSwnjuCqZdnpPJ5MUWN8dy8V/O1uAQgh7gJwV+a1OzI/vxfAe2d5Hz7YBtjvDMdBEUGzJoXbsqHxaDLFcDL15vj1+9EdQSzh62sos8koMFNDzVolNiIKvh2zXJef6+tKTyg0DNdjcxK2XL9XEyn7t3GXj8b1/dayTn4pqCudFBv1LEfgcMQ2oUBfuqtl2bWbRPJspbC2Rrd2oxoTnzp6o/ys5pYjIkjGYKZ3+ICa1pdxBNF3VU89LTjkMXqGCW7Ljqqk3miCwXgazy5IZhiMcL2h8+mCJpgHRI7AxRHsDPHCjcX4Z6VhdK4rh+fMGrN3NVcAbDrqIbX/gAy75SCPrIQ0z5ArLLZkCDyJQuCdIa981NpQxqwaqhtTPH7S10QW+xq1bNdjjdU0DItRz1yz9EpUK4RaxRDxeETnbLtuJ1lsqe93pZMqFXM5rW9wTqOab+bzOUUboTs0KJa2LH0Utka3Vi0gNRRzannjux1HBGmyGDCnXncGeT6vY5HHGE1kg2Y2wl5y9CkoB7GiidQBsEYQ21GEoq6x2mk4OYKLvXRfw2FGFdNuonQEcE0FC+MIAPNwmi7TkCssxs0w42i9uVQvC/uul0kWVyvxDlRfCwA1F+lrIItHY78DqRuMGCcVZTJ+w6iT1lXhZeqAjnfPFbuh1Y9TcOX74x6JXLQTGU9H93PWQHu1kBxO0T1sJ58aMimW2lJDqmooxxHYykcNqaFEm8tQypopBwXcDZs7g/ycZKkhZpeOyTsCe0SgXlcRSjLDwO44FptJhCKrjOxpp+1MF7WKCC4rR0BELyKiDxPR56Ofv5yIfmK2t7Z3UDv13DCYAU8WQsdCM9+lzDXkCrEjiD7E24Mx6lW/cqlt1zsam+UGsqgbykd9OXTAYmAZZLGpXJWTijJHLv6ox7Tjlk6SjJUfgNyp16uUu17fk8KSTjVNavYdEYE8V76hT53D9j6a3vuB5z2UHAGvE9xWAWUji1v1qtGwZ+UoAPl5bdQqccOk7/y27ykgN02LGUdgiwhM4npApGXk6VxWxlpxidsOx6HzjSttWTVkmlQmpevT0Y/iC87uuPWMdgvciODXIBvIRgAghLgXsgroQMBW/78zzO8yfGjX83nJeMANu48gGp8Zfch0OV4fbEbZZZAV6tGoS/3DyknxGFMajPJReT1LtZFlhw5YuAxDWiO/rmpMo/iipbohUhqMJ9LxOsjpbGOfTwXWpPU08KWGTM7NQ9S7UkPZNbYKKOUYTJVA5j4C83yETqNqVEI1pdHc5HJ+02arGrJV4cnUkGUsZ19VMUnjrvobFJeRO34wShn2hWYNU2Fu5IvlMbSIQDWhXVYRAYCOEOJTmdfsUnpXGKwRAVMoTocaTqMjGXDDjAjiD1kSEWR3OzbYjLJNBVSH+qLqJaRjFkeQbwzzNWqp32V3zaPJVPYKOHSRbA7EZ9BNmkicHouGpdHOW1JrcB71Klk1nyR5nnk/GGWxuWfyRgRVo+KnSbHUVzWUTXPZFGxtEZvcOJkcQb4qq2PZsKn7yRLIHctQHVuE7koNKYMfN6xF78u2pSop2wy36BjjmdVJAqLJaRW67CKCM0T0AkSdwUT0nQCemdld7TFMEUGs+BkaERjCUW7Vj4L6QKgPyFZ/FOQIihhJINm960bPJ8Hguqb6nfV6plTURDj5CHmf5qiHo66aT0XZ9ZAUzEbdPHs4XmOIzPojdxmvvI6FI7BKYNjJYtd0uOzO1KZYmojumauGsmmuRhRVZmv9bX8fW9+BqXrMJeEyGOUH/sh5yyEcgSSLTekb5SBUaqhSISxm+n10bPfHuYgAsDiCft4RVCqE9YXGZRcR/DCAXwXwYiJ6CsCPAvjBWd3UXqNjqEboBu7iFRYa1VwVBLfqRyHLEWz107sLF4rmz4Fk16mnQVipoSoZUzX+dWby1puqqVUwzOycbYNYUtczaCmNLCMds+tMu3vXOhOvMJpMrTOV1XVyEUGB1JCPcG8ZJoPZtKGs5bOW6q64yc3webCJ4JlSSaaeC5eEi0kqvR1FPlnDbvtuL7VqGE2EUYHXtGu3ydcDwFYuIjCXqNvODcjKJFel0W6C21D2KID/hYgWAFSEEFuzva29RaNWQb1KqZ38jmXX4IOpiaUoWax2IduDMa5Z4dUSqx2ZDpcKqA6VFkhHBLzqn6lASqKDxy3kjeV4yjDoVTJWN3HIYpOB8nELpgogW847WVPF0JD2cnMm5oY+ckiIm4QChw7RObkmn8qzNQBap61Fhj3LW+nltroRL54a0iICB1lscszNehVC5O8lUULN9xEAUqI6G+1l5x0A8ntqGzaz3R/H51PHAuaIQH3Ps5s9nyzFboJbNfTviGhVCLEjhNgiojUi+ulZ39xeQhpwLSIw/OE5MDWxxOdiksWmqiFuashEbA4DykeBdKmkL80AJA5kZHAg7l2zIaU0ZqRqjKko/7rC3ILpPfVELkan47mWKfJQqStboUCjVs05KXUO1WCYhUwn8Rrr6lWp7WSa2WzUMooH7PC4GFu5qakqS8ldmKqSZERgi2bSx7tSQ0Ay2F7H1mCMZq2SemanIxhYUkOGVFJWME9BVRrtBbipodcLIS6oH4QQ5wF880zuaJ+wmCGWktF34X0EebKYJxqnkM0nBqWGauY0DYcsNg22SRyBm7wF0k1XHAdia4by6SKZDDprZ2+MlopFEr51TWO+313Ga0vrNT2lu1a5Dcv7aIuM1D3oICKjbIQtNdY0bCbkPZkdtW2AvakqqxYNvjER3f3RJJd2szXDdUfm7+OiI+LI5vyBSBPM4AjGkym6w0lKHkN9f00qx7bU0OXoCKpEFIuQEFEbQNNx/BWHTnYYzKAoR1DDYDxNjbzrDeWHmlPLD8gvaqNWiSsS5Icwr2digs2Y8Mji/M4+NhDOChmyr3M5EIukBWdnX4Qsrpu4jIJkse8+6zVTmmcS3o8xdpf+GglwT6VRaDdyq25OJZkjgvymQAhh/ft0LBGBrSrL5JSEEMaIwDZDuRdzBNm+A3t5arYKCFCbR/tAKxNH4KwayqaG9tARcLe7/xXAh4no1yErh74XwG/M7K72AdnxkHFEEJgaij9MowmWoy+imlfMnWsAqLBzhMF4guFkmuo6dMGW/ggji5NdczzkxFHXb3IgLPVRrW9BvTec2QnS2eXJYt97ZO5knrBIZuOMB8+zbU6y+ktuFVhb5OFypjYCvEJw9jhMBVKCZi7FUpsjND2/Ob2oeKb8uU0RtLofsyPIVzzJz1A+DWnTSYo5u6wIXmBEsNg0l5tuRaWmS8zU0FZ/jEa1kiO7V9r1eMCNb8zspYJLFv8sEd0H4JsgB868UwjxoZne2R5joZFuRy9cNaQ+TINEv7074E8nU1CpqmzZmg+NaiX1QRbCPTJSR0IWJ1/MkNSQ7kB8GjlAOhWlvgSc8lFTQxnLgVi4BZNwXPY+szu50SQ/pjG7JrR5zZjy8qWgLMSvL/IApLGtZQy3URDPEnWY3m9TROCLNswNZROrsmu+gsmsBKv+rtmIozeSzj9rXF19ClsGnm6xZR5oZdrhd+pVEJmrhuRIzvyzKl2jrf4oFrmbFdjbXSHEnwD4kxney76i06jhfDfR/w6dV5ycJ69bJGcRhJ1nsVnDZm9krDF2QaY/DAY5ICIYpiICd9UKYK424paPymsIqMfjdfrmUzysXL+FI1guEEn41jUsMxP8fQThRj3LRfichy5OGMnlOzkdI4ltSQ01DAq4rs+CNTU0zvcFAFFEYBtvmjneVvE0GE2NMh+u0Znb/XE8lEZBzQ3RI1p1rPq9QqVCkaowL+0EaFLXvdk7Am7V0LcT0cNEdJGINoloi4g2Z3pne4xstU/RiCBuetE+TF2DMqIP6wtyJqrKEa60+RyBqXrHR6QCspwTyJePuqpWADNHMJpMUa3Yu2j1e8oajUIcAWudmUgv0lnsrRqyOA+fgQ4dM2qSpfA6HJPjdmwYTLLfttSNaViOKyJo16sYT0W+T2OU7wsAFEeQvxfA1NNgKX2dmM+dzDk2axllOYXFZg3jab7voBsXmmQkL5pVoyTFloUD1B3BrMEli38WwN8TQqwIIZaFEEtCiOVZ3thewzowPnAnnzS96E5lgk497DxrCw2c747i0XlrC7wdQbZb1ydRkFpraCjjpFzUFzzrgDipmuy68ZSxrlrBeJruXi1aPsrRYWpaq5tcxtbUUOa+R5sMuM95DMbTVNMUh78AMsbasWu3poYcHIFePjpwnLsddwtntIzGE8v0N0MF08icGmpbOAIbEe0ii3uG6W02NVSbRPeCVfLCLGWz0rn8HMFzQogvzvRO9hnZdvTucIxWPZ9H9CEhnLSIIHD2MQCsd+o4tzOMh3ysM0PDZi4i8JO2CrbyUV80Yas24pRzqmNDrhc7nml6R+uvGjIbdFd5prqekSwtUhrrmeEcGumoe5hoTpHT4wBkjLUrIjCknwaTKRqGXbVpJoarOzqRojakbyxksX02RJYsjspHuaWvtQoqZCaLeyODjLZKJVkGEGUbSG1pMKmTlH8vFce4F46Au009TkT/H4A/BBCrIAkhfn8WN7Uf6DRr6I0mMUMfOp1MwTTbYGc4wdG18IjgYm+E01vJEGwOZENZYhRcuzHTWiCf4uHstIE0t8AxzMYUBSdVY+AWOL0SZlE3XgNbXp5beKt5QktO1XX0nPPQ141cS94LZZc5BDOQIXQdkWOzVsmRnFayWCOi42NdqaGGmdAdGMZRAjLvn530ZeUILKmhgSWaIZJ5fFNE0DcM1oknrNkGUWUdQb1mdjKGEaHA5ZkaWgbQBfB3AHxb9N+3+hYR0euI6EEiOkFE77Ac8xoiuoeI7ieij3BvfLeh1AHVB1LOGA7bxQPJLqA7yHAEgVyD0iN/7MwOqhXil49mGso4KqAKph36kGEojQ4kyKBnHA+jH0BdI7WuCEfAbWALrAAKqbTRryNEenfPqTRSz5Fcx/03axp27T6OIP8s5l21MSJQMxWMEYG5ZNNePmpPDWWVUF2OwKb51G5UU99dQP4NTFViCR+Yfm96lvLUVqOKnkGGum+JCPbSEXDLR98aemIiqgJ4N4DXAjgJ4G4i+oAQ4gvaMasAfhnA64QQTxDRkdDr7BZi4bmoTKx4RJD/YO8Y8os+qAjgkdPbWOvUnbLMOrJ5Zt9MXh221BCnGkcdG7KubnA8o4lA3fOsZk0k3vUUt6DeT5bDsshF+IbnTKYiVQPOcR7q3HpZZ3PBlRoyNPMFVA0phFYN2XbVSbShKfk6IwK7sTapuzYNfQT2iCBKDeWIaLtg4EKzFnceK8SpHsuEtawT61lSQ+16Bc9dNMtim0qRW/UKGtWKVRp7N8GydETUAvB9AF4CIK6hEkJ8r2PZqwCciATrQETvA/AGAF/QjvmHAH5fCPFEdL5TQXe/i1B1vKqRrFvAeANJ2/pOjiMIcyoqInjg2S1cu9Jmr2vYOIKAiCC707bNAVZQO+rsDpNjYPV7VP/mp6LkOiWhzI5cplM0K9X4PrllpykH4uMItGerVpIeCWdDmd6PEWUCuRFBNhpzpclsKTn9d9lrmLgLV0SQSg05mtWSks2scTcba9nlbFNCzezAbYJ5k6m1HLvTyI+3jMlf28zlXOfyBBXKR+GdRs06jc3k9IgIy21ZRj5rcFND/wXA1QD+LoCPADgKwKdAeh2AJ7WfT0av6XgRgDUi+ksi+jQRvcV0IiJ6GxEdJ6Ljp0+fZt5yGOLSMU3oLbSHAJB//FqF4l1CPNcg0KlcHamNbvXHuG6N7wjUlDEFV8dobq2VLGYa5oD0hL4u23zkreLJOBCbhHIWOrcAaNIHTIdlKqv1XSubh+eI9w1SDX3+SiN5bj4/E9/bSHfc9jJjK1luIVwBc0OZybDbdtW2Wv9WzVQ+mlcqBWTtfqNWyfMPltJUQEnNZCKCyEnZIoIs0d2LxnJmS66dYzwt9mG5VcfmHkQEXEfwQiHETwLYEUL8BoBvAfAyzxqT5clOfKgBuC06398F8JNE9KLcIiHuFEIcE0Ic29jYYN5yGJKUjsYRFIgIiAidRjUuE4uH0gSmma5f60BtxLONLC7oKQkgjCOIc++phjThNcxZA6uu6yeL89djlZ1mrseRypbr0tzCZCqlCXwOq5lxBOr9LZKH9wnV5dZYUjDxM1m6wZ33ZpgZ4FIsDeE7zDyTq3zUXnlj7iPIz1t2jQBtGSQ4XLMkOo08oasa2PIcgfneu8NJXBabPT5LLNv4B4Wldv2yigjUnVwgopcCWAFwo2fNSQDXaz8fBfC04ZgPRvLWZwB8FMBXMO9pV5GkdMbx/4twBED6wxQPpQkknpXwHADceGiBva5eS+eMg8pHLbtYXzWOrdrIW8Vj6EINKVdV6zgCd0C+b4HLn2T7KzhT2+zvicuomxyxr3xUpeXS5aO+0lZ1nH5vAL+hzOagKpG4IleJ1lo+6iCLx1OREnW0cQTqeCNZbHUE+YggGXafXuNyYtljgWQam97zYeMfFJZbNWz2Lx9HcCcRrQH4CQAfgMzzv8uz5m4ANxPRTUTUgBx2/4HMMe8H8HVEVCOiDoCvArAv/QpxRDDQOIICVUOAdCpxZBE5AtMOwYfX3no1AOCrbjrEXhN/yTOOgBMREFGusmY85VXj6NcCmJ2+tfw6bvUPkH9Gf7VR+r3hDM8B8qkhznuaJWQnU4GpJ/qwOUaXQzX2cHib3Vy79vy1sv0Nk6mQjX+Wa2RTSU6yuJ43pjY1UcBMANtSQ/J4iyOw6Et1DOWjtgax2BEYHIepgbTdqGIq0pGYjX9QWN6jiIBrnT4czSD4KIDnAwAR3eRaIIQYE9HbAXwIQBXAe4QQ9xPR7dHv7xBCfJGIPgjgXgBTAP+vEOLzBZ/lkpAITo3j/xeNCPRa5FizqECa6Z1veAn+ydfcgJcdXWGvye5eh47SPROyQ1iGE4GOYXeTuqalCsXHsWSrjbgCeVmS2Se7nL+eyKwLa3xj6ShZnAdrTUhqyNLM52qSiwndjAO2KZZm+xsSgTrH4BstVeX6+8Tl1poxVe+ZqcRTLwlVhG/flRoyVRlZ0k6AtANZiQmrI7BEM93RxGjY9ePV9W38g8JecQRcS/d7AF6Ree2/Q+b3rRBC3AXgrsxrd2R+/vcA/j3zPmYGZbR2hmOMJ1MMxtPgvL5CR+tS3ok1i8LPtdpp4LYb1oPWZPPnyW6MWX5qkKgIafBSGHqkEfR7VaE9N1WTdSBcYb2s0eRWVGUNdKLfxE+/JNEHQ3spIEIqUvKreABdxdMVwTWi/obxVDbR+WYiZ6WxXX/XZq0ConRlj61TWH8tffzEev5WPd/NK7WG7KmhrGHvW/oC6lU54jZXbjqcoG2IOPRU0mr0mro3Gx+53K7tfx8BEb0YsmR0hYi+XfvVMrQy0oMAlcPvDiea8mix1NBCsxZ3BKsPVdFzhSKbpnFpyJhgKj/1ORHzrtSv85+r/on+7y1XLWjQsyJ3rry4jmxJJEuaO7PGZzz136ljORFS1vGr9aFksa0vIHtf9WolrmpypoYynwX9XnUQETqZucUx+euMCNL3Xq+aBQ6zqSFX2kldczCepkqFe448fsswYa03mmBjKd8pbFI3dZ0bkBHBcDxF31JiulvwbVNvgewgXoXsJlbYAvADM7qnfYEq+9wZjGNyRml9hKKthZfxmMoCqaEisOazGWQxIL+suXJOj4GtVuRc26zoXCjpyzfomXVjZtWQxfGEdkBzIpBmAWeVrU7iXMfc++ERqrOQxdaJZpojWGjqFUYWYxoQEQD5ucXKcJsjgnxvgKsctFlLz5JQQ2ys9xLPXJ7myGBTiacpgugOx2jXO7ljlSHP6pDpv8tiOeou3uyP9s8RCCHeD+D9RPRqIcTHZ3YXlwFU2Wd3GD4MJosF7cORjLwsdq5QZI1WSEMZkNfjYQ+1qWY7mhmdvpYqHm+5aiyXrdJfE966zO6Zq8OUTw2Fk8UumQWFfHWS38GZCHdfRFCpEGoVYnMR2c3FgJEaMklc2zYG2eE0rtRQKx58k04N2Xb4WfE/2xCb7Pn12n71XTYZ4rYh9dQfTa1OQ507OdbuZADEMy82e2McWTIesivgWQfgjUS0TER1IvowEZ0hon80u9vaHyxEM0jjiIA5AyALKWmdjgiKEs+hyDYYhYjOASay2F/Omazja+SoNYDJwPLKTvOEeGhKSe1UeWqncZTFMOpZR8Da3eeuw+EV0mtc84Gz18qX7dqNaepZHIZavp7u/o3Jdcu402zJpstY21JD9lRPOsJ1lZoC5iomV4ln21JlZDzWEBH0fGSxFhHMElxH8HeEEJuQaaKTkB3B/3xmd7VP2LWIICofFULEEYHN4+82st3B3GarZD0FkZUKMpLI5m2ZOfu4MYyZGrL0StiMQf56Wcfj/ttYid8Qo85wcrbrmOSes2vyzXV+55ZVCPVFBFy+wzQlr14lq15Wq15NTR3z9QXIY9LH20Tk8hGB+7NiEqrrjSaoVsjokNuGBjfbfAFT3wGHIwAw8xJSriNQW+NvBvDbQohzM7qffcVCJDa3dYkcQachJxcNJ1PsDCVpys3RXypyZHEgR5D94nCqhuR18ztM3zWJKHW9YI4gUz7K1hoKbESzpYZYc5wDUnT25+Jfh1tBle0NcEmCZLukfdeoZwsOfOR1hlNQ/86qiQJaaih1vD01ZItOfKmhlLEeTo2SEUC+E3k6FeiPzIJ5bYuTAYCWpUR7pR2lhmZcQsq1Tn9ERA8AOAbgw0S0AaA/u9vaH6iyT+V9i0YE8aSjwQS94XjPogHAzBFUPDOHdeS+xB4phXhdTv7aLxUhr0dxlyi358FGMnOrf7KltaGNaCyOIFcaG0AWBzxXrhIqgPfI5vG5HIGvb8NceWa/n3xEEKWGDBFBTBZzI4KajSPwRQTJGpsonDpel5aO5SgcEYGeGrKVpipcVhGBEOIdAF4N4JgQYgRgB1JJ9EBBDZdOUkPFIoK4S3k0wdZgXKiZrChyO0Tmjl5BNxCyfNEvFaGuG6LFE6/TrsfZaQN62WnGoHMjgkwOnlvmGtJQll/D5xWyu3uOUF1oRGCK/Gx/52zHs6tuX91vVgQvKCJg9BGkj3dFBBmOYOR+f0ypoYFFMgLI6wf1HFWCqts4qHx0jzgCXx/BNwoh/lzvIciER78/qxvbD6x2GvjCM5vYGsgxlUXTOUqaojsYY6s/Lkw6F0GyW06MJJcfUOtHmkHniLIB6V0gR4vHdL3Q1FBoZ3GeIwjUGgrgMuwpGwbxG1AWm5UU4abJGrVqLiJoWQjUHPHtSa/UM3zRcOyODpsZaemkfNTdWawwGNtlpfMRgfvebekbm6FuZ3ogbF3IQJL+yXIEjWrF2NEt71POJNjszTY15Mt9/G0Af450D4GCwAFzBIcWGzi7M8TF7qhwNAAkEcHOcIKt/qhwiqkIktLKxLiGODS9+odbzgkog57pZmY6EN1pqXO5kO1bKN6RbG90Sq2LDeEkWuevNsqXj/qfzVZp5Hr/s1VisTyD973IT7JbtnxOk9SQfH5f+Wi9Suy0k7rXbBWQ7RlMMteD0RSHFmx9BFKkTg0IcqWdAHtqyOoIGumGNVdE0KjKGeipiGA4sTpgQJtJsJ8RgRDiX0f/D55QdiXi0EIDw/EUT5zr4tACb0awCbpu0VZ/jKuX964JO1eLzpgLkFqv7ea4hlkeQ/mUC8OB6NVGY2Y5p7qnoiTzKLPj9pedplVBk5SSPeVXhCyuVcxEv+v+qhXZURsaVdWq2V377pWP5jiCsbsEWRK6Bkdg5Ajka/1saogRzbQbVU2y2ny8qXy0N5w4R1tyq4CIKNd30BvaZxEoLLdmLzznSw39767fCyF+fndvZ3+hpoI9fGoLL756ufB5VDSx1R9hsz/CzUcWd+X+ODAZIE6KJllPBlI0zDCHOpBQgy7vqZIrZ6yxZaizOkzMiCAg9ZUlcTnXIpKDVJQYnM/g6tcqkibjDrPJS2y4IzA9OvSdG5CVOoNMTh4wp29qVakAkCsfdXAEQOIIkojJ31Cm0B9NsNoxbwzb9SpGExGXWdvmFSfnzzsO27EKS+3ZC8/5chaql+0WAK9EIiP9bZBKpAcKhxblH/vM9hCHF4tHBKsd6QgudEfY6o8vKc0UiizhOGD2AcTrC+y01XVVEx1X8kGdu0gXdKNWwTiaTBZPKAvkCEJVS0Pel7g0NnB4Tjo1x490ikhyb48T4+KaYWCrZrJpQkmyOJ12cpPF5ojAljKR+kHp422GPXFiEwB1f9WQofu3N5rgGgdHoI6vVyuxAJ1tl99uVFJd1K6KJIXl1uzHVfpSQz8FAET0pwBeIYTYin7+NwB+d6Z3tg9YX0iEog4t5kWjuFCO4HzkCJbbe8cRJKWOSRqjKFmsDJKNyMqvCyvLVOsSbR2+ZLZemcLp9AWSztZQx6PmNITep05UhiidFtFCih1/gHPjlnhmy259zibblOgblNOqy89Aksd3P4OsBEpHEL6IQJ3TWzVk0DLqOUZJtrQmsaVWPTbyNuPeqddyXcve1FC7jqcu9JzHXCq4FuJ5AIbaz0P4J5RdcdBz+desFM/rt+tVNKoVPHOxh8lU7G1EYEgNBZHFNYNBZ6WGTOkJhjSFbviYDV5Aum9hPOX1SiiNnSLNdrr0BqePAMjwJsz3Uk+V+UjZZE2e4A9ZA7g5gqSaSZHl7ufPOhk/WZwpT40Mu6mBC1AD7JPz98f2QTPZtJavs1ipmGYbyuyGPT1XwFcO2srIafSG/tSQ5Aj2NzWk8F8AfIqI/gCyWuiNAH5jZne1T7hqOYkCQsZDZkFEWOnU8eS5LoDijWlFkJ3LyxWNS9bn0wzcFE/IEBaFRjUZLh6yLh1JTFlRi1qXOLoo4mE02+n9DiEdybneAwYfkauG8jmcWgGOIFfiyW8oG3sionq1gqlAvMN3KZsCeaE3V84fUFVG8jMjhJCDeBzqo+r5AL1ZzXw8EaGVqWLqe6qGAKA7GkfHTlOv546vV3Icgcog2LDSrl8eWkNCiH8L4K0AzgO4AOCtQoj/Z4b3tS8gItx6jSSJvzxgKpgJq+06vnRWOoI1C9E0C1QrBNJKK4eeio0sVFmhEML7hU+vy6dBODvtSyGL9WojtpaSdj2VsrDtPG3XG02mIEYEoqfZuCW1Uu8/7P0wrfFHBFlNKfuGIVc1NXFHYKFNjcqIDzRjbTPUQDIzQF/jkphQ5wTgrRoC0oSuECJKDXmqjKJdvnJQJnkMdXw27eTlCNq1eCbBrMDeqgohPgPgMzO7k8sE733rK/HomR0cucSSz7VOA5/60jkAwOFL4BtCIfPZCUk5nEyxZJifakO9mkyjCs71B3Te6utCd8DZ640mU2/FkEI2B891knW9SikiP30ORK8A4vIYegqKM8zGtsYbrWjvu0+xNCuX4esU1kn5VlRV4/oMJXn8xFhzIwKfIzD1ZlTIHQXqw2xGE9mD4KoCApKUkE87qN2oonc+PDUESJmJWc0k4OcMCoCIXkdEDxLRCSJ6h+H3ryGii0R0T/Tfv5rl/XBwZLmFr37+oUs+z4oW7h26hAqkIsjuXn3lh6m1WtURN68NyPTEeBrWGKauF1KWGV9PU7gMSX9lq5S4/ElT405Gnk5ZhSyvwOExjO8HQ147VGKilnn/5P2ar5MTB/Q8v0ncz50aSjdxDcb2nDyQlo3wpXpyZHGURnI58ZamKOrqFAaSFFDsmEYTENk3M6byUX9EMHuZiZklr4moCuDdAF4LKV19NxF9QAjxhcyhfyWE+NZZ3cd+QR9Vt5cRAZAnKUM5AkB+2UM5guxAFY4DMkpTFOAk6kGiesnw+hAHoshSrgPJGnVeqiy8oS9F8LP7CMKcR/Yz5YrAspIcMtpwNZRlIgKHdhAgjemFrqxd8aV6shHBYDTx/h1kt3ByvHrNhGT8ZEIW25RK1fHZHgV/Q5k00xdnSBjPMiJ4FYATQohHhRBDAO/DARSqs+F568moOlvr/qyQTdOEdhYD8ssbWtdfhCzWq1dCyNt0R/KUlb6S16NUmoD73ujGlptSyhrbUK4lpDks1xvBSSflqrXsa7KlsJzUUCryCowIfKmhcI6Ad25A5vezEYFLawhIRk7aJKj141XVkBB2yWodexERzNIRXAfgSe3nk9FrWbyaiD5HRH9CRC8xnYiI3kZEx4no+OnTp2dxr7uOGzRHwCEjdxPp3XJY1VBD+xJzZaHlOs3oBZaBDjMGI5y8FSznIe8pvY6bNiviQLIEOiuyqulNaFPUKvaBLsm9FRPuC+mwTjsC92fKJFvtJIsz4yf7o4m1CghI5/B9oydNOkm2UlMFXTbCWw6aEanrjyZoOd5HJTEhhPA2zinshRT1LB2B6dMrMj9/BsANQoivAPCLAP7QdCIhxJ1CiGNCiGMbGxu7e5czwrEb17HYrOH7v/amPb92thcgpI9AV8BMKlB4u19VMhgmMaGllMZ80jcdSYTu7JPrhazTnWuwUWdqPqV5hYA1gRFBvVqJ5cI55H4It6LOM9YcmlNiomaICBzG0RQR2IxpVrbaVWqarElSQ655xUB+6lhvNImbzExoNaoQQt63T45CYXkPhtPMMmdxEsD12s9HATytHxCNv1T/vouIfpmIDgshzszwvvYEG0tNfPzHv9EqjztLSB0ZfScaUD5qIItrllmz6Wvm17E4gkxKiZ2zz6aGCvUR8HWY6tVKPKeC60CyEQHnWo0aZXbSzBRURpbCFyHp4z45iqWpCihPakz/LEynwuvQshGBr2pIbyhLOAJfRKCVprI4AiZZXE9zBP3R1Fo6qh/fH03iITZ+iYkrOyK4G8DNRHQTETUAvAmJVhEAgIiupigPQESviu7n7AzvaU+x1KrveVoISO96gzmCavLF4erWyHVJyWBI+WgqxcMcZqOul05/cQ06FXI8+mhHrgNpaI1eXM2n1O5+MnXOK1bQneIg2n17S1sNjts39yCUIxhOpokOVFBE4E4NpctH3bLS2YiAxxFUUqkewE4WVytSKLCnHe8if3V107j5zOMIWvUqGrXKlVk1JIQYE9HbAXwIQBXAe4QQ9xPR7dHv7wDwnQB+kIjGAHoA3iSEyKaPSgQiV1ETUj6qVXwkmjo88haQu+UwByJTSuPJ1Dkly3Sf6SoWfkSwHYnj+coas+tCK7HSJZdMjiC3hldKG/+9maWtemUPx3GnNhc+jkArH+UQ3rmIYGwfkgNkIgJuH8EkiSA4VUNxqmfoN9btejKlrD9yzxfQx1XGs5k9HAEwe5mJmeYthBB3Abgr89od2r9/CcAvzfIe5hHKmMhRk+ESE4A0dmPGTjG/TgQ3hiXrQqp/0lU8nQbvo5zdcXNTd3oKixtlFXEe2XkQnPcjVWk08ZdHAolk91iPCFypoYyDcqaGavnPgnsegTw+XTXkjggmU/l58U0cyzbDDcZ2SWkFnYz2kcXqdzpHsOKYSKh3IqvPE6dJbNbDaWaZGiqxT1DGLqSeX0HX0Q8tAwVUtVFA1ZCeRghsDFPXGU9DUkPpfDqbLE6RuLwmvWylDVfcTh94HyzlPea9F3oKcMCMCPipIUPayVMFBOgcgTuPn0QQU212gfn4eMaDFm14U0ORzLWSlwDsncKAiiAUR+Ami/UGNJ9SqY5ZD6cpHcEBhNq9hnQG62uBxDDL9YEcQUQmcviRZBD9NKj6p54hVdmpoVraaLJTQ7XwktysMB43xZbwOwHX0fgLrsMBMhGcpxKIy+WY+SLHWE+1a1fG1FM1pJdsuqaZKTQzvRkuHSN5/sTRcIy1TA0pjoBHFvc0stjHEQCyl2CWVUOlIziAUNOnQkjbZG0+v8syYPouMLBjN17HNJbqerE43jRQdK4Af5InS3lOjruLNl2HnRqqVTDS5D3CUnK8PH5DK9n0/X31zwKnR6ESEa798QTjiZxL4EsNAVFE4EkNAdJJhDaUATJ9w0oNaVVGfYdAHaBpEw0nMf/Aiwhq2JphRLD3tY0lZo5GrRrcGZyszeTsmTv7rDQFt3ch27cQ4kBE1LcQJDqn7WwHzPJMIMMRcPP92dp75hrVjzEaT9FkOY9EMZbPXyQpQE7vQVasz9VbEupkgKg3YDT1kr9AOiLoe1JD6tq6DDWnfBQA+mPpCOpVcr6nsltYI4tdEYHWd6D6LDgRweHFJp7d7EMIMZNKxDIiOIBQXbBc6eP0WpUamgTX9at1YSkePZLg5/p1hzUOFp3TGp1COILAdY2aXhHFbA7L9HHwBPgSxVguf1HPXEeex52+4ZbdxufW9Kp8zy4rgSYsR6CXhLKOz/QdeKuGUrt2vyicFJJLOAXX8Yk20UTrUfD/vV5wZBHd4QTPbva9xxZB6QgOIJoxR8CbVqWjkfoSh0g3qB2mwDAg957MAxYYTcNTSjonwV2XarYLiFxUJ264URfsNI96jsGYz5noRp0bEaTKhDkRQTVNfLuaDE0cge99TiICf06+qZHLg/HE2zch710ni/2GHVC1/n6ZaCUkN5oITIW95wDI9hEoItofEbzoyCIA4Lc/+YT32CIoU0MHECr9UYwjSJO+XEOZbVBir9M6XMN26Mm6oM7iXOcu9z4zpDZzp66uw+0SbmaMOlfuGghLyyUyELzIMS9b4id/dY7A9z43I44g7hT2SEwAkpj1dSGrcw3Gsst5yIiYsrt2nzqoIov7Y3+aSncy04jXcaWSFG67YQ3f9hXXYmVGQ65KR3AAoapIOPXhWew16Ztax9xp56/n3qHqyIrVhXQyA3JHyZaYSFVghZG4MWcSkOYZTCYYjnk9FbXAXXuq0S2ofJQXlbbqVQxGU82Y8spNOTt8Fc3EUhqeVExHa/riDI5pNyRHoCqMXI5DzmIG+sMJpkJ2JnO+K7VqBb/45q/0HlcUpSM4gFDEXsjweYVUQ1mQ5EN4ekJfF2L4Utcbq0Y0vuNRefuQiCdb5srbqWeinYCderHrSEXL8PLRsD4CHydj6lrmkMWpiIDBEfSjVBInIuiPeBVGQLr7lzM4RvUdqIY41w6fiDQFUilnsR8yNFmUHMEBhCJE1ZeqEFk8Dsu9q1TBUOXDZ53rz+y268yIQF1PEXVcSYusw+IOpknWhPEKwzGfPM/KexQuH/VUDQ1VZZI3Isg3JXIjAo6xzkUEnh2+iggGjNQNoOfxxyyOoF2XVXpKusTnOJQj4KSd9gqlIziAUF86VdIWJkOdGPQgmeaq/ECPAgwloO+0RZBBVwa8P5IhdmiKZ2cwiX4O4whCrpfd3YesGV5KFBFASqdVZl15/yTqkM7GfiwRoVahoC5zJS3tE5FTxwJR1ZCngUseH1UkjXhOSaXWukyOQKWS1NQ0Vx+B/L0cTtMf+dNae4UyNXQAoXaEaocSQhar+bShpG9K1piZp9bvTeWTQ1I8AOJpTyHD64Hw90YdF+JAdJ0bn/FUaGoRQYimUXydUImJ8ZSlWKret95QpjR896V6KLhksdL3GTDSK8nUMVk15IsIFFnMKTUFtFp/Jkegqn7Od0feewdkRNAfTUCgyyYiKB3BAYT60obuehXUPINQ7R9AM+gBstBAsYYyANhRUU/gutBoKX5Ph8qBcKavJVEEx3jq96dKTrn9CnINPyJQx4yjclhf41oj8377JD3kjOOEI/Dm8WsZY+2UmNA5AkbVUKY01bcL16uGOKMklaM4F0UEvnLQdkNWGVWIWD0Ee4HL4y5K7CoSR6CMXRgZpfR4QuvzgSilVKSzOLB8NM71D5WzC7tP5STZ5aq599S/ThnXOPoI5hUKEL+BUYT8O0+891bPPL/v/c7rXTEjAkYePx0R8PoC+uNJ4pQ8xrderaBeJXTjPD6PUzi/EzkCxv3EHAGjq3gvUDqCA4h6xgCp/D0Xqos0VJ0TCBePS3fS8iMQ5dyUYeLPI0iv46aisuk21tS2ApFZvQD3YUrzcK+jSGmfQ0wigsiB+hxH1LjHJYtVRNAf+cniepVQoUR9lB8RROdmvKeqN4CTGlIRxLkdxRFwyOJpNLugdAQlZoTs7pVr7BRkzbgITNVEHMGY3wgFJAamP5JiY6HkdFL9E3a9JKXE+yJm020hO/U4igioAEpSUAFRRCQLEjyYhlF+GxcgMLmVerWCMXMeMiC7hbkRARFFc4UnLDVRRRbHukQM49tp1LAzGLN27cqYn9keAAAWmgyOYDhBlyFfsVcoOYIDCGVwQoxJan2UGgrR/qlWCER6QxmXZFYpnjCnVY8jgjCyOEsyh6eUAjiCagGjno3mAnb3veEUQvCcYmi/QmgBglJ55U65a2UjAk/6pqURwP75ArJ3RH1WOFpM7UYVFyK1T07OHwDObsuIYMFTKNFpVNEdjUGgfZlpbsLlcRcldhVFdq866pGaZYhBJ5IKjcNI357PEUgDsa3SJ6H9AIHOLpfrDpTQCCGni/AKsaMa8PsccvfGuE61ItMr3AbAULJcpYZUdOhrmlK79K1oCpcv768iAlZDWXQuNeGLQ9C269Uk1cMki1VEwKkykk5bxGml/cZMU0NE9DoiepCIThDROxzHvZKIJkT0nbO8n3lBdvca0lCm1o8mcvB4yFpd44htmCvpKp6QeQRAkrMOrVIKJ4vT3EJQaiguceU7j5Dy1lwqMIBwVxITPmOa21wwHEeI3Ii6/mZ/hHqVUPWIHaqIgFOLrwz/xR7PyQBy1x7qCM7uDNFpVFHx3LsqH90Zji+biGBmjoCIqgDeDeD1AG4F8GYiutVy3Lsgh9yX2AVcSh8BgHgISYj2j7wOBYmyAXIoSa1CCaEaWG2UpHjCdvbB5aMRl7ATcL1cRBCQsikSRYSsUfcznvC4oOw1/OWjiZNhVVhFxvRib8wy1KmIwNdHUFPnVo6Alxo6FUk+L7bcxrqtkcWc/pl2vYrtwRj9Eb/fZtaYZUTwKgAnhBCPCiGGAN4H4A2G4/4pgN8DcGqG9zJX0MP4asW/u8oijggm/IEvap0qe+RUZuj32w1N8UT31WXWtevXAsL5k3omIgipANoOKOMtcn9JlVhYlBOXCTOMdXJfvAhMznjmFxy0VETQG7EMdbNeQW8kK804WkPq3PJnXkSgntW3a9erhHxEcZHj9wKzdATXAXhS+/lk9FoMIroOwBsB3OE6ERG9jYiOE9Hx06dP7/qNHjTotfIhgnMKaghLiCibum4voHlKXxeXJQbn+i+tjyA4FRXzLgFGPSTNkzHqhSKPAH6GywXFkRQzuqnXKhgENLglEQHPEbRq1XiH76u8KRIR6Dv1JU9E0NLO5yOKgfT9Lhz01BAA0zdFZH7+jwD+hRBi4jqREOJOIcQxIcSxjY2N3bq/Awv1Qd8ejIPTQoCW4gkoAwWkAQolYeX1KsH57WyKJ1wqIrAjuYBRz6avihDMHKOVLVMN4gjGPMXSXKTiOz6KDrmfIT0i4JRUNusVXIw6edkRQZ//nuq79sVm3XlsrVqJn5Gzw+80Lj9HMMu7OAngeu3nowCezhxzDMD7ooqCwwC+mYjGQog/nOF9HXjUtd1bES0TPcUTyhGEkrDyWAoqy1T3CITl7PX76jIbo/LripeCcpyjqua5tBQU3xFwBwJlIyJ/lVESbYREBBd6I1y13PIfX6vGht0rE62qhnojr6aSwpJmoH0cASAdx7DHy/nr5PPCHFQN3Q3gZiK6iYgaAN4E4AP6AUKIm4QQNwohbgTw3wH8UOkELh16PrdIRNCoVeLO1tDUUJe5Y8xeLyZ9meuUwVTlo+yRmrW00Qx1ICHrEiI8LPpo1CpBmkaqhyNkjTqOq0+UT3P5OYIQ4Ty1S+emhpr1ipYa4nMEnHMDwNpCMgmMU9mjjPtqxx09AOloY6XtP34vMDNHIIQYA3g7ZDXQFwH8jhDifiK6nYhun9V1S6SNVCFHoOXsQ3P9iaxFGMncDRSPA1QqKmxnnxWd4z5ftropyKgXuMeQKiqlGLtd4DrcxsGsbAlXfZQrQKh29ZOpYNX5t2pVTJijHvXSVK7s83qgI1AOYI0xSvLwYnKM7nD2EzNNUAkh7gJwV+Y1IzEshPgns7yXeYJuCEKbyYDIMAematS1Qss51bFFSl0VOQ0ERAQ50rcoB8K/XmiFUkPjTEKkM7hErn5v3KKAZMYFNzWkcQScHb52DJcjMP3bhJZGRB9Z8qedAGBdM+icqrtDkXHnRAQbi8k9rM9oBnEoZpkaKrFPaFxiRFCvEVtcLHvd0MoVeb1KLC0Q5Hh2gWQOnd4WKnudch4BMxOCK4AKrInnToynXs2lbPcyS310zBcg1I0/Z5i7fgw3IhhNhNdpKKwvhhloimpjDi82vcduLCXHLB/01FCJ/UMqIihUPpp8sULJ4iIRgX6PoRGBUpQMlbQILR9V14iyEUHSFNwB7sk9hpfTNgpEVfWI0B0EiM5xIxVdwZaV809FBDyOIPk3r3wU4DkZAHjB4UUAwNG1Nuv4l163AgB4/saC91idIwjt8ZkVLo/apRK7CtVENpnyNX906EYh1DCPI0sZmus3/dt7Pe0+uakhRawOJ1MQhX0R9WcKnWNQeA07IiDsbIdzBNv9MasBUH2muuyqoaSzODgiYKSGUhEBQ6BOgRsRrHTq+LW3HMMtVy2xjn/7N74QNx9ZxFffdIh1/O/94NewP7N7gdIRHFDUqxQRb+HlaSmDF5jiKbSuYCortY55vVgcjymGpkN/ptAO6LA1BRxBtRKTpyEcQUgDYL1K2GZWhdWrFQghZcJnzhEwBOpM1/HhtbdexT52sVnDd9x2lH38bTessY/dC5SpoQMKZQyKDMdOGaLAnL3pHCHXC+UI4nVM1VJ9XbgYn7y3SkAkUcR5FIkIijjvRk1GBNw1jci4AzyOAJCppJDyUYC3aw9JJalBNnLd5VG3f7mhdAQHFOqLWGQmavEdemIcQklm07/91yvmQNSxoWkzdfysnZzufIukk7i73mYgr9DQjKjP8ep9F5xz16qV2FizyOKAVJIaZANcPto+lxtKR3BAob6IRWaiFtnFZo8NM+hFOYlkh84VndOvET6nITyS0KMPbhqqCK9QJJ3UqFWwFSSrLe+/WiGv1LJKHXFE4RTU+8NKDQU6PhVlXC5qn5cbSkdwQFGPI4IiHEFxsjg+R2B9fnyOAtxCYYMeGhFE60K7puU1wxrsAEmA+wxu9t6y/3ahWUtSPSE6SKz5AgU+C4rj4AxrSQm3cYTeoojgctH/v9xQOoIDChUJXGpqqKhBDyHlCnMEtWI8SLKuaGoonMcoQp4XIeqrFQqW5OZeKyQ1pld0hXIxPrVPIP234zhLlRK6XCaCXW4oHcEBhVI1LBIR7EaKp3DVUBHSNzgiKMgRXIJRD+qrqIVfR0VxsyzbDUmN6b0ooe8zZ9ce+rlW34fLRe3zckPpCA4o1I6pUGoolaoJT2kAxQxSSCpEv16RUZyXsq6IUQ9rsCvARRRwHnokFZYa4h8bek8AsNTyd9uGRnOqZr8IZzYPKB3BAUUlIt6KhMK7Udcftmum4GsBSTok1CgU5ghq4Qa6yLUuJfIIeS9SqTxGCjEuSWYcWy9Q+aTASQ2FbnCi5u6UmFyJBKUjOKCYRixgEZnbRsH8btFGNJVGCB+peakpnjBjEhvbACNUiCyuhT/XpTic7L+txwc43kuJCI4s+/V6Qjc46p4PBWoIzQtKR3BAoSowijiC3egjCGmfV4YvtOO+UWAXrF8vdKca92YUyN2H7GCVYyySGgp5L/SdPbehLPRY7vEA8BPf8mV40VWLOLzgdwQ3HFrAi69ewk9+662sc//ga16Aq5abuPWaZdbx84aSOTmg+MrnreGTj51jqSFmoRv/kKojvawySLohul4IPwBcevloqAMpwrs0CqxRjioozVMgytGNNafyKqRKq0gF2fd/3fPx/V/3fNax1Qrhgz/69axjAeAbbjmCT/z4NwV9LucJpSM4oPjhb3gBvuyaJXz50ZXgtfqXOGwnWyzAVIa5Gvgl3WtHUI939+EpmxCSslmIIwhPJxWtGuK8b0ULB2aJ0gnYUTqCA4qlVh1vePl1hdamd4rFDEvQ9bQa+BAkO+fAPoJLJIuLRARF+iqKEPU+FVEdoWJsIZVJlzocqcTeovwLlchBN0Ahu6giQ3D0daGO4FLF44pWGxUx6rUgslitKcARBEQrzUBjHUQWFywcKLE/mOlfiIheR0QPEtEJInqH4fdvIKJ7iegeIjpORF87y/spwUNIhUt6XVFHoMjivU0NdQKbi4pELspohjxbLGURcJ0i0gmhefwQHqKIXlKJ/cPMUkNEVAXwbgCvBXASwN1E9AEhxBe0wz4M4ANCCEFEXw7gdwC8eFb3VIKHogY9dIetkJDFgesKNFEBya55IbAEMb7PAlFSiPNIIgL+GjXyUFWLcVBGBCUUZvkXehWAE0KIR4UQQwDvA/AG/QAhxLYQSvYKCwD4n+ISM0NRg84Z3G1CnBMP9ASKgA2NX5IUT5gjUL0ZITXsygFwRyQCSZ4/xHmoMuHxhP8VCt21h5Tr1lPRRtnNe7ljlo7gOgBPaj+fjF5LgYjeSEQPAPhjAN9rOhERvS1KHR0/ffr0TG62RIJWQWGutU6xZh01wzV056iMX280CVqnehzGAbtnAAg8HEASPUwFf/Fi1Fk7jOYxc1Ckg1x/vzl8RAiJrTuWUvHz8scsHYFpO5P7Nggh/kAI8WIAfx/AO00nEkLcKYQ4JoQ4trGxsbt3WSKHpYJfXNW+/9XPXw9ad92qHBB+7SpvUHj2ekfXOkHrrlpuRddrBa37phcfAQB8Q/R/Dv7WCw+hViF888uuYa/5yuetolmr4LuOXc9e8/wNOWz9+77uJvaa0J16SMezzjN1ymEwlz1m6apPAtA/yUcBPG07WAjxUSJ6AREdFkKcmeF9lfCAiPDet74SV6+EGcqFZg0//11fgVfdFOYIbjjUwU/9vZfgdS+9Omjdq19wCD/xLV+Gf3Ab32ACwHe/Uh4fer2br1rCiX/7+qBqnhsOLeCBd74uaM01K2184f9+XVBq6LrVNh766dcHRVUcTR8dKu3ESSPp1WaceQEl9hez/AvdDeBmIroJwFMA3gTgH+oHENELATwSkcWvANAAcHaG91SCidfcwt/16vj2V/AHeCsQEb7na24MXlevVtidqDqWWvVC64Cwks5LWRNaSguEp9ZCU3kqVRVaTFBkJkaJvcXMHIEQYkxEbwfwIQBVAO8RQtxPRLdHv78DwHcAeAsRjQD0AHy3Rh6XKFFihlDOhut0ltvSXKy0w8xG2dF7+WOmMZsQ4i4Ad2Veu0P797sAvGuW91CiRAk7fubbX4armCnA7z72PDxzsY83MqO+X/3Ht+FCd3gpt1dij0BX2gb82LFj4vjx4/t9GyVKlChxRYGIPi2EOGb6XZm8K1GiRIk5R+kISpQoUWLOUTqCEiVKlJhzlI6gRIkSJeYcpSMoUaJEiTlH6QhKlChRYs5ROoISJUqUmHOUjqBEiRIl5hxXXEMZEZ0G8HjAksMADrqI3UF/xoP+fMDBf8aD/nzA5f+MNwghjPLNV5wjCAURHbd10x0UHPRnPOjPBxz8Zzzozwdc2c9YpoZKlChRYs5ROoISJUqUmHPMgyO4c79vYA9w0J/xoD8fcPCf8aA/H3AFP+OB5whKlChRooQb8xARlChRokQJB0pHUKJEiRJzjivOERDR9UT0F0T0RSK6n4h+JHp9nYj+jIgejv6/Fr3+WiL6NBHdF/3/G7Vz3Ra9foKIfoEuk5l6u/WMRNQhoj8mogei8/zMfj6Xwm7+DbVzfoCIPr/Xz2LDLn9OG0R0JxE9FP0tv2O/nku7p918vjdHr99LRB8kosP79Vw6Cjzjq4jonui/zxHRG7VzXZa2JoYQ4or6D8A1AF4R/XsJwEMAbgXwswDeEb3+DgDviv79lQCujf79UgBPaef6FIBXAyAAfwLg9fv9fLv5jAA6AL4h+ncDwF9dDs+4m3/D6LVvB/DfAHx+v59tRp/TnwLw09G/KwAOH5TngxyXe0o9U7T+3+z38xV8xg6Amrb2lPbzZWlr4mfd7xvYhT/W+wG8FsCDAK7R/ggPGo4lAGcBNKNjHtB+92YAv7rfz7Obz2j43X8C8AP7/Ty7+XwAFgH8dfQFvWwcwS4/45MAFvb7GWbxfADqAE4DuCF6/Q4Ab9vv59mFZ7wJwHORo7vsbc0VlxrSQUQ3Qu40PgngKiHEMwAQ/f+IYcl3APisEGIA4DoAJ7XfnYxeu6xwic+on2cVwLcB+PAs7zcUu/B87wTwcwC6s7/bYriUZ4z+bgDwTiL6DBH9LhFdtQe3zcalPJ8QYgTgBwHcB+BpSIf+n/fivkPAfUYi+ioiuh/yeW4XQoxxBdiaK9YRENEigN8D8KNCiE3G8S8B8C4A/6t6yXDYZVVLuwvPqF6vAfhtAL8ghHh0FvdaBJf6fET0cgAvFEL8wSzv81KwC3/DGoCjAD4mhHgFgI8D+A8zut1g7MLfsA7pCL4SwLUA7gXw4zO74QIIeUYhxCeFEC8B8EoAP05ELVwBtuaKdATRh+f3APyWEOL3o5efI6Jrot+r/Jw6/iiAPwDwFiHEI9HLJyG/YApHIXcklwV26RkV7gTwsBDiP878xpnYped7NYDbiOhLkOmhFxHRX+7NE/ixS894FjLaUc7udwG8Yg9u34tder6XA4AQ4hEh8ya/A+Br9uYJ/Ah9RgUhxBcB7EDyIZe1rQGuQEcQse3/GcAXhRA/r/3qAwC+J/r390Dm81RK5I8B/LgQ4mPq4Cik2yKir47O+Ra1Zr+xW88Y/e6nAawA+NHZ3jUfu/g3/BUhxLVCiBsBfC2Ah4QQr5n5AzCwi88oAPwRgNdEL30TgC/M8t452MXP6FMAbiUipYr5WgBfnOGts1HgGW+Kom8Q0Q0AbgHwpcvZ1sTYb5Ii9D/IL7yADCHvif77ZgCHIPPfD0f/X4+O/wlIz3yP9t+R6HfHAHwewCMAfglRp/V+/7dbzwi58xCQXyz1+vcflOfLnPNGXEZk8S5/Tm8A8NHoXB8G8LwD9ny3R5/ReyGd3qH9fr6Cz/iPAdwfHfcZAH9fO9dlaWvUf6XERIkSJUrMOa641FCJEiVKlNhdlI6gRIkSJeYcpSMoUaJEiTlH6QhKlChRYs5ROoISJUqUmHOUjqBEiRIl5hylIyhRokSJOcf/Dzol/TskyNXpAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Arrange in a dataframe\n", "df = pd.DataFrame(\n", " {\n", " \"epoch\": t,\n", " \"datetime\": cdflib.cdfepoch.to_datetime(t),\n", " \"distance_km\": d,\n", " \"distance_au\": d/AU\n", " }\n", ")\n", "\n", "# Plot the whole file:\n", "fig, ax = plt.subplots()\n", "ax.plot(df['datetime'], df['distance_au'])\n", "ax.set_ylabel('distance (AU)')\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "b8080591-68f4-44cb-8995-20ddcec2d942", "metadata": {}, "source": [ "### Find the intervals, and overplot" ] }, { "cell_type": "code", "execution_count": 11, "id": "e32a3bc2", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2020-05-06 19:55:50 - 2020-05-22 08:55:50\n", "2020-07-09 05:55:50 - 2020-07-24 18:55:50\n", "2020-12-31 03:55:50 - 2021-01-15 01:55:50\n", "2021-03-08 16:55:50 - 2021-03-23 15:55:50\n", "2021-08-12 21:55:50 - 2021-09-02 16:55:50\n", "2021-09-21 23:55:50 - 2021-10-12 18:55:50\n", "2022-02-18 01:55:50 - 2022-02-27 00:55:50\n", "2022-04-22 22:55:50 - 2022-05-01 20:55:50\n", "2022-09-05 22:55:50 - 2022-09-14 23:55:50\n", "2022-11-09 13:55:50 - 2022-11-18 15:55:50\n", "2023-03-04 16:55:50 - 2023-03-13 17:55:50\n", "2023-05-08 07:55:50 - 2023-05-17 09:55:50\n", "2023-08-31 10:55:50 - 2023-09-09 11:55:50\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAD4CAYAAADhNOGaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACBRElEQVR4nO39ebhtaVUeir9jzTlXt/t9mupOdUBRRQGC1AHF7qJcFEwMQU2EPAkGGy5GEr15vL+gMVGveiNJ9N5rbEjlBjGxIfagloCiSESaKpqqogqq6lR/qjv92c1qZvf9/vjmN9uvXWevfc7Ze77PU0+dvfb81vzW2muN8Y13jPEOYoyhRYsWLVrsX3Qu9gZatGjRosXFResIWrRo0WKfo3UELVq0aLHP0TqCFi1atNjnaB1BixYtWuxz+Bd7A644ePAgu+GGGy72Nlq0aNHissJnP/vZU4yxQ7LfXXaO4IYbbsBdd911sbfRokWLFpcViOhx1e9aaqhFixYt9jlaR9CiRYsW+xytI2jRokWLfY7WEbRo0aLFPkfrCFq0aNFin6N1BC1atGixz9E6ghYtWrTY52gdwWWKOEnx3z/1OJ49P5lp/d1PnsNfffnEDu+qRYsWlyNaR3CZ4rc/8wT+zR99ET/yu3c7rx2HCb7jV/8Wb3vfnXjwuc057K5FixaXE+bmCIjovUR0goi+qPg9EdEvEtExIrqHiF4xr73sRXw0O81/6pHTGIeJ09rPP3kWccoHErVRQYsWLeYZEbwPwOs1v38DgJuy/94O4FfnuJc9h7ufPIfVYYA4Zc6n+s8/cQ4AsNT3cc/x83PYXYsWLS4nzM0RMMY+DuCM5pI3AvhvjONTAFaJ6Kp57Wcv4fwowtlRhG99KX+7HnjWzRE8fGILV6/08VU3HsDDJ7fmscUWLVpcRriYOYJrADxZ+vl49lgDRPR2IrqLiO46efLkrmzuUsajp7cBAF/3goPoEHD87Mhp/fGzYxxZH+K69SGePDNCO7e6RYv9jYvpCEjymNQiMcZuZ4wdZYwdPXRIqqK6r/B45ghuOryIK5b7eOqcW+XQ8bMjHFkb4Lr1AbbDBGe2w3lss0WLFpcJLqYjOA7g2tLPRwA8Pc8bfvqR0/jCk+cu+HnObIf4s3ufQRinF76pGXBycwoAOLzcxzWrAzx1zj4iiJIUz25McGR1gGvWhgB4hNCiRYv9i4vpCD4I4K1Z9dBXAzjPGHtmXje75/g5fNftn8Lf/+VPOFMpdfyL3/48fuA3P4df/OhDO7Q7N2yMIxABSz0fV68O8LRDRPDs+QlSBlyzNsDhpR4A4NTWdF5bbdGixWWAeZaP/jaATwK4mYiOE9H3EtE7iOgd2SV3AHgEwDEA/wXAP5vXXgB+Er7x4AIA4I/vnt3fPH1ujL85dgoA8LufffKi8OvnxxGW+wE6HcLBxR5OOxjy0xkNdHCxh4OZIxARRosWLfYn5lk19BbG2FWMsYAxdoQx9l8ZY+9hjL0n+z1jjP0gY+z5jLGXMsbmOnbstuvX8Vc/8hrcfMUS/vbhUzM/zycfPg0AeOurr8dzG1M8fvrCootZcH4cYXnAh8sdWOxiO0wwiex6Cc5sc6O/vtDFgYUugDYiaNFiv2PfdRa//NpV3Pf0xswn+c8/eRZLfR//8ChPb+xEzsEV58cRVgYBAOTG/LRlwvf0Vpit66EfeFju+21E0KLFPse+cwS3XLWEM9vhzMbv4RPbeMHhRdxy5RICj/DARZBoKDuC9cwRnNmycwSiQmhtga8/uNTDKcu1LVq02JvYd47g5iuWAAAPPjdbI9XDJ7fw/EOL8L0Orl0f4tGT2zu5PStsTGIs97OIYDFL+G7bObYz2yG6XgeLPU4tHVzotdRQixb7HPvOEVy7Lkom3bn9jUmEE5tTPP/QIgDgeQcX8Njpi+AIsmQxUFBDLhHB+kIXRLyNY3kQ4Pw4ms9GW7RocVlg3zmCq1b68DqEJ2dwBI+f4mtE9dENBxbw6KltpOnuVg6NwgTDngcAWF/MHIFljuDMdoi1zHkAwMogwEbrCFq02NfYd47A9zq4erWPJ8+4N1E9c56vuXq1DwA4sjbANE5xdrR7HDtjDKMwxkKXUzuLXR9EwObEzpif3g7zKALgjsA1IoiTdKZkezTjumnspq4qMGvDn20FVosWewX7zhEAwJHV4UwRwXMbvHHrymXuCK5c4f9/ZsbhME+dG+M9f/0wzjk4kmmcImXAoMsjgk6HsNTzsTGJrdZvjCOsDIP855VBgO0wQZTYGc1TW1O8+uf+Ev+/37vHes8A8MjJLbz8pz6Cn//Ig07rPvv4Gbz0Jz6CX/vEo07r/vz+5/Din/gQ/uQet2b13/jU43jpT374gkqMW7S43LAvHcHVqwM8N4PxfnZjAr9DeYL2iswhCAfhih/5nbvxc3/2Zbz7Qw9YrxllswcWMkcAcJ7flt7hiWY//3kl60ewXf/n9z+Hk5tT/O5nj2N7aud8AN7Etx0m+M8ffxiJA5X2+597CmE2jc0F7//ME4gShvd/5knzxeV1d/J1v3vXcad1LVpcztiXjuDQUg8nt6bONMUz5yc4vNSD1+GJ1qtWBgC4g3DF6a0pPvUob0770BefsTaOwvgOe4UxX+oH1hHB1jTKK4YA5NGBLT1012Nn83/f+5T9LIO7HueK5FHCnBL1dz7K1z1ychtbDo7n81l/xxeePGf9dx6FMb70DC8Hvvv4Oet7tWhxuWPfOoIoYc7c+HMbE1yR0UEAcHCxiw5hprnBn3rkDBjjHcpnRxGOnbArZxURwbAcEfR9bFjkCKIkxSRKsdSvUkMArB3JY6e3cc0qd4CPOJTOPnJyO0+yP2RZupumDI+fGeH6A7zSy7ZU99woxJntENeuD7A1ja17Rh4/PUKSMjzv0AIePz2aOTfRosXlhn3rCADghGNT2bPnJ7iq5Ah8r4PDS/2ZHMEDz22iQ8A/uI13KH/pmQ2rdaOQG2yRLAbsqaGtzNhXIoKBW0Tw2KltfM3zD6Drd6xLZ6MkxTPnx/iGmw4CAB49Zbfu2Y0JwjjFN958GADwyCk7B/JI9vzflK17zFIGRKiwfsNNh5CkbKaCghYtLkfsS0dweEaxtRMbUxxe6lceO7jUnakh69iJTVy3PsQtVy2h63UcHAE/pQ4qEUGATYsTvbhmqT+bI9icRDi9HeJ5hxZxZG1gTfE8c44rnr746hUMAs+aShM6Tl/9vAP8eSwd7pNn+LpveOGh7HnsHI94Pa+8YR3A7LmfFi0uN+xLR3BoBkcwjRNsTmMcXOxWHl9f6M002OXYiS284PASAq+D6w4MrU/XRbK4nCOwo4Y2p1F+vcBijzuCLQtHIiKoq1b6OLzUs37/hIE9sj7AVSv2EZQo133hFYtY7PnWhlns6yXXrPCfLR318bNjDAIPt169DGA2yq9Fi8sR+9oRnNi0/6KfG3EjWm7GAnhn7xnHPoI0ZXjs1AjPO8Q5c366tqMhBDU0qFUNbU1jY2NbEREUOYKFrDHNpgIoF6xb7OLwUt+aWntusyi7vWK5bx0RiPsdXOrh8HIPJzbs7ndyc4qu18HhpR4We/aies9tTHDlSj8vD56lCKBFi8sR+9IRLPV89PyOk9iaOPWvD6uOYG3YtZZ3yJ9rFCJMUlyd5RuuXRs6OIIsIuhVk8WMAVuh3phvSaihYRZZ2FTkCArs4GIPh5e4YbapyDmzzZ3ogYUernSICE5nukhLPR+Hl3r2EcHWFIeWeiAiXiFm6QiE/Magy1VZ24igxX7BvnQERIS1YRdnHSids7lqZy0icJwHABSUw5VZ+emRtQHOjyMreicvH61RQwCMeQJBDZWTxV6HMAg8y4iAG9QDi10cXu5hHCVWDuTsdgivQ1jq+1hf6Fp3Yp/emua6SC4RyKmtMKfwDi26OwKAi/m5RnotWlyu2JeOAABWhwHOOZSPCqOwvlDPEbhp/QBF0lNUIF2zxh3C0+fMUYGsfFQ4hbF1RBBUHl/o+dgOzY7s5FYIIh4VHcya6k5bREOnt0OsDflEtbVhgFGYWJVmntkOcSAz6OsLXesO7FOb03x/omfEBmdHYR7xrQ4Dp47vFi0uZ+xvR+DwRc8jguGFO4JnsySocASiEsnm5DoKE3S9DgKv+NMVPL/euG5IqCEAWOzZRwRrwy58r+NUbXS2dNJezd4/kXPR4VRp3cqAN83ZNN6d3CocAf87m+/FGKsI8vGIsRXja7E/MFdHQESvJ6IHiOgYEb1L8vs1IvpDIrqHiD5DRC+Z537KWBt2cdbCQAgInnt1WD1NzxoRlKUqBI1h5wjiSqIYKCKCbUNEsDmJEXiEnl/9sy/0fOtk8YGF4sQMwCqqOjMKcwcq/m9DD53Zrhp0wCyFwRjDuVFh0IWonimXsTWNESWs8vraiKDFfsE8h9d7AH4ZwBsA3ArgLUR0a+2yHwPwBcbYVwB4K4D/d177qWN12LU6KQqcHYVY6vuVkzgwY0SwUZWqEFVMNv0IozCp6AwBRSnpyBARCHkJMYsgX9/zrbj+c+MwN8guEUGZe1/L1tuctk9vFQ7E1vFMohRRwvKZziuDAEnKjNSX2I9wIOuOB4UWLS5nzDMieBWAY4yxRxhjIYD3A3hj7ZpbAXwUABhjXwZwAxFdMcc95RAnPlsdmrOjsJEfAIoqIhdHcGorzI0/wJO3/aAze0QgqCFDRDCaJljo+Y3HF7qecS3AIwqRX1gZ8Nd93uLULKeG9OuiJMUoTHKHszqwWyfkuJf7bg7rdDbhbT0b4bm20MU4si8CuOuxM/jrB09aXStw95Pn8FdfPuG05p7j5/CxB9zWfPnZDfz5/c85rWmxvzBPR3ANgLL04/HssTLuBvDtAEBErwJwPYAj9SciorcT0V1EdNfJk25fNhXWhgHilFkLmZ3ZDhv5AcC+YqeMsmEE4FTmOAqbxjyPCAyn3lGYYBB4jccXer4xmgCEIyhO2oDZwDLGcK40Y1nMSjadtgVVld/PMiIQlVfLg5ojMNxPvA7h4PIIxCIqmEQJvvM9n8R3v/cz1l3mYZziTb/yCbztfXdal6kmKcM/+i+fxj/9tTvz7mkTGGN48+2fwvf/t7vw8MnZxrO22PuYpyMgyWP14/fPAVgjoi8A+OcAPg+gYVEZY7czxo4yxo4eOnRoRzbnkrQE1BGBn83/dRGwq08JA7IyRxtqaJpUKoaAUkRgcGqjqLkW4BGJjUPcnBQjMrt+B8OuZ3z/JlGKJGV5JCFO9qb3qy6HIZywyaDXE+K2Dku8/uXaOpuS3rIiq1BLNeHep85B5L0/9chpqzUPPreZ7/NvjtnNS3joxFb+N/qE5ZoW+w/zdATHAVxb+vkIgMqUEMbYBmPsbYyxl4PnCA4BcJtAMiPWHB3BuVFxqq3DVv1ToFymKHDQst59O4wrPQQAMMxO+aaIYBImDVoJsEsWM8YqEQFgN90s713I1vWDDrwOGe9XdwTivTdRQyKZLBzWsq0jEIJ82f1Er4WNg3zoxGb+7y8+bSfN/WBJgfU+yzX3lKSxH3xuU31heW/l+zxlp2fVYv9hno7gTgA3EdGNRNQF8GYAHyxfQESr2e8A4PsAfJwxtiufVhH62zY3bdYGupThMgB+EiUYhUkjIlizTF6Pw+ap3vc66Podc44gajoRQOQIEq1ExThKEJdO9gA3ziaqJu9dyAwrEVlFIILrF/ezNczCgTRO9pYRgbiPcEA2GkzHTmxhue/j2vUBnrBULH3ouS0MAg8vOLyYi+uZcPzsGF6HcMuVS9ay5ULD6pYrl/CEJZ3UYv9hbo6AMRYDeCeADwP4EoDfYYzdR0TvIKJ3ZJe9CMB9RPRl8OqiH5rXfupYcyh/5KfhqNGIJeAyIeysojFtdYEbVVPymkcEklN91zPy/LocAcCpIxVkyqWLFpFE3cCKfxu7oGv36/oddL0Otoy9ErUcwdCO4hH3y2dBCzE+i4jgiTMj3HhwAdeuDfGUpSLr8bMjXLs+wJG1AZ4+b+c8njo7xpXLfTz/0CKesmg+BPgciCuX+7j5yqWZxrO22B+QH3F3CIyxOwDcUXvsPaV/fxLATfPcgwrCqNsMfd8OE6QMeUliHcv9wFqS+YyiMW110EUYpxhHifTULjAK5b8fdn1jRKCihsS0s9E0rhjsMooTevH7hZ5vpGrqlAsgchJuJ3TxHKZ1RUTA/74LllpKW9lr72QlvaJJzyYiOLk5xbXrQ6wOAuvKoRObXNL8mtUB7s6mqZnw1Lkxrl7t4/ByD3/9oL2Q3lWrfVy7NsSf3PMM0pTlr7FFC4F921nsUu1TpynqWBnYzQMAinr1RkRgUaXCGMscgexUbxERKJLFIscwidQD7DdqBhawSzJvKg26GzUE8Ndo7J4eR/A7hH7AP9pCS8mUP9maxBUhv6UsIti0iAhObnKRu6tW+ji5NbXrft6c4vBSD1evDnB2FFmVqXJHMMAVy31sTWOrJsBTWZf1oaUekpRZU6Et9hf2rSMYBB68DllFBBtjuTSDwPLAvmqo0CyqOpU1C0cQJrwCR9YLYBMRjBQRgXhs7EgN2RhmmeIpdyD6dZtTyf26Ng4kxvIgqDTNLfQ864igvKa8fxXiJMWZUYhDiz2sL3TBmDmhzRjjzmO559SQKJzHFcv2E/a4AF8v12w6PcPsjBZ7H/vWERBxNcydiAiW+3weQJyoT9QCKs0iUb9+bqz+oooTvzIi0Jx6k5QhjFNpjkA8pncEshP6jDmCvo8tC86+Loex2PONhnljEjUctu0+F0uvzfc6GASekYo6sx2CMT4zYT2TwzAZ9fPjCGGS4vBS39oRjMME0zjF6rCba1OZZLmTlGUyHV0cWLDvXm+x/7BvHQEAB0dQrUSpQ1Sm2CQWxRe+XopqQw2JZK6U3unqjd1Ys7aX0ShjjSMRUVE5T7KQRSG6BHfuCErv3ZJl1dBSv3qyX+ybo57taVyZ3gaI98YkvxHnlU3l+5n2KU7lhxZ7uU6R6dSdr1myXyMonbVh17r0+dwoRMr48CShZ+Uyg6PF/sH+dgS9wI4aMkUEDro758f8xOrXNItsvtwjySwCgWHX057oi8lmzbWDPEdgXr/QqyaLU6aPJLamMbpeBz2/cEA2J/utWs+CuJ/JMMtyKDbqqluTZqLchsISJ+xDS92CfjEYWxEVrg+7pYhAf1IvHEGQHxrOa6JHvjcxUa6XCxyebiOCFhLsa0ew2PfzJKgOG4aIQDwuTs068H6EpkMpRNXUX+5tySwCgb7vaQ25OO1LqSGLHEE+ByEoG3Sz/HU9CQuIk32iTapuSU72i12zA5HlQfi8BRtqSOIIDAeFQpoisDbq4vNUXWM63ReieEVznUE/qTRIaHUQgAhOw5ha7B/sa0ew7JgjWFZ0FrsocW5KOGwA6Aceen5HHxGE6oigH3S0VT+ygTYCeY5AQw2JOQh+ZQ5CJn+tOW2rDKz4ne5+UgdiONnLGu4WDLQZwP8usojAtt9huR+UJLbtnMfywMdyP4DXIYeIoIth10PgkbEHRtxnddBFp8Mb+WwOPi32H/a1I1jqW1JDY7mOv4AwdKbEIiBPZgqYNPB1yeJ+oKeGxO+kVUMWyeKxRPV0wcKgb06aJ3uxzuR46jSWmKSm64CWdU+bqpsY4+KDTSrKXHZabmALvA76Qcf4mSrLYHSyEZ6maFI4l7Uhz5usDMyd6PVKr+V+4CSF0mL/YJ87AvuIoJ64rDyPqDm3TDyrcg1LWfWRCiJZXD8pA9wRhHGqNJJjCbWTr+3a5AiaXcnCwOtO26Owyb3bOJ5JlGAQVD+eixZy21JRPkNp7TROkbKmkzQ5V4AfEngOhO+VHy7MlU382kLOwpT7OJdROqul+QymHMFGTZKbd8C3EUGLJva9I9ia6qtegKr8sgwLluqf4rlUuQaTY9Ili/uZcZ3GcnpInGx1EYHWEUia0RYsDPM4anL24ueRzqBLOqhzKQxDJNFMFnNqSPV3VjnJYdfTRi1AQfWJQ4LN4WJjzJ2joNlsihbOjSMMux66mcNZHZhHcIp9LOYRgZs4Yov9g33uCPj0KlP4X5ZflkF80WwGwOs0i0yctI7nF520qhNskV9org28DvwO6U/okiRswfXrk9T9moG1q1JSr1MZ5zRlmeOplY/2PKRM3Tmtos2GXV/rrACe+C3njmzol41JVDkMLFkULdSrqFYsHcFC18sn4bloYrXYX9jnjsBOh2bDEBH0fJ68M50EZVLOZSwbcha6ZLHJuE5yYye/9yDwMA71yebmHAT+XBONA+QUT/OkLZ5Tt65+PxOlNInljlI4LFXkIp6v7nj6gadNwAOc768bdXNEEFWchw2dtBXGjdJdcyVU9dCxbHGfFvsT+9wR2AnPqSp9yrDpXp1EaUPKuQyTds92mMDvUE4PlNE3OAJZ+WdlvakPQXLStkoySxxB33Cy55pKcXOdocxVFTEJx6nSYlKV1g67HsIk1XaMb0xkxtYmIig7At+4ZrshgWFTCVU9dPCkdBsRtGhinzuCrP7fokRQZbwFbCSZ60lC2X50JzZZaaSAoIZUJ1hdjgDgRlDfhxA3nIiJjhJ7VuUIVOvCRJ68zaMehQNRGfT8vYlN0ZJbBAIIbSPXiKC5xhSV1jumFy30k+qOYHkQYHMaW4nitdhf2NeOQIT0xrBcI88ssNjzjUqVpn6Exb6PkabRansqHywDAD2D0Zpmj6tKYDk15JaE7fvm/oNJlEq4d/06lUE3GeYiIlBELqr7RYr7GfYJCGqofrq30UNqrtFLdSQNamgS6aOVej5q2ZIK3ZxE+Ge/+Vn87cN2oy3/6ssn8M7f+pxVsUSaMvzEB76I3/r0E1bP3WJ3sK8dQV6FovkAM8a0J/Hyc5kjAr2KqfjSqrpnR1GSzyeuQxixqcoRxCl6fkdZAmuihmQn+06H91aoIok4SREmTaE78bMqR6DSRTJFEqqEuMmBCEOvTE5r3pd6X8hSP8A4ShBpDHT9YCGKFnT34dRQVaYD0BcobE6qzXxF1ZX+c/rBu5/GHfc+ix//oy9qrxP41394L/7knmfw4fueNV77uSfO4tc/+Th+7A/vtRJpbLE72N+OwGJoSZhwXl8m/VyGDTVkEq8Tomebisa0kUR2QSDPESjoD+EIVBgEaoMOIBuYI5enUK0znrQNJ3slVWOKJGrreob8ibpqSO+wkpRhEqWVz4ZwCvreiup7aTMbY3vaTBab7rNRK1UW9zQJ8N17nM9QfuL0SOvQAP6eP32eq6DaDNj5QumaxyxHdLaYP/a3I7CpS9d085bhQg2pG8r0BkE1TwAw5wimcbMcswxdjoAxJk36inXKk7aoxqntuet10KEZqCFD45tKi8lUUTVWrDMnp5sRiMl5xEmKME4r9JWN5EZzXoKNvEeVGhKHCFNE8MBzm3yvKcOThjnHYiZyeZ0OD58srn/I4voWu4O5OgIiej0RPUBEx4joXZLfrxDRHxPR3UR0HxG9bZ77qUN8afWVOtVZtiosWChcyoa7lLFo4HBHYYIFTbIXUBvXaZTmctPS9RpqaBKlYExeesq7b+XOZ5KVo9YNOhFh2PWNDkTJ9TtSQyYnqYpchtYRSLHPgcHYyqTETX+7KEkxjauRh6CJVJ+VPFrpNiMCU9/ME6dHuPmKJQBVQy/DY6f47593aAFPn9PPRwCAR09t4abDiwCA42ft5i63mD/m5giIyAPwy+BD6W8F8BYiurV22Q8CuJ8x9jIArwHw80TUxS6h53fgdcjY4QpAyc0LLPYCqzGKAJTNaaZyVj64fjZqaBInFSnoxnpf7Qh0zWh9TZJZZWDFOmWOIDew1Y9nUXaqMOgSw8x/NuQIFH0EpmSxrCR3aMp/SD5PeXmrYo04YFQdQZD9zj7PMrTIEUTZxLWXHlkBAJzY0IvhPXWOG/NXXr+OZzcmxi79Z89PcMtVy+gHHeNgnRa7h3lGBK8CcIwx9ghjLATwfgBvrF3DACwRz2AuAjgDYNc6XogIC129IFn+JTREBKKcTyeItjmJ4XVISTOJ0F9FDWnLR31BfyiooUifI+h31Q1lenmKDqYK51Nw7837cvkGdeTDn7v6nntZD4WxasixI3kSJiBqVlTlJ2iHvgXTqVu2xiS5UUx5K9YsGCKCYv5EaY1FjuD0Fp+4dutVywD4eEwdTm2FCDzCTVcsIoxTC2nsEAcXu7hyuY9nW0dwyWCejuAaAE+Wfj6ePVbGLwF4EYCnAdwL4IcYYw1rRERvJ6K7iOiukydP7ugmTdU+OlmH+vMAaqMB8C/tsOspK3dM5X3b2XoZejn9MVuyuOd3EBoMuipZbKJOZLkJXW5BVdefrzN0CMvE4wB1tCTyH/W/S75OGfE0ja04daveE/FZKzs5UzmtMNzViECfI5DlWWzKYYXhP7I2wMogMM5FPrM9xfpCF1etDAAAz22qjfskSrA5jXFwsYfDy/02IriEME9HILN29ePytwD4AoCrAbwcwC8R0XJjEWO3M8aOMsaOHjp0aEc3udDztZypLCyXYdGqWkRd9VN+DmVEECW5oamDl4bqHIE+Wdz1O0rBOlUyFbA06IphOC6nZpv7TSP5yV78rDbo8kT40MT3h00DnVfmGJzVQs8+itiSfAZNEuAy552rxWqooROZIT+01MPhpV7+swpntkOsL/SwthDkP+uuBfjozEOLPeOc5ha7h3k6guMAri39fAT85F/G2wD8AeM4BuBRALfMcU8NLHT1HZq603AZJloHyMoGNbmGQcAFwmQ5gjBOESVMmSwmIu2UMnNE4GEap1KOd6Q52eukmk0zENRyGLHyfjyprUhOK3oliIgP7lE6ulQZtfDXYaDMyqduQ45ASw2pqqFyaqjkCAx5hfxvJqWgzBHBoaUe1he6OGuYnHYqo3rEUJ7zGmrodGl05sowsBrk1GJ3ME9HcCeAm4joxiwB/GYAH6xd8wSA1wIAEV0B4GYAj8xxTw2YFCa38/JRQ9WQlTa/vjFNl7NQJULL0E0pm0T6ZLFwEqGkblzkAFSOwCT5MAyae9ZFBBON89Ulp6ea16jrnJ5IpLL5vXiUpaSiNDkCNV3WFA4stJD0fH/5PnmUYyqJDaprTMURYurZejYS06SkyiOCbj5qVTed7VRtdOa5UWRMLrfYHczNETDGYgDvBPBhAF8C8DuMsfuI6B1E9I7ssp8G8DVEdC+AjwL4V4wxu772HcKCYUB5/iU0VQ1ZUkNGh9KTO6aijFVT+aOjTWJ9+agwLDJ6SDwmiygGgac+aed9BDKRPA0VFXFxvcCT3U/d+DaN07xUVLpPTQQio4aIeOe0acaD1Kgb+H5Z+ajaMTbLcE1d3WPJ3njZrr44YmMcwe8QBoGHlYH51C4cgYgIdDO3BRW0PuSOI06ZlXR7i/lDb5UuEIyxOwDcUXvsPaV/Pw3gm+e5BxP4OEKLiEDDrwMlasgQERxY0FfHDrue9MshqwKpQ2fsTFVDovs21DoCt2SxLkfQ9z1ltdEkklM1+f00OQlVRGCisGR7BPhrNlFY5b+JiCJMfQTlNV5m1F0pNt17UdynWQml+7wLoToiwrLBEURJiq1pjLVh12rmdllra3WQOY5RaNTxGocJOh3550+G86MIK0O9SKRAnKSYxKlxD3sd+7qzGMjGGBoigp5fHdoug42Oi040rvw8MoogT0xq1nd9NTVkShZrIwKNYJ0wsLIQX1c11As6mGr2Oos4ni4P0tc4yXGUNrqf831aRQTFWiLCUNsjIS9H1hlose9+zRDq3otx7qRqU966+uKIsiDeyiDAKFTrJm3VGiRNM7e3SrkOYahN5aZxkuLv/Kf/iTf+0ie0pdkCH3vgBF7+0x+xFrX7ofd/AV//7r/U5jb2A/a9I1g0dARv1waCqGDTtWkjXqeKCGSUQh29wJNy/IBd+SggF63LIwIJ7ZKL3UmM5ShKEHhyikckp2UIZzToOvpLmz8Jm/OR831qKKxRKH99A42xFX/H5uwD9RplRKCLclQzFnpqJwVUZbVXMpVc1QyDrVoSe3XQ1eYINif8UNX1O1jNnttEPX3x6Q08cnIbX352Ew+d2NJeCwB/fPczYAz4g88dN157fhThT+99BmdHEf76oZ0tS7/csO8dgZA6UEk/y4ahy5BXZOia08LE6FQWFMlrUbOuKh8FgJ4n7wVgjFkni3U5AnlljTppqbsnP2nrDLp8na7MdRonjVNzvk8NjTKO1NGSjsIaK3ILumY5QUN1OtXKJj3FJs/R6JyiTMoC4Il73cFncxJhqVdEBIDaWAtHICKCxb7huafFfIRlS0fwhSfO5v++5/g57bUA8Pkn+fX3P7NhTEQ/eGIz//eXntkwPvdexr53BIsGSmfbUPsvYJIWAORSzo3n6flSZ5I3FWnWd/2OlOOPU4aUqWcRAAX/KncEemoIkL/uKEml09TEc6nKVXX5DF0kMdHoKelKa6MkRVdB/fU0kQSvAmt+NlRRHV8jbwocaqqoplGS5R4kzkOVL1F0S/c0ZbRAdWiO+L/JEQi5C1NzZnnAk43QHgA8eXacixSaBPDSlOH4mXH+Xp7c0jfDPfDsZr6X/S6AZ3QERHSEiH6EiD5ARHcS0ceJ6FeI6O8Q0WXvSEQ1kK4e21QxBAC+x0NelUMJY67NrzPkADf0soYflcRyGSo+W0ftlNcCcmpIGEKZsRTPKXNAUcw0BtYDY+pyVbUjUHdATzV6Srp5C3qHpY4IRgpp7qHmdD+ayg8DWr5fEbEMNDOVR6G8W3oQeMqZFUB1mM1yrn2lkL7IHhcVc0uGUatbkyh3ADazQABu/G84OMTVqwM8YXAEp7amCJMUX/P8AwCAxw0y10+eHaHrdfDKG9asBPP2MrSGnIh+DcB7AYQA3g3gLQD+GYC/APB6AH9DRN8w703OE6b6f674aVdRsKA51dn0AQAZVyyLCCxUUFURgfjia5PFgY4aSpRDbXSRRJSkCHy5nIaJilIaZg1nr4sk+r663yGMU2keQ+xTldQeKw4JWspGkScadj2MInWyWCXcp3Q4Ciel2xuQ5QgyB2Aqhd2s5Qi4Aq/+ucW1RQe2vnz0+NkxjqwNcWRtYFQrFQJ4t12/DoAL3OlwcnOKQ0s9XLky2PdyFyYL9/OMMdmYoi8C+IOsUey6nd/W7qHQdVd3dZpKPgV0CT/xJTdGBD0eETDGKobXRgW153cUJ2x1H0CxVkMNaQysOPFLHVCiN7DiudFv7ldZzukVlFLdMfE+ArUWkyqRHmoign7gaeZDxNJmuX7gKRuxprFazmIcyk+wY0U5bV/TUzEJ5VFEP1CXqSYpq/D4C3m0rI8Iiuv11NDWNMZ160MA/O/vd8hIDT1zfoyvvG4V58cR7ntaz+OLU/3LMuVUk2CecARXrfRxejvURpR7HdqIQOEEyr8PM2mIyxa5IqNGT8YmWQzoSwDzqh9DsnjY9ZGypkEehVy5VEW1AFkiVXJ6FcZC9yEXhlDVR6BK3haRhCRHEGu499zxyKqUNNSQ6HeQGHWenFY7LFUkESUMgaeOXHRlpzK6racp451EqfTvoIt0JhpqSKfEqtJqUu2tnvwVEYHq1L6VTdITh6nFno+tUD17uTw6k4iMjiNJGc6NIxxY7OHgYg+nDIZd6CK98Mol+B3KO5lVyCOCZX4SMUlu72VorVLW8Vv+qzIApwD8FYD/yBi77OOpYc9EDcVG450/l1ZILav6MTSmiVPY9jSufPm3s+ollXIpkFFDM0cEaoOu5+zVjWha7l1HRWmSvmVKqW5QdSWyKkObpAxJytD13PsIplGC/lKv8Xg/UOcVprG8cow3rukcgaR019BQphokZGqQK1M9gPq7sTWJQVR8phd6PhjL6FTJa9ycRJVZHIuGnML5cQTGgLVhgK5H2JzGSqcIFPIWa8MuDix2jRHBic0pbrt+DetZxH9mO8S1WcSy32CycH9X8tg6gO8G8J8AfP+O72iXITTeVaee7al6KlgdA01EYDvgpszLHig9PrbIVfR8eSLQLlmsqxrSUEMarj/UUkPZ/STGT2bki3W6yEVtJHhpbZNSEs1S6lyG2qir+h30mk8p1hckEYGmnHYSyctitfpJit6IXsCrrtKUNUpY63Mn+r7Hu6QVxnpzGmOx6+fPUx6fWXcEjDHJuE1PW26dS1IsdHM67fR2iGtWB9Lrz49CLPd9eB3CoaWeNiKIkhRntkMcXurnOknn9rEIntayMMYelzz8OIDPE9Hn57Ol3YVO8CtNWTa03TZZ7CuHbRTCYaY+AjlVta0oOyxDGREoOlPL0DaUKeiMyjopNaShXDSUkk0EUnc8vFdCFxEUlFL5tYj3S0Vh6TSRVDkJlUPmaxKpQ+4HnjIpPYlSHFpqSib0NUZ9HCU4uNjMbfVLkVi9eqnehNbp8C5pJTVUonoAXjUEcIrpcO3acZQgZU0pbZ0ktuhSXht2Mezy9+b01lTpCM6OIqxlp/u1ob65TTTJrS0EWB0Wchf7FRdS/nnZl44Cel13mXa8DrqmoEIiwtxHADST1zZlrF2vgyhhjVZ8q4jAVDWkWGuKCLpGBzJjBFIztMKgKxvRFEntKPtZGbloO5nlDktXq696bSLRL5NRGGuoIfGcdYSKyith5GWvSdbBPFSIIAL8OyObkaBTzy0fZkzUkIgI1oaFuqmuAe3sKMw7lpcNyqkbWaJ7uR8UEYGFzMT2NMYvfOSBvELJhM88egbv/ZtHrVRW05ThVz52DPceP2/13DsJU47gFZKH1wD8YwAfn8uOdhk6aYhty1O8wEJXfcIZWSaLF/L9VJ/HRrk0r+lPUvQ7xRfOKlnsGTh7TV2/ah1v1FJTLtr7KUXg5Pcz5UHKjm6pskf+BbVpfKvnZ5SJX9+T0lBijUrOW+yvflLXJYuBrFu5tiZSOGHdtDaZdpJulCsvrS5dqxmfKZtfvdD1tWWbZ0VEsBBge8r/PhtjfU5BqKCuDAKlNAZ/HiGA5+fO46xFRPDrn3wMv/iXx/DQiS386j++TXstYwzf+747sTmN8ZJrVvCqG9e113/0yyfw7z/0AP7b8uP41I+91riXnYSxfLT2MwNwGsDHANw+jw3tNoKsEUxmwEdT94hAqTFjmSzOKzUkEcG6oYy1bMzLhsMmWex7XKteRdWsKe6tSxab6vOB5smeMWaoGlI4gmi2iCA0RQR+B4xxh9H16+WqKppH/nfI1ygiAvH7WRxBHdM4ldJyYm+yyLWghqqy2qqIoN7oVgzLaV6fC+dVHIde8LGc/PU6/N+6U/7ZUYgbDy4A4Cf9jXGzDFtAPM9yP4DvdbDU960igk8+fBoA8OlHzyifW+DRU9t5r8Unjp0yOoJPPcKf+9mNCU5tTXFwsVmIMC+YcgTfqPpdNkjmuR3f0UXAQleetHKOCDJBL9kHxDZZrKrd3p7GuHZNX9GgkpLWyUhX1qsa0mZMFkeaZLEwSHUKxSSHoXI8RdTj5kAEpaTKZRQn9aQSNSQpQ5QwebLYL+iXpiOQRxH5SV1a/ivvq+jlRr1peHkuREcNNe8jm8G8qDHWk6h6OBHrpM8dNmcqiJ4ZFc6OQgQen6EgiBXdKf/cqBoRhEnK3zsJHSsiC6F5tDoMrCIC0ctwZjvEme0QBzTG+v6SftFDJV0jFe59qqCEHnpua1cdgRPPT0QrRPQ9RPQXAD43pz3tOoYKSsdG+rn+PEnKpAnbURjDN/QBiOcAmlVMNsqlPa84VZZRTBjT31snUTFTsjhhWukGQML1Z/fXUTWy+5miHlEe2sgRJAZKKZ8EJo8kVBLb5T0JpClDqBieo3sfVTkCnfNQOeGeIzU01MzrqFcz6SIUGTWkq3oCeDJ6qR/kk/s6pJa7iJMUm5M4F8oTOkmqCKIcEQDAgQXzDOXtaYwz2yG+KjvZHzOooYpO6FfduI5HTm5rrwWAp86O8ernHcie2+w4dhI2WkMDIvouIvoAeEfxLwD4GVTnEV/WUJWxifppk1CcQP5FkA6W4SG/LpQUewGaVUzbFo5A1RSW0ybGiEBeuTLVNGrpyzktqCHH6KVIFtfX6fMgKkObl486Nr5phfhKEUEZeULbISKIkhRJyqQRQf73lhw8QkUzX743HTVUG3ivqhqq5yaEA9Qnoos9DbpF1ZMMvAyVP78YlKMy7CIvIU74JuXUco4A4DOaTX0HwrB/3QsOVn5W4ckzI6wOA7zk6hWjTlKUpHjm/Bi3Xb+Gq1f6+Jtjuzqo0ag19JsAHgSfIvZLAG4AcJYx9jHGmFrC8DKDMSKwzBEsaHoSRlM7zSJRu11/DpvGNpNxVZ2y8/WBvJZdp/NPRMqu3VmSxToDq18nTuh6Css5R6A43U80ztWYxzDkCKr3UetE9RROERA6T+pKI1lEIJth0A88ZR/BOKzmQLQVSZIhRQNNdAIAW7XvzFLfV1JD27WqPHHSV12/MSlGcgLAwUUbR8CN+W3XrwHgDWk6PHl2jGvXhrhiuYdRmGgrpJ49P0HKgGvXB/iGFx7C3z582moQz07BFBG8BMBZ8JnDX2aMJah2GmtBRK8nogeI6BgRvUvy+/+DiL6Q/fdFIkqISJ9RmQMWFMNpxGO21JDo5JRxtjZ9AEBRu13+8oVxiihhxtJTk7EzOoIZqCGxbvbO4tpJO0/6ulFDpsoolZMMbSMCZQRiHxFMYrVRL1cNlTHWOoKiN6KMNMtfSCOC/NQuoS8lg4QGXXUpbF0Mr6+JiMV7UZ3vrK7YA/j3b6nUp7DcD/KyT9m1QFHCao4IipGcAI8IzoxCxAo9KgB4OisZfcHhRSz19BVPAPDs+TGuWunjilzCQn29KEe9enWArziyis1JbF2iuhPQWgbG2MsA/EMAywD+goj+J4AlIrrS9MRE5AH4ZQBvAHArgLcQ0a215/8PjLGXM8ZeDuBHAfw1Y+zMTK/kAqAa3zfOP7yWEYHmg61SqpRh2KuG47bKpSqqIEpSdIjPxjWt16mP6tfJcgSWonOVe+mpIZXstalXQuUkjeWjeVJbkZNw4O61EYGCVhFr9BFBje5K1Y5f5aQA/jmr30c3x6EufRF4XEjONkegcxxAs09huR/kc49l1wJoyFyraK2NSZTTSABwYKELxvTdxWe2RRNaF4eXe7m2ke76A4s9HM5kSJ7TaBmJ/MShpR5uvnIRAPDgLs5IMB51GWNfBvBvAfxbIjoKLkX9GSI6zhj7Gs3SVwE4xhh7BACI6P0A3gjgfsX1bwHw2y6b3yl88/YfI906BuAbKo/3Tn8J7/Zvx8pH/gTwPIAIoA6A7P/i5+yxF21M8S/9MxhvvQTAauW5bjv3Z3jp5HPAH/xGsaa8vvTfv2JP4fi5bwHwUgDAaDrFv/H/O/7X+wA8u1DdfCnn8KKtEP/O30Cy9UJwJRCO5a2H8a+D3wQ+9DcAYwBLAWT/L/38Q6NncTdeDeCV+dokSfFD9H687tgY+L3lYt+g/P8/xZ7Fk+femO8XANI4wo/R+/DqYwB+b6m0T8q2Tfi/u8/g/LnvAfDCfN3M1JAwmMqIQM7159GSq8OyooYUeQVdstgh8lBGOZrXJKghqbGW5KGENpFMmiOW5C5UonYy2km3F4Cf8q87UFTKLfZ95XCaUW2UqyrXJrAxjirRRllvSFWtczaTsAi8Dg4v9bWGPU0Zzo5CrC8EOCwiAo3jyOU0hl1ctcw7px85uY3Xvki5ZEdhx3lkYIzdBeAuIvoR1K1mE9cAeLL083EAXyW7kIiG4PMN3qn4/dsBvB0ArrvuOpctW+HG8EFcE3+m8fjzn/5jHPU/BvbYNWoDmj/GcEUS4V/4Y9xz/HXALddXnuvvbf4PHExPAU9cUV0r+e9N8QbuPXUawFsBAOGpR/G9/p9hcuYgMF0pPWuVpVueTvAW/2nc8+yngZfclD/+smd/H1/V+VPgc0uZIUfNoXHj/prxGRyJnwDww/na6eg83ul/AOPzq0C6Xnsf+B5enz6Ne08nAP5pvi4+/TDe5n8Y2xsHgHQ1f4/K+35j53HceeIwgDcV97NNFjsY2eo6edWQTmtIdz9pNY8vjwh0eQVVo5cuv6MqpRVRjrSPQLE3QMwwqJqEftBBKumjUM1RVg0AyimuUrJYGG1dRLBY2o9u3sdWjRoyKadu1/IPQmr+9FYIXCFdgtPbYe4w1he72vGWm5MYScqwvtDLpT50VUlns9+tDrsIPMJiz8fT53ePGjJ1Fv84gF+p0zWM90v/NRF9E4AhY+xPZMslj6nyC98G4BMqWogxdjuyBrajR4/ueAaFBUP0meRND7dxmq3gwL9UBTFVPPnA3bjht78B8bhZVtZLx/j80mvw6h/+H8bneeBnvxbduCg3m27z+uIHbvtJvOyb/4n6/vffhef/zmuRRtWTSied4iTWcOjHHtPe9/53vw7D6GzlsXAywRDA3c/73/DVb/kx6brH/s+vRCepfsjjcIougL994bvwuu98u3Td2Z+8Dkirobh1jkBxQjee7BWOQKc1JFunc1imiEBbPtqoGmLK/amcW5ETUu/NmhoqOaiyM5pIkr+AiAjkOYIOVV+HbswpkDmC0ql9oGlu265RQ3n+QaUqHMU4VDr5ry3YGWtx3eog0DagnRkJwbwgK4HVS1icGYVY6vn5e3z1aj/PSewGTBHBvQD+mIgm4H0DJ8HHiNwE4OXgk8r+L8Xa46iWmB4B8LTi2jfjItFCAMC6QwwwRZKyCo/eiUcYU1+zsorukHN76bRZM9xjEyS+ncRt6C9gMSzKx6Ix5wr94bJ2nd/le607AkpCxBbBX9rpwkurX4Qw5OEs+ermlqQTNNZFU/O6CAEoru7VRA35HUKH1ElfkwNx7yyWU1G6BradjAhyR+VQaWRaQ6Qq8WwWNPRKlUBlCemJInfRDzrq8una6ExxapftJU1ZQ85aFxEIByGuLxQD1PIY5So8ERGc0TSVndkOcfUq/46tDgOcG4VSwT9+Lf9c865ownI/0OoknRtFWF0o3t+rVwe7Oj7TlCz+AGPsawG8A8B9ADwAGwB+A8CrGGP/O2PspGL5nQBuIqIbs0lmbwbwwfpFRLQC4H8B8IHZX8aFgYIF9CjGaFz1wF48wrRj7wiGmaFOw6YjGGACZukIYn8BvbTYSzzmIWh3sKRdFwR8r0lcjwgiRNRUr6wj7XThs+qHNZryfXQCtUGPKUCnti6JMkdgWlePCAzUEBFxLZ9EcQp2NOihxmjydYokrjZZrFpjwfcrGtfkzWH6xLmMGiIiZQJYNsymr9iXLPkrflZRQ3UaSdeAtp3PRig3t/GiDllZ5ZZEDmZBIws/DpOK3EseEWxpIoJRmHcurw27SBmwpXh+kVg+sMA//8JxqHBmO8T6sOjSvmplsKtVQ1Y5AsbYQwAecnlixlhMRO8E8GFwB/Jexth9RPSO7PfvyS59E4CPMMbMrXdzAvX4SX68vYmlhcJY+/EYYUcueStDf4EbalZ3BEmEAAnQtXMESbCIISuSYlFGNXUNEUHQ5x+6NKp+4DpphNjKEQTwUf1gRyF3KjpHkHYCdNPq6SUW6zQRQUy+xBHwL7Su1LUXdJqVMoYyUGXV0Ix9BLoGNnOvg31EoJPAUAkFmrqlVfMSxmHSSJT2SxFB5VpJg5i4XkoNSWgnreBjZtgXe0Hj+knczGVsT2N0qOqYhl21PMb2tBr9BF4Hiz0f58ZyY80Y48Z6oZCwAIDzo+qwHYE8IlgoGtx0FUk8sVw4giuWeadznKTwDWoEOwGnZLErGGN3ALij9th7aj+/D8D75rkPEzo9Xokz3j6PcqYoSEaYODiCXq+HkHlAbfZsNNlCAACBnSNIgwUMSjkLNuHUUG9B7wi6maNhtYjAS0IrR8C8LoLayT7OTvY6g552uvDiauJMRASeNiLoolOjlHQllgKyBjYdHVJ+vEmj6MtH8/p+ZVmnfSmojk7KjXo9R6BJFquEAk10V19xapeJ5KlKYVWNboOulyc+yxhHzVnNfU1EUIzBrJ7wATGtr+YIwhgLXb9CPS1o5DHGkgluy31fKWExjhJM4zSfXVDMMIhwraTz6VxJMA/IHIEuR7Ad4gWHFvOfhYbRmREfnjNvzN/VXAbw+/wkP9muJnmDdILYs3cERIQJ+qCo6gjG29yQi8jDiN4ShjRFEvMPZRryfQ0WVnSr0O1mRjeufhGJRUgsqaEA1Q+ryDfoDHrS6cJj1S+QcAS6SCKhZm7BaqymZFiMMH6+olfC6xD8DjVzBAbROdWpW5ssVtT36+ikToekfRymwTldr9nMZ3KKqoRuKFEs7Sv6KFSNbn3fgRrKq4aaxncrjwiqyWJArm66PY0bfTpDhTxGlPAGzToNtqTpUxAOYqUkUgdAGUFsZRGKuMfqsKvNEZwfV/saDparmHYBrSMA4Pd5RDAdVxs4eukYsSWvLzCmHjpx1RFMRvy07PXsnot63DFtb50DALApdwT9RUNE0OMnB5ZUIwI/Da0cAbwugoZB58+lS/qyTtDILQh6yvPVp5mk08wt2HRBy4xfmHXS6rScZB3Q+em5oza05esEdHw/UWbUa3kM06Q4PvTejfLSOUX9sJ0mNSRTLFVRQxNJXwCgnqMsq0gqtLkkpay1clBA37C5PW3OSeYaYmrpmKYjUEcE4nERoRQzDNSOY7FXRCi8ykhNO23VuqhFRHBJOQIieiERfZSIvpj9/BVZaemeQDDkhjeuOwI2QerbRwQAMKU+OnE1yTMdcUPuWUYEnSxCGWeOAOE2Quaj29WHiJ7vI2adRkTgsQhJx5IaqkcEoflkzyml6j2FA+l01TMUEgrg13IEpuQtwA2crHvaJKEhO3Hz7meSVn4A/KQeeNS438RAYXW9DqK4mtScaCIC/lzN+cjiOWS6QeL+DSdleA95jsCuE1xVAaVKFvcDT2rY63IUAI/Sun4Ho8hukI2o8pFGBGF1HjKgjghk4npApmVk6FwWxnopywtsaRzHUil3sDLgVUOySWVcur4a/Yh8weltvZ7RTsE2Ivgv4BIQEQAwxu4BrwLaE+j1uYGOavX/A4yRBguyJUqEnQG8WkQQjrLyz76dI/AEVbXJ+wco3MKI+kblUgAIEQC1iIA7Av1QGwBgXg9dSsDS4ssjKpBMjqCeZLanlGpVSon+hA5w4xdJqoZU9E6xzpPSKKqTs0DgdXKuXmAaJ5xu0iSnw0Sho6SJCBoUlIkakjk3QwWVjhqqr1FVQAnHIKsEkvcRyOcjDLueVAlVRqPpk8vN0ldV1VBealrLEXBqSDGWcyKqmLhxF/0NIpfRuH4aVQz7Qs9HyuSNfLk8RikiEE1ol1REAN40Vm+9VUvpXWboZknYeFqKCBhDn4XWCV6BqNOHn1QjgjBzMEHfzqn4Q54LmGyfAwB0om2MyS4yicgHkuqH008jpBbUEPn8migsPnwsM+h+oI5GWKeZZE7jiXFdKqGUoiTlvQIaXaTAazoCG4Mu00TSTVErr6tHBLrxnYCcvprGXNRNpfnU9Tt58lrAZNRlUYQpquoHnlTxU6ZYaqoaqtNcKgVbVcQ2COS9AbKqrGEpWVzHWNIVPVQM1RlpIgKlI8gMft6wlr0vW4qqpHoz3KJmjGddJwnIJqd16JKLCE4R0fORdQYT0XcCeGZuu9plDDJqKJkUZZ/hZBsdYkDXMsGbIfIGCOqllBNR/qnvAxAQZaKhyC1EW5iSZQ8CgkYljm8ZEcDjp/cwLBxZmtFMXqBezyml6gecZev8rr7aqOkIGHzDyT6Q5gjM1FBPYtDDhJkdgdSoy2cP52sklM0kks8IqN5HkSNQSmCok8W66XD1k6lKsbQQ3VM01NVorq7HnVm91l/191H1Hciqx4ShH0uopGnUHPjD5y275Ah4slhG3wgHIaihTofLQKiooa1J3IgIAIUjmDQdQadDWF/oXnIRwQ8C+M8AbiGip8DFaH5gXpvabfSzahyRlAWA8VZGy1gmeAUSb9Csqc+etzewcyq9bD95I1m8hYlnF01EFIBqcg8+IjDPJiLoZvst9i9KUT1dfkLqCLJIQrNO5ghsTuiB30FYOzmrBrGU0fU7EgkH/clerJOd7nXrZHmFKEmVM5XFfRoRwQzUkCnh3pdMBlMplirLZxXVXXmTm4S6U4ngyagkWc/FgiYikEmlD7LIp27YR/kI2mZEECVMqsArO7Wr5OsBYLMREYhZ5HYRAcArk3SVRjsJ24ayRwD8r0S0AKDDGNs9fdRdQDcz0KxU/z8ebWIFQKfrliNI/CG6rOoIhOSEbUQgHFOURQT9ZAuj7gGrtTGaTVoBi5DaRASZI4jCpiMINBQPvC4CSpAkCbxsJGTuCHSUksSBxKmFQfeowdnbJotlBsqUW5BVAKk472KNhzBuGnXdvQJPUt4apyCNhHhP4txCjegcX+NZK5Yqp61lhr2etyqX25aN+OzUUCki0CSLZY65F3hgrLmXQgm12UcAcInqerRXn3cAcMOtGjazNYnz5xPXAvKIQEQbZcfB97N7jsC2auj/IqJVxtg2Y2yTiNaI6GfmvbldQzBACgLCIiKYbnMj7BtkHepIgwH6dUeQUU4DQ0OYwHBxNVvH/W0/3UYY2K3lsg31iCAG88yOQDSNiW5iwI7iYWLdtNQEF0eIWQe9rjoSYV7QqFKKYguqRpIsjiwonplzC5JksSlykTodw71kkYegrlSFAl3fazgp8Rw9T52UblJQ8sa6wOPaTrKZzVIto3zAjl0uRlVuKqvKEnIXsqokHhGoopnq9TpqCCgG25exOY3R8zuV16x1BFMFNSShkuqCeQKi0mg3YEsNvYExdk78wBg7C+Bb57KjiwEiTNADRYUhm2YJXt+2CSwD84fos2qCR2gPDRftnMpwaZU/V5a8XmDbSAK7fcj0ewJEYBblo5TV/MfliCAxUzzInMy0RCkhmSJEoDV8siQzT1iacwSzVA0J/rp+v1kiCdO6npTvl08NE1CVxfY0a3hE4JZXUEVGYg9lEJFUNkJFjYm9yqSxpY5AMcBeVpXlZ4NvZInuSZQ0aDdVM9xIDJyqNaAtaiKOOucPcOMucwRxkmIUJhV5DHHal43EVVFDl6Ij8IgoPxIS0QCA+oh4GWJMA3RK0s9RVvIZWPL6OboLGNIUcVz8wSkaYcp8dDUJ18pT9AdcqmK6BTCGRTZC2rWLCBrduoyhy2wjAn6N6AoGAAiKp6d2BHluIayuC+FrjTPzew1qyHRqBmZPFgceNR3IjMli0z4DX0bzJNo1qoY3VQ9BvsZRgM+1G7kfyKkkeUTQ7MRmjCn/PkNFRKCqypI5JcaYNCJQzVAe5zmCet+Bujy1XgUEcMMtH3GbdUVbUkOy8lEAWL4EHcFvAPgoEX0vEX0PgD8H8Ovz29buI6Q+vFIjWORY6SNAmd7PaFTqSYhGmFDPqg9AYJuGoGgT08kWAkqAviU11AmqtflpwqufbBxBVvOflMpHRSmqLkcguo7jsgNJIkTwtUYMnQA+pWBJ8eXgA+8tHIEkWTyTA4kTqySzzHCaup+b0QfTRjuqyEPnTFUJ8A5B2+OQMlTm8+oUS1WOUPb6ZZ3YxUyF5nMPFNLSMsMOyCueooSBsWbiWqWTNFJ0RetyELKIYLEnLzfdzEpNlyypoc1JjK7XaSS7VwZBPuBm3rByBIyxfw/gZwG8CMCLAfx09tiewdQbVOr/EyH05ugIRHJ5sl3k0zvRCBO4CUdNaAAv2sbWeT4ohvp6nSGBlAJ4JWqIZfX8ojRUh46QsS5RZEhCRMxDoBleT9lzJ9OCEqNkiikCvZH1m+WqNuWjsoYyKweiyi0YqCi5gdbTPCrnoY0iZJSXiYJSJH5NkQdQPbXneQWZcVdEHbLXL4sITNGGvKEsUSq7NiuY5EqwghqqRxzjiDv/egJe16ewOZU4gr6vp3pKJ/xh4IFIXjW0PY0r4noCQtdIpX+0k7BWH2WM/RmAP5vjXi4q4s4AQVJUDcVZgrdvkH6uI1cyHZUcQTLG1GHADQBsd5YQhBuYbJ7hzzGwcwRJpwsvPp//HIZTzuH55ojAy6mhkkHPKJ6hpsGLMsorKq9LQoTM1xokypPTIXpZv5xdp2+T4rHi+hU5guW+/msgiyRM67qKmQnmPgJ3o17PRZicR1mSO5PL1/YeSJPYCmqo6zXHZ+q0j5TUUNzsCwCyiEAxzrPe06CqeJpGqVTmQzc6c2sS50NpBBayPoL6POetSbPCqNMhLHR9aQOajHYCSlLX4yhXO50XbKuGvp2IHiKi80S0QUSbRLQx153tMiJ/WKn/Z1kFkZgxYAsvk5EIS9SQH48QOgy4AYCxv4JBdA6jDR4RdBdWrdalVO3WzXl7C2rIyyqDKo4gjRAi0NJaotqonFuglE9FU5U98htmncylaqPZKZ7ZHMisncXGqiGF8zAZaOkaU7JYZqQN0QpQrfXXdSPLZL9V1I1sWI4uIhgEHuKUNfs0omZfACByBM29ALKeBkXpayJ/7mLOsVzLqJ5TWOz5iNNm34Ggnhbqkhc9TypJwQXqmsUcZUcwb9jmCP49gL/HGFthjC0zxpYYY25H5UsciTeoTAVjUx4dDC1LPgWEjMS0FBH4ycRpwA0AhN1VLCTnMcoigsHyutW6xKs5AouRkQJe0Bx1aTPmUuQW0rIjSMxT0fLcQinJHKd6ygXgjiBOq92rs5aPhhJZhTp6CqOuW9f1ZQ1l+j3Kq4aY0XlM47TSNGWTvwBqxlpzaldSQ7q5CqVT+1Tz3IO8W7imZRQniulvkgqmSE4NDRQ5AlUiWpcsHkumt6nUUFUS3QtKyYumThIArAwvPUfwHGPsS3PdyUVGGgzRK9f/R1sIma+VVpAhyPoORLIZAIJ07DTXAADi3hqW2QbC8ycAAAurVxhWcLBOVQAuzmWkza9D9AqkcTkimBoNepFbKNZ10nAmR2BKjgLFyVJ0wwK2VUNyg64rzxT3kyZLZ4gIdK9NWTVkuA+ASkLRpscBqBlrXUQgoZ+mSYqu5FQtmwSn644upKgl9I0iWayeDVFPFmflo7alr34HHZIni8eRREZbUEmKAUR1LSMVDcZ1kprvpZh8thuOwDZHcBcR/Q8AfwQg/7Yzxv5gHpu6GGDBAgaY5APsKRphRH24MnNCRqLsCLrpBFuOctbp8ACWsY1kg0s6LR24ympdvTY/shg+LyAigvKEs04SIYLeoOe5hbDsCMzjMTsit1BaZ0XVZL+PEgZBw3I6xOxAmqJudg1sjZN6zIzVPK4lp+I+Zc45NHUj+8V7IeyyTYIZqCV0NSWnPb/TSHIqk8WyRLSOGurKE7pTyThKgPP+9UlfyhyBghqaKqIZIs7jyyKCiWSwTj5hreY4VA1rw8CXOxnJiFDg0qSGlgGMAHwzgG/L/vu7pkVE9HoieoCIjhHRuxTXvIaIvkBE9xHRX9tufMfRHWKIaf6B7EQj5wQvfxruCJKSI+ilYySOEUFngUtKeGeOIWIellbsJCZSTx4R6EZNCgRdQQ0V5aOUhojJRA1l6+JqRBAbHEjHF5FELSKw6AcAUOn2nTXJHNo2sDlWALlU2pTvw1j1dG9TaSReR3EfvXPrSU7tphxB87XIT9XSiEDMVJBGBPKSTXX5qJoaqiuh6hyBSvNp0PUwqtE3YqJZo9xU0Fq1Tuexojy13/UwlshQTxQRwW46Alutobe5PjEReQB+GcDrABwHcCcRfZAxdn/pmlUAvwLg9YyxJ4josOt9dgoU8EawE5MQiz0fXjzCxFL6uQxRZST0hQCghylSRzlrf5Eb/qWtR3GelnDQYOQE6rINicXw+fyeGTVUnnDmWZzsfYkj8NIQSUev05QnmSvUEEOgSzCjOAWXjZNN1VA5tyBkrqNZk8WG+wVeB0nK8ghT7NGGsgmTYmB5GKfoLeiooeK5bfemo29sq4ZUp+oi2iiMqT4iUBtrmbprT9JHoI4IMmqokYhWCwYu9Py881ggp3oUE9bqTmysoIYGQQfPnZfLYtefW+y/63WU0tg7CStHQER9AN8L3kOQH5MZY9+jWfYqAMcywToQ0fsBvBHA/aVr/hGAP2CMPZE93wmn3e8gKKv2WfmN1wEe4dbRw3jKu8b5eQZZldHLjv0q8Ku/DzCGFbYF5ugIesvcJ94S3Y9HOjfgoOU61uliiCnwa98KMIYrtrLyUwtHEHS547v14f8K3P5hwAtw4/g+PEpHtOu8bArZCz7/74BHfw3wujgUHsfDnt6vi6E1ScmB2A6KAQpDJiSUbddFaYpeJytzdCg7rTgQU44gp2xSeNm9jA1lYn8xg+AkbSOCsiOIDDSZtGpIo1iqqtLSRQSyHgUZxVWUbNaNu9xY8y5nlRJqc34yIK8aqvcElPdTH2+ZJ39VM5cbncsJOtSk2YZdXzmNTeb0iAjLAx8blxA19N8BXAngWwD8NYAjAEwKpNcAeLL08/HssTJeCGCNiD5GRJ8lorfKnoiI3k5EdxHRXSdPnrTcshvOX/Ma/HnyCoT9w8DyEdwdvAx/sfz3nZ+nu3gAv558C55ZuBlYvR7p2o24I30VHr/yW5yeZ+HGo/hwwv/7i/W3WK97Yv1r8cn0VgAEdDyEvQP4UPJKTNZfZFwbDJfx/8VvwKmlW4DhOhAM8GT3efhY9xu169jajfjN+LU4v3wzMDwA+D083r0Jn+x9rXZdsnglPpi8GmOv6JGwoYZ6ftX4qSSU6yjnFoCS9IFFRADUIxC945ENvTdFH+J1T0uTzWwqjfhzl5LFllVD5Y5knWKpMlmuSLgC8oYymWFXnapVtf59X1Y+2lQqBXjtftfvNPMPitJUIHMEtRzBJHNSqoignugeZ2M56yXX2jGeEmoI4AnjjUslIgDwAsbYPyCiNzLGfp2IfgvAhw1rZEeSeq+0D+A2AK8FMADwSSL6FGPswcoixm4HcDsAHD16dD791odehO+PfgS/802vxqtuXMdP/N8fxw1rbqd4AKBOBz/vfS8evfEIfvLvvRibowjvvPsj+LcHbnV6nmuuvAqvjf8lUgb84yPXWa87sfYK/Gj443j4u78VXodw/yOn8Y7bP4XfGprLTwO/g5+J/wmmL74ZP/iNLwAA/Ox7P4Pz4wg/qFnX7fbwr+PvxX985cvwnbcdydedM5xk0gMvwL+I/jl+ffmF+WOmjl2gfApmlf+b15VyCz3OxTOmHuAi0Cs5gn7g5ZTPLDy8SaiusUZBweSvKd9b2XmklUHojftIZgboFEtd8h2y0lRdQ5mu8kbeR9Cct6wbAdqXSHDoZkkMu35jyLxoYGvmCOR7H4VJXhZbv76eWFblHwSWBsElFRGInZwjopcAWAFwg2HNcQDXln4+AuBpyTUfyuStTwH4OICXWe5pRyGUCEUzyXYYN2aaWj9Xt6gOEM8nayHXoVuSvL3hgP1MBEE95KflXInS/KdWnWJN1ThSeiKxqOKRdKHalI/WqSEd9VBZV48kFNLLqvuJhLFJ3bOypvGe6Iy6+NtVT/f6yCOrLoqr5aOm0lZxXXlvgH1DmcpBdTrUuF6Xf1CWj2qSxXHKpDpJ1n0HiucG5BFBMey+ukbnxOrXAsU0tnLPhyr/ILDc97GxCxITto7gdiJaA/DjAD4IzvO/27DmTgA3EdGNRNQFH3b/wdo1HwDw9UTkE9EQwFcBuCj9CsLoi4qBUZg0ZGptMewVHybhEGQnBBNed+uVAICvutGuYggofclrjsB0WgY4J1mvrIlTu2qc8r0A20ljzXW21T9A8zWaq42q702oMVBl1Kkhm/e0npBNUobUEH2oHKPOococjqlJTpYs1jnTen9DkjLe+Ke4R51K0iaLg6YxVamJAvIEsIoa4tcrHIHEaQDiEGfXIDZQSFKMwwTDoPl9H3Q9pKwaianyDwLLuxQR2Fqnj2YzCD4O4HkAQEQ36hYwxmIieic4heQBeC9j7D4iekf2+/cwxr5ERB8CcA+AFMD/xxj74oyv5YJQCE7F+f9njQjKtciik7Debm6Dn37ji/FPv+Z6vPSInc4Q0Dy9hprSPRnqQ1jChGEoOd1U7qmoQllQJOTK9xLXAtwA2CR9u7WTvUl2uXk/VltnnmNQvl5HddT32HBWDs5D/NtUnVR+frFG1ySXJ3RrDlilWFrvbygE6jSDb0pUle7vM5B054r3TFbiWS4JFQnfiY4aklUZKWgngNuBusSE0hEooplRlEgNe/l6cX9V/kHgUssR/D6AV9Qe+z1wfl8JxtgdAO6oPfae2s//AcB/sNzH3CCM1nYYI05STONU2tBig2FpaPa2QvvcBqvDLm673sztl1Hnz4vTmJ0EdlDjg63KK2v35PfVSyOU9ypCe1uqpu5AdCdO2f0atJltRFDbp02yOHce+b30oyrL14o9ujgcsUafixCTu6rGWnWfbtbfEKe8ic40E7kuja37u/b8DoiqlT2qTuHyY9XrE+Xz94NmNy/XGlJTQ3XDPlH0BQReB4FHzXLTMMFAEnGUqaTV7DGxN1kfAQAsD/yL30dARLeAl4yuENG3l361DDjqKl/iEBz+KEyKU/yM1NBCz8fJTV4SKT5Usz6XK+o0jU5DRoZuTY+HGxX3HIGNzn+j+if7v2/qI5jRoOeGtsb120woA8oOy8Ko19aYjGf5d+Jamwip7vjFem0iW5IsVvUF1PcVeJ28qklLDdU+C+W9lkFEGNbmFufJX21EUN174JFU4LBODeloJ3HPaZxWSoXHGh6/L5mwNo4SHFpqlmvL1E11zw3wiCCMU0wUJaY7BdMx9WbwDuJV8G5igU0A3z+nPV0UdD0+Bm97GufJGaH14YpBKbzcVrSbzwtKPtsiWQzwL2s90WcysF6Hz7WtOhC9BIO4F9A0zK59BKJ00mYeQfk+1hHBDBFIbwZn1av97WzuU/QeVN97Z9E5TSK7LltdVBgpjKlDRAA05xYLwy2PCJq9Abpy0J7fqUwFE0NslHvJZy6njWSwrMRTFkGMwhgDSd+QMORlpydyiCojv5x1F29MoovnCBhjHwDwASJ6NWPsk3PbxSUAIsorBkQnn64ET4eF0odDNKfMSjO5om60bI1dvr6mx2PD2Yvnd+70VVTxmJK+IkIJc/orsVtXOz3rVDGr95M7LJdksU5mQaBZnWR2cLKEuyki6HQIfoescxH1w8XUghqSSVyrDgb14TQ6aqifD76pUkOqE35d/E81xKb+/OXafvFdlhnigYR6mkSp0mmI5y6uVTsZAPnMi41xjMNuivhOsLMOwJuIaJmIAiL6KBGdIqJ/PL9tXRwsZDNI84hgMFtEMOz6pRyBSBbvjiOoNxjZGjuBZrLYXM5ZrLPXyBFrAJmBtSs7bSbEXSklcVK1UzvNoywLo153BFan+8Z9bPIK1TW6+cD1ezXLdtXGtPJaNIaaP17t/s2T6x01L1+hhjTGWkUNqameaoSrKzUF5FVMuhLPgaLKSHqtJCIYm5LFpYhgnrB1BN/MGNsAp4mOg3cE/x9z29VFwo5FBFn5KGMsjwhUHn+nUdfhsW22KtaTU7JSgEcSdd7WkrPPG8MsqSFFr4TKGDTvV3c8+r+NMvHrYtQtnJzqPjK55/qaZnOd2bnVFUJNEYFtvoOXINd7ISjn3OvoB15l6pipL4BfU71eJSLXjAj0nxWZUN04SuB1SOqQB5IGN9V8AVnfgU2OAMDcS0htHYE4Gn8rgN9mjJ2Z034uKhZ6fAbp5gXmCIZdPrkoTFJshzxpasvRXygayWLHHEH9i2NTNcTv2zxhmu5JRJX7OecIauWj1lpDjo1oKmrI5qTuQtGpX5f9fWwrqOq9ATrF0nqXtOke9QFAps9QPacg/l1XEwVK1FDlejU1pIpOTNRQxViHqVQyAqg2jwJc92oSyQXzBgonAwB9RYn2yiCjhuZcQmprnf6YiL4M4CiAjxLRIQATw5rLDqLsU3jfWSOCfNLRNME4jHctGgDkOYIOQT8ysoTGl9ggpZCv86uNaDZSEfx+lHeJ2vY8qJLMttU/9dJa10Y0qxxBozTWIVns8LoalVAOeY86j2+bIzD1bcgrz9T7aUYEGTUkiQjyZLFtROCrcgSmiKBYoxKFE9eXpaVzOQpNRFCmhlSlqQKXVETAGHsXgFcDOMoYiwBsgyuJ7imI4dIFNTRbRJB3KUcJNqfxTM1ks6JxQrQ80QuUDQQvXzRLRYj7inU2Wjz5utL9bE7aQLnstGbQbSOCGgdvW+bq0lDWXGOfV6if7m2E6lwjAlnkp/o71zuedXX7Yr91ETyniMCij6B6vS4iqOUIIv37I6OGpgrJCKCpHzTWVAmKbmOn8tFdyhGY+gi+iTH2l+Ueglp49Afz2tjFwOqwi/uf2cDmNOZa4DPSOUKaYjSNsTmJZ046z4LitFwYSdv8gFgflQy6jSgbUD0F2mjxyO7nSg25dhY3cwSOWkMOuQw1ZWOR+HUoi61LitjSZF3fa0QEfUUCtZH4NtArQS1fFMb66LBXk5Yuykf1ncUC01gtK92MCPR7V9E3KkM9qPVAqLqQgYL+qecIul5H2tHN98lnEmyM50sNmbiP/wXAX6LaQyDAsMccwYHFLk5vhzg/imaOBoAiItgOE2xOopkppllQlFYWxtXFoZWrf2zLOQFh0GvdzJYOpOy0xHPpUO9bmL0jWd3oVFmXG8IkW2euNmqWj5pfm6rSSPf+16vEcnkG43tRp/LSvFRRua/s/TKVjwYeWdNOYq/1KiDVa5DJXE+jFAcWVH0EXKRODAjS0U6AmhpSOoJutWFNFxF0vQ68DlUjgjBROmCgNJPgYkYEjLGfyP7vPKHscsSBhS7COMUTZ0Y4sOA6rbhAWbdocxLjyuXda8Ju1KJbzOStrC+d5mwNM7+GmpSLhQMpVxvFluWcYk+zJpmj2onbXHZaVQUtKCU15TdLstjvyBP9uv15Hd5R6xpV+V791L5z5aONHEGsL0HmCV2JI5DmCPhjkzo1ZBHNDLpeSbJafr2sfHQcJtrRlrZVQETU6DsYh+pZBALL/fkLz5mooX+p+z1j7Bd2djsXF+uZ8X/oxCZuuXJ55ucR0cTmJMLGJMJNhxd3ZH82kBkgG4qmWE+SpKibYXZ1IK4Gne+p0yhn9K1lqOs6TJYRgQP1VU/i2tyLiA9SEWJwJoNbvtcsNJntMJumxIY+AitHh6bnBnilzrTGyQNy+sb3uAJAo3xUkyMACkdQREzmhjKBSZRgdSg/GA4CD1HC8jJr1bzi4vmbjkN1rcDSYP7CcybOQvSy3QzglShkpL8NXIl0T+HAIv9jn9oKcXBx9ohgdcgdwblRhM1JfEE0kyvqCcepZR9Avn6Gk7a4r2iis5V8EM89Sxd01+8gziaT5RPKHHMErqqlLu9LXhrrODynSs3ZRzqzSHJvxYVx0c0wUFUzqTSheLK4Sjvpk8XyiEBFmXD9oOr1KsNeOLEEQGCuGpJ0/46jBFdpcgTi+sDr5AJ0qlP+oNupdFHrKpIElvvzH1dpooZ+CgCI6CMAXsEY28x+/kkAvzvXnV0ErC8UQlEHFs0zflUQjuBs5giWB7uXIyhKHQsaY9ZksTBIqkRWc51bWaZYV2jr2EtmlytTbDp9gaKz1dXxiDkNrvssJypdlE5n0ULKHb+Dc7Mt8ayX3ZqcTb0p0TQopx/wz0DB4+tfA68EqkYQpohAPKexakiiZTTWjJLsl5rElvpBbuRVxn0Y+I2uZSM1NAjw1Lmx9poLha2FuA5AeX5bCPOEsssOZS7/qpXZef1B4KHrdfDM+TGSlO1uRCChhpySxb7EoFtRQzJ6wkKaomz4LBu8gGrfQpza9UoIjZ1Zmu3K0hs2fQRALW9i+V6WqTJTUrZY00zwu6wB9DmCoppJJMv1r7/uZMzJ4lp5ambYZQ1cgBhgXzz/JFYPmqnTWqbOYqFiWm8oUxv26lwBUzlovyanMQ7N1BDPEVxcakjgvwP4DBH9IXi10JsA/PrcdnWRcMVyEQW4jIesg4iwMgzw5JkRgNkb02ZBZS4v7EXjivVNmsGW4nEZwiLQ9Yrh4i7rqpFEahW1iHWFo8siHotmu3K/g0tHcqP3wCIf0aiGMjkcf4YcQaPE076hLDZERIHXQcqQn/B1yqZAU+hNx/kDosqIf2YYY3wQj0Z9VLw+oNysJr+eiNCvVTFNDFVDADCK4uzatPJ44/qg08gRCAZBhZVBcGloDTHGfhbA2wCcBXAOwNsYY/9ujvu6KCAi3HoVTxJ/hcNUMBlWBwEeO80dwZoi0TQPeB0ClUorQ0PFRh2irJAxZvzCV9c1aRCbk/aFJIvL1UbWWkql+wnKQnXyVN0vSlKQRQRSptlsS2q53r/b+yFbY44I6ppS6gNDo2oq0Udgrk2NwohPS8ZaZaiBYmZAeY1OYkI8JwBj1RBQTegyxjJqyFBllJ3yhYOSyWOI6+u0kzFHMPDzmQTzgvVRlTH2OQCfm9tOLhG8722vxCOntnH4Aks+14ZdfOaxMwCAgxeQb3AF57OLJGWYpFiSzE9VIfCKaVTOXL9D5215nesJuH6/KEmNFUMCdQ7e1kkG5SqlLPlpciDlCiDbPEaZgrIZZqNaY4xWSu+7SbG0Lpdh6hQuJ+X7WVWN7jNU8PiFsbaNCEyOQNab0SF9FFgeZhMlvAdBVwUEFJSQSTto0PUwPutODQFcZmJeMwnsOYMZQESvJ6IHiOgYEb1L8vvXENF5IvpC9t+/ned+bHB4uY+vft6BC36elVK4d+ACKpBmQf30aio/rKwtVR3Z8toApyfi1K0xTNzPpSwzv19J4dKF/qpXKdnmT3ql3Elk6JQVqOcVbPIY0vfDQl7bVWLCr71/fL/y+zTEAQ2vXybup6eGqk1c01jNyQNV2QgT1dNIFmc0ks6J90uKorpOYaCggHLHFCUgUh9mZOWj5ohg/jITcyOvicgD8MsAXgcuXX0nEX2QMXZ/7dL/yRj7u/Pax8VCeVTdbkYEQDNJ6ZojAPiX3TVHUB+oYuOApNIUM+QkAidRvWJ4vYsDEclSWwdSN+p2VJl7Q18lwW/dR+DmPOqfKV0EVpfk4NGGrqGsFhFotIMAbkzPjXjtionqqUcE0ygx/h14t3BxvXhMhmL8ZJEsVimViuvrPQrmhjJups/PMWE8z4jgVQCOMcYeYYyFAN6PPShUp8J168WoOlXr/rxQp2lcO4sB/uV1reufJVlcrl5xSd5WO5JTK/qK348qNIHte1M2traUUt3YuuZaXJrDGr0RNnRSo1pLvaZeCmtDDVUiL8eIwEQNuecI7J4b4Px+PSLQaQ0BxchJlQR1+XpRNcSYWrK6jN2ICObpCK4B8GTp5+PZY3W8mojuJqI/I6IXy56IiN5ORHcR0V0nT56cx153HNeXHIFNMnInUT0tu1UNdUtfYltZaL6uZPQcy0DDmsFwT94yK+fB91RdZ0ubzeJA6gl0q8jKLzehpfA76oEuxd5mE+5z6bCuOgL9Z0omW61NFtfGT06iRFkFBFQ5fNPoSZlOkqrUVKAsG2EsB62J1E2iBH3N+ygkJhhjxsY5gd2Qop6nI5B9elnt588BuJ4x9jIA/wnAH8meiDF2O2PsKGPs6KFDh3Z2l3PC0RvWsdjz8X1fd+Ou37veC+DSR1BWwCwqUOxOv6Jk0E1iokQpxfZJ32ok4XqyL+7nsq7sXJ2NuqXmUzWv4LDGMSIIvE4uF26T3HfJrYjniUsOTSsx4UsiAo1xlEUEKmNal63WlZoWawpqSDevGGhOHRtHSd5kJkO/64Exvm+THIXA8i4Mp5knZ3EcwLWln48AeLp8QTb+Uvz7DiL6FSI6yBg7Ncd97QoOLfXwyR/9JqU87jzBdWTKJ1GH8lFJsthXzJqt3rO5zipHUKOUrDn7OjU0Ux+BvQ5T4HXyORW2DqQeEdjcq+tT7SRtSUHVZClMEVJ53KeNYmmlAspAjZU/C2nKjA6tHhGYqobKDWVFjsAUEZRKU61yBJbJ4qCaI5hEqbJ0tHz9JEryITZmiYnLOyK4E8BNRHQjEXUBvBmFVhEAgIiupIwHIKJXZfs5Pcc97SqW+sGu00JA9dTrnCPwii+OrW4NX1eUDLqUj1YoHsthNuJ+VfrL1qDTTI6nPNrR1oF0S41etppPldN9kmrnFQuUneI0O30bS1sljts098A1RxAmaaED5RQR6KmhavmoXla6HhHY5Qg6FaoHUCeLvQ4XChyXrtclf8vqpnnzmcER9AMPXb9zeVYNMcZiInongA8D8AC8lzF2HxG9I/v9ewB8J4AfIKIYwBjAmxljdfqohSMaFTUu5aOlio9CU8cueQvw07KbA+GUUpyk2ilZsn1Wq1jsI4KtTBzPVNZYX+daiVUtubTMETTW2JXS5n9vy9LWcmWPjeOuHC5MOYJS+ahNwrsREcTqITlALSKw7SNIigjCpmoop3pCs7EeBMWUskmkny9QHleZz2Y25AiA+ctMzJW3YIzdAeCO2mPvKf37lwD80jz3sB8hjAkfNekuMQFwYxdbnBSb65hzY1ixzqX6p1rFM+zafZTrJ25b6q5MYdlGWbM4j/o8CJv3o1JplJjLI4FCsjsuRwQ6aqjmoLTUkN/8LOjnEfDrq1VD+oggSfnnxTRxrN4MN43VktIC5WS0KVksflfOEaxoJhKWO5HF58mmSWzew2nmSQ21uEgQxs6lnl+grKPvWgYKiGojh6qhMo3g2Bgm7hOnLtRQlU+3ThZXkrh2TXr1ShtbcbvywHtnKe/Y7r0oU4BTy4jAnhqS0E6GKiCgnCPQ8/hFBJGWZhfIr89nPJSiDSM1lMlcC3kJQN0pDIgIQuQI9MnicgOaSam0jHkPp2kdwR6EOL26dAaX1wKFYebrHXMEWTLRJj9SDKJPnap/glpS1Zoa8qtG05oa8t1LcuvCeLYUW5HfcbhPKX9h63CAWgRnqASyzeXI80WasZ7i1C6MqaFqqFyyqZtmJtCr9WbodIz48xeOxsZYc2pI5AjsksXjUrLYlCMAeC/BPKuGWkewByGmT7kkbYu1TX7XyoCVT4GOHbv5OktjKe6Xi+OljqJzM+RPmslSOydne4qW3ceaGvI7iEryHm6UnB2P3y2VbJr+vuXPgk2PQidLuE7iBHHC5xKYqCEgiwgM1BDAnYRrQxnA6RsraqhUZTTRCNQBJW2iMMnzD3YRgY/NOUYEu1/b2GLu6Pqec2dwsbbG2Vue7OvSFLa9C/W+BRcHwrK+BSfRudLJdmpZngnUcgS2fH+99t5yjejHiOIUPSvnUSjG2ucvCgrQpvegLtan6y1xdTJA1hsQpcbkL1CNCCYGakjcuyxDbVM+CgCTmDuCwCPte8q7hUvJYl1EUOo7EH0WNhHBwcUent2YgDE2l0rENiLYgxBdsLbSx9W1ghpKnOv6xTo3iqccSdhz/WWHFTuLzpUanVxyBI7run65IsqyOazWx2EnwFcoxtrmL4Laffjz6Okb27Lb/LlLelWm184rgRIrR1AuCbW6vtZ3YKwaqpzazaJwXEiuyCnori+0iZJSj4L57/X8w4sYhQme3ZgYr50FrSPYg+jlOQK7aVVldCtfYhfpBnHCZAgduPdiHjBDlLpTSuWchO26SrOdQ+QiOnHdjTqzpnnE65jG9jmTslG3jQgqZcI2EYFXTXzrmgxlOQLT+1xEBGZOvldKLk/jxNg3wfdeThabDTsgav3NMtFCSC5KGFKm7jkA6n0EIhFtjgheeHgRAPDbn37CeO0saKmhPQhBf8yWI6gmfW0NZb1ByXpdqcPV7YRerHPqLG507trus5bUtjypi/vYdgn3akbdVu4acKPlChkIu8ixKVtiTv6WcwSm97mX5QjyTmGDxATAE7OmLmTxXNOYdzmHFhFT/dRuUgcVyeJJbKapyk4mzfI6OipJ4Lbr1/BtL7saK3MactU6gj0IUUViUx9ex24nfSvrLE/azfvpT6hl1MXqXDqZAX6itJaYqFRguSVx85yJA80zTRKEsV1Phe94aq80ujmVj9pFpf3AwzRKS8bUrtzU5oQvoplcSsNAxQxLTV82g2MGXZ4jEBVGOsfBZzEDkzBBynhnss13xfc6+E9v+UrjdbOidQR7ECKx5zJ8XqDSUOYk+eBOT5TXuRi+yv1i0Yhm73gEb+8S8dTLXO1O6rVox+GkPtt9uKKle/moWx+BKScj61q2SRZXIgKLHMEko5JsIoJJZFdhBFS7f20Gx4i+A9EQpzvhE1FJgZTLWVwMGZo62hzBHoRIiIov1UzJ4tiNexdUQSj48Hlz/bXTdmAZEYj7iUSdraRF3WHZDqYp1rjlFcLYPnlel/eYuXzUUDUUisokY0TQbEq0jQhsjHUjIjCc8EVEMLWgboAyjx9b5QgGAa/SE9IlJschHIEN7bRbaB3BHoT40omSNjcZ6sKgO8k0e/wDHTkYSqB80mZOBl0Y8EnEQ2xXimd7mmQ/u+UIXO5XP927rAkvJIpwSEpXVWZ1vH8RdXBno76WiOB3yKnLXEhLm0TkxLVAVjVkaODi12cVSZGdUxLU2sgyRyCoJDE1TddHwH/Ph9NMIjOttVtoqaE9CHEiFCcUl2SxmE/rmvStyBpb8tTlvQk+2YXiAZBPe3IZXg+4vzfiOhcHUta5MRlPgV4pInDRNMrv4yoxEadWiqXifRuHnNIw7Uv0UNgmi4W+z9SCXimmjvGqIVNEIJLFNqWmQKnW3zJHIKp+zo4i494BHhFMogQEumQigtYR7EGIL63rqVdAzDNw1f4BSgbdQRYamK2hDAC2RdTjuM41Wsrf01A4EJvpa0UUYWM8y/sTJae2/Qp8jX1EIK6Js3JYU+Nat/Z+myQ9+IzjIkdg5PH9mrHWSkyUcwQWVUO10lTTKbxcNWQzSlI4ijNZRGAqBx10eZVRh8iqh2A3cGnsosWOonAEwti5JaOEHo9rfT6QUUqzdBY7lo/mXH8onJ3bPoWTtC5Xbbyn5nXCuObRh3NeYYbEr2MUwf/OiXFvQe31m97vpt6VZURgweNXIwK7voBJnBROyWB8A6+DwCOMch7fLqdwdjtzBBb7yXMEFl3Fu4HWEexBBDUDJPh7W4guUld1TsBdPK7aSWsfgQjnJgyT/TyC6jpbKqpOt1lNbZshMgtmyH3IaB7b+4iktMkhFhFB5kBNjiNr3LNNFouIYBKZk8WBR+hQoT5qHxFkz23xnoreABtqSEQQZ7ZFjsAmWZxmswtaR9BiTqifXm2NnQCvGWeOVE2WI4jtG6GAwsBMIi425pqcLqp/3O5XUEp2X8Q63eZyUs+jCIcKoIKCcogiMlkQ58E0FuW3eQGCZW4l8DqILechA7xb2DYiIKJsrnBipSYqksW5LpGF8R12fWxPY6tTuzDmp7amAICFnkWOIEwwspCv2C20OYI9CGFwXIxJZX1GDblo/3gdAlG5ocw2ySwoHjenFeQRgVuyuJ5kdqeUHHIE3gxGvR7NOZzux2EKxuycomu/gmsBglB5tZ1y169HBAb6pl9KAJvnC/DeEfFZsdFiGnQ9nMvUPm04fwA4vcUjggVDocSw62EUxSDQRZlpLsOlsYsWO4pZTq9lBJmapYtBJ+IKjWGmb2+fI+AGYkvQJ679AI7OrsF1O0pouCSnZ8kr5I5qat/n0NibxX28DqdXbBsAXZPlghoS0aGpaUqc0jezKVwm3l9EBFYNZdlziQlfNgnaQeAVVI9lslhEBDZVRtxps5xWutiYKzVERK8nogeI6BgRvUtz3SuJKCGi75znfvYL6qdXl4YysT5K+OBxl7VljSNrw9ypVvG4zCMACs7atUrJPVlczS04UUN5iau983Apb21QgQ4JdyExYTKmjcOFheNwkRsR99+YRAg8gmcQOxQRgU0tvjD858d2Tgbgp3ZXR3B6O8Sw66Fj2LsoH90O40smIpibIyAiD8AvA3gDgFsBvIWIblVc927wIfctdgAX0kcAIB9C4qL9w+9DTqJsAB9K4neoSKg6VhsVFI/byd65fDTLJWw73K8REThQNrNEES5rxH7ixC4XVL+HuXy0cDJWFVaZMT0/jq0MdSUiMPUR+OK5hSOwo4ZOZJLPi329sR6UksU2/TODwMPWNMYksu+3mTfmGRG8CsAxxtgjjLEQwPsBvFFy3T8H8PsATsxxL/sK5TDe65hPV3XkEUFiP/BFrBNljzaVGeX9jlwpnmxfI8u69vK9APf8SVCLCFwqgLYcynhn2V9RJeYW5eRlwhbGutiXXQTGZzzbFxz0RUQwjqwMdS/oYBzxSjMbrSHx3Pxnu4hAvFbTqb1cJWRKFM9y/W5gno7gGgBPln4+nj2Wg4iuAfAmAO/RPRERvZ2I7iKiu06ePLnjG91rKNfKuwjOCYghLC6ibOK+Y4fmqfK6vCzRmeu/sD4CZyoqz7s4GHUXmqdm1GeKPBzyM7a5oDySsoxuAr+DqUODWxER2DmCvu/lJ3xT5c0sEUH5pL5kiAj6peczJYqB6n4X9jo1BED2TWG1n/8fAP+KMZbonogxdjtj7Chj7OihQ4d2an97FuKDvjWNnWkhoETxOJSBAtwAuSZh+f06zvx2neJxl4pw7EiewajX6atZEsw2RqtepuqUI4jtFEsbkYrp+iw6tP0MlSMCm5LKXtDB+ayT1zoimNi/p+VT+2Iv0F7re538Ndqc8IfdS88RzHMXxwFcW/r5CICna9ccBfD+rKLgIIBvJaKYMfZHc9zXnkdQOr3NomVSpnhccwSuSVh+LTmVZYo9Am6cfXlfI8vGqOa62UtBbZyjqOa5MArK3hHYDgSqR0TmKqMi2nCJCM6NI1yx3Ddf73u5YTfKRIuqoXFk1FQSWCoZaFOOAOCOIxzbcf7l5PPCPqgauhPATUR0IxF1AbwZwAfLFzDGbmSM3cAYuwHA7wH4Z60TuHCU+dxZIoKu38k7W12poZHlibF+vzzpa7lOGExRPmo9UtOvGk1XB+KyrkiEu0UfXb/jpGkkejhc1ojrbPWJmjSXOUfgIpwnTum21FAv6JSoIfscgc1zA8DaQjEJzKayRxj31aE+egCq0cbKwHz9bmBujoAxFgN4J3g10JcA/A5j7D4iegcRvWNe921RNVIzOYISZ+/K9ReyFm5J5pGjeBwgqCi3k31ddM729dWrm5yM+gx7dKmiEoqxWzPcx7ZxsC5bYqs+aitAKE71Scqs6vz7vofEctRjuTTVVvZ53dERCAewZjFK8uBicU3Z4VxMzJWgYozdAeCO2mPSxDBj7J/Ocy/7CWVD4NpMBmSG2ZGqEfdyLecU185S6iqS04BDRNBI+s6aA7G/n2uFUreUM3GRzrBN5Jb3ZlsUUMy4sKWGSjkCmxN+6RrbHIHs3zL0S4now0tm2gkA1ksG3abq7kBm3G0igkOLxR7W5zSD2BXzpIZaXCR0LzAiCHyyFher39e1coXfr5NLCzg5nh1IMrtOb3OVva44D4eZCc4VQDOsyedOxKlRc6nevWylPhrbCxCWjb/NMPfyNbYRQZQwo9MQWF90M9CU1cYcXOwZrz20VFyzvNepoRYXD5WIYKby0eKL5ZosniUiKO/RNSIQipKukhau5aPiHhkb4SRNYTvAvdijezltd4aoKsgSulMH0TnbSKWsYGvF+VciArscQfFvu/JRwM7JAMDzDy4CAI6sDayuf8k1KwCA5x1aMF5bzhG49vjMC5dG7VKLHYVoIktSe82fMspGwdUwx5mldOX6Zf823q+0T1tqSCRWwyQFkdsXsfyaXOcYzLzGOiIgbG+55wi2JrFVA6D4TI2sq4aKzmLniMCCGqpEBBYCdQK2EcHKMMB/eetR3HzFktX17/ymF+Cmw4v46hsPWF3/+z/wNdaf2d1A6wj2KAKPssSbe3laxeA5UjwzrZuRyqqss7xfLo5nKYZWRvk1uXZAu62ZwRF4nTx56pIjcGkADDzClmVVWOB1wBiXCZ97jsBCoE52HxNed+sV1tcu9nx8x21HrK+/7fo162t3Ay01tEchjMEsw7ErhsiRs5c9h8v9XHME+TpL1dLyOncxPr63jkMkMYvzmCUimMV5d30eEdiu6WbGHbDLEQCcSnIpHwXsTu0uVJIYZMPXXRp1+5caWkewRyG+iLPMRJ39hF4YB9cks+zf5vvN5kDEta60mbh+3k6u7HxnoZNsT709x7xCt2RETY633Hdh89y+18mNtVWy2IFKEoNsgEtH2+dSQ+sI9ijEF3GWmaiznGLr17oZ9FlzEsUJ3VZ0rnwP9zkN7pFEOfqwpaFmySvMQid1/Q42nWS1+f69DhmllgV1ZCMKJyDeHytqyNHxiSjjUlH7vNTQOoI9iiCPCGbJEcyeLM6fw7E+P3+OGXILMxt014ggW+faNc3v6dZgB/AEuMng1vdW/7cOPb+gelx0kKzmC8zwWRA5DpthLRXhNhuhtywiuFT0/y81tI5gj0JEAhdKDc1q0F2ScjPnCPzZ8iDFulmpIfc8xizJ81kS9V6HnCW5be/lQo2VK7pcczEmtU+g+rezcZaCErpUJoJdamgdwR6FUDWcJSLYCYpn5qqhWZK+zhHBjDmCCzDqTn0Vvvt9RBQ3z7JdF2qs3Ivi+j7bnNpdP9fi+3CpqH1eamgdwR6FODHNRA1VqBp3SgOYzSC5UCHl+80yivNC1s1i1N0a7GbIRczgPMqRlBs1ZH+t654AYKlv7rZ1jeZEzf4sObP9gNYR7FF0ssTbLKHwTtT1u52ayfleQEGHuBqFmXMEvruBnuVeFxJ5uLwXFSrPgkLMS5Itrg1mqHwSsKGGXA84WXN3RUyuRYHWEexRpFkWcBaZ2+6M/O6sjWiCRnAfqXmhFI+bMcmNrYMRmilZ7Lu/rgtxOPV/K693cLwXEhEcXjbr9bgecMSeDzhqCO0XtI5gj0JUYMziCHaij8ClfV4YPteO++4Mp+Dy/VxPqnlvxgzcvcsJVjjGWaghl/eifLK3bShzvdb2egD48b/zIrzwikUcXDA7gusPLOCWK5fwb/7urVbP/QOveT6uWO7h1quWra7fb2gzJ3sUX3ndGj796BkrNcQ6ysbfpeqoXFbpJN2Q3c8lPwBcePmoqwOZJe/SnWGNcFRONM8MUU7ZWNtUXrlUac1SQfZ9X/88fN/XP8/qWq9D+NAPf4PVtQDwjTcfxqd+9LVOn8v9hNYR7FH84Dc+Hy+6aglfcWTFeW35S+x2kp0twBSG2XP8ku62Iwjy0707ZeOSpOzNlCNwp5NmrRqyed9mLRyYJ1onoEbrCPYolvoB3vjya2ZaWz0pzmZYnO5XqoF3QXFyduwjuMBk8SwRwSx9FbMk6k0qomW4irG5VCZd6HCkFruL9i/UooGyAXI5Rc0yBKe8ztURXKh43KzVRrMYdd8pWSzWzJAjcIhWeo7G2ilZPGPhQIuLg7n+hYjo9UT0ABEdI6J3SX7/RiK6h4i+QER3EdHXzXM/LezgUuFSXTerIxDJ4t2lhoaOzUWzRC7CaLq8tlzKwuE+s0gnuPL4LnmIWfSSWlw8zI0aIiIPwC8DeB2A4wDuJKIPMsbuL132UQAfZIwxIvoKAL8D4JZ57amFHWY16K4nbIEiWey4boYmKqA4NS84liDm+5whSnJxHkVEYL9GjDwU1WI2aCOCFgLz/Au9CsAxxtgjjLEQwPsBvLF8AWNsizEhe4UFAPaf4hZzw6wG3WZwtww5J+7oCUQC1jV+KSgeN0cgejNcatiFA7AdkQgUPL+L8xBlwnFi/xVyPbW7lOsGlWij7ea91DFPR3ANgCdLPx/PHquAiN5ERF8G8KcAvkf2RET09ow6uuvkyZNz2WyLAv0ZhbnWhrM164gZrq4nR2H8xlHitE70OMQOp2cAcLwcQBE9pMx+8WLWWRtm85htMEsHefn9tslHuCSxy46lVfy89DFPRyA7zjS+DYyxP2SM3QLg7wP4adkTMcZuZ4wdZYwdPXTo0M7uskUDSzN+cUX7/lc/b91p3TWrfED41at2g8Lr9zuyNnRad8VyP7tf32nda285DAD4xuz/NvjaFxyA3yF860uvsl7zldetoud38A+PXmu95nmH+LD17/36G63XuJ7UXTqey3mmYTsM5pLHPF31cQDlT/IRAE+rLmaMfZyInk9EBxljp+a4rxYGEBHe97ZX4soVN0O50PPxC//wZXjVjW6O4PoDQ/zU33sxXv+SK53Wvfr5B/Djf+dF+Ae32RtMAPiuV/LrXe930xVLOPazb3Cq5rn+wAK+/NOvd1pz1coA9/+fr3eihq5ZHeDBn3mDU1Rlo+lThqCdbGikcrWZzbyAFhcX8/wL3QngJiK6EcBTAN4M4B+VLyCiFwB4OEsWvwJAF8DpOe6phSVec7P9qbeMb3+F/QBvASLCd3/NDc7rAq9j3YlaxlI/mGkd4FbSeSFrXEtpAXdqzZXKE1SVazHBLDMxWuwu5uYIGGMxEb0TwIcBeADeyxi7j4jekf3+PQC+A8BbiSgCMAbwXaXkcYsWLeYI4Wxsnc7ygJuLlYGb2Wg7ei99zDVmY4zdAeCO2mPvKf373QDePc89tGjRQo2f+/aX4gpLCvC7jl6HZ85P8CbLqO8//5PbcG4UXsj2WuwS6HI7gB89epTdddddF3sbLVq0aHFZgYg+yxg7KvtdS961aNGixT5H6whatGjRYp+jdQQtWrRosc/ROoIWLVq02OdoHUGLFi1a7HO0jqBFixYt9jlaR9CiRYsW+xytI2jRokWLfY7LrqGMiE4CeNxhyUEAe13Ebq+/xr3++oC9/xr3+usDLv3XeD1jTCrffNk5AlcQ0V2qbrq9gr3+Gvf66wP2/mvc668PuLxfY0sNtWjRosU+R+sIWrRo0WKfYz84gtsv9gZ2AXv9Ne711wfs/de4118fcBm/xj2fI2jRokWLFnrsh4igRYsWLVpo0DqCFi1atNjnuOwcARFdS0R/RURfIqL7iOiHssfXiejPieih7P9r2eOvI6LPEtG92f+/qfRct2WPHyOiX6RLZKbeTr1GIhoS0Z8S0Zez5/m5i/m6BHbyb1h6zg8S0Rd3+7WosMOf0y4R3U5ED2Z/y++4WK+rtKedfH1vyR6/h4g+REQHL9brKmOG1/gqIvpC9t/dRPSm0nNdkrYmB2PssvoPwFUAXpH9ewnAgwBuBfDvAbwre/xdAN6d/fsrAVyd/fslAJ4qPddnALwaAAH4MwBvuNivbydfI4AhgG/M/t0F8D8vhde4k3/D7LFvB/BbAL54sV/bnD6nPwXgZ7J/dwAc3CuvD3xc7gnxmrL1P3mxX9+Mr3EIwC+tPVH6+ZK0Nflrvdgb2IE/1gcAvA7AAwCuKv0RHpBcSwBOA+hl13y59Lu3APjPF/v17ORrlPzu/wXw/Rf79ezk6wOwCOBvsi/oJeMIdvg1Pglg4WK/hnm8PgABgJMArs8efw+At1/s17MDr/FGAM9lju6StzWXHTVUBhHdAH7S+DSAKxhjzwBA9v/DkiXfAeDzjLEpgGsAHC/97nj22CWFC3yN5edZBfBtAD46z/26Ygde308D+HkAo/nvdjZcyGvM/m4A8NNE9Dki+l0iumIXtm2NC3l9jLEIwA8AuBfA0+AO/b/uxr5dYPsaieiriOg+8NfzDsZYjMvA1ly2joCIFgH8PoAfZoxtWFz/YgDvBvC/iYckl11StbQ78BrF4z6A3wbwi4yxR+ax11lwoa+PiF4O4AWMsT+c5z4vBDvwN/QBHAHwCcbYKwB8EsB/nNN2nbEDf8MA3BF8JYCrAdwD4EfntuEZ4PIaGWOfZoy9GMArAfwoEfVxGdiay9IRZB+e3wfwm4yxP8gefo6Irsp+L/g5cf0RAH8I4K2MsYezh4+Df8EEjoCfSC4J7NBrFLgdwEOMsf9n7hu3xA69vlcDuI2IHgOnh15IRB/bnVdgxg69xtPg0Y5wdr8L4BW7sH0jduj1vRwAGGMPM86b/A6Ar9mdV2CG62sUYIx9CcA2eD7kkrY1wGXoCLJs+38F8CXG2C+UfvVBAN+d/fu7wfk8QYn8KYAfZYx9QlychXSbRPTV2XO+Vay52Nip15j97mcArAD44fnu2h47+Df8VcbY1YyxGwB8HYAHGWOvmfsLsMAOvkYG4I8BvCZ76LUA7p/n3m2wg5/RpwDcSkRCFfN1AL40x61bY4bXeGMWfYOIrgdwM4DHLmVbk+NiJylc/wP/wjPwEPIL2X/fCuAAOP/9UPb/9ez6Hwf3zF8o/Xc4+91RAF8E8DCAX0LWaX2x/9up1wh+8mDgXyzx+PftlddXe84bcAkli3f4c3o9gI9nz/VRANftsdf3juwzeg+40ztwsV/fjK/xnwC4L7vucwD+fum5LklbI/5rJSZatGjRYp/jsqOGWrRo0aLFzqJ1BC1atGixz9E6ghYtWrTY52gdQYsWLVrsc7SOoEWLFi32OVpH0KJFixb7HK0jaNGiRYt9jv8/4MdrZ0j9IQsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Distance you might want\n", "#min_dist = min(df['distance_au'])\n", "\n", "cond1 = df['distance_au'] <= max_dist\n", "cond2 = df['distance_au'] >= min_dist\n", "cond3 = df['datetime'] <= today\n", "\n", "# These are the indices of df that fulfil the criteria, returned as a np.array\n", "idist = np.array(df[cond1 & cond2 & cond3].index)\n", "# Remember that idist[0] is the first index of df that fulfils the criteria')\n", "\n", "# Find the gaps where the criteria weren't fulfilled:\n", "# Take the qualifying indices and take them away from themselves shifted by one\n", "idist_diff = idist[1:] - idist[:-1]\n", "\n", "# Without the [0] it's a tuple\n", "gaps = np.where(idist_diff>1)[0]\n", "\n", "# This will be a list of the indices of qualifying interval starts\n", "starts = idist[([0] + list(gaps+1))]\n", "# This will be a list of the qualifying interval's last element\n", "ends = idist[(list(gaps)+[-1])]\n", "\n", "# List comprehension to extract the relevant slices and store them in a list\n", "# ends+1 to include the end\n", "dfs = [df[start:end+1] for start, end in zip(starts, ends)]\n", "\n", "# Concatenate the list of DataFrames into a single DataFrame\n", "result_df = pd.concat(dfs)#, axis=0, ignore_index=True)\n", "\n", "# List the intervals\n", "time_ranges = [(df['datetime'].iloc[start], df['datetime'].iloc[end]) for start, end in zip(starts, ends)]\n", "for start_time, end_time in time_ranges:\n", " print(f'{start_time} - {end_time}')\n", "\n", "#PLOT THEM\n", "fig, ax = plt.subplots()\n", "ax.plot(cdflib.cdfepoch.to_datetime(t), d/AU)\n", "ax.plot(result_df['datetime'], result_df['distance_au'])\n", "\n", "ax.set_ylabel('distance (AU)')\n", "plt.show()\n", " " ] }, { "cell_type": "code", "execution_count": null, "id": "78ee8d82-9bf5-44f9-8828-197289e05e9b", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.7" } }, "nbformat": 4, "nbformat_minor": 5 }