{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 30 - Huygens - Titan" ] }, { "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 Huygens probe.\n", "https://en.wikipedia.org/wiki/Huygens_(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.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", "vehicle=Vehicle('Huygens', 318.0 , 35.0, 0.00, 5.73, 0.0, 1.25, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(1270.0,0.0,0.0,6.0,0.0,-65.4,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (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": "iVBORw0KGgoAAAANSUhEUgAAAt4AAAHtCAYAAAAuknnBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzde3xcdZ3/8denaZqkF+glaUsLJUABy1UhIMoqCKjsoiwLZRG8FdQg/nRXVxAXES/IroA3BEEiC64gyqL85AdyEW/cBDQFQdACXVoKKYX0RkPTS5p8fn98z6STYZLJJDNzzsm8n4/HeczMOWdmPkk7n3zyzef7PebuiIiIiIhIeY2LOwARERERkWqgwltEREREpAJUeIuIiIiIVIAKbxERERGRClDhLSIiIiJSASq8RUREREQqYHzcAVRCY2OjNzc3F/Wczs5OmpqayhNQGaUx7jTGDOmMO40xQzrjLlXMixcvXu3u6friR2kkORvS9/8kbfFC+mJWvOWXtpjLHW+hnF0VhXdzczPt7e1FPaelpaXo5yRBGuNOY8yQzrjTGDOkM+5SxWxmz5cgnFQZSc6G9P0/SVu8kL6YFW/5pS3mcsdbKGer1UREREREpAJUeIuIiIiIVIAK70G0trbGHcKIpDHuNMYM6Yw7jTFDOuNOY8xpl7bvedrihfTFrHjLL20xxx2vuXusAVRCS0uLp6n/SEQkw8wWu3tL3HFUknK2iKRVoZytEW8RERERkQqoilVNREQqbdu2bQCMH680W27u0NcX7tfUxBuLiMhQNOItIlJip59+OrW1tVx//fVxh1IV/vAHGD8e3v72uCMRERmaCm8RkRKrq6sDYPPmzTFHIiIiSaLCW0SkxBoaGgDYtGlTzJGIiEiSqPAWESkxFd7xqIJFukQk5VR4i4iUmArvyjKLOwIRkeFR4S0iUmL19fWACm8RERlIhbeISIllRrw1ubKy1GoiIklX8cLbzN5pZo+YWXPWvhPNbJmZrTGzy8xsfNaxRjO72swuNrNvmVntcI6JiMRlLLWapCFnq9VERNKiooW3mc0CJgOHZu2bB5wALAQ+A3wU+HTW024GrnX3c4FngK8O85iISCzGSuGtnC0iUloVLbzd/WXg1pzduwJnuPtid/8R8D3gHQBmdhiwt7s/Ep17G/ApM5s01LGyfyEiIkMYK4V32nK2Wk1EJOkq3mri7n05j+93921ZuzqAFdH9o4BlWed2RHdbChwTEYnNWJpcmYacrVYTEUmLJE6uPAT4fnR/LrAm53gXMKfAsQE6OztpaWnp39ra2kocsojIdqOdXNnW1tafr4DG0kVWFsrZIlLVisnZ44c6WGlmNh/odPfHo10O5A4Z1QE9BY4N0NTURHt7e4mjFRHJb7StJq2trbS2tgJgZqtLFliJJS1nq9VEROJQTM5OzIi3mdUAZwKfz9rdAUzLOseAScDKAsdERGIzVnq8h5KknK1WExFJi8QU3sDZwDfcfUvWvnuAvbMeNwNbgfYCx0REYlMNhTfK2SIiRYtjHe/M2IRl7TsPeAxoMLPdzewMM5vv7u3AS2a2X3TqscDl7r51qGMV+lJERPLKTK4cCxfQSVPOVquJiCRdRXu8zWwK8IHo4SIzu4Lwp8oLc05d4u7XRvdPAS4ws+WE9WS/kHXeUMdERGIxVka805Kz1WoiImlhXgVDBC0tLa7JlSJSKRs3bmTy5Mk0NDTQ3d09qtcys8XuXlXLpBabs//4R3jzm+GQQ8J9EZG4FMrZSerxFhEZE7JHvKthcCMp9K0WkaRT4S0iUmLjxo1jwoQJAGzZsqXA2TJaajURkbRQ4S0iUgZjaYKliIiUhgpvEZEyGCsTLNNErSYiknQqvEVEykCFd+Wo1URE0kKFt4hIGajwFhGRXCq8RUTKQIV35anVRESSToW3iEgZZApvTa4sP7WaiEhaqPAWESmDzKomGvEWEZEMFd4iImWgVpPKU6uJiCSdCm8RkTJQ4V05ajURkbRQ4S0iUgYqvEVEJJcKbxGRMtDkyspTq4mIJJ0KbxGRMtDkyspRq4mIpIUKbxGRMlCriYiI5FLhLSJSBiq8K0+tJiKSdCq8RUTKQD3elaNWExFJCxXeIiJloBFvERHJNb7Sb2hm7wS+Bpzi7sujfQ3AJUAXMA041903jOaYiEicxsrkyjTlbLWaiEjSVXTE28xmAZOBQ3MOXQk84O7nAbcCV5fgmIhIbMbCiHdacrZaTUQkLSpaeLv7y4Rk28/M5gCnAbdFu34NHG9mzSM9VtYvQkRkGMZC4a2cLSJSWhVvNXH3Phs4PHEksNrdu6Pj28xsGXAE0DPCY8sr9OWIiOQ1ViZXpilnq9VERJKu4oV3HnOBNTn7uoA5wLYRHhugs7OTlpaW/setra20traOLmoRkSGMZsS7ra2Ntra2zMPG0kVVEonL2Wo1EZE4FZOzk1B4O5D7k6mOMDrCKI71a2pqor29fZRhiogM32gmV2YXmma2uqSBjZ5ytohIlmJydhIK7w7C7PZsU4CVhAQ/kmMiIrEaCz3eg0hszlariYgkXRLW8f4dMCdaZgozqwXmAfeO4piISKzGcOGtnC0iMkIVL7xt+ywdA3D3VcDtwNHR/mOAn7t7x0iPVeDLEBEZ0sSJE4H0F95pyNnq8RaRtKhoq4mZTQE+ED1cZGZXuHsncBZwiZkdCMwGsmfRjPSYiEhsMoV3d3d3zJGMnHK2iEhpVbTwdvcu4Kpoy96/DvjYIM8Z0TERkTiNhcI7bTlbPd4iknRJ6PEWERlz6urqMDO2bNlCb29v3OGMaWo1EZG0UOEtIlIGZjYmRr1FRKR0VHiLiJTJpEmTABXelaJWExFJOhXeIiJlohHvylCriYikhQpvEZEyUeEtIiLZVHiLiJSJCu/KUquJiCSdCm8RkTLJFN4bN26MOZKxTa0mIpIWKrxFRMpEkytFRCSbCm8RkTJRq0llqdVERJKuJIW3mf1DKV5HRGQsSWrhPdZytlpNRCQtirpkvJkdC5wPzGZ70T4OmAU0lDY0EZF0i7vHWzlbRCRZiiq8geuB/wQeA/qifQacVsqgRETGggSMeFdVzlariYgkXbGF93Pu/q3cnWa2pETxiIiMGQmYXFkVOVutJiKSFsX2eJ9lZgvz7D++FMGIiIwlCRjxVs4WEUmQYgvvw4FrzKw3a+sDripDbCIiqZaAwruqcrZaTUQk6YptNfkq0Aq0A73RPgM+WMqgRETGgrgnV1IlOVutJiKSFsUW3ncCt7r7luydZvaD0oUkIjI2JGDEWzlbRCRBii28fwr8u5n9NmufAScC/1qyqERExoAETK6sqpytVhMRSbpiC+8vAnsCi4BMijPCGrGjSuJmtgD4JLA0eo82d/+zmTUAlwBdwDTgXHffED1n0GMiInFLwIh3VeRstZqISFoUW3hfCtzp7l3ZO0t0FbTrgX909w4zmwfcDSwArgTucvebootBXA2cGj1nqGMiIrFKQI+3craISIIUu6rJC7kJPDK7BLEsAKZE9zcBO5rZHMKFHm6L9v8aON7Mmoc6VoJYRERGLQEj3lWVs9VqIiJJV2zh/YncHWZ2MvCdEsTyM8KyV1OADwCfAo4EVrt7N4C7bwOWAUcUOCYiErsEFN5VkbPVaiIiaVFs4T3dzL4IYGYTzewa4BrgTyWI5RPA5ui1XnP3nwNzgTU553UBcwocG6Czs5OWlpb+ra2trQThiogMbaSTK9va2vrzFdA4ihCUs0VEyqyYnF1Uj7e7H2dmh5nZ5cC7gVXAgcCrIw02SwPwE0Jy/o6ZLSNMBtqUc14d0BPdH+pYv6amJtrb20sQoojI8I10xLu1tZXW1lYAzGz1SN+/2nK2Wk1EJA7F5OwhC28ze/sgh54HngW+AewKnAB8puhIB7oBeJ+7rzczJyT0TxNmvmebAqwkJPjBjomIxK7SkyurNWer1URE0qLQiPclhAk064G+PMevJSxNtROjSOJm1ggc6O7ro11fIyx/9Twwx8wa3H2TmdUC84B7CVdh+8Egx0REYldfXw/Ali1b6O3tpaamptxvqZwtIpJghQrvrwLt7v7KUCeZ2cJRxrEW2Gxmc929I9q3BngCuB04Oro9Bvh55hwzG/SYiEjczIyJEyfS3d3Npk2bmDx5crnfsqpztlpNRCTphiy83f2O4byIu/9sNEG4e5+ZnQBcYGaLgVnAOe6+wczOAi4xswMJS2C1Zj11qGMiIrGbNGkS3d3ddHd3l73wrtacrVYTEUmLYi+gUzbu/jhwZp7964CPDfKcQY+JiCRBAi6iUxbK2SIixSt2OUERESlCAtbyrhpqNRGRpFPhLSJSRiq8y0+tJiKSFqMuvM3sTWZ2g5l9VpdrFxEZKGmFt3K2iEh8SjHi/XHgGHf/JvBWM2sowWuKiIwJI716ZRmN2ZytVhMRSbpSFN6fBA6O7v8UqC3Ba4qIjAkJnFw55nK2Wk1EJC2KLrzN7O/N7P3R/UOBPTPrsLp7n7tvKHGMIiKpFXeriXK2iEhyFFV4m9lFwM3APwG4+x+BD5nZaWWITUQk9eIsvKstZ6vVRESSrtgR77cCOwP3Z+27E/h6ySISERlDYh7xroqcrVYTEUmLYgvve9x9PZA9rjDaSw+LiIxZMfd4K2eLiCRIsVeufM7MzgPmmtn7gJMJf8L8VMkjExEZAzKF96ZNm+J4+6rK2Wo1EZGkK2rE291/CvwWqAE+BKwH3uXu3ytDbCIiqdfQEFbri6PwVs4WEUmWYke8cfeHgYczj82sxsyOdvfflDQyEZExIFN4x7WqSTXkbPV4i0haDFl4m9n/AyYXeI1GYDUwZpK4iEipVLLVRDlbRCTZCo14PxGd8zTQB3wMeAxozzpnN2B6WaITEUm5CreaVGXOHh/9JOvpiTcOEZFCChXelwIWzYrHzPZw9wtyTzKzW8sRnIhI2lW41aQqc3ZdXbjdsiXeOEREChmy8Hb3V3N2zTKzGnfvzewwsyOBQ8sQm4hI6lWy1aRac3Z9fbjdvDneOERECil2cuXPgMfN7AGgC1gAvBt43YiKiIjEu6oJVZKzM4W3RrxFJOmKXU7wHuBdhIk5exOWplro7v9ZhthERFIvzlVNqiVnjx8PNTXQ2wvbtsUdjYjI4EaynOBK4PzMYzMbZ2YHu/viUgRkZocBfwcsBR4ANgKXEEZrpgHnuvuG6NyGwY6JiCRBzBfQqZqcXV8PGzeGdpPJhdZ1ERGJSVGFt5ldm2d3E9ALnDDaYMzsdGAPd8/+IXEdcJe732RmxwJXA6dGh68c4piISOzibDWpppytwltE0qCoVhNgT2AF8HzWVgeMeuTEzI4A3g98MWvfHOA04LZo16+B482seahjo41FRKRUYu7xrpqcnenzjukPCyIiw1Jsq8lH3f3p7B1m9gZgUQli+RbwB+ByM5sPXAjsCqx2924Ad99mZsuAI4CeIY4tL0E8IiKjlmk1ienKlVWTs6dNg44OWLsWdtlltK8mIlIexU6ufDrP7i3AR0cThJntDbwJ+IG7f5JwRbW7gF2ANTmndwFzgLlDHBugs7OTlpaW/q2trW004YqIDFv2iLe7D+s5bW1t/fmKcKXJEammnN0YfZfW5L6DiEiZFZOzi+3xXgZk/+SoAWYB1xcbZI59gbXu/kT0+ArgS9F75f7hsI4wckKBY/2amppob2/P3S0iUnbjxo1jwoQJbN26lS1btlCf6YkYQmtrK62trQCY2eqRvnc15ewZM8Lt6hF/t0RERqaYnF1sq8k9wI1sT+R9wMvu/kyxQeaJoz8Wd99kZs8SfkhMyzl3CrAyimGwYyIiiTFx4kS2bt1Kd3f3sArvEqqanJ0Z8VbhLSJJVmzhfYO735e9w8z2NrP57r50FHE8AexoZrPc/eVo3zZCQp5jZg1RYq8F5gH3Embl/2CQYyIiidHQ0MD69evjmGBZNTl7p53CbUdHKV5NRKQ8il3V5Ng8+14Gbh5NEO6+BLgTWAhgZlOBCcCPgduBo6NTjwF+7u4d7r5qsGOjiUVEpNRiXNmkanL2/Pnhdulofp0QESmzYY14m9knga8BU8zs3Dyn3FOCWD4EXGZmEwijIKdGs97PAi4xswOB2UBr1nOGOiYikgiVXtmkGnO2Cm8RSYNhFd7ufoWZ3Qd8GPhuzuGt7v7SaANx99WENWFz968DPjbIcwY9JiKSFLW1tQBsq9D1zKsxZ++5Z7hdsiRcOr6mphzvIiIyOsPu8Xb3J8zs8+7+uhnoZjbP3VeUNjQRERmpasvZ06fDbrvBsmXw1FNwwAFxRyQi8npDFt5mdjZwh7v/Ndr1OTPLHUcw4B3AkaUPT0REhqvac/Zb3hIK74ceUuEtIslUaHLlbsDUrMcHAW8Fdo+OZbbZZYlORESKUdU5++/+Ltz+6lfxxiEiMpghR7zd/f/k7LoIeNLdt2bvNLN9Sh2YiIgUp9pz9nveA5/4BNx1F2zaBNFiMiIiiVHscoL750ngewN7lC4kEREpkarK2bvsAi0t0N0Nt94adzQiIq833OUE30boC3xndAnibNMJs+ZvK3FsIiJjQk9PmN84fnyx1ywbmWrO2aefDu3tcMUV8L73xR2NiMhAw/0psB74EaE38PCcY1uBa0oZlIjIWNLV1QXAlClTKvWWVZuzP/Qh+Pd/hwcfhPvvh7e9Le6IRES2G+463n8xsyOBk919zCZsEZFyqHThXc05e/Jk+Mxn4CtfgXPOCSucmMUdlYhIMOweb3d/dbAEbmZ/V7qQRETGlhhGvKs6Z599NsyeDY88AldeGXc0IiLbFVrH+/8Bkwu8hgH7AjNLFZSIyFjx2muv0dPTQ319PXV1dWV9L+XsYPLk0OO9cGEowt/+dth//7ijEhEp3GryNKFXsAPwQc4x4D2lDEpEZKx44YUXANhll10q8XbK2ZGTToJFi+CHPwzLDD78MOy0U9xRiUi1K1R4Xwx0ufuWoU4yszE5O15EZLRWrAhXZp83b14l3k45O8uVV8LTT4c+73e9K1xYR8W3iMRpyB5vd19dKIFH/qlE8YiIjCnPPfccUJnCWzl7oIaGsJ73ggXw5JOh5eTZZ+OOSkSqWbEX0BnAzOaZ2cXAZSWKR0RkTHn00UcBOPDAA2OOpDpzdlMT3HcfHHQQLF0aLrCji+uISFyKLrzNrMbM/snM7gD+F3g/oadQRERyLF68GICDDjoolvdXzobGRvj978Nkyw0b4IQT4GMfg/VV9V0QkSQYduFtZruZ2UXAC8BPgfnAe4FdgJPLE56ISHqtWbOGxx9/nNraWt70pjdV9L2VsweaMgX+53/gm9+ECRPgmmtg333hhhugry/u6ESkWhQsvM3sZDP7FfAscBLwLWBn4JfufpcHfyhznCIiqXPXXXfR19fHEUccweTJhVb5Kw3l7MGZwb/9G/z5z/CWt8DKlfDBD8LBB8Pdd4MPtg6MiEiJDGfEexPhEsN/AI5392+4eyeDL1UlIiLATTfdBMB73lPR1fuUswtYsCBcTv6662Du3FCIH3tsKMB/8hPYti3uCEVkrCpYeLv77e7+HuADwAfN7HozWwjUZs4xK80Fec1sgpk9Hl3qGDNrMLPLzew/zOwqM9sh69xBj4mIxG3FihX88pe/pLa2llNPPbVi71vJnB29Virzdk1NWOf72Wfh4oth1ix47DE47TTYfXf48pfh+efjik5ExqpiLhm/wt2/CJwe7drZzC40s2bgYyWK5xygOevxlcAD7n4ecCtw9TCPiYjE6tJLL6Wvr4+TTjqJmTMrf5HICuVsSHnebmiAz30Oli+HtjbYay944QX4yldgt93C+t///d+aiCkipWE+iqY2M5sPfAr4mLtPHFUgZm8F3gBcACwCngGWAdPcvdvMxgOvEi51vHWwY+6+PPe1W1pavL29fTThiYgM29KlS1mwYAG9vb08/vjj7D+K65Wb2WJ3bylFXKXM2dHrlSVvx5mz+/rCCijXXAO33AJbolXRa2vhmGPg5JPhuOMght+lRCQFCuXsUa3j7e5L3f1fgY+P5nXMbBKw0N2vzdp9JLDa3buj99pGSNpHFDgmIhKb3t5ezjjjDLZt28aHP/zhURXdpVaqnA1jN2+PGwdHHQU33hgmX151VXjc2wt33glnnBHaUlpa4Pzz4YEH1BMuIsM3qsI7w91/NMqXOBf4es6+ucCanH1dwJwCx16ns7OTlpaW/q2trW2U4YqI5HfRRRdx//33M3v2bC699NIRvUZbW1t/vgIaSxogJcnZUMa8nZScPX06fPzj8JvfwEsvwdVXw7vfDXV1sHgxXHQRvO1tMGMG/MM/wNe/Dg8+uH2UXESqQzE5e3xlQhqcmb0baHf3V3IOOWF2frY6oCe6P9SxAZqamlCriYiU249//GO+9KUvYWZcd911NDaOrGZubW2ltbUVADNbXcoYS6HceTuJOXvmTGhtDVt3d7ga5l13he3pp8No+J13hnPr6+Gww8KShS0tcMghsPPOYTlDERl7isnZsRfewNnAm7Im2U8jTLr5VnQ/2xRgJSG5D3ZMRKTibrjhBhYtWgTAt7/9bY499th4Ayqvqs7bEyeG5Qcz/8QvvhiWJ7zvvrD99a+hT/z3v9/+nJkzQwHe0hK2Aw9UMS5SjUY1ubIkAZg1EUY9Mh4C/g14BPgr0OTum8ysFngN2B3oBZbmO+buHbnvocmVIlIufX19XHzxxZx33nkAnH/++Vx44YUle/1STq4slXLn7bTn7NWrQ+/3n/4E7e3hdt2615+3446w336w//7bt/32g2m5v56ISGoUytmxj3hHF3boZ2a9QKe7rzCz24GjgduBY4CfZxL0UMdERCqhs7OTj3zkI9x2220AXHLJJZxzzjkxR1V+yttDa2yEE04IG4QrYi5btr0QX7wY/vKXUKA/+GDYss2dC3vvHZY23Guv7febm2F87D+1RWQ0kv4RPgu4xMwOBGYDrcM8JiJSNu7Oj3/8Yz796U+zZs0apk2bxvXXX89xxx0Xd2hJoLydwyxclGf33eGUU8I+d3j5ZXjyyVCEZ7annoKOjrD99rcDX6e2FvbYIxThe+4ZXm+33cK2665hTXIRSbbYW00qIe1/thSR5HjggQc499xz+cMf/gDA0UcfzTXXXENzc3NZ3i+JrSblVs05u7c3XMzn2WfDpM1nngnb00+HC/sMZfbs7YV4c/P2+7vtFkbR6+qGfr6IjF7iW01ERJLO3bn33nu59NJLueOOO4Cw8sbFF1/MokWLKOEV2KXK1dSEUe099tg+eTOjuxuWLg1F+LPPhvaVZctCof7887BqVdgeeij/a8+cGSZ07rJLuM3d5s7VqLlIuanwFhEZRHd3Nz/72c+47LLLePTRRwGYNGkSZ599Np/97GeZMmVKzBFKNZk4EQ44IGy5entDe0qmEM8U5ZnHK1fCK6+ELfqvnNeMGaEwnzsXdtopjKJntuzHE0d93VOR6qTCW0Qki7vz0EMPcd1113HTTTfR1dUFhBHuT37yk5x11lk0NTXFHKXIQDU1MG9e2I7Icy3Q3t7QU/7ii2F74YXt9zNbRwesWRO2P/956PfbYYeBRXluYT5zJjQ1ha2+vjxfs0gaqfAWkarX29vLAw88wC233MItt9zCiy++2H/ssMMO46Mf/Sjvf//7qVcFISlVUwNz5oTt0EPzn9PXB52d2wvxTOvKSy9tv595vGFD2J55pvB7T54cCvDsYjx3yz6m0XQZy1R4i0hVWrlyJffccw+/+tWvuOeee+js3L5C3ty5cznttNM4/fTTWbBgQYxRilTOuHEwa1bYDj548PPcYf36oQvzzs7Q1tLZCa+9FrZly4YXx8SJ24vwxkaYPj20wAx1u+OOIX6RpFPhLSJVYdWqVTz44IM88MAD3HPPPTz11FMDjs+fP58TTzyRE088kUMOOYRx+ikukpdZuMjPtGlQ6PdSd3j11VCAZ7ZMQT7Yvu7uMFn0+eeHH9O4cSGeQgV69u20aTBligp2qSwV3iIy5vT09PC3v/2Nhx56iAcffJAHH3yQ5557bsA5kyZN4sgjj+Rd73oX73znO3nDG96g1UlESswMpk4N2557Fj7fHbq6thfhmZ7ztWuHvu3q2n5uMcaNC6PlU6eGQjwTa+Z+oX1aBUaKpcJbRFJt8+bNPPnkkzz66KP92xNPPMGWLVsGnDdp0iQOO+wwDj/8cI466ije8pa3MGHChJiiFpF8zMLEzR12CEsqDldPTyjCM9twivV160ILzLp1YRtuK0y2urrii/ZMob/jjuGiSFJdVHiLSCr09PSwdOlSnnrqKf7617/23y5ZsoRt27a97vw99tiDlpYWDj/8cA4//HAOOOAAxut62yJjUm3t9v70YmzbFlph1q0LfeuZ2+z7g+1btw62bAmrxbz88sjinjRpe3Geu2UK9KGOa+wgffRTSEQSw91Zs2YNS5cuZenSpTz77LMsWbKEp556imeeeYaenp7XPcfMWLBgAQcddFD/9sY3vpGpU6fG8BWISJqMHx/6vWfMKP657rB5c/FF+/r1odhfvx42bgxbR8fI4p84sXCBPlQhr6uZVp4KbxGpqL6+PlatWsXy5cv7C+xMkb106VLWr18/6HObm5vZd9992WefffpvFyxYwOTJkyv4FYiIhLaYhoawzZlT/PPdQ6tLphgfbMsU6fm27u6wvfTSyL6G+vriR9qzz6uvD98HGT4V3iJSUj09PXR0dLB8+XKef/75120rVqxg69atgz5/ypQpzJ8/nz333JP58+ez1157se+++7JgwQImTZpUwa9ERKR8zMKqKlOmhKuFFss9jJYPp0DPV8ivWxdG7DPLQI7EhAkjG2nPnpxabYW7Cm8RGbZNmzbR0dFBR0cHK1eu7L+f2V544QU6Ojro6+sb8nUaGxvZddddmT9//oAie88996SpqUmri4iIFGAWLk40eTLsvHPxz3eHTZtGPtq+bh1s3RqWgnzllZF9DbW1Ixtpz9yfNCl9hbsKbxFh48aNvPzyy/3bqlWrBhTUmSJ73bp1BV/LzJg7dy677rrroJtGrkVE4mUWesQnThx5q8zmzcWNtueeu3nz9qUjR6KmZiAZzgsAACAASURBVGBRPnPm9km2+bYZM+Jft12Ft8gYlV1Mr1q1akBhnbvvtddeG9Zr1tbWMmfOHObOndt/m73tvPPO7LLLLlqmT0RkjMvucd9pp5G9xubNgxftwynmN20qbv32mppwRdRddoHm5rDtumu43W23sNZ8uZd4VOEtkgJ9fX28+uqrdHZ2snr1alavXj3k/ZdffpmNGzcO+/UnTJjArFmzmD17NrNmzWLWrFkDCupMkd3Y2KgrOoqISEnU14et2GUgM7ZsCQX6q6+G9dlfeWX78o75tnXrtve0/+lP+eM54AB485vhQx+ClpbRfX35qPAWqTB3Z8OGDaxdu5a1a9eyZs2agoX0mjVr6O3tLep96urq+ovo7II6374dd9xRfdUiIpIqdXWhvWTmzOGdn+lJX7ECli8P2/PPh9tnnw0XUfrjH8N2xRXw61/DUUeVNmYV3iIj1Nvby6uvvtpfPGcK6ULbunXrii6iAXbYYQeamppobGyksbFx0PuNjY3MmjWLHXbYQcW0iIhUNfew5GL2VUs7O8MSjKtXb19PffPm1z9viAW4RiwxhbeZnQh8E9gBuAH4rLtvM7NG4CJgPVALnOvuPdFzBj0mMhzuzsaNG1m/fj3r1q3rv81sQxXQ69evx91H9L5Tpkxh+vTpTJ8+nWnTpvUXzvmK6aamJmbMmKG+aUkU5WwRKbfMyitdXbBhw/C3V18dWGhv2TK896uvh733hre+FU45BY44ovRfUyIKbzObB5wALAT2Ba4CXgC+AdwMfN7dHzGzjwNfBf49eupQx6RK9Pb2sn79+gFFc24hPdSxfJcbH66pU6f2F9DD3aZNm6YiWlJNOVtE+vrCKPGmTWHL3M/cbtwYRpozV+cczuPcfd3dpYm1vj6saDJ9erhtbIS5cwdu8+aFiZblnsaUiMIb2BU4w923AYvNbD/gHWb2ALC3uz8SnXcb8LSZfQ3Yf7Bj7j78WWUSO3dn8+bNIyqa161bR1dX16jev6GhgalTpzJt2jSmTZvG1KlTX1dQz5gx43UF9NSpU6mpqSnRd0EkVZSzRRKgry+M5ma2zZtff3/z5tcXxYPdL3Q8+3452jDyqa+HHXYIFxraYYfhb9mF9sSJlYl1OBJReLv7/Tm7OoApwFHAsqzzOqKe1Rbg8CGO3VvmkCWHu9PV1dXfw5xbKBcqpIe6kmEhZsaOO+7YXzRnF9D5bnP31dXVlfA7ITL2KWdLNcv0/uYWt0MVv0MVxaM5tyfmRq2GhlAYZ5YVzL4/cWK4wE3mNnvL3TfY44kT4193u9QSUXjncQhwKdAK5K7O2AXMAeYOcWyAzs5OWrLWhGltbaW1tbWU8Y4ZW7Zs6e9tzr3Nty/7diQTBjPq6upGXDjvsMMOWuJOxpS2tjba2toyDxvjjGWYlLOl7LZtG7oIrcS+zJYkdXXbt/r6/Pdzi+KhCubh3q+rS99VI8ulmJyduMLbzOYDne7+uJk5sCnnlDqgBxjq2ABNTU20t7eXI9xEcndee+01Vq9e3b/axlAFc3ZR3T2KhqrJkyf39zBnb8MpoOvr60v4HRBJt+xC08xWxxzOkJSzx7bsVoa4i96+vri/G9tNmDB0oVup+xMmqPhNgmJydqIKbzOrAc4EPh/t6gD2yjpuwCRgZYFjY0pfXx/r16+ns7Ozf33nzP3Bti0j/JV8/Pjx/cVz9mTAQvumTp2qCYMiVUY5u3Iy7Q3d3aG/djS3xRS9cbcyZKupGVh05t6WY1/mcfa+CRPGXvuDVE6iCm/gbOAb7p6pGu8BPp51vBnYCrRHt4MdS4Wenh5WrVrFypUr824vvfQSr7zyCqtXry66jaO+vr5/GboZM2YMq3iePn06kyZN0trPIjJcVZWzh6OnB157LSx/1tU1+P3sx5nVGwoVzCNcvXTUSlXIjuY16upgfNIqFpERSMx/YzM7j5CAG8xsd+BI4D7gJTPbz92fBI4FLnf3rUC7mQ12LBE2btzIsmXLeO655/pvn3vuOV588UVWrlzJK6+8MuzXylw8JXvLrPGcb5s0aVIZvzIRqXZjMWdDKG67usKFNdauDZeYXrt28Pvr1g0sosvZ/1tbGyabZSaujeQ205873AJZrQwipZWIwtvMzgcuzNm9xN2vNbNTgAvMbDkwGfhC1jlDHauYjRs38uSTT/LEE0/0b0uWLClYWI8bN47Zs2czZ84c5syZw0477dR/P/N45syZNDY2auUNEUmMNOfsLVtg6VJ4+ml47jlYuTJsL720/f5o1g4eNy4sezZlCkyePLz7mZUcChXMGvEVST8b6ZX30qSlpcVLOVGnr6+P++67j7vvvpvf/va3LF68OG8ryIQJE2hubmb33Xdn9913Z7fddmP33Xdn3rx5zJ07l5kzZ2odaBEZkpktdveWwmeOHaXM2UuXwq9+BffdB3/6EyxfXniS3sSJMHMmTJsW1gGePn3w+1Onbi+ip0wJo8UaIRapXoVytn5/LkJPTw9XXXUV3/72t1m+fHn//vHjx3PAAQf0b/vvvz/77LMPc+fOVWEtIhKDRx6Bz30uFNzZxo2DPfYIl4WePx923hnmzAnbTjuF2ylTVDyLSHmo8B6mtWvXctxxx/Hwww8DMG/ePE455RSOPvpoDj/8cCZPnhxzhCIiAnD33fCe94R1nydNguOPhyOOgLe+FfbaK/Qui4jEQYX3MLg7H/zgB3n44YfZZZdduPzyy3nve9+ri7aIiCTQOeeEovvMM+HSS8MItohIEqjwHoY77riDO+64g6lTp/LAAw8wb968uEMSEZE8Nm2Cv/wlrMbx3e+GWxGRpNCQbQHuzoUXhsn7X/ziF1V0i4gk2Nq14XbGDBXdIpI8KrwLeOyxx3jkkUeYPn06Z555ZtzhiIjIEDJXWlTRLSJJpMK7gGuvvRaAD37wg7oojYhIwmUKb615LSJJpMJ7CFu3buXGG28E4PTTT485GhERKWTbtnBbWxtvHCIi+ajwHsJvfvMb1q1bx3777ceBBx4YdzgiIlKARrxFJMlUeA/hZz/7GQAnn3xyzJGIiMhwaMRbRJJMhfcgenp6+MUvfgHAwoULY45GRESGQyPeIpJkKrwH8bvf/Y61a9eyzz77sM8++8QdjoiIDINGvEUkyVR4D+Lmm28GNNotIpImfX3hVhcWFpEkUmoaxB133AHASSedFHMkIiIiIjIWqPDOo7Ozk5UrVzJ58mT222+/uMMREZEiuccdgYjI66nwzuPpp58GYJ999mGc/l4pIpIaZnFHICIyOFWVebz44osAzJs3L+ZIRERkJDTiLSJJpMI7j0zhvfPOO8cciYiIFEMj3iKSZKlf6dTMGoBLgC5gGnCuu28YzWs2Nzdz/PHHc/DBB5ciRBERiZQjZ2ebPh2OPx7e8IZSvaKISOmkvvAGrgTucvebzOxY4Grg1NG84MKFC7WMoIhIeZQ8Z2dbsABuvbVUryYiUlqpbjUxsznAacBt0a5fA8ebWXNcMYmISH7K2SJS7VJdeANHAqvdvRvA3bcBy4Aj4gxKRETyOhLlbBGpYmlvNZkLrMnZ1wXMyd7R2dlJS0tL/+PW1lZaW1vLH52IyAi0tbXR1taWedgYZywlppwtImNOMTk77YW3A5ty9tUBPdk7mpqaaG9vr1hQIiKjkV1omtnqmMMpJeVsERlzisnZaW816SDMis82BVg52hfO+s0lVdIYdxpjhnTGncaYIZ1xpzHmCihbzob0fc/TFi+kL2bFW35piznueNNeeP8OmBMtT4WZ1QLzgHtH+8Jx/8OMVBrjTmPMkM640xgzpDPuNMZcAWXL2ZC+73na4oX0xax4yy9tMccdb6oLb3dfBdwOHB3tOgb4ubt3xBeViIjko5wtItXOPOXX1TWzaYSLMSwHZgP/7u6v5ZzTCTxf5Es3AmnsrUxj3GmMGdIZdxpjhnTGXaqYd3X3phK8TiKUMWdD+v6fpC1eSF/Mirf80hZzueMdMmenvvAWEREREUmDVLeaiIiIiIikhQpvEREREZEKUOGdw8wazOxyM/sPM7vKzHaIO6bhMrN3mtkjabn8spmdaGbLzGyNmV1mZqlYV97MjjCzp81snZl9J+54imFmE8zscTM7Mu5YhsvMvmZmHm2Pxx1PMczsMDM728xOMLOxdCGcRElb3k5x7ktj/kjNZ9DMFpjZ98zsM2Z2pZm9Me6Y8slXayT5MzhIvLF9BlPxYa+wK4G73P0mMzsWuBo4NeaYCjKzWcBk4NC4YxkOM5sHnAAsBPYFrgJeAL4RZ1yFmNmOwOHAYcBbgV+Y2e3u/ut4Ixu2c4DmuIMYLjObRJgI885o10gm3MXCzE4H9nD38+OOpQqkJm+nNfdF0pY/0vYZvB74R3fviP6f3A0siDmmAYaoNRL5GcwXb9yfQU2uzGJmc4BlwDR3745+A3oV2Nfdl8ca3DCY2TigF9gt6fGa2duAh9x9W/T4EsL3+bh4IxtatP7wZo8+OGa2GDjb3X8Xb2SFmdlbgTcAFwCL3P338UZUmJl9Orr7fXffHGswRTCzI4AvAu90JdmySlveTnHuS1X+SONn0Mw2Age7+xIzawIed/c5cceVK7fWSPpnME+8sX4G1Woy0JHAanfvBoj+UZYBR8QZ1HC5e1/cMQyXu9+f+U8f6QBWxBXPcLn7pqyiexLwN+D3sQY1DFGsC9392rhjGS4zM+B9wDeBl80s9tGTInyL8H/jcjO7y8wOjzugMexIUpS305j70pg/SOdn8GfANWY2BfgA8KmY48krT61xJAn+DObGG/dnUIX3QHOBNTn7uoDE/cY5Bh0CfD/uIIYr+o35TqAWmBhzOMNxLvD1uIMohgeHEVpNLgd+bGZ/H3NYBZnZ3sCbgB+4+yeB3wB3RX/ylNJLe95OQ+5LVf5I8WfwE8Bm4E/Aa+7+85jjGS59BougwnsgBzbl7KsDemKIpWqY2Xyg093TNHFuKfBDwhX4Lo03lKGZ2buBdnd/Je5YRsLd10U9mv8JfLrQ+QmwL7DW3Z+IHl8BGPCP8YU0pqU2b6ch96U0f6T1M9gA/AS4EfiOmR0TczzDpc9gETS5cqAOYFrOvinAyhhiqQpmVgOcCXw+7liK4e4vAdeaWR9hwlGSnQ28KXRuAOH/+K1mdpG7XxJfWEX7Lilo6yHk1f7c6u6bzOxZYEZ8IY1pqczbKcp9acwfaf0M3gC8z93Xm5kDPzGzZnffGHdgBegzWASNeA/0O2BONIEOM6sF5gH3xhrV2HY28A133xJ3ICPUTkg6SXYa8MasbSXwUZL/5+1cfcCjcQcxDE8AO+b8WXsb8NeY4hnr0pq305L70pg/UvcZjJY6PNDd10e7vgZsIGGrmgxCn8EiqPDO4u6rgNsJ7QMAxwA/d/ekF1ZA/2Q0CH9SSzwzOw94DGgws93N7Izozz6JZWYTzexNWbv+AbgsrniGw9073f3FzEaY3d3p7hvijm0oZjbHzE4zs3HR/+3PAl+IO65C3H0Jof9/IYCZTQUmAL+MM66xKo15O025L435I6WfwbXAZjObm7VvDfBMTPEMKrfWSPpnMF9tFOdnUK0mr3cWcImZHQjMBlpjjmdYsmZBAywysyvcvTPOmIZiZucDF+bsXpKCWfP7Ab80syXAw8Af3T3JyTzNphEmdJ0HPAB8JwlLUw3Th4DLzGwCYeTn1JxZ9FJaqcnbKc59aZOqz6C795nZCcAF0TK1s4BzkvYLzhC1RiI/g/niJbSXxPYZ1DreIiIiIiIVoFYTEREREZEKUOEtIiIiIlIBKrxFRERERCpAhbeIiIiISAWo8BYRERERqQAV3iIiIiIiFaDCW2QEzOxoM/uDmR1Z4Lz9zOx/zOyLJXzvS0v1WiIicTCzA8zs5lLmxkHe53Yze18536NczOwYM3uo0M8ZSRddQEdSzcwWABcRrvp1CHAA8EF3v6HMb/0CcOAwzltBuBBMTSneNLoUb08pXktEBMDMDgH+E3g78CtgGeEiKLXAVe5+dxnedimwIyXKjUO4EniyHC8c5eMTge8TrjD5UHSoDjgJOHSUF/0a7s8ZSREV3pJ2vwD+2d0fN7NxwHWVeFN3f8bM1gzjvA1m9lIJ3/oI4L4Svp6IVDl3/5OZ/RR4g7u/J7PfzI4D/q+Zne3uV5b4PbtLnBsBMLP3An/JFLzufkep3yPD3XuAm8zsU8Dv3f38rDhuHclrmtlM4B3ufpO7P21mq0sUriSEWk0ktcxsOrAXsBHCJXeBC4BKXY61r8TnDcc7gN+V8PVERABedzl1d/8l8GXg22Y2uwzvWdJcbWY7E0a4Ky3fpeh/DXQU8yJmVgf8CGjI2l3Knx+SACq8JbXcfS3wNHCXmR0W7Xse+KOZTTWzC83sj2Z2vJmtMrO/mdlBmeeb2alm9gUz+72Z/UfW/hYz+7KZ3WRmt5jZjtH+ejP7npldYGb/RfgzaVHM7HNRTP9sZoeY2S+i17swivEmM2uK9r9qZl/KeYkJ7r4leq3Pmtm/mdmDZvb1or+BIiKF/Reh5eRkADPbw8y+ambXmNlvomIXM6uLctkXzOx+Mzsl8wJm9gkz+6KZPWxmnxzsjfLlZDNbGOXMRWb2rJl93czmmNl1UT6938yOil7iNGBn4FwzOyGai/NrM/tw1nt8zMwuMrMfmdl/m9kkM9vRzL5kZu1mdpSZPW1mz5lZ80i+YdH35MhoRBwzO9HMLjOzb5rZXWY2OfoaLjOz26KfNy9F3+N9gJPN7N+yXnJa9LOoy8w+Nsh7jjOzi82sNXoPN7P/m/n3kQRxd23aUrsB+wLPE0YFrgVmRvtrCElsLaHXbipwD/AU4RfOtwDnR+dOI4yaHwfsAPw06/UfBC6N7n8baI3uzwJ6Ccm1UIw/JIwazST0Ak6M9tcBdwO3AnMIPzB6gC8BE4GDo8eZ8zP965mv+8ro/hTg83H/W2jTpi29G7AIeHGQYy8D343y6i1ATbT/RuDm6P63gBOj+2cAndH9U4EPRPf3inL1/tHjHwJfju4PlpOnA+uB7wCHAscA3wAuiM49B7gtK1YHmqP7DcDfgEXR43cDN2adeztwNaHt9p+inxcnRl/nXcBFw/ze/T76WfF14FJCv/cxWcdXA/Oi+08RfibVAV8AXgTeDHwkyvu/z8Qbnb8c+GZ07FRgxSAxfBS4Pbo/K/o+vyvu/1faXr+px1tSzd2fMrP9ga8BnwDea2bHeOj57gQ2uPvPAczsy8ADwO6EHzJ9Zvbp6KXuJiT444AZWfv/Fj13MnBW9Fzc/WUze6GIUBuBy4GPuHt39BpbzGwVsMzdV0bvswq410P/42OEHwiNhEmaxxF+SAB0Ax8ys2eAqwi/dIiIlIMTitFDgV2BT5kZwBpgRwsPziQMGgD8N6HVAkKubc/JtTsBf8l5j0XkycnuvtbMXgVudfc/ApjZCmCjme1EGCGenDdo901m9krWrk9Gr5txDfA/0f71hJ8Xt0Tv0R7FOVy/86jH28xuIgz2ZJzs7ivM7HBCwT05yv8rgf9190eAR6Ln5nvt26KfCX8aIqaDgE3R1/2yma1lYMuKJIQKb0ktM6sBdvEwieZfovaPO4A2wghCrqei26nAPOB6d78x2ved6DXPJRTC38l5r4MJCbNrhOHuBfwdsAfweNb+3B7H3v4D7n1REs60hM3w0F6Duy8zs0XAFcBngPcDryAiUkJmNpXw17qnCHlzTZ78OIswIjsBwN17CYMFRM/5irv/IXo84LlZ8ubkiDMwV74InAs8CzwM7DbEl5D9vD2B32Y9/l9CG81MXp+LtzHyFVceJfyCkrHezC4GbibEnqmuc7+uQrYxeN12J3CVmdUT/i02EUbPJWHU4y1p1t93CODujwPnAW8Y5PwJhCS3FHgJWJh90MzeHO0/zswacvZviB4uGGGsfwAuA34cJcaimNkcsibqmNlcwoouexEmW944yFNFREbjw4S/sN1MyI+HRyPNQH9+XA1sBY7N2r+/mU0hJ9ea2fhoICPXYDk5n+8Dyz0sG7uliK9lBWGEvP8tgFcJrTQl4+590eDIvtEvLncAX3X39lK+T8573gbcRvjL7KnAEe7+arneT0ZOhbek3aeiojRjLqE3L2O6mU2I7h9J6O9bD/wE+Cczu8LMjoxGIzYSEuQk4HYze3c0keUNhGL9MeBCM5tgZo2E3uqZZlZjZp83sw8MEqNF2wWEH07fynMs9/zcx+8h9CNm7Amc4O4bgE/leY6ISDFeN5JqZm8HvgKc6e6dhNHll4E7zOy90QTKt0Uj3D8lrH5ygpkdDXzY3bsIufZfzOwrZnYEYQAi06aXnf8Gy8kZ2aPPbwKazGwicDjQYGa7R8d6CJMRMwMw2e/xfeAkM5sWPX4z8AN330aoh/LmYjM7LIqrbpDv3Thy6qnol5PjCH/lnEZoYdwb2AWYmDVxM3dUfWue+AcGlacfJfq32AH4DWFUv5hfSKSC1GoiabcL8JSZ/ZLwZzgDPp51fBvwBTNbTxgd/hSAu99jZv9C+HPlCYTJiU8CmNk/At8j/CC5FjjH3d3M/pmw1NNTwE2EPxm+kXDBiX8mtLAMuHCPmR1KmDQ0j/Cnx18AX4niuRs4DGg2sz0IIzGzgX82syVsH/1ZBExy97acr/06Cxe+2EaYzCQiUrQoj5xGKA5/AHQSJujtCPy9uz8E4O5bo/z4A0Kuu5XQ2w3wr4RBix8R/sKXWUnkvwi57yxCS9xZ7v5KVFgeBuxhZjcMlpPNbCGhr/mjZrbE3TsIEz2/QSi6fwT8I7A/8BxwPWF0/kMWVrvaF/gHM7vD3W8xsz2Bn5vZg4S/mn4hGpk/GZgdfX1PES4mtJOZ7RvF+T5C3/ovsr5vtYTcfyAwNyrMndBbfWL0mn8G7gX+RGifuZOQ028nTOg80MxOcPfM695EuJjRa1GhPhtYaGbLCCPZRM/PvWbFCsLk0ROAemCcmf2aMMGyUkvsyjCY/j1krLJwmd0funtzhd7vHHfX5dxFRMaYqO2lzt0TeQEzMzuHsMLM8uhxPeGvFRdn5gZJMqjVRMayfG0c5Xmj8CfZRyvxXiIiUjnR5NGDElx01wPnA81ZuxuBjSq6k0etJjImRf11HyD86fBkd7+5zG/5v9GfQEVEZGwZ7+5XxR3EYNx9s4XL1n83mni/ktAWc/7Qz5Q4qNVERERERKQC1GoiIiIiIlIBKrxFRERERCpAhbeIiIiISAVUxeTKxsZGb25uLvp5nZ2dNDU1lT6gMlLMlZPGuNMYM6Qz7lLFvHjx4tXunq4vfpSqKWdDeuMGxR6HtMYN6Y29mLgL5eyqKLybm5tpby/+Sq0tLS0jel6cFHPlpDHuNMYM6Yy7VDGb2fMlCCdVqilnQ3rjBsUeh7TGDemNvZi4C+VstZqIiIiIiFSACm8RERERkQpQ4T2E1tbWuEMommKunDTGncaYIZ1xpzHmtEvr9zytcYNij0Na44b0xl7KuKviAjotLS2exp4iEREzW+zuLXHHUUnK2SKSVoVytka8RUREREQqQIW3iIiIiEgFqPDO58474aCD4Lzz4o5EREQK6PpzF+0HtbPk9CVxhyIiMqSqWMe7aOvWwWOPwd57xx2JiIgU0Lexj9cee41x9RpLEpFkU5bKZ1z0benrizcOEREpLPOTbOyvFSAiKafCO5+amnCrwltEJPFsnAHgfaq8RSTZVHjnoxFvEZH0sOhWKVtEEk6Fdz4qvEVE0iNK2dVwXQoRSTcV3vlkCu/e3njjEBGRgjKtJhrxFpGkq3jhbWbvNLNHzKw5a9+JZrbMzNaY2WVmNj7rWKOZXW1mF5vZt8ysdjjHRkU93iIiQEpydlR3q8dbRJKuooW3mc0CJgOHZu2bB5wALAQ+A3wU+HTW024GrnX3c4FngK8O89jIqdVERCQ1Obt/xFt1t4gkXEULb3d/Gbg1Z/euwBnuvtjdfwR8D3gHgJkdBuzt7o9E594GfMrMJg11bNSBqvAWEUlPztbkShFJiYq3mrh7X87j+919W9auDmBFdP8oYFnWuR3R3ZYCx0ZHPd4iIkBacnb/a476pUREyimJkysPAb4f3Z8LrMk53gXMKXBsgM7OTlpaWvq3tra2oSNQj7eIxKitra0/XwGNccdTQOw5W5MrRSROxeTsRF0y3szmA53u/ni0y4FNOafVAT0Fjg3Q1NREe3v78ANRq4mIxKi1tZXW1lYAzGx1zOEMKjE5W5MrRSRGxeTsxIx4m1kNcCbw+azdHcC0rHMMmASsLHBsdFR4i4gMKUk5W5MrRSQtElN4A2cD33D3LVn77gH2znrcDGwF2gscGx0V3iIihSQoZ0e3StkiknBxrOOdmX9uWfvOAx4DGsxsdzM7w8zmu3s78JKZ7RedeixwubtvHerYqIPU5EoRESAlOVutJiKSEhXt8TazKcAHooeLzOwKwp8qL8w5dYm7XxvdPwW4wMyWE9aT/ULWeUMdGzlNrhQRSU3OVquJiKRFRQtvd+8Croq2jK9F22DPeR74SLHHRkWtJiIi6cnZGvEWkZRIUo93cqjwFhFJDS0nKCJpocI7H/V4i4ikR+YnmQa8RSThVHjnox5vEZH0UKuJiKSECu981GoiIpIamlwpImmhwjsfFd4iIukRpWyNeItI0qnwzkc93iIiqdG/1LjGSkQk4VR456MRbxGR9NDkShFJCRXe+WhypYhIamR6vNVqIiJJp8I7H414i4ikR+Zi9krZIpJwKrzzUY+3iEh6ZCZXuka8RSTZVHjnU1sbbnt64o1DREQK0uRKEUkLFd75TJgQbrdujTcOEREpTJMrRSQlVHjnkym8NeItIpJ4mlwpImmhwjsfjXiLiKSHJleKSEqo8M4n0+OtwltEJPk0uVJEUkKFdz41NWAWlhPUyiYiIomWaTVB6VpEEk6Fdz5m6vMWEUkJq4l6vHs14i0iyabCezDq8xYRUl1T8wAAIABJREFUSYesVU00wVJEkmx8pd/QzN4JfA04xd2XR/sagEuALmAacK67bxjNsVFTn7eISCpytplhtYb3ON7jWJ0VfpKISAwqOuJtZrOAycChOYeuBB5w9/OAW4GrS3BsdDTiLSJVLk0522pDsd3Xo6VNRCS5rNKzwM1sHGEKzG7uvtzM5gDLgGnu3m1m44FXgX2BrSM5lhmVyWhpafH29vbiAt11V1ixApYtg+bmkX/BIiKjYGaL3b0lxvdPRc6+f8f76d3Qy+HrDqd2au1ovmQRkRErlLMr3mri7n39l/cNjgRWu3t3dHybmS0DjgB6Rnhs+agD1eRKEZHU5OzMiLf3qMdbRJKr4oV3HnOBNTn7uoA5wLYRHhugs7OTlpbtv3y0trbS2to6dFTq8RaRmLS1tdHW1pZ52BhnLHkkMmePqw2dkyq8RaTSisnZSSi8HdiUs6+OMDrCKI71a2pqouhWE/V4i0hMsgtNM1sdczi5EpmzbXw04r1NhbeIVFYxOTsJhXcHYXZ7tinASkKCH8mx0csU3lu2lOTlRETGiETmbLWaiEgaJGEd798Bc6JlpjCzWmAecO8ojo1eQ0O43ZQ7QCMiUtUSmbO1qomIpEHFC2/bPkvHANx9FXA7cHS0/xjg5+7eMdJjJQl00qRwu3FjSV5ORCSN0pKz+0e81WoiIglW0VYTM5sCfCB6uMjMrnD3TuAs4BIzOxCYDWTPohnpsdFR4S0iVS5NObu/x1utJiKSYBUtvN29C7jq/7d333GS1HX+x1+fyXlnc45sIoMsLOApK0EUOUQOjBjwZMQshp8JMZx3JvQ4A8qq6HkqoqKIKNyhSEZkYdkFZBNhc5zd2cn5+/vjW73T09NheqZ7unr6/Xw86tHdVdVVn66p+cxnvv2tbwVT9PxDwJUJ3jOiZaOmwltEClw+5WyNaiIi+SAMfbzDSYW3iEje0MWVIpIPVHgnosJbRCRvaDhBEckHKrwTUeEtIpI3NKqJiOQDFd6JqPAWEckb6moiIvlAhXciKrxFRPJGUVlwcWW3Cm8RCS8V3olECu/29tzGISIiKRVV+j9nfR19OY5ERCQxFd6J1Nb6x+bm3MYhIiIpRQrv/nb18RaR8FLhncikSf7x0KHcxiEiIikVVxUD0N+hwltEwkuFdyKRwvvgwdzGISIiKamriYjkAxXeiajwFhHJG+pqIiL5ICOFt5ldkInthEp04e10lbyIjB/jMWerq4mI5IOSdFY2s1cB1wAzGCjai4DpQGVmQ8uxykqoqIDOTujogKqqXEckIpKWQsrZR7qatKuriYiEV1qFN/A/wJeBtUCkWcGAN2cyqNCYNAl27YLGRhXeIpKPCiZnq8VbRPJBuoX38865b8bONLMNGYonXKIL77lzcx2NiEi6CiZnH+njrcJbREIs3T7e7zGzS+PMvygTwYTO9On+cc+e3MYhIjIyBZOz1dVERPJBuoX3S4Efmllf1NQPfC8LseXe7Nn+cefO3MYhIjIyBZOzj3Q10agmIhJi6XY1+SLQAKwBIs0KBrw1k0GFxpw5/lGFt4jkp4LJ2cV1vvDube7NcSQiIomlW3jfCfzeOdcVPdPMfpC5kEIk0uK9Y0du4xARGZmCydkl9f7PWe9hFd4iEl7pFt6/BD5lZvdEzTPgEuBDownEzI4G3g9sAZYAq51zT5pZJfA1oAWYCHzCOdccvCfhsoxQVxMRyW8Fk7OPFN5NKrxFJLzSLbw/i0+w7wAid5Ux/Bixo0ri+GGvXuuc22lm84D/BY4GbgDucs7dEoxJeyPwpuA9yZaNXqSriVq8RSQ/FUzOVuEtIvkg3Ysrvw7Mdc4tcM4tDKYF+NaT0ToaqA2edwATzGwWfrzZPwTz/wxcZGYLki3LQCzeokX+ccsW6NcFOyKSdwomZxfXFEMR9Lf109+jfC0i4ZRu4b3dOdcSZ/6MDMTyG/zV97XA5cAHgFXAAedcO4Bzrhd4ATgrxbLMmDgRpk6F9nY/nreISH4pmJxtZurnLSKhl27h/d7YGWZ2GXB9BmJ5L9AJPAa0OuduBWYDjTHrtQCzUiwbZP/+/axYseLItHr16uFHtXSpf9y4cfjvEREZhdWrVx/JV8CUUWyqoHK2upuISC6kk7PT7eM9ycw+65z7NzOrAr4FXIZPvKNVCdyMT87Xm9kL+D6JHTHrlQM9wfNky46YOnUqa9asGVlUy5bBQw/5wvucc0a2DRGRNDQ0NNDQ0ACAmR0YxaYKKmdHCu++w7qJjoiMnXRydlqFt3PuNWZ2upl9Gzgf2AOcCBweYazRfga80TnXZGYOn9A/jL/yPVotsAuf4BMty5zly/3jP/6R0c2KiGRboeXs0smlAHTv787UJkVEMipp4W1mL0+waCuwGbgOmA9cDFw90iDMbApwonOuKZj1JfxV+FuBWWZW6ZzrMLNSYB5wH/5mED9IsCxzTj7ZPz7+eEY3KyKSaYWes8tmlAHQvUeFt4iEU6oW76/hr1xvAuJdJn4TfmiqmYwiiQMHgU4zm+2ciwya3QisB+4AzgkezwVujaxjZgmXZcwpp/jHJ5+Enh4oLc3o5kVEMqigc/aRwnu3Cm8RCadUhfcXgTXOuX3JVjKzS0cThHOu38wuBq41s8eB6cDHnXPNZvYe4GtmdiL+SvyGqLcmW5YZEyfC4sV+SMFnnoGTTsr4LkREMqSgc7ZavEUk7JIW3s65Pw1nI86534w2EOfcOuDdceYfAq5M8J6EyzLq1FN94f3IIyq8RSS0Cj1nl81U4S0i4ZbucIKFadUq//jnP+c0DBERSUxdTUQk7FR4D8d55/nHe+6BXo0PKyISRpEW766dXTmOREQkPhXew7Fwoe/n3dQEj2Vi+FsREcm0ivkVAHRt66K/V7eNF5HwGXXhbWYnm9nPzOyjZrZg9CGF1AUX+MffjLprpIhIzoznnF1cWUzZ7DJcr6Nrm1q9RSR8MtHifRVwrnPuG8CZZlaZgW2Gzxvf6B9vuQX61ZIiInlrXOfsyqP8x+l4LvYmmSIiuZeJwvv9QDDYNb8ExudA16efDvPnw86dcP/9uY5GRGSkxnXOVuEtImGWduFtZq82s7cEz08DlkRugOCc63fONWc4xnAwg7e+1T//3vdyG4uIyDAVWs6uXBwU3ptUeItI+KRVeJvZvwO/Bl4H4Jz7O/A2M3tzFmILn6uugpISuPVW2L4919GIiCRViDm7+rhqANqeastxJCIiQ6Xb4n0mMAd4IGrencBXMhZRmM2eDZdeCn198J//metoRERSKbicXXNSDQCtT7binMtxNCIig6VbeN/tnGsCorPZqG49nHc+8Qn/eMMNsG1bbmMREUmu4HJ2+dxySiaV0HOgR+N5i0jopFt4P29mnwaWmdkbzexW4H3AVzMfWkiddBK86U3Q1QXXXpvraEREkim4nG1mA63ea1tzHI2IyGBpFd7OuV8C9wDFwNuAJuCVzrnvZiG28PriF6GsDP77v+Hee3MdjYhIXIWas+tOqwPg8IOHcxyJiMhgaY9q4pz7m3PuKufcBc65fwX+ambnZCG28Fq8GD7zGf/8yiuhvT238YiIJFCIObt+VT0ATfc15TgSEZHBSpItNLPbgZoU25gCHAD+kqmg8sInPwm//jU8/TS8731w001+yEERkRxRzvbqzqyDYmhZ00JvSy8ltUn/1ImIjJlU2Wh9sM5GoB+4ElgLrIlaZyEwKSvRhVlZGfziF7ByJfzkJ3Dmmb71W0Qkd5SzgZLaEmpX1NLyaAtNf21iykVTch2SiAiQuvD+OmDBVfGY2VHOuSFXFJrZ77MRXOgdfzx8//vw9rfDe98Lc+bAq1+d66hEpHApZwcmXziZlkdbOHDbARXeIhIaSft4O+cORxJ4YLqZFUevY2argNOyEFt+eNvb4OMfh95eP8b3ww/nOiIRKVDK2QOmvm4qAAduP0B/b3+OoxER8dLt+PYbYJ2ZPQi0AEcD5wOFPa7eV78K+/b5UU5e+Ur4/e/hnHF97ZKI5IeCzdlVx1RRubSSjk0dNP2liUnnj+veNSKSJ9IdTvBu4JX4C3OW4YemutQ59+VMBWRmp5vZx8zsYjObYmaVZvZtM/sPM/uemdVFrZtw2Zgygx/+EC6/HNra4IIL/IWXIiI5VMg528yYfvl0AHb/cHe2diMikpaRDCe4yzl3jXPuIufc5cAfzOyUTARjZlcAFzrnrnPO3eacOwDcADzonPs08Hvgxqi3JFs2tkpKfIv3Bz4A3d3w+tfDpz/tby8vIpIjhZyzZ75zJhTDgdsO0L23O5u7EhEZFnPOpV4rsrLZTXFmTwX6nHMXjyoQs7OAzwLnuSAoM5sFvABMdM61m1kJcBg4FuhOtMw592L0tlesWOHWrIm+qD+LnIPrr/f9vvv64Lzz/Kgns2aNzf5FZFwxs8edcytG+N6Cz9lPXfwUjb9vZP5n57Pwiwszsk0RkURS5ex0W7yXANuArVFTOfD4iCMc8E3gWeDbZnaXmb0UWAUccM61AzjnevGJ+6wUy3LHDK6+Gv7v/2DyZLj7bj/6ya9+ldOwRKQgFXzOnvuxuQDs+NYOepp6srkrEZGU0r248l3OuY3RM8xsOfCO0QRhZsuAk4ErnHPrzezjwF3Al4DGmNVbgFlAb5Jlg+zfv58VKwb++WhoaKChoWE0Iad29tmwbh3867/C//4vvOENftzv66+HBQuyu28RyWurV69m9erVkZejGQuv4HN2/T/VU/+Kepr+2sSO63ew8PNq9RaRzEonZ6fV1STuBswWAo8550b8x8HMLgFWR7ZhZpXAfuCLwL8451ZGrfsE8Ivg5WXxljnnrove/ph2NYnlHNx4o+960toKFRXwqU/515WVuYlJRPLGaLqaJNheweXspvubePKsJymqLOK0Z0+jYn5FxrYtIhIto11NzOwFM3s+atqK/6rxd6OMs4So1nfnXAewGSgGJsasWwvsAnYmWRYeZnDVVbBhA7z5zdDZCZ/7HBx1FNxwg78QU0QkC5SzvfqX1zP1DVPp7+hny9Vbsr07EZGE0u3jfTfwTuCKYLocOME5N9p7pa8HJpjZ9Kh5vfiEPCtoTcHMSoF5wH3AX5MsC5/Zs+HnP4d774WXvAR274b3vQ+WLoUf/MAX5CIimaWcHTjquqMoqi7iwO8OsPcXe8dilyIiQ6RbeP/MOXevc+6+YHoAMDNbPJognHMbgDuBS/EbrAfKgJ8DdwCRu9GcC9zqnNvpnNuTaNloYsm6s86CNWvg1lvhmGNg61ZoaPD9vv/936ExtgukiMiIKWcHKuZUsPg//cfe9J5NdLzYMRa7FREZJN3C+1Vx5u0FMnG3mLcBZ5rZ1cDngDcFV72/B3itmX0GuACIvsIm2bLwMoNLLoH1630r+Eknwd69cM01MG+e75qydm2uoxSR/KecHWXmu2Yy5eIp9DX38cwlz9DXpvssiMjYGtbFlWb2fvzV6rUJVrnbORcvwYdCTi+uHA7n4J574Lrr4K67BuavWAHvfje88Y1QU5O7+EQkZ0ZycaVydmI9jT08vvJxOp/rZMrFUzj21mOxIsvKvkSk8GTk4krn3HeAlwPXA4tipjlhTuB5wQzOOQfuvBOeeQY+9CGor/ddUq68EmbOhLe/3Q9L2Nub62hFJOSUsxMrnVzKCX88gZL6Eg7cdoBN796E6x/d6F4iIsOV7p0rS51zQ+5AYGbznHPbMhpZBoW+xTuejg74zW/8UIQPPTQwf9o0Pyb4m98MK1f6ol1Exq1R3rlSOTuBQ/ce4qlXP0V/Zz8zG2ay9HtL1fItIqOWKmcnLbzN7GPAn5xz/whefwY/XNSg1YBXOOdWjT7c7MjLwjva5s1w882+P/imTQPz58yB174WXvc6ePnLobQ0dzGKSFakU3grZ6fn4J8P8vQ/P01/Zz/T3jyN5Tctp6g83UufREQGjLaryUKgPur1S4Az8V9XLoyaZowyTklmyRK49lo/FviaNf6W9LNnw44d8N3vwrnn+pbwt77Vj5bS3JzriEUkN5Sz0zDp3Ekcd/txFFUXse8X+1j3ynX0HNRt5UUke9LtavIS4GnnXHfM/GMiLSxhlPct3vH098Pjj8Ntt/npH1GHv6QEzjgDXvUqOP98OPlkKFIrjkg+GmVXE+XsYWhZ28JTFz5F965uKhZUcMwtx1B3Wt2Y7V9Exo+M3rkSOD5OAl8GHDWS4GQUiorg1FP92N/PPAMbN8JXvwovfakfJeWBB+Azn/Ejo8yYAZdfDv/zP7B9e64jF5Gxo5w9DLUn1/KSR19CzSk1dL7Yydp/Wsv267froksRybjhDif4Mny/wAZgdcziScC3nHPzMh9eZozLFu9kmpr88IR33eVHQtkWcw3VokX+Rj6rVvnH+fNzEqaIpDbC4QSVs0egv6uf5/7fc+z8lr+nT/2qepb9cBmVR1WOeSwikp9S5eySYW6nCfgpvm/gS2OWdQM/HFl4khX19f4GPZdc4lu/N270Rfif/+xbwp9/3k8//rFff8ECX4CfdZbvorJ0qbqmiOQ35ewRKCovYsl/LaH+FfVsathE071NPHb8Yyz8t4XM/uBsikqVF0VkdIbdx9vMJgCXOefyLmEXXIt3Mn19sG4d3Hsv3Hcf3H+/byGPNnGiH6rwjDPg9NP98wkTchKuSKEbaR9v5ezR6T7QzXNXP8fen+0FoHJZJYu/uZhJr56EaRhXEUlgVMMJprGTf3LOPTjqDWVJGJJ4aPX1wVNP+UL8wQfhkUdg167B65jB0Uf7QnzlSjjlFDjuOCgry0nIIoVkNBdXJtmmcvYwNf6pkS0f3kLH5g4AJp4/kUX/sYjalyS6KaiIFLLRjuN9O5DqXuUGHOucmzayELMvTEk8L2zf7gvwv/3NPz7xBHR3D16ntBSOP94X4S95iX88/nioqMhNzCLjVJrjeCtnZ0F/dz87v7uTF7/wIn2H+wCYfNFkFnxugQpwERlktH28N+L7Cu4EElXoBlw4svAklObO9dPrX+9fd3XB2rW+CH/sMV+Ib9rkH594YuB9JSVw7LG+CD/5ZF+IH388TJqUm88hUniUs7OgqKyIuVfPZfpbp7PtK9vYdcMuGm9vpPH2RiZfOJm5H5vLhJdPUBcUEUkpVYv3FKDFOdeVdCNmk51zjZkOLlPC1noyLrS0wJNP+rHEH3/cF+AbNvjxxWPNnj1QhJ9wgn9cvhzKy8c+bpE8k2aLt3L2GOje2822r/sCvL/D57yak2uY8+E5THvDNN39UqSAjVUf73eF+QKesCfxcaOtzV+4+fjj/nH9ej/GeHv70HVLSmDZsoGC/JhjfD/yRYt8NxYRAbLWx1s5OwO693az84ad7PreLnr2+ztels0oY8Y7ZzDznTM1DKFIAcpq4W1m84D3Ae93zlWPeENZli9JfFzq7/dDF65f7y/ifOop/3zLFj/UYazSUli82Bfh0dOyZVAd2lNMJGsyWXgrZ2dHX2cf+36xjx3X76DtqbYj8+tX1TPjX2cw9V+mUlxZnMMIRWSsZLzwNrNi4CLgSuA8YG+wndmjCTSb8i2JF4T2dt8a/tRT8PTT8Oyzftq6NfF75s0bXIgvWeKnOXM07riMW6MtvJWzx45zjsMPHGb3j3az/9f7j3RDKa4rZuolU5n2xmnUn12v8cBFxrGMFd5mthB4F3AFMBnYCnwQ+F/gDOfcw6MPNzvyNYkXpLY2f8OfSCG+YYN/3LwZenriv6eiAo46aqAQX7LEt5ovWQKzZqkol7w2inG8lbNzqPdwL/t+uY/dP9pNy2MtR+aXTC5h6qVBEf6yeqxYF2SKjCejLrzN7DJ8S8nZwBb8Hc/+G/i0c+7qDMaaNeMhiRe8nh7fZSVSkG/a5IvxzZth377E76usHCjCIwX5UUfBwoW+pbxkuDdvFcmNdAtv5ezwadvQxv5b9rPvl/to3zBwzUvplFImXziZyf88mYmvnEhJjfKRSL7LxC3jO/C3GH4YeJdzblOw4dFflRnDzMqAx4APOefuNbNK4GtACzAR+IRzrjlYN+EyGYdKS333kmXL4OKLBy87fNj3Gd+yZaAY37zZv96/f6BveaySEt99ZeHCgWnRooHnU6f6mweJ5Jcxy9nBdpW3U6heXk3156qZf+182p5qY98t+9h3yz46n+tkz0/2sOcne7AyY+LZE5l80WQmv2YyFfN0TwSR8Shl4e2cuwO4I7go50ozWwD8Hjgy9ISZmcvE8CjwcWBB1OsbgLucc7eY2auAG4E3DWOZFJIJE/zY4aecMnRZU9NAQR55fP55eOEFf4fO55/3UzzV1YOL8khhvmABzJ/v9ysSMmOcs0F5e9jMjJoTaqg5oYaFX1pI+7PtHLj9AI23N9L8t2YO3nWQg3cdZDObqVxaycTzJjLpvEnUv6Kekjq1houMByO5uLIEuBh4C/A08CPglc651aMKxOxMYDlwLfAOYBPwAjDROdce7PcwcCy+NSfuMufci7HbHm9fW0qGdHb6izkjhXhkirxuakr+/ro632I+f75/jH0+c6a6ssioZeDiyqzk7GDbWcnbhZizu/d10/hHf1OeQ385RF9L38DCYqhbWcfE8yYy8eyJ1J5WS3GFRkkRCaNMdDUZxDnXC/wG+I2ZLQauxvcnHHESN7Nq4FLn3EfM7Npg9irggHOuPbJfM3sBOAvoSbLsxZHGIQWmomKg+0o8TU2DC/Hoads2aG72I7I8/XT89xcX+37k8YryyOuaVHf3FhmdbORsUN7OtLJpZcy8YiYzr5hJf08/LX9v4eDdBzl09yGaH22m+WE/bf3CVqzMqDutjgkvm8CEl09gwpkT1CIukidG9ZvqnNsCfMjMHh9lHJ8AvhIzbzYQe2e1FmAW0Jtk2RD79+9nxYqBfz4aGhpoaGgYTbxSCOrr4eST/RTLOTh40BfgW7f6x9jne/b418mGSJw40Rfn0dPs2YNf19Wpr3mBWb16NatXH6mLp2RquxnM2ZDFvF3oObuotIgJL53AhJdOYOHnF9J7uJeme5s4ePdBDt9/mLan2zj84GEOP3gYvgwUQc1JNb4QP3MCdSvrKJ9XrlvYi4yRdHJ2Ru5cORpmdj5Q7py7PXj9Iv4ryxXAZc65lVHrPgH8IngZd5lz7rrYfRTi15YSAl1dsH17/KI88rwr6Z29vZqaocV4bIE+ZYqK83EqG3euHK1s523l7OR6DvVw+KHDHL7/MIcfOEzLmhZc7+C/5aXTS6k7rY66lXXUrqyl7tQ6SiaoVVwk2zLe1SQLPgacHPWf+UT8hUDfDJ5HqwV2AS7JMpFwKC/3wxcuXhx/uXN+1JUdO2DnTv8YPe3c6Qv31lY/tvnGjcn3NXv20AJ99mw/lvmsWb7PeVlZdj6rFBrl7RwqnVjKlAunMOVC37DW19ZH86PNHH7gsO+W8mgzPXt7aPxDI41/GPiSoWp5FbUra6k9pZaak2uoObGGktowlAEihSMMLd5TgfKoWY8AHwEeBf4BTHXOdZhZKdAKLAL68OPTDlnmnNsZuw+1nkjecs4PlxhdjMcr0A8dGt72pkwZKMQTTdOn66LQEAlpi3dW87Zy9ug45+h4roOWR1to/rsvxFvXtuK6h/69r1xc6YvwYKo9uZay6foHXWSkQt/i7ZzbH/3azPqA/c65bWZ2B3AOcAdwLnBrJEEnWyYybpj5vub19XDccYnXa2sbKMpji/Pdu/3QiXv2wIEDflq/Pvk+p09PXaBPnaq7ghYo5e1wMzOqFldRtbiK6W+ZDkB/Vz+t61pp/rsvwlvXttL2dBsdWzro2NLB/l8P/EjLZpQdaRGvPq6aqmOrqFpepZFURDIg54V3Cu8BvmZmJwIzgIZhLhMpLNXVsHSpnxLp6/NdW3btGjxFCvPItHevL9L37IEnnki8veJimDEjdYE+ebL6nxcW5e0QKiov8n2+T6s7Mq+/u5+2f7TR+mTrkWK89clWuvd0c/DOgxy882DUBqBySSXVx1ZTfVwwHVtN5ZJKikr1D7jIcOW8q8lY0NeWImno7fXFd2yBHjsdODC87ZWV+f7lkUI80u889rG2NrufK0+FsatJtiln547rd3Q833GkRbztmTbfMr65A/qHrm+lRtXyKqqPrabq6CqqllVRuaySqqVVFFephVwKT+i7mohIyJSUDFyomUxXl28VT1WgNzWlHlYRfOGdrDCfPdu3sOsCUZGssaKBbipcNjC/r7OPjo0dvhiPKsg7X+ik7ak22p5qG7Kt8rnlA4X4Ml+UVy2vonxOOVakb8GkMKnwFpGRKS/3NwGaPz/5eh0dvjvLzp2+EI/3uHMntLSkHr0FYNq01AX65Mnqfy6SQcUVxdSc6Pt9R+tt7aX92Xbanm6jfWM7HRs7/OOWDrq2d9G1vYtDfx588XdRZRGVS4NifGkVFUdVUHlUJZVHVVI2s0zjj8u4psJbRLKrshIWLfJTIs75lvFEhXnkcc8e2LfPT08+mXh7paWD+5nPng1z5w6+g+j06SrORUappKaEulPrqDu1btD8/t5+Ol/oHFSMt29sp31DOz37emhb10bbuqGt5EVVRVQuqhwoxhdXHinKy+eXU1Si31nJbyq8RST3zPxdPCdOhGOPTbxeX99A//NkBfqhQ6m7t5SVDS3Go5/Pnev/aRCRtBWVFFG1pIqqJVVw4eBlPU09Ay3jmzvoeG5g6m3sPdKdZYhiqJhfMagYr1hYQcUCP5VMLFFruYSeCm8RyR/FxQOt2CuSXG/Y0THQxzzSlSVyF9HIXUMPHIDnnvNTItOmwVFHwZIl/kZIS5YMTHV1id8nIgmV1pdSurKUupVDf4d6mnrofK5zoBjf4h87n+uka0cXnc930vl8J4cYeu+C4triI0V4xfyBgvxIYT5JhbnkngpvERl/Kit9wXzUUYnXaWsbWoxHHrdt88si3VoeeWTo+6dN88X40Ud2NbFXAAAer0lEQVTDiSf66YQT/JjrIjIipfWllJ5SSu0pQ0c56uvso/OFzkHFeOfWTjpf7KTzhU76WvoSXugJUFxTPKQYr1hQQfn8cirmVVA6tVSFuWSdCm8RKUzV1bB8uZ/i6evzLebPPQebNw9MW7b4KVKUP/zw4PfNn+9b41et8tMxx6gvuUgGFFcUU310NdVHVw9Z5pyj91CvL8Ij09ao55HCPFE3FqCooojyOeWUz/VTxbyKI8/L55ZTMbeCkgkqm2R0dAaJiMRTXOz7ec+d6wvoaP39vvvK5s3wzDOwbp2fnn56oG/5rbf6defPhyuugKuvVvcUkSwxM0onlVI6qZTalwxtLXfO0dvUm7Ao79reRe+h3iN38kykuK74SBEet0CfU05xpcYvl8R0Ax0RkUzp7YVNm3wr+L33wj33+KEUwd9VdN06qKhIa5O6gY7I2Oht7T0yBGLX9i46t3fStS3q+fYu+tvj3EUoRumUUsrnDbSSl80u8y3pswem4moV5+OVbqAjIjJWSkp815Kjj4YLLoD774ePfMQX35s2wY4dvl+4iIROSU0JJUeXxO3KAlHdWbZ1Di7Qo1/v6KLnQA89B3pofaI18b7qSwYX5NGF+ZxyymaXUTpZfc7HIxXeIiIj1d3tu5VE+n1v2eIL7LVr/bCHEWbwgQ8kv9hTREJtUHeWk4Z2ZwFwfY7uvd2DWsm7dnbRvbObrh3+edfOLnqbeult6qX9mfbE+yu3Qa3kkYI8ulAvm1lGUamuIcknKrxFRKI5B4cPw/79fsjB/fv9RZSxd9rcudMvS6S+Hk45Bc4/H173OrV0ixQAKzbKZ5VTPqs87nCJ4FvOew70xC3II8+7d3b7PunB8ImJdwhl08uGFOSxhXpJrcq9sNBPQkTGp74+aG72N9Npahp4jH5+6BA0NvoCOlJoHzgAPT3D20dxsb8rZmSc78WLfav2iSfCwoW+pVtEJIqZUTa1jLKpZXBS4vX62voGFeOR5907uwee7+k+MrU+nrhrS3Fd8eDCPLgQNHrSDYjGhgpvEQmXvj5obfVFc0uLn1I9b272rdTRBXVz88hjqKmBqVP9NGWKf5w9e+g0bZovvkVEMqy4upiqpVVULa1KuE5/bz/de7qHFOSxhXpfcx/tze20P5u4a0tRVdGQYjy2QFe/89FT4S0iI+Ocv0NkW9vA1No6+HWyZa2t8Qvp9sR/GNJWV+dvQ19fP/AY+3zy5IEiO1JopznyiIhILhSVFFExp4KKOYlzVuSi0CMF+Y6YKbgwtK+1j45NHXRsSjyc4pGxzpMU6KVTSrEiFeeJqPAWGa96enwR29HhH2Ofx1uWTvHc3u6L70wz8y3OtbV+qquL/zz2dWxRPWGCWqNFpOBFXxRac0JNwvV6m3sHCvGY4rxzeyddO7roO9yXcqxzK7NBRXnF/Jg7hc4rp7iicHOzCm+RseQcdHYOrwhO9TzVer292f88FRX+DpA1Nf4x0RRveaS4ji2mq6t1p0cRkTFWUldCyTElVB8TfzhFgN6WXt9qHqc4j8zrPZT6otCymWWDivFCKsxVeIvA0NbhVIXvaJaPleJiX8RWVkJVlZ8SPY+8Tqd4rqpSi7KISAEpqS2hZHkJ1csTF+dHLgqNugnRkbuEvthJ57ZOund30727m+ZH4l+LUzazjIqFFVQuqaRqadXA4+LKvL/5UGgKbzO7BPgGUAf8DPioc67XzKYA/w40AaXAJ5xzPcF7Ei6TcaC/3xeqwy12R1MYj0XrcERFxfCK4HSex1tWWjp2n0kKjnK2iMST6qLQ/t5+und1Dy7Go6au7V0DhfnDQwvzstllVC2ponJppS/Il1dRc3wN5fPK8+LCz1AU3mY2D7gYuBQ4FvgesB24Dvg18Enn3KNmdhXwReBTwVuTLZNscM7fNCTTrcDx5nUmGbs004qLUxe+qQrgVPOqqnzRrW4UkueUs0VkpIpKiqiYV0HFvAp4+dDlrs/RtauLjuc66NjsL/Zs39Tun2/poHtnN907u2m6t2nQ+4rriqk+vpqa42uoPqHaPz+xJnRjmIclmvnAO51zvcDjZnYc8AozexBY5px7NFjvD8BGM/sScHyiZc65trH+AKHU0zNwQVy8qaUl8bJ4U6Qo7u8fu88w2sJ3uMvVOiySDuVsEckKKzYq5lZQMbeCiasmDlrm+hyd2zp9Mb65nY5NHbT9o4229W307O+h+aFmmh+KaiU3qD6umroz6qg7vY66M+qoWlqV01FXQlF4O+ceiJm1E6gFzgZeiFpvZ/A1wgrgpUmW3ZflkLOrv39gXOKmpuE9NjcPLZS7urITX1nZ2BTDFRW6AYlICClni0guWLFRubCSyoWVTDp/0qBl3Xu7aX2qlbb1bbQ91Ubr+lbanmo7Mu1evRuA0umlTDp/EpNeNYlJr5xE6eSxbXgLReEdx6nA14EGoDFmWQswC5idZNkg+/fvZ8WKFUdeNzQ00NDQkMl4E+vuHrgbXvQtqGMfGxsH7qrX0pKZYdqKivwoETU1qadU60UuqKushJKwnjYi48Pq1atZvXp15OWUXMYyTOMnZ4tIXiqbXsak6ZOYdO5AQd7X0Ufr2laaH2mm+W/NHH74MN27utn7073s/eleKIbJr57MjHfMYPI/T6aobGRdQdPJ2aGroMxsMbDfObfOzBwQOwxEOdADJFs2yNSpU1mzZk3mg3UOdu2CzZth+/aBaceOgecHD45s23V1fhziyHjEiR4jzyNDskUXy+XlajEWyUPRhaaZHchxOEnlVc4WkYJSXFnMhDMnMOHMCYC/mVD7P9o5eNdBGu9spOneJhrvaKTxjkYqFlQw/9r5zHj7jLS7oqSTs0NVeJtZMfBu4JPBrJ3A0qjlBlQDu1Isy7y+Pli/Hh54AJ54Ap59FjZsSH1b6uLiwXfGi9x+Ovr5lCl+mjjRF9K1tRqmTURCL9Q5W0QkhplRfWw11cdWM/ejc+ne283eX+xl14276NjYwcZ3bmTvz/dyzM3HUDa1LCsxhKrwBj4GXOeci3ROvhu4Kmr5AqAbWBM8JlqWOffdBzfdBLfdFr/InjQJli+H+fNh7lyYM2fw49SpGsVCRMar8OVsEZFhKptextyr5zLng3PYe/NenvvIczT9pYn156/n5AdPprgq842goSm8zezT+ARcaWaLgFXA/cBuMzvOOfc08Crg2865bmCNmSVaNnqHD8Pb3ga33z4wb8ECeNnL4Iwz4Jhj4OijfWGt7hwiUmBCl7NFREbIio0Zl89g4tkTefKsJ2ld28r267az4NoFGd9XKApvM7sG+LeY2RucczeZ2RuAa83sRaAG+EzUOsmWjc5b3gJ//KPvO/3BD/oifMmSjG1eRCRfhTJni4iMUvmscpbeuJR156xj9w93M/+z8zN+Ux5zmRg9I+RWrFjh0rpQZ8MG35pdWwvr1sHChdkLTkQkCTN73Dm3IvWa40faOVtEJENcv+OhyQ/R29TLGbvPoHxGeVrvT5Wz1fk4nrVr/eO556roFhERESkQVmRULKoAoGtr5u+HosI7nkOH/OP06bmNQ0RERETGVMkE3xO7r60v49tW4R1P5I6P5el9vSAiIiIi+a2oypfHfe0qvMdGpN+7RisRERERKShFFb487u/oz/y2M77F8UCFt4iIiEhBsmJf/7m+zA9AosI7HhXeIiIiIgUpUniT+Z4mKryTUuEtIiIiUlDU4j3WCmBscxERERGJI7hTvArvsaKuJiIiIiIFSS3eY02Ft4iIiEhBUh/vXFHhLSIiIlJQ1OItIiIiIjIW1MdbRERERCT71OItIiIiIjIG1MdbRERERGQsBNWx61eLt4iIiIhI1lhkcI0s3NZFhbeIiIiISEQWB7VT4S0iIiIiEisLLd4lmd/k2DKzSuBrQAswEfiEc655VButqYGZM6G2NgMRiohIRFZytohIBpVMKKFsZhlF1Zlvn877whu4AbjLOXeLmb0KuBF406i2+L73+UlERDIt8zlbRCSD5n9mPvM/Mz8r287rriZmNgt4M/CHYNafgYvMbEGuYhIRkfiUs0Wk0OV14Q2sAg4459oBnHO9wAvAWbkMSkRE4lqFcraIFLB872oyG2iMmdcCzIqesX//flasWHHkdUNDAw0NDdmPTkRkBFavXs3q1asjL6fkMpYMU84WkXEnnZyd74W3Azpi5pUDPdEzpk6dypo1a8YsKBGR0YguNM3sQI7DySTlbBEZd9LJ2fne1WQn/qr4aLXArkxsPOq/l7yhmMdOPsadjzFDfsadjzGPAeXsOPI1blDsuZCvcUP+xp7JuPO98P4rMCsYngozKwXmAfdlYuP5eIIo5rGTj3HnY8yQn3HnY8xjQDk7jnyNGxR7LuRr3JC/savwDjjn9gB3AOcEs84FbnXO7cxdVCIiEo9ytogUOnMuC7flGUNmNhF/M4YXgRnAp5xzrTHr7Ae2jmDzU4B861+pmMdOPsadjzFDfsadqZjnO+emZmA7oaCcHVe+xg2KPRfyNW7I39jTiTtpzs77wltEREREJB/kdVcTEREREZF8ocJbRERERGQMqPAWERERERkDKrzjMLNKM/u2mf2HmX3PzOpyHdNwmNmXzMwF07pcx5OImZ1nZo+a2YKoeaE/5vHiDuaH8rib2SVm9oKZNZrZf5lZSTB/ipndaGZfNbNvBkO6hUKimINlP4s6zr/PZZyxzOwsM9toZofM7Pqo+aE/r8NuuMfQzF5tZt8xsxvM7NLhLsum4cRuZhPM7Ndm1mxma83s9JjlS8ysL+rcPzkMcQfr1QdxR2J7XdSyMB/zV0fFHJn+PpzPNQbxx/07E7NOGM/1pHGH8TyP2vdwjnnmznXnnKaYCfgx8Ibg+auAm3Md0zBirga+jx+e61xgSa5jShDndOB1+DvYLciXY54k7lAed/zYyD8FTgHeBrQBHwuW/RVYGTy/CvhyruMdRszzgG9GHee5uY43Ku4JwKfxN4Z5Df4ujOcGy0J9XufDNJxjCCwFngCKgtf3A8elWhaS2P8DeFNw3v8R2AtURy3/KvDK4Lw/KyxxB8s+HuTFyO9lcZ4c868A/wKcHkzXANek+lxjEHvcvzMx64TuXB9m3KE7z4cbe7JzYiTHfEw+WD5NwCygC6gKXpfgC4CEP5AwTMCHg6ki17EMI9ai6JM8X455bNxhPu7Ay4CSqNdfC5Ld6cCuqPmzgdboBBi2mIPn/wlcGr08LBNQSTBCVPD6ceAV+XJeh3ka7jHE//P7H1GvPwP8ONWykMR+dtTzCfh/3E4NXk/Hj3s+J4THvBy4mziNDXlwzOfFvP4xcEyqzzVGx3/I35nhHttcHfdhxh2q8zzN2DN6rquryVCrgAPOuXYA51wv8AJwVi6DSsbMDHgj8A1gr5m9KcchJeWc64+ZtYo8OOaxcYf5uDvnHgiOY8ROYBtwNv7YRtaL3LhkxRiGF1eimM2sCvhn4FfB63PibiBHnHMdLsi4ZlYNPAvcS56c1yG3iuEdw0HnNbAxeG+qZdm0imHE7py7J+r5YaAZ2BHMehP+n7jtZvbj4PwKRdzAhcBKYJOZ/dHMpkQtC/sx3xZ5bmbFwFLn3D+CWck+V9bF+fsYK4znesq4Q3ieR8eW6phn9FxX4T3UbKAxZl4L/j/pUHLe6fgB3r8N/NzMXp3jsNKRd8cc8u64n4r/zzyfjvWpwPedc+3OucX42P8I3Glmx+c2tKHM7GXAnUApUEV+HeuwGu4xjF0vep1ky7Ip7Z+/mS0F7nXO7QZwzl0P1AEX4b/e/lF2Qh1kWHE7527Ft1yeBSwAbjezSE2RN8cceDm+ewCQ8nOFQRjP9bSE5Dwftkyf62E6mcLCAR0x88rxX4uEmnPukHPuGuDL+O4P+SJvjzmE/7ib2WJgv3NuHXlyrGNiBsA5t9s5dyXwa+D9OQsusS3AT/C3Q/86eXKsQ264xzB2veh1ki3LppH8/N+H70s6sBHn+pxzfwDOBy4xs5kZjXKoYccdND7cj2/1W45vFYy3jTAf84uB3w3aSOLPFQZhPNfTFYbzPC2ZPNdVeA+1E3+RVLRaYFcOYhmpb+EvRssX4+GYQwiPe/A16ruBTwazBh3roLtMNSE61nFijhW64wxH/jG4CfgYvmVkvJzXuTTcYxi7XvQ6yZZlU1o/fzM7F98K+Hy85UFXiL+Q/XM/7fPWObcXuIWB2PLimAdWAI/FWxDnc4VBGM/1YQvReT4imTjXVXgP9VdglplVApgfam0ecF9Oo0pPP/4q23wxHo45hPO4fwy4zjnXFby+G1gWtXwB0A2sGeO4komNOVYYj3O0NfhkPF7O61wa7jGMPa+XBu9NtSybhv3zN7OjgUXOud/FLovRBmzIdKAxRnre9gFPBs9Df8yD5acAT0Su0Ugg+nOFQRjP9WEJ2Xk+GqM611V4x3DO7cFfXRu5gOtc4Naoi9BCx8xmmdmbzawoaMH8KP7K2lAKYgQwyJ9jHht32I+7mX0aWAtUmtkiM3sn0ATsNrPjgtVeBXzbOdedqzijxYvZzI4zs4uC5WXA2/EjnoSCmVXFjDl7AfBf+XJeh1myY2hmHzOzyB+87+O/oo78np6N/2Yk1bKcx25mi4AG4P/MbIGZnWhm7wmWvcPMpgXPXwbcF1yYFoa4Lw1ix8yWAwedcxuD94T6mEe5GPht9IwUnyvrYv/OBPNCfa5H7Q8SxB228zzN2DN7ricb8qRQJ/zXBj/AF1HfBmpyHVOKeI/Fj1jxdHASLM91TElirQXeg+8X9QVgaj4c83hxh/m448eldTHTs8Gy+fiLVz6LHzs1FEP0JYoZ/0d0H/Ao8F/ArFzHGhP3acB+4AF83+7LopaF+rzOhynRMcQP23hp1HqX48d6/ybw2phtJFyWy9iBmUEOiT3vI+NQ/x44ANwMvC0scQfPv4P/R/63wAcJxjEO+zGPWfdeYsboTvW5shx7or+PoT7XU8Ud1vN8uMc80+e6BW8SEREREZEsUlcTEREREZExoMJbRERERGQMqPAWERERERkDKrxFRERERMaACm8RERERkTGgwltEREREZAyo8BYRERERGQMqvEUkb5jZVDO73cy2mtkncx2PiIgkp7w9mG6gIyJ5w8zOB/4K1AEbgKXOuYO5jUpERBJR3h5MhbeI5CUzux842znXm+tYREQkNeVtdTURkTxkZnOB/yvk5C0ikk+Utz0V3jKmzOz/mdl3cx1HRCbiMbNzzewRM1uVobDGFTO7LsPbM+D1wFcyuV0RGUz5evTMbIGZ/cDMbszQ9orN7Eoz2zqMdV9jZo+O9LNmMncrbw9Q4T2OmNnRZvZbM/uhma0zM2dml+c6rhh/Am7OdRBRMhHPduDEDMQSSmZ2gZkdMLPdZnaGmc0NzjNnZt8xs/pgvdPM7Gkz+4aZVQbzpgN7MhzSm4EfOed6g+2L5B3l6xHJx3y9L3gsz9D2ioAmYN4w1l0DHJ1ooZkda2a3BOfeNWY2I2pZpnO38nZAhff4chvwBefcu4CTgZ9GFpjZCWZ2Vi6Cit63c+5p59yDuYgjnnTiSXQMnXMbgQOZjCuTP6/Rbss59yfgBqDXOfeIc247cCXQDTzqnGsK1vs7PtF/3DnXEbz9n4E7RvUBopjZZ4HPA382sw3AskxtW2SMKV+nKR/ztXOuHdiVqX0553qAx4e57l58kZ5o+TPA6uDldc656EI7Y7lbeXswFd7jhJlNApYCbQDOuX7gWsCZ2QR8UrccxJWzfWfSMD5H/xjuKxfbugWYY2anADjnGoG/ABdG7ase2BWcexHLnXMbRrnvI5xz/+acW+KcW+GcW+6cuz9T2xYZK8rX2RXCfJ3pUSzSiT/Vun3BY2y/64zlbuXtwVR4jxPB0DwbgbvM7PRg3lbg7/h+VQuBBjO7AsDMvmhmnw2+ZjrSj8vMVpjZN83sw2bWbmYbzewjZnaxmT1mZq8zs9+ZWbOZvTvqfXG3F7PvK83sJ9F93YKvW79tZl82s7+a2WnB/KT7i3q/mdm7zewZMzvPzNab2V4ze8cw9rEkOp4U+xxyDOOYGHx13GJmV8Yc088Hx+a3QaIe7jEbsi8zu9TM/m5m7zCzzWb2lSTbi/ezjxtPMkHLyDPAZVGzpwKvsaBbCXARcHtUnJVAe9Tr5Wb2FTO7zszuM99l5Q1mtsbMXmtmN5tZk5l9yMxOCT5jo5ldkCo+kXyifF04+TqeRJ8x2X7MbKGZ3WS+9fjrw9lPnP3+ycx+bmYvTbGecnc2Oec0jZMJOBbYiv8P9yZgWtSyF4FVwfPjgOeD51Px/41Pxv/Hvg1YESz7LvBw8LwMOAh8Bz8W55uAncm2F7tvfB+3G4GfBPMr8MVcffD6UmA/UJtsf3E+94xgn58O1v0S0AnMT7GP2HiS7jP6GMaJ4UXgG0BV8L5twfw64JdR6z2ET5rDOmYJ9jUJ//Xh9cBpwLnJthfzs48bzzDPr2uA54LnRwFfBrqAS6LOF4ta/yLgZcHzKuAxoCp4fX9wHOqAw8AXgUrgDcHr1+MbBj5KcA5q0jSeJpSvCyJfB8s/P8zjmOxn/RhwdNR73DDPs8jP8wzgc3GWrwr2UxI1T7k7i5NavMcR51slj8cnorcBz5hZvItINgNvNbNS4GXBvBr8L/pcINI/91n8LxTOuW6gBfiNc64Z/4s4I8X2YuPrYvDFGq8BOlzQRxj4Lf4X9pIU+4vdbmSbvwjW/RK+//F5KfYxKJ509pnAH5zvz/cYMDPqM04OWqQ+jD+mMMxjluDzHsQnuN875/7unPtzGttLFM9w3AIsMt/d5F+AHwD/C7zezGqAFhdk5sBLgYeD5xfh/7hFWlFeB/xbcJybgHuc7xe+Bqhzzv3K+a/f15Lez0AkLyhfF0a+jiPhZ0yyn3OBCc65SDxr0tznqcBbnHNfGOb6yt1ZVJLrACQzzKwYmOucexH4oJn9CH8F+GpgZfS6zrkuM5sNfAL478gmnHP7zGwNcA7+P/KF+GLryFujnvcSdFVKtL0EoUZvYwlQGhVXv5m9CMxJtr9UnHOdZvY8UA9MSWMfI95njF4GfrfmAS84566PXWk4x8zM3oovcCOWOf+VtIuONY2fQcJ4UnHObTaztfjuJhOcc8+b2S34VqjLgD9GxR05NyL9B+fjW6gi22qM3nTU8z4G60dd4mScUb4e9PkKIV9HS3gck/xsTsD/kzFSZwNLzKzGOdeaYJ3+4DMod2eZDsr4UUpU/1vn3Dr8V3nLY1c0s7OB9zrnvuT8CBXR3gecan5Yq/XAV1PtOMX2ktmGb0GNHmbJ8H0fR6sM2JTlfQzHbgb3g8bMVqZxzG4HToqa4l4dn8b24sYz/I/DL4ErgOej4isGPoT/WjZiJb6/asRO4J/MrDpqv0n7GYqMY8rXgxVSvk74GZPspzl4Txkjcx3+24L/SrD8sBu4KF65O8tUeI8vHzCzWVGvZwN3Bc+78ReTLMcPXTXBzMrN7Nxg+TQzm4Lva/gj/NdEDxN8dRkoIua/fDOzFNuL3bdFbeM2oBV4Y7Ct+mB/kZbTRPtLZGawzuzgvXcNYx/R8aTaZ/TniBWv9cPwrVjVwB1mdr75i3hS/QyO7AuY6ZzbEDX1RO2iOOp5su1Fxx03HjMrM7PVQTeSZH6F/4r7NwDOuRbgTuAhN3g0k1fiu6FERI73zebHAn8/EBnLNfZnEPtzzusRFkQSUL6mYPL1cI9j3P3gc2kZ/p8z8N9uYGazhpG7Dd/6/FbgMjN7ffDej5rZNfjW/r9Era/cnW3Z6Diuaewn/AUbDjgE/Az4Cf6rqonB8s8BO4DX4r862oxvYXgNvqXkJ/gkdjv+Qo/uYHvt+F/EVwM9+At45gKfCpa/M9n2Yvb9IeABfL+5Y4NlK/Atpd/A93U8LpifcH8JPr8L3v8hfNeH46OWJdrHguh4Uu0z+hjG7Pu1+IsMvxMci8j7rgiWnwU8HfxsvhEc5+Ees9fG+ayXBvu7GZgdzEv2Mx20rQTxzMHf6OEXwzjXvhXz+g3AK2LmXRfnfS8H/gE0AtcE8y7An2vfBmYRDKkGXIX/w/yD4LNemOvfMU2aMjWhfF1I+Xoe8CC+1f6kFJ8x2c/6NfhvGv+Gv9B9E/De4LPHzd34b1V68C3dk4B78EX/u4B348cz/wa+m0vkPcrdWZ4sOIAimNk04Crn3BeD10X4Vph3O+euyWlwKZiZAxY632dSRsjMPu6cG9FQVSIydpSvJZpyd/5QVxOJ9n58t4NJUfPOJuqiuTCK+mpLX2uNkJkVm9lbgFtzHYuIDIvytSh35yEV3hLt+/i+YM+a2XbgbmCHc+6R3IaVWHCxyYeCl+8wP6ydpK8GuM0593zKNUUkDJSvBZS78466moiIiIiIjAG1eIuIiIiIjAEV3iIiIiIiY0CFt4iIiIjIGFDhLSIiIiIyBlR4i4iIiIiMARXeIiIiIiJjQIW3iIiIiMgYUOEtIiIiIjIGVHiLiIiIiIyB/w9Clw7ydSlarQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(12,8))\n", "plt.rc('font',family='Times New Roman')\n", "params = {'mathtext.default': 'regular' } \n", "plt.rcParams.update(params)\n", "\n", "plt.subplot(2, 2, 1)\n", "plt.plot(vehicle.v_kmsc, vehicle.h_kmc, 'k-', linewidth=2.0)\n", "plt.xlabel('Speed, km/s',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 2)\n", "plt.plot(vehicle.acc_net_g, vehicle.h_kmc, 'b-', linewidth=2.0)\n", "plt.xlabel('Deceleration, Earth g',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 3)\n", "plt.plot(vehicle.q_stag_total, vehicle.h_kmc,'r-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-rate, '+r'$W/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "\n", "plt.subplot(2, 2, 4)\n", "plt.plot(vehicle.heatload/1.0E3, vehicle.h_kmc, 'm-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-load, '+r'$kJ/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.savefig('../plots/huygens-titan.png',bbox_inches='tight')\n", "plt.savefig('../plots/huygens-titan.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/huygens-titan.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Radiative heating is neglected in these calculations." ] } ], "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 }