{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 54 - ExoMars 2016 - Mars" ] }, { "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 Schiaparelli EDM, flown in the ExoMars 2020 mission.\n", "https://en.wikipedia.org/wiki/Schiaparelli_EDM" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"MARS\") \n", "planet.h_skip = 120.0E3\n", "planet.h_trap = 2.0E3\n", "planet.loadAtmosphereModel('../atmdata/Mars/mars-gram-avg.dat', 0 , 1 , 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Schiaparelli', 577.0, 75.0, 0.00, np.pi*2.4**2.0, 0.0, 0.6, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,5.900,0.0,-12.48,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (30*60.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXxcdfX/8dfplpbuGy3pQtpCFwq0YKAshbIUBRFFBFlEQYUgAl/cEAR+iIrK6hcQQSKLAl+k7Fi+LFL8yi7doCylFLrvpHuBbmnO74/PnWaSJmkmM5k7M3k/H4/PY+7cuXPvmVBOTj5z7r3m7oiIiIiISHpaxR2AiIiIiEghUGEtIiIiIpIBKqxFRERERDJAhbWIiIiISAaosBYRERERyQAV1iIiIiIiGdAm7gAypVOnTj58+PC4w0hLRUUFvXv3jjuMtOgz5AZ9htzQ2M8wbdq0le6e3x82Rb169fKSkpKU31cI/y6aQp+7ZdHnzm0N5eyCKazdnalTp8YdRlpKS0v1GXKAPkNuaEmfwcwWZCGcnFJSUtKk/76F8O+iKfS5WxZ97tzWUM5WK4iIiIiISAaosBYRERERyYCCKax79eoVdwhpKysrizuEtOkz5AZ9htxQCJ8h17TUn6k+d8uiz52/zN2ze0CzY4BrgFPdfX607iTgJqAL8ADwU3evjF7rBfwWWAu0BS51962191taWur50JcjIlKbmU1z99K446iLcraISE0N5eyszlibWR+gE3Bg0rqBwInAycCPgXOAHyW97RHgHne/FJgN/DprAYuItGDK2SIiqclqYe3uK4Cnaq3eHfieu09z9/uAPwFHApjZQcAwd38z2nYicJGZdcxWzCIiO+PubN26laqqqrhDyahcytnuUFkZhohIrsp6j7W7V9V6/kriK8TIEmBhtHwUMC9p2yXRYk5+ZSoiLdP+++9Pu3btmDFjRtyhZFyu5OzXX4e2bWHcuHT3JCLSfHLx5MUDgD9Hy/2AVbVe3wAU135TRUUFpaWl20d5eXkzhykiErRpE24JUJnCdGp5efn2fAXk89nXWc3ZWT4tSEQEaHzOzqkbxJjZHkCFuyemfRzYWGuzImCHE2F69+6dFxcVF5HCkyist27dITXVq6ysbPsZ8Ga2slkCa2bZzNlm0QFUWItIDBqbs3NmxtrMWgPnAZclrV4CdE/axoCOwNLsRiciUr+2bdsCqc1Y57ts52wV1iKSD3KmsAZ+Btzo7puT1r0ADEt6XgJsATQ1LSI5oymtIAUgqzlbhbWI5IOsF9bRDAaAJa27HHgL6GBmg83se2a2h7tPBZaZ2d7RpscCf3T3LdmNWkSkfk1pBckXuZKzt0chIpLDstpjbWadgTOjp2eb2W2ErxJ/U2vTWe5+T7R8KnCVmc0nXE/1imzEKiLSWIXaCpKLOVsz1iKSy7JaWLv7BuCOaCRcE4363rMA+H4zhyYi0mSF2gqSSzlbrSAikg9yqcdaRCQvFXIrSK5QYS0i+UCFtYhImgq1FSSXqLAWkXygwlpEJE2asW5+OnlRRPKBCmsRkTQVao91LtKMtYjkMhXWIiJpUitI81MriIjkAxXWIiJpUitI81NhLSL5QIW1iEia1ArS/FRYi0g+UGEtIpImtYI0P528KCL5QIW1iEia1AqSPZqxFpFcpsJaRCRNagVpfmoFEZF8oMJaRCRNagVpfiqsRSQfqLAWEUmTWkGanwprEckHKqxFRNKkGevmp8JaRPKBCmsRkTSpx7r56aogIpIPVFiLiKRJrSDZoxlrEcllKqxFRNKkVpDmp1YQEckHKqxFRNKkVpDmp8JaRPKBCmsRkTSpFaT5qbAWkXygwlpEJE1qBWl+OnlRRPJBm2wf0MyOAa4BTnX3+dG6DsD1wAagO3Cpu6/f2WsiIrmgkGescy1na8ZaRHJZVmeszawP0Ak4sNZLtwOvuvvlwFPAnY18TUQkdoXaY51LOVutICKSD7JaWLv7CkKi3c7MioEzgInRqknAV82spKHXshKwiEgjFGorSC7lbBXWIpIPst4K4u5VVrNZ7ghgpbt/Hr1eaWbzgHHA1gZem5/NuEVE6lPIrSC5krNVWItIPsh6YV2HfsCqWus2AMVAZQOv1VBRUUFpaen252VlZZSVlWU2UhGROjSlFaS8vJzy8vLE016Zj6rZxJKzdfKiiMSpsTk7FwprBzbWWldEmPlgJ69t17t3b6ZOnZr56EREdqIprSDJhaSZrWyWwJpHrDlbM9YiEofG5uxcKKyXEM4cT9YZWEpI4PW9JiKSEwq5FaQOseRstYKISD7IhetY/x9QHF2iCTNrCwwEXtrJayIiOaFQT16sRyw5W4W1iOSDrBfWVn0WjAG4+3LgaeDoaP144DF3X9LQa1kMWUSkQYV6uT3InZytwlpE8kFWW0HMrDNwZvT0bDO7zd0rgPOB681sFNAXSD6DpaHXRERiV6itILmUs3Xyoojkg6wW1u6+AbgjGsnr1wDn1vOeel8TEckFhdoKkos5WzPWIpLLcqHHWkQkrxVyK0iuaN06PG7bFm8cIiINUWEtIpKmQm0FySXt2oXHLVvijUNEpCEqrEVE0lSorSC5RIW1iOQDFdYiImlSK0jzKyoKjyqsRSSXqbAWEUmTWkGaX/SlAJs36wRGEcldKqxFRNKkVpDm17p19QmM+jGLSK5SYS0ikibNWGeH+qxFJNepsBYRSZN6rLNDfdYikutUWIuIpEmtINmRmLHevDneOERE6qPCWkQkTWoFyQ61gohIrlNhLSKSptatW2NmuDtVVVVxh1OwOnQIj59/Hm8cIiL1UWEtIpIB6rNufl26hMf16+ONQ0SkPiqsRUQyQO0gza9r1/C4bl28cYiI1EeFtYhIBugExuanGWsRyXUqrEVEMkCtIM1PM9YikutUWIuIZIBaQZqfCmsRyXUqrEVEMkCtIM1PhbWI5DoV1iIiGaBWkObXu3d4/OSTeOMQEalPRgprM/tyJvYjIpKv8qkVJF9zdnFxeFy6NN44RETq0yaVjc3sWOBKoC/VRXkroA/QIbOhiYjkj1xsBSm0nL3bbuFRhbWI5KqUCmvgfuD3wFtA4vZiBpyRbiBmNgK4EPgY2BMod/e3zawDcD2wAegOXOruutiSiOSUHG0FKaicnZixXrYsE3sTEcm8VAvrue7+h9orzWxWBmK5H/iauy8xs4HA88AI4HbgOXefEM2+3AmcnoHjiYhkTI62ghRUzk7MWC9bBlVV0EpnCYlIjkk1LZ1vZifXsf6rGYhlBNA5Wt4IdDWzYsLMysRo/STgq2ZWkoHjiYhkTI7OWBdUzi4qgl69YNs2zVqLSG5KtbA+FLjLzLYljSrgjgzE8mi0787AmcBFwBHASnf/HMDdK4F5wLgMHE9EpNAVXM4eMiQ8zpmTqT2KiGROqoX1r4EyQj/d4KRxdQZi+SGwCZgCfOrujwH9gFW1ttsAFNd+c0VFBaWlpdtHeXl5BkISEWkcdwfAzBq1fXl5+fZ8BfRqprAKLmfvsUd4/PjjpgcuIpKqxubsVHusnwWecvfNySvN7C+ph7iDDsDfCYn5ZjObBzjhK8ZkRcAOTYy9e/dm6tSpGQhDRCR1qRbWZWVllJWVJd6zspnCKricveee4VGFtYhkU2NzdqqF9UPAL8zsX0nrDDgJuDjVIGt5ADjN3deamRMS9o8IZ5Un6wzoYksiklNSLayzpOBydmLGelYmTr8UEcmwVAvr/0f4SvFswswEhCTdlzSStJn1Aka5+9po1TXRMRYAxWbWwd03mllbYCDwUlOPJSLSHHK0sC64nD16dHicPj1TexQRyZxUe6xvAAa4e4m7D4pGCWH2Ix2rgU1m1i9p3SrgHeBp4Oho3XjgMXdfkubxREQyKkcL64LL2cOHwy67wIIFsLK5GmhERJoo1RnrRe6+oY71fdMJwt2rzOxE4Cozm0a4K9gl7r7ezM4HrjezUdFxytI5lohIc8jRwrrgcnbr1rDffvDaazBtGnzpS5ncu4hIelItrH8IvJG8wsxOAW4G7kknEHefAZxXx/o1wLnp7FtEpLnlaGFdkDl7zJhQWL/yigprEcktqbaC9DCz/wdgZruY2V3AXYTLLYmItFg5WlgXZM4+Omo0mTQp3jhERGpLacba3Y83s4PM7I/Al4DlwChgXXMEJyKSL3KxsC7UnH344dCmDUyZAuvWQdeucUckIhI0WFib2eH1vLQA+Ai4EdgdOBH4cWZDExGRVLSUnN2pExx0ELz6Kvz73/C1r8UdkYhIsLMZ6+uBEcBaoKqO1+8hXLppN/I4SYuIpGvr1nAPlDZtUj11JaNaTM4++uhQWE+apMJaRHLHzn4D/BqY6u6fNLSRmZ2cuZBERPLPpk2bAOjQoUOcYbSYnD1+PPzqV/DMM3DrrZBDHTgi0oI1ePKiuz+zswQdbfdo5kISEck/uVBYt6ScffDBsNtuMHcuvPlm3NGIiASpXhVERETqkCis27dvH3MkLUPr1nD66WH5gQfijUVEJEGFtYhIBqiwzr5vfSs8PvwwVFbGG4uICKiwFhFJ27Zt29iyZQtmRrt27eIOp8XYbz8YOhQqKuBf/4o7GhGRDBTWZrafmT1gZj81s5L0QxIRyS+bN28Gwmx1Ll3Hui6FlLPN4LTTwvLDD8cbi4gIZGbG+gfAeHe/CTjEzGI9JV5EJNsSbSBFRUUxR9IoBZWzv/nN8Pjoo7BhQ7yxiIhkorC+EPhCtPwQ0DYD+xQRyRvr1oUbGXbp0iXmSBqloHL2yJEwdmy4A+Pdd8cdjYi0dCkX1mZ2nJl9K1o+ENjT3ZcAuHuVu6/PcIwiIjlt7dq1AHTv3j3mSHbUEnL2JZeExz/8AaL79IiIxCKlwtrMfgs8AnwdwN0nA98xszOaITYRkbywZs0aIPcK65aSs7/yFRg+HBYtggkT4o5GRFqyVGesDwH6A68krXsWuDZjEYmI5JlEYd2tW7eYI9lBi8jZrVrBz34Wln/3O6iq62buIiJZkGph/YK7rwU8aV3e3xpXRCQdOdwK0mJy9re/DQMHwgcfwGOPxR2NiLRUqRbWc83scmCYmZ1mZo8BFwDXZT40EZH8kKutILSgnN2uHfziF2H5N7/RrLWIxCOlwtrdHwL+BbQGvgOsBb7o7n9qhthERPLCypUrAejRo0fMkdTU0nL2d78LxcXw7rvw9NNxRyMiLVGbVN/g7v8B/pN4bmatzexod38xo5GJiOSJJUuWANCvX7+YI9lRS8rZRUWh1/onP4HLL4cvfxnapPxbTkSk6RpMOWb2D6DTTvbRC1gJZCRJm9lBwFjgY+BV4DPgemAD0B24tBAuDyUihWPp0qUAFBcXxxqHcjacfz7cdhu8/z78+c9w4YXZOrKIyM5nrN+JtvkQqALOBd4CpiZtMwjIyPefZvZdYIi7X5m07l7gOXefYGbHAncCp2fieCIimZBDM9YtPme3bw833QRf/zpcdRWcfjr07Jmto4tIS7ezwvoGwKKzyjGzIe5+Ve2NzOypdAMxs3HAt4BjktYVA2cQTrYBmAQ8ZmYl7j4/3WOKiGRCrsxYo5wNwNe+BkcfDS++CJddBn/5S7aOLCItXYMnL7r7ukSCjvQxs9bJ25jZEcCBGYjlD8AHwB/N7DkzOxQ4Aljp7p9H8VQC84BxGTieiEjaNmzYwIYNG2jfvn3s17FWzg7M4I9/DFcKuesueOGFbB5dRFqyVC+39ygww8z+bGY3mNnTwAvArekEYWbDgP2Av7j7hYTev+eAAcCqWptvAHaYFqqoqKC0tHT7KC8vTyckEZFGmTt3LgC77747Ztbo95WXl2/PV4S+5+bQYnP2iBFwZdSg8o1vwCuvNLy9iEhDGpuzUzpf2t1fMLMvAj8E9iVcuulkd0/3a8WRwGp3fyd6fhvwS8JNDTbW2rYI2Fp7B71792bq1Km1V4uINKsPP/wQgGHDhqX0vrKyMsrKygAws5UZDwzl7F/8AmbOhIcegvHj4e674cwzm/WQIlKgGpuzm3K5vaVA8okqrczsC+4+rSmBJsWxPRZ332hmHxGuvVr7jgudgaVpHEtEJGOaWlhnS0vO2W3awP33Q+/eoTXk29+Gjz6Cq68O7SIiIpmWUmFtZvfUsbo3sA04MY043gG6mlkfd18RraskJONiM+sQJe62wEDgpTSOJSKSMblcWCtnh+L61lth6FC4+GL49a9h+nS4917o1VwNOCLSYqXaY70nsBBYkDSKgHRmPnD3WcCzwMkAZtYNaAf8D/A0cHS06XjgMXdfks7xREQy5e233wZg5MiRMUdSJ+XsyIUXhrsxdusWHvfdVyc1ikjmpdoKco67f5i8wsyGA2dnIJbvALeYWTvCDMfp7l5pZucD15vZKKAvUJaBY4mIpG39+vXMnDmTtm3bMnr06LjDqYtydpLjjoMZM+Bb34JXX4UvfhHOPhtuvFHXuhaRzDB3T28HZoOAKe4e65dqpaWlrpMXRSSbXnzxRcaPH88BBxzA5MmTm7wfM5vm7qUZDK2hY7X4nF1ZCddfH9pCNm8OPdg33xxuJqPeaxHZmYZydkqtIGY2z8zmJo0FhOuYPpGJQEVE8slrr70GwJgxY2KOpG7K2XVr0wYuvxzeeQfGjYOKijCLvd9+8NhjsG1b3BGKSL5KtRXkBeBBwiWVINwyd4W7z85oVCIieeCZZ54B4JhjjtnJlrFRzm7A0KHwr3/BPfeE25/PmAEnnwwDB8J558H3vw99+sQdpYjkk1RPXnzA3f/t7i9F4xXAzGyP5ghORCRXVVRUMHnyZNq1a8dRRx0Vdzj1Uc7eiVat4JxzYO7ccEm+IUNg4UK44goYMCC0hzz/fGgfERHZmVQL62PrWLcCeCQDsYiI5I0nn3wSd2fcuHF06tQp7nDqo5zdSO3bhyuHzJ4Nzz4LX/1qaAl56CE49thQZP/0p2FWW0SkPo0qrM3sQjNbC1xqZtuSB+H2tSt2sgsRkYJy9913A3BmDt7KTzm76Vq1CoX0U0+FWexf/Qr22AOWL4c//AFGj4bhw+Gyy+DNN6GqKu6IRSSXNPqqIGa2L3AWcGutl7a4+7JMB5YqXRVERLLlnXfeYdSoUXTu3Jnly5ezyy67pLW/5rgqiHJ25riHIvr++2HCBFi1qvq14mI48UQ4/vhwImTHjvHFKSLZ0VDObvTJi+7+jpld5u5b6zjAQHdfmE6QIiL54ve//z0AZ511VtpFdXNRzs4cMzjooDBuuQVeeQWeeAKefBIWLYLbbw+jXTsYOzZcH/uLX4RRo8IMuIi0HA3OWJvZz4Bn3H1m9PwKoHXtzYAj3f2I5gqyMfJp9kNE8tcHH3zAyJEjadOmDXPmzGHAgAFp7zNTM9bK2dnlHm6P/uST4QTHqVPDuoTevcMsdmKMHKlCW6QQpDNjPQjolvR8f6AjsJzqyzcZ4e5aIiIFzd25+OKLcXe+973vZaSozjDl7Cwygy98IYzf/Ca0iLz4Yiiy//lPWLwYHn00DIAePeCww6oL7VGjoHXtP3tEJK+ldOdFM9sfeM/dt9Rav1dihiQuhTD7ISK5bcKECZx22ml0796dDz/8kN69e2dkv81150Xl7Pi4hyuMvPQSvPxyeFy8uOY2XbrAoYfCmDFw4IFh6NbqIrkvIz3WkX3cfXqtnQ8DhgCxJmkRkea0YMECzj//fACuvfbajBXVzUw5OyZmMGxYGGVlodCeN6+6yH7ppfD82WfDSBgypLrQHjMmzGp36BDf5xCR1DSqsDazwwhfHx5jZvNqvdyDcNb5xAzHJiKSEzZv3szpp5/OmjVrOP744znnnHPiDqlBytm5xwwGDw7j7LPDukWL4I03YPLkcNWRadNgzpwwHnwwbNOqVbi83+jRNUd+/F0n0vI0dsZ6LXAfoX/v0FqvbQHuymRQIiK5oqqqirPOOos33niD/v3787e//Y1WuX8GmnJ2HhgwIIxvfjM837oV3n8/FNmJYnvWLJg5M4xEsQ3hMn+JInuffcKJkUOHQlFRPJ9FRIJGFdbu/q6ZHQGc4u5KyCLSIlRVVXHxxRczYcIEOnfuzMSJE+mZB02wytn5qW3b6mL5vPPCuo0bQ7H99tvVY8YMWLo0jGeeqX5/69ahlWTkSNhrr+oxbJjaSUSyJZXrWK+jnlkOMxvr7q9mLCoRkZht27aNc889l3vvvZe2bdvyxBNPMHr06LjDajTl7MLQoQOUloaRUFUV7gqZKLTfey/MaM+ZE06YnD07XGc7oVUrGDQoFNhDh4aRWC4u1iUARTKpwcLazP4BdNrJPgwYCeyaqaBEROK0bt06zjzzTJ5++mk6dOjAE088wdFHHx13WDulnN0ytGoVbrO+xx5w8snV6zdtCkX1++9Xt4/MnAkffVTdu508ww2wyy6w557VBXdy0d29e3Y/l0gh2NmM9YeEXr0lVF8DtTYDvpLJoERE4vLBBx/w9a9/nQ8//JDu3bvzj3/8g7Fjx8YdVmMpZ7dg7dvDvvuGkWzzZvj44+rZ7MT48EOoqAitJTNm7Li/Xr1qFtyJonvIELWWiNRnZ4X1dcAGd9/c0EZmprPLRSSvuTu33347l1xyCRs3bmSfffbhiSeeYMiQIXGHlgrlbNlBUVHoux45csfX1qwJM9q1i+7Zs2HlyjBef73me8xg4MCaBfeQIeGKJyUlKrqlZWuwsHb3lY3cz9fRWeYikqfmzp3L+eefzz//+U8Avv3tb3PHHXfQsWPHmCNLjXK2pKp79+qb0yRzDydH1lVwz50LCxaE8cILO+6zb9/Q0z14cHhMjMGDoX9/3W1SCluqN4ipwcwGAhcAF6IkLSJ55vPPP+e6667juuuuY/PmzfTo0YM777yTk5MbVwuIcrY0lhn06xfGkUfWfG3r1nBzm9rF9rx5MH8+LF8exhtv7LjfNm3CbHftonvgwDB2202Ft+S3lAtrM2sNfBU4FzgGWEHo6UubmbUDpgAXu/u/zawDcD2wAegOXOru6zNxLBFpuSorK3nggQe4+uqrWbBgAQBnnnkmN9xwA3379o05usxqzpwd7V95u4Vp27a6BaS2bdtgyZJQZCeK7cSYOxeWLQuPc+fWve/WrcOsdqLQrmt06dK8n08kHY0urM1sEHAO8F2gJ7AAOAF4Hjg4Q/FcApQkPb8deM7dJ5jZscCdwOkZOpaItDDbtm3j4Ycf5uqrr2b27NkA7Lvvvtx2220cdthhMUeXWVnK2aC8LUlat64ugMeN2/H1jRtDC0ly0b1gASxcGMby5dVtJvXp2nXHYrt//zASs+y77NJ8n1GkITstrM3sFMJMx1HAx8AfgL8Bl7v7c9Fmr9fz9kYzs0OAZcCa6HkxcAbha0uAScBjZlbi7vPTPZ6ItByfffYZf/3rX7n55pv5+OOPARgyZAi//OUvOeOMM2hdQN89ZytnR8dS3paUdOgQbtE+fHjdr2/eDIsXVxfatceCBbBuHbz7bhj16dEjFNjJBXft5a5dQ8uLSCY1ZsZ6I+EWuK8D57j7bAAzq+9STikzs47Aye7+EzO7Klp9BLDS3T8HcPdKM5sHjAPmZ+rYIlK45s2bR3l5OXfeeSdr1qwBoKSkhCuvvJLvfOc7tG3bNuYIm0Wz5+xof8rbknFFReEKI/VdjMcdVq/esdhesiQU5IsXh+XVq8NoqPju2HHnxXfv3rqBjqRmp4W1uz8NPB2d9HKumZUATwHbfyOZmbl7Okn7UuDaWuv6AatqrdsAFNe1g4qKCkqTbk1VVlZGWVlZGiGJSD7avHkzTz75JHfddReTJk3avv7ggw/mJz/5CSeeeCJt2qR13nZGlJeXU15ennjaK1P7zVLOhjTztnK2NIUZ9OwZxn771b1NVVW4TGCiyE4uuBPLixfDZ59Vn3xZn7Ztw90pGyq+d9stbCeFrbE521LNrWbWBjgR+BbwHnA38EV3L2/wjfXv70tAkbv/I3o+HzgbKAVOcfcxSdtOBx509xtr76e0tNSnTp3alBBEJM9VVVXxxhtv8NBDD/H3v/+dVatCbde+fXu+8Y1vcMEFF3DwwZlsK84sM5vm7qU737JJ+85ozo72mXbeVs6WOLnD+vV1F9zJz1ev3vm+zMIlBhua/Vbfd2FpKGenPG3j7pXAo8CjZrYH8GNCP19Tk/TPgP2sutGpO2F25Q/RcrLOwNImHkdECoi7M336dB566CEmTJjAokWLtr82atQozj33XM444wy6t/D7MjdDzgblbclzZqHHumvXum+ck7BxY3WhXd/s9/Ll4Wony5ZBQ38rqu+7ZUjr+1B3/xi42MympbGbM4CipOdvAD8B3gQuMbMO7r7RzNoCA4GX0jiWiOSxLVu28PLLLzNx4kQmTpzIvHnztr82YMAATj31VE4//XT2228/TL+ddpChnA3K29JCdOgAe+wRRn22bg3FdUPFt/q+W46MNBq6+31pvLci+bmZbQMq3H2hmT0NHA08DYwHHnP3JWkFKyJ5ZeXKlTz33HNMnDiR5557jvXrqy+J3KdPH0455RROO+00Dj74YFrpt02jpJOzo/crb4tE2raFAQPCqI/6vluO+M/gadj5wPVmNgroC+jMFpEC9/nnn/PKK6/w4osvMmnSJN56660ar++9996ccMIJnHDCCRx44IEFdam8AqG8LVJLq1aw665h7L9/3dvU1/dde3nVqp1f61t93/HJucLa3UuSltcQegFFpEBt2rSJKVOm8PLLLzNp0iRef/11tmzZsv31oqIiDj/8cL7yla9wwgknMGjQoBijlboob4ukL9W+74ZmvjPV992vH3Trpr7vVORcYS0ihW316tW89tprvPrqq7z66qtMnTq1RiFtZpSWljJ+/HjGjx/PIYccQocOHWKMWEQkd2S773uXXaqL7OSCO3m5b99w101RYS0izaiyspL33nuPKVOmMGXKFF5//XXef//9GtuYGfvuuy+HHnooRx99NEceeSQ9evSIKWIRkfyXqb7vJUvg00/ho4/CqE+rVqGvu66iO/l5S2g9UWEtIhlRVVXFxx9/vL2InjJlCtOnT2fTpk01tisqKmLMmDGMHTuWsWPHcvDBB9OtW7eYohYRaZka0/cNNfu+kz0YFsEAACAASURBVC89mPz8k0+qnzekW7e6i+7dd4eSEhg4MMzI5zMV1iKSss2bNzNz5kzefvttZsyYwYwZM3jrrbdYt27dDtsOGTKEAw44gAMOOICDDjqIL3zhCxQVFdWxVxERyTVdusBee4VRny1bQj93fQX44sWwdCmsXRvGe+/Vv68+fUKhnSi2ay937pzpT5hZKqxFpEGffPLJ9uI5UUjPmjWLysrKHbbdbbfdthfRBxxwAKWlpfTs2TOGqEVEJFvatasugOvjHlpPahfdixaFK5zMnx+WV6wIY/LkuvfTowcMGgR77glDh9Z8zIX7gamwFhHcnU8++YSZM2duH++//z4zZ86koqJih+3NjGHDhjFq1KjtY/To0fTr1y+G6EVEJNeZhZvb9O4No0fXvc22bWHme/786mI7cWnBxHLihMtpddzmqlevmoX2sGGw997hRM82Wap4VViLtCDuzrJly/jggw9qFM8zZ85k1apVdb6nU6dO7LvvvtuL51GjRrH33nvTsWPHLEcvIiKFrHXr6sv9jR274+tVVaGfe+7ccDLl7Nk1H1euDOONN2q+r6gotLLss08YY8fCAQc0z5VMVFiLFKB169Yxe/bsOsenn35a53u6dOnCXnvtxV577cXIkSO3Lw8YMEC3BxcRkdi1ahUu7de3LxxySM3X3EMfd3KxPWtWuJTgwoXw1lthJPToARdfDFdemdlbxKuwFslTmzdvZs6cOXUWzytWrKj3fd27d99eNCcX0cXFxSqgRUQkL5lVX2nkyCNrvrZuXThh8t134e234YUXwqz3L38JI0bAKadkLg4V1iI5bMOGDcyZM2f7mDt37vblBQsWUFVVVef72rdvz5577smwYcMYOnRojaGTCUVEpCXp2hUOPTSMTz+FP/8ZLrkkvJbpi1SpsBaJkbuzYsWKGsVz8qjrxMGEVq1aMWTIkB0K56FDh9K/f39aZfK7LRERkTy1YAH861/wz3/CP/4Bn38e1l94IZxwQmaPpcJapJlt3LiRBQsWMH/+fObNm7fDDPRnn31W73uLiooYPHgwQ4YMqTEGDx7MoEGDdD1oERGRiHu4lfv06eGqIYlR+8Y1hx4KV1wBxx2X+RhUWIukacuWLSxcuHB74Vz7cfny5Q2+v3v37jsUzolRXFysmWcREZEkmzaFWejZs+GDD8KYNSs81nGfMrp1gyOOCL3XX/kKDB7cfLGpsBbZicrKSpYsWVJn0Txv3jyWLFmCu9f7/jZt2jBw4EBKSkoYNGjQDjPQ3XPhivYiIiI54rPPqm8kM38+zJsXRmJ56dL639utG4waBV/4QhilpeE61tmao1JhLS1eVVUVy5Ytq7Nwnj9/PosWLarzLoMJrVq1YsCAAdsL59qPxcXFtMnWlelFRERy1JYtUFERCuOlS6vvvpgYiXV1zTona90aBgwIN4EZMQKGDw+PI0bArruGK4TERb/tpeC5O8uXL2f+/Pk1CubE8sKFC9myZUuD+yguLq63cO7fvz/t2rXL0qcRERHJHZ9+Gm7asmJF9WN9y2vWNG6fRUVQXBwunVdSEm5hPmhQ9XL//tm7k2KqcjQskcZzdyoqKuosmufPn8+CBQvYtGlTg/vYdddd6y2cBw4cSPv27bP0aUREROKxbVu4XXjiDoarVlUv1x6JYjlxhY3GaN063NK8b9/qa07361ddRCdGjx7xzjqnQ4W15Dx3Z9WqVXXOOCfG5zv5P7tnz57bi+XkwrmkpITdd99dt+cWEZGCUlUFa9c2XBzXHmvWhCtrpKJ9e+jTJ4xdd214uUeP7PU6x0WFtcTO3Vm7dm29M87z58+v9zbcCd27d9+hYE4833333encuXOWPo2IiEhmbd0aZpJXraoulOtaTn6+enUorlPVowf06tXw6Nmzulju1Cl/Z5ebQ84U1mZ2EnAT0AV4APipu1eaWS/gt8BaoC1wqbtvjS9SaYr169fXWzTPmzeP9evXN/j+Ll261DvjXFJSQteuXbP0SUQElLNFmmrjxroL4YaWd3YyX326dNl5kZw8unfP3d7lfJETPz4zGwicCJwMjATuABYBNwKPAJe5+5tm9gPg18Av4opV6rZx40bmzp3L3Llza7RoJAroNTs5Y6Fjx44MGjSo3uK5W7dumP4kFskJytkioWViw4adzxzXXk6lJznBLMwk9+xZPWPc0HLPnmF7nVeffTlRWAO7A99z90pgmpntDRxpZq8Cw9z9zWi7icCHZnaNu9d/uzppFp999hlz5szh448/3mEsWrSowffusssuO7RoJC/36NFDhbNI/lDOloJSVRX6i1NptVi1KrRopKpt250XxbVf69at8HuTC0VOFNbu/kqtVUuAzsBRwLyk7ZZExVcp8FLWAmxBNmzYsL14/uijj2oUz0sbuCJ7mzZttt/8pPaVNUpKSujdu7cKZ5ECoZwtuWzr1tRbLdasaVo/8i677Lworv2aepILW04U1nU4ALgBKANW1XptA1Bc+w0VFRWUlpZuf15WVkZZWVlzxpi33J2lS5fywQcf1BizZs1q8Pbbbdu2ZfDgweyxxx7sueee7LHHHtuXBw4cqJugiKSgvLyc8vLyxNNeccaSAcrZ0iw+/zz1Inknp+zUq2vX1GaRe/YMV8SQlqGxOTvnKiEz2wOocPcZZubAxlqbFAE7fPnSu3dvpk6dmo0Q80ZlZSXz5s2rs4Cu72TBdu3aMWTIkDqL5wEDBtC6dessfwqRwpRcSJrZypjDaTLlbGkM91DwptqPvLH2v6ZGaNWqZkHcmL7kHj100p40rLE5O6f+GZlZa+A84LJo1RJgaNLrBnQEGrhLfMu0evVq3nnnHWbMmLF9vP/++2zevLnO7Xv06MGIESNqjOHDhzNw4EAVzyLSKMrZLVtlZbhJyLJlYSxfXr28YkXNm4ysWhW2T1W7do0/WS/xvGtX9SNLfHKqsAZ+Btzo7olq8AXgB0mvlwBbgBY7zeHuzJkzh2nTpjFjxoztxfTixYvr3L5///47FNAjRoxQz7OIZIJydgHatAmWLq27YE5e98knqd1MpFOn1GaRe/aEjh3Vjyz5JWcKazO7nJB8O5jZYOAI4GVgmZnt7e7vAccCf3T3LfFFml0rVqxgypQpTJ48efuo69J1HTp0YJ999mHUqFHsu+++2x91fWcRaQ7K2flp69ZQNC9aVP+oqGjcvszCTUJ22y2Mvn1rLiffTKRnTygqat7PJpILcqKwNrMrgd/UWj3L3e8xs1OBq8xsPtAJuCLb8WXLtm3bePfdd3n55Zd59dVXmTx5MgsWLNhhuz59+lBaWsro0aO3F9B77LGHWjhEJCuUs3NXZSUsXAhz58KcOeFx3rzqonnZsp1f/aJNGygurrtgTl63667h0nEiUi0nCmt3vwa4pp7XFgDfz25E2bF161beeustXnrppe3F9Nq1a2ts06lTJ0pLSznwwAM58MADOeCAAxgwYIDaOEQkNi01Z+eKzz6Djz4KhXOieE48LlgA27bV/16zUDQPGFD/6NMHNE8j0jQ5UVi3JIsXL+bZZ5/l2WefZdKkSWzYsKHG6yUlJRx++OEcfvjhHHTQQQwfPlwz0SIiLYx76GOeNWvHsXBhw+/t3x8GD4YhQ8LjoEEwcGAomouLdTc+keakwrqZuTtvvvkmjz/+OM8++yzvvfdejdf33HNPxo0bx7hx4zj88MMZOHBgTJGKiEgc1q+HGTPg7bfDeO+9UEDXdz3mtm1D0bznnjUL6CFDoKRE11YWiZMK62aQKKYffvhhHn300Rq3++7UqRNHH300xx13HMceeyy77757jJGKiEg2rV0L//kPTJ0aiui33gotHHXp1g1GjIDhw2uOQYPU2yySq1RYZ9DChQu59957uffee2ucdNi/f39OPvlkTjjhBMaOHUs7fQ8nIlLw3GH2bHj9dXjjjfA4c+aOl6hr1w722QdGjw5jn31CQd27ty41J5JvVFinadu2bfzjH/+gvLyc559/Ho8yZv/+/TnllFM45ZRTGDNmDK10tXoRkYK3bBlMmgQvvBAely2r+Xq7dlBaCmPGwH77hUJ6+HDNQIsUChXWTbRx40b++te/ctNNNzFnzhwg3A78pJNO4pxzzuHII49UMS0iUuDcYfp0eOwxePppePfdmq/vuiuMHQuHHBLG/vvres4ihUyFdYo2bdrEbbfdxvXXX09FdBX9wYMHc9FFF/Htb3+bnj17xhyhiIg0J/fQG/0//wOPPw7z51e/1qEDjBsHxxwTxt57q51DpCVRYd1I27Zt4/777+eqq67afjJiaWkpP//5zznppJN0STwRkQK3alUopu+5J1zFI6FvXzjpJPj61+GwwzQjLdKSqbBuhHfeeYdzzz2XyZMnA7Dvvvty7bXXcuyxx+pGLSIiBW76dLjpJnj0UdgS3Zy9Z0844ww49VQ4+GBQ55+IgArrBm3dupVf/epXXHfddVRWVtKvXz+uvfZazjjjDPVPi4gUuBdfhN/9Dv71r/DcDI49Fr7/fTjhBM1Mi8iOVFjXY+HChZx66qn85z//wcy44IIL+N3vfkeXLl3iDk1ERJrRe+/BJZfAc8+F5506QVkZ/Nd/gW49ICINUWFdh5deeomTTjqJ1atX079/fx588EEOO+ywuMMSEZFmtGkT/PKXcOONUFUFXbrApZfCD38YbtYiIrIzKqxreeSRRzjzzDPZsmULxx13HPfddx+9evWKOywREWlGb78deqY/+CD0S19wQSiye/eOOzIRySdqFE5y3333ceqpp7JlyxYuvPBCJk6cqKJaRKTAPfooHHpoKKqHDYPXXoPbblNRLSKpU2Edefzxx/nud7+Lu3PNNddw66236hJ6IiIF7s474ZRT4PPP4ayzwvWpDzoo7qhEJF+psAYmT57M6aefTlVVFVdddRVXXHGFLqMnIlLg/vxn+MEPwvLvfw/33htu8CIi0lQtvsf6k08+4Rvf+AZbtmzhBz/4AVdffXXcIYmISDO7+244//yw/N//DT/6UbzxiEhhaNGFtbtzzjnnsHjxYg455BBuueUWzVSLiBS4F1+E884LyzffDBdfHG88IlI4WnQryMMPP8zEiRPp0qULEyZMoF27dnGHJCIizWjVqnD1j23b4LLLVFSLSGblxYy1mXUArgc2AN2BS919fTr7/PTTT/lR9N3fDTfcQP/+/dOOU0REmidnZ8rPfw6ffAKHHw6//W3c0YhIocmXGevbgVfd/XLgKeDOdHd4ww03sHz5csaMGcM555yTdoAiIrJdxnN2JsyYEU5QbNcO/vKXcL1qEZFMyvkZazMrBs4ALohWTQIeM7MSd5/f1P2OGzeO0aNHc9NNN9FK2VVEJCOaK2dnQo8ecOaZ0KsXDB0aZyQiUqhyvrAGjgBWuvvnAO5eaWbzgHHA/Kbu9KijjmL69Ok6WVFEJLOOoBlydiYMGAD33QfucUYhIoUsH6Zq+wGraq3bABSnu2MV1SIiGddsOTtTlPpFpLnkw4y1AxtrrSsCtiavqKiooLS0dPvzsrIyysrKmj86EZEmKC8vp7y8PPG0V5yxZJhytogUnMbm7HworJcQzipP1hlYmryid+/eTJ06NWtBiYikI7mQNLOVMYeTScrZIlJwGpuz86EV5P+A4ujyTZhZW2Ag8FKsUYmISF2Us0Wkxcr5wtrdlwNPA0dHq8YDj7n7kuTtKioqsh1axiV9xZC39Blygz5DbiiEz5CqxubspmqJP1PQ525p9LnzV84X1pHzga+Z2RXAl4EdGvFWrsz/b1IL4R+UPkNu0GfIDYXwGZpopzm7qVrqz1Sfu2XR585f+dBjjbuvAc6NOw4REdk55WwRaanMC+SCnma2Afgw7jjS1AvI96l3fYbcoM+QGxr7GXZ3997NHUwuMbMKYEET3loI/y6aQp+7ZdHnzm315uyCKaxFREREROKULz3WIiIiIiI5TYW1iIiIiEgG5H1hbWYdzOyPZvY7M7vDzLrEHVNTmdkxZvammZXEHUuqzOwkM5tnZqvM7BYzy4sTY5OZ2Tgz+9DM1pjZzXHHkw4za2dmM8zsiLhjaQozu8bMPBoz4o6nqczsIDP7mZmdaGaFdHfFrCukXJ+KQsit6cj3XNZULTF3mNkIM/uTmf3YzG43s9Fxx9QUhfA/6O3Ac+4+wcyOBe4ETo85ppSZWR+gE3Bg3LGkyswGAicCJwMjgTuARcCNccaVCjPrChwKHAQcAjxpZk+7+6R4I2uyS4CSuINoCjPrSDiB5ZhoVVNOcIudmX0XGOLuV8YdS4EoiFyfikLIrRmQt7msqVpw7rgf+Jq7L4n+7T8PjIg5ppTldWFtZsXAGcAF0apJwGNmVuLu82MLrAncfYWZPRV3HE20O/A9d68EppnZ3sCR5Ffy3wL83sPZvP9rZu8A22KOqUnM7BBgGbAm7lia6FxgFvCqu2+KO5imMLNxwLeo/uNA0lBIuT5FhZBbm6wAclnKWnjuGAF0jpY3Al1jjKXJ8r0V5Ahgpbt/DhAln3nAuDiDaip3r4o7hqZw91ein33CEmBhXPE0hbtvjIrqxIzpB8C/Yw2qCaLYT3b3e+KOpSnMzIDTgJuAFWaWrzOSfyD8G/qjmT1nZofGHVCeO4ICyvWNVQi5tanyPZeloSXnjkeBu8ysM3AmcFHM8TRJvhfW/YBVtdZtAIpjiEWqHQD8Oe4gmsLMDgOeBdoCu8QcTlNcClwbdxBN5cFBhFaQPwL/Y2bHxRxWSsxsGLAf8Bd3vxB4EXguaveSplGuD/I2tzZBXueyplDu4IfAJmAK8Km7PxZzPE2S74W1E74uSFYEbI0hFgHMbA+gwt3z9YSzj4G/AkcDN8QbSmrM7EvAVHf/JO5Y0uXua6L+wt8DP4o7nhSNBFa7+zvR89sAA74WX0h5r8Xn+gLIrY1WSLksRS09d3QA/g48CNxsZuNjjqdJ8rrHmvC1WPda6zoDS2OIpcUzs9bAecBlccfSVO6+DLjHzKoIJ83kk58B+4VuCiD8v/GUmf3W3a+PL6y03Er+teS0ISm3uvtGM/sI6BlfSHmvRef6QsitKSrEXNYYLT13PACc5u5rzcyBv0fnUXwWd2CpyPcZ6/8Dis2sA4CZtQUGAi/FGlXL9TPgRnffHHcgGTCV8Ms8n5wBjE4aS4FzyO+vjquA6XEHkaJ3gK61vr6tBGbGFE8haOm5vpBya2MUYi5rjBabO6JLCo5y97XRqmuA9eThVUHyurB29+XA04Sv7QHGA4+5e74VRMD2E7cgfPWTV8zscuAtoIOZDTaz70VfXeYFM9vFzPZLWvVl4Ja44mkKd69w98WJQbiqSYW7r487tsYys2IzO8PMWkX/P/wUuCLuuFLh7rMIffonA5hZN6Ad8L9xxpXPCi3XpyLfc2tTFEIua4oWnjtWA5vMrF/SulXA7JjiabJ8bwUBOB+43sxGAX2BspjjaZKks2ABzjaz29y9Is6YGsvMrgR+U2v1rDw7m3tvwmX2ZgH/ASa7e0tIZrmmO+GEpcuBV4Gb8/Ryat8BbjGzdoSZ1dNrXd1BUlcQuT4VBZJbJTUtMne4e5WZnQhcZWbTgD7AJfn4x5RFVxgTEREREZE05HUriIiIiIhIrlBhLSIiIiKSASqsRUREREQyQIW1iIiIiEgGqLAWEREREckAFdYiIiIiIhmgwlqkHmZ2tJm9bmZH7GS7vc3sYTP7fxk89g2Z2peISBzMbF8zeySTubGe4zxtZqc15zGai5mNN7M3dvZ7RvJHIdwgRgqcmY0Afku4M9MBwL7At939gWY+9CJgVCO2W0i4sUnrTBw0ul3z1kzsS0QEwMwOAH4PHA78E5hHuNFOW+AOd3++GQ77MdCVDOXGBtwOvNccO47y8UmE26nPBt6IXioCvgEcmOZNrBr7e0byhApryQdPAt909xlm1gq4NxsHdffZZraqEdutN7NlGTz0OODlDO5PRFo4d59iZg8Bw939K4n1ZnY88ISZ/czdb8/wMT/PcG4EwMxOAN5NFLTu/kymj5Hg7luBCWZ2EfBvd78yKY6nmrJPM9sVONLdJ7j7h2a2MkPhSg5QK4jkNDPrAQwFPoNw21PgKiBbtwytyvB2jXEk8H8Z3J+ICMAOt8Z29/8Frgb+28z6NsMxM5qrzaw/YYY62+q6rfgkYEkqOzGzIuA+oEPS6kz+/pCYqbCWnObuq4EPgefM7KBo3QJgspl1M7PfmNlkM/uqmS03sw/MbP/E+83sdDO7wsz+bWa/S1pfamZXm9kEM3vczLpG69ub2Z/M7Cozu5vwNWZKzOznUUzfNLMDzOzJaH+/iWKcYGa9o/XrzOyXtXbRzt03R/v6qZn9xMxeM7NrU/4Biojs3N2ElpBTAMxsiJn92szuMrMXo2IWMyuKctkVZvaKmZ2a2IGZ/dDM/p+Z/cfMLqzvQHXlZDM7OcqZZ5vZR2Z2rZkVm9m9UT59xcyOinZxBtAfuNTMTozOhZlkZmclHeNcM/utmd1nZn8zs45m1tXMfmlmU83sKDP70MzmmllJU35g0c/kiGhGGzM7ycxuMbObzOw5M+sUfYZbzGxi9PtmWfQz3gs4xcx+krTL7tHvog1mdm49x2xlZteZWVl0DDezJxL/fSRHuLuGRk4PYCSwgPBX/T3ArtH61oQktZrQ69YNeAF4n/BH48HAldG23Qmz3scDXYCHkvb/GnBDtPzfQFm03AfYRkieO4vxr4RZn10JvXi7ROuLgOeBp4Biwi+ErcAvgV2AL0TPE9sn+scTn/v2aLkzcFnc/y00NDTydwBnA4vreW0FcGuUVx8HWkfrHwQeiZb/AJwULX8PqIiWTwfOjJaHRrl6n+j5X4Gro+X6cnIPYC1wM3AgMB64Ebgq2vYSYGJSrA6URMsdgA+As6PnXwIeTNr2aeBOQuvr16PfFydFn/M54LeN/Nn9O/pdcS1wA6HfenzS6yuBgdHy+4TfSUXAFcBiYAzw/Sjv/zsRb7T9fOCm6LXTgYX1xHAO8HS03Cf6OX8x7n9XGjWHeqwl57n7+2a2D3AN8EPgBDMb76HnugJY7+6PAZjZ1cCrwGDCL5EqM/tRtKvnCQn8eKBn0voPovd2As6P3ou7rzCzRSmE2gv4I/B9d/882sdmM1sOzHP3pdFxlgMveeg/fIuQ8HsRToI8nvBLAOBz4DtmNhu4g/BHhYhIc3BCsXkgsDtwkZkBrAK6WnhyHmFSAOBvhFYICLl2aq1cuxvwbq1jnE0dOdndV5vZOuApd58MYGYLgc/MbDfCDG+nOoN232hmnyStujDab8JdwMPR+rWE3xePR8eYGsXZWP/nUY+1mU0gTOYknOLuC83sUEJB3SnK/0uBOe7+JvBm9N669j0x+p0wpYGY9gc2Rp97hZmtpmZLieQAFdaS08ysNTDAw0kq/xW1ZzwDlBNmAGp7P3rsBgwE7nf3B6N1N0f7vJRQ6N5c61hfICTEDU0MdygwFhgCzEhaX7vHcNv2F9yroiSbaMvq6aH9BXefZ2ZnA7cBPwa+BXyCiEgGmVk3wrdt7xPy5qo68mMfwoxqOwB330aYDCB6z6/c/fXoeY33JqkzJ0ecmrlyMXAp8BHwH2BQAx8h+X17Av9Kej6H0OayKzvm4kqafsWS6YQ/QBLWmtl1wCOE2BPVc+3PtTOV1F+bPQvcYWbtCf8tNhJmvyWHqMdact32vj8Ad58BXA4Mr2f7doQk9jGwDDg5+UUzGxOtP97MOtRavz56OqKJsb4O3AL8T5T4UmJmxSSdCGNm/QhXRBlKOJnxwXreKiKSjrMI35A9QsiPh0YzxcD2/LgS2AIcm7R+HzPrTK1ca2ZtoomK2urLyXX5MzDfw2VVN6fwWRYSZri3HwJYR2h1yRh3r4omP0ZGf5g8A/za3adm8ji1jjkRmEj4ZvV0YJy7r2uu40nTqLCWfHBRVHQm9CP0xiX0MLN20fIRhP66tcDfga+b2W1mdkQ0m/AZIQF2BJ42sy9FJ4oMJxTjbwG/MbN2ZtaL0Nu8q5m1NrPLzOzMemK0aFxF+OXzhzpeq7197edfIfQDJuwJnOju64GL6niPiEgqdpgJNbPDgV8B57l7BWF2eAXwjJmdEJ2geFg0Q/0Q4eohJ5rZ0cBZ7r6BkGv/y8x+ZWbjCBMMiTa65PxXX05OSJ493g/obWa7AIcCHcxscPTaVsLJfokJluRj/Bn4hpl1j56PAf7i7pWEmqfOXGxmB0VxFdXzs2tFrZop+uPjeMK3lN0JLYbDgAHALkknRtaeFd9SR/w1g6qjXyT6b9EFeJEwK5/KHxySJWoFkXwwAHjfzP6X8DWZAT9Ier0SuMLM1hJmdy8CcPcXzOy/CF8nnkg4+e89ADP7GvAnwi+Ke4BL3N3N7JuESyG9D0wgfKU3mnBDhW8SWkxq3JjGzA4knJQzkPDV4JPAr6J4ngcOAkrMbAhhJqUv8E0zm0X17M3ZQEd3L6/12e+1cGOHSsLJQiIiKYvyyBmE4u8vQAXhBLiuwHHu/gaAu2+J8uNfCLnuKUJvNcDFhEmJ+wjf0CWuxHE3IfedT2hZO9/dP4kKx4OAIWb2QH052cxOJvQVn2Nms9x9CeFEyhsJRfV9wNeAfYC5wP2E2fXvWLha1Ejgy2b2jLs/bmZ7Ao+Z2WuEbz2viGbWTwH6Rp/vfcLNcnYzs5FRnKcR+safTPq5tSXk/lFAv6jwdkJv80nRPt8GXgKmENpbniXk9KcJJ0yOMrMT3T2x3wmEm/V8GhXifYGTzWweYSaa6P2179mwkHBy5olAe6CVmU0inMCYrUvQyk6Y/ltIPrNwG9i/untJlo53ibvrduMiIgUmakspcvecvEGXmV1CuELL/Oh5e8K3Ddclzs2R+KkVRPJdXW0WzXOg8JXp9GwcS0REsic6OXP/HC6q2wNXAiVJq3sBn6mozi1qBZG8FfW3nUn4au8Ud3+k1lILsAAAIABJREFUmQ85J/qKUkRECksbd78j7iDq4+6bLNxW/dboxPalhLaVKxt+p2SbWkFERERERDJArSAiIiIiIhmgwlpEREREJANUWIuIiIiIZEDBnLzYq1cvLykpSek9FRUV9O7du3kCaib5GDPkZ9yKOTsUM0ybNm2lu+fXDyFNTcnZkJ//XnKVfpaZo59l5uTDz7KhnF0whXVJSQlTp6Z2J9HS0tKU3xO3fIwZ8jNuxZwdihnMbEHGdpYnmpKzIT//veQq/SwzRz/LzMmHn2VDOVutICIiIiIiGaDCWkREREQkA1p0YV1WVhZ3CCnLx5ghP+NWzNmhmCUV+tlnjn6WmaOfZebk+88y6zeIMbNjgGuAU5Pud38ScBPQBXgA+Km7V0av9QJ+C6wF2gKXuvvW2vstLS31XO/JERGpi5lNc/fSuOOoi3K2iEhNDeXsrM5Ym1kfoBNwYNK6gcCJwMnAj4FzgB8lve0R4B53vxSYDfw6awGLiLRgytkiIqnJamHt7iuAp2qt3h34nrtPc/f7gD8BRwKY2UHAMHd/M9p2InCRmXXMVswiIi2VcraISGqy3mPt7lW1nr+S+AoxsgRYGC0fBcxL2nZJtPj/2bvz+Ljqev/jr0/2fWnTpvtC9401lE1o2WRT8CK44AYqFRfUnxcvCopeUPAqKoKiRAW8qIgCVxAQBIGyQ1OgULrQ0n1P0+x7mu/vj+8kmaRJmjQzmTPJ+/l4nMeZOefMOZ9Mp5/55Jvv+X4j8yfT++6Do4+Gm26KyOlERIaaoOTs6jeqKTmmhDWXrRnoqUREoiaI41gfC/wk9Hg8UNZlfzUwruuLSktLKSrqyN1Lliw5eAf4mhp44w048siBxCsi0m/FxcUUFxe3PS2IZSwDNCg5u7W+lZrXa0hIHdb33ItIjPQ1ZweqsDaz6UCpc25FaJMD6rsclgoccCPMqFGj+j+g+MiRfr13bz8jFREZmPBC0sziMgkNZs62ZPMXaB7cG+5FRKDvOTswv/qbWSLwBeBbYZu3A/lhxxiQCeyIyEXbCuuyrg0sIiLSm8HO2W2FdWtz60GOFBGJncAU1sBVwM3OucawbU8Cs8KeTwGagMiM0aTCWkTkUA1qzk5I9l9XarEWkSAb9MI61IIBYGHbrgHeANLN7DAz+6yZTXfOlQA7zWx+6NCzgducc00RCaYg1EVGXUFERLoVlJytriAiEg8GtY+1mWUDnww9vdTMfon/U+INXQ5d45y7M/T4o8B1ZrYJP57qtRELaMQIvy4vh9ZWSAhSA76ISGwFKWdbkgprEQm+QZ95MVoOeRavvDyorPTdQdoKbRGRQRTkmRejpb85u2FrA69MeoWU8SmcuO3EKEYmItK7wMy8GEht3UH27IltHCIi0iN1BRGReKDCevx4v96+vffjREQkZnTzoojEAxXWEyf69datsY1DRER6pBZrEYkHKqwnTPBrFdYiIoHVXli3qLAWkeBSYd3WYr1tW2zjEBGRHrWNCqIJYkQkyFRYqyuIiEjgtRXW7IehMpqViAw9KqzbCuvNm2Mbh4iI9MjMNJa1iASeCusZM/x6/XrYvz+2sYiISI90A6OIBJ0K66wsfwNjUxNs2hTraEREpAdthbX6WYtIUKmwBpg926/XrIltHCIi0qOE1NBY1o1qsRaRYFJhDSqsRUTiQEKG/8raX69ueyISTCqsoaOwXr06tnGIiEiPEjMSAWitU1cQEQkmFdYACxb49RtvxDYOERHpUXuLdZ1arEUkmFRYAxx9NCQkwFtvQUNDrKMREZFuqMVaRIJOhTX4kUHmzIGWFlixItbRiIhIN9pbrGvVYi0iwaTCus2xx/r1smWxjUNERLqlFmsRCToV1m2OO86vX3wxtnGIiEi31MdaRIJOhXWbU0/166efBqcxUkVEgkYt1iISdCqs28ycCePGwZ49sHJlrKMREZEu1GItIkE36IW1mZ1pZq+a2ZSwbelmdpuZ3WhmvzaznL7si3BgcPrp/vG//x2VS4iIxJsg5ezETN9irZsXRSSoBrWwNrNCIAtY2GXX7cALzrlrgIeAO/q4L7LaCut//jNqlxARiRdBy9lJuUkAtFS0ROqUIiIRNaiFtXNuNz7RtjOzccAlwD9Cm54CzjezKb3ti0qA553nx7N++mkoL4/KJURE4kXQcnZSfqiwLldhLSLBNOhdQZxzXe86WQzsdc7Vhfa3ABuBRQfZF3kFBbBokR/P+pFHonIJEZF4EqScrcJaRIIuCDcvjgfKumyrBsYdZF90fPjDfn3//VG7hIhIHItZzk7OTwZUWItIcAWhsHZAfZdtqUDzQfZ1UlpaSlFRUftSXFx8aNFceKHvDvLPf0Jp6aGdQ0TkIIqLi9vzFVAQ63j6IWY5u63Furn8gNOJiERVX3N20uCF1KPtQH6XbdnADnyS7mlfJ6NGjaKkpGTg0YwdC+ecA48+CvfcA9/4xsDPKSLSxZIlS1iyZAkAZrY3xuH0R8xydntXEN28KCKDrK85Owgt1s8A48wsHcDMkoFJwNKD7Iuez3/er3//e00WIyLSWcxytvpYi0jQxWIca2t7COCc2wU8AoTGuuMM4AHn3Pbe9kU1yPPOg8JCWLUKnnkmqpcSEQmyIOXspJwkSID9VftpbdbsiyISPIM9jnU2cEXo6aVmNir0+IvABWZ2LXAusCTsZb3ti47kZLjySv/4ppuifjkRkSAKWs62BCNldAoATbubInVaEZGIMTdEujoUFRW5iPSxblNRAZMmQXU1vPIKHHdc5M4tIhLGzJY754piHcdgOtScXXJMCTWv13D0q0eTszA6E/GKiPSmt5wdhD7WwZSXB1/5in989dXqay0iEgCp41IBaNzRGONIREQOpMK6N//1XzByJCxdCv/4x8GPFxGRqEoZG+oKslNdQUQkeFRY9yYvD773Pf/4qqugoSG28YiIDHMp40KF9Q4V1iISPCqsD+YLX4A5c2DdOrjhhlhHIyIyrLV3BdmuriAiEjwqrA8mJQV+9zswg//5H3jzzVhHJCIybKVNSQOgfkPXCR5FRGJPhXVfnHiiv5Fx/3749KehXgldRCQW0qenA9DwnrrmiUjwqLDuqxtvhJkz4e23Nc25iEiMpE5KxZKMxm2N7K/fH+twREQ6UWHdV1lZcN99vmvIb34Df/1rrCMSERl2EpISSJvqu4M0bFCrtYgEiwrr/jjySPjpT/3jSy+F5ctjGo6IyHDU1h2kbl1djCMREelMhXV/ffnLcNllvp/1Bz8IW7fGOiIRkWElY3YGALUra2MciYhIZyqs+8vMdwVZtAh27oRzzoG9e2MdlYjIsJF1ZBYANW/WxDgSEZHOVFgfipQUePBBmDsX3nkHzjwT9u2LdVQiIsNC1lEqrEUkmFRYH6oRI+Cpp2DGDD+29VlnQXl5rKMSERnyMmZnYClGw3sNtFS1xDocEZF2KqwHYuxYePppmDoVSko6uoeIiEjUJCQnkDk/E4CaN9RqLSLBocJ6oCZMgKVLYfZsP8b1+94HGzbEOioRkSEt98RcACpfqIxxJCIiHVRYR8LEifD881BU5IvqE0+El1+OdVQiIkNW7im+sK54riLGkYiIdFBhHSkFBb5byBlnwO7dsHgx3HVXrKMSERmSck8OtVi/WElrS2uMoxER8VRYR1J2Njz2GFx5JTQ1wWc/C1/5CjRodjARkUhKHZNK+sx0WmtbqS6pjnU4IiKACuvIS06GW2+F3/3OP/7Vr+C442D16lhHJiIypOSfmQ9A2SNlMY5ERMSLSGFtZudG4jxDyuc+5/tZT58Ob70Fxxzji23nYh2ZiAxzQyVnF3ywAICyf6iwFpFgSOrPwWZ2NvAdYAwdRXkCUAikDyQQM5sDfAVYD8wAip1zb5pZOvBjoBrIB652zlUN5FqD5phj4PXXfdeQP/wBLr8cnngCioshPz/W0YnIEDfUc3be4jwSsxKpfauWhs0NpE1Oi8ZlRET6rF+FNXAPcBPwBtB2t4gBl0QglnuAC5xz281sEvAEMAe4HXjcOXdf6EviDuDjEbje4MjOhrvv9rMzfvGLcP/98Mor8Mc/+nGvRUSiZ0jn7ITUBPLPymfvA3spfaCUid+YGI3LiIj0WX+7gmxwzv3MOfeMc25paHkWuC4CscwBskOP64FcMxuH/wL4R2j7U8D5ZjYlAtcbXJ/4hJ+h8fjjYds2OPVUuPZaaG6OdWQiMnQN+ZxdeEkhALvv2R2tS4iI9Fl/C+svmtlF3Ww/PwKx3A/8zsyygU8CVwKLgb3OuToA51wLsBGIz6beww6D556D734XzODGG/2EMuvXxzoyERmahnzOHnneSJLyk6h5s4aatzULo4jEVn8L65PwiXR/2NIK/DoCsXwJaACWATXOuQeA8UDXu1KqgXFdX1xaWkpRUVH7UlxcHIGQoiA5Ga6/Hp59FiZNgtdeg6OOgttvh/37Yx2diAyS4uLi9nwFFETpMkM+ZyekJjD6o6MB2PWHXYd0DhGRg+lrzjbXj1EqzKwc+AJQArRVgQZ8yjl3w6GHC2ZWAFyAT8xXhx4fCVzsnDsu7LjXgT87524Of31RUZErKSkZSAiDr6ICrrgC7rvPP1+4EH7zG19oi8iwYWbLnXNFUTjvsMjZVa9W8frxr5OUn8QJ204gMSMxIucVEelObzm7vy3W/wQecs5tcM5tDi2bgN8ONEjgj8ADzrnrgR8B9wI78XeVh8sGdkTgerGXlwf33utvaBw/3rdeFxXB//t/vugWERmYYZGzc47LIXthNi3lLez+o/pai0js9Lew/gvwbTM7JWxZBHx7IEGEWj6OcM61VZM/AKqAzcC40PBNmFkyMAlYOpDrBYoZfPjDfgKZr3/db7vlFj/+9W236eZGERmIYZOzJ3xtAgDbbt1Gf/4SKyISSf0trL8LfB34X+APYcsXBhjHPqDBzMaHbSsD3gIeAU4PbTsD30KyfYDXC57sbPj5z6GkBE45BcrK4KtfhXnz4P/+TxPLiMihGDY5e9RFo0gZl0LdO3WaMEZEYqa/hfVPgInOuSnOuamhZQpw4UCCcM61Ah8CrjOzJfgJDb4ZmlTgi8AFZnYtcC6wZCDXCryjjvI3Nv797zBzJqxbBxdeCCeeCP/6lwpsEemPYZOzE1ISmPhNP471xus24lqVK0Vk8PV3gpitzrnqbraPGWggzrkVdNOK4pwrBy4f6PnjihlccAGce66fpfH73/eTypx1Fpx0kn9++un+OBGRng2rnD3uinFs/clWalfUUvpgKaMvGh2LMERkGOtvi/WXum4ws4uBWyITjnSSnAxf/jJs3Ag33QQjRsCLL/pZHBctgmeeiXWEIhJswypnJ6YlMvnayQBsvHYjrU2tB3mFiEhk9bewHmFm3wUwswwz+x3wO/w4phItWVnwrW/Bpk3wwx9Cfj48/zycdhosXgxLh869nCISUcMuZ4/9/FjSZ6ZT/249236xLdbhiMgw06/C2jl3HvCkmd0GvAnMBI4AupvZSyItOxuuucYX2Ndf74frW7rUF9enneaLbRGRkOGYsxNSEphx6wwANl+/mcYdjTGOSESGk14L6y5DNJ1iZqcAKfghldYB3wMmA9dFP1Rpl5Pjp0XfuNH3t87N9d1CTjkFzjjDdxcRkWFHOdsbcdYICj5UwP6a/az/+vpYhyMiw0ivMy+a2SvAHKAC6KmzmgFjnXOpkQ+v7+Jy5sVIqajwQ/XdcgtUVfltZ54J//3fcMIJsY1NRA4qUjMvKmd3qN9Uz7L5y2itbWXufXMZ/RHdyCgikTGQmRevB2Y45yaHDdXUdZkCfCLSQUs/5OX5InrjRvjOd3yXkSef9EP0nX02PPechukTGR6Us0PSp6Qz/afTAXj3S+/SuEtdQkQk+notrJ1zjznn9hzsJM65+yMXkhyyESPghht8gX3NNf6mxyee8COIHHcc/O1vsH9/rKMUkShRzu5s7JKx5L8/n5ayFtZetlZjW4tI1PV3VBCJByNH+tFDNm6E667zz5ctg498BGbMgF/+EmprYx2liEhUmRmzfj+LpJFJ7Ht8H5t/uDnWIYnIEKfCeigrKPBdRLZsgV/9Cg47zBfbV17pC+y//EVdRERkSEubkMbcP80Fg03f28S+J/fFOiQRGcJUWA8HGRnwpS/Bu+/C/ffDMcfAzp3w8Y/7Yfo0VbqIDGEjzhrB5Osmg4NVH1tF3fq6WIckIkPUgAtrMzvKzP5oZv9pZlMGHpJETWIifPjD8Nprfqr0/Hx49lk/VfoRR8Ddd0NDQ6yjFJEoGq45e8p3pzDivBG07Gvh7fPepnlfc6xDEpEhKBIt1lcAZzjnfgqcaGbpETinRFNCAlx+OaxfDzfeCGPGwNtvw2WXwcSJ/sbHLVtiHaWIRMewzNmWaMy9dy6Zh2dS/249Ky9cqSnPRSTiIlFYfwU4JvT4L0ByBM4pg2HECPj2t/1MjnfdBUcfDXv3wk03wdSpvnX72WfVTURkaBm2OTspO4kFjywgZWwKlUsrWf2p1bj9ym8iEjn9LqzN7Bwz+0To8UL8mKnbAZxzrc65qgjHKNGWmgqXXgolJX7Wxo99zLdqP/ggnHqq7yZSXKyRRETikHJ2Z2kT01jwyAISsxMp/Wspaz+vYfhEJHL6VVib2Q+BvwH/AeCcew34tJldEoXYZLCZ+Ull7r0XNm/2Q/UVFvpuIl/4AkyYAN/8pm/hFpHAU87uXvbR2Sx4bAEJGQnsunsX665cR2+zEIuI9FV/W6xPBCYAz4dt+yfwo4hFJMEwblzHUH1/+hMcf7yfOv3mm2HaNLjoIl9wi0iQKWf3IO99ecx/aD6Wauy4fQfvfvFddQsRkQHrb2H9pHOuAgjPPhdFMB4JmpQUuOQSePllP5rIJz/pRxd54AHfReRjH4PVq2MdpYh0Tzm7FyPOGMH8B+eTkJbAzjt2svpTq2lt1g2NInLo+ltYbzCza4BZZvYxM3sA+DLwP5EPTQLn2GPhnnt8V5Arr4TkZLjvPpg3zxfc774b6whFpDPl7IMYee5IFvxzAYlZiey5dw/vfPgd9jfsj3VYIhKn+lVYO+f+AjwNJAKfBiqA9zvnfhWF2CSoxo2DW2+F996DL34RkpJ8d5E5c3wXkRde0EgiIgGgnN03+YvzOeLfR5A0Iomyf5Tx1plv0bS3KdZhiUgcsoHesGFmicBi59y/IxKQ2fHA+4D1wAtALfBjoBrIB67u7i72oqIiV1JSEokQpL82b4Yf/tBPMNMcmnTh6KPha1+Dj37UjzoiIj0ys+XOuaJBupZydg9qVtbw1tlv0bS9ibRpaRz+6OFkzMqIdVgiEjC95exeC2szexjIOsj5C4C9zrnTDj3E9utdBkxzzn0nbNtdwOPOufvM7GzgM865j3d9bRCT9LCzcyf8+td+2bvXbyss9NOpX3EFjB4d2/hEAipShbVy9sA17mjk7Q++Tc3rNSTlJzHvwXnkL86PdVgiEiADKax/ACQBa4FW4HLgDSA8G04FRjjnvjrAIBcB3wXOdKGgzGwcsBHId87VmVkSUAnMc85tCn99UJP0sNTQ4Ifsu+UWeOstv63tJsivfQ2OPDK28YkETAQLa+XsCNhfu59Vl6yi7OEyLMmY9rNpjP/KeMws1qGJSAD0lrOTDvLan+CL74rQiaY5567r5gIPDTxMfga8BNxmZtOBG4DJ+JaVOgDnXIuZbQQWAZsicE2JhrQ0Pz36pZf6mRt/8Qt4+GHfVeTuu2HRIr/vvPNg1KiYhioyxChnR0BiZiLzH5zPhm9tYOvNW1n/1fVUvVzFrN/OIjEzMdbhiUiA9XrzonOusi1BhxSG+ue1M7PFwMKBBGFms4CjgN86574C/Bt4HJgIlHU5vBoYN5DrySAx8zM3/v3vsG4dfP3rkJ0NS5f6wruwEE44wU9E8+yzvqVbRA6ZcnbkWKIx7SfTmHvfXBIyE9hz7x6WH7ec2jWagVZEetbf4fbuB1aY2W/M7Cdm9gjwJHDrAOOYB+xzzoX6DfBLwPBjr9Z3OTYVaO56gtLSUoqKitqX4uLiAYYkETVtGvz857Btm++DfdZZfri+V16BG27wBXh+PpxxBtx4I7z6KrS0xDpqkagpLi5uz1f4fs/RoJw9QKM/Mppjlh1DxuwM6t6pY/kxy9n+6+2aqVFkmOlrzu73qCChPnRfAg4HqoC/OecG9GdFM/sIUOycywvb9gbwV+Ay59zMsO3rgO855/4cfo4g99eTHlRXwzPPwNNP+6XrTI7Z2b7byGmn+WXBAkjo7++CIsEXzVFBlLMjo6W6hXVfWsfuP+4GYMQ5I5j1+1mkjtWoRyLDzSHfvNjHkycARznnlg/gHLOB1cAY59zu0LZl+FaQXwGjnHP1ZpYM1ACHOee2h58j3pK0dGPPHt8lpK3QXreu8/6RI33L9imnwMkn+0I7Uf0dJf4N8nB7ytkDsOeve3j3indpKW8haUQS038+ncJPFerGRpFhJGKFtZnd2c3mUcB+59yHDjG+tnM/BjzqnPuVmeUBS4FjgD8Cf3TOPWJm5wCfcs5d0vX18ZqkpRdbt3Zu0d66tfP+nBw46SR43/t8oX3ssf7GSZE4E63CWjk7Ohp3NLLms2sof6IcgLzT8ph5+0yNeS0yTESysH4ef5NK+ItOAp53zt0wwCALgF/gh4WahL8pZpWZ5eMnG9gEjAG+7Zyr6fr6eE7S0gfO+Zken3kGnn/eL5s2dT4mJQUWLuwotE88EfLyuj2dSJBEsbBWzo4S5xy779nN+m+sp6WsBUsxJn17EpOunkRiuv6SJjKURbKwnuWcW9tl22zgUufctwYW5sDEe5KWQ7Btm58+va3QXrmy81TqZnD44X7kkYUL/TJ7trqPSOBEsbBWzo6y5rJm3vuv99h15y4AUiemMvXGqRReUoglqHuIyFAU7T7WU4Flzrlo3dXeJ0MlScsAlJfDSy91FNrLlnVMsd4mOxuKiuCYY3wf7QULYM4cdSGRmBrkPtbK2VFQ8VwF6766jtoVfji+7KJspv10Gnmn6K9mIkPNQCaI6XqijXT+k2IiUAjcc+jhiURIfr6fdOa88/zz+nooKfFD9732ml9v2eK7kzzzTMfrEhNhxoyOQrttmTpVo5BIXFPOHjx5p+RRtLyIXffsYuO1G6kuqebNRW8y4rwRTPn+FHKKcmIdoogMgv52BSkG/kxHom4Fdjvn3o1CbP0y1Fo/JEp27fJF9ptv+uH93n7bjz7S2nrgsenpvjV77tzOy9SpkNSv30lFehXFriDK2TGwv3Y/W3+6lS0/3kJrrc8tIz8wksnfm6wCW2QIiGQf61Occ8912TYLf4f5+oGFOTBDOUlLlNXXw+rVvsheubKj4N6xo/vjU1Nh1qwDi+7p0/0NlCL9FMXCWjk7hppKm9h681a2/3I7rXW+wB5x9ggmXjWRvNPyNESfSJyKWFcQ4GzguS7bdgPP4Ke3FYk/6elw9NF+CVdR4QvuVas6L1u2wFtv+SVcUpLvUhJebM+Z47dlaBguiQnl7BhKGZXCtP+ZxsSrJrYX2Pse38e+x/eReUQmE78xkdEfG01CirqciQwVfWqxNrOvAD8Asns45Enn3NmRDKy/hkPrhwREdTWsWXNgwb1xY+dRScJNnAgzZ/qW7vD15MkapUQi3mKtnB1MTXub2PGbHWz/5Xaad/sbq1PGpjD+y+MZ89kxmsVRJE5EpCuImR0OfAa4tcuuJufczoGFOHDDMUlLwNTVwdq1HYX2O+/45++9d+DoJG1SUnwXku6K7oICP2SgDHnR6AqinB1crY2t7P7zbrb9bBu1K/0oIiRCwfkFjPvCOPLPzNdQfSIBFsk+1snOuQMqBDOb5JzbMoAYB2w4J2kJuJYWP5nNu+/6Qjt8vX17z6/Lz/cF9syZvjvJ9Ol+PWMG5OYOWvgSfVHsY62cHWDOOcqfKmfHHTsoe6gM1+K/j9OmpDH282Mp/FQhaZM0FKhI0BxyYW1mVwGPOedWhZ5fix+uqdNhwKnOucWRCffQKElLXKqp8aOSdC261671XU56UlDQUWSHF9zTp6vojkORKqyVs+NX485Gdt21i52/3UnDpob27bmLcin8ZCGjLhpFcl5yDCMUkTYDuXlxKhA+uv3RQCawi47hmww/ba2I9FdWFhx1lF/COQe7d/sCe/16X3y3LevXw969fnn55QPPOWpU9wX3jBmQo6G+hjjl7DiVOjaVyddMZtK3JlH+ZDk779xJ2cNlVC6tpHJpJeu+vI6RHxhJ4ScLGXHOCBLTdG+GSBD1tyvI0cBK51xTl+1z21pIYkWtHzJsOOeHAmwrssML7vXr/fCBPRk9uvuCe8YMPyulxEQUu4IoZ8exlqoWSh8sZfcfd1PxdEX7r0aJWYmM/MBICj5cwMhzRpKYqSJbZDBFcri9Bc6517ucfBYwDYhpkhYZNsxg/Hi/LF7ceV9rqy+6uxbcbes9e/zy0ksHnrewsPuie/p0Fd3xSzk7jiXlJDH20rGMvXQsjdsb2X3vbvb8eQ81b9Sw5y972POXPSSkJzDi7BEUXFhAwQcLSMrV5FUisdTX4fZOxv/5cAlQ3GX3COBW59ykyIfXd2r9EDmI1lZ/s2R3Rfd770FDQ8+vLSzsuU93Vtbg/QxDVBSG21POHsLqN9RT+mApex/YS9UrVe3bLdnIPSWXkeeOZOR5I0mfma5JaESiYMCjgpjZAuB/8f33yrvsbgL+5Jy7fqCBDoSStMgAtBXdXQvutqK7sbHn144Z033BraK7z6JQWCtnDxON2xsp/T9fZFc8V+EnrQ9Jm5bGyPN8kZ17Sq76ZYtESKTGsc4FLnbO/S6SwUWKkrRIlLS2wrZtBxbcbUV3U1PPrx07tnOXkrbCe9o0Fd1hojSOtXL2MNNc1sy+J/ZR9mgZ+x7fR8u+lvZ9CRkJ5J+eT/4ZfsmYk6FjPJ3UAAAgAElEQVTWbJFDFLFxrHu5wPuccy8M+EQDoCQtEgP79/dcdG/YcPCiu2tf7rZ1Zubg/QwBEK2bF3u5nnL2EOf2O6peraLs0TLKHi2jdkVtp/0pY1PIPz2fvNPzyD8jn7QJGi9bpK8GMo71w8DBmpUMmOecG33oIQ6ckrRIwOzfD1u3dr55sq9F97hxMGcOzJ0L8+Z1rEeMGLz4B1EEx7FWzpZuNWxroPzJcsr/XU75U+XtU6q3SZ+V7gvtU/PIOzmPlMKUGEUqEnwDGRVkLVABbKdjDNQDzg984NDDE5EhKTERpkzxyxlndN7XtegOL7w3bPAjm+zYAf/+d+fXFRZ2LrTb1iNHDtZPFXTK2dKttAlpjL1sLGMvG4tzjrpVdZQ/5QvtimcrqF9bT/3aenbcvgOA9Jnp5J2SR+7JueSekkva5DR1HRHpg4O1WBcA1c65Xu5cAjMb6Zwri3Rw/aHWD5EhYv9+2LwZVq+Gd96BVav8evVqqK3t/jWjRx9YbM+d6yfLiQMRbLFWzpZ+a21upbqkmvKnyql8rpLKlytprW3tdEzqxFRyT85tL7YzZmdgCSq0ZXgajD7Wn4/1DTJK0iJDXGsrbNniC+22YrvtcU1N968ZNQoOPxyOPNIvRxwBs2dDcrCmho5BH2vlbOlRa3MrNW/UUPFcBZXPV1L5fCUt5S2djknKSyL7uGxyT8gl5/gcshdmk5wfrP9XItEStcLazCYBXwa+4pwb8N1GZpYCLAO+5px71szSgR8D1UA+cLVzrqq71ypJiwxTzvluJeGt220Fd3X1gcenpPgW7bZCu22dl3fgsYNksArrSOfs0DkPKW8rZ8cP1+qoXVVL5XOV7cV2044D75HImJNBzvE57UvmvEwsUa3aMvREcuZFzCwROB+4HDgT2I3v0xcJ3wSmhD2/HXjcOXefmZ0N3AF8PELXEpGhwAwmTfLLOed0bG8ruFesgDff9MuKFX6IwDfe8Eu4qVPh2GNh4UK/HH30kBidJMo5G5S3hzxLMLLmZ5E1P4vxXxoPQMPWBqpeqWpfqpdXU7e6jrrVdey6axfgp17PXphNzsIcso7JIrsoW321Zcjrc2FtZlOBzwOXASOBzcAHgSeAEwYaiJmdCOwkNJmBmY0DLsG3rgA8BTxgZlOcc5sGej0RGeLCC+4PfrBje1UVvP12R6H95pv++caNfvnrX/1xCQm+Zbut2D7hBFiwwJ83DkQ7Z4euobw9TKVNTCNtYhqjL/aDy7Q2tlKzoqaj2H65ioZNDVQ8XUHF0x2/xyWNTCL7mGyyi7Lb16kTU1Vsy5Bx0MLazC7Gt3ScBqwHfgb8AbjGOfd46LCXBhKEmWUCFznnvmFm14U2Lwb2OufqAJxzLWa2EVgEbBrI9URkGMvJgZNO8kublhZ/c+Rrr/ll2TJ46y1fcL/9Ntx5pz9u7Fg4+2y48EK/Tur3H/2ibjBydug6ytvSLiE1gZyFOeQszIGv+m2NuxqpfrWa6pKOpXlvM+X/Kqf8Xx0TgiYXJJNdlN3eqp19TDapE1RsS3zqy7dCPX4K3JeAzzvn3gUws4Hf9djhauBHXbaNB7retV4NjOvuBKWlpRQVdXR3WbJkCUuWLIlgiCIyZCUl+dboBQvgc5/z2+rqfGv2smXw6qvw7LOwcyfcdZdfxoyByy6Dr33NDwPYT8XFxRQXF7c9LYjUj8Lg5GwYYN5Wzh76UsekknpBKgUX+I+3c47GrY1ULw8V2ss7iu19j+9j3+P72l+bNDKJrCOzyDoiy6+PzCJjdgYJyQmx+nFkmOtrzu7PlOaT8K0gU4CHgEXOuStD+8wd4l2QZnYWkOqcezj0fBNwKVCEn473uLBjXwf+7Jy7uet5dCOMiESVc74V+5FH4J57YO1avz09Ha64Ar7/fd8afgiiNKV5VHJ26PUDztvK2QK+2G7Y3EDN8pqOYnt5dafp2NtYipE5L7NTwZ15RCbJeRqNRAZXREcFMbMk4EPAJ4CVwO+B9zvnint9Yc/nexI4KmxTPlCD//PlJ5xzM8OOXQd8zzn3567nUZIWkUHjHLz0EvzkJ/DQQ37b5Mnw6KO+X3Y/RXNUkEjn7NA5B5y3lbOlJ845Grc1UvNmDTUraqh5s4baFbXUr6/v9vjUyamdW7ePyCJtSprG2ZaoieZwe9OBK4HLnXMZh3iOUUBq2KaXgW8ArwKrgFHOuXozS8Yn7sOcc9u7nkdJWkRi4o034PLLYflyP7LI2rX9Hid7EIfbG3DODp1nwHlbOVv6q6W6hdq3atuL7Zo3a6h9u5bWhtYDjk3ITCBzXiaZ8zPJXBBaz88kpTBFfbdlwCI63F4459x64GtmtnwA5ygNf25m+4FS59wWM3sEOB14BDgDeKC7olpEJGaOOgqef973s9640Y+jfeSRsY6qW5HI2aHzKG/LoEvKTiL3pFxyT8pt39ba0kr9uvoDWrebdjVR/Vo11a91Hss+uSD5gGI7c34mSTnBuxFZ4lNEPknOuf+NxHm68UXgx2Z2BDAG0J0tIhIsGzfCt77lJ6PJyvKt1gEXxZwNytsyiBKSEsick0nmnEwKP95xE3FzWTO1K2upXVlLzds17Y+b9zZT8WwFFc92Hso9dVJqp2I7a0HoZslU3Swp/RO4X9Gcc1PCHpfjb74REQmO2lp47DH405/g4Yd9n+v0dH9TY27uwV8/xChvS9Akj0wmb1EeeYs6ZlRt67tdu7KW2rdrO9ara2nc0kjjlkb2PdoxMgmJkDEjg4y5GWTOzWxfp89KJzEtMQY/lcSDwBXWIiKB45zvO/300/DUU/DEE344PvD9qT/2Mbj+epgyJaZhikjPzKx9YpuR54xs397a0krDew0dLduhort+fT11a+qoW1PH3gf3dpwoAdIPSz+g4M6YnUFipgru4U6FtYhIVzt3+psRly+HkhI/lvXu3Z2POe44uPhi+MQn/JjWIhKXEpISyJiVQcasDLioY/v++v3Ura2jblUdtatq29f16+vbl7KHOw/bnjYl7cCCe06G+nAPI/qXFpHhyTnYvt3PuLhqlV/aHpd1neMEGD0aTjvNL+9/vx9eT0SGrMT0RLKPzCb7yOxO21sbW6lb103B/W49DZsaaNjUwL7H9nV6TeqE1G4L7uQRGoN7qFFhLSJDV309bNoEGzb45b33Oh5v2OD3dycnB44+Go45BoqK/Hr6dNAwXSLDXkJqAlnzs8ian9Vpe2tzK/Xv1R9QcNetqaNxWyON2xo7TeUOkDIm5cCCe24GKaNSBvNHkghSYS0i8ampybc4b90K27Z1v96zp/dzjBwJc+f6Zc6cjsfjxqmIFpF+SUhOIHN2JpmzMxl14aj27W6/o35jNwX36jqadjXRtKuJiqc7j1KSXJDsC+15XQpujcMdeCqsRSRY6uth1y7fz7mn9fbtB/Z57k5iou+yMW0aHHZYxzJtmh8WLy/v4OcQERkASzQypmeQMT2DgvML2re7Vkfj1kZq36ntXHCvqqN5bzOVz1VS+Vxlp3Ml5Sd1buGel0nm3ExSxqngDgoV1iISXc754elKS2HvXr+UlvrCuLuiubLy4OcESEjwLcsTJ8KECR3r8MdjxkCS0pyIBI8lGGmT00ibnMbIcztGKXHO0bi9saPQfqej4G4pb6HqxSqqXqzqdK7EnMQDupNkzs0kdWKqpnYfZPrGEZH+qa+Hffv8UlZ2YMEcvm573NjY9/MnJ/uCeOzYjnX44zFjfEGtollEhiAzI21CGmkT0hjx/hHt251zNO1q6tyl5J1aat+ppWVfC1WvVFH1SueCOyEzoduCO21KmgruKNG3kshw5BzU1HQUyP1ZGhr6f720NBg1yi8FBX4pLOy+aB4xQv2bRUS6MDNSx6aSOjaV/NPz27c752gubT6gO0ntO7U072mmelk11cs6T+2ekJ5AxuyMTt1JMuZmkH5YOpao/DsQKqxF4tn+/b7rRG+FcHl599tbWg7tmsnJ/qa/ESP80rVgbnscvs7IiOzPLSIigC+4U0ankDI6hfzF+Z32Ne1tom71gUMDNu1oouaNGmreqOl0fEJaAhnzMsg6PIusI7LIPCKTrMOzNCxgP6iwFgmChobOBXDXYrin5xUVvvX5UGRm+sI4P7+jSO7LkpGhFmURkTiQUpBCyskp5J3c+Ubt5vJm6lbXHdDK3bi1kZrlNdQs71xwp05IJfPwTLKOCBXch2eSPiOdhKSEwfxx4oIKa5FIcQ6qqvpfHO/b1/N4yn2Rm9u/writmE5NjdzPLiIicSM5P5ncE3PJPTG30/bmimZq366lZkUNtW+F1m/Xto/DHT7xTUJaApnzM8k8PJPsY7LJLsom8/BMEtOG97TuKqxFumpq6ih6+1Mcl5dDa+uhXTM5+cCW496e5+d3LLqBT0REIiA5L5m8k/M6tXC7/Y769+o7Fds1K2po3NJIdUk11SXV7LpzFwCWZGQuyCS7KLt9yZyfSULK8GnZ1jeyDE3hN+f1pzjet88PDXeosrL6XhyHP8/MVPcKEREJHEs0MmZmkDEzAy7u2N5c0ewL7TdrqH7dF9h1q+ra+27v/O1O//oUI+uILLKLssk5Pofck3JJOyxtyI67rcJagq2lpfvW474Uy4d6c15i4oGtw30pjvPzfcuziIjIEJecl0zeKXnkndLRut1S0+IL7VBLdnVJNfVr69tHJtnx6x3+tYXJ5J6U275kHZU1ZFq1VVhL9DkHdXWHVhxXVR38/D3JyOh7t4rwfdnZaj0WERHpp6SsJPLel0fe+8KK7coWqt+opvq1aipfqqTqxSqadzez98G97H1wL+CH/8temE3+qfnkn5FP9sJsEpLjs9BWYS3955wfjaK0tO9LfyYICWfWuUW4r8Wxbs4TERGJuaTcJPIX57cPBeico35dPZUvVlL5oi+069bUUbm0ksqllWz6/iYSsxPJW5xH/hn55J+ZT8bsjLjpOqLCWryamo5ppffs6VwUd32+d2//u1mkpvqxj/tbHOfm+qmrRUREJO6ZdfTZHnvZWMCPt135QiUV/66g/Kly6tbUUfaPMsr+UQZA2tQ0Ci4ooOBDBeSclBPoYf5UWA9l+/f7InjXro6ieefOzo/b1v29YS87u2NikL4smZnR+RlFREQkrqUUpDDqQ6MY9aFRADRsa2gvsvf9ax8NGxvYdss2tt2yjaSRSRR8sIDRl4wm/7T8wM0UGZjC2swuBH4K5AB/BP7TOddiZgXAD4EKIBm42jnXHLtIA6KpCbZvh61bOy9btsCOHb5Y3rPHF9d9kZbWMaV0YWHvRXJBgT9eRIYt5WwRiZa0CWmM+cwYxnxmDK7VUfVqFXv/vpe9/7eX+nX17Lp7F7vu3kXqhFQKP1lI4WcKyZwdjAa8QBTWZjYJ+BBwETAP+DWwFbgZ+BvwLefcq2Z2BXA98O1YxTpoamvhvff8smlTR9HcVkDv2tW3GfdGjuwomHtb5+Tohj0R6RPlbBEZLJZg5J6QS+4JuRz2o8OoW1PHnvv2sPue3TRsaGDLj7aw5UdbyD8rn4nfmEj+mfkx7Y8diMIamAx81jnXAiw3s/nAqWb2AjDLOfdq6Lh/AGvN7AfOuQEMNhwQ+/bBu+92FNBty/r1sHt3769NSIBx42DSJJg4sWOZNAnGj/fFcmEhpKQMzs8iIsPJ8MzZIhJTZkbmnEymfn8qU743hcoXKtn1h13suXcP5U+UU/5EOZnzM5lywxQKLiiISYEdiMLaOfd8l03bgWzgNGBj2HHbQ29SEbB00AIcqH37YNUqeOcdv7Q93rWr59ekpMDUqTBtml+3FdBt63HjNOOeiMTEkM/ZIhJ4ZtY+S+S0H09jxx072H7bdmpX1vLOf7xDzkk5TL9lOjlFOYMaV1Ars2OBnwBLgLIu+6qBcYMeUV+VlcGyZVBS0rHesaP7YzMyYOZMmD7dF9Bty/TpvtU5MXFwYxcROTTxm7NFJO4lj0hm8rcnM/E/J7Ljjh1svn4zVS9W8frxrzPp6klMuW4KCamDM5JI4AprM5sOlDrnVpiZA+q7HJIKHHAjTGlpKUVFRe3PlyxZwpIlS6IaK+D7PT/zDDz9NDz/PGzceOAxGRkwZw7MnQvz5vll7lyYMkVDyYkMU8XFxRQXF7c9LYhlLAMRdzlbRIashJQEJlw5gTGfGcOm/97Etp9vY8uNW6h4poL5f59PyuhD7x7b15xtri83wA0SM0sEfgR8xznXaGbfBk51zr0/tN+AhtC2l8JfW1RU5EpKSqIfZHMzLF0KDz4I//qX7xMdLj0djj4aiorg2GP9esYMFdAi0iMzW+6cKzr4kcESFzlbRIatyhcrWfXxVTRubSRtShpHLj2StEkDH9Wst5wdtBbrq4CbnXNt0/Q9CVwRtn8K0AQMbjZ2Dl5+GX7/e/j7332f6TY5ObB4MZx2GixaBPPnq++ziAwXwczZIiJA7km5HP3a0aw8fyXVy6pZceYKjnrxKFIKojewQ2AqQDO7Bp98083sMGAx8Byw08zmO+dWAmcDtznnmgYlqOZm+MMf4NZb4e23O7bPng0f/jCcf75vnVYhLSLDTCBztohIF6ljUjn8X4fz5uI3qV1Ry9rL1jL/4flRGzEkEBWhmX0HuKHL5jXOuTvN7KPAdWa2CcgCro16QM7BfffBtdfChg1+26hR8LnPwac+5ftHi4gMU4HL2SIivUjOS2bBPxZQcngJZY+UsftPuxnzyTFRuVag+lgPRMT66+3cCUuWwCOP+OezZsF3vwsXX6wxoUUkKuK1j/VAqI+1iAy2nXftZO1n15I2JY2FaxeSkHJo97/1lrN1R124FSv8DYePPAK5uVBcDCtXwic+oaJaREREJI6N+fQYMuZk0LCpgb1/3xuVa6iwbrNypb/5cPt2OOkkP4HL5Zer/7SIiIjIEGCJxrgv+mH1d93dyyR9A6DCGqC8HM49Fyor4YIL4Kmn/AQtIiIiIjJkjP7YaADKny5nf/3+iJ9fhTXAVVfB1q2wcCH85S+QNvAxDkVEREQkWFJGpZB1ZBau0VH1UlXEz6/Ces0auOsu34f6D39QUS0iIiIyhOWclANA9evVET+3CutbbvHD6112mR+fWkRERESGrMx5mQDUra6L+LmHd2Hd0gL33+8ff/nLsY1FRERERKIuc26osF6jwjqyXnkFyspg5kw/FbmIiIiIDGmpE1IBaNzRGPFzD+/C+uWX/frUUyFKU1uKiIiISHAkFyYD0Ly7mUhPlDi8C+tly/z6uONiG4eIiIiIDIqkrCQSsxJpbWilpbIlouce3oX1unV+PW9ebOMQERERkUGTNNJPANhSocI6cjZt8uspU2IZhYiIiIgMosTMRAD210R2kpjhW1hXVUFFBaSnw6hRsY5GRERERAZJYpYK68gqK/PrggLduCgiIiIyjKiwjrTycr/Oz49tHCIiIiIyqFRYR5oKaxEREZFhKSHdl8Ctda2RPW9EzxZPKir8Oi8vtnGIiIiIyKBKSPElsGvWONaRUV/v15mZsY1DRERERAaVpfj761qb1WIdGQ0Nfp2WFts4RERERGRQJSSHWqyb1GIdGW2FdXp6bOMQERERkUHV1mId6a4gSRE9W5SYWTrwY6AayAeuds5VDeikarEWEYmKqORsEZEIsuRQV5Cm4dkV5HbgBefcNcBDwB0DPmNaGowbByNGDPhUIiLSSeRztohIBCXnJ5MyNoWEtMiWwuZcZJvAI83MxgEbgXznXJ2ZJQGVwDzn3Ka244qKilxJSUmMohQROXRmttw5VxTrOCJBOVtEhrrecnY8tFgvBvY65+oAnHMt+KS9KJZBiYhItxajnC0iw1Q89LEeD5R12VYNjAvfUFpaSlFRxy8PS5YsYcmSJdGPTkTkEBQXF1NcXNz2tCCWsUSYcraIDDl9zdnxUFg7oL7LtlSgOXzDqFGj0J8VRSRehBeSZrY3xuFEknK2iAw5fc3Z8dAVZDv+rvJw2cCOGMQiIiK9U84WkWErHgrrZ4BxoeGbMLNkYBKwdKAnDmvSjxvxGDPEZ9yKeXAo5iEnajkb9N5Hkt7LyNF7GTnx/l4GvrB2zu0CHgFOD206A3jAObd9oOeOx3+8eIwZ4jNuxTw4FPPQEs2cDXrvI0nvZeTovYyceH8v46GPNcAXgR+b2RHAGEB3uIiIBJdytogMS4Efx7qvzKwU2NzPlxUA8XbTUDzGDPEZt2IeHIoZJjvnRkXwfIF3iDkb4vPzElR6LyNH72XkxMN72WPOHjKFtYiIiIhILAW+j7WIiIiISDxQYS0iIiIiEgEqrEVEREREImBYFtZmlm5mt5nZjWb2azPLiXVMPTGzM83sVTObErYtsPGb2YVmttHMyszsF2aWFNpeYGZ3mNn/mNnPQmPbBoKZLTKztWZWbma3hG0P7PvcxsxSzGyFmS0OPY+HmH9gZi60rAhti4e4jzezq8zsQ6HPc+BjHkr0fkdWd98t0j89fd/JoenpuzjeDMvCGrgdeME5dw3wEHBHjOPplpkVAlnAwi67Ahm/mU0CPgRcBPw/4PPA10O7/wbc6Zy7GngXuD4mQXZhZrnAScDxwCeBL5vZGaHdgXyfu/gmMCXseaBjNrNM/B3fZ4aWi0K7gh73ZcAHnHM3O+f+7pzbS8BjHoL0fkdIL98t0kcH+b6TfjrId3FcGXajgpjZOGAjkO+cqwv9hlkJzHPObYppcN0wswRgPzDVObcpyPGb2cnAy865ltDzHwPzgBuAB51z40LbxwNrgULnXG2s4g3Fkg40uNB/BDNbDlwVii+Q73MbMzsRmA1cB1yK/4Ul6DG3ffH8xjnXENoW2M80+FYU4LvAmWGfk0DHPNTo/Y68rt8tMQ4n7vT0feecOy+2kcWnnr6LnXPPxDay/huOLdaLgb3OuTqA0H+KjcCiWAbVE+dca5dNiwlo/M6559uSTMh2YAtwGj7GtuPaZmArGsTwuuWcqw/7j5wJrAaeJcDvM7THepFz7s6wzYsJdswGfAz4KbDbzD4e2rWYAMcN/Az/ubjNzB43s5MIfsxDzWL0fkdUN98t0g+9fN/JIejluzjuDMfCejxQ1mVbNTAuBrEciniK/1jgN8RBzKHWh38CyUAGwY/5auBHXbYFOmbnHY/vCnIb8CczO4cAx21ms4CjgN86574C/Bt4HJhIQGMeogL7GREJafu+kwHo5rs47gzHwtoB9V22pQLNMYjlUMRF/GY2HSh1zq0gPmJeD9wNnA78hADHbGZnASXOuT1ddgU25nDOuXLn3HeAm/B9EoMc9zxgn3PurdDzXwJGsGMeivR+S2B1+b6Tgen6XRx3huMdrNuB/C7bsoEdMYjlUAQ+fjNLBL4AfCu0aTswM2y/AZkEKGbn3E7gTjNrxd8Q+DzBfZ+vAo7ybyPg43wI32UhqDF351b8n/qC/JlOIixPOufqzWwdkEhwYx6KgvwZkWGsm+87GYBuvovjznBssX4GGBfqKE9o2LdJwNKYRtV38RD/VcDNzrnG0PMngVlh+6cATUDJIMfVFyX4L/Egv8+XAEeGLTvwd6TfRXBj7k4r8DrBfq/fAnJDoyi0acG/50GNeSgK8mdEhreu33cSGW3fxXFn2BXWzrldwCP4PzMAnAE8EHZDXaBYR7OkQfDjN7NrgDeAdDM7zMw+C1QAO81sfuiws4HbnHNNsYqzjZllmNlRYZvOBX4R5PfZOVfqnNvWtuDv7C91zm0hoDGDH9nBzC4xs4TQ5/o/gWsD/l6vwff3uwjAzPKAFOBPBDTmoSjIn5F41fW7Rfqvu++7ULcQ6aeevotjFc9ADLvh9gDMLB/4MbAJGAN82zlXE9OgumFm2fjxHG/Hj/v8S+dcaVDjN7Pv4IfWC7fGOTfHzCbjh4XbhB8/9doud1THhJktBB4F1gCvAK855/4W2hfI97krM9sEXOqcezbIMZvZPHyRWgW8ANwSKlwD/V6bWQE+wZfgW0l/65xbFeSYhyK935HT03dLbKOKL71938UinnjX23dxvBmWhbWIiIiISKQNu64gIiIiIiLRoMJaRERERCQCVFiLiIiIiESACmsRERERkQhQYS0iIiIiEgEqrEVEREREIkCFtYiIiIhIBKiwFpFAMbNRZvawmW02s2/FOh4REemd8nYHTRAjIoFiZmcBzwA5+Fm4Zjrn9sU2KhER6YnydgcV1iISWGb2HHCac64l1rGIiMjBDfe8ra4gIhJIZjYR+NdwTc4iIvFGeVuFtUSBmf2Xmf0q1nG0iUQ8ZnaGmb1sZosjFNaQYmY3R/h8BnwE+FEkzysinSlfD5yZTTGz35rZHRE6X6KZXW5mm/tw7Hlm9uqh/qyRzN3K254K6zhjZnPM7EEz+52ZrTAzZ2afjHVcXTwG3BvrIMJEIp6twBERiCWQzOxcM9trZjvN7AQzmxj6nDkz+6WZ5YWOW2hmK83sp2aWHtpWCOyKcEiXAL93zrWEzi8Sd5SvD0k85us9oXVqhM6XAFQAk/pwbAkwp6edZjbPzO4Lffa+Y2ZjwvZFOncrb6PCOh79Hfhv59zngaOA/23bYWaHm9miWAQVfm3n3Ern3AuxiKM7/Ymnp/fQObcW2BvJuCL57zXQcznnHgNuB1qccy8757YClwNNwKvOuYrQca/hE/k3nXP1oZd/EHhkQD9AGDP7LvB94CkzWwPMitS5RQaZ8nU/xWO+ds7VATsidS3nXDOwvI/H7sYX4T3tfwcoDj292TkXXkhHLHcrb3dQYR1HzGwEMBOoBXDOtQLXAc7McvFJ22IQV8yuHUl9+DlaB/FasTjXfcAEMzsGwDlXBvwb+EDYtfKAHaHPXpvZzrk1A7x2O+fcDc65Gc65IufcbOfcc5E6t8hgUb6OrgDm60iPBNGf+A927P7Qumu/54jlbuXtDiqs40ho6Jq1wONmdnxo22bgNXy/pqnAEjO7DMDMrjez74b+DAjh8jMAAAfOSURBVNTej8rMiszsZ2b2dTOrM7O1ZvYNM/uQmS0zs/8ws/8zsyoz+0LY67o9X5drX25md4f3NQv9OfQ2M7vJzJ4xs4Wh7b1eL+z1ZmZfMLN3zOxMM3vLzHab2aV9uMaM8HgOcs0D3sNu5If+tFttZpd3eU+/H3pvHgwl4r6+Zwdcy8wuMrPXzOxSM1tnZj/q5Xzd/dt3G09vQi0b7wAXh20eBZxnoW4fwPnAw2FxpgN1Yc9nm9mPzOxmM1tqvkvJR82sxMwuMLN7zazCzL5mZseEfsYyMzv3YPGJxBPl6+GTr7vT08/Y23XMbKqZ3Wm+9fcnfblON9d9zMz+ZGYnHeQ45e5occ5piaMFmAdsxv+GeicwOmzfJmBx6PF8YEPo8Sj8b9Mj8b9xbwGKQvt+BbwUepwC7AN+iR+L8uPA9t7O1/Xa+D5mdwB3h7an4Yu1vNDzi4BSILu363Xzc48JXfOa0LE/ABqAyQe5Rtd4er1m+HvYTQybgJ8CGaHXbQltzwH+Enbci/ik2Kf3rIdrjcD/ee8WYCFwRm/n6/Jv3208ffx8fQd4L/R4GnAT0AhcGPZ5sbDjzwdODj3OAJYBGaHnz4XehxygErgeSAc+Gnr+Efwv9/9J6DOoRctQWlC+Hhb5OrT/+318H3v7t14GzAl7jevj56zt3/ME4Hvd7F8cuk5S2Dbl7igtarGOM863Ki7AJ5pPA++YWXc3aawDPmVmycDJoW1Z+P/IE4G2/rGr8f9hcM41AdXA/c65Kvx/tDEHOV/X+BrpfDPEeUC9C/XRBR7E/4e88CDX63retnP+OXTsD/D9f888yDU6xdOfa/bgH873p1sGjA37GUeGWpS+jn9PoY/vWQ8/7z58AnvIOfeac+6pfpyvp3j64j7gMPPdQT4M/BZ4AviImWUB1S6UeUNOAl4KPT4f/+XV1gryH8ANofe5Anja+X7ZJUCOc+6vzv95/A36928gEheUr4dHvu5Gjz9jL9c5A8h1zrXFU9LPax4LfMI59999PF65O0qSYh2A9J2ZJQITnXObgK+a2e/xd1AXA8eFH+ucazSz8cDVwB/aTuGc22NmJcDp+N+op+KLqfaXhj1uIdRdqKfz9RBq+DlmAMlhcbWa2SZgQm/XOxjnXIOZbQDygIJ+XOOQr9lFCx3/fyYBG51zt3Q9qC/vmZl9Cl/Atpnl/J+MXXis/fg36DGeg3HOrTOzN/DdQXKdcxvM7D58K9LFwKNhcbd9Ntr6703GtzC1nass/NRhj/fTWSvqlvb/27t3ELmqOI7j318W1yIBNUgguypaCAEfpFDSiOIiikbYxheIYERIiJEtxEJYEMRGyBY+ClGEFGJEFEIQDIg2PltFVCKk0IiViBosspC/xf+MnhnmzsyOd9xM5veBCzv3dc49d/efM+cVu8A4Xnc93yzE61pjOQ54NzeSXyLGtQRcK2lbRJxpOOdceQbH7glygUyXi6jGv0bEV2RX267eEyUtAQcj4vnIFR5qTwA3K5d9+hp4YVjCQ+43yI9kC2i9DJHIsYf/1TxwcsJpjOIXuschI2nPBsrsOLC72vrOLt/A/frmZ/TH4W1gH3Cqyt8csEJ2m3bsIceLdvwM3CJpa5XuwHF+Zhcwx+tusxSvG59xQDp/lGvmGc9hsrX/xYbjv8e/k84duyfIFevp86SkherzInCi/HyWnKyxi1za6RJJF0u6oxzfIelycqzfG2Q3zueUrsViCz3f0iVpyP1601Z1j2PAGeChcq9LS3qdls+m9JrsLOcslmtPjJBGnZ9hadbP0atf64XIVqitwPuS7lJOkhn2Dv5JC9gZEd9X23qVxFz186D71fnumx9J85JeK8M8BnmH7IJ+FyAi/gQ+AD6L7tVA7iSHiXR0yvuoci3sQ0BnLdPed9D7nqd6hQKzBo7XzEy8HrUc+6ZDxtJ58ssXZO8EkhZGiN0iW48fAe6X9EC59ilJq2Rr/UfV+Y7dkzSJgdveJrOREyIC+A14EzhCdiVdVo4/C5wGlsmunR/IFoK9ZEvHETJIHScnUpwt9/uL/EO7G1gnJ8hcCTxTjj826H49aa8An5Dj1q4rx24iWzrXyLGG15f9jek1PH+U61fIoQk3VMea0ri6zs+wNOsy7El7mZzE90opi851+8rx24BvyrtZK+U8apkt93nW+0p6R4HFsm/QO+26V0N+riD/I4O3Rvhde6nn84PA7T37Dve57lbgW+BXYLXsu4f8XXsZWKAsOQYcIP/hfb08672b/TfmzVtbG47XsxSvrwI+JVvddw95xkHvei/ZU/glOZH8JHCwPHvf2E32iqyTLdXbgY/JSv3jwH5yPe81chhK5xrH7gluKoVnM0LSDuBARDxXPm8hW1H2R8TqpmZuCEkBXBM5ZtHGJOnpiBhrKScz+/84XlvNsXs6eCjI7DlEDgvYXu1bopqUdj6qup7c7TQmSXOSHgbe2+y8mNlIHK/NsXvKuGI9e14lx2J9J+kn4EPgdER8sbnZalYmc6yUj48ql32zjdsGHIuIU0PPNLPzgeO1gWP3VPFQEDMzMzOzFrjF2szMzMysBa5Ym5mZmZm1wBVrMzMzM7MWuGJtZmZmZtYCV6zNzMzMzFrgirWZmZmZWQtcsTYzMzMza4Er1mZmZmZmLXDF2szMzMysBX8DjcV5AL9mNqYAAAAASUVORK5CYII=\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/exomars-2016-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/exomars-2016-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/exomars-2016-mars.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Compare to results from https://doi.org/10.1016/j.actaastro.2018.05.029" ] } ], "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 }