{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 39 - Oceanus - 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 Oceanus mission study\n", "https://doi.org/10.1016/j.asr.2017.02.012" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"SATURN\") \n", "planet.h_skip = 1200.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", "vehicle=Vehicle('Oceanus-Saturn-Probe', 140.0, 127.0, 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", "vehicle.setInitialState(1200.0,0.0,0.0,28.46,0.0,-11.40,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 (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": "iVBORw0KGgoAAAANSUhEUgAAAt0AAAHwCAYAAAB67dOHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZwcZbX/8c+ZLJNlkpCQyTKBkEAg5JLARYZFQHYURdlEVNALF2VwQa/3J/cHAhevIlwFxQUvy4gKCF4R4SeCLILsiJEJEQiSQCD7ZJJJyL5Nkjm/P57qTKczPTM90zXVNfN9v171qu6q6qpT3cnpM08/9ZS5OyIiIiIiEp+ypAMQEREREenpVHSLiIiIiMRMRbeIiIiISMxUdIuIiIiIxExFt4iIiIhIzFR0i4iIiIjETEW3iIiIiEjMur3oNrOTzWy6mU3IWnaWmc0zs5Vm9mMz65vzmrvNzKPpwazlI83sNjP7npndaGb9uu9MRER6PuVsEZHi6Nv+JsVjZqOBCuCwrGXjgTOAs4EDgFuARcD3s9YvB06OXjIna5f3AZe7+3Qz+wLwbeAbMZ+GiEivoJwtIlI81t13pDSzMmA7MNHd55vZB4CX3H1btP564AB3PzV6/kPgReD3mW2i5UcAD7h7VfR8HCG5j3b3DdnHHDlypE+YMKHDMTY2NlJZWdmFs0yeziF5aY8f0n8OaY8fYMaMGevdfUhSx09Dzs5I2+eteOOVtnghfTEr3l3NmDFjhbu3epBubekGcPdmM8t+/nzOJkuAIQBmNgj4GPBvQIOZfdbd/xxtdwIwL2s/S6L9VgPPZu9wwoQJ1NXVdTjG6urqgrYvRTqH5KU9fkj/OaQ9fthR9CYmDTk7I22ft+KNV9rihfTFrHh3ZWYL8q0rxQspDwVuBXD3je4+CRgH/BF41MymRduNA1bmvHYdUJW7w8bGRqqrq3dMtbW18UUvItJFtbW1O/IVCTSOFEg5W0R6tZycPTLfdiWVzM1sEtDo7q9mL3f3pcBFUSvKJcDFgAObcnZRDmzN3W9lZWWq/hITkd6tpqaGmpoaAMxsWzubJ0Y5W0Rkl5y9It92JdPSbWZ9CIn58jY2+wkwPnq8BBie9XoDBgP1XY0l88almc4heWmPH9J/DmmPP9KYdACtKaWcnZG2z1vxxitt8UL6Yla8hen2CykBzMyJLsrJWnYZcIe7L2vjdYcCZ7j7lWZWDdzv7ntF6yYCrwG7u3tT9uuqq6tdrSYikkZmNsPdqxOOQTlbRKQD2srZSYzTnbkix7KWXQHMBAaa2d5mdqGZTTKz/c3stGib/sD5wPUA7l4HLDWzqdFuTgFuyk3eIiLSecrZIiLF0d3jdA8BPhM9vcDMfkr4efKanE1nu/svzOxE4HYzmwf8FbjO3ddkbfdJ4Gozm08YS/bKWE9ARKRA7s7WrVvp168f2aOApEFacrY7NDeHx336FGOPIiLFl0j3ku6mnypFJCnl5eU0NTWxefNmysvLC359KXQv6W6F5uwXX4Sjj4YjjwyPRUSSUlLdS0REepO+fcMPilu37jJIhxRZL2hDEpEUU9EtIhKj/v37A9DUpK7Lccn02lHRLSKlTEW3iEiM+vXrB6ilO04p6yovIr2Uim4RkRhlWrpVdMdPLd0iUspUdIuIxCjT0q3uJfFR9xIRSQMV3SIiMVJLd/zUvURE0kBFt4hIjNTS3X3U0i0ipUxFt4hIjHQhZfzUvURE0kBFt4hIjDRkYPxUdItIGqjoFhGJkVq646c+3SKSBiq6RURipAspu49aukWklKnoFhGJkS6kjJ+6l4hIGqjoFhGJkVq646fuJSKSBiq6RURipJbu7qOWbhEpZSq6RURipAsp46fuJSKSBiq6RURipCED46eiW0TSQEW3iEiM1NIdP/XpFpE0UNEtIhIjtXR3H7V0i0gpU9EtIhIjtXTHT91LRCQNVHSLiMRIQwbGT91LRCQNVHSLiMRIQwZ2H7V0i0gpU9EtIhIjtXTHT91LRCQNVHSLiMRILd3xU9EtImmgoltEJEa6kDJ+6tMtImnQt7sPaGYnA98BPunu86NlA4HrgXXAcOAyd1/blXUiIqUg7UMGpilnq6VbREpZt7Z0m9looAI4LGfVzcAL7n4F8CBwWxHWiYgkLs0t3WnJ2epeIiJp0K1Ft7svIyTaHcysCjgXeCha9CRwmplN6Oy6WE9CRKQAab6QMi05W91LRCQNur17ibs3284Z8jhghbtvjNZvM7N5wLHA1k6um99NpyMi0qa0X0iZppytlm4RKWXdXnS3YhywMmfZOqAK2NbJdTtpbGykurp6x/Oamhpqamq6FrWISAd0pqW7traW2trazNORxY+qS0ouZ6t7iYgkqaM5uxSKbgc25SwrJ7SK0IV1O1RWVlJXV9fFMEVECteZlu7sItPMVsQSWOeVXM5W0S0iSepozi6FonsJ4Sr2bEOAekJy78w6EZGSkOYLKfMouZytPt0ikgalME7300BVNJQUZtYPGA8824V1IiIlIe1DBraiZHO2WrpFpJR1e9FtLVfkGIC7NwAPAydGy08C7nf3JZ1d1w2nISLSIWlv6U5Dzlb3EhFJg27tXmJmQ4DPRE8vMLOfunsj8EXgejM7CBgDZF8x09l1IiKJS/OQgWnJ2epeIiJp0K1Ft7uvA26Jpuzlq4CL8rymU+tEREpBmocMTFvOVku3iJSyUujTLSLSY6W5pTst1L1ERNJARbeISIzS3NKdFiq6RSQNVHSLiMRILd3xU59uEUkDFd0iIjFSS3f3UUu3iJQyFd0iIjFK+5CBaaDuJSKSBiq6RURi1ANvjlNyorcY/V0jIqVMRbeISIzU0h2/AQPCfNOmZOMQEWmLim4RkRjpQsr4ZYruzZuTjUNEpC0qukVEYqQLKeNXXh7mW7aoX7eIlC4V3SIiMVJLd/zKylr6dW/ZkmwsIiL5qOgWEYlR3759Adi2bRuuZtjYqIuJiJQ6Fd0iIjEysx2Ft1q746OiW0RKnYpuEZGYadjA+KnoFpFSp6JbRCRmGjYwfoMGhfmGDcnGISKSj4puEZGY6WLK+I0YEearViUbh4hIPiq6RURipmED45cput97L9k4RETyUdEtIhIztXTHT0W3iJQ6Fd0iIjFTS3f8VHSLSKlT0S0iEjO1dMcvU3SvXJlsHCIi+ajoFhGJmVq641dVFeaLFycbh4hIPiq6RURipiED47fXXmE+f36iYYiI5FWUotvMPlKM/YiI9GRmlnQIQM/M2RMmhLmKbhEpVX0L2djMTgGuAsbQUrCXAaOBgcUNTUSkZ2hubga6v+juTTl7zz3DfPFi2LoVoh8XRERKRkFFN/Ar4L+BmUBztMyAc7saiJlNBw7LWXyquz9iZt8BroyWvebuB0WvGQhcD6wDhgOXufvarsYiIlJM7g5AWVm39+jrNTm7vBwmToR582D2bJg2rRh7FREpnkKL7nfd/cbchWY2uytBmNkU4CXgUiDT6fF24CkzGwyMBE6Oli/IeunNwGPufm/UonMb8OmuxCIiUmxJtXTTy3L2wQeHonvmTBXdIlJ6Cm12+aKZnd3K8tO6GMcGd/+auz/v7n8FlgJvuPtm4CJgNvCCuz/p7m8DmFkVobXmoWgfTwKnmdmELsYiIlJUCbZ096qcffDBYT5zZrH2KCJSPIV+AxwF3G5m27OmZuCWrgTh7gtzFp0J/N5Cs9CngB8Ay8wsu0XkOGCFu2+M9rENmAcc25VYRESKLcGW7l6Vsw8/PMyfe65YexQRKZ5Cu5d8G6gB6oDt0TIDPlvMoICPAGd7aB46wsyGA18H7jGz1e7+KDAOyL0NwjqgKndnjY2NVFdX73heU1NDTU1NkUMWEWldoS3dtbW11NbWZp6O7MKhe1XOPvro0Lf7lVegsREqKzt/QiIiHdXRnF1o0f0o8KC7b8leaGY/K3A/eZnZSGBz9sU17r4KuCpqRflaFIcDm3JeXk5L/8IdKisrqaurK1aIIiIFKbSlO7vINLMVXTh0r8rZAwfCMcfAE0/AY4/BZ4v9p4WISCs6mrML7V7yG+AbZnZM1nQs8I3Oh7qL04A/5Fn3E2B89HgJ4er3bEOA+iLGIiLSZZk7Ufbr/nHsel3OPuOMML/nnmLuVUSk6wpt6f5PYF/gAkKrBYSfKscA/1akmE4j/Bzammbglejx08DPzGygu28ys36E5P5skeIQESmKDRs2AFBRUdHdh+51OfuTn4SvfS20di9ZAuPGFXPvIiKdV2hL9w3Anu4+wd0nRtME4KxiBBMNNTXA3ZdHz6vM7FwzK4t+pvw60div7t4APAycGL38JOB+d19SjFhERIpl/fr1AAwePLi7D93rcvbuu8Ppp0NzM/z0p8Xcs4hI1xRadC9y93WtLB9TjGCADwOPZT0fDnwXeI1wtf0d7j4/a/0XgdPN7ErChTy6OlJESoq772jpTqDo7pU5+9JLw/zmm2H16jiOICJSuEK7l3yJcEOEHczsE8CPgF90NRh3/13O8zdo6Q/Y2varCGPCioiUpM2bN9Pc3Ex5eXkSfbp7Zc4+/HA44QR46im47jq4/vq4jygi0r5CW7pHmNl/ApjZIDO7nXAXspeLHpmISA+wbl1oaE6gPzf04pz9ve+BGfzwh/DGG0lHIyJSYNHt7qcCT5jZTcDfgf2Ag4DW7ngmItLrNTQ0ADB69OhuP3ZvztnV1VBTA9u2wbnnwqbcwQpFRLpZm91LzOyYPKsWAG8D3wf2As4A/r24oYmIpN/SpUsBGDt2bOzHUs7e2fXXhy4mr70GX/gC3HFHaP0WEUlCe326rwemAKsJQz/l+gVh+Kmx9IIELiJSqPr6MAx1VdUuN16Mg3J2lqFD4b774P3vh7vuCneovOEGFd4ikoz2iu5vA3WZ4aDyMbMe/1OliEhnLF68GOi2ols5O8dBB8EDD8Bpp8EPfgAbN8JNN0GfPklHJiK9TZt9ut39kfaSd7Td79rbRkSkN/rHP/4BwOTJk2M/lnJ26045Be69F8rL4ZZbwjje772XdFQi0tsUOnqJiIgU4I1o6IwDDjgg4Uh6tzPPhCefhOHD4Y9/hH/+Z3jxxaSjEpHeREW3iEhMtmzZwpw5cwCYMmVKwtHI0UfDjBlw2GGwaBF84APwla/AmjVJRyYivYGKbhGRmEyfPp2mpiamTp3KkCFDkg5HgIkT4fnn4fLLoaws3Cp+yhT45S/D8IIiInHpctFtZgeb2d1m9nUzm9D1kEREeoZnnnkGgOOOOy7ROLIpZ0P//vDf/w2vvBJGNlm6FC68EA48EO6/H5pbG/dFRKSLitHS/QXgJHf/AXCkmQ0swj5FRFLv4YcfBuCEE05IOJKdKGdHDjwQXngBfvWr0AL+5ptw9tkwdSrcfjts3px0hCLSkxSj6L4EOCR6/BugXxH2KSKSam+//TYvv/wyFRUVnHLKKUmHk005O0tZGXzmMzB7NvzP/8Aee4Ti+6KLYK+94MorYe7cpKMUkZ6g4KLbzD5sZudFjw8D9nX3JQDu3uzua4sco4hI6txxxx0AnHnmmQwcmFxjsnJ2x/TvD1/6Erz7LtxzDxx8MCxfDtddB/vuC8ceG+5ouXp10pGKSFoVVHSb2bXAfcCZAO7+N+BfzOzcGGITEUmldevWcfPNNwNw8cUXJxaHcnbh+vWDc88No5w8/zxccAEMGgTPPQf/+q8wahR8+MOh+8mKFUlHKyJpUmhL95HAHsDzWcseBb5btIhERFLuxz/+MatXr+aoo47iqKOOSjIU5exOMgtDDP7yl9DQEIrs44+H7dvhscdC95MxY8I211wDf/tbWCcikk+hRfcT7r4a8KxlveZ2wiIi7Vm4cCHXXXcdANdcc03C0ShnF8OQIfC5z8FTT4UC/Gc/C3e5LCsLN9i5+mo4/PDQCv6pT8Ftt8E//gHu7e9bRHqPvgVu/66ZXQGMM7NPAZ8g/Gz5laJHJiKSMtu3b+f8889n06ZNnHPOORx//PFJh6ScXWSVlfD5z4dp7dpQiD/+eJjmzQu3m7/33rDt7ruHlvAPfCBM//zPoe+4iPRO5gX+KW5mRwAXAOOBpcD/uvuTxQ+teKqrq72uri7pMESkh7vqqqu49tprGT16NH//+98ZM2ZMl/dpZjPcvboLr1fO7gbuYZSTP/0p9P9+/vkw/ne28nI46CA49NAwVVfD/vtDnz7JxCwixddWzi60pRt3/yvw16yd9zGzE939z12IUUQk1W699VauvfZaysrKuPvuu4tScBeDcnb3MAujnOy7L3z5y6EInzcvFN/PPx/GA58zJ/T9/tvfWl5XUQHve18YLeXAA2HaNDjggHDxpoj0LG0W3Wb2B6CinX2MBFYASuAi0ivdeuutfOlLXwLglltu4aSTTkokDuXs0mEGe+8dpvPPD8tWrw6jotTVwcsvh2nhwtAy/txzO7920qSWIvzAA0MhPnFiGF1FRNKpvZbu16Jt5gDNwEXATCD7d7+JwIhYohMRKWHbt2/n6quv3nHh5He/+11qamqSDEk5u4TtthuceGKYMpYvD0X4q6/C66/Da6+FFvG33w7T/fe3bNu3L+yzT+iSMnlyy7T//qH/uIiUtvaK7hsI/b5XA5jZPu5+de5GZvZgHMGJiJSqhoYGzjvvPJ566inKysq47bbb+PznP590WMrZKTNqFHzkI2HK2LIl3CEzU4S//nq4S+aCBaEgnzNn1/3svjvst18oyjMt7BMnhnlVVRhpRUSS1WbR7e5rchaNNrM+7r5jNFIzOw44LIbYRERKjrtz5513cumll7Jy5UpGjRrFPffck1iXkpzYlLN7gMwFlwcdtPPyjRtD6/ecOaEozxTgc+bAypXw0ktham1/Eya0FOOZgnzPPcM0alTo0iIi8Sr0QsrfAa+a2QvAOmAK8CFgl5YUEZGepq6ujksvvZRnn30WgJNOOom77rqLsWPHJhxZXsrZPcigQa0X4+5QXx8K8nnzwq3ss6fly/O3kEMoyvfYA8aPbynEM48z86FD4z8/kZ6uoKLb3Z8wsw8CXwIOBFYDZ7t7UX6qNLPvAFdGT19z94PMbCBwPeELYzhwmbuvjbbPu05EpFhmzZrFN7/5TR544AEARo4cyY033shnPvMZrISbCJWzewczGDcuTMcdt+v69et3LsbnzQvTokVheu89eOedMOUzdGgovseNg7FjQ5eV3GnMGI1DLtKWzgwZWA9clXluZmVmdoi7z+hKIGY2mHBV/cnRogXR/GbgMXe/18xOAW4DPt2BdSIinebuPPnkk9x444089thjAAwYMICvfvWrXHbZZYwYkY5rEZWzpaIijIIybVrr6zdsCMX3woX552vXwhtvhKktlZUtRXhrxfno0aE7y4ABxT9PkVJX0M1xzOwXrSyuBLa7+xldCsTsa9HDW919c7SsCpgHDHf3jWbWF1gDHAA05Vvn7vOz953GGy2ISDIaGhq4++67+cUvfsGbb74JwMCBA7nwwgu54oorqKqq6tZ4unJzHOVsKQb30Gd80aJww5/6+pYp+3lDAzQ3d2yfQ4eGAjzfNGpUy+OK9gbBFCkhxbw5zr6EsV2zK/U9gOc7GRsAFn6f/RRwKPAtM/uCu/8vcBywwt03Arj7NjObBxwLbG1j3fyuxCMivcvatWt55JFHuOeee3j00UfZvj1cdzh27FguueQSLr74YnZP55hsytnSZWYwcmSYDj44/3bbt4f+47mFeXZxvmxZ2Gbt2jC9/Xb7xx80KH9hPmpUS2yVlWEUF3VxkVJVaNH9eXff6VIMM9ufcIvhTvPQ3H6EmQ0Hvg7cY2argXHAypzN1wFVwLY21u2ksbGR6uqWPzpqamqSHktXRBK2YsUKHnroIR544AH+9Kc/0dTUBEDfvn05/fTTueCCC/jIRz5C/wS+wWtra6mtrc08HdmFXSlnS7fp0yd0KRk7NtxlMx93WLUqFOCZIjzzuLVp48aWfugdMXRoKMCzi/G2Hg8bpiEVpWs6mrML6l7S6g7MJgIvu3tXvhhy93ktUA08AXzC3Q/PWvcK8Ovoaavr3P372fvTT5UisnXrVqZPn87jjz/On/70J15++WUy+c/MOOqoo/j4xz/Oueeey6hRoxKOtkVXupfk2Z9ytqSGO6xbl79Ab2yEFSvC1NgYusFs397+frP16dNSiGcX4yNGhGn33Xd9PHy4WtSldUXrXhL9FJhdpfcBRgO/6nx4rfoJ8AywhHCFe7YhQH0UR751ItLLNTU1MXPmTF544QWee+45nnnmGdaubRkoo3///hx77LF8/OMf5/TTT2fMmDEJRhsP5WxJO7PQcj10KOy7b/vbNzfDmjW7FuPZ89zHa9e2FPGFqKjYuSDPLdBbK9hVrPduhXYveYLQYpFJ4s3AMnd/q6hRhf2+AjwN/MzMBrr7JjPrB4wHngW2t7FORHqZFStWUFdXx4svvsgLL7zA9OnT2bRp007b7L///nzoQx/igx/8IMceeyyDBw9OKNpuo5wtvUpZWShshw8Pd+jsiC1bQgt5bjH+3nthWrmy5XH28/Xrw7RgQfvHyFZR0X6BnjsNH64RX3qCQovuu939uewFZjbZzCa5+9zOBhFd8X4c8BvCl8PXgSvdvcHMHgZOBB4GTgLud/cl0evyrhORnuu9995jxowZ1NXVUVdXx4wZM1jQyjfflClTOProoznqqKM4/vjjGT9+fALRJko5W6Qd5eUtQxp2VHNz6PbSXmHeVrG+cGFhcQ4atHMRnq9Az50GD9YdR0tFoUX3KcBzOcuWEVo32rimuV3Dge8CVwAvAD/KGkLqi8D1ZnYQMAbIvpqmrXUiknJNTU289dZbvP7668yaNYtZs2bx2muvMX/+/F22HTRoEAcffDDvf//7OfrooznyyCOprKzs/qBLi3K2SAzKysIFmMOGwcSJHX9dpo96W4X6qlWtF+wbN4Zp8eLCYu3Xr/3CvLUiXheYFl+HLqQ0s0uA7xD637XmCXc/pZiBFZMuyhEpbZs3b2bu3Lm89dZbvPnmmzuK7Dlz5rBt27Zdth84cCAHH3ww1dXVHHLIIVRXVzN58mT69OmTQPTx6syFlMrZIj2Leyi4c4vxjkwbN3bumGaFtahnF/B9C771Ys/R5Qsp3f2nZvYccD7hgplsTe6+tIsxikgP19zczOLFi5kzZw5vvfXWTvMFCxbQWgOAmbHPPvswbdo0pk6dumOaPHkyfXtzVm+HcrZIz2IWuokMHgx77lnYazdvbmk9b60VPd+0Zk3L40INGVJ4sT5iRM/vt97hby13f83MLnf3rbnrzGy8uxfYO0lEepr169czb9483n333Z3mmce5FzZm9OnTh4kTJzJ58mQmT57M1KlTmTZtGlOmTOkNFzvGQjlbRCAUspnx0wuxbRusXt2xAj23mF+3LkyFXmQ6ZEjLTY8yU/aNkLKXjRiRvu4vbRbdZnYp8Ii7/yNa9H/NLPf3WwOOJ1xUIyI92Pr161m0aBGLFi1iwYIFOxXU8+bNo7Gxsc3Xjx49msmTJ7Pffvux33777Xi89957J3ITmp5GOVtEiqVv35axywuRe5FpR6eVK1uK9Xfeaf84ZWVhTPVMET56NOyxR5jGjWt5PGpUGIu9FLTX0j0R2C3r+fuAwUADLUNQGeGCGBFJsaamJpYsWcKiRYtYuHDhjuI6+/GqVava3Ed5eTkTJkxg7733ZuLEibvMhw0b1k1n02spZ4tIorpykWlmzPTly1um7OfZj7Pvavr66/n326dPGJkmU4xPnBjGfJ80KcyrqrqvxbzNotvdv5yz6Fpglrs3ZS80s38qdmAiUhzuzqpVq6ivr2fp0qXU19fvNC1evJhFixbR0NDQar/qbOXl5ey5557sueeejB8/fpfCesyYMZSl7fe+HkQ5W0TSyqylWO/IGOtNTWFc9Uwx3tAQRnZZsiTMM48bG2HRojC1ZuDAUIC/731w6KFw2GHhcRyt44VeiTTN3V/JXmBmk4F9gH+0/hIRiYO7s3bt2p0K6NaK6vr6erZs2dLu/srKyhg3btyOgjq7uM7MR44ciWnA1zRRzhaRHql//9ByPW5c29tt2QL19aEIX7QI3n0X5s6Ft98O8+XLQ0v566/DnXeG14wdC1/8Ilx5ZXFbwTtUdJvZBwg/SZ4c3VY42wjC1fEPFS8skd7J3Vm3bh3Lli2joaGBZcuWtTlt7OBYUEOHDqWqqmqnaezYsVRVVe0otKuqqjQiSA+hnC0iEpSXhy4l+bq6rF0Lb74JdXXw9NNw//2wdClcfTVMmQJnn128WDr6DbsauIvQX/ConHVNwO3FC0mkZ3F31qxZs0vBnK+o3rx5c4f3PXjwYMaNG7ejgG5tGjt2rEYA6X2Us0VEcriHoRAXL4b581umd94JLd3ZF3D26xf6fBdTR8fpft3MjgM+4e5K1tKrZbp1NDY20tjYyPLly3c8zresI907MgYPHszo0aN3mcaMGbPLsiFD8t37RHoz5WwR6Q0yd/hs7e6emYssGxp2ntr6Ou7XD6ZNg5NPhgsv7Fjf8kIUMk73GvK0jpjZ0e7+QtGiEulGmZbofAVza0V1U1NT+zvOUlFR0WbxnD1VVFTEdKbSmyhni0ipa26G9etDF481a1rm2Y8z89wxwzOF9vbthR2zoiKMWDJxIkyY0DJNnRqK7DhHr21vnO4/AO1VAAYcAIwqVlAiXeHurF69ukMF9PLly1mxYgVbt+5y/5A2DR48mFGjRlFZWbljyn2evWzQoEExna1IC+VsEYmTe2gpXr8+tDCvX7/r1NrytWt3LaIzy9oZNKtdFRUtd7TcffeWx5mxu8eMaZlGjw539UxKey3dcwh9A5fQMsZrLgM+WsygRLJlLi7MFMytzbMfd6aIHjJkSIcL6MrKSgYOHBjT2Yp0iXK2iNDcDBs3hmnDhjBlHre2rK2COXdZoS3L7Rk8uGWowKFD255niurMfPjwcKFkWrRXdH8PWOfubXZINTNdBS8d5u5s2LAhbwHdWiFdaHeOoUOHtllE5z4fMGBATGcr0q2Us0VSYPv2tovg1pa1tz57WQHX4xesvDy0LmemIUN2fp5veWuF9dChpXO3yB+LZwoAACAASURBVO7Q3s1xVnRwP2eiq+F7tU2bNrVaLOebb9q0qaD9Z3fnaGs+atQoRo4cqSJaeiXlbJHiaK8o7sjztgrjAq6t77RBg0Ircmae/Th3WXaB3FYRPXhwuNhQOqdLg/Ka2Xjgy8AlKIH3ONu3b2f58uUsXbp0x9TQ0LDT82XLltHY2Mj69esL2veAAQN2KZbbKqjVJ1qk65SzpadoboZNmwovfjv6PM6WYgh3X+xoQdyZZQMHhmNIaSm46DazPsBpwEXAycAyQh9CSYlt27ZRX1/PkiVLdiqgc4vq5cuX09zc3KF99u/fv0Ot0JnHgwcP1p0NRbqBcrYkxT0UsevWtUzr13e91XjDhlBwx8ms9aK2tef51rX1mvJyFcW9UYeLbjObCHwe+Fdgd2AB8DHgceD9sUQnBXN3VqxYwcKFC1m0aNEu80WLFlFfX9/hYrqyspKxY8cyZswYxo4du8s0evRoRo0axdChQ1VEi5QQ5WzpjC1bdi6Ss6e1a/Ova21av77rI1O0ZeDA/EVtZwvlzDRggIpiKb52i24z+wShheQEYC5wI3AncIW7PxZt9pfYIpRdbNq0iXfffZd33nmHuXPn7pjmz5/PokWL2r2joZntuP13pnhuragePXo0/dR5SyRVlLN7t82bw3jG2dOqVR1btno1FDjwU7sy/YUzU0VF5wvl3O4TZWXFjVUkbh1p6d5EuG3wX4DPu/tbAGYW49+v4u7U19cza9YsZs2axZtvvsncuXN55513WLx4cZuv3W233Rg/fjx77rlnq/Oqqir6xzn6u4gkSTm7h9i4EVasgMbGMM9+nDtftSpMXb1Ar1+/nYvkjk5Dh+66rKKid41MIdKedotud38YeDi6AOciM5sAPAjsaAI1M3OP80eknm3Lli28+uqrvPzyyzuK7FmzZrF6devdLvv27cvEiRPZZ599mDRpEpMmTWKfffZh4sSJjB8/XrcGF+nFlLNLm3u4i159ff6poSEU0hs3Fr7/fv1gt93C+MW77dYy5T5vbdmwYaFbhYjEo5DbwC8E/tPM+gJnAHuY2TXAz4EPArXxhNjz1NfX88wzzzB9+nSmT5/OzJkzWx2HesSIEUydOpWpU6fyT//0T+y7775MmjSJ8ePH07dvlwaeEZEeTjk7Ge6h9XnePHj33TBlHs+bF4rqjt52oLwcKith5Mgw5Xs8cmS4Wchuu2nUCpFSVnDl5u7bgN8BvzOzScC/E/oPKoHnsWnTJp566imeeOIJnnjiCf7xj3/sss3+++/PYYcdxsEHH7yj0B49erQuThSRLlHOjoc7LF4Ms2aF6Y03wnz27DC6RluGDYOqqjCNG9fyODONGROK6sGDVUCL9CRdai5197nAv5nZjCLF02Ns2bKFxx57jN/+9rf84Q9/2Gkc68GDB3PMMcdw5JFHcvjhh3PooYey2267JRitiPQGytmdt3o1/O1v8NJLYZo+PSxrzbBhsPfeYZo4cefH48aFYlpEep+i9FFw97u6ug8zOwv4ATAUuBv4etRCg5ndDZwXbfoHdz89Wj4SuJYw5mw/4DJ3L/K114VZvHgxt9xyCz/72c9obGzcsfyQQw7hwx/+MCeffDJHHHGELmQUkcQoZ7fPHWbOhEcegT/+MRTZub3gR46EqVPDdMABYT5lSujqISKSqyQ6BkcX/JwBnA0cANwCLAK+H61bTripA8CcrJfeB1zu7tPN7AvAt4FvdFvgWZYuXcq3vvUtbr/9drZv3w7AtGnTOPfccznnnHPYe++9kwhLRKToekLOzmfRIrjzTrjjDnjnnZbl/frBIYfA+98fpiOOgD33TCxMEUmhkii6gb2AC6NWkhlmNhU4Hvg+of/hi8AzmVYUADM7Apjs7tOjRQ8Bc8zsO+7eTo+64tm+fTs/+tGPuPrqq9m4cSNlZWWcc845fOUrX+Goo45Sn2wR6YlSm7PzeecduPZauOsuiNpNGDsWPvYx+MhH4MQTwxB4IiKdVRJFt7s/n7NoCTDEzAYR7qD2b0CDmX3W3f8cbXMCMC9rH0uiArcaeDb+qKGhoYFzzjmH558P4Z955plce+21TJkypTsOLyKSiLTm7NZs2wbXXw/f+lYYVaRPHzjnHLjwQjjpJI0zLSLFUxJFdysOBW5w943AJDMbS/gZ8lEzO8TdXwfGAStzXrcOqMrdWWNjI9XV1Tue19TUUFNT06UAZ8+ezSmnnMKCBQsYM2YMt99+O6eeemqX9ikiAlBbW0tt7Y7BRUYmGUsHlXzObs369XD22fD44+H5Zz8L3/wm7LNP0Q8lIj1YR3O2ldr9EaIhrb7s7v/eyrp7gPXufrGZ/RSodPdPZq1/D6hx999lv666utrr6uqKFuOCBQs48sgjqa+v5/DDD+cPf/gDo0aNKtr+RUQyzGyGu1e3v2Uy0pCzW7N9O3z0o/DYY2F4vl//OrRsi4h0RVs5u6y7g2mLmfUBLgYuz7PJT4Dx0eMlwPCs1xowGKiPM8ZNmzbx0Y9+lPr6eo455hieeuopFdwi0iulIWfnc8stoeAeORJefFEFt4jEr6SKbuBS4PvuviXP+mbglejxE8DkrHUTgCYg1uaRyy+/nFmzZrHffvvx4IMPMmjQoDgPJyJSyko+Z7dm27Zw0STAbbfBvvt2dwQi0huVTNFtZlcAM4GBZra3mV1oZlPN7LRofX/gfOB6AHevA5ZGV80DnALc5O4dvMFu4WbOnMlNN91E3759+fWvf60b2ohIr5WGnJ3PCy9AQ0Mots88s7uPLiK9VUlcSGlmVwHX5CyeDVwC3G5m84C/Ate5+5qsbT4JXG1m84EK4Mq4YnR3vv71r+PufPWrX+WQQw6J61AiIiUtDTm7LS+9FOYf/rBusy4i3ackim53/w7wnTyr83aYdvcFwOdiCSrHCy+8wNNPP83w4cO56qqruuOQIiIlKQ05uy1vvBHm06YlG4eI9C4l072k1N1www0AXHLJJQwfPrydrUVEpFTNie6RqVsqiEh3UtHdAXPmzOGhhx5iwIABXHLJJUmHIyIineQOb70VHu+3X7KxiEjvoqK7A26//XYAzjvvPA0PKCKSYo2NsHYtDBsWhgsUEekuKrrb0dTUxJ133gnARRddlHA0IiLSFdmt3LqIUkS6k4rudjz88MM0NjYydepUDjvssKTDERGRLnj77TBX1xIR6W4qutvx85//HIDPfe5zmJpFRERSbfr0MD/ggGTjEJHeR0V3G9avX8+TTz6JmXHeeeclHY6IiHSBOzz+eHis276LSHdT0d2GJ598kqamJo444ggqKyuTDkdERLpg7lyYPx9GjID3vS/paESkt1HR3YY///nPAJxyyikJRyIiIl0VpXROOgn69Ek2FhHpfVR0t6Gurg6AI488MuFIRESkq15/PcwPPzzZOESkd1LRnce2bdv4+9//DsAhhxyScDQiItJVs2eHue5EKSJJUNGdx4IFC9i8eTN77LGHbvsuItID1NeH+fjxycYhIr2Tiu485s6dC8CkSZMSjkRERIph9eowVzuKiCRBRXce77zzDgD77LNPwpGIiEgxrFkT5sOGJRuHiPROKrrzUNEtItJzNDfDpk3h8aBBycYiIr2Tiu48Fi5cCMCECROSDURERLps8+YwHzAAdHNhEUmCiu48li1bBsCYMWMSjkRERLpKrdwikjQV3XksX74cgNGjRycciYiIdNXGjWE+cGCycYhI76WiO49MS7eKbhGR9Mu0dKvoFpGkqOhuRVNTE6tXr6Zv374ao1tEpAdQ0S0iSVPR3Yo10bhSu+22G2VleotERNJORbeIJE0VZSvWrVsHQEVFRcKRiIhIMWzdGub9+iUbh4j0Xiq6W5EpuocMGZJwJCIiUgyZYQLdk41DRHovFd2tUNEtItKzqOgWkaT1TTqArjKzgcD1wDpgOHCZu6/tyj5VdIuIxCOOnN2x44a5im4RSUpPaOm+GXjB3a8AHgRu6+oON2zYAMDgwYO7uisREdlZ0XN2R6joFpGkpbql28yqgHOBL0eLngTuN7MJ7j6/s/vdtm0bAP10xY2ISNHElbMBmpuhqQm2bYPt23eetm2D6H5nKrpFJDGpLrqB44AV7r4RwN23mdk84Fhgfmd3un37dgANFygiUlzHEUPOfu89OPBAWLKk/W1VdItIUtJedI8DVuYsWwdUZS9obGykurp6x/Oamhpqamry7jRTdPfp06dYcYqIdFhtbS21tbWZpyOTjKXIYsnZb77ZUnAPHAh9+0KfPi1T5nm/fnDhhUU6ExGRSEdzdtqLbgc25SwrB7ZmL6isrKSurq7DO21ubgZUdItIMrKLTDNbkXA4xRRLzs448kh48cXOByci0hkdzdlp7z+xhHD1e7YhQH1XdqruJSIisYglZ4uIpEHaq8qngapoCCrMrB8wHni2KztV9xIRkVjEkrNFRNIg1UW3uzcADwMnRotOAu539w5cTpNfpnuJWrpFRIonrpwtIpIGae/TDfBF4HozOwgYA+S/2qaD1NItIhKboudsEZE0SH3R7e6rgIuKuU9dSCkiEo84craISBqo/0QrdCGliIiIiBSTqspWqHuJiIiIiBSTiu5WqHuJiIiIiBSTiu5WZIpuM0s4EhERERHpCVR0t0FFt4iIiIgUg4puEREREZGYqegWEREREYlZ6sfpjsO+++7LaaedxpQpU5IORURE2jFiBJx2Guy/f9KRiIjkp6K7FWeddRZnnXVW0mGIiEgHTJkCDz6YdBQiIm1T9xIRERERkZip6BYRERERiZmKbhERERGRmKnoFhERERGJmYruVtTW1iYdQpfpHJKX9vgh/eeQ9vgjI5MOIC3S9nkr3nilLV5IX8yKtzAquluR9IdSDDqH5KU9fkj/OaQ9/khl0gGkRdo+b8Ubr7TFC+mLWfEWRkW3iIiIiEjMzN2TjiF2ZtYILCjgJSOBFTGF0110DslLe/yQ/nNIe/wAk919SNJBdKdO5OyMtH3eijdeaYsX0hez4t3VXu7e6i+UvaLoFhERERFJkrqXiIiIiIjETEW3iIiIiEjMVHSLiIiIiMRMRXcWMxtoZjeZ2XVmdouZDU06ps4ws++YmUfTq0nH0xFmdrKZTTezCVnLUvV5tHYO0fKS/zzM7Cwzm2dmK83sx2bWN1o+0sxuM7PvmdmNZtYv6VjzyXcO0bq7sz6DB5OMsy1mdqyZzTGzVWb2o6zlqfq/0N3S+P6kJC+kKi+nKQenLeemMb+WYj7t2/4mvcrNwGPufq+ZnQLcBnw64ZgKYmaDCVfnnhwt6swIAN3KzEYDFcBhOatS83nkO4c0fB5mNh44AzgbOAC4BVgEfB+4D7jc3aeb2ReAbwPfSCrWfNo6h2jdclo+gzmJBNkOMxsGHAUcARwJ/N7MHnb3J0nR/4WEpOr9SUleSFVeTlMOTlvOTWN+LdV8qtFLImZWBcwDhrv7xuivuDXAAe4+P9HgCmBmX4se3urumxMNpgBmVgZsBya6+/w0fh655xAtK/nPw8w+ALzk7tui59cTEus1wAPuXhUtH0dIqKPdfUNS8bYm3zm4+6lm9kPgReD3mfWlyMwGAps9SspmNgO4lPCep+r/QndKaa4o+bwA6cvLacnBacu5acyvpZpP1b2kxXHACnffCBD945kHHJtkUIUwMwM+BfwAWGZmibc+dJS7N+csOo6UfR6555CWz8Pdn89JlkuAhcAJhPc8s92S6GF1N4bXIfnOwcwGAR8Dfhs9PzGRADvA3TdlfUEMBt4EniGF/xe62XGk6P1JS16A9OXltOTgtOXcNObXUs2nKrpbjANW5ixbB1QlEEuneHAE4ae0m4B7zOzDCYfVWfo8knMocCvp/gwOJbRsbXT3SYRz+SPwqJlNSza0tkWtSo8C/YBBpPtz6A6pen9SnBdA73Vc0pZzU5NfSy2fquhu4cCmnGXlwNYEYukSd1/l7lcB/w18rb3tS5Q+jwSY2SSg0d1fJaWfQc45AODuS939IkJ/yUsSC65j5gJ3ACcCN5DSz6EbpfL9SVNeyKL3usjSlnNTmF9LKp+q6G6xBBies2wIUJ9ALMXyE2B80kF0kj6PbmZmfYCLgcujRTt9BtFPtYMp4c+glXPIVdKfAez4AvsFof/hsfTM/wvFlPb3p+T/TWbRe11Eacu5acyvpZZPVXS3eBqoijrfEw3TMx54NtGouqYZeCXpIDpJn0f3uxT4vrtviZ4/AUzOWj8BaALqujmuQuSeQ65S/wyy1RG+IHri/4ViSvv7k6Z/k3qviyttOTfN+bUk8qmK7oi7NwAPE36CADgJuD/rQoaSZ2ZVZnaumZVFfyF/Hbgy6bg6IooXwCCdn0fuOaTp8zCzK4CZwEAz29vMLgRWA0vNbGq02SnATe7elFScbWntHMxsqpmdFq3vD5wPXJ9knPmY2SAzOzhr0UeAH6fx/0J3Stv7k7K8kKq8nKYcnLacm7b8Wqr5VEMGZjGz4YR/MPOBMcA33H19okEVwMwOIFwwsBZ4AfiRu89ONqr2mdkQ4DOEsTO/DfzU3RvT9Hm0dg7AKFLweZjZVYShqrLNdvcpZrYXcDXhM6gAriylYaEy8p0DoX/h/xKuTv8r8D13L4mfanOZ2WGEi5FmE2L9m7vfF61Lzf+FJKTp/UlLnk5bXk5TDk5bzk1jfi3VfKqiW0REREQkZupeIiIiIiISMxXdIiIiIiIxU9EtIiIiIhIzFd0iIiIiIjFT0S0iIiIiEjMV3SIiIiIiMVPRLdIJZnaimf3FzI5rZ7upZvZbM/vPIh77hmLtS0QkCWZ2oJndV8zcmOc4D5vZp+I8RlzM7CQze6m97xlJj75JByDSFWY2BbgWeA84FDgQ+Ky73x3zoRcBB3Vgu4XAcKBPMQ4a3bJ2azH2JSICYGaHAv8NHAP8iXCzkzFAP+AWd388hsPOBYZRpNzYhpuBWXHsOMrHZwG3Am8BL0WryoGPA4e5+/wuHKKj3zOSEiq6Je1+D5zj7q+aWRnwy+44qLu/ZWYrO7DdWjNbWsRDHws8V8T9iUgv5+4vm9lvgP3d/aOZ5WZ2KvD/zOxSd7+5yMfcWOTcCICZfQx4PVPsuvsjxT5GhrtvBe41s68Az7j7VVlxPNiZfZrZKOB4d7/X3eeY2YoihSslQN1LJLXMbASwH7ABwN2bCbfP7a7brDYXebuOOB54uoj7ExEB2OVW4+7+R+C/gB+a2ZgYjlnUXG1mexBatrtba7dpfxJYUshOzKwcuAsYmLW4mN8fkjAV3ZJa7v4eMAd4zMyOiJYtAP5mZruZ2TVm9jczO83MGszsTTN7X+b1ZvZpM7vSzJ4xs+uylleb2X+Z2b1m9oCZDYuWDzCz/zGzq83s54SfRgtiZv83iukcMzvUzH4f7e+aKMZ7zawyWr7GzL6Zs4v+7r4l2tfXzez/mNmLZvbdgt9AEZH2/ZzQzeQTAGa2j5l928xuN7M/R4UuZlYe5bIrzex5M/tkZgdm9iUz+08z+6uZXZLvQK3lZDM7O8qZF5jZ22b2XTOrMrNfRvn0eTM7IdrFucAewGVmdkZ07c2TZnZ+1jEuMrNrzewuM7vTzAab2TAz+6aZ1ZnZCWY2x8zeNbMJnXnDovfkuKglHDM7y8x+bGY/MLPHzKwiOocfm9lD0ffN0ug9/ifgE2b2f7J2OTz6LlpnZhflOWaZmX3PzGqiY7iZ/b/M5yMlwt01aUrtBBwALCC0BvwCGBUt70NIYO8R+tbtBjwBvEH4Y/P9wFXRtsMJreWnAkOB32Tt/0XghujxD4Ga6PFoYDshsbYX4x2E1qJRhL5/g6Ll5cDjwINAFeHLYivwTWAQcEj0PLN9pr965rxvjh4PAS5P+rPQpElTeifgAmBxnnXLgJ9EefUBoE+0/NfAfdHjG4GzoscXAo3R408Dn4ke7xfl6mnR8zuA/4oe58vJI4DVwI+Aw4CTgO8DV0fb/gfwUFasDkyIHg8E3gQuiJ5/CPh11rYPA7cRutqeGX1fnBWd52PAtR18756Jviu+C9xA6N99Utb6FcD46PEbhO+kcuBKYDFwOPC5KO8/k4k32n4+8INo3aeBhXli+DzwcPR4dPQ+fzDpf1eadp7Up1tSzd3fMLNpwHeALwEfM7OTPPTxbgTWuvv9AGb2X8ALwN6EL5hmM/tatKvHCcn9VGD3rOVvRq+tAL4YvRZ3X2ZmiwoIdSRwE/A5d98Y7WOLmTUA89y9PjpOA/Csh/6OMwlfBiMJF2SeSviCANgI/IuZvQXcQviDQ0QkDk4oRA8D9gK+YmYAK4FhFp5cTGgwALiT0L0CQq6ty8m1Y4HXc45xAa3kZHd/z8zWAA+6+98AzGwhsMHMxhJahitaDdp9k5ktz1p0SbTfjNuB30bLVxO+Lx6IjlEXxdlRT3vUp9vM7iU09GR8wt0XmtlRhGK7Isr/9cA77j4dmB69trV9PxR9J7zcRkzvAzZF573MzN5j524qUgJUdEtqmVkfYE8PF8x8Nery8QhQS2g5yPVGNN8NGA/8yt1/HS37UbTPywhF8I9yjnUIIVmu62S4+wFHA/sAr2Ytz+3TuH3HCvfmKAFnuoHt7qFLDe4+z8wuAH4K/DtwHrAcEZEiMrPdCL/SvUHImytbyY+jCS2x/QHcfTuhoYDoNd9y979Ez3d6bZZWc3LE2TlXLgYuA94G/gpMbOMUsl+3L/BU1vN3CF1nRrFrLt5G50dWeYXwx0nGajP7HnAfIfZMZZ17Xu3ZRv667VHgFjMbQPgsNhFazaWEqE+3pNmOfoYA7v4qcAWwf57t+xMS3FxgKXB29kozOzxafqqZDcxZvjZ6OqWTsf4F+DFwT5QUC2JmVWRdlGNm4wgjt+xHuLDy13leKiLSFecTflm7j5Afj4pamIEd+XEF0ASckrV8mpkNISfXmlnfqBEjV76c3JpbgfkehobdUsC5LCS0jO84BLCG0H2maNy9OWoYOSD6o+UR4NvuXlfM4+Qc8yHgIcIvsp8GjnX3NXEdTzpHRbek3VeigjRjHKEvXsYIM+sfPT6O0J9vNfC/wJlm9lMzOy5qhdhASI6DgYfN7EPRRSv7Ewr1mcA1ZtbfzEYS+lKPMrM+Zna5mX0mT4wWTVcTvphubGVd7va5zz9K6H+YsS9whruvBb7SymtERAqxSwuqmR0DfAu42N0bCa3Ky4BHzOxj0cWSH4hatn9DGOXkDDM7ETjf3dcRcu1XzexbZnYsofEh0zUvO//ly8kZ2a3OBwOVZjYIOAoYaGZ7R+u2Ei48zDS+ZB/jVuDjZjY8en448DN330aoh1rNxWZ2RBRXeZ73royceir6w+RUwq+bwwndFicDewKDsi7SzG1Nb2ol/p2DaqUPSvRZDAX+TGjNL+SPEekm6l4iabcn8IaZ/ZHw05sBX8havw240sxWE1qFvwLg7k+Y2VcJP1GeQbgQcRaAmZ0O/A/hS+QXwH+4u5vZOYThnN4A7iX8TPjPhJtJnEPotrLTTXnM7DDCBULjCT83/h74VhTP48ARwAQz24fQAjMGOMfMZtPS6nMBMNjda3PO/ZcWbmqxjXDhkohIwaI8ci6hMPwZ0Ei4GG8Y8GF3fwnA3Zui/PgzQq57kNCXG+DfCA0WdxF+2cuMGPJzQu77IqEb3BfdfXlUVB4B7GNmd+fLyWZ2NqEf8+fNbLa7LyFc1Pl9QsF9F3A6MA14F/gVoVX+XyyManUA8BEze8TdHzCzfYH7zexFwq+lV0Yt8p8AxkTn9wbhRkFjzeyAKM5PEfqp/z7rfetHyP0HAeOiotwJfanPivb5d+BZ4GVCl5lHCTn9YcLFmweZ2RnuntnvvYQbFa2PivQxwNlmNo/Qgk30+tx7UiwkXCh6BjAAKDOzJwkXU3bXMLrSDtNnIT2VhVvn3uHuE7rpeP/h7rpFu4hIDxN1dSl395K8OZmZ/QdhJJn50fMBhF8pvpe5FkiSp+4l0pO11nUjngOFn2Ff6Y5jiYhI94kuFH1fCRfcA4CrgAlZi0cCG1RwlxZ1L5EeKepP9xnCz4WfcPf7Yj7kO9HPniIi0rP0dfdbkg4iH3ffbOFW9D+JLrKvJ3SFuartV0p3U/cSEREREZGYdXv3EjM72cymZ125m7lF6jwzW2nhtqh9c15zt4VbmrqZPZi1fKSZ3Wbh1qc3Rhc1iIhIkShni4gUR7d2L4n6RVUQ7mqVWTaecLXt2YSrjG8hDCf0/az1y4GTo5fMydrlfYQrnKeb2ReAbwPfiPk0RER6BeVsEZHi6fbuJWZWRrjr3kR3n29mHwBeisbJxMyuBw5w91Oj5z8EXgR+n9kmWn4E8IC7V0XPxxGS+2h3zx7bk5EjR/qECRMKjrWxsZHKyspOnGXPo/eihd6LFnovWsT1XsyYMWOFuyf2Jitn9056L4tL72dxlfL72VbO7vYLKbNubZ15/nzOJksINx0hGvj+Y4TxPxvM7LPu/udouxOAeVn7WRLtt5owJuYOEyZMoK6u8BtBVVdXd+p1PZHeixZ6L1rovWgR13thZguKvtMCKGf3Tnovi0vvZ3GV8vvZVs4uxSEDDyXcNQp33+jukwh3Gfwj8KiZTYu2GweszHntOqAqZxmNjY1UV1fvmGprc+8xIiJSOmpra3fkK8LQX6VMOVtEerWO5uySGjLQzCYBje7+avZyd18KXBS1olxCuAOWA5tydlFOuAXsTiorK0v2LyIRkVw1NTXU1NQAYGYrEg4nL+VsEZGO5+ySaek2sz6ExHx5G5v9hHA7bQg/aQ7Per0RbkFbX6yYMm+g6L3Ipveihd6LFr3tvVDO7tn0XhaX3s/iSuv7mcg43WbmRBflZC27jHDL7mVtvO5Q4Ax3v9LMqoH73X2vaN1E4DVgd3dvyn5ddXW1q9VERNLIzGa4e3XCMShni4h0QFs5O4lxujNX5FjWsiuAmcBAM9vbzC40s0lmtr+ZnRZtlGhPGwAAIABJREFU0x84H7gewN3rgKVmNjXazSnATbnJW0REOk85W0SkOLp7nO4hhFtzA1xgZj8l/Dx5Tc6ms939F2Z2InC7mc0D/gpc5+5rsrb7JHC1mc0njCV7ZawnICLSiyhni4gUT6+4DXynfqo86SRYtgz+8hcYMiSewERE2lEK3Uu6W6E5e93Mdcz5/BwqDqxg/1/uH2NkIiJtaytnl9ToJSXljTegoQHWrVPRLSJSwpo3NbP+lfWUlZfM2AAiIrtQhsqnoiLM169PNg4REWmT9Q3dzX1bz//lVkTSS0V3Piq6RURSQUW3iKSBiu58MkX3hg3JxiEiIm1S0S0iaaCiOx+1dIuIpIKKbhFJAxXd+ajoFhFJBRXdIpIGKrrzUdEtIpIKKrpFJA1UdOejoltEJBVUdItIGqjozkdFt4hIKqjoFpE0UNGdz9ChYb5qVbJxiIhIm1R0i0gaqOjOZ+TIMF+5Mtk4RESkTSq6RSQNVHTnkym6V6xINg4REWnTjqJ7q4puESldKrrzUdEtIpIKaukWkTRQ0Z2Pim4RkVRQ0S0iaaCiOx8V3SIiqZBddLur8BaR0qSiO58RI8AsjF6ybVvS0YiISB5WZmDRk+ZEQxERyUtFdz59+sCYMeAO9fVJRyMiIm2w/qHqbm5S1S0ipUlFd1vGjw/zBQuSjUNERNrUZ2AfAJo3qegWkdKkorste+0V5gsXJhuHiIi0qWxQ+DrbvnF7wpGIiLRORXdbMi3dKrpFREpa2cDwdaaWbhEpVSq625Jp6Vb3EhGRktZnUNS9ZKOKbhEpTSq62zJpUpjPnp1sHCIi/7+9O4+Tq6zyP/45vW/p7oQsna1JgLAkBFEiiGwRRBEUwR+I4gYOxHGZHzrOzwXQUXFGZdBxdFzIqOMo7ggiKCooO4gEGJaEACEr2TtJdzrd6S19fn88t9LV3dVr3eqq6v6+X6/7qqp7q+49dbty6uSp53muDCrR0n1gv7qXiEhuUtE9mEWLwu3KldmNQ0REBqWWbhHJdSq6BzNnDlRXhwvk7NiR7WhERGQA6tMtIrluzItuMzvbzB41s3lJ68rN7Jtm9q9m9h0zq053W0zBwsKF4f6zz8a6axGRfJAvOfvg7CXqXiIiOWpMi24zmwFUASf22fRt4EF3vxq4Dbgxhm3xOP74cLtiRey7FhHJZfmUsw/O063uJSKSo8a06Hb37YREe5CZzQIuBW6PVt0NnG9m80a7LdagTzkl3D74YKy7FRHJdfmUswurQ9Hdtbcrjt2JiMSuaKwP6O7dZpa8ainQ4O6t0fYuM1sHnAF0jnLb+tgCPu20cPvQQ9DdDQXqBi8iE0e+5OyiyeHrrGu3im4RyU1jXnSnMBvY1WddMzAL6Brltl527tzJkiVLDj5etmwZy5YtG1509fVhQOXLL8OqVXDsscN7nYjIKC1fvpzly5cnHk7NZiwp5GTOLp5SDEDXHhXdIjK2hpuzc6HodmB/n3WlhFYR0th20LRp01gx2j7ZZnDWWfA//wO3366iW0QyLrnINLOGLIfTV07m7KIp4eusc3e/3YmIZNRwc3Yu9JXYDEzus24SsCWNbfG64IJwe+utse9aRCTP5GTOLp4ctXSre4mI5KhcKLrvAWaZWTmAmRUD9cB9aWyL1xveABUV8NhjsGZN7LsXEckjOZmzD7Z071FLt4jkpmzM050YkWMA7r4NuAM4K1r/euDX7r55tNtiD7qiAi65JNzv6bMjIjLu5UvOPjiQcpdaukUkN431PN2TgL+PHl5mZtOi+x8E3mpm1wDnAskjZka7LV4f+EC4/cEPoLU1Y4cREckV+ZSzS2eWAtC+tT2uXYqIxMrcPdsxZNySJUt81AMpE9zhpJNCF5PrroNrr40nOBGRQZjZ4+6+ZOhnjh+jydnuzgMVD9Dd1s2pe0+laFIuzBMgIhPNYDk7F/p05wczuP76cP/LX4Zt27Ibj4iIHGRmlM6NWrtfVmu3iOQeFd0jsXQpnH8+tLTAZz+b7WhERCRJ6RwV3SKSu1R0j9T110NREXzve7o0vIhIDkm0dLdtaMtyJCIi/anoHqmjjoJPfCL08X7ve2Hv3mxHJCIiQMWRFQC0rtZgdxHJPSq6R+Of/xle+UpYtw4uvzwU4CIiklUVC6Oie5WKbhHJPSq6R6OkBH7xC6ipgVtuga98JdsRiYhMeJWLKgFoWdmS5UhERPpT0T1aCxbATTeF+9dcA3/6U3bjERGZ4MoOK8NKjfaN7XTt1UVyRCS3qOhOx5vfHLqadHfDO94Bzz2X7YhERCasgqICqo6vAmDv3zTeRkRyi4rudH32s/DWt8KePfCGN8CmTdmOSERkwqp5bQ0Aex9W0S0iuUVFd7oKCuCnP4VTT4WXXw6F965d2Y5KRGRCqj65GoCmh5qyHImISG8quuNQUQG//S0sXgyrV8Ob3gSNjdmOSkRkwklu6e7u6M5yNCIiPVR0x2XyZPjDH2D+fHjsMXjd62DnzmxHJSIyoZTOLqXy2EoO7DtA0wNq7RaR3KGiO06zZsF994WZTf73f8Nl4zdvznZUIiITypTzpgCw6w519ROR3KGiO25z58L998OiRbBqFZx0EjzxRLajEhGZMA459xAAGm5vwHXxMhHJESq6M6GuLrR4n3ZaaOk+7TT4zW+yHZWIyIRQ/dpqimcU0/ZSG81/a852OCIigIruzDnkELjrLnjf+6C1FS68EK6+Grp0wQYRkUwqKCpgxrtmALDtR9uyHI2ISKCiO5NKS+G//ztcJr6gAL70pdDPW3N5i4hkVN176wDY8bMdHNh/IMvRiIio6M48M/jEJ+Cee8JAy4ceguOPh5/8BNTXUEQkI6peUUXVCVV07eli+4+3ZzscEREV3WPm9NPDjCbnngu7d8O73x2uZKnZTUREMmLuP84FYNPXNuHdauQQkexS0T2Wpk2DO+6A738famrg9tvDLCff+pb6eouIxGzaxdMonVPK/uf303BbQ7bDEZEJTkX3WDOD978fVq6EN78ZmprgIx+BV70K7r0329GJiIwbBcUFzP1EaO1ed806urt0hUoRyZ5Yim4zOzeO/Uwos2eHS8ffcgvMmwfPPBOuYnnxxfDCC9mOTkTGsYmUs2d9YBZlh5XR+lwr236omUxEJHtGVHSb2Tlm9qCZrTGztdGyHvh1ZsIb58zCVIKrVsF110F5Odx8MyxcCH/3d7B+fbYjFJE8ppwNBSUFzP+X+QCsu3odnbs6sxyRiExUI23p/jFwC3AlcHm0XBatl9EqL4drr4Xnn4crrgjrfvADOPJI+OAH4aWXshufiOQr5Wxg+iXTqX1dLZ07O1nz8TXZDkdEJqiRFt1r3f1r7n6Pu98XLfcCn003EDN71My8z3JutO2LSeueSnpNuZl908z+1cy+Y2bV6caRVXPnwn/9F6xeDe96Vxhc+d3vhuL74ovhb3/LdoQikl+Us8NxOfLGIykoK2D7/2yn4Q4NqhSRsTfSovuDZnZRivXnpxOEmR0DPAKcDpwcLSuBv5hZJTAVODtako//beBBd78auA24MZ04csYRR8BNN8Gzz8Jll0FhYeh2ctJJcMYZ8MtfQkdHtqMUkdynnB2pWFDBvOvmAbD6fatp29Q2VocWEQFGXnSfAnzPzA4kLd3Ad9KMo8XdP+ruD7j7X4GtwEp3byP8LLqakKjvdvcXAcxsFnApcHu0j7uB881sXpqx5I6FC8MVLdetCxfYqa6G+++HSy4JreJXXx22iYikppydZO4/zmXKm6bQtbuLVZesortds5mIyNgZadH9BWAZsAA4LGn5XDpBuPvGPqsuBH5jZga8A/gqsN3M3pn0nKVAg7u3RvvoAtYBZ/Td/86dO1myZMnBZfny5emEO/Zmzw6Xkt+0KczpfeyxsGNHuKz84YfDG98YrnDZ0pLtSEUkBsuXLz+YrwitxqOlnJ3ECoyjf3Q0pXNK2fvIXla/fzWuKwOLSJqGm7NtJAnHzH4KXO7u7X3W17l7bHMxmdmfgIvcfW/0eDLwceBq4Dx3v9PM/h/wHnc/Lul1jwC/dfcvJe9vyZIlvmLFirjCyz53eOSR0N/7l7+E9ujPUVUFb3tbuNrlmWeGbikiktfM7HF3XzLK1ypnp9D8ZDP/e/r/cmDfAeo/Xc9h/3pYxo4lIhPLYDl7pC3dPwc+bWanJy1nAJ9OO8qImU0F2hLJG8Dd97j7tcCXgI8mVgP7+7y8FBj/80GZwWtfCz/6UbiM/Le+Ba95DezbF9a94Q1QXw//+I/w8MPQrZ9QRSYo5ewUJr1yEgt/tRAKYeOXNrL+uvVjHYKITEBFI3z+Zwg/U15GSKAABtQBV8UU0/nAbwfY9g3g3uj+ZmByn+2TgC0xxZEfDjkEPvShsLz4YhiAedNNsHYt/Pu/h2XWrDAf+EUXwWmnqQVcZOJQzh7AIeccwtE/PJrV71vN+s+uxw848/55HqGHjIhI/Eba0v1vwFx3n+fu86NlHvC2GGMaLIF3A09E9+8BZplZOYCZFQP1wH0xxpJfFiyAz38e1qwJLdwf+1ho8d6yJbSGv+51MHMmLFsGd9wBra3ZjlhEMks5exB1767jmB8fAwWw4fMbePEfXtSl4kUkY0ZadG9y9+YU6+viCCaaaqrM3XdEj2eZ2aVmVhAN0Pk4cA1A1B/xDuCs6OWvB37t7pvjiCWvmcHJJ8PXvhauavnYY/DJT4apCHfuDHOBv+UtoZX8vPPg29+GDRuyHbWIxE85ewgzLp3Bwp8vxEqMLd/awrNvfZau5q5shiQi49RIB1L+2N3f02fdxcD33T3tixxE88nOcfevR48XAXcCe4EHga+7++qk508GrgfWE75EPu3u+/rud9wNpBwtd3jmGbjlFvjd76DvOVm0CN78ZjjnnFC0l5ZmJ04ROSjNgZTK2cPU+GAjz17wLF27uqhYWMGiXy2icmHlmMYgIvlvsJw90qL7d8Bf3f06M6sg9Ne7GFjh7mcN/ursUdE9gG3b4M47QwH+pz9Bc1KDWFkZnHoqnHVWWF71KvUFF8mCNItu5ewRaF3TyrPnP0vrc60UVBRw5HeOpO69sfwoICITRGxFd7Sz1wDvAt4IbAPeCzS5+550A80UFd3D0NEBDzwAv/893H03PP107+01NbB0aVhOOQWOPx6Ki7MRqciEkk7RHb1eOXsEuvZ18eIHX2T7TdsBmP6O6RzxzSMomVqSlXhEJL+Muug2s9MH2HQi8DrgBsJAmQvc/WPpBpopKrpHYccOuOce+Mtf4M9/hpde6r29ogJOPDEU4KecErqj1NZmJ1aRcWwkRbdydjzcnW0/2MaL//dFulu7KZ5WzIJvLWD6xdOzFpOI5Id0iu6/AscAjYREnfJpwEx3z9kOwNlO4OPChg2hAH/gAXjoIXjhhd7bzUKf8JNOgle/GpYsgcWLoUStQyLpGGHRrZwdo/0v7ef5K56n8d5GAKacO4Uj/v0IKo6syHJkIpKr0im6zyX0/dsxxAEucveb0wszc3IlgY8rO3eGaQkfeigsK1aELirJSkrgFa/oKcKXLIFjjoGikU4PLzJxjbDoVs6OmXc7W27cwtpPruVA8wGs2Jhz1RwOvfZQimqUy0Skt1j7dOejXErg41ZbGzzxRJie8LHHQhH+/PP9n1dREfqDv+IVPcvixVCpWQJEUkm3T3c+ysWc3bG9g7VXr2Xbf28Dh6JDiqj/RD2zPzybwkoNMheRQEV3DibwCaGpqacQX7Ei3K5f3/95ZmEO8eRC/PjjYc6csE1kAlPRnVv2rtjLSx97iaYHmwAonl5M/afrmfWBWRSWq/gWmehUdOdwAp9wGhrgySfhqad6lueeg64UF6OYPBmOOy70FU8sCxfCtGljH7dIlqjozj3uzp4/7WHdZ9bR/FiYarV4WjGzPjSL2R+aTcl0jWURmahUdOd4Ap/w2ttD4Z1ciD/1FOzalfr506b1FOAqxmWcU9Gdu9ydXb/bxfrPrWff4+EaP1Zq1L23jjkfm0PlMeo2JzLRZLToNrNXEi71+yThkr7r09phBuRLApck7rBlS7iC5sqVsGpVz21zqqtaE4ruRCG+cCEcdVRYZs+GgoKxjV8kJnEX3crZ8XN3Gu9r5OWvvsyuO3oaC2pOr2HWsllM/T9TKSxT1xORiSDTRfeNwFvdvc7MLgVudff9ae00ZvmWwGUQ7rBpU08RnijEByvGKypgwYKeIjx5mTRpbOMXGaEMFN3K2RnUsrqFl7/+Mttv2k53S5i1sWhKEXXvraPu/XVULa7KcoQikkmZLrqLgenuvtnMCoAqd9+b1k5jls8JXIapbzH+3HNh9pTnnw/TGw5k5sz+hfgRR8C8eZpjXHJCBopu5ewx0LW3ix0/28GW5VvY98S+g+srj61k+rumM+OdMyg7tCyLEYpIJsR9Gfg3AVPc/SdmdiKwz91XxRBnxoyHBC5p2LMnXMwnUYQnlhdfDP3JUykogLlz4fDDQxF++OG9F7WQyxiJ4TLwytlZ1vx4M1u/t5Udv9xB1+6eQePVp1Qz49IZTL1gKqWzcvZaRSIyArEV3Wb2L8BVwB/c/aJo3ZeBp939p3EEmwnjLYFLTA4cgI0b+xfjL70UWs27B7qgH6H/eHIRnijMDzsMZszQVIcSm3SKbuXs3NLd0c3uP+5mx0930HBbA937e3LMpJMmMe3CaUy9YCoVR+mKlyL5arCcPdLLab0WmAO8L2ndncCPgZxN4CIpFRbC/PlhOeec3ts6OsKc4i+9BGvWhNvEsnZt6LKycyf89a/991tWBoceGrqopFpUlMvYUc7OIQUlBUx9y1SmvmUqXfu6aPhNAzt/tZM9f9pD86PNND/azNpPraXi6AqmXjCVKW+aQvXJ1RQUayC4yHgw0qL7LndvNLPk5vGL4gxIJCeUlMCRR4alr+7uMLNKoghPLsrXrg3dWRKt5qkkivKBCvMZMzTbisRFOTtHFVUVUffuOureXceBlgPs/tNuGm5tYNcdu2hd3crGL29k45c3UjipkNoza5lyzhSmvHEK5fPLsx26iIzSSIvutWZ2NTDbzN4BXAxcCPxD7JGJ5KqCgnC1zDlz4Iwz+m/fuxc2bAgt5clLYt2uXYMX5aWlPUX53LlQXx+WxP25c6FcX7wyLMrZeaCwspBpF05j2oXT6O7spumBJnbdvovdf9xN63Ot7LptF7tuC1MRli8oZ8obp1C7tJaa02somaYB3yL5YjQDKV8DXAbUA1uBn7n73fGHFp/x2j9Q8lRzc+qiPLEMdFGgZFOn9i/Ek4vzmTND9xnJezEMpFTOzmNtG9vY/cfd7P7jbvbcvYcDTQd6ba9YVEHtGbXULq2l9vRaSmaoCBfJpkxPGVgILHX3P6e1owxSApe8kijKN24My6ZNve+//DJ0dg6+j8LCcFGgVK3kiXW1tepbngcyMGWgcnae6u7qpvnRZvb8ZQ+N9zay9+G9dLf1HvBdcXQFNafXUH1yNdWvqabiyAqsQP/ORcbKqItuM/stMNRM/lOBBnc/c/QhZpYSuIwr3d2wfXv/ojz5dvv2ofdTWdlTiCeWRLeZxFJdrcI8y0ZSdCtnTyzd7d3sfWwvTfc10XhvI00PN9Hd2rsIL5pcRPVrQgFefXI11SdWU1Qz0p6lIjJc6cxe8nT0nOeBbuBKwqWDk7PhfGBKDHGKyHAUFITuIzNnwkknpX5OW1toEU9VkCeK9X37wkWEnntu4GNVVfUvxPsuU6aoMM8dytkTSEFpAbWn1lJ7ai2HXnMo3R3dNK9opumhJvb+dS97H9lLx9YOdt+5m9137g4vMqhYWEH1idVMOmESVSdUUXVcFYUV6o4mkmlDFd3/RmgNbwQws8Pd/bN9n2Rmt2UiOBEZpbKyMHf4EUek3u4OTU29C/HNm0Ohnlg2bQqF+erVYRnsWEMV5tOmaUaWsaGcPYEVlBRQ89oaal5bA4C7076pnb2P7GXvX/fS9EgT+57YR+vKVlpXtrLtv7dFL4TKhZVUvaoqFOKvqqLq+CqKqtQiLhKnQf9FuXtTn1UzzKzQ3Q+O5DCzpcCJGYhNRDLFLPTprq2F445L/ZxEYZ5ciKdamprCtIlr1gx8vOLi0Md8sMK8rk6DP9OknC3JzIyy+jLK6suYfsl0AA60HWDfk/toXtHMvif20fx4My2rWmh5NizbfxR1TTOoOKqCquOrqDyuksrFlVQtrqK0vhTTL1siozLS/8beDDxlZg8CzcAxwBuBfi0pIpLnkgvzY48d+HnNzf1byfsuu3b1zM4ykMLC0GVmsMJ81qxQwMtwKWdLL4VlhdScXEPNyTUH1x3Yf4CWp1tofqKZ5sdDMd7yTAutq1tpXd0KP096fXXhwQK8cnHlwaW4Vv8uRYYymikDZwEfAo4D9gK/cvdYfqo0sy8C10QPn3b3V5hZOXA94QtjMvBJd98bPX/Abck0KEcky/bvH7owH87gT7PQIp6YgeXQQ3vuJ5Zx1sc8hikDlbNlxLrbu9n3zD5anm6h5ZmWg/c7d6aeOal0bimViyupOKaCiqMrqDwm3C+eomJcJpY4LwOPu28Brk3aeYGZneDuj6cRI2ZWSRhVf3a0akN0+23gD+7+CzM7B7gReOcwtolIrigvH7yPOUBHR7jS52CF+datPcujj6beT0VF/0I8eZkzJ1yAaIJQzpbRKCgtoHpJNdVLqnut79je0a8Yb13ZSvumdto3tbP797t7Pb94WnEoxPsU46Vz1U1FJp4RtXSb2Q9SrJ4GHHD3C9IKxOyj0d3vuntbtG4WsA6Y7O6tZlYENAGLgI6Btrn7+uR9q9VEZJzo6goF98aNvecyT16a+nZrTqGuLrSSH354/6WuLqdaytNp6VbOlrHgB5z9a/bTsrKF1udCl5SW50L3lO6W7pSvKagsoOKoUIiXLyinYkG4LT+iXK3jktfibOleAPwZSK7U5wAPjDI2ACz8d/cdwKuBz5vZ37v7z4ClhPlkWwHcvcvM1gFnAJ2DbFufTjwikqOKinrmFD/llNTPaWrqPTVi3wJ982bYti0sqVrLKyrgsMNCq3yiED/6aFi4EKZPz6mCfBiUsyXjrNBCAX1URa/13u20b24PhXhyMf5cK507Otn3xD72PbGv3/6KphT1K8TLF4RFfccln4206L7C3Z9PXmFmRxMuMTxqHprbX2Nmk4GPAz8xs0ZgNtD3mtjNwCyga5BtvezcuZMlS3r+07Fs2TKWLVuWTsgikqtqasIy0ODPRGv5unXw0kv9l9274dlnw9LXlCmwaFEowBPL4sUwY0asb2H58uUsX7488XBqGrtSzpassQKjbG4ZZXPLmPKG3lPDd+7uDAM1n29l/4v7w7JmP60vttK1u4vmR5tpfrS53z6Lpxb3LsQPL6dsfhnlh5VTPL1YXVYkK4abs+O4DPx84DF3T+eLoe8+/wVYAtwFXOzuJyVtewL4afQw5TZ3vyF5f/qpUkSGbc+e3kX4mjVhnvKVK2FvvzF/wZw58OpX9ywnnQSTJsUSTgYuA6+cLTnL3enY1tGvEE/c73vFzWQFFQWhAJ9fTtlhZQeL8bL54b7mHZexEFv3kuinwOQqvRCYAfx49OGl9A3gXmAzYYR7sknAliiOgbaJiIzO5MmwZElYkrmHgZ6rVvUsK1fC00/3DPS89dbw3KIiOPlkeOMb4e1vhwULxv59oJwt+cfMKJ1ZSunMUmpPr+21zd3p2NLRqxBvW9fG/rXhtmt318EL/6RSPK2YssP6FOXR/dI5pRQU6wJeklkjHUi5nNBikXhRN7Dd3V+INSizacDXCT9brgGmuft+MysG9gGHAQcG2ubum5P3p1YTEcmY7m544QV47LGwPPoorFgR1iecfTZ84xuhb/gIpTmQUjlbJoyupi72r9tP29q2XsX4/rX7aVvfhrcPUu8UQOnsUsoOLaO0PtyWHVpG6aHR/foyCit18S4ZWpwDKW9y9/v77PwoMzvC3Qe5HN2QAc4iDMD5OeHL4ePANe6+zczuAM4C7gBeD/w6kaAH2yYiMiYKCkIxffTR8J73hHWNjfDnP8Ptt8Mvfwl33RUGfq5bB9XVg+8vXsrZMmEU1RQx6fhJTDq+f9cu73Y6tnYMWJR3bOk4OO3hQIqnFqcuyKNCvfgQ9SmXwY206D4HuL/Puu3APcAr04hjMvBl4GrgQeDrSVNIfRC43sxeAdQByaNpBtsmIpIdtbVw1lnh6pnbt8Mf/hAGaO7cOdZFt3K2CGFQZ+nsUkpnl8Kp/bd3d3TT/nI7bRvaaNvQRvuGnvttG9po39ROZ0MnnQ2dKWdcgTANYll9/xbyxP2SmSUUFKkLy0Q2rKLbzD4CfBGYZGafTPGUu9IJwt1XAvUDbNsDXDnSbSIiY6K7O/TnXrUqzHiyciU88QQ880zoBw7hEvfXXBOmHxwDytkiI1NQUkD5YeWUH1aecrt3Ox3bO1IX5NHjA80HDk6PmFJhny4s9UldWaL7RZM02HM8G9Zf193/08zuB95HGDCTrMPdt8YemYhILti3Lwyg3Lo13G7YELqJrF8fbjdsCFfT7Ku4OAzGPO88uPRSmD9/zEJWzhaJlxX0DPDkNamf09nYmbog39hG+8Z2OrZ10L6xnfaNA3dhKZpcdLAgT1Wcl9SVYAXqwpKvhv1fKnd/2sw+5e6dfbeZWb27b4w3NBGRDHAPF9DZuRMaGnrf7tjRc5n5RKHd3H+u4H6mTQtzdh97bJjHe/FiOOEEKE/dajYWlLNFxlZxbTHFtcVUvaIq5fbu9m7aNoUCvG1j74K8bUMbbRvb6NrTRdeeLlqeakm5Dys2Suf0aS0/NKkwry+jsEIDPnPVoEW3mf0T8Ht3XxWt+oSZ9f1rGvA6wqAaEZGx4R4K4j17Ui+NjaEfdXJhnVi6uoZ/nLIymDkTZs0Kt/X1MG9eaLmePz9cTr4q9ZfsWFPOFsl6SzgzAAAdc0lEQVRdBaUFVBxRQcURFSm3uzudOzsHLMjbN7aH7evCQNCBHBzw2bcgj7qx6CJC2TNUS/d8IHmizFcBlcA2eqagMsKAGBGR4evsDBebaWoKt4kl+XFj48AFdWMjHDgwumNPmhRap6dO7blN3E8U14nb2tp8uvS7crZInjIzSqaXUDK9JFxqKoUDrQdo35S6IB/ugE8rtV4t4327sZTOLaWwTK3lmTBo0e3uH+6z6l+AZ929VwdGM1sYd2AikoPcobU19HNubu65Haxw7vs4cb9t4JaaYausDEXx5MkDL6mK69LS9I+dg5SzRca3wopCKo6qoOKoAVrLkwd8buxdnCced+3uOnjFz4EUzyjuX5AntZpresTRGekw2cXu/kTyCjM7CjgcWJX6JSKSNe3toTDuWySPdt2+fT0zcqSroCBMn1dTE24TS/LjmpqBi+naWigpiSeW8Us5W2QCGc6Az67mrn6t5YmCvG1DG+0vt9O5vZPO7Z00P5Z6TEtBRUGv1vK+3VhKZ5dSUKLpEfsa7pSBpxF+kjw7uqxwsimE0fG3xxybyMThHmbASBS2I1kGK5I7+42hS19ZWeieUVXVc9u3UB6oiE5+XFGRT9028opytogMpGhSEUULi6hcWJlyux9w2re29+5XnijIo/sHmg7QurqV1tUDTI9oUDKzpNcVPsuPKKd8QTnlR5RTOqt0Qs7CMtyW7kbgR4T+gqf02dYBfC/OoERyWnIXiziXkQzuG66iolAYJ4rj5EJ5oPuDba+sDPuUXKecLSKjYoVG2ZwyyuaUUXNKTcrndDV1pe5XHg0Cbd/STseWDjq2dMAj/V9fUF5A+eE9RXjy7XguyIc7T/czZrYUuNjdlawlf3R1QUtLvMVxS0t8XSySlZT0FLjDXSorBy+qx2nfZRmccraIZFJRTRFVi6uoWjzA9Iid3XRs6ThYkLeta2P/mtCPfP+a/XTu7KTl2RZanu0/NWKiIK84poKKhRVULqqkclEl5UeU532XlZHM093EAK0jZnaquz8YW1Qy8aTTvWKwJY7BeqlUVIy8QB6qeFb/ZImRcraIZEtBcUGYovDQspTbOxs7aXupjdYXW0MxnlyQ70hdkFuRUX5kOZULK6lYVEHlsZVMeuUkyg4ry5tBnUPN0/1bYKgJaA1YBEyPKyjJcfnUvaKgIN7iuKoqFNyFmk5Jco9ytojkg+LaYopPKGbSCZP6betq6qL1xVZan2ulZWULrStbaVnVQtu6NlpXtdK6qhVu7nl+UW0RVa+qYtIJk8KyJHcL8aFaup8n9A3cTM8cr30Z8OY4g5IM6uoKU7Y1NvbcJpaBHjc15W73iqGWsjIN1pOJRDlbRPJaUU0R1UuqqV5S3Wv9gZYweLNlVQstK1toebqF5iea6dzeSeNfGmn8S+PB5xZPLabmtBpqTq+h9oxaqo6rwgqzXwsMVXR/BWh29/bBnmRmGgU/ltxD8dvQALt29b5taAhX4RuoiN6XerL8EVP3CpFcpJwtIuNSYWXhwdbsBHenY0sHzY839ywrQiHecGsDDbc2hNfWFDL5rMlMfetUppw7hZKp2ak3hro4TsMw93MhGg2fnrY22Lq197J9e08hnVxU79oV+j+PhlmY37imJtwmL33XJR7X1PQenKfuFSI5STlbRCYSM6N0dimls0uZev5UIBTibWvbaLy/kcb7Gmm6v4m2dW003NJAwy0NUAA1p9VQd1kd0y+eTmHl2NUzac39ZWb1wIeBj6AEPrC2NtiwAdatC8v69bBlS+8Ce8+eke2zogIOOaTnCnvJ96dMGbigrqoK/ZxFZMJRzhaR8c7MwnSEh5cz8/KZALRtaGPXHbtouK2Bxnsbabqviab7mlhz1RrqLq+j/lP1lNZlfravERfdZlYInA9cCZwNbCf0IZzY3OHll+HZZ+GZZ2DlSli7NhTZmzcP/fqiIqirg5kze5a6up5COrmwPuSQUHSLiAxBOVtEJrqyQ8uY/eHZzP7wbLqautj5651s/d5W9j6yl83/sZmty7dS/8l66q+up6A4cw2Twy66zWw+cAVwOXAIsAF4C/BH4OSMRJfL9u6Fhx+GBx6ABx+Ep54KfaZTKSyE+nqYPz8s8+bB7Nm9C+ypU9UCLSKxUc4WEemvqKaIme+fycz3z6T5f5vZ8PkNNPymgfWfW8+u3+1i8e2LKZmRmT7fQxbdZnYxoYXkTGAN8DXgf4Cr3f0P0dMezkh0uWbnTvjVr+Dmm+G++6C7u/f2Qw6BxYvDsmgRLFgQiuy5c3UVPxEZE8rZIiLDM+n4SRx767HsuXcPqy9bTfNjzTx5xpO86pFXUTy5OPbjDacS3E+4bPDDwBXu/gKAmWVgzrgctWULXHcd/PCHPRdbKSqCE0+E006DU08N92fM0PR0IpJtytkiIiMweelkTvjbCTz1+qdoeaaFFz/8Igt/ujD24wxZdLv7HcAd0QCcK81sHnAbcPC/AGZm7pmYuDkH/PrX8Hd/19N15E1vgne+E97yljAwUUQkh0z4nC0iMgol00s49tZjeWzxY+z42Q7qP1VP1XFDXWtsZIbdidjdN7r7Zwj9AwHmmNl1UUK/MtaocsXNN8PFF4eC+7zzYNUq+P3v4T3vUcEtIjltQuZsEZE0lB9eTt3ldQBs/8n22Pc/4pF77t7l7je7+4WEfoIfA74ee2TZtmMHXHFFmJXk85+H22+HY47JdlQiIiMyYXK2iEgMpr1tGgB77h7hVM7DkNboPndfA1xlZo/HFE/uWL48tHCffTZ85jPqqy0ieW9c52wRkRhUvTJ0KWld3Yp3O1YQX/0Xyxx17v6jOPaTU26PrpJ81VUquEVkXBmXOVtEJAbFU4opmlJEd2s3nTs7Y913zkwMbWZvM7N1ZrbLzP7DzIqStt1kZh4ttyWtn2pmN5rZV8zsa2YWz/wuXV3wxBOh2D799Fh2KSIynuRUzhYRiVHxISE1de6Jt+jOicmjo1H2FwAXAYuA7wCbgBuibTsIV1IDeD7ppb8CPuXuj5rZ3wNfAD6ddkBbtoTCe9YsmDQp7d2JiIwnOZezRURiVFQbyuOuxq549xvr3kbvUOD97t4FPG5mxwKvA24gDPp5CLg32g6Amb0GOMrdH41W3Q48b2ZfdPeWtKJJXLZ99uy0diMiMk7lVs4WEYlR4aRCAA40H4h1vznRvcTdH0hOzsBmYKOZVRAuW/zL6PFZSc85E1iXtI+oUmZJ2gHt2xduq6vT3pWIyHiTczlbRCRGVhzG8nlnvJczyImiO4VXA99191Z3PwKYDfwOuNPMFkfPmQ3s6vO6ZmBW353t3LmTJUuWHFyWL18++NFbokaXysp03oOIyKgsX778YL4CpmY7nmHIbs4WEYlRQXEoj71reEX3cHN2rnQvOcjMjgB2uvtTiXXuvpVwZbUK4CPABwAnXO44WSnQr9f7tGnTWLFixfCDaG0NtxUVIwteRCQGy5YtY9myZQCYWUOWwxlUTuRsEZEYjbSle7g5O6daus2skJCcPzXAU74B1Ef3NwOTk15rQCWwJe1AOqPvgJKStHclIjJe5UzOFhGJUaLo7u7sjnW/OVV0A/8E3ODu7QNs7waeiO7fBRyVtG0e0AGk3zxyIOo4X5Brp0dEJKfkRs4WEYnRuO/TbWZXA08C5WZ2mJm938yONbPzo+0lwPuA6wHcfQWwNRo1D3AO8E1370g7mO7ofzaFhWnvSkRkPMqpnC0iEiMrjIruA/EW3TnRp9vMrgWu67N6NaEv4PfMbB3wV+Bf3b0p6TmXAJ81s/VAFXBNLAElim61dIuI9JNzOVtEJE6JC5HHW3PnRtHt7l8EvjjA5umDvG4D8HexB6TuJSIiA8q5nC0iEqMw5ITYi25Vlamoe4mIiIjIxJShlm4V3amopVtERERkYko0dPs4HUiZU9TSLSIiIjIxqaV7DCWK7kSfHhERERGZEKwgqv/inaZbRXdKiZ8T1L1EREREZGJR9xIRERERkQxT9xIRERERkczSlIEiIiIiIpmmlm4RERERkQxTn24RERERkQxTS/cYivl/NiIiIiKSH9SnOxs0T7eIiIjIxKKWbhERERGRDFOfbhERERGRDFNLt4iIiIhIZqlPt4iIiIhIpqmlewxp9hIRERGRiUl9urNAs5eIiIiITCwZKv9UdIuIiIiIZJiKbhERERGRDFPRLSIiIiKSYSq6RUREREQyTEV3Kpq9RERERERiVJTtANJlZuXA9UAzMBn4pLvvjWnnsexGRESCjOZsEZEcNh5aur8NPOjuVwO3ATdmOR4RERmYcraITEh5XXSb2SzgUuD2aNXdwPlmNi9bMYmISGrK2SKSV3RFyl6WAg3u3grg7l3AOuCMbAYlIiIpLUU5W0RynGWoe3G+9+meDezqs64ZmJW8YufOnSxZsuTg42XLlrFs2bLMRyciMgrLly9n+fLliYdTsxlLzJSzRWTcGW7Ozvei24H9fdaVAp3JK6ZNm8aKFSvGLCgRkXQkF5lm1pDlcOKknC0i485wc3a+dy/ZTBj9nmwSsCWtvWrKQBGRTMhMzhYRyQP5XnTfA8yKpqDCzIqBeuC+WPauKQNFROKU2ZwtIpLD8rrodvdtwB3AWdGq1wO/dvfN2YtKRERSUc4WkYks3/t0A3wQuN7MXgHUARptIyKSu5SzRWRCyvui2933AFdmOw4RERmacraI5A3N0y0iIiIikiEZGtKnoltEREREJMNUdIuIiIiIZJiKbhERERGRDFPRLSIiIiKSYSq6RUREREQyTEW3iIiIiEiGqehOpaoKZs6ESZOyHYmIiIiIjKHC6kJKZpZQWFkY637z/uI4GfHhD4dFRERERCaUedfOY96182Lfr1q6RUREREQyTEW3iIiIiEiGqegWEREREckwFd0iIiIiIhmmonsQy5cvz3YIOUPnoofORQ+dix46F9mnv0F8dC7jpfMZr3w9nyq6B5Gvf9RM0LnooXPRQ+eih85F9ulvEB+dy3jpfMYrX8+nim4RERERkQwzd892DBlnZjuBDaN46VSgIeZw8pXORQ+dix46Fz0ydS4OdfdpGdhvzlLOzgk6l/HS+YxXLp/PAXP2hCi6RURERESySd1LREREREQyTEW3iIiIiEiGqegWEREREckwFd0pmFm5mX3TzP7VzL5jZtXZjinTzOyLZubR8lS0bsDzMN7OkZmdbWaPmtm8pHWjev/5fm5SnYtofb/PSLR+XJ4LM3ubma0zs11m9h9mVhStn2pmN5rZV8zsa2ZWnPSaUW2T9OTz5yxXjPQ7QHqL8ztE4v0eyiUqulP7NvCgu18N3AbcmOV4MsrMKgkjgc+OlouiTYOdh3FzjsxsBlAFnNhn02jff96em4HOxSCfERiH58LM6oELCO/zY8AVwEejzb8CfuDunwReAL6Q9NLRbpP05OXnLFeM8jtAIhn4DpnQMvA9lDM0e0kfZjYLWAdMdvfWqHWrCVjk7uuzGlyGmFmimPiuu7dF6wY8D0DHQNvy9RyZWQFwAJjv7utH+/4H25Yv56bvuYjW9fuMROvH5efEzE4DHnH3rujx9YT3dB1wi7vPitbPBp4HZgCLR7PN3VvG8r2NNxMxZ8dtpN8BOq/9xfUdonMbxPU9lGvnUy3d/S0FGty9FSD60l0HnJHNoDLFzAx4B/BVYLuZvTPatJSBz8Ng2/KSu3f3WbWU0b3/wbblhb7nYpDPCIzTc+HuDyQK7shmYCNwJuE9JJ63Obq7JI1tkp6l5OnnLBeM8jtA+ojxO0SI9XsopxRlO4AcNBvY1WddMzArC7FknIefOl5jZpOBjwM/MbNGBj8PXYNsGy9G+/7H3bkZ6DPi7ncycT4nrwb+DVjGwO9psHMxofLKGNO5TcMovwNkaBMlN46JNL6HcoqK7v4c2N9nXSnQmYVYxoy77wGujf43+VHgLgY/D+P9HA31OZhw5ybFZ+RO0jtPecHMjgB2uvtTZjbY+x3tNkmPzm0MRvEdIIMb97kxG0b5PZQz1L2kv83A5D7rJgFbshBLNnwDqGfw8zARztFo3/9EODeJzwiM83NhZoXAB4BPRat6vaco8VeS4v2OYJukJ+8/ZzlmON8BMrRxnRtzwHC/h3KKiu7+7gFmmVk5QDStVz1wX1ajGjvdwBMMfh4mwjka7fufCOcm8RmB8X8u/gm4wd3bo8d3AUclbZ9HGBS1Io1tkp7x8DnLJcP5DpChjffcmG3D/R7KKSq6+3D3bcAdwFnRqtcDv04a+DSumNksM7vUzAqi1rePA9cMdh7G4zmK3juAweCfg/F+bvqei4E+IzD68zRGbyUtZnY18CRQbmaHmdn7gUZgq5kdGz3tHOCb7t7h7itGs23s3tH4lO+fs2wbzXdAtmLNZXF9h4xhyDktru+hMQ57SJoyMIWoo/71wHqgDvi0u+/LalAZYmaLCH2i9gIPAl9399XRtgHPw3g6R2Y2CXg3YZ7PLwD/6e47R/v+8/ncpDoXwHQG+IxErxl358LMriVMD5hstbsfY2aHAp8lvKcqQoGSmFpwVNskPfn6OcsFo/0OkB5xf4dMdHF/D+USFd0iIiIiIhmm7iUiIiIiIhmmoltEREREJMNUdIuIiIiIZJiKbhERERGRDFPRLSIiIiKSYSq6RUREREQyTEW3iIiIiEiGqegWkbxhZtPM7LdmtsHMPpXteEREZHDK2z10cRwRyRtm9kbgHqAaWA0c6e67sxuViIgMRHm7h4puEclLZnY/cKYupS4ikh8met5W9xIRyTtmNhf400RN3CIi+UZ5W0W3jDEz+4SZfSvbcSTEEY+Zvd7MHjGzpTGFNa6Y2Q0x78+AtwNfjnO/ItKb8nX6zGyemf2Xmd0Y0/4KzexKM9swjOeeZ2aPjva9xpm7lbcDFd3jiJkdY2a3mNn3zOwpM3Mze3e24+rj98DPsh1Ekjji2QS8IoZYcpKZnWtmDWa21cxONrO50efMzew/zaw2et6JZvasmX3VzMqjdTOAbTGHdCnwfXfvivYvkneUr0clH/P1jui2NKb9FQCNQP0wnrsCOGagjWa2yMx+EX32rjWzuqRtcedu5W1UdI83vwE+7+5XAK8EfpTYYGbHmdkZ2Qgq+dju/qy7P5iNOFIZSTwDnUN3fx5oiDOuOP9e6e7L3X8PfBvocvdH3H0TcCXQATzq7o3R8/5GSPL/z933Ry9/C3BHWm8giZl9BvgccLeZrQaOimvfImNM+XqE8jFfu3srsCWuY7l7J/D4MJ+7nVCgD7R9JbA8eniDuycX2bHlbuXtHiq6xwkzmwIcCbQAuHs38FnAzayGkNAtC3Fl7dhxGsb76B7DY2VjX78A5pjZCQDuvgv4M/DmpGPVAluiz17C0e6+Os1jH+Tu17n7Andf4u5Hu/v9ce1bZKwoX2dWDubruGesGEn8Qz33QHTbt591bLlbebuHiu5xIpp+53ngD2b2mmjdBuBvhH5U84FlZnY5gJl9wcw+E/20dLDflpktMbOvmdlHzazVzJ43s380swvM7DEzu9DMbjWzvWb2gaTXpdxfn2NfaWY/TO7bFv3E+k0z+5KZ3WNmJ0brBz1e0uvNzD5gZivN7Gwze9rMtpvZZcM4xoLkeIY4Zr9zmMLk6OfiZjO7ss85/Vx0bm6JkvRwz1m/Y5nZRWb2NzO7zMxeNLMvD7K/VH/7lPEMJmoRWQlcnLR6GnCeRV1JgPOB3ybFWQ60Jj0+2sy+bGY3mNl9FrqpXGJmK8zsrWb2MzNrNLOrzOyE6D3uMrNzh4pPJJ8oX0+cfJ3KQO9xsOOY2Xwz+4GFVuN/G85xUhz392b2EzM7ZYjnKXdnirtrGScLsAjYQPif7Q+A6Unb1gNLo/vHAmuj+9MI/ws/hPA/9Y3Akmjbt4CHo/slwG7gPwlzbb4T2DzY/voem9Cn7Ubgh9H6MkIhVxs9vgjYCUwa7Hgp3ndddMyro+d+EWgDDh3iGH3jGfSYyecwRQzrga8CFdHrNkbrq4GfJz3vIULCHNY5G+BYUwg/GX4dOBF4/WD76/O3TxnPMD9f1wIvRfcPB74EtANvS/q8WNLzzwdOi+5XAI8BFdHj+6PzUA00AV8AyoFLosdvJzQKfJzoM6hFy3haUL6eEPk62v65YZ7Hwf7WjwHHJL3Gh/k5S/w9Twb+OcX2pdFxipLWKXdnaFFL9zjioTVyMSEJvRdYaWapBoy8CLzHzIqB06J1VYR/5HOBRH/c5wj/mHD3DqAZuNnd9xL+EdYNsb++8bXTe2DGecB+j/oEA7cQ/rG+bYjj9d1vYp8/jZ77RUJ/47OHOEaveEZyzAHc7qH/3mPAzKT3eEjUEvVRwjmFYZ6zAd7vbkJyu83d/+bud49gfwPFMxy/AA6z0MXk/wD/BfwReLuZVQHNHmXlyCnAw9H98wlfbInWkwuB66Lz3Aj8xUM/8BVAtbv/0sNP7k8ysr+BSF5Qvp4Y+TqFAd/jIMd5PVDj7ol4VozwmK8G3uXunx/m85W7M6Qo2wFIPMysEJjr7uuB/2tm3yeM9F4OnJT8XHdvN7PZwCeB/0nswt13mNkK4CzC/8TnEwqtgy9Nut9F1D1poP0NEGryPhYAxUlxdZvZemDOYMcbiru3mdlaoBaYOoJjjPqYfXTR82+rHljn7l/v+6ThnDMzew+huE04ysPP0J4c6wj+BgPGMxR3f9HMniR0Malx97Vm9gtC69PFwO+S4k58NhL9BQ8ltEwl9rUreddJ9w/QWzfqBifjjPJ1r/c3EfJ1sgHP4yB/m+MI/8EYrTOBBWZW5e77BnhOd/QelLszSCdk/Cgmqb+tuz9F+Pnu6L5PNLMzgQ+5+xc9zESR7MPAqy1MXfU08JWhDjzE/gazkdBymjyVkhH6OqarBHghw8cYjq307veMmZ00gnP2W+D4pCXlKPgR7C9lPMN/O/wcuBxYmxRfIXAV4afYhJMI/VMTNgOnmlll0nEH7VcoMo4pX/c2kfL1gO9xkOPsjV5TwujcQPiV4D8G2N7kPQPglbszSEX3+PIPZjYr6fFs4A/R/Q7CwJGjCdNT1ZhZqZm9Pto+3cymEvoWfp/w09DDRD9XRgro8797M7Mh9tf32Ja0j98A+4B3RPuqjY6XaDEd6HgDmRk9Z3b02j8M4xjJ8Qx1zOT30VeqVg8jtF5VAneY2RstDNgZ6m9w8FjATHdfnbR0Jh2iMOn+YPtLjjtlPGZWYmbLo64jg/kl4WftmwHcvRm4E3jIe89a8gZC15OExPn+mYW5vj8CJOZq7fs36Pt3zuuZFEQGoHzNhMnXwz2PKY9DyKUlhP+YQfhVAzObNYzcbYRW5/cAF5vZ26PXftzMriW08v856fnK3ZmUiY7iWsZ+IQzOcGAPcBPwQ8LPU5Oj7f8MvAy8lfBz0YuEloXzCC0kPyQksN8SBnV0RPtrJfwjfBPQSRisMxf4dLT9/YPtr8+xrwIeIPSTWxRtW0JoIf0qoW/jsdH6AY83wPv36PVXEbo7LE7aNtAx5iXHM9Qxk89hn2O/lTCg8D+jc5F43eXR9jOAZ6O/zVej8zzcc/bWFO/1ouh4PwNmR+sG+5v22tcA8cwhXMThp8P4rH2jz+NLgNf1WXdDitedDqwCdgHXRuvOJXzWvgnMIpo2Dfh7wpfyf0Xv9c3Z/jemRUtcC8rXEylf1wMPElrrjx/iPQ72tz6P8AvjXwmD2l8APhS995S5m/BrSiehhXsK8BdCwX8F8AHCfOVfJXRtSbxGuTuDi0UnTwQzmw78vbt/IXpcQGh9+YC7X5vV4IZgZg7M99BHUkbJzP6fu49qOioRGTvK15JMuTs/qHuJJPsIoavBlKR1Z5I0QC4XJf2cpZ+yRsnMCs3sXcCvsx2LiAyL8rUod+cZFd2S7LuEvl/Pmdkm4C7gZXd/JLthDSwaWHJV9PAyC1PXychVAb9x97VDPlNEcoHytYByd15R9xIRERERkQxTS7eIiIiISIap6BYRERERyTAV3SIiIiIiGaaiW0REREQkw1R0i4iIiIhkmIpuEREREZEMU9EtIiIiIpJhKrpFRERERDJMRbeIiIiISIb9fws49re7Iwx/AAAAAElFTkSuQmCC\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/oceanus-saturn-probe.png',bbox_inches='tight')\n", "plt.savefig('../plots/oceanus-saturn-probe.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/oceanus-saturn-probe.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The max. heat rate is smaller than reported in the paper by about a factor of 3. The difference is attributed to different heating correlations used." ] } ], "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 }