{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 32 - Deep Space 2 - Mars" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from AMAT.planet import Planet\n", "from AMAT.vehicle import Vehicle" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook simulates the atmospheric entry of the Deep Space 2 aeroshells.\n", "https://en.wikipedia.org/wiki/Deep_Space_2" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"MARS\") \n", "planet.h_skip = 128.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('DS2', 3.67, 36.2, 0.00, 0.096, 0.0, 0.0875, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(128.0,0.0,0.0,6.90,0.0,-13.25,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+17YcXAAAgAElEQVR4nOzdeZyT1b3H8c9vhlVGAdlBhgEKiqIiDiqKdcOlVq1rW71dXCpVa3u7WO1Vq621i631arXWjtVbq7alaq11X1q14lYHKW4gqOwDMqDs68z87h/nCYTZM5PkSTLf9+t1XkmePEl+eQK//ObkPOeYuyMiIiIiIu1TFHcAIiIiIiKFQIW1iIiIiEgaqLAWEREREUkDFdYiIiIiImmgwlpEREREJA1UWIuIiIiIpEGnuANIl759+3pZWRnV1dX069cv7nAKho5n+uhYpk+hHcvp06evcPfCeUOtkMjZqSq0zz5B7yu/6H3ll3S/r+ZydsEU1mVlZVRWVlJeXk5lZWXc4RQMHc/00bFMn0I7lma2IO4Ysi2Rs1NVaJ99gt5XftH7yi/pfl/N5WwNBRERERERSQMV1iIiIiIiaVBwhfWUKVPiDqGg6Himj45l+uhYdlyF+tnrfeUXva/8ks33Ze6etRfLpPLyci/EcUEiUvjMbLq7l8cdRzYpZ4tIvmouZxdcj7WIiIiISBwKZlYQEZG4uDu1tbWYGcXFxXGHU5Dcoa4uXNchFpFcpR5rEZF2+ta3vkXnzp351a9+FXcoBeull6BTJ/jkJ+OORESkaSqsRUTaKdFLXVtbG3Mkha9ATgsSkQKlwlpEpJ1UWGeeWbhUYS0iuUyFtYhIO6mwzrxEYS0ikstUWIuItJMK6+xRj7WI5DIV1iIi7aTCOvM0FERE8oEKaxGRdlJhnXkaCiIi+UCFtYhIO6mwzh71WItILlNhLSLSTiqsM09DQUQkH2S9sDazo83sVTMrS9p2qpnNM7OVZnaTmXVKuq+vmf3WzK4zsxvMrHO2YxYRaU4hF9a5krM1FERE8kFWC2szGwCUAAckbSsFTgZOB74FfAX4ZtLD7gPudPfLgDnANVkLWESkFQq1sM7FnK0eaxHJZVktrN39Q+ChepuHAee6+3R3/wPwa+AIADM7CNjd3V+N9n0Y+LqZ9chWzCIiLSnUwjqXcraGgohIPsj6UBB3r6t3+wV3r0natARYGF0/EpiXtO+S6Gp5RoMUEUlBp05hJERNTU0Le+afXMnZGgoiIvkgF09enADcFl0fAqysd/9aYHD9B1VXV1NeXr6tVVRUZDhMEZGgLT3WFRUV2/IV0DczkWVFVnO2eqxFJA6tzdmdmrojDmb2CaDa3WdGmxzYWG+3rsDW+o/t168flZWVGY5QRKShthTWU6ZMYcqUKQCY2YqMBJZh2czZGgoiInFqbc7OmR5rMysGvgp8L2nzEqB30j4G9ACqshudiEjTCnWMdXOynbM1FERE8kHOFNbAJcD17r45advTwO5Jt8uALYC6pkUkZ3TEwpqYcrZ6rEUkl8Uxj3Wi38GStl0OzAC6m9kIMzvXzD7h7pXAUjMbG+16HHCzu2/JbtQiIk0r5MI6V3K2hoKISD7I6hhrM9sZ+EJ082wzu4XwU+KP6u06293vjK5/DrjKzOYT5lO9Ihuxioi0VqEW1rmUszUURETyQVYLa3dfC/wmagnXRq2pxywAzstwaCIibVaohXUu5mz1WItILsulMdYiInmpUAvrXKKhICKSD1RYi4i0kwrrzNNQEBHJByqsRUTaSYV19qjHWkRymQprEZF2UmGdeRoKIiL5QIW1iEg7qbDOPA0FEZF8oMJaRKSdOnUKEyzV1NTEHEnhU4+1iOQyFdYiIu2kHuvM01AQEckHKqxFRNpJhXXmaSiIiOQDFdYiIu2kwjp71GMtIrlMhbWISDupsM48DQURkXygwlpEpJ1UWGeehoKISD5QYS0i0k4qrLNHPdYikstUWIuItJMK68zTUBARyQcqrEVE2kmFdeZpKIiI5AMV1iIi7aTCOnvUYy0iuUyFtYhIO6mwzjwNBRGRfKDCWkSknVRYZ56GgohIPlBhLSLSTiqss0c91iKSy1RYi4i0U6dOnQCoqamJOZLCpaEgIpIPOmX7Bc3saOBa4HPuPj/a1h34ObAW6A1c5u5rWrpPRCQXFHKPda7kbA0FEZF8kNUeazMbAJQAB9S761ZgmrtfDjwE/LaV94mIxK5QC+tczNnqsRaRXJbVwtrdPyQk2m3MbDBwFvBwtOkZ4CQzK2vuvqwELCLSCoVaWOdSztZQEBHJB1kfCuLudbbjb3qHAyvcfUN0f42ZzQMOA7Y2c9/8bMYtItKUQi2sIXdytoaCiEg+yHph3YghwMp629YCg4GaZu4TEckJhVxYNyKWnJ0orOvq2vtMIiKZkwuFtQMb623rSuj5oIX7tqmurqa8vHzb7SlTpjBlypQ0hiki0ri2FNYVFRVUVFQkbvZNf1QZE0vO7tw5XG5t8EwiIpnX2pydC4X1EsKZ48l2BqoICbyp+3bQr18/KisrMxKgiEhz2lJYJxeSZrYiI4FlRiw5u0uXcLllS+sDFRFJl9bm7FyYx/pZYHA0RRNm1hkoBZ5v4T4RkZzQwYaCxJKzE4W1eqxFJJdlvbC27WfBGIC7LwMeAY6Ktk8GHnD3Jc3dl8WQRUSaVVQUUqm74wU2bUWu5Gz1WItIPsjqUBAz2xn4QnTzbDO7xd2rgQuBn5vZvsBAIHmgXXP3iYjEzswoKiqirq6O2trabSsx5rtcytmJMdZbtoQp9zRLiIjkIiuU3pXy8nLXGGsRiUuXLl3YunUrmzZtomvXrik91symu3t5y3sWjrbk7E6doLY2FNeJQltEJNuay9m5MMZaRCTvJXqpO8g461hoOIiI5DoV1iIiaZA4gbGmpibmSAqXCmsRyXUqrEVE0qCDzQwSC80MIiK5ToW1iEgaqLDOvERhvXlzvHGIiDRFhbWISBqosM68kpJwuW5dvHGIiDRFhbWISBqosM68nj3D5Zo18cYhItKUtBTWZnZ8Op5HRCRf5VNhna85e5ddwqUKaxHJVSmtYmBmxwFXEib9TxTlRcAAoHt6QxMRyR+5WFgXWs5WYS0iuS7V5cHuBn4KzADqom0GnJXOoERE8k0uFtYUWM5WYS0iuS7VwvoDd7+h/kYzm52meERE8lKOFtYFlbMThfXq1fHGISLSlFTHWF9oZqc3sv2kdAQjIpKvcrSwLqic3adPuFyxIt44RESakmphfQjwOzOrTWp1wG8yEJuISN7I0cK6oHL2oEHhsqoq3jhERJqS6lCQa4ApQCWQ+PYw4IvpDEpEJN/kaGFdUDl78OBwuXRpvHGIiDQl1cL6ceAhd99h3Sszuz19IYmI5J9OnUI6zbHCuqBytnqsRSTXpVpY/xn4HzP7Z9I2A04F/jttUYmI5JlEj3VNTU3MkeygoHJ2osdahbWI5KpUC+vvA6OAswGPthlhjtS8S9IiIumSo0NBCipn9+8P3bqFkxfXrNk+S4iISK5I9eTFXwBD3b3M3YdHrYzQ+yEi0mEVFYV0WldX18KeWVVQObuoCEaPDtfffTfeWEREGpNqYb3I3dc2sn1gOoIREclXiZ7qRM91jii4nL3HHuFydl7OxC0ihS7Vwvqi+hvM7AzgxvSEIyKSn3K0sC64nL377uFy1qx44xARaUyqhfWuZvZ9ADPbycx+B/wOeC3tkYmI5JHEEJAcK6wLLmePGxcuKyvjjUNEpDEpnbzo7p82s4PM7GbgWGAZsC/Q7gVmzWwMcDHwHuFkmwp3/4+ZdQd+DqwFegOXufua9r6eiEg6JXqsE2Otc0Eh5uwDDwyXr70GdXVh3LWISK5otrA2s082cdcCYC5wPTAMOBn4VjtjuRv4jLsvMbNS4ElgDHAr8IS7TzWz44DfAme287VERNIqF4aCdIScPWRImHavqgrmzt0+NEREJBe01GP9c0KiXAU0dqr7nYSpmwbR/iQ9Btg5ur4R6Glmg4GzgK9F258BHjCzMnef387XExFJm8RQkJh7rDtEzj7wQHjwQXj1VRXWIpJbWvoGuAYY5e7DkqZqqt/KgP9KQyz3A78zs52BLwBfBw4HVrj7BgB3rwHmAYel4fVERNImF3qs6SA5e+LEcPn88+l8VhGR9mu2sHb3x9x9eUtP4u73pyGWi4BNhJNq1rn7A8AQYGW9/dYCg+s/uLq6mvLy8m2toqIiDSGJiLROqicvVlRUbMtXQN90xNBRcvYxx4TLJ54A9+b3FRFJh9bm7FRXXsyk7sCfCIn5RjObR1gpbGO9/boCW+s/uF+/flTqNHERiUmqJy9OmTKFKVOmAGBmKzIWWObElrP32QcGDgzjrN98M9wWEcmk1ubsXDqf+h7gAXe/BvgZIWEvJZxVnmxnoCrLsYmINCtHhoJkU2w52ww+9alw/a9/Teczi4i0T04U1mbWF9jX3VdFm64F1hDOZB8cTd+EmXUGSgGNrBORnLJxY+io7d69e8yRZF4u5Owzo3lG7r5bw0FEJHe0u7A2s/3M7B4z+46ZlbXxaT4CNpnZkKRtK4E3gEeAo6Jtkwk9JEvaGq+ISCZs2LABgJ122inmSJpXKDn7yCPDtHsffAAvvpjuZxcRaZt0jLG+AJjs7l8ws7PM7EN3rz/GrlnuXmdmJwNXmdl0YADwXXdfY2YXAj83s32BgcCUNMQsIpI27p43hTUFkrOLi+FLX4Kf/QxuvBEmTcrEq4iIpMa8nb+hRT/19Y8WCSgCSuJYGbG8vNx18qKIxGHTpk10796drl27smnTppQfb2bT3b08A6E19loFk7OrqmD4cNi6FWbPhtGj0xSciEgzmsvZKQ8FMbNPmdl/RdcPIMyZugRCL4aWGxeRjiaXe6sLOWcPHgxf/GIYY/2jH8UdjYhIioW1mf0YuA84BcDd/w18yczOykBsIiJ5IVcL646Qs6+4Ajp3hnvvhbffjjsaEenoUu2xPhjYDXghadvjhKmWREQ6pPXr1wO5V1jTAXL28OFw/vmh1/qqq+KORkQ6ulQL66ej6ZWSB2afnsZ4RETyzkcffQTArrvuGnMkDXSInH3FFdCtW5jTevr0uKMRkY4s1cL6AzO7HNjdzD5vZg8AXwOuS39oIiL5YcWKsAhX375pWZk8nTpEzh48GC6+OFy/+GKI1uoREcm6lAprd/8z8E+gGPgSsAo4xt1/nYHYRETyQq4W1h0pZ19xRSiwX3klTL8nIhKHlOexdvdXgFcSt82s2MyOcvd/pDUyEZE8kauFNXScnN2rF1RUwAknhCL72GNh7Ni4oxKRjqbZwtrM/g6UtPAcfYEVQEElaRGR1sqVwrqj5+xPfxrOPRfuvDMsef7vf0MHWGFeRHJISz3Wb0T7vAvUAecDM4DkWf2HAzl3xo6ISLYsW7YMgH79+sUciXL2TTfBtGnw1ltw5ZXwy1/GHZGIdCQtFda/IKzOuArAzEa6e4MJjczsoUwEJyKSDxYsWABAaWlpzJEoZ5eUhDmtDzwwFNlf/jLss0/cUYlIR9HsyYvuvjqRoCMDzKw4eR8zOxw4IAOxiYjkhURhXVZWFmscytlBeTlceGGYHeS002DlyrgjEpGOItWTF+8HZprZNGAtMAY4FtC0/CLSIdXU1LB48WIAhg4dGnM0DXTYnH3ddfDSSzBjRiiun3oKunSJOyoRKXSpTrf3NHAM4cSX3QlTN53u7j/NQGwiIjmvqqqKmpoaBg0aRLdu3eIOZwcdOWf36AF//zsMGgTPPw9f+ILmtxaRzGvLdHtVwJWJ22ZWZGb7u7vWuxKRDue9994DYPjw4TFH0riOnLN32w0efRQOPxzuuw923hluvx2KUl0aTUSklVIqrM3szkY29wNqgZPTEpGISB558803ARibg5MmK2fDfvuF4vqYY8I0fJ06wW9+o+JaRDIj1dQyClgILEhqXYGC7/kQEWnMW2+9BcDee+8dcySNUs4GJk2Chx6Cbt3CIjIXXgjucUclIoUo1aEgX3H3d5M3mNkewNlpi0hEJI/kco81ytnbHH00PPwwnHhiKK5LS8MKjSIi6ZTqyYvvNrJ5M/CV9IQjIpI/ampqtvVY52JhrZy9o8mTYepUMAuLx9x7b9wRiUihSXWM9Twg+Qe0YmAAcHc6gxIRyQdvvvkm69evZ+TIkbEvZ94Y5eyGTjoJrr8evvMd+OIXYd06+OpX445KRApFqkNBngb+yPZEXQd86O5z0hWQmR0ETALeA6YB64GfE+Zg7Q1c5u5r0vV6IiJt9dJLLwFw8MEHxxxJk5SzG/Gtb8HGjaHX+oILwgIy//M/oSdbRKQ9Ui2s73H3fyVvMLPdzewT7v5ee4Mxs3OAke6ePDXU/wFPuPtUMzsO+C1wZntfS0SkvV5++WUAJk6cGHMkTVLOboRZGF/dpw9cdFG4vmQJ3HgjdO4cd3Qiks9SnRXkuEa2fQjc195AzOww4L+A7ydtGwycBTwcbXoGOMnMytr7eiIi7eHuPP/880BO91grZzfjggvgT38KKzLeeisce6yWPxeR9mlVYW1mF5vZKuAyM6tNbsBKQqJurxuAWcDNZvaEmR0CHA6scPcNAO5eA8wDDkvD64mItNnbb7/N4sWLGTBgQM5Ntaec3Xqf+xw89xwMGADPPgsTJsB//hN3VCKSr1pVWLv7LcAngRuBEfXabu7eWK9Iq5nZ7sB+wO3ufjHwD+AJYCjhSyDZWmBwe15PRKS9HnvsMQA+9alPUZRjq40oZ6dm4kSorIT994d582D8+HBi4+zZcUcmIvmm1d8G7v4G8D13X1CvLTWz0nbGsRfwUfQaALcARjjhZmO9fbsCW+s/QXV1NeXl5dtaRUVFO0MSEWlacmHdFhUVFdvyFZD2KUWUs1Oz227wwgvwjW+E1RnvuQfGjIETToB//lMLyoh0dK3N2ebNZAszuwR4zN3fiW5fQZiuaYfdgCPc/fC2BmtmnwUq3L1X0rYZwF+Ac9x9dNL2ucDV7v7H5OcoLy/3ysrKtoYgItJqixYtYtiwYXTp0oVly5bRq1evlh/UDDOb7u7l7Y1LOTs95s+Hn/0M7roLNm0K28aNg29/Owwd6dIl1vBEJGbN5eyWeqyHA8nfGOOBgwk/Jw5PagPbGeMbQE8zG5C0rQaoAgabWXcAM+sMlALPt/P1RETa7O6778bd+cxnPtPuojrNlLPToKwMbrsNFi2Ca66B/v3DuOsvfSms2HjppfBuY0vviEiH12yPdYOdzcYDb7n7lnrb90z0kLQ5ELPHgEfd/ddm1ouQiPcH7iFMGfWImX0K+KK7n1X/8bnc+yEihaO2tpbdd9+d999/n8cee6zNQ0GSpavHupHnVc5Og02bwuwhN9wA0UKbABxyCJx3Hpx6KvTsGV98IpJd7emxrm/vRhL07sDItgaX5EvAwWb2LeBq4MzojPILgc9EP2keD0xJw2uJiLTJX//6V95//31GjBjB0UcfHXc4LVHOToNu3eCcc+CNN+Cll0Ix3aMHvPginHtu6NE+5ZSwXPr69XFHKyJxalWPtZkdShiXNwWof4bJrsCv3L29J8O0Sz71fohIfnJ3JkyYwPTp07n11lu58MIL0/K86e6xVs7OvHXr4C9/CSc5Pvfc9pMbd9oJTjwRPv95OOaYcFtECktzObu1Ky+uAv5AGJt3SL37tgC/a3t4IiL5YerUqUyfPp0BAwZw9tlnxx1Oc5SzM6ykJPRWn3suVFXB/ffDn/8ML78ceq6nTg0nOR5yCBx9NEyeHE6A1MqOIoWt1WOszawncIa752RCzvfeDxHJbRs2bGCPPfZg0aJF/O53v+O8885L23NnYoy1cnY85s8PPdn33x/mxk7+it1pp7AAzcEHh7mzJ06EvmmfaFFEMq25nJ3SyYvNvMAkd5/W7idqh0JN0iKSG77xjW9w8803M27cOCorKykurj+LXdtl6uTFZl5POTsLVq4Mc2A//XRY1fG99xruM3o0HHQQlJeHtu++Gj4ikuvaPBTEzP4OlLT0/ITFAvq3LTwRkdz2+OOPc/PNN9O5c2fuuOOOtBbV6aScnVv69IEzzggNoLoaXnklnAD50kvw2mswZ05of/hD2KeoCPbcM6wCWV4eLlVsi+SPlsZYv0sYq7eEsKJWYww4IZ1BiYjkitmzZ3PmmWcCcM011zB+/PiYI2qWcnYO69cvnNh44onh9tatMHMmvPoqTJ8eho68806Y0u+tt8ICNQDFxTBqFOy9945t+PBQiItI7mhp5cW+wFp339zsk5j1cfeV6Q4uFR3hZ0URya6lS5dy6KGH8v7773PKKadw//33U5SBSiaNKy8qZ+e5DRvCtH6VlaHYnj4d3n4b6uoa7rvTTrDXXtsL7b32gj32gCFDVHCLZFKbh4K4+4pWvsYp6CxzESkgVVVVHHHEEbz//vvsv//+3H333RkpqtNJOTv/7bRTGHN90EHbt23cCLNmwZtvbm9vvRVmI3nttdDqP8fo0bD77tvbHnuEbSUtDRQSkXZp7XR7jTKzUuBrwMUoSYtIgZg1axYnnHACH3zwAePGjePJJ5+kR48ecYfVbsrZ+al7dxg/PrRkK1eGAjtRbM+aFZZaX748LMH+n/80fK4hQ0KhPWoUjBwJn/hEuBw5Mix6IyLtk3JhbWbFwEnA+cDRwIeEMX0iInnviSee4HOf+xxr1qyhvLycJ554gj59+sQdVpspZxeuPn3gsMNCS/bxx6HATm6zZ4dZSZYsCe2f/2z4fAMHbi+061/uumt23pNIvmt1YW1mw4GvAOcAfYAFwInAk8DEjEQnIpIlW7Zs4corr+T666/H3Tn99NO566672ClPp2NQzu64evduOJwEoLY2zLP97rvw/vuh0E5czpsHy5aFNq2RiRh7925YcJeVwbBhsNtuWvhGJKHFwtrMziD0dBwJvAfcANwFXO7uT0S7vZSxCEVEMmzGjBmcd955zJgxg+LiYn7wgx9w+eWX5/yY6sYoZ0tTiou3D/uor7YWFi/eXmgnF93vvx96wSsrQ6uvqAgGD95eaDfWunfP+NsTyQmt6bHeSFgC9yXgK+4+B8DM2r+yjIhIjFavXs1VV13FLbfcQl1dHcOHD+fee+9l4sS87tBVzpaUFRdvL4KPPHLH+9zhww8b9nIvWBBaVVUoyhcvbry3G6B//x0L7d1227ENHAid2nXWl0huaPGfsbs/AjwSnfRyvpmVAQ8B2374MTPzdCzhKCKSBZs3b6aiooJrr72W5cuXU1RUxDe/+U1++MMfsssuu8QdXrsoZ0u6mYXCd+BAOOSQhvdv2RKK6vnztxfbyW3RonBC5fLlDWcwSSgqgkGDGhbcyW3wYOjSJaNvVaTdWv33obsvBL5vZp2Ak4HdzOxHwB3AMUBFZkIUEUmPLVu2cM8993DNNdewYMECACZOnMitt97KuHHjYo4uvZSzJVu6dIERI0JrTG0tLF26vdBeuDCcQJno5V68OIztTpxY+eqrTb/WgAE7FtuDBm1vgweHy759NY+3xKfZBWJafLDZJ4CvA+e7e6xn+GixARFpyqpVq6ioqOCmm26iqqoKgL322osf//jHnHTSSZhZrPGla4GYVryOcrbkpC1bQvGdXGzXb1VVjS+UU1+nTqEATy62GyvA+/fX8BNpmzYvENMSd38P+G8zm96e5xERyYQZM2Zw++23c/fdd7Nu3ToAxo4dy2WXXcaZZ55JcXFxzBFml3K25KouXbaPv25KTU0Y651cbC9dur1VVYXLjz7a3vvdHLNQXLdUgA8cqCEo0npp+VvN3f+QjucREWmv1atXM3XqVCoqKpg+fXv9eNRRR3HJJZdw7LHHxt5DHTflbMlHnTqFBW6GDIEDD2x6v82bw9CSRKFdv/BOtOrqUKh/+GHji+kk69On6cI7ueXp7JySRvoRRETy3oYNG3jkkUf405/+xGOPPcaWLVsA6N27N1/84hc5//zzGTt2bMxRikg2dO3acu83wNatoahurvheujQU6StXbl/psjm77NK6AnyXXUKPuRQeFdYikpfWrFnDU089xYMPPshDDz3E+vXrATAzjjjiCM4991xOO+00umsCXRFpROfO20+CbE5tbejdrl9w1y/Gly2DNWtCe/fd5p+ze/fWFeB9+qgAzzcqrEUkb8yfP5+HH36Yhx9+mOeee46tW7duu+/AAw/kzDPP5IwzzmDw4MExRikihaS4ePt0g/vt1/R+7mF8d0sF+NKlsGEDfPBBaM3p3Hn7iZjJbeDAHW8PGKDVL3NFThXWZtYFeA34b3d/zsy6Az8H1gK9gcvcfU2cMYpI9qxbt45//etfPPPMMzz99NO8lfQ7bFFREZMmTeLEE0/k9NNPZ0RTc31JRilviwRmoYe5Tx9obuSZO6xd23QBnijCly2DVau2n6jZkr59my++E7dLStL3nqWhnCqsge8CZUm3bwWecPepZnYc8FvgzDgCE5HM27x5M6+99hr/+Mc/eOaZZ3jllVeoqanZdv/OO+/Msccey4knnsjxxx9P3759Y4xWIsrbIikwC2Osd9kFdt+9+X03bgwFdmKoSf2x34nry5fDihWhvflm889ZUtJy8a1hKG2XM4W1mR0MLAU+jm4PBs4Cvhbt8gzwgJmVufv8WIIUkbT6+OOPeemll5g2bRovvvgi//73v9m8efO2+4uKijjwwAOZPHkykydPZuLEiXTt2jXGiCWZ8rZIZnXvDsOHh9ac2tpQXDdXfCfaunUwd25ozencefsQGA1Dab2cKKzNrAdwurt/28yuijYfDqxw9w0A7l5jZvOAw4D5sQQqIm1WW1vL7Nmzee2113jllVeYNm0ab7/9doP99txzT4488kgmT57MYYcdRq9evWKIVlqivC2SO4qLtxe6LY0DX726+eI7cXvVqrAc/aJFLb9+S8NQhg4NJ4l2hAV5cuUtXgb8rN62IcDKetvWAo2elVRdXU15+fZFcKZMmcKUKVPSGaOItJK7M2/ePF577bVt7fXXX9+2SEtCly5dOOCAAzjkkEOYNGkSBx98MLvuumtMUWdXRUUFFRXbVhXPxzEt7crbytki2WcGvdLFTYkAACAASURBVHqFNmZM8/umexhKp06hwC4r294LP3z49tsDB+b2UvStzdntWtI8HczsWKCru/89uj0fOBsoB85w9wOT9n0d+KO7X1//ebQ8rkg8amtrmTt3LjNnzmTmzJlMnz6dyspKPvroowb7Dhs2jPLycg444AAmTZrE/vvvr6EdZG9J83RJR95WzhYpDC0NQ6mqgoULw2VzEvOPJxfcI0bAqFGh9eiRlbfTKhlb0jxNLgH2S1oJrTfwEHBDdD3ZzkALH42IZMqqVat44403mDlz5rbLt956i40bNzbYt3///kyYMGFbKy8vp3///jFELRmgvC0iQOuHoWzaBAsWwPz5MG9eaMnXV6yAOXNCa8yQITB6dMM2fHhujfHOhcL6LCC5y+pl4NvAq8B3zay7u280s85AKfB8DDGKdCjr1q1j9uzZvPPOO9vaG2+8wYIFCxrdv7S0lH333Zd9992XcePGMWHCBIYOHdrhlw4vYMrbIpKSbt3CLChNzYSybt32Qjtx+f774STL996DJUtCe/bZHR/XqVMY1rLPPrDvvqHts08YWhKH2Atrd69Ovm1mtUC1uy80s0eAo4BHgMnAA+6+JIYwRQrSqlWrmDVrFrNmzdqhiG6qgO7WrRtjx47dVkTvs88+7LPPPvTuXb+TUgqZ8raIpFtJSZj/u7E5wGtqwnCSRI92clu4MIztfvNNuPfe7Y/p3x/Ky2HSJDj00HC9W7fMv4/YC+sWXAj83Mz2BQYCOrNFJEVbt25l/vz5zJkzh7lz5zJnzhzmzJnDrFmzqGpi0FuXLl0YPXo0e+6557Y2duxYRo0aRaeOcFq3tIfytoikVadOYbz1iBFw3HE73rduHbz9NsycCW+8EdrMmWHc92OPhQZhDPdRR8EZZ8Cpp4Z5xDMh9pMX00UnwkhHVldXx+LFi3conhOX8+bN22GRlWTdunVjzJgxOxTQe+65JyNGjFABnUX5dvJiOihni0imuIfhJK+8Ai+8EFrSwr307Alf/zpcfnmYKzxVuX7yooi0wubNm1mwYAEffPDBDm3u3Lm89957bNq0qcnHlpaWMnr0aEaNGrXtcsyYMQwbNozi4uIsvgsREZHMMts+u8iZ0bqvy5fDgw/CPffAtGlw7bXw97+HojudvdcqrEVyhLtTXV3doHBOtMWLF9PcL0z9+/dn9OjROxTQo0ePZuTIkXRvy5/kIiIiBaJ/f/jqV0N76ik49tgwbOQvf4GvfCV9r6PCWiRL3J2PP/6YBQsW7NCSi+f169c3+fiioiKGDRvGiBEjGD58OCNGjGDEiBGMHDmS0aNH07Nnzyy+GxERkfzgHk5ynDYNHn8cHnkkbC8qgnHj0vtaKqxF0qSuro5ly5Y1KJyTW/2VB+vr2bMnI0eO3FY0J7fS0lI659JknSIiIjlm3Tp4912YPRtmzYIZM6CyMgwFSTZxIvzkJ2G2kHRSYS3SSlu2bGHx4sVNFs2LFi1iy5YtzT5HSUkJw4YN26ElF8+atk5ERKRpmzaF3ucFC3ZsibmvFy9u/HG77goHHADHHAPHH9/0fNrtpcJahLAs97Jly1i0aFGTbdmyZc2OcQbo27dvg8I5ufXu3VuLpoiIiNSzZcuOy6BXVTV+fcWK5p+nS5ewIuMee4Q2dixMmBBOZMzG168Kayl4iZMCmyuaq6qqmpySLqGoqIghQ4Y0WTSXlpbSo0ePLL0rERGR3FVbCytXQnV1KIYbu6yuDkM0WlMwJ3TqBEOHwrBhDVtZWWhxTnalwlrymruzatWqZovmxYsXs3nz5hafq3///gwdOrTJNmjQII1xFhGRDmPrVli9Gj7+GFatCi1xvbnLlSvDZSpLpRQVhWXIBw2CwYNDa+x6v37xFs4tUWEtOW39+vUsWrSIhQsXNlk4NzeTRkLv3r2bLZqHDBlCt2ysdSoiIpIlmzaFwnjVqu2XTV1P3pZoLZxv3ywz6NMH+vYNxXDiMvl6375hGrx8KJhbS4W1xCYxi8bChQubbCtXrmzxeUpKSpotmnfbbTdKSkqy8I5ERETSwz0Uts0VwC1db8WPtc0qKgqrFPbuDb16NbxsalufPuFkwY64gG8HfMuSLevWrdvW25zcFixYwMKFC1m8eDFbt25t9jm6dOnC0KFDKS0tbbJw7tmzp04IFBGRnFJT03wx3Joe5Lq69sXQufP2Yrdnzx0vm7qeuOzdG0pKQnEtrafCWtokMYtGc73NH330UYvP079/f0pLS5ts/fr1o0j/q0VEJIvcwzCKVHuIk7e1YpRii3r0SK0Yrn+9W7fszIQh26mwlkYlVgmcN28eH3zwQYPLBQsWtNjb3LVr12aL5qFDh2qpbRERSbu6Oli7tu1DKFavDtO/tUdiGEWqxXDi+i67hB5nyS8qrDuwTZs2MX/+/EaL53nz5rF69epmH9+/f/9t08w11dusIRoiIpKqLVtCcZtKD3Hy9TVrUpuRojFdu6ZeDCdfLylRb3FHpMK6wG3ZsoV58+YxZ84c5s6dy5w5c7a1JUuWNPvYkpISRowYwfDhwxk+fPi26yNGjKCsrIyddtopS+9CRETySV1dKHBXrgzzEyfPNNGaAnnDhvbHsPPOqRfDyZeaKEraQoV1AXB3lixZwrvvvrtD4TxnzhzmzZtHbW1to48rLi5m2LBhDYrmxGWfPn3U4ywi0sHV1oY5iVes2F4ot3T50UftO/GuuLhtY4qTh1EUwtRtkn9UWOeZ6upq3nrrrQZtzZo1je5vZpSVlTF69Ogd2qhRoygtLaVTR5wLR0SkA3MPhfKHHzbeli8PxXGiUE51oY+ExLRrianXUuk53mknDaOQ/KSqKkfV1NQwa9Yspk+fzowZM7YV0MuXL290/z59+jBmzJgGBfTIkSO18ImISAewaRMsWQKLF4clopOL5WXLdiycWzj3vIFdd92+2EdTl8nXO+ocxiL6Z58Dkovo6dOnU1lZycyZM9m4cWODfUtKShg7dmyD1r9/fw3bEBEpUGvWhIJ58eLtxXP9261YT2ubXXaBAQMatoEDw0p4yYVy794qkkVaK2f+q5jZqcAvgV2Ae4DvuHuNmfUFfgysAjoDl7l7in9r55b169fz6quv8sILLzBt2jRefvnlRpflHj58OOXl5YwfP569996bsWPHUlpaqgJaRGLXkXJ2NmzZAgsWwAcfwLx5DS8//rjl5+jUCYYMCW3w4FAkN1U4a6ZTkczIicLazEqBk4HTgb2A3wCLgOuB+4DvufurZnYBcA3wP3HF2habNm3ihRde4Mknn+SFF17g9ddfp6amZod9ysrKKC8vZ//9999WTO+6664xRSwi0rRCz9mZ4h56l2fNCm327HD5/vuhx7m5cczdu8Nuu21vQ4Y0vN6/v1bJE4lbThTWwDDgXHevAaab2VjgCDObBuzu7q9G+z0MvGtm17p7GtY0ypy5c+fy+OOP8+STT/Lss8/uMKyjqKiI8ePHc+ihhzJp0iQmTZrEwIEDY4xWRCQlBZez0+3jj2HGDHj9dXjjje2F9Lp1je9fVATDhsHw4aGNGLHjZf/+OplPJB/kRGHt7i/U27QE2Bk4EpiXtN+SaBhEOfB81gJspXfeeYf77ruP++67j7fffnuH+8aNG8exxx7LkUceycSJE9l5551jilJEpH0KJWeny8cfw8svhyL69ddDQT1/fuP79u0LY8bs2D7xCSgt1Sp7IoUgJwrrRkwAfgFMAeqfjrEWGFz/AdXV1ZSXl2+7PWXKFKZMmZLJGAFYvnw5d911F3fdddcOxXSvXr047rjjOO644zjmmGMYNGhQxmMRkfxRUVFBRUVF4mbfOGNJg7zJ2e3lHormF1+EadNCq9ePAoShG/vuC/vtB+PGwZ57wh57hMJaRPJPa3N2zhXWZvYJoNrdZ5qZA/WnxugKNDgRpl+/flRWVmYjRNydZ599lt/+9rc8+OCDbI3mLerVqxcnn3wyn/3sZznqqKPo0qVLVuIRkfyTXEia2YqYw2mzfMjZ7bVhAzz3HDz+ODz2WDihMFnXrlBeDgccAOPHh2J69901k4ZIIWltzs6p//ZmVgx8FfhetGkJMDrpfgN6AFXZjw7q6ur429/+xk9+8hOmT58OhPHSJ5xwAueffz7HHXecimkR6TByPWe3x8aN8MgjcO+98OSTYY7ohN69YdKk7W3//UNxLSKSU4U1cAlwvbtvjm4/DVyQdH8ZsAXIajeHu/Poo49y6aWXMmvWLCD0tnzta1/jvPPOY7fddstmOCIiuSInc3Z7/PvfcOut8MADO55ouP/+cPzxoU2YoOWyRaRxOVNYm9nlhOTb3cxGAIcD/wKWmtlYd38LOA642d23ZCuud955h29+85s8/fTTAJSWlnLppZdy7rnn0l0TgYpIB5WrObstamvh/vvhxhvhlVe2bz/gADjrLDjjjDAvtIhIS3KisDazK4Ef1ds8293vNLPPAVeZ2XygBLgiGzHV1tZy/fXXc9VVV7FlyxZ69erF1VdfzUUXXaThHiLSoeVizm4Ld3j0Ubj8cnjzzbCtVy84//zQRo2KNz4RyT85UVi7+7XAtU3ctwA4L5vxLF26lDPOOIMXX3wRgPPPP5+f/vSn9OnTJ5thiIjkpFzL2W2xaBF89avhhEQI091973vwpS9Bjx7xxiYi+SsnCutc8tprr3HyySdTVVXF4MGDueOOOzjuuOPiDktERNLk97+Hb3wD1q4NPdRXXw0XXADdusUdmYjkOxXWSZ5++mlOOukkNm3axCc/+Unuv/9++vXrF3dYIiKSBps3h4I6MRXtKafAr38NWmZARNJFhXXkqaee4qSTTmLz5s2cc8453HbbbRpLLSJSIDZsgM98Bp55JkyNd+utcM45WiZcRNJLhTXwxhtvcNppp7F582YuvPBCbrnlFoqKiuIOS0RE0qCmJvROP/MMDBgQTljcf/+4oxKRQtThC+uVK1dy4oknsm7dOs466yx+/etfY+rCEBEpGJdeCk89Bf37w7PPwpgxcUckIoWqQ3fLujsXXHABCxcu5MADD+SOO+5QUS0iUkCefBL+93/D8uJ//auKahHJrA5dWE+dOpX777+fkpIS/vSnP9FNp4SLiBSMTZvgoovC9WuvhUMOiTceESl8HbawXrduHd/5zncAuOGGGxg+fHjMEYmISDr9/vfwwQew557w7W/HHY2IdAQdtrC+7rrrqKqqYsKECZx3Xs6vZSAiIinYuhWuuy5cv/pq6Nw53nhEpGPosIX1UUcdxb777suNN96oGUBERArM1q1w9tkwcSKcdlrc0YhIR2HuHncMaVFeXu6VlZUpPaaurk5FtYjEzsymu3t53HFkU1tydlu4a65qEUmv5nJ2h64qVVSLiBQ2FdUikk2qLEVERERE0kCFtYiIiIhIGqiwFhERERFJg4IrrCsqKuIOoaDoeKaPjmX66Fh2XIX62et95Re9r/ySzfelwlqapeOZPjqW6aNj2XEV6mev95Vf9L7yiwprEREREZE8UzDzWJtZNbAA6AusiDmcQqLjmT46lulTaMdymLv3izuIbErK2akqtM8+Qe8rv+h95Zd0v68mc3bBFNYiIiIiInHSUBARERERkTRQYS0iIiIikgYFU1ibWXczu9nMfmJmvzGzXeKOKV+Z2almNs/MVprZTWbWKe6YCoGZdTGzmWZ2eNyx5DszO8jMLjGzk82sb9zxSHq1Np+b2afM7BYzu9XMTs92nKlobV41s1FmVmtmHrX9sh1rW5jZtUkxz2xin7z5vADM7NWk95RoxzeyXy8zW5O0zylxxNsSMzs6ek9lSdtaXTvl6ufXxPtqdR2T7v9zhVQw3Qo84e5Tzew44LfAmTHHlHfMrBQ4GTgd2Av4DbAIuD7OuArEd4GyuIPId2Z2DjDS3a+MOxbJmBbzuZmNBn4MlLt7nZn9y8xmu/tbMcTbrBTz6leATwF1wFZ3n5GtONvKzHoQTg47OtrU4KTUfPq8AMxsDPAycAmwNdr8O+Cfjex+PvBlYG10+9mMB5giMxsAlAAH1LurVbVTrn5+jb2vNtQxaf0/VxAnL5rZYGAe0NvdN0R/mawG9nL3+bEGl2fM7FDgZXeviW7/nHAcPx1vZPnNzA4G9gCuAs529+fijSg/mdlhwPeBo70Qkpc00Np8bma3AR+5++XR7SuAT7j7OTGE3azW5tWoSLgDuMDdF2c/0rYxs29GV29z901N7JM3nxeE4szdFybdHgb83N0/V2+/rsAjwEXuPjfLYabEzIqAWmC4u89PpXbK5c+vkffV6jomE//nCmUoyOHACnffABAdzHnAYXEGlY/c/YXEP8bIEmBhU/tLy6LenNPd/c64YykANwCzgJvN7AkzOyTugCTtDqd1+fzIaHvCu9Fjc04KefVM4AhgkZn9X5Q7cpqZGfB54JfAh2bW1C/FefN5ASQX1ZFTgL81susJwIHAHDN7NJeHprl7Xb1Nh9P62ilnP7/67yvFOibt/+cKpbAeAqyst20tMDiGWArNBOC2uIPIc5cBP4s7iHxnZrsD+wG3u/vFwD+AJ6IeBykcrc3n9ffLp5zfaF519xuBXYCTgMmEnrSc5sFBhKEgNwP3mtmnGtk1nz8vgOOBR+tvdPcHgJ6EYrQM+HvUg5oPUqmd8vnza7KOycT/uXz58FviwMZ627qyfVyUtIGZfQKodvdGT0SRlpnZsUCluy+PO5YCsBfhp8g3otu3AAZ8Jr6QJANam8/r75cXOb+lvOrute7+MHAscKqZDcpqgG3k7h9H5z38FPhmY7uQh58XQNQLvcnd1zR2f/THxb8Ivbp7EHqw80EqtVNefn6tqWPS/X+uUArrJUDvett2BqpiiKUgmFkx8FXge3HHkucuAe40sxVmtgIYCjxkZpfGHFc+6kTSCdfuvhGYC/SJLSLJhNbm8/r75XzOTyWvuvs7hF9lSjMdV5r9isZjzrvPK8lJwN9b2sndPwSmkj+fWSq1U959fqnWMen6P1cohfWzwGAz6w5gZp0JB+b5WKPKb5cA17v75rgDyXNnAeOSWhXhDGQNr0ndG0DPekM/aoB3YopHMqO1+fxpYPek26PJwdkY6kk1r64HZmcwnkyoA15vZHs+fl4JrSqsI7XAfzIYSzqlUjvl4+fXljqm3f/nCqKwdvdlhLNyj4o2TQYecPcl8UWVv8zscmAG0N3MRpjZudHPKZIid69298WJRki61U39pChNc/fZwOOEKZQws15AFxoZ9yj5q7l8bmHu8sSX+22En24TJ9AdSegtzUlN5NVRZnaNmQ2M9jnbzPpH1w8Fnnf31TGG3SIzG2xmZ5lZUfQ5fAe4Irovbz+vhOhktm7Jw/mS35eZnW5mI6LrexCGq70bT7TNi447hCF0LdZO+fL51X9f0bZG6xgLMvp/riCm2wMws97Az4H5wEDgf9x9XaxB5SEzuxL4Ub3Ns919TBzxFBozm4+m22uzaKzjTUAloWfl9ujnOykgTeVzM5sO/NTd74/2+wIwPnrY8+7+UBzxtqSpvEqI/R3gTHd/xcweAg4h9A4+7u5/yG6kqTOzvQh/8K4BpgE3Rn8Ek6+fVzILC6HsFp3klti27X2Z2S3AFwjzWz8H3NLI7BuxM7OdCXHeClxDiLO6udopHz6/xt4XYfhHo3VM1Duf0f9zBVNYi4iIiIjEqSCGgoiIiIiIxE2FtYiIiIhIGqiwFhERERFJAxXWIiIiIiJpoMJaRERERCQNVFiLiIiIiKSBCmuRJpjZUWb2kpkd3sJ+Y83sL2b2/TS+9i/S9VwiInEws33M7L505sYmXucRM/t8Jl8jU8xsspm93NL3jOSPTnEHINISMxsD/Bj4CJgA7AN80d3vyfBLLwL2bcV+C4HeQHE6XjRaVnZrOp5LRATAzCYAPwU+CTwFzCMsCNIZ+I27P5mBl30P6EmacmMzbgXeysQTR/n4VMLKg3OAl6O7ugKnAQe4+/x2vERrv2ckT6iwlnzwN+Cz7j7TzIqA/8vGi7r7HDNb2Yr91pjZ0jS+9GHAv9L4fCLSwbn7a2b2Z2APdz8hsd3MPg08aGaXuPutaX7NDWnOjQCY2YnAm4mC1t0fS/drJLj7VmCqmX0deM7dr0yKo00rD0ZLaB/h7lPd/V0zW5GmcCUHaCiI5DQz2xUYDawHiJaKvQrI1pKhrV2aNp1L2B4BPJvG5xMRAaipv8HdHwV+APyvmQ3MwGumNVeb2W6EHupsa3DsgGeAJak8iZl1Bf4AdE/anHNLoEvbqbCWnObuHwHvAk+Y2UHRtgXAv82sl5n9yMz+bWYnmdkyM5tlZuMTjzezM83sCjN7zsx+krS93Mx+YGZTzeyvZtYz2t7NzH5tZleZ2R2EnzFTYmaXRjF91swmmNnfouf7URTjVDPrF21fbWZX13uKLu6+OXqu75jZt83sRTP7WcoHUESkZXcQhoScAWBmI83sGjP7nZn9IypmMbOuUS67wsxeMLPPJZ7AzC4ys++b2StmdnFTL9RYTjaz06OcebaZzTWzn5nZYDP7vyifvmBmR0ZPcRawG3CZmZ0cnQvzjJl9Oek1zjezH5vZH8zsLjPrYWY9zexqM6s0syPN7F0z+8DMytpywKJjcnjUo42ZnWpmN5nZL83sCTMrid7DTWb2cPR9szQ6xnsCZ5jZt5Oesnf0XbTWzM5v4jWLzOw6M5sSvYab2YOJz0dyhLurqeV0A/YCFhD+qr8T6B9tLyYkqY8IY916AU8DbxP+aJwIXBnt25vQ6/1pYBfgz0nP/yLwi+j6/wJTousDgFpC8mwpxt8Ten36E8bi7RRt7wo8CTwEDCZ8IWwFrgZ2AvaPbif2T4wfT7zvW6PrOwPfi/uzUFNTy98GnA0sbuK+D4FfRXn1r0BxtP2PwH3R9RuAU6Pr5wLV0fUzgS9E10dHuXrv6PbvgR9E15vKybsCq4AbgQOAycD1wFXRvt8FHk6K1YGy6Hp3YBZwdnT7WOCPSfs+AvyWMPT1lOj74tTofT4B/LiVx+656LviZ8AvCOOtJyfdvwIoja6/TfhO6gpcASwGDgTOi/L+c4l4o/3nA7+M7jsTWNhEDF8BHomuD4iO8zFx/7tS27FpjLXkPHd/28z2Bq4FLgJONLPJHsZcVwNr3P0BADP7ATANGEH4Eqkzs29GT/UkIYF/GuiTtH1W9NgS4MLosbj7h2a2KIVQ+wI3A+e5+4boOTab2TJgnrtXRa+zDHjew/jDGYSE35dwEuSnCV8CABuAL5nZHOA3hD8qREQywQnF5gHAMODrZgawEuhp4cZXCZ0CAHcRhkJAyLWV9XLtIODNeq9xNo3kZHf/yMxWAw+5+78BzGwhsN7MBhF6eEsaDdp9o5ktT9p0cfS8Cb8D/hJtX0X4vvhr9BqVUZyt9axHY6zNbCqhMyfhDHdfaGaHEArqkij/VwHvu/urwKvRYxt77oej74TXmolpPLAxet8fmtlH7DikRHKACmvJaWZWDAz1cJLKN6LhGY8BFYQegPreji57AaXA3e7+x2jbjdFzXkYodG+s91r7ExLi2jaGOxqYBIwEZiZtrz/GsHbbHe51UZJNDMvq42H4C+4+z8zOBm4BvgX8F7AcEZE0MrNehF/b3ibkzZWN5McBhB7VLgDuXkvoDCB6zA/d/aXo9g6PTdJoTo44O+bKxcBlwFzgFWB4M28h+XGjgH8m3X6fMMylPw1zcQ1tn7HkdcIfIAmrzOw64D5C7Inquf77akkNTddmjwO/MbNuhM9iI6H3W3KIxlhLrts27g/A3WcClwN7NLF/F0ISew9YCpyefKeZHRht/7SZda+3fU10c0wbY30JuAm4N0p8KTGzwSSdCGNmQwgzoowmnMz4xyYeKiLSHl8m/EJ2HyE/HhL1FAPb8uMKYAtwXNL2vc1sZ+rlWjPrFHVU1NdUTm7MbcB8D9Oqbk7hvSwk9HBvewlgNWGoS9q4e13U+bFX9IfJY8A17l6Zztep95oPAw8Tflk9EzjM3Vdn6vWkbVRYSz74elR0JgwhjI1L2NXMukTXDyeMr1sF/Ak4xcxuMbPDo96E9YQE2AN4xMyOjU4U2YNQjM8AfmRmXcysL2Fsc38zKzaz75nZF5qI0aJ2FeHL54ZG7qu/f/3bJxDGAyaMAk529zXA1xt5jIhIKhr0hJrZJ4EfAl9192pC7/CHwGNmdmJ0guKhUQ/1nwmzh5xsZkcBX3b3tYRc+w0z+6GZHUboYEgMo0vOf03l5ITk3uP9gH5mthNwCNDdzEZE920lnOyX6GBJfo3bgNPMrHd0+0DgdnevIdQ8jeZiMzsoiqtrE8euiHo1U/THx6cJv1L2Jgwx3B0YCuyUdGJk/V7xLY3Ev2NQjYwXiT6LXYB/EHrlU/mDQ7JEQ0EkHwwF3jazRwk/kxlwQdL9NcAVZraK0Lv7dQB3f9rMvkH4OfFkwsl/bwGY2WeAXxO+KO4EvuvubmafJUyF9DYwlfCT3jjCggqfJQwx2WFhGjM7gHBSTinhp8G/AT+M4nkSOAgoM7ORhJ6UgcBnzWw223tvzgZ6uHtFvff+fxYWdqghnCwkIpKyKI+cRSj+bgeqCSfA9QQ+5e4vA7j7lig/3k7IdQ8RxlYD/DehU+IPhF/oEjNx3EHIfRcShqxd6O7Lo8LxIGCkmd3TVE42s9MJ44q/Ymaz3X0J4UTK6wlF9R+AzwB7Ax8AdxN6179kYbaovYDjzewxd/+rmY0CHjCzFwm/el4R9ayfAQyM3t/bhMVyBpnZXlGcnyeMG/9b0nHrTMj9+wJDosLbCWObT42e8z/A88BrhOEtjxNy+iOEEyb3NbOT3T3xvFMJi/WsiwrxgcDpZjaP0BNN9Pj6azYsJJyceTLQDSgys2cIJzBmawpaaYHps5B8ZmEZ2N+7e1mWXu+77q7lxkVECkw0LKWru+fkAl1m9l3CDC3zo9vdCL82XJc4N0fix4O6EwAAIABJREFUp6Egku8aG2aRmRcKP5m+no3XEhGR7IlOzhyfw0V1N+BKoCxpc19gvYrq3KKhIJK3ovFtXyD8tHeGu9+X4Zd8P/qJUkRECksnd/9N3EE0xd03WVhW/VfRie1VhGErVzb/SMk2DQUREREREUkDDQUREREREUkDFdYiIiIiImmgwlpEREREJA0K5uTFvn37ellZWcqPq66upl+/fukPqADpWLWejlXr6VjB9OnTV7h7hzoIytnpo2OyIx2PhnRMGmrPMWkuZxdMYV1WVkZlZeoriZaXl7fpcR2RjlXr6Vi1no4VmNmCuGPINuXs9NEx2ZGOR0M6Jg2155g0l7M1FEREREREJA1UWIuIiIiIpEGHL6ynTJkSdwh5Q8eq9XSsWk/HSlKhfy8N6ZjsSMejIR2ThjJ1TApmgZjy8nLX+CERyUdmNt3dy+OOI5uUs0UkXzWXszt8j7WIiIiISDqosBYRERERSYOOW1j//vcwfjzcdFPckYiISAvWzlhL5fhKZp8zO+5QRESa1HEL6w8/hBkzoKoq7khERKQFdRvqWDdjHRve3RB3KCIiTeq4hbWIiOSP4nDhdYVxwr2IFCYV1iIikvOsyMKVunjjEBFpjgrrLVugQKYcFBEpWNG3ldcqX4tI7uq4hXVR9NZvvBEuvjjeWEREpFnqsRaRfNBxC+tTT4XDDgvXn3gi3lhERKRZVhwKa42xFpFc1nEL65Ej4YEHwvWPPoo3FhERaV7i20o91iKSwzpuYQ3Qq1e4XL0aamvjjUVERJqUGAqiMdYiksuyXlib2dFm9qqZlSVtO9XM5pnZSjO7ycw6Jd3X18x+a2bXmdkNZtY5bcEUF0PPnuHkxdWr0/a0IiKFImdytqbbE5E8kNXC2swGACXAAUnbSoGTgdOBbwFfAb6Z9LD7gDvd/TJgDnBNWoPadddwqeEgIiI7yKWcrZMXRSQf/D979x1fZ1n/f/z1yd5J26Qj3YPuAqWlIEM2whcExLJ/KirUBe6F4kJUEFQERQzDyZKtqKiILBFKSxld0N3SmabNTjOv3x/XSXNyspOT3Cc57+fjcT/OOfd9n/v+5M7pp59c57qva0ALa+fcbuCJiNUTgY8555Y7534P/BI4CcDMjgZmOOdeCe37F+BqM8uMWlDDhvnH/fujdkgRkaEgpnK2htsTkUFgwLuCOOeaIl6/4JxrCFu1Hdgaen4ysCls3+2hpwujFpBarEVEOhQrObt5VBC1WItILIvFmxePBO4IPR8LlERsrwAKI99UXFzMwoULDy5FRUXdO1tzYV0SeRoRkf5TVFR0MF8B+UHH0wcDkrMPDrfXoBZrERl43c3ZSR1tCIKZTQOKnXNvhFY5oCZit1SgPvK9BQUFLFu2rOcnHTXKP+7a1fP3ioj00pIlS1iyZAkAZrY34HB6ZSBztqX4wrqpTk3WIjLwupuzY6bF2swSgU8AXw9bvR0YFraPAZnAjqideOxY/7gjeocUERnqBjpnJ6T6/66aalVYi0jsipnCGvgycLNzrjZs3b+AGWGvJwF1QC+apjtQGPqGcvv2zvcTEZFwA5qzE1L8f1euTl1BRCR2BTGOtTU/DVv3DWAFkG5mU8zsY2Y2zTm3DNhpZnNDu54B3Oacq4taQGqxFhHpUKzkbEsNdQVRi7WIxLAB7WNtZtnA/wu9vNzMfoH/KvH7Ebuudc7dE3p+EfBtM9uMH0/1m1ENqrnFWoW1iEgrsZSzLdF8ad/kh9w7OEqIiEgMMeeGxtdqCxcudL26ebGiAnJyID0dqqrAlKxFZGCZ2XLnXPSGER0EepOzn09/nqYDTRxfdTyJGYn9FJmISOc6y9mx1Mc6GNnZfpKYmhrYvTvoaEREpAMHu4NoZBARiVEqrAEOOcQ/rlsXbBwiItKhgzcw1g6Nb1pFZOhRYQ0wfbp/fOedYOMQEZEOacg9EYl1KqxBLdYiIoOAJokRkVinwhpaCmu1WIuIxKzmFmt1BRGRWKXCGmBGaD6D1auDjUNERDqkFmsRiXUqrAHmzIGkJN9iXVkZdDQiItIO9bEWkVinwhogNRXmzgXn4PXXg45GRETa0Tx2dVO1CmsRiU0qrJstWOAfly8PNg4REWlXYpYvrBsrGwOORESkfSqsmzUX1r2ZvVFERPqdCmsRiXUqrJsdfbR/fPHFYOMQEZF2JWaHCusKFdYiEptUWDc79FDIy4PNm/0iIiIxRS3WIhLrVFg3S0yE977XP3/uuWBjERGRNlRYi0isU2Ed7sQT/eMzzwQahoiItNXcFaShoiHgSERE2qfCOtypp/rHv/8dGtUiIiISS9RiLSKxToV1uLlzYcoUKC6Gl14KOhoREQmjwlpEYp0K63Bm8IEP+OePPx5sLCIi0kpSThIAjWUqrEUkNqmwjnTeef7xkUegSbN7iYjEiqQRvrCu31cfcCQiIu1TYR3pmGNgwgTYskWjg4iIxJDk/GQA6veqsBaR2KTCOlJCAnzkI/75PfcEG4uIiBykwlpEYt2AF9ZmdpqZvWJmk8LWpZvZbWb2QzP7lZnldGdbv7n8cv/4yCNQWtrvpxMRiVWxlLOThyWDQcP+Bpoa1FVPRGLPgBbWZjYKyAIWRWy6HXjROfcN4Ang193c1j+mTIGTT4aaGrjrrn4/nYhILIq1nG2JRtLwJHC+uBYRiTUDWlg753bjE+1BZlYIXAr8JbTqaeAcM5vU2bZ+D/aLX/SPt9wCdXX9fjoRkVgTizlb3UFEJJYNeFcQ51zk93cnAnudc9Wh7Q3AJuCELrb1rzPPhNmzYft2eOCBfj+diEgsirWcfbCw3qPCWkRiTyzcvDgWKIlYVwEUdrGtleLiYhYuXHhwKSoq6ltUCQnwpS/559//PtQriYtI9BQVFR3MV0B+0PH0QKA5O7UwFYDa7bU9DFtEpPe6m7OTBi6kDjmgJmJdKtBcyXa27aCCggKWLVsW3cg+9CG44QZYtw7uvhs++cnoHl9E4taSJUtYsmQJAGa2N+BweiLQnJ06IVRYb1NhLSIDp7s5OxZarLcDwyLWZQM7utjW/5KT4frr/fPvfQ8qKwfktCIiMSzQnJ02Pg2AA9sOROuQIiJREwuF9X+AQjNLBzCzZGAC8FwX2wbG4sWwcCHs2gXf+c6AnVZEJEYFmrNTx4darLeqxVpEYk8Q41hb81MA59wu4EnglND6U4FHnHPbO9s2YAEnJMAdd/jHW26BaHc3ERGJYbGWsw8W1uoKIiIxaKDHsc4GmjsqX25mBaHnnwLONbNvAv8HLAl7W2fbBsaCBfD5z0NTk+93XVU14CGIiAy0WMzZaRNCXUG2qiuIiMQec84FHUNULFy40EX95sVw1dW+S8iaNfCxj/mbGUVEosDMljvnFgYdx0Dqbc52zvFi7os0VjRyTPExpOSn9EN0IiId6yxnx0If68EhIwP+9CdIS4N77oE//CHoiERE4o6ZkTEzA4DqtdUBRyMi0poK656YOxduvdU/v/JKeOmlYOMREYlDBwvrNSqsRSS2qLDuqSuugE9/Gmpr4dxzYePGoCMSEYkrarEWkVilwrqnzODnP4f3vQ/27oXTT4cdAzOstoiIQMYsX1hXrdKN5CISW1RY90ZSEjz4IBxxBGzYAKecAnv2BB2ViEhcyJqfBUDFsgqGyg34IjI0qLDurdxc+Oc/Yd48WLsWTj0Vdu8OOioRkSEvbWIayQXJNJQ0cGCzht0TkdgRlcLazP4vGscZdEaMgH/9C2bOhLfeguOOg02bgo5KRKRTgz1nmxnZR2YDULG0IuBoRERa9KiwNrMzzOxFM1tvZhtDy2bgkf4JbxAYNQqefRbmz4f16+GYY+D114OOSkRkSOfsnEU5AJQvLQ84EhGRFkk93P8PwI+AFUBTaJ0Bl0YzqEGnubg+7zz4z3/g2GPhN7+BCy8MOjIRiW9DNmfnvMcX1qXPlQYciYhIi54W1hudcz+NXGlma6MUz+CVkwN/+xt84hPw+9/DRRfBihVw/fWQmBh0dCISn4Zszs49LhdLMSpfq6S+pJ7kEclBhyQi0uM+1p8ys8XtrD8nGsEMemlp8Nvfwi23+GL6hhvgtNPg3XeDjkxE4tOQzdmJGYnkHpcLDvY/sz/ocEREgJ4X1scCd5lZY9jSBPyqH2IbnMzgc5/zI4aMHOm7hhx6KDz6aNCRiUj8GdI5e9ipwwDY/08V1iISG3paWF8HLAEOAaaELd+NblhDwMknw5tvwplnwv798MEP+mnQy3WjjYgMmCGds4efORyAvX/ei2vUeNYiEryeFtZ/B55wzm10zm0JLZuBO6Mf2hAwahT89a9w662Qmgp33QVz5sCTTwYdmYjEhyGds7MOyyJtahr1e+opfV43MYpI8HpaWD8AXGNm7w1bTgCu6YfYhgYzuPpqWLYMFi3y/a3f/3645BLYuTPo6ERkaBvSOdvMGHnBSACKHyoOOBoREbCeTAdrZq/iv1IsBZrfaMBo51xa9MPrvoULF7ply5YFGULXGhvhttvgm9+E6mrIyoJrr4XPf963aItIXDKz5c65hf1w3CGfsytWVLD8iOUkjUjimO3HkJCqCYVFpH91lrN7moFuAsY75yY55yaHlknA+X0NMi4kJvoieuVKOPdcqKyEr3/ddw957DHowR85IiLdMORzdtbhWWQemklDSQPFj6rVWkSC1dPCeptzrr35Y0dHI5i4MXkyPP64Hzlk9mzYsAHOP9/P2vjcc0FHJyJDx5DP2WZG4ScLAdjx6x0BRyMi8a6nhfWnI1eY2QXALdEJJ86cdhq88YbvHlJQAC+/DCeeCGecAa+9FnR0IjL4xUXOHnXZKBIyEyh7rozKlZVBhyMicaynhfVwM/sWgJllmNldwF3Aq30NxMxmmdkvzewLZna7mR0eWp9uZreZ2Q/N7FdmltPXc8WUpCS46irfan3ddZCdDf/4ByxYABdfDOvWBR2hiAxecZGzk3KSGPPRMQBs/dHW/j6diEiHelRYO+fOAv5lZrcBrwPTgcOA9mb26qk/AD90zv0MuAG4P7T+duBF59w3gCeAX0fhXLEnOxu+9S3YuBG+9CV/M+ODD8KsWfDJT8IOfcUpIj0TTzl7/FfGY0nGngf2UL2+eiBOKSLSRqeFdcQQTe81s/cCKcAWYB3wHWAi8O0oxDILyA49rwFyzawQuBT4S2j908A5ZjYpCueLTfn5cPPNvqX64x/3NzT++tcwdSp89atQrJtzRKR98Zyz0yakMerDo6AJtly/pb9PJyLSrk6H2zOzl/HJsxRo6mg3YIxzrk/jxZnZ74CpwJnAFcBWIBW4yTk3Nmy/laF1vwt//6AYbq831qzxLdmPPOJfZ2T4FuwvfxnGjAk2NhGJimgNtxfvObtmYw1LZy7FNTgWLF9A9vzsrt8kItJDfRlu7zrgEOfcxLChmiKXScBlUYjz08ABfN+/SufcI8BYoCRivwqgMPLNxcXFLFy48OBSVFQUhZBiwKxZ8PDDsHQpnHWWH//6pz/1I4tcdRVsVX9CkcGoqKjoYL4C8qN02LjO2elT0hl79VhwsOFLG+jJPA0iIp3pbs7u0QQx/cnM8oFz8Yn5a6HnhwMXOOeOCtvvNeA+59zN4e8fsi3WkVasgOuvh0cf9a+TkuDCC+Fzn/MzO4rIoNNfE8T0p1jN2fX76nll2is07G9g9kOzGbl4ZNTPISLxLZoTxPSnPwKPOOeuo+VGmJ3AsIj9soH4vZNv/nzfLWTlSrj0Ut8H+7774Kij4Nhj4aGHoKEh6ChFZOiLyZydPDyZyddPBmDdVeuo31c/UKcWEYmNwjrU8nGYc640tOp6oBx/w02hmaWH9ksGJgCaRWXOHLj3Xj+KyFe/Cnl58NJLvvV6yhT43vdg27agoxSRISjWc3bhJwvJPS6X+t31bPjyhoE8tYjEuZgorIF9wAEzGxu2rgR4E3gSOCW07lR8C8n2AY4vdk2YADfeCO++C7/8JUyf7gvq734XJk2Cs8+GP/9ZrdgiEk0xnbMtwZhx1wws1dj1m13sfWLvQJ5eROJYn/tYm9l84EvACnwC3dzL4xyGvxlmOTAKeN4595yZDQN+DGzGT8N7jXOuzdRacdPHuitNTfCf/0BRETz2GNSHvgYtLISPfhQuvxymTQs0RBFpbSD7WMdTzt72021s+NIGkoYlsfD1haRNSOvX84lIfOgsZ0ejsP41cK5zbrSZXQo85pyr6dNBe0GFdTuKi+F3v4M774R33mlZf/zxvsi+4ALIygouPhEBBrywjpuc7Zocb53zFvv+uo+c9+Rw+LOHk5ASK1/Uishg1d83L14FLAg9fwBIjsIxJRoKCvx412vXwrPPwoc/7MfBfuEF+NjHYPRoX2A//7y/CVJE4kHc5GxLMGb+diYpY1Mo/18573zqHQ3BJyL9qseFtZmdaWaXhZ4vwo+Zuh3AOdfknCuPcozSV2Zwwgm+9XrXLrjrLjjuOKiqgt/+1m+bMcP31d61K+hoRSSK4j1np+SnMPexuSSkJ7Drnl1su1k3dYtI/+lRYW1mPwAeAj4A4JxbCnw49HWiDAbZ2X6q9Bde8N1DvvEN3/963Tr4+tdh3Dg47zx4/HGorQ06WhHpA+VsL+fIHGb+fiYAG7+2kT0P7wk4IhEZqnraYn0MMA54IWzd3/FjmMpgc8gh8IMfwJYt8OSTvqAGeOIJ+MAHfFeRK66AZ56BxsZgYxWR3lDODhm5eCSTfzgZHKy5dA0lf4+cIFJEpO96Wlj/KzRuaXgntcVRjEeCkJTkp0t/7DE/bN+NN8Jhh0FpKdx9N5xyCowfD1/8Irz6qvpjiwweytlhJnx9AuO+MA5X71h1/ir2P7s/6JBEZIjpaWG90cy+Acwws4vN7BHgM8CN0Q9NAjF6tJ9w5vXX/eyO3/ymn3Bm50742c/8tOnTp8O3vw2rVgUdrYh0Tjk7jJkx9SdTGfOJMTQdaOKts99i39P7gg5LRIaQHhXWzrkHgGeARODDQClwunPul/0QmwRtzhy4/npYvx5efhk++1kYNcq//v73Ye5cv8/3vgerVwcdrYhEUM5uy8yYfvt0Rl8+mqaqJt466y2KHysOOiwRGSKiMY51InCic+7f0QmpdzSO9QBpbPQT0DzwgO86si+stWfOHD+l+gUXwKxZwcUoMsgM8DjWytn4Ma7Xf34922/bDgkw484ZjPnYmMDiEZHBo9cTxJjZn4GuZhDJB/Y6507ufYh9F3SSjkv19f7Gxocegkcfhf1h/RWPOAK+8hVYvNj34RaRDkWrsFbO7hnnHJu/u5kt120BfB/syT+YjCVYoHGJSGzrLGd3VfG8GdrnbaAJuBI/DW54NpwMDI9CnDLYJCfD+97nl1/9Cv79b/jTn/xQfa+9BpdcAtdc44f3u+giPwqJiPQn5eweMDMmf28yKWNSWHfVOrbesJXqt6uZ9YdZJGYmBh2eiAxCXbVY54b2KQ29vs459+129nvCOXdu/4XZtVho/ZCQAwfgD3+Am27y42M3O+II31Xkgx+EadOCi08kxkSxxVo5u5f2Pb2PVYtX0VjWSNb8LOY8Oof0SelBhyUiMajXU5o758qaE3TIqFD/vPCDnwgs6nOUMnSkpcGVV8KaNfC3v8GHPuQnpnntNT8JzSGHwPz5fgztt98OOlqRIUM5u/eGnzqcI14+grSpaVSuqGT5/OXs/fPeoMMSkUGmp8PtPQy8YWZ3mNlNZvYk8C/g1uiHJoNeYiKceSb8/vewe7fvh33ZZb7Ifv11uPZamDkT5s3zI4usWqUxskWiSzm7BzJnZrLg1QWMeP8IGkobWHnuSjZ8dQNN9U1BhyYig0RPh9v7F3A6sBeYgR+6abFz7kf9EJsMJenpfjbHP/4R9uyBv/wFPvIRyMvz42V/97t++L5Zs3zB/frrKrJF+kg5u+eShyUz94m5TLlpCiTCtpu2seL4FVS/Ux10aCIyCERjuL0EYL5zbnl0QuqdWOuvJ91UV+dHFnn44bbD902d6kcVWbwYFiwA0536MjQN8HB7ytndVPbfMlZfvJrad2tJSE9g6k1TKfxUoUYNEYlzvR5ur50D3dPO6gKg0Tl3Xi/ji4rBkKSlC/X18Nxzvsh+9FEoDpu0YeJEX2B/8INw9NEqsmVI6a/CWjm77+r317P+s+vZ/cfdAAw7dRgz7p5B2oS0gCMTkaD0+ubFdhwCbAW2hC2pQKAtHzJEJCfDqafCHXf4KdT/8x/4zGf8NOtbtsBPfgLHHOO7i9x6K5SWdn1MkfimnN1HycOSmfWHWcx5eA5JI5LY//R+ls5eyrZbttHUoL7XItJaT1usZzjn3o5YNxO43Dn39WgH1xODpfVDeqGpCV56CR55BB580BfdABkZvhX7kkvglFN8YS4yCPVji7VydhTV7a5j3VXrKH7Yf5uWdUQWM4pmkL0gO+DIRGQgRa3FOjJBh9QCV/QmMJFuSUiA446Dn/3Mt1w//DCcdBJUV/sRR848E8aMgU9+0nclaWwMOmKRmKCcHV0po1KY89Ac5v5lLqkTUql8rZLli5az/gvraShrCDo8EYkBPW2x3gSEvyERGAX8wTl3ZVQCMjsaOA5YD7wIVAE/BiqAYcDXnHPlke8bjK0f0kfvvAMPPAD33w9r17asLyz0E9FcfDEsWqT+2BLz+rHFWjm7nzRUNrD5O5t595Z3oQmSC5KZfP1kxnx8DJaonCMylEXz5sUi4D5aEnUTsNs5906fo/TH/ygw1Tl3bdi63wBPOeceNLMzgI845y6JfO9gT9LSB87BW2/5IvuBB2DTppZt48f7Yf7OP9+3eidqmmKJPf1YWCtn97OKFRWs/+x6yl4sAyDz0Eym3TKNYScNCzgyEekv0Sys3+ucez5i3Qz8Hebr+xjkCcC3gNNcKCgzKwQ2AcOcc9VmlgSUAXOcc5vD3z9UkrT0kXOwdKkvsP/0J9ixo2VbQQGcc44vsk85BVJTg4tTJEw/FtbK2QPAOUfxQ8Vs+MoGarfWApB/Xj6TfziZzFmZAUcnItEWzVFBzmhn3W7goR5H1dZPgTXAbWb2lJkdC5wI7HXOVQM45xrwSfuEKJxPhiIzOOoo3x972zZ4+WX46ldh2jQ/fN/dd8NZZ/ki+5JL4KGHoLIy6KhF+oty9gAwM0ZeOJJFaxcx6fuTSMhIYO/je3l17qusuXwNNZtrgg5RRAZItwprM7vKzEqBr5lZY/gClOATda+FWlDmA3c6564C/g08BYwPHT9cBVDYl/NJnEhI8EX2jTf6/thvveWnTj/sMKio8K3aF14I+fn+ZsivfMW3cm/cqFkfZVBTzg5GYnoik66dxFHrjqLwk34imd2/283S6UtZd/U6anfVBh2iiPSzbncFMbNDgY8At0ZsqnPO7exTEGbnA0XOufzQ63SgGLgO+KBz7qiwfV8D7nPO3Rx+jIkTJ7qCgoKDr5csWcKSJUv6EpYMZRs3+pkeH30U/ve/toX0sGGwcKFfFizwy8SJuhFSoqaoqIiioiIAli9fvsU5Nymax1fODl7Nhho2f3czu+/dDQ4S0hIYc8UYxn95PGkTNcGMyGDS3Zzd0z7Wyc65+nbWT3DObe1lrJjZhfgknRe2bgXwJ+CjzrnpYevXAd9xzt0Xfoyh1F9PBtiePb5f9rJlfnn1Vb8u0vDhcMQRflmwwD9OnapiW/qsH/tYK2fHgMq3Ktn0rU2UPOEb8y3JGHnZSCZ8bYL6YIsMQp3l7KQu3vhl4G/OudWhVV81s8hhFQw4Cd+3rrfeBHLNbJRzrvkrygZgB1BoZunOuRozSwYmAM/14VwirY0cCWef7Rfwrdfbt/sCe9kyeO01WL7c99F++mm/NMvNhfnzWwrtBQvgkEN8NxSRAaacHZuy5mUx7/F5VL5ZydYbt7LngT3s/t1udv9+NyPOGcG4z44j76Q8TH+kiwx6nRbWwGQgL+z1EUAmsIuW4ZsMGN2XIJxza83s78Bi4JdmlgekAPcCZwKnAE8CpwKPOOe29+V8Ip0yg3Hj/PKBD/h1zcX28uV+aS62d+2CZ5/1S7OsrLbF9owZGupPBoJydgzLOjSL2ffOZvL3J7Ptpm3s/M1OSp4ooeSJEjLmZDDu6nGM+n+jSMxUrhAZrHraFeQIYKVzri5i/eywFpLeBWKWD/wcWIZv4bjTObfazIbhJxvYjP/P4BrnXJthHOLpa0WJITt3ti60X3sN3n237X4ZGTBvnr9x8tBDW5bc3IGPWWJOP3YFUc6OYbW7atlZtJMdv9pB3S7/K0rKS2L0x0dTuKSQjOkZAUcoIu2J5jjWH3HO/S5i3QxgunPuL30Ls2+UpCVm7N4NK1a0bt3esqX9fSdObCm2mx+nTlXrdpzpx8JaOXsQaKprovjhYrbftp3yl1smqcw9PpcxHx9DweICtWKLxJA+F9Zmdjz+68MlQFHE5uHArc65CX0NtC+UpCWmlZTAm2/65Y03/OPKlVDbzvBbGRkwd25Lq/bcub7fdmGh+m4PUdEurJWzB6/yV8vZcfsO9vxpD03VTQAkZicy8tKRjPn4GLIXZqsvtkjAolFYzwN+j++/tz9icx1wr3Puur4G2hdK0jLoNDTAunWti+033/QT27QnPd23Zh9yiJ/wJvxRRfeg1g+FtXL2INdQ3sCeB/ew8+6dVLxScXB9+vR0Rl4yklGXjCJjhrqKiAQhKl1BzCwXuMA5d1c0g4sWJWkZMvbt85PZNBfca9b4Ary4uOP3NBfdkQX3tGkwdqyK7hjXH11BlLOHjqpVVey8eye7/7ib+uKW0ROzjshi5CUjGXnRSNLGa1xskYEStT7WnZzgOOfci30+UB8oScuQV1YG69f7IjvysbOiOzUVJk/2RfbUqa2XyZMhJWXgfgZpV3/1se7kfMrZg1BTQxOl/y5l9/272fvYXhrLGw9uyzk6h/zz8sk/L18t2SL9rNeFtZlYuOf+AAAgAElEQVT9Gcjq6vjAHOfcyN6H2HdK0hLX2iu6162DDRvan+ymWUICjB/futgOL8CzswfuZ4hj0SqslbPjR+OBRvb9bR977t9DyZMlNB1oOrgtY2YG+R/wRXb2wmwsQX2yRaKp1xPEAG8DpcB2WsZAbXN84OzehycifZab2zL1eqSKCl9ghy/r1/vHbdv8iCVbtsAzz7R9b0FB22K7+XVBgWadjD3K2XEiMS2RgvMLKDi/gMaqRvb9Yx97H99LyZMlVK+tZuuPtrL1R1tJKUxhxP+NYPiZwxl2yjCScrv6b19E+qKrFut8oMI5187QBa32G+GcK4l2cD2h1g+RXqirg82bWxfbzcvGje2PWtIsK6ttK3fz47hxGjKwB6LYYq2cHeea6psoe6GMvY/vZe/je6ndFvZRSITcY3IZfsZwhp85nKzDstSaLdILA9HH+oqgb5BRkhaJsqYmP9tkey3dGzb47icdSUmBKVPaFt3Tpvmxu9Wvu5UA+lgrZ8cB5xyVKyrZ99Q+9j21j7KXyqClWzbJo5IZfrpvyc47KY+0CboBUqQ7+q2wNrMJwGeAq5xzmb0+UBQoSYsMIOf86CWRxXbz8127On5vQoIvrtsruqdM8eN4x5mBKqyVs+NbQ1kD+5/ef7DQrn239RcbaVPSyDspj2En+0I7dUxqQJGKxLaoFtZmlgicA1wJnAbsDh1nbF8D7QslaZEYUlnpu5I0F9rhj1u3+sK8I4WFvsieMaP1MnkyJA3N/qH9WVgrZ0t7nHNUr65m3z/3UfqfUkqfK201yghA+ox0hp00jLwT88g9LpfUsSq0RaBvNy+GH2QycAXwUWAEsAV4P/AP4D1RiFNEhoqsrJaZIyPV1vp+3e0V3Zs2wY4dfnn++dbvS05uv+CeMQNGjBiQH2swUc6WzpgZmXMyyZyTyfgvjMc1OipWVPgi+5lSSl8opebtGmrermHHHTsASJ2YSu6xueQel0vusblkzsnEEtVHWyRcly3WZnYBvqXjZGA9cBfwO+Abzrkv9HuE3aTWD5EhoLHRj1Sybh28/XbLsnZtxzNSAuTnt19wT53qC/IYF80Wa+VsiYam+iYqlvlCu+yFMspeKmvTop2Yk0jOe3J8sX1sLjlH5ZCYqZuWZejra4t1DX4K3JeAK5xz74QO2ve7HkVEwiUmwqRJfjnttNbbqqraL7jfeQf27vXLf//b9nhTp8LMmTB3Lsyb55fp0wdFwd1LytnSZwnJCeS+J5fc9+QC4BodVauqKPtvmV9eLKN2Sy37/7Gf/f/YH3oTZM7NJGdRDtlHZpO9KJvMuZkkJGnmV4kfPZnSfAK+FWQS8ARwgnPu6tA2c9EYXqQP1PohEqec811Hwgvu5qJ7y5b2+3MnJ/tie9681gX3hAmBjM3dT1OaK2dLv6rdXtuq0K58o7LVqCMACekJZB2R5YvtRdnkLMohbXIapjHwZRCL9s2LScB5wGXASuBu4HTnXFFfA+0LJWkRaaOmxvfdXr0aVq6Et97yy8aN7e+fne0L7COOaJlwZ9asfr9psp9vXlTOlgHRWN1I5euVlC8tp2JpBeVLyzmw4UCb/ZJGJJG9MJvsI7LJmp9F1vws0qeka0xtGTT6c7i9acDVwJXOuUDHyFKSFpFuq6z0xXZzod1cdLc3/Xt6Ohx2GBx5JJxyCpx0EuTkRDWcARxuTzlbBlR9ST0Vyypaiu1Xyqkvrm+zX2JOIlmH+yI7e342WUdkkTEzg4RkdSOR2DMQE8R82Dn3+z4fqA+UpEWkz/bsgTfegNdeg+XL/RLZup2YCO95D1x6KVx8MQwb1ufTBjBBjHK2BMI5x4EtB6hcXknFigoqV1RSuaKSup11bfa1VCNrXhZZR/iCO+vwLDLnZpKUNTSH3ZTBo98L61igJC0i/WLfPl9ov/QS/Otf8L//+dFLAFJTfXF9zTV+FJJeGujCOhYoZ0u42l21B4vsyhW+6G6vGwn4iWwy52WSdWjWwcf0aeka+k8GjAprEZFoKSuDJ5+E3/0Onn7a3xxpBp/6FNx0U69mjlRhLdJWQ1kDla+HWrZfq6TyzUqq11Tj6trWLQlpCWTMzjhYbGcemknWvCxSRqUEELkMdSqsRUT6w8aNcOON8JvfQH09nH46PPVUj0cWUWEt0j1N9U3UvFND5VuVVL1ZRdVbVVS+WUnt1tp2908emexbtef5biQZczLInJVJUq66k0jvDZrC2sxSgFeBzznnnjWzdODHQAUwDPiac668vfcqSYtIYF5/HRYt8sX1qlUwe3aP3j6YC+ve5m3lbImm+tJ6qlZWtSq2q96qorGisd39U8am+JknZ4eK7dmZZMzOIDlvyI5vL1EUlSnNB8hX8GOuNrsdeMo596CZnQH8GrgkiMBERNq1Zw/88pe+qE5JgYKCoCMaaMrbErjkvGTyjssj77i8g+uab5SseitUcK+uonp1NVVrqqjbXkfd9jr2/3N/q+OkFKa0KbYz52SSPEwFt3RPzBTWZnYMsBPYH3pdCFwKfCa0y9PAI2Y2yTm3OZAgRUTAD9f39NPwyCPw0ENQW+vHur777rgqrJW3JZaZGemT0kmflE7++/MPrneNjppNNb7IXlXV8rimmrodddTtqGP/0xEF9+iU1sV26DElX324pbWYKKzNLBNY7Jz7opl9O7T6RGCvc64awDnXYGabgBOAzYEEKiLxp6rKz+K4ciUsXQovv+yH5GseGcQM3v9+39d61qxgYx1AytsyWFmikTEtg4xpGeSf07rgPrDlQEuxvbqqpeDeVUfdrjpK/13a6ljJ+clkzM4gY1ao2A49phSmaHbJOBUThTXwNeCGiHVjgZKIdRVAYXsHKC4uZuHClu4uS5YsYcmSJdGMUUSGqqYmePdd2LDBL+vW+QlkVq2CzZvbTouekABHHQXnngsXXADTpvX4lEVFRRQVHZz8ML+zfWNUn/K2crbEGks00qekkz4lHd7fst41OQ5sPUD1quqW7iShgrt+bz1lz5dR9nxZq2Ml5iT6InuWL7YzZvvnaZPSNCzgINXdnB34zYtm9j4g1Tn359DrzcDlwELgAufcUWH7vgbc55y7OfI4uhFGRDrkHJSWwrZtsHWrH81jwwY/3fmGDbBpE9S1naAC8F08ZszwNyQuWABHH+0fs7KiFt5gu3kxGnlbOVsGO+cctdtrqV5dTfWaUNEdemwoaWj3PQlpCaTPSG/Vup0xK4P0aekkpGiWycEi1m9e/DIwP+wrk2HAE8BPQ8/DZQM7Bi40ERkUKith586Wwrm9x8rKzo8xejRMndqyzJ7tl0MOgWTduBRBeVvinpmRNi6NtHFpDD99eKttdcV1B2+UPFh4h26arHqjiqo3qlofK8lIn5beqnU7Y3YGGTMySMxIHMgfS/ooFgrrS4HUsNf/A74IvAJ8xczSnXM1ZpYMTACeCyBGERlozsH+/b5g7mrpqmgGyMyECRNg/HiYMqV1ET1lSlRboOOA8rZIJ1IKUkg5IYW8E/JarW8oa6B6bevW7eo11RzYdIDqtdVUr62Gx8LeYJA2Ka1V63Zz4a2xuGNT4L8V51xx+GszawSKnXNbzexJ4BTgSeBU4BHn3PYAwhSRaKirg717/RB1xcWtH8Of79wJu3b50Ta6IzUVxozxRfP48S0F9IQJLc/z8no8cYu0T3lbpHeScpPIOSqHnKNyWq1vrGmk+m3fsh3e0l2zroYDmw5wYNMB9v1tX6v3pBSmtPTjDruBMrkgWTdOBijwwroLnwJ+bGaHAaMB3dkiEiuc8yNmlJS0XjormktLuz5uuOxsXzB3tahojiXK2yI9lJieSPbh2WQfnt1qfVN9EzXra1q1blev9i3bzUMDRo5UkjQ8qU3rdsbsDFLHpargHgAxV1g75yaFPd8PXBlcNCJxorER9u1rWyR3tXR0w19HEhMhPx9GjvTjPXf0OHq0L5gzM/vn55WoUt4W6R8JyQlkzsokc1YmBee3jJHvmvzQgG36ca+uomFfA2UvllH2YsRIJVmJZMzMaDM8YPqUdI1UEkUxV1iLSB80NPhW4f37/bJvX+vn4UXx3r0tz3vaktwsPR1GjGi9dFY0Dxvmh6oTEZFeswQjfXI66ZPTGXHWiIPrnXPU7axr1brdXHjXF9dTsayCimUVrY+VamRMj7hpclYGGYdkkJCqfN1TKqxFYk1TE5SVtS2Mu/O6vLx35zTz3Skii+SOlvx8/5ieHt2fXUREes3MSC1MJbUwleGnth6ppL6kvk3rdvWaamq31fpp39+qopiw2ycSIX1qepvJbzJmZpCYqZFKOqLCWqQ/NDVBRUXrore7hXJZWdsJSbqruUAeNqxlGT685bGjQnnYMN9NQ0REhqTkEcnkHZdH3nERI5VU+JFKIruV1GysoeYdv5Q80Xrep9SJqQdbt7MOzSLr8CwyZmeQkKwWbhXWIh1pavItwJHFcXeW0lL//t7KyWlbGHfndW6uulqIiEi3JWUnkXNkDjlHRoxUcqCRmndq2kx+U/NODbVbaqndUsu+p1pGKrEUI3NuJlnzs8ien03W/CwyD80kKSu+Ss34+mkl/oR3q+jpUlbWt+I4K6t1AdzdIjk318/2JyIiEpDEtETfGn1o6zH+mxqaOLDhgG/dXlVN5RuVVK6opGZ9DZWvVVL5WiW72OV3Nkifnk7W4aFie0EWOUfmDOkxuIfuTyZDR2d9jrtTHPe2WwX44d4ii+PuLHl5mq1PRESGnISkBDJm+FkhOa9lfUN5w8Eiu/J1/1i1qoqat2uoebuG4gdD/bcNMmZlkHO0H8875+gcMudkDpmRSVRYy8BobOxby3FfiuPwbhU9LY7VciwiItKlpJwk8o7PI+/4lj7cTbVNVK2u8sX2ikrKXy2nckWl78e9uppd9/iW7YTMBLIXZpNzdA55J+SRe1wuSdmD8//fwRm1BOfAAT9MW/NSUtL94rgvcnN7VxyrW4WIiEggElITyJ6fTfb8lolvmmqbqHy9kvJXyil/uZzyV8o5sPEAZc+VUfZcGdtu3AaJkL0gm7wT88g7cXAV2oMjSukfdXUt4xl3tERur6rq3bnM+lYca8QKERGRQS8hNaFlWvfP+nV1xXW+0P5vOaXPlVLxagUVS/2y7ce+0M59Ty7DzxrOiLNH+K4jMTqLpArrocQ5P4rF7t2wa1fbZfduP610c5FcUdH1MSOlpPgxjJuX4cPb3oDX3pKTo+JYRERE2kgpSCH/7Hzyz84HoKGygfKXyil9tpTSZ32h3Tyb5KZrNpE6IZURZ4+gYHEBeSfkYQmxU2SrsB4MnPPdKbZt88u777YUypHF84ED3T9u8/TS7S3Nk4BELllZvvVZREREpB8kZSUx/PThDD/dT3LTUN7A/n/tp+SvJZT8tYTarbXsuH0HO27fQcrYFEZePJJRl41q1eUkKCqsY0F1NWzaBFu3ti6ew59XV3fvWJmZMHp062XUqJbHkSNbiuTcXBXJIiIiEtOScpIo+GABBR8swDU5KpZVsPfxvey5fw8HNh/g3Z+8y7s/eZfso7IZ99lxFCwuICElmDkdVFgPlP37Yf162LCh5bH5+c6dXb8/JwfGjYPx4/1jYWFLwRxeQGdldX0sERERkUHIEoycRTnkLMph8g8mU/6/cnbft5s99+2h4pUK1ly2hg1f2sCEayZQ+IlCElIHtsBWYR1tFRWwahWsXNl62b274/ckJ8OkSTBxoi+cm5fmQnr8eF9Yi4iIiAgAZkbuMbnkHpPL1B9PZfe9u9l+63aqVlax/nPr2fbTbUz98VQKLigYsJsdVVj3RW0tvP46vPIKvPyyf9y4sf19MzNh6lSYNs0/hj8fP1439omIiIj0UmJGIoVXFjLmijGU/KWEjd/YSPWqalZftJr8B/OZfsd0UgpS+j0OFdY90dgIy5fDP//pl1de8UPWhUtJgVmzYO5cv8yb5x/Hj4eEYPr7iIiIiMQDMyP/nHxGnDWCHXfuYONXNrL30b1Urqhk3pPzyJyd2a/nV2HdlaYmeO45uPdeeOwxP6V2uNmz4aij4Oij/TJ7tiYkEREREQmQJRpjPzmWEWeOYNWFq6hYWsGK41Yw/4X5ZM7pv+JaFWBHysrgjjvgF7/wo3I0mzQJ3vc+v5x0kp/2WkRERERiTtrENA5/9nBWX7Sakr+U8Mb73mDh6wtJye+fbiEqrCPV1cHPfgY/+EHLBCqTJsFll8Gll/puHhqiTkRERGRQSExPZPaDs3njlDco/1856z6zjjkPzumXc8VMp18zO9/MNplZiZn93MySQuvzzezXZnajmf3UzJL7LYi1a+GII+DrX/dF9Uknwd//7m9IvP56381DRbWISGzkbBGRbkpMT2TWvbNIyEyg+E/FlL9S3i/niYnC2swmAOcBi4EvAFcAnw9tfgi4xzn3NeAd4Lp+CeK552DRIj9U3iGHwD/+Ac88A2ecoWJaRCRMTORsEZEeSp+czrirxwGw9cdb++UcsdIVZCLwMedcA7DczOYCJ5nZi8AM59wrof3+ArxtZtc756qidvaVK+H97/et1BdcAL/5jR8eT0RE2hNszhYR6aWxV49l641bKXmyhIbyBpJyolsKx0SLtXPuhVCCbrYd2AqcDGwK22976OnCqJ28vh4+/OGWovr++1VUi4h0ItCcLSLSB6mFqeQel4urc+x/en/Ujx8ThXU7jgTuAMYCJRHbKoDCqJ3p3nthxQp/g+I992iiFhGRnhu4nC0i0kd5J/gR3fqjn3WsdAU5yMymAcXOuTfMzAE1EbukAvWR7ysuLmbhwpZGkSVLlrBkyZKuT/izn/nH730PsrJ6HbeISE8UFRVRVFTU/DI/yFj6YsBztohIH2UvygagckVlt9/T3Zxtzrm+xBZVZpYI3ABc65yrNbNrgJOcc6eHthtwILTupfD3Lly40C1btqxnJ1y3DqZPh9xc2LPHz5ooIjLAzGy5c27QdZcY8JwtIhIFVWuqeHX2q6RPS+eodUf1+P2d5exY6wryZeBm51xt6PW/gBlh2ycBdUB0svG//+0fzzxTRbWISM8NbM4WEYmCtIlpABzYcgDXGN0G5pgprM3sG8AKIN3MppjZx4BSYGfojnOAM4DbnHN1UTnpa6/5x6N6/teKiEg8CyRni4hEQWJGIkl5Sbh6R0NZQ9dv6IGY6GNtZtcC349YvdY5d4+ZXQR828w2A1nAN6N24lWr/OOhh0btkCIiQ11gOVtEJEqS8pJoKG2gobSB5OHRm8cqJgpr59z1wPUdbNsCfLxfTvzuu/5x4sR+ObyIyFAUWM4WEYmSpDxfAjeURrfFOma6ggy4pibYudM/L9RIUCIiIiLxQoV1tJWV+clhsrMhPT3oaERERERkgCRk+hK4sboxuseN6tEGk8rQ2IU5OcHGISIiIiIDKiHFl8CuboiOCjLgKir8oyaFEREREYkrlmIANNU2RfW48VtYN7dYZ2cHG4eIiIiIDKiEVLVYR1dVlX/MzAw2DhEREREZUM1dQZrq1GIdHbWhicJSU4ONQ0REREQGlLqCRFt9vX/UVOYiIiIiceXgzYu16goSHXWhGXaTozfbjoiIiIjEPkvyLdauUYV1dCQn+4lhRowIOhIRERERGUBJeUmkjEkhISO6pXBMTGkeiHPO8YuIiIiIxJWJ35zIxG9OjPpx47fFWkREREQkilRYi4iIiIhEgQprEREREZEoUGEtIiIiIhIFcV9YFxUVBR3CoKFr1X26Vt2nayU9oc9LW7omrel6tKVr0lZ/XRMV1vqwdZuuVffpWnWfrpX0hD4vbematKbr0ZauSVsqrEVEREREYpg5F90ZZ4JiZsXAll68NR/YG+Vwhipdq+7Tteo+XSuY6JwrCDqIgaScHVW6Jq3perSla9JWX65Jhzl7yBTWIiIiIiJBUlcQEREREZEoUGEtIiIiIhIFKqxFRERERKIgbgtrM0s3s9vM7Idm9iszywk6plhiZtebmQstb4TW6ZqFmNlpZvaKmU0KW9fh9Ynna9fetQqtb/MZC62P22slHdPnojUzO9/MNplZiZn93MySgo4pVphZipm9YWYnBh1LrDCzo83sy2Z2npnlBx1PkMxslpn90sy+YGa3m9nh0Tx+PP9DvB14yjn3oJmdAfwauCTgmGKCmWXi75Y9LbSq+c59XTPAzEYBWcCiiE2dXZ+4vHYdXatOPmMQp9dKuqTPRYiZTQDOAxYDc4BfAduAm4OMK4Z8BZgUdBCxwsw+Ckx1zl0bdCwx4g/Auc657aF/S/8AZkXr4HE5KoiZFQKbgGHOuerQX/plwBzn3OZAg4sBZvb50NM7nHMHQut0zcKYWQLQCEx2zm3u7PoAdR1ti4drF3mtQuvafMZC6/U5kzb0uWjNzI4H/uecawi9/jH+WpwVbGTBM7NjgJnAt4HLnXPPBhtRsMzsBOBbwGkuHgu+dphZFbDAObfWzAqAN5xzhdE6frx2BTkR2OucqwYIJadNwAlBBhULzMyAi4GfALvNrLlF6ER0zQ5yzjVFrDqRjq9PZ9uGvMhr1clnDOL8WkmHTkSfi4Occy80F9Uh24GtQcUTK0LfhC12zt0TdCwx5KfAGuA2M3vKzI4NOqAY8DBwl5llA/8PuDqaB4/XwnosUBKxrgKI2l8sg5XzjsZ/TX8bcK+ZnYmuWVc6uz66dmE6+YyBrpW0T5+Lzh0J3BF0EDHga8ANQQcRK8xsBjAfuNM5dxXwb+CpUBe9ePZp4ADwKlDpnHskmgeP18LaATUR61KB+gBiiUnOuf2h/lg/Aj6PrllXOrs+unbtaOczBrpW0j59LjpgZtOAYufcG13uPISZ2fuAZc65PUHHEkPmAPucc2+GXv8CMODc4EKKCenA/cB9wC1mdmo0Dx6vNy9uB4ZFrMsGdgQQS6y7FXgWXbOudHZ9XCfbpOUzBvqcSfv0uWiHmSUCnwC+HnQsMeDLwHzf0wzwn5cnzOwHzrkfBxdWoJIIq/OcczVmtg4YEVxIMeGPwMXOuVIzc8D9ZjbJOVcVjYPHa4v1f4BCM0sHMLNkYALwXKBRxaYm4DV0zbrS2fXRtetc82cMdK2kffpctO/LwM3OudqgA4kBlwKHhy07gCuI7y4ybwK5EV0/GoDVAcUTuNBQg4c550pDq64HyoniqCBxWVg753YBTwKnhFadCjzinNseXFSxwcwKzexSM0sI3WT2JeCbumatWUuziEHnn6l4v3aR16qjzxjo36a0T5+LtszsG8AKIN3MppjZx0LdQuKSc67YOfdu84IfiajYOVcedGxBcc6tBf6OH5YRM8sDUoC/BhlXwPYBB8xsbNi6EuCdaJ0gLofbAzCzYcCPgc3AaOAa51xloEHFADObg/+HWA68CNwS+sepaxYSdifx7cB1wC+cc8WdXZ94vXbtXStgJB18xkLvictrJZ3T56KFmV0LfD9i9VrnXNRa3QY7M9uMhttrbqH9ObAM/y3Pnc65uG2xBjCzw/A3MC4HRgHPO+ei9u1X3BbWIiIiIiLRFJddQUREREREok2FtYiIiIhIFKiwFhERERGJAhXWIiIiIiJRoMJaRERERCQKVFiLiIiIiESBCmsRERERkShQYS0iMcXMCszsz2a2xcy+HnQ8IiLSOeXtFpogRkRiipm9D/gPkAOsBaY75/YFG5WIiHREebuFCmsRiVlm9jxwsnOuIehYRESka/Get9UVRERikpmNB/4Zr8lZRGSwUd5WYS39wMy+ama/DDqOZtGIx8xONbP/mdmJUQprSDGzm6N8PAMuBG6I5nFFpDXl674zs0lmdqeZ/TpKx0s0syvNbEs39j3LzF7p7c8azdytvO2psB5kzGyWmT1qZneZ2Rtm5szs/wUdV4S/AfcHHUSYaMSzDTgsCrHEJDP7PzPba2Y7zew9ZjY+9DlzZvYLM8sL7bfIzFaa2U/MLD20bhSwK8ohXQrc7ZxrCB1fZNBRvu6VwZiv94QeU6N0vASgFJjQjX2XAbM62mhmc8zswdBn71ozGx22Ldq5W3kbFdaD0ePA95xzVwDzgd83bzCzQ83shCCCCj+3c26lc+7FIOJoT0/i6egaOufeBvZGM65o/r76eizn3N+A24EG59z/nHPbgCuBOuAV51xpaL+l+ET+FedcTejt7wee7NMPEMbMvgV8F3jazNYCM6J1bJEBpnzdQ4MxXzvnqoEd0TqXc64eWN7NfXfji/COtq8CikIvb3bOhRfSUcvdytstVFgPImY2HJgOVAE455qAbwPOzHLxSdsCiCuwc0dTN36OpgE8VxDHehAYZ2YLAJxzJcC/gbPDzpUH7Ah99prNdM6t7eO5D3LOfd85d4hzbqFzbqZz7vloHVtkoChf968YzNfRHgmiJ/F3tW9j6DGy33PUcrfydgsV1oNIaOiat4GnzOzo0LotwFJ8v6bJwBIz+yiAmV1nZt8KfQ10sB+VmS00s5+a2efNrNrM3jazL5rZeWb2qpl9wMweM7NyM/tE2PvaPV7Eua80s9+G9zULfR16m5n9yMz+Y2aLQus7PV/Y+83MPmFmq8zsNDN708x2m9nl3TjHIeHxdHHONtewHcNCX+1WmNmVEdf0u6Fr82goEXf3mrU5l5ktNrOlZna5ma0zsxs6OV57v/t24+lMqGVjFXBB2OoC4CwLdfsAzgH+HBZnOlAd9nqmmd1gZjeb2XPmu5RcZGbLzOxcM7vfzErN7HNmtiD0M5aY2f91FZ/IYKJ8HT/5uj0d/YydncfMJpvZPeZbf2/qznnaOe/fzOxeMzu2i/2Uu/uLc07LIFqAOcAW/F+o9wAjw7ZtBk4MPZ8LbAw9L8D/NT0C/xf3VmBhaNsvgZdCz1OAfcAv8GNRXgJs7+x4kefG9zH7NfDb0Po0fLGWF3q9GCgGsjs7Xzs/9+jQOb8R2vd64AAwsYtzRMbT6TnDr2E7MWwGfgJkhN63NehMuL8AAAYxSURBVLQ+B3ggbL//4pNit65ZB+cajv967xZgEXBqZ8eL+N23G083P1/XAhtCz6cCPwJqgfPDPi8Wtv85wPGh5xnAq0BG6PXzoeuQA5QB1wHpwEWh1xfi/7j/EqHPoBYtQ2lB+Tou8nVo+3e7eR07+12/CswKe4/r5ues+ff5HuA77Ww/MXSepLB1yt39tKjFepBxvlVxHj7RfBhYZWbt3aSxDviQmSUDx4fWZeH/IY8HmvvHrsH/g8E5VwdUAA8758rx/9BGd3G8yPhqaX0zxFlAjQv10QUexf+DPL+L80Uet/mY94X2vR7f//e0Ls7RKp6enLMDf3G+P92rwJiwn3FEqEXp8/hrCt28Zh38vPvwCewJ59xS59zTPTheR/F0x4PAFPPdQT4I3An8A7jQzLKAChfKvCHHAi+Fnp+D/8+ruRXkA8D3Q9e5FHjG+X7Zy4Ac59yfnP96fAU9+x2IDArK1/GRr9vR4c/YyXlOBXKdc83xLOvhOY8ELnPOfa+b+yt395OkoAOQ7jOzRGC8c24z8Fkzuxt/B3URcFT4vs65WjMbC3wN+F3zIZxze8xsGXAK/i/qyfhi6uBbw543EOou1NHxOgg1/BiHAMlhcTWZ2WZgXGfn64pz7oCZbQTygPwenKPX54zQQMu/nwnAJufcLZE7deeamdmH8AVssxnOf2XswmPtwe+gw3i64pxbZ2Yr8N1Bcp1zG83sQXwr0gXAX8Pibv5sNPffm4hvYWo+Vkn4ocOeN9JaE+qWJkOM8nWrny8e8nW4Dq9jJ7+bQ/F/RPTWycAhZpblnKvsYJ+m0M+g3N2PdEEGl2TC+r86597Af9U2M3JHMzsZ+LRz7nrnR3gI9xngSPPDPr0J3NjVibs4Xme24ltAw4chMnzfw75KAd7p53N0x05a90PGzI7qwTX7M3B42NLu3eU9OF678XT/x+EB4KPAxrD4EoHP4b82bXYUvr9os+3AcWaWGXbeTvv5iQxhytetxVO+7vBn7OQ85aH3pNA7N+Nb+3/ewfYy13LTuXJ3P1JhPfhcbWaFYa/HAk+Fntfhb9aYiR/aKdfMUs3s1ND2kWaWj+/rdzf+a5yXCH21GJJAxF/pZmZdHC/y3BZ2jMeBSuDi0LHyQudrbvns6HwdGRPaZ2zovU914xzh8XR1zvCfI1J7rReGb4XKBJ40s/eZv0mmq9/BwXMBY5xza8OW+rBTJIY97+x44XG3G4+ZpZhZUaibR2f+hP8K+mEA51wF8Hfgv671aCCn47uJNGu+3vebHwv7KqB5LNPI30Hk73lQj1Ag0gHla+ImX3f3OrZ7HnwuTcH/8QX+2wnMrLAbudvwrccfAi4wswtD7/2SmV2Lb63/d9j+yt39qT86bmvpnwV/Q4QD9gN/BH6L/yppWGj7d4B3gXPxX+2sw7cQnIVv6fgtPkn9GX8jRV3oeNX4f2hnAvX4G2TGA9eEtn+ss+NFnPtzwAv4fmtzQtsW4ls6f4Lvazg3tL7D83Xw87vQ+z+H75owL2xbR+eYFB5PV+cMv4YR5z4XfxPfL0LXovl9Hw1tPwFYGfrd/CR0nbt7zc5t52ddHDrf/cDY0LrOfqetjtVBPOPwExnc143P2q0Rry8CTopYd3M773svsBooAa4Nrfs//GftNqCQ0JBjwCfx//HeGfpZzw7635gWLdFaUL6Op3w9AXgR3+p+eBc/Y2e/67Pw3xS+jL+R/B3g06Gfvd3cjf9WpB7fUj0ceAZf1F8BfAI/nvf/b+cOcSIGojAA/w8sCoHnKlwDgwNBgsKh8HguQYLnABwBh4F7PERLUkMWNgPLZr8vqZimTSYVL39m3vQuUxvK5ztq9y9eNX88dkRVHSW56O7bebyXaRXlvLtvNjq5Faqqkxz31LPImqrqurvX+pUT8HfUa5bU7u2gFWT3XGZqCzhc3DvJ4lDaf7TYerLttKaq2q+q0yQPm54L8C3qNWr3lhGsd899pl6sl6p6S/KU5L27nzc7ra/Nhzmu5uFZTb994+cOkjx29+vKJ4H/QL0mUbu3ilYQAAAYwIo1AAAMIFgDAMAAgjUAAAwgWAMAwACCNQAADCBYAwDAAII1AAAMIFgDAMAAgjUAAAzwAWG7BoZuUIiyAAAAAElFTkSuQmCC\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/deep-space-2-mars.png',bbox_inches='tight')\n", "plt.savefig('../plots/deep-space-2-mars.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/deep-space-2-mars.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }