{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 40 - Oceanus - Uranus (Concept)" ] }, { "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 Uranus entry probe which was part of the Oceanus mission study\n", "https://doi.org/10.1016/j.asr.2017.02.012" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"URANUS\") \n", "planet.h_skip = 1000.0E3\n", "planet.h_trap = -100.0E3\n", "planet.loadAtmosphereModel('../atmdata/Uranus/uranus-ames.dat', 0 , 1 , 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Oceanus-Uranus-Probe', 120.0, 170.0, 0.0, np.pi*0.8**2.0*0.25, 0.0, 0.28, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(1000.0,0.0,0.0,21.47,0.0,-31.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 (100*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": "iVBORw0KGgoAAAANSUhEUgAAAt0AAAHtCAYAAADFpcLCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXwkdZ3/8ddnck1mJnOHGTIwDDfjoHhEUFAZBRRFEVkQQRT8KXFxxdX9reIiui54gReuLi5B8UZYxJVjBX+wAi6CSkauYTllLmeYIZkzcyaZfH5/VPWkk0l3utNVXdXd7+fjUY/urm931bc6yac/+fanvmXujoiIiIiIxGdC0h0QEREREal2SrpFRERERGKmpFtEREREJGZKukVEREREYqakW0REREQkZkq6RURERERiVp90B8ph9uzZvmDBgoKf393dTWtra3wdSrlaP37Qe6DjT8/xL1mypMfd09GZMik2Zmek6edWLfSeRk/vafTS9J7mi9k1kXQvWLCArq6ugp/f3t5e1POrTa0fP+g90PGn5/jNbEXSfSi3YmN2Rpp+btVC72n09J5GL03vab6YrfISEREREZGYKekWEREREYmZku5RdHR0JN2FRNX68YPeAx1/bR9/pdLPLXp6T6On9zR6lfKemruXd4dmJwFfAM5y9+XhumbgSqAXmAFc7O5bSmnL1t7e7mmp9RERKYaZLXH39gT3r5gtIlKgfDG7rCPdZjYHmAIcPaLpauB+d78EuAW4JoI2EREpgWK2iEh0khjpngDsBg509+Vm1gYsA2a4+3Yzqwc2A4uAvvG0ZUZjMjRqIiJJ2L17N4ODgzQ0NIx7GykY6U59zHaHwcHgfl1dKUcrIlKa1Ix0A7j74IhVi4Eed98etg8QBObjS2gTEUnUF7/4Rerr6/n85z+fdFdKUgkx+4EHoL4e3vCGUrckIhKfNMzTPQ9YP2JdL9AGDIyzbZju7m7a24f+6ejo6KiYonsRqUyNjY0A7Nq1q+jXdnZ20tnZmXk4O7peRUIxW0QkS6ExOw1JtwM7RqxrAvrD++Nt26O1tTU1k6aLSG1oamoCxpd0ZyeZZtYTacdKl9qYXeZqSRERoPCYnYakezXBWezZWoA1BMF9PG0iIokqJelOudTFbLNStyAiEr80zNN9D9AWTiWFmTUA84H7SmgTEUlUFSfditkiIuNQ9qTbbM+YhAG4+1rgduCEcP2JwM3uvnq8bWU4DBGRvKol6a6kmK3yEhFJs7KWl5hZC3Bu+PB8M/uOu3cDFwJXmtlRwFwg+4yZ8baJiCQmk3Tv3Lkz4Z6MX6XEbJWXiEglKGvS7e69wHfDJXv9RuCCHK8ZV5uISJImTpwIVPZIt2K2iEh00lDTLSJSdaqlvKSSqLxERNJMSbeISAyUdJePyktEpBIo6RYRiYGSbhERyaakW0QkBkq6y0/lJSKSZkq6RURioKS7fFReIiKVQEm3iEgMqmHKQBERiY6SbhGRGFTDlIGVRuUlIpJmSrpFRGKg8pLyUXmJiFQCJd0iIjFQ0i0iItmUdIuIxCA76XbVPZSF3mYRSTMl3SIiMairq6Ourg53Z2BgIOnuVDWVl4hIJVDSLSISE5WYiIhIhpJuEZGYaNrA8lJ5iYikmZJuEZGYaNrA8lB5iYhUAiXdIiIxUXmJiIhkKOkWEYmJku7yUnmJiKSZkm4RkZgo6S4PlZeISCVQ0i0iEhMl3SIikqGkW0QkJkq6y0vlJSKSZkq6RURioikDy0PlJSJSCZR0i4jERFMGiohIRn3SHcgws4XAR4HngEOBTnd/xMyagSuBXmAGcLG7bwlfk7NNRCRp1VxeksaYrfISEUmz1CTdwE+Ad7r7ajObD/wGWAhcDdzp7jea2cnANcDZ4WvytYmIJKqak25SFLNVXiIilSBN5SULgZbw/g5gmpm1AecAt4Xr7wZONbMF+drK1mMRkTyqPOlOXczWSLeIpFmaku5fAN8zsxbgXOAiYDHQ4+7bAdx9AFgGHD9Gm4hI4qo86VbMFhEpQpqS7o8AO4GHgK3ufjMwD1g/4nm9QNsYbSIiiavypDs1MVvlJSJSCdJU090M/JwgMF9lZssAJ/jaMlsT0B/ez9e2R3d3N+3t7Xsed3R00NHREVG3RURGN94pAzs7O+ns7Mw8nB1tryKTupit8hIRSUKhMTtNSfdPgfe4+yYzc4Jg/nGCM9yztQBrCIJ7rrZhWltb6erqir7HIiJ5jHfKwOwk08x6Iu9YNFITszXSLSJJKjRmp6K8xMxmA0e5+6Zw1ReALcAKoC2cZgozawDmA/cB9+RpExFJXLWWlyhmi4gULxVJN7AB2Glm87LWrQceA24HTgjXnQjc7O6r3X1trrYy9VlEJK9qTbpJacxWeYmIpFkqykvcfdDMTgM+Z2ZLgDnAJ919i5ldCFxpZkcBc4Hswr58bSIiiarWpDttMVvlJSJSCVKRdAO4+6PAh0dZvxG4IMdrcraJiCStWpNuUMwWESlWWspLRESqTjUn3Wmk8hIRSTMl3SIiMRnvlIFSHJWXiEglUNItIhKT8U4ZKCIi1UdJt4hITFReUl4qLxGRNFPSLSISEyXd5aHyEhGpBEq6RURioqRbREQylHSLiMRESXd5qbxERNJMSbeISEyUdJeHyktEpBIo6RYRiYmmDBQRkQwl3SIiMdGUgeWl8hIRSTMl3SIiMVF5SXmovEREKoGSbhGRmCjpFhGRDCXdIiIxUdJdXiovEZE0U9ItIhKT7KTblRHGRuUlIlIJlHSLiMRkwoQJ1NfXA9Df359wb0REJElKukVEYqRpA8tHXyaISJpFknSb2dui2I6ISLVJ47SB1RazVV4iIpWgvpgnm9nJwKXAXIYS9gnAHKA52q6JiFS+JE+mVMwWEUmPopJu4CfAl4GHgcFwnQHnRNkpEZFqkfAMJjUVs1VeIiJpVmzS/by7f2PkSjN7KqL+iIhUlYST7pqI2SovEZFKUGxN94VmdsYo60+NojMiItUm4aRbMVtEJCWKTbqPA75nZruzlkHguzH0TUSk4iWcdNdUzFZ5iYikWbHlJZcBHUAXsDtcZ8D7ouqQmb0GeB3wHHA/sA24EugFZgAXu/uW8LnNudpERNIg4SkDayJmq7xERCpBsUn3HcAt7j5syMbMro2iM2b2AeBgd780a90PgDvd/cbwTPxrgLPD5qvztImIJC7hKQMVs0VEUqLYpPsG4J/M7LdZ6ww4Hfj7UjpiZscD7wVOylrXRnCW/d+Fq+4GbjazBUBfrjZ3X15KX0REopJweUlNxWyVl4hImhWbdH8WOBQ4H8iENyOYA7akAA58A3gA+LaZHQJcDhwA9Lj7dgB3HzCzZcDxQH+etuUl9kVEJBIJJ901EbNVXiIilaDYpPurwB3u3pu9stSrm5nZ4cArgA+4+2Nm9kngTuALwPoRT+8F2oCBPG3DdHd3097evudxR0cHHR0dpXRZRKQg40m6Ozs76ezszDycXcLuaypma6RbRJJQaMwuNuleNTJ4h+YWuZ2RFgEb3P2x8PF3gH8mGJnZMeK5TQQjJozRtkdraytdXV0ldlFEpHjjSbqzk0wz6ylh9zURszXSLSJJKjRmFztl4EdGrjCzM4GritzOSPVk/QPg7juAZ4E6gjPcs7UAa4DVedpERFIh4fISxWwRkZQoNumeaWafBTCzSWb2PeB7wEMl9uMxYJqZzclaN0AQjNvCaaYwswZgPnAfcE+eNhGRVEg46a6pmK3yEhFJs6KSbnc/BbjLzL4NPAIcBhwFjHbFs2K2+xTB1FZnAJjZdKAR+BlwO3BC+NQTgZvdfbW7r83VVkpfRESilJkyMIl5umslZqu8REQqQd6abjN7Q46mFQRfJX6N4Gz104BPlNiX9wPfMrNGgtGPs8Oz2y8ErjSzowjqELPPpsnXJiKSuHKOdCtmi4ik11gnUl4JLAQ2AYOjtF9HMP3UvpQYwN29h2DO15HrNwIX5HhNzjYRkTQoc3lJTcdslZeISJqNlXRfBnS5+4v5nmRmJX1VKSJSrcqcdNdkzFZ5iYhUgrxJt7v/upCNuPsvoumOiEh1KWfSrZgtIpJexc5eIiIiRUh49pKaovISEUkzJd0iIjFS0h0/lZeISCVQ0i0iEqMkpwwUEZH0KDnpNrNXmNlPzez/mtmC0rskIlI90jbSXc0xW+UlIpJmUYx0/y1wort/HTg2c7UxERFJX9JNFcZslZeISCWIIun+KPCq8P4NQEME2xQRqQopTLoVs0VEElB00m1mbzWz94b3jwYOzVzG190H3X1LxH0UEalYSSfdtRSzVV4iImlWVNJtZl8EbgLeBeDufwLeb2bnxNA3EZGKl2TSXSsxW+UlIlIJih3pPhbYD/ifrHV3AF+JrEciIlUk4ZFuxWwRkZQoNum+y903Adlf4lXV5YRFRKKU8JSBNRWzVV4iImmW9zLwo3jezC4B5pnZe4AzCb62vCjynomIVIGER7prImarvEREKkFRI93ufgPwW6AOeD+wCXizu/9bDH0TEal4SSbditkiIulR7Eg37v4H4A+Zx2ZWZ2YnuPt/R9ozEZEqkPTsJbUQsyeEw0cDA8n2Q0Qkn7xJt5ndCkwZYxuzgR6gagK4iEhUypl012rMDsvmSc9U6CIiextrpPux8DlPA4PABcDDQFfWcw4EZsbSOxGRCtfY2AhAX18f7o7FW4BckzG7Obym5o4dyfZDRCSfsZLurwIWnv2OmR3s7p8b+SQzuyWOzomIVLoJEybQ0NBAf38/fX19e0a+Y1KTMTsz0r1jRzCDiU6sFJE0ynsipbtvzgTv0Bwzq8t+jpktBo6OoW8iIlWhXNMG1mrMrquDhvBi9n19yfZFRCSXYk+k/AXwqJndD/QCC4G3AHuNpIiISKCpqYne3t4kTqasmZjd3Az9/cFod7xfJoiIjE+xUwbeBbyZ4CScwwmmnzrD3b8cQ99ERKpCUjOY1FLMzi4xERFJo/FMGbgGuDTz2MwmmNmr3H1JpD0TEakSCc/VXRMxWydTikjaFZV0m9l1o6xuBXYDp5XaGTNrBB4C/t7d7zWzZuBKgq9FZwAXu/uW8Lk520RE0iSppDvumB3uIxVxe9Kk4Hb79ii2JiISvaLKS4BDgZXAiqylCYhqxOSTwIKsx1cD97v7JcAtwDUFtomIpEaCI91xx2xISdyePj243bQp//NERJJSbHnJh9z96ewVZnYEcH6pHTGzY4EXgI3h4zbgHODvwqfcDdxsZguAvlxt7r681L6IiEQp5rm584ktZofbSk3cnjEjuN24sdQtiYjEo9gTKZ8eZfUu4EOldMLMJhOc3JP9VehioMfdt4f7HgCWAceP0SYikiruntR+Y4nZkL64PTO83M+GDVFsTUQkesXWdC8Dsj896oA5wE9K7MfFwFdGrJsHrB+xrhdoAwbytO2lu7ub9vb2PY87Ojro6Ogopb8iIkUrdMS7s7OTzs7OzMPZJewvrpgNMcbt8cTszEi3km4RKbdCY3ax5SV3AdczFMQHgXXu/kzRPQyZ2VuALnd/cUSTAyPPQ28C+sP7+dqGaW1tpaura7QmEZHYFTvSnZ1kmllPCbuOPGaHfYo1bo8nZu+3X3D7l78U9TIRkZIVGrOLTbp/6u6/y15hZoeb2SHu/lzx3QTgH4FXZI0AzSA4weYb4f1sLcAagsCeq01EJFUySXcCtd1xxGxIYdx+5SuD2z//OYqtiYhEr9jZS04eZd064KYS+nAO8PKsZQ1BveEPgLZwiinMrAGYD9wH3JOnTUQklRJIuuOI2ZDCuJ1Juh95RJeCF5F0KijpNrOPmtkm4GIz2529ENTorRtvB9y9293/mlkI5o/tdveVwO3ACeFTTwRudvfV7r42V9t4+yEiEpdyn0gZZ8yGdMbtmTNh0aLg4jgPPhjFFkVEolVQeYm7f8fMfgecB/zriOY+d38h8p4FLgSuNLOjgLlAR4FtIiKpsXXrVgAmT55clv0lGLMhwbh98snwxBNwxx1wvOayEpGUsWJGYMyswd33OunFzOaHIxyp1N7e7jqRUkSSMnv2bNavX8+LL75Ia2trUa81syXu3j72M0d9bU3F7LvvhpNOgqOOCspMRETKLV/MzjvSbWb/CPza3f83XPUpM6sb+TTgjQRzsIqISBZ3Z/PmzQBMnTo11n3Vesx+3euCy8E/+iisWgX77590j0REhoxV030gMD3r8SuBY4GDwrbMMjeW3omIVLgdO3YwMDBAU1PTnsvBx6imY/bEiXDKKcH9n/402b6IiIyUd6Tb3f9uxKovAkvdfdi54Wb2kqg7JiJSDV58MZjKetasWbHvSzEbzj8fbroJfvAD+PSnofwTxoiIjK7YKQNfOkrwPhw4OLouiYhUj1WrVgEwf/78JHZfczH7zW+GtjZ49ln4zW+S7o2IyJCCZi8xs9cT1AGeFF5WONtMgrPjb4u4byIiFW/lyuB8xf3LWGBcyzG7vh4+9rFglPurXw1mNBERSYNCr0i5CfgxQS3gcSPa+oDvRdkpEZFqsWLFCqC8STc1HrM//GH40pfgt7+F//qvoTpvEZEkFTpP9+Nmthg4092rOliLiETp8ccfB2DRokVl22etx+zp0+Hzn4d/+Af4+MfhxBMh/nNYRUTyK7im29035wreZva66LokIlI9HgknjD7qqKPKut9aj9kf/SgsXAjPPQff/GbSvRERGXue7luBKWNsw4BFwD5RdUpEpBps3bqVZ555hrq6urKMdCtmD2logG99Kzix8rLL4B3vCC4TLyKSlLHKS54mqA1cDeS6dKUBb4+yUyIi1eC+++5jcHCQo48+mokTJ5Zjl4rZWU46Cc47D370IzjrLPjjH2Hy5KR7JSK1aqyk+wqg19135XuSmVXlWfAiIqW4++67ATjppJPKtUvF7BG+8x34wx/giSfg7LPhl78MZjgRESm3vDXd7t4zVvAOvSui/oiIVAV357bbgty2XEm3YvbepkyBX/0KZs6E226Dj3wEPNd3ACIiMSr24jjDmNl8M7sC+FZE/RERqQp/+MMf+Mtf/sK+++7L616XjvMWazVmH3EE3HprcJn4a68N5vEeHEy6VyJSa4pOus2szszeZWa/Bv4CvJeghlBEREI//OEPATjnnHOoq6tLrB+K2YHjjoP/+A9obAxKTt73PujrG/t1IiJRKTjpNrMDzeyLwCrgBuAQ4B3A/sCZ8XRPRKTyrFu3jh/96EcAfPCDH0ykD4rZe3vHO+COO4KSk+uvhze9CdatS7pXIlIrxky6zexMM/t/wLPA3wDfAPYD/svd7/TAAzH3U0SkYnzzm99k165dnHrqqSxcuLCs+1bMzu9Nb4L77oN58+D3v4dXvQruvTfpXolILShkpHsHwWWDHwBOdfevuXs3uaejEhGpWc8//zxXXXUVAJdcckkSXVDMHsMrXwldXXDssbB6dZCIf/rTKjcRkXiNmXS7++3u/nbgXOB9ZvYTMzsDaMg8x8wsxj6KiFQEd+cTn/gEu3bt4r3vfS/HHHNMEn1QzC7A3LnBCPdnPwtmcMUV0N4ejIKLiMShmMvAr3T3zwIfCFftZ2aXm9kC4IIY+iYiUlGuu+46br31VlpaWrjiiisS7Yti9tgaGoKrVf7ud3DggfD447B4MbznPbByZdK9E5FqU/TsJe4+4O6/cPd3AT8CPgFcFXnPREQqyOOPP87HPvYxAK6++mrmzZuXcI8CitljO+644OI5l10Gzc1w441w2GHw8Y/D2rVJ905EqkVJ83S7+3Pu/vfA30bUHxGRirNmzRpOOeUUtm/fzrnnnsu5556bdJdGpZidW3NzUGry5JPBSPeuXfCtb8FBB8FHPwpPPZV0D0Wk0pWUdGe4+49L3YaZnW5my8xsvZl9y8zqw/WzzewaM7vCzL5hZtl1iTnbRETKYd26dbzlLW9h1apVHHvssVx77bVJd2lMitm5HXAA/Pzn8Mgj8M53wo4d8G//BgsXwpvfDDfcANu3J91LEalEkSTdpTKz+cBpwBkEX31+CPh42HwTcJ27Xww8A1yW9dJ8bSIisVqxYgVveMMbWLp0KQsXLuSWW25h4sSJSXcrdrUQs486Krh8/KOPwgUXBCPhd90FZ58Nc+bA+98Pt9+uBFxECmfuyc8iZWavBx5094Hw8ZXAIuBy4Jfu3haunwc8DcwBXpqrzd23ZW+/vb3du7q6ynU4IlIDfvvb33LWWWfR09PDUUcdxV133UVra2vk+zGzJe7eHvmGS1CLMXvDBvjZz4Llj38cWj9xIrzxjXDiiXDMMcF0hM3NyfVTRJKVL2bXl7szo3H3/xmxajXQArwJWJb1vNXhTFftwHF52jTpk4jEYufOnVx++eV85StfYXBwkJNPPpnrr7+eGTNmJN21sqnFmD1zJlx0UbA891xQgnLrrcF833fcESwA9fVw5JGwaFFQkvKSlwTLwQcHbSJSu9IaAl4NfBXoANaPaOsF2oB5edpERCLl7txyyy18+tOf5umnn8bMuOSSS7jsssuoq6tLuntJq6mYfcghwUmXn/1scBn5O++E+++HP/0Jli4N6sEfeWT4axoagnrxgw4KlgMPHH6/hv5nE6lZqUu6zewQoNvdHzUzJ7i6WrYmoJ/g6mq52obp7u6mvX1opL+jo4OOjo5I+y0i1am/v59bb72VL3/5yyxZsgSAI444gu9///sce+yxseyzs7OTzs7OzMPZsewkIrUes+fMgfPOCxaArVuD+b6ffBL+93+HlhUrghHy554bfTvTpgWXpt93X2hrG7rNvj9nDkyaFFzMR0TSo9CYnYqa7gwzqwO+Alzq7rvM7J+AN7r7m8N2A3YCbwSOz9Xm7g9kbzeN9YEikl7uzhNPPMGNN97Iddddx5o1awCYO3cun/nMZ7jgggtoamoqS1/SWNOdoZhduO3bYflyeP75YFm2bPj9bdvG3AQATU0we/bwZdasvR/PmBEk8tOmwdSpQe25knWR+KW+pjvLPwJfc/dd4eO7GD6f7AKgD+gKb3O1iYgUZceOHTzwwAPcfffd/Od//idPP/30nrYjjjiCj3zkI3zwgx9k0qRJCfYydRSzCzRp0lB990ju0NMDL7wQLGvWDN1m33/xRdi5E1avDpZiNDQEyXcmCR/rNtdSpv81RapSapJuM7uEIPg2m9lBwGLgd8ALZnakuy8FTga+7e59QJeZ5WoTEclr7dq1PPTQQzz00EP87ne/48EHH6Svbyh8zJo1i3e9612cc845LF68GNMw4TCK2dExg9bWYHnZy/I/d/v2IEHPLOvX7/24uxs2bYItW2Dz5uB2166gbf3IqvoiZZL3QpeWltHXT54ME1IxabFI+aQi6TazSwmmmsr2lLtfZ2ZnAZ8zs+XAFOAzWc/J1yYigruzcuVKli5dyqOPPron0V49YqjQzHjFK17BG9/4Rt761reyePFi6jXdxKgUs5MzaRLMnx8sxdi1a3gSPtZtZuntHbq/eTP090eTvJvtnZBPmxbMEpO9zJq197oZMzQTjFSmVNV0x6Ua6wNFZDh354UXXuDJJ59k6dKlw5atW7fu9fyWlhZe9apX0d7ezmtf+1qOP/54Zs2alUDP80tzTXdcFLPTK5O851uyE/VcS6E17LlMnbp3Mj7akknaZ88ObpWsS9wqqaZbRCQnd6e7u5tnn3122PLMM8/w3HPPsS3HJ/k+++zDkUceyZFHHkl7ezvt7e0cfvjhTND32yJFaWoaKoUpxe7dQ8l5b28wir5pU3ARorGWjRuHkvflywvfp1kwSt7aGiTh2bej3Z8+PbjQUUODTkKVaCjpFpHUcHfWr1/PihUrWLly5Z5lxYoVrFixgmeffZbNmzfnfP2sWbM4/PDD9yTYRx55JIsWLWKfffYp41GIyFjq6oKkdvr04l87OBgk3IUk6Bs2DNW9r18/tC7rPOmC+jpp0vCluXnvdbnWZ9ZNnjy0TJky/HFjoxL7WqCkW0TKYvv27bzwwgu88MILrFmzZs/9zONVq1axcuVKduwYOZXzcFOnTuXQQw/lsMMO49BDDx22zJw5s0xHIyJJmTBhKGE/6KDCX7d7d5Bw9/QEJ5t2d+e+390djMBv2zY0Kt/bG98x1dXlTsjzJeuFtCmhTw8l3SIyLu7Opk2b6O7uHnV58cUXhyXVW7ZsKWi706dPZ/78+cOWAw44gPnz53PIIYfQ2tqqmUREpGh1dUPlIwsXFv66/v5g1pgdO4Lb0Zax2rZt23vZunXo/sDAUMlMHMc93kQ+s2RG6rNvM4uq9AqnpFtE2LlzJxs3bmTjxo1s2LBh2O3GjRtZv3493d3d9PT07Emqe3p6GBgYKHgfjY2N7LvvvnuWtra2YY/3339/5s+fz9SpU2M8UhGR4jQ0DF1oKC59fbkT8nzJeiFtAwNBzXyeyrySTJw4ekKeL1kvpq2aRuqVdItUOHdn586dbNmyhS1btrB58+Y99zOPN23alDepHqukI5eWlhZaW1tHXfbZZ59hSfXMmTM1Qi0iMorGxmCZMSP6bWcS+vEm8tu3D91m39+2LRjF37kzWEqdRjKX7Jr6qBL57Nvm5mAf5aCkWyQB7s6uXbvYunXrsKW3t3fPba4EOvtxZl0xI86jaWhoYMaMGcyYMYOZM2cOu80soyXW5boUuoiIjE+cCf3gYJBwj0zOcyXpY7WNtq6/P/6a+paWodKj1lbYd19405vgzDOjTciVdIuMYawEeax1udaXmihna2pqYurUqUybNo2pU6cOW6ZNm8a0adP2SqSz70+ePFmj0CIiUpQJE4ZGk+OSqakvNEnv7R1+ldbsZffu0feRSeqff35o3bXXBvPSn3dedMeipFuqRn9/P9u2bSt42bp165jtcSTIGY2NjUyZMoWWlhamTJkybGlpadkrgc6VULe0tGjEWUREYuc+VE6ya1ew9PXtfX+0daPd7+sLkurM/WLW5XpOVB/XRx4Jxx4bzbYylHRL2bg7fX19RSXGhSTHmaW/vz+2vo+VII9cN9ZzJ0+eTGNjY2z9FRGR2uAejOCOTEDzLbt2DY0Oj1bbvXXr8CuIbt48dD+GMahImQ3NjZ49b3r2uubmYGJTIeEAACAASURBVMrJGTOCZebM4HbWLDjgANhvv6AkJ2pKumWYTClFKQlwvufszvXdTgTq6uqYPHlyUUsmAR5tUYIsIlLZBgeDJDEzAlrI/Up6bn9/kEC7l+89bWwMktampmBpbBx+W+j9hoahdY2NwePM/VLW1dend7YTJd0VKDsxLjUJHu05g4ODsfW9vr6+6OS30MS5sbFRdckiIkVw3zvZG5n4Vdq67MQ0xo+zVDEbnsBm3x9taWgIRn1Hm6M7s27aNJg6deh26tTghENVM46fku4y27VrF5s2bdozC8VYt6Ot6+3tjTUxbmhoGFdSXMhzNGIsItUok7xm16yWsmTqU+NOUGP88jE1GhqC0c/6+qH7o63LdT/Nz80k0eWa8k5Ko6S7RO5OT08PK1euZNWqVbz44ov09PTkXHojmvMmkxhHkQyPbG9oaIikjyIiaTQwABs3BpcEX78+WDZvzj1Pca5lx47hiXI5v+KPUl1dkMiNTPDSvK6Q19TXB8emL0AlLZR0F8jdWbp0KX/+85955JFHWLp0KStXrmTlypXs3Lmz4O3U1dUxY8aMPdO4jZyVYuS60W6nTJmixFhEJI+VK+H3v4dnnw2mAVu2DP7616EEOw4TJgzVqxa6ZNe4jtYWd9Ka5vpXkWqjpHsMW7du5corr+SHP/whq1atGvU506dPZ/78+ey///7MnTuX2bNn51ymTZumumMRkZjceSd86lPw+OO5n2MWzFwwa1Ywa8GsWUHd6sia1rGW7JPJmpqCBFZEJBeFiDy2bt3KG97wBh5++GEA2traOPbYY3n5y1/Oy172Mg466CD2339/pk6dmnBPRURk6VJ4+9uDOuWpU+H1r4eXvhQOOggOPBD23x9mzw4SbtXAiki5KenO4+tf/zoPP/wwBx98MN///vd5/etfz4QJE5LuloiIjOK224KE+/TT4ec/j2eeXRGR8VLSncPAwADf/e53Afje977H8ccfn3CPREQkn+7u4PY1r1HCLSLpo2HbHO6++27WrVvHYYcdpoRbRKQC7NoV3DY3J9sPEZHRKOnO4c477wTg3e9+t058FBGpAH19wa0u3iEiaaSkO4d7770XgBNOOCHZjoiISEEySbdKS0QkjSq+ptvMmoErgV5gBnCxu28pZZv9/f0sXboUM+OYY46JopsiIkI8MTtDSbeIpFk1jHRfDdzv7pcAtwDXlLrBZcuWsXv3bubPn0+zigNFRKIUeczOUNItImlW0Um3mbUB5wC3havuBk41swWlbPf5558H4JBDDillMyIikiWumJ2ROZFSSbeIpFFFJ93AYqDH3bcDuPsAsAwoabqR9evXA7DPPvuU2D0REcmymBhidkbmipADA1FsTUQkWpVe0z0PWD9iXS/Qlr2iu7ub9vb2PY87Ojro6OjIudENGzYAMHPmzKj6KSJSsM7OTjo7OzMPZyfZl4jFErMzWlqC202bSuyliEgRCo3ZlZ50O7BjxLomoD97RWtrK11dXQVvdPv27QBMmjSpxO6JiBQvO8k0s56EuxOlWGJ2xsKFwe2DD8IHPjC+DoqIFKvQmF3p5SWrCc5+z9YCrIli45qfW0QkUrHG7Le8Jbi9805wj2KLIiLRqfSk+x6gLZyCCjNrAOYD9yXaKxERGU2sMfsVr4DWVli1Cp56KootiohEp6KTbndfC9wOZK5gcyJws7uvTq5XIiIymrhj9oQJcNJJwf1f/zqKLYqIRKeik+7QhcA7zewzwNuAsc+2ERGRpMQas087Lbi96aYotyoiUrpKP5ESd98IXJB0P0REZGxxx+y3vQ0mTYI//hGWL4cFC+Lak4hIcaphpFtERASAyZPh7W8P7v/iF8n2RUQkm5JuERGpKmedFdzeeGOy/RARyaakW0REqspb3xpcKKerC5YuTbo3IiIBJd0iIlJVmpvhfe8L7l9zTbJ9ERHJUNItIiJV58MfDm5//GPYti3ZvoiIgJJuERGpQi97GbzmNbBlC/zsZ0n3RkRESbeIiFSpiy4Kbr/6Vdi9O9m+iIgo6RYRkar07nfDQQfBc89p+kARSZ6SbhERqUr19fCpTwX3v/QlGBxMtj8iUtuUdIuISNU67zyYNw8eewxuvjnp3ohILVPSPYpDDz2UU089lYULFybdFRERKcHEiXDppXDyyXDYYUn3RkRqmbl70n2IXXt7u3d1dSXdDRGRopnZEndvT7of5RR1zHYHs8g2JyKSU76YrZFuERGpakq4RSQNlHSLiIiIiMRMSbeIiIiISMyUdIuIiIiIxExJ9yg6OzuT7kKiav34Qe+Bjr+2j79S6ecWPb2n0dN7Gr1KeU+VdI+iUn54can14we9Bzr+2j7+SqWfW/T0nkZP72n0KuU9VdItIiIiIhKzmpin28y6gRVFvGQ20BNTdypBrR8/6D3Q8afn+A9w99akO1FO44jZGWn6uVULvafR03savTS9pzljdk0k3SIiIiIiSVJ5iYiIiIhIzJR0i4iIiIjETEl3FjNrNrNvm9mXzOy7ZjY16T4lwcy+YGYeLo8m3Z+4mdlJZvZHM1uQta5mfhdGO/5wfdX/HpjZ6Wa2zMzWm9m3zKw+XD/bzK4xsyvM7Btm1pB0XyW3Wvp7jdNof/N6b4tX7GeK3uOxFfs5ldb3VEn3cFcD97v7JcAtwDUJ96fszGwywQkJJ4XLGcn2KF5mNgeYAhw9oqkmfhdyHX8t/B6Y2XzgNIJj+wTwIeDjYfNNwHXufjHwDHBZIp2UQtXE32uc8vzN670twjg/U/Qe5zHOz6lUvqc6kTJkZm3AMmCGu28PR7w2A4vcfXminSsjM8skHf/u7jsT7UyZmNkEYDdwoLsvr7XfhZHHH66r+t8DM3s98KC7D4SPrwQWAZcDv3T3tnD9POBpYI67b0uqvzK6Wvt7jctof/N6b8enmM8UoC9Xm97jIcV8TqX591Yj3UMWAz3uvh0g/CBeBhyfZKfKycwMeA/wdWCdmZ2dcJfKwt0HR6xaTA39Low8/lr5PXD3/8kk3KHVwErgTQQ/78zzVod328vYPSncYmro7zUOef7mF6P3tmhFfqbka5NQkZ9Ti0npe6qke8g8YP2Idb1AWwJ9SYQHXkPwdc23gZ+Z2VsT7lYSavp3oYZ/D14N/Ds1/vOvQPp5lSjP37ze22jkex/1Ho/DGJ9TqX1P65PuQIo4sGPEuiagP4G+JMrdNwKXhv9Jfhy4I+EulZt+F6it3wMzOwTodvdHzUw//8qin1dERvmbvwu9t1EY63dU7/E45ficSm1M0Ej3kNXAjBHrWoA1CfQlLf4VmJ90JxKg34Xhqvr3wMzqgA8Dnw5XDfv5h8F8MrX78087/b1GL/M3r/c2GvneR73H0cj+nErte6qke8g9QJuZNQOEU4TNB+5LtFfJGgT+nHQnEqDfheGq/ffgH4Gvufuu8PFdwOFZ7QsITnbqKnO/pDD6e41e5m9e72008r2Peo+jkf05ldr3VEl3yN3XArcDJ4SrTgRuzjqJquqZWZuZnWNmE8LRvf8LfCbpfsUtPFYAg9r7XRh5/LX0e2BmlwAPA81mdpCZ/R9gE/CCmR0ZPu1k4Nvu3pdUPyW3Wvt7jUOuv3m9t+NTzGeK3uPCFPM5leb3VFMGZjGzGcCVwHJgLvBP7r410U6VkZktIqiH2gLcD1zl7k8l26t4mVkLcC7BnJ6XAd9x9+5a+V0Y7fiBfaiB3wMzu5RgesBsT7n7QjM7APgcwc9/CkECMoCkUq38vcYlX+zXe1uc8Xym6D3ObzyfU2l9T5V0i4iIiIjETOUlIiIiIiIxU9ItIiIiIhIzJd0iIiIiIjFT0i0iIiIiEjMl3SIiIiIiMVPSLSIiIiISMyXdIuNgZieY2QNmtniM5x1pZv9hZp+NcN9fjWpbIiJJMLOXmdlNUcbGHPu53czeE+c+4mJmJ5rZg2N9zkjlqE+6AyKlMLOFwBeBDcCrgZcB73P3n8a861XAUQU8byUwA6iLYqfh5Wz7o9iWiAiAmb0a+DLwBuD/AcsILijSAHzX3X8Tw26fA6YRUWzM42pgaRwbDuPx6cC/A88AD4ZNTcDfAEe7+/ISdlHo54xUCCXdUul+Bbzb3R81swnAD8qxU3d/xszWF/C8LWb2QoS7Ph74XYTbE5Ea5+4PmdkNwBHu/vbMejM7BfhPM/tHd7864n1ujzg2AmBm7wAezyS77v7rqPeR4e79wI1mdhFwr7tfmtWPW8azTTPbB3iju9/o7k+bWU9E3ZUUUHmJVCwzmwkcBmwDcPdBgkt3l+syq4MRP68QbwTuiXB7IiIAAyNXuPt/AZ8Hvmlmc2PYZ6Sx2sz2IxjZLre93jvgbmB1MRsxsybgx0Bz1uooPz8kYUq6pWK5+wbgaeBOM3tNuG4F8Cczm25ml5vZn8zsVDNba2ZPmtkrM683s7PN7DNmdq+ZfSlrfbuZfd7MbjSzX5rZtHD9RDP7NzP7nJl9n+Cr0aKY2afCPr3bzF5tZr8Kt3d52Mcbzaw1XL/ZzP55xCYa3X1XuK3/a2b/YGa/N7OvFP0GioiM7fsEZSZnApjZwWZ2mZl9z8z+O0x0MbOmMJZ9xsz+x8zOymzAzD5iZp81sz+Y2Udz7Wi0mGxmZ4Qx83wze9bMvmJmbWb2gzCe/o+ZvSncxDnAfsDFZnZaeO7N3WZ2XtY+LjCzL5rZj83sR2Y22cymmdk/m1mXmb3JzJ42s+fNbMF43rDwPVkcjoRjZqeb2bfM7OtmdqeZTQmP4Vtmdlv4efNC+B6/BDjTzP4ha5Mzws+iXjO7IMc+J5jZFWbWEe7Dzew/Mz8fSQl316KlYhdgEbCCYDTgOmCfcH0dQQDbQFBbNx24C3iC4J/N1wKXhs+dQTBafgowFbgha/u/B74a3v8m0BHenwPsJgisY/XxhwSjRfsQ1P5NCtc3Ab8BbgHaCD4s+oF/BiYBrwofZ56fqVfPHPfV4f0W4NNJ/yy0aNFSuQtwPvDXHG3rgH8N4+ovgbpw/fXATeH9bwCnh/f/D9Ad3j8bODe8f1gYq18aPv4h8Pnwfq6YPBPYBFwFHA2cCHwN+Fz43E8Ct2X11YEF4f1m4Eng/PDxW4Drs557O3ANQantu8LPi9PD47wT+GKB79294WfFV4CvEtR3n5jV3gPMD+8/QfCZ1AR8BvgrcAzwwTDu35vpb/j85cDXw7azgZU5+vAh4Pbw/pzwfX5z0r9XWoYvqumWiubuT5jZS4EvAB8B3mFmJ3pQ490NbHH3mwHM7PPA/cBBBB8wg2b28XBTvyEI7qcAs7LWPxm+dgpwYfha3H2dma0qoquzgW8DH3T37eE2dpnZWmCZu68J97MWuM+DeseHCT4MZhOckHkKwQcEwHbg/Wb2DPBdgn84RETi4ASJ6NHAAcBFZgawHphmwYMPEwwYAPyIoLwCgljbNSLW7gs8PmIf5zNKTHb3DWa2GbjF3f8EYGYrgW1mti/ByPCUUTvtvsPMXsxa9dFwuxnfA/4jXL+J4PPil+E+usJ+FuoeD2u6zexGgoGejDPdfaWZHUeQbE8J4/8a4C/u/kfgj+FrR9v2beFnwkN5+vRKYEd43OvMbAPDy1QkBZR0S8Uyszpgfw9OmPlYWPLxa6CTYORgpCfC2+nAfOAn7n59uO6qcJsXEyTBV43Y16sIgmXvOLt7GPA64GDg0az1I2sad+9pcB8MA3CmDGyWByU1uPsyMzsf+A7wCeC9wIuIiETIzKYTfEv3BEHcXD9KfJxDMBLbCODuuwkGCghf8y/u/kD4eNhrs4wak0PO8Fj5V+Bi4FngD8CBeQ4h+3WHAr/NevwXgtKZfdg7Fg8w/plV/kzwz0nGJjO7AriJoO+ZzHrkcY1lgNx52x3Ad81sIsHPYgfBqLmkiGq6pZLtqTMEcPdHgUuAI3I8v5EgwD0HvACckd1oZseE608xs+YR67eEDxeOs68PAN8CfhYGxaKYWRtZJ+WY2TyCmVsOIzix8vocLxURKcV5BN+s3UQQH48LR5iBPfGxB+gDTs5a/1Iza2FErDWz+nAQY6RcMXk0/w4s92Bq2F1FHMtKgpHxPbsANhOUz0TG3QfDgZFF4T8tvwYuc/euKPczYp+3AbcRfCN7NnC8u2+Oa38yPkq6pdJdFCakGfMIavEyZppZY3h/MUE93ybg58C7zOw7ZrY4HIXYRhAcJwO3m9lbwpNWjiBI1B8GLjezRjObTVBLvY+Z1ZnZp83s3Bx9tHD5HMEH0zdGaRv5/JGP305Qf5hxKHCau28BLhrlNSIixdhrBNXM3gD8C/Bhd+8mGFVeB/zazN4Rniz5+nBk+waCWU5OM7MTgPPcvZcg1n7MzP7FzI4nGHzIlOZlx79cMTkje9T5FUCrmU0CjgOazeygsK2f4MTDzOBL9j7+HfgbM5sRPj4GuNbdBwjyoVFjsZm9JuxXU473bgIj8qnwH5NTCL7dnEFQtng4sD8wKeskzZGj6X2j9H94p0apQQl/FlOB/yYYzS/mnxEpE5WXSKXbH3jCzP6L4Ks3A/42q30A+IyZbSIYFb4IwN3vMrOPEXxFeRrBiYhLAczsncC/EXyIXAd80t3dzN5NMJ3TE8CNBF8TvpzgYhLvJihbGXZRHjM7muAEofkEXzf+CviXsD+/AV4DLDCzgwlGYOYC7zazpxga9TkfmOzunSOO/QcWXNRigODEJRGRooVx5ByCxPBaoJvgZLxpwFvd/UEAd+8L4+O1BLHuFoJaboC/Jxiw+DHBN3uZGUO+TxD7LiQog7vQ3V8Mk8rXAAeb2U9zxWQzO4OgjvlDZvaUu68mOKnzawQJ94+BdwIvBZ4HfkIwKv9+C2a1WgS8zcx+7e6/NLNDgZvN7PcE35Z+JhyRPxOYGx7fEwQXCtrXzBaF/XwPQZ36r7LetwaC2H8UMC9Myp2glvr0cJuPAPcBDxGUzNxBENNvJzh58ygzO83dM9u9keBCRVvDJH0ucIaZLSMYwSZ8/chrUqwkOFH0NGAiMMHM7iY4mbJc0+jKGEw/C6lWFlw694fuvqBM+/uku+sS7SIiVSYsdWly91RenMzMPkkwk8zy8PFEgm8prsicCyTJU3mJVLPRSjfi2VHwNeyfy7EvEREpn/BE0VemOOGeCFwKLMhaPRvYpoQ7XVReIlUprKc7l+DrwjPd/aaYd/mX8GtPERGpLvXu/t2kO5GLu++04FL0/xqeZL+GoBTm0vyvlHJTeYmIiIiISMxUXiIiIiIiEjMl3SIiIiIiMVPSLSIiIiISs5o4kXL27Nm+YMGCol/X3d1Na2tr9B0qIx1DelTDcVTDMUBlHceSJUt63L0yOhuRWo7ZI1XjMUF1Hlc1HhPouIqVL2bXRNK9YMECurqKv/pqe3v7uF6XJjqG9KiG46iGY4DKOg4zW5F0H8qtlmP2SNV4TFCdx1WNxwQ6rmLli9kqLxERERERiZmSbhERERGRmCnpzqOjoyPpLpRMx5Ae1XAc1XAMUD3HIcNV48+1Go8JqvO4qvGYQMcVpbJfHMfMTgK+AJzl7svDdc3AlUAvMAO42N23lNKWrb293auxHklEqp+ZLXH39gT3r5gtIlKgfDG7rCPdZjYHmAIcPaLpauB+d78EuAW4JoI2EREpgWK2iEh0khjpngDsBg509+Vm1gYsA2a4+3Yzqwc2A4uAvvG0ZUZjMjRqIiKVKgUj3YrZIiIFSs1IN4C7D45YtRjocfftYfsAQWA+voS2aKxdC69+NZx0UmSbFBGpJJUQs7c8tIUlRy/hqQ88VeqmRERik4Z5uucB60es6wXagIFxtg3T3d1Ne/vQPx0dHR2FFdCbQVcXzJ499nNFRCLS2dlJZ2dn5mHaAlD6YrZD70O9+GB5v7kVEYHCY3Yakm4HdoxY1wT0h/fH27ZHa2vr+CZAb24ObnfuLP61IiLjlJ1kmllPwt0ZKXUxu2FWAwD9PXttSkQkdoXG7DQk3asJzmLP1gKsIQju42mLRibp3jHyc0JEpGalLmY3zA6S7oH1A6VuSkQkNmmYp/seoC2cSgozawDmA/eV0BaNhgaoq4Pdu6FfIygiIqQwZtdNrcPqjd1bd7N75+5SNyciEouyJ91mZpm7AO6+FrgdOCFcfyJws7uvHm9bpB2eODG41Wi3iNSgSojZZqbRbhFJvbKWl5hZC3Bu+PB8M/uOu3cDFwJXmtlRwFwg+4yZ8bZFo7kZtm0L6rqnTo188yIiaVVJMbthdgN9a/vo7+mnaV5TVJsVEYlMWZNud+8Fvhsu2es3AhfkeM242iKjum4RqVGVFLPrZwUfZ/3rVQooIumUhprudFPSLSKSepnyEs1gIiJppaR7LKrpFhFJPSXdIpJ2SrrHorm6RURST0m3iKSdku6xqLxERCT19lwgRzXdIpJSSrrHoqRbRCT1NNItImmnpHssqukWEUk9Jd0iknZKuseikW4RkdRT0i0iaaekeyw6kVJEJPX2JN2q6RaRlFLSPRaNdIuIpN6eEym7lXSLSDop6R6LarpFRFKvrqUOazIGtw+ye9vupLsjIrIXJd1j0Ui3iEjqmRmNcxsB6FvXl3BvRET2pqR7LKrpFhGpCI1zwqR7rZJuEUkfJd1j0Ui3iEhF2JN0a6RbRFJISfdYVNMtIlIRVF4iImmmpHssGukWEakIKi8RkTRT0j0W1XSLiFSEzEh3/zpNGygi6aOkeywa6RYRqQgNc4K5ulVeIiJppKR7LKrpFhGpCCovEZE0U9I9Fo10i4hUBJ1IKSJppqR7LEq6RUQqgqYMFJE0U9I9Fp1IKSJSEepa6pgwcQKD2wYZ2DqQdHdERIZR0j0W1XSLiFSE7EvBawYTEUmb+qQ7kGFmC4GPAs8BhwKd7v6ImTUDVwK9wAzgYnffEr4mZ1tkVF4iIrKXtMbshjkN7Fy+k751fTQf3BzlpkVESpKapBv4CfBOd19tZvOB3wALgauBO939RjM7GbgGODt8Tb62aCjpFhEZTSpjtmYwEZG0SlN5yUKgJby/A5hmZm3AOcBt4fq7gVPNbEG+tkh7pZpuEZHRpDJmN+4bJt0vKOkWkXRJU9L9C+B7ZtYCnAtcBCwGetx9O4C7DwDLgOPHaItOYyOYQX8/DOjEHBGRUCpjdtN+TQDs+uuuKDcrIlKyNCXdHwF2Ag8BW939ZmAesH7E83qBtjHahunu7qa9vX3P0tnZWXivzGDKlOD+1q2Fv05EZJw6Ozv3xCtgdtL9ySGVMVtJt4iUW6ExO0013c3AzwkC81Vmtgxwgq8tszUBmdPS87Xt0draSldX1/h7Nm0a9PbC5s0wffr4tyMiUoCOjg46OjoAMLOehLuTSypjtpJuESm3QmN2mpLunwLvcfdNZuYEwfzjBGe4Z2sB1hAE91xt0Zo+Hf761yDpFhERSGnMVtItImmVivISM5sNHOXum8JVXwC2ACuAtnCaKcysAZgP3Afck6ctWtOmBbebNuV/nohIDUhzzM5Out09yk2LiJQkFUk3sAHYaWbzstatBx4DbgdOCNedCNzs7qvdfW2utsh7lykp0Ui3iAikOGbXT6mnfno9gzsH6V+vC+SISHqkorzE3QfN7DTgc2a2BJgDfNLdt5jZhcCVZnYUMBfoyHppvrboaKRbRGSPtMfspv2aGNg0wK6/7qJxdmMcuxARKVoqkm4Ad38U+PAo6zcCF+R4Tc62SGVGupV0i4gA6Y7ZTfs1sW3pNnb9dRctL28Z+wUiImWQlvKSdJsRnvuzcWOy/RARkTE17R/Wda/SyZQikh5Kugsxc2Zwu2FDsv0QEZEx7TmZUkm3iKSIku5CzJoV3K4feV0HERFJm4kHTARg5/KdCfdERGSIku5CaKRbRKRiTDwwTLqXKekWkfRQ0l0IjXSLiFSMTNK9Y9nIC2CKiCRHSXchMkm3RrpFRFKvqa0JazD61/Wze/vupLsjIgIo6S5MprxEI90iIqlndaa6bhFJHSXdhcieMnBwMNm+iIjImFTXLSJpo6S7EPX1wVUp3XWBHBGRCqC6bhFJGyXdhVJdt4hIxdgz0v28RrpFJB2UdBdKdd0iIhWj+aBmQOUlIpIeSroLlRnp7ulJth8iIjImlZeISNoo6S7U/vsHtytWJNsPEREZU3Z5ibsn3BsRESXdhTvggOB21apk+yEiImNqmNVA/Yx6dvfupm9tX9LdERFR0l2w1tbgVuUlIiKpZ2ZMWjgJgO1Pbk+4NyIiSroLp6RbRKSiTDpCSbeIpEckSbeZvS2K7aTa7NnBbXd3sv0QESlRTcRsGBrpfkpJt4gkr76YJ5vZycClwFyGEvYJwBygOdqupUzmRMq//CXZfoiIFKimYzYweeFkALY9uS3hnoiIFJl0Az8Bvgw8DGSuh27AOVF2KpUWLICpU2HtWnjhBdh336R7JCIyltqN2aCabhFJlWKT7ufd/RsjV5rZUxH1J73M4JWvhHvvhYceglNPTbpHIiJjqd2YDUw8YCLWZPSt6WNg8wD104r9yBMRiU6xNd0XmtkZo6yvjQz02GOD29//Ptl+iIgUpqZjttUZkw5XXbeIpEOx//YfB1xuZjdmrTPAgc4oOmRmrwFeBzwH3A9sA64EeoEZwMXuviV8bnOutlhkku4HHohtFyIiEartmA1MPnIy2x7bxtbHtjL1mKlx7kpEJK9ik+7LgA6gC9gdrjPgfVF0xsw+ABzs7pdmrfsBcKe73xieFHQNcHbYfHWetui99rXBbVcX7N4NdXWx7UpEJAK1HbOBlle28OL1L9K7pBcuiHNPIiL5FZt03wHc4u67slea2bWldsTMjgfeC5yUta6N4ISfvwtX3Q3cbGYLgL5cbe6+vNT+jGrmTGhrgzVrYOVKOPDAWHYjIhKR2o7ZwJRXTQFg65Ktce1CRKQgxSbdNwD/ZGa/zVpnwOnA35fYl28ADwDfNrNDgMuBA4Aed98O4O4DZrYMOB7oz9O2BIYWogAAIABJREFUvMS+5HbYYUHS/fTTSrpFJO1qPma3vKIFgK2PbWWwb5AJjbomnIgko9ik+7PAocD5BDWBEATwuZQQwM3scOAVwAfc/TEz+yRwJ/AFYP2Ip/cCbcBAnrb4HH54MIPJM8/AySfHuisRkRLVfMyun1ZP86HN7Hh2B9ue2LYnCRcRKbdi/+X/KrC/uy9w9wPDZQHBqEkpFgEb3P2x8PF3GDrZZ8eI5zYRjJjkaxumu7ub9vb2PUtnZwnnDx12WHD7zDPj34aISB6dnZ174hUwu4RNKWYDLa8KEu3eJb0lbUdEZDSFxuxiR7pXuftoUWtukdsZrR97+uLuO8zsWaCO4Az3bC3AGoIAnqttmNbWVrq6ukrsYuglLwlulyyJZnsiIiN0dHTQ0dEBgJn1lLApxWyCuu4Xb3gxqOv+UGSbFREBCo/ZxY50f2TkCjM7E7iqyO2M9BgwzczmZK0bIAjGbeE0U5hZAzAfuA+4J09bfI47Durr4U9/go0bY92ViEiJFLOBqa8Opgrc/ODm/9/encfJVVf5/3+d3vfuLJ2lQxYgLGEVSNgUiAKiMOLogIrjoM4PgqOOu6OD6LjNd3Dft4yDK+7ihoiCC6AESFhFAsgSIAtJZ+l0J93p9fz++NxKV1eqqqu6q7qWfj8fj/uoqntv33s+VZ2T05/63M/N96lERFLKtuieaWbvBzCzBjP7OvB1YO1kgnD3hwlX2V8UHbsNqAGuBa4Hzo52PQf4qbtvcvdnU22bTCzjam4OhffICPz+93k9lYjIJClnA80nN2PVxt4H9jLYdcBoFhGRKZHV8BJ3v8DMTjWzLwDnAc8CxwO56D64FPicmdUQej8uia5u/zfg42Z2POEr0VVxP5NuW/6cdx7ccgvceCNclOxmbyIihaecHVTWV9K8opnu27vp/ks3sy6YNRWnFREZI23RbWZnptj0FPB34JOEKaL+EXj7ZAJx9+2EOV8T1+8ixS0N0m3LqwsugCuvhOuvDz3eFZqCSkQKTzk7tbYz2+i+vZuu27pUdItIQYzX0/1xYBnQBYwk2X4N4Yr1+UwygZeUY4+FRYvCDXLuuQfC1aoiIoWmnJ1C65mtcDXsvlXjukWkMMYruj8MrHP3bel2MrPpNcbCDM4+G77xDVizRkW3iBQL5ewUWk9vhUroWdvDUPcQVS3ZTt4lIjI5acdFuPsN4yXvaL+f5C6kEnHyyeHxzjsLG4eISEQ5O7Wq1ipaTm3Bh5xdv9fMUyIy9TQYeaJOOSU8rlkD7un3FRGRgpv5opkA7LxxZ4EjEZHpSEX3RB17LLS3wxNPwL33FjoaEREZx6wXhwsod964E1dniYhMMRXdE1VVBa96VXj+ne8UNhYRERlX0wlNVM+ppv/pfnrX9xY6HBGZZiZddJvZCWb2XTN7p5ktmXxIJeQ1rwmP3/8+DA0VNhYRkQxM55xtFcbM88IQkx037ChwNCIy3eSip/sNwDnu/ing9NgtfqeFFSvg8MNh61a4+eZCRyMikonpm7OBWS8JQ0w6f9JZ4EhEZLrJRdH9ZuCk6PkPgOocHLM0mI0OMfnFLwobi4hIZqZvzgZmnT+LioYKeu7soW9DX6HDEZFpJOui28xebGb/HD0/GTjM3TcBuPuIu3fnOMbidsEF4fGaa+BLXwp3qBQRKRLK2WNVNlaqt1tECiKrotvM/hv4MfAyAHe/C7jUzF6dh9hKw4oV8Na3wsAAvPnN8JKXwE5NRyUihaecndycV8wBYNsPxp3SXEQkZ7Lt6T4dOAi4LW7db4CrcxZRqTGDz34WfvxjmDEDbrgBLrxQPd4iUgyUs5OY+eKZVLVVsefuPXSvnVYd/SJSQNkW3Te5excQP8HptLudcFIXXQT33Rfm7v7LX+DuuwsdkYiIcnYSlfWVzL9sPgCbPr+pwNGIyHSRbdH9hJldCRxhZq8ys58CbwI+lvvQStCiRXD00eH5JiVyESk45ewUOt7UARWw7Yfb6H+2v9DhiMg0kFXR7e4/AP4AVAKXAl3AC939S3mIrTTNnRseL70U3vGOcMdKEZECUM5OrX5JPbNfOhsfdPV2i8iUyHr2Ene/w93f4O7nu/v/B/zRzM7OQ2yl6dOfhrPPhp4e+MxnYOlSePnLYcOGQkcmItOQcnZqC/9jIQAbP7+Rgc6BAkcjIuWuKt1GM/sl0DTOMWYD24Hf5yqoktbREW6Uc8898LnPwQ9+AD/7GfzpT2Gs97JlhY5QRMqUcnZ2Wk9tZeb5M9l5w042fGgDh3/x8EKHJCJlLG3RDTwQ7fMIMAJcDtwLrIvb52BgZl6iK2Unngjf+hZcfTW8/vXw29+G6QTvuANmzy50dCJSnpSzs3TI1Yew87c72fzlzcx73TxalrcUOiQRKVPjFd2fACy6+h0zO9TdP5C4k5npdoypzJ8P110HZ5wRer/PPBNuvDFcdCkiklvK2VlqOraJhW9fyDOffIZHr3iUk+46Cau0QoclImUo7Zhud98dS96RuWZWGb+Pma0ETs5DbOWjoQF+9Ss45hhYvx5OPRXuuqvQUYlImVHOnpjF/7WY2oW17LlnD8985plChyMiZSrbCyl/AtxvZl81s0+Y2fXATcDncx9amenogNtug5UrYcuW0ON97bWFjkpEyptydgaqmqo4/CthPPeTVz5Jz709BY5IRMpRtlMG3gS8kHARzhGE6acucvf/yUNs5aetLYztvuIK6O+H17wG3vteGB4udGQiUoaUszM364JZdLyxAx901r96PcO9yssiklvjjek+gLtvBq6KvTazCjM7yd0nfQtGM6sB1gJvdfc/mVk98HGgB5gBvMfdu6N9U24rajU18JWvwHHHwVveAh/7GKxbB9/9LsybV+joRKTM5DNnR8crm7x96CcOpeuPXfSu7+WRyx9h2XeXYabx3SKSG1kV3WZ2TZLV7cAw8I85iOfdwJK4118GbnT3H5rZi4CvAZdksK24mcEb3whHHgmvehX8/vfwnOeE4SZna/pcEcmNKcjZUEZ5u7KhkqN+dBT3nnYv2763jabjm1j0H7roXURyI9sx3YcBTwNPxS21QC56uU8HtgC7otcdwKuBX0W73AxcaGZL0m2bbBxT6gUvgPvvD+O8t26Fc88Nc3u7FzoyESkPecvZUJ55u+mYJpZ9N9xP4Yn3PsH2X24vcEQiUi6yLbovc/cPuvuHYgvwFqBxMkGYWSNhnGF8r8xKYLu79wK4+xDwJHDWONtKy/z54WY6H/hAKLbf9rYwveAPfgB79xY6OhEpbXnJ2VDeeXv2S2ez5CNLwOGhVz5E161d4/6MiMh4sr2Q8pEkq/uByyYZx3uAqxPWLQB2JKzrATrG2XaAzs5Oli9fvn9ZvXr1JMPNscpK+NCH4Hvfg1mzwp0rL7kkPD//fPjqV+EZTWMlMl2sXr16f74i3EFyQvKYsyGPebsYcvbi9y1m/hXzGdk3wl9f8ld67taMJiKSXKY5O9sx3U8C8WMfKoG5wHcmEmR0zPOAde6+LWGTA30J62qBweh5um1jtLe3s27dumSbissll4Qi+9pr4TvfgTvvhN/8JiwQbiH/wheG5ayzoHHSnVUiUoRWrVrFqlWrADCzCY9vyEfOjo6b17xdDDnbzDj8S4cztGuIzh91ct/z7+OYnx/DjBfMKGhcIlJ8Ms3Z2c5echPwPUaT+Aiw1d0fnUiQkXcBJ8RdIT4D+AXw6eh5vGZgc3T+VNtKW2truMjyjW8M47x//etwY52bbw431lm/Poz7rq6G5z53tAg/4QSoyHa0kIiUuXzkbJgmedsqjWXfWQYGnT/s5IEXPcARXz+CeZdqpikRyZ55FhftmdmZ7n5rwrojgGF3f2xCAZi1E3o7YtYA7wDuBB4C2t29z8yqgT3AIYQr7x9Lts3dNyWeY/ny5V7oXpNJGxyEO+6A3/0uLOvWwcjI6PbZs8NY8NNOC8tJJ0F9feHiFZGcMLO73X35BH825zk7OkZe83ax5WwfcR57x2Ns+lwIc/4V81n62aVU1lWO85MiMt2ky9nZdo2+KMm6rcCPs44q4u6d7r4xthASc6e7Pw1cD8Tm0DsH+Km7b3L3Z1Ntm2gcRa+6OhTVH/lIGHbS2Qk//jFcfjksWgTbt8PPfgb/8R9hv5YWWLEC3v52uO462Jb4LbCITAM5z9kw/fK2VRhLP7OUw792OFZrbPnaFu45+R6N8xaRrGTU021mbwY+SvgqMJmb3D1Zcs8+ILMNwOuimyzMINxIYQMwD/hPd98T7ZdyW6Ji6zXJOXd47DG4/XZYsyYsDz44ticcQu/3xReH5ZBDChOriGRlIj3dU5mzo/NtIId5u5hzds89PTz0yofoe6wPKmHRuxex+AOLqaxXr7eIpM/ZGQ8vMbPjgNcCn0/YNODuWyYXYn4VcwLPm56e0CP+5z/DbbeFQrwv7homFeAiJWGiw0uUs/NnuHeYJ9//JBs/sxEc6g6u49BPH8rsl87WHSxFprmcFN3Rgard/YArzc1sUfS1YlEq9gQ+Jfr64MYbw5CUX/0K9sR1Lh19dLgxz7nnwplnQlNT4eIUkTEmOaZbOTuPdq/ZzaOrHmXvg+GeCjPOncHSzy6l8SjNLCUyXU246DazdwE3uPtD0ev3EaacGrMb8Hx3X5mbcHOvVBL4lElXgFdXhwsxY0X48uVhHnERKYhsim7l7Kk3MjTC5q9uZsP7NzDUNQQVMO/SeSz+wGLqD9bF7CLTTbqcPd6UgQcDbXGvTyTcyexZRqegMsLYPCkV9fXwspeFpb8/zIpy001hWbcObr01LO9/P7S1hdvVx4rwQw8tdPQikppy9hSrqKrgoDcfxJxXzWHD+zew5etbePabz7L12q3Mv2w+i9+3mNoFteMfSETKXrbDS04EHnT3gYT1R8V6VopRKfWaFNyuXfCHP4QC/Oab4fHHx24/+OBQfJ9zThiKMnduYeIUmSYmObxEOXuK9T3ex4YPbWDrd7eCg1Ubc18zl4XvXkjjMg07ESl3uZwy8NgkyfsIQN2f5WLGDPinfwq3nn/sMXjiCfja1+Cii8K2J5+E1avhFa+AefPg8MPhssvgW98K+2bxR5yI5J1y9hSrP7SeZd9exooHV9D+inZ82Hn2G8+y9qi1/PWlf2X3X3YXOkQRKZBMpww8g/CV5CpgdcLmmcDn3X1R7sPLjVLuNSkqw8Nwzz2hF/yPfwwzouzdO3afjo7QA37GGWE5+mjdKVNkEiY4ZaBydpHofayXjZ/ayJZvbMH7w/+3zcub6fi3Dua8ag6VDbpmRqScTHr2EjM7Fvg2YbzgroTNA8C17v7hyQaaL+WUwIvK4CDcd18Y/33bbWF6wh07xu4zYwY873mjRfhJJ4WLNUUkIxMsupWzi8zA1gE2fmEjm7+8maFdQwBUtVUx97Vz6XhDB41HauiJSDnI1TzdrcDF7v71XAY3FcoxgRelkRF4+OHRIvzWW2HjxrH7NDTAqaeGQvz00+GUU8LFmiKS1CTm6VbOLkLDfcN0/qiTTV/ZRM+do3e0bD2zlXmvnUf7Re1UtYw3x4GIFKuczdOd5gTPc/c/T/pAeVLOCbyoucNTT4UCPFaEP/LI2H3MYNmyME3haaeFQvyIIzQkRSQymQsp0xxTObsI9NzTw+avbGbr97Yy0hvuIFxRV8Hsl81m7qVzmXHODCqqlAtFSslk5un+JTDenVIMONrd50w8xPyaLgm8JGzbFoah/OUvYUz43XfDwMDYfdraQm94rBA/5RRoaSlMvCIFluU83crZJWioe4jOn3Ty7LefZfctoxda1syvYc4r59B+cTstp7ZgFbrbpUixm0zR/QmgC9jE6ByvB+wG/IO7XzTZQPNluiXwktLfHy7OXLMmLLffDps3j93HDI45JhTgp54KK1aE3nHdtEemgSyLbuXsEte3oY+t39nK1m9vpe+xvv3raxbU0H5RO3MunkPLaSrARYrVZIru2UCPu/ePc4JZ7r4j3T6FNJ0TeMlxh2eeGS3C16wJRfnQ0Nj9GhvDRZknnxyK8BUrYMmSUKCLlJEsi27l7DLh7nTf2U3njzrp/Ekn/c+MfqQ1HTW0/1M7sy6cRduZbVTUaAiKSLGYijHdlxXzxTpK4CWury8MQ1mzBu66C9auDWPFE82eHYrv+EJ8TtF+gy6SkTyN6VbOLiHuTs9dPWz78bZQgD81WoBXtlQy87yZzHrJLGadP4vqWZodSqSQ8lZ0m9ki4E3Am929aOc7UgIvQ9u2heI7VoSvXQvbtx+43+LFYwvxk06C5uapj1dkgnJZdCtnlz53p2dtD9t/vp0dv9rB3gfj7pVQAa2ntzLrJbOY+aKZNB7biOnbP5EpldOi28wqgQuBy4Fzga3RcRZMNtB8UQKfBtxhw4bRAvyuu0LveOLNe2KzpcSK8BNPhOOOC1MZihShyRbdytnlre/JPnZcv4Mdv9pB15+68MHR/9Or51Yz89yZzDh3BjPOnUHt/NoCRioyPeRqnu6DgcuA1wOzgKeAtwC/BU5z99tzE27uKYFPU8PDsH792EL8gQfCTX3iVVSEaQpPOGF0ec5zYNaswsQtEmcS83QrZ08zQ91D7PzdTnb+eic7b9rJwKaxM0M1HtPIjBeGArztzDbdDVMkDyZVdJvZxYQekhcAjwFfB74FXOnub89xrHmhBC777dsH998/Wojfey889FAo0BMtXDi2ED/hhLBOX9fKFMq26FbOFgjDUHrX97Lrpl3s/N1Ouv7UtX8ucACrNppXNNO2so22s9poOb2FqibdlEdkstLl7Ez+hfURbht8O3CZuz8aHXTyV2CKTLW6ujDv9ymnjK7btw8efDAU4LHlgQfCLCrPPAO//OXovjNnhl7w+EL8iCM0faEUE+VswcxoPKqRxqMaOeitBzHSP8LuNbvZ9btQhO+5Zw/dt3fTfXs3T/+/p7Eqo+mkJtrOCkV46/NadWdMkRzLZnjJIkLvyRLgF8BZ7v7v0TbzXEyDkifqNZGsDQ/Do4/CffeNLcZ3JJllrb4+jAuPL8aPPTasF5mkSQwvUc6WlAa7Btn9593svmU3Xbd00XNPD8R/4VcBTSc00frcVlpOa6H19FZqF9bqwkyRceT6Qsoq4B+BfwYeBP4PeKG7r55soPmiBC454Q4bN44W4LGCPNn0hZWVoQf82GPHLosX6xb3kpUcXEipnC3jGuoZYvdf4orwtT340Nj6oKajZn8B3nJaC80nNlNRq3wmEi+fUwYuBf4duNzdi3b6ByVwyaudO8f2iN93X7iAc2TkwH2bmsLdNROLcV20KSnkeMpA5WzJyPDeYbrv6Gb3mt10r+mme003Q7vG3qTMaozmE5tpOb2FlpNbaF7RTN3BdeoNl2ltKm6Oc6m7f3vSB8oTJXCZcr294QLNv/517LJ1a/L9580LPePxy5FHhrtsarz4tJanm+MoZ0tWfMTpfbQ3jANfE4rx3r/1HrBf1Ywqmpc3j1k0LEWmk7wX3blgZi8HPgW0AN8F3unuQ9Ftjf8b6AKqgfe4+2D0Mym3xVMCl6LR2XlgIf7gg6FIT6amBpYuPbAgX7o03IFT/5GVvXwU3bmgnC2DXYP03NnD7jW76VnXQ8/aHga3HfBxUt1ePaYIbzqxidoFKsSlPBV90R1d8PNR4HPA0cBXgP9y90+a2R+B97r7nWb2BmCxu/9n9HMpt8VTApeiNjICTz8Njzxy4LJxY+qfa2qCQw4ZuyxZAgsWQEdHKMo1frzkFWPRrZwtybg7/Zv6QwEeW9b2MLRz6IB9q2ZW0XR8E03HN9F4XCNNxzfRcFQDlXX6Zk9KWykU3WcAa9x9KHr9cUIi/whwnbt3ROsXAI8Ac4FjU21z9zG3IVQCl5K1Z0+YRSW+EH/0UXj8cdi9O/3PVlXB/PmhAO/oCENY2ttHlzlzRp/PmhX2l6JTpEW3crZkxN3Z99S+MUX4nvv2JC3EqYSGIxtoOi4qxo9vpOm4Jmrm16hXXErGZOfpzjt3vy1h1SagmXBzhyfj9tsU/cNbDjw3zbZb8hyyyNRoagq3qj/xxLHr3WHXLnjiibHLhg2wZQts3hwu8IzNNT4eM5gxI3lBHivKZ80K85THHltb1ZM+TSlnS6bMjPol9dQvqWfORXOA0R7xvffvZc/9e9jzwB723r+X3kd76f1bWLZ9f9v+Y1S2VtK4rJGGZQ00LGug8ajwvG5xHVapYlxKR1EU3UmsAD4BrAISJ0buATqABWm2jdHZ2cny5aN/dKxatYpVq1blMl6RqWUWCt+ZM2F5ik7QfftGC/BNm8JFnJ2dyZcdO0KRvnNn6E3PREVFKNQTi/FkjzNmhCK9rS081tTk7r0oQ6tXr2b16v0z+s0uZCwZUs6WjJkZdQfVUXdQHbMuGJ25abh3mL1/C4V4rCDf+9e9DHUN0X1HN913dI85TkVdBfVH1I8W5Ec2UH9YPfWH1uvGPjKlMs3ZRfdbGU1p1enu90d3UOtL2KUWGATSbRujvb0dfVUp005dHRx8cFjGMzwcCu/xivL4x+7u8JjshkGZxNbaOrYQj1+SrUtcX8bDYeKLTDPbXuBw0lLOllypbKikZUULLSta9q9zdwa2DtC7vnf/svehvfSu72VgywB779/L3vv3HnCs6jnV1C+tD0X40rA0HNZA/dJ6qlrLN3dIYWSas4vqN8/MKoErgPdGqzYBh8dtN6AR2DzONhHJRmVlGFIyZ07mPzM4GIa4JBbjscf457t2hTHou3dDV1fohd+3L/UUipmor4eWlrA0N48+j18yWd/QoFlgJkg5W/LNzKidV0vtvFpmPH/GmG2DXYP0PhwV4w/10vv3Xvoe62Pf4/sY3DbI4LZBum/vPuCY1bOrqTukjrrFddQtGX2sXVxL3eI6qpqLqjSSMlJsv1nvAj7p7v3R65uAN8RtXwIMAOuix1TbRCTfqquzL9QhjEfv7R0twuOL8UzWxZa+vrBMpnCHMExmosV7WxvMnRuG0UzP8e3K2VIw1W3VtJ7aSuuprWPW+0gYM973WN/Y5e/hcXD7IIPbB+m5qyfpcatmVo0pxusWh4K89qCw1MypwSr0h7pkr2iKbjO7kpB8683sEGAlcCuwxcyOcfcHgRcBX3D3AWCdmaXaJiLFygwaG8PSccBw3szECvfubujpCY+JS6r1idv6+kJx39U18TbFvimYOzfMErN0KSxbBkcdBaedBrW1Ez92kVLOlmJlFUbdwjrqFtYd0Dvu7gxsGWDfhn1heWrsY/9T/QztHGLPzj3suWdP8uNXGTULavYX4cmWmnk1VFRNyz/EJY2iKLrN7CrCVFPxHnb3a8zslcAHzGwD0AS8L26fdNtEpFzFF+7z50/uWIODo0V4poV6bNm5M/S0d3WFi1a3bDnw+DNmwGtfC+94ByxcOLlYi4RytpQqM6O2o5bajlpaT289YLu7M7ht8MBi/Jl++jeGZbBzkP6n+ul/qj/JGSIVUDMvTWHeUUvN/BoqGzQv+XRSFPN055vmfBWRvOrvh23bQgG+ZUuYAWbtWvjRj0b3aWwM6xcsyOrQxThPd74pZ0sxG943zMDmgf1FeLJl4NmBcOnwOKraqqjpqAlFeEcNNfNHn+9/nF9LRa16zUtF0c/TLSJS0mprQy92rCd70SL44AfH7lNTE4ahiEhJq6yrpP6QeuoPqU+5z8jgCANbUhfmA1sG6N/cz1DXEENdQ/Q+1Jv2nFWzqg4sxhMea+bWUFGt4ryYqegWEckV9zDU5IEHYG/CNGZvelMY7y0iZa+iuoK6RXXULapLuY+7M7RziP7N/aHnPNXjln6GdgwxtGOIvX89cHrE/SxMlRjfQ76/KJ8/2otePbda480LREW3iEg67qGAjk2DuH07PPtsuOlQ7OZDscfNm8NQk2SyHFYiIuXNzKieVU31rGo4NvV+PuJhHPmW9MX5wNYBBrcOMrh1EO5Nd2Kobq8eHcoSK8jnR89jBfu8Gg1ryTEV3SJSvkZGYM+esRdCxl8QmWrd7t3hIsnt20OhPZDFBBvNzXDQQWEGk0MPhSOPhLPOCo8iIlmyCqNmbhg+wnNS7zcyNMLgtsHkRfmWgdFl68D+ecyT3VgoXtXMqvGLc10QmjEV3SJSOMPDoWc4drOcfftCr3LismdP8vWZLLlQXx9uZx9b5s0L0x3On3/gY1NTbs4pIpKFiqqK/TOzpBMrzge2DOzvPd//PCrM+zeHi0GHdg4xtHOI3r+lH3Ne2VI5pjDfX5x3jC3UK5srsWl8MzIV3SLThXuYHq+/P/Tcxi/J1qVaH1sXXygnFs7J1iXbZ2go/+1ubEx+k5vxns+cOVpkNzTkP04RkSkQX5w305xyPx9xBrcPjhbhWwZSFurD3cP0dvfS+3D64ryioSKj4rxqRlVZFucqukVyYWgoedE6XlE70YJ3IvsODhb6XUqurm50qa0NPcWxObhzsWjGEBGRrFmFUTOnhpo5NTQdn/obPHdnaNfQAcNYkhXqI70j++8QmvbctTY6jCXuItDY87qF4W6hlY2lld9VdEtxGxkZv2gthu0jI4V+pzJTXR2mrosttbVjX6dal2x9fLEcXzSne524rro63OhGRERKkplRPbOa6pnVcEzq/dyd4e7hA4exxBfqUXE+3DO8/66h6VTPqaZuSR11B4civP7g+tHXB9cV3RSKKrplVGz4wb594dbY8Y/jrYsNH8h1UTs8XOh3JTNmoZhMVbAmK16rqw/cfzJF8Hj7qsAVEZECMTOqWquoaq2i8cjGtPsO701dnPdv7qf/6X72PbVv/wWhPXf1HHi+aqP+sHoaljXQcGQDjcsaw/NlDVTWF6aHXEV3MRoZSV7UZlIIZ7t/4vNi67GNFbOZFLET3Z6LfTSEQUREJCcqGytpWNpAw9LU19P4iDOwZYC+J/vY92ToFY899j3eR//T/fQ+1HvgjYcqoenfy/a5AAAVMElEQVTYJppPbqbl5BZaTmuhYVnDlIwhV9GdK3190NkZphgbb1qyxMe9e8cWv9lMT5Zr1dXha//6+rGPydbFb4sNG8h1MVylX1EREREZyyqM2gW11C6ohecduH147zC9j/TSu76Xvev30rs+FOC9j/ay57497LlvD1tWbwGgdlEts/5hFnNfM5fW01rzFrMqmkz09sL69fD44/DYY+Fx69ZQZG/bFh5zNTVZTCZF70SK43Tr6urUYysiIiIlr7KxkuYTm2k+cewMLcN7h+m5p4eetT1039VN15+66H+6n81f3szmL2+m5fQWjlh9BI1Hpx8CMxEqulPp74drroHvfQ/uvHP8mR9qamDOHJg9G1pbR6cdy+SxsTEUvrHit6ZGY29FREREcqyysZK2M9poO6MNCMNUeu7uofMnnWz5+ha6b+/m7lPu5sTbT6TpuNzed0FFdzJDQ3DeeXDLLeG1GRx9NBx22Ohd5jo6oL09FNrt7aF4VqEsIiIiUjKswmhZ0ULLihYWX7WYh1/3MNuv284T732C4244LqfnUtGdzB/+EAru2bPh85+HF78Y2toKHZWIiIiI5ElVcxWHf+Vwtl+3na5bunD3nF5gqaI7mb//PTy+/OVwySWFjUVEREREpkTNnBoqGisY2TvCcM8wVS25K5WLa9bwYhGbPUS3fhYRERGZVirqQnk80p/baZRVdCcTm6u6Qm+PiIiIyHRSUauie+qo6BYRERGZlipqQv3n/Z7b4+b0aOXCozdZs5GIiIiITCtWFeo/H1bRLSIiIiKSH9F9AnNddJf87CVmVg98HOgBZgDvcffuwkYlIiLJKGeLSLFTT3dqXwb+7O5XAr8AvlbgeEREJDXlbBEpalYZFd1DKrr3M7MO4NXAr6JVNwMXmtmSQsUkIiLJKWeLSCmIFd0M5/a4JV10AyuB7e7eC+DuQ8CTwFmFDEpERJJaiXK2iBS5/cNL1NM9xgJgR8K6HqCjALGIiEh6ytkiUvT2Dy/RhZRjONCXsK4WGIxf0dnZyfLly/e/XrVqFatWrcp/dCIiE7B69WpWr14dezm7kLHkmHK2iBS9bC+kzDRnl3rRvYlw9Xu8ZmBz/Ir29nbWrVs3ZUGJiExGfJFpZtsLHE4uKWeLSPGLTRmY4fCSTHN2qQ8v+SPQEU1BhZlVA4uAWwoalYiIJKOcLSJFL1/DS0q66Hb3Z4HrgbOjVecAP3X3TYWLSkREklHOFpFSkK8LKUt9eAnAvwEfN7PjgXmABv6JiBQv5WwRKWr5mjKw5Itud98FXF7oOEREZHzK2SJS7HRHShERERGRPNMdKUVERERE8i02e4l6ukVERERE8kN3pBQRERERyTNNGSgiIiIikmexnu5cz16ioltEREREJKILKUVERERE8syqVXSLiIiIiORVbHjJyOBITo+roltEREREJLK/p3tQPd0iIiIiInlRUR3KYw0vmQqNjTB/PrS0FDoSEREREZlCVW1V1MyvoaI+t2VyVU6PVi7e9KawiIiIiMi0svh9i1n8vsU5P656ukVERERE8kxFt4iIiIhInqnoFhERERHJMxXdIiIiIiJ5pqI7jdWrVxc6hElTG4pHObSjHNoA5dMOGascP9dybBOUZ7vKsU2gduWSiu40yuEXTW0oHuXQjnJoA5RPO2Sscvxcy7FNUJ7tKsc2gdqVSyq6RURERETyzNxze7edYmRmncBTE/jR2cD2HIcz1dSG4lEO7SiHNkBptWOxu7cXOoipNM1zdqJybBOUZ7vKsU2gdmUrZc6eFkW3iIiIiEghaXiJiIiIiEieqegWEREREckzFd0iIiIiInmmojsJM6s3sy+Y2f8zs6+YWUuhY0rFzD5qZh4t90frUsZfLG0zs3PN7E4zW5JJbMXapmTtiNYf8LkUYzvM7OVm9qSZ7TCzz5lZVbR+tpl9zcw+ZmafNrPquJ+Z0LapbkO07btxn8MvirUNMjnFktdyIdvcWOwmkmNKgZmdZWaPmNkuM/ts3PqS/axizKzGzO43s5XR63JoU1a1Ut64u5aEBfgG8Mro+YuA7xc6phRxNgJfBc6JlsPGi78Y2gbMBV4GOLAkk9iKsU1p2pH0cym2dgCLgG8DJwGXAnuBd0Xb/gicEj1/A/A/cT83oW0FaMMi4NNxn8PCYmyDlpz8HhQ8r+WoHVnnxmJeJppjin0BWoErgRnABcAgcE4pf1YJ7XsfsBtYWQ5tSvV/ciHaVfA3o9gWoAPoBxqi11VRolhS6NiSxPq2aKnLJP5iahvhW5b9/7FMNO5CtymxHak+l2L8bIAzgKq41x8Hfg2cCmyOW78A2BMlrgltm+o2RM8/A1wUvz1aX1Rt0DLp34GiyWs5ak/GubHQsWbQlqxzTKFjzrBd9USzv0Wv7waeX8qfVVxbTgf+FdgArCyTNmVVK+UzFg0vOdBKYLu79wK4+xDwJHBWIYNKZGYGvAr4FLDVzC6JNq0kdfzptk0pdx9JWLWSicWdblveJbYjzefCOLGm25av2G+LzhOzCXgaeEF07th+m6KnyyexLS9StcHMGoCXAD+KXp8dt09RtUEmbSVFktdyIcvcWNQmmGOKnrv3eVSpmVkjsB74EyX8WcH+tlzk7tfErV5JabdpIrVS3qjoPtACYEfCuh7CX0VFw4NTCZO7fwG41sxeTPr4i7ltE427qNqU5nOB4m/HCsJXcKX8WawAvuruve6+NIrp18BvzOzYaJ9ib4Nkp9w/s3JqXyY5pmSY2RnAb4BqoIHSb9d7gKsT1pV0myZYK+VN1fi7TDsO9CWsqyWM2So67r4LuCr6a+5twE2kj79Y2zbe+15SbUryufyGybUxr8xsKdDp7vebWbo4J7ot7+LbEFvn7luAy6Oe7zcDV4wTZ0n9+xeg/D+zsmhfFjmmlDwGfJMwbOYTwBOUaLvM7DxgnbtvS9hUFp/VBGqlvFBP94E2ES6OiNcMbC5ALNn4POGilXTxF3PbJhp3MbcJRj8XKNJ2mFkloRh9b7RqTCxRkmpMFmcW2/IqSRsSpfwciqUNMmHFngMmq+Tbl2WOKRnuviUaivEuwrCEUv6s3gVcY2bbzWw7sBD4BeFzKdU2JZNJrZQ/hRrYXqwLMI9wQUd99LqaMNh+QaFjGyfuduDadPEXW9sYe7HQhOIuhjbFtyPV5zKZNk5B/O8B5sa9Xg48Fff6YMLXbjUT3TbVbUiyfQXw35Np31T9PmnJ+rMveA7IQ5syyo2FjjOL9mScYwod6wTbdwzwu1L+rKL/qw6KW54BLiYUqCXZpjTtTFsr5fP86ulO4O7PAtcDsQuvzgF+6qMXehQFM+sws1ebWUXUS/BO4H3p4i+mtkUxAxikf9+LuU2J7Uj1ucDE25jn+K8E7gXqzewQM/tXoAvYYmbHRLu9CPiCuw+4+7qJbJvqNpjZMWZ2YbS9Bngt4StgirENMnGFzgG5lk1uLEB4Wcs2xxQqzmyYWYOZnRC36nzgc6X8Wbl7p7tvjC3AMGE40NOUaJtgYrVSXuOJKnyJY2YzCP9BbyD8NfSf7r6noEElMLOjCeOEu4E/A59194ejbSnjL4a2mVkz8Brgy8CHgS+6e+dE4y5Um5K1A5hDis+l2NphZlcBH0lY/bC7LzOzxcAHoliaCElqKPq5CW2byjYQxm9/n3A1+h3Ax9x9c9zPFU0bZPKKIa/lwkRyYzGbaI4pdmZ2MuEC7YcJ+eUud/9xtK0kP6tEZrYBeJ27/6mU2zTRWilv8ajoFhERERHJLw0vERERERHJMxXdIiIiIiJ5pqJbRERERCTPVHSLiIiIiOSZim4RERERkTxT0S0iIiIikmcqukVERERE8kxFt4iUDDNrN7NfmtlTZvbeQscjIiLpKW+P0s1xRKRkmNl5wB+BFsLd4A53952FjUpERFJR3h6loltESpKZ3Qq8oFRuHS0iMt1N97yt4SUiUnLMbCHwu+mauEVESo3ytopumWJm9h9m9qVCxxGTi3jM7BwzW2NmK3MUVlkxs0/m+HgGvAK4OpfHFZGxlK8nz8yWmNn/mtnXcnS8SjO73MyeymDfC8zszom2NZe5W3k7UNFdRsxsmZldZ2ZfN7P7zczN7DWFjivBDcD3Cx1EnFzE8wxwfA5iKUpmdr6ZbTezLWZ2mpktjH7P3My+aGZt0X4nm9mDZvYpM6uP1s0Fns1xSK8G/s/dh6Lji5Qc5esJKcV8vS16rM3R8SqALmBRBvuuA5al2mhmR5vZD6PfvavMbF7ctlznbuVtVHSXm58DH3L3y4ATgG/HNpjZcWZ2ViGCij+3uz/o7n8uRBzJZBNPqvfQ3R8Btucyrlx+XpM9lrvfAHwZGHL3Ne7+DHA5MADc6e5d0X53EZL8u929L/rxlwDXT6oBcczs/cAHgZvN7GHgiFwdW2SKKV9nqRTztbv3AptzdS53HwTuznDfrYQCPdX2vwGro5efdPf4IjtnuVt5e5SK7jJhZjOBw4G9AO4+AnwAcDNrJSR0K0BcBTt3LmXQjpEpPFchjvVD4CAzOwnA3XcAvwf+Ie5cbcDm6Hcv5kh3f3iS597P3T/i7oe5+3J3P9Ldb83VsUWmivJ1fhVhvs71jBXZxD/evsPRY+I465zlbuXtUSq6y0Q0/c4jwI1mdmq07ingLsI4qoOBVWb2egAz+7CZvT/6amn/uC0zW25mnzazt5lZr5k9YmbvMLN/NLO1ZvYyM/uZmXWb2RVxP5f0eAnnvtzMvhk/ti36ivULZvY/ZvZHMzs5Wp/2fHE/b2Z2hZn9zczONbMHzGyrmb0ug3McFh/POOc84D1MYkb0dXGPmV2e8J5+MHpvrouSdKbv2QHnMrOLzOwuM3udmf3dzK5Oc7xkn33SeNKJekT+Blwct7oduMCioSTAhcAv4+KsB3rjXh9pZleb2SfN7BYLw1ReaWbrzOylZvZ9M+sys7ea2UlRG3eY2fnjxSdSSpSvp0++TiZVG9Odx8wONrNrLPQafyKT8yQ57w1mdq2ZPXec/ZS788XdtZTJAhwNPEX4y/YaYE7ctg3Ayuj5McAT0fN2wl/hswh/qT8NLI+2fQm4PXpeA+wEvkiYa/MSYFO64yWemzCm7WvAN6P1dYRCri16fRHQCTSnO1+Sds+LznlltO9HgX3A4nHOkRhP2nPGv4dJYtgAfApoiH7u6Wh9C/CDuP3+QkiYGb1nKc41k/CV4WeBk4Fz0h0v4bNPGk+Gv19XAY9Hzw8F/gfoB14e9/ticftfCJwRPW8A1gIN0etbo/ehBdgNfBioB14ZvX4FoVPgnUS/g1q0lNOC8vW0yNfR9g9m+D6m+6zXAsvifsYz/D2LfZ6nAf+VZPvK6DxVceuUu/O0qKe7jHjojTyWkIQuBf5mZskuGPk78C9mVg2cEa1rIvwjXwjExuOuJ/xjwt0HgB7gJ+7eTfhHOG+c4yXG18/YCzMuAPo8GhMMXEf4x/rycc6XeNzYMb8X7ftRwnjjc8c5x5h4sjlnCr/yMH5vLTA/ro2zop6otxHeU8jwPUvR3p2E5PYLd7/L3W/O4nip4snED4FDLAwx+Sfgf4HfAq8wsyagx6OsHHkucHv0/ELCf2yx3pOXAR+J3ucu4A8exoGvA1rc/UcevnK/l+w+A5GSoHw9PfJ1EinbmOY85wCt7h6LZ12W51wB/LO7fyjD/ZW786Sq0AFIbphZJbDQ3TcAbzGz/yNc6b0aOCV+X3fvN7MFwHuAb8UO4e7bzGwdcDbhL/GDCYXW/h+Nez5ENDwp1fFShBp/jMOA6ri4RsxsA3BQuvONx933mdkTQBswO4tzTPicCYYY/be1CHjS3T+buFMm75mZ/QuhuI05wsPX0B4faxafQcp4xuPufzezewlDTFrd/Qkz+yGh9+li4Ndxccd+N2LjBRcTeqZix9oRf+i458OMNYKGwUmZUb4e077pkK/jpXwf03w2xxH+wJioFwCHmVmTu+9Jsc9I1Abl7jzSG1I+qokbb+vu9xO+vjsycUczewHwRnf/qIeZKOK9CVhhYeqqB4CPjXficY6XztOEntP4qZSMMNZxsmqAR/N8jkxsYey4Z8zslCzes18Cz4lbkl4Fn8XxksaTeXP4AfB64Im4+CqBtxK+io05hTA+NWYT8Dwza4w7b9pxhSJlTPl6rOmUr1O2Mc15uqOfqWFiPkn4luBzKbbv9tEL4JW780hFd3n5dzPriHu9ALgxej5AuHDkSML0VK1mVmtm50Tb55jZbMLYwv8jfDV0O9HXlZEKEv66NzMb53iJ57a4Y/wc2AO8KjpWW3S+WI9pqvOlMj/aZ0H0szdmcI74eMY7Z3w7EiXr9TBC71UjcL2ZnWfhgp3xPoP95wLmu/vDcctg3Ckq456nO1583EnjMbMaM1sdDR1J50eEr7V/AuDuPcBvgL/42FlLXkgYehITe7+/b2Gu7zcDsblaEz+DxM+5pGdSEElB+Zppk68zfR+TnoeQS2sIf5hB+FYDM+vIIHcbodf5X4CLzewV0c++08yuIvTy/z5uf+XufMrHQHEtU78QLs5wYBfwXeCbhK+nZkTb/wvYCLyU8HXR3wk9CxcQeki+SUhgvyRc1DEQHa+X8I/wxcAg4WKdhcB/Rtv/Nd3xEs79VuA2wji5o6Ntywk9pJ8ijG08Jlqf8nwp2u/Rz7+VMNzh2Lhtqc6xJD6e8c4Z/x4mnPulhAsKvxi9F7Gfe320/Szgweiz+VT0Pmf6nr00SVsvis73fWBBtC7dZzrmWCniOYhwE4fvZfC79vmE168Enp+w7pNJfu5M4CFgB3BVtO58wu/aF4AOomnTgDcQ/lP+36it/1Dof2NatORqQfl6OuXrRcCfCb31zxmnjek+6wsI3zDeQbio/VHgjVHbk+Zuwrcpg4Qe7pnAHwgF/2XAFYT5yj9FGNoS+xnl7jwuFr15IpjZHOAN7v7h6HUFofflCne/qqDBjcPMHDjYwxhJmSAze7e7T2g6KhGZOsrXEk+5uzRoeInEezNhqMHMuHUvIO4CuWIU93WWvsqaIDOrNLN/Bn5a6FhEJCPK16LcXWJUdEu8rxLGfq03s2eAm4CN7r6msGGlFl1Y8tbo5essTF0n2WsCfu7uT4y7p4gUA+VrAeXukqLhJSIiIiIieaaebhERERGRPFPRLSIiIiKSZyq6RURERETyTEW3iIiIiEieqegWEREREckzFd0iIiIiInmmoltEREREJM9UdIuIiIiI5JmKbhERERGRPPv/AYLiDZhtOY4KAAAAAElFTkSuQmCC\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/oceanus-uranus-probe.png',bbox_inches='tight')\n", "plt.savefig('../plots/oceanus-uranus-probe.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/oceanus-uranus-probe.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 }