{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 26 - Galileo - Jupiter" ] }, { "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 Galileo entry probe.\n", "https://en.wikipedia.org/wiki/Galileo_Probe" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"JUPITER\") \n", "planet.h_trap = -130.0E3\n", "planet.loadAtmosphereModel('../atmdata/Jupiter/jupiter-galileo-asi.dat', 0 , 1 , 2, 3, heightInKmFlag=True)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Galileo', 335.0 , 256.0, 0.00, 1.26, 0.0, 0.22, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(450.0,0.0,0.0,47.37,0.0,-8.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 (9400.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": "iVBORw0KGgoAAAANSUhEUgAAAt4AAAHtCAYAAAAuknnBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZwcZbX/8c/JZDLZk0kyZIMwrAHZMSxXLibIqiBXEFARZVFGkYteFW+8oFxcWVx+KF7UAVFBcUFQDOASlEVUlkEkLBIIBBKyTvYVkpk5vz+e6kmnM93TPdPdVd39fb9e9arqerqrT1U6p888/VSVuTsiIiIiIlJaA+IOQERERESkFqjwFhEREREpAxXeIiIiIiJloMJbRERERKQMVHiLiIiIiJSBCm8RERERkTIYGHcA5TJu3Dhvbm7O+/nt7e00NTWVLqCEqJX9BO1rNaqV/XziiSdWuHv172iaQnN2Sq18JnLRMdAxSNFxiOcY5MrZNVN4Nzc309bWlvfzp02bVtDzK1Wt7CdoX6tRreynmb0adwzlVmjOTqmVz0QuOgY6Bik6DvEcg1w5W0NNRERERETKQIW3iIiIiEgZqPDOoqWlJe4QyqJW9hO0r9WoVvZT8qfPhI4B6Bik6Dgk7xiYu8cdQ1lMmzbNa32ck4hUJjN7wt2nxR1HOSlni0ilypWz1eMtIiIiIlIGKrxFRErM3ens7KSjoyPuUKqWO3R2QldX3JGIiGSnwltEpMTuvPNOBg4cyHve8564Q6lat90GAwfCBz4QdyQiItmp8BYRKbEBA0KqrZVzauKkQywiSabCW0SkxFKFd5fGQZSMWZir8BaRJFPhLSJSYiq8Sy9VeIuIJJkKbxGRElPhXT7q8RaRJFPhLSJSYiq8S09DTUSkEqjwFhEpMRXepaehJiJSCVR4i4iUmArv8lGPt4gkmQpvEZESU+FdehpqIiKVQIW3iEiJqfAuPQ01EZFKoMJbRKTEVHiXj3q8RSTJVHiLiJSYRd2xKrxLR0NNRKQSDIw7gExmNgh4HPiEuz9gZkOAa4H1QCMw093XRc/N2iYikhTV3uOdhLytoSYiUgmS2OP9GaA57fENwMPufhlwF/D9PNtERBKh2gtvEpS31eMtIkmWqMLbzN4CLAFWR48nAWcDs6Kn3AecambNudrKGbOISG+qufBOSt7WUBMRqQSJKbzNbBhwhrvfnLZ6BrDC3TcBuHsHMB+Y3kubiEhiVGvhnaS8rcJbRCpBksZ4zwSuzlg3GViZsW49MAnoyNG2g/b2dqZNm9b9uKWlhZaWlv7EKyKSl1Th7QVUha2trbS2tqYejit+VEVRsrxdaM7WGG8RiVO+OTsRhbeZnQi0ufvyjCYHNmesawC2Rsu52rbT1NREW1tbf0MVESlYX3q80wtNM1tRksD6odR5u685Wz3eIhKHfHN2Igpv4FLgENvWZdFIOOnmm9FyuhHAYkJyz9YmIpIYVTrUJFF5W0NNRKQSJKXwPpvQ65Hyd+BTwKPAZ8xsiLtvNrN6YArwINAJ3JilTUQkMaq08E5U3tZQExGpBIk4udLd2939tdRESM7t7r4AuBs4NnrqccAd7r7I3Zdmayt3/CIiuVRj4Z3UvK0ebxFJsqT0eOdyEXCtmR0ETABa8mwTEUmEaiy8e1H2vK2hJiJSCRJZeLt7c9ryauDCLM/L2iYikhS1UHjHnbc11EREKkEihpqIiFSz1AmI1Vx4J4V6vEUkyVR4i4iUWC30eMdNQ01EpBKo8BYRKTEV3qWnoSYiUglUeIuIlJgK7/JRj7eIJJkKbxGRElPhXXoaaiIilUCFt4hIianwLj0V3iJSCVR4i4iUmApvEREBFd4iIiWXKrxd3bElox5vEakEKrxFREpMPd6lp8JbRCqBCm8RkRJT4V16upygiFQCFd4iIiWmwrt81OMtIkmmwltEpMRUeJeehpqISCVQ4S0iUmIWVYUqvEtHQ01EpBKo8BYRKTH1eJePerxFJMlUeIuIlJgK79LTUBMRqQQqvEVESkyFd+lpqImIVAIV3iIiJabCu3zU4y0iSabCW0SkxFR4l56GmohIJVDhLSJSYiq8S09DTUSkEqjwFhEpMRXe5aMebxFJskQV3mY23czmmtlqM7subf0QM7vezL5qZt81s5H5tImIJEHqOt7ujldRZZiknB39bYP+thGRJEtM4W1mo4CjgCOBc4CLzey4qPkG4GF3vwy4C/h+2ktztYmIxM7Mtiu+q0HScvagQWG+ZUsxtiYiUhqJKbyBLcBV7r7a3e8B5gCdZjYJOBuYFT3vPuBUM2vO1VbWyEVEepEablIthTcJy9kNDWH+xhv93ZKISOkkpvB2980efSOZ2TDgX8ADwAxghbtvip7XAcwHpvfSJiKSGNV22/ik5exU4f366/3dkohI6SSm8E4xs6OB3wH1wFBgMrAy42nrgUm9tG2nvb2dadOmdU+tra1Fj11EJJtCT7BsbW3tzlfAuNJF1j9JydmDB4e5erxFJA755uyB5Qspb/OAHwHXAl8DXgY2ZzynAdgaLedq69bU1ERbW1tRAxURyVehhXdLSwstLS0AmNmKkgXWf4nI2RpqIiJxyjdnJ67wdvclwM1m1gV8BvgL0JjxtBHAYsBztImIJEa1XlIwKTlbhbeIVILEDTVJ0wYsAu4HJpnZEAAzqwemAA/20iYikhjVWniniTVnpwrvzZn96SIiCZKYwtvMhprZIWmr3gF8y92XAncDx0brjwPucPdFudrKFbeISD6qrfBOWs4eNgzq6sLJlbqkoIgkVZKGmuwP3GNmzwOPAI9Fl6gCuAi41swOAiYALWmvy9UmIpII1VZ4k7CcbQaNjbBiBaxeDePHF2OrIiLFlZjC290fA5qytK0GLiy0TUQkKaqt8E5izh4zJhTeq1ap8BaRZErMUBMRkWpWbYV3Eo0ZE+arVsUbh4hINiq8RUTKIFV4d3Z2xhxJ9VLhLSJJV5TC28zeUYztiIhUq1TBXVdXF3Mk1ZuzJ0wI80U6vV5EEqqgMd5mdhLwOcIJMamifQAwHhhS3NBERKrH1q3hHjH19fVle89ay9m77hrmCxbEG4eISDaFnlx5K3AV8CSQGqhowNnFDEpEpNp0dHQAMHBgWc9pr6mcnSq8X3013jhERLIp9BvgZXf/ZubK6HJSIiKSRRw93tRYzlbhLSJJV+gY74vM7Iwe1p9ajGBERKqRu8dVeNdUzt599zB/4YV44xARyabQwvso4CYz60ybuoDvliA2EZGqkDqxcsCAAd1XNymTmsrZu+wCo0ZBezssXRp3NCIiOyp0qMkXCXcZawNS18Qy4APFDEpEpJqkervLPL4baixnm8GBB8Jf/gJz5my7yomISFIU2vXyO+Aud3/Z3V+NpleAG4sfmohIdUidWFnmYSZQgzn7kEPC/JFH4o1DRKQnhRbePwf+x8zemjZNB/6nBLGJiFSFLVu2ADBo0KByv3XN5exjjgnz2bPjjUNEpCeF/u75eWAv4DzAo3VGuEbsJ4oXlohI9Vi3bh0AI0aMKPdb11zOPuYYqKuDv/8dVq+Gxsa4IxIR2abQHu+vAbu4e7O77xZNzcDpxQ9NRKQ6rF27FoBRo0aV+61rLmePGgUzZkBnJ/zyl3FHIyKyvUIL74Xuvr6H9TqFRUQkixgL75rM2R/8YJjffDO4536uiEg5FVp4fyxzhZmdCVxXnHBERKpPaqhJDIV3Tebsd787DDF57DF48MG4oxER2abQwnuMmX0ewMyGmtlNwE3A40WPTESkSqxZswaIpfCuyZw9bBj813+F5S98Qb3eIpIcBRXe7n4yMNvMrgf+CewNHAT0dGc0EREBFi9eDMDEiRPL+r61nLM//vHQ6/3AA3D77XFHIyIS5LyqiZm9NUvTq8CLwNeBXYF3AZ8sbmgiItVh4cKFAOy8884lfR/l7G1Gj4arroKPfjT0fh97LIwdG3dUIlLreruc4LXAvsAaoKuH9psJl6aaSJUncRGRvkoV3rvsskup30o5O82FF8Itt8Df/gbnngu//S0MKHSApYhIEfVWeH8RaHP35bmeZGZV/7OliEhfLViwAChL4a2cnWbAALjttnA3y3vugcsug6uvjjsqEallOf/2d/d7e0vg0fN+1d9AzOx0M5tvZivN7FtmNjBaP87Mvm9m15jZN82sPu01WdtERJKgo6OD559/HoCpU6eW9L3KmbOhMvL2rruG4nvgQLjmGvjyl3WypYjEJxE/upnZFMKYwzMIP39+GIjOSed24GZ3nwm8QOjRIY82EZHYvfjii7z++us0NzfHcVWTkqmkvH3SSfDjH4MZfP7z8LGPwZYtpX5XEZEdJaLwJpzsc4G7P+HutwD/BxxjZkcCU9390eh5s4BLzGxYrrayRy8iksWTTz4JwAEHHBBzJEVXUXn77LPDnSwbGuB734OjjoIXXyz1u4qIbC8Rhbe7/8XdO9JWLQIWAG8D5qc9b1G0OK2XNhGRRLj//vsBOOqoo2KOpLgqMW+fcUa4vOCuu0JbGxx0EHzpS/D66+V4dxGRIhTeZnaImf3EzD5tZs39DwmAw4DvAZOBlRlt64FJvbTtoL29nWnTpnVPra2tRQpVRKRn7s6f/vQnAI477riCXtva2tqdr4BxxYqpRDkbipy3S5WzjzwS/vlPOOcc2LwZrrgC9t03DEXp6Oj99SIiPck3Z/d2VZN8fBQ4zt3PMbOzzWyZu2/u68bMbE+g3d2fMjMHMrfVAGwFcrXtoKmpiba2tr6GJSJSsH/84x/Mnz+fpqYmDj744IJe29LSQktLCwBmtqKIYRU1Z0Np8nYpc/bo0XDrrfChD8Ell8Azz8B554Xe789+Ft7/fhgypCRvLSJVKt+cXYyhJv8JvDla/jnQ5zPUzawO+Ajw2WjVIqAxrd2AYcDiXtpERGJ3yy23APC+972Purq6mKPpVrScDZWdt2fMgCefDL3de+0FL70Urv29887w3/8Nr7wSR1QiUs0KLrzN7O1m9v5o+XBgr9Q4PXfvcvd1/YjnUuDr7v5G9Hg2kH79rWZgC9DWS5uISKzWrFnDj3/8YwDOPffc2OIocc6GCs/bAwfCBz8Izz0HP/0pHHYYrFoFX/sa7L47nHBC6B3fsCGuCEWkmhRUeJvZVwiXgjoNwN0fAz5oZmf3NxAzuwx4EhhiZrub2QWEu68tMbP9o6edBFzv7lvcvS1bW39jERHpr+uvv561a9dyzDHHcOihh8YSQylzdrT9qsnbAweGK5889hg8+mgYA15fD7Nnh8J8woRw98s//AG29jigUUSkd+YF3EnAzO4nJPBz3f1b0brpwK3uPqXPQZh9DvhSxurn3X1fM9sVuAJ4BRgOXJ46kz5XW6Zp06a5xniLSDksXLiQN73pTWzYsIE///nPHHPMMf3anpk94e4FX/mjVDk72k5J83YScvbq1eEShLfeCn/967b1o0fDO98Jp58eesSHDo0vRhFJnlw5u9CTK2e7+5ro5JmUft962N2/DHw5S9urwIcKbRMRiYO7c/HFF7NhwwZOO+20fhfd/VSSnA21kbcbG+EjHwnTvHnwk5/AHXeEkzFvvTVMQ4fCiSeGm/SceGK4VKGISDaFjvF+OfppcaqZvdfM7gAuBq4pfmgiIpXnuuuuY9asWYwYMYLrr78+7nCUs4tkzz3hyivh6adh7ly4+mo4/HDYtAl+/etQnDc3h0sTfvKT8PvfhzYRkXQFDTUBiO48dh4wBVgC/Mzd7yt+aMWVhJ8tRaS6zZ49m7e//e10dnbyq1/9ine/+91F2W5fh5pEr1XOLqGFC+F3vwuF9n33wfr129rq60Nx/ta3huktb4GRI+OLVUTKI1fOLrjw7mHjdcAMd/9TvzZUYpWSxEWkMj388MOceOKJbNq0iZkzZ3L11VcXbdv9Kbx72JZydols3QqPPBKK8D/8Af7xD0j/ih0wAA45JNyu/rDDQlG+555hvYhUjz6P8Taz3xJOfsllHLACSHQSFxEplXvuuYezzjqLTZs2cf755/PVr341ljiUs+NVXw9HHx2mr3wF1qyBv/0NHnooTI8/Dk88EaaU0aNDEZ4qxA89NFxH3Cy+/RCR0unt5Mo50XPmAl3AhYRLR6V3Q+wGjClJdCIiCebuXH/99Xzyk5+kq6uLc889lxtvvJEB8XVhKmcnyOjR8I53hAnCmO9HHgmXK0xdtnDJknDJwtmzt72usREOOmj76U1vgsGD49kPESme3grvrxGGo6wBMLM93P2KzCeZ2V2lCE5EJKlWrlzJBRdcwG9/+1sArrjiCq688kos3q5K5ewEGzoU3va2MKUsWhSK8NT0z3+GG/g88ECYUurqYOrUUIDvs084iXOffcK6YcPKvSci0lc5C293X5uxaryZ1bl7Z2qFmc0ADi9BbCIiiePu/PznP+dTn/oUS5cuZdSoUdx4442ceeaZcYemnF2BJk+G004LE4Qx4YsWwVNPbZvmzIEXXgh313zuuR23MWXKtiJ8991ht922zYf3NvBIRMqq0Ot4/wp4ysweBtYD+wInEm6GICJS1Z544glmzpzJn/4Uhkf/+7//O7feeivNzc3xBpadcnaFMQtjvHfeGU4+edv6TZtC0f3882H617/C/MUXYcGCMP3xjztur6lp+0J8t91gl122vceoURpPLlJOBRXe7j7bzE4APgYcSLg18Bnurp8tRaRqzZkzhy984QvceeedAIwZM4Zrr72W888/P87x3L1Szq4eQ4fCtGlhStfRAS+/HIrwF16A+fPD4/nzw9TeHqbHHut5u8OGhV73VCGeWp40CcaP3zZpOItIcRTa4427LwY+l3psZgPM7M3u/kSOl4mIVJTOzk5mzZrFt7/9be6//34AhgwZwsUXX8zMmTMZN25czBHmRzm7ug0cCHvvHaZMXV3h5M30Qnz+fHjttW3Txo2hYH/hhdzvM2zY9oV4tmns2HBSaYL/HhWJVUGFt5nd3MPqJqATeFdRIhIRidEzzzzDT37yE2677TYWLlwIwPDhw7ngggv47Gc/y8SJE2OOMH/K2bVtwIDQgz15crjEYSZ3WLcuFOCLFm2bL1wYCvZly7ZNGzeGAv7ll/N738ZGGDcuFOKpKfNx5rr6+uIfA5GkKbTHey/CtV/T77qzM/CXokUkIlJG7s5TTz3F3Xffze23386cOXO623bffXcuueQSzj//fEaNGhVjlH2mnC1ZmYUx3qNGwX77ZX9eqkBPL8R7mpYvh5UrYe3aMF+5srB4hg8PBfvo0WFKLeezbvhwjVWXylBo4f1hd5+bvsLM9iHcjlhEpCKsXr2ahx56iHvvvZd77rmHRYsWdbc1NjZy1lln8f73v5+jjjoq0WO486CcLf2WXqD3NKQl09at4ZKIqeJ75UpYsSL345UrYcOGMEU/NBWkrm5bIZ5emI8aBSNHbj/vad3IkWEcvYp3KbVCT66c28PqN4APA58tSkQiIkW2YsUKHnroIR588EEefPBB5syZg6fdy3vSpEmcfPLJvPOd7+SEE06goaEhxmiLRzlb4lBfv23Md766umD9+nC3z9Wrwzx9ubd1Gzf2rZc9XV1dzwV5rmK9p7mGzEguhY7xns/2P1nWAeOBW4sZlIhIX23cuJEnn3ySxx9/nMcee4zHH3+cl156abvnDBo0iCOOOILjjz+eU045hYMPPjjuG9+UhHK2VIoBA7YVuLvuWvjrt2wJQ1wyC/N168L6zHlP615/Pbxm9er+7cuQIb0X572tGzZMJ6hWq0KHmswGbmNbIu8Clrl7L+dDi4gUl7uzePFinn76aZ5++mnmzJnDU089xbPPPktXV9d2zx08eDBHHnkk06dPZ8aMGRxxxBEMGTIkpsjLSjlbasKgQeGa5U1Nfd/Gli3ZC/VcBXvmfPPmMC1b1vdYzEIRnm+hnm1dlfx4V1UKLbx/4u4Ppa8ws6lmtqe7zytiXCIiQCiwly5dyty5c5k7dy7PPPNMd6G9uoeuqbq6Og466CAOO+wwDjvsMA4//HD2228/6mvz91/lbJE8DRoUrrLSnyuFuodhL30p3NOXN23atq4vY95TGhr6XrSn1g0frt73Yiq08D4JeChj3TLgfuCQokQkIjVp48aNvPjii90Fdmp64YUXWL9+fY+vaWxs5MADD+SAAw7ong466CCGDh1a5ugTSzlbpIzMQqE6fHi4CVFfdXT0r/c9Nb3xRrjazPLl/dunESNCMZ46cbWxEcaM6Xk58/HAgu8YU93yOhxm9p/Al4ERZjazh6fMLmpUIlJ1urq6WLJkCS+//DIvv/wy8+fP32558eLFWV/b2NjI1KlTmTp1Km9605u6i+1JkyZV5djs/lLOFqlsAweG4nXMmL5vwz0MecmnYM9VxG/YEJZT13wv1IgR2Yvy1D6OGxeGCaXmY8aEk12rUV6Ft7t/x8weAs4Fvp3RvMXdlxQ9MhGpKJ2dnSxZsoSFCxeycOFCFixY0F1cz58/n1deeYU33ngj6+vr6+vZc8892XvvvbuL7NRUKXeJTArlbBExC5dIHDoUJkzo+3Y6O8MVZ9JPXk1Nq1Zlf7xqVXj++vVhWrCgsNjHjNm+GE/Ne1q3004weHDf97Gc8v4BwN3nmNln3X1rZpuZTXH3Ag5pdmZ2PKGn5j3u/kq0bghwLbAeaARmuvu63tpEpDjcnZUrV25XVGcuL168mI6Ojpzb2Wmnndhtt93Yfffd2X333bdbnjx5MgP1m2TRKGeLSDGkXyO90CvOpC4TmatIT13XfcUKaG8PU/p14PPV2AgTJ4YhPqkp/fHEiWGKu0DP+S1nZpcC97r7c9Gq/zazzM5/A44BZvQ3GDMbDwwHDs9ougH4vbv/wsxOAr4PvC+PNhHpxYYNG1i8eDFLlixh8eLF202pda+99hqbN2/udVvjx49nl1126Z7SC+vm5maGDx9ehj2qXcrZIpIk6ZeJbG7O/3UdHaH4bm/fVpCnF+aZy8uXbyvon3su97YnToQ99oDdd9827bEHTJ0KY8f2a3fz0lv30m7A6LTHhwLDgKVsuzyVAf34EWMbd19mZnelrzOzScDZwMXRqvuAO8ysGdiSrS3V8yJSi9yddevWsWzZMpYtW5a1oF68eHHWExczjRo1qrugnjJlyg7LO++8c9XceKaCKWeLSMUbODAMH9lpp/ye39UVeseXLIHFi7dNmY+XLg3rliyBhx/ecTu77AKHHAJHHAGnnAIHHFD8u5nmLLzd/eKMVV8BnnH3LekrzexNxQrI3bsyTpaaAaxw901Re0d0U4jpwNYcba8UKyaRJOjq6mL16tXdxXTmtHz58u0e5xpPnW7w4MFMmjSJSZMmMXHixO7l9HWTJ09m5MiRJd5D6S/lbBGpRQMGbBv/feCB2Z/X0QGLFsHLL8NLL4V5avlf/wqXbly4EH77W7j8cjjuOPjNb8INjYql0AGVB7j7P9JXmNlUYA+gl879PpsMZI7yWQ9MAjpytG2nvb2dadOmdT9uaWmhpaWluJGKFGjTpk2sWLGC9vb27nm2gnr58uW9jqFON2zYMMaPH8/48eN3KKbTi+zRo0fryiAJ1NraSmtra+phX88uVc4WEYkMHBjGqe+6KxxzzPZtnZ0wezZ85ztwzz1h3X33QVsbTJ/e+7bzzdn5Xk7waMLPk8dHvRPpxhDOmp+Vz7b6wIHMwaUNhJ4Temnr1tTURFtbW/GjE4l0dnayevXq7Yro3ub5jJtON2rUqO5iOn3aaaeddlg3rJh/okvZpReaZraikNcqZ4uI7Gjr1tDjvWBBmF59FebOheefD1PmyMu3vx0OzzyDJYt8c3a+Pd5rgFsI4wePymjbAtyU53b6YhHhzPd0I4DFhASfrU2kzzo7O1mzZg2rVq1i5cqV3fMVK1Z0T5mF9MqVK3H33jeeZtCgQTQ1NdHU1MS4ceMYN25c1sJ6p512YnDcp2NLpVDOFpGasHlzOKkydZLlsmXbbhqUvrxoURjbnetretw4OPJIeMtb4NRTYb/9ih9vvtfxftrMZgBnunspE3ZP7gduNLMh7r7ZzOqBKcCDQGeONpHukwwzC+jUvKd1q1atYvXq1QUX0RBu9DJu3LjuQrq3+fDhwzXMQ4pOOVtEKkHqGuHr14cb9GQur12b/ZrhqSnP05mAMBZ88uRwEuWUKWG+996wzz6w776h8C61Qq7jvZYsvSRm9u/u3sP5oYWzbVWIRe+71MzuBo4F7gaOA+5w90XR87O2SfVwdzZt2pR34Zze1tnZ2af3HDVqFGPHjmXMmDHbzdN7qNPnY8aMob6+vsh7LtI3ytkiUixbt8KmTbBxY2HzDRu2FdM9FdcFjrbsUUNDuIb32LEwfvy2q6FkLk+YEK7nHffXdG/X8f4t4RqtOZ8G7AfkedGXnO83AjgneniemX3H3duBi4BrzewgwmWw0s+wydUmCfT666/vUChnK57T5/lepSPTsGHDGDt2bI9FdOY8tdzY2KibuUjFUc4WqQ1dXfD666Fw3by55+VC2rMVz6nlrTuchVEcZuGW8iNGwMiR25ZTj0eO3PE285m3nx8ypDSxlUpvlcVcwljBRWy7BmwmA04pRjDuvh74bjSlr18NXJjlNVnbpLS2bNmSd9GcPi/0hMKUhoaGXgvozHVjxozRtaWllihni5RRZ+f2hWyhy/kUyD09d8uW3mMrpgEDwiX1hg0Lt6DPZ55aTi+oM4vrYcOKf53spOut8L4GWO/uObsazaxUZ8dLGXR0dLB69eqCiudVq1axYcOGPr1ffX191uI5V4/00KFDi7znIlVHOVtqTldXGOfbnwK4r8sFXOG16AYPDr29qXm25XzaeyuqBw2qvQK5VHq7gU6+l7A6jdKeJS956ujoYOXKldtd+7m9vT1r8bxq1SrWrl3bp/eqq6vrdehGT8X0sGHDdEKhSAkoZ0tc3EMvbF8K2P4WwH0chVgUZtuK19SU/ri35UIK5PTlhgYVwpWqX4NYzWwK4da//4mSeMls3Lixu4hOL6gzl5cvX86KFSsKvhrHgAEDaGxszKuITi+mR44cqQJapIIoZ1e/rVvL2+ObXjz34UJQRZNP4VtoUZzPciR09JkAACAASURBVH29CmApTMGFt5nVAacSxugdDywjjCmUArk7q1evZtGiRbz22mvdU+bjQnqkzaz7WtCpaz83NTXlHBs9atQoBgwYUMI9FZG4KGcn09at4coOGzZsmzZu3P5xoes2bQpjjuMyaFBhxWuxit9Bg8IYZJFKkHfhbWa7AR8GzgfGAq8C7wT+APxbSaKrAl1dXbz22mvMmzePl156iXnz5nVPL730Ehs3bux1Gw0NDUyYMKH77oSpgjrz8fjx4xk7dqyuxiEiytll8Prr227OsWoVrFmz/bR27Y7rUtOmTaWJaeDA4hS1hb6moQHq6kqzTyLVpNcKzczOJPSUvA2YB3wT+DFwmbv/Pnra30oWYQXZunUrTz75JP/85z+7pzlz5uQsrkeOHMnkyZPZeeedt5vS140ZM0ZDOkQkL8rZ/ffGG+F20q+8Em4pvWDB9nfAS03r1vX9PerqwlUdhg/veRo2LP/1qXVDh4bCW0SSK5//opsJtxj+G/Bhd38BwMxiHM2VHM8++yyzZs3i/vvv569//WuPRfb48ePZa6+92HPPPdljjz3Yc889u5cbGzPvniwi0i/K2Xlavx7mzNk2PfMMvPRSuK10Purrt92gY+xYGD2652nUqB3X1eJl1EQkj8Lb3e8G7o5OyrnQzJqBu4Due/+YmXlf7q9dodauXUtrayu33HILzzzzzHZtU6dO5bDDDuPggw/m4IMP5qCDDmJcOe5BKiKCcnYu7vCPf8Cvfw1/+hM8/njPY6Lr6sKtpJubYdddw62lJ07cVmSn7oQ3apSKZxEpTCG3jF8AfN7MBgLvAnY2sy8BPwBOAFpLE2JydHR08I1vfIOrrrqq+4THxsZGTjvtNI4//nimT5/OxIkTY45SREQ5O11XF/zsZ3DNNfD009vW19XBIYfAgQeG6YADYOrUcFtpDdkQkVIoOLW4ewfwK+BXZrYn8EnCeMKqTuJLly7l1FNP5fHHHwdg+vTpfOpTn+Kkk05i0KBBMUcnItKzWs3ZKcuWwZlnwl/+Eh6PGwfvfS+cdBIcfXS4k56ISLn06296d58HfMLMnihSPIm0cuVKjj76aObNm8eUKVO48cYbOeGEE+IOS0SkILWSs1PWr4djj4Vnn4UJE+ArX4FzzgmXnxMRiUNRfkxz91uKsZ2k+tCHPsS8efM4+OCD+eMf/0hTU1PcIYmI9Fm15+yUr3wlFN377gt//nMovkVE4qRLzvfigQce4K677mLEiBHcddddKrpFRCpAZyfceGNY/uEPVXSLSDKo8O7F9ddfD8CnP/1ppkyZEnM0IiKSj2efDTe1mTIFDj887mhERAIV3jmsW7eOu+++GzPjwgsvjDscERHJ00svhfn+++uSfyKSHCq8c3j44YfZsmULRx55JJMmTYo7HBERydNrr4W5fqgUkSRR4Z3DQw89BIRLB4qISOVYuTLMdVqOiCSJCu8cnngiXHHr3/7t32KORERECrFqVZiPGRNvHCIi6VR45/Dcc88BsP/++8cciYiIFCJVeI8dG28cIiLpVHhnsXbtWhYvXsyQIUPYdddd4w5HREQKoB5vEUmiotxAJ25mNgS4FlgPNAIz3X1df7b56quvAtDc3ExdXV2/YxQRkaAUOTtTaoy3Cm8RSZJq6fG+AXjY3S8D7gK+398NLlmyBEBXMxERKb6i5+xM6vEWkSSq+MLbzCYBZwOzolX3AaeaWXN/trt48WJAhbeISDGVKmdn2rAhzIcPL+ZWRUT6p+ILb2AGsMLdNwG4ewcwH+jXNQDV4y0iUhIzKEHOzqSb5ohIElXDGO/JwMqMdeuB7Srm9vZ2pk2b1v24paWFlpaWrBtdu3YtAKNHjy5WnCIieWttbaW1tTX1cFycsRRZSXJ2plTh3dXVxyhFRAqQb86uhsLbgc0Z6xqArekrmpqaaGtry3ujr7/+OgBDhgzpZ3giIoVLLzTNbEXM4RRTSXJ2pgkTYPFimDsXdtmlz5sREclLvjm7GoaaLCKcFZ9uBLC4PxvdsmULAPX19f3ZjIiIbK8kOTvTsceG+e9/X8ytioj0TzUU3vcDk6LLU2Fm9cAU4MFibNw0UFBEpJhKmrNT3vnOML/jDnAv5pZFRPqu4gtvd18K3A1E/RscB9zh7ovii0pERHpSrpx91FEwcSK88gr0Y8SKiEhRVXzhHbkI+A8zuxx4B5D/GTgiIlJuJc/ZAwbAmWeG5V/+sthbFxHpm2o4uRJ3Xw1cGHccIiLSu3Ll7DPPhG9/G37xC7j6atBNiEUkbtXS4y0iIrKdt7wFdtsNFi6E++6LOxoRERXeIiJSpQYMgA99KCzfdFO8sYiIgApvERGpYuefHwrwu+6C5cvjjkZEap0KbxERqVqTJsHJJ8PWrfCDH8QdjYjUOhXeIiJS1S6+OMyvvx6ie6OJiMRChbeIiFS1E06A/feHJUvgZz+LOxoRqWUqvEVEpKqZwaWXhuVvfEN3shSR+KjwFhGRqve+94U7WT79NNx9d9zRiEitUuEtIiJVb9AgmDkzLF9xBXR1xRuPiNQmFd4iIlITPvIRmDwZ/vlPuPPOuKMRkVqkwltERGrC4MHwuc+F5f/9X+jsjDceEak9KrxFRKRmXHABNDfDc8/Bj34UdzQiUmtUeIuISM0YNAiuuiosX3YZrF0bbzwiUltUeIuISE15z3vgqKPCLeS//OW4oxGRWqLCW0REaooZfOtb2+Zz58YdkYjUChXeIiJSc9785jDee+tW+OhHdVMdESkPFd4iIlKTrrkGmprggQfgppvijkZEaoEKbxERqUljx8K3vx2WL70UFi2KNx4RqX4qvEVEpGa95z1wyimwbh187GMaciIipaXCW0REapYZfPe7MHIk/Pa38IMfxB2RiFSzRBTeFpxlZv/qoW2cmX3fzK4xs2+aWX0+bSIiUjrVlLd33hluuCEsf+ITusqJiJROIgpvYFegEdinh7bbgZvdfSbwAvDFPNtERKR0qipvv//9Ydq0Kcy3bIk7IhGpRokovN39FWB25nozOxKY6u6PRqtmAZeY2bBcbeWIWUSkllVj3v6//wu3k3/iiXBXSxGRYktE4R3p6mHd24D5qQfunjrnfFovbSIiUnpVlbdHjYKf/hQGDoRvfAN+9au4IxKRapOkwrsnk4GVGevWA5N6adtBe3s706ZN655aW1uLHqyISLG0trZ25ytgXNzxFKAoeTuunP2Wt8DXvx6Wzz8f/rXDCHYRkR3lm7MHliMYM7sKOCBL833ufl2WNgc2Z6xrALb20raDpqYm2tra8gtYRCRmLS0ttLS0AGBmK8r9/nHn7Thz9sc/Do8+Cj/7GZx+Ojz2GIwYEUsoIlIh8s3ZZSm83f1/+vjSRcDeqQdmZsAwYHEvbSIi0g+1nLfN4MYb4emn4ZlnwsmWv/411NXFHZmIVLqkDzWZDUxNe9wMbAHaemkTEZF4VEXeHjYM7rwTGhth1qxwZ0sRkf5KUuFt0N0DAoC7twFLzGz/aNVJwPXuviVXWzmDFhGpYVWdt/faK/R019fDdddtu9a3iEhfJaLwNrOdgHOjhxeZ2fC05vcAnzSzzxN6R67Is01EREqkVvL29Olw001h+ZJL4N57441HRCpbWcZ498bdlwNXRlNm26vAh7K8LmubiIiUTi3l7Q9+EF56Cb74RTjrLPjzn+Hww+OOSkQqUSJ6vEVERJLsyivhAx+AjRvhHe/QZQZFpG9UeIuIiPTCDH7wAzjlFFi5Eo4/Hl59Ne6oRKTSqPAWERHJQ309/PKXcPTRsGgRnHACvPFG3FGJSCVR4Z3FgQceyKmnnkpzc3PcoYiISEIMGRIuL3jYYXD55dDQEHdEIlJJEnFyZRJddNFFXHTRRXGHISIiCTNqFPz977qhjogUTj3eIiIiBVLRLSJ9ocJbRERERKQMVHiLiIiIiJSBCm8RERERkTJQ4Z1Fa2tr3CGURa3sJ2hfq1Gt7KfkT58JHQPQMUjRcUjeMVDhnUXS/qFKpVb2E7Sv1ahW9lPyp8+EjgHoGKToOCTvGKjwFhEREREpA3P3uGMoCzNrBwq5we84YEWJwkmSWtlP0L5Wo1rZz13dvSnuIMqpDzk7pVY+E7noGOgYpOg4xHMMsubsmim8RURERETipKEmIiIiIiJloMJbRERERKQMVHhnMLMhZna9mX3VzL5rZiPjjqmYzOx4M3vUzJrT1lXVPpvZ6WY238xWmtm3zGxgtH6cmX3fzK4xs2+aWX3csfaXmU03s7lmttrMrktbX1X/punMbJCZPWVmM6LHVbuvkr9a/hyY2ZfNzKPpqWhdVR+PQr/LqvV49HQcovU7fCai9VV1HPryfR93LaDCe0c3AA+7+2XAXcD3Y46naMxsPDAcODyjqWr22cymAO8CzgA+CXwY+K+o+XbgZnefCbwAfDGWIIvEzEYBRwFHAucAF5vZcVFz1fyb9uAzQHPa42reV8lfTX4OzGwY4eSx46PpjKipao9HH7/Lqu54ZDsOOT4TUEXHoR/f97HWAjq5Mo2ZTQLmA43uvin6y2ktsJ+7vxJrcEViZgOATmA3d3+l2vbZzI4G/u7uHdHja4H9gC8Bd7r7pGj9ZGAuMN7dN8YVb3+Y2RDgdY/+E5vZE8ClhP2qmn/TdGb2FmAf4ArgPELSrMp9lfxVWx4rhJmlCo3vufvr0bqqPx6FfJcBW7K1VfrxyDwO0bodPhPR+qr6XPTl+x44IFtbuWoB9Xhvbwawwt03AUT/mPOB6XEGVUzu3pWxagZVtM/u/pfUf8LIImAB8DbCfqWetyhanFbG8IrK3TenFd3DgH8BD1Bl/6Yp0T6e4e43p62eQRXuqxRsBjX4OTAzA94LfANYZmbvi5pmUOXHo8DvslxtFS3zOOT4TECVHYc+ft/HXguo8N7eZGBlxrr1wKQYYimXat/nw4DvUcX7Gf3V/zugHhhK9e7rTODqjHXVuq9SmJr8HHhwJGFYwfXAT83s7dTm8ci1zzVzPHJ8JqD6j0M+3/exH4OB5XqjCuHA5ox1DcDWGGIpl6rdZzPbE2h396fMrGr3E5gH/Ai4Fvga8DJVtq9mdiLQ5u7LM5qq+d9V8lfTnwN3Xw18Lurt/C9gNrV3PHr7DNTU8ejhM/E7qvj/SQHf97EfA/V4b28R0JixbgSwOIZYyqUq99nM6oCPAJ+NVm23n1EyGkaF7yeAuy+Jhl9cSvjJsBr/TS8FbjazFWa2AtiFcGLQMKpvX6Vw1fiZ74tvA1OozeORa59r8XikpD4TUKXHocDv+9hrARXe27sfmBSdtEZ0iZkpwIOxRlVa1brPlwJfd/c3osezgalp7c2EE27ayhxXKbURkko1/pueDRycNi0mnMH+Q6pvX6Vw1fiZ74su4B/U5vHItc+1eDxSUp8JqN7jUMj3fey1gArvNO6+FLgbODZadRxwR9rg+4oX/XUHYFCd+2xmlwFPAkPMbHczuwBYAywxs/2jp50EXO/uW+KKs7/MbKiZHZK26h3At6rx39Td2939tdREOIu/3d0XUGX7KoWrxs98PsxskpmdbWYDotz+aeDyWjgehXyXVfPxyDwO2T4TUJ3/Twr9vnf3tmxtZYtZlxPcnpk1EsbKvgJMAP7H3TfEGlSRmNkIwvWebyBct/I77t5eTftsZp8jXEoo3fPuvq+Z7Uq4DN0rhGufXp5xRnRFMbPDgXuA54FHgMfc/faorWr+TXtiZq8A57n7A9W+r5KfWvwcmNl+hLG764CHgevc/fmorWqPR1++y6rxePR0HICdyPKZiF5TNcehr9/3cdcCKrxFRERERMpAQ01ERERERMpAhbeIiIiISBmo8BYRERERKQMV3iIiIiIiZaDCW0RERESkDFR4i4iIiIiUgQpvkT4ys2PN7G9mNqOX5+1vZr80s88X8b2/VqxtiYjEwcwONLPbi5kbs7zP3Wb23lK+R6mY2XFm9vfevmekcgyMOwCR/jKzfYGvAKuAw4ADgQ+4+09K/NYLgYPyeN4CoBGoK8abRrf53VqMbYmIAJjZYcBVwFuBPwLzCTdYqQe+6+5/KMHbzgNGUaTcmMMNwDOl2HCUj08Hvge8APw9amoA3g0c7u6v9OMt8v2ekQqhwluqwW+As9z9KTMbAPywHG/q7i+Y2co8nrfOzJYU8a2nAw8VcXsiUuPc/XEz+zmwj7ufklpvZicDvzazS939hiK/56Yi50YAzOydwNOpgtfd7y32e6S4+1bgF2Z2CfCAu38uLY67+rJNM9sJOMbdf+Huc81sRZHClQTQUBOpaGY2Btgb2Ajg7l2EW8GW65asXUV+Xj6OAe4v4vZERAB2uG22u98DXAn8PzObUIL3LGquNrOdCT3c5dbTLcfvAxYVshEzawBuAYakrS7m94fETIW3VDR3XwXMBX5vZkdG614FHjOz0Wb2JTN7zMxONbOlZvYvMzs09Xoze5+ZXW5mD5jZV9PWTzOzK83sF2Z2p5mNitYPNrP/M7MrzOwHhJ9JC2Jm/x3FdJaZHWZmv4m296Uoxl+YWVO0fq2Z/W/GJga5+xvRtj5tZp8ys7+a2dUFH0ARkd79gDDk5EwAM9vDzL5oZjeZ2Z+iYhcza4hy2eVm9hcze09qA2b2MTP7vJk9Ymb/me2NesrJZnZGlDPPM7MXzexqM5tkZj+M8ulfzOxt0SbOBnYGZprZu6Jzce4zs3PT3uNCM/uKmd1iZj82s2FmNsrM/tfM2szsbWY218xeNrPmvhyw6JjMiHrEMbPTzexbZvYNM/u9mQ2P9uFbZjYr+r5ZEh3jNwFnmtmn0jbZGH0XrTezC7O85wAzu8bMWqL3cDP7derfRxLC3TVpqugJ2A94ldArcDOwU7S+jpDEVhHG2o0GZgPPEv7o/Dfgc9FzGwm95icDI4Gfp23/r8DXouX/B7REy+OBTkJy7S3GHxF6jXYijAUcGq1vAP4A3AVMInxhbAX+FxgKvDl6nHp+avx6ar9viJZHAJ+N+99CkyZNlTsB5wGvZWlbBnw7yqt3AnXR+tuA26PlbwKnR8sXAO3R8vuAc6LlvaNcfUD0+EfAldFytpw8BlgDXAccDhwHfB24InruZ4BZabE60BwtDwH+BZwXPT4RuC3tuXcD3ycMvT0t+r44PdrP3wNfyfPYPRB9V1wNfI0w3vu4tPYVwJRo+VnCd1IDcDnwGnAE8KEo7z+Qijd6/ivAN6K29wELssTwYeDuaHl8dJxPiPtzpWn7SWO8peK5+7NmdgDwZeBjwDvN7DgPY77bgXXufgeAmV0JPAzsTviS6TKz/4o29QdCgj8ZGJu2/l/Ra4cDF0Wvxd2XmdnCAkIdB1wPfMjdN0XbeMPMlgLz3X1x9D5LgQc9jH98kvCFMI5wkubJhC8JgE3AB83sBeC7hD86RERKwQnF6OHArsAlZgawEhhl4cFHCJ0GAD8mDLWAkGvbMnLtRODpjPc4jx5ysruvMrO1wF3u/hiAmS0ANprZREIP8fAeg3bfbGbL01b9Z7TdlJuAX0br1xC+L+6M3qMtijNf93s0xtvMfkHo7Ek5090XmNlRhIJ7eJT/FwMvufujwKPRa3va9qzoO+HxHDEdCmyO9nuZma1i+yErkgAqvKWimVkdsIuHk2g+Hg3/uBdoJfQgZHo2mo8GpgC3uvtt0brrom3OJBTC12W815sJCXN9H8PdG/h3YA/gqbT1mWMcO7sb3LuiJJwaFjbWw/Aa3H2+mZ0HfAf4JPB+YDkiIkVkZqMJv9Y9S8ibK3vIj+MJPbKDANy9k9BZQPSaL7j736LH2702TY85OeJsnytfA2YCLwKPALvl2IX01+0F/Dnt8UuEYTQ7sWMu7qDvV1z5B+EPlJQ1ZnYNcDsh9lR1nblfvekge+32O+C7ZjaY8G+xmdB7LgmiMd5S6brHHQK4+1PAZcA+WZ4/iJDk5gFLgDPSG83siGj9yWY2JGP9uujhvn2M9W/At4CfRomxIGY2ibQTdcxsMuGKLnsTTra8LctLRUT641zCL2y3E/LjUVFPM9CdH1cAW4CT0tYfYGYjyMi1ZjYw6sjIlC0n9+R7wCseLhv7RgH7soDQQ979FsBawlCaonH3rqhzZL/oD5d7gS+6e1sx3yfjPWcBswi/zL4PmO7ua0v1ftI3KrylGlwSFaUpkwlj81LGmNmgaHkGYXzfGuBnwGlm9h0zmxH1RmwkJMhhwN1mdmJ0Iss+hGL9SeBLZjbIzMYRxlbvZGZ1ZvZZMzsnS4wWTVcQvpy+2UNb5vMzH59CGI+YshfwLndfB1zSw2tERAqxQ0+qmb0V+ALwEXdvJ/QuLwPuNbN3RidQHh31cP+ccPWTd5nZscC57r6ekGs/bmZfMLPphA6I1DC99PyXLSenpPc+HwI0mdlQ4ChgiJntHrVtJZyMmOqASX+P7wHvNrPG6PERwI3u3kGoiXrMxWZ2ZBRXQ5ZjN4CMmir64+Rkwq+cjYQhjFOBXYChaSduZvaqb+kh/u2D6mE8SvRvMRL4E6FXv5A/SKRMNNREqsEuwLNmdg/hZzgDPprW3gFcbmZrCL3DlwC4+2wz+zjh58p3EU5OfAbAzP4D+D/CF8nNwGfc3c3sLMKlnp4FfkH4yfBgwg0nziIMYdnuxj1mdjjhpKEphJ8efwN8IYrnD8CRQLOZ7UHoiZkAnGVmz7Ot9+c8YJi7t2bs+w8t3Piig3Ayk4hIwaI8cjahOLwRaCecoDcKeLu7/x3A3bdE+fFGQq67izC2G+AThE6LWwi/8KWuJPIDQu67iDAk7iJ3Xx4VlkcCe5jZT7LlZDM7gzCu+cNm9ry7LyKc6Pl1QtF9C/AfwAHAy8CthN75D1q42tV+wDvM7F53v9PM9gLuMLO/En41vTzqmT8TmBDt37OEmwlNNLP9ojjfSxi3/pu041ZPyP0HAZOjwtwJY6tPj7b5T+BB4HHC8JnfEXL63YQTOg8ys3e5e2q7vyDczGhDVKhPAM4ws/mEnmyi12fes2IB4eTRdwGDgQFmdh/hBMtyXWJXemH6t5BqZuE2uz9y9+Yyvd9n3F23cxcRqTLRsJcGd0/kDczM7DOEK8y8Ej0eTPi14prUuUESPw01kWrX0zCO0rxR+En2H+V4LxERKZ/o5NFDE1x0DwY+BzSnrR4HbFTRnSwaaiJVKxpfdw7hp8Mz3f32Er/lS9FPoCIiUl0Guvt34w4iG3d/3cJt678dnXi/mDAs5nO5XynlpqEmIiIiIiJloKEmIiIiIiJloMJbRERERKQMVHiLiIiIiJRBzZxcOW7cOG9ubi7oNe3t7TQ1NZUmoAqm47IjHZMd6Zj0rC/H5Yknnljh7jV1MPuSs6G2Pnfa1+qkfa18uXJ2zRTezc3NtLUVdqfWadOmFfyaWqDjsiMdkx3pmPSsL8fFzF4tUTiJ1ZecDbX1udO+Vifta+XLlbM11EREREREpAxUeIuIiIiIlIEK7xxaWlriDiGRdFx2pGOyIx2Tnum4lFYtHV/ta3XSvla3mrmBzrRp07waxxGJSPUzsyfcfVrccZSTcraIVKpcOVs93iIiIiIiZaDCW0RERESkDFR4Z/O3v8Ghh8LHPhZ3JCIi0ouV966k7dA2Xr785bhDERHJqmau412wDRvgySdh7Ni4IxERkV50rOlgw5MbGDp1aNyhiIhkpR7vbAZGf5N0dsYbh4iI9MrqDQDvqI0LBohIZVLhnU2q8O7oiDcOERHplQ0MhXfX1q6YIxERyU6FdzYqvEVEKkZ3j/dW9XiLSHKp8M6mri7MNdRERCTxBtSHrzMNNRGRJFPhnY16vEVEKkZqqIl6vEUkyVR4Z6PCW0SkYujkShGpBCq8s1HhLSJSMdTjLSKVQIV3NrqcoIhIxUj1eOuqJiKSZCq8s0mdXKkebxGRxOvu8dZQExFJMBXe2WioiYhIxdDlBEWkEiTulvFmNgh4HPiEuz9gZkOAa4H1QCMw093XRc/N2tZvKrxFRPKShLytywmKSCVIYo/3Z4DmtMc3AA+7+2XAXcD382zrH43xFhHJV+x5WydXikglSFThbWZvAZYAq6PHk4CzgVnRU+4DTjWz5lxtRQlGY7xFRHqVlLytoSYiUgkSU3ib2TDgDHe/OW31DGCFu28CcPcOYD4wvZe2/tNQExGRnJKUt3VypYhUgiSN8Z4JXJ2xbjKwMmPdemAS0JGjrf9UeIuI9CYxeVuXExSRSpCIwtvMTgTa3H15RpMDmzPWNQBbo+Vcbdtpb29n2rRp3Y9bWlpoaWnJHpQKbxGJUWtrK62tramH4+KMpSelztuF5mwNNRGROOWbsxNReAOXAoeYWepxI+Gkm29Gy+lGAIsJyT1b2w6amppoa2vLP6JBg8L8jTfAHbbFJiJScumFppmtiDmcnpQ0bxeaszXURETilG/OTsoY77OBg9OmxcCHgR8Ck6LLT2Fm9cAU4EHg/hxt/VdXBw0NoejesqUomxQRqSKJytvdlxNUj7eIJFgiCm93b3f311IT0Am0u/sC4G7g2OipxwF3uPsid1+ara1ogQ0ZEuabNhVtkyIi1SBpeTv9coLuKr5FJJmSMtQkl4uAa83sIGAC0JJnW/8NHQpr1oTCuzHz11EREcmi7Hnb6gyrt1B4b3GsQcMDRSR5Ell4u3tz2vJq4MIsz8vaVhTq8RYRyUsS8vaAIQPo3NpJ5+ZOBjQk4gddEZHtKDPlMnRomG/OPAlfRESSpm5ouPFZ1yZdUlBEkkmFdy6pwls93iIiiTdgSPhK69qswltEkkmFdy4aaiIiUjFShXfn5s6YIxER6ZkK71zU4y0iUjG6h5qox1tEEkqFdy4jRoT5unXxxiEiIr3qHmqiMd4iklAqvHNJXUJw9ep44xARkV5pqImIJJ0K71xUeIuIVAydXCkiSafCO5dU4b1qVbxxiIhIrzTGW0SSToV3LmPGhLl6vEVEEm/A0GioyUYNNRGRZFLhILFW/wAAIABJREFUnYuGmoiIVIyBI8PNmDvXqfAWkWRS4Z2LCm8RkYoxcHQovDvWdMQciYhIz1R45zJ2bJi3t8cbh4iI9CpVeG9dvTXmSEREeqbCO5cJE8J86dJ44xARkV4NbFSPt4gkmwrvXMaNg7q6MNTkjTfijkZERHLQUBMRSToV3rkMGADjx4flZcvijUVERHLqLrxXq/AWkWRS4d0bDTcREakI6vEWkaRT4d2biRPDfMmSeOMQEZGcugvvVSq8RSSZVHj3Rj3eIiIVYdBOgwDY0r4F7/KYoxER2ZEK796o8BYRqQgDGgaEK5t0wtaVuqSgiCSPCu/epArvxYvjjUNERHo1aGLU671kS8yRiIjsKFGFt5lNN7O5ZrbazK5LWz/EzK43s6+a2XfNbGQ+bUXR3Bzm8+cXdbMiIpUuiTm7u/BeqsJbRJInMYW3mY0CjgKOBM4BLjaz46LmG4CH3f0y4C7g+2kvzdXWf3vuGebz5hV1syIilSypOXvQBPV4i0hyJabwBrYAV7n7ane/B5gDdJrZJOBsYFb0vPuAU82sOVdb0aLabTcwg1dfhS1K5CIikUTm7IaJDSE49XiLSAIlpvB2983u7gBmNgz4F/AAMANY4e6boud1APOB6b20FUdDA+y+O3R1wXPPFW2zIiKVLKk5OzXU5I1FutuwiCRPYgrvFDM7GvgdUA8MBSYDKzOeth6Y1Etb8Rx+eJg/+mhRNysiUumSlrMH7zYYgM0vby7WJkVEiiZxhTcwD/gRcCzwNcCBzAzaAGztpW077e3tTJs2rXtqbW3NP6IjjgjzRx7J/zUiIv3Q2trana+AcXHHk0OicvaQPYcAsHmeCm8RKZ98c/bA8oWUH3dfAtxsZl3AZ4C/AI0ZTxsBLCYk8Wxt22lqaqKtra1vQb31rWE+eza4hzHfIiIl1NLSQktLCwD2/9u78zi5qjL/45+nl+o9SYd0ls5GAhiQNRIRR2UTN9wGBXV0HHWEuI/bzA8RxnFcRhR1ZFxG4jLujjPiDIrIDCKCiKMkoFHWQBZCErIvnXTS6/P749xKV1dXVVd139q6vu/X67zq1j1V9557u/rpp0+de67ZrjI3J6tKi9ktS0PifWT9EXzIsXrFaxEpvnxjdiX2eCetBrYAtwPdZtYCYGaNwCLgjnHq4nPGGWE+7y1b4I9/jHXTIiJTREXE7Pq2ehLdCXzAObL5SFybFRGJRcUk3mbWambLU1ZdBFzn7k8CNxG+xgS4ELjB3bfkqou5cfCiF4Xlm2+OddMiItWokmN267JWAHrv741zsyIik1YxiTdwCvC/ZvYrM7sW2BBNUQXwNuDlZnYVIbivTHlfrrr4vOQl4fGHPyzK5kVEqkzFxuz25e0A9NzbE/emRUQmpWLGeLv774CuLHV7gcsLrYvVRRfBtGmwZg08/DAsW1b0XYqIVKpKjtkdT+sA4OC9B4u5GxGRglVSj3dla26GV74yLH/ve+Vti4iIZNX+tKjHe416vEWksijxLsRrXxsev/3tcEMdERGpOK1PaaV+Wj19m/s48rgusBSRyqHEuxDnnw+LF8OGDfA//1Pu1oiISAZWb8w4dwYA+27fV+bWiIiMiCXxNrOL4thOxauvh7e+NSx/8YvlbYuIyATVQsyecUFIvPfevrfMLRERGVFQ4m1mLzSzu8zsUTNbH5WNwA3FaV4FevObIZEI0wquX1/u1oiIZFXLMbvzgnCfnr237sWHvcytEREJCp3V5NvAJ4D7gOQgZwNeG2ejKlpXF/zFX8A3vwnXXgv/+q/lbpGISDY1G7PbTm2jaWETfZv7OPC7A0w/e3q5myQiUnDivd7dP5u+0sweiqk91eGKK+Bb34KvfQ0++EFYuLDcLRIRyaRmY7aZMesVs9hy3RZ23bBLibeIVIRCx3i/zcwuybD+ZXE0pmqcdBK86lUwMACf/GS5WyMikk1Nx+yuV4Zpxnf+cCfuGm4iIuVXaOL9LOCrZjaUUoaB2htvcfXV4Vby118PDz5Y7taIiGRS0zF7+p9Np2lBE0c2HtHsJiJSEQpNvD9CuL3vCcDSlPLheJtVBU45BS67DAYH4R3vAPWmiEjlqemYbfXG3DfPBWDrqq1lbo2ISOGJ98+AG919vbtvispG4CvxN60KfOITcMwxcPvtushSRCpRzcfseZfNgzrY9aNd9G3rK3dzRKTGFZp4/ztwpZmdk1LOBa4sQtsq3zHHwJe+FJbf9z5Yu7a87RERGa3mY3bzgmZm/fksfMDZ/JnN5W6OiNQ4K+SCEzO7h/CV5T4g+UYD5rp7c/zNi8+KFSt89erVxdn45ZfDV78KS5fCb34Ds2cXZz8iUpPMbI27r5jA+xSzgZ57e1hz5hrqWus4e+PZJLoSsWxXRCSTXDG70B7va4GF7n6suy+JyrHAKybbyKp23XVw5pnhhjovexn09pa7RSIioJgNQMfTOph50UyGe4fZ9PFN5W6OiNSwQhPvze7ek2H93DgaU7VaW+Gmm2DxYvjtb+GlL4VDh8rdKhERxezIko8vgTrY8oUtHPzTwXI3R0RqVKGJ99vTV5jZpcDn4mlOFZs7F265JTz+4hdw0UXQk+nvnYhIyShmRzrO6KD7rd0wBOveuU63kReRsig08Z5pZn8PYGatZvZV4KvAPbG3rBqdeCL88pfQ3Q133gnPehY8/ni5WyUitUsxO8WSjy6hsauR/Xfs54nrnih3c0SkBhWUeLv7i4FbzezzwO+BpwCnA5nujFabli2DX/0qJOF//COcdRbcfXe5WyUiNUgxe7TGmY0s+9oyANZ/YD09v9e3kiJSWjkT77QpqM4xs3OABLAJWAf8A7AY+FDxm1pFli4NyfZznwvbt8M558A//RMMDZW7ZSIyhSlmj2/WS2fR/dZuvN+5/+L76d/ZX+4miUgNyTmdoJn9H3ASYSqq4WwvA+a5e1P8zYtPUacTzGZgAD74Qfj0p8Pzc88Nt5hftqy07RCRqpbvdIKK2fkZOjzE78/9PT339DDtWdM4/eenU99cX5R9iUjtmcx0gh8BTnD3xSlTUaWXY4HXxdDIV5jZBjPbbWbXmVlDtH6WmV1vZp80s8+aWWPKe7LWVYTGRrj2Wrj5ZujqgjvugNNOg6uugv37y906EZl6ShazoXrjdn1LPafceAqJ+QkO/PoA919yP8N92f5PERGJT87E291vdvcd423E3X84mUaY2SLgzwnjDt8LXAa8J6r+T+Dr7n4F8AjhDwt51FWOF70IHngA3vxm6O8Pw06WLIFrrtHMJyISm1LFbKj+uN00r4nTfnYaDcc0sOene7j/VfczdETDAUWkuAqd1aRYFgN/7e5r3P1bwBeB883sbGCZu/82et1PgHeZWVuuupK3Ph+zZoW7W/7612HM9969cOWVsGABvOc98Oij5W6hiEghqj5ut5/azuk/P52GzgZ2/3g3a5+3loHdA+VoiojUiIpIvN39V+4+mLJqC/A4cAGwIeV1W6LFFePUVa4/+7Mw5eCtt8JzngMHDoQ7X55wQrgY8+tf1zAUEal4UyVud5zRwRm/PIPE/AT779rPvWffq9lORKRoKiLxzuDpwJeB+cDutLoeoHucuspmBhdeGOb6vvdeeNOboKkp3HjnzW+GOXPgkkvgO9+B3emHKCJSkao2bref1s6Zvz2T9jPaOfzoYe49+162fHELuSYfEBGZiEkn3ma23My+Y2bvN7NjY9je8cBOd/8D4MDhtJc0AQPj1I2xc+dOVqxYcbSsWrVqsk2Nx/LloZd72zb4ylfg/PPDOPAbboDXvx5mz4ZnPzuMB7/vPhjWBUAitWDVqlVH4xUwK67txh2zo23GHrdLHbOb5jex/O7lzFs5D+9z1r1zHX943h/oXddb1P2KyNSQb8zOOZ1gPszseuDl7j7XzF4L/Je7pwfWfLdVD1wDXO3ufWZ2JXC+uz8/qjfgCHA+cG62Oncfc8easkwnOFFPPBES75/+NAxLGUj5m9TZGaYlPP/8UE4+Geoq9YsLEYlDvtMJ5rmt2GJ2tL2ixO1yxuwdP9jBI+94hMHdg1iTsej/LWLh3y6kYVpDWdojItUlV8yOI/FuBGa7+xYzqwPa3f3ABLd1BfANd98ePV8B3ODui6PnS4C1wDHAadnq3H3MHRGqKvFO1dMTxoP/9Kdw222wadPo+s5OeMYzQjn77HCnzJkzy9NWESmKmBPv2GJ2tL2ixO1yx+z+nf089nePsf2b2wFonNXIoqsW0f2WbupbNOe3iGSXK2YX/O+7mb0ImOnu3zWzs4CD7v4AgLsPAxNNuj8IrAZazGwpcB5wJ7DNzE5x9z8BLwQ+HwXo1WaWrW7q6OiAV7wiFIANG0Iv+O23h/LEE3DLLaEkPeUpIQlfsQLOOANOPx2mTStL80WkvIoVs6NtT9m4nehKcNI3TqL78m4eu+IxDvz6AI+99zEe/6fHmf+u+XS/rZvErES5mykiVaagHm8z+zjwbuAWd78kWncNsNbdvzfhRphdDXw0bfVD7n6SmSVvb7wRaAeuSl5Jn6suXbl7T4rCPSTe//d/8Nvfhsc1a+DIkbGvPe64kISnlvnzw4WeIlLRJtrjXayYHW2nqHG7kmK2u7P7p7vZ+A8bOXjvQQDqWuqY87o5zLtsHh1ndWCKpSISiW2oiZndDlwMvMHdr4vWnQt8290XxdHYYqmkIF5U/f2wdm1IxO+7D37/e/jjH8P6dDNmwFOfOrqcfLIScpEKM4nEWzE7Ru7Ovl/uY/O1m9nzsz1H17ed0sbcN89l9mtm0zS3qYwtFJFKEOdQk1vdfZ+ZpWbrl0y8aRK7RCIMMVmR8vMeGICHHgpJeGrZswfuvjuUVB0dIQk/6aQwbOX448M848cfD+3tpT0eEZkMxewYmRmd53fSeX4nhx46xJNfe5Inv/kkh/50iMfe+xiPve8xpp8zndmXzmbWK2cpCReRMQrt8X4NsJQwF+uvgEsJvSnvcvcvFqWFManE3pOycocdO8Kt7O+/Pzwml3ftyv6+uXNHJ+LJxyVLQg+6iMRuEj3eitlFNtw/zO6f7ObJbz7Jnv/Zg/dHf1MNpj1zGjNfNJNjXnQM7cvbsTp9kyhSC2Kd1SS65e8bgUXANuD77v7zyTay2KoliFeEnTtDEv7gg+FW9uvWhcfHHoO+vuzvmzYNFi+GY48d/ZhcnjVLQ1hEJmAys5ooZpfO4P5Bdv1kFzv/cyd7bklJwoHG2Y3MfMFMOp/fyYxzZ9C8sLmMLRWRYir2dIL1wHnuftukNlRk1RjEK87QEGzZMpKIJx8ffRQ2boRDh3K/v7UVFi0Kifj8+bBgQXhMLUrORcaIeTpBxewSGOwZZO9te9nzsz3s+dke+jaP7rRoPraZ6edOZ8Y5M0IivrRZF2iKTBETTrzN7MeEq85zmQXscvcLJt7E4qv2IF7x3MOY8U2bQhK+adPo5Y0bYd++8bfT1ATd3WMT8gULwjCXZOnoUIIuNSPfxFsxuzK5O70P9rLnZ3vYe/te9t+1n6H9Q6Ne0zi7kWlnTaPj6R10nNXBtKdPo/GYxjK1WEQmYzIXV66NXvMwMAxcDtxHmLc1aQmgO7bUOjM45phQnva0zK85cCAk4Y8/HnrOn3giPKaWvXvDXOUbNuTeX0sLzJkzOhlPfZ663NIS//GKVCbF7ApkZrQ9tY22p7ax8P0L8SHn4NqD7L9zP/vu3Mf+O/czsGOA3TftZvdNu4++r3lJMx1nddB+Rjvtp7XTdlobTfOb1DMuUsXGS7yvJfSK7wMws+Pc/UPpLzKzG4vROJlipk2DU08NJZve3rHJeDJB374dnnwStm0Lr9u4MZR89tvVFYax5FM6O6Fed6aTqqSYXQWs3uhY3kHH8g4WvHsB7s6RDUc48LsD9Pyuh557euhZ08ORDUc4suEIO3+w8+h7GzobaDut7Wgi3n5qO60nttIwXbezF6kGOX9T3X1/2qo5Zlbv7ke/IzOz84CzitA2qUWtrWGmlBNOyP26gwdDEp5MxpMl9Xly+cCBUB57LL82mMHMmaH3PjUZnzFj5DHbcnu7hsBI2ShmVyczo2VpCy1LW5jzmjkADA8O0/tALz339HBw7UEOrT3EwbUHGdwzyP479rP/jtE/6sTcBC3LWmhd1nq0tCxrofnYZuoa6spxWCKSQaH/Iv8Q+IOZ3QX0ACcBLyDchUykdNrbwzSGxx+f+3XuYfjKrl2Zy+7dY9ft2RPW794NjzxSWLvq68cm46kJ+owZoQc+WTo6xi63tSl5l7goZlepuoY62k8LQ0yS3J3+rf0jifgfD3Loj4c4/Mhh+p/sp//J/jEJuTUaLce30Ly0mZYlLTQvaT5aWpa0qKdcpMQK+o1z91vN7PnA24HTgH3AJe6ury2lMiV7r2fODDcDysfg4Nhkfd++sC7TY+pyb+9I0j6ZNieT8EyJebbl9vaR0tY2spxITLwtUtUUs6cWM6NpfhNN85s45kXHHF3vw07f5j56H+49Wg4/fJjeh3rpe6KP3gd76X2wN+M2GzobRiflxzaHfSwIpbGrUfOPi8So4H913X0rcHXyuZnVmdmZ7r4m1paJlEtDQxgT3tVV+Hv7+mD//sxJefKxpycMfUk+pi739ITkPbk+ruPJlJBPdlkJfVVQzJ76rM5oXtxM8+JmZj5/9HWzQ4eG6F3Xe3S8+JENRzi84fDR5cG9gxxcc5CDaw5m3najkehOHE3EU5Py5D8BibkJ6po0nEUkHwUl3mb29Qyru4Ah4M9jaZFINWtqgtmzQ5mowcGQgI+XoKevO3gwlEOHRpYPHgzbS/4TEKfGxuzJ+WQem5o01CYmitlS31ZPxxkddJzRMabO3RnYMTCSiK8/wpHNR+h7ou9oGdw9SN+mPvo25bh5GtAwo4HE3ASNcxpJzE2QmBOVuSOPjXMaScxOUJdQki61q9Ae7xOA24DUyb8XEG5FLCJxaGgI48I7O+PZXn//2IQ8juWBgdCLv3dvPO1Mqq+PJ4FPf2xtrcWEXjFbsjKzowny9LOnZ3zN0OEh+rf2j0rG+57oo2/LyHL/jn4G9w0yuG8QHhp/vw0zG47ut3FW49HScEzDqOfJUt9WrykUZcooNPG+zN0fTl1hZicSbkcsIpUokRgZ5x4X95DQJ5PwQ4dCz/uhQ6PXTeSxvz/eoTZJZiH5zpWgT58eLoCdPj33cktLtSTxitkyKfUt9bQc10LLcdnvh+DDzsDuAfq39zOwfSBc6Lk9Kk+OPA5sHwhJ+p5BBvcMZh13ns6aLGNCPqrMbKRhZgONneGxYXqDxqZLRSr04sqHM6zuAy4DPhBLi0Sk8pmFISFNTfEm9BB60uNI4NMfDx8e2e5kNTSMJOJdXSPDi9LLwoWwZEkYllMGitlSClZnJLoSJLoScEru1/pQSpK+c4CBXeOUnQMMHxmmf0s//Vv6C2hUGP7S0NkQkvLOhlGJeaZ1ydfWtdSph12KptAx3hsY/ZVlPTAH+HacjRKRGtbYODL1YpyGhsKFq9kS856ecGFssuzbl/35kSMjM96MNz98fT0sXQqnnQZ/8zdwzjnxHlcOitlSaazeSMxOkJid/8XZQ71D4ybnA3sHQk/63kEG9gwwdGCIwb3h+ZH1RwprY8KyJ+adaes7G46Wxs5GXWQq4yp0qMmtwPcYCeTDwHZ3L3CyYxGREquvD1Mvdoy9yKxgydlr9u0Lyff27bB+PTzwAKxdG0p/1Ds3NATr1oVy441hCE1L9q/tY6aYLVWvvrWe+kX1NC9qzvs9w4PDYdz53jCsJT0xT30c3DN6nff50XnRC1XXWjc6QU9Jysck752je+TrGpW014JCE+/vuPudqSvMbJmZHe/uj8bYLhGR4nIPCXSm3u9My+nTRKaXoaHx9/n610Nz/slDDBSzpSbVNdSRmJUgMavwaU+HDg9lTdbHLO8dvX64d5j+3gKHxSTb3FY3Kik/2uveOTZhH1U3Q0l7NSk08X4hcGfauu3A7cDyWFokIrVteDiMx04vvb2Te55piEk+yXK+WlvDTDTz5sHcuSNl8WJYtiwMNZmeeeaIIlLMFilQfUs99fPraZrfVND73J2hg0Njk/I9oxP0jHX7Bhk+NEzfoTBTTMFtbq/POQQmW13DjAbqGpS0l1JeibeZvRP4GNBhZldkeMmtsbZKRMpraCiMYx6v9PXl97psJVOS3Ff4H50Ja2oaO7NJtllPkuPOM5Xp0yvqhkKK2SKlZ2Y0dDTQ0NEAiwp7r7sz1DM0kpTvGdujnrW3fd8gQweHGDo4RN/mCSTt0+rH9qjPaqSxq5FEV4LGrsajJflcc7FPXF6Jt7t/wczuBN4A/Etadb+7b4urQWb2PMIfjFe7+8ZoXQvwKaAH6ASucPcD49WJVCX3cNObOBPcQrc1OFjec9DSEkpr68hytnXjPU+W9ES6ra1ss40Um2K2SHUxMxqmNdAwrYHmxYUNR/NhZ/DAYOYe9bRkPb23fXD/IEMHhhg6MDTuTZJS1U+rz5mYJ+YkSMxLkOhO0HhMo2aJSZH3UBN3X2tmH3D3gfQ6M1vk7o9PtjFmNgdoB85Kq/oScIu7/8DMXghcD/xFHnUiE5fs9T18eHQPba7lQl6ba3l4uLzHbhaS1ebmwkpTU2GvzZQ0686VsVDMFqkNVmc0zmikcUYjLCnsvT7sDO5PS8r3DB6dLaZ/ZzTt486U57sGjibrRx4bf8YYa7RwB9PuBE3zmo4m5KOWu5to7KqNBN3cPXul2d8CN7v7A9HzqwjTUY16GXC+u58XS4PM6gi3M17i7hvNrBvYAHS6e6+ZNQD7gZOB/mx1yZ6XpBUrVvjq1avjaKKUWvJmLcmp4Hp7R0qcyW/6unL2+tbVTSzxnWwynCwNDUp+K4iZrXH3FXm8TjFbRIrKh53BfYPZE/OdYZ72/q399G8LdzXNhzUZzYubQzl25LFpcVN4nNeE1VfH36VcMXu8Hu8lQOpkuk8D2oAnGZmeyoC5k21kkrsPp/3Hcx6wy917o/rBaG7ac4GBHHUb42qTjGNgYPRdC1MT42SZ7Ppy9ACbhSQ0NQHOtDxefSHbaGoKyw2FXvcsAihmi0iRWV2Y57xxZiOty1rHff3Q4SH6t4UkvG9b39GEfNTy1j4G9wxy+JHDHH7kcOb9JoyW41poPbGV1mWttCxroXVZK60nttLYWT3DBnP+dXf3d6St+jjwJ3cfNU+OmT017oalmA/sTlvXA3QDgznqJBv3MKtDchq05M1DkiXX80x1pbgYLpEIwxBSS3JoQhwJb6blxkb1+kpVUcwWkUpT31JPy9IWWpbmvn/B4MFB+jb1cWTjEY5sisrGkceB7QP0PthL74O9Y96b6E7Q8bQO2pe30/60djqe1kHTwqaKHLpSaLfaqe5+b+oKM1sGHAc8EFurRnMg/d+fJkLPCePUHbVz505WrBjp9V+5ciUrV66MsZll0NcHO3eOlOQcw3v3jl7O9BjnNGrJG5MkL1pLT5BTS676bHXJJFikhqxatYpVq1Yln86a4GYUs0WkKjS0N9BwcgNtJ7dlrB86NETvI730PtzL4YcP0/tQWO59pJf+rf3s3rqb3TeN/F/f2NXI9HOmM+PcGcw4bwZtJ7dhdcVLxPON2flOJ/gcwteTz4u+Fkw1k3DV/E8m0M58bCFc+Z6qA9hKCPDZ6kbp6uqiKsYL9vXB1q2weTM88QTs2BGS6kyPByYxEUBr68g0aMm7+SVLe3th65qb1TMsErPURNPMdhXyXsVsEZlq6tvq6VjeQcfy0Xcf9mHn8PrDHLz3ID339hx9HNg5wK4bdrHrhhA+G2c1MvPFM+m6uIvO53dS35J++cvk5Buz8+3x3gd8izB+8Flpdf3AVyfQxnzdDnzFzFrc/bCZNRJmyLyDcEFPtrrKdOgQPPpouH30+vUjCXbycfv2/LfV0ABdXTB7NsyaBTNnhmS6s3P0Y/q66dPDWGIRmaoUs0WkJlid0Xp8K63HtzL7VbOBMC/64UcPs++Ofey/Yz97b99L/5Z+tn9zO9u/uZ261jpmXTyL7rd0M/3Z00s6JCXfebz/aGbnAZe6ezEDNjZy9Bbt+0kzuwl4LnATcCFwg7tviV6fta6sdu2CP/wBfv97eOihkGivWxd6s3Opq4Publi4EObPD3e+mz17JMFOXZ4xQz3NIjKGYraI1DIzo/WEVlpPaKX7sm7cnd6Hetn137vY9aNd9KzuYcd3d7DjuztofWoriz6wiDmvnVOSWVNyTieY90bMnu3ud8WwnQ7gLwnzvH4E+IK77zSzTsINFzYSrsa/0t0PRu/JWpeqqFNT9fTAb34Dd90F994bku0tWf6OJBKwdCmccAIcdxwsWgQLFoREe8GCkGhrRgsRSZHvdIIFbK+2Y7aI1LTDGw6z7avb2Pa1bQxsD5eYtD61leP/+XhmPn/mpLefK2aPN4/3jwk3R8i5fcIcrLMn3sTiizWIDw3B3XfDjTfCL38ZEu30ixXb2uC00+CMM+Dkk0OifcIJIdGuj3dckYhMbQXM462YLSKSp+GBYbZ/dzsbP7zx6J07u9/azfGfO566proJb3cy83g/TBgruIWROWDHbB94yYRbV03uuw/+9V/hv/4rDCVJamiAs86C5zwnPC5fHnqz6yb+QxMRmQDFbBGRPNU11jHvjfOY89o5bP70ZjZ+eCNbv7yV3nW9nHrjqdS3xd9ROl7i/Umgx91zTtRsZsW6Or4yrFkD73sf3HnnyLrjjoOLL4YXvhDOPjv0cIuIlJditohIgeoSdSz+4GJmvmAma1+8ln237ePhtzzMSd8+KfYLL3N2ybr7rvECeOTimNpTWdzhk5+Epz89JN3TpsF73gNr14YLJa+9Fp77XCXdIlIRaj5mi4hMQseZHZzxizOoa6tjx3d3sPfne2Pfx6TGQpjZIjP7JHBdTO2pLF/+MnzgA2E8A/UAAAAV5UlEQVT5ve8NU/798z/DqadqNhERqTpTPmaLiExS21PbWHzVYgC2fD7+CZcKTrzNrN7MLjazm4HHgNcRxhROLYcOwZVXhuV/+zf47GdDj7eISBWpmZgtIhKTuW+YC8De2/biQ5Of/S9V3om3mS0xs48Dm4F/B44HXgosBC6NtVWV4PbbYf/+MMzkDW8od2tERApSczFbRCQmTd1NJOYmGO4dpm9rPqP38jdu4m1ml5rZ/wLrgFcCnwUWAD9191s8uDvWVlWCBx8Mj898ZnnbISJSgJqN2SIiMUrMTwDQv60/1u3mc6eWw4RbDN8NXObujwCYWbx975VmX/RN7KxZ5W2HiEhhajNmi4jEqK459E0P9w3Hut1xE293vwm4ycwWAZeb2bHAjUBj8jVmZh7HLTArSW9veNSMJSJSRWo2ZouIxOho4n0k3sQ77zHe7v64u/898KZo1QIz+2gU1C+PtVWVIHknSt2+XUSqUM3FbBGRGFlDmL3OB8t0cWWSuw+6+w/d/WLgm8B7gc/F2qpKoM4gEZkCaiZmi4jEKO4b5yRNqjvX3R8F3m1ma2JqT+VIJt6ar1tEpogpHbNFRIoh5n7YSd1AJ8ndvxXHdiqKEm8RmaKmZMwWEYlTkdK/WBLvKU2Jt4iIiEhtqsQe7ylJY7xFREREJEZKvLPRUBMRERGR2qShJmWixFtERESkJsV9ywMl3iIiIiIiqdTjLSIiIiJSQrq4UkRERESkeCryBjqVwsxagE8BPUAncIW7Hyhvq0REJBPFbBGpVVOlx/tLwF3u/kHgRuD6MrdHRESyU8wWkeqgoSajmVk38FrgJ9GqnwMvM7Njy9UmERHJTDFbRKqCLq7M6jxgl7v3Arj7ILABOLecjRIRkYzOQzFbRKqFerzHmA/sTlvXA3SXoS0iIpKbYraIVL4i9XhPhYsrHTictq4JGEhdsXPnTlasWHH0+cqVK1m5cmXxWyciMgGrVq1i1apVyaezytmWmClmi8iUk2/MngqJ9xbCVfGpOoCtqSu6urpYvXp1yRolIjIZqYmmme0qc3PipJgtIlNOvjF7Kgw1uR3ojqanwswagUXAHWVtlYiIZKKYLSI1q+oTb3d/ErgJeG606kLgBnffUr5WiYhIJorZIlLLpsJQE4C3AZ8ys9OBuYAGAoqIVC7FbBGpCu7xTmsyJRJvd98LXF7udoiIyPgUs0Wk4mkebxERERGR6qXEW0RERESkBJR4i4iIiIiUgBJvEREREZFMdMt4EREREZHiMSvO1ZVKvEVERERESkCJt4iIiIhICSjxFhEREREpASXeIiIiIiKZ6OJKEREREZEi0p0rRURERESqlxJvEREREZESUOItIiIiIlICSrxFRERERDLRxZUiIiIiIkWkiytFRERERKqXEm8RERERkRJQ4i0iIiIiUgJKvEVERERESkCJt4iIiIhIBu7xTmvSEOvWJsjMDLgU+Ed3PymtbhbwcWAf0Ahc4e4D49WJiEjxKG6LyJQ2xWc1WQx0AidmqPtP4OvufgXwCPCRPOtERKR4FLdFRApUEYm3u28Ebk1fb2ZnA8vc/bfRqp8A7zKztlx1pWiziEgtU9wWESlcRSTekeEM6y4ANiSfuPuWaHHFOHUiIlJ8itsiIgWopMQ7k/nA7rR1PUD3OHUiIlIeitsiMnXEfMv4klxcaWafAE7NUv1zd/9cljoHDqetawIGxqkbY+fOnaxYMdKpsnLlSlauXDlOy0VEymPVqlWsWrUq+XRWqfdf7ritmC0i5RSuH89fvjG7JIm3u185wbduAZ6SfBJdRd8GbB2nboyuri5Wr149wWaIiJRWaqJpZrtKvf9yx23FbBGpJvnG7EofanIrsCzl+bFAP7B6nDoRESkPxW0RkSwqKfE2ONoDAoC7rwa2mdkp0aoXAp939/5cdaVstIhIDVPcFhEpQKXcQGc28Ibo6dvM7FvufjB6/mrgQ2a2EWgHrkp5a646EREpEsVtEakJ1Xhx5XjcfQfw4aik120C3pzlfVnrRESkeBS3RWRKm+J3rhQRERERmdKUeIuIiIiIlIASbxERERGRElDinc306TBvHrS2lrslIiIiIlJCDZ0NJOYlqGuON1WuiIsrK9I114QiIiIiIjVl2fXL4Pr4t6sebxERERGRElDiLSIiIiJSAkq8RURERERKQIm3iIiIiEgJKPHOYdWqVeVuQkXSeRlL52QsnZPMdF6Kq5bOr451atKxTm1KvHOoxQ9EPnRextI5GUvnJDOdl+KqpfOrY52adKxTmxJvEREREZESMHcvdxtKwsx2ApsKfNssYFcRmlPtdF7G0jkZS+cks4mcl8Xu3lWMxlSqCcZsqK3PnY51atKxVr+sMbtmEm8RERERkXLSUBMRERERkRJQ4i0iIiIiUgJKvEVERERESqCh3A2oRGbWAnwK6AE6gSvc/UB5WxUPM3sF8BlgGvAd4P3uPmhms4CPA/uARsIxD0Tvib2uUplZArgHeLe7/zLXZ6EYdZXIzM4Gng08CtwFHKKGz4mZnQS8k3A+TgBWufvv9Vkpr6l2ruKO1dUgrvhbDeKKq5Us7lg5Zbi7SloB/g14dbT8QuD75W5TTMe1CPgWcCbwV4Rf9L+N6m4HnhEtvxX4RMr7Yq+r1AJcBewHzhvvs1CMukorwJuAj6Wtq/VzshqYHy0vAh7UeSl/mUrnqhixuhpKXPG30kuccbWSS9yxcqqUsjeg0grQDfQBrdHzhijoHVvutsVwbM8BGlKefwr4KXA2sDVl/XzgINBWjLpyn4cc5+fPgL8GNgLn5fosFKOu3Mef4XycC/ycaPajaF1Nn5OofYeAE6PlLmCrzkvZfyZT6lwRc6wu9/HkecyxxN9yH0cexxlbXC33seRxrLHFynIfS5xFY7zHOg/Y5e69AO4+CGwg/LJUNXf/VXQ8SVuAx4ELCMeYfN2WaHFFkeoqjpm1AZe4+9dTVp9H9s9CMeoqzWeBB4HPm9ktZvYsdE4Afgh81cw6gL8E3oXOS7mdxxQ6V0WI1RUt5vhb6eKMq5Uuzlg5ZSjxHms+sDttXQ/hP7Gp5unAl8l9zMWoq0RXANekravZ82Jmy4DlwFfc/Z3AbcAtwEJq9JykeDtwhDAW9aC730ANf1YqxFQ/V5ON1ZUuzvhbsYoQVytdnLFyylDiPZYDh9PWNQFVc4FKPszseGCnu/+B3MdcjLqKYmYvAFa7+460qlo+LycDe9x9bfT8C4BR2+ckqQX4PvA94HNmdiE6L+U2Zc9VTLG6YhUh/layuONqpYszVk4ZmtVkrC2EK2lTdRDGJk0JZlYPvAX4QLRqC/CUlHojjBncWqS6SvO3wPLQRCD8/G8kfCWY7bPgRairJA2kxAd3P2xm64B6avecJH0HeI277zMzJ/xheQ86L+U0JeN2jLG6ksUdfytZ3HG10sUZK6eOcg8yr7QCzCVckNISPW8kDPafX+62xXiMVwBzUp6vADalPF9C+HonUYy6ch9/hvPRBSxIKZuBSwlXYWf8LOT6nEy0rtznIe2cnEgIgqmfk3uAN9TqOYnaNgvYlvLcgMcI04LV7Hkpd5mq54qYYnW5j2OcY4w1/pb7eMY51ljjarmPZ5xjjTVWlvt44iwaapLG3Z8EbgKeG626ELjBRy5UqWpm9kHgPqDFzJaa2V8T5nzdZmanRC97IfB5d+9399Vx15XgMAvi7jvd/YlkAYYIX+0+TpbPQq7PyUTrin2chXD3h4CfAZcAmNkMwh/371Kj5ySyBzhiZvNT1u0G1lLb56WspuK5ijNWl7zxBYg7/pa6/YWIO66Wsu0TEGusLFGbS8Ki/yokhZl1EqZv2kj4D+xKdz9Y1kbFwMyuBj6atvohdz/JzBYDHyIccztwlUdX1RejrpKZ2UbgjR5u4JD1s1CMukoS3YzjOsJcrIsIFwQ9UMvnBMDMTidcNLQGmAPc6e531Pp5KbepdK6KEaurRRzxt5LFHVcrWdyxcqpQ4i0iIiIiUgIaaiIiIiIiUgJKvEVERERESkCJt4iIiIhICSjxFhEREREpASXeIiIiIiIloMRbRERERKQElHiLiIiIiJSAEm8RqRpm1mVmPzazTWb2gXK3R0RExqfYPUI30BGRqmFmLwBuB6YBDwFPcfc95W2ViIjkotg9Qom3iFQlM7sTuKDaboktIlLLaj12a6iJiFQdM1sI/G+tBm4RkWqk2K3EW0rMzP6fmX2x3O1IiqM9Znahmf3GzM6LqVlTipl9OubtGfAq4Jo4tysiYylmT56ZHWtmXzGz62PaXr2ZXW5mm/J47YvN7LcTPdY447did6DEewoxs5PM7Edm9lUz+4OZuZn9ZbnbleZm4PvlbkSKONqzGTg9hrZUJDO7yMx2mdk2M3ummS2MPmduZl8wsxnR684ysz+Z2WfMrCVaNwd4MuYmvRb4mrsPRtsXqUqK2RNSjTF7R/TYFNP26oB9wKI8XrsaOClbpZmdbGY/iD57V5vZ3JS6uOO3YjdKvKea/wb+0d0vA5YD30pWmNlpZnZuORqVum93/5O731WOdmRSSHuynUN3fxjYFWe74vx5TXZb7n4z8CVg0N1/4+6bgcuBfuC37r4vet3vCEH+79z9cPT2lwI3TeoAUpjZ3wMfBn5uZg8By+LatkgZKGYXqBpjtrv3Alvj2pe7DwBr8nztdkKSnq3+fmBV9PTT7p6aaMcWvxW7RyjxniLMbCbwFOAQgLsPAx8C3MymEwK6laFdZdt3nPI4juES7qsc2/oBsMDMzgRw993AbcBLUvY1A9gaffaSTnT3hya576Pc/aPufoK7r3D3E939zri2LVJKitnFVYExO+6ZLApp/3ivHYoe08ddxxa/FbtHKPGeIqJpeR4GbjGzs6N1m4DfEcZULQFWmtmbAMzsI2b299FXTEfHcJnZCjP7rJm9x8x6zexhM3ufmf25md1jZheb2X+Z2QEze0vK+zJuL23fl5vZN1LHuUVftX7ezD5hZreb2VnR+pz7S3m/mdlbzOx+M3uema01s+1m9sY89nFCanvG2eeYc5hBZ/S1cY+ZXZ52Tj8cnZsfRUE633M2Zl9mdomZ/c7M3mhm68zsmhzby/Szz9ieXKJekfuBS1NWdwEvtmhYCfAy4Mcp7WwBelOen2hm15jZp83sDgtDVl5tZqvN7OVm9n0z22dm7zazM6Nj3G1mF43XPpFqo5hdOzE7k2zHmGs/ZrbEzL5uoff42nz2k2G/N5vZd83sWeO8TvG7WNxdZYoU4GRgE+G/268Ds1PqNgLnRcunAOuj5S7Cf+LHEP5bfxxYEdV9Ebg7Wk4Ae4AvEObh/AtgS67tpe+bML7teuAb0fpmQjI3I3p+CbAT6Mi1vwzHPTfa5wej134MOAIsHmcf6e3Juc/Uc5ihDRuBzwCt0fsej9ZPA/495XW/JgTMvM5Zln3NJHx1+DngLODCXNtL+9lnbE+en6+rgcei5eOATwB9wCtSPi+W8vqXAc+JlluBe4DW6Pmd0XmYBuwHPgK0AK+Onr+K0DHwfqLPoIrKVCsoZtdEzI7qP5znecz1s74HOCnlPZ7n5yz583wm8A8Z6s+L9tOQsk7xu0hFPd5TiIdeyVMJQeivgPvNLNMFJOuA15tZI/CcaF074Zd8IZAcn/sg4ZcJd+8HeoAfuvsBwi/h3HG2l96+PkZfqPFi4LBHY4SBHxF+WV8xzv7St5vc5vei136MMP74eePsY1R7CtlnFj/xMJbvHmBeyjEeE/VGvYdwTiHPc5blePcQgtuN7v47d/95AdvL1p58/ABYamG4ySuBrwD/A7zKzNqBHo+icuRZwN3R8ssIf9iSPSgXAx+NzvM+4BcexoWvBqa5+394+Or9Pgr7GYhUDcXs2ojZGWQ9xhz7uRCY7u7J9qwucJ9PB17n7v+Y5+sVv4ukodwNkHiYWT2w0N03An9jZl8jXP29CnhG6mvdvc/M5gNXAN9MbsLdd5jZauC5hP/GlxCSraNvTVkeJBqqlG17WZqauo0TgMaUdg2b2UZgQa79jcfdj5jZemAGMKuAfUx4n2kGGfndWgRscPfPpb8on3NmZq8nJLhJyzx8He2pbS3gZ5C1PeNx93Vmdh9huMl0d19vZj8g9EBdCvw0pd3Jz0Zy7OBiQu9Uclu7UzedsjzEaMNoSJxMQYrZo46vFmJ2qqznMcfP5jTCPxkTdQFwgpm1u/vBLK8Zjo5B8buIdEKmjkZSxt+6+x8IX+OdmP5CM7sAeLu7f8zDDBWp3gE83cKUVmuBT46343G2l8vjhB7U1CmWjDDucbISwCNF3kc+tjF6HDRm9owCztmPgTNSSsYr4wvYXsb25H84/DvwJmB9SvvqgXcTvpJNegZhrGrSFuDZZtaWst+cYwxFpjjF7NFqKWZnPcYc+zkQvSfBxHya8G3BdVnq9/vIhfGK30WkxHtqeZeZdac8nw/cEi33Ey4kOZEwbdV0M2syswuj+tlmNoswzvBrhK+I7ib62jJSR9p/+GZm42wvfd+Wso3/Bg4Cr4m2NSPaX7LnNNv+spkXvWZ+9N5b8thHanvG22fqcaTL1PNhhB6sNuAmM3uBhQt4xvsZHN0XMM/dH0opAym7qE9ZzrW91HZnbI+ZJcxsVTSMJJf/IHy9/UMAd+8Bfgb82kfPZvJ8wjCUpOT5/r6FucDfCSTncU3/GaT/nKt6dgWRHBSzqZmYne95zLgfQjxNEP45g/DtBmbWnUf8NkLv8+uBS83sVdF7329mVxN6+29Leb3idzEVY+C4SukL4WINB/YC3wG+QfiaqjOq/wfgCeDlhK+N1hF6F15M6CX5BiGA/ZhwkUd/tL1ewi/hi4ABwsU7C4Ero/q/zrW9tH2/G/gVYczcyVHdCkJP6WcI4xxPidZn3V+W4/fo/e8mDH04NaUu2z6OTW3PePtMPYdp+3454SLDL0TnIvm+N0X15wJ/in42n4nOc77n7OUZjvWSaH/fB+ZH63L9TEdtK0t7FhBu8vC9PD5r/5L2/NXA+WnrPp3hfecADwC7gaujdRcRPmufB7qJplMD3kr4o/yV6FhfUu7fMRWVOAuK2bUUsxcBdxF67c8Y5xhz/axfTPi28f8IF7s/Arw9OvaM8ZvwrcoAoad7JvALQtJ/GfAWwnzmnyEMc0m+R/G7iMWikyeCmc0G3uruH4me1xF6YN7i7leXtXHjMDMHlngYLykTZGZ/5+4TmqZKREpLMVtSKX5XBw01kVTvJAw7mJmy7gJSLpqrRClfa+krrQkys3ozex1wQ7nbIiJ5U8wWxe8qo8RbUn2ZMA7sQTPbDNwKPOHuvylvs7KLLjR5d/T0jRamtZPCtQP/7e7rx32liFQKxWwBxe+qoqEmIiIiIiIloB5vEREREZESUOItIiIiIlICSrxFREREREpAibeIiIiISAko8RYRERERKQEl3iIiIiIiJaDEW0RERESkBJR4i4iIiIiUgBJvEREREZES+P9jqEN2z7Y5wwAAAABJRU5ErkJggg==\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/galileo-jupiter.png',bbox_inches='tight')\n", "plt.savefig('../plots/galileo-jupiter.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/galileo-jupiter.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Results are off because presumably because of poor validity of heating correlations under extreme entry conditions at Jupiter." ] } ], "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 }