{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 11 - Mass-Benefit Analysis - Venus" ] }, { "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": 3, "metadata": {}, "outputs": [], "source": [ "planet = Planet(\"VENUS\")" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "ISP = 320.0\n", "g0 = 9.80665\n", "\n", "peri_alt_km = 400.0\n", "apoa_alt_km_pp = 400.0 \n", "apoa_alt_km_ab = 60000.0" ] }, { "cell_type": "code", "execution_count": 14, "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": 16, "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": 5, "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, 10.0E3, 11)" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaYAAAGJCAYAAADBi5fgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAByOUlEQVR4nO3deZzN5RfA8c+DGTMY+76EIkpKGSI0lKUQsqQiUmhRCe30Q7tUWlQoRbSICqWEEi1UaKMkobJFluwZ5vz+OHdm7szcWe7MXWY579fr+5q53+0+d9KceZ7vec7jRARjjDEmtygU7gYYY4wx3iwwGWOMyVUsMBljjMlVLDAZY4zJVSwwGWOMyVWKhLsBoVS+fHmpVatWtq/fvXs3FSpUCFyD8pCC/NmhYH/+gvzZoWB//px+9tWrV/8jIn7foEAFplq1arFq1apsXx8bG5uj6/OygvzZoWB//oL82aFgf/6cfnbn3B/Zuc6G8owxxuQqFpiMMcbkKhaY/DB48OBwNyFsCvJnh4L9+QvyZ4eC/fnD9dlduEsSOeeqA3cDscA5QDRQW0S2ZOHaQp5rbwAqA78CD4jIO77Oj42NlYI6VmyMMaHmnFstIrH+Xpcbekx1gCuAfcDnfl77IDAGmAhcCqwEZjvnOgaygcYYY0InN2TlLReRSgDOuYFA+6xc5JyrCNwBPCYiT3h2L3XO1QEeAz4MRmONMcYEV9h7TCKSkM1LOwCRwMxU+2cCDZ1ztXPUMGOMMWER9sCUAw2A/4CNqfav83w9M7TNMcYYEwi5YSgvu8oC+yVt9sZer+Mp7N69m9jY5OdwgwcPLtAZN8YYE2hTpkxhypQpiS/LZ+ceeTkw+a1ChQoFdga3McaEgvcf/M65f7Jzj7w8lLcPKO2cc6n2J/aU9mKMMSbPycuBaR1QFDgt1f7EZ0s/h7Y5xhhjAiEvB6aFQDzQJ9X+vsBaEdkc+iYZY4zJqVzxjMk519PzbWPP10udc7uB3SKyzHPOCWC6iFwPICK7nHNPAfc65w4Ca4DewEVAl5B+AGPyi+7dYa8fo+Bly8K77wavPaZAyhWBCZid6vULnq/LgNae7wt7Nm8jgUPAUJJLEl0hIh8Ep5nG5HN798Jnn2X9/Natg9WSJNOmTWPAgAFJr0uUKMGpp57KoEGDuPHGGylSJLf8Gks2ZswYxo4diz8l31p7fpaf+fPzz6dyxX9REUmdwJClc0TkJPCQZzPG5GOzZ8+mevXqHDhwgNmzZ3Prrbeya9cuHnjggXA3LSBeeOGFzE8qIHJFYDLGmMw0atSIOnXqANC+fXs2btzIM8884zMwxcfHU6RIEdIm7eZeZ55pNQES5eXkB2NMAdakSRMOHDjAN998g3OOF154gbvuuouqVatStGhR9u/fj4gwYcIE6tWrR2RkJFWqVOGWW27hwIEDKe7lnGPkyJE8/PDDVK9enejoaC688EK+//77FOfVqlWLa6+9Nk1bnHOMGTMmw/Y+88wznHHGGURHR1OmTBliY2N57733ko63bt06aThv586dFClShGeffTbNfR5//HEiIiLYvXt30r53332XZs2aUaxYMUqXLk2vXr34888/M/4B5mIWmIwxedLmzZspXLgwJUqUAODhhx9mw4YNTJkyhffee4+oqChGjhzJ8OHDadeuHe+//z533XUX06ZNo1OnTiQkpCzT+dprr/Hhhx8yceJEpk2bxt9//83FF1/MXn+SQdLx+uuvM2LECK666io+/PBDXn/9dXr27JnuvStXrkzbtm2ZOTN1KVCYMWMGl1xyCRUqVABg0qRJ9OjRgzPPPJM5c+YwefJk1q5dS1xcHAcPHsxx28NCRArM1rhxYzGmQILwb9n06quvCiDr16+X+Ph42bt3r0yaNEkKFSokXbt2lc2bNwsg5557riQkJCRdt2fPHomMjJT+/funuN+MGTMEkHnz5nn9eJBy5crJoUOHkvZt3rxZihQpIqNGjUraV7NmzTT3S7x+9OjRSa9Hjx4teH3mIUOGyLnnnpvh54yLi5O4uLik1zNnzkz63Im+++47AWTWrFkiInLw4EEpWbKkDBgwIMW9Nm3aJBERETJhwoQM3zPYgFWSjd/V1mMypiDIaviIi/Mv3Phzfg7Vr1+fiIgIypYty80330yfPn145ZVXko5369YtxTOllStXcvz4cfr27ZviPldeeSVFihRh2bJlKfZ37NiR4sWLJ72uVasWzZo1Y8WKFTlue5MmTfj++++59dZbWbJkCUeOHMn0mssvv5wSJUowY8aMpH0zZsygVKlSdOmiM2JWrFjBgQMH6NOnDydOnEjaatSoQf369Vm+fHmO2x4OFpiMMXnCe++9x7fffsv69es5fPgwr732GmXLJtdqrlKlSorzE4fJUu8vUqQI5cqVSzOMVqlSpTTvWalSJbZt25bjtvfr148XX3yRr7/+mg4dOlC2bFm6d+/Oli1b0r2mWLFi9OjRg9dffx0R4eTJk7z55pv06tWLqKgoAHbt2gVA27ZtiYiISLH99NNP7NmzJ8dtDwfLyjPG5AlnnXVWUlaeL6kz8BKD1s6dO2nQoEHS/hMnTrBnz54UQQ3g77//TnPPv//+m2rVqiW9joqK4vjx4ynOycovf+ccN9xwAzfccAP79u1j0aJFjBgxgt69e/P111+ne90111zD9OnT+eKLLzh69Cg7duzgmmuuSTperlw5QOd6eX/GRDExMZm2LTeywGSMyZeaNWtGZGQkb731FhdffHHS/lmzZnHixImkDLhEH374IYcPH04aztuyZQsrV67knnvuSTqnZs2arF27NsV1CxYs8KtdZcqUSQpIkydPzvDcNm3aUL16dWbMmMHRo0epVasWrVq1Sjp+wQUXEBMTw8aNG+nfv79f7cjNLDAZY/KlsmXLMmLECB599FGKFy9Ox44d+eWXXxg1ahQtW7akU6dOKc6Pjo6mffv23Hnnnfz333+MHj2akiVLMmzYsKRzrrzySq677jqGDRtG586d+eGHH5g2bVqmbRk8eDAxMTE0b96cihUrsmHDBmbMmEH79u0zvK5QoUL06dOHyZMnEx8fz7Bhw1L0DEuWLMn48eMZMmQIu3fv5tJLL6VUqVJs27aNZcuW0bp1a66++mr/fnC5gAUmY0yysmX9KzNUNs16nLnKww8/TIUKFZg0aRIvvPAC5cqVo1+/fjz66KMUKpTyEXu/fv0oXrw4t9xyC//88w9NmjThrbfeSjHk179/f/766y+mTp3K5MmTadWqFe+9916GQ4wALVq04NVXX2XGjBn8+++/VK1alb59+zJ27NhMP8M111zDuHHjkr5P7YYbbqBGjRqMHz+eN954gxMnTlCtWjVatWpFo0aNsvBTyn2cBCBbJq+IjY0VWyjQGJNa4gTbhx6y6maB5JxbLSKxmZ+ZkmXlGWOMyVUsMBljjMlV7BmTMabAK0iPNPIC6zEZY4zJVSwwGWOMyVUsMBljjMlVLDAZY4zJVSwwGWOMyVUsMBljjMlVLF3cGOPbDuAlYAngAAHaAQOBKhlcZ0wOWY/JGJOSAE8AfYGzgE+AZZ6vDTz7n/CcF0KDBg3COZeiqGqotW7dmpYtW4bkvWrVqpVmkcPUxowZk2a5j/zAApMxJqUngT+AxUB3IMKzP8LzerHn+JOha9LRo0d5++23AZIKlRoYOHBgQFbYzW0sMBljku0APgKeIf3fDoU8xz/ynB8Cc+fO5cCBA3Ts2JFdu3axcOHCgNxXRNIs/Bcs//33X8DvWb16dZo1axbw+4abBSZjTLKXgSFk/puhEHAzMDXoLQJg+vTplClThmnTphEdHc306dPTnPPuu+/SrFkzihUrRunSpenVqxd//vlninMSh8deeeUV6tevT2RkZNJCfwsXLqR58+ZER0dTqlQpunXrxq+//uqzPfPmzeOss86iaNGi1K9fP6k3lyhxiG3t2rV06NCBEiVKcMUVVwCwaNEiOnbsSJUqVShWrBhnnXUWTz75JCdPnszwZ3Dy5EkGDx5MyZIlWbJkSYr38eacY9SoUTz77LPUrl2bmJgY4uLiWLduXZr7jRo1KqkdF110EevXr8c5x5gxYzJsS7BZYDLGJFsMXJbFc7t4zg+y7du3s2TJEnr37k2FChXo1q0b77//Pvv27Us6Z9KkSfTo0YMzzzyTOXPmMHnyZNauXUtcXBwHDx5Mcb+lS5fy1FNPMXr0aBYuXMjZZ5/NwoUL6dSpEyVKlGDWrFm8+OKLrF27lpYtW7Jt27YU12/cuJHbbruNESNG8O6771KnTh2uvPJKli5dmqbtXbt2JS4ujvnz5yc9G9u0aRMXX3wxr7zyCgsWLKB///6MGTOGkSNHpvszOHr0KD169GDevHl89tlntG3bNsOf2cyZM1mwYAHPPPMMr776Kn/++Sddu3ZNMQQ6evRoHnnkEfr168e8efNo3749Xbp0yfC+ISMiBWZr3LixGFMghf3/vuw3fdy4cQLIV199JSIiCxcuFEBefPFFERE5ePCglCxZUgYMGJDiuk2bNklERIRMmDAhaV/NmjUlOjpaduzYkeLcxo0bS506dSQ+Pj7F9UWKFJFhw4Yl7YuLixNAVqxYkbTvxIkTUq9ePWnZsmXSvtGjRwsgTz/9dIafLSEhQeLj4+Whhx6S0qVLy8mTJ1O0tU+fPrJ3715p0aKFnHrqqbJx48YU1ye+jzdA6tSpI8ePH0/aN3v2bAHkyy+/FBGRvXv3SvHixeWmm25Kce2TTz4pgIwePTrDdmcVsEqy8a/FekzGFASSxa0VcDyL5x4HLvTj3tk0ffp06tatS/PmzQFo27YtVatWTRrOW7FiBQcOHKBPnz6cOHEiaatRowb169dn+fLlKe7XrFkzKleunPT68OHDrFmzht69e1OkSPIMmtq1a9OiRQuWLVuW4voaNWqkeK5TuHBhevXqxTfffENCQkKKcy+//PI0n2fHjh3ccMMN1KxZk8jISCIiIhg1ahT79+9n165dKc7dvn07LVu25MiRI3z11VecdtppWfqZtWvXjoiIiKTXDRs2BEga2vzpp584fPgwvXr1SnFdz549s3T/YLPAZIxJ1g54P4vnzvecH0SrVq3i559/pnv37uzfv5/9+/dz8OBBunfvzsqVK9mwYUPSL/O2bdsSERGRYvvpp5/Ys2dPintWqZJyEta+ffsQkTT7ASpXrszevXtT7KtUqVKa8ypVqsTx48fZvXt3hu+VkJBAly5d+OCDDxg1ahSffvop3377bdIw3rFjx1Kc/+OPP/Lzzz/Tu3dvn++bnrKplrwvWrRoivvv2KFZKxUrVsz0s4WDTbA1xiQbiM5T6kbGf7YmAC8Arwe3OYm9onHjxjFu3Lg0x1977bWkeUXTpk2jQYMGac6JiYlJ8Tp1skCZMmVwzrFz58401+7cuTPNL/m///47zXl///03kZGRVKhQIcP3+v3331m1ahUzZsxIMUfp/fd9/zVwySWXcM4553D33XcTFRXF0KFDfZ7nr8SAuWvXrhQ/M1+fLRwsMBljklUBLgWGkn7KeAJwm+e8yj6OB8jx48d58803Of/883nsscfSHB82bBgzZszgzjvvJCYmho0bN9K/f3+/36d48eI0btyY2bNnM2bMGAoXLgzAH3/8wVdffcWtt96a4vy//vqLlStXJg3nnTx5ktmzZ9O0aVMKFcp4EOrIkSMAKYbZ4uPjef319CP8nXfeSeHChbn99ttJSEgIyATjhg0bUrx4cWbPnk2bNm2S9s+ePTvH9w4EC0zGmJRGoJNn26Ep4V3QybXx6PDdC2hQGhHcZixYsIA9e/bw5JNP0rp16zTHb7jhBm666SbWrFnD+PHjGTJkCLt37+bSSy+lVKlSbNu2jWXLltG6dWuuvvrqDN/rwQcfpFOnTnTu3Jmbb76ZQ4cOMXr0aEqVKsWIESk/aKVKlejduzdjx46lQoUKvPjii2zYsIEXX3wx0890xhlnULNmTUaOHEnhwoWJiIhgwoQJmV43fPhwChcuzLBhw0hISEjTJn+VKVOG22+/nUceeYSYmBjatm3LmjVrmDpV8/8zC7DBZoHJGJOSA+4A+qDzmp71OtYOmElIauVNnz6dmJiYNA/oE1111VUMHz6c6dOnM23aNGrUqMH48eOTKkNUq1aNVq1a0ahRo0zf65JLLmHBggWMHTuWK664gsjISFq3bs3jjz9O1apVU5xbp04d7rrrLu677z5+++03atWqxZtvvpmi55GeyMhI5s6dyy233EK/fv0oW7Ys1113HaeccgqDBg3K8NqhQ4dSpEgRbr31Vk6ePMldd92V6ftlZOzYsYgIU6dO5dlnn+X8889n2rRptGjRglKlSuXo3jnlpACtdR8bGyurVq0KdzOMMSZXmjNnDr169WL58uW0atUqx/dzzq0WkVh/r7MekzHGFEBff/01CxYs4PzzzycqKorVq1fz2GOP0axZs5AVqk2PBSZjjCmASpQowfLly3n++ec5cOAAFStW5IorruDRRx8Ne8VyC0yZsTVpjDH5UIMGDfjss8/C3QyfbIJtenLpmjTGGJPfWWBKTy5ck8YYYwoCC0y+5NI1aYwxpiCwwORLLl2TxhhjCgILTL74WpPmd+Bz4J9U+4O9Js0O4AG0inOc5+uDWC/NGJNvWVaeL47kZ0qJfgQeB35Bf2pnAPU9X/cCm4GaBC7UC/r86iO093YvyWVh3keTLxLLwoQ3s9MYYwLKeky+CBoAvF0OrAD2AWvRXsx5aE9qK9qTKQE0Aq4CxgJvowHtGP6z5AuTSxw7Btddp1+NCQULTL5ktCaNQysqtwFuAi5Cey1/AX+jz6c6oYuozQKuBEoDpwGd0RpkU4Ev0Z6WL5Z8YXKRl1+GdetgapifpQ4aNAjnXECqa+clY8aM4dNPPw13M0LKauX5sgMdKltM5mvStEPXpMmo/H88sAlYjw4F/uL1fRQphwXre963OZCVxSTf8dxnVBbONcZPx45BXBwsWgTt28OyZRAVFfp2HD16lMqVKydVKNi2bVuK1WbzM+ccI0eO5KGHHgp3U/yW3Vp51mPyxXtNmoR0zvFnTZoIoB7QFbgHmA58DfwLfA/8Dzgb2IhO2p0I9EOHCvugyQ5zgAM+7h3s5AtToL38MvTrB6VKwTXXhK/XNHfuXA4cOEDHjh3ZtWsXCxcuDMh9RYTjx48H5F55yX///RfuJmRMRArM1rhxY8myBBEZLyIXicgcETnu2X/c8/oiz/GErN8yyy4UkX9F5GsRmS4i94jIpSJSTkTuEpGtPs43JsCOHhVp2lS/+nodSh06dJAyZcrIrl27JDo6Wnr27JnmnHfeeUfOP/98iY6OllKlSknPnj3ljz/+SHFOzZo1pU+fPjJ16lSpV6+eFClSRN59910REfnoo4+kWbNmEhUVJSVLlpSuXbvK+vXr07zPu+++KxdccIEUL15cYmJipEmTJjJv3ryk488995w0a9ZMypQpI6VKlZLzzz9fPvjggxT32Lx5swDy/PPPy7Bhw6RChQoSHR0tnTp1ks2bNyedhz7xTrGNHj1aRETi4uIkLi4uTftq1qwp/fv3T3r96quvCiDLli2Tnj17SqlSpeScc84REZH4+Hh55JFHpF69ehIZGSlVqlSR4cOHy9EA/UcGVkk2fldbjyk9iWvSzAR+Btqi6dpt0aGzmZ7jwciIEyAaaIr2nB4FPgRWAf8BDYFrgZ9Im6RhTIAk9pYSh+6iosLTa9q+fTtLliyhd+/eVKhQgW7duvH++++zb9++pHMmTZpEjx49OPPMM5kzZw6TJ09m7dq1xMXFcfDgwRT3W7p0KU899RSjR49m4cKFnH322SxcuJBOnTpRokQJZs2axYsvvsjatWtp2bIl27ZtS7r2ueeeo3v37lSsWJHp06cze/ZsLr/8crZs2ZJ0zpYtWxg4cCCzZ89m1qxZxMbG0rlzZ5+9vEcffZTffvuNV199leeff57Vq1fTvn174uP1f+wVK1YAcO2117JixQpWrFjBwIEDs/Vz7NOnD7Vr12bOnDlJKwL37duXhx56iKuvvpoFCxZw7733MnXqVPr06ZOt9wiY7ESzvLr51WMKpwdE5J0Mju8RkUdEpIqINBKRARKcnpvJNyD8W3aNGzdOAPnqq69ERGThwoUCyIsvvigiIgcPHpSSJUvKgAEDUly3adMmiYiIkAkTJiTtq1mzpkRHR8uOHTtSnNu4cWOpU6eOxMfHp7i+SJEiMmzYMBER+ffff6VEiRJy+eWXZ7ntJ0+elPj4eGnXrp106dIlaX9ij+mMM86QkydPJu3/4osvBJCXX345aR8gI0eOTHNvf3tMt99+e4rzli9fLoBMnz49xf6ZM2cKIN99912WP2d6sB5TPjIQeJ70n2+VRec1/Y6mon8OxAJvAidC0UCT1/gbRp57DiZO9P9YRlt2TZ8+nbp169K8eXMA2rZtS9WqVZk+fTqgvYoDBw7Qp08fTpw4kbTVqFGD+vXrs3z58hT3a9asGZUrJz8YPnz4MGvWrKF3794pEipq165NixYtWLZsGQBfffUVhw4dYvDgwRm2d/Xq1XTu3JlKlSpRpEgRIiIiWLx4Mb/++muac3v27JliGfMWLVpQvXr1pJ5SIF1++eUpXi9cuJDIyEh69uyZ4ufWvn17gDQ/t1CywJQbZTX54k7geuBXdN7UJKAOmkZ+KPjNNPnTsWMwYwZcf73v4wMHwmuvhWZe06pVq/j555/p3r07+/fvZ//+/Rw8eJDu3buzcuVKNmzYwK5duwANWBERESm2n376iT179qS4Z5UqKder2bdvHyKSZj9A5cqV2btX53Uk3qd69erptvevv/7i4osvZu/evTz33HN89dVXfPvtt1xyySUc8/EDq1Spks993sOHgZL68+3atYvjx49TvHjxFD+zihUrAqT5uYVSwci3zItGoJNn26H1+LqQXPlhPvACKSs/dPZs3wDj0Uy+wcCt2LpRxi+pny2l5v2saciQ4LYlsVc0btw4xo0bl+b4a6+9lrTa6rRp02jQoEGac2JiYlK8Tr0IXpkyZXDOsXPnzjTX7ty5k7JlywJQvnx5ALZt28ZZZ53ls70LFy7k33//5e23304RwI4cOeLz/L///tvnvkaNGvk831tUVBQHDqRN1U0MpKml/tzlypUjKiqKzz//3Of5VatWzbQNwRL2HpNzroZzbo5z7l/n3AHn3LvOuVOyeO0pzrnpzrk/nXNHnXMbnHMPOeeKB7vdQZfd5IumwGw0Hf0AunbU9Z5rjMlEZr2lRKHoNR0/fpw333yT888/n6VLl6bZGjVqxIwZM2jevDkxMTFs3LiR2NjYNFu9evUyfJ/ixYvTuHFjZs+ezcmTJ5P2//HHH3z11Ve0bt0agAsuuIASJUowZcqUdO+VGIAiIpJrmm3YsIEvv/zS5/lz5swhISF5WOTLL79k69atScOWAJGRkRw9ejTNtTVr1mTDhg0p0t2XL1+eJtkjPYm9uH///dfnzy2cgSmsyQhAMeA3tMhPN3Smz0/o05PimVxbHNiAVqnrj9ZiuAs4CszydU2eSX4IpH9EkykqiUhnEflMLFHCpOu550QmTgz8udnx7rvvCiDTpk3zefzFF18UQD799FOZNGmSFC5cWG644QaZO3euLF26VGbOnCmDBg2S119/PemaxHTx1D766CMpVKiQXHLJJTJ//nx54403pG7dulK+fHnZtm1b0nnPPfecANK9e3d55513ZNGiRfL444/Ls88+KyIia9eulSJFikj79u3l448/lmnTpknNmjWldu3aUrNmzaT7JCY/VK9eXS677DL54IMP5NVXX5XKlStL3bp15fjx40nnNmrUSOrVqyeLFi2Sb7/9Nqk9n376qQDSp08fWbx4sUyZMkXq168vpUqV8pn88Ntvv6X53FdddZWULl1aHnjgAVm4cKEsWrRIpkyZIt26dZNff/01a/+hMkA2kx9yGljOAnoAdbJ5/VDgpPf1QG30Ef7wTK5tjyZWt0+1/zHP9cVSX1MgA1OiIyIySUTqikgTEXlbROIzvMIUQNddJ3LBBSJxcZlvF1yg5wdL165dJSYmRg4fPuzz+P79+yU6Ojrpl/CCBQukdevWEhMTI9HR0VKnTh0ZMGCArFu3Luma9AKTSNp5TF26dPE5j2n27NnStGlTiYqKkpiYGGnatKm8//77ScdnzZol9erVk6JFi8qZZ54pb775pvTv399nYEqcx1S+fHmJjo6Wjh07yqZNm1K83xdffCHnnXeeFC1aNMU8JhGRSZMmSZ06dSQqKkqaN28uq1atSjcrz1dgOnnypDz99NNy9tlnS9GiRaVkyZJy9tlny5133in79+/3+XPyR3YDU7ZLEjnnbgYeQR+91wOuF5F3/LzHJ0CUiLRItX8ZgIjEZXBtZ7SiXXMRWem1/x5Pu2JE5LD3NVkuSZSfJaDPqMajpZeGAwPQ/mdqO4CXgCXokKGgz7wGYs+tjMmBLVu2ULt2bV566aVsz0vKC8JRkuhm4GwROR9ohf6K81cDdBgvtXXAmZlcuwQdBhznnDvTOVfCOXcR2gublDooGY9C6KDpl+hzqk+BWsD9aBFa0AD0BFov8CzgE2CZ52sDz/4nPOcZY0yA5SQwxYvInwAi8hNaq8BfZdGFJFLbC5TJ6EIROQa0RD/DOuAg+qvzA+AWX9fs3r07xcO9jB5iFggXAO+iQeoftIDsDcDd2JIbxphsmTJlStLvWKB8du6Rk6G8X4ArSM4Lm+X9WkR+zMI9jgNPicg9qfY/BNwjIummszvnotBFH6qiydF/ojlp/wNeF5GbUl9jQ3mZ2IU+oZuIpqLfBbRI59zEyuozsWE9Y4xP2R3Ky8k8pmj0aYW3xNcCnJqFe+zDd88ovZ6Ut+uB1mjixO+efcudc/8CU5xzk0Tkhyy0wSSqiP7XmI7+9Pt79t2B5ksW9jq3EDqYOxVbcsMYE1DZHsoTkVoiUjudLStBCXQILu2MOH2+9HMm1zYE9nkFpUTfeL6ekcU2GG+L0XWgbkbTWkagS8qfCXyR6lxbcsMYEwSZBibn3F7n3Hme719xztUO4PvPB5o555ICmXOuFjqAlLo3ltpOoIxzrk6q/ed7vga+pkdB4Eh+plQYnQywAhiHBqyR6Oq8eJ1njDEBlJUeU3GgqOf7a4EKAXz/l4AtwDznXFfnXBdgHrpQ+eTEk5xzNZ1zJ5xz//O6dhqa8PChc66/c66Nc+5ONF9sNfpI3/hLSLuUhkMz+b4HvkOTJtb7OM8YYwIgK8+Y/gAGOecSg9O5nsQDn0QkyyVpReSwJ8V7AjAD/RX4CXC7iHiXIXXo3++FvK7d4pxrBowBHkKzP/4CpgAPi0h65U9NRtqhs8O6+zhWGViAFottiQartiFrmTGmgMg0K885dx3ae8msd+XQRVcKZ3Je2FhWXhbsQOcpLSbj/+K/oIOmTYDXyXx5eWNMgRO0rDwRecU59xFwOrAUuA0rCZp/eS+58Qy+g1MCul7USHR5jUbAi8DlPs41xhg/ZSldXER2ADucc9OBBSKyObjNMmHl75IblwLXoFObnwZi0tzRGGOyzK90cdFFvC0o5Xf+LrlxAZoYAdp7+ip0TTXG5D+2UKBJXxW0ht79WTg3Bp1s+x6aODHYc52llBtj/BT2hQJNPnM5mlL+LTob7dfwNscYk/dYYDKBVwX4EJ311gJNjLBK5MaYLLLAZILDoYkTX6BDfJeRvKyGMcZkwAKTCa76aDJEI882L5yNMcbkBQELTM65/znnBmVUFcIUUJFobY45wDBgEDr/yRhjfAhkj2kMWiHiD+fc3QG8r8kvWqBp5SeBc4GVYW2NMSaXCmS6eG204GtL4MIA3tfkJyWBV9CVc7uhK+aOwtLKjTFJAtZjEpE/RORnEZkiIn0DdV+TT3VH08q/RntSG8LbHGNM7mHJDyZ8qgAfAf3Q4DSFlGnlO4AH0P53nOfrg579xph8y++hPOdcf+Aq4BQgdaKDiMhpgWiYKSAccAtwMdAHXXLjJbTs0UfAEOBekmv1vY9WP/eu1WeMyVf8CkzOufuBscBa9DH2f0FokymIzkCTIcagdewvJO3SGxHoEGA3tPr5k2jNPmNMvuJvj+l64BkRGRaMxpgCLhK4FfgYWAfcBDyFptR4K4QuydEO7WVVCWEbjTFB5+8zpnLoYIoxwfEyus7TD2h/vBGaIJFaIbSyxNSQtcwYEyL+BqZlwDnBaIgxgA7fXYamlU8DHkPXg3qEtPX2unjON8bkK/4GptuBAc65fs658s65Qqm3ILTRFCSOlHOaegBr0HlP1wHHvY7Z3Cdj8iV/A8kG4CzgVbQkZ3yq7Xj6lxqTBYL+S/JWDe2r/wN0Av717E99njEmX/A3+eEBbAEDE0zt0KeY3VPtL44uQjgUaAUsAL7xnG+MyVf8CkwiMiZI7TBGDUTnKXUjbX++CDARTRNvjvakrFq5MfmOPRMyuUsVdPLsUCDBx3EHDEcHlBNn0xlj8hW/A5Nzropz7gnn3LfOud89Xx93zlUORgNNATQCqIkO071D8rOkeM/rdkBbdL7TtWiKuTEm3/C38sPpwOdAGeBLYCNQGf37tp9zrpWI/BbwVpqCxaEVHfqgQedZr2Pt0HJFiZNqP0d7WJvRNZ+sRJExeZ4TyXoug3PuPXQQpZ2IbPHaXxNYBKwTkdSPrXON2NhYWbVqVbibYQJtNzqn6VR0SY2i4W2OMUY551aLSKy/1/k7lNcGuN87KIEueYFWOWvjbwOMybEKwKfAMaADsC+8zTHG5Iy/gSkSOJjOsYOe48aEXjTwNnAeuoTGlrC2xhiTA/4Gpu+BW1NXeHDOObRy2feBaZYx2VAYLfp6IxqcbNTWmDwpOxNsPwB+cc7NQpdsqwz0Auqi8/KNCa/b0NXCLkWfOV0W3uYYY/zj7wTbhc65zmj+00g0B0qA1UBnEVkU+CYakw3dgKpAV2AUuuCgMSZP8HsFWxFZCCx0zhVD08b3iciRgLfMmJxqik5q6Ig+cxqHTSk3Jg/I9v+mInJERLZZUDK52qnAV+iaTleimXvGmFwt08DknDvpnGvq+T7B8zq97UTwm2yMn8qis+wKAxejVcqNMblWVobyHgC2en1v1cVN3hMFvA7cB1wAfAjUCWuLjDHpyDQwichYr+/HBLU1xgRTIXRF3Nro0hnvolXKjTG5il/PmJxzm5xzPpdWd86d5ZzbFJhmGRNENwBT0TJG74a5LcaYNPzNyqtF+pXIotCa0Mbkfh3R6uRdgD+B28PaGmOMl+xk5aX3jCkW2J/9phgTYueh6eQvofXxT4a3OcYYlWmPyTk3DBjmeSnA+86546lOi0Zzn94KbPOMCbKaaHDqDvQA3gCKhbVFxhR4WRnK2wR84vm+P1qBbHeqc/4DfsaWbDN5UWlgIbqsexvgfaAiWnDrJWAJyTVO2nnOq+LrRsaYQMhKVt48YB6A1mrlQRGxJAeTv0QC09HFW5oBPdFCW0OAe4EIdAXd94G+aB2+EdjChMYEgV/PmERkgAUlk285YCz6tHQi8D90iC/CczzC83ox8AfwZBjaaEwB4G+6+N3OuefSOfasc+7OwDTLmDDZAewB3kNr5vt6aloIeAb4yHO+MSag/M3KGwD8mM6x7z3Hjcm7XkaH7zqgz5buQou/ps5FLYSuQDY1pK0zpkDwNzCdAvyWzrFN2Dwmk9ctJnn9prOBFcCb6Dyn1MGpi+d8Y0xA+RuYjgDV0jlWHc3OMybvciQ/UwL91/4Z8A26Mm6C1zHv84wxAeNvYPocuNM5l6L6g+f1CM9xY/IuQbPvvJVGq5OvB64FEmvopz7PGBMQ/gamMegS6huccw875252zj0MbPDs/1+A22dMaLVDU8JTi0GTHXYCV6NBab7nfGNMQPm7tPoPzrk2wBPA3WhgSwC+AHqIyA+Bb6IxITQQnafUjbR/thVDg9EVaNr4IfT5kzEmoLKztPo3wIXOuWiSl1Y/GvCWGRMOVdDJs0PRlPDUwSkKeBs4C114sGRIW2dMgZCTpdWPish2C0om3xmB5pe2A94h+VlSvOd1J2Aw0MTz/cEwtNGYfMzvHhOAZ02meujfjymIyGs5bZQxYeWAO4A+6LymZ72OtQNmoj2rBDRTrwO6Im7pkLbSmHzLiWR9pXTnXGlgAVpNDJIrhSXdREQK+9UA52oAE9D/5R06rfF2Efkzi9efgS753gYojq6u84KIPJP63NjYWFm1apU/zUty7BjcfDO88AJEpQnHpsASdI7TF2jmXrmwtsaYXMU5t1pEYv29zt+hvEfQ//UuRIPI5cBFwOvoBNum/tzMOVcM+BSoj1YuvwbN7lvqnCuehetjga/RxQsHosu/PYmO/gfUyy/DunUwNUwz/Y8dg+uu068mF3HA0+ifVa2Bv8PZGGPyB38DUwc0OK30vN4qIp+JSD+0pzPUz/sNAk4FuonIXE8l8y7oCP8NGV3onCsEvAZ8IiJdPNcvFZEpIvKUn+3I0LFjMGMGLFoEr70WnuAQzsBoQTETDngUrUgeB2wLb3OMyev8DUxVgE0ichI4hs7uSPQu+ijYH12AlSKyMXGHiGxGl27rmsm1rYEzgIAGIV9efhn69YNSpeCaa0IfHMIdGK23mAUOGA1ch44nbAlra4zJ0/wNTDtJfsT7B9Dc61idbLx/A2Ctj/3rgDMzubal52uUc26lcy7eObfLU+U8Ohtt8SkxKFx/vb4eODD0wSGcgTHcQRHCHxj9che63nMc6VeVNMZkyN+svC/QxIcPgBnAaOdcLbRIS390+qE/ygL7fOzfi86RykhVz9dZ6Oo596Ar6TwA1ECff6Wwe/duYmOTn8MNHjyYwYMHZ/gmiUHh44/1F3R0NBQpApdcAs2a6evsbFFRUCgLfxYkBoZly/T1wIEQF6eBMhRJGL6C4pAhwX/fRN6BsX370H3uHLkFzVdtgyZEZPYnljH5yJQpU5gyZUriy/LZuYe/WXmnAVVF5HPnXATwGNAbnRO/ELhVRPb4cb/jwFMick+q/Q8B94hIuoHTOTcFfUb1nIjc5rX/bk+7zhSRX7yv8Tcr79gxDQLLlsFff8F338HRo3DgADzxhA4vnTih+/zd/vsPihaFYsUyDmB//QUREdC4cXJA+/57iIyEdu3SBruMAmFkJDg/Vlz1/vxRUWlfh8LEidrmIUNSfh9K2c7InAnciZYyahSUphmTq2U3K8/fkkS/A797vo9HpyKO8PdNvezDd88ovZ6Ut8QAmHrhgUVoYDoX+IUcSOwtREVB3bq6JXIuZ78kExL0F15Gwevff2HUKBgxImUArF9fexGlS0N8fMprMrpnQkL6wcvX/t9+g5gYeOSR5H21a8OgQXDZZb6vTf190aL+BUNv4e4tJvIeSvTrv3dftOfUAa2/51fOqjEFV5YDk3MuEh02myAiywP0/uvQ50ypnQn8nIVrM5KQyfEMpf6lmFpOf0kWKqS9pWLF0j9n4kQYNgxuuintsZo1/Q+MJ05kHLi8jx04AEuW6Oc8cQL274edO7W98+fD4cNw/HjKa3x9Hx+fHKwyCmC+guQPP0DVqhoQSpaEKlWgbVt47jm4447sBzx/5HgosScanDqj6UEtMz7dGOP/UN5B4DIR+Swgb+7c7WhB2NNFZJNnXy30sfE9IvJkBteWQxNzXxaRW7z234umtNf1zvYD/4bysjJsFMyhpcyGzYI9rJbRZ/Pnc588qW3NLICl/v7gQXjpJejZUwPjv//Cjh2wdSv8/rs+56tSRQNXlSrJW+rX5cpl7VleVn4OOfrvvRitJPEWOvPPDza52+RV2R3K8zcwLQQ+E5HH/H2jdO5XHPgBOAqMQufRP4imoZ8tIoc859VEhxAfEJEHvK4fDdwPPI5O1I1Fk3Znici1qd8vq4Epq7/0gxkcwhkYwx0UIfPAGB8PXbposErctm9P+/2hQ1CpUuZBrEIFKJxqWnbAn7EtA3oB09Cp4H78LGbM0GHlUD9fMyYnQhWYGgBz0brLc4EdpFpwWkT8GkJzzp1CypJEn6AlibZ4nVML2AyMFZExXvsdmpx7M7rs+w5gOvCg5xlYClkNTP78wg9GcAh3YMxPvcX//tMhSO+g5SuI7d+vwck7aP31l2YjXn01VK6s++fM0WSUbH/ulegMvUn4yBtN/7MmDiWGMvHEmJwKVWBKDDrpXSQZZdKFW1YD0/XXw/r1+gsoM/HxmowQyDk24QyM4Q6KEJ7AGB8Pf/+dHKj+/BMeeww6dIBduzS47dyp54AmwlStmhywKldOu5Utm85zsDXoVPQJwJUZtys3ZCUak12hCkyjMztHRMb624hQyUkR11AKZ2As6L3FROl9toQEePxxfQbWunVywPLeduzQr4cP6zCir6BV+ThUeRgq3wmVb9Vkj8w+YzjS9Y3JiaAFJufcXqCtiKxxzr2CDpNtzmY7wyqvBKZwKsi9xUSBGkr87z/tYfkKWjt3ws4tsPMn2FkIIqNTBq4qVWDTJh1K7NsXzjgDatSA55+3XpPJO4I5j6k4Wr0btLrDJPR5j8mHwl3257vvNDDOnp35uYmBMdC856/5EhWVtSoYRYvCKafolq5NIBfBgUGwo0dy0PrrL3j7bbj4Yhg/Hn7+GY4cgQYNYONG7bk1bgxnnaWp9MFiGYEmHLLSY9qAliJ6Dc18u4kMJq4GcI5TwFmPyWQmLEOJfwEXo3/2jdRdvnqD//wDP/2kvabt2zUVf906KF8ezj4bGjZM/nr66ZpSn1OWEWhyIphDedcBk8m84KtDkx8CvhZSoFhgMpkJ21DiTqAt0BWOjYK41rDsLYiagS4o49CUo3ZwrC/EXalBMTJSh/x+/FGDVuLXrVuhXr3kQJUYtCpXzvrEZMsINDkV1OQH51wV4HRgKXAbGfeYPvG3EaFigclkJqzP2P4B2sHE0uBKw5ADwBDgMiACiEdLGz0PE0uCuxiG3OL7VocP6/Cfd8D68Uc9lrp31aABFPexLKdlBJqcClVW3qvoJNc8+YzJApPJ7Y7tgLhasKwORP2A76fACXBsCMQtgGUbst6LEdHnV949qx9/hF9/hWrVUgarevW0SPHy5ZYRaLIvJIEpr7PAZHK7iQ+DmwZDKqMrnL0M+BocT4CJZ4DrD0Puy9l7xsdrwV7vgPXFF1o1o1EjDVbnnquloKpXh+HDc/Z+puCwwJQFFphMbnd9Y1h/GCIqoEtoRqDrNPt4LhS/E+oXh6mrA9uGxN7R/PmaAfjjj5otuWqVLrlSp45mBJ57Lpx3nn4tVy6wbfBui2UF5l0hWfbCGBNcU4ujZYsigGNodfIItPhr0VQnx6MJEwGWmC5fqZJuLVokH3vmGZ2bVa+eBqsFCzRYlSqlAcp7q1Ej5xXgs73kiMnTrMdkTG4ShxZ7TXQcuAr4D5iDLqGR0fk5lJ3JxQkJsHkzrFmjwSpxO3EibbCqWzdtsdzM2mJZgXmXDeVlgQUmk+tdiJYx9s4KjEeXzDiIrukU7bW/LQENTIGsU7hjR8pA9d132ttKfGaVOAzYoIFORs7ofSwrMG+ywJQFFphMrvcgunRm91T7TwD9gN3APKAY8A46cWNUYN46FJOL9+/XoT/vYPX77zoh2LtnVb8+dO5sdQLzulClixcCConICa99HYCzgE9F5Dt/GxBKFphMrrcDXZJ9MWmntJ8EBgBb0eDUDXgdqByYtw7X5OKjR2Ht2pTBas0aKF0a2rSBCy7QobwlS3TRR+s15R2hSn54Ex3t7ud50xuBFzzH4p1znURkib+NMMZ4VAEuBYaiq555B6fCwKvAQKA+ugpZgIIShK9OYXQ0NGmiG2jv6MILYdIknSS8fDk89ZQ+y/rvPw1YHTtCmTKBeX+T+/jbY/oDuFtE3vK8/h0dER8BTAEqi0ibYDQ0EKzHZPIEAZ4EPkKDTxeSKz/MR/8UPIwGqo+AIBZxDQdfvTERnQg8ahT88osWuW3QQHtSHTpA06aBqQ1oAiu7PabM6t+lVhHY5nnDOkBtYKKIHET/lmvobwOMMak44A5gJvAzmuAQ5/n6i2f/V8A5QAfg3/A0MxiOHdOisddfn3K/c9pDmzkTSpTQwPTQQ8nznCpUgB49YMoU2LIlLE03AeRvYDoAJE6law38IyKeClycJG0yqzEmu6oA96NZd4nbKM/+QsDzQBOgHbAvTG0MsKwuOTJzpi4JMm6cJlP88gt066bDfk2b6jyr226DDz7QChbZdeyYlmY6diz79zD+8zcwfQXc45zrDNwOfOh1rA76WNYYEwoOfQ7VCu1N7Qlvc3Iqvd5SagMHwmuvpQwWlSsnB6ydO+Gtt6BqVXjySV108aKL4LHH9DlaQkLW2+Q9wdeEjr/PmOqiweg0YBO6su0Wz7FPgT9EZEAQ2hkQ9ozJ5EsC3AN8jC6RUT68zcmuYGUFHjoEn30GH3+sk3X374d27fTZVPv2Wt3CF5vgm3MhncfknCsnIntS7WsI7BSR3X7fMEQsMJl8S9Bhv/locKoY3uZkR6iWHNm8WYPNxx/Dp59CrVoapDp00PJLiZN9bYJvztk8piywwGTyNQHGArPRXNkAppLnVydOwNdfJweqn3+GVq10/tSMGbBypaaz2wTf7AlVYJoF/CciPucxAbl6HpMFJlMgPIRm7n0KVA1zW/KYvXvhk09gwgRNqChVSpMfbroJZs2yXpO/QpUu3oyUCQ93oivGlEKreI30twHGmAAbBVyL5s1aOpJfypaFyy6Dkydh+3aYN0+XqT/9dFi9GiZPtgy9ULB5TMbkR/cAN6Dzn/4Mc1vymMSU9ehoOOccnRv1669Qu7bOkTrvPFi8WCf9muCweUzG5FcjgNvQ/1O3hLUleUbqlPWOHTuybt06KlaE//0P/vxTz7n9dq2S/sor1oMKBn+LeCTOYzqBzWMyJvcbiv5f3hpNiDgtrK3J9bwn+K5cuZKPPvqI4sWLM9tTQLB06eSl5evV02dR992nz6Buugkq5sFsyNzI3x7TXWiPaT7aOxrjdaw3sCIwzTLGBMwQ4D6gDfBbmNuSi6XuLY0dO5bKlSuzePFi/vwzeTx04EA9r1Ur+PBDTTnftk0D1cCBOiHX5IxfgUlEfhORukAFEamTOLnWYygauIwxuc1g9M/INsCv4W1KbpW6t/Tzzz9ToUIFunfvzqOPPpp0XmJZpMQ5VGeeqc+hNmyAmjWhbVu45BJNQbfnUNljCwUaU5BMR3tPi4Ezw9yWXMZ7gu+PP15K+fLd2L79RerUmcC6dT1p3Pg7oqJOATKe4HvsGLz5pg7zJSTAsGHQp0/BnP8U6soP5wD18JHsICKv+X3DELHAZAy6uOCdaAkjy6NNY+XKlfTu3ZvffvuNpk2bMm3aNGbNmsX+/ft58cUXs3wfEZ0T9dRTuvDhjTdqJfSC9BwqJPOYnHOlnXNfAt+hiwZO82yvem3GmNysDzABaA/8EOa25EJjx47lvvvuIzIyMmnfiBEjePvtt1M8a8qMczqs9+GHsHQp7Njh33OoglzZ3N/kh0fQ5IdWaG3jy4GL0L/BNgFNA9o6Y0xw9AYmous5rQlzW3KZhg0bMmBAylrU5cuXZ/z48ezfvz9b9zzjDJ2c6/0cqkMHLYOU3qBVQa5s7m9Jot/RalyvoyWImojIas+xF4HiieWKciMbyjMmlbnoRNwP0LWdTAqNGjVi2rRpNGrUKKD3/e8/fQ711FNaZWLYMOjbN/k5VH6pbB6qkkRVgE0ichI4BsR4HXsX6ORvA4wxYdQNmAp0BlaGtykFSdGicO218MMP8Mwz8N572pMaMwb+/js5Q7BUqZQZgAWFv4FpJ1Da8/0fQHOvY3UC0SBjTIh1Rp8UdwG+9OzbATwAXIiWNboQeNCz3wRM4nOoBQt0zajE51Bjx0KzZnqOr4UR8zt/A9MXaCFXgBnAaOfcZOfc88B4NM/HGJPXXIpWJO8G3Az0RRez+QRd0v0ToIFn/xPoEhsmoBKfQ91zjwalzp11PtQvvxS8XpO/JYnGklxIfzyaCNEbKIZWg7g1cE0zxoRUe6AHmlv7AXCx17EIoDsauIYCTwJ3hLh9BcCxYzqst2yZ9qamToVLL9UVd9eu1blWefFZk7/8rfzwu4h87vk+XkRGiEh1ESkrIlenXtXWGJOH7EBLFn0IXIWuhJtaIeAZ4CNsWC8IvKtPFC2q8542bIAaNfRrt25w8GC4Wxl8/g7lAeCcq+Gcu8A5d1HqLdANNMaEyMtoXb02aCrT1cBCH+cVQof7CtDQUiikrtWXqGRJeOQR+O47WLUK6tbVIb8TJ3zfJz/wd4Ltqc65FWgR/c/Rv6mWoAVOEr8aY/KixcBlnu9bAvOAfuiwXmpdsP/bA8y7t+TL6adr1l7//ppqfs45Onk3P1aV8/cZ08vAKeiSF+uB44FukDEmTBz6LClRczQoXYZOxu3ldcz7PJNjib2lZcsyPm/gQJ3f9Nlnuljh8OFak++JJzRQ5Rf+BqYmwLUi8k4wGmOMCSNBp817B52maM/oEuAgcJ1nf3xom5bfZdZbSpRY2fyVV2DIEE2MmDJFJ+F26gQPPgjVqoWmzcHkb2DaivWSjMmf2gHvo9l33s4GPvMcP4hm5c33vDYB8d13Wtncsx5hhhIrm4NWQh8yRKtGPPqorqo7ZAjcdReUKBHcNgeTv4HpEeBu59ynInI4GA0yxoTJQHSeUjfSPn0+HX2q3Bb4F53b9HooG5e/5XSOUqlS8NhjuoruyJH6PGrsWC0CW7hwYNoYSv6mi89A/0lucc6975x7LdU2PTjNNMYEXRV0ou1QIMHH8VPQntMzwEmgUshaZrKoZk2YORPmzdOvjRrBQl+Zlbmcv1l51wL3omWJzkOrjKfejDF51QigJjpM9w7Jz5LiPa+vQQPXEeAmNECZXKdJE02QeOghGDpUK5n/+GO4W5V1/s5jGgu8hy6tXk1EaqfaTg1CG40xoeLQig4zgZ/Robs4z9dfPPv/h5Yo+hVNJ7dEiFzJOejaVStGXHaZVo+4/nrYvj3cLcucv4GpHPCCiOwPQluMMblFFeB+dOA+cRvl2Q+6rsCH6POmnuhaAyZXioiAW26BX3+F8uWhYUOdD3XoUMbXhXOhwuwUcT0jGA0xxuQx0WiFiCi0Qnkmv+hMeJUuDePGwerVWt6oXj1NujiZznBsOBcq9DcwDQUGOef6OOfKOecKpd6C0UhjTC4VCbwB1EKLwO4PZ2NMVtSqBW+8ocVip02Dc8/VBQm9JU74XbQoPEtu+BtIfgEaAq8Bu9DRZe/N5jgZU9AUBl4Czgdao78ZTK7XtCksX65p5UOG6BIba9fqsXAvVOjvPKYHsJVYjDGpOeApND3qQrRaRI2wtshkgXNw+eVaNWLSJLjoIv3+hx/gq6/0nMQySKFccsNJmCsAOudqABPQBFWHFoO9XUT+9PM+9wCPAl+KSEtf58TGxsqqVaty2GJjTIaeRGvrLcbWtc5j9u+Hnj1hxQqtHnHHHVC8OEycqEFsyBD/7uecWy0isf62I6zPhJxzxYBPgfpAf3SWRF1gqXOuuB/3ORXNGbJBBGPCbQQ62zEOWBvmtgBffvkl7du3p2LFisTExHDeeefxyiuvJB3/448/6Nq1KzVr1iQ6Opry5csTFxfHhx9+mOI+R44c4frrr6ds2bKcdtppzJo1K817Pf7445xzzjmcyKNrUkRF6XpPq1driaSRI3V/yJd3F5GwbWgyxUmgjte+2sAJYLgf9/kYmIzOS/8ivfMaN24sxpgQeUNEKonIN+Frwg8//CBRUVHSunVrmTt3rixatEgGDx4sgLzwwgsiIrJ27Vq57rrr5LXXXpNPP/1U5s6dK506dRJA3nnnnaR73X///VKjRg2ZN2+ePPjggxIRESEbNmxIOv7XX39JTEyMfPnllyH/nIHy3HMiEycmvz5xIv1jWQGskuzEhuxcFKgNnab3pY/9y4BlWbzH1cBuoKwFJmNymfkiUkFEPgvP2997770SEREhBw8eTLG/WbNm0qxZs3Svi4+Pl+rVq0vnzp2T9sXGxsq4ceOSXtevXz8puImI9OjRQ6677roAtj60jh4VadpUv2bnuC/ZDUzhTu9ugO/O/jrgzMwuds6VQZ9P3SUiewPcNmNMTl0GvIWu5fRhJucGwfHjx4mIiCA6OjrF/lKlSpGQ4KsgoCpSpAilSpWiSJHk/LDjx4+nuE+xYsU45hnbWrhwIZ999hmPP/54gD9B6KReeuOmm26iadOmSccTl9wISYZedqJZoDY0vfwxH/sfAk5k4fqX0ZrHiUkcn5FBj+mUU06Rxo0bJ22TJ0/Oeug3xmTfChGpKCJvh/Ztf/rpJylWrJjcfPPNsm3bNtm3b59MmTJFihQpIm+++WaKc0+ePCnx8fGyY8cOGTt2rERERMiSJUuSjg8aNEgaN24s27dvl4ULF0qhQoXkq6++kmPHjkmdOnXkpZdeCu2HC6DUvaFDhw5JyZIlJTIyUn755Zd0z/Nl8uTJSb9jgS2SB4fysh2Y0IKxx4GzvPZlGJhsKM+YMPpeRKqIyNTQvu0333wj1apVE3Sqi0RERMjLL7+c5rwRI0YknVOiRIkUz5dERLZu3SoNGzZMOufOO+8UEZGxY8dK8+bNJSEhISSfJxhSPz964okn5LzzzpMmTZrI1VdfneG5GSGbQ3n+zmMKtH1AGR/7y3qOZWQyMBXY6pwr7dlXBCjseX1URP4LUDuNMTl1DrAUrRCRuOBgkP3222/06NGDBg0aMGnSJKKjo5k3bx433ngjUVFR9OnTJ+nc22+/nSuvvJKdO3fy2muvcfXVVzNnzhw6d+4MQLVq1fjhhx/YtGkTpUuXply5cmzatInx48fzxRdfcPToUYYPH857771HsWLFGD58OLfeemvwP2QAeC9UePLkYb7+ejxVqgzmv//+4vvvP+C339ZTrJiuTui9UGGwZDqPyTn3qR/3ExG5OMtvrveOlFTzjpxzn3naFpfBtZlNwBomIk9777B5TMbkAn+g1cr7AyPR2YtB0qtXL9asWcP69euJiEheM75Pnz58/PHH7Nq1i0KFfD9qb926NTt37mT9+vXp3r9Tp06cfvrpTJgwgZEjR/Lxxx/z0UcfsW3bNlq1asXcuXO5+OIs/0rMFZ588klWrlxJ27ZtWbNmDTVr1mTdunW8/rr/K0MGcx5TIfSfTuJWHy08Ugst41jL87oe/v8Tmw8088xDAsA5Vwto4TmWkTY+th/QZIo2wBw/22KMCYWa6JPhWcBdBLWWzE8//cQ555yTIigBNG3alD179rBrV/pTH2NjY9m4cWO6x+fOncv333/PAw88AGgCRP/+/alQoQKNGjWiffv2LMxjq/QdPnyY8ePH87///S9p3y233MKiRYsyDNCBlmlgEpHWItJGRNqga1fGA81F5FQRaS66BlNzz/5n/Hz/l4AtwDznXFfnXBdgHvAXOlQHgHOupnPuhHMu6aclmoCaYkNLSP7reb3Vz7YYY0KlMsnLaQRxwcHKlSvz/fffc/x4yjKeX3/9NVFRUZQtW9bndQkJCXzxxRecdtppPo8fOXKEoUOHMmHCBGJiYpL2Hz58OOn7Q4cOkdmIVG7zyiuv0KpVKxo2bJi0r2TJkgwbNoxx48aFrB3+PmN6ELhfRL723ikiXzvnxqBJC/OyejMROeycuwhN+Z6B9rg+QUsSeRfRd2ipyHCntxtjAqUs+n/7ZeiCg9OAiIwu8N8tt9xCr169uOyyy7j55puJjo5m/vz5vPnmmwwbNozIyEjGjBnD3r17adGiBZUrV2bnzp1MnTqVb775hjfeeMPnfR988EHq1avHFVdckbSvbdu2TJw4kfr167N9+3Y++eQTRowYEdgPFGRt2rShZ8+eafbffvvtrFmzJnQN8SdTAl0O7NJ0jnVEEw7CmumX0WZZecbkQkdEpKOIdBGRxDTk7SIyVkRaiciFnq8PePb76cMPP5S4uDgpX768lChRQs455xx5/vnn5YSnrMG8efOkTZs2UqFCBYmMjJRTTjlFLrvsMvniiy983u+XX36RmJiYFFUfREQOHjwo1157rZQpU0aqVq0q48eP97+xucikSZNk8ODBOboH2czK86uIq3PuF+BXEenm49g84HQRybULCVrygzG51HGgL7AHuAitoDkE7U1FoA8K3geeBy5F6/EFMWnCwOTJk1mzZg2TJ0/O/OR0ZDf5wd+hvLHA6865tWhywd9AJXRx5fpAnwyuNcYY3yKBN9Gn1RPR2i/ej38igO5ANzTN/EngjtA20YSOX89sROQtoAPwL1o/+HnP1/1ABxFJW27XGGOyYhcQA1yB9pp8JcwVQlOsPgJ2hK5pJrT8TiYQkSUi0gJNFa8MRItISxH5JOCtM8YUHC+jw3dPA13R2i5/+TivEHAzOr0+yBo1asT3338f/DcyKWSr8oNz7hx03lKU53XSMRF5LSAtM8YULIuBe9BnR2OBkmhwWoyu0uatC/AsugqbyXf8CkyeUj8L0JFgIfnxo3cGhQUmY4z/HCnTxUcApdCl2t9Ff+skCnBaucld/B3KewQoh/4d44DL0dHg14FNQNP0LzXGmAwImn3nbSA6ZNeFlLVcUp9n8hV/A1MHNDit9LzeKlploR+whJCUZTTG5Evt0JTw1DoCi4DbgSfQADbfc77Jl/x9xlQF2CQiJ51zx9AcmkTvokuCGWOM/waic5m6kfZP5nOBFUAn4HfgV8B3UQaTD/jbY9oJlPZ8/wcpR33rBKJBxpgCqgo6eXYo4Gtx2RrAcvQp9w6gROiaZkLL38D0BdDM8/0MYLRzbrJz7nlgPPBxIBtnjClgRqDVx9sB75D8LCne87oHmlJ+AZoUsT0MbTRBl53KD1U9349HEyF6A8XQUd+8sSqWMSZ3cmhFhz7ovKZnvY61A2aiPSsBHkXHbD4AGmLyEb8Ck4j8jo7wIiLx6N83eat8rjEm96sC3O/ZfHHAfUBt4GI0L9iSIfKNcC+tnjt17w5792b9/LJl4d13g9ceY4xvVwHV0WqdjwLXhbc5JjD8DkzOuSpoLykOLbO4F1gKPCUiOwPbvDDZuxc++yzr57duHayWGGMy0wpNiuiIzqZ8EKs8nsf5W/nhdHRR5DLAl8BGtF7eUKCfc66ViPwW8FYWNNZjM8Y/9dB08i7AZuAVoGhYW2RywN8e0zjgAHC+iGxJ3Omcq4lOgRuHFqc3OWE9NmP8VxEdu+kLtAfeI+XSGSbP8DddvA26tPoW750i8gcwxnPcGGPCIxqYDTRBU8o3hbc5Jnv87TFFAgfTOXbQczx/mjMHZszQYbNy5VJ+3bcPvv8+eV+xYuBskNuYsCiEli6qDbREe07nh7VFxk/+BqbvgVudcx+JSNLcbKfrXtzsOZ4/xcZCkSI6zLZ3L+zZA3/8oV///BOuvVa/37sXTpxIGbh8BTNf+6KjLaAZEyhD0Mm6nYEpaMlpkyf4G5geQKez/eKcm4UWBqkM9EJXTOkU2OblIrVq6eZL69YpnwkdO5YygKX+ftOm5O+9jyckaJDavx/i4tIGrjPOgCZNoEqVoH9cY/KFzsBCNCliC1oI1v72y/X8nWC70DnXGXgIGIn+JxZgNdBZRBYFvol5UFQUVK2qmz+OHtUA1a0bjB2bMnDt3g1Ll8K332rPqkkT3WJjtYdmjPGtMfAV+mfzJnSF3MLhbJDJjN/zmERkIbDQOVcMTRvfJyJHAt6ygig6GqpVg+LF08+0E4HNmzVAffstPPwwrFgBdeokB6smTeC88/Q+/rJUdZMf1UQrffZEh/TeBLLxv4cJjWxXfvAEIwtIoeYcnHqqbr176764OHjxxeRgNWsWrF2r53j3rM4+G4pmMrnDUtVNflUa+BC4AS0P8D5a+sjkOlaSKD9wDs48U7f+/XXf8ePw00/JwWrSJNi4Uc/x7lmdcQYUtnENU0BEopNvH0ILwC4AGoS1RcaHTAOTcy4BfY6UFSIieT/YlS3rX0+gbC6cxRcZCY0b63bjjbrvyBH47jsNVEuWwKOPwo4dcO65yYHq6FEdLrTsQJNfObQ4bG105uVbwEVhbZFJJStB5AGyHpjyh/z6zKRYMWjRQrdE+/bB6tUarN5+G374QbMAY2NT9qyqVQtfu40Jhr5oAdjewONA//A2J7e54YYbcnwPEWlM+vHjbzSrO41MA5OIjMl+s0y2hLLHVqYMtG2rG+j7vvUWrFqlwWrKFBg0SJ9Nde6syRGtW2uPzJi8rjXwGZqxtxkYjfaodgAvodVAr0cTJdqhy7/n0udSW7duZdy4caxatYoffviBo0ePsnnzZmqlmuayefNm7rzzTpYsWUJ8fDxNmzZl/PjxxMbGpnvvt956i6uuuopq1aqxdevWFMeefPJJJkyYQHx8PNdddx0PP/wwhQolFxX6+uuvadu2LT/99FPqtlRK9w1FJMsbmnR5DVDUn+tyy9a4cWMxmYiLS7svIUHk119Fxo0TadZMpEwZkb59Rd55R6Rly5A30ZiA2ykiTUSkr4g8KiIXicg7IuecfY589913Isf1tVwkIuNFJCFcDU3f0qVLpWLFinLppZdK+/btBZDNmzenOOeff/6RqlWrSr169eStt96S+fPnS+vWraVEiRLy888/+7zvvn37pFKlSlK5cmWpVq1aimOffPKJFClSRCZPnixvvPGGxMTEyKuvvpp0/MSJE3LuuefKI488kl6zfcea9A74PFn/tkgA9gBPAfX9uT7cmwWmLPAVmFLbulXk+edF2rYVKVxYpGtXkenTRfbsCXbrjAmewyLSQESqicg/uuucczyBKdFJEblFNDjlMidPnkz6/qWXXvIZmB588EEpXLiwbNy4MWnfoUOHpGLFitKrVy+f9x00aJC0b99e+vfvnyYw3XHHHXLppZcmvb7xxhvliiuuSHr93HPPSf369eX48ePpNdvn72q/iriKSGvgTGA60A9Y55z7zDnX2zkX4c+9TB5WrRrcfDMsXgzNmunw3nvvQe3a0K4dvPACbN8e7lYa459/0QrlPdA1nrZAVFQUUVFRyecUAp4BPkKH+3IR7+Gz9KxcuZK6dety2mmnJe0rXrw4rVq14oMPPuBEqsn6X375JTNnzuT555/3eb/jx48THR2d9LpYsWIcO3YMgL///pv777+fF154gYgI/8KDv9XFEZH1IjIcqAZci86hfgPY6px7zDl3qr/3NHlYkSLQr58Gpu3bNWCtWAFnnQXNm8Pjj8NvtkSXyQNeBm5BA88NwAWw6NFF1K9fP+V5hdDKoFND3cCcK1y4MJE+ng8XLVqUo0eP8vvvvyfti4+PZ/Dgwdx5553UqVPH5/3OP/98lixZwpo1a9i4cSOzZ8+mWbNmANxxxx106tSJNm38X3QiJxNs/wNmOOfWocN6FwJ3AXc4594DbpX8sqKtyZrixeHyy3WLj9eJuu++CxdeCOXLJx9r1Cj9dHSrPGHCZTFwj+f7oUAtKHlFSXgWXcLdWxd0/6jQNS8Q6tWrx+LFi9mzZw/lypUDICEhgW+++QaAvV7/740bN47//vuPe++9N9379e7dm7lz59K4cWMA2rRpw2233cayZcv44IMPWL9+ffYamt4YX0YbuurJdcA3wEngZ+BWNPVvELAN+CQ79w7mZs+YsiArz5j8Pf/kSZEvvxS54w6R004TqVVLZNgwkeXLRU6cCP77G5MVF/rY952InCYiN4vIsSycn0uk94zp999/l6JFi0rHjh1l48aNsn37dhkyZIgULlxYAFm5cqWIiPz2228SFRUlH330UdK1vp4xJdq2bZv88ccfIiJy/PhxadCggUycOFFERJ5++mmpXbu2VKhQQW644QY5cuSI96U+f1f7u7R6Q7ST2wdNoJwH3C0iS71Oe8k5txNdrsvkNcFIVS9UCC64QLfHH9eKFO+9B7feqhN8u3bVntRFNsvRhJEA8YD345BGwCpgADom9DZady8+1I0LjFNPPZXXX3+dIUOGJA3PnXfeeQwbNownnniCKp6VC2677TYuuugimjVrxv79+wF9niQi7N+/n6JFi6Z4tlTVq2D1008/TVRUFDfddBOLFy/m/vvvZ/ny5VSrVo0OHTrwyCOP8OCDD2bc0PQilq8Nzcjbimb7V8ngvDOApf7cOxSb9Zhyod9/F3niCZEWLURKlxapWFHk7bdFDh7M2vXWYzKB8oBoSrgvCSLyuIhUEpGPRGSOiDwYonZlQ3o9pkQnTpyQn3/+OSk778Ybb5QaNWokHa9Zs6agodrnNnToUJ/3/euvv6RkyZLyzTffiIjI8OHDpXv37knHn332WYmNjfW+JOc9JrQ27zwROZlJsPsFW2bdZMWpp8KIEbrt3KnPo6ZOhYEDtTjt5ZfDZZfpMypjgmkgWg2iG2nTwhxwJ7oS7pVAUbRaeR5VuHBhzjjjDAC2b9/OrFmzuPPOO5OOv/XWW0nZdYkee+wxVq9ezezZs6levbrP+95+++306dOHJk2aJO07fPhw0veHDh1K7LxkyN/1mOwpswmeypV1DauFC3WxxAULdMjv9tvh/PM1WHXrZlUnTHBUAS5FEx+ewXfOckvgEmAZOrz3OlAhVA3M3Jw5cwBYvXo1AB999BEVKlSgQoUKxMXFER8fz1133UVcXBwlS5Zk3bp1PProozRo0IARI0Yk3Scxs87btGnTKFq0KK3TGer/+OOP+fzzz/n111+T9rVt25ZnnnmGF154gapVq/Lcc89x7bXXZv5B0utKZbQB5wBXoHOZUmzZuV+oNhvKywN8Dc0dOSLyxhsibdroUN8dd4isX5/++cZkV4Lo5NmLRIfrEueFHve8Tqz8cFxE7hGRGiLyVeibmR7SGXqL8/x/Eh8fL506dZKKFStKZGSknHrqqTJy5Eg5fPhwpvfOKPnh2LFjUqdOHZk+fXqaY+PGjZNq1apJmTJlZMCAAanfy+fvaidZ6FYlcs6VRgvFJ4bTxJzfpJuISK5dQyE2NlZWrVoV7maYjKRepj61337Tob5p06BePdi1C9as0UUWjQmUHei8piVe+9qhdfO8a+XNR4cARwK3Ycu2+8/nT8zfCbaPAOXQ/BSHrgV5Edqh3QQ0zUEDjclc3brw2GPw559w223w999Qo4Z+/9NP4W6dyS+qoEtjLPPaRpG2gGsXYCVaC+dK4GAI25iP+RuYOqDBaaXn9VYR+UxE+qF/WwwNZOOMSVdkJPTooavyrloFpUvDpZdqtYlXXgGvB67GBNWpaHnrUkATYF14m5OH/J3eAX8DUxVgk2hW3jEgxuvYu2jxeGNCq1YteOAB2LIF7rsP5s7VXtRNN+kwnzHBFgVMAe5Fl9KYGdbWBMTkyZNzvCaTc241Orrma/O5FhP4H5h2AqU93/+BLk6cyHcxJWNCpUgRTS2fP1+H9apV0xJHjRvr0vIHDoS7hSa/6w98gi6veiP657vxm7/zmL5AEx8+AGYAo51ztYAT6H+S+QFtnSl4AlV5olo1GDVKe1CLF8NLL8G99+q8qEGDtCq6r3p9VqvP5NTZwLdo0baWwBygVjgblPf4G5jGAom1J8ajiRC9gWJoULo1cE0zBVKgf8kXKgQdOuj2998wfbpWQ4+K0gDVt2/K4LZ3b8ZZgan5E0RNwVEKDUgT0Em5r2APOvzg71Dev+jfAohIvIiMEJHqIlJWRK4WkT2Bb6IxAVKpEtx1F2zYAM8+CytXauWJvn1h2TLwY+qEMZlywHDgHbTC6Ci05LXJVKaByTlX2Dk3xjm3D82iOOCce8czp8mYvMc5aNMG3ngDfv8dYmM1UaJ+ffjrL9i9O9wtNPlJS2A1sAJoD+wKb3Pygqz0mG4E/gd8BzyBVhTvinZSjcnbypXTkkfr1iWnmdetC1dcoc+mEhLC3UKTH1QCFqHpYueRp+vshUJWAtMg4CXRYhx3i0gvYAjQ1zlnRctM/uActGihvaYtW/TZ0Z13Qp068PDDWmDWmJwoDDwETEaXb38Kr5o5xltWAtOppF1baRb6Y64Z8BYZE26lS+sS8d99B7NmaaA64wy49lrdZ0xOdAK+Bt4AegE2iyGNrASmEqT90SUW3ojBmPzKOWjSRFPNf/9dg1OXLrocx3vvwUl7km2yqRY6nFcBiAW8q2ntQOdBXQjEeb4+6NlfQGQ1K6+ac+7UxA3tRaXZ7zlmTP5TtizcfTds2qS9qccf12G+rVvh33/D3TqTF0UBL6I1+S4CpqFP8fsCZ6ETdZd5vjbw7H+CAjH8l9XANAf4zWtb79k/N9X+3/xtgHOuhnNujnPuX+fcAefcu865U7JwXaxzbopzbr1z7ohz7k/n3OvOudr+tsGYLIuIgN69YcUKeOstOHgQatfWIrK/+f3P3xi4BlgK3AW8CrwPdCd5ifcIz+vFaL2dJ8PQxhDLygTbAcF6c+dcMeBT4D+0coSgjweXOufOFpGMKnFeif4d8SxaNrEa+rfHKudcIxH5K1jtNgbQxQvPOANmzoQXXoALLtAiskOHwkUX+a4sYYwv5YAz0aG9xGoRqcefCqELGLYD+pC20nk+kmlgEpHpQXz/QeiPv56IbARwzv2I9rxuQPNW0jNORFJMOHHOfQls9tz3f0FpsTGpVa8OjzyiJZBmzoRbb4XChTUN/eqrba0ok7mX0fWcLkf/1G4GTAUuS3VeIeBmz7FRoWxgaPlbkijQugArE4MSgIhs9gSYrmQQmFIHJc++P5xzu9HekzH+y0mtvmLFYPBgLXW0ZAk8/bTW57vhBp3AW7VqurcxBdxi4B60WsRQdGW73miCxANAUa9zu6DBywJT0DRAJ+ymtg5NpPSLc+4MoCLwSw7bZQqqQNTqcw7atdPt11+1/NFZZ0HHjtqLio3N+XuY/MWR/EwJdCLuanTF3Fj02VPiP5sI8r1wB6aywD4f+/cCZfy5kXOuCDAJ2I12dNPYvXs3sV6/FAYPHszgwYP9eRtj/FOvHjz/PDz0kC4J36OHrhV1++3QrZsu1eHNqpsXTALEkzLoVED/bH8Dnft0HTAanUGai02ZMoUpU6YkviyfnXuEOzAF0kTgAqCTiPgKdlSoUIFVq1aFtlXGAJQpA3fcoQFp7lwd5hsxAm65BQYO1ONg1c0LqnYkZ+N5c2iiw8VovZ3z0Cy+diFtnV+8/+B3zv2TnXv4W1080Pbhu2eUXk/KJ+fcY8Bg4DoRWRSgthkTeEWKQM+e8MUX8M478OOPWuF8yBAd9jMF00DgeSC90oyV0Uy9/wFjgO3A0ZC0LCzCHZjWoc+ZUjsT+DkrN3DOjQTuBm4TkRkBbJsxwRUbCzNmaAHZsmWhVStdeXfRIluCo6CpAlyKJj6kF5wETYa4F9gDNAK+DEXjQi/cgWk+0My7YoRnRdwWZGE1XOfcbei8p5EiMjFYjTQmqKpWhQcfhD/+gPLldYjvrLNgyhQ4ciTcrTOhMgKtPtoOXcMp3rM/3vO6HVrKaDRarfRRNEVsGJDP/pmEOzC9BGwB5jnnujrnuqCP+/5Ca/AC4Jyr6Zw74Zz7n9e+K4GngYXAp865Zl7bmaH8EMYERHQ0VKmiw3vPPQcffAA1a2rK+bZt4W6dCTYH3AHMRMeL2qK18tqiecYzPccT5213R2vs7UKXc18W4vYGUVgDk6eyw0XABmAG8Do6QfYiETnkdapDc1G823uJZ/8l6BJc3tsLQW+8McHinFaOmD8fvvoKDh3SHtQ111h184KgClrDZpnXNgrflR7Kob81nwKuBm4BDvk4L48Jd48JEflTRHqISEkRiRGRbiKyJdU5W0TEicgYr33Xevb52lqH+GMYExx162rvadMmaNgQLrtMV999/31bxNAk6wKsRYPS2Wihtzws7IHJGJMFZcrAXXfB5s1aWWLMGK3Tt327PYcyqgxaoXwiWnn0RvLsWk8WmIzJSyIitP7eqlW6TtTevVCrltbp21GAFuwx6euI9p5OoL2nPDiBxgKTMXmRc3Dhhfrs6csvYf9+OPNMXWX3hx/C3ToTbqXQwrBT0JLWA4E8tGyYBSZj8rq6dWHiRF1lt359rcnXti18+KE9hyro2qOZexHo4oMfhrc5WZWfShIZk/flpLp52bJwzz0wfDi8/TaMHKlzooYN04w+W36jYCqJrpTbC+05XQhMwM9qpKFlgcmY3CQQBVkjI6FvX+jTR+vuTZigz6BuvFGXha9cOefvYfKei4Af0coRDdFJNV3C2qJ0WWAyJr9yTlPL27TROnzPPKOZfJdfrr2ohg1Tnm+VzfO/EsBzaO/pOuBtdFXccl7n7EBLH8xA089/QatODCRkq+ZaYDKmIKhXT5d/f/BBmDwZOnTQxInhw/V756yyeUFyIdp7Gon2niaiq+c+CXyEVjIfBvzgOfY+0Bet5zeC5OoTQWLJD8YUJOXKwX336XyoPn3g7rs1QL38siVKFDTF0GdNs9HhvUbAenQ13e4kr/sU4Xm9GPgDDV5BZoHJmIKoaFHo3x++/15X2H3vPVi5EsaOhV27wt06E0ot0F7SfrRnNAetZJ5aIXTY7yN0uC+ILDAZU5A5BxdfDAsWQKNGWkmiXj2tLvFzllaeMfnB62jvaT663lNP6NikIzfeeGPK8woBN5POGuGBY4HJGKOKFdPnTxs2wCmnaMC69FJbH6ogWAxcBpwPrAHqQY1La3DupnPTntvFc34QWWAyxqRUoQLcf78+h7riCrjzTq0q8cILWunc5D8OfZYEEAU8AizwfJ9ahI99AWaByRjjW1QUDBigz6EmTYJPP9X1oYYNg40bw906E0hC8sKEiWKBTj7OTX1eEFhgMsZkzDmIi4M5c2DNGk2caN7cloHPT9qhiQ9ZMd9zfhBZYDLGZF3NmvDYY/Dnn7oMfOIw3/PPw8GD4W6dya6BwPNAZjMGEtCKEQOD2xwLTMYY/yUuA//995owsXSpLr9x++02zJcXVUEnzw4l/eCUANzmOS/IVa0sMBljsi9x+Y05c3TZ9+hoHebr1Ak+/tgm7eYlI4Ca6DDdOyQ/S4r3vG4H1PKcF2ROCtD4cGxsrKxatSrczTAmdwpUrbyjR+HNN3Xi7rFjcMstOpk3JiZwbTXBswNdy2mJ1752wPX4XSvPObdaRGL9bYIFJmNMcIjA55/Dc8/BJ5/o0hu33KLrR5kCIbuByYq4GmOCI3GY78ILNVnixRehRQto0gRuvRXat4dCXk8TrLq58bAekzEmdLyH+Y4eTR7mK1lSq5X7W93cn/NNyGW3x2TJD8aY0ImOhuuu00SJl1+G5cs1m2/oUA1UxmCByRgTDs5Bq1Ywezb88AMUL67BqmNHWLjQsvkKOAtMxpjwqlEDHnkEmjWDXr3g3nuhfn1NmjhwINytM2FggckYkzsUKqS1+dasgVde0Yy+WrVgyBCdyGsKDAtMxpjcxTlo2RLefht+/BEqVoQuXSA2VqtMWC8q37PAZIzJvapXh9GjdQmOBx/UorE1a8L112uAKkBZxQWJBSZjTO5XuLAuWvjOO/DLL3D66fq1YUN45hnYsyfcLTQBZIHJGJO3VK4Md98N558PEyfCN9/AaafB1VdrMVnL6MvzLDAZY/Ku1q3h9ddh0ybN6rvtNu1NPfYY7NwZ7taZbLLAZIzJ+8qW1aD0448aqDZuhDPOgMsvhw8/hJMnw91C4wcrSWSMyR0CXSvv4EEtf/Tyy9p7uu463U45JTTvb6yIqzEmjwv0L/mYGBg8WLcffoCXXoJzz4WmTWHQILjsMoiISD5/717/a/WZoLChPGNM/nfOOZoosXWrJkk8/bRWnLj7bvjtt3C3zqRigckYU3BER+u6UMuXw7JlmsHXsqX2fv7+Wxc2NGFngckYUzDVqwfjx8Nff+nyG3//rRN6hw6Fn34Kd+sKNAtMxpiCLTISevaEs8+GVaugVCmdzNusGUydaiWQwsACkzHGJKpVCx54ALZsgVGj4P339VlUjx4wZ46tGRUiFpiMMSa1IkWgc2eYO1eDVMeOMGkSVK2qK+4uXGh1+oLIApMxxmSkTBktGrtkCfz8M5x3HowZAytWwM036/IcVgYpoCwwGWNMVlWposkRK1dqgKpeXYNTrVpw1126Cq/1pHLMApMxxmRHVBTcd59m8C1YoMN/3btrKaSxY2HDhnC3MM+yyg/GGJNTDRvq9vDD8PXXWgopLk6fSV11FVx5pfau0mPlkFKwwGSMMaC/7P0pM1S2bNp9zmmaebNm8NRTWuLojTc0Fb1hQw1SPXtC+fIpr7NySClYYDLGGAh8D6RwYbj4Yt1eeEEz+d58U8sgtWihQapbN63pZ1KwZ0zGGBNsRYtC167w1luwbRv06QOzZunw3hVXwD//WDkkLxaYjDEmlEqU0MD0wQe6wGHbtlpctmpVGDAAFi2CEyfC3cqwssBkjDHhUq6cLsvRqJFm9zVsCCNHak/q1lvhyy8L5CKHFpiMMSY3qFYNhg+Hb7/VSbsVK8KNN+rcqWuvhXfe0cUPCwALTMYYk9vUrQv336+9qK+/hsaNYcoUHe5r316fU23ZEu5WBo0FJmOMyc1q19ZhvY8/hu3b4YYbtOfUtCmcdRbcey989VW+GvKzwGSMMXlFTIxWOq9fH3bs0OXindMhv8qVtcDsnDl5fqmOsM9jcs7VACYA7QAHLAFuF5E/s3BtFPAg0BcoDXwP3C0iy4PVXmOMyRUKF4bmzXV75BEd2vvgAw1WAwbo/s6d4bLLtNeVkVxWeSKsgck5Vwz4FPgP6A8I8BCw1Dl3togczuQWU4FOwJ3AJmAI8LFzrrmIfB+0hhtjTG5Tq5auxHvLLTrUt3ixrif10EOaSHHZZRqomjXToOYtl1WeCPdQ3iDgVKCbiMwVkXlAF6AmcENGFzrnzgGuBoaJyEsi8glwBfAn8EAwGjtlypRg3DZPKMifHQr25y/Inx1C9PkTyyFldfNVDslbTIz2gl59FXbuhJdfhkKFtBJ64pDf7NmZDvmF7b+9iIRtAz4BvvSxfxmwLJNr7weOA8VS7R+L9sCKpr6mcePGkhM5vT4vK8ifXaRgf/6C/NlF8uHn37JFZOJEkQ4dREqUEGnbVqROHZHff09zarqfPS4uS28FrJJsxIZw95gaAGt97F8HnJmFazeLyBEf10YCdXLePGOMyWdq1oQhQ7R2344d2os6dEifSTVoAPfcA198EdYsv3AnP5QF9vnYvxcok4NrE4+nsHr16kPOOe9gvBv4JwvtTFTeOefP+flJQf7sULA/f0H+7FAAPn99qLcefmXXLl2ld9y4xEM+P3t9qLfeuV/TuV15oILn+3rZaU+4A1NIiYiV8TXGmFwu3EN5+/DdM0qvN5TVayG552SMMSYPCXdgWoc+K0rtTODnLFxb25Nynvra48DGnDfPGGNMqIU7MM0HmjnnTk3c4ZyrBbTwHMvI+0AE0Mvr2iJAb2CRiPwXiAY652o45+Y45/51zh1wzr3rnDslEPfO7ZxzPZ1z7zjn/nDOHXXO/eqce9Q5VyCHRJ1zC51z4px7KNxtCQXnXEfn3HLn3CHPv/1VzrmLwt2uUHDOtXDOLXLO7XLOHXTOrXHOXRfudgWac666c+4559wK59wRz7/vWj7Oi3LOjXfO7fD8LljhnLswWO0Kd2B6CdgCzHPOdXXOdQHmAX8BkxNPcs7VdM6dcM79L3GfiHwHzAKeds4NdM5dDLwF1AZGB6JxXhOA66MTgK8B6qITgIsH4j1yuTuAk8B9wCXAi8BNwOJUSST5nnPuKuCccLcjVJxzN6D/L64GLkf/AJwNpB6hyHecc2ejFWgi0LmW3YFvganOuZvC2bYgqIPO/9wHfJ7BeVPRn8X/gM7ADrSYQaOgtCo7OeaB3IBTgHeAA8BBYC5QK9U5tdCqEGNS7Y8GngJ2AseAr4HWAWzbUPQXcx2vfbWBE8DwcP/sQvDfpoKPff08/y0uCnf7QvhzKOP5N3aV57M/FO42Bfnz1gKOoqXBwt6eMHz+R9DHASVS7V8BrAh3+wL8WQt5fT/Q8++7VqpzzvHsH+C1rwjwKzA/GO0Ke1aeaE28HpmcswWto5d6/1FguGcLhi7AShFJel4lIpudc18CXdGgmG+JyG4fu7/1fK0WyraE2ThgrYi86Zx7I9yNCYHrgARgUrgbEiaRQDwanL39S+bTWPIUEUnIwmld0J/HLK/rTjjn3gLucc4VlQA9OklUoIZjsiEnE4DzqzjP11/C2ooQcc61RHuJQ8LdlhBqCawHrnTO/e4ZRt/onCsoP4Npnq/POueqOudKO+cGARejBacLmpAXMwh7jymXy8kE4HzHOVcNrUO4RERWhbs9weaci0SfdT4hIulNJsyPqnq28ejzxd/RZ0wTnXNFROSZcDYu2ERkrXOuNfAecLNndzxwo4i8Fa52hZHfxQxyygKTyRLnXAn0YfgJYECYmxMqd6HPMR8Od0NCrBAQA1wrIolrG3zqyda61zn3rHgeNORHzrm66HPvdcCN6JBeV2CSc+6YiLwezvYVBBaYMpaTCcD5hnMuGk3PPxWIE5GtYW5S0HmmBIxEHwgXdc4V9Tpc1DlXGjgoIvln2dBke9Ds08Wp9i9CszOrANtD3agQegTtIXUWkXjPvk+cc+WAZ5xzb2bx2Ux+sQ9d8SG1oBUzsGdMGcvJBOB8wTkXAcwBYoGOIvJTmJsUKqcCUcBM9H/MxA00jX4f0DA8TQu6dZkcz++/lBsCP3gFpUTfAOWAiqFvUliFvJiBBaaM5WQCcJ7nmav0OnARumbWyjA3KZS+B9r42ECDVRvyb3WR9zxfO6TafwmwVUR2hrg9obYTaOR5xujtfHRaSkErdxaSYgbebCgvYy8Bt6ATgEehufwPkmoCcD72PPqP8WHgsHOumdexrfl5SE9E9gOfpd7vnAP4Q0TSHMtHPgSWApOdc+XR1aF7Ae0pGM8XJ6KTid93zr2APmPqgs5jmyAix8PZuEBzzvX0fNvY8/VS59xuYLeILBOR75xzicUMIoDN6ET72kCfoLQpHz/DDAjPs4YJQDt0LtUn6MTDLeFsVyg457bge2wZYKyIjAlda3IH55wAD4vIqHC3JZiccyWBR4Ge6HPW9cBjIlIQ5nHhnLsUuBsdyo9CMxOnAJPz23NFz79pX5aJSGvPOYlJQFcDpYEfgLuD9QeaBSZjjDG5ij1jMsYYk6tYYDLGGJOrWGAyxhiTq1hgMsYYk6tYYDLGGJOrWGAyxhiTq1hgMsYYk6tYYDLGGJOrWGAyeZJz7iXnnDjnArZwm3NuTAaz4IN+b+fcXOfcvlSVzL2PxzjnDjvnpgWzjZ46aEHlnKvinEvwLMRoTAoWmEye4ymPcoXn5dUB/EX6MtA8QPfKjulouZfO6RzvCRTznJfXdQV2A1+FuyEm97HAZPKibkBJtNhoRbTqdYbS64V4HxORrWGuoL4AXQupXzrH+wF/4qO4bB7UDXi/gK1rZLLIApPJi/qj6yFdi1Z+7u990GtI6izn3MfOuUPA21k95vm+l+e8s1O/uXPuQ+fcD16v6zjnZjjnNjvnjjrnNjnnXnTO+VpkMl2eqtVvotWdy6V6z1OAOGBGKFePdc5d4pw75Jyb6Jwr5PXzq+/5+R12zv3pnBvgOf8a59x6zzVLnXOn+bhnSXTZkLme16c7595zzu1yzh3z3G92KIYUTe5kgcnkKc65qkBbYJaI7EZ/uV2WThCYByxDlyxI/Swqo2Oga9D8C/RN9f6V0OUfXvPaXRVdCuV2dA2jB4CL0R6dv6aja99cmWp/X7S6/WtprggS51w/dN2xx0TkllS9m9loD68bsBp4xTn3CLocwj3o8hj1AF/VyDuiC8wt8bxeAFTzXNvBc/1/2O+ngktEbLMtz2zAXei6WM09rzt4Xt/odc4Yz76hPq7P9JjX65eArUAhr323AyeAKhm0sQjQ0vM+5/q6dyafcR3wdap9vwArgvyzTfzZFPH8nOOBgemc089rXxnPz2QPUNJr/22ec2umusdbwBzP9+U953QJ978t23LPZn+RmLymP/CbiKzwvF4CbCfVcJ7Hez72ZeVYotfQv+Qv8tp3DfCJiOxI3OGci3TO3ecZwjqK/kL/3HO4XhbeJ7XpQFPn3Ome+zcF6hO6pIcJwFigp4i8nM45HyV+IyL7gF3AShE54HXOes/XGok7PKvCXopnGA8NZpuAx5xzg5xzdQPyCUyeZoHJ5BnOuVjgTOBd51xp51xpIAZ4F2iW+Ivcyw7Sl9GxRF8AW9BghHPuDOA80g6nPYr2JGYCnYCmQHfPsagsvE9qM4EEkpMg+qFDW7Oyca/suApYS/JQmy/7Ur0+ns4+SPkzuAjNLPwAPF1IXYRzFfpz3OB5RndT9ppu8gMLTCYvSewV3Y3+EkzcbvHsT53NllGSQKYJBJ5fmjOB7s65YmiAOkTa3taVwGsi8pCIfCoi3wL7M7t/Bu+7HVgM9PX0MHqjGWwpfvE75yo6595yzk12zhVzzt3sScxon9339rgYOAX4yDlXIof3Sq0bujLq/sQdIrJJRPoBFYBzgU+BFzyryJoCyAKTyRM8v6CvAr5GM7pSb98D1zjnXIDfegZQAu0B9QHeFZEjqc4phg7feRuQw/edji5r/yj6HMbXMN7dwP+ApWgPpCn6M+rj+Xll1zqgNVCXAAYnz3+bLiQP46Ug6ntguGfXWYF4X5P3WDqmySs6AeWAESLyWeqDzrnJwIvoL9SAEZENzrmvgcfQ502+suIWAv2dcz8BG9EgdkEO33oucAAYhj6/WejjnEIisgEd/rodzZ771zn3LVAZnfOEc64WsBkYKyJjsvLmIvKLc641GvQ+ds5dIiIHc/B5AM4HquAVmDzp+M+gw5QbgcLoNIATaM/JFEDWYzJ5RX/gIJqm7Mub+JjTFCAz0KC0Df1FndqtaFr1w+gv2Bi055JtInIUnV/lgDdE5ISv05xzp3kCyH/AGE9VjFhgp9d5xT1fd+IHEfkVnTtVE1jkmX+UE92A1SKy1WvfTjSADkd/hm+i6fedRWR1Dt/P5FFOh9GNMXmNc64y8DRwEhiClmm6AnhWROZ7nTcYDZo1fQxDhoxzbj0wU0QeClcbTN5ggcmYfM459zqwTkQeCXdbjMkKC0zGGGNyFXvGZIwxJlexwGSMMSZXscBkjDEmV7HAZIwxJlexwGSMMSZXscBkjDEmV7HAZIwxJlf5P03IxblyB7ZfAAAAAElFTkSuQmCC\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([26.0E3, 27.0E3, 28.0E3, 29.0E3, 32.0E3,\n", " 35.0E3, 41.0E3, 48.0E3, 59.0E3, 79.0E3, 108.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\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(fontsize=16)\n", "plt.yticks(fontsize=16)\n", "ax.annotate('',\n", " xy=(vinf_array[8]/1.0E3, fP_prop_ab[8]) ,\n", " xytext=(vinf_array[8]/1.0E3, fP_aero[8]),\n", " arrowprops=dict(arrowstyle=\"<-\"), va=\"center\", ha=\"center\", fontsize=16)\n", "\n", "plt.text(vinf_array[8]/1.0E3, 0.5*(fP_prop_ab[8]+fP_aero[8]), \n", "\t str(int((fP_aero[8] - fP_prop_ab[8])*100/fP_prop_ab[8]))+'%' , \n", "\t fontsize=16, ha='center', va='center',\n", "\t bbox=dict(facecolor='white', edgecolor='white', pad=0.0)) \n", "\n", "ax.annotate('',\n", " xy=(vinf_array[10]/1.0E3, fP_prop_ab[10]) ,\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*(fP_prop_ab[10]+fP_aero[10]), \n", "\t str(int((fP_aero[10] - fP_prop_ab[10])*100/fP_prop_ab[10]))+'%' , \n", "\t fontsize=16, ha='center', va='center',\n", "\t bbox=dict(facecolor='white', edgecolor='white', pad=0.0)) \n", "\n", "\n", "\n", "\n", "plt.savefig('../data/jsr-paper/venus/venus-mb.png', dpi= 300,bbox_inches='tight')\n", "plt.savefig('../data/jsr-paper/venus/venus-mb.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../data/jsr-paper/venus/venus-mb.eps', dpi=300,bbox_inches='tight')\n", "\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 5 }