{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 43 - Ice Giant Pre-Decadal - Uranus Probe (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 described in Sec. A.5.3 of the Ice Giants Pre-Decadal Mission Study. https://www.lpi.usra.edu/icegiants/mission_study/Full-Report.pdf" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"URANUS\") \n", "planet.h_skip = 1015.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", "vehicle1=Vehicle('igpd-uranus', 325.0, 205.0, 0.0, np.pi*1.2**2.0*0.25, 0.0, 0.20, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle1.setInitialState(1015.0,0.0,0.0,23.10,0.0,-30.0,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": "iVBORw0KGgoAAAANSUhEUgAAAuIAAAHtCAYAAABcXRt3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzde3xcdZ3/8dcn16bphZakLSmUcr9ZQYkoilC5r+4CIqggKu5CXFAWFFy8oK4Cv11B+amwIEHR9Q4s/kS6gIJg5bZCuRQEuUkv2EJJekvaJk3TfH5/nDPJJGQmmczlnDPzfj4e85iZc2bO+cwZ+pkP33zO95i7IyIiIiIipVUVdQAiIiIiIpVIhbiIiIiISARUiIuIiIiIRECFuIiIiIhIBFSIi4iIiIhEQIW4iIiIiEgEaqIOoBSampp8/vz54359R0cHzc3NxQsoZvR5y1clfVYoz8/72GOPdbp7eX2oMeSas1PK8fsvNh2z3OmY5a6SjlmuObsiCvH58+ezZMmScb++tbU1p9cnnT5v+aqkzwrl+XnNbEXUMZRarjk7pRy//2LTMcudjlnuKumY5Zqz1ZoiIiIiIhIBFeIiIiIiIhFQIT6Ktra2qEMoKX3e8lVJnxUq7/PKcPr+c6djljsds9zpmGVm7h51DEXX2trqldKbJCLlxcwec/fWqOMoJeVsEUmqXHO2RsRFRERERCJQ8llTzOwY4DLgQ+6+PFzWAFwBdAMzgIvdvSufdSIiURsYGKC/v5+6urqoQ5mwJORsdxgYCB5XV+e7NRGR0inpiLiZzQamAIeMWHUt8IC7fxG4Dbi+AOtERCLz9NNPU11dTWtrcrtKkpKzH3oIamrg8MMLsTURkdIpaSHu7msIku8gM2sBTgduDxfdA5xgZvMnuq6oH0JEZBxSo+B9fX0RRzJxytkiIsVV8tYUdx8ws/RFC4FOd98Sru83s2XAEcC2Ca5bXqKPIyIyqnIoxEE5W0SkmOJwZc25wNoRy7qBFqB/guuG6ejoGPbn4ba2Nk2lIyJFlSrEt27dmvN729vbaW9vTz1tKlxUBRHbnF0Bk4CJSAzlk7PjUIg70DNiWT3B6Al5rBvU3NxcMZdWFZF4yGdEPL3wNLPOggaWv9jl7OED9iIipZVPzo5DIb6K4Oz5dFOB1QQJfyLrREQiVV9fDyS/NWUUytkiIgUSh3nE7wNawmmtMLNaYB6wOI91IiKRKpce8VHENmerNUVEkqbkhbgNnfVjAO7+GrAIOCpcfjRwq7uvmui6EnwMEZGsamtrgaBHPMlXME5CzlZriogkVUlbU8xsKnBG+PRMM7vG3TuAc4ArzOxAYA6QflbORNeJiESmurqa6upqtm/fzvbt26mpiUMnYG6Us0VEiqukvwzu3g1cF97Sl68Hzs7wngmtExGJWl1dHT09PfT19SWyEE9azk7wHx5EpELFoUdcRKQslfEJm7Gi1hQRSSoV4iIiRZLPXOIiIlL+VIiLiBRJGc+cEktqTRGRpFEhLiJSJCrES0OtKSKSVCrERUSKRD3iIiKSjQpxEZEiUY94aak1RUSSRoW4iEiRqDWlNNSaIiJJpUJcRKRIVIiLiEg2KsRFRIpEPeKlpdYUEUkaFeIiIkWiEfHSUGuKiCSVCnERkSLRyZoiIpKNCnERkSLRiHhpqTVFRJJGhbiISJGoEC8NtaaISFKpEBcRKRKdrCkiItmoEBcRKRL1iJeWWlNEJGlUiIuIFIlaU0pDrSkiklQqxEVEikSFuIiIZKNCXESkSNQjXlpqTRGRpFEhLiJSJOoRLw21pohIUqkQFxEpErWmiIhINirERUSKRIV4aak1RUSSpibqAFLMbD/g08BLwF5Au7s/aWYNwBVANzADuNjdu8L3ZFwnIhK1cu4Rj1POVmuKiCRVbApx4CfAie6+yszmAb8F9gOuBe5y95vM7HjgeuC08D3Z1omIRKrMe8SVs0VE8hSn1pT9gKnh4x5gupm1AKcDt4fL7wFOMLP52daVLGIRkSzKvDUldjlbrSkikjRxKsT/G/i+mU0FzgDOAxYCne6+BcDd+4FlwBFjrBMRiVyZF+KxydlqTRGRpIpTIX4u0As8Cmxy91uBucDaEa/rBlrGWCciErkyL8SVs0VE8hSnHvEG4BcEyfrbZrYMcII/eaarB7aFj7OtG9TR0UFra+vg87a2Ntra2goUtojI6CZ6smZ7ezvt7e2pp02FjapgYpez1ZoiIlHIJ2fHqRD/KfBhd99gZk6Q4C8gOLM+3VRgNUHCz7RumObmZpYsWVL4iEVEspjoyZrphaeZdRY8sMKITc5Wa4qIRCmfnB2L1hQzawIOdPcN4aLLgC5gBdASTnmFmdUC84DFwH1Z1omIRK5cW1OUs0VECiMWhTiwDug1s7lpy9YCTwGLgKPCZUcDt7r7Knd/LdO6EsUsIpJVuRbixDRnqzVFRJImFq0p7j5gZicBXzGzx4DZwOfcvcvMzgGuMLMDgTlAeqNgtnUiIpEq1wv6xC1nqzVFRJIqFoU4gLsvBT45yvL1wNkZ3pNxnYhI1Mr5gj7K2SIi+YtLa4qISNkp49aUWFJriogkjQpxEZEiUSFeGmpNEZGkUiEuIlIkKsRFRCQbFeIiIkWSOlmzHHvERUQkfyrERUSKRCPipaUecRFJGhXiIiJFokK8NNQjLiJJpUJcRKRIVIiLiEg2KsRFRIqkurqa6upqBgYG6O/vjzqcsqfWFBFJGhXiIiJFpFHx4lNriogklQpxEZEiUiEuIiKZqBAXESkiFeKlo9YUEUkaFeIiIkWkucSLT60pIpJUKsRFRIpII+IiIpKJCnERkSJSIV46ak0RkaQpSCFuZu8txHZERMpNHAvxcsvZak0RkaSqyeXFZnY8cAkwh6EivgqYDTQUNjQRkeSLshBXzhYRibecCnHgJ8C/A08AA+EyA04vZFAiIuUi4pM1KypnqzVFRJIm10L8ZXe/auRCM3uuQPGIiJSViFtTKiJnqzVFRJIq1x7xc8zslFGWn1CIYEREyk3EhbhytohIjOVaiL8L+L6ZbU+7DQDXFSE2EZHEi7gQr6icrdYUEUmaXFtTvg60AUuA7eEyAz5ayKBERMpFxD3iFZGz1ZoiIkmVayF+J3Cbuw/7RTGzGwoXkohI+Yh4RFw5W0QkxnItxH8JfMHM7k1bZsDJwPmFCMjM3gEcBrwEPABsBq4AuoEZwMXu3hW+tiHTOhGROIi4EK+onK3WFBFJmlwL8S8DewFnAqmUZwRz1Oad1M3sE8Ae7n5J2rIfAne5+03hnLjXA6eFq6/Nsk5EJHIRF+IVkbPVmiIiSZVrIX4lcKe7d6cvLMRV2szsCOAjwDFpy1oI5rv9VLjoHuBWM5sP9GVa5+7L841HRKQQIu4RV84WEYmxXAvxV0Ym9NCcAsRyFfAQcLWZ7QlcCuwKdLr7FgB37zezZcARwLYs65YXIB4RkbxFPCJeUTlbrSkikjS5Tl947sgFZnYq8O18gjCzfYC3ADe4+6eB3wN3AbsAa0e8vBtoAeZmWTdMR0cHra2tg7f29vZ8whURGbeJFOLt7e2D+QpoymP3FZGz1ZoiIlHKJ2fnOiI+08y+7O6Xmtlk4LvAqcCjOW5npAOAde7+VPj8GuCrBD2NPSNeW08wssIY6wY1NzezZMmSPEMUEcndRArxtrY22traADCzzjx2r5wtIlJk+eTsnApxd3+fmb3DzK4GjgNeAw4ENuaynQxxDMbi7j1m9iJQTXBmfbqpwGqChJ9pnYhILKQK8Sh6xCstZ6s1RUSSJmshbmaHZ1i1AngR+CZBT+BJwGfyiOMpYLqZzXb3NeGyfoIE3WJmDWGirwXmAYsJLk5xQ4Z1IiKxkDpZsxQ94pWas9WaIiJJNdaI+BXAfsAGYGCU9TcSTIW1E3kkdXd/zszuBE4B/tPMdgDqgJ8BfwccBSwCjgZudfdVAGa2KNM6EZE4KPHJmsrZIiIJMlYh/nVgibu/nu1FZnZKAWL5GPAdM6sjGCU5LTyr/hzgCjM7kOBM/7a092RbJyISuRIX4hWds9WaIiJJk7UQd/c7xrMRd//vfANx906COWlHLl8PnJ3hPRnXiYjEQSkL8UrN2WpNEZGkynX6QhERyUHEF/QREZEYUyEuIlJEEV/Qp6KoNUVEkkaFuIhIEakQLz61pohIUuVdiJvZW8zsp2Z2oZnNzz8kEZHyEbdCXDlbRCQ+CjEi/s/A0e7+LeCdZtZQgG2KiJSFGPaIl23OVmuKiCRNIQrxTwMHh49/CdQWYJsiImUhbiPilGHOVmuKiCRVzoW4mf2dmX0kfHwIsFfqggzuPuDuXQWOUUQksaIuxJWzRUTiK6dC3MwuB24B3g/g7o8AHzOz04sQm4hI4kVZiFdazlZriogkTa4j4u8EdgbuT1t2J/AfBYtIRKSMpArxiHrEKyJnqzVFRJIq10L8bnffAKSPOxTiUskiImUpdbJmRK0pytkiIjGW9RL3o3jZzL4IzDWzDwOnEvzJ87yCRyYiUgYi7hGvqJyt1hQRSZqcRsTd/ZfAvUA18DFgA3Csu/9nEWITEUm8KAvxSsnZak0RkaTKdUQcd/9f4H9Tz82s2syOcvffFzQyEZEyEHGPuHK2iEiMZS3Ezew3wJQxttEEdAJK6iIiI5SyR7zSc7ZaU0QkacYaEX8qfM3zwABwNvAEsCTtNbsBM4sSnYhIwpW4NaUic3ZteEmi+FwzSURkfMYqxK8ELDzrHjPbw92/MvJFZnZbMYITEUm66upqqqqqGBgYYPv27VRXVxdzdxWZsxsbg/vNm6ONQ0QkV1lP1nT3jamEHpptZsN+RcxsIXBIEWITESkLpeoTr9ScnV6Iqz1FRJIk15M1/xtYamYPAN3AfsBxwBtGXEREJFBXV0dvby99fX1Mnjy5lLuuiJxdUwP19bB1K/T2QkND1BGJiIxPrtMX3g0cS3Cizz4EU2Gd4u7/XoTYRETKQlQX9amknK32FBFJoolMX7gauCT13MyqzOxgd3+soJGJiJSJiOcSr4ic3dgI69bBpk3Q1BR1NCIi45NTIW5mN46yuBnYDpxUkIhERMpMVIV4JeXsKeGkjRoRF5EkyXVEfC+CuWfTT4fZGbi/EMGYWR3wKHC+u//BzBqAKwh6G2cAF7t7V/jajOtEROIkwov6FDVnQ3zytlpTRCSJci3Ez3L359MXmNm+wJkFiudzwPy059cCd7n7TWZ2PHA9cNo41omIxEZUPeIUP2dDTPK2CnERSaJcT9Z8fpTFW4Gz8g3EzN4JvAqsD5+3AKcDt4cvuQc4wczmZ1uXbxwiIoUW1Yh4MXM2xCtvpwrxTZsKsTURkdLItUd8GcP/xFkNzAZ+kk8QZtZIcCb/Z80sNa3WQqDT3bcAuHt/uP8jgG1Z1i3PJxYRkWIxs1Lvryg5O9x2rPK2RsRFJIlybU25G/g5Q4l9AFjj7i/kGcfFwH+MWDYXWDtiWTfQAvRnWfcGHR0dtLa2Dj5va2ujra0tn3hFRMZtYGAAgKqq8f0Rsr29nfb29tTTfOYAKVbOhiLm7YnkbJ2sKSJRySdn51qI/9Td/5i+wMz2MbM93f2lHLeVev9xwBJ3f33EKgd6RiyrJxhVYYx1wzQ3N7NkyZKJhCcikrdcC/H0wtPMOvPYdcFzdriNoubtieTsadOC+7UjS30RkSLLJ2fn1CMOHD/KsjXALTluJ91FwI1m1hkGvwtwG9BIcFZ9uqnAamBVlnUiIrGSayFeQMXI2RDDvL3ffsH9n/9ciK2JiJTGuH4VzOzTZrYBuNjMtqffCP7UuCaPGE4HDkq7rSY4keiHQEs43RVmVgvMAxYD92VZJyISK6UuxIucsyGGefugg4L7J58sxNZEREpjXK0p7n6Nmf0R+Djw3RGr+9z91YkG4O4d6c/DH4oOd19pZouAo4BFwNHAre6+KnxdxnUiInGSKsRLdbJmMXN2uP3Y5e0FC6CmBp55BjZsgB12KMRWRUSKa9w94u7+lJl93t3f0M9nZvPcfWVhQwPgHOAKMzsQmAO0jXOdiEhs9Pf3A1BdXV2yfUaUsyGivD15MrzznfDHP8K998LJJxdqyyIixZO1EDezi4A73P3ZcNG/mtnIXxID3kMwbVXe3H1+2uP1wNkZXpdxnYhInHR3dwMwLXVGYZFEkbMhPnn72GODQvx3v1MhLiLJMNaI+G5A+h/43kpwMs5rDE2HZQQjGyIiMooNGzYAMH369GLvqqJz9rHHwiWXwG9/C+5Q4mnbRURylrUQd/dPjVh0OfBndx92nWYz27/QgYmIlINt27bR09NDdXU1jamrzhRJpefst74VZs6E5cvhr3+FPfeMOiIRkexyPYV/wSgJfR9gj8KFJCJSPjZu3AgEbSmlvrImFZazq6vh6KODx7/+dbSxiIiMx7hO1jSzdxP8OfOY8JLE6WYSnJV/e4FjExFJvFdfDSYomTOndN0glZyzTz8dbr4ZbrgBLrxQ7SkiEm/jnTVlA/Bjgv7Dd41Y1wd8v5BBiYiUi5Urg8lJ5s2bV8rdVmzOft/7YO5ceOEFWLwYFi6MOiIRkczGO4/402a2EDjV3cs2gYuIFNqKFSuA0hbilZyza2rgrLPga1+DK69UIS4i8TbuHnF335gpoZvZYYULSUSkfCxbFnSGlHhEvKJz9qc+BVOmwB13wMMPRx2NiEhmY80j/htgyhjbMOAAYFahghIRKRdLly4FYMGCBUXfl3J2oLkZzj8fLr8cvvQl+P3v1SsuIvE0VmvK8wS9hqsYmoN2JAP+vpBBiYiUA3fnySefBOCggw4qxS6Vs0MXXgjXXgv33Qc/+xmccUbUEYmIvNFYhfg3gG5335rtRWZWlmffi4jk45VXXqGjo4MZM2aUqjVFOTs0YwZ861vwj/8IF1wQXOxnVtn+DUBEkiprj7i7d46V0EPvL1A8IiJl45577gHg3e9+d0nmEFfOHu7MM4N5xdeuhdNOg/7+qCMSERku1wv6DGNm88zsG8B3ChSPiEjZuPvuuwE49thjI44kUGk52wz+679gzhy49174/OejjkhEZLicC3Ezqzaz95vZHcBfgY8Q9CSKiEiot7eXO++8E4i2EK/0nN3SArfcEkxr+K1vwdVXRx2RiMiQcRfiZrabmV0OvAL8EtgT+AdgF+DU4oQnIpJMixYtYuPGjbz1rW9lr732Kvn+lbOHHHYYfO97weN/+ZehxyIiURuzEDezU83sd8CLwAeAq4Cdgf9x97s88FCR4xQRSZQf/OAHAHz0ox8t6X6Vs0f3T/80NBp+zjnwjW+AZ5pXRkSkRMZzZc0egksiPwSc5e4vAJiZUpiIyCieeuop7rrrLhoaGjij9PPmKWdn8OlPB8X3+ecH/eIvvxwU53V1UUcmIpVqzELc3RcBi8xsHnC2mc0HbgNqU68xM3PX2IKICMDll18OwFlnnUVTU1NJ962cnd155wV94x/5CLS3w9KlQQ/5LrtEHZmIVKJcLnG/0t2/DHwiXLSzmV0aJvmzixCbiEjiPPjgg9x8881MmjSJiy66KLI4lLMz+8AH4I9/DIrvP/0JDjoIfvITtaqISOnlPGuKu/e7+3+7+/uB/wI+A3y74JGJiCTM1q1b+dSnPgXARRddVKqL+GSlnD26Qw6Bxx+H44+HdevgYx+D974XVqyIOjIRqSR5zSPu7i+5+/nAPxcoHhGRxPryl7/M0qVL2W233fh8DCetVs4erqkJ7rgDbrwRdtgB7roL9t8fvvpV2LQp6uhEpBLkVYinuPuPC7EdEZGk+s1vfsM3v/lNqqur+dnPfkZjY2PUIWWknD3EDD7xCfjLX+CDH4QtW+DrX4c994QrroANFTPjuohEoSCFeCGY2clmtszM1prZd8ysJlzeZGbXm9k3zOwqM0s/4SjjOhGRUlmyZAmnnXYa7s6ll17KoYceGnVIRVduOXvOHLjpJrj/fnj722HNGrj44qCP/IIL4Omno45QRMpRLArx8Oz+k4BTCPoXzwIuCFffAtzo7hcDLwBfT3trtnUiIkX3yCOPcOyxx7JlyxbOPPPMWLakFFo55+zDDoOHH4Y774SjjgpaVL7zHXjzm6G1Fa65BlatijpKESkXFocZrMzs3cDD7t4fPr8COAC4FPiVu7eEy+cCzwOzgQWZ1rn75vTtt7a2+pIlS0r1cUSkQtx555188IMfZNOmTZx44oncfPPN1BV4Umoze8zdWwu60TxVUs5+4gm44Qb4xS+Gt6kcdFBQqL/jHXDooTB3bnQxikh85Jqzx3NBn6Jz9/tHLFoFTAWOBJalvW6VmQG0Au/Ksm5xkUMWkQrW19fHZZddxmWXXYa7c/rpp/OjH/2I2trYdFoUVSXl7Le8Ba69Fq66Cm67LSjI77kHnnwyuKXMmQMLFsCb3hTcFiyA/faDKVOii11E4i8Whfgo3gZcCbQBa0es6wZagLlZ1g3T0dFBa+vQ/5y0tbXR1tZWyHhFpEIsXryYc889l2effRYz47LLLuMLX/gCVVWF6/Rrb2+nvb099bS0VwSamLLP2ZMmwYc+FNx6e4N5yB98MGhj+dOf4LXXgtvddw9/35w5wYmf6bfdd4edd4ZZs6C6OprPIyKFk0/Ojl0hbmZ7Ah3uvjS8JHPPiJfUA9uAbOuGaW5uJi5/5hSR5HF37r//fi677DLuDiutPffckxtuuIGFCxcWfH/phaeZdRZ8BwVUiTl70iQ49tjgBjAwAMuXw5//PHR7+ml48cWhAv2BB964nZqa4Cqfc+cGhXnqttNO0Nwc3GbNCqZZrJA/togkUj45O1aFuJlVA58EUmc7rQL2TltvQCOweox1IiJ5W7NmDTfffDPXX389zzzzDABTp07lwgsv5F//9V9paGiIOMJoKWcHqqqCUe7dd4cTThhavn17cGLnSy8FRflLLwW3ZcuC5Z2dsHJlcBvLDjsMFefpBfqMGTB9erB+tPtJk4IpGkUknmJViAMXAd90963h87sZfuGJ+UAfsCS8z7RORCRnAwMDPPXUU/z+97/ntttu44EHHiB1Qvvs2bNpa2vjggsuYObMmRFHGhvK2VlUV8O8ecHtyCPfuL63F1avhr/9bfjttdego2Po1tkZnCi6YUNQ0Oeiri57oZ66z7Rs+nS1z4gUU2wKcTP7IkFCbjCz3YGFwB+BV83sTe7+Z+B44Gp37wOWmFmmdSIiY+rt7eXJJ5/kkUceYfHixfzhD39g3bp1g+vr6uo45phjOPPMMznxxBMr5mTM8VDOzt+kSUMj6dkMDMD69UOF+euvDxXoGzcGBXrqfuTjvr6h903UlCm5Fe8jHzc2alReJJNYFOJmdgnBtFfpnnP3G83sQ8BXzGw5MAX4Utprsq0TERnU29vLs88+y9KlS3n00Ud55JFHWLp0Kf39/cNeN2/ePI488kiOO+443vve9zJt2rSIIo4v5ezSqqqCHXcMbvvum9t7e3uHF+ej3Wdb1tUVzKW+adPE50+vqQnaaJqagraa1ONsyyq860sqSCzmES+2OM1JKyLFNTAwwIoVK3jqqad4+umnB+9ffPFFtm/fPuy1Zsb+++/PIYccwqGHHsqRRx7J7rvvjsVo+C6O84gXm3J2fAwMQHf32AV7tmU9I0/RHYfGxvEV7KnlO+6oFhqJh0TOIy4ikqsNGzbwwgsv8Pzzzw/eXnjhBV588UV6Rvnlr6qqYt9992XBggW0trZyyCGHcPDBBzN16tQIohdJhqqqoTaTierthbVrh9pp0m+Zlm3eHNxWrBjfPsyCojx1Euu0aUHM06YN3Robg9vkyWPfqwtNSkWFuIjE0oYNG1ixYgXLly8f9X7t2pFTUg+ZM2cOCxYsYMGCBbz5zW9mwYIF7LfffhU/y4lIFCZNCqZoHO/VR92DUfixivX0x+vW5d8Ln662dnwF+0TvNZuNpKgQF5GS6uvr49VXX2X16tWDt1WrVr3heVdXV9btNDQ0sPfee7PPPvsM3qceT89n+E5EImU2NIo91omsKdu2BUX5668Ho+9dXcNvGzcGI+xbtozvftu2oTabYn3GYhb6kyerVScpVIiLSF7cna6uLl5//XU6Ojro6OgY9fHrr7/O6tWr6RjnkNXkyZOZP38+u+6666j3s2bNKujVLEUkuWprgwsh7bRT/ttyD2abyaVwz/U+tf3Nm/OPN5P6+okX8w0NweOGhqHbyOcNDRrZLwQV4iICBAV1T08P69evZ926dVnv165dO1hod3R00Nc3/hnoqqurmTNnDi0tLcNuc+fOHfZ85syZsTppUkQqg1lQxNbXQ7EuGdDfHxTloxXqhSj2t2yBrVuD2/r1xfkMEByrSZMyF+ujFe/jLfBHLivXwl+FuEiZ6Ovro6uri40bN9LV1TXscaZl69evH1Zgb926dewdjaKxsZFZs2bR3Nw8eD/a45aWFmbNmkW1/mYqIhWspmao/aYY3IOTZCdSyG/eHMx0k7pt2TL8efqyvr6hZaVSX5+5SB9tefq6xsbMc95Pnx5cAKvUVIiLRMTd2bJlC5s2baK7u5tNmzYN3tKfpx53d3ePWlCnHk+0iE5XX1/PjBkzmDlz5pj36YW2ToIUEYkPs6EitJi2bw8K/mzF+nifZ7r19g49To3yb90aTI1ZaDvsALvuGlzE6uST4fzzi99rr0JcZBxSbRubN2/OWCiPt6BOvxVyHv/q6mqmT5/OtGnTBu/TH4+2bGRxrYJaRETGq7p6aFrIUhgYCC4ulbpybPpt/fqh+9Fu69YFRX02qW0BPPhgMEf9Rz9a3M+kQlzKxsDAwLBiefPmzeO+jfX6LVu2FLRoTmloaGDKlClMmTKFqVOnDj7O9Dxbgd3Q0KCeahERiR33YDab556DF16ANWuGpp/cvDkokFMj3WM9HhgoTcxHHQVHH138/agQl5Jxd3p7e9myZctgcZte6Ga7H89ty5YtRY1/0qRJNDY2DhbI4ymcsz2fMmWKeqVFRGTCBgaCEz+3bQtuqcejLctl/Wiv7esb3hoy8vnI5b29w6eQLFQBXVUV9HynTqgd7XGqJ3zSpOGPR1s2eTLMmDH8NmtWsJ9SUCEugwYGBsYsiLMVyuN5bTFGldM1NGwi3aMAACAASURBVDTQ2Ng46m3KlCkZ1431msmTJ6toFhFJEPeg+ButqMz1vhDbmMi+xlpfqtHhQpg2DfbeG/bZJ7i4U1MT7LgjTJ06eiGdqdCuKbPKtcw+Tnnbtm1b0QrkzZs30ztW81QB1NfXM3ny5MHiNpf7sQroyZMna15pEZFxSh9NjXNhms++K0FNTTCPem3t2I/zWV9XN1QQp26jLUsvmlMzw0ydWn4FdKHosJTI1q1b6ezsZP369WNOKZdp2rmeEswPNNEiebSieeSyhoYGavQvUUTKgHvw5/dssz309cWvME2/T9Jo6kRVVb2xyCz2fbG2Mdrj6urym1e70qgqykNvby+vvPIKK1euZOXKlbz66qt0dnYOu3V0dNDZ2cmmTZvy3l9VVdXgyG8hiuWRRbNO9hORctbbC6++GlwGPXXZ840bR3+8aVP2qdR6e4NiPOmSVqDmeq8/kkrcqRAfpxUrVnDvvffy+OOP88QTT/DSSy+xZs2acb+/traWpqYmZsyYwfTp03OeZm769OkqlEVExqGzE+6+G554Ap5+Gv72N1i9Opi+rJDq6rJfVKS+Pt6FqkZTRaKnQnwMjzzyCJ/97Gd58MEH37CupqaGnXfemXnz5jFv3jxaWlpobm6mqalp8D71eOrUqSqiRUSKqL8fLrwQrrtu9P7gmhrYaadgRoT0q+lNnx70saY/nzIl+1X6Jk0q/oU+RKT8qRDP4tFHH+Xwww9n69atTJ48meOOO463v/3tvOUtb2H//fdnp5120kwaIiIxccMN8N3vBu0IxxwD73wnHHQQzJ8PLS3BLA1qVRCROFEhnsXnP/95tm7dyhlnnMF1113HlClTog5JREQyWLw4uP/ud+FTn4o2FhGR8dDYQAbLly/n3nvvZfLkyVxzzTUqwkVEYq6rK7jfbbdo4xARGS8V4hnccsstAJxwwglMnz494mhERGQsqRleJ02KNg4RkfFSIZ7BvffeC8CJJ54YcSQiIjIeqWuSNTREG4eIyHglvkfczBqAK4BuYAZwsbt35bPNgYEBHnroIQAOO+ywvGMUEZFAMXJ2ikbERSRpymFE/FrgAXf/InAbcH2+G3z11Vfp6uqiqamJnXfeOe8ARURkUMFzdkp/f3CvC/iKSFIkuhA3sxbgdOD2cNE9wAlmNj+f7b788ssA7LHHHvlsRkRE0hQrZ6ekZpPdvr0QWxMRKb5EF+LAQqDT3bcAuHs/sAw4Ip+Nrly5EoBdd901z/BERCTNQoqQs1NUiItI0iS9EJ8LrB2xrBtoyWejXeEcWDNmzMhnMyIiMlxRcnZKqhBPtaiIiMRd0gtxB3pGLKsHhl3cuKOjg9bW1sFbe3t71o1u2rQJgMbGxgKGKiIyPu3t7YP5CmiKOp4CKkrOTkmd0vPss/kHKiIyXvnk7KSf0rKK4Kz7dFOB1ekLmpubWbJkybg3OjAwAECNzvgRkQi0tbXR1tYGgJl1RhxOIRUlZ6csXAi//jX84Q/wiU9MNEQRkdzkk7OTPiJ+H9ASToeFmdUC84DFkUYlIiKjKWrOXrgw3Ml94F6ILYqIFFeiC3F3fw1YBBwVLjoauNXdV0UXlYiIjKbYOXvBAthxR3jlFXjxxUJsUUSkuBJdiIfOAU40sy8B7wXaIo5HREQyK1rOrqqC444LHv/mN4XaqohI8SS+Cdrd1wNnRx2HiIiMrdg5+6ST4Oc/h9tug4suKtZeREQKoxxGxEVERAA4/nioq4OHHoLXX486GhGR7FSIi4hI2Zg6FY48EgYGYNGiqKMREclOhbiIiJSVk04K7m+5Jdo4RETGokJcRETKyimnQG0t/O53sEpzaIlIjKkQFxGRsrLjjnDCCUF7yk9/GnU0IiKZqRAXEZGyc+aZwf0Pf6iL+4hIfKkQFxGRsnP88bDTTvD883DvvVFHIyIyOhXiIiJSdmpq4Nxzg8dXXRVtLCIimagQFxGRsvTP/wyTJsEdd8Bf/hJ1NCIib6RCXEREylJTE3z848HjK6+MNhYRkdGoEB/FXnvtxQknnMB+++0XdSgiIpKHiy6CAw6Ao4+OOhIRkTeqiTqAODr55JM5+eSTow5DRETytOee8PTTYBZ1JCIib6QRcRERKWsqwkUkrlSIi4iIiIhEQIW4iIiIiEgEVIiLiIiIiERAhfgo2tvbow6hpPR5y1clfVaovM8rw+n7z52OWe50zHKnY5aZCvFRVNp/MPq85auSPitU3ueV4fT9507HLHc6ZrnTMctMhbiIiIiISATM3aOOoejMrANYkcNbmoDOIoUTR/q85auSPiuU5+fd1d2bow6ilCaQs1PK8fsvNh2z3OmY5a6SjllOObsiCnERERERkbhRa4qIiIiISARUiIuIiIiIRECFeBozazCzq83s/5jZdWY2LeqYis3MLjMzD29Lo46n0MzsGDP7k5nNT1tWtt/zaJ83XF5W37OZnWxmy8xsrZl9x8xqwuVNZna9mX3DzK4ys9qoY5XiK+d/04U0Wh7QsXujXH83dAxz/+3RMRuiQny4a4EH3P2LwG3A9RHHU1Rm1khwAsUx4e2UaCMqLDObDUwBDhmxqiy/50yft9y+ZzObB5xE8Dk+A5wFXBCuvgW40d0vBl4Avh5JkFJqZflvupCy5AEduzQT/N2o6GM4wd+eij5m6XSyZsjMWoBlwAx33xKOsG0EDnD35ZEGVyRmlipevufuvZEGUyRmVgVsB3Zz9+Xl/j2P/LzhsrL6ns3s3cDD7t4fPr8COAC4FPiVu7eEy+cCzwOz3X1zVPFKcZX7v+lCGS0P6NiNLpffDaAv07pKOoa5/Pbov7vhNCI+ZCHQ6e5bAMIf+WXAEVEGVSxmZsCHgW8Ba8zstIhDKgp3HxixaCFl/D2P/Lzl+D27+/2pIjy0ClgJHEnwXaZetyp82FrC8KT0FlLG/6YLIUseWIiO3Rvk+LuRbV3FyPG3ZyE6ZoNUiA+ZC6wdsawbaIkglqLzwDsI/mx0NfAzM/u7iMMqBX3P5fc9vw34HhX23cogfe9jyJIHdOzGJ9tx0jEcxRi/PTpmaWqiDiBGHOgZsawe2BZBLCXj7uuBS8L/e70AuDPikIpN33MZfc9mtifQ4e5Lzawiv1upzH/TEzFKHrgbHbvxGOu/MR3DDDL89ujfbBqNiA9ZBcwYsWwqsDqCWKLwXWBe1EGUgL7nMvmezawa+CTw+XDRsO82TPyNVM53W6kq/d/0RKTygI7d+GQ7TjqG45P+26NjlkaF+JD7gBYzawAIpz2bByyONKrSGQAejzqIEtD3XD7f80XAN919a/j8bmCftPXzCU6kWlLiuKS0Kv3f9ESk8oCO3fhkO046huOT/tujY5ZGhXjI3V8DFgFHhYuOBm5NO+GrrJhZi5mdbmZV4cjhhcCXoo6r0MLPBmBQ/t/zyM9brt+zmX0ReAJoMLPdzewfgQ3Aq2b2pvBlxwNXu3tfVHFK8ZX7v+lCyJQHdOxGl8vvho5hIJffHh2z4TR9YRozmwFcASwH5gBfcPdNkQZVJGZ2AEGvVhfwAPBtd38u2qgKy8ymAmcQzFf6deAad+8o1+95tM8LzKLMvmczu4RgqsJ0z7n7fma2K/AVgu92CkGx0Y+UtXL9N10o2fK9jt1wE/ndqPRjOJHfnko/ZulUiIuIiIiIRECtKSIiIiIiEVAhLiIiIiISARXiIiIiIiIRUCEuIiIiIhIBFeIiIiIiIhFQIS4iIiIiEgEV4iITYGZHmdlDZrZwjNe9ycxuNrMvF3DfVxZqWyIiUTCzN5vZLYXMjRn2s8jMPlzMfRSLmR1tZg+P9TsjyVYTdQAi+TCz/YDLgXXA24A3Ax91958WedevAAeO43UrgRlAdSF2Gl4KeFshtiUiAmBmbwP+HTgc+B2wjOAiK7XAde7+2yLs9iVgOgXKjVlcC/y5GBsO8/HJwPeAF4CHw1X1wAeAQ9x9eR67GO/vjCSYCnFJul8DH3T3pWZWBfywFDt19xfMbO04XtdlZq8WcNdHAH8s4PZEpMK5+6Nm9ktgX3f/+9RyM3sf8P/M7CJ3v7bA+9xS4NwIgJn9A/B0qgB29zsKvY8Ud98G3GRm5wF/cPdL0uK4bSLbNLNZwHvc/SZ3f97MOgsUrsSUWlMkscxsJrA3sBnA3QcILm9eqsvFDhT4dePxHuC+Am5PRASgf+QCd/8f4N+A/2tmc4qwz4LmajPbmWAEvNTecOyAe4BVuWzEzOqBHwMNaYsL+fshMaRCXBLL3dcBzwN3mdk7wmUrgEfMbAczu9TMHjGzE8zsNTP7i5m9NfV+MzvNzL5kZn8ws/+TtrzVzP7NzG4ys1+Z2fRw+SQz+08z+4qZ/YDgz6o5MbN/DWP6oJm9zcx+HW7v0jDGm8ysOVy+0cy+OmITde6+NdzWhWb2WTN70Mz+I+cDKCIyth8QtKicCmBme5jZ183s+2b2+7D4xczqw1z2JTO738w+lNqAmZ1rZl82s/81s09n2tFoOdnMTglz5plm9qKZ/YeZtZjZD8N8er+ZHRlu4nRgZ+BiMzspPJfnHjP7eNo+zjazy83sx2b2X2bWaGbTzeyrZrbEzI40s+fN7GUzmz+RAxYek4XhiDlmdrKZfcfMvmVmd5nZlPAzfMfMbg9/b14Nj/H+wKlm9tm0Tc4If4u6zezsDPusMrNvmFlbuA83s/+X+n4kxtxdN90SewMOAFYQjBrcCMwKl1cTJLV1BL16OwB3A88Q/A/oocAl4WtnEIyqvw+YBvwybfsPAleGj/8v0BY+ng1sJ0i2Y8X4I4JRpVkEvYSTw+X1wG+B24AWgh+QbcBXgcnAweHz1OtT/e+pz31t+Hgq8PmovwvddNMtuTfgTOBvGdatAb4b5tVfAdXh8p8Dt4SPrwJODh//I9ARPj4NOCN8vHeYqxeEz38E/Fv4OFNOnglsAL4NHAIcDXwT+Er42s8Bt6fF6sD88HED8BfgzPD5ccDP0167CLieoE33/eHvxcnh57wLuHycx+4P4W/FfwBXEvSLH522vhOYFz5+huA3qR74EvA34O3AP4V5/w+peMPXLwe+Fa47DViZIYazgEXh49nhcT426v+udBv7ph5xSTR3f8bMFgCXAecC/2BmR3vQM94BdLn7rQBm9m/AA8DuBD86A2Z2Qbip3xIk/PcBO6Yt/0v43inAOeF7cfc1ZvZKDqE2AVcD/+TuW8JtbDWz14Bl7r463M9rwGIP+iefIPiBaCI46fN9BD8aAFuAj5nZC8B1BP8TIiJSDE5QnB4C7AqcZ2YAa4HpFjz5JMEgAsB/EbRmQJBrl4zItTsBT4/Yx5mMkpPdfZ2ZbQRuc/dHAMxsJbDZzHYiGEGeMmrQ7j1m9nraok+H2035PnBzuHwDwe/Fr8J9LAnjHK/7POwRN7ObCAZ/Uk5195Vm9i6CAnxKmP9XA3919z8BfwrfO9q2bw9/Ex7NEtNbgZ7wc68xs3UMb3GRmFIhLollZtXALh6clPMvYbvIHUA7wQjDSM+E9zsA84CfuPvPw2XfDrd5MUFh/O0R+zqYIIF2TzDcvYHDgD2ApWnLR/ZIbh9c4T4QJuVUC9mOHrTj4O7LzOxM4BrgM8BHgNcRESkgM9uB4K95zxDkzbWj5MfZBCO2dQDuvp1g8IDwPV9z94fC58Pem2bUnBxyhufKvwEXAy8C/wvsluUjpL9vL+DetOd/JWi7mcUbc3E/E5/R5XGC/2FJ2WBm3wBuIYg9VW2P/Fxj6Sdz3XYncJ2ZTSL4LnoIRtcl5tQjLkk22LcI4O5LgS8C+2Z4fR1B0nsJeBU4JX2lmb09XP4+M2sYsbwrfLrfBGN9CPgO8LMwUebEzFpIO/HHzOYSzBizN8HJmz/P8FYRkXx8nOAvcLcQ5Md3hSPRwGB+7AT6gOPTli8ws6mMyLVmVhMObIyUKSeP5nvAcg+mqd2aw2dZSTCCPrgLYCNB603BuPtAOFhyQPg/MncAX3f3JYXcz4h93g7cTvCX29OAI9x9Y7H2J4WjQlyS7rywSE2ZS9DblzLTzOrCxwsJ+gM3AL8A3m9m15jZwnC0YjNBwmwEFpnZceGJMfsSFO9PAJeaWZ2ZNRH0Zs8ys2oz+7yZnZEhRgtvXyH4sbpqlHUjXz/y+d8T9DOm7AWc5O5dwHmjvEdEJBdvGGk1s8OBrwGfdPcOgtHnNcAdZvYP4QmZ7w5HwH9JMLvKSWZ2FPBxd+8myLX/YmZfM7MjCAYkUm196fkvU05OSR+dfgvQbGaTgXcBDWa2e7huG8HJjakBmfR9fA/4gJnNCJ+/HbjB3fsJ6qFRc7GZvSOMqz7DsatiRD0V/s/K+wj+CjqDoOVxH2AXYHLaiaAjR937Rol/eFCj9K+E38U04PcEo/65/A+KREitKZJ0uwDPmNn/EPzZzoB/TlvfD3zJzDYQjB6fB+Dud5vZvxD8efMkgpMd/wxgZicC/0nww3Ij8Dl3dzP7IMHUUs8ANxH8ifEgggtgfJCg5WXYhYTM7BCCk5DmEfyp8tfA18J4fgu8A5hvZnsQjNTMAT5oZs8xNDp0JtDo7u0jPvsPLbgQRz/ByVEiIjkL88jpBMXiDUAHwQl/04G/c/eHAdy9L8yPNxDkutsIesMBzicYxPgxwV8AUzOV/IAg951D0EJ3jru/Hhaa7wD2MLOfZsrJZnYKQV/0WWb2nLuvIjhx9JsERfiPgROBBcDLwE8IRu8/ZsFsWgcA7zWzO9z9V2a2F3CrmT1I8FfVL4Uj96cCc8LP9wzBxY12MrMDwjg/TND3/uu041ZLkPsPBOaGhboT9GafHG7zSWAx8ChBu82dBDl9EcEJogea2UnuntruTQQXV9oUFu5zgFPMbBnBSDfh+0deM2MlwcmoJwGTgCozu4fghM1STekrE2D6fqRcWXBZ4B+5+/wS7e9z7q7Lz4uIlJmwTabe3WN5QTUz+xzBDDbLw+eTCP6a8Y3UuUUST2pNkXI2WttHcXYU/An38VLsS0RESic8GfWtMS7CJwGXAPPTFjcBm1WEx59aU6Qshf15ZxD8qfFUd7+lyLv8a/gnUxERKS817n5d1EFk4u69ZnYe8N3wRP7VBG00l0QbmYyHWlNERERERCKg1hQRERERkQioEBcRERERiYAKcRERERGRCFTEyZpNTU0+f/78nN/X0dFBc3Nz4QNKCH1+ff5K/vwQj2Pw2GOPdbp7RX0RytnFoeOTnY5Pdjo+2aWOT645uyIK8fnz57NkSe5Xlm1tbZ3Q+8qFPr8+fyV/fojHMTCzFZEGEAHl7OLQ8clOxyc7HZ/sUscn15yt1hQRERERkQioEBcRERERiYAK8Sza2tqiDiFS+vz6/JVOxyBZ9H1lp+OTnY5Pdjo+2U30+FTEBX1aW1tdfU0ikkRm9pi7t0YdRykpZ4tIUuWaszUiLiIiIiISgZLPmmJmxwCXAR9y9+XhsgbgCqAbmAFc7O5d+awTEZH8KWeLiBRPSUfEzWw2MAU4ZMSqa4EH3P2LwG3A9QVYl79774WDD4bPfKagmxURSYKk5OyuJV08dshjPPePzxVicyIiJVPSQtzd1xAk30Fm1gKcDtweLroHOMHM5k90XcEC3rwZHn8cXnyxYJsUEUmKxOTsAeh+tJvNT2/Oe1MiIqVU8tYUdx8ws/RFC4FOd98Sru83s2XAEcC2Ca5bXpBg6+uD+61bC7I5EZGkSULOrpkZ/JRtW7stn82IiJRcHK6sORdYO2JZN9AC9E9wXWFMmhTc9/YWbJMiIgkXu5xdu2MtANvWqRAXkWSJQyHuQM+IZfUEoyfksW5QR0cHra1DM8m0tbWNb75HjYiLSATa29tpb29PPW2KMpZRxC5n10yvgSrYvnE7A/0DVNVoQjARKZ18cnYcCvFVBGfPp5sKrCZI+BNZN0xzczMTmpNWI+IiEoH0wtPMOiMOZ6TY5WyrMmpm1NC/tp/+9f3UNdeN+70iIvnKJ2fHYdjgPqAlnNYKM6sF5gGL81hXGBoRFxEZKZY5u3Zm0J7Sv66/EJsTESmJkhfiNnTWjwG4+2vAIuCocPnRwK3uvmqi6woWrEbERaTCJSVnD/aJ64RNEUmQkrammNlU4Izw6Zlmdo27dwDnAFeY2YHAHCC9GXCi6/KnEXERqWBJytmDM6fohE0RSZCSFuLu3g1cF97Sl68Hzs7wngmtKwiNiItIBUtSzh5sTVmr1hQRSY449IjHl0bERUQSoWZHjYiLSPKoEM8mVYj39oJ7tLGIiEhGqRFx9YiLSJKoEM+muhpqg+TONiV3EZG4Sp2sqVlTRCRJVIiPJX1UXEREYkmXuReRJFIhPpbUCZvqExcRiS2NiItIEqkQH4tGxEVEYk/TF4pIEqkQH4tGxEVEYk8X9BGRJFIhPhaNiIuIxF5tU1iId6oQF5HkUCE+Fo2Ii4jEXvWUaqomVTGweYDtm7dHHY6IyLioEB+LRsRFRGLPzKhtDkbF+zr6Io5GRGR8VIiPRSPiIiKJUDsrbE95Xe0pIpIMKsTHohFxEZFEqJtVB0Df6xoRF5FkUCE+Fo2Ii4gkgkbERSRpVIiPRSPiIiKJoBFxEUkaFeJj0Yi4iEgiaERcRJJGhfhYUiPiKsRFRGJNI+IikjQqxMeSGhFXa4qISKxpRFxEkkaF+Fg0Ii4ikggaEReRpFEhPhaNiIuIJIJGxEUkaVSIj0Uj4iIiiVDXHIyIb+vYhg94xNGIiIxNhfhYNCIuIpIIVfVVVE+vxvud/g39UYcjIjImFeJjaWwM7jdvjjYOEREZU2pUXH3iIpIENVEHkGJm+wGfBl4C9gLa3f1JM2sArgC6gRnAxe7eFb4n47qCmTo1uO8q7GZFRJIsrjm7dlYtPS/1sG3NNti3kFsWESm82BTiwE+AE919lZnNA34L7AdcC9zl7jeZ2fHA9cBp4XuyrSuMadOCexXiIiLpYpmz6+aEI+KvaURcROIvTq0p+wHh8DM9wHQzawFOB24Pl98DnGBm87OtK2hUqUK8u7ugmxURSbhY5uz6ucEJ9ltX6wR7EYm/OBXi/w1838ymAmcA5wELgU533wLg7v3AMuCIMdYVjlpTRERGE8ucXdcSjIhvXaVCXETiL06F+LlAL/AosMndbwXmAmtHvK4baBlj3TAdHR20trYO3trb28cflVpTRKTE2tvbB/MV0BR1PBnEMmenRsT7Vqk1RURKI5+cHace8QbgFwTJ+ttmtgxwgj95pqsHUldryLZuUHNzM0uWLJlYVCrERaTE2traaGtrA8DMOiMOJ5NY5my1pohIqeWTs+NUiP8U+LC7bzAzJ0jwFxCcWZ9uKrCaIOFnWlc4KsRFREYTy5yt1hQRSZJYtKaYWRNwoLtvCBddBnQBK4CWcMorzKwWmAcsBu7Lsq5w6uuhthb6+nR1TRER4p2zB1tTVvfhrqtriki8xaIQB9YBvWY2N23ZWuApYBFwVLjsaOBWd1/l7q9lWlfQyMw0Ki4iMlxsc3bN1Bqqp1Qz0DOgq2uKSOzFojXF3QfM7CTgK2b2GDAb+Jy7d5nZOcAVZnYgMAdoS3trtnWFM20arF0bFOLNzUXZhYhIUsQ9Z9fNraPn+R62rtpK7YzaYuxCRKQgYlGIA7j7UuCToyxfD5yd4T0Z1xXU9OnB/caNRd+ViEgSxDln18+tp+f5HvpW98Gbir03EZGJi0trSrylWlNUiIuIxF59Szhzik7YFJGYUyE+HjPCE/3Xr482DhERGVP9LmEhvlKFuIjEmwrx8dhxx+B+7chrUYiISNxMmj8JgN4VvRFHIiKSnQrx8Zg5M7hXIS4iEnuDhfhyFeIiEm8qxMcjNSK+bl20cYiIyJgm7apCXESSQYX4eKg1RUQkMernhT3ir2zFt+uiPiISXyrEx0OFuIhIYlQ3VFM7uxbvd7au1gmbIhJfKsTHI9UjrtYUEZFE0AmbIpIEKsTHQyPiIiKJohM2RSQJVIiPhwpxEZFESZ2wuXWFWlNEJL5UiI9HemuK68QfEZG404i4iCSBCvHxmDQJJk+GbduguzvqaEREZAyDUxiqR1xEYkyF+HjNmhXcv/56tHGIiMiYJu0eFOI9L/VEHImISGYqxMdr9uzgXoW4iEjsNezeANVBa8r23u1RhyMiMqqCFOJm9t5CbCfWdHVNESkTlZCzq+qqaNitAVyj4iISXzW5vNjMjgcuAeYwVMRXAbOBhsKGFjOaS1xEEqaiczbQsE8DPS/10PN8D1PeNCXqcERE3iCnQhz4CfDvwBPAQLjMgNMLGVQspXrEV6+ONg4RkfGr3JwNTN5nMuv+Zx1bnt8SdSgiIqPKtRB/2d2vGrnQzJ4rUDzxte++wf2zz0Ybh4jI+FVuziYoxAEV4iISW7n2iJ9jZqeMsvyEQgQTa296U3D/1FPRxiEiMn6Vm7NRIS4i8ZdrIf4u4Ptmtj3tNgBcV4TY4uXNb4aqKvjzn6FHJ/6ISCJUbs4m6BEH6Hm+B9fF2EQkhnJtTfk60AYsAVLzQRnw0UIGFUuNjbD//kEh/uSTcOihUUckIjKWys3ZQN3sOqqnV9O/oZ++NX3Uz6mPOiQRkWFyLcTvBG5z963pC83shkIFZGbvAA4DXgIeADYDVwDdwAzgYnfvCl/bkGldUbztbUEh/uijKsRFJAkqOmebGVMWTGHjAxvZvHSzCnERiZ1cW1N+CXzBzA5Pux0BfKEQwZjZJ4C/d/dvuvuv3b0TuBZ4wN2/CNwGXJ/2lmzrCq+1Nbh/4omi7kZEpEAqO2cDU94STFvY/UR3sXclIpKzXEfEvwzsBZwJpBrujGCO2vPzCST8cfgIcEzashaCabY+FS66B7jVzOYDfZnWufvyfGLJaO+9g/uXXy7K5kVECqyyczZDhfimJzYVaxciIhOWayF+JXCnuw8bWijQVdquAh4CrjazPYFLgV2BTnffAuDu/Wa2DDgC2JZl3fICxPNGe+wR3KsQF5Fku33eYQAAIABJREFUqOycjQpxEYm3XFtTXhmZ0ENz8gnCzPYB3gLc4O6fBn4P3AXsAqwd8fJuoAWYm2XdMB0dHbS2tg7e2tvbJxboLrtAdTWsWgW9vRPbhojIGNrb2wfzFdCUx6YqO2cDjfs3YrVGz4s99Hf3T3g7IiKZ5JOzcx0RPxd4OH2BmZ0KfBu4McdtpTsAWOfuqUm6rwG+SvCn1JFzBdYTjKwwxrpBzc3NLFmyJI/wQjU1sPvu8OKL8MwzcPDB+W9TRGSEtrY22traADCzzjw2Vdk5G6iqq6LxTY1semITm5ZuYofDdijIdkVEUvLJ2bmOiM80sy+HO5psZt8Hvg88muN2Rqoh7X8K3L0HeBGoJjizPt1UYDWwKsu64nn3u4P7++4r6m5ERApAOZu09pTH1Z4iIvGSUyHu7u8D7jazq4Engb2BA4HRrtyWi6eA6WY2O21ZP0GCbgmnvMLMaoF5wGLgvizriuc97wnuVYiLSMwpZwemHTINgK6Hize7rYjIRGQtxEdMeXW4mR0O1AErCEY/vkpwcs5X8gnC3Z8jmO/2lHC/O/z/9u48Tq6qzv//69P7lnRn6ewJCSQkYUciIkH2TVEUB9BRx6/4lbjBqF9RFBmdQeYng+iIiqNxQ4YZVxgQxlEQUFaVMMgeCCRh6ZDQWbrT6X35/P44t9LVla7qru7qvtVd7+fjcR5Vde+te8+9Vfnk06fOOTc6zn8AtwOnRJueCtzk7g3uvjXdutHUZUiJRPzee6F7n19URURio5g9uNrjagFouq9Jd9gUkbwyVB/xq4GVQBPQN8j6HxGmwpoLfGqUdXk/cK2ZlRFaSf42GlX/UeBqMzucMMBoTdJ7Mq0bG/Pnw4oVsH49PPQQHH/8mB9SRGSYFLMHUbWyipK6Eroauuh8qZOK/SrG47AiIkMaKhG/Aljn7q9l2sjMRvszJ9GNIN47yPJdwIVp3pN23Zg644yQiN9xhxJxEcknitmDsCJj6uqp7PzvnTQ/0KxEXETyRsauKe7+m6ECerTdr3JXpQng9NPDo/qJi0geUcxOL9E9pfm+5phrIiLSL9tZUwTg6KPD41//Cr298dZFRESGVHd8mLZw1127Yq6JiEg/JeIjMXMmLFoEbW3w7LNx10ZERIYw5egpFNcW076hnfbNqdOZi4jEQ4n4SB1zTHi844546yEiIkMqKili2ilhGvNdd6pVXETyw6gTcTM70sxuNLNPm9ni0VdpgjjnnPD4y1/GWw8RkSwUbMwGpp8+HYBddygRF5H8kIsW8Y8Ap7r714BjEzdrmPTe+laoqIAHH4RXXom7NiIiw1WYMRuYdnrUIv77XXiv5hMXkfjlIhG/CDgqev4zoDQH+8x/NTXw5jeH5zfdFG9dRESGrzBjNlC5pJLKZZX0NPXQfL9mTxGR+GWdiJvZm83svdHzo4FliTujuXufuxfOPYTPOy88KhEXkTylmD3QzHNmAtB4U2PMNRERyTIRN7N/Bn4JnAPg7n8B3m9m7xmDuuW/t7wFiovh/vs1aFNE8o5i9r7qz60HQiLufeqeIiLxyrZF/FhgAXBf0rL/Aa7KWY0mktpa+OhHwR3OPBMuu0zziotIPlHMTjFl1RTKF5XTtaWL3X8qqB8DRCQPZZuI3+nuTUByM8Kob5U8oX3jG3DFFVBUBF/5CnzgA3HXSEQkQTE7hZntbRXfduO2mGsjIoUu20R8o5ldBiw3s3eb2U3Ax4F/yX3VJojiYviHf4A774TqarjxRnjkkbhrJSICitmDmnPBHAC23rCVnuaemGsjIoUsq0Tc3X8G3A0UA+8HmoDT3f26MajbxHLSSaEAbNoUb11ERFDMTqfmkBrqTqyjr7WPrTdsjbs6IlLASrJ9g7v/CfhT4rWZFZvZKe5+V05rNhGVRrOAfeUrUFcHp5wCZvHWSUQKmmL24OZfNJ+mPzTR8O0G5n98PlakWC0i4y9jIm5mvwZqhtjHTGA7UNBBHYDLL4f77oP//V847TQ46ii49FI491wl5CIy5hSzh2/G22dQvqCc9ufa2f7r7dS/oz7uKolIARqqRfzxaJtngT7gQuBRYF3SNkuA6WNSu4nmda+D55+H664LgzgfeQTOPz8M4PzhD8OAThGRsaOYPUxFJUUs/OxCnv/759n0hU3MfNtMrFgNJiIyvoZKxL8KWDTqHjM7wN2/mLqRmd06FpWbkGprwzSGn/oU/PjH8NnPwvXXh+X/+q9qGReRsaSYnYV5a+bxytdfoe3pNrZ8fwvzPzI/7iqJSIHJ2ETr7s2JgB6ZbWbFyduY2YnA0WNQt4mtshI+9jH49a9D3/Frr4Wvfz3uWonIJKaYnZ2i8iL2v3p/ADZ9fhNd27pirpGIFJps+0r8CnjMzL5rZl81s9uBO4Fv5r5qk8TJJ8NPfhKeX3IJfPWr4QZAIiJjTzF7CPXn1jP9zOn0NPXwwmdeiLs6IlJgsp2+8E7gdMJAn+WEqbDOdfevjEHdJo+//Vv4ZvT/3mc/G1rKezR3rYiMLcXsoZkZy769jKKKIrb9+zZ2/GZH3FUSkQKS9ehBd9/i7pe7+9nu/j7gNjM7agzqNrlcfDH8/OdQXg7f/S68/e2wZ0/ctRKRSU4xe2iVB1Sy+J8WA7D+gvXqoiIi4yarecTN7EeDLK4HeoF3jLYyZlYGPAx8wt3/YGaVwNVACzANuNTdd0fbpl2Xt84/H+bPD0n4b34Dq1fDzTfDAQfEXTMRmYTGOmZHx5gUcXvhJQvZ+dudNN3TxPoL1nPo7YdqbnERGXPZtogvA14CXkwq5UCu7un+GWBx0uvvAPe7+2XArcD3hrkuf61eDQ89BMuWweOPh7nGb7st7lqJyOQ01jEbJknctiJjxQ0rKJlWws7/2cmLV74Yd5VEpACYZzFw0MyWu/uzKctWAB9w98+NqiJmxwIrgC8CHwCeAzYB09y9zcxKgGbgYKAr3Tp335y671WrVvm6detSF8eruTnML37LLeH1NdfApz8da5VEJP+Y2SPuvmqE7x2zmB3ta0zidpwxe8dvdvDEW58Ah4NvPpj6c3SjHxEZvmxjdraDNZ8dZHEn8KFs9pPKzKoJA4iSf0Y9Edju7m3RsXsIQfyEIdZNDLW1oVvKV6IxU5dcErqs3HMP9PXFWzcRmRTGKmbD5I3bM94yg/3/JUxp+MzfPUPLoy0x10hEJrOsEnEz22RmG5PKi8AzwH+Nsh6XAlelLJsPpA5fbwHmDbFuH42NjaxatWpvWbt27SirmyNm8LnPwb//O1RXhznHTz4ZliyBT34S7r4burvjrqWIjLO1a9fujVeEW9KPyBjGbBjDuB13zF54yUJmv282fa19PH7G47Subx3X44vIxDKamJ3VYE3C/LP/CST6s/QB29z9uSz3s5eZnQGsc/fXUlY50J6yrBxIZKaZ1g1QX19P3nVNSfa+98Gpp8J118ENN8BLL4UbAF17bWg5P/NMOPvs8Di94O9MLTLprVmzhjVr1gBgZttHsaucx+yoTmMat+OO2WbG8h8sp6uxi12/28Vjpz7G4XccTvVB1bHVSUTy12hidraDNW909z+4+x+jcl84pi3Ncj/JLgF+ZGbbo8ovJAziqSaMqk82BdgCNGRYNzHNmQNf/jJs2gQPPBBayg86KPQl//nP4b3vhfp6OOYY+NKX4MEHNRe5iAxlLGI2FEDcLiov4pCbD6H2+Fq6Grp4dPWjNN3bNPQbRUSykO1gzf8vGu2evKwOuMfdjxxRBczqCa0iCQ8B/w/4M/A0UO/u7WZWCuwB9idMvfX8YOvcvSH1GHk5WHO4XnghzKpy221w330Du6rU1cEpp8Bpp4XHAw4I3V1EZNIY5WDNnMfsaB9jGrfzKWb3tvfyzHueYfst27EyY+m/LmXeR+dhirUiMogxGaxpZheZWRNwqZn1JhdCn79tI6wv7t7o7q8kCiFYN7r7S8DtwCnRpqcCN7l7g7tvTbdupPXIWwccEPqL33UX7NwZEvKLLgrTHzY1wU03wUc+El7vtx9ccEHoc94w+S6FiAzPWMZsKKy4XVxZzMG/Opj5fz8f73I2fHwDT7/7abqbNH5HREZv2C3iZnYY8H+Ab6as6nL3V3NWIbPNhKm1/mBm0wg3f9gMzAE+7+57ou3SrkuVT60rObVpE/zud2FQ5913w46UcVDLl4eW8ve8B449Vq3lIhPQSFvExytmR8faTA7jdr7G7G0/28ZzFz5H755eyuaWceC/HcjMt494LK2ITELZxuxsu6aUuvs+zQBmtihqCclL+RrUc6qvD554IrSc33UX3Hsv7En6/23lSjjnHHjb2+Doo6Eo2+EBIhKHUXZNUczOsbbn2lj/gfXsfijcELT+XfUs/cZSyueUD/FOESkEOU3EzewS4Dfu/nT0+gtAcepmwEnufmL21R0f+RzUx0x3Nzz8cJgS8frrYVvSL9GzZsFZZ8Gb3xxazDUTi0jeyiaoK2aPD+91Gr7TwMbPb6SvtY/immL2u3w/FnxyAUXlauQQKWS57iO+BKhLev064FjCwJslSWVOlvWUsVZaGrqjXHUVvPxy6MJy8cWweDG89hr8+Mdw/vn9M7F88Ytw//2at1xkYlPMHgdWbCy4eAGvf/L1zDh7Br17etn4uY385eC/0HhzI9n80iwihS3brimvA550966U5QclWmDyUb63rowrd3jqKbj9drjjjn2T76lTw02FTjstlKVL1bdcJEaj7JqimD0Odt65k+c/+TxtT7cBUHNUDUuuWML0N0/X7CoiBWZMb3EPHDpIQF8OHJDlfiQuZnDIIWGe8rvvDjOx3H57aC1fvhx274ZbboGPfxwOPDDc5fPCC+EXv4Dto7mviIjEQDF7HEw/bTqrHlvFsm8vo2xOGXse2cMTZz3Bo6sfZecdO9VCLiJpDatF3MzeROhXuAZIvdfwdOCb7r4o99XLjYnWuhKrF1+EO+8M5a67Bs7EYgZHHtnfWr56NVRUxFdXkQIwkhZxxez49Lb1suXftvDSVS/RvT382lh9eDULL1nIrHfNoqhUfchFJrMxmTXFzA4FbiD0LdyVsroL+A93vyKbio6niRzUY9XXB48+2p+Y338/dCU1rlVUwHHHhQGfJ58Mr3sdlJTEV1+RSWiEibhidsx6WnpouK6Bhmsb6Noa4mb5wnIWfHIBcy6YQ+m00phrKCJjYcymLzSzWuA8d//BSCsXl8kQ1PNCW1u4u2eitfyvfx24fupUeNOb4KST4Pjj4YgjwqBRERmxUcwjrpidB/o6+9h24zZevuZl2taHPuRFlUXMevcs5n10HlNfPzXmGopILo3pPOIZDnqcu98/6h2NkckU1PNKYyPccw/8/vfh8fnnB66vqgpzlq9eHcob3wh1dYPvS0QGNZrBmhn2qZg9zrzP2fHfO2j4ZgO7ft//I0XNUTXM/+h86t9VT0mNflEUmehyPY/4r4GaofYBHOzus4Z70PE2GYN6Xnr5ZfjDH0JSfv/9sGHDwPVmcPDB/Yn5cceF6RQ1q4BIWlnOI66YPQG0bWhjy/e2sPXHW+nZ2QNAUXUR9e+sZ/b7ZzPtpGlYseKiyESU60T8q0AT0ACk29CAt7r7udlUdDxN9qCet157DR58EB54IJR16/adp3zu3P7EfPVqdWcRSZFlIq6YPYH0tvfS+MtGtqzdwu4Hdu9dXja/jNnvnc2c98+h+uDqGGsoItnKdSI+E2hx984hDjrD3Xdk2iZOhRLU815HR0jGE4n5Aw+E6ROTqTuLyABZJuKK2RNU2/NtbLtxG9tu2EbHpo69y6sPqab+vHrqz6uneqWScpF8F1cf8Q/l84CgQgzqE0JfHzz77MDEfKjuLKtXh7nN1Z1FCsQY9RFXzM5T7k7zA81su2Ebr/3iNXqbe/euqzqoKiTl59ZTfXC1bhYkkofGNRE3s0XAx4GL3D1v/1Qv5KA+4aR2Z3nkkYFTJgLMmQPHHguvfz2sWgVHHQXTpsVTX5ExlstEXDF7Yunr7GPXXbto/FUj22/ZTs+unr3rKg+sZMbbZjDzbTOZeuxUzU8ukifGPBE3s2LgbOBC4DRgW7Sf+VntaBwpqE9gqd1ZHnxw4E2GEg44ICTlq1aFBP3II8N0iiIT3GgTccXsyaGvq4+me5p47Zevsf2/tu8d5AlQUlfC9DOnM+NtM5h+5nRKp2ucjUhcxnIe8SXAh4ALgBnAi8DfA78D3ujuD2Zf3fGhoD6JuIfuLH/+c0jQ160L85l3dAzczgyWL+9PzletCgNBq/O2EVBkUKOYR1wxe5Lq6+5j94O72X7bdnbcvoP2Z9v7VxZB7epapp8xnWmnTWPKUVM0A4vIOMp5Im5m5xFaUk4Gngd+APwEuMzdPzWKuo4bBfVJrrsbnn66PzFftw4ee2zfGVqKiuCggwa2nB92WLhDqEieyjqoK2YXnLYNbey4fQc7bt9B873NeE///+sl00qoO7mO6adNZ9rp06hcUhljTUUmv2xj9nDuHtBOuCXyg8CH3P256ECjH+UpkgulpXD44aH83/8blnV2wpNP9ifmDz8cXifK9deH7UpK4NBDB7acH3IIlJXFdjoio6SYXWCqllVR9akqFn5qIT3NPez6/S523rmTXXfsomNTB9tv2s72m7YDUHFARUjKT51G7fG1lNUr1onEKZuuKYsIrSyLgVuBE9z94mideS6mXxkjal0RANrbQ0t5csv5M8+E2VuSlZWFpD6RmB95JKxcqZZzicUouqYoZgvtL7SHpPzOXTTd3URPU8+A9VUrq6g7oY7a42upO6GO8nnlMdVUZHIYj8GaJcA7gPcCTwI/BE5397VZ7WgcKahLWnv2hD7mycn5s8/uu11xMaxYERL0ww7rb4GfM0dTKcqYysFgTcVsAaCvp489j+xh5507abqnid0P7aavfWBDRMUBFdSdUEfd8XXUnlBLxX4VmiZRJAvjPX3hUuBi4EJ3rxrxjsaYgrpkpbkZHn10YH/z557bt+UcYObM/qQ8kaCvXAnlalWS3Mjx9IWK2bJXX1cfLetaaLq3ieY/NtP8QDO9Lb0DtimbV8bUY6Yy9Y1TqX1jLTWvq6G4sjimGovkv7hu6PN+d79h1DsaIwrqMmptbfDUU/D44yExf+yx8Lypad9tS0r2bT0/9FCYO1et55K1Mbqhj2K27KOvp4/Wx1pp+mNTSM7vax4wTSKAlRo1R9Qw9Y1T9yboajUX6RdLIp4LZvZO4GvAVOBG4NPu3hPdsvmfgSagFLjU3buj96Rdl0xBXcaEO7z8cn9SnkjQN2wI61JNnRqmVEwuK1bAsmXqfy5pjUUinguK2ZOfu9P+XDvNDzWz+6Hd7P7TblqfbIXUYTVzypjyhilMOaq/lM3WIFApTBMyEY8GFV0JXAscDPwb8CV3v8bM7gE+5+5/NrOPAPu5++ej96Vdl0xBXcZVW1uYmSU5QX/ySdi1a/DtzWC//QYm6EuXwpIlsGiRurkUuHxMxBWzC1dPSw8tf2kJyfmfdrP7od37tJoDlC8op+aompCYr4qS81lKzmXym6iJ+JuAh9y9J3p9NSG4fxm42d3nRcvnA88Cs4FD061z99bk/SuoS+zcYfv2MBB0/frwmCgvvAC9vYO/zwzmzQtJ+ZIlsHgxLFgQBonOmRO6u8yerekWJ7E8TcQVswWIWs03tNPycAstj4Sy53/30Ltn35iWSM5rjqih5rAaqg+rpnL/SqxI3Vpk8hiLecTHnLvfl7KoAZhCuCHFpqTtGqJ+aKuA1RnW/XGMqyySHTOorw/luOMGruvqgo0bBybnGzfCpk2h60tDQyj3359+/9Onh6Q8kaDPnh0GkiaX+vrwOG1amAVGZIQUsyXBzKg6sIqqA6uY/d7ZAHif0/ZcG3se2TMgOe98pZPOVzrZceuOve8vqiqi+tDqvYl5zWE1VB9aTem00rhOSWRc5UUiPojXA18F1gA7Uta1APOA+RnWiUwcZWWhr/iKFfuu6+6GV14JSfmmTbB5M7z6aihbt4aybRvs3BnKU08NfTyzkLgnEvPURD25TJsWSl1dGIQqMjjFbNnLiozqFdVUr6jeNzn/3z3seXwPrY+3sufxPXQ1dNHy5xZa/twyYB9l88uoPqiaqhVVoawMj2VzyjQwVCaVvPufNZpeq9HdH4vuBNeeskk50A1kWjdAY2Mjq1b1/0qwZs0a1qxZk9N6i4yJ0tL+binp9PbCjh39ifmrr0JjYyjbt/eXxOtdu8L2O1JzoiFMmdKfmA+n1NWFAaq1tVBZqRljsrB27VrWrt07zffMOOsyFMVsGY4Byfl7Zu9d3r2jmz1P9CfmrY+30vpkK10NXXQ1dLHrzoFja4pri/uT8xVVVK+spnJ5JZVLKikqLxrv0xIBRhez86KPeIKZFQNXAZe7e6eZfR44yd1Pj9Yb0AGcBJyQbp27P5i8X/U3FEnS0xNaz9Ml6smvd+0Kpbl58Jlghqu4OCTl6Uptbeb1U6aEUlNTcN1q8rGPeIJitowF73XaN7XTtr4tlGf6H3t27TswFACD8oXlVB5QSeUBlVQcUBGe7x+el9apq4uMjwnZRzzJJcA17t4Zvb4T+EjS+sVAF7Aueky3TkTSKSmBWbNCGa6+vpCM79oV5k5PJOiZSlMTtLTA7t3Q3t6/fLSqqgYm59mW5PdWVKilfnQUsyXnrNioWlpF1dIqeGv/cnenu7F73wT92TY6Xuqg86VOOl/qpOmefe/vUDK9pD9JX1xB+aJyKhb1P5bU5ls6JIUib755ZnYZISBXmtn+wInAvcCrZnaIuz8JnAl8y927gHVmlm6diORSUVF/l5OR6O4OSXlzc0jMM5XUbZqbw3tbWmDPnjA9ZFtb6IYzWsXFA5P02lqYMaO/zJw58PWMGWFQ7IwZBZ/AK2bLeDMzymaVUTarjLrj6was6+vuo+PFDjpe6KD9hfa9JfG6Z2cPLTtbaHm4ZdB9F9cWD0jME48V+4XnZXPLKCpR1xfJvbxIxM3scsK0V8nWu/uPzOxdwBfNbDNQA3whaZtM60QkX5SWhgGi06ePbj99fdDa2p+Ypyu7dw+9TUsLdHaGlvvB7pCaSXV1mEpyyRI46CA45phQ5s4d3flNEIrZkm+KSov6W9FTuDtdW7tCYr6xg44XQ+t5ohW948UOept7aX2ildYnWgfZO1AM5XPLKZtfRvn88lAWhMfkZcVVhdV1TkYvr/qIjxX1NxSRQSVa6hOlqSkMYt2+vX9Aa6IkljU0hER/MIsWwfHHwwc/CCedlJMq5nMf8bGimC3jyd3p3tE9MDlPStI7X+qka+vwfrgpmVYyMDmPkvW9yxaUUzqjVDO/TGITvY+4iMj4GUlLvXtI2DdvDvO9//Wv8NBDobz0Etx4YyivvRamhBSRvGZmlM0so2xmGVNeN2XQbfo6++jc0klnQyddDV10NoTnna8kLdvSSc+uHnp29dD6ZJqWdcDKjfJ5Ka3p80L3l0Qpn1tO8dRiJewFQIm4iEg2zPr7yx95JPzN34Tljz0GRxzRv92cObD//mF++OXL4bzz4A1viKfOIjIqReVFVC4J0ySm431Ry/ormRP2nqYeOjZ10LGpI/MxK4v6E/OkRL187sCkXS3sE5sScRGRXDjsMLjuOvjOd0L3laYmeP75UG6/Hb71rdBiPnv20PsSkQnHioyy+jLK6suYcuTgLesAva29oXU9OWF/tZOuV7v2ls5XO+lr7Qt92jdmTtit1NIm6eVzyymbEyXss0o14DQPKREXEclGe3vodrJtW3955RXYsCEk3a++Ovjgz5UrRz7rjIhMGsXVxVQtq6Jq2b4DS5P1tPTQtWXwJL1rS//z3ubevVM3tjD4rDAAGJTWl4YkfU5S0j6nP3FPLC+u1qDT8aJEXEQKS/IAzdR50Qd7nrqsI3PrFBBmVVm+PHRdOfJIOPpoWLWq4Kc8FJHhK5lSQsnyEqqWZ07Ye9t66dq6b5LetTUpid/aRXdjN92vhdL6WPo+7ADFNcUDkvPU1vW9rewzSrEixbXRUCIuIvnBPSTJ7e2jK21tYb7x5LnHk2dG6RrltNVlZWEQ5uzZ/WXePFi6FJYtC4+zZyvpFpFxUVxVTOX+4S6imfR199H9Wnd/gr61v5V9QCK/tYvePb20b2infUN7xn1aiVE6u3Rgkj5Y8j6njKJydYsZjBJxkULT2xuS0e7u8Jh43tkZSldX//PBXg9nm5HuYzymU02+iU9dXeguknhMfT7YuspKJdkiMuEUlRbtnUpxCun7sLs7PU09+yTqya3riceeXT10NXTR1TB0A0fJ9JIBLerpusaU1JYU1OBTJeIiI5VowR0sqc1m2Xiv7+uL+8qlV1ISEt3RlilToKZm4F0zE691W3sRkbTMjNJppZROK6V6ZXXGbXs7om4xWwe2rO+TvG/romdnDz07e2h7ui3jPosqijK3rkfLJsvgUyXikp/coadnYGtpogU1+THd89GuT30+WHLb0xP3VRoZs9C9oqwszKOdeCwvD6WsrP95umXD2WYky4omflAVESkUxRXFVC6upHJx5m4x3ud0b+/epwvMYF1j+lr76NjcQcfmIcbjGJTOKqViUQUV+1dQuX/lgMeKhRVYcf43uigRl4GSE+COjvSPmdbl6r353HKbUFIyeFKbuiyf1hdrNLyIiIwfKzLKZpVRNqsMDsu8bc+envSt60nLuhu76d4WSsvD+84WYyVGxeIKKg6ooPqgaqoPrqbq4CqqD6qmZGr+pL/5UxMJCfBok91cvCdfEuCSkn1bTBNJZS6eD3fbRGtxaqJbWqouDiIiIjlUUlNCybKSIad37Ovuo2tbF50vdtK+sZ2OjR39jy+00/VqF+3Pt9P+fDu7frdrwHvLF5RTc2QNtatrmbp6KlNWTaG4Ip5GKiXio9XbCzt3QmNjKDt3DpyhYffuzK9bW/uT4XxJgIsmttJMAAAUNklEQVSLQz/a8vKBj4MtG6ttysvVcisiIiKDKiotomJBBRULKqhdXbvP+t62Xjo2d9D2XBttT7XR+nQrrU+10ra+LdxM6ZVOdty2AwArN6afOZ1Z58+i/m/qx3WGFyXiQ3GHZ5+Fxx+HF14IN+zYuDHcxCOReOcqgS4qii/xTTwvLw8t0SIiIiITVHFVceiSclA1vKN/eV9PHx0vdLD74d3sfmA3zQ800/pkKztu3cGOW3ew8fMbOfA7BzLjrBnjUk9lXOls2wZXXgm//GV4nsm0aWFe4fp6mDEjzMwwderAGRvSLaup6U+ElQCLiIiIjJmikiKqlldRtbyKOe+bA0Dnlk4ab2pky3e30PZ0G0+87QkOve3QcUnGlfkNprUVVq8OLeAAc+bAG94QbtSxdCnsvz/MndufeJeWxltfERERERmR8nnlLLh4AfM/Np+Nl23k5atfZtPlm5SIx+a3vw1J+NKloUX88MM1KE9ERERkErNiY8kVS3j5ay+z56976OvsG/P+4pq0dzANDeHxjDPgiCOUhIuIiIgUgKLyIkpnhJ4O3Tu6x/54Y36EiShxoxZ1OREREREpKEUVIT3u6xr72eyUiA/GPTzqLn8iIiIiBcVKQ08I7/ExP5YyzcEkpiNUlxQRERGRgmIlUSLerUQ8HokWcSXiIiIiIgVlbyI+Di3iE37WFDOrBK4GWoBpwKXuvjtHO8/JbkREJBjTmC0ikgNFpaGdWl1Thuc7wP3ufhlwK/C9mOsjIiLpKWaLSF5T15RhMrN5wHuA26JFvwfONrPFcdVJREQGp5gtIhNCcXjwXiXiQzkR2O7ubQDu3gNsAk6Is1IiIjKoE1HMFpE8N559xCd6Ij4f2JGyrAWYF0NdREQkM8VsEcl7Vhwl4uPQIj7RB2s60J6yrBwYcCukxsZGVq1atff1mjVrWLNmzdjXTkRkBNauXcvatWsTL2fGWZccU8wWkbyXaBGnd3jbjyZmT/REvIEw6j7ZFGBL8oL6+nrWrVs3bpUSERmN5MTTzLbHXJ1cUswWkby3t0V8mF1TRhOzJ3rXlHuAedF0WJhZKbAI+GOstRIRkcEoZotI3hvPrikTOhF3963A7cAp0aJTgZvcvSG+WomIyGAUs0VkItg7WFN9xIflo8DVZnY4MAdQR0IRkfylmC0i+S0xfaHurDk0d98FXBh3PUREZGiK2SKS78azRXxCd00REREREcmlRB/x4c6aMhpKxEVEREREItnOmjIaSsRFRERERCLqmiIiIiIiEgNNXygiIiIiEodxnDVFibiIiIiISERdU0REREREYqBZU0REREREYrC3RVxdU0RERERExo8Ga4qIiIiIxECJeNyqq2HuXJg6Ne6aiIiIiMg4KqkroWxuGUWVY58ml4z5ESaij388FBEREREpKPt9YT/2+8J+43IstYiLiIiIiMRAibiIiIiISAyUiIuIiIiIxECJuIiIiIhIDJSIZ7B27dq4qxArnb/Ov9DpGkws+rwy0/XJTNcnM12fzEZ6fZSIZ1DoXzqdv86/0OkaTCz6vDLT9clM1yczXZ/MlIiLiIiIiEwg5j72dw2Km5k1Ai+O4K0zge05rs5EovPX+Rfy+UN+XIP93L0+5jqMK8XsMaPrk5muT2a6Ppklrk9WMbsgEnERERERkXyjrikiIiIiIjFQIi4iIiIiEgMl4iIiIiIiMSiJuwL5yMwqgauBFmAacKm77463VqNnZqcBVwLvcvfN0bK05zrSdfnIzN4JfA2YCtwIfNrde8xsJvDPQBNQSjiP7ug9I1qXr8zsBGAtMAv4ibt/MlpeEN+BBDMrAx4GPuHufyi085+M9FkMNJJ4V4iyiQWFysyOAY4DngfuB1rRNcLMVgIXEa7LMmCtu/91RN8hd1dJKcCPCckqwJnAT+OuUw7OaTZwDuDA4uGc60jX5VsBFgE3AEcB7ycEkkuidfcAb4iefwT4StL7RrQuHwtQC1wWBYazgG7g1EL5DqRciy8AzcCJhXj+k7HosxhwLUYU7wqxZBMLCrEAFwBXpizTNQrnvg6YHz1fBDwz0usT+8nkWwHmAZ1AVfS6JApki+OuWw7OrYikRDzTuY50XdznmOa83wSUJL2+Gvhv4BhgS9Ly+cAeoHqk6+I+1wzXoJJolqTo9SPASYXyHUg672OBDwKbgRML7fwnY9Fnsc/1yDrexV3nmK7TsGNB3HWN6fqcAPw+5f8NXaP+a9EKrIie1wNbRnp91Ed8XycC2929DcDde4BNhC/lhObufSmLTiT9uY50Xd5x9/uiOiY0AC8BJxPqndiuIXq6ahTr8pK7t3sUGcysGngG+AMF8h2Aved9rrv/KGnxiRTI+U9iJ6LPYq8RxruCMoJYUIi+Tvh/4ltm9lszW42uUbJfAT8wsynA+4CLGeH1USK+r/nAjpRlLYS/dCabTOc60nUTweuB71KA529mbwL+h9A/tIrCugaXAlelLCuk85+s9FlkNpx4V2iyjQUFxcyWA0cC33f3i4C7gN8CC9E1SvgY0EEYY7DH3W9ihN8hJeL7cqA9ZVk5oU/tZJPpXEe6Lq+Z2VKg0d0fowDPnzCw5HrgFOCrFMg1MLMzgHXu/lrKqoI4/0lOn0UaWcS7gjHCWFBoDgZ2uvvj0etvA4auUbJK4KfAfwLfMLNTGeH10awp+2ogDGhLNoXQ/2eyyXSuPsJ1ecvMioEPA5+LFjUAByatN0If8C2jWJfX3P1V4Edm1gd8BriPwvgOXAIcGT4qINT9VsLPr4Vw/pNZIcXsYcsy3hWSkcSCQlNCUn7o7u1mtgEoRtco4Ubg3e7eZGZOSMo/yUiuT9wd3vOtAHMIA1gqo9elhM738+OuW47OL3mwZtpzHem6uM9viHO/FJid9HoV8GLS6yWEn5HKRrou7nPM4locAtxRKN8BwmCaBUnlZeA8wmj3SX/+k7nos0h7XYYd7+Ku6zhfl6xjQdx1juEarSDkCsnfn4eB/6Nr5AAzgVeTXhvwAmGax6yvj7qmpHD3rcDthJ/uAU4FbvL+gS0TlvU3ARhkPteRrhuH0xgRM7sMeBSoNLP9zeyDhLl0XzWzQ6LNzgS+5e5d7r5uJOvG74yyY2ZVZnZk0qK3ANcWynfA3Rvd/ZVEAXoJP9m/RAGc/2Smz2Jf2ca7uOoZh5HEgrjqGhd3X08YS3QugJnVERqh/gNdI4CdQIeZzU9atgN4nBFcH4uydkliZtMIUz5tJrS2fN7d98RaqVFKGtn7HeAK4Nvu3pjpXEe6Lt+Y2eXAl1MWr3f3lWa2H/BFwnnUAF/waMaBka7LR2Z2NGEKs/XAn4C/uPsvo3WT/juQysw2Ax/wcBOPgjv/yUafRb+RxrtCNdxYUGiimz9dS5gvexFh4ObTukaBmR1OGLD5COE+Lfe6+x9Hcn2UiIuIiIiIxEBdU0REREREYqBEXEREREQkBkrERURERERioERcRERERCQGSsRFRERERGKgRFxEREREJAZKxEVEREREYqBEXEQmDDOrN7Nfm9mLZva5uOsjIiKZKW5nphv6iMiEYWZnAPcAUwl3CT3Q3XfGWysREUlHcTszJeIiMiGZ2b3AyYV+i24RkYlCcXtf6poiIhOOmS0E7lAwFxGZGBS3B6dEXMaVmX3WzK6Lux4JuaiPmZ1qZg+Z2Yk5qtakYmbX5Hh/BpwPXJXL/YrIQIrXo2dmi83s+2b2vRztr9jMLjSzF4ex7Vlm9ueRnmsuY7fidnpKxCcRM1tpZjeb2Q/M7DEzczN7X9z1SvEb4KdxVyJJLurzMnB4DuqSl8zsLWa23cxeNbM3mtnC6HvmZvZtM6uLtjvazJ40s6+ZWWW0bDawNcdVeg/wQ3fvifYvMuEoXo/IRIzXr0WP5TnaXxHQBCwaxrbrgJXpVprZwWb28+i7d7mZzUlal+vYrbidhhLxyeUW4J/c/UPAkcANiRVmdpiZnRBHpZKP7e5Puvv9cdRjMNnUJ901dPdnge25rFcuP6/R7svdfwN8B+hx94fc/WXgQqAL+LO7N0Xb/YUQ+D/j7u3R298G3D6qE0hiZv8A/CPwezNbDyzP1b5FxpnidZYmYrx29zZgS66O5e7dwCPD3HYbIWlPt/4pYG308hp3T068cxa7FbczUyI+SZjZdOBAoBXA3fuALwJuZrWEIG8x1Cu2Y+fSMM6jbxyPFce+fg4sMLOjANx9B3AX8NakY9UBW6LvXsIKd18/ymPv5e5fdvdl7r7K3Ve4+7252rfIeFG8Hlt5GK9zPStGNvUfatve6DG133bOYrfidmZKxCeJaCqgZ4Hfmtkx0bIXgb8Q+mUtAdaY2QUAZnaFmf1D9LPU3n5gZrbKzL5uZp80szYze9bM/p+ZvcPMHjazc8zsv8xst5l9OOl9g+4v5dgXmtn1yX3lop9nv2VmXzGze8zs6Gh5xuMlvd/M7MNm9pSZnWZmj5vZNjP7wDCOsSy5PkMcc59rOIhp0U/NLWZ2Yco1/cfo2twcBe7hXrN9jmVm55rZX8zsA2a2wcyuyrC/wT77QeuTSdRy8hRwXtLieuAsi7qhAGcDv06qZyXQlvR6hZldZWbXmNkfLXRxeZeZrTOzt5vZT82sycw+YWZHRee4w8zeMlT9RCYSxevCideDSXeOmY5jZkvM7EcWWpe/OpzjDHLc35jZf5jZ6iG2U+weT+6uMkkKcDDwIuEv4B8Bs5LWbQZOjJ4fAmyMntcT/lqfQfiL/iVgVbTuOuDB6HkZsBP4NmEu0L8FGjLtL/XYhD5y3wOuj5ZXEJK7uuj1uUAjMCXT8QY57znRMS+Ltr0S6AD2G+IYqfXJeMzkazhIHTYDXwOqove9FC2fCvwsabsHCEF0WNcszbGmE35u/AZwNHBqpv2lfPaD1meY36/LgRei5wcAXwE6gXcmfV8safuzgTdFz6uAh4Gq6PW90XWYCjQDVwCVwLui1+cTGgo+TfQdVFGZTAXF64KI19H6fxzmdcz0WT8MrEx6jw/ze5b4PN8IfGmQ9SdGxylJWqbYPY5FLeKTiIdWy0MJgen9wFNmNtiglA3A35lZKfCmaFkN4R/+QiDRv/cZwj8w3L0LaAF+5e67Cf8w5wyxv9T6dTJw8MdZQLtHfYyBmwn/gN85xPFS95vY539G215J6L982hDHGFCfbI6Zxm0e+gM+DMxNOscZUYvVJwnXFIZ5zdKc705CwLvV3f/i7r/PYn/p6jMcPwf2t9A95W+A7wO/A843sxqgxaNIHVkNPBg9P5vwn12ileUc4MvRdW4C7vbQr3wdMNXdf+Hh5/pHye4zEJkQFK8LI14PIu05ZjjOqUCtuyfqsy7LY74eeK+7/9Mwt1fsHkclcVdAcsPMioGF7r4Z+Hsz+yFhhPla4A3J27p7p5nNBy4FfpLYhbu/ZmbrgFMIf7EvISRfe9+a9LyHqGtTuv2lqWryPpYBpUn16jOzzcCCTMcbirt3mNlGoA6YmcUxRnzMFD30/9taBGxy92+kbjSca2Zmf0dIeBOWe/gJ25PrmsVnkLY+Q3H3DWb2KKF7Sq27bzSznxNaqc4D/jup3onvRqL/4X6EFqzEvnYk7zrpeS8D9aEudDLJKF4POL9CiNfJ0l7HDJ/NYYQ/OkbqZGCZmdW4+5402/RF56DYPc50kSaPUpL677r7Y4Sf/lakbmhmJwMfc/crPcyAkezjwOstTKP1OPAvQx14iP1l8hKhhTV5Wicj9J0crTLguTE+xnC8ysB+1JjZG7K4Zr8Gjkgqg46+z2J/g9Zn+KfDz4ALgI1J9SsGPkH4GTfhDYT+rgkNwHFmVp103Iz9FEUmMcXrgQopXqc9xwzH2R29p4yRuYbwa8K1adY3e/8ge8XucaZEfHK52MzmJb2eD/w2et5FGJyygjBVVq2ZlZvZqdH6WWY2k9BX8YeEn5UeJPqpM1JESiuAmdkQ+0s9tiXt4xZgD/DuaF910fESLavpjpfO3Gib+dF7fzuMYyTXZ6hjJp9HqsFaR4zQylUN3G5mZ1gYFDTUZ7D3WMBcd1+fVLqTDlGc9DzT/pLrPWh9zKzMzNZG3U4y+QXhJ/FfAbh7C/A/wAM+cLaU0wndVhIS1/unFuYivwhIzCWb+hmkfs4TegYHkTQUrymYeD3c6zjocQixtIzwxxqEXz8ws3nDiN1GaJ3+O+A8Mzs/eu+nzexywq8BdyVtr9g93sai47nK+BfCABAHdgE3AtcTftqaFq3/EvAK8HbCT00bCC0QZxFaUq4nBLVfEwaOdEX7ayP8w3wz0E0YELQQ+Hy0/oOZ9pdy7E8A9xH63R0crVtFaEn9GqGv5CHR8rTHS3P+Hr3/E4SuEocmrUt3jMXJ9RnqmMnXMOXYbycMWvx2dC0S77sgWn8C8GT02Xwtus7DvWZvH+Rcz42O91NgfrQs02c6YF9p6rOAcOOJ/xzGd+2bKa/fBZyUsuyaQd53PPA0sAO4PFr2FsJ37VvAPKIp3ICPEP6j/n50rm+N+9+YikquCorXhRSvFwH3E1r1jxjiHDN91mcRfon8E2Hg/HPAx6JzHzR2E3516Sa0hE8H7ib8EfAh4MOE+dS/RugWk3iPYvc4F4suqAhmNgv4iLtfEb0uIrTSfNjdL4+1ckMwMweWeOhzKSNkZp9x9xFNjSUi40fxWpIpdk9c6poiyS4idFOYnrTsZJIG4eWjpJ/C9DPYCJlZsZm9F7gp7rqIyLAoXoti9ySgRFySfZfQl+wZM3sZuBN4xd0firda6UWDVz4RvfyAhWn0JHs1wC3uvnHILUUkHyheCyh2T3jqmiIiIiIiEgO1iIuIiIiIxECJuIiIiIhIDJSIi4iIiIjEQIm4iIiIiEgMlIiLiIiIiMRAibiIiIiISAyUiIuIiIiIxECJuIiIiIhIDJSIi4iIiIjE4P8HpLBywrW8oZ4AAAAASUVORK5CYII=\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/igpd-uranus.png',bbox_inches='tight')\n", "plt.savefig('../plots/igpd-uranus.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/igpd-uranus.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Discrepancies are attributed to initial state being not properly converted to planet-relative frame, and different heating relations used in the studies." ] } ], "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 }