{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 19 - PAET - 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 PAET flight test.\n", "https://ntrs.nasa.gov/search.jsp?R=19720045287" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \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('PAET', 62.1, 69, 0.0, 0.66, 0.0, 0.46, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(90.0,0.0,0.0,6.56,0.0,-40.8,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 (2400.0,0.1,0.0)\n", "# bank angle = 70 deg arbitrary, actual profile will give more realistic results" ] }, { "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": "iVBORw0KGgoAAAANSUhEUgAAAs8AAAHtCAYAAADm/pMXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXxU1f3/8dcnIYR9DYIgu+ybSNxwAVGoe10Q61Ztaen3669Wa9262Vq1X2utrVtr0UJtqbtVFDdwgYpLlUU2ZZNN2UwCAWQnOb8/zgyZJJNkJszMnUzez8fjPGZy5869n7niJ5+cOfccc84hIiIiIiI1ywo6ABERERGRukLFs4iIiIhIjFQ8i4iIiIjESMWziIiIiEiMVDyLiIiIiMRIxbOIiIiISIwaBB1APPLy8ly3bt1i2regoIB27dolN6B6Rtc08XRNEy9dr+ncuXMLnXPpF1gSxZOzw9L1v1860rWKna5VbHSdylSXs+tU8dytWzfmzJkT0775+fkx7yux0TVNPF3TxEvXa2pma4OOIdXiydlh6frfLx3pWsVO1yo2uk5lqsvZGrYhIiIiIhIjFc8iIiIiIjHK2OJ5woQJQYeQcXRNE0/XNPF0Tes2/feLna5V7HStYqPrFBtzzgUdQ8zy8/OdxuKISF1kZnOdc/lBx5FKytkiUldVl7MztudZRERERCTRVDyLiMSotLSU/fv3Bx1GRisthT17go5CRKRqKp5FRGKwceNGsrOz6dKlS9ChZKxZsyA7G848M+hIRESqpuJZRCQGDRr4afEPHDgQcCSZq0kT//j118HGISJSHRXPIiIxUPGcfM2a+ccdO4KNQ0SkOiqeRURioOI5+cLFs3qeRSSdqXgWEYlBuHjWDYPJo+JZROoCFc8iIjFQz3PyRRbPdWgJAhGpZ1Q8i4jEIFw8l5SUUJcWl6pLcnIgNxdKSmDv3qCjERGJTsWziEgMzIzs7GzAF9CSHBq6ISLpTsWziEiMNHQj+VQ8i0i6U/EsIhIj3TSYfOHiefv2YOMQEamKimcRkRip5zn52rTxj1u2BBuHiEhVVDyLiMQoJycHUPGcTHl5/rGwMNg4RESqouJZRCRG6nlOvnbt/KOKZxFJVyqeRURipOI5+cI9zwUFwcYhIlIVFc8iIjHSDYPJp2EbIpLuVDyLiMRIPc/Jp+JZRNKdimcRkRjphsHkC4953rQp2DhERKqi4llEJEbqeU6+Ll384xdfBBuHiEhVVDyLiMRIxXPyhYvndeugtDTYWEREommQypOZWT/gh8BKoBcw0Tn3iZk1Bu4BdgCtgVucc1pfSkTSSnZ2NlB/iucgcnazZn6hlC1b4KuvoEOHRBxVRCRxUlo8A/8EvumcW29mXYA3gH7An4HXnXNPm9kZwF+BS1Mcm4hITMws6BBSJZCc3bWrL57XrVPxLCLpJ9XDNvoBzUPPdwMtzawjcBnwcmj7m8B5ZtYtxbGJiFTLORd0CKkWSM4OD91YuzZRRxQRSZxUF8/PAY+ZWXPgCuBaYCRQ6JzbBeCcOwCsBkakODYRkZjUo57nQHJ2jx7+ccWKRB1RRCRxUl08XwPsAT4GvnbOPQ90Aooq7LcD6FjxzQUFBeTn5x9sEydOTHrAIiJh8fY8T5w48WC+AvKSElRyBZKzBw70j0uW1D5wEZF4xZqzUz3muTHwJD75/snMVgMO/3VgpFyg0hJe7dq1Y86cOUkPUkSkOrH2PE+YMIEJEyaE31MXl/0IJGcPGOAfFy+O+60iIrUWa85OdfE8BfiWc67YzBw+KV+Pv1s7UnNgQ4pjExGpVj0c8xxIzu7f3z8uXQr790NobRoRkbSQsmEbZpYHDHHOFYc23QlsB9YCHUNTH2FmOUAXYFaqYhMRiUd9GPMcZM5u3hy6dYN9+zTuWUTSTyrHPG8B9phZp4htRcBCYBpwWmjb6cDzzrn1KYxNRKRG9aznOdCcPXSof/z440QeVUTk0KVs2IZzrtTMzgduM7O5QHvgJufcdjP7X+AeMxsCdAAmpCouEZF41Yee56Bz9vDh8MIL8P77cNVViT66iEjtpXTMs3NuAfCDKNu3At9PZSwiIvGqZz3Pgebs4cP943vvJfMsIiLxS/VUdSIidVZ4We4GDVJ9r3X9M2wYNGzop6vbsiXoaEREyqh4FhGJ0d69ewHIzc0NOJLMl5tb1vs8fXqwsYiIRFLxLCISIxXPqXX22f7xlVeCjUNEJJKKZxGRGO3btw9Q8Zwq4eL5tdegpCTYWEREwlQ8i4jEKNzz3LBhw4AjqR/69oXu3aGoyM+6ISKSDlQ8i4jESMM2UssMxo71z6dMCTYWEZEwFc8iIjFS8Zx63/62f3z6adizJ9hYRERAxbOISEwOHDhAaWkpWVlZmqouhQYOhKOPhm3b4OWXg45GRETFs4hITLZv3w5AixYtAo6k/gn3Pj/6aLBxiIiAimcRkZhs27YNgJYtWwYcSf1z5ZXQpAnMmAGLFgUdjYjUdyqeRURioOI5OG3awHe+45//8Y/BxiIiouJZRCQGxcXFgIrnoPz4x372jSlTYNWqoKMRkfpMxbOISAzU8xysnj398I39++G224KORkTqMxXPIiIxUPEcvNtvh4YN4YknYMGCoKMRkfpKxbOISAy2bNkCQKtWrQKOpP7q1g2uuQacg2uv9Y8iIqmm4llEJAYbN24E4PDDDw84kvrtttvgsMPg3Xfh8ceDjkZE6iMVzyIiMQgXzx07dgw4kvqtdWu47z7//MYbobAw2HhEpP5R8SwiEoMNGzYAKp7TwWWXwahRUFQEEyZo+IaIpJaKZxGRGISLZw3bCJ4Z/O1v0KIFvPACTJoUdEQiUp+oeBYRiYGGbaSXbt3g4Yf98+uug2XLAg1HROoRFc8iIjXYtm0bxcXFNGrUiDZt2gQdjoRcfjlceins3AnnnQehdWxERJJKxbOISA1WrFgBQK9evcjKUtpMF2YwcSIMGgTLl/tCuqQk6KhEJNPpt4CISA2WL18O+OJZ0kuzZjB1KrRtC6+/DjfcoBsIRSS5VDyLiNQg3PPcu3fvgCORaLp3h+eeg5wceOAB+NWvgo5IRDKZimcRkRqEe55VPKevkSPh6achOxvuuAPuuSfoiEQkU6l4FhGpwcKFCwHo169fwJFIdS64AP7+dz8W+pZb4He/0xAOEUm8hBTPZnZWIo4jIpJudu3axaeffkp2djZDhgwJOpyEyOScfcUV8MgjvoC+9VY/Brq0NOioRCSTNIhnZzM7A/gF0IGywjsLaA80TmxoIiLBW7RoEaWlpQwcOJDGjetWmquvOXvCBGjVyhfSf/oTbN7se6QbNgw6MhHJBHEVz8A/gf8D5gPhv+UNuCyRQYmIpIs5c+YAMHTo0IAjqZV6m7PHjfMzcJx/Pjz5JHzxhb+psH37oCMTkbou3uJ5lXPuvoobzWxpguIREUkr7777LgDDhw8POJJaqdc5+7TTYNYsOPdcmD0bhg2D55+H444LOjIRqcviHfP8v2Y2Nsr28xIRjIhIOnHOMXPmTABGjhwZaCy1VO9z9tFHw9y5cNJJsH49nHwy3HuvxkGLSO3FWzyfCDxmZiURrRT4SxJiExEJ1LJly9i8eTPt27enT58+QYdTG8rZQIcO8NZb8KMfwf79cNNNcOaZsHFj0JGJSF0Ub/H8G2AC0AvoEdF+ndiwRESC9+abbwIwYsQIzCzgaGpFOTukYUO4/3546SU/Fnr6dOjfHyZN0nR2IhKfeIvn14CpzrlVzrm1obYGeDTxoYmIBOvFF18E4Oyzzw44klpTzq7g3HNh4ULf81xcDOPHw+jRsGpV0JGJSF0R7w2DTwE/NbO3I7YZcCFwXawHMbPjgZOAlcBsYCdwD7ADaA3c4pzbHmdsIiIJs3XrVmbOnEl2djbnnHNO0OHUlnJ2FB07wiuvwBNPwHXX+SEdAwfCb37jh3ZoSjsRqU68xfMv8V//XQ2Ev+gy/ByiMSViM/sO0NM594uIbZOB151zT4fmJf0rcGmcsYmIJMy0adMoKSnhtNNOo02bNkGHU1vK2VUwg8svhzFjfAH95JN+LPTEifD738N55/l9REQqinfYxu+Bzs65bs657qHWDd+LUSMzGwFcjk/o4W0d8XOOvhza9CZwnpl1izM2EZGEeeqppwC44IILAo7kkChn16BdO98D/eqr0KcPrFjh54Y+7TT45JOgoxORdBRv8fyFc25HlO0dYnz/fcBnwINm9rqZnQiMBAqdc7sAnHMHgNXAiDhjExFJiPXr1/P666+Tk5PDJZdcEnQ4h0I5O0ZnngmLFsEDD0CbNvDOO36au+9+F1avDjo6EUkn8RbP11TcYGYXA3+q6Y1m1gcYCjzqnPsh8BbwOtAZKKqw+w6gY8VjFBQUkJ+ff7BNnDgxzvBFRGr2+OOPU1payje/+U3y8vJqfZyJEycezFdA7Q9Ue8rZccjJgWuvhZUr4frrITsbJk+G3r39jYW6qVAks8Was83FMUePmb0CfOicu8PMmgAPABcDc5xzp9Xw3guBic65vNDPjYEC/FRKFznnjovYdx7whHPu3shj5Ofnu/BSuSIiybB//3569OjBl19+yWuvvcYZZ5yRkOOa2VznXH5CDhb7OZWzD8HKlXDnnTBlCpSU+GL6qqvg5z+HHj2Cjk5Ekqm6nB1Xz7Nz7mxghpk9CHwC9AaGANFWsKqoARE3KDrndgMrgGz83dqRmgMb4olNRCQRnnnmGb788kv69evHmDFjgg7nkChnH5ojj4S//x2WLoWrr/bbJk3yPdFXXgnz5gUZnYgEpdri2cxOqdiAhsBafBL9FdAVuC2Gcy0EWppZ+4htB/AJt2OoVwMzywG6ALPi/jQiIofAOce99/rO0xtuuIGsrHhHtgVLOTs5jjzSD99Ytgy+8x2/bcoUGDYMRoyAF1/0PdMiUj9UO2zDzD4E+gHFQGlVuwGHO+dyazyZ2avAK865h82sFT7ZDgOmAFOcc9PM7EzgSufcZRXfX9e/AhSR9Pb8888zduxYOnTowOrVq2nUqFHCjp2KYRvK2amxdi08+CA8+ihsD81u3bOnnyP6O9+B5s2DjU9EDl11Obum4vks/Ni4r2o4wVjn3HMxBJIH3A/MwfdUPOqc+9TMWuMn3F+Dvwv8p865ryu+P1MTsYgE78CBAwwaNIilS5fy8MMPc801le61OyQpKp6Vs1No+3bfI33//WUzcjRvDpddBt//vu+ZFpG6qdbFc7rJ9EQsIsF59NFHmTBhAt27d2fp0qU0TPAyc0HcMBi0+pKzS0pg6lT44x9h9uyy7UOH+iL6ssugZcvg4hOR+CXshkERkUxUWFjIrbfeCsBvf/vbhBfOktmys+HCC+Hdd/1c0ddd5+eKnj8frrkGDj/c33A4ezbUof4qEamCimcRqfduvvlmtmzZwumnn17XF0WRgA0cCH/6E6xf71cuPPVU2L0bHn8cTj7ZT3H305/6IltE6iYVzyJSr73xxhtMnjyZhg0b8vDDD2NmQYckGaBRI7j0Unj7bb/k9623QqdOsGYN3H03DB7sC+3f/laLr4jUNYdcPJvZUDObYmY/MbNuhx6SiEhqFBQUcNVVVwHw61//mt69ewccUfIpZ6fekUfC//0frFsHs2bBD37gh3UsWeIXXOnZE44/3t94uG5d0NGKSE0S0fP8P8Dpzrk/AMPDc3+KiKQz5xzjx49n8+bNjBgxgptvvjnokFJFOTsgWVlwyinwyCOwcSNMmwaXXw5Nm8J//+uXBO/a1c/SceedvrjWGGmR9HPIs22EJsg/zDm33syygGbOue0Jia6C+nLntogk35133skvf/lLWrVqxcKFC+ncuXNSz5cus20oZ6efnTvh5Zfh+efhtdf8z2G9esF558GZZ8JJJ0FujbNzi0giJHS2DTM708wuDz0/FujlnFsP4JwrTVYSFhFJlH//+9/88pe/xMz45z//mfTCOUjK2emvaVP41rfg2WehoABeegm++13Iy/Pjpf/wBzj9dGjbFr75TfjLX/zYaREJRlzFs5ndBTwLXADgnPsI+LaZVVpZSkQkHc2bN48rr7wSgLvvvptzzjkn4IiSRzm77mncGM49F/72Nz+0Y+ZMuPlmGDTI90i/9JKf/q57d+jXD378Y5g+HfbsCTpykfoj3p7n4cARwLsR214D7k5YRCIiSbJkyRLGjBnDrl27uPLKK7npppuCDinZlLPrsAYNYMQI+N3vYOFC+PJLeOwxuOgiaNECli710+J94xv+BsQzzoB77oE5c/zCLSKSHA3i3H+Gc67YzCIHSo9NZEAiIsmwcuVKRo8eTVFREWeddRaPPfZYfZiWTjk7g3TqBOPH+7Z/P3z4Ibz+uh8nPX8+vPGGb+BXNBw5Ek47DUaNgv79IfP/uYukRrzF8yoz+xnQycy+BVyM/zrw2oRHJiKSIMuWLWP06NFs3LiRU089leeee66+rCKonJ2hcnL8oisnnwx33QWbN8M77/h5pd96y88dPXWqbwCHHeaL6FGj/MItPXuqmBaprbhn2zCz44GrgS7ARuBJ59ybiQ+tMt25LSLx+vjjjznrrLMoLCxk+PDhvPHGGzRr1izlcQQ124Zydv20Zk35YnrjxvKvH364L7xPOsk/DhrklxkXEa+6nJ2IqeqygZHOubcO6UAxUCIWkXhMnz6dCy+8kJ07d3LmmWfy7LPP0rRp00BiSaOp6pSz6xnnYNkyX0i//ba/CbGoqPw+LVrAiSeWFdPHHONXSRSpr6rL2dUO2zCzl4CaumjygEIg6YlYRCQWzjkeeOABfvKTn1BSUsIVV1zBpEmTyMnJCTq0pFLOlmjMoG9f3665xhfTS5fCu+/C7Nn+cc0aP3b6tdf8exo29AX0ySfDCSfAccdB+/aBfgyRtFHTmOeFoX2WAaXA94H5QGRXQnegTVKiExGJ065du/jBD37AlClTALj11lu56667yMpKxIKqaU85W2pk5qe569cPJkzw27780hfR4YJ68WJ47z3fwrp08UV0uB19NDRpEsxnEAlSTcXz7/FDO4oBzKync+62ijuZ2dRkBCciEo8VK1ZwySWXMH/+fJo2bcqkSZMYN25c0GGlknK21MoRR8Cll/oGsHWrL5xnz/ZLh8+ZA+vW+fbss36f7Gw/VvqYY2DoUN8GD1ZBLZmv2uLZObetwqb2ZpbtnDs4g6SZjQSOTUJsIiIxcc4xceJEbrjhBnbt2kXPnj158cUXGThwYNChpZRytiRK69Zwzjm+gZ83+rPPfCH90Uf+cdEi+OQT38KysvzwkHAxHW6tWwfzOUSSId6p6p4DFpjZbGAH0A/4BlCpZ0NEJBU2bdrE+PHjefXVVwG4/PLLeeihh2jVqlXAkaUF5WxJiOxsGDjQt/Hj/badO2HuXJg3z88zPW+eL7A//dS3f/2r7P3dupUV0kcf7R8PP1zT5UndFFfx7JybYWZjgGuAwUAxMNY5p68ARSSlSktLmTx5MrfccgtFRUW0atWKRx55hEsuuSTo0NKGcrYkU9OmcMopvoXt3u3HS8+fX1ZQL1zob0hcswZeeKFs37w8P+wjsg0YAAHMJCkSl3h7nnHObQB+Ef7ZzLLMbJhzbm5CIxMRqcInn3zCNddcwwcffADAmDFjmDRpEp06dQo4svSjnC2p1LixHwN9zDFl2w4c8FPlhQvqcCss9HNRv/NO+WP06FFWTA8e7B+PPNIvVy6SDuL6p2hmk6JsbgeUAOcnJCIRkSps3bqV22+/nQcffJDS0lI6dOjAH//4Ry655JL6sNR23JSzJR00aOB7lAcMgCuu8Nuc8zN8LFpUvn32mV8dMbxCYlhurl9ivGJPtYZ+SBDi/TuuF35u0MiVVY4A3k1YRCIiFezatYsHH3yQu+++m+LiYrKysrjuuuu4/fbbadmyZdDhpTPlbElLZtC5s29nnVW2ff9+WL7cD/WILKrXri3rsY7Utq0vogcO9MV1uLVrl9rPI/VLvMXz95xzyyI3mFlf/NKvIiIJdeDAASZPnsyvf/1rNmzYAMCoUaP4wx/+wFFHHRVwdHWCcrbUKTk5Zb3U4WnzALZv92OpIwvqhQv9SokzZ/oWqW3bskK6X7+y5x07qqdaDl28Nwwui7J5L/A94NaERCQi9d6+ffv417/+xd13383y5csBGDp0KHfffTejR4/WEI0YKWdLpmjRAoYP9y3MOVi/3hfSixf7IR/h2T6KisoWfal4nHAxHVlUd+3qp9kTiUW8Y55XU/7rv2ygPfDPRAYlIvXTzp07eeyxx7j33nv58ssvAejZsyd33nkn48aNqy+rBCaMcrZkMjO/uMsRR8CZZ5Ztdw42bPBFdOTUeeGi+r//9S1S48Z+furIorpvX+jZ0y9VLhIp3mEbM4AnKEvGpcBm59zyhEYlIvVKQUEBf/nLX3jggQcoKioCoH///tx6661861vfIicnJ+AI6yzlbKl3zKBTJ99Gjy7/WkFBWSEdWVhv3Bh9THVWFnTvDr17Q58+5R87ddIQkPoq3uJ5inPuP5EbzKyPmR3pnFuZwLhEJMM55/joo4946KGHeOaZZ9i3bx8Axx9/PD/96U8555xz1NN86JSzRSK0awcjRvgWqbi4ci/18uV+burPP/fttdfKv6dJE19Eh1u4qO7dG7RGU2aLt3g+A/hPhW2bgXeAoQmJSEQy2u7du3nqqad4+OGHmTvXTzVsZpxzzjn85Cc/YcSIERrTnDjK2SIxaNUKTjjBt0h79vhp85Yt88X08uVlzwsKKi9PHnbYYZUL6j59/BzWubmp+UySPDEVz2b2Q+BOoLmZ3RJllxkJjUpEMopzjjlz5vD3v/+dJ598kq1btwLQtm1bxo8fz//8z//QvXv3gKPMHMrZIonRqFHZTYUVbd1auaAOt6++8m327PLvycryNyf26lW5devmZxuR9BdT8eyce8jM/gNcBTxQ4eV9zrmNCY9MROq89evXM2XKFB5//HE+++yzg9vz8/P54Q9/yLhx42jcuHGAEWYm5WyR5GvdGo47zrdIpaV+FpBovdVr1sDq1b5Nn17+fQ0a+AI6WmHdtStkZ6fqk0lNzDlX817hnc1ynHP7o2zv4pxbl9DIosjPz3dz5sxJ9mlE5BAUFxfz0ksv8cQTTzBjxgxKS0sBOOyww7j88su56qqrGDJkSMBRpp6ZzXXO5af4nMrZImlk714/DGTFisrtiy+qfl9Ojh/yEa2w7txZ0+wlQ3U5u9qeZzO7EXjVOfdpaNPNZlbxbx8DTgVGHmqgIlI3hQvmZ555hunTp7N/v6/XcnJyuOCCC7jqqqs444wzNGtGkilni6S33Fw/FV6/fpVf273b35gYrbDesMH3Xi+LMnN7bq6fUi9aYd2xowrrZKhp2EZ3IPKe0aOBpsAmyqY+MqBD4kMTkXRWVFTEtGnTePbZZ8sVzFlZWYwcOZJx48Yxbtw42rZtG3Ck9Ypytkgd1bixX2Z84MDKr+3cCStXRi+sN28umyEk2jGPPDJ6Yd2hg6baq61qi2fn3P+rsOkuYLFzbl/kRjOLMpReRDKJc46lS5fy8ssvM23aNN57772DQzKysrI49dRTufjii7nwwgtp3759wNHWT8rZIpmpaVMYMsS3irZvL19YRz4vKChbzryiZs2qLqzbtVNhXZ14p6ob5JybF7nBzPoAPYEof/NUZmYNgY+B65xzM82sMXAPsANoDdzinNseZ1wikgT79u3j3Xff5eWXX+bll19m1apVB1/Lyclh1KhRXHTRRVxwwQUqmNPTIefs0HuUt0XSVIsWcPTRvlVUXFx1j/WWLVVPtdeiRfSiundvaNMm+Z8p3cU6Vd3J+K/6RoeWe43UBn8398sxnvMmoFvEz38GXnfOPW1mZwB/BS6N8VgikkDOOZYtW8aMGTOYPn06M2fO5Ouvvz74el5eHmeddRbnnnsuY8aMoUWLFgFGK1VJcM4G5W2ROqlVK8jP962iLVuiF9UrVsC2bTB3rm8VtWtXtoR5ePx2//5+fHV96a2Otee5GPgHfjzdiRVe2wc8FstBzGw4sBHYGvq5I3AZEP6q8U3geTPr5pxbE2NsInIICgsLefPNN5kxYwYzZszgiwq3fA8YMIBzzz2Xc889l+OOO45szZdUFyQkZ4PytkimatMm+lR7zkFhYfSiOrw4TEEB/KfC8kstWkDfvpWL6m7dMm+avVjneV5kZiOBi51zMSfdSGbWFBjrnLvBzG4LbR4JFDrndoXOcyDUSzICWFOb84hI9bZs2cJ//vMfZs2axaxZs/jkk0+InLKyXbt2jB49+mDr1KlTgNFKbSQiZ4Pytkh9ZOZ7l9u1g+HDy7/mHHz5ZdlS5pGPRUXw0Ue+RWrUyA/3iCyqBwzw2xrEO3g4TcQctnNuG1X0VpjZSc652dFei3ALcHeFbZ2AogrbdgAdox2goKCA/IjvHiZMmMCECRNqOK1I/VZYWHiwWJ45cyaLFi0qVyzn5uZy0kknMWbMGMaMGcPgwYPJ0txGCTFx4kQmTpwY/jEvledOQM6GQ8zbytkimcXMzyvduTOMGVP+tYKCygX1Z5/5BWMWLvQtUm6uL6KHDIHBg8seg5ygKdacXe0iKWb2EtCshnMZMMA5d1g1x/kGkOuceyn08xrgaiAf3zNyXMS+84AnnHP3VjyOJtwXqZ5zjuXLl/PBBx/wwQcf8P7777N48eJy++Tm5nL88cczYsQIRowYwfHHH0+TJk0Cirj+SMUiKYnK2aFjHXLeVs4WkW3bYOnS8kX14sV+tcVoOnUqX1AfdZTvpU51n06tF0kBluHHzq2nbI7QSscHzqnhODcCQ61sJHlrYCpwX+h5pObAhhqOJyLA119/zccff8z777/PBx98wIcffkhRUflOwUaNGnHCCScwYsQIRo4cyXHHHUejRo0CiliSLFE5G5S3RSQBWraMPrZ62zY/hd6CBb4tXOh/Xr/et1dfLdu3RQs49lg4/viyY7Vrl9rPEammnuc8YIdzbm+1BzFr65yr+DVe5OvtgNyITR8ANwD/xU+X1M45t9vMcoCvgR7OufUVj6NeDKnPnHN8/vnnfPjhhwd7lRcuXHhwruWwDpvB9XQAACAASURBVB06MHz4cE444QROOOEE8vPzyc3NreKokiop6nlOSM4O7XPIeVs5W0TiUVLily8PF9MLFsC8eX6cdUU9e8KoUX74yKhRiZ9Cr9Y9z865whjPcQHV3L3tnCuoEFAJUOCcW2dm04DTgGnA6cDz0QpnkfrEOcfatWuZM2fOwTZ37lyKi4vL7Zednc2wYcM44YQTDhbMXbt2xerLfEFSTqJyduhYytsiklLZ2WVzSo8dW7Z9/Xr47399+/BDmDPHL2X++efw6KN+LHZ+Ppx1Fnz729CjR3LjPKT7HM2sC366oh8Sx9RHFfwvcI+ZDcEvGau7SaRecc6xfv36coXynDlzKg2/AGjfvj3HHnvswUI5Pz+fpk2bBhC11EUJytmgvC0iKdSpE1x4oW8ABw7A/PkwY4Zv770HH3/s2+23w4gRMH48fOtbkJOT+HiqHbYR9Q1m2cB5wPeB0cDm0HGSPp+VvgKUuq6kpIQVK1awYMECFixYwCeffMK8efPYvHlzpX3z8vLIz88v1zp27Khe5ToqFcM2qjivcraIZLSvv4ZZs+Cpp+D552H3br+9b1+YMgWGDYv/mIdyw2DkQboD3wO+A7QF1gLnAm8AJ8Qflkhm2759OwsXLjxYKC9YsIBFixaxO/x/dYTWrVsfLJCHDRtGfn4+Xbp0UaEstaacLSL1RbNmcPbZvj30EDzzDNxzj5/lY8QIP7tHt26JO1+NxbOZXYzvsRgFrMTfaf048DPn3Ouh3d5PXEgidUt4fHJkb/KCBQtYtWpV1P07d+7MUUcdxZAhQxgyZAhDhw6lR48eKpQlIZSzRaQ+a9kSvv99P/b56KP99HgvvgjXX5+4c8TS87wbv5zr+8D3nHPLAcwsvvEeInWcc46vvvqKJUuWsHjx4nKP27Ztq7R/w4YNGTBgwMEi+aijjmLw4MG0SfQtwSLlKWeLSL22eDH8/ve+cAY/V3Qi1Vg8O+emAdNCN5p838y64ef6PDgE28zMxTt4WiSNFRUVsWTJkkqFcrSb+MAvaR1ZJA8ZMoS+ffuSk4w7FUSqoZwtIvXN3r3+psHp0+H11/0Ud+CX/77nHhg5MrHni2d57nXAL82sAXA+cISZ3QH8DRgDTKzu/SLpaNu2bQeL5MhCedOmTVH3b9GiBQMGDGDAgAEMHDjw4PMOHTpo2IWkFeVsEclEpaWwciXMnVvWPvoIdu0q26dlS7jsMrjxxuRMWxf3VHXOuQPAc8BzZnYk8GP8+DolYklLzjk2bdrE0qVLWbp0KZ999tnBxy+jzbwONG3alP79+1cqko844ggVyVKnKGeLSF104IBfwvuzz/yNf+Elvpcsge3bK+8/eLBfMGX0aDjlFEjmQrqHNM+zc24lcJ2ZzU1QPCK1tn//flatWlWpSF66dGnUMcngl67u169fpSK5a9euZGVlpfgTiCSXcraIpAvn4KuvfIG8erV/rPh8377o7+3Y0U8/N2yYvynwmGOgQ4eUhX5oxXOYc+4fiTiOSCy2b9/OsmXLKhXIK1euZP/+/VHf07p1a/r160ffvn3p27fvwefdu3cnOzs7xZ9AJFjK2SKSTM7B1q1+ZcANG8o/rl1bVhxHmbm1nM6d/VzNfftCv35lj6kslKNJSPEskmjhVfeiFcnr11e9CnDXrl2jFsnt2rXTcAsREZFDUFICRUWwebPvNd68GTZurFwkb9gAe/bUfLzWraF7dz8Hc/gx/LxrVz9/czpS8SyB2r17NytWrDhYJC9dupRly5axbNkyvv7666jvyc3NpXfv3pWK5N69e9OkSZMUfwIREZG6a8+e8sVwdY+Fhf6GvVi0aOGHV3TqVP6xS5ey4rhly+R+tmRR8SxJ55xj8+bNBwvjyCJ5zZo1VDVjVl5eHn379qVPnz4HC+V+/frRtWtXDbUQERGpYN8+3zMcboWF5R8rbtu8GXbsiO8cbdpA+/Zw2GG+dewYvUhO117jRFDxLAmzb98+Vq5cWa5ADhfJVd2wl52dTc+ePQ8WyeGe5D59+tC2bdsUfwIREZHgOeenXtuyJbZCOPw83kIYICenfDEcfh7tMS/P71/fqXiWuBUWFkbtRV61ahUlJSVR39OqVatyhXH4eY8ePWjYsGGKP4GIiEjy7dvnb5zbssU/Rj6vaVtVM01UJzsb2rYta3l50R/Dz9u390MndEtQfFQ8S1QHDhw4OO1bxSK5qlX2srKy6NGjR9QiWTfsiYhIXVRSAtu21a4A3rmz9ufNzfVDJCoWvNUVwy1agGZZTT4Vz/Xc1q1bDxbHkUXy559/XuW0b82bN680xKJv374ceeSRNErmrOQiIiK1UFICxcVlRW08bft2P4yiNrKz/YwSbdr4x8jnNW1r3Dix10ASR8VzPVBSUsLatWuj9iJv3ry5yvd16dIlai/y4Ycfrl5kERFJqUMpgKu47SZmLVvWrgBu3lxDIjKRiucMsmPHjqi9yCtWrGDv3r1R39OkSRN69+5dqUju1asXTZs2TfEnEBGRTHbgwKH1AB+Kli3LitrqWmTx27q1f58meJJIKp7roMLCQhYuXMiSJUvKzWqxYcOGKt/TqVOnqEMtjjjiCC1DLSIiMauuAI4c7xut1WY2iDCz2Avgik0FsCSSiuc0tmfPHj799FMWLVrEokWLWLhwIYsWLWLTpk1R9w8vHlKxSO7Tpw/NmzdPcfQiIpKu9u+vfQ9wogvgij29VbUWLVQAS3pQ8ZwGnHOsWbPmYHEcflyxYkXUqd+aNWvGwIEDGThwIP379z9YJGvxEBGR+iU8H3BhYflWUFD+56Ki8r3CVSzgGhMzaNWq9j3A+rJT6joVzynmnGPVqlXMmTOHuXPnMnfuXObNm0dxcXGlfbOysujbty+DBg1i8ODBDBo0iEGDBtGtWzcNtRARyUB791YuhGtqe/bEfx6z2hW/4R5g/QqS+kzFcxI551i9ejVz5sw5WCxXVSi3b9+ewYMHlyuS+/XrR2PNVSMiUqft2+eXQQ63TZvKP0Y+j/LroUaNGkG7dn6e36pa27YqgEUSRcVzAu3bt4/58+cze/Zs3nvvPd577z2++uqrSvu1b9+e/Px8hg0bdrB17NhR07+JiNQhBw7Axo3w5Ze+rV9f9rhpU1lBvHVr7MfMzvZLIVdXCFdsTZok7zOKSGUqng/Bnj17mD17Nu+88w6zZ8/mo48+Yk+F78/y8vI49thjGTZs2MGCWYWyiEh6c84PiVi1yrcvvqhcJG/aBKWlNR8rXBB36OCXQ27fvurnbdqoR1gk3al4joNzjkWLFjFjxgymT5/Of/7zn0rFct++fTnxxBMPtl69eqlQFhFJQyUlvjBeubKsSI5sNd1UZwaHHw5HHFG+depUVhB36KCCWCTTqHiuwYEDB5g1axbPP/88L774Ihs3biz3+tChQzn99NM5+eSTGT58OG3btg0oUhERiaakBD7/HJYsgU8/LXtcutTfoFeVli2hZ0/o3h26dq1cIB9+OOTkpO5ziEh6UPEchXOODz74gMmTJ/PCCy9QVFR08LVOnToxevRoxowZw2mnncZhhx0WYKQiIlLRunXw0Udlbe7cqnuRO3eGXr18kdyzJ/ToUdZat05t3CJSN6h4jrB9+3YmT57MxIkT+fTTTw9u79OnDxdddBEXXXQRQ4cO1TAMEZE0snUrvP02zJjh26pVlffp3BkGDoT+/WHAAP/Yr5+fdUJEJB4qnoHi4mLuv/9+7r//fraGbotu3749V199NVdeeSX9+/dXwSwikkZ27YKpU+GJJ+D11/3MF2GtWsHxx8Oxx8Ixx/jWvn1wsYpIZqnXxbNzjn/84x/ceOONFBYWAnDyySdz/fXXc+6555KjwWwiImllyxZ48EHfwiPqsrNhxAgYPdq3YcO0jLOIJE+9LZ6Li4u54ooreOWVVwA45ZRTuOOOOzjllFMCjkxERCpyDp58Eq67zk8hB5CfD1ddBePG+angRERSoV4Wzxs3bmTkyJEsX76cNm3acP/993P55ZdraIaISBoqLYXrr/e9zQCnnAK33+57m5W2RSTV6l3x/PXXXzNmzBiWL1/OoEGDmDp1Kt27dw86LBERqcJtt/nCuWFDeOgh+N73VDSLSHDqXfF88803s3jxYvr06cPbb79NXl5e0CGJiEgV5s+Hu+7yi4y89BJ84xtBRyQi9V1K1zwyswvNbLWZFZnZ/WbWILQ9z8z+ama/M7P7zCwpd+p99tlnPPLII+Tk5PDMM8+ocBYRqUbQORvg97/3j9ddp8JZRNJDyopnM+sCnA+MBX4MfA+4PvTys8Ak59wtwHLgN8mI4b777sM5x/jx4xk8eHAyTiEikhHSIWfv3w/Tpvnn116bjDOIiMQvlcM2ugLfdc4dAOaa2UDgVDObDfRxzv03tN/LwDIzu9M5tzNRJ9+3bx/PPfccAD/60Y8SdVgRkUwVaM4Gv6T2jh1+aWzdmiIi6SJlPc/OuXdDSThsPbAOGAWsjthvfehpfiLPP3PmTIqLixk0aBD9+vVL5KFFRDJO0DkbYOVK/9inT6KPLCJSeykd81zBMcAjQCegqMJrO4COFd9QUFBAfn7+wTZx4sSYT/bOO+8AcOaZZ9Y6YBGReEycOPFgvgLq+k0WKc3Z4BdEAc3hLCKpEWvODmS2DTM7Eihwzi0wMwfsrrBLLrC/4vvatWvHnDlzanXOWbNmATBixIhavV9EJF4TJkxgwoQJAJhZYcDh1FoQORtg2zb/2KpVrQ8hIhKzWHN2ynuezSwb+AFwa2jTeqB1xOsGNAU2JOqce/bs4eOPPyYrK4uTTjopUYcVEcl4QeTssHDx3KJFoo8sIlJ7QQzbuBG41zm3N/TzDCByRFs3YB9Q++6KChYvXsyBAwfo06cPLZSFRUTikfKcHbZnj39s3DjRRxYRqb1Uz/P8M2A+0NjMepjZd4FiYGPoTm6AM4AHnXP7EnXeefPmAXD00Ucn6pAiIhkvqJwdFi6eGzVK9JFFRGovZWOezewXwB0VNi91zk0ys0uA28xsDdAM+Hkizz1//nwAhg4dmsjDiohkrCBzdtju0MhqFc8ikk5SVjw75+4E7qzitbXA+GSd+5NPPgFUPIuIxCrInB2mYRsiko6CnKouZZYvXw5A//79A45ERERipWEbIpKOMr543rJlC1u2bKFZs2a0b98+6HBERCRGKp5FJB1lfPG8MrRE1ZFHHomfUUlEROoCjXkWkXSU8cXzihUrAF88i4hI3aGeZxFJRxlfPK9duxaAHj16BByJiIjEQ8WziKSjjC+eN2zwi1516tQp4EhERCQe+0IzRzdsGGwcIiKRMr54Xr9+PQAdO3YMOBIREYlHSYl/zM4ONg4RkUgZXzyHe55VPIuI1C2lpf5RxbOIpBMVzyIikpbCPc9ZGf+bSkTqkoxOSaWlpWzcuBGAww8/POBoREQkHup5FpF0lNHFc1FRESUlJbRu3Zrc3NygwxERkTio51lE0lFGp6TCwkIA8vLyAo5ERETipZ5nEUlHGV08FxUVASqeRUTqIvU8i0g6yuiUFC6e27ZtG3AkIiISr3DPs4pnEUknGZ2SVDyLiNRdGrYhIulIxbOIiKQlDdsQkXTUIOgAkql79+6cd955HHXUUUGHIiIicRozBoqKoEmToCMRESmT0cXz2LFjGTt2bNBhiIhILUyeHHQEIiKV6cswEREREZEYqXgWEREREYmRimcRERERkRipeBYRERERiVHGFs8TJ04MOoSMo2uaeLqmiadrWrfpv1/sdK1ip2sVG12n2Kh4lpjpmiaermni6ZrWbfrvFztdq9jpWsVG1yk2GVs8i4iIiIgkmjnngo4hZmZWAKyNcfc8oDCJ4dRHuqaJp2uaeOl6Tbs659oFHUQqxZmzw9L1v1860rWKna5VbHSdylSZs+tU8SwiIiIiEiQN2xARERERiZGKZxERERGRGGVc8Wxmjc3sQTP7rZn9xcxaBB1TXWdmF5rZajMrMrP7zaxB0DFlCjNraGYLzGxk0LFkCjM73sxuNLPzzSwv6HgkNsrdVasqB5tZnpn91cx+Z2b3mVlO0LGmi4q5Vf++qlYxZ+pa1Szjimfgz8Bs59zPgKnAXwOOp04zsy7A+cBY4MfA94DrAw0qs9wEdAs6iExhZt8BznHO3euce9E5pxtf6g7l7ihqyMHPApOcc7cAy4HfBBJkeqqYW/XvK4oqcqauVQ0y6oZBM+sIrAZaO+d2hf463wYMcM6tCTS4OsrMTgY+cM4dCP18D/56nh1sZHWfmQ0H+gK3AVc752YGG1HdZmYjgF8Co10mJbZ6QLm7alXlYOAO4N/OuY6h7Z2AZUB759zOoOJNBxVzK/4PC/37qiBaztT/i7HJtJ7nkUChc24XQCjZrAZGBBlUXeacezectEPWA+uCiidTmFlTYKxzblLQsWSQ+4DPgAfN7HUzOzHogCRmI1HujqqaHDwKf43C+60PPc1PYXhpp4rcOhL9+4omWs4cia5VjTKteO4EFFXYtgPoGEAsmeoY4JGgg8gAtwB3Bx1EpjCzPsBQ4FHn3A+Bt4DXzax9sJFJjJS7YxfOwbpm0UXLrbpWFVSVM4HO6FrVKNOKZwfsrrAtF9gfQCwZx8yOBAqccwuCjqUuM7NvAHOcc18FHUsGGQBscc4tDP38EGDAN4MLSeKg3B2DCjlY16yCanKrrlVlVeVMXasYZFrxvB5oXWFbc2BDALFkFDPLBn4A3Bp0LBngRmCSmRWaWSH+L/2pZnZzwHHVZQ1CDQDn3G5gBdA2sIgkHsrdNYiSg8tdMzMzoCn1+5pFza3466J/X+VVlTOz0bWqUaYVz+8AHc2sMUBo2p4uwKxAo8oMNwL3Ouf2Bh1IBrgMOCqibcDfQa/hMLW3EGhZYZjGAeDTgOKR+Ch316xiDp4B9Il4vRuwD5iT4rjSSVW5dTL691VRVTlzA7pWNcqo4tk5twmYBpwW2nQ68HzEjRRSC2b2M2A+0NjMepjZd0NfH0otOOcKnHNfhhtQgv8qdnvQsdVVzrmlwGv46bwws1ZAQ+CVIOOS2Ch3Vy9aDgaKgY1mNjC02xnAg865fUHFGbRqcus69O+rnGpy5r/QtapRRk1VB2BmrYF7gDVAB+CnzrmvAw2qDjOzX+CnRIq01DnXL4h4MpGZrUFT1R2y0IIo9+N73rrgb4RRz3MdodwdXXU52My64qdjWwM0A35eYWaOei0yt+rfV2VV5Uxdq5plXPEsIiIiIpIsGTVsQ0REREQkmVQ8i4iIiIjESMWziIiIiEiMVDyLiIiIiMRIxbOIiIiISIxUPIuIiIiIxEjFs9RrZnaamb1vZiNr2G+gmT1jZr9M4Ll/n6hjiYgEwcwGm9mzicyNVZxnmpl9K5nnSBYzO93MPqjp94zUHQ1q3kUk+cysH3AXsAU4BhgMXOmcm5LkU38BDIlhv3VAayA7EScNLXm6PxHHEhEBMLNjgP8DTgGmA6vxi1zkAH9xzr2RhNOuBFqSoNxYjT8Di5Nx4FA+vhB4BFgOfBB6KRe4CDjWObfmEE4R6+8ZqSNUPEu6eBEY55xbYGZZwORUnNQ5t9zMimLYb7uZbUzgqUcA/0ng8USknnPOfWxmTwF9nXPnhLeb2dnAC2Z2o3Puzwk+564E50YAzOxcYFG4aHXOvZroc4Q55/YDT5vZtcBM59wvIuKYWptjmtlhwKnOuaedc8vMrDBB4Uoa0LANCZyZtQF6AzsBnHOl+CVnU7X8ZWmC94vFqcA7CTyeiAhApeW5nXOvAL8G/mhmHZJwzoTmajM7At/TnGrRljZ/E1gfz0HMLBf4B9A4YnMif39IwFQ8S+Ccc1uAZcDrZnZ8aNta4CMza2Vmd5jZR2Z2npltMrPPzOzo8PvN7FIz+7mZzTSz30ZszzezX5vZ02b2bzNrGdreyMweNrPbzOxv+K8c42JmN4diGmdmx5jZi6Hj3RGK8Wkzaxfavs3MflXhEA2dc3tDx/qJmd1gZu+Z2d1xX0ARkZr9DT9842IAM+tpZr8xs8fM7K1QwYqZ5YZy2c/N7F0zuyR8ADO7xsx+aWYfmtkPqzpRtJxsZmNDOfNqM1thZnebWUczmxzKp++a2ajQIS4DjgBuMbPzQ/emvGlmV0Wc4/tmdpeZ/cPMHjezpmbW0sx+ZWZzzGyUmS0zs1Vm1q02Fyx0TUaGeqYxswvN7H4z+4OZvW5mzUKf4X4zezn0+2Zj6Br3By42sxsiDtk69Ltoh5l9v4pzZpnZ78xsQugczsxeCP/3kTThnFNTC7wBA4C1+L/OJwGHhbZn4xPRFvzYs1bADGAJ/o+/E4BfhPZtje+9PhtoATwVcfz3gN+Hnv8RmBB63h4owSfImmL8O7735jD82Lgmoe25wBvAVKAjPunvB34FNAGGhX4O7x8ezx3+3H8OPW8O3Br0fws1NbW624CrgS+reG0z8EAor/4byA5tfwJ4NvT8PuDC0PPvAgWh55cCV4Se9w7l6kGhn/8O/Dr0vKqc3AYoBv4EHAucDtwL3Bba9ybg5YhYHdAt9Lwx8BlwdejnbwBPROw7DfgrfijqBaHfFxeGPufrwF0xXruZod8VdwO/x49/Pj3i9UKgS+j5EvzvpFzg58CXwHHA+FDenxmON7T/GuAPodcuBdZVEcP3gGmh5+1D13lM0P+u1Mo3jXmWtOCcW2Jmg4A7gWuAc83sdOfHQBcA251zzwOY2a+B2UAP/C+KUjO7PnSoN/BJ+mygbcT2z0LvbQb8b+i9OOc2m9kXcYSaBzwIjHfO7QodY6+ZbQJWO+c2hM6zCZjl/HjA+fiknoe/8fBsfKIH2AV828yWA3/B/+EgIpIMDl9QHgt0Ba41M4AioKX5H36A/8Mf4HH8sAXwuXZOhVx7OLCowjmuJkpOds5tMbNtwFTn3EcAZrYO2Glmh+N7aptFDdq53Wb2VcSmH4aOG/YY8ExoezH+98W/Q+eYE4ozVu+40JhnM3sa32ETdrFzbp2ZnYgvmpuF8v8G4HPn3H+B/4beG+3YL4d+J3xcTUxHA7tDn3uzmW2h/PAPSQMqniVwZpYNdHb+xpAfhYZSvApMxP8lX9GS0GMroAvwT+fcE6Ftfwod8xZ8MfunCucahk96O2oZbm/gJKAnsCBie8UxfyUHX3CuNJRIw8Ok2jo/VAXn3Gozuxp4CPgxcDnwFSIiCWRmrfDfmi3B582iKPmxPb5ntCGAc64E/wc/offc7px7P/RzufdGiJqTQxzlc+WXwC3ACuBDoHs1HyHyfb2AtyN+/hw/JOUwKufiA9R+JpB5+D8yworN7HfAs/jYwxVyxc9VkwNUXX+9BvzFzBrh/1vsxvdiSxrRmGdJBwfH4QE45xYAPwP6VrF/Q3yiWglsBMZGvmhmx4W2n21mjSts3x76sV8tY30fuB/4Vyi5xcXMOhJx84mZdcLPNNIbfwPhE1W8VUTkUFyF/6brWXx+PDHU4wsczI+FwD7gjIjtg8ysORVyrZk1CHVGVFRVTo7mEWCN81OS7o3js6zD91QfPAWwDT8sJWGcc6WhDo4BoT8+XgV+45ybk8jzVDjny8DL+G9ILwVGOOe2Jet8UjsqniVdXBsqLMM64ceqhbUxs4ah5yPx492KgSeBC8zsITMbGeoV2IlPck2BaWb2jdDNGX3xBfd84A4za2hmefixxoeZWbaZ3WpmV1QRo4XabfhfMPdFea3i/hV/Pgc/Pi+sF3C+c247cG2U94iIxKNSj6aZnQLcDvzAOVeA7+XdDLxqZueGbgo8OdTT/BR+Vo7zzew04Crn3A58rv2Rmd1uZiPwnQjhIW+R+a+qnBwW2Qs8FGhnZk2AE4HGZtYj9Np+/A124U6UyHM8AlxkZq1DPx8HPOqcO4Cva6LmYjM7PhRXbhXXLosKdVHoD4yz8d82tsYPB+wDdAaaRNyMWLF3e1+U+MsHFWVsR+i/RQvgLXzvejx/VEiKaNiGpIvOwBIzewX/lZYB/xPx+gHg52ZWjO+lvRbAOTfDzH6E/+rvfPwNd4sBzOybwMP4XwaTgJucc87MxuGnEVoCPI3/+u0o/KIC4/DDQcotzmJmx+JvhOmC/xrvReD2UDxvAMcD3cysJ75HpAMwzsyWUtYLczXQ1Dk3scJnn2x+cYMD+Bt0RETiFsojl+ELvEeBAvxNZy2BM51zHwA45/aF8uOj+Fw3FT/WGeA6fMfDP/DftIVnuPgbPvf9L3542f86574KFYfHAz3NbEpVOdnMxuLH+X7PzJY659bjb168F184/wP4JjAIWAX8E99L/m3zszANAM4ys1edc/82s17A82b2Hv7by5+HesgvBjqEPt8S/IIxh5vZgFCc38KP434x4rrl4HP/EKBTqLh2+LHGF4aO+QkwC/gYPxTlNXxOn4a/SXGImZ3vnAsf92n8gjVfh4rtDsBYM1uN71Em9P6Kaxqsw98QeT7QCMgyszfxNw2mavpWqYHpv4WkO/NLmv7dOdctRee7yTmnpbNFRDJMaAhJrnMuLRepMrOb8DOfrAn93Aj/rcHvwvfKSPA0bEPqgmhDIpJzIv/15rxUnEtERFIndEPk0WlcODcCfgF0i9icB+xU4ZxeNGxD0lpovNkV+K/hLnbOPZvkU34e+jpRREQySwPn3F+CDqIqzrk95pcIfyB0M/kG/BCTX1T/Tkk1DdsQEREREYmRhm2IiIiIiMRIxbOIiIiISIxUPIuIiIiIxKhO3TCYl5fnunXrFvf7CgoKaNeuXeIDqsN0TSrTNalM16Sy2l6TuXPnFjrn6tXFVM5OPV272tO1OzSZdv2qy9l1qnju1q0bc+bEvypmfn5+rd6XyXRNKtM1qUzXpLLaXhMzW5uEcNKacnbq6drVnq7docm061ddztawDRERERGRGKl4FhERERGJUb0onidMmBB0CGlH16QyXZPKdE0q0zVJPl3j2tO1qz1du0NT2ULydAAAIABJREFUn65fnVokJT8/32XSeBoRqT/MbK5zLj/oOFJJOVtE6qrqcna96HkWEREREUkEFc/x+vrroCMQEclodekbURGpf1Q8x+Of/4Tu3WHhwqAjERHJONve38a8E+exfMLyoEMREamSiud4vPkmFBbC+efDli1BRyMiklGymmSx/f3tbHlji3qfRSRtqXiOxyOPwLBhsHo1XHoplJQEHZGISMZoNrgZDdo2YO8Xe9n9+e6gwxERiUrFczwaN4Z//xvatYPp0+FnPws6IhGRjGFZRutTWwNQ/FZxwNGIiESn4jleXbrAM89Adjbccw/8619BRyQikjFandYKgK1vbw04EhGR6FQ818bIkfDHP/rn3/kOvPVWoOGIiGSK1qNCPc9vF+NKNe5ZRNKPiufauvZa+PGPYf9+uOAC+OSToCMSEanzGvdqTG7XXPYX7mf7R9uDDkdEpBIVz4fi3nvhkktgxw4480xYrumVREQOhZmR9808AIqmFgUcjYhIZSqeD0VWFjz+OJx2GmzaBKeeCitXBh2ViEidFi6eC18sDDgSEZHKVDwfqtxcmDoVTjkFNmyAUaP8VHYiIlIrLU9uSYNWDdi1dBe7lu8KOhwRkXJUPCdC06YwbRoMHw5ffAEjRmgIh4hILWXlZNH2nLYAFL6g3mcRSS8qnhOleXN47bWyAvqkk2D+/KCjEhGpk9qNbQfA5n9tDjgSEZHyVDwnUosWfvGUMWOgoMBPaTd7dtBR/f/27js8jurs+/j3VrMkW7a6Zbl3G1ONaIFg02sooSWUBEjoIfCQkAaBBHhSIS8BQsAJJCSUhx5qANMhEMAUU13ADcuWLVmSLXfLOu8fZ2Wt1iorabWz2v19rmsv787Mztx7NLp96+yZMyIifU7hEYVkFGaw7qN1rJ29NuhwRES2UfEca/37w+OPw0knwZo1cMgh/qYqIiIStbSsNEpPKQVgxd3qfRaRxKHiuTf06wf33QfnnQcbN/rp7H71K3Ca8F9EJFqDzxgM+OK5aUtTwNGIiHgqnntLejr8+c9www1gBldc4e9GuGlT0JGJiPQJA/ceSO7kXDZXbda0dSKSMFQ89yYzuOwyeOQRyM31c0JPm+YvKBQRkQ6ZGeUXlgOw7NZlAUcjIuKpeI6H446D116DESPgrbdg6lR4/vmgoxIRSXhlZ5SR1j+N+pfrWffpuqDDERFR8Rw3U6fCe+/BYYdBTY2fkeN//xe2bg06MhGRhJUxKIOyM8oAWHrj0oCjERGJc/FsZpPN7E9m9j9mdquZ7RpanmNmN5vZr8zsz2Y2MJ5xxU1RETz1FFx1lb948Mor4eCDNYxDRBJSouTsYZcOA4Oqu6rYVKnrRkQkWPHuef4n8Cvn3P8DfgPcF1p+K/C6c+5nwGPA7XGOK37S0+GXv/Q3VCkthZdfhp13hgcfDDoyEZFICZGzcyfmUnJCCW6z48s/qLNBRIIV7+J5MpAXer4BGGRm5cCpwBOh5c8Dx5jZqDjHFl+HHw4ffQRHHQX19XDyyXDmmf65iEhiSJicPeKnIwBYdvsytqza0puHEhHpULyL54eAv5pZHnA6cDEwHahxzq0HcM41AguBaXGOLf5KS+GJJ+DWWyE728/GscMO8OijQUcmIgIJlLPzpuZReHghTeuaWPLbJb15KBGRDsW7eL4Q2Ai8A6x1zj0MDAVWRWzXAJTHObZgmMEFF8D778NXvgLLl8PXvw4nnghVVUFHJyKpLaFy9qhrRwGw9KalbFyysbcPJyLSpngXzzn4MXP3Ajea2cGAw38dGK4fsN33ctXV1VRUVGx7zJgxo9cDjptJk/x0djffDAMGwMMPw+TJcPvtmpFDpI+aMWPGtnwFFAcdTzckVM4eWDGQklNKcJscC69a2KN9iYhEijZnm4vjLaPN7BngG865ejP7OfB94FLgaufchLDt5oeW3Rv+/oqKCjdr1qy4xRuYJUvg/PP9RYUAu+3mi+p99w02LhHpNjN71zlXEXQcXZGIOXvDFxt4e/LbuEZHxQcVDNh5QEz3LyICHefsuPU8m1kxsItzrvmKuOuANcBioNzMckLbZQIjgFfiFVvCGTHCT2l3//0wbJgf0rHffnD66VBZGXR0IpICEjVn54zNofyCcnAw/6L5uKb4dQCJiEB8h23UAhvNbGjYslXAh8CTwEGhZQcDDzvnUrtKNPMzcMyZAz//OfTrB/fcAxMnwrXXwtq1QUcoIsktYXP2qF+MIrM0k9Wvr6bqb7o2RETiK27Fs3OuCTgOuMrMzgWuBC53zq0BLgCONbMrgCOBc+MVV8Lr3x+uuQY++wyOPx7WrfM3WRk3Dv70J9i8OegIRSQJJXLOzizIZNz/GwfAFz/6gs3VyoMiEj9xHfPcUykz5rkjL78MP/4xvP22fz1mDFx3HZxyCqTpbusiiaovjnnuqd7M2c45Pjz0Q+qer2Pw6YOZ/M/JvXIcEUlNCTHmWWJk+nT473/9bBwTJ8KCBXDqqTB1KjzyCDQ1BR2hiEivMzPG/3k8adlprLh7BSsfXBl0SCKSIlQ890Vmfi7ojz+Gv/4Vhg6F2bPhhBNgl13ggQc0vZ2IJL3ccbmMvWEsAPPOm8fGpZr7WUR6n4rnviwjA77zHfj8c7jlFj8zx8cf+yEcO+0E996rIlpEklr5BeUUHlVIY10jc749R7NviEivU/GcDLKz4aKLfBF9++0wcqS/wPC00/zNV/78Z1i/PugoRURizsyYdMckMksyqX+xniW/0627RaR3qXhOJv36wbnnwvz5cMcdMHasL6gvvNAX1L/4BVRXBx2liEhMZQ3OYtLfJgGw8IqF1M6sDTgiEUlmKp6TUWYmnH02zJ3rxz/vsQfU1MAvf+lvwHLhhb7AFhFJEkVHFTHy5yOhCT79xqdsWBR5B3ERkdhQ8ZzM0tPhpJPgrbf8FHdHHQUbN/phHBMnwte+Bs89B31oukIRkfaM+sUoCo8spLG2kY+P+5jGhsagQxKRJKTiORWYwbRp8OST8MkncNZZvnf6ySfhsMNghx38DVcaGoKOVESk2yzNmHz3ZHLG57Bu9jo+PflTmho1faeIxJaK51Szww5w553w5Zf+5ipDh/pbgH/ve362jksv1ZAOEemzMgsy2enpncgszqT2mVrmXzSfvnQzMBFJfCqeU1VpKVxxBSxcCPffD/vtB2vWwB//CBMmwAEH+KnuNmreVBHpW3LH5bLj4zuSlp3G8hnLWfSLRUGHJCJJRMVzqsvMhJNPhtdeg/fe80M6cnL8GOnTTvM905de6od7iIj0EYP2GcTkeydDGiy+ZjFLfq8p7EQkNlQ8S4vddvNDOpYvh1tv9a9ra31v9I47wle+An/7m8ZGi0ifUHJ8CZP+7qewW/CjBVTeWhlwRCKSDFQ8y/YGDYILLvA90bNmwXnnQV4evPmmnwJv8GA4/XR49lndwVBEElrZGWWM//N4AOZfNJ/KP6mAFpGeUfEsHdt9d7jtNli2zPdK77cfbNgA99wDhx8Ow4fD5ZfDhx8GHamISJuGnj+UcTeOA2D+9+brLoQi0iMqniU6Awb48dCvvQZffOFvuDJunB/icf31sMsusOuu8Ic/QFVV0NGKiLQy7JJhTLh9Ahgs+PECFl61ULNwiEi3qHiWrhszBq66CubNgzfegPPPh4ICmD0bfvADf5HhEUf43um1a4OOVkQEgPJzy5n0j0n+IsJrFzPv3Hk0bdE80CLSNSqepfvMYJ99/B0Lly+Hhx+GY4/1dzZ85hk/Lrq0FE45BR59VNPeiUjgyk4vY8rDU/w0dn9dzsfH6E6EItI1Kp4lNvr1g69/Hf71Lz8++pZbYN99/fjoBx7w6wYP9kM/nnsOGvWflYgEo+S4EnZ5aZdtN1L5YNoHbKrcFHRYItJHqHiW2Csuhosugtdfh0WL4Le/9eOh16yBv//d3xK8vLxlmyZ9bSoi8TVo70Hs9uZu5IzLYe37a3m34l1Wv7E66LBEpA9Q8Sy9a+RI+NGP4P334bPP4Oqr/R0Mq6v9XNJf/SqMHg0//KGfCk+FtIjESe64XHZ7czfyD8hnc9VmPpj+Acv+sizosEQkwal4lviZNAl+8QuYMwfefdcXzMOGwZIlcMMN/iYsI0bAxRf7OxxqDmkR6WVZxVns/OzODL1kKG6LY96585h7/ly2blT+EZG2qXiW+DODqVPh97+HxYv99HeXXOLnjK6s9OOlDzgAhgyBc87xFx9u3hx01CKSpNIy0xh/43gm/m0i1s9Yfvty3tv7PdbPWx90aCKSgFQ8S7DS0vyNV2680RfSb78NP/6xn0O6uhr++lc/7V1pKZxxhr8gccOGoKMWkSQ05MwhTP3PVLLHZrNu9jpmTZ3FintWBB2WiCQYFc+SOMxgjz3gN7/xc0h/+KEfI73jjrB6Ndx9Nxx/vL8g8YQT/MWHK1cGHbWIJJG83fOoeK+CklNKaFrXxGenf8acs+bQuEYzBImIp+JZEpMZ7LSTHyP90Ucwdy78+tdQUQHr18Mjj/hp78rK/FjpX/3Kb6c7holID2UMzGCH+3ZgwowJpGWnUfX3Kt7Z6R3qXqoLOjQRSQAqnqVvmDABfvITeOcdP/3dLbf4Ke8yM/0sHVdcATvv7O9+ePHFfi7pTZq3VUS6x8woP6ec3d/bnbyKPDYt2cTsA2cz/5L5bF2viwlFUpmKZ+l7Ro70c0Q/8wzU1LT0QpeWti6sNbxDRHqo/+T+7PbGboz6xSgsw6i8qZJZu82i/pX6oEMTkYDEpHg2syNjsR+RLsvL8+Og77zT3yI8vBd67drWwzv23BOuugr+8x/d4VBSmnJ216RlpjHq6lFM/e9UcnfIZcO8DXww/QPmnjOXLXVbgg5PROLMXBfGiJrZ4cCVQBkthXcaMNg5lxP78FqrqKhws2bN6u3DSLJYvBieegqeeAJefLH1dHeDBsHBB/se6sMO8/NLi/QiM3vXOVcR52MqZ8dY06YmlvxmCYt/tRi32ZE5OJPxfxxPycklmFnQ4YlIjHSUs7taPFcDvwbeB5pvBWfAqc65c3saaGeSMRFLnKxb52+88uyz/jFvXuv1kyfD4Yf7Qnr//SGn1+sKSTEBFc/K2b1k3WfrmHfePFa/5m/pXXBoAeNuHEf/yf0DjkxEYiGWxfNbzrm92lhe5pyr6kGMUUnmRCxxtnBhSyH9wgvQ0NCyLjsbpk1r6ZWePNnP/iHSAwEVz8rZvcg1OZbfsZwFP1pAY30jlmEMvXgoo64eRcagjKDDE5Ee6Chnd3XM8wVmdmIby4/pelgiARo9Gs4/Hx59FFatgldegZ/+1N/5cONGX1RfdhlMmeJvIf6tb8Fdd8HSpUFHLtIVytm9yNL8jBx7zt+TIecNwW11LP1/S3lrwlssv3M5rklTZ4oko672PF8MXAvkhS8GnHMuPcaxbSfZezEkQaxYATNn+gL6uee2n6ljwgQ46CD/OOAAKCwMJk7pUwLqeVbOjqOG9xv4/Pufs/p1P5Sj/y79GfvbsRQcWqDx0CJ9TCyHbdQB5wGzgOaJLg04wzl3bU8D7UyqJWJJAM7BJ5/4oR3PP+97qMOHeJj53urmYnq//SA3N7h4JWEFVDwrZ8eZc46V/7eSBT9awKalfq75/IPyGfu7seRNzevk3SKSKGJZPN8LnOWc2xSxvEvj58xsb2A/4HPgdWAd8DugASgAfuycWxP5vlRMxJJgGhv9jVpeeME/3nij9SweWVmwzz4txfQee/gbuUjKC6h4Vs4OyNYNW6m8pZIlv1pCY72fGrP0m6WMvm40OWN0QbJIootl8XwMMBV4MXwx8HXn3CVR7uMsYKxz7sqwZX8DnnHO3R+aWunbzrlvRr43lROxJKj16/280c8/74vp995rfYvwvDx/8WHzEI+ddoI03ZsoFQVUPCtnB2xL7RYW/2oxlTdX4jY7LMMY/O3BjLxyJDmjVESLJKpYFs/vAOOBeqD5jQaUOeeyo3j/NODnwCEudGAzKwcWAgXOufVmlgGsBqY45xaFv1+JWBJeba2fEq+5Z3ru3NbrCwp8MT19un+omE4ZARXPytkJYuPijSy8eiEr/rkCmsAyjLKzyhh5xUiyR3b6oxCROItl8Xwy8G/nXEPE8iOdc09HEwjwBj55j8NfyDIS+L1zbmjYdh+Hlt0V/n4lYulzli71RfSLL8JLL8GXX7Zer2I6ZQRUPCtnJ5j189ez+NrFrLgnVERnGmVnlzHiJyPUEy2SQGI5Vd2XkUk4pCyKICYCuwF/cc59D3gBeAYYDqyK2LwBKO9ibCKJZ9gw+Pa3/TR3ixfDggX+VuLf+hYMHw51dfCvf8Gll8Kuu0JJib/d+B//CLNnQ1NT58cQaZ9ydoLJHZ/L5H9MZs9P96T0tFJco2P57ct5a9xbfHr6p6z9aG3QIYpIJ7paPF8YucDMTgJujOK9U4Ba59yHode3EJoyCdgQsW0/YEvkDqqrq6moqNj2mDFjRpeCFwmUmZ9f+qyz2i+ma2u3L6a//nW46Sb48EMV033MjBkztuUroDiAEJSzE1TuxFx2uHsH9vhkDwafMRiAlfesZNbOs/joax+x+j+rA45QJPVEm7O7OmzjKeC/zrlrzSwXuAk4CZjlnDuok/eeDMxwzuWHLXsfeAB/NfiEsOXzgaudc/eG70NfAUpScw4WLfJjpl9+ue1hHoWFrYd57Lijhnn0EQEN21DO7iM2Lt7Il9d/yfK/Lqdpo/8jedBXBzH8B8MpOroIS9c80SLx1FHO7tL9Q51zR5nZ3mZ2M3AYUAXsgr9YpDMfAoPMbLBzbkVoWSOwDCg3sxzn3AYzywRGAK90JTaRPq+5Z7q5d9o5fxvx5mL65Zd9Mf3oo/4BLcV0c0GtMdMSRjm778gemc34m8cz8ucjWXrTUipvqWT1a6tZ/dpqssdmM+z7wyg7q4yMPN32WyRoHfY8m9n+7azaEzgAuB5oAo5zzv1Ppwczexp4yjn3JzPLxyfb3YG7gbudc0+a2RH4CfxPjXy/ejEkpUUW0y+9tP3twgsKYP/9W3qmd95ZxXSCiEfPs3J28mhc08jyO5ZTeVMlGxdtBCB9YDpDvjuEoRcP1cWFIr2s27NtmNl/gcn4aY7aG2xpwBDnXL8oAikG/oi/29UI/IUon5pZAX7C/UX4C1l+6pzb7qoJJWKRMM75MdOvvOIfbQ3zyM9vKaanTYNddoH0Xr8rs7QhTsWzcnaScVsdNY/VsPTGpax+LfSFQRoUH1tM+YXlFBxYgKVpSIdIrPWkeD4SPzZuZScHONE591DPwuycErFIB8LHTL/yiv938eLW2wwa5Ivp5mEeu+6qYjpO4lQ8K2cnsYZ3G1h641JW/t9KXKP/vztnfA7l55dTdmYZmYW6m6lIrMRsnuegKRGLdNGiRS2F9Msv+9fhBg6Er361ZZjHrrtChsZU9oYgLhgMmnJ279i0fBPL71jO8tuXs2mpv/N6WnYapd8opfzCcgbuMTDgCEX6PhXPIuItXtxSTL/yih/2ES4vr6WYnjYNpk5VMR0jKp4l1poam6h9qpbKWyupe65u2/IBuw6g7OwyBp86mMwi9UaLdIeKZxFp25IlLWOmX34Zvvii9fq8PNhvv5ZhHlOnQqb+M+4OFc/Sm9bPX8+y25dR9bcqGmsbAbAso/i4YoacPYSCgws03Z1IF6h4FpHoLF3aepjH55+3Xj9gAOy7b8swj913VzEdJRXPEg9Nm5qoebyG5Xcs973Rof/i+w3vR9mZZZSdWUbOGM3UIdKZXi2ezWw34AfA+8DDzrlFPdphB5SIReKssrL1MI9581qv79+/pZieNg0qKiArK4hIE16iFM/K2alj45cbqbqriqq/VbFxwcZty/MPyKfsrDKKjy8mY4CGZYm0pbeL59uBY51zZWZ2KvCocy7y1q0xoUQsErBly1oP85g7t/X63FxfTDcP89hjDxXTIQlUPCtnpxjX5Kh/tZ6qO6uofqiapg1+FsO03DSKjytm8KmDKTi0gLRMzQkv0qy3i+dMoNQ5V2lmacAA59yaHu20HUrEIglm+XJ49dWWYR5z5rRen5OzfTHdr9PphZNSAhXPytkprHF1IyvvX0nVXVWseaPlx55ZnEnJySUMPnUwA78yEDONj5bUFtPiOXQ3qULn3D1mtiew1jn3aQzi7JQSsUiCq6pqKaZfeQU+jUgNOTmwzz4twzz22itliumgimflbGnPhoUbWHnfSlbcs4L1n67ftjx7VDalp5Yy+LTB9N+hf4ARigQnZsWzmf0vcAnwjHPuxNCy3wAfOufujUWwHVEiFuljVqzwxXTzMI9PPmm9PjvbF9MHHQSHHOIvQEzSm7YEUTwrZ0s0nHOsnb2WlfeuZOV9K7fNHQ3Qf5f+lH6jlJITS8gdlxtglCLxFcvi+SXgeODbzrk/hpZNA/7pnBsRi2A7okQs0sdVV7ce5vHxx63XFxTAgQfCoYf6Ynr06CCi7BUBFc/K2dIlzeOjV967kuoHq2msb9y2bsCuAyg5sYSSk0rInaBCWpJbRzm7q5fZznTO1ZtZeMV9YvdDE5GUUlICJ5zgHwA1Nb6Ifv55mDnT37Tl4Yf9A2DsWF9EH3ooHHAA5OcHFnofpZwtXWJpRsH0AgqmFzD+5vHUPlNL9UPV1Dxew9oP1rL2g7UsvHIh/Xfq7wvpE0s0tENSTld7nr8BjAGGAq8BJ+F7NS52zv2pVyIMo14MkST3xRe+iJ45E158EerrW9alpcFXvgJf+5p/TJoEfeiipoB6npWzJSaaNjVROzNUSP+rhq2rt25bl7tDbkshvWN/XWwoSSHWFwzuDZwJjACWA/c5557vaZDRUCIWSSGNjfDuu/Dcc76YfvNNv6zZ2LG+iD76aNh//4S/WUuAFwwqZ0tMNW1uou6FOl9IP1pDY13L72XOuByKji2i+NhiBn1lkO5qKH1Wb09Vlw5Md8690KMdRUGJWCSFrVnjC+knnoCnnoJVq1rWDRwIhx/ui+mjjvJjpxNMAk1Vp5wtMdO0pYn6l+q3FdJbarZsW5dZnEnR0UUUHVtE4aGFpOcm58XAkpy6XTyb2ePAgE72XwzUOOcO7H6I0VEiFhEAtm6F//7XF9JPPNF6SryMDH/R4QknwHHHQWlpcHGGiUfxrJwtQWpqbGLNG2uoeayGmsdq2PhFy10N07LTKDikgOJjiyn6WhFZpbp5kiS2nhTP1+EvKpwLNAHn4G/pGp4NR+PnEP1+zCJuhxKxiLRpwQJfRD/+uJ8Wb2toPGZaGnz1q76QPv54GDYssBDjVDwrZ0tCcM6x/tP12wrphrcbWlYaDNxnoC+kjyoid4dcjZOWhNOT4nlQaJv60OtrnHNXtbHdY865Y2MVcHuUiEWkUzU1voh++GE/VnpLy9fI7LVXy2wfY8bENaw4Fc/K2ZKQNi3bxKonVlHzWA11L9ThNrfUHv1G9qPoyCIKjyyk4MACDe+QhBDLeZ5vBy50zm0NWzYdfwHKkJ4G2hklYhHpktWr4cknfSH9zDOwYUPLul13hdNOg9NPh7KyXg8loNk2lLMl4TQ2NFL7bC2rnlxF7b9r2bKy5Q9c62cUHFBA4ZGFFB1VRM6YnAAjlVQWy+L5EOD/Aa8DDcBk4DDgKufcr2MQa4eUiEWk29at8wX0ww/7groh9DVyerq/2PDMM/0Fh710u/CAimflbElorsnR8G4DtU/XsurpVTS80wBhZUnOxByKjiyi6KgiBn11EGlZacEFKykl1lPVlQMXAjsDa4AHnXOP9TjKKCgRi0hMbNwIzz4Ld93lx0o3T4FXVAQXXQQXXwzFxTE9ZIBT1SlnS5+xeeVmap/xhXTds3Wt7nCYPiCd/APyKTi0gMJDC8kZn6Ox0tJrenuqujRgN+fcuz3aURSUiEUk5qqr4d574W9/g9mz/bKcHPjOd+AHP4BRo2JymASaqk45W/qEpsYm1ry5Zluv9LoP17Va329kPwoPLaTg0AIKDiogsyCx53qXviWWwzbubGNxCbDVOXdcN+OLmhKxiPQa5+D11+G3v/XzSANkZcGll8LPfw4DOpsBrmMBDdtQzpaksalyE7Uza6l7ro66mXWt5pQmDfL2yKPw0EIKDyskb8880jI1xEO6L5bF82vAC7QakcS+wGvOuWt7FGUUlIhFJC4+/hh+/WvfIw1+Zo577oG99+72LgMqnpWzJSm5Jsfa99dS+5wvplf/ZzVuS8tpnj4wnYIDC1qGeIzVhYfSNbEsnic65+ZGLJsEnOmc+0nPwuycErGIxNXbb8O55/rhHFlZcPfdcNJJ3dpVQMWzcrakhMa1jax+ZTW1z9VS+2wtG+ZuaLU+e0w2BYcUUHBgAfkH5JNVopu0SMd6e8zzaOAd51xsr65pgxKxiMTd5s1w2WXwpz/5uxe+8ALsv3+Xd5NAY56VsyXpbVy8sWWIx/N1NNY1tlrff+f+vpA+KJ/8/fPJGJgRUKSSqDrK2V06W8xsIa2//ksHBgP/7H54IiIJLCsLbrnFT2H3hz/ABRf4YR194Cp/5WxJVdkjsyn/bjnl3y3HbXU0zGqg7sU66l6oY81/1rDuw3Ws+3AdS29cCukwcI+B5B+YT8FBBQzcZyDpObpRi7Svq8M2ZgD30pKMm4AVzrl5vRDbdtSLISKB2bwZhg+HlSvhk09ghx269PaAhm0oZ4tE2LpxK2veXEP9i/XUvVjHmrfWwNaW9dbPGPSVQRQcVED+gfnk7ZFHWoYuPkw1Met5Bu52zr0asfOJZjbOOfd5tyMUEUl0WVl+7ueVK1tusJL4lLNFIqRnp1NwQAEFBxQw+trRNDY0svq11dS9UEf9i/Ws/WAt9S/VU/9Svd8+L51B+w/yFyAeVED/nfpjaYn/zZP0nq4Wz4cDr0YsWwG8BOwWk4hERBLRgw/Cp59CXh7svHPQ0URLOVukExl5Gf4uhkcT0hjCAAAgAElEQVQWAbC5ZjP1L9f7nukX6tgwbwO1T9VS+1QtAJnFmeRPz9/2yN0hVzdrSTFRFc9m9j3gOiDPzH7cxiYzYxqViEiicM5fLHjZZf71Ndf4m6gkMOVske7LKs6i9MRSSk8sBWDj0o3Uv+QL6foX6tm0dBPVD1VT/VA1AJklmQzaf9C2Yrr/DuqZTnZRFc/OuVvM7FXg28BNEas3O+eWxzwyEZEgOQcvvgg/+5mfsg7g8svhkkuCjSsKytkisZM9LJuyM8ooO6MM5xwbPt/ge6ZDj83LNlPzcA01D9cAvmd60DQV08ks6mEbzrkPzewnzrktkevMbIRzbklsQxMRCcCqVfDAA3DbbfDhh37Z4MF+xo0TTww2ti5QzhaJPTMjd3wuueNzKT+n3BfTX0QU05Wti+mMogzyp7UM8+g/RcV0X9dh8WxmPwSeds59Glr0IzOLnL/FgAOA6bEPT0QkDpYtgyefhMceg5kzYUuo3iwthYsv9rfo7uHtueNBOVskvsyM3HG55I7L9dPihRXTq19ZTd1Ldb6YfqSGmkciiulQQd1/RxXTfU1nPc+jgfyw11OB/kAVLVMfGVAW+9BERHpJVRW88gq8/LJ/zJnTsi4tDQ47DL71LTjhBD+/c9+hnC0SoLaK6Y0LNm7rlW6zmC6M6JlWMZ3wOiyenXMXRSz6X+Bj59zm8IVmFvWEp2aWBbwDXOKce9nMcoDfAQ1AAfBj59yaaPcnItKhLVv8vMyzZsE77/iiee7c1tvk5MDBB8Mxx8DXvuaHafRBvZGzQ9srb4t0g5mRMzaHnLE5DPnOkNbF9Ct+OrxNSzdR82gNNY+GFdP7hxXTmhov4XR1qrqdnHPvhS8ws4nAWODTtt+yncuBUWGvbwWecc7db2aHA7cD3+xiXCIi/kYmc+b4Qvndd/2/s2fDpk2tt8vNhf32g2nTYPp0qKjw8zgnn1jkbFDeFomJNovphS0909uK6X/VUPOvUDFd4Hummy9CHLDzABXTAYt2qrqv4r/qOyR0u9dwhfiruZ+IYj9fAZYDdaHX5cCpQHNvyfPAw2Y2yjm3KJrYRCQFbd0KCxb422R//LHvWf74Y9+j3Ni4/fbjxvkCeffdfdG8++6QmRn/uOMkVjk7tC/lbZFeYmbkjMkhZ0wOQ86OKKabe6a/jCimm4d5HOAf/af01zzTcRZtz3M98A/8eLp9I9ZtBv7a2Q7MrD9wonPuMjO7KrR4OlDjnFsP4JxrDCX6acCiKGMTkWTV0ADz58O8eS2PTz+Fzz6DjRu3394Mxo71xXFzsTx1KuTnb79tcutxzgblbZF4a7OYXhTRM/1l62EemSWhm7YcmE/BAQXkTMhRMd3Lop3n+SMzmw6c5JyLKum24cfAbyKWDQVWRSxrAMq7eQwR6WvWrYOFC+GLL7YvlJd3MB3x8OGw444wZYr/d8cdYfJkPyQjxcUoZ4PytkigzIyc0TnkjM5hyFlhPdMv+YsP61+sZ/PyzVQ/WE31g/6mLVlDsrb1ShccUED2mGwV0zHWlXmeV9NOb4WZ7eece72995rZYcAs59zKyN0CGyKW9QO2m5cUoLq6moqKim2vzz33XM4999wooheRwGzdCpWVfphF82PhwpbnKyPTQph+/fyQi4kTYcIE/5g40RfMgwbF7zN004wZM5gxY0bzy+J4HrsnOTu0TY/ztnK2SGy16pkOjZneMG+DL6Rf8j3Tm5dvZuW9K1l5r//V7Te8X+tiemR2wJ8icUWbs8051946zOxxoLPJTQ2Y4pwr7WA/M4HdwhYVAGuBPwCnOecmhG07H7jaOXdv5H4qKircrFmzOglHROJq9WpYsgS+/LL1v83Ply5tmTe5LVlZMGoUjB7dUiA3P4YPh/TIaYr7JjN71zlX0fmWPTpGTHJ2aF89ztvK2SLx5Zxj/afrW4rpl+tprG19HUj2mOxthXT+Afn0K+9T03HGTUc5u7Oe57n4sXOVtMwRut3+gaM72c+p+J6JZm8ClwFvAZebWY5zboOZZQIjgFc62Z+IxMOmTb7XuK2iuPl5Q0Pn+ykrgzFj/GP06JbnY8ZAebmfW1liIVY5G5S3RfocM6P/lP70n9KfYd8bhmtyrPtoHXUvhorpV+vZuGAjVQuqqLqjCoCcCTkUHFhAwaG+mM7MT96LqWOls+L5t0CDc25TRxuZWYdXbTvnqiO23wpUO+eWmNmTwEHAk8DBwMPOucpOIxeR7tu61Q+XqKz0d9dr71Fd3fm+cnN9D/GIES3/hj8fNkzjkOMnJjkblLdFkoGlGQN2GcCAXQYw/H+G47Y6Gt5v2DbEY/Vrq9kwbwMb5m1g2W3LIA0G7jWQgkMKKDy0kLw980jLVOdGpM5uklIT5X6OJ8qrt9twAfA7M9sFf9crDYgT6a6tW2HVKn+hXUdFcVUVNDV1vr/0dN8zHFkQh/9bWOhnuZDAxSlng/K2SJ9k6cbAioEMrBjIiMtH0LSliYZ3G6h7vo66mXWseWMNa970j8XXLCY9L538A/MpPKSQgkMLyBmnmTygkzHPnb7ZbAR+rs/vOef6xyyqdmj8nKQk56C+3he8VVWwYkXL88jXK1dGVxQDlJTA0KG+OG7vUVqaNGOOgxaPMc9RxKCcLSLtamxopP6Veuqe88X0+jnrW63vN7IfhYcWUnhEIQUHF5CR19V77fUdPRnz3NbO0oFjgHOAQ4AV+DF2ItIVa9e2XwRHPt+8ufP9NSss9GOMy8vbL47LypL1jnoSQTlbRKKVkZdB8dHFFB/tJ5rYuGQjdTPrqJ1ZS93zdWxavInlf1nO8r8sx7KM/Gn5FB1VROFRheSOS53heVEXz2Y2GvgucBZQBCwGvgY8C+zTK9GJ9CXNQyZWruz40VwYr1/f+T6bDRwIgwf7oresrPXz8NelpSqKBVDOFpGeyx6RzZDvDPHT4jU51r6/ltpnaln19CrWvLmGupm+h5pLIWdiDkVHFVF0VBGD9htEWlbyjpXutHg2s5PwPRYHAp/jpym6C/iZc+6Z0GZv9FqEIkFxzvcOd1YMNxfENTX+PdHKzoYhQ9ovipufDx6sC+4kasrZItIbLM3I2z2PvN3zGHnFSDbXbPaF9JOrqHu2jg1zN7B07lKW/mEp6QPTKTy8kJITSig8spCMAck1vCOaT7MBfzvXN4DvOufmAZhZ9wdLiwRlyxY/g0R7BXDksrZuAd2RoiLf+9vZo6wM8vJ0oZ30BuVsEel1WcVZlJ1eRtnpZTQ1NrHmjTWsenIVq55axfpP11P9QDXVD1STlp1G4RG+kC46uoiMQX2/kO70EzjnngSeDF1oco6ZjQIeA7ZNBGhm5npy5aFIdzkHdXXR9Q6vXOm37YqcHN/zG01BXFwMmZofU4KlnC0i8ZaWkUb+/vnk75/P2N+NZcPCDdT8q4bqh6pZ88Yaah6toebRGizLKDy0kJITSyg6pojMgr75f2aXZ9swswzgOOA04GPgDuBQ59yMDt8YA7pyO0Vs2dK6N7j50V7vcGNj5/tslpbmZ5mILHzbK5D79/qEBJIigpptQzlbRIK0qXIT1Y9UU/1QNatfW73t9k2WZRQfU8zgMwZTeERhws0n3VHO7ulUdeOAi4FznHO9PihTibgP27Bh+yK4rcJ4xQqore3avgcNiq5nuLTUz0Shu9lJABJkqjrlbBEJzKaqTb5H+sFq6l+q31ZIZxZnUvrNUsrOLCNval6wQYbEdKq6cM65z4FLzOzdnuxH+iDn/G2ZOyqCw5dFcwvnZs29w829wYMHtzwie4pLSvyFdyLSKeVsEQlSv7J+DD1/KEPPH8rGpRtZec9Kqv5RxfpP11N5cyWVN1eSt0ce5ReWU3pKKek5iXmfgR71PMebejHiYNMmP43asmX+LnXNj+b5hsML4q5cTJeVtX0h3F5xXFSkG3NI0kmEnud4U84Wkc4451j73lqq7qpixT9X0Fjvh2JmFGRQdnYZwy4ZRvbw+HeS9VrPs/QhGza0FMLhhXHk864Mmejfv/NCuPn5oEGaWUJERERaMWuZAm/Mb8aw8v6VLLt1GQ2zGlh6w1Iqb6pk8LcGM+JHI8idkBjTtqp47uuc8zfm+PLLlsfSpf7RXBgvWwarV0e3v4wMP43akCH+UV7u/22ebzi8MNbFdCIiIhIj6bnpDDlrCEPOGsKad9bw5Q1fUv1gNVV3VFF1ZxUlJ5Uw+trRgRfRKp4TmXNQX9+6KA4vkpuXRTN8IjOzdTHc3vPiYl1QJyIiIoEauMdApvzfFNZfu54vf/clVXdVUf1ANTWP1FB+YTmjrhpFZlEwU92peA7a2rWwcCEsWNDyb/PzxYth3brO9zFoEAwf3voxbFhLQVxe7meZ0LAJERER6UNyx+cy8S8TGXn1SBb/cjHL71xO5U2VrPjHCsb8fgxDvjMEi3N9o+K5tznnb9s8Zw7Mm7d9kVxd3fH7BwzYviiOfJ2XGNO6iIiIiPSG7GHZTPzLRIZePJQvfvgFdTPrmHfOPKofqmbiXyeSPSx+FxWqeI6VrVt9UTxnTsvjs8/8vx1dhNevH4waBWPG+Mfo0S3PR47UhXYiIiIiIQN2HsDOz+7MyvtWMv9786l7to53p77LlIemkL9/flxiUPHcHWvXwuzZ8MEH8P77/vHJJ36at7bk5cGkSTBxIowd27pAHjJEY4xFREREomRmDD51MPkH5jPnW3Oom1nH7INmM+H2CQw5e0ivH1/Fc2e2boWPPoL//Mc/3n0X5s/3wzEiDRvmi+RJk2Dy5JbnQ4ao91hEREQkhvqV9WOnp3diwU8WsPSGpcz9zlxck6P8u+W9elwVz5Gcgw8/hKefhpdegv/+d/u742VmwpQpsNtusOuu/t+dd/ZDLEREREQkLtIy0hh3/Tj6De3HF5d9wbxz5pFVlkXx0cW9dkwVzwBbtsCzz8Jjj8G//w2Vla3Xjx4N++7rH3vvDTvs4O+YJyIiIiKBG/4/w9nasJVFVy/is9M/o+KDCnJG5fTKsVK7eP7kE5gxA+67r/WsF0OGwBFHwGGHwX77+aneRERERCRhjbxyJA3vNrDq8VV8/v3P2enxnXrlOKlZPM+aBddeC48/3rJs8mQ49VQ4+mjYZReNURYRERHpQyzNmHDbBN5+6W1WPbGKuhfrKDiwIObHSa1pHlavhvPPhz328IVzdjZccIEvpj/5BK680o9hVuEsIiIi0uf0G9KP4T8cDsDSPy7tlWOkTvH8+eew115w++3+gr8f/hAWLYJbb4Xdd1fBLCIiIpIEys8vxzKNVU+uYnPN5pjvPzWK52XLYPp0mDsXdtrJz8v8+9/D4MFBRyYiIiIiMZRVmkX+9Hxogrpn62K+/+Qvnp2D00/3M2jsu6+fq3nKlKCjEhEREZFeUnCIH+u8+vXVMd938hfP//63n6+5qAgefdTf7U9EREREktaAXQcAsO7jdTHfd/IXz3fc4f+9/HIoKQk2FhERERHpdbmTcwFYP299zPed3MWzc/Dii/75KacEG4uIiIiIxEXWYH8zuy01W3BbXUz3ndzF85IlUF8PpaUwcmTQ0YiIiIhIHKRlppFRmAFNsKV2S2z3HdO9JZrm22yPHq2p6ERERERSSGZhJgCNdY0x3W9yF88rVvh/S0uDjUNERERE4iot25e5TZuaYrvfmO4t0awLXWGpGTZEREREUsq24nmjiufobQmNccnKCjYOEREREYkrFc/dsTl0S8bMzGDjEBEREZG4sn7+eremDX24eDazr5vZQjNbZWZ/NLOM0PJiM7vdzH5rZn8ws9hUuyqeRUS6Le45W0QkhizDF8+uqY9OVWdmI4DjgBOB/wG+C1waWv0gcKdz7sfAPOCamBy0KfSXRkZGTHYnIpIqAsnZIiIxZGmhmda2xna/8awqRwJnO+cagXfNbEfgADN7HZjonHsrtN0TwFwzu845F/t7KoqISDSUs0WkT7P0Pt7z7Jx7LZSEm1UCS4ADgYVh24UmZ6YiBgft8S5ERFJRIDlbRCSWmqvc2A55DvSCwT2A24ChwKqIdQ1AeeQbqqurqaio2PaYMWNGdEfSDVJEJAAzZszYlq+A4qDj6aH45WwRkRhoHrYRbc9ztDk7kMHAZjYOqHbOzTYzB2yI2KQfsN29FEtKSpg1a1b0B1LPs4gE6Nxzz+Xcc88FwMxqAg6n2+KWs0VEYinURey2RlcPRpuz497zbGbpwHnAT0KLKoGCsPUG9AeWxfCgMduViEgqCSRni4jEQPOY52QYtvFD4Hrn3KbQ65nAxLD1o4DNQM+7K9TzLCLSU/HL2SIisdTc89xXLxgEMLOfAe8DOWY2xszOBuqB5aEruQEOB252zm2O4YFjtisRkVQRWM4WEYmBbVPVxbjnOW5jns3sSuDaiMVznHN3mtkpwFVmtggYAFwRr7hERGR7ytki0ud1ccxztOJWPDvnrgOua2fdYuA7vXDQmO9SRCQVBJKzRURiKJnGPMefhm2IiIiIpJSuTlUXreQuntXzLCIiIpKakvAmKfGjnmcRERGRlKKe5+5Qz7OIiIhIauqlCwaTu3hupp5nERERkdTSXP7FuC81NYpnEREREUkp2+Z5VvHcBRq2ISIiIpKammtnjXnuBg3bEBEREUktGrbRDep5FhEREUlJZhq20X3qeRYRERFJLep5FhERERGJUnPHc4xHIqh4FhEREZGko9k2RERERESi1TxsQ7fnFhERERHphIZtiIiIiIhESRcMioiIiIhER1PVdceAATBkCOTlBR2JiIiIiMRR+qB0soZkkd4/Pab7zYjp3hLNRRf5h4iIiIiklFFXjmLUlaNivt/k7nkWEREREYkhFc8iIiIiIlFS8SwiIiIiEiUVzyIiIiIiUUqJ4nnGjBlBh5Bw1CbbU5tsT22yPbVJ71Mbd5/arvvUdj2TSu2n4jlFqU22pzbZntpke2qT3qc27j61Xfep7XomldovJYpnEREREZFYsFjf77s3mVk1sLgbby0GamIcTl+nNtme2mR7apPtdbdNRjrnSmIdTCJTzg6E2q771HY9k2zt127O7lPFs4iIiIhIkDRsQ0REREQkSiqeRURERESipOJZRERERCRKSV08m1mOmd1sZr8ysz+b2cCgYwqCmV1nZi70mB1alnJtY2aHmNlbZjYqbFm77ZAKbdRWm4SWb3fOhJYndZuY2dfNbKGZrTKzP5pZRmh5sZndbma/NbM/mFlm2HvaXSfti/ZcMrMjzOwWM7vVzE6Md5yJqAttl29ma8J+l4+Pd6yJqL28F7GNzrs2RNl2SX/eZQQdQC+7FXjGOXe/mR0O3A58M+CY4srM+uOvgD0ktKj5yveUahszGwwMAPaMWNVROyR1G7XXJh2cM5DEbWJmI4DjgBOBKcCfgS+B64EHgZ84594ys/OBa4Cfht7a0TppX6fnkplNAP4XqHDONZnZq2Y2xzn3cQDxJpJofw/PAb4NNIRevxSn+BJWB/8XhG+j864N0bRdSNKfd0k724aZlQMLgQLn3PpQD9JqYIpzblGgwcWRmV0aenqbc25jaFlKto2ZpQFbgdHOuUUdtQOwub11ydRGkW0SWrbdORNantTnjZl9FXjTOdcYev07/LlwLfCIc648tHwoMBcYDOzU3jrn3Lr4f4q+IdpzycxuA2qdcz8Lvb4CGOecOyuAsBNCF9quH/AkcKFzbn4gwSaotvJexHqdd+2Iou1S4rxL5mEb04Ea59x6gNB/iAuBaUEGFU9mZsA3gBuAFWbW3DMxnRRsG+dcU8Si6bTfDh2tSxqRbdLBOQNJ3ibOudeaC+eQSmAJcCD+czZvVxl6WtHJOmnfdKI7l1q1L/4Pk+lxiC+RTSe6tjsa2AuYZ2ZPmVlxXKNMYG38XxBJ5107omi7lDjvkrl4HgqsiljWAJQHEEsgnLc3/iv4m4F7zOwI1DbNOmqHlGyjDs4ZSL022QO4DZ0nvSHadovcTm0bZds55x4GBuGL6lHA46FeQ+mczrtuSpXzLuk+UBgHbIhY1g/YEkAsgXLO1TnnrgR+DVyK2qZZR+2Q0m3UxjkDKdQmZjYOqHbOzUbnSW+Itt0it1PbduGcC/0x/Cq+J3USvkdQOqfzrgdS4bxL5uK5EiiIWJYHLAsglkRxEzACtU2zjtpBbeQ1nzOQIm1iZunAecBPQotafe7Q0Jb+tHGeRKyT9kV7LkVul3TnWzd0+ffQObcCuJ+W32XpmM67GEjm8y6Zi+eXgHIzywEITR81Angl0KiC1QS8h9qmWUftoDbyms8ZSJ02+SFwvXNuU+j1TGBi2PpR+AtKZ3WyTtoX7bkU2b4TSMIr97uou7+HW4EPejm2ZKHzLnaS8rxL2uLZOVeFv+LzoNCig4GHwy7oSXpmVm5mp5pZWqhH7AfAFanaNqE2ADDo+BxJlTaKbJP2zhlIjd8pM/sZ8D6QY2ZjzOxsoB5YbmY7hjY7HLjZObfZOTervXVxD74P6ehcMrMfmllz4XIbcBhsO1cPxH8bkrKibTszO9HMxoSeT8LPHjE3kKATTGTeCy3TeReFztouVc67pJ2qDsDMCoDfAYuAMuCnzrm1gQYVR2Y2Bfg3sAZ4HbjROTcntC6l2sbM8oDT8fOjXgPc4pyr7qgdkr2N2moToJR2zpnQe5K2TczsSvy0dOHmOOcmm9lI4Cr85x6A/yO0eUq7dtdJ+9o7l8zsXeDXzrmHQtudDkwNve0V59xjQcSbSKJpOzO7Bf/7/SLwMj7ndTZTQtLr4P8CnXediKbtUuW8S+riWUREREQklpJ22IaIiIiISKypeBYRERERiZKKZxERERGRKKl4FhERERGJkopnEREREZEoqXgWEREREYmSimcRERERkSipeBaRuDOzEjN73MwWm9lPgo5HREQ6przdQjdJEZG4M7PDgJeAgcAcYIJzrjbYqEREpD3K2y1UPItIoMzsVeBA3dJaRKRvSPW8rWEbIhIYMxsOPJeqCVhEpK9R3lbxLN1kZj8ysz8FHUezWMRjZgeb2ZtmNj1GYSUVM7s+xvsz4GTgN7Hcr4i0pnzdc2Y2ysz+Yma3x2h/6WZ2jpktjmLbo8zsre5+1ljmbuVtT8VzAjKzyWb2iJn91cxmm5kzs9ODjivC08B9QQcRJhbxfAnsEoNYEpKZHWlmNWa23Mz2MbPhofPMmdktZpYf2m5PM/vYzG4ws5zQssFAVYxDOhW4wznXGNq/SJ+jfN0tfTFfrwz92y9G+0sD6oERUWw7C5jc3kozm2Jm94fOvSvNrCxsXaxzt/I2Kp4T1b+AXzrnvgvsBvyjeYWZ7Wxm04IIKvzYzrmPnXOvBxFHW7oST3tt6JybC9TEMq5Y/rx6ui/n3NPArUCjc+5N59yXwDnAZuAt51x9aLu38cn6cufchtDbvwY82aMPEMbMfg78AnjezOYAE2O1b5E4U77uor6Yr51z64FlsTqWc24L8G6U267AF9rtrf8EmBF6eb1zLrxYjlnuVt5uoeI5wZhZITABWAfgnGsCrgKcmQ3CJ2YLIK7Ajh1LUXyOpjgeK4h93Q8MM7PdAZxzq4AXgKPDjpUPLAude80mOefm9PDY2zjnrnXOjXfOVTjnJjnnXo3VvkXiRfm6dyVgvo71DAtdib+zbbeG/o0chxyz3K283ULFc4IJTfsyF3jGzPYOLVsMvI0fZzQaONfMzgIws2vM7Oehr2y2jWsyswoz+4OZXWpm681srpldZmbHmdk7Zna8mT1qZmvM7Lyw97W5v4hjn2Nmfw8f+xX66vJmM/u1mb1kZnuGlnd4vLD3m5mdZ2afmNkhZvahma0wszOjOMb48Hg6OeZ2bdiGgtDXsA1mdk5Em/4i1DaPhJJttG223bHM7EQze9vMzjSz+Wb2mw7219bPvs14OhLqofgEOClscQlwlIWGaADHAI+HxZkDrA97PcnMfmNm15vZK+aHf5xiZrPM7Fgzu8/M6s3sEjPbPfQZV5nZkZ3FJ9KXKF+nTr5uS3ufsaPjmNloM7vTfC/u76M5ThvHfdrM7jGzfTvZTrm7tzjn9EiwBzAFWIz/S/NOoDRs3SJgeuj5jsCC0PMS/F/FRfi/nJcAFaF1fwLeCD3PAmqBW/BzNX4TqOxof5HHxo/5uh34e2h5Nr4gyw+9PhGoBvI6Ol4bn7ssdMyfhba9DtgIjOzkGJHxdHjM8DZsI4ZFwA1Abuh9S0LLBwL/F7bdf/CJL6o2a+dYhfiv4m4E9gQO7mh/ET/7NuOJ8vy6Evgi9Hws8GtgE/D1sPPFwrY/Bvhq6Hku8A6QG3r9aqgdBgKrgWuAHOCU0OuT8X+k/4DQOaiHHsn0QPk6JfJ1aP0vomzHjn7W7wCTw97jojzPmn+e+wBXt7F+eug4GWHLlLt76aGe5wTkfO/gTvhk8i3gEzNr68KI+cAZZpYJfDW0bAD+l3U40Dxe9TP8LwXOuc1AA/CQc24N/peprJP9Rca3idYXIBwFbHChMbPAI/hfuq93crzI/Tbv897Qttfhx+Me0skxWsXTlWO24wnnx7e9AwwJ+4xFoZ6hS/FtClG2WTuftxafpB5zzr3tnHu+C/trL55o3A+MMT904wTgL8CzwMlmNgBocKHsGrIv8Ebo+TH4/6CaezOOB64NtXM98KLz46RnAQOdcw84/1X2+3TtZyDSJyhfp0a+bkO7n7GD4xwMDHLONcczq4vH3AM4zTn3yyi3V+7uJRlBByCtmVk6MNw5twj4vpndgb8yeQawV/i2zrlNZjYU+DFwV/MunHMrzWwWcBD+L+PR+IJp21vDnjcSGr7T3v7aCTV8H+OBzLC4msxsETCso+N1xjm30cwWAPlAcReO0e1jRmik5XdkBLDQOXdj5EbRtJmZnYEvUptNdP7rXRceaxd+Bu3G0xnn3Hwzex8/dGOQc26Bmd2P7w06CXgqLO7mc6N5PN1IfE9R875Whe867PlWWmtCwyYM3ukAAAPkSURBVMQkyShft/p8qZCvw7Xbjh38bHbG/6HQXQcC481sgHNubTvbNIU+g3J3L1KDJJ5MwsajOudm478WmxS5oZkdCFzonLvO+ZkTwl0E7GF+yqQPgd92duBO9teRJfiezPApfAw/FrCnsoB5vXyMaCyn9bhgzGyvLrTZ48CuYY82r9ruwv7ajCf6j8P/AWcBC8LiSwcuwX/F2Wwv/PjNZpXAfmbWP+y4HY67E0liytetpVK+bvczdnCcNaH3ZNE91+N77f/YzvrVruVCb+XuXqTiOTFdbGblYa+HAs+Enm/GXyAxCT8t0iAz62dmB4fWl5pZMX7s3R34r1zeIPQ1YEgaEX9tm5l1sr/IY1vYPv4FrAW+EdpXfuh4zT2Y7R2vPUNC2wwNvfeZKI4RHk9nxwz/HJHa6oUwfG9Sf+BJMzvM/IUpnf0Mth0LGOKcmxP22BJ2iPSw5x3tLzzuNuMxsywzmxEaktGRB/BfFz8E4JxrAP4N/Me1nmXjUPyQjmbN7X2f+bmivwc0z/UZ+TOI/Dn36Sv/RdqhfE3K5Oto27HN4+BzaRb+Dyzw3zJgZuVR5G7D9wKfAZxkZieH3vsDM7sS3+v+Qtj2yt29qTcGUuvR/Qf+IgQH1AF3A3/Hf+1TEFp/NbAUOBb/Ncx8/F/6R+F7LP6OT0SP4y9e2Bza33r8L9MRwBb8RSnDgZ+G1p/d0f4ijn0J8Bp+HNmU0LoKfI/lDfixfzuGlrd7vHY+vwu9/xL8MIKdwta1d4xR4fF0dszwNow49rH4C+duCbVF8/vOCq2fBnwc+tncEGrnaNvs2DY+64mh490HDA0t6+hn2mpf7cQzDD+Z/71RnGs3Rbw+BTggYtn1bbxvf+BTYBVwZWjZkfhz7WagnNB0XcD5+P9c/xL6rEcH/Tumhx6xeqB8nUr5egTwOr73fNdOPmNHP+uj8N/4/Rd/8fY84MLQZ28zd+O/3diC73EuBF7EF+7fBc7Dz3d9A37ISPN7lLt78WGhxpMkYmalwPnOuWtCr9PwvSHnOeeuDDS4TpiZA0Y7P4ZQusnMLnfOdWsaJBGJH+VrCafc3Tdo2EZy+h7+K/zCsGUHEnYhWCIK+5pIXxF1k5mlm9lpwMNBxyIiUVG+FuXuPkbFc3K6DT826jMz+xKYCSx1zr0ZbFjtC11AcUno5Znmp0yTrhsA/Ms5t6DTLUUkEShfCyh39ykatiEiIiIiEiX1PIuIiIiIREnFs4iIiIhIlFQ8i4iIiIhEScWziIiIiEiUVDyLiIiIiERJxbOIiIiISJRUPIuIiIiIREnFs4iIiIhIlFQ8i4iIiIhE6f8DVfGaf/t1trEAAAAASUVORK5CYII=\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/paet-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/paet-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/paet-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 }