{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 15 - Mass-Benefit Analysis - Uranus" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from matplotlib import rcParams\n", "\n", "from AMAT.planet import Planet" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "planet = Planet(\"URANUS\")" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "ISP = 320.0\n", "g0 = 9.80665\n", "\n", "peri_alt_km = 4000.0\n", "apoa_alt_km_pp = 903323.04470317 \n", "apoa_alt_km_ab = 4962409.1142301" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "energy_pp = -planet.GM / ((planet.RP+apoa_alt_km_pp*1E3 + peri_alt_km*1E3 + planet.RP))\n", "energy_ab = -planet.GM / ((planet.RP+apoa_alt_km_ab*1E3 + peri_alt_km*1E3 + planet.RP))\n", "\n", "v_peri_orb_pp = np.sqrt(2*energy_pp + 2*planet.GM/(peri_alt_km*1E3+planet.RP))\n", "v_peri_orb_ab = np.sqrt(2*energy_ab + 2*planet.GM/(peri_alt_km*1E3+planet.RP))" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "def compute_v_arrival(vinf):\n", " return np.sqrt(vinf**2 + (2.0*planet.GM)/(peri_alt_km*1E3+planet.RP))" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "def compute_DV(vinf):\n", " v_arrival = compute_v_arr(vinf)\n", " DV_pp = v_arr - v_peri_orb_pp\n", " DV_ab = v_arr - v_peri_orb_ab\n", "\n", " return DV_pp, DV_ab" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "def compute_fP_prop_pp(vinf_array):\n", " v_arrival_arr = compute_v_arrival(vinf_array)\n", " DV_arr = v_arrival_arr - v_peri_orb_pp\n", " \n", " return 1 - 1.12*(1.0 - np.exp(-DV_arr/(ISP*g0)))" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "def compute_fP_prop_ab(vinf_array):\n", " v_arrival_arr = compute_v_arrival(vinf_array)\n", " DV_arr = v_arrival_arr - v_peri_orb_ab\n", "\n", " return 1 - 1.12*(1.0 - np.exp(-DV_arr/(ISP*g0)))" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "vinf_array = np.linspace(0.0, 20.0E3, 11)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAGJCAYAAABVbT4SAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACIyUlEQVR4nO2dZ3hURReA3wktIYTeOxoEERUF6UroCEgTRAUElCJgAREbKCBY+BAVkY50RHoTQVCaNBUUBRQRqQICSu+EnO/H7Ca7yaZssskm4bzPM89mp92zN8k9OzOnGBFBURRFUVKSAH8LoCiKotx6qPJRFEVRUhxVPoqiKEqKo8pHURRFSXFU+SiKoigpTkZ/C5CS5M2bV0qWLJno8adOnSJfvny+E+gWQu9d0tD7l3j03iWepN677du3/ysiHie4pZRPyZIl2bZtW6LHV6pUKUnjb2X03iUNvX+JR+9d4knqvTPGHIqtTbfdFEVRlBRHlY+iKIqS4qjy8YJu3br5W4Q0i967pKH3L/HovUs8yXnvjL/D6xhjigKvApWAe4EgoJSIHEzA2ADH2O5AQeAP4G0RWeCpf6VKlUT3fhVFUVIGY8x2EankqS01rHxCgceAM8B3Xo4dAgwCPgUeBrYC84wxjX0poKIoiuJbUoO12wYRKQBgjOkCNEjIIGNMfuBl4H0R+cBRvdYYEwq8D3yVHMIqiqIoScfvKx8RiUjk0IZAZmBmtPqZwN3GmFJJEkxRFEVJNvyufJLAXcA1YF+0+t2O13IpK46iKIqSUFLDtltiyQ2clZgWE6dd2t04deoUlSpFnX1169ZNLWEURVF8yIQJE5gwYYLzbd7Y+qVl5eM1+fLl895b9zgwEfgGMIAA9YEuQCEfC6goipLGcf1Sb4z5N7Z+aXnb7QyQ0xhjotU7VzynSQoCfAC0B8oD3wLrHa93Oeo/cPRTFEVRvCItK5/dQBbg9mj1zrOe35I0+wjgELAaaAVkctRncrxf7WgfkaSrKIqi3JKkZeWzErgBtItW3x7YJSIHEj3zcWAFMJLY71CAo32Fo7+iKIqSYFLFmY8xprXjx4qO14eNMaeAUyKy3tEnHJgmIs8AiMhJY8yHwOvGmAvAT0BboA7QLEkCTQJ6Eb9qDgB6Ap8BA5J0RUVJHbRqBae92LHOnRsWLkw+eZR0S6pQPsC8aO/HOF7XA2GOnzM4iiv9gYvAi0SF13lMRL5MkjSrgdei1a0HdmHVm6v9RjPgE1T5KOmD06dh3bqE9w8LSy5JIpk6dSqdO3eOfJ8tWzZuu+02unbtyrPPPkvGjKnlMRbFoEGDGDx4MN6ELwtz3Mt13tz/NEyq+K2JSHSjgQT1EZGbwFBH8R2GqDMeJ8HARuANoBbQAXgECPTplRVFiYV58+ZRtGhRzp8/z7x583j++ec5efIkb7/9tr9F8wljxoyJv1M6IlUon1SHYE+TXBVQJWA2cB5YCIzHhjNtAZwDIkjbJ2iKksqpUKECoaGhADRo0IB9+/YxcuRIj8rnxo0bZMyYkZjGsKmXcuVuLb94fVx6oj6wLJa27EAnrN/Pr0A48A9QCrsqSpqNnaIoCeSBBx7g/Pnz/PDDDxhjGDNmDK+88gqFCxcmS5YsnD17FhHho48+okyZMmTOnJlChQrx3HPPcf78ebe5jDH079+fd955h6JFixIUFMRDDz3Ejh073PqVLFmSTp06xZDFGMOgQYPilHfkyJHceeedBAUFkStXLipVqsSiRYsi28PCwiK33v755x8yZszIJ598EmOe//3vf2TKlIlTp05F1i1cuJCqVauSNWtWcubMSZs2bTh8+HDcN9DPqPLxRBdgNHY1ExeFgaPADqyyuoFVXBWBj7BKSVGUZOHAgQNkyJCBbNmyAfDOO++wd+9eJkyYwKJFiwgMDKR///689NJL1K9fn2XLlvHKK68wdepUmjRpQkSE+z/49OnT+eqrr/j000+ZOnUqJ06coG7dupz2xgAjFmbNmkXfvn154okn+Oqrr5g1axatW7eOde6CBQtSr149Zs6MHroSZsyYQaNGjciXLx8A48aN49FHH6VcuXLMnz+f8ePHs2vXLmrVqsWFCxeSLHuyISK3TKlYsaIkmOEi8pyI3Iyl/aaI9HL0cyVcRFaLSEcRySkiDUVkpohcTPilFcXngP9LIpkyZYoAsmfPHrlx44acPn1axo0bJwEBAdK8eXM5cOCAAHLfffdJRERE5Lj//vtPMmfOLB07dnSbb8aMGQLIkiVLXG4PkidPHrl4Meof9cCBA5IxY0YZMGBAZF2JEiVizOccP3DgwMj3AwcOFFw+c69eveS+++6L83PWqlVLatWqFfl+5syZkZ/byc8//yyAzJkzR0RELly4INmzZ5fOnTu7zbV//37JlCmTfPTRR3FeM7kBtkksz2Nd+cRGX6AEdiWzALuqwfG6wFFf0tHPlQxAPWAq8DfQEZgFFAGewlrS3fRCjuPA28BDWEOHh7BZjNS3SPGGhKqIWrW8Uyne9E8iZcuWJVOmTOTOnZuePXvSrl07Jk+eHNneokULtzOerVu3cv36ddq3b+82z+OPP07GjBlZv369W33jxo0JDg6OfF+yZEmqVq3Kli1bkiz7Aw88wI4dO3j++ef55ptvuHz5crxjWrZsSbZs2ZgxY0Zk3YwZM8iRIwfNmllvki1btnD+/HnatWtHeHh4ZClWrBhly5Zlw4YNSZY9uVDlExsGmy1oJvYcpx724V8P+N1R/7KjX2wEA09gMwv9gd2OewMo5hi7g9jD82h4H0VxY9GiRfz444/s2bOHS5cuMX36dHLnjoofXKiQe7BF55ZW9PqMGTOSJ0+eGFteBQoUiHHNAgUKcPTo0STL/tRTTzF27Fi+//57GjZsSO7cuWnVqhUHDx6MdUzWrFl59NFHmTVrFiLCzZs3mT17Nm3atCEw0JrZnjx5EoB69eqRKVMmt7Jz507++++/JMueXKi1W3wUAt50lKRQAOuN9CJWec3CWsqFYBVJO6CoS3/X8D6uXxGc4X1aOOYagVVkipLOKV++fKS1myeiW7Y5FdM///zDXXfdFVkfHh7Of//956a4AE6cOBFjzhMnTlCkSJHI94GBgVy/ft2tT0Ie8MYYunfvTvfu3Tlz5gyrVq2ib9++tG3blu+//z7WcR06dGDatGls3LiRK1eucPz4cTp06BDZnidPHsD6Qrl+RichISHxyuYvVPn4gzuxnklvY32HZgL3ABWw/kM1sGF7oiseV5zhfepjFZdG2FYUN6pWrUrmzJn54osvqFu3bmT9nDlzCA8Pj7Qsc/LVV19x6dKlyK23gwcPsnXrVl57LcrjvESJEuzatctt3PLly72SK1euXJFKZ/z48XH2rV27NkWLFmXGjBlcuXKFkiVL8uCDD0a2V69enZCQEPbt20fHjh29ksPfqPLxJwHYM5yHsFESlgMzgB5Yv6KvsPlaozu8uo7X8D6K4pHcuXPTt29f3nvvPYKDg2ncuDG///47AwYMoGbNmjRp0sStf1BQEA0aNKBfv35cu3aNgQMHkj17dvr06RPZ5/HHH+fpp5+mT58+NG3alF9++YWpU6fGK0u3bt0ICQmhWrVq5M+fn7179zJjxgwaNGgQ57iAgADatWvH+PHjuXHjBn369HFb4WXPnp3hw4fTq1cvTp06xcMPP0yOHDk4evQo69evJywsjCeffNK7G5dCqPJJLQQCjzpKNeBx4D3gaWxInw7AA8Q8Y9LwPoovyZ3bu5A50bauUhvvvPMO+fLlY9y4cYwZM4Y8efLw1FNP8d577xEQ4L6t8NRTTxEcHMxzzz3Hv//+ywMPPMAXX3zhtj3XsWNHjhw5wmeffcb48eN58MEHWbRoUZzbgQA1atRgypQpzJgxg3PnzlG4cGHat2/P4MGD4/0MHTp0YNiwYZE/R6d79+4UK1aM4cOH8/nnnxMeHk6RIkV48MEHqVChQgLukn8w4gMrlLRCpUqVxOtkcv6gFta4AOAv7PnQZKA1MJyYCsi1v6IoXuN0Mh061LeRum51jDHbRaSSpza1dkuNOMP7gM1W9BY2Zve3wCu4W7jdQFEUJc2hyic14im8T25sSB9nxG2nAlrq6K8oipKG0DOf1EgXrPl1C9y/HuTBKqC6jvqh2OQTs1JYPkVJZ9xKxw+pBV35pEYKAQ9j/Xiix5fLi91++xKoAjTCZjJSFEVJQ6jySa3EFd5nPZALG9T0HBrlQFGUNIduu6VWnOF92mHTertGVq8PzMHGkauD/S0OSmH5FEVRkoAqn9ROfOF91gC1scpqYEoJpSiKkjRU+aR18hOlgAJIegw6RVGUFECVT3qgAO4KqL9/xVEURYkPVT7phYK4b8G94V9xlHTAcWAi1rzfYA1b6mNdATSQrZJEVPl4olUr8CZ1bu7csHBh8smTUAphFVAYdgX0Wpy9FcUzgk3VsQLoBbyODW57A+v83B7rCtCXuPNZ+ZiuXbsyadIkevfuzUcffZRyF3YhLCyM8PBwNm7cmOzXKlmyJDVr1vSYStvJoEGDGDx4cJr0U1Ll44nTp2HduoT39yYQY3JTGFhL1BbcK/4VR0mDpMJcUleuXGHu3LkAfP755wwfPpyMGfXx1aVLFxo1auRvMRKF+vmkR4pgFdBEbLZTRUkox7ErnpHEn0tqBSmWzn3x4sWcP3+exo0bc/LkSVauXOmTeUUkRnK45OLatWs+n7No0aJUrVrV5/OmBKp80itOBTQO+NDPsihph0nYrbb4ngyuuaRSgGnTppErVy6mTp1KUFAQ06ZNi9Fn4cKFVK1alaxZs5IzZ07atGnD4cOH3fqULFmS9u3bM3nyZMqWLUvmzJkjk8GtXLmSatWqERQURI4cOWjRogV//PGHR3mWLFlC+fLlyZIlC2XLlo1clTkZNGgQxhh27dpFw4YNyZYtG4899hgAq1atonHjxhQqVIisWbNSvnx5RowYwc2bN+O8Bzdv3qRbt25kz56db775xu06rhhjGDBgAJ988gmlSpUiJCSEWrVqsXv37hjzDRgwIFKOOnXqsGfPHowxDBo0KE5ZfIEqn/RMUawCGg187F9RlDTCauCRBPZt5uifzBw7doxvvvmGtm3bki9fPlq0aMGyZcs4c+ZMZJ9x48bx6KOPUq5cOebPn8/48ePZtWsXtWrV4sKFC27zrV27lg8//JCBAweycuVK7rnnHlauXEmTJk3Ili0bc+bMYezYsezatYuaNWty9OhRt/H79u3jhRdeoG/fvixcuJDQ0FAef/xx1q5dG0P25s2bU6tWLZYuXRqZlG7//v3UrVuXyZMns3z5cjp27MigQYPo3z92M9UrV67w6KOPsmTJEtatW0e9evXivGczZ85k+fLljBw5kilTpnD48GGaN29OeHh4ZJ+BAwfy7rvv8tRTT7FkyRIaNGhAs2bN4pzXp4jILVMqVqwoCaJWrZh1Y8aI1K8vMmGCyKlT8fdPTRwSkVIiMtLfgih+w+//fYkXfdiwYQLI5s2bRURk5cqVAsjYsWNFROTChQuSPXt26dy5s9u4/fv3S6ZMmeSjjz6KrCtRooQEBQXJ8ePH3fpWrFhRQkND5caNG27jM2bMKH369Imsq1WrlgCyZcuWyLrw8HApU6aM1KxZM7Ju4MCBAsjHH38c52eLiIiQGzduyNChQyVnzpxy8+ZNN1nbtWsnp0+flho1ashtt90m+/btcxvvvI4rgISGhsr169cj6+bNmyeAbNq0SURETp8+LcHBwdKjRw+3sSNGjBBABg4cGKfcCQXYJrH8RejKJ6E89RR07QqrV8Ptt0ODBjBpEvz3n78li5/i2BXQR8CnfpZF8Q+SwPIgcD2Bfa9jU8AndO5EMm3aNEqXLk21atUAqFevHoULF47cetuyZQvnz5+nXbt2hIeHR5ZixYpRtmxZNmzY4DZf1apVKVgwKhrvpUuX+Omnn2jbtq2bEUOpUqWoUaMG69e7Z2osVqyY2zlLhgwZaNOmDT/88AMREe6RgFu2bBnj8xw/fpzu3btTokQJMmfOTKZMmRgwYABnz57l5MmTbn2PHTtGzZo1uXz5Mps3b+b2229P0D2rX78+mTJlinx/9913A0RuQ+7cuZNLly7Rpk0bt3GtW7dO0Py+QJVPQgkOhjZtYO5cOHbMKqJVq+C22+DXX1O/IiqBVUAfYLfhFMUTnnJJxUYK5JLatm0bv/32G61ateLs2bOcPXuWCxcu0KpVK7Zu3crevXsjH9j16tUjU6ZMbmXnzp38F+3/slAhdyelM2fOICIx6gEKFizI6WhuFwUKFIjRr0CBAly/fp1Tp07Fea2IiAiaNWvGl19+yYABA1izZg0//vhj5Jbb1atX3fr/+uuv/Pbbb7Rt29bjdWMjd7T05lmyZHGb//hxaymSP3/+eD9bcqG2ionBqYjatIFLl6ByZfj6a+jbF6pWtfUtW0KePP6W1J2SWAUUhv3a0cOfwiipkthySUUnghTJJeVc3QwbNoxhw4bFaJ8+fTo1a9YEYOrUqdx1110x+oSEhLi9j35AnytXLowx/PPPPzHG/vPPPzEe5CdOnIjR78SJE2TOnJl8+fLFea2//vqLbdu2MWPGDNq3bx9Zv2yZZ43fqFEj7r33Xl599VUCAwN58cUXPfbzFqdSPHnypNs98/TZkgtVPkklOBjy5YN586wiWr7c/uxURI89Bi1apB5FVAp3P6Du/hVHSWW45pKKzdw6AnjB0S8Zc0ldv36d2bNnU6VKFd5///0Y7X369GHGjBn069ePkJAQ9u3bR8eOHb2+TnBwMBUrVmTevHkMGjSIDBkyAHDo0CE2b97M888/79b/yJEjbN26NXLr7ebNm8ybN4/KlSsTEBD3ZtLly5cB3LbEbty4waxZsWvxfv36kSFDBnr37k1ERESk4UJSuPvuuwkODmbevHnUrl07sn7evHlJnjuhqPLxJcHBVtk89pi7InrpJahWza6IEqKIkjvCwm24h+LplvChyi1AX6wDaX2sOXUzoiIcLMWueJwRDpKR5cuX899//zFixAjCPDhyd+/enR49evDTTz8xfPhwevXqxalTp3j44YfJkSMHR48eZf369YSFhfHkk0/Gea0hQ4bQpEkTmjZtSs+ePbl48SIDBw4kR44c9O3r/kELFChA27ZtGTx4MPny5WPs2LHs3buXsWPHxvuZ7rzzTkqUKEH//v3JkCEDmTJlSlC0hpdeeokMGTLQp08fIiIiYsjkLbly5aJ37968++67hISEUK9ePX766Sc++8zazsenRH1CbJYI6bEkydotKf0vXBCZM0fk0UdFsmcXadhQZNIkkX//TZnrx8afIlJMRCYmbriSzjkmIm+LyEMuZYijPgVo3ry5hISEyKVLlzy2nz17VoKCgqRjx44iIrJ8+XIJCwuTkJAQCQoKktDQUOncubPs3r07cozTgswTK1askKpVq0pgYKBkz55dmjVrJnv27HHrU6tWLalRo4YsWbJE7rrrLsmcObPccccd8sUXX7j1c1qhuVrPOfn555+lRo0aEhQUJEWKFJE333xTJk6cKIAcOHAgTlk//fRTMcbIsGHD3K7jCiD9+/d3qztw4IAAMmXKlMi68PBweeONN6RAgQISGBgotWrVkk2bNiXISi+hEIe1m7HttwaVKlWSbdu2xd8xOVceFy9GrYhWr7YrIufWnHNvOSzM+/A+3vR35U9sQrrBwNOJm0JRlPTB/PnzadOmDRs2bODBBx9M8nzGmO0iUslTm267eSI5g4RmywZt29riqoj69IHq1e3WnIsjWLJTGvgWq4AM0DnlLq0oiv/4/vvvWb58OVWqVCEwMJDt27fz/vvvU7Vq1UgjjuRElY8/8aSI5s6FrVvh4YejVkS5ciWvHHcQpYACAO/PbBVFSWNky5aNDRs2MHr0aM6fP0/+/Pl57LHHeO+992JY6SUHuu2WGnnwQejVy66IvvkGatSwK6InnwSHvb4bSdl2c2UPUBd4H+iA5nNRFCVJxLXtpk6mqZEMGeDxx2HBAjh61EZXmDEDmjSxK6TkoixW0bwCPIH19yiPXRWtd7ze5aj/gCR5rSuKcmujyie1ky2bVUSrV0PJklCvnnfGEN5yJ1a5LMWe/7TCmtlCVD6X1dh8LyOSTwxFUdI3qnzSChkywMSJdkvuoYdsiJ/k4DjwE7AVuwL6wkMfP+RzURQlfaHKJy1hDPzvf9CunVVCf/3l+2s487ncDawC+gBzPPRL4XwuiqKkL1T5pDWMgddfh3797Aro1199O79rPpfyWAXUG/AUdSOF8rkoipL+UFPrtMqzz0LOnFC/PhQp4rt5DVFnPGBXQF9jreBCgftc2lz7KYqieIEqn9RI7tzWfDohFCgAu3fb9A4NGiT92oKN4eWqWO7B5gLqAGwDAh31N5J+OSX1c/Uq9OwJY8ZAYGD8/RUlIei2W2pk4ULrt5OQ8uuv8O230L49zJ+f9GvHls+lHdYU+02XuhTI56L4n0mT7Pebz/x8vte1a1eMMT6J6pyWGDRoEGvWrPG3GD5HlU96oGZNu/J54YWkPyG6YJPNRUSrN8BYbP6WDUTlc+mStMspqZurV62L2apVMH26fe8Prly5wty5cwH4/PPPCU/JEFR+ZvDgwap8lFRMhQp2JTRkCHzwQeLncc3nEl0B5QPGY8PvPEuy53NR/M+kSdbHOUcO6NDBf6ufxYsXc/78eRo3bszJkydZuXKlT+YVEa5fv+6TudIS165d87cImlIh3XHkiEjZsiKvvy4SEZG4OSJEZLiI1BGR+SJy3VF/3fG+kIhUdvRT0i1XrohUrmxfPb1PSRo2bCi5cuWSkydPSlBQkLRu3TpGnwULFkiVKlUkKChIcuTIIa1bt5ZDhw659XGmKfjss8+kTJkykjFjRlm4cKGIxEyp0Lx58xgpFUREFi5cKNWrV5fg4GAJCQmRBx54QJYsWRLZPmrUKKlatarkypVLcuTIIVWqVJEvv/zSbQ5nioPRo0dLnz59JF++fBIUFCRNmjRxS6uAPYV1KwMHDhQRm96hlod0KiVKlIhMMyEiMmXKFAFk/fr10rp1a8mRI4fce++9IiJy48YNeffdd6VMmTKSOXNmKVSokLz00ktyxUe/ZOJIqaArn/RG0aLw3Xd2n6RnT7h50/s5DPAyMBP4DagH1HK8/g6sA04BX/pIZiVV4lz1OI0MAgP9s/o5duwY33zzDW3btiVfvny0aNGCZcuWcebMmcg+48aN49FHH6VcuXLMnz+f8ePHs2vXLmrVqsWFCxfc5lu7di0ffvghAwcOZOXKldxzzz2sXLmSJk2akC1bNubMmcPYsWPZtWsXNWvW5OjRo5FjR40aRatWrcifPz/Tpk1j3rx5tGzZkoMHD0b2OXjwIF26dGHevHnMmTOHSpUq0bRpU4+rtffee48///yTKVOmMHr0aLZv306DBg24ccNa82zZsgWATp06sWXLFrZs2UKXLonb627Xrh2lSpVi/vz5kZlh27dvz9ChQ3nyySdZvnw5r7/+Op999hnt2rVL1DW8IjatlB7LLbHycXLunE0y17atyLVrvp9/vdgV0EnfT634HvB/SSzDhg0TQDZv3iwiIitXrhRAxo4dKyIiFy5ckOzZs0vnzp3dxu3fv18yZcokH330UWRdiRIlJCgoSI4fP+7Wt2LFihIaGuqW/G3//v2SMWNG6dOnj4iInDt3TrJlyyYtW7ZMsOw3b96UGzduSP369aVZs2aR9c6Vz5133ik3b96MrN+4caMAMmnSpMg6PCSHE/F+5dO7d2+3fhs2bBBApk2b5lY/c+ZMAeTnn39O8OeMDXTlcwuSPTusWAGXL9u0DI7c8T7jIawF3LNogNE0gLeqYtQo+PRT79viKoll2rRplC5dmmrVqgFQr149ChcuzLRp0wC7Ojh//jzt2rUjPDw8shQrVoyyZcuyYcMGt/mqVq1KwYJRh5WXLl3ip59+om3btmTMGOV9UqpUKWrUqMH69esB2Lx5MxcvXqRbt7jzzm/fvp2mTZtSoEABMmbMSKZMmVi9ejV//PFHjL6tW7d2S1ldo0YNihYtGrni8SUtW7Z0e79y5UoyZ85M69at3e5bA4fLRvT75mtU+aRngoJsZOw8eawP0Nmzvp1/CPAH1gJOSTc4LdyeecZze5cuKWf5tm3bNn777TdatWrF2bNnOXv2LBcuXKBVq1Zs3bqVvXv3cvLkScAqpUyZMrmVnTt38t9//7nNWaiQez6QM2fOICIx6gEKFizIaUcgX+c8RYsWjVXeI0eOULduXU6fPs2oUaPYvHkzP/74I40aNeKqhxtWoEABj3WuW32+IvrnO3nyJNevXyc4ONjtnuXPnx8gxn3zNepkmt7JlAmmTYPeva3j6tdfW8dUXxAIzAAaYs+EivlmWsW/RD/riY7r2U+vXskri3N1M2zYMIYNGxajffr06ZFZN6dOncpdd90Vo09ISIjb++iJ0nLlyoUxhn/++SfG2H/++YfcjvT2efPmBeDo0aOUL1/eo7wrV67k3LlzzJ07101JXY5l5+HEiRMe6ypUqOCxvyuBgYGcP38+Rv3pWKLeR//cefLkITAwkO+++85j/8KFC8crQ1Lw+8rHGFPMGDPfGHPOGHPeGLPQGFM8gWOLG2OmGWMOG2OuGGP2GmOGGmOCk1vuNEVAAIwcabffHnwQDh3y3dz3Yc2yOxPTNFtJc8S36nGSEquf69evM3v2bKpUqcLatWtjlAoVKjBjxgyqVatGSEgI+/bto1KlSjFKmTJl4rxOcHAwFStWZN68edx0MdA5dOgQmzdvJswRbaR69epky5aNCRMmxDqXU8lkyhQVImTv3r1s2rTJY//58+cTERH1j7Np0yb+/vvvyC1GgMyZM3PlypUYY0uUKMHevXvdTMU3bNgQw8AiNpyrsXPnznm8b8mtfPxqAABkBf4EdgEtgObATuAvIDiescHAXuAA1vOkNjYJwBVgjqcxt5TBQWx8/LFIsWIiv/3muzlviEgVEfnEd1Mq/mHUKJFPP/V938SwcOFCAWTq1Kke28eOHSuArFmzRsaNGycZMmSQ7t27y+LFi2Xt2rUyc+ZM6dq1q8yaNStyjNPUOjorVqyQgIAAadSokSxdulQ+//xzKV26tOTNm1eOHj0a2W/UqFECSKtWrWTBggWyatUq+d///ieffGL/+Hft2iUZM2aUBg0ayNdffy1Tp06VEiVKSKlSpaREiRKR8zgNDooWLSqPPPKIfPnllzJlyhQpWLCglC5dWq5fvx7Zt0KFClKmTBlZtWqV/Pjjj5HyrFmzRgBp166drF69WiZMmCBly5aVHDlyeDQ4+PPPP2N87ieeeEJy5swpb7/9tqxcuVJWrVolEyZMkBYtWsgff/yRsF9UHBCHwUFSlUd54FEgNJHjXwRuuo4HSgHhwEvxjG2APepuEK3+fcf4rNHHqPJxMG2aSIECIj/+6Ls5/xCRPCLyu++mVFKep58WqV7dGkrGV6pXt/2Ti+bNm0tISIhcunTJY/vZs2clKCgo8kG7fPlyCQsLk5CQEAkKCpLQ0FDp3Lmz7N69O3JMbMpHJKafT7NmzTz6+cybN08qV64sgYGBEhISIpUrV5Zly5ZFts+ZM0fKlCkjWbJkkXLlysns2bOlY8eOHpWP088nb968EhQUJI0bN5b9+/e7XW/jxo1y//33S5YsWdz8fERExo0bJ6GhoRIYGCjVqlWTbdu2xWrt5kn53Lx5Uz7++GO55557JEuWLJI9e3a55557pF+/fnL27FmP98kbkkX5YLO5nAW+d7w+mog5vgU2eahfD6yPZ2xTh/KpGq3+NewGUIyVkyofFxYvFsmXT2TtWt/NOVpEHhC7ElIUJVacymfixIn+FiVZiUv5JOXMpydwj4hUAR4EXkrEHHdht9yisxsoF8/Yb7BbdsOMMeWMMdmMMXWwq6lxInIpEfLcOjRvDnPmwGOPwdKlvpmzB5ATeM830ymKkn5JivK5ISKHAURkJxCUiDlyA2c81J8GcsU1UESuAjWxn2E3cAG7kvoSeM7TmFOnTrkdqMV1cHhLULs2LF8O3brZ0+OkYoDJwChge9KnUxQl7TFhwoTIZyyQN7Z+STG1DjTG3I195AAEub4XER+n2HTHGBOITfCcH5tp5jBQGXgLe+bTI/qYfPnysW3btuQUK+3xwAOwZg00bGj9gF54IWnzFQU+xv5GtpO4rySKks4pWbKk85gg3dGtW7dIR1xjzL+x9UuK8gnCZnRxxflegNsSMMcZPK9wYlsRufIMEIY1VvjLUbfBGHMOmGCMGScivyRABqVcORsPrn59OHMG3nrLputOLE8Ai4H+wIe+EVFRlPRFopWPiJT0wfV3Y899olMOG9IyLu4GzrgoHic/OF7vBFT5JJSSJWHjRihTBiZOhNDQhI3Lndsmv3PFmfvnHuARrBG8oiiKC/EqH2PMaaCeiPxkjJkMDBGRAz66/lLgA2PMbSKy33G9kkANrNVaXPwD5DLGhIrIPpf6Ko5X38enSO8UKABOz+0SJawLe8Z4/kRiS/edB5iIdT79FcjuMykVRUkHJMTgIBjI4vi5EzalmK+YCBwElhhjmhtjmgFLgCPYtGUAGGNKGGPCjTFvuYydijUy+MoY09EYU9sY0w/4AHva4NmlWImbjBltOoaTJ+HRR5Pmwt4YG3qnt49kUxQl3ZAQ5XMI6GqMCXO8v88Y81BsxZuLO8yh62AjFczAhqg8ANQRkYsuXQ2QwVVeETkIVAV2AEOBr4CuwASgvohosJfEkjUrLFliA5M2bgwJDNfhkRFYr60lvhJOUZT0QELOfN4nKnmyAGNi6Wcc7Rm8EcBhrv1oPH0OEmVV51r/G/CYN9dTEkjmzDBrlo0cWaeOTc+QN1arydjJBkwD2gDVsLaJiqLc8sSrfERksjFmBXAHsBZ4AZvPUknvZMgAY8fCG29ArVqwZYvNE+QtNYGngO7AQjx8jVAU5VYjQdZuInIcOG6MmQYs96HBgZLaMQbeew9On7bOqLNnJ84M+23gAWA6dg2tKMotjVcRDkSksyqeW5SPP4bff7dm2IkhC/ZU72XsKaKiKLc0fs/no6QRgoJsLLj+/eHXRAavuBfoi+b+URRFlY/iBWXLwocf2mCkFy/G398T/YBr2PhviqLcsqjyUbyjQweoXh169IDExKbKgLV+G4KarSjKLYwqH8V7Ro2Cn36CqVMTNz4U65nVAbjhO7EURUk7JCWwqJIeyZ079pA5rmTNaq3fatVK3HW6Y4OPvgMMStwUiqKkXXymfByhb44DMxy5dpS0SPQgoXExeTKMGAGXL1tl5A3O3D/3AU2wZtiKotwy+HLbbRA2EsIhY8yrPpxXSa107gz33w/PP5+48YWBT7Dbb1d8KJeiKKkeXyqfUkB54E1sugMlvWOMjYCwaRPMnJm4OdpiVz+v+1IwRVFSOz7bdhMRp+vgb9jgnsqtQLZsMHcu1K0LlSpZc2xvGY3N/dMMG2ZWUZR0j1q7KUnnnntg6FBo2xauJGL/LDcwCet8es7HsimKkirxWvk4cuesNMb8ZozZH61Ezyqq3Cp06wZ33gl9+iRufCNs/p8XfCmUoiipFa+23YwxbwKDgV3YPDrXkkEmJS1iDEyYYA0Q5syxqyBvGQ5UwEa+buVb8RRFSV14e+bzDDBSRBL59VZJ12TPbs9/GjaEihUhNNS78dmwUa9bYROpF/C9iIqipA683XbLAyxLDkGUdML998PAgXblcy0RC+PqwNPYnLTHsKkYHgJqOV6HYL3JFEVJ03irfNZjYxMrSuz06gUlS8LLLydu/EBgG1APa7z/LfYv71vgLqA98AE2b66iKGkSb5VPb6CzMeYpY0xeY0xA9JIMMippDWPgs89g+XJYsMD78aOwK52TwP1AJkd9JuyW3GpsTqARPpFWURQ/4K2y2Iv9LjoFOIENC+larvtUOiXtkjMnfPGFjX59wIv8g8eBFcAs4BWgEzFz/wQAIx39dAtOUdIk3hocvI1udigJpXJleP11e/6zcSNkzhz/mElAL6yC6QssBT4GXorWLwDoCXwGDPChzIqipAhGEpOTJY1SqVIl2bZtm7/FuLUQgRYt4PbbbSK6+HgIe7bj3Gr7C6iCPfO5K1rfG9hzofW+ElZRFF9ijNkuIpU8tekZjZK8GANTptizn6VLE9CfKMUDcDvwLjb6wc1ofTOhKEoaJTERDgoZYz4wxvxojPnL8fo/Y0zB5BBQSQfkzg2zZ0PXrnD4cNx9hZgJ5rpgFc2kaPWaiE5R0ixeKR9jzB3YyAYvABeBHxyvLwI7jDGlfS2gkk6oXh369oXHH4cbcWiN+sT0JAsAxmDjpZ9yqV/q6K8oSprD25XPMOA8cIeI1BaRJ0SkNnAHNiTkMF8LqKQjXn7ZWsENiMNCoAs2ynV0C7d7gXbAa473EViF1MXnUiqKkgJ4q3xqA2+KyEHXSkc6hUGOdkXxTEAATJ8On38OK1Z47lMIeBi7lo6ugAYDK4HvsGvvhwHd7FWUNIm3yiczcCGWtguOdkWJnbx5YdYsmwX16FHPffoCJbBbaguIOtsJAh7DRr8u7uinKEqaxFvlswN4PnokA2OMwXpd7PCNWEq65qGHbOrtJ56A8PCY7QZ4GZiJTU1YDxvxoB4298+92K85JqUEVhTF13jl52OMaQR8ifW+mIP1Ly8ItAFKA01EZFUyyOkT1M8nFRERAY0aQZUqMGSId2P3ADWBX4HCySCboig+wWd+PiKyEmiK3WLrjz0aHoC1eGuamhWPksoICIAZM2DyZPjmG+/GlgW6YVdHiqKkSbz28xGRlQ5NFgIUA0JEpLKIfO1z6ZT0TYECVgE99RQc9zJIW39gMzYagqIoaY5ERzgQkcsiclRELvtSIOUWo04dm4K7XTu4GT2EQRwEY4OL9kLD2SpKGiRe5WOMuWmMqez4OcLxPrbi4fRYUeLhzTft6zvveDeuGRCKplZQlDRIQqJavw387fLzrROJVEkZMmSw5tcVK1pLuLCwhI0zwCfAA8CTWPNsRVHSBBrVWkk9rFoFTz8NP/0E+fMnfNwQ4CdgUXIJpihKYvCZtZsxZr8xxmMabWNMeWPM/sQIqCgANGgAHTtChw7WFDuh9AN2AcuTSzBFUXyNt8nkSgJZYmkLRDc+lKTy22/w/fcQGgrFi8ffP3duWLjQGv0/C9TBRkJQFCVV463ygdjPfCoBZxMviqIAZ87Arl1QqRIMHQo1a8bd33k+1ACoCLyPjQGnKEqqJl7lY4zpA/RxvBVgmTEmunFrEDbwyRe+FU+5JSla1DqfPvmkPf/Jmzdh4z4CKgDtsfE2FEVJtSTkzGc/1pXvW6x90TaX986yAKuguiaPmP7n6lV7Fn71qr8luUVo3Njm/unUKeHnP0WxKReeR20yFSWVE+/KR0SWAEsAbPxQhojILWdYMGkS7N4Nn30GvXr5W5pbhHfesabXH35ocwElhBeBqdivQ62TTzRFUZKGt7HdOt+KiufqVRsFZtUqm45GVz8pRKZM8MUXMHw4bN2awDHYJHN9sBEHFUVJlXhrav2qMWZULG2fGGP6+Uas1MWkSTb8WI4c1gr4s89SXoZbdtuvRAkYP96G37l0KWFjHsKmNXw7OQVTFCUpeBvbrTM2kL0ndjja0xXOVc8zz9j3Xbr4Z/Xjuu13y9GiBVSrFhWGJyEMB6YAu5NJJkVRkoS3yqc48GcsbftJh34+zlXPZ59Z15OHHwYRaNYMRo+GJUtg+3Y4ccI7v0hv0G0/YORIuwW3eXPC+hfAJnbviRofKEoqxFs/n8tAkVjaigLXkiZO6sL50F+/3iqcevXg779h/354+227I7Rypa37+284fx4KF7aWwrGVggVtKDNv8LTtd8sZPeTJA6NG2b3HHTsgMDD+Mc8Ck7EZUTskr3iKoniHt5lMFwGlgCoics2lPguwFTgsIs19LqWP8Da226efgjGeH/Se2q5ehaNHo5SRp/LffzaNjSfFVKyYfS1UyJ61O+esVcsqwMDAmO/THWFhsG5d7O1t2tgl6HvvJaz/90BLbDrunL4RUVGUhBFXbDdvlc+92BRe/2K/Tx7FroTaA3mAGiLyS5IlTia8UT7xPeQTqwSuX7d50+JSUCdOWL/KokXh2jXIlcu6vdx3H1StCtOmxa4U0zzxKZMTJ+Dee+HLL20UhPj6A3QHMgMeTWUURUku4lI+Xm27icgvxpjawAfAq9gzowhgI/BoalY83uLc6opNsQQGJm4LLHNmu11XIo7TsfBw+Ocf+Osva+DQvr1VWEOGWIf/0FDbnj27zcVWJLaN0LRI7tzxp1Rw9qlY0W7Hxce7wF1Yc5j7kyyhoig+INEpFYwxQUAu4IyIXPGpVMlEQlc+CV3VJPcWmKetvWvXrAIaNsyePR07BiEhUKOGDYNWowbcdRcEJDpHbUyuXoWePWHMmFSy1ScCzZvD/ffDoEEJGzMZGA9sIQn5exVF8QafpVRwRUSuiMixtKJ4vCG+VY8T19WPr4lu4u0kSxZrdfzFFxAUBEeOwIoVULu2DQbdqpVdDDRpAu++Cxs2wJUk/oZSnZm3MTB2rNWGv8Zm+R+NTkAGYFIyyqUoSoJJ1MrHcfZTBptGwQ0Rme4DuZKFhK58nnkG9uyJOvSPixs3oGxZ3z+Y4zJ2iK/PiROwaRNs3Ghfd+2Ce+6JWhnVqAH58iVMDufqbtUqm24nVRk6TJ5s7d2//x4yJmAH+Rds9OvdQAJjlSqKknh8aXCQE5uyq6qzyvEaOYmIeGVIbIwpho1HXN8x3zdAbxE5nMDxd2J92WsDwcBhYIyIjIzeN61kMvX1tt/ly/DDD1YZbdxoI9UUKGCVkVMhlS5tFVl0XBVcQhRiiiICjRrZ85/XX0/YmN7YsDu6AlKUZMeXymcMNl3XM8B3WCPWc8DTQDXgcRHZ7sV8WbHfR68BA7BKbCiQFbhHROKMp2KMqQSsAdZhd/XPYYPpZxORD6P3TyvKx5uHfGIUws2bdjXkXB1t3GgVmVMR1axpLesiItKAmfehQ9bqbcMGuPPO+PufB+4E5gHVk1k2RbnF8aXy+QubqmsWcAN4wKlsjDFjgWARecqL+V4EPgTKiMg+R10pbBSFVzwpEJexAdjkyX+ISMuEXC+tKB9/bPsdPmyVkVMh7dtn/Y1uvx1efNGeM+XMmQpXP2DPf6ZPt4InxIN3NjAMmxwkMekUFUVJEL5UPpeBBiKy0fFzYxFZ52irD3whIgmwfY2c71sgUERqRKtfDyAiteIYWwebS+ghEfkuIddLK8onNXDihF3ltGxpj1R+/NGaeNeoYc9/1q5NRSbeERFQt66NedSnT/z9BagLNMemYFAUJVnwpbXbP0T5iR/CbrU5CfVeNO7Crl6isxsoF89YZ37lQGPMVmPMDWPMSUd07aBEyKK4MG8ePP+8DSSwZo2NzDB6tFU4GTLYFdE999g+8+fDyZN+FDYgwJrkvfuuXbLFhwFGYzd4jyezbIqieMRb5bORKGODGcBAY8x4Y8xobBzhr72cLzdwxkP9aawPUVwUdrzOAVZhDRb+B3QBPvc04NSpU1SqVCmyTJgwwUtxbw08mXlnzgzVq9tz/Z9/topn9GgoXhymToU77oBy5aBHD2sGftyHD/UEpZO4/Xbo398KnZAIr3di/1ISmKNOUZSEMWHChMhnLHHYlXq77XY7UFhEvjPGZALeB9piDQRWAs+LyH9ezHcd+FBEXotWPxR4TURi3ZE3xkzApu0eJSIvuNS/6pCrnIj87jpGt90SRmLMvG/ehF9+scYI69fb8/+8ee3WnbMUK5Z4eWbMsL5XcZ413bwJDz5oQ0L07Bn/xJewa+/JWDMaRVF8is/OfHyNMeYEsFhEukerHwO0EZFYvVGMMe8BrwHNRGSZS/19wE9AOxFxWwGp8okfX5l5R0TAzp3uyigkxF0ZlSzp2bzb03US7Ge0Z49VQD/+aC8QH4uB17E2l5nj764oSsLxyZmPMSazMWaRMeYh34nGbux3z+iUw8Yhjm9sXCRTdp30ja+iOwQE2PifL7wACxbYM6Fly+CBB2xEhmrVbHy7Dh3sNffts247scmT4CyyZctC377QrZvnCaPTHLgN62mmKEqK4e222wXgEaeFW5IvbkxvbJDSO0Rkv6OuJNbU+jURGRHH2DzYqNqTROQ5l/rXsaEkSzvNt53oyid+UsrMWwT27o1aGTlzJrmujEqUsP6jXvsZhYfb8N89esSMT+SJ/UBl7Hq5uPefRVEUz/jS1HolsE5E3veRYMHYDY8rRDmZDgFCsE6mFx39SgB/AW+LyNsu4wcCb2INDdYAlYCBwBwR6RT9eqp8Ui8iNlCqqzI6fdqadz/zDDz0kA2YOmZMAv2Mdu60Ib937EiYTfjb2ETwC5P8URRFceCzlApAX2CxMeYidrf8ONGSFItIgre7ROSSw1/nI6z1nMH67vR2Kh4HBhsWMvo24dvABWyy5Jcd8gzHKjAlDWGMNVi7/fYoy7aqVa3dwJYt8PHH8O+/dtvut99svwcfhODgWCa8+2547jl49llYujT+w6VXgLuBr4DGPv1oiqJ4wNuVj1OxxDZI4rJQ8ze68kk7eLK4O3kSNm+22bQPHLCOsHfeaR1fq1e3r0WLukxy/brVVv36WQu4+Pga+zVmF6CeYoqSZHy57TYwvj4iMtgL2VIUVT5pg4Rmkf36a5vqwRkWaPNmyJo1KnJ39epwz43tZHiksbUDL1gw/ou3BsoDg9yvl6ryGSlKGiFJyscYcxqoJyI/GWMmA0NE5EAyyJnsqPJJGyQ2nYQI/PlnlDLatAmOHoUqef6kRtDP1Bj5GFWq2AywsXIEuA/YSmTMjgT7GSmK4kZSlc81IExEthhjbgLVROSHZJAz2VHlk/rxdTqJ//6DLeuvs6nrVDYVaMlPh/NFxqhzluLFox0JDceePK6Aq9dScT4jRUnlJFX57MWG1ZmOtSjrAfweW38R2ZB4UZMXVT6pn2RLJ7FlC7RqxfXtO/n5SF631VGGDO7K6N5ykOkBoA98ugjMH9CrCHx6FExZ6DURKJTED6ootwBJVT5PA+OJ3yHVYA0OvEoml5Ko8kn9JKufUd++8M8/MGtWZJXTxHvz5ihldPAgVMoOVf6BRQVg/fdQsBhcvQC1KsL6whDYFGv7GY8RnaLcyiTZ4MAYUwi4A1gLvEDcK59vEylnsqPK5xbn8mUbEfXDD236hVg4OwS2/ggfb4P9V+H4DRupp2ZNO8Xtt8Gbp8CURAOTKkocJNnPR0SOA8eNMdOA5WnV4EC5xcma1S6T2rWzXqs5c8bscxxyroOwZTCwJvx6GDJ8B7+E21x1GzbA0HdgfD6oGQ41r0LNJtatKGOqdTJQlNSHXwOLpjS68lEAe0h05QpMnhyzbQhwF3x6zHGeJMB87JrfscU2apQ1ZLjtP9i4FTZegr//tk6xzlTkVapAtmy+EVdNvZW0ii+TySlK2uf9922GvK89pJ9aDVfru+Qz6oGNoTEzqkvXrjY46mPvwoSsNuLCgQM2sd6lSzBwIBQoAJUqQe/eNtleUvIbTZpk/ZmSmipdUVITuvJRbk1Wr4YuXWwMOFfHn1rwaZtoVnQ/YKNf78amP8TF0m4usD7m9FevwrZt1oBh40Zr0JAzZ9TKqGZNaywREM/XP69TSihKKiLV5vNJaVT5KG506WJTtI4ZE1l1tSbUug7rN0R7yPcCwrF2nziUwkOwPgsEfhf/pSIirBXfxo1RCunsWRuFoWZNWypWjKlYXM3JvTItV5RUgCofB6p8FDfOnoXy5e0eW+3aAHzaFEwhhy+PK+ewWabmAdVt1addwJyAXstIFMePRymijRutcqpQIWp1dP/90KpVIlJKKEoqQZWPA1U+Sgy+/BJefBF+/RWCg3nmCdizAjJV8ND3JHAYqAgYuLEDyjaGzz730DcRXLwI338ftTrasMFu1TVrZpVOWBgsWqSrHyXt4MvAogFAgIiEu9Q1xIZiXCMiPydV2ORElY/ikQ4dIG9e+MiRzvQD4BAwEneTHAEaAXWwMeBKkmx+PlevWmvwkSNtRvB16+yKJ1cuOH8ehg+H+vWhcOHkub6i+AJfKp85wDURecrx/lnAuWF+A2giIt8kUd5kQ5WP4pH//rOOOvPn20MYAUYAK7ApFpoBmbB/4ROA3tj8P0NJtggHns53IiJg1y54+20bQPXvv63ODAuzu4a1akEhDfujpCJ8aWpdFZtuy0k/YBKQA5sDsn+iJFQUf5Inj3Xeefpp6/9jsCuamcBvQD2gluP1DPAS8GvyiXP1qouptwsBATZAw8yZ9sznyBGYOxfKlYMvvrCZXsuWtfnz5syxkYR8IYszuZ+i+BJvVz5XgAYi8p0xJhTYC1QQkV+NMQ2Az0UkbzLJmmR05aPESZs2NkXq+/Fkib8GVADeBVr6XozEppS4edNajq9bB2vX2jOjggXdV0YFCngvi6aTUBKLL7fdTgDdRWSxMaYL8K6I5He01QWWiIiP/Lp9jyofJU5OnLBLiy+/tBlQ42I90B67MgrxnQi+TClx86a1o3Aqo+++s9tyrsoof/74r6E+Rkpi8eW222bgNWNMU+zOt+sWXCjwd6IkVJTUQIEC1ujg6adtCu64qAXUBeLN7esdkybZVUZ8D/nAQGsnEVfUgwwZ4L77oE8fWLoU/v3XbtmFhsL06XDHHXar7rnn7HHXqVOeZcmRI/5rKYq3eLvyKY1VOLcD+7EZTg862tYAh0SkczLI6RN05aPEiwg0b26f2oPjyQh/CmvnuRKb/dQHJGtKiWjcvAk7dtiV0bp1dmVUtKhdFVWvDh98YE2+1cdISSw+9/MxxuQRkf+i1d0N/CMip2IZ5ndU+Sjx0qqV9f7cvt1uwcUXHfR8G8jYC7YAqTaTVcIID49SRlOnwl9/2VVS7dpQrx7s3QtBQXr2oyScJKdUcJnI6efzn0udq59PqlU8ipIgTp+2WU+nTLEWcKtXx70MqVUbInrZsDs9U0zKZCFjRhsMtXx5mDfPbsP9/ruNwTpqlL0txljd/PDDULlywlZoiuIJb898ZgORcegdfj4rsFnvtxpj6vlQNkXxH506Qb58du8pLozAOOzZjw9Mm1MDzrOebNms3cWrr1odfOoUdOxoozA8/7y9Pc2aWcX0++92x1JREor6+SiKJ4yBCRNs1tPffou7711AF6BPSgiWvET3MTp27BhDhw4F7JbbBx/YCAubN1tH13bt4JdfoFEjKFbM6uxZs6zhoK/kUT+j9Im3yic/cBTA4edTCvhURC4AU4C7fSueoviREiVsOIGnn7an83HxJrAVWJUSgiUf0a3tXnvtNd566y2OHDkCuFvZ5csHbdvaMQcPWnPuypVhwQJrCHHPPdC3L6xcafMcJVYezWWUPvFW+ZwH8jh+DgP+FRGnr/dNQO1glPRF9+72K78z7ltsZAVGY899rqSAXMlA9FXPv//+y5IlS8idOzfvuzjeduliTbVdVyPGQOnSNuPqwoV2i27CBBuL7r33rLNr7drwzjvwww/x63JXeVatink9Je2jfj6KEhcBAfZr9/vvWxvouGiMNbl+NyUE8z3RVz0jRoygbt26lChRgi+++MLj6ic2Mma0acUHDLDm2cePQ79+1tfomWfsqunRR2HcONi3z/N5kfoZpW/Uz0dRXAkLs7bG0Rkzxn793rTJem/G1v8ocC+wESibjHImA64+Rjdu/MsPP5ShTJnPOHz4fXLmrMXNmxcpXXo0kHQfo+PH4ZtvbFm9GrJksVG669WDunUhONjdr0j9jNIm6ufjQJWPEi+xKZ+ICPt0bNgQXnkl7v6fYM1v1pJsUa+Tm9dff52zZ8/SsWNHevfuzbJlyyhbtiw7duygWLFiPr2WiLXpcCqiDRvsaufOO62lXY0aVuFoJte0hy/D6wAQXfE46namZsWjKEnCuf02fHj81m+9gIvA9JQQzPf8+++/TJgwgddffz2yLl++fHTp0sXt7MdXGGPD/Lz4og2r9/ffEBJifY4GDLBRj9q0gaxZYfJkPftJL3jlZOrEGHMvUAYPBgYikkb/5RQlHkqWhKFDrT3x5s32YMMTGbC+P00dJY/nbqmVHTt20KNHD4oXL86xY8ci619++WWefvrpZL/+9Ol2deNc4Zw6ZZXS4sXW8u3ee21b8+bWIFFJm3h75pMTWA5Uw6bccm4qRE4iIqk2yIhuuynxEtu2mxMRG+K5Th14/fW4+78AXMZ6wqVRtm7dSu/evdm6dWuKXC++s53//oNq1aBKFWvCXbQotGhhFdG999pVlJJ68Fl4HawdTx7gQeA7bDaTc8DTWIX0eBLkVBT/kzu3VShxce0avPWWtSmO6/xjKFAOa3xQ03cipmfii+qdJw+88IJVMlOm2AXokiXWcu7mTauEmjeHBx/U0D+pHhFJcAH+Ap7CbixEABVd2sYC072ZL6VLxYoVRVF8wsSJIvffL3L9etz95orIXSJyLSWE8j1btmyRKlWqpMi1rlwRqVzZvnrbLyJCZOdOkaFDRSpVEsmdW6R9e5F580QuXEiaTJ07xy+T4hlgm8TyPPbW4KAQsF9EbgJXcU+jtRBokjRVqChpBKezyrBhcfdrDRQHPkwJodI2ScllZIwNiNq/P/z4o02iV706TJwIhQtDkybW6dXb1OIaYSH58Hbb7R8gp+PnQ9ittnWO96G+EUlR0gDG2CfT/ffDI4/YAweP/YBPgcpAW2xAKsUjP/9s/YzmzYu/r9PPKDaKFIEePWw5dw5WrLDbc6++asc1b27PiuKawzXCQoMG9vuG+hj5Dm+Vz0ZscNEvgRnAQGNMSSAc6Ags9al0ipKaKVrUrnw6dbIxY2I7ZLgN6As8h/3P0UNxjyTX6iJHDnj8cVuuX7f2IUuWWIfW4OAog4WqVa1FvRNPERbUx8h3eLvtNhj42vHzcGw0qybAE1jF87zvRFOUNECnTnZf5914Yur0BQ5iN6cVv5E5s13FjB4NR47YCNyZM9sQfoULQ9eu1qz7zBn3OHee4tkpScMr5SMif4nId46fb4hIXxEpKiK5ReRJ8eB8qijpGmfqhdGj7b5RbGTG+v70xobnVfyOMdaRdcgQ2LnTWs6VK2f9iAsXtkaN8+bZFBIJiWeneEeiIhwYY4oZY6obY+pEL74WUFFSPUWK2EQ3HTvafZ3YeBBogE2/oKQ6brsN+vSBr7+2SqhXL5g/H4oXt9tvpUvr6seXeKV8jDG3GWO2YDcQvgO+cZTVLq+KcuvRoYONgDBkSNz9hgFfANtTQCYlUUyaZFM4de1qz4b+/NPalbzyio3A3aoVHDjgbynTPt4aHEzCGo72BvYAcXzNU5RbCGNg/HioUMGeYFes6LlfXqwCehabfC7VxgO5NXFauK1fb98fPXqUIkWK0Lu3jT33/ff21/vAAzZZXqdO1sE1ONiPQqdRvN12ewB4QURGichqEVkfvSSHkIqSJihUyCad69jRHhjERkds8rmxKSWYklBcfY0uXbpEqVKlGDlyJGC/XzhzFL35pk2cN2eONXrs0sVm20hEkoBbFm+Vz9/oakdRYueJJ+zhwODBsfcxWOODwcCx2LspKUv0TK5jx47FGMOoUaOIiIiI7NelC3z+OTRtCsuX2yDnd9xh68uUsYaPf2tazXjxVvm8C7xqjNFFpqJ4whibnnPyZOtqHxt3At2xG9hKqiD6queDDz6gZMmSBAYGsmDBgsh+0S3fChWy50G//WYNEg4dsltyjRrZlZEaKHjGW1PrGcB64KAxZpkxZnq0Mi15xFSUNESBAjBypN1+i+vJ0x9reLAipQRT4uLnn+2KJiwM7rlnLBERtTh6NIhr17rTqdNgatWKICzMts+eDT/95D7euS03fjwcPWoV2aRJ1hiyZ0/7XcSbbbmrV63hQ3pVXt5au3UCXseG2LkfazwavSiK8thj1l534MDY+wRh3bR7YVMvKH7ls8/suc3y5Ze4dOkD1qx5k9BQmDu3JnfdlZXnnlvAunU2QsKmTXH7/AQFwZNP2sysP/9sV0ePPw53322t8hMSYy69x5VLTISDRUA+ESkiIqWilduSQUZFSXsYA2PGwLRpEFcunEZYM553UkowJT7Gjh3LQw89RPny5QEwxjBo0CAGDx7sdvaTUIoXtwYK+/bZP4ndu21MuWbNYNEiz65hrnHl0qtvkbfKJw8wRkTOJoMsipK+yJ8fRo2y9rhXrsTe72NgAhBPdm4lZRg1ahRvvfWWW93DDz9McHAwa9asSfS8xsBDD9k8RH//bf2FPv7YWsv17g2//BLV11NcufSGt8pnI/aoVFGUhNCmjfX9eTOOsAaFgEFY3x/vv1grPmbr1q2Rqx4nxhhWrFjBgw/65mQhWzb7nWT9etiyBUJCbHD0++6DESOsgkrvceW8VT4vAl2NMe2MMXmMMQHRS3IIqShpmk8/tREsN22Kvc+zwBVATXb8TqFChTzW586dmyxZsvj8erffbgNjHDhg48rNnm0t59q3t6kgsmRJn6sfb5XF78DdwHTgJHAjWlEfIEWJTt68drO/c2e4HItlQQZgPPAa8G8KyqakGjJkgJo17euBA1C/Prz+ul0N5c6d/lY/3obXeRtQH15F8ZaWLW2Uyv79bRQET9yPTU7yCjA5BWVTUg3Os56CBW2ah27dbIqHIUNsCohnn7UroAzpICyTV8pHRAb5WgBjTDHgI6A+1vf7G6C3iBz2cp7XgPeATSJS09dyKkqS+eQTa2vbsqU9efbEEKAcsAGIpYuSPokeVw6skcIjj9hoCl9+ac23N260IX7atYs9f2FawK9nNMaYrMAaoCw24lUHoDSw1psoCsaY24AB2K1ARUmd5MkDY8daz8FLlzz3CQFGYs+ADmH3Gh4CajlehwDHU0RaJYVxjbAQHacSevddaNLEbsHdcYd1aI0rjGBqxt8GAl2xSYZbiMhiEVkCNANKYIOPJJSxwCzsmZSipF6aN4fq1e1mfmy0wG5u1wHKA99i44p8C9wFtAc+IEU2wKtWrcrWuPyUFJ8QPa5cbHTtat3GvvoKZs6ExYshNNRa9MdlzZ8a8bfyaQZsFZF9zgoROQBsAponZAJjzJPY3fI4/psVJRUxciQsXGhd5T3xIVAZOAvcCzi3VjIBrbBZsw4BIxJ2ubVr11KzZk2CgoLInTs3HTp04MSJEzH67d69m1atWlG4cGGCg4O56667+OCDDwgPD4/sc/PmTV599VXy589PsWLF+Pjjj2PMM3fuXAoXLsz585qyNaHEtepxxTWuXI0a1hpu4UL45htrNTdiROyL6tSGtwYHvuYuYImH+t1Am/gGG2NyYc+LXhGR08YYH4unKMlArlx2v+Tpp+HXX63Th5Pj2Fhvq4HhwHPAV9jTUCcB2K25+kA7rJ9QLHz33Xc0aNCAhg0bsmDBAv777z8GDBhA3bp12b59e6Tp8LFjxwgLC6NIkSJ8/PHH5M2bl2+//ZZXXnmFU6dOMWzYMACmTZvG+PHjGTduHGfPnuW5556jQoUKhIWFAXDx4kVeeuklRowYQfbs2X1yu24Ffv4Z9uyxabvj48YNGyHByQMP2KR3v/wCQ4fC//5nnVZ79YJU/SsQEb8VrGn2+x7qhwLhCRg/CZtR1TjerwM2xta/ePHiUrFixcgyfvx4URS/0amTSM+e7nVvi8gCx8/XRaS8iMyJZfx8ERkS9yXq1q0rt99+u9y4cSOy7scffxRARo8eHVk3fvx4AeSPP/5wG9+2bVspWLBg5PvWrVtLjx49It83atRIXnnllcj3ffv2lbp168YtVBri3nvvlZ9//tnfYnjF7t0iTz4pkjevyKBBIqdPp+z1x48fH/mMBQ5KLM9jf2+7JRpjzIPAU0APkYTFis2XLx/btm2LLN26dUteIRUlLj76CJYuhW+/japbDTzi+DkTNu9PH+CMh/HNiDdx/datW6lfvz4ZM0ZtclSqVIk8efKwaNGiyLrrjgBj0VcrOXPmdItndv36dYKCgiLfZ82alasO55Ndu3Yxbtw4Ro8eHbdQSrJSrpz1ad682aZ3CA2FN96Af1PIf6xbt26Rz1ji8Frzt/I5A+TyUJ8bz/9urowHPgP+NsbkNMbkxG4jZnC8970rsqL4kpw5YcIEe8p84YKtM0Sd8QDUAB4FehDTwCABZrYZMmQgc+bMMeqzZMnCrl27It+3adOGvHnz8txzz3HgwAHOnz/PokWLmDFjBn379o3sV6VKFRYsWMC+ffvYvn07q1atomrVqgD07NmT3r17U6ZMmQR8eCW5KV3appXavh1On7bWcS+/nLCI2ilCbEsiidraWuNF+Ta++TzMHWObDLt9tj6esRJP6R19TMWKFZNrpakoieeZZ0S6d7c/Pyh2u82VyyJyp4jMjFZ/XUQeinvqBx54QCpXruxWd/DgQTHGSObMmd3q//zzTylXrlzk/5AxRgYPHuzW5/z581KrVq3IPo8//rjcvHlTpk2bJqVKlZLLly8n4AOnHdLitltsHDki8vzzIrly2dcjR+Luf+WKSOfO9jWxANskCdtuAdjvY85SFggDSmIzkpR0vC+D+7FoQlgKVHX46QBgjCmJ/b63NJ6xtT2UX4Bdjp/neymLoviHESOs2dKqVdaIYFm09iBgJnb77ZBL/VJs/zh48cUX+eGHHxgwYAAnT55kz549dOjQgYCAAAICov79T506RatWrQgODmb+/PmsXbuWAQMGMHTo0EhjA4CQkBDWrVvHwYMHOXbsGLNnz+b8+fP069ePUaNGkSVLFvr370/RokUpUqQI/fv3T1QaAsX3FC1q/Zx/+w0yZ7bZVp99Fg4e9Nw/2fMJxaaVPBWsB8IhoEq0+irAQaC5l/MFA/uAnVjT6mZYBbIfyObSrwQQDrwVz3zriMPgQFc+Sqpl1SqR4sVF9pwTqSMiNz30eV9EaolIuKO9jogcj3/qAQMGSGBgYORq5vHHH5dHHnlESpUqFdnn5ZdflpCQEDkd7XT6jTfekEyZMsmpU6dinb9nz57SokULERGZMGGCFC9eXPbv3y/79++XYsWKyaRJk+IXMpWSnlY+0Tl5UuT110Vy57YrnL17o9quXBGpXFnk7Fn7mtjVD0lc+bgyBHhTRL6PpsC+xwaFH+rNZCJyCetKtxeYgXUUPQDUEZGLLl0NNvSiv8+oFCV5qF8fGjWC4S/Bw9j48dEXDC9jN7uGAy9g+xWMf+ohQ4bw77//8uuvv3L8+HFmz57Nn3/+Sc2aUVGodu7cSWhoKLlyuR/BVq5cmRs3brBv377o0wKwfft2ZsyYwciRIwFYuXIlrVu3plSpUpQqVYo2bdqwcuXKhN0DJUXJl89GTNi3D0qUsL7P7drZlVFK5BPy9mFeGjgVS9tJINRbAUTksIg8KiLZRSRERFqIyMFofQ6KiJF4YsuJSJhoXDclrfLBB9Zb8K6Vdq1fH1iAjRcPVhk9AbyFNTbo63kaTwQHB3P33XdToEABVq5cyZ49e3j22Wcj2wsWLMi+ffs4c8bdzuf77+33zCJFisSYMyIigp49ezJgwACKFy8eWX/Jxcvx4sWLzl0JJZWSK5fN9v7XX1C+PISFwVtvWf8hSL58Qt46mR7Ahr1Z4aGtO3brTVGUxBASYr9iduoEO3dCu5zWk+0Tlz71sREQxgJXsedBcfDzzz+zYsUK7r//fgA2btzI8OHDeeWVV6hevXpkv2effZZZs2bRoEED+vXrR548eVi3bh0ffPABLVu2pFixYjHmnjhxIpcuXaJPnz6RdfXq1ePVV1+NTLr2+eefM2JEAkMxKH4le3Yb9SlzZpt6qkULG1l74MCo1U+vXj68YGz7cZ4K8DhwE3uoPwhrADrI8T4caOvNfCld9MxHSfW0bClSuLBIgQIitWp5Lg/VEsn3rUiRebZ/HOzatUtq1KghOXLkkMDAQLnvvvtk8uTJHvtu2bJFHn74YSlYsKBkzZpVypUrJ0OGDPFowXbq1CnJnTu3rF+/3q0+PDxc+vbtK/nz55d8+fLJyy+/LOHh4Ym4EamD9Hzm4wnnWc+VKyKXL4scPhyz3huI48zHGRkgwRhj6gGDgUrYxf8N4EdgoIh8G9dYf1OpUiVxOD4pSuokLMzGzr/nHhstskkTz/3OAPcAufrBr8NTUMBbiwoVKjB16lQqVKjgb1FShE8/tRG0Pa1w4mqLDWPMdhGp5KnN6wN8EflGRGpgF/wFgSARqZnaFY+ipBmyZbPegd27Q7QzmEhyAVOBP16B/1JQNiXdEj2y9oQJE+jloml8ffaTKOsxY8y9QGugEdDeGPOUs/hGLEW5xQkLg1at4MUXY+9TF8i/Brqh+YWVJOMaWTsiIoIhQ4YwceJETp8+DbhH1PYFXikfR9iaTcDPwGzsd6+pwBSXoiiKL3jvPRuga2kc/talPoM/gWnJI8LWrVsjw+co6Zfoq54lS5YQEBDAbbfdxkcuad99ufrxduXzLpAHeBDre9MS66czC+sYWjnpIimKAkBwMEyZAj16wMlYkvQGXLf/ff2w/4GKkgiir3oGDRpE06ZNue+++xgzZkyyrH68NbVuiDU2cKY2/FtEtgPrjDFjsa5xuvWmKL7iwQehY0do396G4MmQIWafu7GpFDtgM576O0uXkuZwzSd06tRiDh/OyOXLd3Pp0s9kztyK++77iFKlhgAx8wklFm9XPoWA/SJyE+tlEOLSthCIxTRHUZRE8/bbcP26zRQWG72BQGBY7F0UJTY++8z69qxZE0HGjIP5/PNBvPyy4ZFHYPPmN7h4cQwLF55m3TrbzxcrH2+Vzz9ATsfPh4BqLm1eRzdQFCUBZMwIs2fb7KfffOO5TwD23Gck1vFBURLB4sWLyZgxI02bNo2sK1WqFK1ateLDDz/06bW8VT4bAefp4wxgoDFmvDFmNDbi1Ne+FC690bVrV4wxbh7htwqDBg1izZo1/hYj7VKokM0Q1qEDHD3quU9RYBTQHrjkuYuixIUxhlGjRmGMe4KCt956iwIFCvj0Wt4qn8FEKZjhwGjsVtsT2ADvz/tOtPTFlStXmDt3LmBDjoSHh/tZopRl8ODBqnySSu3a1sPv8cchtr+ftlizn5dTUjAlvdCyZUu3sEtOihUrxvPP+/bx7pXyEZG/ROQ7x883RKSviBQVkdwi8qSIqLtbLCxevJjz58/TuHFjTp486dNIvyISmQb5VuLatWv+FiHleeMN64Tav3/sfT7FRl9cnlJCKYr3aIqCFGLatGnkypWLqVOnEhQUxLRpnh0zFi5cSNWqVcmaNSs5c+akTZs2HD582K1PyZIlad++PZMnT6Zs2bJkzpyZ5cvtk2blypVUq1aNoKAgcuTIQYsWLfjjjz9iXGfRokXUqFGDbNmykT17dipXrsxSF3+STz/9lGrVqpE7d25y5sxJ1apVI6/h5ODBgxhjGDNmDC+99BL58+cna9asNG3alIMuGaqcS/h33nkHYwzGGAYNGgRAWFgYYWFhMeQrWbIknTp1inw/depUjDFs2LCBNm3akDNnTqpUqQJAeHg47733HmXLliVLliwULlyYvn37ctXXYXhTgty5rYNpbKVOHTh/HkaOhLvvtv2jkwOYDnTFxppXlFSI10aZxphC2GDutYDcwGlgLfChiKSW7OCpimPHjvHNN9/QtWtX8uXLR4sWLVi4cCFnzpxxy58ybtw4evToQefOnXnrrbe4cOECgwYNolatWvz666+EhEQZF65du5YdO3YwcOBA8ufPT8mSJVm5ciVNmjShTp06zJkzh4sXL/LWW29Rs2ZNduzYERkWf9SoUbzwwgu0aNGCadOmkS1bNn766Sc3hXHw4EG6dOlCyZIlCQ8PZ9myZTRt2pQVK1bQqFEjt8/33nvvUaFCBaZMmcLJkyd54403aNCgAbt37yZTpkxs2bKFatWq0alTJ7p37w5A0aJFE3Uv27VrxxNPPMH8+fMjty7bt2/PsmXLePXVV6levTq///47b775JgcPHmTBggWJuo7fWLgwYf22bLFhh5dFT3vq4CGs00MXYAne5xhWlOQmtoijngpwB3ACuI5VOLMdr9ewlnClvZkvpYu/oloPGzZMANm8ebOIiKxcuVIAGTt2bGSfCxcuSPbs2aVz585uY/fv3y+ZMmWSjz76KLKuRIkSEhQUJMePu6exrFixooSGhsqNGzfcxmfMmFH69OkjIiLnzp2TbNmySct4oiG7cvPmTblx44bUr19fmjVrFll/4MABAeTOO++UmzejUm9u3LhRALcMloD0798/xty1atWSWrVqxagvUaKEdOzYMfL9lClTBJDevXu79duwYYMAMm3aNLf6mTNnCpC+IxJ/+KFIpUoiV696br8mIhVEZELiL7FlyxapUqVK4idI49xqUa1dGTdunHTr1i1Jc+DDTKbDgPPAHSJSW0SeEJHaDqV0DvUy8Mi0adMoXbo01apZy/R69epRuHBht623LVu2cP78edq1a0d4eHhkKVasGGXLlmXDhg1uc1atWpWCBaPSWF66dImffvqJtm3bkjFj1IK2VKlS1KhRg/Xr1wOwefNmLl68SLdu3eKUefv27TRt2pQCBQqQMWNGMmXKxOrVqz1u4bVu3ZqAgKg/pRo1alC0aFG2bNnixV1KGC1btnR7v3LlSjJnzkzr1q3d7luDBg0AYty3dEXv3lCsGPSNJatcZmz0gzewIXgUJRXhrfKpjU2jfdC1UkQOYfP61PaNWOmHbdu28dtvv9GqVSvOnj3L2bNnuXDhAq1atWLr1q3s3bsXgJOO8Cn16tUjU6ZMbmXnzp3895+7LUehQoXc3p85cwYRiVEPNkulMzyGc564tr2OHDlC3bp1OX36NKNGjWLz5s38+OOPNGrUyOM5iicTzAIFCnA0NpPgJBD98508eZLr168THBzsds/y588PEOO+pSuMsdGvV66EOXM89ymHzXzanqiMqIqSCvD2zCczcCGWtguOdsUF5+pm2LBhDBsWc2E4ffp0hg4dSp48eQB7sH7XXXfF6Od63gPEsMPPlSsXxhj++Sfmsds///xDbsfBdN68eQE4evQo5cuX9yjzypUrOXfuHHPnznVTUpcvX/bY/8SJEx7rEpIDJTAwkPPnz8eodyrL6ET/3Hny5CEwMJDvvvvOY//ChQvHK0OaJmdOmDsXGjaEChWgTJmYfXphLd+GYp0lFCUV4O3KZwfwvDHGbZyxT4SejnbFwfXr15k9ezZVqlRh7dq1MUqFChWYMWMGIkL16tUJCQlh3759VKpUKUYp4+mh4kJwcDAVK1Zk3rx53Lx5M7L+0KFDbN68OdKirHr16mTLlo0JEybEOpdTyWTKlCmybu/evWzatMlj//nz5xMRERH5ftOmTfz999+R24wAmTNn5sqVKzHGlihRgr1797qZim/YsIELF2L7juOOczV27tw5j/ct3SsfgPvvt6F32rQBT18QAoDJwHjA9zuhipIovF35vA18CfxujJkDHMcmlGsDlEZju7mxfPly/vvvP0aMGOHRnLh79+706NGDdevWUbt2bYYPH06vXr04deoUDz/8MDly5ODo0aOsX7+esLAwnnzyyTivN2TIEJo0aULTpk3p2bMnFy9eZODAgeTIkYO+jnOBkJAQ3nvvPZ5//nkeffRR2rVrR0hICDt27CAwMJDnn3+eevXqkTFjRp566in69u3L8ePHGThwIMWLF3dTMk4uXLhAixYt6N69O6dOneL111+ndOnSPPVUVIzZcuXKsXz5cho1akSuXLkoXLgwhQsX5vHHH2fChAk8/fTTdOrUiQMHDvDhhx+SI0eOBN3jsLAwnnjiCVq3bs1LL71E5cqVCQgI4ODBg3z11VcMGzaMO+64I0FzpWm6dYPvvoPnnrNbcdEpDIzBBh/9GfeojIriD2KzRIitYBPIbQNuAhGO1x+Aht7OldIlpa3dmjdvLiEhIXLp0iWP7WfPnpWgoCA3q67ly5dLWFiYhISESFBQkISGhkrnzp1l9+7dkX1KlCgh7dq18zjnihUrpGrVqhIYGCjZs2eXZs2ayZ49e2L0mzdvnlSuXFkCAwMlJCREKleuLMuWLYtsnzNnjpQpU0ayZMki5cqVk9mzZ0vHjh2lRIkSkX2c1m6jR4+WPn36SN68eSUoKEgaN24s+/fvd7vexo0b5f7775csWbIIIAMHDoxsGzdunISGhkpgYKBUq1ZNtm3bFqu1259//hnjs9y8eVM+/vhjueeeeyRLliySPXt2ueeee6Rfv35y9uxZj/cpoRw4cEBeeeUVOXXqVJLmSREuXBC5806RyZNj79NZRJ5J+JRq7abWbkmBOKzdEv0gB7ICRYCsiZ0jpYu/TK3TK07lM3HiRH+L4nMOHDggXbt2ldy5c0v//v3dzNdTNbt3i+TNK/LLL57bz4vIbSKyKGHTqfJR5ZMU4lI+iY5wICKXReSoiHg+hVaUNMjBgwfp1q0bFStWJH/+/Ozdu5ehQ4e6ma+nasqVg48+suc/Hgw5CMGGBH4Wu2muKH4ijfxHKUrycvHiRRo2bMiOHTuoWrUqzz//PMHBwUz2dH6SQoSEhNC9e/cYFn7x0r49bNhgz4Fmz7Ym2a5UB7oBTwNfodEPFL8Qr/IxxkQAksD5RERUod0ilCxZ0rkFm+a5ceMG//33H8YYRISzZ89y6ZJ/8xJ4sg5MMJ98AtWqwZgxNhJ2dN4EamCNEDw0K0pykxBF8TYJVz6KkibJlSsXe/bs4aeffuLtt99m3rx5vPrqq3Tt2pWgoCB/i+c9gYE2J3L16lC5MjzwgHt7JmAmdhVUB7gz5UVUbm3iPfMRkUEiMjihJSWETu8MGjTI+60WICIigt69e1OoUCECAgJo0aKF74WLB2ek66lTp0bWderUKUGBRJ2Rq10DnKY0999/P4sXL2bZsmWsWbOG0NBQPvnkE48m5qme0FAYNw4eewzOnInZfgfwDjb6wa2XkUPxM14ZHBhjNhtjOhhjsiSXQAp06dIlUXHR5s+fz8iRI+nXrx+bNm3if//7XzJIFzeFChViy5YtNGmStl2+XJXQzz//zKlTp/wtUuJo1cpGv+7YETwp0G5YH6BBKSuWonhr7XYdmyn+mDHmQ2NM2WSQ6ZanaNGiVK1aNf6O0fj9998B6N27N9WqVfOLc2WWLFmoWrUq+fLlS/FrJwf3338/U6ZM8XkK4RRl2DA4dQo++CBmmwEmAVMAzxGKFCVZ8DaTaRg2VOE0bLaQ3caYdcaYtsaYTHEOVhJM9G03YwwDBgzgk08+oVSpUoSEhFCrVi12794d2adkyZKRCdoyZMjgtvV1+fJlXn31VUqVKkXmzJkpVaoU77zzjttW0tWrV+nTpw/ly5cnW7ZsFCxYkEceeYQ9e/a4yfbPP//QsWNHChcuTJYsWShUqBBNmzaNDIzqadvNyebNm3nggQcIDAykZMmSjBo1KkH3Y8KECdx7770EBgaSN29ennnmmVhjvykeyJzZBh798EMbBSE6BYAJ2P/ocyksm3LL4rWfj4jsEZGXsA6mnYAMwOfA38aY940xt/lWRAVg5syZLF++nJEjRzJlyhQOHz5M8+bNIxOqLVq0KDLz55YtWyK3vsLDw2nYsCGTJk3ixRdfZMWKFXTp0oUhQ4bQr1+/yPmvXbvGhQsXGDBgAMuXL2fs2LFcvXqVatWquQUr7dChA1u2bGH48OGsXr2aTz75hKJFi8YadNTJ+fPnadu2LR07dmTx4sWEhYXxwgsveFRSrrz22mv06tWLevXqsXTpUoYPH87KlSt5+OGH3WLYKfFQvDhMmQJPPAEnPaQ3fQRoCLyQ0oKlbnbs2JGgALlKIojN+zShBbgfWIcNtRMBhAPzgIJJndvXJa1EOBg4cKDYX40FkNDQULl+/Xpk3bx58wSQTZs2Rdb179/fbZyIyPTp0wWQ9evXu9UPHTpUMmXKJCdOnPAoQ3h4uFy6dEmyZcsmH374YWR9cHCwjBw5MlbZnVEPpkyZElnXsWNHAWT27NlufevVqyfFixeXiIgIEYkKoXPgwIHIuQICAmTw4MFu45zJ6hYtWhSrHEosvPGGSL16IuHhMdsuikhpEZkTVZUWIxwsX75cHnzwQQkODpaQkBCpWLGifPvttx77du/eXYBYw1X99ttv0rp1a8mTJ48EBgbKHXfcIR9//HFke3h4uLzyyiuSL18+KVq0qFvSRydz5syRQoUKyblz53zy+VKKVBnhwBgTZIx52hjzA/AjkB94EXt02QNrwDkrUdpQ8Uj9+vXdokzffffdABw+fDjOcStXrqREiRJUr149RrK1GzdusHXr1si+c+fOpUqVKuTMmZOMGTMSHBzMxYsX3RLIPfDAAwwfPpyRI0eyc+fOBPv5ZMiQgUcffdSt7vHHH+fw4cOx5v1ZvXo1ERERMRLsValShZCQkPSdKC65GDwYwsNhyJCYbcFY8+vngO1YJ4tngd+wabmHkOqjIowfP57mzZtTsWJFFi1axLx582jTpo3HlfmmTZuYOXMm2bNn9zjXtm3bqFKlCteuXWPSpEl89dVX9O3b123FPW3aNMaPH88nn3xC//79efnll1m3bl1k+8WLF3nppZcYMWJErNe5VfHKIdQYczfQHWiH/VNdArwqImtduk00xvyDXf0oPsKZj8dJlizW4NBTcjdXTp48yaFDh9wUlyvOZGvLli2L3BYbOHAgefPmJSAggMaNG7tdY86cOQwePJj//e9/kWbdzz77LAMGDHDLZhqdXLlyxZDBeYh/9OhRj6bYznOk0NDQOGVXvCBjRhv1oGJF6wPkyPgayQPAfUB97DlQHeBl4FtgGdYs+2GgL6kuMsLBgwfp3bs3w4cPp3fv3pH1DRs2jNH3xo0bdO/enf79+zN+/PgY7RERETz11FPUrVuXRYsWRdbXru2eL3PFihU8+eSTPP744wAsWbKEFStWREaxHzRoEGXLluWJJ57wwSdMX3gbjeAX4BjwMTBBRGL7HrQPzRySKsiTJw+lSpVi7ty5HttLliwJwBdffEFoaKjbGcyNGzdiHOznz5+f0aNHM3r0aP744w+mTZvGwIEDyZcvHz169IhVjjNnznDjxg03BeRMQlekSJFYZQdYtWoVuXLlirVd8ZKCBWHWLHv+s20buN7/EUAocB74G3B+J8gEtAJaYPc4RmCVUipi8uTJBAQE8Oyzz8bbd/jw4dy8eZOXX37Zo/JZt24dv//+u8c2V65fv+7mhJw1a9bIL2u7du1i3LhxbN++3ctPcmvg7bZba6CEWIfSWBfgIvK7iGhK7VRAo0aNOHLkCNmyZfOYbM2Z2fTy5csxgmfOmDEjzkP9MmXK8O6775IrVy527doVpxw3b95kwYIFbnVffPEFxYsXj1X51K9fn4CAAA4fPuxR9lKlSiXkFiieCAuD55+Htm3hhiO/9nFgBTAKu/32DvZrpCsBwEhHv1S2Bbdx40bKli3LF198we23307GjBkJDQ1l9OjRbv327dvH0KFDGTNmTKw7Ahs3bgTszkLVqlUjU7O/8MILbmGPqlSpwoIFC9i3bx/bt29n1apVkW4SPXv2pHfv3vEmgrxV8WrlIyILk0sQJXlo164dU6ZMoW7duvTt25d7772X69ev89dff7F06VIWL15M1qxZadSoEYsXL6ZPnz40bdqUbdu2MWrUKHLmzBk517lz56hXrx7t2rWjbNmyZMqUiSVLlnDmzBkaRN++iUZISAivvPIK//77L6VLl2b27Nl88803kVENPHH77bfz6quv8txzz/HHH39Qq1YtAgMDOXLkCKtXr6ZLly4xtkEUL3jtNdi4Ed54A4YPt/4+vbAK5nZgGNz9/t0Mfi9a4JIAbN7iz4ABKSxzHBw7doxjx47Rr18/3n33XW6//XbmzZvHc889R3h4OC+++CIAPXr0oFWrVnH+7Rw7dgyAtm3b8txzz/H++++zbds23nrrLY4cORK5Fff888+zatUqSpcuDdhzzLZt2zJ9+nT+/vtv+vfvn8yfOu2SqCCgxph7gTJAYPQ2EZmeVKEU35EpUya+/vpr3n//fSZMmMCBAwcIDg7m9ttvp0mTJmTOnBmArl27cuTIESZPnsz48eN54IEHWLZsGS1btoycKzAwkPvvv5+JEydy6NAhAgICKFOmDLNmzaJ58+ZxypE9e3a++OILXnzxRXbu3EmBAgUYOXIkHTt2jHPcu+++y5133hm51WeMoVixYtStWzfyH15JJAEBMGOGTcNdsyasbg6vubR3huAvg2m4tqE1xXalGfAJqUr5REREcOHCBaZOnUqrVq0AqFOnDgcPHuS9997jhRdeYNasWfz4449uRjSxzQXQvn173n77bcBmzb158yavvfYav//+O3feeSchISGsW7eOQ4cOkTlzZgoVKsTZs2fp168fkydPJkuWLPTv359p06YhInTq1IkhQ4bEeT56yxCbGZynAuQENmGzlzozmUa4vL/pzXwpXdKKqbWipChbtojkyydS+XLMtlNik8995mHcQ8ksl5dUrVpVADl//rxb/YcffiiAHDp0SPLlyyfvvvuunDlzJrIUK1ZMHnvsMTlz5kykO8Nrr70mgCxdutRtrp9++kkAmTVrVqxy9OzZU1q0aCEiIhMmTJDixYvL/v37Zf/+/VKsWDGZNGmSjz958pDaTK3fBfJgDS8N0BJrDzML2A9UTpImVBQl5ala1W69/fY7XLzm3pYXm/PndeAbl/obKSdeQrnrrrvibD927BinTp3ijTfeIFeuXJHlyJEjzJ07l1y5crF8+fIEzRXbymX79u3MmDGDkSNHAtbVoXXr1pQqVYpSpUrRpk0bVq5cmYhPl/7wVvk0xCogp3PI3yKyTkSewv5pvuhL4RRFSSFefBGK/QGtp8VsKwPMBZ4EnBGdlmLNsVMRzi3ir7/+2q1+5cqVFC1alAoVKrB27doYpUCBAtSrV4+1a9dSs2ZNAB5++GGyZMnicS6ASpUqxbh+REQEPXv2ZMCAARQvXjyy3jUv1MWLF9NNDqyk4u2ZTyFgv4jcNMZcxSbldbIQ+MJnkimKknIYA4uawH274PMv4MnH3dtrAR8CTYDN2CR0qcyNvHHjxtSuXZvu3bvz77//cttttzFv3jxWrVrFlClTCAwMjPS/cSUwMJACBQq4teXJk4fXX3+dIUOGkD17durUqcO2bdt4++236dixo0ffs4kTJ3Lp0iX69OkTWVevXj1effVVHnzwQQA+//xzRowY4fPPnhbxVvn8gz33ATgEVMOG1gHrHaAoSlrl9U6Q+yFonwk+qgLBHpLoSUco2QKKzIeeK2Bh6jGANcawePFiXn/9dQYOHMiZM2coW7Yss2bN4sknn/R6vrfeeouQkBDGjBnDBx98QKFChejXrx9vvvlmjL7//vsvb7zxBosWLXIz3+7WrRt//fUXL730EiLCs88+yzPPPJOkz5leMN4sAY0x04HDIjLAGPMG8BY2wnU40BFYKiLe/5ZTiEqVKsm2bdv8LYaipE7CwmDtOnhsq/XjmXAftMliHUxvYLfaxmAjXxcFztSB9Wv8Jq6SvIwfP56ffvopXkfbuDDGbBeRmHuUeL/yGYyN3wYwHGt80BbIiv3TfD6xQiqKkgowwNwq0Pp5GH4MxrckMo5OfazzaR7s6e9f8UcSUJTY8Fb5nAOOAojIDWyEp76+FkpRFD9iDEwfBpUrQ8ez8PTTMfssBIpUhdFYx1RF8ZJ4rd2MMRmMMYOMMWeAE8B5Y8wCY0zOZJdOURT/EBwM8+fDq6/CL7/EbM8F3P0qDAW+TGnhUo4KFSqwY8cOf4uRLkmIqfWz2LOdn4EPsJGsmwMfJaNciqL4mzvvhI8/hpYt4biHQG5Bx2ER0Bn7dFAUL0jItltXYKKIdHdWGGO6A58aY7qLyPVkk05RFP/Srh0cPGhTL6xfD9FSe1AVGIsNt7OFqCjYihIPCVn53EbM3DxzsOmzS/hcIkVRUhdvvAENG0LjxnDxYsz21tj0202wqRgUJQEkRPlkI+af1AXHawiKoqRvjLFRr8uXhxYtwFMCw5exXn9tsY4XihIPCQ2vU8QYc5uzYFdDMeodbYqipDeMgfHjIVcum4QuPJqGMcCngGDTcGsEGSUeEqp85gN/upQ9jvrF0er/9FYAY0wxY8x8Y8w5Y8x5Y8xCY0zxBIyrZIyZYIzZY4y5bIw5bIyZZYzRDGOKkhxkyAAzZ8Lly9ClS8z2jNgYcJuxmU4VJQ4SYnDQObkubozJCqwBrmEjJAjWeHOtMeYeEbkUx/DHgbuwWUV2A0WAN4FtxpgKInIkueRWlFuWLFlsSJ0GDeDQIRCxqyIn2YHlQHWgFPCof8RUUj/xKh8R8RDm1md0xW7hlRGRfQDGmF+xK6ju2FCGsTFMRE65VhhjNgEHHPO+lSwSK0p6JXduG2InIRgDFy7A22/DwIHubcWw8U4aYL8SVvWplEoK0b179/g7JYFEZTL1Ic2ArU7FAyAiBxxKpDlxKJ/oisdRd8gYcwr7J68oijd4GyT0xAmbATVnTpuSwZX7gKnYjF+biDolVhQH/s7lehewy0P9bqCct5MZY+4E8gO/J1EuRVHio0AB+OYbGDECpnnYIGmCTbPdBDiTwrIpqR5/r3xy4/nP8jQ2gEeCMcZkBMYBp4DPPPU5deqUWxKobt260a1bN28uoyiKKyVKwNdfQ506kCOHNcV2pRewD2gFfA1kTnEJlRRmwoQJTJgwwfk2b2z9/K18fMmn2GPOJiLi8XtWvnz50JQKiuJj7rwTvvwSHn4YQkKgbl339g+whgddsVtxJsYMSjrC9Uu9Mebf2Pr5e9vtDJ5XOLGtiDxijHkf6AY8LSKrfCSboigJpWJFG4j0iSfg++/d2zJgs57+Bgzxg2xKqsTfymc39twnOuWwf6rxYozpD7wKvCAiM3wom6Io3vDQQzBlCjRrBjt3urcFA8uAydicQEqqZOPGjdx5553kzZuXgQ4rxq+++ooyZcpQtGhRhg0b5rNr+XvbbSnwgTHmNhHZD2CMKQnUAF6Lb7Ax5gWsX1B/Efk0OQVVFCUBNGliI2E3agQbNsDtt0e1FcT6ANXGmmPX8ouESixcu3aNHj16MGXKFHLmzMmjjz7K9evXmTx5MpMmTeLee++lS5cuVKxYkXr16iX5ev5WPhOxwTiWGGMGYJ1MhwBHgMjcrcaYEsBfwNsi8raj7nHgY2AlsMYY4+pNcF5EErRyUhTFxzzxBJw7B/Xrw8aNULhwVNtdwOfAY8AGoIx/RFRisnPnTqpXr07VqvZROnfuXO69914+/fRTHnnkEQDefPNNli5dmvaVj4hcMsbUweYGmoE9ivwW6C0iruFzDXbn2HWbsJGjvpGjuLIeCEsmsRVFiY9nn4UzZ6wC2rAB8uSJaqsHvAc0BrYC+fwjouJOREQEGTNGqYTChQuTNWtWzp+PiiudIUMGIiIifHI9f698EJHDxBOEQ0QOEs1GRkQ6AZ2SSy5FUZLIa69ZBfTww/Dtt9YSzsnT2L2M5tivm0H+EVGJonz58mzYsIFffvmFggUL0rlzZ9q1a8f777/PPffcQ/Xq1Rk2bJjP3FP8bXCgKEp6xRgYNgwqVIDmzWOmYhiCzQjWEfDNl2klCWTNmpWRI0fSsmVLihcvTpEiRfj000+ZNGkSzzzzDIUKFSI0NJQmTZr45HpG5NaJfV6pUiVRPx9FSWFu3oQnn4Rr16w5tsvWDlex23A1gff9JF8cVKhQgalTp1KhQgV/i5KiiAjGuDtkRUREEBDg3XrFGLNdRCp5atOVj6IoyUuGDDBjhlU+Tz8NrmcGgdjELAuw5kdKqiC64gG8VjzxocpHUZTkJ3NmWLAA9u+H3r1tKgYneYGvsAlR1EX8lkGVj6IoKUPWrDYMz3ffwaBB7m2lsSkr2wM7gePA28BDWH+gh7BnRMdTUF4lWVHloyhKypEzJ6xcCbNnW2dUV2piPfceAtoA5bGWcOsdr3dhldMHaJruFGD8+PHJmtNHlY+iKCmLMxXDRx/ZcDyuHAPKApexyegyOeozYSNjrwYOoWm60wGqfBRFSXmKF7epGN54IyqJ3XFgBbARqAA8CdyMNi4AGOnop1twaRpVPoqi+IeyZWH5chsNYfVqmITN/5MBm5nrEvCSh3EBQE9iydqlpBVU+SiK4j/uv99awT35JCw4D4846jNjza+/Ad4h5hlPM+wWnJJm8Xt4HUVRbnEefNCm4W76C1R/CYKDo9qy54X33ocxOyF0FBgXH6Hdo/GckUVJC6jyURTF/zRuDIHb4OhJ2PAthIZGtZ0DHg2FbC1tROyswA0g3w6/iKr4Bt12UxQldVDgF2g20UbCPno0qj4H1gk1G1AX+BebCSzXdn9IqfgIVT6KoqQOCi6HPxtAt2etAvr336i2zNikK2FANaypdaHl/pBS8RGqfBRFSR1kOQ0PA8dehUea21QMLrlkMFjjg8LALuCaJgJKy6jyURQl9dAXm2Zh27uQqws80gquXLFnPAuA+liLuOnAzmF2O05Jk6jyURQl9WCAl4GZBmp0hd9HQaE/oNZ1+B2Y6WhvAZR/wyalU3+fNIkqH0VRUh+FgIEBcPwO6DkXjpeBR36x9U6y/wYbgHeBgWi8tzSGKh9FUVIvGTLAu+/aUq+edUh15Q5gM3b77Rns9pySJlDloyhK6ueJJ2w07D59YOBA94R0BYB1wAnsedAFv0ioeIkqH0VR0gYVK8IPP9iI2K1b2/TcToKBJUBxrDn2P36RUPECVT6KoqQdChaENWsgd274+Wc4cCCqLSMwHmuMUB34wy8SKglEw+soipI6yJ0bwsIS3v/226FaNZuYrnZtW2ew6biLYjOgLgBq+FpQxReo8lEUJXXgzOvjDd9+C48/Dm+9BT17gjG2vjPWGbUFdjXUymdSKj5Ct90URUm71K0LmzfD2LHQvTtcvx7V1hD4GngeGOUn+ZRYUeWjKEra5vbbYcsWOHHCKqOTJ6Pa7gc2AaOBfkCE5ymUlEeVj6IoaZ+QEFi0yJ4ZVa5sjRGclMQqoC1AO+CaPwRUoqPKR1GU9EFAAAwZAv/7HzRoAHPnRrXlwWY+vQ40As76RULFBVU+iqKkLx57DFatgn79YMCAKIfUIGAucA9QEzjiPxEVVT6KoqRH7rsPfvwR1q+Hli2jUjNkAD7GWsNVB371m4S3PKp8FEVJn+TPb02xCxSw/kB//WXrDTZ1w3CgHrDGfyLeyqjyURQl/ZI5M4wfb32Aqle3ysjJ49htuCeAz/0k3y2MKh9FUdI3xkCvXvDFF9CuHXzyCYgj/0IYduXzOjCMqLQMx4G3gX3YaNkPAUMc9YpPUOWjKMqtQe3a1h9o4kTo0gWuOWyu78KmZZgF9AL+B7QHygO3Y5PVfevo1x74AM0d5ANU+SiKcutQqpRVQGfOQJ068I8j/HUR4DusOfY4bITsVtjzIYBMjvergUPAiBSWOx2iykdRlFuLbNlg/nyoX986pG7fbusvYwOSVgPqA/9CQEAAGTJkiBobAIwEVqBbcElEA4sqinLrERAAgwbB3XdDo0b2HGgIYH6DiI1wtAsUf4g5pQtz+3PPRQUsdXLqIahbFn570h/SpwtU+SiKcuvy6KMQGgotWsDxL+DsYxDoWOmMgdJDv4T3gYejjbsB5NuRoqKmN3TbTVGUW5t777UZUsPDoXXzKIfUnsAMrBFCe+CUy5hMKS5lukOVj6IoSr58EJwNipSEqlXhzz9tfV1gJ1AQuBuYibV0u+EnOdMRqnwURVEAcm+Dhp/CCy9AjRo2PhxAMNa8+kvH68PAJCDXdn9Jmi5Q5aMoigJQcLnN+9PtWZg3Dzp2hI8+inJIrQT8iHU4fQEQAzf9Jm2aR5WPoigKQJbTdlXzIvBgLesPNHOmTVC3e7ftkwE4BrwEnL3fmmVrcNJEocpHURTFSV+gBNbPZ3tJ2Pi9jYpdqx40nQK1w21yuveBe/tAN+y50ADgqr+ETpuo8lEURXFigJexhgW/AY0ywvzn4Y5D8E8u+L0i5JkCEgFGoAt25bMHuBfY4D/R0xrq56MoihKdQsCbjgJAZqAF/FgEnnvORsq+cSOq73xgEfAk0BQbpDRHyoqc1tCVj6IoSkJ54AF7FtStG+zcaV///de2tQR2OfrdhVVGSqyo8lEURfGGgAB4+mkbFy4oCMqVg9GjrZNqTmxg0s+B14BH0RhwsaDbboqiKAC5c0NYWML7588PI0fa9AwvvAATJsCnn8KDD1pz7F+AodizoHexeYFMXBPeWqjyURRFAVi4MHHj7r4b1qyBuXPhySfhoYdg+HAoXNgqn8ewhgkzgQnAHT6TOE2j226KoihJxRho2xZ+/x1KlIB77oH//Q+uX4d7gC1AC6A68B4angdVPoqiKL4jWzZ4911rlLB+vV0Vff21dU7tDWwD1gMPOH6+hVHloyiK4mtKl4bly2HECOjVyzqqHjhgHVRXYH2JmmCdWi9FG3sceBt7blTL8TqEdGe44HflY4wpZoyZb4w5Z4w5b4xZaIwpnsCxgcaY4caY48aYK8aYLcaYh5JbZkVRlATRtCns2gWVKtkycCBcuWxTNOwCTmCjZa/CRsv+ANtWHvgWu0r6Fmu63d7RLn74HMmAX5WPMSYrsAYoC3QEOgClgbXGmOAETPEZ0BV4C+vadRz42hhTIVkEVhRF8ZbAQOjfH37+2Z4JlSsHixZBXrFGCKOxYXoeAP4AVgOtiMoZlMnxfjVwCBiR8h8hOfD3yqcrcBvQQkQWi8gSoBk2ulL3uAYaY+7F+hP3EZGJIvIt1q7kMHbR6nMmTJiQHNPeEui9Sxp6/xJPqrl3xYtbi7jJk2HAAGjYEPbsscFM7+gJu07C1NNQ/m2oFWbNvl1LnTD4tTa8sx0aP5MiIv/+++/JN7mI+K1gF5SbPNSvB9bHM/ZN4DqQNVr9YOAakCX6mIoVK0pSSOr4Wxm9d0lD71/iSZX37vp1kY8+EsmTR+Tll0WKTxRZICJbRaS8iDQWkUOxjJ0vIiUnJbuIBw8elLJlyyZpDmCbxPIM9/fK5y6iAlK4shsol4CxB0TksoexmYHQpIunKIqSDGTKBL1721QN//4LR8rA+VlQWWA7NlXD/cBHQPQnXDPgTMVkF7FEiRIEByfk9CNx+NvJNDdwxkP9aSBXEsY6293Yvn37RWOMq8I9BfybADmd5DXGeNNfiULvXdLQ+5d4Uv29u48tFX7u/McOOrd3b3jJUWL0L1PhZ/PHjhQQLTH3Li+Qz/Fzmdg6+Vv5pCgiEuJvGRRFURT/GxycwfMKJ7ZVTULHQtQKSFEURUll+Fv57Mae3USnHDaVU3xjSznMtaOPvQ7sS7p4iqIoSnLgb+WzFKhqjLnNWWGMKQnUcLTFxTKsBXwbl7EZgbbAKhG55gsBk+IEe6tjjAkzxoiHctbfsqUmjDFFjTGjHE7Slx33qKSHfupUHQ0v7p2nv0O5lX0CjTGtjTELjDGHHH9Pfxhj3jPGhETrl8sYM8kY868x5pIx5htjzN1Jvb6/lc9E4CCwxBjT3BjTDFgCHAHGOzsZY0oYY8KNMW8560TkZ2AO8LExposxpi7wBVAKGOgL4XzgBKtYXsDa7zhLPf+Kk+oIxfqonQG+i6OfOlXHJKH3DmAq7n+H1YC9ySlcKudl4CbwBtAIGAv0AFY7DbOMMQb7Rb8R8Dw2Q1Em7DOwaJKuHpsNdkoVoDiwADgPXAAWAyWj9SmJDSoxKFp9EPAh8A9wFfgeCPOhbC86fjmhLnWlgHDgJX/fu9RegDDH762ev2VJzQUIcPm5i+OelYzW515HfWeXuoxYn/il/v4MqfneOdoEGOpveVNTAfJ5qHvKca/qON43d7yv7dInB/ZM/ZOkXN/v1m4ichirTePqcxAPaZhE5AqxGiP6hGbAVhGJPD8SkQPGmE3YX8qHyXRd5RZCRCIS0K0ZNhD/HJdx4caYL4DXjDFZxEdbzWmJBN47xQMicspD9Y+O1yKO12bAMRFZ6zLunDFmGfYZ+EJir+/vbbfUTlKcYJUoZhljbhpj/jPGfK5nZolCnaqTTg9jzDXH2dAaY8yD/hYoFVLL8eqMqxPXM7C4MSZbYi/k95VPKicpTrAKnMOGQVyP3Va9D7u/vMUYc5+InPSncGkMr52qFTdmAl8Cx7CxI/sBa4wx9UVknT8FSy0YY4pg42J+IyLObEO5sefy0XH+3eUCLibmeqp8lGRDrFHIzy5V640xG4AfsMv1AX4RTLnlEJEOLm+/M8YswX6jHwrU9I9UqQfHCmYJ9jy7c0pcU7fd4iYpTrCKB0TkJ6yF0QP+liWNoU7VPkRELgDL0b9DjDFBWIu224CGIvK3S3N8f3eJfg6q8ombpDjBKnGTTlJipRjqVJ083NJ/h8aYTMB8oBLQWER2RusS1zPwsIgkassNVPnER1KcYBUPGGMqYYMN/uBvWdIYKeJUfatgjMmO9ZW6Zf8OHb48s4A62JxqWz10WwoUMcbUchmXHXiEJD4D9cwnbiYCz2GdYAdgvyUNIZoTrOIZY8ws4ADwE3AWa3DwOnAU+MR/kqU+jDGtHT86Y+U/bIw5BZwSkfUi8rMxxulUnQl7X3tg/c7apbzEqYf47p0x5mXsF561RBkcvAwU5Na+d6OxX2beAS4ZY6q6tP3t2H5bCmwBZhpj+mG32V7Hur78L0lX97ejU2ovJMAJVkus9+514Fes1dsNrNKeABTyt2yprWC/2Hgq61z6JKtTdVot8d077Lf0Tdj0KTeA/xwP1cr+lt3P9+1gHPdukEu/3MBk7LniZWwS0HuTen3jmFxRFEVRUgw981EURVFSHFU+iqIoSoqjykdRFEVJcVT5KIqiKCmOKh9FURQlxVHloyiKoqQ4qnwURVGUFEeVj6IoipLiqPJR0iTGmInGGDHGfOTDOQcZY5LF6zohcxtjFhtjzhhjssTSHmKMuWSMmZqcMjpixiUrxphCxpgIY8wtn87gVkWVj5LmcISAf8zx9kkfPiwnAdV8NFdimAbkxAa89ERrIKujX1qnOXAK2OxvQRT/oMpHSYu0ALIDXwH5gUbxDYhtNeHaJiJ/i+fIvinFcmzcsadiaX8KOAysSymBkpEWwDIRifC3IIp/UOWjpEU6YqPrdgKuON5H4rJ9VN4Y87Ux5iIwN6Ftjp/bOPrdE/3ixpivjDG/uLwPNcbMMMYcMMZcMcbsN8aMNcZ4lWpdRK4Ds7FRmfNEu2ZxoBYwQ1IwIKMxppEx5qIx5lNjTIDL/SvruH+XjDGHjTGdHf07GGP2OMasNcbc7mHO7EBtbJBejDF3GGMWGWNOGmOuOuablxLbf4r/UOWjpCmMMYWBesAcETmFfYA9EsuDfgmwHmgGRD8biqsNbP6cc0D7aNcvADQAprtUF8ZG7O4NNATeBupiV2beMg2bt+fxaPXtsWHsp8cYkUwYY57CRn9+X0Sei7ZKmYddqbUAtgOTjTHvYtM8vIZNxVwG+NzD1I2xCfC+cbxfDhRxjG3oGH8NfT6lb/wd1luLFm8K8Ao25Hs1x/uGjvfPuvQZ5Kh70cP4eNtc3k8E/gYCXOp6Y/Pcx5oWApsnq6bjOvd5mjuez7gb+D5a3e/AlmS+t857k9Fxn28AXWLp85RLXS7HPfkPyO5S/4Kjb4loc3wBzHf8nNfRp5m//7a0pGzRbxZKWqMj8KeIbHG8/wabIKyjh76L4pgnrjYn07HfyOu41HUAvhWR484KY0xmY8wbju2mK9iH9neO5jIJuE50pgGVjTF3OOavDJQl5QwNPgIGA61FZFIsfVY4fxCRM8BJYKuInHfps8fxWsxZYYzJDDyMY8sNq7D2A+8bY7oaY0r75BMoqR5VPkqawZGCuxyw0BiT0xiTEwgBFmLTnd8RbchxYieuNicbsQm3OjiufydwPzG3vt7DrghmAk2AykArR1tgAq4TnZlABFGGB09ht6HmJGKuxPAEsIuobTFPnIn2/nosdeB+D+pgLfa+BMdSEOoD27D3ca/jzKxH4kRX0gqqfJS0hHN18yr2Qecszznqo1uJxXUwH++hvePBOBNoZYzJilVCF4m5anocmC4iQ0VkjYj8iE0bnihE5BiwGmjvWCm0xVqGuT3cjTH5jTFfGGPGG2OyGmN6OowhGiT22g7qYjP4rjDGZEviXNFpAawXkbPOChHZLyJPAfmwqdbXAGOMMQ/7+NpKKkKVj5ImcDyEn8Cmjq7toewAOhhjjI8vPQPIhl3JtAMWisjlaH2yYrfaXOmcxOtOA0pgVwN58bzl9irwFrAWu5KojL1H7Rz3K7HsBsKA0vhQATl+N82I2nJzQyw7gJccVeV9cV0ldaKmjEpaoQmQB+grIuuiNxpjxgNjsQ9NnyEie40x3wPvY89/PFmbrQQ6GmN2Avuwiqp6Ei+9GDgP9MGep6z00CdARPZit6p6Y63SzhljfgQKYn2CMMaUBA4Ag0VkUEIuLiK/G2PCsIrta2NMIxG5kITPA1AFKISL8nGYso/EbinuAzJgTejDsSsgJZ2iKx8lrdARuIA18fXEbDz4/PiIGVjFcxT7MI7O81iT5HewD9EQ7Aok0YjIFaz/kQE+F5FwT92MMbc7lMQ1YJAj+kMl4B+XfsGO13/wAhH5A+tbVAJY5fDPSQotgO0i8rdL3T9YJfkS9h7OxpquNxWR7Um8npKKMXZbW1GUtIYxpiDwMXAT6IUNOfQY8ImILHXp1w2rGEt42DJMMYwxe4CZIjLUXzIoqQdVPoqSzjHGzAJ2i8i7/pZFUZyo8lEURVFSHD3zURRFUVIcVT6KoihKiqPKR1EURUlxVPkoiqIoKY4qH0VRFCXFUeWjKIqipDiqfBRFUZQU5/8SD5qdptiB/wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fP_prop_pp = compute_fP_prop_pp(vinf_array)\n", "fP_prop_ab = compute_fP_prop_ab(vinf_array)\n", "\n", "Q_array = np.array([41.89E3, 48.6E3, 64.47E3, 83.6E3, 103.6E3,\n", "\t 126.0E3, 152.5E3,183.5E3, 222.8E3, 275.7E3, \n", "\t 352.0E3])\n", "\n", "tps_mf_array = 0.091*np.power(Q_array, 0.51575)/100.0\n", "ess_mf = 0.23\n", "\n", "fP_aero = 1.0 - ess_mf - tps_mf_array*0.60\n", "\n", "\n", "fig = plt.figure()\n", "fig.set_size_inches([6.25,6.25])\n", "rcParams['font.family'] = 'sans-serif'\n", "rcParams['font.sans-serif'] = ['DejaVu Sans']\n", "params = {'mathtext.default': 'regular' } \n", "plt.rcParams.update(params)\t\n", "\n", "plt.plot(vinf_array/1.0E3, fP_prop_pp, linestyle='-', linewidth=1.00, marker='s',ms=12, markerfacecolor=\"None\", markeredgecolor='red', markeredgewidth=0.75, color='red', label='Propulsive')\n", "plt.plot(vinf_array/1.0E3, fP_prop_ab, linestyle='-', linewidth=1.00, marker='o',ms=12, markerfacecolor=\"None\", markeredgecolor='magenta', markeredgewidth=0.75, color='magenta', label='Aerobraking')\n", "plt.plot(vinf_array/1.0E3, fP_aero, linestyle='-', linewidth=1.00, marker='^',ms=12, markerfacecolor=\"None\", markeredgecolor='blue', markeredgewidth=0.75, color='blue', label='Aerocapture')\n", "\n", "plt.xlabel(\"Arrival \"+r'$V_\\infty$'+r', km/s' ,fontsize=16)\n", "plt.ylabel('Payload mass fraction, '+r'$f_P$' ,fontsize=16)\n", "ax = plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "plt.legend(loc='upper right', fontsize=16, frameon=False)\n", "ax.set_ylim(0,1.03)\n", "plt.xticks(np.linspace(0,20, 5, endpoint=True), fontsize=16)\n", "plt.yticks(fontsize=16)\n", "ax.set_ylim(0,1.03)\n", "\n", "ax.annotate('Aerocapture \\ninfeasible',\n", " xy=(10, 0.3) ,\n", " xytext=(4, 0.3),\n", " arrowprops=dict(arrowstyle=\"<-\"), va=\"center\", ha=\"center\", fontsize=16)\n", "\n", "\n", "ax.annotate('',\n", " xy=(vinf_array[6]/1.0E3, fP_prop_ab[6]) ,\n", " xytext=(vinf_array[6]/1.0E3, fP_aero[6]),\n", " arrowprops=dict(arrowstyle=\"<-\"), va=\"center\", ha=\"center\", fontsize=16)\n", "\n", "\n", "\n", "plt.text(vinf_array[6]/1.0E3, 0.5*(fP_prop_ab[6]+fP_aero[6]), \n", "\t str(int((fP_aero[6] - fP_prop_ab[6])*100/fP_prop_ab[6]))+'%' , \n", "\t fontsize=16, ha='center', va='center',\n", "\t bbox=dict(facecolor='white', edgecolor='white', pad=0.0)) \n", "\n", "ax.annotate('',\n", " xy=(vinf_array[8]/1.0E3, fP_prop_ab[8]) ,\n", " xytext=(vinf_array[8]/1.0E3, fP_aero[8]),\n", " arrowprops=dict(arrowstyle=\"<-\"), va=\"center\", ha=\"center\", fontsize=16)\n", "\n", "plt.text(vinf_array[8]/1.0E3, 0.5*(fP_prop_ab[8]+fP_aero[8]), \n", "\t str(int((fP_aero[8] - fP_prop_ab[8])*100/fP_prop_ab[8]))+'%' , \n", "\t fontsize=16, ha='center', va='center',\n", "\t bbox=dict(facecolor='white', edgecolor='white', pad=0.0)) \n", "\n", "ax.annotate('',\n", " xy=(vinf_array[10]/1.0E3, 0.0) ,\n", " xytext=(vinf_array[10]/1.0E3, fP_aero[10]),\n", " arrowprops=dict(arrowstyle=\"<-\"), va=\"center\", ha=\"center\", fontsize=16)\n", "\n", "plt.text(vinf_array[10]/1.0E3, 0.5*(0.0+fP_aero[10]), \n", "\t r'$\\infty$' , \n", "\t fontsize=16, ha='center', va='center',\n", "\t bbox=dict(facecolor='white', edgecolor='white', pad=0.0)) \n", "\n", "\n", "plt.savefig('../data/jsr-paper/uranus/uranus-mb.png', dpi= 300,bbox_inches='tight')\n", "plt.savefig('../data/jsr-paper/uranus/uranus-mb.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../data/jsr-paper/uranus/uranus-mb.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "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.7.3" } }, "nbformat": 4, "nbformat_minor": 5 }