{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 21 - Viking-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 Viking-2 Lander.\n", "https://en.wikipedia.org/wiki/Viking_2" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"MARS\") \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('Viking-2', 982, 63.6, 0.18, 9.65, 0.0, 0.88, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,4.48,0.0,-17.08,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,20.0)\n", "# bank angle = 20 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": "iVBORw0KGgoAAAANSUhEUgAAAtcAAAHtCAYAAADSh1ILAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZxcZZX/8c9JJ93ZQ/Y96ZBAQAKB0GAAmYRtAAFBBsQAKqhphYFR1IgDiI7gj0VlYEAcW2UQkBExKDsIShh2SYIhAUIIZCNrZ09I6KTT5/fHc6u7utJ7V9et5ft+vZ5X3bp1695zq9OnT5567nPN3RERERERkfbrFHcAIiIiIiL5QsW1iIiIiEiaqLgWEREREUkTFdciIiIiImmi4lpEREREJE1UXIuIiIiIpEnnuANIlwEDBnhpaWncYWREZWUlAwcOjDuM2OlzCPQ5BLn8OcyZM2e9u+dm8G3U1pydyz/nltI55gedY35o6Byby9l5U1yXlpYye/bsuMPIiLKysoI516bocwj0OQS5/DmY2bK4Y8i0tubsXP45t5TOMT/oHPNDQ+fYXM7WsBARERERkTRRcS0iIiIikiYqrnNQeXl53CFkBX0OgT6HQJ9DYSiEn7POMT/oHPNDW87R3L0DQmnigGYnAdcD57n70mjd2cDPgN7AfcC33b06em0A8GNgM9AFuNLdd6fut6yszPN93I+I5Cczm+PuZXHH0RDlbBGR+prL2RntuTazwUBP4MikdaOAs4BzgCuArwLfTHrbg8Bd7n4lsAj4UcYCFhEpYMrZIiKtl9Hi2t3XAg+nrB4NfNnd57j7PcDPgeMAzGwyMN7dX4u2fRS43Mx6ZCpmEZGmuDu7du0i098CZkI25Wx3qK4OTUQkm2V8zLW716Q8fyHxdWJkJbA8Wj4eWJK07cpoMSu/PhWRwjNy5EhKSkpYtWpV3KF0iGzJ2S+9BF26wJQp7d2TiEjHysYLGo8A/jtaHg5sSHl9GzAs9U2VlZWUlZXVtoqKig4OU0QEunTpAkBVVVWr3ldRUVGbr4AB6Y8sYzKSs83CYx5+QSAiOaA1OTurbiJjZuOASnefF61yYGfKZiXAXhfHDBw4MO8nMheR7NO1a1cAPv7441a9r7y8vPYqdDNbn/bAMiCTOVvFtYjEqTU5O2t6rs2sCPga8L2k1SuBvknbGNADyM/vX0Uk5ySK69b2XOe6uHK2imsRyXZZU1wD3wF+6u7Jf6GeAcYnPS8FdgHqohaRrFBSUgK0vuc6D2Q0Z6vnWkRyRcaL66gnA8CS1l0FvAF0M7N9zezLZjbO3WcDq81sQrTpKcDt7r4rs1GLiDSsrcNCckW25GwV1yKSKzI65trMegEXRk8vMrM7CF8rXpey6UJ3vytaPg+41syWEuZbvToTsYqItESi5zofh4VkU86uLfFFRLJcRotrd98G/CJqCddHrbH3LAO+0sGhiYi0ST73XGdTzlbPtYjkimwacy0iknPyuec6m6i4FpFcoeJaRKQd8rnnOpuouBaRXKHiWkSkHdRznRkqrkUkV6i4FhFpB/VcZ4aKaxHJFSquRUTaoVBvIpNpKq5FJFeouBYRaYcCvolMRqm4FpFcoeJaRKQdNCwkM1Rci0iuUHEtItIOuqAxM1Rci0iuUHEtItIO6rnODBXXIpIrVFyLiLSDeq4zQ8W1iOQKFdciIu2gnuvMUHEtIrlCxbWISDuo5zozVFyLSK5QcS0i0g7quc4MFdcikitUXIuItINuIpMZKq5FJFeouBYRaQfdRCYzVFyLSK5QcS0i0g4aFpIZKq5FJFeouBYRaQdd0JgZKq5FJFeouBYRaQf1XGeGimsRyRUqrkVE2kE915mh4lpEcoWKaxGRdlDPdWaouBaRXKHiWkSkHdRznRkqrkUkV3TO9AHN7CTgeuA8d18aresG3AxsA/oCV7r71uZeExGJW773XGdLzlZxLSK5IqM912Y2GOgJHJny0p3Ai+5+FfAw8MsWviYiEqvknmvPs8ovm3K2imsRyRUZLa7dfS0h2dYys2HA+cCj0apngc+YWWlTr2UkYBGRZnTu3JmioiJqamqorq6OO5y0yqacreJaRHJFxsdcu3tNyqqpwHp33xG9Xg0sAaY085qISFbI56Eh2ZKzVVyLSK7IhgsahwMbUtZtA4Y185qISFYosIsaY8nZKq5FJFdk/ILGBjiwM2VdCbA7Wm7qtVqVlZWUlZXVPi8vL6e8vDyNYYqINKwtPdcVFRVUVFQkng5If1QdJpacreJaROLUmpydDcX1SsIV5cl6AasISbyx1+oZOHAgs2fP7pAARUSa0pae6+Ri0szWd0hgHSOWnK3iWkTi1JqcnQ3DQp4DhkXTN2FmXYBRwPPNvCYikhXyecx1A2LJ2SquRSRXZLy4NkukSAzA3dcAjwEnROtPBGa6+8qmXstgyCIiTcrnMdfZkrNVXItIrsjosBAz6wVcGD29yMzucPdK4BLgZjObCAwBkgfeNfWaiEjs8rXnOptytoprEckVGS2u3X0b8IuoJa/fBExv5D2NviYikg0SxXW+9VxnU85WcS0iuSIbxlyLiOS0xLCQfOu5ziZFReFxz5544xARaY6KaxGRdsrXYSHZpEuX8Lh7r0n9RESyi4prEZF2yucLGrNForjOszvMi0geUnEtItJO6rnueJ2jK4TUcy0i2U7FtYhIO6nnuuMlius9e3RRo4hkNxXXIiLtpJ7rjmem3msRyQ0qrkVE2kk915mhixpFJBeouBYRaSf1XGdGoudaFzWKSDZTcS0i0k75ehOZbKOeaxHJBSquRUTaSTeRyQwV1yKSC1Rci4i0k4aFZIaKaxHJBSquRUTaSRc0ZobGXItILlBxLSLSTuq5zgz1XItILlBxLSLSTuq5zozi4vC4a1e8cYiINEXFtYhIO6nnOjO6dw+PH30UbxwiIk1RcS0i0k7quc6Mnj3Do4prEclmaSmuzezT6diPiEguyrWe61zN2T16hEcV1yKSzTq3ZmMzOwW4BhhCXWHeCRgMdEtvaCIiuSFbbyKTbzk70XO9fXu8cYiINKVVxTVwL3AD8AZQE60z4Px0BiUikkuy+CYyeZWz1XMtIrmgtcX1B+5+S+pKM1uYpnhERHJOtvZck2c5Wz3XIpILWjvm+hIzO6eB9Z9JRzAiIrkoi3uu8ypnq+daRHJBa4vrY4Bfm9mepFYD/KIDYhMRyQlZfEFjXuXsXr3C49at8cYhItKU1g4L+RFQDswG9kTrDPhCewMxswOBy4DFwH5Ahbv/w8y6ATcD24C+wJXurtQqIlkji6fiy6uc3b9/eNy4MR17ExHpGK0trp8EHnb3en9BzOxXaYjlXuBMd19pZqOAp4EDgTuBp9z9gejK918C09JwPBGRtMjinuu8ytmJ4nr9+nTsTUSkY7S2uP498O9m9rekdQacDXyjnbEcCERf+rET6GNmwwhXtf9rtP5ZYKaZlbr70nYeT0QkLZJ7rt0dM4s5olp5lbMHDAiPKq5FJJu1trj+PuHrv4sAj9YZYQ7V9ibqPxLGBp4KXAhcDkwF1rv7DgB3rzazJcAUYGk7jycikhZFRUV07tyZ6upqdu/eTXFxcdwhJeRVzlZxLSK5oLUXNP4EGOnupe4+JmqlhF6Q9roU+Bh4Hdju7jOB4cCGlO22AcNS31xZWUlZWVltq6ioSENIIiIt09rp+CoqKmrzFTCgg8LKq5yt4lpE4tKanN3anusV7r6tgfVDWrmfhnQD/peQnG+Nejuc8HVjshJgd+qbBw4cyOzZs9MQhohI65WUlLB9+3Y+/vhjeiWmtWhCeXk55eXlAJhZR5WLeZWz+/ULj5s2wZ49UFTU+qBFRNqiNTm7tT3Xl6auMLNzgVtbuZ+G3AfMdPcfATcSkvZqwtXmyXoBq9JwPBGRtOncOfRVVFdXxxxJPXmVszt3hn32gZoa2Lw5HXsUEUm/1hbX/czs+wBm1t3Mfg38mvC1YJuZ2QBgorsn0uX1wFZgGTAsmtoJM+sCjAKeb8/xRETSrSjqRt2zZ08zW2ZU3uXsQYPC45o16dqjiEh6taq4dvfTgGfM7HbgH8D+wESgoTuAtcZG4GMzG560bgPwJvAYcEK07kRCT8nKdh5PRCStsrG4zsecPXp0eFy6NF17FBFJrybHXJvZPzXy0jLgPeCnwGjgLOCKtgbh7jVmdhZwrZnNAQYDM9x9q5ldAtxsZhMJ4wTL23ocEZGOkg3FdSHk7NLS8LhkSTr3KiKSPs1d0HgzYS7TzUBNA6/fRZjWaSjtSNQA7j4P+FoD6zcB09uzbxGRjpYNxTUFkLPHjAmP6rkWkWzVXHH9I2C2u69raiMza+9XjCIiOS1Liuu8z9mJ4lo91yKSrZosrt39iZbsxN3/mJ5wRERyU6K4rqlpqMM4MwohZ48bFx7ffTfeOEREGtPa2UJERKQBnTqFdJpNFzTmowMPDI+LFsHuvWbPFhGJn4prEZE0yJJhIXmvR48wNGT3bnjvvbijERHZm4prEZE0UHGdOQcdFB7feiveOEREGtLu4trMDjOz+8zs22ZW2v6QRERyT64U1/mQs1Vci0g2S0fP9deBE939Z8DRiTtziYgUklwprsmDnP2JT4RHFdciko3SUVxfBhweLf8e6JKGfYqI5BQzA8DdY46kWTmfsw85JDzOnRtvHCIiDWl1cW1mp5rZBdHykcB+iVvbunuNu29Nc4wiIlkv0WOd6MHOFvmYsydMCBc2fvABrF0bdzQiIvW1qrg2sx8DDwKfBXD3vwNfNLPzOyA2EZGckY3Fdb7m7M6d4cgjw/Irr8Qbi4hIqtb2XB8NjABeSFr3JHBj2iISEclB2Vhck8c5++ijw+PLL8cbh4hIqtYW18+4+2YgeVBhzt5GV0QkXbK0uM7bnJ0ortVzLSLZpsnbnzfgAzO7ChhuZp8HziV83Xh52iMTEckhWVpc523Onjw5PL7+OuzaBcXF8cYjIpLQqp5rd/898DegCPgisBn4Z3f/eQfEJiKSM6qrqwHo3Lm1fRYdJ59zdr9+cMABUFUFb7wRdzQiInVa/VfA3V8FXk08N7MiMzvB3f+a1shERHJIlvZc53XOPvZYWLgQnnkGPvnJuKMREQmaLK7N7BGgZzP7GACsB3I+UYuItFU2FNeFlrNPPx1+9St45BG45pq4oxERCZrruX4z2uZdoAaYDrwBzE7aZgzQr0OiExHJEYlhITH3XBdUzj7xROjaNYy7XrUKhg2LOyIRkeaL658AFl1tjpmNdfdrUzcys4c7IjgRkVyxc+dOALp37x5nGAWVs7t3DwX2Y4/B00/DxRfHHZGISDMXNLr7lkSSjgw2s3rdMmY2FTiyA2ITEckZO3bsAOItrgsxZ590Unj8a84PchGRfNHaCxr/CMwzsxeBbcCBwMnAXj0jIiKFoqamho8//hiAbt26xRxNPXmfsxPF9dNPw549kGXXk4pIAWrtVHzPAP9MuBhmPGFap3Pc/YYOiE1EJCck91qbWczR1CmEnH3AATB2LKxfr7s1ikh2aMtUfKuA2uuyzayTmR3u7nPSGpmISI7IhiEhjcn3nG0GZ54Jt9wCf/xjmJ5PRCROrSquzeyuBlYPBPYAZ6UjIDObDHwKWAy8CHwE3Ez4SrMvcKW7b03HsURE0iFRXPfo0SPmSOorlJz9+c+H4vr3v4ef/Qyy6D4+IlKAWpuC9iPMjepJ60YAL6QjGDO7GBjr7sm9LP8DPOXuD5jZKcAvgWnpOJ6ISDpkcc91QeTssjLYf39YtChc2HjyyR15NBGRprW2uP6qu7+bvMLMDgAuam8gZjYFuAA4KWndMOB84F+jVc8CM82s1N2XtveYIiLpsH37diAri+uCyNlmcOGFcO21cN99Kq5FJF6tvaDx3QZWVwFfTUMstwDvALeb2VNmdgwwFVjv7jui41cDS4ApaTieiEhabNy4EYB+/bLr3iyFlLMvuCA8/vGP4YYyIiJxaVVxbWZLzOyDpLaMkFz/1J4gzGw8cBjwK3e/jPA15lPASGBDyubbgL3uw1VZWUlZWVltq6ioaE9IIiIttmFDSFP9+/dv8XsqKipq8xXhluRpV0g5e9994eyz4eOPYcaMdu1KRGQvrcnZrR0W8gxwP3Xj92qAte6+qNVR1ncQsNHd34ye3wH8IDrOzpRtS4DdqTsYOHAgs2fPTl0tItLh2lJcl5eXU15eDoCZre+QwAosZ990Ezz1FNx/P5x1Fpx7blp3LyIFrDU5u1U918B97j7L3Z+P2gvhGDaujbEmdCap0Hf3ncB7QBHhavNkvQB96SciWaMtxXWGFFTOHjcOfvrTsPz1r8OHH2biqCIi9bW2uD6lgXVrgQfbGcebQB8zG5y0rpqQkIeZWTcAM+sCjAKeb+fxRETSJouL64LL2V//Opx6KmzcWDdMREQkk1pUXJvZZWa2GbjSzPYkN8L4urXtCcLdFwJPAudEx9sHKAZ+BzwGnBBteiIw091Xtud4IiLptHZtSIEDBw6MOZKgkHO2GdxzD4weDa+/DtOnQ01Npo4uItLCMdfufoeZ/R/wJeC/Ul7e5e6r0xDLF4HbzKyY0NMxzd2rzewS4GYzmwgMAcrTcCwRkbRZtmwZAKNGjYo5kqDQc/aAAfDnP8Mxx4Sp+QYNCsNFsujO9CKSx1p8QaO7v2lm33P3vS5MMbNR7r68PYG4+3rCnKmp6zcB09uzbxGRjpQorktLS+MNJEmh5+xDD4WHHoIzzgh3b+zXD66+Ou6oRKQQNFlcm9l3gCfc/e1o1XfNrCh1M+A4wvymIiIFZceOHaxbt44uXbowdOjQWGNRzq7v5JPh3nth2jS45powDvvmm6Eo9RMREUmj5nquxwD7JD2fBPQA1lA3tZMRvvoTESk4yUNCOnVq7TXiaaecneK886CqCr7yldCDPWcO/Pa3YUy2iEhHaLK4dvd/TVn1Y2CBu+9KXmlmn0h3YCIiueDtt0Mn8X777RdzJMrZjfniF2HkyNCD/fzzcPDBcOutcPHFGoctIunX2m6WgxtI0uOBsekLSUQkd8ybNw+AiRMnxhxJg5SzI8cdB/Pnh+n5tm0LPdlTpoR1IiLp1KILGs3sWMJXiSeZ2ZKUl/sRrkZ/NM2xiYhkvWwsrpWzGzZwIPzxj+EOjt/6FrzwAhx2GFx6aRiTPWhQ3BGKSD5o6Wwhm4F7COP5jkl5bRfw63QGJSKSC9y99hbehx56aMzR1KOc3QgzuOACOO20UFDfeSfcfjvcdRdccQV85zvQp0/cUYpILmvRsBB3n0+4svw77j4mpY139x91aJQiIlnovffeY9WqVQwcOJADDjgg7nBqKWc3b5994I474B//gNNPh48+guuvDxc6fu97sCojN2wXkXzU4jHX7r7F3Rvs7TCzT6UvJBGR3PDcc88BMHXqVCzLroxTzm6ZQw6BRx+Fl16CqVNhyxa46SYoLQ0XPP7jH3FHKCK5prl5rh8BejazDwMOAjRaTUQKylNPPQXA8ccfH3MkgXJ22x19NDz3HLz2GvzsZzBzJtx9d2hHHgnl5fD5z0OPHnFHKiLZrrkx1+8Sxu6tpG6O1FQGnJ7OoEREst3WrVt58sknMTM+85nPxB1OgnJ2O33yk/CHP8D778N//VeYE/vvfw/tW9+CCy8MhXYWXb8qIlmmueL6JmCbu1c1tZGZFdxV5yJS2B599FGqqqo49thjGTZsWNzhJChnp8nYsXDbbXDDDfDgg/DLX8Irr4QLIO+8M8yVPW1a6M0eMybuaEUkmzQ55trd1zeXpCOfTVM8IiI54Ve/+hUA5513XsyR1FHOTr/u3eFLX4KXX4Y334TLLoO+fcP82FddBfvuC0cdFQrx6GadIlLg2nWvXjMbZWY3AbelKR4Rkaw3b948nn/+eXr27MkXvvCFuMNpMeXs9jn44DBt35o18Mgjoee6e3d49VX45jfDRZCHHgrXXguzZ4M3NjBHRPJaq4trMysys8+a2RPA+8AFhDF+IiIF4ZZbbgHg4osvpnfv3jFH0zTl7PQrLoYzzgg3o1m3Ljyecw707Anz5sF118ERR4Rbrn/1q2EM94YNcUctIpnS0pvIYGZjgK8CFwP9gWXAGcDTwFEdEp2ISJZ56623uO++++jcuTPf+MY34g6nUcrZmdGjR+jBnjYNqqrCjCOPPBLaypXwm9+EZgaHHw4nnQT//M8weTJ07Rp39CLSEZrtuTazc83sL8B7wL8AtwAjgMfd/SkPXu7gOEVEssJVV11FTU0N5eXljB07Nu5w9qKcHZ+SEjjllHDB44oVMHcu3HgjnHACdOkShorccAMcd1y4ic2xx4Zx208+GebXFpH80JKe652E2+W+DHzV3RcBmJlGk4lIQXnooYd45JFH6NGjB9///vfjDqcxytlZwAwOOyy0K6+EHTvghRfgL3+BZ54JF0S++GJoN9wAnTqFG9oce2xdGzIk7rMQkbZotrh298eAx8xsFDDdzEqBh4EuiW3MzNx16YaI5K+NGzdy6aWXAnDTTTcxJEsrH+Xs7NS9O5x8cmgAmzaFu0K+8EJos2eHu0H+4x/hokkI0wEec0yYjWTyZJgwATq3eDCniMTFWptfzawzcBbhopgFwG+Af3b3ivSH13JlZWU+e/bsOEMQkTxVU1PDGWecwRNPPMGnPvUpnn/+eTp1atdkS/WY2Rx3L0vbDuvvWzk7B+zYEW5Ukyi2X3kFtm+vv02PHuFukUcdVVdwDxgQT7wihay5nN3q4jpl5+OAy4Hp7t69zTtKAyVqEeko1113Hddeey39+vVj7ty5jB49Oq3778jiOuU4ytk5oro6zDzyyit1bcmSvbcbMwYmTQrt8MPDMJRBurG9SIdqLme36wsmd18MfMPM5rRnPyIi2eq+++7j2muvxcy4//77015YZ5Jydu7o3DkUy4cfHm5cA7B2bZhTO1Fsv/56KLiXLIGZM+veO2JEXcGdaMOGhXHgItLx0jJ6y93vScd+zKwYeB34hrvPMrNuwM3ANqAvcKW7b03HsUREmvP4449z8cUXA/CTn/yEkxMDZnNcunI2KG9n0uDBcOaZoUHo3V64MMxKMncuzJkDb7wBH34Y2iOP1L23f/9wE5wJE+oeJ0yALJ+mXSQnZdulETOA0qTndwJPufsDZnYK8EtgWhyBiUhheeihh5g2bRrV1dVceeWVfPvb3447pGylvB2Tzp3riuQvfjGsq6mB996rK7gTbcMGmDUrtGSjRoViO7nwHj8+TCsoIm2TNcW1mR0NrAY2Rc+HAecD/xpt8iww08xK3X1pLEGKSEG46667mD59OjU1NXzjG9/ghhtuiDukrKS8nX06dQrF8fjx4cY2EG7D/uGHsGBBmAIw8fjOO7B8eWiPP15/H2PGwAEHhP0ccEBdGzBAw0tEmpMVxbWZ9QDOcfdvmdm10eqpwHp33wHg7tVmtgSYAiyNJVARyWvV1dXMmDGDW2+9FYAf/vCHteOtpT7l7dxhFm7FPnIknHpq3frqali8eO+i+/3361py0Q3Qt2/9YjtRfO+7b7hRjohkSXENXAncmLJuOLAhZd02YFhGIhKRgrJmzRouvPBC/vrXv9KlSxfuuOMOysvL4w4rmylv57jOneuK5HPOqVtfVRUK64UL69q774bHTZvqLqhM3de4cXv3dI8fHwpykUISe3FtZicDs919XcpLTrjTWLISYHdD+6msrKSsrG5WlPLycv1hFJEW+fOf/8z06dNZv349gwYNYubMmXzqU5/q0GNWVFRQUVE71XROzVacjrytnJ29SkrgE58ILZk7rFmzd8G9cCEsW1a3/PDD9d83cCDst19d23//8DhuHPTsmbnzEmmP1uTsds1znQ5m9gxwWNKqvsB24BbgAnffP2nb94AfuPv9qfvRnKki0lrr169nxowZ3H333QCceOKJ3H333QwfPjyjcWRqnut0SUfeVs7OLzt2hAspk3u7EwX4ztT/biUZOrSu2E4uvseMCXe1FMlGHTrPdZqcT+jZSHgF+BbwGjDDzLq5+04z6wKMAp6PIUYRySM1NTX89re/ZcaMGWzYsIGSkhJuuukmLr/88rTeeTGPKW9LPd27w8SJoSWrqYFVq2DRolB8v/de3fL778Pq1aE938C/kCFDwljuffcNxXZied99w7zd+lWVbBV7ce3ulcnPzWwPUOnuy83sMeAE4DHgRGCmu6+MIUwRyRMvvvgiM2bM4NVXXwXg+OOP584772T8+PExR5Y7lLelpTp1Cje1GTECjj++/mt79oSZShoqvJctC0NQ1qyBl1/ee7/FxVBaWr/gTi7CNX+3xCn24roZlwA3m9lEYAigAXki0iZvvfUW//7v/86jjz4KwODBg7nllluYNm2aZgNJL+VtaZGiolAMjxkDqfdn2rMnTB+4ZAl88EFdSzxfty4U4osWNbzv/v1D8T16dJjLe/ToujZqVHhdv/bSUbKuuHb30qTlTcD0+KIRkVw3Z84cbrjhBh566CHcnR49ejBjxgy+/e1v01NXU6WF8rakW1FRXTE8derer2/fDkuXNlx4f/BBuGnOhg3hrpUN6dGjruhOLr4Ty8OGhRlQRNpC/3REJO+4O7NmzeLGG2/kL3/5CwDFxcVMnz6d73//+wwePDjmCEWkPXr2rLs7ZarErCbLltW15cvrP9+6NdxE5513Gt5/UVEYytJQ4T18eGj77KPeb2mYimsRyRtbt27lvvvu48477+Stt94CoGfPnlxyySVcccUVDB06NOYIRaSjmYVZSIYOhcmTG95my5aGC+/E8urVdeteeKHhfXTvHorsESPqCu7EcuJx8OBQqEthUXEtIjnN3Zk7dy6/+c1vuPfee9m+fTsAQ4YM4ZJLLuGyyy6jX79+MUcpItmkTx845JDQGlJVFcZ8pxbdy5fDypXhte3b6y7EbExRUSjyGyu+hw8PQ1A07WB+UXEtIjlp+fLl/O53v+Pee+/lnaTvdqdMmcKll17KWWedRXFxcYwRikiuKimBsWNDa8zWraHQThTbDT2uWxeWP/yw6eP16ROK7NQ2dGj95a5d03ue0jFUXItIzli9ejUPP/wwf/jDH5g1aymse1YAACAASURBVBaJm2ANGDCAadOmUV5ezoSGBmGKiKRZ796hHXhg49tUVYUhJskFd2oRvnp1GKayZUvjY8AT+vVruPBObkOGhKkKJT4qrkUkqy1ZsoQ//elPPPTQQ7z88su1BXVJSQlnnnkmX/jCFzj55JPp0qVLzJGKiNRXUhKmBCwtbXwbd9i4MdxsZ/Xq8NhQW706bLdxIyxY0PRxBwxouPBOLsoHDwalzY6h4lpEssquXbt46aWXeOqpp3jyySeZP39+7WslJSWcfPLJfPazn+Wzn/0sffr0iTFSEZH2MwvzbvfvDwcf3Ph2NTVhesHGCu/E8po1sH59aG++2fRxBw1qfBhKog0apIsyW0vFtYjEbsmSJTz99NM8+eST/O1vf6u9KBGgV69enHbaaZx99tmceuqpmptaRApSp04wcGBoqbeZT7ZnD1RWNlx4J7e1a+vaG280vr+iojDUZNiwugswk5cTj5qasI6KaxHJKHdnyZIlzJo1i+eff55Zs2axfPnyettMmDCBU089lVNOOYVjjjmGkpKSmKIVEcktiWJ4yBCYNKnx7aqrwwWXTQ1DWbkyFOqJseKvv974/rp1a74AHzYsbJfvVFyLSIeqqalh0aJFvPTSS7UF9YoVK+pt07dvX44//nhOPfVUTj75ZEaMGBFTtCIihaFz57qCtym7dtXv/V65suHHbdvg/fdDa0rfvnUF9777wn771bUxY8I49Vyn4lpE0mrz5s38/e9/55VXXuHVV1/l1VdfZfPmzfW26devH1OmTGHKlClMnTqVgw8+mE6dOsUUsYiINKa4uO4ulU3Ztq35AnzVKti0KbToPl/1dOoU7oSZKLYPPBAOOgg+8Ykw9jtXhp2ouBaRNquqqmL+/PnMnTuX1157jVdffZW33357r+2GDh3KUUcdxdSpU5kyZQoTJkxQMS0ikkd69YLx40NrTPJFmStWhF7uxI143nsv3Khn6dLQnnmm/nv79QtFdmobNiz7im4V1yLSIjt37mTevHnMnTuXuXPnMmfOHBYsWEB1dXW97YqLi5k0aRKTJ0/mqKOOYvLkyYwcORLLtuwnIiIZ1dxFmVVVsGRJKLQXLQrzfr/9dmgbN8KLL4aWrHfvUGRPmABlZXDEEWHWlTinGVRxLSL1uDtr1qxh/vz5zJ8/nzfffJO5c+fyzjvvsGfPnnrbmhkHHHAAkyZNoqysjKOOOorDDjtMFyCKiEirlZTAAQeElsw9jPtOFNqJ9tZboeh+9dXQfv3ruv0cemgotI89Fo47LhT0maLiWqSAbd++nQULFtQW0om2YcOGvbbt1KkTEyZMYNKkSUyaNInDDz+ciRMn0qtXrxgiFxGRQmFWd/HliSfWrXcPs5m8/TbMmxdmM3n99dDr/dprod1xR9h24kQ44QQ44wz4p38KvegdRcW1SAHYuHEj7777LgsXLmThwoW88847LFiwgCVLljS4fZ8+fTj44INr22GHHcYhhxxC9+7dMxy5iIhIwxI3whk0CKZOrVu/eTPMnh2K6+eeC0NJ5s0L7ZZbwkwl5eVwxRVhrHja40rcSjjXlZWV+ezZs+MOQyQ2e/bsYdmyZbUFdHIxvW7dugbf06VLFw488MB6hfTBBx/MiBEjNEY6g8xsjruXxR1HJilni0imfPwxvPwy/OUv8MAD4YJJgHHjwnCS/v1bt7/mcrZ6rkVyyJ49e/jwww9ZvHgx77//PosXL65tixYtoqqqqsH3de/enQMOOKC2jR8/noMOOoj999+fLnFe9SEiItLBunaF448P7YYbYOZMOPdcWLwYHnoIpk9P7/FUXItkmV27drFs2bLaojm5iF6yZAm7du1q9L3Dhg2rV0Qn2vDhwzX1nYiIFKTqapg/H156CZ54Ap59NqwvKYGjjkr/8VRci2TY7t27WbFiBcuWLWPp0qW1LfF8xYoV1NTUNPr+oUOHMnbsWMaNG8e4ceMYO3YsY8eOZfz48fTu3TuDZyIiIpI9amrqZhV56626x3nz4KOP6rbr1Ak+/Wm4/vowhV+6qbgWSbNdu3axYsWKBgvnpUuXsnLlyiaL506dOlFaWlpbQCcX0vvuuy89evTI4NmIiIhkh6oqWLOm7q6PS5bUb0uXhm0aMnYsHH00TJkCn/lMx07Np+JapBWqqqpYuXIlK1as4MMPP+TDDz+st/zhhx+ydu3aJvdhZowcOZLRo0dTWlq6Vxs5ciTFxcUZOiMREZH47N4N69eHKfWS27p1dbdMT7T165vf34AB4S6REyaEW6cfdFC4qUxBznNtZmcDPwN6A/cB33b3ajMbAPwY2Ax0Aa50993xRSr5avv27axevZqVK1c2WDSvWLGCysrKZvfTqVMnhg8f3mDhXFpayogRI1Q8S85TzhaRZO6wbRts2tRw27hx7wK6sjJMm9dSRUUwdGiY73roUBg9GvbdF8aMCa20tGOm1mutrCiuzWwUcBZwDnAQ8AtgBfBT4EHge+7+mpl9HfgR8O9xxSq5xd3ZvHkzq1evrm2rVq2q9zzRtm/f3uz+ioqKGDZsGCNHjmTEiBGMGDGi3vKIESMYMmQInTtnxa+WSIdQzhbJP9XVsHVr423LlsYL502bQpGcchPfFunUKUyFl7gteqINGlR345hEGzgwFNjZLlsqgNHAl929GphjZhOA48zsRWC8u78Wbfco8K6ZXe/uHzW2M8l/O3bsYN26dfXamjVr9iqe16xZw8cff9yifZaUlDB06NB6xXNqET148GCKcuE3W6RjKWeLxMw9jC/evj1crLd9+97LqcVxU8Xzjh3tj6lHD+jXD/r2bbglCucBA+qW+/Xr2LslxiErimt3fyFl1UqgF3A8sCRpu5XRjS3KgOczFqB0uD179rBhw4baQnnt2rV7Fc/JrSW9zAk9e/asLZqHDh3aaNtnn3104xSRFlDOFmkZ9zCmeOfOULwmt4YK4saK5Ma2bUtPcWM6dYLevRtuffqE4RaNFc19+8I++4BGPAZZUVw34AjgJ0A5sCHltW3AsNQ3VFZWUlZWd7Oc8vJyysvLOzJGaYS7s2XLFjZs2FDb1q9fv9fy+vXra4vl9evX05q7hRYXFzN48GAGDRpU2wYPHtxg0dyzZ88OPFuRtqmoqKCioiLxdECcsaSBcrbknN276wrdhorf5Nae19NZAKcqLoaePUPr0WPv5T59Gi+YE0VzYrl793A7cWlYa3J21hXXZjYOqHT3eWbmwM6UTUqAvS6OGThwILqVbvpVVVWxefNmNm3a1GiRnPp848aNVFdXt/pY/fv336tgbqz17t1bvcyS05KLSTNrwTXw2Uk5W9or0bv78cehSN25s2OXE0VvG/5MtUlRUShck1u3bk0XxS1Z7tEDdIPdzGlNzs6q4trMioCvAd+LVq0E9k963YAewKrMR5eb3J0dO3awadMmNm3aVFsoN/Y8dd3Onal/J1umV69e9O/fv7YNGDCgweVEMT1gwABdBCiSY5Sz8497KEITvbEdXegmHpuY+r/DdOoUCtTkgje1CE4tiNvymgrgwpNt1cx3gJ+6e2IK8GeArye9XgrsAgqquyO1QN64cWODyw21zZs3s3t322fBKioqom/fvvTt27fZQjmx3K9fP0pKStL4CYhIllLO7mB79tQVpK1tyQVyS1sLr/9Ou86doWvXUKR261a33NC69i4nF736AlQ6QtYU12Z2FSEBdzOzfYGpwP8Bq81sgrsvAE4Bbnf3XfFF2j47duxg9erVrF27trYwbqxYTl7etavtp9y1a9faAnmfffapXW7J8x49emj4hYjspVBydkvV1IQ5frdtCzMvJJabe759e9PFcTv6RtqspKSuEE1uHVHoJpb1xaXkk6z452xm1wDXpaxe6O53mdl5wLVmthToCVyd6fhaavfu3SxbtozFixezePFili5dutdcylu2bGnTvhMFcr9+/eoVv6nPGyqUu3btmuYzFZFCli85uyW2bIFly+raunXhxhfr19fdVS6x3BEXrpk1XOi2pSWGLzTVunbNv2nRRDItK4prd78euL6R15YBX8lsRM1bu3Ytc+bMqW0LFixg6dKl7GkmuxYXFzN06FCGDBlC//79Gy2UU4vmbt26ZejMRESalos5uzlbtsDrr8Obb8K8eTB/PnzwQVjfUj17hlkXevWqe0y05OfJyz17Nl0AFxdr6IJIrsmK4joXfPTRRzz77LM888wzPPPMMyxatKjB7UaOHMm4cePYb7/9KC0tZdiwYfXmV+7Xr5+GWYiIZIF58+APf4Bnn4XZsxu+qK5bt3CL5UQbOrTuBhjJj/37h+EUIiIqrptQU1PDs88+y7333suf/vQnPvqo7gZjPXv25PDDD2fSpEkcfvjhTJw4kbFjx6qHWUQki1VXw733wm23heI6oXNnOPJIOPRQmDgRDjkE9t8/FM3qDxGR1lBx3YDq6mp+97vfcfPNN/P222/Xrv/kJz/JqaeeykknncSRRx6pqeNERHLIK6/ARRdB4ovHfv3g85+H00+HY48NQzRERNpL1WGKWbNmcfnll7NgwQIARowYwfTp07ngggsYO3ZszNGJiEhb/PKXcOmlYejHuHFw7bXwuc9pKIeIpJ+K60hVVRUzZszg9ttvB2DMmDH84Ac/YNq0aRQXF8ccnYiItNU998DXo9m3Z8yA665TUS0iHUfFNbBmzRpOP/105syZQ5cuXbjmmmv47ne/qynsRERy3OLFdYX1f/4nfPOb8cYjIvmv4IvrpUuXcuKJJ/L+++8zZswYHnjgAY444oi4wxIRkTS4+upwQ5bzz1dhLSKZUdDFdWVlJSeddBLvv/8+kyZN4qmnnmLgwIFxhyUiImmwahU8+GCYK/qmm+KORkQKRcEW11VVVZx55pksXryYww47jOeee47evXvHHZaIiKTJQw+BO3z60zBiRNzRiEihKNibnH73u9/llVdeYdSoUTz++OMqrEVE8szf/hYeTz893jhEpLAUbHH9pS99iQMPPJCZM2cydOjQuMMREZE0mzwZjj4apk6NOxIRKSQFOyxk0qRJzJ8/n6KiorhDERGRDvDd74YmIpJJBdtzDaiwFhEREZG0KujiWkREREQknVRci4iIiIikiYprEREREZE0UXGdgyoqKuIOISvocwj0OQT6HApDIfycdY75QeeYH9pyjiquc1Ah/GNuCX0OgT6HQJ9DYSiEn7POMT/oHPODimsRERERkRiZu8cdQ1qYWSWwLO44MmQAsD7uILKAPodAn0OQy5/DaHcfGHcQmdSOnJ3LP+eW0jnmB51jfmjoHJvM2XlTXIuIiIiIxE3DQkRERERE0kTFtYiIiIhImqi4zjFmdpKZvWZmpXHHEhczO9vMlpjZBjO7zcw6xx1THMxsipm9a2abzOzWuOOJk5kVm9k8M5sadyzSMcysm5ndbmb/z8x+YWa9444p3QoltxXK76uZTTaz75jZWWY2IO540snMDjSzn5vZFWZ2p5kdGndM6dBQjdWW3JOXv7j5yswGAz2BI+OOJS5mNgo4CzgHOAj4BbAC+GmccWWamfUBjgEmA0cDfzazx9z92Xgji80MoDTuIKRD3Qk85e4PmNkpwC+BaTHHlDYFltvy/vfVzC4Gxrr7NXHH0kHuBc5095XRv92ngQNjjqldmqixWp171HOdQ9x9LfBw3HHEbDTwZXef4+73AD8Hjos5pjjsAm5w903u/jjwJrAn5phiYWZHA6uBTXHHIh3DzIYB5wOPRqueBT6TZ9/gFURuK4TfVzObAlwAfD/uWDrQgUCvaHkn0CfGWNKioRqrrblHxXWOcfeauGOIk7u/4O7VSatWAsvjiicu7r7To6l+zKwH8A4wK9agYhCd+znuflfcsUiHmgqsd/cdAFEOWAJMiTOodCqE3FZAv6+3EHLy7Wb2lJkdE3dAHeCPwK/NrBdwIXB5zPGkRQM11lTakHtUXEuuOwL477iDiIuZHQs8CXQBusccThyuBG6MOwjpcMOBDSnrtgHDYoglU/Ixt+X976uZjQcOA37l7pcBfwWeioYc5JNLgY+B14Ht7j4z5ng6Sptyj4pryVlmNg6odPd5cccSo8XA3cAJwE/iDSWzzOxkYLa7r4s7FulwTvjqOVkJsDuGWDpcPua2Avp9PQjY6O5vRs/vAAw4M76QOkQ34H+B+4FbzezEmOPpKG3KPbqgUXKSmRUBXwO+F3cscXL31cBdZlZDuEiokHwHOMzMEs/7Ag+b2Y/d/eb4wpIOsJLw803WC1gVQywdKo9zW6H8vnYmqbZy951m9h7QP76QOsR9wOfdfbOZOfC/Zlbq7h/FHViatSn3qOdactV3gJ+6e1XcgWSJ2YQkUEjOBw5NaquAr5J/X6ULPAcMM7NuAGbWBRgFPB9rVB0jX3Nbofy+vgn0SRkGUg28HVM8aRdNKzjR3TdHq64HtpLjs4U0ok25R8V1jrG6//ZbkxvmMTO7CngD6GZm+5rZl6OvUQuGmXU3s8OSVn0auC2ueOLg7pXu/mGiEWZLqXT3rXHHJunl7muAxwjDnwBOBGa6e179hzKfc1uh/L66+0LCdTDnAJjZPkAx8HiccaXZRuBjMxuetG4DsCimeNImtcZqa+7RsJAcknRVLsBFZnaHu1fGGVOmmdk1wHUpqxcWwNXnqSYAj5vZQuBV4O/RlHwi+eoS4GYzmwgMAcpjjietlNvyyheB28ysmNDLOS1lJpic5u41ZnYWcK2ZzQEGAzNy/T9KTdRYrc49Fs3mJSIiIiIi7aRhISIiIiIiaaLiWkREREQkTVRci4iIiIikiYprEREREZE0UXEtIiIiIpImKq5FRERERNJExbVII8zsBDN72cymNrPdBDP7g5l9P43H/km69iUiEgczO8TMHkxnbmzkOI+Z2ec78hgdxcxONLNXmvs7I7lFN5GRrGdmBwI/JtwV6gjgEOAL7n5fBx96BTCxBdstB/oCRek4aHR71d3p2JeICICZHQHcAPwT8BdgCeGGGF2AX7j70x1w2MVAH9KUG5twJ7CgI3Yc5eOzCbdpXwS8Er1UAvwLcKS7L23HIVr6d0ZyiIpryQV/Bj7n7vPMrBPwP5k4qLsvMrMNLdhuq5mtTuOhpwD/l8b9iUiBc/fXzez3wAHufnpivZmdBvzJzL7j7nem+Zg70pwbATCzM4D5iaLW3Z9I9zES3H038ICZXQ7McvdrkuJ4uC37NLNBwHHu/oC7v2tm69MUrmQJDQuRrGZm/YD9gY8g3HYVuBbI1K1Fa9K8XUscBzyXxv2JiADsdQtud38c+CHwn2Y2pAOOmdZcbWYjCD3VmdbQ7cufBVa2ZidmVgLcA3RLWp3Ovx+SBVRcS1Zz943Au8BTZjY5WrcM+LuZ7WNm15nZ383sM2a2xszeMbNJifeb2TQzu9rMZpnZ/0taX2ZmPzSzB8zsITPrE63vamY/N7Nrzew3hK80W8XMvhvF9DkzO8LM/hzt77ooxgfMbGC0fouZ/SBlF8XuXhXt69tm9i0ze8nMbmz1Bygi0rzfEIaHnAtgZmPN7Edm9msz+2tU0GJmJVEuu9rMXjCz8xI7MLNLzez7ZvaqmV3W2IEayslmdk6UMy8ys/fM7EYzG2Zm/xPl0xfM7PhoF+cDI4Arzeys6NqYZ83sS0nHmG5mPzaze8zst2bWw8z6mNkPzGy2mR1vZu+a2QdmVtqWDyz6TKZGPduY2dlmdpuZ/czMnjKzntE53GZmj0Z/b1ZHn/EngHPN7FtJu+wb/S3aZmbTGzlmJzO7yczKo2O4mf0p8fORLOLuampZ3YCDgGWE/93fBQyK1hcREtVGwti3fYBngLcI/3E8Crgm2rYvoff7NKA38Puk/b8E/CRa/k+gPFoeDOwhJNDmYryb0PsziDA2r3u0vgR4GngYGEb4o7Ab+AHQHTg8ep7YPjGePHHed0bLvYDvxf2zUFNTy90GXAR82Mhra4H/ivLqQ0BRtP5+4MFo+Rbg7Gj5y0BltDwNuDBa3j/K1QdHz+8GfhgtN5aT+wGbgVuBI4ETgZ8C10bbzgAeTYrVgdJouRvwDnBR9Pxk4P6kbR8DfkkYBvvZ6O/F2dF5PgX8uIWf3azob8WNwE8I469PTHp9PTAqWn6L8DepBLga+BD4JPCVKO/PSsQbbb8U+Fn02jRgeSMxfBV4LFoeHH3O/xz3vyu1vZvGXEvWc/e3zOxg4HrgUuAMMzvRwxjsSmCru88EMLMfAi8C+xL+kNSY2TejXT1NSOKnAf2T1r8TvbcncEn0Xtx9rZmtaEWoA4Dbga+4+45oH1VmtgZY4u6rouOsAZ73MB7xDULSH0C4MPI0wh8CgB3AF81sEfALwn8sREQ6ghMKziOB0cDlZgawAehj4cnXCB0DAL8lDIuAkGtnp+TaocD8lGNcRAM52d03mtkW4GF3/zuAmS0HPjKzoYSe3p4NBu2+08zWJa26LNpvwq+BP0TrNxP+XjwUHWN2FGdLPefRmGsze4DQoZNwrrsvN7NjCEV1zyj/rwLed/fXgNei9za070ejvwmvNxHTJGBndN5rzWwj9YeXSJZQcS1ZzcyKgJEeLlz5t2ioxhNABaEnINVb0eM+wCjgXne/P1p3a7TPKwnF7q0pxzqckBS3tTHc/YFPAWOBeUnrU8cc7ql9wb0mSrSJIVr9PQyFwd2XmNlFwB3AFcAFwDpERNLIzPYhfOv2FiFvbmggPw4m9KwWA7j7HkKHANF7/sPdX46e13tvkgZzcsSpnys/BK4E3gNeBcY0cQrJ79sP+FvS8/cJQ14GsXcurqbtM5nMJfwnJGGzmd0EPEiIPVFBp55Xc6ppvDZ7EviFmXUl/Cx2EnrBJctozLVku9pxgADuPg+4Cjigke2LCYlsMbAaOCf5RTP7ZLT+NDPrlrJ+a/T0wDbG+jJwG/C7KPm1ipkNI+niGDMbTpgpZX/CBY73N/JWEZH2+BLhm7IHCfnxmKjHGKjNj+uBXcApSesPNrNepORaM+scdVakaiwnN+S/gaUeplytasW5LCf0dNceAthCGPaSNu5eE3WAHBT95+QJ4EfuPjudx0k55qPAo4RvWKcBU9x9S0cdT9pOxbXkgsujwjNhOGGsXEI/MyuOlqcSxtttBv4X+KyZ3WFmU6NehY8ISbAH8JiZnRxdPHIAoSB/A7jOzIrNbABhrPMgMysys++Z2YWNxGhRu5bwB+iWBl5L3T71+emE8YEJ+wFnuftW4PIG3iMi0hp79Yia2T8B/wF8zd0rCb3Ea4EnzOyM6KLFY6Oe6t8TZhU5y8xOAL7k7tsIufbfzOw/zGwKoZMhMaQuOf81lpMTknuRDwMGmll34Bigm5ntG722m3ABYKKTJfkY/w38i5n1jZ5/EviVu1cTap4Gc7GZTY7iKmnks+tESs0U/QfkNMK3lX0Jww3HAyOB7kkXS6b2ju9qIP76QTUwdiT6WfQG/kronW/NfzokgzQsRHLBSOAtM3uc8JWZAV9Per0auNrMNhN6eS8HcPdnzOzfCF8tnkW4IHABgJmdCfyc8MfiLmCGu7uZfY4wTdJbwAOEr/cOJdx04XOE4Sb1bl5jZkcSLtQZRfia8M/Af0TxPA1MBkrNbCyhR2UI8DkzW0hdL85FQA93r0g59/+xcPOHasIFRCIirRblkfMJBeCvgErCRXF9gFPd/RUAd98V5cdfEXLdw4Sx1gDfIHRM3EP4pi4xQ8dvCLnvEsLwtUvcfV1UPE4GxprZfY3lZDM7hzDO+KtmttDdVxIurvwpobC+BzgTOBj4ALiX0Mv+RQuzSB0EfNrMnnD3h8xsP2Cmmb1E+Pbz6qiH/VxgSHR+bxFuqDPUzA6K4vw8YRz5n5M+ty6E3D8RGB4V304Y63x2tM9/AM8DrxOGujxJyOmPES6inGhmZ7l7Yr8PEG7osz0qxocA55jZEkKPNNH7U+/psJxwweZZQFegk5k9S7ioMVPT00oLmH4eksss3DL2bncvzdDxZri7bk0uIpJnoiEqJe6elTfxMrMZhJlblkbPuxK+dbgpca2OZAcNC5Fc19CQi445UPj6dG4mjiUiIpkTXbA5KYsL667ANUBp0uoBwEcqrLOPhoVIzorGu11I+JrvXHd/sIMP+X70daWIiOSXzu7+i7iDaIy7f2zhFuz/FV3svoowhOWapt8pcdCwEBERERGRNNGwEBERERGRNFFxLSIiIiKSJiquRURERETSJG8uaBwwYICXlpa26j2VlZUMHDiwYwLKAMUfL8Ufv1w/h0T8c+bMWe/uuXsibdCWnA25/zOH/DgH0Hlkk3w4B8id82guZ+dNcV1aWsrs2a2762hZWVmr35NNFH+8FH/8cv0cEvGb2bK4Y8m0tuRsyP2fOeTHOYDOI5vkwzlA7pxHczlbw0JERERERNJExbWIiIiISJoUdHFdXl4edwjtovjjpfjjl+vnkOvxxyEfPrN8OAfQeWSTfDgHyJ/zyPhNZMzsJOB64Dx3XxqtOxv4GdAbuA/4trtXR68NAH4MbAa6AFe6++7U/ZaVlXkujNMREUllZnPcvSzuOBqinC0iUl9zOTujPddmNhjoCRyZtG4UcBZwDnAF8FXgm0lvexC4y92vBBYBP8pYwCIiBUw5W0Sk9TJaXLv7WuDhlNWjgS+7+xx3vwf4OXAcgJlNBsa7+2vRto8Cl5tZj0zFLCJSqJSzRURaL+Njrt29JuX5C4mvEyMrgeXR8vHAkqRtV0aL6fn69JRT4KCDYNOmtOxORCTfZEvO3vbGNmYfPpuFFy9s765ERDpUNs5zfQTwk2h5OLAh5fVtwLDUN1VWVlJWVpe/y8vLmx8Y/847sHw5bNkCffu2J2YRkVapqKigoqIi8XRAnLG0U0Zyds3OGrbP3U6n4oK+Dl9EYtKanJ1VxbWZjQMq3X1etMqBnSmblQB7XRwzcODA1k883rt3eNy6tZWRioi0T3IxSUOZGAAAIABJREFUaWbrYw6nTTKZs62LAVCzu6aZLUVE0q81OTtrugDMrAj4GvC9pNUrgb5J2xjQA1iVloP26RMet2xJy+5ERApFpnO2dQ7FtVdndoYrEZHWypriGvgO8FN3r0pa9wwwPul5KbALSM/8Teq5FhFpq4zm7ETPte9WcS0i2S3jxXXUkwFgSeuuAt4AupnZvmb2ZTMb5+6zgdVmNiHa9BTgdnfflZZgEsW1eq5FRBqULTm7U5fw50rFtYhku4yOuTazXsCF0dOLzOwOwteK16VsutDd74qWzwOuNbOlhPlWr05bQIlhIeq5FhHZSzblbA0LEZFckdHi2t23Ab+IWsL1UWvsPcuAr3RIQBoWIiLSqGzK2RoWIiK5IpvGXGeehoWIiOQEzRYiIrmisItrzRYiIpITNCxERHJFYRfXiRvHbNwYbxwiItIkDQsRkVxR2MX14MHhcd26eOMQEZEm1c4Wop5rEclyhV1cDxoUHlVci4hktdphIeq5FpEsp+IaVFyLiGS55GEh7iqwRSR7FXZxPWBAeKyshBpdgS4ikq2sk9X9xVK6FpEsVtjFdXFxuKixpkYXNYqIZLnE0BBNxyci2aywi2uoGxqydm28cYiISJM0Y4iI5AIV18OGhceVK+ONQ0REmqS5rkUkF6i4Hj06PC5bFm8cIiLSpE7F4U9WTZWGhYhI9lJxnSiuly6NNQwREWlap+5Rcb1TxbWIZC8V16Wl4VE91yIiWa2oWxGg4lpEspuKa/Vci4jkBPVci0gu+P/s3XecXGW9x/HPb3vf7CabXoEkQAKhJFQlAQIoKCLlImJBr0RFFLuCXMWGBbhXQfQalauiIgIqRZCi9J4AgXRCCpCy2bStyZbkuX88Z3ZnJ1tmdmf3zOx836/Xec3MOWfO+Z0pv/3tM895jorrSMu1imsRkZSWVej/ZO1t2htyJCIi3VNxPWEC5Ob60UIaGsKORkREupFdpG4hIpL6VFzn5MD06f7+ypXhxiIiIt1Sy7WIpAMV1wCHHupvly8PNw4REelWpLhWy7WIpDIV1wAzZvjbZcvCjUNERLrV3i2kScW1iKQuFdfQ0XKt4lpEJGW1dwvZrW4hIpK6VFwDzJzpb5csCTcOERHplk5oFJF0oOIaYNo0KCuDt9+GTZvCjkZERLoQGed6b6NarkUkdQ16cW1mp5nZ82Y2OWpeoZndZGbXmtkvzKwsnmVJk5UFc+b4+y+8kPTNi4ikq1TK2TllOQDsrVNxLSKpa1CLazMbBZQAx8Qs+jnwlHPuKuBu4JdxLkueY4/1t88/PyCbFxFJN6mWs3OG+eK6bVdbsjYpIpJ0g1pcO+eq8cm2nZmNBT4I3BvMegQ428wm97Qs6cEdE/ztUHEtIgKkXs7OLvd9rttqVVyLSOrKGewdOuf2mVn0rHnANudcU7C8zczWAXOB1h6WrU9qYMcf72+ffRb27IGCgqRuXkQkHaVSzlbLtYikg1Q4oXEcsD1mXj0wtpdlndTU1DB79uz2aeHChYlFMXIkzJrlC+tnnknsuSIiCVq4cGF7vgJGhB1PAkLL2SquRSQsieTsQW+57oIDdsfMy8e3gNDLsnZVVVUsWrSof5HMn++H43v4YTjllP5tS0SkBwsWLGDBggUAmNm2kMNJRGg5O6dcxbWIhCORnJ0KLdcbgYqYeaXApl6WJd/8+f724YcHZPMiIkNAaDlbLdcikg5Sobh+FBhrZoUAZpYLTAQe72VZ8p10ku9rvXixH/NaRERihZaz21uua9twziVjkyIiSRfGONeRM2MMwDm3BbgPODWYPx+4yzm3sadlAxJcURG8+93+/l13DcguRETSSSrl7Ky8LLKKs2CvRgwRkdQ12ONclwKfCh5eYmZVwf1PA+8zs28AZwILop7W07Lku+ACf3vHHQO6GxGRVJeKOTtvdB4ArdX7deMWEUkJg3pCo3OuHvhFMEXP3wlc2s1zul02IN7zHt815OmnYd06mDJl0HYtIpJKUjFn543KY88be2ipbqFoetFA7UZEpM9Soc91aikthfPO8/d/85twYxERkU4iLdctW1pCjkREpGsqrrsSDLXCLbdAq356FBFJFXmjguK6WsW1iKQmFdddeec74eCDYfNmuPvu3tcXEZFBoZZrEUl1Kq67Ygaf+Yy//8MfgoZ8EhFJCe0t1yquRSRFqbjuzsc/DlVVfszrRx4JOxoREQHyxvniuvnt5pAjERHpmorr7hQVwec/7+9fe224sYiICACFUwoB2LNuT8iRiIh0TcV1Ty67DMrL4bHH4KGHwo5GRCTjFUwuAGDPhj24feqyJyKpR8V1T4YNg6uu8ve/8hXYuzfceEREMlx2cTa5Vbm4FkfzJnUNEZHUo+K6N5/9LEyYAK++Cr//fdjRiIhkvIIpQeu1uoaISApScd2bwsKOPtdf/Sps2xZuPCIiGa69uF6r4lpEUo+K63hcfDGcfLIvrL/4xbCjERHJaJHLnjetbAo5EhGR/am4jocZ/PKXUFAAt94KDzwQdkQiIhmreGYxAI1LG0OORERkfyqu4zV1Knz72/7+JZfAli2hhiMikqnai+tlKq5FJPWouE7El77ku4ds3Qof/jDs2xd2RCIiGafwoEIsz9izbg9tDW1hhyMi0omK60RkZ8Mf/gAjRvirNv7wh2FHJCKScbJysyg6OOh3vUz9rkUktai4TtTYsfC73/n7V18N994bbjwiIhmo5IgSAOoX1YcciYhIZ0kprs3szGRsJ22ceSZ873vgHHzwg7B0adgRiYjEbSjk7LLjygCoe64u5EhERDrLSWRlM3sXcDUwmo7CPAsYBRQmN7QUd9VVvqj+85/hrLPgqaf8xWZERFLEUM7ZKq5FJFUlVFwDtwI/AF4GImfzGfDBZAaVFszglltg/Xp47jk47TR44gkYOTLsyEREIoZszi4+rJisoix2r9lNS00LeVV5YYckIgIkXlyvdc79d+xMM1uZpHjSS2Eh3H8/zJ0Lr70GZ5wB//oXVFaGHZmICAzhnJ2Vk0XpnFJqH6+l7pk6RrxvRNghiYgAife5/rSZnd/F/LOTEUxaqqiAhx7y42C/8gqceqoukS4iqWJI5+xh84YBsPORnSFHIiLSIdHi+kTg12a2N2raB/yiv4GY2SFmdrOZfcHMfm5mRwTzC83sJjO71sx+YWZl/d1X0o0eDf/+d0eBfcopfixsEZFwDemcXXmG/5Vwx0M7BmoXIiIJS7RbyHeABcAiYG8wz4APJyGWW4H3Oec2mtlE4EHgEODnwD+dc7cHJ+f8ErgoCftLrvHj4fHHfWH92mswb57vIjJmTNiRiUjmGtI5u3ROKdnl2exevZvd63dTODmtz9EUkSEi0ZbrB4C7nXNrnXMbgmk98KskxHIIUBrc3w2Um9lY/Ik3kcGkHwHONrPJSdhf8o0ZA489BjNmwIoVvsDeuDHsqEQkcw3pnJ2Vk0XFqRUA7HxQXUNEJDUkWlz/GbjSzE6KmuYCVyYhljvxP1+WAh8CPgvMA7Y555oAnHNtwDpgbhL2NzBGjYJHH4VZs2D1an+y45tvhh2ViGSmIZ+zh585HICav9UM1C5ERBKSaHH9X8Dngd8Dv4uaPpmEWC4D9gAvAg3OubuAccD2mPXqgbFJ2N/AqaryfbCPPhreeMMX2OvXhx2ViGSeIZ+zh79vOGTDrn/tonVH60DtRkQkbokW19cBE5xzk51zU4JpMnBuEmIpBG4D/gT8xMzmAw7/c2O0fGC/DFpTU8Ps2bPbp4ULFyYhpH6orIRHHoFjj/WF9UknwbJl4cYkIilj4cKF7fkKGKhx5IZ8zs4bkUfFyRW4Nse2uzVSk4gMjERydqInNL7lnKvvYv7oBLfTlT8AH3DO7TIzh0/anwcqYtYrBTbFPrmqqopFixYlIYwkGjbMD9N35pnw9NNw/PFw223+io4iktEWLFjAggULADCzgaoKMyJnV11Qxc5HdlLzlxrGfEwnkYtI8iWSsxNtub4sdoaZXQD8JMHtxG5jBDDLObcrmPU9oA7YAIw1s8JgvVxgIvB4f/Y3qMrKfIF94YVQXw/vfS/ccAM4F3ZkIjL0ZUTOHvH+EViuseOhHTRvbB7IXYmI9CrR4rrSzP4LwMyKzOzXwK/xfe76Ywewx8zGRc3bDrwK3AecGsybD9zlnEuvITiKinyL9be/7YvqL38ZLrgAdu3q/bkiIn2XETk7ryrPX6FxH2z57ZaB3JWISK8S6hbinDvLzI4zs5uAM4AtwCygtj9BOOf2mdk5wDfNbDEwCviKc67OzD4N/NjMZuF/ylzQn32Fxgy++U0/TN/HPgZ33QUvvQR/+Qv4/jsiIkmVSTl7zCfGUHNnDZt/s5mJV07Esmwwdisisp8ei2szO6mbRRuA14HrgUnAOcAX+hOIc24JXZzB7pzbCVzan22nlPPO88P0XXihL65POAF+9CO44grISvSHBBGRDpmcsyvmV5A/MZ896/aw8+Gd7VdvFBEZbOZ66PtrZs/hLxSwC9jX3WrAGOdcfvLDi9/s2bNdyp3Q2JPmZt895Gc/849PPx1++1td0VEkA5nZYudcv3/CyvScveHaDaz7xjoqTqtg1kOzkrptEZGI3nJ2b02l3wGmOucmRQ3jFDtNBi5OZtAZIT8fbroJ/vY3GD7cn/R4+OFw7729P1dEpGsZnbPHfmosWcVZ7Hx4Jw1LGsIOR0QyVI/FtXPufufc1t424py7M3khZZhzzoFXX4X582HbNjj7bN9FpFlnvItIYjI9Z+dW5jLmP/2vf29d/1bI0YhIplIn31Qwdiw8+CBcdx3k5sKNN/oxsVesCDsyEZG0Mv4L4yEbqm+rpml1U9jhiEgGUnGdKrKyfB/sp5+GAw6Al1+GI46AH/wAWnVJXxGReBROLmT0JaNhL6z/1vqwwxGRDKTiOtXMmeML60svhZYWuOoqOO44WLIk7MhERNLC5G9NxvKMrX/eqr7XIjLo+l1cm9mRZvYHM/uSmU3uf0hCWRksXAgPPwyTJ/sh+2bPhm99S32xRaRfMiFnF0woYNxl/vo2a69cG3I0IpJpktFy/SlgvnPuBuCEyGVvJQnmz4fXXoPPfhba2uA73/FdRZ54IuzIRCR9ZUTOnnjlRLJLs9nxwA623bct7HBEJIMko7i+HDg6uP9nIDcJ25SIkhJ/guMTT8D06bByJcyd67uN7NwZdnQikn4yImfnjcxj8ncmA7DmijXs3bM33IBEJGMkXFyb2bvN7OLg/jH4MVU3gr8krnOuLskxCsA73+n7XX/rW5CXB7/+NRx8MNx+e9iRiUgKy+ScPe7ycRTPLGbP2j289SMNzScigyOh4trMvg/cAbwfwDn3AvARM/vgAMQmsfLz4Zpr4JVXfLG9dSt84APwhS/4biMiIlEyPWdn5WQx9eapgL96Y+OyxpAjEpFMkGjL9QnAeODJqHkPAD9MWkTSu0MOgcceg5tv9uNi/+QnfkSRxYvDjkxEUkvG5+xhJw1jzCfG4FocKy9Zyb627q4KLyKSHIkW1w8753YBLmre+UmMR+KVlQWXXQaPPALjx/vC+phj4DOfge3bw45ORFKDcjZw4A0Hkj8hn/pF9bx1nbqHiMjASrS4XmtmVwHTzewDZnYX8BngR8kPTeJy0kn+So5f+hKYwc9/DtOm+VZtdRURyXTK2UBOWQ7Tfz0dgPXXrKfhVY19LSIDJ6Hi2jn3Z+DfQDbwEWAXcLpz7uYBiE3iVVIC11/v+2Kfeirs2AGXXw5HHgn//nfY0YlISJSzO1SeXsmYBb57yPILl9PWoMYHERkYCY8W4px7zjn3Kefcmc65/wQeNbNTByA2SdTMmf7CM3/7G0yZAkuX+mL7vPNg3bqwoxOREChndzjofw6iaEYRTSubeP2y13HO9f4kEZEE5fS00MzuAUp62cYIYBvwr2QFJf1gBuecA+96F/zP/8D3vw9//Sv84x/w+c/DlVdCeXnYUYrIAFDO7ll2UTYz/jKDxXMWU31rNcNOHsaYj40JOywRGWJ6a7l+FXgBuBX4HZAHLAvuR6a/AksHMEbpi4ICX0ivWgUf+pC/bPqPfgRTp8K994YdnYgMDOXsXhQfWsy0n08DYPWnV1P3wpAd5ltEQtJjyzVwHWDB2eaY2YHOuW/GrmRmdw9EcJIE48bBrbfC5z4HX/wiPPUUnH22PwHy+9/3Y2eLyFChnB2H0R8dTe2ztWz+5WaWnrOUo188mvxxyoUikhw9tlw752ojSTowysyyo9cxs3nAMQMQmyTTnDnw+OPw4x9DdjbccAMcdRQ891zYkYlIkihnx2/qTVMZNm8YLZtbWHrOUvY26fLoIpIciZ7QeCewxMz+18yuM7P7gIeBG5MfmiRdVhZ85Svw5JN+uL7ly+H44+E//9Nf7VFEhhrl7G5k5WYx484ZFBxQQP2ielZ+bCVun05wFJH+S3QovoeB0/Enw0zHD+t0vnPuBwMQmwyU44/3w/Z9/ev+Co+33OKL7Z/8BFpbw45ORJJEObtnucNzOeyew8guzabmLzWs+eIajSAiIv1m/U0kZpYFHOmcS8q1t83sOOAdwBrgKaAR+DFQD1QAX3PO7XcGyuzZs92iRYuSEUJmWb3ajyLywAP+8aGHwo03+iH8RGRQmNli59zsQdqXcnaMnf/ayatnvoprcUy5dgqTrpwUdkgiksJ6y9m9ndAYu7FbuphdBewFzkkwtq62/zHgQOfc1VHz/g/4p3PudjN7F/BL4KL+7ksC06b5YfoiQ/UtXw7z58O55/p+2ZMnhx2hiPSRcnZ8Kk6t4JBbD2H5B5az7qp15I3MY8x/aog+EembRPtcTwXeBDZETflAv1tAzGwucDHwX1HzxgIfBCJjxz0CnG1mk/u7P4liBu95DyxbBtdeC0VFfmzsQw6Bb38bdu8OO0IR6Rvl7DiN/I+RTL1pKgCrFqyi+rbqkCMSkXSVUMs18Ann3KroGWZ2MHBJEmL5b+AZ4CYzOwj4LjAJ2OacawJwzrWZ2TpgLrA+CfuUaPn5fmzsD38YvvpVuO02uOYa+MUv4FOfgk98AsaPDztKEYmfcnYCxn1mHK3bW1n/rfWs+NAKMBj1gVFhhyUiaSbRExpXdTG7GfhEf4Iws+nAkcCvnHOX468c9k9gArA9ZvV6YGzsNmpqapg9e3b7tHDhwv6ElNnGj4c//QmeeMIP11dd7VuwJ06EefNg4ULYHvu2iEgiFi5c2J6v8FdNTDrl7MRN/uZkJn1rEuyDFRevoPrPasEWkcRydkInNAYtENFPyAZGAbc65y7tQ6yR7Z4LLHTOjQgeFwI1wHeA85xzx0at+xLwJ+fc9dHbSKWTY4YU5/z42Dff7K/s2Nzs5+fk+EusX3QRvPe9UFoabpwiaWygTmhUzu67ddesY8O3N0AWHHzLwYz+6OiwQxKRFJHUExrx46P+iY5kvQ+ods6t7mN80XG0x+Kc221mr+P/EFTErFsKbOrn/iReZr61et48qK2Fv//dt2o/8gjcd5+fCgrgzDPhggt83+2SkrCjFhFPObuPplwzBYAN397AyktW0rypmYlfn4iZhRyZiKS6RE9o/INz7jHn3OPB9CRgQX+7/ngVKDez6M5tbfiEPDZoFcHMcoGJwOP93J/0RXk5fPSj8OCDsGkT3HQTvOMdsGePPwHyoougqgrOOw9uvx0aGsKOWCTTKWf3w5RrpnDQjQeBwbqr1vH6Z1/H7dU42CLSs0SL63d1Ma8auKM/QTjnVgIPAOcDmNkwIA/4I3AfEBl0eT5wl3NuY3/2J0kwahRcfrm/2uNbb/kL0Jx4Ykeh/YEPwMiR8P73w//9H9TUhB2xSCZSzu6n8Z8dz6G3H4rlGZtu3sTSc5fSVtcWdlgiksLi6nNtZpcD38P/vNeVh51zXSXx+AMxGwH8FFiEb+n4lXNuuZlV4C9IsB4YDVzpnNuvSTTV++9ljLffhjvvhDvugGee6ZhvBiecAGef7afp0/08EUl6n2vl7OTb9fgulp6zlLZdbRQdUsTMu2dSNLUo7LBEJAS95ey4T2g0s8OBjwI3xixqcc5t7nuIyZFuiTojvP22Pwnynnvg3/+GlpaOZVOn+iL7fe/zl2PPSbT7v8jQMRAnNCpnJ1/TmiaWvm8pTcubyC7P5tDbDmX4u4eHHZaIDLKkFdfBxnKdc61dzJ/onHuzjzEmRTom6oxSXw8PPQR33+2vBrljR8ey4cPhtNPgpJNg7lx/8Rq1aksGGcDRQpSzk6ytvo2VH1nJtr9vA4MJX53AlO9OISs30V6WIpKu+lVcm9mXgfudc8uDx9/Anw3eaTXgZOfcvP6H23fpmqgzUlub7zJyzz2+2F6zpvPyESN8oX3CCXDccXD00X5EEpEhKlnFtXL24HD7HG/+4E3WfXMd7IPSY0s59LZDKZxSGHZoIjII+jsU3xRgWNTjo4BiYAsdQzsZvl+dSHxycnzxfNJJcN11sGoVPPqov2jN44/D5s3+pMi//tWvn5sLRxzhC+3INGWKWrdF9qecPQgsy5j0jUmUzy1nxQdXUP98PYuOWMRBPz2I0R8dreH6RDJcot1CjgKWOudaYuYfGmkpCUs6t4JIFOfgjTd8of3cc35autTPj1ZV5Yvs44/3t7Nn60I2krYGsFuIcvYAa93RyqpPrGLb37YBUHFGBdMXTqdgon5tExmqkn0RmcOccy/F7GA6cCAQaqKWIcIMDjrITx//uJ9XXw8vvthRbD/7rB/a7957/QSQlQUzZ8Kxx8KcOb4rycyZkJcX3rGIhE85e4DlVuYy464ZVP+hmjVXrGHngzt5ccaLHPCjAxj7ybFYtlqxRTJNvEPxvRP/U+ICYGHM4krgRufcxOSHF7+h0goicXAO1q7tKLafew5eecX35Y6WlweHHeYL7aOO8reHHQb5+eHELdKNARiKTzk7BM1bmnn98tfZdpdvxS45uoSpN02l/PjykCMTkWRKVsv1LuD3+P58J8YsawF+3bfwRPrADA480E8XX+zn7d4NL70Ezz8Pixf7adWqjvsRubm+Rfvoozumww7TCZMy1ChnhyB/dD4z75xJzV01rPn8GhoWN/DyCS8z6iOjOOBHB5A/Wv/Yi2SCRMa5LgcucM6lZFIeiq0g0k91dfDyyx0F9uLFsHr1/v23s7P9RW2OOAJmzeqYRuucLxkcAzTOtXJ2iPY27mXDtRt46/q3cC2O7NJsJnx1AuM/P56cEo3rL5LOkjrOdQ87eYdz7ql+b6gfhnqiliSpr/cF90svdW7h3rdv/3VHjuwotCOF9/TpvvVbJIkG6oTGHvannD1ImtY08cYX3mD7fdsByB2Zy6T/msTYBWPJytPY2CLpqL/jXN8DlPS2D2CGc25k30JMjkxJ1DIAmppg2TJYsqTzVFe3/7p5eTBjRucW7lmzoLJy8OOWISOJ41wrZ6eonY/uZO2Va6l/vh6AgskFTL5mMiM/OFIXoBFJM/3tc70K33dvIx1jpO63D+A9fQtPJAUUFfkRRubM6ZjnHKxfv3/BvXatb/l++eXO25gwYf+C+8ADfZcTkcGjnJ2iKk6u4Khnj2Lb3dtY9411NC1vYuUlK1l/zXomfG0Coy8ZTXaB8oXIUNBby/UIoN4519zjRsyGO+e2Jzu4RGRaK4iEpK4OXnvNj04SKbhfe82fUBkrP99fyn3GDH8S5cyZ/v6kSX7oQJFAEluulbPTgNvrqP5jNRu+v4Hdq33uyBuTx4QvTWDMJ8eoT7ZIihusPtefCPukmUxO1BKyvXv9JdyjW7hfew3efLPr9YuL4dBDO4rtSOE9dqyuOpmhQuhzrZydAtxeR81fa9jw/Q00LmkEIGdYDmMuHcO4y8fpQjQiKWpAi2szmwh8BrjcOVfc5w0lgRK1pJzaWli+3PfnXrq043bLlq7XLy/vXHDPmOGL8FGjVHQPcYNVXCtnpybnHDvu38GGazdQ90xwrkcWVJ1bxbgrxlF+YrkuqS6SQpJeXJtZNnA2cClwGlAdbGdcfwLtLyVqSRvbt+9fcC9dCjt2dL3+sGFw8MG+i0lkOvhgmDJFfbqHiIEsrpWz00vd83W8/dO3qbmjBtfm/z6XHF3CuE+Po+rCKnUZEUkBSSuuzWwK8AngY8BwYAPwOeBB4Hjn3DP9D7fvlKglrTkH1dX7F90rVsCuXV0/Jz8fpk3bv/CeNg0KCwc3fumXARrnWjk7jTVvbGbjLzay6X830bbdX302uySbkR8cyZhLx1B6dKlas0VC0u/i2swuwLd4nAKswV/Z63fAVc65LyQx1n5RopYhyTnYutUX2ZFp5Up/+/bbXT/HDCZP7tzKfeihvptJuS7DnIqSWVwrZw8te3fvZevtW9n8q80dXUaAkiNKGHPpGEZeNJLcCo29LzKYknH58934y+U+A3zCObc62HD/z4QUkZ6Z+T7Xo0bBvHmdl9XX+0I7UmxHpjVrYN06P91/f+fnjBvX0Zd7xoyO+yq6hxLl7CEkuzCbMZeMYcwlY2hc1sjmX29my++30PBKA69/5nXWfGENw98znFEfGsXwM4eTla+RiETClki3kIn41pDJwN3AXOfcZ4Nl5pIx7Eg/qBVEJNDSAm+80dHKvXy5n1asgD17un7O+PGdi+3IbVnZ4MaeoQaoW4hy9hC1d89etv19G1t+s4Wd/9rZPqJ5zrAcqv6jilEfGuVPgsxStxGRgTAQJzTmAOcAFwNLgd8ApzvnFvYn0P5Sohbpxd69/sI4y5Z1nlau7L7onjCho+iOFNyHHgqlpYMa+lA3wCc0KmcPYc0bm9n6561U/6Gahlca2ufnT8pn5IUjqTq/itLZ6p8tkkwDPRTfQcBngUudc0V93lASKFGL9NHevf7Kk8uWdQwdGCm6m7u5FsnEiV0X3SW9XXlbujKIQ/EpZw9hDUsb2PrHrVT/sZrmtzq+u/mT8qlXEitgAAAgAElEQVQ6r4qqC6ooO6ZMLdoi/TRYF5H5iHPu9/3eUD8oUYskWVtbR9EdmZYv90V3S0vXz5k0qaPgjozVfcgh/hLz0q0QLiKjnD2EuX2O2qdqqbmzhpq7amjZ1PF9zR+fz4jzRlB1XhVlx5eRlaM+2iKJGpTiOlnMLA94EbjCOfeYmRUCPwbqgQrga865uq6eq0QtMkja2nyf7tjuJatWQWvr/uubwYEHdr4E/MyZfsjAvLzBjz8FDXZxnUx9zdvK2YPD7XPUPVvnC+07a2h+u6NFO2d4DsPPHM7w9w6n8oxKcso0hrZIPJIxWshg+gr+5JuInwP/dM7dbmbvAn4JXBRGYCISyMmB6dP9dO65HfNbW/cvupcu9UX3mjV++vvf999ObNF9wAG6OE56Ud5OYZZllJ9YTvmJ5Rx4w4HUv1jP1ju2sv3u7exes5vqW6upvrUayzWGzR3G8Pf6YrtwisbKF+mrlGm5NrMTgIOBbwKXAKuBdUCFc64pOCmnFpjhnFsf+3y1goikqOZmWL16/ytSrl3rx/GOVVDg+29HF9wzZ/qTK4foSVnp2nLdn7ytnB0u5xxNq5rYfu92tt+zndpnamFfx/KiGUUMP9O3aJe/o1xD/IlESYuWazMrBs53zn3RzL4ZzJ4HbHPONQE459rMbB0wF1gfSqAikrj8fDjsMD9Fa2z0wwNGF9xLl/qL47z0kp+ilZZ2FNyHHw5z5sARR/hiXAad8nZ6MzOKDy6m+OBiJn5lIq3bW9l+/3a237udHf/cQdOyJpqWNfHWdW+RVZTFsLnDqDyjkoozKiiaXqTRR0R6kBLFNfA14Icx88YB22Pm1QNju9pATU0Ns2d3/BOxYMECFixYkMwYRSSZioth9mw/Rdu1y584GV1wL10KNTXw7LN+isjJ8YX2Mcf4YvuYY/wJlGnQrWThwoUsXNg+Gt6IMGPpo37lbeXs1JI7PJfRHx7N6A+PZl/LPmqfrGXHgzvY8eAOGl9tZMcDO9jxwA4A8ifk+0L79Aoq5lfoCpGSERLJ2aF3CzGzM4B859w9weP1+J8XZwMXOOeOjVr3JeBPzrnrY7ejnxhFhritW30r92uvwSuvwIsv+sexOay4GI46qnPBPXlySncpSbduIcnI28rZ6aN5czM7H97Jjgd3sPPhnbTWRJ24nAWlR5Uy7ORhDDtlGOXvKCenJFXa7UQGRjp0C/kycGTUT0wV+KuJ/XdwP1opsGnwQhORlDFypJ9OPrljXkOD7z7ywgu+2H7hBX+hnCef9FPE6NFw4okd05FHQq5a2/pBeTuD5I/JZ/RHRjP6I6Nx+xwNrzSw46Ed7HxwJ7VP11K/qJ76RfW8dd1bWI5ROicotk8eRvkJ5WQXpf4vSSLJlAot11VAftSsZ4EvAs8Dy4Eq59xuM8sFGoADnHMbY7ejVhARAXz3kRdf7Ci2n38etsf0VCgshGOP7Si2TzgBysvDiZe0bLnud95Wzh4a2hraqHu6jp2P7mTXo7uoX1Tf6cRIyzXKjitrL7bLjisju0DFtqS3lG+5ds7VRD82s71AjXPuTTO7DzgVuA+YD9zVVWEtItKuqgrOPNNP4LuNrF4NTz8NTz3lb1evhsce8xNAVpbvQjJ/vp+OP96fiCldUt6WiJySHCrPqKTyjEoA2uraqH2ytr3Ybni5gdona6l9spYN39mA5Rmls0spf0e5n04sJ7dSvyLJ0BJ6y3WsSN+94GIEFfiLEawHRgNXOucaunqeWkFEJG41NfDMMx0F94sv+ovjRBQWwkknwemnw/vfD1OmDGg46dZyHasveVs5OzO07mxl1+O72PWonxqXNkJM2VE0o6ij2H5HOQWTCjQaiaS0tLpCY38oUYtIn9XX+z7ajzzip9de67z8iCP8BXPOPdcPBZhk6V5c94VydmZq3dlK3bN11D5VS+1TtdS9UIdr7lyH5I3L69SyXXxYsS7TLilFxbWISKK2bIF//Qvuu89PDVENr4cfDh/9KFx8MYwalZTdqbiWTLWveR/1i+vbi+3ap2pp29nWaZ2soixKjy6l7LgyPx1bRv44dduS8Ki4FhHpjz17fKH917/C3/4GO3f6+dnZcNZZ8NWv+pMi+0HFtYjn9jmaVjb5QvvJWmqfrWXPG3v2Wy9/fD5lx5VReqwvukuPKtWoJDJoVFyLiCRLc7Nvyf7d7+D++2HvXj9//nz4xS/goIP6tFkV1yLda6lpof6Feuqeq/PTC3XsrdvbeaVsKJlV0t6yXTqnlKJpRVi2+m5L8qm4FhEZCNXVcPPNcOONUFvrL17zj3/A3LkJb0rFtUj8Iq3bdc/VUfe8L7gblzZ2GgIQILskm5IjSyidXdo+FR5UiGWp4Jb+UXEtIjKQduyAT34S7rwTRoyADRugqCihTai4FumftoY26hf51u365+upX1xP81vN+62XXZZN6dGlfgoK7oIDNDqJJCblx7kWEUlrlZW+m8gDD8C2bbBypb/8uogMmpySHCrmVVAxr+MCoS3VLdQv9oV25CqSLZta2ocFbH/usBxKZ5dScnTQyn1kKQVTCtTCLX2m4lpEpD+eew4uuwwaG2H8eDj00LAjEhEgb1Qew88czvAzh7fPa97c3KnYrl9UT2t1Kzsf2cnOR3a2r5ddlk3JrBJKjiih5Eh/WzyjmKw8DQkovVNxLSKSqOpq37/697+Hxx/386ZMgbvvhoKCcGMTkW7lj8kn/z35jHjPCACcc7RsaukothfX0/BKAy2bW9qvLBlhuUbRoUWUHlnaUXTPKiGnXKWUdKZPhIhITxoa/EVlliyBV16BJ56AFSs6lhcVwRVXwJVXQmlpeHGKSMLMjPxx+eSPy2fE+0a0z2+pbqHhlQbqX/bFdsPLDex+fTeNSxppXNLYaRsFBxR0auEuPbKUvLF56sedwVRci4g4B5s2werVflq1yt+uXAlr1/rl0YqK/NjWF14I558P5eXhxC0iAyJvVB6VZ1RSeUZl+7y2hjYaX2uk4eWG9oK74bUG9qzdw561e9j2123t6+aOyKV4VrHvWnJ4CcWHF1N0SBHZBRqLOxOouBaRzNDSAm++CevXw7p1/nbt2o5CurGx6+fl5vp+1LNm+enYY2HOHMjLG8zoRSRkOSU5lB9fTvnxHf9M72vbR9PKpo5iO7ht3dbKrn/tYte/Ok6cJBuKDi5qL7ZLDvfdStTKPfSouBaRoaGxETZuhLff7lxERwrpjRth377unz9iBEyb5qfp0zvuT5umQlpEupSVk0XJzBJKZpbAh/w85xzNbzbT8FoDjUsaaXi1gcZXG2la3UTTMj9xW8c2cipzOgruWf62eEYx2YVq5U5XKq5FJLU5B9u3dxTO3d3W1va8nawsmDjRn3g4ebK/nTKlo4CurOz5+SIicTAzCiYVUDCpoP3ESYC9TXtpXN5I46tBwR0U3m072tj12C52PRbVyp0FhVML21u3Iy3d+RPz1cqdBlRci0g4Ghr8qBtbtvipu/tbtvguHb3Jz/dD4Y0b528jxXOkkJ4wwXfxEBEJQXZRNmWzyyibXdY+LzJaScOrDTQsaWgvvJtWNrF71W52r9pNzR01Hdsoz+7UraT48GJKDishu1it3KlExbWIJEdbm29hrqnxF1Opqem4HymWo4vm7vo4d2XYsI6iubvbykpQi46IpJHo0UqGv7tjPO69e/bStKKpvdhuWOJbulu3te43RCAGhQcWdhTcs/xtwWRdCCcsKq5FpGuNjV0Xyt3N27mz921Gy8+H0aP3n0aN2v9xgpcTFxFJZ9kF2ZQe6a8WGeGco6W6xRfc0a3cK5rYvWY3u9fs7jRiSXZpNsWHBSOWRPpyH1ZMTolKv4GmV1hkKHMOmppgx479p+3bu54fWbZnT2L7MvOtx1VVfhoxouN21Kj9i+ayMrU0i4jEyczIH51P/uh8Kk/vOEdkX0swYklw4mSk8G7Z0kLdM3XUPVPXaTsFBxZ0KrhLZgWt3MrHSaPiWiQd7N0LdXW+dXjXLn8bT5G8Ywc0N/dtnwUF+xfJXRXOkdvKSshWvz8RkcGUlZflT3w8vKTT/JatLR0t3EsafCv38ib2vLGHPW/s2b+V+/DOrdzqy913Kq5FBkOkBTlSHEcK5Nj73c2rq+t9H90pLPSFb/Q0fPj+82KnoiK1LIuIpKm8kXlUnlZJ5WndtHIv6Si6W6tbqXu6jrqno/7WRPpyR18MZ1YxBZPUyt0bFdci8Wpt7b0Q7m75rl3++X1l5q8COGyYnyoqei6QI8sqKnxxLSIiGa9TK/eHOua3VLe0F9qRIQKblkf15b4rqpW7LLvjxMlI0T2zWK3cUVRcS+ZobfVjIUcXvLFTbEEcPTU19W//hYW+2I0ukKNve5pXWqouFyIiMiDyRuVReXrl/n25VzR1HiZwSQOtW1upfaqW2qdiRiw5qNAX20eUUHJUCaVHlZI3KjMvwJUyxbWZnQvcAJQBfwC+5JxrM7MRwPeBXUAu8DXnXD+aACVtxVMc9zQlMvRbV7Kz+1YYR+7n5yfndRBJAcrZIkNbVl5Wex9sPtwxv3lLc0c/7qDoblrRxO7Xd7P79d3U3NkxLnfeuDxKjyql9OjSjoI7Ay73nhLFtZlNBM4BzgdmAL8A3gKuB+4Avu6ce97MPgV8B7gyrFilH8IujrOyOgrd3qboojgyFRerD7IIytkimazLEUuag77cSxpoeLmB+pfqaXi5gZaNLWzfuJ3t925vXzd3VG7ngvvoUvInDK0rT6ZEcQ1MAj7unGsDFpvZTOBkM3sKmO6cez5Y715glZl9zznXz0pLEtbW1r/iuKGhf/tPpDjuaiopUXEskhzK2SLSLis/qpX7I36e2+fYvWY39YvraXipgfrF9dS/VE9rdSs7HtjBjgd2tD8/d0Rue8t26bGllB1bRv6Y9P21NyWKa+fckzGzNgKlwCnAuqj1Ngb/2cwGHh+0AIeKVCiOY0/KU3EsknaUs0WkN5ZlFE0romhaEaMuGgX4C+HsWbuH+pfqOxXdrdta2fnQTnY+1HExsvyJ+ZQdV+anY8soOaqE7IL0OPcoJYrrLswBrgMWANtjltUDYwc9olTQ1uaHZIv3BLxkF8dm/W85zspKzmshIqlEOVtEemVmFB5YSOGBhYy8YCTgC+7mN5t9wf1iPXXP11H/Yj3NbzZT82YNNX/xfbgt1yiZVdJRcB9XRsEBqTksYMoV12Z2EFDjnFtiZg7YHbNKPrDfyTE1NTXMnj27/fGCBQtYsGDBgMbaZ875Ijn2UtLbtnW+H1sc19f3b78qjkVSxsKFC1m4cGHk4YgwY+mPjMjZIjJgzIyCSQUUTCqg6v1VALi9jsYVjdQ9V0f98/XUPVdH47JG6hfVU7+ono0/2whA3tg8hs0dxrC5wyifW07R9KIBK7YTydnmnBuQIPrCzLKBHwJXO+eazexK4GTn3OnBcgP2BPOeiX7u7Nmz3aJFiwY95na7d8PmzX7atAmqq7svmrdt69uYx7FjHScyVVSoOBZJUWa22Dk3u/c1U0ta52wRSSttdW3UL/It23XP1VH7dC1t29s6rZM7KpdhJ/liu+K0CoqmFQ1ILL3l7FRruf4ycL1zLnK95oeBT0Utnwy0AIObkffsgXXrYM0aWL/eF8+bNnUU0ps2+e4ZiSgp6Xzp6Nj7kctJRxfIpaUqjkUklaRmzhaRISenLIeKUyqoOKUC8CdMNq1oYtfju9qn1upWau6ooeYO35Wk8KBCKs+qZPhZwxl20jCy8genhkqZ4trMrsIn4EIzOwCYBzwBbDazmc65pcC7gJuccy0DEoRz8MYb8Pzzflq61BfUb7/tl/UkNxfGjPHT2LEwalRHsRx7O2IEFBQMyCGIiAyGlMjZIpKxLMsonlFM8Yxixl02Duccu1fv9oX2o7vY8dAOdq/ZzcafbmTjTzeSU5HD6I+OZuxlYymaOjAt2hEpUVyb2dXAd2Nmr3TO3WJmFwLfNLP1QAnwjaQH0NAAN94Iv/89rFq1//LsbJg8GQ48EKZMgfHjO4royDR8uFqVRSQjhJ6zRURimBlF04soml7E2AVjcXsddc/Vsf0ffpztxqWNvP2Tt9n4s41MunoSk745acD6Z6dUn+v+6HP/veXL4V3vgrfe8o9HjIATT4RjjoGjjoKpU2HiRN8yLSIyANK1z3V/qM+1iAym+sX1vH3T21T/rhqASVdPYsp3p/RpW+nW53pwNTbCWWf5wvqoo+D734f58yEns18WERERkaGk9OhSDvntIYy8YCSvvfc1NvxgA6M/NprCAwqTvq/M7sfwu9/5ExRnzYKnnvIt2CqsRURERIak4WcNZ+RFI2EvVP+xekD2kdnF9V13+duvfAUKk/+fi4iIiIiklqrz/HjatU/UDsj2M7e43rcPXnzR3z/llHBjEREREZFBUTqnFID6l/p5cb5uZG5xvXWrv+Lh8OF+5A8RERERGfLyx+dj+UbbjjbaGtp6f0KCMre43ugvncm4ceHGISIiIiKDxszIH58PQPNbzb2snbjMLa63bPG3arUWERERySj543xx3bIl+de4ytziujboxD5sWLhxiIiIiMigyqnwo8O17VK3kOSpDzqxl5aGG4eIiIiIDKqcYSquky9SXJeUhBuHiIiIiAwqFdcDIS/P97euqgo7EhEREREZRHkj88gbk4dlW9K3nbmXI/zc5/wkIiIiIhll0lWTmHTVpAHZdua2XIuIiIiIJJmKaxERERGRJFFxLSIiIiKSJCquRURERESSJKOL64ULF4YdQr8o/nAp/vCl+zGke/xhGAqv2VA4BtBxpJKhcAwwdI5DxXUaU/zhUvzhS/djSPf4wzAUXrOhcAyg40glQ+EYYOgcR0YX1yIiIiIiyWTOubBjSAozqwE2JPi0EcC2AQhnsCj+cCn+8KX7MUTin+Scy6grWvUxZ0P6v+cwNI4BdBypZCgcA6TPcfSYs4dMcS0iIiIiEjZ1CxERERERSRIV1yIiIiIiSaLiWkREREQkSTKyuDazQjO7ycyuNbNfmFlZ2DElysy+Z2YumJaEHU88zOw0M3vezCZHzUub96Kr+IP5Kf9emNm5ZrbOzLab2U/NLCeYP8LMfmlmPzKz/zaz3LBj7Up38QfL/hD1+t8dZpzdMbO5ZrbKzHaa2U+i5qfN538wxfu6mNm7zexnZvZzMzt/sOPsSTzHYGblZnaHmdWZ2ctmdlwYsfYk0c+omV1oZo8NUnhxSeQYgnU/b2YfTbX3I4HvxVVmdnVw+z0zs8GOtTfd/T2NWSdlv9+9cs5l3AT8H3BhcP9dwG1hx5Rg/MXA/wLzg2lq2DHFEfMo4P2AAyan23vRQ/wp/14AE4HfA0cDHwEagS8Hyx4Fjg3ufwr4QdjxJhj/ROC/o17/CWHH20X85cBVQAVwFtAKzA+WpcXnP4TXrNfXBZgGvARkBY+fAGaGHXuCx3AtcFHw2f4HUA0Uhx17oscRte5Y4FXgsbDj7ssxBN/Re4ADwo65H5+p9wC/inp8S+Q5qTJ19/c0Zp2U/n73eoxhBxDCmzoWaAaKgsc5wR/rLt/gVJyAzwdTQdixJBh3VvSXKd3ei9j40+W9AN4J5EQ9/nHwh/w4YFPU/HFAQwr+ce8y/uD+/wDnRy9PtQkoJBiZKXi8GDg53T7/g/h6xfW64P+pvTbq8TeA/ws7/gSP4ZSo++X4f7zmhB1/oscRtf5PgctIoeI6kWMAHgZODjvmfn6mvgL8OerxzcCCsOPv4nj2+3saszxlv9/xTJnYLWQesM051wTgnGsD1gFzwwwqXsHPOx8AbgCqzeyikEOKm3NuX8yseaTRexEbf7q8F865J4PXNmIj8CZwCv71jqy3Mbg7exDD61V38ZtZEfBe4C/B41NDCbAXzrndLvjrYGbFwArgMdLs8z+I5hHf69Lp8wusCp6bCuYRxzE45/4ddb8WqAPeHrwwezWPOD+jZvafwJ+ApsEMMA7ziOMYzOy9wFTgKDO738x+bGZ5gx1sD+YR33txH3C2mZ1vZhX4caNvHcxA49FFPRArlb/fvcrE4nocsD1mXj3+v8KU57zj8F+Ym4A/mtm7Qw6rr/RehGMOvlUgXV//OcD/OueanHMH4Y/jH8ADZnZYuKF1z8zeCTwA5AJFpO/rP9DifV1i10ul1y7h99bMpuFbfDcPZGAJius4zOxAoNI59/xgBZaAeN+LDwLPAjcCH8Y3nHx/wKOLX1zH4Zxbge9q9Dv8Pzsfc87tHpQIkyuVv9+9ysTi2gGxH7R8/M9xacM5t9M5dzXwA3y3hHSk92KQmdlBQI1zbglp+PrHxA+Ac26zc+5S4A7g8tCC690a4LfAqcB1pOHrP0jifV1i10ul164v7+1n8D/pp5Jej8PMsvCx/88gxpWIeN+LGcDTzrlW59x24FfAxYMQX7wS+UzlAucABwG3p+qJ6r1I5e93rzKxuN6IP2khWimwKYRYkuFG/Eld6UjvxSAys2zgk8DXg1mdXv+gm0sxKfr6dxF/rJR+/YN/Am4Bvoz/KXeoff6TJd7XJXa9VHrtEnpvzWw+vtV67UAHlqB4juNEYAGwxcy24X/FOzG4nwrifS9yginiVaByAONKVFzHYWZz8f32H8bnmSOALw5KhMmVyt/vXmVicf0oMNbMCgGC/+gmAo+HGlXf7cOfUZuO9F4Mri8D1zvnmoPHDwPTo5ZPBlqARYMcV7xi44+V6q9/xCL8H46h9vlPlnhfl9jP77Tguakg7vfWzA7Bj07xt8ENMS7xHMcLwMH4Iu4I4Jv4z/gRgxtqt+J9L16l8+epDX9+RKqI9zguANYCOOc2AVfjTwpPN6n8/e5VxhXXzrkt+A7/kZOf5gN3RZ3MldLMbKyZfdDMsoKWxi/hz6JNeVFjbRqk33sRG386vRdmdhXwMlBoZgeY2ceBXcBmM5sZrPYu4CbnXEtYcXanq/jNbKaZnR0szwM+ih9JJKWYWZGZHRk160zgp+n2+R8sPb0uZvZlM4v8wf1f4Axo/26egv/1InTxHoOZHYBv9X3IzCab2Swz+3Q4Ue8vnuNwzjU7596OTMBOoDm4H7oEPk8/Bd5tZvnB4xNIkc8TJHQcrwDR+WYv/h+glBL79zSYlxbf73hYcBJ7RgnOoP0xsB4YDVzpnGsINag4mdkM/ElRdcBTwE+ccyvDjap3ZlYKfAj4OfAd4GfOuZp0eS+6ih8YSRq8F2Z2NfDdmNkrnXOHmNkkfEvTeqAE+EbMyByh6y5+fP/q2/BnlD8H/ChoqUkpZnYM/oTLlfg4X3DO3REsS4vP/2Dr7nUxs8X4sdjvDNb7EHBU8LTHnXMpcxGh3o4BeBp4HpgQ89QPOOduH8xYexLvexG1/iXAJc65eYMcarcS/Dy9Ez8yRaFzLpVOaIzrOII+8NfgR23Zhv98XdvDL36Drod6IG2+373JyOJaRERERGQgZFy3EBERERGRgaLiWkREREQkSVRci4iIiIgkiYprEREREZEkUXEtIiIiIpIkKq5FRERERJJExbWIiIiISJKouBaRlGJmVWZ2j5ltMLOvhx2PiIj0THm7M11ERkRSipmdATwKlOGvajjNObcj3KhERKQ7ytudqbgWkZRlZk8Ap6TaJdlFRKRrytvqFiIiKcrMJgAPZXKCFhFJJ8rbnoprSToz+6qZ3Rx2HBHJiMfM5pvZs2Y2L0lhDSlmdn2St2fAfwA/TOZ2RaQz5ev+M7PJZvYrM/tlkraXbWaXmtmGONY9y8ye7+uxJjN3K293UHGdZszsEDP7q5n92syWmJkzsw+FHVeM+4Hbwg4iSjLieQuYlYRYUpKZnWlm28xss5kdb2YTgs+ZM7OfmdmwYL1jzGypmd1gZoXBvFHAliSH9EHgN865tmD7ImlH+bpP0jFfbw1u85O0vSxgFzAxjnUXAYd0t9DMZpjZ7cFn72ozGx21LNm5W3k7oOI6/fwd+LZz7hPAkcDvIwvM7HAzmxtGUNH7ds4tdc49FUYcXUkknu5eQ+fcKmBbMuNK5vvV32055+4Hfg60Oeeedc69BVwKtADPO+d2Beu9gE/mX3HO7Q6e/l7gvn4dQBQz+y/gGuARM1sJTE/WtkUGmfJ1gtIxXzvnmoBNydqXc64VWBznutX4Qry75cuAhcHD651z0cV00nK38nZnKq7TiJlVAtOARgDn3D7gm4Azs3J84rYQ4gpt38kUx3HsG8R9hbGt24HxZnY0gHNuO/Av4D1R+xoGbAo+exEHO+dW9nPf7Zxz33XOTXXOzXbOHeyceyJZ2xYZLMrXAysF83WyR4dIJP7e1t0b3Mb2g05a7lbe7kzFdRoJhrVZBfzTzI4L5m0AXsD3c5oCLDCzjwGY2XfM7L+Cn4Ta+1WZ2Wwz+28z+7yZNZnZKjP7opmdY2Yvmtn7zexvZlZnZp+Mel6X24vZ96Vm9tvovmfBT6M3mdkPzOxRMzsmmN/j/qKeb2b2STNbZmanmdmrZlZtZpfEsY+p0fH0ss/9XsMuVAQ/89ab2aUxr+k1wWvz1yAZx/ua7bcvMzvfzF4ws0vM7HUz+2EP2+vqve8ynp4ELRzLgAuiZlcBZ1nQBQQ4G7gnKs5CoCnq8cFm9kMzu97MHjffveRCM1tkZu8zs9vMbJeZXWFmRwfHuN3MzuwtPpF0onydOfm6K90dY0/7MbMpZnaL+Vbg6+LZTxf7vd/M/mhmJ/aynnL3QHLOaUqjCZgBbMD/p3oLMDJq2XpgXnB/JrA2uF+F/696OP4/7zeB2cGym4Fngvt5wA7gZ/ixKi8CNva0vdh94/uc/RL4bTC/AF+wDQsenw/UAKU97a+L4x4d7POqYN3vAXuASb3sIzaeHvcZ/Rp2EcN64AagKHjem8H8MuDPUes9jU+Mcb1m3eyrEv9T30+AY4D5PW0v5r3vMp44P19XA28E9w8EfgA0A+dGfV4sav2zgQEtW6IAAAXLSURBVHcG94uAF4Gi4PETwetQBtQC3wEKgQuDx/+B/wf/SwSfQU2ahtKE8nVG5Otg+TVxvo49vdcvAodEPcfF+TmLvJ/HA9/qYvm8YD85UfOUuwdwUst1mnG+dfEwfLL5CLDMzLo6ceN14MNmlgu8M5hXgv8yTwAi/WVX4L80OOdagHrgTudcHf7LNrqX7cXG10znEyTOAna7oM8u8Ff8l/LcXvYXu93INv8UrPs9fH/g03rZR6d4EtlnN+51vn/di8CYqGMcHrQsfR7/mkKcr1k3x7sDn8Tuds694Jx7JIHtdRdPPG4HDjDfNeQ84FfAg8B/mFkJUO+C7Bs4EXgmuH82/g9YpDXk/cB3g9d5F/Bv5/tpLwLKnHN/cf6n8pdJ7D0QSQvK15mRr7vQ7TH2sJ/5QLlzLhLPogT3OQe42Dn37TjXV+4eQDlhByDxM7NsYIJzbj3wOTP7Df7M6oXAsdHrOueazWwc8DXgd5FNOOe2mtki4FT8f9ZT8AVV+1Oj7rcRdB3qbnvdhBq9jalAblRc+8xsPTC+p/31xjm3x8zWAsOAEQnso8/7jNFGx/dnIrDOOfeT2JXiec3M7MP4IjZiuvM/H7voWBN4D7qNpzfOudfN7GV815By59xaM7sd35p0AfCPqLgjn41If75J+JamyLa2R2866v5eOtuHuqjJEKN83en4MiFfR+v2dezhvTkc/49EX50CTDWzEudcQzfr7AuOQbl7gOlFSS+5RPWHdc4twf/sdnDsimZ2CnCZc+57zo/8EO0zwBzzQ0K9Cvyotx33sr2evIlvCY0eosjwfRH7Kw9YPcD7iMdmOvdLxsyOTeA1uwc4Imrq8qzzBLbXZTzxHw5/Bj4GrI2KLxu4Av8TasSx+P6jERuBd5hZcdR+e+z3JzKEKV93lkn5uttj7GE/dcFz8uib6/Gt/j/tZnmt6zgRXbl7gKm4Tj+fNbOxUY/HAf8M7rfgT+A4GD/sU7mZ5ZvZ/GD5SDMbge/79xv8TzrPEPzMGMgi5r91M7Nethe7b4vaxt+BBuADwbaGBfuLtIB2t7/ujAnWGRc8959x7CM6nt72GX0csbpqxTB8a1QxcJ+ZnWH+xJne3oP2fQFjnHMro6bWqF1kR93vaXvRcXcZj5nlmdnCoMtHT/6C/zn6TgDnXD3wAPC06zxKyOn4LiMRkdf7NvNjZV8ORMY6jX0PYt/ntB65QKQbytdkTL6O93Xscj/4XJqH/wcM/K8UmNnYOHK34VuRPwxcYGb/ETz3S2Z2Nb7V/l9R6yt3D7SB6MitaWAm/EkSDtgJ/AH4Lf5npYpg+beAt4H34X/meR3fUnAWvsXjt/hEdQ/+5IqWYHtN+C/bu4FW/EkzE4Arg+Uf72l7Mfu+AngS349tRrBsNr7F8wZ838OZwfxu99fN8bvg+VfguykcFrWsu31Mjo6nt31Gv4Yx+34f/sS+nwWvReR5HwuWzwWWBu/NDcHrHO9r9r4ujvX8YH+3AeOCeT29p5221U084/EXO/hTHJ+1G2MeXwicHDPv+i6edxKwHNgOXB3MOxP/WbsJGEswHBnwKfwf318Fx/qesL9jmjQla0L5OpPy9UTgKXzr+xG9HGNP7/VZ+F8Mn8OfXL4auCw49i5zN/7XkVZ8i3Ul8G98Yf8J4JP48b5vwHdJiTxHuXuAJwteQMkQZjYS+JRz7jvB4yx8a8onnXNXhxpcL8zMAVOc78MofWRmX3HO9WmYJxEZPMrXEk25O32oW0jmuRzfRaAyat4pRJ2oloqifobST1B9ZGbZZnYxcFfYsYhIXJSvRbk7Dam4zjz/i++btcLM3gIeBt52zj0bbljdC07wuCJ4eIn5IeEkcSXA351za3tdU0RSgfK1gHJ32lG3EBERERGRJFHLtYj8f7t1LAAAAAAwyN96EjuLIgBgItcAADCRawAAmMg1AABM5BoAACZyDQAAE7kGAICJXAMAwESuAQBgEgyyv0E6S8fEAAAAAElFTkSuQmCC\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/viking-2-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/viking-2-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/viking-2-mars.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "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 }