{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 42 - Dragonfly - Titan (Planned)" ] }, { "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 Dragonfly probe at Titan. https://en.wikipedia.org/wiki/Dragonfly_(spacecraft)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"TITAN\") \n", "planet.h_skip = 1270.0E3\n", "planet.h_trap = 0.0E3\n", "planet.loadAtmosphereModel('../atmdata/Titan/titan-gram-avg.dat', 0 , 1 , 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle1=Vehicle('Dragonfly', 1000.0, 140.0, 0.0, np.pi*3.7**2.0*0.25, 0.0, 0.43, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle1.setInitialState(1270.0,0.0,0.0,7.3,0.0,-49.7,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle1.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle1.propogateEntry (120*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+17YcXAAAgAElEQVR4nOzdeXxc5XX/8c+RLO+7Jdt4w4DZAsYBFKAhCU4MhYaGOBQCAZpSAkpJkyZp4QclhCYsbSBLISRQlJS0kJJQQlhMIGzZoGyRQwwhmCXYeMO2ZFveLcnW+f1x70ijsUbSaJZ778z3/Xrd18zcZ+bOGdk+On7m3OeauyMiIiIiIsVTFXUAIiIiIiLlTkW3iIiIiEiRqegWERERESkyFd0iIiIiIkWmoltEREREpMhUdIuIiIiIFNmQqAMohdraWp89e3ZOr2lubqaurq44AZUZ/awGTj+r3OjnBYsXL25x94r6IQwmZ0N5/H1J+mdIevygzxAXSf0MfeXsiii6Z8+eTVNTU06vqa+vz/k1lUo/q4HTzyo3+nmBmb0ddQylNpicDeXx9yXpnyHp8YM+Q1wk9TP0lbPVXiIiIiIiUmQqukVEREREikxFdxYNDQ1Rh5AY+lkNnH5WudHPS3JRDn9fkv4Zkh4/6DPERTl8hkzm7lHHUHT19fWexL4gEREzW+zu9VHHUUrK2SKSVH3lbM10i4iIiIgUWUWsXiIiEoXdu3cDMGSIUm0xuUNnZ3C/ujraWEREstFMt4hIEfzzP/8zNTU1fPOb34w6lLL3zDMwZAh84ANRRyIikp2KbhGRIqipqQGgvb094khERCQOVHSLiBTB0KFDARXdIiISUNEtIlIEqaK7ra0t4kgqRwUsxiUiCaaiW0SkCDTTXTpmUUcgItI/Fd0iIkWgoltERNKp6BYRKYJhw4YBKrpLSe0lIhJnJS+6zewkM3vezGan7TvdzJaZ2QYzu8nMhqSN1ZrZbWZ2vZl9y8xqBjImIhKlcunpTkLOVnuJiCRBSYtuM5sCjAaOSds3C1gInAF8EbgQ+ELay+4Bbnf3y4DXgasHOCYiEplyaC9RzhYRKZySFt3uvg54IGP3vsAF7r7Y3e8Avgt8EMDMjgMOdvfnw+cuAj5nZqP6Giv6BxER6Uc5FN1Jy9lqLxGROCt5e4m7d2Y8fsrdd6ftWg2sCO9/CFiW9tzV4d36fsZERCJVLj3dScjZai8RkSSI44mU7wH+I7w/HdiQMb4VmNbPWA/Nzc3U19d3bY2NjQUOWUSkp3x6uhsbG7vyFVBb2MgKTjlbRCraQHP2kGwDUTCzOUCzuy8JdzmwM+Npw4COfsZ6qKuro6mpqcDRiohkl097SUNDAw0NDQCYWUtBAyuguOVstZeISBQGmrNjM9NtZtXAp4HL03avBiakPceAUcCafsZERCJVDj3dfYlTzlZ7iYgkQWyKbuAS4Bvunv5d7OPAwWmPZwPtQFM/YyIikSr3ohvlbBGRnESxTndqTsLS9l0BvAiMMLP9zewCM5vj7k3AO2Z2ePjUU4Cb3b29r7ESfRQRkazK5URK5WwRkcIoaU+3mY0Bzgsfnm9m3yH4evKajKcudffbw/tnAVeZ2XKC9WK/lPa8vsZERCJTDhfHSVrOVk+3iMRZSYtud98K3BpuKdeGW7bXvA18KtcxEZEolUN7SVJytnq6RSQJ4tTTLSJSNsqh6BYRkcJR0S0iUgQquktP7SUiEmcqukVEiiB1ImWSe7qTQu0lIpIEKrpFRIpAM90iIpJORbeISBHU1NQA0NHRgavvoST0YxaROFPRLSJSBGbWo/CW4lF7iYgkgYpuEZEiUV+3iIikqOgWESkS9XWXltpLRCTOVHSLiBSJiu7SUHuJiCSBim4RkSJR0S0iIikqukVEikRFd2mpvURE4kxFt4hIkehEytJQe4mIJIGKbhGRItFMt4iIpKjoFhEpEhXdIiKSoqJbRKRIVHSXlnq6RSTOVHSLiBRJquhWT3dxqadbRJJARbeISJGkTqTUTLeIiAwp9Rua2UnAtcBZ7r483DcCuAHYCkwALnP3LfmMiYhErRzaS5KUs9VeIiJxVtKZbjObAowGjskYugV42t2vAB4AbivAmIhIpJJedCclZ6u9RESSoKRFt7uvI0i0XcxsGnAOsCjc9QRwmpnNHuxYUT+EiMgAJb2nWzlbRKRwSt5e4u6d1nNaYj7Q4u47wvHdZrYMOAHoGOTY8hJ9HBGRrMqhpztJOVvtJSISZyUvunsxHdiQsW8rMA3YPcixHpqbm6mvr+963NDQQENDQ35Ri4j0Y7DtJY2NjTQ2NqYe1hY2qrzFLmervUREojTQnB2HotuBnRn7hhHMipDHWJe6ujqampryDFNEJDeDLbrTi0wzayl4YPlRzhYRSTPQnB2Hons1wVns6cYAawiS+2DGREQil/Se7ixim7PVXiIicRaHdbp/CUwLl5LCzGqAWcCv8xgTEYlc0lcvySJ2OVvtJSKSBCUvuq37jBwDcPe1wEPAgnD/icC97r56sGMl+BgiIv0qhxMplbNFRAqjpO0lZjYGOC98eL6Zfcfdm4GLgRvMbB4wFUg/Y2awYyIikUr6THfScrbaS0QkzkpadLv7VuDWcEvfvwm4KMtrBjUmIhK1pBfdScnZai8RkSSIQ0+3iEhZKtMTKUVEZBBUdIuIFEnSZ7pFRKRwVHSLiBRJOZxImSTq6RaROFPRLSJSJJrpLg31dItIEqjoFhEpEvV0i4hIiopuEZEiUdFdWmovEZE4U9EtIlIkqZ7ujo6OiCMpb2ovEZEkUNEtIlIkmukWEZEUFd0iIkWiEylLS+0lIhJnKrpFRIpESwaWhtpLRCQJClJ0m9mHC3EcEZFyEtf2EuVsEZHSG5LLk83sFOBKYCrdBXsVMAUYUdjQRESSLer2kkrL2WovEZE4y6noBu4E/g14EegM9xlwTiGDEhEpB6n2kghnuisiZ6u9RESSINei+y13/1bmTjNbWqB4RETKRtQz3Shni4jERq493Reb2Rm97D+tEMGIiJSTGJxIWVE5W+0lIhJnuRbdxwPfN7M9aVsncGsRYhMRSbQYnEhZETlb7SUikgS5tpdcDTQATcCecJ8Bf13IoEREykEM2kuUs0VEYiLXovsR4AF37zFtY2bfyzcQMzsU+CzwJnAg0OjuvzezEcANwFZgAnCZu28JX5N1TEQkajU1NQDs3r2bzs5OqqpKfmkE5WwRkZjItej+MfDPZvaLtH0GnA58Ps9Y7gQ+6u6rzWwW8ChwKHAL8HN3vztc/uo24BPha/oaExGJlJkxdOhQ2tvbaW9vZ/jw4aUOoaJytnq6RSTOzHPIUmb2W4IZjVYg9UIDprp7Xr9NzGw7cLS7LzWzOmAJUA8sAya4+w4zGwJsBg4D2rONufvy9GPX19d7U1NTPuGJiAzKmDFj2LZtG5s3b2bs2LE5v97MFrt7/WDeu1Jy9muvwSGHwEEHBfdFRKLSV87Odab768Aj7r414w0KcXWznxCc8PMXwHnA54D5QIu77wBw991mtgw4AejoY2x5AeIREcnbsGHD2LZtW1QnUypni4jERK4Nhiszk3doagFi+QywC/gtsM3d7wWmAxsynrcVmNbPWA/Nzc3U19d3bY2NjQUIV0Skf4M5mbKxsbErXwG1ebx9ReVstZeISBQGmrNznen+DPBs+g4zOxO4Ebg91yAzjAB+RJCYbwxnQBzYmfG8YQQzJvQz1qWurg61l4hIFAZTdDc0NNDQ0ACAmbXk8fYVkbO1ZKCIRGmgOTvXonuimX3Z3a8xs5HAt4EzCWY68vVD4Gx3bzUzJ0jmXyA4wz3dGGANQXLPNiYiEgsRXwpeOVtEJCZyai9x91OBx83sZuD3wEHAPKC3K54NmJnVAvPcvTXcdS2wBXgbmBYuM4WZ1QCzgF8Dv+xjTEQkFqJcq7vScrbaS0Qkzvqc6TazD2QZeht4A/gGsC+wEPhiHnFsBHaZ2XR3Xx3u2wC8BDwELAhvTwTuTT3HzLKOiYjEQSkvBV+pOVvtJSKSBP21l9xAsO5qK9DZy/jtBMtP7UMeCdzdO81sIXCVmS0GpgCXuvsWM7sYuMHM5hGc/NOQ9tK+xkREIlfiS8ErZ4uIxFR/RffVQJO7r+/rSWaW11eVAO6+BPh0L/s3ARdleU3WMRGROChxe0lF52y1l4hInPVZdLv7wwM5iLv/pDDhiIiUl1KeSFmpOVvtJSKSBLmu0y0iIjmI8kRKERGJDxXdIiJFpKK7dNReIiJxpqJbRKSIIl6nuyKovUREkiDvotvMjjSzH5rZP5nZ7PxDEhEpH3Gb6VbOFhGJRiFmuv8OONHdvwm8N3XhAxERiV/RjXK2iEgkClF0fxY4Orz/Y6CmAMcUESkLMWwvKducrZ5uEYmznItuM/sLMzs3vH8McGDqimLu3unuWwoco4hIYkU9010JOVs93SKSBDkV3WZ2HXAP8DEAd38B+KSZnVOE2EREEi/KmW7lbBGR+Mh1pvu9wAzgqbR9jwBfK1hEIiJlJOKZ7orK2WovEZE4y7XoftzdW4H01Jb35YRFRMpVxEV3ReRstZeISBL0eRn4XrxlZlcA083sbOBMgq8tP1fwyEREykDEJ1IqZ4uIxEROM93u/mPgF0A18EmgFfhzd/9uEWITEUm8KGe6Ky1nq71EROIs15lu3P054LnUYzOrNrMF7v5kQSMTESkDUa9eUgk5W+0lIpIEfRbdZvYgMLqfY9QCLUDZJHARkUIpZXuJcraISHz1N9P9Uvic14BO4CLgRaAp7Tn7AROLEp2ISMKVeKa7onO22ktEJM76K7q/Dlh49jtmdoC7X5X5JDN7oBjBiYgkXaroLtGJlBWZs9VeIiJJ0GfR7e6bM3ZNMbNqd9+T2mFm84FjChWQmR0HvA94E3ga2A7cAGwFJgCXpa6gZmYjso2JiMRBqr2kFDPdytkiIvGV64mUPwGWmNnTBEnzUOBkYK+ZlMEws78FDnD3K9P2/QD4ubvfbWanALcBnwiHb+ljTEQkchGfSFlROVvtJSISZ7kuGfg48OcEJ+EcTLD81Bnu/m/5BmJmJwDnAl9O2zcNOAdYFO56AjjNzGb3NZZvLCIihRLlOt2VkrNraoLbjo58jyQiUjyDWTJwDZA+q1FlZke7++I8Y/kW8Axws5nNAa4B9gVa3H1H+N67zWwZcALQ0cfY8jxjEREpiBgsGVj2OXv48OB21658jiIiUlw5Fd1mdnsvu+uAPcDCwQZhZgcDRwJ/6+4vmdmlwM+Ba4ENGU/fCkwDdvcxJiISC1EW3ZWSs1V0i0gS5NReAhwIrADeTtuGAfnOmBwGbHT3l8LH3wEMcGBnxnOHEcyY9DXWQ3NzM/X19V1bY2NjnuGKiAzMYNpLGhsbu/IVwbrag1UROTu96FZft4iU2kBzdq7tJRe6+2vpO8zsEOD8nCPcO46uWNx9p5m9QXDp4gkZzx0DrCFI4NnGeqirq6OpqSlzt4hI0Q1mpruhoYGGhgYAzKwlj7eviJxdVQVDh0J7O7S1dRfhIiKlMNCcneuJlK/1srsNuDCn6Pb2EjDOzKak7dtNkIynhctMYWY1wCzg18Av+xgTEYmFEq/T3UMl5Wy1mIhI3OVUdJvZMjN7K217G3gVuC+fINx9KfAIcEb4PuOBocD/AA8BC8Knngjc6+6r3X1ttrF8YhERKaRSrtOdqZJy9ogRwe3OzAYWEZGYyLW95HHgLoKvCSG4zPA6d3+9ALF8ErjJzIYSzH58Ijy7/WLgBjObB0wFGtJe09eYiEjkIl69pGJy9ujRsG4dbNkC++xTqKOKiBROrkX3D939N+k7zOxgM5vj7m/mE4i7txCs+Zq5fxNwUZbXZB0TEYmDKNtLqKCcXVsLf/oTtLTAwQcX4x1ERPKT6+olp/Sybx1wTwFiEREpO1G2l1BBObuuLrhtyee0UxGRIhpQ0W1mnzWzVuAyM9uTvhGsu7quqFGKiCRUTXi5xI6ODrxE69lVYs6uDRfpam6ONg4RkWwG1F7i7t8xs98AfwN8O2O43d3fKXhkIiJloKqqipqaGjo6Omhvb++a+S6mSszZqZnu9eujjUNEJJsB93SHVx273N33upCBmc1y9xWFDU1EpDwMHTq0pEU3VF7Onj07uF22LNIwRESy6rPoNrNLgIfd/Y/hrv9nZtWZTwM+CMwvfHgiIsk3dOhQtm/fTltbG2PGjCna+1Ryzp4zJ7h9M6/TQ0VEiqe/me79gPFpj48CRgFr6V6CygiWfhIRkV6U8GTKis3ZqaL7jTeijUNEJJs+i253//uMXdcBf3D3Hr85zOxdhQ5MRKRclGqt7krO2fvuG1wgZ/Vq2LABJk2KOiIRkZ5yXTJwbi/J+2DggMKFJCJSXiJcq7ticnZ1NdTXB/efey7aWEREejOgEynN7P0EX0meZGaZp6lMJDg7flGBYxMRKQulXqu7UnP2n/0ZPPUUPPMMnHpq1NGIiPQ00NVLWoE7CPoFj88Yawe+X8igRETKSQQz3RWZs084AW64AR55BK67LupoRER6Gug63S+b2XzgTHcvy2QtIlIsu3fvBrovlFNslZqzP/QhGD0aXnwR3n476PMWEYmLAfd0u/vmbMnbzN5XuJBERMrLjh07ABg5cmTJ3rMSc/bw4fDhDwf377wz2lhERDL1t073g8Dofo5hwGHA5EIFJSJSTkpVdCtnw4UXwv/+L9x2G1x+OQwZ8CXgRESKq7909BpBb+Bqutd4zWTAXxYyKBGRclLCme6Kz9kLFsDBB8Nrr8Edd8AFF0QdkYhIoL+i+3pgq7v3efaPmZXdWfAiIoXg7qUsuis+Z1dVwVVXwbnnBrdnnw0l7OoREcmqz55ud2/pL3mHPlageEREykpraysdHR2MHj26a+nAYlHODpx9Nhx9dHChnMsvjzoaEZFArhfH6cHMZpnZ9cBNBYpHRKSsrF27FoB99tkn4kgqJ2dXVcH3vhf0c998M/zsZ1FHJCIyiKLbzKrN7GNm9jDwJ+Bcgh5CERHJ8M477wDRFd2VmrOPPBKuvjq4f/bZsGRJtPGIiAy46Daz/czsOmAl8GNgDvARYCZwZiGCMbOhZrYkXF8WMxthZjeb2b+a2a1mNjbtuVnHRETi4u233wZgxowZJX3fUuTs8H1im7cvvxzOOQe2bYM//3N46aVivpuISN/6LbrN7Ewzewx4A/gr4FvADOBn7v5zDzxToHguBWanPb4FeNrdrwAeAG4b4JiISCz84Q9/AOCwww4ryfuVOGdDjPO2Gdx+O5x0EqxfH1yx8je/KeY7iohkN5CZ7p0Elw1+BjjN3b/h7s1kX45qUMzsvcA7wKbw8TTgHCB1lv0TwGlmNruvsULGJCKSr5dffhmAuXPnluotS5KzIRl5e9gwWLQIPvYxaG0Nrlr5zW+CF/ynISLSt36Lbnd/yN3/EjgP+Gszu9PMzgC6rmdsZpZPEGY2CjjD3W9P2z0faHH3HWEcu4FlwAn9jImIxMKePXt44YUXAJg3b15J3rMUOTs8RmLy9rBhwQVzLrsM9uyBSy6BE0+EP/2p2O8sItItl8vAr3D3LwN/G+6aYWbXhLMUF+UZx2XA1zL2TQc2ZOzbCkzrZ2wvzc3N1NfXd22NjY15hisi0r/FixezefNmDjjgAGbNmjXg1zU2NnblK6B2MO9d5JwNRczbxcjZQ4bA174G990HtbXwi1/A3Lnwla/A1q15H15EKthAc3bOF8gNZyd+AvzEzOYAXyRI4IPKimZ2MtDk7usz34rga9J0w4CO8H5fYz3U1dXR1NQ0mPBERAbtZ+FadSeeeGJOr2toaKChoQEAM2vJJ4ZC5+wwpqLm7WLm7IUL4fjj4YtfhP/5H/jqV+GWW+DKK4NLyOtCOiKSq4Hm7LzW6Xb3N93988Df5XGYS4DbzawlDHQmwQk2o4AJGc8dA6whuMRxtjERkch1dnZyxx13AHDmmQVbLCQvBcrZkPC8XVcHP/xhcFLln/0ZNDfD5z8PM2cGxXe4tLqISEHlVXSnuPsdebz8HODdadsa4ELgB8A0MxsBYGY1wCzg18Av+xgTEYncY489xvLly5k5cyYf/OAHow6nhzxzNpRJ3n7/++H//g/uvx+OPRY2boTrrguK74UL4cEHoaPX709FRHJXkKI7H+7e7O6rUhuwB2h29xXAQ8CC8KknAve6+2p3X5ttrNTxi4hkcneuDq/M8tnPfpaqqshTbUGVU942g49+FJ59Fp5+OljlxB0eeCDYP3Mm/P3fw5NPqgAXkfzk3NNdYhcDN5jZPGAq0DDAMRGRyPzoRz/i2Wefpba2ls985jNRh1NqiczbZkGv9/HHB+0ld94JP/gBvPpq0PN9yy0wcSJ85CPw4Q/DggUwaVLUUYtIkphXwGKl9fX1rhMpRaQU1q5dy9y5c2lpaeF73/seF154YV7HM7PF7l5foPASIS452x1+97tgxZOf/jQowFPM4OijgwvvnHRS0J6ikzBFpK+craJbRKRA2tramD9/Ps899xwLFizg8ccfJ98lsVV0x8errwZ93o89FrSitLd3jw0ZAu9+d/ds+XvfC9OnRxeriESjr5wd9/YSEZFEaGtr4+Mf/zjPPfccM2fO5K677sq74JZ4OfTQYLvsMtixI1j95PHH4Ve/gt//Hpqagu2mm4LnT5sGRx0VbEceGdzOnBnMkotI5VHRLSKSp82bN3PWWWfx6KOPMmHCBBYtWsTkyZOjDkuKaORIOOWUYAPYtg2efx6eeSZYEeXZZ2HNmmB76KHu102cGBTfhx8O73pX9zYhczFFESk7KrpFRPLwyiuvcPrpp/P6668zadIknnjiiZJd8l3iY/To4OTKBeHaLJ2dwWXmf/e7YHvxxeB2wwZ44olgS7fPPt0F+KGHwoEHwgEHwKxZUF1d+s8jIoWnoltEZBA6Ojq44YYbuPrqq2lvb+eII47gvvvuY//99486NImBqqqgcD7wQDjrrGCfO6xaFRTgf/xj9/bqq/DOO8H25JM9j1NTA/vtFxTgc+Z0386ZA/vuC8OHl/6zicjgqOgWEcmBu/Pggw9yxRVX8Mc//hGACy+8kBtvvJFRo0ZFHJ3EmVnQ0z1zJpx2Wvf+zk5YsQJeeaW7CP/Tn+DNN4P2lNdfD7beTJkSFN+zZnVv6Y8nTVIPuUhcqOgWERmAzs5OFi1axPXXX8+zzz4LwAEHHMBtt93GglRPgcggVFXB7NnBduqpPcd27IC33goK8FQhnrpdsQLWrQu2F17o/dgjRwbF94wZwYmdvW1Tp8KwYcX+lCKioltEpA8bNmzgrrvu4tvf/jZvvvkmAHV1dVx11VU0NDQwdOjQiCOUcjZyZHDS5eGH7z22Z0/QkrJiBbz9dnCbeX/zZli6NNj6MmlS78X4lCkweXL3NmFC8J8EEcmdim4RkQxtbW08+uij/Pd//zeLFi2iI7z+9+zZs/nCF77ABRdcwJgxYyKOUipddXUwgz1jRrAueG82bw6K79Wru1dTSW3vvNN9u2FDsL38ct/vOWQI1NX1LMQnT967OE9tI0YU/nOLJJWKbhERoLW1lYcffpj777+fRx55hG3btgFQVVXFySefzIUXXsjChQsZMkRpU5Jj3DiYOzfYstmzB1pa9i7K166F9et7bq2t3Sd9DsSIEcEsel9bbW3Px+PGaTZdypN+e4hIRWpra+PZZ5/lySef5Be/+AUvvPACu3fv7ho/4ogjOPfccznvvPOYNm1ahJGKFFd1dTBTPWVKcBGfvrS1QXPz3sX4unV771u/HnbuDFZsWbVq4PFUVQXrmWcW5xMnBu0t48cHt+n3U7dazUXiTEW3iFSElpYWnn/+eZ577jmeffZZnnnmGXbu3Nk1Xl1dzfz581m4cCGnnXYa++23X4TRisTTsGHdLS39cYft27tbV9K3lpbe92/YAFu2BOMtLbnHN3z43oV4b8V5b/fHjNEMuxSXim4RKTsbN27kpZdeYsmSJTQ1NfHcc891nQSZbu7cuXzoQx9iwYIFfOADH2DcuHERRCtSnsyCiwaNHh0sYzhQHR2wcePeBfqmTcHW2trzNv3+rl1BW8zatYOLd8yYoL0ltY0dm/t9nVst2ajoFpHEamtr48033+SVV15hyZIlXduqXr7LHjFiBO95z3s47rjjOPbYYzn++OOZMmVKBFGLSF9qarrbXXLhHrSzDKQ4723ftm3BLPuWLbBy5eDjHz6878I821j67ahRmnUvRyq6RSTW3J3m5maWLl3Ka6+91uN22bJldHZ27vWaESNGMHfuXObNm8eRRx7Jcccdx+GHH05NTU0En0BESsEsWGJx5EiYPj331+/ZExTcmzd336a29Mf9je3aFWzr1uX3WcaO7VmID+b+yJG6OFKcqOgWkcjt2LGD5cuXs2zZMt566y2WLVvW4/7WrVt7fV1VVRUHHHAAhxxyCPPmzeva5syZQ3V1dYk/hYgkWXV1d4/3YLkHFzTKtXDfvBm2bu0eSx1j8+b8Zt2rq/cuxgdTxA8bpuK9EFR0i0hR7dmzh3Xr1rFy5UpWrVrFqlWruu6vWLGCZcuWsbafBsxx48Zx8MEHc8ghh/S4nTNnDsN0KT0RiQmzoDVk1KjgAkODtXt3d6tL+ux7rvd37epuo8lHTc3AivT+xiu93z02RbeZnQ58ExgL/BD4J3ffbWa1wHVAK1ADXObuHeFrso6JSHG5O62traxdu7bHliqsU8X1mjVr2LNnT5/HGjJkCPvuuy/7778/++23X9dt6v7EiRMxTbPEinK2SPEMGRIskThxYn7HaW/Pv3DfvDk4uTV1Yms+hg0bfLtM+r6kXi4hFmGb2SxgIXAGcBhwK7AS+AZwD3C5uz9vZn8HXA38c/jSvsZEZBC2b99Oc3Mz69atY+3atbzzzjt7Fdapra2tbUDHnDx5MjNmzGDGjBnMnDmz6/6MGTPYf//9mT59utpBEkQ5WyQZhg4NLj5UW5vfcXbtKkzx3tbWvYZ7PkaO7FmUT5jQ/Tnr6rrvp2+TJgUz9lGKRdEN7Atc4O67gcVmdjjwQTN7GjjY3Z8Pn7cIeM3MrgXmZhtz9+2l/gAicY1JPT4AACAASURBVNTZ2cmmTZtobm6mpaWF5ubmfu+nr13dnzFjxrDPPvswdepUpk6dypQpU3oU1DNmzGD69OlqASk/ytkiFWT48GCbPHnwx0itLjPYgj293WbHjmDLdWnI8eODonz69GCt+Zkzu29nzoSDDgoK+mKJRdHt7k9l7FoNjAE+BCxLe97q8CvmeuD4PsZ+XeSQRUrK3dmxYwcbN27sd9uwYUNXAb1hw4Z+WzsyDRs2jLq6uq5COrWlF9epAnvUqFFF+sQSZ8rZIpKr9NVlpk4d/HFSF11KPxl106buCyq1tARXTU1/nFrrvbU12N54I3uMBx4YXJn1nHPgtNMGH2dvYlF09+I9wNeBBiCzg2grMA2Y3sdYD83NzdTX13c9bmhooKGhoZDxigxIe3s7mzdvprW1ldbWVjZt2jSgQnrjxo10dAyu9XX8+PHU1dVRW1tLXV1dj/u97Rs1apT6pyPW2NhIY2Nj6mGeXwyXhHK2iJRE+kWXcjlZtbMzKM7Xr4dVq4Jt5crg9rXX4De/CQr6118PtrvvhnvvhdNP7//YA83ZsSu6zWwO0OzuS8zMgczvuocBHUBfYz3U1dXR1NRUjHClgrg7u3btorW1tUfh3NvjbM/ZsWPHoN9/5MiRTJw4sWubMGFCj8fpW6qQnjRpktamTqD0ItPMBnEx7NJRzhaRuNi1K5jJTs2AZ95vbQ1mwZubu3vL168PxjMNHx60ngzEQHN2rIpuM6sGPg1cHu5aDRyUNm7AKGBNP2Mie+ns7GTr1q1dBXCqKM6831fx3N7enlcMQ4YMYfz48YwbN47x48fvVSxnK6QnTJjA8OHDC/STECkM5WwRyUVHR1AY79zZvaU/3r49uDLo9u3dW/rjvsa2bAlWaxmMoUNh331hv/1gzhyor4dTToF99ins549V0Q1cAnzD3VNLIjwO/F3a+GygHWgKb7ONSZlxd3bu3NlvwdzX/S1btuDuecUxdOhQJkyY0KNwTt/62zdy5Ei1bkg5Uc4WSQj3YP3vXbuCVUQGuqU/v7dCub/H6fdzPMUoZ0OHBidLjhsXbL3dr6sLTgidPLn7/vjxUFVV3NggRkW3mV1BkHxHmNn+wHzgN8A7Zna4u/8BOAW42d3bgSYzyzYmMdPR0bHXbHKuxfPu3bvzjmPMmDGMGzeuqxDOvJ9ZMGfu12yzSEA5W6Rve/b0LF7b2/cuaHvbN9iCeCBbnvNOeauqghEjgm348O77qcejR3dfXKi/++n7Ro0Klg+M+6/oWBTdZnYlcE3G7qXufruZnQVcZWbLgdHAl9Ke09eYFEhmW8ZgCuZclqHLZvjw4X0WzJnFc+a+sWPHai1okQJQzpY46ezsWbwOpLgdaAGcz7GKPas7GNXVwQVq0rfhw/fe19uWKpJ7K5b7KqTTHw8ZUtmXk7d8v25Pgvr6eq/Uk3L6assYaMG8devWvNsyqqqq+iyIB1I8a61nqURmttjd6/t/Zvmo5JwdZ+5BT27cittBLuxUdGbdBevQoXsXsf3tG2gxnEsBrXmn4usrZ8diplv61tnZSWtra9cFTJqbm/c6ya+v4rlYbRm5FM9ahk5EJHednUEvbOpiIKmtUIVsroVyXKWK1cEUt/m8rq9jVfqsruxNRXdE2traWLNmDatWrWLdunU9Cur169f3eNzS0pLzBU7S9deW0V/BrLYMEZHs2tp6Xi1v69a9r5zXW+Hc31jqJLQ4GTIk2kK2t31Dh6q4lWRQ0V0E7s6GDRt4/fXXefPNN1m5ciWrVq1i1apVrF69mlWrVtHc3JzTMceOHdu19nJdXR0TJ05UW4aISIF0dARXrEtdyS7ztqVl74I69TjPlUT7NWJE95X8Ro7s7pUtZiGbrbjV/IvI4KnoztOOHTv43e9+xwsvvMCLL77I66+/zhtvvMGmTZv6fF11dTXTpk1j+vTpTJ06tauYnjx5co/iOnWlQBXOIiL52bABliyBV16Bt9/uua1fP/jjDhkSLEU2diyMGRPcpu6PGtVdKKcXzplFdLax4cM1iytSLlR0D8LKlSv58Y9/zH333ccLL7zQa+vHmDFjOOiggzjwwAPZd999mT59OjNmzOjaJk+erJYNEZEi6uiABx8Mtl/8IrjcczZVVTBpEtTWBmv3pm5T9ydN6i6s07cxY4JZYBXGItIfFd052Lx5M5deeik/+MEPuk5OrK6uZt68eRxzzDHU19dz6KGHcuCBBzJlyhSdOCgiEpHnn4fzzoM33+zeN2oUzJ0bbPvvH1yBbt99Ydas4MpzmgcRkWJS0T1ALS0tvP/972fp0qVUVVVxxhlncNZZZ3HyySczZsyYqMMTEZFQSwucfDJs3gwHHggXXggf/jC8612lueqciEhvVHQPgLtz7rnnsnTpUg477DDuueceDj300KjDEhGRXjz4YFBwv+99QVtJTU3UEYmIqOgekLvvvpvHHnuMCRMm8NhjjzFt2rSoQxIRkSyWLQtuFyxQwS0i8aEv2vqxZcsW/vEf/xGA66+/XgW3iEjMtbYGtxMnRhuHiEg6Fd39uPrqq3nnnXc49thj+dSnPhV1OCIi0o9U0T1+fLRxiIikU9Hdhz/+8Y/ceOONmBnf/e53qdIZOCIisZcquseNizYOEZF0qiKzcHc+//nPs2fPHhoaGjj66KOjDklERAZg8+bgVjPdIhInKrqzeOqpp3jiiScYP3481157bdThiIjIAG3cGNyq6BaROFHRncX9998PwEUXXURtbW3E0YiIyEBt2hTc6kRKEYkTFd1ZPP300wCccsopEUciIiK5cA9udYVJEYkTFd292L17Ny+//DIARx11VMTRiIhILlJFt4hInCT+4jhmNgK4AdgKTAAuc/ct+Rxz5cqV7Nq1ixkzZjBeTYEiIgVTjJyd/b2KcVQRkcEph5nuW4Cn3f0K4AHgtnwPuGbNGgBmzJiR76FERKSngufsbDTjLSJxkuii28ymAecAi8JdTwCnmdnsfI6bKrp19UkRkcIpVs7OlOrl7uws5FFFRPKT6KIbmA+0uPsOAHffDSwDTsjnoBvD9aYmTZqUZ3giIpJmPkXI2ZmGhI2Tu3cX8qgiIvlJetE9HdiQsW8rkNcU9ebwygpjx47N5zAiItJTUXJ2JhXdIhJHST+R0oGdGfuGAR3pO5qbm6mvr+963NDQQENDQ9aD7r///px22mkcccQRBQxVRGRgGhsbaWxsTD0spwsFFCVnZzrxRJg7F0aMyCNSEZEBGmjONk/wmSZm9gngq+5+UNq+N4B/cfe7Uvvq6+u9qakpihBFRPJiZovdvb7/Z8afcraIlLu+cnbS20t+CUwLl6DCzGqAWcCvI41KRER6o5wtIhUr0UW3u68FHgIWhLtOBO5199XRRSUiIr1RzhaRSpb0nm6Ai4EbzGweMBUYeOOfiIiUmnK2iFSkxBfd7r4JuCjqOEREpH/K2SJSqRLdXlJMaWehSj/0sxo4/axyo5+X5KIc/r4k/TMkPX7QZ4iLcvgMmVR0Z1GOf9jFop/VwOlnlRv9vCQX5fD3JemfIenxgz5DXJTDZ8ikoltEREREpMgSvU73QJlZM/B2ji+rBVqKEE450s9q4PSzyo1+XrCvu9dFHUQpDTJnQ3n8fUn6Z0h6/KDPEBdJ/QxZc3ZFFN0iIiIiIlFSe4mIiIiISJGp6BYRERERKTIV3RnMbISZ3Wxm/2pmt5rZ2KhjiiszO93MlpnZBjO7ycwSv+57sZnZUDNbYmbzo44lCczsODO7xMwWmllt1PFIPCU9b5dTLi2HHJfkvGNmh5rZd83si2Z2i5m9O+qYBsLMTjKz581sdtq+RP+77k1i/2EX0S3Az939bjM7BbgN+ETEMcWOmc0CFgJnAIcBtwIrgW9EGVcCXArMjjqIJDCzvwUOcPcro45FYi+xebsMc2mic1wZ5J07gY+6++rw79ajwKERx9QnM5sCjAaOyRhK7L/rbHQiZRozmwYsAya4+45wtmEzcJi7L480uJgxs/cDz7r77vDxDQQ/p1OjjSy+zOy9wCHAVcD57v6raCOKLzM7AfgycJIrSUkfkp63yymXJj3HlUPeMbPtwNHuvtTM6oAl7j4t6rj6Y2ZVwB5gP3dfnvR/19movaSn+UCLu+8ACJPgMuCEKIOKI3d/KvVLIrQaWBFVPHFnZqOAM9z99qhjSYhvAa8CN5vZz83s+KgDktiaT4Lzdrnk0jLJceWQd34CfN/MxgDnAZ+LOJ4BcffOjF3zSfC/62xUdPc0HdiQsW8rEPv/JcbAe4D/iDqIGLsM+FrUQSSBmR0MHAl8z90/CzwJ/Dz8ClIkU7nl7aTm0kTnuDLKO58BdgG/Bba5+70RxzNY5fbvGlDRncmBnRn7hgEdEcSSGGY2B2h29yVRxxJHZnYy0OTu66OOJSEOAza6+0vh4+8ABnw0upAkxsombyc1l5ZJjiuXvDMC+BFwF3CjmZ0YcTyDVTb/rtPpRMqeVgMTMvaNAdZEEEsimFk18Gng8qhjibFLgCPNLPV4AvCAmV3n7jdEF1ZsDSEtN7n7TjN7A5gUXUgSY2WRtxOeS8shx5VL3vkhcLa7t5qZAz8ys9nuvj3qwHJUFv+uM2mmu6dfAtPMbASAmdUAs4BfRxpVvF0CfMPd26IOJMbOAd6dtq0BLiSZXyGXwkvAuIyvdXcDf4woHom3csnbSc6l5ZDjEp93wuUN57l7a7jrWmALMV+9JIty+Xfdg4ruNO6+FngIWBDuOhG4191XRxdVfJnZFcCLwAgz29/MLgi/HpU07t7s7qtSG8EZ2s3uviXq2OLI3ZcCjxAsoYaZjQeGAj+LMi6Jp3LI20nPpeWQ48ok72wEdpnZ9LR9G4DXI4pnwKz7axKD8vh33Ru1l+ztYuAGM5sHTAUaIo4nlszsSuCajN1LE37musTHJ4GbzGwowezGJzJWeBBJl9i8rVwaK4nOO+7eaWYLgavMbDEwBbg07v/5SVtpBeB8M/uOuzeT4H/X2WidbhERERGRIlN7iYiIiIhIkanoFhEREREpMhXdIiIiIiJFpqJbRERERKTIVHSLiIiIiBSZim4RERERkSJT0S0yCGa2wMyeMbP5/TzvcDP7XzP7cgHf++uFOpaISBTM7Agzu6eQuTHL+zxkZmcX8z2KxcxONLNn+/s9I8mhi+NIopnZocB1BFfieg9wBPDX7v7DIr/1SmDeAJ63ApgAVBfiTcNL4XYU4lgiIgBm9h7g34APAI8BywguRlID3Orujxbhbd8ExlGg3NiHW4A/FOPAYT4+neBy968Dz4ZDw4C/Ao5x9+V5vMVAf89IQqjolqS7H/i4uy8xsyrgB6V4U3d/3cw2DOB5W8zsnQK+9QnAbwp4PBGpcO7+WzP7MXCIu/9lar+ZnQrcZ2aXuPstBX7PHQXOjQCY2UeAl1PFrrs/XOj3SHH3DuBuM/sc8Ct3vzItjgcGc0wzmwx80N3vdvfXzKylQOFKDKi9RBLLzCYCBwHbIbgELnAVUKrLrHYW+HkD8UHglwU8nogIwF6XO3f3nwFfAf7dzKYW4T0LmqvNbAbBzHap9Xap+CeA1bkcxMyGAXcAI9J2F/L3h0RMRbcklrtvBF4Dfm5mx4X73gZeMLPxZnaNmb1gZqeZ2Voze9XMjkq93sw+YWZfMrNfmdm/pu2vN7OvmNndZvZTMxsX7h9uZt81s6vM7D8JvhrNiZn9vzCmj5vZe8zs/vB414Qx3m1mdeH+zWb2LxmHGOrubeGx/snM/tHM/s/MvpbzD1BEpH//SdBmciaAmR1gZleb2ffN7Mmw0MXMhoW57Etm9pSZnZU6gJl9xsy+bGbPmdlns71RbznZzM4Ic+b5ZvaGmX3NzKaZ2Q/CfPqUmX0oPMQ5wAzgMjNbGJ5784SZ/U3ae1xkZteZ2R1m9t9mNsrMxpnZv5hZk5l9yMxeM7O3zGz2YH5g4c9kfjgTjpmdbmY3mdk3zeznZjY6/Aw3mdmi8PfNO+HP+F3AmWb2j2mHnBD+LtpqZhdlec8qM7vezBrC93Azuy/15yMx4e7atCV2Aw4D3iaYDbgdmBzuryZIYBsJeuvGA48DrxD8Z/PPgCvD504gmC0/FRgL/Djt+P8HfD28/+9AQ3h/CrCHILH2F+N/EcwWTSbo/RsZ7h8GPAo8AEwj+GXRAfwLMBI4Onycen6qXz31uW8J748BLo/6z0KbNm3J3YDzgVVZxtYB3w7z6k+B6nD/XcA94f1vAaeH9y8AmsP7nwDOC+8fFObqueHj/wK+Et7PlpMnAq3AjcAxwInAN4CrwudeCixKi9WB2eH9EcCrwPnh45OBu9Ke+xBwG0Gr7cfC3xenh5/z58B1A/zZ/Sr8XfE14OsE/d0npo23ALPC+68Q/E4aBnwJWAUcC3wqzPu/SsUbPn858M1w7BPAiiwxXAg8FN6fEv6c/zzqv1faem7q6ZZEc/dXzGwucC3wGeAjZnaiBz3ezcAWd78XwMy+AjwN7E/wC6bTzL4QHupRguR+KjApbf+r4WtHAxeHr8Xd15nZyhxCrQVuBj7l7jvCY7SZ2VpgmbuvCd9nLfBrD/odXyT4ZVBLcELmqQS/IAB2AJ80s9eBWwn+wyEiUgxOUIgeA+wLfM7MADYA4yx48GmCCQOA/yZor4Ag1zZl5Np9gJcz3uN8esnJ7r7RzDYDD7j7CwBmtgLYbmb7EMwMj+41aPedZrY+bddnw+OmfB/433B/K8Hvi5+G79EUxjlQv/Swp9vM7iaY6Ek5091XmNnxBMX26DD/rwH+5O7PA8+Hr+3t2IvC3wm/7SOmo4Cd4edeZ2Yb6dmmIjGgolsSy8yqgZkenDDzD2HLx8NAI8HMQaZXwtvxwCzgTne/K9x3Y3jMywiK4Bsz3utogmS5dZDhHgS8DzgAWJK2P7OncU/XgHtnmIBTbWCTPGipwd2Xmdn5wHeALwLnAusRESkgMxtP8C3dKwR5c0Mv+XEKwUzsUAB330MwUUD4mq+6+zPh4x6vTdNrTg45PXPlKuAy4A3gOWC/Pj5C+usOBH6R9vhPBK0zk9k7F+9m8Cur/I7gPycprWZ2PXAPQeypyjrzc/VnN9nrtkeAW81sOMGfxU6CWXOJEfV0S5J19RkCuPsS4ArgkCzPH0qQ4N4E3gHOSB80s2PD/aea2YiM/VvCh4cOMtZngJuA/wmTYk7MbBppJ+WY2XSClVsOIjix8q4sLxURycffEHyzdg9Bfjw+nGEGuvJjC9AOnJK2f66ZjSEj15rZkHASI1O2nNyb/wCWe7A0bFsOn2UFwcx411sAmwnaZwrG3TvDiZHDwv+0PAxc7e5NhXyfjPdcBCwi+Eb2E8AJ7r65WO8ng6OiW5Luc2FBmjKdoBcvZaKZDQ3vzyfo52sFfgR8zMy+Y2bzw1mI7QTJcRTwkJmdHJ60cghBof4icI2ZDTWzWoJe6slmVm1ml5vZeVlitHC7iuAX07d6Gct8fubjvyToP0w5EFjo7luAz/XyGhGRXOw1g2pmHwC+Cnza3ZsJZpXXAQ+b2UfCkyXfH85s/5hglZOFZrYA+Bt330qQa//BzL5qZicQTD6kWvPS81+2nJySPut8JFBnZiOB44ERZrZ/ONZBcOJhavIl/T3+A/grM5sQPj4W+J677yaoh3rNxWZ2XBjXsCw/uyoy6qnwPyanEny7OYGgbfFgYCYwMu0kzczZ9PZe4u8ZVC89KOGfxVjgSYLZ/Fz+MyIlovYSSbqZwCtm9jOCr94M+Lu08d3Al8yslWBW+HMA7v64mf0DwVeUCwlORPwDgJl9FPguwS+R24FL3d3N7OMEyzm9AtxN8DXhuwkuJvFxgraVHhflMbNjCE4QmkXwdeP9wFfDeB4FjgNmm9kBBDMwU4GPm9lSumd9zgdGuXtjxmf/gQUXtdhNcOKSiEjOwjxyDkFh+D2gmeBkvHHAX7j7swDu3h7mx+8R5LoHCHq5AT5PMGFxB8E3e6kVQ/6TIPddTNAGd7G7rw+LyuOAA8zsh9lyspmdQdDHfKGZLXX31QQndX6DoOC+A/goMBd4C7iTYFb+kxasanUY8GEze9jdf2pmBwL3mtn/EXxb+qVwRv5MYGr4+V4huFDQPmZ2WBjn2QR96ven/dxqCHL/PGB6WJQ7QS/16eExfw/8GvgtQcvMIwQ5/SGCkzfnmdlCd08d926CCxVtC4v0qcAZZraMYAab8PWZ16RYQXCi6EJgOFBlZk8QnExZqmV0pR+mPwspVxZcOve/3H12id7vUnfXJdpFRMpM2OoyzN1jeXEyM7uUYCWZ5eHj4QTfUlyfOhdIoqf2EilnvbVuFOeNgq9hf1eK9xIRkdIJTxQ9KsYF93DgSmB22u5aYLsK7nhRe4mUpbCf7jyCrwvPdPd7ivyWfwq/9hQRkfIyxN1vjTqIbNx9lwWXov92eJL9GoJWmCv7fqWUmtpLRERERESKTO0lIiIiIiJFpqJbRERERKTIVHSLiIiIiBRZRZxIWVtb67Nnz875dc3NzdTV1RU+oBJQ7NFQ7NEo59gXL17c4u7J/HCDVIk5ezD0ecubPm8y9ZWzK6Lonj17Nk1NuV99tb6+flCviwPFHg3FHo1yjt3M3i5hOLFQiTl7MPR5y5s+bzL1lbPVXiIiIiIiUmQqukVEREREikxFdx8aGhqiDmHQFHs0FHs0FLtA5f0s9XnLmz5v+amIi+PU19d7OfQJiUjlMbPF7l4fdRylpJwtIknVV87WTLeIiIiISJGp6BYRERERKTIV3b356U/hqKPgmmuijkRERPqx9cWtNB3ZxNK/XRp1KCIiWVXEOt0527gRXnwxKLxFRCTW9mzfw7bfb6NqpOaRRCS+lKF6YxbcdnZGG4eIiPTLqsKcXf7rAohIgqno7k1V+GOpgJVdREQSL1Vzdypni0h8qejuTaro1ky3iEj8hUW3ZrpFJM5UdPdG7SUiIolhpvYSEYk/Fd29UXuJiEhyaKZbRBKg5EW3mZ1kZs+b2ey0faeb2TIz22BmN5nZkLSxWjO7zcyuN7NvmVnNQMbyovYSEREgITk7NdGtiRIRibGSFt1mNgUYDRyTtm8WsBA4A/gicCHwhbSX3QPc7u6XAa8DVw9wLJ9Ag1sV3SJSwZKTs8Nb1dwiEmMlLbrdfR3wQMbufYEL3H2xu98BfBf4IICZHQcc7O7Ph89dBHzOzEb1NZZ3oGovERFJTM5WT7eIJEHJ20vcvTPj8VPuvjtt12pgRXj/Q8CytOeuDu/W9zOWH7WXiIgACcnZmukWkQSI44mU7wH+I7w/HdiQMb4VmNbPWA/Nzc3U19d3bY2NjX1HoPYSEYlQY2NjV74CaqOOpx8xyNnBjXq6RSQKA83ZsboMvJnNAZrdfUm4y4GdGU8bBnT0M9ZDXV0dTU1NAw9E7SUiEqGGhgYaGhoAMLOWiMPJKjY5WzPdIhKhgebs2Mx0m1k18Gng8rTdq4EJac8xYBSwpp+x/Ki9RESkT3HK2erpFpEkiE3RDVwCfMPd29L2PQ4cnPZ4NtAONPUzlh+1l4iI9CdGOTu8VdEtIjEWxTrdqfRoafuuAF4ERpjZ/mZ2gZnNcfcm4B0zOzx86inAze7e3tdY3kGqvUREBEhIzlZPt4gkQEl7us1sDHBe+PB8M/sOwdeT12Q8dam73x7ePwu4ysyWE6wX+6W05/U1NnhqLxERSUzOtiq1l4hI/JW06Hb3rcCt4ZZybbhle83bwKdyHcuL2ktERBKUs8NbpWwRibE49XTHh9pLRESSQ+0lIpIAKrp7o/YSEZHk0ImUIpIAKrp7o/YSEZHE0JKBIpIEKrp7o/YSEZHk0Ey3iCSAiu7eqL1ERCQ51NMtIgmgors3ai8REUkOzXSLSAKo6O6N2ktERBJDPd0ikgQqunuj9hIRkeTQTLeIJICK7t6ovUREJDnU0y0iCaCiuzdqLxERSQ7NdItIAqjo7o3aS0REEkM93SKSBCq6e6P2EhGR5NBMt4gkgIru3qi9REQkOVIpu1M5W0TiS0V3b1JF95490cYhIiL9UnuJiCSBiu7eVFcHtyq6RUTiT+0lIpIAKrp7U1MT3HZ0RBuHiIj0T0sGikgCqOjujYpuEZHk0Ey3iCTAkFK/oZmdBFwLnOXuy8N9I4AbgK3ABOAyd9+Sz1heVHSLiADJyNnq6RaRJCjpTLeZTQFGA8dkDN0CPO3uVwAPALcVYGzwVHSLiCQnZ2umW0QSwErdA2dmVcAeYD93X25m04BlwAR332FmQ4DNwGFA+2DGUrMxKfX19d7U1DTwINevhylToLYWmpvz/MQiIoNnZovdvT7C9499zm5vbueZyc8wZNIQ3tfyvnw/sojIoPWVs0veXuLunV1fBQbmAy3uviMc321my4ATgI5Bji3PK0jNdIuIAAnJ2ZrpFpEEKHnR3YvpwIaMfVuBacDuQY710NzcTH199386GhoaaGhoyB6Rim4RiVBjYyONjY2ph7VRxtKL2OVs9XSLSJQGmrPjUHQ7sDNj3zCCWRHyGOtSV1dHTu0lKrpFJELpRaaZtUQcTqb45WzNdItIhAaas+NQdK8mOIs93RhgDUEKHcxYftKLbnfo+dWqiEgli1/O1jrdIpIAcVin+5fAtHApKcysBpgF/DqPsfxUVelS8CIivYtdzrbqsOpWuhaRGCt50W3dZ+QYgLuvBR4CFoT7TwTudffVgx0rSKBqMRERSUTOtiFBiL5HM90iEl8lbS8xszHAeeHD883sO+7eDFwM3GBm84CpQPoZM4Mdy09NDbS1BUX3iBEFO6yISFIkJWd3Fd27VXSLSHyVtOh2963AreGWvn8TcFGW1wxqLG+a6RaRCpeUnK2iW0SSIA493fGkoltEJBGsyoLmFwfvVOEtIvGkojsbFd0iIomh2W4RiTsV3dkMHx7ctrVFG4eIiPRLRbeIxJ2K7mxGjgxud+yINg4REemXim4RUzRXGwAAIABJREFUiTsV3dmo6BYRSQwV3SISdyq6s1HRLSKSGCq6RSTuVHRno6JbRCQxVHSLSNyp6M5GRbeISGKo6BaRuFPRnY2KbhGRxFDRLSJxp6I7m1TRvX17tHGIiEi/VHSLSNyp6M5GM90iIomholtE4q4gRbeZfbgQx4kVFd0iUqbKMWer6BaRuBuSy5PN7BTgSmAq3QV7FTAFGFHY0CI2alRwu21btHGIiAxSJeVsFd0iEnc5Fd3AncC/AS8CneE+A84pZFCxMGFCcNvaGm0cIiKDVzE5W0W3iMRdrkX3W+7+rcydZra0QPHER6ro3rgx2jhERAavYnK21QRFd2d7Zz/PFBGJRq493Reb2Rm97D+tEMHEysSJwe2mTdHGISIyeBWTs6tGBL/OOneq6BaReMq16D4e+L6Z7UnbOoFbixBbtFIz3Sq6RSS5KiZnVw1X0S0i8ZZre8nVQAPQBOwJ9xnw1/kGYmaHAp8F3gQOBBrd/fdmNgK4AdgKTAAuc/ct4WuyjuVN7SUiknwVk7OrR1QD0LlLRbeIxFOuRfcjwAPu3pa+0/5/e3ceJ1dZ5v3/c/W+pJNO0lk7CQmQzbAJQRFUIuCAOgPIDzccFZ8HIi4zOoqDC+oM6owiOjogalwfN0RBERmNyii4gwmYgCQkmI10EuhO0klv6fX6/XGf6q6uVFVvVV3VXd/363VeVXXOqXOuOnVy5eq77nMfsy9nIJZvAZe5e4OZLQJ+DqwEbgfWu/ud0ZX4XwJeF70n3bKxUfcSEZn4CiZnq3uJiOS7kRbd3wPeb2a/iptnwBXAO8cYy0qgJnreAUwzs/mEq+zfHs2/H7jbzBYDXamWufuuMcYCNTVQVAQtLdDdDaWlY96kiMg4K5icHSu6ezt6h1hTRCQ3Rlp0f4jwM+LVQGxcJiOMATvWBH4Xoe/hy4B/BP4JWAM0uXs7gLv3mNlO4HygO82yXWOMJRTcs2fDgQNhWrhwzJsUERlnBZOz1dItIvlupBdSfgpY6O6L3X1JNC0mtJqM1duAY8CfgVZ3vxuoBw4mrNcCzB9iWWbU14fHhoaMbVJEZBwVTM7u79OtoltE8tRIW7qfdveWJPPnZiCWSuAOQmL+bNQC4oSfLeOVE1pMGGJZv8bGRlavXt3/eu3ataxdu3boiBYsgI0bVXSLyLhat24d69ati72sG8OmCiZnq6VbRHJluDl7pEX324A/xs8ws1cBnwW+NsJtJfo28Fp3bzYzJyTzdxGucI9XA+wjJPdUywaZNWsWGzZsGHlEaukWkRyILzLNrGkMmyqYnN1fdGv0EhEZZ8PN2SPtXjLDzD4UbbTKzL4CfIXw8+KomVkdcLq7x+65/jHgKLAbmB8NM4WZlQKLgAeBX6dZlhkqukVkYiuYnB0bp1sXUopIvhpR0e3urwB+aWa3An8BlgGnA8nueDYSh4BjZlYfN+8gsBm4D7gwmncRcLe7N7j7gVTLxhjLgAULwuPu3RnbpIjIeCmknF1cHfXpblNLt4jkp7TdS8zsxSkW7Qa2A7cAJwCXA/8y2iDcvc/MLgc+bGYbgTnAe939qJm9FbjZzE4n9EOM79iXbtnYLV0aHrdty+hmRUSyoZBzdsm08N9Zz5GeTG1SRCSjzN1TLzT7E2Es1mYgVfOBAfPcvTzz4WXG6tWrfVR9ug8dgpkzobo6jNdtlvngRETSMLON7r566DULO2cf/t/DbLpoE7UvqeWMX52RpchERNJLl7OHupDyJmCDuz87xA7G+lNlfpoxA+rqoKkp9OuOdTcREclPBZuzS2qjlu5mtXSLSH5K26fb3X86VPKO1rsrcyHlmeXLw+PWrbmNQ0RkCIWcs/uLbnUvEZE8NdLRSwrPqaeGx7/8JbdxiIhISsXTwoWUaukWkXylonsoz3teeHz44dzGISIiKfVfSNncQ7prlUREckVF91DOPjs8qugWEclbRaVFFFUXQR/0tmqsbhHJP2Muus3suWb2bTN7j5ktHntIeWblyjB6ye7dsH9/rqMRERmTyZyzS2eWAtDdeNyd5UVEci4TLd3XARe5+6eBc2N3G5s0iovhhS8Mz++/P7exiIiM3aTN2WXzygDo2t+V40hERI6XiaL7HcBZ0fPvAaUZ2GZ+ueSS8Pjzn+c2DhGRsZu0Obt8Xhh6vHN/Z44jERE53oiLbjN7mZm9Pnr+PGBp7Da+7t7n7kczHGPuXXxxePzFL6BXfQVFZOIopJytlm4RyWcjKrrN7OPAD4BXArj7w8AbzeyqLMSWP1asgBNPhMZGePDBXEcjIjIshZazVXSLSD4baUv3ucAC4Ldx834GfCJjEeUjM7gq+j/q29/ObSwiIsNXUDk7VnR3Nqh7iYjkn5EW3b9092YgfhDUSXc74aRe//rweNdd0Nqa21hERIanoHJ2xQkVABzbeSzHkYiIHG+kRfcOM/sAsNzMXmtmdwNvBz6Z+dDyzIoVcO650NICX/96rqMRERmOgsrZlUvDQCwd2ztyHImIyPFGVHS7+/eAXwHFwBuBZuDv3P3zWYgt/7z73eHxs5/VBZUikvcKLWdXLKzAyoyuA130tOp28CKSX0pG+gZ3/xPwp9hrMys2swvd/X8zGlk+uvzycEHljh3wrW/B1VfnOiIRkbQKKWdbsVF5UiXtW9rpeKqDmjNqch2SiEi/tEW3md0LTBliG3VAEzDpEvhxiovh3/8d3vAG+OAH4dWvhqqqXEclIgIoZwNULa+ifUs77VvaVXSLSF4ZqqV7c7TOk0AfcC3wKLAhbp0lwIysRJePrroqdC/ZuBE+/vEwiYjkh4LP2VPOmELTPU20PtLKnNfNyXU4IiL9hiq6PwVYdPU7ZnaSu384cSUz+3E2gstLRUXwuc/Bi14En/xk6HJy9tm5jkpEBJSzmXJWaOhv2dCS40hERAZLW3S7+5GEWXPMrNjd+68iNLM1wPMyFZCZnQO8EHgK+B3QBtwMtADTgRtid1Azs8pUy7LqvPPgXe+C//oveOMb4eGHoUY/Y4pIbilnQ81ZIRe3PNKC9zlWZJncvIjIqI30Qsq7gE1m9jtC0lwJXAwc15IyGmb2ZuAkd78xbt7XgfXufqeZXQJ8CXhdtPj2NMuy6+Mfh/XrYcuWcEHlD34QWsFFRPJHweXs8nnllNWX0dXQRdtf25hy6lBd3EVExsdIhwz8JfB3hItwlhOGn7rS3f9zrIGY2fnA64EPxc2bD1wF/CSadT9wqZktTrdsrLEMS2Ul3HMPTJ0KP/whfDgj/4eJiGRMoebs6RdMB+DwLw9netMiIqM24qZZd9/n7je6+6Xu/o/AT8zsrAzE8hlgC3Crma03s/OANUCTu7dH++4BdgLnD7FsfCxbBnfcEVq4P/5xuPnmcdu1iMhwFGLOnv7SqOi+X0W3iOSPEXUvMbOvJZk9C+gFLh9tEGa2HHgu8GZ332xm7wXWAx8DDias3gLMB3rSLBuksbGR1atX979eu3Yta9euHW24g7385fCNb8Cb3gQ33BDmvfe9YOpHKCKjs27dOtatWxd7WTfa7RRqzp5+USi6mx9opretl+Lq4mG/V0RkpIabs0fap3spYWxXj5u3APjtCLeTaBVwyN03R69vAz4S7Sfxfr7lQHf0PN2yfrNmzWLDhg2JszPnDW+Ajg54y1tC4d3QAJ/5TBjXW0RkhOKLTDNrGsOmCjJnl88rp+b5NbQ81MLB+w4y+zWzR70tEZGhDDdnj7Tovsbdn4yfYWYrgKtHGmCSOPpjcfcOM9tOuHXx9IR1a4B9hOSeatn4W7s29O9+05vgv/873LXym9+E6YkhioiMm4LN2XNeN4eWh1p45o5nVHSLSF4Y6YWUTyaZ3QlcM8Y4NgPTzCz+TgY9hGQ8PxpmCjMrBRYBDwK/TrMsN1772jCiSW0t3HcfnHkmZLOFXUQkjULO2bNePQuK4NBPD9F5oDMbuxARGZERFd1mttPMdsRNuwkX0vxoLEG4+1bgZ8CV0X5qgTLgO8B9wIXRqhcBd7t7g7sfSLVsLLGM2UteAo88AmedBbt2wbnnwkc/Ct3H/YIqIpJVhZyzy+eVU3dpHd7t7Ptibn4AFRGJN9LuJb8EvstA/8A+4Bl335aBWN4IfM7MygitH69z9x4zeytws5mdDswF4q+mSbcsd5Ysgd//Hq6/Hm67LQwn+KMfwde+BmeckevoRKRwFHTOrn9nPU33NLHvC/tYdMMiiit1nY2I5I65+9BrxVY2e7G7/yZh3nKg192fynRwmbJ69WrP6oWU6fzqV/B//29o9S4qguuug5tugpkzcxOPiEwoZrbR3VcPvWbS9xZ0znZ3Nq7eSOsjrZz4qRNZdP2iDEQnIpJaupw90nG6L0ky7xngByOOqlBccAE89hj88z+HYQRvvx2WLoVbb4VO9TMUkawq6JxtZiz5+BIA9vzHHrqb1c1PRHJnWEW3mb3DzJqBG8ysN34ijLv6TFajnOimTIHPfQ7+8he48EI4fDgU4cuWwVe+ov7eIpJRytkDZlw8g9o1tfQc7mHnB3fmOhwRKWDDKrrd/TbgxcBngRMTpgXunqw1RRKdcgr88pfh9vGrVsGePXDttbBiBXz1q2r5FpGMUM4eYGac/LmTsRJj3+37aP5tc65DEpECNezuJdFNEN7n7rsTpv1mpo5yw2UGl10GmzeHW8gvXx7G9L7mGli8ONxO/mDiTdtEREZGOXvAlNOmsOh94SNvffNWeo705DgiESlEaS+kNLPrgZ+6+xPR6w8Sbn4waDXgJe6+JltBjlVOL6QcSk8P3HknfOpTsGlTmFdVFW6yc911cNppuY1PRHJqJBdSKmen1tfZx8bnbaRtcxt1r6xj1d2rMLOM7kNEZCwXUi4BauNenwmcS/iJckncNDcDcRamkhJ4/evh0Ufh/vvhkkugvR2+8AU4/XQ455ww1GBbW64jFZH8p5ydQlF5EavuXkXx1GKaftTE7o/tznVIIlJgRjpk4JnA4+7elTD/ObGWlXyU1y3dyfz1r/DFL8K3vgVHjoR5U6fCa14DV10FL35xGH5QRCa9MQ4ZqJydoOneJh6//HFwWLZuGfOvnZ+V/YhIYcrkkIGnJkney4GTRhucJLFqVRhScN8++MY3wl0tjx6FL3853PHyhBPgX/81jIYygj+aRKTgKGcnqLu0jqW3LwVg23XbeOa7BTOQi4jk2LDuSGlmLyL0A3ypmSWOuTQD+G/gJxmOTWJ9u9/0JnjiCfjOd+C73w032vnUp8K0ciVccQVcfnm49bz6KIoUPOXs9Oqvq6e7sZtdH97Fln/cQm9LL/PfohZvEcmuYXUvMbNTgW8S+gIeTljcBXzH3W/KfHiZMeG6l6TjDn/4QyjAv//9wSOdLFgQiu9XvjJ0QSkZ1t9UIpLHRtO9RDl7eHZ/Yjc73x/+Jjnhwyew+COLsSI1XIjI6KXL2cPu021m04BXuftXMhnceJhURXe87m544IEw7vc994TuKDHTp8NLXxouzLz4YpivVhyRiWi0fbqVs4en4QsNbH/HduiDusvrWPHNFZTUqMFCREYnI0X3EDt4obv/bswbypJJW3TH6+uDP/85FN8/+hE8+eTg5aedFgrwSy6B886DsrLcxCkiIzKWCynTbFM5O87B9Qd54rVP0Hukl8rllTznO8+h5qyacdu/iEweoy66zexeYMpQ2wdWufvs0YeYXQVRdCd66in4+c9h/Xr41a/CMIQxlZWh8F6zJkxnn60iXCRPjXCcbuXsUWrf1s7jVzxO+1/bsRJj8UcXs+i9i7BidTcRkeFLl7OH+g3tSaAZaABSVecG/P3ow5OsOPnkML397eH28r/7XSjA16+Hxx8PY4Lff39YV0W4yGShnD1KVcuqOOvPZ7Hjhh003NrAzvfvpOmHTSz74jJqzlSrt4iM3VAt3XVAi7t3pt2I2Ux3z9t7lxdkS3c6zz4LDz4Y+oM/8EAYGSVeeXkYCeUFLxiY1CdcJCdG2NKtnJ0BB9cfZNu12+jc2wlFUP+Oehb/+2JKa0tzFpOITAzj0af7mny+WCfXCTzvDVWEAyxaNFCAn3NOuFtmRcV4RypScLLUp1s5ewg9LT3s+rdd7P3cXuiFkhklnPCBE5j/9vkUVxTnNDYRyV9ZK7rNbBHwduAd7l496g1lWT4k8Anl8GF46CH44x/D9NBD4eY88UpKwk18zjprYDrttNBVRUQyJpNFt3L2yLVuamX7O7dz5MFwd+DyheUs/shi5rxhDkVlujOwiAyW0aLbzIqBS4FrgZcCz0Tbqc9AoGXAn4F3uvsDZlYJ3Ay0ANOBG9z9aLRuymWJ8imBT0i9vbBly+AifOvWMGJKvOLiwYX4mWfCKadAjfpDiozWWIvubObsaPsZz9v5lrPdnUPrD7HjfTto29wGhOJ74XsWMu+aeRRXq+VbRIKxXEgZv5ElwDXAm4GZwG7gH4CfAy/IQJwA7wUWx72+HVjv7nea2SXAl4DXDWOZZFJxcSieTzkFrr02zGtrC7eh37hxYNqyBTZvDtPXvz7w/iVL4NRTQ0t47PHkk3XzHpEsGqecDQWQt82MmS+byYyLZ/DsHc+y+z920/5EO0+96yl2fXQX9W+vZ/7a+ZTXl+c6VBHJY0O2dJvZqwgtJBcATwFfAf4f8AF3/5eMBWJ2LrAC+DBwNbAN2AlMd/d2MysBjgCrCHdUS7rM3XclbjvfWk0mrbY22LRpoAjftCn0D+/qOn7d8vLQKh4rwletCre0X7hQt7IXiTPSlu7xytnRvrKSt/M9Z3ufc/AnB9n9n7tpeaglzCyGusvqqH9bPbUX1GLKYyIFaawt3R2EZPkH4Bp33xZtdOxXYA4EWA1c6e7vNrMPR7PXAE3u3g7g7j1mthM4H+hOs2xXpuKSEaquhnPPDVNMdzds3x5avx97bOBx92545JEwJW5jxYowrVw5MJ10koYxFBmerOfsaHsFm7etyKi7rI6Zl87kyG+P0PD5Bpp+2NQ/VS6rZO7Vc5nz+jlULNIF5yISDFl0u/t9wH3RBTjXmtli4MdA/9hJZmY+tmFQbgA+kTCvHkgc0qoFmA/0pFl2nMbGRlavHvijY+3ataxdu3YM4cqwlZbCc54Tpte+dmD+kSNhvPBYd5Qnngj9xJ99dqClPF5JSSi8Y0X48uWwdGnopjJrllrHZVJZt24d69ati72sG8l7xylnQxbz9kTJ2WZG7YtrqX1xLZ37O9n/lf3s+9I+OrZ1sPMDO9n5gZ3UrqllzhvmMOvKWZRMVZc6kclouDl7NBdSlgCXA68HHge+Cvydu69L+8bU27sYKHf3e6PXuwg/U64GXuXuz49b9xHgu9HLpMvc/ZbEfeT7T5US5+DBUHxv2TIwbd0Ku3ZBqnN16tSBAnzp0sHP6+pUkMuEloELKTOas6NtZjVvT+Sc3dfTx+GfH+bAtw7QdE8T3hnylpUbMy6ewawrZjHzH2ZSOkNjfotMRhm5kDLG3XuAu4C7zOxk4F8I/QdHm8CvB54b1/9tOqFV5jPR83g1wD7CndZSLZOJbObMcHfM884bPL+9HbZtGyjEt20Lt7rfvj0MZ5isdRxg2rSBAvzkk0Nr+eLF4eLOBQvCRaIik1gWcjYob6dUVFLEzFfMZOYrZtJzpIfGuxo58K0DHHnwCAfvPcjBew9iJUbtS2qpu6KOukvrKJ+vCzBFCkGmbo7zRnf/5ijfOwuIzzh/BN4NPAQ8Acxy9w4zKwVagROBXsIFQsctc/eGxH1M5FYTGYI7NDWF4nv79oFCPL4gT6WkJNz0Z8mSwVOsKJ8zR63kknNZujnOqHN29P6s5u3JmLM793XS9OMmGu9upPmB5nA0ItWnVTPjZTOYcckMpp07TeN/i0xgGW3pTmYsydvdG+Nfm1kv0Ojue8zsPuBC4D7gIuDuWHJOt0wKiFno0z1r1uALOGFwQR4rwnfuHJj274cdO8KUTGXlQAEeK8YXLQrTwoUwd65aymVCGkvOjt6vvD1C5fPLqX9rPfVvraf7YDdNP2mi6UdNHL7/MG2b22jb3MbTn3ya4ppipl84nekXT2f6BdOpXFqpkVBEJol8v6rjrcDNZnY6MBdYO8xlIukLcoCOjjCKSqwI37VrcFF+6NBAd5ZkSkqgvn5wIR7/uGhR6N6i/zClsChvD6F0Zinzrp7HvKvn0dfZR/Nvmzm0/hCHfnaI9ifaabqniaZ7mgAom1dG7Zra/klFuMjElZHuJfluMv5UKePg6NHBBfmuXfD007BnT5iefXbobUyZcnxRvmBBKNZj09SpKswlpWx0L8l3hZyzj+05xqGfH+LwLw7T/GAz3Y3dg5aXzSuj9vxapp43lannTGXK6VMoKlV3FJF8kdHbwE9EhZzAJYuOHYO9e0MBHl+Mxz9vaxt6O9XVg4vwxGn+fJg3T3fwLFAquguXu9O+pZ3mB5r7p8QivKiyiJrVNUw9ZypTXxAK8fJ5ujBTJFey3qdbpCBVVIQRUU4+Oflyd2huPr4Q37sXGhoGpra2MBrLtm2p92UWLuxMV5zPmwe1tWo1F5kkzIzq51RT/Zxq6t9WH4rwre00P9jM0T8d5egfj9KxrYMjvz3Ckd8e6X9f+aJyalbXUHNmDVPOmkLNmTWUzdbNxURyTUW3SLaYwfTpYTr99OTruIcbBcUX4fv2DX7d0ADPPAMHDoQp2dCIMeXlofieOzc8pno+e7ZazkUmGDOjemU11Surqb+uHoDug90cfSgU4Ef/dJSjDx2lc08nnXs6afphU/97y+rLQhF+5pT+x/L6cvUPFxlH+l9XJJfMQut0bS2sWpV6ve7uUHCnK8z374eWloH+50Ptd/bsoYvzefOgqiqTn1hEMqh0ZikzXz6TmS+fCYD3htbwlkdaaH2klZaNLbQ+2kpXQxcHGw5y8CcDNwUtmVFC9SnVVJ9aTfUp1Uw5dQpVq6oordWNe0SyQUW3yERQWhouxFy4MP16bW2hON+/P0zxz+NfNzaG1vNnnoFNm9Jvc+rUoYvzuXPDjY3UaiaSU1ZsVK+qpnpVNbwhzPM+p+OpjsGF+COt9Bzq4chvjnDkN0cGbaN8YfmgYrx6VTWVyyopmaKSQWQs9C9IZDKprg533TzppPTrdXeH0VeGU6AfPRqmdH3OIfxhMGfO8cV4YoE+Z07oBiMi48KKjKplVVQtq2LOa+cA4SLNrn1dtD7WStvjbbQ91kbb4220P9FO59OddD7dyaGfHRq0nfKF5VQtr6JqRRWVyyupWlFF1fIqyheom4rIcKjoFilEpaUDF2Cm4w6HD6cuymPzDhwIF43u3RumocyYcXxhnuxRF4aKZIWZUV5fTnl9OTMvmdk/33udjr919BfirY+10r6lnY7tHf3F+OH7Dw/aVlF1USjqoyK8clkllSeFqWRGiQpykYiKbhFJzSwUyDNmpO9zDuFmQ4mFeLLHZ54JNx46dAieeCL9NsvLUxfliReGlqofqshYWfFAq/isK2b1z+/r6ePYrmO0b22n48kO2re20/5kO+1b2+lu7Kb10VZaH209bnsltSVUnFTRX4RXnlRJ5cmVVJxUQfn8cqxIBbkUDhXdIpIZlZWwZEmY0unrg6amoYvzAwfChaG7d4cpHTOoqxte63lNjVrPRUaoqKSIqpOrqDq5Cv5+8LLuQ939BXj71nY6nuqg428dHPvbMXqae2jd2ErrxuML8qKKIiqWVAwU4osrqFhcQfkJ5VQsrtAFnTLpqOgWkfFVVBRapmfPhtNOS79u/IWh6R6ffTZcHNrYCI89ln6bVVUDRfiCBQN3Co1dqLpoUSjgVZiLDEvpjFKmvWAa014wbdB8d6e7sbu/CI8V4rHX3Y3dtG9pp31Le9LtFk8tpuKEUIjHHmMFecUJFZTWlarrikwoKrpFJH8N98LQnp7Qej5Ucb5/P7S3w44dYUqlomKgCI8vyhcvhqVLw/Pi4ox+VJHJxswom11G2ewypp077bjlPUd76NgxUIgf230sTLvCY+/R3nCB52PJ7+xbVFU0qCgvX1gepgUDU3Gl/p1K/lDRLSITX0lJaL2eOzf9eu7Q2jpwMejevYPvGBp7bG6G7dvDlExZGZx4Yrgb6dKl4XH58nATpLq6zH8+kUmoZGoJNWfUUHNGzXHL3J3ug9107u7sL8IHPe4KRXm6lnKAkpkllC8op2JhRSjE44vyheFC0uIqFeYyPlR0i0jhMAt9umtqYNmy1Ou1tIQCPDbt2ROmnTtDIb5/P2zdGqZECxbAGWfAc58L558PL3yhhkgUGSEzo6yujLK6MmrOOr4oB+huHlyUd+4No6t07o2mhk56DvbQc7CHtk3JW8sh3CQosYW8fH45ZfPL+h9LZ6ori4yduXuuY8i61atX+4YNG3IdhohMFm1t8NRTA9P27bBlS7jRUFvCf+5VVXDppfCRj8CKFSPelZltdPfVGYp8QlDOlkzwPqfr2a7ji/H453s78e6h6yArM8rmlh1XjJfNGzyvZLqGSCx06XK2WrpFREaqujp0JTn99MHze3vhb3+DRx+FDRvgF7+AzZvhe9+Dn/0sdGeZMiU3MYsUGCsyyueWUz63HFL82ep94WLPzr2dHHv6WH8h3rW/i659XXTu66RrXxc9zT107umkc09n+n2WG+Xz0hfmZfPKKKlVcV6IVHSLiGRKX1+4wLKiAqZOhRNOCEU3hL7kR46o6BbJI1ZklM0po2xO6m4sAL0dvXTtHyjCY4+J83qP9vb3OU+nqKKIsrllA0X5vFCM9xfo0fPSmaUay3wSUdEtIjKU2AWYBw6Em/vEP+7ZA7t2hWnv3lB4Jzr9dLjppqHvACoieam4spjKEyupPLEy7Xq9bb107h9cmHfui2s5bwjPe1uHV5xbSejWUjYvFOjxxfmgAn12GVas4jzf5U3RbWZXAJ8GpgLfBt7j7j1mVgd8HGgGSoEb3L07ek/KZSJ4ZEqfAAAZ/klEQVQix3EPd848ciTc3j52Z8x008GDocDu6Bh6+2bhQsrnPCeMQX7aabBmTRhicJJRzhY5XnF18cBNhNLoae0Jhfj+uKI89nr/QCt6T3NPf5eXtIoIwzOmKs6jLi9lc8ooKivK4CeWkciLotvMFgGXA1cCq4AvAE8DtwA/AN7n7g+Z2XXATcD7o7emWyYik4F7GFu7rS1Mra3hsaUFjh5NPh05knpZb+/o4qisDEMSzpkzMDzhnDkD43cvXhyel5Vl8tPnJeVskbEpmVJCydISqpamL857O3rpOtCVtEDv3D/wvLuxO6x3oAseTb/v0rrSpF1ZEru5aIzzzMuLohs4Afg/7t4DbDSzU4CXmNnvgOXu/lC03k+AJ83sY8CpqZa5e+qxgUQks9yhuzu0BLe3H/8YK5bjC+bY43Cet7eHfWRKZWUYMnD6dJgxY3jT7Nm6ffxgytki46C4spjKJZVULknfraWvq4+uZ5K3lscX513PdNHd1E13U3fKmw7FlNSWJO/KklCcl9TkSymZ//LiSLn7bxNmNQA1wAXAzrj1GqKrfVcD56VZ9mCWQxbJb+7Q1RUK31TFcLLH0a472tbj4aqsDCOGTJkSHmPPp00LFyzGHoeaamoKoiU625SzRfJLUVkRFQsrqFhYkXY97w3DKB7XWr4voTg/ELq29DT3pL35EEDxlOLkXVkSi3ON2JIfRXcSZwOfAtYCBxOWtQDzgfo0y0TyQ08PHDsWCtNjx8b3+XiOwV9SEsajrqw8/jFZsTyS51VVuuV6/lPOFpkArDgMaVg+L/0Nu7zP6T7UnbS1PLEPem9rLx3bO+jYnv66l6KKouRdWRL6oE/mEVvyrug2s5OBRnffZGYOJH6L5UA3kG7ZII2NjaxePTBI59q1a1m7dm1G45Y81tc3UJCOpYgdzfuy3QKcTrpCeKSPQy0rLc3d55yE1q1bx7p162Iv8/q+8srZIpOPFQ3cEZTTUq/n7vQc6UneWp7QzaW3tZdjO49xbOcQI7aUhmEcyxeUU7Gkon+qXFJJxZIKyheWU1SaXxeDDjdn59UdKc2sGPgEcKO7d5rZ+4GXuPvfRcsNOAa8BDg/1TJ3/0P8dnV3szzS1zfQNSFWmCZOI5k/nHW7unL3eYuKQlFaURGm8XpeXh6Kbpnw8vmOlMrZIjJcQ47YEhXuPc096TdUBOULy/uL8IolFVQtq6JqZRWVSytzfgHoRLoj5fXALe4eGxvnl8B1ccsXA13Ahugx1TIZia6uwRe7xU/JCtqxFMy5KoBzUfhWVITCt8D7sMmkppwtIsMy0hFbOp/u5NjOY3Ts7OhvIT+28xidDZ107g4TDyS82aBicQVVK6uoWhEK8SmnTaH6tGqKK3LfTTFvim4z+wAh+Vaa2YnAGuA3wH4zO8XdHwcuAW519y5gg5mlWjY5dXYOHg6tpSV1sZxsNIhUU88Qf1VmWqxLQqxATTalWjaa95SVqfAVyTDlbBHJhkEjtrz4+OV9nX0c2zNQhHfs6KD9yXbat7bT8dRAgX7op4f632MlRvUp1Uw5awo1q2uYdt40qk+pHvcLO/Oi6DazG4GPJsze6u5fM7PXAB82s13AFOCDceukW5Z/OjqgqWnwFCugh/OYrVbi4uKBi9YSp/h+u6MpeBOXlZerABaZ4AomZ4tI3ikqL6JqaVXSFvO+rj46/tZB+5ZQhLc90UbrX1pp39JO619aaf1LKwe+egCA0jmlTL9wOjMunkHd5XWUTM1+SZxXfbqzJav9A7u7Yd++cPvnvXvh6afD47PPHl9gD+eOdumUlISh0WLDo9XUpC6WRzKpJVgkb+Vzn+5sUZ9uEcmkntaeUHRvbOXow0dpfqCZrn0DDZlFFUXMvGwmC69fyNTVU8e0r4nUpzt/9fTAI4/Apk3wxBPw17/Cli3Q0DD8odnKyqCubmCaOTPcoCN+nOFk4w7HnldUqDgWERERGYGSKSXUvrCW2hfWAmHUlfat7Rz+5WEaf9jIkQeP0HhnI413NjLr1bNY9sVllE7P/KhcKrrT6emBH/8YvvlNeOCB0M0jkRnMnw8LFoRbQC9YEKa5cwcX2HV1oVVZRbOIiIhIzpgZ1SurqV5ZzYJ/XsCxPcdouK2BhtsaaPx+I62bWjnz92dSOjOzhbeK7lR27IArrggt2zFLl8Lznw+rVoVp5Uo44QSNUSwiIiIyQVUsquCkm09i/lvn8/ilj9P2eBt/u/5vrPj6iozuR0V3Mn19cPnl8Nhjoah+97vhla8MLdkiIiIiMulULqlk1Y9W8fCyh3nmO8+w9PNLKa7K3FCDKrqTeeihUHDX14eW7mnTch2RiIiIiGRZ1clVVJ9aTdvmNlo3tzLtnMzVgPl1H818sW1beDz/fBXcIiIiIgWkfEE5AN2N3RndroruZFpbw2NtbW7jEBEREZFxFRuzu7elN6PbVdGdTOwOjcW5v2WoiIiIiIwfKwkjzXlPZu9lo6I7md7oLxsV3SIiIiIFxUpVdI+fvr7wqKJbREREpKD0t3R3q+gWEREREckKdS8REREREckyFd3jyTN7kEVERERkggg1N57helBFdzpmuY5ARERERMZTrPzLcBusim4RERERkYhlqdFVRXcy6l4iIiIiUtjU0j2O1L1EREREpLCoe4mIiIiISJZl6ULKkoxuLQfMrBK4GWgBpgM3uPvR3EYlIiLJKGeLSN7LUkeHydDSfTvwO3f/APBj4Etj3qL6dIuIZEvmc7aISDaoe8kAM5sPXAX8JJp1P3CpmS3O0A4yshkRERmHnC0ikgH9o5eo6B5kDdDk7u0A7t4D7ATOz2VQIiKS1BqUs0Uk32XpQsqJ3qe7HjiYMK8FmB8/o7GxkdWrV/e/Xrt2LWvXrk291SlTYN48qKnJXKQiIsO0bt061q1bF3tZl8tYMiw7OVtEJINKppVQNq+M4inFw1p/uDnbMn1l5ngys+uBV7n78+PmPQJ8191vic1bvXq1b9iwIRchioiMiZltdPfVQ6+Z/5SzRWSyS5ezJ3r3kgbC1e/xaoB9OYhFRETSU84WkYI10YvuXwPzoyGoMLNSYBHwYE6jEhGRZJSzRaRgTeii290PAPcBF0azLgLudveG3EUlIiLJKGeLSCGb6BdSArwVuNnMTgfmArraRkQkfylni0hBmtAt3QDuftjdr3X3j7v7P7l7a6a2HXcl6oSj2HNDseeGYp84lLMzR593ctPnnXwmfNGdTRP5BFDsuaHYc0OxCxTesdTnndz0eScfFd0iIiIiIlk2ocfpHi4zawR2j+KtdUBThsMZL4o9NxR7bkzm2E9w91njFUw+KNCcPRr6vJObPu/ElDJnF0TRLSIiIiKSS+peIiIiIiKSZSq6RURERESyTEW3iIiIiEiWqehOwswqzexWM/sPM/uCmU3NdUypmNkVZrbTzA6a2efMrCRu2bfNzKPpx7mMMxUz+1hcjJuieXl//M3sobi4Y9PLo2XHfaZ8YGYvjeJeHDcv5bHOp+8hRewpz/1oeV6c/8lij+YnPU/y6bhPFIV2zIY69ycrMyszs01mtibXsYwHMzvHzK43s8vNrC7X8WSLma00s8+b2b+Y2e1mdkauY8qWgviHOgq3A+vd/U4zuwT4EvC6HMd0HDNbBFwOXAmsAr4APA3cEi17FnhptPqTOQkyDTOrJlytHIsxNlpBXh9/M1sJ/BG4HuiOZn8F+FWaz5RTZjYHmAI8L2FRumOdF99DstjTnftxy3N+/qc67kOcJ3lx3CeYgjlmQ537k9x7gcW5DmI8mNmbgZPc/cZcxzIOvgVc5u4N0fn9c2BljmPKCo1eksDM5gM7genu3h61IBwBVrn7rpwGl8DMXgT80d17otc3E+J8hZn9F/B74J7Y8nxjZu+Knn7R3Y9F8/L++JvZInffE/f6BOBmd39Nss+UL8ysCOgFlrj7rnTHGuhKtSwX30OS2FOe+9HrvDn/E2OP5iU9TybC+Z9vCu2YDXXuT1Zmdi6wAvgwcLW7P5DbiLLHzM4HPgS81AugSDOzNuAsd99qZrOATe4+P9dxZYO6lxxvDdDk7u0AUWLbCZyfy6CScfffJhQUDcAeM6sC/gH4fvT6wpwEmIaZGfBa4NPAM2YWa5VaQ54f//iCO/JK4J40nykvuHtfwqw1pD7W6ZaNu8TYU537APl2/ifGPsR5soY8Ou4TxBoK6JilO/cnq+iXoSvd/Wu5jmWcfAbYAtxqZuvN7LxcB5RldwFfMbMa4B+Bf8pxPFmjovt49cDBhHktwET4q+tsQstZu7ufTPgs/wP8zMxOzW1og3lwDuEn9luB75jZy5iYx//lwP+k+Uz5Kt2xnmjfw9nAFwHy/fwf4jyZaMc9HxT6Mes/9yexG4BP5DqI8WBmy4HnAl9293cA/wusj7qqTVZvA44BfwZa3f3uHMeTNSq6j+dAR8K8cgb67uYlMzsZaHT3/guy3H2/u18L/AB4R86CS8PdD0d91v4TeBcT7PhHF7ccc/ejsXlJPlO+SnesJ8z3kOzch/w//1OcJxPmuOeRgj1mqc79ycTMLgY2uPuzuY5lnKwCDrn75uj1bYABl+UupKyrBO4Avgt81swuynE8WaOi+3gNwPSEeTXAvhzEMixmVgy8BXhfilX+G1g0fhGNSizGiXb8LwXuTbEs3497umM9Ib6HYZz7kP/fQ3x8E+K455mCPGbDPPcng+uBr5lZk5k1AQuBH5vZv+Y4rmwpIW6QC3fvALYDM3MWUfZ9G7jb3W8i/KJxR9SlaNJR0X28XwPzzawSwMxKCf8hPpjTqNK7HrjF3TtTLO8DHhnHeEYjFuNEO/7piu58P+7pjvVE+R6GOvch/7+H+PgmynHPJ4V6zIZz7k8GVwFnxE37gGuYvF1qNgPTErqT9ABP5CierIp+LT7d3ZujWR8DjjJJRy9R0Z3A3Q8A9wGxi68uIvwF1pC7qFIzsw8AjwKVZnaimf0fMzvFzC6NlpcBbwJuzmWcicxsvpldZWZF0YVl7wE+OJGOf/SXeEXsZ89UnymnQcaJYoLwU2Xacz3fvofE2KN5yc79k81sRT6d/4mxpztP8u24TwSFeMxSnfu5jisb3L3R3ffGJsJIQI3xXfomE3ffCvyMMCQkZlYLlBGuT5mMDgHHzKw+bt5BYFuO4skqDRmYhJlNJ/wnvQuYC7zf3VtzGlQSZnYj8NGE2VsJ/VfvIFzB/yfgk+6eVz+1mtkqQmI5CvwO+GyUbCbS8b8SWODun41ep/xMuRZ3VfjtwE3Abe7emO5Y58v3kCx2ws/qx5377r4yGq0kL87/FLHPJs15ki/HfSIppGOWKu+7+6RsGUxkZruY/EMG1gGfAzYQfrX5srtPypZuADM7nXAx5UZgDvAbd5+Uv1Sp6BYRERERyTJ1LxERERERyTIV3SIiIiIiWaaiW0REREQky1R0i4iIiIhkmYpuEREREZEsU9EtIiIiIpJlKrpFRERERLJMRbeITBhmNsvM7jWz3Wb2vlzHIyIi6SlvD9DNcURkwjCzi4FfA1MJd19d5u6HchuViIikorw9QEW3iExIZvYb4AJ378l1LCIiMrRCz9vqXiIiE46ZLQR+UaiJW0RkolHeVtEt48zM/tXMPp/rOGIyEY+ZXWRmfzSzNRkKa1Ixs1syvD0DXg18IpPbFZHBlK/HzswWm9mXzexLGdpesZlda2a7h7HuK8zsodF+1kzmbuXtQEX3JGJmK83sh2b2FTPbZGZuZv+Y67gS/BS4I9dBxMlEPE8Dp2cglrxkZi83syYz229mLzCzhdF55mZ2m5nVRus9z8weN7NPm1llNG8OcCDDIV0FfNXde6Lti0w4ytejMhHz9bPRY3mGtlcENAOLhrHuBmBlqoVmtsrM7ozOvRvNbG7cskznbuVtVHRPNvcA/+7u1wDPBb4ZW2Bmp5nZ+bkIKn7f7v64u/8uF3EkM5J4Uh1Dd38SaMpkXJn8vsa6LXf/KXA70OPuf3T3p4FrgS7gIXdvjtZ7mJDk3+vuHdHb/wG4b0wfII6ZfQj4N+B+M9sKLM/UtkXGmfL1CE3EfO3u7cC+TO3L3buBjcNc9xlCgZ5q+V+BddHLW9w9vsjOWO5W3h6gonuSMLMZwDKgDcDd+4APA25m0wgJ3XIQV872nUnD+Bx947ivXGzrTmCBmZ0F4O4Hgf8F/j5uX7XAvujci1nh7lvHuO9+7v5Rd1/q7qvdfYW7/yZT2xYZL8rX2ZWH+TrTI1aMJP6h1u2NHhP7WWcsdytvD1DRPUlEw+88Caw3s3OiebuBhwn9qJYAa83szQBmdpOZfSj6aam/35aZrTazz5jZu8ys3cyeNLN3m9nlZvZnM3ulmf3IzI6a2Vvi3pd0ewn7vtbMvhHfty36ifVWM/tPM/u1mT0vmp92f3HvNzN7i5n91cxeamabzewZM7t6GPtYGh/PEPs87hgmMT36ubjFzK5NOKb/Fh2bH0ZJerjH7Lh9mdmVZvawmV1tZtvN7BNptpfsu08aTzpRi8hfgVfFzZ4FvMKiriTApcC9cXFWAu1xr1eY2SfM7BYze9BCN5XXmNkGM7vMzO4ws2Yze6eZnRV9xoNm9vKh4hOZSJSvCydfJ5PqM6bbj5ktMbOvWWg1/tRw9pNkvz81s++Y2XlDrKfcnS3urmmSTMAqYDfhL9uvAbPjlu0C1kTPTwF2RM9nEf4Kn0n4S30PsDpa9nngD9HzMuAQcBthrM3XAQ3ptpe4b0Kfti8B34jmVxAKudro9ZVAI1CTbn9JPvfcaJ8fiNb9GHAMOGGIfSTGk3af8ccwSQy7gE8DVdH79kTzpwLfi1vv94SEOaxjlmJfMwg/GX4WeB5wUbrtJXz3SeMZ5vl1I/C36PlJwH8CncAVceeLxa1/KfCi6HkV8GegKnr9m+g4TAWOADcBlcBrotevJjQKvIfoHNSkaTJNKF8XRL6Olv/bMI9juu/6z8DKuPf4MM+z2Pf5AuAjSZavifZTEjdPuTtLk1q6JxEPrZGnEpLQG4G/mlmyC0a2A28ws1LgRdG8KYR/5AuBWH/cLYR/TLh7F9AC3OXuRwn/COcOsb3E+DoZfGHGK4AOj/oEAz8k/GO9Yoj9JW43ts3vRut+jNDf+KVD7GNQPCPZZwo/8dB/78/AvLjPODNqiXoX4ZjCMI9Zis97iJDcfuzuD7v7/SPYXqp4huNO4EQLXUz+P+DLwM+BV5vZFKDFo6wcOQ/4Q/T8UsJ/bLHWk1cCH42OczPwKw/9wDcAU939+x5+cn+UkX0HIhOC8nVh5OskUn7GNPu5CJjm7rF4Noxwn2cDr3f3fx/m+srdWVKS6wAkM8ysGFjo7ruAfzazrxKu9F4HPD9+XXfvNLN64Abg/8U24e7PmtkG4ELCX+JLCIVW/1vjnvcQdU9Ktb0UocZvYylQGhdXn5ntAhak299Q3P2Yme0AaoG6Eexj1PtM0MPAv61FwE53/2ziSsM5Zmb2BkJxG7Pcw8/QHh/rCL6DlPEMxd23m9mjhC4m09x9h5ndSWh9ehXwP3Fxx86NWH/BEwgtU7FtHYzfdNzzXgbrQ93gZJJRvh70+QohX8dLeRzTfDenEf7AGK0LgKVmNsXdW1Os0xd9BuXuLNIBmTxKietv6+6bCD/frUhc0cwuAN7m7h/zMBJFvLcDZ1sYumoz8MmhdjzE9tLZQ2g5jR9KyQh9HceqDNiW5X0Mx34G93vGzJ4/gmN2L3BG3JT0KvgRbC9pPMP/OHwPeDOwIy6+YuCdhJ9iY55P6J8a0wC80Myq4/abtl+hyCSmfD1YIeXrlJ8xzX6ORu8pY3RuIfxK8LkUy4/4wAXwyt1ZpKJ7cvknM5sf97oeWB897yJcOLKCMDzVNDMrN7OLouWzzayO0Lfwq4Sfhv5A9HNlpIiEv+7NzIbYXuK+LW4b9wCtwGujbdVG+4u1mKbaXyrzonXqo/euH8Y+4uMZap/xnyNRslYPI7ReVQP3mdnFFi7YGeo76N8XMM/dt8ZN3XG7KI57nm578XEnjcfMysxsXdR1JJ3vE37WvgvA3VuAnwG/98GjlvwdoetJTOx432FhrO93ALGxWhO/g8TveUKPpCCSgvI1BZOvh3sck+6HkEvLCH+YQfhVAzObP4zcbYRW5zcArzKzV0fvfY+Z3Uho5f/fuPWVu7MpGx3FNY3/RLg4w4HDwLeBbxB+npoeLf8IsBe4jPBz0XZCy8IrCC0k3yAksHsJF3V0RdtrJ/wjfBnQTbhYZyHw/mj5/0m3vYR9vxP4LaGf3Kpo2WpCC+mnCX0bT4nmp9xfis/v0fvfSejucGrcslT7WBwfz1D7jD+GCfu+jHBB4W3RsYi9783R8vOBx6Pv5tPRcR7uMbssyWe9MtrfHUB9NC/ddzpoWyniWUC4icN3h3Gu/XfC69cAL0mYd0uS970YeAI4CNwYzXs54Vy7FZhPNGwacB3hP+UvR5/173P9b0yTpkxNKF8XUr5eBPyO0Fp/xhCfMd13/QrCL4x/IlzUvg14W/TZk+Zuwq8p3YQW7hnArwgF/zXAWwjjlX+a0LUl9h7l7ixOFh08EcxsNnCdu98UvS4itL68xd1vzGlwQzAzB5Z46CMpo2Rm73X3UQ1HJSLjR/la4il3TwzqXiLx3kHoajAjbt4FxF0gl4/ifs7ST1mjZGbFZvZ64O5cxyIiw6J8LcrdE4yKbon3RULfry1m9jTwS2Cvu/8xt2GlFl1Y8s7o5dUWhq6TkZsC3OPuO4ZcU0TygfK1gHL3hKLuJSIiIiIiWaaWbhERERGRLFPRLSIiIiKSZSq6RURERESyTEW3iIiIiEiWqegWEREREckyFd0iIiIiIlmmoltEREREJMtUdIuIiIiIZJmKbhERERGRLPv/Afx/cLZn9wk7AAAAAElFTkSuQmCC\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(vehicle1.v_kmsc, vehicle1.h_kmc, 'k-', linewidth=2.0)\n", "\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(vehicle1.acc_net_g, vehicle1.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(vehicle1.q_stag_total, vehicle1.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(vehicle1.heatload/1.0E3, vehicle1.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/dragnofly-titan.png',bbox_inches='tight')\n", "plt.savefig('../plots/dragnofly-titan.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/dragnofly-titan.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Vehicle geometry and ballistic coefficient data is not publicly available in the literature. Hence, results are only approximate." ] } ], "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 }