{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 18 - Apollo-6 - Earth" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from AMAT.planet import Planet\n", "from AMAT.vehicle import Vehicle" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook simulates the atmospheric entry of the Apollo-6 flight test.\n", "https://en.wikipedia.org/wiki/Apollo_6" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \n", "planet.loadAtmosphereModel('../atmdata/Earth/earth-gram-avg.dat', 0 , 1 ,2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Apollo-6', 5400.0, 400.0, 0.3, 12.0, 0.0, 3.0, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,9.60,0.0,-5.90,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 (2400.0,0.1,60.0)\n", "# bank angle = 70 deg arbitrary, actual profile will give more realistic results" ] }, { "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+17YcXAAAgAElEQVR4nOzdeXiU1f3//+c7YZF9XwSEsIkIRcGAKHVBQHEBkeLe1qpIv7Ufrbb60+7WpVpa9z2urVpEwZZFRcUVq7K5oIIiCChbSFjDFkhyfn+cmWQSkpBlZu7JzOtxXeeayT2Te94T4snLM+c+x5xziIiIiIhI7aQFXYCIiIiISDJQsBYRERERiQIFaxERERGRKFCwFhERERGJAgVrEREREZEoULAWEREREYmCekEXEC1t27Z1GRkZlT4nJyeHdu3axaeggKXKe9X7TD6p8l4j3+fixYtznXPJ/6YjVKXPLk+y/H7ofSQWvY/Ek8jvpbI+O2mCdUZGBosWLar0OZmZmQd9TrJIlfeq95l8UuW9Rr5PM1sTcDlxV5U+uzzJ8vuh95FY9D4STyK/l8r6bE0FERERERGJAgVrEREREZEoSKlgPWnSpKBLiJtUea96n8knVd5rqrzPaEuWn5veR2LR+0g8dfW9mHMuvi9oNgq4FTjfObc6dGw8cCfQHHgW+I1zriD0WFvgNmAbUB+4wTm3v+x5MzMzXaLOxRERqYyZLXbOZQZdR3nUZ4uIlFZZnx3XEWsz6wA0BYZEHOsKjAMmANcCE4FrIr7tReBJ59wNwHLg5rgVLCKSwtRni4hUT1yDtXMuG5hR5nA34DLn3GLn3L+AB4HhAGY2FOjjnJsfeu4s4CozaxKvmkVEDsY5R0FBAYWFhUGXElWJ1mcXFUFBQTTOJCISG3GfY+2cKyrz9bzwR4gh64DvQvdPAVZFPHdd6G5CfmQqIqnpqquuon79+jz00ENBlxJ1idJnz5sH6elw8sm1PZOISOwk4sWLg4FHQvc7A5vLPJ4HdCr7TTk5OWRmZha3rKysGJcpIuLVq+e3BCioxnBqVlZWcX8FtI1NZXERlz47Pd3fJtmHAiJSR1S1z06oDWLMrBeQ45z7LHTIAXvKPK0hcMCFMO3atUvYhcRFJLmFg3V1poJMmjSp+Kp3M8uNSWExFs8+O/QjVrAWkUBUtc9OmBFrM0sHfg7cGHF4HdAq4jkGNAHWx7c6EZGK1WTEuq6Ld58dHrFOoR+xiNRBCROsgeuAfzjn8iOOvQH0ifg6A9gHaGhaRBJGKgZr4txnayqIiNQFcQ/WoREMAIs49jvgE6CRmfUws8vMrJdzbhGwwcz6h546GrjfObcvvlWLiFQsmYN1ovTZmgoiInVBXOdYm1kz4MehL39mZg/gP0q8pcxTv3LOPRm6fz7wJzNbjV9P9ffxqFVEpKqSNVgnUp+tqSAiUhfENVg75/KAh0Mt7NZQq+h71gCXx7g0EZEaS9ZgnUh9tqaCiEhdkEhzrEVE6qRkDdaJRFNBRKQuULAWEaklBevY01QQEakLFKxFRGpJwTr2NBVEROoCBWsRkVpSsI698FQQ/YhFJJEpWIuI1JKCdexpxFpE6gIFaxGRWlKwjj1dvCgidYGCtYhILSlYx54uXhSRukDBWkSklhSsY09TQUSkLlCwFhGpJQXr2NNUEBGpCxSsRURqKT00nKpgHTuaCiIidYGCtYhILWnEOvbCwbqoCJwLthYRkYooWIuI1FI4WBdqnkLMmGnUWkQSn4K1iEgtacQ6Pg45xN/u2RNsHSIiFVGwFhGpJQXr+GjUyN/u3RtsHSIiFVGwFhGpJQXr+AgHa41Yi0iiUrAWEaklBev4ULAWkUSnYC0iUksK1vGhYC0iiU7BWkSklhSs40PBWkQSnYK1iEgtKVjHh4K1iCQ6BWsRkVpSsI4PBWsRSXQK1iIitaRgHR8K1iKS6OrF+wXNbBRwK3C+c2516FgjYDKQB7QCbnDO7TjYYyIiiSCZg3Ui9dlNm/rbnTujcTYRkeiL64i1mXUAmgJDyjz0EPC+c+53wAzg0So+JiISuGQN1onWZ7ds6W+3bYvWGUVEoiuuwdo5l43vaIuZWSfgImBW6NBcYKyZZVT2WFwKFhGpgmQN1onWZytYi0iii/tUEOdckZlFHjoZyHXO7Q49XmBmq4CTgP2VPLY6nnWLiFQkWYM1JFaf3aqVv926tbZnEhGJjbgH63J0BjaXOZYHdAIKKnlMRCQhJHOwLkdgfbZGrEUk0SVCsHZA2Wu8G+JHPjjIY8VycnLIzMws/nrSpElMmjQpimWKiJQvPT0dqF6wzsrKIisrK/xl2+hXFTOB9dkasRaRoFS1z06EYL0Of+V4pGbAenwHXtFjpbRr145FixbFpEARkcrUZMQ6MkiaWW5MCouNwPpsjViLSFCq2mcnwjrWbwOdQks0YWb1ga7Auwd5TEQkIYSDdWFhYcCVxEVgfXabNv5206ZonE1EJPriHqyt5CoYA3DObQRmAyNCx0cC051z6yp7LI4li4hUKpnnWCdSn90pNFN7w4ZonE1EJPriOhXEzJoBPw59+TMze8A5lwP8AphsZkcBHYHIiXaVPSYiEri0tDTMDOccRUVFpKUlwoeBtZdofXbLln73xbw835o1i9aZRUSiI67B2jmXBzwcapHHtwJXVPA9FT4mIpIo6tWrx/79+ykoKKBBgwZBlxMVidZnm/lR65Ur/ai1grWIJJrkGFYREQlYMk8HSSTh6SDrNCFQRBKQgrWISBQoWMdH587+9rvvgq1DRKQ8CtYiIlGgYB0fvXv72+XLg61DRKQ8CtYiIlGgYB0fffr426+/DrYOEZHyKFiLiESBgnV8KFiLSCJTsBYRiQIF6/gIB+tvvoH9B2yULiISLAVrEZEoULCOj2bNoGdPyM+HL78MuhoRkdIUrEVEokDBOn6OPdbfzp8fbB0iImUpWIuIRIGCdfwoWItIolKwFhGJAgXr+Dn+eH/7zjvgXKCliIiUomAtIhIFCtbxM3AgtG4Nq1bBihVBVyMiUkLBWkQkChSs4yc9HUaN8vdfey3YWkREIilYi4hEQVqa704LCwsDriQ1jB7tb2fNCrYOEZFICtYiIlFQVFQElARsia0xY6B+fZg7FzZsCLoaERFPfwFERKLAha6iU7COjzZt4MwzoagI/v3voKsREfH0F0BEJAo0Yh1/P/2pv330UdAMHBFJBFH5C2BmZ0TjPCIidVU4WJtZwJUcXLL02WPGQEaG39585sygqxERqWawNrPRZva+ma0ws29DbTUwPTbliYjUDYk4FSTZ++x69eDXv/b3//pXPy1ERCRI9ar5/GeA24FPgHAXZsBF0SxKRKSuSdAR66Tvsy+7DG6/HRYtgilT4OKLg65IRFJZdYP1t865u8oeNLOvolSPiEidlKBzrJO+z27SBG67zQfsG2+EsWOhWbOgqxKRVFXdvwC/MLMJ5RwfG41iRETqqkScCkKK9NmXXALHHANr18J11wVdjYiksur+BRgGPG5mhRGtCHg4BrWJiNQZCToVJCX67LQ0ePppaNAAsrJg9uygKxKRVFXdqSA3A5OARUB4cSMDflLbQsysL/B/wAqgN5DlnPvUzBoBk4E8oBVwg3NuR21fT0QkmhJ0KkjK9Nn9+/spIddf7+dZz58PRxwR61cVESmtusH6VWCGcy4/8qCZPRaFWp4BznbOrTOzrsBrQF/gIWCOc26qmY0GHgUujMLriYhEzd69ewE45JBDAq6klJTqs3/9a/joI5g+3c+1fv99aN8+Hq8sIuJVd2jleeC3ZnZiRDsJ+G0UaukLhC852QO0MLNO+KvXZ4WOzwXGmllGFF5PRCRqdu/eDUDjxo0DrqSUlOqz09Lgn/+Eo4/2a1uPHAmbN8fjlUVEvOqOWP8R/5HfzwAXOmZAR+BXtaxlGn4u4OnAj4GrgJOBXOfcbgDnXIGZrQJOAlbX8vVERKJmz549QMIF65Trs5s0gTlz4OST4fPP4ZRT4JVXoHPneLy6iKS66o5Y/x04zDmX4ZzrHmoZwPgo1HIlsBdYCOx0zk0HOgNlxxvygE5lvzknJ4fMzMzilpWVFYWSREQOzjlXPGLdqFGjKn1PVlZWcX8FtI1RaSnZZ3foAG+9BYcfDkuWwNCh8NlnUTu9iKSgqvbZ1R2x/t45l1fO8Y7VPE95GgFT8B3zPaFRDof/iDFSQ2B/2W9u164dixYtikIZIiLVk5+fj3OOBg0akJ6eXqXvmTRpEpMmTQLAzHJjVFrK9tmHHgoffADjxvm51sOGwUMPwU9/GrOXFJEkVtU+u7oj1leWPWBm5wL3VPM85XkWmO6cuxm4A99hb8BfVR6pGbA+Cq8nIhIVCTq/GlK8z27TBt54A378Y9i1y693fcklsH17ENWISCqobrBubWZ/BDCzxmb2OPA4/qPAGjOztsBRzrltoUO3AjuANUCn0PJNmFl9oCvwbm1eT0QkmvLy/KBw06ZNA67kACnfZx9yCPzrX/DEE9Cokb/fty+89BI4d/DvFxGpjmoFa+fcmcAbZnY/8ClwOHAUUN7OXtWxBdhrZpGXl2wGlgCzgRGhYyPxIyTravl6IiJRk5OTA/jpDYlEfbZn5rc8X7zYz7fesAF+9CMYMwaWLQuyMhFJNpXOsTazEyt4aA3wDfAPoBswDri2pkU454rMbBzwJzNbDHQArnfO7TCzXwCTzewo/LzASTV9HRGRWNi0aRMQfLBWn125vn3hf/+DRx6BG2+El1+GV1+FiRPhppv8vGwRkdo42MWLk/FrlW4Disp5/En80k2HUotOGsA59xnw83KObwWuqM25RURiKYFGrNVnH0RaGlx5JYwfDzff7LdAz8ryU0Quv9zv3NitW9BVikhddbCpIDcDvZ1z3SKWairbMoCLY1+qiEhiSqBgrT67ijp29KuEfPmlD9l798KDD0LPnn7lEC3PJyI1UWmwds694pzbdLCTOOemRa8kEZG6Ze3atQAcGvBcAvXZ1denj98C/fPP/eohAM8843dvHDYMnn3Wh24Rkaqo7qogIiJSxsqVKwHo1atXwJVITfXv7wP1N9/A1VdD8+Z+Heyf/AQOOwyuu06j2CJycArWIiK1tGLFCgB69uwZcCVSW927w733wvr18NhjfuQ6NxfuvNPfHzAAJk+G0IcUIiKlKFiLiNRCUVERq1atAhSsk0mTJn61kI8/ho8+8hc8tm7tp4zccAN07QonnAB33QXffht0tSKSKGodrM1soJk9a2a/MbOM2pckIlJ3LF++nPz8fLp06ZKIG8QcQH129ZjBscf6Cxs3bIAZM+Dcc6FBA79V+m9+4y94POoo+POfYdEiKCpvPRYRSQnRGLH+f8BI59ydwPHhHbdERFLBwoV+E8PMzMyAK6ky9dk11KABjB0LL7wAmzbB1KlwwQXQrBksWeKX7xs8GNq3hwsvhKeegnXazkwkpUQjWP8fcEzo/vNA/SicU0SkTli0aBEAgwcPDriSKlOfHQXNm8N558GUKZCT4zea+fnP/RrYmzfD88/73R67dPEXRv7qV34b9dzcoCsXkVg62AYxBzCz04HWzrnnzGwIsNM5txT8blzAjijXKCKSsObNmwfAkCFDAq6kfOqzY69hQxg92jfn/Moir7/u21tv+bWyv/wS7rvPP//II+Gkk3w78UTt+CiSTMw5V/Unm90G/AqY45ybEDp2B7DEOffv2JRYNZmZmS48ciQiEg8bN27k0EMPpVGjRmzZsoVDDjmkRucxs8XOuajPJVGfHbx9++DDD+Gdd+C99/z9PXtKP6d3b79m9tChvvXrB/WqPewlIvFSWZ9d3f90jwe6AJdEHHsVeAYItJMWEYm3OXPmADB8+PAah+oYU58dsAYNSkanwQfthQt9yH73Xfjf//wI9zffwNNP++c0buznaoeD9tChfqdIEUl81Q3WbzjntplZ5DD3hGgWJCJSV0ydOhWAM888M+BKKqQ+O8E0aOBHp4cNg9/+FgoK4JNP/JJ+H30E8+fDypU+dL/7bsn3de0KgwbBwIElt506+VVLRCRxVHcqyAVAD6AzMA84FzgHuMo592BMKqyiVPlYUUQSw/r16znssMNIT09nw4YNtGnTpsbniuFUEPXZdVBOjg/Y4bC9YAHk5R34vPbtSwftgQOhRw9I0w4VIjEVtakgzrnnzWwo8DPgp8AG4FTn3NxaVykiUoc8+eSTFBUVMW7cuFqF6lhSn103tWsHZ53lG0BhISxf7ke2P/645HbTJnjtNd/CGjf2c7T79y+57d9fo9si8VKtEetyT2CWDpzsnHszOiXVjEY/RCRedu/eTbdu3cjNzeWNN95g5MiRtTpfrEasK3gt9dlJwDlYs6YkaIfb+vXlP79ly9Jhu18/6NPHr0iiwC1SPTUesTazmcDBthJrC+QCgXbSIiLx8uijj5Kbm8vgwYMZMWJE0OUUU5+dOswgI8O38eNLjm/Z4pf2++KLktvPP/fH33/ft0jNmsHhh/vWp09J690b6sBGoiIJ52BTQZaEnvM1UARcAXwCRA4zdAdax6Q6EZEEk5OTw8033wzAH//4RyyxhvvUZ6e41q3hhBN8C3MOsrNLh+0vv/TTSzZvhsWLfSurc2cfssOhu1cvP4e7e3dopP06Rcp1sGD9d/x0kW0AZtbTOfensk8ysxmxKE5EJNHccMMNbNu2jVNPPZWzwpNgE4f6bDmAmV+ur2NHKDtrafNm+Ppr35YvL7m/YoXfjn3dOr/JTVmdOvmQ3aMH9OxZ+n779ppeIqmr0mDtnNte5lAHM0t3zhWGD5jZyUBibjkmIhJFL730Ek899RQNGjTgnnvuSbTRavXZUm1t2sDxx/sWqbDQz+GODN0rV8K338Lq1X4u9/r1B04tAX8BZThkZ2T4bd67dvW33bpB27YK3pK8qruO9TTgMzN7H8gD+gKnAQeMiIiIJJOVK1cyceJEACZPnkzfvn0DrqhK1GdLjaSnl4xCn3566ccKCmDtWh+yv/22JHCH72/d6qebfPFF+edu1KgkaJe97dbNT0GpXz/271EkFqq73N4bZnYqcCUwANgGTHDO6WNFEUlamzdv5owzzmDr1q2ceeaZXH311UGXVCXqsyUW6tUruXDylFMOfHzrVli1yofsNWtK2nff+dtt20pGwsuTluanmnTtCl26+Na5c+nbQw/1m+2IJJrqjljjnFsP/CH8tZmlmdkxzrlyLn0QEanbtmzZwumnn87y5cs56qijmDJlSsJNAamM+myJt1atfBs0qPzHd+woCdmRgTt8f/16PyK+dm3lr9Ohw4GBu+ytVjaReKtWsDazJ8s53A4oBMZFo6DQZgY/BFYA7wO7gMn4jzFbATc453ZE47VERCqTnZ3NqaeeypIlS+jevTsvv/wyzZo1C7qsKlOfLYmoefOSjWvKs2+fD9XffVdyAeXataVvN2zwK51kZ/u1vCt7rfAId9nWsWPJ/WbNNO9boqO6I9a98WufRu4q0wW/VW6tmdmlQE/nXOToylPAHOfcVDMbDTwKXBiN1xMRqciiRYs455xzWLt2LX369OHNN9+kc+fOQZdVXeqzpc5p0KBkfndFCgp8qC4buMuG8B07YOlS3yrTuHH5gbtsEG/bVlvGS+WqG6wnOudKzYoysyPw2+XWipmdBFwMjIo41gm4CPhl6NBcYLqZZTjnVtf2NUVqY8uWLSxevJgvv/yS7OxsduzYQcOGDWnatCkZGRn07t2bAQMG0KJFi6BLlWpwzvHII49w7bXXkp+fz/HHH89LL71Ehw4dgi6tJtRnJ7i5c2HOHNi/3y9T16mTb507+9aypUZSy1OvXsnPqCLO+fnea9f6Ee7y2saN/nb3bj8nfOXKg79uhw4lYbtDB//v1qHDgfdbtVIIT0XVvXixvEsN8oGJwI21rOUu4APgfjPrBdwCdANynXO7Q69fYGargJOA1bV8PZFq27RpE//617+YPn06H3300UGfb2YMGDCAE044gdGjRzNixAgOOeSQOFQqNbFmzRomTpzI3LlzAZg0aRL3338/DeroVVLqsxPb7t1wxhk+VFekWTO/IUtGRvm3zZvHqdg6yMxvmNO6NQwYUPHznIO8vIpDd2TburVkZPxg6tWDdu0qDt7h++3b+6aVUJJDdedYr6L0R4rpQAfgmdoUYWZ9gIHApc65JWZ2PTAHuBXYXObpeUCnsufIyckhM7Nk2/ZJkyYxadKk2pQlUuyrr77ijjvuYMqUKezbtw+Ahg0bMmjQII4++mg6depEixYt2LdvH9u3b+fbb7/l66+/5rPPPituDzzwAE2bNuX000/nnHPO4YwzztBodoLYtWsXkydP5u9//zt79uyhTZs2PPTQQ5x33nkxe82srCyysrLCX7aNxWuoz05sjRrBWWfBf/5T+njr1j5wff+9D3xLlvhWnlatKg7eGRnQpEls30MyMPP/g9K8ud9hsjJ79/opKOGgnZ0NmzaVvg3f37at5HlV0bp11UN4kyb6JCPeqtpnm3OuoscOfLJZFvBvSjrqIiDbObe8hnWGzzseyHLOtQ193QjIAW4GfuScOzbiuR8D/3bO/SPyHJmZmW7Roshde0Vqb926ddx000089dRTFBYWYmacddZZXHLJJYwePZomB/mrtXv3bhYuXMjbb7/NjBkz+PTTT4sfq1+/PqNGjeK8887j7LPPpmXLlrF+O1LGnj17eOqpp/jrX//KutAQ1Hnnncd9990X16kfZrbYOZd58GdW+7zqsxNcURHMmAF33AELFpQc79cPzj3X75R4yCF+U5ZVqw683bOn8vO3bg2HHeYv4KvoVtuTx0Z+vg/YFQXvyPs5Of53oaoaNSoJ2ZGBu7zWrp0fPZfoqazPrm6wPtE5916ZY32AQufciloUeB6+k24ZcewT4AX8iMjhEce/Af7snPt35DnUSUs0FRQU8MADD/DHP/6RnTt3kp6ezuWXX84NN9xAj8quqDmI1atX89///pf//ve/zJs3j6JQT9qgQQNOO+00zjvvPMaOHUtzfb4bU1u3buXxxx/nzjvvJDs7G4BjjjmGu+++mxNOOCHu9cQwWKvPrkM++wyeeAKeew62bCk53r+/D9ljx8JRR5WMVDrnA1l5gXvVKr98XegDtkq1aVM6bHfuXHIBX3gr9PbtFc5iqbDQ/5tXFLzDt+H7e/dW7/xt2pQfussL5c2bazT8YKIZrP/qnPtdmWMtgbedcwNrUeARwDKgo3MuO3RsIfAA8CDQzjm3x8zqAzuBHs65UjOc1ElLtHz88cdcccUVfBxaw+nss8/mb3/7G30O9hlhNW3atImXXnqJF154gXfeeYfwf4sNGzZk9OjRnHvuuZx++um0bt06qq+bqpxzfPDBBzz66KO8+OKL7A39ZRo4cCC///3vOeecc0gL6EqjGAZr9dl10L598Oab8OKLfprItm0lj3Xu7KePnHWW35ylceOKz1NUBLm5fkrJ2rX+NvJ+eK3oyuZ4h5n5kc/IwB0ZvMP327eHFi0UzGLJOdi168DAXV7Lzva/A9WIejRoUPXR8PbtU3OjnloHazP7P/zcuYoWcH3DOTe65iWCmb0CvOycezDU8b8LHAM8CzzrnJttZqcDP3HOXVT2+9VJS23t37+f2267jVtvvZXCwkK6du3KAw88wJgxY2L+2hs3bmT69Om88MILzJs3rzhkp6enM2zYMMaMGcOYMWOiHu6TnXOOzz77jBdeeIEXX3yRFStKBmlHjRrFtddey+jRowPf8CXawVp9dvIIh+yXXoKXXy49X/eQQ2DECBg92k8Z6dOn+oG2qMiPekcG7vXrSy7eC1/Al5NT9XBWr55flq5tWx/G27UruV/ebdu2qRnO4qWwEDZvLj90lxfGd+6s3vlbtqw4dJcN5i1bJsdKKVEZsTazAcAlwH1lHtrnnKvi1PxKz98WuBdYBHQFHnPOLTWzVvjNBlYDHYHfOucO+GdXJy21sWzZMn7yk5+weLHfjO5Xv/oVt956K00D2LZr/fr1TJ8+nRkzZvDuu+9SUFBQ/FivXr0YOXIkw4cP5+STT6Z9+/Zxry/R7dq1i3fffZfXX3+dl19+uVSY7tChA5deeilXXHFFrab0RFssRqzVZyefoiL49FOYPdu3hQtLP965sw/YI0f6wH3oodF77YICH7oiw3bkbWQAz8ur/vmbN/dBu3Xrkp0bW7U68OuyTRu7RN/u3QcfBQ/fz8nxwb2qIldKqcqIeKLO/4/mVJD6zrkDPjQys67Oue9qUWOtqZOWmigqKuK+++7jxhtvJD8/n65du/L0008zfPjwoEsDYPv27bz22mvMmjWLV155hS2REy+B/v37M3z4cI477jiGDBlCjx49Ah99jbcdO3Ywf/58PvjgA+bNm8e8efOKV24BaN++PePHj+fcc8/lxBNPpF4CThSN4VQQ9dlJbONGeOUVvxb23Lk+5ETq188H7JNOgh/+0AeVeMjP99MPcnN9TVW5rU44i5Se7kdBywbuli1LVvo4WNMKGzVXVOSXIDzYKHj4sR3V3IO1WbOqj4a3bu1/H+KhxsHazK4DXnHOLQ19/Xv8ck2lngYMd86dHJ1ya0adtFTX6tWrufTSS3nnnXcAuPTSS7nnnnsS9sLBgoICFi1axNtvv81bb73F//73P/aUWRKgdevWDBkyhCFDhjBgwAD69etHr169EjJM1sTmzZv5/PPPWbJkCUuWLGHBggV88cUXRPZjZsYxxxzDaaedxmmnncZxxx2X8O8/WsFafXbqKiqCL74oCdnvvutHHiP16QMnnFDSMjISI1A65+eR5+T4C/i2bj2wVXR8167av35amg9wlYXvpk19a9KkdCvvWJMmfmpLIvxsE83evf7fuarTUqoy/z8sLc1PKwqH7q5dD1yCslOn6ExFqU2wfhB4zjn3Qejr6UATYCMlyzcZMNQ5d0TtS605ddJSVc45nnjiCa699lp27txJ+/bteeyxxxg7dmzQpVVLfn4+8+fP57333mPBggXMnz+fTZs2HfC8+vXr06dPH/r160efPn3o3r07GRkZZGRk0KVLl4QLnbt27WLlypUHtKVLlxYviSyqxmQAACAASURBVBepfv36DBo0iOOOO47jjz+e4cOH07ZtTJaFjpkoBmv12QL4udkffeTnZ8+b5++XXZqvc2cfsI8/Ho491q840rBhMPXW1L59PpSXDd/btvkpKTt2HLwdbMnCmqhXr/zAXTaIN27s58o3alRyG3m/KscSrAuPGudg+/aqjYRv2uT/3Q+mQQPo2RMGDYJjjoExY6BXr+rXFs2pIIOAL5xz+8ocPzI8QhIUddJSFevXr+eKK67glVdeAWDChAk8/PDDdS6Ilcc5x/fff8+CBQtYuHAhX3zxBV9++SVr1qyp8HvS09Pp0qULhx56KO3bt6dDhw60b9+e9u3b065dO5o3b07Tpk1p1qxZ8W2TJk2oV68e6enpxS08/cQ5R0FBAfv27WPfvn3s37+f/Px88vLy2L59O9u2bSt1m52dzcaNG9mwYUNx21HJZ4WNGzemf//+DBgwgAEDBjBo0CCOOeaYOr+bZQyngqjPFsAH0I8/9iF73jx4//0Dg0iDBjBwoA/Z4dajR/KPvO7fX3EIjzy+a1f5befOA49FXBoTc/XqlR/AGzb0/6bltdo8Vq+en3IReVvesao8Fjl6HP49q+rvW1GRn0JUUODb7t1+etT69f5C3BUr4MsvfVu9uvxzpKXBW2/56VLVEc1gfYlz7p9ljvUBDnfOzapeWdGlTloq45xj6tSpXHnllWzdupWWLVvy4IMPcuGFFyb9nOSdO3eydOlSvvzyS1auXMnq1atZvXo1q1atYv369VF5jbS0NNLS0kpdaFlTDRs2JCMjg549e5Zqffr0oWfPnoEtiRdLMQzW6rOlXEVFsHSpD9nz5/v21VcHPq9NGxgyxIfswYP9KF8c906qs/btqzh0Rx7fs8e3vXvLvz3YsT17qreUXl1VNnQ7F533bQavv+4v+K3e99V+ub0T8B8fTgKyyjzcGrjPOde1emVFlzppqUhubi5XXnklL774IgCnn346jz/+OJ06HbDLcsrJz8/n+++/Jzs7m+zsbDZt2sSmTZvIzs4mNzeXvLw8du7cSV5eXvH9Xbt2UVBQQGFhIYWFhcWb3ITVq1ePBg0aFLf69evTrFkzWrRoQYsWLWjZsmXx/Xbt2nHooYeWaq1atUr6/9kpKwbL7anPlmrbts2vNLJgQUnYLmd2GZ07+4Ad2Tp2jH+94sPl/v3lB/B9+0q3/PwDj1X3sfz8klHiyNHisseq+lg4gtY0JIdHxMuOjIdbeLpN5PSbdu38/xx26wbDhkV/KkhVZ+ZsA/4FdAeGlXlsH/B49csSiS3nHC+++CJXX3012dnZNG3alLvuuouJEyemXHCrSMOGDenVqxe9atKzhDjnigN2vXr1knJEuQ5Sny3V1rIljBrlG/iws2ZNSdBetAg++QTWrfNt5syS7+3U6cCwHc3l/qR8ZiXTNFq0CLqa6CobusuG73itAFJdVQrWzrnPzexk4FznnDpkSXirVq3il7/8Ja+++ioAJ554Ik8//TTdu3cPuLLkY2YJdwFkqlOfLdFg5ldSyMiA887zx4qK4JtvYPFi3z7+2Lf1632bFTHBqGPHA8N2p07JP2dboqO6c64TRbXmWFd4ErMfOufej0I9NaaPFQX87ol33303N910E3v27KFFixZMnjyZiRMnaiRVElas5lhX8nrqsyVqiopg5cqSsB0O3Nu3H/jcNm1gwIDS7cgjK9+aXSTR1HgqiJnNBA629ZwB/QBtASeBev/99/nlL3/JkiVLALjwwgu566676KjJf5Ii1GdLENLSoHdv3y64wB8rKoJvvz0wbG/eDG+/7Vvk9/fqVTps/+AHfqRc4yFS1xzs89uv8XP11lGyBmpZBpwVzaJEqmPVqlXccMMNxRcndu/enYcffpjTTjst4MpE4k59tiSEcFju1QvOP98fc87PzV6ypHT76itYvty3adNKztGoERxxhB/R7tvX3x55pF+HWLPPJFEd7Ffzb0Cecy6/sieZWaDLNklqysvL4/bbb+euu+4iPz+fRo0acf3113PDDTfQWJ8rSmpSny0Jywy6dPHtjDNKjufnw7JlpcP255/7NYk/+cS3SPXrw+GHlw7bffv6Y3V8SXtJApUGa+dcbhXPcw66ylziZN++fTz22GPceuutbNy4EYCLL76Y22+/ncMOOyzg6kSCoz5b6qKGDeHoo32LtHWrD9zLlvk1t8O3a9aUbPwRKS3Nb2gTOcLdt68f9W7WLH7vR1JbrT5MMbOuwC+B/0OdtMRYQUEBzzzzDH/5y1+KdxM89thjueeeexg6dGjA1YkkPvXZUpe0auW3Wz/++NLHd+6Er78uHbaXLfM77YVb5FKA4EfJw0E7MnQnwaa7kmCqHazNLB0YC1wBjAKy8XP6RGKisLCQadOm8ec//5mvv/4agCOPPJKbb76Z8ePHa01qkUqoz5Zk07RpyfJ9kfLz/VKAS5eWhO1ly3wIX7vWt9dfL/09bduWDtrh1qVL3VvmTRJDlYO1mXUHJgKXAm2ANcAY4DXguJhUJylt3759PPvss0yePLk4UPfo0YO//OUvXHjhhaQn6urwIglAfbakmoYNoX9/3yIVFsKqVaXDdvh+bq7f1n3evNLf06yZn0JSdoS7e3ddOCmVO+ivh5mdix/pOAVYAdwF/BP4nXNuTuhpH8SsQkk5u3bt4rHHHuPOO+9k7dq1AHTr1o3f/va3XHbZZdSvXz/gCkUSl/pskdLS00tWKBk7tuR4eJWSyKAdbjk5fnv3hQtLn6tBA3+RZNkRbl04KWFV+f+uPfgtcD8AJjrnlgOYWe13lhGJsGHDBh555BEefPBBNm/eDPgpHzfeeCMXXHCBArVI1ajPFqmCyFVKwtu4h+Xmlg7a4eD9/ffwxRe+RQpfOFl2hPuII6B58/i9JwneQYO1c242MDt00csVZpYBzACKU46ZmYvGFo6SkubPn899993HCy+8QEFBAeAvSvztb3/LmDFjtGOiSDWozxapvbZt4YQTfIu0c6dfd7vsCHfkhZOzyixm2bnzgSPcRx4J7drF7/1I/FR5ppBz7jvgj2ZWDxgHdDGzW4AngFOBrNiUKMlo7969TJs2jfvvv58FCxYAkJaWxo9+9COuvvpqTjjhBF2UKFIL6rNFoq9pU8jM9C1S+MLJstNKvv7aTzdZtw7eeKP097Rpc+AId9++cNhhunCyLrPaDFqYWS/gKuAK51ygO3JkZma6RYsWBVmCVMHnn3/O448/zjPPPMPWrVsBaN26NVdccQW/+MUv6NatW8AVisSfmS12zmUe/Jm1fh312SJxFL5wsrxpJXl55X9P06blXzjZo4cunEwUlfXZtfoncs6tAH5lZotrcx5Jbnl5eTz//PM8/vjjxaPTAAMHDuTKK6/koosu0k6JInGgPlskviIvnBwzpuS4c7B+/YFTSpYu9RdOLlrkW6QGDaB379Jhu39/H8IVuBNHVP4pnHP/isZ5zKwBsBD4lXPuHTNrBEwG8oBWwA3OuR3ReC2JrYKCAt566y2ee+45pk+fzq5duwBo0aIFF198MZdffjmDBg0KuEqR1BStPhvUb4vUhJmfe92584EXTm7efODo9rJl8N135e842agRHHWUn54SXt+7b1+F7aAk2o/9eiAj4uuHgDnOualmNhp4FLgwiMLk4JxzLF68mOeee47nn3++eLtxgBNPPJGJEyfyox/9SKPTIslF/bZIFLVpAz/8oW+RwhdORobuJUv8VJOPPvItrFEjGDQITjzRX4B5/PHQokV830eqSphgbWbHAxuAraGvOwEX4bffBZgLTDezDOfc6kCKlHKtWLGCKVOm8NxzzxVv5ALQu3dvLr74Yi666CJ69+4dYIUiEgvqt0Xip6ILJ7dsgY8/9lNHFi/2bdUq+N//fLv9dr8c4FFH+aB92mlw8sk+fEv0JUSwNrMmwATn3K/N7E+hwycDuc653QDOuQIzWwWcBKwOpFAptmzZMqZPn860adP47LPPio+3b9+eCy64gB//+MdkZmZqZQ+RJKV+WyQxtG4NI0f6FrZlC3z4Ibz3nt9VcuFC+OQT3+6914fqU06BM8+Ec86Bjh2Dqz/ZJESwBm4A7ihzrDOwucyxPKBTXCqSUpxzfPHFF0ybNo1p06axdOnS4seaN2/O2WefzcUXX8yIESOop4ldIqlA/bZIgmrd2ofmM8/0X+/eDfPnw1tvwSuv+BHul1/27aqr4NRT4ZJLYPx40F5stRN4AjKz04BFzrlNZR5y+B3EIjUE9pd3npycHDIjPh+ZNGkSkyZNimapKaeoqIiFCxcyY8YMpk+fzvLly4sfa9WqFePGjWPChAmMGDGChg0bBlipSN2TlZVFVlbxUtJtg6yluqLRb6vPFomfxo1h+HDfbrkFNmyAV1+FGTN80H71Vd+6d4dbb4ULL9Ra2mVVtc+u1TrW0WBmbwADIw61AnYCdwEXO+cOj3juN8CfnXP/LnserYkaHbt372bu3LnMnDmT2bNnk52dXfxY27ZtOeecc5gwYQLDhw/XFuMiURKvdayjJRr9tvpskcSQmwvPPw/33w/h8bNrroG77w62rkRWWZ+dCMG6HX5EI+xD4NfAfGAp0M45t8fM6uM77h7OuXVlz6NOuuY2btzI7NmzmTlzJm+88QZ79+4tfqxbt26MHTuWs88+m5NOOknTPERioA4G61r32+qzRRJLYSE88QT84hdQVORXHjniiKCrSkwx2yAmGpxzOZFfm1khkOOc+87MZgMjgNnASGB6eaFaqsc5x+eff86sWbOYOXNmqU1bAAYPHszYsWMZO3YsP/jBD3QBooiUon5bJPkUFMCeiIlcaWnB1VKXBR6sD+IXwGQzOwroCGgCXg3l5+fz7rvvMmvWLGbNmsWaNWuKH2vYsCEjR45k7NixnHXWWXTqpOuMRKTG1G+L1BHOwaefwpQp8Oyzfu41wG9+A4cfXvn3SvkSLlg75zIi7m8FrgiumrotNzeXV155hZkzZ/Laa6+xc+fO4sfat2/PmWeeydixYxk1ahRNmjQJsFIRqcvUb4vUDc75Na7ffx/eeAPmzoWIvdwYMAD+8hcYNy64Guu6hAvWUnPOOZYtW1Y8Kv3hhx9SVFRU/Hj//v0ZO3YsY8aMYciQIaTpcx4REZGk5JwPzZ995texnj8fFiyAnJzSz+vY0S+zd+GFMGyYVgOpLQXrOm7//v3MmzePmTNnMmvWLL799tvix+rXr8+IESMYM2YMY8aMISMjI7hCRUREJCby8vzFhp9/7rc5D99uLruqPNC2LQwdCiNG+E1l+vVTmI4mBes6aMuWLbz66qvMmjWLOXPmsH379uLH2rRpw5lnnsmYMWM49dRTad68eYCVioiISDQUFMDq1fD1174tX15yPzw3uqyWLeEHP4BBg+DYY33r3l1BOpYUrOuI5cuXF0/xeP/99yksLCx+rG/fvsWj0scddxzp6ekBVioiIiI1kZ8Pa9bAypXw7bclt8uXw4oVsL/cLfKgYUN/seEPfuDbgAH+tksXheh4U7BOUM45FixYwLRp05g5c2apXQ/r1avHKaecUhyme/bsGWClIiIiUhXOwZYtpUNz+Pbbb+H77/1zKtKlC/Tpc2A77DDQmFpiULBOIM45Pv30U6ZOncrUqVNZvXp18WOtWrXi9NNPZ8yYMYwePZqWLVsGV6iIiIiUa/9++O67krAcbuEAHTF78wBpadCtG/TsCT16lNz26uVHpLWAV+JTsE4AS5cu5fnnn2fq1KmlRqY7derEueeeyznnnMOwYcO066GIiEjAnPMXBYYD86pVpQP0d9/5nQsr0rTpgcE5fNutG9SvH7/3ItGnpBaQ7du3M2XKFJ588kkWLlxYfLxdu3ZMmDCBCy64gB/+8IdaEk9ERCTOwnOdy446h1teXsXfa+anZvToUdK6d/fhuWdPvyqH5j0nLwXrOCoqKuK9997jiSeeYNq0aezduxeA5s2bc+6553L++eczfPhwjUyLiIjEkHN+PeeKgvPatZXPdW7WrGSUuXv30iG6Wzd/MaGkJiW4ONi1axfPPPMM9957L1999VXx8eHDh3P55Zdzzjnn0Lhx4wArFBERSS579/rl6SoKz7t2Vfy94bnOkYE5srVurVFnKZ+CdQytXbuWBx54gKysLLZu3QpA586dufTSS7n00kvp0aNHwBWKiIjUXXl5JWs5r1hROjivX1/597ZsWXFw7tpVc52lZhSsY2DNmjXcfvvtPPnkk+wPLTo5dOhQrrnmGsaPH099/dcqIiJSJYWFfr5zOEBHtsrCc716FY86d+8OrVrF7z1I6lCwjqLvvvuOv/71r8WBOi0tjfPPP59rrrmGoUOHBl2eiIhIwtq7F776Cr74wm/PHQ7P33zjLyYsT8OG0Lu3X8u5d++Sec89evg1n3XJksSbfuWiYNeuXfztb3/j73//O3v37iUtLY2LL76YP/zhDxxxxBFBlyciIpIwCgr8tI0vvijdvvmm4mXqOnU6cFOUI47wUza0MYokEgXrWnDO8fzzz3P99dezbt06AM4//3xuuukmBWoREUl52dnwySfw6aclAXrZMti378DnpqX5sNy/P/Tt6+/36eM3RmnWLP61i9SEgnUNrV+/np///OfMnj0bgGOOOYZ7772XYcOGBVyZiIhIfDnnt+P++GMfpD/+2LeK5kB36+YDdP/+8IMf+Ns+feCQQ+Jbt0i0KVjXwJQpU7jyyivZtm0bLVq04B//+AeXXXaZNnMREZGUsGEDfPghzJ9fEqY3bz7weU2bwsCBcPTRMGCAD9BHHgnNm8e/ZpF4ULCuhn379nHttdfy0EMPAXDGGWeQlZVF586dA65MREQkNgoKYMkSH6Q/+MC31asPfF6bNjBokA/Sgwb51rOnn+IhkioUrKsoOzubcePG8dFHH9GgQQPuu+8+Jk2ahGmFeBERSSLOwdKl8MYbMHcuvPsu7NxZ+jlNm8LQob5lZvowfdhh2jRFRMG6Cr799ltOPfVUVq5cyWGHHcb06dMZPHhw0GWJiIhERV4evPoqzJ7tw/SGDaUf79kTjj++pPXrp9U4RMqjYH0QX375JSNGjCA7O5tjjjmGV155hfbt2wddloiISK3k5MDMmfCf//gwHblWdIcOMHIkjBrlbzXjUaRqFKwrsXLlSkaNGkV2djYjRozgP//5D8205o+IiNRRu3fDf/8LzzwDr79esm60mR+JHjcORo/2FxlqWodI9SVMsDaz8cCdQHPgWeA3zrkCM2sL3AZsA+oDNzjn9se6no0bNzJq1Cg2bNjA8OHDmT17NodoHSARESDx+mypmHPw3nvwz3/CtGl+2gf4XQlHjYJzzoGzz4aOHYOtUyQZJESwNrOuwDhgAtAPeBj4HvgH8CJwo3Nuvpn9P+Bm4LexrGf//v2cd955rFq1isGDBzNjxgyFahGRkETrs6V8u3fDc8/Bfff5jVnCjj0WfvpTOP98v5KHiERPQgRroBtwmXOuAFhsZv2B4Wb2PtDHOTc/9LxZwNdmdqtzblesirnuuuuYN28enTp1YtasWZr+ISJSWkL12VLazp1w771w552wdas/1rEjXH65D9SHHx5sfSLJLCGCtXNuXplD64BmwCnAqojnrQstb5cJvBuLWl5//XXuu+8+6tevz/Tp0+nQoUMsXkZEpM5KpD5bShQWwkMPwS23+AsTAYYMgV/9CiZMgAYNgq1PJBUk6rLtg4FHgM5A2b2c8oBOZb8hJyeHzMzM4paVlVXtF92xYwcTJ04E4Oabb2bo0KHVPoeISFVkZWUV91dA26DrqaVA+mwp8cknforH1Vf7UD10KLz5pt8Z8aKLFKpFaquqfXZCjFhHMrNeQI5z7jMzc8CeMk9pCBxwIUy7du1YtGhRrV77L3/5C99//z2ZmZlcd911tTqXiEhlJk2axKRJkwAws9yAy6mxIPtsgX374KabYPJkP2J92GF+TvXZZ2tVD5FoqmqfnVDB2szSgZ8DN4YOrQMOj3jcgCbA+mi/9qpVq3jggQcAePTRR6lXL6F+NCIiCSfIPlv8BYk/+Ql8+qkP0ddc46eBNG0adGUiqSvRpoJcB/zDORdepv4NoE/E4xnAPiDqwxx/+tOf2LdvHz/+8Y8ZNGhQtE8vIpKMAuuzU1lREdx1l99K/NNPoUcPmDcP7r5boVokaAkTrM3sd8AnQCMz62Fml+HXQd0QuuIcYDRwv3NuXzRf+/vvv2fKlCmkp6dzyy23RPPUIiJJKcg+O5WtXQsjRsBvfuN3SrziCh+uhw0LujIRgQSZCmJmfwDKJtqvnHNPmtn5wJ/MbDXQFPh9tF///vvvp7CwkAsuuICMjIxon15EJKkE3Wenqvfeg3PPhU2b/Jbjjz8OZ50VdFUiEsmcc0HXEBWZmZmuJhfCFBQUcOihh5Kbm8tHH33EscceG4PqREQqZmaLnXOZQdcRTzXts1PVI4/AVVdBQYEfsZ4yBdq1C7oqkdRUWZ+dMFNBgvLmm2+Sm5vLEUccwZAhQ4IuR0REpJR//AN+8Qsfqq+/HubMUagWSVQJMRUkSFOnTgXgwgsvxLQ2kYiIJJCnn/Zh2syPWodW+xKRBJXSI9bOOebMmQPA+PHjA65GRESkxOefw//7f/7+Qw8pVIvUBSkdrJctW8aGDRvo2LEj/fr1C7ocERERAJyDK6/0K39cfnlJwBaRxJbSwXru3LkAjBw5UtNAREQkYTz/PLz/PrRv7+dYi0jdkNLB+oMPPgDgpJNOCrgSERERLz8ffh9apPC226Bly2DrEZGqS+lgHV7qSauBiIhIonjkEVi1Co48En72s6CrEZHqSNlgvWXLFlauXEmjRo048sgjgy5HRESE7dshvAHwHXdAvZRfu0ukbknZYL148WIABg4cSD31XCIikgAmT4bNm+HEE7WrokhdlLLBesmSJQAcffTRAVciIiIC69bB3Xf7+5Mn+7WrRaRuSdmh2m7dunHWWWcxbNiwoEsRERHhkEP80npbtsCxxwZdjYjURMoG6wkTJjBhwoSgyxAREQGgTRu4/36/hrWI1E0pOxVEREQkEWkKiEjdpWAtIiIiIhIFCtYiIiIiIlGgYC0iIiIiEgUpFayzsrKCLiFuUuW96n0mn1R5r6nyPqMtWX5ueh+JRe8j8dTV96JgnaRS5b3qfSafVHmvqfI+oy1Zfm56H4lF7yPx1NX3klLBWkREREQkVswlyYKZZpYDrDnI09oCuXEoJxGkynvV+0w+qfJeI99nN+dcuyCLibcq9tnlSZbfD72PxKL3kXgS+b1U2GcnTbAWEREREQmSpoKIiIiIiESBgrWIiIiISBSkRLA2s0Zmdr+Z/dXMHjaz5kHXFCtmNt7MVpnZZjO718zqBV1TLJlZAzP7zMxODrqWWDOzoWZ2nZmNM7O2QdcTbWbW18weNLNrzewhMzs66JqiycxGmdl8M8uIOJYyfVNtJdPPqrzfhbomWf7WmNlJZva1mW01s3uCrqe2kuVvopndamYu1D4Lup7qqJP/IdTAQ8Ac59xUMxsNPApcGHBNUWdmXYFxwASgH/Aw8D3wjyDrirHrgYygi4g1M7sU6Omc+0PQtcTQM8DZzrl1od/l14C+AdcUFWbWAWgKDCnzUEr0TVGSFD+rSn4X6oxk+VtjZi2AYcBQ4Hjgv2Y22zk3N9jKaqXO/000syb4CxdHhQ7V5CLnwCT9xYtm1glYBbRyzu0O/V/1dqCfc251oMVFmZmdAHzonCsIfT0Z/z7PDLay2DCz44EjgD8BP3POvRNsRbFhZicBfwRGuST+D9bMdgHHOOe+MrN2wGfOuU5B1xUtZpYGFALdnXOrU6lvqq1k+1mV/V0IuJxqS5a/NWbWCNgb7lfNbDFwnXPu7WArq5lk+ZtoZteE7j7inNsbaDE1kApTQU4Gcp1zuwFCHcEq4KQgi4oF59y8cEcXsg74Lqh6Yin0f7QTnHNPBl1LHNwFLAPuN7M5ZjYs6IJiZBrwuJk1A34MXBVwPVHlnCsqc+hkUqRvioKTSaKfVTm/C3VKsvytcc7tiQjVTfD97DuBFlVDyfI30cwMuAC4E8g2szr3qVQqBOvOwOYyx/KApBkJq8Rg4JGgi4iRG4A7gi4i1sysDzAQeMw593/Am8Cc0MfJyeZKYC+wENjpnJsecD2xlsp9U3XpZ5XY6vTfmtAI/KtAfaBxwOXUVFL8TXTeUPxUkPuB58zs9IDLqpZUCNYO2FPmWENgfwC1xI2Z9QJynHN1atJ/VZjZacAi59ymoGuJg37AFufcktDXDwAGnB1cSTHTCJgC/Bu4x8xGBlxPrKVk31RD+lklqCT5W7MCeBoYAfw92FKqLxn/JjrntoauKboduOZgz08kqXDx4jqgVZljzYD1AdQSF2aWDvwcuDHoWmLkOmCg/8QI8P++M8zsNufc5ODKiol6RPx36pzbY2bfAG2CKylmngUucM5tMzMHTDGzDOfcrqALi5GU65tqQT+rBJQsf2uccxuAJ82sCH/xX12TzH8T76OOTc9JhRHrt4FOoYsUMLP6QFfg3UCriq3rgH845/KDLiRGLgKOjmjrgYnU4Y8iK7EEaFFm6kcBsDSgemIitHzgUc65baFDtwI7SJJVQSqQin1TTelnlZiS7W/NIvz/xNU1yfw3sQj4OOgiqiPpg7VzbiMwG/8RD8BIYLpzri7+x3NQZvY74BOgkZn1MLPLQh/VJQ3nXI5zbm244a+uz3HO7Qi6tmhzzn2Fn/s3AcDMWgINgJeDrCsGtgB7zaxzxLHNwPKA6ok6KxlOMki9vqk2ku1nVfZ3oS5Khr81ZtbYzAZGHDoDuDeoemoqmf4mmlknM7vIzNJC/538Bvh90HVVRypMBQH4BTDZzI4COgKTAq4nJszsD8AtZQ5/VdevEhZ+CtxrZg3wo3QXlrkiv85zzhWZ2TjgT6ElrzoA19fF0rMsAwAAIABJREFUPwzliVjpBOBnZvaAcy6HFOmboiQpflaV/C7UGUn0t6Y/8LKZfQV8BCxwziXboEVd0wp/EebvgPeBe+rakpRJv461iIiIiEg8JP1UEBERERGReFCwFhERERGJAgVrEREREZEoULAWEREREYkCBWsRERERkShQsBYRERERiQIFa5EKmNkIM/vAzE4+yPP6m9kLZvbHKL7236N1LhGRIJjZADN7MZp9YwWvM9vMLojla8SKmY00sw8P9ndG6o5U2SBG6jAz6wvcht+dbzAwAPiJc+7ZGL/098BRVXjed/hF7dOj8aKh7Zr3R+NcIiIAZjYYuB04EXgdWIXfaKc+8LBz7rUYvOwKoAVR6hsr8RDwRSxOHOqPx+O3B18OfBh6qCHwI2BILTcwqerfGakjFKylLvgvcJ5z7jMzSwOeiseLOueWm9nmKjxvh5ltiOJLnwS8F8XziUiKc84tNLPngSOcc2eFj5vZmcB/zOw659xDUX7N3VHuGwEwszHA5+FA65x7JdqvEeac2w9MNbOrgHecc3+IqGNGTc5pZu2B4c65qc65r80sN0rlSgLQVBBJaGbWGjgc2AV+62vgT0C8tgwtivLzqmI48HYUzyciAlBQ9kBoC++bgLvNrGMMXjOqfbWZdcGPUMfbAT87YC6wrjonMbOGwL+ARhGHo/n3QwKmYC0JzTm3BfgamGNmQ0PH1gALzKylmd1iZgvMbKyZbTSzZWY2KPz9Znahmf3ezN4xs79GHM80s5vMbKqZvWRmLULHDzGzB83sT2b2BP5jzGoxs/8vVNN5ZjbYzP4bOt8toRqnmlm70PHtZvbnMqdo4JzLD53rN2b2azP7n5ndUe0foIjIwT2BnxJyLoCZ9TSzm83scTN7MxRmMbOGob7s92Y2z8zOD5/AzK40sz+a2Udm9n8VvVB5fbKZTQj1mT8zs2/M7A4z62RmT4X603lmdkroFBcBXYAbzGxc6FqYuWZ2ScRrXGFmt5nZv8zsn2bWxMxamNmfzWyRmZ1iZl+b2bdmllGTH1joZ3JyaEQbMxtvZvea2Z1mNsfMmobew71mNiv092ZD6Gd8JHCumf064pStQn+L8szsigpeM83M/mZmk0Kv4czsP+F/H0kQzjk1tYRuQD9gDf7/6p8E2oeOp+M7qS34uW4tgTeAL/H/03gc8IfQc1vhR73PBJoDz0ec/3/A30P37/7/2bvv+Lbqe//jr4/lve3E2XESMklISEjKHmGV0MFogRY64LaQtnT/OriFQm+BbuDS0l5KyqWDltIBtxQKtFDCKqMkhCxCBomzE+wk3vH+/v74SrHs2I6HLMnS+/l4nIekI+no42P5o4+/+pzvARYHr48EWvHJ80gx/go/6jMC34uXHVyfAfwdeAQYg/9AaAa+BWQD84O3Q48P9Y+Hfu7/CV7PA/4z1r8LLVq0DN0FuArY0c19e4GfBPPqw0AguP4B4E/B63cAHwhe/wRQHrx+OfDR4PVpwVw9O3j7V8B/Ba93l5OLgUrgTuB44BzgNuCm4GO/BjwaFqsDJgavZwHrgKuCt88DHgh77GPAPfjW14uDnxcfCP6cTwLf6eW+ezb4WfF94Ef4futzwu6vAEqD19fiP5MygBuAHcAJwCeDef/ZULzBx5cBtwfvuxzY1k0MVwOPBa+PDO7nd8f6faWl46Iea4l7zrm1ZjYbuBW4Fni/mZ3jfM91OVDtnHsIwMz+C3gROAr/IdJmZl8Kburv+AT+XmBY2Pp1wefmAp8JPhfn3F4z296HUIcDdwGfdM7VB7fRaGZ7gC3OuV3B19kDPOd8/+EKfMIfjj8I8r34DwGAeuDjZrYBuBv/T4WIyGBw+GLzeGAC8HkzA9gHFJi/8Sn8oADAr/GtEOBz7bJOuXY0sLrTa1xFFznZObffzKqAR5xz/wYws21AnZmNxo/w5nYZtHMHzeydsFWfC2435F7gj8H1lfjPi4eDr7EsGGdvLXXBHmsz+wN+MCfkUufcNjM7BV9Q5wbz/y7gbefcq8Crwed2te1Hg58Jr/UQ03HAweDPvdfM9tOxpUTigApriWtmFgDGO3+QyheC7RmPA0vwIwCdrQ1eFgKlwP3OuQeC6+4MbvM6fKF7Z6fXmo9PiDX9DHcacCowGVgZtr5zj2HroTucawsm2VBb1jDn219wzm0xs6uAnwJfBj4CvIOISASZWSH+27a1+Ly5r4v8OBI/opoO4JxrxQ8GEHzOt51zLwVvd3humC5zcpCjY67cAVwHbAReASb18COEP28q8EzY7bfxbS4jODwXt9D/GUtex/8DElJpZj8A/oSPPVQ9d/65jqSF7muzJ4C7zSwT/7s4iB/9ljiiHmuJd4f6/gCccyuB64EZ3Tw+HZ/ENgG7gUvC7zSzE4Lr32tmWZ3WVwdvHt3PWF8Cfgz8Lpj4+sTMxhB2IIyZjcXPiDINfzDjA908VURkIK7Ef0P2J3x+PCU4Ugwcyo8VQBOwKGz9bDPLo1OuNbPU4EBFZ93l5K78HChzflrVxj78LNvwI9yHXgKowre6RIxzri04+DEr+I/J48DNzrllkXydTq/5KPAo/pvVy4EznHNVg/V60j8qrGUo+Hyw6AwZi++NCyk2s/Tg9YX4/rpK4PfAxWb2UzNbGBxNqMMnwBzgMTM7L3igyAx8Mb4CuMXM0s1sOL63eYSZBczsP83so93EaMHlJvyHzx1d3Nf58Z1vvw/fDxgyFbjIOVcNfL6L54iI9MVhI6FmdjrwbeBTzrly/OjwXuBxM3t/8ADF04Ij1A/iZw+5yMzOBq50ztXgc+0XzOzbZnYGfoAh1EYXnv+6y8kh4aPH84ASM8sGTgGyzOyo4H3N+IP9QgMs4a/xc+CDZlYUvH0C8AvnXAu+5ukyF5vZicG4MrrZdyl0qpmC/3y8F/8tZRG+xXA6MB7IDjswsvOoeFMX8XcMqot+keDvIh/4J35Uvi//cEiUqBVEhoLxwFoz+xv+azIDPh12fwtwg5lV4kd3Pw/gnHvKzL6A/zrxIvzBf2sAzOxC4Gf4D4r7gK8555yZXYafCmkt8Af8V3pz8SdUuAzfYtLhxDRmdjz+oJxS/FeDfwG+HYzn78CJwEQzm4wfSRkFXGZmb9E+enMVkOOcW9LpZ/+l+RM7tOAPFhIR6bNgHrkCX/z9AijHHwBXAJzvnHsZwDnXFMyPv8DnukfwvdUAX8QPSvwG/w1daCaO/8Xnvs/gW9Y+45x7J1g4nghMNrPfdpeTzewSfF/x1Wb2lnNuJ/5AytvwRfVvgAuB2cBm4H786PrHzc8WNQt4j5k97px72MymAg+Z2b/w33reEBxZvxQYFfz51uJPljPazGYF4/wwvm/8L2H7LQ2f+48FxgYLb4fvbf5AcJtvAM8Br+HbW57A5/TH8AdMHmtmFznnQtv9A/5kPbXBQnwUcImZbcGPRBN8fudzNmzDH5x5EZAJpJjZ0/gDGKM1Ba0cgel3IUOZ+dPA/so5NzFKr/c155xONy4ikmCCbSkZzrm4PEGXmX0NP0NLWfB2Jv7bhh+Ejs2R2FMriAx1XbVZDM4L+a9MX4/Ga4mISPQED848Lo6L6kzgm8DEsNXDgToV1fFFrSAyZAX72z6K/2rvUufcnwb5Jd8OfkUpIiKJJdU5d3esg+iOc67B/GnVfxI8sH0Xvm3lmz0/U6JNrSAiIiIiIhGgVhARERERkQhQYS0iIiIiEgEqrEVEREREIiBhDl4cPny4mzhxYp+fV15eTklJSeQDShLafwOnfThwQ30fLl++vMI5N3R/gH5Qzo4u7be+0z7rn2TYbz3l7IQprCdOnMiyZX0/k+iCBQv69TzxtP8GTvtw4Ib6PjSzrbGOIdqUs6NL+63vtM/6Jxn2W085W60gIiIiIiIRoMJaRERERCQCkr6wXrx4caxDGNK0/wZO+3DgtA+Th37X/aP91nfaZ/2T7Pst6ieIMbNzgVuBD4Wd7/4DwO1APvBb4CvOuZbgfcOB7wCVQBpwnXOuufN2FyxY4BK9p0dEEpOZLXfOLYh1HF1RzhYR6ainnB3VEWszGwnkAseHrSsFLgIuAb4MXA18KexpfwLuc85dB2wAbo5awCIiSUw5W0Skb6JaWDvn9gKPdFo9AfiEc265c+43wM+AMwHM7ERgunPu1eBjHwU+b2Y50YpZRCRZKWeLiPRN1HusnXNtnW6/EPoKMWgnsC14/SxgS9hjdwavRu4r0zfegBNPhI9+NGKbFBFJFPGSs2tW1LBs/jLe+o+3BropEZFBE4/zWL8L+FHw+lhgX6f7a4AxnZ9UXl7OggXtuXvx4sW9a6APBODVV6Gmpr/xioj02ZIlS1iyZEno5vBYxjJAUcnZbQfbqH29lpT0pD/mXkRioLc5O64KazObApQ751YGVzngYKeHZQCHHQhTUlLSvwnJR4/2l7t29f25IiL9FF5ImllFjMPpl2jmbEszANqa247wSBGRyOttzo6bf/3NLAB8CvjPsNU7gaKwxxiQA0SuCh42DNLSoLISDnb+PBARka5EO2eHCmvXHN2ZrERE+iJuCmvgq8BtzrnGsHVPAdPDbk8EmoDIzdFk1j5qvXt3xDYrIpLgopqzU9L8x5UKaxGJZ1EvrIMjGAAWtu56YAWQZWZHmdknzGyKc24ZsNvMjgk+dBFwl3OuKaJBjQm2/6mwFhHpIF5ytqUHW0Ga1AoiIvErqj3WZpYHhKbfuMrMfor/KvGWTg99yzl3X/D6h4CbzKwMP5/qDREPTH3WIiKHiaecrVYQERkKolpYO+dqgLuDS8itwaW752wFPjmogYVGrFVYi4gcEk85W60gIjIUxFOPdeyUlvrLsrKYhiEiIl3TiLWIDAUqrAGmTPGXmzbFNg4REenSoen21GMtInFMhTWosBYRiXOhE8NoxFpE4pkKa4DJk/3l5s3Q2hrbWERE5DBqBRGRoUCFNUBOjp8ZpKkJtm+PdTQiItLJocK6xeGcimsRiU8qrEOOPtpfrlkT2zhEROQwZgYBf921qLAWkfikwjrkuOP85euvxzYOERHp0qE+6yYV1iISn1RYh6iwFhGJa4dmBmnWzCAiEp9UWIeECuvly2Mbh4iIdEkHMIpIvFNhHTJ1KhQWwo4dOlGMiEgc0tkXRSTeqbAOSUmBM87w15cujW0sIiJymJRM/5HV1qBWEBGJTyqswy1c6C+ffTaWUYiISBdSsv1HVmu9zjcgIvFJhXW4M8/0l888A5onVUQkrgRy/Hx7bXUasRaR+KTCOtzs2TBmjO+z1kGMIiJxRSPWIhLvVFiHS0mBD3zAX//zn2Mbi4iIdBAasW6tU2EtIvFJhXVnH/ygv/zzn9UOIiISRwLZwVaQerWCiEh8UmHd2Wmn+XaQt9+G556LdTQiIhKUkhNsBdGItYjEKRXWnQUCcM01/vrdd8c2FhEROUQj1iIS71RYd+Xqq32/9cMPw+7dsY5GRERQj7WIxD8V1l0ZNw4uughaWuCOO2IdjYiIoFlBRCT+qbDuzvXX+8u774aKitjGIiIiBPKCI9bVKqxFJD6psO7O/Plw/vlQVwe33RbraEREkl5acRoALQdaYhyJiEjXol5Ym9m5ZvaqmU0MW5dlZneZ2XfN7G4zy+/NfYPuW9/yl//9336WEBGRJBNPOTu1OBWA5v3NkdqkiEhERbWwNrORQC5wfKe7/gd40Tl3PfAIcE8v7xtcJ5wAH/84NDXBl7+sea1FJKnEW84+NGK9XyPWIhKfolpYO+f24hPtIWY2BrgCeDS46mngAjOb2NN9UQkY4Ac/gLw8ePRR+MMfovayIiKxFm85WyPWIhLvot4K4pzrPAHpQqDCOVcfvL8F2AKccYT7omPUKLj9dn/92ms1/Z6IJJV4ytkasRaReBcPBy+OBfZ1WlcDjDnCfR2Ul5ezYMGCQ8uSJUsiF+HVV8OiRXDgAHzyk9CmkxOIyMAsWbLkUL4Chsc6nj6IWc5OLWofsXZqzRORKOptzk6NXkjdcsDBTusygNB3fT3dd0hJSQnLli2LfHQAZnDvvTB7NjzxBHz/++3T8YmI9MPixYtZvHgxAGY2lOb0jFnODmQFSMlOoa2+jdaaVlLz4+EjTESSQW9zdjyMWO8EijqtywN2HeG+6Bo7Fn77W3/9xhvh6aejHoKISByIac5OH50OQOOuxkhtUkQkYuKhsF4KjDGzLAAzSwNKgeeOcF/0vec9vqhua4PLLoM334xJGCIiMRTTnJ0xNgOApp1NkdqkiEjExGIeawtdBXDO7QEeA84Orj8HeMg5t7On+6IYckff+hZceKHvt160CHbsiFkoIiKDLd5ydsY4X1g37tSItYjEn2jPY50HfDp48yozKwle/wxwoZndALwHWBz2tJ7ui75AAH7/ezj5ZNi+Hd79bs0UIiIJKR5zdmjEWoW1iMSjqB754ZyrAe4OLuHrDwDXdPOcbu+LmawsP6/16afD2rWwcCE884zvwxYRSRDxmLMPFdY7VFiLSPyJhx7roam4GJYuhWOPhQ0bfJG9aVOsoxIRSWgZ44OF9TYV1iISf1RYD0RJiR+pXrAANm+Gk06Cf/0r1lGJiCSsrKlZANRvqI9xJCIih1NhPVDFxb64Pv98qKiAs8/2PdgiIhJxWVOywKBhcwNtzTpZl4jEFxXWkZCXB3/9K3zmM9DYCFdcAV/6EjRpOigRkUgKZAXIKM3AtTgatjTEOhwRkQ5UWEdKair87Gfwk5/46z/+MZx2GmzdGuvIREQSSvb0bADq16sdRETiiwrrSDKDz38eXngBxo+Hf//bH9z461+Dc7GOTkQkIWTP8IV13dq6GEciItKRCuvBcOKJsGIFvP/9UFUFV10FF1wAu6J/JnYRkUSTd1weALXLa2MciYhIRyqsB8uwYfDII/CrX0FBATz2GBxzDNx7rz8luoiI9Evu/FwAapbVxDgSEZGOVFgPJjO48kp/Epnzz/enQb/mGjj1VFi5MtbRiYgMSdkzsknJTqGhrIHmfc2xDkdE5BAV1tEwdiz87W/wwAMwahS8/DIcdxx8+ctQXR3r6EREhpSU1BRy5/lR6+p/K4eKSPxQYR0tZnD55fDWW/CFL/h1d94J06fDPfdAS0ts4xMRGUIKTysEoHJpZYwjERFpp8I62goK/FR8y5b5gxz37IFPfxpmz4a//EWzh4iI9ELhWSqsRST+qLCOlXnz4KWX4I9/hMmT/Uj2xRf7/ut//EMFtohIDwpOKcDSjJrXa2g+oD5rEYkPKqxjyQwuvRTefBPuuguGD/fF9nnnwQkn+LM5qsAWETlMIDtA/kn50AYHnj4Q63BERAAV1vEhPR0+9znYvBm+/30oKYHXXoMLL4S5c/1Bjzo9uohIB8MvGA5AxcMVMY5ERMRTYR1P8vLguuugrMwf2DhmDKxaBR/5CEycCLfeCuXlsY5SRCQuDL/YF9b7HttHa0NrjKMREYlQYW1m74nEdiQoOxu++EV4+20/Y8jMmbB7N9x4oz9V+ic/6c/sKCLSD4mSs7OOyiJ3bi6tta0ceErtICISe30qrM1skZm9aGabzGxzcCkDHhqc8JJcZiYsXgxr1vgDGt/7XmhshPvu8/Ngz58Pd98NlToqXkQOlww5u+RDJQDs+fWeGEciItL3Eev7gYeBa4D/CC5XBdfLYDGDc8/1p0Vfv97Pg11UBK+/DtdeC6NHw8c/Ds89p9Oli0i4hM/Zoz4+CgKw75F9NO3VsSgiElt9Law3O+fucM4tdc49F1yeBW4ahNikK9Om+Xmwd+3yBzWedRY0NMD998PChb4X+7rr4I03NKOIiCR8zs4Yk8Gw9w7DtTiNWotIzPW1sP6MmV3SxfoLIhGM9EFmpj+T4z//6Xuxr78eSkth+3b44Q/9PNmzZsEtt8CmTbGOVkRiIyly9pjFYwDYeddO2pr0rZ2IxE5fC+tTgHvNrDVsaQPuHmggZna0mf3MzL5sZv9jZnOD67PM7C4z+66Z3W1m+QN9rYRz1FHwne/Ali3wwgu+PWT4cFi3Dm66CaZOhWOPhW99yx/0qJFskWSRFDm7+Pxismdm07ijkb2/2zvYLyci0i1zfSiyzOwA8ClgGRCa28iAjznnbhlQIGbLgAudczvNrBT4u3PuaDP7JfCkc+4PZrYIuNI5d3nn5y9YsMAtW7ZsICEkluZmP5r9wAP+VOk1Ne33TZgAF13kl1NPhdTU2MUpIpjZcufcgkHYbtLk7D2/2cNbV75F1vQsjl97PBawiG1bRCRcTzm7r4X1A8B/OOcaO60f5ZwbUHObmdUB851zb5lZCbASWABsAYqcc/VmlgpUAbOcc2Xhz1dh3YPGRli61BfYjzwCe8J+VYWF/sDI887zy7hxsYtTJEkNYmGdNDm7rbmNf0/7Nw1lDUz7xTTGXD0mYtsWEQnXU87uayvIg8A3zOz0sOUM4BsDjhL+jP/KMg/4KPB5YCFQ4ZyrB3DOteCT9hkReL3kkZEBixbBz38OO3fCyy/7AxynTfNT9f3pT3D11X6O7Nmz4atfhaef9gdFishQljQ5OyUthUnfnQRA2Y1ltNS0DPZLiogcpq+F9Y3Al4DfAL8OWz4VgViuBRqA14Ba59xDwFhgX6fH1QAaiuivlBQ48UR/6vT16/2BjT/7GVxwAeTk+Dmzb7/dj2IXF8M55/gzPr7wgh/5FpGhJKly9ogPjyDvhDya9jSx9Ttbo/GSIiId9LW59kfAE865mvCVETqLVxbwe3xivtPMtgAOONjpcRlAc+cnl5eXs2BB+6j84sWLWbx4cQTCSnCTJ/uDHa+9Fpqa4KWX4Mkn4e9/91P2/fOffgE/E8nJJ/tp/c44A044wY+Gi0ifLVmyhCVLloRuDh+kl0mqnG1mTLlzCitOXsH227Yz4tIR5M3PG9A2RUSg9zm7rz3WJznnXu5i/Secc/f1OcqO23gS+LBzrtLMbgS+gB9p+ZZzblrY4zYG1z0Q/nz1WA+C8nJ4/nl49lm/rFnT8f6MDFiwAE46qX0ZPToWkYoMaYPYY52UOXvjlzay88c7yZmTw/zX5pOS3tcvZ0VEuhfJgxfvd859rNO6S4H/dc71e0olMxsOrHbOjQ7eNmATcCXwJFDinDtoZmlALXCUc25n+DZUWEdBRUXHQnv16sMfM2FCe5F98sl+mr+0tGhHKjKkDGJhnZQ5u7WuldfmvEbD5gbGfnEsU++cOiivIyLJKZIHLxYHRyYws2wzuxe4F99jNxD7gQYzGxu2bh+wCngMODu47hzgoc4JWqJk+HD4wAfgJz+BVatg3z54/HG48Ubfi52bC1u3woMPwhe/CO96F+Tn+57uz34W7rvPP69FBxWJRElS5uxAToCZD8zE0oydP97JO398J5ovLyJJrE8j1gBmdiLwEeA8YA/wcaDKOXdgQIGYHYs/GGY5MBJ43jn3nJkVAT8EyoBRwDecc7Wdn68R6zjQ2gpr1/pZR0LLhg2HPy4z049kz5/fvsycqZFtSVqDNWId3HbS5uwdd+1g0xc2EcgNMPeFueTNVb+1iAxcv1tBzOz0bu46HjgTuA1oAy5yzn15oIEOhArrOHXgALz+Oixf3r68/fbhj8vM9KdgnzOn4zJ8sI7pEokfkSqslbM7cs6x7qPreOeBd0gflc68l+eRNTFrUF9TRBLfQArrV4CjgUp8Mu7yYcBo51xMp4dQYT2E9LbYBn8wZOdie8YMSE+PbswigyiChbVydidtjW2sWrSKymcryZqexdxn55IxSrMZiUj/9ZSzjzTd3s3AMudcjw1qZnZJf4OTJFRUBGef7ZeQyko/68iqVe3L6tWwe7df/v739sempvrieubM9uXoo2HqVE3/J8lOObuTlIwUZv3fLN44/Q3qVtfxxsI3mPvMXDLGKFeISOT1ucc6XmnEOgG1tUFZWcdie9Uqf1Kbrt63gQBMmeKL7PCCe8YMyM6OevgivTWYPdbxKto5u6m8iZXnrKRuVR1ZU7KY89QctYWISL9EbLq9eKbCOonU1fmDJNet88ubb/pl8+auC24zmDixvdCeNs0vU6f6VhOzqP8IIuFUWEdH875mVp67ktoVtaSNSOOYR46h4MSCqMYgIkPfQFpBROJPTg4cf7xfwh086GchCRXaoaJ740bYssUvf/vb4duaOrVjsR26XlwcvZ9JRAZd2rA05i6dy9pL13LgqQO8sfANZvxyBiMvHxnr0EQkQaiwlsSRleWn8Tv22I7rm5t9+0io2N640S8bNvi5uN94wy+dFRd3LLanTIFJk+Coo/xsJRrpFhlyUgtSmf232Wz8/EZ237ObdVeso+qFKibfMZlAZiDW4YnIEDfgwtrM5gFfAVbgTwRQNtBtikRUWppvATn66MPv27+/vcjesKH9+saN/r5XXvFLZ7m5vsAOX0JF98SJfvpAkTiknA0paSlMu3saubNz2fT/NrHr7l1UvVTFzAdnkjMjJ9bhicgQFokR608D5zjnPmpmV5jZXufcwQhsV2TwFRfDCSf4JZxzsGdPx2J782a/vP02VFe3H0zZlbFjuy66J02CUaMgpa8nPRWJGOVswMwY+9mx5J+Uz9rL1lK3so5lc5cx6duTGPeVcaSk6m9URPpuwAcvmlkaMMI5t9PMUoBc51x1RKLrAx28KFHjnJ+Le/Nm37cdKrhDy9at/iyU3UlPh9JSmDDBj26HX06Y4IvyVHVpJZNoHryonH24luoWNn1xE3t+tQeA3ONymX7vdPLm6UyNInK4iB68aGbnA8XOud+Z2fFArXPuTQDnXBsQ9QQtElVmfqS7uBgWdPF31dIC27d3XXSXlUF5ue/53rSp6+0HAjBu3OFFd+hy/HidIEd6TTn7yFLzU5nxyxmM+PAI1i9eT+3rtSyfv5zRV49m0q2TSB+hvzcR6Z0+jVib2XeALwJPOucuCa77PrDKOffA4ITYO/E0+iHSo7o62Lbf53JcAAAgAElEQVTNj2yXlfkldH3rVn9CnJ6YwZgx3RfepaX+QE4ZMgZrxFo5u+9aaloo+1YZO+/aiWtxBPIDTPjmBMZ+biyBLB3cKCKRHbE+GRgHXBm27gngfiCmSVpkyMjJ6f5gSoCGBj/i3bngDl3fubN9+de/ut7GyJG+wB4/vuvLkSPV550clLP7KDUvlSl3TGH04tG8/ZW32f/4fjZ/fTM7/nsHpdeXMvrq0Zo9RES61dfC+innXKWZhQ9zJ82pcUWiIjPTT/E3dWrX9zc3w44dHYvu8Mvt22HvXr+89lrX20hL8+0mXRXeoesFBZpScOhTzu6nnBk5zPnbHPY9uY8t39hC7Ru1bPr8Jrb/YDul15cy6qpRGsEWkcP0tbDebGbXA2PN7MPApcDFwOcjHpmIdC0tzc8uMmlS1/e3tsKuXb7A3rat42XoekVF+0lzupOb2/Oo97hxmlYw/ilnD9CwRcMoPq+Yir9UUPZfZdStqmPjtRspu6mMMdeOYexnx6oHW0QO6fOsIGZ2InAVUArsBn7vnHs68qH1Tbz264nEpfp6P+rdVdEduqyrO/J2RozoerQ7dDlqlD8YU3o0mLOCKGdHjmtzlD9czrbvb6N2eS0AlmGM+tgoxn5uLLnH5sY4QhGJhp5ydiSm2wsAC51z/xzQhgZoKCZpkbjlHFRWdj/ivX27L8xbWnreTmqqnz6wu1Hv8eOhqCjpW06iPN2ecvYAOeeoer6K7bdvZ9+j+w6tzzshjzGLxzDiQyMI5OgfSpFE1e+DF83sr8CR/gUfDlQAMU3SIhJBZr7gLSo6/BTxIa2tvo+7p+J77972Ay+7k5PTc+E9frxmOekl5ezoMDMKzyik8IxC6tfXs+OuHez97V5qXq1h/avr2fTlTYz8yEhGXTmKvOPzsCT/x1EkmRypx3pV8DHrgTbgGvxpcMOHGSYBxYMSnYjEr0DAT/s3ZgyceGLXj2lo8LOXdC6+wy9rauCtt/zSneHDey68R4/WSXU85ewoy56ezbSfTmPyDyfzzh/fYfeS3VS/XM2uu3ex6+5dZE7OZORHRjLyipFkT8+OdbgiMsh6bAUxs4LgYyqDt292zt3UxeMecc5dOHhhHtlQ/lpRJKlVVfVceO/Y4WdC6UmoyO+u+J4wIa5bTiLVCqKcHR9qV9ey55d7eOfBd2ja3XRofe78XEZ8eAQlF5eQNVnfwogMVf1uBXHOVXVaNdLMAs65Q+drNrOFwPEDjlJEklNBAcye7ZeutLX5lpLuCu/t22HPnvY2lJ5eZ/JkOOoov4RfLy1NiBFv5ez4kDs7lyl3TGHyjyZT+Wwle3+3l/KHyqldXkvt8lo2f20zObNzGH7xcIZfPJzcY3PVLiKSIPr6SfJnYKWZvQjUAEcD5wGHjYiIiERESopv9Rg9Go7vph5sbOy55aSszI+Mv/66XzoLBPyodqjQnjIFZs70y4QJQ/lkOsrZMWQBo+jsIorOLmLqz6ay//H9lD9Uzr6/7aNudR11q+vYevNWMidlMvzC4RS/p5jC0wtJyRiy7zeRpNef6fbGANcCc4Bq4E/OuUcGIbY+SeSvFUVkgJzzc3dv3gxvv+0vQ8vbb/uivLtcmJXlz5IZKrTnzIETTvB93xEyyNPtKWfHmbbGNg4sPUDF/1VQ8UgFzXvbW51SslMoOquI4kXFFJ9fTNZRahkRiTeDPd1eCjDPObd8QBtq396JwKnAJuBFoA74IX60pQi4zjlX3fl5yZykRWSAGhr8zCWhQnvDBli3Dtauhd27u37OUUf5AnvhQrjgAj9ndz9Febo95ew44lod1a9UU/FoBfuf2E/dqo7zx2dNzaJ4UTFF5xRRcHoBaYVpMYpUREIiVlib2X1drC4BWp1zF/UzvvDt/wcw2Tn3zbB1vwSedM79wcwWAVc65y7v/FwlaREZFAcO+CL7zTd9ob18uV/q69sfYwZnnw3f/z7Mn9/nlxiswlo5e+hp3NXI/if3++Uf+2mtam2/MwXyjsuj8KxCis4qouDUAs2XLRIDkSysX8DPfRr+pFOAF5xztwwwyDOAG4FzXTCo4FeYW4Ai51y9maUCVcAs51xZ+PMTKkm3tsKKFX707OBBP5vB6NH+AKsRI+J2ZgORpNHSAmvWwMsvwxNPwGOP+VaSnBzf113ct9nsBrGwVs4ewtpa2qh5tYb9f99P5dJKql+pxrW0/yot1cg7IY+iM4soPLOQ/BPzCWSr0JbIc87RtKuJpr1NNL3TRMv+Flybfy+mZKSQPjKd9FHpZE7KJCUt8Y8R6PesIF242jm3vtPGZ+BPlztQdwAvAXeZ2RTgFmACUOGcqwdwzrWY2RbgDKAsAq8Zf1auhMsu819Fd6WkxPd4zpnjT9wxZ47v/8zMjG6cIsmqqcn/07t9O+zf79eFBigOHvRtJfFDOXsIS0lNoeCUAgpOKQCgpbaF6n9Vc2DpASqfqaRmeQ3V/6qm+l/VbL11K5Zq5M7PpfC0QgpO889LG6bWEekb1+aoW1PHgWcOULOshvp19dSvr6etru2Iz7V0I2dWDvkn5FO8qJjCswpJzRv6My71RSR6rCcBrznn+n0kj5lNB9YBc51zq8zsa/ij1m8FPuKcmxP22JeBvzrnvhe+jQkTJriSkpJDtxcvXszixYv7G1JstLTA1Kl+BoNx43z/Zmam//Devdt/mFcf1qroZzSYOROOO84v8+bB3LmQlxf1H0FkyGpq8n9r+/b5Ax137fJ/d7t2tS/btvle7LYuPmDmzoVbb4X3vrdXL7dkyRKWLFkCwPLly7c65yZG8KfplnJ24mipaqHy+Uoql1ZS+WwltStr/WmBwmTPzKbgtAJfbJ9aQOYEDcJIR845Dm44yIFn/D9sB5YeoGVfy2GPSytJI31MOukj0kktTsVS/bfnbfVtNO1tonFnI41bGzs8xzKMkotLGPWJURSdXYSlDN1v3Hubs/vaCrKFjl8pBoCRwP3OuWv6G6yZfQBYEkr0ZpYFlAM3Ax90zp0Q9tjXgQecc7eFbyMhvlZ86SU45RR/UNTatYePQjvnP9hXrfIj26tW+WXjxsM/6M18kR4qtEOXw4ZF7+cRiYTVq32Pc1GRv+1c75bQYxsb/VR7VVVQWekv9+1rL6JDS01N7+JJSfFT8E2dCjNm+CkATz4ZJk3q9484iK0gytlJpKW6heqXq6l8oZKqF6qofrUa19jxMz5jfIYfzT6tgIJTC8iZmTOkix3pn4NlB30R/cwBKpdW0rSrqcP9GeMyKDzL/zOWMyuH7BnZpBUf+duPluoWalfVUrm0kv1P7Kf6lepDGSh7ZjYTbphAyWUlpKQO7XaRSLaCPAU8QHuibgP2Oue66VvotdTwWJxzB81sI/5DoKjTY/OAXQN8vfi0K/hjzZnTdWuHmf9AnzAB3v/+9vV1db7AXrGifZ7eNWt8O8mGDfDgg+2PHT8eZs2CY45pX44+GrJ1ql2JQ6tX+7+HaAgEfG/0sGF+CZ2uPXwZOxYmToSMjOjENHDK2UkkNT+V4vOKKT7P9/i3NbZRs7yGqheqqHrRL43bG3nngXd454F3/HOKUw+1m+SflE/e/DwdEJmAGnc1Urm0vZBu2NKxZS2tJO3QQbGFZxWSNTmrXyctSs1PpfDUQgpPLWTijRNp2NbAnl/vYfeS3dS/Wc+6j6yj7OYyptw5hWGLEnOgr6+F9W+dc8+HrzCz6WY2xTm3aQBxrAIKzGykc25vcF0LPhmPMbOsYOJOA0qB5wbwWvErVNwePNi35+XkwEkn+SWkqcmPeocK7ddf96PcobPTPflk+2PN/Cj5Mcf4ovvoo/1o3NSpfT4ISySiSkt9IVtWdvh9U6b4syeadb8ApKf7sy4WFvrLgoKOBXRoyc8fyieC6Y5ydhJLyUih4OQCCk4ugOuCvbNr69oL7ReqaNzRyL5H97Hv0X3+SQHInZNL/on55J+UT/6J+WRN6V+RJbHTVNFE5bPBNqFnKql/q77D/amFqRQuLKTwzEIKzyokZ1bOoPyOM0szmXjjREqvK2Xv/XvZ+r2tHFx/kNXnr2bY+4Yx9adTE649qa+tIN91zl3faV0hsNQ5N29AgZg9DvzNOfez4DafA+YDv8V/ODxmZucDH3POXdH5+QnxtWJodG7q1O4PXhyI1lY/R+/atX5EO7Rs2OD7u7sybFh7kT1tWvv1KVN8ISISDfv2wV//Cn/4A/z9735ddjbcey9cfthMbkPOILaCKGdLt5xzNGxt8G0jL1dT/Uo1tatqobXj41KLU8mdl0vu3Fzy5uWROy+X7OnZWEDFdrxoqmii6vkqX0w/W0nd6o7zoafkpFB4WuGhUencubkx+f21NbWx4yc72HrzVlprWgnkBZjy4ymMumrUkPrnbcDT7ZnZ5/AHpXR3NNxTzrlF/Q8RzGw48GNgGX6E4xfOuTfNrAh/soEyYBTwDedcbefnJ0SSbmjwo8/g58iN1tfNTU2+uA4vtDdu9EtdXffPKyxsb03paikp0dSAEnnr18N//ZdvcUpP92dNjOBZEGMh0oW1crb0V2tdKzXLaqh+xRfaVS9XdTgzZEhKVgo5s3MOFdw5M3PIPjqb9JL0GESdXJxzNG5vpPrV6kPFdN2ajp/VlmEUnFxwqJDOe1deXE2D17inkY2f3UjFwxUADLtgGDN+OaNXfdzxICLzWJvZHOBK4Ced7mpyznVzarLoSZgkPXUqbNrkR6+POSa2sTjnZ0UIFdnhBXdoju2eZGX52U1Gj/b9qaNHd309P18FuPROayu8+CL86ld+Adixw/c+D2GDMWKtnC2R4JyjcUcjtStqqX2jltoVtdSsqDls9oeQtOFpZB+dTfbM7EPFdvaMbDLGZuggyX5qqWrx/+y8Wk31q9XU/LuGpj0dDzZMyUwh/6R8Cs8opHBhIXkn5BHIjO9eeecce3+3l42f20hrVSuZkzKZ9dAs8ubF/4xmkTxBTJpz7rB/Xc2s1Dm3bQAxDljCJOkLL/RfeT/4IHzoQ7GOpnvO+SnJtm5tX8rKOt6urOzdtrKyfJE9cqQfeSwp8ZehJfx2SYkK8WRRW+tHp0NnPfz3v+HVV/168Acbfu978LWvxTbOCBjEVhDlbBkUzQeaDxXadavrqFtXR/2b9bTWtHb5eMswsiZlkTk5k6yjssia3H49c1Imgaz4LgKjoa25jfr19X5/Bpfa1bVd/hOTWpRK3vF5FJxSQOHCQvKPzyclI35GpPuiYWsDay9ZS82yGlIyU5h2zzRGfXxUrMPqUb9nBTGzrwKPO+feDK76upl1fvcbcCawcKCBCv7gwb/+1fdBxzMzX+SWlMCCbuqB6mr/NX1oPuDQnMCdr9fV+RHwzZt799qpqb7ILirqeFBaby7z8ny7TWpyTVgfN5zz74uKivb5oisq/Ptg+3Y/+hy6fOedrrcxdSp88IPwyU/6Xn85RDlboiWtKI2iM4soOrN9EhjnHI07G/0JRd6sP1Rs16+vp/mdZurfqj/sILpD2xuRRsbYDDLGhS1ht9NHpRPIDwypPtzOnHO0VrfSuKORg28f9Msmf9nwdgMNZQ0dzqwZYhlG7txc8o/PJ/+EfPKOz0uoA0ozJ2Qy94W5bPrCJnb/YjdvXfkWBzcdZOK3Jw7Jn/FI1cUkoDDs9nFADrCH9umbDN9HJ5Ewa5a/jPfCujfy8/1y9NE9P66mxhfZ5eXthVZFRfe3a2pgzx6/9FdGhi+wQ0tubs+3c3L8FIgZGe2XPV3vfDsQGFqj7M3Nvue/N0t9vf+dhC/V1YffDs0b3d2Bsp2lp/sDZmfM8Mtxx/mZb0Yp3fRAOVtixszIHJdJ5rhMis/tOKNUS20LDZsbfBEZvAxdbyhroPmdZprfaaZ2xWHt+O3bTzfShqf5E5WUpJNWkta+DE8jNT+VQH6g42WevxyM0dy25jba6ttoPtBMy4EWWva3dLxe0exPmrKrkaZd/gQqbfU9nL3QIHNyJrmzc8mZnXNoyZqSNeTnfT6SQGaA6Uumk3tcLhs/t5Gtt2yloayB6fdOJyV9aP3sPRbWzrnPdlr1HWCNc65Dc4+ZzYx0YEkrkQrr3srLg+nT/dIbjY2+wA6d7CP8sqt14ffV1PgR8sZGv4ROSR0NgYAfKQ9dhl/vaV1KSv+L8rY2X8i2tvrLzkt36wd4RtYe5eX52WbC231GjPBzrI8f7/vyx43zBXRAXw/3hXK2xKvU3FRy5+SSOyf3sPtcq2s/c9+OrpemvU201bXRtKuJpl1N1NHDgfVdsDQjJSMFS/eXna8TwP/rGVycc+3X2xxtB9sOLa0HW2lraDts9pTeSMlJIWNshm+BmZxJ1hTfFpM1JdgSE+d90YNt7KfHklmaydrL1rL3/r007mzkmL8cM6ROi97XSGc7514PXxE8te1k4M2unyJ9Mn26L6Q2bfKF39A5EUX0ZGT4g9X6e8Cac36ktbbWF9ldLV3d19jonxcqyru63t39bW2+iG3tRyaOhZQU3/uemXn40tX6vDy/5Oe3X+98OzR3tN7T0aScLXHPAkbGmAwyxmTAu7p/XOvBVprLm/1S0UxTeVP77X3NtFa30lLT4i+rwy6rWnHNjtbmCOffAASyA6QWppJanEpaURqpRWHXh6X6dpYxGaSPTSdjTAap+UOnQIyVYe8Zxrzn5rH6faupfKaSVe9exewnZpNWODRmDOnVb9jMTsN/fXhu8BS54YrxR50/GuHYklNWlj+Ib/du3/owblysI0o8Zn4/Z2X5HvFoCC+sO48UH2ndQIpxs/aR8PAlfIS88/rQfTJkKWdLIgpkBQiUBsgs7dsJRZxzuCZHW2MbbU1tuMZO15vafG+z+XYWjEOLmUGKn3UjJSuFQFaAlCx/PZ6mr0s0efPzmPfiPN446w2qX6lm5VkrmfOPOaQPj//pHHv76VkJ/Abfv3dKp/uagHsjGVTSy8rylw0NPT9Oho6UFL+kDY3/uGXIU84WCTIzLMOG7KwZySprchbzXvDFde2KWlaeuZJjnz6W9JHxXVz3qrB2zq02s4XApc45JeTB1hicWic9vt88IhKflLNFJBFklmYy7/l5rDx7JXVr6lh5zkrmPjuXtGHxO0jV63/fnHNV3SVoMzs1ciEludD80DDkzyYnIrGjnC0iiSBjTAZzn5tL9oxs6tbUsWrRKlqqezm7VAwcaR7rvwKHH8Lb6WHALGBEpIJKauXlfsS6oACys2MdjYgMIcrZIpKI0kekc+zTx7LitBXULKth9ftWM+fJOQSy428WlSO1gqzH9+rtpH0O1M4MeF8kg0pqa9b4y5maDUtE+kw5W0QSUsbYDI79py+uq16oYs3Fa5j96Oy4m+f6SIX1D4Aa59zh59MMY2Y6ujxSXn3VX86eHds4RGQoUs4WkYSVNSmLuf+cy4rTV3DgHwdYf/V6Zvx6RlydobHHMt85V3GkBB10cYTikf/7P3+5aFFs4xCRIUc5W0QSXfb0bOY8PoeUnBT23r+XLTd2nlE0tgY0fm5mpWb2A+DHEYonub38Mrz2GhQWwrvfHetoRCTBKGeLSCLIm5/HrD/OggBs+842di3ZFeuQDulzYW1mATO72MweB94GPoLv6ZOBcA5uuMFf/8xnICcntvGISEJQzhaRRDTsPcOYdvc0ADZcu4F9f9sX44i8XhfWZjbJzL4DbAceBKYA7wfGA5cOTnhJ5P77YelSf8rnr3wl1tGIyBCnnC0iiW7MNWMovaEUWuHND79J3Zt1sQ7pyIW1mV1qZv8ANgIfBO4AxgF/c8496byXBjnOxLZuHXz2s/76bbf54lpEpB+Us0UkmUy6ZRIll5XQWtvKmovW0FzZHNN4ejNifRB/CtyXgAucc7c558rpfion6YuKCrj4YqithQ9/GK68MtYRicjQppwtIknDzJhx3wxy5uRwcONB1l2xDtcau3R3xMLaOfeYc+59wEeBj5nZ/WZ2CXDofJIWT/OcDCU1NXD++bB+PRxzDPziF6BdKSIDoJwtIskmkBPgmL8cQ+qwVPY/sT+mM4X05ZTm25xzNwL/EVw1zsxuMbOJwDWDEFtiKy+Hs8+GZcvgqKPgH/+A3COdME1EpHeUs0UkmWRNymLWH4IzhXxvG/ueiM3BjH2eFcQ51+Kc+7Nz7mLg18CXgTsjHlkie/ttOPVUP7XexInw1FMwenSsoxKRBKScLSLJoujsIibdMgmAt658i8bdvZnWP7IGNI+1c26Tc+6LwKcjFE/ie+QRmD8fNmyAY4+Fl17yI9YiIoNMOVtEEl3pdaUUnl1Ic3kz6z62DtcW3X7riJxg3Tn3m0hsx8zSzWylmS0M3s4ys7vM7LtmdreZ5UfidWLi4EH4f/8PLroIqqr8AYvPP6+RahGJukjlbEjwvC0iQ46lGEfffzRpJWlU/rOS7bdtj+rrR6SwjqCvARPDbv8P8KJz7nrgEeCeWAQ1YK+8AvPmwX//NwQCfkq9hx6CfH3eiMiQl5h5W0SGrIzRGcz49QwAtty0hfr19VF77bgprM3sZGA3cCB4ewxwBfBo8CFPAxcED7wZGurq4Otfh1NO8TN/zJzpT1v+la9o9g8RGfISMm+LSEIYdv4wRl01CtfoWH/N+qi1hMRFYW1mOcAlzrn7wlYvBCqcc/XgD8ABtgBnRD/CfnjkEV9I/+hH/vbXvw7Ll8O73hXbuEREIiAh87aIJJTJt08mbWQaVS9UsWvJrqi8ZlwU1sB1wPc7rRsLdJ4rpQYY09UGysvLWbBgwaFlyZIlgxBmL5SVwQUX+F7qbdtg7lx/gOIPfgCZmbGJSUTizpIlSw7lK2B4rOPphwHl7bjJ2SKSsNKK05h611QAtnxzy4DOytjbnJ3a71eIEDM7D1jmnHun010OfwaxcBlAl3ulpKSEZcuWDUKEvdTYCLffDrfe6g9UzMvz16+9FlJjvptFJM4sXryYxYsXA2BmFTEOp08ikbdjnrNFJCmUXFJCwekFVD1fxbbvbWPyDyb3azu9zdnxMGL9VeA+M6sIBjoef8BLDlDU6bF5QHTG8vvi6adhzhy44QZfVH/oQ/DWW/CFL6ioFpFENPTztogkBTNj8u2+mN5x5w4atjYM6uvFQ2F9BTA3bNkFXA38EhhjZlkAZpYGlALPxSjOw+3YAZddBuee6+elnjbNn0HxwQdhTJcdKyIiiWDo5m0RSTr5C/IZ8eERuCbH9jsGd/q9mBfWzrly59yO0AK0AuXOuW3AY8DZwYeeAzzknNsZq1gPaW72U+bNmAF/+hNkZcF3vwurVvkiW0QkgQ3JvC0iSa30G6UA7L53N837+t9rfSTx3qfwGeCHZnYsMApYHON44Nln4bOfhTff9LcvvhjuvBNKS2MalohInIi/vC0iSS93Ti7Fi4rZ/+R+dv1iFxP+c8KgvE7cFdbOuYlh1w8A18QumjB79sBXvwq/+52/PXky3HUXnH9+bOMSEYmxuM3bIiJhxn5uLPuf3M/e3+yl9LpSbBDOKRLzVpAh4S9/gVmzfFGdmQnf/jasWaOiWkRERGSIKHp3EWkladSvq6f29dpBeQ0V1j2pr4dPf9q3e+zfD+9+N6xdCzfdpDmpRURERIaQlLQURnxoBAAVfxmcWU5VWHdn1y44/XS45x5IT/d91E88AUcdFevIRERERKQfis8vBmD/U/sHZftx12MdF1atgve8B3bu9IX0ww/DscfGOioRERERGYDCMwqxNKPmtRpaqlpILYhsKawR687WrIGzz/ZF9amnwquvqqgWERERSQCBnAC5c3OhDWrfiHyftQrrcNu2wTnnQEWFH7F++mkY3u3p4EVERERkiMmdmwuosB5cjY1wySWwdy+cdRY89BBkZMQ6KhERERGJoJw5OQDUrlZhPXhuvhleew0mTvRnU9SsHyIiIiIJJ+uoLAAatzZGfNsqrAHWr4cf/chf/93voLg4tvGIiIiIyKDIGO87Ehp3qLAeHLfeCs3N8IlPwMknxzoaERERERkkocK6YXtDxLetwnr7dvj97yEQgBtvjHU0IiIiIjKIUgtSsVSjra6Ntsa2iG5bhfWDD0Jrqz+74sSJsY5GRERERAaRmREoCADQUtUS0W2rsP7zn/3l5ZfHNg4RERERiYrQiWFUWEdSVZWfCSQtDRYtinU0IiIiIhIFocK6tbo1ottN7sL6lVfAOTjuOMjOjnU0IiIiIhIFgVzfCtJaq8I6cl5/3V+eeGJs4xARERGRqLE0A6CtWQcvRs6GDf5yxozYxiEiIiIiUZOS7ktg1+Qiu92Ibm2o2bjRX06bFts4RERERCRqQiPWrlmFdeTs2OEvJ0yIbRwiIiIiEjWWHmwFaVIrSOSUl/vLESNiG4eIiIiIRE1KWrAVRCPWEVJf75f0dMjNjXU0IiIiIhIlCT9ibWYfMLMtZrbPzH5sZqnB9cPN7B4z+4GZ3WFmaRF5wX37/OXw4WAWkU2KiCSLqOdsEZEIskCwx7olAUeszawUuAi4BPgycDXwpeDdfwLuc85dB2wAbo7Ii9bU+Mv8/IhsTkQkWcQkZ4uIRFKoAo5sXU1qZDfXbxOATzjnWoDlZnYMcKaZvQhMd869Gnzco8B6M7vVOVc3oFc8eNBfZmUNaDMiIkko+jlbRCSCLNStEOHCOi5GrJ1zLwQTdMhOYBtwFrAl7HE7g1cXDPhF6+v9pc64KCLSJzHJ2SIikRTqAk7EwroL7wJ+DowF9nW6rwYYM+BXCI1Yq7AWERmowc/ZIiKRFBqwdpGtrOOlFeQQM5sClDvnVpqZAw52ekgG0Nz5eeXl5SxY0D4osnjxYhYvXtz9C4VGrNUKIiIxsGTJEpYsWRK6OTyWsQxE1HK2iEgk9XHEurc5O64KazMLAJ8C/jO4aicwLex+A3KAXZ2fW1JSwrJly3r/Yg0N/jIjo7/hioj0W3ghaWYVMQ6nX6Kas0VEIqivPda9zdnx1gryVQlpc30AAA95SURBVOA251xj8PZTwPSw+ycCTcDAs3Frq79Mjav/LUREhpLo5WwRkUhK9B5rM7seWAFkmdlRZvYJoBLYHTziHGARcJdzrmnALxgqrAOBAW9KRCTZRD1ni4hEUiL3WJvZN4FbOq1+yzl3n5l9CLjJzMqAXOCGiLyoCmsRkX6JSc4WEYmkQRqxjovC2jl3K3BrN/dtBT4Z8RdVYS0i0i8xydkiIhGU0PNYx4QKaxEREZHkNEhnXlRhrcJaREREJLmEBqzbIltZq7BWYS0iIiKSVNQKEmkqrEVERESSU6JPtxd1KqxFREREktMgTbeXvIV1aEeGvgoQERERkeSgEesIU2EtIiIikpRskOo/FdYqrEVEREQkAlRYq7AWERERSU5qBYkQFdYiIiIiyWmQyr/kLaxzc2H0aMjPj3UkIiIiIhJFgfwA6aPTCeREdna41IhubSj53Of8IiIiIiJJZeI3JzLxmxMjvt3kHbEWEREREYkgFdYiIiIiIhGgwlpEREREJAJUWIuIiIiIREDSF9ZLliyJdQhDmvbfwGkfDpz2YfLQ77p/tN/6Tvusf5J9v6mwTvI3wEBp/w2c9uHAaR8mD/2u+0f7re+0z/on2fdb0hfWIiIiIiKRYM5F+FyOMWJm5cDWfjx1OFAR4XCSifbfwGkfDtxQ34cTnHMlsQ4impSzo077re+0z/onGfZbtzk7YQprEREREZFYUiuIiIiIiEgEqLAWEREREYkAFdYiIiIiIhGQtIW1mWWZ2V1m9l0zu9vM8mMd01BgZreamQsuK4PrtC+PwMzONbNXzWxi2Lpu95v26eG62ofB9Ye9J4PrtQ8TiH6f/dfd34h01Nc8LV5fc3OiS411ADH0P8CTzrk/mNki4B7g8hjHFNfMLAd/tO+5wVWhI/q1L3tgZiOBXOD4Tnf1tN+0T8N0tw97eE+C9mGi0e+zH47wNyJB/czTSa+fuTmhJeWsIGY2BtgCFDnn6s0sFagCZjnnymIaXBwzsy8Fr/7cOdcQXKd92QtmlgK0ApOcc2U97Tegqbv7knmfdt6HwXWHvSeD6/W+TCD6ffZfd38jcri+5Gm979r1JTcng2RtBVkIVDjn6gGccy34P54zYhlUPDMzAz4M3A7sNbPQf+wL0b48IudcW6dVC+l+v/V0X9LqvA97eE+C9mGiWYh+n312hL8R6aSPeVqC+pibE16ytoKMBfZ1WlcDjIlBLEOC819tnGhmRcBXgN+ZWSXal/3V035r6eE+CeruPemcewK9LxONfp/9cIS/ETkyve/6Idnfd8laWDvgYKd1GUBzDGIZUpxzB4BvBv8j/RLwFNqX/XGk96D2aS918Z58Av2NJxr9Pgegm78ROTK97wYgWd93ydoKshMo6rQuD9gVg1iGqp8ApWhf9ldP+037tH9C70nQPkw0+n1GRvjfiByZ3neRkVTvu2QtrJcCY8wsC8DM0vC/9OdiGtXQ0ga8jvZlf/W037RP+yf0ngTtw0Sj32dkhP+NyJHpfRcZSfW+S8rC2jm3B3gMODu46hzgIefczthFFd/MbIyZXWFmKcGvdb4C3KB92TvBfQZg0PN7UPu0a533YXfvSdDfeKLR77N/evobkcP1JU/HILy41ZfcnAyScro9gGBT/Q+BMmAU8A3nXG1Mg4pjZjYL3x9VDbwI3Omceyt4n/ZlD8wsD/gofj7Um4GfOufKe9pv2qcddbUPgRF0854MPkf7MIHo99l3PeVt6ag/eVr6l5sTXdIW1iIiIiIikZSUrSAiIiIiIpGmwlpEREREJAJUWIuIiIiIRIAKaxERERGRCFBhLSIiIiISASqsRUREREQiQIW1iIiIiEgEqLAWkbhiZiVm9lcz22pm/xnreEREpGfK2+10ghgRiStmdh6wFMgH3gKmOef2xzYqERHpjvJ2OxXWIhK3zOx54CznXEusYxERkSNL9rytVhARiUtmNh74R7ImZxGRoUZ5W4W1DAIz+7qZ/SzWcYREIh4zO8fMXjazhREKK6GY2W0R3p4BlwHfj+R2RaQj5euBM7OJZvYLM7snQtsLmNk1Zra1F499r5m92t+fNZK5W3nbU2E9xJjZ0Wb2sJnda2YrzcyZ2UdjHVcnjwO/j3UQYSIRz3bg2AjEEpfM7D1mVmFmu83sJDMbH3yfOTP7qZkVBh93vJmtMbPbzSwruG4ksCfCIV0B/K9zriW4fZEhR/m6X4Zivn4neJkRoe2lAJVAaS8e+//bO/uYrasyjn++IpRIgKiMV40aiUOdTZSamUZYSwzSeKk5GjQJZxYu51oNE5FNW9AstRVOR61EnBmRGS21ZYoFT3M5SZBFiC/NCkMwlBe5+uM6N5zn9v7dLw/3zcPNc322M36/33m5rnPO776e8zvnOocO4PSiSEljJa1I7958SUOyuGbb7rDbxMC6HVkJ3GRmVwIfBH5SipB0lqQLu0OpXLaZPWtmT3SHHpVoRJ+iNjSzjcB/mqlXM/vrUMsys4eBHwD7zOwpM3sRmAPsAf5sZttTurW4Ib/ezN5M2T8NPHRIFciQdAOwAHhE0gbgtGaVHQSHmbDXDdKO9trMdgGvNEuWme0F/lJn2lfxQXhR/HpgabpdbGb5QLpptjvs9kFiYN1GSBoEfAD4H4CZ7Qe+BZikAbjRVjfo1W2ym0kd9dh/GGV1R1krgBGSzgEws23Ao8ClmayBwCvp3Ssxxsw2HKLsA5jZzWY22szGmdkYM3u8WWUHweEi7HVrOQLtdbNPgmhE/1pp307/lvs9N812h90+SAys24h0dM1GYLWkD6VnLwBrcb+mUcCXJM0GkLRQ0g1pGeiAH5WkcZK+K+laSbskbZT0NUmfkbRO0mWSfiFph6S5Wb6K5ZXJniNpWe5rlpZDb5d0i6TfSzovPa8qL8svSXMlrZd0saRnJL0qaVYdMkbn+tSQ+Y42rMAJaWl3p6Q5ZW26ILXNg8kQ19tm75AlaaqktZJmSdok6dYq5VXq+4r6VCPNbKwHpmWPTwYmKbl9AJOBVZmexwG7svsxkm6VtFjSH+QuJTMkdUiaImm5pO2S5kk6J9Vxm6RLaukXBO1E2OueY68rUVTHanIkjZJ0j3z29zv1yKkg92FJP5N0fo10YbtbhZlFaKMAjAVewL9Q7wEGZ3FbgIvS9RnA5nR9Mv41fSL+xb0VGJfi7gTWpOs+wGvAHfhZlJ8HXq5WXrls3MfsR8Cy9Pzd+GBtYLqfCvwbeE81eRXqPSTJ/GZKuwh4Czi1hoxyfarKzNuwgg5bgCVA35Rva3reH7gvS/ckbhTrarMCWYPw5b3bgPOAidXKK+v7ivrU+X7NB/6ert8P3ALsBi7P3hdl6ScDF6TrvsA6oG+6fzy1Q3/gdWAhcBwwI91Pxz/uryO9gxEiHE2BsNc9wl6n+AV1tmO1vl4HnJ7lsTrfs1J/fhi4sUL8RUnOsdmzsN0tCjFj3WaYzyqeiRuaLwDrJVXapLEJmCmpN3BBetYP/yGPBEr+sc/hPxjMbA+wE3jAzHbgP7QhNcor1283nTdDTALetOSjCzyI/yAvryGvvNxSmfemtItw/9+La8jopE8jMgv4lbk/3TpgaFbHE9OM0rV4m0KdbVZQ39dwA/ZLM1trZo80UF6RPvWwAnif3B3ks8BdwG+B6ZL6ATstWd7E+cCadD0Z/+NVmgW5DLg5tfN24DFzv+wOoL+Z3W++PP40jfVBELQFYa97hr2uQGEdq8iZCAwws5I+HQ3KPBe4wsxuqjN92O4WcWx3KxDUj6RewEgz2wJ8VdLd+A7qpcD4PK2Z7ZY0HPg68ONSEWb2L0kdwMfxL+pR+GDqQNbseh/JXaiovAJV8zJGA70zvfZL2gKMqCavFmb2lqTNwEDgpAZkdFlmGfs4+Ps5BfiHmd1WnqieNpM0Ex/AljjNfMnYcl0b6INCfWphZpskPY27gwwws82SVuCzSNOAX2d6l96Nkv/eqfgMU6msbXnR2fXbdGY/4ZYWHGWEve5Uv55gr3MK27FK35yFf0R0lQnAaEn9zOyNgjT7Ux3CdreQaJD2ojeZ/6uZ/RVfahtTnlDSBOBqM1tkfsJDzpeBc+XHPj0DfLuW4BrlVWMrPgOaH0Mk3PfwUOkDPN9iGfXwTzr7ISNpfANttgo4OwsVd5c3UF5FfeqvDvcBs4HNmX69gHn4smmJ8bi/aImXgY9IOj6TW9XPLwiOYsJed6Yn2evCOlaRsyPl6UPXWIzP9n+vIP51O7jpPGx3C4mBdfvxFUnDsvvhwOp0vQffrDEGP9ppgKR3SZqY4gdLOgn39bsbX8ZZQ1paTBxD2Ve6JNUor1y2sjJWAm8An0tlDUzySjOfRfKKGJrSDE95V9chI9enlsy8HuVUmr0QPgt1PPCQpE/KN8nU6oMDsoChZrYhC3szEb2y62rl5XpX1EdSH0lLk5tHNe7Hl6AfADCzncBvgCet82kgn8DdREqU2nu5/Czsa4DSWablfVDez219QkEQFBD2mh5jr+ttx4pycFvaB//4Al+dQNKwOmy38NnjmcA0SdNT3uskzcdn6x/N0oftbiWtcNyO0JqAb4gw4L/AT4Fl+FLSCSn+RuAlYAq+tLMJnyGYhM90LMON1Cp8I8WeVN4u/If2KWAvvkFmJPCNFP/FauWVyZ4H/BH3Wxub4sbhM51LcF/DM9LzQnkF9beUfx7umnBmFlck4725PrVk5m1YJnsKvonvjtQWpXyzU/yFwLOpb5akdq63zaZUqOvUJG85MDw9q9anncoq0GcE/h8Z3FvHu/b9svsZwMfKni2ukO+jwN+AbcD89OwS/F27HRhGOnIMuAr/w3tXquul3f0bixChWYGw1z3JXp8CPIHPup9do47V+noSvlL4J3wj+fPA1anuFW03viqyF5+pHgQ8hg/qrwTm4ud5L8HdUEp5wna3MCg1XtBDkDQYuMrMFqb7Y/BZlLlmNr9blauBJANGmfssBl1E0vVm1qWjnIIgOHyEvQ5ywna3B+EK0vO4BncLGJQ9m0C2Ke1IJFt6imWnLiKpl6QrgJ93ty5BENRF2OsgbHebEQPrnscPcV+s5yS9CPwOeMnMnupetYpJmznmpdtZ8mPfgsbpB6w0s801UwZBcCQQ9jqAsN1tRbiCBEEQBEEQBEETiBnrIAiCIAiCIGgCMbAOgiAIgiAIgiYQA+sgCIIgCIIgaAIxsA6CIAiCIAiCJhAD6yAIgiAIgiBoAjGwDoIgCIIgCIImEAPrIAiCIAiCIGgCMbAOgiAIgiAIgiYQA+sgCIIgCIIgaAL/BxHk6nac6iVDAAAAAElFTkSuQmCC\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/apollo-6-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/apollo-6-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/apollo-6-earth.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }