{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 38 - Curiosity - 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 Curiosity aeroshell (MSL)\n", "https://en.wikipedia.org/wiki/Curiosity_(rover)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"MARS\") \n", "planet.h_skip = 126.0E3\n", "planet.h_trap = 2.0E3\n", "planet.loadAtmosphereModel('../atmdata/Mars/mars-gram-avg.dat', 0 , 1 , 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Curiosity', 3257.0, 146.0, 0.00, np.pi*4.5**2.0, 0.0, 1.125, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(125.0,0.0,0.0,6.08,0.0,-15.48,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (30*60.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZxcZZX/8c/pTtIJSUgC3WQPnQWSQNibEImyGkBRtl9QQUREaWTRkRkZHFEcBWZUkBFRkZZBB5FFZBRBDMu4sGkwIUBiCEkgnUBCQichC2Qh3X1+fzy3uqsr1Wttt6q/79freVXVrVv3nluVOn3y1HOfa+6OiIiIiIhkpqzQAYiIiIiIlAIV1iIiIiIiWaDCWkREREQkC1RYi4iIiIhkgQprEREREZEsUGEtIiIiIpIFfQodQLZUVlZ6dXV1XvfZ0NBAVVVVXveZKzqW+Cql49GxpDd//vz17l4ab0wX9TRn699QPOlY4knHkhsd5eySKayrq6uZN29eXvdZU1OT933mio4lvkrpeHQs6ZnZyqxsqIj0NGfr31A86VjiSceSGx3lbA0FERERERHJAhXWIiIiIiJZoMI6A7W1tYUOIWt0LPFVSsejY5FMldL7rmOJJx1LPBXLsZi7FzqGrKipqfG4jL0REekOM5vv7jWFjiOflLNFpFh1lLPVYy0iIiIikgUqrEVEMtTc3MyuXbsKHUZJc4emptBEROJKhbWISAaWLVtGeXk5BxxwQKFDKWnPPgt9+sAxxxQ6EhGR9qmwFhHJQJ8+4XIA6rHOjxI5LUhESpQKaxGRDPTt2xdQYZ1rZoWOQESkcyqsRUQykCisGxsbCxxJ76AeaxGJs7wX1mY2y8zmmll10rKzzGyFmW0ws5vNrE/Sc5VmdpuZfcfMbjKzvvmOWUSkPaU+FCQuOVs91iJSDPJaWJvZcGAQMD1p2TjgDGA2cAXwOeBLSS+7H7jD3a8ClgLfylvAIiKdKOWhIHHM2eqxFpE4y2th7e7rgAdTFu8LXOju8939TuBHwPEAZjYDmOzuc6N1HwK+YGYD8xWziEhHSrmwjlPOVo+1iBSDvA8FcffmlMdPuXvy4MTVwKro/gnAiqR1V0d3e9UVykQkvhJDQUp1jHXccrZ6rEUkzuJ48uKRwE+i+6OBDSnPbwVGpb6ooaGBmpqallZXV5fjMEVEWgvrpqYmvBtVX11dXUu+AipzE11e5CVnq8daRAqpqzm7T3tPFIKZTQIa3P3FaJED21NWqwB2+821qqqKefPm5ThCEZG2zIw+ffrQ2NjIrl276NevX5deV1tbS21tbWIb63MZY64UImerx1pECqGrOTs2PdZmVg5cDHwlafFqYFjSOgYMBNbkNzoRkfb1xin38p2zEz3WKqxFJM5iU1gDXwZudPedScseByYnPa4G3gPUNS0isVHqU+61I685W0NBRKQYFGIe60R6tKRlXwUWAAPMbIKZXWhmk9x9HvCmmU2LVj0FuMXd38tv1CIi7SvlmUHilrPVYy0icZbXMdZmNhg4L3p4gZn9kPBT4rUpqy5x9zui+x8HrjGzesJ8qlfnI1YRka4q1cI6TjlbPdYiUgzyWli7+1bg1qglXBe19l6zEvhsjkMTEemxUp1yL445Wz3WIhJncRpjLSJSlEq1xzpO1GMtIsVAhbWISIZUWOePeqxFJM5UWIuIZKg3TreXb+qxFpFioMJaRCRDvXS6vYJQj7WIxJkKaxGRDGkoSO7pAjEiUgxUWIuIZEiFde5pKIiIFAMV1iIiGSrV6fbiSD3WIhJnKqxFRDKkHuvcU4+1iBQDFdYiIhlSYZ0/6rEWkThTYS0ikiFNt5d76rEWkWKgwlpEJEOabi9/1GMtInGmwlpEJEMaCpJ76rEWkWKgwlpEJEMqrPNHPdYiEmcqrEVEMqTp9nJPF4gRkWKgwlpEJEPqsc49DQURkWKgwlpEJEMqrPNHPdYiEmcqrEVEMqTp9nJPPdYiUgxUWIuIZEjT7eWPeqxFJM765HuHZjYLuA74uLvXR8sGAN8FtgLDgKvcfUtnz4mIxEEpDwWJS85Wj7WIFIO89lib2XBgEDA95akfA0+7+1eBB4HbuviciEjBlepQkDjmbPVYi0ic5bWwdvd1hETbwsxGAecCD0WLngBOM7Pqjp7LS8AiIl1QqkNB4pSz1WMtIsUg70NB3L3Z2mbI44D17r4ter7RzFYAxwK7OniuPp9xi4i0p5SHgsQtZ6vHWkTiLO+FdRqjgQ0py7YCo4DGDp4TEYmFUi6s0yhIztYFYkSkGMShsHZge8qyCkLPB50816KhoYGampqWx7W1tdTW1mYxTBGR9Hoyxrquro66urrEw8rsR5UzBcnZGgoiIoXU1Zwdh8J6NeHM8WSDgTWEBN7ec21UVVUxb968nAQoItKRnoyxTi4kzWx9TgLLjYLmbPVYi0ghdDVnx2Ee6z8Bo6IpmjCzvsA44C+dPCciEgu9bChIQXK2eqxFpBjkvbC21rNgDMDd1wIPAydGyz8IPODuqzt6Lo8hi4h0qFSn24P45Wz1WItInOV1KIiZDQbOix5eYGY/dPcG4BLgu2Z2CDACSB5o19FzIiIFV6rT7cUpZ6vHWkSKQV4La3ffCtwateTlbwMXtfOadp8TEYmDUh0KEsecrR5rEYmzOIyxFhEpaqVaWMeJeqxFpBiosBYRyVBiKEgpjrGOG/VYi0icqbAWEcmQeqxzTxeIEZFioMJaRCRDKqxzT0NBRKQYqLAWEclQKU+3FzfqsRaROFNhLSKSoVKdbi9O1GMtIsVAhbWISIY0FCR/1GMtInGmwlpEJEMqrHNPPdYiUgxUWIuIZEjT7eVeWfTXqqmpsHGIiHREhbWISIbUY517/fqFW73FIhJnKqxFRDKkwjr3EoX1e+8VNg4RkY6osBYRyZCm28s9FdYiUgyyUlib2YezsR0RkWJUbNPtFWPOVmEtIsWgT3dWNrNTgK8BI2gtysuA4cCA7IYmIlIc4joUpJRydvR/Fxobobm59WRGEZE46VZhDfwC+E9gAdAcLTPg3GwGJSJSTOJaWFNCOdsMKipg587Qa92/f6EjEhHZXXcL69fc/abUhWa2JEvxiIgUnRhPt1dSObtfPxXWIhJv3f0x7RIzm51m+WnZCEZEpBjFuMe6pHK2xlmLSNx1t7CeCdxuZk1JrRm4NQexiYgUhRgX1iWVs1VYi0jcdXcoyLeAWmAekLj+lQGfyjQQM5sKXA4sB/YD6tz9BTMbAHwX2AoMA65y9y2Z7k9EJFuSh4K4Oxaf62+XVM6uqAi3O3dmY2siItnX3cL6D8CD7t4mrZnZT7MQyy+A0919tZmNAx4FpgI/Bua4+33RGe63AedkYX8iIllRVlZGWVkZzc3NNDU1tRTaMVBSOXvQoHD7zjvZ2JqISPZ1N/vfC/ybmf0xaZkBZwH/lGEsU4HB0f3twBAzG0U4e/2yaPkTwANmVu3u9RnuT0Qka/r27cvOnTvZtWtXnArrksrZg6O9bdFvliISU93N/l8n/OR3AeDRMiPMkZppkv41YSzgh4DzgC8AxwHr3X0bgLs3mtkK4FigPsP9iYhkTXJhPWBAbKaILqmcnSist27NdEsiIrnR3ZMXbwDGunu1u4+PWjWh9yNTlwI7gL8D77j7A8BoYEPKeluBUVnYn4hI1sR0yr2Sytl77hlu1WMtInHV3R7r1909XV/BiCzEMgC4h5CYvx/1cjjhJ8ZkFcBup943NDRQU1PT8ri2tpba2toshCUi0rnuzgxSV1dHXV1d4mFlbqIqrZytHmsRKZSu5uzuFtaXAn9NXmBmZwPfB+7o5rZS3QV8wt03mZkTEvaXCGeVJxsMrEl9cVVVFfPmzcswBBGRnuluYZ1cSJrZ+hyFVVI5W2OsRaRQupqzuzsUZC8z+3q00T3M7HbgdsJPgT1mZpXAIe6+KVp0HbAFWAmMiqZvwsz6AuOAv2SyPxGRbEsU1jEbClJSOTsxFEQ91iISV90qrN39VOBxM7sFeAHYHzgESHdlr+7YCOwws9FJyzYALwEPAydGyz4IPODuqzPcn4hIViXGWMfpIjGllrPVYy0icdfhUBAzO6adp1YCy4AbgX2BM4ArehqEuzeb2RnANWY2HxgOXOnuW8zsEuC7ZnYIYVygBk6LSOzE4eqLpZ6z99or3G5IPT1SRCQmOhtj/V3CXKWbgOY0z99BmLppJBkkaQB3fxG4OM3yt4GLMtm2iEiuxWQoSEnn7BHRKZfr1uVi6yIimeussP4WMM/d3+poJTPL9GdFEZGiFpOhICWds4cPD7dr1xY2DhGR9nRYWLv7I13ZiLv/OjvhiIgUp7KycMqKu3eyZu6Ues5O9FirsBaRuOrurCAiIiIFsc8+4fatt6A53UAXEZECU2EtIpIFheyp7i0qKsIJjE1NOoFRROJJhbWISBaZWaFDKGkjR4bbN94obBwiIulkXFib2WFmdpeZ/YuZVWcekoiI5Eqx5+zx48PtihWFjUNEJJ1s9Fh/Hvigu38PODpxxS0Rkd6kiIaCFHXOVmEtInGWjcL6cuCI6P69QN8sbFNERHKjqHP2hAnh9rXXChuHiEg63S6szexDZvbJ6P50YL/E5WrdvdnddbFZEem14jbGutRydqKwVo+1iMRRtwprM7seuB84E8DdnwPON7NzcxCbiIhkoBRzdmIoiHqsRSSOuttjfTQwBngqadkfgG9nLSIRkSIU0zHWJZezE4V1fb3mshaR+OluYf24u28Ckv+CFOWlcUVEciFmQ0FKLmcPGgSjRsHOnRoOIiLx093C+jUz+yow2cw+YWYPAJcB38l+aCIikqGSzNkHHxxuX3qpsHGIiKTqVmHt7vcCfwTKgfOBTcBJ7v6jHMQmIlI0mpqaACgri891t0o1Z6uwFpG46tPdF7j734C/JR6bWbmZneju/5fVyEREish7770HQEVFRYEjaasUc7YKaxGJqw4LazP7HTCok21UAuuBok3SIiKZ2rlzJwD9+vUrWAy9JWersBaRuOqsx/qlaJ1XgGbgImABMC9pnfHAXjmJTkSkSMSkx7pX5OzJk6FvX3j1VXjnnXBCo4hIHHRWWN8AWHRWOWY20d2vSV3JzB7MRXAiIsUiUVgXsseaXpKz+/WDqVNDj/WiRTBjRqEjEhEJOjzLxt03JxJ0ZLiZlSevY2bHAdNzEJuISNGIw1CQ3pSzDz003M6b1/F6IiL51N2TF38NvGhmTwNbganAycBuPSI9ZWYzgPcDy4GngXeB70b7GwZcVWyX4BWR0heToSCpSjZnH3003HknPPMMXH55trcuItIz3Sqs3f1xMzsJuBQ4mDB102x3z8rPimb2GWCiu38tadnPgDnufp+ZnQLcBpyTjf2JiGSDu7cU1n379i1wNK1KOWfPnBlun30221sWEem5nky3twZITqJlZnaEu8/PJBAzOxb4JDAradko4FzCBQ0AngAeMLNqd6/PZH8iItnS1NSEu1NeXk55eXnnL8ijUs3ZBxwAQ4bAqlXwxhswZkw2ty4i0jPdKqzN7I40i6uAJuCMDGO5CXgWuMXMJgHXAvsC6919G4C7N5rZCuBYoD7D/YmIZMWOHTuA2A0DKemcXVYG73sfzJkDTz4J556bza2LiPRMdy8Rth+wCliZ1CqATHs+JgOHAT9198sJ86vOAcYCG1JW3wqMSt1GQ0MDNTU1La2uri6TkEREumzr1q0A7Lnnnl1+TV1dXUu+IswtnQslnbNPOSXc/uIX3X6piEi3dDVnd3coyOfc/ZXkBWY2Bbig2xG2dSCw0d0T0/3/EPgG4MD2lHUrgF2pG6iqqmKeTg8XkQLYvHkz0L3Cura2ltraWgDMbH1OAivxnH3eeXDVVfDoo/Dyy2EKPhGRXOhqzu5Wj3Vqgo7sBD7Xreh214ekIt/dtwPLgHLCWeXJBgNrMtyfiEjWbNkSJr3oTmGdD6Wes/feGz7zGXCHb387F3sQEemebhXWZrbCzF5LaiuBl4HfZBjHS8AQMxuetKyRkIxHmdmAaP99gXHAXzLcn4hI1sS1sO4NOftf/xXKy+GXv4QVK3K1FxGRrunuUJDHgbsJP/dBuGTuOndfmkkQ7r7EzP4AzAZ+ZGZDgX7AL4EPAScCDwMfBB5w99WZ7E9EJJt6MhQkT0o+Z48fH05c/MUv4Oqr4e67c7UnEZHOdbewvsvdn0xeYGaTzWySuy/PMJbzgZvNrB+hh+Oc6IzyS4DvmtkhwAigNsP9iIhkVYwL616Rs6+9Fu6/H+65By6+GI49Ntd7FBFJr7uF9SnAkynL1gF/Ipwh3mPuvp4wJ2rq8reBizLZtohILq1btw6A4cOHd7Jm3vWKnL3vvvDVr8I118Cll8Jzz8HAgfmMQEQk6NIYazO73Mw2AVeZWVNyI0yttC6nUYqIxNiaNeHcvFGjdptVriB6Y86+8kqYPBkWL4YLLoCmpkJHJCK9UZcKa3f/IXAM8H1gQkob4+6n5CxCEZGYi1th3Rtzdv/+8JvfwODB8Otfw0UXQXNzoaMSkd6my0NB3P0lM/uKu+82H6mZjXP3VdkNTUSkOMStsIbembOnToWHHw4XjvnZz8JsIbfeCn26O+hRRKSHOkw3ZvZl4BF3Xxwt+lczK09dDTgeOC774YmIxF99fT0AY8aMKWgcytlwzDGhuD71VLj9dlizBu69N/Rki4jkWmf/jx8PDE16fDgwEFhL6/RNRjjzW0Sk19m0aRNr165lwIABjBs3rtDhKGcDJ5wATzwBp58OjzwCM2fCr34FU6YUOjIRKXUdFtbuflnKouuBRe7+XvJCMzsg24GJiBSDl19+GYApU6ZQVtata25lnXJ2q5kz4W9/Cz3XCxfCEUfAD34AF14IZoWOTkRKVXf/ChyUJkFPBiZmLyQRkeKxeHEYdTF16tQCR5JWr87ZkybB3/8O550H27bB5z4HH/oQvPZaoSMTkVLVpVM6zOwDhJ8PZ5lZ6kVj9wJ+ADyU5dhERGJv7ty5ABx66KEFjqSVcnarPfcMV2U8+WT44hfh0UfhwAPhG9+Af/kX6Nu30BGKSCnp6rnSm4A7CeP3ZqY89x5wezaDEhEpFs888wwA73//+wscSRvK2SnOOw9OOgn++Z/hl7+Ef/u3MHPIf/wHnHWWhoeISHZ0qbB294Vmdhxwtrv3uoQsIpLOxo0bWbx4MRUVFRx++OGFDqeFcnZ6++wDd90F558Pl10GS5fC7Nlw1FHwrW/BrFkqsEUkM10eY+3um9tL0GYWq64aEZF8eOSRR4DQW11RUVHgaNpSzm7fSSeFKzT++McwfDjMnRuGihxxBNx3HzQ2FjpCESlWnc1j/TtgUCfbMOBAYJ9sBSUiUgx++9vfAnD66acXOJJAObvr+vaFSy6BT30KbrkFbr4ZFiyAT3wCJkwIz33601BVVehIRaSYdDYU5BXCWL3VtM6BmsqAj2QzKBGRuNu6dStz5swB4lNYo5zdbYMGhfHWV1wB//M/cMMN8OqrcOWVcPXVYfz1RRfBccdBgWdTFJEi0Flh/R1gq7vv7GglM+sVZ5eLiCTcc889vPvuu8ycOTMOF4ZJUM7uof794eKLw5R8v/891NWFi8vce29oY8bAOefAuefCIYdoLLaIpNfh/7/dfX1nCTpyZpbiERGJPXfn1ltvBeDzn/98gaNppZydufJyOO20cFn0+nq45hoYNw7eeCP0Zh92WJiu77rrQs+2iEiyjH7YMrNxZvYd4OYsxSMiEntz5szhhRdeoLKyktmzZxc6nC5Tzu6ecePgm9+EFSvgqafg0kuhshJefhm+/vVwAZoZM8L47Pr6QkcrInHQ7cLazMrN7EwzewR4FfgkYUyfiEjJc3euueYaAK666ir69+9f4Ig6ppydubIyeP/74Uc/gjVrwhCR886DgQPDjCJf+hKMHw8HHwxf+xo89xw0Nxc6ahEphC4X1mY23syuB14H7gUmAR8FxgJn5yY8EZF4ueeee5g3bx7Dhw/n0ksvLXQ47VLOzo2+fcNl0X/xC1i3Du65B84+O5wEuXAhXH99mBd79Gi48MLw/FtvFTpqEcmXTgtrMzvbzB4DlgH/D7gJGAP83t3nePBsNoIxs35m9mJ0YQPMbICZ3WJm/2Fmt5rZntnYj4hIT2zcuJEvfelLAFx//fXsscceBY5od/nM2dH+em3eHjgwTM/3q1/B+vXhcumXXQZjx8LateHKjueeG+bKPvTQMNPIo4/Ctm2FjlxEcqUrV17cTrgE7rPA59x9KYCZtTeVUyauBKqTHv8YmOPu95nZKcBtwDk52K+ISIfcncsvv5yGhgaOOeYYLrzwwkKH1J585mxQ3gagoiJceOakk8K82C+9BI89Bo8/HsZnv/hiaDfeGHq9jzgiDC95//th5swwdltEip+5dy3Xmtk44CJCAn0QONbdvxA9Z97VDbW//aOBKcA1wAXAUmAFMMzdt5lZH2AzcKC716e+vqamxufNm5dJCCIi7aqrq+Piiy9m4MCBzJ8/n8mTJ2dt22Y2391rsrZBcp+zo+30OG/3ppy9Ywc880wosh9/PFyIJvXdnzIlFNkzZoShJFOnhhlKRCR+OsrZ3bmk+Sp3/zrwmWjRGDO71syqCck7kwAHArPd/Y6kxccB6919W7T/RkLCPjaTfYmIdNeTTz7JF7/4RQBuu+22rBbVuZLLnA3K293Rvz+ceCJ8+9swfz68/Tb84Q/hAjTHHhueX7IEbr89zKN90EEwdCgcfzx85Svwm9+EkyZFJP66MhSkjShR/hr4tZlNAq4gJOm6DOK4Cvh2yrLRwIaUZVuBUek20NDQQE1N638eamtrqa2tzSAkERFYuHAhp512Gjt37uSyyy7jk5/8ZFa2W1dXR11dS9rM2UCAHOVsyDBv9+acPWQInHJKaADvvRd6sZ9+OswyMncurFoFf/5zaAljxsD06aFH+6ijwnCSQZ1dwF5EsqKrObvLQ0E6Ymbnu/udPXztyUCFu/8uelxP+EmxBjjb3Y9KWvd54G53vzF1O73pZ0URyY9FixYxa9Ys1q5dy5lnnsn9999PeQ5+n8/FUJBO9tfjnB29PuO8rZzdsbVrw7R9iUL773+HLVvarlNWFi5Wkyi0a2rggAOgX7/CxCzSW3SUs7vdY51OJgka+DJwmLVeH3YYYTzgTdH9ZIMB/SAmIjk3b948Tj75ZDZu3MgJJ5zAL3/5y5wU1YWQYc4G5e2cGzEiXAHytNPC4+ZmeOWV1kL7uefCCZILF4Z2++1hvT59wvjsQw4J7dBDw21VVeGORaQ3yUqPdUYBmFUBFUmL/gr8MzAXWAxUuft2M+sLvANMcPfVqdtR74eIZMsDDzzA+eefz7Zt2/jIRz7C/fffn9MLweS7xzpT2cjbytmZ274dnn++tWd7wQJYtmz3EyMBRo5sLbanTQs921OmQAxnjBSJvZz3WGfC3RuSH5tZE9Dg7qvM7GHgROBh4IPAA+mKahGRbGhqauKb3/wm1157LQDnn38+t99+O3379i1wZPGivB0PAwaEqfpmzmxd9u67oQc7Mb3fiy+Gnu033wxtzpzWdc2gujr0cB9wQGubOhX2LNnZx0Vyq+CFdScuAb5rZocAI4DecWaLiORdfX09559/Pk899RRlZWXccMMNXHHFFSQNd5CuUd4uoIEDw5R9M2a0LmtuhhUrWovsxYtDW7o0LF+xIlymPdmYMa1F9v77w377hduxY8PYbhFJr+BDQbJFPyuKSE+4Oz//+c/5p3/6J7Zu3cqIESO48847mTVrVt5iKLahINmgnF14u3bBq6+2FtqJtmQJ7NyZ/jUVFTBxYttie7/9Qhs5MvSCi5S6WA8FEREplIULF3LZZZfx1FNPAXDWWWdx2223UanL4Ekv0LdvGGc9ZQqcdVbr8qam0Iu9eDG8/HIYt71sWejhXru2tQBPtcceYWjJ+PGhJd8fPz7MzS1S6lRYi0ivs27dOq677jpuvfVWmpqa2Gefffje977HJz/5SQ39kF6vvBwmTQotMStJwtatsHx5KLKTC+6lS2HjxvaLbgiFdbqCO7FMJ1JKKVBhLSK9xqZNm7jxxhv5/ve/z7vvvktZWRmXXXYZ1113HUPVnSbSqcGD4bDDQku1eXPrmO3kVl8fbjdtCjOXLFiQftuVlTBuXBjHne525Ehd5l3iT4W1iJS8NWvW8IMf/ICf/OQnbN68GYCPfvSjXHfddRx88MEFjk6kNAwZEubNPvTQ3Z9zh4aG9gvvlSth/frQnn8+/fbLy2H06I6L72HDNM5bCkuFtYiUrEWLFnHTTTdx1113sWvXLgCOP/54rrvuOo4++ugCRyfSe5jBPvuEdtRRuz/f1ATr1sHrr4fLuae7Xbcu3F+1qv397LFHKLJHjQpF+OjRbe+PHh0uvqMZNCVXVFiLSEnZtm0b999/P3V1dTz77LMAlJWVMXv2bK688kqmT59e4AhFJFV5eSiAR41KX3gD7NgBq1e3X3ivWhXGgL/ySmjtSRT5qUV3agGu3m/pCRXWIlL03J25c+dy1113cdddd7UM9xg8eDCf+tSnuOKKK5g0aVKBoxSRTPTvH6b6mzix/XU2bQrF9+rVsGZN6/3kx+vWtbb2hp0k9pdcaKdrI0dCv37ZP1YpXiqsRaRoLVq0iLvvvpt7772XFStWtCw/6qijqK2t5WMf+xiDBg0qYIQikk9Dh4Z24IHtr9PYGIrqdAV48v0tW8I836++2vE+E73fiTZmTBiOkrgdOzZcJVN6BxXWIlI03J0XX3yRBx98kF//+tcsWrSo5bmRI0fy8Y9/nAsuuIBDDjmkgFGKSJz16dNaBHfknXfa9nintjfeCPN6v/VWaO3NdgKw996tRXZq0T12bIiloiK7xymFocJaRGJt165dPPnkkzz44IP87ne/Y+XKlS3PDRs2jNmzZ3POOedwzDHHUK65uEQkSwYNgsmTQ2tPau93ouB+/fW2txs2hPbCC+1va/jw3Qvu5DZypE66LAYqrEUkdlatWsVjjz3GY489xuOPP86mTZtanhs+fDinnXYap59+OrNmzaKfBjiKSIF0pfe7uTn0aL/+emtLFN2Jtgk5mUEAACAASURBVGZN67jv+fPTb6esLMxoklpwJxfjI0Zoru9CU2EtIgX37rvv8uc//7mlmF6yZEmb5w844ICWYnr69OmUlZUVKFIRke5JFMQjRsCRR6Zfp7ExDCtJV3Qn2tq1oQBfswbmzk2/nT59wuwm7fV6jx0LVVWa7SSXVFiLSN5t2bKFZ599lieffJInn3yS5557rmWeaQizeZx44omcdNJJnHTSSUzsaBoAEZEi16dP6HkeM6b9dd57D958M33RnWgNDZ3P9V1R0fGQk7FjwwmgKr57RoW1iOTcxo0beeqpp1oK6eeff57m5uaW58vKyjjqqKM46aSTOPnkk5k+fTp9NZhQRKRFv36w776htWfHjt3HeKe2t9/ufLaTgQM7L74HD87+MZYCFdYiklW7du1i4cKFzJ07l7/97W/MnTuXV1Ku1tCnTx+mT5/OMcccw7HHHsvMmTMZMmRIgSIWESkN/fvDpEmhteedd9ovuhPtnXc6v9DOkCEdF95jxvTOaQZVWItIj7k7b7zxBnPnzm0ppOfPn8/27dvbrFdRUcGMGTM45phjOOaYY3jf+97HwIEDCxS1iEjvNWgQTJkSWjrusHlzxydbvv56WGfzZkia9XQ3lZWhyB4/HiZMCBf3mTAhtH33Lc1ZTlRYi0iXuDsrVqzg+eefZ8GCBTz//PM8//zzvPXWW7utO2nSJGbMmMFRRx3FjBkzOPjggzV7h4hIETBrvdDOQQelX8c9TB/YUa/36tWwfn1o6eb4LiuDceN2L7gTj4cNy+1x5ooKaxHZTVNTE0uXLm0pnhcsWMCCBQvaTHuXMHToUKZPn95SSE+fPp3KysoCRC0iIvlgFnqjKyvhsMPSr9PcHKYPXLUKVqyA114L7dVXw+3rr0N9fWh//OPurx86NAxpSfSuJ9qkSfG+mE5sCmszOwv4HrAncBfwL+7eaGaVwPXAJqAvcJW772p/SyLSVe7O2rVrWbhwYUtbtGgRixcv3m04B8A+++zD4Ycf3tIOO+wwxo8fj+n08V5HOVtEOlJWFi5qM3IkHHXU7s/v3BmK7kShnVp4b9oE8+aFlrrdCRPaFtsHHwzTpsVjTHcsCmszGwecAcwGDgRuBV4HbgTuB77i7nPN7PPAt4B/K1SsIsVqy5YtLFq0qKV4ThTSGzduTLv+2LFjOeKIIzjssMNaCumRI0eqiBblbBHJWEUF7LdfaKncw9SBy5bBkiVt22uvwfLloT38cOtryspCkX3ooW1bVVX+jgliUlgD+wIXunsjMN/MpgHHm9nTwGR3T0yF/hDwipld5+7vFipYkbhyd9566y2WLFnCkiVLePnll1tuV7UzsemQIUM46KCDWtq0adOYNm0aw4p1gJvkg3K2iOSMGeyzT2gzZ7Z9bseOUFQnCu3Fi+Gll1rvL14Md9/duv6++8LRR4ftHH10GDfeJ4fVbywKa3d/KmXRamAwcAKwImm91VFvWQ3wl7wFKBIzjY2NrFixYrcCesmSJbz99ttpX1NRUcEBBxzAtGnTWgrogw46iNGjR6sXWrpFOVtECqV//zDsY9q0tsu3b4d//ANeeKG1vfgirFwZ2j33hPX22gvOOgs++1mYMSP78cWisE7jSOAGoBbYkPLcVmBU6gsaGhqoqalpeVxbW0ttbW0uYxTJKXdn3bp1LFu2rE1bsmQJy5Yt47333kv7uiFDhjB16lSmTJnSpk2cOJE+ufxvunRLXV0ddXV1iYfFfrancraIFNSAAVBTE1pCU1PowX7mGXj2WXj66XAi5e23h3bppfDDH3btKpNdzdnm7pkdSZaZ2STgMne/wsx+CFS5+8eTnt8I1Lr7r5NfV1NT4/NSR7iLxJy709DQsFvxvGzZMpYvX84777zT7mvHjh3bpoBO3B8+fLh6oIuMmc1395rO14wf5WwRKSaLFsGdd8LNN4fLxD/0EHzkI93bRkc5O1bdV2ZWDlwMfCVatBrYP+l5AwYCa/IfnUjPJMY9v/baa7sVzsuWLWPLli3tvnbYsGHst99+bdrkyZOZPHkygwYNyuNRiOxOOVtEis20aXDOOWEc9urVsHZtdrcfq8Ia+DJwo7vvjB4/Dnw+6flq4D1A3RwSK9u3b6e+vp7XXnstbdu2bVu7rx0yZMhuxfN+++3HpEmT2HvvvfN4FCLdppwtIrHnDkuXwu9/Dw88EIaFAOy/P8yend19xaawNrOvEpLvADObABwHPAm8aWbT3H0RcApwi7unH1wqkiOJ+Z7bK5zXrOm4Q26vvfZi/PjxTJo0abfiubKyUkM3pOgoZ4tIXG3cGE5cfP751rHVyRcJHjQILroIvvENGDIku/uORWFtZl8Drk1ZvMTd7zCzjwPXmFk9MAi4Ot/xSe+wadMm6uvrW9qKFSvaFM87duxo97V9+vShurqaCRMm7NbGjx/P0KFD83gkIrmlnC0icbB1a5jretmycJJiYjaQdLPL7rMPHH88nHkmfPjDMHhwbmKKRWHt7tcB17Xz3Ergs/mNSEqNu+9WONfX17Ny5cqW+5s3b+5wG5WVlWkL5wkTJjB69GjNuCG9hnK2iOTDzp1hHPSqVeES6InLoC9dGtqbb6Z/3YAB4WqMhxwSptR7//vDpdDz8eOwKgEpCe7O22+/vVvhnFxAd3SSIMDAgQOprq6murqafffdd7de5z333DNPRyMiIlK6du2CdevCiYNvvhluE/fffBPeeCMU0evWdbydiopQMO+/P0yeHArpQw8NV3MsL8/PsaRSYS1Fwd3ZuHFju4VzfX19h1PTAQwaNKilcE4Uz8mP9957b411FhER6YFdu2D9+jCWuaEh3KYWz4nb9eu7ts3ychg9GsaODW3cuND23z+0MWMKV0C3R4W1xIK7s2HDhg4L53ff7fiKyIMHD25TKKcW0HvttZcKZxERkS5obIQNG9oWyh3dtnPR37TKysKY55EjYcSI1ttEGzMmFNIjR8avcO6MCmvJi2wUznvuuWfawjlRPA8bNkyFs4iISIrGxlD4rl8fiuWOWqJQ3rgxTFPXVWVlUFkJVVWhaK6qguHDdy+eR44M6xVbwdxVKqwlK7JVOI8fP363gjmxTDNriIhIb7dtW/qCuKOiedOm7u/HLBTAiSK5s9u99grFdW+nwlq6JFeFc3JT4SwiIr2Be5gqbuPG1vb22x3fTxTJHcz82i6zUPjuvXfnLVEo77136fYq55IKawFUOIuIiHTXe++1Fr6dFcbJ999+G5qaerbPiorWIriysmvF8tChKpLzRYV1L6HCWUREZHfu8M476QvhzorkTiaj6tCgQaEXea+9YNiw1vupj5Pv77037LFHfuZjlp5RYV1CmpubWbNmDa+++irLly/f7bazeZxVOIuISLFqaoLNm3cvjLvSGht7ts/y8q4XxcmPhw2Dvn2ze/wSDyqsi8yuXbtYuXJl2uK5s8tuDx48mAkTJqhwFhGR2Nq1q2fF8aZN3ZvFItkee7QtgLtaIA8erN5jaUuFdQy5O6tXr2bJkiUsWbKEV155haVLl7J8+XJWrlxJUwcDs6qqqpg0aRITJ05k4sSJLfcnTZpEZWWlpqMTEZG8ee+9cMJdQ0OYtSK5JS/bsKG1QN66tef7GzJk9wK5szZsWBi3LJINKqwLaPv27SxbtqxNAZ24bW+8s5kxduzY3YrmRCGty26LiEguNDeHXuHOiuTkx52MQEyrrKx7RXHi/tCh0EdVjRSY/gnmwfbt21m8eDELFy5k4cKFLF68mCVLlrBy5Uq8nd+tqqqqmDx5MlOmTGHKlCnsv//+TJo0ifHjx9O/f/88H4GIiJSixsZQBCcuPZ18m3z/rbdCr3Jzc/e2X14eZq5ItKqq3R8nZrdIFMiDB2s+ZCleKqyzqKmpiddee62lgE605cuX05wmG5WXlzNp0qQ2BfTkyZOZPHkye++9dwGOQERESsHWrbB6dWhr1qQvlteuDb3K3RmXPHRo54Vy8uMhQ1QkS++iwrqHGhsbWbx4MfPmzWtpixYtYvv27butW15eztSpUznooIM46KCDmDZtGlOnTmXChAn01WnBIiLSRY2NoShOFM3Jbc2a1vtdHads1nrp6REj2r9NFMz6kyXSMRXWXeDuvPLKKzz33HMtRfQLL7yQtogePXp0SwGdaFOnTqVCZ0aIiEgH3MOwjPp6WLEi3K5atXvPc1eGY/TvD6NHhzZqVCiO2yuYNS5ZJHv0dUqjsbGRF154gaeeeoonn3ySp59+mvXr1++23oQJE6ipqaGmpoYjjjiCQw89lL322qsAEYuISDHYsQOWL4elS1uL5xUrWu9v29b5NvbZp7VoTi6ekx8PG6Zp4EQKQYU1oUd64cKFPProozzxxBM8++yzvJNyOaURI0bwvve9jyOPPLKlkFYRLSIiqdzDyX5LlsArr4TbxP0VKzoe0zx0KFRXw/jxoY0b17ZwHjkS+vXL26GISDcVRWFtZgOA7wJbgWHAVe7eg0l8Wm3YsIEnnniCOXPm8Nhjj7FmzZo2z0+aNIkPfOADLW3ixImaA1pEpAtykbPjbNEieP55ePFFeOmlcNvQkH7d8nKYMAH23z/cjh/fWkhXV4fCWkSKV1EU1sCPgTnufp+ZnQLcBpyTyQavvfZabr755pbHI0aM4JRTTuHkk0/m2GOPZeTIkZlFLCLSe2U9Z8fZxRfDs8+2XbbnnjB1KkyZApMnh9spU2DiRPU4i5Sy2BfWZjYKOBe4LFr0BPCAmVW7e31Pt3vqqaeycOFCTj75ZE455RQOOugg9UiLiGQoVzk7zk45JQzVOPhgOOSQ0MaN0xhnkd4o9oU1cByw3t23Abh7o5mtAI4F6nu60VmzZjFr1qysBCgiIi2OIwc5O86+/vVCRyAicVEM07aPBjakLNsKjCpALCIi0jHlbBHptYqhx9qB1AmjK4BdyQsaGhqoqalpeVxbW0ttbW3uoxMR6YG6ujrq6uoSDysLGUuWKWeLSMnpas4uhsJ6NeGs8mSDgTbTeFRVVTFv3ry8BSUikonkQtLMdp8ov3gpZ4tIyelqzi6GoSB/AkZF0zdhZn2BccBfChqViIiko5wtIr1W7Atrd18LPAycGC36IPCAu68uXFRB0k8CRU/HEl+ldDw6ltKX65xdSu+7jiWedCzxVCzHEvvCOnIJcLqZXQ18GIjFQLxi+ZC7QscSX6V0PDqWXiNnObuU3ncdSzzpWOKpWI6lGMZY4+5vAxcVOg4REemccraI9Fbm7oWOISvMrAFYmefdVgKlctKRjiW+Sul4dCzp7evuVVnaVlHIIGfr31A86VjiSceSG+3m7JIprEVERERECqlYxliLiIiIiMSaCmsRERERkSxQYd1DZjbLzOaaWXWhY8mEmZ1lZivMbIOZ3WxmRXFCa3vM7Fgze8XM3jaz7xc6nkyZWT8ze9HMjit0LJkys+vMzKP2YqHjyZSZzTCzL5vZGWZWSldOLDgzG2Bmt5jZf5jZrWa2ZzvrfcjMfmhmPzaz2fmOszNdza9mtp+ZNSV9Pw7Ld6xd1ZXvcdw/F4Do77entA+nWW+omW1JWufMQsSbKl0N0tXvTbRubD6jdo6ly7VJ3L4/RV1EFYqZDQcGAdMLHUsmzGwccAYwGzgQuBV4HbixkHH1lJkNAWYCM4Cjgd+a2cPu/kRhI8vIlUB1oYPIlJkNJJx4MitalO8TjbPKzD4DTHT3rxU6lhL1Y2COu99nZqcAtwHnJK9gZvsD1wM17t5sZk+a2RJ3X1SAeHfTzfz6OeBDQDOwy90X5CvO7ujK9zjunwuAmU0F/gp8GdgVLb4d+GOa1S8CPg1sjR7/KecBdqKDGqTT7030+th8RumOpQe1Say+Pzp5sYfMrAxoAsa7e32Bw+kRM/sA8Fd3b4wefxc40N1PLWxkPRNd6W2HR/+ozWw+8GV3L3gi7AkzOxqYAlwDXODufy5sRD1nZl+K7v7E3XcUNJgMmdmxwNeBWa4EmnVmNgpYAQxz921RT9VmQm6qT1rvJ8BGd/9q9PhqYJK7f6YAYe+mq/k1Kiz+G/i8u7+R/0i7rivf47h/LhAKN3dflfR4X+C77v7xlPUqCBc7utTdl+U5zA6l1iBd/d5Er43VZ5TmWLpcm8Tx+6OhID3k7s2FjiFT7v5U4h9uZDWwqr31487dtycV1QOBl4E/FzSoHorin+3udxQ6lkyZmQGfAL4HrDOz3XpQisxNhH9bt5jZHDObWeiASsxxwHp33wYQ5agVwLEp650QLU94JXptLHQjv54DHA+8bmY/i777sdON73GsPxeA5KI6cibw2zSrfgQ4ClhqZr+3GA35SlODHEfXvjcQs88o9Vi6WZvE7vujwlqSHQn8pNBBZCr63+4fgL7AHgUOp6euAr5d6CCywYMZhJ+QbwF+aWYfKnBYPWJmk4HDgJ+6++XA/wFzol4TyY7RwIaUZVuBUZ2sl26dOEmbX939+8CewGmEy7//d57j6pJufI+L7XOBcHXQ36cudPcHgCGE4rQa+F3UuxpHXf3epFs37p9Ru7VJHL8/cf0HInlmZpOABncv+pPKgOXAz4ETgRsKG0r3mdnJwDx3f6vQsWSTu78djUn+T+BLna0fUwcSfkJ9KXr8Q8CA0wsXUslxYHvKsgpax8K2t166dWKhs/zq7k3u/hBwMnCWmY3Ma4Dd0IXvcdF8LgBRL/QOd9+S7vnoPxRPEnp5pxB6sOOoq9+bdOvG9jPqSm0St++PCmvBzMqBi4GvFDqWbHD3N6MhFF8m/c9gcfdl4A4zW29m64GxwINm9q8FjitbfgCMK3QQPdSHpJO+3X07sAzYu2ARlZ7VwLCUZYOBNZ2sl26dgutOfnX3xYRfQYrh+9He97goPpckpwG/62wld18H3Ed8P5uufm/SrRvLz6i7tUlcvj8qrAVCIXeju+8sdCBZNo+QQIrNucChSW0N4aznoh+mE2kGni90ED30EjAkZehHI7C4QPGUoj8Bo6KTkTGzvoQ/lH9JWe9xYHLS4/2JwYwNaXQ3v74LLMlhPNnS3ve4WD6XhC4V1pEm4IUcxpKJrn5voHg+o57UJgX//qiw7qHoRA4IPwMXLTP7KrAAGGBmE8zswuinl6JjZnukzF/5YeDmQsXTU+7e4O5vJBohmTe091Nl3JnZKDM718zKou/NvwBXFzqunnD3JYTx+7MhzHEL9CPN+EzpGXdfS5iJ4cRo0QeBB9x9tYV5wxMFwU8IP/0m8vEJhF7U2Ggnv+5nZt8ysxHROheY2T7R/Q8Af3H3zQUMO62OvsfF9rkkRCe69U8edpd8LGY228wmRPenEIaBvVKYaNtKrUE6+t5E68f2M0pXT7VXm1gQ7++Pu6t1sxF+NrmEME7pm0BVoWPq4XF8LTqG5PZyoePK4HimAw3AU4Sx1WcXOqYsHVc9cFyh48gg/gMJZ3QvIiT0KYWOKcPjqQR+CVwB/BdwQKFjKrVG+Jn6p4TC7RZgULR8PmG2nMR65xFmabkJOL3QcaccQ9r8CgwgzMgwI1rvQWA9cA9wfqHj7uB42v0eF9PnknJMs4EvpSxrORbCORSbgP8FvgiUFTrmKK60NUh735s4f0bpjqWj2qQYvj+ax1pEREREJAs0FEREREREJAtUWIuIiIiIZIEKaxERERGRLFBhLSIiIiKSBSqsRURERESyQIW1iIiIiEgWqLAWaYeZnWhmz5rZcZ2sN83MfmVmX8/ivm/I1rZERArBzA42s/uzmRvb2c/DZvaJXO4jV8zsg2b2187+zkjx6FPoAEQ6Y2ZTgeuBjcCRwMHAp9z9rhzv+nXgkC6st4owMX95NnYaXYp2Vza2JSICYGZHAv8JHAM8RrjIxgigL3Cruz+ag90uB4aQpdzYgR8TLlyTdVE+PotwUZylwF+jpyqA/wdMd/f6DHbR1b8zUiRUWEsx+C3wMXd/0czKgJ/lY6fuvtTMNnRhvS1m9mYWd30s8GQWtycivZy7/93M7iVcMfEjieVmdirwGzP7srv/OMv73Jbl3AiAmX0UWJgoaN39kWzvI8HddwH3mdkXgD+7+9eS4niwJ9uMLsF9vLvf5+6vmNn6LIUrMaChIBJrZrYXsD/wLoC7NwPXEC5xmg/NWV6vK44H/pTF7YmIADSmLnD33wP/DvyXmY3IwT6zmqvNbAyhhzrfdnvvgCeA1d3ZiJlVAHcSLs2dkM2/H1JgKqwl1tx9I/AKMMfMZkTLVgLPmdlQM7vWzJ4zs9PMbK2ZvWxmhydeb2bnmNnVZvZnM/uPpOU1ZvbvZnafmf2vmQ2Jlvc3sx+Z2TVm9t+EnzG7xcz+NYrpY2Z2pJn9NtretVGM95lZVbR8s5l9I2UT/dx9Z7StfzGzfzazZ8zs291+A0VEOvffhCEhZwOY2UQz+5aZ3W5m/xcVs5hZRZTLrjazp8zs44kNmNmlZvZ1M/ubmV3e3o7S5WQzmx3lzAvMbJmZfdvMRpnZz6J8+pSZnRBt4lxgDHCVmZ0RnQvzhJl9OmkfF5nZ9WZ2p5n9j5kNNLMhZvYNM5tnZieY2Stm9pqZVffkDYvek+OiHm3M7Cwzu9nMvmdmc8xsUHQMN5vZQ9Hfmzej9/gA4Gwz++ekTQ6L/hZtNbOL2tlnmZl9x8xqo324mf0m8flITLi7mlqsG3AgsJLwv/o7gH2i5eWEJLWRMNZtKPA48A/CfxrfB3wtWncYodf7VGBP4N6k7T8D3BDd/y+gNro/HGgiJM/OYvw5oddnH8JYvD2i5RXAo8CDwCjCH4RdwDeAPYAjoseJ9RPjxxPH/ePo/mDgK4X+LNTU1Iq3ARcAb7Tz3DrgB1Fe/V+gPFp+N3B/dP8m4Kzo/oVAQ3T/HOC86P7+Ua4+KHr8c+Dfo/vt5eS9gE3A94HpwAeBG4FronWvBB5KitWB6uj+AOBl4ILo8cnA3UnrPgzcRhj6emb09+Ks6DjnANd38b37c/S34tvADYTx1h9Men49MC66/w/C36QK4GrgDeAo4LNR3v9zIt5o/Xrge9Fz5wCr2onhc8DD0f3h0ft8UqH/Xam1bRpjLbHn7v8ws4OA64BLgY+a2Qc9jLluALa4+wMAZvbvwNPABMIfkWYz+1K0qUcJCfxUYO+k5S9Hrx0EXBK9FndfZ2avdyPUSuAW4LPuvi3axk4zWwuscPc10X7WAn/xMP5wASHhVxJOgjyV8EcAYBtwvpktBW4l/KdCRCQXnFBsTgf2Bb5gZgAbgCEWHlxM6BQA+B/CUAgIuXZeSq4dCSxM2ccFpMnJ7r7RzDYDD7r7cwBmtgp418xGEnp4B6UN2n27mb2VtOjyaLsJtwO/ipZvIvy9+N9oH/OiOLvqTx6NsTaz+widOQlnu/sqM5tJKKgHRfl/DfCqu88F5kavTbfth6K/CX/vIKbDge3Rca8zs420HVIiMaDCWmLNzMqBsR5OUvliNDzjEaCO0AOQ6h/R7VBgHPALd787Wvb9aJtXEQrd76fs6whCQtzaw3D3B94PTAReTFqeOsawqeUJ9+YoySaGZe3tYfgL7r7CzC4AfghcAXwSeAsRkSwys6GEX9v+QcibG9Lkx+GEHtV+AO7eROgMIHrNN9392ehxm9cmSZuTI07bXPkGcBWwDPgbML6DQ0h+3X7AH5Mev0oY5rIPu+fiRno+Y8nzhP+AJGwys+8A9xNiT1TPqcfVmUbar83+ANxqZv0Jn8V2Qu+3xIjGWEvctYz7A3D3F4GvAlPaWb8fIYktB94EZic/aWZHRctPNbMBKcu3RA+n9jDWZ4GbgV9Gia9bzGwUSSfCmNlowowo+xNOZry7nZeKiGTi04RfyO4n5MeZUU8x0JIf1wPvAackLT/IzAaTkmvNrE/UUZGqvZyczk+Aeg/Tqu7sxrGsIvRwt+wC2EwY6pI17t4cdX4cGP3H5BHgW+4+L5v7SdnnQ8BDhF9WzwGOdffNudqf9IwKaykGX4iKzoTRhLFxCXuZWb/o/nGE8XWbgHuAM83sh2Z2XNSb8C4hAQ4EHjazk6MTRaYQivEFwLVm1s/MKgljm/cxs3Iz+4qZnddOjBa1awh/fG5K81zq+qmPP0IYD5iwH3CGu28BvpDmNSIi3bFbT6iZHQN8E7jY3RsIvcPrgEfM7KPRCYofiHqo7yXMHnKGmZ0IfNrdtxJy7RfN7JtmdiyhgyExjC45/7WXkxOSe48PA6rMbA9gJjDAzCZEz+0inOyX6GBJ3sdPgP9nZsOix0cBP3X3RkLNkzYXm9mMKK6Kdt67MlJqpug/H6cSfqUcRhhiOBkYC+yRdGJkaq/4e2nibxtUmvEi0WexJ/B/hF757vyHQ/JEQ0GkGIwF/mFmvyf8TGbA55OebwSuNrNNhN7dLwC4++Nm9kXCz4lnEE7+WwRgZqcDPyL8obgDuNLd3cw+RpgK6R/AfYSf9A4lXFDhY4QhJm0uTGNm0wkn5Ywj/DT4W+CbUTyPAjOAajObSOhJGQF8zMyW0Np7cwEw0N3rUo79ZxYu7NBIOFlIRKTbojxyLqH4+ynQQDgBbgjwIXf/K4C7vxflx58Sct2DhLHVAP9E6JS4k/ALXWImjv8m5L5LCEPWLnH3t6LCcQYw0czuai8nm9lswrjiz5nZEndfTTiR8kZCUX0ncDpwEPAa8AtC7/r5FmaLOhD4sJk94u7/a2b7AQ+Y2TOEXz2vjnrWzwZGRMf3D8LFckaa2YFRnJ8gjBv/bdL71peQ+w8BRkeFtxPGNp8VbfMF4C/A3wnDW/5AyOkPE06YPMTMznD3xHbvI1ys552oEB8BzDazFYSeaKLXp16zYRXh5MwzgP5AmZk9QTiBMV9T0EonTJ+FFDMLl4H9ubtX52l/V7q7LjcuIlJiomEpFe4eywt0mdmVhBla6qPH/Qm/NnwncW6OFJ6GgkixSzfMIjc7Cj+ZPp+PfYmISP5EJ2ceHuOiuj/wNaA6aXEl8K6K6njRUBApbjrWyQAAIABJREFUWtH4tvMIP+2d7e7353iXr0Y/UYqISGnp4+63FjqI9rj7DguXVf9BdGL7GsKwla91/ErJNw0FERERERHJAg0FERERERHJAhXWIiIiIiJZoMJaRERERCQLSubkxcrKSq+uru7WaxoaGqiqqspNQEVI70dbej9a6b1oK9vvx/z589e7e696g3uSs6E4/y0WY8xQnHEr5vwpxrizFXNHObtkCuvq6mrmzevelURramq6/ZpSpvejLb0frfRetJXt98PMVmZtY0WiJzkbivPfYjHGDMUZt2LOn2KMO1sxd5SzNRRERERERCQLVFiLiIiIiGRBry6sa2trCx1CrOj9aEvvRyu9F23p/SicYnzvizFmKM64FXP+FGPc+Yi5ZC4QU1NT48U21kdEBMDM5rt7TaHjyCflbBEpVh3l7F7dYy0iIiIiki0qrEVEREREsqD3FtZ33AGHHw633FLoSEREpBNbX9jKvMPnseQzSwodiohIu3pvYb1+PSxYACt73fSxIiJFp3lbM+8seIdtr2wrdCgiIu3qvYX1gAHhdvv2wsYhIiKdsnIDwJtK44R7ESlNKqx37ChsHCIi0rnycKPCWkTiLO+FtZnNMrO5ZladtOwsM1thZhvM7GYz65P0XKWZ3WZm3zGzm8ysb1YCUY+1iEin4pKzrSz0WNOUja2JiORGXgtrMxsODAKmJy0bB5wBzAauAD4HfCnpZfcDd7j7VcBS4FtZCUaFtYhIh+KUs1uGgjSrx1pE4iuvhbW7rwMeTFm8L3Chu8939zuBHwHHA5jZDGCyu8+N1n0I+IKZDcw4GBXWIiIdilXO1lAQESkCeR8K4u7NKY+fcvfGpEWrgVXR/ROAFUnrro7uZn6FMhXWIiKdikvOTvRYayiIiMRZHE9ePBL4SXR/NLAh5fmtwKiM99K/f7jVyYsiIpnIS87WrCAiUgz6dL5K/pjZJKDB3V+MFjmQ2qVcAf+fvTuPj6uu9z/++mTfk6ZJmu4LUOgCZQmbbGUVhVtFEWRR0StRfuh1w6vgghdQFkHxsklQ9HJdQERlUZCiXMouLaXQQqGl+56kTZM0SbN9f398J8lksieTWTLv5+NxHjNz5sw5n/l2+plPvvM930NL6GsrKyspK+vqFCkvL6e8vLzvg6nHWkSiqKKigoqKio6HRdGMZbgimbNVWItINA02Z8dMYW1mycAXgG8Hrd4KzA7axoBsYFvo64uLi1m6dOngD6jCWkSiKLiQNLOqKIczZBHP2RpjLSJRNNicHUtDQa4CbnXO7Q9atxg4OOjxDKAZGEI27oMKaxGRkYhoztYYaxGJB9GYx9o67gatuwZYDmSa2Swz+5yZHeicWwpsN7P5gU3PBu5wzjWPOJCsLH/boMvjioj0JVZydsc81uqxFpFYFtGhIGaWC1waeHiZmd2J/ynx+pBNVzvn7g/cvxD4vpltwM+n+p2wBJOX52/37gXnoPO7Q0REIMZytoaCiEgcMOfGRpIqKytzQxqvB5CTA/v2QW0t5OaOTmAiIgMws2XOuZFPIxpHhpqzmyubeankJVLGp3Bi1YmjGJmISP/6y9mxNMY68goK/G1NTXTjEBGRfmmMtYjEAxXWoMJaRCTGabo9EYkHKqzBj7MWEZHYpTHWIhIHVFiDeqxFRGJc51CQ9v63ExGJJhXWoMJaRCTGaSiIiMSDxC6sx43zt1Vxd9EzEZHEEvi2UmEtIrEssQvr0lJ/u2NHdOMQEZF+BQ8FGSvTxIrI2JPYhfXEif52+/boxiEiIv0ys65rP2qctYjEqMQurNVjLSISNyzNV9btzaqsRSQ2JXZhrR5rEZG4kZTuv7La96uwFpHYpMIaVFiLiMSBpIxAYd2kwlpEYlNiF9bFxZCU5GcFaW6OdjQiItKPjh5rt18nL4pIbErswjo5GSZP9vc3bYpuLCIi0i8NBRGRWJfYhTXAAQf42/ffj24cIiLSr86hICqsRSRGqbA+8EB/u3ZtdOMQEZF+WXpgVhAV1iISo1RYq8daRCQudA4F0cmLIhKjVFirsBYRiQs6eVFEYp0K64MO8rerV0c3DhER6ZfGWItIrIt4YW1mZ5rZq2Y2I2hdppndYWY/MrN7zCxvMM+FxZw5kJICa9bAvn1h3bWISLyLpZytWUFEJNZFtLA2swlADnBMyFN3Ay84564BHgXuHeRzI5eeDoccAs7BypVh3bWISDyLtZzdefKixliLSIyKaGHtnNuJT7SdzGwScDHweGDVM8AiM5vR33NhDWzBAn+7YkVYdysiEs9iLWcnZyYD0N6gwlpEYlPEh4I450Iz4kKgyjnXEHi+FVgPnDLAc+GjwlpEpFexlLOT83xh3VrXGo7diYiEXUq0AwAmA9Uh6+qASUBrP8+Fz5FH+ttXXw3rbkVExqCo5ezkXF9Yt9W1hWN3IiJhFwuzgjigMWRdOtAywHPdVFZWUlZW1rlUVFQMPoJjj/WXN1++HOrqhhS8iMhwVFRUdOYroCja8QxB1HJ2Sp7vC1JhLSKRNticHQs91luBcSHrcoFt+CTd13PdFBcXs3Tp0uFFkJPje61few1eeQXOPHN4+xERGaTy8nLKy8sBMLOqKIczFFHL2R091q21GgoiIpE12JwdCz3WzwKTzCwTwMxSgWnAcwM8F14nnuhvX3gh7LsWERlDopaz1WMtIrEuGvNYW8ddAOfcDuAJ4PTA+jOAR5xzW/t7LuyBnRI4t2bx4rDvWkQkXsVSztYYaxGJdREdCmJmucClgYeXmdmdzrlK4ArgFjNbAJQC5UEv6++58Dn9dEhL80NBdu2CkpJROYyISLyItZytoSAiEusiWlg75+qAewJL8Po9wOV9vKbP58IqJwdOPRX+/nf429/gsstG/ZAiIrEs1nJ251CQWvVYi0hsioUx1rFj0SJ/+8gj0Y1DRER6SBnnC+uW3T0mGRERiQkqrIN9/ON+2r2nnvLDQUREJGaklqQC0FKpwlpEYpMK62ATJsDZZ0NrKzz4YLSjERGRIMnZySRlJNHe0E7bPg0HEZHYo8I61Kc/7W/vvx+ci24sIiLSycxILfa91s2VzVGORkSkJxXWoRYtguJiWLECnn8+2tGIiEiQzuEguzQcRERijwrrUBkZcMUV/v7tt0c3FhER6SatOA3QOGsRiU0qrHtzxRV+Tuu//AXefz/a0YiISEDnUJBdGgoiIrFHhXVvSkvh4ov9GOsbboh2NCIiEpA+OR2A/Zv3RzkSEZGeVFj35bvfhZQUeOABePvtaEcjIiJAxowMAJo2NkU5EhGRnlRY9+WAA+Dzn4f2dl9ki4hI1KVP9z3WKqxFJBapsO7P974HmZnw5z/D009HOxoRkYSXMd33WO/fqKEgIhJ7VFj3Z9IkuPZaf//KK6FJPSQiItGUMS0wFGRTE65d1xoQkdiiwnogX/86zJsHa9fCD38Y7WhERBJacnYyqcWpuGbH/m3qtRaR2BKWwtrMPhyO/cSk1FS45x5//8Yb4dVXoxuPiMgIxXvOzjokC4CGtxuiHImISHdDKqzN7Gwze8HM1prZusCyAXhkdMKLESed5Huu29rg0kuhvj7aEYmIDGis5uzsedkA7Fu1L8qRiIh0lzLE7f8XuBFYDrQH1hlwcTiDikk//CEsXgxvveXHW//612AW7ahERPozJnN21jzfY63CWkRizVAL63XOuZ+ErjSz1WGKJ3ZlZMBvfwvHHuvntj7ySPjKV6IdlYhIf8Zkzu7osW5YpaEgIhJbhjrG+gozO7+X9YvCEUzMO/RQ31MN8I1vwDPPRDUcEZEBjMmcnT0/MBRk5T5cm2YGEZHYMdTC+gTgF2bWFrS0A/eMQmyx6YIL4Jpr/Hjrj38cVqyIdkQiIn0Zkzk7rTiN9GnptNW30bBavdYiEjuGOhTkOqAcWAq0BdYZ8KmRBmJmc4AvAWuBg4AK59wbZpYJ3ALUAeOAbznnakd6vBG5/npYswYefhg+9CF4+WWYPj2qIYmI9GLM5uy8Y/Oo3FRJ7Su1nUNDRESibaiF9ZPAo865bpOHmtl9YYjlf4GPOOe2mtk04O/AHOBu4Cnn3ENmdjZwL3BRGI43fElJfpz1rl3w3HNw9tnwwgswfnxUwxIRCTFmc3becXlUPuwL64n/PjHcuxcRGZahDgV5ELjazE4OWk4Brg5DLHOA3MD9RiDfzCbhz15/PLD+GWCRmc0Iw/FGJiMD/vIXP+569WpYtAgaG6MdlYhIsDGbs/OOywNg78t7w71rEZFhG2qP9ffwP/ldBnScMWJAKTDSKTL+iB8L+CHgUuDLwEKgyjnXAOCcazWz9cApwIYRHm/kCgrgySfh+OPhpZfgoovgkUcgOTnakYmIwBjO2TlH5pCUkUTDqgaadzaTNiEtnLsXERmWofZY/xiY6pyb4ZybGVhmAB8LQyz/D2gCXgPqnXOPAJOB6pDt6oBJoS+urKykrKysc6moqAhDSIMweTI89ZQvsh99FL70JXA6S11E+ldRUdGZr4CiUTrMmM3ZyRnJ5J+cD8CeZ/YMI3wRkcEbbM4eao/1ZudcXS/rS4e4n95kAr/HJ+bbA70cDv8TY7B0oCX0xcXFxSxdujQMYQzD3Lnw+ONwxhnw859Dfr6//LkuICMifSgvL6e8vBwAM6sapcOM6ZxdeFYhe57ew+7Fu5lwyYQR7UtEpD+DzdlD7bH+f6ErzOwTwO1D3E9vfgM84py7DrgJn7C3488qD5YLbAvD8cLrxBPh97/3w0BuvtlfnbG9feDXiYiMnjGds8ed6Q+1Z/EenH4pFJEYMNTCutDMvgdgZllm9gvgF/ifAofNzIqABc65msCqG4BaYCMwKTB9E2aWCkwDnhvJ8UbNeefBn/8M6elwzz1w6aWwf//ArxMRGR1jOmdnH5pN2sQ0mrc1U/96/WgcQkRkSIZUWDvnzgEWm9kdwBvAbGAB0NuVvYZiN9BkZpOD1lUDbwJPAKcH1p2B7yHZOsLjjZ5/+zc/5jo31/dgn366n5ZPRCTCxnrONjOKPuaHOu76g/KsiERfv2OszezkPp7aCKwBbgWmAx8FvjbcIJxz7Wb2UeD7ZrYMmAB80zlXa2ZXALeY2QL8uMDy4R4nYhYu9PNbL1oEL74IRx/tx2Afdli0IxORMSwRc3bJJ0rYdtc2Kh+uZNZNszCd2yIiUWT9jUszs1fwc5XWAH0NGDZgonMuPfzhDV5ZWZmL2smLfdmxww8PeeUVyMmBBx+Ec86JdlQiEmPMbJlzriwM+0m4nO3aHC9NfomWnS0c+dqR5JXlhSE6EZG+9ZezBxoKch1wkHNuetBUTaHLDOCScAc9JpSWwrPPwsUXQ32978G+/XZNxycioyXhcrYlGyUXlACw49c7ohyNiCS6fgtr59zfnHMDDlxzzv0xfCGNMRkZ8JvfwA9+4GcJ+drXoLwcmpujHZmIjDGJmrM7Lmm+8zc7aWtoi3I0IpLIhjoriAyHGVx7rT+ZMSMDfvELf1Ljli3RjkxEJO7lLMgh99hc2va2UflwZbTDEZEEpsI6kj75SXj+eZg0CV54AQ4/HP7612hHJSIS9yZd7i/uuO3nsXeZAxFJHCqsI62sDJYvh7PPhupqOPdcuOoqDQ0RERmBkk+WkDIuhdpXaql5oWbgF4iIjIIRF9ZmdoSZ/cbMvmFmM0YeUgIoKfE91Tff7K/UeNtt/sqNq1dHOzIRGePGas5Ozk5m8pV+Wu3NN2+OcjQikqjC0WP9ReAM59xtwAc6rrglA0hKgv/8T1iyBKZNg9degyOO0KwhIjLaxmzOnvzlySRlJFH9RDX1K3UlRhGJvHAU1l8CjgrcfxBIDcM+E8cHPgBvvgmf+Qw0NflZQy691N8XEQm/MZuz00rSKP1cKQAbr9sY5WhEJBENubA2sw+Z2SWB+8fg50zdCv5qXM652jDHOPbl58Ovfw1/+pO/kMzvfgdnnAFVVdGOTETiXKLl7GlXTyMpI4nKhyupXTqm3pqIxIEhFdZm9kPgYeA8AOfcv4BPm9nFoxBb4jnvPD9byOTJ/lLoxx0H//pXtKMSkTiViDk7Y0oGk7/sx1qv+/a6KEcjIolmqD3WHwCmAM8HrXsSuClsESW6BQt8MX3kkfD++3D88fCtb2loiIgMR0Lm7GnfnkZKQQo1/6hh99O7ox2OiCSQoRbWi51zNUDw2XXnhzEega55rr/5Tf/4llv8iY2vvBLduEQk3iRkzk4tTGXa1dMAWPMfa2hvbo9yRCKSKIZaWK8zs2uAg83sk2b2CHAlcHP4Q0twmZm+oH7xRTjkED8V3wkn+GK7sTHa0YlIfEjYnD3lK1PIPDiTxncb2Xybpt8TkcgYUmHtnHsQ+CeQDHwaqAHOcs7dNQqxCfhx1suX++EgALfe6q/Y+Oqr0Y1LRGJeIufspPQkZt81G4CN12+kcYM6JERk9A15VhDn3CvOuS865z7snPt34FkzO30UYpMOGRlw003w8sswdy68957vvb7pJs15LSL9SuScPe70cZR8soT2xnbeu/w9XLvypYiMrpT+njSzx4CcAfZRBFQB/whXUNKHY46B11+Ha66Bn/wErr7a92bffz9kZ0c7OhGJMuXsng68/UD2PLOHPc/sYds92zqvzigiMhr6LayBNwPbvAu0A5cDy4GlQdvMBApHJTrpKT3dXwL9lFP8hWT+8Ad491347W9h3rxoRyci0aWcHSJtQhqzfz6bVeev4v1vvs+4s8aRdVBWtMMSkTFqoML6x4AFzirHzA5wzn0/dCMze3Q0gpN+LFrkZwn5yEdgxQo/Pd9118E3vgEpA/2zisgYpZzdi+KPF1NySQm7fruL1Z9ezeFLDicpNRwXHhYR6a7fzOKc29uRoAMmmFly8DZmthA4JlwBmdlxZnaVmX3UzIrMLNPM7jCzH5nZPWaWF65jxb25c2HZMigvh+Zm+Pa34cQT/QwiIpJwlLP7dtAdB5E+JZ3aV2pZd7UuHCMio2Oof7L/EVhhZj83sx+b2RPAYuC/wxGMmX0WONc5d6tz7i/OuSrgbuAF59w1wKPAveE41piRlwf33gtPPQVTpvjZQg4/3A8XaWuLdnQiEl3K2QGp41KZ++BcLMXYctsWKv9SGe2QRGQMGup0e4uBs/AnvhyMn7rpfOfcjSMNxMxOAS4Bvhe0bhJwMfB4YNUzwCIzmzHS4405H/wgvPUWfPazsH8/XHUVnHyyn0FERBKScnZ3+SfkM+vmWQCsvmw1je9rCj4RCa/hTLe3zTn3XefcIufcpcDjZnZUGGL5CfAOcIeZPWVmJwALgSrnXEPg2K3AeuCUMBxv7Cko8DOEPPEETJwIL70Ehx0GV14J6/TTp0giUs7ubsrXplB0XhFte9t46yNv0bq3NdohicgYMqSz3Mzs/l5WFwNtwEeHG4SZHQwcAXzWOfemmX0TeAq4AagO2bwOmBS6j8rKSsrKyjofl5eXU15ePtyQ4ts558CqVfDVr8IDD8Ddd8M998CZZ/rx2IsWQWpqtKMUSWgVFRVUVFR0PCwajWMoZ/dkZhzyq0N4/d3XaVjVwKoLVnHoXw8lKUUnM4pI3wabs80N4QIjZvY8fu7T4BedADzvnLt+GHF27PdjQIVzrijwOBOoBK4DPu6cOzZo29eB3znnbg3eR1lZmVu6NHhGKQF8gX3LLfDQQ36ICMCECXDZZfD5z8OBB0Y1PBEBM1vmnCsbeMsh71c5uw+N6xt5/djXaalsYdIXJ3HQ3QdhZlGLR0TiR385e6jzsn3eOfduyM4PAS4bZmzBcXTG4pxrNLM1+MvwjgvZNhfYNsLjJY558+B//gd++lP4zW+gosIX2zff7JdTT4XPfAY+9jHIzY12tCISXsrZfcicmcn8R+fzxqlvsO3n20gtSmXm9TOjHZaIxLmhnrz4bi+r9wOfH2EcbwL5ZjYhaF0rPhlPCvSGYGapwDTguREeL/EUFsJ//Ic/wfHFF32PdWYmPPusv19aCp/6FCxerNlERMYI5ez+5R+fz9zfz4Vk2HjDRjb+cGO0QxKRODekwtrM1pvZuqBlI/7klT+PJAjn3GrgSeD8wHEKgDTgt8ATwOmBTc8AHnHObR3J8RKaGXzgA/CrX8H27XDffXDSSdDQ4Hu0zzoLpk+Hr3/dT903hKFCIhJblLMHVnxeMXMemAMG67+7nk0/3hTtkEQkjg11jHUF8Du6xuu1AzudcyOe083MioCf4S+9Ow24zzn3tpmNA24BNgClwNXOufrQ10d7vF7ce/99X1g/8ED3GUSmT4cLLoALL/RXd9QYRJGwG8Ux1srZg7T9V9t593O+g3/ad6Yx8/qZGnMtIr3qL2cPtbA+2Tm3JGTdwUCbc27tyMIcmVhL0nHLOX+p9Icegocfhm1BQyMPOMAX2eecA8cco5lFRMJkFAtr5ewh2PE/O1j976uhDSZ+YSKz75qNJau4FpHu+svZQ51f6Oxe1u0EHh5yVBKbzOD44+H222HzZliyxM+DPWGC79W+8UZ/2fTCQjj3XL/dW29Be3u0IxeRnpSzh6D0M6XM/9N8LN3Yfu92Vn58Ja11mudaRAZvUIW1mX3JzGqAb5lZW/CCn7N056hGKdGRlOTHX995J2zdCv/8J3zpS3DIIVBfD3/9K3zta/4iNCUlcN55fvaR11/XCZAiUaScPXxFi4pY8PQCUgpSqH60mtePf53GdbpCo4gMzqCHgpjZYcBngP8OearZObc93IENVSz+rDimbd0K//gHPPOML7i3hpyblJcHJ5zgL6t+8slQVgZpadGJVSTGjcZQEOXskWl4r4GVH1lJw+oGUgpTmPu7uRR+sDDaYYlIDAjnGOtU51xLL+unOeeieip1rCfpMc05WL/eDxtZsgSee67nJdQzM/0QkxNPhOOOg2OP9cNJRGQ0x1grZ49A695W3r7kbXb/dTcAU6+ayswfziQpTVdpFElkwy6szewq4G/OubcDj7+DvwBAt82AU51zC8MT7vDEQ5JOKFu2wPPPdxXbb7/dc5vZs32RfdxxvuiePx9ShnrNIpH4F67CWjk7/FybY9NNm1h/7Xpog5yjcpj7u7lkzc6KdmgiEiUjufLiTKAg6PGRQDawg67pmww/pZJIlylT4KKL/AJQWekL7Zdf9rOOLF0K773nlwce8NtkZcHRR/tC++ij/fCRadM0xZ/I4Clnh5klG9O/M52C0wp45+J3qF9Wz9IFS5nxXzOY8vUpJKWo91pEugx1KMiRwErnXHPI+rkdPSTREi+9HxLQ0gIrVvgiu2N5//2e240fD0cd5Yvso47yi4ptGWNGcSiIcnYYte5tZc2X17Dzf/25nzlH5HDwLw8m94jcKEcmIpE0kh7rUIc6514P2fnBwAFAVJO0xJnUVF8sl5X5mUbA92q/+mpXj/bSpVBdDU8/7ZcORUVdRbaKbZH+KGeHUUp+CnMemEPJxSW894X3qF9ez7KjlzH5ysnMuHYGqYWa218k0Q2qx9rMTsL/fFgOVIQ8XQj8t3NuWvjDG7x47P2QATgHmzbBsmV+WbrU31ZX99w2tNguK4OpU1VsS1wId4+1cvboa61vZf131rP1zq3QDimFKcy8biYTvzBRw0NExrgRzwpiZocCD+DH7+0JeboZ+K1z7rqRBjoS8Z6kZZA6iu2OIrtj6avYPvxwP8/2ggX+ds4cSE+PfNwi/RiFwlo5O0LqV9Sz9qtrqfm/GgCy5mUx68ZZjD93vC6JLjJGhWW6PTPLBz7hnPtFOIMLl7GSpGUYnIONG7sX2kuXwu7dPbdNSfEXuOkotDtuS0vVuy1RM0rzWCtnR4hzjqq/VPH+Ve/TtK4JgNxjcpl5/UzGnTlOBbbIGBO2eaz7OcCJzrkXRryjERhLSVrCoKNne8UKePPNrts1a/xzoYqKuorsww6DefN873ZOTuRjl4QzWicv9nM85exR0L6/na33bGXTjZto2eWnD88/KZ8Z/zWDgoUFKrBFxoiRzGP9GDBQZWHAPOdcyfBDHLmxmKRlFOzbB6tWdS+233wTamp63376dJg71xfaHbdz5kCuZgGQ8AnjPNbK2TGgbV8bW+/ayqabN9G6uxWA3KNzmfqfUyk+rxhLVoEtEs9GMivIu0ANsJWuOVB77B84d/jhiURQdjYcc4xfOjgHmzd3FdlvveWL73ff9UNMNm6EJ5/svp+pU7sX23Pn+iUvL7LvR6Q75ewYkJydzLT/nMakL05iy8+2sPW/t1L3Wh1vf+JtMg7IYOpVUyn9TCnJmaHX7hGReDdQj3URUOec29/vTszGO+d6OXsscsZy74dESWurn1v77bd9od1xu3o1NDf3/prJk32Pdkeh3XG/qCiysUtcCWOPtXJ2DGpraGPHr3ew+bbNnWOwUwpTmPi5iUy6YhKZszKjHKGIDEUkxlh/PtonyCRSkpYoa22Fdet8od1RbHf0cDc19f6aoqKexfacOTBpkk6alGiMsVbOjoL21naq/lTF5h9vpm5pnV9pUHh2IZOvnEzh2YUaJiISB0atsDazacCVwJecc9nD3lEYJGKSlhjT1gYbNsA773QV3R336+t7f01eXu893NOnQ5Lmwk0UkSqslbNjR+2/atl691Z2PbgLt99/D6dPT6f0slJKP1NK5kz1YovEqrAW1maWDCwCLgfOBHYG9jN5pIGORKInaYlhzsHWrT2L7bff7n1KQIDMTD8tYGgP9wEH+KtWypgymoW1cnZsa65qZsevdrDtnm00re/6xatgYQGlny2l+OPFJGdrLLZILAnXPNYzgc8DnwXGAxuB/wD+DhzvnHspDIGmAa8BX3HO/Z+ZZQK3AHXAOOBbzrna3l6rJC1xxzl/GffgYrvjdvv23l+TmgoHHdSzh3v2bMjIiGz8EjajNI/1qOfswHGGlbeVs7tz7Y6aZ2vY/qvtVD1SRXtTOwDJOckUX1DMhIsn+Cn7NFREJOqMr5oZAAAgAElEQVRGMisIZvYJfE/HacBa4CfA/wDXOOeeCmwWlgQNfBOYEfT4buAp59xDZnY2cC9wUZiOJRJdZlBS4peFC7s/V1PTs9h++20/Q0nH/WBJSTBrVs8ebs3FnXAinLNBeTssLMkYd/o4xp0+jta7Wtn1h13s+PUOal+qZcf9O9hx/w5SJ6RSfH4xJReWkH9CPpakIlsk1gxYWAON+EvgvgR83jn3HoCZjfysxyBm9gFgO4HL75rZJOBi/HhAgGeAR8xshnNuQziPLRJzCgrg+OP9EmzfPj8rSWgv99q1Xctjj3V/zdSpflhJ6DJxok6cHJsikrMD+1TeHgUp+SlMunwSky6fRMO7Dez43x1UPlRJ49pGtt21jW13bSNtcholF5RQcmEJuUfnqsgWiREDFtbOuSeAJwInvVxuZjOAR4HOgZ5mZm4EZ0GaWTZwvnPu62b2/cDqhUCVc64hEEerma0HTgE2DPdYInEtOxuOOsovwfbvh/fe69nL/d57fo7uzZth8eLur8nN7b3gPuAASE+P3HuSsIpEzg7sQ3k7ArIOzmLWDbOYef1M6pfXs+vBXex6aBf7N+1ny0+3sOWnW0iblEbRR4oo+kgRBacWkJSmE59FomU4Jy+mAB8FLgFWAr8EznLOVQw7CLPrgDudc7vMbANwGXA08Cnn3GFB270MPOacuzF0H9OnT3fFxcWdj8vLyykvLx9uSCJjQ2urn6lk9eruyzvv9H3iZHKyH1bSUWjPnQvz5/thJdlRnUhiTKmoqKCiwqfNZcuWbXTOzRiN44xGzg7sd0R5Wzl7+Jxz1L5ay64Hd1H5x0qat3bNq5+cm0zhhwop+mgRhR8qJLVAJzuLhMNgc/ZIp9s7EPgycLlzLmuY+/ggkO6ceyzweAM+QZcBn3DOHRu07evA75xzt4buRyfCiAxRVVXPgnv1ali/Htrbe25vBjNm+CJ73jy/zJ/vi2+dODkiEZxub8Q5O7CfEedt5ezwcM5R/3o9VX+pourRKva9ta/zOUsxChYWUPjhQgrPLiTrkCxMw79ERiwSF4j5tHPugWG+djFwRNCqcUA9/oSbS5xzs4O2XQNc65z7Xeh+lKRFwqSpyY/V7hjLHXwBnNbWntsnJfnhI6EF9+zZkJYW+fjjUBQuEDPsnB14/YjztnL26Ghc10jVY1VUP1pNzZIaCPobOX16OoVnFzL+Q+MpOK2AlNzBnGYlIqFGvbAeCTMrBoIHdL4MfB14FXgbKHbONZpZKj5xz3LObQ3dj5K0yChrboY1a7oK7ZUr/e2aNb33cKek+KkBQwvuAw/0z0mnSBfWIxWOvK2cPfpaqluofrKa3U/tZs/f99BS1dL5nKUa+SfmU/ihQgrPKiT70GydACkySCOabm+0Oecqgx+bWRtQ6ZzbZGZPAKcDTwBnAI/0VlSLSASkpXUVyMGamnxvdmjBvW6dH8v9zjvw8MPd93PwwT0L7pkz/fhuiXnK2/EhdXwqpZeWUnppKa7dUbesjt1P7Wb3k7upfbWWmmdrqHm2hnX/uY7UolQKTi2g4LQCxp0+jswDMzVsRGQYol5YD+AK4BYzWwCUAjqzRSTWZGTAggV+CdbQ4Ivq0IJ740Z46y2/hO5nzpyeBfe0abq8e3xR3o5BlmTkHZ1H3tF5zPjeDFp2t7DnmT3sfnI3e57Zw/4t+6l8uJLKh/3fTOlT0juL7ILTCsiYovMoRAYj6kNBwkU/K4rEibo6X3B3FNodRffWPjo1s7O7ZiYJLrgnTx4z83DH21CQcFDOjh3OORrXNlLzzxr2/GMPNc/WdBs2ApB5UCYFCwvIPymf/BPzyZiRoR5tSVgxPcY6XJSkReJcTY0/WTK04N65s/ft8/N7L7gnTIi7gluFtcQS1+7Yt3KfL7L/WUPNczW01bV12yZtchoFJwUK7ZPyyZ6nMdqSOFRYi0j8qq7uPpSk4351de/bFxZ2FdnBBXdRUWTjHgIV1hLL2lvbqV9WT82SGvY+v5e9L+yldU/3GYJSClLIOyHPF9sn5pNblktSuoZwydikwlpExhbnYNeunuO3V63yPd+9KSnpXnDPnw+HHeavQBllKqwlnrh2R8M7DdQ831Vo79+0v9s2lmbkHJFD3nF5fjk2T8NHZMxQYS0iicE52Lat94K7vr7n9mZ+hpKysq7l8MMjfoVJFdYS75o2NXUW2TXP19CwqqHHNqklqd0K7dyjczWXtsQlFdYikticg82buxfaHTOTtHQ/SYukJD92+6ij4Ljj4KST/GwlozgziQprGWta97ZS+1otta90La3VIReYMsien+2L7GNzyTs6j6y5WSSlagiJxDYV1iIivdm/3xfXS5d2LStXQlv3E7UoLIQTT/RF9kknwZFHQmpq2MJQYS1jnXOOpnVNXYX2q7XUL6/HtXavQSzdyDk8h9yjcsktyyX3qFxfbKeo2JbYocJaRGSwGhthxQp47TV48UV4/nk/vCRYVpbvzT7rLPjEJ2DWrBEdUoW1JKK2xjbql9dT+2qg0F5WT+Paxh7bJWUkkXN4DjlH5XQV23NUbEv0qLAWERku52D9enjhBV9kP/+8v9JksJNPhttvhyOOGNYhVFiLeC01LdS/Xk/d0jrqltVRt7SOpnVNPbZLykzq7NnOOTyHnCNyyJ6XrZlIJCJUWIuIhNOuXbBkCfz5z/DYY/7EyIwMePppP1RkiFRYi/StZXcLda/XUb+sq+BuWt+z2LYUI2tOli+0g5bUwvAN2xIBFdYiIqOntha+/GV44AE/JGT16iGPv1ZhLTI0LdUt1C2ro355PfVv+KXh3QbopaRJn5pOzhHdi21N/Scj0V/O1jw3IiIjkZcHv/yl77let85fPXLBgmhHJTKmpY5PpfCsQgrPKuxc17avjX0r91G3vK6z2N735j72b97P/s37qX6s66JSyfnJ5CzI6RxGknN4Dtlzs0lK01ASGRkV1iIiI9HUBN//vr8wTUYGTJ8e7YhEElJydjJ5x/o5sju4NkfDmobOQrt+uV9aKlvYu2Qve5fs7dy2YyhJ9mHZ5Bya428PyyFtUpp6t2XQVFiLiAxHTY0f/vHTn8KGDZCcDPfdBwUF0Y5MRAIs2cg+JJvsQ7KZ8MkJgJ/6r3lHc1exHVga1zSy76197HtrH7vY1bmPlHEpnUV2Z9E9P5vk7ORovS2JYSqsRUQGa+tWeOopePRRWLzY91aDvzz63XcP68RFEYksMyN9YjrpE9MZ/6Hxnetb61tpWNVA/Vt+CEn9m/62dU8re5/by97n9gbtBDIPyCT70OxuRXfmrEwsSb3biUyFtYhIb1pa/Hjp117rmmZv/fru25x+OlxxBXzkI5CidCoSz1JyUnoOJXGO5m3NnUV2R9Hd8E4DjWsbaVzbSNWfqzq3T8pKInt+V6Gdfajv4U4dr5lJEoW+CUQksTkHlZXwzjv+KozLl/tl1Spobu6+bV6e75VetAjOPRcmTYpOzCISEWZG+uR00id3791ub26n4d2Grp7tt/xt89Zm6v5VR92/6rrtJ21yWrdx29mHZZN1cJZOlhyDVFiLSGJobvZjodeu9UV0x7J6Neze3ftrDjzQX778hBN8QX3YYX4stYgktKS0JHIOzSHn0BwmXDKhc31LdYvv1X5rX1fRvXIfzVub2b11N7uf6so1lhqYd/uw7gV3WqlOloxnKqxFZGxwDqqq/JR3vS2bN/ttepObC3PmwNy5/uqJRxzhp8zLy+t9exGRXqSOT2XcwnGMWziuc51rdzSua+zeu72inqZ1Tex70xfg3fZRlNr9ZMnDcsiam0Vypv6ojwcxU1ib2ceA24A84DfAN5xzrWZWBPwQqAFSgW8551qiF6mIREVDgy+O+1o2bfJXQOxLUpKfCm/WLF9EH3KIv50zByZOBPUQDYlytsjgWJKRdWAWWQdmUfyx4s71nSdLvtn9ZMmWqhZq/llDzT9runaSBFmzs3oU3OnT0tW7HWNiorA2s2nAR4HzgXnAPcBm4FbgYeDbzrlXzeyLwHXA1dGKVUTCrL3dj3Hevt0v27Z13QYXzn0N1wiWlwcHHOCL546l4/G0aUO+IqL0TjlbZOT6Olly/5b93Qrt+jf9VSUbVvul8g+Vndsn5yX3GEqSPT+blNyYKO8SUqy0/HTgc865VmCZmc0HTjWzF4CDnXOvBrZ7HHjXzG5wzu3ra2ciEgPa2mDXrq5CObhoDr6/cye0tg68v7Q0mDIFpk7texk3Tj3PkaGcLTIKzIyMqRlkTM1g/DldJ0u2NbXR8E5D94J7ReBCNy/sZe8Le7vtJ2NWRo+CO3NWJpas/DjaYqKwds49H7JqK5ALnAasD9pua+AnjzLguYgFKCJ+fHJ9ve9d3rWr9yX0ufb2we27sNDPsDFxYtftxIndi+biYj+cQ6JOOVskspIzksk9IpfcI3K7rW/e2dxjKMm+t/fRtK6JpnVNVP2l76kAcw7LIXtBNqkF+iUvnGKisO7F0cCPgXKgOuS5OkBzXImMlHNQWwvV1f6kv46ltwK5Y+m4IMpgFRf3LJZD75eWQnr66LxHiRTlbJEoSJuQRuGZhRSeWdi5rr2lncb3GnsU3Pu37O91KsCMAzLIPTKXnCNz/O0ROaQVp0X6rYwZMVdYm9mBQKVzboWZOaAxZJN0oMeJMJWVlZSVlXU+Li8vp7y8fFRjFYkZHb3JwQVyaMHc27rBDMEIlpkJEyb4grmkpPsSvK7jfpqSc18qKiqoqKjoeFgUzVhGQjlbJLYkpSaRPS+b7HnZcFHX+pbdLZ3zbXcMJal/s56m95toer+Jyoe7xm6nT03vKrQDt2kTE3sawMHmbHN9TT8VBWaWDNwEfNc5t9/MrgZOdc6dFXjegKbAupeCX1tWVuaWLl0a8ZhFws452LdvaAVydXXPi5kMRk4OFBV1LePH910sl5RAdnb4369gZsucc2UDbxlblLNF4lt7SzsNqxuof72eutfr/O3yOtr39RzGlzohldwjugrt3GNzyZiSEYWoo6+/nB1rPdZXAbc65/YHHi8Gvhj0/AygGVA2lvjgHNTVdRXCg73dv3/gfYfKzu5eIIcWzMGPO9ZpCIaMjHK2SBxLSu260E3pZ0oBcG2OxrWNXYV24LZlZwu7n+p+kZu0SWmdM5vkHptLblkuKTmxVlpGVsy8ezO7Bp98M81sFrAQWAJsN7P5zrmVwNnAHc65YXTNiYyQc7B379AK5OpqaBnGFL5ZWX0Xw72tGz/eD9MQiRDlbJGxyZKNrIOzyDo4iwkX+atKOudo2tDUWWjXLfVjtZu3NVP15yqq/hw4STIJsudldxbaecfmkT03O6FmI4mJwtrMvgtcH7J6tXPufjO7EPi+mW0AcoDvRDo+GYPa22HPnr6L4b7WtbUN/VgdPckdRXFvt8EF8/jxvrAWiVHK2SKJxczInJlJ5sxMij/uL3Lj2h2NaxqpfbW2c9m3Yp+/nPtb+9j+i+0AJOcmk/eBPApOLiD/5Hzyjs4jKX3szvAUU2OsR0Lj9RJYa6u/eMhQepF37+778tb9ycvrv0DurWDOSMwxaDJ48TrGeiSUs0XGnrbGNuqX13cW2nWv1tG0oftsUpZu5B0bVGgfnxd3w0fiaYy1JLrm5q7e4cEWyjU1A++3NwUFgy+Qi4r8XMua5UJERKRXyZnJ5H8gn/wP5Heu279tP3uf30vNkhr2LtnLvpX72LtkL3uXBC5qkwx5x+ZReFYh484aR+7RuSSlxG+PtgprGT1NTUM/aa+ubuD9hjLzRe9QepILCyFFH38REZHRlD4pnZILSyi5sASAluoW9r7YVWjXvV5H7Uu11L5Uy4YfbCClIIWC0wsoPKuQwg8WkjE9vn71VWUhg9fe7gvgnTv7XjouLlJd7aeMG6rkZF/0DqUnuaDAv05ERERiWur4VIoWFVG0yE8F3VrbSs3/1bD76d3s+fseGtc2UvVIFVWP+BMicw7Poei8Ioo+WkT2odkxP5e2CutE19rqC+GBiuWdO/12Qzl5LyVlaAXy+PGQn6/LVouIiCSIlLyUboV247pG9izew+6/72bP4j3Uv1FP/Rv1bLh2AxmzMig6r4gJF00g58icmCyyVViPVc75E/S2boVt2/zScT943c6dvid6sMaN81feC15KSro/Li72hXJurh+mISIiIjIImbMyyfxCJpO+MIm2pjZq/lFD1V+qqHq0iqZ1TWy5bQtbbttC1rwsSj9dyoRLJpA+OXauyaDCOh4553uRN270y6ZNsGVLz6J5MBcZMesqjEML5NDiWZeoFhERkQhJzkhm/DnjGX/OeGb/fDZ7X9pL5R8q2fXgLhpWNbDuW+tY9+11jD9nPFO+OoWC0wqi3outwjoWtbX5ArmjcO5taWoaeD/5+TBpEkye3PfthAmQmjr670lERERkmCzZKDipgIKTCjjgtgPY/dRudjywg+rHq6l+wi/Z87OZ8rUpTPjUBJJSozOsVIV1tLS3+17mNWt6LuvW+Wnn+jNuHEyf3rVMndq9aJ40yV+YRERERGQMSUpL6hyX3byrmW33bmPb3dvYt3If7/77u2z80UZm/tdMSj5ZEvGrPqqwHm1tbbB2Lbz1ll9WroT33vPr+ut1Li3tXjh3LDNmwLRp/kIlIiIiIgksrSSNGd+bwbRvTWPXH3ax8fqNNL7XyDuXvsPm2zZz8H0Hk3tUbsTiUWEdTvv3w/Ll8Oqr8MYbvpBetarvAnrCBDjooJ7LgQeqt1lERERkkJLSkii9tJSST5aw8393suHaDdQvr2fZMcuY+o2pzPzhzIgMD1FhPRK7dsE//wmvvOKX5ct7H8IxdSocemjXMmeOL57V6ywiIiISNkkpSUz87ERKLihh/bXr2fLTLWz+8Wb2vrSXeQ/PI33i6M4gosJ6KJqbYckSePppv6xY0XObuXPhuOPgqKN8ET1/vh8PLSIiIiIRkZydzIG3Hkjxx4tZ9YlV1L5Yy/ITlnP4s4eP6tUcVVgPpK3N90o/9BD86U+wZ0/XcxkZcPLJcNJJvpg++mg/E4eIiIiIRF3+8fmULSvjrX97i7rX6nhj4Rsc+cqRpE0YnemDVVj3paYGfvlLuPNO2LCha/28eXDOOXDmmXDiib64FhEREZGYlDYhjQWLF7DizBXUvVbHqgtXseCZBSSlhH/MtQrrUI2N8LOfwY03Qm2tXzdrFlxyCVx4oS+sRURERCRupOSnMP/R+Sw7chl7n9vLtru2MeUrU8J+nOjMnh2rXnrJF85XX+2L6lNPhccf93NLX3edimoRERGROJU+MZ3ZFbMBWP/99bTsaQn7MVRYd6io8OOl16/3Jx0uXuzHVp97LiSpmURERETi3fhzx1NwWgFttW1s/+X2sO9fFSPAfffBF77gT1S86ipYuhTOOCPaUYmIiIhIGJkZU77qh4DsuH9H2PevwnrZMrjySn//jjvgxz+GtNE5U1REREREoqvw7EKS85NpeKeBxnWNYd13XBTWZpZpZneY2Y/M7B4zC8+VVdraoLwcWlp8cf2lL4VltyIiiWzUcraISBgkpSZReGYhAHsW7xlg6yHuO6x7Gz13Ay84564BHgXuDctea2th8mSYMgVuvjksuxQRkVHK2SIiYZJ/Uj6pxam0N7WHdb/mnAvrDsPNzCYB64FxzrkGM0sB9gLznHMbOrYrKytzS5cuHd5Bdu2CkpIwRCsiMnRmtsw5VxbtOMIhIjlbRGSE2lvasRTDzIb82v5ydjz0WC8EqpxzDQDOuVZ80j4lbEdQUS0iEi4LGe2cLSIyQkmpScMqqgcSDxeImQxUh6yrAyYFr6isrKSsrOuPh/LycsrLy0c/OhGRYaioqKCioqLjYVE0Ywkz5WwRGXMGm7PjobB2QOgpm+lAt1m9i4uL0c+KIhIvggtJM6uKcjjhpJwtImPOYHN2PAwF2QqMC1mXC2yLQiwiItI/5WwRSVjxUFg/C0wys0wAM0sFpgHPjXTHQV36gtojlNqji9qiO7VHv0YtZ0N8tn08xgzxGbdijpx4jDsSMcd8Ye2c2wE8AZweWHUG8IhzbutI9x2PH4rRpPboTu3RRW3Rndqjb6OZsyE+2z4eY4b4jFsxR048xh2JmONhjDXAFcAtZrYAKAV0houISOxSzhaRhBTz81gPlplVAhuH+LIiYCydNDRSao/u1B5d1Bbdhbs9pjvnisO4v5g3zJwN8flZjMeYIT7jVsyRE49xhyvmPnP2mCmsRURERESiKebHWIuIiIiIxAMV1iIiIiIiYaDCWkREREQkDBKysDazTDO7w8x+ZGb3mFletGOKJDP7mJmtN7NqM/uZmaUE1heZ2b1mdrOZ/SQw/2zCMLM0M1thZgsDjxP9c3KcmV1lZh8NfDYSsj3MbI6Z3WVmXzOzu83s8MD6hGyPSBhs25rZh8zszsC/y/mRjjMklsHGXGBmtWbmAst5kY41JJ4zzexVM5vRzzYx084dBhl3zLR1X9+7vWx3mZn91Mx+YWYnRzrOkFgGG/NBZtYW1M5HRDrWkHhOMbN3zWyPmd3ez3aj87l2ziXcAvwKuDBw/2zg99GOKYLvfRrwAHAU8GlgH3BV4LlngWMD978I3BjteCPcNt8B9gIL9Tnhs8ANIesSsj2ApcDkwP1pwDuJ3B4RavMB2xaYDbwOJAUeLwHmx3LMgee+CZyHn9/7DCA5ijFPCMTigBl9bBNT7TzYuGOprfv73g3Z7lTgb4H7qcCbwIRYjjmw7c3AWYE2PiXKn4184Br81V/PAVqAM3rZbtQ+11F781Fs9EnAfiAr8Dgl8IGZEe3YIvT+TwJSgh7fAvwVOA7YFrR+MlAPZEc75gi1yweAzwEbgIWJ/DkBTgGeITBrUGBdIrfHPuCQwP1i/KW5E7Y9ItDeg2pb4OfAj4Iefwf4VYzHnA4sBg6KdjsHxZRE/4V1zLTzEOOOmbbu63u3l+2eAsqDHt8HXBvjMU/AXxBqSrTbORBPZsh31zLg1F62G7XPdSIOBVkIVDnnGgCcc63AenwxMeY5554PvOcOW4FNwGn4dujYruMqaWURDC8qzCwbON85d3/Q6oUk7ufkJ8A7wB1m9pSZnUBit8cfgV+YWS5wKfBlErs9RttCBte23XIW8G7gtdGwkMHFfC5wLPCemf3VzIoiGmUvnHPtA2wSS+3caRBxx0xb9/O928nMkvGfl5ho68HEHHARvqd9s5n9KvB9GjXOuUYXqJQDsbwD/F8vm47a5zoRC+vJQHXIujp8j0MiOhr/l1sit8u3gJtC1iVke5jZwcARwH3OuS8B/8D3okwlAdsj4P8BTcBrQL1z7hES9PMRIYNt29Dtotn+g4o58NnJxxdQM4DHzCzWv4djqZ0HLcbbuuN7N1ghkEHstnVvMeOcux3IAxbhh4L8MsJx9crMTgKexA+pyeplk1H7XMfKhyySHNAYsi4dPw4noZjZgUClc24FCdouZvZBYKlzblfIUwnZHsA8YLdz7s3A4zsBI3HbA/xPi78HfgfcbmZnkNjtMdoG27ah20Wz/Qf9eXDeEnyP2SH4XtVYFkvtPCSx2NYh37vBOq7WF3Nt3U/MADjn2pxzjwMfBD5mZhMjGmDv1gK/Bk4HftzL86P2ue71DM8xbit+UHuwXPy4yYQR+NnpC8C3A6u24gfzdzxvQDZjv12uAo7wbxfwn41H8cMhEvFzkkJQXnDONZrZGiCZxGwPgN8An3TO1ZiZwxfZXyVx22O0DTZHh24XzfYf8veKc26nmT2EP0ns5VGMbaRiqZ2HJVbaupfv3WDV+HH6MdXWA8TcjXPubTP7B76dt492bAPEsh2438za8Sexhhq1z3Ui9lg/C0wys0yAwJRy04DnohpV5F0F3Oqc2x94vBg4OOj5GUAzfkaEsexi4PCgZRvwefwZ/on4OXkTyDezCUHrWgmcsJdo7REYl7nAOVcTWHUDUAtsJAHbI0IGm6NDc9bswGujYbjfK23AG6Mc20jFUjuPRCy0dej3bqfAuOBniL227jPmPuwDVo9iPEO1FF9Ehxq1z3XCFdbOuR34M1hPD6w6A3gk6GS9Mc/MrgGWA5lmNsvMPgfUANvNbH5gs7OBO5xzzdGKMxKcc5XOuS0dCz75VjrnNpGAnxPn3Gr8uLTzwc8DC6QBvyUB2wPYDTSZ2eSgddX4P0ASsT1GXX852vy86h1fhj/H//Tc8QvbacB/RzpeGHzMZna+mc0K3D8EP+zq3WjE3MG6fq6zoHUx2c7BBoo71tq6t+/dwPzP15lZaWCzO/FTxGFmafip7n4RnYgHF3Ng3u2SwP2TgOecc3ujGHNWyDzaHwZ+FnguIp9rC5w8mVDMbBx+6pgNQClwtXOuPqpBRYiZfRe4PmT1aufcHDObDnwf3y45wHdCzgoe88xsA3CZc+7/EvVzEuil/Rn+L/1p+BMZ307g9liAP4FxGX5qqSXOuecStT0ioa+2NbNl+Pn1/xjY7lLgyMDLnnPOPRqNeAOxDBizmd2Jn1nmn/iZCu4cxOwWoyZoppu7gesC8VTGcjsH4hkw7lhq676+d/Ft+jZwkXPulcC23wDG44cmPOicezGSsXYYbMxm9ihwAr4H+Enn3AORjbQ7MzsGP4XwauAV4F/OuYcDz0Xkc52QhbWIiIiISLgl3FAQEREREZHRoMJaRERERCQMVFiLiIiIiISBCmsRERERkTBQYS0iIiIiEgYqrEVEREREwkCFtYiIiIhIGKiwFpGYYmbFZvaYmW00s29HOx4REemf8nYXXSBGRGKKmX0QeBbIw189a7Zzbnd0oxIRkb4ob3dRYS0iMcvMlgCnOedaox2LiIgMLNHztoaCiEhMMrOpwNOJmpxFROKN8rYKaxkFZvafZnZXtOPoEI54zOwMM3vZzBaGKawxxcxuDfP+DLgAuCmc+xWR7pSvR87MZpjZfWZ2b5j2l2xml5vZxkFse+e7NaIAAAmZSURBVI6ZvTrc9xrO3K287amwjjNmNsfM/mRmvzCzFWbmzOzSaMcV4m/A76MdRJBwxLMZWBCGWGKSmX3YzKrMbLuZHW9mUwOfM2dmd5pZQWC7Y8xspZndZmaZgXUTgB1hDuli4JfOudbA/kXijvL1sMRjvt4VuE0P0/6SgBpg2iC2XQrM6etJM5tnZg8FPnvfNbPSoOfCnbuVt1FhHY/+AvyXc+7zwBHAAx1PmNlhZnZKNIIKPrZzbqVz7oVoxNGbocTTVxs6594FqsIZVzj/vUa6L+fc34C7gVbn3MvOuc3A5UAz8Kpzriaw3b/wifybzrnGwMv/DXhiRG8giJl9D/gB8IyZrQYODte+RSJM+XqI4jFfO+cagG3hOpZzrgVYNshtd+KL8L6eXwVUBB7e6pwLLqTDlruVt7uosI4jZlYIzAb2ATjn2oHvA87M8vFJ26IQV9SOHU6DeB/tETxWNPb1EDDFzI4CcM5VA/8Azg06VgGwLfDZ63CIc271CI/dyTl3vXPuIOdcmXPuEOfcknDtWyRSlK9HVwzm63DPBDGU+Afati1wGzruOWy5W3m7iwrrOBKYuuZd4CkzOy6wbiPwL/y4pplAuZl9FsDMrjOz7wV+BuocR2VmZWb2EzP7qpk1mNm7ZvZ1M/uomb1mZueZ2Z/NrNbMvhD0ul73F3Lsy83s18FjzQI/h95hZjea2bNmdkxgfb/HC3q9mdkXzGyVmZ1pZm+a2U4zu2wQxzgoOJ4BjtmjDXsxLvDTbp2ZXR7Spj8ItM2fAol4sG3W41hmdr6Z/cvMLjOzNWZ2Uz/76+3fvtd4+hPo2VgFfCJodTFwjgWGfQCLgMeC4swEGoIeH2JmN5nZrWb2nPkhJRea2VIz+4iZ/d7MaszsK2Z2VOA9VpvZhweKTySeKF8nTr7uTV/vsb/jmNlMM7vffO/vjwdznF6O+zcz+62ZnTDAdsrdo8U5pyWOFmAesBH/F+r9QEnQcxuAhYH784F1gfvF+L+mx+P/4t4ElAWeuwt4KXA/DdgN3Imfi/IiYGt/+ws9Nn6M2b3ArwPrM/DFWkHg8flAJZDb3/F6ed+lgWNeE9j2BqAJmD7AMULj6feYwW3YSwwbgNuArMDrNgXW5wEPBm33Ij4pDqrN+jhWIf7nvduBY4Az+ttfyL99r/EM8vP1XeD9wP0DgBuB/cDHgj4vFrT9IuCkwP0s4DUgK/B4SaAd8oC9wHVAJnBh4PEF+D/uv0HgM6hFy1haUL5OiHwdeP4Hg2zH/v6tXwPmBL3GDfJz1vHveTxwbS/PLwwcJyVonXL3KC3qsY4zzvcqHopPNJ8GVplZbydprAE+ZWapwEmBdTn4/8hTgY7xse/g/8PgnGsG6oA/Oudq8f/RSgfYX2h8++l+MsQ5QKMLjNEF/oT/D/mxAY4Xut+Off4usO0N+PG/Zw5wjG7xDOWYfXjc+fF0rwETg97j+ECP0lfxbQqDbLM+3u9ufAJ71Dn3L+fcM0PYX1/xDMZDwCzzw0E+DtwH/B24wMxygDoXyLwBJwAvBe4vwn95dfSCnAdcH2jnGuCfzo/LXgrkOef+4PzP48sZ2r+BSFxQvk6MfN2LPt9jP8c5A8h3znXEs3SIxzwauMQ591+D3F65e5SkRDsAGTwzSwamOuc2AP9hZr/En0FdARwbvK1zbr+ZTQa+BfxPxy6cc7vMbClwOv4v6pn4YqrzpUH3WwkMF+prf32EGryPg4DUoLjazWwDMKW/4w3EOddkZuuAAqBoCMcY9jFDtNL1/2casN45d3voRoNpMzP7FL6A7XCw8z8Zu+BYh/Bv0Gc8A3HOrTGz5fjhIPnOuXVm9hC+F+kTwF+D4u74bHSM35uO72Hq2Fd18K6D7rfRXTsaliZjjPJ1t/eXCPk6WJ/t2M+/zWH4PyKG6zTgIDPLcc7V97FNe+A9KHePIjVIfEklaPyrc24F/qe2Q0I3NLPTgP/nnLvB+Rkegl0JHG1+2qc3gZsHOvAA++vPJnwPaPA0RP+/vXt5kaOKAjD+nQyMiwhqkEBmVHQnqJCFkI0oiigaYTa+QAQjQoJEZiEuhAHBdWbhYyGKkIUYEYUQBAOiG59/gKhEyMaIKxEfuMgsrotzW+80Xd09Y7WTsb8fFHRVddW9VdV95vZ9TZB9D/+tReDcjNOYxk9s7odMRBzawj07AxxslpGjy7dwvpH5mf5yeAc4Apxv8rcArJLNpgOHyP6iAz8Ct0XE3ibdsf38pP8x4/Vm8xSvO69xTDq/1WMW2Z4TZG3/Sx37fy3/DDo3ds+QBevd55mIWGrWl4Gz9fVFcrDGjeTUTldExGURcXfdvz8irib7+r1JNuN8QW1arPYw9Cs9ImLC+YbTjuYcp4E/gEfrua6s6Q1qPrvS63Kgvme5Hnt2ijTa/ExKs72OYaNqL4KshdoLfBAR90YOkpn0DP5OCzhQSvmuWTaaJBaa1+PO1+Z7ZH4iYjEiXq/dPMZ5l2yCfg+glPI78CHwedk8G8g9ZDeRgcH9PhU5F/ZxYDCX6fAzGH7Ou3qGAqmD8Zq5idfT3seR6ZCxdJH88QXZOkFELE0Ru4OsPX4ceCgiHq7HPhsRa2Rt/cfN+43dszSLjtsus1nIAREF+AV4CzhJNiVdVfe/AFwAVsimne/JGoLDZE3HSTJInSEHUlys5/uT/KLdB2yQA2SuBZ6v+58cd76htFeBT8l+azfVfbeSNZ3rZF/Dm+v2zvQ6rr/U41fJrgm3NPu60ri+zc+kNNt7OJT2CjmI79V6LwbHHan77wC+rs9mvd7nae/ZyohrfbCmdwpYrtvGPdNN5+rIzzXkPzJ4e4rP2stD648Adw5tOzHiuNuBb4CfgbW67X7ys/YKsESdcgw4Rv7hfaNe6wM7/R1zcelrwXg9T/H6OuAzstb94IRrHPesD5MthV+RA8nPAU/Xax8Zu8lWkQ2ypnof8AlZqH8KOErO571OdkMZHGPsnuES9eZpTkTEfuBYKeXFur6HrEU5WkpZ29HMTRARBbihZJ9FbVNEPFdK2dZUTpL+O8ZrtYzdu4NdQebPcbJbwL5m2100g9IuRU3Tk81O2xQRCxHxGPD+TudF0lSM1zJ27zIWrOfPa2RfrG8j4gfgI+BCKeXLnc1WtzqYY7WuPhE57Zu27nLgdCnl/MR3SroUGK8Fxu5dxa4gkiRJUg+ssZYkSZJ6YMFakiRJ6oEFa0mSJKkHFqwlSZKkHliwliRJknpgwVqSJEnqgQVrSZIkqQcWrCVJkqQeWLCWJEmSevAX+kHjivnhXK4AAAAASUVORK5CYII=\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/curiosity-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/curiosity-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/curiosity-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 }