{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 16 - Mass-Benefit Analysis - Neptune" ] }, { "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(\"NEPTUNE\")" ] }, { "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 = 400000.0 \n", "apoa_alt_km_ab = 5248815.88" ] }, { "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/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACJLUlEQVR4nO2dZ3hURReA3wmkQULvLaEJSlUQAogBpQlIE0QMgnSxUkT0A6UoCioiRbpSpQhSRRCQ3lSaAgqIVCmC0msScr4fs5vsJpuyyW42Zd7nmWf3zp1y9hL27MycokQEg8FgMBhSEy9PC2AwGAyGzIdRPgaDwWBIdYzyMRgMBkOqY5SPwWAwGFIdo3wMBoPBkOpk9bQAqUm+fPkkODg42f0vXbpE/vz5XSdQJsI8u5Rhnl/yMc8u+aT02e3Zs+dfEXE4QKZSPsHBwezevTvZ/atXr56i/pkZ8+xShnl+ycc8u+ST0menlDoV3z2z7WYwGAyGVMcoH4PBYDCkOkb5OEHPnj09LUK6xTy7lGGeX/Ixzy75uPPZKU+H11FKFQMGAtWBKoA/UFJETiahr5elby+gEHAEGC4i3zpqX716dTF7vwaDwZA6KKX2iEh1R/fSwsqnDPAscAXY6mTf94GhwATgKWAXsEgp1dSVAhoMBoPBtaQFa7ctIlIQQCnVHWiUlE5KqQLAm8BIEfnUUr1RKVUGGAl87w5hDQaDwZByPL7yEZGoZHZtDPgAc2PVzwUqKaVKpkgwg8FgMLgNjyufFFABuAcci1V/yPL6UOqKYzAYDIakkha23ZJLHuCqxLWYuGxz345Lly5RvXrM2VfPnj2NJYzBYDC4kKlTpzJ16lTrZb742qVn5eM0+fPnd95b9zwwDVgPKECAhkB3oLCLBTQYDIZ0ju2PeqXUv/G1S8/bbleAXEopFaveuuK5TEoQ4FOgI1AR+BHYbHmtYKn/1NLOYDAYDE6RnpXPIcAXKB2r3nrW83uKRh8NnALWAW0Ab0u9t+V6neX+6BTNYjAYDJmS9Kx81gARQFis+o7AQRE5keyRzwOrgbHE/4S8LPdXW9obDAaDIcmkiTMfpVRby9tqltenlFKXgEsistnSJhKYJSLdAETkolLqM+AdpdQNYC/QHngCaJEigaYDr5C4avYCXga+BAanaEaDIW3Qpg1cdmLHOk8eWLLEffIYMixpQvkAi2JdT7S8bgbqWd5nsRRbBgE3gTeICa/zrIh8lyJp1gFvx6r7BW3U3QzIYVPfAhiHUT6GjMHly7BpU9Lb16vnLkmimTlzJl26dIm+DggIoFSpUvTo0YOXXnqJrFnTytdYDEOHDmXYsGE4E76snuVZbnLm+adj0sS/mojENhpIUhsRuQ98YCmuQxFzxmPlDtp9tRdQH2gLPA3kcunMBoMhHhYtWkSxYsW4fv06ixYt4rXXXuPixYsMHz7c06K5hIkTJybeKAORJpRPmkPQp0m2CuhxS7kKrAQWo7fm6gAX0LZ1cTyLDAaDq6hatSplypQBoFGjRhw7doyxY8c6VD4RERFkzZqVuMawaZeHHspcfvHp2eDAfTREKxhH5AJeAJYDfwPl0So8GB3wZxpwye0SGgyZnkcffZTr16/z888/o5Ri4sSJvPXWWxQpUgRfX1+uXr2KiDBmzBjKlSuHj48PhQsX5tVXX+X69et2YymlGDRoECNGjKBYsWL4+/vz+OOPs3//frt2wcHBvPjii3FkUUoxdOjQBOUdO3YsDz74IP7+/uTOnZvq1auzdOnS6Pv16tWL3nq7cOECWbNmZdy4cXHG+fjjj/H29ubSpZgvmiVLlhASEkK2bNnIlSsX7dq14/Tp0wk/QA9jlI8jugNfAIlFnQsAfkP7/pyz9FuPjtP9JDAJvSoyGAwu58SJE2TJkoWAgAAARowYwdGjR5k6dSpLly7Fz8+PQYMG0a9fPxo2bMjKlSt56623mDlzJs2aNSMqyv4/+OzZs/n++++ZMGECM2fO5J9//uHJJ5/ksjMGGPHw9ddf079/fzp06MD333/P119/Tdu2beMdu1ChQjRo0IC5c2OHroQ5c+bQpEkT8ufPD8DkyZN55plneOihh1i8eDFTpkzh4MGDhIaGcuPGjRTL7jZEJNOUatWqSZL5REReFZH78dy/LyKvWNrF5paILBGR50Ukp4jUFZGxInIm6dNHc05EhlnGeNzyOtxSbzAkFfB8SSYzZswQQA4fPiwRERFy+fJlmTx5snh5eUnLli3lxIkTAsjDDz8sUVFR0f3+++8/8fHxkc6dO9uNN2fOHAFk+fLlNo8HyZs3r9y8eTO67sSJE5I1a1YZPHhwdF1QUFCc8az9hwwZEn09ZMgQweYzv/LKK/Lwww8n+DlDQ0MlNDQ0+nru3LnRn9vKvn37BJCFCxeKiMiNGzckR44c0qVLF7uxjh8/Lt7e3jJmzJgE53Q3wG6J5/vYrHzioz8QhN6C+xZ9BoTl9VtLfbClXWyyAa2Br9ErnwHAHqAyUBv4DO2gmhAmwoLBlSRVRYSGOqdSnGmfQsqXL4+3tzd58uTh5ZdfJiwsjK+++ir6fqtWrezOeHbt2kV4eDgdO3a0G+e5554ja9asbN682a6+adOmZM+ePfo6ODiYkJAQdu7cmWLZH330Ufbv389rr73G+vXruX37dqJ9WrduTUBAAHPmzImumzNnDjlz5qRFC+1NsnPnTq5fv05YWBiRkZHRpXjx4pQvX54tW7akWHZ3YZRPfCh0tqC56FgJDYBQy+sflvo3Le0Swg9tFTcLrYjeRcdmqAbUAD4G/nLQz0RYMBjsWLp0Kb/88guHDx/m1q1bzJ49mzx5Yqx8Che2D7Zo3dKKXZ81a1by5s0bZ8urYMGCceYsWLAgZ8+eTbHsnTp1YtKkSfz00080btyYPHny0KZNG06ePBlvn2zZsvHMM8/w9ddfIyLcv3+f+fPn065dO/z8/AC4ePEiAA0aNMDb29uuHDhwgP/++y/FsrsLY+2WGIXRCuNdF4zlg863+hQwGb2SWQzUAoqhzbfbAoHoyAnrSDzCQkN0jAcT5NSQwalYsWK0tZsjYlu2WRXThQsXqFChQnR9ZGQk//33n53iAvjnn3/ijPnPP/9QtGjR6Gs/Pz/Cw8Pt2iTlC14pRa9evejVqxdXrlxh7dq19O/fn/bt2/PTTz/F2++FF15g1qxZbNu2jTt37nD+/HleeOGF6Pt58+YFtC+U7We0EhgYmKhsnsIoH0/hjV5FNUAnAd+GVkShaEOH+sBhEs5KZCIsGAzxEhISgo+PDwsWLODJJ5+Mrl+4cCGRkZHRlmVWvv/+e27duhW99Xby5El27drF22/HeJwHBQVx8OBBu36rVq1ySq7cuXNHK50pU6Yk2LZ+/foUK1aMOXPmcOfOHYKDg6lbt270/dq1axMYGMixY8fo3LmzU3J4GqN80gJZ0XEc6qFXM9WA/OiE4jmIWRFVdtDXRFgwGBySJ08e+vfvz0cffUT27Nlp2rQpf/zxB4MHD+axxx6jWbNmdu39/f1p1KgRAwYM4N69ewwZMoQcOXLQt2/f6DbPPfccXbt2pW/fvjRv3pxff/2VmTNnJipLz549CQwMpFatWhQoUICjR48yZ84cGjVqlGA/Ly8vwsLCmDJlChEREfTt29duhZcjRw4++eQTXnnlFS5dusRTTz1Fzpw5OXv2LJs3b6ZevXo8//zzzj24VMIon7RGFiAnMB6tiH5Gr4iaAu3QZzy2W3GxIzEYDCkhTx7nQubkSdue1SNGjCB//vxMnjyZiRMnkjdvXjp16sRHH32El5f9nnanTp3Inj07r776Kv/++y+PPvooCxYssNue69y5M2fOnOHLL79kypQp1K1bl6VLlya4HQhQp04dZsyYwZw5c7h27RpFihShY8eODBs2LNHP8MILLzBq1Kjo97Hp1asXxYsX55NPPmHevHlERkZStGhR6tatS9WqVZPwlDyDEhdYoaQXqlevLk4nk/MEj6Ot2mwVyxW04UJJ4CubexHorTt7wx2DweAEVifTDz5wbaSuzI5Sao+IVHd0z1i7pUUcRVjIDaxFh/FpDVgtNVdY2hsMBkM6wiiftEh8ERayAcvQiqgxWhFNtLQ3GAyGdIQ580mLFEabY79B3IR23mifoT7AA+jgpoVSWT6DIYORmY4f0gpm5ZNWSSjCwlK0o2p1YA5w3BMCGgwGQ/IxyietkpQIC2vQoXvqogOcGgwGQzrBbLuldRKLsNAbyEvMCumxVJLLYDAYUoBZ+WQEnkVvv7UGnHO2NhgMBo9glE9GoRHaPLsbekvOYDAY0jBm2y0jEQJsAJqgzbBf96w4BoPBEB9G+WQ0HgK2os+A/gWGkXjaB4PBEefRaeHXo/+GBP131R0TRd2QYsy2W0YkCB0l+zu0H9B9z4pjSGek0USGPXr0QCllF+gztalXrx6PPZY6Vj3BwcFxEuHFZujQoXFSSaQXjPLJqBQANqHNtMOA8ARbGwwxpMFEhnfu3OGbb74BiA6eaYDu3bu7JNOqJzDKJyOTA+0LdA8dlPSWZ8UxpAPOoxMZxo6sYYs1keFqS/tUYNmyZVy/fp2mTZty8eJF1qxZ45JxRSROcjh3ce/ePZePWaxYMUJCQlw+bmpglE9Gxw9YBBQFngTSblZdQ1pgOnqrNrFvBttEhqnArFmzyJ07NzNnzsTf359Zs2bFabNkyRJCQkLIli0buXLlol27dpw+fdqujXUr66uvvqJ8+fL4+PhEJ4Nbs2YNtWrVwt/fn5w5c9KqVSuOHDniUJ7ly5dTsWJFfH19KV++fPSqzIp1O+zgwYM0btyYgIAAnn32WQDWrl1L06ZNKVy4MNmyZaNixYqMHj2a+/cT3h+/f/8+PXv2JEeOHKxfv95uHluUUgwePJhx48ZRsmRJAgMDCQ0N5dChQ3HGGzx4cLQcTzzxBIcPH0YpxdChQxOUxSWISKYp1apVk0xLlIgMEJGHRORvD8tiSLvUFZHwJLYNF5HH3SiLhbNnz0qWLFnkpZdeEhGRDh06iK+vr1y+fDm6zaRJkwSQLl26yKpVq2TBggVSvnx5CQ4OluvXr0e3CwoKkiJFikiFChVk3rx5sn79ejl27JisXr1avLy8pEGDBrJ8+XL5+uuvpXTp0pIvXz75+++Y/zChoaFSsGBBKVGihHz11Vfy3XffSbNmzUQpJRs2bIhuN2TIEAGkVKlSMmLECPnxxx9l48aN0bJ++umn8v3338uGDRvk448/loCAABk4cKDd5w4KCpKwsDAREbl9+7a0bNlSChQoIHv27Ikzjy2ABAUFSaNGjWT58uWyaNEiCQ4OltKlS0tERER0u0GDBolSSt566y1Zu3atfPTRR1K2bFkBZMiQIcn817IH2C3xfB97XCGkZsnUysfKKBEJFpEjnhbEkKp4/H9f8kUfNWqUALJjxw4REVmzZo0AMmnSJBERuXHjhuTIkUO6dOli1+/48ePi7e0tY8aMia4LCgoSf39/OX/+vF3batWqSZkyZey+nI8fPy5Zs2aVvn37RteFhoYKIDt37oyui4yMlHLlysljjz0WXWdVCp9//nmCny0qKkoiIiLkgw8+kFy5csn9+/ftZA0LC5PLly9LnTp1pFSpUnLs2DG7/vEpnzJlykh4eMyviEWLFgkg27dvFxGRy5cvS/bs2aV37952fUePHp1qysdsu2U23kKn3K4H7PWsKIZURJJY6qKNU5LSNhyd+DCpYyeTWbNmUbZsWWrVqgVAgwYNKFKkSPTW286dO7l+/TphYWFERkZGl+LFi1O+fHm2bNliN15ISAiFCsWEgr916xZ79+6lffv2ZM0a431SsmRJ6tSpw+bN9pkaixcvbnfOkiVLFtq1a8fPP/9MVJR9HpTWrVvH+Tznz5+nV69eBAUF4ePjg7e3N4MHD+bq1atcvHjRru25c+d47LHHuH37Njt27KB06dJJemYNGzbE2zsmG2WlSpUAorchDxw4wK1bt2jXrp1dv7Zt2yZpfFdglE9mpBs6X1ATtEWcwWDFUSLD+EiFRIa7d+/m999/p02bNly9epWrV69y48YN2rRpw65duzh69Gj0F3aDBg3w9va2KwcOHOC//+wPOgsXtndSunLlCiISpx6gUKFCXL582a6uYMGCcdoVLFiQ8PBwLl26lOBcUVFRtGjRgu+++47BgwezYcMGfvnlFwYNGgTA3bt37dr/9ttv/P7777Rv397hvPGRJ1Z6c19fX7vxz5/XliIFChRI9LO5C+NkmllpDeRCx4WbCrTypDCGNEN3tB9PKxL+aRqFTmT4tXvFsa5uRo0axahRo+Lcnz17drTfzcyZM6lQoUKcNoGBgXbXsQ/oc+fOjVKKCxcuxOl74cKFOF/k//zzT5x2//zzDz4+PuTPnz/Buf766y92797NnDlz7Hx4Vq50rPGbNGlClSpVGDhwIH5+frzxxhsO2zmLVSlevHjR7pk5+mzuwiifzEx9tLlsc+AK0MWz4hjSAAklMrQShQ7d9BRuTWQYHh7O/PnzqVmzJiNHjoxzv2/fvsyZM4cBAwYQGBjIsWPH6Ny5s9PzZM+enWrVqrFo0SKGDh1KlixZADh16hQ7duzgtddes2t/5swZdu3aFb31dv/+fRYtWkSNGjXw8kp4M+n27dsAdltiERERfP11/Fp8wIABZMmShT59+hAVFeUSJ9tKlSqRPXt2Fi1aRP369aPrFy1alOKxk4pRPpmdauitt8bocDwDPCqNIS3QH+1A2hBtTt0C7WAagd5qm4hWPP3dK8aqVav477//GD16NPXq1Ytzv1evXvTu3Zu9e/fyySef8Morr3Dp0iWeeuopcubMydmzZ9m8eTP16tXj+eefT3Cu999/n2bNmtG8eXNefvllbt68yZAhQ8iZMyf9+9t/0IIFC9K+fXuGDRtG/vz5mTRpEkePHmXSpEmJfqYHH3yQoKAgBg0aRJYsWfD29mbMmDGJ9uvXrx9ZsmShb9++REVFxZHJWXLnzk2fPn348MMPCQwMpEGDBuzdu5cvv9S284kpUZcQnyVCRizG2i0BzojIgyLylmizbIPhnIgMF21ObS3vW+pTgZYtW0pgYKDcunXL4f2rV6+Kv7+/dO7cWUREVq1aJfXq1ZPAwEDx9/eXMmXKSJcuXeTQoUPRfWzNl2OzevVqCQkJET8/P8mRI4e0aNFCDh8+bNcmNDRU6tSpI8uXL5cKFSqIj4+PPPDAA7JgwQK7dlYrNFvrOSv79u2TOnXqiL+/vxQtWlTeffddmTZtmgBy4sSJBGWdMGGCKKVk1KhRdvPYAsigQYPs6k6cOCGAzJgxI7ouMjJS/ve//0nBggXFz89PQkNDZfv27Umy0ksqJGDtpvT9zEH16tVl9+7dnhYj7fIf0BSoBEwGLmECSxoMmYjFixfTrl07tmzZQt26dVM8nlJqj4hUd3TPbLsZYsiLDh7ZGngYyAe8BrxDzLbLSvSBtHXbJX3GNDQYMj0//fQTq1atombNmvj5+bFnzx5GjhxJSEhIqgRPNabWBnsC0GF4rBafDUgTgSUNBoNrCQgIYMuWLXTq1IkmTZowduxYnn32Wb7//vtUiZRtVj4Ge86jlcvv6FWP1SLO1h3AGliyITpittmCMxjSHRUqVGDTpk0em9+sfAz2WANLegOT0GdAddErHVtSObCkwWDIWBjlY7BnHTr9AujznPfRSuZx4GKsti0s7Q0Gg8FJjPIx2KOIOeOx8gZ6e60T2sHQSux2BoPBkESM8jHYI2irttgMB24An9jUOWpnMBgMScAoH4M98QWWzArMBz4DdljqUiGwpMFgyJgY5WOwpzs64nWUg3sl0E6nHdCheCZa2hsMBoOTGOVjsMc2sKQjBdQC7YRaA52SwY2BJQ1pg7t3oWtX/WowuAqjfAxx6Q8EobfUviXmbCfCcv2r5b2PR6QzpDLTp8OhQ/Clh83qe/TogVLKJVGd0xNDhw5lw4YNnhbD5RjlY4iLAt4E5qKdTRsAoZbXP4B56EjYIwATKi9Dc/cuzJkDa9fC7NmeW/3cuXOHb775BoB58+YRGRnpGUE8wLBhw4zyMWQyCgPvApttymBLfWlgAvAccM1TAhrczfTp0KkT5MwJL7zgudXPsmXLuH79Ok2bNuXixYusWbPGJeOKCOHh4S4ZKz1x7949T4tglI8hBTyLXg31RJtoGzIU1lVPt276unt3z61+Zs2aRe7cuZk5cyb+/v7RGU5tWbJkCSEhIWTLlo1cuXLRrl07Tp8+bdcmODiYjh078tVXX1G+fHl8fHxYtWoVAGvWrKFWrVr4+/uTM2dOWrVqxZEjR+LMs3TpUurUqUNAQAA5cuSgRo0arFixIvr+hAkTqFWrFnny5CFXrlyEhIREz2Hl5MmTKKWYOHEi/fr1o0CBAmTLlo3mzZtz8uTJ6HbWGGsjRoxAKYVSiqFDhwJQr149h3mOgoODefHFF6OvZ86ciVKKLVu20K5dO3LlykXNmjUBiIyM5KOPPqJ8+fL4+vpSpEgR+vfvHyedtzswyseQMsagt+KmeVoQg6uxrnr8/PS1n59nVj/nzp1j/fr1tG/fnvz589OqVStWrlzJlStXottMnjyZZ555hoceeojFixczZcoUDh48SGhoKDdu3LAbb+PGjXz22WcMGTKENWvWULlyZdasWUOzZs0ICAhg4cKFTJo0iYMHD/LYY49x9uzZ6L7jx4+nTZs2FChQgFmzZrFo0SJat25tpzBOnjxJ9+7dWbRoEQsXLqR69eo0b97c4Wrto48+4s8//2TGjBl88cUX7Nmzh0aNGhERoQ9ad+7cCcCLL77Izp072blzJ927J8/ENCwsjJIlS7J48eLozLAdO3bkgw8+4Pnnn2fVqlW88847fPnll4SFhSVrDqeIL9FPRiwmmZyb+ENE8onIb54WxBAf4PmSXEaNGiWA7NixQ0RE1qxZI4BMmjRJRERu3LghOXLkkC5dutj1O378uHh7e8uYMWOi64KCgsTf31/Onz9v17ZatWpSpkwZu+Rvx48fl6xZs0rfvn1FROTatWsSEBAgrVu3TrLs9+/fl4iICGnYsKG0aNEiut6a3O3BBx+U+/fvR9dv27ZNAJk+fXp0HQ6Sw4noxHahoaFx6oOCgqIT7ImIzJgxQwDp06ePXbstW7YIILNmzbKrnzt3rgCyb9++JH/O+CCBZHJm5WNIOeXR6RWeBW55WBaDQ5xVFePHw4QJzt9LqCSXWbNmUbZsWWrVqgVAgwYNKFKkSPTW286dO7l+/TphYWFERkZGl+LFi1O+fHm2bNliN15ISAiFCsX4CNy6dYu9e/fSvn17smaNCfRfsmRJ6tSpw+bNmwHYsWMHN2/epGfPngnKu2fPHpo3b07BggXJmjUr3t7erFu3zuEWXtu2be1SVtepU4dixYpFr3hcSevWre2u16xZg4+PD23btrV7bo0aNQKI89xcjVE+BtfQCagJvOppQQwpJfZZT2xS8+xn9+7d/P7777Rp04arV69y9epVbty4QZs2bdi1axdHjx7l4kUd8bZBgwZ4e3vblQMHDvDff//ZjVm4sH0OkCtXriAiceoBChUqxOXLlwGixylWrFi88p45c4Ynn3ySy5cvM378eHbs2MEvv/xCkyZNHJ6jFCxY0GGd7Vafq4j9+S5evEh4eDjZs2e3e2YFCuj8KbGfm6sx+XwMrmMC8CgwG62MDOmS2Gc9sbE9+3nlFffKYl3djBo1ilGjRsW5P3v27OismzNnzqRChQpx2gQGBtpdx06Uljt3bpRSXLhwIU7fCxcukCdPHgDy5csHwNmzZ6lYsaJDedesWcO1a9f45ptv7JTU7du3Hbb/559/HNZVrVrVYXtb/Pz8uH79epx6q7KMTezPnTdvXvz8/Ni6davD9kWKFElUhpTg8ZWPUqq4UmqxUuqaUuq6UmqJUqpEEvuWUErNUkqdVkrdUUodVUp9oJTK7m65DQ4IABainVQPe1gWQ7JIbNVjJTVWP+Hh4cyfP5+aNWuycePGOKVq1arMmTOHWrVqERgYyLFjx6hevXqcUq5cuQTnyZ49O9WqVWPRokXcv38/uv7UqVPs2LEj2qKsdu3aBAQEMHXq1HjHsioZb++YkO9Hjx5l+/btDtsvXryYqKiYUCLbt2/n77//jt5iBPDx8eHOnTtx+gYFBXH06FE7U/EtW7bEMbCID+tq7Nq1aw6fm7uVj0cNAIBswJ/AQaAV0BI4APwFZE+kb3bgKHAC6IzOufkWcAdY6KiPMThIJSaLSGURue1pQQzOMn68yIQJrm+bHJYsWSKAzJw50+H9SZMmCSAbNmyQyZMnS5YsWaRXr16ybNky2bhxo8ydO1d69OghX3/9dXSfoKAgCQsLizPW6tWrxcvLS5o0aSIrVqyQefPmSdmyZSVfvnxy9uzZ6Hbjx48XQNq0aSPffvutrF27Vj7++GMZN26ciIgcPHhQsmbNKo0aNZIffvhBZs6cKUFBQVKyZEkJCgqKHsdqcFCsWDF5+umn5bvvvpMZM2ZIoUKFpGzZshIeHh7dtmrVqlKuXDlZu3at/PLLL9HybNiwQQAJCwuTdevWydSpU6V8+fKSM2dOhwYHf/75Z5zP3aFDB8mVK5cMHz5c1qxZI2vXrpWpU6dKq1at5MiRI0n7h0oAEjA4SKnyqAg8A5RJZv83gPu2/YGSQCTQL5G+jdDeJY1i1Y+09M8Wu49RPqlElIg8KyIveVoQg7N07SpSu7ZIaGjipXZt3d5dtGzZUgIDA+XWrVsO71+9elX8/f2jv2hXrVol9erVk8DAQPH395cyZcpIly5d5NChQ9F94lM+IloBhYSEiJ+fn+TIkUNatGghhw8fjtNu0aJFUqNGDfHz85PAwECpUaOGrFy5Mvr+woULpVy5cuLr6ysPPfSQzJ8/Xzp37uxQ+XzxxRfSt29fyZcvn/j7+0vTpk3l+PHjdvNt27ZNHnnkEfH19RVAhgwZEn1v8uTJUqZMGfHz85NatWrJ7t2747V2c6R87t+/L59//rlUrlxZfH19JUeOHFK5cmUZMGCAXL161eFzcoaElI+SZJqhKKVeBj4EjgDlgG4i8q2TY/wI+IlInVj1mwFEJDSBvs3Rwf9ricgum/q3LXIFioid7VX16tVl924TDyZVuAY8gv4p0M7DshgMaYyTJ09SsmRJpk2blmy/nfSAUmqPiFR3dC8lZz4vA5VFpCZQF+iXjDEqoLfcYnMIeCiRvuvRW3ajlFIPKaUClFJPoFdTk2MrHkMqkxN9/vMKcNzDshgMhjRHSpRPhIicBhCRA4B/MsbIA1xxUH8ZyJ1QRxG5CzyG/gyH0Hk2fwS+Ix6D30uXLtkdqCV0cGhwAdWBQUB7IPOFzzIYMiVTp06N/o4F8sXXLiWm1n5KqUroGMgA/rbXIvJbCsZOFKWUH/q3dQHgBeA0OsvMe+gzn96x++TPnx+z7ZbKvA5sAN5GZ0E1GAwEBweT3COPtE7Pnj2jHXGVUv/G1y4lyscfnUjZFuu1AKWSMMYVHK9w4lsR2dINqIc2VvjLUrdFKXUNmKqUmiwivyZBhri0aQPx2Mo7JE8eWLIkWVNleBQwA3gY/a/VwqPSGAyGNEKylY+IBLtg/kPoc5/YPITOJJMQlYArNorHys+W1wfRac+c5/Jl2LQp6e0dRJY12JAHmI/OgFoVnY7bYDBkahI981FKXVZKPWJ5/5VSqqQL518BhCiloldJSqlgoA5xV1WxuQDkVkqViVVf0/Lq+vgUhuRTG+gLdCAmM6rBYMi0JMXgIDvga3n/IpDfhfNPA04Cy5VSLZVSLYDlwBlgirWRUipIKRWplHrPpu9MtJHB90qpzkqp+kqpAcCnwB7AsUuxwXO8hY6CMMTTghgMBk+TlG23U0APpZRVAT1sOex3iIgkORSqiNyymEePAeagTwh+BPqIyE2bpgrIgo2yFJGTSqkQYCjwAdqq4gwwFRghIlEY0hZe6H9l6/lPI49KYzAYPEhSlM9I9CqkM9qQYGI87ZTlfhZnBLCYaz+TSJuTxFjV2db/jg7k735+/x3On4c6deKPuGhInALAXCAMvT6NG0jYYDBkAhJVPiLylVJqNfAAsBFtPPuHuwVLc5w6Be+/DwcOQK1a0LAhNGgAVap4WrL0R3106u2OwFqc/LliMBgyAkmydhOR88B5pdQsYJWInHCvWGmQp57S5do1bQm3fj106AD//QdK6Tj0DRtCUJCnJU0fvAs0QAdCetfDshgMhlTHqQgHItIlUyoeW3LmhJYtdTrHw4dh717t57NxI9SoAWXLQu/e2u/nSmKuSpmYLMDX6E3czR6WxWAwpDommVxKKV4cChWCr7+GqCg4eBDWrYNp0+DFF+HBB/X2XIMGULs2+PomOmSmcXItgnZADQP24Vo7SoPBkKYxyseVeHlB5cq69O8P9+7Bzp16i+7tt7XRQp06MedFlSrpPrHJTE6uTdDKpzM6Kp/H0xsaDIbUwPxXdye+vloxfPAB/PQTnD4NvXrBX39Bu3ZQuDA8/zx89RWcOeNpaT3HB+hgSqM9LYjBYEgtzMrHEXnyOLeasOR4T5TcuaF1a11AW9CtXw9r18LAgXqchg3h33/h6lXIlctJwdMp3sACdFjYukCIZ8UxGAzuxygfR6TW+UlQEHTrpktUFPz2mz4vOndOnyVVrKi35xo2hJAQ8PFJHbk8QRDam+w59PlPggk1DAZDesdl225KqfeUUj0Sin5gSAAvL6haFQYM0GdGly7BiBEQGanPj/Lnh48/1koqo9IKaImOV54xo80bDAYLrjzzGYr+7XpKKTXQheNmTvz84Ikn4KOP4Jdf4NdfYelSaNFC+xZlVD5GB3T6AjgPDAceB0Itr+9b6g0GQ7rGlcqnJFAR7TJYyYXjGgCCg2HLFm26/fDDsGOHpyVyD77o85+B6FVQRXS0v82W1wroyAifYlZHBkM6xmXKR0ROicjvIjJVRDq6alyDDd7e8Mkn8MUX2mjhk08y5jbccnSC9MvoKAjelnpvoA2wDr06MtZxBkO6xZhap0eefhp+/hm+/VZvw0VkoAQ554HVlvIE8BJxVzhewFhLG7MFZzCkS5xWPpbcOWuUUr8rpY7HKrGzihrcRVCQ3oYrXx727NHOrBmB6cAr6L/Mz4EDwFcO2nkBLwNfpppkBoPBhTilfJRS76IDohQB9qN34m1LknP5GFyAjw98+qmOJ9eqlX6f3rfh1gFPW95nAxYCbwMHHbRtYWlvMBjSHc76+XQDxopIX3cIY7DgrJNr2bLw3XfQvj1s3gwzZ0LevO6Szr0oYs54AB4CRqFD8OyOdc/2vcFgSFc4q3zyAivdIYjBhuQ6uW7ZAu+8A488AgsW6LxD6Q0BIrBXLF2Ab9AWbu/Y1Gegoy6DIbPh7JnPZsBkT0ur+PjA6NEwbpzehhs9GiSd2SM3JO7PGwVMRlu3HbWpX2FpbzAY0h3OKp8+QBelVCelVD6llFfs4gYZDc7SsqUOZPrNN/q9M+kZPE13tINp7KOrYLQHWU/LvSh0LqDuqSmcwWBwFc4qi6Not78ZwD/ojQ/bEu5S6QzJJzgYtm7V50EPPwy7dnlaoqRRGHgKeIO4CuhV4A4wDZ3M/SmgUKpKZzAYXISzZz7DMX7l6QfrNtzjj2t/oIEDoV8/nfY7LdMfvcXWEG1O3QJ9BhQFtAdeQ1vA9feUgAaDIaUoSW9nAimgevXqsnv3bk+L4RlOntTWcAULamu4pKaB8CTn0X4/623qGgL/AqeBdJi81WDITCil9ohIdUf3zBlNZsG6DVemjLaGSw/bcIXR5zy2nmSDgZHA7xjlYzCkY5IT4aCwUupTpdQvSqm/LK8fK6XM7ntax8cHPvsMPv9cGyJ89ln6s4YD8EOf+7wGXPWsKAaDIXk4G+HgAXRkg9eBm8DPltc3gP1KqbKuFtDgBlq10tZwCxfq9+nJGs5KXfRZ0FueFsRgMCQHZ1c+o4DrwAMiUl9EOohIfeAB4JrlviE9YN2GK11ab8P99JOnJXKekcD3wCYPy2EwGJzGWWu3+sBLInLStlJETimlhqI9LwzpBes2nNUa7u23tUJyZiWUJ0/qpR2PTU60T1BP4FfA3zNiGAwG53FW+fgAN+K5d8Ny35DeaNUKqlTR1nDHj8Off0Lu3Enr60wMOnfQEpiLznD6oWdFMRgMScfZbbf9wGuxIxkopRTaI2O/a8QypDolS8K2bTp9d3rbhhuPNsn+1dOCGAyGpOKs8hmOzi35h1JquFKqt1JqGHAI7YExzNUCGlIRHx9tij1mjE5Y9/nn6cMarhD6/Kc7EOlhWQwGQ5JwSvmIyBqgOXqLbRB6x30w2uKtuYisdbmEhtTHag03b55O133liqclSpwuQA5gnKcFMRgMScFpPx8RWWPxWA0EigOBIlJDRH5wuXQGz2HdhgsO1ttwP//saYkSRgFT0Oc+xz0si8FgSJRkRzgQkdsiclZEbrtSIEMawsdHb7199hk0awabNnlaooQpg/b76YWJQGgwpHESVT5KqftKqRqW91GW6/iK2XHPiLRurdMzPPss7N/vaWkSph/wHzDb04IYDIaESIqp9XDgb5v35jdlZqR+fZg4Ua+AtmzRzqlpkaxoy7cmllLQs+IYDAbHJKp8RGSYzfuhbpXGkLZp2xb+/RcaN4bt23WE7LTII8CL6NSH8z0qicFgiAdnY7sdV0o5TKOtlKqolDJHvRmdl16CF16Ap56C69c9LU38DEVHHvzOw3IYDAaHOBvhIBjwjeeeHxCUImkM6YP33oOLF7VJdlTsdKNphGzAVLQJdijaNtNgMKQZnFU+EP+ZT3VMgPv0T548SQuZIwK//64t4u7fhyxZ3C6a0zyJdn3+HzoKgsFgSDMkqnyUUn2BvpZLAVYqpcJjNfMH8gALXCueIdVxJkjovXvQtCm8+qo2RkiL6bk/ASoCHYDaHpbFYDBEk5SVz3HgR8v7zsBu4FKsNvfQuSWnu040Q5rH1xeWLtWWcMOGwdChnpYoLnmAsejQO/uIf9PYYDCkKkmxdlsOLAfQ8UN5X0SMYYFBkyMHfP89PPaYtn7r3dvTEsWlLTry9UhgiIdlMRgMgJNnPiLSxV2CGNIxBQvCDz9A3bqQLx+0a+dpiexR6CiED6MVUQXPimMwGJw3tR6olHJ4dKuUGqeUGuAasQzpjlKl9ArolVdgwwZPSxOXYmgX6e7AfQ/LYjAYnI7t1gX4LZ57+y33DZmVKlVg0SJ47jnYu9fT0sSlF5AFmORpQQwGg7PKpwTwZzz3jpOB/Xzu3oWuXfWrIQFCQ2HKFGjeHI4d87Q09ngB09AOqKc9K4rBkNlxVvncBorGc68Y2uotQzJ9Ohw6BF9+6WlJ0gGtW2vrt8aN4fx5T0tjz4PAG0BvTJRCg8GDOKt8tgIDlFJ2BquW6/6W+xmOu3dhzhxYuxZmz/bM6ifdrbx69IAuXXQYnmvXPC2NPQPRK5+FnhbEYMi8OKt8hgJlgaNKqRFKqZeVUiOAo5b691wsX5pg+nTo1Aly5tRhzTyx+kmXK69Bg7QFXMuWaUtr+qA90vqi0y8YDIZUR4k4t/dgye3zKdpf3AuIArYBb4rIbpdL6EKqV68uu3c7J+Ldu/oYY/Nm8POLe50aWOdcuxYaNUrduVNMVBQ8/zyEh2tjhLQUhqcPOiDUTI9KYTBkWJRSeyyZr+Pec1b52AzqD+QGrojInRTIl2okR/lMmKCjxjz5pF55BARolxY/P72zFBCgi5+f+6LLWGV45RX79+mGe/e0AUKpUjB5ctoJw3MTHXpnGjoGnMFgcCluUT7pEWeVj+0qZ/16+OoruHlTH2EcOAD588OtW7ouMlIrocBAx68J3XPUJiAAsmZNOyuvl1/W4duSPeeNG/DEE/oMaPhwl8qXIlYDrwAHgOwelsVgyGC4XPlYcvqUQ6dRsENE0mwCY2eVT0KrjNj3wsO1IrpxQysj66vt+4TqHN3z8dEKyNcXihaFIkUgOFgbkOXNC716QcmSOqiAOxcTEyZog4tOnVK44rp0CerUgddf18FI0wod0RlPR3taEIMhY+Ey5aOUygWsAkKsVZbX6EFExKlNfaVUcWAMeuNDAeuBPiKSJE8MpdSDaN/1+ujfrqeBiSIyNnZbZ5RPYisMd69ARODKFWjQAObO1TtX587BiRPafWb2bAgKgtOn9b3g4JhSsqT9+9y5k6+cXH7edPKkjgP32Wfw7LMpGMiFXEJvv30HPOphWQyGDERCysfZfD4fAnmBx9Fm1a2Ba0BXoBbwnJOCZQM2oP2DOqOV2AfARqVUZRG5lUj/6pb+m9CBU66hre4CnJHDEVYLt/i+aP38Yizf3HH+ohTMmwfdusFDD+m6hx+OuV+mTMzK6/p1/Z1+8qRWTidPwtatMdci9koptnLKmTN+ORxZ+qXo8wYH6zA8DRvq3EENGqRgMBeRH73q6Y6O2e7tWXEMhsyAsyufv4BhwNdABPCoiOyx3JsEZBeRTk6M9wbwGVBORI5Z6kqioyi8JSKfJdDXCzgIHBGR1kmZL6krn6Suaty5+nHlyuvKFXvFZPv+xAnw9naslAoX1lmzt21zw3nT1q3wzDNaEVV3+MModRHgKXTW03c8LIvBkEFw5bbbbaCRiGyzvG8qIpss9xoCC0QkrxPj/Qj4iUidWPWbAUQkNIG+T6DzDD0uIklybk2q8nHGosxd1mdJGdcVc4vAf/85Vk5798LlyzprQsmSUKGCNrYoUgT+9z/9mqKzphUrtHbbvBnKlo23mUsMHpLCSXQ+3h3AA26cx2DIJLhy2+0CkMvy/hR6q22T5bpMMmSrgCVXUCwOAYnF5X/M8uqnlNoFVAOuoLOpDkyJ+fe+fXD4sHZLSYyICChfPrkzOcYaUWHz5oTbde+uVyLduiX/S1kpbbCQL5/9AsS6yjlxQm/r/fWXNjXfv1+fQS1apD97xYq6VKoU8z5PniRO3qKFNkJo3Fgvr4oUcdjM1sHWrSbmwcBgoCd6M9dZF2yDwZBknFU+29DGBt8Bc4AhSqlgIBJ9ZrPCyfHyoBVGbC6jfYgSwvpNtRCYALyN/t06HCiOPo+y49KlS1S3+Ybt2bMnPXv2jDOwp6MIJHbeZMWd505WGfz9dSlYEGpb0lA/+KBWWu3awcGDuliV0sGD2kw8tkJ66CHI7siUuVs3uHgRmjSBLVsgVy6727ahjRo1SpmiTRKvAfOBL4EebpzHYMigTJ06lalTp1ov88XXztltt9JAERHZqpTyRueGbA9kA9YAr4lIkgOWKKXCgc9E5O1Y9R8Ab4tIvMpRKTUV/fUwXkRet6kfaJHrIRH5w7ZPcpxMPUG3bnrl5Z2Eg2/rysuVCjMl500i2gLPqpSs5fBhvbCxVUgVK8IDD4CPt0Dfvnqf74cftLaz4BEH2wPAE8BaYCXa/lKhz4Uaog0TCrtZBoMhA5BmnUyVUv8Ay0SkV6z6iUA7EcmfQN+P0KudFiKy0qb+YWAvECYi82z7pBfl42nccd4UGalNxGMrpVOnoHRpqFRRqPjHIir6HaPi7LcoWTYr4eEecrC1KpnfgMnA02gLuAi0MvoCbZzQnxhnA4PBEAeXnPkopXzQW1xjRGSLi2Q7hOOkxg8Bvyehb0JEJUuiTI67zpuyZtUrtPLloW3bmPo7d/Sq6OBBxcFizzB1/AYOPPQPl8lDXu9rFPH5j4kVNlE54DiPBBzlhcsN+LKS8ErRZXqAPHlgyZJkf16HjEYb7J+xXHvbvLYBWqHTMowG3nTt1AZDZiHJykdEwpVSDYA4zpspYAXwqVKqlIgcB7CcIdVBr2oSYjXaP6gx+veolSaWV7PESQapfd7k76/9lx5+GHghC/w0HMLD+SekBU+sH0THlwtx5EgFlu3X50q5c+sIEP+90IeQEKg2rAVJNq9MCufRf1nrgO3Ac3C3Nrz8PxtrOy/0/4KGQBhmC85gSAbOGhxsRxscbHLR/NOAV4HlSqnB6A2P99G/OadYGymlgoC/gOEiMhxARP6zbL29q5S6jrZPqo5O6zDL6jdkcA5PW/qRJQusWsWiCpN4+fFNvPJKvehbUVFw9Ch89BFs2gQbNsC+n+aROxiqVbMv+eI95kyE6ehYb15AXaAFTH8GDkXGUrRewMtow4TByZzLYMjMiEiSC3qL7E+0wigGZEH/N4wuzoxnGbME8C1wHbgBLAOCY7UJRiumobHqFdAPOAaEo82/hwPejuaqVq2aGNI4oaFy545IjSp35U7R0iLz5sVpcueOSI0a+vX+4/XkyBHdrH9/kXr1RHLkEClRQqR1a5EPPhBZvVrk4sUkzl9XRMJt5rogUsNb5OrKmDmjCReRx1PwWQ2GDA6wW+L57nfW2s16jhJfJ5EELNQ8jTE4SAfUq8eEtpu0MUPoQZ3LYs4cbWdtQ7TBw6J6ehlkQ1SU9kvasyem7N2rnWVjr5AKFIg1fyhgc941YQKoA/DKRpjQG5RPrG3GWO0NBkMMrnQyHeYCeQyGeLkb5RNj8OBXURsTtG4N330HNWpEt4s2ePD1iRNa3ctLB0woWxaes0QbjIqC48djlNGnn2qFlD17LIV0DwpFAN6xjC86QvcLELrJxsgiIlUeicGQIUlU+SilLgMNRGQvEAS8LyIn3C6ZIVMy/XwzOvWxMXioU0cftrRpA7/8ogPOYWPwMLYpSbF38PLSwVjLlIH27XWdiL1CGjMG9v4G/oWgWh1tHv7IIzo2XuHx4FcZXuhuc/azApOEzmBIJklZ+WQHfC3vO6M9H4zyMbiFfTfLcHhebIOHp8G7EpS/AFULgtJxbyIioPzN+GPCJYZS2seodOmY7A5yDk60hV3Pwptv6XQVFStCtmxQKwiqzYXZgRDWAXJNRIfYNRgMTpMU5XMK6KGU8kUf8D+slIrXEFdc5wNkyIR8We4T2NQs7o2oEtDqdR1ye9y4mPp6nwAO2icTVQRKtYHvp8Gg/8Err+oV0l9/wc7tsOsduH4OCuSHR4pByEioVQtCQqBECfck9Uu1wKoGQyqSFOUzEm32bM23MzGedtYAJE4lkzMYkoSXl86g9+ijOohcx45um+ruKzBnHGxeDBQG1cKyZRcEL1yD0f2gbgEYPk2bps+fr5OzZsmilZBVGVWvbhcpKNmkWmBVgyEVSVT5iMhXSqnV6CDzG4HXgT8S7mUwuIFcuWDpUqhfXweJq1LFLdNM/xI6DQS/Nmi/H5uFFg3B703ovBL+PAoDB+pqEZ2KYtcu2LlTbxseOqSDsNaqFaOQSpZ0bnWU6oFVDYZUwllT6xloR890eeZjTK3TAfXqxTGdjsOCBTBokDZAaNMm8fZOkKT4ceFw92EIvQebD8bf7s4dbVG3c2eMUoqMtF8dPfpoPNG+LXgksKrB4CLSbGDR1MYon3RAmzY6g11i/PUX3L4Njz+uV0MuIslf8D/BhIagBsMrbyVtbBH4+297ZfTbbzqyt1UZ1aoVkyI9tiJMtcCqBoOLcKWfj8HgXpIaJDQiAho0gMqVXTq9U+GFAqD8NCCJykcpKF5cF6t13b17es6dO2HVKhg8WOvUkBB9v2ZN/VH9/Nybv8lgSG3MyseQfvnnH32qP3EiPP106s9/A6iIju/WwHXDnj0LW7fqFEelSsGvv+oYevXqabenESN04lez+jGkdRJa+ZhEwYb0S8GCeonSrZtOFpTaBKJtP3sBt103bNGi8O+/ehW0fTv89x98/jnkzKm3BQ8cgHLlYMAAvVq6ft11cxsMqYVRPob0TUgIDB+uQ/DcupX68zcDagBDXTek1cKtWzd97esLjz0G774LP/6os45nz67NuD/7TCurGjXgrbfg++9dq4zu3oWuXfWrweBKjPIxpH969dLbb92761P91GYsMAudP9cFJJZTKWdO7XRasKBWRpcu6Vh1AQH6tUgRrYwGDoTVq3X+o5TIYvUxMhhcibOm1ta0CZE2dY3RO98bRGSf60V0HebMJwNz545eHrzwAvTpk/rzz0T7A/1Misx4kmrRllC7u3fh55+1BfrGjdoivUIFfWZUr55+TIGBSZfF6mNkrOwMzuIyU2ul1ELgnoh0sly/REzEgwigmYisT6G8bsMonwzOyZN6G27hQv2tmZoIOshoE1KUWtsZX56ktr17F376SSujTZu0MqpYMUYZ1anjWBkZHyNDSnGl8jkFDBSRBZbrv4Afgf7AVKCQiNRPucjuwSifTMDatfDii/obtmjR1J37L6AmevVTKnlDdOumTb29vRNva80k6+yW2N272s/Iqox2745RRvXra2WUNavxMTKkHFcqnztAIxHZqpQqAxwFqorIb0qpRsA8EUluAmO3Y5RPJuHDD2HlSv3N6uubaHOX8jGwDliLjnaYDrhzx35ltHs35M+v8yH176+VUUCAWf0YnMeVptbXgbyW9/WAf0XkN8v1fYiT18tgSH3efhsKFdKOMqlNP+BfYE7qT51c/P31qmfoUK18zpzRq5tHH9V6vHBhaNhQK6mpU43lm8E1OKt8dgBvK6WaA32A723ulQH+dpFcBkPy8fKCWbO0KdisWak7d1Z0MNIBwMXUndpVfP01vPqqdmbdvBnOndPWdUeO6GO1EiW0Xl+/HsLDPS2tIb3i7LZbWbTCKQ0cR2c4PWm5twE4JSJd3CCnSzDbbpmM33/XBxU//KBTkqYmbwLnSXfJ5hI727lzR6+I2rTRx2uHD8MTT0CzZtC0aXSiWYMBcOG2m4j8KSJlgfwiUsaqeCy8QZKjXBkMqcBDD+nQO888o8MEpCbDgJ3A6tSdNqXE9jFasGABxYoVw/oj1d8fXnpJ+xjt2gV//gmtWmlF9NBDUK0avPeePkOKivLc5zCkfZxSPkopL6VUVhH5z6ausVKqP5BVRC65XEKDISW0awdt28Lzz8P9+6k3b3Z0wvnewM3UmzYlxI6sICKMGjWK8+fPs2LFiuh23bvrvH5372rDhE6dtHX7xYs64oI1KkKhQtrwcNEiuHYt+TKZCAsZE2fPfOYDX1kvLH4+q4FPgF1KKReGVzQYXMRHH2m75PfeS915GwGPA++m7rTJJfaqZ9OmTfz777+ULl2aoUOHRq9+bKNr2+LtrbfsPv5YR0X4+We9RTdjho7kXb++jsDwxx9JD0RhIixkXIyfjyFzcPGiDsEzfjy0bJl68/6Ljv+xAh0DLg1j62MkIvz6az1y5gzl8uUfgHCCgoaSL59+ds76GN2+DRs26ECoq1ZpP6JmzXSpV8/x+ZKJsJD+SejMBxFJcgHuAHUt78sAUUBly3UjtOm1U2OmZqlWrZoYMjE//SSSP7/I4cOpO+9cEakkIuGpO21K2LBhg5QpU0a2bt0qNWvWlGXLlknVqlUlKioqxWNHRYn89pvIRx+JPPaYSGCgyNNPi0yeLHLmTEy78eNFJkyI+96QfgB2Szzfx8bPx5B5qFEDPvhAm2rdTMWDmOeBIsCnqTdlShARhg4dyrvvvkvWrDpQXYsWLQDszn6Si1JQqZJ2x9q6VZtvd+ig31etClWq6AjdEyfqMyOwP2cyZAyMn48hc9Gjh85V3bVr6kXAVmjjg9HAn6kzZUrYvn07586d4/nnn4+uU0oxdOhQRowY4fL58uTRymfuXJ0fcNIknUDv+nUICtLbgT/9BB07mrOfjISz8XffQiucFWg/n6E299qjjUsNhrSLUjpOTN26MHo0vJmCKKDOEAwMAnoCG0jToXeCg4P55ptvolc9Vlq0aEGWLFncOneWLNol6+pVnR/w0iVtSffqqzo1RGSk9it68EG3imFIBYyfjyHz4ecHixdr06uNG1Nv3tfRZtdfJdbQsxQrVoyHH344Tr1SiubNm7t9fluru+LF9e+D336DJUugTBm9e/r443oVZLK4pl+SlUxObPx8bOoOiPHzMaQXgoL0Ps/zz+tgZqlBFnTonXeAC6kzZXojtq+RFaX0imjNGr3qee01bTVXogSEhWmLuNR04zKknGSlvVJKVQHK4cDAQERmp1QogyFVaNBAJ55r2xa2bEmdCNhVgG7oVdA37p8uvZFYFlc/P33/4kW9Evr3X1iwAP73P7hwQZ8Lde5stuXSA85GOMillNoO7EM7nM60lBk2xWBIP7z1lt7bef311JvzPfT/oJQbjmUo4lv1xMbW8i1fPn0etHu3DuEXFQVPPqm35r74IvWjKhmSjrPbbh+iTa3roo9MWwNPoMMnHifNu9EZDLFQSrvgb9mSeqZU/miX7FfQzgsGIPFVj5X4IixUqKCjK5w+DcOHw7ZtULq0Du23YoV2jHUWE97HfTi77dYYHTJxl+X6bxHZA2xSSk1CGx10cqF8BoP7CQyEpUu1BVyVKjoSgrupj3bL/h8wwf3TpQf27dMRFhYtSrytNcKCI7JmhSZNdLl2Db75RiulHj20SXfnztqfSCXB4tA2vI9JoudanFU+hYHjInJfKXUXsM38vgRY4DLJDIbUpHx5mDJF/0y2pvJ0N5+gQ+88D9R2/3RpHXcsPHPm1EqnRw9tuj17NrRuDTlyaCUUFqYDoDrCug1oDe/TrZsJ7+NKnN12uwDksrw/BdSyuVfGFQIZDB6jTRv907hDB+1Q4m7yAJ8DPYB77p8us1OmjN6OO34cxo2Dgwe1YUKzZnp1FHtrzboNmDOn420+Q8pwduWzDQgBvkMnCh6ilAoGIoHOmCNUQ3rnyBHYswdKldIlMfLk0WZXyaUd+n/SKLQhgsHteHnpYKb16uk4s0uX6vTgvXvDs8/qFVGVKnrVs3mz7tO9uw5yalY/rsNZ5TMMHaUK9KZBXnRkg2xoxfOa60QzGDzAlStw9KjOBfDaa3obLiHq1UvZfAqYCDyMVkTGRDhVCQjQq5oXXtCGCnPm6NXO1asQEqIjLBQvbm/kYM5+XIOzEQ7+EpGtlvcRItJfRIqJSB4Red6R86nBkO7In19HQHjpJZ18xt0URweq6oGOE2/wCCVKwKBBOq5c/vxQoIA2TGjbVi+GTXBT15KsCAdKqeJKqdpKqSdiF1cLaDB4hOrVYeRIfQ6UGjFceqPjwk91/1SGhPnyS3j5ZX3mc/o0PPaYTgHVsqVeDU2f7mkJMwZObbsppUqhfXqs/jxWY0WxvBd0EBGDIf3TrZtOx9mli14JJcU2N7lYQ+/UA54GirpvKkP8WC3crGc92bPrIBi9e+toTB9+qKMrFCmirebc+SeR0XF25TMdKIFOp9AE7a1QH+1oan01GDIO48bB339rRxF3UwG9AnrV/VMZHGPr6CoilCxZknXr1uHrq3+LHD2qt+Fef11vyc2fnzqGkRkRZ5XPo8DrIjJeRNaJyObYxR1CGgwew9dXr3o+/xx+/NH98w0CDqO95gypSuzwPsuXL+fUqVN8+OGH0W2yZNH5hooWhWHDdMK78uW1tdw9Yy7vFM4qn7+BcHcIYjCkWYoXh6+/1lErT59271y+wDS03ehhYDjwOBBqeX0fOO9eETIrsVc9w4YNo3jx4hw8eJB9+/ZFt7Navp09q7OvzpgBy5Zpy/zRo1M3SW56Jjmx3QYqpbK7QxiDIc3yxBPQt6/ec3H3T9w66M3tJ9EREH4ENlteKwAd0Sm5UykRa2Zh3z6YN09bz1eqtJxjxxT//psLP78uNGw4LNo3qF49vd22d6/uV7cufP89fPedPiIsWRKGDjVBTRPDWVPrOej/BieVUiuVUrNjlVnuEdNgSAMMGKBXQW+84d55RgOVLO/zAt6W995AG2AdOr7IaPeKkdn48kvYvh02bhS8vYfx9ddDKVtWsWhRG3x8fmbMmH1s2gSbNul2sSMePPywzrq6fbs+JixbVifCO3fOE58m7eNsSoUX0amwcgGPoKNbxy4GQ8bEGgF70yb96g7OA6uByeiAoz2B2H4lXsBYSzuzBedyli9fjlKKp59+GgA/Pz8GDhzIsGHDktT/gQf0Ft5vv+kEdxUrQq9e8Ndf7pQ6/eHsttswYCk6jXZRESkZqyQhHonBkI7JkUOH03nrrZh9F1cyHZ1qwQudsKQCMMJBOy/gZcDEG3M577//PkOHDkXZ2FH37NmTn3/+mV9//TXJ4xQrBmPGaAu5ggWhZk2dOPe339whdfrD2fA6eYGJInLVDbIYDJ4nT56khcwpUABq19bZUF3JOuBtm+sJ6OynrdF7Dba0AMYBg10rgi0hISHs2rUr8YYZiH79+kWveqz4+/uzYMECcufO7fR4+fLpgKZvvqkDpzdurH2Y33lH/wnFx9272tl14sQMGk9ORJJcgDXAq870SUulWrVqYjC4jP79RRo3FomMdN2Yjzuomysi5UXkVhLbxyI0NFTQ5glxSuPGjaPb7du3Txo3bizZs2eXwMBAefrpp+XPP/+0G+vWrVvStWtXyZ07t5QqVUoWLFgQZ75Ro0ZJ5cqVJSIiInHh0jhVqlSRffv2uXTMO3dEJk4UCQ4WCQ0VWbNGJCoqbrvx40Vq1BCZMMGl06cqwG6JT5/Ed8NhYygH/AqEoVdBXrGLM+OldjHKx+BSIiJE6tUTGTzYdWPWFZHwWHVRIvKciLwaqz5ckqR8Dh06JDt37rQrn332mQDyxRdfiIjI0aNHJSAgQEJCQmTZsmWyePFiqVSpkhQuXFj++eef6LHeffddKV68uCxfvlzef/998fb2lqNHj0bfP3PmjAQGBsr27duT8eHTHu5QPlYiIkTmzhWpUEHkkUdEFi8WuX9f37tzRyueq1f16507bhHB7bhS+URZyv14SqQz46V2McrH4HIuXBApVkxkxQrXjDdcRL51UH9ZRIqLyPc2dYtF5P3kTdO1a1fx8fGR//77T0REunXrJjlz5pQrV65Etzlz5oz4+vrKgAEDouuqV68uo0aNir4uX768TJw4Mfr6mWeeka5duyZPqDSIO5WPlfv3RZYt00qmXDmRr74SGTMmZsUzfnz6Xf0kpHycPfMZjvEuMBhiKFhQZyJr2RJ27NAZy1JCd7QfTyvszYFyA7Ms9/ZjOX1FR1p0ktu3b7No0SKefvpp8uTJA8CuXbuoVasWuXLlim5XrFgxKlasyNKlS/nYEl4oPDwcf3//6DbZsmXjriXM85o1a9i0aRNHjhxxXqhMjJeX/vNp0QI2boQPPtDm2h99BLdvZ9xcQk4pHxEZ6moBlFLFgTFAQ3Rw0vVAHxFxypVcKfU28BGwXUQec7WcBkO81KqlvQrbtIGdO3U0yuRSGHgKeANtTm2rgOqjU273QAcefQqIJwV0QixdupQbN27QuXPn6LosWbLg4+MTp62vry9//fUXd+/exc/Pj5o1azJr1izatm3Lb7/9xv79+5kwYQL37t3jtddeY+TIkeTNm9d5oQwopX2Zf/9dGyRs2waffqoT3mXIXELxLYlSo6CT0P0JHET/1msJHAD+ArI7MU4p4CbwD7AtvnZm283gNqKiRF54QeT55x2fHjs1loh8IiJPiN5as54BhYvIfBEJEJF2lnbJoFGjRlKgQAE7g4B27dpJ0aJFJTw85sDp+vXrkjNnTgHk3LlzIiLy999/S6VKlaINFqxbcsOGDZNatWpJVEo/exojNbbdbLGe9VjPeLZt01txrVuLPPxw+jv7wVVnPq4u6N9394EyNnUl0Wm5+zkxzg/AFGCTUT4Gj3HrlkiVKiLjxrlmvHOiz4Aetynvi8gGEckrIn/G3zU+zp49K15eXtK3b1+7+q1btwogXbt2lb///ltOnjwpzzzzjGTJkkUAuXDhQnTbqKgoOXbsmPz7778iIvLXX39JQECA7N+/X27duiW9evWSAgUKSHBwsIxz1bPwEKmtfByd79y5I/L22yIBASKdOqX8t01qkpaVz4/obbLY9ZuBzUkc43ngEpDHKB+Dxzl2TCR/fpGtW907z+ciUlN0PmEnGDVqlACyf//+OPe++OILyZEjR/SqpkGDBtGGCbYrotg0bdpU+vTpIyIi//vf/6RatWpy8eJF2bdvnwQEBMj69eudEzINkZrKJ/aqJzbbt4v4+4s0aSJy5kyqiJRiElI+ycpk6kIqoLfcYnMIeCixzkqp3OjzordE5LKLZTMYnKd0aZg5E9q3h/NujH3zGpADx9EPEmDWrFlUqVKFKlWqxLn38ssvc/HiRQ4ePMjp06dZt24d586do2bNmnh7ezsYDZYtW8b+/fsZPnw4oI0OOnfuTP78+alatSqNGjVizZo1Tn64zIltVG1H1K6tk9l5eek4ctOmgaRj8y9PK588wBUH9ZfR9j2J8QlwFJiZlMkuXbpE9erVo8vUqSZnscENNG0KPXpoBRQR4Z45vNB/9ROBJAYg2L17N7///rudoUFsfH19qVChAsWLF+fAgQOsX7+e3r17O2x7+/Zt3njjDcaMGUNgYGB0/a1bt6Lf37x507pDYUiA2LmE4uOll+Dff2HNGq18GjaEEydSR8akMnXq1OjvWCBffO2cNbVOMyil6gKdgEckiX/d+fPnZ/fu3e4VzGAAeO89+OUXHQNuzBj3zFEE+IIY8+uAhJvPnj2brFmzEhYWFufe33//zaRJk6hduza+vr7s3r2bjz76iDZt2tChQweH473//vuUK1eOZ599NrquQYMGTJgwgfLly3Pu3Dl+/PFH+vfvn9xPmGlIbNVjxZpLaNcubdk/Zgw8+qj+c3v1Vb0q8jQ9e/akZ8+eACil/o2vnaeVzxUcr3DiWxHZMgUdVvFvpVQuS11WIIvl+o6ImNyCBs/g5aV/yj76qI4o+dxz7pmnLbAK6ItOQhcPERERzJ8/nyZNmlCgQIE49729vfnpp5+YMmUKN27coHTp0rz33nu8EU/6iMOHD/PFF1+wZ88eu/p3332Xixcv0rVrV/z9/Rk5ciSNGjVK/ufLJOzbB4cPw6JFibeNiNDZU7Nm1Vk+WrbUK6ZvvtHm2OXKuV9eV6ASWzQopTY4MZ6IyJNJnlyP7SOx/HKUUpsssoUm0Dex1U5fEfnctqJ69epiVj6GVGX/fr03snGjjq3vDq4DD6Pz+7RyzxSZlapVqzJz5kyqVq3qaVESJCpKByAdOlQrpP79tXLyNEqpPSJS3dG9pCzSvNDOn9ZSHqgHBAP+ltd66LhvykH/hFgBhCilolMxKKWC0bkcVyTSt76D8ivagKE+sNhJWQwG11O1qs6t3KYNXLvmnjlyAHOAlzD5fTIpXl562233bli/HkJC0n7qhkSVj4jUE5H6IlIf7XMdAdQSkVIiUkt0Dp9alvqxTs4/DTgJLFdKtVRKtQCWA2fQ22oAKKWClFKRSqn3bOTaFLsAV4Frluu/nZTFYHAPnTrp1U/nzvonqjuojU481xUTACsTExwMa9dC794628fQoRAe7mmpHOPs8dT7wLsi8pNtpeV6KPCBM4OJyC3gCbTF2hx0pKoTwBMictOmqQKyJENegyFtMGYM/PMPjBrlvjneBf5DW8C5iF27dhESEuK6AQ1uRyl9BrRvn853WK2atn1Jazi7K1gW7dDpiIuA01EVRcdweyaRNidJwpaeiNRzdn6DIVXw8dGnyTVq6MBdDRu6fg5vYC560/oJ4EHXT2FIPxQtCsuXw4IF0Ly5XngPGwY2cWE9irMriRNAr3ju9UJvoRkMBkcUKwbz5mlb2VOn3DPHA2jH0zAgjW63GFIPpaBDBzhwQP/JVamiA5amBZxVPsOAp5VSB5VSQ5VSvS2vB4Fm6K03g8EQH/XqaXOkZ57RnoXuoAdQDBjinuEN6Y8CBWDhQr3r2749vP463LyZeD934pTyEZEFQGPgGvAO2sXtHfRBf2MRWehqAQ2GDEe/flCyJLz2mnvGV8B0dP6fLe6ZwpA+ad1ar4KuX4dKlbRlXHzcvQtdu7rvN5LTB/gisl5E6qDNrAsB/iLymIj86HLpDIaMiFLw1Vc6Y9j06e6ZowDalrQT+qeiwWAhTx4dfnDiRG2Y0KOHYy+A6dPh0CHtuOoOkmU9ppSqgvatbgJ0VEp1shaXSmcwZFQCA2HJEnjnHfeZIjUDmgKvumd4Q/rmqaf0KihrVu3//N13MfessebWroXZs92z+nFK+SilcimltgP7gPno0IYzgRk2xWAwJIXy5WHyZGjXTkeLdAefAr8AC9wzvCF9kyMHTJqkFU2fPtCxo/5TtMaay5kzJouqq3F25fMhOnt8XfTOcmu0UefXwHGghkulMxgyOs88o0+AO3SA+/ddP342tPn162jXbYPBAfXqwa+/asOEihVh7Fh93gPQvbt7Vj/OKp/GaAVkDeL+tyWaQCdgPTozqcFgcIYRI7Tiee+9xNsmh+pAH6Az4KYAC4b0T/bs8NlnEBYGN27o1wsXYiJpu3r146zyKQwcF5H7wF0g0ObeEvQus8FgcIasWbUn4Jw52ivQHQxE+/24KbuDIWNw9672Azp8GB58EMaN0/XuWP04q3wuALks70+hY7pZcTq6gcFgsFCggI6A0KMHHD3q+vGzoANYjUSH3zUYHGA968mVSy/IR1gy5bpj9eNseJ1tQAjwHfpPeYglCnUkelGfWCRqg8EQHzVr6pgoDz8MjzwCWbIk3idPHm01lxRKotMuhAG7gUQSlxkyF1YLt82b9bWIICIopdco3btDaKg2z04s6V1ScFb5DEPnTwSdwjov0B59rLkCnVneYDAkl5w54dln4c4dmD9f+wQlRL16zo3/Avqn4zuYLTiDHbGzqT733HPs37+fI0eOAParn1deSfl8zkY4+EtEtlreR4hIfxEpJiJ5ROR5Efkv5SIZDJmciRP11ttYZzOUJAEFTEZnu1rn+uEN6RPrqqdbN3199epVVq1axenTp9m7d290O1ee/ZgUBQZDWsPfH779Fj76CLZudf34edAeeV3QKRgMmZ7Yq56xY8dSuXJlHnnkEYYNGxbdzpVnP04nWlVKFQb6A6HoP+PLwEbgMxG5kHKRDAYDJUvCrFnw3HM6AkKRIon3cYYGwLPoWPSLcD4HsSFDsW+ftnBbtAgiI6/y00/jKVr0De7dO8nu3WuoVm0vgYGPABARof2jU4pTykcp9QCwFcgNbAeOoeO7vQF0UkrVFZE/Uy6WwWCgSRN46SUdAWHjRp0TyJV8CDwKzEabCxkyLbYrmWHDxlK27NOEhBRg796/qVBhID/+OIzlLnYDcHbbbRRwHXjAklq7gyW99gPo8IVuTNNoMGRCBg3SZtjdu4O4OD+2Hzo2yZvo+CSGTM/Vq1cZP348gwYNiq7r0aMHu3fvtjv7cQXOKp/66DTaJ20rReQUOpdPfdeIZTAYAPDygrlz4Y8/YPhw149fGW351gntMGHI1EybNo2nn36aMmVi3Db9/f0ZOHAgn376qUvncvbMxwe4Ec+9G5b7BoPBlWTPDitXQkgIlCqlT3xdSR/ge/S+xaCEmxoyNu3atSN37txx6nv37k2DBg1cOpezK5/9wGvK6nVkQSmlgJct9w0Gg6spVAhWrYL+/WHTJteO7YWOTT8OHQHbkGkJDg4mZ86cceq9vb156KGHXDqXsyuf4WgXtT+UUguB82iDg3ZAWUxsN4PBfVSooB1P27fXbuiuMDmyUgyYAHQE9gLZXTe0weAIZ51M1wDN0Vtsg9BptAcDN4HmIrLW5RIaDIYYnnwSRo6EZs3g0iXXjt0OHTyrv2uHNRgc4bSfj0UBrVFKZUObXF8Rkdsul8xgyIzkyZO0kDmRkVC6NNR3sY3PeKAKMAsdOngJ2hLucaAh0B0d295gSCFOKx8rFoVjlI7B4EqSGiRURCdciYiAqChtFecKAoGngJ7AJPRW3JvAj8BK9LbcU+jVkXFMNaQAE17HYEiPKAUzZuhsX++847pxR6PTL/RHr3qsgbW9gTboeHCnLO0MhhSQ6MpHKRUFJNW7TUQk2aspg8HgBL6+sGwZ1Kqlt+B69kzZeOeB1WgFEwnUhuzfZyd7dhvrAy9gLHoLLgyzBWdINklRFMNJuvIxGAypSd688P338NhjUKKEDsmTXKYDr6AVjA8wFyrVrcT3G763b+eFdqz4Em1ulIHZv3+/p0XIsCS67SYiQ0VkWFJLaghtMBhsKFNGR8Hu1Al++y3546wDnra5Lg8MB99OvnArVtsWpIuUDE2aNEEpxeDBMVryxx9/pGPHjpQuXRp/f39Kly5N7969uXjxYoJjjRw5EqUUjz32mF39/fv3GThwIAUKFKB48eJ8/vnncfp+8803FClShOvXr7vkc2UEnDrzUUrtUEq9oJTydZdABoMhGdSpA+PHQ/PmcO5c8sZQ6LMdW15Ch+DpBETZ1MdulwaZP38+v/4aN2f45MmT+e+//xg8eDBr1qzhnXfeYcWKFYSEhHDz5k2HYx0/fpwPPviAAgUKxLk3a9YspkyZwrhx4xg0aBBvvvkmm2wcgW/evEm/fv0YPXo0OXLkcNnnS+84ez4TjjbC/FwpNQuYKiKHXS+WwWBwmvbt4fhxrYC2bIGAAOf6CxCBvWJRwFR0Cob3gA8s9REpF9edXLlyhb59+zJmzBief/55u3sTJ04kf/780dehoaE88MADhIaG8s0339C1a9c44/Xu3ZuwsDCOHDlCZKR9ELzVq1fz/PPP89xzzwGwfPlyVq9eTT2LyfzQoUMpX748HTp0cPGnTN8462RaD3gIrYA6AYeUUpuUUu2VUungt5DBkMF5+2145BGdByjSyUihDdHm1LHxRVu+zQPmWupWWNqnUQYOHEjFihUdfuHbKh4rjz76KABnz56Nc2/evHns3buXjz76yOFc4eHh+Pv7R19ny5aNu5ZUnwcPHmTy5Ml88cUXyfocGRmnTa1F5LCI9AOKAi+ijTHnAX8rpUYqpUq5VkSDwZBklIJJk+DePejTx7k0DN3RMUuiHNzLj1ZM/YBtwERL+zTItm3bmD17tlNf+Js3bwbgwQcftKu3rqA+/vhj8uTJ47BvzZo1+fbbbzl27Bh79uxh7dq1hISEAPDyyy/Tp08fypUrl8xPk3FJtp+PiNwTkTnoRHJb0X+ebwFHlVKLlFKFXCSjwWBwBm9vWLxYByB1cPgdL4XRDqRv4FgBVUCn324M1ERHdUxjhIeH06tXL958880kf+HfuHGDPn368OCDD9KqVSu7ewMGDOCBBx7gxRdfjLf/a6+9RnBwMGXLlqV69eo0b96c9u3bM3v2bP7++2+73DiGGJKlfJRS/kqprkqpn9FxcAug/2SLAL2B2ug0VQaDwRPkzKmjYH/6KSxdmvR+/YEg9Jbat8Sc7URYrj9Dn/+sQKeVTGN8/PHH3LlzJ8lf+JGRkXTo0IGzZ8+yYMECsmaNOQbfunUrs2fPZtKkSejA/Y4JDAxk06ZNnDx5knPnzjF//nyuX7/OgAEDGD9+PL6+vgwaNIhixYpRtGhRBg0aRFSUI+2euXA2jXYldNb3MHTc2+XAQBHZaNNsmlLqAjozvMFg8BRBQbB8OTz1FBQtCjVqJN5HocPphKH9fsbZ3GuIPvMphP6J+Tz6GyALaYLTp08zYsQIpk+fzr1797h37170vXv37nH16lUCAwPJkkULHBUVRefOnVm/fj2rVq2icuXKduP16tWLbt26UaxYMa5evQpoZXX//n2uXr2Kv78/vr4xhr9BQUHR7wcNGkTt2rVp1qwZ06ZNY+7cuWzduhXQBg6lSpWiW7du7noU6QMRSXJBL8b/BoYAhRNo9yCw0ZmxU6NUq1ZNDIZMx7JlIoULi5w44boxw0XkCRHp57ohU8rGjRsFbbMXb9m3b190+x49ekiWLFlk6dKlDsdLbKwxY8Y47Ld7924JDAyUU6dOiYhImzZtpF+/mAfVr18/adu2ras+ttuYPHmy9OzZM0VjALslnu9jZ02t2wLLReR+IgrtD0xKbYMhbdCyJZw8CU2bwo4dkCtXysf0Ru9thKB/aqYB44OqVauycePGOPX169enY8eOdOvWLTo9dP/+/Zk+fTqzZs2Kc85jxdFYffr04f79+4wfP94u1bSVqKgoXn75ZQYPHkyJEiWi62/divHSvXnzpvVHeqbGKeUjIkkMuWswGNIUb7wBf/0FzzwDq1eDjwsy3udBp5asC5TG4z83c+XKFe1bE5ugoKDoe6NGjeKzzz6ja9eulC1bll27dkW3y58/P6VLlwZwOFauXLmIjIyMd55p06Zx69Yt+vbtG13XoEEDBg4cSN26dQFtuj16tInMmqwgoEqpKkA5wC/2PRGZnVKhDAaDGxgzBlq3hpdegi+/1GbZKeUBYD7wHNoEu2zKh3Q3q1evBuCrr77iq6++srvXuXNnZs6cmaxx//33X/73v/+xdOlSvL1j3B579uzJX3/9Rb9+/RARXnrpJXPeAyhnln9KqVzAKvRiG2IyekQPIiJp5PgxLtWrV5fdu3d7WgyDwXPcvAmhodCmDbjSBHgKMAbYiU4xaUj3TJkyhb179zJlypRkj6GU2iMi1R3dc9bU+kMgLzqvoQJaA0+gzaqPA0kwpzEYDB4jIAC++w6mToX58103bi+gCfAsaT70jiFt4KzyaYxWQNZN0r9FZJOIdALWo319DAZDWqZwYa2A3ngDLOa/LmE02hDhDUwSFkOiOKt8CgPHLdZud9FJd60sAZq5SjCDweBGKlWCuXOhXTs4etQ1Y2YBFgBb0Om3MwBVq1Y1OX3chLPK5wKQy/L+FFDL5l5cu0ODHT169EApZWcJk1kYOnQoGzZs8LQYBlsaNYIPPoBmzeDff10zZg50DLgPgTWuGdKQMXFW+WwjxthgDjBEKTVFKfUF8AnwgyuFy0jcuXOHb775BtCmlrHDsmd0hg0bZpRPWqR7d2jbVvsCWSIxp5iSaB+gTsDvrhnSkPFwVvkMI0bBfIKOgdsM6ICO9vSa60TLWCxbtozr16/TtGlTLl68yJo1rvtZKCKEh4e7bLz0gm34FEMKGDECiheHF18EV8Ucewz4FJ0Z1UWLKkPGwlnlcw0dSBQRiRCR/iJSTETyiMjzIvKf60XMGMyaNYvcuXMzc+ZM/P39mTVrlsN2S5YsISQkhGzZspErVy7atWvH6dOn7doEBwfTsWNHvvrqK8qXL4+Pjw+rVq0CYM2aNdSqVQt/f39y5sxJq1atOHLkSJx5li5dSp06dQgICCBHjhzUqFGDFStWRN+fMGECtWrVIk+ePOTKlYuQkJDoOaycPHkSpRQTJ06kX79+FChQgGzZstG8eXNOnjwZ3c4alHHEiBEopVBKMXToUEA78jly2AsODraLJDxz5kyUUmzZsoV27dqRK1cuatasCeh4Wx999BHly5fH19eXIkWK0L9//+icKoZE8PKCmTPhzBmwSTedYjqhrd/aAOZ3giEWiTqZKqWyAO+ibVhyAPeVUiuBbiJy1b3iZQzOnTvH+vXr6dGjB/nz56dVq1YsWbKEK1eukDt3jFPE5MmT6d27N126dOG9997jxo0bDB06lNDQUH777TcCA2PsOzZu3Mj+/fsZMmQIBQoUIDg4mDVr1tCsWTOeeOIJFi5cyM2bN3nvvfd47LHH2L9/P0WLFgVg/PjxvP7667Rq1YpZs2YREBDA3r177RTGyZMn6d69O8HBwURGRrJy5UqaN2/O6tWradKkid3n++ijj6hatSozZszg4sWL/O9//6NRo0YcOnQIb29vdu7cSa1atXjxxRfp1asXAMWKFUvWswwLC6NDhw4sXrw4euuyY8eOrFy5koEDB1K7dm3++OMP3n33XU6ePMm3336brHkyHX5+kDu3dkT99lttEZcYefLAkkSCnoxAB+XqDXxJjGegwRBf0DdrAV5BBxTdAIxC7+ZGAjMS65vWiqcCi44aNUoA2bFjh4iIrFmzRgCZNGlSdJsbN25Ijhw5pEuXLnZ9jx8/Lt7e3nZBDIOCgsTf31/Onz9v17ZatWpSpkwZiYiIsOufNWtW6du3r4iIXLt2TQICAqR169ZJlv/+/fsSEREhDRs2lBYtWkTXnzhxQgB58MEH5f79+9H127ZtE0CmT58eXQfIoEGD4owdGhoqoaGhceqDgoKkc+fO0dczZswQQPr06WPXbsuWLQLIrFmz7Ornzp0bJ5CkIRFCQ0WOHBEpWFBk7dqktU8KN0XkYRH5OPmieYoqVapk2r8hdwcWTcq2Ww9gmugYtgNFpJ1FIXVUSrkgQFTGZ9asWZQtW5ZatbRxYIMGDShSpIjd1tvOnTu5fv06YWFhREZGRpfixYtTvnx5tmzZYjdmSEgIhQrFZPO6desWe/fupX379nY5SUqWLEmdOnWiMzXu2LGDmzdv0rNnzwRl3rNnD82bN6dgwYJkzZoVb29v1q1b53ALr23btnh5xfwp1alTh2LFirFz504nnlLSaN26td31mjVr8PHxoW3btnbPrVGjRgBxnpshER54ABYtgrAwOHjQNWNmR58Ij7W8Ggwk7cynFHFz8yxEW/UHxW1usGX37t38/vvvtGnThqtXr3L16lVu3LhBmzZt2LVrF0ctPhYXL14EtGLy9va2KwcOHOC//+yP0wrH2ha5cuUKIhKnHqBQoUJcvnwZIHqchLa9zpw5w5NPPsnly5cZP348O3bs4JdffqFJkyYOz1EKFizosO7s2bMJPZpkEfvzXbx4kfDwcLJnz273zAoUKAAQ57kZkkDdujoDarNmcP68a8YsBiwFugH7XTOkIX2TlMCiAcTNWXjD8hqIIUGsq5tRo0YxatSoOPdnz57NBx98QN68eQF9sF6hQoU47WzPe4A4mRVz586NUooLFy7E6XvhwoXo/PP58uUD4OzZs1SsWNGhzGvWrOHatWt88803dkrq9u3bDtv/888/DuuqVq3qsL0tfn5+XL8eNyWmVVnGJvbnzps3L35+ftGJumJTpEiRRGUwOOD55+H4cXj6adi8GbJnT/mYj6LtY1sCP5Em03AbUo+kWrsVVUqVshb0aihOveWewUJ4eDjz58+nZs2abNy4MU6pWrUqc+bMQUSoXbs2gYGBHDt2jOrVq8cpieWjz549O9WqVWPRokXcvx+TbunUqVPs2LEj2qKsdu3aBAQEMHXq1HjHsioZ28i8R48eZfv27Q7bL1682C4t8Pbt2/n777+jtxkBfHx8uHPnTpy+QUFBHD161M5UfMuWLdy4cSNOW0dYV2PXrl1z+NyM8kkBgwbpSAjPPguuMmt/Fr36aQXE/XMwZCbiOwyyFrSxwX0HxWF9YuM5GL84sBhtxn0dHaanRBL6VQemAoeB28BpdIDTkvH1SW2DgyVLlgggM2fOdHh/0qRJAsiGDRtERB/wZcmSRXr16iXLli2TjRs3yty5c6VHjx7y9ddfR/cLCgqSsLCwOOOtXr1avLy8pEmTJrJixQqZN2+elC1bVvLlyydnz56Nbjd+/HgBpE2bNvLtt9/K2rVr5eOPP5Zx48aJiMjBgwcla9as0qhRI/nhhx9k5syZEhQUJCVLlpSgoKDocawGB8WKFZOnn35avvvuO5kxY4YUKlRIypYtK+Hh4dFtq1atKuXKlZO1a9fKL7/8Ei3Phg0bBJCwsDBZt26dTJ06VcqXLy85c+Z0aHDw559/xvncHTp0kFy5csnw4cNlzZo1snbtWpk6daq0atVKjhw5koR/KYOIODYguHdP5JlnRJ56SuT27cTbJ4UoEXlORDpY3qdhjMGB+wwOkqIcOjtTEhsv1tjZgD+Bg+jfQi2BA8BfQPZE+n4KbAdeBkLRGeX/AP4Dijvqk9rKp2XLlhIYGCi3bt1yeP/q1avi7+9v9yW7atUqqVevngQGBoq/v7+UKVNGunTpIocOHYpuE5/yEdEKKCQkRPz8/CRHjhzSokULOXz4cJx2ixYtkho1aoifn58EBgZKjRo1ZOXKldH3Fy5cKOXKlRNfX1956KGHZP78+dK5c2eHyueLL76Qvn37Sr58+cTf31+aNm0qx48ft5tv27Zt8sgjj4ivr68AMmTIkOh7kydPljJlyoifn5/UqlVLdu/eHa+1myPlc//+ffn888+lcuXK4uvrKzly5JDKlSvLgAED5OrVqw6fU1I5duyY9OvXTy5evJiicdIF8SmTiAiRDh1EGjQQuXkz8fZJ4baI1BSRYckfIjUwyseDysedBe07dB8oY1NXEm3K3S+Rvvkd1AVZVmTDHfXxlKl1RsWqfKZNm+ZpUVzOsWPHpEuXLpI3b15577337MzXMywJKZPISJHOnUUef1zk+vXE2yeF8yJSQkQWpmwYd2KUj2dNrd1JC2CXiByzVojICfSKpmVCHUXkkoO6U8AloKiL5TRkEv766y+6du1KzZo1KVGiBH/++SfDhg2zM1/PlGTJAl99BeXKQePGcO1ayscshDa9fgX4OeXDGdIXnv4fVQFY7qD+ENDO2cGUUg8CBdDbbwZDkrl58yZPPvkkBw8eJCQkhJdffhlfX18mT57sMZkCAwN55ZVX4lj4eQwvL5gyBV5/HRo0AF/flI9ZBR35oA06S1jyAl8Y0iGeVj55gCsO6i/jZDJepVRWYDJ65fOlozaXLl2ievWYjK49e/ZM1NnSED/BwcHW7c50T2RkJDdv3kQphYhw+/ZtIiI8m5IzTT5bpWDcOBgwACZPhkuXIH/+lI3ZAjhied2Kdko1pFumTp1qa02bL752ypN/4EqpcOAzEXk7Vv0HwNsikmTlqJSajDbibCYiax21qV69uuzevTslIhsyOL/99hvDhw9n+/btDBgwgJdeeols2bJ5WqzUoU0biMe/yiEXLkDWrLB+PRRKodOOoP/3XgG+Bf4BpqHzIyvL/YZAd3RKy1SiatWqzJw5M0k+a4a4KKX2iEh1hzfjOwxKjYL+E5vioH4icMmJcUaiDQ1eSKhdejE4GDJkiOh/Gue4f/++vPHGG1KoUCFRSknLli1dL1wiWI0QZsyYEV3XuXNnKVq0aKJ9rRZtJ06ccJ+ASeTXX3+VZ555RgoVKiSjR4+2i11nsGH4cJFy5UT+/jvlY90TkcdEpL7oYF7fiojVWj/ccv2EiHwiqWainZkNDlwBadjg4BD63Cc2D5HENFRKqUHAQOB1EZnjQtk8Rvfu3ZMVF23x4sWMHTuWAQMGsH37dj7++GM3SJcwhQsXZufOnTRrlr4zqleuXJnFixfzww8/8Mcff/CvqzJ9ZjTefRe6doXQUDh1KmVj+QANgD3AC+hzIKufs7fleh06h/LolE1l8DyePvNZAXyqlColIscBlFLBQB3g7YQ6Wtq+DnwADBKRDJI1XsddS07KgT/+0HYWffr0sQv0mZr4+voSEhKSeMN0QuXKlZk2bZqnxUjbvPWWNj4IDYUff4TSpZM3znlgC9rW9QmgDDopnS1e6AClDYEwUnULzuBaPL3ymQacBJYrpVoqpVqgrd/OAFOsjZRSQUqpSKXUezZ1zwGfozPFb1BKhdiUh1LzQ7iaoUOH2lk4KaUYPHgw48aNo2TJkgQGBhIaGsqhQ4ei2wQHB0cnaMuSJQtKKWbOnAnocDkDBw6kZMmS+Pj4ULJkSUaMGGEXEufu3bv07duXihUrEhAQQKFChXj66ac5fPiwnWwXLlygc+fOFClSBF9fXwoXLkzz5s2jA6NaE8xZ57Zlx44dPProo/j5+REcHMz48eOT9DymTp1KlSpV8PPzI1++fHTr1i3e2G8GD/HGG/D221CvHjiIfJ4kpqPNrisCc9D2ricctPNCu5Y7NCsypBc8uvIRkVtKqSeAMeg/NwX8CPQRkZs2TRU6iratsmxiqW9iKbZsBuq5SWyPMHfuXMqVK8fYsWMJDw9nwIABtGzZksOHD5M1a1aWLl3KuHHjmDlzZvSWXenSpYmMjKRx48b8/vvvvPvuu1SqVIldu3bx/vvvc/nyZUaP1vsX9+7d48aNGwwePJjChQtz+fJlJk6cSK1atfjjjz+i0ze88MILnDp1ik8++YTixYvzzz//8OOPP8YbdNTK9evXad++PQMHDqRMmTIsWLCA119/ncDAQLuMpbF5++23GT16NK+//jqffPIJZ8+eZfDgwRw8eJAdO3aQJUsW1zxgQ8p56SW9AnriCVi7FhwEyE2QdcTsdzQG/odOw70DncbSlhbAOMCFiVcNqUx8h0EZsaRXgwNAypQpYxcrbdGiRQLI9u3bo+sGDRoUx1Bh9uzZAsjmzZvt6j/44APx9vaWf/75x6EMkZGRcuvWLQkICJDPPvssuj579uwyduzYeGWPz+AAkPnz59u1bdCggZQoUUKiovTpcWyDgxMnToiXl5cMG2Yfg8WarG7p0qXxymHwIHPnihQqJLJ/v3P9Ho91HSUir4hIHRG5koT2biCzGRxs3bpVypcvHx3ZQ0SH/HrggQekaNGiMnLkSKfGIw0bHBiSSMOGDe2iTFeqVAmA06dPJ9hvzZo1BAUFUbt27TjJ1iIiIti1a1d022+++YaaNWuSK1cusmbNSvbs2bl586ZdArlHH32UTz75hLFjx3LgwIEk+6JkyZKFZ555xq7uueee4/Tp0/Hm/Vm3bh1RUVFxEuzVrFmTwMBAkygurRIWpn2BGjUCZ1wbBLB1rVLo1c0jQH3gos09z7pgZUju3btH7969mTFjBtu2bWPx4sW88847dOnShU8//ZQdO3bw448/sn79epfMZ5RPOsGaj8eKr8W73FFyN1suXrzIqVOn4iSoq1GjBhCTbG3lypW0b9+eBx98kHnz5vHTTz/xyy+/kD9/frs5Fi5cSIsWLfj444+pXLkyRYsWZfjw4XbnR47InTu3nfKEmCR08Skf6zlSmTJl4sh/48YNkyguLdOuHUydCk2bgs0PnARpCKyMVWc1MGgB1EWfBoM2VWroGlENmgMHDlC7dm1CQkIoX74833zzDZ988gnDhg3j6aefpkSJErz77rv88MMPLpnP09ZuBjeTN29eSpYsyTfffOPwfnBwMAALFiygTJkydoYCERERcQ72CxQowBdffMEXX3zBkSNHmDVrFkOGDCF//vz07t07XjmuXLlCRESEnQKyJqErWtRxKD5rgr21a9eSO3fcgBfW+4Y0SsuW4OMDLVrAt9/qDKkJ0R3oiI5vb/uzWAHDgJxoBbQG7Qn4tRtkzsRERUXZxTAsUqQI2bJls0v2mCVLlkR/aCYVo3wyOE2aNOHbb78lICCA8uXLx9vu9u3bcYJnzpkzxy4xXWzKlSvHhx9+yOTJkzl48GCCcty/f59vv/2W5557LrpuwYIFlChRIl7l07BhQ7y8vDh9+jQNG5qfuemSp56C+fN19IQFC+DJJ+NvWxh4Ch3rfixx92X6oXMnVwN6YTKhupiKFSuyZcsWfv31VwoVKkSXLl0ICwtj5MiRVK5cmdq1azNq1CiXhSQzyieDExYWxowZM3jyySfp378/VapUITw8nL/++osVK1awbNkysmXLRpMmTVi2bBl9+/alefPm7N69m/Hjx5MrV67osa5du0aDBg0ICwujfPnyeHt7s3z5cq5cuUKjRo0SlCMwMJC33nqLf//9l7JlyzJ//nzWr1/PzJkz4w2cWbp0aQYOHMirr77KkSNHCA0Nxc/PjzNnzrBu3Tq6d+9O/fr1Xfm4DO7gySf1yqdtW5g9G5rENk61oT/agbQh2py6BdrBNAK91bYA7Ww6F2gL1Hav6JmJbNmyMXbsWFq3bs3Zs2fp1KkTEyZMoGHDhnTr1o2rV6/y0ksvucyB3CifDI63tzc//PADI0eOZOrUqZw4cYLs2bNTunRpmjVrho+PDwA9evTgzJkzfPXVV0yZMoVHH32UlStX0rp16+ix/Pz8eOSRR5g2bRqnTp3Cy8uLcuXK8fXXX9OyZYIZMMiRIwcLFizgjTfe4MCBAxQsWJCxY8fSuXPnBPt9+OGHPPjgg9FbfUopihcvzpNPPknZsmVT/oAMqcPjj8Py5Xorbvp0vRXnCAW8iXYgnY42OLDSEK10CqO33lqht97MothlPPHEExw/fhwRif5R2KZNG9q0aUNUVJRLndc9Glg0tTGBRQ0GD7N7NzRrBl98oVdCKWEb8AwwCb0acgMmsGjKSCiwqFn5GAyG1KN6dfjhB30WFB4Ozz+f/LEeQ6+AmgHXgRddIqEhlTDKx2AwpC5Vq8K6dTojang4JBDhIlEeBjYCjYBraGMFQ7rA+PkYDIbUp2JF2LBBR8WOSTyWPMqhA5J+AQxHO6saUsyUKVPo1auX28Y3Kx+DweAZypWDTZu0Ndzduzo9d3IJQmdBbQRcRVvMpZHs4wbHmJWPwWDwHKVLw+bNMHYsfPJJysYqCGwCdqGzokamWDqDGzErH4PB4FmCgrQCevBBmDhRXyeFPHlgyRL7utzo6NitgefQpti+LpXW4CKM8jEYDJ6nWDGoVAmuXdNheIYPh3icj6OpV89xfXZ0jLjn0U6qSyx1hjSF2XYzGAxpAx8ffQa0YgUMHAgp8UH0BRYCRYg5BzKkKYzyMRgMaYf8+bUV3IYN0KdPyhRQVnS20+rETclg8DhG+RgMhrRF3rywfj38/DP07g0piaLsBXwOtERHxE44/ZUhFTHKx2AwpD1y5dKpuH//HZ59FmzC+juNAoYCvdEK6KgrBDSkFKN8DAZD2iQwUCugvHl1WJ7ffkvZeH3QSqgesD+FshlSjFE+BoMh7eLnB1Om6EgITz4JNskOk0UXYDzQGNiecvEMyccoH4PBkPZ54QVtCTdqFHTrBnfuJH+sZ4DZaF+gta4Rz+A8RvkYDIb0QYUK8MsvcPs21KqVMgXUGFgKvAB86yL5DE5hlI/BYEg/BATAvHnQsyfs2xc3woEz1AF+AF4DZrhIPkOSMcrHYDCkL5SCl1/WERH694e+fXVqhuRQFZ2SYSgw1qb+PDpC9jF0nLjHgfct9QaXYMLrGAyGtEGePPGHzHFEUJDOC9Spk+63cCEUL+78vOXQEbEbApeBAHSSuleA0mhH1QrokD0dgaeA/pio2SnEKB+DwZA2SO4W2ooV8PHH8OijMGuWTlLnLCXQOYEeBgoAu9HfjsMt973RqbpboRPWjQbeTJ64Bo3ZdjMYDOkbLy94+21YsAC6doWhQ+H+fefHiQLKAtmAHjhOyeCF3p5bjdmCSyFG+RgMhoxBvXqwe7c2yW7SBC5dcq7/dPSqZh1wDmgPr730GqVKlbJv5wW8jN6OMyQbo3wMBkPGoXBhHRfu0UfhkUdguxOepOuAp9HpF1YACrot6UYOyRG3bQtLe0OyMcrHYDBkLLJmhQ8/hEmToE0b+OyzpEXHVuizHdApGRYADwBV0OdBtnhjSCFG+RgMhoxJ8+bw008wfz488wxcvZpwewEibK6zAhMs5TlgIHDPci8CQwoxysdgMGRcgoNh2zYoUkQHJ92/P/62DdHm1LFpjg5EehioCRxCb8s1dLGsmQyjfAwGQ8bG1xcmTIAPPoCGDWH6dMfbcN2BL9BWb7EpACwDXgVCgXeArm6TOFNglI/BYMgcPPccbNkCY8ZAly46RpwthdEOpG/gWAEptMJ5CggHOgF/u1XiDI1RPgaDIfPw4IM6Q2pUFNSsCUeO2N/vDwSht9S+JeZsJ8Jy3RBtgPAnOi/QI8DCVJE8w2GUj8FgyFxkz64jIbz+Ojz2GHzzTcw9hY5cMBf4HWiA3mZrAPxhqX8Tbe02GFgFvIcOu3M19T5CRsCE1zEYDJkPpaBHD6hWDdq1g61b4cyZuBZxtvHb1luKlTx5dEigvcAA9IpoJlDfrZJnGIzyMRgMmZdHHoE9e+DFF2HDBjhwQAcsTQrWIKjZgYnA90CYpXyA9hUyxIvZdjMYDJmbXLlg6VLInx9q1IDVq5M3TlPgV3QahhrAAZdJmCExysdgMBiU0ukYFi/W23GDBycvOGl+YAnQB3gC+AzHlnMGo3wMBoMhmrp1Ye9e2LULGjWCf/5xfgwFdAF+QiuiBsAZl0qZITDKx2AwGGwpUAB++AHq1NEGCYsXJy02XGxKAZvRyqcaMM+1YqZ3jPIxGAyG2GTJAsOHw9df68gIoaHaMMHpcYD/ofP/vA90AK64VNJ0i1E+BoPBEB9WpdOpkw5U2qULnDvn/DjVgD1APrRJ9gbXipkeMcrHYDAYEiJLFujeXUdDKFgQKlWCESN0lARnyAaMB6YBLwD9gLsulzbdYJSPwWAwJIUcOWDkSPjlFx0d++efYeFC58+DGgO/AaeBR9Hm2bacB4YDj6OjKzyO3rLLYGm7jfIxGAwGZyhVChYt0nHiRo3SIXp+/tm5MfICi9ChehoAnwCRwKfoUD0VgR/RBgs/AhUs9Z+i8w5lAIzyMRgMhuSQM6deBXXvDq1a6XOhs2eT3l8BnYGf0fmBygEH0em52xCTLdXbcr0OOAWMdtUH8CwmvI7BYDCAjtVmDZmT1PZZsmgjhLZt9ZZclSo6YOmbb0K2bEkbpyQwHx0l+zu0SXYY9nHlQC8VxqIja4ehU0CkY4zyMRgMBtBBQpNLYKA2QujRA95+G8qXh48+gg4dwCsJG0wzgI/RiigMnVF1EpAnVjsv4GXgS3RU7XSM2XYzGAwGVxEcDAsWwLx58PnnULu2jpaQGOuAp4GH0SbZhYDKwA8O2rawtE/nGOVjMBgMruaxx+Cnn+Dll/WWXFiYTtkQH4qYMx5/9PbaDKA32iBhm01bbzIERvkYDAaDO/Dy0kYIhw9D6dJQtSoMGQK3bsVtK8RkTbXSEDgCPI9O2d0A2OqgXTrF48pHKVVcKbVYKXVNKXVdKbVEKVUiiX39lFKfKKXOK6XuKKV2KqUed7fMBoPBkGQCAnSonn374NgxKFcO5syxd1JtiD7niY030JUYJdQZnbq7rPvFdjceNThQSmVDB5q4h36sgk7DtFEpVVlEHPxEsONLoBk6j+Bx4BXgB6VULRHZ7zbBDQaDwVlKlNCx4nbuhL59Yfz4mHOhnd3gs5dh7JugEnDkKZoV9k6Hv7LBiQIwBO2Emh4REY8V4A3gPlDGpq4k2t2qXyJ9q6CVVRebuqzo3wgrHPWpVq2apIQpU6akqH9mxjy7lGGeX/JJk8/u/n2RuXNFihUTad9epGZNkU9E5FURuR9fHxF5RXS7uk+IfCUipUSkvohsdr2IkydPlrp166ZoDGC3xPMd7ulttxbALhE5Zq0QkRPAdqBlEvpGAAtt+kYCC4DGSimXJ7GdOnWqq4fMNJhnlzLM80s+afLZeXlpI4TDh3WkhD174MpgKHRPb8F9S8zZToTluiEQDPQHvO7rnEGH0edBXdDJ6za7TsQmTZpw6dIl1w0YC08rnwpon97YHAIeSkLfEyJy20FfH6BMysUzGAwGN5I9uzZCqF4dzpyGiaWg+QI4FKUNDELRr38Ac9HheGydT72BF4lRQt2A+rhECQUFBZE9e/aUDxQPSpKTJMlVkysVDnwmIm/Hqv8AeFtE4j2TUkqtBXKISEis+gZoK/jHRWRrrHs3sFe4l4B/nRA5n5PtDTGYZ5cyzPNLPmn+2ZWHcof1kYFb2qeA5Dy7fOiE4gBRIhLoqFGminAQ30MwGAwGQ+ri6W23K0BuB/V5SDzfX0J9AS6nQC6DwWAwuBFPK59D6LOb2DwE/J6EviUt5tqx+4YDx+J2MRgMBkNawNPKZwUQopQqZa1QSgUDdSz3EmIl+ritnU3frEB7YK2I3HOFgClxgs3sKKXqKaXEQbnqadnSEkqpYkqp8RYn6duWZxTsoJ1xqo6FE8/O0d+hKKWqpr7UaQOlVFul1LdKqVOWv6cjSqmPlFKBsdrlVkpNV0r9q5S6pZRar5SqlNL5Pa18pgEngeVKqZZKqRbAcuAMMMXaSCkVpJSKVEq9Z60TkX1oM+vPlVLdlVJPos2sS6Jdr1KMjRNsebQT7Ato3+KNSin3mYFkPF4HatmUBp4VJ81RBngWvZW8NYF2XwI9gPeA5ujclj9k5i9Qkv7sAGZi/3dYCzjqTuHSOG+i/Sz/BzRBx9HuDaxTSnkBKKUU+od+E+A14Bn0j/6NSqliKZo9Pgeg1CpACbQV+3XgBrAMCI7VJhjtUDo0Vr0/8BlwAZ0N/SegngtlS7YTrCkCOkOJAA08LUtaLoCXzfvulmcWHKuN007VmaEk5dlZ7gnwgaflTUsFyO+grpPlWT1huW5pua5v0yYn+kx9XErm97i1m4icRmvThNqcJG5qJUTkDtDPUtyBQydYpZTVCfYzN81ryESISFTirRw7VSulFgBvK6V8xUVbzemJJD47gwNExJEH6S+W16KW1xbAORHZaNPvmlJqJfo78PXkzu/pbbe0TkqcYA0xfK2Uuq+U+k8pNc+cmSUL41Sdcnorpe5ZzoY2KKXqelqgNEio5fUPy2tC34EllFIByZ3I4yufNE58Jt+XcWzmbbDnGjrj/Gb0turD6P3lnUqph0XkoieFS2ck9LdovW+In7noJNXngCB0MOINSqmGIrLJk4KlFZRSRYHhwHoR2W2pzoM+l4+N9e8uN3AzOfMZ5WNwG6KNQvbZVG1WSm0BfkYv19N5ImBDekFEXrC53KqUWo7+Rf8B8JhnpEo7WFYwy9Hn2V1SY06z7ZYwKXGCNThARPaiLYwe9bQs6QzjVO1CROQGsArzd4hSyh9t0VYKaCwif9vcTuzvLtnfg0b5JExKnGANCeO5oILpE+NU7R4y9d+hUsobWAxUB5qKyIFYTRL6DjwtIsnacgOjfBIjJU6wBgcopaoD5dBbb4akkypO1ZkFpVQOtK9Upv07tPjyfI1OxtBKRHY5aLYCKKqUCrXplwN4mhR+B5ozn4SZBryKdoIdjP6V9D6xnGANjlFKfQ2cAPYCV9EGB+8AZ4FxnpMs7aGUamt5W83y+pRS6hJwSXSqsH1KKatTtTf6ufZG+52Fpb7EaYfEnp1S6k30D56NxBgcvAkUInM/uy/QP2ZGALeUUrYZAv62bL+tAHYCc5VSA9DbbO+gXV8+TtHsnnZ0SuuFJDjBmhLvs3sH+A1t9RaBVtpTgcKeli2tFfQPG0dlk00btzpVp9eS2LND/0rfjk4NEAH8Z/lSreFp2T383E4m8OyG2rTLA3yFPle8DfwIVEnp/B7N52MwGAyGzIk58zEYDAZDqmOUj8FgMBhSHaN8DAaDwZDqGOVjMBgMhlTHKB+DwWAwpDpG+RgMBoMh1THKx2AwGAypjlE+BoPBYEh1jPIxpEuUUtOUUqKUGuPCMYcqpdzidZ2UsZVSy5RSV5RSvvHcD1RK3VJKzXSnjJaYcW5FKVVYKRWllMr06QwyK0b5GNIdlhDwz1oun3fhl+V0oJaLxkoOs4Bc6ICXjmgLZLO0S++0BC4BOzwtiMEzGOVjSI+0AnIA3wMFgCaJdYhvNWF7T0T+FseRfVOLVei4Y53iud8JOA1sSi2B3EgrYKWIRHlaEINnMMrHkB7pjI6u+yJwx3Idjc32UUWl1A9KqZvAN0m9Z3nfztKucuzJlVLfK6V+tbkuo5Sao5Q6oZS6o5Q6rpSapJRyKtW6iIQD89FRmfPGmrMEEArMkVQMyKiUaqKUuqmUmqCU8rJ5fuUtz++WUuq0UqqLpf0LSqnDlj4blVKlHYyZA6iPDtKLUuoBpdRSpdRFpdRdy3iLUmP7z+A5jPIxpCuUUkWABsBCEbmE/gJ7Op4v+uXAZqAFEPtsKKF7oPPnXAM6xpq/INAImG1TXQQdsbsP0BgYDjyJXpk5yyx03p7nYtV3RIexnx2nh5tQSnVCR38eKSKvxlqlLEKv1FoBe4CvlFIfotM8vI1OxVwOmOdg6KboBHjrLdergKKWvo0t/e9hvp8yNp4O622KKc4U4C10yPdaluvGluuXbNoMtdS94aB/ovdsrqcBfwNeNnV90Hnu400Lgc6T9ZhlnocdjZ3IZzwE/BSr7g9gp5ufrfXZZLU85wigezxtOtnU5bY8k/+AHDb1r1vaBsUaYwGw2PI+n6VNC0//bZmSusX8sjCkNzoDf4rITsv1enSCsM4O2i5NYJyE7lmZjf5F/oRN3QvAjyJy3lqhlPJRSv3Pst10B/2lvdVyu1wS5onNLKCGUuoBy/g1gPKknqHBGGAY0FZEpsfTZrX1jYhcAS4Cu0Tkuk2bw5bX4tYKpZQP8BSWLTe0wjoOjFRK9VBKlXXJJzCkeYzyMaQbLCm4HwKWKKVyKaVyAYHAEnS68wdidTlP/CR0z8o2dMKtFyzzPwg8Qtytr4/QK4K5QDOgBtDGcs8vCfPEZi4QRYzhQSf0NtTCZIyVHDoAB4nZFnPElVjX4fHUgf0zeAJtsfcdWJaC0BDYjX6ORy1nZr2TJ7ohvWCUjyE9YV3dDER/0VnLq5b62FZiCR3MJ3pob/linAu0UUplQyuhm8RdNT0HzBaRD0Rkg4j8gk4bnixE5BywDuhoWSm0R1uG2X25K6UKKKUWKKWmKKWyKaVethhDNEru3BaeRGfwXa2UCkjhWLFpBWwWkavWChE5LiKdgPzoVOsbgIlKqadcPLchDWGUjyFdYPkS7oBOHV3fQdkPvKCUUi6eeg4QgF7JhAFLROR2rDbZ0FtttnRJ4byzgCD0aiAfjrfcBgLvARvRK4ka6GcUZnleyeUQUA8oiwsVkOXfpgUxW252iGY/0M9SVdEV8xrSJsaU0ZBeaAbkBfqLyKbYN5VSU4BJ6C9NlyEiR5VSPwEj0ec/jqzN1gCdlVIHgGNoRVU7hVMvA64DfdHnKWsctPESkaPorao+aKu0a0qpX4BCaJ8glFLBwAlgmIgMTcrkIvKHUqoeWrH9oJRqIiI3UvB5AGoChbFRPhZT9rHoLcVjQBa0CX0kegVkyKCYlY8hvdAZuIE28XXEfBz4/LiIOWjFcxb9ZRyb19AmySPQX6KB6BVIshGRO2j/IwXME5FIR82UUqUtSuIeMNQS/aE6cMGmXXbL6wWcQESOoH2LgoC1Fv+clNAK2CMif9vUXUAryX7oZzgfbbreXET2pHA+QxpG6W1tg8GQ3lBKFQI+B+4Dr6BDDj0LjBORFTbteqIVY5CDLcNUQyl1GJgrIh94SgZD2sEoH4Mhg6OU+ho4JCIfeloWg8GKUT4Gg8FgSHXMmY/BYDAYUh2jfAwGg8GQ6hjlYzAYDIZUxygfg8FgMKQ6RvkYDAaDIdUxysdgMBgMqY5RPgaDwWBIdf4PLkGtF1BFsKYAAAAASUVORK5CYII=\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([52.86E3, 56.12E3, 65.05E3, 78.26E3, 94.96E3,\n", "\t 115.3E3, 140.7E3, 173.5E3, 219.8E3, 286.7E3, \n", "\t 393.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.2) ,\n", " xytext=(4, 0.2),\n", " arrowprops=dict(arrowstyle=\"<-\"), va=\"center\", ha=\"center\", fontsize=16)\n", "\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", "\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/neptune/neptune-mb.png', dpi= 300,bbox_inches='tight')\n", "plt.savefig('../data/jsr-paper/neptune/neptune-mb.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../data/jsr-paper/neptune/neptune-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 }