Skip to content

Instantly share code, notes, and snippets.

@TTY6335
Last active December 4, 2021 13:01
Show Gist options
  • Save TTY6335/0410d588bd72d9b997e461a935835f1c to your computer and use it in GitHub Desktop.
Save TTY6335/0410d588bd72d9b997e461a935835f1c to your computer and use it in GitHub Desktop.
Didymosを2Dで表示する
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 46,
"id": "ee00adcb",
"metadata": {},
"outputs": [],
"source": [
"from astropy.utils.data import conf\n",
"conf.remote_timeout = 10000"
]
},
{
"cell_type": "code",
"execution_count": 47,
"id": "894cf19c",
"metadata": {},
"outputs": [],
"source": [
"#astropyより単位と時間のモジュールを呼び出す\n",
"from astropy import units as u\n",
"from astropy.time import Time\n",
"\n",
"#poliastroのclassで定義されている惑星を呼び出す\n",
"from poliastro.bodies import Sun, Earth, Mars, Jupiter\n",
"from poliastro.ephem import Ephem\n",
"from poliastro.frames import Planes\n",
"from poliastro.twobody import Orbit\n",
"from poliastro.plotting import StaticOrbitPlotter\n",
"from poliastro.plotting.misc import plot_solar_system"
]
},
{
"cell_type": "code",
"execution_count": 48,
"id": "1e2d8804",
"metadata": {},
"outputs": [],
"source": [
"#計算する時刻\n",
"EPOCH = Time(\"2021-11-24 00:00:00\", scale=\"tdb\")"
]
},
{
"cell_type": "code",
"execution_count": 49,
"id": "b5a04e76",
"metadata": {},
"outputs": [],
"source": [
"#NAS/JPLよりdidymosの軌道データを呼び出して位置を計算させる\n",
"didymos = Ephem.from_horizons(\"65803\", EPOCH, plane=Planes.EARTH_ECLIPTIC)\n",
"didymos\n",
"\n",
"#NAS/JPLよりdidymosの軌道データを呼び出して軌道を計算させる\n",
"didymos_orbit= Orbit.from_sbdb(\"65803\")"
]
},
{
"cell_type": "code",
"execution_count": 50,
"id": "dc565bcd",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAGBCAYAAABvrRvrAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACW4UlEQVR4nOzdd1iV9f/H8ec5h71BZMgQGe6tuBX3yJmZudPKkebISsuG9W3YttS01Fyp5chy771FcSAgypK99+Zwzu8Pkl+kKOtwH+DzuK77Ss45932/Dia8z2fKADWCIAiCIAi1kFzqAIIgCIIgCJoiCh1BEARBEGotUegIgiAIglBriUJHEARBEIRaSxQ6giAIgiDUWqLQEQRBEASh1qrxhc6vv/5KXFwcvr6+z3ytk5MTp06dwsfHh9u3bzNkyJBqSCgIgiAIglRk1PB1dHr27ElmZiZbtmyhVatWT33tL7/8ws2bN/n5559p1qwZhw4dolGjRtWUVBAEQRCE6lbjW3TOnz9PcnJyicdcXV05fPgw169f59y5czRp0gQAtVqNmZkZAObm5kRHR1d7XkEQBEEQqpe6ph8NGzZU+/r6Fn994sQJtbu7uxpQd+rUSX3y5Ek1oLazs1PfuXNHHRERoU5OTla3b99e8uziEIc4xCEOcYhDc4cOtYyxsTHdunVj165dxY/p6+sDMH78eDZt2sT3339Ply5d+O2332jZsiVqtVqquIIgCIIgaFCtK3Tkcjmpqam0a9fusedeffVVBg8eDMCVK1cwMDDA2tqahISE6o4pCIIgCEI1qPFjdP4rIyOD0NBQxowZU/xY69atAQgPD6dfv34ANG3aFAMDA1HkCIIgCEItJ3n/WWWO7du3q6Ojo9X5+fnqiIgI9SuvvKJ2cXFRHz58WH3r1i21n5+f+sMPP1QD6mbNmqkvXLigvnXrlvrmzZvqAQMGSJ5fHOIQhzjEIQ5xaO6o8dPLBUEQBEEQSlPruq4EQRAEQRAeEYWOIAiCIAi1Vo2edRUfH8/Dhw+ljiEIgqA1GjZsiI2NTZVe09LSko8//hh3d3fkcvH5WNAuKpWKoKAgPv74Y1JSUh57vkaP0fH29sbT01PqGIIgCFpDEz8Xf/zxRzp16oSOTo3+bCzUYkqlkmvXrjF//vzHnhOluSAIgvBU7u7uosgRtJqOjg7u7u5PfE4UOoIgCMJTie4qoSYo7f9T8X+vIAiCUCN5enqyfPny4q9/++031q5dK2Giyrtx4wZvvvlmtd937dq1/Pbbb9V+3+ogCh1BEAShRtLT0+PMmTOkpqZKHUXQYqLTVRAEQaiRFAoFo0aNYvv27cyePbvEc+fOnWPDhg0UFBRgbm7Op59+Sr169Vi7di3R0dEkJiYSHh7OggULuHv3LpcuXaJ+/fosX74cHR0dAgICWL58OTk5OVhYWLB06VKsra35448/2LNnDwqFgkaNGvHFF1+UuG90dDRLly4lJycHgHfeeYc2bdpw48YN1q5di4WFBcHBwTRt2pRPP/0UmUzGpUuX+P7777GwsKBp06ZPfK+FhYWsWrWKGzduUFBQwIsvvsjo0aO5ceMGv/zyC+bm5jx8+JB27dqxePFi5HI5R48eZePGjajVanr06MHcuXMBuHTpEqtXr0alUmFubs6aNWsACA0NZebMmcTGxjJ+/HjGjRtHTk4O7733HvHx8RQWFvLqq68ycODAqv6r1ChR6AiCIAg11osvvsiECROYMmVKicfbtm3Lxo0bkclk/P3332zZsqW4SygyMpKff/6ZkJAQXnnlFb766ivmzZvHO++8w4ULF+jRowfffPMN3333HZaWlhw7dozVq1fz0UcfsXnzZvbu3Yuenh4ZGRmP5bGysmLVqlXo6+sTHh7OBx98wJYtWwAIDAxkx44d1K9fn9dee43bt2/TrFkzPv/8c9asWYOTkxNLlix54vvcu3cvJiYmbNmyhfz8fF577TU6d+4MgJ+fHzt27MDe3p558+Zx+vRpWrduzcqVK/ntt98wNTVl7ty5nDlzhjZt2vD555+zdu1aHBwcSEtLK75HWFgYP//8M9nZ2YwZM4YxY8Zw+fJlrK2t+eGHHwDIzMys9N9ZdROFjiAIglBjmZiY8Nxzz/HHH3+gr69f/Hh8fDxLliwhMTGRgoICGjRoUPxct27dimfpqFQqunXrBoCbmxvR0dGEhYUREhLCnDlzgKJ1WqytrYGiGWgffvghXl5e9O7d+7E8SqWSr7/+mvv37yOXywkPDy9+rkWLFtja2gLQuHFjoqOjMTQ0xMHBAWdnZwAGDx7M33///dh1r169SlBQECdPngQgKyuLiIgIdHV1adGiBY6OjgAMGjSIW7duoaOjQ4cOHbC0tCy+7s2bN5HL5bRr1w4HBwcAzM3Ni+/Ro0cP9PT00NPTw9LSkqSkJNzc3Pjxxx9ZuXIlPXr0oF27duX429EOotARBEEQarTx48czefJkhg8fXvzYN998w4QJE/Dy8iruNnpEV1cXKJqlo6Ojg0wmK/66sLAQAFdXVzZs2PDYvX744Qdu3rzJuXPn+PXXX9mxY0eJqffbt2/HysqK7du3o1Kp6NGjR/Fzenp6xX/+973KQq1W8/bbb9O1a9cSj9+4ceOx18pkMtTqJy+Rp1ari9/vfz36vvw7X8OGDdmyZQsXL17kp59+onPnzkyfPr3MubWBGIwsCIIg1Gjm5ub079+fvXv3Fj+WmZlZvEL0gQMHynW9hg0bkpKSwp07d4CiVprg4GBUKhVxcXF07NiRefPmkZmZWTwW59/3tba2Ri6Xc+jQoWcWMy4uLkRHRxMZGQnAsWPHnvi6Ll268Oeff6JUKgF4+PBh8b39/f2JiopCpVJx/Phx2rZtS8uWLfHx8SE1NZXCwkKOHj1K+/btad26NT4+PkRFRQGU6Lp6koSEBAwMDHjuueeYNGkSgYGBz/juaR/RoiMIgiDUeBMnTmTnzp3FX8+YMYN3330XGxsbWrZsSXR0dJmvpaury5dffsl3331HZmYmSqWS8ePH07BhQz766CMyMzNRq9WMHz8eU1PTEueOGTOGxYsXc/LkSTp06IChoeFT76Wvr8+SJUtYsGABFhYWtGnThpCQkMdeN2rUKGJiYpg0aRJqtRpLS0u+/fZbAFq1asWqVasIDg6mXbt29O7dG7lczpw5c5g1axZqtZru3bvj5eUFwJIlS1i0aFHxdX766adS8wUFBbFixQpkMhk6Ojq8++67Zf4+aguxBYQgCEItoomfi4cPHy4eoyJolxs3brB169YS6wnVVYmJiQwZMuSxx0WLjiA8g7GBLob6OujpKtDX1UFfV4GejgJ9XUVxX7gaUKnVoC76b15BITl5BeTmK8nJU5KTryS/oOz98YIgCELV0KpCZ8GCBbz22muo1Wp8fX2ZNm0aeXl5UscSaiFdhRxrCyNsLIyxsfz//1qaGGBmrI+5sT5mxvqYGemjo6iaoWyFKhXpWfmkZuaSkpFDSmbeP3/OJT41i5ikTGKTs0hKz6aUcYSCIAgldOjQgQ4dOkgdQ6tpTaHToEED5s2bR/PmzcnNzWXHjh2MGzeOzZs3Sx1NqMGszQ1xtjHH2cYMJ1uz4j/bWBo/9trUzFyS03NIz84nPC6dtKw80rPySMvOIyevgPyCQvL+OYr+rESlLur/lcv/mbUhkyGTydDTlWOor4uhng6G+jrFfzY30cfCxABLEwOaOFlhYWKAiaFeiRz5ykLikrOITc7kYVwaYbFphMWmEhabRlZuQXV82wRBEGoNrSl0oGj3UUNDQwoKCjAyMirX4DFBaFDPhMZOVjR2rEdjJys8HKwwM/7/dTWycosKmFtBcUQlZRCfkkV8SjbxqVnEp2ZL1rWkp6vA1sIYu3rG2FmZYG9lgn29omNoFw8M9f//n2lCahahMWkERiRxLzyJgPBEUjJyJcktCIJQE2hNoRMdHc23335LeHg4OTk5HDt2jOPHjz/2uunTpzNjxgwAMTiuDtPVkdPUuR6tXW1o7WpLs4b1MDUqKmoKlIWExqRy7k44wVEpPIxLIzw+naT0nGdcVRr5BYVEJKQTkZD+2HMyGdhaGuNiZ4GLnTmN7CxwbWDBhH4tUPzTpRaXksW98EQCHiZxOziO+5HJqFSi70sQBAG0qNCxsLBg5MiRNGrUiNTUVHbt2sXEiRPZtm1bidetW7eOdevWAUWzC4S6QSGX0bJRfdo3tqeNqw3NGlqjp6sAIDQmlbO3wwmMSOJ+RDKhMakUFKokTlw11GqITc4iNjmLK/5RxY/r6yrwcLSiqXM9mjpb09S5Hl5tGgJFLVe+IQncfBDLraA4gqJSigZKC4Ig1EFaU+j079+f0NBQEhMTAdizZw/dunV7rNAR6g5rcyM6N2tAp6YNaN/YDhNDPQoLVTyISuHvC4HcCYnHNzSB9Ky6N2A9r6CQu6EJ3A1NKH7M0sSAtu62tPWwpa27HV2aFy3xnpaVx7WAaC77R+J9L4bMnHypYguCIFQ7rSl0wsPD6dKlC4aGhuTk5NCvXz+uX78udSyhmrnYmePVpiE9Wzvh1qBoj5b4lCzO3HrItYBofB7EigG5pUjJzOX0rYecvvUQgHpmhrR1t8WzqT2dmzkwoGMjCgtV+IYmcNkvkgt3I4hOrHkb9AnaT6UG3zgZYWkyXMzVtLJVI3/yrgNlEhsby8cff0xSUhIymYznn3+e8ePHA0Ur+y5ZsoSYmBjs7e1ZtmwZZmZmXL16lVWrVlFQUICuri7z5s0rXl9o9erVHDx4kIyMDM6dO1fqfUt7nY+PD99//z1BQUF8/vnn9OvX74nnb9u2jb1796JQKLCwsOCjjz7C3t6++PnMzEzGjh1L7969WbRo0WPnl/beADZu3Mi+ffuQy+VP3BqiKs4HWLx4MXPnzsXR0ZERI0ZgZGSEXF7Ubd6+fXvefvvtUr9//3Xjxg10dHRo06YNAB9//DE9e/Z87PuXkpLCRx99xMqVK8t87afRmkLn2rVr7N69Gx8fH5RKJTdv3iyxN4lQe7k1sMSrjTNebZxxtjVHpVLjGxrPmr03uHYvmrDYpy9RLjxZUnoOJ33COOkThlwmo6lzPbq2cKBLcwdeH9mB10d24H5EEqduPuTMrYfEpWRJHVmoBVRqWHFVh5ux/78sQzs7FfM6Kytc7Ojo6LBgwQKaNm1KVlYWU6ZMoXPnzri6urJ582Y8PT2ZOnUqmzZtYvPmzcydOxcLCwu+//576tevT1BQEPPmzePQoUMA9OzZk7FjxzJ69Oin3re019nZ2bF06VK2bt361PObNGnCli1bMDAwYPfu3axYsYJly5YVP//zzz/Tvn37Us8v7b2FhIRw/PhxduzYQUJCAnPmzOHPP/9EoVBU6fmPtr14tGHoo8wWFhZPfd9PolQquXHjBoaGhsWFTmksLS2xtrbm9u3bz3xtWWhNoQNF1d3HH38sdQyhGtS3MGJgx0YM9HTD2caMQpWK20Fx/HnuHud9I8RMoiqmUqvxf5iI/8NEfj10GxsLI7zaNKRPu4bMGtGeWSPaczc0gTO3igqj1My61x0oVA3fOFmJIgfgZqwc3zgZbewqNlbM2tq6ePKJsbExLi4uJCQk4OrqytmzZ/nll18AGDZsGDNnzmTu3Lk0adKk+Hw3Nzfy8/PJz89HT0+PVq1alem+pb3u0U7opW2O+UjHjh1LXOvw4cPFXwcEBJCcnEzXrl0JCAh44vmlvbezZ88yYMAA9PT0cHBwwMnJCT8/P1q3bl2l5x85coRevXo99T0C/PXXX/z1118olUocHR353//+h4GBAR9//DHm5uYEBgZiZmbG7du3USgUHD58mHfeeQeAmzdvsm3bNpKSkpg3b15x607v3r05fPhw7St0hNpNX1dBz9bODPJ0pb2HHXK5jNtBcew648/5OxGk1cGxNlKJT81m19kAdp0NwL6eCb3bNqRvu4a88bwns0Z04NLdSA5dDcL7XowYyCyUS1jak3/5P0yreKHzb9HR0QQGBtKiRQsAkpOTi4sga2trUlJSHjvn1KlTNG7cuMTu4dVt7969dOvWDQCVSsUPP/zAJ5988tRJNaW9t4SEBFq2bFn8OhsbGxISisbrffbZZ4wePZrmzZtX6Px/u337NoMGDSrx2KxZs4q7roYNG8aECRPo06cPzz//PABr1qxh7969vPTSS0DRsJSffvoJhULB2rVrMTQ0ZPLkycXfk8TERNavX09YWBhvvfVWcaHTrFkz1qxZ8+xvbBmIQkfQOBc7c0Z1b0z/jo0wNtAjOimDLcd8OXY9hJgkMUZEajFJmfx+0o/fT/rhYmfOkE5uDPR0pVcbZxJSszhyLYRDV4OITRZdW8KzuZg/uZhpWMrj5ZGdnc3ixYtZuHAhJiYmZTonODiYlStXsmrVqkrfv6IOHTpEQEBAcevK7t276d69O3Z2dhW6nvoJHz4etS598MEHlTr/35KSkh7rpnpS11VwcDA///wzGRkZ5OTk0KVLl+Ln+vXr91iX2L95eXkhl8txdXUlOTm5+HErK6viyUmVJQodQSPkchk9WjoyqkcT2nnYkV9QyOlbYRy+GsydkHixxYGWCotNY80+H9YdvEW3Fg4M6ezOhP4tmNi/JZf8Ivnz3D1uBcVJHVPQYq1s1bSzUz02RqeVbeX+0SuVShYvXszgwYPp27dv8eOPfiFaW1uTmJiIpaVl8XNxcXEsWrSITz75pMQ4kycpLCwsbmno1asXs2bNKnfG1atXc+HCBQC2b98OwNWrV9m4cSO//PJLcYvSnTt3uHXrFrt37yY7OxulUomhoSFz584tcb3S3puNjQ1xcf//7zA+Pv6J68pV9nx9fX3y8589S/N///sf33zzDY0bN2b//v3cuHGj+Lln7d7+71a2fxdgeXl56OvrP+mUchOFjlClTAz1GNndgxHdGmNjaUxscia/7Pfh0NXgOjkNvKZSFqo4dyeCc3cisDY3Yng3D0Z086BHKyeCo1LYc/4eJ3zCxEalwmPkMpjXWYlvnIyHaTIaVsGsK7VazaeffoqLiwsTJ04s8VyvXr04cOAAU6dO5cCBA3h5eQGQkZHBm2++yZw5c8o0zkOhUBQXJxU1e/ZsZs+eXfx1YGAgy5YtY8WKFVhZWRU//tlnnxX/ef/+/QQEBDxW5EDp761Xr158+OGHTJw4kYSEBMLDw4u78qryfBcXFyIiIorHJJUmKysLa2trlEolR44coX79+k98nZGREVlZZWsZDg8Px9XVtUyvfRZR6AhVwtrckDFezRje1QMjA128A6P58U9vrvhHiTEeNVxiWjYbD99m24m79Gvnwgu9mvLOuK5MH9aOPecC2XP+npjyL5Qgl0EbO3WVjMmBorEihw4dwt3dnQkTJgAwZ84cunfvzssvv8x7773Hvn37sLW15csvvwRg586dREREsH79etavXw/AqlWrsLKyYsWKFRw9epTc3FyGDh3KyJEji1fc/7fSXufn58eiRYtIT0/nwoUL/PLLL+zcufOx83/88UdycnJ49913gaLZWt9//32Z33dp783NzY3+/fszduxYFAoFixYtKu4e+vcYnYqc/289evTAx8eHzp07Fz/27zE6Hh4efPLJJ8yaNYtp06ZhZ2eHu7t7qcVMz549effddzl79mzxYOTSXL9+nR49epT5e/U0MqDG/hby9vYuXhdBkIZTfTNe6tucgR0bIZfLOH3zIb+f9CMkJlXqaIIGtXW35cXezejWwpHMnHz2XrjPrrMBYkC5FtDEz8XDhw+LLXfqoNzcXF5//XXWr1//1HE2mjBjxgy+/fbb4nV/yiIxMZEhQ4Y89rho0REqpIG1CdMGt6FvOxcKlIUcvBLEzjMBYnBxHXErKI5bQXG4NbBk0oCWjO/XgtG9mnLg8gN2nPbX2n3FBEEoOwMDA2bMmEFCQkKFB05XREpKChMmTChXkfM0otARyqW+hRFTBrZiSCc3CgpV7Djtz64zAaRkinVv6qLg6BQ+2XweJxszJvZvweieTRjezYM95+/x+0l/sd2EINRwpa2YrEmWlpb07t27yq4nCh2hTCxM9Jk0oBXDu3kAsPfifbaeuCsW9hMAiIhP58vtl9l81Jepg1szrk8Lhnf1YPtJP/acDxSDlgVBkIwodISn0lXIecGrKZMGtMRAV4cj3iFsOXqH+NRsqaMJWigmKZNl2y6x87Q/rw1tx8zh7RndsymbjtzhyLVgMTBdEIRqJwodoVQ9Wzkxc0R7HKxNuXg3kl/2+RCRkC51LKEGCI5O5b11p2ntasOM4e14Z1wXRvVozIo93iV2XBcEQdA0UegIj3FrYMGcUR1p52FHaEwqb685wY37sVLHEmqgOyHxvPHjUXq3deb1ER1YOW8Qx6+HsvaAD4lpYsCyIAiaJ3/2S4S6wkBPh9dHtOeXt57D1d6C5buu8tq3B0WRI1TamVvhvPzlPrYc88WrjTNb3hvBuL7N0VGIH0G1kloF4QHgc6Lov2pVpS4XGxvLrFmzePHFFxk7diy///578XNpaWnMmTOH0aNHM2fOHNLTi1qdr169yuTJkxk3bhyTJ08usafU6tWrGTp06DM3rCztdT4+PkyaNIkuXbpw8uTJUs/ftm0bY8eOZfz48bz++uvExMSUeD4zM5PnnnuOr7/++onnl/beADZu3Mjzzz/PCy+8wOXLlzVyPsDixYuJjIwEYMSIEUyfPr3E8xMmTCje16qylixZQnh4eJVc69/ETxkBgK4tHNj07nDG9mnOwStBTPpiH/suPUClEmMqhKqRm1/IxsO3mfrVfm7cj2Xm8Pb8snAIzZzrSR1NqEpqFRzbjOzIr8iuH0F25Fc4trlSxY6Ojg4LFixg165dbNy4kd27dxMSEgLA5s2b8fT0ZM+ePXh6erJ582YALCws+P777/njjz9YunQpS5cuLb5ez549i1/3NKW9zs7OjqVLlz624eV/NWnShC1btvD777/Tr18/VqxYUeL5n3/+mfbt25d6fmnvLSQkhOPHj7Njxw5WrFjBV199RWHh4wP+K3t+cHAwKpWqxPYZWVlZxMYWffgNDQ196vt/kifd55EXXniBLVu2lPuazyIKnTrO2tyQj6f25IvX+pCdW8AbPx5h+a5rYlqwoDExSZl8uOEsS9afxtRIj1XzBzN7VAcM9Kp3QTJBQyICkT30K/GQ7KEfRARW+JLW1tY0bdoUAGNjY1xcXIp32z579izDhg0DinbTPnPmDFBUZDzaisDNzY38/PzifZtatWpVpgUQS3tdgwYN8PDweOJGmP/WsWNHDAwMiq8VHx9f/FxAQADJycklVh3+r9Le29mzZxkwYAB6eno4ODjg5OSEn59flZ9/5MiRx1qzBgwYwPHjxwE4evRoiWIvOjqa6dOnM2nSJCZNmsTt27cBuHHjBrNmzeKDDz5g/Pjx5OTksGDBguLWoGPHjgHQrl07vL29USqVpX9TK0AUOnXYIE9XNi4eTpdmDqw9cJMZ3x3CL6xqdosVhGe57BfF1C8PsP/SfV70asaGRcPp2MRe6lhCZSVGle/xcoqOjiYwMLB4b6bk5OTiYsTa2pqUlJTHzjl16hSNGzcusYFkddu7dy/dunUDQKVS8cMPPzBv3rynnlPae0tISMDW1rb4dTY2NsWF32effYa/v3+Fz/+327dv06xZsxKP9e3bl9OnTwNw/vx5evbsWfyclZUVq1atYuvWrXzxxRd89913xc/5+fnx+uuvs3PnTi5fvoy1tTXbt29nx44dxd8XuVyOo6MjDx48eOr3pbzEYOQ6yMJEn4UvdqZna2duB8fx9e+XiRYrGgsSyM4r4Ic/vTnhE8Y7L3Xhm1n92H/pAav3Xic3X6y9UyNZO5Tv8XLIzs5m8eLFLFy4EBMTkzKdExwczMqVK1m1alWl719Rhw4dIiAggF9++QWA3bt307179wqvNqx+wjINj1qXPvjgg0qd/29JSUlYWFiUeMzc3BwzMzOOHTtGo0aNilusoGiH+a+//pr79+8jl8tLjLdp0aIFDg5F/w+4ubnx448/snLlSnr06EG7du2KX2dpaUlCQsJjBVZliEKnjunRyomFL3bG2FCX1Xtv8OfZe2JtE0Fyd0MTeO3bg0wb3IaX+jSnrbstX2y7yL3wJKmjCeXl1AR1wxYluq/UDVuAU5NKXVapVLJ48WIGDx5M3759ix+3srIiMTERa2trEhMTsbS0LH4uLi6ORYsW8cknn5QYZ/IkhYWFTJ48GSja3XvWrFnlzrh69WouXLgAULwT+tWrV9m4cSO//PJLcYvSnTt3uHXrFrt37yY7OxulUomhoeFjO5iX9t5sbGyIi4srfl18fPwTu9gqe76+vn5xd9+/DRgwgK+//pqPPvqoxOPbt2/HysqK7du3o1KpSmzKaWhoWPznhg0bsmXLFi5evMhPP/1E586diwc55+fno6+v/9g9K0MUOnWEob4O81/wZJCnG/cjkli4+hJhsWlSxxKEYgVKFWsP3ORqQDTvTezGqnmD2HLMl60n7opB8TWJTA4DX0YdEVjUXWXtUFTkyCo+UkKtVvPpp5/i4uLCxIkTSzzXq1cvDhw4wNSpUzlw4ABeXl4AZGRk8OabbzJnzhzatGnzzHsoFIri4qSiZs+ezezZs4u/DgwMZNmyZaxYsQIrK6vixz/77LPiP+/fv5+AgIDHihwo/b316tWLDz/8kIkTJ5KQkEB4eHhxV15Vnu/i4kJERAQNGjQo8Xjv3r1JTEyka9euJbq8MjMzsbGxQS6Xc+DAgVIHHickJGBmZsZzzz2HkZERBw4cKH4uPDwcNze3J55XUWKMTh3g1sCCnxc+R/8OjdhyzJfZPxwRRY6gtW4Hx/Hq1wc4fesh04a04cc3BmJraSx1LKE8ZHJwbgbt+xf9txJFDhSNFTl06BDXr19nwoQJTJgwgYsXLwLw8ssvc/XqVUaPHs3Vq1d5+eWXAdi5cycRERGsX7+++Jzk5GQAVqxYwdChQ8nNzWXo0KGsXbv2ifct7XV+fn4MHTqUkydPsmzZMsaOHfvE83/88UdycnJ49913mTBhAgsXLizX+y7tvbm5udG/f3/Gjh3LvHnzWLRoUfHu4v8eo1OR8/+tR48e+Pj4PPa4sbExL7/8Mrq6uiUeHzNmDAcPHmTatGmEh4eXaMX5t6CgIKZOncqECRPYsGEDr7zyClDUVaavr1+mgeLlIQNq7Eclb29vPD09pY6h1YZ39eCN5zuSnp3HZ79d5HZw3LNPEgQt0a+9CwvGdEKlUrNs+yWu+FfNgNbaTBM/Fw8fPlzlv3wE7Zebm8vrr7/O+vXrn1gIVbXt27djbGzMyJEjK3R+YmIiQ4YMeexx0aJTSxnp6/Lh5B4sHNuZ28FxTP/2oChyhBrnpE8YM78/RHxqFsum9+HV59oglz99Sq8gCFXDwMCAGTNmPHFGliaYmJgwdOjQKr+uGKNTCznZmPHZq71xqGfCugM3+f2UH2K8sVBTRSdmMufHo8wb7cmkAa1o7lKfz7ZcICUzV+poglDrde3atdruNWLECI1cV7To1DJdmjuwesFgTA31WLj6BNtPiiJHqPnyCwr5dscVvtx+ieYNrVn79nM0FSsqC4JQBqLQqUUm9GvB56/2Jjoxg1nfH+JOSPyzTxKEGuSodwhzfjhCvrKQH+YMoE+7hlJHEgRBy4lCpxbQ11Xw4eQeTB/WjtM3w5i38hjxqdlSxxIEjQiJSWX28iPci0jioyk9mTakDc9YiV8QhDpMFDo1nIWJPsvnDKB324b8st+Hz7ZeJK9ArCgr1G5pWXm8veYkB68EMWVgK5a+3FPslVWLid3LtWf38tTU1FJf+zRnz55l06ZNAJw5c6Z4U9byOH/+fPHq0uUhCp0azLG+KavmD6aRvQUfbjjLH6f8pY4kCNVGWaji2x1X+Onv6/Ro5cTyOQOxMKnaFVWFilGr1SiTQ8kPv4YyOfSJWw6Uh9i9XHt2L68oLy8vpk6dChQVOuXd+VypVNKjRw/OnTtHbm75JiKIQqeGat7QmlXzBmFsoMvCn45zyS9S6kiCIIndZ+/x4YazuNiZs3LeIOzrlW0PJEEz1Go1uQEHyPXfR374ZXL995EbcKBSxY7YvVx7di+Hoo1VX3rppeKvf/vtt+LFFGfOnMl3333HK6+8wksvvVR83f379/P1119z+/Ztzp8/z4oVK5gwYQKRkZFERkYyd+5cJk+ezPTp0wkLCwPg448/Zvny5cyaNYuVK1cik8no0KED58+fL/V79iRaVeiYm5uza9cuAgIC8Pf3p0uXLlJH0krdWzry/ez+ZOTkM+fHIwSI/YCEOu6yXxRvrT6BmZE+q+YNwq2B5bNPEjSiMCWMwuSS3RKFySEUpoRVyfXF7uXS7l5eFjk5OWzYsIHFixfzv//9r8Rzbdq0oWfPnsybN4/t27fj6OjI559/zjvvvMNvv/3G/Pnz+eqrr4pfHx4ezk8//cSbb74JQLNmzbh161a58mjVOjo//vgjR44c4cUXX0RXVxcjIyOpI2mdQZ6uvDOuC4HhSSxZf4a0rDypIwmCVvB/mMi8lcf4emZffpw7gA9/PcvNILFIZnVTZT55cTlVVgJYNarUtcXu5UWk3L28LB516bVv356srCwyMjJKfW12dja+vr68++67xY8VFBQU/7lfv34lVmW2srIq9wKGWlPomJqa0qtXr+I+vIKCAtLSxH5M/zaye2MWjOmEd2A0H204S26+GHQsCP/2MC6NOT8e5ZtZfflyZl8+2XRedOtWM7lJ/Sc/bvzkx8tK7F6uPbuXKxSKEsXSf1/z36LpaV18KpUKExOTUjdU/e9+WXl5eeXe3Vxruq5cXV1JSEhg48aN+Pj4sG7dOtGi8y9jezdjwZhOXLwbyfvrz4giRxBKkZiWzbyVxwiOSuGTab3o0cpJ6kh1isLSBYWVa8nHrFxRWLpU+Jpl2b0cqJLdy7dv316hIgeKdi9/dA34/93Lv/vuu8d2Lz9w4AD79u1j/vz5PPfcc0/dvfy/761Xr14cP36c/Px8oqKinrl7eUXPf7R7+X/Vq1eP5ORkUlNTyc/PLy7uHjl+/DgAt27dwsTE5LHWN2NjY7KysoCibR8aNGjAiRMngKK/6/v37z92z0cqsru51hQ6Ojo6tG/fnjVr1hQ3d/27KeuR6dOn4+3tjbe3d53ZZG7KoFa8PrIDp3zCWLrxLAVKldSRBEGrZWTn887PJ7kXnsTHL/fEq42z1JHqDJlMhkGzYRg0H4Few64YNB+BQbNhzxy4+zRi93Lt2L28sLAQXV1ddHR0eO2115g2bRpvvvkmDRuWXLjT1NSUV155hWXLlvHhhx8+dt2BAweydetWJk6cSGRkJJ9++il79+5lwoQJvPTSS5w9e7bU78mNGzfo0aNHub6PWrN7ua2tLVeuXKFRo6I+3B49evDuu+8Wjxh/krqwe/m0wa2ZMqg1h68G8+2OK6jEfg6CUGaG+jp8OaMvLRpa88W2S5y6GSZ1JI0Tu5cLVeXfu5enp6czceJEDh069NRzZs6cyfz582nevHmV50lKSuKDDz5gzZo1T3xe63cvj4uLIyIigsaNGwNFA5AeVaV11YR+LZgyqDUHrwTxzY7LosgRhHLKyVOy+JdT3AmJZ8mkbvRv7yJ1JEGoMR7tXn7kyBGmT5/OnDlzJM0TGxvLggULyn2e1gxGBpg7dy7btm1DT0+PkJAQpk2bJnUkybzQqynTh7Xj+PVQvt95VWzMKQgVlJuv5L11p/lieh/endCN7DylGKAsCGX0aPfyoUOHlun1FVm5uKyeNI6oLLSmRQeK+mE9PT1p06YNzz//fIWXmq7phnfz4I3nO3L29kO+/P2SaMkRhErKKyjkg1/PcD8ymaUv96S9R8Wm9QqCUPNoVaEjQP8OjVj4Ymcu+UXy2W8XUalEkSMIVSEnT8nitaeITEjns1e9aOZcT+pIgiBUA1HoaJGOTexZPL4rPvdj+HjTOZSFYnaVIFSljOx83vnlFMnpuXw5oy+N7C2kjiQIgoaJQkdLuDtY8snUXjyMTeOjjefEFHJB0JDk9BzeWnOCvAIlX83og7W5WK9LEGozUehoAVtLY76c0Zf07DwWrz1FVm7Bs08SBKHC4lKyeHftaYwMdFk2vTeG+lo1L0P4j9jYWGbNmsWLL77I2LFj+f3334ufS0tLY86cOYwePZo5c+aQnp4OFK1IPHnyZMaNG8fkyZPx9vYGiqZML1iwgDFjxjB27FhWrlxZ6n0DAgIYN24czz//PN9++23xasA+Pj5MmjSJLl26cPLkyVLP37ZtG2PHjmX8+PG8/vrrxMTEFD83d+5c+vTpU7yH05OU9t4ANm7cyPPPP88LL7zA5cuXy3V+amoqs2bNolevXnz99del3r+2EIWOxEyN9PhqZl/0dOS8u/Y0Sek5UkcShDohJCaVjzedp5GdBUtf7olcXvEF7QTN0tHRYcGCBezatYuNGzeye/duQkKKNg7dvHkznp6e7NmzB09PTzZv3gyAhYUF33//PX/88QdLly5l6dKlxdebNGkSu3fvZtu2bdy5c6d48cH/+vLLL1myZAl79uwhPDycS5cuAWBnZ8fSpUuL93QqTZMmTdiyZQu///47/fr1Y8WKFcXPTZ48mU8++eSp55f23kJCQjh+/Dg7duxgxYoVfPXVVxQWPr5afmnn6+vrM2vWLObPn//U+9cWotCRkEIu45OpvbC3MuH9X8/wME7s7SUI1el6YAzf77pG52YOzH+hdi8+WpNZW1vTtGlToGj7ABcXl+KNHc+ePVu8sOywYcM4c+YMUFRk1K9ftL+Wm5sb+fn55OfnY2BgQMeOHQHQ1dWlSZMmxMfHP3bPxMREsrKyaN26NTKZjKFDhxav2NugQQM8PDyeudpzx44dMTAwAKBVq1Yl7tOpUyeMjY2fen5p7+3s2bMMGDAAPT09HBwccHJyws/Pr8znGxoa0rZtW0l3c69OotCR0BvPd6Sdhx3f/HEZ35Dy7cYqCELVOHQ1iG3H7zKiW2PG9m4mdRzhGaKjowkMDCxeUyU5Obl41WZra2tSUlIeO+fUqVM0btz4sV/sGRkZnD9//okrScfHx2NjY1P8tY2NTbl3zf63vXv30q1bt3KdU9p7S0hIwNbW9onZ/r0FRFm+N3WB6JiWyLCu7ozq0YTfT/pxwidM6jiCUKf9evgWDvVNmTG8HcHRKdy4Hyt1JOEJsrOzWbx4MQsXLnxso8jSBAcHs3LlSlatWlXicaVSyfvvv89LL730xJ3N1VW4ftmhQ4cICAiossX0npTtUevSBx98UCX3qE1Ei44EWrnWZ/4LnbjiH8X6g7ekjiMIdZ5aDV/9fomw2DQ+nNITO6undykI1U+pVLJ48WIGDx5M3759ix+3srIiMTERKOpusrS0LH4uLi6ORYsW8cknnzxWzHzxxRc4OzszYcIEoGjDykebf/7888/Y2tqW6GqKj48v7gorzerVq4uv8cjVq1fZuHEj3333Xbm7ikp7bzY2NsTFxZXI9qS9yJ72valLRKFTzWwsjPhkqhfRiRl89tsFseqxIGiJ3PxCPtpwFrkMPn3FC33dx3dzFqShVqv59NNPcXFxYeLEiSWe69WrFwcOHADgwIEDeHl5AUXdUm+++SZz5syhTZs2Jc5Zs2YNmZmZJXYTVygUbN++ne3btzNr1iysra0xMjLC19cXtVrNwYMHi69dmtmzZxdfAyAwMJBly5bx3XffYWVlVe73Xdp769WrF8ePHyc/P5+oqCjCw8OfuD1CaefXNVqze3lF1LTdy3UUcn58YwDOtubMXn6EiIT0Z58kCEK16tS0Acum9+HUzTA+3/rk2TjarDbuXn7r1i2mT5+Ou7t7cRfNnDlz6N69O6mpqbz33nvExcVha2vLl19+ibm5Ob/++iubNm3Cycmp+DqrVq2ioKCAYcOG4eLigq6uLgBjx45l1KhRj93X39+fTz75hLy8PLp168Y777yDTCbDz8+PRYsWkZ6ejr6+PlZWVuzcufOx82fPnk1wcDD16hWtwm1nZ8f3338PwPTp0wkLCyMnJwdzc3M++OCD4n2lHintvQFs2LCBffv2oVAoWLhwId27dweKxuiMHj2a5s2bP/X8ESNGkJWVRUFBAaampqxcuRJXV9fK/DVJrrTdy0WhU41mj+zAi72bsXTjWc7diZA6jiAIpZjYvwWvDW3Hij3e/HU+UOo45VIbCx1BKIvSCh3RdVVNurd05MXezdhz7p4ocgRBy2074cclv0hmjWiPu0PdHNcgCLWFKHSqga2lMYvHdyUwIomf9/lIHUcQhDL46vfLpGXm8eGUnhjoiQmqglBTiUJHwxRyGUtf7olMJuOTzecpEBt1CkKNkJ6VxxfbLuJobcq80TWni1wQhJJEoaNhUwa2ollDa77ZcYWYpEyp4wiCUA63guLYeuIuQzq70a+9i9RxBEGoAFHoaFAz53pM7N+So97BnLsdLnUcQRAqYPPRO/iGxPPmi52wtRTr6whCTSMKHQ3R11Xw7sTuJKbnsHLPdanjCIJQQSqVms+3XkSGjHfGdZE6jiAI5SQKHQ2ZMawdzjZmfLX9Elm5BVLHEQShEuJSsvh5vw8dGtszvJuH1HEEYO3atfz222/8/PPPXL169bHnb9y4wZtvvilBMkHbiKkEGtDew47RvZqy+2wAN4Pinn2CIAhab/+lB/Rq7cys4e25FhBNXEqW1JG0mBoz83yMjArIztYlPU2PomXbqt6sWbM0cl2h9hCFThXT11Xw1kudCY9PZ53Yx0oQapVvd1zh10VDWTSuC2//fBKxg8uTqHFzS8XCMq/4kdQUfYKDLahssbNhwwYOHjyIra0tlpaWNG3alI8//piePXvSr18/Ll26xPfff4+FhQVNmzYFQKVSMWbMGH799VcsLS1RqVS88MILbNy4kR9++AEDAwPCwsKIiYnho48+4uDBg/j6+tKiRQs+/vhjAI4ePcrGjRtRq9X06NGDuXPnUlhYyKeffkpAQAAymYwRI0aU2ONK0B6i0KliUwa1okE9UxasOk5+QaHUcQRBqEJxKVms2evD2y91YVgXD/ZffiB1JK1jZp5fosgBsLDMw8w8n/Q0/QpfNyAggGPHjrFt2zaUSiWTJ08uLmYA8vLy+Pzzz1mzZg1OTk4sWbIEALlczpAhQzh8+DATJkzg2rVreHh4YGFhAUB6ejpr1qzh3LlzvPXWW6xfvx5XV1defvllAgMDsbKyYuXKlfz222+Ympoyd+5czpw5g62tLQkJCezYsQMo2ltL0E5ijE4VcrW34KXezTl8NYjbwaLLShBqo4NXgrj5IJbpw9piYVLxX9y1lZHRk8cklvZ4Wd28eZPevXtjYGCAiYkJvXr1KvF8WFgYDg4OODs7I5PJGDx4cPFzw4cP59ChQwDs27eP4cOHFz/Xs2fROmdubm5YWVnh7u6OXC7H1dWVmJgY/P396dChA5aWlujo6DB48GBu3ryJg4MDUVFRfPPNN1y6dAljYzEjT1uJQqeKyGUy3hrbmfTsfLH6sSDUcj/svoaBng4zh7eXOorWyc7WLdfj5fFoQ8/ysrOzw8rKCm9vb/z8/OjWrVvxc3p6ekBRy8+jTT4f3auwsBB1Kf2TZmZmbN++nQ4dOrBr1y4+++yzCmUTNE8UOlVkRDcPmrvUZ/Xf10nPzpc6jiAIGhQen87O0wEM7uRGa1cbqeNolfQ0PVJTSrZ0pabo/zMgueLat2/PmTNnyM3NJSsri/Pnz5d43sXFhejoaCIjIwE4duxYiedHjRrFRx99RP/+/VEoFGW+b8uWLfHx8SE1NZXCwkKOHj1K+/btSU1NRaVS0bdvX2bNmsW9e/cq9f4EzRFjdKqAubE+rw5ty/XAGE74hEkdRxCEavDbcV/6dXDhzRc78do3BylUiZHJRWQEB1tU+ayrpk2bMmDAACZOnIi9vT1t27Yt8by+vj5LlixhwYIFWFhY0KZNG0JCQoqf79WrF//73/9KdFuVhbW1NXPmzGHWrFmo1Wq6d++Ol5cX9+/f53//+x8qVdG2PnPmzKnU+xM0RwbU2H+d3t7eeHpKvwfNgjGdGNbFnVe/OcjDuDSp4wiCUE26tnDgi9f68PM+H3ac9pc6DqCZn4uHDx/G2tq6Sq9Z3fz9/Vm+fDnr1q2TOoqgIYmJiQwZMuSxx0XXVSW52lswrKs7ey/eF0WOINQxl/2iuOwXyaQBLTE3FgOTtdWmTZtYvHixaHWpo0ShU0mzR3UgM6eATUfvSB1FEAQJ/LzPB0M9HV4e1FrqKEIppk6dyv79+x/r7hLqBq0rdORyOT4+Puzfv1/qKM/UvaUjHRrbs+nIbTLEAGRBqJPC49M5cPkBI7p54FTfTOo4giD8h9YNRp4/fz4BAQGYmWn3DwyFXMbM4e0Ji01l3yWxaFhNJFcoUBUWLero1rEd7p07Yl7fGmMrCwxNTdEzMuSHl6YBMOKdeXQcPoTCwkJUSiXK/AIyU1JYOWkGAO2GDMDKsQHp8QkkR8WQ8DCC9IREyd6bUL02Hb1D/46NmDm8HR9sOCt1nCr3aMCtIGiz0v4/1apCx8HBgaFDh/L555+zcOFCqeM81SBPN5xszHh//RlUYraF1tM3MsK5VXMcmjXBoVljGjTxwNrZkY97DyUnPYPGXTvR97UpZCYlk5mcQnZ6Bmmxcch1FKiUhUT43UNHTw+5QoFCRwcdvZJrgrTo3YN2zw0s8VhCWDhfDn8JAI8unuRmZhF1LxCVUqyYXdukZuax/YQf04e1o42bba1bMDQoKAgrKyt0dLTqV4YgFFMqlQQFBT3xOa2adbVr1y6WLVuGqakpb7/99jOnAUo160pXIee390eSnJ7D7B+OVPv9hWfTMzTEvVMHHt72JSs1ja5jn2fMh4sASI6OISYwiNjgUM5u3k5Wahq6BvoUFiiLW3gqQkdfH/P61tRzaoC1sxNyhZwL23cDsGjv79i6upCfk8vDO3cJ9bmN35kLRPqLtTdqCz1dBb+9N4L41GzmrjgqWQ5N/Fy0tLTk448/Ll41WBC0iUqlIigoiI8//piUlJTHntea8nzo0KHEx8fj4+ODl5dXqa+bPn06M2YUdRdINd1xeDcPbC2N+fqPy5LcX3gyC1sbWvXvTYs+PWnUvg06urr8/v6nXN93CL/T50mKiCLS/x7ZaemPnVuQm/eEK5aPMi+PpMgokiKjuH/Zu8Rza2cuwLl1Cxq1bU2j9q3pP2MqFva27PjwcwBa9vXiwVVv8rKyK51DkEZ+QSFbT9xl4Yud8Wxij3dgjNSRqkxKSgrz58+XOoYgVIjWtOh88cUXTJ48GaVSiYGBAWZmZuzZs4fJkyeXeo4ULToGegq2fTCKh7FpLFx9olrvLZTOrL41S08VDWCPeRDMvfOXCbx0lRCf2xQWVG6PHU0wNDNF39CQ1Lh47DzceGfPVpT5+QRevMr1A0fwP3MBZb4Y4F7T6Cjk/LZkBMnpOcz5UZpWHW1ZX0wQtIla2w4vLy/1/v37n/k6b2/vas82vl8L9enlk9QtXKwl/z7V1UMmk6mb9uyqnvbjl+qJX31S/HjnF0aorRs6SZ6vIu/HpW1r9Yh35qk/PLFX/Z3vZfWnF4+qnVo2lzybOMp/DOvqrj69fJK6U9MGktxfip+L4hCHNh9a03VVE+jrKnjRqxnXAqLxCxMzaqqbgYkxnUePoNtLo7F2diQ9MYlre/5/GYKrf+6TMF3FqdVqwm7dIezWHfZ/twqPzh1o99xAYoOCAWjVzwuZXI7vybOoxewXrXfkWggT+7dk6uDWXLsXLXUcQajztLLQOXv2LGfPat8UzSGd3LA0NWDbibtSR6mT+kybRP8ZUwm5cYvDK37G9+RZCpVKqWNVKbVKxf3L3iXG+HQd+zxNunUm4WEEpzdu5fq+w1rZHScUURaq+O3YXd4Z14VOTRuIYkcQJCaGz5eRXC5jbJ/m3A1N4E5IvNRx6gRjSwuGvz0Xjy5F4w3Ob9vJ92Nf5qepr3Pr6MlaV+SUZt3rC9n05nvkZmYy9uP3WHJ4N636eUkdS3iKY9dDSEjN4qU+zaSOIgh1nla26GijPm0bYl/PhFV/XZc6Sq2nZ2hIn2kT6TVlHHoGBmSlpPHgijeZySlkJqdIHa/aqVUqfE+cwffEGTy6eDJ4zvTiafCP1vkRtIuyUMWf5wKZNaI9Ho5WPIhMljqSINRZotApo/H9WhAWm8pl/0ipo9RqrQf0YdTiNzG3rc/tY6c4smot8aEPpY6lNR5c8ebBlf/v1how8xXcPdvz91fLiQq4L2Ey4b8OXH7A5IEtGdu7GZ9vvSh1HEGos0TXVRl4NrHHrYElf5zyR62WOk3tZlrPirSEBFZOmsGWt94XRc4zJEVEUd/FmQV/bGT0+29jYGIsdSThH1m5BRy8HESftg2pb2EkdRxBqLNEoVMGo3o2ITk9h5M+YVJHqXV09PQY/MYMOgwbDMClnX+xYsJrhN32lThZzXB93yG+HDaWC9t20fXFUSz6+3fcPNtLHUv4x5/nila+fqFXU4mTCELdJQqdZ7CzMqZLMwcOXAlCWSim9lalBk08WPDHBgbMnIZji6JfBGqVCrVoNiuX3Mws9n79AysmTic9MZHstDSpIwn/iE/N5sztcIZ2ccdATyF1HEGok0Sh8wwjujVGrVazX+xQXmVkMhl9XpnE/N9/xcjcjHWvv8ner36QOlaNF+EXwA/jXiHmftH6O0PmzsS1YzuJUwl7LwZiYqhHn7YuUkcRhDpJFDpPoasj57nObly8G0limtiDqKq4ebZn2Jtz8Dt9nm9HT+LehStSR6p1DEyMaTOwL6+vX0m/115GJpNJHanO8g1JICw2leHdPKSOIgh1kih0nqJvOxfMTQz4+0Kg1FFqBUMzUwCCrt1g1ZSZbHnr/SdusClUXm5mFstfmsbtoyd5bv4sXv3pWwzNzKSOVWftv/SAZg2tcXewlDqKINQ5otB5iqFd3AmPS+NmUJzUUWo8z5HP8f6RPTi1KFpALfTmHYkT1X552dlsXbyU3Z9+jUcXT2au+1G07Ejk2PVQcvOVolVHECQg1tEpRQNrE1q52vDLfh+po9RoMrmcYW/OoffUCdy/4k1iRJTUkeqcyzv/IuZ+MIZmpmKgt0Qyc/I5ffMh/ds34ud9PuTk1Y1VvQVBG4gWnVIM6uhKoUrFiRuhUkepsRS6ukz+5lN6T53Ahd93s+71N8lJF11VUgi7dYeAc0WL1nV+YQTthw2SOFHdc/DKA4wMdOnZyknqKIJQp4gWnSeQyWCgpys+92NJTMuROk6N1eWFEbQZ2Jd936zg7JbfpY4jUDTjre2gfjTu2gkTK0vObflD6kh1hl9YItFJGfTv0Ihj18UHKEGoLqJF5wlau9piZ2XCMe8QqaPUaJd27OGnabNFkaNF1Go16+e8za2jJxn5znwGzHpF6kh1yskbYbRvbIelqYHUUQShzhCFzhMM8nQlKzef874RUkepcYzMzZj245dYNrBDrVYTcv2m1JGE/ygsKGDroo/w3nuQwXOmM2j2a1JHqjNO3AhFIZfTt52L1FEEoc4QXVf/oaOQ07O1E+fvRJBXIHaFLg8DUxNmrVuJjWtDrJ2dSImOlTqSUAq1SsWOj75AVaiiIC9P6jh1Rnh8OvcjkujfoVHx9hCCIGiWKHT+o72HHSaGepy9HS51lBpFz9CA1376Dlv3RmyYu6jEDtuCdlKrVOxc+kXx14ZmpuSkZ0iYqG44cSOM2aM64FjflMgE8f0WBE0TXVf/0auNM1m5+dwIjJE6So0h11Hw8vfLaNi6BdsWLyXwoljpuKaxc3dlyaHdtHtuoNRRar3Ttx4C0EPMvhKEaiEKnX+Ry2X0aOnIFb8oCsQGnmWmb2SMgakxuz75ijvHT0sdR6iAxPBIou8HMe7T92nUvo3UcWq1xLRsAiOS6N5SFDqCUB1EofMvrV1tMDcx4Nwd0W1VHjnp6fz08utc+2u/1FGEClLm57NpwXskR8Uw7YcvqefkKHWkWu2ibyTNG1qL2VeCUA1EofMvvVo7k5uv5Nq9aKmj1AjNvXrw2urvMDA1QVUoBm7XdDnp6ayf/RYAr6z4Cn0jI4kT1V4X7kYgl8vo1kIUlIKgaaLQ+ZfOzRtw434sufnil/azWDs7MmHZUkzqWaHML5A6jlBFkiKj2PL2B0T43UOlEv8ONCU0JpXopAy6txSFjiBomph19Q8Ha1Ma1DNl5+kAqaNoPV0Dfab+8CUqpZLNC95DKaYn1ypB124QdO2G1DFqvYt3IxnZrTEGegrx4UoQNEgUOv/o1KwBAN6i2+qZRrwzH1u3Rqyb9SYpMXVzrRwduQwbY33qGephZaiLhYEuOnIZcpkMuQxylSrS85Rk5CtJzsknMj2H/MKataGmtbMjk775lF1LlxF1777UcWqdq/5RvOjVjFauNnjfE7M8BUFTRKHzD88m9kQmpBOdlCl1FK1maGZGc6/unNm0jfuXr0kdp1royGW0tDGjvZ057lbGuFkZ09DcEB152Xt+VWo1cZl5PEzLxj8hg9tx6dyNTydNi3exzk5Lx6xePSZ9/T+WvzSV/JxcqSPVKndDE8hXFtKhsb0odARBg0ShA+gq5LR1t+PItWCpo2i9nPR0vnthMnlZ2VJH0ah6hrr0d61PT+d6tLUzx0BHAUBkeg7ByVmcCUskKj2XpJx8ErPzScsroKBQjUpddBjoKDDT18FUT4f6xno4mRnhbG5IIwsjprRxRkcuAyAgIYOzD5M48zCRoOQsKd/yY7LT0tm+5BNmrlvB0AWz+WvZ91JHqlXyCgq5G5pAh8Z2UkcRhFpNFDpAy0b1MdTXEbOtnqHtoH7cOXmG7LR0qaNUmgzo5mRFU2sT7iVmcikiGV2FjIGuNjznYUsHewsUchkhKVn8dS+G69Gp3IxNI72MLTBpeUrisp48dslAR05za1Pa2pnT3dmKGR0aMqujC2Gp2fx1L4aDD+JIzdWOAd5B125wYfsuek16iVtHTxLqc1vqSLWKz/0YXhvaDgsTfVIzxVg3QdAEUegAbdxtKVSpuBMcL3UUrdV2UD8mf/sZOz9extU/90kdp1JkwLcDWuDlYl38WFhqNmb6OlgZ6vEwLZtfbz7keEgCoalV33KVq1ThE5uGT2waG26FY2WoSy/negxrbMebXdx4w7MRx0MS2HQ7nJAU6VvODq/4mRa9e9Bz4lhR6FSxG/djeW0otHO3K14xWRCEqiUKHaBVo/oER6eSnacdn6K1jYGJMSMXLyD8rj/efx+UOk6ldXOyKlHkALhYGOEbl8aSUwFcj06t1jzJOQX8HRjL34GxuFoaMaqJPaOa2jPY3YbToYms9XlIcIp03Vr5Obn8Mn0+KbF1c+C5Jt2PSCYzJ5/2jUWhIwiaojXr6Dg6OnLq1Cn8/f25e/cu8+bNq5b7KuQymjWsz90Q0ZpTmoGzX8OknhV/fvpNrVgYsHl90yc+fj48udqLnP8KScnm+yvBDP/9ChtuhtPZ0ZLtozuwuLs75vrSfS5JioxCpSzEwMQYU+t6kuWobVRqNXdDE2jZqL7UUQSh1tKaQkepVPLWW2/RvHlzunTpwpw5c2jWrJnG7+vuYIWhvg6+oaLQeRI7d1d6jB/DlV1/E+l/T+o4ldbB3oJRTZ48+DNQi2bcpeUp+flGGMN/v8pO/yieb9qAP8d2YmQp2auDXKHgzR2beOGDdyTLUBv5hyXiYmeBsYGu1FEEoVbSmkInNjaWmzdvApCZmUlAQAAODg4av28r16JPUndDEzR+r5pIrlAQdPU6h1f+InWUStGVy5jf2ZVfhrUhv1CNb1zJAdVnwxK5FJEsUbrSZeQr+e5yMBP33CA4JYsPezXhuwEtsJTgl6KqsJBrfx2gVT8vXNq2rvb711Z+D4t+9jR1Fi1lgqApam07GjZsqH748KHa1NT0qa/z9vau9L0+mdpLve2DkZK/Z3Fo7nCxMFJvH91BfX26l3pxd3e1vkKuloG6u5OV+pW2zuruTlZqmRbkfNYhA/WElg7qi9N6qo9O7Kru1MCi2jPoGuirPzq5T/3Gll8k/37UlsNIX1d98ruJ6ikDW1XJ9ari56I4xFGbDq1p0XnE2NiYP//8kwULFpCRkfHY89OnT8fb2xtvb2+sra2fcIXyaexkRcDDpEpfpzbymjIe03pWUseolG6OVmwa2Q5rIz3mH/Hlq4tB5BWqUAMXI5LZcCucixHJqKUOWgZqYPvdKKb8fYPU3AJWDGnNSy003+r5bwW5eRxdvZ5G7VrT3KtHtd67tsrOKyAsNpVmDSv/80wQhMdpVaGjo6PDn3/+ybZt2/jrr7+e+Jp169bh6emJp6cniYmJlbqfmZEedlYmBEVpX5eF1Jr26MKId+bR3Ku71FEqbHxLB5YPaklURi6T//LhohZ2TVVEcEo20/bd5EJ4Eu90c+f9no1RyKrv/t57D5IUGYVHl47Vd9Nazv9hIs1FoSMIGqFV08t//fVXAgICWL58ebXcz92hqLUiKCqlWu5Xkwyc9SrJUTFc33dY6igVMsezEdPaOnM6LJGPTgeQo1RJHalKZRcU8s5xP2Z1dOHVdg0x09Ph/dMBKFWab5tSKQv5YdwrtWLhSG1xPzKZYV09sLEwIl4DazcJQl2mNS063bt3Z8qUKfTt25ebN29y8+ZNhgwZotF7ujlYAvAgsnZ80q8qbp7tadimJad+/Y1CpfbuxVSaeZ1cmdbWmd3+0Sw67lfripxH1MCa62F8dzmIfq71+XZAC/QV1fNP+lGRY2hmVi33q+1CY1IBcG1gKW0QQaiFtKZF5+LFi8hk1dj+Dng4WJGQmkVaKUv111VeU8aTkZSM996atzjg/M6uTG7txE6/KL6+FCR1nGrx+90ocpUq3uvhwed9m7H4hB/VsVF6s57dmPrDMn4Y/yox9+vG91pTHhU6jewtuOIfJW0YQahltKZFRwoejpY8EN1WJcgVCvKzs7mwfRfK/Hyp45TLuBYOda7IeeSvezF8eymI3i7WLOruUS33DLvtS6GykF6TXqqW+9VmWbkFxCZn4mZvIXUUQah16myhI5fLcKxvVvxJSiiiKixk6+KlnFi7Seoo5eLVsB4Lu7pxKjSBb+pYkfPITv9oNt4K54VmDXi5jZPG75eTnsH1fYdo99wA0YVVBUJiUmkkCh1BqHJ1ttCxszJGRyEnIl4MqHxER08Pm0YNpY5Rbm6WRnzWpxkBCRl8ePpejZgqrik/eYdyNCie2R0b0bGBhcbvd2X3XnT19Wk/dKDG71XbhUSn4Gxrjk41jbMShLqizv6Lcqpf9Ak0KvHxtXrqqlb9vFi87w+cWzWXOkqZ6SvkfNG3OVkFhSw85kdeoeYHHssNTDGwb46xazfMmg3CvOUwzJoNxMTDCwP75iiMpV176LPzgTxMy+bzPs2wNtLT6L2iAx8Q4X+PTqOGafQ+dUF4fDo6Cjl2VsZSRxGEWkVrBiNXN8d/Ch3RovP/PEc+R3JUDBF3A6SOUmbzO7viZmXMG4fvkJSjgTFFMjmGDq0xdumMkWMbDOybo2P07JkxhTnp5MYGkBPtS2bwRbKj7oCqemaw5ShVLD7hz5ZR7fnEqylzDt/R6P32fPYNWSlpGr1HXRDzz15rDeqZEpkgPoAJQlWpw4WOKRnZeWLG1T9MrCxx79yR0xu3oVbXjM6fbo5WjG3hwNY7EVyJrNpB5fq2TbBs9wJmTfqjY1IPtVpFXnwQGYGnyEsKoyAlAmVWMqq8TNSFBcgUush0DdA1tUHXzB79+m4Y2DejXpeXse7+GoW56aQHHCf11l/kRPtWadYnCU3N5oerwbzXozHDG9uy/36cxu4V7uuvsWvXJdHFhY6JxEkEoXapw4WOmfjU9C+t+vdGoaPDrSPHpY5SJvoKOYu6uxOSksVP3qFVdl0Tt+5Yd5+OkVM7VAU5ZDw4R8a9k2SGXEKV9+z/X3Jj/Ep8LdczxrhRF0wb98a8xXNYtnuBnBh/Ei+uIyPwNGhwRNGegBgGudnwZhc3LkWkaKbF6x+uHdrSZmBf/lr2vcbuUdslp+eQm6/EXhQ6glCl6uwYHXsrk+JPUAK07t+buJAwYu4HSx2lTKa0ccLRzJCvLwVRUAWrARvYN8dlyiacx61Gx9SG2OPfcH/FAKL+WkR6wNEyFTlPosrPIiPwJNH7P+T+j32JPvgJcj1jnMYsp9Grf2Do2K7S2UujBj4/fx99hYK5nRpp7D4ANq4u9JjwInburhq9T20Xk5RJA2tTqWMIQq1SZ1t0rM0NSbwrllp/ZMvbH2JhZyN1jDKxN9FnahtnjgbFcz06tVLXkukaYtN7LlYdx6HMTibm8Gek3PpLI+NpVPnZpN7aQ+rtvzFrPhjbPvNo9PImUm/vJfb416jyqr7wfpiWww6/KCa1duS3O5EEp2RV+T0AAs5dBKBZz67EBoVo5B51QXRSpmjREYQqVidbdEyN9NDX0yExTRQ6j+Skp9eY1W1faVc0Bf7Ha5VrfdK3dsN12jasPMeT4rOL4J9HkeKzS/ODhtUq0v0OEfTL8yRcXId5q6G4vrYTQ4fWGrndptvhZOYrNdqqkxaXQNS9+zTrVXM3gdUGcSmZ2FqKWVeCUJXqZKFjbW4EIAqdf/SY8CJdxz4vdYwysTPRZ5iHLXsDY4jPqviYE9MmfWn0yjYUhuaEb59J7NFlGmlReRp1QQ4JZ1YRtmUaqNW4TN6Aeauqn6adnqdk8+0IejjXo7kGu0UCzl/GpW0r9I2MNHaP2i4lIxcTQz10derkj2ZB0Ihy/2syMjJCLq/Z/witzQ0BSEzLkTiJdug1+SWadOssdYwyebTi7+bbERW+hmWHsTi+8B258Q8IXj+WrLBrVRWvQnKi7hDy6ziyw31wGPE51j1mVPk9dvtHk5mvZEIrhyq/9iNB126QGB5ZY7pAtVFyei4AliYGEicRhNrjmRWLTCZj/PjxHDhwgLi4OO7du0dMTAx3797l66+/xt3dvTpyVqn6/7ToJKSKFh1LezvqOToQdO2G1FGeyVRPh+GN7TjwII64Ci4LYNlxHPaD3yfzwVkebp1OYVZSFaesGFVeBg//mE3qnX3YeM3BusfMKr1+VkEhf9+Lob9rfWyN9av02o88uOLN1yPHExcSppHr1wUpmUUfvixNDSVOIgi1xzMLndOnT+Pm5sZ7772HnZ0dzs7O2Nra0rNnT65cucKXX37JxIkTqyNrlXn0QyQ5Q7TouHZoC0Dw9ZvSBimDwe42GOgo2O0fXaHzzVsOw37Qe6QHniLiz7dRK3OrOGElqZRE7/+I1Nt7sfGajZXnhCq9/A6/KGTIGNXUvkqvK1SdlIyi/yetzESLjiBUlWfOuurfvz9K5eODM1NSUtizZw979uxBR6dmTd4yMdQjN19JgVLz2wVoO+dWzcnLzq4RM2VGNrHjXmIGgRVYFsDQoTX2Q5eSFXqVqL8WVdsqxeWnJvrgx8j1TbDt/zZ5SWFkhVyqkivHZOZxPTqVIe42/HIjrEqu+V/dXhpNn2mT+OK5MahV4t9XeSWn/9OiYyJadAShqjyzRedRkdOhQwf27NnDjRs3uH37Nnfu3OH27dslXlNTmBjqkqnBxdNqEh09PcJu3tH6X0qNrYxpam3K3sDYcp+rMLLC8YXvUGbEEbHnbdSFBRpIWIXUKqL2vU9ewgMcn/8KHTO7Krv04aA4HM0MaVFfM4OS87JzsHKwx8bFWSPXr+1S/+mStTDRTPeiINRFZW6K2bZtG++88w6+vr6otPyX4rOYGOqJQucfuz75UuoIZeLlYo1KreZ4SEK5z7Uf8j4KQwtCN05AlVsz9jZTF+QQsXshrq/twmH4pzzcNoOqWEX5dFgi7ylVDHKzwU8DK4NH+t8DwKF5EzFWpwLyCwpRFqow0teVOoog1BplLnQSEhLYv3+/JrNUm6JCR8s/1Qsl9HCy4m58Bqm55ft7M2s2ELOm/Yk7uZy8+AcaSqcZBalRxB3/mgbDPsGy/ZiiNX4qKaugkBsxqXR2fPbGpBWR8DCcwgIlto1cNHL9uiAnrwBD/Zo1HEAQtFmZ54kvXbqUdevWMW7cOJ5//vnioyYSLTpFmnTvwtzf1mLZoOq6RjTBylCXFjZmXIwo3wwpmY4+tv0WkhMbQNLVLRpKp1mpt/8mK+wa9b3mIDcwq5JrXotKwc3SGGsjvSq53r+plIUkRUZh06hhlV+7rsjOU2IoWnQEocqU+WPDtGnTaNq0Kbq6usVdV2q1mr/++ktj4TTF2ECX6ESxoadDUw9c2rYiO027u3M8G1gAcDEiuVznWXlOQNfcnqh9H4C65na3xh7/BtdX/6B+99eIO1n5TTOvRRXt9N7ZwZKDD6p+V/MbB4+izBMfJCoqJ68AIwPRoiMIVaXM/5ratGlD69aaWaK+uunqyMlXFkodQ3JWDg3ISEomL0u71xNqZm1KrrKQB+WYbSXT0ade5ylkBl8gO/y6BtNpXl78fdL8jmDZ/kUSLq6v9DijB8lZZOUraVHfVCOFzolfNlb5NeuSnDylGKMjCFWozF1XV65coVmzZprMUm10FHKUhTX3E35VsbS3IyWm/LOYqlsza1MeJGVRWI6xuOathqFjbEXipdrxSzfpykbkekZYdRhb6WupgaDkLDw0uHmkQkcHWQ1fQV0qYoyOIFStMv8k6tGjB7du3eLevXuPTS+vaUShU8TC3pbUmKr/RF+VZEATaxMCytnVaNl2NLlx92p8a84jefEPyAq9ikWbUVVyvcCkTDysNLN5ZIvePfjyxlns3DW3iWhtVqhSI5PJpI4hCLVGmT82DB48WJM5qpUodIrEBD4g0j9Q6hhPVd9YDxM9HYKSs8p8jl49FwwbtCT2+DcaTFb9Un334TDicwwd25ITeatS1wpOycJET4f6RnokZFfteJqcjEzkcjkmlpqZ2VXbqdRq5KLQEYQqU+YWHWtra8LDw0scNXXMTlGhU/k1SWq6rYuXcmbzdqljPJXNP/sylWdvK7Om/QFI9z+ikUxSSb93EpUyr/j9Vcajnd/rG1X9wnSPBrcbWZhX+bXrArUa0aIjCFWozIXOunXraNmyZfHX48aN44MPPtBIKE3TUcgpFC06NYJtBQod40ZdyIm9hzIzUVOxJKEuyCE74ibGjbpU+lrx/3w/6xtX/RTzrNQ0AIzMq2Y6fF2jVquRizpHEKpMmQudMWPGsHnzZpo2bcprr73G7NmzGThwoCazaYxKrUZWx8dJWjaw4/2je2jZt5fUUZ7qUYtOfBkLHZlCDyPHtmSFXtFkLMlkhV7BwMYDhZFVpa6TmP1PoaOBFp387KL9mvSNjKr82nWBSg2IFh1BqDJlHqMTGhrKuHHj+Pvvv4mIiGDgwIHk5mrZ7s9lpFSq0KnjM0L0jYywamCPXMs3ZDXVK8qXnle2/dT067shU+iSE+2ryViSyYn2A8DAtnGlirmsgqLlFQx1q/7fQUFeHmc2bSfCL6DKr10XiBYdQahaz/wtd+fOHdTq/x/PYmVlhUKh4OrVq0DR+jo1jbJQhY6ibhc6uvpFn+QLcsveJSQFHbkMZTn2VtO3aQxQbds9GBlAj7ZFH8Av3IKsHM3eLy/+PgAGNpUrdJSqon/Tuhoo+FWFhez/bmWVX7euEGMIBaFqPbPQGTZsWHXkAGDQoEH8+OOPKBQK1q9fz1dffaWR+ygLVejo1O1CR/FPS45Ky3ee15HLi38pl4WehQNqtYr8lEgNpiryzhRYMg0s/tkIPD0Tvv4NPv9Vc/cszEmlMCcNXQuHSl3n0fdUR0NNB4ZmphTk5qHMFyskl5e+roK8Au3+dykINckzC53w8PDqyIFcLuenn35iwIABREZG4u3tzb59+wgIqPrmb2Wh6LpSP9oJW8vHAhS16JS90NExqUdhVjKoNbvy9ewX4et5JR8zM4HPXofcPPhuq+buXZCZiI6JdaWu8ehvXVOze/537jAnf93CkZVrNXL92kxfV0FWOTevFQShdM/8bX/69GneeOMNnJycSjyuq6tLnz592LRpEy+//HKlg3Tq1ImgoCBCQ0MpKCjgjz/+YOTIkZW+7pMUiBYdslLTuH3sFJlJ5ds/qrrlFarQL0c3o8LQEmVOmgYTgVwO7z7lf/m3J4GuBoc+FWYnozCq3Bo1RroKADLzNdNyIFcoUIltVipEX1eHvALxvROEqvLMH8eDBw/mlVde4ffff6dRo0akpqZiYGCAQqHg2LFjLF++vEpWSHZwcCAiIqL468jISDp37lzp6z5JXr4SAz3tHoSraYkPI9jy1vtSx3im9LwCdBVyDHXk5CifPVZHJleASrOfhhs7g9NTNny3s4YWbnBLQ2sxqguVyHUNKnUNk3/+/8/SQKHzqFu0UHS/VIierkIUOoJQhZ752z4vL481a9awZs0adHR0sLa2Jicnh7S0qv3U/KQm9H8Pgn5k+vTpzJgxAyhaxLAiMnMLMDEQm+bVBI9mW5np65KjLMPAaZn8if/fVKW8MtRRGh/jXcn1ER4VOpka+IVqYFK0tURuVtlXsxb+n76ugnxRJApClSnXT0ulUklsbGyVFzlQ1ILz7+4xR0dHoqOjH3vdunXr8PT0xNPTk8TEii0Il5mTj4kG1g+pSUysLPnf+SN0Hj1c6ihPlfZPVWFRxsJUlZ+FXFez67eERoHPvdKf9wuGe2Gau79czxBVfuV2nG9gWtQiFJtZ9UtE6D8qdDJFoVMRRga6ZJdxOQVBEJ5NawaqeHt74+HhgYuLC7q6uowbN459+/Zp5F6Z2fmYGNbtFp38nFyMLcy1fvXamIyiphFHs7J11RTmpKEw1PzWA4tXwpM+dCuVRc9pksLQgsKc1Epd49H3MzK96gudvMwsDq34magA7d5HTRsp5DJMDPVIy9TuZR8EoSbRmkKnsLCQN954g6NHjxIQEMDOnTvx9/fXyL0yc/IxNaz6pe9rkvycHPJzcjGxqtwKu5oWlpqNSq2mkUXZdtpWZiagY2SBTNdQo7lOXIUBc+DMjf9/7PxNGDwPDl7Q6K3RNbVBmZVUqWs4mhmSma8kVQOze7JS0zi5bjOxQSFVfu3aztykqABNy6qZi7EKgjYq84jc48eP89Zbb3Hnzh2NhTl8+DCHDx/W2PUfyczJx9hAD5msaAO9uiojKQlTa+0udPIKVURn5OJqWbbuqLykMAD0rRqSG/eU/qUqcPYG9JkJxoZFM7EyqqGnRtfMHrmeEXmJoZW6jruVMWGplev+Ko2xhTm6+vqkxSdofLxUbWP+z5YnaeXY200QhKcrc4vOokWLWL58ORs2bMDO7ilTTmqA9Ow85HJZnW/VyUhMxqyCA7qrU2hqNm6WZWvRyUssakXQt/HQZKQSsnKqp8gB0LdxByA/seKtJQoZtKhvim98elXFKqHbuBf48MRe5AqFRq5fm4lCRxCqXpkLnZs3b9KvXz8OHDjAkSNH+OijjzAwqNwUV6kkphWt01/PvG5vOnj7+CkCL1+VOsYz+cal42ZljLn+sxsg8xNDKcxJw8i5fTUkq35GTh1QFxaQE1vxhTTdrUww0FHgG6eZQqeeowOpsXEUavmq29pIFDqCUPXKPUYnMDCQNWvWMHfuXB48eMCkSZM0kUujEtOKmuytzTU7jkPbndvyB6c3aHAJ3yriHZ0KQMcGFmV4tZrsiJsYN/TUZCTJGLt4khPli7qg4ptqtbUrGqytqRYda2dHEsM1vwVHbWT5z2y4VA3MhhOEuqrMhc758+eJiopi+fLlODg4MHXqVHr37k2nTp345ZdfNJmxyj1q0alvXrbukNpMz9AAHT3t7sLzT8ggK19ZxkIHMoLOo2fpVK3dV9VBx9QGA/vmZIZcqtR1ejnXIzQlixgNzewRhU7F2Voak68sFIWOIFShMg9GnjVrFn5+fo89Pm/ePI3NjtKUpPSiQsfaom636Dg0a8yCPzayacG7+J0+L3WcUhWq1fjEptHVsWwDpzMCT6EevASzZgNJqKZdzKuDWbNByGRy0vyPVPgaJnoKOjQwZ+sdzRQipvWsMK1nRWxw5QZL11W2lsbEp2TV6UkSglDVytyi86Qi55GhQ4dWSZjqoixUkZyeQ/06PkYnISwCuVyOvYeb1FGe6WRIAo5mhrSyefa6P4XZyWSGXMKizSiQ156tPizajCIn+i4FKRHPfnEpujtZoSOXc+5h5aanlyY/J5dt7y4l4NxFjVy/trO1NCYuRSy0KAhVqUrW0QkNrXmf3hLSsqlvUbcLnfycHBLDI2nQRPu7eE6HJZKrLGSIu02ZXp9y/Q90TW0wa9pfw8mqh7FrNwxs3En2/r1S1xne2I7YzFzuJmhmfE5edjY+B4+JrqsKsrUyIS5ZFDqCUJW0ZsHA6hadmIGjtXavClwdIv3v4dSimdQxnimroJBzD5MY4FofxRP2RfuvzOCL5CWGYN1jRqX3hdIG1t1foyAjvlLdVg6mBnRxtOLve7GoNNQ10rRHF2waNdTMxWs5XYWcemaGokVHEKpYzf8NUEER8enY1TNGR1FnvwUAhN2+i5WDPWb1tX89nYMP4rA01KO/a1myqok/+xMG9d0wb1mzulb/y8TDC2PnDiReWAuqik/ZHtXUnkKVmr2BMVWYrqRxn31Iv+kva+z6tZnNP2tFiUJHEKpW7RnAUE4RCRko5HLs65kQoaFptjVBwPlLqFWFFOTlSx3lmS5FJBOSksWUNs4cDU545usz7p0gJ8oX274LyHhwBlVuRjWkrFoyHX1s+71FXmIoKbf+qvB1jHUVPN/UnnPhSSRka+bvup6TI6b1rAi76auR69d2DW2Lpv2Hx1X9psmCUJfV2eaMR8WNcxkGt9ZmiQ8juLB9Nznp2l/sqYHf7kTQpJ4JXRwty3RO9OFPURhZYNvvLc2G0xAbrzno12tI7NFllWrNGdvCAQsDXTbcfFiF6Ury6NwBgJAbNzV2j9qs4T/rG4XX4Q9egqAJdbbQifxnMKZj/bpd6ACYWFnSqp+X1DHK5HBQPHGZebzS1rlMr8+LCyTp8mYs2z6PWbNBGk5XVWSYuPWgwbBPsOo8meQbO8kKq/gK1ka6Cia2cuR8eBIBiZlVmLOkxl07kRobR1xImMbuUZs1tDUjITWLLA1stCoIdVmdLXSycgtIzsip8y06AO2HDmLqD19iYWcrdZRnUqrUbLj1kPb2FvRrVLZxRfHnVpMdcRP7oUvRs3bVcMLKkuH04nKcx/2ERZtRyGRydM3sgGcPwC7NlNZOWBjost5Hc605AK4d2hJ4Ufu3FNFWDW0tCBPdVoJQ5epsoQMQGpOKa4OydYHUZvcuXAaguVd3iZOUzd/3YniQlMn8zm7ol2UwuUpJ5F+LUBfk0HDcanRMyzZFXQombt0xbdynxGOmHr0wcavY342TmSFT2jhx6EEcfgmaHaP05fCXOPbzBo3eo7aSyYpadB7GikJHEKpanS50HkSm4GpvgUJe8U/LtUF86EPiQx/Ssm8vqaOUSaEavr0cRANTAya3dirTOcqMeML/mIPcwBTn8T+jY6Kds8wMHds+8XEDu6YVut6i7u7kF6r48WpwJVKVTW5GJqmxcRq/T21kY2GMob4u4RraaFUQ6rI6Xugko6erKJ7tUJfdPXUWd88OGJqZSh2lTG7EpHEsOJ5X2jnjYVW2Pcty4+4RsXMeeub2uEzeiK65g4ZTlo/CuB7mLYY88bnc2Hvlvt5QD1u6Olqx5nooSTmaG/chk8mY+sOXNPfqobF71HaN/9neJCgqWeIkglD71O1C558fKh5l3EOpNrtz/AwKXR3cOraXOkqZfX0piPS8Aj7r26xsXVhAdvgNHm6bgcLQnEbTtmKkJbucG9i3wPWV39Exrkd25K0Sz2XcP01mcPm2VHA2N2Rxdw9uRKeyyz+6CpM+rlGHtrTq54WBSd1eabwymjjXQ1moIig6ReooglDr1OlCJyohg5w8pSh0gAi/AD4f8gJ3T52VOkqZpeYWsPRMIG6WxszvXPZBxjnRvoRumkxhdgoNJ/yCdffpEu6JJcOq0yRcpmxErVISunkKYZunEv7HHOLPrCT8jzlE7HqTosn1ZaMrl/FF32bkF6r48EyAxlZBfqTj8CHkZmVx99Q5zd6oFmvmXI/g6BQKlCqpowhCrVOnCx2VWk1wdHJxs3Fdlxyp2U/+mnA1KoVtvpGMbeHAc2XcBwsgP/khoZsmke5/FJveb9Bo2lYM7JtrMOnj9Ou703Dyr9gNeIes0CuEbphAXvx9QE1m8AUSL64nM/gC5SlyABZ396CptSmfngskPkuzC0HqGujTZmBf7hw/TX5OrkbvVVvJZNDYqR73wjWz0aog1HV1utABuBeeRGNHqzq/FQSATC5n0lef0H/mNKmjlMvKayFcj07lw15N6GBf9vFWqvxsova+R8Tuheia2OD6yu84jFyGrmXZBjhXlI6ZHfZDPsT1tZ0Y1Pcgat8HROycR2FOaqWvPbWNE6Oa2rPe5yFnNbRD+b+1GzwAAxNjvPce0vi9aitHazNMDPUIFIWOIGhEnf/t7hsSj76ejmjVAdQqFfomxnQb+zxyHYXUccpMqVLzznE/ItNz+GZACxqaG5br/IzAkwStGU7ChbWYNumD++v7cByzHGOXTlW4IagMQ8d2OIz8Ao85B7FoM5Lk678TtGYYab77q+QOA13r80YnVw4HxfHzjbAqueazpCclcePAEUKui9WQK6qJcz0AAsITJU4iCLVTnd3r6hHfkKI9k1q52uD/UPygufTHHl5b/R0t+3px59gpqeOUWUa+kvlH77JxRDt+eq41rx+8Q0R6TpnPV+VnkXD2J1Ju7MCyw0tYtn8RsyZ9KciIJ/3eCbJCr5AdfgNVXtlXFpbpGmLk1A6TRl0wbdofPQsHCvOySL62nSTvbSjTYyvyVp+oXyNr/tenKTdiUvnf2cAqu+6z3Dt/mXvnL1fb/Wqj1q42ZObki6nlgqAhdb7QScnMJTw+nVauNuw47S91HMndu3iFxPBIek8ZX6MKHYDojFzeOHyH1c+1Zt3wNsw+dIeQlOxyXUOZmUjC2Z9IvLgeE/demLd8Dst2L1DPcwJqVSH5KZHkJQZTkBKJMjuFwtw0UKtBJkOhb4qOsRW65vbo13dHz8oZmVwHlTKf7IfeJJxdTcb9U6jyy5fpWQa52fC/3k25E5/OwqN3KdD06ON/tB7QhwdXr5OTXvM2S9Umbd1tuRMSj0pdPX9vglDX1PlCB4q6r3q2ckImK/qdVZepVSrObN7OmA8X4dqhLSE3bkkdqVweJGcx48BtVj/XmrXD2jLn0B0Ck8q/v5NamUfGveNk3DuOTKGHoUNrjBp2wKC+B/r13TBp1BW53uNdZKqCXAoy4shLCCb93gmyw33IjriJWqmZgbovNLMvmkYek8rCo3fJqaZZO/VdnJn87WecWLuJoz+tq5Z71kbW5oY42Zix/9IDqaMIQq0lCh2KCp2hXdxpaGtOmFiCHe+9hzCzrkdCWLjUUSokNDWb6Qdusea5Nvw6oi3/OxvIsX+6KCtCXZhPdvh1ssOvl3hcpmuIwuD/F1gszM1AXVD27rLKUMhkvNPNnTHNG3A+PIl3T/iTV1h9U5P7TJtEYX4BF3/fXW33rI3auhXtL3crWKwoLQiaUucHIwPcCir6IdOhsZ3ESbSDMi+Po6vXk5FUc1dpjUzPZepeHwISM/miX3Pmd3ZFUcU7fagLclBmxBcf1VXkWBnqsvq51oxp3oDNt8N569jdai1yrJ0d6ThiCJd3/01msljgrjLaetiRkZ1HcJT4PgqCpohCB4hLySI8Lg3Ppg2kjqJV3DzbM/ztuVLHqLCknAJeP3ibnX5RTG7txOqhbXAwNZA6VqX0cbFmxwsdaWljxoenA1h5LVTjCwL+16DZr6HML+Dk+s3Ve+NaqK27LbeDxfgcQdAkUej84+q9aNq62aKnW3OmVWuaQ7PG9H55Am6eNWdbiP9SqtR8fSmIj04H0KSeCTvGdGRiK0dq2j6uZvo6fNK7Cd8MaEFsZh6T/rrB4aD4as8hk8vRNzLiwvZdZCaJVojKaGBtgoO1KTcfVN3sO0EQHicKnX94B0Sjr6dT3GcuwKUdf5EcHcPId+Yjk9fs/1UOBcXz4i5vrkWl8GYXN34d3o7WNmZSx3omhUzGSy0c2DO2E4PcbFl7I4ype28Smlq1M7fKSq1SsWHeIo6sWivJ/WuTrs0dAbjsHyVxEkGo3Wr2b68qdCs4jtx8JZ2aie6rR5R5eRz8/iccmjXGc+RQqeNUWkJ2PguP+fH+KX/sTQ3YMLId3w9siXsZdz+vTnJZUTfVzjEdeaebO4GJmUz66wZrfR5SKFE3h0ubVlg5Fv37UBUWSpKhNunWwoGw2FRiKjArUBCEstOKQufrr78mICCA27dvs2fPHszNy76Mf1UpUKq4FRRHJzFOp4RbR08SevMOz82fhaGZ9reAlMXR4ARG7bjKqmshtLUzY/voDnzdvzkd7C2kjoa+Qs7opvbsetGTbwa0QA0sOOrLnMN3CErOkiyXnqEhE7/+hMnffCpZhtrE2ECX1m62XPKLlDqKINR6WlHoHD9+nJYtW9KmTRvu37/Pe++9J0mOqwFRONmY4VQDujSq05+ffcPJdZvJy5LuF21Vy1Wq2HQ7gpF/XGPTrXDa21vwy7A27BrTkZdaOGBpoFuteVrUN+Wtrm4cGN+ZJT0bk5mv5N2T/ry025sL4dLPfhv8xnSsGtiz98sfpI5SK3g2tUdHIeeyn+i2EgRN04p1dI4fP1785ytXrjBmzBhJclzwjWD+C53wauPM1uN3JcmgjWLuBxFzP0jqGBqRka9k9fUwfr0ZTn/X+rzYvAHvdHPnra5u3IlL5+zDRM49TOJhWtVOHTfQkdPG1hzPBhb0d62Po5kheUoVF8KT2OEXhY8Wrefk1KIZPSeO5eIffxJ221fqOLVC1+aOpGXm4h8mtp0RBE3TikLn31555RV27NhR6vPTp09nxowZAFhbW1fpvRPTcvANiad3m4ai0HmCZr26M2DmNH5+bS75OdWzZkx1yStUcfBBHAcfxOFhZUwfF2u8GtZjfmc35nd2IzknH/+EDPwTMghIzCQ6I4fYzDyyCp49VsVMX4dGFka4WBjRyMKIFjZmtKxviq5CjlKl4kZ0Gut9HnI6LLFM16tOOvr6jP/iI9ITkzj04xqp49QKOgo5XZo7cNk/UkwrF4RqUG2FzvHjx7Gze3xBvvfff599+/YBsGTJEpRKJdu2bSv1OuvWrWPduqIl5729vas859nb4bzxfEcc65sSmSD28Pm33MxMnFo2Y/hbb/DnZ99IHUdjHiRn8SA5i7U+D7Ez0ae7kxUt6pvSor4Z3ZyskMv+f256Zr6ShKw88lVqlIWq4n2mTPV1MNfXwVRPF32d/+8hzlUWEpScxTbfSK7HpHI7Nq3atm2oCJkMgq7d4O6ps+Rm1p6uSyl5NrXHzFif0zcfSh1FEOqEait0BgwY8NTnp0yZwrBhw+jXr181JXqyR4WOVxtntp3wkzSLtgn1uc3Zzb/TZ9pE/M9dIuDcRakjaVxsZh5/BsTwZ0AMAEa6CtwsjbEz0cfWWB87E33qGemhK5ejK5eho5AjAx6mZpOWpyQjT0lSTj5hqdmEpmYTk5FLTfoMX5Cbx57Pv5U6Rq3St50LaVl5XA+MkTqKINQZaqmPQYMGqf38/NTW1tblOs/b21sjeVbOG6Re+9Zzkn9ftPHQ0dNTL9y5Wf2/80fUFna2kucRh2YOc9v66rlb16rtG7tJnqU2Hfq6CvWhL19SL3yxs8buoamfi+IQR009tGLW1apVqzA1NeX48ePcvHmTNWukHQtw9vZDPBytcBazrx6jzM9ny9vvo9DRod1zT2+lE2omuY6CSV/9Dzt3V5T5BVLHqVW6tXDEUF+Xkz5hUkcRhDpDKwYje3h4SB2hhJM3wpg1vD1DOrvxy/6bUsfROonhkXw7ehIpMWLp+tpo1OI3ce3Qlq2LPqqxO9hrqz7tGpKQmo1vSPVv3yEIdZVWtOhom5TMXK74RzGgoyvymrYpUjV5VOTYubvSfuhAidMIVaXr2OfpPu4FTm34jZuHjz/7BKHMzI316dLcgTO3HorZVoJQjUShU4rD14KpZ2YoVkp+hv4zpjLusw9p3NVT6ihCFWjVtxf+Zy9y6MefpY5S6wz0dEVXR8Ghq7VzTSpB0Fai0CnFFf8oUjJyGdLJTeooWm3XJ18SFxzKy8uX4dCssdRxhEr6de4iti76CLVKe6e811TDurhzNzSeMC1aDFIQ6gJR6JSiUKXm+PUQurZwwMxYX+o4WisvK5t1ry8kOy2d11Z/Tz1HB6kjCeVk3dCJqT98iaGZGYUFBeRlS7Mzem3W2tUGZ1tzDlwWrTmCUN1EofMUh68Fo6uj4LnOolXnadITElk3600UCgX9pr8sdRyhHOo5OvD6r6twadsKY0tzqePUWkO7upOZk8/Z22KRQEGoblox60pbhcWmcfNBLCO7N2bn6QAxgPAp4kMfsnLKTFKixUysmsKygR2zfl2Jrp4eq199g8SHEVJHqpVMDPXwau3M4WvB5OZr1xYfglAXiBadZ9hzPhA7KxO6thBdMs+SEBaOMj8fA1MTXl6+DCtHMZBbW1k3dGLOpjUYmBjzy4z5xD4IljpSrTW0izv6ejocuPxA6iiCUCeJQucZLvlFEpucyeheTaWOUmNY2tvh7tmeNzb9jK2ri9RxhCcozC8gMzmF1dPmEHXvvtRxai2FXMbonk3weRBLcHSq1HEEoU4Shc4zqFRq9l68T3sPOxrZW0gdp0aIuR/ET9NmI1PImbP5Zxq1ay11JOEfjs2bIJPLSYmJ5YdxrxBzXwyO1SSvNg2xsTRm15kAqaMIQp0lCp0yOHgliNx8JaN7NpE6So0R+yCYVZNnkpWSyqz1K2nas6vUkeq8zqOHM2/berymjJc6Sp3xYu9mhMenczUgSuooglBniUKnDDKy8zl2PYSBHV2xMjOUOk6NkRQZxYpJM/A9eZboQNFyIBW5joIRi+Yz9pMlPLhyncu7/pI6Up3QyrU+TZ3rsftsAGIegyBIRxQ6ZbTjlD8KhYwXvZpJHaVGyUlPZ+uij0iPT0Aml9N/5jQMTE2kjlVnmFhZMnPtCrwmj+Pc1h38Ovdt8rLEOjnVYWzv5qRl5XHMO0TqKIJQp4lCp4yikzI55fOQkd09MDPSkzpOjeTUshkDZ77Cwp2bcGrZXOo4dYK5TX0cmniw7b2P2fvVD6iUYnpzdXC1t6BHKyf2Xggkr0B8zwVBSqLQKYftJ+9iqK/LC2IGVoWE3/Fj9bTZyORy5m75hd5TJyKTiU1Tq5pcR0GLPj0BiLp3n88Gj8bnwFGJU9UtUwa1IjMnn93n7kkdRRDqPFHolENYbBrn7oTzfM8mGOnrSh2nRgq77cv3L77M3dPnGP7WG7z06ftSR6pV7Bu7MX/br7yy4mscmxcV5LkZmRKnqlsa2Vvg1aYhe84HkpGdL3UcQajzxMrI5bTt+F16tXbm+Z6N2XbCT+o4NVJOegZb3nofz1FDi1fj1dHXp7CgQGwmWUEKXV36vTqZfjOmkpOewaY33yPSX7QmSGHKwFZk5eaz+6yYUi4I2kAUOuV0PzKZy36RvNSnOXsvPiAzR3xiqyjvvw8W//m5eTNx69iePz/7mnBffwlT1UyzN/yES9tW3Dx0jL+WfU9WqtghWwqN7C3o3bYhW475itYcQdASouuqAtYfvIWxgR4T+rWQOkqtEXbLF9N6Vszduo6xnyzBrL611JG0Xj1HB+QKBQBnf/uDtTMXsHXxUlHkSOjlQaI1RxC0jSh0KiAkJpUTPqGM7tkEa3MjqePUCneOn+arEeM4u/l3OgwbxLsHdtJ2UD+pY2klY0sLRr37Jov3/YHnqKEA3Dl2isBLVyVOVre1cLHGq01Ddp4JEK05gqBFRKFTQRsO3UYml/HyoFZSR6k18rKzOfD9Kr4aOR6/M+eJCizaBNHMpj4GJsYSp5OeST1Lhi18g/eP7KHbS6O5+td+/M9ckDqW8I9ZIzqQmJbNztOi61UQtIkYo1NBcSlZ7Lt4n+d7NmHnmQAi4tOljlRrJEdGs23x0uKvn3/3Tdw7d+DSjr+4+Ptu0hMSJUxX/eo5OdJr0li6jBmJXKHA59AxTqzdREJYuNTRhH/0auNMy0b1+eaPy+Tmi3VzBEGbiEKnErYev8uQzm7MHNaODzaclTpOrXVy/WaQyej76mR6T53ArcMnOL9tZ62eVSSTyWjctRPdXnqe5r17olIqubH/CKc2bi2eqSZoBx2FnBnD2hESncKRa2IVZEHQNqLQqYS0rDy2Hr/LzOHt6dS0AdfuRUsdqVaK9A9k85vvUc/RgZ6TxtLp+WGkxMYS6X8PuY4CIzMzMpNTpI5ZJfQMDcnPyUGuo8PELz8G4OS6zVz8fTcZScnShhOeaGT3xjhYm7Lol5OoxKZWgqB1ZECN/Zfp7e2Np6enpBl0FHJ+XTQMGfDq1wcoKBTrwGiavpERch0FOekZNOvVnWk/fknQ1evcPHwc/7MXtWrWkUwmY/Dg9rRv74aPTzBHjvig/s8vQwtbG1r260W7IQMxtbbiiyFjALBv7E586EMKCwqkiC6UgaWJAZvfG869iCQW/XxK6jiAdvxcFARtIlp0KklZqGLlHm++mdWPF7ya8scpMRBR0/Ky/39TyriQUE5v2Eq75wYw7rMPURUWEnbLl01vvkdWSqp0ISkqcvb8tYSRI7sUP7Z37xVGP/8FarWaln17MXDWqzg0awxA9P0gruzeh0JXl8KCAmLuix3ftd3MEe0w0NNh5Z7rUkcRBKEUotCpAtcDYzh/J5wpA1tx4kYYiWlid+jqkhwZzeGVv3B45S84Nm9Ki949cG7VorjIGf7WXOo5OfDwti8Pff2J9LtHfk5OtWQbPLh9iSIHYOTILkx4dRjb1u9HrVaTn5PDge9X4XfmAvGhD6sll1A12rjZMMjTja3HfcVkBEHQYqLQqSKr995g0+LhvD6iPZ/+Jqb8SiHS/95jA5QL8vOwb+xGq35eAKgKCwk4d4kN8xYB0KxnN/Kys0mOiiE9IRFVYcVmzCh0dbGwtcHKsQEpMbEkPoygV/+OT3xt+w7ubFsPfqfP43f6fIXuJ0hLRyFnwZhOxCRlsvX4XanjCILwFKLQqSKxyVlsP+nHtCFtOH4jlCv+UVJHEoAjK9dyZOVajC3McWrVHOcWzcj51yaXL336Pqb1rIq/zkpN4/r+w+z7+kcAJn/zKWqVqnggm0wmI/DSVbz/Poi+kRFv7tyEiaUFhmamxdc4tuZXjq5ez5VLfrBw2GOZTuw7p5k3K1SbMV5NcbGz4L11p8krENPJBUGbiUKnCm0/6YdXG2cWvtiJaV8dICtXDCLVFlmpadw7f5l75y+XeHzVlJlYOTTAyrEBZtb1MLGyJPZ+cPHz9V2c0TMwKP5arVYT86Do+fycHCL8AshMTiErJZW0+ASSI6OJCwkDYN+ei+zde+WxMTpHjvho8J0KmmZnZcyUga05fydcfKARhBpCrS3HW2+9pVar1ep69eqV6fXe3t6SZ/7v0cSpnvrEdxPUC1/sLHkWcUh/yGQy9ZAhHdRLloxVDxnSQS2TySTPJI7K/H2i/n52f/WBZWPV9S2MJM/zpEMbfy6KQxxSHlrTouPo6MiAAQN4+LBmD8gMjEhi15l7jOvbnFM3w7gVFCd1JEFCarWaw4dvcPjwDamjCFVgZPfGtPOw45s/LpOQKiYdCEJNoDV7XS1fvpxFixY9tsZITbTpyG0iE9J5+6UuGOgppI4jCEIVaGBtwoxh7bniH8Whq8HPPkEQBK2gFYXO8OHDiYqK4s6dO8987fTp0/H29sbb2xtra+tqSFd+eQWFfPPHFRysTZk5vL3UcQRBqCS5TMa747tRqFLx3c4rUscRBKEcqq3r6vjx49jZ2T32+Pvvv8+SJUsYOHBgma6zbt061q1bBxStAKqt7oTEs/O0P2P7NOfavWgu+4lBi4JQU73QqymtXG1Ytu0iiWnVsw6TIAhVo9oKnQEDBjzx8ZYtW9KoUSNu374NFI3V8fHxoVOnTsTF1ezxLesP3qKdhx2LxnXl1W8OkpwufkAKQk3j4WjFa8PactE3gmPXQ6WOIwhCOUnedXX37l1sbW1p1KgRjRo1IjIykvbt29f4IgegoFDFp79dwEBPh3cndEUmkzqRIAjlYaivw0dTepCWmcvXO0SXlSDURJIXOrVdRHw6q/++gWeTBozp1UzqOIIglMObYzpjX8+ET3+7QHpWntRxBEGoAK0rdBo1akRSUpLUMarU/ssPuOAbwWvD2tLEqZ7UcQRBKIPBnVwZ0LERm4/44huSIHUcQRAqSOsKndrqmz8uk5yewyfTemFurC91HEEQnsLZxox5ozvhcz+GbSfEXlaCUJOJQqeapGfn89HGc1iaGPDhlB7I5WLAjiBoIyN9Xf73ihe5+Uo+33YJVS1Y20sQ6jJR6FSjB5HJLN99lQ6N7XllSBup4wiC8B8yGSyZ1B1Ha1M+2XxezJQUhFpAFDrV7Mi1EPZdus/E/i3p0cpJ6jiCIPzL1EGt6d7SkVV/X+d2cM2f+SkIgih0JLFqz3UCHiby7oSuuNiZSx1HEASgZysnpgxqzeGrQfx94b7UcQRBqCKi0JFAQaGKjzaeIydPybLpfbA0NZA6kiDUaS525rw3sRv+YQks33VN6jiCIFQhUehIJDEtm/fXn8Hc2IDPX+2Nvq7Y/FMQpFDPzJBl0/uQnafko43nKChUSR1JEIQqJAodCd2PTObzrRdo4lSP9yZ2EysnC0I1M9TXYdn0Ppgb67Nk3WmSxOBjQah1RKEjsYt3I1mz7wZebRoyfWg7qeMIQp2hkMv4+OVeuNpb8PGm89yPTJY6kiAIGlBtm3oKpdt99h6O1qaM79eC+NQsMRBSEKrBW2M706lZA7754zLX7kVLHUcQBA0RhY6WWPHXdeqZGzH/hU5k5hRw4obYJVkQNGXq4NYM6ezOpiN3OHQ1WOo4giBokOi60hIqlZr/bTmPz/0Y3h3flW4tHKWOJAi10otezXj5n2nkm4/ekTqOIAgaJgodLVKgVPHBhrPcj0xm6cs9aetuK3UkQahVRnTzYPaoDpy59ZBvd16VOo4gCNVAFDpaJidPyeK1p4hKzODzV3vT1Fnsdi4IVWFwJ1fefLEzF+9G8vnWi6hUYg8rQagLRKGjhTKy83nn55OkZObyzax+otgRhErq286Fd17qindgNJ9sPodSrJUjCHWGKHS0VFJ6Dgt/Ok5aVh7fzOpH84bWUkcShBqpVxtnlkzsxp2QeD789SwFSlHkCEJdIgodLRafms2CVcdJyShq2WnlWl/qSIJQowzydOWjKT3wf5jIkvWnySsolDqSIAjVTBQ6Wi4xLZs3fzpOQlo2X83oRxs3MUBZEMpiZPfGvDuhGzcfxLHol5Pk5CmljiQIggREoVMDJKXn8Oaq48QmZ/LljD54NrWXOpIgaLVxfZuzYEwnLvpGsGT9aXLzRUuOINRVotCpIVIyc3nzp+OEx6fzxWt9GNixkdSRBEErTRvShpnD23PSJ4ylm86JMTmCUMeJQqcGScvK481Vx7kdHMd7E7szvl8LqSMJgtZQyGW8NbYzUwa24uCVIL7YepFCMYVcEOo8UejUMNl5Bby79jQnboQyY1g75j7fEbnY9lyo44z0dfnitT4M6+rBlmO+fLvjCiq1KHIEQRB7XdVIykIVX2y7SFJ6Di/1aY6VmSHLtl8iX8woEeoga3MjvpzeBxc7c77547LYu0oQhBJEoVNDqdXw8z4fktJymD2qA7aWxny44SxJ6TlSRxOEauNqb8GXM/pgZKDLu+tOcz0wRupIgiBoGdF1VcPtOhvAB7+ewcXOnJ8XDhGrKAt1RveWjqyYNxCA+SuPiSJHEIQnEoVOLXDxbiRzfjxKgVLFj28MpH97F6kjCYLGyGRFM6s+e7U34XHpzP7hCMHRqVLHEgRBS4lCp5YIjUll1vLD+D9M5P3JPZg+rK0YpCzUOsYGRYOOpwxsxaErQcxfdYzENNFdKwhC6cQYnVokPSuPt9ecYO5oTyb0a0lTZ2s+33oBDwcrPByteBCZzLV70YjJKEJN5GJnzqeveGFracz3u66y/9IDqSMJglADiEKnlilUqflh9zUCHiayYEwnfv/gefR0FcXPX/SN4MONZ0WxI9Qogzxdmf+CJ9l5ShauPsHd0ASpIwmCUENoTdfVG2+8wb1797h79y5fffWV1HFqvKPeIazZe6NEkQPQvZUTnZo2kCiVIJSPkb4u70/qzrsTuhEYkczM7w6JIkcQhHLRihad3r17M3LkSFq3bk1+fj7164tduquCmbH+Ex93d7DiakB0NacRhPJp6lyPDyb3wM7SmF8P3WL7CT+xCKAgCOWmFYXO66+/zpdffkl+fj4ACQniE1tVeBCZ/MTHbS2NkMkQ3VeCVpLJ4KU+zXn1ubYkpmUzf9Ux/MISpY4lCEINpRVdV40bN6Znz55cuXKFM2fO0LFjx1JfO336dLy9vfH29sba2roaU9Y81+5Fc9E3osRjyRk5DO/WmK9m9sXa3FCiZILwZA2sTVg+ZwAzh7fngm8Er31zUBQ5giBUigyols/1x48fx87O7rHH33//fT7//HNOnTrF/Pnz8fT0ZMeOHbi6uj7zmt7e3nh6emoibq0hk0Gnpg1wd7AiKKpo1tWwLh68PrIDBcpCVu7x5oRPmNQxhTpOLpMxulcTXn2uLcpCFav+us5R7xCpY9VI4ueiIDxOLfVx+PBhtZeXV/HXQUFBamtr62ee5+3tLXn2mno41jdVr5o/SH16+ST1VzP6qm0tjSXPJI66eTjVN1OvnDdQfXr5JPXnr/VW1zMzlDxTTT7Ez0VxiKPkoRVdV3///Td9+/YFwMPDAz09PRITRXO1JkUmZDBvxTFW7PGmZaP6bFw8nBe9miGXi0UGheqhkMsY17c5695+Dmcbcz7fepH3158R+7UJglCltGIw8oYNG9iwYQO+vr7k5+fz8ssvSx2pTlCp1fx1PpALvhEsGNOJ2aM60K+DC9/uuEJQVIrU8YRarJ27LfNe8MTFzoLzd8L54U9vkkWBIwiCBlTbGB1NEH3RVatP24a8Mboj5sb6HLgcxMbDt0nLypM6llCLWJsb8vqIDvRt70J0UgYr91znin+U1LFqFfFzURBK0ooWHUE7nL71EO/AGKYObs2o7o3p264hm4/e4e8L9ylU1dh6WNACOgo5L/RqysuDWiGXydh45DZ/nPInv6BQ6miCINRyotARSsjMyWfVX9fZf+kBc0Z14I3nPRnerTE//X0d73sxUscTahiZDHq3acirz7XFob4pF+9G8tPf14lJypQ6miAIdYTouhKeqktzB+aM6oBjfTO8A6P59eBtAiOSpI4l1ADt3G2ZMbw9TZ3rERydwtr9N7l2T6zIrWni56IglCRadISnuuIfxY3AGEb1aMLE/i34eeEQzt8JZ+ORO4TGpEodT9BCbg0smD6sHZ2bORCXksWybRc5cSNMbN8gCIIkRKEjPFNBoYpdZwM4eCWIMV5NebF3M7q3dOLUzTA2HrlNdKLohhDA3cGSSQNa4tWmIelZeazZe4O/LgRSoFRJHU0QhDpMFDpCmWXnFbDlmC9/XQhkXN/mjO7ZlD5tG3LyZhh/nPIXLTx1VLOG1kwa0JJuLRzJzMnnt2O+7DwTQGZOvtTRBEEQRKEjlF9Gdj7rDtxi99l7jO/bgmFd3RnY0ZUr/lH8ftKPOyHxUkcUqkEbN1smDWhJxyb2pGXmsv7gLf6+EEhWboHU0QRBEIqJwchCpZka6TGqexNG92qChYkBfmEJ/H7Sj8t+UWJcRi2jq5DTt70Lo3s1pbGjFcnpOew47c++Sw/IzVdKHU9A/FwUhP8ShY5QZfR1FQzp5MbYPs2xr2dCTFIm+y8/4PDVIFIzxcKDNZmVmSEju3kwrJsHVqaGhMWmsudcIEevh4i1cLSM+LkoCCWJQkeocnK5jJ6tnBjRvTHtPezIVxZy7nY4ey/e525ogtTxhDKSyYq6p4Z2ccerjTMKuZzL/lHsOX8Pn/uxUscTSiF+LgpCSWKMjlDlVCo1Z2+Hc/Z2OA1tzRnRzYOBnq7079CI4OgUjlwL5qRPGCkZuVJHFZ7A1tKYQZ1cGeTpSoN6pmTm5LP34n3+uhAoZtgJglDjiBYdoVoY6Cno264Rw7t50NS5HoUqFTcCYzl+PYQLdyPIzRfdH1Iy0Pu/9u49tqmC7wP4t7d1F3aBDejYsOseHMx3MLlKuIiAgBMF/8AYMEKCAWOChoUEESNCSAwiCfAmxkdhkYsQlTCM4AZOYNyUy1w3OtaxDQod666Uld1v/b1/DPo6mA/wADtd9/0k37SnO2zfNmv347TnHC0mJUTjlRf+hZFDDAAAc3E50s9dwSlLCd+e6kH4ukjUGbfoULdoamlH2rlipJ0rxjMDQvDyaBNmjDHhk3cmobG5FacsJThuvo7sonL+Ue0mAXotJvxPNF5MfAYvDBsEvZ8Wjpu12HnkIo5cuIqKW/VKVyQiemzcokOKUamABFN/zBgdi5eefwbBgXo0NrfifIEDZ/Ju4Gx+KWobeCyWJyk40A8vxEdhSuIzGDdsEPx0GlTVNODURTsyc68jz1YF7ijXs/F1kagzbtEhxYgAlqtVsFytwv+mXsDzQwZiUkI0JiQMxpREI9rb3ci9UoE/80vxV2E5D0j4X1CpgLjocLwQPwhjhw1CvDEcGrUaVTX1+OWPQpzItePSNQ43ROS7uEWHvI5KBQwdHI5JCYMxcXg0YgxhAADn7UZkF5Uju7AcfxWWobKmQdmiXmpQRB+MiB2A0XGRGDM0EmF9/OF2Cy6X3MT5AgfOWR0osFdzuPFRfF0k6oxbdMjriAAF9psosN/E9rQcDAgLxKg4A0Y9G4nRcQa8PNoEACitqsWla1WwXq9G/vVqXHHcQru7d/31VqtVGBLVF8NNAzA8tj+GmwagX0gAAMBZ24jzVgfOFziQdbkMrnoey4iIeh8OOuT1KmsacPj8VRw+fxUAEGMIxeg4AxL/NRCjh0Zi5thYAEBzSxsKbzhhvV6NotJbsJXVwF7p8pmTSup1GpgiwzAkqi+GRPXDkKi+iI3siwB9x9PYcbMWWZfLYLFVwnK1CvZKF7faEFGvx0GHepxr5S5cK3dh/8nLAIABYYF4LiYC8c9E4LmYCLwxaSj8dBoAQLvbjdKqWlwrd8FWXoPr5S6UOetQ7qzz2qM19wsJQHREMKIighHVv+MyxhCKwQNCoFGrAQB1jS0oLr2FX88WIf9aNSy2KlS7+FYeEdG9OOhQj1dZ04DKHDsyc+wAAI1ahej+ITBFhsJkCEOMIQwxkWGYODzaMygAQFNLG8qddSh31qPCWQdnbRNq6ppQU98MV10zXPUdy7cbWuB+zLfE/P20CPLXIVCvQ0iQH8JDAxEeEoDw4AD0CwnouB4agEHhfRCg13n+XVu7G46bdSipvI0TuXYUl95CcakT5U7u+k1E9DA46JDPaXcLrle4cL3ChUzYPbf76TSIigiGoV8QDP36wND3zmW/IMQbIxAapP/H79na1o7m1rtpQ8ud6263QKVSQa3CnUsV1GoVNGoVAvRaBOg7hhu1WtV113Y3nLVNcNY2otxZD3NROUqra1FaVYsb1bWouFX/2EMWEVFvxkGHeo2W1nbYymr+cTd1jVqF0D7+CA3SIyxIj7A710OC9PDTaaDXaeDvp4WfTgN/XcelWqWCWwQiArdb7lwH3CJoaGpFQ3Prncs2z/LthmY4bzfi5u1GuOqb+TkaIqKniIMO0R3tboHzdiOctxuVrkJERE+I+sGrEBEREfVMHHSIiIjIZ3HQISIiIp/FQYeIiIh8FgcdIiIi8lkcdIiIiMhnecWgk5iYiD///BNms5ln3iUiIqInxisGnY0bN2LdunUYOXIk1qxZg40bNypdiYiIiHyAVww6IoKQkBAAQGhoKBwOh8KNiIiIyBeoACh+APphw4bhyJEjHecKUqsxYcIE2O32LtddsmQJli5dCgCIiIiAyWTqzqpERF6Nb/8TddZtg05GRgYMBsN9t3/yySeYPn06Tpw4gdTUVLz55ptYunQpZsyY8cDvySc0EVFnfF0k6swrtujU1NQgLCzMs+xyuRAaGvrAf8cnNBFRZ3xdJOrMKz6j43A4MGXKFADAtGnTUFRUpHAjIiIi8gVecfbyJUuWYOvWrdBqtWhqavJ8BudBjEYjLly48FDrRkREoLq6+nFqPhXs9WjY69Gw16PxhV5Go/EptyHqeaQ35MKFC4p3YC/2Yi/28qVeDNMT4hVvXRERERE9DRx0iIiIyGf1mkHn22+/VbpCl9jr0bDXo2GvR8NeRL7HK3YvJyIiInoaes0WHSIiIup9es2g481nSF+2bBkKCgqQl5eHL774Quk691mxYgVEBOHh4UpXAdBxElir1Yrc3FykpqY+1MEln5ZZs2ahoKAARUVF+OijjxTrca/o6GgcO3YM+fn5yMvLw4cffqh0JQ+1Wo3s7GwcPHhQ6SqdhIaGYt++fbBarcjPz8f48eOVrgQAWL58OfLy8mCxWLB3717o9XqlKxH1OIrv+tUdOXLkiLzyyisCQJKSkuT48eOKdwIgL730kmRkZIifn58AkP79+yve6e+Jjo6Ww4cPy7Vr1yQ8PFzxPgBkxowZotFoBIBs2LBBNmzYoEgPtVotxcXFYjKZRKfTSU5OjsTHxyv++AAQg8EgI0eOFADSp08fuXz5std0S05Olj179sjBgwcV7/L37NixQ959910BIDqdTkJDQxXvNGjQILl69ar4+/sLAPnxxx9l0aJFivdimJ6UXrNFx1vPkP7+++9jw4YNaGlpAQBUVVUp3KizzZs3Y+XKlRARpat4ZGRkoL29HQBw9uxZREdHK9Jj3LhxKC4uhs1mQ2trK3744QfMnTtXkS73Ki8vh9lsBgDU1dXBarUiKipK4VZAVFQUZs+eje3btytdpZPg4GC8+OKLSElJAQC0trbC5XIp3KqDVqtFQEAANBoNAgMDvea1i6in6DWDzvLly/Hll1/Cbrdj06ZN+Pjjj5WuBACIi4vD5MmTcfbsWWRmZmLMmDFKV/J4/fXXUVpaiosXLypd5R8tXrwY6enpivzsqKgolJSUeJZv3LjhFcPEvYxGI0aOHIlz584pXQVbtmzBypUr4Xa7la7SSWxsLKqqqvDdd98hOzsb27ZtQ2BgoNK14HA4sGnTJtjtdpSVlcHlciEjI0PpWkQ9ilecAuJJedAZ0pOTkz1nSE9JSXmoM6Q/7V5arRZ9+/bF+PHjMXbsWPz000+IjY3tll4P6rZ69WrMnDmz27r83X/q9csvvwAAVq9ejba2NuzZs6e76wEAVCrVfbd505YvAAgKCsL+/fuxfPly1NbWKtpl9uzZqKysRHZ2tufcdt5Cq9Vi1KhR+OCDD3D+/Hls2bIFq1atwpo1axTtFRYWhrlz58JkMqGmpgb79u3D22+/rdjvPFFPpfj7Z92RmpqaTssul0vxTgAkPT1dpkyZ4lkuLi6WiIgIxXslJCRIRUWF2Gw2sdls0traKtevX5eBAwcq3g2ALFy4UP744w8JCAhQrMP48ePl8OHDnuVVq1bJqlWrFH9s7kar1crhw4clOTlZ8S4A5PPPP5eSkhKx2WxSVlYm9fX1snv3bsV7AZCBAweKzWbzLE+aNEkOHTqkeK958+bJ9u3bPcvvvPOOfPXVV4r3YpgeFsULdEvy8/M9A8W0adMkKytL8U4A5L333pN169YJAHn22WfFbrcr3qmr2Gw2r/kw8qxZs+TSpUuKD4QajUauXLkiMTExng8jP/fcc4o/Pnezc+dO2bx5s+I9usqUKVO87sPIJ0+elLi4OAEgn332mWzcuFHxTuPGjZO8vDzPQL9jxw5ZtmyZ4r0YpodF8QLdkokTJ0pWVpbk5OTI2bNnZdSoUYp3Ajr27ti9e7dYLBb566+/ZOrUqYp36ireNOgUFRWJ3W4Xs9ksZrNZvv76a8W6JCUlyeXLl6W4uFhWr16t+GNzNxMnThQRkdzcXM/jlJSUpHivu/HGQScxMVEuXLggubm5cuDAAQkLC1O8EwBZu3atWK1WsVgssmvXLs8emgzDPFx4ZGQiIiLyWb1mrysiIiLqfTjoEBERkc/ioENEREQ+i4MOERER+SwOOkREXiAlJQUVFRWwWCwPXHfw4ME4duwYsrOzkZubi6SkpG5oSNQzca8rIiIvMHnyZNTV1WHXrl0YPnz4f1z3m2++gdlsxr///W/Ex8cjLS0NJpOpm5oS9SzcokNE5AVOnToFp9PZ6bbY2Fikp6cjKysLJ0+exNChQwHAa09STOStFD+YD8P0lPj7+0tmZqao1WoxGo1isVj+q++j0+nkxIkTotFoFL9PjPfk3t+p33//XYYMGSJAx1GSjx49KgDEYDDIxYsXpaSkRJxOp9ccAJVhvDHcokP0CBYvXozU1NTHPvt2a2srjh49irfeeusJNSNfExQUhAkTJmDfvn0wm8345ptvEBkZCQCYP38+duzYgcGDB+PVV1/F7t27uzzJLBF1UHzaYhhvyLFjx+Tll18WALJ+/XrZunXrfeucOXNGjEajAJ3/920ymSQ7O1vGjBkjRqNRrFarbNu2TSwWi3z//fcyffp0OX36tBQWFsrYsWMFgIwYMUJ+/fVXxe834z35++9UcHCwOByOLtfLy8uT6Ohoz/KVK1ekf//+ivdnGC+N4gUYxisyefJkOX78uCxYsEAOHTokarW609d1Op2UlZV5lu/+UYqLi5Ps7GxJTEz03N7a2ioJCQmiUqkkKytLUlJSBIDMmTNHDhw4IABErVZLZWWl4veb8Z7c+9bVmTNnZN68eZ7lESNGCABJS0uTRYsWCQAZNmyYlJaWKt6dYbw4ihdgGK9JZmamZGVlSZ8+fe77WmRkpFitVs+y0WiU8vJysVqtnc5abjQapbCw0LO8c+dOWbBggQAdW37MZrPnazdu3OjyZzG9L3v37hWHwyEtLS1SUlIiixcvlpiYGElPT5ecnBy5dOmSfPrppwJA4uPj5fTp05KTkyNms1lmzJiheH+G8dZoQUQAgISEBERGRqK6uhp1dXX3fb2xsRH+/v6dbnO5XCgpKcHEiRORn5/vub25udlz3e12e5bdbje02v9/2un1ejQ1NT3pu0I90IIFC7q8vatj5FitVkyaNOlpVyLyCfwwMhEAg8GAPXv2YO7cuaivr8fMmTPvW6empgYajQZ6vd5zW0tLC9544w0sXLgQ8+fPf6Sf2a9fP1RVVaGtre2x+xMRUdc46FCvFxAQgNTUVKxYsQIFBQVYv3491q5d2+W6v/32233/k25oaMBrr72G5ORkzJkz56F/7tSpU5GWlvY41YmI6CEo/v4Zw/SUPP/887Jr164n8r32798vcXFxit8nhmEYXw636BA9gpycHBw/fhxq9eM9dXQ6HX7++WcUFhY+oWZERNQVnuuKiIiIfBa36BAREZHP4qBDREREPouDDhEREfksDjpERETkszjoEBERkc/ioENEREQ+i4MOERER+az/A7AN3y5xZwmhAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 1008x432 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"#グラフを書く\n",
"frame = StaticOrbitPlotter(plane=Planes.EARTH_ECLIPTIC)\n",
"#惑星を描写\n",
"frame.plot_body_orbit(Earth, EPOCH, label=\"Earth\")\n",
"frame.plot_body_orbit(Mars, EPOCH,label=\"Mars\")\n",
"frame.plot_body_orbit(Jupiter, EPOCH, label=\"Jupiter\")\n",
"#didymosを描写する\n",
"didymos_orbit_plot_t,didymos_orbit_plot_p=frame.plot(didymos_orbit)\n",
"didymos_orbit_plot_p.set_marker(\"\")\n",
"frame.plot_ephem(didymos, EPOCH, label=\"didymos\")"
]
}
],
"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.8.0"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment