{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 37 - Opportunity - Mars" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from AMAT.planet import Planet\n", "from AMAT.vehicle import Vehicle" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook simulates the atmospheric entry of the Opportunity aeroshell (MER-B)\n", "https://en.wikipedia.org/wiki/Opportunity_(rover)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"MARS\") \n", "planet.h_skip = 125.0E3\n", "planet.loadAtmosphereModel('../atmdata/Mars/mars-gram-avg.dat', 0 , 1 , 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Opportunity', 836, 88.0, 0.00, 5.52, 0.0, 0.66, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(125.0,0.0,0.0,5.55,0.0,-11.5,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (30*60.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZhcVbX38e/qIXOAQDqJnakTQhIgiCFhvkiYEWQQQQW5F0RoQPHFi3gB9YKAXgFHQBGjBCdQLqAgYbqAgCJjhwAhkIRA5oF0AhnI3On1/rFPpas7PXdVnVNVv8/z7KeqTp06tU53smr1rr3PNndHRERERES6piTuAERERERECoEKaxERERGRDFBhLSIiIiKSASqsRUREREQyQIW1iIiIiEgGqLAWEREREcmAsrgDyJT+/ft7VVVV3GE0q7a2loqKirjDyAqdW37SuSXLtGnTVrp7fgXdRZ3N2fn4+22OziNZCuE8CuEcID/Oo7WcXTCFdVVVFTU1NXGH0ayJEycmNrau0rnlJ51bspjZgrhjyLXO5ux8/P02R+eRLIVwHoVwDpAf59FaztZQEBERERGRDFBhLSIiIiKSASqsc6C6ujruELJG55afdG6Srwrl96vzSJZCOI9COAfI//Mwd487hoyYOHGiJ31MjohIc8xsmrtPjDuOXFLOFpF81VrOVo+1iIiIiEgGqLAWEemi+vp6tm7dGncYBc0dtm0LTUQkqVRYi4h0wcsvv0xpaSmHHHJI3KEUtOefh7Iy+OQn445ERKRlKqxFRLqgvLwcQD3WOVIg04JEpECpsBYR6QIV1rlhFm5VWItIkqmwFhHpAhXWuZEqrEVEkiznhbWZHWNmL5lZVdq208xsnpmtMrObzaws7bn+ZvYrM7vRzH5iZuW5jllEpCWFXlgnLWerx1pEkiynhbWZDQT6AAekbRsGnAqcDvwncD7w9bSX3QtMcfcrgDnAdTkLWESkDd26dQNgy5YtMUeSeUnK2RoKIiL5IKeFtbu/DzzYZPNw4Dx3n+buvwd+ARwBYGYHAWPc/aVo34eAr5lZ71zFLCLSmkLusU5SztZQEBHJBzkfCuLu9U0e/9Pd69I2LQEWRvePBOal7bskultUK5SJSHIVcmENycvZ6rEWkSRL4uTF/YHbo/uDgVVNnl8HVOY0IhGRFhR6Yd0OOcnZGgoiIvmgrO1dcsfMRgG17v56tMmBjU126w7s8AlWW1vLxIkNnSLV1dVUV1dnK1QREaDzhfXkyZOZPHly6mH/zEaVG7nM2RoKIiJxam/OTkxhbWalwIXAlWmblwCj0/YxoDewtOnrKyoqqKmpyXaYIiKNpArrLVu24O5YOyvA9ELSzFZmLcAsiStnq8daROLQ3pydpKEglwM/cvfNadueAMakPa4CtgCqoEUkEUpLSykpCal027ZtMUeTUznN2RoKIiL5II7rWKe6cyxt27eA6UBPMxtpZueZ2Sh3rwGWmdm4aNfjgVvdvfCuayUieauQx1knJWdrKIiI5IOcDgUxs77A2dHDc83s54SvEq9vsussd58S3f88cLWZzSdcT/XbuYhVRKS9ysvL2bx5M1u3bqVnz55xh5MxSczZ6rEWkSTLaWHt7uuAX0Yt5XtRa+k1C4AvZzk0EZFOK9Qe6yTlbA0FEZF8kKQx1iIieamQV19MCg0FEZF8oMJaRKSLCrXHOonUYy0iSabCWkSki1RYZ5+GgohIPlBhLSLSRSqss09DQUQkH6iwFhHpIhXWuaMeaxFJMhXWIiJdpMmL2aehICKSD1RYi4h0kXqss09DQUQkH6iwFhHpIhXWuaMeaxFJMhXWIiJdpMI6+zQURETygQprEZEuUmGdfRoKIiL5QIW1iEgXafJi7qjHWkSSTIW1iEgXqcc6+zQURETygQprEZEuUmGdfRoKIiL5QIW1iEgXqbDOHfVYi0iSqbAWEekiFdbZp6EgIpIPVFiLiHSRJi9mn4aCiEg+UGEtItJF6rHOHfVYi0iSqbAWEekiFdbZp6EgIpIPVFiLiHSRCuvs01AQEckHKqxFRLpIhXXuqMdaRJJMhbWISBdp8mL2aSiIiOSDsly/oZkdA3wP+Ly7z4+29QRuAtYB/YAr3H1tW8+JiCRBIfdYJyVnayiIiOSDnPZYm9lAoA9wQJOnbgOec/dvAQ8Cv2rncyIisSvUwjqJOVs91iKSZDktrN39fUKi3c7MKoGzgIeiTU8CJ5tZVWvP5SRgEZF2KNTCOkk5W0NBRCQf5HwoiLvXW+Pv9CYBK919Q/R8nZnNAw4Htrby3Pxcxi0i0pJCLawhOTlbQ0FEJB/kvLBuxmBgVZNt64BKoK6V5xqpra1l4sSJ2x9XV1dTXV2d2UhFRJrRmcmLkydPZvLkyamH/TMfVdbEmrPVYy0icWhvzk5CYe3AxibbuhN6Pmjjue0qKiqoqanJfHQiIm3oTI91eiFpZiuzElh2xJKzNRREROLU3pydhMJ6CWHmeLq+wFJCAm/pORGRRCjkoSDNiCVnayiIiOSDJFzH+mmgMrpEE2ZWDgwDnm3jORGRRCiywjrWnK0eaxFJspwX1tYwC8YA3H05MBU4Ktp+NHC/uy9p7bkchiwi0qrUGOtCLKyTkrM1FERE8kFOh4KYWV/g7OjhuWb2c3evBS4GbjKzfYFBQPoMltaeExGJXarHutBWXkxSzi6LPq22bcvE0UREsiOnhbW7rwN+GbX07R8CF7TwmhafExFJgkIdCpKknB19KUCB/e0iIgUmCWOsRUTyWqEW1kmiwlpE8oEKaxGRLlJhnX0qrEUkH6iwFhHpokKevJgU0d8ubNmiCYwiklwqrEVEuqhQJy8mSWlpaO6awCgiyaXCWkSkizQUJDc0HEREkk6FtYhIF6mwzg0V1iKSdCqsRUS6SIV1bqiwFpGkU2EtItJFmryYGyqsRSTpVFiLiHSRJi/mRo8e4XbDhnjjEBFpSUYKazM7IRPHERHJR/k2FCRfc/ZOO4XbdevijUNEpCUdWtLczI4HvgMMoqEoLwEGAj0zG5qISH5IamFdaDk7VVivXRtvHCIiLelQYQ38AfgBMB2oj7YZcFYmgxIRySdJLawpsJytwlpEkq6jhfV77v6TphvNbFaG4hERyTsJnrxYUDl7553DrQprEUmqjo6xvtjMTm9m+8mZCEZEJB+VlpZiZtTX17MtWcsCFlTOVo+1iCRdRwvrQ4HfmNm2tFYP/DILsYmI5I2EDgcpqJydKqzXrIk3DhGRlnR0KMh1QDVQA6S6ZQz490wGJSKSb8rLy9myZQtbt26lR+q6cPErqJy9227htrY23jhERFrS0cL6UeBBd9+cvtHMfp25kERE8k9Ce6wLKmd/7GPhdtmyeOMQEWlJRwvrPwNXmdnf07YZcBpwacaiEhHJMwmdwFhQOTtVWC9dGm8cIiIt6Whh/d/AHsC5gEfbjHCN1LxL0iIimZLQ1RcLKmerx1pEkq6jhfUPgUfdvdG6V/m6ipeISKYkdChIQeXsyspwu3QpuINZvPGIiDTV0cJ6UdMEHRnU1UDMbE/gEmAuoYdlsru/ZmY9gZuAdUA/4Ap318WWRCRRElpYF1TO3mmncC3rNWtgxQoYODATRxURyZyOXm7vK003mNkZwM8yEMsfgP9x958CNwB/irbfBjzn7t8CHgR+lYH3EhHJqIQW1gWVs81gjz3C/blzM3VUEZHM6WhhvauZ/TeAmfUys98AvwFeyUAsewJ9o/sbgZ3NrJKw9O5D0fYngZPNrCoD7ycikjGpyYsJG2NdcDk7VVi/806mjigikjkdKqzd/UTgCTO7FXgNGA3sCzS3sldH3UdYyKAvcDbwNWASsNLdN0TvXwfMAw7PwPuJiGRMSUlIp+7exp65U4g5e9SocKseaxFJolYLazP7ZNMGdAMWAO8A1wDDgaszEMtXgE2EnpSP3P1+YDCwqsl+64DKpi+ura1l4sSJ29vkyZMzEJKISHZMnjx5e74C+mfimMWQs9VjLSJxaG/Obmvy4k2Er/tWA/XNPD+FcOmmjwH/2blQt+tJGKM3GPiZmc0jXB5qY5P9ugM7DGKsqKigpqamiyGIiORGdXU11dXVAJjZygwdtuBzdqqwnjOnwy8VEem09ubstgrr64Aad1/R2k5mlomvFf8IfMHdV5uZExL21wmzytP1BbQ8gIgkSkKGgBR8zh43DkpK4M03YdMmSM7q8SIibQwFcfdH2krQ0X73dSUIM+sP7Ovuq6NN3wPWEr6+rIwu34SZlQPDgGe78n4iItliMV5cuRhydp8+sNdeUFcH06dn6qgiIpnR0auCZMsHwCYzG5y2bRXwBjAVOCradjRwv7svyXF8IiLSINacfcAB4fbllzN5VBGRruvoAjFZ4e71ZnYqcLWZTQMGAt9097VmdjFwk5ntS1jUoDrOWEVEil3cOfugg2DKFHj2Wbg07xZmF5FClojCGsDdXwcubGb7h8AFuY9IRKT9EjLGOmfizNlHHx1un3oKtm6FaG0eEZHYdXkoiJmNN7M/mtk3tHCLiBS7OMdYt0ch5OwRI2DMGFi7Fl56Ke5oREQaZGKM9UXA0e7+Y+CQ1KQVERFJpILI2ccfH26nTo03DhGRdJkorC8BJkT3/wzoSzkRkeQqiJx98snh9sEH441DRCRdhwtrM/uUmX0xun8AsEdqxre717v72gzHKCKSeNu2bQMaljZPikLN2Z/8JPTrB7NmhSYikgQd+gQws+8D9wKfAXD3l4H/MLOzshCbiEje2LRpEwA9ErRiSSHn7LIyOOmkcP/ee+ONRUQkpaNdK4cAQ4B/pm17FLghYxGJiOShzZs3A9C9e/eYI2mkoHP22WeH29/8BqIvDEREYtXRwvqJaKWt9OtKZWJpXBGRvJbQwrqgc/ZRR8HIkbBwITz2WNzRiIh0vLB+z8y+BYwxsy+Y2f3AV4EbMx+aiEj+SOJQEAo8Z5eUQHW0/Mxtt8Ubi4gIdLCwdvc/A38HSoH/AFYDx7r7L7IQm4hI3khij3Ux5OzzzoPu3eGRR+Cdd+KORkSKXYdXXnT3F4EXU4/NrNTMjnL3pzIamYhInnD37T3WSSqsofBzdkUFnHUW3Hkn3H47/PjHcUckIsWs1cLazP4G9GnjGP2BlUBBJGkRkY7aunUrAGVlZZSWlsYWR7Hm7AsvDIX1H/8IP/gBdOsWd0QiUqza6rF+I9pnNlAPXABMB2rS9hkB7JqV6ERE8sD69esB6NWrV8yRFGfOPuAA2GcfmDED7rgDLr447ohEpFi1VVj/ELBoVjlmtru7X910JzPT2lciUrRWr14NQL9+/WKOpDhzthlccw2cfjpcfz2ccw7E/zeOiBSjVicvuvuaVIKODDSzRt9zmtkk4IAsxCYikhdShfUuu+wSaxzFnLNPOw0mTIBly+DWW+OORkSKVUcnL94HvG5mzwHrgD2B44AdekRERIpFUgrrZhRNzjYL46uPPTb0Wp95JgwbFndUIlJsOnq5vSeAYwkTX8YQLt10urv/IAuxiYjkhaQW1sWWs485JgwHWb8evva1uKMRkWLUmcvtLQW+k3psZiVmNsHdp2U0MhGRPJHUwhqKL2fffDM8/jj87W/w1FNhdUYRkVzpUGFtZlOa2VwBbANOzUhEIiJ55oMPPgCSV1gXY86urISrroJvfQv+67/glVfCCo0iIrnQ0XSzB7AQWJDWugMF2fMhItIey5YtA6CysjLmSHZQlDn70ktDgf3qq3DPPXFHIyLFpKNDQc5399npG8xsLHBuxiISEckzS5cuBRJZWBdlzu7VC667Ds4/H77+dTj66LBCo4hItnV08uLsZjZvBs7PTDhgZgeZ2eVmdqqZ9TeznmZ2q5n9j5n90sx2ytR7iYhkQlIL62LO2V/6EhxxBKxYAV/5CrjHEYWIFJuOjrGeB6Snp1JgIPCHTARjZl8Cdnf39Ik2dwKPufs9ZnY88CvgzEy8n4hIJixZsgRIXmFdzDm7pASmTAkrMt53H/zpT3DWWbmOQkSKTUeHgjwB3E1Doq4H3nf3OV0NxMwOB74IHJO2rRI4C/hqtOlJ4H4zq3L3+V19TxGRrnL3xPZYU+Q5u6oKfvITqK4Obfx42HPPXEchIsWko5MX/+juz7j7s1H7J2BmNioDsfwEeBu41cweM7NDgUnASnffAODudcA84PAMvJ+ISJetWLGCDRs2sPPOO9O3b9+4w2mq6HP2+eeHxWLWr4fPfAbWro0rEhEpBh0trI9vZtv7wL1dCcLMxgDjgV+7+yXAU8BjwFBgVZPd1wGJ6xYSkeI0Z07o/B0zZgxmFnM0Oyj6nG0Gv/41jBsHs2fDccdB9AWDiEjGtauwNrNLzGw1cIWZbUtvhCT6fhfj2Bv4wN3fiB7/HDDC15cbm+zbHdja9AC1tbVMnDhxe5s8eXIXQxIRaVuqsB49enSHXjd58uTt+Qron8mYlLMb690b/vIXGDoUXnwxDAl55JGsvZ2IFKD25ux2jbF295+b2T+Ac4Bbmjy9xd2XdTrShji2x+LuG83sHcJEm35N9u0L7NDfUFFRQU1NTRfDEBHpmPQe646orq6muroaADNbmcmYlLN3tMceUFMThoX8/e9w4olw0UXwox+FwltEpDXtzdntHgoS9Uxc6e4LmrRlZjasi/G+AexsZgPTttURknGlmfUEMLNyYBjwbBffT0QkI2bOnAl0vMc625SzdzRgADzxBPzwh1BeDrffDvvtF3qxRUQyodUeazO7HHjE3d+KNv2XmZU23Q04gjBppVPcfZaZPQqcDvzCzHYBugF3AZ8CjgKmAkcD97v7ks6+l4hIprg7L730EkDq68FYKWe3raQELr8cjjkGvvhFmDkTDj4YzjsPfvCDUHyLiHRWW0NBRgC7pD3eD+gNLKfh8k0GDMpALP8B3Gxm3Qg9HGe6e52ZXQzcZGb7Ru9TnYH3EhHpsnnz5rFy5Ur69+/PiBEj4g4HlLPbbd99w9CQ7343XJJvyhS4/364+uqwoEyPHnFHKCL5yLwDy1GZ2X7Am+6+pcn2vdJ6SGIxceJE1xhrEcmlu+++my9+8YuceOKJTJ06tdPHMbNp7p7xLm/l7PaZPTssff7YY+Hx0KFwzTVwzjlQ1tHVHkSk4LWWszt6ub19mknQY4DdOxuciEi++r//+z8APvnJT8YcSYuUs9thzJhwlZCpU+HjH4dFi8L1r/feG267DdasiTtCEckX7b3c3mFm9kngGDP7ZHoD9gR+kdUoRUQSpr6+nkcffRSAE088MeZoGlPO7jizcKWQ6dPh7rth991hzhz46lehshK+9KUw8bGuLu5IRSTJ2vsl12rg94Txe4c2eW4L8JtMBiUiknQvv/wyK1asYNiwYey1115xh9OUcnYnlZSES/KdfnoYc/2rX8Ezz8BvfxvagAFwxhnwhS/AIYeE/UVEUtp7HesZZjYJOMPdlZBFpOj97ne/A+Czn/1s4lZcVM7uuvLyUDx/4QthDPbdd8Of/xx6sX/xi9AGD4aTTgrtyCM14VFEOjh5scWDmP2buz+XgXg6LUkTYUSksG3YsIHKykrWrFnDjBkzGDduXJeOl63Ji628n3J2J7jDa6+FAvvPf4aFCxue6907XMLvpJPg2GNhyJD44hSR7GotZ7d1Heu/AX3aOj5heVtd/VNEisIdd9zBmjVrOPDAA7tcVGeScnZ2mYXl0MePhxtugGnT4KGHQps+HR54IDQIEyKPOiq0I46Afk3XoxSRgtTWUJDZhLF6S2i4BmpTBnw6k0GJiCTV5s2bufHGGwG48sorY45mB8rZOWIGEyeGdu21sHhxuKrIww+HMdmzZ4d2221h3wkT4PDD4dBDw9jsgQPbfAsRyUOtDgUxs/7AOnff3OpBzHZz91WZDq4j8vFrRRHJPzfccANXXXUV48aN4/XXX6ckA7PXMjUURDk7GbZuhVdegaeeCu3558O2dLvvHgrsVKG9115Q2nSNTBFJpNZydqbGWJ8f9wSZQk7SIpIMCxYsYK+99mLDhg088cQTHH300Rk5bgxjrJWzc2jDBnjuudCefx5efBHWr2+8T58+oVd7//1DL/j++8OIEaG3W0SSpdNjrNtx4GHAV4FL0OWbRKSAbdu2jXPOOYcNGzbwuc99LmNFdS4pZ8ejV68wofHYY8PjujqYMQP+9a9QaP/rX2Ei5LPPhpay224NRXbqtrIynnMQkfbpcGFtZqXAycAFwDHA+4QxfSIiBeu6667j2WefZdCgQdxyyy1xh9NuytnJU1bWMAnykkvCtvffh5qaMIQk1Wpr4fHHQ0sZNCisDpnexo6F7t3jORcRaazdhbWZjQDOB74E7AYsAE4CHgcOzkp0IiIJ8Ic//IHrrrsOM+Ouu+5iYB7MPFPOzi8DB4aVH1OLeLqHXuz0YrumBpYvD+3//q/htWVl4SokqUJ7n33CmO3hw7WAjUiutVlYm9kZhJ6OI4G5wE+A3wHfcvfHot2ez1qEIiIxeuSRRzjvvPMA+OlPf8qRRx4Zc0StU84uDGahMB4+HD772bCtvh4WLIA33mjc3nkHZs4M7U9/ajhGjx6h4N5zz9Crveeeoe2xhxazEcmW9vRYbyQsgfs8cL67zwEws67PehQRSbB7772XL37xi9TV1fGNb3yDSy+9NO6Q2kM5u0CVlIQJjSNGwCmnNGzfsCEU1alC+8034e23YdkyeP310Jo7TqrgHjMmXKVk1KiwmqR6uUU6r83C2t2nAlOjSS8XmFkV8CBQntrHzMwzcXkREZEEcHd+9rOfcfnll1NfX89ll13GD3/4w7jDahfl7OLTq1eY2Lj//o23r14drqX99tuhzZoVbt99t6FNndr4Nd27hyI7VWin2u67h97zsi5d8kCk8LX7v4i7LwT+28zKgFOBIWZ2PXAHcCwwOTshiojkzvr167nwwgu56667gDBp8Tvf+Q6WZ9c9U86WXXaBAw8MLd3mzWH4SKrQnju3oa1YAW+9FVpTZWWhuE71mo8YAVVVDfcHDNDlAUW6dB1rMxsFfA24wN17ZSyqTiima6KKSHY8//zznHPOOcydO5fevXvz29/+ltNPPz3r75ur61grZ0tb1q0LPdlz5zbcpu4vWtT6a3v2bFxop98fMSIU+iq8pRBk7TrW7j4XuNTMpnXlOCIicVq7di3XXnstP/3pT3F39tlnH+6++27GjRsXd2gZpZwtbenbFz7xidCa2rgR5s+HefMablNt/nz44IOGYSfN2Wmnlnu7q6rCIjki+S4jo6Xc/feZOI6ISC7V19fzhz/8gSuuuIL333+fkpISrrzySq655hq6F/CFgZWzpTN69my4skhz1qxpXHA3Lb7Xrm1+MmVK//4t93YPG6YrmUh+0DQEESk67s7DDz/MNddcw6uvvgrAQQcdxC233ML+TWeAiUi77Lwz7LtvaE25w6pVO/Zyp99fuTK0V15p/vgDB8LQoTBkSLhter+yEsrLm3+tSK6osBaRolFfX8+jjz7KtddeyyvRp/egQYO48cYbOfvssynRdcZEssIs9Ej377/j1UsgXKN7+fKWC+9Fi8LqlKkVKlt6j0GDmi+6U48/9jFd2USyK1H/vMysG/AKcKm7P2NmPYGbgHVAP+AKd18bZ4wikn8++ugjfv/733PLLbcwe/ZsAAYMGMAVV1zBRRddRK9esc7jy2vK25IJJSWhx7myEg49dMfn6+pC4b1oUWiLF+94f9myhvbyy62/z5AhoaXes2nbaSdNtJTOSVRhDXwTqEp7fBvwmLvfY2bHA78CzowjMBHJPzNmzODOO+9kypQprFmzBoAhQ4bw9a9/nYsuuojevXvHHGFBUN6WrCsrayiGDz64+X22bg1FdWvF9/Ll4f7ixa2/X69eLRfdlZWh57uyUhMuZUeJKazN7BBgGfBh9LgSOAv4arTLk8D9Zlbl7vNjCVJEEm/VqlX86U9/4s4779w+fhrgkEMO4dJLL+Uzn/kM5RqImRHK25Ik5eVhkuOwYS3vs2ULLF0aiuwlS8L95tr69Q2XGmxN376NC+7Bg3dsGn5SXBLxqzaz3sDp7n6ZmV0dbZ4ErHT3DQDuXmdm84DDgfmxBCoiibRq1SoeeOAB7rvvPp588knq6uoA2GWXXTjzzDP58pe/zIQJE2KOsrAob0s+6tYtXHGkqqr1/data7noTm/r1oXVLaMRZs0yCxMvU4X2kCHNF+A77ZTJM5W4JKKwBq4AbmiybTCwqsm2dUBlcweora1l4sSGa3VXV1dTXV2dyRhFJEEWL17Mww8/zP3338/f//53tm3bBkBpaSnHHXccX/rSlzjllFPokdBrdE2ePJnJk7cvftg/zlg6qUt5WzlbkqxvXxgzJrSWuIdl41NF9pIlzbf33w9DUJYvh2mtXEG+T5+GwjvVUpMuU/f79dPY77i0N2d3aeXFTDCz44Du7v636PF84FxgInCGux+Ytu+rwN3u/qOmx9EqXiKFbevWrTz//PM8+uijPPLII8yYMWP7c6WlpRx11FGcccYZnHLKKVRUVMQYacflauXFTMlE3lbOlmKxdWsoqlsqvFNtw4a2j9WzZ/MFd/r9XXdV8Z1tWVt5MUMuB8Zbw7+CfsCDwE+i++n6AktzF5qIxMXdefvtt3n66af5+9//zpNPPsnatQ0Xl+jduzdHH300J598Mqeccgq77bZbjNEWHeVtkXYqL2+45F9L3MMCO4sXhyI7NeEyNekydX/tWnjnndBa0qNH68X3kCHhsocqvrMjCYX1WUD6EmcvAJcBLwHfNLOe7r7RzMqBYcCzMcQoIlnm7syZM4enn36ap59+mmeeeYYVK1Y02mfs2LGccMIJnHDCCfzbv/1bQa+OmHDK2yIZZAa77BLauHEt77d2bUOR3bToTt1fs6btiZfdu7fc4526X1Gh4rszYi+s3b02/bGZbQNq3X2hmU0FjgKmAkcD97v7khjCFJEM27BhA6+88govvPDC9lZb2ygdMGjQII444giOOOIIjjrqKEaOHBlTtJJOeVskHjvtBHvtFVpL1q1rXHw3V4CvXg3vvhtaS7p1a7nHO/W4oiJcG1waxF5Yt+Fi4CYz2xcYBGhmi0geqq+vZ31qKEQAACAASURBVM6cOUybNm17Ef36669vn3CYMmDAACZNmsQRRxzBpEmTGDNmDKYuk3yjvC0So759Yc89Q2vJRx+1PuRk0SL48EN4773QWtKtW+MJl831fg8YUFzFd+IKa3evSrv/IXBBfNGISEfV1dUxa9YsXn31VaZNm8arr77Ka6+9xkcffdRov5KSEsaPH8/BBx+8vY0cOVKFdB5S3hbJL336tH3Vk/XrQ/Hd0pCTxYth1aqGZedbUl7eevE9YgTstlvhDDtJXGEtIvnjww8/5M033+TNN99kxowZTJ8+nddff52NGzfusO+QIUPYb7/9OPDAAzn44IPZf//96aNly0REEql3bxg9OrSWbNjQfPGd/njlSpg/P7SW9OkTCuwRI2DkyIb7qZZPi+SqsBaRNm3atIm3336bGTNmMGPGjO2F9JIlzQ+draqqYsKECey3335MmDCB8ePHM2DAgBxHLSIi2dSrF+yxR2gt2bixofhuWoAvWBB6u9euhRkzQmtORUXjwnv33WHs2DDcZddds3NunaXCWkS2W7VqFbNnz2bWrFnbb2fNmsXcuXOpr6/fYf+ePXuy1157sc8++zBu3Dj23Xdf9ttvP3ZNWqYTEZFY9OwJo0aF1hz3MJ47NaTkvfca7s+bF3q6a2tDe/nlHV8/YEBDkZ1q++4bVruMgwprkSKzdetW5s2bt0MBPXv2bFauXNnsa0pLSxk7duz2Ajp1O3LkSEpLS3N8BiIiUijMQq/zrrvChAk7Pl9fD8uWNS6433kHZs0KbcWK0P7xj8avGzIkHG/ixNAOPTRM7Mw2FdYiBWjDhg289957vPvuu7z77rvMnTt3+/358+fvcDWOlN69ezN27FjGjh3LmDFjtt+OHj06sUuDi4hI4SopCZMfBw+Gww5r/Fx9fRhS8vbboc2aBTNnwmuvNQw3efDBsG95ORxyCJxwAvz7v8PHPpadeFVYi+Qhd2f16tXbi+X04nnu3LksXdr6QnfDhg3boXgeO3YslZWVuiqHiIjkhZISGDYstOOOa9heXw9z5kBNDUybBi+8AK+8As8+G9q3vx2K6xtvDOO3M0mFtUgCpQrn+fPnt9jSl/duqqysjBEjRrD77rszatQodt999+33R4wYod5nEREpWCUlYdz12LFw9tlh24cfwlNPwd13w9/+BnfeCU88EQruQYMy994qrEVi4O58+OGHOxTLCxYsaFfhDGHYRqpgblpADx06lLIy/fcWEREB6NcPjj8+TGrcdVe4444wVGTqVDj//My9jz55RbJgy5YtLFmyhIULF7Jo0SIWLly4Q1u3bl2rx+jTpw9VVVUttl133VXDNkRERNKkJju++y7Mnh3GXc+eHdp774XnU6qqQrGdSSqsRTrI3Vm5cuX2Arlp4bxo0SKWLVuGu7d6HBXOIiIiHeMOH3zQ+JJ86W3BAti8ufnXlpWFS/FNnBjGZJ94ImR6ZKQKa5EmNmzYsL1Ybq63edGiRWzatKnVY5SUlDB48GCGDh3KsGHDtrf0xyqcRUREGriHZdJTC8o0d7toUVhQpjWpBWXGjAnjrFPLt48alflCuikV1lJU6uvrWb58ebNDM1JFc0vXck63yy677FAopz+urKykvLw8B2ckIiKSbHV1YYGX5csbt/ffD7dLl4aiecmSlnub0/Xtu+Oy56lWVRWWSI+LCmspKOvXr9/ey7xgwYJGRfOCBQtYvHgxW7dubfUY5eXlDB06tMWieejQoey00045OiMREZHkSa2Y2Fyh3LTV1ob922PnncPiLoMHN387ZEiYfJjUL3xVWEveqK+vZ8WKFTsUy+mP29PbXFFRsb1AHj58+A49zwMHDqSkpCQHZyQiIpIM7vDRR2EVw9ra5m+bbmujn2o7szA8Y9CglltqEZg4e5szQYW1JMamTZt2GJqRXjgvWrSIzW18R9StW7dGRXKqcE7vde7Vq1eOzkhERCQ+Gza0XBQ3t62N6UM72Gmn1ovlVKuoCBMHi0GRnKYkwdq1a5k3bx7z5s1j/vz5OxTOK1asaPMYu+22W6NCuWnhrN5mEREpVJs2hQK4pR7lprfr13fs+D17woABoVVU7Hi/6TatNbYjFdaSMVu3bmXRokW89957zJs3j/fee6/R/VWrVrX6+rKyslZ7m4cNG0bv3r1zdDYiIiLZtWULrFzZ/uEXbSx/sIPu3Vsuipvbpo/YrlNhLe3m7qxatYp33323UcGcul20aBHbtm1r8fU9evRgxIgRjBw5kqqqqh0K50GDBlFaWprDMxIREcmcurpQKLd3jPLq1R07fllZ84VxS7d9+yZ3kl+hUmEtO/joo4945513mDNnDnPmzGl0/8MPP2zxdWbGkCFDGDlyJCNHjtxeRKduNUxDRETyiXvoJV6+PKzm9/77rY9Z/uCDjh2/tBT6929fkTxgQLhihgrlZFNhXaTq6up49913mT179vaiOVVEL126tMXX9e3bl1GjRjUqmFP3hw8fTvfu3XN4FiIiIh23bVvoWV62LLRU4dzc/Q0b2n9cs1Aot7dHuV8/UH9TYUlMYW1mpwE/BnYC/gh8w93rzKw/8H1gNVAOXOHu7bzAi6QK6JkzZ/LWW28xc+ZMZs6cyezZs9myZUuzr+nWrRujRo1i9OjR29see+zB6NGjGThwoFYLFBHlbEmk+vrQe7x4cUNrrmhesSIU1+3Rsyd87GOhDRzY8jjligrYbbfQCy3FKxGFtZkNA04FTgf2Bn4JLAJ+BNwLXOnuL5nZRcB1wFVxxZpU7s6iRYuYPn06b775ZqMCuqVL1A0bNow999xzhwJ62LBhGussIi1SzpY41NeHgjhVMC9atOP9JUvChMD22G23hoJ50KCW72ucsnREIgprYDhwnrvXAdPMbBxwhJk9B4xx95ei/R4CZpvZ99y9gxeRKRzbtm1j9uzZTJ8+fXt77bXX+KCFwV3Dhg1jr732Yu+9997e9txzT/r27ZvjyEWkQChnS8Zt3hwK5PnzQ1uwoOF+qmhuz4Iku+4KQ4c2rNZXWdlQLKduBw6Ebt2yez5SnBJRWLv7P5tsWgL0BY4E5qXttyQahjAReDZnAcaovr6e2bNn8+KLL/Lyyy8zffp03njjDTZu3LjDvrvtthvjx4/n4x//eKMCWstvi0gmKWdLZ2zeDAsXNhTLTYvnpUvbXva6f/+GZa1TxXP6/cGDQWuASZwSUVg3Y3/gh0A10PTix+uAyqYvqK2tZeLEidsfV1dXU11dnc0Ys+KDDz7g5Zdf5sUXX+SFF17gpZdeYs2aNTvsN3z4cMaPH8/48eP5xCc+wfjx4xkyZIjGP4vkicmTJzN58uTUw/5xxpIBRZuzpYE7rFoF77zT0N57r3Hh3JqSEhg2DKqqGrfhw0PhPHhwGO8sEof25uzEFdZmNgqodffXzcyBpl2z3YEdvgyqqKigpqYmFyFm1NKlS3n22Wd5+umn+cc//sHs2bN32Gfw4MEcfPDBHHDAAUyYMIFPfOIT7LrrrjFEKyKZkl5ImtnKmMPptGLL2QJr18KcOaGlF9HvvAOtXJGV0tJQIKeK5aYF9ODBUF6ek1MQ6bD25uxEFdZmVgpcCFwZbVoCjE573oDeQBt/9ybXsmXLthfSzzzzDHPmzGn0fI8ePZgwYQIHHXTQ9jZkyJCYohURaVkx5OxitnYtvPUWzJzZ0N56K4x3bkmfPrDHHg1t1CgYMaKhcC5LVNUhknlJ+yd+OfAjd09dxuIJ4KK056uALUDedHPU1dXxwgsv8Mgjj/DII4/wxhtvNHq+d+/eHHbYYRxxxBEcfvjhjB8/nm6aUSEi+aHgcnYxqquDWbNg+nR4/fWGInrRoub3794dRo9uKJ7T7w8cqCtoSHFLTGFtZt8iJN+eZjYSmAT8A1hmZuPc/U3geOBWd2/nxXTisXLlSqZOncqjjz7K448/3miMdK9evTjssMOYNGkSkyZNYsKECZTruy8RyTOFlLOLyebN8Oab8OqroaWK6U2bdty3WzcYOxb23rtxGzlS12oWaUkiCmsz+w5wfZPNs9x9ipl9HrjazOYDfYBv5zq+9lixYgV//etfue+++3j66afZlnbl+dGjR3PCCSdwwgkncNhhh9GjR48YIxUR6ZpCyNnFwD1cdeP55xvajBmhh7qpESNg/Hj4xCdg3LiGAlpDN0Q6JhH/Zdz9e8D3WnhuAfDl3EbUPqtXr+aee+7hf//3f3nmmWeor68HoKysjGOPPZaTTjqJT33qU+y+++4xRyoikjn5mrML3datMG1a40J62bLG+5iFXuj99gtt/PjQ+vWLJ2aRQpOIwjqfbNu2jaeeeoo777yTv/71r9tXNSwvL+e4447jjDPO4JRTTtFVO0REJKtShfQzz4T23HOwvskyPLvuCocc0tAmTAgTDEUkO1RYt9Py5cu5/fbbueOOO1icNiX6qKOO4uyzz+aUU06hn/7kFxGRLKmrg5qa1gvpMWPgsMMaCunRozWZUCSXVFi3oaamhptvvpl77rmHrdFaqiNGjODcc8/lnHPOYfjw4TFHKCIiherDD+Gxx2DqVHj00R2vEz12LEyaFNrhh4clu0UkPiqsm+HuPPHEE1x//fU899xzAJSUlPCZz3yGSy65hEmTJlFSUhJzlCIiUogWL4Z774UHHoB//QvS5sIzahQcfbQKaZGkUmGdxt15/PHHufbaa3nxxRcB2HnnnTn//PO55JJLqKqqijdAEREpSMuXw333wT33hCEeKWVlcOSR8OlPw4knhqEdIpJcKqwjL730EpdddhnPP/88AP379+fyyy/nK1/5Cn379o05OhERKTRbtsDf/ga//jU8+SREF5aiR49QRJ9xBhx/POy8c7xxikj7FX1hvWjRIq666iruuusuACoqKvjmN7/JxRdfTB9NnRYRkQybOzcU03feCbW1YVt5eSimP/95OPlkUH+OSH4q2sJ669at/PjHP+a6665j48aNdO/encsuu4yrrrpKPdQiIpJR7qFX+pZb4OGHw2MIi7FccAGcfXa4NJ6I5LeiLawvuugipkyZAsAZZ5zBTTfdpDHUIiKSFX/5C5x+erjfvTuceSZceCEceKAuhydSSIq2sL7ssst44YUXuPnmmznmmGPiDkdERArYpz8dVjo87TSoroaKirgjEpFsKNrCeu+99+bNN9/UZfNERCTruncPi7uod1qksBV1VamiWkREckVFtUjhU2UpIiIiIpIBKqxFRERERDJAhbWIiIiISAaosM6ByZMnxx1C1ujc8pPOTfJVofx+dR7JUgjnUQjnAPl/HiqscyDf/5G0RueWn3Rukq8K5fer80iWQjiPQjgHyP/zUGEtIiIiIpIB5ql1VfOcmdUCC+KOowX9gZVxB5ElOrf8pHNLluHuXlRLhnQhZ+fj77c5Oo9kKYTzKIRzgPw4jxZzdsEU1iIiIiIicdJQEBERERGRDFBhLSIiIiKSASqss8zMjjGzl8ysKu5YMsnMTjOzeWa2ysxuNrOyuGPKFDM73Mxmm9mHZvazuOPJNDPrZmavm9mkuGPJNDP7npl51F6POx7JHDPraWa3mtn/mNkvzWynuGPqrHz/XCik/F9I+b5Qcnu+5/G8/c+QD8xsINAHOCDuWDLJzIYBpwKnA3sDvwQWAT+KM65MMLOdgUOBg4BDgAfMbKq7PxlvZBn1TaAq7iAyzcx6Eya9HBNtSupkZumc24DH3P0eMzse+BVwZswxdVi+fy4UUv4vwHyf97m9EPK4Cusscvf3zezBuOPIguHAee5eB0wzs3HAEeRhYm3GFuAHHmb1PmxmbwDbYo4pY8zsEGAZ8GHcsWTBBcAs4Dl33xR3MJI5ZlYJnAV8Ndr0JHC/mVW5+/zYAuuEAvhcKKT8XzD5voBye97ncQ0FyTJ3r487hkxz939GSTVlCbAwrngyyd03Rkk29Zfz28AzsQaVIdH5nO7uU+KOJdPMzIAvAD8G3jezvOvJlFZNAla6+waAKP/MAw6PM6jOyufPhULK/4WS7wsltxdKHldhLZmwP3B73EFkkpkdBjwKlAO9Yg4nU64Abog7iGzw4CDCV4i3AneZ2adiDksyZzCwqsm2dUBlDLFIY3mf/wsg3xdEbi+UPK7CWrrEzEYBte6edxMM2jAX+C1wFPDDeEPpOjM7Dqhx9xVxx5JN7v6hu38H+AHw9bjjkYxxYGOTbd2BrTHEIpECyv95m+8LMbfnex7XGGvpNDMrBS4Erow7lkxz92XAFDOrJ0wIyXeXA+PDN20A9AMeNLPvu/tN8YWVNbeQh1/pSouWEP7NpusLLI0hFqGw8n+e5/tCzu15mcdVWEtXXA78yN03xx1IFtUQPtTz3VmEHr6UF4DLgMfjCSfr6oFX4w5CMuZp4Ndm1tPdN5pZOTAMeDbmuIpZIeb/fMz3hZzb8zKPayhIllnDn5HW6o55xsy+BUwHeprZSDM7L/paMK+ZWS8zG5+26QTg5rjiyRR3r3X3xalGmPle6+5r444tE8ys0szOMrOS6P/cN4Bvxx2XZIa7LwemEr6qBzgauN/d860IAvL/c6FQ8n8h5PtCyu2FksfVY51FZtYXODt6eK6Z/dzda+OMKRPM7DvA9U02z8r3GcmRcYTLLs0CXgRedveHY45J2taPMHnnW8BzwM/y7TJs0qaLgZvMbF9gEFAdczydku+fCwWW/5Xvk6Ug8rhFV5oREREREZEu0FAQEREREZEMUGEtIiIiIpIBKqxFRERERDJAhbWIiIiISAaosBYRERERyQAV1iIiIiIiGaDCWqQFZnaUmT1vZpPa2G+cmf2vmf13Bt/7h5k6lohIHMzs42Z2byZzYwvvM9XMvpDN98gWMzvazF5o63NG8ocWiJHEM7M9ge8DHwD7Ax8H/t3d/5jlt14E7NuO/RYSLmxfmok3jZZr3pqJY4mIAJjZ/sAPgE8C/wfMIyy0Uw780t2zsQT2XGBnMpQbW3Eb8GY2Dhzl49OA24E5hCXDISwj/lnggC4uYtLezxnJEyqsJR88AHzO3V83sxLgzly8qbvPMbNV7dhvrZkty+BbHw78I4PHE5Ei5+6vmNmfgbHu/unUdjM7EfirmV3u7rdl+D03ZDg3AmBmJwEzUgWtuz+S6fdIcfetwD1m9jXgGXf/TlocD3bmmGY2ADjC3e9x99lmtjJD4UoCaCiIJJqZ7QqMBtYDuHs9cDWQqyVD6zO8X3scATydweOJiADUNd0QLeH9XeCnZjYoC++Z0VxtZkMIPdS5tsPPDngSWNKRg5hZd+D3QM+0zZn8/JCYqbCWRHP3D4DZwGNmdlC0bQHwspntYmbXm9nLZnaymS03s7fNbL/U683sTDP7tpk9Y2b/k7Z9opl918zuMbO/mNnO0fYeZvYLM7vazO4gfI3ZIWb2X1FMnzOz/c3sgeh410cx3mNmFdH2NWZ2TZNDdHP3zdGxvmFml5nZv8zshg7/AEVE2nYHYUjIGQBmtruZXWdmvzGzp6JiFjPrHuWyb5vZP83s86kDmNlXzOy/zexFM7ukpTdqLieb2elRzjzXzN4xsxvMrNLM7ozy6T/N7MjoEGcBQ4ArzOzUaC7Mk2Z2Ttp7XGBm3zez35vZ78yst5ntbGbXmFmNmR1pZrPN7D0zq+rMDyz6mUyKerQxs9PM7GYz+7GZPWZmfaJzuNnMHoo+b5ZFP+O9gDPM7LK0Q/aLPovWmdkFLbxniZndaGbV0Xu4mf019fuRhHB3NbVEN2BvYAHhr/opwIBoeykhSX1AGOu2C/AEMJPwR+PBwHeiffsRer1PBHYC/px2/H8BP4zu/xSoju4PBLYRkmdbMf6W0OszgDAWr1e0vTvwOPAgUEn4QNgKXAP0AiZEj1P7p8aPp877tuh+X+DKuH8Xampq+duAc4HFLTz3PnBLlFf/ApRG2+8G7o3u/wQ4Lbp/HlAb3T8TODu6PzrK1ftEj38LfDe631JO3hVYDfwMOAA4GvgRcHW07zeBh9JidaAqut8TeBs4N3p8HHB32r5TgV8Rhr5+Jvq8OC06z8eA77fzZ/dM9FlxA/BDwnjro9OeXwkMi+7PJHwmdQe+DSwGDgS+HOX9Z1LxRvvPB34cPXcmsLCFGM4Hpkb3B0Y/52Pj/nel1rhpjLUknrvPNLN9gO8BXwFOMrOjPYy5rgXWuvv9AGb2XeA5YCThQ6TezL4eHepxQgI/Edgtbfvb0Wv7ABdHr8Xd3zezRR0ItT9wK/Bld98QHWOzmS0H5rn70uh9lgPPehh/OJ2Q8PsTJkGeSPgQANgA/IeZzQF+SfijQkQkG5xQbB4ADAe+ZmYAq4CdLTy4kNApAPA7wlAICLm2pkmu/Rgwo8l7nEszOdndPzCzNcCD7v4ygJktBNab2ccIPbx9mg3afaOZrUjbdEl03JTfAP8bbV9N+Lz4S/QeNVGc7fW0R2OszeweQmdOyhnuvtDMDiUU1H2i/L8UeNfdXwJeil7b3LEfij4TXmklpv2AjdF5v29mH9B4SIkkgAprSTQzKwWGepik8v+i4RmPAJMJPQBNzYxudwGGAX9w97ujbT+LjnkFodD9WZP3mkBIiOs6Ge5o4N+A3YHX07Y3HWO4bfsT7vVRkk0Ny9rNw/AX3H2emZ0L/Bz4T+CLwApERDLIzHYhfNs2k5A3VzWTHwcSelS7Abj7NkJnANFrrnX356PHjV6bptmcHHEa58rFwBXAO8CLwIhWTiH9dXsAf097/C5hmMsAdszFdXT+iiWvEv4ASVltZjcC9xJiT1XPTc+rLXW0XJs9CvzSzHoQfhcbCb3fkiAaYy1Jt33cH4C7vw58Cxjbwv7dCElsLrAMOD39STM7MNp+opn1bLJ9bfRwz07G+jxwM3BXlPg6xMwqSZsIY2aDCVdEGU2YzHh3Cy8VEemKcwjfkN1LyI+HRj3FwPb8uBLYAhyftn0fM+tLk1xrZmVRR0VTLeXk5twOzPdwWdXNHTiXhYQe7u1vAawhDHXJGHevjzo/9o7+MHkEuM7dazL5Pk3e8yHgIcI3q2cCh7v7mmy9n3SOCmvJB1+Lis6UwYSxcSm7mlm36P4kwvi61cCfgM+Y2c/NbFLUm7CekAB7A1PN7LhooshYQjE+HbjezLqZWX/C2OYBZlZqZlea2dktxGhRu5rw4fOTZp5run/Tx58mjAdM2QM41d3XAl9r5jUiIh2xQ0+omX0SuBa40N1rCb3D7wOPmNlJ0QTFw6Ie6j8Trh5yqpkdBZzj7usIufb/mdm1ZnY4oYMhNYwuPf+1lJNT0nuPxwMVZtYLOBToaWYjo+e2Eib7pTpY0t/jduCzZtYvenwg8Gt3ryPUPM3mYjM7KIqrews/uxKa1EzRHx8nEr6l7EcYYjgGGAr0SpsY2bRXfEsz8TcOqpnxItHvYifgKUKvfEf+4JAc0VAQyQdDgZlm9jDhazIDLkp7vg74tpmtJvTufg3A3Z8ws/9H+DrxVMLkvzcBzOwU4BeED4opwDfd3c3sc4RLIc0E7iF8pfcJwoIKnyMMMWm0MI2ZHUCYlDOM8NXgA8C1UTyPAwcBVWa2O6EnZRDwOTObRUPvzblAb3ef3OTc77SwsEMdYbKQiEiHRXnkLELx92ugljABbmfgU+7+AoC7b4ny468Jue5BwthqgEsJnRK/J3xDl7oSxx2E3HcxYcjaxe6+IiocDwJ2N7M/tpSTzex0wrji881slrsvIUyk/BGhqP49cAqwD/Ae8AdC7/p/WLha1N7ACWb2iLv/xcz2AO43s38RvvX8dtSzfgYwKDq/mYTFcj5mZntHcX6BMG78gbSfWzkh9+8LDI4KbyeMbT4tOuZrwLPAK4ThLY8ScvpUwoTJfc3sVHdPHfcewmI9H0WF+CDgdDObR+iJJnp90zUbFhImZ54K9ABKzOxJwgTGXF2CVtpg+l1IPrOwDOxv3b0qR+/3TXfXcuMiIgUmGpbS3d0TuUCXmX2TcIWW+dHjHoRvG25Mzc2R+GkoiOS75oZZZOeNwlemr+bivUREJHeiyZn7Jbio7gF8B6hK29wfWK+iOlk0FETyVjS+7WzCV3tnuPu9WX7Ld6OvKEVEpLCUufsv4w6iJe6+ycKy6rdEE9uXEoatfKf1V0quaSiIiIiIiEgGaCiIiIiIiEgGqLAWEREREckAFdYiIiIiIhlQMJMX+/fv71VVVR16TW1tLRUVFdkJKIGK7Xyh+M652M4XCuOcp02bttLd8/skOqgzORsK4/cdB/3cOk8/u84r1J9dazm7YArrqqoqamo6tpLoxIkTO/yafFZs5wvFd87Fdr5QGOdsZgvijiHXOpOzoTB+33HQz63z9LPrvEL92bWWszUUREREREQkA1RYi4iIiIhkQFEX1tXV1XGHkFPFdr5QfOdcbOcLxXnOxUy/787Rz63z9LPrvGL82RXMAjETJ070QhzHIyKFz8ymufvEuOPIJeVsEclXreXsou6xFhERERHJFBXWIiIiIiIZULyF9ZQpsN9+cOutcUciIiJtWDd9HTUTapj1pVlxhyIi0qLiLaw/+ACmT4f58+OORERE2lC/oZ6PXv2IDbM3xB2KiEiLirew3mmncLtuXbxxiIhIm6zMAPBthTHhXkQKU/EW1n37htu1a+ONQ0RE2lYabrxOhbWIJFfxFtbqsRYRyRtWqh5rEUm+nBfWZnaMmb1kZlVp204zs3lmtsrMbjazsrTn+pvZr8zsRjP7iZmVZyQQ9ViLiLQpKTk7VVizLRNHExHJjpwW1mY2EOgDHJC2bRhwKnA68J/A+cDX0152LzDF3a8A5gDXZSQY9ViLiLQqSTlbY6xFJB/ktLB29/eBB5tsHg6c5+7T3P33wC+AIwDM7CBgjLu/FO37EPA1M+vd5WDUYy0i0qok5WwNBRGRfJDzoSDuXt/k8T/dvS5t0xJgYXT/SGBe2r5LortdX/pXPdYiIm1KTM7W5EURyQNJadZzywAAIABJREFUnLy4P3B7dH8wsKrJ8+uAyqYvqq2tZeLEidvb5MmTW38X9ViLSIwmT568PV8B/eOOpwtykrPVYy0icWpvzi5r6Yk4mNkooNbdX482ObCxyW7dga1NX1tRUUFNTU3736x7d+jWDbZsgU2boEePTkYtItJx1dXVVFdXA2BmK2MOp1NymbNTY6w1eVFE4tDenJ2YHmszKwUuBK5M27wE6Je2jwG9gaUZeEPoH/3BUVvb5cOJiBSTXOds9ViLSD5ITGENXA78yN03p217AhiT9rgK2AJ0oGu6FRUV4VaFtYhIR+U0Z28vrDXGWkQSLI7rWFvqbtq2bwHTgZ5mNtLMzjOzUe5eAywzs3HRrscDt7r7lowEM2BAuFVhLSLSrMTk7NTkRfVYi0iC5XSMtZn1Bc6OHp5rZj8nfJV4fZNdZ7n7lOj+54GrzWw+4Xqq385YQKke6xUrMnZIEZFCkaScretYi0g+yGlh7e7rgF9GLeV7UWvpNQuAL2clIA0FERFpUZJytlZeFJF8kKQx1rmXGgqiHmsRkUTT5EURyQfFXVgPHBhuly2LNw4REWnV9qEgmrwoIglW3IX1sGHhduHC1vcTEZFYqcdaRPJBcRfWw4eH2wUL4o1DRERal/q0qgd3FdcikkzFXVgPHRpuFy2CbZoRIyKSVGbWqLgWEUmi4i6se/YMExjr6jTOWkQk4TTOWkSSrrgLa9BwEBGRPGHlobCu36ouaxFJJhXWI0eG27lz441DRERaVdIjfGTVb1JhLSLJpMJ6r73C7VtvxRuHiIi0qqR7+MjyzRoKIiLJpMJ6773D7cyZ8cYhIiKtUo+1iCSdCmv1WIuI5IVUj7UKaxFJKhXWo0ZBWRnMnw/r18cdjYj8//buO06q6v7/+OuzfbYvywIC0gQBRQEFo7FhwRhN0BhNvjHF8lVSvtHEmMQaNUaNv0TT1JgQo6bHGE3sRoy9oWBFaUoHwV1g2WV7Ob8/zuzubG8zO+39fDzuY2bunbn3c2fgM589c+45It1obbGuU2EtIrFJhXV6OkydCs7B8uXRjkZERLqhriAiEutUWAMccoi/XbIkunGIiEi31BVERGKdCmuAj33M36qwFhGJWeoKIiKxToU1tBXWr7wS3ThERKRblhmcIEYt1iISo1RYA8yYAdnZsHYtbN8e7WhERKQL6mMtIrFOhTX4UUEOP9zff+qp6MYiIiJdau1jra4gIhKjVFi3mD/f3y5eHN04RESkS6nZqQA0V6mwFpHYpMK6xQkn+NvFi/3QeyIiElNS831h3VjZGOVIRES6NuSFtZnNN7MlZjYhZF3AzG4xsxvM7HYzy+/LtrA64AAYMQI2b9YsjCIiQbGUs1PzfGHdVNEUrl2KiITVkBbWZjYSyAUO6bDp18ALzrnLgQeA3/ZxW/ikpMCnPuXv33dfRA4hIhJPYi1np+WnAdBUqcJaRGLTkBbWzrnt+ETbysxGA2cCDwVXPQksMLMJPW2LSICnn+5v7703IrsXEYknsZazW7uCVKgriIjEpiHvCuKc63jVyTygzDlXHdzeCKwDju5lW/gddxwUFvqpzVeujMghRETiSSzl7LS8YIu1uoKISIyKhYsXxwA7OqyrBEb3sq2d0tJS5syZ07osWrSo/5FkZMCpp/r7f/lL/18vItJHixYtas1XwPBox9MPUcvZunhRRKKlrzk7behC6pYDajqsywQagvd72taqpKSEpUuXDj6as86Cu++GO++Eq6/2Y1yLiITZwoULWbhwIQBmVhblcPojajm7tY+1WqxFZIj1NWfHQov1FqCow7o8YGsv2yLj6KNh6lTYuhUefjhihxERiVNRy9kto4Koj7WIxKpYKKyfBkabWQDAzNKBccCzvWyLDDMI/kXC7bdH7DAiInEqajk7bZhvsW7cqcJaRGJTNMaxtpa7AM65bcDDwHHB9ccD9znntvS0LaJBnnUWZGfDE0/A669H9FAiIrEslnJ2enE6AA1lDbgmTeQlIrFnqMexzgO+Fnx4tpmVBO9/HTjFzK4ATgIWhrysp22RUVwMX/+6v3/ddRE/nIhILIq1nJ2SnuJbrR007OzUbVtEJOrMJcj03XPmzHFhuXixxbZtMHEi1NbCa6+BvwpURCTszGyZcy6pksxAc/aSaUuoWVXD3HfnkrNfTgQiExHpWU85Oxb6WMemUaPgwgv9/YsuggT5A0REJJ5ljMgAoP6j+ihHIiLSmQrrnlxxBYwYAS+8AP/8Z7SjERFJeukjgv2sS9UVRERijwrrnuTnt/WxvugiKC+PbjwiIkkuoyTYYr1dLdYiEntUWPfm3HPh0ENhyxb4zneiHY2ISFLLGO0L67rNdVGORESkMxXWvUlNhbvugsxMf/voo9GOSEQkaWWNzwKgboMKaxGJPSqs+2LaNPjRj/z9s87yrdciIjLkWgrr2g21UY5ERKQzFdZ99Z3vwPz5UFYGn/88NOjCGRGRoZY5PhNQYS0isUmFdV+lpsJf/gJjxsCLL8Ill0Q7IhGRpJM5JhNSoP7Deprrm6MdjohIOyqs+6OkBO65B9LS4Oc/h9tui3ZEIiJJJSU9hcyxmeDUai0isUeFdX8dfjjccYe/f+GF8MAD0Y1HRCTJZE/NBqB6RXWUIxERaS8shbWZnRSO/cSNs86Ca6+F5mb4n/+BJ56IdkQiIn0W7zk7Z38/lXnVu1VRjkREpL20/jzZzE4ErgRG0VaUpwAjgUB4Q4txV14JH34It98Op5ziW65POCHaUYmItErUnJ29f7DF+j21WItIbOlXYQ38Cfgx8AbQctWIAWeGM6i4YAa33grOwW9+AwsWwP33w0lx3RAkIoklIXN2zn5qsRaR2NTfwnqtc+5nHVea2cowxRNfUlL8BYxmvuV6wQLf//rss6MdmYgIJGjObukKUr2imuaGZlLSdbmQiMSG/majr5vZ6V2sXxCOYOJSS3F92WXQ1ATnnAPXX+9bskVEoishc3ZaQRqBKQGaa5upeket1iISO/pbWB8O3GFmTSFLM3B7BGKLH2Zwww1wyy3+/pVXwnnnQa2GghKRqErYnJ3/sXwAKpZURDkSEZE2/S2srwUWAlOASSHLNeENK05985tw772QlQV33glHHQWbNkU7KhFJXgmbs/MPDRbWr6iwFpHY0d8+1o8BDzjn6kJXmtnvwhdSnPvsZ2HSJDjtNHjtNTjoIPjzn+ETn4h2ZCKSfBI2Z+d9LA9QYS0isaW/LdZ/By4zs6NClqOByyIQW/yaPRuWLvXD75WVwYknwgUXQE1NtCMTkeSSsDk7d2YuqXmp1KyuoXaTut2JSGzob2H9A+DbwB+BP4QsXw1zXPGvuBgefRR+/GM/Bfqtt/rW69dfj3ZkIpI8EjZnp6SnUHhMIQC7Fu+KcjQiIl5/C+ufAns75yY45yYGlwnAaYMNxMymm9ltZnaRmf3azGYF1wfM7BYzu8HMbjez/MEea8ikpsKll8KSJTBtGqxcCR/7GFxzDdTXRzs6EUl8CZ2zh50wDICdT+yM1CFERPqlv4X1JudcZRfrR4Uhlj8BNzjnfg7cCPwtuP7XwAvOucuBB4DfhuFYQ+ugg2DZMt8dpLERfvhDX2C/9Va0IxORxJbQObtofhHgW6ybG5t7ebaISOT1t7D+RscVZnYG8IswxDIdyAverwEKzGw0foawh4LrnwQWmNmEMBxvaGVnw69+BU8/DRMmwJtvwty58KMfQUNDtKMTkcSU0Dk7MCVAYN8AjTsbKX+mPBKHEBHpl/4W1sPM7AcAZpZtZncAdwCvhSGWf+LHW80DvgRcAMwDypxz1QDOuUZgHXB0GI4XHfPmwTvvwNe/7gvqq66CQw+F5cujHZmIJJ6EztlmxojPjQCg9B+lkTiEiEi/9Kuwds6dDCw2s1uAN4F9gZlAVzN79dc3gFp8wt/jnLsPGAPs6PC8SmB0GI4XPbm58Otfw5NPwrhx/oLGgw/2Fzo2NkY7OhFJEMmQs0vOKAGg7F9l6g4iIlHX4zjWZnZUN5s2AGuAm4DxwKnARYOMJYDvozcG+IWZrQMc/ifGUJlAp74TpaWlzJkzp/XxwoULWbhw4SBDirDjjvOt19/7HixaBJdfDv/6F9x9N+y3X7SjE5EIWrRoEYsWLWp5ODwc+0zGnJ1zQA6BqQFqVtWw64ldFJ9U3O99iIj0pq8525xz3e7EzF7B96MrB7prCjBgL+dc5sBCbT3W48D/OOfKgz9dXogfJupq59y+Ic9bE1z319DXz5kzxy1dunQwIUTXE0/4adA3bYL0dD+L4w9+AEVF0Y5MRCLMzJY55+b0/sxe95OUOXvDjzew7vJ1DD91ODP+NSMs+xQR6U5PObu3riDXAlOcc+NDhmrquEwAvjjIAIcDM51zLVefXAdU4FtZRptZIPi8dGAc8OxgjheTTjjBt16ff77vDvLzn8Pkyf6CR13cKCJ9k5Q5e9Q5o7A0o+yhMuq21vX+AhGRCOmxsHbOPeqc+6i3nTjn/jnIOHYCtWY2JmTdDuBt4GHguOC644H7nHNbBnm82FRQ4LuELF3qL3LcuRO+9S2YMQP+8Q9oVv9BEelesubszFGZFC8ohib48PcfRvJQIiI96u+oIBHhnGvG9/m7yswWAlcC33POVQBfB04xsyuAk4AY7zgdBgcdBE89BQ88AFOmwOrV8PnPw8yZcN99KrBFJKpiMWeP+Yav8bfcuoWmmqahOKSISCc99rGOJ3Hfx7o7DQ1w111w3XW+/zX4Avvqq+HUU8EsuvGJyKCFq491PAl3znbOsWzuMvYs28OUX09hzNfH9P4iEZEBGEwfa4m29HRYuBDWrPFD9I0Z42dsPO0037L9739DgvxxJCIyUGbGuEvHAbDpp5tobtAveyIy9AZdWJvZbDP7s5ldHJczIsaLzEw/qcz778Mtt8Do0X72xs98xo+B/fDDKrBFpFeJnLNLPlNCYN8Atetq1ddaRKIiHC3WXwOOd87dDHy85WpwiZCsLD8U3wcf+BFD9toL3ngDPv1pOOwwWLxYBbaI9CRhc7alGhOvnwjA+qvX01ipCbdEZGiFo7D+JnBw8P7fgfQw7FN6k5UFF1zgC+yf/QxKSmDJEj9s3+GH+y4iushRRDpL6Jxd8tkS8g/Np+GjBjb9dFO0wxGRJNPvwtrMPmlmXwzePwQ/ZuoW8FeKB68Kl6ESCMBFF8HatX5K9GHD4OWXfReR/faDO+6A2tpoRykiUZJsOdvM2OemfQDYdNMmatZ2nAhSRCRy+lVYm9n1wL3AZwCcc68CXzGzMyMQm/RHbi5ceils2AC/+AWMGwerVvkJZyZM8EX3rl3RjlJEhlCy5uyCwwsY8YURNNc0s/obq0mU0a9EJPb1t8X648BY4PmQdY8BN4YtIhmc3Fw/qcz778Nf/gKzZsH27XD55b7Y/ta3fOu2iCSDpM3Zk38+mbSiNHb9Zxcf/a3XOXNERMKiv4X14uAUtqF//p8exngkXNLT4cwz4fXX4Ykn4PjjYc8ef8HjlClw+unw0ku60FEksSVtzs4YmcE+P/VdQtZcuIa6LZrqXEQir7+F9VozuxyYamb/Y2b3Af8H/L/whyZhYQbz5/vRQt58E846C1JT/QyOhx/uRxL5xz+gUVfPiySgpM7Zo84dRdEnimjc0ciKr6zANashQUQiq1+FtXPu78BTQCrwFaAcOME5d1sEYpNwmzkT7r7b98O+4gp/oeOSJX669AMPhCefjHaEIhJGyZ6zzYxpd08jfUQ65U+Vs+kmjRIiIpHV71FBnHOvOOe+5pw7yTn3v8DTZnZcBGKTSNlrLz9F+saNfjbHiRNhxQrfsj1vHjzwADQ1RTtKEQmDZM/ZmaMymXbXNADWXbGO8ufKoxyRiCQy6+lqaTN7EMjtZR/DgTLn3LHhDKy/5syZ45YuXRrNEOJXXR38/Odwww1QWenXTZoEX/0qnHOOHyNbRCLGzJY55+aEYT/K2d344PsfsOmnm0gvSefgZQeTtXfWkB1bRBJLTzm7txbrt4FXgT8BfwAygHeD91uW+4HlYYtWhl5mph+qb/NmX2BPnOhHDrnkEhg7Fr74RXj+eV3oKBL7lLO7MfGGiRTNL6KhtIHlpy6nqUa/yolI+PXWYl0QfE558PG1zrmrunjeA865UyIXZu/UYh1GTU3w+ONw++3w6KNtBfV++8HXvgZf/jIUFkY3RpEEEsYWa+XsHjTsbGDZ3GXUrq1l+GeHs/89+2OpNqQxiEj8G3CLtXNud0uCDhppZqkddj4POGTQUUrsSE2Fk0+Ghx+Gdev8hY4jR8J778GFF8Lo0fC//wuvvaZWbJEYopzds/Rh6cx4YAapBamU3VfGmm+u0eQxIhJW/b148Z/AW2b2GzP7qZk9DCwGfhX+0CQmjB/vL3TctAnuvReOOw5qauDOO+GQQ2DOHLjrLk2bLhKblLM7yJ2RywEPHoBlGlt/s5X1P1wf7ZBEJIH0d7i9xcAJQBkwFT900+nOuR9HIDaJJenpflKZJ5/0U6VffLEfru/11+Hcc30B/qMfwe7d0Y5URIKUs7tWeFQh+/19P0iBDT/cwMafbIx2SCKSIHrsY92nHZilALOdc8vCE9LAqI91FNTWwj33wC9+4SefASgqgu98x3cZyc+PbnwicSJcfaz7eCzl7KAP7/yQVeetAgcTfzyR8ZeOj2o8IhIfesrZaf3c0Z1drC4BmoBTBxCbxLOsLD+T41e+Ak8/DT/8ITz3HPzgB/Czn/kC++tfh+LiaEcqkpSUs3u217l7AbDqvFWsu2wdOBh/mYprERm4/vaxngJsBDaELJlAVFs+JMrM4Nhj4Zln4Kmn4MgjYdcuX2CPHeu7iizTPxGRKFDO7sVe5+7F1DungsG6y9ex9rK1uqBRRAasXy3WwHnOuVWhK8xsGnB2uAIys0OBI4D3gReAKuAnQCVQBFzinKsI1/EkjMzgmGP87I1PPQU33wyPPeYvbrzrLjj0UPjmN31f7czMaEcrkgyUs/tgr7P3wtKMlWevZOONG6nbWsfUO6aSkt7vyYlFJMn19+LFVV2srgPOC0cwZnYO8Cnn3E3OuX8758qAXwMvOOcuBx4AfhuOY0kEmfnRQx59FNasgYsu8uNev/IKfOlLsPfecOWVfqQREYkY5ey+G/WlURzw0AGkZKew/Y/beefT79C4pzHaYYlInOlXYW1m68xsbciyAVgB/GuwgZjZ0cAXgR+ErBsNnAk8FFz1JLDAzCYM9ngyRCZP9v2tN2+GRYvgwAOhtBSuv97P8Hj66b5/tn56FQk75ez+Kf5kMbOemUV6STq7/rOLNw5/g5p1NdEOS0TiSH9/51oMnAucE1y+BBzonDs/DLH8DJ/wbzGzx83scGAeUOacqwZwzjUC64Cjw3A8GUo5OXD++X70kOeeg89/3rds33ef7589Ywbceits2RLtSEUSiXJ2P+XPzWf2i7MJ7Bug6u0qls1dxq6nd0U7LBGJE/0trP/snHvGOfdscHkeMDObPJggzGwqMBv4nXPum8B/gceBvYEdHZ5eCYwezPEkisz8xY1//zts2ABXXw2jRvlZHS+4wF/sePDBcNll8J//wJ490Y5YJJ4pZw9A9pRsDlpyEMM+OYzGHY28Nf8tNt+yWRc1ikiv+ltYn9jFuu3AvYOMY39gp3Pu7eDjWwEDHNDxd7hMoKHjDkpLS5kzZ07rsmjRokGGJBE3ejRcc40vsP/2N1iwAAIBP+nMjTfCiSf6cbE//nG4/HJYvBiqqqIdtUhYLFq0qDVfAcMjdBjl7AFKL0zngIcOYO/v7w1N8P6F7/Pu6e/SsLPTqYhIEuhrzu7TBDFm9k3gOiCvm6csds51lcD7xMw+ByxyzhWGrHsD+AdwjnNu35D1a4CrnXN/Dd1HLEw2IGFQU+OH7Xv6aX+7bBk0N7dtT0vzU6nPmwdHHeWL7rzu/lmKxIdwTxCjnB1eH93zEavOX0VTZROZe2cy/a/TKTyisPcXikhC6iln93nmRTM7EDgL+FWHTfXOuQ8HGeA0fF+9Uc657cF1r+FbQW4DSpxzNWaWDuwBJjnn2nXGjackLf2weze88IIvsp95xrdmhxbaKSkwa5Yvso88Eo44AkaMiFa0IgMSiZkXlbPDq2ZtDe994T0qX62EFNj74r2Z8MMJpAZSox2aiAyxsBTWwR2lO+c6/Q5mZuOccxsHESNm9ijwiHPuNjMrBJ4FDgb+jO8n+LCZfRL4snPuzI6vj7ckLQNUXt5WaD//vC+0GzsMiTVtmi+yjzzSt2hPmuT7dovEqEhNaa6cHV7NDc2sv2o9G3+yEZohMCXA1N9PpfBItV6LJJMBF9Zm9l3gUefce8HHVwAd/zw34Bjn3LxBBjkc+CWwFBiHvyjmPTMrwk82sB4YBVzmnOt0RVs8JmkJg6oqPz7288/75eWXfXeSUMOH+8lpDjvM386dq+4jElPCVVgrZw+NiiUVrDx3JdXvVQMw+mujmXj9RNKHpUc5MhEZCoMprG8D/uKceyn4+D4gB9iGv0gFfJI+1Dk3LaxR91M8J2kJo/p634rdUmi/8oofNzuUmR/e79BD25Zp03y3EpEoCGNhrZw9RJrrmtlw/QY2/ngjrtGRVpzGxOsmMvr80ViqfiETSWTh7ApyELDcOVffYf1+LS0k0RLvSVoixDlYv963ZL/yil/eeKNz95GCAn9RZEuhPXOmH7VEXUhkCESwK4hydoTtWb6H9y98n/KnywHImZnDlF9OofBodQ8RSVQ95ey0fu7rAOfc6x12PhXYB4hqkhbpkpmf4XHiRDgz2M2zpsYX1y2F9iuv+OnVFy/2S4uCAt+yPWMG7L9/2/2Skuici0j/KWdHWO6MXGb+dyal95XywcUfUPVWFW/Oe5NhnxzGxOsnkjdb3c5EkkmfCmszOxL/8+F8M1vXYfMw/FXnD3V6oUgsCgT8RY0f/3jbui1bYMkSX2S/+iosXw47dsCLL/olVEmJL7SnTYPp09uWMWPUwi0xQTl7aJkZI04fQfFJxWy6aRObfrqJnY/tZOdjOyk5o4QJ104gZ1pOtMMUkSHQ13GsDwD+CEwEOs7tWo/v03dt+MPru0T5WVFihHPw0Ue+wG5Z3n3X31ZWdv2a3NzOxfb06X5UknRd1CTdi8A41srZUVRfVs/GGzey5dYtuDoHBsNPG864S8aRPzc/2uGJyCCFaxzrAuAM59wd4QwuXBI5SUsMcc53G3nvPVixAlau9LcrVkBZWdevSU+HyZM7F9xTp0KOWrEkYuNYK2dHWe3mWjZct4Ftd23D1fvv2sLjChl36TiKjivC9AuXSFwK28WLPRzgCOfcC4Pe0SAkQ5KWGFdW1rnYXrHCT9nenXHjfJE9ZQrsu2/b7fjxkKqJJ5JFpC5e7OF4ytlDqG5rHZt/sZmtv9lKU2UTADkH5DDm/8Yw4osjSMvt7+VOIhJNgxlu70Egt7f9A/s756I63V0yJWmJM1VVsHp1+2J7xQpYswYaOs3d4aWnwz77dC64p0zxo5VoaMCEEsbh9pSzY1hDeQNbb9/K5l9upmG7/7+fmp/KqLNHMeYbY8iemh3lCEWkLwYzKsgqoBzYQtsYqJ32D3xq4OGJJLicHJg92y+hGhth7Vrfwr1mjS++W263bPHrV67svL9AwBfYXRXdJSW6gDK5KWfHsPTCdMZfNp69L96b0vtK2XLrFipeqmDLr7aw5VdbKDi6gFFnj6Lk9BK1YovEqd5arIcDlc65uh53YlbsnNsR7uD6IxlbPySBVVXB++93LrjXrOk84U2oggLfn7ul8A5diouHLn7plzC2WCtnx5nKNyvZettWtv9lO801zQCk5KRQcnoJo84eReFRhViK/lgWiSVD0cf6vGhfIKMkLUmjvLzrgnv1aqio6P51RUVdF9xTpkChJrOIpij0sVbOjjGNFY2U3lvKtru3sfuF3a3rM8dnMuLzIxjxuRHkHpSrCx5FYkDECmszGwf8H/BN51xUhzdQkpak55xvzV6zpuulqqr71w4f3n3RnacJLiJtqApr5ez4UP1+Ndv/sJ1tf9hG3aa2Hx+yJmUx4nMjKPlcCbmzVGSLREtYC2szSwUWAOcD84Htwf2MGWygg6EkLdID52Dbtq4L7vff97NRdmfkyK4L7smTNVxgmESysFbOjl+u2bH7xd2U/qOU0n+WUr+tbWb6wOQAwz8znOIFxRQcVoClqsgWGSrhGsd6InAecA5QDGwALgT+AxzmnHspPOEOjJK0yAA1N8PWrW1Fdseiu66H7rqjR7cvtEPvBwJDdw5xLkLjWCtnJxDX5Nj9wm4+uucjSu8rpeGjthGF0orTKD6pmOIFxQw7YRhp+brwUSSSBlVYm9kZ+JaOY4H3gTuAPwCXO+cuCnOsA6YkLRIBzc2weXPXLd0ffND9cIEAY8d23dK9zz6QmTl05xAHwllYK2cnvubGZnY/v5sdD+2g7MEyaj+obd1m6UbhvEKGnTSMYScMI3t6trqMiITZYIbbA6jBT4H7EnCec251cKeDv+pRRGJbSoqfxGbcODjuuPbbmppg48aui+5163xBvnkzPP10+9eZ+f11VXRPnAgZGUN3folJOTvBpaSlUHRMEUXHFLHPzftQvbK6tciueLmCXYt3sWvxLj7gAzLGZDBs/jCKTiii6PgiMkr0/0skkvrTFWQcvhVkAvAAcLRz7oLgNnPhGF5kENT6IRJDGhr8jJNdFd3r1/uW8K6kpvpZJ7squidMgLTE/Ik7Ql1BlLOTUH1pPTsf28nOJ3aya/Gudl1GAHJn5/oi+5gi8g/P13jZIgMQ7osX04BTgS8Cy4HfAyc45xYNNtDBUJIWiRP19b5Fu6uie+NGf6FlV9LSfIt2V0X3uHFxPQV8hC9eVM5OUq7ZUfVOVWuRXf5cOa4u5P9XKuQdnEfh0YUUHl1IwREFpBVCkGD/AAAgAElEQVSo0BbpTSSH25sMXACc75yL6lysStIiCaC21s9G2VXRvXlz96/LyIDp0+HAA9svI0fGxUyUQzjcnnJ2EmuqaWL3C7vZ9cQuyp8tp/L1SmgKeUIK5M7M9UX20QUUHFFAxnB1HRHpaCgmiPmKc+6Pg97RIChJiyS46urui+6tW7t+TUmJL7BnzmwrtqdPh6ysoY29F1GYIEY5W2isbKTipQrKny33hfZrlbiG9jVBYHKA/MPyW5ecGTmkpKVEKWKR2BDxwjoWKEmLJLHKSli+HN5+u/3S1UyUqakwdaovsmfPhrlz4eCDIT9/6OMOGurCOhYoZ8eepuomKl6uoPy5YKH9amXrNOstUnJSyD+krdDOPzRfrdqSdFRYi0jycc732Q4ttN96y7dwd7x40swX23Pm+EJ77lyYNWvIxuJWYS2xqLmhmaq3q9j98m4qXq6g4uUKatfVdnpeYHKAvDl5rUvu7FyNpS0JLW4KazPLAF4DvuWce8bMAsBPgEqgCLjEOddFE5SStIj0UXU1vPeeL7KXLYPXXvP3O47JnZbmW7KPPhqOOgqOOAIKCiISUjwX1gPN28rZ8al+e327QrtyaedWbQwC+3YotmflagQSSRjxVFhfAXwfOCWYoO8CHnfO3WNmJwJnOee+0NVrlaRFZMDq6uCdd3yR3bK89177lm0z34p91FEwfz4ce2zYWrTjvLAeUN5Wzk4MzQ3NVL1bReXSSiqXVrJn2R72vLWnU19tDLKnZ7e2aOfOyiV3Zi7pRenRCVxkEOKisDazjwPTgKuAs4HVwDqgyDlXHRwyajewv3NufcfXK0mLSFhVVMBLL8Fzz8Gzz/piO7RVOxDwxfWnPw1f+MKg+mjHa2E9mLytnJ24muuaqVreVmxXLquk6p0qXGPneiNz78zWIjtnZg65s3IJTApgKbE/mo8kr8HOvBhxZpYDnO6c+46ZXRVcPQ8oc85VAzjnGs1sHXA0sD4qgYpI8sjPhxNP9Av4LiRLlviZJB97DJYuhUce8ct3vwvf/jZcdRWkJ0cLnPK2dCclM4W8g/PIOzgPvurXNdU2UfW2L7b3vOlbtaveqaJuUx11m+rY8dCO1ten5qaSc0BO+4L7gFxSc+J3rHpJHjFRWAOXADd2WDcG2NFhXSUwuqsdlJaWMmdO2x8PCxcuZOHCheGMUUSSWXY2HHOMX669Fj780BfYf/oTPPMMXHcdfPAB/PWvfdrdokWLWLSodY6W4ZEKO4IGlbeVs5NLalaqH03kkLZfdlyTo3pNNVVvVbUW23ve2kP91vrWPtyhsiZmkTMjxy/7+9vA1ACpWSq4JfL6mrOj3hXEzD4BZDrnHgw+Xo//SXEOcIZz7mMhz30d+Ktz7qaO+9HPiiISNc8+61u2Wya4mTixXy+Pt64g4cjbytnSnfrSet+i/VaVL7bf3EP1yurO/bYBUiAwJdBaaLcW3FMCpKRrvG2JjFjvCvJdYLa1zY5WBDwA/Cx4P1Qe0M1MECIiURI6nXrH0UUSk/K2RExGSQbDjh/GsOOHta5rbmimZk0NVcurqHq3yt8ur6Lm/RpqVvml7P6y1udbupE9NZucGTlk759Nzv45ZE/PJrCPCm6JrFgorM8EMkMevwx8B1gCfM/MAs65GjNLB8YBz0YhRhGR9pqafH/r226Df//brzvjDNh33+jGNTSUt2VIpaSnkLNfDjn75bRb31TbRPXKaqqWV1H9bnVr4V27rra1+A5laUZgcoDsadlkT89ud5uWFwslkcS7qP8rcs6Vhj42syag1Dm30cweBo4DHgaOB+5zzm2JQpgiIr5f9fPPwxNPwIMPQmkwfWVk+AsYr7kmquENFeVtiRWpWankzcojb1Zeu/WNexqpfq+6tXW7+r1qqldWU7u+luqV/j7/br+vjDEZ5EzP6VR0Z4zKIOTXGZEeRb2w7sXXgZ+Y2UxgFKArW0RkaFRUwJtv+uWNN+DFF/2sjaH22Qe+/GVYuBD22is6ccYe5W2JurTctE4XS4Kftr16dTXVK3xx3Xq7upr6LfXUb6ln15O72r0mtSDVF9nTstsK72nZZE3KUrcS6STmCmvn3ISQ+7uA86MXjYgkvN27YeVKWLHC365cCcuX+xE+OsrJgcMP95PELFgAM2b4iWOSnPK2xIvU7K5buF2To3Z9LVUrqtoX3CuqadzVSOWSSiqXVHbYGQQmBQjsGyB732x/OzWb7H2zyRitVu5kFXOFtYhIWDkHO3b40TrWrfPL2rXw/vu+mN62revXZWT4wnn2bL/Mnetvk2ScapFkYqlGYJ8AgX0C8Km29c45Gkob2hXbVSuqqFlVQ+2GWmrW1FCzpoadj+xst7+UnJS2YnvfbLKntt1PK1Dplcj06YpIfKurg61bYfNmv2zZ4m/Xr28rpvfs6f71WVkwdSpMnw7Tpvll+nS/qIgWSWpmRsaIDDJGZFB4VGG7bU21TX5UktU1VK+upmZV8HZ1DQ1lDex5Yw973uice9JHpLcrtFvuByYFSMlU15J4p8JaRGJTfT189JFvUd6+3d9++KEvnFuK582b2y4g7El+Pkya5MeXnjjR3580yRfP48ZBir7MRKR/UrNSyZ2RS+6M3E7bGnY2tBbZ7YruNTU0fNTA7o92s/v53e1flAJZE7K6LLozx2Rqmvc4ocJaRIZOXZ3vltFSMIcWzR1vd+7sfX/gx5AePRrGjm1bxoyB8ePbCumiIvWFFpEhkz4snYJDCyg4tKDdetfsqNtc11Z0r2orvmvX11K71i881n5/KYEUAlMCbUX31OCoJVOzSctXKRdL9GmIyMA0N8OuXVBW5luNy8raLx3XlZZCZWXv+22RmgojRsCoUTByZNttaPE8dqx/TqqmNBaR2GcpRta4LLLGZfnBKEM01zVT80HnbiXVq6ppKG2g6u0qqt6u6rTPjFEZBKYGi+3gEpgaIGtCFilp+jVuqKmwFkl2zkF1tW8h7mrZsaPtfmiRvHOnL677IzUVhg+HkhI/PF1owTxqVPv7xcXqoiEiSSMls+tJcAAadjVQs8YX2dWrgoX3Kt+1pH5bPfXb6tn9bPuuJZYenAwnWGiHtnKnD9P1I5GiwlokUTjnW4S7K5C7K5Z37vT9mQeisNAXyi3Fclf3Qx8XFKhLhohIP6UXpZN+SHqncblds6N2Y21rod26rPTjclev8KOZdNrf8PQuW7k15fvgqbAWiSVNTX5ikvJy382ivLxt6erxrl1thfKuXdDYOLDjBgIwbFjPS1FR+4K5uFijZoiIRJGlGIEJAQITAgz7xLB22xr3NLZ2JQlt5a5eXU1DWQMNZQ1UvFjRfofBsbmzp2aTvV82OfvnkDPDT4qTmq0ud32hwloknJzzQ7v1pSgOfdyyrqKi92P0JDe39wK5qyUQCM/5i4hITEjLTSPvoDzyDuowGY5z1G2p69TK3XFs7h0P72h7kUHWpCxyZuT4Yju4ZE/L1hCBHaiwFgnlHNTU9K0I7u5xf/sdhzLz3SUKC9uWoqLuH3dsUc7ICN97ISIiCcfMyBqbRdbYLIqOK2q3rammqa0v93vVVC2vourdKj9qyQe11H5Qy44HQgruVAhMDrS2bOcemEvu7FyyJmYl7cyTKqwlMdXXt3WT6M9tefnA+xu3yMnpuRju+Dj0fl6eRrgQEZGoSA2k+uL4wPZjczfXN1O92hfa1e9WU/VuFVXLq6j5oIaaVX4pu7+s9flphWnkzvJFdsuSPS07KUYpUWEtsau52XeNaLnArj8FclXnIYn6LCPDF7t9LYZDHxcUqN+xiIgklJSMlC4nw2mqaaJ6VbBle3kVe970s002fNRA+TPllD9T3raPrBRyDsgh7+A88g/NJ/+wfAJTAgnXsq3CWoZGfb2/wK6ncY5bLsALbT0eaLeK1FRf7LZ0kQi97WpdSyFdVOSnuBYREZEepQZSyZuVR96stn7czjnqP6xnzxt7qHyjsnVq99p1tVS+Vknla5Vs/c1WANKK03yRfWg+BYcVkHdIHml58V2axnf0Eh29TQzSVeE80Ivy8vL6VhR3vM3L07BuIiIiQ8zMyBydSeboTIpPLm5d31DewJ4391D5WiUVL1dQ8XIF9dvq2fnITnY+EpxpNwXy5uRRdGwRhccVUnB4AamB+OoeqcJavOZm31Lc0xTTLbelpQObGKS4uPtxjocP73whXmGhulWIiIgkgPTCdIrmFVE0z18w6ZyjdkNta5Fd8XKFL7xfraTy1Uo23rgRyzQKPl5A0fwihi8YTvZ+2THfdUSFdaJrbvaF8ObNbcuWLb5ADi2WP/qof2MgFxV1Lox7KpoLCjSLnoiIiAC+ZbtlDO6RXxgJ+LG3dz+/m13/3UX5f8vZ8+Yeyp8up/zpctZdvo7A5ADDTx1O8SnFFHy8AEuJvSJbhXU8c84XxRs2tC+cQwvoLVugoaFv+yssbD+ldHfTTY8YoZZkERERCau03DSKP1lM8Sd9F5L6snrKnyln5yM7KXuojJr3a9h00yY23bSJzHGZjPrKKEaeNZLsydlRjryNCutYt2cPrFsHa9f6peP92tre91FcDGPHti1jxrQVyqHFsi7aExERkRiRMTyDEaePYMTpI2hubKbipQrK/l1G6f2l1G2oY8N1G9hw3QYKjipg7EVjGb5geNRbsVVYx4LmZti0CVasaL+sXu27aPRk+HAYPx723rutaO5YRGtWPREREYljKWkpFB5VSOFRhexz0z6UP1fOtru3UfrPUnY/t5vdz+0mMCXAuEvGMfKskVEbM1uF9VDbvRvefBNefx3eeAPefRdWroTq6q6fn5EBEyfCpEltS8vjiRMhP39o4xcRERGJIkux1gshp9wyhW13bWPTzzZRs6aGVeetYtPPNzH55skM+8SwIY8tZgprMzsNuBnIB/4MXOycazSz4cD1QDmQDlzinOtjp+Eoq6+HpUvhhRdg2TJfTL//ftfPHTECpk9vv0yb5lucddGfiMSYhMzZIhJ30vLSGHvhWEZ/YzSl/yhl3ZXrqH63mrdPfJuRXx7J5F9NJr1w6K4Li4nC2szGAacCpwP7A7cDm4CbgHuBS51zS8zsa8C1wGXRirVH9fXw/PPw7LP+9pVXOveBzsiAAw+Egw6C2bPhgAN8ET1s6P+qEhEZiITJ2SKSMFLSUhh55khKPlvC5l9tZv3V69n+p+2UP1vOAQ8d0Gma9kiJicIaGA+c65xrBJaZ2QzgGDN7AZjqnFsSfN5DwCozu845N4g5q8No+3Z49FF4+GF44gl/sWGo/faDI4+EQw7xxfR++/niWkQkfsVvzhaRhJaSmcK4741j+ILhrPjSCiqXVvLGEW+w//37M+z4yDdixkRh7Zx7vsOqLUAecCywLuR5W4IDg88Bnh2yADvaswfuvx/+9Cf473/9sHctZsyA+fPhqKPgiCP8xYUiIgkk7nK2iCSd7KnZzHp+FivPXknpPaUsX7Ccmf+dScFhBRE9bkwU1l2YC/wUWAjs6LCtEhjd8QWlpaXMmTOn9fHChQtZuHBheKNauRJ+/nP485/bLjbMyIDjj4dPfQpOPhnGjQvvMUUkIS1atIhFixa1PIz3v8BjM2eLSFJLzUplv7/ux6rAKrbdvY13Pv0Oc9+eS+bozH7vq685O+YKazObDJQ6594yMwfUdHhKJtDpQpiSkhKWLl0amaDeeQeuuAIeeqht3RFHwJe/DGec4WchFBHph9BC0szKohzOgMVkzhYRCbIUY9/f7Uvdljp2Ld7FyrNWcuATB/Z7avS+5uyYGm7CzFKBrwKXBldtAYpCthuQA2wdkoB27IDzz4dZs3xRnZUFCxf6Maaff97fV1EtIkkq5nK2iEgXUtJSmPbHaaQPT2fXk7so+1fk2jJiqrAGvgvc5JyrCz5eDEwN2T4BqAci38zx+OO+v/Qdd/jh7i64wE8d/tvf+mHwREQkdnK2iEgPMkdlMuHaCQCsvXwtrtn1/IIBipnC2swuB94AAmY2yczOxY+D+mHwinOAE4FbnHP1EQvEObjhBvjkJ2HbNj+ix/Ll8Ktf+bGmRUQkdnK2iEgf7XXeXmSOy6RmVQ27/rsrIseIiT7WZnYl8KMOq1c65+40s88DV5nZeiAXuCJigTgHF14It94KZnD99fD970NqasQOKSISb2ImZ4uI9ENKegp7nb8X63+wng/v+JBh88M//F5MFNbOueuA67rZtgH43yEJ5NprfVGdmQl//SucdtqQHFZEJJ7ETM4WEemnkV8ayfofrGfn4ztpbmgmJT28nTdipitI1D3yCFxzje9Pfc89KqpFREREEkxgQoDs/bJpqmhi94u7w75/FdbgR/845xx//8Yb4ZRTohuPiIiIiERE0XF+8KKKlyvCvm8V1gDXXQelpXD00XDxxdGORkREREQiJG9OHgCVyyrDvm8V1h9+CLfd5i9W/OUvfVcQEREREUlIuQflArDnzT1h37eqyN/9Dhoa4NRTYebMaEcjIiIiIhEUmBwAoG5DHa4pvONZJ3dh7Rzceae//41vRDcWEREREYm41KxUMvbKwDU66jbX9f6Cfkjuwvqdd/xsiiNHwrHHRjsaERERERkCWROzAKhZVxPW/SZ3Yf3II/725JPVt1pEREQkSWSOyQSgflt4J4ZN7mryuef87fz50Y1DRERERIZM+oh0ABo+agjrfpO7sE5P97MsHnFEtCMRERERkSGSOSaTjL0ywMK735iY0jxqHnwQ6ushIyPakYiIiIjIEBl/2XjGXzY+7PtN7hZrUFEtIiIiImGhwlpEREREJAxUWIuIiIiIhIEKaxERERGRMEjqwnrRokXRDmFIJdv5QvKdc7KdLyTnOSczfd4Do/dt4PTeDVwyvncqrJNIsp0vJN85J9v5QnKeczLT5z0wet8GTu/dwCXje5fUhbWIiIiISLiYcy7aMYSFmZUCG/r5suFAWQTCiVXJdr6QfOecbOcLiXHO451zJdEOYigNMGdDYnze0aD3beD03g1cor533ebshCmsRURERESiSV1BRERERETCQIW1iIiIiEgYqLAWEREREQmDpCyszSxgZreY2Q1mdruZ5Uc7pkgws/lmtsTMJoSsS8hzN7PTzGydme0ws1+aWVpw/XAz+62Z/T8z+5mZpUc71nAxs6PNbJWZ7TKzX4SsT8jPuIWZZZjZW2Y2L/g4oc9XPH3Og9PV94H0rLvvFeldd99PySApC2vg18ALzrnLgQeA30Y5nrAzs5FALnBIh00Jd+5mNg44FTgduAg4D/h2cPO9wJ3OuUuA1cC1UQkyzMysADgcOBT4EvB/ZnZ8cHPCfcYdfA+YEPI40c9XPH3OA9TD94F0o5fvFelBL99PCS/pRgUxs9HAOqDIOVcd/At0N7C/c259VIMLMzNLAZqAic659Yl67mZ2JPCyc64x+PgnwP7Aj4D7nXOjg+vHAKuAkc65qmjFGw5mFgBqXfA/sJktA76LP7+E+4xbmNnHgWnAVcDZ+D+WEvZ8xUvU3DWUOn4fRDmcmNfd94pz7uToRhb7uvt+cs49Hd3IhkYytljPA8qcc9UAwf8064CjoxlUJDjnmjusmkcCnrtz7vmW5Be0BdgIHIs/v5bnbQnenTOE4UWEc64mJGnlACuAZ0jQzxhaz/N059ydIavnkaDnK+3MQ5/zoHTxfSA96OF7RXrRw/dTUkjGwnoMsKPDukpgdBRiGWrJcu5zgd+QBOcbbFV5DEgHsknsc74EuLHDukQ+X2mjz1mireV7Rfqoi++npJCMhbUDajqsywQaohDLUEv4czezyUCpc+4tkuB8gfeBu4HjgJ+SoOdsZp8AljrnPuqwKSHPVzrR5yxR0+F7Rfqu4/dTUkjGK1y3AEUd1uUBW6MQy1BL6HM3s1Tgq8ClwVVbgH1DthuQQ4KcL4Bz7kPgTjNrxl/U9zyJ+Rl/F5jtP0LAn+MDwM9IzPOV9hI6d0ns6uJ7Rfqoi++npJCMLdZPA6ODnesJDr82Dng2qlENjUQ/9+8CNznn6oKPFwNTQ7ZPAOqBpUMc11BYii8+EvUzPhOYFbJsxV+lfxeJeb7SXqL+u5bY1/F7Rfqv5fspKSRdYe2c2wY8jP9pAuB44L6QC9sShrU17xkk9rmb2eXAG0DAzCaZ2blAOfChmc0IPu1E4BbnXH204gwXM8s2s9khq04Cfpmon7FzrtQ5t7llwY9uUOqc20gCnq+0l6j/rodSx+8D6V1X3yvBbiHSg+6+n6IVz1BLuuH2AMysCPgJsB4YBVzmnNsT1aDCzMzy8ONH/ho/dvOtzrnSRDx3M7sSP7ReqJXOuelmNh4/NNt6/DiuV3S40jsumdkhwCPASuAV4FXn3L3BbQn3GXdkZuuBs51zzyTD+Upy/LuOlO6+D6IbVWzr6XslGvHEk56+n5JBUhbWIiIiIiLhlnRdQUREREREIkGFtYiIiIhIGKiwFhEREREJAxXWIiIiIiJhoMJaRERERCQMVFiLiIiIiISBCmsRERERkTBQYS0iMcXMSszsQTPbYGaXRjseERHpmfJ2G00QIyIxxcw+ATwN5ONn7trXObczulGJiEh3lLfbqLAWkZhlZs8BxybCNPQiIskg2fO2uoKISEwys72BJ5I1OYuIxBvlbRXWEgFm9n0zuy3acbQIRzxmdryZvWxm88IUVkIxs5vCvD8DPgfcGM79ikh7yteDZ2YTzOx3ZvbbMO0v1czON7MNfXjuyWa2ZKDnGs7crbztqbCOM2Y23czuN7M7zOwtM3Nm9qVox9XBo8Dfoh1EiHDEswmYGYZYYpKZnWRmZWb2oZkdZmZ7B/+dOTO71cwKg887xMyWm9nNZhYIrhsJbAtzSGcCv3fONQb3LxJ3lK8HJB7z9UfB28ww7S8FKAfG9eG5S4Hp3W00s/3N7J7gv70rzWxUyLZw527lbVRYx6N/Az90zp0HzAb+2LLBzA40s6OjEVTosZ1zy51zL0Qjjq70J57u3kPn3CqgLJxxhfPzGuy+nHOPAr8GGp1zLzvnNgHnA/XAEudcefB5r+IT+fecczXBl38aeHhQJxDCzH4AXAM8aWYrganh2rfIEFO+7qd4zNfOuWpga7iO5ZxrAJb18bnb8UV4d9vfBRYFH97knAstpMOWu5W326iwjiNmNgzYF6gCcM41A1cBzswK8EnbohBX1I4dTn04j+YhPFY09nUPMNbMDgZwzu0A/gt8KuRYhcDW4L+9FtOccysHeexWzrkfOeemOOfmOOemOeeeC9e+RYaK8nVkxWC+DvdIEP2Jv7fnNgVvO/Z7DlvuVt5uo8I6jgSHrlkFPG5mhwbXbQBexfdrmggsNLNzAMzsWjP7QfBnoNZ+VGY2x8x+ZmbfNrNqM1tlZt8xs1PN7DUz+4yZ/cvMKszsqyGv63J/HY59vpndHdrXLPhz6C1m9mMze9rMDgmu7/F4Ia83M/uqmb1rZvPN7G0z225mZ/fhGFNC4+nlmJ3ewy4UBX/arTSz8zu8p9cE35v7g4m4r+9Zp2OZ2elm9qqZnW1ma8zsxh7219Vn32U8PQm2bLwLnBGyugQ42YLdPoAFwIMhcQaA6pDH08zsRjO7ycyeNd+l5PNmttTMTjGzv5lZuZl9y8wODp7jDjM7qbf4ROKJ8nXy5OuudHeOPR3HzCaa2Z3mW39/2pfjdHHcR83sL2Z2eC/PU+6OFOecljhagP2BDfi/UO8ERoRsWw/MC96fAawN3i/B/zVdjP+LeyMwJ7jtNuCl4P0MYCdwK34syi8AW3raX8dj4/uY/Ra4O7g+C1+sFQYfnw6UAnk9Ha+L8x4VPOblwedeB9QC43s5Rsd4ejxm6HvYRQzrgZuB7ODrNgbX5wN/D3nei/ik2Kf3rJtjDcP/vPcL4BDg+J721+Gz7zKePv77uhL4IHh/H+DHQB1wWsi/Fwt5/gLgyOD9bOA1IDv4+Lng+5AP7AauBQLA54OPP4f/4/5igv8GtWhJpAXl66TI18Ht1/Txfezps34NmB7yGtfHf2ctn+dhwNVdbJ8XPE5ayDrl7ggtarGOM863Kh6ATzRfAd41s64u0lgDfNnM0oEjg+ty8f+R9wZa+seuwP+HwTlXD1QC/3TOVeD/o43qZX8d46uj/cUQJwM1LthHF7gf/x/ytF6O13G/Lfv8a/C51+H7/87v5Rjt4unPMbvxkPP96V4D9go5x+Jgi9K38e8p9PE96+Z8d+IT2APOuVedc0/2Y3/dxdMX9wCTzHcH+SzwO+A/wOfMLBeodMHMG3Q48FLw/gL8l1dLK8hngB8F3+dy4Cnn+2UvBfKdc/9w/ufxN+jfZyASF5SvkyNfd6Hbc+zhOMcDBc65lniW9vOYc4EvOud+2MfnK3dHSFq0A5C+M7NUYG/n3HrgQjP7Pf4K6kXAx0Kf65yrM7MxwCXAH1p24Zz7yMyWAsfh/6KeiC+mWl8acr+RYHeh7vbXTaih+5gCpIfE1Wxm64GxPR2vN865WjNbCxQCw/txjAEfs4NG2v7/jAPWOed+0fFJfXnPzOzL+AK2xVTnfzJ2obH24zPoNp7eOOfWmNkb+O4gBc65tWZ2D74V6QzgkZC4W/5ttPTfG49vYWrZ147QXYfcb6K9ZtQtTRKM8nW780uGfB2q2/exh8/mQPwfEQN1LDDFzHKdc3u6eU5z8ByUuyNIb0h8SSek/6tz7i38T23TOj7RzI4FvuGcu875ER5C/R8w1/ywT28D/6+3A/eyv55sxLeAhg5DZPi+h4OVAayO8DH64kPa90PGzD7Wj/fsQWBWyNLl1eX92F+X8fT9dPg7cA6wNiS+VOBb+J9NW3wM31+0xRbgCDPLCTluj/38RBKY8nV7yZSvuz3HHo5TEXxNBgNzE761/5fdbN/t2i46V+6OIBXW8ecCMxsd8ngM8Hjwfj3+Yo1p+KGdCsws08yOD24fYWbD8X39fo//Geclgj8tBuYg5REAAALsSURBVKXQ4a90M7Ne9tfx2Bayj38De4D/Ce6rMHi8lpbP7o7Xnb2CzxkTfO3jfThGaDy9HTP0PDrqqvXC8K1QOcDDZvYJ8xfJ9PYZtB4L2Ms5tzJkaQg5RGrI/Z72Fxp3l/GYWYaZLQp28+jJP/A/Qf8TwDlXCTwGvOjajwZyAr6bSIuW9/tv5sfC/ibQMpZpx8+g4+cc1yMUiHRD+Zqkydd9fR+7PA4+l2bg//gC/+sEZja6D7nb8K3HXwbOMLPPBV97sZldiW+t/2/I85W7IykSHbe1RGbBXxDhgF3An4G78T8lFQW3Xw1sBk7B/7SzBt9CcDK+peNufJJ6EH8hRX1wf9X4/2ifBBrwF8jsDVwW3H5uT/vrcOxvAc/j+63tH9w2B9/SeTO+r+GM4Ppuj9fN+bvg67+F75pwQMi27o4xITSe3o4Z+h52OPYp+Iv4bg2+Fy2vOye4/WhgefCzuTn4Pvf1PTuli3M9PXi8vwFjgut6+kzb7aubeMbiJzL4ax/+rf2qw+PPA8d0WHdTF687CngP2AFcGVx3Ev7f2i3AaIJDjgFfw3/x/i54rp+K9v8xLVrCtaB8nUz5ehzwAr7VfVYv59jTZ30y/pfCV/AXkq8GvhE89y5zN/5XkQZ8S/Uw4Cl8UX8e8FX8eN4347uhtLxGuTuCiwXfPEkSZjYC+Jpz7trg4xR8K8pXnXNXRjW4XpiZAyY632dRBsjMvuecG9BQTiIydJSvJZRyd3xQV5Dk8018t4BhIeuOJeSitFgU8tOTfnYaIDNLNbMvAvdFOxYR6RPla1HujjMqrJPPb/B9sVaY2SZgMbDZOfdydMPqXvBijm8FH55tftg36b9c4N/OubW9PlNEYoHytYByd1xRVxARERERkTBQi7WIiIiISBiosBYRERERCQMV1iIiIiIiYaDCWkREREQkDFRYi4iIiIiEgQprEREREZEwUGEtIiIiIhIGKqxFRERERMJAhbWIiIiISBj8f9YVHCBCenAsAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(12,8))\n", "plt.rc('font',family='Times New Roman')\n", "params = {'mathtext.default': 'regular' } \n", "plt.rcParams.update(params)\n", "\n", "plt.subplot(2, 2, 1)\n", "plt.plot(vehicle.v_kmsc, vehicle.h_kmc, 'k-', linewidth=2.0)\n", "plt.xlabel('Speed, km/s',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\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", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 2)\n", "plt.plot(vehicle.acc_net_g, vehicle.h_kmc, 'b-', linewidth=2.0)\n", "plt.xlabel('Deceleration, Earth g',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\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", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 3)\n", "plt.plot(vehicle.q_stag_total, vehicle.h_kmc,'r-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-rate, '+r'$W/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\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", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "\n", "plt.subplot(2, 2, 4)\n", "plt.plot(vehicle.heatload/1.0E3, vehicle.h_kmc, 'm-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-load, '+r'$kJ/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\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", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.savefig('../plots/opportunity-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/opportunity-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/opportunity-mars.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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": 2 }