{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 48 - Crew Module Atmospheric Re-entry Experiment" ] }, { "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 Crew Module Atmospheric Re-entry Experiment. https://www.isro.gov.in/sites/default/files/flipping_book/LVM-3/CARE/files/assets/common/downloads/LVM3-brochure.pdf" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \n", "planet.h_skip = 126.0E3\n", "planet.h_trap = 0.0E3\n", "planet.loadAtmosphereModel('../atmdata/Earth/earth-gram-avg.dat', 0, 1, 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle1=Vehicle('care', 3735.0, 350.0, 0.0, np.pi*3.1**2.0*0.25, 0.0, 3.0, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle1.setInitialState(126.0,0.0,0.0,5.326,0.0,0.0,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle1.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle1.propogateEntry (120.0*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+17YcXAAAgAElEQVR4nOzdd5xcdfX/8ddJ7430uikkgQQisJAAAkHA0EFApFiwsApi+wqCgqjYEIGfIFJWREVFKaEYEJTeAyYgECAhIT0hYRNSNj3ZPb8/Pnc2k8m22Z2dO+X9fDw+j5m5c+feczOZM2c/87mfa+6OiIiIiIg0T6u4AxARERERKQQqrEVEREREMkCFtYiIiIhIBqiwFhERERHJABXWIiIiIiIZoMJaRERERCQD2sQdQKb07t3bS0pKGlyvoqKCPn36tHxAWaLjyW2FdDyFdCyQW8czc+bMVe6eG8FkSWNzdqpcet8yodCOB3RM+aDQjgeye0z15eyCKaxLSkqYMWNGg+uVlpY2ar18oePJbYV0PIV0LJBbx2Nmi+KOIdsam7NT5dL7lgmFdjygY8oHhXY8kN1jqi9nayiIiIiIiEgGqLAWEREREcmAoiusy8rK4g4ho3Q8ua2QjqeQjgUK73iKRaG9b4V2PKBjygeFdjyQO8dk7h53DBlRWlrqhTZeSESKg5nNdPfSuOPIJuVsEclX9eXsouuxFhERERFpCQUzK4iISLZt27YNgHbt2sUcSeFzh+rqcL9163hjERGpi3qsRUSa4Pzzz6d9+/b87W9/izuUovDSS9CmDRx+eNyRiIjUTYW1iEgTdO7cGYC1a9fGHElxKZDTgkSkQKmwFhFpgh49egAqrLPFLNyqsBaRXKbCWkSkCbp37w6osM6WRGEtIpLLVFiLiDSBeqzjoR5rEcllWS+szewYM3vFzEqSlp1mZgvMbLWZ3WBmbZKe621mt5nZr8zsejNrm+2YRURSFUthnSs5W0NBRCQfZLWwNrN+QBfgoKRlQ4FTgTOA7wBfAb6d9LJ7gTvc/VLgPeCqrAUsIlKHRGG9bt26mCNpObmUszUURETyQVYLa3dfCTyUsngY8CV3n+nudwK/A44EMLNJwBh3fyVadxrwDTPrnK2YRURqUww91rmYs9VjLSK5LOtDQdy9OuXx8+6+I2nRMmBxdP8TwIKkdZdFd4vq0r8iknuKobCG3MnZGgoiIvkgF09ePBC4Nbo/CFid8nwlMDCrEYmIpCiWwroRspKzNRRERPJBTl3S3MxGARXu/ka0yIHNKau1B7anvraiooLS0p2dImVlZZSVlbVUqCJS5Lp16wbA+vXrqa6uplWr9PopysvLKS8vTzzsndnosiOOnK0eaxGJQ2Nzds4U1mbWGvgqcFnS4mXA6KR1DOgMLE99fZ8+fZgxY0ZLhykiAkDr1q3p1q0b69evZ/369TU92I2VXEia2aqWiLElZTtnayiIiMSpsTk7l4aCXAxc6+5bk5Y9DoxJelwCbANUQYtI7Ip8OEhWc7aGgohIPohjHutEerSkZT8AXgc6mtkIM/uSmY1y9xnAB2Y2Plr1WOC37r4tu1GLiOyuGArrXMvZ6rEWkVyW1aEgZtYV+Gz08Dwzu4nwU+JPU1ad7e53RPc/A1xpZgsJ86leno1YRUQaUuiFdS7lbA0FEZF8kNXC2t0rgVuilvCzqNX1mkXAl1s4NBGRtBV6YZ1LOVtDQUQkH+TSGGsRkbxS6IV1LlKPtYjkMhXWIiJNpMI6ezQURETygQprEZEmUmGdPSqsRSQfqLAWEWkiFdbZozHWIpIPVFiLiDSRCuvsU4+1iOQyFdYiIk2kwjp7NBRERPKBCmsRkSZSYZ09GgoiIvlAhbWISBOpsM4+9ViLSC5TYS0i0kQqrLNHQ0FEJB+osBYRaSIV1tmjoSAikg9UWIuINFG3bt0AWL9+PdXV1TFHUxzUYy0iuUyFtYhIE7Vu3Zpu3brh7qxfvz7ucAqahoKISD5QYS0i0gwaDpIdGgoiIvlAhbWISDOosM4u9ViLSC5TYS0i0gzdu3cHVFi3NA0FEZF8oMJaRKQZ1GOdHRoKIiL5QIW1iEgzqLDOLvVYi0guU2EtItIMKqyzQ0NBRCQftMn2Ds3sGOBnwGfcfWG0rCNwDVAJ9AQudff1DT0nIhK3Qi+scyVnayiIiOSDrPZYm1k/oAtwUMpTNwMvuPsPgIeA2xr5nIhIrAq5sM7FnK0eaxHJZVktrN19JSHR1jCzgcA5wLRo0RPAyWZWUt9zWQlYRKQBhVxY51LO1lAQEckHWR8K4u7VtutvepOBVe6+KXp+h5ktAI4Attfz3MJsxi0iUptCLqwhd3K2CmsRyQdZL6xrMQhYnbKsEhgI7KjnuV1UVFRQWlpa87isrIyysrLMRioikqKphXV5eTnl5eWJh70zG1WLUs4WkaLT2JydC4W1A5tTlrUn9HzQwHM1+vTpw4wZMzIfnYhIPZpaWCcXkma2KuOBtZxYcrZ6rEUkTo3N2blQWC8jnDmerCuwnJDA63pORCR2hT4UpBax5GwV1iKSD3JhHuungYHRFE2YWVtgKPBsA8+JiMSuCAvrWHK2ptsTkXyQ9cLadp4FYwDuvgJ4GDgqWn40MNXdl9X3XBZDFhGpU7du3QBYv3491dXVMUeTebmWs9VjLSK5LKtDQcysK/DZ6OF5ZnaTu1cAFwDXmNkEoD+QfAZLfc+JiMSqTZs2dO3alcrKStavX1/Tg10IcilnayiIiOSDrBbW7l4J3BK15OVrgPPreE2dz4mI5IIePXpQWVnJ2rVrC6qwzqWc3bZtuN2+22mQIiK5IxfGWIuI5LUiHGedde3bh9utW+ONQ0SkPiqsRUSaSYV1y1NhLSL5QIW1iEgzJQrrdevWxRxJ4VJhLSL5QIW1iEgzqce65bVpA61aQXU17NgRdzQiIrVTYS0i0kwqrLNDvdYikutUWIuINJMK6+xQYS0iuU6FtYhIM6mwzg4V1iKS61RYi4g0kwrr7EgU1lu2xBuHiEhdVFiLiDRTorBes2ZNzJEUtq5dw+2GDfHGISJSFxXWIiLN1DWq+CorK2OOpLB17x5u9cOAiOQqFdYiIs3UuXNnADZt2hRzJIUtcbV4TRcuIrkqI4W1mR2fie2IiOSjTp06AbBx48aYI2mcfM3ZicJaPdYikqvapLOymR0LXAH0Z2dR3groB3TMbGgiIvkhV3usCy1nayiIiOS6tApr4C/AL4HXgepomQHnZDIoEZF8kiisc7DHuqBytnqsRSTXpVtYz3f361MXmtnsDMUjIpJ3cngoSEHl7D59wu2HH8Ybh4hIXdIdY32BmZ1Ry/KTMxGMiEg+Sh4K4u4xR7OLgsrZgweH26VL441DRKQu6RbWhwK3m1lVUqsGbmmB2ERE8kLbtm1p06YNVVVVbNu2Le5wkhVUzh40KNwuWxZvHCIidUl3KMhVQBkwA6iKlhnwueYGYmZ7ARcB84A9gXJ3/5+ZdQSuASqBnsCl7r6+ufsTEcmkzp07s27dOjZt2kT7xCUC41dQOVs91iKS69ItrB8FHnL3rckLzez3GYjlL8Ap7r7MzIYC/wb2Am4GHnP3u6Mz3G8Dzs7A/kREMiZRWG/cuJGePXvGHU5CQeXsAQOgVStYsQK2bYN27TKxVRGRzEl3KMg/gO+b2eFJ7Qjg+xmIZS8gumAtm4HuZjaQcPb6tGj5E8DJZlaSgf2JiGRM4gTGHJtyr6Bydtu2MGIEVFfDe+9lYosiIpmVbo/1Dwk/+Z0HJM7QMcIcqd9qZiz3EcYCHgd8FvgGMBlY5e6bANx9h5ktAI4AFjZzfyIiGdOmTUinVVVVDayZVQWXs/feG+bNg3fegfHjM7FFEZHMSbfH+tfAEHcvcffhUSsBTstALBcCW4D/AhvcfSowCFidsl4lMDAD+xMRyZhWrUI6zbHCuuBy9t57h9u3387UFkVEMifdHusl7l5Zy/L+GYilI/B3QmL+TdTL4YSfGJO1B7anvriiooLS0tKax2VlZZSVlWUgLBGRhrVu3RpofGFdXl5OeXl54mHvlomq8HL2vvuG25kzmxi1iEgTNDZnp1tYXwi8nLzAzD4N/Aa4I81tpforcJa7rzUzJyTsbxPOKk/WFVie+uI+ffowY8aMZoYgItI0icK6urq6gTWD5ELSzFa1UFgFl7MPPTTcvvhiGGvdKt3fXUVEmqCxOTvdlNTLzH4YbbSTmd0O3E74KbDJzKw3MMHdExeq/RmwHlgEDIymb8LM2gJDgWebsz8RkUzL0aEgBZezhw6FIUPCZc01HEREck1ahbW7nwA8bma/Bf4HjAYmALVd2SsdHwFbzGxQ0rLVwJvAw8BR0bKjganurssDiEhOSXcoSDYUas4+/PBw+5//ZHKrIiLNV+9QEDM7vI6nFgFzgWuBYcCpwHeaGoS7V5vZqcCVZjYT6Adc4u7rzewC4Bozm0AYF6iB0yIitSiWnH3KKfC3v8EDD8B3v5vprYuINF1DY6yvIcxVuhaobeDgHYSpmwbQjCQN4O5vAF+tZfka4PzmbFtEpKVt3RquwdKhQ4c4wyiKnH3ccdC+Pbz0EixfDgM1T5SI5IiGCuurgBnu/mF9K5lZc39WFBHJa5s3h8kwYi6siyJnd+kCJ54IU6fCHXfAFVfEHZGISFDvGGt3/1dDCTpa777MhSQikn+2bNkCQMeOHWOLoZhy9lejvvLycsihYe0iUuQ0UZGISAYkCuuYe6yLxlFHwahRsGQJ3Htv3NGIiAQqrEVEMmDDhg0AdOrUKeZIikOrVnDJJeH+T36iXmsRyQ0qrEVEmmnjxo1s3ryZDh060Llz57jDKRpf/CIMHw6zZ8Mf/xh3NCIiGSiszWw/M/urmX3XzEqaH5KISH758MMwrLlv376YWczR1K+QcnbbtvDzn4f73/sefNjg6HIRkZaViR7rrwFHu/t1wCGJK26JiBSLRGHdp0+fmCNplILK2WedBZ/8JKxZAxddBO5xRyQixSwThfVFwAHR/X8AbTOwTRGRvLF8+XIA+vXrF3MkjVJQOdsMbrkFOncOJzGWl8cdkYgUs7QLazM7zszOje4fBOyZuFytu1e7+/oMxygiktPmzp0LwJ577hlzJLsrhpw9YsTOgvpb34KZM+ONR0SKV1qFtZn9HLgX+BSAu78KfN7MzmmB2ERE8kKuFtbFlLPPOQfKymDr1nDxmEWL4o5IRIpRuj3WhwCDgeeTlj0KXJ2xiERE8sx7770HwOjRo2OOZDdFlbNvvBEmT4YVK8Jlzysq4o5IRIpNuoX14+6+Fkg+PSSvL40rItIcVVVVvP766wCMHz8+5mh2U1Q5u317eOABGD8e3n03FNkffBB3VCJSTNItrOeb2Q+AMWZ2lplNBb4O/CrzoYmI5L7Zs2dTWVnJ0KFDGTBgQNzhpCq6nN2jBzz+OIwbB++8A4cfDtFIHRGRFpdWYe3u/wCeAloDnwfWAp9099+1QGwiIjlv+vTpAEycODHmSHZXrDm7f3945hnYf3+YNw8mToQnn4w7KhEpBm3SfYG7TwemJx6bWWszO8rdlbZEpOg8/vjjAHz84x+POZLaFWvO7t07FNfnngvTpsGUKfCzn4XLoLduHXd0IlKo6i2szeyfQJcGttEbWAUUdJIWEUm1fft2HnvsMQBOOOGEmKNRzk7VtWsYc33FFXD11fD978O//w133glDhsQdnYgUooZ6rN+M1pkDVAPnA68DM5LWGQ70apHoRERy2PPPP8+6desYO3YsI0eOjDscUM7eTevW8MtfhrHWX/xi6MUeNy70Xl94IbRJ+3dbEZG6NZRSfg1YdFY5ZjbS3a9MXcnMHmqJ4EREctmf/vQnAE4//fR4A9lJObsOxx0Hb74JX/0qPPhguJDMn/4Ev/sdHHxw3NGJSKGo9+RFd1+XSNCRfma2y+g0M5sMHJSpgMxskpldbGanmllvM+toZr81s1+Y2S1m1i1T+xIRaaq1a9dy3333AfClL30p5mgC5ez69e0bhoY8+GAYCvL663DIIXDaaTB7dtzRiUghSHe6vfuAN8zsVjP7tZk9DDwO3JiJYMzsi8CJ7n6tuz/o7quAm4EX3P0HwEPAbZnYl4hIc9x+++1s3ryZI488khEjRsQdTl2Us2txyilhnuvLL4eOHUOxPW4cfPnLmppPRJon3en2Hgc+STjxZQxh6qYz3P2XzQ3EzI4AzgV+mLRsIHAOMC1a9ARwspmVNHd/IiJNtWnTJq699loALrnkkpijqZtydt06dw7jrOfNC8NDzOCOO2DsWPjMZ+B//4s7QhHJR+n2WOPuy939Cnc/2d0/C0wzswMyEMv1wLvAb83sMTM7FJgMrHL3TdG+dwALgCMysD8RkSa5+eabWblyJaWlpRx77LFxh1Mv5ez6DRwIt94aerC//OVwsuM998B++8Exx4RhIzt2xB2liOSLtM6HNrM7alncB6gCTm1qEGY2BtgP+KK7v2lmlwCPAT8DVqesXgkMTN1GRUUFpaWlNY/LysooKytrakgiIrVavnw5V111FQBXXXUVZtak7ZSXl1NeXp542Dsz0e1KObvx9twTbr8dfvxjuP56uO02eOKJ0AYPDr3aX/lKuPiMiBSfxuZsc/dGb9TMnifMfZr8okOB5939p02IM7Hd04Byd+8dPe4IVABXAae7+8SkdV8D7nL3a5O3UVpa6jNmJM8oJSKSWe7OWWedxT333MPJJ5/MQw9lZnINM5vp7qUNr5n2dpWzm2jNGvjzn+Hmm3eOu27dGo49Fj7/eTj5ZOjQId4YRSQe9eXsdGfw/Iq7z0nZ+FjgvCbGlhxHTSzuvtnM5hIuw9szZd2uwPJm7k9EJG1/+ctfuOeee+jUqRO/+c1v4g6nMZSzm6hnT/j2t+Gb34SnngoF9rRp8MgjoXXvHsZin3suHHqoruYoIkG6Jy/OqWXxVuArzYzjTaC7mfVLWraDkIwHRr0hmFlbYCjwbDP3JyKSltmzZ3PhhRcCcNNNNzF8+PCYI2qYcnbztWoFRx8N998Py5bBDTdAaSmsWwfl5XDEEWGoyIUXwpNPajy2SLFLq7A2swVmNj+pLSKcvPJAc4Jw99nAo8AZ0X56AO2AvwEPA0dFqx4NTHX3Zc3Zn4hIOioqKjjhhBPYuHEjZ511Fuedd17cITWKcnZm9e0berD/+194+2247DIYMQJWrIBbbgkFeP/+4STIhx6CDRvijlhEsi3dMdblwF3sHK9XDax09/eaHYhZb+AGwqV3hwK/d/d3zKwncA2wEOgPfN/dd0tXuT5eT0Ty04YNGzjmmGOYPn06BxxwAM8++yydO3fO6D5acIy1cnYLcw9T802dGlryhWbatoXDDoPjjw9XftxrrzCtn4jkt/pydrqF9eHu/lzKsjFAlbvPa16YzVMoSVpEckdlZSXHH388L7zwAkOGDOGVV15hwIABGd9PCxbWytlZ9s47YdjIv/4Fr7wC1dU7nxs6FKZMgSOPhMmToQX+K4lIFtSXs9Odx7q2CVtXAvemHZWISA5btWoVU6ZM4YUXXmDQoEE8+eSTLVJUtzDl7Czbe2+44gp46SX48EP4+9/hc5+DPn1g8WL4/e/hnHPC/Nljx8LXvgb/+Ad88EHckYtIJjRqVhAzu4gwP2lXM7u0llUez2hUIiIxmjNnDieccALvv/8+Q4YM4emnn2bkyJFxh9Voytm5YY894KyzQquuhtdeCyc4PvMMPP88zJkT2m3RRd9HjYKDDw5t0iTYZx9ok+7cXSISq0YPBTGzfYEvADemPLXN3WP/W7tQf1YUkex6+OGH+dznPsfatWvZb7/9mDZtGoMGDWrRfbbEUBDl7Ny2fXsotJ9+OhTaL7wAGzfuuk6nTnDggTuL7YkToV+/WjcnIlmUyTHWbd19ey3Lh7r74mbE2GzFnqRFpHm2bNnCpZdeyo03hjr01FNP5a9//WvGT1SsTQuOsVbOzhPbt8Nbb8HLL4c2fTq8//7u6w0eDAccsGtTsS2SXU2+QIyZXQz8y93fiRZ9z8xSp8E34EhgcnMDFRGJw/Tp0zn//POZNWsWbdq04eqrr+Y73/kOrVqlexpKvJSz81fbtrD//qF9/eth2YcfhhMgE8X2jBmwdGloyRf9HDQoFNilpSq2ReLW0Oit4UCPpMf7A52BFeycvskIUyqJiOSV9evXc/nll/O73/0Od2fkyJH8/e9/58ADD4w7tKZSzi4gffvCSSeFBmGc9nvvwcyZoc2YAa+/Hi5cs2wZ/POfu752n312bePGheElItJy0h0Ksj8wy923pSzfO6mHJBb6WVFEGmvHjh384Q9/4Ec/+hErV66kdevWXHzxxVx55ZV0iqHyaMGhIMrZBS612J45M4zdru3iNGYwcuTOQnvvvWHMGBg9WgW3SDqaPBSkFvu4+2spGx8DjARiTdIiIg1xd6ZNm8Zll13Gu+++C8CkSZO49dZbmTBhQszRtQjl7ALXqlWYtm/sWDj33LCsuhoWLQpjtpPbnDkwb15oD6Rce3Po0FBkjxkTtpW4HTRIF7URSUdjp9s7jPDz4TFmtiDl6V6Es86nZTg2EZGMqKqqYurUqfziF7/gjTfeAGDEiBFcffXVnHHGGViBVQ7K2cWtVSsYPjy0k0/euXzr1nBlyLfeglmzwv3Zs8NJkosXh/Z4ykSMnTuHHu2RI0MbNWrn/cGDw75EZKfG9livBe4kjN87NOW5bcDtmQxKRCQTNm7cyF133cV1113HnDlzABgwYACXXXYZX/va12jXrl3MEbYY5WzZTfv2MGFCaMm2b4cFC0KP9uzZO+fXnj0bVq0K47hff3337bVrByNG1F50l5SE/YkUm0YV1u7+lplNBj7t7krIIpLT5s2bx80338wf//hH1q5dC0BJSQmXXnop5513Hh06dIg5wpalnC3paNs29EqPHr3zRMmEjz4KY7jff39nmzcv3K5cubPXO5VZGF4ycmQovktKQhs2LNwOGACtU+erESkAjR5j7e7rqKOXw8w+7u4vZCwqEZE0bdy4kQceeIA///nPPPHEEzXLJ02axEUXXcSZZ55J27ZtY4wwu5SzJRN69QpXgZw0affnNmzYteBOLroXLw7jvBctgqee2v21bduGwjtRcCcX3SUl4ZLvKrwlHzU0j/U/gS4NbMOAcUDfTAUlItIY1dXVPPfcc/z5z3/mvvvuY0M0FUL79u05++yz+frXv05pacYn28hZytmSTV261D60BGDbNli4MBTZCxfu3j78cGcxXps2bXYW3sOGhfHcQ4bsvB0yBLp104mVknsa6rGeQxirt4ydc6CmMuDETAYlIlKXHTt28PzzzzN16lQeeOABli9fXvPcpEmT+MIXvsCZZ55Jr169YowyNsrZkhPatds5vKQ2mzaFXu3kYnvRop33V6yA+fNDq0uXLrsW28lFd+J+t24ZPzSRejVUWP8KqHT3rfWtZGY6u1xEWszmzZt5+umnuf/++3nooYdYtWpVzXPDhg3js5/9LJ///OcZXde3ePFQzpa80KnTzmkCa7N5867DSZYsCVecXLJkZ9uwAd59N7S6dO26s9AeODC0AQN2vT9gQPhDQCQT6i2s3X1Vfc8n+RQ6y1xEMsTdmTt3Lo8++iiPPfYYzzzzDFu2bKl5fs899+T000/n9NNP54ADDii46fKaSjlbCkXHjjvn1a6NO6xdu3vBnbifuK2shHfeCa0+vXvXXnQn7g8cCP37h7HhIvVJ9wIxuzCzocDXgYtQkhaRZqioqOC5557jqaee4tFHH2XBgl2nXz7ggAM48cQTOf300xk/fryK6SZQzpZCYQY9e4a27761r+MeZjVZujS0Dz6A5ctDS76/cmWYVnDVKnjzzfr326fPrj3dAwaEgjvR+vULt127avx3sUq7sDaz1sDJwPnAMcBKwpi+ZjOzdsB/gW+5+zNm1hG4BqgEegKXuvv6TOxLROKVKKSfeeYZnnnmGWbNmrXL87169WLKlCkce+yxTJkyhX79+sUUaX5ryZwdbV95W3KSGeyxR2j1XVi1qgoqKnYtuGsrwlesCOtVVEB0nak6dey4a8FdV+vXT/N9F5pGF9ZmNhz4CvBFYA9gEXAS8G/g4AzFcwlQkvT4ZuAxd7/bzI4FbgPOztC+RCRLqqurmTNnDtOnT2f69Om8+OKLvP3227us06FDBw455BCOOOIIpkyZQmlpKa0131aTZSlng/K25LnWrXcWuvVJLsATBffKlaHgTm4ffBBOzlywILSG9OzZuCJ8jz00BWE+aLCwNrNPE3o6PgHMA64H/gz8wN0fi1Z7qbmBmNkhwAfAmujxQOAcws+WAE8AU82sxN0XNnd/ItJy1qxZw6uvvsr06dN5+eWXeeWVV2ou1JLQoUMHDj30UI444ggmT57MQQcdRHt13TRbtnJ2tC/lbSkayQX4/vvXv+6GDbsX3LW1lSthzZrQ6jsJM7H/vn0bV4RrKEp8GtNjvZlwCdyXgK+4+3sAZlbXVE5pM7POwBnu/n9mdmW0eDKwyt03Abj7DjNbABwBLMzUvkWkeVavXs3rr7/Oa6+9VtPmzp2723oDBw7k4IMP5uCDD2bSpEmUlpaqkG4ZLZ6zo+0pb4vUoUuXcIn3UaPqX6+6OowDb0wRvnp16A3/4IOG99+x487x3g0NRSnwC9FmXYOFtbs/DDwcnfRyvpmVAA8BNefGmpm5e3OS9qXA1SnLBgGrU5ZVAgObsR8RaSJ3Z/ny5bzxxhu7FNGLFi3abd127dpxwAEHMGnSpJpCesiQITFEXXyylLNBeVuk2Vq1CjOS9O4N48fXv+62beHCOg0V4ImhKIk5wRtS11CU1BMz99gjxCv1S+eS5ouBH5pZG+BUYLCZ/RT4A/BJoLwpAZjZFGCGu3+YuktCz0uy9sD22rZTUVGxyxXWysrKKCsra0pIIkVvzZo1zJo1i1mzZvHWW2/V3F+zZs1u63bs2JGPfexj7L///uy///7st99+jBs3jnaaGLZe5eXllJfXpM3emd5+S+VsyEzeVs4WSU+7dmE+7sGDG163oaEoifHh6Q5FSfSCpxbdqQV5586ZOeZc0ticbc3ptLMuW30AACAASURBVDCzUcA3gPPdvVMTt/E4sF/Sop7ABsK4wHPdfXTSunOBH7n7XanbKS0t9RkzZjQlBJGitWbNGubMmcPs2bN5++23a4roZcuW1bp+r1692GeffWqK6P33358xY8boJMNmMrOZ7t7i117PRM6OttPsvK2cLRK/5KEoH3xQfzH+0UeN326XLg33gg8YEMaM5+PXR305u1nzWLv7POBbZjazGZs5h9CjkfAy8H/AK8AlZtbR3TebWVtgKPBsM/YlUnSqqqpYuHAhs2fPrimiE7cffpja4Rh07NiRcePGsc8++zB+/Pia2/79+2v+6DyWoZwNytsiBSGdoShbt9Y+C0pyL3ji/oYNMG9eaPVp0ybMC564Ombq5ekHDw695PlUfDersE5w9zub8dqK5MdmVgVUuPtiM3sYOAp4GDgamOrutXeliRQxd2fFihW8//77zJs3j7lz59YUz3PnzmXbtm21vq5Tp06MHj2aMWPGMH78+Jo2fPhw9UIXsObk7Oj1ytsiRaZ9exg6NLT6uMP69Q0X38uXh+kLFy8OrS6pxffw4TByZGgjRoRlufR1lZHCugVdAFxjZhOA/oAG4EnRqqqqYsmSJcybN6+mgE6+3bRpU52vHTRoEGPHjmXs2LGMGTOm5nbw4MG00tkoklnK2yJFzAy6dw+trkvSJ2zZEgrsui5Jv3Rpw8V327ZQUrKz2N5zz9D7Pn586O3OtpwrrN29JOn+GsJ8rCJFYd26dSxcuJCFCxeyYMEC3n///ZrieeHChWzfXuu5u0AY/zxy5EhGjRrFqFGjaorn0aNH07Vr1ywehRQb5W0RaYoOHUKv84gRda+zZQssWxYK7cWLw0V33n8f5s8PtytWwNy5oaXq02dnkT1xIhx+eOj5bkk5V1iLFLL169fXFM61tdpm3Ug2cODAmuI5+XbkyJH07NkzS0chIiKSHR067OyNrs3GjaHYThTac+bAW2/BrFmht/vpp0P77W/D+sOGwXHHwbnnwqGHZv5COiqsRTLE3Vm3bh1Lliyps3D+qIHTqjt27EhJSUlNSxTNo0aNYsSIEXTq1OSJHERERApO5847e6WTuYde7lmz4H//gxdfDG3RIrj11tCOOw7uuSfMYpIpKqxFGmnLli0sXbqUxYsXs2TJklpvN2zYUO82OnToQElJCcOHD9+lgE60Pn36aNYNERGRZjLbebLl8ceHZVVV8MYbcN11cNdd8OijcO+98MUvZm6/KqxFCCcGrlixYrdiOfl+XVPTJevcuTNDhw5l6NChtRbPffv2VeEsIiKSJdu3h/HXb70F06fD88/DzGjC0bZt4WMfy+z+VFhLwXN3Vq9eXVMoJ7dE0bxs2TJ27NhR73batGnDoEGDGDp0KEOGDKn1tkePHiqcRUREsmjLlnD59gULdt4uWBDGW8+eHYrrZB06wEknwQ9+oMJaZDeJcc11taVLl7J5c+pVlnfXt2/fOgvmIUOG0L9/f83tLCIikiVbt8KHH4b5r5cvr/u2oR+Uhw8PY7D33x8OOwwmTWq5y66rsJactnHjxt2K5NTCubKyssHtdOvWjSFDhuzSkgvnwYMH06FDhywckYiISPFJXDimogJWrQq3yfdrW9aIr3cgXCBm6NBQQCe3UaNg770ze3JiQ1RYS2y2bt1aa6Gc3Bqafg7C1QNTi+YhQ4YwePDgmvvdunXLwhGJiIgULvdQ7K5ZE9pHH+28X1/76CNYuxaqq9PbX5s2YS7qAQPC1Rfruu3bN6ybC3IkDCk0O3bsYPny5fUWzY05GbBdu3a7FMi1tZ49e2pcs4iISD3cYdMmWLduZ1u/ftfHdS1bty4UyGvXhpk1mqpLl1Ao9+kDvXvXfpt8v3v3zM8z3dJUWEvaqqurWbFiRb29zR988AHVDfxp2rp1awYOHFhv0dynTx9dcltERIpaVVUoeBOtKcXx+vXNK4oTOneGnj2hV69wW1ur7bkePcIsHIVOhbXswt1ZtWpVvT3NjZlBw8zo379/vUXzgAEDdDKgiIgUrOrqMHQiuSBuyv2NGzMTT4cOoRc40bp12/VxfcsTxXG7dpmJpVCpsC4yVVVVLF++nPfff5/58+ezYMGCXeZsXrp0KVu2bGlwO7179653iMagQYNop0+fiIjkIfdQzCb3/jalIG7syXcNMYOuXUPBm2i1FcR1FcWJZfpabnkqrAvQhg0bWLBgAfPnz68poBP3Fy5cyLZt2+p9fffu3evtaR48eDAdO3bM0tGIiIg0jjts3tz0Qjhxv7Iy/RPt6tKly67FcFPud+kCGhWZH1RY56mtW7cyb9485syZw+zZs5k9ezbz5s1j/vz5rFy5st7X9uvXjxEjRtS0xLRzida1a9csHYWIiEgoiLdubX5BnKlxxACdOjWtGE5+3LVrmApOiocK6xy3evVq3n333ZriOVFIz58/v86TA9u1a8fw4cNrCueRI0fW3B8+fDhdsjmho4iIFLTt23c9Qa6ugreh51OvjtdUHTo0r3c40XJl+jbJL/pvkyN27NjBe++9xxtvvMGbb75Zc7ts2bJa12/VqhUjR45k7NixjBkzhrFjx7LnnnsycuRIBg4cqJMCRUQkLdXVochdvTq0jz6q/X7q40yNI27bdtcxwU0tiDWOWOKkwjoGVVVVvPvuu0yfPp1XXnmF1157jbfffputW7futm6nTp3Ye++9GTt27C5t1KhRtG/fPoboRUQkH1RVhQJ45cpwyefEbfL9Vat2Fshr1jRtXHGrVmG2iOYUxN27g77SpBCosM6CdevW8fzzz/Pyyy8zffp0/vvf/9Z6Ge5hw4YxYcIEJkyYwL777suECRMYOXKk5nEWEZEamzbBsmWwdGm4XbYsFMqpBXRFRfqFcrdusMceofXqVfv91MfduunEOpGEnCmszew04DqgG/BX4LvuvsPMegM/B9YCbYFL3T1DI7FaxpYtW3jppZd48sknefLJJ5kxYwZVKWdTDB06lEmTJjFx4kQOPPBA9t13X7p37x5TxCIi6SmknJ1LtmyBhQthwYKdhXPq7Zo1jd9er17hcs/9+tV+26fPziK5Z8/iuICHSEvKicLazIYCpwJnAOOAW4AlwLXAvcBl7v6KmX0NuAr4flyx1mXt2rU88sgj3H///Tz22GNs2rSp5rk2bdpwyCGHcPjhhzNx4kQmTpzIgAEDYoxWRKTpCiFnx2nVKnjvPZg/f/dWx2k1u2jXDgYOhMGDYdCg0Pr3r71o1nhjkezKicIaGAZ8yd13ADPNbDxwpJm9AIxx91ei9aYBc8zsZ+6eoesQNd327dt55JFHuP322/nPf/7D9qRTmidMmMBRRx3FUUcdxWGHHaYp7ESkkORlzs62NWvgrbfg7bd3bRUVdb+mTRsYNgyGD4chQ0LRnFxADx4MvXuHC4aISO7JicLa3Z9PWbQM6Ap8AliQtN4yC9mkFHg2awGmWLlyJTfccAN33HFHzZzRrVq1YvLkyZx22mmceuqpDBkyJK7wRERaVL7l7GzYvh3+9z945RV49dVw+957ta/bpQuMGQMjR8KIETtvR4wIhbOmeRPJX7n68T0Q+DVQBqxOea4SGJj6goqKCkpLS2sel5WVUVZWltGgVqxYwS9/+UvKy8trLvu91157cf7553PuuefSt2/fjO5PRApXeXk55eXliYe944wlA3IyZ7ck91A4P/44/Oc/8PTTsGHDruu0bw/77APjxu3ahgzRyX4i+aaxOTvnCmszGwVUuPsbZubA5pRV2gO7nQjTp08fZsyY0SIx7dixg5tuuokrr7yyZjaPU045hYsvvphDDz0U029yIpKm5ELSzFbFHE6T5WLObknz5sHf/x7au+/u+tzo0TBpEkycGNq+++pkQJFC0dicnVOFtZm1Br4KXBYtWgaMTnregM7A8mzFtGTJEj7zmc/w8ssvA3DiiSfy85//nH333TdbIYiI5KRczNktoboaHnsMrrsOnnpq5/JevWDKFDjmmNAGD44vRhHJDTlVWAMXA9e6e+JKKY8DX0t6vgTYBmSlm+PVV1/lhBNOYNWqVQwaNIhbbrmFk046KRu7FhHJBzmVs1vC9Olw0UUwc2Z43LkznHYanH02HH20eqRFZFc5U1ib2Q8IybejmY0AJgPPAR+Y2Xh3nwUcC/zW3be1dDwvv/wyU6ZMobKykk9+8pP87W9/o3fvfB8GKSKSGbmWszOtqgouvxx+9avweMAA+Pa3oawsXGVQRKQ2OVFYm9kVwE9TFs929zvM7DPAlWa2EOgCXN7S8SxevJhTTjmFyspKzjrrLO68807aqltCRATIvZydaTt2wFlnwdSp0Lo1fO97ocju3DnuyEQk1+VEYe3uPwN+Vsdzi4AvZyuW6upqzj77bCoqKjjmmGP4y1/+QhvNfSQiUiOXcnZLuOyyUFT36AEPPACTJ8cdkYjkC1WMKW677TZeeuklBgwYwD/+8Q8V1SIiReS11+D660NP9T//CYcdFndEIpJPNJNmksrKSn74wx8CcOONN9KrV6+YIxIRkWy69towR/U3vqGiWkTSp8I6yU033cTq1as55JBDOP300+MOR0REsmjt2jAEpFUr+M534o5GRPKRCuvI1q1buf766wH4yU9+oou+iIgUmX/+E7ZtgyOOgKFD445GRPKRCuvIfffdx6pVq/jYxz7GUUcdFXc4IiKSZXffHW7PPDPeOEQkf6mwjtxyyy0AXHDBBeqtFhEpMmvWwH/+E4aBaCSgiDSVCmtg/vz5vPjii3Tp0oVzzjkn7nBERCTLHnwwzF/9iU9Anz5xRyMi+UqFNTB16lQATjrpJLp06RJzNCIikm0PPRRuzzgj3jhEJL+psGZnYa2ZQEREis+2bfDkk+H+8cfHG4uI5LeiL6w/+ugjXn31Vdq1a8exxx4bdzgiIpJl06fDhg2w994wZEjc0YhIPiv6wvq5557D3Zk0aRKdO3eOOxwREcmy118Ptx//eLxxiEj+K/rC+umnnwbgyCOPjDkSERGJw5tvhtt99403DhHJf0VfWD/33HOACmsRkWI1e3a43XvveOMQkfxX1IX1li1bmDVrFq1ateLAAw+MOxwREYnBkiXhdtiweOMQkfxX1IX1m2++yY4dO9hrr73o1KlT3OGIiEiWVVXB8uXh/qBB8cYiIvmvqAvrmTNnAnDAAQfEHImIiMRhxYpQXPftC+3bxx2NiOQ7FdaosBYRKVarV4dbXW1RRDKhTdwBNIaZdQSuASqBnsCl7r6+udt9MzoVfL/99mvupkREJNJSObslrFkTbnv2jDcOESkM+dJjfTPwgrv/AHgIuK25G3R33nnnHQD21qngIiKZlPGc3VLWrg23KqxFJBNyvrA2s4HAOcC0aNETwMlmVtKc7S5evJiNGzfSr18/9thjj+YFKSIiQMvl7JaS6LHu0SPeOESkMOR8YQ1MBla5+yYAd98BLACOaM5G1VstItIiJtMCObulqMdaRDIpHwrrQcDqlGWVwMDmbFSFtYhIi2iRnN1S1GMtIpmUDycvOrA5ZVl7YHvygoqKCkpLS2sel5WVUVZWVudGhw0bxkknncQhhxySwVBFRBqnvLyc8vLyxMPeccaSYS2Ss1vKnnvCSSfB+PFZ37WI5JHG5mxz9+xE1ERmdjbwE3cfnbRsLvAjd78rsay0tNRnzJgRR4giIs1iZjPdvbThNXOfcraIFLr6cnY+DAV5GhgYTd+EmbUFhgLPxhqViIjURjlbRIpWzhfW7r4CeBg4Klp0NDDV3ZfFF5WIiNRGOVtEilk+jLEGuAC4xswmAP2B7A/EExGRxlLOFpGilBeFtbuvAc6POw4REWmYcraIFKucHwqSaUlndBYEHU9uK6TjKaRjgcI7nmJRaO9boR0P6JjyQaEdD+TOMamwznM6ntxWSMdTSMcChXc8xaLQ3rdCOx7QMeWDQjseyJ1jKrrCWkRERESkJeT8PNaNZWYVwKJGrNobWNXC4WSTjie3FdLxFNKxQG4dzzB37xN3ENmURs5OlUvvWyYU2vGAjikfFNrxQHaPqc6cXTCFtYiIiIhInDQUREREREQkA1RYi4iIiIhkQNEU1mbW0cx+a2a/MLNbzKxb3DE1l5kdY2avmFlJ3LE0h5mdZmYLzGy1md1gZnkxv3p9zOwIM5tjZmvM7Ddxx5MJZtbOzN4ws8lxx5IJZvYzM/OovRF3PLJTY/O1mR1nZjeZ2c1mdka242ysxuY4M9vTzKqS/l/ul+1Y09GYz1C+vEcA0fepp7Tja1mvh5mtT1rnU3HEW5faaoN0aqBcfM/qOKZG1w7Z/GzlfQGThpuBx9z9bjM7FrgNODvmmJrMzPoBXYCD4o6lOcxsKHAqcAYwDrgFWAJcG2dczWFm3YFDgUnAIcCDZvawuz8Rb2TNdglQEncQmWBmnQknuhwTLWrKSXTSchrM12Y2Gvg5UOru1Wb2nJnNdvdZMcRbpzRz3FeA44BqYLu7v56tONPVmM9QvrxHAGa2F/AycDGwPVp8O/BULaufD3wBqIweP93iATZSPbVBo2qgXHzPajumJtQOWftsFcXJi2Y2EFgA9HT3TdFfNeuAce6+MNbgmsHMWgFVwPB8PQ4zOwx42d13RI+vIbwvJ8QbWdOZWUdgi0cfLjObCVzs7jmTfNNlZocAY4ErgfPc/Zl4I2oeM/t2dPdWd98SazCyi8bmazO7FfjI3X8QPb4cGOXuX4wh7Do1NsdFxcMfgK+5+9LsR5qexnyG8uU9glCoufvipMfDgGvc/TMp67UHHgYudPe5WQ6zUVJrg3RqoFx9z2o5pkbXDtn+bBXLUJDJwCp33wQQvRELgCPiDKq53L067hiay92fT3wwIsuAxXWtnw/cfXNSUd0ZeBd4JtagmiE6hjPc/Y64Y8kEMzPgLOA6YKWZ5e0vVwVqMo3L15+IlifMiV6bU9LIcWcDRwJLzOyP0ecuJ6XxGcqL9wgguaiOfAp4sJZVTwQmAu+Z2SNm1rvFg0tTLbXBZBpfA+Xke5Z6TGnWDln9bBVLYT0IWJ2yrBIYGEMsUr8DgVvjDiITor+oHwXaAp1iDqc5LgWujjuITPFgEuFn7N8CfzOz42IOS3ZqbL5OXS9fcnqtOc7dfwN0A04Gjib0sOWkND5D+foeARwPPJK60N2nAt0JRWkJ8M+oNzWXpVMD5et7VmftkO3PVq7/Z8gUBzanLGvPznFUkgPMbBRQ4e6FciLZPOBPwFHAr+MNpWnMbAoww90/jDuWTHP3Ne5+BfBL4NsNrS9Z09h8nbpezuf0hnKcu1e5+zRgCnCamQ3IaoBpasRnKO/eI4CoF3qLu6+v7fnoD4vnCL27Ywk92LksnRoo796zxtQO2fxsFUthvQzombKsK7A8hlikFmbWGvgqcFncsWSKu38QDZ+4mPwddnQxcIeZrTKzVcAQ4CEz+17McWXSjcDQuIOQGo3N16nr5XROTyfHufs7wJPkz//Luj5DefUeJTkZ+GdDK7n7SuBucv99SqcGyqv3LN3aIRufrWIprJ8GBkYnlWFmbQn/qM/GGpUkuxi41t23xh1IC5hBSFb56BzgY0ltOeHs6oIYrhOpBl6LOwip0dh8/TgwJunxaHJodoZapJvjNgKzWzCeTKrrM5Rv71FCowrrSBXwvxaMJRPSqYHy7T1rSu3Qop+toiis3X0F4Szeo6JFRwNT3T1fix2g5gQSAKt3xRxnZj8AXgc6mtkIM/tS9NNOXjKzTilzZB4P3BBXPM3h7hXuvjTRCF8iFXX9RJoPzGygmZ1jZq2iz9B3gcvjjkuC+vK1mV1sZokv/VsJP+smcuEnCD2nOaeOHLenmV1lZv2jdc4zs77R/cOAZ919XYxh16m+z1C+vkcJ0YltHZKHvyUfk5mdYWYjovtjCTNozIkn2tql1gYN1UD58J7VVu/UVTtYEN9ny92LohF+2vg94cP/W6BL3DE183i6AhcQxkP9BOgTd0xNPI4romNIbu/GHVczj+kgoAJ4njC2+tNxx5TBY1sITI47jmYewzjC2eOzCF8iY+OOSW2396jWfA3MJMxQk1jvs8D1UTsl7rjrOJZacxzQkTD7wqRovYeAVcDfgc/HHXcDx1TnZygf36OUYzsD+HbKsppjAm4C1gL3A98EWsUdc0qstdYG9dVAuf6e1XZM9dUOcX+2imIeaxERERGRllYUQ0FERERERFqaCmsRERERkQxQYS0iIiIikgEqrEVEREREMkCFtYiIiIhIBqiwFhERERHJABXWInUws6PM7CUzm9zAeuPN7B4z+2EG9/3rTG1LRCQOZravmd2bydxYx34eNrOzWnIfLcXMjjazlxv6npH80SbuAEQaYmZ7AT8HPgIOBPYFPufuf23hXS8BJjRivcWEyfdbZ2Kn0eVmt2diWyIiAGZ2IPBL4HDgP4QLaPQH2gK3uPu/W2C384DuZCg31uNmwsVqMi7Kx6cRLoTzHvBy9FR74HTgIHdf2IxdNPZ7RvKECmvJBw8CZ7r7G2bWCvhjNnbq7u+Z2epGrLfezD7I4K6PAJ7L4PZEpMi5+3/N7B+EqySemFhuZicAD5jZxe5+c4b3uSnDuREAMzsJeCtR0Lr7vzK9jwR33w7cbWbfAJ5x9yuS4nioKduMLq99pLvf7e5zzGxVhsKVHKChIJLTzKwXMBrYCODu1cCVhMuXZkN1htdrjCOBpzO4PRERgB2pC9z9EeDHwP8zs/4tsM+M5mozG0zooc623f7tgCeAZelsxMzaA3cSLrudkMnvD4mZCmvJae7+ETAHeMzMJkXLFgGvmlkPM/upmb1qZieb2Qoze9fM9k+83szONrPLzewZM/tF0vJSM/uxmd1tZvebWfdoeQcz+52ZXWlmfyD8jJkWM/teFNOZZnagmT0Ybe+nUYx3m1mfaPk6M/tRyibaufvWaFvfNbP/M7MXzezqtP8BRUQa9gfCkJBPA5jZSDO7ysxuN7Mno2IWM2sf5bLLzex5M/tMYgNmdqGZ/dDMppvZRXXtqLacbGZnRDnzPDOba2ZXm9lAM/tjlE+fN7NPRJs4BxgMXGpmp0bnwjxhZl9I2sf5ZvZzM7vTzP5sZp3NrLuZ/cjMZpjZJ8xsjpnNN7OSpvyDRf8mk6MebczsNDO7wcyuM7PHzKxLdAw3mNm06Pvmg+jfeG/g02b2f0mb7Bl9F1Wa2fl17LOVmf3KzMqifbiZPZB4fyRHuLuaWk43YBywiPBX/R1A32h5a0KS+ogw1q0H8DjwNuGPxoOBK6J1exJ6vU8AugH/SNr+i8Cvo/v/DyiL7vcDqgjJs6EY/0To9elLGIvXKVreHvg38BAwkPCFsB34EdAJOCB6nFg/MX48cdw3R/e7ApfF/V6oqanlbwPOA5bW8dxK4MYor94PtI6W3wXcG92/Hjgtuv8loCK6fzbw2ej+6ChX7xM9/hPw4+h+XTm5F7AW+A1wEHA0cC1wZbTuJcC0pFgdKInudwTeBc6LHk8B7kpa92HgNsLQ109F3xenRcf5GPDzRv7bPRN9V1wN/Jow3vropOdXAUOj+28TvpPaA5cDS4GJwJejvP9MIt5o/YXAddFzZwOL64jhK8DD0f1+0b/zJ+P+f6W2a9MYa8l57v62me0D/Ay4EDjJzI72MOa6Aljv7lMBzOzHwAvACMKXSLWZfTva1L8JCfwEYI+k5e9Gr+0CXBC9FndfaWZL0gi1N/Bb4MvuvinaxlYzWwEscPfl0X5WAM96GH/4OiHh9yacBHkC4UsAYBPweTN7D7iF8EeFiEhLcEKxeRAwDPiGmQGsBrpbePBVQqcAwJ8JQyEg5NoZKbl2APBWyj7Oo5ac7O4fmdk64CF3fxXAzBYDG81sAKGHt0utQbtvNrMPkxZdFG034Xbgnmj5WsL3xf3RPmZEcTbW0x6NsTazuwmdOQmfdvfFZnYooaDuEuX/5cD77v4K8Er02tq2PS36TvhvPTHtD2yOjnulmX3ErkNKJAeosJacZmatgSEeTlL5ZjQ8419AOaEHINXb0W0PYCjwF3e/K1r2m2iblxIK3d+k7OsAQkKsbGK4o4GPAyOBN5KWp44xrKp5wr06SrKJYVl7eBj+grsvMLPzgJuA7wDnAh8iIpJBZtaD8Gvb24S8ubqW/NiP0KPaDsDdqwidAUSv+Ym7vxQ93uW1SWrNyRFn11y5FLgUmAtMB4bXcwjJr9sTeCrp8fuEYS592T0X76DpM5a8RvgDJGGtmf0KuJcQe6J6Tj2uhuyg7trsUeAWM+tAeC82E3q/JYdojLXkuppxfwDu/gbwA2BsHeu3IySxecAHwBnJT5rZxGj5CWbWMWX5+ujhXk2M9SXgBuBvUeJLi5kNJOlEGDMbRJgRZTThZMa76nipiEhzfIHwC9m9hPx4aNRTDNTkx1XANuDYpOX7mFlXUnKtmbWJOipS1ZWTa3MrsNDDtKpb0ziWxYQe7ppdAOsIQ10yxt2ro86PcdEfJv8CrnL3GZncT8o+pwHTCL+sng0c4e7rWmp/0jQqrCUffCMqOhMGEcbGJfQys3bR/cmE8XVrgb8DnzKzm8xsctSbsJGQADsDD5vZlOhEkbGEYvx14Kdm1s7MehPGNvc1s9ZmdpmZfbaOGC1qVxK+fK6v5bnU9VMfn0gYD5iwJ3Cqu68HvlHLa0RE0rFbT6iZHQ78BPiqu1cQeodXAv8ys5OiExQPi3qo/0GYPeRUMzsK+IK7VxJy7TfN7CdmdgShgyExjC45/9WVkxOSe4/3A/qYWSfgUKCjmY2InttOONkv0cGSvI9bgdPNrGf0eCLwe3ffQah5as3FZjYpiqt9Hf92rUipmaI/Pk4g/ErZkzDEcAwwBOiUdGJkaq/4tlri3zWoWsaLRO9FN+BJQq98On9wSJZoKIjkgyHA22b2COFnMgO+lvT8DuByM1tL6N39BoC7P25m3yT8nHgq4eS/WQBmgxXRlQAAIABJREFUdgrwO8IXxR3AJe7uZnYmYSqkt4G7CT/pfYxwQYUzCUNMdrkwjZkdRDgpZyjhp8EHgZ9E8fwbmASUmNlIQk9Kf+BMM5vNzt6b84DO7l6ecux/tHBhhx2Ek4VERNIW5ZFzCMXf74EKwglw3YHj3P1lAHffFuXH3xNy3UOEsdUA3yJ0StxJ+IUuMRPHHwi57wLCkLUL3P3DqHCcBIw0s7/WlZPN7AzCuOKvmNlsd19GOJHyWkJRfSdwCrAPMB/4C6F3/fMWZosaBxxvZv9y9/vNbE9gqpm9SPjV8/KoZ/3TQP/o+N4mXCxngJmNi+I8izBu/MGkf7e2hNw/ARgUFd5OGNt8WrTN/wHPAv8lDG95lJDTHyacMDnBzE5198R27yZcrGdDVIj3B84wswWEnmii16des2Ex4eTMU4EOQCsze4JwAmO2pqCVBpjeC8lnFi4D+yd3L8nS/i5xd11uXESkwETDUtq7e05eoMvMLiHM0LIwetyB8GvDrxLn5kj8NBRE8l1twyxaZkfhJ9PXsrEvERHJnujkzP1zuKjuAFwBlCQt7g1sVFGdWzQURPJWNL7ts4Sf9j7t7ve28C7fj36iFBGRwtLG3W+JO4i6uPsWC5dVvzE6sX05YdjKFfW/UrJNQ0FERERERDJAQ0FERERERDJAhbWIiIiISAaosBYRERERyYCCOXmxd+/eXlJSkvbrKioq6NOnT+YDyjE6zsKi4ywsM2fO3ODuXeOOI5uUsxtHx1vYdLz5aebMmavcvdYDKZjCuqSkhBkz0r+SaGlpaZNel290nIVFx1lYzKzofj1Uzm4cHW9h0/HmJzNbVNdzRZfMRURERERaggprEREREZEMKPrCuqzs/7N33/FxlVf+xz9HktWLZUvu3cbYGEwzpiQBA6YlhCSUDSFlSXMSNpuyYTe/kCybBdKATdmQgjaw6SQbCCEQkgAJoYRqU4wBY4x7l2TZapbVzu+POyONZFXPHc1o5vt+ve5rZu7cuffcGeno6Jnnee6KZIcwInSe6UXnmXaqkx3AaJFBPxOAzjfd6XzTT9pcIGbJkiWeDv12RCTzmNkqd1+S7DhGknK2iIxWA+XsjG+xFhEREREJgwprEREZFdob25MdgojIgDK7sL7tNjjhBLjllmRHIiIiA2jd3crKxSvZ+l9bkx2KiEi/Mruwrq2F55+Hzf1ORygiIilg36P7aNnYwoYvbqCjuSPZ4YiI9CmzC+vi4uC2sTG5cYiIyIAmXDaBoqOL8INO08tNyQ5HRKRPKqxBhbWIyChQdHQRAM2vNCc5EhGRvqmwBmhS64eISKrLm5kHwMFtB5MciYhI30a8sDazc8zsaTObFbPuYjPbaGa1ZvYdM8uJea7CzG41s2+Y2TfNbExowajFWkRkQKmUs/OmqbAWkdQ2ooW1mU0EioGlMetmAO8ELgU+C3wE+EzMy34D3O7unwfWAdeFFpAKaxGRfqVazo4W1i1bW8LapYhIqEa0sHb33cA9vVbPBD7k7qvc/afA94AzAczsFOBId386su29wD+bWVEoAamwFhHpV6rl7K4W661qsRaR1DTiXUHcvbPX48fcPXbW/+3Alsj9s4CNMdtuj9wN59K/KqxFRAaUSjk7b0pQWLfubA1jdyIioUvFwYsnAT+M3J8K1PZ6vgGYEsqRVFiLiMRrxHL2mAljIAvaqtvobO0c/AUiIiMspQprM5sHVLv7i5FVDhzotVke0Nb7tdXV1SxZsqRrqaqqGvyA0cK6vh7c44hcRGR4qqqquvIVUJHseA7HSOfsrJwscifmAtC6S63WIjJyhpqzc/p7YqSZWTbwMeD/xazeDsyP2caAImBH79dXVlaycuXK4R20oADy86GlBZqboSicrtsiIoNZsWIFK1asAMDMapIczrAlJWcDuVNyad3ZSuvOVvJn5A/79SIih2OoOTuVWqyvBm5299hRKQ8CR8Y8ngW0AsPPxn0xg4rIPx01o+7vmohIMo18zqa7n/XBHRrAKCKpJxnzWFv0bsy6a4DngQIzm2NmHzKzee6+EthpZkdHNj0f+K67h/cdYLSwru3dLVBERFItZ+dOjnQF2aGuICKSeka0K4iZlQDvizy80sxuIfgq8fpem65199sj998NXGtmmwjmU/1iqEGNHx/cqsVaRKSHVMzZarEWkVQ2ooW1uzcAP4gsUTdElv5esxn4cMKCUlcQEZE+pWLOzp0SabHWlHsikoJSqY91cqiwFhEZNdRiLSKpTIV1ZWVwu3t3cuMQEZFBdbVYq4+1iKQgFdbTpgW3W7cmNw4RERlUdPCiWqxFJBWpsJ45M7jdsmXg7UREJOlyK3MhG9pr2+k8qKsvikhqUWE9Y0Zwu3lzcuMQEZFBWbaRO0lXXxSR1KTCOlpYb9sGHR3JjUVERAalAYwikqpUWOfnw6RJ0N6uftYiIqOABjCKSKpSYQ2wcGFw+/LLyY1DREQGpRZrEUlVKqwBjjkmuF2zJrlxiIjIoLoua66LxIhIilFhDXD00cGtCmsRkZSnFmsRSVUqrKG7xXr16uTGISIig1IfaxFJVSqsARYvhpycoMW6vj7Z0YiIyADUYi0iqUqFNUBhIZx4InR2whNPJDsaEREZgFqsRSRVqbCOestbgtvHHktuHCIiMqAxFWPIKsqiva6dttq2ZIcjItJFhXXUsmXB7QMPJDUMEREZmJlRuKAQgObXmpMcjYhINxXWUWeeCQUFsHIlbN+e7GhERGQAhUdGCuu1KqxFJHWosI4qLIRzzw3u//73yY1FREQGpBZrEUlFKqxjveMdwe2ddyY3DhERGVBXYa0WaxFJISNeWJvZOWb2tJnNillXYGbfNbOvmtkPzKx0KM+F7l3vgvx8+OtfYePGhB1GRGS0SNWcXbgwKKyb1jQlYvciIodlRAtrM5sIFANLez31feBxd78GuAe4dYjPhWvsWLjkkuD+//5vwg4jIjIapHLOLlxQSFZ+Fi0bWmir08wgIpIaRrSwdvfdBIm2i5lNAa4A7o2segi4yMxmDfRcwoL8yEeC26oqOHAgYYcREUl1qZyzs3KyKD6uGIDG5xrD3r2IyGEZ8a4g7t7Za9UyoMbdmyPPtwMbgTMGeS4xzjgDjj8edu+GH/84YYcRERkNUjlnlywpAaBhZUMidi8iMmypMHhxKlDba10DMGWQ53qorq5myZIlXUtVVdXhRWMGX/hCcP/GG6FNXzGKSPiqqqq68hVQkex4hiFlcnbxiUGLdf2z9cN+rYjIcAw1Z+eMXEj9cqB3n4s8IFrRDvRcl8rKSlauXBlORBdfDPPnw7p1cNtt8PGPh7NfEZGIFStWsGLFCgDMrCbJ4QxHyuTs0lOCcZH7H9+Pu2Nmce1PRKQ/Q83ZqdBivR0o77WuBNgxyHOJk50NN9wQ3L/2WqhXa4iISETK5OzCIwvJnZRL2+42TbsnIikhFQrrh4EpZlYAYGZjgBnAI4M8l1iXXgqnnQbV1fDVryb8cCIio0TK5GwzY+yZYwHY99d9iTiEiMiwJGMe6+h3dQbg7ruA+4CzI+uXA3e5+/aBnhuBQOGb3wxub74Znnsu4YcUEUk1qZ6zo4V13V/rEnUIEZEhG+l5rEuAaIflK82sMnL/E8A7zOyLwFuBFTEvG+i5xDr5ZPjUp6CjAz74QWhtHbFDi4gk22jI2eXLg54ndQ/W0Xmw9wQmIiIjy9w92TGEYsmSJR7a4MVYTU1w7LHwxhvwmc/At74V/jFEJKOZ2Sp3X5LsOEZSmDn72WOfpWl1E8fcfwzjLxgfyj5FRPozUM5OhT7Wqa2oCH7+cxgzBr79bfi//0t2RCIiEqPiXcHMVzV3j6bJVUQkHamwHopTTgn6WwN86EOwenVy4xERkS6V7wp6qNTcU0Nnu7qDiEjyqLAeqn/6J7jiiqBryAUXwJYtyY5IRESAosVFFBxZQNueNvb+aW+ywxGRDKbCeqjMgovFnH467NgRFNd7lcBFRJLNzJj8wckA7Lp9V5KjEZFMpsJ6OPLz4Xe/g6OOgldegbPPhhr16RMRSbaJH5gI2VB7by2tuzWDk4gkhwrr4SovhwcegCOOgBdegLPOgj17kh2ViEhGy5ucx/i3jsfbnR1Vib04r4hIf1RYH46pU+GRR2DBAnjppaB7yMaNyY5KRCSjTfvsNAC2f3c7HQc6khyNiGQiFdaHa/Jk+NvfYPFieO01OPVUXZ1RRCSJxi4bS/GJxbRVt7HrJ+prLSIjT4V1PCZOhEcfDfpa794dtFz/6U/JjkpEJCOZGTP+bQYAW7+xVVdiFJERp8I6XmVlcP/98N73BlPxve1tcNNNkCZXtBQRGU0qL6mkcFEhLZta2HGr+lqLyMhSYR2G3Fz46U/hi1+Ezk74t3+Dyy6DhoZkRyYiklEs25jztTkAbL5+M+317UmOSEQyiQrrsGRlwQ03wN13Q2kp3HUXnHgiPP10siMTEcko4y8cT+mbSmmraWPzDZuTHY6IZJBQCmsze2sY+0kL73wnPPssHHMMvP46vOlN8OUvQ1tbsiMTEQHSP2ebGfO+NQ8Mtn5zK42rG5MdkohkiGEV1mZ2vpk9bmbrzWxDZNkE3JWY8Eap+fPhmWfgc58Luob8538GBfbq1cmOTEQySCbn7NKTSply1RTogHUfW4d3atyLiCTecFusfwb8Fvgo8MHIcmVkvcTKz4ebb4a//AWmTw9asU84Ieh/3dSU7OhEJDNkdM6e85U55E7Opf6perZ9a1uywxGRDDDcwnqDu3/T3R9290ciy9+AaxMQW3o488zgIjL/9E9B6/VNN8GiRXDPPZo5REQSLaNzdk5ZDvOr5gOw4QsbaHheA8pFJLGGW1h/wswu7WP9RWEEk7bKyuCWW+Cpp+C442Dz5qAv9tln66IyIpJIGZ+zKy6sYMpVU/A259UrXqWjSVdkFJHEGW5h/SbgR2bWEbN0Aj9IQGzpZ+nSoEvId74D48bBww8HM4d84AOwZUuyoxOR9KOcDcy9aS6FCwtpXtvM2g+uxfVtoYgkyHAL6+uAFcARwJyY5cvxBmJmC83se2b2WTP7vpkdF1lfYGbfNbOvmtkPzKw03mMlVU4OfOpTsH59MLgxNxd+9jM44gi46irYujXZEYpI+lDOBrILs1l05yKyS7Kp/k01W76mhgwRSYzhFtZ/BO5x9w3uvjmybAL+J4RYfgZ81d2/BXwduCOy/vvA4+5+DXAPcGsIx0q+8vJgcOOrr8LllwfT8f3gBzB3LnziE2rBFpEwKGdHFB1VxMJfLASDjV/aSPXd1ckOSUTS0HAL618BXzCz02OWM4AvhBDLQqAkcv8AUGZmU4ArgHsj6x8CLjKzWSEcLzXMmQN33AFr1gQFdns7/PCHMG8e/OM/wgsvJDtCERm9lLNjVLy9gtlfmQ0Or1z+CnUP1yU7JBFJMzacvmZm9izBV4r7gOgLDZjk7vlxBWL2E2AucAHwEWALkAfc5O5TY7ZbE1n3k9jXL1myxFeuXBlPCKnh1VeDKzj+6lfBLCIQzCzy2c/C294WXOFRRNKKma1y9yUJ2K9ydi/uzuv//Do7vreD7OJsjn34WEqXJL23ioiMIgPl7OFWaTcB0919lrvPjiyzgIvjDRK4CmgBngUa3f0uYCpQ22u7BmBK7xdXV1ezZMmSrqWqqiqEkJJg4UL4xS+CPtif+QwUFweDHC+6CBYsCLqPVOsrTJHRrqqqqitfARUJOoxydi9mxhH/fQQT3jOBjsYOVp+7mvqV9SNybBEZvYaas4fbYn2quz/Zx/oPufvthxVp9z4qgHcQJObPR+4fB1zm7ifHbPcc8Et3vzn29WnTYt3b/v1w223BTCLRftdjxsDFF8NHPxq0ZqsVW2RUS2CLtXJ2PzpbO3n5spep/X0t2aXZLL5/MWVvKktaPCIyeoTZYn1VHzu/DPj24QTWy8+Bu9z9OroHwuwEynttVwLsCOF4o0NZGfzLv8AbbwQXlXnb26CjA379a1i+PLh8+vXXw4YNyY5URFKPcnY/snKzWHTnIir/oZKO+g5ePPdF9j6wN9lhicgoN9zCepyZ/TuAmRWa2Y+AHxF8FXjYIi0fx7r7vsiqG4B6YDMwxcwKItuNAWYAj8RzvFEpJyfoDnLffbBpE3z5yzBtWlBwX3ttMJvIaafB978PNTXJjlZEUoNy9gCyxmRx1C+PYuI/TqSzuZPVb13NjltT7n8AERlFhlVYu/vbgAfN7LvAC8B84Figryt7DcdeoMXMpsasqwVWA/cBZ0fWLSdoIdke5/FGt+nT4T/+Iyiw//hHeN/7oLAQnnwyuHT65Mnw9rcHfbX37092tCKSJMrZg7NsY8HtC5hxzQzogHUfX8f6z63HO3QRGREZvgH7WJvZ6f08tRQ4E7gZ6ATe6e6fjSsQs2MJvrZcBUwEHnX3R8ysHLgR2ARMAr7g7o29X5/s/npJ19gYdBX5xS/ggQeC7iIQ9MdevhwuuQTe8Q6oSNQYKRE5XGH1sVbOjs/OH+9k3Yp1eJtTfm45C3++kNzK3GSHJSIpZqCcPVhh/RTBXKX7CJJxn5sBk909L95A45GKSTppdu+G3/wG7rwTHnuse9q+7Gw444ygyH7nO2HKIQP1RSQJQiyslbPjVPe3Ol75h1doq24jd2ouR/3qKMa+eWyywxKRFBJPYf1WYKW77xnkAJe6+53xhRmfVE3SSbdnD/zud/Db38Jf/hJcgCbqxBPhwguDAZEnnqjZRUSSJMTCWjk7BAe3H+Tld79M/d/rIRtmXzeb6f82nawc5UgRiaOwHk1SOUmnjLo6uPdeuOsuePBBOHCg+7mJE4MC+8ILg64jJSX970dEQpWo6fZSWarn7M62TjZ+cSNbb9oKQOkppSz4yQIK5xcmOTIRSbYwp9uT0ay8HD7wgaAvdm0t/OEP8IlPBIMhd++G228P5scePx7OOSe4GM0LL3R3JRERyRBZY7KYe+NcFv95MblTc6l/qp6Vx61k239v08BGEemXWqwF3GHNmmAqv/vug6ee6llMV1bC2WcHLdnnnAMzZiQvVpE0pBbr1Na2r431n1rP7p/tBqBkaQnzfzCfkhP0zZ5IJlJXEBmemppgZpGHHgq6jGzb1vP5I44ICuwzz4TTT4cJE5ITp0iaUGE9OlT/rprXP/k6rdtbIQumfnIqs6+bTU5ZTrJDE5ERlNDC2syOBz4HPE8wX+mmuHZ4mEZjkh4V3GHduu4i++GHob6+5zZHHRXMNhJdJk1KTqwio9RIFtbK2fFpb2hn05c3se0726ADciflMuv6WUz+4GQs25IdnoiMgEQX1rcC73D3SWZ2BXC3ux8Y7HVhG61JetRpb4eVK4Mi+5FH4Ikneg6ChOAy68uWdRfaU6f2uSsRCYxwYa2cHYLGFxtZ94l11D8ZNDQULipk7k1zGXf+OMxUYIuks0QX1mOACe6+3cyygGJ3rx/sdWEb7Ul61GpthWefDYrsRx6Bv/8dmpp6bjN3bneRffrpMHMm6A+PSJcRLqyVs0Pinc6e/9vDxi9spGVTCwBjzx7LnK/PoXRJaZKjE5FECbWwNrMLgHHu/gszWwo0uvsrIcQZl3RI0mmhrQ2eew7+9reg0H78cWho6LnNlClw2mnwpjcFt8cfH1whUiRDJbKwVs5OvM6DnWy/ZTubb9hM+77gWgHjLxzPzP+YqQJbJA2FVlib2VeATwN/cvdLI+u+Dqx291+GEezhSqcknVba24Mp+2JbtPfu7blNQQGcdFJ3oX3qqcGUfyIZIlGFtXL2yGqrbWPLN7aw/Xvb6WwOZlZSgS2SfsIsrB8G3gX8o7t/J7LuDOBn7p7UOdjSMUmnpc5OeO21oG/23/8e3L722qHbLVgQFNqnnAJLlwYDJHM08l7SUwILa+XsJGjd08rWm7f2KLDLzy1n+tXTKV9erj7YIqPcQDl7uJXKg+6+z8xiq/FLDz80yThZWbBwYbB8+MPBupoaePLJ7kL72Wdh7dpgue22YJvCQjjhhKDIPumk4Hb2bPXVFhmYcnYS5E7IZe6Nc5l+9fSuArvugTrqHqij6Ngipl89nQnvnkDWGF2jTSTdDLfF+nJgDjAVeAy4jKA15J/d/XsJiXCI0rn1I+O0tsLzzweF9jPPBMvGjYduN358z0L7pJM0p7aMSglssVbOTgFte9vY8cMdbPvvbbTtbgMgb1oeUz89lckfmsyYcRpjIjKahD148RTgSmAGsBO4w90fijfIeGVSks5INTVBS/Yzz3TfVlcfut3MmUHL9gknBIMiTzgBJk8e+XhFhiHBgxeVs1NER0sHe36xh603b6V5bTMAWflZTLh8AlOumkLpSeqHLTIaJHq6vWxgmbv/Ja4dxSkTk3RGc4fNm7uL7GeegVWrDp3qD2DixJ6F9vHHqxuJpJQRnm5POTvJvNOpvb+W7bdsp+7PdV3ri08sZupVU5lw+QSyC7OTGKGIDOSwC2sz+z1QPMj+K4Aadz/r8EOMXyYnaYno6Aj6ZT//fDDlX/S295UiAcaOheOO6y60jz8+uLCNpv2TJAirsFbOHn2a1zez89ad7Lx9J+17g6n6ssuymXD5BCZdOYnSk0s12FEkxcRTWN9AMMDxNaAT+CjBZXBjs+FsgjlSPxVaxIdBSVr61NkZ9M/uXWzv2XPotrm5wewjixf3XCZOHPm4JaOEWFgrZ49SHQc6qP6/arb/YDsNT3fP/V+4oJBJV05i4vsnkjclL4kRikhUPIV1WWSbfZHH17n7tX1sd4+7vyOkYE8B3gysBx4HmoAbgQagHPh8X1cJU5KWIXOHnTt7FtovvACbNvW9/YQJhxbbCxdCfv6Ihi3pK8TCWjk7DTSuaWTXj3ex+2e7adsTDHYkC8adN46J759IxUUVZBepq4hIshz2dHvuvr/Xqolmlu3uHTE7XwYsjTvKYF8fBOa6+5di1v0vwcUNfm1m5wO3Au8J43iSocyCqz9OmQIXXti9vr4e1qyB1at7Lnv2wEMPBUtUdjYceSQcc0zPgnv6dPXdlqRRzk4PxUcXM+/mecz52hz2/nkvu/53F7X31rL3j3vZ+8e9ZBVmMf7t45lw+QTGnT+O7HwV2SKpYrjT7Z0DfIugVaIBWAicB1zr7l+LK5DgogX/DpzjkaDMbAqwESh392YzywH2A4vcfVPs69X6IQnhDlu2HFpsr1sXdDPprbQ06E6yaFGwHH10cDt5sgpu6VcCp9tTzk4TrTWt7LljD3t+tYf6J7q/AMguzaby4komXD6BsWeN1dzYIiMg7On2pgBXAYuBeuA37n5PGEECTwAGzAOuB2YCN7n71Jjt1kTW/ST29UrSMqIOHIBXXoGXXupZcPc1BSAEgyV7F9uLFgXdTFRwZ7wET7ennJ1mWja3sOf/giK78bnGrvU543IY//bxVL6rkvJzy8kuUEu2SCIkerq9LOB4d18Vxz6OBF4FjnP31Wb2r8C1wA3Ae919ccy2TwK/793aMnPmTK+srOx6vGLFClasWHG4IYkcnj174OWXey5r1kBdXd/bjx9/aLG9aBFUVIxs3DLiqqqqqKqqAmDVqlWb3X3WSBxXOTu9NK9rZs+vgiK7+dXmrvVZhVmMO38cFe+sYPyF4xlTrhmPROIx1Jw93K4gt/exuhLocPd3Hkac0f1eDFS5e0XkcQFQDVwHXOLuJ8ds+xzwS3e/OXYfav2QlOUOu3YdWmy//HLfUwFC0JLdu9heuDAoxCXtJLAriHJ2Bml6tYmau2uo+V0NDc92zyxiOcbYZWO7iuz8mRp4LRKPwx682IcjgL8AsdX4NIJL5cYjJzYWdz9gZq8D2QSjymOVADviPJ7IyDEL+lhPngzLl3evd4ft2w8ttl95JWj5/utfgyXWhAlBH+6FC3veTpqkLiXSF+XsDFK0sIiihUXMvGYmLVtbqLmnhpq7a9j3yD7qHqqj7qE6Xv/k6xQeVcj4t45n3FvHUfamMrJy1S9bJCzDLaw/4u6vxa4wswUEl8uNx2qgzMwmuvvuyLp2gmQ8xcwKIol7DMFleR+J83giyWcG06YFy3nnda/v7IStWw8ttl99NSi49+yBv/2t577KynoW2tH7M2ZAlv5oZjDl7AyVPz2faZ+cxrRPTqNtbxu199VSc08NdQ/W0fxKM82vNLP15q1kl2RTfk55UGhfME5zZYvEKYw+1rOBZ6NfCcaxn/uBP7j798xsLEEiPhH4OfBzd7/PzC4A3u/uV/R+vb5WlLTX2QnbtnUX2bG3/fXhLiyEBQsObeWeOxdyhvt/tSTKCF/SXDk7g3W2drL/if3svX8vtffX0vxyc4/ni48rpvy8csqXl1P2pjINgBTpQ2iDF81sIz2/UswGJgI/c/ePxhlkBfAdgiuEzQD+x91fMbNygosNbAImAV9w98ber1eSlozlDrt3H1psv/JKsL4vublwxBE9W7ePOiq4rHueWqxGWgL7WCtny4BaNrdQe38te+/fS91f6ug80D2NqOUZZW8uo3x5UGiXHF+CZavLmUiYhXUV8Eu6E3UnsNvd18UdZZyUpEX6sHdvUGj3Lrq3bOl7+6ysoDW7dx/uBQuguHhkY88gCSyslbNlyDpaOtj/6P6u/tiNz/f8fyhnXA7lZ5V3Fdr5c/Ixje2QDBRmYX26uz/aa92RBCPM18cXZnyUpEWGobER1q49tFvJG2/0feEbCPpr9zVwsrz3WDUZrgQW1srZcthaq1vZ9/A+6h6sY++Dezm4+WCP5/Nm5jH2jLFdiwptyRRhFtZfdfdreq0bCzzs7sfHF2Z8lKRFQtDSAq+/3t2VJFp0r1sHbW19v2bSpO4uJbEFd2WlZioZogQW1srZEgp3p2VDS1drdt1f62jf295jm9ypuYw9PSiyy84oo/DIQhXakpbiLqzN7JMEE/+X9LPJg+5+/uGHGD8laZEEam8PWrN7dyl59dUmktE2AAAgAElEQVTgKpR9GT++Z8EdXXR590OEXVgrZ0uieYfT+FIj+x/Zz75H9rHv0X201/YstMdMHNOj0C46qgjL0u++jH6htFib2WLgH4H/7vVUq7vvjC/E+ClJiyRBZ2fQXzvawh27NDT0/Zro1IC9pwecPj1jpwZMRIu1craMJO90ml9tDorsyNK2u+e3XDljcyg9tZTS00opO62MkqUl5BRrdiIZfcLsCjLG3Q/5PtjMZrh7P6OhRoaStEgKiV78pneXkpdf7n9qwKKi7kL7mGNg8eJgmTgx7Vu4E9gVRDlbksLdObDuQFdr9v5H93Nwa88+2mRD8bHFlJ1W1lVs583IU/cRSXmHXVib2dXA/e7+SuTxFwmma+qxGXCmuy8LJ9zDoyQtMgq4Bxe46d26/eqr/U8NWFnZXWRHl6OOgvz0uSxzWIW1crakspatLdQ/Uc/+J/ZT/0Q9Dc83QEfPbXKn5vYotIuPLSYrLzO/yZLUFc8lzWcDY2MenwAUAbvonr7JCOYqFREZmFnQAj1xIpx5Zs/namuDAnvNGnjpJVi9Oliqq+EvfwmWqKwsOPLInsX2MccEM5dkdmuXcrakrPzp+eS/O58J754AQEdTB/XP1vcotlu3t1L9m2qqf1MNgOUaxccVU7q0lJKlJZQuLaXgiAL11ZaUNdyuICcAa9y9tdf6o6ItJMmi1g+RNOQe9OGOLbRXr4bXXut7WsCysp7dSBYvhqOPhpL+xvClhgR2BVHOllHDO53mdc3U/z1SaD9ZT/Pa5p6XOCLoq11yUklXoV2ytIS8SbqwlYyceFqsezvG3Z/rtfMjgblAUpO0iKQhM5g5M1guvLB7/YEDQet2bLH94otQUwOPPx4ssebM6S60lywJlsmTR/ZckkM5W0YNyzKKFhRRtKCIyR8Ofj/b97fTsKqB+mfqaXi6gfqn62nd2Urdg3XUPdg9XiNvRl6PVu3iE4rJKdHASBl5Q/qpM7O3EHx9eE7kErmxxhGMOr835NhERPpWUAAnnBAsUdFLu8cW26tXB324N2wIlt/9rnv7qVO7i+yTToITT4SKipE/lwRQzpZ0kVMWudrjWd0Xojq4/WBQaD8TKbifbeDgloNUb6mm+s6gCwkGBfMLKDmhhJITSyg+oZiSE0rIKVOxLYk11J+wfcBPCfrvvanXc63Aj8IMSkRk2MyCi9VMmgTnntu9vq0tuMDN6tXw/POwciWsWhXMWrJ9O9xzT/e28+fD6ad3LzNnjvx5hEM5W9JW3tQ8Kt9VSeW7KoFgTu3m15p7FNtNq5s48NoBDrx2gD137Ol6bcG8gqDIjim2x4wbk6xTkTQ0nHmsy4DL3D0lE7L664nIkHV2wvr1QZH97LPB7XPPQXNzz+1mzgwK7PPPD5Zx4xISToLmsVbOlozV2dpJ05omGp5roHFVY3D7YiN+8NCaJ39WPsUnBsV2yQklFJ9YTG5FbhKiltEitHmsBzjAm9398cG3TBwlaRGJS3t70KL9yCPw6KPw2GOwb1/381lZcOqpQV/vd78bZs8O7dCJGrw4wPGUsyXjdLZ10vxKMw3PNdCwKii4G19spPPAoQOh86bnUXx8McXHFlN8XLDkz87XHNsCxDeP9e+B4sH2Dyxy9wmHH2L8lKRFJFSdncHUf3/5C/zhD0Gx3RZzrZXTT4cPfACuuCLo8x2HEOexVs4WGYbO9k6a1zbT+FwjDasiBfcLjXQ2HVpsZ5dmU7y4u9AuOraIoqOLyM7vPVW8pLt4CuubCPrqbeeQCW+6NwMudPdL4w00HkrSIpJQ9fXw4IPw29/C3XcHM5NAMCf35z4Hn/oU5B3elF8hFtbK2SJx8o5g2r/GFxtpfKGRphebaHyhkdZdrYdunA2FCwqDYjumdTu3Ul1J0lk8hXUF0ODuB/vdKNhuvLvXxhdmfJSkRWTENDTAXXfBd78b9M0GWLQoKLiPOGLYuwuxsFbOFkmQ1t2tXcV24wtBN5Lmtc3Qx5T6uVNyexTaxccWUzCvAMtWV5J0MBJ9rD+S7AEyStIiMuLc4c9/hk9/Oph5ZM4cWLsWxgxvloEk9LFWzhYJQceBDprWNPUouJtebKKjseOQbbMKsyg6pigouCNL0eIizbc9CoV5gZjeO54B/BPwSTR9k4hkGrNgtpDHHw+m+duwIZg3+9hjkx1Zn5SzRcKVXZBN6UmllJ5U2rXOO52WjS09WrYbX2jk4NaDNDzdQMPTDT32kT8nv7vQPjYovPNnaaDkaDXswtrMsoGLgI8C5wC7Cfr0xc3McoFngU+7+9/MrAC4EWgAyoHPu3t9GMcSEYlbdTX84hfwrW8Fgx0nTz6sriCJlMicHdm/8rZIDMsyCuYWUDC3gMpLKrvWt9W20bg6KLSbXgxauZtebqJlQwstG1qoubuma9vs0uxDiu2io4vILtBAyVQ35MLazGYDHwE+CIwHNgNvB/4MnBpSPP8KzIp5/H3gT+7+azM7H7gVeE9IxxIRGZ7W1mBKvgcfhD/+EZ56KiioIbhc+h13QGFhcmOMGKGcDcrbIkMyZvwYys8sp/zM7qtIdrZ10vxac1ehHV3adrex/7H97H9sf/cOsqBwfmFXoR1dcqfkqnU7hQxaWJvZZQQtHWcB64FvAj8BrnH3P0U2eyLeQMzsNGAnUBd5PAW4guBrS4CHgLvMbJa7b4r3eCIiA4pesfGll+CZZ4Ii+rnn4GDMuMCcnKAryIc/DO94B2QnvzVppHJ25FjK2yJxyBqTRfHRxRQfXczE907sWn9w18FDiu3mtc1dS/Wvq7u2zRmf07Pf9rFFFB1VRFZuVjJOKeMNpcX6AMElcJ8APuLu6wDMLP5RjxFmVgRc6u7/YmbXRlYvA2rcvRnA3dvNbCNwBrAprGOLSIZraoI33giuxLh2bTB39Zo1wf3YeaujFi6EN78ZLrgAzj4bSksP3Sa5Ep6zI/tT3hZJkLxJeeRNymPced1Xe+1o6aD5leauQjtaeLfXtrPvr/vY99fuHl6WYxQuLOzuShKZmUTTACbeoIW1u98H3BcZ9PJRM5sF3AN0DXs3M/P4phf5PPD1XuumAr2ng2oApvS1g+rqapYs6R6guWLFClasWBFHSCKSFlpbYft22Lo1WDZu7C6k16+HXbv6f+2cOXD00XDiiXDKKbB0KYwdG0pYVVVVVFVVRR9WhLJTRixnQ5x5WzlbZHiy87MpOSG47HqUu3Nw28EehXbji40ceP0ATS810fRSE/y8ex+5k3N79NsuPraYgvkFZOWodXswQ83Zw55uz8xygHcC7wXWALcB57p71YAv7H9/5wF57v77yONNwJXAEuAydz85ZtvngF+6+82996Opm0QyUGNjUBjv3Bncbt8OW7YEBXT0dteuYFq8/uTmBgX03Lkwf35QSB9zTNAyXTzYRQzDkcjp9sLO2ZF9xp23lbNFEqejqYPGl3oW202r+5kGMD+LwkWFh3QnGTN2eNOGZpJQp9tz93bgTuBOM5sHfJagP9/hJumrgeNjOt6XE7SufDNyP1YJsOMwjyMio8HBg1BTA3v2BEVxbOHc+35j4+D7y8qCqVNh+vRgmTkT5s0LlrlzYdq0lOgbnSgJyNmgvC2S0rKLsik7pYyyU8q61nVNAxjTb7vpxSZaNrXQuKqRxlU982nejLyeM5McV0zBnAIsSwMlBxLXPNbuvh74tJmtimM3VwCx1wF+EvgX4GngX82swN0PmNkYYAbwSBzHEpGR5A779weFcnV1923s/d63DQ2D7zcqPz+Y4m7SpGCZMgVmzAgK6OjtlCnBIEMJK2eD8rbIqNNjGsCLu6cBbN/ffug0gC81cXDLQQ5uOUjtvd29u7KLs7suctM1DeAxReQUK8dGhfJOuPtP43htdexjM+sAqt19i5ndB5wN3AcsB+5y9+1xBSsiw+cetA7v3dtzqas7dF10iRbL7e3DO1ZODlRUBMvkyT0L5+j96G1paXCRFhmWeHJ25PXK2yJpIqcsh7FvGcvYt3SPH+ls7+TA6wcO6bvduqOV+ifrqX8yZmp6g4K5BRQfV0zJ0hJKTy6l+ITijC22U/2sPwHcaGbHApMAjWwRiUdLC+zb1730Loz7K5Tr6oZfIEcVF0NlZVAoV1b2vN/7trISyspULI9uytsio1xWThZFC4soWlgEl3evb61pPXQawFeaObD+AAfWH6D6zsj/3FlQtKgoKLSXllKytISio4syYpBkyhXW7j4r5n4dQV9AEXGH5uagIN6/v2eBHF36Wx99LnYO5uEqKoLychg3rufS37powZyfH957IClJeVskM+RW5JJ7di7lZ8dc5Ka1k+a1zTQ8F1yuvf6ZeppWN3XNSrLrtmDmpayCLIpPKKZ0aSmlp5RS9uYy8qbk9XeoUSvlCmuRtNXZGXSniC2AByuEe6873FbjqDFjgqJ37NigZbi/wrj3+vJyyEu/BCgiIvHJys2ieHExxYuLmXzlZAA6DnTQ+EIjDc8EhXbDMw0cWH+A+r/XU//37m4k+XPyKXtLGWPfMpayt5RRcETBqL+KpAprkaFqbe1Z7PYujvu7jd7fv3/gad+GoqAgKIjHjj106W997HP5+epmISIiCZVdkE3ZqWWUndo9K0lbbRsNKxuof7qe/U/sp/6Jelo2tNCyoYXdP9kNwJgJYyh7cxljzxhL+TnlFC4oHHWFtgpryQzRwXcDFb6D3R44EH8cxcXdRe5QCuHe69RqLCIio9CY8WMYd964rqtJeofT+GIj+x/bz77H9rH/sf207Wmj5rc11Py2BoC8aXmUn1vOuPPGUX52OWPGp/7c2iqsZXSIthYPt5U4trW4szO+GHJyDi2KYwvg/p6L3paWato3ERERwLKt60qS0z49DXfnwOsH2P/YfuoerqPugToObjvIrtt3sev2XWBQsqSE8ReNp/KSymBgZQrSX3lJPHdoaoqvtbi5Of44ioqGVgD391xhobpRiIiIJICZUTi/kML5hUz+8GS802lc3Ujdn+vY+8Be9j++n4ZnG2h4toFN/76JwgWFVFxSQeUllRQfV5wyXUZUWMvgooVxbW331GvDbS3uOPQyqsOSnT20ori/4ri0NBi4JyIiIinPsoyS40ooOa6EGZ+fQUdTB3V/rQu6itxTQ/PaZrZ8ZQtbvrKFomOKmPzhyUx47wRyK3KTGrcK60zT1tZdINfW9lwGWtfaGt9xi4qGXhT3VRyrtVhERCRjZRdlU/H2CireXkFnWyf7HtlHzV017PnNHppeamL9Z9bzxr++QeUllUz/t+mUHF+SlDhVWI92Bw/Cnj2we/eht7t3B1e+iy2Uh3O56FgFBTB+fLBEp2sbanGs1mIREREJSdaYLMYtH8e45eOY95151N5by87bdrL3z3vZ86s97PnVHsrPLWfWf8yi7LSywXcYIhXWqaijIyiOt28Plh07+i6c9+wJulkMR1ZWMC9xtEiOLn2ti32uoCAx5yoiIiJymLJys6i8pJLKSypp2dLCtu9sY8etO6h7IBgAOeE9E5jzjTnkTx+Zi5WpsB5pnZ1BobxxI2zbFhTO0dvo/Z07h34hkJwcmDAhWCZO7Hk7YUL31e+ihXNZWVBci4iIiKSR/Bn5zPuvecz84ky2fWsbW2/eyp479lD7h1rm3zqfiZdPTHgMKqwToaUFXn8dNmzoubzxBmzaNLTLSldUwLRpMHUqTJkCkyb1XTyPHatCWURERCRizLgxzL5+NpM/MpnXP/06tffU8up7XqXhmQbm3jwXy0rcmC0V1vFob4e1a2HNGnj55e5l/fqB50yeMAFmz4bp04PCOVpAxxbS+SPzlYWIiIhIOsqfmc/Rdx/Njlt3sP5T69n2rW10NHUw/4fzEzY9nwrr4dizB/7+d3jqqWBZubLv+ZWzsmD+fJg3D+bM6bnMnh1cfU9EREREEsrMmPrxqRTMKWDNO9ews2onRUcXMe2fpyXkeCqsB9LZGRTQ998Pf/wjPPfcodvMng2LF8OiRd3LkUeqxVlEREQkRYw7dxwL/ncBr1z+Chs+v4GKd1SQPyP8Wk2FdV927IAf/Qh+/ONgkGFUfj68+c1wyilw8snBUlmZtDBFREREZGgmvHsC1XdXU/3rajZdt4kFP1oQ+jFUWMeqqYHrr4dbb+0eYDhtGlxyCVxwAZx+uqadExERERmlZl83m+pfV7Pnjj3M++Y8ckrDLYVVWEf96U9w5ZXB/NAQFNMf/ziceWZwOW0RERERGdUK5xdSelop9U/Us+9v+6i4qCLU/afMPG1mdrGZbTSzWjP7jpnlRNZXmNmtZvYNM/ummYV/Cb9f/hIuvDAoqk8/HV58Ee68E5YvV1EtItKHpOZsEZE4lC8vB2DfI/tC33dKFNZmNgN4J3Ap8FngI8BnIk//Brjd3T8PrAOuC/Xgq1YFLdUdHXDNNfDww8FgRBER6VNSc7aISJxKTigBoPnVPmZ2i1OqdAWZCXzI3duBVWZ2NHCmmT0OHOnuT0e2uxd4zcxucPemuI/qDp/8JLS1wVVXwVe+EvcuRUQyQHJytohICArmBePlDrx+IPR9p0SLtbs/FknQUduBLcBZwMaY7bZH7i4J5cBPPhlMp1dRATfeGMouRUTSXdJytohICPKm5wHQuqs19H2nRGHdh5OAHwJTgdpezzUAU3q/oLq6miVLlnQtVVVVgx/ld78Lbt//figqijNkEZGhq6qq6spXQLijZ0beyORsEZEQZBdnQxZ0NHbQ2TbAlbJjDDVnp0pXkC5mNg+odvcXzcyB3u30eUBb79dVVlaycuXK4R3sqaeC2/POO5xQRUQO24oVK1ixYgUAZlaT5HAO24jmbBGREFiWkTM2h/a97bTvaye3MnfQ1ww1Z6dUi7WZZQMfA/5fZNV2oDzmeQOKgB2hHHDNmuBWgxVFRIZtxHO2iEhIcsqCtuWO+o5Q95tShTVwNXCzu0euzsKDwJExz88CWoH4mzmam6GuDvLyYNKkuHcnIpKBRi5ni4iEKCsvKIGH2hVkyPsNdW9xMLNrgOeBAjObY2YfAvYBOyMjzgHOB77r7vH3Nq+uDm4rK8Es7t2JiGSSEc/ZIiIhstyg9vODHup+U6KPtZl9Cbi+1+q17n67mb0buNbMNgHFwBdDOWhtZHxNxWgfMyQiMrKSkrNFREKUlRtpsW4Nt8U6JQprd78BuKGf5zYDHw79oAci42sKC0PftYhIOktKzhYRCVFXi3VruC3WKdMVZMS1tAS3eXnJjUNERERERlSiWqwzt7A+GBlrk5+f3DhEREREZERFW6w7D6qwDodarEVEREQykmVHJq4Id7a9DC6sWyOD1HMHnxRcRERERNJHtLD2TvWxDkdH5F+UnJQYvykiIiIiIyVaAYfbEySDC+vOyDuZlblvgYiIiEgmsqxIi3WHWqzDocJaREREJDNlBzfqChIWFdYiIiIiGSnaYq2uIGFRYS0iIiKSmSLln7qChCU6eFGFtYiIiEhG6ZpuTy3WIYm2WGdnJzcOERERERlRGrwYNnUFEREREclMGrwYMhXWIiIiIhlJgxfDpsJaREREJDNp8GLIVFiLiIiIZCQNXgybCmsRERGRjNQ1eFF9rEMSLazNkhuHiIiIiIwsdQUJmUfeSBXWIiIiIhklUV1BcsLdXWKYWQFwI9AAlAOfd/f6uHaqwlpEJCESkrNFRMKU4S3W3wced/drgHuAW0PbswprEZGwJS5ni4iEIGMHL5rZFOAK4N7IqoeAi8xsVlw79nD/QxERkQTmbBGREGXylReXATXu3gzg7u3ARuCMuPaqriAiIomwjETkbBGRMEXLv5DbWUdDH+upQG2vdQ3AlNgV1dXVLFmypOvxihUrWLFiRf97LS6GyZOhtDS8SEVEhqiqqoqqqqrow4pkxhKyxORsEZEQ5ZTlkDs5l6yiobUxDzVnm6d4lwgzuxq4zN1Pjln3HPBLd785um7JkiW+cuXKZIQoIhIXM1vl7ksG3zL1KWeLSLobKGePhq4g2wlGlccqAXYkIRYRERmYcraIZKzRUFg/DEyJTN+EmY0BZgCPJDUqERHpi3K2iGSslC+s3X0XcB9wdmTVcuAud9+evKhERKQvytkikslGw+BFgE8AN5rZscAkQCNcRERSl3K2iGSklG+xBnD3Onf/qLt/xd3/2d0bw9p3zAjPtKbzTC86z7STTrOCKGeHSOeb3nS+6WdUFNaJlAkfMug8043OM+1UJjuA0SKDfiYAnW+60/mmn4wvrEVEREREwpDy81gPlZlVA5sP46UVQE3I4aQinWd60XmmlyPdvSTZQYwk5ewh0/mmN53v6DTT3fv8pjFtCmsRERERkWRSVxARERERkRCosBYRERERCYEKaxERERGREGRsYW1mBWb2XTP7qpn9wMxKkx1TWMzsHDN72sxmxaxLq/M1s4vNbKOZ1ZrZd8wsJ7K+wsxuNbNvmNk3I5dTHrXM7Awze83M6szs2zHr0+rzjDKzXDN70cyWRR6n3Xma2Q1m5pHlxci6tDvPeAz1/TCzC8zsFjP7vpldOtJxhmUY5zvWzOpjfn7eNdKxhqWvv1N9bJMWny8M+XxH/efb39/mPra70sy+ZWY/MrPTRzrOhHL3jFyA/wXeHbl/PnBHsmMK6bwmAu8CHJiVjucLzAB+CpwIfABoAq6OPPcwcHLk/seBryU73jjOswy4BigH3ga0AcvT7fPsdc5fBPYDy9LxPIEi4IcEl/leDhyRjucZwvs06PsBzAeeA7Iijx8Fjk527Ik638hz/xrJ79Gfn+xkx36Y59vn36k0/nwHPd90+HwH+tvca7szgfsj98cAq4GJyY4/rCUjZwUxsynARqDc3Zsj/1HtBxa5+6akBhcCM8sCOoDZ7r4p3c7XzN4CPOnu7ZHHNwKLgOuB37r7lMj6qcBrBL+wTcmK93CZWQHQ4pFfUjNbBVxNcE5p83lGmdlpwALgWuBKYB1pdp5m9pnI3R+6e0tkXVr9fsZrqO+Hmf0Q2Ovu10QefxGY5+4fTELYh20Y55sH3Adc5e6vJyXYEPX+O9XH82nx+UYN4XxH/efb399md39br+3+RPC3uiry+H+Abe7+nyMdcyJkaleQZUCNuzcDRH4INgJnJDOosLh7Z69Vy0ij83X3x6K/uBHbgS3AWQTnFd1ue+TukhEMLzTufiCmqC4CXgX+Rpp9ntB1fpe6++0xq5eRRudpZgZcDvwXsNvM3hN5ahlpdJ4hWMbQ3o8ev+8E/3AuG4H4wraMoZ3vhcDJwDoz+4OZVYxolCHr4+9Ub+ny+QJDOt9R//kO8Le5i5llE/xsp81n21umFtZTgdpe6xqAKUmIZSSk+/meRPD1elqeZ6QV4I8EX5kVkp7n+Xng673WpdV5euAUggskfBf4hZldQJqdZwiG+n703m60vmdDOl93v4uge9gZwCzg95FW0HSVLp/vkKTp5xv92xxrHJBPGn+2o/1DO1wOHOi1Lo+gD2s6StvzNbN5QLW7v0j6nud64MfA2cBNpNl5mtl5wEp339PrqbQ6zyh3r3P3LwFfAz5Dmp5nHIb6fvTebrS+Z0P+/CP/nD1K0Jq7gKCFM12ly+c7ZOn0+fb62xwr2v84bT/bTC2stxMMCItVAuxIQiwjIS3PN/KV0seA/xdZ1eM8I1+9FzHKz9Pdd0a6SFxN0JqRbp/n1cDtZlZjZjXAdOAegs8unc6zt/8mGOyTbp9nvIb6fvTebrS+Z8P+/N19N/Brgp+fdJUun++wjfbPt4+/zbFqgYOk8WebqYX1w8CUyOAwIlOyzQAeSWpUiZOu53s1cLO7H4w8fhA4Mub5WUArsHKE40qUlQR/bNLt87wCOC5m2QF8hGCmhHQ6z946CWY9SLfPM15DfT96/77Pj7x2tDncz78DeCHBsSVTuny+h2s0f769/zZ3iYwbeog0/mwzsrB2910Eo2/PjqxaDtwVM9htVIu01AIYpOf5mtk1wPNAgZnNMbMPAfuAnWZ2dGSz84HvuntrsuKMh5kVmtnxMaveCnwn3T5Pd692923RheAPSrW7byGNztPMppjZFWaWFfkd/RzwxXT7POM10PthZlebWfQP8g+B86Ar551F8C3AqDLU8zWzS81sTuT+AoIZM15LStAh6P13KrIu7T7fqMHON10+377+NpvZEWZ2nZlNimx2C8EUsphZLsH0fD9KTsThy8jp9gDMrBy4EdgETAK+4O6NSQ0qBGZWArwP+D5wHXCLu1en0/ma2ZcIptaLtdbdF5rZTILp2jYBxQSFSzujkJktBf4ArAWeAp5x999Enkubz7M3M9sEXOnuf0un8zSzRQSDUOuBx4Fvu/vayHNpc55h6O/9iEw5+TV3vzOy3fuAEyIve8Td70lGvPEayvma2S0Euf2vBLMD3TKEmSZS0gB/p9L18x30fNPh8+3vbzPBZ/gK8B53fyqy7eeA8QTdQH7l7n8fyVgTKWMLaxERERGRMGVkVxARERERkbCpsBYRERERCYEKaxERERGREKiwFhEREREJgQprEREREZEQqLAWEREREQmBCmsRERERkRCosBaRlGJmlWb2ezPbbGb/L9nxiIjIwJS3u+kCMSKSUszsPOBhoJTgql3z3X1vcqMSEZH+KG93U2EtIinLzB4Fzhqtl6UXEck0mZ631RVERFKSmU0HHsjU5CwiMtoob6uwlgQws38zs+8lO46oMOIxs+Vm9qSZLQsprLRiZjeHvD8D/gH4epj7FZGelK/jZ2azzOx/zOzWkPaXbWYfNbPNQ9j2bWb29OGea5i5W3k7oMJ6lDGzhWb2WzP7kZm9aGZuZu9Ldly93A/ckewgYoQRz1bg2BBiSUlm9lYzqzGznWZ2qplNj/ycuZndYmZjI9stNbM1ZvZfZlYQWTcR2BVySFcAt7l7e2T/IqOO8vVhGY35ek/kNi+k/WUB+4AZQ9h2JbCwvyfNbJGZ/fidzBMAAAkNSURBVDrys/clM5sU81zYuVt5GxXWo9HvgP90948AxwM/jT5hZovN7IxkBBV7bHdf4+6PJyOOvgwnnv7eQ3d/DagJM64wP6949+Xu9wPfB9rd/Ul33wp8FGgFnnb3fZHtniFI5P/q7gciL387cF9cJxDDzP4d+DLwkJmtBY4Ma98iI0z5ephGY75292ZgR1jHcvc2YNUQt91NUIT39/zLQFXk4c3uHltIh5a7lbe7qbAeRcxsHDAfaAJw907gWsDNrIwgaVsS4krascM0hPPoHMFjJWNfvwammdmJAO5eC/wFuDDmWGOBHZGfvagF7r42zmN3cffr3f0Id1/i7gvc/dGw9i0yUpSvEysF83XYM0EMJ/7Btu2I3Pbu9xxa7lbe7qbCehSJTF3zGvAnMzslsm4z8AxBv6bZwAoz+yCAmV1nZv8e+Rqoqx+VmS0xs2+a2WfMrNnMXjOzfzGzd5rZs2b2LjO728zqzexjMa/rc3+9jv1RM/txbF+zyNeh3zWzr5nZw2a2NLJ+wOPFvN7M7GNm9rKZnWNmq81st5ldOYRjHBEbzyDHPOQ97EN55KvdBjP7aK/39MuR9+a3kUQ81PfskGOZ2aVm9oyZXWlmr5vZ1wfYX1+ffZ/xDCTSsvEycFnM6krgbRbp9gFcBPw+Js4CoDnm8QIz+7qZ3Wxmj1jQpeTdZrbSzN5hZneY2T4z+7SZnRg5x1oze+tg8YmMJsrXmZOv+9LfOQ50HDObbWa3W9D6e9NQjtPHce83s1+Y2ZsG2U65O1HcXcsoWoBFwGaC/1BvBybEPLcJWBa5fzSwIXK/kuC/6fEE/3FvAZZEnvse8ETkfi6wF7iFYC7K9wDbB9pf72MT9DG7FfhxZH0+QbE2NvL4UqAaKBnoeH2c96TIMa+JbHsD0ALMHOQYveMZ8Jix72EfMWwC/gsojLxuS2R9KfCrmO3+TpAUh/Se9XOscQRf730bWAosH2h/vT77PuMZ4s/Xl4A3IvfnAl8DDgIXx/y8WMz2FwFvidwvBJ4FCiOPH428D6XAfuA6oAB4d+TxPxD8c/85Ij+DWrSk04LydUbk68jzXx7i+zjQZ/0ssDDmNT7En7Po53kq8B99PL8scpycmHXK3Qla1GI9ynjQqngMQaL5APCymfU1SON14P1mNgZ4S2RdMcEv8nQg2j/2VYJfGNy9FWgA7nT3eoJftEmD7K93fAfpORjibcABj/TRBX5L8At58SDH673f6D5/Gdn2BoL+v+cMcowe8QznmP2414P+dM8Ck2POcXykRekzBO8pDPE96+d89xIksHvc/Rl3f2gY++svnqH4NTDHgu4glwD/A/wZ+AczKwYaPJJ5I94EPBG5fxHBH69oK8i7gOsj7/M+4K8e9MteCZS6+/958PX48wzvMxAZFZSvMyNf96HfcxzgOMuBMnePxrNymMc8CXivu//nELdX7k6QnGQHIENnZtnAdHffBHzKzG4jGEFdBZwcu627HzSzqcDngZ9Ed+Hue8xsJXA2wX/UswmKqa6XxtxvJ9JdqL/99RNq7D6OAMbExNVpZpuAaQMdbzDu3mJmG4CxQMUwjnHYx+ylne7fnxnARnf/du+NhvKemdn7CQrYqCM9+MrYY2MdxmfQbzyDcffXzex5gu4gZe6+wcx+TdCKdBnwh5i4oz8b0f57MwlamKL7qo3ddcz9DnrqRN3SJM0oX/c4v0zI17H6fR8H+GwWE/wTcbjOAo4ws2J3b+xnm87IOSh3J5DekNFlDDH9X939RYKv2hb03tDMzgKucvcbPJjhIdY/ASdZMO3TauAbgx14kP0NZAtBC2jsNERG0PcwXrnAugQfYyh20rMfMmZ28jDes98Dx8UsfY4uH8b++oxn6KfD/2/vXl7kqKIAjH8nA+MighokkBkV3QkquBDciKKEiEaYjS8QwYiQIJFZiAthIOA6s1CzkIiQhRgRhSCCAdGNzz9AVBSyMS5FfOAis7guzm2903R1t0O1k7a/HxR0VXXVvVXVfeb2fc3bwBHgQpO/JWCdbDYduJPsLzrwE3BXROxt0h3bz0/6HzNeb7dI8brzGsek81s9ZpmdOUnW9r/csf/X8s+gc2P3DFmwnj/PRcRKs74KnK+vL5GDNW4mp3a6KiKuiIiDdf/+iLiW7Ov3BtmM8wW1abHaw9Cv9IiICecbTjuac5wD/gAer+e6uqY3qPnsSq/Lgfqe1Xrs+SnSaPMzKc32OoaNqr0IshZqL/BBRNwfOUhm0jP4Oy3gQCnlu2bZapJYal6PO1+b75H5iYjliDhdu3mM8w7ZBP0uQCnld+BD4POyfTaQQ2Q3kYHB/T4bORf2cWAwl+nwMxh+znM9Q4HUwXjNwsTrae/jyHTIWLpM/viCbJ0gIlamiN1B1h4/CTwSEY/WY5+PiA2ytv7j5v3G7lmaRcdtl9ks5ICIAvwCvAmcIZuSrqn7TwAXgTWyaecHsobgMFnTcYYMUu+TAyku1fP9SX7RHgC2yAEy1wMv1v1PjzvfUNrrwKdkv7Vb6r47yJrOTbKv4a11e2d6Hddf6vHrZNeE25p9XWnc2OZnUprtPRxKe40cxHeq3ovBcUfq/nuAr+uz2az3edp7tjbiWh+u6Z0FVuu2cc9027k68nMd+Y8M3pris/bK0PpjwL1D206OOO5u4BvgZ2CjbnuQ/Ky9CqxQpxwDjpF/eF+v1/rQbn/HXFz6WjBeL1K8vgH4jKx1v33CNY571ofJlsKvyIHk3wPP1msfGbvJVpEtsqZ6H/AJWah/BjhKzue9SXZDGRxj7J7hEvXmaUFExH7gWCnlpbq+h6xFOVpK2djVzE0QEQW4qWSfRe1QRLxQStnRVE6S/jvGa7WM3fPBriCL5zjZLWBfs+0+mkFpl6Om6clmpx2KiKWIeAJ4b7fzImkqxmsZu+eMBevF8xrZF+vbiPgR+Ai4WEr5cnez1a0O5livq09FTvumf+9K4Fwp5cLEd0q6HBivBcbuuWJXEEmSJKkH1lhLkiRJPbBgLUmSJPXAgrUkSZLUAwvWkiRJUg8sWEuSJEk9sGAtSZIk9cCCtSRJktQDC9aSJElSDyxYS5IkST34C2IBYTB91FAsAAAAAElFTkSuQmCC\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(vehicle1.v_kmsc, vehicle1.h_kmc, 'k-', linewidth=2.0)\n", "\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(vehicle1.acc_net_g, vehicle1.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(vehicle1.q_stag_total, vehicle1.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(vehicle1.heatload/1.0E3, vehicle1.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/crew-module-atmospheric-re-entry-experiment-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAF8CAYAAAA6gkcXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdd5gUddLA8W+RkaAgSZAgwYTKoYuiiICgYk4Yz3yy6usZznCceMYTPdOhhxFFMUdUlMOABA9UVIIJTyUnFcmgxN2t94/qcYd1F3Z3eqYn1Od5+umentnuslm35pdFVXHOOedKUyXqAJxzzqUvTxLOOefK5EnCOedcmTxJOOecK5MnCeecc2XyJOGcc65MniScc86VKe2ShIjcJiIabF8E52qLyBARuV1EHhaR+lHH6ZxzuaBa1AHEE5E6QCPgsODU/GD/EPCOqr4kIn2BR4EzIgjROedyiqTTiGsRuTI4fERVNwTnmgNzgQaquk5EqgGrgY6qOi+aSJ1zLjekTUlCRAQ4HegC3CIiF6vqC0BPYJmqrgNQ1QIRmQv0AObFX6NRo0bapk2bSt1/6dKlNG7cuNLxZxt/HsX8WRTzZ7GlbHkeU6dOXaaqpf6HpE2SUCvSdBWRBsDVwHMisgpoASwv8fG1QPOS16hTp84Wr/Pz88nPzy/X/fPy8pgyZUolIs9O/jyK+bMo5s9iS5n8PIYOHcrQoUNjL38t63NpkyRiVHUl8PegZHElMAZYX+JjNYHNJX+2cePGGfsP5pxzqRT/JVpElpX1ubTr3RTn30ArYDHQoMR79YAfUh6Rc87lmHROEkXANGA80FxEagOISHUseXwQ5s3KWy2VK/x5FPNnUcyfxZZy4XmkTe+moBdTT+BFQIE7sF5O80TkReBZVR0lIkcCZ6vqmSWvkZeXp17d5JxzFSMiU1U1r7T30qlNogHwT2AgMAm4L66L6yXAXSLSCWgGZH/6ds65NJA2SUJVZ2DVSKW9txLon9qInHPOpXObhHPOuYh5knDOOVcmTxLOOefK5EnCOedcmTxJOOecK1Pa9G5yziVAFQoLoaAAqlWDqlVBJOqoXBbwJOFculm7FhYtgsWLi7dFi+Dnn2HVquJt9WpYvx42bYLNJaYyE4EaNaB6dahdG3bYARo2hAYNbN+4MbRqVby1bAnNmnlicb/jScK5KBQVwZw58M038N138P33tv/uO0sGlVGlipUiCgtt27jRtl9+gaVLt/3zdevCnnva1rGjbfvtB02aVC4elxU8STiXbEVF8O23MG1a8TZ9OqxZU/rna9WCnXeGFi223Jo1s5LADjvYtv32sN12xSWGKlW2vOfmzVbKWLfOSh4rVhRvS5bAggXF2/z5dv7TT22L17YtHHAAdO0KBx4InTtbMnI5wf+lnQvb5s0wdSpMnGjbpEmwcuXvP7fTTrD33rDbbrbtuqvtd955yz/4lVGlCtSsaVu9etC06bZ/ZulSK9l88w3MmAFffWX/HXPm2PbCC/a57beHXr2gTx/o3dti9mqqrJU2E/yFwSf4c5FQhZkz4e23bZs40b69x2vRAvbfH/bd17bOnS1JpLuCAksYn3wCkyfDf/8Ls2dv+ZlWreD442075BAr1biMsrUJ/jxJOFcZmzfDuHHw1luWGObM2fL93XazP5jdu9vWunX2fNuePx/GjoX337ctvr1jhx3g6KOhXz848kgrybi050nCuTBs2mR/FF95BUaO3LIKqWFDOPxw+8N4+OHWfpALiorgs8/sebzxBvzvf8Xv7bADnHIKnHmmJcxEq9Bc0niScK6yVK1NYfhwGDHCup3GdOwIJ54IRx1lVUlVq0YWZtqYORNefx2efx6++KL4/M47w3nnwZ/+BG3aRBWdK4MnCecqasECePppSw7xdfD77GNVKf36wR57RBZeRpgxw5LF88/DvHl2TsRKWvn5cOyx3n6RJjxJOFceRUXw7rswZAi8846VIsAanc85x7bdd482xkykag3ejz0Gr75qYzfAquQuuQQuvtjHYkTMk4RzW7N6tZUYHnzQqkvAGlxPOAHOP9+6enpVUjhWrIBnnoGhQ62rLdiz/uMf4corrUuwSzlPEs6VZvFi+Ne/4NFH4ddf7VyrVvB//wcXXgg77hhtfNlMFSZMgMGDYdSo4lJbnz5w/fXQo0f29AbLAFtLEt7dwOWeWbOsTrxtW0sSv/5qg8Nee83aHwYM8ASRbCL2zN9806YiuewyqFPHeo/16gUHHwyjRxcnDxcZTxIud8ycadUau+1m9eObN1sXzWnTbMzDiSf6dBNR6NAB/v1vWLgQbr3VuhN/9JGNt9hvP+te68kiMp4kXPZbtMhKDnvsYT1tqlSBCy6wPv0vv2yjn130GjSAG26wwXr33GMN29OnW9tQ1642gM+lnCcJl72WL4err4b27a3kANbWMHs2DBtmJQqXfurWtX+3uXOthNG0qU062KePbZ98EnWEOcWThMs+mzdbN9YOHazNYeNGOPVU67f/2GPWOO3SX61a1lYxezYMGmQTC44da6WKU08tHnvhksqThMsuY8bAH/4Al19u02b07m1tDi+95CWHTFWnDgwcaCWLv/3NFlF65RUbs/L3v9t6GS5pPEm47LBwoTU8H3649b9v29bmEhozxtscskWDBnDHHbY2xxlnWAlx0CBL/s88443bSeJJwmW2wkKrWtpzT0sKdevCP/9pieL4472vfTZq1co6IHz4IXTpAj/8YKPhe/Wy7rQuVJ4kXOb68ks46CCrWvrlFzjpJPuWOWCAT1GdCw46yNa4GD7c1uz+4AObW+vmm4un/nAJ8yThMs/mzdaffr/9rNdLixZWihgxwo5d7qhSBc49174cXHihTed+yy3QqZPNF+US5knCZZbvvoNu3eCmm2zVtD//ubhqyeWuhg2t59oHH1iD9nff2dQef/kLrF8fdXQZzZOEywxFRdb28Ic/2CI3rVrZKOkhQ6B+/aijc+nikEPg88/tS0TVqnDffbZcrM/pVmmeJFz6W7YMjjnG2h42bLDFa7780hoqnSupZk1rl5g82UbZf/utja246SarqnQV4knCpbeJE6308PbbVqXw2mvw5JM2sMq5rcnLg6lTrcqpqMjasXr0sAWlXLl5knDpqagIbr/dSguLF1tPls8/t7EQzpVX7do26n7cOFtC9eOP7UvHW29FHVnG8CTh0s+yZXDkkbauQGGhdWmdMAFatow6Mpepeva0yQKPOspG4h93HFxzjVc/lYMnCZdePv/cqgnee8/WdBg92gbH+VrILlGNGlkJ4q67rFH73nut+unHH6OOLK15knDp4+WXrXvr/Pk2kvbzz61E4VxYqlSBa6+1MRQtW1r1U16e9ZhzpfIk4aJXWGgTuJ12GqxbZ1Ms/Pe/VofsXDIcdJB1i+3e3ab1OOQQm+rD/Y4nCRetX36xRWXuuMOqAAYPtmkWatWKOjKX7Zo0seVS+/e3rtV//KPNMltYGHVkacWThItO7BvcqFHWvfXdd+HKK31SPpc6NWrAo4/CAw/Yl5Q774TTT7ek4QBPEi4qX39tA5ymT7eV4z75xNZ+cC7VRODSS62zxPbbw6uvwhFHwKpVUUeWFjxJuNQbO9YaqBcuhAMPtMbD9u2jjsrlukMPtcGbzZtbm1j37rY+eo7zJOFS66mnoG9fWLMG+vWzhNGoUdRROWf23tu+tOyxh5V2DzzQJpDMYZ4kXOrceafNu1RQYAOZXnrJRsQ6l05atYJJk6y0u2iRjaX48suoo4qMJwmXfKpw3XXWc0TEZm69+27rs+5cOmrY0Ja+PfJImwGgVy9rP8tB/n+pS67CQrjkEhs1XbUqPPusrQHhXLqrXRtef91mIF6xwtoscnDKcU8SLnk2bbK+548+auMe3ngDzjwz6qicK7+aNW3FwxNOsN5OffpYT7wc4knCJceGDTZj60svQb168M479o3MuUxTo4ZNGXPyybB6tXW8yKE2Ck8SLnzr19tyoqNH2yR948db459zmap6dXjxRfvis2oVHH44zJoVdVQp4UnChSuWIN57Dxo3tim+99sv6qicS1y1aja/U+/esGQJHHaYrXWS5TxJuPDEEsSYMTYvzvjxsNdeUUflXHhq1bLG7P33h3nzrESxfHnUUSWVJwkXjvXrbSGX+ATRsWPUUTkXvnr1rCq1Y0cbaHfMMfb7n6XSMkmISA0R+UJEegava4vIEBG5XUQeFpH6EYfo4m3caHW1778PTZtagthzz6ijci55dtzRqlRbtYLJk+H8823J3SyUlkkCuBZoE/f6IWCSqg4ERgKPRhGUK0VBgXVzffdda4MYN84ThMsNzZvbDMb16lkvvptuijqipEi7JCEiBwE/AiuD182BM4HYyuXvA8eJSJso4nNxioogP9/6kW+/vSUKTxAul+y9t3WPrVIFbrsNnn466ohCl1ZJQkTqAP1U9Ym40z2BZaq6DkBVC4C5gPepjJIqXHUVPPmkjUz9z3+gc+eoo3Iu9fr2talmAC680GaSzSJplSSAAcA/S5xrAZTsPrAWaF7yh5cuXUpeXt5v29ChQ5MUpuPWW+H++63/+Ouv22RozuWq//s/uOIK2LwZTj0Vfvwx6oi2aejQob/9rQTKnIpZVDV1UW2FiBwB1FTVN4PX84DzgDzgFFU9IO6z04DnVfWe+Gvk5eXplBycWyXlHn/clnysUqV4JKpzua6gwMZOTJhgX5rGj7cvURlARKaqal5p76VTSeIa4AkRWSYiy4CWWCN1HaBBic/WA35IcXwObHqNiy+244ce8gThXEy1ajYqu0UL+PBD+Otfo44oFOmUJM4E/hC3/QBcCDwJNBeR2gAiUh1oBXwQUZy5a/p0OOUUm9n1uuvgoouijsi59NK0KbzyipUg7rvPkkaGS5skoapLVXVRbAMKgaWqugAYBcQWQO4DjFDV7B8Pn07mz4ejjoJffrGZXAcNijoi59LTgQdaggBryP7++2jjSVDaJIltuAQ4XkSuB44C8iOOJ7esWQNHHw0//QQ9e8ITT9jiQc650l1yCZxxBvz6q32p2rQp6ogqrVxJQkR2EJEqca+biMjByQsLVLWNqk4Ijleqan9VHaSql6nqL8m8t4tTVGSD5WbMsHV/X3vN5th3zpVNBB5+GNq0galTM3qg3TaThIgMAGYC40RkoIhUUdWfAe9fmgtuuMFGlTZoAG+9ZXvn3LZtv72txFiliq3vPmFC1BFVSnlKEhcDu6tqT+Bp4AoR2R74NZmBuTTw0ktw++227OjLL0O7dlFH5Fxm6dYN/v53G3x69tm2DGqGKU+S+BJYDRA0Kg8GjgZ8kr1sNm2aTVoGcO+9tmyjc67ibrgBunaFRYtswF2GKU+S+BPQK/6Eqj4PPJOUiFz0li2zNX3Xr7dEcfnlUUfkXOaqVg2eecamr3n2WZtmPINsM0mo6jJVHVPK+duSE5KLVFERnHUWLFxo334efth7MjmXqPbt4R//sOOLLrIegxkiaV1gReSoZF3bJdEdd9hsrjvuaO0Q3pPJuXBccQV06WLVTn/7W9TRlFvCSUJE+orIJBGZJSJzgm0eMCLx8FxKjR8PN95oJYdnn4WWLaOOyLnsUa2ajTGqXt1K6B9kxqQRYZQkngFeA/oD5wfbeXibRWb58Ucb/FNUBNdfb9MfO+fCtddeMHCgHV98cUYMsqsWwjXmqOq/Sp4UkW9DuLZLhYICSxBLlkCvXnDzzVFH5Fz2uu46eP55+PZbW4fi6qujjmirwihJXCIi/Uo5f1wI13ap8I9/WNG3WTP75a1aNeqInMteNWvaWiwAt9yS9mtPhJEkugGPi0hh3FYEPBzCtV2yffihLbsoYgmiWbOoI3Iu+x15JBx7LKxdm/aN2GEkiVuxCfc6AG3jtptDuLZLpjVrrLtrUZHNfd+r17Z/xjkXjsGDrVTx9NPw0UdRR1OmMJLE28BIVZ2jqvODbR7wWAjXdsn05z/DvHmw7762HKlzLnXatYNrrrHjK6+0qTvSUBhJ4kXgOhE5JG7rAVwXwrVdsrz4YvEo0Oeegxo1oo7Iudxz3XW2UNFnn9la8Wko4TWuReQzrKppFRC7mADNVLVWYuFVjK9xXU4LFsA++8Dq1fDII77CnHNReughuPRS2H13+OorG0+RYsle4/puoGWw/sMuwdYGODGEa7uwqcIFF1iCOO44yPf1m5yL1IUXQtu21iX2qaeijuZ3wkgSHwF/F5G6ACJykIicoKpvh3BtF7bHHoOxY6FRIzv2eZmci1aNGsXzOt18s02smUbCSBLPAgcA1QFU9SOghYgMDOHaLkzz5xcP3HnwQWjSJNp4nHPm9NOhUyeb1+mRR6KOZgthJInPVLWnqq6MOzcV+EsI13ZhUYX+/eGXX+Dkk+GUU6KOyDkXU6VKcWninntg48Zo44kTRpIoEJHfpgoVkdpYz6YlIVzbhWXYMBgzxmZ3ffBBr2ZyLt0cc4x1KPnhh7RqmwgjSQwHJorIKyIyElgEdMeWPXXpYMECuOoqO37gAety55xLLyLFk//deafNqZYGEk4Sqvod0AN4FfgQuBxoq6qTEr22C4GqdXFduxZOPBFOOy3qiJxzZenXzxYomjPH1phPAxVOEiJyt4hcJSK7xs6p6npVfUlV71LV51R1Vbhhukp76SV45x1o0MD6Y3s1k3Ppq2rV4rmc7rjDpsyJWIWThKpei/Vo6h4kjJtFpJuI//VJO6tW2XB/sOKrT97nXPo7+2xo0QJmzID33486mspVN6nqz6o6LEgY/wQaAneIyB0icqKIbBdqlK5yBg60NSK6dYM//SnqaJxz5VGjBlxyiR0PGRJtLFRyWg4Rqamqv+ujFZQmDgAOB+qq6l8TD7H8fFqOOJ98AgceaMXX6dNtRSznXGb4+WdbPnjzZpg1y0ZkJ1EypuW4rbSTaiar6q2pThAuTkGBNVar2uA5TxDOZZYmTWyAnap1WY9QZZNEXxF5QkRuE5H9Q43IJe7+++GLL6BNG7jxxqijcc5VxmWX2f6JJ+DXXyMLo7JJ4gJVvQCb3G83EXlARG4XkW4hxuYq44cf4Kab7PjBB2E7bx5yLiPl5UHXrtYB5cUXIwujsg3XnwX71ar6jKr+GRgEdBCRaSIyJFhTwqXa3/5m3zpOOAGOOirqaJxziYjN0jx8eGQhVCpJiEj9uOOaInICthLdEGB3oDngs8el2uTJtpBQzZpw771RR+OcS1S/flYbMGmSNWBHoNIN1yLyRxF5AVgKvADUAS4Cmqjqyar6SlhBunIoKoLLL7fjq69Oem8I51wK1KtniQJsLewIVDZJ/Bl4HKgPfAsMBq5Q1edV9ZewgnMV8NRTtgRi8+a2JKJzLjucd57tn3oqkhHYlV0n73Wgv6quABCRtsApItIGmAO8qqrzQ4nQbduaNcWJ4a67oG7daONxzoWnRw9o3drWg5kwAQ49NKW3r2xJ4opYggBQ1TnA+8Aa4FpgtojcF0J8rjwGDbKR1QcdBGeeGXU0zrkwValiU3VAJJP+VTZJtAAQkc5B19eZwBTgMOBeoJ2qXhlSjG5rFiywcRFge59Cy7nsE1sk7PXXobAwpbeudHWTiGwGWgIzsDUlXlLVaJrfc9kNN9gqVmecYf2qnXPZZ++9bQrxWbNg4kTo2TNlt65sSaI6lhj2UtW9VXWQJ4gIfP65dXmtXt2qnJxz2UnElh0GGDEipbeubJI4R1VvUtVvQo3GVcyAATa3y6WXwi67RB2Ncy6ZYknitddS2supsiOu3w47EFdBY8bAe+9B/fpw/fVRR+OcS7a8PGjVyqbe+fTTlN02jDWuXaoVFVkpAqzra6NG0cbjnEs+ETj2WDt+552U3daTRCZ68UVbI6JFC7jiiqijcc6lyhFH2P7dd1N2S08SmaagoHiW11tugdq1o43HOZc6vXpZR5VPP4WVK1Nyy6QkiWD8xLMicnUwCtuF5dlnrRtc+/Zw7rlRR+OcS6W6dW3QbFFRyta/TlZJ4mKgj6reCxwkIv51NwybN8Ott9rxTTdBtcoOc3HOZawUVzklK0n8GdgvOH4RG1fhEjV8OMydC7vtZoPnnHO55/DDbT9uXEpul3CSEJF6InKniNwZvN4fOF1VFwOoapGqrkn0Pjlv40a4LVha/OaboWrVSMNxzkWkUyerdpo7F378Mem3C6MkMRzoBjQAUNVPgZoicncI13Yxw4bZPE0dO8Kpp0YdjXMuKtWq2bKmAB99lPTbhZEkNqvqwdgcTjELgAtCuLYD2LABbr/djm+5xWaFdM7lrm7dbP/hh0m/VRh/bf4X7BVsOVPgr8DiEK7twBYbWbwY9tkHTjwx6micc1E76CDbp6AkEUb3mFHBMqYNRWRP4BhsKdPjQri2KyiwhYTARld7KcI517WrjcCeNg3Wr0/qeKmE/+Ko6lQgH3gSW5XuOqC9qk5M9NoOePllmDPHxkXE5pR3zuW2+vWtfXLzZvjqq6TeKozeTeeq6lpVfVFV71LVZ4BGInJsJa7VQ0S+E5GV8SvbiUhtERkSLHD0sIjUTzTujFBUBHfcYccDBniPJudcsU6dbP/FF0m9TaWThIh0F5FDgMNE5JD4DdgDeLCC19se6yXVFTgLuFRE+gRvPwRMUtWBwEjg0crGnVH+8x/4+mto3rx4+ULnnANrowT48suk3iaRNolVwNPALtgf93ibgMcreL1NwB2qqsB/RORLoFBEmgNnApcGn3sfGCEibVR1XmWDT3uqxT2arrkGataMNh7nXHpJUUmi0klCVb8SkR7Aqaq6RUIQkR1VdXkFr7c+7ufrYL2mJgBnAMtUdV3wuQIRmQv0AObFX2Pp0qXkxS3hmZ+fT35+fkXCSB8ffACTJ0PDhtC/f9TROOfSTXxJQrXC69sPHTqUoUOHxl6Wud5AQr2bgpHUj4tIA6AuEIuyiYjcraoVHishIt2BQcCPwHZAC6BkwlkLNC/5s40bN2bKlCkVvWV6uvNO219xhY2udM65eM2aQePGsHSpDbRt3bpCPx7/JVpElpX1uTAaru/G/ojPi9s+Afat5CVnYaO4ewN3Y+Mv1pf4TE1gcyWvn/6++cYWFald25Ymdc65kkRg993teNaspN0mjE73LYEmQD+glapWAY7AusJWmKr+qKpPANdgVUqLCab8iFMP+KHSEae7+4KOXeedBzvuGGkozrk01q6d7WfPTtotwkgSX6vqMuAt4Jzg3H+BwQledwqWIMYDzWPTjYtIdaAV8EGC109PS5fC00/bsa8655zbmrZtbT9nTtJuEUaSqC4is4AuwEwRGQ98XNFri8h2ItI57tRRwP2q+hMwCqt+AugDjIjNMpt1HnnEZnw9+mibEtw558qSgiSR8LQcqnqTiLwMzFTVj4MGkM7AqxW81F5Y19dvgcnAp6r6n+C9S4C7RKQT0Awb4Z19Nm6EB4PhJVddFW0szrn0F6tuSuckISK3A++r6gwAVR2PVRFVSDDFeOMy3lsJZH8/0BdegCVLrGtbr15RR+OcS3e77GL7NK9u6gf8blEhEdkphGvnDlUYHDTj/OUvFe7z7JzLQY0b23Q9K1fCpk1JuUUYs8AOAs4MBsBpcK4qcCHwxxCunxsmTbJBMU2a+NKkzrnyqVLFEsVPP8HPP8POO4d+izCSxP8BuwInUZwkwNoOXHk99JDt8/N9Cg7nXPk1bZr2SeIuYFzQbvAbETk5hGvnhiVLYMQI+1aQqdOIOOei0aSJ7X/+OSmXD2M9iRGxBCEiLUWkaux8otfOGcOG2bzwxxwDLVtGHY1zLpM0bWr7JUuScvmwlzn7DmgX8jWzW2EhPBrMfP5//xdtLM65zNMomJtvWZnTLyUk7CThXXIqavRom5yrXTs47LCoo3HOZZrYBKC//pqUy/uCyVF7+GHbX3yxr1/tnKu47baz/bp1Sbl8KMuXxr3sD/woIrtVZvnSnDNnjs32WrMmnH9+1NE45zJRnTq2T1KSqHTvpmDdB8GWL50bnF6ATcnREPg3NumfK8uwYTaI7tRTfbZX51zlxEoSSapuSqflS3NLYSEMH27HvvKcc66yklzdlJTlS105vPce/PADtG8PBx8cdTTOuUwVG3ybpGk5EmqTUNU1pSUIEWmYyHVzwpNP2v78832eJudc5RUV2T5Jf0dC6U4jIg2CgXStRKQV0FZEngjj2llp2TJ44w3rzXTOOdv+vHPOlUWD2ZCSlCTCmCr8buBqbN6mWJQKfJXotbPW88/bCOsjj0zKXCvOuRwSSxJJ6kIfxtxNsTWuuwOfqeoiEekDVA/h2tlH1Xo1gXd7dc4lLgOqm5K1xnV2mj7dpgRv2BCOOy7qaJxzmS7dq5soXuP6bIrXuK6Pj+YuXazB+qyzfEpw51zi1q+3fa1aSbl8mGtcf5/gGtfZb/NmePFFOz733K1/1jnnymPtWtvXr5+Uy4fRcF2P4mqmAcCvwFJVXZDotbPOmDHWs2mPPaBz56ijcc5lg1iSqFcvKZcPo0poODbiugGAqn4K1Ax6Pbl4zz1n+zPP9LERzrlwZECS2KyqBwMz4s4tAC4I4drZ49dfbWwEWJJwzrkwrFlj+zROEv8L9gogIjWBvwKLQ7h29hg50uZWOfBAaNs26micc9kiyUkijN5No0TkBaChiOwJHAPUAbx/Z7z4qibnnAtLbNnS2DKmIQujd9NUEckHjgZaAdcBo1V1eaLXzhpLl8K770LVqjYtuHPOheXHH22/005JuXwYJQlUdS3wYvw5ETlXVZ8K4/oZ75VXbGrwI4+EJk2ijsY5l02SnCTCWJnuTyIyX0Q2ikhhsBUBPsFfzEsv2f6MM6KNwzmXXdavh1WroHr1pC1cFkZJYjCQD3wGFATnBDg9hGtnvp9+gokToUYNOP74qKNxzmWTWHtEs2ZpPS3HGOANVd0Qf1JEngvh2pnv9ddtbpUjjkjaiEjnXI5aEIxZbtEiabeocJIQkUNKnHoa+IeIlFzPuh9weWUDyxojRti+X79o43DOZZ+ZM23fvn3SblGZksQ/gU7AMqAo7nz8X0EBmpLrSWLpUpgwAapVg2OPjToa51y2iSWJDh2SdovKJIkbgIWq+v3WPiQiPuXX9pkAACAASURBVE5i5Ejr1dS3LzRoEHU0zrlsM2uW7dOpJKGqY8v5uTcrHk6WeTWYCNermpxzyZCCkoSv+ZAsK1bA2LE2gM57NTnnwlZUlJKShCeJZHnzTSgogF69oFGjqKNxzmWb2bNtPriddkpqdXYYg+l+t3qOiOwmIrndUjtypO1POinaOJxz2Wn6dNsneW2aSo+TEJHuWC+mw0Rkbom3GwL/xta9zj0bNsB779mx92pyziXD55/bPl2TBLAKGyOxC7boULxNwOMJXDuzjRtnxcDOnWHnnaOOxjmXjdK9JKGqX4lID+BUVc3dhFCat4IC1HHeC9g5lyQpShIJtUmo6pqyEoSInJfItTOWanGS8Kom51wyLFxo8zZtvz20aZPUW4XRcN1fRBaJyOYSs8AOCyG+zDN9OixeDM2bw777Rh2Ncy4bTZxo+27doEpyO6mGMcHfPcBF+CywJr4UkaRZGZ1zOS6WJA4+OOm3CiNJjAVGqeov8SdF5JkQrp153gwGmntVk3MuWWJJonv3pN8qjCTxNHC7iLwad06wCf8uC+H6mWPxYpg2DWrXhkMPjToa51w2WrECZsyAmjWhS5ek3y6MJDEQaAccB2hwrgrQjFxLEu++a/vevS1ROOdc2CZNsv3++1uiSLIwksQgYLSqbo4/KSJ9Q7h2Zoklib6595/unEuRd96xfe/eKbldwklCVUeW8db8RK+dUQoLYcwYOz7iiGhjcc5lJ1V4+207PvLIlNyyMivTDQFeUNWPgtcPAyXLPALsD3RMOMJMMWUKrFwJbdsmdUZG51wO++47mDfPJg3Ny0vJLStTkliATckRUwvYCPxEcZuEYO0UuSNW1eSlCOdcssRKEUcckfTxETGVWXTo7hKnBgELVHVT/EkReTKRwDKOJwnnXLKluKoJwmmTmCUidUTkNGBn4DvgP6qaO20Sq1bBJ5/YWta9ekUdjXMuG61YAePH20JmKfwyGsa0HF2AucBg4GTgemCaiOxZiWudJCJzRWS5iNwvItWC841E5FERuVNE/iUi1RONO1Rjx1rD9UEHQf36UUfjnMtGb7xhC5kdemhKFzILo1LrTuASoLGq5qnqftjU4RdX5CIi0go4ARuE9xfgQuDK4O1XgCdUdQDwPXBrCHGHx6uanHPJ9vLLtj/llJTeNoxxEpNVdUT8CVVdJSIbK3id1sAFqloATBWRvYBeIjIJ2E1VPwk+9xbwnYjcpqq/Jhx9GGJdXw8/PNo4nHPZaflyq7GoWhVOPDGltw6jJFFXROrFXohIdRE5G+hakYuo6sQgQcQsxnpSHYpVZ8U+tzg4/F3/r6VLl5KXl/fbNnTo0IqEUDnz5tm2ww5Jn9fdOZejYlVNvXuHVtU0dOjQ3/5WAmVeNIySxEPAxyKyDqgOtAV+ARKd4a4LcDeQDywv8d5aoHnJH2jcuDFTpkxJ8LYVNH687Xv0sCzvnHNhe/ZZ2596amiXzM/PJz8/HwARWVbW58Lo3fStiOwHHA10wL79v6Gq6yt7TRFpDyxV1S9ERIGS16oJbP79T0ZgwgTb9+wZZRTOuWw1e7b9naldO+XtERBOdRNAE2AH7A/3WhL4Ay4iVbH1Kf4WnFoMNIh7X4A6wA+VvUdoVItLEt711TmXDE8GQ85OOSWS3pNhdIE9F+txNBg4G+vt9KWIVHZKjmuAe1Q11vA9Btgt7v02wCYgxfVKpZgzx5YR3HFH2HvvqKNxzmWbwkIYPtyOL7ggkhDC6gJ7D9BEVTurakfgMODPFb2QiAwEpgO1RaStiFyATQHyY9DbCaAvMKTkCO9IxLdHpGiIvHMuh4wZY+vUtG8PhxwSSQhhNFwvAe6M++aPqi4Wkf/FXovIdqq6bmsXEZG/A/8ocfpbVX0iGM19o4jMA+piA/ai5+0Rzrlkevhh259/fmTLIYuqbvtTW7uASB+sO+rzcad3Av6I9U6qApynqrckdKNyyMvL05T1blKFnXeGH36Ar76Cvfba9s8451x5zZ4NHTpA9eqwYAE0bZq0W4nIVFUtdVrZMEoS9wF7AreX8l6sykmBpCeJlJo1yxJE48bQMXdmRHfOpcgDD9iX0TPOSGqC2JYwKtL/DeykqlVK24CqwFUh3Ce9xC9EHlEx0DmXpdasgWHD7PiKKyINJYxxElsd1qyqKiIPJXqftPPhh7bv1i3aOJxz2Wf4cFi71hqrI57JISVdckquf50VPvrI9p4knHNh2rwZBg+244hLEZCiJJF1li+Hb7+FWrUiz/LOuSzz3HM2H9xuu8Hxx0cdjSeJSomVIrp0gRo1oo3FOZc9Cgvh9qAP0PXXp8V8cElLEiJyXrKuHTmvanLOJcPLL8PMmbDLLtarKQ2EMS1HfxFZJCKbRaQw2IqAYSHEl55ijdYHHRRtHM657FFUBIMG2fF119lyyGkgjCjuwSbk+wyIrQchwOkhXDv9bNoEn31mx54knHNhefFFmDEDWraEc8+NOprfhJEkxgKjVPWX+JMi8kwI104/06fDhg2w++42sZ9zziVq40ZrgwC4+ea0ausMI0k8DdwuIq/GnRNsrerLQrh+eom1R3gpwjkXlkcftR5Ne+4J55wTdTRbCCNJDATaAcdh02+AtXU0IxuTxKef2r5rhVZndc650q1ZA/8I5ja94460aYuICSOaQcDokgPmRKRvCNdOP7EJBPNKnQvLOecq5u67Ydky6y15bKKrPocvjGk5RopIHRE5HdgZ+BZLGu8kHF26WbnSJvarVctnfXXOJW7uXEsSAHfemZbzwCWcJESkC/AfrIppHtYecZuInKKq3yR6/bQSK0X84Q82fa9zziXiqqus0frMM9N23FVYK9NdAjRW1TxV3Q/oBlwcwrXTSyxJdOkSbRzOucz37rvwxhtQt25xaSINhdEmMVlVR8SfUNVVIrKxrB/IWLHxEd4e4ZxLxKZNcPnldnzDDdC8ebTxbEUYJYm6IlIv9kJEqovI2UD2df/xkoRzLgz33APffw+77gpXXhl1NFsVRkniIeBjEVkHVAfaAr8A6ddMn4glS2DhQqhXz2ZndM65yvj2W7j1Vjt+8MG0GjhXmjB6N30rIvsBRwMdgAXAG6q6PtFrp5VYVdN++0EVnzzXOVcJRUXQv781Vp9/PvTpE3VE2xTKqA1V3Qi8Fn9ORPZQ1f+Fcf204OMjnHOJeuQRmDQJmjWDe++NOppyqXCSEJEhwAuq+lHw+mGgZsmPAfsDHROOMF18/rnt99032jicc5lpwQIYMMCOH3gAGjSINp5yqkxJYgGwKu51LWAj8BPF03IINlVH9vjiC9t36hRtHM65zFNYaHMy/fILnHginHxy1BGVW4WThKqW7ND7IfBU/LQcIrIbsDzB2NLH6tU2+VbNmtYbwTnnKuJf/4IPPoCmTW0yvwxS6TYJEemOlRh6At+JbDGevAFwLTAkoejSxVdf2b5jx7SbfMs5l+a++KJ4GvBhw6Bx42jjqaBE/uKtwqYJ3wUbYR1vE/B4AtdOL17V5JyrjA0b4I9/hM2b4ZJL4Oijo46owiqdJFT1KxHpAZyqqtmTEErjScI5VxlXXmmrze26qw2gy0AJdfhX1TWlJQgRaZjIddOOJwnnXEU995y1P9SsCS+8ANttF3VElRLKqDARaSAiLUWklYi0AtqKyBNhXDtyhYXFbRL77BNtLM65zPDtt3DRRXZ8//0Z3XU+jKnC7wauxrq/xhqvFfgq0WunhVmzYP162HlnaJhdBSTnXBKsWwf9+sGvv9oU4Pn5UUeUkDC66rQEmgDdgc9UdZGI9MHmccp8XtXknCsvVbjwQmuH2H13q25Kw4WEKiKM6qavVXUZ8BYQW8H7v8DgEK4dvRkzbL/33tHG4ZxLf3feae0PdevCq6/aPsOFUZKoLiKzgLOBmSIyHqhPSO0dkftfMP3UHntEG4dzLr2NGgUDB9rxs8/auKosEMYssDeJyMvATFX9WESWA38AXk04unQQSxJ77hltHM659PW//1n7gyr84x9w/PFRRxSaUL7tq+oMVd0UHI9T1X8Bh4Zx7UgVFNjCIGD1i845V9LPP8Mxx8DatXDKKcWjq7NEwklCRPqLyCIR2SwihbENGBZCfNGaO9eWGWzZMivqFp1zIVu3Do49FubMsbVmnnwy4xuqSwqjTeIe4CLgM6AgOCfA6SFcO1rffGN7b49wzpVUWGhVTJ9+Cq1bW5tEnTpRRxW6MJLEWGCUqv4Sf1JEngnh2tHyRmvnXGlUbcqNkSNhhx3g7bdtIaEsFEaSeBq4XUTiG6oF6AdcFsL1o+ON1s650gwaZAsH1ahhiSKLv0iGkSQGYgsMHUfxokNVgGZkS5LI4l8A51wF3X8/3HCDrXX/zDNwyCFRR5RUYSSJQcDo+EWHAESkbwjXjo6qzb8CniScc+aJJ6yaCeCxx+DUU6ONJwXC6AK7S8kEESgQkVtF5IAQ7pF6P/1kXdoaNoRGjaKOxjkXtZdfhv797fi+++CCC6KNJ0XCSBIXi8gkEZkqIn8FEJHawBvAo8BuItIrhPuk1qxZtu/QIdo4nHPRGzXKFg8qKrLBcldcEXVEKRNGkmgJfA48ATQRkbOBFsB2qroYeAZrt8gss2fbvn37aONwzkVr5Eg46SQbXHvttVk3WG5bQmm4VtX7Yy9E5EygBraEKaqqIrJrCPdJrVhJol27aONwzkXnlVdsLERBAfzlLzaBX5YNltuWMJJEKxHpDGwAegNdgEXAarAFiYCaIdwntWJJwksSzuWm55+Hs8+2KqYBA+COO3IuQUA4SeJ5rP2hJTbq+s9Y1dNkEfk7sDswJoT7pJZXNzmXu556Cs4/33o53ngj3HxzTiYICGcW2KlAaxFprKpLAUSkp6ouF5HGWNJ4ONH7pJQqzJxpx17d5FxuGTLEGqZjM7r+/e9RRxSp0Na4BmrFrXG9i4g8oapLVfUmVf0pjPukzIoVsHo11KsHjRtHHY1zLhVUbT2Iyy+347vvzvkEAb7Gdeniq5pytIjpXE4pKICLLrLBclWrwrBhcO65UUeVFnyN69J4zybncse6dXDaaTYWonZt69F09NFRR5U20m6NaxE5TEQ+EZE2cedqi8gQEbldRB4WkfqJBr1VsZKEJwnnstvPP0OfPpYgGjaEceM8QZQQRpKIrXHdheI1rj+uzLVFpClQF9i/xFsPAZNUdSAwEhvJnTzz59t+l12SehvnXIS+/hr23x8+/tgWFps0Cbp2jTqqtJOMNa6XAZ2pxBrXqrpEREbGnxOR5sCZwKXBqfeBESLSRlXnJRZ9GRYssH2rVkm5vHMuYqNHw+mn2/xs++9vo6qzdD2IRIWxfOm5Jda4Hg/8B+hUmeupalGJUz2BZaq6Lni/AJgL9Kh00NviScK57KRqU30fe6wliNNOgwkTPEFsRaVLEiLSHevNdJiIzC3xdgNgCNZOkagWwPIS59YCzUt+cOnSpeTl5f32Oj8/n/z8/IrdTbU4SbRsWbGfdc6lr40b4bLLbIpvgJtusi1HezAOHTqUoUOHxl6WOdV1ItVNq7BV6XYBupV4bxPweALXjqfA+hLnagK/m568cePGTJkyJbG7LV8O69fbkoT1k9s+7pxLkYUL4eST4bPPoGZNGD7cqptyWPyX6KCZoFSVThKq+pWI9ABOVdWwEkJpFmMlk3j1gB+ScjevanIuu4wbZ9VKy5ZB69bw2muw775RR5UxKtUmISJnBrO9Hg2MD87tLSIfi8hKEXleRLYPKcbxQPNgjQpEpDrQCvggpOtvyZOEc9lBFe66Cw47zBLE4YfD1KmeICqosg3XzwKHYd1SZwdzNI3FBtZdgY22vq8yFxb5rYJQAIIpPUZhM8wC9AFGBGtVhC/W/dWThHOZa/lyOPFEm721qMjWgBg9GnbcMerIMk5lq5u+VNXz417/E9gB2E9VvwIQkScrelERqQecFbw8T0QeCCYNvAS4S0Q6Ac2ACrZGV4CXJJzLbBMn2hoQixbB9tvbjK7HHx91VBmrskni09hBsIb1ecCQWIIIVK3oRVV1LTZj7MMlzq8E+lcq0oryJOFcZioshNtvt2m9i4psYNwLL0CbNlFHltEqmyR2FJG9gHVY1dMC4LfpEkVkJ6y9IvMsXGh77/7qXOZYvBjOOsvGPIjA3/4Gt94K1TN3Crl0Udk2ieuxEdWzgALgaFX9RURai8g1wGSs+inz/BTMar7TTtHG4ZzbNlV47jnYay9LEE2bwrvv2ipyniBCUamShKp+C+wuIg2CqqCYhdgguiFhBJdyqvDjj3bsScK59LZ0KVxyCYwYYa+PPtqm+G7aNNq4skxCczeVSBCxKTU2JhRRlFauhE2bbBDddttFHY1zrixvvgn9+9ssrnXrwn33wQUX5Ozo6WQKYz2J7BGravJ5XJxLT6tWwV/+YiOmAXr0sGNvnE6aUJYvzRpe1eRcelK1aqU997SkUKsWDB5so6k9QSSVlyTieZJwLv0sWgSXXmpVTAAHHmhtD3vsEW1cOcJLEvG8usm59FFYCA88YKWHN9+EevXgoYdscSBPECnjJYl4XpJwLj188QVcfDFMnmyvTzjBEkaLFtHGlYO8JBHPk4Rz0VqxwqqW9t3XEkTz5jZr6+uve4KIiCeJeF7d5Fw0Cgvh0Udh112tSkkErrgCvvnGJupzkfHqpng//2z7Jk2ijcO5XPLRR7Zi3LRp9rpXL/j3v20UtYuclyTiLQ9WSW1U5kp+zrmwzJkDZ5wB3bpZgmjZEl5+GcaO9QSRRrwkEaNanCR8znnnkmfFCrjtNmuI3rzZlhO95hq47jqoUyfq6FwJniRi1q61X9g6dWygjnMuXBs2wIMPWoJYtcraHc4+21771Pxpy5NEjJcinEuOggJ49lm45RaYN8/O9e4Nd98NnTtHGprbNk8SMZ4knAtXYSG89JIlh++/t3N77WXrTvft65PxZQhPEjHLltneG62dS0xRkY1tuOkm68IK0K6dvT7zTKha4UUrXYQ8ScR4ScK5xKjCW2/BjTfaiGmA1q3hhhvgnHN8EaAM5UkixpOEc5VTUACvvGKrwX0VLHPfogVcfz386U9Qo0a08bmEeJKI8eom5ypmwwZ46ilrY5gzx841bw5//StcdJH3EswSniRivCThXPmsXQuPPAL/+lfxVDbt28OAAdaltWbNaONzofIkEbNmje132CHaOJxLV/Pn2wC4xx6D1avtXKdONgiuXz9vkM5SniRiYkmiXr1o43AunajCxx/bGtKvvWbdWgEOPhgGDvSurDnAk0RMLEnUrx9tHM6lg82bbbnQwYPh00/tXLVq1oX1yiuhS5do43Mp40kixpOEc7BggVUnPf54cXtDw4bWEH3ppb6mQw7yJBHjScLlqsJCeOcda4wePdoGw4EtEXrFFdYYvd120cboIuNJIsaThMs1P/0ETzwBQ4daozTYgLfTTrOlQ7t39/YG50niN2vX2t6ThMtmmzbBqFEwfLiVGmIN0W3bWpXSeef5oltuC54kwBrp1q+3Lny1a0cdjXPhUoXPP4cnn4Tnny8eE1S1KpxwgpUaDjsMqvgaZO73PEnAlqUIL167bPHDDzYL6/Dh8OWXxef33hvOP996KjVtGll4LjN4kgBvj3DZY/ly67r64oswYYKVIsB6KP3xj1ad1Lmzfxly5eZJAopLEnXrRhuHc5Wxdi2MHGmJ4d13bcI9sIn1jjzSZmA9+mifLsNViicJsPYI8G5+LnOsXm0Nz6+9Zg3RGzbY+apV4fDD4YwzrL3Bp5lxCfIkAcVJwmetdOlsyRIrMbz+Oowdax0uYg4+2BJDv37eO8mFypMEFH8L855NLt3MnWtJ4bXX4KOPitsYqlSBQw6BE0+Ek06CVq2ijdNlLU8S4CUJlz42b4YPP7SqpNGjYcaM4vdq1LCuqieeCMcdB40bRxenyxmeJMBLEi5aP/xg02KMHg3vvVfckQJsVuKjjrLSwpFH+izFLuU8SYCXJFxqbdxoVUfvvw9vvw3Tp2/5/p57WmI46ijo1s2X/3SR8iQBxUnCSxIuGQoLbcTz++9bg/OkScW/c2C96nr3tqRw5JHQunV0sTpXgicJKK5u8pKEC4MqzJxpCWHsWBg3Dlau3PIz++xjiaFvX2uA9t89l6Y8SYCXJFxiCgtt2ouJE+G//7X9zz9v+ZnWra3RuXdvOPRQ76bqMoYnCfCShKuYjRvhs88sGUycaL2RYlO7xDRpAj16QJ8+lhjatvWpMFxG8iQB9j89+LQF7vdUba2FTz6xbfJkmDat+HcmZpddbP2F7t2t+qhDB08KLit4koDiOfWrVo02Dhe9tWutlDB5cnFiWLLk95/ba6/ipNC9O+y8c+pjdS4FPElA8YRo1fxx5JQ1a+CLL6xkMH06TJkC33xTPKo5pmFD6NoVDjjAtv33hwYNoonZuRTzv4rgJYlcsGyZJYJp04qTwsyZv/9c9erwhz8UJ4SuXaFdO686cjnLkwR4ksgmGzbAt9/C11/bNmOGlRYWLvz9Z6tXtwV4OneGffe1fefO3oHBuTieJKA4SXh1U+YoKIBZs4qTQWybOROKin7/+Tp1oFOn4mSw7742stlHMzu3Vf5XEYrbJLwkkV5UrZro++/hu+9six3Png2bNv3+Z6pUgd12s4bl+K1DB//3da4SPEmAVzdFbf16KxXEJ4HYccmRyvFat/59Mth9d68uci5EGZMkRKQ2cBewFmgADFDVNVv/qXLyJJF8BQU23uD773+/LVhQ9s/Vq2clg912g113Ld536ODLzTqXAhmTJICHgHdU9SUR6Qs8CpwRypW9C2x4iorg44+tK2l8Ipg9e8uV1OJVq2Y9iDp0KE4IsWTQtKn3LHIuQhnxV1FEmgNnApcGp94HRohIG1Wdl/ANvCQRnltugVtvLd9nO3SAjh2hfXs49VTo0iW5sTnnKiwjkgTQE1imqusAVLVAROYCPYB5CV99772t7rtZs4QvlfOaNy//Z2fOLB6rsPvuniScS0OZkiRaAMtLnFsLbPEXaenSpeTl5f32Oj8/n/z8/G1f/cYbE4/QmYsugrPOsnaGlSvh119tW7fOqpsKCkrfx/27OeeSb+jQoQwdOjT2slFZnxMtOQVBGhKRa4BTVPWAuHPTgOdV9Z7Yuby8PJ0yZUoUITrnXMYSkamqWuo3tSqpDqaSFmM9muLVA36IIBbnnMsZmZIkxgPNg26wiEh1oBXwQaRROedclsuIJKGqPwGjgN7BqT7ACFVdHF1UzjmX/TKl4RrgEuAuEekENAPK0SLtnHMuERmTJFR1JdA/6jiccy6XZER1k3POuWh4knDOOVcmTxLOOefK5EnCOedcmTxJBOKGpzv8ecTzZ1HMn8WWcuF5eJII5MI/dkX48yjmz6KYP4st5cLz8CThnHOuTBkxwV95ichSYH4lf7wRsCzEcDKdP49i/iyK+bPYUrY8j9aq2ri0N7IqSTjnnAuXVzc555wrkycJ55xzZfIk4ZxzrkyeJJxzzpUpY2aBDUOwaNFd2PrYDYABqrqmlM8dCRyNJdFxqvpqSgNNgQo8i4bAtUB9Vb00tVGmTnmeh4hsDzwOHAHMBi5R1cmpjjXZyvksqgGPAqcA84BzVXV6ikNNifL+vxL3+dOw342eqYkwuXKtJPEQMElVBwIjsV/yLYjIrsAg4HJV/T/gchHZK7VhpsQ2n0WgHdARqJOqwCJSnucxAHgN6IUtnTtSRLLxuZTnWZwVnG8HLASGpC68lCvv/yuISHPg+lQFlgo50wU2+MebCzRQ1XXBN6HVQEdVnRf3uUeAFcEvBCJyPdBeVc+PIOykKO+ziPv8rUArVT0vpYGmSAV+Nw5V1XHB8fZY//iDVPWzCMJOigo8i+1UdV1wfCxwjar2iCLmZKrE/yv3A98Bp3pJIvP0BJbFfrFVtQD7xy/5i31ocD7mu+Bns0lPyvcsYopSFFdUelKO5xFLEMHxamANsCh1YaZET8r3LNbFvdwdGJiqAFOsJ+X8f0VE/gQ8D6wr+V4my6Uk0QJYXuLcWqD5Nj5X2mcyXXmfRa6o8PMIqiUnqOqPyQwsAuV+FiKyo4hcC/wZ2D4FsUWhXM9DRNoBDVX1k1QFliq5lCQUWF/iXE1g8zY+V9pnMl15n0WuqMzzuBRr0M82FXkWa4EJwGfA6yLSMrmhRWKbz0NEqmC/D4NTGFfK5FLvpsVYz4R49bAGyK19rrTPZLryPotcUaHnISJ9sFLEnGQHFoFyPwtV3QR8FvTm+R7oijViZ5PyPI9uQD5wjoiAJZFaIrJMVRulJMokyqWSxHigedCdDRGpDrQCPijxuTHAbnGvdw1+NpuU91nkinI/DxHZA2irqq+nNsSUqfDvhqoWAtOxP6jZpjzP41OsXeYPwXYjMCU4zng5kyRU9SdgFNA7ONUHGKGqi0XkGhGJJYZHsH7wiH0tOBT4d6rjTaYKPIsYCbasVN7nISJtsW+M74lIGxHpJCKXRBN1clTgWewlIjsExzsAVYGPo4g5mcrzPFR1o6ouim3ASmBjcJzxcqm6CeAS4C4R6QQ0w/6HBzgDGxD0nap+LSIPiMi/gvfuUdUZqQ816bb5LABEJA8bF9BYRHqoaraWNrb6PERkDVb/3hK4Mu7nTk9lkClSnt+N64DDRGQksATI1+ztT1+u/1eyVc6Mk3DOOVdxOVPd5JxzruI8STjnnCuTJwnnnHNl8iThnHOuTJ4knHPOlcmThHPOuTJ5knA5TUT+KiIPRh2Hc+nKk4RLChHZQ0ReE5HHReQLEVEROSvquEoxGngh6iCcS1eeJFyyvAHcoqoXAp2Bp+PfFJHLoghKRJoEE9IBoKpfq+qkKGLJJlH9e7rk8yThQhesi70r8CuAqhZhk55p8P4+wG0RxFUTS1a1U33vbCYi5wEnRR2HSw5PEi50qroCm8/mHRHpGpybD3wa/KE+D6gv1lpd1gAABf1JREFUIv8UkR4i0lxEngzaByaKyKGxa4lIXxG5XURuFJFCEZkuImeJyAUi8q2IHCIi40RkpYgcE/xMtWD+rb+JyCgRuSq43LHAnsApInKViHQQkeEi8mjc/Q4WkftE5B4RGSMiHYLzZd4vnoh0EZERInKTiAwVkbUi8qKI1BKRXYP/zkdFZJiIzA7O54nIzSLyUlBFt31wrZNF5GoReU5EJgfn6orIrSJynYgsDJ7PASLyvYjcHHzmJBHZKCI9K3rPrQmue7+I3Csi7wSx7AL0A9rF/Xs+Fdzv0iDG8SLSMHjWK0RkWMV+o1ykVNU330LfgI7AfGzp0yeAJnHvtbFfvd9e3wPcGBxfC7wVHDfAFrZpFLz+D/B8cNwYK5lch5UMrgM+DN47BhgXHHcB1sbdawJwXnBcE1vUfnjwuinwJVAteH0N8A02w2mZ9yvx310NmzV0PNAB2A8rUf0NqAUMDa65N9AfW9Htxbif/xC4Ozj+IO78oGB/KbZ+MkB3oG9w/BRwc9zn52NLb1bontv4N12GrXUOMAM4OTg+D1tfI/6ZTsJKk/WxCQAfBhoCOwGbgF2j/h31rXyblyRcUqjNnLs38ABwDjAjmEWzNEOBYSKyE/ZNv25wvkNwHFsZ7H8UVxUtC/YvqOp6bHW0ZsG5icA1IrIdcGDc9UrGuBH4Ke7UWcD3ausYAzwO7AEcvI37xV+zIPjsB6o6U1WnYuseH6OqG4Afgc9U9StVfQw4CthRRK4UkSuD/8aYliLyiIg0Ae4Pzv2CzUh6LPaHeFrs1iVDCeKp6D235hRVXSAi3bBk8LvnGvdMZ6nq96q6BluQ6BNVXaG23OtPpT07l55ybapwlwIiUhVoqarzgMuD6oXRWDI4oJQfWQQMAGYCk4FdgvNfAwuwufzfDM6/BFYMEdliiYsC7Bs/qrpaRA7AFqv/7zbCjf/j2gGo/tsbqqtEZAWw89buVw4zsGQVu1/8PVsBc1X1vlJ+7gRgODAHK2E9DDwDtAVewRa2OaUc96/IPbdmlYjcGdx7EWWvMVIyYRWWeF2EV3VnDP+HcslQnbg/Xqr6BTAQW72rNI8A81T1WWBj3M+tA64G+ojIOcCzqvritm4uIhcAXVV1ML9fxH5rFmAlmS0uR+LrBdTAvk2X5kfgaAlWPgMIEhzAGiAPq9p6QETaYdU9NwVx1gTuDD5bRPmT1tbuWSqxhYVGA7eq6pRy3sdlAU8SLlkuE5Hmca9bAO8Ex5sARKSBiOyOdZFtHFQPdQNqi0jb4PVtWHXNVOBLEakR/Gzsdzf+22zsuDNWnVKNYEUxEWknIrWCe8fuG/uZ2M89BbQQkR7Bz+wKzFTVKdu4X2l2ijvuCTwU9zr+j/looA4wSkSOEJH+FCfTa1S1SFWHYMuDApwmIjurra99a1wMS4AuQaP0SVhbQJOgVFfue4o19pc2nqUd1ka0o9jqdC2B7USkDcXPtFbwurSVDLf12qWrqBtFfMu+DWssVWwZx2exKpOngAbB+wK8B3yBNW72B1YDI4GTgZ+B47Fvyh8DK7DqHcVKBp2Cn1Gsa+3OwH1YKaQv9u37R+wP66FY1chtwb3/FFy/P9aAPhGrk+8YvN83uOc/sfaUnYPzZd6vlP/+4Vh7wdXAP4DLg/O7YY3EC4EecZ/vgVWtrQTuBaoE54uC5zYAuDI4dzMwGyuZDQZaB+dbYNVSM4P/5mnA7Vj1XkXuOQ2YU8p/U1XgXSwZXRc8m0+BJkBzbIW2N7BkMwlrKN8Ha8/5GRgBtAZODZ7bI0C9qH9Xfdv25ivTubQVfNvvo6oPBK+rYr2meqtVJaUlERmOVZ/dHHEolSIi16r+f3t3jMIgEAVheAZLT+Et0pjGW+UANt4lRfqQk4h4kpdCQZA8SLcr/F+33XTDvn2wMZXOgTowbkLNHpI62+1+brR9RP8qF+kvv8Ytl2D7rmNjCqAkULVR28bRanuR9JT0iYi5bKyc7V7STdJg+/wIfgVzRLxLh0A9GDcBAFLcJAAAKUoCAJCiJAAAKUoCAJCiJAAAqS896X2z26GRnwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(6,6))\n", "plt.plot(vehicle1.stag_pres_atm[0:2700], vehicle1.q_stag_total[0:2700], 'r-', linewidth=2.0)\n", "plt.xlabel('Stagnation pressure, atm',fontsize=14)\n", "plt.ylabel('Stagnation point heat-rate, '+r'$W/cm^2$',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/crew-module-atmospheric-re-entry-experiment-earth-thermal.png',bbox_inches='tight')\n", "plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth-thermal.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/crew-module-atmospheric-re-entry-experiment-earth-thermal.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 }