{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 35 - Hayabusa - 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 Hayabusa sample return aeroshell.\n", "https://en.wikipedia.org/wiki/Hayabusa" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \n", "planet.h_skip = 120.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('Hayabusa', 16.27, 114.0, 0.00, 0.128, 0.0, 0.202, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,11.3,0.0,-13.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 (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+ZLJN9IwtknSyQhABJcEjwQiRA2IwsYhBB1AgYUeGi98pFcUEBFRH8uWDQEfGiCAZE5SYgSJAliCzBEAiQkEA2Ekgm+77MzPn98VRnejo9S890T3XPfN+vV72qu6q66nTNzOkzTz/1lLk7IiIiIiLSNEVxByAiIiIi0hKosBYRERERyQIV1iIiIiIiWaDCWkREREQkC1RYi4iIiIhkgQprEREREZEsaBt3ANnSu3dvLykpafD25eXl9OnTJ3cBNQO9h/zQEt4DtIz3Uajv4eWXX97g7oUXeBNkmrMTCvVnnI90LrNH5zJ7CuFc1pWzW0xhXVJSwvz58xu8fWlpaUbb5yO9h/zQEt4DtIz3UajvwcxWxh1Dc8s0ZycU6s84H+lcZo/OZfYUwrmsK2erK4iIiIiISBaosBYRERERyYJWW1jPmDEj7hCaTO8hP7SE9wAt4320hPcgddPPOHt0LrNH5zJ7Cv1cmrs37wHNTgNuAi509xXRsvOB24BuwD3Af7t7RbSuN/A9YAvQDrjW3fen7re0tNTzvU+OiEg6Zvayu5fGHUc6ytkiIjXVlbObtcXazPoBXYAJScsGA+cB04CvAJcDX0562QPAXe5+LfAWcEOzBSwi0oopZ4uIZKZZC2t3Xwc8lLJ4CHCpu7/s7r8DfgGcDGBmxwMj3f2FaNvZwFVm1rm5YhYRqY+7U1lZSUVFRdyhZFU+5Wx3qKgIk4hIvmr2PtbuXpXyfF7iK8TIGmBV9PgUYHnStmuih3n5lamItE6PPfYYbdu25SMf+UjcoWRdvuTsf/4T2rWDk05q6p5ERHInHy9ePA74ZfR4ALAxZf12oH/qi8rLyyktLT0wlZWV5ThMEZGgqCik0srKyga/pqys7EC+AnrnJrJm0Sw52yzMm/myIBERoOE5O69uEGNmI4Byd18YLXJgd8pmxcBBF8L06dMn7wcUF5GWqU2bNgBUVVXVs2W1GTNmHLj63cw25CSwHGvOnK3CWkTi1NCcnTct1mbWBvg88LWkxWuAnknbGNAZWNu80YmI1K4xLdaFrrlztgprESkEeVNYA18FbnX3vUnLHgdGJj0vAfYBapoWkbyRaLFuTYU1zZyzVViLSCFo9sI6asEAsKRl1wELgI5mNszMLjWzEe4+H3jPzI6KNj0T+Lm772veqEVEateYriCFIl9ytgprESkEzdrH2sy6ApdET6eb2e2ErxJvTNl0sbvfFT2+EPi2ma0gjKf6jeaIVUSkoVpqV5B8ytkqrEWkEDRrYe3u24E7oinhpmiq7TUrgctyHJqISKO11K4g+ZSzVViLSCHIpz7WIiIFqSV3BckXKqxFpBCosBYRaaKW2hUkn6iwFpFCoMJaRKSJWmpXkHyiwlpECoEKaxGRJlJXkNxTYS0ihUCFtYhIE6krSO6psBaRQqDCWkSkidRinXsqrEWkEKiwFhFpIrVY554KaxEpBCqsRUSaSBcv5p4KaxEpBCqsRUSaSF1Bck+FtYgUAhXWIiJNpK4guafCWkQKgQprEZEmUleQ3FNhLSKFQIW1iEgTqStI7qmwFpFCoMJaRKSJ1BUk91RYi0ghUGEtItJE6gqSeyqsRaQQqLAWEWkidQXJPRXWIlIIVFiLiDSRuoLkngprESkEKqxFRJpILda5p8JaRAqBCmsRkSZSi3XuqbAWkUKgwlpEpIl08WLuqbAWkUKgwlpEpInUFST3VFiLSCFQYS0i0kSWqPpQcZ0rKqxFpBA0e2FtZqeZ2QtmVpK0rKOZ/dzMvm9md5hZt4asExHJFy21O0i+5GwV1iJSCJq1sDazfkAXYELKqpnAs+5+HfAQ8KsGrhMRyQstsTtIPuVsFdYiUgiatbB293WERHuAmfUHLgZmR4vmAueYWUld65olYBGRBmqJI4PkU85WYS0ihaBtcx/Q3auS+yMCk4EN7r4rWl9hZsuBk4D9daxb0Zxxi4jUpaV2BcmXnJ0IoQV9ISAiLVCzF9ZpDAA2pizbDvQHKupYJyKSN1piV5BaxJKz27cP8/37m7onEZHcyYfC2oHdKcuKCS0f1LPugPLyckpLSw88nzFjBjNmzMhimCIitcu0K0hZWRllZWWJp71zE1VOxJKzE4X13r2Zhisi0nQNzdn5UFivAXqmLOsKrCUk8NrW1dCnTx/mz5+fkwBFROqTaYt1ciFpZhtyFlj2xZKzi4vDfN++hgcqIpItDc3Z+TCO9ZNAfzPrCGBm7YDBwNP1rBMRyRst8eLFWsSSsxOF9d69uoBRRPJXHONYJ66CMQB3fx+YA5waLZ8CPOjua+pa14whi4jUq6VevJgvObuoCNpG37Gqn7WI5Ktm7QpiZl2BS6Kn083sdncvB74A3GJmY4FDgeSOdnWtExHJCy3x4sV8y9nFxVBREVqtE32uRUTySbMW1u6+HbgjmpKXbwY+V8tral0nIpIvWmJXkHzL2cXFsHNnKKy7ds3FEUREmiYf+liLiBS8ltoVJJ9oZBARyXcqrEVEsqAldgXJN8kXMIqI5CMV1iIiWdASu4Lkm86dw3znznjjEBGpjQprEZEsUFeQ3OvRI8y3bo03DhGR2qiwFhHJAnUFyb3u3cN8y5Z44xARqY0KaxGRLFBXkNxLtFirsBaRfKXCWkQkC9RinXsqrEUk36mwFhHJArVY5576WItIvlNhLSKSBbp4Mfd69gzzDRvijUNEpDYqrEVEskBdQXKvf/8wX7s23jhERGqjwlpEJIvcPe4QWqxEYb1mTbxxiIjURoW1iEgWJLqAJFquJfsGDAhzFdYikq9UWIuIZIEK69xLFNZr14K+GBCRfKTCWkQkCxJ9q1VY507nzuECxn37YN26uKMRETmYCmsRkSxItFgnht2T3DjiiDBfsiTeOERE0tEngIhIFqgrSPMYPTrMFy+ONw4RkXSyUlib2YezsR8RkUJVSIV1IefsUaPC/M03441DRCSdtplsbGZnAt8EDqW6KC8C+gEdsxuaiEjhyMfCuiXm7DFjwvyVV+KNQ0QknYwKa+D3wA+ABUDiLggGXJzNoERECk3i4sU862Pd4nL2hAlh/tJLUFEBbTP9FBMRyaFMU9I77v7j1IVmpt5uItKq5WOLNS0wZ/ftC8OGwTvvwOuvw9ixcUckIlIt06aVL5jZtDTLz8lGMCIihSpPC+sWmbOPPz7Mn3km3jhERFJlWlifANxpZpVJUxVwRw5iExEpGHv27AGguLg45khqaJE5+7TTwvxvf4s3DhGRVJl2BbkBmAHMByqjZQZ8qqmBmNlo4EpgGXA4UObur5hZR+AWYDvQE7jW3bc19XgiItm0a9cuADp37hxzJDW0yJx91llh/uSTsGsXdOqUzb2LiDRepoX134CH3H1v8kIz+3UWYvk9cK67rzGzwcBjwGhgJvCou8+KrnD/FXBRFo4nIpIVVVVVBwrrTvlV5bXInN2vHxx3XLiA8ZFHYFq6zi4iIjHItCvIH4Gvm9mHkqaTgK9nIZbRQNfo8W6gu5n1J1y9PjtaPhc4x8xKsnA8EZGs2L17NwAdOnTIt1FBWmzOvigq1X//+2zvWUSk8TJtsf4W4Su/6YBHy4wwRurVTYzlT4S+gGcBlwBXAZOBDe6+C8DdK8xsOXASsKKJxxMRyYqdO3cCedcNBFpwzr74YrjmmtBi/d57cNhh2dy7iEjjZNq08iNgkLuXuPvQaCoBzs9CLF8E9gAvATvc/UFgALAxZbvtQP/UF5eXl1NaWnpgKisry0JIIiL1a0w3kLKysgP5Cuidm8habs7u1w/OPTeMZf2znzXuDYiINFRDc3amLdar3X17muWHZrifdDoC9xES80+iVg4nfMWYrBjYn/riPn36MH/+/CyEISKSmca0WM+YMYMZM2YAYGYbchJYC8/Z//M/8Oc/w8yZ8LWvQffuTdqdiEitGpqzM22x/mLqAjO7APhJhvtJ5x7gQXe/AbiZkLDfI1xVnqwrsDYLxxMRyYrNmzcD0KNHj5gjOUiLztkTJ8LkybBtG9x8cy6OICKSmUwL615m9i0AM+tkZncCdxK+Cmw0M+sNjHX3LdGim4BtwEqgfzR8E2bWDhgMPN2U44mIZNP69euB0AqbZ1p8zk4U1D/+cbgbo4hInDIqrN19KvC4mf0ceAU4AhgLNHWwo03AHjMbkLRsI/AqMAc4NVo2hdBCsqaJxxMRyZry8nIA+vbtG3MkNbWGnD1xIlxyCezbB5dfDlVVuTqSiEj96uxjbWYfqmXVSmApcCswBDgP+Epjg3D3KjM7D/i2mb0M9AOucfdtZvYF4BYzG0voFzijsccREcmFRGEdd4t1a83Zt90Gjz0Wbhjzs5/Bl7+c6yOKiKRX38WLtxDGKt0CpGsHuIswdNNhNCFJA7j7QuDzaZZvBj7XlH2LiORSvhTWtNKc3bcv/PrXcN554YLG0lI48cTmjEBEJKivsL4BmO/u6+vayMx03ysRabVWrVoFwKBBg2KOpPXm7HPPhauvhp/+FM4/P9yVcciQuKMSkdamzj7W7v5IfQk62u5P2QtJRKSwvBNdNTds2LBY42jtOfvWW2HKFCgvD/P33os7IhFpbfLq3rsiIoXG3fOmsG7t2raFBx6AceNg2bJQXK+v998MEZHsUWEtItIEGzZsYMeOHfTo0YOePVOHcJbm1qMHPP44jBkDb7wBJ5wAb78dd1Qi0lqosBYRaYLXXnsNgFGjRsUciST07g1PPAHjx4eW6w9+EF54Ie6oRKQ1aHJhbWbjzeweM/tvMytpekgiIoXj3//+NwDjx4+POZKGaS05u18/ePppOP300Of6Qx8Ktz53jzsyEWnJstFifQUwxd1vA/4jccctEZHWYMGCBUDhFNa0opzdtSvMmQNf/GK4gcyXvgQXXhgKbRGRXMhGYX0l8IHo8R+BdlnYp4hIQXjxxReBgiqsW1XObtcOfvELuO8+6NIlXNw4Zgz8qUWOiyIiccu4sDazs8zsk9HjCcDhidvVunuVu2/LcowiInlp1apVLFu2jK5duzJu3Li4w0lLOTv4xCdg4UI4+eTQYn3BBWG865Ur445MRFqSjAprM/se8ADwUQB3fxH4tJldnIPYRETy2hNPPAHAySefTNu29d1vq/kpZ9c0bBjMnQt33BFar//yFxg1Cr7zHdi1K+7oRKQlyLTF+j+AgcC8pGV/A27OWkQiIgXiscceA2DKlCkxR1Ir5ewURUVwxRXw5ptw0UWwZw9897swejT87/9CRUXcEYpIIcu0sH7c3bcAyddVt7hb44qI1GfXrl3MmTMHgKlTp8YcTa2Us2sxcCDcey8880y4ocyqVfDZz4b+1/fdB1VVcUcoIoUo08L6HTO7DhhpZp8wsweBLwE/zH5oIiL56+GHH2bnzp1MmDAhn++4qJxdj0mTYP58+P3vYfhweOstuPhiOOYYuPvuMJqIiEhDZVRYu/sfgX8AbYBPA1uA0939FzmITUQkb919990AXHjhhTFHUjvl7IZp0wYuuSR0D7nzThg8GF5/HaZPh5IS+MEPYNOmuKMUkUJg3sTR8s2sDTDZ3Z/ITkiNU1pa6vPnz48zBBFpJZYuXcoRRxxBcXExq1evpk+fPk3an5m97O6lWQqvvmMpZ9dj797QHeS222DRorCsUyf45Cfhc5+D0lIwizdGEYlPXTm7zsvYzez/gC717L83sAGINUmLiDSXn//85wB88pOfbHJRnU3K2dlRXBxaqz/zGXj88VBg//3v8Otfh2nsWLj88lBo9+wZd7Qikk/qGx/q1WibJUAV8DlgAZDczDAU6JWT6ERE8szatWu58847Abj66qtjjuYgytlZZBZuiX766dXdRO6+O4yHfdVVcM01cM45YYzss86CDh3ijlhE4lZfYf0jQneRLQBmNtzdv526kZk9lIvgRETyzXe/+112797Nxz72MY455pi4w0mlnJ0jo0eHluvvfx8eeigU2Y8/DvffH6auXeG880KRPWUKtG8fd8QiEoc6L150962JBB3pF/XPO8DMJgMTchCbiEheWbhwIb/5zW8oKiripptuijucgyhn515xMXz846FryIoVcMstcOyxsH17GFlk6lQ49FD41KdCwb11a9wRi0hzyvRWYX8CFprZs8B2YDRwBnBQi4iISEtSUVHBZZddRmVlJVdeeSWjRo2KO6SGUM7OoSFDQneQa66BpUth1qxw0eMbb8A994SpbVs46SQ4++ww5e/IjCKSDRmPCmJm/YEvAscA24AH3D32rxXz+QpzESl8N998M1//+tcZNGgQr7/+Ol27ds3avnM5KohydvNbvBhmzw7TP/9Z82Yzo0aFriJTpsDkydC9e2xhikgj1ZWzszHcXhEw3t1fbtKOmqglJ2kRidfTTz/NqaeeSmVlJY888ghnnXVWVvffzMPtKWc3o40b4W9/C0X2o4/Ctm3V69q0gQkTqgvt449X32yRQtDo4fbS7OiuNIv7AJXAeY2ILd0xjgdOBJYBzwI7gVsIX2P2BK51922170FEJHvWrl3LhRdeSGVlJddee23Wi+pcUs6O3yGHhJvPXHIJ7N8PL74Ic+eGCx+ffx7+9a8w3XgjdOwYiutJk8L0wQ9C585xvwMRyUSmfawPJ4x9mtzMPRCYl41gzOyzwHB3/2bSst8Cj7r7LDM7E/gVcFE2jiciUpctW7Zw5plnsm7dOk4++eS8vGCxHsrZeaRdOzjhhDBdf31ovX7mmepC+4034MknwwShRfvYY6sL7RNPhN69430PIlK3jLqCmNlId1+SsmwUMN3dv9akQMxOAr4FnOZRUFHfwOVAT3ffZWZtga3AGHdfkfz61vK1oog0j507d3LWWWcxb948Ro0axbx58+ido6omV11BlLMLy/r18OyzMG9emC9YAJWVNbcZPhwmTqyexo0LI5WISPPJdR/rocBL7t6kTxwzexl4DjBgBHAjMAT4kbsPSNpuUbTs7uTXK0mLSLZs2bKFqVOn8txzzzFgwACee+45Bg8enLPjNXMfa+XsArF9e+guMm9emF54AXbvrrlN+/ahuJ4wIRTaEybAiBFQVOdguiLSFNnsY72cml8ptgH6Ab9vfHihVQUYD3zW3V81s2uAR4GbgI0pm28H+qfuo7y8nNLS6vc4Y8YMZsyY0ZSwRKQVWrt2LR/5yEdYsGABgwYN4oknnshJUV1WVkZZWVniaU6awpWzC1vXrnDaaWECqKiARYtCgZ2Y3nwz9Nt+8UW4/fbq140dC+PHh+nYY+HII0NXFBFpnIbm7Ey7gpQB91KdqKuAde7+ViPjTOz3fKAs0YJiZh2BcuAG4GPuPjFp238D97r7rcn7UOuHiDTV/PnzOffcc1m7di0jRoxg7ty5DBkyJOfHzWFXEOXsFm7bNpg/v7rQnj8f1qw5eLv27eGoo6qL7fHjQ/GtiyNFMpe1FmvgHnd/JmXnI81shLsva3SEIY4Dsbj7bjNbSmhd6ZmybVdgbROOJSJSg7vz29/+li996Uvs2bOHSZMm8eCDD9KnT5+4Q2sq5ewWrls3OOWUMCWsXx/6ZydPS5fCv/8dpgQzOPxwOProUHQn5iNGhAsnRSRzmRbWZwLPpCxbBzxJ+FqwsV4FuptZP3dfFy2rICTj/mbWMUrc7YDBwNNNOJaIyAFbtmzh85//PPfffz8Al112GTNnzqR9yxhQWDm7FerbF844I0wJ27bBwoWhyP73v8P8jTfgrbfC9OCD1dt26BC6jiQX20cfDf37h2JcRGrXoMLazK4k9J3rambXptnk8aYE4e6LzexvwDTgF2bWA2gP/AE4CzgVmANMAR509zRfdImIZObhhx/mC1/4AqtXr6ZLly7MnDmTT33qU3GH1WTK2ZKqW7fqYfsS9uwJd4lctAhee616vnr1wa3bAD17hiL7yCPDHSRHjw7zQYN0saRIQoMKa3e/3cyeAT4D/Cxl9T53fy8LsXwa+KmZtSe0cFzk7hVm9gXgFjMbCxwK6OoWEWmSdevWcfXVVzNr1iwAjjvuOO69915GjBgRc2TZoZwtDdGhQxhRZNy4msu3bg1FdmrBvWlT9QglyTp1gpEjQ5GdXHAffng4hkhrkunFi+3cfX+a5YPdfVVWI8uQLoQRkfrs2bOHn/zkJ3z/+99n+/btdOrUiRtvvJH//M//pG3bTHvGZU8OL15UzpascIf33w8F9uLFYTSSxYvD9P776V9TVARDh1YX3EccEYrtESNgwAC1ckvhavTFi2b2VeARd38jWvQ/ZpZ6SYMBJwOTmxqoiEguVFZWMmvWLK677jpWrlwJwNSpU7n99tspKSmJN7gsUs6WXDGDww4L0+mn11y3eTMsWXJwwf3229XTww/XfE2HDuFmNyNGVBfbifnAgSq6pXDV10QzFOiR9PxYoDPwPtXDNxnh6z4RkbxSWVnJ/fffz4033sibb74JwNFHH81tt93GaYnBgVsW5Wxpdj17wvHHhynZ3r2hqE4U3MuWhdFJli2Ddevg9dfDlKq4+OCie/jw0Po9eHAYOlAkX9VZWLv7l1IWfQ9Y5O77khea2ZHZDkxEpLH27t3Lfffdxw9/+EMWL14MwODBg/n2t7/N9OnTadNCxxJTzpZ8UlwcLnQ8Ms1v27ZtoehOLrYT8/ffDyOWvPHGwa8rKgrdSIYOTT/176/WbolXpp0Kj3b3GtcJR3fgGg6k+RMQEWk+69ev54477mDmzJmsX78egJKSEr7xjW/w6U9/uqUMoZcJ5WzJS926Vd+oJtX27TWL7qVLYflyeOcdePfdMGrJ6tXwTOpAkoTW7CFDQpE9bNjBhXevXhoyUHKrocPtTSJ8fXhadIvcZL0IV53PznJsIiL1cneefvppfvOb33D//fezb19onB07dixf+cpXuPjii2nXyu7lrJwthaxr1/SjlQDs3w+rVoVCO920fn11MV7bvktKwjRkSPWUeN6njwpvaZqGtlhvAX5H6L93Qsq6fcCd2QxKRKQ+a9eu5e677+auu+5i2bJwE0Ez45xzzuHLX/4ykydPxlrvJ6RytrRI7dqF/tbDh6dfv2MHrFhRe+G9fXsY2eS119K/vmPH0I87udhOLr4PO0x3pZS6NXQc69fMbDJwgbsrIYtILLZt28b//d//8cc//pFHH32UyspKAAYMGMD06dO59NJLGTZsWMxRxk85W1qrLl3CTWyOOurgde5hLO4VK2DlyjAlP165snqEkyVL0u+/bdtwQ5zUlu7ENGiQLq5s7Rrcx9rdt1JLK4eZnejuz2YtKhGRyK5du5gzZw6zZs3i4YcfZu/evQC0bduW888/n8suu4wzzjijxV6Q2FjK2SI1mcEhh4TpAx9Iv822bTUL7dTCe9266tbv2o7Rv3/6biaJwrtLl1y9Q8kH9Y1j/X9Afb8CBowB+mYrKBFp3TZu3MgjjzzC7NmzeeSRR9i5cycQunpMmjSJT3ziE0ybNo2+fZV2kilnizRNt25w9NFhSmf37tDHu7bie82a6um559Lvo0ePMFb3wIGh0E48Tn7etWvO3qLkWH0t1ksIffXWUD0GaioDPpLNoESkdXF3Fi9ezOzZs5k9ezbPPfccVVVVB9Yff/zxXHjhhVxwwQUMGDAgxkjznnK2SA517Bhu3z5yZPr1+/eHkUuSC+/k4vvdd2HLljAtWlT7cbp1q7vwHjgwbNN6LyPJX/UV1j8Etrv73ro2MjNdXS4iGdm8eTNPPvkkc+fO5e9//ztvv/32gXXt2rXj1FNP5eyzz+acc85hyJAhMUZaUJSzRWLUrl310H7puMPGjWG4wHffrZ6Sn69eHbqk1DaWd0KXLukL78SyQYOge/fcvE+pXX03iNnQwP18FF1lLiJ12Lt3L8899xxz585l7ty5zJ8/v0ar9CGHHMLUqVM5++yzOf300+nWrVuM0RYm5WyR/GYGvXuHKd0Y3lB9kWVthXfi+Y4d4Y6W0U1l0+ratWahnVp4q8939mV6g5gazGww8CXgSpSkRSTJjh07eP7555k3bx7z5s3j+eefZ/fu3QfWt2vXjkmTJjFlyhROPfVUJkyYoAsQc0w5WyT/JV9kOXZs+m3cQ3eS5EI7UXwnz7dvr7/lO9Hnu7bCe+BA6NQpN++1Jcq4sDazNsA5wOeA04B1hD59ItKKlZeX8+yzzx4opBcsWHBgOLyEY445hilTpnDaaacxadIkOnfuHFO0rYdytkjLYwY9e4aptgst3cPwgcmFdroCvCF9vnv1OrjoHjy4+mY7Gt+7WoMLazMbClwOfBY4BFgJnA08BnwwJ9GJSF7avXs3CxYs4IUXXuDFF1/khRdeYHnK+FNt2rShtLSUSZMmMWnSJE488UT69OkTU8Stj3K2SOtmFgriXr3qbvlO9PlOLbqTu59s2hSmhQvT76ddu5qFdmKYwcTUv3/rKbzrLazN7AJCS8cpwDLgx8DdwHXu/mi0WS2DyohIoauqqmLx4sUHCugXX3yRV199lYqKihrbdezYkYkTJx4opD/4wQ/SRZ33mp1ytog0VEP6fFdVQXn5wUX3qlVhtJMVK8L43m+/HaZ02ratWXinFt8DBrScwrshLda7CbfAfQ643N3fAjCz2oZyEpECtXfvXhYtWsSCBQt45ZVXWLBgAQsXLjwwjnRCUVERRx99NBMnTmTChAlMmDCBMWPG0LZtky7bkOxQzhaRrCkqgn79wlRamn6bxPjeiUI7dXr/fXjnnTCl064dDBsGI0YcPA0ZEtYXino/Bd19DjAnuujlc2ZWAjwEHHibZmburqQtUkC2bNnCK6+8cqCAXrBgAW+++eZBLdEAgwYNOlBET5w4kWOPPVat0XlKOVtEmlt943vv2VN34f3ee7XfSr5t29CqnVxsjxwJRx4Z+nrn21jemdzSfBXwLTNrC5wHDDSzG4HfAKcDZbkJUUSaYufOnbzxxhssWhSecEIAACAASURBVLSIRYsW8frrr7No0SLWrFlz0LZFRUWMHj2acePGMX78eMaPH8+4cePo3bt3DJFLUyhni0i+6NABjjgiTOns3Blas5ctO3havbr6caouXWD06FBkjxkT5kceGVq5i4py+55qY01ptDCzEcBVwOfcPdbBWEpLS33+/PlxhiASqz179rBkyZIaxfOiRYsOuqgwoWPHjowZM+ZAAT1+/HiOPvpojdQRAzN72d1r+ZI1q8dRzhaRgrJnDyxfHgrrpUvDtHhxGEJw/fr0r+naFT7wgdB15bjjwnzo0Oy1bteVs5vUIdLdlwFXm9nLTdmPiDSMu7N+/XqWLFnC4sWLWbJkyYHHy5cvr3HDlYR27doxcuRIjjrqKI466ijGjBnDUUcdxdChQzVudCujnC0ihaZDh9AqPXr0wes2bAg3yEmM1f3GG/D666FryVNPhSmhXz847TQ4/XSYOjWMlpILTWqxzjYzaw+8BFzt7k+ZWUfgFmA70BO41t23pXutWj+kJdm7dy/Lli1LW0Bv3bo17WuKiooYMWLEgcI5MR1++OG0K6QrP1qh5mqxzoXG5m3lbBHJlfffh5dfhpdegvnzwzy5dbu4GC68EK6/Plw0mamctVjnwDVASdLzmcCj7j7LzM4EfgVcFEdgItm2d+9eli9fzrJly1i2bBlvv/02y5YtY+nSpbW2PgN0796dkSNHMmrUKEaOHHng8fDhw+nQoUMzvwsR5W0RyS+HHhpapadODc/dQ/eRv/8d5syBJ56A3/0OHnoIFiwI3USyJW8KazP7D+A9YHP0vD9wMeH2uwBzgQfNrMTdV8QSpEiGduzYwdtvv32gaE4uoFevXk1t3xgVFRUxfPjwGoVzYt63b18s3y6DllZJeVtECoFZdXeSq6+Gp5+GyZNh69ZQZF9+efaOlReFtZl1Bqa5+3+Z2bejxZOBDe6+C8DdK8xsOXASsCKWQEVSJPo8L1++nOXLlx8omhPz999/v9bXFhUVMXToUIYPH86IESMYMWJEjcfFxcXN+E5EMqO8LSKFwj30v37yydBK/cQTYXnPnnDqqdk9Vl4U1sC1wM0pywYAG1OWbQf6N0tEIoTCefPmzaxYseJA8bx8+fIDz1esWMHu3btrfX379u0ZNmxY2sJ5yJAh6vsshUx5W0Tyzt69YeSQhQvhlVfCtGBBuHV7Qvv2cP75cNNN2e0GAnlQWJvZGcB8d08dNMUJdxBLVgzsT7ef8vJySpNuCTRjxgxmzJiRzVClhdqxY0eNYjm1eN62Le31sgf06tWLoUOHUlJSclDr84ABAzTyhqRVVlZGWdmBoaQLaqDwbORt5WwRaYyqqnAhYuLW6slD8S1bFm5Ek66X5WGHwcknhxbqj340tFZnoqE5O/ZRQczscSD5DvU9gR3Aj4FPuvsRSdsuBa5393tT96MrzCUdd2fTpk2sXLmSlStXsmrVqhrzFStWsGHDhjr30aVLF4YOHXqgeE48Tjzv1q1bM70baakKbVSQbORt5WwRSeYOW7bAunU1pzVrqovo1avD8/1pm1iDNm3CnRqPOQbGjQvT2LEweHABjGOdJRcTWjQS/gX8F/ACcI2ZdXT33WbWDhgMPB1DjJKnKioqWLNmTdqiOTHftWtXnfsoLi6uUTCnFs+9evXSxYIiNSlvi0idqqrCxYEbN9acNmw4uHhety60QtdVMCc75JBwO/NBg8JdFg8/PNzq/PDDw/P27XP73uoSe2Ht7uXJz82sEih391VmNgc4FZgDTAEedPeD78MsLdaOHTvqLJrXrFlT67B0Cd27d2fw4MEMGTKkxnzw4MEMHTqUfv36URTXvU9FCpDytkjrUVkZCuQtW6qnTZsOLphTp82bQ3GdiW7dwo1cElPfvjBgAAwcWF1IDxwIHTvm5r1mQ+yFdT2+ANxiZmOBQwF1wGtB3J1169YdVCwnP968eXOd+zAz+vfvf1DRPGTIkAOPu3fv3kzvSERQ3hbJK1VVsG1bzcI4ddq8ufZ19VxmVKdu3ULrcuqUXDwnF9H5XDA3VN4V1u5ekvR4M/C5+KKRpti/fz/vvvvugWI5XfG8d+/eOvfRoUOHA63L6YrmgQMH0j7O73xERHlbJAcqK2HHjlDYppu2bq19XerUlMvpzKB7d+jRI1zw16NHmNIVzMlTz57QGge+yrvCWgrHzp07axTNqYXz2rVr6+2mccghh6RtbU7M+/Tpo/7NIiJSMPbvh+3bQ0GbmGdSBCe23bEjezF161ZdECdPyYVybcu7dgX1lmw4FdaSlruzcePGOgvnjRtTh6utycwYOHDgQS3NyS3OXbp0aaZ3JCIikl5VFezcWbMYTkyZPN++Heq4tUHGunYNRXG6qXv32tclb9OlC7RVtddsdKpbqcrKStauXVtn4dyQ0TSSu2mkTgMHDtQNUEREJCfcw81AaitwMymOd+xoWneJZGbVhW3XrrUXwPUVxmopLkwqrFuwbdu28fbbb/POO+8cNF+1ahUVFRV1vr5bt25pC+bE1LdvX42mISIiGamoOLirRENbg1PXNXR4tobo1Ongwra253Wt69Qpe+MlS+FRYV3AqqqqWLt2ba3Fc303PunXr1+dhbNG0xAREQitudnqKlHPl6EZadcufUtvpgWxuktItujXqABs2rSJxYsXs2TJkgPzJUuWsHz58jpH1ejYsSPDhg1j2LBhDB8+vMa8pKSEDh06NOO7EBGR5uYeCtnNm6unxnaVyHRM4tokuko0tkU4+Xlxcf3HE2lOKqzzRGVlJcuXL2fx4sUHFdHl5eW1vq5v3741iubkx4ceeqhG1BARaQF2765ZHG/eHG7SUd+yTZuy110i0VWiqQVx587qKiEtlwrrGGzZsoXXXnuNhQsX8uqrr7Jw4UIWLVpU68WCnTt3ZuTIkYwaNYpRo0YxcuRIRo4cyfDhwzWqhohIgamoCLd1Li+ve9qwobpQ3rOn8cfr0CEMn5aYEhfNZVIQq6uESMPozyTHtm7dyksvvcTzzz/Piy++yKuvvsrKlSvTbjtgwACOPPLIg4roAQMGqOVZRCSPuYfbOK9ZA2vXhinxeN26mgVzPTeUTatdu1AU9+pVs0iu73nPni3jbnYihUKFdRa5O2+//TZPPvkkzz//PM8//zxvvvkmnjKGT4cOHTjqqKM45phjGDt2LMcccwxHH300hxxySEyRi4hIXbZsgeXLq6dVq2oWz++9B/v2NWxfZuHOdH361D317l19BzuNNCFSGFRYN9GGDRuYO3fugSm1Nbpdu3Yce+yxHH/88UycOJHx48czYsQI2uo7NRGRvLJ9O7z5Jrz+epi/8051Id2QVuaePaF//+ppwAA47DA49NCaBXOvXtCmTe7fj4g0P1V3jbBu3Tr+8pe/8Kc//YmnnnqKysrKA+t69erFKaecwoknnsjEiRMZN26cRt8QEckza9fCCy+E6dVXQzG9alXt23fqBEOHwrBhYT5kSCickwtpdbkQERXWDVRVVcXcuXOZOXMms2fPpioad6ht27ZMmTKF0047jSlTpjBu3DjdNEVEJM+8+y48+ig8/jj861+wevXB27RvDyNHwpgxMHo0jBhRXUj37auuGCJSPxXW9aiqquKBBx7g+uuvZ8mSJUDo3vHhD3+YadOmcc4559CzZ8+YoxQRkVSrV8Pvfgd//CMsWlRzXbduMGECTJwIxx4biunhwzXyhYg0jVJIHRYsWMBll13GggULABg8eDBXXHEFl156Kf369Ys5OhERSWf+fLjhBpgzJ4zWAWHs5FNPhbPOgg99CEaNAn25KCLZpsI6DXfn1ltv5brrrqOiooIBAwZw/fXXM336dNq1axd3eCIiksamTXDVVXDvveF5+/Zw3nnwmc+Eolp36RORXFNhnaKqqoorr7ySO+64A4CrrrqKH/zgB3Tu3DnmyEREpDZvvAFnnBH6UnfsCFdeCddcE0bhEBFpLiqsU1x//fXccccdFBcX84c//IGPfexjcYckIiJ1WL0aJk8ON1+ZOBHuuSdceCgi0txUWCf5xz/+wU033URRURF/+ctfOOuss+IOSURE6uAO06eHovrUU2H2bA17JyLx0aUbkYqKCq644goAvvOd76ioFhEpAM8+C//4R7g5y333qagWkXipsI488MADLF26lBEjRnDttdfGHY6IiDTArFlhfsUV6k8tIvFTYR2ZOXMmAF/96ldp3759zNGIiEhDzJsX5lOnxhuHiAjkUWFtZueb2XIz22hmPzWzttHy3mb2KzP7oZn92MyyPt7du+++y7PPPkuHDh24+OKLs717EZEWJ86cneAOS5eGx0cdlaujiIg0XF4U1mY2GDgPmAZ8Bbgc+HK0+gHgLne/FngLuCHbx//rX/8KwIc//GG6du2a7d2LiLQocefshO3bYfdu6NIFunfP1VFERBouX0YFGQJc6u4VwMtmdhRwspk9C4x09xei7WYDS8zsJnffma2D/+Mf/wBCYS0iIvWKNWcn7NgR5moPEZF8kRct1u4+L0rQCWuAVcApwPKk7dZED0uzdeyqqiqeeuopAE455ZRs7VZEpMWKM2cn2749zLt0ycXeRUQylxeFdRrHAb8EBgAbU9ZtB/qnvqC8vJzS0tIDU1lZWYMO9Oqrr7J582aGDBnC0KFDmxq3iEiDlJWVHchXQO+442miZsvZydRiLSLNpaE5O1+6ghxgZiOAcndfaGYO7E7ZpBjYn/q6Pn36MH/+/IyP98IL4RvLE044IfNgRUQaacaMGcyYMQMAM9sQcziN1tw5O1misFaLtYjkWkNzdl61WJtZG+DzwNeiRWuAnknrDegMrM3WMV988UUAJk6cmK1dioi0CnHk7GQqrEUk3+RVYQ18FbjV3fdGzx8HRiatLwH2AU1r5kiSKKwnTJiQrV2KiLQWzZ6zk6mPtYjkm7wprM3sOmAB0NHMhpnZpcAW4L3oinOAM4Gfu/u+bBxzx44dvPHGG7Rt25Zx48ZlY5ciIq1CHDk71bZtYa4+1iKSL/Kij7WZfRO4MWXxYne/y8wuBL5tZiuALsA3snXchQsXUlVVxdixY+nQoUO2disi0qLFlbNTbd4c5j171r2diEhzyYvC2t1vAm6qZd1K4LJcHHfJkiUAjBkzJhe7FxFpkeLK2alUWItIvsmbriBxeOuttwA44ogjYo5EREQypcJaRPKNCmvg8MMPjzkSERHJ1KZNYa7CWkTyRasurN955x0ARowYEXMkIiKSqUSLda9e8cYhIpLQqgvrtWvD0KoDBgyIORIREcnUli1h3qNHvHGIiCS02sJ6//79lJeXU1RURN++feMOR0REMrRnT5h37BhvHCIiCa22sH7//fcB6NevH23atIk5GhERydS+aHTs4uJ44xARSWi1hfV7770HwGGHHRZzJCIi0hh7o/s9tm8fbxwiIgmttrBOtFirsBYRKUyJFmsV1iKSL1ptYb1hwwYAevfuHXMkIiLSGCqsRSTftNrCeuPGjQAccsghMUciIiKNoT7WIpJv8uKW5nEYOnQo55xzDuPGjYs7FBERaYQPfxh27VJhLSL5o9UW1tOmTWPatGlxhyEiIo00a1bcEYiI1NRqu4KIiIiIiGSTCmsRERERkSxQYS0iIiIikgUqrEVEREREsqDVFtZlZWVxh9Bkeg/5oSW8B2gZ76MlvAepm37G2aNzmT06l9lT6OdShXUB03vIDy3hPUDLeB8t4T1I3fQzzh6dy+zRucyeQj+XrbawFhERERHJJnP3uGPICjMrB1Zm8JLewIYchdNc9B7yQ0t4D9Ay3kehvoch7t4n7iCaUyNydkKh/ozzkc5l9uhcZk8hnMtac3aLKaxFREREROKkriAiIiIiIlmgwlpEREREJAtaXWFtZh3N7Odm9n0zu8PMusUdU6bM7HwzW25mG83sp2bWNu6YGsvM2pvZQjObHHcsjWVmx5vZV83sPDPrHXc8mTKz0Wb2CzP7ipnNNLNxccfUEGZ2mpm9YGYlScsK/u9b0tPPtvFq+8wws95m9isz+6GZ/djM2sUda6FI/ezS72fTpH6OFvL5LNiCrAlmAo+6+ywzOxP4FXBRzDE1mJkNBs4DpgFjgDuA1cCtccbVBNcAJXEH0Vhm9llguLt/M+5YmuD3wLnuvib6/XoMGB1zTHUys35AF2BCyqqC/vuWOuln2wj1fGY8AHzN3V8wsyuAG4CvxxVrgUn97NLvZyOl+xw1s99SoOezVV28aGb9geVAT3ffFf3XvhUY4+4rYg2ugcxsEvAvd6+Int9CiH9qvJFlzsz+AxgFfBuY7u5PxRtRZszsJOBbwGlewH9IZrYT+IC7LzazPsBCd+8fd1z1MbMioBIY6u4rWsLft6Snn23j1faZAdwI/Dnxt25mA4AlQD933xlXvIUg9bMLeAv9fjZKus/RQv97b21dQSYDG9x9F0CUaJYDJ8UZVCbcfV4iQUbWAKviiqexzKwzMM3d74o7lib4MfAm8HMze9TMTog7oEb6E3CnmXUFLgGuijmeBnH3qpRFkynwv2+p1WT0s22UOj4zTiGcw8R2a6KHpc0YXsGp5bNrMvr9bKx0n6OTKeDz2doK6wHAxpRl24G8b52rw3HAL+MOohGuBW6OO4jGMrORwHjg1+5+JfAE8GjURaHQfBHYA7wE7HD3B2OOp7Fa4t+3BPrZZk/iM0PntHHSfXbpXDZCbZ+jwCAK+Hy2tsLagd0py4qB/THE0mRmNgIod/eFcceSCTM7A5jv7uvjjqUJxgCb3P3V6PntgAHnxhdSo3UE7gPuBX5iZlNijqexWtTft9Sgn20WpHxm6JxmqI7PLp3Lxqntc7Sgz2dru3hxDdAzZVlXYG0MsTSJmbUBPg98Le5YGuGrwHgzSzzvCTxkZt9z91viCysjbUn6+3H33Wa2FDgkvpAa7R7gE+6+xcwcuM/MSgqwn2WL+fuWg+hn20RpPjPWAEckrTegMzqndUn72UXozqDfz8zV9jnahgI+n62txfpJoL+ZdQSIhhYaDDwda1SN81XgVnffG3cgjXAxMC5pWgtcTmF1aXkV6J7S9aMCeCOmeBolGh5wrLtviRbdBGwjz0cFqUVL+vuWmvSzbbrUz4zHgZFJ60uAfcD8Zo6rkNT22fVb9PvZGLV9jq6lgM9nqyqs3f19YA5warRoCvBg0kUbBcHMrgMWAB3NbJiZXRp9xVcQ3L3c3d9NTISRHcrdfVvcsTWUuy8G/kYYwgoz6wG0Bx6OM65G2ATsiUYESNhIuMo9r1l1s5FBy/n7loPpZ9s06T4zgC3Ae2Z2VLTZmcDP3X1fXHHmuzo+u1ah38+M1fE5+gcK+Hy2quH2AMysJ3ALsAI4FPi6u++INagMmNk3CcMkJVvs7oXYwgiAma2gMIfb6w38lNDCM5hwAUZBtVgDmNlYwgWMLwP9gGfcPa9bBpJGMJlJGHv3dncvL/S/b6mdfraNU9dnhpkNIQwZt4IwLvw3UkYQkTokf3bp97NxavscLeTz2eoKaxERERGRXGhVXUFERERERHJFhbWIiIiISBaosBYRERERyQIV1iIiIiIiWaDCWkREREQkC1RYi4iIiIhkgQprkVqY2alm9pyZTa5nu6PM7H4z+1YWj/2jbO1LRCQOZnaMmT2QzdxYy3HmmNkncnmMXDGzKWb2r/o+Z6RwtK1/E5F4mdlo4HuEuwQeBxwDfMrd78nxoVcDYxuw3SqgJ9AmGweNbt+6Pxv7EhEBMLPjgB8AHwL+Diwn3HijHXCHuz+Wg8MuA7qTpdxYh5nAolzsOMrH5wO/JNyR9l/RqmLgY8AEd1/RhEM09HNGCoQKaykEfwU+7u4LzawI+G1zHNTd3zKzjQ3YbpuZvZfFQ58EPJPF/YlIK+fuL5nZH4FR7v6RxHIzmwr8xcy+6u4zs3zMXVnOjQCY2dnAa4mC1t0fyfYxEtx9PzDLzK4CnnL3bybF8VBj9mlmfYGT3X2Wuy8xsw1ZClfygLqCSF4zs17AEcBOAHevItyCt7luGVqV5e0a4mTgySzuT0QE4KDblbv7w8B3gP9nZofm4JhZzdVmNpDQQt3c0t3qfS6wJpOdmFkx8DugY9LibH5+SMxUWEtec/dNwBLgUTM7Plq2EnjRzHqY2Y1m9qKZnWNm75vZm2Z2bOL1ZnaRmX3DzJ4ys+8nLS81s++Y2Swz+7OZdY+WdzCzX5jZt83sN4SvMTNiZv8TxfRxMzvOzP4a7e/GKMZZZtYnWr7VzK5P2UV7d98b7eu/zey/zOyfZnZzxidQRKR+vyF0CbkAwMyGm9kNZnanmT0RFbOYWXGUy75hZvPM7MLEDszsi2b2LTN73syurO1A6XKymU2LcuZ0M1tqZjebWX8z+22UT+eZ2SnRLi4GBgLXmtl50bUwc83sM0nH+JyZfc/Mfmdmd5tZZzPrbmbXm9l8MzvFzJaY2TtmVtKYExadk8lRizZmdr6Z/dTMbjOzR82sS/Qefmpms6PPm/eic3wkcIGZ/VfSLntGn0XbzexztRyzyMx+aGYzomO4mf0l8fORPOHumjTl9QSMAVYS/qu/C+gbLW9DSFKbCH3degCPA68T/mn8IPDNaNuehFbvqUA34I9J+/8n8KPo8f8DZkSP+wGVhORZX4z/S2j16Uvoi9cpWl4MPAY8BPQnfCDsB64HOgEfiJ4ntk/0H0+875nR467A1+L+WWjSpKlwJ2A68G4t69YBP4vy6p+BNtHye4EHosc/Bs6PHl8KlEePLwIuiR4fEeXqo6Pn/wt8J3pcW07uBWwBfgJMAKYAtwLfjra9BpidFKsDJdHjjsCbwPTo+RnAvUnbzgF+Rej6+tHo8+L86H0+Cnyvgefuqeiz4mbgR4T+1lOS1m8ABkePXyd8JhUD3wDeBSYCl0V5/6lEvNH2K4DbonUXAatqieFyYE70uF90nk+P+/dKU81Jfawl77n762Z2NHAT8EXgbDOb4qHPdTmwzd0fBDCz7wDPAsMIHyJVZvblaFePERL4VOCQpOVvRq/tAnwhei3uvs7MVmcQam/g58Bl7r4r2sdeM3sfWO7ua6PjvA887aH/4QJCwu9NuAhyKuFDAGAX8Gkzewu4g/BPhYhILjih2JwADAGuMjOAjUB3C08+T2gUALib0BUCQq6dn5JrDwNeSznGdNLkZHffZGZbgYfc/UUAM1sF7DSzwwgtvF3SBu2+28zWJy26Mtpvwp3A/dHyLYTPiz9Hx5gfxdlQT3rUx9rMZhEacxIucPdVZnYCoaDuEuX/tcDb7v4C8EL02nT7nh19JrxUR0zHAruj973OzDZRs0uJ5AEV1pLXzKwNMMjDRSr/GXXPeAQoI7QApHo9mvcABgO/d/d7o2U/ifZ5LaHQ/UnKsT5ASIjbGxnuEcCJwHBgYdLy1D6GlQdWuFdFSTbRLesQD91fcPflZjYduB34CvBJYD0iIllkZj0I37a9TsibG9Pkx36EFtX2AO5eSWgMIHrNd939ueh5jdcmSZuTI07NXPkucC2wFHgeGFrHW0h+3eHAP5Kev03o5tKXg3NxBY0fseTfhH9AEraY2Q+BBwixJ6rn1PdVnwpqr83+BtxhZh0IP4vdhNZvySPqYy357kC/PwB3XwhcB4yqZfv2hCS2DHgPmJa80swmRsunmlnHlOXboqejGxnrc8BPgT9EiS8jZtafpAthzGwAYUSUIwgXM95by0tFRJriM4RvyB4g5McTopZi4EB+3ADsA85MWn60mXUlJdeaWduooSJVbTk5nV8CKzwMq7o3g/eyitDCfeAQwFZCV5escfeqqPFjTPSPySPADe4+P5vHSTnmbGA24ZvVi4CT3H1rro4njaPCWgrBVVHRmTCA0DcuoZeZtY8eTyb0r9sC3Ad81MxuN7PJUWvCTkIC7AzMMbMzogtFRhGK8QXAjWbW3sx6E/o29zWzNmb2NTO7pJYYLZq+Tfjw+XGadanbpz7/CKE/YMLhwHnuvg24Ks1rREQycVBLqJl9CPgu8Hl3Lye0Dq8DHjGzs6MLFCdFLdR/JIwecp6ZnQp8xt23E3Ltf5rZd83sJEIDQ6IbXXL+qy0nJyS3Ho8H+phZJ+AEoKOZDYvW7Sdc7JdoYEk+xi+Bj5lZz+j5RODX7l5BqHnS5mIzOz6Kq7iWc1dESs0U/fMxlfAtZU9CF8ORwCCgU9KFkamt4vvSxF8zqDT9RaKfRTfgCUKrfCb/cEgzUVcQKQSDgNfN7GHC12QGXJG0vgL4hpltIbTuXgXg7o+b2X8Svk48j3Dx3yIAMzsX+AXhg+Iu4Bp3dzP7OGEopNeBWYSv9MYRbqjwcUIXkxo3pjGzCYSLcgYTvhr8K/DdKJ7HgOOBEjMbTmhJORT4uJktprr1ZjrQ2d3LUt77by3c2KGCcLGQiEjGojxyMaH4+zVQTrgArjtwlrv/C8Dd90X58deEXPcQoW81wNWERonfEb6hS4zE8RtC7vsCocvaF9x9fVQ4Hg8MN7N7asvJZjaN0K/4cjNb7O5rCBdS3kooqn8HnAscDbwD/J7Quv5pC6NFjQE+bGaPuPufzexw4EEz+yfhW89vRC3rFwCHRu/vdcLNcg4zszFRnJ8g9Bv/a9J5a0fI/WOBAVHh7YS+zedH+3wFeBp4idC95W+EnD6HcMHkWDM7z90T+51FuFnPjqgQPxSYZmbLCS3RRK9PvWfDKsLFmecBHYAiM5tLuICxuYaglXqYfhZSyCzcBvZ/3b2kmY53jbvrduMiIi1M1C2l2N3z8gZdZnYNYYSWFdHzDoRvG36YuDZH4qeuIFLo0nWzyM2Bwlem/26OY4mISPOJLs48No+L6g7AN4GSpMW9gZ0qqvOLuoJIwYr6t11C+GrvAnd/IMeHfDv6ilJERFqWtu5+R9xB1Mbd91i4rfrPogvb1xK6rXyz7ldKc1NXEBERERGRLFBXEBERERGRLFBhLSIip3rd3gAAIABJREFUIiKSBSqsRURERESyoMVcvNi7d28vKSnJ+HXl5eX06dMn+wHlkGJuHoUYMxRm3K095pdffnmDuxfWCWii1pSz46Dz1DA6Tw2j81RTXTm7xRTWJSUlzJ+f+Z1ES0tLG/W6OCnm5lGIMUNhxt3aYzazlVnZUQFpTTk7DjpPDaPz1DA6TzXVlbPVFUREREREJAtUWIuIiIiIZEGrL6xnzJgRdwgZU8zNoxBjhsKMWzFLQ+m8N4zOU8PoPDWMzlPDNfsNYszsNOAm4MKk+92fD9wGdAPuAf7b3Suidb2B7wFbgHbAte6+P3W/paWlrv4/IlKIzOxldy+NO450lLNFRGqqK2c3a4u1mfUDugATkpYNBs4DpgFfAS4Hvpz0sgeAu9z9WuAt4IZmC1hEpBVTzhYRyUyzFtbuvg54KGXxEOBSd3/Z3X8H/AI4GcDMjgdGuvsL0bazgavMrHNzxSwi0lopZ4uIZKbZ+1i7e1XK83mJrxAja4BV0eNTgOVJ266JHmbvK9O774Zjj4Wf/jRruxQRaSnyJWdv//d2Xj7uZRZfuripuxIRyZl8HMf6OOBH0eMBwMaU9duB/qkvKi8vp7S0OnfPmDGjYZ3tt22DBQvggx9sbLwiIhkrKyujrKws8bR3nLE0UbPk7Kp9VWyfvx3aNDleEZGMNTRn51VhbWYjgHJ3XxgtcmB3ymbFwEEXwvTp06dxg5f37Rvm69dn/loRkUZKLiTNbEPM4TRKc+bsNp1CRV21q6qeLUVEsq+hOTtvhtszszbA54GvJS1eA/RM2saAzsDarB24X78wX7cua7sUEWnpmjtnF3UKH1eVuyqbuisRkZzJm8Ia+Cpwq/v/Z+++4+Qq6/7/vz7bN9uz2WwKqRBCEggJBEJvgjRvmiAIFkSJivpVb/GnFFEBbxtyY+OWWG8UkFuwICoKSBWEVEIoKZAN6dmU7X33+v1xndmdnWzf2Sk77+fjcR4zc87MOZ85mfnks9dc57pcc9i6x4HZYY+nAy1A9MZoUou1iMhQxDRnq8VaRJJBzAvroAUDwMLW3QisAnLNbKaZXWNmhzjnlgM7zOzw4KnnAD90zrVELSC1WIuI9CpRcrZarEUkGcS0j7WZFQAfCB5ebWY/wv+UeFvEU990zv0iuH85cIuZVeDHU70pqkGVlEB6OlRVQUsLZGVFdfciIskqkXK2WqxFJBnEfObFkTKsWbwmTYIdO2DLFjjooOgGJiLSj0SeeXGkDDZnO+d4JvMZaIdTWk4hLTORejKKSCpJmJkXE9aECf5227a+nyciInFhZqTnBq3WjWq1FpHEpMIaYPp0f1tREc8oRESkD+pnLSKJToU1wIwZ/laFtYhIwlI/axFJdCqsQS3WIiJJQC3WIpLoVFhDV2G9aVNcwxARkd6pxVpEEp0Ka4CZM/3txo3xjUNERHrV2WJdrxZrEUlMKqwBZs2CjAx4+22or493NCIi0oP0At9i3V6rwlpEEpMKa/CTwsyZA87Ba6/FOxoREelBZkkmAG1VbXGORESkZyqsQ+bP97evvhrfOEREpEcZxX6y4Lb9KqxFJDGpsA454gh/u3p1fOMQEZEedRbWarEWkQSlwjrkmGP87b//Hd84RESkRxklKqxFJLGpsA5ZvBjS02HVKqiri3c0IiISQS3WIpLoVFiH5OXBUUdBezu89FK8oxERkQgqrEUk0amwDnfSSf726afjGoaIiBwo1BWkdX9rnCMREemZCutwZ53lb//yl/jGISIiB+hssd6nFmsRSUwqrMOdfrrvErJqFWzZEu9oREQkTFZ5FgAtu1riHImISM9UWIfLyYF3v9vf/+Mf4xuLiIh0kzk+EwxaK1vpaOuIdzgiIgdQYR3p8sv97b33xjcOERHpJi0jjcyyTHDQulv9rEUk8aiwjnThhVBUBMuXw9q18Y5GRETCZE0MuoPsUHcQEUk8Kqwj5eTAFVf4+7/6VVxDERGR7rImBIX1ThXWIpJ4VFj35Jpr/O0vfqHJYkREEkj2xGwAmnc0xzkSEZEDqbDuybHHwvHHw/798MtfxjsaEREJZB8UFNabVViLSOKJeWFtZmeZ2UtmNj1sXa6Z/dDM/svM/sfMCgeybURdf72//d73oFkJXERSU6Ll7JyDcwBofKsxmrsVEYmKmBbWZlYO5APHRmy6G3jeOXcj8CfgngFuGzkXXgjz5sHmzXD33TE5pIhIIknEnJ17cC6gwlpEElNMC2vn3C58ou1kZpOAK4E/B6ueAC4ws+l9bRvxYNPT4Vvf8vdvvx2qqkb8kCIiiSQRc7YKaxFJZDHvCuKcixzV/zRgj3OuIdjeBmwCTu1n28g7/3w49VTYtw++/vWYHFJEJJEkWs7OmphFWm4abXvbaKvW1OYiklgS4eLFycDeiHW1wKR+tnVTWVnJokWLOpelS5cOPzIzuPNOSEuDH/wAVqwY/j5FRIClS5d25itgXLzjGYS45mwzI2em72fdsKFhkKGLiAzNQHN2RuxC6pUDIn/TywZC02r1ta1TWVkZy5cvj350Rx0Fn/ucL7CvvRZefhkyEuG0iUgyW7JkCUuWLAHAzPbEOZzBiHvOzjs8j4bXGqhfU0/hothczy4iqW2gOTsRWqy3ASUR6wqA7f1si52vfx2mTYNVq+Cb34zpoUVEEkzcc3b+gnwA6l7RPAMiklgSobB+CphkZrkAZpYJTAWe6Wdb7OTnw89+5u9/7Wvw3HMxPbyISAKJe87OPzIorFersBaRxBKPcawtdBfAObcTeBR4V7D+TOBh59y2vrbFMOTgyGfCl78MHR1w5ZVQWRnzEEREYi0Rc3Z4i7VzLpq7FhEZlliPY10AfCJ4eLWZlQX3PwlcaGY3AecBS8Je1te22Lr1VjjhBNi6FS66CJqa4haKiMhIS9ScnTUhi6wJWbRXt9O4QcPuiUjisNHy1/6iRYvciFy8GGnHDli8GLZsgfe/H+67z48eIiIyRGa2wjm3KN5xxNJwc/baS9ey5+E9zP7ZbCZ+dGIUIxMR6VtfOTsR+lgnl4kT4dFHoaAAHnjAT30+Sv44ERFJFsWnFANQ9awm7xKRxKHCeijmz4ff/Q4yM/0wfDfeqOJaRCSGik4uAqD6ueo4RyIi0kWF9VCdfbYvrjMy/NTnt9yi4lpEJEby5+eTXpRO06YmGjepn7WIJAYV1sNx4YW+O0h6Otx+O1x3HbS3xzsqEZFRz9KNse8eC8DeRyMnexQRiQ8V1sN16aXw0EOQkwM/+Qlcdhk0qvVERGSklV5QCsDeR1RYi0hiUGEdDRddBI8/DsXF8Ic/wOmnw7bYD7UtIpJKSs8thTSoerqKtuq2eIcjIqLCOmpOOgmefx6mToWXXoKjj9YMjSIiIyizNJOik4twbY7KP2jSLhGJPxXW0TRvHixfDmecAbt2+du77tJFjSIiI6T8qnIAdv3vrjhHIiKiwjr6ysrg73/341u3tcHnPw8XXKAp0EVERsD4940nLSeNqqeraKzQ9S0iEl8qrEdCRgZ897v+osbiYj+hzJFHwhNPxDsyEZFRJaMog3EXjQNg5y92xjkaEUl1KqxH0nvfC2vWwCmn+KnQ3/1u+OIXNWqIiEgUTfrEJAC2/2Q77Y0a8lRE4keF9UibMgX++U+49VZIS4M77vAzNz79dLwjExEZFYpOKSL/qHxaK1vZdZ/6WotI/KiwjoX0dPjKV+Bf//IXOG7c6IfkW7IEqqriHZ2ISFIzM6b85xQAtt65FdeuC8ZFJD5UWMfS4sWwciV8/euQmQk//SnMng0/+5lmbBQRGYay95WRPTWbhjca2P3g7niHIyIpSoV1rGVlwS23wOrVcOKJsHs3XHstLFqk7iEiIkOUlpnG9K9OB2DTLZvoaO2Ib0AikpJUWMfL3Ll+Apnf/tZPKrN6te8ecskl8Prr8Y5ORCTplH+onNzZuTS91cSOpTviHY6IpKCoFNZmdl409pNyzODyy+HNN+H22yEvz0+Jfvjh8IEPwIYN8Y5QREah0Zqz0zLSmPlfMwHYdPMmWna3xDkiEUk1gyqszewcM3vezDaa2dvBUgE8PDLhpYjcXLjpJli/Hj75ST8O9n33wZw5cM01sGlTvCMUkSSUijl73MXjKDm7hLaqNt7+0tvxDkdEUsxgW6x/DfweuBb4SLBcHayX4Zo0Ce6+2xfYH/2oX/fLX8KsWb4Fe/Xq+MYnIskm5XK2mTHrR7OwbGPnr3ay/8n98Q5JRFLIYAvrt51zdzrnnnLOPRMsTwO3jEBsqWv6dD9SyJtvwgc/6Nfddx8sXAhnnQX/+Ac4DSclIv1KyZw95pAxTLt5GgBvfPgNWve1xjkiEUkVgy2sP2lml/aw/oJoBCMRDjkE7r0X3noLPvtZ3wf7iSfg7LN9kf3LX0JDQ7yjFJHElbI5e+qXp1J4fCEt21pY/4n1ODVGiEgMDLawPhH4mZm1hy0dwP8MNxAzm2NmPzazz5vZ3Wa2IFifa2Y/NLP/MrP/MbPC4R4r6UybBnfdBVu2wDe+AeXl8Morvv/15Mnwuc/BG2/EO0oRSTwpm7PTMtKY85s5pOenU/m7Srbfsz0eYYhIihlsYX0rsASYBcwMW74WhVh+DfyXc+6/gW8BDwTr7waed87dCPwJuCcKx0pOJSVw441QUeFbqxcv9jM3fv/7fvi+00+HBx+E5uZ4RyoiiSGlc3buzFwO/cmhAGz8zEaqntNMtyIysgZbWP8N+JNz7m3n3OZgqQB+GoVY5gAFwf1GoMjMJgFXAn8O1j8BXGBm06NwvOSVkwNXXw3//refyXHJEt9N5Omn4Yor/EWQn/oUvPSS+mKLpLaUz9nlV5Vz0OcPwrU5Xrv0NZo2N8UrFBFJAYMtrH8L3GBmp4QtpwI3RCGWh/A/WRYAHwA+A5wG7HHONQA459qATcCpUTje6LBwIdxzD2zbBj/+McyfD/v2+dFFjjsODjvMdx/ZvDnekYpI7ClnAzO/M5OSM0to3d3KK+9+hZZKjW8tIiNjsIX1V4DPAfcC/xu2fDwKsVwHNAHLgDrn3MPAZGBvxPNqgUlRON7oUlQE113n+16vXg3/+Z8wYYIfuu/mm/1II8cfD3feCe+8E+9oRSQ2lLPx/a3n/m4uefPzaFzfyKvnvUpbbVs8QxKRUWqwhfV3gSnOuenOuRnBMh24JAqx5OL76N0P3GVmZwIO/xNjuGzggLGTKisrWbRoUeeydOnSKISUpI48Er73PX+x49/+Bu9/v5+E5t//hi98wV8MuXgx3HGH768tIjG3dOnSznwFjBuhwyhnBzKLM5n/2HxyZuZQu7zWF9c1Kq5FZGAGmrNtMEMQmdnxzrkXe1h/jXPuF0OKtGsfjwFXOOeqzOwrwP/Dt7R81Tl3aNjzNgTr7g9//aJFi9zy5cuHE8LoVl8Pf/0rPPQQPPpo92H6FiyA//gPvxx9NKRFZaZ7ERkgM1vhnFs0AvtVzo7Q+FYjq09bTfPWZgqOLWD+Y/PJLMmMeRwikrz6ytmDraCu62HnlwF3DSWwsH2MA450zoUu2b4dqAE2A5PMLDd4XiYwFXhmOMdLSXl5cNllftSQykpfYF9+uV+/ejXcdhsce6wfvu/aa+GRR3wxLiLJTDk7Qu7BuSx4dgE503OofbmWV854hebtGklJRKJjsIX12KBlAjMbY2Y/A36G72M3HPuAJjObHLZuL7AGeBR4V7DuTOBh59y2YR4vtY0ZA+99L/z2t7Bnj+8uct11MGUK7NzpZ3288EIoLfUzPX73u77vdkdHvCMXkcFRzu5B7gxfXOfOyqVudR0rF6+kbk1dvMMSkVFgUF1BAMzsOOAq4GxgJ/AhoNo5t39YgZgdiW9dWQGUA886554xsxLgO0AFMAG4wTl3QAZUV5AocM4X0H/+s1+WRfzfO368L7RDyyRdQyoSDSPVFSTYt3J2L1oqW1h70VpqXqghPT+duQ/OpfS80rjGJCKJr6+c3WdhbWan9LLpWOB04A6gA7jIOff54QY6HImQpEed3bvhySfhH/+Axx/3Q/qFmzfPF9innw4nn+wnsBGRQYtWYa2cPXjtTe2s+8g6dv92NxhMu2Ua078yHUu3eIcmIglqOIX1v/GTAFThk3GPTwMmOueyhxvocCRKkh61nPPTpj/+uC+0n366+wWQZn4M7VNP9cspp8C4kRroQGR0iWJhrZw9BM45Nt++mYqvVoCD4ncVM/e+uWSVZ8U7NBFJQMMprM8DljvndvdzgEudcw8NL8zhSaQknRKam+HFF32L9jPP+FkeWyImXZg3r3uhPWFCfGIVSXBRLKyVs4dh3xP7eOOqN2jd3UrWhCxm/3y2uoaIyAGGXFgnk0RM0imlsdEX188845cXX4SmiKmDZ86EE07wy/HHwxFHQHp6fOIVSSAj2cc6USVqzm7e3szrV75O9TPVAEy4ZgKH3HkIGUUZcY5MRBKFCmuJveZmf/FjqNB+4YUDh+/Lz/cT1YSK7eOOg+Li+MQrEkcqrBOLa3ds+e8tbLp5E67ZkX1QNocuPZTSc9V6LSIqrCURtLXB2rW+wA4tmzYd+Ly5c32BfcwxsGiR77edpX6OMrqpsE5M9W/W8+aH36T25VoAxl08jkP++xBypuXEOTIRiScV1pKYdu70XUZChfby5Qf2087K8sX1okVdxfbcuZChn2Vl9FBhnbg62jrYetdWKr5WQUd9B2m5aUy9cSpTrp9Ceo66somkohEtrM1sIfAFYBV+IoCKYe1wiJIlSUsfmpth5Up4+WVfZC9bBuvWHfi83FxYuNAX2YsW+WnYDz1UxbYkrVgW1srZQ9O8rZm3rn/LD8sH5MzIYcZtMxj//vFYmobmE0klI11Y3wNc6JybYGZXAn9wzjUOa6dDkGxJWgaouhpWrfJF9vLlfnn77QOfl5PjL4ZcuBAWLPDLEUf4ftwiCS7GhbVy9jDsf2o/Gz6zgYbX/HCjefPzmPlfMxl73ljMVGCLpIKRLqwzgfHOuW1mlgbkO+dqhrXTIUjWJC1DsHcvrFjR1aq9ahVs3nzg88xg1qyuQju0TJwY+5hF+hDjwlo5e5g62jrY9etdVHy1guYtzQAUnVzE9K9Np/j0YhXYIqNcVAtrMzsXGOucu8/MjgXqnHOvRyHOYUnmJC1RsH+/n4599equ5bXX/EWTkcrLDyy2Z83S0H8SNyM8pbly9ghpb2pn+93b2fyNzbTt87mm8LhCpt44ldL3lKrAFhmlolZYm9k3gM8CjznnLg3WfQtY45y7PxrBDtVoSNISZc3NfrbI1at9q3ao4K7poXEuN9dfFHnEEd2X8nLf8i0ygkaqsFbOjo226ja2/mArW7+/lba9vsDOm5/H1BumUnZpGWkZaXGOUESiKZqF9VPAxcCHnXPfD9adCvzaOTc1GsEO1WhK0jKCnIOKiu4t26tWwZYtPT9/3LgDi+1589R3W6JqBAtr5ewYaqtrY8dPd7Dlji20bPcjHGVPzWbypyYz8WMTyRybGecIRSQa+srZgx1G4XHnXJWZhVfjlw49NJEYM4MZM/xy8cVd6/fv9+Nsv/pq92XPHnjqKb+EmznzwIJ71iyNTCKJRjk7hjLyM5jy+SlMvm4yO+/dyZbvbqFxQyNvf+ltKr5WQfkHypn8mcnkH6E/zEVGq8G2WF8BzAQmA88Bl+FbQz7jnPvxiEQ4QKOx9UPizDnfkh1ZbL/5JrS2Hvj87GyYM8d3KQktc+bAwQdDplqqpHcj2GKtnB1HrsOx7+/72PaDbex7bF/n+qJTi5h07STGXTKO9Fxd2yGSbKJ98eJxwNXAVGAH8IBz7onhBjlcqZCkJUG0tMD69QcW3D2NTAK+qD700O5F95w5fl2OZnCTEb94UTk7ATSsa2Dbj7ax81c7aa9rByCjOIPxV41n4kcnUrCwIM4RishAjfRwe+nAac65J4e1o2FKtSQtCaimxo9E8sYb8PrrXUtvBXdamm/NjmzlPuwwyMuLbewSVzEebk85O47aatrY/cBudvx8B7XLajvX5y/MZ+JHJ1J2eRlZ47LiGKGI9GfIhbWZPQL01xlsHLDHOXfG0EMcvlRN0pIE6ut995HXX+9edL/1FnR09PyaadO6iuzZs/1y2GEapWSUilZhrZydXOrW1LHj5zvY9etdtO33o4lYhlFydgnlV5ZTekEpGfm6bkMk0Qzn4sU1wXPWAR3AtfhpcMOz4QxgbBTiFBmd8vL8tOtHH919fVMTbNhwYMG9fr1v5d68Gf72t+6vKSzsXmiH7s+apW4lAsrZSSV/fj6zvj+Lmd+eyZ4/7mHXvbvY94997PuLX9LGpDHuwnGMv3I8Y88eS1qmhu0TSXT9FdbfxbdqVwGY2cHOuVsin2RmfxqJ4ERGtdA07Ecc0X19a6uftv31131L97p1fnnzTaiq8rNNLlvW/TVmMH16V6EdXnhPnKhW7tShnJ2E0nPSKb+inPIrymmpbKHyd5Xsun8XNf+qYfcDu9n9wG4yijMovaCUskvKKHl3iS56FElQfRbWzrnqiFXlZpbunGsPrTCz04BjRyA2kdSUmdlVHIdzDioruwrtULG9bp0vxDdt8stjj3V/XUFB94I7VHTPmuUnxpFRQzk7+WWVZTH5uslMvm4yjZsa2f3b3ey+fzf1a+vZde8udt27i7S8NErPLWXce8dRel4pGYXqLiKSKAb7bXwIeMXMngdqgTnA2cABLSIiEmVmMH68X04+ufu2lhZfXIe3cIcK7337YPlyv0Tub+rUrq4khx7qb2fN8q3fGpN7NFDOTmK5M3KZdsM0pt0wjYZ1DVT+vpI9v99D7fJaKh+qpPKhSizLKDmzhNL3lFJ6fik5U9UlTCSehjLc3iTgOmA+UAP8zjkX958VdSGMSC/27DmwhXvdOn/xZFtbz6/JyPCT4IQK7dBy6KEwZYof0USiZoSH21POHmWaNjdR+QdfZFc/Xw1h/43nHZ7H2PPHUnp+KYXHF2o6dZERMNLD7aUBC51zK4a1o679HQecBGwEngfqge/gW1tKgC8552oiX6ckLTJIra2+68i6df4iyg0b/IWTGzb0PsU7+IlwDj64e7Eduj9pkvpzD0GMh9tTzh5Fmnc2s+8v+9j7l73sf3x/5xjZABklGYw9eyxjzxtLyVklZE/IjmOkIqNH1AprM/tFD6vLgHbn3EVDjC98/x8BDnbO3Ry27pfAY865B83sHODDzrn3R75WSVokihobfYt2qNAOL7x37uz9dWPG9NzKPWsWlJWp6O7FCM68qJydQjqaO6h6rqqz0G7c0Nhte94ReZScWULJWSUUn1JMep4ugBQZimgW1s8BT9LthydOBJ5zzt02zCBPBb4CnOWCoIKfMDcBJc65BjPLAKqBec65ivDXK0mLxEhtLWzc2L2FO7Ts2dP76woLe27lnjULxqb26G8jWFgrZ6ewhg0NviX77/upeqaKjsaucfMtyyg6oYiSs0ooObOEgqMLsHT94SsyENEsrGc759ZFrDsMuNo59+XhBgm8ABhwCHAbMA34rnNuctjz1gbr/jf89UrSIglg//7uhXZ48V0dOWBFmLFjDyy2Q48LRv9UzyNYWCtnC+Bbs6tfrGb/4/vZ//h+apfXdvtzK70wnaKTiig+tZjiU4vJPypf42aL9GKk+1jPAJY558YNYx+zgTeABc65NWb2RfxV67cDVznn5oc990XgEefcN8P3MW3aNFdWVtb5eMmSJSxZsmSoIYlINDnnW7N7auXesMHPTtmb8vKeW7oPOcR3PUlSS5cuZenSpQCsWLFis3NueiyOq5wtAK37Wql6qop9j+9j/xP7aXqrqdv2tLw0ik7sKrQLjikgLUuFtqSugebswbZYb6L7T4rpQDnwa+fctUMN1swuAZaGEr2Z5QKVwK3Ae51zi8OeuxK43zl3R/g+1PohkqScgx07em7l3rgRmpt7f+3kyT23dM+c6S+yTBIj2GKtnC0D0rS1iepnq6l6poqqZ6poXNe9f3ZabhqFxxdSfGoxRScWUXBsARkFGpJTUtNwpjSP9DhwP12JugPY5ZxbP4z4QnF0xuKcazSzDfj/BEoinlsAbB/m8UQkUZj50UQmTYJTT+2+raMDtm7tuZX7rbdg2za/PPVU99elpfkxug89FObNg8MP9zNczp3rp5hPHcrZMiA5B+WQc2UO5VeWA360kfBCu+G1Bqr+WUXVP6v8C9L8lOyFJxRSdEIRhScUkjM9B9MFypLiBltY/8Y592z4CjObbWaHOOc2DiOONUCRmZU753YF69rwyXiSmeUGiTsTmAo8M4xjiUiyCBXIU6fCmWd239bWBps3H9jKvWEDVFR0Lf/4R9drzHxr9hFHdBXbhx/uC/DROSGOcrYMSfaEbMa/bzzj3zcegJbKFqqfq6bq2SpqXqihblUddav9sv1u/3dT1oSsboV2wVEFpGWr+4iklsH+T3IO8GzEul3AU8DCoQbhnHvTzP4GXAr82MyKgSzgPuBc4F3Ao8CZwMPOuW1DPZaIjBIZGX487YMPhnPO6b4tNBPlunXw2mvw6qt+CU2M89Zb8Mc/dj0/K8u3Zi9aBMcc45fDD/fTyyc35WyJiqyyLMouKaPsEt8vvr2hndrltdS8WEP1C9XUvFBDy84W9vx+D3t+70cHsiyj4KgCCo4toPDYQgqOLSD3kFy1asuoNqA+1mb2afxFKb1dnv+4c+6cXrYNLBCzccD3geX4Fo6fOudeN7MS/GQDFcAE4AbnXF3k69VfT0T61dLii+u1a32hHbqtqDjwuTk5sGCBL7IXL/bdVA46aETCinYfa+VsiTXnHI0bG6l5oavQrn+tvnsPfyCjOIOCY8KK7WMKyJ6YPNdDiECURgUxs/nAh4EfRGxqcc7tGF6Iw6ckLSJDVlsLr7wCy5Z1LRt76CkxaxacfjqccQacfTYUF0fl8CNx8aJytsSQALFRAAAgAElEQVRba1UrtctrqX25lpqXa6h9qZaWnS0HPC/7oOxurdoFRxeQUTgqu2bJKBHNcawznXOtPayf6px7ZxgxDpuStIhE1f79sHy5L7JfeAGefdYX4CFZWb64ft/74KKLID9/yIcawVFBlLMlYTjnaN7WTO3LtdQuC4rtZbW017Z3f6LBmNljyD8qn4KjC/ztwgIyilRsS2IYcmFtZtcDf3XOvR48vgl/1Xe3pwGnO+dOi064Q6MkLSIjqq0NVqzwI5D8/e++0O4IZrIrLYXrr/fLEC6CjFZhrZwtycZ1OBrWNXS1ai+rpW51Ha71wNokd1ZuZ7FdcJQvuDNLkv46CElCwxlubwYQ/lvnUUAesJOunlOG70cnIjJ6ZWT4vtaLF8OXvww7d8LDD8O998LLL8MNN8C//gV/+pMfzSQ+lLMlqViakTcnj7w5eUz4sP9YdjR3UPdqHXUr66hdUUvdyjrq1tTRuKGRxg2NVD5Y2fn6nBk5Xa3awW3WuKx4vR2RQXcFOQpY65xriVg/N9RCEi9q/RCRuHAOHngArrrKP16zxg/jNwgj2BVEOVtGhY6WDupfq+8stmtX1lL/Sj0dTR0HPDd7ana3YrvgqAKyylVsS/REc4KYI5xzKyN2Phs4GIhrkhYRibmNG+ErX/Gt1ADjx8P06XENKYJytowKaVlpFCwsoGBhARM/OhGAjrYOGt5o6GzVrl3hu5E0v9NM8zvN7PnDns7XZ03KIn9BPvkL88lf4Pts58zIwdI09J9E14AKazM7Gf/z4VnBFLnhxuKvOv9zlGMTEUlsF1/sh+wDOPdcuOsuKOhthLvYUc6WVJCWkUb+EfnkH5EPV/t1rj3os72ylroVdf52ZR0t21vYt30f+/66r/P16YXp5B/ZVWznL8wnb24eaVma1EaGbqAt1lXAvfj+eydGbGsBfhbNoEREksLVV/sLFgF+/nOYODGu4YRRzpaUZOlG3tw88ubmwQf8OtfhaHyr0c8Uuaquc9bIlp1+Nsnq56q7Xp9p5M3L61Zs5x+Zr+H/ZMAG9Elxzr1qZqcBlznnlJBFRDo6fNePkMbG+MUSQTlbpIulGWNmjWHMrDGMv6zrO9u8s/mAYrtxQ2PnVO3hcg7O6exCEiq4syZmaRZJOcCA/wRzzlXTSyuHmZ3knHs+alGJiCSqigr43e/gV7+C14NuytddBzNnxjOqAyhni/Qte0I22edkU3pOaee6tto26tfUU7uqtrPYrl9bT9NbTTS91cSeh7v6bWeOzzyg33buIblYuortVNZnYW1mjwD9zXpgwDxgfD/PExFJLu3tfgr0ZcvgpZfgiSdgw4au7ZMnw9e/DtdcE78YwyhniwxPRkEGRScWUXRiUee6jpYOGt5soG5VnS+4gxbt1t2t7P/Hfvb/Y3/nc9Py0sifH9Fv+/A80nMih5OX0aq/Fut1+L562+gaAzWSAe+JZlAiIjFVV+dbojdtgjfe8C3RoaW+vvtzi4rg3e/2w+ude66fgTFxKGeLRFlaVlAsz8/vHGvbOUdTRVNnq3ao6G7Z1kLNizXUvFjTtYN0yJsT0W97gSa3Ga36K6y/DdQ655r7epKZ6epyEUlMbW2waxfs2AHbt8OWLb6IDi2bNsHevb2/fsoUOPZYOOYYOOUUfzuE2RVjRDlbJAbMjNwZueTOyKXskrLO9S2VLZ0t2qGiu2FdA/Vr66lfW8+uX+/qfG72tOyufttBsZ09JVv9tpNcn/87OOf29LU9zMXoKnMRiYetW+HFF+Gtt2DPHqis9MvOnb6Y3rXLT+LSl+xsP/70tGlw2GEwd65f5syBceNi8jaiQTlbJL6yyrIYe9ZYxp41tnNde0M79a/Wd3UjWVVH/Zp6mjc307y5mb1/6vrDPmNsRrc+2/kL8smdnUtahoYATBbDanYxs6nAp4BPoyQtIiNtzx545RW/vPwyvPCCb4HuixmUl/uh8CZN8v2iZ8zwhXRoKS+P5zTkMaOcLRJ76WPSKVxcSOHiws51HW0dNK5v7GzVDl0s2bavjaonq6h6sqrzuWk5aeQdEXQlCXUnmZ9P+hj1205Egy6szSwduAC4FjgL2IXv0yciEh319f4iwTfe6Cqk16zxXTkiFRXBccfB4Yf74e/GjYOyMpgwwRfT5eWQmbp9GZWzRRJPWkZa53jb5VeVA77fdvPW5m7dSOpW1dFU0UTtslpql9WG7QDGzB7TNW370b51W+Ntx9+A/wXMbAbwMeAjQCmwGfgP4O/A8SMSnYiMXq2tvrV5/Xo/8sa6dV33t27t+TV5eTB/Phx5JBx1FJxwgu+ukQKtzYOlnC2SXMyMnCk55EzJYdx/dHVBa61qPbDf9usNNLzhl9337Q52ALmzcrsK7aPyKTiqgIwiFdux1O/ZNrPL8C0dZwAbgTuB/wVudM49FjzthRGLUESSU3u7b2EOXSAYug3d37LFT7LSk8xMOOQQmD27q5CeP9+PFa0iuk/K2SKjS2ZxJiWnlVByWknnuvamdurX1lO3so7aFbXUrqil/tV6Gtc30ri+kd0P7O58bu4hud1bto/SiCQjaSB/xjTip8B9AfiYc249gJn1czWQiIxazsH+/b5lOXIJjbqxebNvle6NGRx0EMya5QvoQw/1t7Nn+4sIE3fkjUSnnC0yyqXnpFO4qJDCRWH9tls6qH+tntoVtdStqKN2ZS11r9TRuLGRxo2NVP5fZedzc2bk+CL7aN+qXXB0AZmlKrajod//uZxzjwKPBhe9XGtm04E/AZ3/AmZmzvV32b2IJIWWFj+iRmhUje3bey6gGxr631d5edeFguG3M2b4Yeyys0f63aQc5WyR1JSWlUbBwgIKFhb4TmBAR2sHDa83ULuytrPgrnuljqZNTTRtaqLyoa5iO3tadmeRHSq4s8Yn1Dj9ScEGm1vNLAO4CLgKWAv8HHi3c25p9MMbuEWLFrnly5fHMwSRxOUcVFf7QjlUMPd2u2/fwPZZWOhbnHtaQkPXjRkzom9rtDCzFc65RSO0b+VsEenU0RbMJLki6Eay0o9I0tFwYNe87IOyfZG9qIDCYwspOLaAzGK1bPeVswddWEfs+BDgM8C1zrm4/g+qJC0pp6Wla8zm3bsPvA3dD7U+NzUNbL/p6V3D04VG1pgypatonjLFD1lXWNj/vmRARrKwjjiOcraIHMC1OxrW+ZbtUMFdt6qO9rr2A56bOzu3c/jAwsWF5M3PIy0zta596StnD6sTo3NuI/BZM1sxnP2ICL4/8t69vRfIkeuqqwe3//z87sVy+P3w23HjdIHgKKWcLSI9sXTrHP6PD/h1rsPRuKGR2hW11LxcQ+3LvnW7cV0jjesa2XWvn0UyLSeN/KPyKVzsW7QLFxeSMz0nZWeQjMrVQc65e6OxHzPLApYBn3XOPW1mucB3gFqgBPiSc64mGscSGVEtLb5I3rPnwNue1u3dO/hCOS3Nj9c8fvyBt6H7ofGcJ0zwhbUI0cvZoLwtMlpZmjFm9hjGzB5D+ZV+rO2Olg7q1tRR+1ItNS/VUPNSDY3rG6l5oYaaF7q+5pnjMyk8tpDCEwopOqmIgmMKSM9JjQltEu2y+y8C08Me3w085px70MzOAe4B3h+PwCSFNTUNvkiure1/v5HS0mDs2J4L5ciCefx4KClRy7IkAuVtkRSRlpXWORrJ5E9NBqB1fyu1L/tW7ZqXaqh9qZbW3a3sfXQvex/107VbllFwTAHFJxdTdHIRhScUjtq+2glTWJvZCcAOYH/weBJwJX76XYAngIfNbLpzriIuQUpyc86PZBFeAPdWMIff1tcP/ljp6VBa6rtVDPS2uFiFsiQV5W0RySzJZOzZYxl79ljAzyDZtKmJmn/XUP2vaqqfq6Z+bT01/6qh5l818C3AIO+IPIpOKqL4lGKKzygmq2x0jECSEIW1meUBlzrn/tPMbglWnwbscc41ADjn2sxsE3AqUBGXQCVxOAc1Nb0Xx72tG+gFfOEyMwdXIJeW+mm2U7R/maQG5W0R6YmZkTszl9yZuZ1dSFr3t1LzQg1Vz1VR/Vw1tctqqV9TT/2aerbfvR2A/AX5lJxZQslZJRSdVET6mOTsOpIQhTXwJfzfMOEmA3sj1tUCk3raQWVlJYsWdV2guWTJEpYsWRLNGGWkdHT4yUYGUiSH329rG/yxcnK6it+BFsr5+SqSJeqWLl3K0qWdI96N6+u5CWpYeVs5WyR1ZJZkUnp+KaXnlwLQ3thO7bJaqp+rpurpKqqeq+qctn3LHVuwLKPoxCJKziph7DljyV+QH/eLIQeas4c13F40mNnZQLZz7pHgcQVwNbAIuMw5tzjsuSuB+51zd0TuR0M3JYi2Nj8O8kBbkPfs8UV1b1Nb9yU/v+diuK91GldZElCshtuLlmjkbeVsEQlpb2yn+l/V7H9iP/uf2E/dyjoIK0+zp2Yz7oJxlF5YSvGpxXEf3m/EhtuLkuuBhWF/iZTgZwm7M7gfrgDYHrvQUlxz88CL49DtYEe2CCku7r9IjryvWftE4kV5W0SiJj03nbFnjmXsmb6fdsueFqqeqmL/P/az99G9NL/TzLYfbWPbj7aRXpRO6XmllF1WRul5paRlJ9a1SYnQYl0GhFdILwL/CbwEvA6UOecazSwTqANmOue2Re5HrR8D4JxvHQ6Ng7xrV9f9ntYNpUgOjWwx0Bbk0lL//IxE+BtPJD6SsMV62HlbOVtEBsJ1OGqX1bLnkT3s+dMeGl5r6NyWUZJB2fvKKP9AOUUnFsWsu0hCt1g75yrDH5tZO1DpnHvHzB4F3gU8CpwJPNxTUZ3SnPOtxdu3+2Xnzt6L5t27B9cvOSNj4MVx6FYjW4iMesrbIhIrlmadszzO/MZMGt9qpPIPlez6zS7qX6lnxz072HHPDnJm5jD5U5OZcM2EuA7lF/fCuh+fBL5jZkcCE4DUubIlNOpFqGAOX7Zt67q/Y4efjGSgioq6xkMeP95PXd3bYxXJIjJ4qZu3RWTE5R6cy9TrpzL1+qnUvVrHrt/sYtd9u2h6u4m3vvAWm27ZxIQPT2DKF6eQOz035vHFvStItCTdz4qtrbB1K1RUwObNXUtFBWzZ4ovmhob+9uIVF8OkSX6ZMKGrOI4smsePV79kkQSUbF1BoiHpcraIJCzX7tj76F62/mArVf+sAsAyjYlLJjLtpmlkT4xu7ZPQXUFGtZoa2LAB1q2D9eth48auQnr79v5HwhgzBiZP7iqaQ0v4uokTNdKFiIiIpCxLN8ZdOI5xF46jbm0dW769hV337WL7j7ez695dzPzmTCZ9chKWNvJ9sFVYR0NlJaxaBa++2lVEr1/vu2n0xswXyNOnw7RpfgndnzrVbyso0PjJIiIiIgOUf3g+c349h6lfnsrbN7zN3j/vZcOnN7D7wd3MfXBu1FuvI6mwHqzdu+H552HlSl9Mr17tW597kp0Ns2bBoYf6ZdYsmDHDF88HHQRZo2P6ThEREZFEkjcvjyMeOYLK31ey4VMbqH6umhWLVjD/r/PJPzJ/xI6rwro/+/bB3/4GTz8Nzz3nW6Qj5efDkUf65bDDfBE9ezZMmQLpyTklp4iIiEiyK7ukjKITi3jtsteofq6aV858hQXPLCBvbt6IHE+FdU+qquD+++Hhh+GZZ6C9vWvbmDFw/PGweDEsWOCXgw/W6BkiIiIiCSirPIsjHz+StRevZd/f9rH2orUcvfJoMvKjXwarsA63eTN84xtw331dI3Kkp8O73gVnnw2nnAJHHQWZ8RsfUUREREQGJy07jXkPz2PlsSupX1tPxdcqOOSOQ6J+HBXW4CdNue02+Pa3/TTe4IvpD38Yzj/fzwwoIiIiIkkrPTedw351GCsWrWD7j7cz9YtTySqP7vVu6r/Q0ADnngu33uqL6iuvhDfegCeegA9+UEW1iIiIyChRcHQBpe8ppaOpg93/tzvq+0/twto53yr9xBN+MpWnnvLdQA47LN6RiYiIiMgIGH/FeAD2Pro36vtO7a4gjz0GDz0EhYW+qJ4zJ94RiYiIiMgIKjq5CIDaFbU457AozhmS2i3Wd9/tb2+6SUW1iIiISArInpJNemE6bXvbaN3bGtV9p25h3d4OTz7p73/oQ/GNRURERERiwszImugvWmzdpcI6OjZuhMZGP334hAnxjkZEREREYiRrvC+sW3a3RHW/qVtY79jhb6dNi28cIiIiIhJT6Xl+ZuyOho6o7jd1C+u9wZWgpaXxjUNEREREYsqy/QWLHc0qrKOjrs7f5ufHNw4RERERiam0bF8Cq7COlrY2f6vpyUVERERSigrraFNhLSIiIpKSLCMYu7o9uvtVYZ2R2nPkiIiIiKScoK52zkV1tyqs09PjG4eIiIiIxFZossXo1tUpXFh3BH1qVFiLiIiIpJTOacxHa2FtZpeY2SYz22tm3zezjGD9ODO7x8y+bWZ3mpk6RYuIxJlytogktdHcYm1mU4GLgEuBzwMfAz4XbP4d8Avn3JeA9cCtcQlSREQA5WwRGQWCCjjafawT5cq9acA1zrk2YIWZHQ6cbmbPA7Odcy8Fz/szsM7MbnfO1Q/riFE+kSIiKST2OVtEJIo6u4JEd7S9xGixds49FyTokG3AO8AZwKaw520L7i6K2sFDJ1ZERAYkrjlbRCQaRnNXkB4cA/wEmAzsjdhWC0yKeUQiItIb5WwRSS4jNNxeonQF6WRmhwCVzrlXzMwBjRFPyQZaI19XWVnJokVdjSJLlixhyZIlvR9IXUFEJI6WLl3K0qVLQw/HxTOW4YhZzhYRiaZBtlgPNGcnVGFtZunAx4EvB6u2AYeGbTcgD9ge+dqysjKWL18+lIMOJVQRkWEJLyTNbE+cwxmSuORsEZEoGOxwewPN2YnWFeR64A7nXHPw+HFgdtj26UALoGwsIhJ/ytkikpxCFfBo7WNtZjcCq4BcM5tpZtcAVcCO4IpzgHOAHzrnWoZ9QHUFEREZspjnbBGRaAo1WHeMwj7WZnYzcFvE6jedc78ws8uBW8ysAsgHborywaO6OxGR0S6uOVtEJApGaubFhCisnXO3A7f3sm0z8NHYRiQiIr1RzhaRpJdiw+2NPHUFEREREUlNI9RhIXUL6xB1BRERERFJSdEex1qFtYiIiIikFrVYi4iIiIgkrtQtrNXHWkRERCS16eLFKFMfaxEREZGUYiNU/6mwFhEREZHUpBbrKFFXEBEREZHUpIsXR4i6goiIiIikJrVYi4iIiIgMg1qsoyw/HyZOhIKCeEciIiIiIjGUUZhB1sQs0vKiWwpnRHVvyeRTn/KLiIiIiKSUaTdNY9pN06K+39RtsRYRERERiSIV1iIiIiIiUaDCWkREREQkClRYi4iIiIhEQcoX1kuXLo13CIOmmGMjGWOG5IxbMctA6bwPjM7TwOg8DYzO08CpsE7CD4tijo1kjBmSM27FLAOl8z4wOk8Do/M0MDpPA5fyhbWIiIiISDSYc1GeyzFOzKwS2DyEl44D9kQ5nJGmmGMjGWOG5Iw71WOe5pwri9K+kkKK5ex40HkaGJ2ngdF56q7XnD1qCmsRERERkXhSVxARERERkShQYS0iIiIiEgUqrEVEREREoiAj3gHEi5nlAt8BaoES4EvOuZr4RgVmdgnwPaAQ+A3wBedcm5mNA74BVAGZ+Hhbg9f0ui3GsWcBy4DPOuee7uscJ9L5N7PjgJOAjcDzQH2ixm1mc4BPB7HOApY651Yn4rk2s7OA24HLnXMV/cWSCO+hl5h7/E4G2xL+ezlaxPu7l0zM7HbgpuDhGufckfGMJ5EMNi+lqp7OU7Ben63+OOdScgF+if/AAJwDPJAAMU0F7gWOBj6EL/CuD7Y9BSwO7n8C+GbY63rdFuP4bwKqgdP6O8eJcv6BjwC3D/SzEe+4geXA5LDPyxuJGDNQDlwMOGD6cM9tLN5DTzH39Z0Mtif893K0LPH+7iXLAuQBPwHODJZZ8Y4pUZah5KVUXPo4T/psDeT8xTuAuLxpmAQ0A2OCxxnBf5jT4xzXyUBG2OPvAH8BjgO2h62fDNQFH/Jet8U49hOAa4AK4LS+znGinH/gVOAJgtFx+vtsJELcwfEOC+6XAdsTNWZ8V7PwInVIccbyPfQQc4/fyeB+wn8vR8uSCN+9ZFmAzwVLTrxjScRlMHkp3rEm0nkK1umzNYAlVftYnwbscc41ADj/s+4mfKEVN86554JYQrYB7wBn4OMLPW9bcHdRP9tiwszygEudc78IW30avZ/jvrbF0p3AG8APzewxMzsxCeJ+CPiZmRUAHwA+k6gxO+c6IlYNNc6YvYfImPv4TkKCfy9HmdOI/3cv4ZmZAVfguy7tMrP3xzmkhDPIvJSyIs+TPlsDl6qF9WRgb8S6WvxfronkGPzPLn3Fmwjv5UvAtyLWJXTMZjYbWAj81Dn3aeBJ4DFgSh+xxT1u4DqgCd+Xvc4593A/cSVCzCFDjTOR3kPoOwnJE/NooPM5AM47Dj+Zxw+B+8zs3DiHlej02RoAfbYGLlUvXnRAY8S6bCBhLiwys0OASufcK2bWV7xxfS9mdjaw3Dm3O2JTf3HF+/zPA/Y559YEj38EfJXEjzsXeAD/n8FdZraJxI85ZDhxxv09hH8ng1UJ+70chXQ+B8E5tx+4OWhl/BzwtziHlMj02RoEfbb6l6qF9Tb8lb/hCvD9VePOzNKBjwNfDlZtAw4N2274fpzb+9kWC9cDC/1hAX9e/4TvZtHbOXZ9bIuVDMI+/865RjPbAKT3EVsixP0b4ArnXFXwB9cD+OSWyDGH9PW96yvOuL+HHr6TkNjfy9EmoXN2AvsB8HS8g0hw+mwNjT5bvUjVriBPAZOCIXYws0z81f/PxDWqLtcDdzjnmoPHjwOzw7ZPB1rwI0T0tS0WrgQWhC3bgY/hr7Lu7RwnwvlfAxSZWXnYujaCiwETMe5g+LYjnXNVwarbgRpgc6LGHKGvWIa6LVYiv5OQ2N/L0SYRPgPJqANYGe8gEpw+W0Ojz1YvUrKwds7tBB4F3hWsOhN4OOwCo7gxsxuBVUCumc00s2vw4+DuMLPDg6edA/zQOdfinFve27ZYxOucq3TObQ0tQDv+5/J36OUcJ8L5d869if8J61IAMysGsoD7EjjufUCTmU0OW7cX/0dCwsVsXT9jGPT9vRvqtpGOOVh3wHfSzA7p67sX7+/laJMA372kYGaTzOxKM0sLPstfoGvMYWFweSkO4SWMyPOkz9bAmfNDqKQcMyvBD51VAUwAbnDO1cU5ppuB2yJWv+mcm2Nm04Bb8PHmAze5rkkqet0Wa2ZWAVzt/AQxvZ7jRDj/QQvw9/GtiFPxFzK+nshxm9mR+AsYV+DHGn3WOfdMosUcNmrJ3cCtwI+cc5VDjTMW76GnmPHdP3r8TgavSYrv5WgQ7+9eMjCzefgGgxr8ZFd3BY0IwtDyUirqJReOR5+tAUnZwlpEREREJJpSsiuIiIiIiEi0qbAWEREREYkCFdYiIiIiIlGgwlpEREREJApUWIuIiIiIRIEKaxERERGRKFBhLSIiIiISBSqsRSShmFmZmT1iZpvN7MvxjkdERPqmvN1FE8SISEIxs7OBp4BC4E3gUOfcvvhGJSIivVHe7qLCWkQSlpk9C5yhqcBFRJJDqudtdQURkYRkZlOAf6RqchYRSTbK2yqsZQSY2f9nZj+Odxwh0YjHzM40sxfN7LQohTWqmNkdUd6fAe8DvhXN/YpId8rXw2dm083sp2Z2T5T2l25m15rZ5gE893wze2mo7zWauVt521NhnWTMbI6Z/d7MfmZmr5iZM7MPxDuuCH8FHoh3EGGiEc8W4MgoxJKQzOw8M9tjZjvM7HgzmxJ8zpyZ/cjMioPnHWtma83se2aWG6wrB3ZGOaQrgZ8759qC/YskHeXrIUnGfL07uM2O0v7SgCpg6gCeuxyY09tGM5tnZg8Gn72bzWxC2LZo527lbVRYJ6M/Al93zn0MWAjcG9pgZvPN7NR4BBV+bOfcWufc8/GIoyeDiae3c+icWwfsiWZc0fz3Gu6+nHN/Be4G2pxzLzrntgDXAi3AS865quB5L+MT+Redc43By/8DeHRYbyCMmX0F+BrwhJm9CcyO1r5FYkz5epCSMV875xqA7dE6lnOuFVgxwOfuwhfhvW1/DVgaPLzDORdeSEctdytvd1FhnUTMbCxwKFAP4JzrAG4BnJkV4ZO2xSGuuB07mgbwPjpieKx47OtB4CAzOxrAObcXeBJ4T9ixioHtwWcv5DDn3JvDPHYn59xtzrlZzrlFzrnDnHPPRmvfIrGifD2yEjBfR3skiMHE399z24PbyH7PUcvdyttdVFgnkWDomnXAY2Z2XLBuM/Ayvl/TDGCJmX0EwMxuNbOvBD8DdfajMrNFZnanmX3OzBrMbJ2Z/aeZXWRmy8zsYjP7g5nVmNnHw17X4/4ijn2tmf0qvK9Z8HPoD83sm2b2lJkdG6zv83hhrzcz+7iZvWZmZ5nZGjPbZWZXD+AYs8Lj6eeYB5zDHpQEP+3Wmtm1Eef0a8G5+X2QiAd6zg44lpldamYvm9nVZrbBzL7Vx/56+rfvMZ6+BC0brwGXha0uA863oNsHcAHwSFicuUBD2OPDzOxbZnaHmT1jvkvJ5Wa23MwuNLMHzKzKzD5rZkcH73GvmZ3XX3wiyUT5OnXydU96e499HcfMZpjZL8y3/n53IMfp4bh/NbP7zOzEfp6n3D1SnHNakmgB5gGb8X+h/gIYH7atAjgtuH848HZwvwz/13Qp/i/ud4BFwbYfAy8E97OAfcCP8GNRvh/Y1tf+Io+N72N2D/CrYH0OvlgrDh5fClQCBX0dr4f3PSE45o3Bc28HmoBp/RwjMp4+jxl+DnuIoQL4HjAmeN07wfpC4Ldhz/sXPikO6Jz1cqyx+J/37gKOBc7sa38R//Y9xjPAz9fNwFvB/YOBbwLNwCVhnxcLexcfVPgAAAXxSURBVP4FwMnB/THAMmBM8PjZ4DwUAtXArUAucHnw+H34P+6/QPAZ1KJlNC0oX6dEvg62f22A57Gvf+tlwJyw17gBfs5C/57HA1/tYftpwXEywtYpd4/QohbrJON8q+IR+ETzIeA1M+vpIo0NwAfNLBM4OViXj/8iTwFC/WPfwH9hcM61ALXAQ865GvwXbUI/+4uMr5nuF0OcDzS6oI8u8Hv8F/KSfo4Xud/QPu8Pnns7vv/vWf0co1s8gzlmL/7sfH+6ZcDEsPdYGrQofQ5/TmGA56yX97sPn8D+5Jx72Tn3xCD211s8A/EgMNN8d5D3Aj8F/g68z8zygVoXZN7AicALwf0L8P95hVpBLgZuC85zFfBP5/tlLwcKnXP/5/zP46sY3L+BSFJQvk6NfN2DXt9jH8c5EyhyzoXiWT7IYx4DXOWc+/oAn6/cPUIy4h2ADJyZpQNTnHMVwP8zs5/jr6BeCiwOf65zrtnMJgNfAv43tAvn3G4zWw68C/8X9Qx8MdX50rD7bQTdhXrbXy+hhu9jFpAZFleHmVUAB/V1vP4455rM7G2gGBg3iGMM+ZgR2uj6/kwFNjnn7op80kDOmZl9EF/Ahsx2/idjFx7rIP4Neo2nP865DWa2Ct8dpMg597aZPYhvRboM+EtY3KHPRqj/3jR8C1NoX3vDdx12v53uOlC3NBlllK+7vb9UyNfhej2PffzbzMf/ETFUZwCzzCzfOVfXy3M6gveg3D2CdEKSSyZh/V+dc6/gf2o7LPKJZnYGcJ1z7nbnR3gI9yngGPPDPq0Bvt3fgfvZX1/ewbeAhg9DZPi+h8OVBawf4WMMxA6690PGzBYP4pw9AiwIW3q8unwQ++sxnoG/HX4LfAR4Oyy+dOCz+J9NQxbj+4uGbANOMrO8sOP22c9PZBRTvu4ulfJ1r++xj+PUBK/JYmjuwLf2f7+X7dWu66Jz5e4RpMI6+XzGzCaFPZ4MPBbcb8FfrHEYfminIjPLNrMzg+3jzWwcvq/fz/E/47xA8NNiII2Iv9LNzPrZX+SxLWwffwTqgCuCfRUHxwu1fPZ2vN5MDJ4zOXjtYwM4Rng8/R0z/H1E6qn1wvCtUHnAo2Z2tvmLZPr7N+g8FjDROfdm2NIadoj0sPt97S887h7jMbMsM1sadPPoy//hf4J+CMA5Vwv8DfiX6z4ayLvx3URCQuf7AfNjYX8aCI1lGvlvEPnvnNQjFIj0QvmalMnXAz2PPR4Hn0uz8H98gf91AjObNIDcbfjW4w8Cl5nZ+4LXfsHMbsa31j8Z9nzl7pE0Eh23tYzMgr8gwgH7gd8Av8L/lFQSbP8qsBW4EP/TzgZ8C8H5+JaOX+GT1CP4Cylagv014L9o5wKt+AtkpgA3BNuv6Wt/Ecf+LPAcvt/avGDbInxL5/fwfQ0PD9b3erxe3r8LXv9ZfNeEI8K29XaM6eHx9HfM8HMYcewL8Rfx/Sg4F6HXfSTYfiqwNvi3+V5wngd6zi7s4b1eGhzvAWBysK6vf9Nu++olnoPwExncP4DP2g8iHl8OnB6x7o4eXncK8DqwF7g5WHce/rP2Q2ASwZBjwCfw//H+NHiv74n3d0yLlmgtKF+nUr6eCjyPb3Vf0M977Ovf+nz8L4X/xl9Ivh64LnjvPeZu/K8irfiW6rHAP/FF/ceAj+PH8/4evhtK6DXK3SO4WHDyJEWY2XjgE865W4PHafhWlI87526Oa3D9MDMHzHC+z6IMkZl90Tk3pKGcRCR2lK8lnHJ3clBXkNTzaXy3gLFh684g7KK0RBT205N+dhoiM0s3s6uAh+Mdi4gMiPK1KHcnGRXWqecn+L5Yb5jZFuBxYKtz7sX4htW74GKOzwYPrzY/7JsMXj7wR+fc2/0+U0QSgfK1gHJ3UlFXEBERERGRKFCLtYiIiIhIFKiwFhERERGJAhXWIiIiIiJRoMJaRERERCQKVFiLiIiIiESBCmsRERERkShQYS0iIiIiEgUqrEVEREREouD/b7eOBQAAAAAG+VtPYmdRJNYAADAIOluU8zifFPEAAAAASUVORK5CYII=\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/hayabusa-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/hayabusa-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/hayabusa-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 }