{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 29 - Mirka - Earth" ] }, { "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 Mirka entry probe.\n", "http://adsabs.harvard.edu/full/2006ESASP.631E..26B" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \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", "vehicle=Vehicle('Mirka', 154.0 , 214.9, 0.00, 0.785, 0.0, 0.5, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,7.60,0.0,-2.51,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (2400.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXxU1f3/8deHJWxBQMK+i7IICEpYFQEVcQFxrUu1tVVTrVr71VqXWtva2lZtq/6sS6O1rUVaVGwtuCAqsggKQRRFVmXfDJvsAZLz++PMJJNJQiZkJncmeT8fj/OYmTt37v3kip98cubcc8w5h4iIiIiIVE6toAMQEREREakOVFiLiIiIiMSBCmsRERERkThQYS0iIiIiEgcqrEVERERE4kCFtYiIiIhIHNQJOoB4ycjIcJ07dy53v9zcXFq0aJH4gFKYrlFsdJ3Kp2sUmwULFuxxzjUOOo6qFGvOjqZ/U7HTtaoYXa/Y1fRrtWDBgq3OuVIvQLUprDt37kxOTk65+2VmZsa0X02maxQbXafy6RrFxsxq3LeHsebsaPo3FTtdq4rR9YpdTb9WZramrPdqXDIXEREREUkEFdYiIiIiInFQ4wrrrKysoENIerpGsdF1Kp+uUcxygw4gVejfVOx0rSpG1yt2ulZlM+dc1Z7QbBTwG+By59zq0LaLgT8CxwDjgTucc4dD72UADwI7gbrAXc65Q9HHzczMdDV5vI+IpC4zW+Ccyww6jtIoZ4uIFHeknF2lPdZm1gpIBwZGbOsIXAhcCvwfcD3w44iPvQw875y7C1gOPFBlAYuI1GDK2SIiFVOlhbVzbgvwWtTmTsD3nXMLnHMvAE8CIwHMbDDQ3Tn3UWjfycCtZtaoqmIWESmPc45Dhw5RUFAQdChxlUw52znIz/dNRCRZVfkYa+dcQdTrWeGvEEM2AGtDz88AVkXsuyH0NCm/MhWRmunkk08mLS2Nzz77LOhQ4i5ZcvYHH0CdOnD66ZU9kohI4iTjzYsDgGdCz9sB26Le3w20jf5Qbm4umZmZhS07OzvBYYqIeHXr1gXg4MGDMX8mOzu7MF8BGYmJrEpUSc42849VfFuQiAgQe85OqgVizOx4INc592lokwP2R+1WDyhxI0yLFi1q9GTlIhKctLQ0APLy8mL+TFZWVuGd9Wa2NSGBJVhV5mwV1iISpFhzdtL0WJtZbeAHwN0RmzcAzSL2MaARsLFqoxMRKVu9evWAivVYp7qqztm1Qr+tVFiLSDJLmsIa+AnwB+dcZJfPNKB7xOvOwEFAXdMikjSOpse6GqjSnB3usa5m94eKSDVT5YV1qAcDwCK23QssBBqY2XFm9n0zO945lwNsMrPeoV3PAZ5wztWcbiERSXrhwro69lgnS87WUBARSQVVOsbazBoDV4deXmtmf8Z/lfjrqF2XOueeDz2/HLjfzFbj51P9WVXEKiISq+o6FCSZcraGgohIKqjSwto5txt4OtTCfhNqZX1mDXBdgkMTETlq1XUoSDLlbA0FEZFUkExjrEVEUlJ17bFOJhoKIiKpQIW1iEglVecx1slCQ0FEJBWosBYRqaTqOhQkmWgoiIikAhXWIiKVpKEgiaehICKSClRYi4hUknqsE09DQUQkFaiwFhGpJPVYJ56GgohIKlBhLSJSSbp5MfE0FEREUoEKaxGRStJQkMTTUBARSQUqrEVEKklDQRJPQ0FEJBWosBYRqSQNBUk8DQURkVSgwlpEpJLCPdYaCpI4GgoiIqlAhbWISCWpxzrxNBRERFKBCmsRkUrSzYuJp6EgIpIKVFiLiFSSbl5MPA0FEZFUoMJaRKSSNBQk8TQURERSgQprEZFK0s2LiaehICKSClRYi4hUknqsE09DQUQkFaiwFhGpJBXWiaehICKSClRYi4hUkoaCJJ6GgohIKlBhLSJSSeqxTrw6dfzj4cPBxiEiciQqrEVEKkk91olXv75/1CUWkWRWp6pPaGajgN8AlzvnVoe2NQAeBnYDzYC7nHO7yntPRCQZVOce62TJ2aG/XThwoLJHEhFJnCrtsTazVkA6MDDqraeA2c65e4HXgL/E+J6ISOCqa2GdTDk7dIk5dEg3MIpI8qrSwto5twWfaAuZWVvgKmByaNM7wAVm1vlI71VJwCIiMaiuQ0GSKWebFfVaV7PLLCLVSJUPBXHOFVj49m5vBLDVObcv9P5hM1sFDAcOHeG91VUZt4hIWaprjzUkV86uX98X1Xl50KBBZY8mIhJ/VV5Yl6IdsC1q226gLXD4CO+JiCSFOnXqYGYcPnyYgoICatWq1veFB5az69eHb77ROGsRSV7JkP0dsD9qWz18z8eR3ismNzeXzMzMwpadnZ2QYEVEoplZhXuts7OzC/MVkJG46OIusJytoSAiEpRYc3Yy9FhvwN85HqkxsBGfpMt6r5gWLVqQk5OTkABFRMqTlpZGXl4eBw8epH54brgjyMrKIisrCwAz25ro+OIosJwdvqz7o0t3EZEEizVnJ0OP9XSgbWiKJsysLtARmFHOeyIiSSN8A2N1HGcdJbCcfcwx/nGXJlwVkSRV5YW1Fd0FYwDOuc3AFODM0PazgEnOuQ1Heq8KQxYRKVd4KEh1mxkkmXJ206b+cceOeBxNRCT+qnQoiJk1Bq4OvbzWzP7snMsFbgIeNrO+QGsgK+JjR3pPRCQpVMce62TL2c1Cg0xUWItIsqrSwto5txt4OtQit+8AbijjM2W+JyKSLKpjj3Wy5exwYb1zZyKOLiJSeckwxlpEJOVV57msk4V6rEUk2amwFhGJg/Dc1c65gCOpvsKF9bbombJFRJKECmsRkThSYZ047dr5x/Xrg41DRKQsKqxFROIgatlvSYAOHfzjunXBxiEiUhYV1iIicaQe68Tp2NE/rl0bbBwiImVRYS0iIimhXTswg40b4fDhoKMRESlJhbWISByEh4Koxzpx0tL8cJCCAvjyy6CjEREpSYW1iIikjN69/eNnnwUbh4hIaVRYi4jEgXqsq0afPv5RhbWIJCMV1iIikjLChfWnnwYbh4hIaVRYi4jEQX5+PlC0UIwkxsCB/nHOHNCXAyKSbPQbQEQkDvbv3w9AgwYNAo6kejv+eGjVCnJzYfnyoKMRESlOhbWISBwcOHAAUGGdaGYwbJh/PmtWsLGIiERTYS0iEgfhwrp+/foBR1L9jRzpH994I9g4RESiqbAWEYkDFdZV54IL/ONbb8G+fcHGIiISSYW1iEgchMdYq7BOvPbt/U2M+/fDm28GHY2ISJG4FNZmdl48jiMikooOHjzIoUOHqF27NmlpaUGHU67qkLOvuMI/PvdcsHGIiESqUGFtZueY2WwzW2lmX4XaamBSYsITEUl+27ZtAyAjI6NwoZhkUJ1z9ne+A/XqwdSp8NVXQUcjIuLVqeD+/wR+BywECkLbDLgqnkGJiKSSrVu3Ar6wTjLVNmc3bw6XXw4vvAB/+hP8+c9BRyQiUvHC+ivn3J+iN5rZ0jjFIyKScsKFdfPmzQOOpIRqnbPvvBP++U949lm46y7o0CHoiESkpqvoGOubzOzSUrZfEI9gRERSURL3WFfrnN27t++1PngQfvGLoKMREal4YX0q8JyZ5Ue0AuDpBMQmIpIS1q9fD0CbNm0CjqSEap+zf/UrSEuDv/1NC8aISPAqOhTkASALyAHyQ9sMuKaygZhZT+AWYCVwApDtnPvEzBoADwO7gWbAXc65XZU9n4hIvHz55ZcAdO3aNeBISqj2ObtbN7j7bnjgAcjKggULoGHDRJ1NROTIKlpYvwm85pzLi9xoZs/GIZZ/AuOccxvMrCMwFegJPAW85ZybaGbnAH8BrozD+URE4uKr0LQUxx13XMCRlFAjcvY998BLL8HSpXDLLfD884k8m4hI2So6FOTfwD1mdnpEGw7cE4dYegKNQ8/3A03MrC3+7vXJoe3vABeYWec4nE9EJC7CPdZJWFjXiJxdvz5MnOgf//Y3yM5O5NlERMpW0R7rn+O/8rsWcKFtBrQGbqtkLK/gxwKeC1wN3AqMALY65/YBOOcOm9kqYDiwupLnExGptD179rBixQrq1q1Lt27dgg4nWo3J2SedBE89Bd//Ptx0E7RpA2PHJvKMIiIlVbTH+hGgg3Ous3OuS6h1Bi6OQyw/BA4A84E9zrlJQDtgW9R+u4G20R/Ozc0lMzOzsGWry0JEqsDChQtxztGnTx/q1asX8+eys7ML8xWQqOlEalTO/t734L77oKDAzxYybVqlDykiAsSesyvaY73OObe7lO2tK3ic0jQA/oVPzI+Fejkc/ivGSPWAQ9EfbtGiBTk5OXEIQ0QkdvPnzwcIJ9uYZWVlkZWVBYCZbY17YF6Ny9kPPABbtvi5rceMgX//Gy66KO6nEZEaJtacXdEe6x9GbzCzy4DHKnic0owHJjnnHgB+j0/Ym/B3lUdqDGyMw/lERCpt+vTpAAwePDjgSEpV43K2GTzzDNx6q5/f+tJL4ZFHwLnyPysiUlkVLayPNbOfA5hZQzN7DngO/1XgUTOzDKCvc25naNNvgF3AGqBtaPomzKwu0BGYUZnziYjEQ15eHu+99x4Ao0ePDjiaUtXInF2rFjz+uJ/juqAAfvpTuOoq2LOnKqMQkZqoQoW1c+58YJqZPQF8AnQD+gKlrexVEduBA2bWLmLbNmARMAU4M7TtLHwPyYZKnk9EpNJmzpzJvn37OOmkk2jbtsQw4sDV5JxtBvffD6++CunpfkhI377wwQdVHYmI1CRHHGNtZqeX8dYaYAXwB6ATcCHwf0cbhHOuwMwuBO43swVAK+BO59wuM7sJeNjM+uLHBWYd7XlEROJp/PjxAIwbNy7gSDzl7JIuugg++giuvBIWLYJhw+COO/wS6OnpQUYmItWRuSMMPDOzD/Fzle4ECsraDWjjnIv9dvgEyMzMdLp5UUSqyu7du2ndujX79u1j5cqVlVp10cwWOOcqdvdj6cdRzi5DXp4fGvLQQ354SLt28Mc/wre+5Xu3RURidaScXd5QkAeAE5xznSKmaopunYFvxztoEZFk9ve//519+/Zx2mmnJdNS5srZZahXD377W5g7FzIzYcMGuOIKGDFCw0NEJH6OWFg7595wzn1d3kGcc6/ELyQRkeR28OBBHnnkEQBuv/32gKMpopxdvoED4cMP/eqMzZvDzJlw2mlw3nnw8cdBRyciqa6is4KIiNR4zz//POvWraNnz55JM75aYle7NtxwA6xcCT//uR9r/eab0L+/L7CnT9f0fCJydFRYi4hUwPbt2/nZz34GwAMPPECtWkqjqappU7+gzFdf+RsaGzTwBfYZZ/ie7YkT4VCJpW1ERMqm3wgiIhXw05/+lO3bt3PGGWdwySWXBB2OxEGLFvCHP8DatfDLX0JGBuTk+DHYnTr5Xu3Vq4OOUkRSQaULazM72czGm9kdZta58iGJiCSnSZMm8de//pV69erxxBNPYCk4nYRydtkyMvw0fGvWwFNPQffusGkT/OY3cNxxcO658MorsD960XYRkZB49FjfCJzlnPsjMDS84paISHWyYsUKbrjhBgAeeeQRTjzxxIAjOmrK2eVo2BBuugmWLIH33/erNtatC2+9BZddBq1awbXXwttvw+HDQUcrIskkHoX1LUD/0PN/A3XjcEwRkaSxbds2zj//fHbs2MHYsWO55ZZbgg6pMpSzY2QGw4fDiy/66fkefdRP1bd7N/zjHzB6tJ8P++ab4Z13NB5bRI6isDazc83s26HnA/Fzpm4AvxqXc25XnGMUEQnMrl27GDt2LCtWrKBv3768+OKLKTUERDk7PjIy4Mc/hvnzYdkyPxa7Wzf4+ms/bGTUKGjZEq65xi+jvndv0BGLSBAqVFib2YPAy8BFAM65ecB3zOyqBMQmIhKob775htGjRzN37lw6dOjAlClTaNy4cdBhxUw5OzG6dfNjsZcuhQUL4N57oWdP2LkTxo+HSy7xhfi4cfDcc7B+fdARi0hVqWiP9VCgPTArYtubwO/jFpGISBLYsGEDI0eO5MMPP6RTp07MmDGD9u3bBx1WRSlnJ5AZnHIKPPggfPGFL7R//3sYNAgOHID//c/Pl92hA5x0Etx1lx+zffBg0JGLSKJUtLCe5pzbCUROnX9pHOMREQncwoULGThwIAsXLqRr1668//77dOnSJeiwjoZydhXq3t0Xzx9+6MdkP/00XHABNGoEn30GDz8MI0f63uyLLoK//AVWrNBiNCLVSZ0K7v+Vmd0LtDOzK4DL8F8x3hr3yEREqphzjueee47bbruN/fv3M2zYMF599VUyMjKCDu1oKWcHpG1buPFG3/LyYPZsv/jMm2/63u3//tc3gPbtfcE9cqRfnKZTp2BjF5GjZ66Cfyqb2WDgWqAjsAn4l3PunfiHVjGZmZkuJycn6DBEJEXt3LmTH/zgB7z00ksAXHfddTz55JPUq1cv4ec2swXOucwEHVs5O8msXeun7ps2zQ8N2bq1+Ptduvgie8QIOO006NzZDzsRkeRwpJxd4cK6lIPXBkY4596t1IEqqSYnaRGpnP/+97/88Ic/ZNOmTaSnp/PMM8/w7W9/u8rOn8jCupRzKWcnkYICWLwY3nsPpk+HGTP8TZCRWreGU0/1behQOPlkSEsLJl4ROXLOPuJQEDP7H5BezvEzgK1AoElaRKSiNm7cyG233cYrr7wCwODBg3nhhRc44YQTAo7s6Chnp55ataBPH99uuw3y8+HTT32hPWsWfPABbN4Mkyb5BlC/Pgwc6IvsoUP981atgv05RMQrb4z1otA+y4AC4AZgIRDZzdAFODYh0YmIJMC+ffv44x//yEMPPcTevXtp1KgRv/3tb7n55pupXbt20OFVhnJ2iqtd2880csop8JOf+Bsbly/3BfacOf5x6VKYOdO3sA4dYMCAopaZCU2aBPdziNRU5RXWj+CHi+wEMLOuzrn7o3cys9cSEZyISDzl5+czYcIE7r33XtaHJhceN24cjz/+OJ2qxx1jytnVjJmfbaR7d/j+9/22rVth7lxfaM+d6+fSXrfOt1dfLfpst27Fi+1+/fxy7SKSOEcsrJ1z30RtamVmtZ1z+eENZjYCGJiA2ERE4iI/P5+JEyfywAMPsGzZMgD69evHn/70J0aOHBlwdPGjnF0zZGTA2LG+gR8+smyZXxUy3D75xPd0L1/ul2QHP+ykWzdfYIdb375+DLeIxEdFp9t7BfjUzGYDu4GewGigRI+IiEjQDh48yMSJE3nwwQcLC+ouXbpw//33c80116T6sI9YKGfXALVrw4kn+vbd7/ptBw/6ubPDhXZOTtEiNkuXwr//XfT5Vq2KF9v9+sEJJ/jjikjFHM10e22BHwInAbuAl51zgX+tqDvMRSRs+/btZGdn88QTT7Bx40bAF9T33Xcf11xzDXXr1g04wuISPN2ecrYAfj7tL77wvdmRbdeukvs2bAi9ekHv3sVbmzaa+k/kqGcFKY1zbiNwX8TBa5lZf+fcgkrEKCJSaV988QVPPvkkf//739m3bx8AvXr14o477uDqq69OuoK6KihnS1i9en6qvpNPLtrmHKxeXbLYXru2qLc7UrNmJYvt3r3hWN0OKwJUsLA2s+dL2dwCyAcujEdAocUMTgNWArOBvcDD+K8xmwF3OedK+ftaRGqi/fv38/LLL/Pss88ye/bswu2jR4/m9ttvZ9SoUVgN7WJTzpbymPkFabp08cush23f7ufX/vzzovbZZ7Bjh58GcNas4sdp06Z4oX3iidCzp2YmkZqnoj3WJ+DnPo0cP9IemFX67hVjZt8DujrnIntX/ga85ZybaGbnAH8BrozH+UQkNTnnWLRoEX/961/55z//yc7Qihrp6elcddVV3HrrrfTu3TvgKJOCcrYclWOPhWHDfAtzzs+pHVlsf/65L8A3bfJt2rTix2nTxhfYPXr4x3DTkBKpripaWF/vnFsWucHMeuCXy60UMxsOfBsYFbGtLXAVcHNo0zvAJDPr7JxbXdlzikhqWbt2LRMmTGD8+PEsXry4cPuAAQPIysriiiuuID29vPVRahTlbIkbM18Qt2kDo0YVbS8ogDVrivdsL1niZyoJF9zvvVf8WMccU7LY7tnT95zXqfAgVZHkUaF/vtEJOiQPuB64u5Kx/AmYAzxhZscDvwY6AVudc/tC5z9sZquA4cDqSp5PRFLAjh07eOWVV3jxxReZMWNG4fbmzZtz+eWXc/3113Ny5KBRKaScLVWhVq2i4SThKQChqOBessS3pUuLnm/fDvPm+RYpLc3PSBIutHv08K1bN9DfzJIKKjrGehXFv1KsDbQC/lmZIMysO3Ay8D3n3CIzuxN4C/gNsC1q991A2+hj5ObmkplZdINmVlYWWVlZlQlLRAKydetWXnvtNV555RXeffddDh06BED9+vW54IILuPrqqxk9ejRpaWkBR3r0srOzyc7ODr/MSMQ5lLMlSJEF93nnFW13DnJzSxbbS5b4RW4WL/YtWrt2fqGcbt2KFs3p1g06d9bUgJJ4sebsCk23Z2bZwASKEnUBsMU5t/wo4wwf92Ig2zmXEXrdAMgFHgAucc4Nitj3Y2CCc+4PkcfQ1E0iqW3Tpk385z//YdKkScyYMYP8fL+mSa1atRgxYgTXXHMNF198Mcccc0zAkcZfoqbbU86WVLNnjx9CEllsL1sGK1f6ublLk5YGxx9fetGdkZA/WaWmi+d0e+OdczOjDt7dzI53zq086gh9HIWxOOf2m9kKfO9Ks6h9GwMbK3EuEUkCzjmWLFnClClTmDx5Mh988AHhP/Tr1KnD6NGjueSSS7jwwgtp0aJFwNGmLOVsSSnp6dC/v2+R8vP9sJJly4ra8uX+ccMGPz/3F1+UPN6xx5ZecB9/PNSvXzU/k9QsFS2szwFmRm3bAkzHfy14tBYBTcyslXNuS2jbYXwybmtmDUKJuy7QEZhR1oFEJHnl5eUxY8YMpkyZwpQpU1i1alXhe/Xq1ePss8/mkksu4YILLqBZs+j6TI6CcrZUC7Vrw3HH+XbuucXf27MHVqwoWXAvW+bHcs+d61skMz+EJLLgDhfd7dr5YSwiRyOmwtrMbsGPnWtsZneVssu0UrbFzDm31MzeBC4FnjSzpkAa8CJwLnAmMAU4C5jknNtQmfOJSNXZtGkTb775JlOmTOHtt99m7969he9lZGRw/vnnM2bMGM4+++xqOcwjCMrZUpOkp5dc+AaKpgcsreBetaqoTZ1a/HMNG/oCO7qXu3t3P5uJyJHEVFg75/5sZjOB7wL/L+rtg865TXGI5TvA42aWhu/huDJ0R/lNwMNm1hdoDejuFpEktn//fmbNmsXbb7/N1KlT+fzzz4u9f9JJJzFmzBjGjBnDwIEDqa27juJOOVuk+PSAI0YUf+/gQfjqq5IF9/Ll8PXXRStQRmvduqjIjiy8u3SBGriwq5Siojcv1nXOHSple0fn3Nq4RlZBuhFGJBjOOT7//HOmTp3K22+/zcyZM8nLyyt8v2HDhowYMYIxY8Zw/vnn07FjxwCjTU4JvHlROVukgnbs8AV2dMG9fDkcOFD6Z+rU8cNUShvP3aqVFsOpbo765kUz+wnwhnMufEvAT80sunvJgJHAiMoGKiKpYf369UyfPp133nmHt99+m82bNxd7/5RTTuHss8/m7LPPZujQodSrVy+gSGsW5WyRymvWDAYN8i1SQYGfDjC64F62DNauLSq+ox1zTPGCO/x4wgnQqFHV/ExSdcobCtIFaBrx+hSgEbCZoumbDP91n4hUU5s3b2b69OmFbeXK4hNKtG3btrCQPvPMM2nZsmVAkdZ4ytkiCVKrFnTq5FvkypMA+/f7KQGjh5YsWwY7d8L8+b5Fa9++5Dju7t2hY0fNzZ2qjlhYO+dujtr0IPC5c67YbJJmdmK8AxOR4OTm5vL+++8XFtJLly4t9n7jxo05/fTTOeOMMzj77LPp1asXpu86A6ecLRKMBg2gTx/fIjkHW7eWPpZ75UpYv963d98t/rl69UrOzR1+bN686n4uqbiKTrfXxzn3ceSG0ApcXYFSZpAUkVSwZcsWZs+ezYwZM5g+fXqJGw4bNWrEaaedxsiRIxk5ciSnnHIKdepUNH1IAJSzRQJkBi1a+HbaacXfO3y4+NzckYX3xo1lr0DZvHnpM5Z07aq5uZNBrNPtDcN/fTgqtERupGPxd51PjnNsIpIAzjm+/PJLZs+ezaxZs5g1axYrVqwotk/9+vU59dRTCwvpAQMGUFe3vKcM5WyR5Fenji+Gu3YtvuQ7+Lm5w4V29OO2baXPzR0eqlLaDZSam7vqxNrltBN4AT9+79So9w4Cz8UzKBGJn/z8fBYtWsSsWbMKi+nomw0bNmzIkCFDGDZsGCNHjmTQoEG64TC1KWeLpLD0dDjlFN8iOQebNpV+A2Xk3NxvvVX8c9Fzc0cW3pqbO75incf6MzMbAVzmnFNCFkli+/fvZ/78+YW90XPmzGH37t3F9snIyGDYsGGcdtppDBs2jH79+qlHuhpRzhapnsygbVvf4jU3d6tWxXu3e/SA3r39DZTq5a64mAdJOue+oYxeDjM7zTk3O25RiUjM1q1bx9y5c5k7dy5z5sxh4cKFHDpUfOriLl26MGzYsMJiunv37rrZsJpTzhapWdLSfFHco0fJ93buLL2Xe/ly2LLFt5kzi38mPR169fJFdmTTvNxHVt481v8D0ss5hgG9AM2vJZJgeXl5LFy4sFghvWFD8dWizYyTTjqpWCHdrl27gCKWqqScLSKladoUBg70LVJBgZ+VJLLg/uILf9Pk5s3w0Ue+RWrevGSx3auXn/9byu+xXoYfq7eBojlQoxkwJp5BiYi3cePGwiJ67ty5LFiwoNiqhgBNmjRhyJAhhW3gwIE0adIkoIglYMrZIhKzWrX8kI+OHUvOzb11qy+wP/+8qH32mb95csYM3yK1a1ey4O7Zs+YtglNeYf0QsNs5l3ekncxMd5eLVNLBgwdZtGhRYU/03LlzWbNmTYn9evbsyZAhQxg6dChDhgyhR48e1NJAOPGUs0UkLjIyYPhw38Kc81MBRhbbn3/uC/ANG3ybOrVofzO/1Ht0wd2tmx+6Uh2Vt6WSee8AACAASURBVEDM1hiPcxG6y1wkZuEp7+bNm8dHH33EvHnzWLhwYYne6PT0dAYPHlzYGz148GCa6fs2KYNytogkkpnvmW7XDkaPLtqenw+rV5csuJcuhS+/9O2114r2r1PH3ywZXXB36ZL6K05WaoUHM+sI3AzcgpK0SJlyc3OZN29esbZ9+/YS+3Xr1o3BgwcX9kb36tWL2qmeZSRpKGeLSCLUrl00J/e4cUXbDx7047ajC+6vvipaAGfixKL9GzSAE0+Evn1h8GDfTjwxtYrtChfWZlYbuAC4ARgFbMGP6RMR/HR3H3/8cbHe6FWrotfogJYtWzJo0CAGDhzIwIEDGTBggHqjJe6Us0UkKGlpRb3RkfbuhSVLShbcGzbAggW+Pf+83zc93d90GS60hwzxw1SSVcyFtZl1Aa4Hvgc0B9YAY4GpwJCERCeS5PLz81m6dGmxInrRokXk5+cX269hw4b079+fgQMHFhbTHTt21JR3kjDK2SKSrBo1gsxM3yLt2OEL7I8/hg8/9G31anjvPd/C+vSBM86AkSP9GPCmTas0/CMy58q6cTy0g9ll+J6OM4CV+K8P/wHc65z7v4RHGKPMzEyXk5MTdBhSjTnnWLt2LfPnzycnJ4d58+aRk5NTYvGVWrVq0atXr2K90b169aJOnUqNvJJqzMwWOOcyy98zpmMpZ4tItRGe9u/DD/0y7h99BAcOFL1v5leoHDUKxo6FQYMSP3TkSDk7lt/0+/FL4M4BrnfOLQ8d9MgVuUiK27x5c2ERHX7Mzc0tsV+HDh2KFdH9+/cnPb28qYRFEkY5W0Sqjdat/bjt8NjtAwd8cf3eezB9ui+4w8NHfv97P0zk/PN9kT16tB9KUpXK7bEu3NHf9HID0Bl4DRjunLs19J65WA+UIOr9kMrYvn07CxYsYP78+YVF9Pr160vsd+yxx5KZmcmAAQMKx0W3adMmgIilOolnj3XEMZWzRaTa27sXPvgAXn8dJk+GyFuaGjSACy+Eq6/2Pdp168bnnEfK2TEX1hEHqwNcCHwb+Bz4K3C2cy67soFWhpK0xGr37t0sXLiwWBH95ZdfltgvPT2d/v37M2DAAAYMGEBmZiZdunTRuGiJu0QU1hHHVs4WkRrBOb9y5OTJ8L//+aEjYRkZcO21cOutfkGcyohrYR114OOBW4EbnHMNj/pAcaAkLaU5cOAAn376abEiesmSJUT/u69fvz79+vUrVkR3795dC69IlUhkYR11HuVsEakxVq+GCRNg/Hg/Cwn48ddXXOGHjbRvf3THTVhhHXGC7zjnXqj0gSpBSVoOHTrE4sWLi42L/uyzzzh8+HCx/erUqUOfPn2KFdG9evWibry+IxKpoKoqrCPOp5wtIjWGczB/Pjz2GLz0kl/QplEj+Mc/4JJLKn68yt68WK54JWgzSwPmA7c55943swbAw8BuoBlwl3NuVzzOJamtoKCAZcuWFSuiP/nkEw5E3ioMmBknnnhiYQE9YMAA+vbtS/369QOKXCR48SyqlbdFJNmZ+bmwJ0yA3/4W7rgDXn3V91wvWgQ9e8bvXMk2/9ed+Bttwp4C3nLOTTSzc4C/AFcGEZgExznHqlWrihXRCxYsYM+ePSX27dq1a7Ei+pRTTtEMHSKJpbwtIimjc2d45RUYMQJmzoR3362mhbWZDQU2ATtCr9sCV+GX3wV4B5hkZp2dc6sDCVKqxKZNm5g3b16xcdGlLf/dvn37YkV0//79OfbYYwOIWKRmUt4WkVRz+DA8+ijMnu1fDxoU3+MnRWFtZo2AS51zt5vZ/aHNI4Ctzrl9AM65w2a2ChgOrA4kUIm7Xbt2sWDBAubNm1fYSpvmrkWLFsWK6MzMTFq3bh1AxCICytsiklry8/1MIfff71d3BH8D44AB8T1PUhTWwF3A76O2tQO2RW3bDbStkogk7g4ePMhnn31WrIgubYaOxo0bF95YGC6itfy3SNJR3haRpLdmDfzrX/Dss/DVV35b587wxBMwZkz8zxd4YW1mo4Ec59zXUW85/ApikeoBh0o7Tm5uLpkRi85nZWWRlZUVz1ClApxzrFy5slgRvXDhQvLy8ortV7duXfr27Vu4auHAgQM1zZ3UCNnZ2WRnF04lnRFkLBUVj7ytnC0iiVBQAB9/DFOnwhtvwJw5Re916QI//jFkZUFF5zCINWfHZbq9yjCzacDJEZuaAXuAPwHfds51i9h3BfAL59yE6ONo6qZghZf/DhfR8+fPZ8eOHSX269atW7EiWjN0iFT9dHuVFY+8rZwtIvGQl+cL6TlzfJs5E7ZuLXq/QQO/HPpVV8F55/l5rCsr4dPtVdJV+B6NsLnA7cBHwJ1m1sA5t9/M6gIdgRkBxCgR9u7dS05OTrHe6LVr15bYr1WrVgwaNKiwiM7MzKRZs2YBRCwicaa8LSJVbs8ePz3ep5/69sknsHAhHDxYfL+OHeGcc3wbNQqqcnKwwAtr51xu5GszywdynXNrzWwKcCYwBTgLmOSc2xBAmDWWc441a9YwZ84c5syZw9y5c/n000/Jz88vtl96ejqZmZnFeqPbt2+vcdEi1ZDytogk0u7dsGyZXy1x6VLfFi2CL7/0i71E69ULhg6FIUPg1FPhhBP83NVBCLywLsdNwMNm1hdoDWgAXoIdOHCAjz/+mLlz5xYW05s3by62T+3atTn55JOL9Ub36NGD2vH4fkVEUp3ytoiUKy/P31j45Ze+hQvopUthQxl/itetCyeeCH37Qr9+/vHkkyGZvgxPusLaOdc54vkO4Ibgoqn+Nm7cWFhEz507lwULFnAw6juVY489liFDhjB06FCGDBnCgAEDtOiKiBRS3haR0mzf7mfiCBfP4edffQXr1pXe+wxQrx506wY9ehS1Xr38Qi5paVX7M1RU0hXWkjjOOVasWMHMmTOZMWMGs2bNYs2aNSX269WrV2ERPXToULp166YhHSIiIlIoPx82bYK1a31bs6bk465dZX++Vi3o1Am6doXjjoPu3YuK6M6d43OTYRBUWFdjBQUFLFmyhBkzZhQW09HDOho3bsygQYMYOnQoQ4cOZdCgQTRt2jSgiEVERCRozvmieP163yIL5vDz9ev9KoZH0qiRL5zDxXPkY8eOyd/7fDRUWFcj+fn5LFq0qFiP9NbIOWfwKxiefvrpDB8+nGHDhtGnTx+NjRYREakh8vPh6699Ybxhg2+lPd+7t/xjtWzpe507diz52LEjNG8e3E2EQVFhneK++uorpk2bxjvvvMO7775bYu7otm3bMnz48MJiukePHhrWISIiUg3t3euHZ4SL5NKK5k2bfHFdnoYNoX17aNeu9OK5fXs/R7QUp8I6xezYsYP33nuPadOmMW3aNL4Kr88Z0qlTJ4YPH15YTHft2lWFtIiISArbs8cXxBs3Fn+M3nakMc2RWrTwBXO4cA63yNdNmtS83uZ4UGGd5AoKCsjJyWHKlClMnTqVnJwcCgoKCt9v2rQpZ555JmeddRajRo2ia9euAUYrIiIisdq9u2SxXFrRvHt3bMdLS4M2bUoWyZHP27b1s25IYqiwTkJ79uxh2rRpTJkyhddff50tW7YUvle3bl2GDRvGqFGjGDVqFP3799cYaRERkSQRvvGvtGI5umiOZRwzQP36vmBu27b0x/DzZs3Uyxw0FdZJYtOmTfznP/9h8uTJvPfee8Xmku7YsSNjx47lvPPOY/jw4TRq1CjASEVERGqmXbt8QRzdoovmfftiO16DBiWL49IemzZVwZwqVFgHaMuWLUyaNImXXnqJmTNn4kIzpZsZQ4YMYcyYMYwdO5bevXtrnLSIiEiC7N9fVBxv3Ohv9CutgN6zJ7bjNWxYdpEc+XjMMSqYqxsV1lVsx44dvPTSS0ycOJEZM2YUjpdOS0vjnHPO4aKLLuK8886jZcuWAUcqIiKS2g4dgi1bjlwwb9gAURNqlSncwxweqxwukKOL5saNVTDXVCqsq0B+fj7vvvsuf/vb3/jPf/5DXl4e4MdLn3feeVx++eWMHTuWJk2aBBypiIhI8nMOtm4tKpTLKpi//rrsZbMj1alTVCBHtnABHW6aKUPKo8I6gdavX092djZ///vfWbduHeCHeZx11llcffXVjBs3TqscioiIRPnmG1i3rqitXVv89fr1cOBA+ccxg9atyy6Uwy0jwy+xLVJZKqwT4KOPPuKxxx7j5ZdfJj80C3uXLl343ve+x3e/+106duwYcIQiIiLB2L+/eJFcWgEdy/RyTZsWn4O5tN7mVq18b7RIVdE/tzhxzvHaa6/x0EMP8eGHHwJQu3ZtvvWtb3HjjTcyfPhwaunPYRERqeYOHvQF8qpVpbfc3PKP0aCBX+GvQ4eiFv06PT3xP4tIRamwriTnHJMnT+aXv/wlCxcuBKBZs2ZkZWVx880306FDh4AjFBERiZ+CAj+GuazCecMGv09Z6tTxC5YcqWg+9liNZZbUpMK6EubOncttt93G/PnzAWjTpg1333031113neaaFhGRlOWcn01jxQpYvty38PMvvzzy+OZatXyh3KULdO7sHyNbmzagdc2kulJhfRQ2b97MXXfdxQsvvABAq1atuOeee8jKyqJBgwYBRyciIhKbAwdgyRJYvLh48bxixZHHObdsWbJgDrcOHfzS2iI1kQrrCnDOMWHCBG6++Wa++eYb0tLSuPPOO7nnnnvUQy0iIkkrP9/3NH/2GXz+edHjihVlD9to2hS6dStqJ5zgH48/3i9sIiIlqbCO0c6dO8nKyuLll18G4Nxzz+WJJ56ga9euAUcmIiJSpKAAli2D+fOL2qeflj58o1Yt6N4devf2j5FFdPPmGucsUlEqrGOwbNkyLrjgApYvX056ejqPPvoo1113nZYZFxGRwO3ZA3PmwPvvw9y5sGBB6cM4OnSAPn18Ed27t3/eowfUr1/lIYtUWyqsyzF9+nQuuugivvnmG0466SReffVV9VKLiEhgDhyAmTNh+nRfTOfkwOHDxfdp3x4GDChq/ftDs2aBhCtSo6iwPoK3336bcePGceDAAS6++GL+8Y9/kK6JM0VEpIpt3gxvvAGTJ8O0abB3b9F7tWr54nn4cBg2DAYO9KsNikjVS5rC2swuBv4IHAOMB+5wzh02swzgQWAnUBe4yzl3KNHxzJw5kwsuuIC8vDyysrJ4+umntcCLiEhIsuXs6mj7dnj5ZRg/HmbPLv5ev34wahSMHAmnnqqbCUWSRVIU1mbWEbgQuBToBTwNrAP+ALwM3O2c+8jMbgQeAO5JZDwrVqzgoosuIi8vjxtvvJEnn3xSRbWISEiy5ezqxDk/xOOJJ+D11+FQ6E+SevXgzDNh7Fg4/3w/XlpEkk9SFNZAJ+D7zrnDwAIz6w2MNLPZQHfn3Eeh/SYDy8zsN865vWUdrDL279/PRRddxPbt2xkzZgx//vOfVVSLiBSXNDm7uigogIkT4eGH4ZNP/LZateDss+Hqq+HCC6Fx42BjFJHyJUVh7ZybFbVpA9AYOANYFbHfhtBMHJnAjETEcu+997J48WK6d+/OhAkTqK3loUREikmmnF0dzJ4NP/6xn80DoFUruOUWuO46v0qhiKSOZO2KHQA8A7QDtkW9txtoG/2B3NxcMjMzC1t2dnaFTzp37lwee+wx6tSpw/jx42ms7gERSZDs7OzCfAVkBB1PJQWSs1PdgQNw++3+hsMFC6BtW3j2WVizBu67T0W1SDKJNWcnRY91JDM7Hsh1zn1qZg7YH7VLPaDEjTAtWrQgJyfnqM9bUFDAj3/8YwDuvPPO8IUTEUmIrKwssrKyADCzrQGHc9SCytmpbuVKuPhivwJi7dpwzz1w992gRXxFklOsOTupCmszqw38ALg7tGkD0C3ifQMaARvjfe6XXnqJefPm0aZNG+699954H15EpNoJMmensjffhKuugp07/QqHL77op8sTkdSXbENBfgL8wTmXF3o9Dege8X5n4CAQ124O5xy/+93vAPjVr36luapFRGITSM5OVc7BQw/5WT127oRx4/ziLiqqRaqPpCmszexeYCHQwMyOM7Pv4+dB3RS64xzgHOAJ59zBeJ576tSpLFq0iDZt2vCd73wnnocWEamWgszZqejwYbjpJj/cwzn41a/g1Vc1/7RIdZMUQ0HM7D7g11Gblzrnnjezy4H7zWw1kA78LN7nf+aZZwD40Y9+RL169eJ9eBGRaiXonJ1q9u2DK6+E//3Pz0c9YYIfXy0i1Y8554KOIS4yMzPd0dwIs3XrVtq2bUt+fj4bNmygtdaBFZEqZmYLnHM16o7po83ZqWbfPjj3XJg5E5o188X1aacFHZWIVMaRcnZS9FgH6eWXX+bQoUOce+65KqpFRCRuDhzwC7vMnAnt2sG0adCzZ9BRiUgi1fjC+vXXXwfgsssuCzgSERGpLpyD66/3xXTLlvDuu9C9e/mfE5HUljQ3LwbhwIEDvPfeewCcc845AUcjIiLVxTPP+Gn0GjXyxbWKapGaoUYX1h988AH79++nX79+tNESVyIiEgerVsEdd/jnzz4LJ50UbDwiUnVqdGE9Z84cAEaMGBFsICIiUi04B7fcAvv3wxVX+NlARKTmqNGF9YcffgjAoEGDAo5ERESqg3ffhTfegCZN4NFHg45GRKpajS2snXN89NFHgAprERGpPOfg5z/3z+++GzTRlEjNU2ML63Xr1rFt2zYyMjLo3Llz0OGIiEiKmzEDPvwQMjL8cBARqXlqbGG9ePFiAHr37o2ZBRyNiIikuuef94833gjp6cHGIiLBqPGFda9evQKOREREUt3evfDKK/75974XbCwiEpwau0BMp06dGDNmDKeeemrQoYiISIrbuxeuvRY2b4bjjgs6GhEJSo0trC+77DKttigiInHRsiU89VTQUYhI0GrsUBARERERkXhSYS0iIiIiEgcqrEVERERE4kCFtYiIiIhIHNS4wjo7OzvoEJKerlFsdJ3Kp2sUs4ygA0gV+jcVO12ritH1ip2uVdlUWEsJukax0XUqn65RzFoEHUCq0L+p2OlaVYyuV+x0rcpW4wprEREREZFEMOdc0DHEhZnlAmti2DUD2JrgcFKdrlFsdJ3Kp2sUm+7OucZBB1GVKpCzo+nfVOx0rSpG1yt2Nf1adXLOlfpNY7UprEVEREREgqShICIiIiIicaDCWkREREQkDmpMYW1mDczsCTP7rZk9bWbHBB1TsjGzi81slZltM7PHzaxO0DElMzNLM7NPzWxE0LEkKzMbbGY/MbMLzUxTykUxs55m9qSZ/Z+ZPWVm/YKOKVkph8dOubzilM8rRrm9bDXpf7angLeccxPN7BzgL8CVAceUNMysI3AhcCnQC3gaWAf8Ici4ktydQOegg0hWZvY9oKtz7r6gY0li/wTGOec2hP4fnAr0DDimZKUcHgPl8qOmfB4j5fYjqxGFtZm1Ba4Cbg5tegeYZGadnXOrAwssuXQCvu+cOwwsMLPewEiUjEtlZkOBTcCOoGNJRmY2HPg2MCroWJJcTyA8G8h+oEmAsSQt5fAKUS6vIOXz2Cm3l6+mDAUZAWx1zu0DCCWcVcDwIINKJs65WaHrErYBWBtUPMnMzBoBlzrnng86liT2J2AJ8ISZvWVmpwYdUJJ6BXjOzBoDVwO3BhxPshqBcnhMlMsrRvm8wpTby1FTCut2wLaobbuBtgHEkioGAM8EHUSSugv4fdBBJCsz6w6cDDzrnLsFeBd4y8xaBRtZUvohcACYD+xxzk0KOJ5kpRx+9JTLj0z5PEbK7bGpKYW1w3/NGqkecCiAWJKemR0P5DrnPg06lmRjZqOBHOfc10HHksR6Adudc4tCr/8MGDAuuJCSVgPgX8AE4DEzOyvgeJKVcvhRUC4/MuXzClNuj0GNGGON/yqsWdS2xsDGAGJJamZWG/gBcHfQsSSpnwAnm1n4dTPgNTN70Dn3cHBhJZU6ROQW59x+M1sBNA8upKQ1HrjCObfTzBzwr9C44b1BB5ZklMMrSLk8JsrnFaPcHoOa0mM9HWhrZg0AzKwu0BGYEWhUyeknwB+cc3lBB5KkrgL6RbSNwPXoq9ZIi4AmUV8PHga+CCiepBSaoqqvc25naNNvgF1oVpDSKIdXnHJ5+ZTPK0a5PQY1orB2zm0GpgBnhjadBUxyzm0ILqrkY2b3AguBBmZ2nJl9P/RVooQ453Kdc+vDDcjHf9W6K+jYkoVzbinwJn66L8ysKZAGvB5kXEloO3DAzNpFbNsGLA8onqSlHF4xyuWxUT6vGOX22NSUoSAANwEPm1lfoDWQFXA8ScXM7gN+HbV5qe6UlqP0HeBxM0vD9yxeGTVTQY3nnCswswuB+81sAdAKuFO/1MukHB4D5XJJMOX2cphzLugYRERERERSXo0YCiIiIiIikmgqrEVERERE4kCFtYiIiIhIHKiwFhERERGJAxXWIiIiIiJxoMJaRERERCQOVFiLlMHMzjSzOWY2opz9epvZS2b28zie+5F4HUtEJAhmdpKZvRzP3FjGeaaY2RWJPEeimNlZZja3vN8zkjpq0gIxkqLMrCfwIH6lugHAScA1zrnxCT71OqBvDPutBZoBteNx0tByzYficSwREQAzGwD8DjgdeBtYhV9opy7wtHNuagJOuxJoQpxy4xE8BXyeiAOH8vHF+GXOlwNzQ2/VAy4BBjrnVlfiFLH+npEUocJaUsF/gW855z41s1rA36ripM655Wa2LYb9dpnZpjieejgwM47HE5Eazjk338z+DfRwzo0Jbzez84H/mNlPnHNPxfmc++KcGwEws7HAZ+GC1jn3RrzPEeacOwRMNLNbgfedc/dFxPHa0RzTzFoCI51zE51zy8xsa5zClSSgoSCS1MzsWKAbsBf8MtDA/UBVLRlaEOf9YjESmB7H44mIAJRYeto59zrwS+BRM2udgHPGNVebWXt8D3VVK23Z7neADRU5iJnVA14AGkRsjufvDwmYCmtJas657cAy4C0zGxzatgaYZ2ZNzezXZjbPzC4ws81mtsTMTgl/3syuNLOfmdn7ZvbbiO2ZZvZLM5toZq+aWZPQ9vpm9qSZ3W9mf8V/jVkhZvbTUEzfMrMBZvbf0PF+HYpxopm1CG3/xsx+EXWINOdcXuhYd5jZ7Wb2gZn9vsIXUESkfH/FDwm5DMDMuprZA2b2nJm9GypmMbN6oVz2MzObZWaXhw9gZj80s5+b2YdmdktZJyotJ5vZpaGcea2ZrTCz35tZWzP7WyifzjKzM0KHuApoD9xlZheG7oV5x8y+G3GOG8zsQTN7wcz+YWaNzKyJmf3CzHLM7AwzW2ZmX5lZ56O5YKFrMiLUo42ZXWxmj5vZH83sLTNLD/0Mj5vZ5NDvm02ha3wicJmZ3R5xyGah30W7zeyGMs5Zy8weMrOs0Dmcmf0n/N9HkoRzTk0tqRvQC1iD/6v+eaBlaHttfJLajh/r1hSYBizG/9E4BLgvtG8zfK/3+cAxwL8jjv8B8Ejo+aNAVuh5KyAfnzzLi/Hv+F6flvixeA1D2+sBU4HXgLb4XwiHgF8ADYH+odfh/cPjx8M/91Oh542Bu4P+b6Gmppa6DbgWWF/Ge1uA/xfKq68CtUPbJwAvh57/Cbg49Pz7QG7o+ZXA1aHn3UK5uk/o9d+BX4ael5WTjwV2Ao8BA4GzgD8A94f2vROYHBGrAzqHnjcAlgDXhl6PBiZE7DsF+At+6OtFod8XF4d+zreAB2O8du+Hflf8HngEP976rIj3twIdQ88X438n1QN+BqwHBgHXhfL+++F4Q/uvBv4Yeu9KYG0ZMVwPTAk9bxW6zmcH/e9KrXjTGGtJes65xWbWB/gN8ENgrJmd5fyY61xgl3NuEoCZ/RKYDRyH/yVSYGY/Dh1qKj6Bnw80j9i+JPTZdOCm0Gdxzm0xs3UVCDUDeAK4zjm3L3SMPDPbDKxyzm0MnWczMMP58YcL8Qk/A38T5Pn4XwIA+4DvmNly4Gn8HxUiIong8MXmQKATcKuZAWwDmph/8QN8pwDAP/BDIcDn2pyoXNsG+CzqHNdSSk52zm03s2+A15xz8wDMbC2w18za4Ht400sN2rn9ZvZ1xKZbQscNew54KbR9J/73xauhc+SE4ozVdBcaY21mE/GdOWGXOefWmtmp+II6PZT/NwJfOuc+Aj4Kfba0Y08O/U6Yf4SYTgH2h37uLWa2neJDSiQJqLCWpGZmtYEOzt+k8qPQ8Iw3gGx8D0C0xaHHpkBH4J/OuQmhbY+FjnkXvtB9LOpc/fEJcfdRhtsNOA3oCnwasT16jGF+4RvOFYSSbHhYVnPnh7/gnFtlZtcCfwb+D/g28DUiInFkZk3x37YtxufNbaXkx1b4HtU0AOdcPr4zgNBnfuWcmxN6XeyzEUrNySGO4rlyPXAXsAL4EOhyhB8h8nMnAO9FvP4SP8ylJSVz8WGOfsaSj/F/gITtNLOHgJfxsYer5+ifqzyHKbs2exN42szq4/9b7Mf3fksS0RhrSXaF4/4AnHOfAvcCPcrYPw2fxFYCm4BLI980s0Gh7eebWYOo7btCL3seZaxzgMeBF0OJr0LMrC0RN8KYWTv8jCjd8DczTijjoyIilfFd/DdkL+Pz46mhnmKgMD9uBQ4C50Rs72NmjYnKtWZWJ9RREa2snFyaZ4DVzk+rmleBn2Utvoe78BTAN/ihLnHjnCsIdX70Cv1h8gbwgHMuJ57niTrnZGAy/pvVK4HhzrlvEnU+OToqrCUV3BoqOsPa4cfGhR1rZmmh5yPw4+t2Av8CLjKzP5vZiFBvwl58AmwETDGz0aEbRXrgi/GFwK/NLM3MMvBjm1uaWW0zu9vMri4jRgu1+/G/fP5UynvR+0e/HoMfDxh2AnChc24XcGspnxERqYgSPaFmdjrwK+AHzrlcfO/wFuANMxsbukFxWKiH+t/42UMuNLMzge8653bjc+2PzOxXZjYc38EQHkYXCBA3+gAAIABJREFUmf/Kyslhkb3HJwMtzKwhcCrQwMyOC713CH+zX7iDJfIczwCXmFmz0OtBwLPOucP4mqfUXGxmg0Nx1Svj2tUiqmYK/fFxPv5bymb4IYbdgQ5Aw4gbI6N7xQ+WEn/xoEoZLxL6b3EM8C6+V74if3BIFdFQEEkFHYDFZvY6/msyA26MeP8w8DMz24nv3b0VwDk3zcx+hP868UL8zX+fA5jZOOBJ/C+K54E7nXPOzL6FnwppMTAR/5VeP/yCCt/CDzEptjCNmQ3E35TTEf/V4H+BX4XimQoMBjqbWVd8T0pr4FtmtpSi3ptrgUbOueyon/1v5hd2OIy/WUhEpMJCeeQqfPH3LJCLvwGuCXCuc24ugHPuYCg/PovPda/hx1YD3IbvlHgB/w1deCaOv+Jz3034IWs3Oee+DhWOg4GuZja+rJxsZpfixxVfb2ZLnXMb8DdS/gFfVL8AjAP6AF8B/8T3rn/H/GxRvYDzzOwN59yrZnYCMMnMPsB/6/mzUM/6ZUDr0M+3GL9YThsz6xWK8wr8uPH/Rly3uvjc3xdoFyq8HX5s88WhY34CzADm44e3vInP6VPwN0z2NbMLnXPh407EL9azJ1SItwYuNbNV+J5oQp+PXrNhLf7mzAuB/8/encfXVdf5H399sjZtlqZNum90oZRiS2kpu5SlA46ogDAKo4gKVWdk0HFBEVABR2XQURFH64i7iIgDyiiL/Ni3rrR0BbrTNV2ytk2a5PP743vT3KZpmjTnbsn7+Xicx80959xzPjm999NPvvd7vt8+QJaZ/Z1wA2OyhqCVozD9W0gmszAN7C/cfUySzvcFd9d04yIiPUysW0q+u6flBF1m9gXCCC3rY8/7EL5t+HbLvTmSeuoKIpmuvW4WiTlR+Mp0UTLOJSIiyRO7OfOUNC6q+wC3AGPiVpcBdSqq04u6gkjGivVv+xDhq70r3f3BBJ9yTewrShER6Vly3P2/Ux3Ekbj7fgvTqv8gdmP7FkK3lVs6fqUkm7qCiIiIiIhEQF1BREREREQioMJaRERERCQCKqxFRERERCLQY25eLCsr8zFjxnTpNRUVFZSXlycmoARRzImXafFC5sWcafFCYmNeuHDhTnfPrAvSTceSsyEz3zvpTtc0Wrqe0Uu3a9pRzu4xhfWYMWNYsKBrM4nOmDGjy69JNcWceJkWL2RezJkWLyQ2ZjPbkJADp7FjydmQme+ddKdrGi1dz+il2zXtKGerK4iIiIiISARUWIuIiIiIRKBXF9Zz5sxJdQhdppgTL9PihcyLOdPihcyMuSfSv0P0dE2jpesZvUy6pkmfIMbMZgN3Ah+Im+/+cuA7QDHwG+Bz7t4Y21YGfAOoBHKBm9z9QNvjzpgxw9Op/42ISGeZ2UJ3n5HqONqjnC0icqiOcnZSW6zNbDBQCMyMWzcKuBS4AvgscB3wmbiXPQjc5+43AW8AtyctYBGRXkw5W0Ska5JaWLv7duCRNqtHAx9z94Xu/ivgXuA8ADM7HZjo7q/G9v0LcIOZ9UtWzCIivZVytohI1yS9j7W7N7d5/nzLV4gxm4GNsZ/PB9bF7bs59mP0X5nedBNMnAgvvRT5oUVEMlW65OzqBdUsnLmQVR9d1d1DiYgkTDrevHgq8OPYz8OBXW221wDD2r6ooqKCGTNmHFzmzp3btbNu3QpvvAHLlx9DyCIiXTN37tyD+QooS3U83ZCcnN0ENfNrqFtWF0HIIiJd09mcnVYTxJjZeKDC3ZfEVjmwr81u+cBhN8KUl5d3b/DwSZPC48qVx34MEZFOmjNnzsE73c1sZ4rDOSbJzNk5peG/qwN7DjuUiEjCdTZnp02LtZllA58AvhS3ejNQGrePAf2ALZEHoMJaRKTTkp2zWwrrxsrGo+wpIpI6aVNYA58H7nb3+rh1TwIT456PARqA6MdoUmEtItIVSc3ZOf1bC+tkDxMrItJZSS+sYy0YABa37mZgMVBgZmPN7GNmNt7dFwBbzeyk2K4XA/e4e0PkgY0dC7m5sGED1KkPn4gIpE/OzsrNIqtfFjRBU01Tdw8nIpIQSe1jbWZFwIdiT681sx8Svkq8o82uq9z9vtjPHwBuM7P1hPFUv5KQ4HJzYcIEWLEiLKeempDTiIhkinTL2bmludTX1dO4p5Gc4rS6RUhEBEhyYe3uNcB/x5YWd8aWI71mA/DxBIcWzJgRiur581VYi0ivl245O6c0h/q36zmw5wB9RvdJxClERLolnfpYp97M2ORi8+alNg4RETmMbmAUkXSnwjpeSyu1CmsRkbRzsLDeo8JaRNKTCut4U6eGvtarVkFlZaqjERGROLmluQA07lZhLSLpSYV1vPx8OO00cIdnnkl1NCIiEid3cCisG7ZHPzCUiEgUVFi3NXt2eHziidTGISIih8gbkgdAw1YV1iKSnlRYt9VSWD/5ZGrjEBGRQ+QPzQegYZsKaxFJTyqs2zr1VOjfH956Kwy9JyIiaaGlxbp+a/1R9hQRSQ0V1m3l5MDll4ef778/tbGIiMhBB7uCqMVaRNKUCuv2XH11eLz//nAjo4iIpFze0NY+1q7cLCJpSIV1e2bNgiFDYM0aeP75VEcjIiJAdlE2WQVZNO9tpqmmKdXhiIgcRoV1e7Kz4brrws/33pvaWEREBAAzOziV+f71+1McjYjI4VRYH8knPhEK7D/9CbZsSXU0IiIC9BkbCut9a/elOBIRkcOpsD6SESPgssugsRG+851URyMiIkDBuAIA9q9Ri7WIpB8V1h358pfD449+pFZrEZE0oBZrEUlnKqw7csop8P73w/798PWvpzoaEZFer2BsrMV6rVqsRST9qLA+mttvD2Nb//SnMG9eqqMREenVCsaHwnrv6r0pjkRE5HAqrI/mxBPh3/89jGc9Z07ocy0iIilRMKEAyzP2r9tPY7XysYikFxXWnXHbbTBmDCxZoi4hIiIplJWbRb8T+wFQt6wuxdGIiBxKhXVn9OsHv/gFmME3vgHPPJPqiEREeq1+U0NhXbukNsWRiIgcSoV1Z517LtxyS+gS8s//DFu3pjoiEZFeqXBKIaDCWkTSjwrrrrjtNjjnnDD03qWXwj4N9yQikmxF04sAqH6lOsWRiIgcSoV1V+TkwB//GPpbz5sHH/kINDWlOioRkV6laGYRlmvULa3jQOWBVIcjInJQ0gtrM5ttZq+a2Zi4dQVmdo+Z/YeZ/beZFXdmW0oMGgR/+QsUFcGDD8LHPw7NzSkNSUQkUdIxZ2cXZFM0owgcql9Wq7WIpI+kFtZmNhgoBGa22fQj4AV3vxl4BPhJJ7elxkknwf/9H/TtC7/8JVx/vVquRaTHSeecXXJOCQBVz1Ul4vAiIsckqYW1u28nJNqDzGwYcDXwl9iqvwPvNbMxHW1LSsAdOeecUFwXFMB998EVV6jPtYj0KOmcs/uf1x+A3Y/vjvrQIiLHLOldQdy9bb+JWcBOd98b294IrAPOPcq21Js1C554Avr3h4cfhtmzYbeSvIj0HOmas/vP6k9W3yxqF9ey/21Nby4i6SEdbl4cDuxqs64GGHaUbYeoqKhgxowZB5e5c+cmJNjDnH02vPACjBgBL74IZ54Jq1Yl59wikrHmzp17MF8BZamOpwvSImdn98mmdHYpALv/Tw0aIpJYnc3ZOckL6YgcaNuHIh9oudW7o20HlZeXs2DBguij64zJk+Gll+Dd74bXX4eZM+HXv4b3vS818YhI2pszZw5z5swBwMx2pjicrkibnF32njJ2PbKLiocqGPaJw2p3EZHIdDZnp0OL9WagtM26ImDLUball5EjQ3H9T/8ENTVhnOubb4YDGgpKRHqUtMnZZZeVYXnGnqf2UL+5PhGnEBHpknQorJ8GhplZAYCZ5QKjgGePsi39FBbC738P3/42ZGXBN78ZuoasXp3qyEREopI2OTt3QC4D3zMQmmH7b7Yn4hQiIl2SinGsreVHAHffBjwKXBBbfyHwkLtv7mhbEkPuGjP44hfhmWdg9GhYsACmTYN77tGQfCKScdI9Zw/5yBAAtv58K97siTqNiEinJHsc6yLgk7Gn15pZeeznTwHvM7OvAP8IzIl7WUfb0tc558CSJfDhD4dh+P7t3+CMM2Dx4lRHJiLSKZmQswdcPID8EfnsW72P3Y/pJkYRSS1z7xl/4c+YMcNTdvPi0fzv/8INN8DmzaGLyL/9G3z1q2GYPhHp9cxsobvPSHUcyRRlzt74nxtZ+8W19L+gPyf//eRIjikiciQd5ex06GPd8112GaxcCTfeGJ5/73swfnzoHqKbG0VEumXo9UPJLsym8qlKqudrinMRSR0V1slSVBQK6vnz4Z3vhF27Qsv15MmhRbuHfHMgIpJsuf1zGfYvYbi9tV9em+JoRKQ3U2GdbKecEm5sfPhhOP54ePNNuPxymD4dHnlEBbaIyDEYddMocvrnUPlUJbufVF9rEUkNFdapYBYmj1m2DH74Qxg6NNzUeOmlofB++GFobjuLsIiIHEnugFxGfWkUAG999i2aG5RDRST5VFinUm4u/Ou/wpo18P3vhwL7tddCn+yTToL/+Z8wooiIiBzV8H8bTsH4AvYu38umuzelOhwR6YVUWKeDgoLQ33rNGvjBD2D48HCz4/XXh7Gwv/Y12LEj1VGKiKS17IJsjv/x8QCsv309e9/Ym+KIRKS3UWGdTgoKwrB869bBb38buoVUVMDXvw6jRoUxsV98Uf2wRUSOoPSCUgZfMxivd1ZctYLmenUJEZHkUWGdjnJz4eqrw6yNzz4L730vNDTAb34DZ58N73hH6JtdVZXqSEVE0s6EH0ygz3F9qF1Uy5qb1qQ6HBHpRVRYpzOzMDTfI4/A2rVw880weDAsXx5atocNg+uug3nz1IotIhKTU5LDib8/EcsxNn9/M9vv357qkESkl1BhnSnGjIFvfAM2boQ//AHOPx/27oWf/QxOOy2Mh33XXbBlS6ojFRFJueKZxYz7zjgAVn10FVUv6Rs+EUk8FdaZJi8PrrwSnnoKVq+Gz30OBg0KNzvedBOMHAnvehc88ADs35/qaEVEUmb4DcMZ9qlheL2z7NJl7H1TNzOKSGKpsM5kxx8Pd98Nb78Nf/5zmGgmOxseeww++EEYMgQ+/nF44glNnS4ivY6ZMf4H4ym9qJQDFQdYcv4S9q3TEKYikjgqrHuC3Fx4z3vgoYdg61a4554wk2NVFdx3H1x0URgje86c0NLd2JjqiEVEkiIrJ4vJf5xM8VnF1L9dz5Lzl7B/o77NE5HEUGHd0wwcCJ/+dBhRZPlyuO02mDgRdu2Cn/4ULrww3PT4yU+Glu36+lRHLCKSUDmFOUz56xSKZhaxf/1+Fp+9mLqVdakOS0R6IBXWPdmJJ4YxsFeuhKVL4ZZbYMKEMDb2T34S+mKXlYU+27/+dSi+RUR6oJziHKY8NoXiM4up31TP4rMXU/WybmgUkWipsO4NzMLY13fcEW54XLwYbr0Vpk6F2lr44x/hmmvCUH6zZsF3vwtvvKEh/ESkR8ktzWXqk1MZeMlAGnc3suSCJWz/vYbiE5HoRFJYm9k/RnEcSQIzOPlkuP12eO01WL8+TKN+4YVh27PPhpFGJk6EcePgU5+Chx+G6upURy4iEenNOTu7bzaT/3cyQ68bSvO+ZlZetZI1X1hDc6NmaBSR7jPvQqukmV0M3AIMobUozwIGu3tB9OF13owZM3zBggWpDCHzVVXB3/4WRhh5/HHYvbt1W04OnHlmuBHyootg2jTI0hceIlEws4XuPiMBx1XOPgJ3Z/O9m1nz2TV4o9P//P5M+vUk8oflpywmEckMHeXsrhbWFcA3gcVAy5/3Blzt7nO6G2h3pDpJ9zhNTbBwYbjB8fHH4ZVXoDmuRae8HM47r3U5/vjQ4i0iXZbAwlo5+ygqn6tk+ZXLObDjADkDc5j404mUX1ae6rBEJI1FWVi/6u6ntbN+iLtv60aM3ZYuSbrHqqwMQ/U9/ngotjdtOnT70KGHFtpjx6rQFumkBBbWytmdUL+1nlUfXcWex/cAMPS6oYy7exw5JTkpjkxE0lGUhfUpwFh3/2Ob9XPcfW73wuyedErSPZ57uLnx6adbl4qKQ/cZOTLcCDlrFpx1llq0RTqQwMJaObuTvNnZ/MPNrPniGrzeyRuWx4QfTlDrtYgcJsrC+gbgDqAofjXg7p7dzSAnAZ8G3gImAHPd/TUzKwDuAmqAUuAmdz/sTrp0S9K9ijusWNFaZD/zzKH9syEM63fmmXD22aHQnj4d8tWXUQQSWlgrZ3dR3fI6Vl+3mupXQsgD3zeQCfdMoM/IPimOTETSRZSF9R7gE8ACoKllNfBhd7+jm0EuAN7n7pvNbBTwuLtPMrOfA4+5+wOxG3E+4u5XtX19uibpXqm5GV5/PRTZL7wQlu1thrTKz4cZM1oL7TPPDJPbiPRCCSyslbOPgTc7W368hbVfWktTTRNZBVmM/PxIRn5xJDmF6h4i0ttFWVj/Dviou9e3Wd/t/npmVgdMd/dVZlYOLAFmAOuAUnffa2Y5QBUw2d3Xx78+nZN0r+cOa9fCiy+2LsuXH77f+PEwc2brcvLJUJDSgQtEkiKBhbVydjfUb67nrc+9RcUDoatb3pA8jrvzOIZcOwTLVtc2kd4qysL6vcApwP+LXw1c7u43djPIXwLjgHcB1wEbgXzgP919eNx+y2Lrfhn/+kxI0hJn9254+eXWQnvePNi//9B9cnLCxDbxxfakSZDdrW+wRdJOAgtr5ewIVL1UxVuffYuaeTUA9J3Ul9G3jWbQlYNUYIv0QlEW1vMJfekqgZYXGjDE3bvVAc3M+gGPACOA77j7T83sC4SvLKfE7fcy8Gd3/2b86zMpSUs7DhwIrdjz5rUuy5cfOsQfQL9+oX/29OlwyilhPO2JE0MRLpKhElhYK2dHxJudHQ/sYN3N69i/PjQC9D2xL2O+OobyK8qxLBXYIr1FRzm7q9XIfwJ/c/eaNieIYhavAuB+YDjwPTNbR/iPYF+b/fKBA21fXFFRwYwZrb/jnDlzmDMnpcO0Slfk5oauHyefDC3/brW1sGgRzJ/fWmyvXw/PPReWFn36wJQpocieNi0U3O94R1gvkqbmzp3L3LkHB+YoS9BplLMjYlnG4KsGU/7+crb9chsb7tzA3hV7WfGBFfQ9sS8jPzeSwf88mKx8TZwl0hN1Nmd3tcX6DHd/uZ31H3P3+7oc5aHHeAz4oLtXmtmtwL8BnwG+6u7Hx+33Zmzd7+Jfn2mtH3KMduwIhfaiRbB4cVjWrz98v+zs0G2kpdieNi0U3wMGJD1kkaNJYIu1cnaCNDc0s+3n29jwjQ3Ubwpd2POG5DH8huEM++QwcgfkpjhCEUmUKLuC/NrdP9xm3ZXAz9y9uBsBlgGvu/vQ2HMjDOH0EeAxoNzd95lZLlBLGJd1c/wxMj1JSzfs2QOvvdZaaC9aBKtWHd6NBGDYsNCaHb9MmqTWbUmpBBbWytkJ1tzQzI4/7GDT3ZuoW1IHQFbfLIZcM4RhnxxG4dTCFEcoIlGLsivIADO71d3vMLO+wA+AK4H53YxxN7DfzIbHJd9dwFLgUeCC2OOFwENtE7T0cqWlrTM+tti7Nwz511Jsv/YaLFsGW7aE5fHHW/fNzg4T2LQtuMeMgSx9rSsZTTk7wbLyshjyoSEM/ufBVP6/SjbdvYndj+1my4+3sOXHWyg+vZhhnxxG+T+Vk12gG69FeroutVgDmNnpwD8DFwHbgGuAKnff061AzKYC/wIsBAYDz7n7s2ZWSphsYD0wBPiyu9e2fX1Paf2QBGpuhnXrYOnSUHS3LG++2X7rdmEhTJ4ciuzJk+HEE8MyfLhmkZRIJarFOnZs5ewkq1tRx5afbGHbL7fRVBWGD8/pn8PgawYz5NohFJ5ciCmHiGSsY+4KYmbvPMKmmcB5wN1AM3Cpu3+2u4F2R09O0pJg+/bBypWhyI4vurcdYZjfoqLQfaSl0G5ZRo9WC7cck6gKa+Xs9NK0t4kdD+xgy4+3HByqD6Dv5L4MuSa0cucP1wy0IpmmO4X1K8AkwlBN7TTphd2Aoe6e0uzQG5K0JNnOna1F9sqVYdr25cth16729y8oaL/gPu44DQcoHYqwsFbOTlM1i2rY9vNtbL9/O427GsNKg9ILShn84cGUva+MnBLlCZFM0J3C+h+BBe6+4ygnuMLd/9i9MLuntyVpSaGKilBkr1jRWnCvWAFbt7a/f35+6MPdUmifcEJYJkzQzJICRFpYK2enueaGZnY/tpttv9rGrr/swhvC/8GWa5T+QynlV5RT9r4ycks1qohIuopsVJB01luTtKSRPXsOLbRblk2b2t/fLHQfaSm0J05s/XnwYPXj7kUS2cc6XSlnw4E9B6j4QwU7fr+DyucqD37HYDlG6YWllF8ZK7IHqsgWSScqrEVSqbo6DP+3cmXoSrJ6dXi+Zg00NbX/mpKSQwvtlsJ7/HjIy0tu/JJwKqylYXsDFf9bQcUfK6h8urXIJgtKziph4HsGMvCSgfQ9oa9ufBRJMRXWIumooQHWrg1F9qpVrQX3qlVQWdn+a7KzYezY9ovuskRN3ieJpsJa4jVUNLDz4Z1UPBiKbG9s/X+6z7g+DLxkIGXvKaPknBKy8nTDtEiyqbAWySTuoR93ewX3unVhe3sGDjy8S8nEiaEQ182TaU2FtRxJY1Ujux/fza5Hd7Hrr7tab3wEsouzKZ1dyoB/GEDp7FIKjtM9GyLJkNDC2symAZ8DFhMmAljfrQMeIyVp6RX274e33mq/6K49bKjgIDcXxo1rLbgnTmz9WVO8p4VkFtbK2ZnLm5zqV6rZ+Zed7Hp0F3uX7z1ke8H4Akr/IRTa/c/rT06x/qAWSYREF9Y/Ad7n7kPM7Grgf919X7cOegyUpKVXcw8zSsYX2i2F98aNR35dWdnhxfbEiWGIwFzdMJUsSS6slbN7iH3r9rHniT3sfmI3e57ac3AyGgCyofj04lBkn9+f4lOLycpXtxGRKCS6sM4FBrn7ZjPLAgrdvbpbBz0GStIiR1BXF2aXbCm6V69uXerq2n9NTk64UbJtwT1xYuhyIpFKcmGtnN0DNTc2U7Og5mChXf1KNcTV2VkFWRSfUUz/Wf3pf25/ik9ToS1yrCItrM3sXcAAd/+tmc0Eat19RQRxdouStEgXucPmzYcW3J1t5W6vW4lauY9Zgqc0V87uhRqrGql8ppLdT+6m6tkq6pYd+kd0Vp9YoX1uf/rP6k/RaUVk98lOUbQimSWywtrMvgHcCDzm7lfE1n0LWOruv4si2GOlJC0SofhW7rYt3R21co8bd+gNlC2T4hQVJTf+DJOowlo5W1o0VDRQ9VwVlc9WUvlMJXWvH/o5tnyj+NRiis8qpuSsEorPKCavTEN7irSno5zd1TsbzgRGAB+JW/c34NdASpO0iESoXz84+eSwxGuvlbvl540bW5+3NXo0TJ4clpNOCo+TJkHfvsn5fXov5WwBIK88j/L3l1P+/nIAGnY2UPV8FZXPxArtpXVUvVBF1QtVbCJMalUwsYCSs0pCoX1mMX0nagxtkaPpamH9pLtXmll8M/cVUQYkImnMDEaMCMsFFxy6be/e0MrdUnC3zEK5ahVs2BCWv/710GONHdtacLcU3RMnQp8+yf29ei7lbGlXXlke5ZeVU35ZKLQP7D5A9cvVVL1YRdWLVdTMq2Hf6n3sW72PbfdtAyBnYA4lZ4Yiu+SsEoqmF5HdV91HROJ1tSvIB4GxwHDgeeBK4DLgBne/NyERdpK+VhRJU42NYYjA5cth2bLwuHw5vPFG2NZWVla4cbKlZXvKFJg6NXQzyeqZN1slsCuIcrYck+aGZmpfq6XqxSqqXwoFd8PWhkN3yobCKYUUzSyi+LRiimcWh5khs9WqLT1b1Dcvng5cC4wCtgL3u/vfuxtkdylJi2SYhoZQXMcX28uWhanem5sP379fv1BgT53a2k3lpJN6RHeSBN+8qJwt3ebu7F+//2CRXfVi7IbINh/V7KJsimaEQrul4M4flp+aoEUSJNHD7WUDs9z9qW4dqJuUpEV6iH37QleSZcvCsnQpLFkSxuluKysLjj/+0GL75JNhyJDkx90NSR5uTzlbItFU10TNwhqq51VT82p4rN9Yf9h+ecPzDrZoF51WRNH0InKKNHmNZK5jLqzN7M9A4VGOXwbsdPfzjz3E7lOSFunhKipCgf3aa62PK1dCU9Ph+w4bBjNnwqmnhscZM6B//+TH3ElRFdbK2ZJq9VvrqZlfQ/Wr1dTMC8V2U3Wbz6hB34l9KZxeSNH0UGgXTitUsS0ZozujgiyN7bOa8IXP9YRpcOOz4XGA5kUWkcQqL4cLLwxLi/37ww2S8cX2a6+F1u2HHw5Li+OPD4V2S7E9bVpPvElSOVtSKn9oPvnvzafsvWUAeLOz9429B1u0q1+tpm5JHXtX7WXvqr3s+O2O8EIV29JDHK3FuiS2T2Xs+e3ufls7+z3i7u9LXJhHp9YPEQFC/+w334T582HevPC4eDHUt/mKOi8vFNnvfCeccw6ceSaUlKQk5AhbrJWzJe011zdT+3ottQtrqVlYQ82CGuqW1eEH2tQjKrYlTR1zi7W7V7VZNdjMst394Pc6ZjYLmNntKEVEopCV1Toj5Ic+FNY1NIT+2i2F9quvhpbuF18Myze/GV43dWooss87D84/H4qLU/u7dJFytmSCrPwsimcUUzyj9fN1WLG9sIa619tv2S44voCiGUWtxfbJheQUq9iW9NDV4fZmA/8FvADUAJOAi4Db3P2bCYmwk9T6ISJdsmdPKKqffz4sCxbAgQOt23Ny4Oyz4eKLwzJlShh7OwESONyecrZkrCMV24e1bAMFEwpCkX1KIUWnhMfc0twURC29QdTD7Q0D/gWYAlQDD7r7I92OspuUpEWkW/buDS3Zzz0HTz4JL78tcL0NAAAgAElEQVR86LB/o0bBBz4AH/xg6J8dYZGd4OH2lLOlx2iub6ZuWd3BLiQ1i2LFdsPhtUyf4/qEQnt60cFiO69c07RL9yV6uL0sYJq7L+zWgVqPdzpwNvAWoZWlDriL0NpSCtzk7tVtX6ckLSKR2rMH/v53eOwx+NvfYOvW1m3HHw+f/CR87GOR9MtO8nB7ytnSozQ3NFO3oi60bC+qoXZRLbWv1dK8//Dx8PNH5Ic+27FCu2h6EflDNc62dE1khbWZ3dfO6nKgyd0vPcb44o//UWCcu98St+7nwGPu/oCZXQx8xN2vavtaJWkRSZjmZnjpJfj97+EPfwhD/wEUFcGnPgVf+Uq3+mMnsCuIcrb0Ss2NzexdtZfaRbFie2EtNYtraK47vNjOG5J3sMhu6UqSPzIfS1DXL8l8URbWzwNPAfEvOgt43t3v6GaQ5wK3ArM9FlTsK8x1QKm77zWzHKAKmOzu6+NfryQtIknR2AiPPgrf/z4880xYN3gw/OIXoS/2MUhgYa2cLRLjTc6+t/aFbiSxlu2aRTU0VR0+Fn5uWW5rf+1YC3ef4/qo2Bage+NYt3Wdu69uc/ATCNPldtd3gZeAe8xsPHAHMJowkcFeAHdvNLN1wLnA+gjOKSLSNTk5cOmlYZk3Dz772dCa/e53wwMPwBVXpDrCeMrZIjGWbfSd2Je+E/sy+OrBQBhne/+6/dQsCjdH1i4KN0oe2HmAPU/sYc8Tew6+Pqd/ziE3RxZNL6JgfAGWpWJbWkXRx/o4YL67l3XjGBOBlcDJ7r7UzL4A3AbcCfyzu0+J2/dl4M9t72gfPXq0l5eXH3w+Z84c5syZc6whiYh0TnMz3Hor/Md/hNkdt2yBgoKjvmzu3LnMnTsXgIULF25w9zEJjhRQzhY5GnenfmN9a6t2bESSAzsOHLZvdmE2hdMO7UZSMLGArJysFEQuidTZnN3VriDrOPQrxWxgMPBrd7/+WIM1s8uBuS2J3swKgArgduD97n5a3L6LgN+5+93xx9DXiiKSMu5h1JC334aFC+GUU7r08gR2BVHOFomAu9OwtaG1VTvWwt2wueGwfbMKsig8ufBgoV00vYi+J/YlK1fFdk8RZVeQJ4Hf0Zqom4Ht7v5GN+JrieNgLO6+z8zeJPwnUNpm3yJgSzfPJyISDXe4995QVBcUwNixqY4onnK2SATMjPxh+eQPy6fsPa1f9jRsb6BmcWvLdu2iWvav30/1y9VUv9w6GI7lG4XvKDxkRJLCdxSSla9iu6fpamH9G3d/Ln6FmU00s/Hu/lY34lgKlJjZYHffHlvXSEjGw8ysIJa4c4FRwLPdOJeISDRefx1uuikMxwdw112hO0j6UM4WSaC8wXkMvHggAy8eeHDdgd0HDrk5snZhbbhpckEYe3srYehOyzH6ndTvkBFJCqcUkt03O1W/jkSgq4X1xcBzbdZtB54Gph1rEO6+ysz+BlwB3Gtm/YE84LfAu4ALgEeBC4GH3H3zsZ5LRKRb6uvD2Nb33hsmkoEw7N6PftQ6hXr6UM4WSbLcAbkMuHAAAy4ccHBdY1Ujta/VHjIiyd5Ve6l9LYy5ve2+bWHHLOh3Yr9DRiQpnFpITpGmbM8UnfqXMrNPE25KKTKzm9rZ5ckIYrkG+L6Z5RFaOK6K3VH+KeAuM5sKDAF0d4uIJNf+/fDss/Dgg/DQQ1BZGdb37Qsf/3i4eTHuRrxUU84WSS85JTn0P7c//c9t/UarsbaRuiV1h4xIUreijrplYdn+q9iXQQYFxxccMoNk4bRCcvtryvZ01OmbF81sCvAR4AdtNjW4+9Z2XpJUuhFGRCLjDm+9BY8/Hrp5PP007NvXun3KFLjmmjDzYmnbLsVdl4ibF5WzRTJP074m6pbWHTIiSd2yOvxAO1O2j+tD0fQiik8vpviMYoqmFanPdpJEcvNibEilL7n7YePNmNkod9/YnSBFRFKmuTn0l37++bA89xxs23boPtOmwSWXwAc/CCeemJo4u0A5WyTzZBdkU3xaMcWntc7k2lzfTN3yutYZJBfVULuklv1r9rN/zX4q/hBmgrV8o+iUIorPCIV2yRkl5A/XdO3J1mFhbWafB/7q7itiq75oZm171RtwHjAr+vBERBJg82ZYsADmzw/LvHmt3TtalJfDeefBu94FF10EQ4emJtYuUM4W6Xmy8rPCsH2nFMF1YV3zgWb2rtxLzfwaql6uovrlavau2HvYaCT5o/IPFtnFZxRTeHIhWXlq1U6ko7VYHwfE3+J+CtAP2Ebr8E1G6EcnIpJe3MMweK+/DosWtRbSW9vpCTF6NJxzDrzzneFx4kTIvOmLlbNFeoGs3CwKp4RRRIZ+PPzRf2DPAapfrT5YXFe/Wk39xnoqNlZQ8UBo1c7qk0XRjCJKzimh/7n9KT6zWDdGRqyrE8ScAixz94Y260+MayFJCfXXE+nlamth2TJYujQU0kuXhqVtSzRASQnMmAGnnhoeZ86EkSOTH3NMAieIUc4W6aW8yalbWReK7JeqqXq5in2r9x26UzYUnVJEyTtL6P/O/pScU0JuqW6KPJooJ4h5h7svanPwicA4IKVJWkR6AXeoqIBVq1qX1ath5UpYt6791wwYAFOnhhsOTz01LOPHQ1av+DpUOVukl7Jso/CkQgpPKmTY9cOAMMZ21UtVVD1XReVzlWFs7flhefs7b4NBv3f0o/+5/UOxfW5/8srzUvybZJbODrd3DuHrw9mxKXLjDSDcdf6XiGMTkd6qrg7Wrg0jc7zxxqGFdHst0AC5uTBpUiigW5Z3vCP0jc68Lh3dopwtIu3JHZBL2SVllF0SZo9srG2k+uVqKp+tpOq5KqpfraZuaR11S+vYfE8Yfr7f1H4MmD2A0tmllJxTQnaBJrDpSGdbrCuBXxH6753VZlsD8D9RBiUivcDu3bBmTSie4x/XrGm/D3SL4uJQQE+cCCecEJaJE2HChFBcCyhni0gn5BTmMGD2AAbMDpPZNO1romZeDZXPVVL5bCXVL1ZTt6SOuiV1bLp7E5ZvlJxdcrDQLjy5EMvqXQ0XR9OpwtrdXzezWcCV7q6ELCJH19QEW7aELhotBXN8Ab1nz5Ffm5sLY8fCuHGhYG4poE84AQYP7nUt0F2lnC0ixyK7ILt1IptboWl/E1UvVLHn73vY8+QeahfXUvlUJZVPVcKXILcsl9LZpQy8ZCADLh5A7gA1bnRlHOsqjtDKYWZnu/sLkUUlIunPHXbuDIVze8uGDXDgsCGUW/XrFwrn8eMPfxwxArL1dWN3KGeLSHdl98lunZ79W9Cws4HKpyrZ8/c97H5yN/Ub6tlx/w523L8DsqHkrBIGXjKQgZcMpO8JfbFe2AhytHGs/wwUHuUYBkwGBkUVlIikiZqaIxfO69aFvtAdGTQIjjsuFMxti+dBg9TyHDHlbBFJpLyyPAZ9YBCDPjAId2ffW/vY/dfd7PzLTqqeDTdFVj1XxdovrqXPuD4MvGQg5ZeVU3J2CZbdO/L90VqsVxP66m2mdQzUtgy4JMqgRCRJ6upg48awtFc479rV8etLSkLhPGZMeIxfxowJrdKSTMrZIpIUZkbfCX3pe2NfRtw4gsaqRnY/sZtdf9nFrr/uYv+a/Wz+/mY2f38zeUPyKHt/GYOuHNTji+wOx7E2szKgxt3rOzyI2UB3P8r/wImlMVFF2mhuhu3bQ5eMluI5ftmwIdxA2JE+fdovmluW0tKk/Co9XVTjWCtni0g68Can+pVqdj6yk4o/VrB/3f6D2/KG5FF2eRmDPhArsjPw5sdjHsfa3Xd28hyXobvMRZKrrg42bTq0UI4vnDdt6riPM0BeHowaFZb2CujBg3vLeM89gnK2iKQDyzZKziqh5KwSxn57LLWLatnx4A4qHqxg/9r9bPnRFrb8aAt9xvRh8EcGM+SaIRSMLUh12JHo0syLh73YbBTwr8Cn3T2l3/mq9UN6lKamMORcS4Ec/9iyHK2bBkB5eWvh3LKMHt36c3m5Cuc0kKiZF9s5j3K2iKSMu1O7uJYdf9jBjt/toH5T65drJeeWMOTaIZRfUU5OYXpPsx7lzIuYWTbwXuB6YDawndCnT0Q6wz0MNXekonnTJti8ORTXHcnLC9Nwty2WW5aRI6Fv3+T8TpK2lLNFJF2YGUWnFFF0ShFj/2MslU9Xsu0X26h4qCLc/PhsFW995i2GfnQow/5lGH0nZN7/YZ0urM3sOOA64KPAQGAD8B7gceCMhEQnkon27Tu8YG77uHfv0Y8zeHBr4Rz/2FJEDxqk1mY5IuVsEUlnlmWUXlBK6QWlTLh3AhUPVrD1Z1upfrmat7/3Nm9/720GvGsAwz89nAEXD8iYvthHLazN7EpCS8f5wFvAd4FfAje7+2Ox3V5KWIQi6eRoXTQ2bQpjOx9NUdHhBXP844gRkJ+f+N9HehzlbBHJNDnFOQz9+FCGfnwoNYtr2PzDzez43Q52/203u/+2m34n9WPUzaMY9E+D0n5Ekc60WO8jTIH7EnCdu78BYGbH3jlbJB3Fd9E4UmtzZ7po5OaGwrilSG6vcC4pSc7vJL2RcraIZKyiaUWc8LMTGHfXOLb+bCub79lM3bI6Vl69kvW3rWfUl0Yx+JrBZOWm5ze2Ry2s3f1R4NHYTS/Xm9kY4BHg4LyVZmbenbsgRZIhyi4aHbU2ayQNSSHlbBHpCXIH5jLqi6MY8ZkRbPvVNjZ+cyP73trH6utWs/GujYz99ljK3leWdrM7dmVK843ArWaWA1wKjDCzO4CfAf8AzE1MiCKd0F4XjbaFc1e6aBypcFYXDckQytki0hNk5WUx7LphDLl2CBUPVLD+6+vZ98Y+ll+2nJJzShj33XEUzyhOdZgHdXe4vfHADcD17p7SWzc1dFMP5h4mMumotbkrXTTaK5pbflYXDUmBJA63p5wtIhmt+UAzW36yhfVfW0/jrkYwGH7DcI77xnFJG6Yv0uH24rn7W8CNZrawO8eRXq6jLhotP3emi8aQIUfunjFypLpoSK+nnC0imS4rN4sRnx7BkA8PYcOdG9j0X5vY/IPN7Hx4JxPnTmTARQNSGl8kpb27/yqK45hZHjAfuNHdnzGzAuAuoAYoBW5y9+ooziVJ4h6m1V63LswM2F4BfaxdNOJ/Hj5cXTREOimqnA3K2yKSGjklOYz7z3EMunoQq69fTe3CWpZevJSRXxjJcd84LmU3N6bb1DZfAMbEPf8R8Ji7P2BmFwM/Aa5KRWDSgerqUDjHL2vXtv68b1/Hr8/NDQVyR63N6qIhkq6Ut0UkZYqmFXHKK6ew6a5NrLttHZv+cxNVL1Qx+aHJ5A9NfoNb2hTWZnYmsBXYE3s+DLiaMP0uwN+Bh8xsjLuvT0mQvVVDQ2hdblswt/x8tKm1BwyA444LE5uMHn144ayJTkQykvK2iKSDrJwsRt88mpJ3lrDyqpVUv1zNojMWMeVvU+g3qV9SY0mLwtrM+gFXuPu/m9ltsdWzgJ3uvhfA3RvNbB1wLrA+JYH2ZA0NoUh+883Dl02boLn5yK/t0ycUzi3L2LGHPldrs0iPo7wtIumm/9n9mb5oOsveu4zqV6pZfOZipjw+heKZyRs1JC0Ka+Am4Ftt1g0H2jaF1gDD2jtARUUFM2a03qA5Z84c5syZE2WMma+xMfRzbls4v/FGWH+kUTWyskLLctuiueVn3RQo0mVz585l7tyDI96VpTKWY9StvK2cLSKJkFeex9SnprLi6hXsemQXSy9eysnPnEzhlMJuHbezObtbw+1FwcwuAvLd/c+x5+uBa4EZwJXuflrcvouA37n73W2Po6Gb4tTVwapVsGJF67JqVWiRPnCg/deYhW4aEyYcvowZA3l5Sf0VRHqTZA23F5Uo8rZytogkUvOBZpZfuZxdj+wid1Au0xdMp8/IPpEcO2HD7UXk88C0uJlzSgmzhH039nO8ImBL8kJLc1VVsHLloQX0ihWh9flIRoxov3geOzZ06RAROTrlbRFJa1m5WUx+YDJL/3Eplf+vkuVXLGfac9PIyk/sN+zpUFhfDcTftvky8O/Aq8AXzKzA3feZWS4wCng2BTGmVlNT6K6xZAm89lp4XLoUthzh/6rcXDj+eDjxxNblhBNg/Hjom9I5IUSkZ1DeFpG0l5WfxYkPnMjC6QupmVfD2pvXMv474xN6zpQX1u5eEf/czJqACnffaGaPAhcAjwIXAg+5++YUhJk8tbWhaG4poF97DV5/vf0h6/r0gUmTwhJfRI8dG4prEZEEUN4WkUyRV5bH5Acns+iMRbz9vbcZdNWghE6BnvLC+ig+BdxlZlOBIUDPurOlqSl03Xj1VXjllfC4fHmYVKWt0aNh6lQ4+eSwTJkS+j5nZyc9bBGRDvTsvC0iGad4ZjEjPjOCt7/7Nm/d8BbTXppGXFe2SKVdYe3uY+J+3gNcn7poIlZdDS+8AM8/HwrpBQtCC3W8nBw46aRQPLcU0lOnQmnbbosiIumhR+dtEekRxnx9DNt/tZ3qV6rZ88SehE19nnaFdY9SUxMK6aefhmeegYULDx8PeswYOO20sJx+OkybppsIRURERCKUU5jDyC+OZO0X17LxWxtVWGcE99C149FHw/Lyy4eODZ2TE4rnc8+FM86AmTPDGNAiIiIiklDDPjGM9V9dT+Uzlexbu4+CsQWRn0OFdXc1N4euHX/8Yyim169v3ZadHQrp886DWbPgzDOhsHsDlIuIiIhI1+UU51B+RTnbf72d7b/dzphbx0R/jsiP2FssWQK//S3cfz+8/Xbr+vJyePe74ZJLYPZsKE7eNJoiIiIicmQthfXuv+1WYZ1y9fXwhz/APffA/Pmt60ePhquugksvhVNP1fTeIiIiImmo/6z+WI5R/Wo1ByoPkNs/2uGJVVh3Rl0d/OAH8L3vwY4dYV1paSimr746dPFI0LAtIiIiIhKNnOIcCqcXUvNqDbULaym9INpR11RYd+TAAfjv/4ZvfKO1oJ4yBW64IRTUmsVQREREJKMUTSsKhfVrKqyTZ+FCuO66MPMhhBE87rwTLrxQrdMiIiIiGarf1H4A1C6tPcqeXafOwG25w3e+Ewrp114L40w/8kiY0GX2bBXVIiIiIhmsYFwYZq9+Y33kx1aLdbzGRvjYx+DXvw7PP/OZ0Erdr19q4xIRERGRSOSPzAdg/6b9kR9bhXWLpib4yEfgd78LhfSvfgWXX57qqEREREQkQvkjQmFd/3Y97o5F2BtBhXWL228PRXVhITzxRJgZUURERER6lJzCHLL6ZNG8v5nmvc1k98uO7NjqYw1h5sQ77gj9p//0JxXVIiIiIj1YdkkophurGyM9rgrr5ma48cZw0+KXvxxuUBQRERGRHiunOHTaaKxSYR2t//s/WLwYRoyAr3wl1dGIiIiISIJlF4cW66bqpkiPq8L6xz8OjzfeqAlfRERERHqB7MJYYV2nwjo6NTXhRsWsLLj22lRHIyIiIiJJkJUbSmA/4NEeN9KjZZqnnw5jV59+OpSVpToaEREREUkCyw1D7KmwjtL8+eHxne9MbRwiIiIikjQthXXzgeZIj9u7C+vXXw+PU6akNg4RERERSZqDLdYNarGOzptvhsdJk1Ibh4iIiIgkTVZeD+9jbWaXm9k6M9tlZt83s5zY+jIz+4mZfdvMvmtmuZGddOvW8Dh8eGSHFBHpDVKSs0VEItKju4KY2SjgUuAK4LPAdcBnYpsfBO5z95uAN4DbIznpvn2wZw/k5sLAgZEcUkSkN0hJzhYRiZBlh8KaaEfbIyfawx2z0cDH3L0RWGhmJwHnmdkLwER3fzW231+A1WZ2p7vXdeuMO3eGx/LyMNyeiIh0VvJztohIlGKln3sP7Ari7s/HEnSLzcBG4HxgXdx+m2M/zuj2SWtrw2NRUbcPJSLSm6QkZ4uIRMgs1mIdbU+Q9Cis23Eq8GNgOLCrzbYaYFi3z1AXazzp16/bhxIR6eUSn7NFRKIUq6uJtsE6bbqCHGRm44EKd19iZg7sa7NLPnCg7esqKiqYMaO1UWTOnDnMmTPnyCdqabFWYS0iKTB37lzmzp3b8jRjZ6hKWs4WEYlSrLDubFeQzubstCqszSwb+ATwpdiqzcDxcdsN6Adsafva8vJyFixY0PmTqcVaRFIovpA0s50pDueYJDVni4hEqYst1p3N2enWFeTzwN3uXh97/iQwMW77GKAB6H42ro+doqCg24cSEemlkpezRUQidLCPdcRdQdKmsDazm4HFQIGZjTWzjwGVwNbYHecAFwP3uHtDt0/YGLvvJju724cSEeltkp6zRUSi1JP7WJvZLcAdbVavcvf7zOwDwG1mth4oBL4SyUlbCuuctLgEIiIZIyU5W0QkSl3sY91ZaVFVuvudwJ1H2LYB+HjkJ1VhLSJyTFKSs0VEopSgFuu06QqSdE2xqXZUWIuIiIj0Kj2+j3XSqcVaREREpHdSi3XEVFiLiIiI9E4J6mOtwlqFtYiIiEjvohbriGm4PREREZFeybLUxzpaLU3/LZ3XRURERKR3aKmrm9UVJBothXVW770EIiIiIr2SuoJErLk5PKrFWkRERKRX0XB7UVNXEBEREZHeSS3WEVNhLSIiItI7Jaj8672FdWEhDB0KJSWpjkREREREkiinOIe8oXlk9Yu2FO69gzj/67+GRURERER6ldFfGc3or4yO/Li9t8VaRERERCRCKqxFRERERCKgwlpEREREJAIqrEVEREREItCrC+u5c+emOoQuU8yJl2nxQubFnGnxQmbG3BPp3yF6uqbR0vWMXiZdUxXWGUYxJ16mxQuZF3OmxQuZGXNPpH+H6OmaRkvXM3qZdE17dWEtIiIiIhIVc494LscUMbMKYEMXX1YG7ExAOImkmBMv0+KFzIs50+KFxMY82t3LE3TstHSMORsy872T7nRNo6XrGb10u6ZHzNk9prAWEREREUkldQUREREREYmACmsRERERkQiosBYRERERiUCvLKzNrMDM7jGz/zCz/zaz4lTH1JaZXW5m68xsl5l938xyYuvLzOwnZvZtM/uumeWmOta2zCzPzJaY2azY80y43qeb2efN7NLYNU7bmM1skpnda2afNbMfmdnJsfVpFbOZzTazV81sTNy6I8aY6viPEG+7n8PYtrT/LPYkqX5/9FRmdqeZeWxZkup4MlFXc50cXXvXNLY+7d+vvbKwBn4EvODuNwOPAD9JcTyHMLNRwKXAFcBngeuAz8Q2Pwjc5+43AW8At6ckyI59ARgT9zzdr/dHgUvc/W53f9jdd5LeMf8a+A93/y/gW8D9sfVpE7OZDQYKgZltNnUUY8riby/eo3wOITM+iz1J2ry/ewoz60cYbWF2bLkitRFlnmPMddKBI13TTHm/9rpRQcxsGLAOKHX3vbEWqCpgsruvT2lwMWZ2DvCyuzfGnt8FTAbuAP7k7sNi64cDq4HB7l6XqnjjmdmZwAnAbcC1hIIjba+3mZ0L3ArM9tiHId3fI2ZWB0x391VmVg4sAWaQZjGbWRbQBBzn7us7uq5AQ6rjbyfedj+H7v5uMzudNP8s9iTp/pnMVGbW8ofij919f0qDyWBdyXV6v3ZO22saW5cR79fe2GI9C9jp7nsBYv9prgPOTWVQ8dz9+Zb/zGM2AxuB8wmxtuy3OfbjjCSGd0SxvyavcPf74lbPIr2v93eBlcA9ZvaYmZ1F+sf8R+B/zKwI+BBwA2kYs7s3t1k1iyPH2NG2lMTbwecQ0vyz2APNIs3e35nOzAz4IPAdYLuZXZXikDJWF3OddELba5pJ79feWFgPB3a1WVcDDEtBLJ11KvBj0j/2mwhdE+KlbcxmNhGYBvzU3T8NPAU8BowkTWOO+RdgPzAfqHX3h0jj6xynoxgzIf6WzyFkRrw9ia53xDw4nfDV+j3Ab83sXSkOq6fQ+zVimfR+zTn6Lj2OA/varMsHDqQglqMys/FAhbsvMbO0jd3MLgIWuPuONpvSNmZCF4Td7r409vyHwFdJ75gBCgj9qocD3zOzdaR/zHD0GNM2/vjPYWxVJlzvnkTXO0HcfQ9wS6xF8DPA31IcUk+g92uCZML7tTcW1puB0jbrioAtKYilQ2aWDXwC+FJs1Wbg+LjtBvQjPWL/PDAthASEa/wIoatFul7vHOI+A+6+z8zeBLJJ35gBfgN80N0rY39s3U9IMOkcM3T82fMOtqVUO59DSO/PYk+UMXk7g/0AeCbVQfQQer8mXtq+X3tjV5CngWFmVgAQGyJrFPBsSqNq3+eBu929Pvb8SWBi3PYxhJu+FiQ5rvZcDZwct2whjKLwc9L3ei8FSmJ3ILdoJMSeljGbWRkw1d0rY6vuBKqBDaRpzHE6+uyl8+ey7ecQ0vuz2BOl8/ujp2gGFqU6iB5C79fES9v3a68rrN19G/AocEFs1YXAQ3E3H6UFM7sZWAwUmNlYM/sYUAlsNbOTYrtdDNzj7g2pirOFu1e4+9stC+Fu3gp330iaXm93X0X4GukKADPrD+QBvyVNYwZ2A/tjo1C02EX4IyGtYrbWry8MOv7spcPnsm28sXWHfQ7NbLy7LyBNP4s9UTq8P3oaMxtmZlebWVbsvf854CupjisTdSXXpSC8jNT2mmbS+7XXDbcHYGalwF3AemAI8GV3r01pUHHM7BbC0HrxVrn7JDMbTRjKbj1hnMevtBm5IC2Y2XrgWnd/Jp2vd6wF+PuElsZRhBsZV6R5zFMJNzAuBAYDz7n7s+kUc9yIJT8ijO/8Q3ev6CjGVMbfXryE7h/tfg5jr8mIz2JPkU7v757AzCYTGhaqgReA78UaG6QLjiXXSceOkI8HkSHv115ZWIuIiIiIRK3XdQUREREREUkEFdYiIiIiIhFQYS0iIiIiEgEV1iIiIiIiEVBhLSIiIiISARXWIiIiIiIRUGEtIiIiIhIBFdYiklbMrNzM/mxmG8zsS6mOR0REOqa83S2+SUoAAAoRSURBVEoTxIhIWjGzi4CngWJgFXC8u+9ObVQiInIkytutVFiLSNoys+eA8zVVuIhIZujteVtdQUQkLZnZSOCJ3pqcRUQyjfK2CmtJADP7opndm+o4WkQRj5ldaGYvm9msiMLqUczs7oiPZ8A/Ad+K8rgicijl6+4zszFm9lMz+0lEx8s2s+vNbEMn9n23mb16rL9rlLlbeTtQYZ1hzGySmf3JzP7HzJaYmZvZh1IdVxt/Be5PdRBxoohnEzA1gljSkpn9o5ntNLOtZnaGmY2Mvc/czH5oZv1j+800s2Vm9h0zK4itGwxsizikq4GfuXtj7PgiGUf5+phkYr7eEXvMj+h4WUAlMKoT+y4AJh1po5lNNrMHYu+9W8xsSNy2qHO38jYqrDPRw8DX3f06YBrwq5YNZjbFzM5NRVDx53b3Ze7+QiriaE9X4jnSNXT31cDOKOOK8t+ru8dy978CPwIa3f1ld98EXA80AK+6e2Vsv3mERP4Fd98Xe/l7gEe79QvEMbNbga8BfzezVcDEqI4tkmTK112Uifna3fcCW6I6l7sfABZ2ct/thCL8SNuXA3NjT+929/hCOrLcrbzdSoV1BjGzAcDxQB2AuzcDtwFuZiWEpG0piCtl545SJ36P5iSeKxXHegAYYWbTAdx9F/AUcEncufoDW2LvvRYnuPuqbp77IHe/w90nuPsMdz/B3Z+L6tgiyaJ8nVhpmK+jHgmiK/Efbd+m2GPbfs+R5W7l7VYqrDNIbOia1cBjZnZ6bN0GYB6hX9NxwBwz+yiAmd1uZrfGvgY62I/KzGaY2XfN7DNmttfMVpvZv5vZpWY238wuM7P/NbNqM/tE3OvaPV6bc19vZr+I72sW+zr0HjP7ppk9bWYzY+s7PF/c683MPmFmy81stpktNbPtZnZtJ84xIT6eo5zzsGvYjtLYV7s1ZnZ9m2v6tdi1+VMsEXf2mh12LjO7wszmmdm1ZvammX2rg+O192/fbjwdibVsLAeujFtdDrzbYt0+gPcCf46LswDYG/f8BDP7lpndbWbPWuhS8gEzW2Bm7zOz+82s0sxuNLPpsd9xl5n949HiE8kkyte9J1+350i/Y0fnMbPjzOw+C62//9mZ87Rz3r+a2W/N7Kyj7KfcnSjuriWDFmAysIHwF+p9wKC4beuBWbGfTwLWxn4uJ/w1PZDwF/dGYEZs273AS7Gf84DdwA8JY1FeBWzu6Hhtz03oY/YT4Bex9X0IxVr/2PMrgAqgqKPztfN7D4md8+bYvncC+4HRRzlH23g6PGf8NWwnhvXAd4C+sddtjK0vBn4ft9+LhKTYqWt2hHMNIHy99z1gJnBhR8dr82/fbjydfH/dAqyJ/TwO+CZQD1we936xuP3fC5wT+7kvMB/oG3v+XOw6FANVwO1AAfCB2PN/Ivxx/zli70EtWnrSgvJ1r8jXse1f6+R17Ojfej4wKe413sn3Wcu/5xnAV9vZPit2npy4dcrdCVrUYp1hPLQqvoOQaK4BlptZezdpvAl82MxygXNi6/5/e/cWakUVx3H8+/ekPXhIK5H0qGUgCV3wIfGhKDIryki6aIEIGYFi1gkiIhCi8iVQ6PYQimAQmmJhVmRUvlQWJUgRGQonKSMI7KJSecl/D/81uvZmz76c9u543L8PDOzLzFpr1pz5O3ut/4y9xIk8GSjyY/cQJwzufgw4DGxx90PEiXZRg/Kq23eUypsh5gJ/ecrRBd4kTsi7GtRXXW5R5oa07koi//emBnVUtKeVOku87ZFP9yUwIdvHC9OI0qNEn0KTfVayv78SAewtd//C3T9sobyy9jRjE3CpRTrI3cBa4H1ggZn1Aoc9Rd7kGmBnen0H8Y9XMQpyJ/Bs6uffgR0eedm7gPPcfbPH9PhuWjsGIsOC4nV3xOsaSvexTj1zgDHuXrRnV4t1zgQWuvvTTa6v2N0h5wx1A6R5ZtYDTHb3/cAjZraOuIN6DTArX9fdj5pZH/AE8GpRhLv/Yma7gBuJX9RTiYupU5tmr0+Q0oXKyitpal7GNGBk1q6TZrYfmFSvvkbc/W8zGwDGAuNaqGPQdVY5wenzZwrwvbs/X71SM31mZouIC9jCZR5Txp63tYVjUNqeRtx9n5ntJtJBxrj7gJltIkaR5gPvZu0u/jaK/L2LiRGmoqyDedHZ63+odBKlpclZRvG6Yv+6IV7nSvuxzrG5ivgRMVizgWlm1uvuR0rWOZn2QbG7g9Qhw8tIsvxXd/+KmGqbXr2imc0Glrn7So8nPOQeAmZaPPbpa+C5RhU3KK+eH4gR0PwxREbkHv5Xo4C9Ha6jGT9TmYeMmc1qoc+2ATOypebd5S2UV7M9ze8OrwOLgYGsfT1APzFtWphF5IsWfgKuNbPRWb118/xEzmKK15W6KV6X7mOdeg6lbUYxOKuI0f4XSr7/w0/fdK7Y3UG6sB5+Hjazidn7PmB7en2MuFljOvFopzFmdq6ZzUnfjzezcUSu3zpiGmcnaWoxGUHVr3QzswblVddtWRlbgSPAfamssam+YuSzrL4yE9I6fWnb7U3UkbenUZ35flSrNXphxCjUaOAdM7vF4iaZRsfgVF3ABHf/LluOZ1X0ZK/rlZe3u2Z7zGyUma1JaR71bCamoLcAuPth4D3gU698GsjNRJpIoejvjRbPwl4OFM8yrT4G1cd5WD+hQKSE4jVdE6+b7cea9RCxdBTx4wtidgIzm9hE7DZi9HgRMN/MFqRtHzOzFcRo/UfZ+ordndSJxG0tnVmIGyIc+A14DVhPTCWdn75/CjgAzCOmdvYRIwRziZGO9USQ2kbcSHEslfcncaLdChwnbpCZDDyZvn+gXnlVdfcDHxN5a5en764mRjpXE7mGV6TPS+sr2X9P2/cTqQlXZt+V1XFJ3p5GdeZ9WFX3POImvpdTXxTbLU7fXw98k47N6tTPzfbZvBr7ek+qbyPQlz6rd0wryippzyTiPzLY0MTf2otV7+8Fbqj6bFWN7a4DvgUOAivSZ7cRf2svARNJjxwDlhL/8K5N+3r7UJ9jWrS0a0Hxupvi9RTgE2LUfUaDfax3rOcSM4WfEzeS7wWWpX2vGbuJWZHjxEj1BcAO4qL+QWAJ8Tzv1UQaSrGNYncHF0udJ13CzMYDS939mfR+BDGKssTdVwxp4xowMwemeuQsyiCZ2ePuPqhHOYnI/0fxWnKK3cODUkG6z3IiLeCC7LPZZDelnYmyqSdNOw2SmfWY2ULgjaFui4g0RfFaFLuHGV1Yd59XiFysPWb2I/ABcMDdPxvaZpVLN3P0p7f3Wzz2TVrXC2x194GGa4rImUDxWkCxe1hRKoiIiIiISBtoxFpEREREpA10YS0iIiIi0ga6sBYRERERaQNdWIuIiIiItIEurEVERERE2kAX1iIiIiIibaALaxERERGRNtCFtYiIiIhIG+jCWkRERESkDf4F5IJ+IBK7x2wAAAAASUVORK5CYII=\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/mirka-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/mirka-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/mirka-earth.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 }