{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 44 - Ice Giant Pre-Decadal - Neptune 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 Neptune 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(\"NEPTUNE\") \n", "planet.h_skip = 1000.0E3\n", "planet.h_trap = -100.0E3\n", "planet.loadAtmosphereModel('../atmdata/Neptune/neptune-gram-avg.dat', 0 , 7 ,6, 5 , heightInKmFlag=True)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle1=Vehicle('igpd-neptune', 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(1000.0,0.0,0.0,25.73,0.0,-20.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.0*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+17YcXAAAgAElEQVR4nOzdeZhcdZX/8ffpNZ0NsjSEBEIjIsQomw0i+hMUUFyIwBBHQBH9SSPuM8qPEdFRwAVUxpWlBZxRdERgRgZUGFBkBlEkkYCigEAWCCR09pB00unu8/vje6tTXemq7upa7r1Vn9fz3Keq7rfq1rnV9Zw6/a1zb5m7IyIiIiIildMQdwAiIiIiIrVORbeIiIiISIWp6BYRERERqTAV3SIiIiIiFaaiW0RERESkwlR0i4iIiIhUWFPcAVTDzJkzvaOjY8z37+npob29vXIBxaQW96sW9wlqc79qcZ+g8vu1ePHiNe5eey9cAcXm7Iy0vscUd3Up7uqqt7gL5ey6KLo7OjpYtGjRmO/f2dlZ1P3Tohb3qxb3CWpzv2pxn6Dy+2Vmyyu28YQqNmdnpPU9prirS3FXV73FXShnq71ERERERKTCVHSLiIiIiFSYiu4RdHV1xR1CRdTiftXiPkFt7lct7hPU7n6lUVr/Foq7uhR3dSnunczdy77Rgk9odgJwKfD37r4sWtcGXA5sBqYBF7j7plLGsnV2dnoa+4lERMxssbt3xvj8ytkiImNUKGdXdabbzPYEJgNH5gxdCdzn7hcCtwLXlGFMRERKoJwtIlI+ccx0NwADwH7uvszMZgNLgWnuvtXMmoCNwHygbzxjmdmYDM2aiEgcBgcHGRgYoLm5edzbSMBMd+JztjsMDobrjY2l7K2ISGkSM9MN4O6DOauOBda4+9ZovJ+QmI8pYUxEJFY//vGPaWxs5L3vfW/coZQkDTn7/vuhqQle//pStyQiUjlJOE/3HGBtzrrNwGygf5xjw/T09NDZufOfjq6urtQ29otIOrS0tADQ19dX9GO7u7vp7u7O3JxZvqjKQjlbRCTLWHN2EopuB3pz1rUCO6Lr4x0b0t7ensoTs4tIemWK7u3btxf92Owi08zWlDWw0ilni4hkGWvOTkLRvZJwFHu2KcBzhOQ+njERkViVMtOdcInN2VU+RElEpChJOE/3PcDs6FRSmFkzMBe4t4QxEZFYtba2AjVZdCcuZ5uVugURkcqretFtNpQeDcDdVwG3A8dF648HbnH3leMdq8JuiIgUVEp7SZIoZ4uIlEdV20vMbArw7ujm2Wb2HXfvAc4DLjezQ4BZQPYRM+MdExGJTS3MdKctZ6u9RESSrKpFt7tvBq6Kluz164Fz8jxmXGMiInGqhZnutORstZeISBokoadbRKTm1PCBlCIiMg4qukVEKqAW2kvSRu0lIpJkKrpFRCqgFtpL0kLtJSKSBiq6RUQqQO0lIiKSTUW3iEgFqL2k+tReIiJJpqJbRKQC1F5SPWovEZE0UNEtIlIB2e0lrilYEZG6p6JbRKQCGhoaaGoKP4XQ398fczQiIhI3Fd0iIhWiFpPq0hcKIpJkKrpFRCpEZzCpDvV0i0gaqOgWEakQncFEREQyVHSLiFSI2kuqS+0lIpJkKrpFRCpE7SXVofYSEUkDFd0iIhWi9hIREclQ0S0iUiFqL6kutZeISJKp6BYRqRDNdFeH2ktEJA1UdIuIVIhmukVEJENFt4hIhehAyupSe4mIJJmKbhGRClF7SXWovURE0kBFt4hIhai9REREMlR0i4hUiNpLqkvtJSKSZE1xB5BhZvOAjwBPAgcA3e6+xMzagMuBzcA04AJ33xQ9Ju+YiEjcarm9JEk5W+0lIpIGiSm6gR8C73D3lWY2F7gTmAdcCdzh7jea2YnANcDp0WMKjYmIxKrG20uUs0VEipCk9pJ5wJToei+wm5nNBs4AbovW3w0sMLOOQmNVi1hEpIAaby9RzhYRKUKSiu6bgWvNbArwbuCjwLHAGnffCuDu/cBS4JhRxkREYlfL7SUkMGerp1tEkixJRfeHgG3Ag8CL7n4LMAdYm3O/zcDsUcZERGJX4+0licnZ6ukWkTRIUk93G/DvhMT8DTNbCjjha8tsrcCO6HqhsSE9PT10dnYO3e7q6qKrq6tMYYuIjGy87SXd3d10d3dnbs4sb1Rlo5wtIsLYc3aSiu4bgHe5+wYzc0Iy/wThCPdsU4DnCMk939gw7e3tLFq0qPwRi4gUMN72kuwi08zWlD2w8khczlZ7iYjEYaw5OxHtJWY2EzjE3TdEqy4FNgHLgdnRaaYws2ZgLnAvcE+BMRGR2NVqe0nScrbaS0QkDRJRdAPrgG1mNidr3VrgEeB24Lho3fHALe6+0t1X5RurUswiIgXV8IGUytkiIkVKRHuJuw+a2cnA58xsMbAncL67bzKz84DLzewQYBaQ3dhXaExEJFa1esrApOZstZeISJIlougGcPeHgXNHWL8eOCfPY/KOiYjErVbbSyBZOVvtJSKSBklpLxERqTk13F4iIiJFUtEtIlIhtTzTnURqLxGRJFPRLSJSIbXa0500ai8RkTRQ0S0iUiFqLxERkQwV3SIiFaL2kupSe4mIJJmKbhGRClF7SXWovURE0kBFt4hIhai9REREMlR0i4hUiNpLREQkQ0W3iEiFqL2kutTTLSJJpqJbRKRC1F5SHerpFpE0UNEtIlIhai8REZEMFd0iIhWi9pLqUnuJiCSZim4RkQpRe0l1qL1ERNJARbeISIWovURERDJUdIuIVIhmuqtL7SUikmQqukVEKqSxsZGGhgYGBwcZGBiIO5yapfYSEUmDshTdZvbWcmxHRKTWJLHFRDlbRKT6moq5s5mdCFwEzGJnwd4A7Am0lTc0EZH0a21tZdu2bfT19TFx4sSqPne95Wy1l4hIkhVVdAM/BL4MPAQMRusMOKOcQYmI1IqYZ7rrImervURE0qDYovtpd78id6WZPVameEREakrM5+pWzhYRSYhie7rPM7PTRli/oBzBiIjUmpjPYFJXOVvtJSKSZMUW3a8FrjWzgaxlELiqArGJiKRezO0ldZGz1V4iImlQbHvJxUAXsAjInP/KgPeUKyAzOwp4HfAkcB+wBbgc2AxMAy5w903RfdvyjYmIJEHM7SXK2SIiCVFs0f1L4FZ3HzZlY2bfK0cwZvY+YH93vyhr3feBO9z9xuhI/GuA06PhKwuMiYjELub2EuVsEZGEKLbo/gnwaTP7ddY6A04FPl5KIGZ2DHAmcELWutmEo+w/HK26G7jFzDqAvnxj7r6slFhERMol5vaSusrZ6ukWkSQrtuj+LHAAcDaQSW9GOAdsSQkcuAK4H/i2mb0UuATYF1jj7lsB3L3fzJYCxwA7CowtKzEWEZGyiLm9pC5ytnq6RSQNii26vwr80t03Z68s9dfNzOxA4DDgfe7+iJmdD9wBXAqszbn7ZmA20F9gbJienh46OzuHbnd1ddHV1VVKyCIiYzKe9pLu7m66u7szN2eW8PTK2SIiFTbWnF1s0f1MbvKOzCpyO7nmA+vc/ZHo9neAfybMzPTm3LeVMGPCKGND2tvbWbRoUYkhiogUbzztJdlFppmtKeHp6ypnq71EROIw1pxd7CkDP5S7wswWAt8ocju5msj6B8Dde4G/AY2EI9yzTQGeA1YWGBMRSYSYD6Ssi5yt9hIRSYNii+7pZvZZADObaGbXAtcCD5YYxyPAbma2Z9a6fkIynh2dZgozawbmAvcC9xQYExFJhJh7upWzRUQSoqii293fBtxlZt8GlgAvAw4BRvrFs2K2+xjh1FanAZjZ7kAL8CPgduC46K7HA7e4+0p3X5VvrJRYRETKKc6zl9RbzlZ7iYgkWcGebjN7fZ6h5YSvEr9GOFr9ZOAfSozlLOCbZtZCmP04PTq6/TzgcjM7hNCHmH00TaExEZHYVbO9pF5zttpLRCQNRjuQ8nJgHrABGBxh/HrC6af2osQE7u5rCOd8zV2/Hjgnz2PyjomIJEGV20uUs0VEEmq0ovtiYJG7v1DoTmZW0leVIiK1qsrtJXWds9VeIiJJVrDodvdfjGUj7n5zecIREakt1WwvqdecrfYSEUmDYs9eIiIiRYj57CUiIpIQKrpFRCoozrOX1Bu1l4hIkqnoFhGpoJh/HEdERBKi5KLbzA4zsxvM7JNm1lF6SCIitSNpM921mLPV0y0iaVCOme4PAse7+9eBozO/NiYiIons6VbOFhGJQTmK7o8Ar4qu/wRoLsM2RURqQgLbS2o2Z6unW0SSrOii28zeYmZnRtePBA7I/Iyvuw+6+6Yyxygiklpxt5fUQ85We4mIpEFRRbeZfRG4CTgFwN3/AJxlZmdUIDYRkdSLs71EOVtEJDmKnek+Gtgb+N+sdb8EvlK2iEREakjM7SV1lbPVXiIiSVZs0X2Xu28AslNbTf6csIhIOcTcXlIXOVvtJSKSBgV/Bn4ET5vZhcAcM3sXsJDwteVHyx6ZiEgNiHmmWzlbRCQhiprpdvefAL8GGoGzgA3Am9z9uxWITUQk9eLs6a63nK32EhFJsmJnunH33wO/z9w2s0YzO87df1XWyEREakDcZy+ph5yt9hIRSYOCRbeZ/RcweZRtzATWADWTwEVEyqWa7SXK2SIiyTXaTPcj0X0eBwaBc4CHgEVZ99kPmF6R6EREUq7K7SV1nbPVXiIiSTZa0f1VwKKj3zGz/d39c7l3MrNbKxGciEjaVbm9pC5zttpLRCQNChbd7r4xZ9WeZtbo7gOZFWZ2LHBkBWITEUm9araXKGeLiCRXsQdS3gw8bGb3AZuBecCbgV1mUkREJN6zl6CcLSKSGMWeMvAu4E2Eg3AOJJx+6jR3/3IFYhMRSb04z15SLzm7sTFc9vfHG4eISCHjOWXgc8BFmdtm1mBmr3L3xWWNTESkBsT84zh1kbOjl5iYXmIRkTEpqug2s+tHWN0ODAAnlxqMmbUADwIfd/ffmFkbcDnha9FpwAXuvim6b94xEZGkaGoKaba/v5/BwUEaGor6grEklc7Z0XPEnrejLxOI6VToIiJjUmz2PwBYASzPWlqBcs2YnA90ZN2+ErjP3S8EbgWuGeOYiEgimFmcfd2VztmQgLytmW4RSYNi20s+4O6PZ68ws4OAs0sNxMyOBp4H1ke3ZwNnAB+O7nI3cIuZdQB9+cbcfVmpsYiIlFNrayt9fX309fUxYcKEaj51xXJ2tK1E5O3smW53nUJQRJKp2AMpHx9h9XbgA6UEYWaTCAf3ZH8Veiywxt23Rs/dDywFjhllTEQkUeI6mLJSORuSlbcbGiDq4tHBlCKSWMX2dC8Fsn/zqxHYE/hhiXFcAHwlZ90cYG3Ous3AbKC/wNguenp66OzsHLrd1dVFV1dXKfGKiIxZse0l3d3ddHd3Z27OHO/zVjBnQwXz9nhydmtrKLi3b4fm5jHFLyJSFmPN2cW2l9wF/JidSXwQWO3uTxQdYcTM3gwscvcXcoYc6M1Z1wrsiK4XGhumvb2dRYsWjTQkIlJxxZ7BJLvINLM1JTx12XN2FFNF8/Z4cnZLC2zZor5uEam+sebsYovuG9z9f7JXmNmBZvZSd3+y+DAB+BRwmO1swptGOMDmiuh6tinAc4TEnm9MRCRRYjxXdyVyNiQwb2cOptQZTEQkqYo9e8mJI6xbDdxUQgxnAIdmLc8R+g2/D8yOTjGFmTUDc4F7gXsKjImIJEqM5+quRM6GBOZtnTZQRJJuTEW3mX3EzDYAF5jZQPZC6NFbPd4A3L3H3Z/NLITzx/a4+wrgduC46K7HA7e4+0p3X5VvbLxxiIhUSrVPGVjJnA3JzNtTp4bLTfq1BhFJqDG1l7j7d8zsf4D3At/KGe5z9+fLHllwHnC5mR0CzAK6xjgmIpIY2T+QUw0x5myIKW9PixpX1q8v1xZFRMprzD3d7v6Imf2Tu+9y0IuZzY1mOErm7h1Z19cD5+S5X94xEZEkGRwcBKjqr1FWK2dHz9WRdT2WvL377uFSRbeIJFXBotvMPgX8wt3/Eq36f2bWmHs34A2Ec7CKiEiOahXd9ZyzMzPdGzbEG4eISD6jzXTvB+yedftwYBKwip2noDLC14QiIjKCKs50123O1ky3iCRdwaLb3T+cs+qLwJ/dfdjRQGb28nIHJiJSKwYGBoDKF931nLNnRj9H8ULumcNFRBKi2E+AV46QvA8E9i9fSCIitSWOnu5I3eTsffYJl888E28cIiL5jOlASjP7P4SvJE+IflY423TC0fG3lTk2EZGakCm6Gxtz26srox5z9ty54XJF2Q4PFREpr7GevWQD8ANCv+Brc8b6gGvLGZSISC3JFN1Zv+BYaXWXs1V0i0jSjfU83X8ys2OBhe5ec8laRKSSent7AZgwYUJVnq8ec/Y++0BjY2gv2boVJk6MOyIRkeHG3GDo7hvzJW8ze135QhIRqS2bN28GYMqUKVV7znrL2a2tcNBBMDgIf/pT3NGIiOxqtPN0/xcweZRtGDAf2KNcQYmI1JIXX3wRqHzRXe85+7DD4NFHYckSePWr445GRGS40dpLHif0Bq5k5zlecxnw9nIGJSJSK/r6+ujr66OpqYnW1tZKP11d5+xDD4UbboAHH4Rzz407GhGR4UYrui8DNrv79kJ3MrOaOgpeRKRcsltLqnAgZV3n7GOOCZd33w3uUL3jVkVERlewp9vd14yWvCOnlCkeEZGasmrVKgD22KPy3Rz1nrMPOwxmzIDly+GJJ+KORkRkuJJ+qcHM5prZZcA3yxSPiEhNWblyJQBz5syJOZLaz9mNjXDiieH6T38abywiIrmKLrrNrNHMTjGzXwBPAWcSeghFRCRH3EV3veXs97wnXP7rv4YzmYiIJMWYi24z28/Mvgg8A/wEeClwErAPsLAy4YmIpNvy5csB2Hvvvav6vPWas48/Ppyz++mn4fbb445GRGSnUYtuM1toZv8N/A34O+AKYG/g5+5+hwf3VzhOEZFU+lN00uj58+dX5fnqPWc3NsInPxmuf+ELmu0WkeQYy0x3L+Fng+8HFrj719y9h/ynoxIRkcjDDz8MwMEHH1ytp6z7nN3VBbNmwR//CNdfH3c0IiLBqEW3u9/u7m8H3g28x8x+aGanAc2Z+1gVzoMlIpI269at46mnnqKlpYUDDzywKs+pnA1tbfD1r4fr558Pzz4bbzwiIlDcz8CvcPfPAu+LVu1tZpeYWQdwTgViExFJtV/96lcAvPa1r6WlpaWqz13vOfv00+Gtb4UNG+CUU6C3N+6IRKTeFX32Enfvd/eb3f0U4N+AfwC+UfbIRERS7s477wTghBNOiC2Ges3ZZvBv/wYdHbBoEZx6KmzbFndUIlLPSjpPt7s/6e4fBz5YpnhERGpCb28vN998MwALFiyIOZqg3nL2zJnhDCYzZ8Idd8BJJ8H69XFHJSL1qqSiO8Pdf1DqNszsVDNbamZrzeybZtYUrZ9pZteY2WVmdoWZZfcl5h0TEYnTzTffzMaNGzniiCOqduaSsaqnnD1/PtxzD7S3h5+Hf/Wr4S9/qfSziojsqixFd6nMbC5wMnAa4avPDwCfiIZvAq539wuAJ4CLsx5aaExEJBYDAwN86UtfAuDcc8+NOZryS1vOfsUr4A9/gIMPhr/9DQ4/PBxoOTBQjWcXEQkSUXQD+wLvd/fF0QzMd4E3mNlRwIHu/kB0v9uAj5rZpEJjVY9eRCTL9ddfz2OPPUZHRwdnnXVW3OFUQupydkcH3H8/nH02bN8On/pUmPX+5S/B6+ZkiiISp0QU3e7+v+7en7VqJbACeCOwNOt+K6OrnaOMiYjE4tlnn+VTn/oUAF/60pdobq69rre05uxJk+D734fbboPZs2Hx4nCGkyOOgBtugL6+akUiIvUoEUX3CI4ArgbmAGtzxjYDs0cZExGput7eXk499VQ2bdrESSedxLve9a64Q6qWVOXst789tJl89avhIMvFi+E974F994WLLoLHHqt2RCJSDxJXdJvZS4Eed3+Y8AtquWdXbQV2jDI2TE9PD52dnUNLd3d3BSIXkXq2bds2Fi5cyIMPPkhHRwfXXXcd4/0Nmu7u7qF8Bcwsa6BlltacPXFiaDFZsQKuvTb0fa9aBV/8IsybF2a/v/ENWL687E8tIjVmrDnbPEHNbGbWCHwFuMjdt5vZp4E3uPubonEDtgFvAI7JN+bu92dvt7Oz0xctWlTFPRGRerJ27VoWLlzIPffcw/Tp07n33nt5xSteUZZtm9lid09k21wt5Wx3uPde+OEP4eabYdOmnWOvfGWYHT/pJDjySGhsrGpoIpIihXJ20ma6PwV8zd23R7fvArJ/O7kD6AMWjTImIlIVDzzwAIcffjj33HMPs2bNKmvBnQI1k7PN4Nhj4brrwoz3T38Kp50GU6bAn/4EX/4yHH00zJgBCxbAFVfAQw/pDCgiMnaJKbrN7ELgIaDNzF5iZu8HNgDPm1nmE+xE4Nvu3ufui/KNVT14Eak7mzdv5hOf+ARHH300K1as4Mgjj+T3v/993RTctZyz29pg4UK46SZYswbuugs+9jF4yUtg48ZwIOYnPxlOPThzZijCv/AF+K//gmee0dlQRGRkiWgvMbOLgEtyVj/m7vPMbF/gc8AyYDLwmcxR84XGsqm9RETKZevWrVx99dVcdtllvPDCCzQ2NvKP//iPXHLJJbS2tpb9+ZLYXlLPOXv5cvjNb8Jyzz0j93xPmQIHHrjrcsABoZdcRGpXoZydiKK70pKcwEUkHZYtW0Z3dzfXXnstPT09ABx11FFcddVVHHrooRV73iQW3ZWWppy9bBn89reh1eShh2DJEli3Lv/9586F/fYbedlrL2hIzPfPIjIehXJ2U7WDERFJi7Vr1/Kf//mf3HTTTdx1111kJile9apXcfHFF/OWt7xl3GcokdrQ0RGWM88Mt91DS8rjj++6PPVUOFvKihXhoM1cra3htIX77gt77x2WOXN2Xt97b5g+PfSfi0j6qOgWEYm4O48++ij//d//zR133MGvf/1rBqIj5VpaWli4cCHnnXceRx99tIptGZEZtLeH5XWvGz62Y0doR1m6dNdl2TJ44QV44omw5DNhws5CfM6cMDs+a9auy/TpmjUXSRoV3SJStwYHB/nLX/7C/fffz3333cfdd9/N888/PzTe2NjIm9/8ZhYuXMjJJ5/MjBkzYoxW0q65GV760rCMZMuWUHyvWAHPPhuWlSuHX9+wIcyYP/VU4edqaoI999xZhGdfz10mT9bsuUg1qOgWkbrg7ixfvpwlS5bw0EMP8cADD/D73/+ejRs3DrvfXnvtxQknnMCb3vQmTjzxRBXaUjWTJsH8+WHJ58UXdxbiK1eG0xuOtKxfH8ZXrhz9eSdO3Dk7P9oycyZMnaoiXWQ8VHSLSM3ZsGEDjz32GH/96195+OGHWbJkCUuWLNmlwAaYO3cuRx99NK95zWt44xvfyPz589U6Iok1efLOs6EUsn07rF49vBDPvb1qFTz/PGzdGtpexvrrmy0tofgeqSifMSO0tkyfvvP6jBmaTRcBFd0iklLbt29n+fLlPP300zz22GPDltWrV4/4mPb2dg477DAOPfRQOjs7Ofroo5kzZ06VIxepvNbWcKaUuXML388dNm+Gnp6dy5o1w2/nLlu2wHPPhWWsmpt3FuO5BflIl5nrEyeqWJfaoaJbRBJpx44dPP/886xYsYKlS5fy9NNPs3Tp0qHrK1euJN8pT9va2jjwwAM56KCDOPjggzn00EM59NBDmTVrlmaxRbKYhXaRqVNh//3H9pje3vzF+bp1YVm7dvjlli1hpj3P/8N5tbSEAny33UZedt89/1hmaW4u/nURqQQV3SJSdVu2bOG5557j2WefZeXKlTz77LO7XF+9enXeohqgoaGBfffdl/3222+owM4s++yzDw06dYNIRbS1jW0WPdv27cOL8JEK85Eut23b2QpTSryjFeajFfAq3KUcVHSLSMncnY0bN/LCCy+wevVqVq9ePex67rJly5ZRt2lmzJ49m7333puXvOQl7Lfffuy3335D1/fZZx+a9UkokgqtreH0hnvtVdzjentDAb5xY/5lw4bC4729YSm1cJ86NRzsOnFiWNraRr4+0tiECeFytKW5We00tUxFt4gM6evrY+PGjWzcuJH169ezbt26oWXt2rV5b69fv37ofNZj0dLSMlRQZ5Y5c+YMuz1r1iyampSiROpZW1s4H/l4D71wD60t4y3Ycwv3SmtoGHuBXkwxX+i+SrPVo5dapAbs2LGDLVu28OKLL/Liiy+yadOmoeI5s4y0Lnds27Zt445hypQptLe3s+eeew5b9thjj13WTZ06Vb3VIlJxZuHMKZMnl6dw7+0NZ3vJXI7leqZgH2nZtm347R07dj6+WpqaRi/YMzP2hdZNmBCW1tbCl5nrTU31N6uvolukCgYGBujt7R1atm7dOux2Zt2LL744rHge7Xrmsq+vryxxNjU1sdtuu7Hbbrsxbdo0pk+fzvTp05kxY8bQ9ZFuT5s2Ta0eIlKTsgv3ShsYGHuBXkwxX2jp7w9nsNm8ufL7l62hYfQCPV/BXs7HVLP4V9EtdcPd6e/vZ/v27Wzbtm3oMvt6obFCxfJoBXW5iuJ8GhsbmTx5MpMmTWLSpElDhXNmmTp16i7rRhpva2vTDLSISEwaG6tX4EOYxe/vH7kYz52lL3R769ZwsOy2bWO/zP4HI2laW+FDH4JLLw0z+uWioluqYnBwcFghm1vYjrXwLbZQzr0cHByMZf/NjLa2NiZOnEhbW9vQknt7ypQpTJo0aaiAnjx58rDrI62bPHkyLS0tKpZFRGRU7qHgHRgIBXfmemNjmP1tbg4HjOaO594uNDaWx/b1hWJ9y5bwS6uZ65nlxRd3Xfr7q/Mabd8O//Iv8KpXwZlnlm+7KrrrQG7Bm1vsFiqAx3s7d92OHTvifhmAMCM8YcIEJkyYQGtr67DL0dYVKpZHK6hVFIuIjI87DA4OL+DGuuQWfuVeKrn9cm07dzsxzT0lTqYPPbu9ZMeOUHD39cGFF8I73lHe51TRHZO+vj62bNky5mXr1q1s2bKF3t7eogvgpBS8uXHuVPcAACAASURBVMVs5nZrayttbW1jLn7Hu661tVVnwxCRqhscDB/imWXHjmQWZ0ktLFUkll9DQ5jZbmoKl7nXS71djsc2NYVZ9+bm4ddzb4/1fs3NOwvtuE7NqAqkBNu2bRv6MY+enh7Wr1/P+vXr2bBhw9D1zO0NGzYMHfi2ZcsW+qv1HUkkt9Ct1u3Mdc30ikgS9PWFX07csGH419j5rm/bNrxgziyZ2bCxrK9yuq9JZsOLtGKW3AKv3Eslt1/qtkcqbhsa6u+sIUmhonsM3J2HH36Ye++9l8WLF/PnP/+ZZ555hjVr1ox7m01NTUMHvY11mThxIhMnTiy6CG5ublbBKyJ1wx0WLYL77oNHHoEnn4QXXgjLhg3xxNTSsnNpbk5mgZbkwlUfYVILVHSP4ne/+x3nnXceDz/88C5jTU1NQz/oscceezBt2rShZffdd9/ldvbZJVpaWmLYGxGR2vbss3DKKaHoHkljI7S3h5/8njw5HDCWucws2bfb2oYXzK2tw2/nW5e9vh7PRywiu1LRXcCSJUs47rjj6O3tZebMmSxYsIAjjjiCQw45hI6ODvbYYw8aGxvjDlNERCKf/nQouNvb4eST4ZBDYN48mDUL9tgDpk8PX6+LiFSbiu483J0Pf/jD9Pb2cuaZZ3LdddfR2toad1giIlLAn/8cLm+9FV7zmnhjERHJpv/38/j1r3/N/fffT3t7O1deeaUKbhGRFOjpCZfj/clvEZFKSf1Mt5m1AZcDm4FpwAXuvqnU7V599dUAfOxjH2Pq1Kmlbk5ERKhczs5YuzZczphRri2KiJRHLcx0Xwnc5+4XArcC15S6wf7+fu68804AzjrrrFI3JyIiO5U9Z2f09YVT/DU1lfenm0VEyiHVRbeZzQbOAG6LVt0NLDCzjlK2++ijj7J582b2339/5s6dW1qQIiICVC5nZ2zcGC53201nCxGR5El10Q0cC6xx960A7t4PLAWOKWWjTz75JAAvf/nLSwxPRESyHEsFcnaGe7hUwS0iSZT2nu45wNqcdZuB2dkrenp66OzsHLrd1dVFV1dX3o2uWLECgI6OjjKFKSIydt3d3XR3d2duzowzljKrSM7OmDQpXG7aFH4ZUse/i0g1jDVnp73odqA3Z10rsCN7RXt7O4vy/VLCCPr6+gBoa2srMTwRkeJlF5lmNv6fvk2eiuTsjEmT4NBDYckSuO02OO208QcqIjJWY83ZaW8vWUk4+j3bFOC5GGIREZHCKp6z3/e+cHnFFTvbTUREkiDtRfc9wOzoFFSYWTMwF7g31qhERGQkFc/Z739/OF3g734Hd91Vrq2KiJQu1UW3u68CbgeOi1YdD9zi7ivji0pEREZSjZw9eTKcf364/vnPa7ZbRJIj1UV35DzgHWb2GeCtwOhH24iISFwqnrM//GGYOTPMdt90U7m3LiIyPmk/kBJ3Xw+cE3ccIiIyumrk7MmT4YtfhHPPhY9/HE44AabldpKLiFRZLcx0i4iIDPOBD8DRR8OqVfDe98LgYNwRiUi9U9EtIiI1p6EBbrgBdt89nD7wq1+NOyIRqXcqukVEpCbttx/84Afh+oUXwu23xxuPiNQ3Fd0iIlKzTjoJLrootJcsXAi/+U3cEYlIvVLRLSIiNe3ii+GDH4Rt2+Dtb4ef/zzuiESkHqnoFhGRmmYG3/0unH02bNkCCxbAt76lc3iLSHWl/pSBIiJS3wYG4LOfhaefDtcHB3deZl8fGAj3HxwMpxJ85BG49tp4YxeR+qGiW0REUu2Pf4Qvf7n4x91+eyjAG/Sdr4hUgYpuERFJtb6+cHnAAXDppaGIbmwcfpm7zgzmz1fBLSLVo6JbRERqQns7vPOdcUchIjIyFd0jOOCAA1iwYAHz5s2LOxQRERnF9Onh4MiDDoo7EhGR/FR0j+DUU0/l1FNPjTsMEREZg3nz4NZb445CRKQwdbOJiIiIiFSYim4RERERkQpT0S0iIiIiUmEqukVEREREKkxF9wi6u7vjDqEianG/anGfoDb3qxb3CWp3v9IorX8LxV1diru6FPdOKrpHkNY3yGhqcb9qcZ+gNverFvcJane/0iitfwvFXV2Ku7oU904qukVEREREKszcPe4YKs7MeoDlRTxkJrCmQuHEqRb3qxb3CWpzv2pxn6Dy+7Wvu7dXcPuJM46cnZHW95jiri7FXV31FnfenF0XRbeIiIiISJzUXiIiIiIiUmEqukVEREREKkxFdxYzazOzb5vZl8zsKjObGndM5WJml5qZR8vDccczXmZ2gpk9YGYdWetS/3cbab+i9an8u5nZqWa21MzWmtk3zawpWj/TzK4xs8vM7Aoza4471rHKt0/R2A1Zf6db44yzHqUlB6T5PTRSLkr66x7lVM9Z3hqNJS63Fvv5lpTXP0/ced/r0Xjs7/diP3fL8Xo3jX6XunIlcIe732hmJwLXAKfHHFPJzGwS4YCAE6JV4zlAKXZmticwGTgyZyjVf7d8+5XWv5uZzQVOBk4D5gNXAc8AXwNuAv7J3R8wsw8CFwOfjivWsSq0T9HYC+z8Oz0eS5D1LfE5IM3voQK5KLGvu5nNA34HfArYEa2+Fvh1EnPrOD/fYn/9R4p7lM8AkvB+H+fnbsmvtw6kjJjZbGApMM3dt0b/lW0E5rv7sliDK5GZfSK6erW7b4s1mBKZWQMwAOzn7stq5e+Wu1/RulT+3czs/wC/c/f+6PblhMR7CfAf7j47Wj+HkGz3dPctccU7Fvn2yd3fZmb/AvwW+FlmXKonLTkgze+hkXJR0l93M5vr7iuybu8LXO7uf5/U3FrM5xvQl2+s2q//CHHnfa9HtxPxfi/mc7dc73e1l+x0LLDG3bcCRG+EpcAxcQZVKjMz4F3A14HVZpaIWYjxcvfBnFXHUgN/t9z9SvPfzd3/NyeRrgRWAG8k/G0y91sZXe2sYnjjkm+fzGwicBLw0+j2cbEEWN+OJQU5IK3voQK56FgS/LpnF9yRU4CfJTm3Fvn5VmisqnLjLvAZQJLe70V+7h5LGV5vFd07zQHW5qzbDMyOIZay8eAowtcl3wZ+ZGZviTmsctLfLfmOAK6mtv5WRxBmQra6+0sJ+/Zz4Jdm9sp4Q6s7aX1fpeI9VCAXpe11fyvw85Tl1kKvcZpe/8xnAEl+v4/y3ijL662e7p0c6M1Z18rOXrBUc/f1wEXRf3KfAH4Zc0jlor9bgpnZS4Eed3/YzGrib5W9T5l17v48cE40i/MR4Ny44qtDqXtfpfE9NEIuuouUvO5mNhPY5u6bMutSkltHe28n/vUf6b0OyX6/53lvlCXPaKZ7p5XAtJx1U4DnYoilkr4FzI07iDLS3y2hzKyRkEj/KVo17G8VJbRJpOhvNcI+5Urd36kGpCoH1MB7KBNfml73BcB/5RlL8utd6DVO/Os/hvc6JPv1z46tLK+3iu6d7gFmm1kbgIVTmc0F7o01qvIbBP4YdxBlpL9bcn0K+Jq7b49u3wUcmDXeQTgYaFGV4ypF7j7lSuPfKe3SlgPS/h7KxJem171Q0Z3k17vQa5yG13+09zok+/XPjq0sr7eK7oi7rwJuBzJN/ccDt2Qd7JVKZjbbzM4ws4ZoZvGTwGfijmu8on0AMKidv1vufqX972ZmFwIPAW1m9hIzez+wAXjezF4R3e1E4Nvu3hdXnMUYaZ/M7BVmtiAabwHeC1weZ5z1Jk05II3voXy5KC2vu4VTwE1w9xei24nNrcV8viXp9c+NO1o30nv9pWZ2UFLe78V87pbt9XZ3LdFC+Orge9GL/G1gctwxlWGf5hOOGv4z4UCGg+KOqYR9mQKcR+it+gLQXgt/t5H2K81/N+CiaF+yl79GY/sC1wGfBS4DmuKOt5R9ihLwC8ADwDeB2XHHWo9LGnJAWt9DhXJRSl7304BPjGV/Yo6z6M+3JLz+eT6/Cn0GJOL9Pp7P3XK83jpPt4iIiIhIham9RERERESkwlR0i4iIiIhUmIpuEREREZEKU9EtIiIiIlJhKrpFRERERCpMRbeIiIiISIWp6BYZBzM7zszuN7NjR7nfK8zsp2b22TI+91fLtS0RkTiY2cFmdlM5c2Oe57ndzN5VyeeoFDM73sx+N9rnjKRHU9wBiJTCzOYBXwTWAUcABwPvcfcbKvzUzwCHjOF+Kwgn1G8sx5NGPz27oxzbEhEBMLMjgC8Drwf+G1gKzAKagavc/c4KPO2TwG6UKTcWcCXhx07KLsrHpxJ+SOUJ4HfRUCvwd8CR7r6shKcY6+eMpISKbkm7nwHvdPeHzawB+H41ntTdnzCztWO43yYze76MT30M8D9l3J6I1Dl3f9DMfkL4Bb63Z9ab2duA/zSzT7n7lWV+zq1lzo0AmNlJwJ8yxa67/6Lcz5Hh7juAG83so8Bv3P2irDhuHc82zWwP4A3ufqO7P25ma8oUriSA2ksktcxsOvAyYAuAuw8CnyP8rGs1DJb5fmPxBuCeMm5PRASgP3eFu/8c+DzwL2Y2qwLPWdZcbWZ7E2a2q22X1w64G1hZzEbMrBX4AdCWtbqcnx8SMxXdklruvg54HLjDzI6K1i0H/mBmu5vZJWb2BzNbYGarzOyvZnZ45vFmdrqZfcbMfmNmX8pa32lmnzezG83sP8xst2j9BDP7rpl9zsyuI3w1WhQz+39RTO80syPM7GfR9i6JYrzRzNqj9RvN7J9zNtHi7tujbX3SzP7RzH5rZl8p+gUUERnddYQ2k4UAZra/mV1sZtea2a+iQhcza41y2WfM7H/N7O8zGzCzD5nZZ83s92b2kXxPNFJONrPTopx5tpn9zcy+Ymazzez7UT79XzN7Y7SJM4C9gQvM7OTo2Ju7zey9Wc9xjpl90cx+YGb/ZmaTzGw3M/tnM1tkZm80s8fN7Gkz6xjPCxa9JsdGM+GY2alm9k0z+7qZ3WFmk6N9+KaZ3RZ93jwfvcYvBxaa2T9mbXJa9Fm02czOyfOcDWZ2mZl1Rc/hZvafmb+PJIS7a9GS2gWYDywnzAZcD+wRrW8kJLB1hN663YG7gEcJ/2y+Brgouu80wmz524CpwE+ytv9b4KvR9X8BuqLrewIDhMQ6Woz/Spgt2oPQ+zcxWt8K3AncCswmfFjsAP4ZmAi8KrqduX+mXz2z31dG16cA/xT330KLFi3pXYCzgWfzjK0GvhXl1f8AGqP1PwZuiq5fAZwaXX8/0BNdPx14d3T9ZVGufmV0+1+Bz0fX8+Xk6cAG4BvAkcDxwNeAz0X3PR+4LStWBzqi623AX4Gzo9tvBn6cdd/bgWsIrbanRJ8Xp0b7eQfwxTG+dr+JPiu+AnyV0N99fNb4GmBudP1RwmdSK/AZ4Fng1cD/jfL+bzLxRvdfBnw9GjsdWJEnhg8At0fX94xe5zfF/b7SMnxRT7ekmrs/amavBC4FPgScZGbHe+jx7gE2ufstAGb2eeA+4CWED5hBM/tEtKk7Ccn9bcCMrPV/jR47GTgveizuvtrMniki1JnAt4H/6+5bo21sN7NVwFJ3fy56nlXAvR76HR8ifBjMJByQ+TbCBwTAVuAsM3sCuIrwD4eISCU4oRA9EtgX+KiZAawFdrNw41zChAHAvxHaKyDk2kU5uXYv4E85z3E2I+Rkd19nZhuBW939DwBmtgLYYmZ7EWaGJ48YtHuvmb2Qteoj0XYzrgV+Gq3fQPi8+I/oORZFcY7VPR71dJvZjYSJnoyF7r7CzF5LKLYnR/n/OeApd38AeCB67Ejbvi36THiwQEyHA73Rfq82s3UMb1ORBFDRLallZo3APh4OmPlY1PLxC6CbMHOQ69HocndgLvBDd/9xtO4b0TYvIBTB38h5rlcRkuXmcYb7MuB1wP7Aw1nrc3saB4YG3AejBJxpA5vhoaUGd19qZmcD3wH+ATgTeAERkTIys90J39I9Ssiba0fIj3sSZmJbANx9gDBRQPSYL7j7/dHtYY/NMmJOjjjDc+WzwAXA34DfA/sV2IXsxx0A/Drr9lOE1pk92DUX9zP+M6v8kfDPScYGM7sMuIkQe6ayzt2v0fSTv277JXCVmU0g/C16CbPmkiDq6ZY0G+ozBHD3h4ELgYPy3L+FkOCeBJ4HTsseNLNXR+vfZmZtOes3RTfnjTPW+4FvAj+KkmJRzGw2WQflmNkcwplbXkY4sPLHeR4qIlKK9xK+WbuJkB9fG80wA0P5cQ3QB5yYtf6VZjaFnFxrZk3RJEaufDl5JFcDyzycGnZ7EfuygjAzPvQUwEZC+0zZuPtgNDEyP/qn5RfAxe6+qJzPk/OctwG3Eb6RPR04xt03Vur5ZHxUdEvafTQqSDPmEHrxMqabWUt0/VhCP98G4N+BU8zsO2Z2bDQLsYWQHCcBt5vZm6ODVg4iFOoPAZeYWYuZzST0Uu9hZo1m9k9m9u48MVq0fI7wwXTFCGO598+9/XZC/2HGAcDJ7r4J+OgIjxERKcYuM6hm9nrgC8C57t5DmFVeDfzCzE6KDpb8P9HM9k8IZzk52cyOA97r7psJufZjZvYFMzuGMPmQac3Lzn/5cnJG9qzzYUC7mU0EXgu0mdlLorEdhAMPM5Mv2c9xNfB3ZjYtuv1q4Hvu3k+oh0bMxWZ2VBRXa57XroGceir6x+RthG83pxHaFg8E9gEmZh2kmTub3jdC/MODGqEHJfpbTAV+RZjNL+afEakStZdI2u0DPGpmPyd89WbAB7PG+4HPmNkGwqzwRwHc/S4z+xjhK8qTCQci/hnAzN4BfJfwIXI9cL67u5m9k3A6p0eBGwlfEx5K+DGJdxLaVob9KI+ZHUk4QGgu4evGnwFfiOK5EzgK6DCz/QkzMLOAd5rZY+yc9TkbmOTu3Tn7/n0LP2rRTzhwSUSkaFEeOYNQGH4P6CEcjLcb8BZ3/x2Au/dF+fF7hFx3K6GXG+DjhAmLHxC+2cucMeQ6Qu47j9AGd567vxAVlUcB+5vZDflyspmdRuhj/oCZPebuKwkHdX6NUHD/AHgH8ErgaeCHhFn5syyc1Wo+8FYz+4W7/4eZHQDcYma/JXxb+ploRn4hMCvav0cJPxS0l5nNj+J8F6FP/WdZr1szIfcfAsyJinIn9FKfGm1zCXAv8CChZeaXhJx+O+HgzUPM7GR3z2z3RsIPFb0YFemzgNPMbClhBpvo8bm/SbGCcKDoycAEoMHM7iYcTFmt0+jKKEx/C6lVFn4691/dvaNKz3e+u+sn2kVEakzU6tLq7on8cTIzO59wJpll0e0JhG8pLsscCyTxU3uJ1LKRWjcq80Tha9g/VuO5RESkeqIDRQ9PcME9AbgI6MhaPRPYooI7WdReIjUp6qd7N+HrwoXuflOFn/Kp6GtPERGpLU3uflXcQeTj7tss/BT9t6KD7J8jtMJcVPiRUm1qLxERERERqTC1l4iIiIiIVJiKbhERERGRClPRLSIiIiJSYXVxIOXMmTO9o6Oj6Mf19PTQ3t5e/oBSoJ73HbT/2v/k7P/ixYvXuHsygqmSesrZirl60hi3Yq6ecsVdKGfXRdHd0dHBokXF//pqZ2fnuB5XC+p530H7r/1Pzv6b2fK4Y6i2esrZirl60hi3Yq6ecsVdKGervUREREREpMJUdIuIiIiIVJiK7gK6urriDiE29bzvoP3X/tf3/qdVGv9uirl60hi3Yq6easRd9R/HMbMTgEuBv3f3ZdG6NuByYDMwDbjA3TeVMpats7PT09hfJCJiZovdvTPG51fOFhEZo0I5u6oz3Wa2JzAZODJn6ErgPne/ELgVuKYMYyIiUgLlbBGR8oljprsBGAD2c/dlZjYbWApMc/etZtYEbATmA33jGcvMxmRo1kRE0ioBM93K2SIiY1QoZ1f9lIHuPmhm2auOBda4+9ZovN/MlgLHADvGObaspCAHBuDVr4bBQfjjH0valIhImqUhZ29+aDOPf+BxJh88mYO+f1ApmxIRqZgknKd7DrA2Z91mYDbQP86xYXp6eujs3PlPR1dXV+GG+YYGWLw4XHeH4R84IiIV1d3dTXd3d+bmzDhjGUHicvZg7yAv/vFFGlp1bgARqb6x5uwkFN0O9OasayXMilDC2JD29vbiTnhuBo2NYca7vx+am8f+WBGREmUXmWa2JuZwciUuZ1tTmBjx/uq2S4qIwNhzdhKK7pWEo9izTQGeIyT38YyVrrlZRbeIyK4Sl7NVdItIGiThu7h7gNnRqaQws2ZgLnBvCWOla4r+H+nvL8vmRERqROJytopuEUmDqhfdtvOIHANw91XA7cBx0frjgVvcfeV4x8oSqIpuEZFU5OyhonuHim4RSa6qtpeY2RTg3dHNs83sO+7eA5wHXG5mhwCzgOwjZsY7VppMS8mOXdoNRUTqQlpytma6RSQNqlp0u/tm4KpoyV6/Hjgnz2PGNVYyzXSLSJ1LS862ZhXdIpJ8SejpTiYV3SIiqaCZbhFJAxXd+ai9REQkFVR0i0gaqOjORzPdIiKpoKJbRNJARXc+KrpFRFJBRbeIpIGK7nxUdIuIpEKm6B7cMRhzJCIi+anozkc93SIiqaCZbhFJAxXd+WimW0QkFVR0i0gaqOjOR0W3iEgqZIpuBsBdhbeIJJOK7nzUXiIikgpmBo3hug+o6BaRZFLRnY9mukVEUkMtJiKSdCq681HRLSKSGiq6RSTpVHTnk2kvUdEtIpJ4KrpFJOlUdOeTmelWT7eISOINFd07VHSLSDKp6M5HB1KKiKSGZrpFJOlUdOczYUK43L493jhERGRUDc3h40xFt4gklYrufDJF97Zt8cYhIiKj0ky3iCSdiu58MkV3b2+8cYiIyKjU0y0iSaeiOx/NdIuIpEbDhPBxNrh9MOZIRERGpqI7HxXdIiKpYa1hpntwm4puEUkmFd35qOgWEUkNzXSLSNKp6M5HRbeISGo0tEZFt2a6RSShVHTn09YWLlV0i4gk3tBMt4puEUmoprgDyDCzecBHgCeBA4Bud19iZm3A5cBmYBpwgbtvih6Td6xkmukWEckraTk7M9Pt23X2EhFJpsQU3cAPgXe4+0ozmwvcCcwDrgTucPcbzexE4Brg9OgxhcZKo6JbRKSQROVszXSLSNIlqb1kHjAlut4L7GZms4EzgNui9XcDC8yso9BYWaLRebpFRApJVM4e6unWgZQiklBJKrpvBq41synAu4GPAscCa9x9K4C79wNLgWNGGSudZrpFRApJVM7WTLeIJF2Siu4PAduAB4EX3f0WYA6wNud+m4HZo4wN09PTQ2dn59DS3d09ejQqukUkJt3d3UP5CpgZdzx5JCpnD52nWzPdIlJlY83ZSerpbgP+nZCYv2FmSwEnfG2ZrRXYEV0vNDakvb2dRYsWFReNim4RiUlXVxddXV0AmNmamMPJJ1E5WzPdIhKXsebsJBXdNwDvcvcNZuaEZP4JwhHu2aYAzxGSe76x0mVOGbh1a1k2JyJSYxKVs9XTLSJJl4j2EjObCRzi7huiVZcCm4DlwOzoNFOYWTMwF7gXuKfAWOmmRMcHbd5cls2JiNSKJOZszXSLSNIlougG1gHbzGxO1rq1wCPA7cBx0brjgVvcfaW7r8o3VpaIpk4Nl5vKc9pvEZEakricrfN0i0jSJaK9xN0Hzexk4HNmthjYEzjf3TeZ2XnA5WZ2CDAL6Mp6aKGx0mRmujdtAncwK9umRUTSLIk5WzPdIpJ0iSi6Adz9YeDcEdavB87J85i8YyVrbQ3L9u1hyRxYKSIiicvZDW2h6B7YOlCJzYuIlCwp7SXJlD3bLSIiidU4uRGAgS0qukUkmVR0F6K+bhGRVBgqul9U0S0iyaSiuxAV3SIiqaCiW0SSTkV3ITptoIhIKjROiYruzSq6RSSZVHQXopluEZFU0Ey3iCSdiu5CMkX3xo3xxiEiIgWp6BaRpFPRXciMGeFy7dp44xARkYJUdItI0qnoLqS9PVz29MQbh4iIFNQ4sREMBrcO4gP6VUoRSR4V3YXMnBkuVXSLiCSaNRiNk3SubhFJLhXdhWRmutesiTcOEREZlVpMRCTJVHQXovYSEZHU0GkDRSTJVHQXoqJbRCQ1mnZvAqB/Q3/MkYiI7EpFdyGZnm61l4iIJF7T9FB071i3I+ZIRER2paK7kOxTBg7o60oRkSRrnt4MwI61KrpFJHlUdBfS1BRmu93VYiIiknDNM0LR3b9O7SUikjwqukezzz7h8pln4o1DREQKGmov0Uy3iCRQWYpuM3trObaTSHPnhssVK+KNQ0SkTGo1Z2faSzTTLSJJ1FTMnc3sROAiYBY7C/YGYE+grbyhJYSKbhFJqXrL2U0zdCCliCRXUUU38EPgy8BDwGC0zoAzyhlUomTaS1R0i0j61FXO1oGUIpJkxRbdT7v7FbkrzeyxMsWTPJrpFpH0qqucnenp7l+r9hIRSZ5ie7rPM7PTRli/oBzBJNK++4bLp5+ONw4RkeLVVc5umdUCQN/qvpgjERHZVbEz3a8FLjGzG7PWGeBAdzkCMrOjgNcBTwL3AVuAy4HNwDTgAnffFN23Ld9Y2bzsZeHyiSdgcBAadMIXEUmNusrZQ0X38334oGMNVq5Ni4iUrNii+2KgC1gEZH4txoD3lCMYM3sfsL+7X5S17vvAHe5+Y3RQ0DXA6dHwlQXGymP69PBz8D09sHLlzh5vEZHkq6uc3TihkabpTfSv62fHmh207NFSrk2LiJSs2KL7l8Ct7r49e6WZfa/UQMzsGOBM4ISsdbMJB/x8OFp1N3CLmXUAffnG3H1ZqfEMc9BBoeh+7DEV3SKSJnWXs1tnt9K/rp/tz21X0S0iiVJsr8RPgE+b2euzlmOAT5chliuAvwLfNrM7zOy1wLHAGnffCuDu/cBS4JhRxsrrwAPD5eOPl33TIiIVVHc5u2V21GLynPq6RSRZip3p/ixwAHA2oScQwleVs4CPjzcIMzsQOAx4n7s/YmbnA3cAlwJrc+6+GZgN9BcYK69588Lle5+OyQAAIABJREFUn/9c9k2LiFRQ3eXs1tmtAGx/bvso9xQRqa5iZ7q/Cuzj7h3uvl+0dACnlhjHfGCduz8S3f4OOw/26c25byuwY5SxYXp6eujs7BxauruLPH7o8MPD5eLFxT1ORGQcuru7h/IVMLOETdVdzh6a6V6pmW4RqY6x5uxiZ7qfcffNI6yfVeR2RopjKBZ37zWzvwGNhCPcs00BniMk8Hxjw7S3t7No0aLxR5cpuh95BPr6oEV9giJSOV1dXXR1dQFgZmtK2FTd5ezWfcJM97YV28b1eBGRYo01Zxc70/2h3BVmthD4RpHbyfUIsJuZ7Zm1rp+QjGdHp5nCzJqBucC9wD0Fxspr6tTQ193XpxYTEUmTusvZbfuHX7fvfTJ3Ul1EJF7FFt3TzeyzAGY20cyuBa4FHiwlCHd/jHCU/WnRtncHWoAfAbcDx0V3PR64xd1XuvuqfGOlxJJX+MoASpkxFxGprrrL2UNF91MqukUkWYpqL3H3t5nZUWb2beDNwCrgEGBjGWI5C/immbUQZj9Od/d+MzsPuNzMDiF8JdqV9ZhCY+V1xBHwox/B/fdDV+WeRkSkXOoxZ7fObcWajL6VfQz0DtDY1ljupxARGZeCRbeZvT7P0HLgb8DXgH2Bk4F/KCUQd19DOOdr7vr1wDl5HpN3rOze+MZweddd4A6mXzoTkWRRzoaGpgZa921l21Pb2LZ0G5NePqmSTyciMmajzXRfDswDNgCDI4xf///bu/M4ucoq/+Of0/uWpOlsnT0EshEwIAEUQQKyCQguqIw4jDoQGDecEX+g8nP3pSMy4zpKUHSUETf4DcimoMgiKATZlxCykaSzr5303nV+fzy3uqs71Ut1V3VVdX3fr9fzqrpL3XtuVdfJya3nPpdwxfoUhpnAc96RR0J9PTQ0wIsvwqJF2Y5IRKQ35WxCF5OW1S00v9qsoltEcsZARfeXgBXuvq2/lczswvSFlKPM4Iwz4Oc/h7vuUtEtIrlIORuoWlDF7j/s5sCLB5hw/nBGXBQRSZ9+L6R097sHSt7Rer9NX0g57B3vCI+/LYzDFZH8opwd1LyuBoADzxzIciQiIt1SHb2ksJ19NlRXwxNPwPr12Y5GRESSqF4cupTsf2Z/liMREemmojsVlZVw3nnh+W9+k91YREQkqepF1VAETSub6GzpzHY4IiKAiu7Uvfe94fGnPw2jmIiISE4priymal4VxODA8+piIiK5YdhFt5kdY2Y3m9knzWz28EPKceedB5MnwwsvwGOPZTsaEZGUFErOHrNkDAD7/rovy5GIiATpONN9BXC6u18PnBi/xe+oVVoKH/xgeH7DDdmNRUQkdQWRs8edNA6AvY+k4z5AIiLDl46i+6PAsdHzXwKladhmbrvsMigqgltugXXrsh2NiEgqCiJnJxbdrq6AIpIDUi66zeytZnZx9Px4YK67bwJw95i7j/7f8ubMgYsvhvZ2+OIXsx2NiEifCjVnVy2souSQEto2tdGyviXb4YiIpFZ0m9lXgd8A7wBw98eBS8zsfRmILbd94QtQUgI/+1no3y0ikmMKOWdbkTHu5HC2e/f9u7McjYhI6me6TwSmAw8nzLsH+HraIsoXc+bAsmUQi8FHPqKRTEQkFxV0zh5/zngAdt21K8uRiIikXnTf5+57gMQKc1TfTrhfX/4yTJwIDz4IN9+c7WhERHor6Jxdd04dALvu20WsNZblaESk0KVadK8xs88A883sIjO7FfgI8O/pDy0P1NXBddeF5//6r7BlS3bjERHpqaBzdsWMCqoXVxM7EGP3A+piIiLZlVLR7e6/BP4EFAOXAHuAM939+xmILT9ccgmccQbs3BmGElQ3ExHJEcrZMPGdEwHY9ottWY5ERApdyqOXuPtf3f0Kdz/H3f8ZeMDM3pKB2PKDWbg7ZV0d3HsvfL9g/i0TkTxQ6Dl78sWTAdh+23Y6D+iW8CKSPSX9LTSzO4CaAbYxAdgB/DFdQeWdqVNh+XK48EL41Kfg1FNh0aJsRyUiBUY5+2CVh1Uy9sSx7Ht0H9tv2079P9ZnOyQRKVD9Ft3As9E6K4EYcBnwFLAiYZ1DgbqMRJdP3vWu0L3kJz8Jz//6V6itzXZUIlJYlLOTqL+knn2P7qPhBw0qukUkawYquq8DLLr6HTM7zN0/13slM7s9E8Hlne98B1asgOeeg4sugjvvDGN5i4iMDOXsJCa/fzJrrlnDvsf2se/xfYw9fmy2QxKRAtRvn2533xtP3pHJZlacuI6ZLQWOz0Bs+aemBu64AyZMgN//Hq64IozjLSIyApSzkyuuLmbKZVMA2PitjVmORkQKVaoXUv4WeMbMfmhm15nZncB9wHfSH1qemj0bbrsNKirgxz9W4S0i2aScHZn20WlQDNt+vY2mV5qyHY6IFKBUhwy8DziTcBHOfMLwUxe6+9cyEFv+Ovlk+N3vQuF9443w4Q+r8BaREaec3a1iZgX1H6iHTlj3+XXZDkdEClDKHY7dvQG4Nj5tZkVmdqy7PzncYMysDHgCuNLd/2xmlcA3gEbgEOBqd98Xrdvnspxw+umhq8n558MNN4ShBb//fShKeZRGEZEhy2TOjraXN3l79udms/XnW9n2y23MvGYmNYsHGuhFRCR9Uiq6zeymJLMnAp3A29MQz6eA2QnT/wXc6+6/MrOzgRuAfxjEstxwxhlw++2h8P7hD6G1NTyWlWU7MhEpACOQsyGP8nbFzAqm/stUNn17E6s+voqj/3w0ZjZSuxeRApfqade5wGvA+oRWDqTjLPeJwGZgdzQ9FXgf8LtolfuB881sdn/LhhtH2p15ZjjjXVERhhM8+2zYtSvbUYlIYchYzob8zNuzPz+b0oml7H1oL1t/vnUkdy0iBS7V7iWXuvvKxBlmtgD4wHCCMLNqQj/DfzOz+PBWS4Ed7t4E4O4dZrYWOAVo72fZuuHEkhFnngkPPRTOeD/wALzhDeFiyyOPzHZkIjK6ZSRnR9vJy7xdekgph113GC9/4GVWX7WaunPqKJugXx9FJPNSvZByZZLZrcClw4zjauDrveZNA3b2mtcITB1g2UG2b9/OkiVLutry5cuHGe4QHHccPP44LF4Mq1aF6RtvBPeRj0VEctry5cu78hXhDpJDksGcDRnM25nO2ZMvmUzt0lrat7ez8p9X4srDIjIMg83ZqfbpXgskZqdiYDLw86EEGW3zLGCFu2/rtciB5l7zyglnSxhgWQ8TJ05kxYoVyRaNrBkz4C9/gY9/HG66CZYtg/vvh+99DyZOzHZ0IpIjli1bxrJlywAwsx1D3U4mcna03Yzm7UznbDNjwU8X8MTiJ9h5x042L9/M1MuTnrMRERnQYHN2qt1L7gN+QXcSjwFb3f2VoQQZuQo4JuFilkOA24H/iJ4nGgM0RPvva1luq64O43efdloYw/vXv4Y//Qm+9S143/vCKCciIumRiZwNoyBvV8yqYN4P5/HSP7zEq594lZrX1zD2ON2pUkQyJ9ULKW929z+7+4NRexgwMzt8GDG8Dzg6oTUQfvr8CTA1GmIKMysFZgIPAg/0syw/XHwxPP10KL537ID3vx/OPRfWr892ZCIyemQiZ8MoyduTL5rMlMunEGuJ8fwFz9Pa0JqtUESkAKRadJ+dZN5W4DdDDcDdt7v7xngjDGW13d1fA+4E3hKtejpwq7tvcvctfS0bahxZcdhhoXvJj38MtbVwzz2waBH8+79DS0u2oxOR/Jf2nA2jK2/P/c5cxr15HG2b23ju/OfoaOzIZjgiMooNqug2s4+a2R7gajPrTGyEC2MyNe7SvwAXmNlngXOAZYNclj/M4EMfgpdeggsvhAMH4JprYMECuOUWXWgpIinLYs6GPMvbRWVFLLp1ERVzKtj/5H6ee9tzdDZ1ZjssERmFbLBXbZvZ64B/Ar7Ta1Gbu29Od2DptGTJEs+JCykH4w9/gE9+Ep5/Pkwfdxx88YthfG/19xYpOGb2pLsvGcLrlLNT0LymmadOfoq2hjbq3lrHotsWUVxRPKIxiEj+6y9nD7p7ibs/C1zj7ut7tc1mNjNt0Ra6M88Mfb1vvBHq6+GJJ+Ccc+D448NNdnTmW0QGQTk7NZVzKll8/2JKJ5Sy655dPHfOc3TsU1cTEUmffs90m9lVwN3u/mI0/VnCkFM9VgNOdfelmQpyuPLqTHei/fvDbeOvuw62RSNzHXkk/Nu/hZFOysuzG5+IZFwqZ7qVs4dv/3P7efasZ2nb3EbN62t43T2vo2ySbp4jIoMznDPdhwK1CdOvB04E5kTL4q0+DXFKbzU1cNVVsHYt/Od/wpQpodvJhz4Es2fDV78KO3vfa0JECphy9jDVHFXDMY8cQ8VhFez/+36ePP5JGp9uzHZYIjIKDLpPN4CZvR543t3bes0/In5mJRfl7Znu3lpb4Ve/guuvh2efDfPKy+Hd7w432jnpJPX7FhllhtqnO3qtcvYQtW0No5k0Pt5IUWUR82+az+SLJmc1JhHJfWnp0x05Kknyng8cNtTgJAXl5XDJJaHP9/33h77ebW1w883w5jfDEUeEM+I6+y0igXL2EJVNLuPoB4+m/gP1xJpjvPQPL/HKv7yikU1EZMgGO2TgyWb2ZuAMM3tzYgMWAt/PaJTSkxm85S1w112wejV89rPhosuXXw79vadMgQsugF/+MgxBKCIFRTk7PYoripl/03wO/+7hWJnR8MMGVrx+BY1/V3cTEUndoLqXmNlRwM8IfQF391rcBvyPu38p/eGlRy78VJlx7e1w552wfHkYdjAWC/Orq+Htb4f3vCeMjFJRkd04RSQlQ+leopydfo1PN/LS+16i6aUmrNSYefVMZn5mJsWVGlZQRLr1l7NTGad7HPBud/9ROoMbCbmYwDNqyxb49a/DzXX++tfu+TU1cN558K53wVvfGgpyEclpwxinWzk7zTqbOln9f1bT8P0GACoPr2TuD+ZSd3pdliMTkVyRlqJ7gB2c5O6PDHtDGZKrCXxErFkTupn89rfw1FPd8ysrww133va2UIDXazADkVw0nAsp+9mmcvYw7P3LXlYuW0nTi00ATHzPROZ8bQ6VcyqzHJmIZNuQi24zuwOoGWj7wCJ3nzT0EDMr1xP4iFmzBm67DW69tecZcIBjjw0XZp5zTrgLZrF+MhXJBSmO062cPUJibTE2XL+B9V9aT6wlhpUa0z42jVmfnUVpXWm2wxORLBlO0X0dsAfYBPS1ogHnufuFww00U/IhgY+4jRvh9tvh7rvhT3+ClpbuZRMmwFlnwemnhws2Z8zIXpwiBS7Fols5e4S1bGhh7bVr2frzreBQUlvC9H+dzrSPT6O0VsW3SKEZTtE9AWh099YBdjDe3XN2nLp8SuBZ0dwMDzwQCvC77oJ163ounzsXTjstFOCnnhqKchEZESkW3crZWdL4VCOrr1rNnj/tAaB4XDHTr5zO9Cun68y3SAEZiT7dl+byxTr5mMCzxh1WroR77w1nwP/8Z2jsNTzW4sWwdCm86U1w4okwbVo2IhUpCBnq062cnSF7HtzDui+t6y6+a4qp/1A90z42jarDq7IcnYhkWsaKbjObCXwE+Ki75+xQGPmcwLOuowOefBL++MdQhD/ySLgzZqJZs0IBHi/CjzpKfcJF0iSdRbdy9sjZ88ge1n95Pbv/EI3YaDD+vPFM/8R0ak+txXT3YJFRKa1Ft5kVA+cDlwFnAFuj7eTs6c7RkMBzRksLPPYYPPwwPPpoeL5vX891qqvDhZnHHx8uyjz++FCY6x8ZkZQNt+hWzs6u/c/uZ+O3N7L1f7bireHf26oFVdT/cz31l9RTNqksyxGKSDqla5zuQ4FLgQ8C44H1wMeB3wNvdPdH0xNu+o2mBJ5zOjvhhRfgL3/pbr37hEPoB55YhC9ZApNydvAEkZwxjHG6lbNzSNu2NhpuaKDhBw20bW4DwEqM8eePZ8o/T6HurDqsWCcmRPLdsIpuM3s34QzJacCrwI+A/wY+4+7/muZYM2I0JvCctm0bPPFEaI8/HtrOJNdsTZ0KRx8NxxwTHo8+GubMgaKikY9ZJEelWnQrZ+e2WHuMXffsYvOPNrPz7p3QGeaXTi5l0nsmMemiSYx9w1isSAW4SD7qL2eXDOL1zYTbBj8KXOrur0QbHf4VmDI6TZoE554bGoSLM9etC8V3vBD/+9+hoSG0u+/ufm1NTbhQM16IL14MRxwBVboASWSQlLNzWFFpERPOn8CE8yfQ2tDKlv/ewpafbKF5VTObvruJTd/dRPnMciZdNIlJ751EzTE16v8tMkqk0r1kJuHsyWzgduAUd/9YtMw8HcOgZMhoPmuSt2IxWL0ann66Z2toOHhdM5g9OxTfiW3hQhgzZsRDFxlJw+heopydJ9yd/X/fz9ZbtrL9V9tp3dh9sXr5rHImnD+B8eePp/bNtRSV6ZdAkVyW7gspS4C3AxcDzwM/Bs509+XDDTRTCi2B57Vt2+CZZ3oW4qtWQXt78vVnzOguwhcsgHnzQpsyRRduyqiQhgsplbPziMecvY/uZdst29h+63bat3bnvuJxxYx/63jGnz+eurPqNP63SA7K5JCBhwMfAy5z95z9/b+QE/io0N4Or74KL77Ys738MrS1JX9NTU24qc+8eTB/fncxPncu1NaObPwiw5DmIQOVs/OIx5zGJxrZcfsOdtyxg6YXmroXGoxZMoZDzjiEujPrGPvGsToLLpIDRuLmOJe4+8+GvaEMUQIfpTo6YO3aUIC/8AK88kp3S3bhZtykSd1FeLwQnzMntLFjRy5+kUHI0M1xlLPzUPPqZnbcsYOdv9vJ3r/sxdu6//0uqi6i9pRa6s6so/bUWqqPrNbFmCJZkPGiOx3M7J3A9cBY4Gbgk+7eEd3W+KvAHqAUuNrd26PX9LkskRJ4Adq5M3RLSSzE4625ue/XjR/fXYDH22GHhcfp03XTHxlxmSi600E5O7s6D3Sy56E97L5vN7v+sKvnWXCgpLaEcSeNY9ybx1H75lpqXl9DUanOhItkWs4X3dEFP18Bvg0sAn4AfN7dv2lmDwDXuPvfzOwKYJa7fzp6XZ/LEimBS5dYLFysmViEr1oFa9aE1tLS92tLSsIFnfFi/NBDQ5/y6dNDmzoVystH7FCkMORi0a2cnXtaG1rZff9udt+3mz0P7aH1tZ53Di6qKmLcieMYd/I4xr5hLGOOH0NprfqEi6RbPhTdJwOPuXtHNP0NQiL/MnCbu0+N5k8DVgKTgaP6WubuBxK3rwQug+IOW7Z0F+C9W7KRVXqbNKm7CE9s8eJ82jSorMz8sciokaNFt3J2jmtZ38Keh/ew96G97HloD80rD/6Fr3J+JWNPGMvYE0IRXvO6GvULFxmm4Y7TnXHu/nCvWZuAMYSbO6xNWG9TNF7pEuBN/Sx7MMMhy2hkFkY9mTIF3vSmg5c3N4fxxuNF+Nq1sHFjd2toCKOvbNsWxiHvy/jxyQvzxFZTk7HDFBku5ezcVzGrgvpZ9dS/vx6Atq1toQh/ZC+Nf2uk8alGmlc207yyma0/2wqAlRtjXj+GMceOoeaYGmqOqaH6iGqKylWIi6RDThTdSRwHXAcsA3pfEdcITAWm9bOsh+3bt7NkSfd/OpYtW8ayZcvSGa8UgsrKMDb4woXJl3d2wtatPQvxxLZhA2zaFPqb79wZhkbsy5gxobvKlCnhMd4Sp6dMgerqzByrZNXy5ctZvrxrRL8J2YxlkJSzc1zZ5DImXTiJSRdOAiDWFmP/M/vZ97d9ND7eyL6/7aP5lWb2PbaPfY/t63qdlRpVR1Qx5pgx1BwdCvGaxTWUjMvV8kFk5A02Z+fctyYa0mq7uz8T3UGt929i5UA70N+yHiZOnIh+qpSMKy7uLoiPPz75OrEYbN/ed2Eeb42NsHJlaP0ZO7bvojw+PWWK7uiZZxKLTDPbkeVw+qWcnZ+KyooYe9xYxh7XPWJT+652Glc0sv+p/TQ+1cj+p/fT/EozB545wIFnevQAomJOBdVHVlO9KLSqRVVULaiiuEIXm0vhGWzOzqmi28yKgcuBa6JZm4B5CcsNqAYaBlgmkpuKimDy5NCOPTb5Ou6wd2/orhJvmzcnn963L7SXX+5/v7W1yc+U19eHNnlyeBw3TjcVkkFTzh5dSutKqTuzjroz67rmdezv4MCzB9j/9P6uYvzAcwdoWdNCy5oWdt6R8ONFEVQeVtlVhHcV5POr1EVFhBwruoGrgG+6e/yy6/uAKxKWzwbagBXRY1/LRPKXWSiSa2vDnTb74g67dx9clCcr0PfsCe3FF/vfd3l5dwEef+zrufqdi3L2qFdSUxJGPTlxXNe8WHuMppVNNL3QxIEXDnS15lXNXY3/TdhIEVQcWkHVvCqq5ldROa+y67F8Wjmm/+hLgciZotvMPkNIvpVmNgdYCjwEbDazI939eeBs4Lvu3gasMLO+lomMfmZQVxfakUf2vZ477NqVvCjfujW0LVtCa2yE114LbSDV1aEAnzChZxs/Pvl0XR2Uaoiy0UI5u3AVlRZRc2QNNUf2/I93Z0snzSubexTiTS800by6mZbVLbSsbmHXPbt6bquqiKp5USE+r4rK+ZVUza2i4rAKSseXqiCXUSVXhgy8ljDUVKKX3X2hmc0CPgesA2qAzyYMU9XnskQafkpkkJqauovwxGI82XR/Y5r3pba2Z1HeV4Eeb3V1YXz0ApajQwYqZ8ugdbZ00rK6haZXmmh+pTmcJY+et28/qEt/l+IxxVTMqaDy0MrwOCc8VhxaQcXsCvUfl5yU8+N0Z5oSuEiauYez4lu2hJFYduwIra/nO3aEs+2xWOr7qq0dXIGeeEZ9FBXquVh0Z5pyduFo39VO86pmml5pomllVJS/0kTL2hY693X2+9qyaWVdhXjloZVUzK6gfGY55TPKKZ9erqJcsiLnx+kWkTxjFkZOGTt24HXjYrHQr3wwBXp8eufO7v7or746+H3FC/VkBfr48TBxYncf9fp6qKhI/T0QkWErrSul9IRSxp7QM5e4Ox27O2he00zLmpaux5a10fP1LbRtaqNtUxt7H96bfNuTS6mYGQrxihnR48zux9KJpViRuq/IyFHRLSIjo6iouw/6vHkDrw9h7PM9e/ouyvs6o55qoV5b27MIT2zxEV9mzQrjp4tIxplZKMjrShm75OD/3Mc6YrRuaO0uwte00PJaC62vtYbHja20b22nfWs7jU80Jt9HmVE+o7sQL59WTvnUcsqmlXU9L51cSlGJRl6R9FDRLSK5q7g4nJkePz71Qr2/An3bttBHffPm0EUmXqQPNPRiXV0ovmfPDo/xdthhIT6dMRcZEUUlRVQeWknloZUcctohBy33Tqd1cyutG1q7C/HExw0tdOzs6LrAs+8dQVl9KMLLpkbFeK/n5dPKKR5brIs+ZUAqukVkdEks1OfPH3j9WCwMvZh4keiWLd0F+ebN4YZF69eHs+i7dsFTTx28naIimDOn+66lxx4LJ58czpSLyIiyYqNiegUV0yvgjcnX6TzQScuG7mK8raGN1k2ttG5q7Xrevq2dtoY22hr6H2SnqKqoZ0E+tZyyKWWU1Uctel5SW6LivICp6BaRwlZU1F2kL1rU93qxWDhDvn59d1u3LjyuWhW6ssTb737X/bojjoArr4RLLw37EpGcUFxdTPWCaqoXVPe5TqwtRtvmNlobomJ8U1SYNyQ839RKrCnWPUZ5P6zMehTh8VY+pbzHdFl9mW4oNAqp6BYRGYyiou5+3ieccPDy1tZQfL/4IrzwAjz2WGgvvgiXXw4PPgg336w7forkkaKyIipmVVAxq++uY+5O577OnmfJG1pp29LW3TaHx87GTlpfC11eBlJySMlBZ8p7TE8qo2xyGSXjS9TvPE9oyEARkXSKxWDTpnDG+8UX4YtfhO3bw7KXXoIFC1LanIYMFBk9Og900ra1uwiPt9bNPYv09q3teMcg6zOD0vGllE4upWxSGaWTSimbXNbjeemk0q4ivbhaQylmkoYMFBFJt5YWePppePJJeOUVWL06tLVrw1nv3hYvDn2+RaRgFVcXUzmnkso5lf2u5zGnfWd70rPlXc+3tdG+rZ32ne207wit6YWmAWMoqioKBfkgivTS8aVYsX6dSxcV3SIinZ3hbpwHDsD+/aHt3t3ddu0Kj9u2wYYNoa1aBR0H3UwxmDQpjGgyf344s33ssXDKKVBaOrLHJSJ5yYqMsolllE0sg6P6XzfWEQtF97Z22raGQrxtW1uP5+1bu+fFmmK0rGuhZd0g7ipcBKUTSvs8a95VnEfLiqt0Fr0/KrpFJDe5Q3s7NDen3pqaDp534EDPtn9/9/Oh3NLeLFx4edxx4WLJww4Lbc4cjectIiOmqKSI8vpyyuvLB1zX3els7Ow6Sz5Qkd6xqyOcTd/WPrhYqou6i/BkZ9ETHycU3ll0Fd0i0i1e6La1dT/GW+L0YJ63tYVidrBFcbI2lNvGD4UZVFVBdXVoNTVwyCGh1dV1P58wAWbMCG3OnLCeiEieMDNKxpZQMrYEDh94/VhbOIueWIjHi/W2rW20b+85L3YgRsvacOfQgYMJZ9EHVaBPLqW4Ov/HQlfRLTJY7qEIbG8P3Qo6Ovp+3t+ybD2PT/dXKPfVXSJbSkqgsnLgVlU18DrxgjpeVCdOV1ZqVBERkV6KyooonxruzjmQrrPoiWfMez0mLuvY2UH79nbatw+yL3pl0UHF+EFdXCaF0V1KJ5ZiRbmX01V0y9DFYskLu8HOG+rrsrH9+HQhKC0NrawstFSfx6dLSwdXMPdXOJcoRYmI5IMeZ9HnDrx+j77ovc+k937c2kasOUbr+lZa1w883KKVWbhZ0YyoTS+nYkZFeD6znKq5VVkZxUX/omVLLBZGOGhtDWcY4897Tw92Wfws5UgVnx0d4cxvoSkqCoVgSUkoKkfqeTq31V/RXFqqM7449DfZAAAUtUlEQVQiIpJxKfdFP9DZZ0GeON3a0ErHro4Bu7mUzyqnakEV1QurqV5cTe3JtVTMqchoFxYV3YPlDo2NYbzdeNuxI4xokHhBVuLz3tNNTaGPa2trGC0h35kdXBgmKxZHw7zi4jCtOwqKiIiMKDOjpKaEkpqSAYdbBOhs6qR1YyutG1pp2dBC64bW7um1LTS/2tx11nz373d3va5sahmT3juJqR+eStXhVWk/DhXdfWlpgdtug9//Ptzg4qWXQuGcTuXl3a2sLLXp3vPiZylHsjBVASoiIiI5priqmKp5VVTNS144x9pjtKxpoenlJg68eIB9f9vH3kf20tbQxsb/3EjDDxpYePNCJr5rYlrjUtGdTEsLnHRSuOlFoqoqmDixZ6urS35hVrLpqiqoqAhFckmJfsYXERERGWFFpUVUza+ian4VEy6YAIQuLI0rGtn4HxvZ9sttrLx0JePfNp6isvSdYFTRncy994aCe/Jk+PSnYckSWLgwFNgiIiIiMqqYGWOPG8vCXyxk3+P7wpnwlU3UHJW+oWFVdCezcWN4fOc74corsxuLiIiIiIwIM6OsvoyWNS107EnvqGXqlJtMfGi4srLsxiEiIiIiIyrepcTb0ztKm4ruZOJ3wdOFgiIiIiKFJRrC22PpLbrzvnuJmVUC3wAagUOAq91937A2qqJbRCQjMpKzRUTSqOtulmke3Xk0VJX/BTzi7p8BbgduGPYW40W3RhcREUm39OdsEZE0suJQ/3mnupd0MbOpwPuA30Wz7gfON7PZw9pw/E6LOtMtIpI2GcvZIiLpFJV/6e5eku9V5VJgh7s3Abh7B7AWOCWbQYmISFJLUc4WkRwXP9Ot7iU9TQN29prXCEzNQiwiItI/5WwRyXld3Ut0IWUPDjT3mlcOtCfO2L59O0uWLOmaXrZsGcuWLct8dCIiQ7B8+XKWL18en5yQzVjSTDlbRHJfvHvJIPt0DzZn53vRvYlw9XuiMUBD4oyJEyeyYsWKEQtKRGQ4EotMM9uR5XDSSTlbRHJeqt1LBpuz8717yQPA1GgIKsysFJgJPJjVqEREJBnlbBHJffEz3a4LKbu4+xbgTuAt0azTgVvdfVP2ohIRkWSUs0UkH2RqnO58714C8C/AN8xsMVAPqOOfiEjuUs4WkZymCyn74O67gcuyHYeIiAxMOVtEcl68H0gsM5sVERERESl48e4lujmOiIiIiEimpDhkYIqbFRERERGRriED1b1ERERERCRD4me61b1ERERERCQzuoYM1JluEREREZEMUZ9uEREREZHMUp9uEREREZEM05CBIiIiIiKZppvjiIiIiIhkVteZbvXpHgE1NTBlCowZk+1IRERERGQEFY8rpmxqGcXVxWndbklatzZafOQjoYmIiIhIQZl97WxmXzs77dvVmW4RERERkQxT0S0iIiIikmEqukVEREREMkxFt4iIiIhIhqno7sfy5cuzHULWFPKxg45fx1/Yx5+v8vFzU8wjJx/jVswjZyTiVtHdj3z9w0mHQj520PHr+Av7+PNVPn5uinnk5GPcinnkqOgWERERERkFzD29d9vJRWa2HVg/hJdOAHakOZx8UcjHDjp+HX/uHP8sd5+Y7SBGUoHlbMU8cvIxbsU8ctIVd585uyCKbhERERGRbFL3EhERERGRDFPRLSIiIiKSYSq6RUREREQyrCTbAeQiM6sEvgE0AocAV7v7vuxGNXxmdgbwFeC97r4umtfnsQ51Wa4xs3cC1wNjgZuBT7p7h5lNAL4K7AFKCcfQHr1mSMtylZmdAiwHJgH/7e6fiOaP+s8/zszKgCeAK939z4V07KNdPnwmQ8lDuSSV70+uMLM3ACcBrwKPAAfI4ZjNbCHwUUK8c4Hl7v50Lr7X6awnshxz0u9ltCz93013V+vVgJ9EHwrA2cAt2Y4pDcc0GXgH4MDswRzrUJflUgNmAj8DjgUuISTdq6JlDwAnRM+vAL6W8LohLcvFBowDPkNIdOcC7cDphfD593ofPgvsBZYW2rGP9pbrn8lQ81AutVS+P7nQgA8CX8mzv5MVwLSEv5mXcjFu0lxPZCvm/r6X0fK0fzez/keWaw2YCrQCVdF0SfRBzM5mXGk6tqJef3B9HutQl2X7GJMc88lAScL0N4C7gDcADQnzpwH7geqhLsv2sfbzHlQSjVQUTT8JnFoIn3/CMZ8IfAhYBywtpGMf7S0fPpOh5KFsx9wr/kF/f7IdaxTPKcD9vfJeTsccxXQAWBA9nwg05GrcpKmeyHLMSb+X0fOMfDfVp/tgS4Ed7t4E4OFnhrWEL3Fec/dYr1lL6ftYh7osp7j7w1F8cZuA14DTCDHH19sUPV0yjGU5yd2bPcoaZlYNvAT8mQL4/KHrmC9095sSZi+lAI69QCwlxz+TIeahnDCE708u+A9Cnvuumd1rZm8i92MG+C3wIzMbA7wf+Bg5Gnca64kR0zvmfr6XkKHvporug00Ddvaa10j4n9po09+xDnVZrjsO+CEFeOxmdjJwD6FvWhWF8x5cDXy917xCOfZCkI+fyWDyUK5I9fuTVWY2HzgGuNHdPwr8EbgXmEGOxpzgw0ALoe/8fne/lRx+r3sZDXkz/r2EDMWsovtgDjT3mldO6Ac72vR3rENdlrPM7HBgu7s/Q4Ede+RV4KfAW4DrKID3wMzOAla4+7Zei0b9sReQvPpMUshDWTfE70+2LQJ2ufuz0fT3ACO3Y46rBG4BfgF8y8xOJz/ihjzPm72+l5ChmDV6ycE2ES44SzSG0LdqtOnvWH2Iy3KSmRUDlwPXRLM2AfMSlhuhz3bDMJblNHffDNxkZjHgU8DDjP7P/yrgmPAxASHu2wk/P4/2Yy8UeZOzU8xDuWAo359sKyGhtnH3ZjNbBRSTuzHH3Qxc5O57zMwJBfgnyP24Yej1RNYl+V5Cpr6bI9mJPR8aUE/oLF8ZTZcSLgCYlu3Y0nR8iRcR9HmsQ12W7ePr57ivBiYnTC8B1idMH0r46ahsqMuyfYwpvBdHAn8ohM+fcDHS9IS2AXg34ar1UX3shdLy6TNJJQ9lO9YonpS/PzkQ8wLCv3OJ7/MTwD/lasxRPBOAzQnTBqwmDHmYk3GThnoimzEnzOvxvYzmZeS7qe4lvbj7FuBOwk/wAKcDt3p3J/q8Zd2nKwz6P9ahLhuBw0iZmX0GeAqoNLM5ZvYhwribm83syGi1s4Hvunubu68YyrKRO6LUmFmVmR2TMOsc4NuF8Pm7+3Z33xhvQCfhJ8TXGOXHXijy5TNJNQ9lK85EQ/n+ZCvWOHd/mXDtyoUAZlZLOGHyP+RozJFdQIuZTUuYtxN4lhyMO131xAiGfFDM0byDvpdmdnim/q23qIKXBGZ2CGHomHWE/6F92t33ZzWoYUq4Gvq/gC8B33P37f0d61CX5RIzuxb4cq/ZL7v7QjObBXyOcAw1wGe9e1D8IS3LRWZ2PGF4speBvwKPu/tvomWj+vPvzczWAR/wcHOPgjr20SzXP5Oh5qFcM9jvT7ZFNzX5NmHc65mEiypfzOWYAcxsMeFiyicJ40o/5O4P5lrc6a4nshUzoUtJ0u9l9Jq0fzdVdIuIiIiIZJi6l4iIiIiIZJiKbhERERGRDFPRLSIiIiKSYSq6RUREREQyTEW3iIiIiEiGqegWEREREckwFd0iIiIiIhmmoltE8oaZTTSzO8xsvZldk+14RESkf8rb3XRzHBHJG2Z2FvAAMJZwh8157r4ru1GJiEhflLe7qegWkbxkZg8Bp+XqLbNFRKSnQs/b6l4iInnHzGYAfyjUxC0ikm+Ut1V0ywgzs/9jZt/Pdhxx6YjHzE43s8fMbGmawhpVzOybad6eAe8Bvp7O7YpIT8rXw2dms83sRjO7IU3bKzazy8xs/SDWPdfM/jbUY01n7lbeDlR0jyJmttDMbjOzH5nZM2bmZvb+bMfVy93ALdkOIkE64tkALE5DLDnJzM4xsx1mttnM3mhmM6K/Mzez75lZbbTe8Wb2vJldb2aV0bzJwJY0h/Q+4Mfu3hFtXyTvKF8PST7m623RY3matlcE7AFmDmLdFcDCvhaa2SIz+1X0t3etmdUnLEt37lbeRkX3aPO/wBfd/VLgGOBn8QVm9jozOyUbQSXu292fd/dHshFHMqnE09d76O4rgR3pjCudn9dwt+XudwP/BXS4+2PuvgG4DGgD/ubue6L1Hick+U+5e3P08rcBdw7rABKY2f8FvgDcb2YvA/PTtW2REaZ8naJ8zNfu3gQ0pGtf7t4OPDnIdbcSCvS+lr8ALI8mv+nuiUV22nK38nY3Fd2jhJnVAfOAAwDuHgM+B7iZjSMkdMtCXFnbdzoN4jhiI7ivbGzrV8B0MzsWwN13An8EzkvYVy3QEP3txS1w95eHue8u7v5ld5/r7kvcfYG7P5SubYuMFOXrzMrBfJ3uEStSiX+gdTujx979rNOWu5W3u6noHiWi4XdWAvea2RuieeuBxwn9qA4FlpnZBwHM7Etm9n+jn5a6+m2Z2RIz+w8z+4SZNZnZSjP7NzN7u5k9YWbvMLP/Z2b7zOzyhNcl3V6vfV9mZj9N7NsW/cT6XTP7mpk9YGbHR/P73V/C683MLjezF8zsDDN71sy2mtkHBrGPuYnxDLDPg97DJA6Jfi5uNLPLer2nX4jem9uiJD3Y9+ygfZnZhWb2uJl9wMxWmdnX+9less8+aTz9ic6IvAC8O2H2ROBci7qSAOcDdyTEWQk0JUwvMLOvm9k3zexBC91U3mtmK8zsAjO7xcz2mNmVZnZsdIw7zeycgeITySfK14WTr5Pp6xj724+ZHWpmN1k4a3zdYPaTZL93m9n/mNmbBlhPuTtT3F1tlDRgEbCe8D/bm4BJCcvWAUuj50cCa6LnEwn/Cx9P+J/6a8CSaNn3gUej52XALuB7hLE2/wHY1N/2eu+b0KftBuCn0fwKQiFXG01fCGwHxvS3vyTHXR/t8zPRul8BWoBZA+yjdzz97jPxPUwSwzrgeqAqet1r0fyxwC8T1vsLIWEO6j3rY191hJ8MvwUcD5ze3/Z6ffZJ4xnk39e1wOro+WHA14BW4J0Jfy+WsP75wMnR8yrgCaAqmn4oeh/GAnuBLwGVwHuj6fcQTgp8kuhvUE1tNDWUrwsiX0fLvzDI97G/z/oJYGHCa3yQf2fxz/ONwOeTLF8a7ackYZ5yd4aaznSPIh7ORh5FSEKXAC+YWbILRlYB/2hmpcDJ0bwawpd8BhDvj/sS4cuEu7cBjcBv3X0f4UtYP8D2esfXSs8LM84Fmj3qEwzcRviyvnOA/fXebnybv4jW/Qqhv/EZA+yjRzyp7LMPv/PQf+8JYErCMY6PzkR9gvCewiDfsz6Odxchud3u7o+7+/0pbK+veAbjV8AcC11M3gXcCPweeI+Z1QCNHmXlyJuAR6Pn5xP+YYufPXkH8OXofd4D/MlDP/AVwFh3/7WHn9yfIrXPQCQvKF8XRr5Oos9j7Gc/pwPj3D0ez4oU93kccLG7f3GQ6yt3Z0hJtgOQ9DCzYmCGu68DPm5mPyZc6b0cOCFxXXdvNbNpwNXAf8c34e7bzGwF8BbC/8QPJRRaXS9NeN5B1D2pr+31EWriNuYCpQlxxcxsHTC9v/0NxN1bzGwNUAtMSGEfQ95nLx10f7dmAmvd/Vu9VxrMe2Zm/0gobuPme/gZ2hNjTeEz6DOegbj7KjN7itDFZJy7rzGzXxHOPr0buCsh7vjfRry/4CzCman4tnYmbjrheSc9xVA3OBlllK97HF8h5OtEfb6P/Xw2ryP8B2OoTgPmmlmNu+/vY51YdAzK3RmkN2T0KCWhv627P0P4+W5B7xXN7DTgw+7+FQ8jUST6CHCchaGrngX+faAdD7C9/rxGOHOaOJSSEfo6DlcZ8EqG9zEYm+nZ7xkzOyGF9+wO4OiElvQq+BS2lzSewR8OvwQ+CKxJiK8YuJLwU2zcCYT+qXGbgJPMrDphv/32KxQZxZSveyqkfN3nMfazn33Ra8oYmm8SfiX4dh/L93r3BfDK3Rmkont0+ZiZTU2YngbcGz1vI1w4soAwPNU4Mys3s9Oj5ZPMbAKhb+GPCT8NPUr0c2WkiF7/uzczG2B7vfdtCdv4X2A/cFG0rdpof/Ezpn3try9TonWmRa+9dxD7SIxnoH0mHkdvyc56GOHsVTVwp5mdZeGCnYE+g659AVPc/eWE1p6wi+KE5/1tLzHupPGYWZmZLY+6jvTn14SftX8L4O6NwD3AX7znqCVnErqexMXf71ssjPX9USA+Vmvvz6D355zXIymI9EH5moLJ14N9H5Puh5BLywj/MYPwqwZmNnUQudsIZ53/EXi3mb0neu0nzexawln+Pyasr9ydSZnoKK428o1wcYYDu4GbgZ8Sfp46JFr+eWAjcAHh56JVhDML5xLOkPyUkMDuIFzU0RZtr4nwJXwr0E64WGcG8Olo+Yf6216vfV8JPEzoJ7coWraEcIb0ekLfxiOj+X3ur4/j9+j1VxK6OxyVsKyvfcxOjGegfSa+h732fQHhgsLvRe9F/HUfjJafAjwffTbXR+/zYN+zC5Ic64XR/m4BpkXz+vtMe2yrj3imE27i8ItB/K19p9f0e4FTe837ZpLXvRl4EdgJXBvNO4fwt/ZdYCrRsGnAFYR/lG+MjvW8bH/H1NTS1VC+LqR8PRN4hHC2/ugBjrG/z/pcwi+MfyVc1P4K8OHo2JPmbsKvKe2EM9x1wJ8IBf+lwOWE8cqvJ3Rtib9GuTuDzaI3TwQzmwRc4e5fiqaLCGdfLnf3a7Ma3ADMzIFDPfSRlCEys0+5+5CGoxKRkaN8LYmUu/ODupdIoo8SuhrUJcw7jYQL5HJRws9Z+ilriMys2MwuBm7NdiwiMijK16LcnWdUdEuiHxL6fr1kZhuA+4CN7v5YdsPqW3RhyZXR5AcsDF0nqasB/tfd1wy4pojkAuVrAeXuvKLuJSIiIiIiGaYz3SIiIiIiGaaiW0REREQkw1R0i4iIiIhkmIpuEREREZEMU9EtIiIiIpJhKrpFRERERDJMRbeIiIiISIap6BYRERERyTAV3SIiIiIiGfb/Af+zKmg8NwyTAAAAAElFTkSuQmCC\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-neptune.png',bbox_inches='tight')\n", "plt.savefig('../plots/igpd-neptune.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/igpd-neptune.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 }