{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 41 - Hera - Saturn (Concept)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from AMAT.planet import Planet\n", "from AMAT.vehicle import Vehicle" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook simulates the atmospheric entry of the Saturn entry probe which was part of the Hera mission concept study.\n", "http://dx.doi.org/10.1016/j.pss.2015.06.020" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"SATURN\") \n", "planet.h_skip = 1000.0E3\n", "planet.h_trap = -100.0E3\n", "planet.loadAtmosphereModel('../atmdata/Saturn/saturn-nominal.dat', 0 , 1 , 2, 3, heightInKmFlag=True)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle1=Vehicle('Hera-Saturn-Probe-steep', 220, 269, 0.0, np.pi*1.0**2.0*0.25, 0.0, 0.18, planet)\n", "vehicle2=Vehicle('Hera-Saturn-Probe-shallow', 220, 269, 0.0, np.pi*1.0**2.0*0.25, 0.0, 0.18, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle1.setInitialState(1000.0,0.0,0.0,26.3,0.0,-22.0,0.0,0.0)\n", "vehicle2.setInitialState(1000.0,0.0,0.0,26.3,0.0,-9.0,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle1.setSolverParams(1E-6)\n", "vehicle2.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle1.propogateEntry (100*60.0,0.1,0.0)\n", "vehicle2.propogateEntry (100*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": "iVBORw0KGgoAAAANSUhEUgAAAucAAAHtCAYAAAC6dNAzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXxU1f3/8ddJQsIWCAlbEnYIsskaluICVqm4IVi+rrVYq9Stdan9+W21ta2tVdGvtipqXOpSwQ2tiIpbFawLCgqCQABJWAKBBCQsgazn98edmUz2bSZ3ZvJ+Ph73cWfuvXPncydw8smZcz/HWGsRERERERH3RbkdgIiIiIiIOJSci4iIiIiECCXnIiIiIiIhQsm5iIiIiEiIUHIuIiIiIhIilJyLiIiIiISIGLcDaAldu3a1/fr1a/DxeXl5dOvWLXgBuSDSrinSrgd0TeHAjetZtWpVvrU2cj7EBmhsm+0Vaf/eAkGfSc30udRMn0t1jf1MAtFmt4rkvF+/fqxcubLBx6enpzfq+HAQadcUadcDuqZw4Mb1GGO2tegbhoDGttlekfbvLRD0mdRMn0vN9LlU19jPJBBttoa1iIiIiIiECCXnIiIiIiIhQsl5DebOnet2CAEXadcUadcDuqZwEGnXE2n086lOn0nN9LnUTJ9LdW58JsZa27JvaMw04C/ABdbabM+2dsA9wCGgC3CLtfZgc/b5S09PtxpDJSLhyBizylqb7uL7q80WEWmgQLTZLdpzbozpAXQEJlTZNR/4r7X2d8DrwGMB2CciIs2gNltEpOW50XMeBZQB/a212caYFCAL6GKtLTTGxAAFwHCguCn7vL07XuqFERG3WGspLS0lOjqaqKjG94eEQM95yLfZ1kJZmfM4plXUIBORUBV2PecA1tryKpumAvnW2kLP/lKcBnxKM/aJiISEcePGERsby6pVq9wOpUnCoc3+5BNo0wZOOqm5ZxIRcV8o3BCaCuyrsu0QkNKMfZXk5eWRnp7uWzIyMgISuIhIfTZt2gTAmjVrGvyajIwMX3sFdA1OZE0Wcm32t98663XrGnYBIiKBFOg2OxS+ALTA0Srb4oASz+Om7vPp1q2biuqLiCu8Q1mOHTvW4NfMnTvXVyHAGJMflMCaLuTabGMafKiISMAFus0OhZ7zHJy79v3FA7uasS8iLFy4kN/85jdcc801dOzYkZtvvtntkESkkaKjo4HGJechLuTabG9y3sK3UIlIECj3CY3k/EMgxVNiC2NMG6APsKwZ+8Lenj17WLBgAfPmzWP+/Pm8+eabvh64J554wuXoRKShvMl5UVGRy5EETMi12U24z1ZEQpByH0eLN2nG+L6ANADW2lxgCXCqZ/tpwCJrbU5T97XAZQRdbm4umZmZFBQUADBlyhQmTJjAe++9x8KFC12OTkQaKtx7zsOhzVbPuUhkUO7jaOk65/HAVZ6nlxljunkeXw2ca4y5FTgT8J+Oqan7wtrw4cNJSEhg7NixvPPOOwDMmDGDN954g+zsbO69914A3njjDebNm8fpp5/Oq6++CsDu3bu5/fbbuf7667nssssoKiriL3/5C+eeey633XYbycnJ3HTTTa5dm0hrEs7Jebi02eo5F4kMyn08rLURv4wbN842BM6NTkFfGmrfvn32/PPPt4D98Y9/bPft22c//PBDO2XKFGuttXv37rUXXHCBtdbab7/91nbp0sUeO3bMXnTRRXbPnj3WWmvHjx9vX3rpJfvOO+/YtLQ0u2fPHrtlyxYbGxtrP//88wbHIiJN06tXLwvYX/3qV016PbDShkA72pJLQ9tsr6efthasbd++US8TafWc75uCvzRGuOc+gWiz1d8QwhITE3nxxRdZsmQJn332GRdeeGGl/StWrODw4cM8/fTTfPbZZ5x88sl8//33fPHFF7z11ls8/fTTjB8/npKSEmJjY0lJSaF79+4MHDiQKVOmqIKNSAsYMmQIAH369HE5ksiVmuqse/RwNw4RaT7lPqFRSjFkOH/whIb169fTtm1bBgwYwFlnncVrr73GpEmTuOGGG3zHlJaWEh8fz2WXXQbA5ZdfTmlpKQUFBcyePZuOHTv6vtr57LPPKp0/MTGRTp06teQlibRKvXv3BiAhIcHlSCJXYqKz7tzZ3ThEwk0IpT2Ach8v9ZyHqPLych5//HHf8x49ejBw4EDi4+MpLi6msLCQMWPGsGTJEv71r3+xZ88e/vGPf1BSUsLkyZO54ooryMnJ4eOPP+bDDz8E4MiRI77zZWdnc+qpp1Z7XxEJrDZt2gBQUlKtnLcEiG4IFYkMyn0cSs5D2F133cXZZ5/Nrbfeyu9+9ztefvllhg8fzt69e/njH/9I3759efzxx/nf//1f0tPTGTRoEO3bt+fBBx9k165dDB06lAULFnD66acDsGvXLubPn8+dd97JTTfdREpKtYn5RCTA9u1zJsXMzs52N5AIduCAs/7+e3fjEJHmU+6jYS0ha8SIEbUOs9myZYvv8YUXXlhtPFafPn1Yvnx5tdelpaVxzTXXBDZQEanT1q1bAdiwYYPLkUSuvDxnvXevu3GISPMo93Go51xEJIhiYpw+kOLiYpcjiVzeUooa1iIikUDJeStw5MgR3nzzTTZv3szatWvdDkekVfGOOVdyHjy+aZJERDzCOfdRct4KdOjQgXnz5pGTk8Pxxx/vdjgirUpsbCygG0KDSZMQiUhV4Zz7qEkTEQkiVWsJPg1rEZFIouRcRCSI1HMefBrWIiKRRMm5iEgQeXvOy8rKXI4kcmlYi4hEEjVpIiJBNG3aNAAmTpzociSRa9AgZ92vn6thiIgEhJJzEZEgiouLA1StJZg0Q6iIRBIl5yFs4cKF/OY3v+Gaa66hY8eO3Hjjjfzxj3/k7LPPbtDrv/vuO84++2yeeeYZiouL+fOf/8w555wT5KhFxJ9uCA0+JecikUO5j5LzkLVnzx4WLFjAvHnzmD9/Pm+++SZt2rRh0qRJHD58uEHnGDhwIB06dMBaS2xsLOPHj+fQoUNBjlxE/H377bcArFy50uVIItfu3c561y534xCR5lHu44hxOwCpWW5uLpmZmRQUFNC5c2emTJlCXl4ebdu2bdR52rVrV+NjEWkZ3h7zcPvlEE68X0po5JBIeFPu41DPeYgaPnw4CQkJjB07lnfeeQeA2bNn+/bffffdDBw4kEcffRSANWvWMHfuXG6//XYuvfTSBr3HvHnzuP/++5kzZw6LFy+mvLycOXPmMGLECA4cOMDNN99Mr169yMzMZP/+/fzwhz8kMzMz8BcrEsG8Y85LS0tdjiRyqZSiSGRQ7uNQcl6FMabWJSMjw3dcRkZGncf6GzduXI3b6xITE8PSpUtJT09n+vTpzJ49m/379wOwe/dufvGLX/DMM8/w2GOPAc4YrTPPPJM//elPvPvuu+z2fs9bi4ULF3L06FFuvPFG7r//fi677DK2b9/Ovffey/79+4mPj+eOO+7g8OHDJCcnEx8fz6xZszjuuOMafA0iUpGcl5eXuxxJ5FIpRZGmM6b2xS/tISOj7mP9jRtX8/b6KPdxqEkLYYmJibz44ossWbKEzz77jAsvvBCA5ORkEhISSE1NpaCgAIC77rqLvn378tRTT1FWVkZRUVGd537llVfo37+/733S09N555136NatG6NHj+aDDz7g66+/ZtCgQSxatIi3336b6dOnB/eCRSKQ9+tY9ZwHj2YIFYkcyn2UnFdjra11mTt3ru+4uXPn1nmsv1WrVtW4vS7r169n69atAJx11lm89tprvP/++xQWFvqOMcb4euMeeeQRPvzwQy6//HI6duzYoOvcs2eP73nXrl19VSUuueQSFi5cyIoVK5g3bx7PPfccW7ZsIS0trcHxi4jDm5yr5zx41HMu0nTW1r74pT3MnVv3sf5Wrap5e32U+zjUpIWo8vJyHn/8cd/zHj16MHDgQNq3b1/j8fPnz2fs2LEUFBRQWFjI4cOH65yR8JxzzuG1117zPc/NzeWss84CYObMmbzzzjskJyczdepUvvvuO+Lj4wN0ZSKti5Lz4NOYc5HIoNzHETLVWowxQ4HrgC1AGpBhrV1tjGkH3AMcAroAt1hrD3peU+u+SHDXXXexdu1aRo0aRXZ2Ni+88AIvvfQSW7du5dtvv2X58uXk5eWxcuVKfvzjHzNnzhyuuuoqhg0bxgsvvMBVV13FN998Q1RUFLNmzWLp0qV89913rF+/njlz5vDVV1/xs5/9jOOOO46bbrqJHj16ANChQwfOP/98Zs2ahTGGn/zkJ8ycOdPlT0MkPCUnJwPQqVMnlyMJrFBqs72/P2v5/S0iYUS5D5jGDLUIJmPMSuBca22OMaYP8I61dqgx5p/AUmvti8aY6cAca+1FntfUus9fenq6VY1hEXHDN998w6hRoxgxYgRr165t9OuNMaustelBCK1ZQqnN3rED+vSB1FTYuTMglyci0iSBaLNDaVjLUMD7/cFRoLMxJgW4GHjDs/19YIYxpl9d+1osYhGResTGxgIROUNoyLTZmiFURCJJKCXnrwBPGGPigZ8AvwSmAvnW2kIAa20pkAVMqWefiEhI8I41j8BJiEKmzfYm5ceONfdMIiLuC6Xk/BrgGPAlcNhauwhIBfZVOe4QkFLPPhGRkLBvn9NM+VcIiBAh02Z7k/Lvv2/umURE3BcyN4QC7YCFOA34A8aYLMDifF3qLw7wfj9c1z6fvLw80tMrhv/MnTu3UllEEZFg8U4d3Zj7ezIyMvwnPesa+KgCImTa7A4dnLWGtYiIGwLdZodScv4v4EJr7QFjjMVp9G/AuaPfXzywC+eXQG37KunWrRu6IVRE3NCU5Nw/GTXG5AclsOYLmTa7AeWNRUSCJtBtdkgMazHGdAVGWWsPeDb9BTgIbANSPOW3MMa0AfoAy4AP69gnIhISvPV5Q6UyViCEWpvtX0JR5eRFJNyFRHIO7AeOGWNS/bbtA74BlgCneradBiyy1uZYa3Nr29dCMYuI1Ku2yTPCXEi12TF+3wHrplARCXchMazFWltujJkJ/MEYswroAfzGWnvQGHM1cI8xZhTQE/AfeFjXPhER10Vich7Kbfbhw5qMSETCW0gk5wDW2jXAL2rY/j1wZS2vqXWfiEgo8NY5BygrKyM6OtrFaAIn1NpsY5wbQg8fhu7dg/EOIiItI1SGtYiIRKQ2bdrQpk0bAIqLi12OJnIlJDjruDh34xARaS4l5yIiQRQVFUWcJ2OMwFlCQ4Z3KItuCBWRcKfkXEQkyLw950rOg8fbY64bQkUk3Ck5FxEJopKSEo4edebe0bCW4CksdNYbN7obh4hIcyk5FxEJorKyMo55unPVcx483uQ8K8vdOEREmkvJuYhIC1FyHjzeYS35oTqfqohIAyk5FxFpIRrWEjzt2jnrffvcjUNEpLmUnIuItBD1nAdPhw7O+vvv3Y1DRKS5lJyLiLSQctX5C5qOHZ31gQPuxiEi0lxKzkVEgqi0tNTtEFqFTp2ctZJzEQl3Ss5FRIKorKyM9p4Zcqy1LkcTuQYMcNb6ckJEwp2ScxGRIOrcuTODBw92O4yIN2OGs+7Sxd04RESaS8m5iIiEveRkZ52b624cIiLNpeRcRCSIysrKfDeCalhL8PTs6ax37nQ3DhGR5lJyLiISRF9++SXffPON22FEvK5dnXV+vsopikh4U3IuIhJEhd555SWo4uKgTRvn8YcfuhuLiEhzBCQ5N8acGYjziIhEmn1+U1aGyrCWSG2zExOd9SefuBuHiEhzxDTmYGPMdOA2oCcViX0U0ANoF9jQRETCX15enmvv3dra7N69Yc8e0CgiEQlnjUrOgeeAvwFfA95qsga4OJBBiYhEivz8fN9jF3rOW1WbPXIkrFwJGza4HYmISNM1Njnfaq39v6objTEbAxSPiEhE8U/OY2NjW/rtW1WbPXMmPPWUU07RWjDG7YhERBqvsWPOrzbGzK5h+4xABCMiEml2+tX2i4uLa+m3b1Vt9hlnOOuyMqcHXUQkHDU2OT8BeMIYU+a3lAOPBCE2EZGw99vf/pYunmkr27Zt29Jv36ra7JgYGDvWebwxIr8bEJHWoLHDWv4MzAVWAmWebQa4NFABGWMmAScCW4D/AkeAe4BDQBfgFmvtQc+x7WrbJyISCsaPH08bT40/F3rOW12b/ZOfwFdfwdKlcGnArlJEpOU0Njl/G3jdWlvkv9EY83gggjHG/AwYaK29zW/bP4Gl1toXPZUHHgMu8uyeX8c+EZGQcOzYMcCVnvNW12affTbcdBO8/TaUljq96SIi4aSxzdYLwG+NMf/x22aA84DrmxOIMWYKcAkwzW9bCk5VgWs9m94HFhlj+gHFte2z1mY3JxYRkUD4/PPPefXVV30TEbnQc97q2uy0NOjWDfLy4MEH4cYbA3VmEZGW0djk/PdAGnAZ4K0JZnBq6DaroQf+D/gUeNAYMwi4A+gL5FtrCwGstaXGmCxgClBSx77sZsYiItJsS5cuZd68eb7nLiTnrbLNHjJEybmIhK/G3hA6D+htre1nre3vWfrh9MI0mTHmOGAM8Li19jrgA2Ap0BvYV+XwQ0AKkFrHvkry8vJIT0/3LRkZGc0JV0SkQVasWOF73LlzZ6KiGtbkZmRk+NoroGszQmiVbfYddzjrrCz47rtGvVREpNEC2GYDje8532GtPVTD9p7NjGM4sN9a653X7SHgdpyenqNVjo3D6YGhnn0+3bp1Y6XqaolICzp27BjLly/3PU/0zi3fAHPnzmXu3LkAGGPy6zm8Lq2yzZ4yBZKSYN8++PnP4aOPmnwqEZF6BbDNBhrfc35N1Q3GmP8BHmhmHDH4/aFgrT0KbAaice7o9xcP7AJy6tgnIuKqZcuWUVhYyODBgwFISkpyI4xW22b/4Q/OetkyyM4OxjuIiARHY5PzRGPM7wGMMe2NMU8ATwBfNjOOb4DOxpgefttKcRrtFE/5LYwxbYA+wDLgwzr2iYi46tVXXwVgzJgxQON6zgOo1bbZv/wldO7sPJ45MxjvICISHI1Kzq21ZwHvGWMeBFYDg4FRQE0z0DXmvBtxSn7NBjDGJACxwPPAEuBUz6GnAYustTnW2tza9jUnFhGR5jpy5AgLFy4EYNSoUYA7Peetuc02BubPdx6vWaOhLSISPuocc26MObmWXdtwvsK8F+fu/JlAc++J/ynwd2NMLE5vykWeu/mvBu4xxozCGSc51+81de0TEXFFUVERV199NZmZmZSVOXP/9OnTJ+jvqza7sosvhiVLYOFCuOwy+Ppr6FJ1YI2ISIgx1tradxrzOTAUOACU13YYkGytbfEaYQ2Vnp5udUOoiLjhiiuu4Mknn+SRRx7hqquuavTrjTGrrLXpDTxWbXYVxcVwwgmwcqUzvGXRImhg0RwRkUZrTJtdm/qqtfwZWGmt3VtPIM36ilREJFKUl5dXKpmY7bkbsV+/fi3x9mqzq4iNhRdfhNGj4d//hvPOc9YiIqGqzv4Da+1b9TXynuNeCVxIIiLhKT8/n+HDh/PMM8/g/VbyO0+h7f79+wf9/dVm12zAALjkEufx66/DNdVq2IiIhA59uSciEgDWWq6//no2btzIs88+C8CBAwfIzs6mbdu2DBw40OUIW7f585365wCPPAJz5kAdozpFRFyj5FxEJAAee+wxFixYQPv27Xn00UcxxrB69WoARo4cSUxMY+d8k0AyBt57D8aPd54/+yxMmwaFhe7GJSJSlZJzEZFmev/997n++usBZxrntLQ0AL766iugota5uKtNG/j004oe9A8+gLQ02LTJ3bhERPw1Ozk3xowxxvzLGPNrY0y/5ockIhI+PvvsM2bOnElxcTE33HADl3gHNwMffvghAJMmTXIrvGpae5sdE+Mk5Zdf7jzftcu5WfTvf4fy2urbiIi0oED0nF8FnGatvQ+Y7J39TUQk0pWWlnLppZdy5MgRfvrTn3Lffff59hUXF/uS82nTprkVYk1afZsdHQ1PPgkvv+zUQj96FG64wSm5+PXXbkcnIq1dIJLz64BxnscvAG0CcE4RkZAXExPDSy+9xC9+8QuefPLJSiUUP/nkE44cOcKwYcNITU11Mcpq1GZ7zJ4Nzz/vVHDp3h0+/xzGjoWf/hRyc92OTkRaq0Yn58aYM4wxl3geTwDSvNMvW2vLrbUHAxyjiEjIOHDgAM8884zv+dixY3n00Uer3fC5YMECAGbMmNGi8VWlNrt+M2bAX//q3DQK8Nxz0Lcv3Hwz7Nnjbmwi0vo0Kjk3xvwVeBmYBWCt/QL4qTHm4iDEJiISUl577TWGDx/OZZddxmuvvVbrcceOHePll18G4NJLL22p8KpRm91wV1wBX34JI0Y4z4uL4b77oE8fuPFG2LbN3fhEpPVobM/5ZKAX8LHftreBuwIWkYhIiNm0aRPnnnsu5513Hrt27WLSpEmMHDmy1uMXLFhAQUEB48aNY9iwYS0YaTVqsxth3DhYvRoefhg6d3a2FRfDAw84ExnNmgX/+Y/qo4tIcDU2OX/PWnsA8G+aWs000CLSuuzevZurrrqKYcOGsXjxYuLj43n44Yf55JNPap1UqKysjLvvvhuAG2+8sSXDrYna7EaKjnZmEN22DW69Fdq2hVNOcbb/+99w6qkwdCjceSds3+52tCISiRqbnG81xvwOOM4Yc6ExZhFwLXB34EMTEXHX448/zmOPPYa1liuvvJKNGzdyzTXXVLrxs6oFCxawadMm+vXrxwUXXNCC0dZIbXYTde4Mf/kL7NjhlF7cvh3+/GcnWc/MdBL3vn2dxP3JJyE/3+2IRSRSNCo5t9a+APwHiAZ+ChwAfmStfTgIsYmItJjy8nLeffdd3njjDd+2a6+9lgsvvJB169aRkZFBSkpKnec4dOgQ/+///T8Abr/9dtdnBVWb3Xxduzo3ivbsCddfD+3bV97/0UfOePUePWDqVKdeunrURaQ5jG3m4DljTDQw1Vr7QWBCCrz09HS7cuVKt8MQkRC0c+dO/vnPf/Lkk0+ybds2evXqxdatW2nTpvEVBq+77joefvhhJk6cyKefflpnD3tDGWNWWWvTm32iivOpzW6GggJ46imYPx+2bHG2eau8+P86HTUKfvQjOP10OPFEiItr+VhFpOUFos2us1vHGLMY6FjPOboC+UDINvQiIv7y8/N54YUXePHFF/nvf//r296vXz9+/vOfU1RU1OjkfPHixTz88MO0adOGRx55JCCJeWOpzQ6+zp2d6i3XXw/vv+/cPLpkibPvH/+AZcvg7bdhzRpnmTcP2rVzetV/9CM4+WQncY+OdvUyRCSE1fed6zeeYzKBcuBK4GvAv0ujP5AYlOhERALAWsuRI0fo2NHJW9euXcsvf/lLANq2bcs555zDlVdeyamnntqkpHrDhg3MmTMHgL/97W+MGTMmcME3jtrsFhIV5STbP/qRc/Poxx/DT34C113nzDjqvV84KgpycpyE/e23nW2dOjmzkZ50kpOsp6erZ11EKtSXnM/DGfpyAMAYM9Ba+4eqBxljXg9GcCIiTXX06FGWL1/OW2+9xZtvvsmAAQN49913ATj55JO59NJLmT59Oueccw7x8fFNfp/du3dz5plncuDAAWbOnOl2hRa12S7o29dZvHJy4PBhOHTIeR4dDcOGOb3uO3dCdnblZL1tWxg/HiZMqFj69q0YLiMirUudybm1tqDKph7GmGhrbZl3gzFmKjAhCLGJiDTKunXreOmll/joo49YsWIFxcXFvn1Hjx6lpKSENm3aEB0dzbPPPtvs99uxYwennnoq2dnZTJgwgeeff96V4SxearNDw6BBkJvrlF58/nl4911Yu9bZFxfnPN+3D5Yvd3rc161z1h/7VaPv1q1ysj52LHTv7s71iEjLamwpgVeANcaY/wKHgKHA6UC1nhkRkWDau3cvK1asoE+fPowaNQqAFStWcMcddwBgjGHs2LGcccYZnHnmmUycOJHoAA703bBhA2eeeSbZ2dmMHj2aJUuW0L5qKQ/3qc12Sfv2cPHFzpKfD4sWwcKFThnGU06BmBi48EK4/35ISXGO/+or+OILZ8nLgzffdBav5GRnvPro0RXrtDSNXxeJNI2u1mKMSQGuAUYCB4GXrbUh/RVpKN/5LyL1Kyoq4ptvvuHzzz/3LVu3bgXg+uuv54EHHgBg27ZtPPDAA0ydOpWTTz6ZLl26BCWexYsX85Of/IRDhw4xYcIEli5dGrT3au6d/2qzQ8vhw+C59YH9+50SjSUlkJQE//M/cNFFznj07dsrEvUvvnBuLvUOk/HXrh0cfzyMHAnDhzvDZ4YOhV69NCxGxA2BqNYSiFKKUcAYa+2qZp0oiCK5oReJNAcOHGD16tWcdNJJvp7uE044gU8//bTSce3bt2fChAlccMEFXHXVVS0S27Fjx/j973/PvffeC8D555/PU089RYcOHYL2nkEopag2O0Ts3+9Ue1m4EDZsqNjety9ccgn86ldO/XSA8nLIynKS9NWrK9a11VSPj3eS9GHDKhL2YcOgXz/nJlURCY4WT86NMU/VsLkbUGatndmcQDznjwW+BK631n5kjGkH3IPzdWwX4BZr7UHPsbXuq6q1NPQi4aS0tJTNmzfz7bffsn79etasWcPXX39NVlYW4AwbGTJkCODUD//ggw+YNGmSbxk+fHiLTvKzcuVKLrvsMr799luioqL461//yi233IIJcvdkcxr6YLfZnvcIeLvd2tpsa+Gbb5wkfeFCJ+E2xqkC07u3c0xJCdRU3XP/fue133zjJPjr1ztLbTOWxsU5lWTS0pxl0KCKx6mpStxFmivodc5rkIZTG9c/o+8FfFzz4Y32G6Cf3/P5wFJr7YvGmOnAY8BFDdgnIiGipKTEl4S3b9+es846C3CS75EjR1Y7Pi4ujpEjR3LwYEXO9uCDDwY9Ca5NXl4et956K0888QTWWgYPHsyzzz7LxIkTXYmnkYLdZoPa7WYzxhlDPmoU3Hmnc6Pol19WJObWOvsGDYIrr4Qzz6wYZ56Y6NRQnzq18jnz8ion6xs2wLffwu7dFduqatu2cuI+cCD07+/0tvfp4+wXkeBrbM/5cdbazI3hpCEAACAASURBVCrbhgCXWWv/t1mBGDMZGIJzo9JlwCYgC+hirS00xsQABcBwoLi2fdba7Krnbm29MCJuWrJkCf/5z3/YvHkzmzZtYuvWrZSWlgIwZcoUPvroI8AZRz58+HCGDBnC8OHDGTFiBGPGjGHIkCGuT3sPcPjwYR566CHuvvtuDhw4QExMDDfccAN/+tOfWvTGz2b2nAetzfacKyjtttrsyjZscMaVl3lq7gwYANdeC5dfDgkJjTvXoUPOzKabN1esvY/37Kn7tSkpTqLuTdj9H/fpU3PPvkhr0+I951UbeY8i4AqgyQ29MaYDMNtae5MxxltFYCqQb60t9Lx3qTEmC5gClNSxL7upcYhI7crKyti5cydZWVls2bKFzZs3+xLwF154gREjRgDw1ltv8cgjj1R6bf/+/Rk+fDiTJ0/2bYuLi2OLd/7zEHL48GEeffRR7r77bvI9YwOmTZvG3//+d4YOHepydI0TrDYb1G63pKFDndrpzz0H8+fD1q3w61/D7bfDzTc7S0Nve4iPhzFjnKWqgwcrJ+xZWc6Sne0Mtdm1y1mq3P4BOMNhUlOdRL1375qXpCTdpCrSEI1Kzj0NqX9XezTQA3iumXHcAtxVZVsqsK/KtkNAClBax75q8vLySE+v+CNm7ty5zJ07tznxikQcay179+4lKyuLrKwsOnbsyDnnnANAdnY2aWlpvh7wqjZu3OhLzmfNmkWfPn1IS0tj8ODBDBw4MBRLDFazbds2HnzwQZ544gkKCpxy4ZMmTeLPf/4zp512WosOq8nIyCAjI8P7tGtTzxPENhuC2G6rza6uRw8nCb/xRqe84j/+AR98AI884iTqgdCpk1NPfezY6vtKS50/ELzJetX1zp2wY4ez1KZdO6eKTO/eFeuqS+fOSuAl/ASqzfZq7LCWDGABFY19ObDHWrupyQEYczoQZ61d7HmejfP1aDrwP9baiX7HfuV5f2rbZ629t+p76CtSEafnOzc3l65duxLnmSs8IyODxYsXk5WVRXZ2NoWFhb7jTzzxRD72zIpSXFxM+/bt6datG/3792fAgAEMHjyYwYMHk5aWxtChQ8MiAa+qpKSEpUuX8s9//pPXX3+d8vJywKkOc+uttzJ9+nTXxrp7NXNYS8DbbM95g9puq81umOXLncmMZs1yOxIoLnYS8+zsiiTdf9m5EwqqTpFVg44dneEzdS3JyU5deJFQ5MYNof+y1i6vEsRxxphB1tqmfj99MzDG7xdgF+B14P88j/3FA7twftHUtk+kVdu/fz8vvPACO3bsYPv27Wzfvp0dO3aQk5NDaWkpK1asYMIEZ4LIdevW8abfLCcJCQn079+f/v37M27cON/22NhYDh8+TNsIuSNs3bp1PP300/zrX/9ij2egbUxMDBdffDHXX399pV7bMBeMNhvUboeEk0+ueFxWBi+95Iz9PuGElo8lNta5gXTgwNqPOXSo5qTd//nhw7Bpk7PUJSHBSdLrS+I9/RAiYaWxyfl0YHmVbXuAD4EaRrA1yMWA/3+fz4CbgBXAb4wx7ay1R40xbYA+wDKgDHi8ln0iESkzM5MtW7aQk5PjW7yJ9+TJk3nyyScBKCgo4Nprr63xHN27d+eQ30wmc+bMYerUqb6EPKGOu8vCOTG31rJ27VpeeeUVFi1axHq/UhVDhgzhZz/7GZdeeinJyckuRhkUwWizQe12yHngAWfYy89/7k5y3hDx8RV112tiLRw44FSU8Y5vr2nZvds57sCByvXha5KY6Ez01LOnMzSoR4/Kj73Pu3XTDa0SOhqUnBtjrgP+AsQbY26p4ZD3mhqAtTavynuVAXnW2u3GmCXAqcAS4DRgkbU2x3NcrftEwkVxcTHZ2dns2rXLl3D7P160aJEvYfztb3/La6+9VuN5EhMTfY9TU1O58sor6dOnD7179/ate/XqVS3BHjduXKUe8khSUlLCZ599xltvvcWiRYsq3XyakJDA+eefz89+9jMmTpzo+tCVQAtmmw1qt0NNSQl4iiDRubOroTSLMdCli7PUlsCDk8Tv3193Au9N4vfvd5aaSkdWlZRUPXlXIi9uaFBybq19yBizHJgD/KPK7mJr7e6AR+a4GrjHGDMK6AnMbeA+EVeVlJSwZ88edu3aVS3x/sEPfsAvfvELAL744gtOOumkWs+zY8cOX3I+ceJECgsLSU1NJTU1lZSUlEoJuFdsbKz/jSmtyvbt21m6dClLly7l/fffr/QtQbdu3Zg5cyazZ8/mlFNOoU0E/3Z1sc0Gtdst6osv4LrrnLro8fHO40hnjJNIJyU5JSZrU17u1HvPzXXKRHqXmp7n5zvj9/fta3gi3727k6h71/6L/7akJAiB6rASRhp7Q2gba21JDdv7WGtrmUTYfbq5SAKlvLycvLw8X6Ltv8yfP99Xn/sHP/gBn3/+eY3nOP/883nxxRcBp0LID3/4Q1JSUnwJt3/yPWbMGOLj41vs+sLNzp07WbZsmW/ZVGWg6tChQ5k+fTrnnHMOJ510UkjUT2+sZt4QqjY7Qm3b5tQ5/89/nOe9esGrr8L48e7GFa7KypwEvbbk3f95Xp6T+DeU9xuBqkl8bUl9165K5sNZ0G8INcbcDLxlrfX+Hfn/jDHRVQ8DTsGpbysSlqy17N+/v1rCPWnSJE499VQAFi9ezHnnnUeZdyaQKm6//XZSU1MBZ2hJz549SUlJITk52Zdwp6am+koOAvTt25fvvvsu+BcYAcrLy8nMzGTFihUsX76cZcuWsXXr1krHdOzYkdNOO43p06dz+umn069fP3eCdYna7Mh17JgzvtpbnzwxET77zOktv+YauPVW57E0TXR0xdCVGiYursSbyO/d6yTq/kvVbXv3Vgyt2b8fNm5sWDxdujSsV75bNyeZj41t/mcgoaO+v836A/53iI0FOgC5VJTmMjhfT4qEHGstBw8eZPfu3b6E+8CBA1zn993vCSecwKpVqygqKqr2+ptuusmXnCcmJlJWVkbXrl1JSUmptnTwmwXk5ZdfjrhxzC1t7969rFixwrd8+eWXvvrjXvHx8Zx44olMmTKFKVOmMG7cuIgertIAarMjhLXOREAffghvvQXvv+/0pubnO+Od4+PhlVdg8uTGzxIqzeOfyDdEWZkzXKa2BL7qtn374PvvnaW+qjVenTvXnszXlNSrik1oqzM5t9ZWLfnwV2CdtbbYf6Mxpo5bN0SCo7CwsFIv9+jRoxkyZAgAL730Erfeeiu7du2qVLsbICoqiquvvproaKdDsaSkhKKiIjp37uxLtL3DSqZMmeJ73aRJkzh27JivRnhdlJg3nLWWnJwcVq9e7VtWrVpFdnZ2tWN79erFxIkTmTx5MlOmTGHUqFFhOVQlWNRmh7/16+HOO52kfFeVIpNjxzrb+vZ1np95ZsvHJ40XHe0kxt27N+z4sjInMW9Iz7x3KShwloZOuhwfX3+PvP+iuvItq7G/1Y631n7lv8EYcxwwEGjALRQi9SsqKiI3N5ecnByOHj3q67m21jJ9+nR27NjBrl27qvWi3nfffb7k3Frrq87Rvn17X7LtXYqKinyT5ixevJj4+PhKPd81iYmJUSLYTMXFxWRmZrJmzZpKyfi+fVUnjnSGqKSnpzNx4kTfkpJS4yTAUju12SHKWsjMhGXLnETp4osr9j3/vLPu1g1OOQVOO81JxD2j5iTCRUc7Q1W6NnCeyfJyJ5mvKWmvLZk/dMhZqowMrFWHDrUn7jUl9h06aKbX5mhoKcWTcL4KneaZDtpfIk41gDcCHJtEGP8KJgMGDKCrp+V56aWX+Oc//+nrAc/Pz/e9JjU1lZ07dwJOb/S6devY5elOio2NrZRwDxgwwPe6008/nfXr15OSkkKnTp3q7Mnu2VPf8Afa0aNHyczMZP369ZWWLVu21DhmPykpidGjR1dahg4d6vt2QxpHbXboKS+HtWudWT29y969zr7x4yuS86FD4dFH4cQTnXKCSnCkPlFRFdVrPP1TdbLW6WWvLXGvKbE/csRZavhCs0Zt29bdI+8tS9mzp7NfY+Yra2g34AHgWZzxjFWnNygGnghkUBJevBVMCgsL6d+/PwBHjhzhpptuqjTsZM+ePXirA73wwgtccMEFgFNxY+nSpb7zRUdHk5yc7CsVaK31JdevvPIKHTt2JCUlhcTExFqT7oSEhDon1JHms9ayd+9eNm/ezObNmysl41u3bqWmSlDGGAYNGsSoUaMqJeKpqakaChRYarNdVlrqDE/wjoK76Sb4+98rH9OzJ0yZ4vSOexkDnkqrIkFhjHOfQkICDB5c//HWOr3sDemR9+47dqxi1teGSEpyZnT1Juy1LYmJreMP1obWOV9rjJkK/I+1Vo16K/bBBx/w0UcfkZWVRVZWFtu3byc3N5fS0lImTJjAihUrAIiLi+OJJ56g3K/elDGGHj16kJKSUmnc9rnnnstxxx3n6wHv2rVrrT2mP/jBD4J7gVKJtZb8/HxfAr5lyxbf482bN1eqI+4vOjqatLQ0hg0bVmkZPHgw7dq1a+GraH3UZre8oiJYudLpEV+2DD75BB5+GH76U2f/pEnw2mtOMn7yyc560KDWkWhIeDMGOnVyloED6z/eWqeXva7e+L17nUmicnOdx94a8+vW1X3uNm0qJ+upqdC7t7P06lWxDuMJrYFGjDm31hZQS2+LMeZEa+1/AxaVtChrLbt372bLli1kZWWRnZ3tS76zs7N56aWXmDhxIgBvvvkm999/f7VzJCUl0dlvarqYmBgef/xxEhMTfUl3jx49aqykMXDgQAY25H+8BEVRURHbtm3z/dy96++++47NmzdXG9vvLyEhgbS0NNLS0hg8eLAvCU9LSyNW31O6Sm12y/jb3+C995yyhseOVd737bcVj88/Hy68sGVjE3GDMdCxo7N4vkyvk7c0ZW5uxeJN3KsuBQUN65Hv1q0iWe/dG/r1c/4YHjQIBgwI/Rtc66tzvhjoWM85DDAcaOB9yOKWAwcOsGnTJjZv3kxcXByzZ88GnIlw+tfxP+i7777zJednnHEG8fHx9OvXj/79+9O3b1+Sk5OrTQsPcPnllwfnQqRRSkpK2LlzZ6XE23+9a9euGoegeHXq1MmXgA8aNMj3OC0tjaSkJA1HCSFqs4Nr715491245JKKHu8333R6ycEZI+7tFT/ppMo3cEZFtXy8IuHAvzTlqFF1H3v0aOUEPienIlnfsQN27nS2eXvqv/665vOkpkJampOsjxrlVEIaNcq5kTUU1NdznokzdjGHihq5VRng7EAGJYGxePFi3njjDTZs2MCmTZvIy8vz7Rs/frwvOe/Tpw+pqan07t2b/v37079/f1/y3a9fv0pTw0+bNo1p06a1+LVIzay15OXlsWPHjlqXnJycWidOAmcIivdn7/9zHzBgAGlpaXTr1k0JePhQmx1gubmwYIEz++annzpf2Q8b5vwyB/jtb50hLSed5PTWiUjwtGvn9MbX1SNfVubM5OpN1nfsgKwsp8zkli1OhZqcHGf56KOK10VFOTfUnnYazJjh3Avi1h/V9SXndwOHrLXVZ2fxY4zRXf8uKC0tZd26daxevZq1a9eydu1a7rnnHkaPHg3Axx9/zBNPVHyr3a5dOwYPHszgwYMZ451mDqfut7ciioQOay0HDhyoM/HeuXNnjZMn+TPGkJqaWi359v4h1qtXL5WIjBxqswOgqMhJxp991ukp9946ExsLp57q/PL3Oussd2IUkZpFR0NKirN4vvSvpLQUtm93EvWNG2H1avjqK2cY2vr1zvKPfzj3ibzxRsNLWgZSfZMQ5de1388sdPd/izh8+DC33XYbX375JV9//TVHjx6ttH/VqlW+5HzWrFn06tWLESNG+G64jNJ3qyHh6NGjlSrZeJecnJxKz48cOVLvuRISEujdu3edS0MmTpLwpzY7MEpK4Jpr4MAB5wa0GTOcUofTpzs1yUUkfMXEOOPOBwyAH/2oYvuxY7BoEdxwgzMG/vPPYfFicGOEbrO6y4wxfYBrgetQQx9QZWVlrF69mo8++oh9+/Zx5513Ak7v95NPPsnhw4cBGDRoEGPHjmXkyJEcf/zxTJo0yXeOyZMnM3nyZFfib62Ki4vJzc2tMfH2T74PHDjQoPN17NixzqS7V69edOxY3xBjEYfa7Jrt2QMPPAC33+5UeejY0Zmls7wcLrjAnZ4zEQmO4mLYtq1imIu3B33VKqdijFdysvNNmRsanZwbY6KBGcCVwDRgD84YR2mmzZs3s2TJEj744AM+/vhjDh48CDhlCf/whz/Qtm1boqOjeeihh0hOTiY9PZ3ExESXo4581loOHjxIbm5urcvu3bvZvXs3e72zitSjTZs2lSZQSklJqTaLaXJyMp07d9Z4b2kWtdm1KyqCu++Ge+5xSr917Qq//rWz7+qr3Y1NRBrOWjh4sHrZxl27KsaXexe/2++qSUqCE06As8+Giy5y/lB3Q4OTc2NMf+AK4GdAErANOAd4B1Dx6WZ67rnn+Km3IK7HgAEDmDp1KlOnTq1UTWPOnDktHV5EOnbsGHv27Kkz6fYux6rWSKtFVFQUPXv2rJZ4V03AExMTNcRIgkptdt1WroTLLqsodzhjBpxxhqshibR65eXOhEfff+8MKztwoOKx/zo/v3Iinp/vDEdriKgop8yit1rLoEHO49GjoU+f0Jh7oN7k3BjzPzg9Lj8EtgD/BzwD/M5a653W8dOgRRhhrLWsWrWK559/nl69evFrTzfNlClT6NKlC9OnT+eMM85g6tSp9O7d2+Vow09ZWRn5+fkNSrgbOrQEnOElPXv2rHXxTq7UvXt33VwprlKbXb/HH3fGlJeWOr+UH3/cKX8oIk1jrTNm++BBJ7luyPrgweqJd0FBxQ3YjdWxo1MxqWtXZ92tm1MyMSXFWXuXHj2cm0ZDWUOyiKM40z1/Clxhrd0EYIypvTCyVFNQUMDTTz/No48+ysaNGwHo27cvN954I1FRUfTp04e8vLxaZ8Zszay1FBQU1Jpk+/d+7927t9KspHWJiYmpM+H2T7w1rlvCiNrsOrz8Msyd6zy+7jpnWEuoT0giEmjFxc5QrsOHa15X3Xb4cP3Jdh0VexslPh4SEqBLl5rXCQmVE3BvQh7us4L6qzc5t9YuAZZ4biS60hjTD3gd8E31aIwxtq5ZTFqx7Oxs7r33Xp555hnfTZzdu3fnwgsv5JJLLqk0nri1JeaFhYU1Jtg1LcXFxQ0+b2JiIsnJyXUm3MnJyXTp0kVDSyTiqM2u249/DNdeC+npzrAWkVBlrTPpTk3JcmMS65rWDR0C0hhxcdCpk5Nc17f2Pu7SpXLy3bmzU02ltWvwR2Ct3Q783hgTA8wEehlj7gCeBH4EZAQnxPC2Zs0aHn74YQCmTp3KL3/5S2bMmBGxQx9KSkrYu3dvvT3cubm5HDp0qMHnjY+Pr9abXVPS3b17d9q0aVP/CUUinNrsmkVFwUMPuR2FRIrycigsrEiK61v8E+j6lsJCJ0EPhuhoZxhIx47OrJgNWdeXbMfGBifW1qjRGaK1thR4BXjFGDMIuBFnfGOrbOirOnbsGJ988gmneurvzJgxgz/84Q/Mnj2b448/3uXomqa8vJx9+/Y1qJd7n38donrExcXVmmT7J+A9evSgQ6jMqSsSZtRmw4YNThWWp56Cnj3djkZaWllZwxPixibWVaYaCbi4OCc59ibIjUmm61rHxobGjY9Ss2Z131prtwDXG2NWBSiesLZ8+XKuuOIKsrKyWL16NcOHD8cYw5/+9Ce3Q6vRkSNHyMnJqbeHe8+ePXVO/+4vKiqK7t2715po+y8qEyjSslpjm20t/OpX8P77cNddTj1zCT2lpc3rZa5raWCxrSZr164iga5p8SbEjV3at9cQj9YqID92a+2zzT2HMeY84D6gE/Av4NfW2lJjTFfgrzh1edsAt1hrSzyvqXVfSyorK+POO+/kj3/8I+Xl5QwdOrTazJ0tyVrL/v37ycnJYefOnb61/+OcnJxGVStJTEyscziJd1/Xrl1b3dh5kXDTmtrs//zHScy7dIHbbgvmO0U+702EgRy64V0acVtRoxnjJLoNTYobk0y3b+8MlRIJpJD4m8xz49JMYDYwHHgE2AHcC7wM/K+1doUx5irgz8BvPS+ta1+LOHLkCOeffz5vvfUWxhh+97vfcfvttxMbxMFXBQUFZGdnk5WVxY4dO6ol3Tt37mxQXe64uDhSU1NrvHnSPwnv3r27pn8XEZ9warPvv99Z33hj5M/0aa0zsVKgh254l9LS4MUeFdW03uWGJNPt2mkIh4SXkEjOgb7A5Z6xkauMMSOAU4wx/wWOs9au8Bz3BpBpjPkLcHxt+6y1R1oi6Pz8fM466yy++OILkpKSWLhwIdOmTWv2eY8cOeJLvmtaf//99/Weo1OnTvTq1YvU1FR69epV4+OkpCQNKxGRpgiLNnvPHnjrLWd87VVXBeMdAqe83EmMCwqcxVvzubbHVbd5k+pAlbOrSUxM0xPo+hLpuDgl0CJeIZGcW2s/rrIpB4jHmUQjy++4HE8ymQ6cUMe+ZUEOGYCtW7eSmZlJ3759effddxk8eHCDX1taWkpWVhaZmZls3LiRzMxM31LfFPDt2rWjX79+9O/fnz59+lRLvlNTU4mPj2/u5YmI1Chc2uzFi53e5GnTnFrILcFap+5zfj7s21f3ev/+iiT74MGmT77iLza2+Ql0bcm0qnGItIyQSM5rMB6YB8wFqpb/OASkAKl17GsREyZMYPny5SQlJZGamlrjMdZadu/ezddff+1bNmzYwJYtWyippdBobGwsffv2pX///r4k3H/dvXt39XiLSCgJyTb7vfec9YwZzT/XsWOwezfs2lWx9n+cm+sk3Pn5TR/+0aGDU++5c2dnqe+x/7b4eOf1qiQrEv5CLjn3lPrKs9au8cxoV/XOyjigBKhrXyV5eXmkp6f7ns+dO5e53inimmD37t0kJycDMHLkyEr7jh49yooVK1i2bBmff/45X331Va094b179+a4445jyJAhldapqamaHEekFcvIyCAjw1fpMKRHSodym33jjTB2rNNzXh9rYedO2LIFtm6F776rvN6/v+Hv26GDM749KanudWJi5QRblTlEwlOg22wTSpPEGWOigbuA26y1RcaY3wKnWGt/5NlvgGPAKcCU2vZZaz/1P296erpduXJlQGL86KOPmD59OnfddRc33HADABs3buTVV19l6dKlrFixotpslp07d2bMmDG+ZcSIEQwePFi1u0WkXsaYVdba9PqPbHnh0GbXJisLli+H1asrlroKWMXEQHIypKRUrP0fJyc7Q2eSkpzx0yLSOgWizQ61v9NvBu611hZ5nr8H+N/G0w8oBlZ61rXtC4qioiKuuuoqioqKyMnJ4b777uOpp55i/fr1vmOMMYwaNYqpU6dy4oknMm7cOPr166dhKCISiUK6zfZXXg4ffwwLFjjDXbKyqh+TlATHHQcDB8KAARXrAQOgRw+VzBORlhEyybkx5nc4jXQ7Y8wAYCqwHNhtjBlhrV0HTAcetNYWAyuNMbXtC4r77ruPzMxMEhMTefjhh321zBMSEpgxYwYzZ85k6tSpdOnSJVghiIiEhFBvs9euhbffhnHjYMcOuOMOZ3iKV0ICTJkCEybA6NHOkpysiiEi4r6QSM6NMbcBd1TZvNFa+5Qx5gLgD8aYbKAjcKvfMXXtC6h9+/bx17/+FYD9nsGH06ZN47rrruOMM86gje7CEZFWIhza7E8/hVtugZ49nZs1AXr3hksvhZkznbHomi9NREJRSCTn1tq/AH+pZd824OeN3Rdov/rVrygsLARg/PjxPPTQQ0yYMKEl3lpEJKSEQ5u9Y4ezzs11esTvuQcuukgJuYiEvpBIzkPdpk2bWLBgAQAzZszglVdeUU+5iEgIe/99Z921K3z1ldODLiISDpScN8Cf/vQnAJKSknj11VeJVteLiEhI27DBWV95pRJzEQkvuve8HgcPHuTf//43AO+//74ScxGREHfsmDPjJsCpp7obi4hIYyk5r8frr79OYWEhJ598MqNHj3Y7HBERqYf3BlBwbgIVEQknSs7r8cQTTwBw7rnnuhyJiIg0hHcyIWOc+uQiIuFEyXkdrLV88sknAHTv3t3laEREpCGOHHHWEydCp07uxiIi0lhKzuvwySefUFZWhjGGCy+80O1wRESkATxVb+nQQZMKiUj4UXJehxdffBGArl27EhOjwjYiIuHAm5y3b+9uHCIiTaHkvA4ff/wxACNHjnQ5EhERaajly531xo3uxiEi0hRKzuuQlZUFwCmnnOJyJCIi0lD79jlra92NQ0SkKZSc16K8vJxDhw4BcMYZZ7gcjYiINJT3htC2bd2NQ0SkKZSc12Lz5s1YazHGqL65iEgYOXrUWcfFuRuHiEhTKDmvRU5ODgDjx48nKkofk4hIuDh2zFkrOReRcKSssxabN28GYNiwYS5HIiIijVFU5KyVnItIOFJyXgtvz3lvzf0sIhJWioudtZJzEQlHSs5r4a1xvs9727+IiISF5GRnnZbmbhwiIk2h5LwW+fn5APTo0cPlSEREpDG8zbZGJYpIOFJyXosjnlpcQ4cOdTkSERFpDO+wlthYd+MQEWkKJee1KPa07mn6XlREJKzk5jrrw4fdjUNEpCmUnNfAWov1TC3Xt29fl6MREZHG2LDBWWdmuhuHiEhTKDmvwf79+32PExISXIxERERERFoTJec12LFjBwDGGIwxLkcjIiIiIq1FjNsBNJcxph1wD3AI6ALcYq092Jxzeoe0dO7cudnxiYhIhWC02SIikSQSes7nA/+11v4OeB14rLknjI6OBqBXr17NPZWIiFQW8DZbRCSShHVyboxJAS4G3vBseh+YYYzp15zzHjzodOJ06tSpOacRERE/wWqzRUQiSbgPIkFioQAAIABJREFUa5kK5FtrCwGstaXGmCxgCpDd1JNu2bIlIMGJiEglUwlCm11eDpdc4jxOSADPHHIiImEp3JPzVGBflW2HgBT/DXl5eaSnp/uez507l7lz59Z60i+//BKAnTt3BipOEZEGy8jIICMjw/u0q5uxBFhQ2uyXX4YXXqi+ffr0ZkQqItJAgW6zwz05t8DRKtvigBL/Dd26dWPlypUNPmlBQQEA7dq1a2Z4IiKN55+MGmMiqR84KG22t/ptVBQ8+CC0bQu9esFppzUvWBGRhgh0mx3uyXkOzt3+/uKBXc056aFDhwAl5yIiARaUNturQwe45ppAnElExD1hfUMo8CGQ4inNhTGmDdAHWNackxYXFwMQGxvb3PhERKRCUNpsEZFIEtbJubU2F1gCnOrZdBqwyFqb415UIiJSE7XZIiL1C/dhLQBXA/cYY0YBPYHa7xoSERG3qc0WEalD2Cfn1trvgSvdjkNEROqnNltEpG5hPawlWCZPngzAD37wA5cjERGR+gwYUHktIhLOlJzXwHsjaFxcnMuRiIhIfTp2dNYdOrgbh4hIICg5FxEREREJEUrOa/DNN99UWouISOjatavyWkQknCk5r8F+z3Rz3rWIiISufM98fPv2uRuHiEggKDkXEREREQkRSs5FREREREKEknMRERERkRCh5FxEREREJEQoOa9B165dAUhKSnI5EhERqU+nTs46Pt7dOEREAiHG7QBC0axZszh48CCzZs1yOxQREanH2LEwYwYMGeJ2JCIizafkvAbnnXce5513ntthiIhIAwwdCq+/7nYUIiKBoWEtIiIiIiIhQsm5iIiIiEiIUHIuIiIiIhIilJyLiIiIiIQIJec1yMjIcDuEgIu0a4q06wFdUziItOuJNPr5VKfPpGb6XGqmz6U6Nz4TJec1iMR/nJF2TZF2PaBrCgeRdj2RRj+f6vSZ1EyfS830uVSn5FxEREREpBUz1lq3Ywg6Y0wesK0RL+kK5AcpHLdE2jVF2vWArikcuHE9fa213Vr4PV3VhDbbK9L+vQWCPpOa6XOpmT6X6hr7mTS7zW4VybmIiIiISDjQsBYRERERkRCh5FxEREREJEQoOfdjjGlnjHnQGHOnMeYRY0wnt2MKBGPMX4wx1rOscTuepjDGTDPGrDDG9PPbFtY/r5quybM97H5expjzjDFZxph9xpi/G2NiPNu7GmMeM8bcbYz5P2NMG7djbajarsmz719+P6PX3YxTwr8tCJRI/H8YKMaYWGPMGmPMVM9z/ZsBjDGTjDE3G2Nmev6dtOrPxRgz1BjzsDHmRmPMfGPMaM/2Fv1cYuo/pFWZDyy11r5ojJkOPAZc5HJMzWKM6YBzM8M0z6am3GTlKmNMD6AjMKHKrrD9edV2TeH48zLG9AFmArOB4cAjwA7gXuBl4H+ttSuMMVcBfwZ+61asDVXXNXn27aXiZ5TpSpDiL2zbgkCJxP+HAfYboJ/fc/2bMeZnwEBr7W1+2/5J6/5cngPOtdbmeP5PvQMMpYX/veiGUA9jTAqQBXSx1hZ6ehwKgOHW2mxXg2sGY8wNnoePWmuPuRpMMxhjooAyoL+1NjsSfl5Vr8mzLex+XsaYk4DPrLWlnuf34CQHdwCvWmtTPNtTcRLZHtbaI27F2xC1XZO19ixjzP3AJ8C/vfvFPZHQFgRCJP4/DBRjzGRgCPAH4DJgE63834wxZgrwe2Ca9SSC+r8ExpgjwDhr7UZjTDdgDZBOC38uGtZSYSqQb60tBPA0cFnAFDeDag5jjIH/z959h0dVpo0f/97pvZEQSOjFUKUF1LWxlteCXXRtu6Cr6PLKFn+67qq76uq6uviudcEF6+qKDRtgd7GgWJAmvYaSQEiB9J7n98dzJhlCAiSZTEnuz3Wda2bOOXPOMzPJPfc85z7P4Qrg/4BcEQnYX7/GmPomsyYS4J9X09cUqJ+XMebLJklqNrATOA37mbjWy3buZnqxeW3S0msSkSjgfOA15/HpPmmgcjeRAI8FntAZ/w89wTkaOdkY86zb7Ino38w/gPXAEyLygYiciL4vAG8AT4tILHANMAMfvC+anDdKBwqazCsB0nzQFo8w1vHYMokngP+IyDk+bpan6Oflv8YDT9G5PqPx2KMZ5caYQdjXtgh4X0RG+rZpXV5n+jvzpM74f9gWtwMPNpnXpd8TEckAxgBzjTE3A58CHwC96cLvi2M6UAl8D5QaY+bjg78XrTlvZICKJvPCgRoftMWjjDH7gbucntnfAu/7uEmeoJ+XHxKRQUCeMWaViHSKz8j9NbnmGWP2ADc4Pek3Azf6qn2q88aCtuqM/4dtISJnAcuMMfuaLOqy74ljOFBojFntPH4SuBt9XwAigXnYhPxREdmOD94X7TlvlA0kNpkXC+T4oC0d5XGgj68b4SH6efkZEQnGJql/cGYd9Bk5PzaiCaDPqJnX1FRAfUadVFeIBUetM/4ftsOtwLMiki8i+die4Xewr78r/82E4NY5a4ypADYDwXTt9wXgJWC+MeYv2CMu84A9ePl90eS80WIgTUQiAZyhpvoAn/u0VZ5VDyz3dSM8RD8v/3Mr8LAxpsp5/DGQ4ba8H1ANLPNyu9qj6WtqKtA+o86oK8SC1uiM/4dtdRUw2m3KAa4HnqNr/82sBuKdUcNcarHvT5d9X0QkGRhljDngzLofKMaOmubV90WTc4cxZi+wEHCd4HUG9tdTdsvP8m8ikiYiV4lIkNNb8v+AO33drrZw2g8g0Dk+r6avKZA/LxG5A1gBRIrIABG5DjgA7BGREc5qZwNPGGOqfdXO1mjuNYnICBG5wFkeBkwB/u7LdnZ1nSEWeEpn/D9sD2NMnjFmt2vCjo6VZ4zZSRf+mzHGbMCWS04GEJEEIAz4D134fQEKgUpnRCOXAuyPGa++LzqUohsRScR+0WYBPYA/GmNKfdqodhCR4dh/wGJgCfCo808ZUNzOmp6FHZ/3SWNMXiB/Xs29JqA7Afh5ichd2OHa3G0wxgwVkb7Y4cuysOO63xkIww+29Jqw9eXzsGfqfwM8ZIzpSod8/VIgxwJP6Yz/h54mIlnAVGPMZ139b8bpJX4MewSlD/bk0HX6vsgo7EmhPwCpwBfGmM+9/b5ocq6UUkoppZSf0LIWpZRSSiml/IQm50oppZRSSvkJTc6VUkoppZTyE5qcK6WUUkop5Sc0OVdKKaWUUspPaHKulFJKKaWUn9DkXKk2EJHTReRrEZl4hPVGiMhrIvInD+57pqe2pZRSviAix4rI656MjS3sZ6GIXNGR++goInKGiCw90veM6nxCfN0ApdpDRIYCf8Ve2Ws8cCzwc2PMSx28613AqKNYbyeQCAR7YqfOZYNrPLEtpZQCEJHxwN+AU4CPsBfZ6gGEArONMR92wG63APF4KDYexixgTUds2InHlwBPAZuApc6icOBSYIIxJqsduzja7xnVyWhyrgLd28DlxphVIhIEPOeNnRpjNolIwVGsVywiezy461OBLzy4PaVUF2eM+V5EXgGGGGPOc80XkUnAWyJyqzFmlof3We7h2AiAiJwP/OhKio0x73l6Hy7GmBrgVRGZAXxmjLnLrR3vtGWbItId+Kkx5lVjzEYRyfdQc1UA0bIWFbBEJAk4BigDMMbUYy9R7a3L3tZ7eL2j8VNgsQe3p5RSALVNZxhjFgH3AI+ISI8O2KdHY7WI9ML2lHvbIe8d8AmQ3ZqNiEg48G8g0m22J78/VIDQ5FwFLGNMIbAR+EBEjnfm7QC+E5EEEblPRL4TkQtEZK+IrBeRsa7ni8iVInKniHwmIg+4zc8UkXtE5FUReVNE4p35ESLyTxH5s4g8gz0k2yoi8nunTZeLyHgRedvZ3n1OG18VkRRnfpGI3N1kE2HGmCpnW/9PRG4Rka9E5MFWv4FKKXVkz2DLWy4DEJGBIvIXEXlaRD51EmJEJNyJZXeKyJci8jPXBkRkuoj8SUS+EZGbW9pRczFZRCY7MXOqiGwWkQdFJE1EnnPi6ZcicpqziauAXsDtInKRc27QJyIyxW0fN4jIX0Xk3yLygohEi0i8iNwtIstE5DQR2Sgi20SkX1veMOc9mej0rCMil4jIYyLyfyLygYjEOK/hMRFZ4Hzf7HHe42HAZSJyi9smE53vohIRuaGFfQaJyEMiMs3ZhxGRt1yfjwowxhiddArYCRgO7MD2LjwLdHfmB2MDXSG29i8B+BhYi/1RegJwl7NuIrb3fRIQB7zitv2vgJnO/UeAac79VKAOG4CP1Mbnsb1P3bG1iVHO/HDgQ+AdIA37pVID3A1EAeOcx671XfX0rtc9y7kfC/zB15+FTjrpFLgTMBXY3cKyXOBxJ66+CQQ7818GXnfu/wO4xLl/HZDn3L8SuMa5f4wTq0c6j58H7nHutxSTk4ADwKPABOAM4GHgz866twEL3NpqgH7O/UhgPTDVeXwW8LLbuguBf2FLfC92vi8ucV7nB8Bfj/K9+8z5rngQmImtPz/DbXk+0Me5vxb7nRQO3AnsBo4DfunE/c9c7XXWzwL+z1l2JbCzhTZcDyx07qc67/P/+PrvSqe2TVpzrgKaMWatiIwE7gemA+eLyBnG1qDnAcXGmPkAInIPsAQYgP0iqheR3zqb+hD7JTAJ6OY2f73z3BjgV85zMcbkisiuVjQ1GXgC+KUxptzZRpWI7AW2G2NynP3sBT43th5zBfZLIxl7Yukk7BcJQDnwCxHZBMzG/jBRSqmOYLAJ6wSgLzBDRAAKgHixD27EdiwAvIAt6wAba5c1ibU9gR+b7GMqzcRkY0yhiBQB7xhjvgMQkZ1AmYj0xPY0xzTbaGMqRGSf26ybne26PA285sw/gP2+eNPZxzKnnUdrsXFqzkXkVWyHkMtlxpidInIiNimPceJ/DrDVGPMt8K3z3Oa2vcD5Tvj+MG0aC1Q4rztXRAo5uDxGBRBNzlXAEpFgoLexJ/782ik1eQ+Yg+2JaGqtc5sA9AFeNMa87Mx71Nnm7dhk+dEm+xqHDaolbWzuMcBJwEBgldv8pjWXdQ0LjKl3ArWr/KybsaU8GGO2i8hU4Engd8DVwD6UUsqDRCQBe9RvLTZuFjQTH1OxPbthAMaYOmyHAs5z7jXGfO08Pui5bpqNyQ7DwbFyN3A7sBn4Buh/mJfg/rzBwH/dHm/Flux059BYXEvbR5JZjv0R43JARB4CXse23ZWBN31dR1JLy3nb+8BsEYnAfhYV2F54FYC05lwFsoY6SABjzCrgDmBIC+uHYQPhFmAPMNl9oYgc58yfJCKRTeYXOw+HtrGtXwOPAf9xgmeriEgabicXiUg6dqSaY7AniL7cwlOVUqo9pmCP1L2OjY8nOj3WQEN8zAeqgbPd5o8UkViaxFoRCXE6O5pqKSY35ykgy9ghc6ta8Vp2YnvaG3YBFGHLdjzGGFPvdKAMd37cvAf8xRizzJP7abLPBcAC7BHeK4FTjTFFHbU/1bE0OVeBboaTuLqkY2sFXZJEJMy5PxFbb3gAmAdcLCJPishEp1ejDBtEo4GFInKWc/LNEGxCvwK4T0TCRCQZW+vdXUSCReQPInJNC20UZ/oz9gvsH80sa7p+08fnYesjXQYDFxljioEZzTxHKaVa45AeWRE5BbgXuNEYk4ftpc4F3hOR852TPk92espfwY7qcpGInA5MMcaUYGPtr0XkXhE5FdtJ4SoJdI9/LcVkF/de7DFAiohEAScCkSIywFlWgz2B0tVJ476Pp4BLRSTReXwcMNcYU4vNh5qNxSJyvNOu8BbeuyCa5FPOD5hJ2KOlidhyyQygNxDldrJp09756mbaf3Cjmql9cT6LOOBT7NGB1vxoUX5Gy1pUoOsNrBWRRdhDfgLc5La8FrhTRA5ge5lnABhjPhaRX2MPjV6EPaFyDYCIXAj8E/tl8yxwmzHGiMjl2GGu1gKvYg9PjsZetONybLnMQRc/EpEJ2BOd+mAPc74N3Ou050PgeKCfiAzE9uj0AC4XkQ009iJNBaKNMXOavPbnxF48pBZ7ApZSSrWaE0euwiaQc4E87EmF8cA5xpilAMaYaic+zsXGunewteYAv8F2bPwbe6TQNULKM9jY9yts+d2vjDH7nOTzeGCgiLzUUkwWkcnYOuvrRWSDMSYbe3Lqw9jE/N/AhcBIYBvwIraX/xdiR/EaDpwrIu8ZY94UkcHAfBH5Cnv09U6nh/8yoIfz+tZiL8jUU0SGO+28AltH/7bb+xaKjf2jgHQneTfYWu9LnG2uBD4HvseW6ryPjekLsSehjhKRi4wxru2+ir0gVKmTzPcAJovIdmyPOM7zm17TYyf2hNeLgAggSEQ+wZ4U6q3hhZWHiH5mqrMSe8nj540x/by0v9uMMTO9sS+llFLe45TYhBtj/PIicCJyG3bknCzncQT2qMdDrnOVVODQshbVmTVXMtIxO7KHf5d7Y19KKaW8xznhdawfJ+YRwF1AP7fZyUCZJuaBSctaVKfk1Ptdgz1MeZkx5vUO3uVW53CrUkqpziXEGDPb141oiTGmUkRmAI87gwXkYEtw7vJty1RbaVmLUkoppZRSfkLLWpRSSimllPITmpwrpZRSSinlJzQ5V0oppZRSyk90iRNCk5OTTb9+/Vr9vLy8PFJSUjzfoACl78eh9D05lL4nh2rPe/LDDz/kG2O61BuqMfvw9HV2Lvo6O5cffvih1BgT255tdInkvF+/fixb1vqr5mZmZrbpeZ2Vvh+H0vfkUPqeHKo974mI7PBwc/yexuzD09fZuejr7FxEpN1VKVrWopRSSimllJ/Q5FwppZRSSik/ocn5YUybNs3XTfAr+n4cSt+TQ+l7cih9T7yjq7zP+jo7F32dnU5eezfg9YsQiciZwP3Az4wxWc68SODvQAmQCNxujCluzzJ3mZmZpivUOSmlOh8R+cEYk+nD/WvMVkqpo+SJmO3VnnMRSQVigAlNFs0Clhhj7gDeAf7lgWVKKaXaQWO2Ukp5ny96zoOAOqC/MSZLRNKA7UCiMaZcREKAImA4UN2WZa7eHRfthVFKBSo/6DnXmK2UUkcp4HrOAYwx9U1mTQTyjTHlzvJabAA/tR3L2m/XLjjxRDjzTI9sTimlAlEgxOySFSUsG7eMDddtaO+mlFLK5/xhnPN0oKDJvBIgDaht47KD5OXlkZnZ+CNm2rRpRz4xITgYvv4agvScWaWUd82ZM4c5c+a4Hib7si3N8LuYXb6hnNLlpVTuqDzKl6CUUp7j6ZjtD8m5ASqazAsHapz7bV3WICUlpfUD3/foYW/r62H3bujVq3XPV0qpNnJPRkUk38fNacrvYrYECwD1ZU07+ZVSquN5Omb7Q7dwNvasfXexQE47lrVfUBCEhtr7K1d6ZJNKKe+bN28et912G9OnTycmJoZbb73V100KdH4Xs4PjgwEwdd49h0op5Xkas/0jOV8MpDlDbCEioUAf4PN2LPOM6Gh7u3atxzaplPKe3NxcXn75ZWbOnMmsWbNYtGgRQU6p2tNPP+3j1gUsv4vZIXH2ILAm50oFNo3ZlteTcxER110AY8xeYCFwujP/DGC+MSa7rcs81thEp5Nn82aPbVIp5T179+5l48aNFBUVAXDqqacyYcIEPv74Y+bNm+fj1gWGQIjZIQlOhaZWtSgV0DRmW94e5zwWuMl5OFVEUpz7vwIuFJE7gXMB9zN/2rqs/VJT7e22bR7drFLKO4YPH05CQgJjx47lww8/BOCCCy5gwYIFZGVl8fDDDwOwYMECZs6cyVlnncWbb74JwJ49e7j77rv5zW9+w9SpU6mqquL+++/nwgsv5K677qJnz57ccsstPntt3hAoMTsopvGrrL5aM3SlApXGbMvr45z7QpvHzJ06FV54Afr1g+3bPd0spTqfhk7WDtaKuFVYWMivfvUrXnvtNS699FLmzJnD6tWrueeee/jss8/Iy8tjxowZvPLKK6xbt46TTjqJPXv2cO211/Loo4/SvXt3JkyYwG233UZ8fDw333wzS5YsoaSkhGHDhvHFF19w3HHHddhL9fU4577Q2phdva+ar1O/BuDEwhMJTQztqKYp1al8Jp95ZT8TzcSjXldjtn/UnPuv006zt95KOJRSHpeUlMSrr77KwoULWbp0KVdcccVBy7/99ltKS0t5/vnnWbp0Kaeccgr79+/nu+++47333uP5559n/Pjx1NTUEBYWRlpaGt27d2fgwIGceuqprR8JSnmchEnDrQRpvFYqkGnM9o+hFP3XpEn2dt8+21OnSbpSh+dnR+LWrVtHREQEAwYMYNKkSbz11lscf/zx/Pa3v21Yp7a2ltjYWKZOnQrAddddR21tLUVFRUyePJmYmJiGQ6RLly49aPtJSUnExcV58yWpZgSF2X4mCRZC4vVrTamj1ZoebW/QmG1pz/nhdOtm687LyrSsRakAVF9fz9y5cxsep6amMnDgQGJjY6murqa8vJwxY8awcOFCXnrpJXJzc3n88cepqanhJz/5Cddffz3Z2dl8+eWXLF68GICysrKG7WVlZXH66acfsl/lXa6ec603Vyqwacy2NDk/kpEj7W0XOktYqc7kwQcf5LzzzuPOO+/kjjvu4PXXX2f48OHs27ePe+65h759+zJ37lz+8Ic/kJmZyaBBg4iKiuKJJ54gJyeHoUOH8vLLL3PWWWcBkJOTw6xZs3jggQe45ZZbSEs75AKXystcFyGiDip36VVClQpkGrP1hNAjO+MM+PRTGDVKL0akVBf32WefNZyU5C16QujR+Sz4M6iHQU8Ootf/6hWdlVKBG7O15/xIzjzT3m7Z4tt2KKWUapGr97xmb42PW6KUUu2jyfmRXHmlvS0rsyeGKqW6pLKyMhYtWsTmzZv58ccffd0c1YSr7rx6X7WPW6KU8geBHLM1OT+SPn0gOtrenzXLt21RSvlMdHQ0M2fOJDs7m5Guc1GU3wiODAageq8m50qpwI7ZmpwfjTFj7O3bb/u2HUoppZoVHOMk5zmanCulApsm50fDVdqybp1v26GUUqpZwQk2Oa/KrvJxS5RSqn00OT8a115rL0BUU6MJulJK+aHQpFDA1pzX1+p450qpwKXJ+dGIjITzz7f3v/zSt21RSil1iJAEe2XQjDkZjeOeK6VUANLk/GhNmmRvFy3ybTuUUkodwlVzLsGCiCbnSqnApcn50brgAnv7/vuQne3btiiljtq8efO47bbbmD59OjExMfzud7/jnnvu4bzzzjuq52/dupXzzjuPF154gerqav7yl79wvutImvIbruS8tqTWxy1RSrWHxmxNzo9ejx6QkgK1tXD//b5ujVLqKOTm5vLyyy8zc+ZMZs2axaJFiwgNDeX444+ntLT0qLYxcOBAoqOjMcYQFhbG+PHjKSkp6eCWq9ZyJec5s3LY+fedPm6NUqotNGZbmpy3xsSJ9laHVFQqIOzdu5eNGzdSVFQEwKmnnsqECROIiIho1XYiIyObva/8R0iirTkvX1/O/v/u93FrlFJtoTHb0uS8Ne68097u3Qs7dvi2LUqpIxo+fDgJCQmMHTuWDz/8EIDJkyc3LH/ooYcYOHAgTz31FACrVq1i2rRp3H333fz85z8/qn3MnDmTRx55hClTpvDuu+9SX1/PlClTGDFiBAcOHODWW2+lV69ebNy4kcLCQk477TQ2btzo+RfbxYUmhzbcL19b7sOWKKXaSmO2pcl5a4waBUlJ9v5tt/m2LUr5K5GWpzlzGtebM+fw67obN675+UcQEhLCBx98QGZmJmeffTaTJ0+msLAQgD179nDjjTfywgsv8K9//QuwtY7nnnsu9957Lx999BF79uw57PbnzZtHRUUFv/vd73jkkUeYOnUqO3fu5OGHH6awsJDY2Fjuu+8+SktL6dmzJ7GxsVx88cVkZGS06nWoI2tIzoOgancVNQdqfNsgpQLEZ/JZi1POnJyG9XLm5Bx2XXfLxi1rdv6RaMy2NDlvrV/8wt4uWAD1OpauUv4uKSmJV199lYULF7J06VKuuOIKAHr27ElCQgLp6ekNh1AffPBB+vbty7PPPktdXR1VVYe/oM0bb7xB//79G/aTmZnJhx9+SEpKCqNHj+bTTz9lxYoVDBo0iPnz5/P+++9z9tlnd+wL7qJCu9nkPCjSfq2VLj+6+lSllH/RmK3Jeevdd5/tvaushGef9XVrlPI/xrQ8TZvWuN60aYdf190PPzQ//wjWrVvHtm3bAJg0aRJvvfUWn3zyCeXljWUPIkK980N79uzZLF68mOuuu46YmJijeKmG3NzchsfJycmEhtok8eqrr2bevHl8++23zJw5kxdffJEtW7YwePDgVr0GdXRcPecSao+uFC8t9mVzlAoYE83EFqe0aWkN66VNSzvsuu4yf8hsdv6RaMy2NDlvrZgYGD/e3v/3v33bFqXUYdXX1zN37tyGx6mpqQwcOJCoqKhm1581axZjx46lqKiI8vJySktLqaura3H7559/Pm+99VbD47179zLJuSbCRRddxIcffkjPnj2ZOHEiW7duJTY21kOvTDXlSs5Nrf0BV7S0yJfNUUq1gcZsK8Qnew10Tz8NY8bAV1/Btm0wYICvW6SUasGDDz7Ijz/+yKhRo8jKyuKVV17htddeY9u2baxdu5YvvviCvLw8li1bxqWXXsqUKVO46aabGDZsGK+88go33XQTq1evJigoiIsvvpgPPviArVu3sm7dOqZMmcLy5cu59tprycjI4JZbbiE1NRWA6OhoLr/8ci6++GJEhGuuuYaLLrrIx+9G5xWSGAIC9aX1pExOIfF/En3dJKVUG2jMBjGtPEzcUURkKHAzsAUYDMwxxqwUkUjg70AJkAjcbowpdp7T4jJ3mZmZZtmyZZ5t8LXXwvPPw003wezZnt22Uko5ROQHY0ymr9vRlD/G7CXJS6gtqOUne39CWGpYG1+ZUkq1nSditj+VtbwIPGCMeQSaxRHcAAAgAElEQVR4EJjnzJ8FLDHG3AG8A/zL7TmHW9axfv97ezt3Lnz/vdd2q5RSfsLvYnZ4r3DAjtailFKByp+S86GAq7inAogXkTTgKmCBM/8T4AIR6Xe4Zd5p7VDo2xfq6uCqq7yyS6WU8iN+F7NdyXn5tnIKPy4ke3a2pzatlFJe40/J+RvA0yISC1wDzAAmAvnGmHIAY0wtsB049QjLvOOZZ+ztli3wxhte261SSvkBv4vZEb3tVQQrt1ay+qzVbPntFmpLaz21eaWU8gp/Ss6nA5XA90CpMWY+kA4UNFmvBEg7wrKD5OXlkZmZ2TDNcb8QSnucfrq9MBHAjTfquOdKKY+YM2dOQ7wCkn3dnhb4Xcx29ZzX7q8l7vg4TLWhcFFhK16SUkq1nqdjtj+N1hKJrVlMBx4Vke2AwR4udRcOuC79drhlDVJSUvD4CaEur7wCw4ZBYSHMmAH//GfH7Ecp1WVMmzaNac6Y8CKS7+PmtMTvYnZ4b6fmfFcV3X/WneKlxex7ZR/df9a91dtSSqmj5emY7U895y8B840xf6Hx5KI92DP63cUCOUD2YZZ5z5AhMGWKvT97ti1xUUqpzs/vYrYrOa/cWUnK5SkQBAXvFVBz4JD8Xyml/JZfJOcikgyMMsYccGbdDxQDO4A0Z/gtRCQU6AN8Diw+zDLveuYZSEiwVy+86qpWX8VQKaUCib/G7MhBkQBUbKogvGc4CRMTMNWGvNfyPLULpZTqcH6RnAOFQKWIpLvNKwBWAwuB0515Z2B7arKNMXtbWualNjcKCoLXXoPwcDus4nPPeb0JSinlRX4Zs8PTwwmKCqImr4aa/TX0vK4nANmzsvGXa3oopdSR+EXNuTGmXkQuAv4sIj8AqcBtxphiEfkV8HcRGQX0AKa5PfVwy7zrzDNhzhxb4nLzzXD88bYWXSmlOhl/jdkSJEQdE0XpylIqNlWQMjmFnQ/tJPmCZEyNQcLEk7tTSqkO4RfJOYAxZhVwYzPz9wM3tPCcFpf5xC9+AZ98Ai++COPHw5o10L+/r1ullFIe568xO/KYSEpXllK+sZy44+LIXJWJiCblSqnA4S9lLZ3HP/8JUVFQXg7jxkFF08EJlFJKdZSojCgAyjeWA2hirpQKOJqce1psrO09DwqC/fth9Gh7FVGllFIdLnpENAClq0ob5tVV1LHrkV1sunmTr5qllFJHTZPzjnDCCbb+HGDTJsjM1AsUKaWUF8SMjQGgdHljcl5bWMu2P2wj5585lCwv8VXTlFLqqGhy3lF++Uu4/357f+VKOOUUHWJRKaU6WOTASILjg6neU03VnirAjuKSPsMOLLN5xmYduUUp5dc0Oe9Id94Jv/+9vf/VV3DttZqgK6VUBxIRYsfEAlDyQ2Mveb8/9yM0NZTir4vJfSnXV81TSqkj0uS8oz30EMyYYWvQX3gBbrhBa9CVUqoDxYyzpS0l3zcm5yFxIQz8+0AAtt6ylercap+0TSmljkSTc294/HFYuBAiI+3VRM85B4qKfN0qpZTqlOJPjAeg6IuD42zqNakknpFITX4NG6dt1PIWpZRf0uTcW845Bz76CGJi4OOPoXdvWLfO161SSqlOJ+GUBACKlhZRV9l4pFKChIxnMwiOD6amoIa6Ej2KqZTyP5qce9NJJ8HTT0NwMJSUwKhR8Oabvm6VUkp1KqHdQok+NhpTZSj57uDRWSJ6RzBmyRjGfD6GkDi/uQ6fUko10OTc2372M1i2DKKjobYWLr0Upk/XE0WVUsqDEiba3vMDiw8csixmRAwSbC9OZOoMtUW1Xm2bUkodjibnvjB6NOzcCX372sezZ0NGBuzb59t2KaVUJ5F4eiIABYsKWlynpqCG1eesZvW5q6kr1xIXpZR/0OTcV5KSYPt2uPpq+3jzZhg7FhYv9m27lFKqE0g8I5GgyCBKvi+hKqeq2XXqq+op31BO8dfFrL18LfU1erE4pZTvaXLuSyLw0kswf77tRc/OhtNOs8Mtai+6Ukq1WXBUMIlnOr3nC5rvPQ9PC+fYj44lpFsIhYsK2TB1A/W1mqArpXxLk3N/cMklsGkT3HcfhIXZk0bT0+GPf4R6/aJQSqm2SL4wGYC8t/JaXCd6SDTHvncsQdFB7Ht5H+uuWEd9tcZdpZTvaHLuL8LC4K67YPlyiI21J4s++CAkJ8Pbb/u6dUopFXC6XdANCRH2f7Kfqj3Nl7YAxE2IY9THowiODyZ/fj5rLl5DfZUm6Eop39Dk3N8MHw5799pRXERg/364+GIYMQJWrPB165RSKmCEJYfR7bxuUAe5/8k97LrxJ8QzevFoQpNDCY4JRkLFS61USqmDeSQ5F5FzPbEd5YiKgjfegPXrbbIOsHatPWH00kvtSC9KKdVGXSlm95jaA4DcF3KPeEXQ2DGxjPl6DEOeG4IENQ61qJRS3tSq5FxEzhaRJSKyRUS2OVMWML9jmtfFZWTAmjU2Ue9hv2B4800YPBhmzICsLJ82Tynl3zRmQ9K5SYSmhFK2pozir4uPuH7U4CiCo4IBqKuoY/nxy8n+Z/YRE3ullPKU1vacvwi8CdwAXOtMU535qqNceink5MDXX8OVV0JNDTz5JPTvD+PGwXff+bqFSin/1OVjdlBoED1v6AnArkd2teq5+W/lU7KshM03b2btpWup2V/TEU1USqmDSGt6A0TkW2PMcc3M72GM2evRlnlQZmamWbZsma+b4Tk//gi/+AWsXNk4b/BgeOQROPdcW6uulOoUROQHY0xmG5+rMRuo2lPFN32/wdQZjtt8HJEDIo/6ufte38fG6zdSV1xHWFoYxzx1DMnnJ3usbUqpzqU9MdultT3nvxKRyc3Mv6A9jVCtNHKkHdXlmWegp+0RYvNmOO88SE2Fhx6CykrftlEp5Q80ZgPhPcPpflV3qIddD7eu97z7Zd3JXJlJ3PFxVOdUs+aCNay7ah3VedUd1FqlVFfX2uT8ROBpEalzm+qB2Z5qkIgcLyK3ishFIpIsIpEi8oSIPCAis0Ukzm3dFpd1eiJw3XW23OXdd2HQIDs/Lw/+8Afo3RvuuENPHlWqa9OY7ejz+z4gsOfpPVRsr2jVcyP7RzJmyRgG/mMgQZFB7Ju3jwOLD3RQS5VSXV1rk/O/ANOAwcAAt+keTzRGRK4FzjPGPGyMedsYkw/MApYYY+4A3gH+5faUwy3rOs4/3/acf/stnHyyLXHJz4e//Q369bM97S+9BNXa06NUF6Mx2xE9LJrUn6diagxZd2e1+vkSLPT+XW/GrxlPv3v6kXJZSsOyiq2tS/aVUupwWltz/jJwrTGmqsn8dtcvisipwJ+AM43TKBFJA7YDicaYchEJAYqA4UB1S8uMMVnu2+50NedHYgx88w088QS88op9DBAeDhdcYC92dOyxvm2jUuqotLPmXGO2m4rtFXyX8R2m1jBu+ThiR8e2e5tla8v4ftT3dDuvG/3v60/MyBgPtFQpFah8UXP+CvBHETnFbToV+GN7GuH4B7AeeEJEPhCRE4GJQL4xphzAGFOLDe6nHmFZ1yYCJ5wAL78Mc+bAwIF2flUVvP46jBoFvXrBPffArtbVXyqlAorGbDeR/SNJ/990MLB5+mZMffuHRyz9sZSg8CAK3ilg2ahlrL1iLWUbyjzQWqVUV9Xa5PxPwG+BfwMvuE03tqcRIpIBjAHmGmNuBj4FPgB6AwVNVi8B0oD0wyxTLtdfD1u2wKZNMHWqvcARQHY23Hsv9OkDJ51ke9n37PFpU5VSHqcxu4l+9/YjrGcYxUuL2fNs+2Ne6hWpHLf1ONJ/nY6ECnmv5vH98O9Zd/U6SleVtr/BSqkup7XJ+UygtzGmnzGmvzP1Ay5pZzuGA4XGmNXO4ycBAQzQtJgvHKg5wrKD5OXlkZmZ2TDNmTOnnc0NQIMHw3PPQVERzJ8PmZlw9tkQEQFffQW//jWkpcGIEXZIRj2RVCmfmTNnTkO8Atozbp/G7CZC4kIY+A97NHHb77dRlV11hGccWXiPcAY/NpjjthxHzxt7IkHCvpf3kf1kdru3rZTyfx6M2QCEtHL9XcaYkmbm9/BAOxraYoypEJHNQDCQ2GTdWCAHG+hbWnaQlJQUulTN+eGEhMAll9gJoKQEFiyAG2+E0lJYuxZuucVO/frZix5Nngxjxuj46Up5ybRp05g2bRoAIpLfjk1pzG5G9591J/ffuRS+X8iGqRs49sNjkaD2x7eI3hFkPJVB3z/2Zdcju0i7qfGgQOHHhVTuqCT1ylSCo4PbvS+llP/wYMwGWt9zPr3pDBG5DHi0ne1YDcSLSKrbvFps0E4TkUhnX6FAH+BzYPFhlqmjFRsLV10F69bBffdB376Ny7Ky7Igv48bZoRmnT4f33oMyradUKkBozG6GiJDxbAahyaHs/2Q/ux/b7dHtR/SNYPCjg4keEg2AMYbtd2xn0w2b+DrtazbP2EzpGi15UUo1r7XJeZKI/AlARKJE5GngaeD79jTCGLMBeB+Y7Gw7AQgD/gMsBE53Vj0DmG+MyXZGGmh2WXva0mX17m1HccnKgg0b4C9/saUwAElJtkZ99myYNAni4+2QjTNnwqpVjaPBKKX8jcbsFoT3CCfjmQwAtt2+jaKvijpuZwbSZ6QTd0IcdcV1ZD+ZzbKRy1hx8gr2PLOH2qLajtu3UirgtGooRbAXnACuBs4C9gK/AIqMMfvb1RCRZOAxYBm2N2WuMWadiCQCfweysIdi/2iMKXWe0+Iyd11uKEVP2rkTUlJsucu778KTT8L+Jh91t25w7rm2hv3MM+36SimPaO+wXBqzD2/zbzeT/Vg2oamhjFs2joheER26v9JVpeT8K4fcF3OpK60DIOPZDHpe27ND96uU8g5PDKV42ORcRE5pYdEE4KfAw0A9cJEx5nftaUhH0uTcg158EV54Ab74AmoOOY/LGjkSJk6En/4UTjnFJu9KqTZpTaDXmN169bX1rD5rNQf+e4CYcTGM+XyMV2rCa0tqyXs9j32v7mP468MJibMl/Nv/tJ2a/TWkXJpC/MnxBIW09gC3UsqXvJGcfwMMBQ5gA3qzqwE9jTHh7WlIR9LkvAOUlcF//2tPJn3nHdi3DwYNgt27obLy4HVHjoTTTrMJ+ymn2DIZpdRRaWVyrjG7DWoKavhh/A9Ubq8k8X8SGblgJEFh3k+K66vr+brn19QW2jKX0ORQki9OJuXSFBJOSyAoVBN1pfydN5Lzc4Flxph9R2jIZGPMG+1pSEfS5LyDGWNrz7t3t4n3d9/B3XfDZ581v/7IkfCTnzROAwfqSDBKtaCVybnG7DYq31TOipNWUJNXQ8rlKQx7eRgS7N24ZIyhdEUpea/nkfdGHhVbGkeeDEkIIeOZDFIu0bJBpfxZhyfnnYUm5z7w3nu2BObjj6Gg6XVHmkhJOThZHzcOIiO9006l/JwnAn2g8VXMLllRwsqJK6krriP1F6lkPJPhs7ISYwxlP5aRNz+PvPl5lK8tZ+y3Y4mbEAfAvlf3Ubmjkm6TuhE1LArRDg6l/IIm50dJk3Mfqq+HFStskv7RR/aCRxMnwhlnwNdfw5IlkN9kSNDQUBg9GiZMsBdLysyEoUMhWMcGVl2PJufedeCLA6w+ZzX15fWkTE5h6H+G+qTEpamyDWVEHRPVMB77iokrKPrcjjAT3jecpLOTSDw9kYSfJhCWHObLpirVpWlyfpQ0OfcjZWX2KqVpzsU55s+3FzlyCQ6GurpDnxcVBWPH2kR9/Hh7O2gQBPn+S1OpjqTJufcVfVXE6nNXU1dcR9I5SQx7bRghMa29Zl/Hynszj/x38yl8v5CafW4n5wv0ub0PA/42wHeNU6oL80TM9q9oozq/6Gg7uZxyii1/WbzYnmCalXXw+uefD6tXw44dtpd9yZLGZfHxtgRm7FgYNcr2tmdk2J53pZRqo/gT4xm9eDSr/mcVhe8XsvLklYxYMKLDh1lsjZRLUki5JAVTbyhZVsL+T/az/9P9FH1VROTgxrLA/IX57HpoFwkTE4g/KZ64E+IaRoZRSvmndveci8gY4P8BK7AXlMjyQLs8yte9MKoVsrJsor54sS13ee89O3/fPhgxwo4EExcHpaW2B76p8HC7nitZHz0ajj3WJvJKBSBP95xrzD565RvL+fG8H6nYUkFYjzBGvDuCuPFxvm7WYdVV1IGB4ChbBrh5xmayn3S7zlMQxIyKIf6keBJOTSDlUj3BVClP8ouyFhH5F3ChMaaHiFwFvGWMqTjS87zJXwK9aoe8PJt072syCEVyMqSn21FisrJg+/bmn9+/v03UXUn7scdC375aFqP8Xgck5xqzW6GmoIa1k9dy4LMDSLgw6NFBpN2YFjAnYNYcqKHoiyKKviyiaEkRJctKMLX2ez9mTAyZy+2fljGGPXP2EDMuhphjY/yizl6pQOQvyXko0N0Yky0iQUCMMaa4XRv1MH8K9KodjIH16235y+LFdqjGwkK7bPFie6JpURE895y9omllpV1/zRqoqjp0e9HRMGyYTfpHjIDhw+1tWpoO7aj8Rgck5xqzW6m+up4tv9lCzlM5AKT8LIWMORkBWR5SV15HyfclFC0pIqRbCOk3pQNQvqWc7wZ/B4CECTFjYoibEEfshFjiJsQROSiy4WRUpVTLfJKci8g5QJIx5j8iMgEoNcasa08jOpq/BXrlIfX18OOPNlmfPt2WtACcey68/77tFc/MtBdAGjzYJtzr18PKlTZh37On+e0mJDQm6u6Je4oe/lXe195ArzHbc3Ln5bJp2ibqSuuIGBjB0H8PJf4nnaNkrnxTOTv+uoOS70oo31B+yHL3YRxLfywlOCqYiAERAXMEQSlv8XpyLiJ/BX4DfGCMmezMexBYbYx5uT0N6Uj+GuhVB3n4YXj3XfjmG6hxG8UgIgJuvx3uucc+LiiwPexr19pk3TW5euOb6t7dJurDhtkTT4cMsVN6uva0qw7TnkCvMdvzyjeVs/bytZStKgOBXrf0ov99/QmO7DxDvdYW1VKyrITi74op+a6E0pWlTNgwgaBwW+qy8oyVHPj0AMFxwcSMjiFmTAyxY2KJGRND1NAovZKp6tJ8MVrLT4BewBS3ee8DLwJ+G+hVF3PrrXYqLYUvv4RPP4VPPrFXMXXv/d64EZ5/3va0P/CAPWnUGMjNbUzUXYn72rW23v2//7WTu+joxmTdPWkfPFgvpqR8TWO2h0UdE8W4b8eRdW8WOx/aye7/203BggIyns4g4eQEXzfPI0LiQ0g8PZHE0xObXR6WGkZoaig1uU49+xeNJ+f3vL4nGXMzAKjZX0PZmjKiR0QTmqijaCl1tFqbnH9sjDkgIu7d7ZNbXFspX4qJgXPOsRPYk0rdh1l8801bn/7ccxASAiedZBP1SZPg9NPthZJcjIFdu2yivn69Tew3bLBTXh4sX24ndyL2pFP3pD0jwybtaWl6MqryBo3ZHSAoPIgBDwwg+aJkNkzdQPn6claespLUa1IZMHMA4T3Cfd3EDjXsP8MAqNpbRemK0oapZEUJMaNjGtbb/+l+1l1mK6jC0sOIHhF90BQ7JhYJ1qOOSjXV2rKWK4ABQDrwJXAZcDEwwxjzzw5poQf48yFS5UMbNsDbb9vhGr/++uCLH/30p4f2kLeksPDgZH3DBvt4y5bmL6gEtkd94EB7IaXBg+2t6356uibuqkE7y1o0Znewuso6dv5tJzsf2ompMgTHBdPv3n6k/296lyzvMMY01KHnvZnHzr/tpGxtGfUV9QevGAwnl55McIQtB9rzzB6CY4KJHhFN5OBIHS1GBSxfnRB6PDAV6APsAeYZYz5pTyM6WiAFeuUj+/fDRx/ZRH3RIrj6anjsMbssPx9mzYLrroNevY5+m9XVsG3bwYm7K2nPy2v5eeHhNnF3T9pdiXuvXvYqqqrL8MAJoRqzvaBiawWbf7OZwkX2nJXIQZH0f6A/KZNTuvxJk6beULm9krI1ZQ1TXVkdI98d2bDOV92/oibPOUcoGCIHRBI1JIqoIVEkX5xM/Amd48Rb1fn5y1CKwcBEY8yn7dpQBwrEQK98qLYWysoaL1w0dy5Mm2aT4vPPh//9X1v20p4v3AMHYOtWm6hv2QKbNzfez81t+XlhYTBggE3WBw609/v3b7yNimp7m5Rf6oChFDVmd6D8BflsvXUrFZvs0PGx42MZ8NAAEn/afP22gvraerLuzrKJ+9oyKrdXgltH+6AnBtHrZtsxkr8gn10P72pI3F1TRN8IHepR+YUOT85F5F0gpsUVrGQg3xhzWnsa0pECOdArP/D993YEmDfftIk7wCmnwEMPwfHHe35/xcUtJ+4tDf/o0qOHTdTdJ1fyrnXuAak1gV5jtn+or6ln73N7yboni+o91QAkTEyg75/6kvDThC7fk34kdZV1VGypoHxDOeUbykm+MJmYkfbPOuveLLLuyTrkOUERQUSPjGbst2Mb3t+SFSWE9QgjrEeYvufKa7yRnN+PPWl0I/Z37A3YSz67R83+2DF0f92ehnSkQA/0yk/s3QtPPw2PPGLrzIOC4N574c47vTeUYmlpY+K+fbstm3FNWVkHDx3ZVHg49Ot3aPLuSuBjY73zGlSrtDI515jtR+rK6tj92G52/n0ndUX2/JO4E+Loe1dfks5J0oSxDar3VVO6qrQhcS/fUE75+nKq91QTNSSKCesnNKy7JGkJtftrCYoOInJQJFGDo4gcHEnk4EgSTkkgcqCOpqU8zxvJebyzzgHn8V+MMX9uZr13jDEXtqchHamzBHrlJ4qK4K9/hccfhyefhOuv93WLrLo6yM62iXrTxH3bNjsU5OEkJ7ecuPfqZUe0UV7XyuRcY7Yfqi2qJfvJbHY9sovaAnv0LWZMDL1+14vuP+uuJz96QG1RLdX7qokabEv76irqWPnTlVRsrqC2sPaQ9QfPHtxwddSC9wrY+/zehsQ9clAkkQMibY+7lsqoVvLFRYj+BUw3xtS5zZuIPcGoZ3sa0pE6W6BXfqKgALp183Urjl5p6cFJe9MEvqqq5eeGhNhhIVsqmUnUetqO0s7RWjRm+5Ha0lpynsph18O7qMm1R7nCeoSR9r9ppN2YRlhKmI9b2DnVFNZQsaWCis0Vtlxmczm9ftOLuPH2iqfb/7ydHfftOOR5QRFBRA2PInNZ47/f/s/2E5YSRkT/CIKj9OR8dShfXIToDWCViCwBSoChwFnAIT0zSnVqZWW2zCWQkvOYGBg50k5N1dfb19O0t337dltGs2ePvd26tfltJyS03Ovep489kVX5gsZsPxISE0KfW/uQfnM6+/6zj92P7qZsTRlZf8pix/07SL0mlfRfpRM7TkvMPCk0KZTQCaHETYhrdnnqNalEHhNpk/fNFVRsq6ByWyU1eTXUVzaemWqMYc0Fa6grsb91w3qEETEggsgBkUQMiCD5wmRix+pnp9qvLUMppgHTgWOBYuB1Y8w7HdA2j+msvTDKBwoK4IUX4B//sCdufvRRx5wU6m8qKmxNe3PJ+7Zt9sdKS4KCoHfv5nvcBwyw5TRae9siDwylqDHbTxljOPDfA+x+dDcFCwsa5seMjSFtWhrdr+xOSJyWk/lKbUktNQU1RPaztel1FXWsvWStTd63V2JqDs6fjpl7DGnXpwGw96W97PzrTiIGNibvDbf9IwmO1l73zspfhlIMAsYYY35o14bstsKA74HfGGM+E5FI4O/YHp9E4HZjTLGzbovLmuoqgV51kF274OOP7WgtH37YOGLL2LH2BNExY3zbPl8zxo7b3lLivmuXXaclMTEtl8v06wcREV57Kf6oA4ZS9FjMdrbn8bjdFWN2+aZycmbnsPffextqpIOiguh+ZXfSbkgjdkKsnkDqR0ydoSqnisptlQ097Sk/SyFmhB1VZttd29j5153NPjc4JpiTik9q+Dx3P7ab4JhgIgZEENE/gvBe4QSF6HkIgcoXNefPNjM7BagzxlzUnoY4278T+D1woRPknwM+MMa8KiJnA1OMMVc667a4rKmuGOhVG9XX25Mqe/e2jysrbcmGqx47KAjOOAN+/Ws491zt8T0aVVWwc2fzifvWrfYIxOGkp7fc696jR6f/DNpZc96hMdvZh8fjdleO2XWVdeS/mU/OnByKPi9qmB95TCSp16SSenUqkQN0lBF/V1dWR8XWxhKZim0VVG6tpGJ7BSFxIYz7bhxgj558GfMl9eWN5TMSIoT3CSdyQCTpN6eTfGEyYHvyTbUhJClEf6j5MV8k518CnwLuTzoR+NIYc1+7GiLyE2AIthZyKrAJ2A4kGmPKRSQEKAKGA9UtLTPGZDXddlcO9OowqqrsVTtXr4Yff4RVq+yY5lVVdkQW1+gkl15qe8vPOsveT031bbs7E2Ps1VmbS9y3bYMdO+woNC2JjDw4WW96Pzrae6+lg7QzOe+wmO1sv0PitsZsq3xjOTlzc8h9KbfhBFKAuBPjSL0mle6XdSe0W6gPW6jawhjTkFzXV9ez44EdVG6rpHK7TeKrc6ob1s14NoOe19pzt3Pm5rBp2iaC44JtiUz/iMZymf4RJJ2tw3P6A1+cEHq9MWZjk0YMwQblNhORaGCyMeYWEXGdqDQRe6GMcgBjTK2IbAdOBWoOsyyrPW1RndD+/fZCPvHxkJFh573zjk20m0v8evaE3bttSQXA/Plea2qXIwJJSXbKbCaW1dbaspiWRpgpKIB16+zUnNTU5hP3gQNtj3zn/yLrkJjtbEfjdgeLyohi0MODGPDgAA58eoC9L+4l/618ir8qpvirYrb8egtJZyeRclkK3c7vRmiCJuqBwD2BDgoLov89/Q9aXldRR+WOSiq3VRJ9bGMHQ11ZHcExwdQV11G6spTSlaUNy4Ljgzn5wMkNj9f/fD0EcUgSr8NDBoZWJedNg7yjCrge+EM72nE78GCTeelAQZN5JUAaUHuYZYfIy8sj0+2Lf9q0aUybNmeDppMAACAASURBVK0dzVV+a+FC2xO+aZNNyDdtgvx8u+zmm+GJJ+z9fv1sr+0xxzSOYDJyJIwfb8f07vxJW2AICbFJdf/+cPrphy4vKmpM2Ju7KFNurp2WLj30udHR9sfakCGN09ChMGiQz+vc58yZw5w5c1wPk9u6nQ6M2dCBcVtj9sGCQoJIOiuJpLOSqC2tJf/tfHJfzGX/J/spWFBAwYICJFRIPDORlMtSSL4gmdAkTdQDVXBkMNFDookecvCRv96/7U2v3/SipqDmoJ72yu2V4Faibowh7428g0aacZFwYeDMgfSa0QuAiu0VlK4sJaK/Td71BOS28VTMdmnVp+D0crgfHg0GUoEX29oAETkLWGaMaXqFFANUNJkXju194QjLDpKSkoIeIg1g9fX2Ajo7dtiEq+m0YkVjMjVzJnzxxcHPj4qCwYNtfbLLiBF23O9Ird0MaPHxMHq0nZqqq4OcnOYT9y1b7Emsy5fbyV1QkP0x4ErYhw2z2x8+3F5l1Qvck1ERyW/rdjoiZjvb7dC4rTG7ZSExIfS4pgc9rulB1d4q8t/MJ++NPA58foDC9wopfK+QTSGbSDg9gZTJKSSfn0xYqg5l2lmICGHJYYQlh7U4NCT1MHLRyIbEvaHmfbsdHjIkvjH1K/ygkM3TNzc8DukWQmT/xtFl+t/fHwm2nVWm3mivews8FbNdWvsT6WPgZRqDfT2Qa4zZ1I423AqMcTvMkwi8A/zDue8uFshx9t/SMhVI6uttkrR7t5127bK3J5wAFzoXMFy4sPF+c3butD3gAJddZnu/Bw+28445BtLSDu0JDw7WxLyzCw62J/b27g2nnnro8sJC2LjRnnfgPrnGc9+6FRYtalw/JMT2rI8ebUfocf0o8O8LMHVEzAaN234hvEc46dPTSZ+eTvW+avLfzifv9Tz2L97P/g/ttEk2EXdcHN0u6Ea387sRPTxa65I7OQkWEk9LJPG0Q2NTbWntQQl2eFo4SZOSGnriawtqKSkooWRZCSFJIQz424CGdb8f+T315fUNo8q4SmYiB0QSmRGpZVUe1NoTQk8xxnzRZF4G9sz/LW1qgEgKtvfEZSlwC/AtsA5IMcZUiEgoUAoMAOqALc0tM8ZkN92HnlzkI64eb1finZsLN97YuPy00+Crr6C6+tDn3nQTzJ5t7y9fbkdI6d/flqM0nTIy9CI3ynOqq21ivmEDrF9vTxZeudIm8s3Fy6FD4cQT4aST7O3AgR4ti2rnCaEej9nONjo0bmvMbp/q/GoK3ikgb34e+/+7H1PV+Hcb0S+iIVFPOCWBoDAdsk9Zxhiq91bbnvbtFdSX15M2zVadmXrDl9FfNlsqA9Dvvn70u6sfAMXLitn3n31EZkQSlRFF1JAoW+veRX4U+uKE0LOBJjUD5AKLgTYN9myMyXN/LCJ1QJ4xZqeILAROBxYCZwDzXUH8cMuUF9TW2qtG7t5te6b79rXzP/4Y7r3Xzs/JgZomR6x//nNbZgI20amuticD9uplezhdt+4X9hk71vZyKuUNYWE24R46FC6+uHF+WRmsWWMT9ZUrbTnVqlU2gV+/3o55D/Zv+Lzz7HTaab4+QuPxmA0at/1dWHIYPX/Zk56/7EltaS37P3Zq0xcWUJlVSfbj2WQ/nk1wXDBJZyeRdI6tZw/v6Z2yLeWfRITwnuGE9wwn/sT4g5cFCScdOInKnQeXybiGiowe2lgfX/xNMbsf3X3Q84PjghsS9YxnMggKtT8KTZ1pKJtRjY6q51xEbgbuxx6CbM7HxpizPdIgkSxgqjNebiL2ghVZQA/gj8aYUme9Fpc1pb0wrVRZaUtNXGN9A9x+u63TdfWC791re8YB7r8f7rzT3l+wAC64oPF53brZUTFc5QUPPNBYBpCbC7Gxjcm6UoGmutoe2fnqK1iyxN7mueWt8fFw7bUwY4YdJaYN2tIL482Y7ewvCw/GbY3ZHcPUGYq/K6ZgQQH57+ZTvrb8oOXRo6Jtsn5WEvEnxmuvumqT0lWlFLxfQMXGCso3llO+oZza/fbCWqHdQzkx98SGdb8d8i2mxhA1JIroYdFEHxtNzLExRA2JIig8MP/+vDrOuYgcC0wBHm+yqNoYs6c9jehoGujdGNN4yH35cnjvPXvRHVfSvXu3Hd0kKckOU+cyeLBNzl1E7AmWvXrBddfZMhSwz1mzxs5PS/N1r6FS3lVfb3vUFy6Ed99tPNk0IgJmzbKJeiu1NdBrzFZHUrGtgoJFBRR+UMiBxQeor2gsWQiOCSbhtISGZF0vfKTayhhDTX6NTdIP1JJ8vh3MpL6mni9jvzyo7MpFQoSBjwyk1812VJmaghrqyusI7xXu9+UxvrgIUagx5pAz60WkjzGm+evU+oEuFejr6+0l5psm3K7pscdgyhS77uzZMH36odsICbG93Zs2NdZyv/aave3Vy049e0Konvyh1GEtX25HEHrlFft/tXOn/d9phXbWnGvMVkelrrKOoiVFFH5QSOEHhYf0qkcOjmzsVT81npAYHXJPtV99VT0VWyso31BO2Zqy/9/evcfJUdV5H//8JpnJJDPJzCSZkMkdngDhEgUTQGWByE0FbyCXXXe9S/BZeIRnxcULuqywj5cF1MVrXFlwFdZVVkEUQQREUMAgoiCXEAjknkkyQ2Zyn8zv+eNUZ2p6uufWPVPV3d/361Wv7q6qrjrnVPfpX586dYrtf95O55862bliJ0fcegTNZzUDsOYra3j+0ucZO2Usk46dxMRjJu5/rJmWrmvORrzPuZldBvzc3TN3+PhHMxuTvRrwBsLNJ2SkdHeHriTZd1JcvRrq60MrHYQW7bPPDl1Tclkb69553HGhu0om4M5M06aF4eTizjtvZPIlUq66ukI/9cx3ad++vtdgFJnqbBmuMbVjmHzqZCafOhmugV2rd9F2d1sI1n+5lZ0rdrJ2xVrWXr8WqzYa/qqBptObmHz6ZOqPqtcQezIsVeOqQneWw+toPrt5//x9O/b1Gru9e283Y6eMpWtLF1vv3MrWO3uuQ6tfVM/i5T2xcPwOrKVqoL++BwKNsdevAeqADfQMzWWEvoNSqN27Q9eRlStD4H3GGaE7CcCnPx36a+cyeXLPc7MQSFdV9Q26Z83qve5rXhMmESncjh2hpfyRR8L0q1/1XMhcXQ3XXQdz5ox0KlRnS1HUzq7df1Fpd1c3HY90hED97q10/L6D9vvaab+vnRc/8SLVzdU0nRYC9abTmhg3QxeWSmHGTOjdpjDnsjnM/uhsdq/ezbZHt9HxaEd4XN5B9dSes/j7du3j0UMeZdJxk2h8Q+iWVYpdsvoNzt39oqxZ/wI86e69xr4zs8OLnbCKsGsXXHllz2gPK1f2XGQJ4WKyTHA+bx40N/fcfjxzO/I5c0IXlLibbhqtHIhUFvdwTcbKlT3f28z04ou9v78Q7jb6zneGazLmzRuF5KnOluKrGltFw/ENNBzfwIFXHcjerXtp+1UbW+/aSttdbexes5tNN29i083hnlR1C+v2t6o3nNDAmPHZJ29Ehs7MqJ1TS+2cWqadMw0IFznvbes5I9nx+w52r95N6+pWWn8ULs6vW1jH1LdPZcrbpzBx0cSSaFUfap/z97r7TVnzDgUOcfefFjtxxZJo/8XWVnjssZ7JDG69NSxzD6OVbN8eXmfuTHjwwSH4fte7wrjJmXVL4AMlUtL27Andx9avD/3Dc92RdseO3O8dMwYWLoRjjw3TCSf03ByrAAX2OVedLSPK3dnx7A7a7mpj691bab+/ne4dPX9Sq2qraDixIbSqv7FJN0GSEeXu7Hx+J+33tYc/kHduZV/Hvv3LF/9xMfWvrh/RNIzaOOdmdgLhVOhp0e2g4yYTRgNIbUU/6u6/H771rTCs2urVvZfV1oa+qGPHhmD7mmvC0IKHHRZ+yDO3oc+mykxkeNzDH+DW1jCtXx/G4F+3rvfzdet6D4OYT2Nj+BO9YEHPeOiHHRZaycel43S+6mwZLWZG3YI66hbUMeuSWXTv7uaVh15h691babu7jc7HO2m7u422u9vgMqhpqQmt6m+cTNOpTdQ0p+tiPiltZsaEgycw4eAJzFg6g+7d3bTf387m2zaz68VdvQLzNdevYfKbJzNhfvqGcx7s5dbtwHcJ/RmPz1q2B/j3YiaqpOzcGW68M28evOpVYd7LL4fRGQDq6sKtvhct6pniF1tmhiAUkcFxh46OMGzn5s3hLrStrX0f48937hzctquqwhChLS3hGo3su9LOnRuC8/RTnS2JqBpX1XPr+M/Dno17aLunbX+wvmf9HjbetJGNN20Eg/rX1IdW9dObaHi9xlaX4qoaV8XkN4ZRhuI6/tjB8x95HqpgxtIZHPT5gxjbkJ4RiIYyznkDcK67l1ylXvRTpO7hhiPf+EYYJWX79t63m1+/Hn7yEzjxxNC6Nkb97UT62LcPXnklXDQ51GnfvoG3H1dbG0Yham4OwfeMGT1TS0vP82nTUvd9LWCcc9XZkiruzvY/b98fqLc/0N5rjOuquiqa3tC0v7/6+EPGqwuMjIhdL+9i1T+vYuN3N+JdTk1LDYf/4HAaTyi88WXUxznvJyF/5e4PFryhEVK0it4dfvzjcBHnn//cM3/RonAjnlxjhouUq66uEFy3t/c8xp/399jWFp4Pt/6ZODF0B2tu7pkywXeu53V1Jds1rBgVfY5tVkadLam2b8c+XvlN6AKz9a6+Y6uPmztuf6t60ylNVDfp3hpSXJ1PdvLc0ufY9rttWLWx4KYFHPA3BxS0zREPzs3sdmCgnvMGHOHu0wpJyEgqWkV/6aXhJj4QfvAvuAA++MFw6lukVLiHYTu3bQtTR0fP88w0mAA7cyFzIRobw/CeQ5mamnpujlUBhlLRq86WUrZ77W62/jK0qrf9so29m2P3BaiCScdO2t9ffeKxE6kaqy4wUjjf56z82ErWfGkNNtZY+POFTD5t8sBvzGM0Lgh9ltB3cS09Y+T2SQfwlkISUTIuvTTc7v7SS0NQnpKLv6RCdHX1BNLZAXV/r3MtK8bNcKqqwnCfjY1De2xoCEF2Y2PqupCUAdXZUrLGzRxHy/taaHlfC97tdD7eub9VfdtD29j2cJhe+uxLjJk0hqZTmpj85slMfcdUXVgqw2ZjjPnXzceqjfb72qlpSf6zNFDL+VSgw91397sRsynuvqXYiSuWglphtmyBKVN6XmtIQxmK7u7QwhwPkvMFzgMF2PmG8BuO6uoQJE+aFLqITJrUM02cOLhAu75e34VRMMSWc9XZUpa6Orpo/3X7/iEbdz4Xu8i7ChqXNNJ8bjPNZzVTc0DywZWUHu928BCsF2LEW87dffMgt3MW5Xj1/9q18OpXw7vfDV/8YghoFIyUv8zQe9nBdK5AeqB1OjuH3686m1n+YHqor3XWpyxVfJ0tZWvsxLFMfctUpr5lKgA7V+2k7e42Nv9kM233tNF+bzvt97az4qIVNJ7YSPP5zUw7f5r6qcugWVV64ruCxo0xsznARcDFlGNFf+WVoeX8mWd0+j3t3MNweYMJmAdap7Oz750eCzFhQt9AOfsxXyAdn1fCFzVKOpR9nS0VY/y88YxfOp4ZS2ewt20vW27fwqYfbgqjwNzfTvv97Tx/6fM0n93M9PdPp+mUplQFX5JeW+/eysbvbaTlQy00npjM0LlDDs7NbAzwNuAC4DRgI6GPY3nZvBluvDH0q/3KV3qPTS7Fke/CxOEG10MdXq8/48f3H0zHg+j+1qmvDzecEklIxdTZUrGqm6qZ/t7pTH/vdPa2h0B9439upO1XbWy6ZRObbtnEuDnjmPn3M2lZ2qLWdOlX+/3tbPzPjdS01KQ/ODezA4EPAe8HpgAvAW8F7gJeNyKpS9Jtt4UL8E4/vSi34C4r7qH/c2Y4vMzoHcMJrru6ipeu2trCg+nMowJqKXEVV2eLANWN1Ux/z3Smv2c6u17axYabNrDhPzawa9UuXvj4C6y6ahUtH2xh9j/MpnZunjtyS0WbuHgiANv/XIQRyYZpwAjEzM4ltLicDDwPXAfcBHzS3X8RrfbbEUthUh54IDy+9a3JpmOk7NnTO7gezPP462IF1TU1hQXR8XWr1RoiUrF1tkiW2rm1zPvMPOZeMZetd29lzbVraLunjbX/tpZ131zHzItmMvdTc6meot8O6VF3RB0A2/+S4uAc2Em43fNvgQ+5+3MAZlakq9xS6rHHwuOxxyabjoF0d4eAOXMr8/iUPS8eXA/2dub5jB8fxpuOD4832KA681wXJoqMhMqss0XysCpjypumMOVNU+h8opOXv/Aym27ZxJovrWH9DeuZf+18pn9guu5GKgCMmxXikj0b9uDuiXwuBgzO3f0O4I7oQqILzGwecBuw/6+mmZkX41ajaeEOK1eG54cdNvr77+yE9ethw4bej62tfQPvrVuHd/Hi2LE9wXXmMTPFX+da1tCgoFokpSqyzhYZpPpX13P4zYcz+2OzeeHyF2j7ZRvPfuhZWm9t5bD/PEyt6MKYujFU1VXRvb2bfR37GDtp9Lu5DnqP7v4y8GkzGwu8A5hlZlcB3wFOB5aNTBIT8oc/wLp1oXW3WLq7YeNGeOklWLUqPK5b1zsAX79+6HdebGyEqVN7T1Om9H2dufFLY2MYQUStBCJlq+LqbJEhmHj0RF5116vYdMsmVly8gq13buUPr/0DC3++kAkHT0g6eZKwmuYadm3fxd7WvekOzjPcvQv4EfAjM5sP/F9C/8byqejNQov5cFrNu7rghRfg6afDtHJlTzD+8sthdJKB1NZCSwtMn977cdq0vkH45Mm6eFFE8qqIOltkGMyMA951AA0nNvDkW5+k84+dPHHqE7zmd69h3AydHa5kk143iXFzx0FCo2gXFNW5+/PAJWb2WJHSU1o6OmD5cnj44dDS/vTTsGJFuNgynylTYN48mDs3TLNmhcA7HoRPmqRWbREpuoqvs0VyqJ1Vy1G/OYo/nfYntj28jafOe4qjHzha46JXsMNvPjzR/RelydXdv1voNszsbOBaYBLwPeCj7t4V3Y76Xwjj8lYDl7v73ug9eZcV5Kmn4Lrr4Oij4eKLe+bv2QMPPgh33AH33ANPPpn77o+zZ8Phh4eW94MP7h2M19cXnDwRkUKUXZ0tUqCx9WM58qdHsnzhcrY9tI0N/7GBlg+2JJ0sqVCp6A8RXbj0DuAc4AjgG8Bq4Brgh8DH3f0RM/sw8FngE9Fb+1s2fC+8ADfcAGeeGYLzlSvha18LNyVqa+tZb+xYOOooOO44OOYYOOIIWLBAAbiIlLXU1dkiRVAztYaD/vUgnnn3M6y+ZrVGcKlg3V3ddO/spmpcFVU1o38TylQE58Bc4ANR38jHzOxI4A1m9iBwqLs/Eq33U+BZM7saWJhvmbsXNjjljh3hsbYWLr88tKJnxvU+/PAw9vkZZ4SAfPz4gnYlIlKC0lVnixTJtPOnsfKylex4Zgfbn9pO/ZFqbKtET73zKbbcvoUjfnwEze9oHvX9pyI4d/ffZM1aC0wk3ETjxdh6a6N/sYuB4/tZ9uuCEpQJxB96CG69NfT/fs974CMfgUWLCtq0iEipS12dLVIkVdVVTD5tMhu/t5FXHnhFwXmF2n+9wTBGqi6G0W+rH5xjgG8CM4EtWcs6gBkDLOultbWVxYsX75+WLRtgkIJ9+8Ljhg1hhJSHHoKbblJgLiKjYtmyZfvrK2Bq0ukZhGTrbJEiqjsy3CFy54oCb9YnpSsapcX3De52EMWus1PRch4XDfXV6u5PRHe0y/52jAP2Av0t66W5uZnly5cPPhHx0VZuvhle97rBv1dEpEBLly5l6dKlAJjZ5oST069U1NkiRTRubhhGcdfqXQmnRJIy1JbzYtfZqWo5N7MxwIXAx6NZa4Gm2HID6oB1AywrzNq14bGuDk4+ueDNiYiUo9TU2SJFlLnpzL7OfQmnRJJiY0JwPtiW82JLVXAOXAZc4+6ZO/X8Ejg0tnwesAdYPsCywuyK/i3PmKHxxkVE8ktHnS1SRGPqQp+G7h0JdTiW5GW6tXRXeHBuZp8EHgfGm9lBZvYBwli466ORAADeBFzv7nvcfXm+ZQUnJnNn0GOPLXhTIiLlKFV1tkgRVY0PoVH3TgXnlWp/t5aETp6kos+5mV0BXJU1+xl3v8HMzgc+Y2argHrgU7F1+ls2fNujUb00XrmISB+pq7NFiilqtkyq1VSSN+PCGTSd3kTD6xsS2X8qgnN3vxq4Os+yl4APDnVZQTKjtVSl5sSCiEhqpK7OFimi/TceUmxesRqOb6Dh+GQCc0hRt5ZUuffe8PjQQ8mmQ0REREaXWs4lYQrOc/HoC6mWcxERkcqSGQdCsXnFaru3jdVfWk3nE52J7F/RZy6ub6SIiEglSvrukJK81h+2svIfVvLKb19JZP8KznPJBOcaRlFERKSyZLq1qKFOEqLgPBcF5yIiIhVp/wWhajmXhP6fKTgXERERychERmo4r1wJt80qOBcRERHJyIykqNFaRC3nKTJ/fng89ND+1xMREZGyYurSKmo5T6EZM3o/ioiIiEhFqKqpompCFTY2mSg9FXcIFRERERFJg/lfms/8L81PbP9qOc9l3brejyIiIlJZ1OVcEqLgPJcVK8Ljc88lmw4REREZXepyLglTcC4iIiIiEll19SoeWfAIG763IZH9KzjPRXcFExERqWwKBSrW3o172fnsTrrauhLZv4JzERERkQx1a5EMjXOeQhrrVERERKSyaJxzEREREREBBeciIiIifbiuPxN1a0mRk04Kj0uWJJoMERERGWXq0SoJfwZ0h1ARERERkUjjkkasyqg/qj6R/Ss4FxERERGJNJ/VTPNZzYntv+SDczMbD3wR6ACagMvdfVtBG12+PDw+9liBqRMRkbgRqbNFRoK6nEtCyqHP+deBB939k8BtwLcK3mJnZ+9HEREpluLX2SJFZBpGueLtXLmT9gfa2bVmVyL7L+ng3MxmAO8CfhrNugd4m5nNK2jDukJbRKToRqzOFhEporVfXcsfT/ojrT9sTWT/JR2cA0uAze6+A8Ddu4AXgZOSTJSIiOS0BNXZUirUTicJKfXgfCawJWteBzCjKFvXqS0RkWIa2TpbpBj00y8ZCf1BK/ULQh3YmTVvHLA3PqO1tZXFixfvf7106VKWLl068qkTERmGZcuWsWzZsszLqUmmpchUZ4tI+g3xD1qx6+xSD87XEq72j5sIrIvPaG5uZnlmBBYRkZSLB6Nmtjnh5BST6mwRKR2DbDkvdp1d6t1a7gNmRENzYWbVwBzg1wVttaWl96OIiBTDyNTZIiNBfc4rV8Jdm0o6OHf3DcAdwCnRrFOBW919bUEbPvTQ8HjIIQVtRkREeoxYnS1STOpzLgkr9W4tAP8b+KKZvRqYDqhjoohIeqnOFpFUm/OPc2i5oIWa5ppE9l/ywbm7twEXFHWj26Kb1XV0FHWzIiKVbkTqbBGRIqo5oIaaA5IJzKHEu7WMmD/8ofejiIiIVBTXDQklIQrOc9EXUkREpDKpz3nF23jzRv7yd39hy8+yb8swOhSci4iIiIhEOh/vZNP3N7H9L9sT2b+CcxEREZFsOokuCX0GFJyLiIiIRMzUr6XiaZxzEREREZGUUcu5iIiIiEjC1HKeQosW9X4UERGRyqI+55IQBee5TJzY+1FEREQqg7qcV7zaA2tpOKGBcbPGJbL/kr9DqIiIiIhIscz88ExmfnhmYvtXy3ku9fXQ0qKWcxERkQpjY4yalhqqp1UnnRSpUGo5z+Wii8IkIiIiFWXczHG8ft3rk06GVDC1nIuIiIiIpISCcxERERGRlFBwLiIiIiKSEgrORURERERSQsF5P5YtW5Z0ElJF5dGXyqQvlUlfKpPRUSnlrHyWF+Wz7EwtdAMKzvtRQR+kQVF59KUy6Utl0pfKZHRUSjkrn+VF+Sw7zYVuQMG5iIiIiEhKmLsnnYYRZ2atwEvDeOtUYHORk1PKVB59qUz6Upn0VUiZzHX3gltiSonq7AEpn+VF+Swvh7p7QXexrIjgXERERESkFKhbi4iIiIhISig4FxERERFJCQXnIiIiIiIpoeA8BzMbb2bXm9n/M7NvmNmkpNM0UszsNDN7xMzmxeblzf9wl5UKMzvbzF40sy1m9hUzGxvNn2pm3zKzL5jZdWZWHXvPsJaVCjM7ycyeNbM2M/tybH7Ffk4AzKzGzJ4wsyXR64oujySVcxma2dVm5tH0RDSvLPJbzN+fNMuVz2h+n2MbzS+5fBb7tzOt8uUzWva92PG8LTZ/6Pl0d01ZE/AfwPnR8zcBtySdphHK5wHAWYAD8waT/+EuK4UJmAN8F1gEvAfYDlwWLbsPOC56/mHgc7H3DWtZKUxAA/BJoAk4E9gLnFrJn5NYPj4FvAIsUXkkfizKsgyBOuCbwKnRdHC55LfYvz9pnfrJZ85jW4r5HInfzjROA+RzDnBd7HjOLiSfiWc2bRMwA9gNTIhej40OwLyk0zZC+a2KVxr95X+4y5LO4xDK4gRgbOz1F4GfAa8F1sXmzwQ6o8p1WMuSzusQymQ80ahO0evHgDdU8uckSvfrgQ8Aq4AllV4eCR+Lsi1D4NJoqi3H/FKk35+k8zHUfOY7tqV6fCnyb2fS+RlqPqPnXwLOiS+P5g8rn+rW0tcSYLO77wBw9y7gReCkJBM1Uty9O2vWEvLnf7jLSoK7/yZKd8Za4GXgZEJeMuutjZ4uLmBZSXD3nR7VKGZWBzwN3E8Ff06icjjH3W+IzV5ChZZHCiyhDMvQzAz4a+BaYKOZ/U20aAllkt8i/v6kWnY++zm2UIL5HIHfzlTKl08zmwC8Ffjv6PUpsXWGlU8F533NBLZkzesg/JutBP3lf7jLStUxhNOOFV8mZnYCcCdQDUygssvkcuDzWfMquTySVpZl6MFrCTduuR74vpm9mTLNb6Qiviv9WU7D+AAADP1JREFUHFsoj3wW+ttZKo4BvunuO9x9PiFPPwPuNLOF0TrDyqeC874c2Jk1bxyhr20l6C//w11WcsxsPtDq7k+gMgF4HrgROAX4Vyq0TMzsjcByd9+UtagiyyMlyroM3b3N3a8APkfoBlHO+a2o70qOYwslns8i/XamXlY+AXD39e5+AfBD4OLMbIaRTwXnfa0lXPwWNxFYl0BaktBf/oe7rKSY2RjgQuDj0axeeYtOSdaRI99DWFZSokrnBuAywunVSv2cXAbcYGabzWwzMBu4jXBcK7E80qBSyvDfCBedlXN+K/W7kjm2UML5LOJvZ6rlyGe2vMdz0PlMuoN92iZgOqGz/vjodTXh4oyZSadtBPMcvyAnb/6Huyzp/A2jPC4HDoi9Xgy8FHt9IOG0VM1wlyWdxwLK5kjg7kr9nADNwKzYtBo4N6qIK6480jBVShlGn73vl1t+i/H7k3QehprPfMd2oDJIOg+DyGNRfjuTzsdQ85lj+THAvwxUBv3tQy3nWdx9A3AH4fQ9hCFxbvWeTvxlJfoXB2DQf/6Hu2wUslE0ZvZJ4HFgvJkdZGYfANqB9WZ2ZLTam4Dr3X2Puy8fzrLRy1FhzGyCmR0dm3UG8JVK/Zy4e6u7r8lMwD7Cqc2XqcDySINyLUMzm2Fm7zKzqqie/ijwqXLKb7F+f0YxycOSnc98xxZK9/NczN/OUU/8EOTKp5kdaWZvi5bXAO8ljOTCcPNpUSQvMWbWRCjYVYR/sZ9w985EEzUCzGwi8HfA14HPAl9199b+8j/cZaXAzK4Arsqa/Yy7H2Zmc4HPEPJWT/iR7IreN6xlpcDMjiVc4PIM8DDwqLv/MFpWkZ+TODNbBbzP3e9XeSSnHMvQzI4gXIS9DXgQ+LK7PxMtK/n8Fvv3J61y5ROYRp5jG72npPI5Er+daZQvn4T+5bcQRmV5GPiCu6+LvW/I+VRwLiIiIiKSEurWIiIiIiKSEgrORURERERSQsG5iIiIiEhKKDgXEREREUkJBeciIiIiIimh4FxEREREJCUUnIuIiIiIpISCcxEpGWbWbGa3m9lLZvbxpNMjIiL9U709dLoJkYiUDDN7I3AfMIlwZ7ZD3H1rsqkSEZF8VG8PnYJzESlJZvYAcHKab/csIiI9VG8Pjrq1iEjJMbPZwN2q4EVESoPq7cFTcC6jysz+0cy+lnQ6MoqRHjM71cx+Z2ZLipSssmJm1xR5ewacB3y+mNsVkd5UXxfOzOaZ2bfN7FtF2t4YM7vAzF4axLpnmtkjw81rMetu1dtDo+C8jJjZYWb2P2b272b2hJm5mf1d0unK8nPglqQTEVOM9KwGXl2EtKSSmZ1hZpvNbL2Zvc7MZkefMzezr5pZY7TesWb2pJlda2bjo3kHABuKnKR3Ad9x965o+yIlR/X1sJRifb0pehxXpO1VAe3AnEGsuxw4LN9CMzvCzH4QffauMLPpsWXFrrtVbw+BgvPy8hPgn939Q8DRwHczC8zsVWZ2UhKJiu/b3Z909weTSEcuQ0lPvjJ092eBzcVMVzGPV6HbcvefA18Hutz9d+6+GrgA2AM84u7t0XqPEn4MPubuO6O3vxW4o6AMxJjZp4ErgXvM7Bng0GJtW2SUqb4eolKsr919B7CuWPty973AY4NcdyMhkM+3/ClgWfTyGnePB+NFq7tVbw+dgvMyYWaTgUOA7QDu3g18BnAzayBU/JZAuhLbdzENIh/do7ivJLb1A2CWmS0CcPctwK+At8T21Qisiz57GQvc/ZkC972fu1/l7ge7+2J3X+DuDxRr2yKjRfX1yEphfV3skTeGkv6B1t0XPWb3Ay9a3a16e+gUnJeJaFiiZ4FfmNlro3kvAY8S+nkdCCw1s/cDmNlnzezT0Smt/f3KzGyxmV1nZpea2Q4ze9bM/sHM3mFmvzezs8zsx2a2zcwujL0v5/ay9n2Bmd0Y73sXndq93sw+Z2b3mdmx0fx+9xd7v5nZhWb2lJmdZmZ/MrONZva+Qezj4Hh6BthnnzLMoSk6Td1hZhdklemVUdn8T1SZD7bM+uzLzM4xs0fN7H1mtsLMPt/P9nId+5zp6U/UwvIUcG5sdjNwpkVdWIC3AbfH0jke2BF7vcDMPm9m15jZry10jznfzJab2dvN7BYzazezS8xsUZTHLWZ2xkDpEyklqq8rp77OJV8e+9uPmR1oZjdYaIX+18HsJ8d+f25m3zez4wdYT3V30txdU5lMwBHAS4R/yjcA02LLVgFLoudHAi9Ez5sJ/+qnEP75vwwsjpZ9Dfht9LwG2Ap8lTBW6d8Aa/vbXva+CX3uvgXcGM2vJQR8jdHrc4BWYGJ/+8uR7+nRPj8ZrXs1sAuYO8A+stPT7z7jZZgjDauAa4EJ0ftejuZPAv4rtt5DhIp1UGWWZ1+TCacqvwwcC5za3/ayjn3O9Azy83UFsDJ6/r+AzwG7gbNjnxeLrf824ITo+QTg98CE6PUDUTlMAl4BPguMB86PXp9HaDz4KNFnUJOmcppQfV0R9XW0/MpBlmN/x/r3wGGx9/ggP2eZ4/k64J9yLF8S7WdsbJ7q7oQntZyXEQ+tmwsJldV7gKfMLNeFLyuAd5tZNXBCNK+eUBnMBjL9hZ8mfOlw9z1AB/Ajd99G+LJOH2B72enbTe8LTM4EdnrUZxn4H8KX+uwB9pe93cw2b47WvZrQH/q0AfbRKz1D2WceP/XQv/D3QEssj1Oilq1LCWUKgyyzPPndSqgEb3P3R939niFsL196BuMHwEEWura8E/g2cBdwnpnVAx0e1d6R44HfRs/fRvgBzLTGnAVcFZVzO3Cvh37qy4FJ7v7fHk71P87QjoFISVB9XRn1dQ5589jPfk4FGtw9k57lQ9znMcDfuvs/D3J91d0JG5t0AqQ4zGwMMNvdVwEfMbPvEK5sXwYcF1/X3Xeb2UzgcuCmzCbcfZOZLQdOIfyzP5AQkO1/a+x5F1G3qHzby5PU+DYOBqpj6eo2s1XArP72NxB332VmLwCNwNQh7GPY+8zSRc93aw7wort/OXulwZSZmb2bEARnHOrh9LfH0zqEY5A3PQNx9xVm9jiha0uDu79gZj8gtGadC/wslu7MZyPTn3EuoaUrs60t8U3Hnu+jt27U/U7KjOrrXvmrhPo6Lm859nNsXkX4IzJcJwMHm1m9u3fmWac7yoPq7hRQwZWPamL9gd39CcJpwwXZK5rZycDfu/vVHkbeiLsIOMbCkF5/Ar4w0I4H2F5/Xia0xMaHmDJCX8xC1QDPjfA+BmM9vftlY2bHDaHMbgeOik05r/ofwvZypmfw2eG/gPcDL8TSNwa4hHAKOOM4Qv/ZjLXAX5lZXWy//fZ7FCljqq97q6T6Om8e+9nPtug9NQzPNYSzDl/Js/wV77mQX3V3Cig4Ly//x8xmxF7PBH4RPd9DuABmAWHYrgYzG2dmp0bLp5nZVELfx+8QTkn9lug0aaSKrNYCM7MBtpe9b4tt4ydAJ/DX0bYao/1lWmDz7S+flmidmdF7fzGIfcTTM9A+4/nIlqsVxQitYXXAHWb2RgsXHg10DPbvC2hx92di097YLsbEnve3vXi6c6bHzGrMbFnUZaU//004nf4jAHfvAO4EHvLeo7ScTujykpEp71ssjJV+MZAZ6zb7GGQf55IeOUIkD9XXVEx9PdhyzLkfQl1aQ/gDB+EsCWY2YxB1txFasd8NnGtm50Xv/aiZXUE4a/Cr2Pqqu9NgJDqyaxr9iXCRiQNtwPeAGwmnxZqi5f8ErAHeTjhNtYLQUnEmocXlRkJFdzvh4pQ90fZ2EL6sbwb2Ei46mg18Ilr+gf62l7XvS4DfEPrxHREtW0xocb2W0PfyyGh+3v3lyb9H77+E0M1iYWxZvn3Mi6dnoH3GyzBr328nXBj51agsMu97f7T8JODJ6NhcG5XzYMvs7Tnyek60v1uAmdG8/o5pr23lSc8sws0ybh7EZ+3fsl6fD7wha941Od53IvAXYAtwRTTvDMJn7XpgBtFwcsCHCT/e347y+pakv2OaNBVrQvV1JdXXc4AHCa3/Rw2Qx/6O9ZmEM5YPEy7Ofw74+yjvOetuwtmZvYQW88nAvYQ/Bh8CLiSM934toUtN5j2qu1MwWVTIIpjZNODD7v7Z6HUVoTXnQne/ItHEDcDMHDjQQx9OGSYz+5i7D2uYLhEZPaqvJU51d3lRtxaJu5jQxWFybN7JxC70S6PYaTSdQhsmMxtjZn8L3Jp0WkRkUFRfi+ruMqXgXOK+Seib9rSZrQZ+Caxx998lm6z8ogtkLolevs/CkH4ydPXAT9z9hQHXFJE0UH0toLq7LKlbi4iIiIhISqjlXEREREQkJRSci4iIiIikhIJzEREREZGUUHAuIiIiIpISCs5FRERERFJCwbmIiIiISEooOBcRERERSQkF5yIiIiIiKaHgXEREREQkJf4/A81/6xsXhroAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(12,8))\n", "plt.rc('font',family='Times New Roman')\n", "params = {'mathtext.default': 'regular' } \n", "plt.rcParams.update(params)\n", "\n", "plt.subplot(2, 2, 1)\n", "plt.plot(vehicle1.v_kmsc, vehicle1.h_kmc, 'k-', linewidth=2.0, label='Steep')\n", "plt.plot(vehicle2.v_kmsc, vehicle2.h_kmc, 'k--', linewidth=2.0, label='Shallow')\n", "plt.xlabel('Speed, km/s',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "plt.legend(loc='upper left', fontsize=12, frameon=False)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 2)\n", "plt.plot(vehicle1.acc_net_g, vehicle1.h_kmc, 'b-', linewidth=2.0, label='Steep')\n", "plt.plot(vehicle2.acc_net_g, vehicle2.h_kmc, 'b--', linewidth=2.0, label='Shallow')\n", "plt.xlabel('Deceleration, Earth g',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "plt.legend(loc='upper right', fontsize=12, frameon=False)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 3)\n", "plt.plot(vehicle1.q_stag_total, vehicle1.h_kmc,'r-', linewidth=2.0, label='Steep')\n", "plt.plot(vehicle2.q_stag_total, vehicle2.h_kmc,'r--', linewidth=2.0, label='Shallow')\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", "plt.legend(loc='upper right', fontsize=12, frameon=False)\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "\n", "plt.subplot(2, 2, 4)\n", "plt.plot(vehicle1.heatload/1.0E3, vehicle1.h_kmc, 'm-', linewidth=2.0, label='Steep')\n", "plt.plot(vehicle2.heatload/1.0E3, vehicle2.h_kmc, 'm--', linewidth=2.0, label='Shallow')\n", "plt.xlabel('Stagnation point heat-load, '+r'$kJ/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "plt.legend(loc='upper right', fontsize=12, frameon=False)\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/hera-saturn-probe.png',bbox_inches='tight')\n", "plt.savefig('../plots/hera-saturn-probe.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/hera-saturn-probe.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 }