{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 46 - Uranus Probe - Decadal Study - 2010 (Concept)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from AMAT.planet import Planet\n", "from AMAT.vehicle import Vehicle" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook simulates the atmospheric entry of the Uranus entry probe studied during the 2010 Decadal Mission Concept Study (page 42). https://sites.nationalacademies.org/cs/groups/ssbsite/documents/webpage/ssb_059323.pdf" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"URANUS\") \n", "planet.h_skip = 500.0E3\n", "planet.h_trap = -100.0E3\n", "planet.loadAtmosphereModel('../atmdata/Uranus/uranus-ames.dat', 0, 1, 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle1=Vehicle('uranus-probe', 88.0, 190.0, 0.0, 0.46, 0.0, 0.19, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle1.setInitialState(500.0,0.0,0.0,22.3,0.0,-68,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": "iVBORw0KGgoAAAANSUhEUgAAAt4AAAHtCAYAAAAuknnBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzde3xcdZ3/8dcnvbdJadombdPSlmtBKAUJBURoERAXAcGfoLCI7ipxEXVdkR8suCiKqyKCioLEywqyiLKoSEX5gQsoIJcAcpOrtLSkt6SX9EqbNJ/fH98zzWSa2yQzc87MvJ+Px3nMzPmeOfM5M+1nPvnO93yPuTsiIiIiIpJfFXEHICIiIiJSDlR4i4iIiIgUgApvEREREZECUOEtIiIiIlIAKrxFRERERApAhbeIiIiISAEMjzuAQpk8ebLPnj17wNu3tLRQU1OTv4ASqlyPG8r32HXcyffkk0+2untxBJsj2ebslGL6XJNG793Q6P0bvFJ77/rK2WVTeM+ePZumpqYBb19fX5/V9qWiXI8byvfYddzJZ2ZvxB1DoWWbs1OK6XNNGr13Q6P3b/BK7b3rK2drqImIiIiISAGo8BYRERERKQAV3r1oaGiIO4RYlOtxQ/keu45bSok+18HTezc0ev8Gr5zeO3P3uGPoxsyuBC6LHj7r7vPMbAxwFbARqAYudvcN0fa9tqWrr6/3Uho/JCLlw8yedPf6uOPoiXK2iEh3feXsRPV4m9k4YDJwQrR8IGq6HnjI3S8F7gRuTHtaX20iIpInytkiItlJVI+3mX02uvsDd38rWlcHLAaq3X2LmQ0H2oADgO29tbn7kvR9q/dEROLi7nR2duLuDB+e/WRSSe3xTlLOdofOznB/2LChHJWIyNAURY+3mRnwIeBbwCozOytqWgi0uvsWAHfvICTuBf20iYgkwq9+9SuGDx/OmWeeGXcoOZO0nH3LLTB8OHzkI0Pdk4hI/iRmHm8PXe9HmFk1cCHw32a2HpgOrMnYfCNQB3T00dZNS0sL9fVdf3w0NDSU1WB+EYlPRUXo49ixY8eAn9PY2EhjY2Pq4eTcRzU0ScvZ0Vu8s9dbRKSQBpqzE1N4p7j7OuALUW/KZ4F7ga0Zm40C2qP7fbXtVFNTU1KTs4tI8RgWjX3ozKIqTC80zaw1L4HlQFJydqrwTtDoSREpIwPN2YkZatKD7wIzgWbCme/pqoDl/bSJiCRCqvDOpse7CMWas9XjLSLFIMmFdyfwFHA/UBdNQYWZjSAk9wf7aRMRSYQyKbxjzdkqvEWkGCSm8DazOjM728wqop8sLwQuc/eVwCLguGjT44E73L25r7ZCxy8i0ptSLLyTlrNVeItIMUjSGO9q4OvApcBDwLfTppc6H7jKzOYBU4H0M2z6ahMRiV0pFt4kLGer8BaRYpCYwtvdXyD85NhT2zrgvGzbRESSoBQL76TlbBXeIlIMEjPURESkVJVi4Z00KrxFpBio8BYRyTMV3vmnwltEioEKbxGRPFPhnX8qvEWkGKjwFhHJMxXe+afCW0SKgQpvEZE8U+Gdfyq8RaQYqPAWEckzFd75p8JbRIqBCm8RkTxT4Z1/KrxFpBio8BYRyTMV3vmnwltEioEKbxGRPFPhnX8qvEWkGKjwFhHJMxXe+afCW0SKgQpvEZE8U+Gdfyq8RaQYqPAWEckzFd75p8JbRIqBCm8RkTxT4Z1/KrxFpBio8BYRyTMV3vmnwltEioEKbxGRPFPhnX8qvEWkGKjwFhHJMxXe+afCW0SKgQpvEZE8U+Gdfyq8RaQYqPAWEckzFd75p8JbRIqBCm8RkTxT4Z1/KrxFpBgMjzuATGY2EngC+Fd3f8DMxgBXARuBauBid98Qbdtrm4hIUpR64Z2EvK3CW0SKQRJ7vC8CZqc9vh54yN0vBe4Ebhxgm4hIIlREVaG74+4xR5MXsedtFd4iUgwSVXib2TuAFcC66HEdcDZwV7TJfcCpZja7r7ZCxiwi0h8z21l8l1qvd1LytgpvESkGiSm8zWwc8AF3/0na6oVAq7tvAXD3DmAxsKCfNhGRRCnF4SZJytsqvEWkGCRpjPfFwNcz1k0H1mSs2wjUAR19tO2ipaWF+vr6nY8bGhpoaGgYSrwiIgM2bNgw2tvbB1x4NzY20tjYmHo4OW+BDU3e8na2OTv6u4b29gHFLSKSUwPN2YkovM3sRKDJ3VdnNDmwNWPdKCCVWvtq66ampoampqahhioiMijZ9ninF5pm1pq3wAYp33k725w9bly43bJlwE8REcmZgebsRBTewOeBQ8ws9biacNLNNdH9dFXAckJy761NRCRRSnCoSaLydmVluN20aah7EhHJn6SM8T4bODhtWQ58HPgvoC6afgozGwHMBB4E7u+jTUQkUUqw8E5U3h45Mgw3aW+H7duHujcRkfxIROHt7i3u/mZqAXYALe6+FFgEHBdtejxwh7s3u/vK3toKHb+ISH9KrfBOWt426+r13rx5qHsTEcmPpAw16cv5wFVmNg+YCjQMsE1EJDFKrfDuRyx5e9w4aGsLhXd15oAWEZEESGTh7e6z0+6vA87rZbte20REpHCSkLc1zltEki4RQ01EREpd6oqVqQvpSO5VVYXbtrZ44xAR6Y2+AURECqAzurJL2iwgkmOTo5lz12TOFC4ikhAqvEVECkA93vk3aVK4VeEtIkmlbwARkQJI9Xir8M4fFd4iknT6BhARKQANNcm/1FCT1sRd51NEJFDhLSJSABpqkn/q8RaRpNM3gIhIAajHO//U4y0iSafCW0SkADTGO//U4y0iSadvABGRAtBQk/xLFd7q8RaRpNI3gIhIAWioSf5pHm8RSToV3iIiBdDR0QHA8OHDY46kdKX3eEc/MIiIJIoKbxGRPOvo6GDHjh1UVFSo8M6jsWNhwgTYvl3DTUQkmVR4i4jk2bZt2wAYNWqUhprk2e67h9tly+KNQ0SkJyq8RUTyLFV4jx49OuZISt/MmeF26dJ44xAR6YkKbxGRPHvrrbeA0OMt+aUebxFJMhXeIiJ5ph7vwlGPt4gkmQpvEZE8U4934aQKb/V4i0gSqfAWEckz9XgXTqrwXrIk1jBERHqUk8LbzE7KxX5EREpR+qwmSVDKOXuffcLtK69oLm8RSZ6sJpQ1s/cAXwCm0lW0VwBTgDFDDcbMFgCNQC1wk7t/Nlo/BrgK2AhUAxe7+4b+2kREkmDr1q1A4Xu8yzFnT5kC48fDunVhLu+amlztWURk6LK9ksPPgK8BTwOd0ToDzh5qIGa2G3AUcATwDuA3ZrbI3e8Drgf+4O6/iL5IbgTOip7aV5uISOza2toA2G233Qr90mWXs81gzhx44gl46SUV3iKSLNkW3q+7+zWZK83spRzEsh34mrs78DszexbYYWZ1hC+JC6Lt7gPuMLPZ0XN6bHP3JTmISURkyNavXw/AhAkTCv3SZZmz99svFN4vvghHH52rvYqIDF22Y7zPN7MP9LD+1KEG4u5bowSOmY0DXgQeABYCre6+JdquA1gMLOinTUQkEWLs8S7LnH3ggeH22WdzuVcRkaHLtvA+CviRme1IWzqBG3IVkJkdDfweGAGMBaYDazI22wjU9dPWTUtLC/X19TuXxsbGXIUsItKnwfR4NzY27sxXwORBvnRZ5uxDDgm3Tz894KeIiAzJQHN2tkNNvgw0AE3AjmidAR8eTJC9eA34KeHkm28CrwNbM7YZBbRH9/tq26mmpoampqYchikiMjCD6fFuaGigoaEBADNrHeRLl2XOThXezzwDO3bAsGGD2o2IyIANNGdnW3j/HrjT3belrzSzH2YdYS/cfQXwk6hX5iLgz4Qz39NVAcsB76NNRCQRYhzjXZY5e/LkcOn4Zcvg1VfDmG8RkSTIdqjJbcC/m9kxacsC4N/zEFsT0AzcD9RFU1BhZiOAmcCD/bSJiCRCjGO8yzZnH3pouH3ssVzvWURk8LItvP8D+CxwM3BT2vKJoQZiZmPN7JC0VScB33H3lcAi4Lho/fHAHe7e3FfbUOMREcmV1atXA2H4RIGVbc5OzWby5z/nes8iIoOX7VCTbwK/d/eN6StzdBW0AwlTUr0EPAo87u6/i9rOB64ys3mEC0E0pD2vrzYRkditWLECgGnTphX6pcs2Zx9zTLj905/ysXcRkcHJtvBelpnAI1OHGoi7Pw702B3k7uuA87JtExFJglThPXXqkFNltso2Zx98MFRWhjHeK1ZA4f/mERHZVbZDTT6ZucLMzgC+nZtwRERKy8aNG9m8eTNjxoxh/PjxhX75ss3Zw4fDO94R7mu4iYgkRbaF90Qz+w/YOb7vR8CPgCdyHpmISAlYuXIlEIaZmFmhX76sc/bCheH2nntiDUNEZKesCm93fy9wr5ldB/wV2BeYB/R0ZTQRkbIX4/juss/ZJ58cbhctgs7OeGMREYF+xnib2TG9NL0BvApcDcwCTgP+LbehiYgUvzfeeAOAGTNm5P21lLO7O/BAmDUL3ngDHn8cjjgi7ohEpNz1d3LlVcD+wHqgp/6CnxCugjaNMkjiIiLZeu211wDYe++9C/FyytlpzODUU+G66+Cuu1R4i0j8+iu8vww0ufvqvjYys7L42VJEJFt///vfAdhrr70K8XLK2RlOOSUU3nfcAVdeGYpxEZG49DnG293v7i+BR9v9T+5CEhEpHake70IU3srZuzr2WKithZdfhqamuKMRkXKX7awmIiKShQL3eEuG4cPh7LPD/Z/9LN5YRERUeIuI5MmqVatobW2lqqoqlllNJDj33HD785/D9u3xxiIi5U2Ft4hInjzzzDMAzJs3j4oKpdu4HHwwHHQQtLbCL38ZdzQiUs6G/E1gZoeY2S1mdqGZzR56SCIipeGvf/0rEArvpCjHnG0G//qv4f4114B7vPGISPnKRRfMvwDHu/u3gHeY2Zgc7FNEpOil93gnSFnm7LPPDidZPv00PPBA3NGISLnKReH9KeDQ6P5twIgc7FNEpOilerwPPvjgmCPppixz9ujRcMEF4f43vhFvLCJSvrIuvM3sH8zsH6P784F93L0ZwN073X1DjmMUESk6a9eu5W9/+xujRo1i7ty5scWhnN3lk5+Eykq45x7405/ijkZEylFWhbeZfRW4HTgdwN0fB841s7PzEJuISNF66KGHADj88MMZPXp0LDEoZ3c3eTJ8/vPh/sUXa6y3iBRetj3e7wBmAH9OW/d74Os5i0hEpAT8KepSPeaYY+IMQzk7w+c+F8Z6P/oo/OhHcUcjIuUm28L7XndfD6T3E5TNpYdFRAYqIYW3cnaGqir4znfC/QsvhGXL4o1HRMpLtoX362Z2KTDHzD5kZncAFwA6VUVEJLJ69WqampoYOXIkRx55ZJyhKGf34IMfhNNOg40b4SMfgY6OuCMSkXKRVeHt7rcB/wsMA84F1gPvdvfv5yE2EZGidNddd+HuHH/88VRWVsYWh3J2z8zghhtgyhS4/37493+POyIRKRfDs32Cuz8KPJp6bGbDzOw4d//jUAIxs/cD3wLGA7cAF7p7h5lNBr5K+MIYAVzs7u3Rc3ptExGJy29+8xsATjvttJgjyV/OjvZVtHl76lS4/XZ417vg6qthxoyui+yIiORLn4W3mf0W6K+7ZjLQCgw6iZvZTOA0wtjDA4AbgGXA1YQz8i9x98fM7F+ALwOp/om+2kRECq6trY17770XM+OUU04p6GsXKmdHr1X0efvoo+EHP4CPfxw++1lYsQKuvBKGZ90lJSIyMP2ll2ejbV4GOoHzgKeBprRt9gAmDjGOWcA/u3sH8KSZHQgca2YPAXPc/bFou7uAl83sSmBub23uvnmI8YiIDMqtt97Ktm3bOPbYY5k6dWqhX75QORtKJG9/7GNhjPcFF4QL6zz8MDQ2wv77xxGNiJS6/grvbwIWnRWPme3l7pdnbmRmdw4lCHf/c8aqZqAKeBewOG27ZjMDqAeO6qPtwaHEIyIyGO7OD3/4QwDOO++8OEIoSM6G0srbn/gEzJkDH/oQPPQQzJsXph28+GKoro4rKhEpRX2eXOnubakEHpliZsPStzGzhcD8HMd1GPADYDqwJqNtI1DXT9suWlpaqK+v37k0NjbmOGQRKXf33nsvTz/9NJMnT+b0008f0r4aGxt35ivC8JB+xZizIcd5u9A5e+FC+Nvf4LzzoL099H7vuSf853/C+vX9Pl1EytxAc3a2I9n+B3gm+ilxI7A/cCKwS4/KYJnZ3kCLuz9jZg5szdhkFNBOmJe2t7Zd1NTU0NTU1FOTiMiQuTuXXx5S4UUXXTTkq1U2NDTQ0NAAgJm1DnI3ec/ZkJ+8HUfOnjgxDDP5+MfhkkvCjCeXXQZf+1ooyD/5Sdh774KGJCJFYqA5O9vpBO8F3k04MWcO4az0D7j71wYfapeoZ+YTwCXRqmagOq3dgHHA8n7aREQK6o477uCxxx6jpqaGCy64IO5wgPznbCjNvD1/Pvzxj3DvvXDccbBpE1x7LeyzD7zjHWEqwjWZ/fYiIgOQ7QV0cPfl7v4Fdz/V3c8B7jKzQ3MUz+eBq919W/T4XsKXRcpsYDvhRKG+2kRECmb9+vV8+tOfBuCKK65g3LhxMUfUJc85G0o0b5vB8cfDfffBU0+FC+2MGwd/+Uvo+a6thWOOgauuCkNU3Pvfp4hIVkNNzOwnPayuAXYQppUatOjqak3AGDPbE1gI/AlYYWYHuvvzwHuA69x9O9BkZr21iYgURGdnJ+eeey4rV67kqKOO4hOf+ETcIe2Uz5wd7b8s8vYhh8BPfwrf/z78+tdwyy2hR/zPfw7LxRfDzJnwznfCUUeF5cADYdiwfnctImUm2zHe+xDmfk3/234GkHl2e1bM7AvAVzJWv+TuPzGzDwKXm9kSwvy0l6Vt01ebiEjeXX755dx1111UV1dz0003UVGR9Q+J+ZSXnA3lmbfHjYNzzglLW1sYirJoEfzud7B0Kdx6a1gAxo+Hww8PRfvBB4eZUvbdV3OEi5Q78yx+HzOzOe7+csa6/YCPuvslvTwtEerr610nV4pIrrg7X/rSl/jyl79MRUUFd999NyeeeGJeXsvMnnT3+kE8Tzm7AHbsgBdeCHOAp5YlS3bdbvTo0BOeKsQPPhgOOigU6SJSOvrK2VkV3r3sfA/gCXcf0HRXcSmmJC4iydbe3s6FF17IddddR0VFBT/72c84++yz8/Z6gy28e9mXcnYBLF8Ojz0GzzwDf/1ruO2pGIcwbWGqEE/dzpwZxpmLSPHpK2dnO8Z7Md1/shwGTAF+NvjwRESKx/LlyznzzDN5+OGHGTFiBLfccgtnnnlm3GH1SDk7PnV1cPrpYUlZvz4U4Kli/K9/DT3lr78ell//umvbCRNCEZ4qxA8+GA44AEaOLPyxiEjuZDva7F7gVroSeSewyt1fyWlUIiIJ4+7cdtttfOYzn6G1tZXp06dz++23c+SRR8YdWl+UsxNkwgRYsCAsKe3t8NJL3XvG//pXaG2FBx8MS8qoUWHM+Pz5Yfz4/Pmw117qGRcpJtkW3re4+5/SV5jZHDPb291fy2FcIiKJ8eabb3LBBRfw29/+FoATTjiBW265hdra2pgj65dydsKNGAFz54blnHPCOndYsaJ7If700/Dqq/Doo2FJqa4OBXj6kvx/liLlK9vC+z2EqaLSrQLuBw7JSUQiIgmxefNmrrrqKr75zW+ydetWxo8fz7e+9S0+9rGPYcXRzaicXYTMwlCVujo46aSu9evWQVMTPP54WB57DFatgnvuCUvKrFldPeLz58Pb3x5mZBGR+A2o8DazTwFXAlVmdnEPm9yb06hERGLU3t7OTTfdxBe/+EWWLw8XVTzjjDO49tprmT59eszR9U85uzRVV8MJJ4QFQs/4smVdhfjjj4fC/I03wvLLX4btKirCbCpHHAFHHx2WWbPiOw6RcjbgWU3M7CDgI8B3M5q2u/uKXAeWa8V+hryI5F97ezs333wzV155JUuiKSjq6+u59tpreec73xlbXIOZ1UQ5uzzt2AEvvti9GH/22bA+3cyZoQA/9FCYMycss2froj8iuZCz6QTNbIS7t/ewfqa7Lx1CjHmnJC4ivWltbeXHP/4x119/PUuXhlS23377cfnll/PBD34w9oviDGEeb+VsYcuWMEb8kUfgT3+Chx4KM6xkGjkSpk8Py4wZu96fMQOmTQvj0kWkd4OeTtDMPg/c7e5/i1b9XzPL/HvYgGMJlwoWESkK7k5TUxPXX389P//5z9m2bRvQVXCfeeaZDCuy7j/lbOnJ2LFdl7K/6CLo7ITnnw8F+PPPw8svh6W5GRYvDktvzMLJm6mivK6u59uJEzXbikhP+hvjvQcwIe3x24FxwEq6pqcyYGruQxMRyb1XXnmFn//859x666288kqYVc/MOOmkk/jUpz7FiSeeGHsP9xAoZ0u/KirCFTMPOqj7+s2bQ/GdWt58c9f7K1eGEzpXrYKnnur9NUaO7L0oT7/VSZ9SbvosvN39goxVXwWed/ft6SvN7G25DkxEJBfcnRdffJFFixZx++23kz58oba2lnPOOYfzzz+fvffeO8Yoc0M5W4Zi3DjYd9+w9KajIxTfzc3h6pzLl3fdT79tawtX6uztap0p48f3XpSn7k+dquEtUjqynU5wrrt3+xvXzOYAewF/6/kpIiKFtW3bNh588EEWLVrEokWLWJz223lVVRXvf//7+cd//EeOPfZYhg/PNg0WFeVsyanhw8NY7xkz+t5u8+auwry34nz5ctiwISwvvtj7vlLDW9KL8mnToKYGJk/edRk1KrfHLJJLA51O8GjCz5MnRJcgTjeRcNb8XTmOTURkQLZv305TUxMPPvggDzzwAA8//DCbN2/e2T558mTe+973csopp3DSSScxZsyYGKPNP+Vsidu4cbDPPmHpjTusXdt/cZ4a2rJqVThJtD9VVbsW4z0V6al11dWazUUKZ6BdPeuBmwnjB4/KaNsO/CiXQYmI9GXDhg08+eSTPPLIIzz44IM8/PDDbNmypds28+bN4+STT+bkk0/msMMOK7oTJYdIOVsSzwwmTQrL3Lm9b9fREYru9GJ8xQpYswZaWqC1tfuycWNY+jpJNDOOiRN7L9InTQrF+YQJYUndr6rSCaSSvQEV3u7+nJktBM5wdyVsESmYbdu28cwzz/DEE0/w+OOP88QTT/DSSy+RORXq/vvvz4IFC1i4cCELFixg6tTyPX9QOVtKyfDhXbOoHHZY39u6h6Erra09F+WpJb1t7dpQxK9ZE2Z3GaiKiq5iPL0g7+nxbruF8eyZt6U90k16MuCP3N3b6KWXxMze6e4P5SwqESk77k5zczPPP//8zuW5557jueeeo729+1TUI0aM4KCDDmL+/PksXLiQY445pqwL7Z4oZ0s5MgtF7W67wV57Dew5HR2wbl3fxXpbW9hm/fqwrFsHmzaFon3t2sHHO3ZsKMB7Ksr7u03dr6wMfwRIcehvHu/fApX97MOAA4DaXAUlIqWrs7OTZcuW8eqrr/LSSy91K7Tb2tp22d7M2H///TnssMOYP38+hx12GAcddBCjR4+OIfpkU84Wyd7w4WF4SU0N7L//wJ/X0dFViKeK8d7up04ibWvrfrtlS1hWrhx8/GZh2MtAivS+bseM0dCZQuivx/tlwljBZrrmgM1kwMm5DEpEipu7s2LFCl599VVeffVVXnnllZ33X3vttZ0Xq8k0adIk5s6dy4EHHrhzmTt3LuPHjy/wERQt5WyRAhk+vGsc+GC4h9lfeirIU7d9taVuN23q2vbNN4d2PKlCfOLE3pdJk7o/rq4O87bLwPRXeH8D2OjuPX9LRswsZ2fHm9kJwJXAB919SbRuDHAVsBGoBi529w39tYlI/rS1tbF48WKWLFnC4sWLd7mfPqtIpqlTp7Lvvvuy7777diuya2trMXW5DIVytkiRMAvDRCorwxSJg7VjRziZNJtivacC/623uobODPTE1JTKyr4L9dra8ItCbW3X/XL90bK/C+i0DnA/p5ODs+TNbArhZ9L5GU3XA39w91+Y2XuAG4GzBtAmIoPQ2dnJqlWrWLZsWbclVVgvWbKEdevW9bmPiRMnsu+++7LPPvt0u917772pqqoq0JGUF+VskfIzbFjXSZxDsW1bKOBTw2NSRXjq5NP0x5nLpk1hWbp04K83fnz3Yry2NswPv+eeYXz+/vsP/ZiSyDJnBsjqyWYzgQuAT7l7Ti78amYVwA5gD3dfYmZ1wGKg2t23mNlwoI0wRnF7b22pnpeU+vp6T79inUi5cndaWlpYtmwZb7755i7F9bJly2hubqajo6PP/YwZM4Y99tiD2bNnd7tN3Z84cWKBjqj0mdmT7l6fg/0oZ4tITrmHgr23ojx10urq1d2Xfr5iAPjIR+DGG4vvokh95eysJ7Ixs2HAqcB5wAnAKsKYwpxw986Mn5oXAq3uviVq74guCLEAaO+jbUmuYhIpFu3t7axcuZLly5ezYsUKli9fTnNzc7ei+s033+x1jHW6mpoaZsyYwe67775zmTVr1s7iuqamRsNCioBytojkk1nXCZyzZw/sOe5hmEtmMf7GG/Doo/DAA2G7m24KPd9vf3vXEJUJE4r7RNABF95mtgfwceCfgEnAG8ApwD3AkXmJLpgOrMlYtxGoAzr6aBMpCe7Ohg0bWL169c4lVVSnF9grVqygpaVlQPusrq7eWUxnFtepdZo1pLgpZ4tIUpl1DY/Zd9/ubS+/DPPmhaEvAJdcsuvzhw0L48qrqrpmdEndz1z6a6usLOx86v2+lJmdQegpeRfwGnANcBNwqbv/IdrskbxFGM7M35qxbhSh54R+2nZqaWmhvr6r17+hoYGGhoYchikycJs3b2b16tW0tLR0K6jTH6fut7S0sH379gHtt6KigqlTpzJt2jTq6up23mYW1ZWV/c04J3FrbGyksbEx9XDA8yYoZ4tIMZszB159FRYtgr/9rWuYSuq2rS0U5W1tYcmFMWPCSaBTpoSe9SlTYL/94MQT4eCDB7aPgebsgdT4Wwnj8h4BPu7urwCY2eAHh2enmXDme7oqYDkhwffW1k1NTQ0aLyi55u5s3LiR1tZW1gV36K4AACAASURBVKxZM6ClpaVll8ub96eqqora2lpqamqora3dpbBO3dbU1JTbpdFLVnqhaWYDPWkSlLNFpMjtvjucf37v7e3tYVx55rJhQ8/r+2vbujVMxZg5HeMll8Bf/gJHHNF/zAPN2f0W3u6+CFgUnZRznpnNBu4ERqS2MTPzoZyl2bf7gR+a2Rh332pmI4CZwIOEE3p6axPJSkdHB2vXrmXNmjVZFdL9nYTYk9GjR1NbW9utmO7pcU1NDTU1NYwZMyYPRyylSDlbRErdiBFd0xX2Ztu2rgsZrVvXdXXSFSu6L8uXw9//3vM+pk0b+BVQByqbS8YvBf4jOgv9NGCGmX0F+DHwbqCxr+cPlHWdpWPR6640s0XAccAi4HjgDndvjrbvtU3Kz44dO2hra2PdunWsW7eOtWvX7ryf+Tjz/qZNmwb1mpWVlUyaNGnAS21tLePGjdOJiZJXytkiUizcwzziqWkJN2/u+35P61JTIaYK7a2Zg9r6MWVKGG++zz5hmT8fjjwyDEPJpayHk7t7B/A/wP+Y2d7AvxHGEw45iZtZFXBO9PCjZvY9d28BzgeuMrN5wFQgfaBfX21ShFLDN3orkvsqoNva2hhsR15FRQUTJ07MqoieNGkSo4ptniMpK8rZIjJQ7mGav23bQiH81luhgN26Nf/3c/0b3PDh4aqaqWXChDCOe9q0XZfp08PJloUwpHm8d+7E7Fx3vzkH8eSN5oQtLHdny5YtWRXNqcfr169nx44dg37tCRMmUF1dvXOZOHFij/czH1dVVVFRUZHDd0EkN3I1j3fa/pSzRRKiszMUuqnlrbe6Py7UutT6vA1C68fo0WGGkXHjuq7ombrf07qe2idM6Cq0x46Nb8rBnM7j3ZOkJ3AZPHdn06ZNO8c+r127dpf7PT1eu3Yt7e27TFQwYJWVlX0Wyb217bbbbjq5UKQfytlSztJ7dXsrRNMf99WWiwJ4EKcJ5dWwYeGCNaNGhWEWqWX06P7vD3S7np5TLl/dBZy5UOI0mAI6dX8wJw9CuLJhf73MPbVNmDCBESNG9P8CIiJSFNxh+/bBF7O5ftzZGfc70t2oUaH4TBW86UtP6/O17ahRhZ3Tuhzp7S1SHR0dtLa2dpvrOf1+6jY188ZQCuhx48btHPs8ceLEXe739Li6uloXYBERSYAdO7rG0fa0pIrSfBe7STJ8eN9F6FAfZ7N+xIjivQqjZE+Fd4K0t7ezatUqmpubd16VsKeLqaxevZq1a9dmvf9sC+hJkyZRXV2tkwdFRHLMPRS9qZkZUrMzbNnSe4HcV1tf7QO8/lbejRyZ32I3m8flMqxBkkeFd4Fs376dpUuX8vrrr7Ns2TKWL1++s8BO3a5atWrAM3KY2c6p6VLzPvd0O3nyZBXQIiI5kCqWU3MCp5b168NUZqkCOv1+X+sKNdzBrPt42rFjdx1fO5jiNZttR44EnbsuosI7p7Zt28ZLL73Eiy++yOuvv95tWbZsGZ39ZFkzY+rUqdTV1TF9+nSmTZvG1KlTeyyoJ02apJMIRUSGwB3WrIGVK7suppG639Kya4G9bl1ue49Hjw5TmKXP0JBeEPdUJA+mbeRIDWUQSQoV3oO0ZcsWHnvsMf7yl7/w7LPP8txzz/HKK6/0Oo7azJg5cyZ77rknM2fOZPr06TsL7NTtlClTdFKhiEiOrV4NTzwBr7wCr70Gr74alubmcOnpbIwa1X1u4NSSKqDTC+m+Ho8bp5PYRMqR/ttnYfHixdx2220sWrSIJ554Ypfp8syMffbZhwMPPJC9996bPffckz333JM99tiDWbNmMXLkyJgiFxEpH52dcP/98POfwwMP9H45aAjz/qYuojF1atdtbW3PBXaur2InIuVFhfcAvPDCC1xyySUsWrRo5zoz4+CDD+boo4/mkEMOYe7cubztbW9j7NixMUYqIlLeXnwRPvxhePLJrnVjx8Jhh8Hb3tZ1Oeh99oFZs8JwDxGRQlHh3Y977rmH008/na1btzJmzBhOO+00zjzzTBYuXMiECRPiDk9ERCIbNsDxx8Py5eES0B/7GLzvfXDQQRrWISLJoFTUh82bN/PhD3+YrVu38uEPf5irr76a2trauMMSEZEe/PKXoeg+9FB48MEwjlpEJElUePfh7rvvpqWlhUMPPZSf/vSnVGguJBGRxHrhhXB75pkqukUkmVRJ9uGxxx4D4NRTT1XRLSKScJs2hdvq6njjEBHpjarJPjz77LMAzJs3L+ZIRESkP6nLkmsCKRFJKhXefViyZAkAe++9d7yBiIhIv1IX/tUPlCKSVEpPvejs7GTp0qUAzJo1K+ZoRERERKTYqfDuxerVq9m2bRsTJ06ksrIy7nBEREREpMip8O7FypUrAZg+fXrMkYiIiIhIKVDh3Yu1a9cCMHHixJgjEREREZFSUBLzeJvZGOAqYCNQDVzs7huGsk8V3iIi+ZGPnC0iUgxKpcf7euAhd78UuBO4cag7VOEtIpI3Oc/ZIiLFoOgLbzOrA84G7opW3Qecamazh7LfVOFdrSsxiIjkTL5ytohIMSj6whtYCLS6+xYAd+8AFgMLhrLTLVu2AFBVVTXE8EREJM1C8pCzRUSKQSkU3tOBNRnrNgJ1Q9nptugSaCN1CTQRkVzKS84WESkGpXBypQNbM9aNAtrTV7S0tFBfX7/zcUNDAw0NDb3udPv27YAKbxGJR2NjI42NjamHk+OMJcfykrOh64qVO3bkIEoRkSwMNGeXQuHdTDgrPl0VsDx9RU1NDU1NTQPeqQpvEYlTeqFpZq0xh5NLecnZALvtFm7XrRt8cCIigzHQnF0KQ03uB+qi6akwsxHATODBoey0s7MTgGHDhg01PhER6ZKXnA0waVK4jc6NFxFJnKIvvN19JbAIOC5adTxwh7s3xxeViIj0JJ85O1V4r8kcQS4ikhClMNQE4HzgKjObB0wF+h4IKCIiccpLzlbhLSJJVxKFt7uvA86LOw4REelfvnK2Cm8RSbqiH2oiIiICkLrQsMZ4i0hSqfAWEZGSoB5vEUk6Fd4iIlISVHiLSNKp8BYRkZJQVQXDh8PmzRBdfFhEJFFUeIuISEkwg8nR9eJaS+mSQyJSMlR4i4hIyZgyJdyuWhVvHCIiPVHhLSIiJWPq1HC7cmW8cYiI9ESFt4iIlIxU4b1iRbxxiIj0RIW3iIiUjGnTwq16vEUkiVR4i4hIydBQExFJMhXeIiJSMlR4i0iSqfAWEZGSocJbRJJMhbeIiJQMFd4ikmQqvEVEpGSo8BaRJFPhLSIiJWP8eBg9GjZtgo0b445GRKQ7Fd4iIlIyzGD33cP9ZcvijUVEJJMKbxERKSmzZoXbN96INw4RkUwqvEVEpKTMnBluly6NNw4RkUwqvEVEpKSo8BaRpBoedwAAZmbAGcAV7r5/Rttk4KvAemAEcLG7t/fXJiIi+ZPkvK2hJiKSVEnp8Z4FVAP79dB2O/ATd78YeAX48gDbREQkfxKbt9XjLSJJlYjC292XAPdmrjezI4A57v5YtOou4NNmNq6vtkLELCJSzpKct9XjLSJJlYjCO9LZw7p3AYtTD9y9Obpb30+biIjkXyLz9owZ4ba5GTo6crlnEZGhSVLh3ZPpwJqMdRuBun7aREQkHrHn7VGjoK4OduyAN9/M5Z5FRIamICdXmtnXgLm9NN/n7t/upc2BrRnrRgHt/bTtoqWlhfr6rk6VhoYGGhoaeo35oIMO4tRTT2X27Nm9biMiki+NjY00NjamHk4u9OvHnbezzdmZTj4ZNmyAzp765EVEcmygOdvcvTAR9cPMZgOL3d3S1v07cKy7vzt6bMBbwLHAgt7a3P2RzP3X19d7U1NTvg9DRCTnzOxJd0/cMLp85m3lbBEpVn3l7KQPNbkXmJP2eDawHWjqp01EROKhvC0i0oskFd4GO3tAAHD3JmCFmR0YrXoPcJ27b++rrZBBi4iUMeVtEZEsJOUCOrXAR6KH55vZze6+KXr8QeByM1sCVAKXpT21rzYREckT5W0RkewlZox3vmm8oIgUq6SO8c4n5WwRKVbFPMY7NmlnppaVcj1uKN9j13FLKdHnOnh674ZG79/gldN7p8K7F+X0jyBduR43lO+x67illOhzHTy9d0Oj92/wyum9U+EtIiIiIlIAZTPG28xagDeyeMpkoDVP4SRZuR43lO+x67iTb5a718QdRCENImenFNPnmjR674ZG79/gldp712vOLpvCW0REREQkThpqIiIiIiJSACq8RUREREQKQIV3BjMbY2bXmdl/mtkNZjY+7pgKycyuNDOPlmfijidfzOwEM3vMzGanrSuLz76nY4/Wl+Rnb2bvN7PFZrbGzL5jZsOj9ZPN7EYz+4aZXWNmI+KOVQavXP7/DkW2eU/vaTCYHKL80sXMFpjZy2a2zsy+nba+LP/tqfDe1fXAQ+5+KXAncGPM8RSMmY0jnOBwQrR8IN6I8sPMphCumDc/o6nkP/vejr1UP3szmwmcRjiefwM+Dnw2ar4d+Im7Xwy8Anw5liAlV0r+/+9QDDLvlf17OoQcovwCmNluwFHAEcA5wAVmdnzUXJb/9nRyZRozqwMWA9XuviX6q7YNOMDdl8QaXAGYWSqZ/MDd34o1mDwzswpgB7CHuy8pp88+89ijdSX52ZvZ0cBf3L0jenwVcADwFeBX7l4XrZ8OvAxMcffNccUrg1NO/3+HIpu8B2zvra2c3tPB5BBgbm9t5ZZfzGwM8JZHxaaZPQl8nvB+lOW/PfV4d7cQaHX3LQDRf7TFwII4gyoEMzPgQ8C3gFVmdlbMIeWVu3dmrFpImXz2mcdeyp+9u/859YUZaQaWAu8ifL6p7Zqju2V1WfYSspAy+f87FFnmvb7aysYgc4jyS8Tdt6YV3eOAF4EHKON/eyq8u5sOrMlYtxGoiyGWgvLgCMJwg+uA/zazf4g5rELSZ18en/1hwA8o48+7ROnzHJy+3je9pz0bSA7Re5ch+uXg98AIYCxl/P4NjzuAhHFga8a6UUB7DLHEwt3XAV+IekE/S/iPUg702Zf4Z29mewMt7v6MmZX9511i9HkOTn/vm97TNFnkEP173NVrwE+Bq4BvAq9Tpv/21OPdXTNQnbGuClgeQyxx+y4wM+4gCkiffZeS++zNbBjwCeCSaFW3zzv6Y2Mc5fl5lwL9/x2cvt43vadpsswhyi8Z3H2Fu/+EML57AWX8b0+Fd3f3A3XRyQBE0//MBB6MNap4dAJPxR1EAemz71KKn/3ngavdfVv0+F5gTlr7bMIJPU0FjktyQ/9/B6ev903vaXfZ5BDll941EQrrsv23p8I7jbuvBBYBx0WrjgfuSDsxomSZWZ2ZnW1mFdFf5xcCl8UdV75ExwhgUF6ffeaxl/pnb2aXAk8DY8xsTzP7Z2A9sMLMDow2ew9wnbtvjytOGbxy+v87FNnkPb2nXbLNIe7e1FtbwYOPmZmNNbND0ladBHynnP/taTrBDGZWTRiDtASYCvy7u2+KNagCMLMDCGN6NwAPAd9295fijSo/zKyKMJ/o9YS5Vb/n7i3l8Nn3dOxALSX62ZvZFwjTfqV7yd33N7NZwOWEz7sSuCxj9gIpIuXw/3coBpP39J4OPocovwRmNh/4HfAS8CjwuLvfHrWV5b89Fd4iIiIiIgWgoSYiIiIiIgWgwltEREREpABUeIuIiIiIFIAKbxERERGRAlDhLSIiIiJSACq8RUREREQKQIW3yCCZ2XFm9oiZLexnuwPN7Jdm9h85fO1v5mpfIiJxMLODzOz2XObGXl5nkZl9KJ+vkS9mdryZ/aW/7xkpHsPjDkBkqMxsf+CrwFrgMOAg4MPufkueX3oZMG8A2y0FqoFhuXjR6PK57bnYl4gIgJkdBnwNOAb4f8BiwoVLRgA3uPs9eXjZ14DdyFFu7MP1wPP52HGUj98P/AB4BfhL1DQK+D/AfHdfMoSXGOj3jBQJFd5SCn4DnOnuz5hZBfBfhXhRd3/FzNYMYLsNZrYihy+9APhTDvcnImXO3Z8ws9uA/dz95NR6M3sv8Gsz+7y7X5/j19yS49wIgJmdAjyXKnjd/e5cv0aKu7cDvzCzTwMPuPsX0uK4czD7NLNa4Fh3/4W7v2xmrTkKVxJAQ02kqJnZRGBfYDOAu3cSLtNbqEuyduZ4u4E4Frg/h/sTEQHY5ZLm7v474EvAtWY2NQ+vmdNcbWYzCD3chdbT5eDvA5qz2YmZjQJuBsakrc7l94fETIW3FDV3Xwu8DPzBzI6I1r0BPG5mE8zsK2b2uJmdamYrzexFM3t76vlmdpaZXWZmD5jZf6atrzezL5nZL8zsV2a2W7R+tJl938wuN7MfE34mzYqZ/d8opjPN7DAz+020v69EMf7CzGqi9W1m9sWMXYx0923Rvi40s8+Z2cNm9vWs30ARkf79mDDk5AwAM9vLzL5sZj8ysz9GxS5mNirKZZeZ2Z/N7IOpHZjZJ83sP8zsUTP7VG8v1FNONrMPRDnzo2b2qpl93czqzOy/onz6ZzN7V7SLs4EZwMVmdlp0Ls59ZvaRtNc4z8y+amY3m9lNZjbOzHYzsy+aWZOZvcvMXjaz181s9mDesOg9WRj1iGNm7zez75jZt8zsD2ZWGR3Dd8zsruj7ZkX0Hr8NOMPMPpe2y+rou2ijmZ3Xy2tWmNk3zKwheg03s1+nPh9JCHfXoqWoF+AA4A1Cr8BPgNpo/TBCEltLGGs3AbgXeIHwR+eRwBeibasJvebvBcYDt6Xt/2Hgm9H9a4GG6P4UYAchufYX408JvUa1hLGAY6P1o4B7gDuBOsIXRjvwRWAscGj0OLV9avx66rivj+5XAZfE/Vlo0aKleBfgo8CbvbStAr4b5dVfAcOi9bcCt0f3rwHeH93/Z6Alun8WcE50f98oV8+NHv8U+FJ0v7ecPBFYD3wbmA8cD1wNXB5texFwV1qsDsyO7o8BXgQ+Gj0+Ebg1bdtFwI2EobenR98X74+O8w/AVwf43j0QfVd8HfgmYbz38WntrcDM6P4LhO+kUcBlwJvA4cDHorz/QCreaPslwLeitrOApb3E8HFgUXR/SvQ+vzvuf1daui8a4y1Fz91fMLO5wJXAJ4FTzOx4D2O+W4AN7n4HgJl9CXgI2JPwJdNpZp+NdnUPIcG/F5iUtv7F6LmVwPnRc3H3VWa2LItQJwPXAR9z9y3RPraZ2Upgsbsvj15nJfCgh/GPTxO+ECYTTtJ8L+FLAmALcK6ZvQLcQPijQ0QkH5xQjM4HZgGfNjOANcBuFh58gtBpAHATYagFhFzblJFrpwHPZbzGR+khJ7v7WjNrA+5098cBzGwpsNnMphF6iCt7DNp9q5mtTlv1qWi/KT8CfhmtX0/4vvhV9BpNUZwDdb9HY7zN7BeEzp6UM9x9qZkdRSi4K6P8vxz4u7s/BjwWPbenfd8VfSc80UdMbwe2Rse9yszW0n3IiiSACm8pamY2DNjdw0k0n4mGf9wNNBJ6EDK9EN1OAGYCP3P3W6N13472eTGhEP52xmsdSkiYGwcZ7r7AO4G9gGfS1meOcdyxs8G9M0rCqWFhkzwMr8HdF5vZR4HvAf8G/COwGhGRHDKzCYRf614g5M01PeTHKYQe2ZEA7r6D0FlA9Jwr3P2R6HG356bpMSdHnO658k3gYuBV4FFgjz4OIf15+wD/m/b474RhNLXsmos7GPyMK08R/kBJWW9m3wBuJ8Seqq4zj6s/HfReu/0euMHMRhM+i62E3nNJEI3xlmK3c9whgLs/A1wK7NfL9iMJSe41YAXwgfRGMzs8Wv9eMxuTsX5D9HD/Qcb6CPAd4L+jxJgVM6sj7UQdM5tOmNFlX8LJlrf28lQRkaH4COEXttsJ+fGoqKcZ2JkfW4HtwHvS1s81syoycq2ZDY86MjL1lpN78gNgiYdpY7dlcSxLCT3kO18CaCMMpckZd++MOkcOiP5wuRv4srs35fJ1Ml7zLuAuwi+zZwEL3L0tX68ng6PCW0rBp6OiNGU6YWxeykQzGxndX0gY37ce+Dlwupl9z8wWRr0RmwkJchywyMxOjE5k2Y9QrD8NfMXMRprZZMLY6lozG2Zml5jZOb3EaNFyOeHL6Zoe2jK3z3x8MmE8Yso+wGnuvgH4dA/PERHJxi49qWZ2DHAF8Al3byH0Lq8C7jazU6ITKI+OerhvI8x+cpqZHQd8xN03EnLtZ8zsCjNbQOiASA3TS89/veXklPTe50OAGjMbCxwFjDGzPaO2dsLJiKkOmPTX+AHwf8ysOnp8OPBDd+8g1EQ95mIzOyKKa1Qv710FGTVV9MfJewm/clYThjDOAXYHxqaduJnZq769h/i7B9XDeJTosxgP/JHQq5/NHyRSIBpqIqVgd+AFM/sd4Wc4A/4lrb0DuMzM1hN6hz8N4O73mtlnCD9XnkY4OfF5ADN7H/B9whfJT4CL3N3N7EzCVE8vAL8g/GR4MOGCE2cShrB0u3CPmc0nnDQ0k/DT42+AK6J47gGOAGab2V6EnpipwJlm9hJdvT8fBca5e2PGsf+XhQtfdBBOZhIRyVqUR84mFIc/BFoIJ+jtBvyDu/8FwN23R/nxh4RcdydhbDfAvxI6LW4m/MKXmknkx4Tcdz5hSNz57r46KiyPAPYys1t6y8lm9gHCuOaPm9lL7t5MONHzakLRfTPwPmAu8DrwM0Lv/LkWZrs6ADjJzO5291+Z2T7AHWb2MOFX08uinvkzgKnR8b1AuJjQNDM7IIrzQ4Rx679Je99GEHL/PGB6VJg7YWz1+6N9/hV4EHiCMHzm94ScvohwQuc8MzvN3VP7/QXhYkabokJ9KvABM1tM6Mkmen7mNSuWEk4ePQ0YDVSY2X2EEywLNcWu9MP0WUgps3CZ3Z+6++wCvd5F7q7LuYuIlJho2Msod0/kBczM7CLCDDNLosejCb9WfCN1bpDET0NNpNT1NIwjPy8UfpJ9qhCvJSIihROdPPr2BBfdo4EvALPTVk8GNqvoThYNNZGSFY2vO4fw0+EZ7n57nl/y79FPoCIiUlqGu/sNcQfRG3d/y8Jl678bnXi/nDAs5gt9P1MKTUNNREREREQKQENNREREREQKQIW3iIiIiEgBqPAWERERESmAsjm5cvLkyT579uysntPS0kJNTU1+AopBKR2PjiWZdCz58eSTT7a6ezKCKZDB5GxI1udWCDre0lZuxwulccx95eyyKbxnz55NU1N2V2qtr6/P+jlJVkrHo2NJJh1LfpjZG3HHUGiDydmQrM+tEHS8pa3cjhdK45j7ytkaaiIiIiIiUgAqvEVERERECkCFdx8aGhriDiGnSul4dCzJpGORuJXb56bjLW3ldrxQ+secuAvomNmVwGXRw2fdfZ6ZjQGuAjYC1cDF7r4h2r7XtnT19fVe7GOGRKQ8mdmT7l4fdxw9Uc4WEemur5ydqB5vMxsHTAZOiJYPRE3XAw+5+6XAncCNaU/rq01ERPJEOVtEJDuJ6vE2s89Gd3/g7m9F6+qAxUC1u28xs+FAG3AAsL23Nndfkr7vvPSevPUWjB6d232KiGRIao930eVsEZECKIoebzMz4EPAt4BVZnZW1LQQaHX3LQDu3kFI3Av6acsfd7j8cth9d2huzutLiYgkUdJyduudrTTVN7H4i4uHuisRkbxJzDzeHrrejzCzauBC4L/NbD0wHViTsflGoA7o6KOtm5aWFurru/74aGhoGPwAfjN48UVobYUrroDGxsHtR0SkF42NjTR25ZbJccbSk6TlbN/hbHpyEyNrRw7iaEREhmagOTtRQ03SmdlXgXrgXuAMdz88re0p4NboYY9t7n51+v5y/rPlyy/DAQeE+y+8AHPm5G7fIiJpkjrUJF3cOXvLK1t4fM7jjJoxiiOXHTmEIxERGZqiGGrSg+8CM4Fmwpnv6aqA5f205decOfBP/wQ7dsCVV+b95UREEi7WnD1mrzFUjKlg25vbaF/fPtTdiYjkRZIL707gKeB+oC6aggozG0FI7g/205Z/l10Whp388pewcWNBXlJEJKFizdk2zBj7trEAbH5+81B3JyKSF4kpvM2szszONrOK6KSdC4HL3H0lsAg4Ltr0eOAOd2/uq60gQc+eDYcfDtu3wyOPFOQlRUSSIIk5e9yB4wDY/JwKbxFJpsScXEn4+fHrwKXAQ8C306aXOh+4yszmAVOB9DNs+mrLv/32g0cfhWXLCvqyIiIxS1zOrjyoklWsYtMzm3K1SxGRnEpM4e3uLxB+cuypbR1wXrZtBZGax3v79thCEBEptCTm7Kr6KgA2Nmnon4gkU2KGmhSt1KwwZvHGISJS5ioPqQSDzc9upnNbZ9zhiIjsQoX3UI0fH25feCHeOEREytzwquGMnTMWb3c2PafhJiKSPCq8h+qss0Jv9w03wD33xB2NiEhZ03ATEUkyFd5Ddcgh4fLxnZ3woQ/Bq6/GHZGISNlKFd6bnlSPt4gkjwrvXLj8cjjtNFi/Ht73PtiwIe6IRETKUqrw3vCE8rCIJI8K71yoqICbbw6XkH/xRbjqqrgjEhEpS5WHVGIjjc3PbqZ9ra5gKSLJosI7V6qq4PvfD/dvuikMPRERkYIaNnYY448cDw7rH1gfdzgiIt2o8M6lY46B6dPhzTfhqafijkZEpCxVH1cNwLo/ros5EhGR7lR455IZnHpquH/nnfHGIiJSplR4i0hSqfDOtdNPD7e33AI7dsQbi4hIGao6rIphVcPY+vJWtjVvizscEZGdVHjn2nHHwZ57wpIl8LvfxR2NiEjZqRhRwYQFEwBY+//WxhyNiEgXFd65VlEBF1wQ7l93XbyxiIiUqYnvnQjAmrvWxByJiEgXFd758E//BGPHwn33wd/+Fnc0IiJlZ9LJkwBYe89adrylYX8ikgwqvPOhuhrOPTfcv+aaeGMRESlDo2eMpvLQSjq3dLL+fzWtoIgkgwrvfPnc58Is6lPBsAAAIABJREFUJz/7GaxYEXc0IiJlZ/KpkwFo/W1rzJGIiAQqvPNln33CZeS3b4fvfS/uaEREys6kU8NwkzW/XYN3eszRiIio8M6viy4Kt9dfD21t8cYiIlJmKudVMmrWKLav2E7bw8rBIhI/Fd75dOSRsGABrF8P3/lO3NGIiJQVM2PKWVMAWH3r6pijERFR4Z1/V1wRbq+5BtbpKmoiIoVUe3YtAKt/uZrO7Z0xRyMi5S5xhbeZjTSzZ8xsYfR4jJldZ2b/aWY3mNn4tG17bUuMBQvgXe8KQ02uvTbuaEREci7JebtybiXjDhxHx9oOXUxHRGKXuMIbuAiYnfb4euAhd78UuBO4cYBtyZHq9f7ud2Hz5nhjERHJvUTn7Z293v+t4SYiEq9EFd5m9g5gBbAuelwHnA3cFW1yH3Cqmc3uq62QMQ/IO98Jhx8eer1vvTXuaEREcqYY8vaUs6eAQcuvW2hf057PlxIR6VNiCm8zGwd8wN1/krZ6IdDq7lsA3L0DWAws6KcteT71qXB76aXQ3BxvLCIiOVAseXv0rNFUv7sa3+asvHllPl9KRKRPiSm8gYuBr2esmw6syVi3Eajrp20XLS0t1NfX71waGxtzEHIWzj4bTjgBWlvh2GNhyZLCvr6IFJXGxsad+QqYHHc8vchb3s51zq77RHiJ5Tcux11zeotIbg00Zw8vXEi9M7MTgSZ3zxyA58DWjHWj4P+3d99xctX1/sdfn63Znt4rJQWIkGTpAvkhKKKiYENU9KpEFMu9NsQuYrkoioJ6jYrlqoiIAnIBK1UIIUACIaQ30usm25LsZj+/P75ns7OT7Tt938/H4zxm5pyZcz5nJvnMZ7/z/X4Prb8VdrWtnREjRrBo0aL+htp3eXnw29/Cq18NixeHrid33AFz56YvJhHJWPPmzWPevHkAmFnGXXox2Xk70Tl72OuHUTSmiMYVjex7bB+Dzx2csH2LiPQ0Z2dKi/engNvMbFcU7ATCoJsyYEjccyuALcDmLrZlphEj4OGH4VWvgh07wu2NN4JaX0Qk+2RV3s4rzGPM+8cAsOl7m5J9OBGRDmVK4X0FcErMsgX4APALYKyZlQCYWSEwEXgEeKiLbZmrqgr++tfQ17ulBa69Fi67TFe2FJFsk3V5e+w1Y7FiY9fdu6hbWpeKQ4qItJMRhbe773T3Ta0LcBjY6e4bgfuAV0VPvQC4y903u/u2zralOv5ey8+Hr38d7r03FOJ33w2zZ8NTT6U7MhGRHsnGvF08upgxHwit3hu/uTEVhxQRaScjCu9ufAh4o5l9HrgYmNfDbZnvDW+AZ56BWbNg7Vo4++xQkB8+nO7IRET6I2Pz9sTPTMQKjB2/36FWbxFJORsoo7urq6s9rYMru3LwYOh68t3vhsfnnQe33QbHHJPeuEQkI5jZM+5ene44UimZOXvVR1ex+dbNDHnNEE5+8OSkHENEBq6ucnY2tHjnvuJiuOmm0Pd79Gh45BE48UT45jfh0KF0RyciklMmfXkS+VX57P3rXnY/GD+7oYhI8qjwziSvfjU8/zy8611w4EBoBZ85Ex58MN2RiYjkjKLhRUz+4mQAVn9sNYcb1b1PRFJDhXemGTEC/vd/4e9/h2nTYOVKeO1r4ZJLdMVLEZEEGfeRcZSeUErjqkbWf3l9usMRkQFChXemuuCC0Pr9ne9ARQX85S9wzjm64qWISALkFecx/bbpkAcv3/Qy+5/en+6QRGQAUOGdyYqK4JOfhBUr4LTTYN06OPdcWLMm3ZGJiGS9ytMrmfCJCdACy69czuF6dTkRkeRS4Z0NxowJXU/OPhtefjm0fP/73+mOSkQk602+fjKlM0ppWN7A6v9cne5wRCTHqfDOFpWV8MADYarBrVth7tww/eAAmQ5SRCQZ8kvyOeGOE7BiY+vPtrLjjh3pDklEcpgK72xSURFavj/5SWhuDreXXAI7d6Y7MhGRrFU+s5zjvnccACvmraBxbWOaIxKRXKXCO9sUFoYBl3/6EwweDPfdByefDP/4R7ojExHJWmOvHsvwS4dzeP9hll66lOa65nSHJCI5SIV3trr0UliyJPT33roVLrwQPv3pcBVMERHpFTNj2m3TKJlaQv3z9Sy/cjneoq58IpJYKryz2cSJ8NBDcP31kJ8fWsJPPx1efDHdkYmIZJ3CwYXMvHcm+VX57PrzLtZ/dX26QxKRHKPCO9vl58MXvwiPPQbHHBNawaur4ZZbNPBSRKSXSqeVcsLtJ0AebLh+A1t/uTXdIYlIDlHhnSvOPBMWL4b3vS9cbv5jH4PXv14DL0VEemnYa4dx3M3RYMsPrGD3/+1Oc0QikitUeOeSigr4+c/hrrtgyBC4//4w8PLPf053ZCIiWWX8R8cz8bqJcBhefOuL7FuwL90hiUgOUOGdiy67rP3Ay8suC4MxN21Kd2QiIlljytenMPp9o2lpbOGF175A7TO16Q5JRLKcCu9cNWFCGHh5662hJfzuu+GEE+CHP4TDuiyyiEh3zIypP5nK8EuH01zTzJILl1D7rIpvEek7Fd65LD8frrkGli2DN70JamvhIx+BV74SXngh3dGJiGS8vII8Tvj9CQx74zCa9zaz5IIl1D6n4ltE+kaF90Awfnzo5/2nP8HYsbBgAcyeDZ//PDTqCm0iIl3JK8rjxD+cyLA3RMX3+UvY96T6fItI76nwHkguvTS0fn/4w6G7yTe+Aa94Bfzzn+mOTEQko+UV5XHinScy/E1Rt5MLlrD7Qc12IiK9k5DC28wuTtB+zjOzFWa218xujllfYma3mNk3zOzHZlbZk23Sgaqq0M/78cfhxBNh9Wq44AJ473th1650RyciKaCc3Td5xXmccOcJjH7vaFoaWlj6hqVs//32dIclIlmkV4W3mV1kZo+b2WozWxst64G7+huImVUBZwNnAO8CrjGzC6LNPwIed/fPAfcAP4l5aVfbpDNnnQXPPgs33ADFxfCrX8H06TB/vgZfiuQI5ezEyyvIY9pt05jwqQl4s/PSFS+x8TsbcV2wTER6wHqTLMxsJ/BN4DmgpXU1cIW7z+tXIGYlwAGPAjKzZ4BPASuAdcAQd28wswJgH3AicKizbe6+Pnb/1dXVvmjRov6EmLtWroSrrw6zoEC48uWtt4bLz4tI2pnZM+5e3YfXKWcn0cYbN7L22rUAjPnAGI7/4fHkFakHp8hA11XOLujlvta6+3c7OMDyPkUWw92PjPIzszLgJeBh4B3ALndviJ7XbGbrgPOApi62re9vTAPG1Kmhn/edd8InPgGLFsEZZ4SrYH7zmzByZLojFJG+Uc5OoomfmcigYwax/N3L2fqzrTSuaeTEu06kcEhhukMTkQzV2z/NP2Rmb+lg/SWJCAbAzM4BHgAKgVJgHBA/gqUWGNvNtnZ27txJdXX1kWX+/PmJCjk3mMHb3gbLl8N110FhIdx2WyjKb7kFmpvTHaHIgDJ//vwj+QoY3sfdKGcn2ci3jOSUR0+haHQRNQ/V8Ozpz1K3tC7dYYlIivU0Z/e2q8lHga8BFbGrAXf3/D7GGn+MMcBrgRuBPwBrgbe6++kxz3kW+F30sMNt7v6d2P1mw8+WGWXlSvj4x+HBB8PjmTPh29+GV786FOkikjL96GqinJ0iBzYe4IVLXqB+ST15pXlM++k0Rl0xKt1hiUgadJWze9vifT0wDzgeOCZm+Up/Aozl7lvd/TZCX8HzgM3AkLinVQBbutkm/TF1Ktx/f7ji5eTJ4YI7F10UZkDJoi9DkQFOOTtFBk0cxOwnZjPqylG0NLTw0jtfYtXHVtFyqKX7F4vIgNHbwvsB4B53X+vuG6JlPfDTxIfGIkKSfggYGw3kwcwKgYnAI91sk/4ygze+EV56KbR2DxkC//oXnHoqXH45rFmT7ghFpGvK2SmUX5rP9F9O5/gfH48VGptv2czi8xbTuE4XKhORoLeF9++B68zs3JjlPOC6/gZiZqVmNitm1cXA9919G3Af8Kpo/QXAXe6+uatt/Y1HYgwaBJ/6VCi0P/OZMP3gHXeE6Qc/+lHYsSPdEYpIx5SzU8zMGHf1OGY9NoviCcXsX7CfRacsYvvvNN+3iPS+j/fThJ8sa4DWFxow2t0H9SsQs9OA/wOWAwuAhe5+Z7RtCKH/4HpgNHCdu9d1ty1WtvUXzGgvvwxf/jL88pfgDuXl8OlPhxlRysvTHZ1IzulHH2/l7DRq2tPEiqtWsOtP4eJko949iuNvPZ6Cyt5OKCYi2aSrnN3bwvttwAPuXhu3/mJ3v79/YSZXLiTxjLN0aZgB5b77wuORI0NBftVVYVYUEUmIfhTeytlp5u5s/flWVn98NS0NLQw6ZhDTfzGdwecOTndoIpIkiRxc+XJ8Ao+M7n1YkvVOOgn+8hd45JFwsZ0dO+Caa+CEE+D223UFTJH0U85OMzNj7AfGMueZOZSfUs6BtQdYfN5iVn1sFYfrlSNFBpreFt4fjl9hZm8Fbk5MOJKVzj0XnnwS/vjHMBvK6tVwxRVhCsLf/14FuEj6KGdniLLpZcx+ajaTvjQJKwgDL59+xdPUPFKT7tBEJIV6W3gPNbMvwpGBNT8DfgY8nfDIJLuYwZvfHLqfzJ8PkyaF2VDe8Q4V4CLpo5ydQfKK8pjy1SnMXjibspPLQuv33MWs/PBKmmqa0h2eiKRArwpvd38d8HczuwVYDEwFTgY6ujKaDESFhaGP98qVKsBF0kw5OzNVzKpgzsI5TP7KZKzA2PLjLSycvpDtt2+nN+OuRCT7dFl4x01Bda6ZnQsUARuAVcCXgUnAl5IfqmSVoiIV4CIpppydPfKK8pj85cnMeXYOlWdX0rS9iZeueIklFy6hYWVDusMTkSTpclYTM1sAzCBMRdXZ5bcMGOPuxYkPL3FyZYR81jp0CH71K/j612HDhrBuxgz4/Ofh7W+HAk2vJdKZns5qopydnbzF2faLbaz5zBqa9zRjRcbEz0xk4mcnkl+Wn+7wRKSX+jOryfXA8e4+yd2ndLJMBt6Z6KAlx3TWAv6ud8G0aWHdwYPpjlIk2ylnZyHLM8a8fwynrTiN0f8xGj/kbLhhA09Ne4ptv9mGt6j7iUiu6LLwdvf73b3byxK6+x8TF5LktNgC/Gc/g+OOg7Vr4YMfhGOPhZtvhvr6dEcpkpWUs7Nb0fAipt82nVmPz6J8TjmHNh9i+buX8+xZz7Jvwb50hyciCdDbWU1EEqOoCN7//tDqffvtod/35s3wX/8F48fDtdfCxo3pjlJEJOWqzq5izsI5TP/ldIrGFFH7VC3Pnfkcy965jMZ1jekOT0T6QYW3pFdBAVx+OSxeDPfeC2eeCTU1cOONcMwxof/3k0+mO0oRkZSyPGP0e0Zz2srTmPj5iVixseN3O1g4bSGrPraKQ9sPpTtEEekDFd6SGfLy4A1vgCeegAULwuwnZvCHP8BZZ8EZZ4SZUJo0162IDBwF5QUcc8MxnL7idEZdOQpvdjbfspkFxy5g3ZfW0byvOd0hikgv9LvwNrNZZvYbM/ukmU3uf0gy4J1+Ovzud7BuHVx3HQwdCk89FYrxKVPgW9+C3bvTHaVIVlLOzk6DJg1ixq9mUP18NcMuGUZLfQsbvraBBccuYON3Nury8yJZIhEt3lcDF7j7TcBZZlaSgH2KhL7e3/gGvPwy/M//wPTpoR/4ddfBhAlw9dXw4ovpjlIk2yhnZ7Hyk8qZec9MZv17FlXnVNG8u5m1n17LgikL2PjtjTTXqQVcJJMlovD+CDAnuv97oDAB+xRpU1oaZj158UV48EG46CJobISf/AROOgnOOQd+8xs4cCDdkYpkA+XsHFB1VhWnPHIKM++fScWpFTTtbGLtZ9by1JSn2PCtDTTXqgAXyUS9LrzN7LVm9s7o/mmEOWM3A7h7i7vvT3CMIkFeHrzmNfDAA7BsWWjxLi+Hxx+Hd78bxo2DT3wCnn8+3ZGKZAzl7NxlZgx77TBmPzWbmQ/MpPKMSpp2NbHuunUsmLyA9Tesp2mvxsWIZJJeFd5m9nXgTuBSAHdfCFxpZlckITaRzs2YAT/+MWzZElq+Z8+GPXvge9+Dk0+GV7wizIyyaVO6IxVJG+XsgcHMGHbRMGY9MYtX/O0VVJ5dSfOeZtZ/cT1PTniS1f+1mgMb9IugSCbobYv3WcB44LGYdQ8A30pYRCK9UVEB8+bBM8/A00+HVvAhQ+CFF8Jc4BMnwqteBb/4BexXw54MOMrZA4iZMfTCocx6bBYn//Nkhlw4hJb6FjbdvIkFxy5g2TuXUbu4Nt1higxovS28/+7uNUDs9WvfksB4RPquujq0gm/bBnffDW95S7hQz7/+Be97H4wcCW96E/zv/4a5wkVyn3L2AGRmDDl/CCf/7WTmPDuHkVeMBGDH73bwzKxnWPLqJez56x5dil4kDcy95//xzOxy4BhgHKEF5a2EnzA/6u4/TEqECVJdXe2LFi1KdxiSajU1cNddodh+9FFo/fdeWAgXXghvfjNccgkMH57eOEW6YGbPuHt1H16nnC0AHNhwgE03b2LLT7fQUt8CQMnUEsZdM47R7xlNQVVBmiMUyR1d5exeFd7Rzs4A3gtMBLYCt7v7PxIQ5GXATUAl8Bvgk+7ebGbDga8DNYTR99e6e1P0mk63xVMSF7ZsgT//ORTijzwCLeHLB7PQWn7RRWHw5umnhytqimSIvhbe0WuTkrOjfSctbytnJ0fT3ia2/M8WtvxoCwc3HQQgryyP0VeOZtw14yg7sSzNEYpkv4QW3h3sPB+Y6+7/7Mc+JgI3AN8HTgR+DHzZ3b9jZg8Bn3X3p8zsamCSu18Xva7TbfGUxKWdHTtCd5S77oKHH4ZDMZdfrqgIxfdZZ4VL2J9xBgwenLZQRfpTeHewr37n7Gg/Sc3bytnJ1dLcwu57d7P51s3UPNTW9W7w3MGMvWYswy8ZTl6RLm4t0hd9LrzN7F6gvJv9Dwd2ufv5/QjwHOBJd2+OHt9ISORfA/7k7mOj9eOAFcAoYGZn29y9Pv4YSuLSqfr60AL+4INhWbXq6OfMmAFz5oQZU045JdyOGJH6WGVA6mnhnaqcHR0rqXlbOTt16l+sZ/OPNrPtV9uOdEMpHF7IqHePYvT7RlN+Unf/pEQkVlc5u7vf05+PnrMCaAGuAp4DYrPhFGBofwJ098fiVm0GKoDzgXUxz9tsZgDVwNldbHukP/HIAFNWBhdfHBYIXVKefDIsTzwRZkx56aWw/OY3ba8bM6Z9IX7yyTB1KuTnp+c8RFKUs0F5O5eUnVjG1B9O5ZhvHMO2X29j6/yt1C+tZ9P3NrHpe5uoOLWC0e8bzcjLR1I4WNdbEumP7grvbxNaxWsAzOxYd/9S/JPM7J4Ex3VqdOx5wO64bbXAWMJgoc62HWXnzp1UV7f98TFv3jzmzZuXqHgll4wdGwZdvvnN4fHBg7B4cViWLAnL88/D1q1hefDBttcOGhQubX/CCe2XY49Vv3Hplfnz5zN//vzWhz0d/ZuunA0JztvK2alXUFXA+I+OZ9xHxlH7TC3bbtvG9t9tp/bpWmqfrmXNf61h2BuGMfLykQy9eCj5g9TIINKqpzm7y0rA3ffFrRplZvnufrh1hZnNBU7rY5xHMbPjgJ3uvsTMHGiMe0ox0ESYHquzbUcZMWIE+tlS+qS4OPT5Pv30tnUtLbBuXftifMkS2LChrUiPVVgYWsNbC/EZM8LtccdBSUlqz0eyQmyhaWa7evKadOTsaJ8Jz9vK2eljZlRWV1JZXcmxNx3Lrj/vYuttW6n5Zw0779zJzjt3kl+Rz/BLhzPy8pEMuWAIeYXqDy4DW09zdm+b4P4ILDGzxwmtFDOA1wBHtaj0RTTo54PAZ6NVm4GpMdsNKAO2dLNNJLny8kIr9rHHtrWMA+zbF7qkLFvWdrtsGaxfDy++GJZ448aF/Rx3XPvbY4/VoE7pr6TmbFDeznX5JfmMumIUo64YxYGXD7Djjh3s+P0O6p6pY/uvt7P919spGFbAiLeMYMRlIxg8d7AGZYp0oS/TCY4FPgy8AtgP3OnuCfnZ0syuBX7p7tujx9XAXe4+KXo8hdCHcVh0/A63ufuh+H1roI6kVX09LF/evhhftiy0mjc3d/66YcNg8mSYNKn9bev9qqrUxC9p1c/pBJOWs6P9JyVvK2dntoaVDaEIv30HDS81HFmfX5HP0NcOZfgbhzP0tUMpHKI+4TLw9Gdw5VHcfQvwhZid55nZHHd/ph8xYmafIwwAKjGzY4C5wKPAVjM7yd2XAhcBt0QJepGZdbZNJLOUlYVZUebMab++uRlefhnWrIHVq9vfrlkDu3eH5ZlO/nsNHhwK8IkTQ8t56zJ+fNv9ysowV7kMSMnK2dG+lLcHqNKppUz+4mQmfWES9S/Us/POney6Z1e4/4ed7PzDTqzAqDq3iuFvHM6wS4ZRMlnd6kR6e+XK2zpYPQI47O5v6nMQZl8gTEEVa7m7zzCzSYSfRdcTpsn6fMz0VZ1ui6fWE8k67rB9e+g3vn59WOLvNzR0vQ8IRX9sUT5uXBhAGruMGRP6sktG6seVK5OSs6N9JzVvK2dnp8a1jez+y2523bOLmkdr4HDbtpJpJQx99VCGvHoIg+cOpqBcA84lNyXsAjpm9hjwT8IAmVZnA4+5e3wCzihK4pJz3GHXrlCEv/wybN7c8VJ/1LT2HRs69OhiPL5AHz0aioqSelpytH4U3srZkjZNe5rY88Aedt2ziz1/3cPh/W1VuBUalWdVhkL8VUMon12uAZqSMxJZeE9z9xVx66YD73X3z3bysoygJC4DknsY8BlfjG/dGuYrj73tqq95rOHDuy7Ox48PBXqevkQTpR+Ft3K2ZISWphZqF9ay52972Pu3vexfuD/MNB/JK8uj6uwqBp87mKrzqqg8tZK8YuUQyU7JvmT8FOBpd+/pPLNpoSQu0oWWltB63lqIxy6x67Ztg8OHu99fUVHoez5lShgIGn87cqT6nfdCgi8Zr5wtade0t4maf9Ww5697qHmkhsaV7WeZzBuUR+UZlVSdV8Xg8wZTeXol+aWaN1yyQ8IGV5rZOtr/ZJlPuAzw//Y9PBFJu7y8UAyPHBmuwNmZw4dDgd5ZYb5lS+j2smMHrFoVlo6UlYV5zGfObL+MHJmc8xuglLMlUxUOKWTEm0cw4s0jADi49SD7Ht1HzaM11DxSQ8OLDdQ8XEPNwzVsYANWaJTPKqfyzEqqzqyi8sxKiicUY/oDXrJMb7uazAd+R1sibwG2u/vKJMSWUGo9EUmh+vq2QaDr1rXdtt7fu7fj140cCaecAqedFpZTTw3dVga4fnQ1Uc6WrHRo5yH2Pb6Pmkdq2PfIPuqW1LX/ExIoGlvUrhAvn12uq2lKRkhkH+9z3f3RuHXTCCPkV/cvzORSEhfJILt3w9Kl8MILbcvSpVBbe/RzJ0xoK8RPOy1MyVhRkfqY06gfhbdytuSE5n3N7F+4n/1PRsuC/TTXtB+XYkWhVby1EK88s5Li8WoVl9RL5DzeFxHmaI21HXgImNWH2ERkIBo2DM47Lyyt3GHjRnj2WVi4MCxPPx26rrz8Mtx1V3ieWeimctppcO65cOmlupBQ55SzJScUVBUw9MKhDL1wKADe4jSsaGD/grZivP7FemqfqqX2qVq4ObyuaHQRFadWUFFdEW5PraBouGZmkvTpUYu3mX0EuAHorJnp7+5+USIDSzS1nohkoZYWWLGirRBfuBCWLIGmprbnDBoEl1wCb3976JoyfnzODdzsbYu3crYMRO1axRdEreJ7j56tqXhSMZWnVh4pxCtmV1BQpTnFJXES0tXEzF4BvAf4QdymQ+6+tX8hJp+SuEiOOHAgFN8LFsDdd8PDD7ffPmRIGKh57LFhFpXWZfLkMP1hFk5z2JeuJsrZMtC5O41rGql9ujYsi2qpfaaWloaWo55bMq0kFONRy3j5yeXkl6m/uPRNIvt4F7p7UwfrJ7r7xn7EmHRK4iI5auNG+O1v4R//CAX57t2dP7e4uG2aw/jlmGPCRYQyUD/6eCtni8Tww079S/XtivG6JXX4obhayKBkagnlp5RTMauC8lPKKT+lnKJR6qYi3etz4W1mnwLud/dl0ePPE6ajavc04P+5+9zEhJscSuIiA4B7mNLwxRfbz6LSen/nzq5fP3ZsmFUldjn22LS3kve08FbOFum9loMt1L1Q19YqvqiWhmUNePPR9VHRmKIjRXj5KeWUzyqn5NgSLC+3urdJ//RncOUUYHDM49lAGbCNtol9DNB8XyKSfmYwblxYOlJXd/TUhq3316xpm4v8/vvbXlNeDmecAa98JZx9drhfXp6Ks+kL5WyRXsorzqOyupLK6soj6w4fOEzDsgbqFteF5bk66pbUcWjrIfZs3cOeB/a0vb4sj7KTyiifWU7ZSWWUzQxL0Qi1jsvRetvVZDaw1N0Pxa0/obWFJVOp9UREutTSEgrwxYvhuefC7eLFsHlz++fl54eLDLUW4mef3XmhnyD96GqinC2SIN7iHFh3gLrFddQ+V3ukKD+0+VCHzy8cWUjZzKggn1kWivITy9R3fABIZB/v97j7r+LWTQOmuvtf+hdmcimJi0ifbNsG//53WB5/PEx3ePhw++dMntxWiM+dC9OnJzSEfhTeytkiSXZo1yHql9ZT/0LMsrSew3WHj36ywaBjBrW1kM8MxXjJ8SXkFWXfwG/pWL8LbzM7h/Dz5DxgftzmocAP3H1ifwNNJiVxEUmI+vowreHjj4di/Iknjr7wzwknwHveA+9/f5izvJ/6MJ2gcrZIGrk7BzYcOFKEtxbkDcs77jtOPpQcV0LZjDJKTyildEZpuD+9VC3kWSgRhfdM4NeE/oPx13o+BPzW3a/vb6DJpCQuIklx+HC46ubjj4flr3+FvVGaHDQIrrgCrrkGZs/u8yH6UHgrZ4tkoJZDLTSsbGjfOr6sngN9VcTrAAAe6klEQVTrDrSNwohTPKk4FOEzSik9ofTI/cKhhakNXnosUfN4VwFvdfefJTK4VFESF5GUaGqCBx+EH/8YHnigbf1ZZ8H3vw/Vve4x0td5vJWzRbLE4cbDNK5spH5ZPQ0vNdDwUgP1y+ppXNWIN3VcpxWOLKTshNAqXjKthNKppZROK6V4UjF5Beq2kk4J6+PdxQFe6e6P93tHSaQkLiIpt2oV/Pd/w89/Hh4PHw6bNoX5xHuhr328u9ifcrZIFmhpauHA2gPUv1RPw7KoIH8pFOcdXQgIwIqMkmNLKJ1WSsnU9reFwwuxHLuybybq83SCZnYv0N28WQacCIzsW3giIlnu8GHYvj1czGf16lBwr1wZBmKuWtX2vMmToSh5U4wpZ4vklrzCPEqnhZZs3tS23lucgy8fDEX48gYaVzbSsCLcHtx08EirebyCwQXtWsePFObHlZBfqr7kqdDdPN4rgBpgM532PsKA1ycyKBGRjNDcDLt2hWXnzjC1YOyyaVO43bbt6JlOWhUWwpw5cNllcNVVYa7x5FHOFhkALM8YNGkQgyYNYthF7QdwN9c107i6kcYVjTSsbDhSkDesaKC5ppnap2qpfar2qH0WTyymdGrpUYX5oImDsHy1kidKd4X3fwO17n6wqyeZWcKmpTKzC4EbgLe7+/poXQlwI1ALDAGudff93W0TEQFCAb1/P+zbF5aamrbb1qI6/nbnzrC9p0aOhPHjw5Uujz8+LDNnwkkn9bprST8oZ4sMcAXlBVScUkHFKRXt1rs7h7Yfatc63rCigYaVDRxYc4CDGw9ycONB9v6j/XhsKzZKj4/rtjK9lNLppRQO1gDP3uqy8Hb3XT3cz6VAvwfwmNkows+kp8Vt+hHwoLvfYWYXAT8B3tGDbSKSzZqawlR9rcv+/e0fx65vLaZjl9Z19fV9O75ZmA5wxIjQP3vs2HCxnPHj266QOW4cjBmTyuK6U8rZItIZM6N4dDHFo4sZfO7gdttamlo4sP5Au4K89fbQ1mie8qVH59Gi0UVHivDYpXhCMZanVvKO9GtwpZlNBK4BPuLuZQkJyCwPOAxMcff1ZjYWWAcMcfcGMysA9hH6KB7qbFtry0srDdQRSZHDh9sK4dbCN75g7qyAjl8Odtlw23NmUFkJVVUweHC4bV2GDw+FdWtxHXs7ZEi4UmWaJWpwpXK2iPRWc21zKMJbu62saKRhebjf0tjxAM+80qhvelxBXnJ8Cfkl6c+pydbnwZWd7CwfuAS4CrgQ2E7oU5gQ7t4SN+J2LrDL3Rui7c1mtg44D2jqYtv6RMUkMqAcOnR0d4zYIjr2cUfr4i8m0x/5+VBR0f0SW0i3LrEFdnk55A3M6bWUs0WkPwoqCqiYU0HFnLiuK9EAz4blDe2W+pfqadreRN1zddQ9V9d+ZwaDpgxqK8ZnlFJ2UhllJ5RRUNnrkjQr9fgszWwK8AHgP4BhwAbgDcBfgTOTEl0wDtgdt64WGAs0d7GtnZ07d1IdM3/uvHnzmDdvXmIjFUkH91As19dDQ0Pbbfz92tquC+bWxw1Hj4TvtdbCt7X4razsvnju6DmDBiV7MGLGmj9/PvPnH7no5PDevl45W0SSKXaA59DXDG23rWlvU+g/HluUv9RA45pGDqw9wIG1B9hz/552rymeWByK8BPLjtyWzijNmtlWepqzuy28zeythJaS84HVwHeBXwGfc/cHo6c90a9ou+ZAY9y6YkLLCd1sO2LEiBHoZ0tJi6amjgvh7u735rktHf/c1yd5eW1Fc2zx3NN1FRUZ0T0j28UWmmbW077bytkiknaFQwqpOqOKqjOq2q1vOdRC45rGI4V4/bLQd7xhecORwZ3tCnKDQccMOrogn1ZKXnFm/YrZ05zdkxbvRkK/vCeAD7j7ymin/b/yTs9sJox8j1UBbCEk+M62iXTNPfQhbi1iOypsu1p6+rymo2qKxCsshNJSKCtrfxt7v7y8Z8VzWdmAbWXOEcrZIpKR8oryKJtRRtmMsjDEO9LS3MKBNQeof7H+yEDO+hfraVzZyIE1Bziw5gC772n7scwKjLKTyiifXU7FnArKZ5dTfnJ5VvQf77bwdvf7gPuiQTlXmdlk4B7gyBwyZmaeiEtgduwh4KdmVuLujWZWCEwEHiEM6Olsm2SzlpaeFbX9LZST9s82Rl5ex4Vwd4VyT++XlobCWwTlbBHJPnkFbRcKGnHZiCPrWw610LCyIbSKv9hwpChvXNNI3eI66hbXse22beHJ+VA2IxTjlWdUUvXKKspOLMu42VV63Mfb3TcCX4xGob8JGG9mXwN+DrwamN/V63vK2kbpWHTcbWZ2H/Aq4D7gAuAud98cPb/TbZIksV0nktVSfOBAas6lqKh9ARtf0Ha29OZ5hYVqQZaUU84WkWyXV5RH+UnllJ/U/oK8zXXNofB+to7aZ2upe6buSLeV+qX1bP/1dgDyq/KpOquKqldWUXVuFZVnVJJXkN4uKv2dTvA44KPAVe5e2u9gzCqAdxHmeb0euNXdd5rZEMIFF9YDo4Hr3L0uek2n22INyKmpDh+Gurqup2uLn9ato6I4fl1zc2ri70lh29MiuKPnlJRAwcAYRS3ZLYHTCSpni0hOOtxwmPoX6qldVMu+J/ax7/F9HNzYfkra/Kp8hl44lKEXDWXY64dRNKooKbF0lbP7VXjHHOBKd/91v3eURFmXxA8dgr17j1727Am3XV1IpPV+Iman6Eh814lktBYPGjRgp38TiZeowjtmf8rZIpLzDrx8gH3/3se+x/ax9597aVwRM7Y7H4ZcMIRR7xrFiEtHkF+WuP7hCZ3HuyOZnsDTrqkpXH56+/aw7NjRdn/37qOL6r17E1c0l5f3bgq3srKjC+OO+hOr64RI1lLOFpGBYNCEQQy6fBCjLh8FQOO6RvY8uIfd/7ebvX/by96/hmX10NWM/9h4xn10HIVDkztmSr+z91dLC2zdCuvXw4YN7W83bWorrnsrPz9cNa91GTq0/ePWadu6WsrK1GosIiIiApRMKWHch8Yx7kPjaNrdxI47d7DtF9uoXVjL+q+sZ9MPNnHczccx6l2jsCQ1MKrw7o2WFliyBP79b3jhBXj+eVi6NPSj7kpeXrj89MiRMGpU+2X48PYFdetSXq5WZREREZEkKBxWyLirxzH2g2PZ9+g+1n9lPTUP17D8yuXULqzluO8fl5QZUVR4d8cdnngCfvITuP/+jluvhw+HKVNg8mSYNCncTp4MEybA6NEwbJguKCIiIiKSYcyMwecN5uR/ncy2X2xj5YdWsvnWzRSNLWLSdZMSfjwV3l3ZtQve855QcLeaMAHOPx9mzYKZM8MyYkTn+xARERGRjGZmjHnfGApHFrL0DUtZ/9X1jHn/GIpGJnbmExXenWlpgUsugSefDFfz+/CH4corYepUdQERERERyUHDXz+cYW8Yxu6/7GbHH3Yw/iPjE7p/Fd6d+ec/Q9E9ZgwsXAjjE/vGi4iIiEjmGXrxUHb/ZTe1C2sTvm9NedGZp58Ot5dfrqJbREREZIAoOaYEgENbDyV83yq8O7N3b7gdPTq9cYiIiIhIyhQMDR1CmvY0JXzfKrw70xhd3WjQoPTGISIiIiIpk1ccymNv6v/V3Y/ad8L3mCtaWsJtgbrBi4iIiAwUVhAm0VDhLSIiIiKSREcK72YV3iIiIiIiSdNaeLc0tSR83yq8RUREREQiavEWEREREUkBy48ulHg48ftW4S0iIiIiEmktvP2wWrxFRERERJInP9yo8BYRERERSSK1eKeDJ/7NFhEREZHMlsw+3jlxdRgzKwFuBGqBIcC17r4/QTtPyG5ERCRIas4WEekntXh370fA4+7+OeAe4CdpjkdERDqnnC0imUt9vDtnZmOBK4C/RKv+AVxiZpPTFZOIiHRMOVtEMp1avLs2F9jl7g0A7t4MrAPOS2dQIiLSobkoZ4tIBlMf766NA3bHrasFxsau2LlzJ9XV1Ucez5s3j3nz5iU/OhGRPpg/fz7z589vfTg8nbEkmHK2iGQ0y2sb3+ct3u5xZ3qas3Oh8HagMW5dMdAUu2LEiBEsWrQoZUGJiPRHbKFpZrvSHE4iKWeLSObLA1p6Xnj3NGfnQleTzYRR8bEqgC1piEVERLqmnC0iGS9Z3U1yofB+CBgbTU+FmRUCE4FH0hqViIh0RDlbRDJesgZYZn3h7e7bgPuAV0WrLgDucvfN6YtKREQ6opwtIlkhSVMK5kIfb4APATea2cnAaEAjcEREMpdytohktGS1eOdE4e3ue4Gr0h2HiIh0TzlbRDKd+niLiIiIiKSA+niLiIiIiKRCkvp4q/AWEREREYmhFm8RERERkRRQH28RERERkRRQi7eIiIiISCpEFbK3qPAWEREREUkadTUREREREUkBdTUREREREUkFTScoIiIiIpJ8avEWEREREUkB9fEWEREREUkBtXiLiIiIiKSC+niLiIiIiCSfWrxFRERERFJAfbxFRERERFJALd4iIiIiIqmgS8aLiIiIiCSfupqIiIiIiKRAsrqaFCR0b31kZga8Ffiqu8+I2zYc+DpQAxQC17p7U3fbREQkeZS3RSSn5fh0gpOAIcD0DrbdCdzm7tcCK4Hre7hNRESSR3lbRHJWTg+udPf1wN/j15vZGcA0d38qWvUX4KNmVtbVtlTELCIykClvi0guGwh9vFs6WHc+sK71gbtvju5Wd7NNRESST3lbRHJSTvfx7sI4YHfculpgbDfbjrJz506qq9ty+7x585g3b17nR66qgjFjoLS091GLiPTT/PnzmT9/fuvD4emMpZcSkrd7nbNFRBKoYFgBRWOKyCvuWRt1T3O2uSe2ku/wIGbfBGZ2svkf7n6zmU0G1rm7xbzuVmCEu789Zt0eYB4wt7Nt7v7H+INUV1f7okWLEnA2IiKpZWbPuHtKW4XTnbeVs0UkW3WVs1PS4u3u1/XxpZuBqa0PolH0ZcCWbraJiEg/KG+LiCReJvXx7sjfgWkxjycDh4BF3WwTEZH0UN4WEelEJhXeBkdaQABw90XAVjM7KVp1EXCLux/qalsqgxYRGcCUt0VEeiEjBlea2UjgPdHDD5nZr929Lnr8duBLZrYeKAc+H/PSrraJiEiSKG+LiPReRrR4u/sOd/+Ku5u7/ygmeePuG9z9/e7+NXe/1t2be7ItEWJGp+aEXDofnUtm0rkMHMrbmUHnm9sG2vlC7p9zRhTemSrXPvxcOh+dS2bSuUi6DbTPTeeb2wba+ULun7MKbxERERGRFEjJPN6ZwMx2Aht6+bLhwK4khJMuuXQ+OpfMpHNJjknuPiLdQaRSH3M2ZNbnlgo639w20M4XcuOcO83ZA6bwFhERERFJJ3U1ERERERFJARXeIiIiIiIpoMJbRERERCQFMuICOpnGzEqAG4FaYAhwrbvvT29URzOzC4EbgLe7+/poXaex93VbCs7jMuAmoBL4DfBJd282s+HA14EaoDCKqSl6TZ+2peBczgPmAyOBX7n7f0brs+5ziTmnIuBp4OPu/nCWn8sNtF2w5Xl3Pzmbz0faDJTPo7d5P1v15Xsh2/Xl+yMX9OY7Jie4u5a4BfgFIalBuKTx7emOqYMYRwGXAg5M7knsfd2W5POYCPwamANcCdQDn4q2PQScHt2/GvhmzOv6tC3J51IFfI6QKF4HNAEXZOPnEndenwf2AXOz+VyAMuB/gAui5fhsPh8tR32+Of959CXvZ+PS1++FbF76+v2RC0tvvmNyYUl7AJm2AGOBg0Bp9Lgg+k8/Od2xdRBrXmwC7ir2vm5LwTmcAxTEPL4R+D/gDGBLzPpxQF1UPPVpWwrOpYRopqDo8TPA/8vGzyXmHM4C3gesB+Zm+bn8Z7QMilmXteejpd1nO2A+D3qR99Mdaz/OsdffC+mOOQHn3Ovvj3THnKDz7vF3TLpjTdSiPt5HmwvscvcGAA+XM14HnJfOoDri7i1xq+bSeex93Zbsc3jM218yejOwETg/iqH1eZuju9X92JZU7t7oUaYwszLgJeBhsvBziTmHt7j7bTGrs/VcDLic8NP1djN7RzafjxxlLgPk8+hl3s9KffxeyGp9/P7Ian34jskJKryPNg7YHbeulvBXWKbrKva+bku1UwndAbL2XMzsHOABQv/D0m5iyuRzuRb4Vty6rDwXD84gXJjhFuC3ZvbabuLK2PORowzkz2MgnHtPvhdyQi+/P7Jdb79jcoIK76M50Bi3rpjQ3yrTdRV7X7eljJkdB+x09yXdxJTp57Ia+CXwKuDb3cSUkediZq8BFrn7jrhNWXcusdx9r7t/AfgmodtJVp+PHDGQP4+cPvdefC/kit58f2StPn7H5ATNanK0zYTBDbEqgC1piKW3uord+7gtJcwsH/gg8Nlo1WZgasx2I/Th3tKPbSnh7luB28ysBfg08BjZ97l8CpgV3j6I4rgH+G4XMWXquXTkB4SfcbP2/4y0k815u79y9tx7+b2QE3r5/ZHN+vIdkxvS3ck80xZgNGGwRkn0uJDQ0X9cumPrJN7YQTadxt7XbSk8j2uBUTGPq4ENMY+nEH5uKurrtjR8NicBf8vGzwUYAYyPWV4G3kqYbSCrzqWL8/ttNn42Wjr8PAfU50EP836640zAefb4eyHdsSbh3Lv9/kh3jP08v15/x6Q75kQt6moSx923AfcRfuaBMO3YXd42iCNjWNufigZdx97XbSk4Dczsc8BzQImZHWNm7yPM0brVzE6KnnYRcIu7H3L3RX3ZloLzKDWzWTGrLga+n42fi7vvdPdNrQtwmPBz78ZsOxcAMxtrZleYWV70/+aTwOez8bORow2kz6M3eT8N4SVMb78X0hVnovTl+yPVMSZSX75j0hVroln0F4XEMLMhhOmL1hP+2rzO3evSGlQcM6sA3gX8CLgeuNXdd3YVe1+3Jfk8vgB8LW71cnefYWaTgC9FMZUTCqXm6HV92pbkczmNMOXVcmABsNDd74y2ZdXnEs/M1gPv9XBxg6w7FzM7kTBgaT/wOHCzuy/vT8yZ8tlIMBA+j77k/WzU1++FbNbX749c0dPvmFygwltEREREJAXU1UREREREJAVUeIuIiIiIpIAKbxERERGRFFDhLSIiIiKSAiq8RURERERSQIW3iIiIiEgKqPAWEREREUkBFd4ikjXMbISZ3WtmG8zss+mOR0REuqfc3UYX0BGRrGFmrwEeAioJV3ib6u570huViIh0Rbm7jQpvEclKZvYocH4uXC5aRGSgGOi5W11NRCTrmNkE4G8DNXGLiGQj5W4V3pJiZvYZM/thuuNolYh4zOwCM3vSzOYmKKycYmbfSfD+DHgb8K1E7ldEjqac3X9mNtnMfmpmP0nQ/vLN7Coz29CD577OzJ7q67kmMn8rdwcqvHOImc0wsz+Z2c/MbImZuZm9K91xxbkfuD3dQcRIRDwvAycnIJaMZGYXm9kuM9tqZmea2YTo35mb2a1mNjh63mlmttTMbjKzkmjdKGBbgkO6Avi5uzdH+xfJSsrZfZKNOXtHdFucoP3lATXAxB48dxEwo7ONZnaimd0R/dv7gpmNjtmW6Pyt3I0K71xzN/BVd/8AMAv4desGM3uFmZ2XjqBij+3uS9398XTE0ZHexNPZe+juK4BdiYwrkZ9Xf/fl7vcDPwKa3f1Jd38ZuAo4BDzl7jXR8xYSkvyn3b0xevkbgPv6dQIxzOyLwFeAf5jZcmBaovYtkgbK2b2UjTnb3RuALYk6lrs3Ac/08LnbCUV6Z9tfBOZHD7/j7rGFdsLyt3J3GxXeOcLMhgJTgXoAd28BvgS4mVURErqlIa60HTuRenAeLSk8Vjr2dQcw3szmALj7buCfwOtjjjUY2BL922s13d2X9/PYR7j719z9eHevdvfp7v5oovYtkkrK2cmVgTk70TNZ9Cb+7p57OLqN73edsPyt3N1GhXeOiKblWQE8aGZnROs2AAsJfaqmAPPM7D8AzOx6M/ti9BPTkT5cZlZtZt81s/80swYzW2FmnzCzN5nZ02Z2qZn92cz2m9kHY17X4f7ijn2Vmf0ytp9b9FPrLWb2TTN7yMxOi9Z3ebyY15uZfdDMXjSzC83seTPbbmbv7cExjo+Np5tjHvUedmBI9LNxrZldFfeefiV6b/4UJemevmdHHcvM3mJmC83svWa2ysy+1cX+OvrsO4ynK1GryIvAW2NWjwBeZ1G3EuAS4N6YOEuAhpjH083sW2b2HTN7xEKXlbeb2SIze6OZ3W5mNWb2cTObE53jbjO7uLv4RLKNcvbAydkd6ewcuzqOmU0xs9sstB5/uyfH6eC495vZb83s7G6ep/ydLO6uJUcW4ERgA+Gv29uAkTHb1gNzo/snAWuj+yMIf4kPI/y1vhGojrb9EHgiul8E7AFuJczD+Q5gc1f7iz82oX/bT4BfRusHEYq5wdHjtwA7gYqujtfBeY+Ojvm56Lk3AAeASd0cIz6eLo8Z+x52EMN64CagNHrdxmh9JfD7mOf9m5Awe/SedXKsoYSfDm8GTgMu6Gp/cZ99h/H08N/XF4A10f1jgW8CB4HLYv69WMzzLwHOie6XAk8DpdHjR6P3oRLYB1wPlABvjx6/jdAw8Emif4NatOTagnL2gMjZ0fav9PB97OqzfhqYEfMa7+G/s9bP80zgyx1snxsdpyBmnfJ3kha1eOcQD62SMwlJ6ErgRTPraADJKuDdZlYInBOtKyf8J58AtPbPfYnwnwl3PwTUAn909/2E/4Sju9lffHwHaT9Q43VAo0d9hIE/Ef6zXtbN8eL327rP30XPvYHQ//jCbo7RLp7eHLMTf/HQl+9pYEzMOQ6LWqP+k/CeQg/fs07Odw8hud3j7gvd/R+92F9n8fTEHcAxFrqbvBn4KfBX4G1mVg7UepSVI2cDT0T3LyF8sbW2oFwKfC16n2uAf3noF74IqHT3P3j46f05evcZiGQN5eyBkbM70Ok5dnGcC4Aqd2+NZ1Evj3kq8E53/2oPn6/8nSQF6Q5AEsPM8oEJ7r4e+JiZ/Zww+ns+cHrsc939oJmNA64FftW6C3ffYWaLgFcR/hqfQii2jrw05n4zUVelzvbXSaix+zgeKIyJq8XM1gPjuzped9z9gJmtBQYDw3txjD4fM04zbf+3JgLr3P3m+Cf15D0zs3cTCtxW0zz8HO2xsfbiM+g0nu64+yoze47Q3aTK3dea2R2EFqi3Av8XE3frv43WvoOTCK1TrfvaHbvrmPuHaa8FdYmTHKSc3e78BkLOjtXp+9jFZ/MKwh8ZfXU+cLyZlbt7XSfPaYnOQfk7ifSG5I5CYvrfuvsSws940+OfaGbnAx929xs8zFAR6xrgVAtTWj0P/Hd3B+5mf13ZSGhBjZ1iyQj9HvurCFiZ5GP0xFba94PGzE7vxXt2L3BKzNLhyPhe7K/DeHp+Ovwe+A9gbUx8+cDHCT/Jtjqd0Fe11WbglWZWFnPcLvsYiuQ45ez2BlLO7vQcuzjO/ug1RfTNdwi/Fny/k+37vG1gvPJ3Eqnwzi0fNbOxMY/HAQ9G9w8RBpJMJ0xbVWVmxWZ2QbR9pJkNJ/Qz/DnhJ6IniH62jOQR9xe+mVk3+4s/tsXs426gDrg82tfg6HitLaedHa8zY6LnjIte+2APjhEbT3fHjD2PeB21fBihBasMuM/MXmNhAE93n8GRYwFj3H15zNIUc4j8mPtd7S827g7jMbMiM5sfdSPpyh8IP2//EcDda4EHgH97+9lMXk3ohtKq9f2+3cJc4B8BWudxjf8M4j/nrJ5dQaQLytkMmJzd0/exw+MQ8mkR4Y8zCL9uYGZje5C/jdD6/G7grWb2tui1nzSzLxBa+/8Z83zl72RKRsdxLalfCIM1HNgL/Ab4JeFnqiHR9i8Dm4A3En42WkVoXXgdoZXkl4QEdi9hkMehaH8NhP+ErwWaCIN3JgDXRdvf19X+4o79ceAxQp+5E6Nt1YSW0psI/RxPitZ3erxOzt+j13+c0PVhZsy2zo4xOTae7o4Z+x7GHfuNhEGGt0bvRevr/iPafh6wNPpsbore556+Z2/s4FzfEh3vdmBctK6rz7TdvjqJZzzhIg+/68G/tR/EPX478P/i1n2ng9edCywDdgNfiNZdTPi3dgswlmg6NeBqwpfyT6NzfX26/49p0ZLIBeXsgZSzJwKPE1rtT+nmHLv6rF9H+LVxAWGw+0rgw9G5d5i/Cb+qNBFauocC/yIU/R8APkiYz/wmQjeX1tcofydxsejNE8HMRgJXu/v10eM8QgvMB939C2kNrhtm5sAUD/0lpY/M7NPu3qdpqkQktZSzJZbyd3ZQVxOJ9RFCt4OhMevOJ2bQXCaK+VlLP2n1kZnlm9k7gbvSHYuI9Jhytih/ZxkV3hLrfwj9wF4ys5eBvwOb3P3J9IbVuWigycejh++1MK2d9F45cLe7r+32mSKSKZSzBZS/s4q6moiIiIiIpIBavEVEREREUkCFt4iIiIhICqjwFhERERFJARXeIiIiIiIpoMJbRERERCQFVHiLiIiIiKSACm8RERERkRRQ4S0iIiIikgIqvEVEREREUuD/A76+N1YamOc8AAAAAElFTkSuQmCC\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/decadal-2010-uranus-probe.png',bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-uranus-probe.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-uranus-probe.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAF4CAYAAACYdLRaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdd3gU5fbA8e+hd0UJIiAgSFGs19hFEBuK7aqgotiNeu1dsRfsXnshKhZUbPwUReXaUEQsFAtFAaUjYuggJQk5vz/OrlliCGR3srPlfJ5nn5md2bxzCJCTt4uq4pxzzgWpWtgBOOecyzyeXJxzzgXOk4tzzrnAeXJxzjkXOE8uzjnnAufJxTnnXOBqhB1AKmjSpIm2adMmrq8tKCggJycn2IDcJvHvfXj8ex+eVPrejxs3bqGqlhuMJxegTZs2jB07Nq6vzc3NjftrXWL8ex8e/96HJ5W+9yIya0P3vFnMOedc4Dy5OOecC5wnlwTl5eWFHULW8u99ePx7H550+d6Lry0Gubm5miptmM45ly5EZJyq5pZ3z2suzjnnAufJxTnnXOA8uTjnnAucJxfnnHOB8+TinHMucJ5cnHPOBc6Ti3POucB5cnHOVb3iYnu5rOELVzrnNl1REcydCzNnwqxZMGcOLFkCK1ZU/Fq92r6+Th1o2BAaNbJjeectWsCOO9qrWTMQCfWP7OLjycU5V2rdOksc06ZZ8ij7+v13KCmpfLki9lqzxl4FBZv2dVtsYUlmp51KE86OO8Lmm1c+BpdUnlycy0aLFsHUqTBlyvqvX3+FwsINf50ItGwJrVpB69Z2bNKktOYR+4qtkdSrZ1+/ejUsX75+rSb2/fLlMGMGTJgAEyfC4sUwcqS9Ym27LRxxBBx9NHTtCrVrV933ysXF1xbD1xZzGaqwEH777Z8JZMoUSy4b0qIFdOhgP8Bbt17/1bIl1KyZnPhVraY0cWJpspk4ESZNstpPVMOG0KOHJZrDD4ctt0xOfK7CtcU8ueDJxaW5lSth8uTSH7y//GIJZMaMDTdh1a8PHTuu/+rQwV4NGiQ3/spatw7GjYN337XXhAml96pVg/33t0TTuzdss014cWYBTy4b4cnFpYWSEksc339f+lv8xInWR1IeEWjT5p9JpGNHaN48czrKZ86E996zRPP556Wj0qpXtwRz5ZWw++5hRpixPLlshCcXl5JWrYLvvoPRo+Grr+Drr21kVlk1a0KnTtbR3bmznXfsCNttZ6OzssmyZfC//8Gbb8Lbb1stB6xf5soroWdPq924QHhy2QhPLi4lFBZax/UHH8CoUVZDKTs3pGVL2HPP9UdObbdd8vpB0sns2fDoo5Cfb4MFwJLu5ZfDaadB3brhxpcBPLlshCcXF5qlS+HDD2HoUDsuX156r1o12GUX2G8/2HdfO7ZqFV6s6Wr5cnj2WXjkEUs4YCPcrr0WLr3UE3MCPLlshCcXl1Rz51qTzdCh8MUX69dOOne2zuju3WGvvWwklAtGcTG89RY8+CBE/79vvz08/rh9v12leXLZCE8ursotXw5DhsCgQdbpHP1/V706dOliCeXoo6Fdu1DDzAqqMHw4XHKJzesBOPFESzotWoQbW5rx5LIRnlxclSguho8+soTyzjulczNq17aO5WOPtYmAPi8jHGvXwgMPQP/+Nrmzfn245RZrKqtVK+zo0oInl43w5OICNWmStfG/+ir8+Wfp9QMOgL594YQTfPmSVDJrlnXyv/22ve/UCZ56Crp1CzWsdFBRcvHlX5wLwsqV8MYb8Mwz8M03pdc7dbKE0qePzTlxqad1a/i//7OmsosvtrlE3btDv35w661Qw39MxsNrLnjNxcVJ1TqGn3kGBg+2BAPWCX/KKXDWWZCbmzmTFbPB2rVw111w5502afWAA6wG6n0x5fKai3NBWrXKfuA8+aTNRYnabz845xzo1cva7136qV0bbrsNDjwQTj7Z5h3tuiu8/DIcdljY0aUVn6rq3KaaOtXa5ps3h3PPtcSy5ZZwxRXWzzJqFJxxhieWTNCtG/zwAxxyCCxcaAtjXn+9b3hWCZ5cnKtISQkMG2Y/ZDp2hIcftiVG9toLXnrJ5qw8+CDssEPYkbqgbbWV9cPceadNaL3nHks6c+eGHVla8OTiXHnWroWBA215laOOgk8+seVCzj7b+lm++cY66rNt7a5sU60a3HADjBhhNdavvoK9915/JWZXrpTrcxGRO4EbIm9/UtVdRKQucB+wAmgMXKuqyyOfj+uec+VasgQGDLA1qebPt2stWsBll1liadw43PhcOA44wJrJjj8evvzSJr4OHWoLYrpypVTNRUTqA02AQyKvEyK3ngRGqWo/YCgwIObL4r3nXKlly2wCXevW1rY+fz7svLM1fU2fDldd5Ykl2+Xk2KTY44+3fy+HHmrLybhypdRQZBG5LHL6tKquiVxrDswAGqvqKhGpASwDOgOF8dxT1Zmxz/WhyFls1Sp44glrT1+82K4ddBBcc431s/gwYlfWunU2i/+JJ+zfx6OPwkUXhR1VKCoaipwyNRcREeAk4EFggYicHLnVDVioqqsAVLUYSxpdE7jnsl1hoc3C3m47SySLF1tTx5dfWv/KoYd6YnHlq14dHnvMlo1RtYmX/fqVrhfngBTqc1GrQu0tIo2BK4FXRGQp0AIou+H3CqA5UBznPZfNPvsM/vMf2woY4F//solznlDcphKxhNK8uc1tuvtuW+rnmWf831BEytRcolR1iareCNwNXAYosLrMx2oDRQncW09BQQG5ubl/v/Lz8xP/g7jUs2ABnHqqNXtNmQLt21ub+dixNkHOfyi4yjrjDNtiuV49eO45uP32sCOqcvn5+X//rMT6yMuVMjWXcjwKfA7Mw0Z6xWoI/I4lkHjurScnJwfvc8lg69bZCLB+/awjtk4dG1569dU2I9u5RBx+uK0rd/TRthZZhw42uz9D5eXlkZeXB4CILNzQ51Ku5hKjBBgPjACaR4YVIyI1gVbAFwncc9li8mSbl3DhhZZYDj/cZtPfeKMnFhecnj3hv/+18zPPhK+/DjeeFJAyyUVEmotIHxGpFuncvxK4QVX/AIYBB0U+ejAwRFXnxXsvWX8mFyJVePpp2H13a/Zq0cKawN5/H9q2DTs6l4kuuQTOP98m4B5zDMycGXZEoUqZocgi0hn4EFgOjAIeVtVfIvcaY5MhZwLNgOtVdWUi92L5UOQMs2iRdbK+8469P+MMGy7qWwa7qlZUZLWYjz+2LatHj4ZGjcKOqsr4ZmEb4cklg3z+uXXaz5tn/6kHDICTTgo7KpdNli6FffeFn3+2BS/fey9j94RJi3kuziWkuNg66bt3t8Syzz7w44+eWFzybb65LXbapIktfPn882FHFApPLi79LVliTRF33WXDiW+6yfbh8J0fXVjatrWJlmD/Lov+MQMi43lyceltyhQbDfbRR7b206ef2lyDDG2GcGmkVy/b5nrmTBg0KOxoks6Ti0tfH31k+6pMnWqLTI4ZY/ttOJcKqle3Ie9gS8Vk2UZjnlxc+lGFRx6xOSvLlsG//237bLRuHXZkzq3vxBNtJYjp0+GVV8KOJqk8ubj0UlhoWwxfdpntEnnTTTZ/pUGDsCNz7p9q1Mja2osnF5c+li61oZ3PPWe7Qr7+uvWvVPN/xi6F9ekD7drBtGnw2mthR5M0/r/SpYdZs2C//Wy72WbNbGn83r3Djsq5jatRw4bJA9x/f7ixJJEnF5f6xo+3EWGTJ8MOO9j+9bvvHnZUzm26U06xDv4JE2DNmrCjSQpPLi61ffCB7V/+xx9w4IHece/SU61aNvdFFX77LexoksKTi0tdTz8NRx0Ff/0FffvabOfNNw87Kufi06GDHaOb1GU4Ty4u9RQV2RL5F1xQOiLsxRfttz/n0lU0uUydGm4cSeLTmF1qWbjQZjZ//rklk2eegdNOCzsq5xLXsaMdPbk4l2QTJthufjNn2oiwt9+2jnznMoE3izmXZEVFcM89sMcellhyc22DL08sLpNkWbOYJxcXrvHjYc894frrbQe/s86yFY1btAg7MueC1bw51K9vTb+LF4cdTZXz5OLCsXo1XHutJZYffrDl8f/3v9LZ985lGhHb4wWyIrl4n4tLvgkT4PjjbTkMEVsn7I47fH0wl9lUYcECO99qq3BjSQJPLi65hg61bYhXrrQ9xp97zpbNdy7TLVtms/MbNICGDcOOpsp5s5hLDlW4+25bHn/lSlvMb8wYTywue8yfb8ettw43jiTx5OKqnipccQX062fv774bXn7Z+1Zcdsmy5OLNYq5qlZTAf/4DAwZAzZq25Phxx4UdlXPJ58nFuYCsWwdnn21Lt9SpY5Mie/QIOyrnwuHJxbkAlJTAmWfCoEE2tv/dd6F797Cjci48nlycS5CqLToZTSzDh8P++4cdlXPhyrLkskkd+iKyuYhUi3nfVET8p4X7J1W4/HLIz7emsPfe88TiHNieRGDr5mWBjSYXEbkWmAZ8JiL9RKSaqv4J5Fd5dC69qNp2ro88Yp33b79tG3w55+DPP+3YtGm4cSTJpjSLnQ90UtVFItISuFREBgJ/VW1oLu3072/DjKtXhzfe8M5752IVFNgxJyfcOJJkU5rFfgKWAajqXFV9COgJNKrKwFyaefBB29SrWjV45RU49tiwI3IudZSUwKJFdh5dXyzDbUpyORtYr21DVV8FBlVJRC79PPssXHWVnQ8cCCeeGG48zqWaJUtsaP5mm2XNjqobbRZT1YXAx+Vcv7NKInLp5YMP4Pzz7fyJJ+D008ONx7lUFG0Sy5L+FqjC5V9E5IiqKtuliHHjoHdv+43sxhttJr5z7p+yrL8FApjnIiI9gBuBZpQmq2rAVoAvHpWpZs2Cnj3hr79sj/vbbw87IudSlyeXuAwC7ga+B0oi1wToE0DZLhWtXm3rgy1YAAcdBM88Y/uyOOfKFx2G7MmlUqar6n/LXhSRXwIo26Wa6Oz78eOhbVt4882s6aB0Lm7r1tlRNdw4kiiIPpcLROSEcq4fHUDZLtU8+aQtRFmvnk2SbNw47IicS30772zHH34IN44kCiK57Ac8KyLrYl4lwFMBlO1SyejRtiUx2A6S0f8wzrmK7bqrHSdMgMLCcGNJkiCSy+1AHtAeaBvzujWAsl2qWLoUTj4Ziott46+TTgo7IufSR8OG0KGDJZbJk8OOJimCSC4fAkNVdbqqzoq8ZgLPBFC2SwWqcN55MHs27LEH3HNP2BE5l37+9S87jh8fbhxJEkRyeQ24XkQOiHl1Ba4PoGyXCp5/3tYKa9AABg+2RSmdc5WTZckliNFiN2FNYmcA0aEQgs17uTSA8l2YpkyBiy+28yefhHbtwo3HuXQVTS7jxoUbR5IEkVzuBz5U1RWxF0Xk8ADKdmEqKoJTToFVq+zYt2/YETmXvnbbzY4//miTj+vXDzeeKhZEs9ho4EYRaQAgIvuKyLGq+mEAZbsw3Xef/ZbVqpWtG+aci98WW1jtZfVqyMvL+DkvQSSXl4G9gJoAqjoaaCEi/QIo24Vl4kS47TY7f+45W83VOZeYl16yGsurr9qmehksiOQyRlW7qeqSmGvjgMvjLVBEaonIjyLSLfK+rog8JiJ3ichTItIo5rNx3XMVKCqCM86w43nnwcEHhx2Rc5mhc2d44QU7v+oq+PzzMKOpUkEkl2IRqR19IyJ1sZFiCxIo82qgTcz7J4FRqtoPGAoMCOCe25DY5rD77gs7GucyywknwLXX2pIwvXvD3LlhR1QlRBNs9xORjtjilbOAWsD+2Gixo1V1VBzl7Qt0Am7GRqBNBWYAjVV1lYjUwHbG7AwUxnMvMg/nb7m5uTp27NjKhpqZJk2yjseiIvj4Y6+1OFcViovh8MPhk09gzz1h5EioXXvjX5diRGScquaWdy/hmouqTgG6Am8BXwGXAG3jTCz1gRNUdWDM5W7AQlVdFXleMZY0uiZwz5WnpMQ2/ioqgnPP9cTiXFWpUQNeew1at4bvvoOLLsq4Dv5KD0UWkfuB+cAwVZ0KoKqrgdcDiOdaoOz07xbAojLXVgDNgeI4762noKCA3NzS5JuXl0deXl5lY09/L7wAo0bZbnneHOZc1dpyS1v8dd99bavwRYsgPx+aNAk7sgrl5+eTn58ffbvBYCudXFT1ahFpChwlIucCf2HbII/WBNrYROQwYKyq/ln2kcDqMtdqA0WR83jv/S0nJ4esbxYrKICrr7bzhx6CzTcPNx7nssFuu9mqF6edZolm9GgYOBCOSN2NfGN/+RaRhRv6XFzNYqr6p6o+p6pXYzWNLYC7ReRuEfm3iNSLo9irgIEisjAS8DZYJ3x9oOy67g2B34F5cd5zZV1zDSxebE1hJ58cdjTOZY9jj4WffoIuXWwDvp49bcvwv/4KO7KExNWhLyK1VXVtOdcFm/NyKNBAVa+pRJk5WM0i6mvgCuBbYDKQo6qrRaQmsBJbeXkd8Gtl76nqvNhnZ32H/hdfQLdu1qE4YQK0bx92RM5ln3Xr4MEH4cYbrd+zQwcYNMg6/FNUVXTo31neRTXfqOrtlUkska8tUNW50ReWHApUdTYwDDgo8tGDgSGqOk9V/4jnXmXiynhFRbazJEC/fp5YnAtL9erWgvDddzYfZupU64+5/XYbXZZm4k0uPURkoIjcKSLJSKsXAMeIyA3AEdj+MYnec2CLUf78M2y3nY29d86Fa9ddYexY2zdp3Tq45RbYf39LNmkk3maxPVR1jIhshm1nvBewHHhfVb8KOMYql7XNYgsXWk1l6VJ491046qiwI3LOxfrsMzj9dJtoWbcu3HsvXHghVAti/nviAm8WU9UxkeMyVR2kqhcB/YH2IjI+suSKzydJdTffbInl0EPhyCPDjsY5V1b37tbZ37evLXh5ySU26GbmzLAj26i4kkuZ9btqi8ix2M6Tj2Gz65sDTQOJ0FWNn36CAQOsnfehh0Ak7Iicc+Vp3NgWvPy//4OcHBgxAnbayebGpPDEy7g79EXkFBEZDBQAg7Ehw+cBTVX1eFV9M6ggXcBU4bLLbEb+hRfCDjuEHZFzbmP+/W9bnun442HlSltFo2dPmJeaY5TiTS4XAc8CjYBfgIeAS1X1VVVdGVRwroq884799rPllnDrrWFH45zbVDk58OabtmR/48bw4Yew447wyispV4uJN7m8DbRQ1Z6quieWaHqJyBMicqWItA4uRBeooiK47jo7v+02+wfqnEsfIjbReeJEm8m/dCmceqqttrx8edjR/S3e5HKpqi6OvlHV6cAn2Iixq4HfROThAOJzQXvuORvS2L697YbnnEtPzZvDsGHW99KwofXJdO9uSzmlgHiTSwsAEdktshHXNGAscAjwINBOVS8LKEYXlJUrS5vB7roLatYMNRznXIJE4Oyz4fvvoW1b24epSxeYMyfsyCq/cGXE2yJShK3/NQl4AXhdVX8NKjBXBR5+2NYu2nNP6xR0zmWGdu1sRfPDDrMlnPbbz/Zj6tgxtJDirbnUxBLKjqq6k6r298SS4goKSpfRv+8+H3rsXKbZemtbJ3Cffazmsv/+MH58aOHEm1xOU9VbVHVyoNG4qtO/P6xYYR2AXX1+q3MZqXFjq7H06GErcHTrBp9/Hkoo8c7Q/zDoQFwVmjnT1hATgbvvDjsa51xVql8fhg6FE0+0Xyh79LDlnZIsNRaocVWrf38bgtynD+y8c9jROOeqWq1aNvfl/PNh7Vo47rikJxhPLpnut9/g+edtmZdbbgk7GudcslSvbi0W111nqysneQMyTy6Z7s477R9W376+V4tz2UbEWi52392Wibn33qQ9ukqSS2T+y8uR2fptquIZbhNMm2YL3lWvDjfdFHY0zrkwVKsGjz5q5/ffn7QVlauq5nI+cLCqPgjsKyJ1q+g5riK3326LU555pk2wcs5lp333tT7XNWtst8skiGuzsI0WavvVN1XVeSJSDWigqqmz6E0ZGblZ2C+/2Fap1apZDaZNm7Ajcs6Fac4cm1S5erUNTw5gSkLgm4WVKbyhiNwrIvdG3u8JnBTdq15VS1I5sWSs/v2t1nLWWZ5YnHOwzTali9Zeeqn1xVahIJrFXgD2AxoDqOp3QG0RuT+Asl08pk+HwYOtr+X668OOxjmXKq66Clq1gh9/hIEDq/RRQSSXIlXdH1tjLGo2cFYAZbt43Huv/VZy6qlea3HOlapXzzr1AW64webAVJEgksvPkaOCbXsMXAOk5vZomW7ePHjhBRuC6LUW51xZvXrZtISCAtsTpooEkVyGRbY77ikiTwO/AbsDFwZQtqusBx6AwkLbOCjEFVGdcylKBHIjffA//FBlj0k4uajqOCAPeB6YDlwPbKeqXyZatqukggIYMMDOvdbinNuQXXe1YxUml3j3c/mbiJyuqi8Cr8Vc6ygi+6rqe4mW7yrh0UdtmOERR8Buu4UdjXMuVaVychGRLoAAh4jIjDK3twAeBTy5JMvKlfDEE3butRbnXEV22cWOP/5oUxaqBT+fPpGay1LgJWBbbChyrELg2QTKdpU1cCAsWQJ772270Dnn3IZstZVtLjZ/PsyYYTtZBizu5KKqE0SkK9BbVddLJCKypaouSjg6t2mKi+Ghh+z86qt9l0nn3Mbtuqsllx9+qJLkklBdSFWXq+qzItJYRLYRkVYi0grYVkSqdoaOKzVkiC1G1749HHNM2NE459JBFfe7BNGhfz9wJTbPJforswITEi3bbQLV0klRV15ps/Kdc25jdtzRjj//XPHn4pRwcgG2AZoCXYAxqjpXRA4GagZQttuYL7+EceMgJwdOOy3saJxz6SL6i2gVdOZDMJMoJ6rqQmxkWPSn20jgoQDKdhsT3afhggugru9s4JzbRKtX27GKfm4EUXOpKSK/An2BaSIyAmiE73JZ9ebMgXfegRo14Lzzwo7GOZdOUj25qOotIvIGME1VvxaRhcBuwFsJR+cq9tRTtkDlSSdB8+ZhR+OcSyepnlxE5C7gE1WdBKCqI4ARiZbrNmL1asjPt/NLLgk3Fudc+qni5BJE09UJwD82AxORrQMo223Ia6/BokWw++42cdI55yoj1WsuQH+gj4jUJ7LsPlAdOAc4JYDyXVmq8Nhjdn7xxT5p0jlXeWmQXP4DdACOozS5ADQLoGxXntGj4fvvoUkTOPHEsKNxzqWjNEgu9wGfqeqS2IsicnwAZbvyRIcf5+VBnTrhxuKcS0/LI70Z9etXSfFB7OcyJJpYIkvAVI9eT7RsV45582y5l+rVbW6Lc87FY+pUO263XZUUH/RclClA8CuguVIDBtjw4+OOg5Ytw47GOZeOSkpg8mQ732GHKnlE0MnFe5arUnExPPecnV/ou0g75+I0ezb89Rc0awZbbFElj/BZ9Olk+HD4/Xdb/fiAA8KOxjmXrqq41gIBJBcROT3m7bnA/Mg2x0clWrYr49nItjnnnOPDj51z8Zs0yY6dO1fZI4Le5ng2tvSLb3MctPnzYdgwW0fMVz92ziUiCTWXlNrmOLKzZT62hP+LqnpZ5HpdbMjzCqAxcK2qLk/kXtp58UXryP/3v62d1Dnn4hWtuaRicqlom+N4iMhmWJLaG9gXeEdEhqnqJ8CTwHBVfV1EegADgJMjXxrvvfShun6TmHPOxUu1tOZShc1ioqob/1RlCxXZQlUXV/Jr6gJrNBKQiIwDrsKGN88AGqvqKhGpASwDOmM1pErfU9WZsc/Ozc3VsWPHxv8Hrmqffw4HHggtWsCsWb7bpHMuftOnQ7t20LQpLFiQUFEiMk5Vc8u7F8QMfUSkMdCA0qHITUXkAVU9a1PLUNXVMeXVB34GPsdqGgtVdVXkc8WRPp6uQFGc92Ym8MdNvmit5ayzPLE45xIzapQdq3jB2yCW3L8fuBJbVyyaXBSYEGd5XbDFMOcD9YAWwKIyH1sBNAeK47y3noKCAnJzS5NvXl4eeXl58YQfvCVL4K23bHTYWZucq51zrnwjR9oxzukM+fn55Ee3+4AmG/pcEDWXbbAO+C7AGFWdKyIHAzXjLO9X4AWsI/5+YDqwusxnamM1ExK497ecnBxStlns9ddh7Vo4+GBo0ybsaJxz6e7LL+3YpUtcXx77y3dkc8hyBTGJcqKqLsSGHUfHyI4EHoqnMFWdr6oDsf6WrsA8bKRXrIbA7wncSx8vv2zH00+v+HPOObcxf/xha4rVrw+77ValjwoiudQUkV+BPYBpIjIC+DqAssdiCWIE0DzS4Y+I1ARaAV8kcC89TJ8OX30F9erBsceGHY1zLt1Fay377gs1421c2jRBrIp8C3AMMFZV3wRuB14BDq5MOSJST0RiU+kRwCOq+gcwDDgocv1gYIiqzov3XqX/kGF55RU7HnccNGgQbizOufSXYH9LZQTRod+Q0uawa4G/gAJVnV3JonYE3heRX4BvgO9U9f3IvQuA+0RkF2wTstje9njvpTZVGDTIzk89NdxYnHOZIcH+lspIeJ6LiAwBtgImq2pe5No5QEdVvTrxEKteSs5z+e472Gsvm40/Z44t++Kcc/FasgS23NKaw5YuDWQHyormuQTR51KkqvsDk2KuzQZ83Gwioh35J5/sicU5l7iRI61FZI89qmxr41hBJJefI8fozPrawDVYZ7yLR1ERvPaanfftG24szrnM8MEHdjzssKQ8LohfiYeJyGBgCxHZATgSqA8cHUDZ2enTT6GgALbfHnbdNexonHPpTrU0ufTsmZRHJpxcVHWciOQBPbGhvtcDH6hq2dnxblO98YYdTzrJ921xziVuwgSYO9f6cJP0C2sgjfmqugJ4LfaaiJyuqi8GUX5WKSyEt9+28169wo3FOZcZ3o8MvD3iCKiWnA2Ig9iJ8mwRmSUia0VkXeRVAgwMIL7s8+mnNpJjxx2tWcw55xIVbRI74oikPTKImstD2PyRMdhikWALWJ4UQNnZJ9ok1rt3uHE45zLD4sUwerSNOj3kkKQ9Nojk8jHwjqquib0oIq8EUHZ2KSyEd96xc28Sc84F4aOPoKQEunaFRo2S9thKJxcRKbtuwEvAHSLyXpnrJwCXxBtYVvrkE2sS22kn6NQp7Gicc5kg2t+SpFFiUfHUXB9qpNQAACAASURBVO4BdgEWAiUx10+IORds1r4nl8p48007eq3FOReENWvgvcjv/WmQXG4C5qjq1Io+JCI+z6Uyiovh3Xft/IQTKv6sc85tig8/hGXLbPhxkltDKp1cVPXTTfzcu5UPJ4t99ZV1vHXo4KPEnHPBiK6sfsopSX90cgY8u42LduQfc0y4cTjnMsOyZTBsmE3EPin5g3c9uaQCVRg61M59UzDnXBCGDLEt0rt2hZYtk/74ICZR/mP/XRHpKCJHJVp21pg4EWbMgKZNbZl955xLVIhNYpDAPBcR6YKNCjtERGaUub0F8ChQdniyK0+0Sezoo6F69XBjcc6lv3nzYMQIqFUrtAFCiUyiXIrNcdkW2K/MvULg2QTKzi7RUWLe3+KcC8Lrr1tze8+esPnmoYQQd3JR1Qki0hXoraqeSOL1xx8wdizUqQMHHRR2NM65dKcKL0bWDA6pSQwS7HNR1eUbSiwickYiZWeN//3Pjt27J2V3OOdchvviC/jpJ+vDPfLI0MIIokP/XBGZKyJFZVZFfi6A+DJfdLXSww8PNw7nXGZ4+GE7XnAB1K4dWhhBLFz5AHAevipy5RUX26Jy4MnFOZe4336zPtxatSy5hCiI5PIpMExVV8ZeFJFBAZSd2b791haqbN8e2rULOxrnXLp77DHrc+nTB7baKtRQgkguLwF3ichbMdcEW8jy4gDKz1whbODjnMtQy5bBc5HeiEsvDTcWgkku/YB2wNGARq5VA5rhyaViH35oR28Sc84lauBAWLkSunWzhSpDFkRy6Q98oKpFsRdFpEcAZWeu+fPh++9thFjXrmFH45xLZ+vWwaOP2vnll4cbS0TCyUVVh27g1qxEy85ow4fb8cADbY6Lc87F6+23YeZM67tN8r4tGxLPTpSPAYNVdXTk/VNA2fFuAuwJdE44wkwVbRLz/hbnXCKKi+Hmm+388stTZgmpeGous7GlX6LqAGuBPyjtcxGsH8aVZ90629IYoIe3HjrnEvDSS/Dzz7DttnDuuWFH87d4Ngu7v8yl/sBsVS2MvSgizycSWEb74QdYsgTatPEhyM65+K1eDbfcYud33GHzW1JEEH0uv4pIfRE5EWgJTAHeV1Xvc9mQzz6zo68l5pxLxBNPwNy5sMsucPLJYUeznoSTi4jsAbyPDT+eiTWJ3SEivVR1cqLlZ6RPIztFd+8ebhzOufS1dCncdZed33MPVEutvR+DiOZe4AIgR1VzVXV3bAn+8wMoO/MUFsKXX9q5JxfnXLzuvdea17t1g8MOCzuafwhinss3qjok9oKqLhWRtQGUnXm+/RZWrYLOnaFZs7Cjcc6lo99/h0cesfN77gGRcOMpRxA1lwYi0jD6RkRqikhfYO8Ays483iTmnEvUzTdbZ/5xx6Xs1uhB1FyeBL4WkVVATaAtsBI4KoCyM080uXhnvnMuHiNH2hpiNWtC//5hR7NBQYwW+0VEdgd6Au2xeTDvqOrqRMvOOKtXW7OYiC/54pyrvDVrSueyXH89dOoUbjwVCGp4QVNgc6AIWBE5urLGjIGiIth559D2tXbOpbE77oCpU2H77aFfv7CjqVAQQ5FPB54GCoHpQC3gvshQ5EmJlp9RRo2y4/77hxuHcy79/Pgj3HeftXw8+2you0xuiqCGIj8ANFXV3VS1M3AIcFEAZWeW6BBkTy7OucpYtw7OOcfWEbvwQth337Aj2qggkssC4F5V/XvosarOA36OvheRegE8J72tWwejR9u5JxfnXGU88giMHQstW5ZOnExxQSSXK4GLRKRVzGsvYDsR2UZEWgNXB/Cc9DZxIixfbuuJtWwZdjTOuXQxfTrceKOdP/00NGxY8edTRBBDkR8GdgDKS6fRpjEFbgvgWenLm8Scc5W1bh2cfbaNND355JTZq2VTBFFzeRTYWlWrlfcCqgNXBPCc9Oad+c65yrrtNvj8c9hqK3j44bCjqZQg5rnkb+S+isiTm1KWiBwHPAg0Al4GrlTVYhFpgi3tvxSbqHltdFvleO8lXTS5dOkSyuOdc2nmo4/gzjttQcrBg6Fp07AjqpSkLKO5KT/QRaQVcCxwAnA5cA5wWeT2m8BAVb0WmArcHvOl8d5Lnnnz7LXZZik96ck5lyJ+/x1OPRVU4dZbbTv0NBNEn0tQWgNnqWoxME5EdgQOFJFRQEdV/TbyufeAKSJyJ7BTPPdU9a9k/aEA+O47O+65Z8oti+2cSzHFxXDSSVBQAIcckvKTJTckZX7SqeqXkcQSNQ9bSqY7MCPmc/Mip7kJ3EuuaHLZY4+kP9o5l2ZuvtkGADVvDi+/DNWrhx1RXKosuYjIGQkWsQc2878FsKjMvRVA8wTuJVdszcU55zbkww/h7rvTtp8lVsLJRUTOFZG5IlIkIusirxLguQTK3A4oUNUfsWHMZRfBrI2tXxbvvfUUFBSQm5v79ys/v8IxCpVTUmKTn8CTi3Nuw2bPhr597fzOO+GAA8KNZwPy8/P//lkJNNnQ54Loc3kAOA8YA0SbtQQ4KZ7CRKR6pLzrIpfmAR1i7gtQH/g9gXvrycnJYWw0AQRt6lSbPNmyJWy9ddU8wzmX3pYvhyOPhEWLoEcPuPbasCPaoLy8PPLy8gAQkYUb+lwQyeVTYJiqroy9KCKD4izvKuCBmOVkPmb9LZPbYItkjo0c47mXPN7f4pyrSFER9OoFEybYaNJXX82IgT9B/AleAu4SkQNiXl0prXlsMhHpB3wP1BWRtiJyFjZHZX5k9BhAD+AxVS1U1bHx3Iv/jxqHMWPs6E1izrmyVOE//7E5LTk58P770Lhx2FEFIoiaSz+gHXA01s8BlrSaARdvaiEiciNwR5nLv6jqQBE5EbhZRGYCDYAbYj4T773k8M5859yG3HuvLZ9fpw689x60bRt2RIERVd34pyoqQOQY4IOyEyVFpIeqDk+o8CTJzc3VKulzKS6GBg1g7VpYssQ3CHPOlXrtNVsvTATeeguOOy7siCpNRMaparnTO4JY/mWoiNQXkZOAlsAvWLJJi8RSpaZOtcTSpo0nFudcqVGj4PTT7fyBB9IysWxMEDtR7gG8jzWFzcRGit0Z2YlycqLlp7UffrDjrruGG4dzLnVMnQrHHAOFhbbx1+WXhx1RlQhqJ8oLgBxVzVXV3YH9WH+kVnb68Uc77rJLuHE451LDzJlw8MGweLENPX74YWsWy0BBdOh/o6pDYi+o6lIRWbuhL8ga0eTiNRfn3Jw5tgDlnDm2TfHgwVAjlZZ3DFYQNZcGIvL31mgiUlNE+gJ7B1B2eos2i3nNxbnsNm+eJZaZM23k6Acf2GCfDBZE2nwS+FpEVmF7prQFVgJHBVB2+vrjD1iwABo1sg5951x2+uMPOOgg+O03+Ne/YPhw234jwwUxWuwXEdkd6Am0x1YyfkdVy67rlV1i+1sytE3VObcRBQWWWKZMgZ13tsmSGTJJcmMCafCLLNXyf7HXRGR7Vf05iPLT0k8/2dGbxJzLTosWWef95Mmwww7wySew5ZZhR5U0lU4uIvIYMFhVR0feP4WtNrzex4A9gc4JR5iuJkdGYe+4Y8Wfc85lnsWL4dBD7ZfMDh3g009teZcsEk/NZTa23ldUHWAt8Aely78ItiRM9ooml+23DzcO51xyzZsHhx0GkyZBu3bw2WfQrFnYUSVdpZOLqt5f5tJXwIuxy7+ISEf+uVFX9lCFnyMtgjvsEG4szrnkmTbNaiwzZ9ovlh99BC1ahB1VKOLucxGRLlgNpRu2N31sr3Vj4GrgsYSiS1fz5sGKFdCkib2cc5nv++9tL5Y//ywdbpxFfSxlJdKhvxRbbn9bbEZ+rELg2QTKTm9ea3Euu4wcCUcdZZt+HXII/N//Zfw8lo2JO7mo6oTIvi29VTV7E0l5vL/Fuezx7rtw4omwZg307g0vvQS1y45xyj4JzdBX1eXlJRYR2SKRctNetObiycW5zPbii7ai8Zo1cP75toukJxYgmOVfEJHGIrKNiLQSkVZAWxEZGETZaSlac/FmMecykyrcdReccQasWwc33QRPPgnVq4cdWcoIYsn9+4ErsWHI0U59BSYkWnba8pqLc5lrzRo45xx45RVbfeOhh+DSS8OOKuUEMUN/G6Ap0AUYo6pzReRgbJ2x7LN4MSxcCPXrZ+0QROcy1oIFcOyx8M039n/81Vfh6KPDjiolBdEsNlFVFwLvAadFro0EHgqg7PTz22923G47X1PMuUzy44+wxx6WWFq1gtGjPbFUIIjkUlNEfgX2AKaJyAjg64DKTj/R5NIuuxcocC6jvPsu7Lef7cWyzz7w3Xe2EKXboCBWRb5FRN4Apqnq1yKyCNgVeCvh6NLRr7/acbvtwo3DOZc4Vbj/frjuOjs/5RR49lmoUyfsyFJeILULVZ2kqoWR889U9b9A9yDKTjvR5OI1F+fS26pVcPrpcO21llj694dBgzyxbKIgRoudC9wCbMU/k9ULiZafdmL7XJxz6WnqVDjhBJgwAerVs6Ry3HFhR5VWghgt9gBwHjAGKI5cE+CkAMpOP94s5lx6e+stOOssWx+wQwcYMsS3zohDEMnlU2CYqq6MvSgigwIoO7389ZdtaVqrlg9Ddi7dFBZaE9jDD9v73r2tf6Vhw3DjSlNBJJeXgLtEJLYDX4ATgIsDKD99RJvE2rb1mbrOpZO5cy2ZfP011KwJDz4IF13k0wkSEERy6YdtDHY0pZuFVQOakW3JxTvznUs/H38MffrY5OdttoE33oC99w47qrQXRHLpD3wQu1kYgIj0CKDs9DJzph233TbUMJxzm6CwEG69Fe65x0aDHXqoLeniezAFIoihyNuWTSwRxSJyu4jsFcAz0sPs2XZs3TrcOJxzFZsyBfbdF+6+25q+br3VNvfyxBKYIJLL+SIySkTGicg1ACJSF3gHGAB0FJEDA3hO6pszx46tWoUbh3OufKowYAD8618wbpz9IvjFF3DLLd5PGrAgkss2wA/AQKCpiPQFWgD1VHUeMAjrl8l80ZrLNtuEG4dz7p8KCmzRyfPPtwmSp55q64Xtv3/YkWWkQDr0VfWR6BsR6QPUwrY6RlVVRDoE8JzUF00uXnNxLrUMHw5nnmlTBTbbDJ5+Gk7Kzql4yRJEcmklIrsBa4CDsAUs5wLLwDYSAzJ/a7Y1a+DPP6FGDWjWLOxonHNgc8+uuw4ef9zed+1q2xD7L4BVLojk8irWv7INNkv/IqyJ7BsRuRHoBHwcwHNS29y5dmzRwttunUsFI0bYpl7Tp9svfXfcAVdf7f8/kySIVZHHAa1FJEdVCwBEpJuqLhKRHCzZPJXoc1KeN4k5lxpWrIBrrrGmL7Cl8Z9/3jrxXdIEsipypOmrjoi0EpFWwLYiMlBVC1T1FlX9I4jnpDRPLs6F73//s3XAnn7aZtrfdhuMGeOJJQRBrIp8P3AlNjs/ulaCAhMSLTut+DBk58KzdClceSUMHGjvc3PtfKedwo0riwXR57IN0BToAoxR1bkicjBQM4Cy00c0ubRsGW4czmWboUPhP/+B33+H2rWttnLlldbP4kITRLPYRFVdCLwHnBa5NhJ4KICy08f8+XZs3jzcOJzLFrNmwTHH2NyV33+37Yd/+MFWNvbEErogkktNEfkVG4I8TURGAF8HVHb6+CPSrbT11uHG4VymKyqyrYd32MH2tm/YEB55BL78Ejp1Cjs6FxHEaLFbROQNYJqqfi0iC4HdgLc28qWZJZpcfI6Lc1Xnq69shv3Eifa+d2946CFvMUhBCdcuROR0VZ2kqtEZ+SOA94FdEi07bZSUwIIFdr7VVuHG4lwmWrTI5qzsv78llnbtbNb96697YklRcddcRKQLNjrsEBGZUeZ2Y+AxrB8m8y1ebFX1zTeHOnXCjsa5zLFunY36uv56SzC1almfyvXXQ926YUfnKpBIs9hSbBfKbYH9ytwrBJ5NoOz04v0tzgVv5Ei49FLrpAfo3h2efBI6dgw3LrdJ4k4uqjpBRLoCvVU1sEQiIocAdwInqurMyLW6wH3ACqxWdK2qLk/kXqC8v8W54Myebcu0vPGGvW/Vyjrwe/XybYfTSFx9LiLSJ7L6cU9gROTaTiLytYgsEZFXRWSzOMrdCmgA7Fnm1pPAKFXtBwzF9olJ9F5wosOQPbk4F79Vq2zTro4dLbHUrWtzVn7+2TruPbGklXg79F8GDsF+cP8WWUPsU2xC5aXY7PyHK1uoqi7AksDfRKQ50IfS/ptPgKNFpE289yob10Z5zcW5+KnC4ME2jPi222yF8ZNPtt0ib74Z6tULO0IXh3ibxX5S1TNj3t8DbA7srqoTAETk+XgKVtUSWf83lG7AQlVdFblfHBlA0BUoivPezHhi2yDvc3EuPl98YU1gY8bY+912g0cf9Q28MkC8yeW76ImI7AWcATwWTSwRQa1r3QJYVObaCqA5UBznvWD9+acdc3ICL9q5jDRpku2zMmyYvW/WzJbEP/NMXxI/Q8SbXLYUkR2BVVgT2WzgxuhNEdka648JggKry1yrjdVMSODe3woKCsjNzf37fV5eHnl5eZse4aJIDttyy03/Guey0bx5tl/988/b/LAGDWxo8eWXQ/36YUfnNkF+fj75+fnRt0029Ll4k8sN2AZhHYApwDGqulJEWgO9gIuxZrIgzMNGesVqCPyOJZ547q0nJyeHsWPHxh+hJxfnKrZ8Odx3H/z3v7B6ta39dcEF1qfStGnY0blKiP3lO7IiS7ni6tBX1V9UtROwpapur6qTI7fmYJMnOwBB9cKNAJpHhhUjIjWBVsAXCdwLVjS5NNlgEncuO/31F9x7L2y7LfTvb4nl+OOtWezxxz2xZLCEln9R1SVl3peo6troK54ypbQ3XyJl/gEMAw6KXD8YGKKq8+K9F09cFfKai3PrW7vWOubbtbO+lcWLoUsXGD0a3noLOnQIO0JXxVJqXWoRaQicGnl7hog8Htk6+QLgPhHZBWgGxHaIxHsvGMXFtlGRiC3/4lw2KyqCF1+E228v3eMoN9dqLYcc4nNVsoioatgxhC43N1fj7nMpKLCq/RZblNZgnMs2xcXw2ms2T+XXX+3aTjvZCLCjj/akkqFEZJyq5pZ3L6VqLmnJm8RcNisqgldesZpJNKm0b281l969oVp2bevkSnlySZQnF5eNCgvhpZfgrrtgRmRR9Hbt4IYboG9f3wnSeXJJmCcXl03WrrU5KnffbQtMgnXO33ijLdniScVF+L+ERHlycdlgxQp45hmbpzIvMuBy++3hppus+ctn1bsyPLkkatkyO/pIMZeJCgpsSPHjj9uoSLCO+ptusvkq3qfiNsCTS6KWR7aHadQo3DicC9LMmfDgg/DcczbxEWwxyeuugyOO8NFfbqM8uSTKk4vLJOPHW1J5/XXbYhjgyCNt/S9fqdhVgieXRHlycemupMRWJ/7vf20JfLA+lL594ZprYMcdw43PpSVPLony5OLS1V9/2Wz6hx+GadPsWqNGcO65cPHF0Lp1uPG5tObJJVGeXFy6mT0bnnoK8vNtzS+wRHLZZXDWWf5v2QXCk0uiPLm4dKAKn38Ojz0GQ4daUxjA3nvDFVfAv//tc1RcoPxfU6I8ubhUtnIlvPyyDSWeNMmu1awJJ50EF10E++wTbnwuY3lySZQnF5eKfvzRJj2+/HLpXKxmzeD88yEvD7beOtz4XMbz5JIoTy4uVaxcaUOI8/Phu+9Kr++zj3XQH3881KoVXnwuq3hySYRqaXJp2DDcWFz2Gj/eEsqrr9oyLQCbbQannmojv3bZJdz4XFby5JKINWtsolmtWv4boUuu5cth8GBLKuPHl17fbz9LKL16Qb2gdhp3rvI8uSRi1So7+n9ilwyq1tz1zDOWWKL//ho3htNPh3POgc6dw43RuQhPLomIrrlUt264cbjMNneubcg1aFDpiC+Arl2tlnL88VCnTnjxOVcOTy6J8OTiqsry5TBkiI32GjHCai0ATZrAGWdYLaVjx1BDdK4inlwSEU0u3izmglBUBB99ZAnlnXesTw+gdm3bh/7UU6FHD+/fc2nBk0sivObiEqUK48ZZk9fgwbZ/StQBB9jikSec4PsFubTjySURnlxcvCZNgjfftHkpv/xSer1TJ0soffpAmzahhedcojy5JMKTi6uMSZPgjTcsqfz8c+n1nBzbf75vX9h9d9+Iy2UETy6J8OTiKqJaWkMpm1AaN7bFInv1goMOsvW+nMsgnlwSEZ1n4MnFRanCxImlCSW2yWuLLUoTSvfunlBcRvPkkgivuTiwUV6jRtlS9u++CzNmlN7bcsvShHLggZ5QXNbw5JIITy7Za/lyGD7cEsoHH8DSpaX3cnLgmGOgd2/o1s0TistKnlwSsXatHWvXDjcOlxxz5ljN5N13bWJjUVHpvU6dbC7KMcfAXnvZHvTOZTFPLokoLraj/2aamdauteau4cPhww/XX3qlWjXo0sUSytFHQ4cO4cXpXAry5JIITy6ZZ8YMSyTDh8Nnn8Fff5Xea9AADjnEaic9e9pSLM65cnlySUS0WcT3Hk9fy5fDyJHw8ceWUKZOXf/+zjvbkis9ethy9r70inObxH8qJiJac/Hkkj7WroVvvoFPPoFPP7Ul7NetK72/2WZWOzn8cDjsMGjRIrxYnUtj/lMxEd4slvqKi+GHH6yJ69NP4csvS0f5gXW877OPTWTs0cM64/2XBecS5v+LEuHNYqlnzRqrjXz5pTV3jR5te8vH2nFHSyYHH2yLQzZqFE6szmUw/6mYCK+5hG/5cksgI0daQvnuOygsXP8z7drZxloHHWQz45s1CydW57KIJ5dEeM0luVRh1iz49lv4+mtLJj/8ACUlpZ8RsU74Ll1KX82bhxezc1nKfyomwjv0q9aKFTB2rHXAf/utHRcsWP8zNWrAnntaEjngABvR1bhxOPE65/7mPxUT4c1iwSkshMmTLZl8+629Jk1av1YCtvjj3ntbx/v++9uxfv1wYnbObZAnl0R4s1h8Vq+GCRNg/PjS14QJ/+wrqVED/vUvSyB7722vdu18vxPn0oD/VEyE11w2bvly+Omn9RPJ5Mnrzy2Jat8edtutNJnstpsvCupcmvLkkgjvcym1erXtXTJx4vqv2bP/+dlq1aBzZ6uVRF+77upDgp3LIP5TMRFvvWW/gWdTM01hIfz6q/WHxCaRX3/9Z/8I2IrRZRPJTjtBvXrJj905lzSeXBIhkpm1lpIS+P13W2dryhQ7Rl8zZpTfpFW9Omy/vU1QjH21bZuZ3yPnXIX8f322Ki6GuXMtWcyYAdOnlyaQadNKt3AuSwS23fafSaRjR9/Xxjn3t4xPLiJSF7gPWAE0Bq5V1eXhRpUE69bZnJCZMy15RI/R15w5pX1G5Wna1PYoib46drRj27ZQp06y/hTOuTSV8ckFeBIYrqqvi0gPYABwcsgxJWbNGpg3r/Q1d+4/j/Pnl998Fat5c6uFRF/RRNK+vU9EdM4lJKOTi4g0B/oAF0YufQIMEZE2qjoztMBilZTYTPRly2DJEigosNeff274uGzZppXdpAm0br1+Aom+Wrf2GohzrspkdHIBugELVXUVgKoWi8gMoCsws8qfXlICV10FDz0UbLk1atg+I9FXy5b/PG/e3PtAnHOhyfTk0gJYVObaCmC9lQwLCgrIzc39+31eXh55eXmJP/233yqfWLp1s/6Opk0hJ2f9Y/R8881trohzziVZfn4++fn50bcb3OtbVDU5EYVARK4CeqnqXjHXxgOvquoD0Wu5ubk6duzY4ANQhccfhwcftHkdm21W+tp8czs2amTDeMGaqvr0CT4O55yrAiIyTlVzy7uX6TWXedgIsVgNgd+T8nQRuPhieznnXBbJ9LaVEUDzyHBkRKQm0Ar4ItSonHMuw2V0clHVP4BhwEGRSwcDQ1R1XnhROedc5sv0ZjGAC4D7RGQXoBkQQE+9c865imR8clHVJcC5YcfhnHPZJKObxZxzzoXDk4tzzrnAeXJxzjkXOE8uzjnnAufJxTnnXOA8uSQoZo0dl2T+vQ+Pf+/Dky7fe08uCUqXv+hM5N/78Pj3Pjzp8r335OKccy5wGb0q8qYSkQJgVpxf3gRYGGA4btP59z48/r0PTyp971urak55Nzy5OOecC5w3iznnnAucJxfnnHOB8+TinHMucJ5c4iQidUXkMRG5S0SeEpFGYceULUTkOBGZISKLROQREcn41b1TjYjUEpEfRaRb2LFkGxHZW0SuEpFjRWSDe9iHzf9Txu9JYLiqvi4iPYABwMkhx5TxRKQVcCxwAtAZeAqYAzwQZlxZ6GqgTdhBZBsRORNop6o3hh3LxnhyiYOINAf6ABdGLn0CDBGRNqo6M7TAskNr4CxVLQbGiciOwIF4ckkaEdkXmA8sCTuWbCIiXYFTgEPCjmVTeLNYfLoBC1V1FUDkB90MoGuYQWUDVf0y8v2OmgfMDiuebCMi9YETVHVg2LFkof8CPwOPichwEdkv7IAq4sklPi2ARWWurQCahxBLttsDeDrsILLItcA9YQeRbUSkI7Ab8IyqXgR8CgwXka3CjWzDPLnER4HVZa7VBopCiCVrich2QIGq/hh2LNlARA4Dxqrqn2HHkoU6A4tV9afI+8cBAY4JL6SKeXKJzzygcZlrDYHfQ4glK4lIdeA84LqwY8kiVwEDRWShiCwEtgGGisg1IceVDWoQ00euqquBacCWoUW0Ed6hH58RwDMiUldVV4tITaAV8EXIcWWTq4AHVHVt2IFkkT5YDT3qa+AK4H/hhJNVfgI2E5GtVHVB5FoxMDnEmCrkNZc4qOofwDDgoMilg4EhqjovvKiyh4j0A74H6opIWxE5K9JE5qqQqhao6tzoC1iHNUsuDzu2TKeqvwAfYkPwEZHNgVrA+2HGVRGvucTvAuA+EdkFaAbkhRxPVhCRG4E7ylz+xUcvuSxwGvCIiNTCWkpONwsUagAABwFJREFULjNyMqX4qsjOOecC581izjnnAufJxTnnXOA8uTjnnAucJxfnnHOB8+TinHMucJ5cnHPOBc6Ti3NxEJFrROSJsONwLlV5cnEpRUS2F5H/E5FnIzsdqoicGnZc5fgAGBx2EM6lKk8uLtW8A9ymqudgS4y/FHtTRC4OIygRaSoiJ0bfq+pEVR0VRiyZJKy/T1f1PLm4lCEiWwAdgL8AVLUEuBnb4gAR2Rm4M4S4amNJrm6yn53JROQM4Liw43BVw5OLSxmquhiYgm2CtHfk2izgu8gP+DOARiJyj4h0FZHmIvJ8pP/jSxHpHi1LRHqIyF0icrOIrBOR70Xk1Mgil7+IyAEi8pmILBGRIyNfU0NEHheR60RkmIhcESnuKGAHoJeIXCEi7UXkBREZEPO8/UXkYRF5QEQ+FpH2kesbfF4sEdlDRIaIyC0iki8iK0TkNRGpIyIdIn/OASLynIj8FrmeKyK3isjrkabEzSJlHS8iV4rIKyLyTeRaAxG5XUSuF5E5ke/PXiIyVURujXzmOBFZKyLdKvvMikTKfUREHozsoNhARLbFFmFsF/P3+WLkeRdGYhwhIltEvteLReS5yv2LcqFSVX/5K2Ve2KZIs4ASYCDQNOZeG/sn+/f7B4CbI+dXA+9FzhtjO4M2ibx/H3g1cp6D1YSux2oi1wNfRe4dCXwWOd8DWBHzrM+BMyLntYEBwAuR91thS6LXiLy/ClsKvXpFzyvz566BrbQ9AmgP7I7V4K4D6gD5kTJ3As4FNgNei/n6r4D7I+dfxFzvHzleCPSOnHcBekTOXwRujfn8LGwb70o9cyN/pwuBVpHzScDxkfMzgM/LfE9HYbXXRsAC4ClgC2BroBDoEPa/UX9t2strLi6lqOok7IfZ49gqsJMiK0+XJx94TkS2xmoWDSLX20fOo7uF/kxpk9bCyHGw2oZLY7BVrQG+BK4SkXrAPjHllY1xLfBHzKVTgalaukLts8D2wP4beV5smcWRz36hqtNUdRzwKnCkqq4B5gNjVHWCqj4DHAFsKSKXichlkT9j1DYi8rSINAUeiVxbia3ifRT2A3x89NFlQ4nEU9lnVqSXqs4W2/O9NuV8X2O+p7+q6lS1ZfynAt+q6mJVnR+5/4/vnUtNvuS+Sxliu0tuo6ozgUsizSAfYElkr3K+ZC62p/s04Btg28j1icBsbL+ddyPXXwer9ohIbBnFWA0DVV0mInsBXYGRGwk39odye6Dm3zdUl4rIYqBlRc/bBJOwJBd9XuwzWwEzVPXhcr7uWOAFYDpWo3sKGAS0Bd4ExgK9NuH5lXlmRZaKyL2RZ8/Ftufd0PNirSvzvgRvyk8b/hflUklNYn7oqeqPQD+g0wY+/zQwU1VfBtbGfN0q4ErgYBE5DXhZVV/b2MNF5Cxgb1V9CFhUibhnYzWn9YrD+o8SUQv77b0884GeIvL3IINIYgRYDuRiTXCPi0g7rFnqlkictYF7I58tYdOTXUXPLJfYplYfALer6thNfI7LAJ5cXKq5WESax7xvAQyPnBcCiEhjEemEDVXOiTRj7UfpzpT1sFFlrwLjgJ/ENlhCRKL/5mN/e46e74Y1+9QgssuoiLQTkTqRZ0efG/2a6Ne9CLQQka6Rr+kATFPVsRt5Xnm2jjnvBjwZ8z42CXwA1AeGichh8v/t3U9LVFEYx/Hv02yicCGBixpRcWEQ5CpoU0LOwp1QYC8gZitCRtRKIiSIKKhFSwO3BS4NwkULwYXQpk2MFS1EF4lvwKfF88gMMmLFhXsdfh84i3vmzzn3MMy553nO5Zo1aU/C8+5+4O6viSd2Atw1s7q7bwFPOvqwA1zLZP1tItcxkKvIv27TYhNEt/uRRokc2AUzGwMGgXNmNkx7TM/mceeYHjrpWKqq7KSPisphIZLIDuwBy0Ro5x3Qn68b8BH4QiR9m8A+sALcAXaBaeLKfB34TYShnFiJjOdnnNjiXAdeEaueKeJqf5v4Q75FhHCeZtv38vubxMaCz0TO4Uq+PpVtPiPyRfWsP7a9Lue/RORD7hNP25zN+jEief4LmOh4/wQRAtwDXgBnsv4gx+0hMJd1C0CLWAm+BIay/hIRPvuW57wJLBJhyH9pcxPY6nJONWCVmMQe5dhsAAPAReAHcW/T5Tz3r8BVIl+1C7wHhoCZHLe3QF/Zv1WVk4ueRCk9J1cXDXd/k8c1YhfapEfIq5LMbIkI8y2U3JX/YmYP3P152f2QalBYTHrRY2DEzM7ncQ1oEMn9KusWFjoVzOwm7R1oIppcpCctEju4fprZd+ADsOburXK7dTwzuwFcBybN7OjmgNOg5e6fyu6EVIfCYiIiUjitXEREpHCaXEREpHCaXEREpHCaXEREpHCaXEREpHCaXEREpHB/AKTCgt4g3UN8AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(6,6))\n", "plt.plot(vehicle1.stag_pres_atm[0:1700], vehicle1.q_stag_total[0:1700], 'r-', linewidth=2.0)\n", "plt.xlabel('Stagnation pressure, atm',fontsize=14)\n", "plt.ylabel('Stagnation point heat-rate, '+r'$W/cm^2$',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/decadal-2010-uranus-probe-thermal.png',bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-uranus-probe-thermal.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-uranus-probe-thermal.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }