{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 20 - Viking-1 - 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-1 Lander.\n", "https://en.wikipedia.org/wiki/Viking_1" ] }, { "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-1', 980, 63.0, 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.42,0.0,-16.99,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": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZwcZbX/8c/JNtnXyR6SISEJgUBABoIisiNe8Bq4KAr8BJWM4jVuF4kCboheRUC5IF5GRERFvewS2VEiiywJGEhCEghZIAlhspHJnsmc3x9P9UxPZ9ZMT1d19/f9etWrq6urq071JGfOPP089Zi7IyIiIiIi7dMp7gBERERERAqBCmsRERERkSxQYS0iIiIikgUqrEVEREREskCFtYiIiIhIFqiwFhERERHJgi5xB5AtpaWlXlZWFncYOVNVVcXgwYPjDiMR9FkE+hzq5dtnMXfu3HXunj8BZ8G+5ux8+9nuq2K5Tiieay2W64TCv9bmcnbBFNZlZWXMmTMn7jBypry8vKiutzn6LAJ9DvXy7bMwsxVxx5Br+5qz8+1nu6+K5TqheK61WK4TCv9am8vZ6goiIiIiIpIFKqxFRERERLJAhXWeqqioiDuExNBnEehzqKfPonAVy8+2WK4Tiudai+U6obiuNZO5e25PaHYKcBVwjrsvj7adBVwL9AV+D/yXu9dEr5UCPwQ2AV2Bme6+O/O45eXlXsj9eUSkcJnZXHcvjzuOxihni4g01FzOzmmLtZkNBXoDR6VtGw1MA84GvgZcBHw17W13Are6+0xgCXBlzgIWESliytkiIm2T08La3dcC92dsHgN81t3nuvvtwC+AEwDM7Ghgors/H+37ADDDzHrlKmYRkea4Ozt37iTX3/7lQpJytjvU1IRFRCSpct7H2t1rM54/lfoKMbIKWBmtnwgsS9t3VbSayK9MRaT49OnTh+7du7Nt27a4Q+kQScnZzzwDXbvCcce190giIh0niYMXjwT+N1ofCazPeL0aGJH5pqqqKsrLy+uWysrKDg5TRAS6desGwI4dO9r0vsrKyrp8BZRmP7KcyUnONguPBfjFgIjkgdbm7ERNEGNmBwBV7j4v2uTA9ozdSoC9BsIMHjy4oG9GLiLJ1L17dwB27tzZpvdVVFTUjZw3s3VZDywHcpmzVViLSJxam7MT02JtZp2BzwPfTNu8ChiQto8BvYDVuY1ORKRxJSUlQNtbrPNdrnO2CmsRyQeJKayBS4Br3D292ecxYGLa8zJgF6CmaRFJhFSLdbEV1uQ4Z6uwFpF8kPPCOmrBALC0bZcBLwM9zGysmX3WzA5w9znAGjObHO16GnCDu+/KbdQiIo0r9MI6KTlbhbWI5IOc9rE2sz7A+dHTC83sRsJXiT/I2HWRu98arZ8DfMfMlhPup3p5LmIVEWmNfe1jnQ+SlLNVWItIPshpYe3u1cAvoyXlqmhp6j0rgM91cGgiIvukkFusk5SzVViLSD5IUh9rEZG8U6yDF3NNhbWI5AMV1iIi7VDILdZJosJaRPKBCmsRkXZQYZ0bKqxFJB+osBYRaYdCHryYJCqsRSQfqLAWEWkH9bHOLRXWIpJkKqxFRNpBXUFyQy3WIpIPVFiLiLSDCuvcUGEtIvlAhbWISDuoj3VuqLAWkXygwlpEpB3UYp0bKqxFJB+osBYRaQcNXswNFdYikg9UWIuItINarHNDhbWI5AMV1iIi7aDCOjdUWItIPlBhLSLSDhq8mBsqrEUkH6iwFhFpB7VY54YKaxHJByqsRUTaQYMXc0OFtYjkAxXWIiLtoBbr3FBhLSL5QIW1iEg7qI91bqiwFpF8oMJaRKQd1GKdGyqsRSQfqLAWEWkH9bHODRXWIpIPVFiLiLSDWqxzQ4W1iOQDFdYiIu2gwjo3VFiLSD7okusTmtkpwFXAOe6+PNrWA7gaqAYGADPdfXNLr4mIxK3QBy8mJWersBaRfJDTFmszGwr0Bo7KeOkm4Gl3vwy4H7i5la+JiMSqkFusk5SzVViLSD7IaWHt7msJibaOmY0AzgUeiDY9Dvy7mZU191pOAhYRaUEhD15MUs5WYS0i+SDnXUHcvdZSGTI4Hljn7tui12vMbBlwHLC7mdeW5zJuEZHGdOnShU6dOrFnzx5qamro0iXnabVDJSVnq7AWkXyQhN8AI4H1GduqgRFATTOviYjEzszo3r0727ZtY+fOnQVXWDcilpytwlpE8kES7griwPaMbSWElo/mXmugqqqK8vLyuqWysrJDghURybQv/awrKyvr8hVQ2jGRdYhYcrYKaxGJU2tzdhKaVlYRRo6n6wOsJiTppl5rYPDgwcyZM6dDAhQRac6+9LOuqKigoqICADNb1yGBdYxYcrYKaxGJU2tzdhJarP8OjIhu0YSZdQVGA7NbeE1EJBEK+c4gjYglZ6uwFpF8kPPC2upHwRiAu78DzAJOirafDNzt7quaey2HIYuINKuQC+uk5GwV1iKSD3LaFcTM+gDnR08vNLMb3b0KuBi42symAMOAirS3NfeaiEjsCnWSmCTlbBXWIpIPclpYu3s18MtoSd++EZjexHuafE1EJAkKtcU6STlbhbWI5IMk9LEWEclrhTxJTFJ07hwea2vjjUNEpDkqrEVE2qlQW6yTpGvX8Lh7rxv3iYgkhwprEZF2KtQ+1kmiwlpE8oEKaxGRdlKLdcdTYS0i+UCFtYhIO6mPdcdL72OtftYiklQqrEVE2kkt1h3PrL7VuqYm3lhERJqiwlpEpJ1UWOdGl+gGseoOIiJJpcJaRKSdNHgxN9TPWkSSToW1iEg7qcU6N1RYi0jSqbAWEWknDV7MDfWxFpGkU2EtItJOarHODbVYi0jSqbAWEWkn9bHODRXWIpJ0KqxFRNpJLda5ocJaRJJOhbWISDupj3VudOsWHnftijcOEZGmqLAWEWkntVjnRs+e4XHr1njjEBFpigprEZF2UmGdG717h0cV1iKSVFkprM3s37JxHBGRfJRvgxfzNWf36hUet2yJNw4RkaZ0acvOZnYacAUwjPqivBMwFOiR3dBERPJDUlusCy1nq8VaRJKuTYU18Dvgv4GXgdpomwHnZjMoEZF8kuDBiwWVs9ViLSJJ19bC+k13vy5zo5ktylI8IiJ5J6kt1hRYzlaLtYgkXVv7WF9sZmc3sv3fsxGMiEg+SnAf64LK2WqxFpGka2thfQxwi5ntSVtqgV92QGwiInkhwS3WBZWz+/QJj9XV8cYhItKUtnYFuRKoAOYAe6JtBvy/9gZiZpOALwFvAOOBSnf/l5n1AK4GqoEBwEx339ze84mIZEuC+1gXVM4eNCg8rl+fjaOJiGRfWwvrh4D73b3B951m9qssxPI74GPuvsrMRgOPAJOAm4CH3f3P0Qj3m4FPZeF8IiJZkeAW64LK2aWl4XHdumwcTUQk+9paWP8J+JaZ/S1tmwFnAV9pZyyTgOiLPrYD/cxsBGH0+n9G2x8H7jazMndf3s7ziYhkRYIL64LK2SqsRSTp2lpYf5vwld+FgEfbjHCP1PYm6bsIfQE/ApwPzACOB9a5+zYAd68xs2XAccDydp5PRCQrunTpQufOndmzZw81NTV06dLW1NphCipnq7AWkaRr6+DFnwL7uXuZu+8fLWWE1o/2+iKwA3gR2OLudwMjgczedNXAiMw3V1VVUV5eXrdUVlZmISQRkdZpa6t1ZWVlXb4CSjsorILK2SqsRSQurc3ZbW1WecvdGxuPPayNx2lMD+CPhMT886iVwwlfMaYrAXZnvnnw4MHMmTMnC2GIiLRdSUkJW7duZceOHfRO3XC5GRUVFVRUVABgZh1VKhZUzu7fHzp1go0boaYGkvPFgIgUutbm7La2WH8xc4OZfRz4eRuP05jfA3e7+5XAjwkJew1hVHm6PsDqLJxPRCRrUt0/ampqYo6kgYLK2Z07w8CBYX3DhmwcUUQku9paWA80s28DmFlPM7sFuIXwVeA+M7NSYIq7b4o2XQVsBlYAI6LbN2FmXYHRwOz2nE9EJNs6d+4MwJ49e1rYM6cKLmcPGRIeV6t5RUQSqE2FtbufDjxmZjcA/wImAFOAxmb2aosNwA4zG5m2bT3wCjALOCnadjKhhWRVO88nIpJVSSysCzFnl5WFx+XLs3VEEZHsabaHmpl9qImXVgCvA9cAY4BpwNf2NQh3rzWzacB3zGwuMBT4hrtvNrOLgavNbAqhX2DFvp5HRKSjJKGwLoacvf/+4XHZsmweVUQkO1oa+nE14V6lm4DaRl6/lXDrpuG0I0kDuPs84PONbN8ITG/PsUVEOloSCmuKIGerxVpEkqylwvpKYI67v9vcTmbW3q8VRUTyWkIK64LP2akW6zffjDcOEZHGNFtYu/uDrTmIu9+VnXBERPJTEgrrYsjZ48eHx8WL441DRKQxbb0riIiINCIJhXUxmDgx3HZv6VLYnnnHbBGRmKmwFhHJAhXWuVFSElqta2th0aK4oxERaUiFtYhIFqiwzp2DDw6PCxbEG4eISKZ2F9ZmdriZ/d7M/svMytofkohI/smXwroQcrYKaxFJqmy0WH8BONndrwU+kJpxS0SkmORLYU0B5OzJk8OjCmsRSZpsFNZfAo6I1v8EdM3CMUVE8oqZAeDuMUfSorzP2YceGh7nzoXkf9wiUkzaXFib2UfM7Lxo/ShgfGq6WnevdffNWY5RRCTxUi3VqZbrpCjEnD1+PAwcCKtXw8qVcUcjIlKvTYW1mf0QuBM4E8DdXwA+bWbndkBsIiJ5I4mFdaHm7E6d4P3vD+vPPhtvLCIi6draYv0BYBTwVNq2h4AfZy0iEZE8lMTCmgLO2R/4QHh85pl44xARSdfWwvoxd98EpPdqy9upcUVEsiWhhXXB5uxjjgmParEWkSRpdkrzRrxpZpcBI83sk8DHCV8xzsh6ZCIieSShhXXB5uwjjwwzMM6bB9XV0KdP3BGJiLSxxdrd/wT8DegMfBrYBJzq7r/ogNhERPJGEgvrQs7ZPXvC4YeHGRhfeCHuaEREgra2WOPuzwHPpZ6bWWczO8ndn8hqZCIieSSJhTUUds4+9liYMwceewxOOinuaEREWiiszewvQO8WjlEKrAPyPkmLiOyrJBTWxZazzzgDfvYzuP9++HHeD8cUkULQUov1K9E+i4FaYDrwMjAnbZ/9gYEdEp2ISJ6oqakBoEuXNn8RmE1FlbOPPRb694dFi2DJEpgwIe6IRKTYtfQb4KeARaPKMbNx7v6dzJ3M7P6OCE5EJF9s27YNgJ49e8YZRlHl7K5d4d/+De64A/76VxXWIhK/Zgcvuvt7qQQdGWpmDb7nNLPjgaM6IDYRkbyRhMK6GHP2hz8cHp/I+44tIlII2vqd5V3APDN7GqgGJgEfBvZqERERKRa1tbXs2LEDgO7du8ccTQMFn7NTgxaffBK2b4cePWINR0SKXFtvt/cYcCph4MtEwq2bznb3/+6A2ERE8sL27duB0FrdqVNb593qOMWQs0eOhCOOgK1b4fHH445GRIrdvtxubzVwReq5mXUysyPcfW5WIxMRyRNJ6AbSlGLI2WedBXPnwl13wUc/Gnc0IlLM2lRYm9mtjWweDOwBpmUjIDM7Gvgg8AbwNLAVuJrwNeYAYKa7b87GuUREsiGphXWx5Oyzz4bLL4d77oGbboJevTrybCIiTWvrd5bjgZXAirSlBMhKy4eZfQY4w92vcff73H0dcBPwtLtfBtwP3JyNc4mIZEtSC2uKJGdPmABHHw1btsC993b02UREmtbWriAXufvi9A1mdiBwYXsDMbPjgPOAU9K2jQDOBf4z2vQ4cLeZlbn78vaeU0QkG6qrqwHolbym0qLJ2RdcAM89B7fdBuef35FnEhFpWlsHLy5uZPNO4KIsxHId8Bpwg5k9bGbHAMcD69x9W3T+GmAZcFwWzicikhXr168HYNCgQTFH0lAx5exzzoGSEvjb38KEMSIicWhTYW1my8zszbRlBSGxtuvLNzObCBwO/Mrdv0SYavdhYD9gfcbu1cCIzGNUVVVRXl5et1RWVrYnJBGRVksV1gMHtn5Cw8rKyrp8RZhmPOuKKWcPGACf/jS4w4wZUFvbrsOJiDTQ2pzd1q4gjwF3AB49rwXWuvuSfYqy3sHABnd/JXp+I/Dd6DzbM/YtAXZnHmDw4MHMmTMnc7OISIfbsGED0LYW64qKCioqKgAws3UdEliR5eyrroK77w633fuf/4GvfjWrhxeRItbanN3WwYu/d/cn3X12tDwVjm8HtCNWCAV+XZHv7tuB14HOhFHl6foAq9t5PhGRrElqVxCKLGcPGQK//nVYnzkT/vWvXJxVRKReWwvr0xrZtha4s51xvAL0M7OhadtqCMl4hJn1ADCzrsBoYHY7zycikjUJLqyLLmdPmwZf+ALs2gWf+AREXyaIiOREqwprM/uSmW0CZprZnvSF0J9ubXuCcPdFwEPA2dH5+gPdgD8As4Bo0lpOBu5291XtOZ+ISDa98847QOjekATFnrOvvRamTIHXX4f/+I9QZIuI5EKr+li7+41m9g/gAuB/Ml7e5e5rshDLp4HrzawboYXjU+5eY2YXA1eb2RRgGFCRhXOJiGTNihUrABgzZkzMkQTFnrN79oRZs+Coo+DJJ+HCC+F3v4POnXMdiYgUm1YPXnT3V8zsm+6+1yAUMxvt7ivbE0g0scB5jWzfCExvz7FFRDpSqrAuKyuLN5A0xZ6zR40KxfXxx8Mf/wjdusGtt0KntnaAFBFpg2YLazO7BHjQ3RdGmy41s8y/+Q04gXD/UhGRorJ161aqqqro1q0bw4YNizUW5eyG3vc+ePBBOO00+O1vYceOMLgxefP4iEihaKnFen+gf9rz9wG9gHeov32TEb7uExEpOqnW6v32249O8TeHKmdn+OAH4YEH4KMfhT//GebPh9tvD0W3iEi2NVtYu/t/Zmz6ITDf3RsMBTGzg7IdmIhIPpg/fz4AkyZNijkS5eymnHACvPACnHUWLFgQ+l5/61twxRVhtkYRkWxpa/PKIY0k6InAuOyFJCKSP+bNmwfAlClTYo6kUcrZkYMOgrlz4StfCbMyXnUVHH44PPZY3JGJSCFp1eBFMzuW8PXhKWa2LOPlgYRR5w9kOTYRkcRLFdaHHnpozJHUU85uXK9e8POfh1vwfe5z8NprcOqpoZvINdfAhAlxRygi+a61dwXZBNxO6L93TMZru4BbshmUiEg+cHdeeOEFAA4//PCYo2lAObsZxx4Lr74K118PP/hB6IP94INw/vmhe8gB7Z2XUkSKVqu6grj7q4QR5Je4+/4Zy0R3v7JDoxQRSaAFCxZQVVXF8OHDOSBB1ZhydstKSuDSS8MkMhddFLb99rdw4IFwwQUQfREhItImre5j7e7vuXujrRxm9sHshSQikh/+/ve/A3DiiSdiZjFH05BydusMGwa/+hUsWRK6h0C4a8hhh4VBj/fdB3v2xBujiOSPlu5j/RegdwvHMOBgYEi2ghIRyQcPPfQQEArrJFDO3ndjx8Itt8Bll4UuIrfeGmZtfPJJGDMmFN2f+UyYeEZEpCnm7k2/aPZTQl+9VdTfA3Wv3YAz3P3s7IfXeuXl5T5nzpw4QxCRIrJ+/XqGDRuGu7NmzRoGDx68z8cys7nuXt7emJSzs+e99+A3v4EbboA33wzbOnUKk81cdBGccQZ07RpvjCISj+ZydkuDF38CVLv7zhZOUHSjy0WkuN17773U1NRwyimntKuozjLl7Czp1w+++lX48pfhiSdCa/a994ZBjg8+CEOGwNlnwznnhElo4p8bSESSoNlU4O7rWkrQkTOzFI+ISOK5O5WVlQB86lOfijmaesrZ2depE5xySpi1cdUquO66cE/sd9+Fm26C446D0aPha1+DZ59Vf2yRYteuv7HNbLSZ/QS4PkvxiIgk3nPPPceLL77IwIEDOeecc+IOp9WUs9tn8OBQQM+fHyabufTS0P961apwf+xjjoHhw+HCC+Guu2Dz5rgjFpFca3NhbWadzexMM3sQWAqcR+jTJyJSFK655hoApk+fTs+ePWOOpnnK2dlnBu97H/zkJ7BsGTz3XCi4998fqqrCbfs+/nEoLQ2t3T//eSjGmxnSJCIForUTxGBm+wMXAZ8BBgErgI8CjwDv75DoREQS5oUXXuCee+6he/fuzJgxI+5wmqScnRtmMHVqWK69NszmOGtWmHTm2Wfh8cfDAjB0KJx4Ipx0UljKymINXUQ6QIst1mb2cTN7FHgd+A/gOmAU8Fd3f9iDZzs4ThGR2Lk7l156KQBf+cpXGDlyZMwR7U05Oz5mof/1pZfCU0+Ffti/+12Y0XH4cFi7Fv74x3BXkf33D7f4u+ACqKyEhQuhtjbuKxCR9mpNi/V2whS4zwIXufsSADPTl1oiUlR+85vfMHv2bAYNGsTMmTPjDqcpytkJMWhQKKrPPz90A1m0KNxh5Iknwv2xly0Ly+23h/0HDAj9tI85JtxppLwcuneP9RJEpI1aLKzdfRYwy8xGA9PNrAy4H6i7g6eZmTd3Q2wRkTz39ttv8/Wvfx2A66+/ngEDBsQcUeOUs5PJDCZNCsuXvhTuHjJvHjzzTFiefjoMgpw1KywA3brBEUeEIjtVcJeWxnsdItK8ZieIafQNZl2AaYQBMPOBXwOnuntl9sNrvaRPNiAi+WvXrl0cf/zx/POf/+SMM87gL3/5S1anMM/WBDFNHFs5Ow+4w8qVocBOFdqNDXicOBGOOgqOPDIsU6ZAjx7xxCxSrJrL2W0urDMOfAAwA5ju7rEOjVeSFpGOMmPGDG688UZGjRrFSy+9lPUJYTqysM44j3J2Htm0Cf75z/pC+4UXYPv2hvt06QKTJ9cX2kceCQcfrFkhRTpSe2ZebJa7vwF8xczmtuc4IiJJdfXVV3PjjTfSrVs37r777iTNsthmytn5pX9/+MhHwgKwa1foPjJnDrz4YlgWLoR//Sssv/pV2K9bt1BcT5nScBk4ML5rESkW7SqsU9z99mwcx8y6AS8CX3H3J82sB3A1UA0MAGa6u265LyI5ccsttzBz5kzMjNtuu42jjjoq7pCyIls5G5S3c6lbt/pW6YsvDtu2bIGXX25YbL/xRtj28ssN3z9qFBx2WMNie+zY0OotItmRtP9O3wDK0p7fBDzs7n82s9OAm4HkzB8sIgXrhhtu4Mtf/jIAN954Y6KmLk8Y5e0Y9e4Nxx4blpTNm+HVV0Pr9rx5oTX71Vfh7bfDkhocCaFYnzgx3CYwtUyaBOPHh9dEpG0SU1ib2QeANcDG6PkI4FzgP6NdHgfuNrMyd18eS5AiUvDcne9973tceeWVQJhl8Ytf/GLMUSWT8nYy9e1bfxeRlD17Qkt2qtCeNw9eeSUU2q++GpZ0nTuH4jpVaE+aBBMmhKVfv9xej0g+SURhbWa9gLPd/etm9p1o8/HAOnffBuDuNWa2DDgOWB5LoCJS0Kqrq7nwwgu555576NSpE5WVlXzuc5+LO6xEUt7OL507h5bpiRPhE5+o3755c7i/9sKFYdbIhQvDsmxZ2L5o0d7HGjIkHCdVaKfWx46FkpLcXZNIEiWisAZmAj/O2DYSWJ+xrRoYkZOIRKSovPrqq3zyk59k4cKF9O3blzvuuIPTTz897rCSTHm7APTtG27flzl8YPt2WLy4vtBevBiWLIHXXw8zSr77bphdMl2nTmGa9saK7pEjw+sihS72wtrMPgzMcfd3M15ywgxi6UqA3Y0dp6qqivLy+jufVFRUUFFRkc1QRaQA1dbW8rOf/YzLLruMXbt2cdBBB3HvvfcyYcKEDj1vZWUllZV1t5LOq2k/spG3lbOTrUePMNDxsMMabq+tDd1HliypX1JF9/Ll8OabYXnoob2PN358w6I7VXgndK4lkQZam7PbdR/rbDCzx4DD0zYNALYA1wHnufuEtH1fB77r7ndkHkf3RBWRtlq0aBFf+MIXmD17NhCKu2uvvZbevXvnNI5c3cc6W7KRt5WzC8/OnbB0aeNF97uZf4KlKS1t2LqdWsaN0+Q3kkwddh/rLDmX0KKR8k/g68DzwDfMrIe7bzezrsBoYHYMMYpIAdm6dStXXXUV1157Lbt372bo0KHccsstnHHGGXGHli+Ut2UvJSX1dxbJtGlT6EaSKrTTl3XrwvLss3u/b/jw0Hd7//33fhwxIvQdF0mS2Atrd69Kf25me4Aqd19pZrOAk4BZwMnA3e6+KoYwRaQA7Nmzhz/84Q9cccUVvPXWW5gZ06dP50c/+hGlpXnVGyNWytvSVv3719+DO507rF7dsHU7tb58OaxZE5Znntn7mN26wZgx9YV2ZvGtLiYSh9gL6xZcDFxtZlOAYYA64IlIm7k79913H1dccQULFy4E4IgjjuAXv/gFU6dOjTm6gqO8La1mFgY2jhwJJ5zQ8LWaGli1KvTZXras/jG1vnZtaAV//fXGj92vXyiyx4wJy+jRDR+HDAnnF8mmxBXW7l6Wtr4RmB5fNCKSz/bs2cM999zDT37yE+bODbN4l5WV8f3vf5/zzjuPzvoeOSuUt6UjdOlSXxRnFt0AW7eGVu30ojv98b336qd7b0xJyd7Fdvr6qFGaJEfaLnGFtYhIe23bto3bbruN6667jqVLlwIwdOhQvv3tbzN9+nS66belSN7r1QsOPjgsmdxDv+1ly2DlSlixIiyp9ZUrYcOG5lu8zUIf78yie/ToUHSPHBkGXuo2gpJOhbWIFIzFixdTWVnJbbfdxoYNGwAYO3Ysl1xyCRdccAE9e/aMOUIRyQUzGDw4LJn36E6prg4FdnqxnV6Ar1oV+n+vXg3PPdf4Mbp2re/KMnJkfcGdvj5ihFq+i4kKaxHJazt27OC+++7j5ptv5sknn6zbfuSRR3LppZdy5plnqsuHiOylT5+mW7wBdu8ORXVma/dbb4Wi++23YePG0B1l+fLmzzVkSMvFd79+6vNdCFRYi0je2bNnD7Nnz+YPf/gDd999N++99x4APXv25Nxzz6WiooLy8nJMv6VEZB917VrfBaQp27aFIju1vP323utr1tTPVvnyy00fq0eP0PVkxIi9l/TtffqoAE8yFdYikhdqa2t5/vnnuUgxj2oAACAASURBVOuuu/jTn/7E6tWr6147/PDDmT59Oueddx59+/aNMUoRKSY9e4YZJcePb3qfmppwB5OmCu/Vq8Pj1q31M1c2p1evpovu9G05nudKIiqsRSSxduzYwd/+9jfuu+8+/vKXv7B27dq618aOHcu5557Lueeey6RJk2KMUkSkaV261Hf9aKq/N4Q+36k+3atXh5bu9OepZevW5gddpvTp03jrd3oxPnx4+ONAskeFtYgkyltvvcWjjz7KQw89xMMPP8zWrVvrXhszZgzTpk3jk5/8JFOnTlVXDxEpGH36hGndJ05seh932Ly58YI7sxCvrg4T7Sxe3Px5+/dvvvU7tb2kpPnjSKDCWkRitWXLFmbPns2jjz7Ko48+yqJFixq8fthhhzFt2jQ+9rGPMWXKFBXTIlK0zMIgx379oLkv6tzDwMqWWr/XrAnTzW/aBNHcWU0aOLB+oGWqBT61nnocMkS3H1RhLSI5VV1dzbPPPstTTz3FP/7xD5577jl2795d93rv3r058cQTOfXUUznjjDMY09zIIRER2YtZKIQHDoTJk5ver7Y23M+7pdbvNWvCfhs2wKuvNn28Ll1g2LCWC/BCHgqjwlpEOtS6det4+umn6wrpl19+mT179tS93qlTJ6ZOncqpp57KqaeeytSpU+natWuMEYuIFIdOncIkN6WlcOihTe9XWxvuapIqtFMDMDPX168PgzLffrv58/buXV9kjxtXPwB0woTwvHv37F5nLqmwFpGs2blzJ/PmzeP555+vW954440G+3Tu3JmpU6dy7LHH8qEPfYgPfvCDDBgwIKaIRUSkJZ06hZboYcPgfe9rer8dO0LrdmbhnV6Ar1oFW7bAkiVh+fvfGx7DDPbbLxTaEyfCQQeFe40fdFDoapJ0KqxFZJ/U1tby5ptv8sILL9QV0S+//DK7du1qsF+PHj2YOnUqH/rQhzj22GM5+uij6a37QImIFJzu3WH//cPSFPfQp3v16jDZzhtv1N/l5PXX66ehX7kSnnii4XtLS+sn9UlfSks79rraQoW1iLRo9+7dLFy4kJdffrlumTdvHps3b95r30mTJjF16tS6ZfLkyeraISIiQGiRHjAgLI3Nerl7d5jJcsmScEeTBQvCsnAhrFsHs2eHJd2QIaFF+5BDwi0Np06FAw6IZyIdFdYi0sCmTZtYsGAB//rXv+qK6Pnz5+/VEg0wbNgwysvL64roI488kv79+8cQtYiIFIKuXev7XJ9+ev1299B3O1VopxfcqZktn3yyfv+BA0OR/f73w0knhWK7Sw6qXhXWIkVq69atLFy4kPnz57NgwQLmz5/P/PnzWbVqVaP7jxs3jsMPP7zBMmzYsBxHLSIixSjV93q//eC00+q3u4cuJQsWhCnjn38+LGvXwsMPh+W73w13IjnpJPjIR+Dss0OLeUdQYS1S4DZt2sTixYtZsmQJr732Wl0BvWzZskb37969OwcddBCHHHJIXQE9ZcoU+vXrl+PIRUREmmcGo0eH5SMfCdvcQx/t55+Hp56CRx8NXUvuvTcsM2bAOefAVVeFQj2bVFiLFIDdu3ezbNkyFi9evNfy7rvvNvqeLl26cOCBBzJ58mQOPvhgJk+ezOTJk9l///3p3Llzjq9AREQkO8xgzJiwfOITYdvy5aHAvvPOMCjy9tth1qzQyj16dPbOrcJaJE/U1NSwcuVK3nzzTZYuXcrSpUvriuelS5dSU1PT6Pt69OjB+PHjmThxYoNCevz48XTr1i3HVyEiIpJ7ZWVQURGWe++Fs84KE9489BB8/vPZO48Ka5EEqa6ublA4p6+vWLGiwcQqmUaPHs3EiRP3WkaNGkWnYp9jVkREilKqW8icOfCPf8Ajj4S7jQAMGtSwv3Y2qLAWyaFt27axcuVKVqxYUfe4bNmyugK6qqqq2fePGjWKsWPHMm7cOMaOHcuECROYOHEi48ePp2fPnjm6ChERkWRJ3TVk8eKwLFkCixbBSy+F2/Sl698fPv1pmDkzzACZTSqsRbLE3dm4cSMrVqyoW1LFc2ppqXDu3r07Y8eObVA8jxs3jnHjxlFWVkb3fJ7nVUREZB+4h9ka33knFM8rVoRW6NRjatmxo/H3DxoERxwRbr936qlw9NHhtn4dQYW1SCu4O5s2beLtt99m1apVrFq1qm797bffriugt2zZ0uxxunbtyujRoxkzZkyDJVVEDx8+XN02RESk4NXWhhkY168PLcrr1oVb5L3zTsNlzZrwuG1by8csLQ3ToB94YHicOBEOPTQMYszVZDGJKazN7CzgWqAv8Hvgv9y9xsxKgR8Cm4CuwEx33x1fpFJo9uzZw9q1a/cqljMft2/f3uKx+vTps1fRPGbMmLpietiwYSqcpSAoZ4sIhAK5ujoUyenL+vX1RXNjj+vXh/e2Vo8eMHx46LoxZkz9LfbS1/v06bjrbK1EFNZmNhqYBpwNHAz8EngLuAa4E/imuz9vZl8ArgS+FVeskj927NjBO++8U7esWbNmr/U1a9awZs2aZgcFpvTp04dRo0YxcuTIvR73228/xowZQ//+/bE45lAVySHlbJHC4B66T1RXh2Xz5saL5Mxl48b69ffeC8fZF/36hW4apaXhcdiwppc+feKZorytElFYA2OAz7p7DTDXzCYDJ5jZ08BEd38+2u8BYLGZXeXuW+MKVuKzZ88eNmzYwNq1axstlNPXN23a1OrjDhkypNGCOfU4cuRI+vbt24FXJpJXlLNFYlJTA1u3hiWzIM5cb+n55s3QinalFvXpEwYEpi8DBoSCOVU0pz+WloYpxzuqn3OcElFYu/tTGZtWAX2AE4FlafutiloDy4HZOQtQOkxtbS0bNmygqqqKd999t9HH9PX169dT28rvjrp06cKwYcMYPnw4w4YNa3J9+PDhlJSUdPCVihQO5WyR5u3aVV/8tmbZsqX1++7ald1Yu3UL03336VO/DBiwd6HcWOHcv394b5dEVJPJkNSP4kjgp0AFsD7jtWpgr5ujVFVVUV5eXve8oqKCioqKjoxRMrg7W7duZf369WzYsIH169fXLevWrWu0aF63bl2rC+WUAQMG1BXHzRXLAwYMUH9mSazKykoqKytTT0vjjCULlLMl8WpqYPv2sGzb1rrHtuybety6NZyro5hBr15h6dNn76I4/Xlr1jVPWOu0NmcnrrA2swOAKnefZ2YOZI4YKwH2GggzePBg5syZk4sQi8KuXbsaFMeZhXJTz3ftw5/S/fv3Z8iQIQwePHivx8xtgwYNomshfnckRSe9kDSzdS3snljK2bKv3EPra2YB25rn+1L47s7hENouXeqL32wv3bvnR1/jQtPanJ2owtrMOgOfB74ZbVoFTEh73YBewOrcR5efdu7cycaNG9mwYUNd8Ztab2xJFcot3TauKd27d2fQoEEMGjSIgQMH1q2XlpY2WiyXlpZqWm2RPKWcXXgaa9XtyOf7OuhtX3TqFO4s0bNn44/Nvdbaxx49QvGrX2vFK1GFNXAJcI2774yePwZ8Ie31MmAXUHTNHDt27GhVcZz52tat+zZeqHPnznWFcXqBnL7e2PMePXpk+cpFJMGUsztYbW24a0NLRWq2Ct2O7MLQmK5d6wvSzOK2qeepbW0tert1U0uvdLzEFNZmdhkh+fYws7HA8cA/gDVmNtnd5wOnATe4e5a77ueGu7N9+/ZWtRpnbmvNPZQbk1kgN7WkXh8wYACDBg2ib9++6p8sIk0qhpzdFu4Nb1OWfjuyjRvDLcma6srQ3POdO1s+dzalWnVbU+A2VfC25bkGvUmhScQ/aTO7AvhBxuZF7n6rmZ0DfMfMlgO9gctzHV9rbdmyhTfeeIPly5fvtbz77rts2LCBnfuYJbt27dpscdzUa3369NF9lUUkqwolZ7fGjh2wdGmYOnnNmvpZ4FLra9fWF9BtHIfdat27t71g3dciuGtXteqKtEciCmt3vwq4qonXVgCfy21ELauqquKZZ57hpZde4tVXX+WVV17hzTffbPF9JSUlTRbBzRXOvXr1UoEsIomQjzm7JbW1sGQJvPACvPgiLFoUnr/1Vuv7Affu3fA2ZOm3JevXr/H+uC0910A1kfySiMI6H+zcuZMnnniCBx54gNmzZ/Paa6/ttU+XLl0YP348Y8eOpaysrG4ZM2YMw4cPVx9kEZEE2bgRZs2CBx6Axx4Lrc6ZOneG/fcPy4gRYQa44cPrl6FDw0QX/furW4OIqLBulrvz9NNPc/PNN/OXv/yF6urqute6d+/O+9//fqZOncohhxzCIYccwsSJE3WHCxGRhFuwAK67Du64I3T1SBk1Co46Co48Eg45BMaPDwW17vApIq2lwroRtbW13HnnnfzoRz/ilVdeqdt+2GGHceaZZ3LKKadwxBFHqIgWEckj1dUwcybcfHN9f+gTToBp0+CMM2Ds2HjjE5H8p8I6w7PPPsuMGTN46aWXABgyZAjTp0/nM5/5DOPGjYs5OhER2RfLloXieeHC0L3j4ovhq1+FCRNafq+ISGupsI7s3LmTyy+/nOuuuw53Z8SIEXz3u9/lggsuoKSkJO7wRERkH61dCyefDG++CZMmwf/9H0yeHHdUIlKIVFgT7vBx1lln8fTTT9OpUye++c1v8u1vf1sDDUVE8pw7XHRRKKrLy+GJJ6Bv37ijEpFCVfSF9apVqzjuuONYunQpI0eO5K677uLoo4+OOywREcmCRx4Jd/7o1w/uu09FtYh0rKIurNeuXcuJJ57I0qVLOfzww5k1axYjRoyIOywREcmS664Lj9/6FowcGW8sIlL4iraw3rlzJ2eddRZLlizh0EMP5fHHH2fgwIFxhyUiIlmycSP87W9hsGJFRdzRiEgx6BR3AHG55JJLePbZZxk1ahSPPPKIimoRkQLz97/Dnj1wzDFh9kMRkY5WtIX1RRddxIEHHsi9997LsGHD4g5HRESybOBAOPPMcJ9qEZFcKNquIFOmTGH+/Pl07tw57lBERKQDHH98WEREcqVoW6wBFdUiIiIikjVFXViLiIiIiGSLCmsRERERkSxQYS0iIiIikgUqrPNUZWVl3CEkhj6LQJ9DPX0WhatYfrbFcp1QPNdaLNcJxXWtmVRY56li/kebSZ9FoM+hnj6LwlUsP9tiuU4onmstluuE4rrWTCqsRURERESywNw97hiywsyqgBVxx5FDpcC6uINICH0WgT6Hevn2WYxx98FxB5FL7cjZ+faz3VfFcp1QPNdaLNcJhX+tTebsgimsRURERETipK4gIiIiIiJZoMJaRERERCQLVFjnITM7xcyeN7OyuGOJi5mdZWbLzGy9mV1vZl3ijikuZnacmS02s41m9vO440kCM+tmZvPM7Pi4Y5HsMLMeZnaDmf3IzH5pZn3jjqkjFGNuK5b/r2Z2tJldYmbTzKw07ng6gplNMrNfmNnXzOwmMzss7phyreD/wxYaMxsK9AaOijuWuJjZaGAacDZwMPBL4C3gmjjjioOZ9QOOAY4GPgDcZ2az3P3xeCOL3TeAsriDkKy6CXjY3f9sZqcBNwOfijmmrCri3Fbw/1/N7DPAOHe/Iu5YOtjvgI+5+6ro3/MjwKSYY8opFdZ5xt3Xmtn9cccRszHAZ929BphrZpOBEyj8Xz6N2QX8t4dRyH81s1eAPTHHFCsz+wCwBtgYdyySHWY2AjgX+M9o0+PA3WZW5u7LYwss+4outxXD/1czOw44Dzgl7lhyYBLQJ1rfDvSLMZZYqCtIHnL32rhjiJO7PxX94klZBayMK544ufv2qKjGzHoBrwFPxhpUjKLP4Gx3vzXuWCSrjgfWufs2gOj//zLguDiDyrZiy21F9P/1OkJuvsHMHjazY+IOqAPdBdxiZn2A84EZMceTcyqspRAcCfxv3EHEycyOBR4CugI9Yw4nTjOBH8cdhGTdSGB9xrZqYEQMseRSoee2gv//amYTgcOBX7n7l4AngIejbp2F6IvADuBFYIu73x1zPDmnwlrympkdAFS5+7y4Y4nZG8BtwEnAT+MNJR5m9mFgjru/G3csknVO+Fo5XQmwO4ZYcqLQc1sR/X89GNjg7q9Ez28EDPhYfCF1qB7AH4E7gJ+b2ckxx5Nz6mMtecvMOgOfB74Zdyxxc/c1wK1mVksYCFSMLgEON7PU8wHA/Wb2Q3e/Or6wJAtWEX6e6foAq2OIpcMVSW4rlv+vXUirtdx9u5m9DgyKL6QO9Xvgk+6+ycwc+GM0FmJr3IHligpryWeXANe4+864A0mQOYQipBidS2jFTPkn8HXCqHTJb38HfmVmPaLCpCswGpgdc1wdpRhyW7H8f30F6GdmQ919bbStBlgYY0wdIrqF4BR33xRtugq4kDCgcU5cceWauoLkIav/E9+a3bGAmdllwMtADzMba2afjb46LSpm1tPMDk/b9G/A9XHFEyd3r3L3t1ML4e4oVe6+Oe7YpH3c/R1gFqGrE8DJwN3uXnB/RBZLbiuW/6/uvogw/uVsADPrD3QD/hpnXB1kA7DDzEambVsPLIkpnlioxTrPpI20BbjQzG5096o4Y8o1M7sC+EHG5kVFMLK8MZMJt9lbBDwHvODuhZiwRS4GrjazKcAwoCLmeLJOua1gfRq43sy6Eb5p+VTG3V8KgrvXmtk04DtmNhcYCnyj0P5YaolFd+oSEREREZF2UFcQEREREZEsUGEtIiIiIpIFKqxFRERERLJAhbWIiIiISBaosBYRERERyQIV1iIiIiIiWaDCWqQJZnaSmT1rZse3sN9kM/s/M/t2Fs/902wdS0QkDmZ2qJndmc3c2MR5ZpnZJzvyHB3FzE42s3+29HtG8ocmiJHEM7NJwA8JszodCRwK/D93/30Hn/otYEor9lsJDAA6Z+Ok0XTNu7NxLBERADM7Evhv4EPAo8AywkQ7XYFfuntHTCX+BtCPLOXGZtwEzO+IA0f5+CzgfwkzCP4zeqkE+A/gKHdf3o5TtPb3jOQJFdaSD+4DPuHu88ysE/CbXJzU3ZeY2fpW7LfZzNZk8dTHAf/I4vFEpMi5+4tm9ifgQHc/I7XdzE4H7jWzS9z9piyfc1uWcyMAZvZR4NVUQevuD2b7HCnuvhv4s5nNAJ509yvS4rh/X45pZkOAE9z9z+6+2MzWZSlcSQB1BZFEM7OBwARgK4QpU4HvALmaMrQ2y/u1xgnA37N4PBERgL2m0Xb3vwLfA35mZsM64JxZzdVmNorQQp1rjU1B/jiwqi0HMbMS4HagR9rmbP7+kJipsJZEc/cNwGLgYTM7Otq2AnjBzPqb2Q/M7AUz+3cze8fMXjOz96Xeb2afMrPLzexJM/tR2vZyM/uemf3ZzO4xs37R9u5m9gsz+46Z/ZrwNWabmNmlUUyfMLMjzey+6Hg/iGL8s5kNjra/Z2bfzThEN3ffGR3rv8zs62b2jJn9uM0foIhIy35N6BLycQAzG2dmV5rZLWb2RFTMYmYlUS673MyeMrNzUgcwsy+a2bfN7Dkz+1JTJ2osJ5vZ2VHOvNDMXjezH5vZCDP7TZRPnzKzE6NDnAuMAmaa2bRoLMzjZnZB2jmmm9kPzex2M/utmfUys35m9l0zm2NmJ5rZYjN708zK9uUDiz6T46MWbczsLDO73syuNbOHzax3dA3Xm9kD0e+bNdFnfBDwcTP7etohB0S/i6rNbHoT5+xkZj8xs4roHG5m96Z+PpIQ7q5FS6IX4GBgBeGv+luBIdH2zoQktYHQ160/8BiwgPBH4/uBK6J9BxBavU8H+gJ/Sjv+M8BPo/WfARXR+lBgDyF5thTjbYRWnyGEvng9o+0lwCPA/cAIwi+E3cB3gZ7AEdHz1P6p/uOp674pWu8DfDPun4UWLVrydwEuBN5u4rW1wP9EefUeoHO0/Q7gzmj9OuCsaP2zQFW0/ing/Gh9QpSrD4me3wZ8L1pvKicPBDYBPweOAk4GrgG+E+37DeCBtFgdKIvWewCvARdGzz8M3JG27yzgZkLX1zOj3xdnRdf5MPDDVn52T0a/K34M/JTQ3/rktNfXAaOj9QWE30klwOXA28BU4HNR3n8yFW+0/3Lg2ui1TwErm4jhImBWtD40+pxPjfvflZaGi/pYS+K5+wIzOwS4Cvgi8FEzO9lDn+sqYLO73w1gZt8DngbGEn6J1JrZV6NDPUJI4KcDg9K2vxa9tzdwcfRe3H2tmb3VhlBLgRuAz7n7tugYO83sHWCZu6+OzvMOMNtD/8OXCQm/lDAI8nTCLwGAbcCnzWwJ8EvCHxUiIh3BCcXmUcAYYIaZAawH+ll48nlCowDAbwldISDk2jkZuXY48GrGOS6kkZzs7hvM7D3gfnd/AcDMVgJbzWw4oYW3d6NBu283s3fTNn0pOm7KLcD/Rds3EX5f3BOdY04UZ2v93aM+1mb2Z0JjTsrH3X2lmR1DKKh7R/l/NbDU3Z8Hno/e29ixH4h+J7zYTEzvA7ZH173WzDbQsEuJJIAKa0k0M+sM7OdhkMqXo+4ZDwKVhBaATAuix/7AaOB37n5HtO3n0TFnEgrdn2ec6whCQqzex3AnAB8ExgHz0rZn9jHcU/eCe22UZFPdsgZ56P6Cuy8zswuBG4GvAecB7yIikkVm1p/wbdsCQt5c30h+HEpoUe0G4O57CI0BRO/5vrs/Gz1v8N40jebkiNMwV74NzAReB54D9m/mEtLfNx74W9rzpYRuLkPYOxfXsO93LHmJ8AdIyiYz+wlwJyH2VPWceV0tqaHp2uwh4Jdm1p3ws9hOaP2WBFEfa0m6un5/AO4+D7gMOLCJ/bsRktgbwBrg7PQXzWxqtP10M+uRsX1z9HTSPsb6LHA98Ico8bWJmY0gbSCMmY0k3BFlAmEw4x1NvFVEpD0uIHxDdichPx4TtRQDdflxHbALOC1t+yFm1oeMXGtmXaKGikxN5eTG/C+w3MNtVXe24VpWElq4604BvEfo6pI17l4bNX4cHP1h8iBwpbvPyeZ5Ms75APAA4ZvVTwHHuft7HXU+2TcqrCUfzIiKzpSRhL5xKQPNrFu0fjyhf90m4I/AmWZ2o5kdH7UmbCUkwF7ALDP7cDRQ5EBCMf4y8AMz62ZmpYS+zUPMrLOZfdPMzm8iRouW7xB++VzXyGuZ+2c+P4PQHzBlPDDN3TcDMxp5j4hIW+zVEmpmHwK+D3ze3asIrcNrgQfN7KPRAMVjoxbqPxHuHjLNzE4CLnD3akKu/bKZfd/MjiM0MKS60aXnv6Zyckp66/HhwGAz6wkcA/Qws7HRa7sJg/1SDSzp5/hf4D/MbED0fCrwK3evIdQ8jeZiMzs6iqukic+uExk1U/THx+mEbykHELoYTgT2A3qmDYzMbBXf1Uj8DYNqpL9I9LPoCzxBaJVvyx8ckiPqCiL5YD9ggZn9lfA1mQFfSHu9BrjczDYRWndnALj7Y2b2ZcLXidMIg//mA5jZx4BfEH5R3Ap8w93dzD5BuBXSAuDPhK/0DiNMqPAJQheTBhPTmNlRhEE5owlfDd4HfD+K5xHgaKDMzMYRWlKGAZ8ws0XUt95cCPRy98qMa/+NhYkdagiDhURE2izKI+cSir9fAVWEAXD9gI+4+z8B3H1XlB9/Rch19xP6VgN8hdAocTvhG7rUnTh+Tch9FxO6rF3s7u9GhePRwDgz+31TOdnMzib0K77IzBa5+yrCQMprCEX17cDHgEOAN4HfEVrXP23hblEHA/9mZg+6+z1mNh6428yeIXzreXnUsv5xYFh0fQsIk+UMN7ODozg/Seg3fl/a59aVkPunACOjwtsJfZvPio75L2A28CKhe8tDhJw+izBgcoqZTXP31HH/TJisZ0tUiA8DzjazZYSWaKL3Z87ZsJIwOHMa0B3oZGaPEwYw5uoWtNIC089C8pmFaWBvc/eyHJ3vG+6u6cZFRApM1C2lxN0TOUGXmX2DcIeW5dHz7oRvG36SGpsj8VNXEMl3jXWz6JgTha9MX8rFuUREJHeiwZnvS3BR3R24AihL21wKbFVRnSzqCiJ5K+rfdj7hq72Pu/udHXzKpdFXlCIiUli6uPsv4w6iKe6+w8K06v8TDWxfTei2ckXz75RcU1cQEREREZEsUFcQEREREZEsUGEtIiIiIpIFKqxFRERERLKgYAYvlpaWellZWZveU1VVxeDBgzsmoBwqlOuAwrkWXUfyJPla5s6du87dkxlcB9mXnA3J/jm2la4lmXQtyZSka2kuZxdMYV1WVsacOW2bSbS8vLzN70miQrkOKJxr0XUkT5KvxcxWxB1Dru1LzoZk/xzbSteSTLqWZErStTSXs9UVREREREQkC1RYi4iIiIhkQVEX1hUVFXGHkBWFch1QONei60ieQrqWYlZIP0ddSzLpWpIpX64l5xPEmNkpwFXAOWnz3Z8FXAv0BX4P/Je710SvlQI/BDYBXYGZ7r4787jl5eWelL43IiJtYWZz3b087jgao5wtItJQczk7py3WZjYU6A0clbZtNDANOBv4GnAR8NW0t90J3OruM4ElwJU5C1hEpIgpZ4uItE1OC2t3Xwvcn7F5DPBZd5/r7rcDvwBOADCzo4GJ7v58tO8DwAwz65WrmEVEipVytohI2+S8j7W712Y8fyr1FWJkFbAyWj8RWJa276poNTtfmZ56Khx8MGzenJXDiYgUmqTk7OqXq5lzxBwWfWZRew8lItJhkngf6yOBn0brI4H1Ga9XAyMy31RVVUV5eX3urqioaLmj+8KFsGoVvPce9O3bnphFRNqksrKSysrK1NPSOGNpp5zk7NrttWx5aQuduhX1mHsRiUlrc3aiCmszOwCocvd50SYHtmfsVgLsNRBm8ODBbb9xeL9+9YX1fvvtQ8QiIvsmvZA0s3Uxh7NPcpmzrasBULu7toU9RUSyr7U5OzF/+ptZZ+DzwDfTNq8CBqTtY0AvYHVWTppqpVZXEBGRNsl1zk4V1r47t3eyEhFpi8QUfKjEtgAAIABJREFU1sAlwDXuvjNt22PAxLTnZcAuIDv3aOrXLzy+915WDiciUkRymrOtS1RY16iwFpHkynlhHbVgAFjatsuAl4EeZjbWzD5rZge4+xxgjZlNjnY9DbjB3XdlJRi1WIuINCspObtT1/DrSi3WIpJkOe1jbWZ9gPOjpxea2Y2ErxJ/kLHrIne/NVo/B/iOmS0n3E/18qwFlCqs1WItIrKXJOVsdQURkXyQ08La3auBX0ZLylXR0tR7VgCf65CAUl1B1GItIrKXJOVsDV4UkXyQpD7WuacWaxGRvKA+1iKSD4q7sNbgRRGRvKCuICKSD4q7sB4Q3RVqw4Z44xARkWZp8KKI5IPiLqyHDg2P774bbxwiItIsdQURkXxQ3IX1kCHhce3aeOMQEZFmqSuIiOQDFdagFmsRkYSrK6xrHHcV1yKSTMVdWA8eHB6rqqBWt3ASEUkqM4POYV3dQUQkqYq7sC4pgf79Yc8eDWAUEUk49bMWkaQr7sIa1B1ERCRP6M4gIpJ0KqyHDw+Pb78dbxwiItIsDWAUkaRTYV1WFh5XrIg1DBERaZ51i6Y136UxMSKSTCqsVViLiOSFzj3C6MXa7SqsRSSZVFiPGRMely+PNQwREWlepx7hV5YKaxFJKhXWqRZrFdYiIonWqWf4lbVn256YIxERaZwKaxXWIiJ5QV1BRCTpVFiPGgVdu8KqVbB1a9zRiIhIE9RiLSJJp8K6a1eYMCGsv/ZavLGIiEiTOvdUi7X8//buPD6usuz/+OfK1uxNmu77SltaKJSwyGIrOwIiPqCAoqJQcVfEDZCfDyoKAiq4FnlwQQQVkE3RIjuylaUFChRtS1dokjbN3iTN/fvjPkkm0ywzySRnJvN9v17nNTPnnDnnOjmTK1fuuc99RJKbCmuABQv845o14cYhIiI96rh4sUGFtYgkJxXWAPvu6x9ffTXcOEREpEcdXUEa1RVERJKTCmvobLFWYS0ikrQ6Ll5Ui7WIJCkV1gD77ecfX3op3DhERKRH7S3W6mMtIslKhTXAnDkwcqQfGWTr1rCjERGRbrRfvLinXl1BRCQ5DXlhbWbHmdkzZjY9Yl6emd1gZlea2S/MrDiWZQmTkQEHH+yfP/NMwjcvIpKqkilnZxYHhXWtCmsRSU5DWlib2TigEDgkatHPgSecc5cAdwO/inFZ4hx6qH98+ulB2byISKpJtpydNTILgNbq1kRtUkQkoYa0sHbOvYNPtB3MbCJwDnBvMOtB4H1mNr23ZQkPrr2wVou1iAiQfDk7q0SFtYgkt6yh3qFzrs3MImctBSqdcw3B8lYzWw8sAVp6WbYhoYG9613+8emnobER8vISunkRkVSUTDlbhbWIJLtkuHhxElAVNa8WmNjHsi4qKiooLy/vmJYvXx5fFKNHw4EHwu7d8MQT8b1XRCROy5cv78hXwOiw44lDaDlbhbWIhCXWnD3kLdbdcEBj1LwR+JYP+ljWYcyYMaxcuXJgkRx7LLz4Ijz4IBx33MC2JSLSi2XLlrFs2TIAzKwy5HDiEVrOVmEtImGJNWcnQ4v1FqA0al4RsLWPZYnXXkyvWDEomxcRGQZCy9kqrEUk2SVDYf0wMNHM8gDMLBuYCjzax7LEO/JIyM/3rdYbNw7KLkREUlxoObt9uL3WXa24NpeITYqIJFQY41i3XwVjAM65t4H7gGOC+ccCdzjntvS2bFCCy8uDk0/2z//yl0HZhYhIKkmmnJ2RleGLa+eLaxGRZDPU41gXARcGLz9uZmOC558GTjOzS4H3Assi3tbbssQ780z/+Oc/D+puRESSXTLm7JzxOQA0v9OcyM2KiCTEkF686JyrBX4RTJHzdwIX9PCeHpcNive+13cHefppWLsW9tlnyHYtIpJMkjFn54zPoXFtI81vN1Mwr2CwdiMi0i/J0Mc6uRQUwFln+ee//nW4sYiISBc544IW621qsRaR5KPCujvBcCrcfLMf11pERJJCR1eQt1VYi0jyUWHdnUMOgUWLoLJSFzGKiCQRFdYiksxUWHfHDD77Wf/8qqvAaVgnEZFkoMJaRJKZCuuefPSjMHEivPwy3H9/2NGIiAiQM9EX1rs3q5ueiCQfFdY9GTECLr7YP//e99RqLSKSBPJm5AHQtL4p5EhERPamwro3F1wAo0f7offuvjvsaERE0t6IaSMAaNrYRFtrW8jRiIh0pcK6N4WFcPnl/vnXvw4tLeHGIyKS5jJzM8mZlAN7YPcmdQcRkeSiwrovn/oUzJ7tbxazfHnY0YiIpL28meoOIiLJSYV1X3Jy/MggAJdeClu3hhuPiEiay52RC0DjusaQIxER6UqFdSxOPx1OOQV27YLPfS7saERE0lrebN9i3fiGCmsRSS4qrGNhBj//ue9zfddd8Kc/hR2RiEjaKlhYAEDdy3UhRyIi0pUK61hNmQJXX+2fL1sGGzaEGo6ISLoq2M8X1vUv14cciYhIVyqs43HhhXDaab5LyNlna5QQEZEQ5M3MIyM/g+atzbTsUB4WkeShwjoeZnDTTTB5sh/b+pJLwo5IRCTtWIZRsCBotX5FrdYikjxUWMerrAxuvRUyM+Gaa+CWW8KOSEQk7RQuKgSgdmVtyJGIiHRKSGFtZu9NxHZSxlFHwfXX++fnnw/PPBNuPCIicRgOObv4sGIAap6qCTkSEZFOWfGsbGYnApcB4+ksyjOAcUBeYkNLcp/5DLz8Mvzyl3DqqfDYYzBvXthRiYh0GM45u/hwX1jvempXyJGIiHSKq7AGfg98H3gRaAvmGXBOIoNKGddf70cHeeABOO44eOIJmDYt7KhERNoN25ydPzefrJIsmrc007SpidwpuWGHJCISd2G9zjl3XfRMM3s9QfGkluxsuOMOOOEEX1Qfeyw88ghMmhR2ZCIiMIxztmUYxYcVs+OBHex6che5Z6mwFpHwxdvH+tNmdkY389+XiGBSUn4+3HcfLF4M//kPLFkCGzeGHZWICAzznD1yyUgAdq7YGXIkIiJevIX1EcCvzWxPxNQG/GKggZjZfDP7mZl92cx+bmYHBPPzzOwGM7vSzH5hZsUD3VfCjRwJ//ynL67/+19497th3bqwoxIRGdY5e9QJowDY8Y8dOOcGazciIjGLtyvIFcAyYCWwJ5hnwLkJiOX3wGnOuS1mNhX4BzAf+DnwgHPu9uBCnF8BZydgf4lVVgb/+heceKIfJWTJEv96n33CjkxE0tewztmFiwrJHpdN85ZmGtY0dIxtLSISlnhbrP8O3O2cW+eceyuYNgA3JiCW+UBR8LwRGGlmE/EX2dwbzH8QeJ+ZTU/A/hKvpARWrIAjj4TNm31xvWZN2FGJSPoa1jnbMoxRxwet1g/sGIxdiIjEJd7C+jbgm2b27ohpCfDNBMTyF/xXlkXAR4DPA0uBSudcA4BzrhVYDyxJwP4GR1GRHyXk6KPh7bdh6VJYvTrsqEQkPQ37nF12chkAFXdVDNYuRERiFm9h/S3gS8DvgN9GTJ9KQCyfAZqA54A659wdwCSgKmq9WmBiAvY3eAoK/AWNJ5wAFRXwnvfACy+EHZWIpJ9hn7NHnTyKjNwMap6soWlz02DtRkQkJvEW1j8EpjjnpjvnZgTTdOADCYglD/gjcCvwYzM7FnD4rxgjjQBaot9cUVFBeXl5x7R8+fIEhDQAeXnw17/CKafAjh2+Bfuxx8KNSUSSxvLlyzvyFTB6kHYz7HN2VmEWo97ru4NU3lHZr22IiPQl1pwd78WLm5xztd3MHx/ndrpzC3CWc67azBw+YX8JKI1arwjYGv3mMWPGsHLlygSEkUC5uX6c67PPhjvvhGOOgZ//HC64IOzIRCRky5YtY9myZQCY2WBVhGmRs8d+cCyVd1ay/fbtTP7i5IRsU0QkUqw5O94W689EzzCzM4Efx7md6G2MBhY556qDWd8FaoC3gIlmlheslw1MBR4dyP6GVE4O3H47fPnL0NoKy5bBF7/on4uIDK60yNmjTh5FRn4GNU/V0LC2YTB3JSLSq3gL61Fm9i0AM8s3s18Dv8b3sRuIHUCTmUXesrAKWA3cBxwTzDsWuMM5t2WA+xtaWVlw3XVw003+bo3XX++7hmzeHHZkIjK8pUXOzirMYuxZYwHY9uttg7krEZFexdUVxDl3spkdZmY3ACcAbwOLgF0DCcI512Zm7wcuN7PngXHAV51zNWb2aeBqM1uE//py2UD2FapPfALmzoUzz4THH4cDDoDf/hZOPjnsyERkGEqnnD3hggm8/X9v8/Zv3mbGd2eQkRNvu5GIyMD1Wlib2bt7WPQW8CZwDTANeD/w5YEE4pxbRTdXqjvndgLDp1PyEUfASy/BRz8K//iHv7jxoovgyithxIiwoxORFJbOObv40GIK9iug/uV6Ku+qZOyHxg51CCIiWG+3gTWzp/E3AagG2npaDZjgnAu1KiwvL3dJd/Fib9ra4Jpr4JJLYM8eWLQIbr0V9t037MhEZIiZ2fPOufIEbCetc/aWX2zhzc+8SVF5EYufXYyZJXT7IiLQe87u67uyK4A5zrlpEUM1RU/TgQ8nOuhhLyMDvvY1ePJJmDkTVq2Cgw7yo4b08s+OiEgv0jpnj//YeLJHZ1O7spbqR6r7foOISIL1Wlg75/7mnNve10acc39JXEhp5tBDfdeQj38cmprgs5+F978fqqLvsSAi0rt0z9mZ+ZlM+oK/nnLjVRtDjkZE0pGu7kgGRUVw881+WL6RI+Gee3zXkIceCjsyEZGUMukzk8goyGDnP3ZS80xN2OGISJpRYZ1MPvhB3yXk8MNhyxZ/Q5lPfQp2DegCfhGRtJFdls3kz/ubxKy7ZF3I0YhIulFhnWymTYNHH4UrrvBjXi9fDgsWwH33hR2ZiEhKmPK1KWSVZFH9UDU7HtwRdjgikkYGXFib2YFmdouZfcXMpg88JCErC771LXjxRd8He8sWOPVU+PCHoaIi7OhEJIWlQ87OLs1mytenALDuG+twbbogXESGRiJarC8EjnXOXQsc3n4rW0mABQv8qCHXXQd5eX44vvnz4Te/0cghItJfaZGzJ39hMjkTc6h7vo5tN+lujCIyNBJRWH8OOCh4fhuQnYBtSrvMTPjyl+Hll32f66oqOO88eM974I03wo5ORFJPWuTszPxMZl07C/Ct1i1VLSFHJCLpIO7C2sxOMrMPB88PwY+ZugX8bW6dc7oMezDMmgUrVsDvfw+jR/t+2PvvD9/+NrToD4aIdC+dc/bYD42l5OgSWne0su6bupBRRAZfXIW1mX0P+DNwOoBz7lngo2Z2ziDEJtHM4CMfgddfh09+Epqb4X//F046CXboAh0R6Srdc7aZMeenc7BsY9uN29j58M6wQxKRYS7eFuvDgcnA4xHz/g78IGERSd/KyuDXv4aHH4Zx4+Bf/4KFC+FPf1LfaxGJlPY5u2B+AdMumwbAG594g9ba1pAjEpHhLN7CeoVzrhqIrN7OSGA8Eo+lS+G55/y419u2wYc+BMcd5/tji4goZwMw9ZtTKTywkKYNTaz7mrqEiMjgibewXmdmlwBzzewsM7sD+CxwVeJDk5hMmQKPPw6/+hWUlvrW6wMOgM98RkPziYhyNpCRncG8387Dso2tv9zKjn+o65yIDI64Cmvn3G3AQ0Am8FGgGjjeOfezQYhNYpWRAcuWwZtvwuc+5/ti/+IXMGcO/OhHvi+2iKQd5exOhfsVMv1/pwPw2rmvsXvr7lDjEZHhKe5RQZxzTzvnLnTOvdc590ngYTM7ZhBik3iVlcENN/jboh9/vL8V+kUXwX77wf33q/+1SBpSzu409WtTKT22lJaKFtactYa21rawQxKRYSart4Vmdg9Q2Mc2RgOVwL8SFZQM0IIF8MADvpi+6CJYuxZOOQVOPBF+/GOYOzfsCEVkEChn984yjfl/mM/KA1ay6/FdbLh8AzOvnBl2WCIyjPRaWAOrg3XeANqAC4AXgZUR68wARg1KdNJ/Zr6YPv54+OlP4YorfLG9337wla/48a9HjAg7ShFJLOXsPuSMzWHfP+7LS0e/xMbvb6RwcSFjzxgbdlgiMkz0VVj/ELDgqnLMbJZz7vLolczs7sEIThIgJ8e3Wp97LlxyCdx0E/zgB741+5Zb/E1mRGS4UM6OQcmSEmZeNZN1X13H6x99nbwZeRQdVBR2WCIyDPTax9o5t6s9QQfGmVlm5DpmthQ4ZBBik0QaMwZuvBGefNJf1Pjyy3DwwfDDH8KePWFHJyIJoJwduylfmcL488bT1tjGy6e9rIsZRSQh4r148S/AKjP7pZn90MzuA1YA1yc+NBkU73oXvPgiXHihHy3ka1+Dww6Dp54KOzIRSTzl7B6YGfv8Yh9GHjWS5i3NvHLaK+ypVyODiAxMvMPtrQCOx1/4Mhc/dNMZzrnvD0JsMlgKCvxwfPffD5MmwcqV/iYz554LW7eGHZ2IJIhydu8yRmSw4I4F5M7IpXZlLa+e8SptzRopRET6z9wAh2AzswzgQOfc8wkJyOww4EjgP8ATQD1wNVALlAJfd87VRL+vvLzcrVy5Mnq29KWuzve5vuYa2L3bF92XXgpf/jLk5oYdnUhaMLPnnXPlQ7Qv5ewoDWsbePGIF2mpbGHsOWOZ//v5WIaFHZaIJKnecnZchbWZ/V83s8cAe5xz7+9nfJHbPw+Y5Zy7LGLezcADzrnbzexE4GPOubOj35tMSTolrVsHF18Md93lX8+cCddeC6ed5kcYEZFBM1iFtXJ27GpW1rDqPavYU7eHSV+YxOwfz8aU+0SkG73l7Hj7WM8BNgJvRUwjgAG3fJjZEuDDwLci5k0EzgHuDWY9CLzPzKYPdH8SZeZMuPNOePBBPw72unVw+ulwwgmwZk3Y0YlI/yhnx6i4vJiFf12I5Rhbrt/C+svWM9BvdEUk/fQ13F60851zb0TOMLN5wMcTEMt1wL+BG8xsNvAdYBpQ6ZxrAHDOtZrZemAJsCEB+5RoxxwDL73k+2BffjmsWOGH5DvzTPj0p+Goo9SCLZI6lLPjUHpMKfveui+vfuhVNl65ERzM+N4MtVyLSMzivXjxjW5m7wbOH0gQZjYXOBC40Tn3OfwdwR4ApgBVUavXAhOjt1FRUUF5eXnHtHz58oGElN6ysuDzn4c33/SjhzgHt90GS5bAjBnwjW/426arNUek35YvX96Rr/B3Q0w45ez4jfmfMex7276QCRu/v5H1l6rlWkRiz9nx9rFeD0S+IRMYB/zeOXdBP2PFzD4ALHfOjQ5e5wEVwBXA/zjnDo1Y9wXgVufcNZHbSLb+esPKxo1+DOybb4YtWzrnz58PZ5/tp9mzw4tPJMUNYh9r5ex+qrijgjVnrcG1OiZ/aTKzrp2lCxpFBEhsH+sVwCeA84LpI8D+A0nQgSwiuqU45xqBN/F/BEqj1i0CNCbcUJo6Fb7zHV9gP/qob8UuK4PXXvPdRebM6bzZzLp1YUcrIp2Us/tpzP+MYd/b98Wyjc0/3syac9bQtltD8YlI7+ItrG9xzj3inHs0mB4HLOhfNxCrgZFmNi5iXis+GU8MWkMws2xgKvDoAPcn/ZGRAe9+t+9/vW0b/O1vfuzroiI/FvbXvgazZsHixfC978Eb3X0LLSJDSDl7AMZ8YAz7/31/Mosyqbi9gtUnrqZ1V2vYYYlIEou3sD6xm3nvAH8eSBDOudeBvwNnAJhZCZAD/AG4DzgmWPVY4A7n3JbutiNDKDsbTjoJfvc7eOcduOMOOOccX2S/+CJcdhnMmwcLF/pW7eeegza19ogMMeXsASo9ppQDHjuAnPE5VD9SzQvveoGGtQ1hhyUiSSqmPtZm9jngu/iv9LqzwjnXXQKPPRCz0cBPgJX4Fo4bnXNrzKwUf7OBDcB44JvOubro9yd7f7200dTkRxK54w64+26oru5cNn48nHIKnHoqHHss5OeHF6dIEkl0H2vl7MRr3NDIyye/TMOaBjJHZrLvrftS9t6ysMMSkRAk5AYxZrY/8DHg+qhFzc65bQMLceBSLUmnheZmeOghuOceuPde2Ly5c1lurh/a79RTfbE9aVJ4cYqEbDAuXlTOTrzW2lZe/9jrVN5VCQYzvjuDqd+YqosaRdJMIu+8mO2ca+lm/lTn3MYBxDhgqZik04pzfoi+e+/103PPdV2+eDGceKIfJ/vww6G4OJw4RUIwiKOCKGcnmGtzvHXlW2z41gYARp04inm/mUfOuJxwAxORIdPvwtrMLgb+5pxbE7y+FH/Vd5fVgPc455YmJtz+SdUknba2bYP77/dF9ooV0NjYuSwjw9+U5qij4Mgj4V3vgsmTdWMaGbYSVVgrZw+dyvsqef1jr9O6o5XssdnM++08yk5U1xCRdNBbzu7rzoszgJKI14uBAuBtOsdGNXw/OpHYTZgA55/vp8ZGePhheOQRePxxP8LISy/56YYbOtc/7DA49FA/lZdDYWGohyCShJSzh8joU0Zz8OqDee3c16h+uJqXT3qZSV+cxMzvzSSzIPp/GRFJF/F2BVkMvOKca46av297C0lYUr31QyI0NMCzz8ITT/jp2Wdh586u62Rk+BFH2gvtww7zo5Bk6g+apJ5B7AqinD3I3B7Hxqs3sv5b62EP5M7MZe6Ncyk9Ono4bxEZLgbSYh1tP+fcC1EbnwvMAkJN0jKM5OfD0qV+Aj9M35tvwjPP+Onpp2H16s7pxhv9ekVF/kY1ka3aEyeqC4mkM+XsQWaZxrRvTqP0uFLe+OQb1K+uZ9Uxqxj/yfHM+uEsskuzww5RRIZQrMPtHYX/+nAZsDxq8Sjgeufc1MSHF7vh0vohMWpshBde8EV2e7G9adPe640d6y+MPOgg/7h4MUybpmJbksogDLennB2CtpY2Nl29iQ1XbMA1O7JHZzPjyhlM+MQELFM5R2S4GPCoIGa2H/A7fP+9qO/kaQb+4Jy7YqCBDsRwTNISp23bOovsZ5/1N6qJHEe73ahRnUV2e9E9c6bvXiISgkEorJWzQ1T/Wj1rL1zLrsd2AVB4YCFzbpjDyCNGhhyZiCRCosaxHgmc6Zz7dSKDS5ThnKSln5yD9et9y3b79PzzUFm597rFxbBokR+NZNEiPy1YAAUFQx+3pJ1BGsdaOTtEzjkq/lTBf7/6X3Zv2g3AmA+OYcZ3ZpC/j26OJZLKEjaOdS87ONI598SANzQAwz1JS4I4529UE11sb+vmfhlmMGdO12J7//1h6lR1JZGEGqyLF3vZn3L2ENnTsIeNV21k09WbaGtqg0yYcN4Epv2/aeROzg07PBHph4GMY30P0NeYZgYscM6N7X+IA5cuSVoGydtv+xvYrF7d+fjaa9Dauve6JSW+wG4vuPff349Qolu0Sz8lcBxr5ewk1bSpibeueIttN2+DPWAjjEmfncSUr05hxPgRYYcnInEYyKggbwDVwBY6x0Dda/vAKf0PTyQJjB/vpxNO6Jy3e7cvrtuL7fapshIee8xP7dpbtyNbthctgilT1LotQ0k5O0nlTvHD8E356hTWX76eitsr2HzdZrb8bAsTPjmBKV+dQt70vLDDFJEB6qvFejRQ65zb3etGzMqcc1WJDi4e6db6ISFxzrduR7Zsr1oFr7/ee+v2fvv5PtvtU5nu0CadEthirZydImpfrOWt77xF5V3BNR+ZMO7D45j69akU7KtrO0SS2VD0sT4/7Atk0j1JS8jaW7cji+321u3ujBvnC+x99+1acI8aNbRxS1IIoY+1cnaSqH+1no1XbeSdW9+BPX7eqJNGMfmLkyk9vhTTN14iSWfQCmszmwp8Fviccy7Uf7GVpCXptLdur1oFr7wCr77qpzVroL6++/eMH99ZZEcW3aW6i9twNlSFtXJ28mpc38imH27i7Zvf9hc5Avnz8pn0hUmM/+h43SZdJIkktLA2s0zgfcAFwHHAO8F2Jg000IFQkpaU0dbmb2bTXmhHFtwNDd2/Z8KErgV3+6QW7mFhMAtr5ezU0lzZzLYbt7HlZ1to3uLvRJ9VksX488Yz4fwJ6iYikgQSNY71DOB84DygDHgL+ALwD+Bdzrl/Jybc/lGSlpTX1gZvveUL7MiC+7XXei64x43rWmi3T2PG6KLJFDJI41grZ6ewtpY2Ku+qZPNPNlPz75qO+cWHFzPhggmMPXOsWrFFQjKgwtrMzsS3dBwN/Af4NfBb4BLn3JcTHGu/KUnLsNXWBhs2dC2016zxjz11KSkr8wX2/PldC+6JE1VwJ6FEFtbK2cNP7fO1bF2+le23bmdPne+InVmcybhzxjH+E+MpKi9SX2yRITSQ4fYAGvG3wP03cL5zbm2w0YFf9SgifcvI8LdcnzkTTj21c357l5I1a/aeqqrg8cf9FKm4uPsW7ilTdEv34UM5e5gpOqiIub+ay6xrZ1FxewVbb9xK7TO1bP3lVrb+cit5++Qx7iPjGPfhceTN1JB9ImGKpyvIVHwryHTgbmCJc+7zwTJziRheZADU+iEScM7fSTK62H71Vdixo/v3FBTs3bq9774wfTpk6uvmwTZIXUGUs4exupfr2HbTNrb/cTst21s65hcfXsy4j4xj7AfHkl2WHWKEIsNXoi9ezALeD3wYeAW4CTjeObd8oIEOhJK0SB+cg4qK7lu433mn+/fk5sK8eXt3K5k1C7L1RztRBvniReXsYayttY2dD+7knVveofKuStoa/Igilm2UHlfKmDPGMPq00WSP0u+rSKIM5nB7s4HPAxc450K9n7OStMgAVFV19t2OnLZs6X797GzYZ5+9W7jnzIERuj1zvIZwuD3l7GGsta6VqrureOeWd9jxzx3ga2wsyyg5usQX2e8fTc6YnHADFUlxQ3GDmI8653434A0NgJK0yCDYtcsX3NFF94YN3a+fmQmzZ/thARcu7HycM0ct3L0I4QYxytnDXPP2Zir/WknFXyrY+dDOjpvPkAElS0s6WrJHTNQ/wiLxGvTX7Ya2AAAgAElEQVTCOlHMLAd4Dviic+4RM8sDrgZqgVLg6865mu7eqyQtMoTq6/1t3KNbuP/7X9/lJFp2Nsyd64vs9mnBApgxQ324GfrCOpH6m7eVs4dOS1ULlXcHRfaDO3Etnb+jhQcVMvrU0ZSdWkbhgYUaXUQkBgMdFWQofRV/oU27nwMPOOduN7MTgV8BZ4cRmIhEKCiAgw7yU6TGRnjjjc6hAV95xU/r13c+j5SX13mHyciie/JkDQuYOpS3k1x2WTYTPjGBCZ+YQMvOFqruraLijgp2rthJ3fN11D1fx4ZvbyBnYg5lp5RRdmoZpUeXkpmvf3pF4pU0LdZmdjgwD7gc+DiwFlgPlDrnGoILcHYBC5xzG6Lfr9YPkSRWV+e7k7QX1+1Fd099uIuL9+5OsnAhjB07LAvuVG2xHkjeVs4O357GPVQ/VE3lvZVU3VfVcadHgIy8DEqOLmHUiaMYdfwo8ubkqTVbJJD0LdZmVgCc4Zy7yMwuD2YvBSqdcw0AzrlWM1sPLAE2hBKoiPRPYSEcfLCfIu3c6buQRBbcL78MlZXw1FN+ilRW1llkL14MhxziRytRd5Ihp7yd+jLzMik7uYyyk8twzlH3Uh1V91ZRdW8VtStr2XH/Dnbc74fozJ2eS+nxpYw6fhQlx5SQXaJrJkS6kxSFNfB14AdR8yYBVVHzaoGJ3W2goqKC8vLOfx6WLVvGsmXLEhmjiCRaaSkccYSfIm3f3rVlu/15VRU8+qif2hUUQHm5L9oPOcRPU6cmfcv28uXLWb68Y8S70WHG0k8DytvK2cnFzCg6sIiiA4uYfvl0dm/bzY4HdrDznzvZsWIHTRua2LZ8G9uWb4MMKD60uKPQLjqkiIws3WBKhrdYc3boXUHM7ARghHPunuD1BvxXiuXAmc65QyPWfQG41Tl3TfR29LWiyDDnnO868sorvlV75Up49tnuRygZO7ZroX3oob6IT1Kp1hUkEXlbOTt1uDZH7Qu1vsj+5w5qnqzBtXbWDpnFmYw8aiSl7yml5D0lFC4qxDKT+x9bkYFI9q4gFwMHRvTdKsXfJey64HmkImDr0IUmIknDzF/UOHkynHhi5/zt2+G553yR3f64fTvcf7+f2i1Y0Nk6fsQR/hbxSd6qncSUt9OIZRjF5cUUlxcz7ZJptNa2Uv1ItS+0/7GDxjcbu3QbySrNYuS7OwvtgoUFWIZ+1yQ9JEOL9RggciDNp4CLgGeANcAY51yjmWUDdcBM59xeVzyp9UNEAN+yvX69L7CffRaeeQaefx527+663rhxvsA+8kg46ig48MDQ+mqnYIv1gPO2cvbw0bSpiepHqql+2E9NG5q6LM8enc3IJUGhvbSE/Pn5KrQlpSV1i7VzriLytZntASqccxvN7D7gGOA+4Fjgju6KahGRDma+NXrmTDjrLD9v925fXD/5ZOf0zjtw551+AigpgaOPhmOP9dPs2WrR7oHytkTKnZLL+HPHM/7c8QA0bmjsKLKrH65m9+bdVN5RSeUdlUDQon3ESEYe6aei8iIyRqiPtgwPobdYR2vvqxfcaKAUf6OBDcB44JvOubru3qfWDxGJmXOwdm1nkf3II7BuXdd1pk6FE06A00/3Bfcg3qo91Vqso/UnbytnpwfnHI3/7Sy0dz2+i92bu357ZCOM4oOLOwrt4sOLyS7VqCOSvFLmzosDoSQtIgOybh3861/w4IP+sSpicIviYjjlFPjAB+DkkyE3N6G7TvXCuj+Us9OTc47dG3ez64ldHVP9K/V7rVewsMAX2UcUU3xYMXmzNI62JA8V1iIi8Whrg5degvvu811FVq3qXFZSAmefDeed54f5S8AfexXWks5adrRQ81RNR6Fd82wNrrlrbZJVlkXxob7ILj6smOJDiskaGXpvVklTKqxFRAbiP/+Bu+6C226DF17onL9wIVx0EXzkI5Dd/6+uVViLdNrTtIfalbW+yH6qhpqna2jZ3tJ1JYP8efmdhfahxeQvyNd42jIkVFiLiCTKyy/DzTfD73/v7xAJMGsW3HQTLFnSr02qsBbpmXOOpg1N1Dzji+yap2uoe7Fur1btjIIMig8upujQIorK/ZQ7LVddSCThVFiLiCRaczPcfjt85zvw5pu+S8gtt8A558S9KRXWIvFp291G3Ut1vtAOCu6m9U17rZdVltVRZBcd5B9HTB6hYlsGRIW1iMhgaWmBb30LrroK8vNh61YYOTKuTaiwFhm45u3Nvsh+poa65+uoXVlLS2XLXutlj83uLLaDgnvExMEb9UeGn6Qex1pEJKVlZ8PFF8NPfwr19f7mNAccEHZUImknZ2wOo08dzehTRwOdI5DUPl9L7crOqWV7Czv+toMdf9vR+d4JOR1FduHiQgoPLGTEJLVsS/xUWIuI9FdLi+/+ccklvqheuNBPIhI6MyN3Wi6503IZ84ExQGd/7chCu/b5Wpq3NVN1bxVV93YOs5k9OpvCAwo7pwMLydsnTxdISq9UWIuIxKO11Y8McuedcOutsGmTn3/kkX7UkCylVZFkZWbkzcgjb0YeY88cC4BrczSua+wotOteqqPuxTpaKlvY+eBOdj64s+P9GbkZFOxf0FFoFx5QSOF+hWQWZIZ1SJJk9BdARKQnra3+wsTVq/303HPw73/71ul2c+fCpZf6ixYz9cdVJNVYhpE/O5/82fmMO2scEHQj2bybuhfrOgrtupfqfGv3s7XUPlvbuYEMyN8nv2uxfUAhOWNzQjoiCZMKaxFJb875YfPefNPf5nzt2s7nb7wBu3fv/Z45c+C44+Css+CIIyBDXw2LDCdmRu6UXHKn5DL6faM75rfsbKFuVV2Xgrt+TT0NrzfQ8HoD22/b3rFuzoQc37q9f6F/XFRI/tx8MnKUL4YzFdYiMvw1N8PGjbBhQ+e0fr2/8cvatVBd3fN7p0+H/ff306JFvpCeMGFo4haRpJJdmk3p0lJKl5Z2zNvTtIeGNQ0dxXbti7XUr6qneVszzdua2fmPzq4klm3kz8+ncFFhl6J7xHiNSjJcqLAWkdTmHOzcCVu2+GnzZnjrra4F9Natfr2eFBf7Vuh99vFT5PM4h84TkfSSmZtJ0eIiihYXdcxzbY6m9U3Ura6jfnU9davrqFtVR9N/m6hfXU/96vou28gem93Zsr1/IQWLCiiYX0DGCLVupxoV1iKSvJqb4Z13Ogvm9uI5sojeuhUaG3vfTkYGTJ3qW58jp9mzfRE9dqy/wYuISAJYhpE3K4+8WXmMOX1Mx/zWulbqX/GFdd2qzqK7ZfveF0qS6W/bXrh/YZcW7pyJORoGMImpsBaRodXaChUVvmB+++3Ox8jn7Y87dvS9PYCiIpg0qXOaNg1mzOgsoCdP9uNNi4iEKKswi5GHjWTkYZ3fhLWPt123qq6zhXtVHY1vNtLwagMNrzaw/Y+dfbezyrK6tG4XLiokf998MvN08XQyUGEtIgPjnO+jXFnpC+bKys7n7QV0ZPFcWdl7t4xIGRm+Nbm9YJ48uWsB3f66qKjvbYmIJKHI8bYjL5Tc07CH+lc7W7fbi+7WqlaqH66m+uGIa0OCkUnaL5Jsf9Tt24eeCmsR6aqpqfsiuad5lZWwZ0/s2zeDMWNg/HgYN673x7IyDWEnImkpMz+T4oOLKT64uGOec47dW3ZTv6q+S//thjcaOkYmqfhTRcf6WaVZHUV2R9/tBQVq3R5EKqxFhiPn/FjLO3bEP/XVX7k7xcUwerQvmEeP7vo8ulgeM0Y3URER6QczI3dyLrmTcyk7uaxjfsfIJKvrOoruulV1tFa1suvRXex6dFfnRtpbtxd1Lbh1C/fE0F83kWTW0gK7dvlRL6qr/WOsBXJLS//2mZ29d3HcXcHc/lhWBiM0VJSISFi6HZnEOZq3NncU2e0Fd5fW7dsjWrdHZXXpRtLRdztXrdvxUGEtMpja2qCmxhfF7YVxd897WlZf3/c+epKXB6NGxT8VFGiEDBGRFGdmjJg0ghGTRlB2UlTr9qsNXQvuVXW07uim73Ym5M/N36vgzpmgkUl6osJapDfOQUND/AVx+7RrV+wX6nUnIwNKSjqn0tLYiuPSUl9Yi4iIRMjMzaTooCKKDurauh3Zd7u94G5Y20DDGj/xx85tZJVldRTZHRdM7qtxt0GFtQx37YVxZLEb79TaOrAYiov3Lo57eh29rKhIrcciIjKoeuy73RiMTBJVcLdWtVL9UDXVD0W1bs/Lp/CAQl+4Ly6i8MBCsorTq9RMmqM1sw8A1wLFwC3AV5xzrWY2GvgeUA1kA193zvWz86iknPaL8LoreHftGprCOC8v/oK4/XVxsS7Uk2FJOVtk+MvMy6S4vJji8qiRSTbv7tJvu25VHY1rI8bd/kPnuNt5c/IoOqiIwsWFvtheXEh26fC9r0BS/MU3s6nA+4EzgAXAL4BNwDXAn4FvOOeeMbMLgSuAb4YVq8Spt8I41imeody6014YxzuNHOkfdWGeSBfK2SLpy8zInZJL7pRcRp8SNe72K/XUvVRH7fO11L3gx95ufLORxjcb2X5bZ7GdOyO3s9gOHnNG54RxOAmXFIU1MA34hHOuFXjezBYC7zGzJ4C5zrlngvXuBd4ws+865wZwVZfEzDmoqxtYi/FAC+P8/L2L3XiKYxXGIommnC0iXWTmZ1J8SDHFh3S2brc1t1G/pp665+uofaGW2udrqV9VT9P6JprWN1Hxl85RSUZMGUFReRHFhxZTfFgxReVFZBak3ogkSVFYO+cej5q1BSgCjgbWR6y3JbgKtRx4dMgCTGW9FcaxTLt2JbYw7k+rcc7w+C9WZLhQzhaRWGTkZFB0QBFFBxQx4ZMTAGhrbaPhtQbqXvAt27Uv1FL3Yh27N+1m96bdVN5VGbwZCvYroPiw4o5iO39uPpaR3NcdJUVh3Y2DgR8Cy4CqqGW1wMQhjyhMDQ1QVdX9iBOxFMdtbQPbf0FB/1qLVRiLpAvlbBGJSUZWBoX7FVK4XyHjPzYeALfH0bC2gdpna6l5poaap2s6bnZTv6qebb/aBkDmyKBV/LBiig8vZuQRI8kqSq5SNrmiAcxsNlDhnFtlZg6Ivg3cCGCvC2EqKiooLy/veL1s2TKWLVs2qLH2S/vtotunqqreX1dW9u9OeJHaC+P+thhnD9+LDETCsnz5cpYvX97+cnRv6yazYZ+zRWTQWaZRML+AgvkFHcX2noY91D5fS83TNR3FdvOWZnau2MnOFTv9GzOhaHERJUtKKFlawsgjR5I1cnBK21hztrmBjLGbYGaWCfwAuMw5t9vMvgm8xzl3fLDcgKZg3r8j31teXu5Wrlw55DED/uK8bdvg7bf9Y/TziorOgrk/N/zIyfF3uIseeUKFsciwYGbPO+fK+14zuaRszhaRlNS0uYnaZ3yxXf14NbUrayGyt2oGFB5QSMnSEkadOIqSd5cMytjaveXsZGuxvhi4xjm3O3i9ArgwYvl0oBkYumy8Zw9s2QLr18O6dZ2PmzZ1Fs+1tbFvL/J20e1TWVnv83QnPBFJTsmXs0Vk2Gofa3vM/4wBoLW2lZp/11D9aDXVj1ZT+1wwGskLdWy+bjOZhZmUHldK2clljH7/aLLLBr+hMWkKazO7BJ9888xsJrAUeAzYZmYLnXOvACcCNzjnmgclCOfghRfgySf94wsvwOuvQ0sfQ7COGAHjx8OECZ1T++vx42HcuM4iubBQRbKIpLykyNkiktayirIYdcIoRp0wCvDdR2qeqmHngzupur+K+pfrqbyrksq7Kln7mbWM+cAYplw8pctdJxMtKbqCmNllwHeiZr/unJtvZtOAy4ENQCFwaTDEUxcD+lpx1y740Y/gpptg8+a9l48fDzNnwowZfpo5E6ZO7SyiS0pULItIv6VaV5DQc7aISAyaNjZR9bcqKu+q9P2yHWAw+aLJzLp6Vr9HGOktZydFYZ0I/U7Sa9bASSfBxo3+9aRJcMIJUF4OixfDwoW+K4aIyCBJtcI6EVRYi8hQanqric0/2cyWG7bgWh1TvzGVmd+f2a9tpVIf66HV2AinneaL6vJy+OEP4d3vhozEd3QXERERkXDkTstl9nWzGfXeUaw+YTUbr97IhE9NIG96XkL3k94V5C23wH/+AwsWwGOPwdKlKqpFREREhqlRx45i7FljoQ22/2F732+IU3pXkXfe6R8vvhjyEvsfi4iIiIgkn9Gn+2Goqx+rTvi207ewdg6eeso/P+64cGMRERERkSFRfGgxAHUv1iV82+lbWFdV+dFARo6EibrbroiIiEg6GDFpBJZjtFS0sKd+T99viEP6FtYbNvjH6dM1VJ6IiIhImrAMI3dqLuBHC0mk9C2sq6r849ix4cYhIiIiIkMqe5y/C2NLZR83AYxT+hbW1UGH9ZKScOMQERERkSGVXeoL69ade92/akDSt7Detcs/jhwZbhwiIiIiMqSySv2tXFp2qsU6MWpq/GNxcbhxiIiIiMiQai+s1WKdKCNGwIQJUFYWdiQiIiIiMoRyxuWQMz4Hy0zsABbpe0vzz3/eTyIiIiKSVqZdMo1pl0xL+HbTt8VaRERERCSBVFiLiIiIiCSACmsRERERkQRQYS0iIiIikgBpXVgvX7487BASYrgcBwyfY9FxJJ/hdCzpbDidRx1LctKxJKdUORYV1sPAcDkOGD7HouNIPsPpWNLZcDqPOpbkpGNJTqlyLGldWIuIiIiIJIo558KOISHMrAJ4K863jQYqByGcoTZcjgOGz7HoOJJPMh/LNOfcmLCDGEr9zNmQ3OcxXjqW5KRjSU7JdCw95uxhU1iLiIiIiIRJXUFERERERBJAhbWIiIiISAKosBYRERERSYC0LKzNLM/MbjCzK83sF2ZWHHZMA2Fm3zUzF0yrwo4nVmZ2nJk9Y2bTI+al5Lnp7liC+SlzbszsA2a23syqzOwnZpYVzB9tZr8ys6vM7Dozyw471t70dBzBslsizsfdYcYpPYs1D5jZSWb2UzP7uZmdMdRxxiKWYzGzkWb2ZzOrMbMXzeywMGLtS7z52cw+ZGaPDFF4cYnnWIJ1v2RmH0vGcxPH78slZnZZ8PhdM7OhjjUWPf09jVoneX/3nXNpNwE3Ax8Knp8I/DHsmAZwLAXAL4Fjg2lO2DHFGPc44HTAAdNT+dz0ciwpc26AqcDvgIOAjwL1wMXBsoeBQ4PnFwLfDzvefh7HVOC6iPMxJex4NfV4HvvMA8A+wAtARvD6MWBh2LH381iuBM4OPrf3A+8ABWHH3p9jiVh3IrAaeCTsuAdyLEApcA8wM+yYB3IswCnAjRGv/6/9Pck09fT3NGqdpP7dT7tRQcxsIrAeKHXONQStWbuABc65DaEG1w9m9qXg6S+dc02hBhMnM8sA9gAznHMbUvncRB9LMC9lzo2ZHQU85ZxrDV5fDSwAvgPc6ZybGMyfBLwBjHPO1YcVb096Og7n3Mlm9iPgSeCv7csl+cSaB8zsl8AO59wlwetLgdnOufNCCLtbcRzL0c65h4LnI/FDih3unHsuhLC7FW9+NrOf4HPFB51zS4cy1r7EcyxmtgK40jn38NBH2rc4PmNfBQ5yzp0VvP4ZsMo5l3R3Xenu72nU8qT+3U/HriBLgUrnXANA8Ad2PbAkzKD6I/ga5yzgWuAdMzs75JDi4pxri5q1lBQ9N9HHkmrnxjn3eFSxuQXYCByNPwft620JnpYPYXgx6+k4zCwfOBX4U/D6mFAClFgsJbY80OWziS/ilg5BfPFYSgzH0l5UB893ATXA5qELMyZLiTE/m9kngVuBhqEMMA5LieFYzOxUYA6w2Mz+ZmZXm1nOUAfbh6XEdl7uA95nZmeYWSl+TOjfD2WgseqmNoiW1L/76VhYTwKqoubV4r+2SinOOwz/C3ID8AczOynksAZC5yZ5HIzvxpLq5+Rg/DcGDc652fjjuR/4u5ntF25o0oNYP3PR6yXj5zLu3x8z2wfffWLbYAbWDzEdi5nNAkY5554ZqsD6Idbzcg7wFHA9cC6+seR7gx5dfGI6Fufca/juRr/F/9NznnOucUgiTLyk/t1Px8LaAdEfphFASwixJIRzbqdz7jLg+8CX+lo/iencJAEzmw1UOOdWkcLnJOo4AHDObXPOXQD8GfhcaMFJb2L9zEWvl4yfy/78/nwW+OqgRdR/fR5L8BX+Z4EfDWFc/RHreVkAPOmca3HOVQE3Ah8egvjiEc9nLBt4PzAbuD3ZL0TvRVL/7qdjYb0FfzFCpCJgawixJNr1+Au0UpXOTcjMLBP4FPCNYFaXcxJ0cSkgyc9JN8cRLSXOR5qKNQ9Er5eMuSKunGZmx+Jbq9cNdmD9EMuxHAEsA942s0r8t3VHBM+TSaznJSuY2q0GRg1iXP0R07GY2RLgYOfcCnw3kQOAi4YkwsRL6t/9dCysHwYmmlkeQPAf21Tg0VCjSow2/JWyqUrnJnwXA9c453YHr1cAcyOWTweagZVDHFe8oo8jWqqcj3QUax6I/mzuE7w3mcSc08xsPn7kibuGNsSYxXIszwLz8EXbAcDl+FxxwNCG2qdYz8tqun7GWoHXhiTC2MV6LGcC6wCcc1uBy4CjhjDORErq3/20K6ydc2/jO/G3X7x0LHBHxEVZKcPMJprZOWaWEbQkfgW4NOy4YhUxhqZBap+b6GNJxXNjZpcALwJ5ZjbTzD4BVAPbzGxhsNqJwA3Oueaw4uxLd8dhZgvN7H3B8hzgY8DVYcYp3estD5jZxWbW/gf1l8AJ0PH7dzT+m4ikEeuxmNlMfEvvP81supktMrNPhxN192I5Fufcbufc5vYJ2AnsDp4njTg+Yz8BTjKzEcHrw0nRzxjwEnBgxFv34P8RSjrRf0+DeSnzu592w+0BBFfEXg1sAMYD33TO1YUaVD+Y2QLg7/gryJ8Afuycez3cqGJjZkXAR4CfA1cAP3XOVaTiuenuWICxpNC5MbPL8EPrRXrdOTffzKbhW542AIXApck6XF1Px4HvT/1H/JXkTwNXBa02koR6ygNm9jx+HPW/BOt9BFgcvO1R51zS3fSnr2PBDwH5DDAl6q1nOeduH8pY+xLreYlY/+PAx5NtuD2I+zN2FH7kiTznXLJdvBjTsQT937+NH6mlEv95u7KXb/VC0UttkDK/+2lZWIuIiIiIJFradQURERERERkMKqxFRERERBJAhbWIiIiISAKosBYRERERSQAV1iIiIiIiCaDCWkREREQkAVRYi4iIiIgkgAprEUkqZjbGzO4xs7fM7BthxyMiIr1T3u6kG8SISFIxsxOAh4Fi/F0T93HO7Qg3KhER6YnydicV1iKStMzsMeDoZL2FuoiIdJXueVtdQUQkKZnZFOCf6ZqcRURSjfK2CmsZBGb2NTP7WdhxtEtEPGZ2rJk9ZWZLExTWsGJm1yR4ewZ8EPhBIrcrIl0pXw+cmU03sxvN7FcJ2l6mmV1gZm/FsO7JZvZMf481kblbedtTYZ1izGy+md1pZr82s1Vm5szsI2HHFeVvwB/DDiJCIuLZBCxKQCxJyczea2aVZrbNzN5lZlOCz5kzs5+aWUmw3iFm9oqZXWtmecG8ccDbCQ7pHOAm51xrsH2RlKN83S+pmK+3B48jErS9DKAamBrDuiuB+T0tNLMFZnZ78Nm7zMzGRyxLdO5W3kaFdSr6K/C/zrnzgQOB37UvMLP9zWxJGEFF7ts594pz7okw4uhOPPH09DN0zr0BVCYyrkSer4Fuyzn3N+DnQKtz7inn3CbgAqAZeMY5Vx2s9yw+kX/VOdcYvP1U4L4BHUAEM/sW8G3gQTN7HZibqG2LDDHl6zilYr52zjUAWxO1L+dcC/B8jOu+gy/Ce1r+KrA8eHmNcy6ykE5Y7lbe7qTCOoWY2ShgH6AewDnXBlwOODMbiU/aFkJcoe07kWI4jrYh3FcY27odmGxmBwE456qAfwGnROyrBNgafPbazXPOvT7AfXdwzn3HOTfHOVfunJvnnHssUdsWGSrK14MrCfN1okeCiCf+vtbdEzxG93tOWO5W3u6kwjqFBEPXvAE8YGaHBfPeAp7F92uaASwzs/MAzOwKM/tW8DVQRz8qMys3s+vM7Etm1mBmb5jZRWb2fjN7zsxON7O7zKzGzD4V8b5utxe17wvM7DeRfc2Cr0NvMLPvm9nDZnZIML/X/UW838zsU2b2qpkdZ2arzewdM/t4DPuYExlPH/vc62fYjdLgq91aM7sg6mf67eBnc2eQiGP9me21LzM7w8yeNbOPm9mbZvaDXrbX3bnvNp7eBC0brwJnRsweA5xsQbcP4H3APRFx5gENEa/nmdkPzOwaM3vUfJeSD5nZSjM7zcz+aGbVZvZFMzsoOMYqM3tvX/GJpBLl6/TJ193p6Rh724+ZzTCz/zPf+vvDWPbTzX7/ZmZ/MLMj+lhPuXuwOOc0pdAELADewv+H+n/A2IhlG4ClwfOFwLrg+Rj8f9Nl+P+4NwLlwbKfAf8OnucAO4Cf4seiPBvY0tv2oveN72P2K+A3wfxcfLFWErw+A6gAinrbXzfHPT7Y5yXBut8FmoBpfewjOp5e9xn5M+wmhg3AtUB+8L6Nwfxi4LaI9Z7EJ8WYfmY97GsU/uu9HwOHAMf2tr2oc99tPDF+vi4D/hs8nwV8H9gNfCDi82IR678POCp4ng88B+QHrx8Lfg7FwC7gCiAP+FDw+oP4f+6/QvAZ1KRpOE0oX6dFvg6WfzvGn2Nv5/o5YH7Ee1yMn7P28/ku4P91s3xpsJ+siHnK3YM0qcU6xTjfqrgfPtF8FHjVzLq7SONN4FwzywaOCuYV4n+RpwDt/WNfw//C4JxrBmqBvzjnavC/aOP72F50fLvpejHEyUCjC/roAnfifyE/0Mf+orfbvs1bg3W/i+//e1wf++gSTzz77MG9zvenew6YEHGMZUGL0pfwP1OI8WfWw/HuwCewu51zzzrnHoxjez3FE4vbgZnmu4Zg+R4AAAT4SURBVIP8D3Aj8A/gg2ZWCNS6IPMGjgD+HTx/H/6PV3sryOnAd4KfczXwkPP9slcCxc65Pzn/9fiLxHcORFKC8nV65Otu9HiMveznWGCkc649npVx7vNg4MPOuf+NcX3l7kGSFXYAEjszywSmOOc2AF8ws5vwV1AvBw6NXNc5t9vMJgFfB37bvgnn3HYzWwkcg/+Pega+mOp4a8TzVoLuQj1tr4dQI7cxB8iOiKvNzDYAk3vbX1+cc01mtg4oAUbHsY9+7zNKK52/P1OB9c65H0evFMvPzMzOxRew7eY6/5Wxi4w1jnPQYzx9cc69aWYv4ruDjHTOrTOz2/GtSGcC90fE3f7ZaO+/Nw3fwtS+rarITUc830NXbahbmgwzytddji8d8nWkHn+OvZyb/fH/RPTX0cAcMyt0ztX1sE5bcAzK3YNIP5DUkk1E/1fn3Cr8V23zolc0s6OBzzjnvuv8CA+RPgscbH7Yp9XAVX3tuI/t9WYjvgU0chgiw/c9HKgcYO0g7yMW2+jaDxkzOzSOn9k9wAERU7dXl8exvW7jif1wuA04D1gXEV8m8EX816btDsX3F223BTjSzAoi9ttrPz+RYUz5uqt0ytc9HmMv+6kJ3pND/1yDb+3/SQ/Ld7nOi86VuweRCuvU83kzmxjxehLwQPC8GX+xxjz80E4jzWyEmR0bLB9rZqPxff1uwn+N82+CrxYDGUT9l25m1sf2ovdtEdv4K1AHnBVsqyTYX3vLZ0/768mEYJ1JwXsfiGEfkfH0tc/I44jWXeuF4VuhCoD7zOwE8xfJ9HUOOvYFTHDOvR4xtUTsIjPieW/bi4y723jMLMfMlgfdPHrzJ/xX0H8BcM7VAn8HnnRdRwM5Ht9NpF37z/uP5sfC/hzQPpZp9DmIPs8pPUKBSA+Ur0mbfB3rz7Hb/eBzaQ7+ny/w305gZhNjyN2Gbz0+FzjTzD4YvPcrZnYZvrX+XxHrK3cPpsHouK1pcCb8BREO2AncAvwG/1VSabD8/wGbgdPwX+28iW8hOBnf0vEbfJK6B38hRXOwvQb8L9pJQAv+ApkpwDeD5Z/obXtR+/4i8Di+39qCYFk5vqXzWnxfw4XB/B7318Pxu+D9X8R3TdgvYllP+5geGU9f+4z8GUbt+zT8RXw/DX4W7e87L1i+BHglODfXBj/nWH9mp3VzrGcE+/sjMCmY19s57bKtHuKZjL+Rwa0xfNauj3r9IeA9UfOu6eZ97wbWAFXAZcG89+I/azcAEwmGHAMuxP/hvTE41lPC/h3TpClRE8rX6ZSvpwJP4FvdD+jjGHs71yfjvyl8Gn8h+VrgM8Gxd5u78d+KtOBbqkcBD+GL+vOBT+HH874W3w2l/T3K3YM4WfDDkzRhZmOBC51zVwSvM/CtKJ9yzl0WanB9MDMHzHC+z6L0k5l91TnXr6GcRGToKF9LJOXu1KCuIOnnc/huAaMi5h1NxEVpySjiqyd97dRPZpZpZh8G7gg7FhGJifK1KHenGBXW6eeX+L5Yr5nZJmAFsNk591S4YfUsuJjji8HLj5sf9k3iVwj81Tm3rs81RSQZKF8LKHenFHUFERERERFJALVYi4iIiIgkgAprEREREZEEUGEtIiIiIpIAKqxFRERERBJAhbWIiIiISAKosBYRERERSQAV1iIiIiIiCaDCWkREREQkAVRYi4iIiIgkwP8HXA4QO52ELGcAAAAASUVORK5CYII=\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-1-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/viking-1-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/viking-1-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 }