{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 28 - Pathfinder - 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 Pathfinder mission.\n", "https://en.wikipedia.org/wiki/Mars_Pathfinder" ] }, { "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('Pathfinder', 586 , 62.0, 0.00, 5.52, 0.0, 0.66, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,7.48,0.0,-14.06,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,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZyVdfn/8dc1MwwgIOuIgMKwo+CCIpJ7ikuaRqa5/szMpqxss776VVMr/VZuWWrWZGa2mmlu5VqpuAeCCsQmw9KgOCCb7DDX74/PfYYzw+xzlvuc834+HvfjnHOf+9z3NbfjNRefc92f29wdERERERHpmKJsByAiIiIikg9UWIuIiIiIpIAKaxERERGRFFBhLSIiIiKSAiqsRURERERSQIW1iIiIiEgKlGQ7gFTp16+fl5eXZzuMWKupqaGsrCzbYeQknbv207lr2fTp01e6e0GdpPbmbP0+tY/OW/vp3LVfvp675nJ23hTW5eXlTJs2LdthxNqECRN0jtpJ5679dO5aZmZLsh1DprU3Z+v3qX103tpP56798vXcNZez1QoiIiIiIpICKqxFRERERFJAhXUBqaioyHYIOUvnrv107iSV9PvUPjpv7adz136FeO7M3TN7QLPjgeuBs9x9cbTudOAWYHfgd8Bl7r49eq8fcAOwBugEXO7u2xrud8KECZ6PfTwikv/MbLq7T8h2HI1RzhYRqa+5nJ3REWsz6w90ByYmrRsMTAHOAL4BXAx8PeljDwD3uPvlwHzgexkLWESkgClni4i0TUYLa3dfATzSYPUQ4CJ3n+7u9wF3Ah8FMLNJwGh3fy3a9jHgUjPrlqmYRURa4u7s2LGDHTt2ZDuUlIpTznaHHTvCIiISVxnvsXb32gavpya+QoxUA0uj58cCVUnbVkdPY/mVqYgUpvvuu4+SkhIuuuiibIeScnHJ2S+/DCUlcNRRHd2TiEj6xPHixUOAn0fPBwGrGry/HhjY8EM1NTVMmDChbqmsrExzmCIiQVFRSKW1tbUtbLlTZWVlXb4C+qUnsozIaM7O8GVBIiJA63N2rG4QY2YjgBp3fzNa5cCmBpt1Bna5EKasrCwvJyEXkfhLFNZtaQWpqKiou2LezFamJbA0y2TONosOoMJaRLKgtTk7NiPWZlYMfAG4Iml1NdA7aRsDugHLMxudiEjTiouLgbaNWOe6TOdsFdYikgtiU1gD3wJudvctSeueAUYnvS4HtgIamhaR2GjPiHUeyGjOThTWIiJxlvHCOhrBALCkdVcCM4CuZjbMzC4ysxHuPg1418zGRZueBNzu7lszG7WISNPyecQ6bjlbI9YiEmcZ7bE2sx7A+dHLC83sDsJXid9vsOlcd78nen4WcI2ZLSbMp3pVJmIVEWmtfB2xjlPOViuIiOSCjBbW7r4euCtaEq6PlqY+swT4XJpDExFpt3wdsY5TzlZhLSK5IE491iIiOak90+1J26jHWkRygQprEZEOytdWkDjSiLWIxJkKaxGRDsrXVpA4USuIiOQCFdYiIh2kEev0U2EtIrlAhbWISAdpxDr91GMtIrlAhbWISAdpxDpzNGItInGmwlpEpIM0Yp1+agURkVygwlpEpIM03V76qbAWkVygwlpEpIPUCpJ+6rEWkVygwlpEpIPUCpI5GrEWkThTYS0i0kEasU4/tYKISC5QYS0i0kEasU4/FdYikgtUWIuIdJBGrNNPPdYikgtUWIuIdJBGrDNHI9YiEmcqrEVEOkjT7aWfWkFEJBeosBYR6SC1gqSfCmsRyQUqrEVEOkitIOmnHmsRyQUqrEVEOkgj1pmjEWsRiTMV1iIiHaQR6/RTK4iI5AIV1iIiHaQR6/RTYS0iuUCFtYhIB2nEOv3UYy0iuaAk0wc0s+OB64Gz3H1xtK4rcCOwHugNXO7u61p6T0QkDvJ5ur245WyNWItInGV0xNrM+gPdgYkN3voZ8KK7Xwk8Avyile+JiGRdvraCxClnqxVERHJBRgtrd19BSLR1zGwgcC7wWLTqWeA0Mytv7r2MBCwi0gr52goSp5ytwlpEckHGW0HcvdbqN8sdA6x0943R+9vNrAo4GtjWzHuLMxm3iEhT8nXEGuKTs9VjLSK5IOOFdSMGAasarFsPDAS2N/OeiEgs5OuIdROymrM1Yi0icRaHwtqBTQ3WdSaMfNDCe3VqamqYMGFC3euKigoqKipSGKaISOPaM2JdWVlJZWVl4mW/1EeVNlnJ2WoFEZFsam3OjkNhXU24cjxZD2A5IYE39V49ZWVlTJs2LS0Biog0pz0j1smFpJmtTEtg6ZGVnK3CWkSyqbU5Ow7zWP8LGBhN0YSZdQIGA8+38J6ISCzk83R7jchKzlaPtYjkgowX1rbzKhgDcPf3gMeB46L1k4EH3b26ufcyGLKISLPy+eLFuOVsjViLSJxltBXEzHoA50cvLzSzO9y9BrgEuNHMDgD2BJIb7Zp7T0Qk6/L14sU45Wy1gohILshoYe3u64G7oiV5/Wrg8018psn3RETiILkVxN2xPOlbiFPOVmEtIrkgDj3WIiI5zczqimlX5ZcWefJvFRHJcyqsRURSIJ/7rOMg6rZh+/bsxiEi0hwV1iIiKZCvfdZxUVoaHrftMiO2iEh8qLAWEUmBAptyL+MShfXWrdmNQ0SkOSqsRURSQK0g6dWpU3hUYS0icabCWkQkBdQKkl4asRaRXKDCWkQkBTRinV7FxVBUFKbb0ykWkbhSYS0ikgIasU4/jVqLSNypsBYRSQGNWKefCmsRiTsV1iIiKaAR6/RTYS0icafCWkQkBTTdXvppZhARiTsV1iIiKaBWkPRLjFhv2ZLdOEREmqLCWkQkBdQKkn7duoXHjRuzG4eISFNUWIuIpIBGrNOvR4/wuG5dduMQEWmKCmsRkRTQiHX67b57eFy/PrtxiIg0RYW1iEgKaMQ6/RKFtUasRSSuUlJYm9nJqdiPiEiuyqUR61zN2YlWEI1Yi0hclbRlYzM7Cbga2JOdRXkR0B/omtrQRERyRxyn28u3nK0RaxGJuzYV1sBvgR8AM4DEXw8Dzk1lUCIiuSamrSB5lbNVWItI3LW1sF7k7rc2XGlmc1MUj4hIToppK0he5ex+/cJjTU124xARaUpbe6wvMbMzGll/WiqCERHJVTEdsc6rnN2/f3hcsSK7cYiINKWthfXhwN1mtiNpqQXuSkNsIiI5I6Yj1nmVs/fcMzy+91524xARaUpbW0G+B1QA04DEsIwB/6+jgZjZPsBXgIXASKDS3WeaWVfgRmA90Bu43N3VYScisRLTEeu8ytmJEWsV1iISV20trJ8AHnH3LckrzeyXKYjlt8An3L3azAYDTwH7AD8DnnT3+6Mr3H8BnJOC44mIpExixDpmhXVe5ezEiLVaQUQkrtpaWP8J+F8z+2fSOgNOB77WwVj2AaJZStkE9DSzgYSr178crX8WeNDMyt19cQePJyKSMokWkMTIdUzkVc7u1QtKS8OsIBs3wm67dXSPIiKp1dbC+juEr/wuBDxaZ4Q5UjuapP9C6AX8GHA+cClwDLDS3TcCuPt2M6sCjgYWd/B4IiIpkyisEyPXMZFXOdsM9t4b3nkHliyBffbp6B5FRFKrrUMrNwF7u3u5uw+NlnLC6EdHfQnYDPwb+NDdHwQGAasabLceGNjwwzU1NUyYMKFuqaysTEFIIiKtk2gBae2IdWVlZV2+AvqlKay8y9lDh4bHqqr2By4i0latzdltHbFe5u6N3Ux2zzbupzFdgT8SEvNt0SiHE75iTNYZ2Nbww2VlZUybNi0FYYiItF1bR6wrKiqoqKgAwMxWpimsvMvZicJ60aI2f1REpN1am7PbOmL9pYYrzOxM4LY27qcxvwMedPfvAT8kJOx3CVeVJ+sBLE/B8UREUqatI9YZknc5e9iw8KgRaxGJo7b+BehjZt8BMLPdzOxu4G7CV4HtZmb9gAPcfU206npgHbAEGBhN34SZdQIGA8935HgiIqkW0x7rvMvZI0aEx7k5ee9IEcl3bSqs3f0U4Bkzux2YCYwCDgAau7NXW3wAbDazQUnrVgFvAY8Dx0XrJhNGSKo7eDwRkZSK44h1Pubs/fcPj2++mao9ioikTrM91mZ2VBNvLQEWADcDQ4ApwDfaG4S715rZFOAaM5sO9Ae+7e7rzOwS4EYzO4DQF1jR3uOIiKRLHEasCyFnDx8eptmrroZVq6Bv31TuXUSkY1q6ePFGwlyla4DG7tN7D2HqpgF0IEkDuPubwBcaWb8a+HxH9i0ikm4xGbHO+5xdXAz77QevvRZGrY89Nl1HEhFpu5YK6+8B09z9/eY2MrOOfq0oIpLT4jBiTYHk7IMPDoX1q6+qsBaReGl2aMXd/95Sgo62+0vqQhIRyT1xGLEulJx9VNTw8rwuYxeRmInPVTYiIjksJiPWBSFRWL/0Emzfnt1YRESSqbAWEUmBbdvCPVBKStp63y1pqwEDYORI2LABZszIdjQiIjupsBYRSYFNm8INB7t27ZrlSArD0UeHR7WDiEicdLiwNrPxZvY7M7vMzMo7HpKISO7JlcI6X3K2+qxFJI5SMWL9RWCyu98CHJa445aISKHYsWMH27Ztw8woLS3NdjgtyYucnZgN5Nln4cMPsxuLiEhCKgrrrwAHR8//BHRKwT5FRHJG8mi1mWU5mhblRc4eNAgOOww2b4a//S3b0YiIBG0urM3sY2Z2XvR8IjAycbtad69193UpjlFEJNbi3AaSzzn7zDPD4wMPZDcOEZGENhXWZnYD8ADwSQB3fx24wMzOTUNsIiI5YePGjQDstttuWY6kvnzP2aefHh6ffDKMXIuIZFtbR6wPA/YCpiatewL4YcoiEhHJMTEesc7rnD14MIwfH6bd+8c/sh2NiEjbC+tn3H0N4EnrcvrWuCIiHRXXEWsKIGefdlp4fPTR7MYhIgJtL6wXmdmVwGgzO9vMHgS+DPwo9aGJiOSG1atXA9CrV68sR7KLvM/Zn/xkeHzgAYj+fSMikjVtKqzd/U/AP4Fi4AJgDXCCu9+ZhthERHJCorDu3bt3liOprxBy9gEHwMSJsHo1/PGP2Y5GRApdm++96+6vAq8mXptZsZkd5+7qcBORghTXwhoKI2d/+cvw+utw551w0UUQ/xkPRSRfNVtYm9mjQPcW9tEPWAnkTZIWEWmLuBTWhZqzP/1puOwymDEDXnklzG8tIpINLY1YvxVtMw+oBT4PzACmJW0zFOiTluhERHJAXAprCjRnd+kCF18MP/wh/PjHKqxFJHtaKqxvAiy6qhwzG+7u1zTcyMweSUdwIiK5IEaFdcHm7EsvhVtugYcegnfegeHDsx2RiBSiZi9edPe1iQQd6W9mxcnbmNkxwMQ0xCYikhNWrVoFQJ8+2R0ILuScPXAgnHce1NaGUWsRkWxo63R7fwHeNLOfm9lNZvY48Azw09SHJiKSG959910ABg4cmOVIdlFQOfub3wyPv/kNrMvZG7WLSC5r63R7zwAnEC58GU2YuukMd/9BGmITEckJ1dXVQPwK60LL2fvtB0cfDR9+CHfdle1oRKQQtWe6veXA1YnXZlZkZge7+/SURiYikgPcneXLlwMwYMCALEezq0LL2VdeCc8/Dz/6EXzhCxC/e/aISD5rU2FtZvc0sroM2AFMSUVAZjYJOAJYCLwIbABuBNYDvYHL3V1f8olILHzwwQds3bqVnj170q1bt2yHU08h5uzjj4djjoHnnoObboIbbsjUkUVE2t5jPRJYCixJWjoDKRn5MLPPAh9395vd/WF3Xwn8DHjR3a8EHgF+kYpjiYikQmK0etCgQVmOpFEFl7PN4AdRo8ttt0HU/i4ikhFtbQW52N3nJa8wszHAhR0NxMyOBs4Djk9aNxA4F/hytOpZ4EEzK3f3xR09pohIRyX6q+PYBkKB5uxJk2DKFHj4Yfjf/4V7783UkUWk0LX14sV5jazeAlycglhuBf4D3G5mT5rZ4cAxwEp33xgdfztQBRydguOJiHTYwoULARg2bFiWI9lVIefsm26C0tIwQ8hTT2X66CJSqNpUWJtZlZktSlqWEBLrXzsShJmNBsYDv3T3rxButfsksDewqsHm64FdLr2vqalhwoQJdUtlZWVHQhIRaZUFCxYAMGrUqDZ9rrKysi5fEW4znnKFnLNHjICrrgrPP/UpmDo1pbsXkQLT2pzd1laQZ4A/AB69rgVWuPv8dkW501jgA3d/K3p9B3BtdJxNDbbtDGxruIOysjKmTZvWcLWISFrNnx/S38iRI9v0uYqKCioqKgAws5UpDywo6Jx99dWwcCH89rdwwgnwl7/AKaek9ZAikqdam7PbWlj/zt1fSF5hZqPNbIS7L2x7mPXiqIvF3TeZ2QKgmHBVebIewPIOHEtEJGXaO2KdIQWds4uK4Ne/hi5d4Je/DH3X994b7tAoIpIObZ0V5KRG1q0AHuhgHG8BPc2sf9K67YRkPNDMugKYWSdgMPB8B48nItJhGzdupKqqiuLi4lj2WKOcTXEx/OIXcMUVsH07nH8+fPe74dbnIiKp1qrC2sy+YmZrgMvNbEfyQuinW9GRINx9LvAEcEZ0vF5AKfB74HHguGjTycCD7l7dkeOJiKTCzJkzqa2tZezYsXTu3Dnb4dRRzq4vMQXfzTeH59ddB6eeCqtXZzMqEclHrWoFcfc7zOwF4DPATxu8vdXdUzFT6AXAT8yslDDCcY67bzezS4AbzewAYE+gIgXHEhHpsH//+98AHHLIIVmOpD7l7MZddhmMGwfnngt//zsccECYNeSjH812ZCKSL1rdY+3ub5nZFe6+y0UoZjbY3Zd2JJDoxgK7dL65+2rg8x3Zt4hIOiQK6+gq8VhRzm7ciSfC9Olw1lnw+utw3HHwzW+GOzTG6EsHEclRzRbWZvYt4O/uPida9T9mVtxwM+CjhPlLRUQKgrvzwgvhusBJkyZlOZpAObt1ysvhxRdDMX399XDLLfDXv4Z2kTPPDO0iIiLt0VKP9VCgV9Lrg4DDgGHRe4llz7REJyISU/PmzWPZsmWUlZWx//77ZzucBOXsVurUKfRav/gi7LMPLFoURrH33x/uuw+2bs12hCKSi5odsXb3LzdYdQMwy93rpRwz2zfVgYmIxNnTTz8NwPHHH09RUVsnWEoP5ey2mzQJ3noL7rknzBYyaxZ85jNw5ZXwta/BRRdB377ZjlJEckVb/xrs10iCHg0MT11IIiLx99BDDwFw4oknZjmSZilnt0JJCVRUQFVVmPd6332huhr+539g0CC44AJ45RVwb3lfIlLYWjvd3pFmdhRwvJkdlbwA+wB3pjVKEZEYWbZsGc8//zxdunThE5/4RLbD2YVydvuUlsKFF8Lbb8Pjj8PHPhZaQn77WzjssNAm8qMfwbJl2Y5UROKqtbOCrAHuI/TmHd7gva3A3akMSkQkzn73u98BcOqpp9KzZ88sR9Mo5ewOKCoKtz4/5ZTQe11ZCb/6VWgTueKKsBx1VLiD4xlnQJ8+2Y5YROLCvJXfbZlZT+BMd49lQp4wYYJPmzYt22GISJ7bunUrw4YNo7q6mieeeIKTTmrs5oZtY2bT3T2lc/YpZ6fW1q3w5JPw+9/Do4/C5s1hfUlJKLI/8YmwDBmS3ThFJP2ay9mt7rF297VNJWgzO6K9wYmI5JI///nPVFdXs++++8a6v1o5O7VKS+G00+D++2HFitCLPXly6Lv+5z/DhY7l5XDggXDNNfDyy+EW6iJSWFqax/pRoHsL+zBgLLBHqoISEYmjrVu3ct111wFw2WWXYTGb8Fg5OzN23z30Yl94IXzwQbiL4yOPhBHtN98My/e/Dz17hhvQnHBCuDFNeXmWAxeRtGupx3oeoVevGmiqZ8SAj6cyKBGROLrzzjt55513GDNmDBdccEG2w2mMcnaG9ekD558fls2b4V//CoX200/D/Pnw0ENhARg5MhTYxx0HRx6pafxE8lGzPdZm1g9Y7+5bmt2JWV93X5Xq4Noi1/r1RCS3LFu2jHHjxrFu3Toee+wxPv7x1NWmqeqxVs6Ol6qqUGA//TT84x+wdm3998eNg6OPDj3aRx0Fexb8bXtEckNzObulG8SsbOUxPomuMheRPOXuXHzxxaxbt47TTjuNU045JdshNUo5O16GDoUvfCEs27fD66/DU0/B88/Dq6+GWUZmzYI7o8kPR42Cww+HQw8Ny7hx4eJIEckdHfpf1swGA18GvoKStIjkqZtuuomnn36aPn368Itf/CJ2vdWtpZydPSUlYS7sww4Lr7dsCYX2Cy+E5aWXQuvI/PnhwkiA3XaDgw/eWWhPmgR77ZW9n0FEWtbmwtrMioHTgM8DxwMrCD19IiJ558knn+SKK64A4J577mHPHPu+Xjk7njp3Dn3WRx4JV10F27bBjBnhDo+vvRaWRYtg6tSwJPTvD+PH11+GDQtzb4tI9rW6sDazocDFwGeBvsAS4FTgKeAjaYlORCSLpk2bxllnnYW7c+2118byLotNUc7OLZ06wcSJYUmoqQmj2olC+/XXw1R/Tz4ZloQePcI0f4lCe7/9YJ99woi3iGRWi4W1mZ1JGOk4FlgI3Ar8BrjS3RP/a7+ctghFRLLgrbfe4sQTT2TdunWcddZZXHPNNdkOqVWUs/NHWdnOO0BCmDO7qiqMbM+YAW+8ER7fe2/XkW2zMJI9duzOZdw4GD0aunTJzs8jUghaM2K9iXAL3JeBi919PoCZte6WjSIiOWbq1KmcdtpprFmzhtNOO43f/va3FOXOd+3K2XkqUSwPGwaf+tTO9e+9FwrsmTPD4+zZoVf7nXfC8uijO7ctKoIRI+oX3PvuGy6cVMEt0nEtFtbu/jjweHTRy+fNrBx4BOiU2MbMzFt7b3QRkRj785//zAUXXMCWLVuYMmUKf/zjH+nUqVPLH4wJ5ezCs+ee8LGPhSVh69ZQXM+eHZZZs8LjwoU7L5L86193bl9UFAr2ffYJy777hscxY8INcUSkdVrdY+3uS4HvmFkJMAXYy8y+D/wKOAGoTE+IIiLpt3XrVi6//HJuu+02AL70pS/x05/+lOLi4ixH1j7K2YWttDS0fowbV3/95s0wd+7Ognv2bPjPf8LI9sKFYXnssfqfGTRoZ6GdXHSXlWXu5xHJFc3eIKbFD5uNAC4FPu/uWb1MohBuNiAi6VFVVcV5553HK6+8QklJCTfffDNf/epXMzatXqpuENOK4yhnS6M2b4YFC0KRnVjmzAkj21uauN1Q375hRHvEiHBXyeRHjXJLPmv3DWJa4u4Lga+Z2fSO7EdEJBtqa2u56667uPzyy9mwYQN77bUXf/7zn/nIR/Jz0gzlbGlKly5hNpH99qu/fseOcMFkcrGdeL5qVZh/+6WXdt1fWdmuxXbisWfPzPxMItmQkns6uft9qdiPmZUC/wa+5u7PmVlX4EZgPdAbuNzd16XiWCJS2GbOnMlXv/pVpkZTKZx11lnccccd9OvXL8uRpV+qcjYob+e74uJQDI8YAaeeunO9OyxfHka0FywILSSJx4ULw1SBNTXwciPzz/TrF4rsYcOgvDzcoTKx7L237jYpuS1uv77fBsqTXv8MeNLd7zezk4BfAOdkIzARyQ/vv/8+V199NXfffTfuzh577MFdd93F6aefnu3QcpXydgEyC73XgwbBRz9a/73a2lB0Nyy2E89XrgzLK6/sut/i4nB3yeRiO7n4HjBAN8OReItNYW1mhwHvAquj1wOBcwm33wV4FnjQzMrdfXFWghSRnLV69Wp+/OMfc9ttt7F+/XpKSkq49NJLueaaa+jVq1e2w8tJytvSmKKiUBzvtRccc0z99xIj3QsXhhaT5GXxYqiuhiVLwvLcc7vuu7QUhgwJRfaQITB48M7HwYNDoV9amoEfUqQJsSiszawbcIa7f9PMEndhOAZY6e4bAdx9u5lVAUcDi7MSqIjknNWrV3Pbbbdx2223sW5d6Eg4+eSTufXWWxk9enSWo8tdytvSHskj3Ucfvev7W7bA0qU7C+2Ghff774eR7wULmt7/gAG7FtyJZcgQ6NUrbCeSDrEorIHLgR82WDcIWNVg3XpgYEYiEpGctmDBAn7yk59w7733smHDBgAmT57MtddeyxFHHJHl6PKC8rakXOfOof965MjG39+wYWfBvXTprkt1dRgRX74cXn218X107954wZ086p1DU9dLzGS9sDazE4Fp7v5+g7eccAexZJ2BbY3tp6amhgkTds58UlFRQUVFRSpDFZGYq62t5dlnn+X222/nb3/7G4npRE844QS+853vxKqgrqyspLKybirpnLpiMhV5Wzlb2qNbt513jGzMtm2hqG5YcC9ZsvPxww/D7CZz5jS+DzMYOLD5Ue+ePTXqXWham7M7NI91KpjZM8D4pFW9gQ+BW4Hz3H1U0rYLgGvd/Q8N96M5UUUK19KlS/n1r3/Nr3/9a5YsWQJA586dOf/88/n617/OuIZ3yYiZTM1jnSqpyNvK2ZIN7rB27a4Fd/KyfHm4ALM5PXo0P+o9cKBGvfNZ2uaxTpFzCSMaCa8A3wReA75tZl3dfZOZdQIGA89nIUYRiZkPP/yQxx57jN/85jc8/fTTdaPTQ4YM4eKLL6aiooI99tgjy1HmLeVtyUlmoce6Vy/Yf//Gt9m2LbSUNDfqvX79zjtXNqaoKBTXjY14J496S/7JemHt7jXJr81sB1Dj7kvN7HHgOOBxYDLwoLtXZyFMEYmBTZs28cQTT/CnP/2Jxx9/nE2bQtdBaWkpp59+Op/73Oc49thjKdJ8XGmlvC35rFOnMMVfeXnj77vDmjUtj3r/979haewGOhDuTtnSqLfm9M49cf9Pdglwo5kdAOwJqAFPpMCsXbuWp556ikcffZRHH32U9evX17132GGHcfbZZ3PuuefSt2/fLEYpSZS3Ja+ZQe/eYTnggMa32bq15VHvdetg1qywNKaoKFxI2Vyvt24dHz+xK6zdvTzp+Wrg89mLRkSyoaqqiscee4zHHnuM559/nm3bdl77dvDBB3P22Wfz6U9/msGDB2cxSklQ3hapr7R0501tGuMOq1c3P+r97ruwbFlYmhr17tmz8VYTjXpnj063iGTdxqtG6MIAACAASURBVI0bmTp1Ks888wxPPfUUs5KGcIqKijjyyCM59dRTmTJlCiObmodLRCRHmEGfPmE58MDGt9m6NbSSNDXqvXRpuBDz7bfD0phEr/fee4dCu7HHfv00w0kqqbAWkYzbsWMHb7zxBs8++yzPPPMML730Elu3bq17v0ePHpx00kmceuqpnHzyyWrzEJGCU1oKw4aFpTHu8MEHjRfcy5btHPVO9Ho3dgt5gC5ddhbaTRXf3bql7+fMNyqsRSTttm7dyvTp03nxxReZOnUqL774IqtXr65738w4+OCDOf7445k8eTJHHnkkpbovsYhIk8ygb9+wjB/f+DZbt+6c1ztRbDd8XLOm+btZQhhZb2rUWy0n9ek0iEjKffjhh7zyyitMnTqVqVOn8tprr9XN4JEwdOjQukL62GOP1ai0iEiKlZY2P8MJhKkDmyq6Ez3eH3wQljffbHwfajnZSYW1iHTI9u3bmT17Nq+//nrdMmvWLGob3GFhzJgxHHnkkXVLeXOZXkREMqJHD9h337A0xh1qanZtM0l+VMvJTiqsRaTV3J0lS5bUK6KnT5/Oxo0b621XXFzMhAkT6oroI444grKysixFLSIi7WUGe+wRlglN3B+2sZaThsV3KlpOBg2C4uL0/JyposJaRBq1bds25s6dy8yZM5kxYwYzZ85k5syZ9XqjE4YNG8bEiRPrlvHjx7PbbrtlIWoREcm09racNCy+W2o56dwZRoyAkSNh1Kidj6NGQf/+8Wg1UWEtIqxdu5ZZs2bVK6BnzZrFli1bdtm2X79+HHLIIRx66KFMnDiRQw45hH79+mUhahERyRUttZzU1oaWk6b6vRcvhhUrmr6VfPfu9Yvt0aNh7FgYMya0oWSKCmuRArJhwwbmzJnD7NmzmT17NrNmzWL27NksW7as0e2HDx/OgQceyIEHHsj48eM58MADGThwIBaHYQEREckbRUVh1Ll//6ZbTtav39lOMn9+WBLPV6+GN94IS8P9jhwJ48aF5SMfgSOOSF8/twprkTy0adMm5s2bV694njVrFosXL8bdd9m+c+fO7LvvvvWK6P3335+ePXtmIXoREZFd9egBBx0UloZWrapfaP/nP2Fke8ECmDcvLA8+GLbt1AlOOAGuu67pIr69VFiL5Ch357333mPevHnMnTuXuXPn1j1fsmRJowV0p06dGD16NGPHjmXcuHGMHTuWsWPHMnz4cIrjfkWIiIhIE/r2DaPRH/lI/fWbN8PcuTBrFsycCc8/H0a1//Y3+Mc/YPr0pttT2kOFtUjMbdmyhYULF+5SPM+bN49169Y1+pni4mJGjBhRVzwnHkeOHEmnTp0y/BOIiIhkR5cu4bbxBx4I558f1tXUwOTJ8NZb8OyzKqxF8s727dtZsmQJCxcuZOHChSxYsIAFCxYwb948qqqqdpkTOqF3796MGTOGMWPGMHr06Lrnw4YNUwEtIiLSwJtvwh13wNtvh9cTJ6Z2/yqsRTJk27ZtLFmyhAULFtQroBcuXEhVVRXbt29v9HNFRUWMGDGiXuGceN6vXz9dSCgiItKE9eth6lR45hl48snQFgJhPuwf/AAmTUrt8VRYi6TQtm3bqKqq2qVwThTPO3bsaPKze+21FyNGjGDkyJGMGDGirpgeMWIEnTt3zuBPISIikns2bgwXLU6bBq+/Dq+9BnPmhLtHJvTpA2efDV/9apiSL9VUWIu0gbtTU1PDokWLqKqqYtGiRfWeL1u2rMm2DTNj8ODBdUVzcgE9fPhwunbtmuGfRkREJLe4h/msFy2ChQtD4TxnTpgBpKqqfhENUFISZhE5/vgwE8ikSeGGNumiwlqkgU2bNlFVVVWvcE4unjds2NDkZ82MIUOG1CuaE8+HDh2q4llERKQZiRvF/Pe/YVmyJBTRycumTY1/tqQkzFk9fjwcemjonz7wQN0gRiStamtrWb58+S6jzonX7777brOf79WrF8OGDatbhg4dWvd88ODBlKbzn8IiIiI56sMPw2jze+9BdXVYEgV04vny5bBtW/P76dsXhg2DoUNhn33CHRbHjg23O8/2n2AV1pKXVq9eXVc0J0afE68XL17M1q1bm/xsp06dGDJkSKPF89ChQ+ndu3cGfxIREZH42rgxFMorVuwsmpMfk58384VvPX37wl57hWXvvWH48FBIJ4rpON+7TIW15KTNmzezePHiRgvnqqoq1q5d2+zn+/fvv8toc+L5oEGDdLMUEREpOO6wZg2sXBmWmprGHxPP338/jEK3VpcuO29bPmhQKJwTj4ll4EDI5a5JFdYSSzt27KC6urrJUeeW2jW6d+/O0KFD65ZE4ZxYunXrlqGfREREJPPcw1RzH3wQllWrmi6YE89XrYImZn5tUmkp7LlnKJYTj0093313yPcZYlVYS1a4O6tWrWqycF66dCnbmmmyKikpYciQIU0WzprfWURE8kFjBXLiecPXDd9ra5EMofjt1w/KysJj8vOG68rKQluG/tzuFJvC2sxOB24Bdgd+B1zm7tvNrB9wA7AG6ARc7u4ttLVLHGzYsIHFixfvUjwnnn/YwvdHAwYMaHLUedCgQZSUxObXV6TgKGeLtE2mC2SA7t3DvM2JpaWCuW9f0G0TOiYWlYmZDQamAGcAY4G7gGXAzcADwBXu/pqZfRH4HvC/2YpVdtq+fTvLli1rctT5/fffb/bzu++++y4jzYnX5eXlmppOJKaUs6WQ1dbC2rX1C9/mlkSRvHp1+wvkbt1CYdy3784iuaXnvXurSM6GWBTWwBDgInffDkw3s3HAR83sRWC0u78WbfcYMM/Mrnf3Vl5bKu3l7rz//vu7jDQn3wyluTsJlpaW1s2u0dioc+/evdWuIZKblLMl523fHi7Ua22BnCiS16wJxXV7dOsWCt6+fVtfJKtAzi2xKKzdfWqDVdVAD+BYoCppu+qoEJsAPJ+xAPPY+vXrGy2cE8vGjRub/KyZMWjQoCZHnQcOHEhRUVEGfxoRyQTlbImTbdvaVhwnCuQWJo9q1u6712+xaFgQN7aoQC4MsSisG3EIcBNQAaxq8N56YGDDD9TU1DBhwoS61xUVFVRUVKQzxpyxdu1a5s+fz7x585g/fz7z58+vuyHKqlUNT299vXv3bvTiwGHDhjFkyBA6K0uItEtlZSWVlZWJl/2yGUsKKGdLh23e3PYC+YMPQt9ye5hBr16tL4wTS69e0KlTan92ib/W5uzYFdZmNgKocfc3zcyBhjeu7AzsciFMWVkZ06ZNy0SIsbRlyxYWLVpUr4BOPDbX69ylSxfKy8t3KZ4Tz3vGeRZ2kRyWXEia2cosh9NuytmSzD3cbrqxC/BaWpr5grRZRUVtGzlOLD17gm5ZIK3V2pwdq8LazIqBLwBXRKuqgVFJ7xvQDVie+ejiYc2aNcyePZvZs2czZ86cuhHoqqoqapto+uratSujRo2qtwwfPpxhw4bRv39/tWuISLsoZ+e/LVvC/MaJm4EkLzU1jRfIW7a071idOrWtME4sPXqE4lokDmJVWAPfAm5298T/ls8AX0x6vxzYCuT9MMf69euZM2dOXRE9a9YsZs+eTXV1daPbFxUVMXz48HrF8+jRoxk1ahSDBg1S8Swi6aCcnWN27AjFb8MiOblYTn7dnj7kLl3aXhz36RMu7NP17JLrYlNYm9mVhOTb1cyGAccALwDvmtk4d58FnATc7u5bsxdpark7S5Ys4Y033uCNN95g5syZzJo1iyVLljS6fdeuXdlnn30YO3YsY8eOZcyYMYwaNYphw4ap31lEMqZQc3Zcbd8OK1bAu+/C8uU7Hxs+r6lp24wWxcWwxx5hKSvb+TzxurECWjOlSiGLRWFtZlcD32+weq6732NmZwHXmNlioDtwVabjS5Xa2loWLFhQV0S/8cYbzJgxg9WrV++ybWlpKWPGjKkroMeNG8fYsWMZOnQoxWoKE5EsKpScHRebN8OyZbBkCSxdGpZEoZwonFesCP3NrdGnT/0CuWGxnPy6Vy+1WYi0RSwKa3e/Hri+ifeWAJ/LbESpsWbNGl599VVefvllXnnlFV577TXWN3L5cllZGQcddBAHHXQQ48ePZ7/99mPEiBG6s6CIxFK+5uxsWbcO3nknFMxLluwsoBOPK1a0vA+zUAgPHBiWAQN2fT5gAPTvrxktRNJJlVsK/fe//+Wf//wnL774Ii+//DJz5szBGwwh7L333nUFdKKYHjhwoG6UIiKSx9auhYULYcGC+o8LF4Ze5uaUlMBee8GQITB4cFgGDapfOKtgFokHFdYdsHbtWp577jmeffZZnn32WebOnVvv/dLSUg4++GAOO+wwDjvsMD7ykY8wYMCALEUrIiLpVlMDs2bVX+bPh5XNTKjYpQsMGwbl5TuL5+THAQM0LZxIrlBh3UaLFy/m4Ycf5uGHH+bFF1+sd0vv7t27c8wxx3D00Udz+OGHc9BBB+mCQhGRPFRbG0adX38dpk2Dt98ORXRNTePbd+0KI0bsXEaO3Pk4cKD6mEXyhQrrVli0aBG///3veeihh5g5c2bd+pKSEo444ggmT57M5MmTmThxIp30XZyISN5Ztw5eeAFeeSUU0//+d+NT0fXoAePG7VzGjoUxY8Kos4pnkfynwroJ69ev5y9/+Qv33nsvL7zwQt367t27c/LJJzNlyhROPvlk3ZlQRCQPbdkCU6fCP/8ZlmnTwhzQyQYMgEMPhUMOgQMPDIX03ntrLmaRQqbCuoHFixfz05/+lLvvvrtuBo+uXbvyqU99inPOOYfjjjtO7R0iInlo/Xr4+9/hr38Nj8mTOJWUwGGHwZFHhmJ64sRwAaGISDIV1pE5c+Zw3XXX8eCDD9bdGvyII47gs5/9LGeccQa77757liMUEZFUq62F556Du++Ghx6qfzvu/feHE06AY4+FI44IbR4iIs0p+MJ62bJlXHvttfzmN7+htraWkpISzjvvPL7xjW8wfvz4bIcnIiJpsGkT3HMP3HorLFoU1pmFAvqTn4QpU8JMHSIibVGwhfW2bdu49dZb+e53v8umTZsoKSnhi1/8IldeeSWD9P2eiEhe2rwZfvYzuPHGnTde2XtvuOgi+Oxnw/R2IiLtVbCF9SWXXMKvfvUrAM4880z+7//+jxEjRmQ5KhERSacLL4T77w/Px4+Hq64Ko9OaJ1pEUqFgC+tvfvObvPzyy9x2222ccMIJ2Q5HREQy4Otfh7lz4YYb4OSTNYOHiKRWwRbW++67L7NmzaJIE4uKiBSMSZNgxgwV1CKSHgVdVaqoFhEpPCqqRSRdVFmKiIiIiKSACmsRERERkRRQYS0iIiIikgIqrAtIZWVltkPIWTp37adzJ6mk36f20XlrP5279ivEc6fCuoAU4i94qujctZ/OnaSSfp/aR+et/XTu2q8Qz50KaxERERGRFDB3z3YMKWFmNcCSbMcRc/2AldkOIkfp3LWfzl3Lhrh7WbaDyKQO5Gz9PrWPzlv76dy1X76euyZzdt4U1iIiIiIi2aRWEBERERGRFFBhLSIiIiKSAiqsC4CZnW5mVWa2ysx+YmYl2Y4p15hZqZm9aWbHZDuWXGNmk8zsW2Y2xcz6ZTseyU1m1tXMbjez/zOzu8xs92zHlEvM7Hoz82h5M9vxxJmZHW9mr5lZedI6/f61oLHzFq0vqN89FVh5zswGA1OAM4CxwF3AMuDmbMaVg74NlGc7iFxjZp8Fhrv71dmORXLez4An3f1+MzsJ+AVwTpZjyglm1o1wEdnx0Spd6N8EM+sPdAcmNnhLv3/NaOq8FeLvni5ezHNmdiTwirtvj17fCIx191OyG1nuMLPDgDHANcCF7v5cdiPKDWZ2NPAd4HhXopEOMLOBQBXQ2903Rt+6rSXkssVZDS4HmNnXo6c/d/fNWQ0mB5hZEbADGOrui/X71zoNz1u0ruB+99QKkufcfWqiqI5UA0uzFU+uif61fYa735PtWHLQrcB/gNvN7EkzOzzbAUnOOgZY6e4bAaKcVgUcnc2gcoGZGXA2cAuwwsw0ytoCd69tsOoY9PvXoobnrVB/91RYF55DgJ9nO4gccjnww2wHkWvMbDQwHvilu38F+AfwZPR1oUhbDQJWNVi3HhiYhVhyigeTCF/H3w783sw+luWwco1+/9qhUH/3VFgXEDMbAdS4e95fPJAKZnYiMM3d3892LDloLPCBu78Vvb4DMOAT2QtJcpgDmxqs6wxsy0IsOcndV0fXOvwA+HpL20s9+v3rgEL73VNhXSDMrBj4AnBFtmPJId8C7jGzlWa2EtgbeMTM/ifLceWCEpIujnb3TcACoG/WIpJcVg30brCuB7A8C7Hkup8Cg7MdRI7R719qFMTvngrrwvEt4GZ335LtQHLIucCBScty4GLUStMabwE9G7R+bAfmZCkeyW3/AgaaWVcAM+tE+AP9fFajyk21wBvZDiLH6PcvNQrid0+FdQEwsyuBGUBXMxtmZhdFbSHSDHevcff/JhbC1c417r4u27HFnbvPBZ4gTPOImfUCSoG/ZTMuyU3u/h7wOHBctGoy8KC7V2cvqtxgZgPN7FwzK4ouJrsMuCrbccVZdJ4gtK/p96+VGp63Qv3d0zzWec7Mrga+32D1XM1yIRlwAfATMysljO6c02CGGpG2uAS40cwOAPYEKrIcT67oTbgA+0rgReA2TRHXNDPrAZwfvbzQzO5w9xr0+9esxs4bBfq7p3msRURERERSQK0gIiIiIiIpoMJaRERERCQFVFiLiIiIiKSACmsRERERkRRQYS0iIiIikgIqrEVEREREUkCFtUgTzOw4M3vZzI5pYbtxZvZnM/tOCo99U6r2JSKSDWa2v5k9kMrc2MRxHjezs9N5jHQxs8lm9kpLf2ckd+gGMRJ7ZrYPcAPwAXAIsD/w/9z9d2k+9DLggFZst5QwEX5xKg4a3S53Wyr2JSICYGaHAD8AjgKeBqoINzrpBNzl7k+l4bALgZ6kKDc242fArHTsOMrHpwM/B+YDr0RvdQY+BUzs4E1PWvt3RnKECmvJBQ8Dn3b3N82sCPh1Jg7q7vPNbFUrtltnZu+m8NBHAy+kcH8iUuDc/d9m9idgjLt/PLHezE4B/mpm33L3n6X4mBtTnBsBMLNTgbcTBa27/z3Vx0hw923A/WZ2KfCcu1+dFMcj7dmnme0BfNTd73f3eWa2MkXhSgyoFURizcz6AKOADQDuXgtcA2TqlqG1Kd6uNT4K/CuF+xMRAdjecIW7/w24Dvixme2ZhmOmNFeb2V6EEepM2+XcAc8C1W3ZiZl1Bu4DuiatTuXfD8kyFdYSa+7+ATAPeNLMJkXrlgCvm1kvM/u+mb1uZqeZ2Xtm9h8zOyjxeTM7x8yuMrPnzOz/ktZPMLPrzOx+M3vIzHpG67uY2Z1mdo2Z/YrwNWabmNn/RDF92swOMbOHo/19P4rxfjMri9avNbNrG+yi1N23RPu6zMy+aWYvmdkP23wCRURa9itCS8iZAGY23My+Z2Z3m9k/omIWM+sc5bKrzGyqmZ2V2IGZfcnMvmNmr5rZV5o6UGM52czOiHLmhWa2wMx+aGYDzezXUT6dambHRrs4F9gLuNzMpkTXwjxrZp9JOsbnzewGM7vPzH5jZt3MrKeZXWtm08zsWDObZ2aLzKy8PScsOifHRCPamNnpZvYTM7vFzJ40s+7Rz/ATM3ss+nvzbnSO9wXONLNvJu2yd/S3aL2Zfb6JYxaZ2Y/MrCI6hpvZXxP/fSQm3F2LllgvwFhgCeFf9fcAe0TriwlJ6gNCr1sv4BlgNuEfjR8Bro627U0Y9T4F2B34U9L+XwJuip7/GKiInvcHdhCSZ0sx3ksY9dmD0Iu3W7S+M/AU8AgwkPAHYRtwLbAbcHD0OrF9on888XP/LHreA7gi2/8ttGjRkrsLcCHw3ybeWwH8NMqrDwHF0fo/AA9Ez28FTo+eXwTURM/PAc6Pno+KcvV+0et7geui503l5D7AGuA2YCIwGbgZuCba9tvAY0mxOlAePe8K/Ae4MHp9IvCHpG0fB35BaH39ZPT34vTo53wSuKGV5+656G/FD4GbCP3Wk5PeXwkMjp7PJvxN6gxcBfwXOBT4XJT3n0vEG22/GLgleu8cYGkTMVwMPB497x+d5xOy/Xulpf6iHmuJPXefbWb7AdcDXwJONbPJHnqua4B17v4ggJldB7wIDCP8Eak1s69Hu3qKkMBPAfomrf9P9NnuwCXRZ3H3FWa2rA2h9gNuBz7n7hujfWwxs/eAKndfHh3nPeB5D/2HMwgJvx/hIshTCH8EADYCF5jZfOAuwj8qRETSwQnF5kRgCHCpmQGsAnpaePEFwqAAwG8IrRAQcu20Brl2APB2g2NcSCM52d0/MLO1wCPu/jqAmS0FNpjZAMIIb/dGg3bfZGbvJ636SrTfhLuBP0fr1xD+XjwUHWNaFGdr/cujHmszu58wmJNwprsvNbPDCQV19yj/LwfecffXgNeizza278eivwn/biamg4BN0c+9wsw+oH5LicSACmuJNTMrBvb2cJHKV6P2jL8DlYQRgIZmR4+9gMHAb939D9G626J9Xk4odG9rcKyDCQlxfTvDHQUcAQwH3kxa37DHcEfdG+61UZJNtGX19dD+grtXmdmFwB3AN4DzgPcREUkhM+tF+LZtNiFvrmokP/YnjKiWArj7DsJgANFnvuvuL0ev6302SaM5OeLUz5X/BS4HFgCvAkOb+RGSPzcS+GfS63cIbS57sGsu3k77Zyx5g/APkIQ1ZvYj4AFC7InqueHP1ZLtNF2bPQHcZWZdCP8tNhFGvyVG1GMtcVfX9wfg7m8CVwJjmti+lJDEFgLvAmckv2lmh0brTzGzrg3Wr4te7tPOWF8GfgL8Pkp8bWJmA0m6EMbMBhFmRBlFuJjxD018VESkIz5D+IbsAUJ+PDwaKQbq8uNKYCtwUtL6/cysBw1yrZmVRAMVDTWVkxvzc2Cxh2lVt7ThZ1lKGOGuOwSwltDqkjLuXhsNfoyN/mHyd+B77j4tlcdpcMzHgMcI36yeAxzt7mvTdTxpHxXWkgsujYrOhEGE3riEPmZWGj0/htBftwb4I/BJM7vDzI6JRhM2EBJgN+BxMzsxulBkDKEYnwF838xKzawfobd5DzMrNrMrzOz8JmK0aLmG8Mfn1kbea7h9w9cfJ/QDJowEprj7OuDSRj4jItIWu4yEmtlRwHeBL7h7DWF0eAXwdzM7NbpA8chohPpPhNlDppjZccBn3H09Idd+1cy+a2ZHEwYYEm10yfmvqZyckDx6PB4oM7PdgMOBrmY2LHpvG+Fiv8QAS/Ixfg58ysx6R68PBX7p7tsJNU+judjMJkVxdW7i3BXRoGaK/vFxCuFbyt6EFsPRwN7AbkkXRjYcFd/aSPz1g2qkXyT6b7E78A/CqHxb/sEhGaJWEMkFewOzzexvhK/JDPhi0vvbgavMbA1hdPdSAHd/xsy+Svg6cQrh4r9ZAGb2CeBOwh+Ke4Bvu7ub2acJUyHNBu4nfKV3IOGGCp8mtJjUuzGNmU0kXJQzmPDV4MPAd6N4ngImAeVmNpwwkrIn8Gkzm8vO0ZsLgW7uXtngZ/+1hRs7bCdcLCQi0mZRHjmXUPz9EqghXADXE/iYu78C4O5bo/z4S0Kue4TQWw3wNcKgxH2Eb+gSM3H8ipD7LiG0rF3i7u9HheMkYLiZ/a6pnGxmZxD6ii82s7nuXk24kPJmQlF9H/AJYD9gEfBbwuj6BRZmixoLnGxmf3f3h8xsJPCgmb1E+Nbzqmhk/Uxgz+jnm024Wc4AMxsbxXk2oW/84aTz1omQ+w8ABkWFtxN6m0+P9jkTeB74N6G95QlCTn+ccMHkAWY2xd0T+72fcLOeD6NCfE/gDDOrIoxEE32+4T0blhIuzpwCdAGKzOxZwgWMmZqCVlpg+m8huczCbWDvdffyDB3v2+6u242LiOSZqC2ls7vH8gZdZvZtwgwti6PXXQjfNvwocW2OZJ9aQSTXNdZmkZ4Dha9M38jEsUREJHOiizMPinFR3QW4GihPWt0P2KCiOl7UCiI5K+pvO5/w1d6Z7v5Amg/5TvQVpYiI5JcSd78r20E0xd03W7it+k+jC9uXE9pWrm7+k5JpagUREREREUkBtYKIiIiIiKSACmsRERERkRRQYS0iIiIikgJ5c/Fiv379vLy8vE2fqampoaysLD0BpZliz7xcjRsUe7a0Nvbp06evdPfc/CHbqT05G3L79yGVdB4CnYdA5yHI1HloLmfnTWFdXl7OtGltu5PohAkT2vyZuFDsmZercYNiz5bWxm5mSzIQTqy0J2dDbv8+pJLOQ6DzEOg8BJk6D83lbLWCiIiIiIikgAprEREREZEUKOjCuqKiItshtJtiz7xcjRsUe7bkcuxxpXMa6DwEOg+BzkMQh/OQ8RvEmNnxwPXAWUn3uz8duAXYHfgdcJm7b4/e6wfcAKwBOgGXu/u2hvudMGGCq79IRHKRmU139wnZjqMxytkiIvU1l7MzOmJtZv2B7sDEpHWDgSnAGcA3gIuBryd97AHgHne/HJgPfC9jAYuIFDDlbBGRtsloYe3uK4BHGqweAlzk7tPd/T7gTuCjAGY2CRjt7q9F2z4GXGpm3TIVs4hIoVLOFhFpm4z3WLt7bYPXUxNfIUaqgaXR82OBqqRtq6OnqfnK9IorYP/94R//SMnuRETyTVxy9voZ65l28DTmfnZuR3clIpI2cZzH+hDgpuj5IGBVg/fXAwMbfqimpoYJE3bm7oqKipab2Kur4e23w6OISAZVVlZSWVmZeNkvm7F0UEZydu2mWj5840OKSgv6mnsRyZLW5uxYFdZmNgKocfc3o1UObGqwWWdglwthysrK2j4peJ8+4XFVw78DIiLplVxImtnKLIfTLpnM2dbJAKjdVtvCliIiqdfanB2bf/qbWTHwBeCKpNXVQO+kbQzoBixPGMEVywAAIABJREFUyUEThfUHH6RkdyIihSLTOTtRWPu2zM5kJSLSFrEprIFvATe7+5akdc8Ao5NelwNbgdTM0aTCWkSkvTKas4s6hT9XKqxFJM4yXlhHIxgAlrTuSmAG0NXMhpnZRWY2wt2nAe+a2bho05OA2919a0qC6ds3PKqwFhFpVFxytkasRSQXZLTH2sx6AOdHLy80szsIXyV+v8Gmc939nuj5WcA1ZraYMJ/qVSkLSCPWIiJNilPOtpKosN6uwlpE4iujhbW7rwfuipaE66Olqc8sAT6XloBUWIuINClOOVsXL4pILohTj3XmaVYQEZGcoFYQEckFhV1Yl5WFx/ffB1eyFhGJK128KCK5oLAL6913h912gw0bYN26bEcjIiJN0Ii1iOSCwi6szWDQoPBcd18UEYkt9ViLSC4o7MIaVFiLiOQAzQoiIrlAhbUKaxGR2EsU1uwA1zUxIhJTKqwThfXy1NwlXUREUs/Mdo5aq89aRGJKhfVee4XHpUuzG4eIiDRLFzCKSNypsB4+PDwuXJjdOEREpFm6gFFE4k6F9ciR4XHBguzGISIizdKItYjEnQrr8nIoKYFly2Dz5mxHIyIiTSjqEv5k1W7WiLWIxJMK606dYOjQcOfFd97JdjQiItKE4t2KAajdqMJaROJJhTWoHUREJAcUdQ1/snZs2pHlSEREGqfCGmDMmPA4a1Z24xARkSZpxFpE4k6FNcCBB4bHGTOyG4eIiDSpaLdoxHqjRqxFJJ5UWAOMHx8eVViLiMRWohWkdpNGrEUknlRYQ2gF6dKF/8/encfHVdf7H399sjb70nRfaQstbWkpTakspWUVUBYRFQW56JUKXjfuz5VFvVwUBcQFRQnqFUQQAVlErYCyr12gpTulK21o0yVbkzRJ8/398Z00kzR7ZnJmMu/n4/F9ZOacmXM+czLzySff+Z7vYdMm2Lcv6GhERKQdGgoiIrFOhTX46fZmzPC31WstIhKTNBRERGKdCutmxx/vf77ySrBxiIhIu5IzQj3WGgoiIjFKhXWzU07xP59/Ptg4RESkXeqxFpFYp8K62bx5/ucrr0BDQ7CxiIjIYQ6Nsd6vHmsRiU0qrJsNHw6TJ0NNDSxbFnQ0IiLSRnKuL6wbKxsDjkREpH39Xlib2Zlm9rqZjQ9blmFmd5jZD8zsV2aW2511Ebdggf/59NNR24WISDyJpZydkp8CQGO5CmsRiU39Wlib2TAgGzi+zao7gZecc9cCjwN3dXNdZJ1zjv/55JNR24WISLyItZytwlpEYl2/FtbOuZ34RHuImY0EPgX8NbToGeB8Mxvf2bqoBHjGGZCeDm+8ATt3RmUXIiLxItZydkqBCmsRiW39PhTEOdf2rJMFwG7nXE1ofSOwCZjfxbrIy8qC004D5+Bvf4vKLkRE4kks5Wz1WItIrIuFkxdHAXvaLKsCRnaxrpWysjKKi4sPtZKSkt5Fc/75/udDD/Xu+SIi3VBSUnIoXwFFQcfTA4Hl7EOF9T4V1iLSv7qbs1P6L6QOOaC2zbJ0oHnOu87WHTJkyBCWLFnS92guvhi+9CV46ik/HGTYsL5vU0SkjYULF7Jw4UIAzGx3wOH0RGA5O7UgFVCPtYj0v+7m7Fjosd4OFLRZlgPs6GJddBQVwbnnQlMT/OlPUduNiEicCixnJ+eFpturaMQ1uUhsUkQkomKhsH4WGGlmGQBmlgqMBZ7vYl30XHaZ//n73/vx1iIi0iywnJ2UkuTnsm5Sr7WIxKYg5rG25psAzrn3gSeB00PLzwAecc5t72xdVIM87zwYOhTeegtefDGquxIRiWWxlrPThqcBUP9+faQ2KSISMf09j3UOcFXo7hVmNiR0+2rgAjO7DjgXWBj2tM7WRcegQXD11f72T34S9d2JiMSiWMzZaSNUWItI7DI3QIY6FBcXu4icvNhs504YOxYaGmD1apgyJXLbFhEJY2ZLnXPFQcfRn3qbs1ddsoqyB8s4+r6jGXapTi4Xkf7XWc6OhTHWsWnYMPjMZ/wY6+99L+hoREQEDQURkdimwroz110HaWnw4IOwYkXQ0YiIJLz0EemACmsRiU0qrDszZgxcFRpe+I1vaIYQEZGANfdYHyg9EHAkIiKHU2Hdleuug7w8+Oc/4fHHg45GRCShpY/xPdYHtqiwFpHYo8K6K0OHwk03+dtf+QrU1AQbj4hIAhs0YRAAtRvbXuBRRCR4Kqy746qrYOZM2LoVbrgh6GhERBLWoDGDsBSjfkc9B2sPBh2OiEgrKqy7IyUF7r4bkpP9vNbPR/fCjyIi0j5LNgaN973WdZvqAo5GRKQ1FdbdNWcOXHutP4HxiiugsjLoiEREEpKGg4hIrFJh3RPXXw/HHQebN8N//qdmCRERCUDGxAwAatersBaR2KLCuifS0uCBByA3Fx5+GG67LeiIREQSTta0LAD2r9wfcCQiIq2psO6po46Ce+/1t7/1LT8Nn4iI9JusGaHC+m0V1iISW1RY98YFF/j5rZua4KMfhaVLg45IRCRhZE1v6bF2BzUkT0Rihwrr3rrxRrj0Uti/H849F959N+iIREQSQmpBKumj02mqa6L2XY2zFpHYocK6t5KS4He/gzPPhF274IwzYMuWoKMSEUkIWTN9r3XVsqqAIxERaaHCui/S0uCRR2DuXD9TyIIFKq5FRPpB7gdyAah8RVOfikjsUGHdVzk5/gTG44/3xfWpp6q4FhGJsryT8gCoeKUi4EhERFpEpLA2s3MjsZ24lZfni+s5c2DTJjjhBHjrraCjEhFp10DI2TlzciAZqt+q5uB+XdpcRGJDjwprMzvbzF4ysw1mtjHUNgOPRCe8OJKfD089BfPnQ2kpnHIKPP100FGJSAIbyDk7JTuF7JnZcBAqX9NwEBGJDSk9fPwfgJuBN4Gm0DIDPhXJoOJWfr7vub7iCvjTn/xsIT//OVx1FZgFHZ2IJJ4BnbMLTiugelk1e/+5l4LTC4IOR0Skx4X1Rufc7W0XmtnaCMUT/9LT4Y9/hDFj4NZb4QtfgNdfh1/9CjIygo5ORBLLgM7ZhecUsu22bexdtJeJt0wMOhwRkR6Psb7azC5uZ/n5kQhmwEhKgltu8QV2Rgbccw+ceCJs3Bh0ZCKSWAZ0zs47OY+krCT2v72fuvfqgg5HRKTHhfVJwG/M7GBYawJ+1ddAzOxoM/ulmV1jZnea2bGh5RlmdoeZ/cDMfmVmuX3dV7/51Kfgtddg4kR/MuNxx8Ff/hJ0VCKSOAZ0zk5KSzo0BGTv3/dGazciIt3W08L6RmAhcCQwIax9LwKx/AH4gXPuJ8APgQdCy+8EXnLOXQs8DtwVgX31nxkzYMkSfxn0igp/CfSvfAXq64OOTEQGvgGfs4suLAJg14O7orkbEZFu6Wlh/Q/gcefcRufcllDbDNwdgViOBnJCt2uBPDMbiT/J5q+h5c8A55vZ+Ajsr//k58Ojj8JPfgIpKf6ExpNP9lPziYhEz4DP2UUfKcLSjPJnyzlQeiBauxER6ZaeFtZ/Ar5tZqeEtfnAtyMQy8P4ryxzgMuALwELgN3OuRoA51wjsAmYH4H99S8z+OpX4aWXYNw4WLzYDw157LGgIxORgWvA5+zU/FQKzykEB2UPl0VrNyIi3dLTwvoG4KvAvcA9Ye3zEYjlC0AdsBiods49AowC9rR5XBUwMgL7C8bcubBsGZx/PpSXw0c+Atdco6EhIhINCZGzh14yFICdf9gZzd2IiHSpp4X1rcAY59x459wRoTYeuCgCsWTgx+jdD/zUzM4AHP4rxnDpQEPbJ5eVlVFcXHyolZSURCCkKCks9D3Vt9/uh4b89Kd+aMjmzUFHJiL9pKSk5FC+AoqitJuEyNlF5xeRkp9C1eIqqpZV9Xo7IiId6W7O7uk81tucc+1lreE93E577gMucc6Vm5nDJ+yvAm1n/c8BdrR98pAhQ1iyZEkEwugnZr6n+sQT4eMf90NDZs6EH/0IFi70U/aJyIC1cOFCFi5cCICZ7Y7SbhIiZydnJjPsP4ax/Wfb2fGrHUy+e3JEtisi0qy7Obun1dsX2i4ws48BP+3hdtpuowiY6ZwrDy26CagEtgAjzSwj9LhUYCzwfF/2F1PmzoU334QLL4TKSrj6apg3D1atCjoyEYl/CZOzR17lR5vsvH8nDeWHdZCLiPSLnhbWhWZ2A4CZZZrZb4Df4MfY9cVeoM7MRoUt2wOsAJ4ETg8tOwN4xDm3vY/7iy2FhX5+6z//GYYNg1degVmz4LrroLo66OhEJH4lTM7OmpJF/mn5NNU0UVpSGu3diYi0y5xzPXuC2QeAS4EPAu8DlwMVzrl9fQrEbCa+d2UpMAx4wTn3vJkVALcAm/FfX37bOXdYtVlcXOziaihIR8rL4VvfgrtCU7+OGAE/+AFcfrmGh4gMUGa21DlXHKVtJ0zO3rNoD2+f8zapQ1P5wOYPkJyRHNHti4hA5zm708LazE7pYNXxwKnAbUATcKFz7pq+BtoXA6awbvbKK356vsWhjqXZs/082PPmBRuXiERcpArrRM/ZzjmWzllK9dJqJv18EqO/NDqi2xcRgb4V1q/hLwJQjk/G7T4MGOGcS+9roH0x4AprgKYmuP9+34O9PfRN6kc+4k9wPPLIYGMTkYiJYGGd8Dm77NEyVl20irRRacxdP5fkTPVai0hkdZazuxpbcCNwpHNuXNhUTW3bePzXjBJpSUlw2WWwbh1873uQmemv4Dhtmp9RZF+fvskVkYEn4XN20QVFZB+bTf32et776XtBhyMiCabTwto593fn3K6uNuKcezhyIclhsrLgu9+Fd96Bz34WGhv93NdHHunHYh88GHSEIhIDlLPBkoyJP54IwNabt1K/UxffEpH+o7Ph4snIkfDb3/rp+RYsgD174Kqr/Pjrp56CHp6IKiIyEBWcVsDgDw/mYPVBNt2wKehwRCSBqLCORzNnwr//DQ89BGPHwvLl8MEPwvz58OKLQUcnIhK4CbdOwFKN0rtLqXi5IuhwRCRBqLCOV2Zw8cWwdq0/mbGw0BfVp5wCZ58NA+1EThGRHsiaksXYb44FYN2V62g60NG5nCIikdPnwtrMZpnZfWb2/8xsfN9Dkh7JyIBvfAM2bvTjsHNy4J//hDlz/Awii/t6HQgRGUgSKWePvW4sGZMzqFlTw5YfbAk6HBFJAJHosb4KOMM592PgxOZL2Uo/y8vzM4ds2uQL7YwMeOwxOP54Px77b3/TGGwRgQTK2cmDkplcMhmALd/fQuXrlQFHJCIDXSQK6y8Cs0O3/wSkRmCb0luDB/uhIe++C1//OuTmwvPPw4c/DMceCw8+qFlERBJbQuXs/FPyGX3NaDgIqz+5msbKxqBDEpEBrMeFtZmdY2aXhm4fj58zdTuAc67JOacugVgwYgTccgts3Qq33urvr1gBl1ziZxFZsSLoCEWkHyhnw4SbJ5A9K5u6TXWsv3o9nV0YTUSkL3pUWJvZ94GHgI8AOOfeAC43s09FITaJhLw8+NrX/BCRu+6CceP8LCLFxX5Mdm1t0BGKSJQoZ3tJ6UlMfWAqSZlJ7Lp/Fzvu3BF0SCIyQPW0x/pEYDQQPqfbP4AfRiwiiY70dFi4EFauhKuvhoYGuPFGmDwZfvELFdgiA5Nydkjm5Ewm/8aPt97w1Q2UP18ecEQiMhD1tLB+2jlXDoR/j3ZxBOORaMvOhjvv9FPzzZwJ27bBl74E48f7sdmVA/5bYZFEopwdZtgnhzHm62NwjY5VF6+ibmtd0CGJyADT08J6o5ldC0w2s0vM7BHgv4AfRT40iaqTT4Zly+CRR+C442DXLvjWt/xQkRtu8PdFJN4pZ7cx4eYJFJxVQMPuBt4+/20aK3Qyo4hETo8Ka+fcn4B/A8nA5UA5cJZz7pdRiE2iLSkJLrrIX0xm0SKYNw/Ky+Gmm3yBffXVsGFD0FGKSC8pZx/Oko2pD0wl46gM9i/fz8oLV+riMSISMT2eFcQ595pz7irn3LnOuf8EnjWz06MQm/QXM39J9Bde8ENEzjsP6urg17+Go46Cj30MXn896ChFpBeUsw+XWpjKjEUzSBueRvlz5az59BrcQc0UIiJ9l9LZSjN7AsjuYhtFwG7gX5EKSgJ08sm+rV4NP/4x/OEP8PDDvp1wAlxzjb+iY0qnbx0RCYBydvdlHJHBjEUzePOUNyl7qIz1hes56s6jsCQLOjQRiWNdVUcrQo9ZBzQBVwJvAkvCHnMEUBiV6CQ4U6fCb3/rZw75+c+hpARefdW3ceP8iY4f/7jv7RaRWKGc3QPZM7OZ/vh0Vpy9gtK7SjEzjvzlkSquRaTXuiqsbwUsdFY5ZjbROfedtg8ys8ejEZzEgFGjfBH9ne/APffAz34G69f7C8385jfwk5/A9OlBRykinnJ2DxUsKOCYJ45h5QUr2fHrHTjn1HMtIr3W6Rhr51xFc4IOGWZmyeGPMbMFwPFRiE1iSVYWfOELsGaN770uLIRnnoEZM+CKK2DjxqAjFEl4ytm9U3hWIdOfmE7SoCRK7ypl/VXrNeZaRHqlpycvPgwsN7Nfm9mtZvYk8DTw88iHJjEpKQmuvBLWrYMvftGPtb7nHjjySLjsMn8BGhGJFcrZ3VR4ZlhxfXcpqz+5WrOFiEiP9XS6vaeBs/AnvkzGT910sXPu5ijEJrGsqAjuuAPWrvU91klJ8Mc/wjHHwAUXaBYRkRignN0zhWcWMmPRDJJzkyl7qIwV566gsVLzXItI95lzffu6y8ySgFnOuaURCcjsA8DJwAbgJWA/cAtQBRQA33TOHXZ5wOLiYrdkyZK2i6W/bNkCt93mx13Xha5mdtppcO21/qdOchTpkJktdc4V99O+lLO7UL28mhVnr6D+/XqyZ2Uz4x8zSBuWFnRYIhIjOsvZPSqszex37SweAhx0zl3Yy/jCt/8ZYKJz7vqwZf8HLHLOPWhmZwP/4Zz7ZNvnxnKSTig7d/oTHH/5y5bLox9/vC+wzzvP92yLSCvRKqyVs3uvdmMtKz64gtoNtaSPTeeYJ44he2ZXMxmKSCLoLGf3tMo5EtgKbAlr6UCfez7MbD5wKXBD2LKRwKeAv4YWPQOcb2bj+7o/iZJhw+AHP/A92N//vh8y8sYbcOGF/kTH++6DhoagoxRJFMrZvZQxIYNZL88i9wO5HNh6gGUnLqPs0bKgwxKRGNfTwvpzzrnvOef+p7kBXwayIhDL7cAa4A4zW2RmJwELgN3OuRoA51wjsAmYH4H9STTl5/te6i1bfA/26NGwahV8+tMwYoS/XPpzz6nIFoku5ew+SBuaxsxnZzLs08Noqmli1UWr2Py/m+nrEEoRGbh6evLiunYWHwA+15cgzGwyMAu42zn3RfwVwRYBY4A9bR5eBYxsu42ysjKKi4sPtZKSkr6EJJGSmQlf/jK8+66/4Mz06bBnj79c+qmnwuDBcNFFfgq/rVuDjlak35SUlBzKV/irIUaccnbfJQ9KZso9U5hwywQw2Pydzay6aBUN+9QpIJJIupuzezrGehMQ/oRkYBjwB+fclb2MFTO7CChxzhWF7mcAZcCNwEedc3PDHrsMuN85d1v4NmJ9vJ6EOAdvv+1nEPnrX/282OGmToWzz4azzoJ583xhLjLARXGMtXJ2BO1+cjdrLlvDwYqDpI9LZ9qD08idmxt0WCLSzzrL2V1debGtp4H7aUnUTcBO59z6PsTXHMehWJxztWb2Dv6PQEGbx+YAO/q4PwmKmR9rPWOGv6Ljli2waJFvzzwDq1f7dvvtkJbmi+uzzoIzz4SZM3Xyo0jPKGdHUNGHiyheVszqS1ZTtbiKN09+kyN+cARj/nsMlqyZj0Sk52Os73POPeecez7UXgTMzCb1MY4VQJ6ZDQtb1ohPxiNDvSGYWSowFni+j/uTWDFuHHz+8/Doo36IyLPPwre/DbNn+/HX//oXfPObcNxxMHw4XHqpvyDNjgHxd1ok2pSzIyxjQgazXprF6GtG4xodG7+xkTfnv0nNOzVBhyYiMaCnhfXZ7SzbCTzUlyCcc2uBfwAXA5hZPpAG/BF4Ejg99NAzgEecc9v7sj+JUWlpsGCBn1VkyRI/dd8DD8BnP+tPfiwrg/vv9xekGTUKpk2D//oveOgh2LUr6OhFYpFydhQkpSUx6fZJHPPkMaSNSKPy5UqWzFzCez9/D9ekExtFElm3xlib2ReBm/Bf6bXnaedcewm8+4GYFQE/A5bgezjuds6tNrMC/MUGNgPDgW8756rbPj8ex+tJDzjnr/L49NPw1FO+Z7umTQ/R1Km+MF+wAObPh6FDg4hUpMciPcZaObv/NOxtYMNXNrDzvp0A5J6Yy1G/OorsGZrzWmSgisgFYsxsBvAfwM/brKp3zpX2LcS+GyhJWrqpvt7Pj/3cc/D88/Dyy1Bb2/oxU6f6Aru50B42rL0tiQQuGicvKmf3r7JHy1h/9XoadjZAMoz+8mjG/894UnJ6eiqTiMS6SF55MdU5d9gcQ2Y21jkX6FxpAy1JSw/V18Pixb7Qfu659gvto4+GuXP9+O3Zs/3JkJp1RGJAFGcFUc7uRw3lDWy+YTPb79wOTZA2Mo2Jt05k6CVDsSSd3CgyUPS6sDazrwF/d86tDt2/Dn/Wd6uHAac65xZEJtzeGYhJWvogvNBu7tFuO3QkOdkX27Nnw6xZLbOVDB4cSMiSuCJVWCtnx4aqZVWsv3o9VW9UAZBTnMPE2yaSPz8/4MhEJBL6Mt3eEUB4JjgOf8Wu92mZvsnw4+hEYkdaGpx0km/XXecL7WXL/EmRS5f6tno1rFzp2z33tDx31KiWIru5TZ4MqanBvR6R7lHOjgE5x+Vw3KvH8f7/vc+mGzZRtaSKtxa8xeDzBjPh5glkTYvEhS9FJBb1dCjIccBK51x9m+VTm3tIgjKQez8kSmpqYMUKX2QvX+5vv/324T3b4IvqqVN9kT1zZkvBrXHbEgFRHAqinB2wg/sPsu32bWz90Vaa9jeBwZCPD2H8DeNVYIvEqUiOsf4P59w9bZZNBo5yzv21b2H2TaIkaYmypibYuNEX2c1t+XK/rD1Dh7bu2Z4+HaZMgSz9wZTui2JhrZwdI+p31rP5xs2U3l2Ka3C+wL54CONuGEf2MZpBRCSe9LmwNrN5+K8PFwIlbVYXAj93zo3ta6B9kWhJWvpZVZUfMhJecK9YAZWV7T9+/Hjfwz1tmv85daofz53T0exnksiiMN2ecnaMqttWx9YfbfUFdr3/+1t4biFj/nsM+aflY6aTHEViXSQK62OAe/Hj9/a1WV0P/NE5d2NfA+2LRE3SEiDnYOvWlmEkK1b4cdvr1/urRrZnzJjWxXZzy8vr39glpkShsFbOjnEHth84VGA31TUBkDUzizH/PYahlwwlKa2n128Tkf4SqXms84CPOed+E8ngIiXRk7TEkIYG2LDBF9nNbdUqWLfOn0TZnlGjWorsKVNa2rBhoB6sAS9K81grZ8eB+t317Pj1Drb/YrufAxtIG5HGyKtGMuI/R5A+Kj3gCEWkrYiNse5kByc7517q84b6QElaYl5jox+r3VxoNxfda9dCXV37z8nNbSmyJ09uuT1pkp/5RAaEaI2x7mR/ytkxpulAEzvv38l7t7/H/pX7/cJkGPzhwYz8/EgKzyrEkvVPtkgs6Ms81k8AXZ1VYcA051yg149Wkpa4dfAgbN7cumd77Vrfysvbf05yMhxxROve7ebCu6ioX8OXvovgPNbK2XHOOUf5v8vZcdcOdj+6G9fo/0anj0tn5JUjGXb5MAaNGRRwlCKJrS/zWK8DyoHttMyBetj2gQ/3PjyRBJecDBMn+nbeeS3LnYOyspYiO7zg3rzZDzfZsAGefLL19gYPbt273Vx0T5gAKbq88gCnnB3nzIyC0wsoOL2A+p31lP5fKaUlpdRtqmPT9ZvYdMMm8hfkM+zTwxjy0SGk5OozLRJLuuqxLgKqnHMHOt2I2WDn3J5IB9cT6v2QhFJX54vq9oru6ur2n5Oa6oeQtFd05+uKcEGKYI+1cvYA5Joc+57ZR+lvStn9xG7cAf93OykjiaILihj26WEUnFVAUopOeBTpD73usXbO7e7mPj4CxOQJMiID0qBBfs7s6dNbL3cOSkvbL7i3boU1a3xra9iww8dxT54M48b5HnWJC8rZA5MlGYVnFVJ4ViEN5Q2UPVzGzj/spOKFCnb9aRe7/rSLlMEpDPnIEIZ8bAj5p+aTlKoiWyQIfTp50czGAv8FfNE5F+gVMdT7IdKF/fvhnXcOL7rXrYPa2vafk54ORx3VftGdrYtaREp/nbyonD2w1G6uZdcfd7Hzvp3UrG25YmxKYQpFFxYx5GNDKDitQFP3iURYRGcFMbNk4HzgSuBMYGdoO6P6GmhfKEmL9FJTE2zb1rp3u/n2jh0dP2/UqMNPnJw8GUaPhiT9Ie+JaBbWytkDn3OO/Sv3U/ZwGWUPlVGzJqzILkih8NxCis4rovDsQlLyNCZbpK8iNY/1EcDngM8Ag4EtwJeBfwInOOdeiUy4vaMkLRIFlZX+gjfNBXdz0b1+fcdzcmdmHt7LPXmyX6ZLvbcrSvNYK2cnqP2r9rProV2UPVxGzaqWIttSjLxT8hh83mCKzisiY2JGgFGKxK8+FdZm9jF8T8dpwAb8uLx7gGudc9dEONZeU5IW6UfNUwSG93I3F927dnX8vDFjDi+4m3u5E/hCOJEsrJWzJVzNuhp2/3U3e/66h4qXKqCpZV3m0ZkMPm8whecUkndinoaMiHRTX6bbA6jFXwL3FeBzzrn1oY32/coyIhKfwqcIPPfc1uv27fMFdvgY7rVr/Swm27b59vTTrZ+TldW6l3vqVDjmGDjySJ082XPK2XJI5uRMxk4ey9ivjaVhbwN7/7GX3X/dzd5Fe6lZU0PNmhq23bKNpKwk8hfkU/hBf5JkxlEZWAL/syvSWz0ZCjIW3wurND5gAAAgAElEQVQyHngcmO+c+1JonblIXMKxD9T7IRLjGhth06bDC+516/x83e0ZNAimTYMZM1raMcfAkCH9G3uURWkoiHK2dKipoYmKFyvY87c97HtqX8vVHkPSx6ZTeFYhBWf5ObVTC1MDilQk9kT65MUU4ELgUmAl8FvgLOdcSV8D7QslaZE4tndv617uVatgxQo/RWB7hg8/vNg++mg/i0kcivLJi8rZ0qUD2w+w9+m97HtqH/ue3kfD7oaWlUmQPSub/AX5FJxaQN7JeToJUhJaRAvrNhueBHwJuNI5l9nrDUWAkrTIAFReDitX+iK7ub39dvsXwUlO9kNJwovtGTPiYvx2P063p5wtXXJNjuq3qtn7T19oV7xcgWsIqxWSIOe4HPIX5JN/ar4vtHUFSEkgUSusw3ZwuXPu3j5vqA+UpEUSRFOTP3EyvNBescLP0d1ePhs8GGbP9q242P8cOzamiu3+KqzD9qecLd12cP9BKl6poPy5csqfK6fqjSpcY5tCe3ao0F6QT+4JuaQWaOiIDFxRL6wjxczSgMXAV5xzz5lZBnALUAUUAN90zlW291wlaZEEt38/rF7dund7xQo/zKStoqLDi+0xYwIrtvu7sI6k3uZt5ez41WWhDWROzSTvpDzyTsoj96RcMibqZEgZOPo6K0h/+jr+RJtmdwKLnHMPmtnZwF3AJ4MITERiXFYWzJnjWzPnYMsWWLq0pS1ZArt3wz//6VuzoUPhhBPgxBP9z+JiyNA8v92gvJ1gkrOSKTyzkMIzCwForG6k8pVKyp8rp+LFCioXV1Kzuoaa1TWU3l0KQOqwVPJObCm0c47L0fR+MiDFTI+1mZ0ITAG+A1wBrAc2AQXOuZrQCTgVwDTn3Oa2z1fvh4h0S9tie8kS/7Ntz3ZKCsya5QvtE0+EM86AwsKohBSvPdZ9ydvK2QNX04EmqpZVUfFyBZUvV1LxcgUNZQ2tHmPpRu6cXHJPzCX3+Fxyjs8hfXS6erUlLsR8j7WZZQEXO+f+28y+E1q8ANjtnKsBcM41mtkmYD6wOZBARST+mcH48b599KN+mXPw7rvwyivw6qv+58qVsHixbz/7mb9M+8knw4UXwqWX+h7uBKa8LR1JSk8i74Q88k7Ig6/5S67XbqhtVWjXrKmh4qUKf9GakLQRaeQcn0Pu8bnkzs0lpzhHs49I3ImVd+w3gR+2WTYK2NNmWRUwsr0NlJWVUVzc8s/DwoULWbhwYSRjFJGBygwmTfLt8sv9sqoqeOMNX2Q//7xvL7zg2ze/CRddBNdcA3Pn9mqXJSUllJQcmvGuKCKvo3/1KW8rZycOMyPzyEwyj8xkxBUjAGjY00DFqxVUvVFF5euVVL1RRX1pPXse38Oex1veQplTMn2xPdf3amfPyNYQEglEd3N24ENBzOyDQLpz7onQ/c34rxSLgY855+aGPXYZcL9z7ra229HXiiISVRUVsGgR3Hcf/P3vfnYSgLPOgptuaj22u4fibShIJPK2craEa+7Vbi6yK9+opPrNalx96xrF0o2cWTlkz84mZ3YOObNzyDw6k6RUFdvSf2J9KMjXgFlh46oK8FcJuz10O1wOsKP/QhMRCcnLg098wrdt2+DOO+EXv4CnnoJnnvG92Dfd5IeMDHzK2xJR4b3awy8bDvix2tUrqlsV27Xraql8rZLK11ommkkalETWjCxyZueQfZwvuLOmZalnWwIRCz3WQ4Dwy6W9Cvw38DqwGhjinKs1s1SgGpjgnNvedjvq/RCRfrd3L9x8M9x+u+/BvuYaf7uH4rDHus95WzlbeqOhvIGqxVVUL6umamkVVUurqNtYd9jjLM3IOibrUK929nHZZB+TTVK6im3pu5jusXbOlYXfN7ODQJlzbquZPQmcDjwJnAE80l5RLSISiMJCuPVWPxzknHP8SY5f/zqMGBF0ZFGlvC1BSc1PbTXVH0DDvgaq3/SFdnPBXftOLdVLq6leWk0pfso/SzWypmf5Xu3jcsg+NpusY7JIyQm8FJIBJNbfTVcDt5jZTGA4oDNbRCT2zJwJ2dl+HPbevQO+sO6C8rb0q9SCVApOK6DgtJZRSI2VjYeK7aplVVQvraZmXQ3Vb1ZT/WY17//2/UOPzZiUQfax2b7QnplF9rHZpI/S1H/SOzFXWDvnxofd3gdcGVw0IiKdaGyEe++F667zRfWcOTB1atBR9TvlbYk1Kbkp5M/PJ39+/qFljVWNVC/3vdjVb/m2f9V+ajfUUruhlrKHW76ISRmcQvbM7EMFd/bMbJ0kKd0Sc4W1iEjM27oV7rkHSkrgvff8snnz4M9/Duyy6CLSuZScFPJPzif/5JZiu6m+iZo1Nb7gfqv60M/GPY2U/7uc8n+XH3qspRlZ07IOFdrNPdyp+alBvByJUSqsRUS60tAAy5b52T8efdRfqbHZUUfBd78Ll1ySKDOCiAwYSWlJvkiemQ2hKeydcxx470BLsf1WNfuX+57t5qEk4dLHpftie4bfTtbMLDImZGBJ+ic7EamwFhFpq7S05XLnr70GL70E+/e3rM/Kgg99CK68Ek47TQW1yABiZgwaM4hBYwZR9OGW64A0VjWyf8X+1gX32/s5sOUAB7YcaHVhm+TsZLKOyTpUaGfPDJ0oma2ya6DTb1hEEldNDaxdC6tXw6pV/jLmy5bBjnamXT7qKJg/Hz78YTjzTMjI6P94RSQwKTkp5J2UR95JeYeWNTU2Ubu+lurl1S1F9/Jq6nfUU/lqJZWvVrbaxqCJgw71kDcX3YPGDdKJkgOICmsRGdgaGmDLFnj3Xd82bPBt1SrYtAnam8s/Nxdmz/atuBhOOSXRZ/oQkXYkpSSRNTWLrKlZ8MmW5fVl9a0K7erl1dSsrqHu3Trq3q1j9192H3pscl5yq2Ek2TOzyZqWRXJmcgCvSPpKhbWIxLfGRt/DvHWrb9u2+YK5uZDeuhUOHmz/uSkpvid66lSYNs3/nDULJk7U8A4R6bW0IWmknZ5GwektUwA21TdRs7aG6hV+zHZzwd2wq4GKFyuoeLGiZQNJkHFkxmG925oGMPapsBaR2OUc7Nnji+Xmojm8gN661RfVTU0db8MMxo71xXJzmzTJF9GTJkFaWv+9HhFJWElpSb5nekY2XNay/MD7B1oV2tXLq6lZW0Ptulpq19VS9uewaQAL/TSAWTOyDhXcmVMzSR6k3u1YocJaRPpfUxPs3u2L4tJS38JvN99//32or+98W2Z+mMbYsb6NGQPjxrUU0OPHQ3p659sQEQlI+vB00oenU/jBlqtJHqw7SM3qw3u3G/c2Uv5sOeXPtkwDSDJkTsls3bs9I4u04Wnq3Q6ACmsR6Tvn/AVSdu2CsrLWre2yXbt8a2zs3rbz8nyxHF44h98eNUq9ziIyoCQPSibnuBxyjss5tMw5x4Hth/du175TS82qGmpW1bDr/l2HHp86JLXVuO3smdlkTskkKU3D3KJJhbWIHK6mxvco79nT0sLvty2Yd+/2Jwn2RGGh72keOdL/bO/28OGQmRmd1ygiEkfMjEGjBzFo9CAGf2jwoeUHaw6yf5UvtluN3S5rYN8z+9j3zL6WbaQamVMzDztZMm2IOiciRYW1yEDW3JPcXnHcWeFcV9fzfWVnw9ChMGRIS+vo/tChMGhQ5F+viEiCSc5MJndOLrlzcg8tc85Rt6WuVaG9f4W/yM3+5fvZv3w/O/+w89Dj00akHda7nXFUBkkp6t3uKRXWIvGisRH27et+cdzcOpoRozPp6TB4cEsrKmp9v73CWYWyiEhMMDMyxmeQMT6DogvCLnJT3cj+t9v0bq+opr60nr2le9m7aG/LNtJDl3Cf2bp3O7VAl3DvjAprkf7mnL+KX3Phu3fv4cVwe628vOtttyc7+/DCuKv7WVn+pEARERkwUrJTyDshj7wTWi5y45octRtrDxu7fWDLAaqXVVO9rM0l3Mekk31cNrlzc8k9PpecOTmk5KqcbKYjIdIXjY2dF8Ydretqpov2mEFBQfeK4+bbhYWaEUNERDpkSUbmpEwyJ2Uy5KNDDi1vrGikekV1q97t/Sv3c2DbAQ5sC7uEu0Hm0Zm+yJ6bQ+7cXLKmZ5GUmpjDSFRYi4DvRa6q6n5h3NwqK7vednsGDWpdFDe3wsL2lw8e7IvqZM1VKiIi0ZeSl0L+vHzy5+UfWuYOOmreqaFqcRVVb1RR+Xol1W/5q0rWrK7h/d+/D0BSRlJLr/YHcsk7OY/0EYnRyaPCWgae+vruD7FofszevT2f1QL81fnCe5G7Ko6b12mmCxERiTOWbGRNySJrShbDPz0cgKYDTVS/VU3l65VUvlFJ1etV1G6opfLlSipfbul8ypiUQd4peeTNyyN/Xj6DJgwakPNsq7CW2NV2Rovu9iZXV3e97fZkZXVdGLddn5+vS1+LiEjCSkpP8j3Tc1tmJWnY00DlYl9kV7xSQcXLFdRuqKV2Qy3v/873aqeNSCPvFF9k55+eT+bkzAFRaKuwlv5RV9fzYRb79vVuRovk5O71GrddplktRERE+ix1cCqDzx7M4LP9fNtNjb5Xu+LFCipeqKD8xXLqS+spe7CMsgf9JdvTR6dTcFYBhWcVkn96PmlF8Tm3tgpr6Z0DBw6/kl5Hbfduf8GR3sjJ6flY5NxczWghIiISI5JSksgtziW3OJcx14zBOUfNmhoqXqyg/Ply9j2zjwPvHeD9373ve7QNcmbnUPDBAoouKCKnOCduerNVWIvX1OR7icML4s6K5oqKnm0/NbVnxXHzyXq6VLWIiMiAYmZkTc0ia2oWIz8/EtfkqF5Rzb6n9rH3qb1UvFRB1ZIqqpZUsfX7W0kblUbRhUUUXVhE/vz8mJ5xRIX1QHfwoC+ES0t927Gj/dvvv++njuuu5GR/cZDw1nzBkLatqMjPpRwn/22KiIhI/7EkI+fYHHKOzWHsN8ZysOYg5S+Us/dve9n92G4OvHeAHb/cwY5f7iClIIUhFw9h+H8MJ/fE3JjryY6ZwtrMLgJ+DOQC9wH/zznXaGZFwPeBciAV+KZzrhfTNwxA1dWwbRts3erbe+8dXjDv3Ol7o7ujoACGDeu6UB46VCftiSQ45WwRiZbkzORDY7Qn/XwSVUuq2P3YbnY/upuaNTWU3l1K6d2lDJo4iOGXD2f4Z4YzaExsnCcVE4W1mY0FLgQuBqYBvwK2AbcBDwHfcs69bmZXATcC3w4q1n5z8CBs3+4L5vDiOfz2vn3d29aQITBiREsbOfLw28OH6+Q9EekW5WwR6S9mRu6cXHLn5DLh+xPYv2o/79/7Pjvv20ndu3Vs/u5mNt+4mSEXD2H0V0eT94G8rjcazXidc4EGAGBm84BXnXONofu34JP1/wJ/cc6NDC0fBawDhjnn9odvo7i42C1ZsqR/A++rAwdg82Z4913YsMG35tubNnU9r3J6OowdC2PG+J+jR/tCObxwHjZM45RFYpyZLXXOFQcdR3clbM4WkZjhDjr2/Xsfpb8tpezhMghNIpZ7Qi5H3HQEBacVRG3fneXsmOixds692GbRdiAHOA3YFPa47aGxNMXA8/0WYF9VVMCqVa3bO+/43ufOhmkMHw7jxrUunsNvDxmiccsi0u8GfM4WkZhnyUbhmYUUnllI3Xt1fgx2yQ4qX61k+enLKTijgAm3TiDn2Jx+jSsmCut2zAFuBRYCe9qsqwJG9ntE3dHU5HubFy+GZctaiuj33mv/8UlJcMQRMHEiTJrU+ueECf6CJSIisS8+c7aIDAiDRg9iws0TGHf9ON772XtsvWUr+57Zx9LipYz9+ljGfWccyRnJ/RJLzBXWZjYJKHPOLTczB9S2eUg6cNgYibKyMoqLW3rlFy5cyMKFC6MaK+Xl8MIL8NprvphevLj9aegGDYKjj4Zp01ra5MkwfryGaYgkqJKSEkpKSprvFgUZS1/EVc4WkQEtOSuZcdeOY+RVI9n8P5vZfsd2tv5wK3ue3MP0x6aTMTGj19vubs6OiTHWzcwsGfghcL1z7oCZfRs41Tl3Vmi9AXWhZa+EP7dfxus1NPhC+pln4F//gqVLDx/KMWIEzJkDxcVwzDG+iJ4wwU9PJyLSjngbY90s5nO2iCS0ilcrWHvFWmrX15JSkML0J6aTf3J+n7cb82Osw3wNuM05dyB0/2ngqrD144F6oP+ycUODL6Qfeggee6z1TBwpKXDiiTBvHhx/vC+oR43qt9BERAIWezlbRCQk74Q8Zr8xmzWfXsOev+7h7XPeZsZTM8g7IXozh8RMYW1m1+KTb4aZTQAWAC8ApWY23Tm3EjgbuMM5Vx/1gHbuhJIS+PWv/bzQzY4+Gs49F04/3RfU2dlRD0VEJNbEXM4WEWlHSl4K0x+dzprL17Dr/l2svGAlxW8Vkz4yPTr7i8pWe8jMrsdP0xRurXPud2b2CeA7ZrYZyAaui2owFRXwwx/CT38KdXV+2eTJ8KlPwcUXw9SpUd29iEisi6mcLSLSBUs2ptwzhfqd9ZT/q5y1V6xlxj9nROWqjTE1xrovIjJeb9EiuOIK31sNcP758OUvw2mnaVo7EYmaeB1j3RcaYy0i/a1+Zz1vHP0Gjfsamf7EdIrO6915453lbF2TutmPfgTnnOOL6hNP9DN9PP64H/KholpEREQkrqUNS2PcDeMA2PqjrVHZhwprgFtugW99yxfQ3/++n/lj7tygoxIRERGRCBpx5QiSc5KpfLmS/Wv2d/2EHlJh/dxz8O1v+9v33gvXXqup8UREREQGoJTsFIou8kNA9jzZ9npWfZfYhfXBg34MdVMTXHcdXHZZ0BGJiIiISBQVnlUIQPm/yyO+7cQurBctgrffhrFj4frrg45GRERERKIs9wO5AFSvqI74thO7sL7vPv/zqqv8ZcdFREREZEAbNH4QSRlJ1O+op6G8IaLbTtzC2jl/WXKAj3402FhEREREpF9YkpExMQOAus11Ed124hbWGzdCWRkMHQpHHhl0NCIiIiLST1KHpQLQsFM91pGxfr3/OWOG5qkWERERSSBpQ9MAqN9VH9HtJm5h/e67/ueECcHGISIiIiL9KnVIqMd6t3qsI6OpCYYPh0mTgo5ERERERPpR2rA00kakYSmRHbWQEtGtxZMvf9k354KORERERET60bhrxzHu2nER327i9lg30/hqEREREYkAFdYiIiIiIhGgwlpEREREJAJUWIuIiIiIREBCF9YlJSVBh9Brir3/xWvcoNiDEs+xxyodU0/HwdNx8HQcvFg4Diqs45Ri73/xGjco9qDEc+yxSsfU03HwdBw8HQcvFo5DQhfWIiIiIiKRYm6AzONsZmXAlh4+rQjYHYVw+oNi73/xGjco9qB0N/Zxzrkh0Q4mlvQyZ0N8vx8iScfB03HwdBy8/joOHebsAVNYi4iIiIgESUNBREREREQiQIW1iIiIiEgEqLAWEREREYmAhCyszSzDzO4wsx+Y2a/MLDfomDpiZheZ2SYz22NmPzOzlNDyIjO7y8x+ZGa3m1lq0LF2xMzSzGy5mS0I3Y+n4/8BM/uamV0YOuYxH7uZHW1mvzSza8zsTjM7NrQ8JmM3szPN7HUzGx+2rMNYY+l1dBB7u5/Z0Lq4+dzGmlj6vQetvfddounsc5ZIzGy+ma0zs31m9tOg4wlS21ojMM65hGvA/wGfCN0+G3gg6Jg6iHMscC8wG7gc2A98LbTuWWBu6PZVwM1Bx9vJ67gOqAAWxNnx/wxwU7y9d4AlwKiw99CaWI0dGAZ8BHDA+O4c51h5He3F3tlnNrQ+bj63sdZi5fcedOvoM5NIravPWaI0IA+4FigAPgQ0AGcEHVeAx6NVrRFYHEEfiAAO/EjgAJAZup8S+lCODzq2dmKdB6SE3b8F+BvwAWBH2PJRQDWQFXTM7byGE4HPApuBBfFy/IH5wDOEZs6Jp/dOKKYpodtDgB2xHDv+m7Pw4rTDWGPtdbQTe7uf2dDtuPncxlqLtd970K3t+y7RWmefs0RqQEabv1FLgVODjiugY9Gq1ggylkQcCrIA2O2cqwFwzjUCm/CFVExxzr0Yiq/ZdmArcBo+5ubHbQ/dLO7H8LpkZlnAxc6534UtXkB8HP/bgTXAHWa2yMxOIn5ifxj4jZnlAJcBXyKGY3fONbVZtICOY+1sXb9rG3snn1mIk89tjFpADP3eg9bOZyahdPE5SxjOuVoXqipDf2/XAM8FGlQAOqg1ApOIhfUoYE+bZVX4HpFYNwf4NfHzGr4J/LDNspiP3cwmA7OAu51zXwT+BSwCxhDjsYd8AagDFgPVzrlHiIPjHqazWOPpdUDLZxbiL/ZYomMnnQn/nCUcM5sH/ANIBTIDDicI7dUagUnEwtoBtW2WpePHJsUsM5sElDnnlhMHr8HMPggscc7tarMq5mMHpgF7nXMrQvd/ARjxETv4rwcfAO4HfmpmZxA/sUPnscbN62jzmYU4ij0G6dhJu9r5nCWiDcDvgdOBW4MNpX91UmsEJhHPot2OH+gfLgc/DjUmmVky8HngW6FF24GjwtYbkEVsvYavAbN8aIA/5o/jh1jE+vFPIeyz4ZyrNbN3gGRiP3aA+4BLnHPlZubwRfZXiY/YofPPqOtkXcxo5zML8fG5jVVxl7cl+jr4nCUc51wp8DszawK+HnQ8/azdWsPMvu+cuyWIgBKxx/pZYKSZZQCEprsaCzwfaFSd+xpwm3PuQOj+08DksPXjgXr8bBCx4lPAsWFtB/A5/Jn9sX78VwB5ZjYsbFkjoZMAYzl2MysCZjrnykOLbgIqgS3EeOxhOvuMxsvnt+1nFuLjcxur4uX3Lv2rvc9ZIluC/yc0kXRUawQ2NCjhCmvn3PvAk/ivTADOAB4JO5EoppjZtcCbQIaZTTCzzwLlQKmZTQ897GzgDudcfVBxtuWcK3POvdfcgIP4r+u2EuPH3zm3Fj9e7WIAM8sH0oA/EuOxA3uBOjMbFbZsD/6fhZiM3Vq6Ggw6/4zG2ue3beyhZYd9Zs1sknNuCTH+uY1VsfZ7D1p777tE09HnLOi4+pOZZZrZrLBF5wI/CyqeIHRSa1QGFVMiDgUBuBq4xcxmAsOBhQHH0y4zux743zaL1zrnfmdmnwC+Y2abgWz8/I3xIh6O/+XAz8wsDd8z9knnXKOZxXTszrkmM7sQ/95Yip/z9uvOucpYjD1s5hKAK8zsF865Mjp/j8TE62gvdvzX0u1+ZkO34/lzG7SY+L0HrZPPTMLo7G9jEPEEaDrwNzNbC7wGvOGc+1vAMSU8C83UIiIiIiIifZBwQ0FERERERKJBhbWIiIiISASosBYRERERiQAV1iIiIiIiEaDCWkREREQkAlRYi4iIiIhEgAprEREREZEIUGEtIjHFzIaY2RNmtsXMvhV0PCIi0jnl7Ra6QIyIxBQz+yDwLJALrAWOcs7tDTYqERHpiPJ2CxXWIhKzzOwF4DTnXGPQsYiISNcSPW9rKIiIxCQzGwM8lajJWUQk3ihvq7CWKDCzb5jZL4OOo1kk4jGzM8zsVTNbEKGwBhQzuy3C2zPg48API7ldEWlN+brvzGy8md1tZndFaHvJZnalmW3pxmM/ZGav9/a1RjJ3K297KqzjjJkdbWZ/MbPfmNlyM3NmdlnQcbXxd+CBoIMIE4l4tgEzIxBLTDKzc81st5mVmtkJZjYm9D5zZvYLM8sPPe54M1tpZj82s4zQsmHA+xEO6VPAb51zjaHti8Qd5eteicd8vSv0Mz1C20sCyoGx3XjsEuDojlaa2TQzezD03rvezIaHrYt07lbeRoV1PHoM+B/n3OeAWcC9zSvMbIaZzQ8iqPB9O+dWOudeCiKO9vQkno6OoXNuHbA7knFF8vfV12055/4O3Ak0Oudedc5tA64E6oHXnXPloce9gU/kX3fO1Yaefh7wZJ9eQBgzuwH4HvCMma0FJkdq2yL9TPm6h+IxXzvnaoAdkdqXc64BWNrNx+7EF+EdrV8FlITu3uacCy+kI5a7lbdbqLCOI2ZWCBwF7AdwzjUB3wGcmeXhk7YFEFdg+46kbryOpn7cVxDbehAYbWazAZxze4B/AR8O21c+sCP03ms2xTm3to/7PsQ597/OuSOdc8XOuSnOuRcitW2R/qJ8HV0xmK8jPRNET+Lv6rEHQz/bjnuOWO5W3m6hwjqOhKauWQcsMrMPhJZtAd7Aj2s6AlhoZp8BMLMbzeyG0NdAh8ZRmVmxmd1uZl81sxozW2dm/21mF5rZYjP7iJk9amaVZvb5sOe1u702+77SzH4fPtYs9HXoHWZ2s5k9a2bHh5Z3ur+w55uZfd7MVpnZmWa2wsx2mtkV3djHkeHxdLHPw45hOwpCX+1WmdmVbY7p90LH5i+hRNzdY3bYvszsYjN7w8yuMLN3zOyHnWyvvd99u/F0JtSzsQr4WNjiIcCHLDTsAzgfeCIszgygJuz+FDP7oZndZmbPmx9S8gkzW2JmF5jZA2ZWbmZfMbPZode4x8zO7So+kXiifJ04+bo9Hb3GzvZjZkeY2e/M9/7e2p39tLPfv5vZH83spC4ep9wdLc45tThqwDRgC/4/1N8BQ8PWbQYWhG5PBzaGbg/B/zc9GP8f91agOLTul8ArodtpwF7gF/i5KD8JbO9se233jR9jdhfw+9DyQfhiLT90/2KgDMjpbH/tvO7hoX1eG3rsTUAdMK6LfbSNp9N9hh/DdmLYDPwYyAw9b2toeS7wp7DHvYxPit06Zh3sqxD/9d5PgeOBMzrbXpvffbvxdPP9dT3wbuj2ROBm4ABwUdj7xcIefz4wL3Q7E1gMZIbuvxA6DrlABXAjkAF8InT/4/h/7v8fofegmtpAaihfJ0S+Dq3/XjePY2e/68XA0WHPcd18nzX/Pk8AvtvO+gWh/aSELVPujlJTj3Wccb5X8Rh8orkcWGVm7Z2k8Q7waTNLBeaFlmXjP8hjgObxsWvwHxicc/VAFfCwc64S/2Hfzf4AAAVQSURBVEEb3sX22sZ3gNYnQ3wIqHWhMbrAX/AfyIu62F/b7TZv8/7QY2/Cj/89s4t9tIqnJ/vswF+dH0+3GBgR9hoHh3qUvoo/ptDNY9bB692LT2CPO+fecM4904PtdRRPdzwITDA/HOSjwN3AP4GPm1k2UOVCmTfkJOCV0O3z8X+8mntBPgL8b+g4lwP/dn5c9hIg1zn3Z+e/Hn+Tnv0OROKC8nVi5Ot2dPgaO9nPGUCec645niU93Occ4FLn3P908/HK3VGSEnQA0n1mlgyMcc5tBr5sZr/Fn0FdAswNf6xz7oCZjQK+CdzTvAnn3C4zWwKcjv+P+gh8MXXoqWG3GwkNF+poex2EGr6NI4HUsLiazGwzMLqz/XXFOVdnZhuBfKCoB/vo9T7baKTl8zMW2OSc+2nbB3XnmJnZp/EFbLPJzn9l7MJj7cHvoMN4uuKce8fM3sQPB8lzzm00swfxvUgfA/4WFnfze6N5/N44fA9T87b2hG867PZBWmtCw9JkgFG+bvX6EiFfh+vwOHbyu5mB/yeit04DjjSzbOdcdQePaQq9BuXuKNIBiS+phI1/dc4tx3/VNqXtA83sNOALzrmbnJ/hIdx/AXPMT/u0AvhRVzvuYnud2YrvAQ2fhsjwYw/7Kg1YH+V9dEcprcchY2Zze3DMngCODWvtnl3eg+21G0/3Xw5/Aj4DbAyLLxn4Cv5r02Zz8eNFm20HTjazrLD9djrOT2QAU75uLZHydYevsZP9VIaek0bv3Ibv7f9ZB+srXMtJ58rdUaTCOv58ycxGht0fBSwK3a7Hn6wxBT+1U56ZpZvZGaH1Q82sCD/W77f4r3FeIfTVYkgSbf5LNzPrYntt921h23gMqAYuCW0rP7S/5p7PjvbXkRGhx4wKPXdRN/YRHk9X+wx/HW2113th+F6oLOBJM/ug+ZNkuvodHNoXMMI5tzasNYTtIjnsdmfbC4+73XjMLM3MSkLDPDrzZ/xX0A8DOOeqgH8AL7vWs4GchR8m0qz5eD9gfi7sLwLNc5m2/R20/T3H9QwFIh1QviZh8nV3j2O7+8Hn0jT8P1/gv53AzEZ2I3cbvvf408DHzOzjoef+PzO7Ht9b/6+wxyt3R1M0Bm6rRafhT4hwwD7gPuD3+K+SCkLrvwu8B1yA/2rnHXwPwYfwPR2/xyepJ/AnUtSHtleD/6CdAzTgT5AZA3w7tP6znW2vzb6/AryIH7c2LbSuGN/T+WP8WMPpoeUd7q+D1+9Cz/8KfmjCMWHrOtrH+PB4utpn+DFss+8L8Cfx/SJ0LJqf95nQ+vnAytDv5seh49zdY3ZBO6/14tD+HgBGhZZ19jttta0O4hmNv5DB/d14r/28zf1PAKe2WXZbO887BVgN7AGuDy07F/9euwMYSWjKMeAq/B/eu0Ov9cNBf8bU1CLVUL5OpHw9FngJ3+t+bBevsbPf9Yfw3xS+hj+RfD3whdBrbzd3478VacD3VBcC/8YX9Z8DPo+fz/vH+GEozc9R7o5is9DBkwRhZkOBq5xzN4buJ+F7UT7vnLs+0OC6YGYOOML5MYvSS2b2dedcr6ZyEpH+o3wt4ZS744OGgiSeL+KHBRSGLTuNsJPSYlHYV0/62qmXzCzZzC4FHgk6FhHpFuVrUe6OMyqsE8+v8WOx1pjZNuBp4D3n3KvBhtWx0MkcXwndvcL8tG/Sc9nAY865jV0+UkRigfK1gHJ3XNFQEBERERGR/9+OHdMAAAAACOrf2hKekMI5cKwBAGAgrAEAYCCsAQBgIKwBAGAgrAEAYCCsAQBgIKwBAGAgrAEAYCCsAQBgEMrSFeLcYueEAAAAAElFTkSuQmCC\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/pathfinder-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/pathfinder-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/pathfinder-mars.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }