{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 34 - Genesis - Earth" ] }, { "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 Genesis sample return aeroshell.\n", "https://en.wikipedia.org/wiki/Genesis_(spacecraft)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \n", "planet.h_skip = 125.0E3\n", "planet.loadAtmosphereModel('../atmdata/Earth/earth-gram-avg.dat', 0 , 1 , 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Genesis', 210, 80, 0.00, 1.78, 0.0, 0.43, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(125.0,0.0,0.0,10.8,0.0,-8.0,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (30*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+17YcXAAAgAElEQVR4nOzdeZwcZbX/8c+ZyTYkIQzZ98lKIIEkEBLAyyarIoteROGngAKDKIgLKBcBl8hVFFEuroMXvYooshmJgrIvApGEBMhCQshGFpJJICF7MjPn98dTnemZzJqp7uru+b5fr3pVd1V31anp5MyZp596HnN3RERERESkbYqSDkBEREREpBCosBYRERERiYEKaxERERGRGKiwFhERERGJgQprEREREZEYqLAWEREREYlBh6QDiEuvXr28rKysxa+vrKykd+/emQsoiwrpWqCwrqeQrgUK63py6VpmzZq13t1zI5gsaW3OTsmlzy0uhXhNoOvKJ4V4TZC562oqZxdMYV1WVsbMmTNb/PpJkya16vW5rJCuBQrregrpWqCwrieXrsXMlicdQ7a1Nmen5NLnFpdCvCbQdeWTQrwmyNx1NZWz1RVERERERCQGKqxFRERERGLQbgvr8vLypEOITSFdCxTW9RTStUBhXU8hXUt7UoifWyFeE+i68kkhXhMkc13m7lk/aSZMmjTJC7F/kIgUPjOb5e6Tko4jm5SzRSRfNZWz222LtYiIiIhInApmVBARkaS4O9XV1QB06KC0mgnuEP2I0Y9YRHKVWqxFRNrovvvuo2PHjpx//vlJh1Kw/vUv6NgRjj8+6UhERBqnwlpEpI2Ki4sB9rRaS/zMwrpAbgsSkQKlwlpEpI1UWGeeCmsRyQcqrEVE2kiFdeapsBaRfJD1wtrMTjGzGWZWlrbtY2a21Mw2mNntZtYhbV8vM/uVmd1iZreZWcdsxywi0pRCLqxzJWersBaRfJDVwtrM+gLdgMlp24YA5wDnAl8GLgW+lPa2+4C73P3rwCLgO1kLWESkBQq1sM6lnF0U/bZSYS0iuSyrhbW7rwWm1ds8FPisu89y998BPwNOBDCzo4CD3H1G9NqHgavMrGu2YhYRaU6hFta5lLNTLdY1NW09kohI5mS9K4i719R7/py7V6VtWgWsiB5/EFia9tpV0cN2NUOZiOS2Qi2sIXdytrqCiEg+yMWbF48Efhk9HghsqLd/MzAgqxGJiDShkAvrFshKzlZhLSL5IKfmrzKzkUClu78abXJge72XdQZ2139vZWUlkybVNoqUl5dTXl6eqVBFRPbYl8K6oqKCioqK1NNe8UeVednM2SqsRSRJLc3ZOVNYm1kxcDlwXdrmVcDotNcY0BVYXf/9vXv3ZubMmZkOU0RkL/tSWKcXkma2PiOBZVC2c7YKaxFJUktzdi51BbkGuNXdd6Zteww4KO15GbALUAUtIjmjQ4fQRtHOuoJkNWdrVBARyQdJjGNtqYdp264HZgMlZjbczD5rZiPdfSawxszGRS89HbjD3XdlN2oRkcalWqyrqqqaeWX+yZWcrVFBRCQfZLUriJl1Bz4VPb3YzH5K+Cpxar2XvuHud0WPPwHcZGbLCOOpfiMbsYqItFSh3ryYSzlbXUFEJB9ktbB2983AL6Il5bvR0th7lgOXZDg0EZF9VqiFdS7lbBXWIpIPcqmPtYhIXirUwjqXqLAWkXygwlpEpI1UWGeeCmsRyQcqrEVE2kiFdeZpVBARyQcqrEVE2kiFdeZpVBARyQcqrEVE2kiFdeapK4iI5AMV1iIibaTCOvNUWItIPlBhLSLSRiqsM0+FtYjkAxXWIiJtpMI681RYi0g+UGEtItJGHTqEubZUWGdOalQQ3bwoIrlMhbWISBulWqyrqqoSjqRwqcVaRPKBCmsRkTZSV5DMU2EtIvlAhbWISBupsM48FdYikg9UWIuItJEK68xTYS0i+UCFtYhIGxVFd9a5O67KLyM0pbmI5AMV1iIibWRme4prtVpnhqY0F5F8oMJaRCQG6g6SWR07hvXu3cnGISLSFBXWIiIxUGGdWV26hPWOHcnGISLSFBXWIiIxUGGdWZ07h/XOncnGISLSFBXWIiIxUGGdWR07hn7WVVWgH7GI5CoV1iIiMVBhnVlmarUWkdynwlpEJAYdOnQAVFhnkvpZi0iu65DtE5rZKcB3gU+4+7JoWwnwA2AzUAp83d3fb26fiEiuSLVYV1VVJRxJvHIpZ6darFVYi0iuymqLtZn1BboBk+vt+jnwvLtfD0wDftXCfSIiOaEQu4LkWs5OtVirK4iI5KqsFtbuvpaQaPcwswHABcDD0abHgbPMrKypfVkJWESkhQqxsM61nL3ffmG9dWscRxMRiV/Wu4K4e42lptAKTgDWu/u2aH+VmS0Fjgd2N7FvWTbjFhFpSiEW1pBbOfuAA8L6vffaeiQRkczIemHdgIHAhnrbNgMDgKom9tVRWVnJpEmT9jwvLy+nvLw83khFRBrR2sK6oqKCioqK1NNemYkqIxLL2QceGNYqrEUk21qas3OhsHZge71tnQktHzSzb4/evXszc+bM+KMTEWmB1hbW6YWkma3PWGDxSyxnl5aGtQprEcm2lubsXCisVxHuHE/XHVhNSOCN7RMRyRmF2hWkAYnl7FSL9bvvxnE0EZH45cI41k8BA6IhmjCzjsAQ4Jlm9omI5Ix2VFgnlrNTLdYqrEUkV2W9sLbau2AMwN3fAaYDJ0XbTwYecPdVTe3LYsgiIs0q1MI6l3J2v35hvVrfWYpIjspqVxAz6w58Knp6sZn91N0rgSuAH5jZeKAfkH4HS1P7RERyQiEW1rmWs4cODesVK+I6oohIvLJaWLv7ZuAX0ZK+/T3gskbe0+g+EZFcUYiFda7l7CFDwnr58kwcXUSk7XKhj7WISN5LFdY1NTUJR1K40lus9WMWkVykwlpEJAapluqiIqXVTOnWLfSz3rkTli1LOhoRkb3pN4CISAxShXWq5VoyY/z4sH711WTjEBFpiAprEZEYqLDOjsMOC+vXXks2DhGRhqiwFhGJgQrr7Jg4MaxnzEg2DhGRhqiwFhGJQeqmRfWxzqzjjgvr556D3XtNlC4ikiz9BhARiYFarLNj4EA46CDYsgVmzkw6GhGRulRYi4jEQIV19px8clg//HCycYiI1KfCWkQkBiqss+c//zOs770X3JONRUQkXSyFtZl9OI7jiIjkq3wqrPM9Zx93XBjPeskSmDUr6WhERGq1qrA2s9PN7HkzW2xmS6JlGfBAZsITEckP27ZtA6CkpCThSGoVas4uLoZPfCI8vvPOZGMREUnXoZWv/z3wPWA2kJpQ1oAL4gxKRCTfbN26FYBu3bolHEkdBZuzr7gCbr8dfvc7+O//hp49k45IRKT1hfUSd7+t/kYzeyOmeERE8k5NTc2eFuv99tsv4WjqKNicfdBB8KEPwSOPwE9+AlOnJh2RiEjr+1hfYWbnNrD9rDiCERHJR+lFdY6NY13QOfuGG8L6tttg9epkYxERgdYX1h8Afm1m1WlLDfCLDMQmIpIXUt1AunbtmnAkeynonH3MMfDRj8K2bXDttUlHIyLS+sL6O0A5MAoYnrZ8K96wRETyx5YtW4CcLKwLPmf/8IdQUgL33APTpiUdjYi0d60trB8Bprn7EndfHi3LAN2XLSLt1vr16wHomXt30BV8zh4xAr73vfC4vBxWrkw2HhFp31p78+KfgP8ysyfTthnwMeDq2KISEckja9euBaBv374JR7KXdpGzr7oqzML4xBNh8phnnoEuXZKOSkTao9YW1jcSvlK8GEjNd2VAPwooSYuItEYOF9btImcXFcGf/gSTJsG//w0XXAB//jN0aO1vOBGRNmptV5AfAoPdvczdh0VLGaH1Q0SkXcrhwrrd5OxeveCvf4UDDoCHHoKLLoKqqqSjEpH2prV/z7/t7psb2N6vrYGY2cHAlcBiQgtLhbvPMbMS4AfAZqAU+Lq7v9/W84mIxGV1NNZbv35tToVxa1c5+7DD4NFH4eSTw82M778fWrJz755SESlUrW2x/nz9DWb2ceAnMcTye+C/3f3HwPeBP0bbfw487+7XA9OAX8VwLhGR2CxevBiAESNGJBzJXtpdzp4yBf7xDzjwQJg+HU48EVatymYEItKetbawPtDMbgQws/3M7NfAr4GXY4jlYKB79Hg70MPMBhCm3n042v44cJaZlcVwPhGRWCxatAiA0aNHJxzJXtplzj7mGHjhBRg2DF5+GSZMCC3ZIiKZ1qrC2t3PAB4zszuAOcBoYDzQ0MxerXU/YSKD7sCngKuAE4D17r4tOn8VsBQ4PobziYi02Y4dO1ixYgXFxcUMGzYs6XDqaM85+6CDYMYMOPVUWL8+TH9+1VWwuaGOMSIiMWmysDaz4+ovQCdgOfAm8E1gKHBTDLF8HthBaEnZ4u4PAAOBDfVetxkYUP/NlZWVTJo0ac9SUVERQ0giIk2bO3cu7s6oUaPo2LFji99XUVGxJ18BveKIRTm7rt694ZFH4OabwwghP/0pHHJIuMnRvfn3i4iktDRnN3fz4g8IX/dtBGoa2H8XYeim/sCX9y3UPUoIffQGAj8xs6WE4aG213tdZ2B3/Tf37t2bmTNntjEEEZHWeemllwCYPHlyq95XXl5OeXk5AGa2PqZwlLPrKSqC66+HM86Ayy4LXUPOPjv0vf7BD8IQfSIizWlpzm6usP4OMNPd1zX1IjOL42vFu4FPuvtGM3NCwv4S4a7ydN2B1TGcT0SkzWbMmAHAlClTEo4EUM5u1Pjx8OKLodX629+Gp56CI4+E884Lhff48UlHKCKFoMmuIO7+9+YSdPS6+9sShJn1Asa7+8Zo03eB9wlfXw6Ihm/CzDoCQ4Bn2nI+EZE4uDvPPvssAEcddVTC0ShnN6e4GK6+Gt56C669Fjp3DhPJTJgQWrSff15dRESkbVo7KkimvAvsMLOBads2AK8B04GTom0nAw+4uwZPEpHEzZ07lxUrVtCnTx8mTJiQdDjZlNc5u7Q0dAN580344hehpAT+/nc49lg44gi4807YsiXpKEUkH+VEYe3uNcA5wE1mVg7cAFwbTSpwBXC2mX0D+DBQnlykIiK1pk+fDsAZZ5xBUVFOpNOsKJScPXgw3H47rFgBN90EPXvC7NlQXg4DB8IXvhCmSFcrtoi0lHmBZIxJkya5bl4UkWxxd8aOHcuCBQuYNm0aZ5111j4fy8xmuXu7uo0uF3P2jh1w//3wy1/Cv/5Vu33UKPh//y8sI0cmF5+I5Iamcnabm1jMbKKZ3W1mX9XELSLSXrz44ossWLCAvn378qEPfSjpcFpMObtxXbrApz4V+lq/+ip86UvQr1/oMvKtb4UC+4gjYOpUmDtXLdkisrc4vrv8HHCyu/8IOCZ104qISCG74447ALjoootaNX51DlDOboHDDoMf/xjefhv++U+48ELo1g1eeSV0Gzn00FBoX3MNPPccVFUlHbGI5II2dwWJ7vru4+6rzKwI6Bb1s8uqXPxaUUQK04IFCxg7diwdOnTgrbfeYvDgwW06Xja7gihn77vt2+GJJ+AvfwmTzFRW1u7bf3846SQ47bQw22OOTcIpIjFqKmc3N451Qwf7EHCgu//BzCYTZtyaD3tuaMl6ghYRyaYbb7wRd+eSSy5pc1GdacrZ8SkpgY98JCzV1fDCC/DQQ/C3v8GiReHxQw+F144aFQrsE04Io4307Zto6CKSJa1qsTazm4GrgUfd/dxo2/eB19z9nsyE2DL52PohIvnn0Ucf5UMf+hBdu3ZlwYIFsRTWmWqxVs7OnmXLQpeRf/wjtGpv2lR3/0EHwXHHhSL7uONg6NBEwhSRGMTZYn0MMAi4KG3bI8DvgUSTtIhIpm3atIkrrrgCgG9961s531qNcnbWlJWFYfrKy0N/63//Gx5/PPS/fuEFWLgwLHfeGV4/eDAcdRRMnhyWww8PfbhFJL+1trB+LG362pQ4psYVEclp7s6ll17KsmXLOPzww7n66quTDqkllLMT0KEDHHNMWAB27w43PT73HDz7bFi//XZY7rsvvKaoCMaOrS20jzwSxo2D/LovVkRaW1gvMbPrgYFm9kng48BHgatij0xEJIfceuut3H///XTv3p177703X0YCUc7OAR07wpQpYbnmGqipgfnz4eWXQ8v2v/8Nr70Gr78elv/93/C+zp3D6CMTJtQuhx0G3bsnez0i0rhWFdbu/iczOwq4GLgQWAOc6u6PZyA2EZGccM899/C1r30NgLvuuouReTJLiHJ2bioqCq3R48bBZz4Ttm3fDnPm1Bba//43LF4MM2eGJcUsTFKTXmxPmAD9+4d9IpKsVo8K4u4vAS+lnptZsZmd5O5PxBqZiEgO+OMf/8iFF14IwI9+9CPOPTe/elIoZ+eHkhI4+uiwpGzaFCaqmTOndpk7N0xY8+abtd1IAHr3DoX62LG167FjobQ0+9ci0p41WVib2V+B5m6n6AWsB5SkRaSg3HnnnVx++eW4O9dffz1f+cpXkg6pScrZhaVHjzCCyHHH1W7btQsWLKhbbM+ZE8bUfuqpsKQbMKC2yE4V3IccEsbdFpH4Nddi/Vr0moVADXAZMBtIHyNpGHBgRqITEUlAVVUV1157LT/5yU8AuPnmm7n++usTjqpFlLMLXKdOMH58WC6KxnpxhxUrYN68sMydG9bz58Pq1WF57LG6xxk4MAwBWH8ZMgSKi7N/XSKFornC+oeEsa43ApjZCHe/qf6LzGxaJoITEcm21atXc+GFF/LEE0/QsWNHfv7zn3PppZcmHVZLKWe3Q2ZhXOyhQ+HDH67dXl0dxtdOL7bnzQst3qtWheXJJ+seq3PnMLnNQQfB6NF1i251KxFpXpOFtbvXG+KevmZW7O7VqQ1mdgIwOQOxiYhk1UMPPcRll13Ghg0b6N27Nw888ADHHnts0mG1mHK2pCsuhhEjwnLWWbXbq6pg+fLasbVTy6JFoXV77tyw1FdaWnu8+suAAeGmTJH2rrU3L94PvGpmzwObgYOB04C9WkRERPLF6tWr+fKXv8yf//xnAE477TR+85vf0L9//4QjazPlbNlLhw61BXF6CzfA5s2hwE4vtlPr997be5SSlC5dYNiwhovusrLQEi7SHrR2uL3HzOxU4PPAYcBG4Fx319eKIpJ3du3axc9+9jO++c1vsnnzZvbbbz++973vceWVV1JUAM1vytnSWt27wxFHhCWdO6xbB2+91fCybl3oYrJgQcPH7d8/FNhDh4Z1+uMhQ2C//TJ7XSLZsi/D7a0Gbkg9N7MiMzvC3WfFGpmISIZUV1fzxz/+kZtuuomlS5cCcPbZZ3P77bczdOjQhKOLl3K2xMEM+vYNS2pGyXSbN8OSJQ0X3StWwJo1YXnxxYaP36dPbaHd0FqT4ki+aFVhbWZ3NbC5N1ANnBNLRCIiGVJdXc2DDz7I1KlTef311wE4+OCDueWWWzjzzDMTji5+ytmSLd27145WUl9VVei7vWxZ6Nudvl62LBTe69aF5eWXGz5+aSkMHhyWQYP2fjxoUBgLXCRprW2xHkUY+9TTtg0CnostIhGRmG3bto3f/va3/OhHP2LJkiUADBkyhG9/+9t8+tOfprhwxxdTzpbEdegQunsMGdLw/pqa0JqdXmzXL8Dfey8sr73W+Hl69qwtuBsqwAcNUl9vybzWFtaXuvvC9A1mNoYwXa6ISE5ZvHgxd955J3fddRfr168HYMSIEXz1q1/lM5/5DF26dEk4woxTzpacV1QUxtUeOLDhbiY1NWECnJUr4e23w1L/8apVsGFDWObMafxcffrUbeUeODCMaJK+3n9/TQ8v+661Ny8ubGDzTuBS4Lo4AjKzo4D/ABYDzwNbgR8Q7mgvBb7u7u/HcS4RKTy7du1i2rRp/OpXv+KJJ2onFzzyyCP52te+xkc/+tFCbqGuQzlbCkFRUW3/7vo3VabU1MDatc0X36kuJ6+80vj5unbdu9iuv+7fX63f0rDW9rFeSt2vFIuBvsDv4wjGzD4DjHD39BttfgM86u73mtnpwK+A8+M4n4gUhpqaGp599lnuuece7r//ft577z0ASkpK+MQnPsHll1/OlClTsHbWDKWcLe1FUVEodvv3hyOPbPg11dXwzju1BffKlaHv96pVtetVq2DrVnjzzbA0pVevpovvAQOgd2+N793etLYryGPAPdQm6hpgrbsvamsgZnY88P+AU9K2DQAuAL4QbXoceMDMytx9WVvPKSL5q6amhpdffpn77ruPP/3pT6xatWrPvsMOO4xLL72UT3/60xxwwAEJRpk45WyRSHFxbZeTKVMafo07vP9+3WK7oeL7nXdg/fqwvPpq4+fs2DEU+821gGvUk8LR2sL6bnd/Nn2DmR1kZiPdfXEbY7kNeAG4w8xGAlOBocB6d98G4O5VUQvM8cCyNp5PRPLMjh07ePLJJ5k2bRoPP/wwa9as2bNv2LBhXHDBBZx//vmMHTs2wShzinK2SCuYQY8eYTnkkMZfV10dupQ0VHSnb3v33TDqyYoVTZ+3W7fmW7/794dOneK9Xolfawvr04Fn621bCzwFTNzXIMzsoOj9n3H318zsWuBR4LvAhnov3wwM2NdziUh+WbFiBf/85z959NFHefTRR9m6deuefUOGDOHss8/mggsuaJddPVpAOVskA4qLa7ueNNbvG2D79lBgN1aAp55v2VI722VTevduvvW7Vy91P0lSiwprM7uSkDC7m9nXG3jJY22MYyzwrrunBtL5KfBNwteX2+u9tjOwu/4BKisrmTRp0p7n5eXllJeXtzEsEcm2zZs38/TTT/PYY4/xz3/+k4X1ftNMmDCBs88+m3POOYfx48fnbTFdUVFBRUVF6mmvOI+tnC2SG0pKaqd2b4w7bNy4d2t3/WL8nXfC6CiVlep+koSW5mxz98b21X2h2WHARcD/1Nu1y93XNPCWFjOz84AKdz8gbdts4M+EFpHRadvfBL7p7vekH2PSpEk+c+bMtoQhIgnYvHkzL774Is8++yzPPvssL774IlVVVXv2d+/enQ9+8IOceuqpnHHGGQU3MyKAmc1y90nNv7JVx1TOFikg1dVh5JPG+n6nivF3323Z8bp1a9noJ+p+sremcnaLu4JEX/dd5+57tTyY2RB3b6YHUZNeA3qYWV93XxttqwJWAwPMrMTdt5tZR2AI8EwbziUiCdqwYQPPP//8nkJ69uzZVFdX79lfVFTE0UcfzSmnnMKpp57K5MmT6dixY4IR5yflbJHCUlwcCt4BA2BSE3+Gb98eJtxp7gbMLVtg0aKwNKWx7iepSX+GDAlFugRNFtZmdg3wd3efH236mpnVHwDWgBOBE/Y1CHd/w8weAc4FfmZmBwCdgD8AHwJOAqYDJwMPuPuqRg8mIjmjqqqK119/nRkzZuxZFixYUOc1xcXFTJ48meOOO45jjz2WY489ltLS0oQizm/K2SJSUgLDh4elMe6waVPTxffq1S3vflJaWrfQrr/07x/+MGgPmmuxHgakj1V1ONAVeIfa4ZsM6BdDLBcCt5tZJ0ILx/nRHeVXAD8ws/HRedQJTyQHuTtvv/12nSJ61qxZbN9et8ttly5dOOqoozj22GM57rjjOOqoo+im5o64KGeLSLPM4IADwtLUIEqNjX6SGgs8NeJJasr5xorvDh1CS3d6sV1WFor/ESPCTJgdWjucRo5qcR9rADM7HJjr7rvqbT8krYUkEeqvJ5I9NTU1LFmyhFdeeYXZs2fvWdatW7fXa0eOHMnkyZOZMmUKU6ZMYcKECXTWlGV1ZKKPdXRc5WwRySj30KKdKrIbWtaubfoYHTrA0KG1hXb6evjwMM18Lomlj3XkUHevMxFoNOzSCCDRJC0imbF7927mz59fp4CeM2cOmzdv3uu1paWldYroyZMn06tXrANeSOsoZ4tIRplBnz5haazv944doZU7VWgvXw7LlsFbb8GSJaEV/K23wvJYA2MW9eoFI0fCmDFw8MG162HDcq+lu6XD7R1L+PrwlGiw/3QHEu46fzjm2EQki2pqali+fDlz585l7ty5vP7668ydO5c33niD3bv3uv+N/v37M3HiRCZOnMjhhx/OxIkTKSsry9vh7wqJcraI5JIuXUJhPHJkw/u3b69baC9ZUvdxapbLl16q+75OnWDUqNpCe/x4mDgxtHIn9auopXX+RuB3hP57H6i3bxfw6ziDEpHMcXfWrVtXp3ieO3cu8+bNY8uWLQ2+Z8SIEXsV0X379s1y5NIKytkikjdKSkJhfPDBe++rqQmjnLz5JrzxBixYULt++22YNy8s6Xr0gAkTQpF9+OFwzDHZK7ZbVFi7++tmdgLwcXdXQhbJA9XV1SxdupQ33nhjr2XDhvqT4wX9+vVj3LhxjBs3jkMPPZRx48ZxyCGH6ObCPKOcLSKFoqgo3Pg4cCCccELdfakhAxcsCMX1nDkwe3YYzeSZZ8KSMmhQeP+pp8LZZ2eu33ZrxrHeRCOtHGb2H+7+fGxRiUiLbdmyhYULF+5VPC9atIhdu3Y1+J7999+/TvGcWtQfunAoZ4tIoevWLbRIH3543e1r1oQCe/ZsmDkTnnsu9PG+++6wdOkC550HU6eGEUri1Nw41n8FmmuqMsL0tn3iCkpE6tq5cydLlixh8eLFvPnmm3uWhQsXsnLlykbfN2jQIMaMGbNnOfjggxkzZgz9+/dXX+gCpJwtIhLGze7fHz784fC8pia0aD/5JDz4IDz7LPzud/DXv8Irr4SbIOPSXIv1QkJfvVXUjoFanwEfiS8kkfZp165dDRbPixcvZsWKFdTU1DT4vk6dOjFq1Ki9CujRo0fTvXv3LF+FJEw5W0SknqIiOPTQsFx9NcyYAUcdBRs3wj/+AZ/7XHznaq6wvgXY7O47m3qRmenucpEW2LlzJ8uWLatTPKceL1++vNHiuaioiGHDhjFq1ChGjRrFyJEjGTVqFKNHj2bYsGF0yLXxhiQpytkiIg3YsSO0VN9zD9x3X9jWo0focx2nJn8bu/v6Fh7no+gucxHcncrKSpYsWdLgsnLlShqblKmoqIiysrI9xXN6AT1s2DA6deqU5UCk+EwAACAASURBVKuRfKOcLSISbNsGs2aF1umnn4anngrbUs44A370o6anft8XbWrmMrMhwBeAK1GSlnZix44dLFu2rNHieevWrY2+t7i4mMGDB+8pmOsXz5qRUDJJOVtECo173WH3UqODvP56mJI93fjxcOaZcOGFYfzrTGh1YW1mxcBZwGXAKcBaQp8+kYLg7qxdu7bRwnnVqlVNvr+0tJThw4c3uAwePJiOHTtm6UpElLNFJP/t2hVmbGxo8pjFi8Owe/UVF4exrKdMCeNYn3JKuKEx01pcWJvZMOBS4DNAT2A5cCbwD+DojEQnkiFbtmxh6dKlLFmypMH19u3bG31vhw4dGDp0aIOF87BhwygtLc3ilYg0TDlbRPLBtm1hSvOVK8M6/XFq/c47YWSPxvTpA2PHwiGHhPWhh4bJYbp2zd51pDRbWJvZxwktHR8EFgO3Af8HXO/uj0YveyFjEYrsg6qqKt5+++1GC+fKysom33/ggQcyYsSIBovnQYMG6WZByVnK2SKStO3bYe3avZd16+o+f+edMDJHc4qKYOjQ0B86fRkxIqx79sz8NbVUS6qD7YQpcF8ALnX3RQBm1thQTiIZ5+6sX7++0cJ5xYoVVNfvXJWmc+fODBs2bE8rc/31/pmakkkk85SzRSQ27rB1K6xf3/hSWVm3YN68ueXH79SpdmbFQYMaXvfvD/nSi7LZwtrdpwPTo5teLjOzMmAasOcSzcy8saEORPbRtm3bWLp0aaNdNpq6SdDMGDRo0J5iuX7h3K9fP4qKirJ4NSLZoZwtIk3ZsaPpIrmhZWeTA3jurWPH0D2jXz/o27fppWdPKKT5ylozpfkK4EYz6wCcAwwys6nA/wKnAhWZCVEKVXV1NatWrWq01fmdd95p8v09evRosGgePnw4Q4cO1Qgb0q4pZ4sUtu3b4d13YcOGsG5oaWhf+pBzLdWlC/TuDb16Nbz07Fm3WC4tLaxiuTVa3VHU3auA+4H7zWwk8GVCfz4laanD3XnvvfcaLZyXL1/O7t27G31/x44dKSsra7TVWTcJijRPOVskd7mHQrepQrixfTt27Ns5O3VqvEBOL5TTn++3X7zXXcjadAeWuy8GrjazWTHFI3kmNaZzY9013n///Sbf379//0b7Og8YMIDi4uIsXYlI4VPOFskM9zDkW3OtxQ3t27Vr387ZqVMogA88sHZdf2loe9eu7bc1ORtiGdrA3X8Xx3EkN+3evZulS5eyaNEiFi1axMKFC1m0aBFvvvlms2M6d+vWrdHuGmVlZZSUlGTpKkQkRTlbpGGNtSDXXze0rYkvYJtUUtLyojh9e0mJCuRcpDHDBAjdNtasWbNX8bxo0SKWLFlCVVVVg+8rLi5m6NChjbY69+zZE9P/fBERybIdO1pXGKfWrb1RL6WkpOkW5MYKZbUvFRYV1u1MVVUVb731FvPmzWPevHnMnz9/TwG9paGpiwgjbAwdOpTRo0dz0EEHMXr0aEaPHs2oUaMYMmSIxnQWEZGM2bmz8WK4qQK5iXm+mtS5894FcnpR3NC+0lIVyBLkVEVkZp2Al4Gr3f1pMysBfgBsBkqBr7t70512BYCamhqWLl3K3Llz9xTR8+bN44033mBnI3+O9+zZs07hnCqkR4wYoS4bItIg5W1pjd27a8c9bmjZsGHvArmJkVWb1LFjywvj9H3qYiFtkVOFNXAtUJb2/OfAo+5+r5mdDvwKOD+JwHLZ1q1befXVV5k9ezazZ89mzpw5zJ8/v9FpuYcMGcLYsWP3LGPGjGHUqFH0zKWpi0QkXyhvt2M7dzZeJKcv69aFdUtm2auvuLhlhXH91+gmPUlCzhTWZnYMsAZ4L3o+ALgA+EL0kseBB8yszN2XJRJkDnj33Xf3FNCvvPIKs2fPZuHChTQ018PAgQP3FM/jxo1j7NixHHLIIXTv3j2ByEWk0ChvF6bq6lAEr1kTppyuv049XreudTPsQZiaumfPMCZyQ0vPnrVLqkDu3l0FsuSPnCiszawrcK67f8XMboo2nwCsd/dtEMZiNbOlwPHAskQCzbLq6mrmz5/Piy++uGdZuHDhXq/r0KEDY8eO5fDDD2fixIlMnDiRcePGccABByQQtYi0B8rb+ae6Okw3/fbbYVmzpuHiubISampadswOHcI4x40Vyn361H1eWhpaoEUKVU4U1sDXge/X2zYQ2FBv22ZgQEMHqKysZNKkSXuel5eXU15eHmeMGbdjxw5eeuklnnrqKV544QVmzJjB5nrNAZ07d95TPKcK6XHjxmmWQZE8U1FRQUXFnjlaeiUZyz5qU94uhJydS2pqQgtyqmheuXLvx6tXQyMDPO2ld+8wHXX//nXX6Y/79IEDDlBrsrQPLc3ZiRfWZnYaMNPd19Xb5UD9TsKdgQZHiuzduzczZ87MQISZU1VVxSuvvMITTzzBk08+yfPPP8+OelMpDRkyhKOPPpqjjz6aY445hvHjx9OpU6eEIhaRuKQXkma2PuFwWiWOvJ2POTtJ7uFGviVLYOnSsE5//PbbLRtHuU8fGDQIBg+GAQNCkVy/eO7TJ9z4JyK1WpqzEy+sgWuAiWljHZcC04DbosfpugOrsxda/DZu3MgjjzzCX//6Vx555BE2bdpUZ/9hhx3GiSeeyLHHHsvRRx/NgAENNtCLiCSpXeXtbHEP3TDeeAMWLIA336xbPDczkS09e4aCObWkCujU44EDoUuX7FyLSHuVC4X1BYQWjZQXga8AM4BrzazE3bebWUdgCPBMAjG2yZo1a7jvvvuYNm0azz77bJ3JVkaMGMFJJ53ESSedxAknnECfPn0SjFREpEUKPm9nUnV1KJTfeKN2WbAgrN97r/H3de8Ow4eHZdiw2sfDh4fieb/9sncNItKwxAtrd69Mf25m1UClu68ws+nAScB04GTgAXdveg7tHLFlyxYeeugh7r77bh5//HFqojtBiouLOfHEEznrrLM488wzGTFiRMKRioi0TqHm7UzYtg1efx3mzKldXnstbG/I/vvDwQfDmDEwejSMGFFbSPfsqf7MIrku8cK6GVcAPzCz8UA/IOfvbJk7dy533HEHd999N9uizNmxY0fOPPNMPv7xj/PhD3+Y0tL635SKiBSMvMvbcamuhvnz4YUX4MUXYcYMWLSo4RE2Bg2CQw4JBfSYMbXFdN++Kp5F8lnOFdbuXpb2+D3gsuSiaRl3Z/r06dx22208/fTTe7Z/4AMf4FOf+hTnnXceBx54YHIBiohkUD7m7Tjs3h2K5yeegH/9Kzyu3w+6uBgOPRQmTKhdxo8Prc8iUnhyrrDOJ+7O3//+d775zW8ya9YsALp27cpFF13ElVdeycEHH5xwhCIiEqfFi+GRR+Cxx+Dpp/eeIGXoUDj66Nrl0EN1w6BIe6LCeh+9/vrrfOELX+C5554DoF+/flx77bVccskl9OjRI+HoREQkLvPnw/33h+X11+vuGzMGTjkFjj8+FNIayEmkfVNh3Urbtm3jxhtv5Pbbb6e6uprevXtz3XXX8bnPfY79dEu2iEhB2LABfv97+PWvYd682u09esDpp8Npp8HJJ4fROEREUlRYt8Jrr73G+eefz/z58ykqKuLKK69k6tSpmjpcRKRAzJkDt94aWqd37gzbDjwQzjkH/vM/4aSTQBPdikhjVFi30D333MNnP/tZdu7cyZgxY7j77rs54ogjkg5LRERi8PLL8J3vwPTp4blZaJUuL4czz9RMhCLSMiqsm+Hu3HLLLfzXf/0XAJdccgm33347Xbt2TTgyERFpq7Vr4brr4Le/Dc9LSuDyy+Hqq6GsLMnIRCQfqbBuxs0338yNN96ImXHbbbfxpS99KemQREQkBg88AJdeChs3QqdOoZi+5hrQBLgisq9UWDfhF7/4BTfeeCNFRUX84Q9/4JOf/GTSIYmISBtVVcGXvww//Wl4ftppcMcdMGpUsnGJSP5TYd2IF154gS9+8YsAVFRUqKgWESkAO3fCBRfAgw+GVupbb4Urr9RshyISDxXWDdi0aROf+MQnqKqq4qtf/SqXXHJJ0iGJiEgbuYeuHw8+GIbNmz4d/uM/ko5KRAqJCusG3HDDDaxcuZIpU6bwve99L+lwREQkBv/zP3D33dC1Kzz5JBx+eNIRiUihKUo6gFzz2muv8bOf/Yzi4mIqKiroqDGWRETy3rJlcP314fH//Z+KahHJDBXW9UydOhV35/Of/zyHHXZY0uGIiEgMpk6FbdvgvPPCRC8iIpmgwjrNggULuP/+++ncuTPXXXdd0uGIiEgM3nknTE9uBjffnHQ0IlLIVFin+eUvfwnAxRdfzIABAxKORkRE4nDffbB7N3zkIzByZNLRiEghU2Ed2bFjB7///e8BuPzyyxOORkRE4vLgg2GtUVNFJNNUWEf+9re/8d577zFx4kQmTpyYdDgiIhKD3bvhpZfC49NPTzYWESl8KqwjDz30EADnn39+wpGIiEhc5s6FHTtCF5ADD0w6GhEpdCqsgV27djF9+nQAzjnnnISjERGRuMyeHdYaXk9EskGFNfDMM8+wadMmxo4dy6hRo5IOR0REYvLqq2E9YUKycYhI+5AzhbWZfczMlprZBjO73cw6RNt7mdmvzOwWM7vNzGKfseWf//wnAB/5yEfiPrSISEFKMme3xpw5Ya3CWkSyIScKazMbApwDnAt8GbgU+FK0+z7gLnf/OrAI+E7c53/iiScAOPnkk+M+tIhIwUk6Z7eUe22L9fjxSUUhIu1Jh6QDiAwFPuvuVcAsMxsHnGhmzwMHufuM6HUPAwvN7LvuvjWOE2/YsIE5c+bQuXNnPvCBD8RxSBGRQpdYzm6NFStg0ybo3Rv698/22UWkPcqJFmt3fy5K0CmrgBXAB4Glaa9bFT2cFNe5n3rqKdydY445hpKSkrgOKyJSsJLM2a3xyithPX58mHVRRCTTcqKwbsCRwC+BgcCGevs2A3tNi1hZWcmkSZP2LBUVFS060XPPPQfAiSee2KaARURao6KiYk++AnolHU8bZS1nt8bLL0fBHRn7oUWknWlpzs6VriB7mNlIoNLdXzUzB7bXe0lnYHf99/Xu3ZuZM2e2+nwvR5l3ypQprQ9WRGQflZeXU15eDoCZrU84nH2W7ZzdGiqsRSQuLc3ZOdVibWbFwOXAddGmVUBp2n4DugKr4zjf7t27mR0Nchr9BSIiIi2U7ZzdGu6QqtuV3kUkW3KqsAauAW51953R88eAg9L2lwG7gFiaOebPn8+OHTsYPnw4B2pKLhGR1spqzm6Nt96CjRuhb18YNCjbZxeR9ipnCmszux6YDZSY2XAz+yywEVgT3XEOcDpwh7vviuOc8+bNA2DixIlxHE5EpN1IIme3Rqq1+sgjdeOiiGRPTvSxNrMbgKn1Nr/h7neZ2SeAm8xsGdAN+EZc5128eDGAZlsUEWmFpHJ2a6T6V6sbiIhkU04U1u7+XeC7jexbDlySifOmCuuRI0dm4vAiIgUpqZzdGrNmhbUKaxHJppzpCpKEt956C4ARI0YkHImIiMRp5cqw1heSIpJN7bqwXr063Kg+ePDghCMREZE4rVsX1n36JBuHiLQv7bawdnfeeecdAPr27ZtwNCIiEpcdO2DzZujYEXr0SDoaEWlP2m1h/f7777Njxw66du1Kt27dkg5HRERisnlzWO+/v0YEEZHsareFdaq1ul+/fglHIiIicdq6Nay7dk02DhFpf1RYq7AWESkoKqxFJCnttrBevz5M896rV6+EIxERkTilCuv99ks2DhFpf9ptYf3ee+8BaCpzEZECoxZrEUlKuy2s3333XQBKS0sTjkREROK0bVtYq7AWkWzLiZkXkzB8+HDOOussxo8fn3QoIiISo1694Jxz4PDDk45ERNqbdltYn3vuuZx77rlJhyEiIjGbMgUeeijpKESkPWq3XUFEREREROKkwlpEREREJAYqrEVEREREYqDCWkREREQkBu22sK6oqEg6hNgU0rVAYV1PIV0LFNb1FNK1tCeF+LkV4jWBriufFOI1QTLXpcK6ABTStUBhXU8hXQsU1vUU0rW0J4X4uRXiNYGuK58U4jWBCmsRERERkbxl7p50DLEws0pgeSve0gtYn6Fwsq2QrgUK63oK6VqgsK4nl65lqLv3TjqIbNqHnJ2SS59bXArxmkDXlU8K8Zogc9fVaM4umMJaRERERCRJ6goiIiIiIhIDFdYiIiIiIjFod4W1mZWY2R1m9t9m9gsz2z/pmPaVmX3MzJaa2QYzu93MOiQdUxzMrJOZvWpmJyQdS1uZ2VFmdo2ZnWNmvZKOZ1+Z2cFm9jMz+7KZ/dzMJiQdU2uZ2SlmNsPMytK2FUw+KHSF/FmZ2XfNzKPl1aTj2VeF+H+soWuKtuftZ9ZY7WBmvczsV2Z2i5ndZmYdk461pZqqh8zs7rTPalqmYymIQqyVfg486u73mtnpwK+A8xOOqdXMbAhwDnAuMBb4BfA2cGuSccXkWqAs6SDaysw+A4xw9xuSjiUGvwfOdvdV0b+9fwAHJxxTi5lZX6AbMLneroLIB+1EQX5WZtaVcIPVKdGmfbmhM3GF+H+ssWvK58+smdrhPuA6d59hZp8DvgP8V1KxtlRT1xTtW0ftZ7Uw4/G0p5sXzWwAsBQodfdt0V80m4Cx7r4s0eBaycyOBV5096ro+Q8I13FGspG1jZkdA4wBbgIudvenk41o35jZ8cCNwCleAP/JzGwrcIS7v2FmvYFX3X1A0nG1hpkVAdXAMHdfVkj5oNAV8mdlZl+KHv7S3XckGkwbFeL/sfrXFG3L28+ssdoBmAo8mMrrZjaQUIT2dfetScXbEk3VQ2b2Y+BfwF9S+zOtvXUFOQFY7+7bAKIf8lLg+CSD2hfu/ly9fySrgBVJxROHqBXgXHe/K+lYYnAbsAC4w8weNbMPJB1QG90P/NrMugOfAq5KOJ5Wc/eaeptOoEDyQTtwAgX4WZmZAZ8EfgSsNbO8aMltTCH+H6t/Tfn+mTVRO3yQ8NmkXrcqejgpi+Htk8auycz2A84E/hw9Pykb8bS3wnogsKHets1AXrW8NeJI4JdJB9FGXwe+n3QQbWVmBwETgTvd/UrgCeDR6GvFfPV5YAfwMrDF3R9IOJ44FHI+KDQF+Vl5cBShW8EdwB/M7EMJhxWngvvcCvAzS9UOhfRZHUn4NmGbu48kXNvfgEfM7NBMn7y9FdYObK+3rTOwO4FYYmNmI4FKd8+rGyjSmdlpwEx3X5d0LDEYC7zr7q9Fz38KGHB2ciG1WQnwR+Ae4CdmdnLC8cShIPNBgSroz8rd34vuxfge8KXmXp9HCvZzK4TPrF7tUBCfVUP1kLuvcffLCH3Ir8x0DO2tsF4FlNbb1h1YnUAssTCzYuBy4LqkY2mja4C7zGy9ma0HBgPTzOxrCce1LzqQdmOwu28H3gR6JhZR290NPODu3yF8q/DHqOtOPiu4fFDA2stn9T/AkKSDiFF7+Nzy8jNroHao81lFXV66kkefVQvqoax8Vu2tsH4KGGBmJQDRUDJDgGcSjaptrgFudfedSQfSRhcAE9KW1cCl5Gf3lteAHvW6flQB8xOKp00sDBM43t03Rpu+C7xPHo0K0ohCzAeFqr18VjXAK0kHEaP28Lnl62dWv3Z4DDgobX8ZsAuYmeW42qK5eigrn1W7Kqzd/R1gOpDqwH4yoRVuVePvyl1mdj0wGygxs+Fm9tnoa5C84+6V7r4ytRDuwq509/eTjq213P0N4BHC0D+Y2QFAJ0Ifr3z0LrAjuks8ZQOwKKF49knUAgOhW07B5YNCVqiflZkNMLMLzKwo+vf5VeAbSce1rwrx/1j9ayqEz6yh2gHYCKwxs3HRy04H7nD3XUnF2RqN1EPjzOysaH8n4CLgBxmPpQBGAmsVMysl/GCXAf2A/3L3LYkGtQ/M7AbC8Djp3nD3fG9FBMDMlpHfw+31Am4n/LU/hHAjY162WAOY2XjCDYyzgL7As+6eN61OaaOZ/JwwNutP3b2yUPJBe1CIn5WZjSX8Ef4+8Dzwk+gP87xTiP/HGromoA95/Jk1VTuY2VDCULfLCON3fyNbQ9S1RWPXROhP/UfCaCcvAbe4e8a7trS7wlpEREREJBPaVVcQEREREZFMUWEtIiIiIhIDFdYiIiIiIjFQYS0iIiIiEgMV1iIiIiIiMVBhLSIiIiISAxXWIo0ws5PM7AUzO6GZ140zsz+b2Y0xnvuHcR1LRCQJZnaYmd0XZ25s5DzTzeyTmTxHppjZyWb2YnO/ZyR/dEg6AJHmmNnBwM2EGQCPBA4DPu3ud2f41G8D41vwuhVAKVAcx0mjaX93x3EsEREAMzsS+B5wHPBPwqQZ/YCOwC/c/R8ZOO1ioAcx5cYm/ByYm4kDR/n4Y8AvCbPNvhjt6gz8JzDZ3Ze14RQt/T0jeUKFteSDvwDnufurZlYE/CYbJ3X3RWa2oQWve9/M1sR46uOBZ2M8noi0c+7+spn9CRjj7h9JbTezM4CHzOwad/95zOfcFnNuBMDMzgReTxW07v73uM+R4u67gXvN7CrgaXe/IS2OaftyTDPrA5zo7ve6+0IzWx9TuJID1BVEcpqZHQiMBrYCuHsNYcrVbE0ZWhPz61riROCpGI8nIgKw1/TU7v434FvAj82sXwbOGWuuNrNBhBbqbGtoau/HgVWtOYiZdQZ+B5SkbY7z94ckTIW15DR3fxdYCDxqZkdF25YD/zazA8xsqpn928zOMrN3zGyBmR2eer+ZnW9m3zCzp83sv9O2TzKzb5nZvWb2oJn1iLZ3MbOfmdlNZva/hK8xW8XMvhbFdJ6ZHWlmf4mONzWK8V4z6x1t32Rm36x3iE7uvjM61lfN7Ctm9i8z+36rf4AiIs37X0KXkI8DmNkIM/uOmf3azJ6IilnMrHOUy75hZs+Z2SdSBzCzz5vZjWb2kpld2diJGsrJZnZulDMvNrM3zez7ZjbAzH4T5dPnzOyD0SEuAAYBXzezc6J7YR43s4vSznGZmd1sZr8zs/8zs65m1sPMvmlmM83sg2a20MyWmFnZvvzAop/JCVGLNmb2MTO73cx+ZGaPmlm36BpuN7OHo983a6Kf8SHAx83sK2mHLI1+F202s8saOWeRmd1iZuXROdzMHkp9PpIj3F2LlpxegLHAcsJf9XcBfaLtxYQk9S6hr9sBwGPAPMIfjUcDN0SvLSW0ep8B7A/8Ke34/wJ+GD3+MVAePe4LVBOSZ3Mx/pbQ6tOH0Bdvv2h7Z+AfwDRgAOEXwm7gm8B+wBHR89TrU/3HU9f98+hxd+C6pD8LLVq05O8CXAysbGTfWuB/orz6IFAcbb8HuC96fBvwsejxZ4HK6PH5wKeix6OjXH1o9Py3wLeix43l5AOBjcBPgMnAycCtwE3Ra68FHk6L1YGy6HEJsAC4OHp+GnBP2munA78idH39aPT74mPRdT4K3NzCn93T0e+K7wM/JPS3Pjlt/3pgSPR4HuF3UmfgG8BKYApwSZT3n07FG71+GfCjaN/5wIpGYrgUmB497hv9nE9N+t+VlrqL+lhLznP3eWZ2KPBd4PPAmWZ2soc+15XA++7+AICZfQt4HhhO+CVSY2Zfig71D0ICPwPombZ9QfTebsAV0Xtx97Vm9nYrQu0F3AFc4u7bomPsNLN3gKXuvjo6zzvAMx76H84mJPxehJsgzyD8EgDYBlxoZouAXxD+qBARyQQnFJuTgaHAVWYGsAHoYeHJ5YRGAYD/I3SFgJBrZ9bLtf2B1+ud42IayMnu/q6ZbQKmufu/AcxsBbDVzPoTWni7NRi0+3YzW5e26crouCm/Bv4cbd9I+H3xYHSOmVGcLfWUR32szexeQmNOysfdfYWZfYBQUHeL8v9q4C13nwHMiN7b0LEfjn4nvNxETIcD26PrXmtm71K3S4nkABXWktPMrBgY7OEmlS9G3TP+DlQQWgDqmxetDwCGAL9393uibT+Jjvl1QqH7k3rnOoKQEDfvY7ijgf8ARgCvpm2v38ewes8O95ooyaa6ZfX00P0Fd19qZhcDPwW+DPw/YB0iIjEyswMI37bNI+TNDQ3kx76EFtVOAO5eTWgMIHrPt939heh5nfemaTAnR5y6uXIl8HXgTeAlYFgTl5D+vlHAk2nP3yJ0c+nD3rm4in0fseQVwh8gKRvN7BbgPkLsqeq5/nU1p4rGa7NHgF+YWRfCZ7Gd0PotOUR9rCXX7en3B+DurwLXA2MaeX0nQhJbDKwBzk3faWZTou1nmFlJve3vR08P3sdYXwBuB/4QJb5WMbMBpN0IY2YDCSOijCbczHhPI28VEWmLiwjfkN1HyI8fiFqKgT35cT2wCzg9bfuhZtadernWzDpEDRX1NZaTG/JLYJmHYVV3tuJaVhBauPecAthE6OoSG3eviRo/xkZ/mPwd+I67z4zzPPXO+TDwMOGb1fOB4919U6bOJ/tGhbXkg6uiojNlIKFvXMqBZtYpenwCoX/dRuCPwEfN7KdmdkLUmrCVkAC7AtPN7LToRpExhGJ8NjDVzDqZWS9C3+Y+ZlZsZteZ2acaidGi5SbCL5/bGthX//X1n3+E0B8wZRRwjru/D1zVwHtERFpjr5ZQMzsO+DZwubtXElqH1wJ/N7MzoxsUj41aqP9EGD3kHDM7CbjI3TcTcu0XzezbZnY8oYEh1Y0uPf81lpNT0luPJwK9zWw/4ANAiZkNj/btJtzsl2pgST/HL4H/NLPS6PkU4E53ryLUPA3mYjM7KoqrcyM/uyLq1UzRHx9nEL6lLCV0MTwIGAzsl3ZjZP1W8V0NxF83qAb6i0Sfxf7AE4RW+db8wSFZoq4gkg8GA/PM7G+Er8kM+Fza/irgG2a2kdC6exWAuz9mZl8kfJ14DuHmv7kAZnY28DPCL4q7gGvd3c3sPMJQSPOAewlf6U0gTKhwHqGLSZ2JacxsMuGmnCGErwb/Anw7iucf/1K+2wAAIABJREFUwFFAmZmNILSk9APOM7M3qG29uRjo6u4V9a79NxYmdqgi3CwkItJqUR65gFD83QlUEm6A6wF8yN1fBHD3XVF+vJOQ66YR+lYDXE1olPgd4Ru61Egc/0vIfVcQuqxd4e7rosLxKGCEmd3dWE42s3MJ/YovNbM33H0V4UbKWwlF9e+As4FDgSXA7wmt6xdaGC1qLPBhM/u7uz9oZqOAB8zsX4RvPb8Rtax/HOgXXd88wmQ5/c1sbBTnJwn9xv+S9nPrSMj944GBUeHthL7NH4uOOQd4BniZ0L3lEUJOn064YXK8mZ3j7qnj3kuYrGdLVIj3A841s6WElmii99efs2EF4ebMc4AuQJGZPU64gTFbQ9BKM0yfheQzC9PA/tbdy7J0vmvdXdONi4gUmKhbSmd3z8kJuszsWsIILcui510I3zbckro3R5KnriCS7xrqZpGZE4WvTF/JxrlERCR7opszD8/horoLcANQlra5F7BVRXVuUVcQyVtR/7ZPEb7a+7i735fhU74VfUUpIiKFpYO7/yLpIBrj7jssTKv+P9GN7asJ3VZuaPqdkm3qCiIiIiIiEgN1BRERERERiYEKaxERERGRGKiwFhERERGJQcHcvNirVy8vKytr9fsqKyvp3bt3/AHlIf0saulnUUs/i1qZ+lnMmjVrvbu3qx9ye8/Zuo7couvILbl+HU3l7IIprMvKypg5s/UziU6aNGmf3leI9LOopZ9FLf0samXqZ2Fmy2M/aI5r7zlb15FbdB25Jdevo6mcra4gIiIiIiIxUGEtIiIiIhKDdl9Yl5eXJx1CztDPopZ+FrX0s6iln0XyCuUz0HXkFl1Hbsnn6yiYCWImTZrkudwfR0SkMWY2y90nJR1HNilni0i+aipnt/sWaxERERGROKiwFhERERGJgQprgMsvh0MOgddeSzoSkf/P3p3Hx1XX+x9/fbI2e9qmO7SlpWUpBQplFSyyKKCyXcQLV6+IUvcfetWLV3H3XkVxBRTqvu8LgogiOyJLS4EW6AJ0oS1d0jRNmqRZP78/vmeSyTR7ZjKTzPv5eJzHzJxz5pzPOZN88sl3vud7RKQH9SvqWb54OWvevibdoYiI9EqFNcDmzfD88+FRREQyTkdLB/tW7KPhuYZ0hyIi0isV1gAzZoTHrVvTG4eIiPQotzQXgPaG9jRHIiLSOxXWANOnh8dt29Ibh4iI9Ci3JBTWHQ0daY5ERKR3KqxBLdYiIhkuVli371OLtYhkrhEvrM3sHDN7zMxmx827xMw2mNluM/ummeXFLasys1vN7Hoz+5qZ5Sc9KLVYi4j0KFNytrqCiMhoMKKFtZlNAUqBE+PmzQQuAi4FPgS8E/hg3Nt+C/zA3a8F1gGfS3pgarEWETlAJuXsnKLw56qjqQPvGBs3NhORsWdEC2t33wHcljB7FnCVu69w958ANwOvATCzk4HD3P2xaN3bgQ+YWUlSA4sV1i+/nNTNioiMZpmUsy3HyCkOf7LaG9VqLSKZacS7grh7R8Lrh9y9LW7WViA27t2ZwIa4dWNNysm99e/kyVBaCnv2QE1NUjctIjKaZVLOVj9rEcl0mXjx4gnALdHzGcDuhOX1wPTEN+3atYvFixd3TsuWLRv4Hs3g0EPD8/XrhxCyiMjgLFu2rDNfAVXpjmcYRixnx/pZa2QQERlpA83Zeb0tSAczOxTY5e5PR7McaEpYrRBoTXzvpEmTWL58+dB3Pm8ePPVUKKxPOmno2xERGYClS5eydOlSAMysOs3hDMlI5+zOFmtdwCgiI2ygOTtjWqzNLBd4F/CxuNlbgfFx6xhQAiR/+I5588KjWqxFRPqVjpydUxL1sVZhLSIZKmMKa+AjwA3u3hw3727gsLjXs4EWYBhN072IFdbr1iV90yIiY9CI5+zOIffUx1pEMlQ6xrG22NO4eR8HVgJFZjbHzK4ys0PdfTnwipkdFa16LnCju7ckPbAjjwyPq1YlfdMiIqNVJuXsvLLQe7G9XoW1iGSmEe1jbWZlwFuil1ea2U2ErxI/n7DqGnf/QfT8zcCnzGwjYTzVT6QkuIULITcX1qyBpiYoKkrJbkRERotMy9l548OfrLY9bf2sKSKSHiNaWLt7PfCdaIr5QjT19p5NwDtSHFoopA8/HJ59NrRan3hi/+8RERnDMi1nxwrr1poDroUUEckImdTHOv2OOy48PvlkeuMQEZED5E8Id0dXi7WIZCoV1vFihfUTT6Q3DhEROYBarEUk06mwjnfaaeHxgQfSG4eIiBwgb4L6WItIZlNhHW/RIigvhxdfhJdfTnc0IiISJ3981BWkRoW1iGQmFdbxcnPh1a8Oz++7L72xiIhIN7EW69Y96goiIplJhXWic84Jj7ffnt44RESkm87h9tRiLSIZSoV1oosvDo933gmNjemNRUREOhVMLgCgZWcL7p7maEREDqTCOtHBB8NJJ4Wi+q9/TXc0IiISyS3LJac4h46GDt19UUQykgrrnlx2WXj84Q/TG4eIiHQyMwqnFwLQvK05zdGIiBxIhXVP/vM/oaAgtFhv3pzuaEREJFIwPeoOsq0lzZGIiBxIhXVPqqrgkkugowO+9710RyMiIhG1WItIJlNh3Zv3vCc83nwz1NenNxYREQGgYJparEUkc6mw7s3pp8OrXgU1NXDTTemORkREgMIZUYv1FrVYi0jmUWHdGzP4zGfC8xtuCAW2iIik1bg54wBoeqEpzZGIiBxIhXVfzjoLXvOaUFR/4hPpjkZEJOsVzSsCVFiLSGZSYd0XM7jxRsjLg1tvhcceS3dEIiJZrWhuKKz3b9hPR1tHmqMREelOhXV/FiyAD30I3MMwfA0N6Y5IRCRr5RblUnhQId7mNG9SP2sRySwqrAfis58NBfa6dfDhD6c7GhGRrBbrDtK4tjHNkYiIdKfCeiCKiuAXvwg3jbn11vBcRETSomRhCQD7nt6X5khERLpTYT1QRx8NX/taeP6Od8ATT6Q3HhGRLFV2XBkA+55UYS0imWXEC2szO8fMHjOz2XHziszsRjP7PzP7jpmVD2TZiHvve+Hqq2H/frjgAli/Pm2hiIiMhEzM2aWLSgGoX6mbd4lIZhnRwtrMpgClwIkJi74NPOzuHwduA24d4LKRZRZuFnPmmbB9exiOb9OmtIUjIpJKmZqzi48oxgqN/S/up7W2NdmbFxEZshEtrN19ByHRdjKz6cAVwO3RrH8AF5jZ7L6WjUjAPSkogNtug1NPhZdfhjPOgBdeSFs4IiKpkqk5Oyc/h7JFoTtI3b/qkrlpEZFhGfGuIO6eOPDoGUC1uzdGy9uADcCSfpalT2kp3HknnHgibNwYbn2+cmVaQxIRSYVMzdmVZ1QCUHtfbbI3LSIyZJlw8eIMYHfCvHpgej/Lutm1axeLFy/unJYtW5aSYDtVVMA//gFnnw07d8KSJXDPPandp4iMGcuWLevMV0BVuuMZhIzI2ZWviQrr+1VYi0jqDTRn541cSL1yIPHetIVArONcX8s6TZo0ieXLlyc/ur6UlcEdd4Qbx/zmN/C614WRQz7wgdAfW0SkF0uXLmXp0qUAmFl1msMZjIzI2eWnlmN5Rv2Kelr3tJI/Pn/I2xIR6c9Ac3YmtFhvBcYnzCsDtvWzLDMUFoZxra+9Ftrb4Zpr4KqrwsghIiJjT0bk7LzSPCpOr4AOqLmzJtmbFxEZkkworO8DpptZEYCZ5QMzgQf6WZY5cnPhS1+CX/0q3EzmRz8K/a9XrUp3ZCIiyZYxObvq4vBt7K4/7ErF5kVEBi0d41jH+kgYgLtvB+4Azormnw383t239rVsBEMeuDe/GR55BObNC0X14sXw9a9DR+K1PyIio0Mm5+yqi0JhXfPXGtob21OxCxGRQRnpcazLgHdHL680s0nR8/cAF5rZJ4DzgaVxb+trWeY59tgwQsjSpdDSAv/1X3DaabB6dbojExEZlEzP2eMOHkfZiWV0NHVQ/afR1E1dRMYqc/d0x5AUixcv9hG/eLE/t90G7353uJlMXh585CPwyU9CcXG6IxORDGJmK9x9cbrjGEnJytlbb9nK+vesp/I1lRx777FJiExEpG995exM6GM9dl14ITz/fLgVent76Id92GHwk5+oe4iISBJMuXwKOcU51N5XS+P6xnSHIyJZToV1qlVWws03h77XixbBli3wtreF/tf33Zfu6ERERrW8ijwmv3kyAFu/lZmX34hI9lBhPVJOPhmWL4cf/xhmzAj9sM88M4x9/c9/pjs6EZFR66D/OgiAV773Ci07W9IcjYhkMxXWIyknJ9xMZt06+N//DTeY+fvfw8WNZ58NDz6Y7ghFREad0qNKmfjGiXTs72DLN7ekOxwRyWIqrNOhuBg+/nHYuDFczFheHm6HvmRJmG67LfTJFhGRAZn58ZkAbPnGFpq3Nqc5GhHJViqs02nCBPjc50KB/ZnPhP7YDz4IF10E8+fDN78J9fXpjlJEJONVnFxB1b9V0dHYwUufeCnd4YhIllJhnQnGj4dPfxo2bQo3lDnkEHjpJfjgB+Ggg8Ljc8+lO0oRkYw29/q5WIGx48c7qHusLt3hiEgWUmGdScrLQxG9fj388Y+hW0hdXWi5XrAATj8dfvpTaGpKd6QiIhmnaG4RB//XwQCsuWoN7fvVpU5ERpYK60yUmxu6g9x/Pzz5JLzrXVBaCg8/HC5+nD4drrlGd3MUEUkw61OzKDqsiMbnGtn02U3pDkdEsowK60y3aBHccgu88gp897twwglQWwvf+hYsXBiG8bv5ZqjW7XxFRHKLcjn8h4dDDmy+fjM1d9ekOyQRySJJKazN7PxkbEf6UFoK73wnPP54aMV+z3vCcH2PPQbvfz9MmwYXXAC//S3s35/uaEUkg431nF1xSgWzPjkLHJ6/4nn2v6ycKCIjY1CFtZmda2YPm9kLZvZSNG0Efp+a8KRHixbBt78dWrF//nM47zxwh9tvh8sug6lT4e1vh7/8BZo17JRItsrmnD37k7MZ/7rxtFa38uy/PUt7o/pbi0jqmbsPfGWzXcAXgZVAR2w2cIW7L01+eAO3ePFiX758eTpDSK8dO+CXv4Sf/QxWrOiaX14Ob3gDXHppuMtjcXH6YhSRHpnZCndfnILtZnXObt3dyorFK9i/cT8T3ziRBX9YQE6eekCKyPD0lbMHW1g/5u4n9TB/qrtvH0aMw5b1hXW8NWvgd78L09NPd80vLobzzw+F9nnnweTJ6YtRRDqlsLDO+pzd8HwDK1+1krY9bUy7ehrzb52PmaV8vyIydvWVswf7r/t7zOzSHuZfMPiwJGUOPxyuuw6eeioM3Xf99eGix8bGUGxfeWXoLnLyyfD5z4c+24P4B0tERo2sz9klR5Sw8PaF5IzL4ZXvvsL6967HO5TvRCQ1BltYvwr4npm1x00dwHdSEJskw6GHwn//d7jocdMmuPFGOPdcKCgIFz5+6lNw/PEwYwZcfTX84Q+wZ0+6oxaR5FDOBipeVcFRfzoKKzS23bKNde9ep+JaRFJisF1B9gDvApYDsStBDHiru38++eENnLqCDFJDA9xzT7jA8S9/ga1bu5bl5IRi++yzw3TqqTBuXPpiFRnjUtgVRDk7Ts3dNay+YDUd+zuouqSKI356BLnFuSMag4iMfsnsY/0L4O3u3pwwP2v6641J7qEv9l/+AnffDY88Aq2tXcvHjQt3fTzrrFBoH3tsuImNiCRFCgtr5ewEe+7fw+qLVtO+t52yk8pY+OeFFEwuGPE4RGT0SmZhfQFwHHBv/GzgEne/ZlhRDpMK6yRqaICHHoJ//CNM8RdAQhhp5NRTQ7H96leH/tuFhemJVWQMSGFhrZzdg4bnGnjm/Gdo3tRM4cxCFvx2AeUnlqclFhEZfZJZWD8BzANqgdgbDZjq7mntK6DCOoV27oR77w1dR+65BzZs6L68sBBOPLGr0D7llFB8i8iApLCwVs7uRfP2ZlZftJr6x+qxfGPu1+Yy430zNGKIiPQrmYX1ZcBf3b0+Yf757n7nMIM8Ang/8ALhD8Eyd3/KzIqALwP1wHjgWnevS3x/upN0VtmyJbRox6bVq7svz8mBo4+Gk07qmg4/PMwXkQOksLBWzu5DR0sHL37kRbbeGK4xqfq3KuZ/Zz4Fk9Q1RER611fOzhvktl5OTNCRqYMP6wA/BS50961mNhP4G3AE8G3gLnf/tZmdC9wKXJ6E/clQHXQQXH55mABqauCf/wxF9oMPhhvUPPVUmG69NaxTXh5ateOLbY2jLZJqytl9yCnIYd635lHxqgrWvnMt1b+vZu9De5l/y3wmXTwp3eGJyCg02Bbrn7r7WxPmvQn4vrsP67t/M2sAjnf3NWY2CXgaWAxsAMa7e6OZ5QF7gQXuvjH+/ZnQ+iGRxsZQXD/2WJgefTS0cieaPTsU28cfH27TvmgRVFWNeLgi6ZbCFmvl7AFq2tDE2qvWUnt/LQCTr5jM3K/OpXCqrh8Rke6S2WI9wcw+6e6fN7Ni4FvAm4Anhhsk8DvCeKvnAW8BPgCcAVS7eyOAu7eZ2QZgCbAxCfuUVCguDv2tTz+9a962bd0L7eXLYePGMP3mN13rzZwZCuzjjuuapk0D9XsUGQrl7AEqOqSIY+45hq03b+Wla19i5y92svuO3cz+7GxmvG8GOfnqyiYi/RtUYe3urzezk83sRuB1wHbgGEKLxHC9F7iNkPC/6u6/N7OPArsT1qsHpidhfzKSpk+Hiy8OE0BbGzz3XCiwn3wSVq4MXUc2bw7Tbbd1vXfy5FBgL1oECxeG6bDDID8/PcciMkooZw+O5RgHfeAgJr5+Ii9c8wK779jNix96kVe+9wqHfuNQJpw9Id0hikiG67OwNrNX97JoE7AeuAGYBVwEfGiYsRQBvwRmAN+IWjkcaEpYrxBoTZjHrl27WLy4q1V+6dKlLF26dJghScrk5YULHI8+Gq66Ksxrb4d160KRHSu2n3wyjEpy111hisnPDxdExgrt2HTwwWrdllFh2bJlLFu2LPYyKX2glLOTo2hOEQtvX0j1HdW8cM0LND7byDPnPEPlWZXM+eIcyk/QqEci2WagObvPPtZm9ijhYpRaoKO31YBp7j6sjmhmdhfw7+5ea2afBP4f8EHg0+4+P2699dG8X8S/P9P660mSuIfuIitXhmnVqjC99FLP61dUwFFHhWnhQliwIBTgU6ao4JaMlaw+1srZyde+v50tX9/C5us307433Lyy6pIqDvncIZQsKElzdCKSDkMebs/MzgeWu/vOfnZwqbv/bhgBVgGr3H1a9NoIQzi9DbgLmOTuTWaWD+wD5rj71vhtjJYkLUmybx88+2woslev7iq4d+3qef3KSjjiiDAdfnjX89mzdRdJSbskFtbK2SnSWtPK5i9vZuu3ttLRFP5nmXjhRGZ+bCYVJ1ekOToRGUlJG8c6VcwsB3gROC0ausmAx4CzgWXAz9z9jugimbe6+xWJ2xhtSVpSZMeOriJ71Sp4/vkw7e2lS2lhIcyff2DBPX8+FBWNbOyStVI1KkiqZHPObt7WzKb/28T272+nY38osCuWVDDz2plMOHeCbjAjkgUyvrAGMLNjCBfDrACmAA+6+wNmNp5ws4GNhLFX/8fd9yW+f7QmaRkB7qHgjhXZa9Z0Pd+6tef3mIXxuufPh3nzuqb58+GQQ6BAN5CQ5BlthTUoZ7fsaGHLt7aw9eatnV1EiuYXMeO9M5jytinkV+riapGxalQU1sM12pO0pEl9fVehHV9wv/hiGLmkJzk5oQtJrNCOL7xnzQoXZooMwmgsrIdrrOTstro2tt2yjS3f2kLL1hYAcopzmPKWKcx47wxKjylNc4QikmwqrEUGq7UVNm0Ko5SsX981rVsX5vf2e5OfD3PmdBXac+bA3LlhmjUrdD0RSaDCevTraOtg9593s/XmrdTeW9s5v2xxGVPeNoUpl08hf6JasUXGgpQW1ma2CPgwsBL4feLdtUbKWEvSksGam8OoJLFCO77w7ukOkzFmYTjA+GI7/vn48SN3DJJRRrKwVs5OvYbnG9j27W1s/8l22utCNxHLNyZeMJGpb5vKhHMn6IYzIqNYqgvrW4EL3X2qmV0B/NHdE8cxTbmxnKRlFGlshBdeCEX2Cy+ELiUvvRQeN2+Gjt5GQCOMXJJYcMceDzpIo5eMYSNcWCtnj5D2pnaqb6tm+4+2s+fuPZ0DIOZX5VN1SRWTL5tMxZIKcvJUZIuMJqkurPOBydGV4TlAqbvXDWujQ5ANSVpGuVj3kvhiO/55Q0Pv783PD/26Y4X3IYeEafbsMI0fr3G6R7ERLqyVs9OgeVszO362g+0/2k7j842d8/MnxRXZr1aRLTIaJLWwjoZPmuDuPzezE4F97v5cEuIclmxL0jLGuIc7TPZUcL/4Imzf3vf7y8u7iuz4gjv2vELj7GayVBbWytmZxd1pWN3Art/sYudvdtK0ruvLgvyqfCacP4GqC6oY/9rx5JXpQmiRTJS0wtrM/he4BrjL3S+N5n0JeCbxrlojLVuTtGSJhgbYsKGr2N6wIdyRcsOGMPXV2g2hm0lPBXfsealGLkinVBXWytmZzd1pWNXAzt/sZNdvdtG0vqvItgKj8oxKJr5xIlVvrGLcrHFpjFRE4iWzsL4PuBh4m7t/M5q3BPipu89MRrBDpSQtWcsdamq6iu1YwR3/vKmfLrQTJ/ZccMde62Y5KZXCwlo5e5Rwdxqfb2T37bupvr2aun/VdbspfcnCEsa/djwTXjuBitMryC3SNRci6dJXzh7s90x3u3utmcVX45cOPTQRGTazUBhPnAiLe/g9dw+3eu+p4I493707TCtW9LyPqVO7Cu3E6eCDNXZ35lLOHiXMjJIjSyg5soSZ186kpbqFmjtr2H37bmr+VkPDqgYaVjWw5atbsEKj4rQKJrx2AuPPGU/pMaVYjq6xEMkEg22x/ndgDjADeAh4E6E15APufnNKIhwgtX6IDFFHR7gzZWLBHetmsmlT7zfLgTBaycEH9154T52qCyv7kcIWa+XsMaCjpYO9D+9lz917qPl7DftW7oO4P935k/IZf9Z4KpZUULmkkuLDi3VrdZEUSvbFiycDVwIzgVeAX7r7P4Yb5HApSYukSHt7uPV7rNBOnLZt6/2GOQDjxoWb4yQW3HPmwKGHhgsvs1yKL15Uzh5jWqpbqL2nlpq/17Dn7j00v9zcbXn+5HwqX11JxatDoV1yVIlatEWSKNXD7eUCZ7j7PcPa0DApSYukSXNzaNXurfDevbvv90+eHArs2N0qY8+zqOge4eH2lLPHEHencW0jtffWUvtAmFp3tHZbJ298HhWnhyK74lUVlC4qJadAw/qJDNWQC2sz+zPQ33ABVUC1u5859BCHT0laJEPV13fvWhKbYkMJ7t/f+3uzpOhOVmGtnC3uTtP6ps4ie+8De2ne0r1F2wqNsuPLKD+lnIpTKyg/pZzCaYVpilhk9BnOxYvPROusJVyffDXhNrjx2fAQYEIS4hSRsaisDBYuDFOijo7QlSR2p8rYreFfeCFMO3eG6ZFHDnzv1Klw5JFwxBHhMTZNmpTNfbqVs7OcmVE8v5ji+cVMv3o67s7+jftDkf3gXur+VUfjmkbqHqmj7pE6tnx1CwCFswqpOCUU2eWnllN6TKluuy4yBP0V1l8htGrXApjZXHf/VOJKZnZbKoITkTEuJyfcrv2gg+A1r+m+rL+ie/v2MN17b/f3TZjQvdA+5hg49tgwlvfYp5wt3ZgZRYcUUXRIEdOunAZAa00rdY/VUfevaHqsjuZNzezctJOdv9oJQM64HEqPLaVscRllJ5RRtriM4sOKsdys/adVZED6LKzdfW/CrClmluvu7bEZZnYGcGIKYhORbNZf0b15Mzz/PDz3XPeppgYefjhM8ebMgeOOC9OiRXDCCWGIwjFEOVsGIn9CPhPPm8jE88LPv7c7Dc81dBbaex/ZS9O6JuoeraPu0a673eeW5lJ6XFRsRwV30dwijUAiEmewg8/+DnjazB4G6oEjgNcBB7SIiIikTE5O141szjuva747vPJKV5G9ejU89RQ880y4a+VLL8Hvfte1/lFHwZIlXdPkySN9JKmmnC39slyjdGEppQtLmb50OgCte1rZ9+Q+6p6oo355PfVP1NO8uZm9D+5l74Nd/7/lVeZRenwp5SeUU3p8KaXHllI0p0ijkEjWGspwe9OB9wJHA3XAb9097V8r6kIYEelVayusWQNPPgkrV4Yb4Sxf3v3CSTM47TS47DK44orQpWSEpHi4PeVsSYqWnS2hyI5NT9TTsr3lgPVyy3IpPaaU0kXRdGwpJQtKNBKJjBmpHm4vB1jk7r3csm1kKEmLyKA0N8Pjj8MDD4TpoYfCPIDiYrjwQjj1VDj5ZDj6aCgoSFkoIzzcnnK2JIW707KtJbRqP1HPvpX72LdyX4/FtuUbxUcWU7aorLPYLj2mlLwK3bVVRp+kFdZm9oMeZk8C2t39oiHGlxRK0iIyLHV1cMcd8KMfwd13d182bhwcfzyccgqcfTacfnoovpMkhXdeVM6WEde8vZl9T+0L08rw2LS+qdvdImPGzRnXrdAuOaqEcbPGqSuJZLRkFtYPAffQ/dfjVcBD7v75YUU5TErSIpI0a9bA/ffDo4+Gae3a7ssLC0O3kde/Hi6/PAz9NwwpLKyVsyUjtNW30bCqobPQrl9ZT8OqBrzlwBoktzSXkqNKKFkYNx1VQkFV6r41EhmMZBbWh7n72oR5hwNXuvvHhhdm5/ZOBk4DXgAeBhqALxMuvBkPXOvudYnvU5IWkZSpqenqNnL33aGvdix35uaGCyg//WlYPLTaOIWFtXK2ZKyO1g4a1zR2FtsNqxrYt2rfAXeOjCmYWnBAsV1yZAm5xbkjHLlku1T3sT4EeMLdq4a1obCttwNz3f26uHk/BO5y91+b2bmTiFR/AAAgAElEQVTA29z98sT3KkmLyIiprg4F9q9+BXfeCW1tYf4XvgCf+MSgNzfCfayVsyWjtexqoWFVQ5hWh2K7YXUDHQ0dB65sUHRoUVexvSBMRfOKdIMbSZlktlhvoPtXirnAFOCn7n71MINcAnwSOMejoKKr2TcA49290czygL3AAnffGP9+JWkRSYtdu+ArXwkThNu1z549qE2ksMVaOVvGBO8Id5CMtWo3rA6Fd+PaRmg/cH3LM4rmF1GyoITiI4tDwX1kVHBrdBIZpuHc0jzR3cAv6ErUHcAOd183jPhivgY8AtxoZocCnwdmAdXu3gjg7m3RH4olwMYk7FNEZOgaGkJf7CefDK/z8sKFjplDOVvGBMsxiuYUUTSniKoLu75s6WiOupOsCl1JGp9rpOHZBvZv2E/jc400PtfYfTt5RtG87gV38ZHhFvAquCUZBltY/8zdH4yfYWaHmdmh7v7CUIMws8OARcDb3f0ZM/socBfwBWB3wur1wPSh7ktEZEja28Ot1FetCv2tH3oojIfdGvUHraiAb3972BcyJplytoxpOYU5YczsY0q7zW9vaKdxTSMNzzXQ8GxCwf18I43Pdy+4yYXiecUULyim5MiS7gV3oQpuGbjBFtbnAg8mzNsB3EdIskO1AKhx92ei1zcBnya0sjQlrFsIHHBlw65du1gcd+HQ0qVLWbp06TBCEpGs094OW7d23aXxxRfD47p14U6O8TeUgXBTmZNPDjeVufJKGD9+wLtatmwZy5Yti70cdn/nXihnS1bKLcml7Pgyyo4v6za/vbGr4G58tqvw3v/SfhrXNNK4ppHq31fHbSjqwx11JYkV3sWHqeDONgPN2QPqY21m7ye0RJT1ssrd7n7uYIOM2/5lwDJ3r4ybtxL4DaFFZH7c/PXAp939F/HbUH89ERmwNWvgT38Krc51daFLR20tbNnS1QLdk4MPhoULYdGiMNzeKaeEluphSnYfa+VskcFpbwoFd6xlO1Z4N73Y8/jb5ELR3IQ+3AtKKJpfRO44jVIy1g27j7W732RmDwJvA76VsLjF3V8ZZozPABVmNsXdd0Tz2oBtwHQzK3L3JjPLB2YCDwxzfyKSLdxDF46HHoIHHwzThg29rz91KsyZA3PnhsfY8wULoLKy9/dlEOVskcHJLcqlbFEZZYsSWrib2mlc21Vwxx6bXmyiaV2Y+GPcG3KiFu5Y6/aCMCygupRkjwF3BYn60X3M3Q9ozjGzme6+eahBuPsaM/srcClws5lVAgXAz4HzgLOAO4Czgd+7+9ah7ktExrDW1tBtY+XKrumpp2DPnu7rTZgAb3gDnH8+TJsGJSVQXg4zZiT1jorppJwtMny5RbmUHVtG2bEJBff+dprWNh3Qh7vphbiC+0/xG4rrw72gpFsLt4YFHFv6LKzN7CPAne7+XDTrv80s8TsOA14DnDHMWP4T+KaZFRBaOC6Prih/D/BlMzsGmAqoE55Ittu7N9wNcc2arun550PLdGxM6XhTpoSuG6efHqajjw4jeIwxytkiIyN3XG7PF03ub6dpXVPoTvJsGIc71qWkpz7clmcUHVbUrdguXlBM0aFF5OSp4B6N+vvLcggQ/93ncUAJsJ2uXkdGSJ7D4u7VwH/0MH8PMKzxVkVklGlrC/2dYxcRbtjQ/fmuXb2/95BDQh/o+GnatHCh4dinnC2SRrnjcik9upTSoxMK7qgPd8Oz0UWTz8ZdNPlsI43PNrKLrrxmBUbxYcWdhXbJUVEL95wiLDcrctmo1Wdh7e7vS5j1v8Bqd2+Jn2lmRyY7MBEZo9zDLcJffrlr2rKl++vNm3tueY4ZNw4OOwwOP7z7NH/+mOnKMRTK2SKZqdc+3A3tNDzfvdhueLaB5k3NnXefjJczLofiww/sUjLukHFYjgruTDDY70IXuvuT8TOi8UznAs/1/BYRyRrNzfDKKwdO8cXzli3QlDgiWw9mzAitz7ELCOOfT50KOfqadACUs0UyWG5JLuWLyylfXN5tflt9G43PNx7QpaR5SzP7ntrHvqf2dVs/pyin2+gksZbucTNVcI+0ARXWZnY64evDc6K7aMWbQLjq/PYkxyYimaK+vueCOXFKvEiwN+XlYei62HTQQd1fz5qVaXcwHFWUs0VGt7yyPMpPLKf8xISCe29b1wWTca3cLdta2LdiH/tWJBTcJTmdN7wpOaprpJLCgwqx7OgeN+IG2mJdC/yE0H/vVQnLWoDvJTMoEUmx1laoroadOwc2NTb2v02A3NzQmjxtWtfjtGkHFtDl5f1vS4ZDOVtkDMqryKPilAoqTuk+fn7rntYDiu2GZxto3dFK/RP11D9R32393PLc7kMCRlPB9AIV3MM00HGsV5nZGcCb3F0JWSTTtLTA7t1hqq4O065dvRfKNTWD2/64cV1Fcl9TVZW6aGQA5WyR7JI/Pp/K0yqpPK37WPutu1u7FdqxLiWt1a3UPVpH3aN13dbPG58XRjs5tmsqPqKYnALl9YEazDjWe+mllcPMTnP3h5MWlUg2a23tKpDjC+W+XtfV9b/deDk5oQiePHlgU2lptoyqMWYoZ4tI/sR8Kl9dSeWruxfcLTtbOovt+Fbutpo2au+vpfb+2s51Ld8oPrK4W7Fdekwp+ePzR/pwRoX+xrH+M1Da1zqEfnwLgMnJCkpkzGhtDa3D/RXG8a/37h38fnJzYeLEUCxXVXU9nzKl50J5woTwHhlTlLNFZCAKJhdQMLmA8a8Z3znP3UNf7af3dV4gue+pfTS90ETD0w00PN3Ajh/v6Fy/cFZhZ6FddlwZZSeUUTitMB2Hk1H6a7FeS+irt5WuMVATGfCGZAYlkpGam7t3txjI41CL5AkTuork+EK5t9fl5eqCIaCcLSJDZGYUziikcEYhE8+f2Dm/bV8bDasauhXbDc+EIQGbNzWz+7bdnesWzCig/IRyyk4IhXbZ4rKsa9nur7C+Hqh39+a+VjIzXV0uo0tT04F9kvsrlPft63+7iXJyQhGcWAj3VSxXVKhIlqFSzhaRpMorPfCCyY62DprWN4XW7ZX7qF9eT/3yelq2tlC9tZrqP3XdXbJoXlFnoV1+Ujllx5eN6T7b/d0gprqv5XEuRleZSzq4hxErBtOKXF098FEu4uXldRW/A31UkSwjSDlbREZCTl4OJUeUUHJECVP+fQoA3uE0rmvsHIWk/ol66lfW07S+iab1Tez8xc7w3nE5lJ1YRsVpFVScHgr2vIrB3lYlcw3rSMxsJvA+4P0oSctwuYdW4cEUyLt3w/79g99XQcGBRXB/hXJ5uS7gk1FNOVtEUsVyjJLDSyg5vISpb50KQEdrBw2rG6h/op66x+uo+1cdjc81svfBvex9MOoqaVBydAkVp1WECy3PrKSgqiCNRzI8gy6szSwXuAC4GjgH2EHo0yfSxT2MVDHYIrmlpf9tJxo3bnCtyBMnapQLyRrK2SKSLjn5OZ23cp++dDoQhgDc+8he9j60l70P76V+eX3nxZHbbt4GBqXHlTLhtRMYf854Kk6tIKdw9HzzO+DC2swOAd4JvB2YCGwC3gj8DTglJdFJZtm/P4yBvGNH13jI8c/jX1dXQ1vb4PdRXNxzMdxXoVxcnPxjFRnllLNFJBPlT8yn6o1VVL2xCoD2pnbqH6+n9qFaau+rZe/DezvvIrn5i5vJKc6hckklVRdWMfHCiRROzeyRR/otrM3sTYSWjjOBF4CvAT8GPu7ud0WrPZKyCCV13MMtqHfs6L1Ajn9eX9//NuOVlQ2uFXniRCgqSs2ximQJ5WwRGU1yi3KpXFJJ5ZJKuA7aG9upfbCWPXfvYc/de2hY1UDNX2uo+WsNvAfKTymn6uIqJl0yiaI5mVczDKTFuolwC9xHgHe6+zoAM+ttKCdJN/cwzNsrr8C2bQdO8fOb+xw8oLv8/DAGcuLYyPGvY8+rqqAws/+rFBmjlLNFZNTKLc5l4rkTmXhuGPKveXszNXfWUP3HamrurqHukTrqHqnjpY++RMVpFUy9ciqT3jSJvPLMuACy3yjc/Q7gjuiil6vNbDZwG9A5MKGZmbsraY8E99DNYtMm2Lw5PL788oHFc1PTwLZXVhZuRd1XkRx7XlGhfskiGU45W0TGksKphUy7ahrTrppGW30bNXeFIrv6z9XsfTj0017/gfVMetMkDvrQQZQdW5bWeG2wudXM8oCLgP8AVgPfB17r7suSH97ALV682JcvX57OEJKjowO2bIGNG0PRHF9Ax54PpGguKYHp03uepk3reizt7yZtIpJqZrbC3RenaNvK2SIy5rTta2PX73ax/Ufb2ftA183YKs+sZOa1Mxl/zngsRY2BfeXsQRfWCRs+FPgAcLW7p/UKslGVpDs6YOtWWL+++/TCC/Dii/13z6iogFmzwjRzZphmzOhePJel9z82ERm4VBbWCftRzhaRMafpxSa23rSVbd/dRkdDBxAK7LlfmUvZccmvh/rK2cPqkOLuLwDXmNmK4WxnzHIPrc+rVnVNq1fD2rV9j708ZQrMmRMK5vgCOvZYUdH7e0VEeqGcLSJjUdHcIg79+qHM+vQstt2yjZevf5nae2tZcfwKpr9vOnOvn0tuSe6IxJKUnt7u/pNkbGdUcw+tzY89FqYnnwxF9N69Pa8/eTLMm9c1HXpo16Nam0UkhZSzRWQsyq/MZ9bHZjF96XQ2/e8mtn5rK9tu3saev+1hwe8WUHpM6ru/ZsYllKNRWxusWAH/+Ac8/DA8/jjU1By4XlUVLFwYpqOOCo9HHhnu4iciIiIiSZU/IZ9Dv3ooU94yhTVvW0PDqgZWnraSI397ZOdoI6mSUYW1mRUATwDXuPv9ZlYEfBmoB8YD17p7XdoC3LkT/vhH+Nvf4N57D2yNnjwZTjopTCecAEcfHbp1aCQNERmjMj5vi0jWKltUxnGPH8fat69l5692svqi1Rzzt2PCmNkpklGFNfBRYHbc628Dd7n7r83sXOBW4PIRjaihAX71qzDde2+48DBm7lw45xw44ww4+eTQ/1lFtIhkl8zL2yIikdxxuRzx8yPIG5/Htu9sY/VFq1n8zGLGHTwuJfvLmMLazE4FXgH2RK+nA1cA74tW+QfwezOb7e4bUx7Qtm3w9a/D974HtbVhXl4evO51cNFFoaA+5JCUhyEikqkyLm+LiPTAcox5N85j/+b91PylhrXvXMvRdx2dkuH4cpK+xSEwsxLgUnf/QdzsM4Bqd28EcPc2YAOwJKXBNDbCpz8dLiS84YZQVJ98ciiwd+yAO++EpUtVVItIVsuovC0i0g/LNQ773mHkjc9jz9/3UHtfbUr2kxGFNXAt8KWEeTOA3Qnz6oHpPW1g165dLF68uHNatmwI9z54+mk4/nj43OdCgX3xxWGEj3/9C97xDpgwYfDbFBHpwbJlyzrzFVCV7niGYFh5Oyk5W0RkEAqnFnLwhw8GYPOXNw/qvQPN2WnvCmJmrwOWu/vOhEUOJN5isBBo7Wk7kyZNYlg3G7j3XrjggtCn+ogjYNkyOO20oW9PRKQPS5cuZenSpQCYWXWawxmUZOTtYedsEZEhmP7u6Wz8zEb2/GMPLbtaKJhUMKD3DTRnZ0KL9UeAH5hZdRTowcBtQAnhivJ4ZcC2pEewciW8/vWhqP6P/4Dly1VUi4j0Lv15W0RkCPIn5lN5ViW0w56/70n69jOhsL4CODZu2ga8E/ghMD0augkzywdmAg8kde9NTfCmN4U7Ib7tbfCTn0BxWu/0KyKS6dKbt0VEhqHy1WG4vbrHkj8SaNq7grj7rvjXZtYO7HL3zWZ2B3AWcAdwNvB7d9+a1ABuvDHcMfGoo+CWWyAnE/7XEBHJXGnP2yIiw1B2QrjD9b6n9iV922kvrPvxHuDLZnYMMBVYmtStt7fDN78Znt9wA4xLzZiGIiJZJLV5W0RkmIrmFAGwf9P+pG874wprd58d93wPcHXKdnb//WG86rlz4bWvTdluRETGshHN2yIiw1R4cCEYNG9ppqO1g5z85PVWyO5+D/ffHx4vukh3TBQRERHJAjkFOeRNyIMOaKttS+62k7q10ebRR8OjRgARERERyRp5FaHTRtteFdbJs359eFywIL1xiIiIiMiIySsPhXV7XXtSt5u9hXVLC7z8chgFZNasdEcjIiIiIiMktyIXgLY6tVgnx+7d0NEBkyZBwcDuuiMiIiIio19OQSiBvcWTu92kbm00qa0Nj+MTbxImIiIiImOZ5YdBK7xVhXVy7IluY1lZmd44RERERGRExQrrjtaOpG43ewvruug2lmVl6Y1DREREREaU5anFOrlaWsJjYWF64xARERGRERW7KYwK62RpbQ2P+fnpjUNERERERpT6WCdbrMVaI4KIiIiIZBXLjQrrdhXWyaEWaxEREZHsZNFjcutqFdYqrEVERESyTFRYu6vFOjnaojvt5OWlNw4RERERGVlqsU6yjmjcwpzsPQUiIiIi2cgs1mSd3O1mb1UZa/qPnVgRERERyQ5qsU4yFdYiIiIi2Ul9rFNEhbWIiIhIdlGLdZIl+T8UERERERkd1Mc62dQVRERERESSKGMKazO7xMw2mNluM/ummeVF86vM7FYzu97MvmZmyRl4WoW1iMiQjXjOFhFJprHcFcTMZgIXAZcCHwLeCXwwWvxb4Afufi2wDvhcUnaqwlpEZEjSkrNFRJIpRRcvZsrdUWYBV7l7G7DCzI4CXmNmDwOHuftj0Xq3A2vN7Avu3jCsPaqwFhEZqpHP2SIiyTSWW6zd/aEoQcdsBTYDZwIb4tbbGj1dnISdhkcV1iIig5KWnC0ikkTZdvHiCcAtwAxgd8KyemB64ht27drF4sWLO6dly5YNbE8qrEUkDZYtW9aZr4CqdMczTCOXs0VEkmGQLdYDzdmZ0hWkk5kdCuxy96fNzIGmhFUKgdbE902aNInly5cPfEcabk9E0mjp0qUsXboUADOrTnM4QzZiOVtEJJkG2cd6oDk7o1qszSwXeBfwsWjWVmB83HIDSoBtw96ZuoKIiAzLiOZsEZFkGst9rON8BLjB3Zuj13cDh8Utnw20AMNv5lBhLSIyXCOXs0VEkmjM97E2s48DK4EiM5tjZlcBtcAr0RXnAOcCN7p7y7B3qMJaRGTIRjxni4gkU4parDOij7WZXQd8PmH2Gnf/gZm9GfiUmW0ESoFPJGWnKqxFRIYkLTlbRCSZxvI41u7+BeALvSzbBLwjBTsNjyqsRUQGJS05W0RkFMiYriBpo8JaREREJLukqPxTYS0iIiIikgQqrEVEREQkO43VUUFEREREREaCpagrsAprEREREZEkUGEtIiIiIpIEKqxFREREJDupj7WIiIiIyDBouD0RERERkcyVvYV1aSlMmwZlZemORERERERGUG55LgXTCsgtyU3qdjPiluZp8b73hUlEREREssrs62Yz+7rZSd9u9rZYi4iIiIgkkQprEREREZEkUGEtIiIiIpIEKqxFRERERJIg6wvrZcuWpTuEjKFz0UXnoovORRedi/QbK5+BjiOz6Dgyy2g+DhXWo/jDSzadiy46F110LrroXKTfWPkMdByZRceRWUbzcWR9YS0iIiIikgzmnuSbpKeJme0CNg3hrVVAdZLDGa10LrroXHTRueiSqnMxy90npWC7GUs5W8eRYXQcmSXTj6PXnD1mCmsRERERkXRSVxARERERkSRQYS0iIiIikgQqrEVEREREkiBrC2szKzKzG83s/8zsO2ZWnu6YUs3MzjGzx8xsdty8Xs/DWDxHZnaJmW0ws91m9k0zy4vmV5nZrWZ2vZl9zczy497T67LRzsyWmNlaM9tjZt+Im59VPxcxZlZgZk+b2RnR66w8D+k20HNrZueZ2U1m9m0zu3Sk4+xLb7mmh/XmmVm7mXk0LRrpWAfCzL4QF+PTvayTsZ8HQPT3zxOm83tYr9LM6uLWuTgd8SbENKi/3z28PyM+m16OY0C/K9G6mf/74u5ZOQE/BN4cPT8X+GW6Y0rx8U4BLgYcmD2Q8zDWzhEwE/gJcDzwn0AD8JFo2X3ASdHzdwNfjHtfr8tG8wRUAB8HxgOvB1qBs7Pt5yLhnHwC2Auckc3nId3TQM4tMB94EsiJXj8IHJXu2KNYes01Pax7PfBa4GxgSbpj7yXGEuCWKMazgXmj6fOI4jkC+AZwOnByNK0GxvWw7kejv5ex481Nc+yD/vudiZ9NT8cxmN+VaP3M/31JdwBpOWiYDjQDxdHrvOjDnJ3u2FJ83DkJP9C9noexeI6ihJoX9/rLwF+iBLstbv4MYF/0x6TXZek+niScjyKikYGi1yuA12Tbz0Xc8Z8KXAVsBM7I1vOQ7mmg55ZQ6P1f3OtPAD9Md/xRLD3mmh7WmwLcARyU7pj7OZ4PRtMBReho+DyieGYmvJ4F/LqH9QqBu+nhn4c0xz/gv9+Z/Nn0cBwD+l2Jlo2K35ds7QpyBlDt7o0A7t4GbACWpDOoVHP3joRZZ9D7eehr2ajk7g9FxxGzFdgMnEk4tth6W6Oni/tZNqq5e5NH2crMSoDngfvJsp8L6Dz+S939B3GzzyDLzkOGOIOBndtuv5vA2ui9addHrkl0OeGf2ZfN7IfRz2FGMTMD/h34KrDDzC7vZdWM/TwA3D3x/F8M/KmHVd8AnASsM7O/mFlVyoMbgEH+/U6UMZ9N4nEM4ncFRsHvC2RvH+sZwO6EefWE/wCzSV/nIRvO0QmE/+Sz+jyY2enAX4F8oJjsPB/XAl9KmJeN5yETDPTcJq6Xyec/lmu6cfdvAOXABYSvtr8/wnH1y4OTCTfsuBH4uZmd18Oqo+nzADif8I1lN+7+e0I3uSWEb6f+bGaZWCsNJgeNps+mx98VGB2/L5C9hbUDTQnzCgl9TLNJX+dhTJ8jMzsU2OXuT5PF5yHyAvAj4CzgK2TZ+TCz1wHL3X1nwqKsOg8ZZKDnNnG9jDz/CbnmAO7e7u63A68DLjGzaSMa4AC5+x53vw74IqFbyAGrMAo+DwgXpAP73b2up+XRPxMPElp6Dye0YGeaweSgUfHZ9Pe7AqPj9yVbC+uthAu24pUB29IQSzr1dR7G7Dkys1zgXcDHolndjjX66rOEHs5DwrIxwd1fibpAfITQSpNtPxcfAX5gZtVmVg0cDNxG+Jyz6TxkioGe28T1Mu7895BreuXuzwH3EC7mymTfoucYM/7ziHMB8Of+VnL3HcCvyczPZDA5KOM/m8H8rkBm/75ka2F9HzDdzIoAouHTZgIPpDWqkdfXeRjL5+gjwA3u3hy9vhs4LG75bKAFWN7PsrFmOSEBZ9vPxRXAsXHTNuCdhCvus+k8ZIqBntvE38350XszSWKu6U8DsCaF8SRDB2GEiUSj4fOIGVBhHWkHnkphLEM1mBw0Gj6bwf6uQIb+vmRlYe3u2wlXlp4VzTob+H3chWljUtTaCmDQ93kYq+fIzD4OrASKzGyOmV0F1AKvmNlR0WrnAje6e4u7L+9t2YgHn2RmVpwwBuj5wDez7efC3Xe5+5bYRPhDuiu62ClrzkOm6OvcmtlHzCxWINxC+Do4ltvOJLSmZoSeck00Bu/nzGxqtM6VZjY5en468IC7701j2Acws+lmdoWZ5UTn+cOEUSUYTZ9HTHTB27j4rl/xx2Fml5rZnOj54UCNu69NT7RdBvP3O1o/Iz+bxOOI5vX0u3KoBaPq9wWiobaykZmNJwzrshGYCvyPu+9La1ApZGZlwFuAbwOfA25y9119nYexdo7M7Drg8wmz17j7EWY2C/gU4VhLgU/ErlTua9loZmYnEi7eWQM8Cjzu7r+NlmXNz0UiM9sIXOnu92fzeUin3s6tma0gjCP/u2i9twDHRW97wN1vS0e8iXrLNYRYnwMud/dHzew24FWEFsW/uvtPRjbS/pnZAsLFzXXAw8A33H1NtGxUfB7xLNwc5aDoQrjYvM7jMLObCH8r7yWMknRTDyNyjKgh/v3OuM+mp+MgdP/o7e9yEaPs9wWyuLAWEREREUmmrOwKIiIiIiKSbCqsRURERESSQIW1iIiIiEgSqLAWEREREUkCFdYiIiIiIkmgwlpEREREJAlUWIuIiIiIJIEKaxHJKGY2ycz+bGabzOxj6Y5HRET6przdRTeIEZGMYmavA+4Dygl3q5vv7jXpjUpERHqjvN1FhbWIZCwzexA4cyzcQl5EJBtke95WVxARyUhmdjDw92xNziIio43ytgprSQEz+28zuzndccQkIx4zO9vM/mVmZyQprDHFzG5I8vYMuAz4UjK3KyLdKV8Pn5nNNrPvmtmtSdperpldbWabBrDu683ssaEeazJzt/J2oMJ6lDGzI8zsD2b2PTN72szczN6S7rgS3An8Mt1BxElGPC8DxyQhloxkZuebWbWZvWJmp5jZwdHPmZvZTWZWGa13opmtNrOvmllRNG8KsD3JIV0BfN/d26Lti4w6ytdDMhrz9c7osTBJ28sBaoGZA1h3OXBEbwvNbIGZ/Tr62bvOzKbGLUt27lbeRoX1aPQn4LPu/k5gEfCT2AIzO9rMlqQjqPh9u/tqd384HXH0ZDDx9HYO3X0tUJ3MuJL5eQ13W+5+J/BtoM3d/+XuLwNXAy3AY+5eG633OCGRf9Tdm6K3vxG4Y1gHEMfMPgl8BviHma0BDkvWtkVGmPL1II3GfO3ujcC2ZO3L3VuBFQNcdwehCO9t+bPAsujlDe4eX0gnLXcrb3dRYT2KmNkEYD7QAODuHcCnADezCkLStjTElbZ9J9MAjqNjBPeVjm39GjjIzI4HcPfdwD3AG+L2VQlsi372Yg539zXD3Hcnd/+8u89z98Xufri7P5isbYuMFOXr1MrAfJ3skSAGE39/67ZHj4n9npOWu5W3u6iwHkWioWvWAneZ2cnRvE3A44R+TYcAS83s7QBm9jkz+2T0NVBnPyozW2xmXzOzD5pZo5mtNbP/MrOLzOwJM8T1GwoAAAd5SURBVLvYzP5oZnVm9q649/W4vYR9X21mP4rvaxZ9HXqjmX3RzO4zsxOj+X3uL+79ZmbvMrNnzewcM3vGzHaY2ZUD2Me8+Hj62ecB57AH46OvduvN7OqEc/qZ6Nz8IUrEAz1nB+zLzC41s8fN7EozW29mX+pjez199j3G05eoZeNZ4E1xsycBr7eo2wdwAfDnuDiLgMa414eb2ZfM7AYze8BCl5I3m9lyM7vQzH5pZrVmdo2ZHR8d424zO7+/+ERGE+Xr7MnXPentGPvaj5kdYmY/sND6+5WB7KeH/d5pZj83s1f1s55yd6q4u6ZRNAELgE2E/1B/AEyOW7YROCN6fhTwUvR8EuG/6YmE/7g3A4ujZTcDj0TPC4Aa4CbCWJSXA1v72l7ivgl9zG4FfhTNH0co1iqj15cCu4CyvvbXw3FPjfb58WjdLwD7gVn97CMxnj73GX8Oe4hhI/BVoDh63+Zofjnwq7j1/klIigM6Z73sawLh671vACcCZ/e1vYTPvsd4BvjzdR3wYvR8LvBFoBm4JO7nxeLWvwA4PXpeDDwBFEevH4zOQzmwF/gcUAS8OXp9GeGf+w8T/Qxq0jSWJpSvsyJfR8s/M8Dz2Ndn/QRwRNx7fIA/Z7HP8xTg0z0sPyPaT17cPOXuFE1qsR5lPLQqLiQkmv8EnjWzni7SWA+81czygdOjeaWEX+SDgVj/2OcJvzC4ewtQD/zO3esIv2hT+9leYnzNdL8Y4vVAk0d9dIE/EH4hL+lnf4nbjW3zF9G6XyD0/z2nn310i2cw++zF7R760z0BTIs7xolRi9IHCecUBnjOejneGkICu83dH3f3fwxie73FMxC/BuZY6A7yb8B3gb8Bl5lZKVDvUeaNvAp4JHp+AeGPV6wV5GLg89F5rgXu9dAvezlQ7u6/8fD1+EoG9xmIjArK19mRr3vQ6zH2sZ+zgQp3j8WzfJD7PAH4D3f/7ADXV+5Okbx0ByADZ2a5wMHuvhH4f2b2fcIV1MuAk+LXdfdmM5sBXAv8OLYJd99pZsuBswj/UR9CKKY63xr3vI2ou1Bv2+sl1PhtzAPy4+LqMLONwEF97a8/7r7fzF4CKoGqQexjyPtM0EbX789MYIO7fyNxpYGcMzN7K6GAjTnMw1fGHh/rID6DXuPpj7uvN7OVhO4gFe7+kpn9mtCK9CbgL3Fxx342Yv33ZhFamGLb2h2/6bjn7XTXgbqlyRijfN3t+LIhX8fr9Tz28dkcTfgnYqjOBOaZWam77+tlnY7oGJS7U0gnZHTJJ67/q7s/Tfiq7fDEFc3sTOC97v4FDyM8xHsfcIKFYZ+eAa7vb8f9bK8vmwktoPHDEBmh7+FwFQDrUryPgXiF7v2QMbOTBnHO/gwcGzf1eHX5ILbXYzwDPxx+BbwdeCkuvlzgGsLXpjEnEfqLxmwFTjOzkrj99tnPT2QMU77uLpvyda/H2Md+6qL3FDA0NxBa+7/Zy/K93nXRuXJ3CqmwHn0+YGbT417PAO6KnrcQLtY4nDC0U4WZFZrZ2dHyyWZWRejr933C1ziPEH21GMkh4b90M7N+tpe4b4vbxp+AfcC/R9uqjPYXa/nsbX+9mRatMyN6710D2Ed8PP3tM/44EvXUemGEVqgS4A4ze52Fi2T6+ww69wVMc/c1cVNr3C5y4573tb34uHuMx8wKzGxZ1M2jL78hfAX9OwB3rwf+CvzTu48G8lpCN5GY2Pn+pYWxsN8PxMYyTfwMEj/nUT1CgUgvlK/Jmnw90PPY434IubSA8M8XhG8nMLPpA8jdRmg9fivwJjO7LHrvh83sOkJr/T1x6yt3p1IqOm5rSs1EuCDCgT3Az4AfEb5KGh8t/zSwBbiQ8NXOekILwesJLR0/IiSpPxMupGiJttdI+EU7D2glXCBzMPA/0fKr+tpewr6vAR4i9FtbEC1bTGjp/Cqhr+FR0fxe99fL8Xv0/msIXRMWxi3rbR+z4+Ppb5/x5zBh3xcSLuK7KToXsfe9PVq+BFgdfTZfjc7zQM/ZhT0c66XR/n4JzIjm9fWZdttWL/EcRLiRwS8G8LP2rYTXbwZekzDvhh7e92rgOWA3cF0073zCz9qNwHSiIceAdxP+8H43OtY3pPt3TJOmZE0oX2dTvp4JPExodT+2n2Ps67N+PeGbwkcJF5KvA94bHXuPuZvwrUgroaV6AnAvoah/J/AuwnjeXyV0Q4m9R7k7hZNFJ0+yhJlNBt7t7p+LXucQWlHe5e7XpTW4fpiZw/9v545tIgaCKID+ESkRHdAKbZBQABIRvdAEEvkVQAlkJNDHEPgCJwg4jbiz7r3MloONvr68s5vrXmYWOVBVPXb3QVc5Af9HXrMmu7fBKMj5uc8yFnC1eneT1aG0U7TaerLtdKCquqiq2yTPx14L8CvyGtm9MYr1+XnKMov1VlUfSXZJPrv79bjL+t7+MMfD/vGulmvf+LvLJC/d/f7jl8ApkNcksntTjIIAAMAAf6wBAGCAYg0AAAMUawAAGKBYAwDAAMUaAAAGKNYAADBAsQYAgAGKNQAADFCsAQBgwBckvFs/4nJn4wAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(12,8))\n", "plt.rc('font',family='Times New Roman')\n", "params = {'mathtext.default': 'regular' } \n", "plt.rcParams.update(params)\n", "\n", "plt.subplot(2, 2, 1)\n", "plt.plot(vehicle.v_kmsc, vehicle.h_kmc, 'k-', linewidth=2.0)\n", "plt.xlabel('Speed, km/s',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 2)\n", "plt.plot(vehicle.acc_net_g, vehicle.h_kmc, 'b-', linewidth=2.0)\n", "plt.xlabel('Deceleration, Earth g',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 3)\n", "plt.plot(vehicle.q_stag_total, vehicle.h_kmc,'r-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-rate, '+r'$W/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "\n", "plt.subplot(2, 2, 4)\n", "plt.plot(vehicle.heatload/1.0E3, vehicle.h_kmc, 'm-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-load, '+r'$kJ/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.savefig('../plots/genesis-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/genesis-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/genesis-earth.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 }