{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 45 - ADEPT ViTaL - Venus (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 ADEPT ViTaL entry system concept. https://ieeexplore.ieee.org/document/6497176" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"VENUS\") \n", "planet.h_skip = 200.0E3\n", "planet.h_trap = 20.0E3\n", "planet.loadAtmosphereModel('../atmdata/Venus/venus-gram-avg.dat', 0, 1, 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle1=Vehicle('adept-vital', 1620.0, 38.0, 0.0, np.pi*6.0**2.0*0.25, 0.0, 3.0, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle1.setInitialState(200.0,0.0,0.0,10.8,0.0,-8.25,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": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZxcZZX/8c/pJZ1O0tlDQjbCEhABEWgQYRRkGSMw7LKqMAoRFX7iigo4M4iCgIogog2DDqDIIAqCgIIDAiJgRwybLIEQIEDS2Ruydvf5/fHc6lQqvVXXcu+t+r5fr/uqqnur7j1VnZw+/dS5zzV3R0REREREClMTdwAiIiIiIpVAhbWIiIiISBGosBYRERERKQIV1iIiIiIiRaDCWkRERESkCFRYi4iIiIgUQV3cARTL+PHjfcaMGXm9pq2tjQkTJpQmoDLTe0kmvZdkStp7mTNnzhJ3T05AZTCYnA3J+9kNRBpjhnTGrZjLJ41xFyvmvnJ2xRTWM2bMoLW1Na/XNDc35/2apNJ7SSa9l2RK2nsxswVxx1Bug8nZkLyf3UCkMWZIZ9yKuXzSGHexYu4rZ6sVRERERESkCFRYi4iIiIgUQVUX1rNnz447hKLRe0kmvZdkqqT3Um3S+LNLY8yQzrgVc/mkMe5yxGzuXvKDdB/M7Gjge8BI4EbgS+7eYWbjgW8DK4B64Bx33xC9ptdt2Zqbmz1tvT4iIgBmNsfdm+OOI5dytojI5vrK2WUbsTaz6cCRwLHAF4DTgLOjzbcA17n7OcALwAVZL+1rm4iIlIBytohI/so5K8hWwCfdvQOYY2Y7Ax8ys4eBHdz9seh5dwDPm9mFwC69bXP3d8oYu4hIrzo7O3F36uoqZqIlSFjOdoeurnC/traQPYmIlE7ZRqzd/aEoQWcsBF4FDgDmZz1vYXS3uZ9tIiKxu/nmm6mrq+Pkk0+OO5SiSlrOfuQRqKuDD36w0D2JiJROnCcv7gn8BJgCLM3Z1g5M7mfbJtra2mhubu5eWlpaShCyiMimMqPUGzZs1kbcp5aWlu58BYwvfmRFp5wtIlVroDk7lu8tzWw7oM3d55qZA2tyntIAbAD62raJCRMmpG6ichFJv/r6eiD/wnr27NndZ6ib2ZKiB1ZEScrZZTzfXkSk20BzdtlHrM2sFvg08LVo1UJgTNZ2A4YDb/SzTUQkdpnCuqOjo59nplNScrZZoXsQESm9OFpBvgxc5u7rosf3AjtkbZ8BrAda+9kmIhK7wbaCpIhytojIAJW1sDazbwBPAI1mto2ZfZIw1+mb0RnnALOAK919vbu39ratnHGLiPRmsK0gaZDEnK1WEBFJsrL1WJvZecC3clY/5+7XmdnxwDfN7BVgBHBu1nP62iYiEqtKbQVJWs5WK4iIpEHZCmt3vxC4sJdtC4BP5btNRCRuldoKopwtIpK/OKfbExFJvUpuBUkitYKISJKpsBYRKYAK6/JQK4iIpIEKaxGRAmRaQSqtxzqpNGItIkmmwlpEpAAasRYRkQwV1iIiBVBhXR5qBRGRNFBhLSJSgEqdbi+p1AoiIkmmwlpEpACVOt1e0mjEWkTSQIW1iEgB1AoiIiIZKqxFRAqgVpDyUiuIiCSZCmsRkQKoFaQ81AoiImmgwlpEpABqBRERkQwV1iIiBaitrcXMcHe6urriDqfiqRVERJJMhbWISIHUDlJ6agURkTRQYS0iUiC1g4iICKiwFhEpmArr8lEriIgkmQprEZECZVpBNOVe6agVRETSQIW1iEiBNGJdPhqxFpEkU2EtIlIgFdalpxFrEUkDFdYiIgXS1RdFRARUWIuIFEzT7ZWPWkFEJMlUWIuIFEitIKWnVhARSQMV1iIiBVIriIiIgAprEZGCqRWkfNQKIiJJpsJaRKRAagUpPbWCiEgaqLAWESmQWkFERARUWIuIFEytIOWjVhARSTIV1iIiBVIrSOmpFURE0kCFtYhIgVRYl49GrEUkyVRYi4gUSD3WpacRaxFJAxXWIiIFUo+1iIiACmsRkYKpFaR81AoiIkmmwlpEpEBqBSk9tYKISBqosBYRKZBaQUREBFRYi4gUTK0g5aNWEBFJMhXWIiIFUitI6akVRETSQIW1iEiB1AoiIiKgwlpEpGBqBSkftYKISJKpsBYRKZBaQUpPrSAikgYqrEVECqRWkPLRiLWIJFnZC2szO9jMHjOzGTnrLzQzj5a5WesbzexKM/uOmV1tZiPLHbOISF8quRUkKTlbI9YikgZ15TyYmU0ERgB75awfDowHDo5WLcja/GPgHne/2cxmAT8FTixDuCIiA1KphbVytohIfspaWLv7IjO7vYdNpwPPAQ+7+9rMSjObDJwEfC5adR9wq5nNcPdXSh2viMhAVGqPdRJztlpBRCTJyt4K4u5d2Y/NzIATgO8Bi8wse2Rjf2CJu6+OXtsBzAf2K0+0IiL9q+Qe66TkbLWCiEgaxH7yogd7E75WvBL4hZl9JNo8BVia85J2YHLuftra2mhubu5eWlpaShq3iEjGYFpBWlpauvMVIf+lgnK2iFSjgebssraC9MXdlwPnRaMhZwN3Aw6syXlqA7DZb68JEybQ2tpa8jhFRHINphVk9uzZzJ49GwAzW1KSwEoorpytVhARicNAc3bsI9Y9uAKYHt1fCIzJ2d4EvFHWiERE+lDJrSADUJacrVYQEUmDJBbWXcDfo/v3A5PNrBHAzOoJCfzPMcUmIrKZSp0VZICUs0VEInHMY50Zd7Do8WQzO8nMaqJtXwLOBXD3t4A7gQOj1xwE3OruC8sctohIryp1VhBIXs5WK4iIJFlRCmszO2SAz2sCzogenmpmEwhfG14MPAlcDfw8Z1qmzwBHmNm5wCHA7GLELCJSLGlrBUljzlYriIikQV4nL0aT/Z8HTGJjUV4DTAQa+3u9u7cTEvHVWavb2Nif19NrlhPmTBURSaSktoJUYs7WiLWIJFm+s4LcAFwEPEHoq4Pw9eBJxQxKRCRNEtwKUjE5WyPWIpIG+RbWL7v793NXmtlzRYpHRCR1EtwKopwtIlJG+fZYf8bMju1h/eHFCEZEJI2S2gpCBeZstYKISJLlW1jvC1xrZp1ZSxeb9t+JiFSVBBfWFZOz1QoiImmQbyvIBYQzvFuBzmidAR8vZlAiImmS4B5r5WwRkTLKt7C+G7jd3ddlrzSza4oXkohIuiS4x7ricrZaQUQkyfItrH8FfN3M/i9rnQFHA58vWlQiIimS4FaQisnZagURkTTIt7A+H5gJnApkxg2MMEdqqpK0iEixJLgVRDlbRKSM8i2sLwXuji4a0G2gV/ESEalECW4FqbicrVYQEUmyfGcFeS03QUcmFSMYEZE0SnArSMXkbLWCiEga5FtYfzZ3hZl9FLi8OOGIiKRPgltBKi5na8RaRJIs38J6rJmdD2Bmw8zsWuBa4G9Fj0xEJCUS3ApSMTlbI9YikgZ59Vi7+6FmtreZXQl8GHgL2BVYWYrgRETSIKmtIMrZIiLl1WdhbWYf7GXTAuBF4DJgK+BI4AvFDU1EJB2S0gpSDTlbrSAikmT9jVhfAuwIrAC6eth+HWHqpi1JaZIWESlUgkasKzZnqxVERNKgv8L6AqDV3Rf39SQzO7Z4IYmIpEumx7qjowN3x+KrApWzRURi1Gdh7e53DWQn7v7r4oQjIpI+ZkZtbS2dnZ10dHR0j2CXWzXkbLWCiEiS5TsriIiI9CApfdaVSq0gIpIGKqxFRIogwVPuiYhImaiwFhEpggSdwFiRhgwJt2vXxhuHiEhfCi6szWw3M7vRzL5kZjMKD0lEJH3S0gqS1pw9Zky4XbEi3jhERPpSjBHrM4CD3P17wD5m1liEfYqIpEqKWkFSmbOHDg2j1uvWwZo1cUcjItKzYhTWZwJ7RPd/BcRzOryISIxS1AqSypxttnHUeunSeGMREelN3oW1mX3EzE6O7u8FzHT3hQDu3uXuq4oco4hI4iW1FaSScvaMGeH2pZdiDUNEpFd5FdZm9m3gFuAoAHd/HPiEmZ1UgthERFIjia0glZazd9wx3D77bLxxiIj0Jt8R632AqcBDWevuBi4uWkQiIimU0FaQisrZzc3h9v77441DRKQ3+RbW97r7CiD72le6NK6IVL2EFtYVlbMPOSTc3nMPrEpNA4uIVJN8C+uXzewbwA5mdoKZ3Qp8Dvhu8UMTEUmPNdFUFY2NiZpko6Jy9tZbw377QXs7XHFF3NGIiGwur8La3X8F/B9QC3wCWAH8q7tfVYLYRERSY1U0hDpy5MiYI9moEnP2N78Zbi+8EJ54It5YRERy1eX7And/FHg089jMas3sQHf/U1EjExFJkSQW1lB5OfuAA+C00+Daa2HWLLj7bth997ijEhEJ+iyszex3wIh+9jEeWAKkMkmLiBRq1apVtLe309DQQFNTU2xxVEvO/tGPYMECuPde2Hdf+M534MwzoT4VM3KLSCXrb8T6yeg5zwNdwOnAE0Br1nO2BsaWJDoRkRR44YUXAJg5cyY1NcW47tagVUXObmiAO+4IxfS118IXvwg/+Ql85SvwsY+FqzSKiMShv8L6UsCis8oxs23d/Zu5TzKz20sRnIhIGvz1r38F4D3veU/MkVRPzm5ogGuugSOOgLPPhhdegNNPh3PPhU98Av793+Hd7447ShGpNn0Orbj7ykyCjkw0s9rs55jZ/sBeJYhNRCQVfvvb3wIwa9asWOOoxpx92GHwz3/CL34Bu+4KixfDZZfBTjuFx+eeC488Agm7IKaIVKh8T178NTDXzB4G2oEdgQ8Dm42IiIhUg9bWVu6//35GjBjBYYcdFnc4uaoiZ9fXw0knwYknwqOPws9+BjffDE8+GZbvfAeGD4e99w492f/yL7DHHjA21Q0xIpJEeRXW7n6vmf0r8FngPYSpm45199R/rSgikq81a9Zw+umnA3DGGWcwZsyYmCPaVLXlbDN4//vDcuWV8OCDcOedcNddMG8e/OlPYcmYOjWMameWHXaAbbeFEf2d/iki0gtz9/6f1dcOzGqA3dx9TnFCGpzm5mZvbW3t/4kiIkXwzjvvcPzxx/P73/+erbfemieffJIRg6zIzGyOuzcXOcTejlWVOfutt+AvfwnLI4+Ekezomj6b2XJL2G47mDkz3G69NWy1FUyfHrbFe36qiMStr5yd14i1mV3Xw+oJQCdw5CBiExFJnX/84x+ccsopPPnkk4wdO5Y777xz0EV1KSlnbzRpEhxzTFgAOjvDKPbcuWF56il48UV4+WV4882wPPTQ5vuprw8j3ZlCe6utNt6fPDkcZ+zYMHouItUn3x7rmYS5T7OHuacCPaQfEZHKMn/+fC655BJaWlro6upi++2357bbbmPHHXeMO7TeKGf3orY2tH7ssAMcd9zG9Z2d8Nproeh+8cVw+8or8OqrYe7stjaYPz8svamvDwX2lluG2+z7EyfCuHEwfnxYxo6Furwv1SYiSZXvf+fT3P357BVm9i7g1IHuwMwOBi4Ejnf3V6J1jcAlhJNrxgDnuPuq/raJiJTa+vXr+cMf/sD111/Pb3/7Wzo7O6mpqeHzn/88F1xwQeKutJhDOTtPtbUwY0ZYDjpo8+1r1oQiO1NoZ9+++WZoOVm5MhTnr702sGOOHr2x0M4uujP3x44Nz8leRo1SQS6SRPmevPh8D6vXAacBX+vv9WY2kXBVsNypnn4M3OPuN5vZLOCnwIkD2CYiUnRLlizhj3/8I3fffTd33303S5cuBaCuro5PfOITnHPOObw7BZMkK2cXX2PjxpHu3qxZEwrst97aWGxn7re1wZIlG5dly2DFirDMm5dfLCNGbF5s5xbgmfVNTZsvI0fCkCGFfR4isqm8Tl40s/ls+pViLTARuMHdTx/gPmoI/X1bu/srZjYZmA+McffVZlYHrAR2Atb3ti0zcpKhkxdFZDA6Ozt54YUX+Otf/9q9PPvss2Tnxp122omPf/zjnHzyyUydOrXoMZTq5EXl7OTr7AxF9ZIlsHTppkX30qWhEM8U3tnLypXQ1VX48YcM6bnozi3Ac9cNGxamMBw+fPP7urS8VLqinbwI3Av8ko2JugtY5O4vDHQH7t5lm57VsT+wxN1XR9s7ol8G+wEb+tj2Sp6xi0gV6+rqYuHChTz//PM89dRT3cszzzzDmpzpIYYMGcJ+++3HrFmzmDVrFjvuuCOWzrPRlLMTrrY2tHyMG5ff67q64O23Ny+2eyvC29th1apwm72sXx8K+OhLmaKor9+04O6rCO/t/rBh4duB3pYhQ3SCqCRTvoX1je7+YPYKM9vBzLZz9zy/xOo2Bcj9L90OTAY6+ti2iba2NpqbN/7xMHv2bGbPnj3IkEQkbdyd5cuX8/rrr/P6668zf/585s2bx0svvcS8efN4+eWXWbduXY+vnTJlCnvvvTfvf//7ef/738/uu+/O0KFDSxpvS0sLLS0tmYfjS3QY5ewKVVMTRpJHjgwzkgyGO6xbt3mxnVl6KsQzy+rV8M47G2+z72/YsLGoL5WaGhg6tP8CvFjL0KGaZrHaDTRn51tYzwIezFm3CLgf2C3PfWU4kDubaANh5IN+tnWbMGEC+lpRpPK4OytXrmTx4sUsWrSo+/aNN97g9ddf57XXXusuplevXt3nviZOnMh2223HzjvvzC677MIuu+zCzjvvzNgYLsGXXUia2ZISHUY5W3plFgrGoUNhwoTi7NM9jIL3VHDn3u9r2+rVoVe9t6WjIzynn//yRTVkSPisMoV29v1Srquv1+h8Egw0Zw+osDazMwlnhTeZ2Tk9POXewQQZWUg4czxbE/AGIYH3tk1EUsbdWbNmDcuWLetxWbp06WYF9OLFi1m/fv2A9j9ixAimTZvG1KlTmT59OjNnzmS77bZj2223Zdttt6WpqanE7zAZlLMlLmbQ0BCWUl6ItKOj78J7oMvatQN/7vr1YVlV5jluMqPzPRXgpSzuNUo/OAMqrN39R2b2IHAKcEXO5vXu/mYBMdwPXGNmje6+xszqgenAnwknzPS2TUTKzN1ZvXo1q1atYuXKlaxcubL7fu665cuX91g8D7RIztbU1MQWW2zBxIkTu28nTZrUXURPnTqVadOmJX3qu7JRzpZKV1e38UTKcsiMxGcK8bVrN73f17p8n5+7Lo7R+YzcUfpyFfRpHqUfcCuIuz9pZl9z982+0jOz6e7+6kD2YxvPgrFov2+Z2Z3AgcCdwEHAre6+MHp+r9tEpH8bNmzg7bffHvDS3t7ea7G8cuVKOjs7C4qnoaGBcePGMXbs2B6X3AJ6iy22oLGxsUifRvVQzhYpnuyR+FGjynvsjo7QC1+OIj53XRJH6QeybsSIMP/7mDFhydwfOzZsL2XR3mdhbWZfBu5y92ejVV81s9rcpwEfIpwp3iczawI+Fj081cx+5O5twGeAS8xsV2ASkH0GS1/bRCpCR0cHq1evZvXq1bzzzjv93s+nUB7MCHFfhg4dyqhRoxg1ahQjR47c5Db3fk8FtIrk0lHOFqk8dXVhGT68vMfNnNxaioK9vz8OOjtLN0o/dGi4ANTo0XDEEfDFLxZ3Pvf+Rqy3BkZnPd4dGA68xcbpm4yQPPvl7u3A1dGSvX450OOcqn1tEymHrq6u7oJ2oIVvvveLXfxmq62tpampiREjRgx46a1gHjlyJEN0RYkkU84WkaLIPrl19Oj+n19MHR09F92rV2+ccWbZso3zvWfmgM/cvvVWeH5P1q6F554L9x99FLbYAj75yeLF3mdh7e6fy1n1beBpd9+kCjCz5F+CTCqSu7N27do+C9f+Ctv+tq/t7X9nEZkZw4cPZ9iwYQwbNqzP+42NjXkVyg0NDWmdg1nypJwtIsXS1RXaQNat29gS0tPS1/bsbZnR77Vre7/f37YNmzW2Feaww+CQQ4q7z3yn29vF3f+evcLMdgC2BZ7t+SVSzTZs2FCSYjf7fj5XDx2sTGHbX9Gb7/My91X8SokoZ4skiHsYjc0Umxs29F209rXkvrZYBXBmKfB0mpLJ9LpnX1QodxkxYuP90aM39ldnlqlTSzcKP9Dp9j5A+Prw4OgqWtnGEs46v6PIsUmZrVu3jvb29u5l1apVmzzubX1fhXFHR0fJ425oaNisUM2nqO3tfmZ0uLGxkRrNOSQpopwt1cQ9FJnFLFJL+dq0MAu9xw0N4bavZSDPGTo0PC/TXtLb/b62peGKmwMdsV4BXE/o39s3Z9t64NpiBiX527BhAytWrGD58uU9LpltK1eu7LVI3lDs71gI/b3Dhw/Pu6Ad6HMbGxupq8v3ixeRiqecLSXR2bnpV/PZJ6ENZH0pCtwS/Ooqqbq6/ovQzFJfP/Dn5Ra3hRbEtbXJL2KTaKDzWD9lZvsDH3V3JeQS6+rqYvny5SxevHiTi2RkL8uWLdukcH777bcLPm59fT1NTU2bLCNHjtxsXe763grn4cOHU19fX4RPRETyoZxdXbq6Nr2a4UCXgRbD2evL8CXkoORTrBZaxBayr/p6XXSl0uUzj/VKehnlMLN/cfeHixZVBVu7dm33ZZizl1dffZU33niDRYsW0dbWlvdcwTU1NYwePZoxY8b0uWRmduipUG5oaCjRuxaRclPOTrauLnj77TA/8MqVG5e+Hq9aFV6TKYwz98twfnW3zEwR2fMGZy99rc/0xha7iK2rU7EqydHfPNa/A0b0sw8DdgK2KFZQabd+/Xpefvllnn/++U2WefPmsWjRogHtY/To0d0Xx8heJk6cyIQJExg3bhxjxozpLqabmprUByxS5ZSz49HVBW1tsHhxmOore9qvnpYVK0KRXMzzrgd6IldmGTas54tr9Fckp/mKeCLl0N+I9fOEXr2FbJwDNZcBhxUzqDTZsGEDc+fOZc6cOcyZM4fW1laefvrpXvuVa2trmTJlCtOnT2fatGlMmzat+/7UqVO7C2fNFSwig6CcXWTuoWCeNw8WLICFC8Py+usb77/55uBaJIYPD1fxGzUKRo7ceL+3dSNH9lwkNzZqxFYkKforrL8LtLv7ur6eZGZVc3a5u/PMM8/wxz/+kT/96U88+OCDm/U3mxkzZsxghx122GSZOXMmU6ZMobY290JoIiJFoZw9SBs2wD//CU88ES4eMW/exmUgp7CMGwcTJ8L48RuXceM2fZxZN2ZMKJJ13rVI5envAjFLBrifo6jws8wXLFjADTfcwE033cSzz246/ev222/PnnvuyR577MEee+zBbrvtRlNTU0yRiki1Us4euMWL4YEHwvK3v8FTT4W5fHsyejRst124DPLUqTBlSlgy9ydPDm0SIiIF/b1sZtOBzwFnUqFJ+sknn+Siiy7illtu6T6hcNy4cRx66KEcdNBBHHDAAUyZMiXmKEVE+lcNObsvTz0FN90Ev/sdPPPM5tu33RZ22w123hlmzgzF9HbbhQtKiIgMRN6FtZnVAocDpwMHA4sIPX0Vpb29nXPPPZerrrqKrq4u6urqOPHEE/nYxz7GwQcfrGnkRCQVqiVn96arC267DS66CFpbN65vbIR994UPfQj22ScU1KNGxReniFSGARfWZrY1cBrw78A4YAHwb8AfgPeXJLqYvPTSSxx++OE8++yz1NbWcuaZZ/LVr36VadOmxR2aiMiAVFPO7s2zz8InPwmPPRYejx0Lxx4Lxx8fimrNMCoixdZvYW1mHyWMdBwAzAO+D/wP8A13vyd62iMli7DMFixYwH777cfChQvZcccduemmm9h1113jDktEZECqLWf35uGHYdasMM/z5Mnw9a/Dpz4VRqpFREplICPWawiXwH0EOM3dXwAwsyLOwJkM69ev56ijjmLhwoV84AMf4M4772TkyJFxhyUiko+qydm9WboUjjkmFNXHHQctLWrzEJHy6Lewdvc7gTujk15ON7MZwO1Ad5OxmZl7Mae6j8ell17KE088wdZbb83vfvc7FdUikjrVlLN7c+WVYdaPD34QfvlL0AynIlIuA55S3t1fdffzCf16AFPN7FtR0j69BLGV1bJly/jud78LwDXXXMPo0aNjjkhEZPAqPWf35dZbw+3556uoFpHyyvtaTe7e4e6/dvejCH17XwAuL3pkZXbVVVfR3t7OwQcfzIEHHhh3OCIiRVGpObs3S5fC00+HeaU/8IG4oxGRalPQPNbuPg/4vJnNKVI8sejs7OSaa64B4Ktf/WrM0YiIlEal5Oy+PPhguH3f+zTrh4iUX94j1j1x9+uLsZ+43Hvvvbz22mtss802HHDAAXGHIyJSUmnP2X25775wqy8eRSQORSms0+6WW24B4NRTT6WmRh+JiEhaZQrrgw6KNw4RqU5VX0V2dXXx+9//HoAjjjgi5mhERGSwFi6EF16ApibYc8+4oxGRalT1hfXcuXNZtGgR06dPZ5dddok7HBERGaTMJcvf9z6oK+gMIhGRwan6wvrRRx8F4IMf/CBmFnM0IiIyWE88EW533z3eOESkelV9Yf34448D8L73vS/mSEREpBDz54fb7bePNw4RqV5VX1jPnTsXgN01xCEikmpvvhluJ0+ONw4RqV5VXVi7Oy+88AIA73rXu2KORkRECvHGG+F2yy3jjUNEqldVF9ZvvfUW77zzDmPHjmXs2LFxhyMiIgVYuTLcKp2LSFyqurB+8cUXAZg5c2bMkYiISKHefjvcjhgRbxwiUr2qurBeuHAhANOnT485EhERKdQ774Tb4cPjjUNEqldVF9aLFi0CYOLEiTFHIiIihejogA0boKYGhgyJOxoRqVYqrIFJkybFHImIiBSioyPc1teDLkkgInGp6sL6rbfeAjRiLSKSdpnCurY23jhEpLpVdWHd1tYGwIQJE2KORERECtHZGW51KXMRiVNVF9bt7e0AjBw5MuZIRESkEJkRaxXWIhInFdZAU1NTzJGIiEghVFiLSBJUdWH9djTpqQprEZF0U4+1iCRBVRfWGrEWEakMmR5rFdYiEqdEFdZmdqGZebTMjdY1mtmVZvYdM7vazIrWEJ0prEfoMl0iInkrd87ui3u4rUnUbzURqTaJ6UYzs+HAeODgaNWC6PbHwD3ufrOZzQJ+CpxY6PHcndWrVwMwXJfpEhHJS7lztohIGiSmsAZOB3NXOqsAACAASURBVJ4DHnb3tQBmNhk4Cfhc9Jz7gFvNbIa7v1LIwTo7O3F3amtrqdV3hyIi+Sprzh6ozMi1iEgcEvGlmZkZcALwPWCRmWVGN/YHlrj7agB37wDmA/sVesx169YBMETXvhURyUscObv/mEp9BBGR/iVixNrdHdjbzMYAXwJ+YWYrgCnA0pyntwOTCz3m+vXrARXWIiL5iiNnDzy2ch1JRGRziSisM9x9OXBeNBpyNnAvsCbnaQ3AhtzXtrW10dzc3P149uzZzJ49u9djZUasGxoaCo5bRCRfLS0ttLS0ZB6OjzOWwSpnzu6PRqxFpJQGmrMTVVhnuQJ4AFgIjMnZ1gS8kfuCCRMm0NraOuADaMRaROKUXUia2ZKYwylUyXO2iEicBpqzE9Fj3YMu4O/A/cBkM2sEMLN6YDrw50IPoBFrEZGiKXnOHii1gohInBJRWJvZZDM7ycxqoq8UvwSc6+5vAXcCB0ZPPQi41d0XFnpMjViLiAxOHDm7/5hKfQQRkf4lpRVkDHAx8A3gYeDyrKmZPgNcYma7ApOAwTfhZckU1hqxFhHJW9lz9kBpxFpE4pSIwtrdnyF8XdjTtuWE+VKLStPtiYgMThw5uz8asRaRJEhEK0gcNGItIlJ5NGItInGq+sK6vr4+5khERKRQGrEWkSSo2sK6q6sLQJczFxEREZGiqNrCurOzE4Camqr9CEREKo5aQUQkTlVbVWYKa41Yi4ikn1pBRCQJVFirsBYRqRgasRaROFVtYa0eaxGRyqERaxFJgqotrNVjLSJSeTRiLSJxqtqqUq0gIiKVQyPWIpIEVVtYqxVERERERIqpagtrtYKIiFQetYKISJyqtqpUK4iISOVQK4iIJIEKaxXWIiIVQyPWIhKnqi2s1WMtIlI5NGItIklQtYW1eqxFREREpJiqtqpUK4iISOVRK4iIxEmFtQprEZHUUyuIiCRB1RbWmR5rtYKIiFQOjViLSJyqtqr0KPuqsBYRST+NWItIElRtVZkZsTZlYxGRiqERaxGJU9UW1hqxFhGpHBojEZEkqNqqUiPWIiIiIlJMVVtYa8RaRKTyqBVEROJUtVWlRqxFRCqHUrmIJEFd3AHEZebMmRx++OG8+93vjjsUEREpUH09HH44DB8edyQiUs2qtrA+5phjOOaYY+IOQ0REimDECLj99rijEJFqV7WtICIiIiIixaTCWkRERESkCFRYi4iIiIgUgQprEREREZEiqOrCuqWlJe4QikbvJZn0XpKpkt5LtUnjzy6NMUM641bM5ZPGuMsRswrrCqH3kkx6L8lUSe+l2qTxZ5fGmCGdcSvm8klj3CqsRURERERSwrxCrv9qZm3AgjxfNh5YUoJw4qD3kkx6L8mUtPeylbtPiDuIchpkzobk/ewGIo0xQzrjVszlk8a4ixVzrzm7YgprEREREZE4qRVERERERKQIVFiLiIiIiBRBVRbWZtZoZlea2XfM7GozGxl3TINlZkeb2XwzW2pmPzSzurhjKpSZDTGzuWa2f9yxFMrM9jazL5vZkWY2Pu54BsvMdjSzq8zsC2b2YzN7b9wx5cPMDjazx8xsRta6iskD1SCtPy8zu9DMPFrmxh1Pb9L6f6SnuKP1ifzce/udbWbjzeynZvZdM/u+mdXHHWtGX3WGmd2Y9TnfHmecucxsPzN73syWm9nlWetL+u869UXYIP0YuMfdbzazWcBPgRNjjilvZjYdOBI4FtgJuBp4DbgszriK4CvAjLiDKJSZ/TuwrbufF3csRXADcIS7L4z+3f0B2DHmmAbEzCYCI4C9cjZVRB6oIqn7eZnZcMLJUgdHqwZzsmbJpfX/SG9xJ/Vz7+d39i3A19z9MTM7A7gA+HpcsWb0FXO0bTEbP+fnYwmyB2Y2CtgX2BvYB7jNzO509/so8b/rqjt50cwmA/OBMe6+OvrLayWwk7u/EmtweTKzDwB/dfeO6PElhPdxaLyRDZ6Z7QO8C/gmcKq7PxBvRINjZvsB5wMHewX8JzOzd4A93P05M5sAzHX3yXHHNVBmVgN0Alu7+yuVlAeqQVp/XmZ2dnT3J+6+NtZg+pHW/yO5cUfrEvm59/Y7G/gW8JtMTjWzKYQidaK7vxNXvFEsvdYZZvYD4C/AbZntSWFmjcDazO9fM5sDfJnwuZb033U1toLsDyxx99UA0T+G+cB+cQY1GO7+UM4/5oXAq3HFU6holOFYd78u7liK4PvAP4ErzeweM9s37oAK9GvgWjNrAj4GnBVzPHlx966cVftTIXmgSuxPyn5eZmbACcD3gEVmlpiR3p6k9f9IbtxJ/tz7+J19AOGzzTxvYXS3uYzh9ai3mM1sGPBvwP9Gjw+MJcBeuPuarKJ6OOH38QOU4d91NRbWU4ClOevagdSMvvVhT+AncQdRgHOAi+MOolBmtgOwG3CNu58J/Am4J/raMq0+C6wF/ga87e63xhxPoSo5D1Si1P28PNib0JJwJfALM/tIzGHlI3WfOaTuc8/8zk7TZ70n4ZuA1e6+HSH23wN3m9ku8Ya2uWjE/W6gHhhGGT7raiysHViTs64B2BBDLEVjZtsBbe6emJM08mFmHwZa3X1x3LEUwU7AMnd/Mnr8I8CAI+ILqWCNwE3AL4HLzeygmOMpVEXmgQqW2p+Xuy+PzrO4CDi7v+cnSGo/c0j+557zOzsVn3VPdYa7v+nupxN6xM+MLbjezQN+DhwIXEoZPutqLKwXAmNy1jUBb8QQS1GYWS3waeBrccdSgC8D15nZEjNbAkwDbjezr8Yc12DUkXVisLuvAV4ExsUWUeFuBG519wsI3yrcFH29llYVlwcqXCX8vK4ApscdRB4q4TOHBH7uPfzO3uSzjtpZhpOgz3oAdUbiPmfoLvyvI9QY+1GGf9fVWFjfD0yOGtuJprSZDvw51qgK82XgMndfF3cgBTgJeG/W8gZwGulsbXkSGJXT+tEBPBtTPAWxME3gru6+Ilp1IbCKlMwK0otKzAOVrBJ+Xl3A3+MOIg+V8JlDMj/33N/Z9wI7ZG2fAawHWsscV1/6qzOS+DlnayUU1SX/d111hbW7vwXcSfhaAOAgwkjcwt5flVxm9g3gCaDRzLYxs09GX9ekiru3ufvrmYVwlnebu6+KO7Z8uftzhJ6uYwHMbDQwhNCHlkbLgLXRmeoZS4EXYoonb9EIEISWnIrLA5UujT8vM5tsZieZWU307+9LwLlxx9WbtP4fyY076Z97T7+zgRXAm2a2c/S0WcCV7r4+rjiz9VJn7Gxmh0fbhwCnAJfEGWc2MxtmZrtlrToE+GE5/l1X3XR7AGY2hvAP4BVgEvB1d3871qAGwczOI0zTk+05d0/zSCIAZvYK6Z5ubzzwQ8JfydMJJzKmcsQawMx2JZzAOAeYCDzo7qkYucqayeTHhLlhf+TubZWSB6pF2n5eZrYT4Q/sVcDDwOXRH92Jk9b/Iz3FDWxBQj/3vn5nm9lWhGlmXyHMzX1uEqaw6y1mQj/1TYQZNR4FvuvuSWpd2YswmPUcIb7H3f2WaFtJ/11XZWEtIiIiIlJsVdcKIiIiIiJSCiqsRURERESKQIW1iIiIiEgRqLAWERERESkCFdYiIiIiIkWgwlpEREREpAhUWIv0wswONLNHzGz/fp63s5n9r5mdX8RjX1qsfYmIxMHM3mNmtxQzN/ZynDvN7IRSHqNUzOwgM/trf79nJD3q4g5ApD9mtiPwbcIVAPcE3gN83N1vLPGhXwN2HcDzXgXGALXFOGh0idUNxdiXiAiAme0JXAR8EPgj4cIek4B64Gp3/0MJDjsPGEWRcmMffgw8XYodR/n4aOAnhKvN/jXa1AAcA+zl7q8UcIiB/p6RlFBhLWlwG3Ccu881sxrgZ+U4qLu/YGZLB/C8VWb2ZhEPvR/wYBH3JyJVzt3/Zma/At7l7odl1pvZocBvzezL7v7jIh9zdZFzIwBm9m/AU5mC1t3vKvYxMtx9A3CzmZ0FPODu52XFcftg9mlmWwAfcveb3f15M1tSpHAlAdQKIolmZmOB7YF3ANy9i3DZ13JdMrSryM8biA8B9xdxfyIiAJtdItvdfw/8J/ADM5tUgmMWNVeb2VTCCHW59XR58fuAhfnsxMwagOuBxqzVxfz9ITFTYS2J5u7LgOeBe8xs72jdAuBxMxttZt8ys8fN7HAze8vM/mlmu2deb2Ynmtm5ZvaAmX0na32zmf2nmd1sZr8xs1HR+qFmdpWZfdPM/pvwNWZezOyrUUzHmdmeZnZbtL9vRTHebGYTovUrzew/cnYxxN3XRfv6kpl90cz+YmYX5/0Bioj0778JLSEfBTCzbc3sAjO71sz+FBWzmFlDlMvONbOHzOz4zA7M7LNmdr6ZPWpmZ/Z2oJ5yspkdG+XMU83sRTO72Mwmm9nPonz6kJkdEO3iJGAqcI6ZHRmdC3OfmZ2SdYzTzezbZna9mf2PmQ03s1Fm9h9m1mpmB5jZ82b2spnNGMwHFn0m+0cj2pjZ0Wb2QzP7npndY2YjovfwQzO7I/p982b0Gb8b+KiZfTFrl2Oi30XtZnZ6L8esMbPvmtns6BhuZr/N/HwkIdxdi5ZEL8BOwALCX/XXAVtE62sJSWoZoddtNHAv8Azhj8b3A+dFzx1DGPU+FBgJ/Cpr/38BLo3u/wCYHd2fCHQSkmd/Mf6cMOqzBaEXb1i0vgH4A3A7MJnwC2ED8B/AMGCP6HHm+Zn+8cz7/nF0vwn4Wtw/Cy1atKR3AU4FXu9l2yLgiiiv/gaojdb/Ergluv994Ojo/ieBtuj+icDHovvbR7l6l+jxz4H/jO73lpPHAiuAy4G9gIOAy4BvRs/9CnBHVqwOzIjuNwL/BE6NHn8Y+GXWc+8EfkpofT0q+n1xdPQ+7wG+PcDP7oHod8XFwKWEfuuDsrYvAaZH958h/E5qAM4FXgfeB3wqyvsPZOKNnv8K8L1o24nAq73EcBpwZ3R/YvQ5/2vc/660bLqox1oSz92fMbNdgAuBzwL/ZmYHeei5bgNWufutAGb2n8DDwDaEXyJdZnZ2tKs/EBL4ocC4rPX/jF47AvhM9FrcfZGZvZZHqOOBK4FPufvqaB/rzOwtYL67vxEd5y3gzx76D58gJPzxhJMgDyX8EgBYDXzCzF4Arib8USEiUgpOKDb3ArYCzjIzgKXAKAsPPk0YFAD4H0IrBIRc25qTa7cEnso5xqn0kJPdfZmZrQRud/fHAczsVeAdM9uSMMI7oseg3deY2eKsVWdG+824FvjfaP0Kwu+L30THaI3iHKj7PeqxNrObCYM5GR9191fNbF9CQT0iyv9vAC+5+2PAY9Fre9r3HdHvhL/1EdPuwJrofS8ys2Vs2lIiCaDCWhLNzGqBaR5OUvl/UXvGXUALYQQg1zPR7WhgOnCDu/8yWnd5tM9zCIXu5TnH2oOQENsHGe72wL8A2wJzs9bn9hh2dm9w74qSbKYta5yH9hfcfb6ZnQr8CPgCcDKwGBGRIjKz0YRv254h5M2lPeTHiYQR1SEA7t5JGAwges1/ufsj0eNNXpulx5wccTbNla8D5wAvAo8CW/fxFrJfNxP4v6zHLxHaXLZg81zcweBnLPk74Q+QjBVm9l3gFkLsmeo59331p4Pea7O7gavNbCjhZ7GGMPotCaIea0m67r4/AHefC3wDeFcvzx9CSGLzgDeBY7M3mtn7ovWHmlljzvpV0cMdBxnrI8APgV9EiS8vZjaZrBNhzGwKYUaU7QknM/6yl5eKiBTiFMI3ZLcQ8uO+0Ugx0J0flwDrgVlZ63cxsyZycq2Z1UUDFbl6y8k9+QnwiodpVdfl8V5eJYxwdx8CWElodSkad++KBj92iv4wuQu4wN1bi3mcnGPeAdxB+Gb1RGA/d19ZquPJ4KiwljQ4Kyo6M6YQeuMyxprZkOj+/oT+uhXATcBRZvYjM9s/Gk14h5AAhwN3mtmHoxNF3kUoxp8AvmVmQ8xsPKG3eQszqzWzr5nZx3qJ0aLlm4RfPt/vYVvu83MfH0boB8yYCRzp7quAs3p4jYhIPjYbCTWzDwL/BXza3dsIo8OLgLvM7N+iExQ/EI1Q/4owe8iRZnYgcIq7txNy7f8zs/8ys/0IAwyZNrrs/NdbTs7IHj3eDZhgZsOAfYFGM9sm2raBcLJfZoAl+xg/AY4xszHR4/cB17h7B6Hm6TEXm9neUVwNvXx2NeTUTNEfH4cSvqUcQ2gx3AGYBgzLOjEyd1R8fQ/xbxpUD/0i0c9iJPAnwqh8Pn9wSJmoFUTSYBrwjJn9nvA1mQFnZG3vAM41sxWE0d2zANz9XjP7f4SvE48knPz3NICZHQFcRfhFcR3wFXd3MzuOMBXSM8DNhK/03ku4oMJxhBaTTS5MY2Z7EU7KmU74avA24L+ieP4A7A3MMLNtCSMpk4DjzOw5No7enAoMd/eWnPf+MwsXduggnCwkIpK3KI+cRCj+rgHaCCfAjQI+4u5/BXD39VF+vIaQ624n9FYDfJ4wKHE94Ru6zEwc/03IfZ8htKx9xt0XR4Xj3sC2ZnZjbznZzI4l9BWfZmbPuftCwomUlxGK6uuBI4BdgJeBGwij65+wMFvUTsAhZnaXu//GzGYCt5rZXwjfep4bjax/FJgUvb9nCBfL2dLMdoriPIHQN35b1udWT8j9uwJTosLbCb3NR0f7/AfwZ+BvhPaWuwk5/U7CCZO7mtmR7p7Z782Ei/W8HRXik4BjzWw+YSSa6PW512x4lXBy5pHAUKDGzO4jnMBYrilopR+mn4WkmYXLwP7c3WeU6XhfcXddblxEpMJEbSkN7p7IC3SZ2VcIM7S8Ej0eSvi24buZc3MkfmoFkbTrqc2iNAcKX5n+vRzHEhGR8olOztw9wUX1UOA8YEbW6vHAOyqqk0WtIJJaUX/bxwhf7X3U3W8p8SFfir6iFBGRylLn7lfHHURv3H2thcuqXxGd2P4GoW3lvL5fKeWmVhARERERkSJQK4iIiIiISBGosBYRERERKQIV1iIiIiIiRVAxJy+OHz/eZ8yYkffr2tramDBhQvEDqnL6XEtDn2tpxP25zpkzZ4m7V9UPVjk7f3rveu/VJqnvva+cXTGF9YwZM2htzf9Kos3NzYN6nfRNn2tp6HMtjbg/VzNbENvBY6KcnT+9d733apPU995XzlYriIiIiIhIEaiwFhEREREpgqovrGfPnh13CBVJn2tp6HMtDX2u6VHNPyu99+qk954uZb1AjJkdDXwPGAncCHzJ3TvMbDzwbWAFUA+c4+4botf0ui1bc3OzJ7EPR0SkP2Y2x92b444jl3K2iMjm+srZZRuxNrPpwJHAscAXgNOAs6PNtwDXufs5wAvABVkv7WubiIiUgHK2iEj+yjkryFbAJ929A5hjZjsDHzKzh4Ed3P2x6Hl3AM+b2YXALr1tc/d3yhi7iEi1Uc4WEclT2Uas3f2hKEFnLAReBQ4A5mc9b2F0t7mfbYW5/nrYfXf4wQ8K3pWISKVJWs5u/0c7rbu38ty/P1forkRESibOkxf3BH4CTAGW5mxrByb3s20TbW1tNDc3dy8tLS19H33xYnjiCXj99UGGLyIyeC0tLd35ChgfdzwDEGvO7nqni7efeJvVz68ebPwiIoM20JwdywVizGw7oM3d55qZA2tyntIAbAD62raJCRMm5DeJuFk+IYuIFNXs2bO7z3g3syUxh9OnZOTsfCIWESmugebsso9Ym1kt8Gnga9GqhcCYrO0GDAfe6GdbcZRxVhQRkbRJXs4u2p5ERIoujlaQLwOXufu66PG9wA5Z22cA64HWfrYVJjNircJaRKQvCcnZBe9BRKTkylpYm9k3gCeARjPbxsw+SZjr9M3ojHOAWcCV7r7e3Vt721aEYArehYhIJUtUzo6U89oLIiL5KluPtZmdB3wrZ/Vz7n6dmR0PfNPMXgFGAOdmPaevbYVTkhYR2UzicrbGQkQkBcpWWLv7hcCFvWxbAHwq320FUSuIiEivEpezuw9Qsj2LiBQszun24qVWEBGR1LDuwZB44xAR6YsKa41Yi4gkn8ZCRCQFqrewzlBhLSKSHkrZIpJg1VtYqxVERCQ9lLJFJAVUWGvEWkQkNTTdnogkWfUW1hlK0iIiyacRaxFJgeotrNUKIiKSPhoLEZEEU2GtEWsRkcTTdHsikgbVW1hnqLAWEUk+fckoIilQvYW1WkFERNJHYyEikmAqrDViLSKSfBoLEZEUqN7COkOFtYhIami6PRFJsuotrNUKIiKSHkrZIpICKqw1+iEikh5K2SKSYNVbWGeosBYRSTxNtyciaVC9hbVaQURE0kMpW0RSQIW1RqxFRNJDKVtEEqx6C+sMFdYiIsmnEWsRSYHqLazVCiIikjqabk9EkkyFtZK0iEjyaSxERFJAhbUKaxGR9FDKFpEEq97CWkREUkPT7YlIGlRvYa0RaxGR9FAriIikgAprFdYiIumhlC0iCVa9hbWIiKSHRqxFJAWqt7DWiLWISOpouj0RSTIV1krSIiLJpxFrEUmB6i2sRUQkfTQWIiIJVr2FtUasRURSQ9PtiUgaqLBWYS0iknxqBRGRFKjewlpERNJHYyEikmDVW1hrxFpEJD00Yi0iKaDCWoW1iEhqaLo9EUmy6i2sRUQkPTIj1qqrRSTBqrew1oi1iEhqdM8KIiKSYCqsVViLiKSHUraIJJgKaxERST6lbBFJgeotrDM0Yi0ikh5K2SKSYNVbWKsVREQkPTRiLSIpUPbC2swONrPHzGxGzvoLzcyjZW7W+kYzu9LMvmNmV5vZyCIFUpTdiIhUssTk7Iim2xORJKsr58HMbCIwAtgrZ/1wYDxwcLRqQdbmHwP3uPvNZjYL+ClwYsHB1ER/U3R2FrwrEZFKlKicren2RCQFylpYu/siM7u9h02nA88BD7v72sxKM5sMnAR8Llp1H3Crmc1w91cKCmbIkHC7fn1BuxERqVRJytmabk9E0qDsrSDu3pX92EK2PAH4HrDIzLJHNvYHlrj76ui1HcB8YL+CA1FhLSLSr8Tk7O6AirYnEZGii/3kRQ/2JnyteCXwCzP7SLR5CrA05yXtwOTc/bS1tdHc3Ny9tLS09H3ghoZwu25dQfGLiAxGS0tLd74i5L9UiC1na8BaRGI00Jxd1laQvrj7cuC8aDTkbOBuwtjEmpynNgAbcl8/YcIEWltbB37AzIi1CmsRicHs2bOZPXs2AGa2JOZw8lb2nN194PxfIiJSqIHm7NhHrHtwBTA9ur8QGJOzvQl4o+CjaMRaRKQYypOzNWItIimQxMK6C/h7dP9+YLKZNQKYWT0hgf+54KOosBYRKYby5OyIptsTkSSLYx7rzLiDRY8nm9lJZlYTbfsScC6Au78F3AkcGL3mIOBWd19YcCA6eVFEpF+Jydmabk9EUqCshbWZNQFnRA9PNbMJhK8NLwaeBK4Gfp4zLdNngCPM7FzgEGB2UYLRiLWISJ+SlLM13Z6IpEFRTl40s0Pc/a7+nufu7YREfHXW6jY29uf19JrlhDlTi2v48HD79ttF37WISJKlMmdn6uquPp8lIhKrvArr6Cpa5wGT2DjaXQNMBBqLG1qJjRoVbleuhK6ujVdiFBGpEJWUs60+VNbeoV4QEUmufEesbwAuAp5g47iBEa60lS51ddDUBO3tYckU2iIilaNicnbNkPB3Qdc6DVmLSHLlW1i/7O7fz11pZs8VKZ7yGj06FNUrVqiwFpFKVDE52xrCiHXXehXWIpJc+fY/fMbMju1h/eHFCKbsRo8OtytWxBuHiEhpVEzOzoxY+zrXlHsiklj5Ftb7AteaWWfW0sWmJ7akx7hx4XbRonjjEBEpjYrJ2VZjWF3UZ71BhbWIJFO+rSAXEKZOagU6o3UGfLyYQZXNVluF2wUL4o1DRKQ0Kipn2xDDO5yu9V3dI9giIkmSb2F9N3C7u28y+bOZXVO8kMpoxoxwq8JaRCpTReXsmoYaulZ34escRsQdjYjI5vItrH8FfN3M/i9rnQFHA58vWlTlkims58+PNQwRkRKpqJxtQ3QCo4gkW76F9fnATOBUNl5Y1ghzpKYuSfOud4XbuXPjjUNEpDQqKmfXNERT7q1VYS0iyZRvYX0pcHd0Na5uZnZI8UIqo/e+N8xn/eyzYdq9pqa4IxIRKaaKytl1o+tY9+o6OlZ0xB2KiEiP8j3747XcBB2ZVIxgym7o0FBcu8Mjj8QdjYhIsVVUzq4bE8aCOparsBaRZMq3sP5s7goz+yhweXHCicGHPxxu77gj3jhERIqvonJ2/Zh6QIW1iCRXvoX1WDM7H8DMhpnZtcC1wN+KHlm5HB5dJ+G226Czs+/nioikS0Xl7LqxYcR6w7INMUciItKzvAprdz8UuNfMrgT+AWwP7Ar0dGWvdGhuhm23hYULNWotIhWl0nJ2/bgwYr2hTYW1iCRTnycvmtkHe9m0AHgRuAzYCjgS+EJxQyuTmho46yw4+2z49rfDCHaNLjwgIulT6Tl76IyhAKydvzbmSEREetbfrCCXADsCK4Ce5je6jjB105akMEl3O+00uOQSaG2FG26AU06JOyIRkcGo6Jw9dOtQWK+ZvybmSEREetZfYX0B0Orui/t6kpml8mvFbsOHw0UXhYL685+H/fbbePEYEZH0qOic3bhNIwBrXlRhLSLJ1GfPg7vf1V+Cjp736+KFFJOPfzy0gaxcCcccA6tWxR2RiEheKj1nD912KDXDalj36jrWt62POxwRkc2omTjDDH72M9hmG/j73+Gww0KRLSIiiVBTV0PTHuFCXu2P9zQ9t4hIvFRYZxs7Fu67D6ZMgYcegn32gZdfjjsqERGJjNxnJADL/7Q85khERDanwjrX1lvDww/DTjuFS52/971w3XXh6owiIhKr8YePB2DJb5fgyssikjAFF9ZmtpuZ3WhmXzKzGYWHlAAzZoRLnB99NLS3w6c+BQccAHPmG40dvQAAIABJREFUxB2ZiEhB0p6zR+49kiGTh7D2lbWseGBF3OGIiGyiGCPWZwAHufv3gH3MrLEI+4zfyJHw61/DjTeGFpEHHggXkzn+eBXYIpJmqc7ZVmNM/vRkAF677LWYoxER2VQxCuszgT2i+78C6ouwz2Qwg5NPhpdegi9/GYYMgf/931Bgf+hDcPPNsFYXKhCRVEl9zp78mcnUDKth2V3LWHbfsrjDERHplndhbWYfMbOTo/t7ATPdfSGAu3e5e+XNUzd6NFx6KcybFwrspqYwgn3CCTBpEpxxRnjc0RF3pCIim6jEnD1kwhC2On8rAF749AtsWKFLnItIMuRVWJvZt4FbgKMA3P1x4BNmdlIJYkueadNCgf3aa3DFFbDHHmFKvp/+NIxgb7FFmA/7lltgmUZRRCRelZyzp31hGiN2G8Hal9fyz5P/SdeGni40KSJSXvmOWO8DTAUeylp3N3Bx0SJKg1Gj4KyzwiXQn3oKzjkHdtgBli8PPdnHHQfjx4cZRc4+G267DZYujTtqEak+FZuzaxpq2OnWnagbU8eyu5bx7AnP0rVOxbWIxCvfwvped18BZM9xlMpL4xbNzjvDxRfDc8+F5dJLYf/9Qz/23Lnwwx/CUUeFQnvmTDjpJLj88jDryBpdlldESqqic3bj1o2854/voXZULUt+s4R/7P8P1i1cF3dYIlLFLJ95QM3sBGAbYAphBOSjhK8Yz3L3q0oS4QA1Nzd7a2trnCFsau1aePTR0Hv9wAPw2GObn+hYWxsK8/e8J9xmlmnTwomTIlIVzGyOuzeXYL9VkbPbn2jn6SOeZt1r66gbU8e239+WSadMwpRHRaQE+srZeRXW0c72Bk4FpgNvAje5+32FBlmoxBXWuTZsgGeegccfD8vf/gZPPw1dPXx1+f/bu/M4yar67uOfX1X1Nr3M9PRsDDDOAKODSJSIqA/xYURQ3KNx1/hoHhnj8qhRsxmMcYlbIC9UYnSMvNTEBaNGDCpuISpBgxhfICgwwAwDAww9a/dM79W/549za/p2Ta3dt7q27/v1Oq+quvfWvadO1fzm1+eee25/f0iwt2yB004LZfPm8Njfv/R1F5GaqlViHe27LWL21PAUt7/6dg5cG65vGXjyAJs+uInBrYOJHUNEBBJOrAvsPA1sdfcfLWpHi9TwiXUhR4+G4SK33RbGat96a3jct6/4e9aunUu0N22CDRtCOfnkULq7l67+IpKIWibWBY7VsjHb3dn7L3u5++13M70vzBSy/LzlnPSWkxh63hCpjG42LCKLt+DE2sy+BfSV2f8qYJ+7n7/wKi5eUybWxTz8cEiyd+wI5a67wuPdd8NkmfGDa9aEBDuXbK9fH8oJJ8w9rlihoSYiDSSpxFoxO5gZneH+y+/nvkvvIzuSBaDr5C7Wvnota162hr7HlGsiEZHiSsXsTJn33hJtcwcwC1wM/AqIR8NNwMoE6ik5a9aEW6ifn/f/3uws7Nkzl3Dfey/s3h2m/9u9G+6/PyTlDz9c+u6Q3d0hwc6VXMIdf75+fbjjpBJwkWaimA1k+jNsfPdGTnrLSTz0+YfYc8UexneMs/tvd7P7b3ez7IxlrH7halZetJKBJw5gacU5EUlGuR7r5dE2h6LX73P3vy6w3dXu/vzaVbO8luqxXqhsFvbuDUl2LuF+8EF44IH5j6Ojle2vszPcAGf9+vC4ahWsXh1K7nl82bJltf18Ii0qwR5rxewCfNY59ONDPPyVhxn+2jAzB+Zu5pUZzDD49EEGzx9k+bnLWXb6MiylRFtEiltwj7W7H85btNbM0u6eje18K3DOomspi5dOzw39eNKTim935EhIsHMlP/HOPT98eC5Jr0RPz/HJ9uBgKCtXFn4+OKhx4SIJUcwuzFLG4FMHGXzqIJuv2MzBHx0Mt0P/7gHG7xpn+Kphhq8aBkKivfzc5QycO0D/2f30n9VPx1DT3fVdROqk3FCQfF8Dbjaz64FR4HTgGcBxPSLSwPr6wsWPmzeX3m5sDB56KCTZe/fC8HC4sHJ4+Pjnw8NhXu5qEvGcnp7CCXf+64GBUJYvn3s+MBB61kWkEMXsPKmOFEMXDTF00RAAY3eNcfD7Bzn808Mc+ukhpvZMsf+a/ey/Zu6mXl0nd9F3Vl8oj+uj99G9dJ/SrYshReQ4C5lubz3wRuB3gBHgX9396hrUrSoaClJn7qEnPD/xPngw3N794MHjn+dez8yU338p3d3zE+38xLvY6/7++aWnR2PKpS5qPN2eYnaF3J3J3ZMcvv4wh284zJFfHeHIzUeYHTt+WlTrMHo297Bsy7K58qhldJ/STcdQh+bQFmlhtZ5uLwWc5e4lrparvUYM0lIB9zDtYLGkO/56dBRGRsIQlZGRuefZbPnjVCKVmp9o9/Udn3yXK/H3dHUpUZeKLPF0e4rZVfCsM7ZjLCTZvzrCkVuOMHb7GJP3Fp+hKd2fpntTN92buuk5pWfu+aYeuh7RRaav2pPFItJIFjMrSP6OriyweDWQBX5/AXWTdmcWktG+vjA9YLXcwxCUXKKdn3gXep1bNjo6v0xMhOWH84epLlAmU30yXqp0aJynVEcxe/EsbfRu6aV3Sy9rX7722PLs0Sxjd44xdvsYY7+NHu8cY2LnBNmRLEdvOcrRW44W3Gd6eZquk7pCObFr7nnsdWYwo15vkSZU7Z/Nm4EfAfFu7pMIt8qtiJldCHwAeKm774qW9QAfJYwBHAT+3N1Hyq0TwSzMRrJsWZi5ZDGmp8NwlvyEe3S0+PJSZXp6rrc9CV1d85PxXEK+0EcNfWkHitk1ku5N039WuLgxzt2ZOTDD+M5xJnZOMLFzgvF75p5P3jdJ9nCWscNjjN02VnT/qZ4UnWs76VjbQefazmPl2Ot1c8vSA2kl4SINotrE+nXufkd8gZltIdwutywzW0u4eUH+FemfBK5196vM7CLg08DLK1gnkpyOjrkLJZMwNVV9Ml4qgZ+cDKXUnTmrkTtbsNDkvNAy9ao3GsXsJWZmdAx10DHUwcDZA8etzyXek/dPMrlnMjzmSux1diTLxK4JJnZNlD9ml80l3kMddKzqIDOUoWNVx7HX8cfMUIZ0d7oWH1+k7SUxxnoT8At3X1Xh9inCachN7r4rurBmJzDo7mNmlgEOA2cAU8XW5XpOcpplvJ7IgriHoSrx5LvUYyXbjI8nX8+uroUl5StWzP1Rs2JFKJn2GYe6xGOsFbObwMyRGab3TjP10BRTe+dKoWWzR4+/uLKcVG/q+MR7qIPMygyZFRkyg+GxY7Dj2PPMYIZ0n3rHRZIcY72T+acU08Ba4J8r3Ye7z+b9o9xKuL3uWLR+JjrOecB0iXW7qqm7SFMzC0M3enrCnTmTMDMTLhxNIknP71Xfv7/88csplHAXel6odHUt/vgtQDG7eWX6MmT6MvSc2lN22+zR7FzivX+a6X3TzOyfYXrfdHgdLYuvmz06y+TRyZIXYRaUZi7hziXgscR73vIVUVmeIT2QJjOQIbUspcRcWlq1XUI/AL7EXKCeBfa6+52LqMOJQP7/wqPAemCmxLp5hoeHOfvsuT8etm3bxrZt2xZRLZEWl8mEqQeXL09mf7le9WqT85EROHQojEWPP+aS9fvuq74uXV3zp1gsVfKnY8yVwcGazpG+fft2tm/fnntZUe/xAihmt4F0b5qeU3roOaV8Eg5hOEp2JDs/4c4l3QdnmDk0M+9x+uD0seezY7PM7J9hZv8Cp0lNQ2ZgLtHOPcaT7/RA+vjXA0rQpb4qjdnVJtb/4u4/iS8ws0eZ2WnufleV+8pxIP+cdBeh54My645ZvXo1Oq0oUkfxXvXVqxe3r9nZkFTHp1vMJd35zwuVycm5+dQXo68PhoaKlxNOgE2bQlm5sqqLQeOJpJklNHD+OIrZchwzI7M8JKqVJuM5s1Oz8xLveNKdn5Tn1mVHssyMhMfZ8dlj205SZW95XBoy/RnS/WnSfbGS/7rIskx/5rhlqS7d8EeKqzRmV5tYXwT8JG/ZXuA64Kwq95Wzh3DleFw/8AAhgBdbJyKtKpWa603fuLG69+Z6zvOnWKymHD4c5lHP9bbfe29lx37DG+DjH2+k8eGK2ZKoVGeKzjWddK5Z2Nmc2enZeYn2zMgMM4fnJ9/x18XWzU5ECf6hRd5gLMY6rKKkPN2fJt1bYHnf8ctTPepZbzcVXbxoZm8mTLfUX2STH7j7RRUf1MyZuxBmHXAXsNrdx82sAzgCnEK4YKbgOnffE9+nLoQRkYplszA2FsaYFytHjsCDD4ahKPfcA3ffDbt3l9/3jh1w2mlVVSfpixcVs6XVzU7NMjMSxopnj2SZGZ0heyQ7v4xmK142MzoTfr1JMwom4aneVNmkvNTyVKd61+tp0RcvuvsVZvYT4P8AH89bPeXuD1ZRmdyfbhbt+yEzuwZ4GnANcAHw9VwQLrVORJaQe0hIZ2ZCmZ6u7vlC3lPpfqemwvCPiYm5Cyhzz/MfJyfD+5PS0wMnnhiGgmzbBqeemty+F0gxW1pdqjNF56rOxK5OcHd8yo8l2yUT9dEs2aPR86PHJ+q5ZD97JPSs554nyTLH965XnKwXWZbqTZHKKGFfrIrPV7r7LWb2F+5+3P9IZrbB3ct25ZhZP/Cq6OVrzOwKdx8G3gB81MweC6wD4lewlFonUl/ujZFo1jJpjb9uJb294cZCvb3FS2597mLG3GwkucfVq8NjA57qVcwWqZyZYV1GqitFx1By8/HPzsyGRLtAAl7psniinis+44kPhQFIdaeSS9b70mQGM6Q62itZLzkUxMzeCXzH3X8Tvf4rwnRN8zYDnuruW2tVyUrotGITmZ0Nydr0dOhpLFYWun4pk9ZsLc4dNrBUKtwEJpOZK/HXST1fyHu6u0Pp6gol97zQss7OhkqGkxoKopgt0h5mp2YTS9Lj21F+dHDVcgl2ZmU0HePKaGrGlR3zHo+tH8qEO4r2NO5NjBYzFGQTsCL2+neBXuAh5prfCL0S0oiy2XAKPH46vFgptz7JxLeVxJPBeieXtdxvOh0Sa2lkitkibSDVmSK1MkXHyuR6192d2fHZxBL1mZHQo55bPnlfdbPApPvT4Y6i6zrpWNtx7O6iuWVdJ3fRvbGbzGCmoS4QLZlYu/ub8hb9LXCru0/FF5rZo5OuWMsrdPFUuYupCq0fHy+dGDdqEtvREXoN80ux5dWsX8qkNZ1uqJ5PaW+K2SKyUGZGelma9LLkeordnexolukD0RzpB6JpGA9Ec6TnlseeTw9PM7V3iuxolvHRccbvKn2X4HR/mu6N3XQ/opvujd30PLKH3jN76TuzL9FhPZWqdk6oM939f+ILzOxRwKnAbxKrVbPIZsM8uXv3hnLgQJim69ChUOLP81+PjS1NHc3mTpHHS+60eLGSvz73OqnEV8moyFJQzBaRujGzcHOfgQxsrPx97mEM+dTeKab3hkR76qGpY3cYnXpoisndk0zsnCA7muXor49y9NdHj9tP5wmd9D++nxXnr2DFU1fQ99i+mvduV5RYm9lTCKcPL4xuTxu3knDV+b8nXLf6O3gQ7rorlLvvDuX+++cS6eHhcPHaQhW6cKrcxVT52+VuyFEsMVYSK9J22jZmi0hLMDM6BjvoGOyALcW3c3dmDs4wsWsilJ0THP1tSLKP3nqUqQen2H/NfvZfE24I23NaD+v+7zpOfPOJZPpqc7+BSvd6CPgCYfzeuXnrpoB/SrJSdbNzJ1x9NfzsZ3DjjbBrV/n3rFoFa9eGMjQUZgtYvjw8lnre26vxqiJSK+0Rs0WkrZkZHSs76FjZQf/vzp+232ed8XvGGblhhEPXHeLA9w4wftc4O/9yJ3s+sYdHX/VoVvzeiiJ7XrhK57H+tZltBV7s7q0XkK+9Ft77Xvj5z+cv7+mBRz4yzEt76qnhpg8bNswl0qtXhx5hEZEG0vIxW0SkDEsZy05bxrLTlrHu1evwrHPg+wfY9Z5djP5ilFuefguP/+Xj6T29N9HjVjOP9WGK9HKY2e+5+/WJ1WqpuMPb3w6XXx5e9/XBc58LF1wA55wDp58eLk4TEWkyLRmzRUQWyNLG0DOHWPn0lfzmlb9h+Kphdr5rJ4/5t8ckepySibWZfQvoK7MPA84A1iRVqSVz5ZUhqe7shPe/H9785jB2WUSkCbV8zBYRWSRLG6ddfhrDXx1m/7f3k53Iku5OrhO1XI/1HYSxensoPm24Ac9JrEZLxR0+8pHw/DOfgVe/ur71ERFZvNaN2SIiCela10XP5h7G7xxn/I5x+h5brj+icuUS648Ao+5eclZvM2u+q8v37IEdO8KFhK98Zb1rIyKShNaN2SIiCeo6uYvxO8eZeniq/MZVKDkthbvvKxegIy9IqD5L5847w+OZZ2octYi0hJaO2SIiCcosD33LM4eTvZHeouZ7M7MNZvYR4GMJ1Wfp7A9zGrJ6dX3rISKyRJo6ZouIJCjdGzpVZ4/OJrrfqmfHNrM08DzgYuBCYC9hTF9zGRkJj8uX17ceIiI11DIxW0QkQZYON8/z2UXc6K+AihNrM9sEvA54LTAE3As8F/ge8OREa7UUZqKuf81DLSItqOVitohIkqIxG55NNrEuOxTEzF5sZt8HdgB/APw9cBLwbXe/1oMbEq3VUsgl1hpfLSItpGVjtohIgnI91iQ7EqSiHutxwi1wbwBe5+53AphZsin+Ustmw2OmNveKFxGpk9aM2SIiSapRj3XZrNLdrwGuMbMNwMVmthG4Gjg2hsLMzN2bK2irx1pEWlDLxmwRkQQdG2O91ENBctx9t7u/mzBeD+AkM3t/FLQvTrRWS2E26vtPLWpiFBGRhtRyMVtEJEGWqs1QkKqzSnefcfevufsLgM8DfwJcnmy1lpBZvWsgIlIzLRezRUSSUKP0b1EDjN39LuCtZvbLhOojIiI1opgtIlJbiYyDcPcvJLEfERGpPcVsEZHa0ABjEREREWlPCV/GrcRaRERERNpLjcZYK7EWEREREUlA+ybWmsJVRERERBLUvol1jqbbExEREZEEKLEWERERkbaU9E1olViLiIiISHvRxYsiIiIiIo1LibWIiIiISAKUWIuIiIhIe9INYkREREREFs5qNCtc+ybWmsdaRERERBLUvol1juaxFhEREZEEKLEWERERkfakMdYiIiIiIougeaxFRERERBqXEmsRERERkQQosRYRERERSUBDJdZm9gEz86jcHC3rMbNPmNkHzewfzWyg3vUUERHFbBFpAQlfvJhJdncLZ2a9wCrgwmjRvdHjJ4Fr3f0qM7sI+DTw8kUfUPNYi4gs2JLHbBGRJNXo4sWGSayBi4HbgevdfQLAzNYDrwDeFG3zQ+DrZrbR3XclclTNYy0ishD1idkiIg2sIYaCWLiv5MuAy4C9Zpbr3dgK7HP3MQB3nwF2AufVo54iIqKYLSJSTEP0WLu7A08ys0HgHcAXzewQcCKwP2/zUWB9/j6Gh4c5++yzj73etm0b27Ztq12lRUQWYfv27Wzfvj33clU961ItxWwRaRVe4dDgSmN2QyTWOe5+ELgk6g15G/ADYDxvsy5gOv+9q1ev5qabbqp9JUVEEhBPJM1sX52rsyCK2SLStKocCVxpzG6IoSAFfBzYAOwBBvPW9QMPLHmNRESkGMVsEREaN7GeBf4HuA5Yb2Y9AGbWQQjeP65j3UREZD7FbBERGiSxNrP1ZvYKM0tFpxTfAfyVuz8EXAM8Ldr0AuDr7r6nXnUVEWl3itki0jJadB7rQeDDwLuA64HLY1MzvQH4qJk9FlgH6OoWEZH6UswWkaZmNZpuuSESa3e/jXC6sNC6g4T5UkVEpAEoZouIFNYQQ0FERERERJqdEmsRERERkQQosRYRERGR9pTwxYtKrEVERESkvdTm2kUl1iIiIiIiSVBiLSIiIiKSACXWIiIiItKeNMZaRERERGQRNMZaRERERKRxKbEWEREREUmAEmsRERERkQQosRYRERERSYASaxERERGRBCixFhERERFJgBJrEREREZEEKLEWEREREUmAEmsRERERkQQosRYRERERSYASaxERERGRBCixFhERERFJgBJrEREREZEEKLEWEREREUmAEmsRERERkQS0b2Ld1wcnnAD9/fWuiYiIiIgsoXR/ms4TOkn3pRPdbybRvTWTN70pFBERERFpKxsv2cjGSzYmvt/27bEWEREREUmQEmsRERERkQQosRYRERERSYASaxERERGRBLR9Yr19+/Z6V6ElqV1rQ+1aG2rX5tHO35U+e3vSZ28uSqyb8EtrBmrX2lC71obatXm083elz96e9NmbS9sn1iIiIiIiSTB3r3cdEmFmw8C9C3jrKmBfwtURtWutqF1ro97t+gh3X13H4y85xewF0WdvT/rsjadozG6ZxFpEREREpJ40FEREREREJAFKrEVEREREEqDEWkREREQkAW2bWJtZj5l9wsw+aGb/aGYD9a5TszKzD5iZR+XmaJnadwHM7EIz+28z2xhbVrQt1c6VKdSu0fLjfrvRcrVrg2n376TYb7jVmdkLzWynme03s4+ZWabedVoqZnaemd1hZgfN7PJ616cezKzTzG42s631rkul2jaxBj4JXO/u7wKuBj5d5/o0JTPrJVy1e2FUXhStUvtWyczWAn3AOXmrSrWl2rmMYu1a4rcLatdG1LbfSYnY0NLMbAPw+4R/m38CvA54W10rtUTMbDlwLvAk4FXAm8zsgvrWqi7+FNhY70pUoy1nBTGz9cBOYNDdx6K/gA8DZ7j7rrpWrsmYWS7IfcrdJ6Jlat8FMrMUkAU2ufuuUm0JTBVbp3aeL79do2XH/Xaj5fr9Nhh9J4V/w63OzJ4C/MzdZ6LXHyV858+ub81qz8x6gAmPkjQz+yXwTne/rr41Wzpm9r+ALcBfA69x9/+sb40q06491luBfe4+BhD9o90JnFfPSjUbMzPgZcBlwF4ze3m0aitq3wVx99m8RVsp3pal1klMfruW+O2C2rURbaXNv5MCsaHluftPc0l1ZA+wu171WUruPh5LqnuB3wL/WddKLaHoM7/I3a+sd12q1a6J9YnA/rxlo8D6OtSlaXnwJMLp9E8AXzSzZ6L2TVKptlQ7L1CJ3y6oXRuRvhMBeALwqXpXYilFvfbfBTqAZXWuzlL6c+DD9a7EQrRrYu3AeN6yLmC6DnVpeu5+0N0vAT5EGP+m9k1OqbZUOy9Sgd8uqF0bkb6TNmdmpwHD7n5z2Y1by13A54CnAX9X36osDTN7BnCTuz9c77osRLsm1nuAwbxl/cADdahLK/k4sAG1b5JKtaXaOTm53y6oXRuRvpM2ZmZp4PXAX9S7LkvN3R+MhkO8k/YZ+vRO4Eoz22dm+4CTgavN7M/qXK+KtGtifR2wPro4ADPrIPyn+uO61qr5zQL/g9o3SaXaUu2cnNxvF9SujUjfSXt7J3Cpu0/WuyJ1dBPhD8x28ArgcbHyAGFGmKYYBtSWibW7PwRcQzi1AnAB8HV3b5cfbSLMbL2ZvcLMUtHFYO8A/krtu3BROwIYlP6tqp0rl9+uxX67oPjQiPSdHP8bbhdm9i7gV0CPmZ1iZn8UDQtpaWa2zMzOii16FvCxetVnKbn7sLvfnyuE2XCG3X2k3nWrRFtOtwdgZoPAR4FdwDrgL939SF0r1WTM7AzCRRUjwPXA5e5+e7RO7VslM+snzFf6SeB9wBXuPlyqLdXO5RVqV2ANRX670XvUrg2mnb+TYrGhvrWqPTO7BHh/3uLb3f30etRnKZnZOcC3gduBnwM3uvu/1rdW9WFmu2ii6fbaNrEWEREREUlSWw4FERERERFJmhJrEREREZEEKLEWEREREUmAEmsRERERkQQosRYRERERSYASaxERERGRBCixFhERERFJgBJrEWkoZrbazL5lZvea2V/Uuz4iIlKa4vYc3SBGRBqKmT0DuA4YINx17JHufqC+tRIRkWIUt+cosRaRhmVmPwHOd/eZetdFRETKa/e4raEgItKQzOxk4PvtGpxFRJqN4rYSa6kBM/szM/uHetcjJ4n6mNkFZvYzM9uaULVaipldmvD+DHgJ8OEk9ysi8yleL56ZbTSzz5jZpxPaX9rMLjazeyvY9tlm9t8L/axJxm7F7UCJdZMxs9PN7Btm9k9mdrOZuZm9qt71yvMd4Mv1rkRMEvW5D3hsAnVpSGb2LDPbZ2YPmtmTzezk6HfmZnaFma2ItjvHzG41s8vMrCdathZ4KOEqvQL4rLvPRPsXaTqK1wvSjPH64eixK6H9pYBDwIYKtr0JOL3YSjM7w8yuin57l5jZuti6pGO34jZKrJvRN4H3uvvrgLOAL+RWmNnvmNl59ahU/Njufqu7X1+PehRSTX2KtaG73wHsS7JeSX5fi92Xu38H+CQw4+4/c/f7gIuBKeC/3f1QtN2NhED+p+4+Hr39ucA1i/oAMWb2buBvgB+a2e3Ao5Lat8gSU7yuUjPGa3cfAx5I6ljuPg38ssJt9xKS8GLrbwO2Ry8vdfd4Ip1Y7FbcnqPEuomY2UrgkcBRAHefBf4acDNbTgjaVod61e3YSargc8wu4bHqsa+rgJPM7PEA7r4f+BHwnNixVgAPRL+9nC3ufvsij32Mu7/f3Te7+9nuvsXdf5LUvkWWiuJ1bTVgvE56Johq6l9u22z0mD/uObHYrbg9R4l1E4mmrrkDuNbMnhQtuxe4kTCuaROwzcxeC2Bm7zOzd0engY6NozKzs83s783sbWY2ZmZ3mNnbzez3zewXZvYCM/s3Mxsxs9fH3ldwf3nHvtjMPhcfaxadDv2EmX3IzK4zs3Oi5SWPF3u/mdnrzew2M7vQzG4xs71m9poKjrE5Xp8yxzyuDQsYjE7tjprZxXlt+jdR23wjCsSVttlxxzKzF5nZjWb2GjPbYWYfLrG/Qt99wfqUEvVs3Aa8OLZ4NfBsi4bbWxLXAAAG4ElEQVR9AM8DvhWrZw8wFnu9xcw+bGaXmtmPLQwpeamZ3WRmzzezL5vZITN7q5k9PvqM+83sWeXqJ9JMFK/bJ14XUuwzljqOmW0ysyst9P7+XSXHKXDc75jZF83s3DLbKXbXirurNFEBzgDuJfyFeiWwJrZuF7A1ev4Y4J7o+WrCX9NDhL+4dwNnR+v+Abghet4JHACuIMxF+XJgT6n95R+bMMbs08DnouXdhGRtRfT6RcAw0F/qeAU+97romO+Ktv0AMAE8oswx8utT8pjxNixQh13AZcCy6H27o+UDwFdi2/0XIShW1GZFjrWScHrvcuAc4IJS+8v77gvWp8Lf1yXA3dHzU4EPAZPAC2O/F4tt/zzgKdHzZcAvgGXR659E7TAAHAbeB/QAL41ev4Twx/07iH6DKiqtVFC8bot4Ha3/mwrbsdR3/Qvg9Nh7vMLfWe77fDLwngLrt0bHycSWKXbXqKjHusl46FU8kxBoXg3cZmaFLtLYAfyhmXUAT4mW9RH+IZ8M5MbH/pbwDwZ3nwJGga+5+wjhH9q6MvvLr98k8y+GeDYw7tEYXeAbhH+QLyxzvPz95vb5pWjbDxDG/15Y5hjz6lPNMYv4dw/j6X4BnBD7jENRj9LbCG0KFbZZkc97gBDArnb3G939h1Xsr1h9KnEVcIqF4SB/AHwG+B7wEjPrA0Y9iryRc4EboufPI/znlesFeQHw/qidDwH/4WFc9k3AgLt/1cPp8V9R3Xcg0hQUr9sjXhdQ9DOWOM4FwHJ3z9XnpiqP+QTgle7+3gq3V+yukUy9KyCVM7M0cLK77wLeYmafJVxBvR14Ynxbd580sxOBPwc+n9uFuz9sZjcBTyP8Rb2JkEwde2vs+QzRcKFi+ytS1fg+NgMdsXrNmtku4KRSxyvH3SfM7B5gBbCqimMs+Jh5Zpj797MB2Onul+dvVEmbmdkfEhLYnEd5OGXs8bpW8R0UrU857r7DzH5FGA6y3N3vMbOrCL1ILwa+Hat37reRG7/3CEIPU25f++O7jj3PMt8sGpYmLUbxet7na4d4HVe0HUt8N79D+CNioc4HNptZn7sfKbLNbPQZFLtrSA3SXDqIjX9195sJp9q25G9oZucDb3T3D3iY4SHuTcATLEz7dAvwkXIHLrO/UnYTekDj0xAZYezhYnUCd9b4GJV4kPnjkDGzJ1bRZt8CHhcrBa8ur2J/BetT+cfhK8BrgXti9UsDbyWcNs15ImG8aM4e4PfMrDd23JLj/ERamOL1fO0Ur4t+xhLHGYne08nCXEro7f9YkfWHfe6ic8XuGlJi3Xz+n5mtj70+Ebg2ej5FuFhjC2Fqp+Vm1mVmF0Tr15jZKsJYv88STuPcQHRqMZIi7690M7My+8s/tsX28U3gCPCyaF8rouPlej6LHa+YE6JtTozee20Fx4jXp9wx458jX6HeCyP0QvUC15jZMyxcJFPuOzh2LOAEd789VqZjh0jHnpfaX7zeBetjZp1mtj0a5lHKVwmnoL8G4O6jwHeB//L5s4E8nTBMJCfX3l+2MBf2m4HcXKb530H+99zUMxSIFKF4TdvE60rbseBxCLG0k/DHF4SzE5jZ+gpitxF6j/8QeLGZvSR67zvM7BJCb/2PYtsrdtdSLQZuq9SmEC6IcOAg8C/A5winkgaj9e8B7geeTzi1s4PQQ/BsQk/H5whB6luECymmov2NEf6hPROYJlwgczLwl9H6Pyq1v7xjvxX4KWHc2hnRurMJPZ2XEcYaPiZaXvR4RT6/R+9/K2FowpmxdcWOsTFen3LHjLdh3rGfT7iI74qoLXLve220/jzg1ui7uSxq50rb7PkFPuuLouN9GTgxWlbqO523ryL1OYlwI4MvVfBb+3je65cCT81bdmmB9/1v4DfAfuCSaNmzCL+1TwDriaYcA/6Y8B/vZ6LP+px6/xtTUUmqoHjdTvF6A3A9odf9cWU+Y6nv+tmEM4U/J1xIfifwxuizF4zdhLMi04Se6pXAfxCS+tcBryfM530ZYRhK7j2K3TUsFjWetAkzWwP8sbu/L3qdIvSivN7dL6lr5cowMwc2eRizKAtkZn/q7guayklElo7itcQpdjcHDQVpP28mDAtYGVt2PrGL0hpR7NSTTjstkJmlzeyVwNfrXRcRqYjitSh2Nxkl1u3nU4SxWL81s/uAHwD3u/vP6lut4qKLOd4avXyNhWnfpHp9wDfd/Z6yW4pII1C8FlDsbioaCiIiIiIikgD1WIuIiIiIJECJtYiIiIhIApRYi4iIiIgkQIm1iIiIiEgClFiLiIiIiCRAibWIiIiISAKUWIuIiIiIJECJtYiIiIhIApRYi4iIiIgk4P8DCa3RbZ0QT5wAAAAASUVORK5CYII=\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/adept-vital-venus.png',bbox_inches='tight')\n", "plt.savefig('../plots/adept-vital-venus.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/adept-vital-venus.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAm8AAAF4CAYAAAACDR42AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdd3hUdfbH8fehqoCCAioKgn3tuhHbKliwrr33HmXF3lbXtrZVsa4FjawFUSy4imL7KfYOWGDFgkpRsAQUAakh5/fHmTFjpGZucmeSz+t57nNn7kzuHCOEk285x9wdERERESkOjdIOQEREREQWn5I3ERERkSKi5E1ERESkiCh5ExERESkiSt5EREREioiSNxEREZEi0iTtAOpK27ZtvXPnzrX6GeXl5bRr165WP0OCvtd1Q9/nuqPvdd3Q97nu6Hudn+HDh09y9/l+AxtM8ta5c2eGDRtWq59RUlJS658hQd/ruqHvc93R97pu6Ptcd/S9zo+ZjVvQa5o2FRERESkiSt5EREREioiStwSVlpamHUKDoe913dD3ue7oe1039H2uO/pe1x5rKL1NS0pKXHPvIiIiUgzMbLi7l8zvNY28iYiIiBQRJW8iIiIiRUTJm4iIiEgRUfImIiIiUkSUvImIiIgUESVvIiIiIkWk4JI3M+thZu+ZWef5vLalmZ1jZvuYWdvMtaXN7FYzu9rM+pjZsnUds4iIiEhdKajepma2ItAS6Dqf144F1nD3i6q9dAfwvLs/Yma7AncBh9Z6sCIiUlzcYd68ONyhWTNoVHBjGCKLVFDJm7v/YGaDql83s27A4UCPatc7AIcBp2QuvQQ8bmad3X1sLYcrIiK1xR1+/RV+/BHKy6uOH3+En36C6dNh2rQ45x7Za3PmVCVqFRVVCVt1jRtHEtesGTRvXvW4WTNYdllo3TqO5Zarepw92raFDh3iaNMGzOr++yQNUkElbwDuXml//AtwI/A2cKuZrQlc4e5vAd2BSe4+I/O1FWY2BugGjK2zoEVEZPFVVkYSNn78H49vvoEffohEbdas5D+7USNokvmnL5vgzZwZRz6aN4eVV45ELnvu2BHWXDOONdaAZZbJP34RCjB5q87M1gE2BY519xFmdi7wfCaJWwWYXO1LpgEdqt+nvLyckpKqLhOlpaXquyYiUlsqKmDcOPjiizg+/zzOY8dGgjZnzqLvsfTS0L49tGv3+2OFFWJUrGXL3x+tWsW5RYtIpho3jkStceM4GjX6/ehYdhp19uyIJ3vMnh3H1Knwyy8wZcrvj19+gZ9/jgT0u+9g4sR479ixcSxIhw5Vydyaa8I668BGG8Hqq2v6VgAoKyujrKws+7Ttgt5XkL1NzcyBLu4+1sz2A8rc/bcNCkA5cBawLHCgu2+R87UfAA+5+/W591RvUxGRWjBrFowaBSNGwCefVCVrX30Fc+cu+OvatoVOnf54dOwIK60USVqLFnX335GvX3+tSuSyx7hx8OWXcYwZs+DvR8uWsOGGkchtvHGcN9wwElRpsBbW27TgR96IGH+L091nmtloYAViarRNtfe3AibWWXQiIg2BO3z7bSRpucfnn8fo1fysumqMLq29dtWxxhqRoNW3KcQWLapG1OZn3rwYccwmc6NHw6efwscfR6L3zjtx5Fp7bdh66zi22grWW08jdAIUx8jbusCnwEru/kPm9aHAlcB7wJdAu0xS1xSYDqzu7hNy76mRNxGRJfDdd/D++3EMHQrDhsVUYXWNGkWSsdFGsMEGkayts04kMcU0cpamSZMiEf744ziyo5jVp5aXWw623DISuWxSp+9xvVVUI29WtVvBANz9MzN7DjgAuN3MWgPNgGcyGxQGAzsCg4GdgMerJ24iIrIQU6dGcpabrH377R/ft/zyVdN62WP99WNtmtRc27awww5xZM2dG4nc22/H8c47saHjhRfigNgRu8020KMH7LwzbLqpRuYaiIIaeTOzVsARRO22y4Hb3L08U5D3FmAY0Am4291HZb6mDXAdMYW6EnCBu0+vfm+NvImIZHz3Hbz5JrzxRhwjRsQO0FytWsHmm0PXrnGUlMQ0qMphpOfbb6umV994A4YP/335kxVWgB13jGRul11ielqK1sJG3goqeatNSt5EpEFyjzVW2UTtjTdiM0Gupk1hk02qErWuXWMqVKM4hW3yZHj5ZXjxRfi//4sNErk22wz23juOjTZS4l1klLyh5E1EGpAJE2DIkDheeikWxOdq2TLWTW27bRxdu9a/DQQNTTZJzyZyL70UO2CzOneGffeFQw6JEVUlcgVPyRtK3kSkHpsyBV59tSpZ++yz37/erl1VorbttrFurUnBLXmWJM2aFX8ennwSnnoqatJldekSSdwhh0RJEiVyBUnJG0reRKQeqayM9U7PPgvPPRcbDHLXrLVsCd27x/qnnXaKTQX6B7rhmjcP3n0XHnsMHn001jxmbbABHHssHH44rLhiejHKHyh5Q8mbiBS5n3+O6bBswlZeXvVakyYxDZpN1rp2jXVsItXNmxfrHh95JJK5yZkmRU2awO67wzHHwF//qj8/BUDJG0reRKTIuMP//geDB0fC9vbbvx9d69QJ9tgDdtsNtt8+RttElsScOfDMM3DffXHOFlteeWU48cQ4Vl011RAbMiVvKHkTkSIwb14kaU8+GcfXX1e91qRJrFfbffc4/vQnTYVKcn74AR56CO6+Ozo/QPSD3WsvOOWUqEGnP291SskbSt5EpEDNmhWbDLILy3OnQ9u1iymsv/41pkPV61Jqmzu89hr06QP//S9UVMT1jTeGs86KTQ7NmqUbYwOh5A0lbyJSQGbMiGmqRx+N9Wu5JR1WXz1KOuyzT6xja9w4vTilYfv++xiJu/32GJkD6NABTj0VTj4ZWrdON756TskbSt5EJGUzZ8Lzz8dC8aefjgQua7PNIlnbZ5/Y/afpKSkks2fHlOqNN8Y6TIg+q6eeCmecEZ0dJHFK3lDyJiIpmD07+lA+8khMiU7P6dy3xRZw0EGw//6w2mrpxSiyuNxjx/O118Irr8S1Fi3gb3+Dc86B9u3Tja+eUfKGkjcRqSOVlVGK4YEHYOBA+OWXqtf+/Gc4+GA48MCoeC9SrN56C668MkaTIZK4M86Ac8+NUTnJm5I3lLyJSC0bNQr694cHH4Tx46uub7JJjLAddBCssUZ68YnUhqFD4YorYikAQJs28Pe/Q69earmWJyVvKHkTkVrw/fcwYEAkbR98UHW9Uyc44og4/vSn9OITqSvvvAMXXBA7VQFWWQX+9a/o3NCoUbqxFSklbyh5E5GEzJ0bO0X/85/YKZotbLrccjEdeuSR8Je/6B8saXiya+IuuAA+/DCubb453HwzbL11urEVoYUlb/rpIiKyOD7/HM47Dzp2jFIegwfHrtC99oqSH9myCtttp8RNGiYz2GUXGDYsujasvHJMq26zTdSHmzAh7QjrDf2EERFZkOnT4d57YyRt3XWhd++od5V9/O23MGhQjLgttVTa0YoUhkaN4Oij4Ysv4KKL4u/GI4/E35tbbqkq/Cs1pmlTEZHqRo6MCvP9+8O0aXGtRYvYKXr88VE8V7XYRBbPuHGxE/XJJ+P5ppvCnXdC167pxlXgNG0qIrIo2UKk224LG20Uydu0aZGo9e0L330X69y23lqJm8iSWG01eOKJqHW42mqxHm7LLaPdVm6xallsSt5EpGEbOzYWWHfsGDvj3nwTWraMwqMjRkSj+OOPh1at0o5UpLjtuSd88kmsHW3UCG66KUrpvPlm2pEVHSVvItLwuEcz+D33jF6i11wTDeGzI24TJ0Y/xw03TDtSkfqlRYvo0PDuu9EKbvTo2ORz5pkwa1ba0RUNJW8i0nDMnBlToBtuCD16xI7Rpk2rRtw++igabmuUTaR2lZTErtSLLopRuJtvjrIi2d6pslBK3kSk/ps4Mf6R6NgRTjwxpm5WXjkqw3/zTWxM2GYbrWUTqUvNm8ffwXffhbXXjsStpARuuy1Gx2WBlLyJSP01fHh0OVhtNbjqKpg8OfqL9u8fa90uukjNtEXSVlISHUpOOCE2Dp16Kuy9N/z8c9qRFSwlbyJSv2SrvO+4Y/yj8OCD0Sz+gANianTo0JgmbdYs7UhFJKtFiyhyPXAgtG4dvVL//OeqTg3yO0reRKR+qKiIPqObbRZV3l9+OdaunXUWfP01PPaYpkZFCt3++8co3GabwZgxUZrn3nvTjqrgKHkTkeI2Y0askVlrLTjssNh0sOKK0RR7/Hi44YaYNhWR4tClC7z1VkyjzpoFxx0Xf8flNwWXvJlZDzN7z8w6z+e1Zmb2sZl1z7m2tJndamZXm1kfM1u2DsMVkbSMHAnnnhuJ2amnxhq2tdaCsrJ4/Pe/x/SLiBSfpZaKadQ+feL5uedGf2EBoEnaAeQysxWBlsCCemacC3Sudu0O4Hl3f8TMdgXuAg6ttSBFJB1TpsCrr8YxZMjvSwpsvjmcfz7ssw80bpxWhCKStJNPhnfegX794Jhj4I03oElBpS6pKKjvgLv/YGaD5veamW0NfAf8nHOtA3AYcErm0kvA42bW2d3H1nK4IlKb5s2LzQUvvBDHe+/FxoOs1q3hkEPgqKOi1Y7WsonUT7fcEr+wvfsuXH99jKo3cAWVvAG4e6VV+yFsZi2AA9z9LDO7JOel7sAkd5+R+doKMxsDdAPG1k3EIpKYb76pStZeeilG27KaNInFyzvsANtvHwnbUkulF6uI1I3WreGee2Ij0qWXwh57NPjuJwWXvC3A+cA187m+CjC52rVpQIfqbywvL6ekpOS356WlpZSWliYZo4gsKfcomPvkk3EMH/7719dcE3beOY7tt4dltaRVpEHaeWc46SS46y44+mh4//16OX1aVlZGWVlZ9mnbBb2v4P/LzWwXYJi7/ziflx2YWe1ac2Bu9Te2a9eOYcOG1UKEIrJE5s2LNSzZhO2rr6pea9ECdtopfsPeZZfoOyoiAjFl+vzzUfvtllvg7LPTjihxuQNLZjZpQe8r+OQNOAfYNGcqtQ0wyMyuAr7JPM/VCphYd+GJyCJVVkbC9tBDUYTzx5zfxdq1g732is0GO+4ISy+dXpwiUrhatozdp7vvDpdcEjXhOndOO6pUFEPydhgxmpb1DnAW8AKwDHC3mS3t7jPNrCnQCXit7sMUkd9xhxEjonDugAFRcy2rSxfYd99I2LbeWjtERWTx7LYbHHQQPPoonHIKDB7cIDcrFVzyZlVDbAbg7uXVXp8HlLv7VGCqmQ0GdgQGAzsBj7v7hDoMWURyjR0bvUMHDIBRo6qur7oqHHpoHJts0iB/4IpIAm65JTY1PftsjOQfeGDaEdW5gkrezKwVcETm6TFmdlv15G0+egLXmdnGwEqAdiGI1LVZs+CJJ2JH2JAhMeoGsMIK8YP1sMOiNVWjgqsLLiLFZqWV4NprowbcaadBjx4NriC3efaHbD1XUlLi2rAgkrAPP4yE7cEH4edMCcbmzWG//eCII+KHatOm6cYoIvVPZSVst1200Tr55KpODPWImQ1395L5vVZQI28iUgSmT49p0bKySN6yNtsMjj8+pkXbVN9HJCKSoEaNomzIJpvAnXfCkUfG+tkGQnMYIrJ4Pv8cTj8dVlkFevaMxK1Nm+gr+uGHUaPtb39T4iYidWP99eG88+LxySdDRUW68dQhJW8ismDz5sFTT0WBzHXXhX//G6ZOjfVrDz4IEyfGtU02STtSEWmILrooyoWMHBmzAQ2EkjcR+aOpU+HGG2GNNWDvveHFF6P+2gknxCjbm2/GJgS1pxKRNC29NNxwQzy++GKYXL3pUv2k5E1EqkyYENMQHTtG9fJx46LLwQ03xGt3361RNhEpLPvuGz2Pf/opivc2AEreRCSK6R59dEw/9O4dI2/dusWU6ejRcNZZWssmIoXJLGq/NW4cmxdGjEg7olqn5E2kIXvjDdh1V9h4Y+jXL7bfH3xwNH1+9VXYc0/VZhORwrfBBrFhqrIyNlbV8zJo+qks0tC4RyHd7t2jTtILL0RD+NNOgy+/hIcfhs03TztKEZEl889/RmHwV1+Fxx9PO5papeRNpKFwh+eei52iO+0Er70WVckvvTT6jt5yS/QcFREpRm3awJVXxuOzz4YZM9KNpxYtVvJmZq3NrFHO8/Zm9pfaC0tEEuMezZu7doXdd4d33onfTq+6KvqQXnYZLL982lGKiOTvxBNjGcj48VHGqJ5aZPJmZucDo4GXzexCM2vk7j8CDaegikixevXVGGnbc08YNgzat4frrouk7cILYbnl0o5QRCQ5jRvHpiuAa66JHaj10OKMvJ0MrOvu3YF+wOlmthzwa20GJiJ5GDYMdtkFtt8+Rtrato1yH2PGwLnnQsuWaUcoIlI7evSAHXeEX36JBK4eWpzkbQTwC4C7f+vuNwF7AMvWZmAiUgOjRsH++8eGg//7P1h2Wbj8cvj66yj3scwyaUcoIlL7sknbrbfCt9+mG0stWJzk7Xhg+9wL7v4Q8ECtRCQiS+6HH6K334Ybwn//G1XHzzsvkraLL4ZWrdKOUESk7pSUwEEHwaxZsa63njGv57VQskpKSnzYsGFphyGSrJkz4aab4rfMadNivceJJ0bC1qFD2tGJiKRn9GhYb72o/TZyZDwuImY23N1L5vdarZUKMbPda+veIg1eZWU0hl9nHfjHPyJx22OP+AHVp48SNxGRtdaKfsyVlfFzsh7JO3kzs13N7E0z+9LMvs4cY4H6XSFPJC3vvw9bbQVHHAHffBPb4l96KcqB/OlPaUcnIlI4Lrkk1vo++SS8+27a0SQmiZG3B4D/AicCx2aOY9CaOJFklZfHlOiWW0YCt/LKcM89MHx47KwSEZHfW3nlaJcFsXmrnsh7zZuZvefuW8zn+kru/n1eN0+Q1rxJ0Zo3D+66Cy66CH7+GZo2hTPPjHVtKvkhIrJwkyZB587w66/xi2+RtP+r7TVvPc3sgPlc3yuBe4s0bO+8Ez9oTjklErcePWDECLj2WiVuIiKLo23baFoPVe2zilwSyds2QF8zm5dzVAJ9Eri3SMP0yy/xw2brreHDD6FjRxg4MJrIr7tu2tGJiBSXs8+OEkpPPQUff5x2NHlLInm7HCgF1gJWzzkuS+DeIg3PE0/ElvY+faBJE7jgAvj00yi+a5Z2dCIixWfFFaG0NB7Xg9G3JJK354BB7v61u4/LHGOBuxO4t0jDMWEC7LdfHBMnxsaEDz+Eq6+GFi3Sjk5EpLidey40axazGJ98knY0eUkieXsYuMDMtss5ugEXJHBvkfqvshLuvDNG2554Iroh3HYbvPkmbLBB2tGJiNQPq6wSdd8Arroq3VjylMRu06HElOkUIHszA1Zy96XyCy852m0qBWncODj+eBgyJJ7vtRfcfjusumq6cYmI1Efjx8Oaa8Yu/i++gDXWSDuiBart3aa9gY7u3tndu2SOzsC+CdxbpH5yh759oxfpkCGxG+rRR6OQpBI3EZHa0akTHHZYzHjcemva0dRYEsnb28BFZtYSwMy2NrN93P25mtzMzHqY2Xtm1jnn2n5mNsbMJpvZLWbWJOe1tmZ2l5lda2Y3mlnTPP97RGrXt9/CbrtFwd1p02KN2yefwIEHakOCiEhtO+OMOP/nP7Gzvwglkbz1B7YAmgK4+9vAKmZ24ZLeyMxWBFoCXXOudQL2AQ4AzgROAM7I+bLHgHvc/XzgC2L3q0jhcYf+/WMd2wsvQJs28NBDsXi2ffu0oxMRaRg22QS6d4fp0yOBK0JJJG9D3b27u/+cc204kWgtEXf/ARhU7fJqwHHuPtzd+wG3A9sDmNmWwDru/l7mvU8Dp5qZtuZJYZkyBQ4/HI48Mn7T23PPGG079FCNtomI1LUzMynKv/8NFRXpxlIDSSRvFWbWPPvEzJYmdpr+UJObuXtltedvuHvud3YCMD7zeAdgTM57J2QezneBn0gq3ngjmscPGBANkvv2hUGDoueeiIjUvb/+NTYujBsXP4+LTBLJ233AG2b2mJkNAr4FtgVOTuDe87M5cGfm8SrA5GqvTwM6VP+i8vJySkpKfjvKyspqKTyRjLlzox9p9+6xw2nzzeGjj2J3qUbbRETS06hRVcP6m25KN5YcZWVlv+UpQNsFvS/vUiHw22jbXsQU5wTgGXefksf9HOiSKfabe31N4BR3PzPz/DagnbsfnPOen4BSdx+Y+7UqFSJ1auxYOPjgaIJsFl0SLrssmsqLiEj6pk+P1oNTphRkw/pES4WYWW8zO8vM1s5ec/eZ7v6Iu1/n7g/mk7gt5HMbAycBf8+5PAFok/MeA1oAE5P+fJHFNmgQbLpp/DDo1AlefTUKQipxExEpHC1bVhXtvfPOhb+3wCxx8ubu5xI7TLfNJHKXmdk2mcSpNp0DXO/us3OuvQisk/O8MzAH0BCb1L25c6P58T77xG9ye+0V7a222y7tyEREZH6y/U4ffrioyobUaM2bu//o7v/JJHLXAMsD/zKzf5nZvma2TE0DykkCLefahcCHwNJmtrqZHWdma7r7MOA7M8v2ENoVuNXd59T080VqZPz4SNJuvDGayd9wQxTcXX75tCMTEZEFWWst2GEHmDEDHnww7WgWW5NFv+WPzKx5dgTM3WcRJTqeziReWwDnmFlLdz9vCe/bCjgi8/SYzJq2k4Arqr31M3e/J/P4YOASMxtL1Ij7R03+m0Rq7Jln4Kij4KefYv3EI4/AVlulHZWIiCyO0lJ4+WW46y7o2bMoNpTVaMOCmfXOjLoVDW1YkMRVVMRu0muvjee77w79+sEKK6Qbl4iILL7Zs6Mt4aRJ8N570LXror+mDtRGb9NdzeweM7vSzArjv1KkLk2aBLvuGolb48ZxfvppJW4iIsWmefOYPQG4//50Y1lMNR1529zdh5rZckSJkC2AqUSJkLcSjjERGnmTxHz8cWxKGDs22lo99pg2JYiIFLMRI6KY+vLLw8SJkdClLPGRN3cfmjn/4u4PuHsv4CpgLTP7wMxuNbNuNQ9ZpEA9/HCsZxs7NmoCDR+uxE1EpNhttFEkbz/9BM8+m3Y0i1Sj5M3Mls153NzM9gHuBm4F1iU6HKjTttQf8+bBeedFL9KZM+GYY+D112OdhIiIFL8jj4xzv37pxrEYarrm7UozO9zMBgDlwACiOO5JQHt339/dH0sqSJFU/fQT7LYb9O4d69tuvRXuuQeWWirtyEREJCmHHRZts555BiZX77xZWGqavPUC+gLLAp8BNwGnu/tD7j49qeBEUvfZZzE9+uKL0K4dDBkCvXoVxVZyERFZAiuvDDvuGAXXn3oq7WgWqqbJ2xPAKu6+h7t3JRK5A83sdjM728xWSy5EkZQMGQJbbglffw2bbQbDhkE3LeUUEam39t8/zo8/nm4ci1DT5O10d/8p+8TdvwZeInacngt8ZWY3JxCfSDr69o1SIL/8AvvuG+vbOnVKOyoREalN++wTMysvvghTp6YdzQLVNHlbBcDMNjWzq81sNNFPtAdwA7CGu5+RUIwidaeyMjYmnHhiFOE97zwYOBBatEg7MhERqW0rrgjbbgtz5sDgwWlHs0A1njY1s3HAcGBP4D5gHXcvcffe7j4uqQBF6syvv8IBB8TGhCZN4O67o/huo5r+NRERkaJTBFOnNf1XqSmRsG3g7hu6+1Xu/mVyYYnUse++i/VsTzwByy0Hzz8PJ5yQdlQiIlLX9tsvzs89F7/UF6CaJm9Huful7j4q0WhE0vD551F4d/hwWH11ePfd2HEkIiINz6qrwhZbRE3P555LO5r5qmmHhcL8rxFZUu++C9tsA+PGxV/Wd9+FdddNOyoREUnTvvvGuUC7LWgxjzRcgwfDDjtEMcY99ojSIO3apR2ViIikbZdd4vzCC1CDHvC1TcmbNEz33BNbwmfOhOOOgyef1I5SEREJG20UO08nToRRhbdCTMmbNCzucOWVcPzx0a/0oouipluTJmlHJiIihaJRI9h553j8wgvpxjIftZK8Zeq/9c90W+hcG58hssQqK+HUU+Hii6MI4+23wxVXqNWViIj8Ue7UaYGpreGGk4Gd3P0IMzvMzH5w95m19Fkii1ZRAcceC/37Q/Pm8OCDVbV8REREquvRI86vvx5LbJZeOt14ctTWtGkv4M+Zxw8TdeFE0jF7Nhx0UCRuLVrE7iElbiIisjDt28Omm8KsWfDGG2lH8zt5J29m1srMrjWzazPPuwKHuPsEAHevdPfCbRAm9duvv8Kee0bx3dat4aWXYoepiIjIomRH3157Ld04qkli5O0+YBugDYC7vw80N7PeCdxbpOamTIk1Cy++GL9BvfoqbLll2lGJiEix2GabOL/9drpxVJNE8jbX3f8CfJJzbTxwXAL3FqmZ8vIYYXvrraiW/frrsPHGaUclIiLFZKut4vzeezB3brqx5Egiefs0c3YAM2sOnAdMSODeIkvuhx+ge3f48ENYYw14801YZ520oxIRkWLTrh2svXZsWPjoo7Sj+U0SydtgMxsA7GFmdwJfEZsVTkng3iJL5vvvI3EbNQrWWy8Wma62WtpRiYhIsdp66zgX0NRp3smbuw8HSoF7ga+BC4A13b2wtmZI/ffdd7D99vDZZ7DBBvDKK7DyymlHJSIixSy77u2tt9KNI0fedd7M7Gh3v58oCZK9to6Zbe3uT+d7f5HFkk3cPv8cNtxQfUpFRCQZBTjyVuPkzcy2BQzoYWZjqr28PPBvYImTNzPrAVwJHOzuYzPXlgauA6YRu1rPz5YfWdhr0kBMnBiJ2xdfRD+6IUOgbdu0oxIRkfpg3XWhVSuYMAF+/DGqF6Qsn5G3KUA/oAtRKiTXHKDvkt7QzFYEWgJdq710B/C8uz9iZrsCdwGHLsZrUt9NmBCJ2+jRsZv0pZeUuImISHIaNYoZnbffhhEjYKed0o6o5mve3H0k0A04x9275B7A1u5+eQ3u+QMwKPeamXUADqNqFO8lYC8z67yw12ry3yRF5ocfohzI6NGwySYacRMRkdqRLTU1YkS6cWTktWHB3ae6e18za2NmHc2sk5l1ArqY2T01vGdltUvdgUnuPiPzegUwhkgcF/aa1Gc//RSVr9IxxMcAACAASURBVHOnSldYIe2oRESkPtpoozh//HG6cWQksWGhN3A2UefNMpcdGJnvvTNWASZXuzYN6ABULOS13ykvL6ekpOS356WlpZSWliYUotSpqVNh111h5Mio3/bii7D88mlHJSIi9VU2eavlkbeysjLKysqyTxc4lZR38gZ0BNoD2wJD3f1bM9uJ5JrROzCz2rXmQLbU8cJe+027du0YNmxYQiFJan79FfbYA4YOhS5dYsStABaPiohIPbbhhnEeNSo6LTRNKsX5vdyBJTObtKD3JVGk93/uPolYd3ZU5trrwE0J3BuiU0ObatdaARMX8ZrUN7Nmwb77RseEVVaJxG2VVdKOSkRE6rtWrWD11WHOnFhnnbIkkremZvYlsDkw2sxeAd5J6N4ArwAdMiVBMLOmQCfgtUW8JvXJ3Llw8MFVTeaHDImRNxERkbqw1lpx/uqrdOMggWlTd7/UzB4FRrv7O5lhvk2BgTW5n5ll181Z5v7fm9lgYEdgMLAT8Li7T8i8f4GvST0xbx4ceSQ89RS0aRMJnHqViohIXVp99Th//XW6cZDMhoWrgZfc/RMAd3+FGBGryb1aAUdknh5jZre5eznQE7jOzDYGViLacWUt7DUpdu7wt7/BI4/EsPULL1QtHBUREakr9Sl5Aw4A/lv9opmt7O7fLcmN3H0a0Cdz5F7/GThxAV+zwNekHrjsMigrg6WWgsGDYfPN045IREQaomzyNqZ6U6m6l0TydhVwmJm1IHaGAjQGTgAOT+D+0lDdcQdcfnlUt37kEdhuu7QjEhGRhqqejbz9DVgb2I+q5A1iClOkZgYOhF694vFdd8Fee6Ubj4iINGzZTXJffx1Len5bol/3kkjergNezkxf/sbM9k/g3tIQvfoqHH54/OW48ko44YS0IxIRkYZuueXi+OUX+PnnVIvD513Ow90fzyZumRZZjbPX8723NEAffQR77x21dHr1ggsvTDsiERGR0K5dnMvLUw0jqVpsWZ8DayR8T2kovv4adtst2l8ddBDcfHOqw9IiIiK/U0+TN/1LKzUzeXL0K/3+e9hhB+jXDxo3TjsqERGRKvU0eRNZcrNnwz77RMuRjTeGJ56A5s3TjkpEROT36kvyZmZH5zw9EfjOzNYxsz3zvbc0AO5w3HFV/UqfeQaWXTbtqERERP4om7z9+GOqYdR4t6mZbUtMk/Yws2zFuvFEa6zlgX8TzepFFuzSS+Ghh6Bly0jc1GheREQKVdu2cZ48OdUw8ikVMgXoB3QBtqn22hygbx73lobgvvvgiiuqivBuvHHaEYmIiCzYMsvEedasVMOocfLm7iPNrBtwkLsrUZMl8/LLcGKmq9ltt8Huu6cbj4iIyKIstVScU07e8lrz5u5T55e4mVl6leuk8H36Key/P1RUwFlnQc+eaUckIiKyaNnkbebMVMNIZLepmbXJFOjtZGadgNXN7J4k7i31zKRJsMceMGUK7LsvXHdd2hGJiIgsngIZecu7PZaZ9QbOJvqaZuu8OTAy33tLPTN3Lhx4IIwZA3/+M/Tvr1puIiJSPOpL8gZ0BNoD2wJD3f1bM9sJaJrAvaU+OeOM6Fu60krw5JNVCz9FRESKQbNmcZ49O9Uwkpg2/Z+7TyLKghyVufY6cFMC95b64s474Y47ovjuE0/AqqumHZGIiMiSmTcvzinPGiUx8tbUzL4EjgRGm9krwLKoe4NkvfYanHpqPC4rgy23TDceERGRmpg7N85N051czDt5c/dLzexR4At3f8fMJhGFegfmHZ0UvzFjqnaWnnMOHHXUor9GRESkENWX5M3MWlE1XXo+8CtQ7u7j8723FLnp02Hvvauazl9zTdoRiYiI1Fw2ecuufUtJElOb9xEdFtoAuPv7QPPMLlRpqNzh2GNh5EhYZx0YMCD1NQIiIiJ5KZCRtySSt7nu/hfgk5xr44HjEri3FKsbboCBA6PJ/FNPQevWaUckIiKSn2yJkHow8vZp5uwAZtYcOA+YkMC9pRi98gqcf3487tcP1l473XhERESS8NNPcV4+3UZSSew2HWxmA4DlzWw94K9AC2CvBO4txeabb+Dgg6GyEi68MNa8iYiI1AeTJ8d5hRVSDSOJ3abDzawU2APoBFwAPOvuk/O9txSZ2bPhgAOgvBx23hkuvzztiERERJKTTd7atk01jCRG3nD3acDDudfM7Gh3vz+J+0uROP10eP99WG01eOghbVAQEZH6pUBG3vJe82Zmx5vZODObbWbzMkcloMb0Dcl998Fdd0UHhf/+N/U/2CIiIokrkOQtiZG3m4BSYChQkblmwCEJ3Ps3ZvYnoBfwJbAWUObuH5nZ0sB1wDSiXMn57j41yc+WRRg5Enr2jMd33AGbbZZuPCIiIrXhu+/i3L59qmEkkby9CDzp7rNyL5rZgwncO9cDwN7uPsHMOgEvAH8C7gCed/dHzGxX4C7g0IQ/WxZk+nQ46KDYPn3ccXGIiIjUN5WVMHZsPO7cOc1Iljx5M7Ptql3qB1xhZk9Xu34AcFpNA5uPPwGtMo9nAsuZWQfgMOCUzPWXgMfNrLO7j03ws2V+3GPE7bPPYP314dZb045IRESkdkycCHPmxKhby5aphlKTkbdrgI2BSUBlzvUDch4bsCLJJm8Dgb5mthtwBHAq0B2Y5O4zANy9wszGAN2AsQl+tszPvfdC//6wzDLw2GNxFhERqY/GjIlzly7pxkHNkreLgW/c/YuFvcnMkq7z9jdgELG27gZ3f9zMzgWqlySZBnSo/sXl5eWUlJT89ry0tJTS0tKEQ2xARo6EUzIDnn36wJ/+lG48IiIitenrr+O8+uq19hFlZWWUlZVlny6wHskSJ2/uPmQx3/fUkt57EZYGBgCrADdnRticmELN1RyYW/2L27Vrx7BhwxIOqYHKXed27LFw1FFpRyQiIlK7sslbLY685Q4smdmkBb0vkTpvdaQ/cIi7TzEzJxK5M4gdprlaARPrOrgGpVevqnVut92WdjQiIiK17+OP47z++unGQTK9TWudmbUFNnb3KZlLVwJTgXFAh0y5EMysKdHl4bVUAm0IHn4Y7r8fll4aHn1U69xERKRh+PDDOG+6abpxkEyR3qPnc20dM9sz33vn+AmYZWar5FybDIwABgM7Zq7tBDzu7hMS/GzJGjcOTj45Ht90E6y3XrrxiIiI1IXJk2H8+BiwWHvttKOp+bSpmW1L7CrtkVl/lmt54N9A9fIhNeLulWa2D3CJmQ0ndrKe6+5TzawncJ2ZbQysRBQMlqTNmwdHHAG//BLN5rXZQ0REGoqPPorzRhsVROvHfNa8TSFqvHUBtqn22hygbx73/gN3/xg4aT7XfwZOTPKzZD7+9S94801YeWXo2xfM0o5IRESkbhTQlCnkkby5+0gz6wYc5O6JJmpSYN59Fy67LB7ffz+0XeDuZRERkfrn7bfjnFNyLE15rXlz96kLStzM7Jh87i0FYupUOPzwmDY96yzo0SPtiEREROqOO7z+ejzu1i3dWDKS2LBwopl9a2ZzzWxe5qgE/pNAfJK2M8+M2jabbAJXX512NCIiInVr1KjYsLDKKrVaoHdJJFHn7XpiLdpQoCJzzYBDEri3pGnwYLjnHmjeHB58MM4iIiINycsvx7lbt4JZ751E8jYEGOzu03MvmtkDCdxb0jJ5MpyY2Qdy1VUqCyIiIg3Ts8/Geeed040jRxLJWz/gajMbmHPNiEb1pyZwf0nDqafC99/DX/4CZ5yRdjQiIiJ179df4ZVXYsRtt93SjuY3SSRvFwJrAHsRvUYh1tKthJK34jRwIAwYEMUI77uvIGraiIiI1LmXX4bZs6FrV2jfPu1ofpNE8nYV8Ky7/64ZvJntmsC9pa798AP07BmPe/eGNdZINx4REZG0DMxMKu6ZZNOo/OWdvLn7oAW8NC7fe0sdc4/2V5MmwU47VbXCEhERaWhmzYInn4zHBx+cbizVLHHyZma3AgPc/e3M8z5A9W2IBnQF1s87Qqk7/fvHH9Rll4X//Aca5V1JRkREpDg9/3zUOt10U1hrrbSj+Z2ajLyNJ1pjZS0FzAa+p2rNmxHr4KRYTJwIp50Wj2+5BTp1SjceERGRND2QKZpRYKNuUIPkzd17V7t0FTDe3efkXjSze/MJTOpYr14wZQrssQccfXTa0YiIiKTnhx/gqadiw96RR6YdzR8ksebtSzNrYWYHA6sCnwPPuLvWvBWLxx+HJ56AVq2gT5+CKUIoIiKSin79oKIC9toLOnRIO5o/yDt5M7PNgWeI8iBjiSnTK8zsQHcfle/9pZb9/HOMugFccw107JhuPCIiImmaNw/uuiseH398urEsQBIr0q8FegLt3L3E3f8MbANoq2IxOOecqmK82l0qIiIN3eDB8NVX0LlzLCUqQEnUeXvX3R/PveDuU8xsdgL3lto0ZEj0Lm3WDO6+W7tLRUREbropzqedVrBF6pP417qlmbXKPjGzpmZ2JLBlAveW2jJjBpSWxuNLLoF11003HhERkbQNGwavvRZrwAt0yhSSGXm7A3jHzGYATYHVgelAYZUjlt+75BL4+mvYaCM477y0oxEREUnf5ZfH+eSTo+ZpgUpit+lnZvZnYA9gLaIO3JPuPjPfe0st+eCDGBZu1Aj69oWmTdOOSEREJF0ffABPPx19vc85J+1oFiqpRU7tgdbAXGBa5iyFaN68+I2isjLm8zffPO2IRERE0nfppXHu2bOgmtDPTxKlQo4G7gTmAF8DzYDrMqVCPsn3/pKwu++GoUOjbk12eFhERKQhe+212GXasiWce27a0SxSUqVCrgfau/um7r4+0APolcC9JUk//ggXXBCPb745FmSKiIg0ZJWVVdOk550HK66YbjyLIYnk7QfgWnf/rTSIu08APs0+N7NlEvgcyde550YLrF12gQMOSDsaERGR9D34YOwyXXllOOustKNZLEkkb2cDvcysU86xBbCmmXU0s9WAwh+DrO9eey3afTRvDrfdphZYIiIiP/0EZ58dj6++Glq0SDeexZREqZCbgfWAq+fzWnbq1IF/JvBZUhNz5sDf/haPL7gA1lwz3XhEREQKwd//DuXlsN12cPTRaUez2JJI3v4NDHL3H+b3opkZcFoCn5O935bAX4AvgTeBX4HriF2ubYDz3X1qUp9XL9x0E4waFUnb+eenHY2IiEj63nwzNvE1bRq9TItoRiqJOm9li3jdzeyOfD8HwMyOBdZw94tyrt0LPO/uj5jZrsBdwKFJfF69MG5c1a7S22+HpZZKNx4REZG0zZlT1c/7738vui5DddLM0t3zrvtmZt2Aw4GLc651AA4Dns5cegnYy8w65/t59cY550QrrIMOgp13TjsaERGR9F18MXzyScxIXXhh2tEssWLqRH4jsYP1VjN73sy2AboDk9x9BoC7VwBjgG6pRVlIXn0VBg6MatE33JB2NCIiIul7+WXo3TuazvfrV5QzUkWRvJnZOsCmwN3u3gsYAjwPdAQmV3v7NKBD3UZYgCoq4PTT4/EFF8Cqq6Ybj4iISNp++gmOOgrcY/Rtq63SjqhGai15M7NjErzd+sBP7j4i8/w2wIhdrNV7qDZnPu25ysvLKSkp+e0oK1voUr3i17cvjBgBq61WtQ1aRESkoXKH0lKYMAG23hr+8Y+0I/qDsrKy3/IUoO2C3pdEe6wTgUuBFalKBrOJ1X353j+jCTmxuvtMMxsNNCZ2mOZqBUysfoN27doxbNiwhMIpcD//DBdl9nRcfz0svXS68YiIiKTt1lvh8ceju1D//tAkiYIbySotLaW0tBQAM5u0oPclEfn1wEnAUKAic82AQxK4d9YIYDkzWzGnJEkFkaR1MLOlMwldU6AT8FqCn118/vlPmDwZunWD/fdPOxoREZF0vfJKVfeEvn2hS5d048lTEtOmQ4DB7v6Vu4/LHGOBBxK4NwDu/hnwHHAAgJm1BpoBDwKDgR0zb90JeDzTnqthGjUqOig0ahT9S4uobo2IiEjixo6FAw+EefOiLMhBB6UdUd6SGHnrB1xtZgNzrhmRaJ2awP2zjgJuMbNmxOjaoe5eYWY9gevMbGNgJaA0wc8sLu5w5pnxB/Skk2CTTdKOSEREJD0zZsC++8Zs1K67wpVXph1RIszd87uB2fvAGsQuz+zNGgEruXvz/MJLTklJidf7NW/PPQe77w7LLQejR0O7dmlHJCIiko5582LE7Yknop7b++9Dm+rL5AuXmQ1395L5vZbEyNtVwLPVC/Fmuh1IXZk3D847Lx5ffLESNxERabjcoVevSNxat4ZBg4oqcVuUJNpjDTKzFmZ2CLAq8BmRzD2fd3Sy+Pr1g//9L0qDnHJK2tGIiIik58or4c47oXlzeOopWG+9tCNKVBKlQjYHniGmSscS692uNLMD3X1UvveXxTBjRlVpkKuuKspq0SIiIom4+2645JLYuDdgAGy7bdoRJS6J3abXAj2Bdu5e4u5/BrYBTk7g3rI4br4ZJk6EzTaDQw9NOxoREZF0DBhQ1XD+9ttjs0I9lMSat3fd/fHcC+4+xcxmJ3BvWZTycrjmmnjcu3f8piEiItLQDBgARxwBlZVR7/Tk+juGlMS/9C3NrFX2iZk1NbMjgS0TuLcsyhVXwLRpsNtusMMOaUcjIiJS9x5+uCpxu/TSmDatx5IYebsDeMfMZgBNgdWB6cCeCdxbFubLL6FPnxhtu/batKMRERGpew8/DIcfHonbJZfAZZelHVGtS2K36Wdm9mdgD2AtYDzwpLtXbxgvSbvsMqiogGOOgQ03TDsaERGRunXXXdCzZ5QGaSCJGyQz8oa7zwb+m3vNzP7k7p8mcX+Zj08+gYcegqZNY4hYRESkoXCPZUPZf/+uuAL+8Y8G0xJyiZM3M7sVGODub2ee9wGqd1IwoCuwft4Ryvxdemn84T3xROjcOe1oRERE6sa8eXD66bGbtFEjuOOOaAnZgNRk5G08MCXn+VLAbOB7qtpjGdEyS2rDBx/A449HPbd//CPtaEREROrGzJmxVOjRR6FZs9hhut9+aUdV55Y4eXP33tUuvQXcn9sey8zWASbnGZssSHYXzSmnQIcO6cYiIiJSF77/HvbZB957D1q1ipZX22+fdlSpqPGaNzPblhhh6w58bva7ieY2wLnArXlFJ3/0zjvwzDPQogWcf37a0YiIiNS+jz6CvfaCb76BTp3g6adho43Sjio1+WxYmAL0A7oQHRVyzQH65nFvWZCLL47zGWeo+byIiNR/gwZFKZBff4Wttopm8yuumHZUqapx8ubuI82sG3CQuytRqwuvvAJDhsByy8HZZ6cdjYiISO2prIxdpP/8Z2zQO+KI6Fuq/t35lQpx96nMZ4TNzJZ395/yubfMx+WXx/nss6FNm3RjERERqS2TJ0ey9vzzUf7j6qvh739vMKVAFiWROm9m1gZoSayBA2hvZte7+3FJ3F+At96CV1+NUbfTTks7GhERkdoxdCgccACMHw8rrBA7Snv0SDuqgpJ38mZmvYGziTIh2eTNgZH53ltyXHVVnE87LRI4ERGR+sQd7rwz1nTPmQNdu8Jjj8UGBfmdJEbeOgLtgW2Boe7+rZntRPQ5lSQMHw7PPRc7TE8/Pe1oREREkjVpEpxwQmxOgCiFdcMN0Lx6DwABaJTAPf7n7pOAp4GjMtdeB25K4N4CVaNuf/tbDCGLiIjUF0OGwMYbR+K23HLRaP6225S4LUQSyVtTM/sS2BwYbWavAO8kdG/53/9iW/RSS8FZZ6UdjYiISDLmzIlNCD16wMSJsM028PHHcPDBaUdW8PKeNnX3S83sUWC0u79jZpOBTYCBeUcnscMGoofpSiulG4uIiEgSPvoIjj4aRoyI/qSXXQYXXghNEtlHWe8lMjrm7p+4+5zM45fd/UZghyTu3aB9+SU88gg0bQrnnpt2NCIiIvmZOzdqt22+eSRuq68Or78ebR+VuC22vJM3MzvRzL41s7lmNi97AP9JIL6G7aabokjhEUdAx45pRyMiIlJz//sfbLllJGoVFdCrVyRw21Rv0iSLkkSaez1wEjAUqMhcM+CQBO7dcE2eDPfeG4/VTUFERIrVrFmxBOiaa2LkbbXV4J57YAdN0NVUEsnbEGCwu0/PvWhmDyRw74arTx+YORN23x3WXz/taERERJbckCHQsyeMHh3PTzoJeveGVq3SjavIJZG89QOuNrPcDQoGHACcmsD9f8fMmhGjfKe7+6tmtjRwHTANaAOcn2nbVbxmzYJbb43H55yTbiwiIiJLqrw8/v3q1y+er7ce3HUX/OUv6cZVTySRvF0IrAHsRXRWgFhLtxK1kLwB5wKdc57fATzv7o+Y2a7AXcChtfC5dad/f/jxR9hsM+jePe1oREREFo873HdfJG4//RS12i6+ODbdNWuWdnT1RhLJ21XAs+4+N/diJpFKlJltDXwH/Jx53gE4DDgl85aXgMfNrLO7j0368+tEZWVUlYb4w68mvCIiUgw++CC6AL35Zjzv0QPuuAPWXDPduOqhJEqFdKmeuGVUmNnlZrZFAp+BmbUADnD3e3IudwcmufsMAHevAMYA3ZL4zFQ8+yx89ln0cjvggLSjERERWbgff4xapCUlkbi1bw8PPggvvKDErZYkkbydbGZvmtlwMzsPILMO7UliCnMdM9s+gc85H7im2rVVgMnVrk0DOiTweem4KdNV7PTTo76biIhIIZozB268EdZaC/r2hcaNozrCF1/AYYdp5qgWJdWY/iXgU6CLmR1JtMdaxt0nZHad/h/wSk0/wMx2AYa5+4/VXnJgZrVrzYE/jASWl5dTUlLy2/PS0lJKS0trGlLtGDUKXn45GtAfd1za0YiIiMzfs8/CmWdGogZRGeHGG2GdddKNq8iVlZVRVlaWfdp2Qe9LZMOCu9+SfWJmhwHNgGzHBTeztfP8jHOATa0qi28DDAJuzDzO1QqYWP0G7dq1Y9iwYXmGUcvuuCPORx4JrVunG4uIiEh1H3wQ/UhffDGer712zBjtvnu6cdUTuQNLZjZpQe9LInnrZGabArOAHYkG9d8Cv2Q+vA0xGpaPw6rd4x3gLOA94FwzW9rdZ5pZU6AT8Fqen1f3pk6F+++Px6ecsvD3ioiI1KWvvoKLLoKHH47nyy0XnRJ69dIu0hQkkbw9RKxv60jUX+sF3AO8a2YXAesCL+bzAe5envs8036r3N3Hm9lgImkcDOwEPO7uE/L5vFT06wfTp0O3brDBBmlHIyIiEpsRrrgC7rwzWlo1bw6nngoXXADLL592dA1W3smbuw8HVjOzdtkky8y6u/tkM2tHJHN98v2chegJXGdmGxO15QpsIdticIfbbovHvXqlG4uIiMjUqTEdev31MbBgBsccA//8Z1RDkFQlMfKWnRpdysyy/0fbm1lvdz8OuDSJz8jl7p1zHv8MnJj0Z9SpIUPg889hlVVg773TjkZERBqqadOiw88NN0SRXYA994zepJoVKhh5J29m1hs4m9j5md1R4MDIfO/dYPTJDEyedJLKg4iISN2bPj1mgK6/HiZnKnBtuy1cdVWcpaAkVSqkPbAtMNTdvzWznQBlIYvjxx/hqaegUSM4/vi0oxERkYbk11/h9tujWfykzObGbbaJ6dEddlCttgKVRJHe/7n7JOBp4KjMtdeBmxK4d/3Xr18sAt1jD+hQvLWFRUSkiPzyC/zrX9ClC5x/fiRuW20F//d/8MYbsOOOStwKWBIjb03N7EvgSGC0mb0CLEsyiWH95g7/+U881qibiIjUth9+iI0IffrEpgSALbeMkbYePZSwFYkkdpteamaPAqPd/Z1MUblNgYF5R1ffvf129DFdaSUVOBQRkdozZkxMjd5zD8yeHdd22CEK7u60k5K2IpPEhoWj3f3+7HN3f8XMJgIbA+PzvX+9lh11O/pobVQQEZHkffxxJG0PPwzz5sW1ffaJpG2LLdKNTWqsxsmbmW1L7C7tYWZjqr3cBriVWAcn8zN1KjzySDzWlKmIiCSlshKeeSamR1/JtBVv3BiOOgrOOw/WXz/d+CRv+Yy8TQH6AV2Abaq9Ngfom8e967+BA2HGDNhuO1hrrbSjERGRYvfrr3DffXDLLTB6dFxr2RKOOy6ayHfunGZ0kqAaJ2/uPtLMugEHubsStSX14INxPuqohb9PRERkYb75Jmq0lZXBlClxrVMnOO00OOGE6EMq9UqNkjczOyzz0IFXMtc2BMqIXqbPAT3d/Zckgqx3JkyIoexmzWD//dOORkREio17bHq77TZ47LGq9WxbbRWjbPvuC00SaaIkBaim/2f7A/cDl7j7N5kepkOI6dLTgVWAm4FjE4myvnn44fiL99e/QuvWaUcjIiLFYtq0mLm54w4YmWlk1LgxHHxwJG3ahNAg1DR5G+HuuYnZNUBr4M/uPhLAzO7NN7h6Kztlevjh6cYhIiLFYeTIqM32wAPRygqgXbvY8Nazp5rFNzA1Td7ezz4wsy2AY4Bbs4lbRuM84qq/Pv0UPvww1iCotpuIiCzI7Nmxua1PH3jrrarr224bCdt++0Hz5unFJ6mpafK2gpltAMwgplDHAxdlXzSzlYE98g+vHsqOuh1wACy1VLqxiIhI4fnkkyim269fVb/RVq1ig9vJJ8MGG6Qbn6SupsnbP4AngbWBz4G93X26ma0GHAicSkyjSi53GDAgHmvKVEREsn75JdZD33MPvP9+1fWNN45RtsMPj7IfItQweXP3z4B1zayNu/+c89I3RHHeW5MIrt75+GP4+mtYccWo7yYiIg1XZSW8+irce29Mj86aFdeXXRYOPRSOPRa6dlXrKvmDvPYRV0vccPdKYHZeEdVn//1vnPfZJ3YHiYhIwzNmTGw8uPdeGDu26voOO0RB3X33hWWWSS08KXwqAlOXssnbfvulG4eIiNStSZPg0Udj3fPbb1dd79QJjjkmji5d0opOV67XHgAAGB9JREFUioySt7ry+eexCLV1a+jePe1oRESkts2YAU89FQnb889DRUVcX2aZmIE59tgYbWvUKN04pegoeasrTzwR5z33jM4KIiJS/1RUwMsvQ//+8XM/W5OtcWPYbbfYeLD33tp8IHlR8lZXslOmaoclIlK/VFTExoPHHouErby86rUttoiE7eCDoX371EKU+kXJW12YMAGGDo2h8p13TjsaERHJV0VF9KjOJmzZemwAa68dCdthh8Gaa6YXo9RbSt7qwvPPx3nHHWHppdONRUREambu3N8nbJMnV722zjpw4IFxbLihyntIrVLyVheefTbOu+2WbhwiIrJkpk2DF16AQYPgmWfg55wKWdmE7aCDouuBEjapI0reatvcufDSS/FYyZuISOGbOBGefjoStiFDYM6cqtfWXbdqhE0Jm6REyVtte/ttmDo1/sJ37px2NCIiUp17lHIaNChKe+S2pzKDrbeOHaJ77x2jbSIpK5rkzcz2A24AlgX6A2e7e4WZtQWuAqYATYHz3X1uepFW89xzcdaom4hI4Zg9G15/PX5GP/UUfPVV1WtLLQU9ekSy9te/RktDkQJSFMmb/X979x4tV1mfcfz75B4ICQmEhhBCICThHoQEIgiJkgooVheoVGtdVGuKZelyKRZBq9SqVbT1RhUiClbLpVbbKgotKAQEwQQCQW7hkgQTEjghQEJuh5Pz6x/v3pzJnHvOnpmzT57PWu/a19n7PW8mM8+8+yZNBt4BvBM4Evgu6TmqXwN+AnwqIu6VdD7weeDiRtW1HYc3M7P+YcWK9Jl8003pXmxbtrQt23ffFNTe/vYU3Pbcs3H1NOtGKcIbcBDwgYhoAe6TdBTwRkm/BWZExL3Zer8AHpf0hYjY3KjKvqapCZYtS7/iTjml0bUxM9u9bNvW1rt2003pSTeVZs5MP6zf8pZ0aNTPnLaSKEV4i4g7q2atAfYC3gSsqFhvjdLJo7OARXWrYGcWZVU46aQU4MzMrHYiYPnydJHYTTel23pU9q6NHp3utXnmmXDGGTBxYuPqatYHpQhvHZgNfBVYALxQtWwT0D/+R95+exr6WaZmZrWxdm26IvTWW9Nw9eqdl+e9a2eeCa9/PQwd2ph6mhWodOFN0qFAU0Q8KCmArVWrDAfaXbDQ1NTErFmzXptesGABCxYsqGldHd7MzAq2cWM6qpGHtYcf3nn5+PHpYe9vfrN716x0Fi5cyMKFC/PJfTtbTxFRnxoVQNJg4MvAZyJiu6SLgTdGxJuz5QK2ZfPurnztrFmzYsmSJfWr7PPPpyuURo5MN3UcPrx++zYzGyg2b4Z77kk/hn/963Qbjx072pbvsQfMnZueYDN/fnq6waBBDauuWVEk3RcRszpaVraetwuBr0XE9mz6FuD8iuVTgGagjimtE3fckYYnneTgZmbWUxs3wl13pd61O+5Iz4VuaWlbPnhw+lydPz8FtjlzYNiwxtXXrAFKE94kXUIKZSMlHQLMA+4A1ko6KiL+AJwBfDsimjvfUp3cdVca+ipTM7PObdgAd96ZgtqiRbB0KbS2ti0fNAiOPx5OPTUdDj311HThgdlurBThTdJngH+smv1YRPxA0rnAZyWtBEYBn653/Tp0b3b3kjlzGlsPM7P+ZM2a9OSZPKw99NDOy4cMgRNPTCFt7tzUyzZmTGPqatZPleqct76o6zlvzc3pl+H27elX5dix9dmvmVl/8uqr8OCDKazl5Y9/3Hmd4cNTWJs7N5U5c3yDXDMG1jlv5bBsWQpu06c7uJnZ7qOpCX73u1Tuvjudr7a16oYAo0enW3acfHIKayec4PtgmvWSw1st5IdMTzyxsfUwM6uV5ub0Q3Xx4vSZd/fd8MQT7debPj0d+szL4Yf7alCzPnJ4q4Xf/z4NHd7MbCBobU2Pllq8OH2+LV4MDzyQAlylkSNTT1oe1ObMSc8MNbNCObzVwtKlaTirw0PVZmb9V0Q6L60yqC1ZAps2tV93+vQU1k44IR0KnTnTTzAwqwOHt6I1N8Ojj4IERx3V6NqYmXUuD2pLl6Zy330psD3/fPt1J02C2bNTOeGEdPuOvfeuf53NzOGtcI8/nm4oeeihvmLKzPqPHTvSOWl5UMvLC9WPhyZdaFUZ1GbPhv33r3+dzaxDDm9FW7YsDY85prH1MLPd1/bt6ZmfeUC7//702bR5c/t1x42D170uleOOS0Ft6tR09MDM+iWHt6LlN5w8+ujG1sPMdg/PPZc+d5YtS8OlS+GRR9I91qpNmtQW1PKwduCBDmpmJePwVjSHNzOrhS1bUm/aQw/tHNaamtqvK6WLCY47ri2oHXssjB9f/3qbWeEc3oqW3+fo8MMbWw8zK6cdO2DFirZwlg+ffDJdYFBt9Oh0cdQxx6QfjTNnpjJqVP3rbmZ14fBWpJaW9KELcPDBja2LmfVvr74KTz2VDnE+8ki6Sv2RR+Cxx2DbtvbrDxkCM2akgJYHtaOPhsmTfdjTbDfj8FakZ55JAW7SpHSzSjOzbdtg+fL2Ie2JJzo+Lw3SZ0gezvKwNmNGeg6ome32HN6KlB8yPfTQxtbDzOorAtavTyFt+fJ0y6A8pD39dHpCQUcOPhiOOCKdZnHEEakcdhiMGVPf+ptZqTi8FenJJ9PQ4c1sYNq8Of1Iy0NaHtSWL4eXXur4NYMHp4sH8nCWh7UZM3wvSDPbJQ5vRcrPdzvkkMbWw8x2XXMzrFqVQloezPKyenXnrxs9OoW0vOS9adOm+XCnmRXK4a1Izz6bhpMmNbYeZta1jRvTxQKV5emn0/CZZzo/zDl0aOpZrwxpM2ak4X77+cIBM6sLh7ci5eFt4sTG1sNsd9faCuvWtQ9oeUhbv77z10rpCs5p09qCWV4OOihd9Wlm1kD+FCqSw5tZfbS2wtq1sHJlOsS5alXbeD7s6HYbuREj0ukNU6e2Lwcd5MOcZtavObwVJcLhzawoLS3p/LKOQtmqVenQZme32cjts0/H4WzqVJgwAQYNqsdfYmZWOIe3omzalK5E22OPdOKymXXs1VfTD53Vqzsvzz7b+Xlnuf32S71kU6akYfW4/x+a2QDl8FaUDRvScN99fdKy7b62bYM1a7oOZs891/FjnipJcMABnQezyZPTDyUzs92Qw1tR8ns87b13Y+thVrSI9P5euzaVdet2HlaOd3avs0qDBqVTCyZN6rxMnAjDhtX+bzMzKyGHt6I4vFmZRMArr0BTUyrPP995KFu3DrZv79l2hwxJPWZdBbMJE3zFpplZH/gTtCh5ePNjbaxRNm9uC2I9GXZ1NWa10aNT6Np//7Zh5Xg+HDfOFwKYmdWYw1tRXn45Dd3zZn0VkW4i+8ILqWzY0Pkw7zlraoItW3q3n5Ej00n/48en0lk4mzDBj3EyM+tHHN6K8soraThqVGPrYf1DROrZevnl1CubD196KYWurgLZhg2wY0fv9zliRApheSCrHlbPcyAzMyulARHeJI0ELgM2AWOBiyJiY10r0dychj7Juvwi0jlemzalUL5p087hq6Pxjubl74ldsdde6T5l48Z1Pdxnn7ZANmqUr3Q2M9sNDIjwBnwHuDkibpB0BnAl8J661iC/YajDW321tKTDhVu2wNatbeNbtrQFr+phT+a1tPS9bsOGwdix6TzIvfdOZcyYFLy6CmXjxqVnaJqZmXWg9OFN0kTgvcAF2axbgZ9KmhIRK+tWkbyXZXf80m1tTT1VlWXbtvbzuipbt7YPXz2Z7u4u+7tq2LDU+zVqVBpWB7CejI8YUZu6mZnZbq304Q2YB6yPiC0AEdEiaQUwF1hZt1rkIaJM4e3WW+Haa1Pdm5t7Pmxu3jl4FdFLtasGDUo3ax05Mg3zMnJkW/CqDGGV410tdw+qmZn1UwMhvB0AvFA1bxOw0wNGm5qamDVr1mvTCxYsYMGCBcXVooznvD36KFx9dTHbGj686zJiRPfLK8NXZ4GsenrYMJ/nZWZmA8LChQtZuHBhPrlvZ+sNhPAWwNaqecOBnY6njR8/niVLltSuFqefnu6FNXdu7fZRtNNOg6uuSgFo6NDeDSuD19ChDlBmZmZ9VNmxJGl9Z+sNhPC2hnSFaaW9gGfrWot581IpkyOOSMXMzMxKYyDcCv02YGJ2uxAkDQUmA4saWiszMzOzGih9eIuIdcCNwGnZrPnATyNiTeNqZWZmZlYbA+GwKcCHgcskzQQmAAVeiWBmZmbWfwyI8BYRLwIfanQ9zMzMzGqt9IdNzczMzHYnDm9mZmZmJeLwZmZmZlYiDm9mZmZmJeLwZmZmZlYiDm8FqngemdWY27o+3M7147auD7dz/bita8fhrUB+o9aP27o+3M7147auD7dz/bita8fhzczMzKxEFBGNrkNdSGoCVtV4N/sC62u8D0vc1vXhdq4ft3V9uJ3rx23dNwdFxPiOFuw24c3MzMxsIPBhUzMzM7MScXgzMzMzKxGHNzMzM7MSGdLoCpSBpJHAZcAmYCxwUURs7GC9M4G3kkLxbyLiP3uyzJK+trOkMcBVwOnAU8CHI+KeOlW/VIp4T1escy6prefVtNIlVFQ7Z9v5G+BF4HG/r3dW0Gf0Jdn8VmAP4O/DJ4W304u2Hgd8EhgdERdULfP3YV9FhEs3BbgaODcbPwO4roN1pgP3A4Oy6TuAo7pb5lJoO38JeA9wPPBL4Dlgz0b/Xf2x9LWtK9aZCCwDbm/039QfSxHtTPqC/DlwSKP/nv5aCvjsOAv4XsW6P8i359L7ts6Wzc7et9f09N/BpefFh027IWki8F7gF9msW4E/kzSlatWPAzdHRGs2/b/AJ3qwzCisnW+NiOsi4r5sW+OAI2pZ7zIqqK1zFwFX1Kam5VZgO/8H8PWIeLp2tS2vgtr5cGCvinW3AmNqUd8y60VbExGLgQc62Iy/Dwvg8Na9ecD6iNgCEBEtwApgbtV6b8rm5x7PXtvdMkvm0cd2jojf5DMj4mVgI7C6ZjUur3n0/T2NpA8C1wJbaljXMptHH9tZ0tuAacBxkn4l6TJJw2pc77KZR9/fzzeSQsg7JY0l3Z/sRzWsc1nNo2dtnWvtYJ6/Dwvgc966dwDwQtW8TaTDRV2tV7lOV8ssKaKdXyNpOulQ3toiKzlA9LmtJU0FxkXEvZIOr1VFS66I9/R7gd8B3wKuAZYCIp1LZEmf2zkiHpX0HtKPkTuAcyJia22qW2o9bevebMPfh7vAPW/dC1IXeqXhwKvdrFe5TlfLLCminStdgL/gOtOntpY0iNS+X69ZDQeGIt7TRwJ3RcSrEfEC8D3gL2pQ1zIr6rNjKPAO4FDgBklDC67nQNDTtu7NNvx9uAsc3rq3hnTCcKW9gGe7Wa9yna6WWVJEOwMgaT6p183nCHWsr219MrAAWCdpPfBt4ORs3NoU8Z4ews5HSJaRzuW0Nn1uZ0lzgdkRcQvpEOCxpHOzbGc9bevebMPfh7vA4a17twETs8ujyX6NTQYWVa13CzCjYnp69trulllSRDuTHcI7JCL+q7bVLbW+tvXvgcNIX3DHAp8FlmTj1qaI9/SyqmUtwKM1qW15FdHO7wKeBoiIZ4HPAKfUsM5l1dO27oq/Dwvg8NaNiFhHOpn1tGzWfOCnEbFG0oWS8jfhFaT7iyFJpJMyv9WDZUYx7SzpEFKP0P9JmiJppqQP1/PvKIO+tnVEbI+I1Xkh3XtsezZumYI+O74JnClpeDZ9Ev7s2ElB7fwA8LqKze4g/UixCr1o65yyUsnfhwXwg+l7ILv66DJgJTABuDgiXpF0H/BP0XaT2PcBx2UvWxQR/1OxjU6XWdKXdpa0P3AvcGDVZv88Im6oyx9QIkW8pyu2dR5wXvgmve0U+NlxCumqvJER8cU6/gml0Nd2zs7jvJR05fR60ufIlyJiez3/jjLoRVvPAv4FGA+cHxGLKrbh78M+cngzMzMzKxEfNjUzMzMrEYc3MzMzsxJxeDMzMzMrEYc3MzMzsxJxeDMzMzMrEYc3MzMzsxJxeDOzAUnS30n610bXw8ysaA5vZtYrkg6X9DNJV0l6UFJkN93sb34FXNfoSpiZFc3hzcx667+Bf4iIvyY9UujfKhdK+kgjKiVpP0nn5tMR8YeI+G0j6jKQNOrf08w65/BmZj0maRzpQdKbASKilfRg+siWHwN8oQH1Gk4KkSPrve+BLHv02dmNroeZ7czhzcx6LCI2kJ6xebOkOdm8VcDvswB1HjBa0pclzZU0UdLV2flnd0p6U74tSWdI+pKkz0raIWmppPdJ+oCkxySdKuk3kl6UdFb2miGSLpf0KUk3Svp4trm3AUcA75L0cUnTJF0j6cqK/b1B0jckfU3SLZKmZfM73V8lSbMl/VTS5yQtlLRJ0vWSRkianv2dV0r6vqSnsvmzJF0q6YbsUPOYbFvnSPqEpH+XdE82b5Skz0u6WNIfs/Y5UdJySZdm65wtabukeb3dZ1ey7X5T0j9Lujmry8HAO4GpFf+eP8z2d0FWx9skjcvaeoOk7/fuHWVmuyQiXFxcXHpcgCOBVUAr8ANgv4plU9LHymvTXwM+m41/EvhFNj4W2ATsm03/Erg2Gx9P6sm7mNSTdjFwV7bsLOA32fhsYFPFvm4HzsvGhwNXAtdk038CLAOGZNMXAo8Ag7vaX9XfPQS4EbgNmAYcT+qB/BQwAliYbfNo4EPAGOD6itffBXw1G19UMf+L2fAC4N3Z+CnAGdn4D4FLK9ZfBczr7T67+TddD0zOxh8GzsnGzwNur2rT35J6X0cDzwHfBcYB+wPNwPRGv0ddXAZ6cc+bmfVKRDxMCguXA+8HHpY0s5PVFwLfl7Q/qWdsVDZ/Wja+NZt+lLZDnuuz4XURsRVYDEzI5t0JXChpD+D1FdurruN2YF3FrPcByyOiJZu+CjgceEM3+6vcZku27qKIeCIi7gOuBc6KiG3AWmBxRDwUEd8D3gLsI+ljkj6W/Y25AyVdIWk/4JvZvFeAyyS9jRSQ7s93XV2VrD693WdX3hURz0g6mRTS2rVrRZs+GRHLI2IjsBy4NyI2RMTabHm7tjOzYg1pdAXMrDwkDQYOjIiVwEezw2S/IoW0Ezt4yWrgIuAJ4B7g4Gz+H4BngNOAn2fzb4DUbSepchstpB4yIuJlSScCc4E7uqluZeiZBgx9bUHES5I2AJO62l8PPEwKkfn+Kvc5GVgREd/o4HXvAK4Bnib1SH4X+BFwCPATYAnwrh7svzf77MpLkr6S7Xs1oE7Wqw6SO6qmW/HpOGY15/9kZtYbQ6kIFRHxIHAJcFgn618BrIyIHwPbK163BfgEMF/S+4EfR8T13e1c0geAORHxdeCFXtT7GVLP306bI52/1xfDSL1PHVkLvFXSaxdRZMETYCMwi3SI9nJJU0mHLT+X1XM48JVs3VZ6Hia72meHJO1NCuCfj4glPdyPmTWQw5uZ9dZHJE2smD4AuDkbbwaQNFbSYaRbiYzPDnOeDIyUdEg2/QXSYcf7gGWShmWvzT+XKnt/8vHXkQ4LDiH12iFpqqQR2b7z/eavyV/3Q+AASXOz10wHnoiIJd3sryP7V4zPA75TMV0Zsn4F7AncKOl0SR+iLeReGBGtEfFtYGk271xJkyLiaeDzFXV4DpidXYxwNulcs/2yXtAe71PpIo+O7sc3lXQO4j6SZgAHAntImkJbm47IpivbNNfdtJkVzOHNzHrrQNJ5bj+WdA0wAzg/W7YWuIV08UAr8C3g06Sb5d5MOjx6NOlw24uksPEg8CSwNjt37oPZtv5S0iTSIcYJks4ghbDjSeelrQDWAH+Vnf91A6kn65QsaLwJOEHSkRGxhnTLiy9L+jLwUeCcbD9d7a8jR2ZXiv4jcHNE3JqFnj8F5uUBMSLWA28nXSxxPSlE/Sjbxt9mV25eROp1fIp0zt8iSZdkdf9Mtu7lpAsEHgJeAlYCx5J67nqzz3eTQmG1B4BFWZueDdxEulBhC+nfcUy2rRGkwDxb0jGS3kDqJTxL0kGS3k0KtudK2quTtjOzAiii+hQGM7PaynrH5kfE5dn0YNJVrKdlh0T7pSysroyISxtclV0i6ZMR8dVG18PM+sY9b2bWCJcAB0vaM5seDMwnXbzQn3V02LAUJJ1K2xWsZlZiDm9m1ghfIl0BukrSCuBnwG3Z4cN+SdIpwBzgNEnVFz+UwVMR8etGV8LM+s6HTc3MzMxKxD1vZmZmZiXi8GZmZmZWIg5vZmZmZiXi8GZmZmZWIg5vZmZmZiXi8GZmZmZWIv8PkZuigPS4gh0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(10,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/adept-vital-venus-thermal.png',bbox_inches='tight')\n", "plt.savefig('../plots/adept-vital-venus.thermal.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/adept-vital-venus.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 }