{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 36 - Beagle-2 - 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 Beagle 2 aeroshell\n", "https://en.wikipedia.org/wiki/Beagle_2" ] }, { "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.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('Beagle-2', 68.46, 70.0, 0.00, 0.67, 0.0, 0.417, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,5.40,0.0,-15.8,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (30*60.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXxcdb3/8dcnTRu609J0X9JSuligLGERlB0ERUUEBVzAhShXcbmioCBXBK4I6lVBwYi48VMQQarIcgEpyy0iKW2RtpTSpi2kpU03mha65vP743smmUwn+yxnJu/n4/F9zMyZM+d8TibzySff+Z7vMXdHRERERES6pyTfAYiIiIiIFAMV1iIiIiIiGaDCWkREREQkA1RYi4iIiIhkgAprEREREZEMUGEtIiIiIpIBpfkOIFOGDRvmFRUVWdl2fX095eXlWdl2HPW044Wed8w97Xgh3sc8d+7c9e4ez+CypKs5O87vY1fpmAqDjqkw5OKY2srZRVNYV1RUUFNTk5VtV1ZWZm3bcdTTjhd63jH3tOOFeB+zma3Mdwy51tWcHef3sat0TIVBx1QYcnFMbeVsDQUREREREckAFdYiIiIiIhmgwroDqqqq8h1CTvW044Wed8w97XihZx5zMSrG91HHVBh0TIUh38dk7p7bHZqdClwHfNTdV0TLzgZ+CAwC7gS+5u67o+eGAdcDm4HewOXuvit1u5WVlV5s44REpGcws7nuXpnvONJRzhYRaamtnJ3THmszGwEMAI5MWjYeOAs4B/gq8FngK0kvuwe4w90vB14BvpuzgEVEejDlbBGRzslpYe3ua4FZKYsnAJ9297nu/jvgZ8CJAGZ2NDDV3Z+L1v0bcKmZ9c9VzCIi7WlsbGTXrr06ZQtenHK2O+zZE5qISFzlfIy1uzemPH468RVipA5YFd0/CahNWrcuuhvLr0xFpOd55JFH6NWrF+9///vzHUpWxCVnz5kDpaVw3HHd3ZKISPbE8eTFI4DbovtjgA0pzzcAo1NfVF9fT2VlZVOrrq7OcpgiItC7d2+ATvdYV1dXN+UrYFjmI8uZnObsHJ8WJCICdDxnx+oCMWY2Gah39wXRIgfeTlmtDNjrL1h5eXnRTXIuIvGXKKx37tzZqddVVVU1nb1uZuszHlgO5DJnm3U1ShGR7utozo5Nj7WZ9QI+B1yRtLgOGJK0jgH9gdW5jU5EJL0+ffoAnS+sC12+crZ6rEUkzmJTWAOXAT9w9x1Jyx4FpiY9rgB2AuqaFpFYSBTWxXjyYjtymrMTPdYqrEUkznJeWEc9GACWtOxbwDygr5lNMrNPm9lkd68B1pjZgdGqpwM3u3vP6hoSkdjq6lCQQhGXnK2hICJSCHI6xtrMBgIfjx5eZGa3EL5KvDZl1Zfd/Y7o/keBq81sBWE+1StzEauISEcU81CQOOZs9ViLSJzltLB29wbg1qglXBe11l6zEvhMlkMTEemSYh4KEqecraEgIlII4jTGWkSk4BT7UJC40FAQESkEKqxFRLqhmIeCxJF6rEUkzlRYi4h0Q1cvECOdo6EgIlIIVFiLiHSDeqxzQ0NBRKQQqLAWEekGFda5pR5rEYkzFdYiIt3Qq1cvAPbs2UNjY2OeoyleGgoiIoVAhbWISDeYWVFPuRcXGgoiIoVAhbWISDepsM4d9ViLSJypsBYR6SbNZZ19GgoiIoVAhbWISDfpBMbs01AQESkEKqxFRLpJQ0FyRz3WIhJnKqxFRLpJQ0GyT0NBRKQQqLAWEekmDQXJPg0FEZFCoMJaRKSbdFnz3FGPtYjEmQprEZFuUo919mkoiIgUAhXWIiLdpMI6+zQUREQKgQprEZFu0lCQ3FGPtYjEmQprEZFuUo919mkoiIgUAhXWIiLdpMI6+zQUREQKgQprEZFu0lCQ3FGPtYjEmQprEZFuUo919mkoiIgUgtJc79DMTgWuAz7q7iuiZX2BG4EGYAhwubtvae85EZE4KObCOi45W0NBRKQQ5LTH2sxGAAOAI1Oe+jnwjLt/C5gF/KKDz4mI5F2xDgWJY85Wj7WIxFlOC2t3X0tItE3MbDRwAfC3aNFjwAfMrKKt53ISsIhIBxRrj3WccraGgohIIcj5UBB3b7SW3+mdAKx397ei53ebWS1wPLCrjedW5DJuEZHWJHqsi62whvjkbA0FEZFCkPPCOo0xwIaUZQ3AaGB3G8+JiMRCose62IaCtCKvOVs91iISZ3EorB14O2VZGaHng3aea1JfX09lZWXT46qqKqqqqjIYpohIel0ZClJdXU11dXXi4bDMR5U1ecnZGgoiIvnU0Zwdh8K6jnDmeLKBwGpCAm/tuRbKy8upqanJSoAiIm3pysmLyYWkma3PSmDZkZecraEgIpJPHc3ZcZjH+glgdDRFE2bWGxgPPNnOcyIisVCsJy+2Iq85Wz3WIhJnOS+srfksGANw9zeAB4CTo+WnAPe6e11bz+UwZBGRNhVzYR2XnK2hICJSCHI6FMTMBgIfjx5eZGa3uHs9cAlwo5nNBEYCyQPt2npORCTvinge69jkbA0FEZFCkNPC2t0bgFujlrx8E3BxK69p9TkRkTgo1h7rOOZs9ViLSJzFYYy1iEhBK9bCOk40FERECoEKaxGRbirWoSBxoqEgIlIIVFiLiHSTeqxzRz3WIhJnKqxFRLpJhXX2aSiIiBQCFdYiIt2koSDZp6EgIlIIVFiLiHSTeqyzL1FYNzbmNw4RkbaosBYR6aZEj7UK6+yJ/ndBXwqISJypsBYR6aZEj7WGgmRP9L8L+t9FROJMhbWISDdpKEj2JXqs9SMWkThTYS0i0k0aCpJ9GgoiIoVAhbWISDdpKEj2JYaC7NqlKfdEJL5UWIuIdJOGgmSfGZSWhvv6/0VE4kqFtYhIN2koSG5onLWIxJ0KaxGRbtJQkNxQYS0icafCWkSkmzQUJDd0AqOIxJ0KaxGRbtJQkNzQXNYiEncZKazN7L2Z2I6ISCEqtKEghZqzy8rC7fbt+Y1DRKQ1pZ1Z2cxOB64CRtJclJcAI4C+mQ1NRKQwxLXHuthy9oAB4Xbr1vzGISLSmk4V1sDvge8B84DGaJkBF2QyKBGRQpIorHft2oW7Y2Z5jqhJUeXsgQPDbUNDfuMQEWlNZwvr5e7+o9SFZvZyhuIRESk4JSUllJaWsnv3bnbv3t1UaMdAUeVsFdYiEnedHWN9iZmdk2b5BzIRjIhIoYrpcJCiytmDBoVbFdYiEledLayPBW43sz1JrRG4NQuxiYgUjJiewFhUOVs91iISd50dCvJdoAqoAfZEywz4RHcDMbPpwBeBV4EDgGp3n29mfYEbgQZgCHC5u2/p7v5ERDIppnNZF1XOThTWW/QXQERiqrOF9UPALHffkbzQzH6ZgVh+D3zQ3evMbDzwCDAd+DnwsLvfHZ3h/gvg/AzsT0QkY2I6FKSocnZiKIgKaxGJq84W1ncB3zSzfyQtM+Bs4MvdjGU6EPVH8DYw2MxGE85e/0K0/DHgXjOrcPcV3dyfiEjGxHQoSFHl7KFDw+3Gjd3dkohIdnS2sP424Su/iwCPlhlhjtTuJuk/E8YCngF8HLgUOAFY7+5vAbj7bjOrBY4HVnRzfyIiGRPToSBFlbOHDQu369d3d0siItnR2ZMXbwLGuXuFu0+MWgWh96O7/gPYDjwPbHX3e4ExwIaU9RqA0akvrq+vp7KysqlVV1dnICQRkY4pLQ39FB3tsa6urm7KV8CwLIVVVDlbhbWI5EtHc3Zne6xfc/d052OP7OR20ukL/JGQmH8c9XI44SvGZGXAXn+5ysvLqampyUAYIiKdl7gojLu3s2ZQVVVFVVVV4rXZKhWLKmcnCuv6+s4HKyLSHR3N2Z3tsf6P1AVmdi7w405uJ507gXvd/bvADYSEvYZwVnmygcDqDOxPRCRjOltY50hR5ezy8nCrHmsRiavOFtZDzezbAGbWz8xuB24nfBXYZWY2DJjp7pujRdcBW4CVwOho+ibMrDcwHniyO/sTEcm0GF3GPFlR5ezkHuvGxrbXFRHJh04V1u7+PuBRM7sZmA9MAWYC6a7s1Rkbge1mNiZp2QbgReAB4ORo2SmEHpK6bu5PRCQr4tRjXWw5u18/2G8/2LkT1q3LxBZFRDKrzTHWZnZcK0+tBJYCPwAmAGcBX+1qEO7eaGZnAVeb2VxgBPB1d99iZpcAN5rZTMK4wKqu7kdEJFviMBSkJ+TsigrYsAFWrICRmRgpLiKSQe2dvHgjYa7SzUC6L97uIEzdNIpuJGkAd18AfC7N8k3Axd3ZtohItsWhsKYH5OwJE2Du3FBYH310tvYiItI17RXW3wVq3L3NL93MrLtfK4qIFLSYFNZFn7MrKsLtypV5DUNEJK02C2t3f7AjG3H3P2cmHBGRwhSHkxd7Qs6eODHcvvpqfuMQEUmns7OCiIhIG+J08mIxmjYt3C5enN84RETSUWEtIpIBMRkKUvSmTw+3ixeDftQiEjcqrEVEMkCFdW6MHg2DBsHGjboCo4jET7cLazM71MzuNLOvmVlF90MSESk8hVJYF3rONmvutV60KL+xiIikykSP9eeBU9z9h8AxiStuiYj0JLt37wagtLS9yZbyruBzdvJwEBGROMlEYf1F4PDo/l1A7wxsU0SkoOzatQuA3r1jnwILPmcffHC4feGF/MYhIpKq04W1mZ1hZh+L7h8JHJC4XK27N7r7lgzHKCISezt37gSgT58+eY6kpWLM2UceGW7/9a/8xiEikqpThbWZXQ/cA3wIwN3/BXzSzC7IQmwiIgUjjj3WxZqzDz0UevWCl16CrVvzHY2ISLPO9lgfA4wFnk5a9hBwQ8YiEhEpQHEsrCnSnN2vXxgO0tio4SAiEi+dLawfdffNQPJp7wV7aVwRkUyJ6VCQos3ZRx0Vbp97Lr9xiIgk62xhvdzMvgVMNbPzzOxe4AvA9zMfmohI4Yhpj3XR5myNsxaROOpUYe3udwH/AHoBnwQ2A6e5+8+yEJuISMGIY2FdzDk70WM9Z46uwCgi8dHpCVfd/Z/APxOPzayXmZ3s7o9nNDIRkQKSGAoSp8IaijdnT5sGI0bA6tWwcCEceGC+IxIRaaewNrO/AgPa2cYwYD1Q0ElaRKSrdu3axa5duygpKaGsrCxvcfSknF1SAqefDr/9LTz0kAprEYmH9nqsX4zWWQI0AhcD84CapHUmAkOzEp2ISAHYGs35NnDgwKZLm+dJj8rZZ5wRCusHH4Svfz3f0YiItF9Y3wRYdFY5Zra/u1+dupKZzcpGcCIihaChoQEIhXWe9aicfdppoef6//4Ptm2D/v3zHZGI9HRtnrzo7m8mEnRkhJn1Sl7HzE4AjsxCbCIiBSFRWA8Y0N4ojOzqaTl7yBA47DDYtSsU1yIi+dbZkxf/DCwws2eABmA68B5grx4REZGeInkoSMwUfc4++WSoqYFHHgk92CIi+dTZ6fYeBU4jnPgylTB10znu/r0sxCYiUhDi0mOdqifk7Pe+N9z+/e/5jUNEBLo23d5q4KrEYzMrMbPD3X1uRiMTESkQMe6xLvqcfcwxYUjIkiXw8sthGj4RkXzpVGFtZnekWVwO7AHOykRAZnY08C7gVeAZYBtwI+FrzCHA5e6+JRP7EhHJhE2bNgEwePDgPEfSUk/I2aWlcNZZ8Otfw113wXe+k609iYi0r7OXND8AWAWsTGplQEZ6PszsU8CZ7v4Dd7/f3dcDPweecfdvAbOAX2RiXyIimVJfXw9AeXl5niPZS4/I2eedF27vuktXYRSR/OrsUJDPuvuS5AVmNg24qLuBmNnxwMeAU5OWjQYuAL4QLXoMuNfMKtx9RXf3KSKSCTEurHtEzj7pJCgvD8NB5s+HQw/N1p5ERNrW2ZMXl6RZvAP4bAZi+RGwGLjZzB42s2OBE4D17v5WtP/dQC1wfAb2JyKSEevWrQNg+PDheY6kpZ6Ss0tL4dxzw/1f6DtNEcmjThXWZlZrZsuT2kpCYv1Ld4Iws6nAocAv3f2LhEvtPgyMAzakrN4AjE7dRn19PZWVlU2turq6OyGJiHRYV3qsq6urm/IV4TLjGdeTcvYXvxguFvOrX8FLL3VrUyIie+lozu7sUJBHgT8AiVFsjcBad3+lS1E2mwFsdPcXo8e3AP8V7eftlHXLgF2pGygvL6empiZ1sYhI1nWlsK6qqqKqqgoAM1uflcB6UM6ePh0uvjj0WF94Ifzzn9C7d8Y2LyI9XEdzdmdPXrzT3We7+5NRezps3yZ3I1YIBX5Tke/ubwNLgV6Es8qTDQRWd3N/IiIZ88YbbwDxGwpCD8vZN94I48fDCy/A94pmpm4RKSSdLaxPT7NsLXBPN+N4ERhsZiOSlu0mJOPRZtYXwMx6A+OBJ7u5PxGRjNixYwd1dXWUlJQwbty4fIeTqkfl7EGDwrR7AN/9LjypvxQikmMdKqzN7Itmthm43Mz2JDfCeLq13QnC3V8GHgLOifa3L9AH+H/AA8DJ0aqnAPe6e1139icikikrV67E3Rk/fjy9YzL2oCfn7JNOgq9/HfbsCSc0vvZarvYsItLBMdbufouZPQVcCPw05emd7r4mA7F8EviJmfUh9HCc7+67zewS4EYzmwmMBKoysC8RkYxYvnw5ABMnTsxzJM16es7+7/+GefPgscfgfe8LPddDUgeoiIhkQYdPXnT3F83sCnff6yQUMxvv7qu6E0h0YYGPpVm+Cbi4O9sWEcmW2tpaACZNmpTnSFrqyTm7tDRcLOZd74J//xvOPBMefRT69ctnVCLSE7RZWJvZZcCD7r4oWvQNM+uVuhpwImH+UhGRHuWlaG63KVOm5DkS5exk++0H//u/cOyxMGdOKK7//GcYOjTfkYlIMWuvx3oisG/S48OA/sAbNE/fZISv+0REepz58+cDcGg8LvennJ1k3LhQXB9/PDzxBBx5ZOjJDtPQiohkXpuFtbt/IWXR9cBL7r4zeaGZvSPTgYmIxN2ePXtYsGABEI/CWjl7b9Omwb/+BWedFS53/s53wtVXwze/GYaMiIhkUmen2zsoTYKeCuyfuZBERArD0qVL2bZtG2PHjmXYsKxcPLG7lLOBCRPCcJAvfxl27w6F9bHHhvmuRUQyqUP/r5vZuwlfH55qZrUpTw8lnHX+twzHJiISa0888QQAxxxzTJ4jaUk5e299+8KPfwzvfz9cdFHoxa6shM9/Hq67TmOvRSQzOvpF2Gbgd4Txe8emPLcTuD2TQYmIFILHHnsMgFNPPTXPkexFObsVJ58MCxfCNdfAT34Ct94Kd98N3/gGfPGL0L9/viMUkUJm7t7+WoCZDQbOdfdYJuTKykqvqanJdxgi0kPs2rWL4cOHs3nzZmpra6moqOjytsxsrrtn9JQ65ez2LVwIl14aTmwEGD4cvvUtuPhiTc0nIq1rK2d3eIy1u7/ZWoI2s3d1NTgRkUL02GOPsXnzZqZPn96tojpblLPbN2MGPP54mDnkyCNh3Tr4yldg/Hi46ipYk4nL6IhIj9LePNZ/BQa0sw0DZgDDMxWUiEjc/eEPfwDgggsuyHMkzZSzO88MTj0VTjkFHngArr0Wnn8err8ebrwRzj8/nPR42GH5jlRECkF7Y6yXEMbq1dE8B2oqA87MZFAiInG2efNm/vKXvwBw3nnn5TmaFpSzu8gsnNh45plhBpEf/Qj+8hf43e9CO/xwqKoKhfbAgfmOVkTiqr3C+vtAg7vvaGslM+tRZ5eLSM92++23s23bNk466SQmT56c73CSKWd3k1mYiu/YY2H5crj5Zvjtb2HuXPjc5+A//xMuuAAuvDDMiV3S2UlrRaSotZkS3H19ewk68qEMxSMiEms7d+7k5ptvBuCrX/1qnqNpSTk7syZNgv/5H6irgzvvhOOOg23b4Je/hHe9CyZODLOJvPACdHAeABEpct36X9vMxpvZ94GfZCgeEZFYq66uZtWqVUybNo33vve9+Q6nU5Szu6ZvX/jYx+DJJ2HxYrjssnC59FWr4KabwjCRadPC1RznzIE9e/IdsYjkS6cLazPrZWYfMrMHgWXAxwhj+kREitqWLVu45pprAPje975HSQGMA1DOzqxp00IxvWIFPPMMfOELYZq+V16BG24IQ0hGjYJPfQruuw+2bs13xCKSSx3+q2BmE83seuA14C5gMvB+YBxwbnbCExGJjyuuuIL169dz7LHH8sEPfjDf4bRJOTu7SkpCEX3LLWGoyOOPh9lDJk6E+nr4zW/gwx+G/faD006D738famrUmy1S7NotrM3sXDP7X2Ap8GHgR8BY4O/u/rAHc7Icp4hIXs2ePZtbb72V3r17c+utt2Jm+Q4pLeXs3CsthZNOCpdMX7YMXnoJvve9cHLjrl3w6KNwxRVwxBFQXh4K7p//HBYt0thskWLTkUuav024BO4c4LPu/gqAmSkdiEiPsHbtWj72sY8BcNVVV3HQQQflOaI2KWfnkVm48MyMGaGYXrcO/vGP0KP92GNhCMl994UGoUf73e9uboceGgp1ESlM7X583f0B4AEzGw9cbGYVwCygd2IdMzPv6LXRRUQKyK5du/joRz/K6tWrefe73803v/nNfIfUJuXseBk+HM47LzQIU/g9/nhoTz8Nq1fD/feHBjBgQOjpfve7wywkRx4ZTp4UkcJgnc2tZlYKnEU4AeYl4FfAae5enfnwOq6ystJramryGYKIFJnGxkYuvPBC7rzzTkaNGsULL7zAyJEjM74fM5vr7pUZ3zDK2XHmHgrtp59ubkuXtlynT58whCTRo33ssTB4cH7iFZGgrZzd6cI6ZcOTgUuBi929X5c3lAFK0iKSSe7OV7/6VX7yk5/Qv39//vGPf3DkkUdmZV/ZLKxT9qOcHXNr1oTZRp5+Gp56Cl58seU4bDOYPj0U24k2cyaUleUvZpGepq2c3a2RXO7+KvBlM5vbne2IiMRJY2MjX/rSl/jZz35Gnz59uP/++7NWVOeScnb8jRoF554bGsDmzWFu7KeeCsX288+Hkx4XLQpXhATo3RsOPrhlsT1tWlguIrnVrR7rTDOzPsDzwJfdfbaZ9QVuBBqAIcDl7r4l3WvV+yEimbBz504uuugi/vjHP1JWVsaf/vQnPvCBD2R1n7nqsc6GruZt5eyu2b4dFiwIU/c9/3xoixfvPbtIWVk4gfKQQ0KP9iGHhOJ7333zE7dIMclaj3UWfB2oSHr8c+Bhd7/bzE4HfgGcn4/ARKT4rV69mnPOOYdnn32WgQMHMmvWLE488cR8hxV3yts5tM8+cNRRoSU0NITLqicK7ZqaMHb7hRdCS1ZR0Vxsz5wZiu9JkzQTiUimxOajZGbHAGuATdHj0cAFwBeiVR4D7jWzCndfkZcgRaRoPfPMM5x77rm88cYbjB07llmzZnHYYYflO6xYU96Oh4ED4fjjQ0vYsiWMz54/P/Rwz58f5tdesSK0xCwkEE6QnDoV3vGOlu2AAzScRKSzYlFYm1l/4Bx3/08zuzpafAKw3t3fAnD33WZWCxwPrMhLoCJSdHbs2MF3vvMdbrzxRhobGzn++OP505/+xPDhw/MdWqwpb8fboEHwrneFlrB7d7j0eqLQfvHFMIxk5Ur4979DS1ZaClOmtCy2p00LyzQFoEh6sSisgcuBG1KWjQE2pCxrAEbnJCIRKXrz5s3joosu4sUXX6SkpITLL7+ca6+9lt7qpusI5e0CU1raXCCfnzQ4Z+vWUGAnTopMtNra5vvJzMKl26dNCzOUJG6nT4ehQ3N7TCJxk/fC2szeA9S4+7qUp5xwBbFkZcCudNupr6+nsrJ5HHlVVRVVVVWZDFVEisSmTZv49re/za233kpjYyP7778/v/3tbzn22GNzFkN1dTXV1U1TSQ/L2Y4zIBN5Wzk7PgYMaJ5NJNlbb8GSJS2L7cWL4dVXwxju5cvhwQdbvqa8vLnITi64x46FkpLcHZNIpnU0Z+d9VhAzexQ4NGnREGAr8CPgY+4+JWndpcB/ufsfUrejM8xFpD27d+/m17/+NVdeeSX19fX06tWLSy+9lGuvvZYBAwbkLa5CmxUkE3lbObtw7dwJy5aFInvxYnj55ebbbdvSv6Zfv1BoJxfbU6fC/vtrWIkUnrjPCnIBoUcj4VngP4HngK+bWV93f9vMegPjgSfzEKOIFLDGxkb+/Oc/8+1vf5tXXnkFgOOOO45bbrmFgw46KM/RFSTl7R6sT5/m4jhZYyPU1e1dcC9eDOvWpZ+lBEJv9uTJe7f99w+96SKFJO+FtbvXJz82sz1AvbuvMrMHgJOBB4BTgHvdvS4PYYpIAWpsbGTWrFlce+21zJs3D4D999+fa6+9lvPOOw8zy3OEhUl5W9IpKYFx40I77bSWz23c2LJne/HiMMxkxQp4/fXQZs/ee5sjR7ZedGtObomjvBfW7bgEuNHMZgIjAQ3AE5F2bd++nTvvvJObbrqpqYd6zJgxXH311XzqU5/SyYnZpbwtexk6FI45JrRku3fDqlVh3HZqW7YM3ngjtGeeSb/NiRObW0VF8/0JEzTERPIjdoW1u1ck3d8EXJy/aESkkNTV1XH77bdz22238cYbbwAwYcIEvva1r/HZz36WvvpLmxXK29JVpaXhAjWTJu3dy71nTxhaklpsJ+5v3Bja3Lnptz1yZPqie+LE0Kuu/68lG2JXWIuIdEZjYyOPP/44t912G7NmzWLPnj0AzJw5k2984xt85CMfoVSXlRMpOL16wfjxoZ10Usvn3GHt2jAl4IoV4Ta5rVrV3Nv97LN7b7ukJIztTi68J0wI+xo3Ljyn/8OlK/TXRkQK0rJly7jzzjv5/e9/z7JlywDo1asXH/7wh7nkkks46aSTNIZapEiZhR7pkSPhne/c+/lEb3dysZ1cgNfVheJ71Sp4spVTa8vLQ5GdKLZT748apUvBy970KyEiBWPDhg386U9/4ve//z3PJnVDjR07lqqqKj7zmc8werSuRSLS0yX3didf6j1hx45QVCcX26tWwWuvhdu6OqivDy3dTCaJfYwevXfBPWZM6PEeMyYU/r16ZfVQJWZUWItIrK1Zs4b777+f++67jyeeeKJpqEf//v05++yz+cQnPsFJJ51EL/31EpEOKm87RcgAACAASURBVCuDAw4ILZ09e8JQk0SxnSi4k++vXdv8eM6c9NspKQk928nFdvL9xK2GnRQPFdYiEjvLli1rKqafffZZEheyKi0t5T3veQ+f+MQnOOuss+jfv3+eIxWRYpTojR49Go4+Ov06O3Y0DylJ7ul+/fXm23Xrwv26OvjXv1rf35AhexfbqQX4kCFhCIzEmwprEcm7t99+m9mzZ/Pwww/z0EMPsXTp0qbn9tlnH97znvdw9tlnc+aZZzJ06NA8RioiEpSVNc9o0pqdO2HNmpbFdqLQTr6/aVNo//5369vaZ5/0vd2pQ0807ju/9OMXkZxzdxYtWsTjjz/OQw89xOzZs9m+fXvT8/vuuy+nn346Z599NmeccUZeLzcuItJVffqE2UYmTGh9ncZGWL++9cI7cbtlS5huMDpXO62SklBcjx3behs9OvxTINmhwlpEss7dWbJkCbNnz+aJJ55g9uzZrFu3rsU6hx9+OGeccQann346Rx11lKbIE5EeoaQEhg8P7dBDW1+voSF9b3fq0JPVq0Nra+hJeXnrhXeiB1wj7bpGf7lEJOP27NnDwoULmTNnDk899RSzZ89mzZo1LdYZNWoUJ554IqeffjqnnXYaI0aMyFO0IiLxN3AgTJsWWmsSQ08ShXa6tnp184wn8+a1vq19922/+B48WOO+U6mwFpFu27x5M8899xxz5sxhzpw5PPfcczQ0NLRYZ/jw4Zx44omccMIJnHjiiUyZMkXzTIuIZFBHhp7s2RN6tlsrvBNt8+bQXnqp9W3179960T1+fJh1pV+/zB9nnKmwFpFO2b17N4sWLaKmpqapmF64cGHTzB0JFRUVHHPMMRxzzDGceOKJTJ8+XYW0iEie9eoVpgAcNQqOOCL9Ou6wYUP7xfe2bbBkSWitGTcOpkyBqVPDbeL+hAnFOce3CmsRaVVjYyOvvPIKNTU1PP/889TU1DBv3jzefvvtFuv16dOHww8/nGOOOYZ3vvOdvPOd79SFWkRECpQZDBsW2iGHpF/HPZxQ2VrRXVsLy5c3z/X9+OMtX9+nD+y/f3PBnVx4l5cX7hATFdYiAoRx0UuXLmXBggXU1NRQU1PD3Llz9xrSATBx4kSOOOIIjjjiCI455hgOO+ww9tlnnzxELSIi+WAWxlgPHgwzZqRfZ/fucHXLJUvglVdCS9yvq4PFi0NLNWRI2OaMGXDggc1t2LCsHlJGqLAW6YEaGhp48cUXmT9/PgsWLGDBggX8+9//3qsnGsLlwisrK6msrOSII47g8MMPZ7/99stD1CIiUkhKS2Hy5NDe976Wz23dCkuXtiy2E/c3bYJnngkt2YgRzcX2QQfB4YeHx3365O6Y2qPCWqSIuTurVq1iwYIFLYroZa1MhDpu3DgOOeQQDj300KYietSoUTmOWkREit2AAWF6wdQpBt3DzCYLF4b20kuhLVwYLiO/di384x/N6/fp01xkH344HHZYeJyvubpVWIsUAXfn9ddfZ+HChS3aokWL2Lp1617r9+nThxkzZjBz5kxmzpzJIYccwsEHH6yrGoqISF6ZNV9O/tRTm5e7h8vGJwrtBQtg7tzQyz13bmgJvXvDzJlw/PFwwgnwrneF6QNzQYW1SAFxd+rq6pqK5uQCesuWLWlfM3z4cA4++OCmAnrmzJlMmzaN3r175zh6ERGRrjFrnkoweVjJli1hPu65c+GFF8LtkiVQUxPaD38YXnv44fDhD8O554aTJrNFhbVIDO3cuZNXX32VJUuWsGTJEl5++WWWLFnC4sWLefPNN9O+ZtiwYcyYMWOvNqwQzvYQERHpgkGDQs/08cc3L2togOeeg9mz4cknw/1Eof3Nb8Jxx8G3vw2nnJL5eFRYi+SJu1NfX9+icE7cr62tZc+ePWlfN3To0LQF9PDhw3N8BCIiIvEzcGAomhOF81tvwaOPwj33wF//Ck89FYaZfPGL8NOfZnZqPxXWIlm2bds2li1b1qIHOlFAb968Oe1rzIxJkyYxdepUpk6dyrRp05puR4wYoQutiIiIdFC/fvDBD4a2ZQvccgtcc024Pf30vWcs6Q4V1iIZ0NDQ0FQ8L126lFdffbWprV69utXXDR48eK/ieerUqUyePFnzQouIiGTYoEFhrHV1NaxcGWYgySQV1iIdtGXLlrSF89KlS1m7dm2rr+vTpw+TJk1i8uTJHHDAAS0KaPU+i4iIZN9bb8Ejj8Cdd8L990NjY5hf+5xzMruf2BTWZnY28ENgEHAn8DV3321mw4Drgc1Ab+Byd9+Vv0ilWLk7b7zxBsuXL2f58uXU1taybNmypkK6vr6+1deWlZWx//77NxXPkydPbmrjxo2jV69eOTwSkexTzhaRONu5M5ys+OSToT31FCSugda7N3z60/D972d+Gr5YFNZmNh44CzgHmAHcCrwG/AC4B7jC3Z8zs88D3wW+ma9YpbBt3bqV2tpaamtrmwroRBFdW1ub9sqDCX379mX//fffq3A+4IADGDNmDCUlJTk8EpH8Uc4WkTjZvj3MbZ2Ycu+FF+DFF0Nxneyoo+BDH4JPfCLMk50NsSisgQnAp919NzDXzA4ETjSzZ4Cp7v5ctN7fgCVmdp27b8tXsBJfe/bsoa6ubq+iOXF/3bp1bb5+2LBhTJw4kUmTJjW1RCE9atQoFc8igXK2iOSUO2zcGC6Dvnhxy1ZbG4Z2pJo+PUzDd9xx4TZbxXSyWBTW7v50yqI6YCBwElCbtF5dNB61EngyZwFKrGzatKmpWE7teV65ciW7drX+rXOfPn1aFM6p9wcNGpTDIxEpTMrZIpJpjY2wbh3U1YUrLK5YEVptbfP9hob0ry0pgWnTWl7W/NBDw4mKuRaLwjqNI4CbgCpgQ8pzDcBe/3PU19dTWVnZ9LiqqoqqqqpsxihZsn37dlasWNFiiEZyId3aBVISRo4c2aLHObmAHj16tHqdJRaqq6uprq5OPCz0q/goZ4tIq7ZuDQXz6tXhNtGSH69ZA7t3t72dAQPCCYfTpze3adPggAOgrCy7x9DRnB27wtrMJgP17r7AzBxIHfRaBuzVJVleXk5NTU0uQpRuSh6uka5wXtPO3Df9+/dn4sSJew3ZmDRpEhUVFfTr1y9HRyLSdcmFpJmtz3M4XaacLdIzbd8Oa9fCG2+0f7t1a8e2ud9+YbjG+PEwcSJUVDS3iRNhyJDMXsylMzqas2NVWJtZL+BzwBXRojpgStLzBvQHWp8YWPLO3dmwYUOrhfOqVavaHK5RWlrK+PHjmwrnRBGdaOXl5ZqiTiQGlLNFiod7KIDr65vbunWhME5XLLfz5XELZWWhYB4zpmVLXjZ6NBTD5RtiVVgDlwE/cPcd0eNHgc8nPV8B7ATUzZFn27Zta1E0JxfOtbW1bG3n39ORI0emLZwnTZrEmDFjKC2N26+miKShnC0SU42NsHlzy0I5XVu/vvn+jh3tbzehtBRGjAht5Mjm2+T7idvBg/PX05xrsalezOxbhOTb18wmAScATwFrzOxAd38JOB242d13tr4lyYRdu3bx2muvtVo4tze7xqBBg1otnCsqKujbt2+OjkREskE5WyS33noLNmwIhfCGDc33WyuS16+HPXs6t4++faG8PLRhw/YukJNvhw4NJw1KS7EorM3sKuDalMUvu/sdZvZR4GozWwEMAK7MdXzFyN1Zu3Zt2qK5traW1157jT1tfCL79OlDRUXFXsM0EoX0kCFDNFxDpEgpZ4t0nXsYRpEojtMVy+me27698/saPDgUyIliObWlPte/f+aPt6eJRWHt7tcB17Xy3ErgM7mNqDhs2bKl1cK5vYuhmBljx45NWzRPnDhRs2uI9GDK2SLB7t2waVPbBXHq/Y0b25/9Ip2yslAI77dfaIn7rRXOw4ZBnz6ZP2ZpWywKa+manTt3snLlyrSF8/Lly9m4cWObrx86dGirhfOECRMoy/bcNSIiIjGxY0fHe48TbdOmru1r4MC9C+TU+6mP+/XrOeOUC5kK6xhrbGxkzZo1rRbOdXV1uHurr99nn33SFs2JNnjw4BwejYiISPYlZrfo7FCLbV24NqhZmAKurYI49bn99lNPcjFTYZ1n7k59fT2vvPLKXu3VV19lRxun6JaUlLQ6Ld2kSZMYMWKExjmLiEjBSsxs0V7PcerjnV04Xba0dO8hFu0Vy/vuC716Zf64pXCpsM6RhoYGli5dmraAbutKguXl5a0WzuPGjaN37945PAoREZGuSRTJ69aF1pFiedOm8LrO6tev/eI49fHAgRpqId2nwjqD3J3Vq1ezaNEiFi5cyKJFi5qK57auJjh48GCmTp3KlClTWrTJkyczcODAHB6BiIhIx7iH4ROJQjm5JS4ukrqss9O/QfPMFh0pjhNNM7pKvqiw7gJ3p66urql4TtwuWrSo1d7nsrIyJk+ezJQpU/YqoocNG6YhGyIikneNjaGneM2a5qvspSuSE4VyG5NLpbXvvjB8ePOsFe0Vy0OHhiEaIoVCv67tePPNN1mwYAHz589nwYIFTUV0Q0ND2vWHDh3KjBkzeMc73sE73vEOpk2bxpQpUxg3bhy9NBBLRETyYNeuUCSvWdN2W7u2c1PB9e0bLhYyfHhzwZy4n7qsvFwn7UnxU2EdSfRCz5s3j/nz5ze15cuXp11/v/32ayqgk2+HDx+u3mcREcmZrVvhtdfg9dfDbaLV1YViefXq0AvdxiRSLQwdCqNGhZZcNKcrlnVBEZGWemxhvWbNGubMmcPzzz9PTU0N8+fPZ8OGDXutV1ZWxoEHHsghhxzCzJkzOfDAA5kxYwbl5eUqoEVEJOsaG2H27PTF82uvhRMC22MWiuTRo5uL5nRt5MhwIRIR6ZoeW1j/9Kc/5YYbbmixbOjQoRx66KEccsghTW3q1KmaeUNERPLGDN773nABk3TKymDsWBg3rmUbO7a5YB4+XGOVRXKhx37MjjvuOObNm8cRRxxBZWUlhx12GGPHjlUvtIiIxIoZnHUWlJTsXTyPGxdO9NOfLpF46LGF9RlnnMEZZ5yR7zBERETadddd+Y5ARDqiJN8BiIiIiIgUAxXWIiIiIiIZoMJaRERERCQDVFiLiIiIiGSACusOqK6uzncIOdXTjhd63jH3tOOFnnnMxagY30cdU2HQMRWGfB+TCusOyPeblGs97Xih5x1zTzte6JnHXIyK8X3UMRUGHVNhyPcxqbAWEREREckAc/d8x5ARZlYPrMzS5ocB67O07TjqaccLPe+Ye9rxQryPeYK7l+c7iFzqRs6O8/vYVTqmwqBjKgy5OKZWc3bRFNYiIiIiIvmkoSAiIiIiIhmgwlpEREREJANUWLfDzE41s+fMrCLfsWSbmZ1tZrVmtsHMfmJmpfmOKdvM7HgzW2Jmm8zsx/mOJ1fMrI+ZLTCzE/IdSy6Y2XVm5lFbkO94pG1m1tfMbjaz/zazW81sUCvrnWFmt5jZz83snFzH2VEdza1mdoCZ7Un6XT0017F2Vkc+W4XyPgFEf+89pb03zXr7mtmWpHU+lI94W5Ouduno5ypaN3bvWSvH1OG6JVefr6IvnLrDzEYAA4Aj8x1LtpnZeOAs4BxgBnAr8Brwg3zGlU1mNhg4FjgaOAa438wecPfH8htZTnwdqMh3ELlgZv0JJ7OcGi3K1knOkjk/Bx5297vN7HTgF8D5ySuY2RTgeqDS3RvN7Ckze9ndX8pDvK3qZG79LHAG0Ajscvd5uYqzKzry2SqU9wnAzKYDzwKXAbuixbcD/0iz+sXAhUBD9PiJrAfYQW3ULu1+rqLXx+49S3dMXahbcvL50smL7TCzEmAPMNHdV+Q5nKwxs3cDz7r77ujxjcAMd39ffiPLHjPrC2z36ENgZnOBy9w9NgkyG8zsGGAacDVwkbvPzm9E2WVmX4nu3ubu2/MajLTLzEYDtcAQd38r6oF6k5CPViStdxuw0d2/FT2+Epjs7p/KQ9it6mhujQqHXwGfd/fXcx9p53Xks1Uo7xOEQs3dVyU9ngDc6O4fTVmvDHgA+A93X5rjMDsktXbp6Ocqem0s37M0x9ThuiWXny8NBWmHuzfmO4ZccPenE7+ckTpgVWvrFwN3fzupqO4PLAZm5zWoLIuO8xx3vyPfseSCmRlwHvBDYK2Z7dU7I7FzArDe3d8CiPJSLXB8ynonRcsTlkSvjZVO5NbzgROB18zs19FnNbY68dkqiPcJILmojnwIuD/NqmcCRwGvmNnfzWxY1oPrpDS1ywl07HMFMX3PUo+pk3VLzj5fKqylNUcAt+U7iFyI/ut9COgN9MtzONl2OXBDvoPIFQ+OJnxdfTPw/8zsjDyHJW0bA2xIWdYAjG5nvXTrxFHa3OruPwYGAR8ATiH0rsVWJz5bhfo+AbwX+HvqQne/FxhMKEorgL9Gvalx1tHPVbp1C+U9a7VuyeXnK+6/CJIHZjYZqHf3nnKS16vAb4CTgZvyG0r2mNl7gBp3X5fvWHLN3Te5+1XA94CvtLe+5JUDb6csK6N5zGtr66VbJ1bay63uvsfd/wa8BzjbzEblNMAu6MBnq+DeJ4CoF3q7u29J93z0j8VThN7daYQe7Djr6Ocq3bqxf886Urfk6vOlwlpaMLNewOeAK/IdS664+5poaMRlpP9arFhcBtxhZuvNbD0wDphlZt/Ic1y59FNgfL6DkDbVAUNSlg0EVrezXrp1YqMzudXdFwGPU1i/q619tgrqfUryAeCv7a3k7muBu4n/e9XRz1W6dWP9nnW2bsn250uFtaS6DPiBu+/IdyB5UENIKMXqAuCQpLaacJZ0jxjyE2kEXsh3ENKmJ4DR0cnFmFlvwh/AJ1PWexSYmvR4CjGamSGNzubWbcDLWYwn01r7bBXa+5TQocI6sgeYn8VYMqGjnysovPesK3VL1j5fKqzbEZ2gAWBtrlgEzOxbwDygr5lNMrNPR1+vFCUz65cyj+V7gZ/kK55sc/d6d3890Qh/DOpb+6qzGJjZaDO7wMxKos/y14Ar8x2XtM7d3yDMuHBytOgU4F53rzOzy8ws8Qf/NsJXuok8fRKh1zR2WsmtB5jZd81sZLTORWY2PLr/buBJd38zj2G3qa3PVqG+TwnRiW37JA+bSz4mMzvHzCZF96cRZtBYkp9o00utXdr6XEXrx/49S1ePtVa3WJCfz5e7q7XSCF9/XEIYb3QNUJ7vmLJ4rFdFx5ncFuc7riwf85FAPfA0YWz1ufmOKcfHvwI4Id9xZPkYZxDOEn+J8MdiWr5jUuvQ+zYE+CWhULsZGBAtn0uY1Sax3seBH0Xtg/mOu5VjSZtbgb6EmReOjtabBawH/gh8Mt9xd+C4Wv1sFeL7lHJs5wBfSVnWdEzALcBm4D7gS0BJvmNOiTVt7dLa56oQ3rN0x9RW3ZLPz5fmsRYRERERyQANBRERERERyQAV1iIiIiIiGaDCWkREREQkA1RYi4iIiIhkgAprEREREZEMUGEtIiIiIpIBKqxFWmFmJ5vZHDM7oZ31DjSzP5nZtzO475sytS0RkXwws4PN7J5M5sZW9vOAmZ2XzX1ki5mdYmbPtvd3RgpHab4DEGmPmU0Hrgc2AkcABwOfcPc7s7zr14CZHVhvFWHi/V6Z2Gl0qdldmdiWiAiAmR0BfA84DvhfwsUzRgK9gVvd/ZEs7PZVYDAZyo1t+DnhQjUZF+XjswkXwXkFeDZ6qgz4MHCku6/oxi46+ndGCoQKaykE9wMfcfcFZlYC/DoXO3X3V8xsQwfW22JmazK46+OBpzK4PRHp4dz9eTO7i3CFxDMTy83sfcBfzOwyd/95hvf5VoZzIwBm9n7g34mC1t0fzPQ+Etx9F3C3mV0KzHb3q5LimNWVbUaX1j7R3e929yVmtj5D4UoMaCiIxJqZDQWmANsA3L0RuJpw6dJcaMzweh1xIvBEBrcnIgKwO3WBu/8d+A7wP2Y2Mgv7zGiuNrOxhB7qXNvrZwc8BtR1ZiNmVgb8jnDJ7YRM/v2QPFNhLbHm7huBJcDDZnZ0tGwl8C8z29fMrjWzf5nZB8zsDTNbbGaHJV5vZueb2ZVmNtvM/jtpeaWZfcfM7jaz+8xscLR8HzP7mZldbWa/InyN2Slm9o0opo+Y2RFmdn+0vWujGO82s/Jo+Ztm9l8pm+jj7juibX3NzP7TzP7PzG7o9A9QRKR9vyIMCTkXwMz2N7PvmtntZvZ4VMxiZmVRLrvSzJ42s48mNmBm/2Fm3zazf5rZF1vbUbqcbGbnRDnzIjNbamY3mNloM/t1lE+fNrOTok1cAIwFLjezs6JzYR4zswuT9nGxmV1vZr8zs9+aWX8zG2xm/2VmNWZ2kpktMbPlZlbRlR9Y9DM5IerRxszONrOfmNkPzexhMxsQHcNPzOxv0d+bNdHP+B3AuWb2n0mbHBL9LWows4tb2WeJmX3fzKqifbiZ/SXx/khMuLuaWqwbMANYSfiv/g5geLS8FyFJbSSMddsXeBRYSPin8Z3AVdG6Qwi93u8DBgF3JW3//4Cbovv/A1RF90cAewjJs70Yf0Po9RlOGIvXL1peBjwCzAJGE/4g7AL+C+gHHB49TqyfGD+eOO6fR/cHAlfk+71QU1Mr3AZcBLzeynNrgZ9GefU+oFe0/A/APdH9HwFnR/c/DdRH988HPh7dnxLl6oOix78BvhPdby0nDwU2Az8GjgROAX4AXB2t+3Xgb0mxOlAR3e8LLAYuih6/B/hD0roPAL8gDH39UPT34uzoOB8Gru/gz2529LfiBuAmwnjrU5KeXw+Mj+4vJPxNKgOuBF4HjgI+E+X92Yl4o/VXAD+MnjsfWNVKDJ8FHojuj4h+zqfl+/dKrWXTGGuJPXdfaGYHAdcB/wG838xO8TDmuh7Y4u73ApjZd4BngEmEPyKNZvaVaFOPEBL4+4D9kpYvjl47ALgkei3uvtbMXutEqMOAm4HPuPtb0TZ2mNkbQK27r4728wbwpIfxh/MICX8Y4STI9xH+CAC8BXzSzF4BbiX8UyEikg1OKDaPBCYAl5oZwAZgsIUHnyN0CgD8ljAUAkKurUnJtaOAf6fs4yLS5GR332hmbwKz3P1fAGa2CthmZqMIPbwD0gbt/raZrUta9MVouwm3A3+Klm8m/L24L9pHTRRnRz3h0RhrM7ub0JmTcK67rzKzYwkF9YAo/68Glrn7c8Bz0WvTbftv0d+E59uI6TDg7ei415rZRloOKZEYUGEtsWZmvYBxHk5S+VI0PONBoJrQA5BqYXS7LzAe+L27/yFa9uNom5cTCt0fp+zrcEJCbOhiuFOAdwH7AwuSlqeOMdzT9IR7Y5RkE8Oy9vMw/AV3rzWzi4BbgK8CHwPWISKSQWa2L+HbtoWEvLkhTX4cQehR7QPg7nsInQFEr7nG3edEj1u8NknanBxxWubK14HLgaXAP4GJbRxC8usOAP6R9HgZYZjLcPbOxbvp+owlLxD+AUnYbGbfB+4hxJ6onlOPqz27ab02ewi41cz2IbwXbxN6vyVGNMZa4q5p3B+Auy8AvgVMa2X9PoQk9iqwBjgn+UkzOypa/j4z65uyfEv0cHoXY50D/AT4f1Hi6xQzG03SiTBmNoYwI8oUwsmMf2jlpSIi3XEh4Ruyewj58diopxhoyo/rgZ3A6UnLDzKzgaTkWjMrjToqUrWWk9O5DVjhYVrVHZ04llWEHu6mXQBvEoa6ZIy7N0adHzOif0weBL7r7jWZ3E/KPv8G/I3wzer5wPHu/ma29iddo8JaCsGlUdGZMIYwNi5hqJn1ie6fQBhftxn4I/AhM7vFzE6IehO2ERJgf+ABM3tPdKLINEIxPg+41sz6mNkwwtjm4WbWy8yuMLOPtxKjRe1qwh+fH6V5LnX91MdnEsYDJhwAnOXuW4BL07xGRKQz9uoJNbPjgGuAz7l7PaF3eC3woJm9PzpB8d1RD/VdhNlDzjKzk4EL3b2BkGu/ZGbXmNnxhA6GxDC65PzXWk5OSO49PhQoN7N+wLFAXzObFD23i3CyX6KDJXkftwEfNrMh0eOjgF+6+25CzZM2F5vZ0VFcZa387EpIqZmifz7eR/iWcghhiOFUYBzQL+nEyNRe8Z1p4m8ZVJrxItF7MQh4nNAr35l/OCRHNBRECsE4YKGZ/Z3wNZkBn096fjdwpZltJvTuXgrg7o+a2ZcIXyeeRTj57yUAM/sg8DPCH4o7gK+7u5vZRwhTIS0E7iZ8pXcI4YIKHyEMMWlxYRozO5JwUs54wleD9wPXRPE8AhwNVJjZ/oSelJHAR8zsZZp7by4C+rt7dcqx/9rChR12E04WEhHptCiPXEAo/n4J1BNOgBsMnOHuzwK4+84oP/6SkOtmEcZWA3yZ0CnxO8I3dImZOH5FyH2XEIasXeLu66LC8WhgfzO7s7WcbGbnEMYVf9bMXnb3OsKJlD8gFNW/Az4IHAQsB35P6F3/pIXZomYA7zWzB939PjM7ALjXzP6P8K3nlVHP+rnAyOj4FhIuljPKzGZEcZ5HGDd+f9LPrTch988ExkSFtxPGNp8dbXM+8CTwPGF4y0OEnP4A4YTJmWZ2lrsntns34WI9W6NCfCRwjpnVEnqiiV6fes2GVYSTM88C9gFKzOwxwgmMuZqCVtphei+kkFm4DOxv3L0iR/v7urvrcuMiIkUmGpZS5u6xvECXmX2dMEPLiujxPoRvG76fODdH8k9DQaTQpRtmkZ0dha9MX8jFvkREJHeikzMPi3FRvQ9wFVCRtHgYsE1FdbxoKIgUrGh828cJX+2d6+73ZHmXy6KvKEVEpLiUuvut+Q6iNe6+3cJl1X8andi+mjBs5aq2Xym5pqEgIiIiIiIZoKEgIiIiIiIZoMJaRERERCQDVFiLiIiIiGRA0Zy8OGzYMK+oqOjUa+rr6ykvL89OQEVMP7eu0c+ta3rCz23u3Lnr3b24DzJFV3I2FPfvQ7EeW7EeF+jYClV3j62tnF00hXVFRQU1OLvv2wAAIABJREFUNZ27kmhlZWWnXyP6uXWVfm5d0xN+bma2Mt8x5FpXcjYU9+9DsR5bsR4X6NgKVXePra2craEgIiIiIiIZoMJaRERERCQDenRhXVVVle8QCpJ+bl2jn1vX6OcmyYr596FYj61Yjwt0bIUqm8eW8wvEmNmpwHXAR5Oud3828ENgEHAn8DV33x09Nwy4HtgM9AYud/ddqdutrKz0Yh0LJCLFzczmuntlvuNIRzlbRKSltnJ2TnuszWwEMAA4MmnZeOAs4Bzgq8Bnga8kvewe4A53vxx4BfhuzgIWEenBlLNFRDonp4W1u68FZqUsngB82t3nuvvvgJ8BJwKY2dHAVHd/Llr3b8ClZtY/VzGLiPRUytkiIp2T8zHW7t6Y8vjpxFeIkTpgVXT/JKA2ad266G73vzK94w449FC4+eZub0pEpFjFJWc3zGug5vAaXv7Uy93dlIhI1sRxHusjgJui+2OADSnPNwCjU19UX19PZWVz7q6qqmp7cPrGjTB/PtTWtr6OiEiWVFdXU11dnXg4LJ+xdFNOcnbj241sfWErJX169Dn3IpInHc3ZsSqszWwyUO/uC6JFDrydsloZsNeJMOXl5Z2b7Hvw4HC7ZUsXIhUR6Z7kQtLM1uc5nC7JZc623gZA467GdtYUEcm8jubs2Pzrb2a9gM8BVyQtrgOGJK1jQH9gdbd3mCis33yz25sSEelpcp2zE4W1787tTFYiIp0Rm8IauAz4gbvvSFr2KDA16XEFsBPo/hxNgwaFWxXWIiJdkdOcbaVRYb1LhbWIxFfOC+uoBwPAkpZ9C5gH9DWzSWb2aTOb7O41wBozOzBa9XTgZnff2e1ANBRERKRdccnZJb3DnysV1iISZzkdY21mA4GPRw8vMrNbCF8lXpuy6svufkd0/6PA1Wa2gjCf6pUZCUZDQURE2hSnnN00FESFtYjEWE4La3dvAG6NWsJ1UWvtNSuBz2Q8mMRQEPVYi4ikFaec3TQURGOsRSTG4jTGOrfUYy0iUjA0K4j8f/buPD6uut7/+OuTfV+arfteWiiltBQKlKVCq2wqsgkVFa9QUPGqv4sXRVAvesWLyEUF1IKoIAJXAVE2WWQRKgXaQjfKUpruS9o0adJsbfr9/fGdSSZplkkyycxk3s/H4zxm5syZcz7nzOQzn3zne75HJB4kbmGdkwNJSbBvHxw40P3yIiISNeoKIiLxIHELazMYMsTf3xWXQ8iKiCQMdQURkXiQuIU1QFmZv925M7pxiIhIlzQqiIjEg8QurEtL/e2OHdGNQ0REuqSuICISDxK7sA62WKuwFhGJaaFdQZxTcS0isSmxC+tgi7W6goiIxDRLspZvLNeswlpEYlNiF9ZqsRYRiRvqDiIisS6xC2v1sRYRiRs6gVFEYl1iF9YjR/rbTZuiG4eIiHRLQ+6JSKxL7MJ67Fh/W14ezShERCQMlha4+mKjrr4oIrEpsQvrMWP87caN0Nwc3VhERKRLyVnJABysV2EtIrEpsQvrzEwYOtRf0nzr1mhHIyIiXUjK9F9ZKqxFJFYldmENrd1B1q+PahgiItK1pCz/ldVcp18YRSQ2qbBWYS0iEheSM9UVRERimwrrSZP87bvvRjcOERHpklqsRSTWqbCeOtXfrl4d3ThERKRLOnlRRGKdCmsV1iIicaHl5MU6FdYiEptUWB92GKSkwIcfQl1dtKMREZFOtHQFqVdXEBGJTSqs09J8ce0crFkT7WhERKQTLScvqsVaRGKUCmuA6dP97bJl0Y1DREQ6FWyxVh9rEYlVKqwBZs/2t0uWRDcOERHpVPDkxeZ96goiIrEpZaA3aGbzgR8Cn3bOlQfmZQI3AzVAIXCtc25vd89FTLCwfu21iK5WRCTexVLOTinwX1kHqg9EYnUiIhE3oC3WZlYG5ADHtXvqTuAV59x1wGPAr8N8LjKOPhpSU+Gdd2BvZGt2EZF4FWs5Oznft1gfqFJhLSKxaUALa+fcDnyibWFmw4EFwN8Cs54DPmFmY7t6LqKBZWTAjBn+BEa1WouIALGXs1tarFVYi0iMGvA+1s659medzAV2OefqAs8fANYDp3bzXGSdGljlc89FfNUiIvEqlnK2CmsRiXWxcPLiCGB3u3k1wPBunmujoqKCWbNmtUyLFi3qWRTz5/tbFdYiMgAWLVrUkq+A4mjH0wNRy9kqrEUkWsLN2QN+8mIHHFDfbl46sD9wv6vnWpSUlPDmm2/2PoqTToL0dFi+HCoqoKSk9+sSEenGwoULWbhwIQBmtivK4fRE1HK2CmsRiZZwc3YstFhvwZ85HioX2NrNc5GVmQknn+zvP/NMxFcvIjJIRC1na1QQEYl1sVBYvwAMDwzRhJmlAqOBl7p5LvLOOcffPvxwv6xeRGQQiFrOTslLAYPm6mZcs4vEKkVEImrAC2szs+BdAOfcduBx4PTA/HnAw865LV091y/BXXCBv33qKait7ZdNiIjEk1jK2ZZkpBT6Vuv9uw/pXSIiEnUDPY51LnBV4OFlZhbsyPwl4JNm9h3gLGBhyMu6ei6yRoyAOXOgoQEef7zfNiMiEg9iMWenDU0DoGlHUyRXKyISEQN68qJzrgb4ZWAKnb8HuKKT13T6XL+4+GJ49VX47W/9fRGRBBWLOTttaBp1a+po2t4E0/prKyIivRMLfaxjy2c+4y8Y88wz8OGH0Y5GRERCpJUFWqy3q8VaRGKPCuv2Cgvh05/293s6FraIiPQrdQURkVimwrojV17pb+++G/bti24sIiLSoqWwVou1iMQgFdYdOf54mD0bdu+Gu+6KdjQiIhLQ0hVkmwprEYk9Kqw7Ygbf+Y6/f8st0NgY3XhERASA9NHpADSUN0Q5EhGRQ6mw7szZZ8O0abBli+8SIiIiUZc5PhOAhvUqrEUk9qiw7kxSEnzve/7+f/0X7N0b3XhERIT0kelYitG0rYnm+uZohyMi0oYK666cdx6ceCJUVMDNN0c7GhGRhGfJpu4gIhKzVFh3xcz3sQb46U81rrWISAxQdxARiVUqrLtzwglwySX+MudXXQXORTsiEZGEljE+A4D6D+qjHImISFsqrMNx220wZAg8+yzce2+0oxERSWjZR2QDsG+1rjMgIrFFhXU4Skt9cQ3w9a/Dhg3RjUdEJIFlHxkorFeqsBaR2KLCOlyXXgof/zhUVcGCBbB/f7QjEhFJSNnTAoX1qn04dc8TkRiiwjpcZnDPPTBiBCxe3DoUn4iIDKi00jRSS1JprmmmcZMu4CUisUOFdU8UF8Mf/+jHuL7pJnj00WhHJCKSkIKt1rVv10Y5EhGRViqse+qUU+BHP/L3P/tZWLEiuvGIiCSg3GNzAdj7mi7eJSKxIyKFtZmdFYn1xI3//E/f53rfPt/vevv2aEckIhK2wZCz80/IB2Dvv1RYi0jsSOnJwmZ2BnA9MJTWojwJKAMyIxtaDDODu+6C99+HJUtg/nx46SU/JJ+ISIwYzDk774Q8APa+vpeDBw6SlKIfYEUk+npUWAP3ATcBy4GDgXkGLIhkUHEhIwP+9jffNWTVKjjzTHjuOcjNjXZkIiJBgzZnp5WmkTEhg4Z1DexbuY/cGcq9IhJ9PS2sP3TO3dp+ppmtjVA88aWkxF805qST4PXX4Zxz4PHHVVyLSKwY1Dk7f04+DesaqHqxSoW1iMSEnv529iUzu6CD+Z+IRDBxaeRIeP55GD4cXn4ZPvpRP9a1iEj0DeqcPeSjvvtd5d8roxyJiIjX08J6DnC3mTWHTAeBX/Y1EDM73MzuMLNvmNmdZnZ0YH6mmf3CzH5kZr80s7y+biviJkzwRfXo0fDaazBvHuzeHe2oREQGdc4unF8IQPVL1TTXN/fXZkREwtbTriA3AguBN4FgFjPgsxGI5T7gk865LWY2Gvg7cDhwJ/C0c+6hwIk4vwYuicD2ImvCBPjnP+G002DpUvjIR3yf69LSaEcmIolrUOfstNI0cmbmULuslqqXqig6o6g/NiMiEraetlg/BTzmnPvQObchMJUDd0UglsOBYCe5eiDfzIbjT7L5W2D+c8AnzGxsBLYXeaNH+5brKVNg5Uo49VTYujXaUYlI4hr0OXvImb47yO6/6ldCEYm+nhbWDwLfNrNTQqZTgW9HIJY/43+yzAUuBb4KzAV2OefqAJxzB4D1wKkR2F7/GD7cD703bRqsXeuL640box2ViCSmQZ+zS84vAaDi4Qpcs+uvzYiIhKWnhfUNwNeBe4Hfh0xXRiCWLwMNwBtArXPuYWAE0L4ZogYYHoHt9Z/SUnjhBZg5Ez74wA/J98EH0Y5KRBLPoM/ZOUfnkDEhg/0791P1sk4cF5Ho6mlh/RNglHNurHNuXGAaC5wXgVgygQeAPwK3mdk8wOF/YgyVDuxv/+KKigpmzZrVMi1atCgCIfVBUZEfLeT442HDBpg927dki4gAixYtaslXQHE/bWbQ52wzo/Qify5Lxf9V9GodIiLdCTdn9/TkxU3OuZoO5g/t4Xo68gfgYudclZk5fML+OlDYbrlc4JCOyyUlJbz55psRCCOCCgrgmWdgwQI/vvW8efCrX8EXvxjtyEQkyhYuXMjChQsBMLNd/bSZhMjZJReVsPGmjez8004m3jaRpHRdhVFEIivcnN3T7PPl9jPM7ELgth6up/06ioHpzrng73g/BPYCG4DhZpYZWC4VGA3ET9Nvbi785S9wzTVw4ABcfjlcfTU0NkY7MhEZ/BIiZ+dMzyHn6BwO7D5AxaNqtRaR6OlpYT3EzG4AMLMsM7sbuBvfx64vKoEGMxsRMm83sAJ4HDg9MG8e8LBzbksftzewkpPhJz+Bu++G1FS44w6YMwfWrYt2ZCIyuCVEzjYzhi0cBsC2Rdv6c1MiIl3qUVcQ59zZZna8mf0C+BiwHZgOVPclCOfcQTM7F/iumS0FyoBvOuf2mtmXgJvNbDr+58uFfdlWVH3xizB9Olx0kR/reuZM+M1v4IKOLowmItI3iZSzyxaUse6adVS9UEXde3VkHZY1EJsVEWmjy8LazE7p5KkNwPvALcAY4FzgG30JxDn3Nh2cqe6c2wNc0Zd1x5RZs2DZMl9kP/IIXHghXHYZ3HYb5OdHOzoRiWOJnLNT8lMoW1DGtru3sfl/N3PYLw8b6BBERDDnOh/308xew18EoAo42NliwDDnXHrkwwvfrFmzXMydvNgV53yXkG9+ExoaYNQouOcef4KjiCQUM1vqnJsVgfUkdM7e984+3jjiDSzdOGHDCaSVpUV0/SIi0HXO7q6P9Y3AJOfcmJChmtpPY4HPRDroQc/Mn8S4fDkcdxxs2gTz58NXvgK1tdGOTkTiU0Ln7OzDsyn6ZBGu0bH5Z5ujHY6IJKAuC2vn3JPOuZ3drcQ59+fIhZRgpkyBV1+F//5vf2LjnXfC0UfDK69EOzIRiTPK2TD62tEAbLljC027mqIcjYgkGg32GQtSUuC66+CNN/zJjevW+as1XnON7yYiIiJhyT8hn8KPFtK8t5mNP94Y7XBEJMGosI4l06fD66/D9ddDUhL89Kd+5JA3+joylohI4hh/03gAtty+hYZNapwQkYGjwjrWpKXBD34Aixf7biLvvAMnnADf/S406WdNEZHu5M7MpeTTJbhGR/n3yqMdjogkkD4X1mY2w8z+YGb/YWZj+x6SAP6ExmXL4BvfgIMHfbE9ezbE08gnIhJzEiVnj/vBOCzV2P7b7VS/1qdhu0VEwhaJFuurgHnOuZ8CJwYvZSsRkJkJt94KL74I48bBW2/BscfC5z4Hm3XGu4j0SkLk7KxJWYy6ZhQA73/pfQ4e6Gz0QRGRyIlEYX01cEzg/oNAagTWKaFOOQVWrPBjXqelwX33weTJcNNN0NgY7ehEJL4kTM4ec/0Y0sekU/tWLVvv2BrtcEQkAfS4sDazM83sM4H7x+HHTN0C/jK3zrm9EY5RAHJy4OabfZ/r88+Hujo/kshxx8GGDdGOTkRiVCLn7OSsZCb9fBIA669fT/2H9VGOSEQGux4V1mb238CfgE8BOOdeBz5nZgv6ITbpyPjx8Oc/w7PPwsSJviX7uON8dxERkRDK2VD8iWJKLiqhubaZtZ9fi2vu/GrDIiJ91dMW6xOBkcA/Q+Y9Bfw4YhFJeObN88PwzZsHO3fCaafBv/87VFVFOzIRiR3K2cBhdx5G2tA0ql+pZtOtm6IdjogMYj0trJ91zlUBof/yXxDBeKQnCgrgqafghhv8JdJ/8Qs47DC45x4/koiIJDrlbCC1KJXJ90wGfJeQ2hW1UY5IRAarnhbWH5rZdcBkM7vYzB4GvgL8T+RDk7CkpMCNN/qh+U4+GSoq4Itf9GNfv/ZatKMTkehSzg4oOrOIYVcOwzU5Vl+wmgN7D0Q7JBEZhHpUWDvnHgT+ASQDnwOqgI865+7oh9ikJ6ZPh5degvvvh+HD/RUcTzgBFiyAjbqsr0giUs5ua+KtE8melk39+/Ws/be1OKf+1iISWT0eFcQ595pz7irn3FnOuS8CL5jZ6f0Qm/SUmS+k1671I4akp8MDD/ih+a6/Hmpqoh2hiAww5exWyVnJTH14Ksl5yex6eBebb9P1AEQksqyr/9jN7K9ATjfrKAZ2OedOi2RgPTVr1iz3pq5K2FZ5OXzrW/DQQ/5xWRl8//tw+eW+C4mIxAQzW+qcmxWB9Shnh6HikQpWn78akuGoJ45iyMeGRCUOEYlPXeXs7lqsVwCvA/cBvwfSgNWB+8HpEWBVxKKVyBk7Fh58EF59FY4/HnbsgC99CaZN891GRGSwUc4OQ8l5JYy+bjQ0w+oLV1O7SiczikhkdNds+RN8q3YVgJlNcM59t/1CZvZYfwQnEXLiibB4MTz8sG/BXrsW5s6Fr30NfvQjyMqKdoQiEhnK2WEa94Nx1K+rp+KhClaevZKZS2aSPjQ92mGJSJzrssXaOVcdTNABZWaWHLqMmc0FjuuH2CSSzOCCC2DNGvje9yA5GX72Mzj6aPjXv6IdnYhEgHJ2+CzJmPLbKeQdn0fjxkZWfXwVB2o0UoiI9E1PT178M/C2mf3KzH5iZo8DzwI/j3xo0i/S0nw/6yVLYOpUeP99OOkk33q9Z0+0oxORyFLO7kJyZjJHPnYkGeMyqHmzhlWfXEVzQ3O0wxKRONbT4faeBT4K7AIm44duusA5d1M/xCb96ZhjYOlSuPZa//jnP4dJk+DOO+GAWm1EBgPl7O6llaYx/dnppA1No+qFKtZ8eg0H9+sCWyLSO12OChLWCsySgBnOuaURCcjseOAk4APgFWAfcDNQAxQC1zrn9rZ/nUYF6YO33oKvf731hMYjj4TbboPTE3JELpEBF6lRQcLclnJ2B2pX1fLWKW9xYM8BSj9TyuH3Ho4lWbTDEpEY1FXO7tGYa2Z2TwezS4Bm4NxexNZ+/V8AJjjnrg+Z91vgaefcQ2Z2BvBr4JK+bktCHH00vPACPPIIXHMNrFoF8+bBJz8Jt9wCEydGO0IR6QXl7PDlHJnDUU8fxdunv83O+3diycaUe6ZgySquRSR8Pe1jPQnYCGwImdKBPrd8mNmpwGeAG0LmDQcWAH8LzHoO+ISZje3r9qQdMzj/fHjnHT9SSHY2PPaY74d97bWw95AGJxGJfcrZPZB3XB7THp9GUnYSO+7dwTuXvqNuISLSIz29Ssjlzrl3Q2eY2RTgsgjEciuwGPiFmU0EfgCMwV/IoA7AOXfAzNYDpwLlEdimtJeRAd/+Nnz+8/7qjb//Pdx8M9xzD3zxi3DFFTBhQrSjFJHwKGf3UMGpBUz/+3RWnLmCnQ/u5GDjQY548AiS0np8oWIRSUA9PXnx3Q5mNwKX9yUIM5sMzADucs5dDTwPPA2MAna3W7wGGN5+HRUVFcyaNatlWrRoUV9CkuHD4Xe/g9df9+Ng79oF//M/vlvI8cf7PthbtkQ7SpG4tWjRopZ8hb8aYsQpZ/dO/px8pj83nZSCFHY9uotVn9RQfCKJLtyc3aOTFwMtD6EvSAbKgPucc1f0MlbM7DxgkXOuOPA4E6gAbgTOd87NDll2GfBH59wtoeuItRNhBhXn4LXX4Fe/8heZ2bfPzzeDk0+Giy/23UhKS6Mbp0ic6q+TF5Wz+6ZmeQ0rPrqC/bv2kzMjh2mPTyN9uC4iI5LoInbyIn780z/SmqgPAjucc+/1Ib5gHC2xOOfqzex9/JdAYbtlc4Gtfdye9IQZnHCCn375S3jiCX+p9CeegJdf9tNXvwqnnQYXXuhPelSRLRILlLP7IHdGLjP+NYOVZ66kdnkty2YvY9qT08iZlhPt0EQkRvW009gfnHMvOudeCkz/BCzQv64vVgD5ZlYWMu8APhkPD7SGYGapwGjgpT5uT3orK8sXzw8/DDt3wn33wdln++L72Wdh4UIYNgxOPdVf2XHjxmhHLJLIlLP7KGtiFjP+NYO8E/No3NzI8jnLqXi0ItphiUiM6mlhfUYH83YAf+pLEM65tcBTwAUAZlYApAH3A48DwQGV5wEPO+fUuTcW5OXBpZfC44/Djh1w991w1ln+cukvv+zHxh4zBo49Fm66yV9OvY/jpotIjyhnR0BacRrTn59O6cWlNNc0s/q81Xz47Q85eEAjhohIW2H1sTazq4Ef4n/S68izzrmOEnj4gZgVAz8D3sS3cNzlnFtjZoX4iw2UA0OBbzvnatu/Ph766yWM6mp48kk/LvaTT0JdXetz48b5Fu5zzvGt2hkZ0YtTJEZEuo+1cnb/cM6x+dbNrLt2HTRDwekFHPHAEaSVpEU7NBEZQF3l7LBPXjSzo4DPAz9v91STc25b30Lsu3hM0gmhvh6eeQYefdT3yd61q/W57GyYP98X2Wed5buQiCSg/jh5UTm7/+x5cQ9rPr2G/Tv3kzY0jcm/nUzRGUXRDktEBkhECuvAilKdc/s7mD/aORfVzrTxnKQTRnMzvPGG7zry+OPw9tttnz/mGDjjDD/SyPHHQ35+dOIUGWD9OCqIcnY/adzSyJoFa6h+uRqAEf8+gvE/Hk9yZnKUIxOR/tbrwtrMrgGedM6tCTz+Dv6s7zaLAR9xzs2NTLi9E+9JOiFt2uS7ijz+ODz3HDQ0tD5nBkcdBXPmwEkn+bG0R4/280UGmUgV1srZA8s1Ozb+ZCPlN5TjDjiypmYx5Z4p5B2XF+3QRKQf9WW4vXFAQcjjmUA2sJ3W4ZsM349OpGdGjYIrr/RTXR288IKfXnkFli3zLdpvvw133umXLyuD446D2bP97bHHQkFB19sQSSzK2QPIko0x3xrDkPlDWPOZNdStrmPZCcsY8dURjPvhOFJyejqirYjEu552BZkJrHLONbWbf0SwhSRaBkPrh4Sor/fdRl55xU9LlkBl5aHLTZ7cWmjPnu1budN0IpHEl37sCqKcPUCa65op/345m27dBM2QPjqdSXdMovicfrmopohEUSQvEDPNObes3conAxOAqCZpGWQyM+GUU/wEfpi+dev8JdaXLPG3y5fDu+/66d57/XLp6TBjRtuW7QkT1IVEEpVy9gBJzkpmws0TKL2klHcvf5faZbWs+vgqhpwxhAm3TiD78OxohygiAyCswtrMTsb/fDg/cIncUEPwZ53/LcKxibQyg4kT/bRggZ/X1AQrVrQW2kuW+CL7tdf8FFRQAEcfDTNn+qJ7xgzf0p2in2llcFLOjp7cGbnMXDKTLb/YQvn3y6l8upI9z+1h+FeGM/Z7Y0ktTI12iCLSj8Idx3oacC++/96edk83Afc7526MfHjhG2w/K0ovVVX5LiShLds7dhy6XGam7zYSLLRnzoQjj9S42hIV/TCOtXJ2DGja2cT6G9az7a5t4CClMIVR/zmKkV8dSXK2Rg8RiVeRGsc6H7jQOXd3JIOLlERI0tILzsG2bf5kyOXL/bRsGWzYcOiyKSlw+OG+yD76aF94T5sGJSUDH7cklH4ax1o5O0bUvl3LB9/4gKoXqgBILUtlzHVjGH7lcJLSe3oBZBGJtoiNY93FBk5yzr3S5xX1QSIlaYmAykp46622BffatR1fcn3oUF9kB6dp03wBnp4+8HHLoNRfJy92sT3l7AHmnGPPs3tYf/16at6oASB9VDqjvzWaoV8YqvGvReJIX8ax/iuQ0936ganOudLeh9h3iZakpR/s2+f7bC9f7of5W7nST7WHXI0ZkpNhypS2xfZRR8HIkTpRUnosguNYK2fHOOccu/+6m/XXr2ffqn0ApJakMvJrIxn+5eHqgy0SB/pSWP8EqAK20DoG6iGLAec45y7oa6B9kahJWvrZwYNQXu4L7hUrfKG9YgW8/37HrdsFBYcW20ceCTnd1TqSyCJYWCtnxwl30FHxSAUbf7yR2qX+n/fknGSGXzWckV8fSfoI/SImEqv6UlgXAzXOucZuNlDknNvdtzD7JtGTtAywujpYvbptsf322x2PtQ1+yL9goR2cxo/3Ld+S8CJYWCtnxxnnHFX/qGLjjzey5zl/nqmlGMXnFTPi6hHkn5SP6VcwkZgyEH2sL4/2CTJK0hJ1wRMlQ4vtFSvgnXdg//5Dl8/K8n21jzii7TRunAruBBOFPtbK2TFo75t72XTzJioeqYBmPy97ejYjrh5B2YIykrOUF0RiQb8V1mY2GvgKcLVzLqqj3ytJS8xqaoL33msttIPTli0dL5+e7vtvty+4J0yAVPW/HIwGqrBWzo4PDZsb2PbrbWxdtJX9O/0/5SmFKQz9/FCG/ttQcqapa5lINEW0sDazZOATwBXAfGBHYD0j+hpoXyhJS9yprPSt2WvWtJ02b+54+dRUOOywQwvuSZM0Qkmc68/CWjk7fh1sPMjOP+1ky+1bqFlS0zI/d1YuQ784lNKLS0kt0D/bIgMtUuNYjwMuB74AFAEbgH8H/g6c4JxbHJlwe0dJWgaN6mo/9F/7gru8vOPlk5P9FSmD3UqCt5MnQ7YuoxwP+mkca+XsQaRmaQ1der6jAAAgAElEQVTbfrONHX/cQXO17yeSlJFE8fnFDPvCMArmFmDJ6ostMhD6VFib2YX4lo7TgA+Au4HfA9c5574R4Vh7TUlaBr19+zouuD/80I9e0pGxY9sW3MGpsHBAQ5euRbKwVs4e3Jrrm9n16C623bONquerWuanDUuj5KISyhaUkXtsrk54FOlHXeXslDBeX4+/BO5i4HLn3HuBlfb9rEcRCV92NhxzjJ9CNTT4PtyrV/uuJcHuJe+/71u5y8vhqafavmbo0Lat28Hb0lKNwx3/lLMHseTMZMoWlFG2oIz69fVs/912dvxhBw0fNrDlZ1vY8rMtZEzIoPTiUsoWlJF9hH61EhlIPekKMhrfCjIWeAw41Tn31cBz5iIxvEgfqPVDpJ39+2Hdurb9uN95x7d619d3/JohQ2DqVD/29tSprZMu696v+qkriHJ2gnDOUfNGDTsf2MnOh3bStK2p5bnsadkUn1dMyXklZE/LVku2SARE+uTFFOBc4DPAKuA3wEedc4v6GmhfKEmLhOngQdiwobXQDr3du7fj15SUdFxwDxkysLEPUv188qJydgJxzY6ql6vY+cedVDxcwYE9B1qey5iQQcmnSig+r5i82XlYkopskd7oz+H2JgJfBa5wzmX1ekURoCQt0kfO+SEAV6/206pV/nbNmo4v6w6+S0lHBXd+/sDGHucGcLg95ewEcrDpIFUvVFHxSAW7/rKrZeg+8H2yiz9VTPG5xRScWkBSWlIUIxWJLwNxgZjPOefu7fOK+kBJWqSfOAcbN7YW3MFpzRp/BcqOjBhxaMF9xBGQmzuwsceJKFwgRjk7wbhmR/XianY9souKRyto3NB6cc7knGQK5xdSdHYRQ84aQvowDd8p0pV+L6wjxczSgDeArznnXjSzTOBmoAYoBK51znX4W7WStMgAO3jQnxjZvuB+5x1/QmVHRo9uLbSDRffhhyf8sIADXVhHUm/ztnJ29DjnqF1eS8UjFez+6272rdzX5vmcmTkUnV1E0dlFfoQRdRkRaaOvo4IMpG/iT7QJuhN42jn3kJmdAfwauCQagYlIO0lJMH68nz7+8db5zc2wfn3b7iSrV/uTJjdu9FPoKCVmfljA0IJ72jR/9Uld+CYeKG/HGTMjd2YuuTNzGf/D8TRsbGD3k7upfKKSPc/voXZZLbXLatnwgw2klqQy5GNDKJxfSOG8QtKH629SpCsx02JtZicCU4DvApcB7wHrgULnXF3gBJxqYKpzrrz969X6IRLjDhzwo5S0b+F+910/gkl7KSn+IjfTprWdxowZdEMCxmuLdV/ytnJ2bGqub6bqhSp2P7Gb3U/sbtNlBCBrahZD5vtCO/+UfFJyYq19TqT/xXyLtZllAxc45/6fmX03MHsusMs5VwfgnDtgZuuBU4HyqAQqIr0XLJQnT4bzzmudv3+/H3M79KTJlSvbznvwwdblc3PhqKPguONap3HjBl2xHeuUtwen5Mxkis4qouisItztjrp36qh8ppI9z+6h6sUq6lbXUbe6js23bcZSjbwT8iicX8iQ+UPInZWrqz9KwouJwhq4Fvhxu3kjgN3t5tUAwztaQUVFBbNmtf7zsHDhQhYuXBjJGEWkP6Sm+hMbjzgCLrywdX5dnT9BcuVKWLHC365cCTt3wquv+imoqKhtoX3ssTE/9vaiRYtYtKhlxLviaMbSS33K28rZsc/MyD4im+wjshn19VEcbDrI3n/tpfJZX2jXvFlD9cvVVL9cTfkN5STnJZN/cj4FcwsomFtAztE5JKVotBEZHMLN2VHvCmJmHwPSnXN/DTwux/+kOAu40Dk3O2TZZcAfnXO3tF+PflYUSRA7d8KyZfDGG/D667BkCVRUHLrcxIlw4okwZ46fDj/c9wuPQfHWFSQSeVs5O/7t37Ofqn9UtRTaDR+2PWk5ObddoT1DhbYMDrHeFeQaYEbI1aAK8VcJuzVwP1QusHXgQhORmFNaCmec4SfwwwFu2OCL7OC0dCl88IGf7g2MKldQACec0FpoH3ccZEV1KOd4prwtpBamUnJ+CSXn+1+HGjY1UPVSFVUv+qlhXQOVT1ZS+WQlECi0T2pXaKeq0JbBJRZarEuA0NOM/wX8P2AJsAYocc7Vm1kqUAuMd85tab8etX6ISIv9+333kWCXkVdf9Re/CZWSAjNmwGmnwemnw0knQWZmVMKNwxbrPudt5ezBr2FzA9UvVbcU2vUf1Ld5PikribzZeeTPySdvTh75J+STkh8L7X0iXYubcayh9SfFwHioDwJ/cM49bmZnAp91zi3o6HVK0iLSpY0bfYG9eLG/ffttPxZ3UFqab8k+/XQ/zZrli+8BEG+FdXu9ydvK2YmncUtjmxbt+vfbFtoYZE/LJn9OfkuxnTEmA9OJyRJj4rmwLsRfaKAcGAp82znX4bWVlaRFpEdqanyR/fzzflq+3HcrCRoyxI/P/alPwfz5/dptZJAV1mHlbeVsadrRRPXiaqpfrWbvq3upWVqD29+2JkkbkdZSaOfPySd7erb6aUvUxVVh3VtK0iLSJ7t3wwsv+CL72Wf9mNtBWVlw5pnw2c/CWWf5kUwiKN4L695Qzpb2muub/UgjrwSK7cV7ObDnQJtlkrKSyD0ml7zj88ibnUfu7FwyRmZEKWJJVCqsRUR66p134NFH/RSaW0pKYMECuOwyOProiGxKhbXIodxBR93aOqpf9YV29SvVNKxrOGS5tOFp5M3Oay22Z+WSnJ0chYglUaiwFhHpi02b/EVqfvc7P7Z20IknwrXXwjnn9GkoPxXWIuFpqmii5vUa9i7Z2zI1Vze3XSgJso/Mbim082bnkXV4FpakvtoSGSqsRUQiwTk/lN/vfgf33w9VVX7+8cfD3XfD1Km9Wq0Ka5HecQcdde/VUbOktdiufbsW2tXaybnJ5B6bS+6s1iljrE6MlN5RYS0iEmm1tXDXXXDzzbB9u++H/cQTMHduj1elwlokcprrmqlZVtOm2G7c2HjIcilDUnyRfUxrsZ0+Kl3FtnRLhbWISH+proarrvJdRUpLYf36Ho8gosJapH81bmuk5vUaapbWUPOmn/ZX7D9kudSS1EOK7bThaSq2pY1Yv/KiiEj8ys+H++6Dp57yl1tfuxZmzox2VCISIn1YOumfTKf4k8UAOOdo3NzYUmSHFtuVT1VS+VRly2vThqa1FNk5x+T4lu2h6Z1tShKcCmsRkb6oqICFC33LdXExTJ4c7YhEpBtmRsaoDDJGZVDyKX9JduccDRsaDim2m7Y3sfvx3ex+fHfL69OGp/lW7WMCxfYxuaQPU7EtKqxFRHpn3Tr47W/h9tt9UV1QAA8/DNnZ0Y5MRHrBzMgcm0nm2ExKLygFAsX2h+2K7aU1NG1tYvfW3ez+W0ixPSytTaGde0wu6cNVbCcaFdYiIt1pboZ33/XjWb/+ur+AzHvvtT7/sY/BHXfAhAnRi1FEIs7MyJyQSeaETEo/HSi2DzrqP6j3/bWX1lC7tJaaZTU0beugZXtoGjkzc9q0bqeP0AmSg5kKaxGRULt2werVrdPKlf5y57Xtrsqdn+/Hr164EE4+GfRFKZIQLMnIOiyLrMOyKLukDAgU2+vqWwvtpTW+2N7eROWTlVQ+2dpnO7U09dBuJCNVbA8WKqxFJPE4B9u2wfvv+ysshhbSO3d2/JrRo2HWLDjmGDj1VJg9G1KUQkUkUGxPyiJrUhZlFweKbeeL7dpltW1at/fvPPQEydSS1EO7kWjov7ikbwURGZycgx07fPEcnD74oPV2376OX5edDUcc4S/2EpyOOcYPpSciEiYzI2tiFlkTsyi9KKTP9vqGtt1IlgZGI3m6ksqnQ4rt4tRDupFkjNFFbWKdCmsRiU/Bwrm83E8bNhx6v76+89cXFcGkSX4Uj9AietSoPl2eXESkM2ZG5vhMMsdnUnphSLFd3tC2G8nSGvbv2s+eZ/aw55k9La9PKUohd2bbbiS6gmRsUWEtIrGppga2bIHNm1unTZtai+eNG6Ghoet1FBb64jl0mjjR3xYWDsReiIh0yczIHJdJ5ri2o5E0bmxsKbJbupHs2s+eZ/ew59mQYnuIL7ZbupHMzCVjvIrtaFFhLSIDyzmoqmpbMHc07d3b/bqKimDs2NZpzJi29/Py+nVXRET6g5mRMSaDjDEZlJzXOs5246bGDruR7HluD3ueCym2C1IO6UaSOSFTxfYAUGEtIpFx4IA/8W/bNti+ve1t+3ndtTQDZGT4bhkjR8KIEf525Mi2hXNOTn/vlYhITDAzMkZnkDG67UVtGjc3HtqNZOd+qv5RRdU/qlpen5yffEg3kswJmViSiu1IUmEtIp1zDvbs8QXzzp2+KG5fNAdvKyr88uHIyWktmjubCgs1hJ2ISBfaXEHy3NZiu2lr0yHdSJq2N1H1QhVVL4QU23m+2A5t3c6cpGK7L1RYiySahobWQnnnTn8CYOjj9tOBA+Gt18yPnDF0KAwb1nobej94q5ZmEZF+YWakj0gnfUQ6xZ8obpnfuPXQbiRN25qoerGKqhdDiu3cZHJmtO1GknVYlortMKmwFol3zc1QWdl5Ydy+cK6p6dn68/J8wVxW5m87K5ZLSiA1tX/2UURE+iR9eDrpw9Mp/nhIsb0t0I0kZKztpi1NVL9cTfXL1S3LJecmkzc7j7w5eeSfmE/e7DxS8lVCdkRHRSTWNDT4bhW7dvnb7u7v3h1+FwzwxW9paesULJg7mkpKfF9nEREZdNKHpZN+TjrF57QW2007Du1G0ri5se0JkgbZR2aTPyefvBN9sa2RSDwV1iL9yTmoru5Zodz+0tnhKCw8tEDurGDOz1ffZRER6VBaWRpFZxVRdFZRy7zGbY3sXbyX6sXV7F28l5qlNexbuY99K/ex9VdbAUgtS/Wt2SfmkX9SPrnH5JKUmnjXBIiZwtrMzgN+CuQBfwD+wzl3wMyKgf8GqoBU4Frn3P7oRSoJbf9+30Lck0I53D7KQampUFzsW4tLSrq/X1SkS2vLgFPOFkkc6cPSKTm/hJLz/QmSzQ3N1C6tpXpxNdWv+mJ7/4797Hp0F7se3QVAUnYS+SflU/iRQgrmFpBzTA5JKYO/0I6Jb2MzGw2cC1wATAV+CWwCbgH+BHzLObfEzK4CbgS+Ha1YZRBxDurqwiuOg/erqrpfb3u5ueEVyMHbvDy1KEtMU84WSWzJGcnkz8knf04+fNOPRFK/rt63ar/q+2fXra1jz9/3sOfvvvtIck4y+SfnUzC3gIKPFJAzY3AW2jFRWANjgH9zzh0AlprZkcBHzOwVYLJzbklgub8B75rZD51z+6IVrMSwfftai+CdO7u+rajo+pLXHUlK8i3E4RTKwdZk9VGWwUc5W0RamBlZE7PImpjF0M8NBaBxe6MfceQFP+pI/Xv1VD5VSeVTlYC/iE3h/EKGnDGEIR8bQvqI9GjuQsTERGHtnPtnu1lbgFzgNGB9yHJbAh3jZwEvDViAEj319R0XxJ0Vy3V1PVt/Rkb4XS5KSqCgAJKT+2dfReKEcraIdCd9aDplF5dRdnEZAI1bGql6yRfae17YQ8O6Bir+VEHFnyoAyD4q2xfZZwwhf04+SWnx2ZodE4V1B44FfgIsBHa3e64GGD7gEUlkBEe86KoVOXTevh42cqWnt45m0dlt6P3sbHW7EOk75WwR6VL6iHTKFpRRtsAX2vXr6qn8u2/B3vOPPexbsY99K/ax6eZNJOcmU3ROESXnlzDkjCEkZ8dPg1bMFdZmNhGocM69bWYOaP9bfTpwyIkwFRUVzJo1q+XxwoULWbhwYb/GKgE1NYdekW/7dj9+cvuiuadjKKeldV0kt7/NyVGhLHFh0aJFLFq0KPiwuKtlY5lytoj0RuaETEZ8eQQjvjyCg40HqX6lmsqnK9n91G7qVtex84Gd7HxgJ0mZSQw5Ywgl55dQdE5R1MbPDjdnm+vJ+Lf9zMySgR8D1zvnGs3s28BHnHMfDTxvQENg3uLQ186aNcu9+eabAx7zoHXggC+EQwvlzu73pFU5NfXQVuOuWpZ1Ip8kADNb6pyb1f2SsUU5W0T6Q/26eioerWDXw7vY+9relvmWbhSfW8zQzw+lcH5h1E5+7Cpnx1qL9TXALc65xsDjZ4GrQp4fCzQBysa95Zy/St/mzW2nLVtg69bWormiIvyLjmRkHHolvqFDW8dRDi2WNYayyGCinC0iEZc5IZPR14xm9DWjadzS2FJkV71URcVDFVQ8VEHasDTKLi1j2BXDyJqUFe2QW8RMYW1m1+GTb6aZjQfmAi8D28zsSOfcKuAM4BfOuaboRRrDDh703S+2bDm0cA4toBsawltfScmhl60OFs2h99WqLJJwlLNFZCCkj0hn5NUjGXn1SBo2NrDjvh1s//126t+vZ9NPNrHplk0MOWsII782ksJ5hVG/+mNMFNZmdj3wg3az1zrn7jGzTwPfNbNyIAf4zkDHFzOChfP69VBefujtxo3+AibdycuDkSPbTiNG+ClYLJeW+m4bIiLtKGeLSDRkjM5gzHfGMPq60ex9bS/b7t7Gjvt3UPlEJZVPVJJ1RBZjbhhD6UWlWFJ0CuyY6mPdF4Omv15jI3z4Ibz7Lrz3ni+ag4Xzhg3dtzYXF7ctljsqoHNzB2RXRCQ88drHui8GTc4Wkahqqmhi26JtbLlzC01b/Y9j2UdmM/bGsRSfW9wvLdjx1Mc6MTjnu2QEi+fgbbCQPniw89cWFcG4cTB27KG3Y8ZAVuz0MxIRERHpT2klaYz5zhhGfXMU2+/dzoYbN7Bv1T5Wn7eagtMKOOzOw8iaPHC1kQrr/lZTA6tWwYoVflq50t9WV3e8fFISjB8PkyfDpEkwcWJr8TxmjFqbRURERNpJSkti+OXDGfrZoWy9ayvl3y+n6h9VvHHUG4y5fgxjrhuDJfd/9xAV1pFUXw/Ll8Nrr/lp6VLfraMjRUUwZYovoA87rPV2wgR/kRMRERER6ZGk9CRGXj2SskvKWHftOrb/Zjvl3y2n6sUqDr//cNKH9m+NpcK6L6qr4aWX4PnnYfFieOstP/5zqLQ0OPxwOOqo1mnaNH+CoEbSEBEREYm41KJUptw9hdJPl/LOpe9Q9Y8qls1exvRnp5N1WP91DVFh3RPOwbJl8Nhj8Nxz8Prr0Nzc+nxSki+ajz/eT8ce61ulNbqGiIiIyIAbMn8Is96axerzVrP3tb0sP2k5R790NNmHZ/fL9lRYh+Ott+CBB+DPf27btSM5GebMgXnz4JRTfCGtPtAiIiIiMSN9WDrTn5vOqvNWseeZPaw8ayUzl8wkrTQt4ttSYd2Zxkb405/gjjt8f+mgsjI47zw46yxfTOflRS9GEREREelWcnYyRz56JG/NfYuaN2pY+/m1THtyWsSH41Nh3V5zM9x3H3zve/6CK+Avw/2Zz8BFF8FJJ/mWahERERGJG8lZyRz5lyN548g3qHy6kh337WDo54ZGdBtJEV1bvFu1CmbPhi98wRfVU6fCXXf5MafvuANOPVVFtYiIiEicSh+ezsT/nQjA+hvWc7Cpi2uH9IIK66D/+z+YOdMPkTdqlG+1fvttuPxyyO6fDu4iIiIiMrDKPltG1hFZNG5sZMf9OyK6bhXWQccfD5mZcOWVsHo1XHqpWqdFREREBhlLMkb/52hSClJorm3u/gU9oD7WQaNH+0uKl5VFOxIRERER6Uell5RScmEJyVmRbURVYR1KRbWIiIjIoJeUlgSRH21PXUFERERERCJBhbWIiIiISASosBYRERERiYCELqwXLVoU7RDiko5b7+i49Y6Om4QazJ+Hwbpvg3W/QPsWr/pz31RYS4/puPWOjlvv6LhJqMH8eRis+zZY9wu0b/FKhbWIiIiISIwz51y0Y4gIM6sANvTwZcXArn4IZ7DTcesdHbfeSYTjNsY5VxLtIAZSL3M2DO7Pw2Ddt8G6X6B9i1d93bdOc/agKaxFRERERKJJXUFERERERCJAhbWIiIiISASosBYRERERiYCELKzNLNPMfmFmPzKzX5pZXrRjikVmdp6ZrTez3Wb2MzNLCcwvNrNfm9n/mNmtZpYa7VhjkZmlmdnbZjY38FifuzCY2fFmdo2ZnRv4rOm4JZhw33MzO9PMbjezO83sgoGOs6d6sF8FZrbXzFxg+tRAx9obZjbfzJaY2dgulomr9ywozH2Lq/ets+/4Dpa7zMz+18zuNrNTBjrO3ujBvk0ys+aQ92xGX7fd4YYSwJ3A0865h8zsDODXwCVRjimmmNlo4FzgAmAq8EtgE3AL8CfgW865JWZ2FXAj8O1oxRrDvgmMDXmsz103zOwLwATn3PUh836Ljlui6fZvxcwOA/4bmOWcO2hmL5vZWufcqijEG65wc8AVwOeBmsDjFwYovl4zszIgBziui2Xi8T0La98C4uZ96+Y7PnS5jwAXOefOCjSiLTWz+c65HQMdc7jC3beAy4EzgYPAfufc8j5vP9FGBTGz4cB6oNA5Vxf4L6YamOqcK49qcDHEzE4G/uWcOxB4fDP+A/oD4BHn3PDA/BHAu0CZc25ftOKNNWZ2IjAF+C5wGfAe+tx1ycxOBW4A5rtAYtLfa+IJ9z03s18Blc656wKPvwNMdM59IQphd6sH+5UOPA582Tn3flSC7SUzSwKagXEd/X3G23sWKox9i6v3rbPveOfc2e2Wexr/nb8o8PguYLNz7r8GOuZw9WDfyoDfAFc55zZHavuJ2BVkLrDLOVcHEDjw64FToxlUrHHO/TP4oQzYAmwETsMfr+ByWwJ3Zw1geDHNzLKBC5xz94TMnos+d925FXgH+IWZPW1mc9BxS0RzCe89b5OL8P/gzx2A+HprLuHt1znAbOA9M3vCzIoHNMo+cM4d7GaReHvPWoSxb3H1vnXxHd/CzJLxn8+4es/C2beAS4CPAJvM7LeB7+4+S8TCegSwu928GmB4FGKJJ8cCv0LHLxzXAj9uN0/HrQtmNhmYAdzlnLsaeB54GhiFjluiCfdvpf1ysf65CGu/nHMPA/n4gmYs8NdAa+lgEG/vWdgGwfsW/I4PNQTIIP7fs472DefcbUAe8AlgHr71us/i6U2PFAfUt5uXDuyPQixxwcwmAhXOubfR8euSmX0MeNM5t7PdUzpuXZuK/4l4ReDx7YCh45aIwn3P2y8X65+LsD/LznsZ38I7Bd8SOhjE23vWI/H6vrX7jg8V7Csct+9ZF/sGgHOu2Tn3N+BjwHlmNqyv20zEwnoLUNhuXi6wNQqxxLzAT0FXAt8KzGpz/MzMgGx0/IKuAe4xs11mtgvf4voY/hjpc9e5FEJOpnbO1QPvA8nouCWacHN0++Vi/XPR4++ewAliDwGj+zGugRRv71mvxNP71sF3fKjdQCNx+p51s29tOOfW4H8p7fN7loiF9QvAcDPLBAic5ToaeCmqUcWua4BbnHONgcfPApNDnh8LNAFvDnBcsWoBcHTItBV/1vFv0eeuKyuA/MDJJEEH8MdPxy2xhJuj2+eiw4jhURjo/XdPM/BWP8c2UOLtPeuLeHnf2n/HtwicRP4c8fuedbpvndgHrO3rRhOusHbObcefuXt6YNY84OGQk/AkwMyuA5YDmWY23sz+DagCtpnZkYHFzgB+4ZxrilacscQ5V+Gc2xyc8Mm1wjm3EX3uOuWcWws8hR8eCTMrANKA+9FxSyhd5Wjz45sHv+R/hf/5NvjL2WnAzwc63nCFu19mdoGZjQ/cn4LvIvVuVILuocD7AL4bV3Be3L5nobrbt3h83zr6jg+M63yjmQ0NLHY7cHZg+TTgGODu6EQcvnD2zfz43KWB+ycDLznnqvu87UQbbg/AzAqBm4FyYCjwbedcbVSDijFmdj1+aL1Qa51zh5vZGPwwcuX4sT2/0+4MXAkws3LgMufci/rcdS1wFv3P8L9+jMafyLhGxy3xdPaem9lS4Cbn3J8Dy10KzAy87CXn3GPRiDdc4eyXmd0OXAr8A3gRuD2MESmizsxy8XHfib+2we3OuYp4f88gvH2Lt/ets+94/HuzBrjEOfdaYNn/AIrw3UAedM69OpCx9lS4+2ZmjwFz8L+kPOWcuzci20/EwlpEREREJNISriuIiIiIiEh/UGEtIiIiIhIBKqxFRERERCJAhbWIiIiISASosBYRERERiQAV1iIiIiIiEaDCWkREREQkAlRYi0hMMbMSM/urmW0ws29FOx4REema8nYrXSBGRGKKmX0MeAHIw18t6zDnXGV0oxIRkc4ob7dSYS0iMcvMXgZOc84diHYsIiLSvUTP2+oKIiIxycxGAc8kanIWEYk3ytsqrKUfmNl/mtkd0Y4jKBLxmNk8M/uXmc2NUFiDipndEuH1GXAR8ONIrldE2lK+7jszG2tmd5nZryO0vmQzu8LMNoSx7NlmtqS3+xrJ3K287amwjjNmdriZPWJmd5vZ22bmzOzSaMfVzpPAA9EOIkQk4tkETI9ALDHJzM4ys11mts3MTjCzUYHPmTOz282sILDccWa2ysx+amaZgXllwPYIh7QA+I1z7kBg/SJxR/m6V+IxX+8M3KZHaH1JQBUwOoxl3wQO7+xJM5tqZg8FPnvXm9nQkOcinbuVt1FhHY/+AvyXc+5yYAZwb/AJMzvKzE6NRlCh23bOrXLOvRKNODrSk3g6O4bOuXeBXZGMK5LvV1/X5Zx7ErgTOOCc+5dzbhNwBdAELHHOVQWWex2fyL/pnKsPvPzjwON92oEQZnYD8H3gOTNbC0yO1LpFBpjydQ/FY752ztUBWyO1LefcfmBpmMvuwBfhnT2/GlgUeHiLcy60kI5Y7lbebqXCOo6Y2RDgMGAfgHPuIPBdwJlZPj5pWxTiitq2IymM/Tg4gNuKxroeAkaa2TEAzrndwPPAOSHbKgC2Bj57QVOcc2v7uO0WzrkfOOcmOedmOeemOOdejtS6RQaK8nX/isF8HemRIJERdhIAAAfzSURBVHoSf3fLNgdu2/d7jljuVt5upcI6jgSGrnkXeNrMjg/M2wC8ju/XNA5YaGZfADCzG83shsDPQC39qMxslpndamZfN7M6M3vXzP6fmZ1rZm+Y2afM7FEz22tmV4a8rsP1tdv2FWb2u9C+ZoGfQ39hZjeZ2QtmdlxgfpfbC3m9mdmVZrbazOab2Qoz22Fml4WxjUmh8XSzzUOOYQcKAz/t1pjZFe2O6fcDx+aRQCIO95gdsi0zu8DMXjezy8zsfTP7cRfr6+i97zCergRaNlYDF4bMLgHOtkC3D+ATwF9D4swE6kIeTzGzH5vZLWb2kvkuJZ82szfN7JNm9oCZVZnZ18zsmMA+7jazs7qLTySeKF8nTr7uSGf72NV2zGycmd1jvvX3J+Fsp4PtPmlm95vZnG6WU+7uL845TXE0AVOBDfj/UO8BSkOeKwfmBu4fCXwYuF+C/2+6CP8f90ZgVuC5O4DFgftpQCVwO34sykuALV2tr/228X3Mfg38LjA/A1+sFQQeXwBUALldba+D/R4a2OZ1gWV/CDQAY7rZRvt4utxm6DHsIIZy4KdAVuB1GwPz84AHQ5Z7FZ8UwzpmnWxrCP7nvduA44B5Xa2v3XvfYTxhfr6uB9YF7k8AbgIagfNCPi8WsvwngJMD97OAN4CswOOXA8chD6gGbgQygU8HHl+E/+f+Pwh8BjVpGkwTytcJka8Dz38/zOPY1Xv9BnB4yGtcmJ+z4Pt5AvC9Dp6fG9hOSsg85e5+mtRiHWecb1Wchk80nwNWm1lHJ2m8D3zWzFKBkwPzcvB/yKOAYP/Yd/B/MDjnmoAa4M/Oub34P7Sh3ayvfXyNtD0Z4myg3gX66AKP4P8gz+tme+3XG1znHwPL/hDf/3d+N9toE09PttmJvznfn+4NYFjIPhYFWpS+jj+mEOYx62R/K/EJ7DHn3OvOued6sL7O4gnHQ8B4891BzgfuAv4OXGRmOUCNC2TegDnA4sD9T+C/vIKtIJ8CfhA4zlXAP5zvl/0mkOec+z/nfx5fTs/eA5G4oHydGPm6A53uYxfbmQfkO+eC8bzZw20eC3zGOfdfYS6v3N1PUqIdgITPzJKBUc65cuDfzew3+DOoFwGzQ5d1zjWa2QjgWuD3wVU453aa2ZvA6fj/qMfhi6mWl4bcP0Cgu1Bn6+sk1NB1TAJSQ+I6aGblwMiuttcd51yDmX0IFADFPdhGr7fZzgFa/35GA+udc7e1XyicY2Zmn8UXsEGTnf/J2IXG2oP3oNN4uuOce9/MluO7g+Q75z40s4fwrUgXAk+ExB38bAT7743BtzAF17U7dNUh95tp6yDqliaDjPJ1m/1LhHwdqtPj2MV7cxT+n4jeOg2YZGY5zrnaTpY5GNgH5e5+pAMSX1IJ6f/qnHsb/1PblPYLmtlpwJedcz90foSHUF8BjjU/7NMK4H+623A36+vKRnwLaOgwRIbve9hXacB7/byNcGyjbT9kzGx2D47ZX4GjQ6YOzy7vwfo6jCf83eFB4AvAhyHxJQNfw/9sGjQb3180aAtwkpllh2y3y35+IoOY8nVbiZSvO93HLrazN/CaNHrnFnxr/886eb7atZ50rtzdj1RYx5+vmtnwkMcjgKcD95vwJ2tMwQ/tlG9m6WY2L/B8qZkV4/v6/Qb/M85iAj8tBiTR7r90M7Nu1td+2xayjr8AtcDFgXUVBLYXbPnsbHudGRZYZkTgtU+HsY3QeLrbZuh+tNdR64XhW6GygcfN7GPmT5Lp7j1o2RYwzDm3NmTaH7KJ5JD7Xa0vNO4O4zGzNDNbFOjm0ZX/w/8E/WcA51wN8BTwqms7GshH8d1EgoLH+wHzY2FfDQTHMm3/HrR/n+N6hAKRTihfkzD5Otzj2OF28Lk0Df/PF/hfJzCz4WHk7v/f3t27RhGEARh/Xm1tTGHhF1qrYGGtaKcR0vhRiKAiJIhwhVgIAcHaNGoh2KQQrQRJK9oo6B8gCgo2phU7C1OMxTuB9bjLnceel3DPD6bY2dud3dm7N8N8bILsPb4MnI+IC/XYWxGxSPbWv2583tg9TuOYuG0aTyIXRBTgJ/AUWCaHknbW/XeBVWCOHNr5SvYQzJI9HctkkFohF1L8ruf7Rf7QTgNr5AKZfcCduv/aRufrKrsDvCXnrR2q+46RPZ1L5FzDwzW/b3l97r/U4zvk1IQjjX39yjjQvJ5BZTbrsKvsOXIR36NaF+vHXa37TwAf67NZqvU8bJ3N9bjXc7W858CemrfRM/3rXH2uZy/5jwyeDfFde9C1fRE42ZV3v8dxx4FPwA9gseadIb9rD4Hd1FeOAQvkH94n9V7PTvo3ZjK1lTBeT1O83g+8I3vdjw64x42e9Sw5UviBXEj+BbhR771n7CZHRdbInuoZ4A3ZqL8OzJPv814ip6GsH2PsHmOKWnmaEhGxC1gopdyr29vIXpT5UsriRC9ugIgowMGScxY1ooi4XUoZ6VVOkv4f47WajN1bg1NBps9NclrATCPvFI1FaZtRY+jJYacRRcT2iLgEvJj0tUgaivFaxu4txob19HlMzsX6HBHfgVfAainl/WQvq7+6mKNTN69EvvZN/24H8LKU8m3gJyVtBsZrgbF7S3EqiCRJktQCe6wlSZKkFtiwliRJklpgw1qSJElqgQ1rSZIkqQU2rCVJkqQW2LCWJEmSWmDDWpIkSWqBDWtJkiSpBTasJUmSpBb8ASE+ufoH34BqAAAAAElFTkSuQmCC\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/beagle-2-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/beagle-2-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/beagle-2-mars.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }