{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 16 - Apollo-AS-202 - 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 Apollo AS-202 suborbital flight test.\n", "https://en.wikipedia.org/wiki/AS-202" ] }, { "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('Apollo-AS-202', 5400.0, 400.0, 0.3, 12.0, 0.0, 3.0, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,8.29,0.0,-3.53,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,70.0)\n", "# bank angle = 70 deg arbitrary, actual profile will give more realistic results" ] }, { "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+17YcXAAAgAElEQVR4nOzdd5xU1f3/8ddnF5aqgIIUkd5BsSwiTZoKYu+oMT+jYf3aYvSrX2MjEZPYe4urMVGJCbFFRQU1gDSlC4JSRGQFKUvvsMue3x9nZhvLsmVm75T38/E4j5m5c+fezwzwmQ9nzj3HnHOIiIiIiEjlpAQdgIiIiIhIIlBhLSIiIiISASqsRUREREQiQIW1iIiIiEgEqLAWEREREYkAFdYiIiIiIhFQLegAIqVhw4auVatWpe6TnZ1No0aNqiagGKfPoih9HkXp8yhQFZ/FnDlzNjjnkuoDL0vOLkmi/93U+4tven/xryzvsbScnTCFdatWrZg9e3ap+6Snpx9yn2Shz6IofR5F6fMoUBWfhZmtjOoJYlBZcnZJEv3vpt5ffNP7i39leY+l5WwNBRERERERiQAV1iIiIiIiEZBUhXVGRkbQIcQMfRZF6fMoSp9HAX0WsSXR/zz0/uKb3l/8q+x7NOdchEIp4wnNTgf+CFzmnPsxtO1C4HHgcGA08L/OudzQcw2BPwFbgOrAnc65nOLHTU9Pd4k+7kdEEpOZzXHOpQcdR0mUs0VEiiotZ1dpj7WZNQbqAicX2tYCOB+4GLgV+DXw20Ivewt41Tl3J7AUGFVlAYuIJDHlbBGR8qnSwto5tw54v9jmlsA1zrk5zrnXgeeBgQBmdgrQ0Tk3I7Tvh8DNZlanqmIWETkU5xw5OTnk5eUFHUpExVrOzsuD3NxIHElEJDqqfIy1cy6v2OMp4Z8QQ1YDWaH7g4AVhfZdHbobkz+Zikhy6t27N2lpacycOTPoUCIuVnL2lCmQmgoDBlT2SCIi0ROLFy/2AP4Sun80sLHY89uBZsVflJ2dTXp6en7LzMyMcpgiIl61an5JgJycA4YSH1RmZmZ+vgIaRieyKlElOTs11d/u31/ZcEVEyq+sOTumFogxs3ZAtnNufmiTA3YX260GcMC3V6NGjRJ+0nIRiU1paWlA+QrrjIyM/KvPzWxDVAKLsqrM2aH/u6iwFpFAlDVnx0yPtZmlAtcBvyu0eTXQoNA+BtQBfq7a6EREDq569eoA7Nu3L+BIqk5V5+xwj7XGWItILIuZwhq4HXjMObe30LbPgI6FHrcC9gHqmhaRmBEurMvTY50AqjRnayiIiMSDKi+sQz0YAFZo293APKCWmbUxs2vMrJ1zbjawxsy6hXYdCjzrnEuebiERiXnhoSCJ2GMdKzlbQ0FEJB5U6RhrMzsM+EXo4dVm9hz+p8QHiu262Dn3auj+ZcBIM/sRP5/qPVURq4hIWSVqj3Us5WwNBRGReFClhbVzbjvwYqiF/THUDvaalcC1UQ5NRKTCErXHOpZytoaCiEg8iKUx1iIicSlRe6xjiYaCiEg8UGEtIlJJKqyjT0NBRCQeqLAWEamkRB0KEks0FERE4oEKaxGRSlKPdfSFh4Kox1pEYpkKaxGRSlKPdfSpx1pE4oEKaxGRSlKPdfTp4kURiQcqrEVEKkmFdfTp4kURiQcqrEVEKklDQaJPQ0FEJB6osBYRqST1WEefhoKISDxQYS0iUknhwlo91tGjoSAiEg9UWIuIVFJ4KIh6rKMnXFjn5YFzwcYiInIwKqxFRCpJQ0Giz0y91iIS+1RYi4hUki5erBo1a/rbPXuCjUNE5GBUWIuIVJJ6rKtGrVr+dvfuYOMQETkYFdYiIpWkixerhgprEYl1KqxFRCpJFy9WDRXWIhLrVFiLiFSSeqyrhgprEYl1KqxFRCpJPdZVI1xY6+JFEYlVKqxFRCpJFy9WjfCsIOqxFpFYpcJaRKSSNBSkamgoiIjEOhXWIiKVpKEgVUOFtYjEumpVfUIzOx34I3CZc+7H0LZawCPAdqABcKdzbtuhnhMRiQWJ3GMdSzm7Th1/u3NnJI4mIhJ5VdpjbWaNgbrAycWeegGY6py7G3gfeKmMz4mIBC5Re6xjLWfXr+9vt26N1BFFRCKrSgtr59w6fKLNZ2bNgCuAD0ObPgfONbNWpT1XJQGLiJRBol68GGs5O1xYb94ciaOJiERelQ8Fcc7lmVnhTQOADc65XaHnc81sBdAfyCnluR+rMm4RkYNJ5KEgsZSzw4X1li2VPZKISHRUeWFdgqOBjcW2bQeaAbmlPCciEhMSdSjIQQSWs1VYi0isi4XC2gHFr/Guge/54BDP5cvOziY9PT3/cUZGBhkZGREMU0SkZBXpsc7MzCQzMzP8sGHko4qawHK2CmsRCUpZc3YsFNar8VeOF3YY8DM+gR/suSIaNWrE7NmzoxKgiEhpKtJjXbiQNLMNUQksOgLL2SqsRSQoZc3ZsTCP9USgWWiKJsysOtAC+OIQz4mIxIREvXjxIALL2Q1CJfumTZE4mohI5FV5YW0FV8EYgHNuLTAWGBzafhrwjnNudWnPVWHIIiKlKlxYO+cCjiayYilnN27sb9eujcTRREQir0qHgpjZYcAvQg+vNrPnnHPZwPXAI2bWHWgCFB5oV9pzIiKBMzOqVatGbm4uubm5+YV2vIu1nH3UUWAG2dmQkwMJ8jGLSAKxROldSU9PdxpjLSJBqV27Nrt372bHjh3UCS8RWEZmNsc5l37oPRNHRXN206a+x/qnn6B58ygEJiJyCKXl7FgYYy0iEveSbJx1YJo29bdr1gQbh4hISVRYi4hEQJLNZR2YZqEZsX8+YJ4REZHgqbAWEYmARF59MZYcc4y/Xbky2DhEREqiwlpEJALUY1012rf3t8uWBRuHiEhJVFiLiESAeqyrRocO/nbp0mDjEBEpiQprEZEI0MWLVSPcY63CWkRikQprEZEICA8FUY91dLVuDampfoz1rl1BRyMiUpQKaxGRCKhWza+3lZubG3AkiS0tDbp2Befg66+DjkZEpCgV1iIiERBebKtgBXCJlvTQsgyzZgUbh4hIcSqsRUQiSIV19PXo4W+12K6IxBoV1iIiERDusZboO/lkf/vll8HGISJSnAprEZEI0FCQqnPccVCvHixfDj/+GHQ0IiIFVFiLiESQCuvoq1YNBg3y9z//PNhYREQKU2EtIhIBGgpStU47zd+OHx9sHCIihamwFhGJIPVYV41hw/ztxx/Dzp3BxiIiEqbCWkQkAtRjXbVatYJevfwiMR98EHQ0IiJeRAprMxsWieOIiMSreLp4MVFy9uWX+9vXXw82DhGRsHIV1mY21Mymmtn3ZvZDqP0IvBOd8ERE4sOePXsAqFmzZsCRFEj0nH355VCzJowbB0uWBB2NiAhUK+f+bwAPAvOAvNA2A66IZFAiIvFm9+7dANSuXTvgSIpI6JzdsCFcdRW8/DI8/TS88ELQEYlIsitvYf2Dc+6J4hvNbHGE4hERiUu7du0CoFatWgFHUkTC5+zf/tYX1n/9K9x5J7RsGXREIpLMyjvG+nozu7iE7edGIhgRkXgVLqxjrMc64XN2ly5wxRWwbx+MHBl0NCKS7MpbWPcBXjGz/YVaHvBiFGITEYkLeXl5MTnGmiTJ2Q88ANWrwxtvwPTpQUcjIsmsvENBRgEZwGxgf2ibAVdVNhAz6wzcBHwPtAcynXNfm1kt4BFgO9AAuNM5t62y5xMRiZTw+OqaNWuSkhJTs5gmRc5u0wZuvx0efBCuuQbmzYPYGpEjIsmivIX1J8D7zrm9hTea2csRiOUN4Dzn3GozawGMBzoDLwDjnHNjzGwo8BJweQTOJyISERs2bACgYcOGAUdygKTJ2SNHwn/+A999B3fcAc89F+0ziogcqLxdK/8C7jKzUwu1/sBdEYilM3BY6P5uoJ6ZNcNfvf5haPvnwLlm1ioC5xMRiYj169cDcNRRRwUcyQGSJmfXrAmvvQZpafD885rbWkSCUd4e6/vwP/ldDYSXGTOgCXBLJWN5Gz8W8EzgF8DNwABgg3NuF4BzLtfMVgD9gR8reT4RkYiI4cI6qXJ2jx7w7LNw3XWQkQEtWsCAAdE+q4hIgfL2WD8KHOOca+Wcax1qrYALIxDLDcAeYBawwzn3DnA0sLHYftuBZsVfnJ2dTXp6en7LzMyMQEgiIoe2bt06oHyFdWZmZn6+AqI1hiTpcvaIEXDDDbB3L5xzDsyaFZHDikiSK2vOLm+P9U/Oue0lbG9SzuOUpBbwT3xifirUy+HwPzEWVgPIKf7iRo0aMXv27AiEISJSPitXrgTg6KOPLvNrMjIyyMjIAMDMNkQlsCTM2WbwzDOweTP8859w2mnw/vvquRaRyilrzi5vj/UNxTeY2SXAU+U8TklGA+8450YBD+ET9hr8VeWFHQb8HIHziYhExNKlSwHo2LFjwJEcIClzdmqqH299ySWwbRsMGQL//ndVRiAiyaq8hfURZnYfgJnVNrNXgFfwPwVWmJk1BLo757aENv0R2AasBJqFpm/CzKoDLYAvKnM+EZFIWrJkCQAdOnQIOJIDJG3Orl7d91jfeKNfPOayy/yUfDkH9J2LiEROuQpr59xZwGdm9izwNdAB6A6UtLJXeWwC9phZ4d9RNwILgLHA4NC20/A9JKsreT4RkYjIzc2N2cI62XN2aqq/mPHJJ6FaNXj8cRg4EH74IYhoRCQZlDrG2sxOPchTK4FlwGNAS+B84NaKBuGcyzOz84GRZjYHaAzc4ZzbZmbXA4+YWXf8uMCMip5HRCTSFixYwK5du2jXrh1HHnlkoLEoZx/IDH77Wz9jyKWXwrRpcOyx8NBDvjc7ttbzEZF4d6iLFx/Bz1W6Bcgr4flX8VM3NaUSSRrAOTcfuK6E7ZuBEZU5tohItEwPraHdu3fvgCMBlLMPqk8fmD8fbr4Z/vUv+M1v4M034emn4eSTg45ORBLFoQrrUcBs59z60nYys8r+rCgiEpfGjx8PQL9+/QKOBFDOLlXDhn7c9aWX+in5vvoKevaEX/4S/vxnKMekLiIiJSr1RzDn3MeHStCh/d6OXEgiIvFh27ZtfPrpp5gZ55xzTtDhKGeX0QUXwJIlcOedfqXG11+Hdu38kJG1a4OOTkTimUaXiYhU0DvvvMO+ffvo27cvjRs3DjocKYfDD/fjrL/9Fi6+GPbs8cNCWreG226D1bpEXkQqQIW1iEgFOOd4+umnAbj22msDjkYqqm1beOstmDcPzj/fF9hPPgmtWsEvfgFz5gQdoYjEExXWIiIV8PHHHzN//nyOOuoohg8fHnQ4UknHHw/vvQdz5/qFZfLy4B//gPR06N/fX/C4Z0/QUYpIrKt0YW1mJ5jZaDP7XzNrVfmQRERiW05ODrfffjsAd955JzVq1Ag4orJTzi7dCSf4VRqXL/dDQg4/HCZPhssvh2bN/GwiCxYEHaWIxKpI9Fj/D3Cac+5xoHd4xS0RkUQ1cuRIFi9eTLt27bjpppuCDqe8lLPLoFUrv6DMTz/Bc8/5gnvzZr/gTPfucNJJ8NhjkJUVdKQiEksiUVjfBJwUuv8voHoEjikiEpM++OADHn74YVJSUnj11VdJS0sLOqTyUs4uh8MP9wvJzJ3r2403Qr16/v4dd0DLltC3ry++NaOIiJS7sDazM83sytD9k4H24eVqnXN5zrltEY5RRCQmTJo0icsuuwznHKNGjYqVuatLpZwdOSecUFBAv/uunw+7Vi2/muPNN/t5sAcM8LOL/Phj0NGKSBDKVVib2Z+At4ALAJxzM4FfmtkVUYhNRCRmvPnmmwwZMoQ9e/Zw3XXXcffddwcd0iEpZ0dHzZp+LuwxY2D9er+C47nnQrVq8MUXfj7s1q19IX7//X5MtnNBRy0iVaG8Pda9gebAlELbPgEeilhEIiIxZNu2bVxzzTVceeWV7Nu3jxtvvJHnn38eMws6tLJQzo6yunX9hY3vv++L7PDKjnXrwtdfwx/+4Mdkt23rL4acPBn27w86ahGJlkMtaV7cZ865LWZW+P/eSbk0roiULicnh9mzZzN58mQWLlzI8uXLyc7OZvfu3eTm5lKvXj3q169P8+bNOe644zjuuOM49thjadOmDSkpwc8Eum/fPl555RXuv/9+1q9fT40aNXj88ce58cYbgw6tPJSzo2D3br+ozOrVfujHscdCt27QpQsMH+7b3r0wYQL85z++6F6xws+P/eSTfmn1c8/182afdpofTiIiicFcOX6fMrPhQBvgaHwPyCX4nxhvds49H5UIyyg9Pd3Nnj07yBBEkp5zjokTJ/KPf/yDd955h61bt5b7GHXq1KFbt275xXa44G7QoEEUIj7QypUr+dvf/kZmZiZr1qwBoFevXvz1r3+lc+fOUTmnmc1xzqVH4bjK2VGwbBl06FDyc61a+SI7XGx36wbt2/sFaN57zxfa339fsH+dOjB0qC+yzzoLquivuYhUQmk5u1yFdehgpwBXAy2ANcA/nXOfVzbIyornJC0S73JycvjnP//JY489xjfffJO/vVOnTvTv35/09HTatWtHs2bNqFWrFqmpqWzbto3Nmzfzww8/8M0337BgwQIWLFjA6oOsJd28eXOOPfZYunbtmt+6dOlCnTp1KhX7li1bmDFjBtOmTWPs2LHMmzcv/7muXbsyatQoLrjggqgO/YhWYR06tnJ2FDz1FNx6a9n2TU31xXWXLr6BL86//RYK/XOhWjW/GM0FF8B550Hz5pGPW0QqL6KFdQkHTwUGOOf+W6kDVVK8J2mReOSc4+233+auu+5i+fLlADRt2pRrr72WK664okI9vBs3bixSaC9YsICFCxeye/fuEvdv1KgRLVq0oGXLljRv3pwGDRrkDzOpUaMG4RyXk5PD5s2b2bRpE+vXr2f58uUsW7aMrGITEdetW5ezzjqL6667jgEDBlTJWOpoFtYlnEs5O0I2b/azhDz9NGzcWLC9TRvo1Qt27YKFC30PdUlftWalX9Q4YoQ/fvzN6CiS2CpcWJvZB0DdQxy/IbDBOTeo4iFWXiIkaZF4MnPmTH7zm98wY8YMADp06MBdd93F5ZdfHvGVCPfv38/y5ctZtGgRCxcuZNGiRSxatIglS5aQk5NTqWOnpaVx0kkn0atXLwYNGsTgwYOpWbNmhCIvm0gV1srZwdi5E157DV580RfSYf37w5VX+iEe69fDokW+lzp8+8MPh76QcdGigl5uEYkNpeXsQ128uCC0zxIgDxgBzAMKZ8PWwBERiFNE4sCWLVu45557ePHFF3HO0bhxY+6//36uvfZaqlUr7/XQZZOamkqHDh3o0KEDF1xwQf72/fv3s27dOrKyssjKymL16tVs2bKFrVu3smXLFvbt2weAmZGamsoRRxzBEUccwZFHHkmbNm1o3749LVu2pHr1hFkjRTk7AHXqwA03wPXXw5dfwksv+WXRv/jCt7Q0GDbMF9n33ltwsWJOjr8AcsUKP+91Vpbv+d64EXbsgD59IErD+kUkSg7VY10vtM+W0ONRzrmRJez3vnPuvOiFeWiJ1PshEoucc4wZM4Zbb72VtWvXUq1aNW677Tbuu+8+6tY9VCeplCaCPdbK2TFi61a/iMybb/rZQfLy/PbatYterHiE/osjEndKy9mlzmnlnNsaTtAhjUPj8woffABwcqWjFJGYtWLFCoYNG8bll1/O2rVr6d27N3PnzuXhhx9WUR1DlLNjR7168KtfwWefwapV8MQT0KOHH3f97rvwy1/CUUf56faee84PCxGR+FfeyWLfBuab2V/M7FEzGwt8BjwT+dBEJGi5ubk8/vjjdOvWjXHjxlG/fn0yMzOZMmUKxx57bNDhyaEpZ8eApk39DCIzZ8JPP/lCevBg/9x//+uXQ2/b1k/h95vfwEcf+XHbIhJ/KjLdXjPgBuA4YBvwlnPu/SjEVi6J/rOiSFWbM2cOI0aMyJ9+bvjw4Tz55JM0adIk4MgST5Sn21POjlGbNvki+oMPfM924Wnf09Lg1FNhyBAYNMiv3piaevBjiUjVifZ0eynACc65OZU6UCUpSYtExo4dOxg5ciRPP/00eXl5tGzZkhdeeIFhw4YFHVrCquLp9pSzY1Buru/RHjfOt9mzi07FV7++L7QHDvSrPR53HMTAAqUiSakys4IUP9CrJWxuBOwHzq9AbCWd4xSgL/A9MBXYCTwCbAcaAHc657ZF4lwiUtTHH3/M9ddfT1ZWFikpKdx2223cf//9Gkcdp5Sz40e1atC7t2+jRsGGDb4X+9NPYdIkP2vIBx/4Bn6Fxv79fZF96ql+pccoTcojIuVQ3n+G7YH/AoW7uZvjl8qtNDP7FdDWOXdvoW1/A8Y558aY2VDgJeDySJxPRLy1a9fy29/+ljFjxgBwwgkn8PLLL3PSSScFHJlUknJ2nGrYEC6/3DfwhfWkSb5NnOin5vvPf3wDP+Vfz54Fxfkpp2h5dJEglGsoiJl1dM4tKbatE3C1c+53lQrErD9wH3C6CwUVGhu4AmjgnNtlZtWArUBX59yPhV+vnxVFyi8vL49XX32VO+64gy1btlC7dm1GjRrFLbfcErU5qeVA0RoKopydmJwrKLQnToTp0yG08GkRXboUFNq9e/tl1TV8RKTyoj3GujUwyznXsJLHmQNMBwxoBzwAtAQedc4dXWi/haFtrxV+vZK0SPksXryY6667jsmTJwMwdOhQXnjhBVq3bh1wZMmnisdYK2cnoHXr/OI006fDtGl+jHZofaR89erBSSf5af969ID0dGjRwi+tLiJlF8kx1iso+pNiKtAYeKPi4fleFeAE4FfOuQVmdgcwDvgjsLHY7tuBZsWPkZ2dTXp6wXvMyMggIyOjMmGJJKQdO3bw0EMP8eijj7Jv3z4aNWrE008/zfDhwzF9w1aZzMxMMjMzww8rVeQejHJ28mjc2C86c35o5PzevTB3ri+0p0/3RfeaNX6xmgkTCl7XqFFBkR0uuBs3DuY9iMSysubs8g4FyQTepCBR5wHrnHNLKxhn+LgXApnhHhQzqwVkA6OAi5xzPQvtOxd40zn3WOFjqPdDpHT79+/n73//O/feey9r164F4JprruHRRx/lCC3/FqgoDgVRzpZ8P/8Ms2b53uxZs3zbtOnA/Zo3hxNO8O34431r1Uo92yJhEeuxBkY75yYXO3hHM2vnnPu+whH6OPJjcc7tNrNl+N6V4pdfHAb8XIlziSSdCRMmcNtttzF//nwATj75ZJ588kl69+4dcGQSZcrZkq9ZMzjvPN+gYKx2uMiePRvmzPErRa5aBR9+WPDaevUKiuzjj/dFd+fOfr5tESlQ3sJ6KDC52LZ1wET8z4IVtQCoZ2aNnXPrQtty8cm4mZnVCiXu6kAL4ItKnEskacyfP597772XsWPHAnDMMcfw8MMPc9lll5Giq5iSgXK2HJQZtG7t26WX+m15ebB0KcyfD/Pmwddf+9v16+GLL3wLq14dunYt6Nk+7jjo1s3PaCKSrMpUWJvZTfixc4eZ2Z0l7PJZZYJwzi02s0+Ai4Hnzaw+kAb8AzgTGAyMBU4D3nHOra7M+UQS3eLFi/nDH/6QP31e3bp1ueuuu7j11lupVatWwNFJtClnS0WlpECnTr5ddlnB9jVrfJEdLrS//hqWLSvYVliTJr7A7tbNz6/drZufoUTT4UsyKFNh7Zx7zswmA/8PeKbY0/ucc2siEMsvgafNLA3fw3G5cy7XzK4HHjGz7kATQFe3iBzEihUrGDVqFK+//jp5eXmkpaVx/fXXc9ddd9FYVyQlDeVsibSmTX0788yCbdu3w4IFBcX2woW+rV3r2+efFz1GmzYFBXe46O7QQcNJJLGU9+LF6s65nBK2t3DOZUU0snLShTCSzJYuXcrDDz/MG2+8QU5ODtWqVePaa6/lnnvu4Zhjjgk6PDmEKF68qJwtVSovD1auLCiyv/nG3y5eDDkH/E30q0V27OjHa3fu7HvKO3f2BXedOlUfv0hZVPjiRTO7HfjYOfdtaNP/mVlq8d2AgcCAygYqIuUzb948HnzwQd5++22cc6SkpPDLX/6S3//+97Rp0ybo8KSKKWdL0FJSCsZtn3NOwfacHD90JFxoh9vy5bBokW/FtWxZUGiHh6d07uynCNQMJRKrDjUUpDVQv9DjE4E6wFoKpm8y/M99IlIFnHNMnjyZBx98kPHjxwOQlpbG1VdfzR133EG7du0CjlACpJwtMal6dT/OukuXomO3d+6E777zPdqLFxfcX7bM93yvXAmhNJevQYMDe7fbt/dDTWrUqNr3JVJcqYW1c+7GYpv+BCx0zhVZz8nMukQ6MBEpas+ePfzrX//imWeeYd68eQDUqVOH//mf/+G2226jWbMD1uCQJKOcLfGmTh2/OE16sR/Vc3JgxYqCQrvw7ebNBQvfFJaS4leSbN/et3DB3b69n4e7evUqe1uSxMo73d6xzrm5hTeEVuBqC3xb8ktEpDLWrFnDiy++yF/+8heys7MBaNSoETfeeCM33XQTRx55ZMARSgxTzpa4VL26L4w7dCiYdxv83Ntr1xYU2d9953u3ly3zc3KH22fF5r2pVs0X18UL7vbtfTGeWnzAlEgFlXW6vX74nw9PDy2RW9gR+KvOPzzghSJSIc45pk6dyksvvcS///1vckJX/Rx//PHccsstDB8+nJo1awYcpcQq5WxJVGYFM5QMHFj0uX37fC/3smV+Lu5wwb1sGfz0E3z/vW+ffFL0dWlpfhhJ4WK7bVto1w6OOcYX5SJlVda/LluA1/Hj9/oUe24f8EokgxJJVhs2bOD111/n5ZdfZvHixQCkpKRw4YUXcsstt9CvXz9MV+3IoSlnS9JJS/MzjHTseOBzu3fDDz8cWHAvW+aXeg+P8S6uenXf092unW/hgrttW3+BpsZ0S3Flncf6GzMbAFzinFNCFomgvLw8Jk6cyMsvv8x7773Hvn1+OGyTJk245pprGDFiBK1atQo2SIkrytkiRdWq5VeJ7Nr1wOd27PA92eFCe/ly/3j5cli9umB7cWZ+GEm42C5eeGu6wORU5h84nHNbOUgvh5n1dc5NjVhUIklg6dKljB49mtGjR7NihX7VV5YAACAASURBVP+13swYNmwYI0aM4KyzzqK6rraRClLOFimbunX9kuzHH3/gc7t2+Z7uwsV2eEhJeNaSlSthwoQDX9ukyYEFd/h+gwbRf18SjEPNY/0BcKhFSA3oChwVqaBEElV2djZjxozhjTfeYObMmfnbmzdvzrXXXss111xDixYtAoxQ4plytkhk1a5dsFJkcfv2+aK6eMG9fLkvxsMrUE4t4b+wRxxRcsHdrh0cdZTm6Y5nh+qxXoIfq7eagjlQizPg7EgGJZJIduzYwUcffcTo0aMZN24cubm5ANStW5eLL76Yq666iv79+5Oqy9Kl8pSzRapIWlrBxY7F7d/vh5EULrYL39+0CWbO9K242rX9xZTh1rp1wf1WrfzzErsOVVg/DGx3zu0tbScz09XlIoVs376djz76iLfeeotPPvmE3bt3A5CamsqwYcO46qqrOPfcc6mtDCmRpZwtEgNSU/346xYtYNCgos85B+vWlVxwf/+9n6c7vDJlSZo2PbDgDremTf183hKcQy0Qs6GMx7kAXWUuSW7btm2MHTuWt956i3HjxrFnz57853r16sVll13G8OHDady4cYBRSiJTzhaJfWZ+/HWTJtC374HPb9nipw384YeCFn7844+wZo1v06Yd+NoaNXyv9sF6vA87LNrvTio1O6OZtQBuBG5CSVqS0Jo1axg7diwffPABn332GXv3FnQU9unTh0suuYSLLrqI5s2bBxiliKecLRL76teHE07wrbjwEJPiBXe4rV8PS5b4VpKGDUsuuNu0gebNNWd3JJT7IzSzVOBcYARwOrAOP6ZPJOE551iwYAEffvghH3zwAbNmzcp/zszo168fl1xyCRdeeCFHH310gJGKeMrZIomj8BCTAQMOfH7HDl9sFy+4w0X4hg2+lTS2u1o1aNmy5CEmbdpoJpOyKnNhbWatgV8DvwKOBFYC5wDjgV5RiU4kBuzdu5cvvvgiv5jOysrKf65mzZqcfvrpnHPOOZx99tk0bdo0wEhFCihniySfunXh2GN9Ky4vz4/tLl5wh4vu1av9OO/ly0s+dr16RQvt8JLzHTpA48aaySTskIW1mV2C7+kYBHwPPAG8BtztnBsX2m161CIUqWLOOZYtW8b48eMZN24ckyZNYteuXfnPN27cmHPOOYdzzz2XwYMH6wJEiSnK2SJSkpSUguXg+xRfjxXYs8eP4S5pbPcPP8DWrTBvnm/FHXZY0UI73Nq39wV5MilLj/Vu/BK404FfO+eWApjZwaZyEok727dvZ8KECYwbN47x48fnL9gSdtxxx3HOOedwzjnn0KNHD1J02bXELuVsESm3mjWhUyffinPODyEJF9vh1SiXLvXjubdsgTlzfCuuceOSi+62bRNzSfhDFtbOubHA2NBFLyPMrBXwPpC/JJyZmXNOSVvixv79+5k7dy6ff/4548ePZ9q0afnzSwMceeSRnH766QwZMoQzzjiDZs2aBRitSNkpZ4tIpJlBo0a+9exZ9DnnYONGX2QXb8uW+eEn69bBlClFX5eS4mcw6dKlaOvc2Q9piVflWdI8C7jPzKoB5wPNzewB4K/AGUBmdEIUqby8vDwWLVrEhAkTmDBhAl988QVbt27Nfz4lJYXevXszdOhQhgwZwkknnaQFWySuKWeLSFUw87ONNGwIvXsXfS4vD1atKrnoLjzMZOzYoq9r0eLAgrtLl/gYVmKV6bQws3bAzcAI51ygA03T09Pd7NmzgwxBYohzju+//z6/kJ44cSLZ2dlF9mnXrh0DBw7kjDPOYPDgwTTQJc8SEDOb45xLr4LzKGeLSEzYt88viPPtt0XbkiX+uZIcfTR07eqnIjzxRH/btm3VL4pTWs6u1IyFzrnvgVvMrIRRNSJVxznH8uXLmTJlCpMmTWLChAmsWrWqyD7NmjVj8ODBDBo0iIEDB9KyZcuAohUJhnK2iMSKtLSCnujCcnN9L3bxgvu77/zMJatXw6efFux/2GEF83736AH9+vke76BEZCpw59zrkTiOmaUBs4BbnHOTzKwW8AiwHWgA3Omc2xaJc0l8Cw/tmDx5MpMnT2bKlCmsWbOmyD4NGzZk4MCBDBo0iEGDBtG+fXtM8wGJRCxng/K2iERWtWoFFzief37B9v37YeVKmD+/YHaSuXPh559h8mTfwlq2hFNP9W3YMKjKy6RibY2dO4BWhR6/AIxzzo0xs6HAS8DlQQQmwcrJyWHu3LlMmTKFyZMnM3XqVDZv3lxkn4YNG9KvXz9OPfVUBg0aRLdu3TR7h0j0KW+LSNSlphbMoX3BBQXb160rKLK//NJfJLlyJbzxhm8Ap5ziX3PVVX66wWiKmcLazHoDa4DNocfNgCvwy+8CfA68Y2atnHM/BhKkVJndu3czY8aM/N7o6dOnF5lLGqB58+b0798/v5ju1KmTeqRFqpDytogErXFjGDrUN/A92wsX+h7szz7z7auvfLvnHrjoIrjrLujePTrxxERhbWZ1gIudc7eZ2cjQ5gHABufcLgDnXK6ZrQD6Az8GEqhEzZo1a5g2bRrTp09n2rRpzJ07t8j0dwAdOnTg1FNP5dRTT6Vfv360bNlShbRIQJS3RSQWpab6orl7d7j5Zr/M+6ef+t7rDz6AMWPg3/+G666DJ5/083dHUkwU1sCdwEPFth0NbCy2bTugCYXj3P79+/n222+ZNm1afiu+IEtKSgrdu3fPL6T79u1LkyZNAopYREqgvC0iMa9uXbjwQt9++gkefRRefBH+8hdYvNj3aFeLYDUceGFtZkOA2c659cWecvgVxAqrAeSUdJzs7GzS0wtmPsnIyCAjIyOSoUoF7dy5kxkzZuQX0V9++SXbthW9lqlu3bqccsop9OnThz59+tCzZ08OP/zwgCIWib7MzEwyM/Onkm4YZCzlFYm8rZwtIlXtmGPgmWfgmmv80JFJk+Ddd+HSSw/92rLm7MALa+B24IRCP+k3wK8S9kTofmGHAT+XdJBGjRqhOVFjw6pVq4r0Rs+fP5/9+/cX2adFixb5RXSfPn3o1q0b1SL5X0aRGFe4kDSzDQGHU16VztvK2SISlPr1/TR969bB9u1le01Zc3YsVDJX4Hs0wr4EbgNmAHeYWS3n3G4zqw60AL4IIEY5iNzcXL755pv8Inr69OlkZWUV2Sc1NZWTTjopv4ju3bs3zZs3DyhiEYkA5W0RiSvOwbRp8Pe/w+uvQ04OtG8Pl1wS2fMEXlg754osh2dm+4Fs51yWmY0FBgNjgdOAd5xzqwMIU0J27NjBjBkzmDJlCtOmTeOrr75ix44dRfapV68evXr1yi+ke/ToQd26dQOKWEQiTXlbROLBli3w3//C+PEwbpwfYw1+GfYrroCnn4ZIjzoNvLA+hOuBR8ysO9AE0AC8KrZ+/XqmTp3K1KlTmTJlCvPmzTtgWEebNm2KDOvo0qWL5o8WSV7K2yJS5fLyYOlSP5f19On+9ttvfU912NFHw5VX+jHWHTtGJ46YK6ydc60K3d8MjAgumuQSXhY8XERPnTqVpUuXFtknNTWVHj160Ldv3/xCWrN1iCQ35W0RqUo7dsCCBX4VxnBbsACKLXdBWhr07AlnnukvVuzeHaLd7xdzhbVUndzcXBYsWJBfRE+dOpW1a9cW2ad27dr06tWLvn370q9fP3r27KlhHSIiIhJVzsH69fDdd35avMJt5cqSX9O8uS+ke/eGXr3gxBOhRo2S940WFdZJJCcnh1mzZjFhwgQmT57Ml19+ecD46EaNGuUX0X379uX444+nevXqAUUsIiIiico52LABli8v2pYu9QX0li0lv656dejSBY4/vmAxmO7d4cgjqzb+kqiwTmD79+9n3rx5TJw4kQkTJjBlyhR27txZZJ+2bdvmF9H9+vWjffv2Ws1QREREIiI3F1atOrB4DrfSprurVw86d4ZOnYq2Nm18cR2LVFgnkLy8PBYtWsSECROYOHEikyZNYuvWrUX26dy5MwMHDmTAgAH07duXpk2bBhStiIiIxLu9e/1sGz/+6IdohFv48apVUGzOgyIOPxzati3a2rf3BfVRR/kZPOKJCus4t379esaPH8+4ceP47LPPyM4uMgsWbdq0YeDAgQwaNIiBAweqkBYREZEy27nzwGK58OO1a4vOvFGSpk0PLJ7D7cgj4694Lo0K6ziTm5vLzJkz+eSTTxg3btwBK5c1a9aMQYMG5RfSrVq1CiZQERERiWnOQXY2ZGX5Qrn47cqVsHFj6cdISfFLhbdqBS1bFrTw42OOgZo1q+LdxAYV1nFg06ZNjB07lo8++ojPPvuMzZs35z9Xo0YNBgwYwNChQxk6dCgdO3bUGGkRERFh3z4/FONghXNWFuzZU/ox0tKgRYsDC+bw/aOPhmqqJvPpo4hRq1at4v333+e9995j0qRJRRZlad++PUOHDuXMM8+kf//+1K5dO8BIRUREJAhbtx68YF65EtasOfQwjfr1Cwrnkm6bNIn+3M+JRIV1DFm+fDlvvfUW7777LrNmzcrfnpqayuDBgznvvPMYNmwYbdu2DTBKERERiba8PF8YFx6WUbxw3rat9GOkpPge5dIK50gv6Z3sVFgHLDs7mzFjxvCPf/yDr776Kn97rVq1GDJkCBdccAFnn302RxxxRIBRioiISCTt3u0L5IP1Nq9aBTk5pR+jVq2DF8wtW/qiOlanpUtUKqwDsHv3bt577z1Gjx7Np59+mj/Mo06dOpx33nlcfPHFDBkyREM8RERE4pBz/qK/0oZpFJvEq0SNGpVeOCfajBqJQIV1FVq8eDEvvfQSr732Wv4FiKmpqQwbNowrr7yS8847jzp16gQcpYiIiJQmJwdWrz544ZyVBbt2lX6MatX8jBkHK5xbtPA90hJfVFhHWU5ODm+//TZ/+ctfmDx5cv729PR0rr76ai699FIaNWoUYIQiIiJS2LZtpc+k8fPPfgx0aerVO/RFgampVfN+pOqosI6S7du388orr/Dkk0/y008/AX6ox5VXXsl1113HiSeeGHCEIiIiyWn7dr+c9vffw4oVBxbOW7aU/nqzQ18UWK9e1bwXiS0qrCNsw4YNPPHEE7z44otsCf3L7NSpE7/5zW+48sorOVyX34qIiERVeIxzuHgO34bvr19f+utr1jz0RYFpaVXzXiS+qLCOkG3btvHEE0/w+OOPs2PHDgD69u3L//3f/3HWWWeRokkgRUREIiY8HV3xojl8u3XrwV9bo0bBktpt2hQseBIunBs21EWBUjEqrCspJyeHF154gQceeICNoXU/zzzzTO677z569eoVcHQiIiLxbedOWLoUFi8u2pYt81PWHcxhh0G7dr54Lnzbrh00a6ZFTyQ6VFhXwtSpU7n++utZuHAhAH369OHBBx+kX79+AUcmIiISX9avh2+/LSicv/vO32ZlHfw1DRuWXDy3beunqlOvs1Q1FdYVsHXrVm699Vb+9re/AdCmTRueeuopzj77bEz/ikVERA4qJweWLIH58wvaggWwdm3J+1erBu3bQ6dORVuHDn45bpFYosK6nCZPnsxVV11FVlYWNWrU4He/+x133nkntTTZpIiISBH79vmieeZMmDULvv7a90rv23fgvnXrQrdu0Llz0QK6dWutHijxQ4V1GeXl5XH//ffzwAMP4JwjPT2d0aNH07Fjx6BDExERCZxz/sLBmTNhxgx/O29eyUV0mzbQvbtvxx3nb1u10rhniX8qrMtg+/btXHXVVbz//vukpKRw99138/vf/57q+i+0iIgkqZwcXzhPnuzbtGmwadOB+3XqBCefDD16wIknwrHH+gsLRRJRzBTWZnYh8DhwODAa+F/nXK6ZNQT+BGwBqgN3Oudyqiqun3/+mSFDhrBw4ULq16/PmDFjOOOMM6rq9CIiMSlWc7ZEz549vid6yhRfSE+f7mfsKKxJE+jZ0xfSPXtCeroWSpHkEhOFtZm1AM4HLga6Ai8CPwGPAW8Bv3POzTCz/wFGAXdVRVwrV65k8ODBLF++nI4dO/Lhhx/Svn37qji1iEjMitWcLZH3/ffwyScwbhxMnHjg9HYdOsCpp/rWr5+fA1rX8Esyi4nCGmgJXOOcywXmmFk3YKCZTQU6OudmhPb7EFhiZn90zu082MEiYdWqVZx66qlkZWVx4okn8umnn3LkkUdG85QiIvEi5nK2REZenu+Jfust+Ogjv9BKYcceC/37FxTSTZoEE6dIrIqJwto5N6XYptXAYcAgYEWh/VaHprNLB76IVjxbt25l2LBhZGVlccopp/DJJ59QX3P6iIgAsZezpXKc8xca/utfvqBevbrguQYN4IwzYOhQGDIEmjYNLk6ReBAThXUJegCPAhnAxmLPbQeaFX9BdnY26enp+Y8zMjLIyMgo94n379/PpZdeyjfffEOnTp346KOPVFSLSMRlZmaSmZkZftgwyFgiILCcLRW3eTOMHg2ZmRBa5wzwwzkuvRQuuMCPlU5NDS5GkVhR1pwdc4W1mbUDsp1z883MAcUXLK0BHHAhTKNGjZg9e3alz//QQw/x6aef0qhRIz7++GOOOOKISh9TRKS4woWkmW0IOJwKCzpnS/nNmQPPPgtjxvgLEgGOOgquvBIuu8wX0xonLVJUWXN2TBXWZpYKXAf8LrRpNdCh0PMG1AF+jsb5v/zyS0aOHAnAG2+8QevWraNxGhGRhBB0zpayy8uDsWPhiSfgi0KDck47Da67Ds49F9LSgotPJFHEVGEN3A485pzbG3r8GfA/hZ5vBewDIt7NkZOTQ0ZGBnl5edxxxx0MGTIk0qcQEUk0geVsKZvdu+Fvf4Mnn/QzfICfQ/rXv4Ybb4S2bYONTyTRxExhbWZ345NvLTNrAwwAJgNrzKybc24hMBR41jlXwjpOlfPcc8+xcOFCWrduzf333x/pw4uIJJSgc7aUbu9eeOUV+NOfYM0av61lS7jlFrj2Wjj88GDjE0lUMVFYm9m9wAPFNi92zr1qZpcBI83sR6AucE+kz79z507+/Oc/A/DMM89Qq1atSJ9CRCRhBJ2z5eCcgzfegPvug6wsv+2EE+Cuu/zFiNVi4ltfJHHFxD8x59wfgT8e5LmVwLXRPP9LL73Ehg0b6NmzJ2eddVY0TyUiEveCztlSskWL4Prr/cqIAF27wqhRvqDWxYgiVSMmCusg5eXl8eyzzwJwzz33YMo+IiISR5yD55+H226DnBw/w8cjj8AvfqGp8kSqWtIX1pMmTeLHH3+kRYsW6q0WEZG4kpvrL0R87TX/OCMDHnrIL+wiIlUv6Qvr10LZ6Fe/+hUpKSkBRyMiIlI2eXnw//4fvPkm1KkDr77qF3YRkeAkdWG9f/9+xo4dC8Dw4cMDjkZERKTs/vxnX1QfdhiMGwe9ewcdkYgkdRftV199xaZNm2jXrh0dO3YMOhwREZEymTABRo70FyWOGaOiWiRWJHVh/fHHHwNw9tln66JFERGJC5s3+yEgzvlp9c48M+iIRCQsqQvryZMnA3DaaacFHImIiEjZ3HgjrFoFPXv6wlpEYkfSFtZ79+5l1qxZAPTWb2giIhIH/vlP32rX9gvBaMEXkdiStIX1nDlz2Lt3L127dqWB5iUSEZEY99NPcMMN/v6TT0L79sHGIyIHStrC+ssvvwTUWy0iIrHPOfjVr2DLFjj7bBgxIuiIRKQkSfsjUqtWrTjrrLMYNGhQ0KGIiIiUyswvV756NbzyipYoF4lVSVtYX3TRRVx00UVBhyEiIlImF10E55+vZcpFYlnSDgURERGJNyqqRWKbCmsRERERkQhQYS0iIiIiEgEqrEVEREREIiCpCuvMzMygQ4gZ+iyK0udRlD6PAvosYkui/3no/cU3vb/4V9n3qMI6SemzKEqfR1H6PAros4gtif7nofcX3/T+4p8KaxERERGRGGDOuaBjiAgzywZWHmK3hsCGKggnHuizKEqfR1H6PApUxWfR0jnXKMrniCllzNklSfS/m3p/8U3vL/6V5T0eNGcnTGEtIiIiIhIkDQUREREREYkAFdYiIiIiIhGQFIW1mdUys2fN7M9m9qKZHR50TEExswvNbIWZbTSzp82sWtAxxQIzSzOz+WY2IOhYYoGZnWJmt5vZ+WbWMOh4gmJmnc3seTO71cxeMLPjg44pWSVDHjez081shpm1CjqWSEv07x4z629mS8xss5k9FXQ80ZLo35Vm9kczc6E2vyLHSKi/2KV4ARjnnBtjZkOBl4DLA46pyplZC+B84GKgK/Ai8BPwWJBxxYg7gFZBBxELzOxXQFvn3L1BxxID3gDOc86tDv37GQ90DjimZJXQedzMGgN1gZODjiXSEv27x8zqAX2AU4DewH/MbKxz7vNgI4uKhP2uNLM6+AsXTw9tqsjF1YlfWJtZM+AK4MbQps+Bd8yslXPux8ACC0ZL4BrnXC4wx8y6AQNJkORWUWbWG1gDbA46lqCZWX/gSgoSS7LrDBwWur8bqBdgLEkrGfK4c26dmb0fdBxRkujfPfuAB52fDeIjM1sA7A84pohLgu/KEcBiYKpzbk9FD5IMQ0EGABucc7sAQv+wVwD9gwwqCM65KaH3H7YayAoqnlgQ+h/qxc65V4OOJUY8AXwHPGtm48ysT9ABBext4BUzOwz4BXBzwPEkqwEkQR53zuUFHUM0JPp3j3Nud6ioDn+nfAdMCjSoCEv070ozM2A48Diwzswq/GtYMhTWRwMbi23bDjQLIJZY0wP4S9BBBOxO4KGgg4gFZtYROAF42Tl3E/BfYFzoJ+pkdQOwB5gF7HDOvRNwPMlKeTyxJOR3j5n1Az4BqgO1Aw4n0hL6u9J5p+CHgjwL/MPMzqzIsZKhsHb4n3ALqwHkBBBLzDCzdkC2c65Cg/MTgZkNAWY759YHHUuM6Apscs4tCD1+DjDgvOBCClwt4J/Am8BTZnZawPEkK+XxBJHg3z3fA38HBgOPBhtK5CTTd6VzbnPo+qIHgd9W5BgJP8Ya/5NTg2LbDgN+DiCWmGBmqcB1wO+CjiVgtwMn+F+AAP/35H0z+5Nz7pHgwgpMNQrlBOfcbjNbBhwZXEiBGw0Md85tMTMH/DM0rndn0IElGeXxBJDo3z3OuTXAq2aWh7/IL1Ek43flM1RwOE8y9FhPBJqZWS0AM6sOtAC+CDSqYN0OPOac2xt0IAG7Aji+UPsZ+DUJ+BNlGS0A6hUb+pELfBtQPIEKTTPY3Tm3JbTpj8A2NCtIEJTHE0OyfPfMxv9nMFEk43dlHjC3Ii9M+MLaObcWGIv/aQbgNOAd51wi/aUvMzO7G5gH1DKzNmZ2TeinuaTjnMt2zq0KN/xV3NnOuW1BxxYE59xi/PjAiwHMrD6QBnwUZFwB2gTsMbOjC23bCCwNKJ6klSx53Aq6BK3UHeNQIn/3mFltMzuh0KZhwNNBxRNpyfBdaWbNzOwKM0sJ/Tv8X+CeihwrGYaCAFwPPGJm3YEmQEbA8QTCzO4FHii2eXGiXuUrFfJL4GkzS8P3CF5e7Gr+pOGcyzOz84GRZjYHaAzckUhfJnEmofN4oZlnAK42s+ecc9lBxhQpSfDd0w0/zd5i4CtgpnMuWTsk4lUD/MWZdwNTgacqOpWnhWaIERERERGRSkj4oSAiIiIiIlVBhbWIiIiISASosBYRERERiQAV1iIiIiIiEaDCWkREREQkAlRYi4iIiIhEgAprkYMws8FmNt3MBhxiv25m9m8zuy+C5340UscSEQmCmR1nZm9FMjce5DxjzWx4NM8RLWZ2mpl9eajvGYkfybJAjMQxM+sM/Am/El4P4DjgKufc6Cif+iegexn2y8JPLp8aiZOGlmvOicSxREQAzKwH8CBwKvApsAK/0E514EXn3PgonPZ7oB4Ryo2leAFYGI0Dh/Lxhfjlu5cCX4aeqgFcBJxc0YVEQsr6PSNxQoW1xIP/AJc65+abWQrwt6o4qXNuqZltLMN+28xsTQRP3R+YHMHjiUiSc87NMrN/AZ2cc2eHt5vZWcB7Zna7c+6FCJ9zV4RzIwBmdg7wTbigdc59HOlzhDnncoAxZnYzMMk5d2+hON6vyDHN7ChgoHNujHNuiZltiFC4EgM0FERimpkdAXQAdoJfZhoYCVTVkqF5Ed6vLAYCEyN4PBERgNziG0JLb/8BeNLMmkThnBHN1WbWHN9DXdUO+OyAz4HV5TmImdUAXgdqFdocye8PCZgKa4lpzrlNwBJgnJmdEtq2EphpZvXN7AEzm2lm55rZWjP7zsxODL/ezC43s3vMbJKZ/bnQ9nQz+4OZjTGzd82sXmh7TTN73sxGmtlf8T9jlouZ/V8opkvNrIeZ/Sd0vAdCMY4xs0ah7VvN7PfFDpHmnNsbOtb/mtltZjbNzB4q9wcoInJof8UPCbkEwMzamtkoM3vFzP4bKmYxsxqhXHaPmU0xs8vCBzCzG8zsPjP7ysxuOtiJSsrJZnZxKGdebWbLzOwhM2tmZn8L5dMpZjYodIgrgObAnWZ2fuhamM/N7P8VOscIM/uTmb1uZq+ZWR0zq2dmvzez2WY2yMyWmNkPZtaqIh9Y6DMZEOrRxswuNLOnzexxMxtnZnVD7+FpM/sw9H2zJvQZdwEuMbPbCh2yQei7aLuZjTjIOVPM7GEzywidw5nZe+E/H4kRzjk1tZhuQFdgJf5/9a8CR4W2p+KT1Cb8WLf6wGfAIvx/GnsB94b2bYDv9T4LOBz4V6HjTwMeDd1/EsgI3W8M7Mcnz0PF+Hd8r89R+LF4tUPbawDjgfeBZvgvhBzg90Bt4KTQ4/D+4fHj4ff9Quj+YcDvgv6zUFNTi98GXA2sOshz64BnQnn1XSA1tP1N4K3Q/SeAC0P3rwGyQ/cvB34Rut8hlKuPDT3+O/CH0P2D5eQjgC3AU8DJwGnAY8DI0L53AB8WitUBrUL3awHfAVeHHg8B3iy071jgJfzQ1wtC3xcXht7nOOBPZfzsJoW+bT8izAAAIABJREFUKx4CHsWPtz6t0PMbgBah+4vw30k1gHuAVUBP4NpQ3p8Ujje0/4/A46HnLgeyDhLDr4GxofuNQ5/zGUH/vVIr2jTGWmKec26RmR0L/BG4ATjHzE5zfsx1NrDNOfcOgJn9AZgKtMF/ieSZ2W9DhxqPT+BnAUcW2v5d6LV1getDr8U5t87MfipHqA2BZ4FrnXO7QsfYa2ZrgRXOuZ9D51kLfOH8+MN5+ITfEH8R5Fn4LwGAXcAvzWwp8CL+PxUiItHg8MXmyUBL4GYzA9gI1DP/4Dp8pwDAa/ihEOBz7exiubYp8E2xc1xNCTnZObfJzLYC7zvnZgKYWRaw08ya4nt465YYtHO7zWx9oU03hY4b9grw79D2Lfjvi3dD55gdirOsJrrQGGszG4PvzAm7xDmXZWZ98AV13VD+/xlY7pybAcwIvbakY38Y+k6YVUpMJwK7Q+97nZltouiQEokBKqwlpplZKnCM8xep/CY0PONjIBPfA1DcotBtfaAF8IZz7s3QtqdCx7wTX+g+VexcJ+ET4vYKhtsB6Au0BeYX2l58jOH+/Cecywsl2fCwrCOdH/6Cc26FmV0NPAfcClwJrEdEJILMrD7+17ZF+Ly5sYT82Bjfo5oG4Jzbj+8MIPSa+51z00OPi7y2kBJzcoijaK5cBdwJLAO+AlqX8hYKv649MKHQ4+X4YS5HcWAuzqXiM5bMxf8HJGyLmT0MvIWPPVw9F39fh5LLwWuzT4AXzawm/s9iN773W2KIxlj/f/buPD6uqv7/+OszSZqkadOmTbo33feyFEvZRJAdRUBWARf0C/26sSmIAvpVREVEBFGUoPhDBUQEZVFBVGSRtaULLXQButA2bdM1bZo0y5zfH2emmUyTNstN7izv5+NxH3fm3pk7nzuTfPLJmXPPkVS3p98fgHNuAXAdMLmNx/fCJ7F3gErgnMSdZnZYbPtHzawwaXt17O6UTsb6EnAHcH8s8XWImQ0j4UIYMxuOHxFlIv5ixgfaeKqISFd8Bv8N2cP4/HhUrKUY2JMfNwH1wCkJ2w8ws74k5Vozy401VCRrKye35pfASueHVd3dgXNZjW/h3vMSwHZ8V5fAOOeiscaPabF/TP4G3OicmxPk6yS95hPAE/hvVi8AjnHObe+u15POUWEt6eCyWNEZNxzfNy5ugJn1it0+Ft+/bhvwIPBxM/uZmR0ba02owSfAIuBJMzs5dqHIZHwxPg/4rpn1MrNSfN/mQWaWY2ZfN7NPthGjxZZv4f/43NbKvuTHJ98/Dd8fMG4CcKZzrhq4rJXniIh0xF4toWb2IeA7wP8656rwrcMbgL+Z2cdiFygeHWuh/gN+9JAzzex44DPOuR34XHu5mX3HzI7BNzDEu9El5r+2cnJcYuvxDKDMzHoDRwGFZjY2tq8Bf7FfvIEl8TV+CZxtZiWx+4cB9zjnGvE1T6u52MwOj8WV38Z7FyGpZor98/FR/LeUJfguhpOAkUDvhAsjk1vF61uJv2VQrfQXiX0WxcC/8K3yHfmHQ3qIuoJIOhgJLDazv+K/JjPg8wn7G4HrzWwbvnX3MgDn3DNmdjn+68Qz8Rf/LQIwszOAn+P/UNwLXOOcc2Z2Hn4opMXAQ/iv9A7GT6hwHr6LSYuJacxsFv6inHL8V4N/Ab4Ti+dp4HBgtJmNw7ekDAHOM7MlNLfeXAwUOecqks79N+YndmjEXywkItJhsTxyIb74uweowl8A1w841Tn3MoBzrj6WH+/B57rH8H2rAa7AN0r8Fv8NXXwkjl/jc98X8F3WvuCc2xgrHA8HxpnZ79vKyWZ2Dr5f8SVmtsQ5txZ/IeWt+KL6t8AZwAHAe8Dv8K3rnzY/WtQ04CNm9jfn3KNmNgF4xMz+i//W8/pYy/q5wJDY+S3GT5Yz1MymxeL8BL7f+F8S3rc8fO4/CBgeK7wdvm/zWbFjzgeeA17Hd2/5Oz6nP4m/YPIgMzvTORc/7kP4yXp2xgrxIcA5ZrYC3xJN7PnJczasxl+ceSZQAETM7J/4Cxh7agha2Q/TZyHpzPw0sP/POTe6h17vGuecphsXEckwsW4p+c65lJygy8yuwY/QsjJ2vwD/bcMP49fmSPjUFUTSXWvdLLrnhfxXpm/0xGuJiEjPiV2ceUgKF9UFwA3A6ITNpUCNiurUoq4gkrZi/ds+if9q71zn3MPd/JLvxr6iFBGRzJLrnPtF2EG0xTlXZ35a9Z/GLmxfh++2csO+nyk9TV1BREREREQCoK4gIiIiIiIBUGEtIiIiIhIAFdYiIiIiIgHImIsXS0tL3ejRozv0nKqqKsrKyronoB6WKeei80gtmXIekNrnMnfu3E3OudQMrpt0JmdDan+OqUzvW+fpveucTH7f9pWzM6awHj16NHPmdGwm0ZkzZ3b4OakqU85F55FaMuU8ILXPxcxWhR1DT+tMzobU/hxTmd63ztN71zmZ/L7tK2erK4iIiIiISABUWIuIiIiIBCCrC+vZs2eHHUJgMuVcdB6pJVPOAzLrXLKZPsfO0fvWeXrvOidb37cenyDGzE4EbgLOT5jv/izgx0Ax8Hvgq865xti+UuB7wDYgD7jWOdeQfNyZM2e6TO3LIyKZzczmOudmhh1Ha5SzRURa2lfO7tEWazMbDPQBZiVsKwfOBM4BrgIuAa5MeNrDwL3OuWuBZcCNPRawiEgWU84WEemYHi2snXMbgMeSNo8CPuecm+uc+y3wc+DDAGZ2ODDJOfdq7LFPAJeZWVFPxSwikq2Us0VEOqbH+1g756JJ91+If4UYsxZYHbt9HLAi4bFrYzeD/8p0+3Y46SSYPh16uHuMiEiqSpWcveONHcz5wByWfHZJVw8lItJtUnEc60OBH8VuDwc2J+3fAQxLflJVVRUzZzbn7tmzZ3es43xxMbz2mi+w16+HoUM7GreISLtVVFRQUVERv1saZixd1CM5O1oXZecbO4n0yupr7kUkJO3N2SlVWJvZeKDKObcgtskBtUkPywf2uhCmrKysawORm8HBB8Nzz8H8+SqsRaRbJRaSZrYp5HA6pSdzdiTfF9TR+uh+HikiErz25uyU+dffzHKA/wW+nrB5LVCS8BgDioB13RLEwQf79bx53XJ4EZFM0dM523oZANHdKqxFJHWlTGENXA3c6pzbnbDtGWBSwv3RQD3QPWM0xQvr+fO75fAiIhmkR3N2vAuIq9c1MCKSunq8sI61YABYwrbrgHlAoZmNNbPPmdl459wcoNLMpsceegpwp3OuvluCmzHDr9ViLSICpE7OtvxYi7W6gohICuvRPtZm1hf4ZOzuxWb2M/xXid9NeugS59y9sdvnA98ys5X48VSv77YAp0yBwkJ45x2oqoKysm57KRGRVJdKOVst1iKSDnp85sXuEtgsXscdB88+C48+Ch//eNePJyKyH6k882J36WjOrt9Uz0tlL5E7IJcPbv5gN0YmIrJvKTPzYlr40If8+vnnw41DRET2iI8KohZrEUllKqyTHX20Xz/3XLhxiIjIHvGuIOpjLSKpTIV1siOOgPx8fwHj+vVhRyMiIoDl+YsXXb0jU7owikjmUWGdrHdvOOEEf/vJJ8ONRUREALCINRfXDSqsRSQ1qbBuzRln+PXjj4cbh4iI7LFnkhh1BxGRFKXCujWnnebX//gHbN0abiwiIgJoyD0RSX0qrFszdKjvDrJ7NzzwQNjRiIgIzSODaFpzEUlVKqzb8j//49e//nW4cYiICNDcFUQt1iKSqlRYt+XMM2HAAD86yIsvhh2NiEjWixT4P1lNtU0hRyIi0joV1m0pKIAvf9nfvvnmcGMRERFy+uQAEK1RVxARSU0qrPflssugsBD++ldYuDDsaEREslq8sG7aqRZrEUlNKqz3pbQULr3U377uunBjERHJciqsRSTVqbDen+uug759fav1M8+EHY2ISNZSYS0iqU6F9f4MHgzXX+9vX3EF1NWFG4+ISJZSYS0iqU6FdXtccQVMnAhvvw3/939hRyMikpVUWItIqlNh3R4FBXDffRCJwK23wvPPhx2RiEjWUWEtIqlOhXV7HX44XHstRKNw3nmwdm3YEYmIZBUV1iKS6lRYd8SNN8Jxx8GGDXD22epvLSLSg+KFdWN1Y8iRiIi0ToV1R+Tmwh/+AOXl8OqrcNFF0KSWExGRnpA3MA+Axi0qrEUkNamw7qiyMnjySejXDx59FL74RXAu7KhERDJeXqkvrBs2NYQciYhI61RYd8YBB8ATT/iLGisq4CtfUXEtItLN4i3WKqxFJFWpsO6so4+GP/4R8vLg9tv9DI3qFiIi0m3UYi0iqa7HC2szO9HMXjWz0QnbCs3sTjP7vpn9wsyK27MvdB/7GDz+OBQWwq9/DWedBTt2hB2ViEhgUilnJxbWTt8SikgK6tHC2swGA32AWUm77gJedM5dBzwG3N3OfeE75RR4+mno398X2UccAe+9F3ZUIiJdlmo5O9I7QqQgQrQuSnRXNKjDiogEpkcLa+fcBnyi3cPMhgEXAk/ENv0TON3MRu9rX48E3F5HHw2vvQaTJ8PixXDoob4PtohIGku1nG1m5JX5Vuv6DfVBHFJEJFA93hXEOZfczHAssMk5tyu2vxFYARyzn32pZcIEeOUV+OhHYcsWOP10uPJK2L077MhERDot1XJ2/sh8AHa/r9wqIqknFS5eHA5sTtq2Axi2n30tVFVVMXPmzD1LRUVFtwS7T/36+e4gP/qRH/P6jjt815C33+75WEQkpVVUVOzJV0Bp2PF0QKg5u6C8AIC61ZqgS0R6Tntzdm7PhdQmB9QmbcsH4pd972vfHmVlZcyZMyf46DoqEoGrr4ZjjoFPfALmzYODD4ZvfQu+9jU/ioiIZL3Zs2cze/ZsAMxsU8jhdESoOTu/PNZivVot1iLSc9qbs1OhxXotUJK0rS+wbj/7Utuhh/qi+pJLoL4ebrjBb0uF4l9EpPNCzdkFo9RiLSKpKxUK62eBYWZWCGBmeUA58Nx+9qW+4mK45x745z9hzBhYsABmzYL/+R9Yvz7s6EREOiPUnB1vsa5bpcJaRFJPGONYW/wmgHNuPfAkcHxs+wnAI865tfva14Mhd93xx8Obb/ouIrm5cO+9MHEi/PCHUKc/DiKSulItZxeOKwSgdllyjxMRkfD19DjWfYHPx+5ebGZlsdtfAM4ws+uBjwCzE562r33po6jIX9S4aBGcdpqfSObrX4fx4+EXv/DdRUREUkgq5uzC8YVYnlG3oo6mGs12KyKpxTJl9qqZM2e6lLh4sb2efhquvdZ3DwEYPRq++U341Kd0gaNIljGzuc65mWHH0ZO6krNfm/4auxbv4gNzPkDfD/QNODIRkX3bV85OhT7W2enkk+GNN+Dhh2HqVFi50ve9HjcObrtNU6OLiLShaGoRADVv1YQciYhISyqswxSJwDnnwMKFcP/9MGUKvP8+fPWrMHIkfOMbUFkZdpQiIiml99TeANQsVGEtIqlFhXUqyMmBCy/0/a+feMJPkb59O9x8M4wa5fe98AJkSLcdEZGuKD60GIDqV6tDjkREpCUV1qkkEvEXNj7/PLz8Mpx1FjQ1wYMPwoc+BAceCHfdBdX6YyIi2avvYb5f9Y45O4g2JM+4LiISHhXWqerww+GRR2DFCrj+ehg0yLdof+lLMHw4fO5z8NxzENUfFRHJLr1Ke1E4oZBobZSaN9UdRERShwrrVFdeDjfd5Pte/+EPfqr0nTvhN7+BY4+FsWP9aCLLl4cdqYhIjyk+3HcH2f7C9pAjERFppsI6XfTqBeefD//5DyxZ4luxy8th1SpfeE+cCEccAbffDmvTa/4cEZGOKjnez5y+5ektIUciItJMhXU6mjTJF9MrVsC//w0XXwx9+sArr8BVV8GIEf4CyDvvhHXrwo5WRCRwJSf7wnrbs9toqtVEMSKSGlRYp7NIBD78Yd8tZP1631XkrLOgoABefBEuv9wX2cccAz//uX+MiEgGyB+ST58ZfYjWRdn23LawwxERAVRYZ46iIt9V5JFHYONGeOABOPNM34Xk+efhy1+GYcPgqKPglltg6dKwIxYR6ZKBpw0EoOpPVSFHIiLiqbDORH37wgUXwJ//7Ivs3/8eTj/dT5X+0kt+KvXJk/1y7bV+m0YXEZE0M+iCQYAvrKO7lcNEJHwqrDNdcTFcdBE89hhs3gx/+hN86lNQUuJbrW+5xbdiDx3qp1R/+GHYujXsqEVE9qtoShF9ZvShaXsTm/+2OexwRERUWGeVPn3g7LPht7/1Ldn//jdccQWMHu3v33svnHcelJbCkUfCd77jL4hs0oVBIpKaBl80GIDKX1WGHImIiArr7JWb6y98vP12eO89WLAAfvADPzZ2JOJnfvz2t/0QfmVlvuD+1a/g3Xc1tbqIpIzBnxlMpCDClr9tYdeyXWGHIyJZToW1gJmfLv3rX4dnn4UtW3zXkS9+EcaN811DHn4YLr0Uxo+HUaPg05/2LdwrVoQdvYhksV6lvRh0ke9rvfZOjeEvIuEylyGtjzNnznRz5swJO4zM9O678PTT8K9/+QlqtiRNyDBqlG/pji+jRvliXUTaxczmOudmhh1HTwoyZ+98cydzDpxDpDDCYe8eRv7Q/ECOKyLSmn3l7EBarM3sI0EcR1LUuHG+9fqRR6Cqyncbuf12P5xfSYmf/fG+++Czn4UxY2DkSD/03x13wJw50NAQ9hmISIJMy9l9DuhD6cdLidZGWf391WGHIyJZrEMt1mZ2CnADMITmojwCDHbOFQYfXvupxTok0SgsXOi7kPznP37M7G1JkzX07g2HHeYviDzqKN9vu3//UMIVSUXd1WKdTTl75yLfam25xqwlsygcG+rpiUgG21fOzu3gsX4H/ACYB8QHDTXgws6HJ2ktEoGDD/bLVVf5QnvJEvjvf/3y0kuwfLkvvJ991j/HDKZOhVmz4NBD/fqAA/xkNiISpKzJ2X2m92HwRYPZ8PsNLL98OQc8cQCmLmki0sM62mL9qnPusFa2D3HOhTpftlqsU9jGjb7AfuklX2zPmQP19S0fk5/vi/N4oX3ooTBxoi/cRTJcN7ZYZ1XO3l25m9cmv0ZTdRPTHp1G2cfLAj2+iAjsO2d3tLA+BBjrnPtT0vbZzrmKroXZNSqs00hdHcyfD6+/Dq+95tetTbFeXAwzZ/oie8YMX3hPmKBiWzJONxbWWZez1/58Lcu/vJxeQ3sxc8FMepXpmzARCVaQhfVlwHeBvombAeecy+likFOALwPvABOACufcfDMrBG4BdgAlwLXOuerk56uwTnPbtsHcuc2F9uuvw5o1ez+uqAgOOsgX2TNm+GXaNCgo6PmYRQLSjYV11uVs1+SY/+H5bH9hOwNPG8j0x6erS4iIBCrIPtY3ArOBOUB8Oj4DPtX58Pb4HXCGc26tmZUDTwNTgLuAp5xzD8UuxLkbuCCA15NU0r8/HH+8X+IqK32BPWeOb+GeN88X2/FuJXG5uTBlSnOr9owZvs/2wIE9fx4iqSXrcrblGFN+P4U5B81h85ObWXPbGkZ+dWRPvbyIZLmOtlg/AHzWObc7aXuX++uZWQ3wAefcEjMrAxYAM4EVQIlzbpeZ5QLbgWnOuZWJz1eLdZaID/c3b55f5s/33Uii0b0fO3QoTJ/ecpk61U/tLpJCurHFOmtzdtWfq1h81mIwmP74dEpPK+221xKR7BJkV5DTgUOAfyduBs5yzl3RxSDvA8YBpwKXAKuBfOBHzrnhCY9bFNt2X+LzVVhnsZoaePPN5lbt+fNh8WK/vTVjxvgW7cSCe9IkjUoioenGwjqrc/bKG1ey8v9WEimKMOP5GfQ9pO/+nyQish9BdgX5Jr4v3cVAvCI3/BipXUrSwBeBx4DXgR875x4xs2uAzUmP2wEM6+JrSSYpKoLDD/dLXDQKK1fCokUtlyVL/DTsK1bA4483Pz4310/XPnmyX6ZM8etJk6Bfvx4/JZGAZHXOHvXNUexauouND2xkwUkLOPjZg+lzgL6xEpHu09HC+kfA351zOxI3BjSLVyHwIDAcuN3MVuD/ENQmPS4f2Gsqv6qqKmbObP7nYfbs2cyePTuAsCQtRSIwdqxfTj+9eXt9vR9XO7HYfvNNeO89X3QvWbL3sYYObS60Ewvv4cM1dbt0WkVFBRUVewbm6K5+Clmds82MyfdOpnF7I1v+uoUFxy/goH8dpOJaRDqsvTm7o11BjnDOvdzK9s855+7tcJQtj/EU8Ann3DYz+yZwOXAl8H/OuYkJj1se2/ZA4vPVFUS6ZNcuWLbMF9Zvv91cZC9dCrt3t/6coqLmQnvChJaLZpaUDujGriDK2UBTXROLzlzE1qe3ktMvh+l/nk7Jh0t65LVFJPME2cf6d865TyVtOxf4tXOuuAsBlgJvOueGxu4bfginzwBPAWXOuVozywN24sdlXZt4DBXW0i2ammD16uZCO7Horqpq+3kDBzYX2ePHt7ytoluSdGNhrZwd01TbxNufeptNj2zC8oxJv57EkE8N6bHXF5HMEWQf6wFm9k3n3HfNrDfwU+BcfB+7rtgC1JnZ8ITkuxlYCDwJHB9bnwA8kpygRbpNTo6/2HHMGDj11Jb7Nm/2LdpLlvjuJcuXwzvv+PXmzX555ZW9j1la2nrBPW4clKgVTQKlnB2TU5jDtIem8e7V77Lm9jUs+fQSql+uZtxt48gp6NKQ3iIie3SoxRrAzA4HLgJOBtYDnwa2O+e2dikQs4PwF8PMBQYDzzvnnjOzEvxkAyvxF9x8wzm3M/n5arGWlOGcH4M7XmQnFt3vvOO7nbSlf//mvuFjxjTfHjsWyss1ckmG6q4W69ixlbOTrLt7HcsvX46rd/Q5pA9TH5hK70m9Q4lFRNJPp7uCmNmH2tg1C/gwcCsQBc50zl3V1UC7QoW1pAXnYN26lkV3/PZ777U9RCD4CzJHjGhZbCcW32VlupgyTQVVWCtnt9+OuTtYfO5i6lbUESmIMPrG0Yz8ykgsR79DIrJvXSmsX8HPpLUNn4xbfRgw1DmX39VAuyLsJC3SZc75ftsrVvgiO3FZsQLef7/1iXDiiopaFtrxLiyjRsHo0VDc6S610s0CLKyVszugYVsD71z5Dhvu2wBA31l9mfCzCRQfqt8VEWlbV/pY3wjMcc5t3M8LnNPZ4EQkxgwGDfLLYYftvb++3l9IGS+0k4vvbduahxBsTUlJc5EdXxLv66LKTKCc3QF5/fOY8v+mMOj8QSy9dCk7XtvBG7PeYPCnBzP2B2PJHxbq/x4ikoY63Mc6VaVC64dIqLZu3bvgXrXKT5SzahXUJg8vnKS4uPWCO35/wAB1Nekm3dnHOlWlWs5urG5k1fdXseYna3D1jkjvCMO/NJyRV4+k1yBd2yAizQIbbi+VpVqSFkkp8W4m8SJ75cq9b++rfzdAnz5tF92jR/vRTlR4d4oK69RR+24t717zLpv+vAmASO8Iw74wjJFfGakWbBEBVFiLyP4454cHbKvoXrkSduzY5yHo3btl0R2/HV8PHqzCuw0qrFNP9ZxqVt24is1P+BnaLdcoO6+MEZePoPgw9cEWyWYqrEWka5zzfbhbK7rjt7dt2/cx8vN9kZ1ccMfXQ4f6ccOzkArr1LVj7g5W/WCVb8GOXQ7ad1Zfhl4ylEHnDSK3X0engxCRdNethbWZzQC+CszDTwSwsksH7KR0SdIiGWv79tYL7vh68+Z9Pz8314/V3VbhPXw45OV180mEoycLa+XszqlbVcfau9ZSeU8ljVsbAYgURCg9s5TBnxlMyQklRHIjIUcpIj2huwvru4EznHNDzOxC4M/Ouf1cJRW8dEvSIlln505fZCcX3PH1hg37fn58HO+2Cu+RI32reBrq4cJaObsLmmqaqPpTFevvW8+2Z5u/pckdkEvp6aWUfryUkhNLyCnMzm9fRLJBdxfWecAg59xaM4sAfZxz1V06aCeka5IWkZjaWj+cYFuF97p1vktKW8x8d5K2Cu/yct8PPAX1cGGtnB2QulV1rP/dejbev5FdS5pnVI0URRhw8gAGnDSAkhNLKBxbGGKUIhK0QAtrMzsVGOCcu9/MZgE7nXNvBRBnl2RCkhaRfaiv95PktFV4r1kDTU37PsagQW0X3uPGQWE4BVA3T2munN0Dat6uYdOfN1H1aBU757acwb1gbAElJ5ZQclwJ/Y7qR/7w9PxmRUS8wAprM/secAXwlHPunNi2m4GFzrkHggi2szItSYtIBzU2wtq1bRfeq1dDQ0PbzzfzrdqTJrVcJk/2/bu7cUST7iqslbPDUbe6ji1PbWHrM1vZ+s+tNG5rbLE/f1Q+/Y7sR/GRxfQ7oh9FBxQR6aX+2SLpIsjC+lng48BnnHN3xLYdA/zOOVceRLCdlclJWtJQfT28+CLMnev7Dm/eDNXVzV0Z4mszPxV5UZEfJzq+JN9vbVthoYav64hoFCorWy+84xPrNDa2/tyiIpg4EQ48EGbM8MvBBwc2TXw3FtbK2SFzTY4dc3ew9ZmtbHt+G9WvVNNU3fKbFcszek/tTZ+D+zQvB/Yhb0BmXqybSqINURo2N9C4uZGGzQ00bGrw97c3Et0VJVobpam2iWhttHlpiLY4hsXysPUycopyyCnKIVIUIacoh9x+ueQNyqPXoF70GtyLvEF55A3MwyLK3emsK1OaJ3vGObfNzBKrcU2NK5JoyxY46ihYsqR7XydelHekGG/P/fz8zCzYIxHf8jx8OBx55N77Gxp8gb10qf/sli5tXqqqYN48v9x3X/Nzxo+HD34QTjoJTjzRT5KTWpSzQ2Y5RvGsYopnFTPq+lG4JkfNWzVUv1TN9pe2U/1KNbXLa6lZUEPNgho23Nd8EW9eWR69J/WmcGIhvSf13nO7YFQBOb11cWQi5xzRXb5I3rNsSiiYE4vmhG3J/+T0BOtlFIwuoHBcIYXjCikYV0DRlCKKDixEI360AAAgAElEQVSi15Beewp1SU8dbbH+BDAWGA68AJyLbw25zDn3826JsJ2ypfVD0sDTT8Mpp/jbs2bB2WfDwIG+dTMSaS5azXwr6q5dfsSM+FJT0777+5uivLNycpoL7sR1V7cVFKTvcHlbt/pie/58eOMNX2C/+ab/ZiIuJwe++lW4+eYO/2PSjS3WytlpoHFnIzVv1rBz/s49S82iGqK7om0+J680j/zyfApGFfh1eQH5w/PJG+xbR/MG5ZE3ID1bRqMNURq3NNKwJVYIb2lssW7YklQwx4rlaF3b71ebIpA3II/cgbnklfrW5LyBeeT2zyXSO0JOYQ6RwkiL25ZnEH9bXfM6Wh8lWhOlqaaJpp1NNNU00bi9kYaNDdRvrKdhg183bmnjmzH8P1N9DvLfWhQfXkzxUcXkD1Gf/FQT9MWLhwMXA+VAJfCgc+6fXQ2yq5SkJWXs2gVHH+0LMIAzz4TrroNDDw32dZqafJHd3kK8vfcTi8Wg5eT4LiwFBX7dntuJ23r18sV5fJ18u619OTm+2E1e4v/oxBfw3UEaGpqX+P36ej/7ZHzZsgX++U94/vmW5/j2275fdgd088WLytlpyEUdu9fuZtfSXdQuq2XX0l3+9vJadr+/G9fQjr/dEV+o9RrUi7zSPHKKc8gtzt2zzu3nb+f0zSGnMAfrZUTyI0TyI1i+EekVu90roTh3SWt8a7Hb7Yjuju5Zku9Ha6M07WiiaUcTjdWNzbd3NLbY3rjF3+8MyzdfGCcWyANz97ktt19uj//z0VTTRO17tdS9V0ftu7XUvlNLzeIadi7YSdP2vc+9YFwB/Y7qR/8P9afkpBIKRhb0aLyyt+4ebi8HONY5968uHaiLlKQlpezaBTfeCHfcAXV1ftusWfD5z8P556fssG+ALyITC/Z9rdv7mJoa38Ie7USLUjr55CfhN7/xk910QA8Pt6ecneZc1FG/oZ7dq3dTt6qOutV17F69m91rd9NQFWsd3diwZyKbtBNvRR7gC+D4ek/LcsK+xKI50juS1t0onHPsXr2bnQt2smPuDqpfqvZ98ne2LLZ7T+3NgFMGMODkAfT7UD9yCtQtqKd1urA2s8eBPvs5fimwyTl3XOdD7DolaUlJ69fDbbdBRYWfmRCgXz+48EJfhB1xRGb2Z26Nc771t7bWL3V1e99ubVvi7fr65tbjxFblxPut3W5q8q/f2hKNNt8GXxQntoDH7/fq5bu29O3bvJSV+UlrRo70FzP22V+6bF1QhbVytiSKNkR9v+KNsf7G1Y00VTc1r7c37rkdrYsSrW/Z0uzqY7froy0LVktaQ+st3fnNLeCRwgg5fWMt5n19K3lu34Tbse25A3LJLe75VuRUFW2MUrOwhu0vbmfrv7ey7V/bWhTaOX1yGHjaQErPLmXgqQPJKVKR3RO6UljfhL/AcSkQBS7FT4ObmA3H4MdIvTywiDtBSVpS2q5d8Mc/wi9/Ca++2rx97FhfYH/ykzBhQnjxSagCLKyVs0UyWLQ+SvXL1Wx5agtbntrCzvnNY6ZHCiMMOHUAg84fxMDTB6oluxt1pbDuF3vMttj9G51z32rlcY85584IKuDOUJKWtLFwIfzud/DAA342wbjDDvMF9nnn+YlMJGsEWFgrZ4tkkdoVtVQ9UsWmRzZR/UrzBKo5/XIYdP4ghnxmCMVHFKd1F5lUtK+cvc8R6Z1z2+MJOmZwrH9e4sGPBWZ1OUqRbHHggfCjH/kJS555Bj7zGd+F4NVX4bLLYNgw+MhH4P77fT9lkXZSzhbJLoVjCim/upxDXj6Ew1cfzrifjKPPIX1o2t5EZUUl846ax2uTXmPV91axu3J32OFmhY4Ot3ci8BPgRWAHMAU4GfiWc+4H3RJhO6n1Q9Larl3w2GO+mH766eaJSnr39qOKXHSRHyc5XYerk33qxuH2lLNFstDORTvZcN8GNvx+A/Xr/UhPlmuUfryUYV8cRv9j+qsVuwuCHm5vGPBF4ECgGnjYOfdYl6PsIiVpyRhVVfDww77Ifuml5u2lpX5EkYsugsMPz56LHrNANw+3p5wtkqWijVG2PrOVyl9VsumxTRC77rH3lN4M+8Iwhnx6CLn9OjpXoHT3cHsRYIZzbm6XDtR8vMOBDwLv4FtZaoBb8K0tJcC1zrnq5OcpSUtGWrHC98W+/34/PnLc2LF+ZJGLLurwmMmSenp4uD3lbJEsVLemjsp7Kqm8p5L6St+KndMnh6GXDGX4FcMpHF0YcoTpI7DC2szubWVzGdDknDuzk/ElHv+zwDjn3A0J234DPOWce8jMTgE+45y7IPm5StKS0Zzzs/7dfz88+GDLix4POcQX2J/4hO+fLWmnG7uCKGeLSAvRhiibHtvEurvWse3Z2CUZESg7t4yRXx1J8aHF4QaYBoIsrF8A/kWLOZc4CnjBOffdLgZ5DPBN4EQXCyr2FeYKoMQ5t8vMcoHtwDTn3MrE5ytJS9ZoaoLnnvNF9p/+BNWxxkAzOOkk+MIX4KMf7fAkJRKebiyslbNFpE075u9gzY/XsPEPG3GNPk30O7ofI68eycDTBmo88TYEWVhPcs4tTdo2GbjYOff1rgYJvIQfcn488F1gFPAj59zwhMctim27L/H5StKSlerq4K9/9UX2X//aPB35iBEwezZccgkMHRpujLJf3VhYK2eLyH7Vralj7U/Xsu7udTRV+47YRdOLKL++nEHnDsJyVGAn6vRwe8mSE3TMbuCSzgQWZ2aTgBnAPc65L+NbWJ4CRgKbkx6+A9jr++6qqipmzpy5Z6moqOhKSCLpoaAAzj4bHn0UKivhxz/2E82sWQPf+haUl8MFF8CCBWFHKkkqKir25Cv8bIiBU84WkfYoGFHAuFvGccSaIxh32zjyR+RTs6iGty94m9emvsb6+9YTbYiGHWao2puzO9pivYKWXynmAIOB3znnLu1krJjZWUCFc640dr8QqAJuBM52zh2W8Ng3gAecc7cmHkOtHyIx0Sj861/wi1/A44/7riPgx8a+7jo46qhw45O9dGOLtXK2iHRYdHeU9fetZ/UPVlO3sg6AgtEFlH+9nCEXDyGS36F22Yyzr5zd0U6YzwAP0Jyoo8AG59yyLsQXj2NPLM65WjNbjv8jUJL02L7AOkSkdZGIH/P6xBP9JDS33QYVFfC3v/nltNPg1lth0qSwI5Xup5wtIh0WyY8wbPYwhnx2CBsf2Miq76+idlktyz6/jFU3rWLUDaMY8rkhRPKyu8BuTUffkd875/7jnHsutrwAmJmN72IcC4F+ZjY4YVsjPhkPi7WGYGZ5QDnwXBdfTyQ7lJfD7bfDqlVw/fV+hscnn4Tp0+Gaa3wfbclkytki0mmRvAhDPjOEWW/NYsqDU+g9rTe71+xm2eeX8drk11j/+/W4pq4N25xpOlpYn9LKtg3Aw10Jwjm3BPg7cA6AmfUHegH3A08Cx8ceegLwiHNubVdeTyTrlJXBTTfB8uX+gsamJt9qPXMmLFwYdnTSfZSzRaTLLMcY/InBHLrwUKY+NJXCSYXUvVfHkk8t4fUDX6fq0Sq6Oi9KpmhXYW1mXzazbcC1ZtaUuOAvVNkQQCyfBo40s6uA/wMucM41Al8AzjCz64GPALMDeC2R7DRkCNxzD7zyCkycCIsX+z7XzzwTdmQSIOVsEekOFjEGnTeIQxcdyqR7J5E/Kp9db+1i8dmLeWPWG2x5ekvWF9jtvnjRzA4EPgP8NGlXvXOuMujAOkoXwoh0UE0NXHqpn3AmLw+efVYXNoakOy5eVM4Wke4W3R2l8leVrLppFfXr/XCv/Y/rz7hbx9F3Rt+Qo+s+QY5jneeca2hle7lzbnUXYuwyJWmRTohG4YtfhLvv9uNdL10KfTM3GaaqbhwVRDlbRLpd064m1v58Lat/sJrGrY1gMPhTgxlz0xgKRhaEHV7gOj0qiJldDfzNOfdWbNPXzCwn+WHAh4FjuxqoiPSwSATuvBPmzYPXXoO77oJrrw07Kukk5WwRCUNO7xzKryln6CVDWfW9Vay9cy0bfruBqj9WMeIrIyi/tpzc4uyYDXh/fazHAP0T7h8CHAmMje2LL0O6JToR6X55efCNb/jbTzwRbizSVcrZIhKavJI8xt86nllvz6Ls/DKidVFWf381r45/lbW/WEu0MfMnmdnnvw/OuS8lbfoesMg5V5+40cymBh2YiPSg8bHR16qqwo1DukQ5W0RSQeHYQqb9YRrbr9zOu1e/S/V/q1n+xeWs/dlaJvx0AiXHJw93nzk6OtzeAa0k6EnAuOBCEpEe98c/+vWUKeHGIUFTzhaR0PQ7vB8zXpjBtEemUTCugF1v7WLBCQtYfN5i6lZn5jwK7erwYmZH4/vlnRibIjfRAPxV5/oOWSTdOAd33OHHuAa48spw45FAKGeLSKowM8rOKmPARwaw5rY1rPreKqoermLzXzcz6rpRjPjqCHIKki8FSV/t7Um+Dfgtvm9e8nhc9cCvggxKRHrA4sVw2WV+mD2AW26BY48NNSQJjHK2iKSUnIIcRl03isGfGsy7V79L1R+rWHHDCip/U8n428dTelpp2CEGol2FtXPuTTM7FjjXOaeELJLOFi2C733Pd/+IRmHgQPjlL+Gcc8KOTAKinC0iqapgZAHTHprG1v/dyvLLlrPrrV0s+tgiSj9eyoSfTSB/WH7YIXZJu/tYO+e2t5WgzeyDwYUkIoFraoLHH4dTT4UDDoA//AFycvwY1suWqajOQMrZIpLKSo4rYeb8mYz7yThy+uaw6c+beG3Ka6z95VpcNH1nb9zfONaPA332cwwDpgGDggpKRAKyZIkvou+9F95/328rKIBLLoFrroHy8nDjk0ApZ4tIOonkRRh55UjKzilj+ReXs/mJzSz/wnI23r+RiRUTKZpSFHaIHba/riBL8X311gJt/ftgwGlBBiUineQcvPWWb51+6CFYsKB53/jx8PnPw8UX++4fkomUs0Uk7RSMKGD6Y9Op+lMVyy9bzvYXtzPn4DmMum4U5deVE8nr6CB24dlfYf1DYIdzbve+HmRmurpcJCzbtsF//gN//7tf4i3TAP36wdlnw4UXwoc/7GdalEymnC0iacnMGHTuIEpOKOG9r71H5a8qWfntlWx6fBNTfjeFoqnp0Xq9vwliNrXzOB9HV5mL9IwNG+CFF+D55/2ycKFvqY4bNAhOOcUX1CefDPnpfSGItJ9ytoiku7ySPCbdM4lBFw1iycVL2PnGTuYcMoexN49lxOUjsIiFHeI+dWnidjMrB74EfBklaZHg1dTAvHnw+uvNyzvvtHxMXh7MmuWL6VNPhRkz1DItrVLOFpF0UXJsCYcuPJR3rnqH9feu592r3mXz45uZ/P8mU1BeEHZ4bepwYW1mOcDpwKXAicAGfJ8+EemKHTv8UHjz58OcOb6IXrzYD4mXqHdvOPJIOPpo+NCH4LDDoLAwnJgl5Slni0i6yi3OZfKvJ1N6eilLL13Ktme38fqBrzPlvimUnpGa4163u7A2szHAJcBngYHAKuBjwNPAEd0SnUgmikZh5UrfhWPBAr8sXAjvvrv3Y3Ny4KCD4NBDm5fp030rtcg+KGeLSKYoPaOU4iOKWXrpUjY/vplFZy5ixFdGMPbmsSl3YeN+C2szOxff0nEc8A5wG3AfcJ1z7qnYw17qtghF0pVzsGaNH6Xj7bf9+q23fBG9Y8fej8/Lg6lT4cAD4QMf8EX0wQf7FmqRdlLOFpFM1GtQL6b/ZTprfrKG9659jzW3raH65Wqm/mFqSnUNaU+LdS1+CtyXgEucc8sAzCx9R+8WCVJTE6xY0Vw8J6537mz9OUOG+JboAw9sXk+erJZoCYJytohkJDNj5FdGUnxEMW+d9xbVL1czd+Zcpv95Ov2O6hd2eEA7Cmvn3JPAk7GLXi41s9HAY8CeCsDMzDmnpC2ZbfNmP0vhsmWwfLlfL13ql91tjG5WVuZboadM8eupU/3Mh4M0N4d0D+VsEcl0/Y7ox8z5M3nrwrfY+o+tzD9uPpN+NYkhnxoSdmjt72PtnFsNfNPMcoEzgRFm9l3g18BJQEX3hCjSg3bs8EVzvHCOr5ctg61b237eiBHNxXPiujQ1L66QzKecLSKZLG9gHgf89QDevepd1v5sLUs+vYRdS3Yx5qYxmIU3JF+HRwVxzjUCfwL+ZGbjgavw/fmUpCX1OQdbtsB77zUv777bXExXVrb93L59YcIEmDix5XrKFCgu7rlzEOkA5WwRyVSR3AgT7pxA78m9WX7FclZ/fzUNVQ1M/MVELCec4rpL41g7594BrjCzuQHFI9J1DQ2wenVz0ZxcRFdXt/3c/Hw/9Xe8cE4sogcPhhD/CxbpKuVsEclEw780nIIxBSw+ezGV91TStLOJyfdNDmXEkC4V1nHOud8GcRwz6wW8DlzhnPuPmRUCtwA7gBLgWufcPqoiyRpbt7YslhOL59Wr/QWFbenTB8aNg7FjW64nTfJdOnJyeu48REIQVM4G5W0RSQ0DPzKQA586kDdPe5OND24kWh9l2kPTerzlOpDCOkDXAKMT7t8FPOWce8jMTgHuBi4IIzDpYYmtzq0t2/Yxv4UZjBzZsnCOL+PGwcCBankWCY7ytoikhP7H9Oegfx3EgpMWsOmRTSz7/DImVkzs0T7XKVNYm9mRQCWwNXZ/GHAhfvpdgH8Cj5jZaOfcylCClOA450fZaK1oXrHCF9XJMw4mKiqCMWP2bnkeOxZGj/ZdOkSkWylvi0iqKZ5VzAFPHsDCExdS+atK8gbnMfamsT32+ilRWJtZEXCOc+4rZvat2OZjgU3OuV3gL8AxsxXAMcDKUAKVjtm9G1atarvVubVJUuISW51bW8rK1OosEiLlbRFJVf0/2J9pf5rGm2e8yervrabPQX0YdG7PDHObEoU1cC1wc9K24cDmpG07gGGtHaCqqoqZM2fuuT979mxmz54dZIySzDnYuLFlS3Ni4bxmjX9MW/r23burRnwpL1ers2S0iooKKir2DMyRjuMydilvK2eLSHca+NGBjP/xeN658h2WfHYJRdOKKJpa1OnjtTdnh15Ym9nJwBzn3MakXQ4/g1iifKChteOUlZUxZ86cbogwy9XWwsqVbbc679rV9nNzcvbd6jxggFqdJWslFpJmtinkcDokiLytnC0i3W345cOpfq2ajQ9s5O1Pvs0hrx7S6ZFC2puzQy+sgauBGQkdy0vws4TdFrudqC+wrudCyxLV1S0nQkksnNft5+0uKWm7cB45UlN0i2Qm5W0RSXlmxsS7J1L9UjU75+3k/R+9z6jrRnXra6ZCYX0hvkUj7mXgK8CrwDVmVuicqzWzPKAceC6EGNPf7t2+UI4Xz4nL+vVtPy83118MOGbM3oXzmDG+sBaRbKO8LSJpIbdPLhPvmcjCExey8jsrGXzRYApGFXTf63XbkdvJOVeVeN/MmoAq59xqM3sSOB54EjgBeMQ5tzaEMNODc76F+a239i6eV65se5SN/PyWMwlOmNBcPA8f7otrEZEY5W0RSScDThjAoAsGsfHBjaz89kom/2Zyt71WqldMXwBuMbODgCGArmwBX0CvXw+LF++9bN/e+nMiEV8ox2cTTFxGjvT7RUS6TnlbRFLOmO+OoerhKtb/dj3l15XTe0LvbnmdlCusnXOjE25vBS4NL5oUsGMHLFgA8+bBokW+eH7rLT/zYGsGDIBp02Dy5JbF85gxGmVDRLqF8raIpLrCcYUM/vRg1t+7nnV3rWP8T8Z3y+ukXGGd1TZv9gX0G280r5cvb33Iuv79fQGdvAwerJE2RERERJIM/9Jw1t+7nsrfVDLm+2PIKcwJ/DVUWIdl925fPL/0kl9ef93PNpgsLw+mT4cZM+CAA5oL6KFDVUCLiIiItFPfQ/rS55A+7HxjJ1uf2Urp6cFPIaDCuqds3gzPP99cSM+d64vrRL17w0EHwSGH+EL6kEN8Ed2rVzgxi4iIiGSQsrPK2PnGTjb9ZZMK67RSWwv//S/885/wzDO+dTq5S8fUqXDkkX457DCYNMlPqiIiIiIigRt4xkBW3LCCLU9v6Zbjq7AOUmUlPPYY/OUv8J//tGyR7tULjjoKPvjB5kJaY0CLiIiI9JiiqUXk9s+lfl09de/XUTAy2DGtVVh31bp18OCD8Oij8PLLLVulZ8yAE07wywc/6Lt6iIiIiEgoLGL0PawvW5/eSvUr1SqsU8Lu3fDEE/Cb38BTTzVPvJKfDyedBB//OJx2GpSVhRuniIiIiLTQ5+A+bH16K7ve2hX4sVVYd8TWrfDLX8Idd8CGDX5bXp4vpD/xCTjlFOjTJ9wYRURERKRN8clhat+pDfzYKqzbo7oafvAD+NnPYOdOv236dLjkErjoIigN/qpSEREREQle4YRCAHYtV4t1z4pG4Ve/ghtugKoqv+2EE+BrX/NrjSMtIiIiklbyR/qZqOsr6wM/tgrrtrz/Pnz2s/Cvf/n7Rx0FP/6xH81DRERERNJSXmkeAA1VDYEfW4V1a55/3veb3rLFd/O48044/3y1UIuIiIikuZw+OVi+Ea2N0lTTRE5RcHOIRAI7UqZ4+GHfzWPLFjj1VFi0yF+YqKJaREREJO2ZGb3K/KzWDZuCbbVWYZ3oqafgwguhoQEuv9wPqTd4cNhRiYiIiEiAcvr4VuqmnU2BHleFddw778C550JjI1xzjR9ST9OLi4iIiGScSKEvgZtqVVgHr6kJPvlJP5TeOefAD38YdkQiIiIi0k3ihXW0NhrscQM9Wrq6/3549VUYPhwqKtSfWkRERCSD5fT2vRJUWAetqQm+/W1/+/vfh5KSUMMRERERke6lFuvu8o9/wIoVMGaMn0VRRERERDJaJD9WWO9WYR2s3/3Ory+9VBcrioiIiGSDWMnnoi7Qw2Z3Yd3U5FuswV+0KCIiIiIZzyKx6+mCHRQkdQprMzvLzFaY2WYzu8PMcmPbS83sbjP7oZndZmZ5gb3oggWwebPvBjJ+fGCHFRHJdKHkbBGRgFiOL6wzssXazMqBM4FzgKuAS4ArY7sfBu51zl0LLANuDOyF583z68MP10ggIiLtFFrOFhEJSqwCdk3BFta5gR6t80YBn3PONQJzzWw68GEzexGY5Jx7Nfa4J4ClZnaTc66my686f75fH3xwlw8lIpJFwsnZIiIBibdYZ2RXEOfcC7EEHbcWWA0cB6xIeNza2M2ZgbzwsmV+PXVqIIcTEckGoeVsEZGgxFusM7ErSCsOBX4JDAc2J+3bAQwL5FXWrPHr8vJADicikqV6JmeLiARkTx/rDO0KsoeZjQeqnHMLzMwBtUkPyQcakp9XVVXFzJnNjSKzZ89m9uzZ+36x99/36xEjuhSziEhHVVRUUFFREb9bGmYsXdGjOVtEJCB7uoK0cxjr9ubslCqszSwH+F/g67FNa4GJCfsNKALWJT+3rKyMOXPmtP/Fampgxw7Iz9dsiyLS4xILSTPbFHI4ndKjOVtEJEgdvHixvTk71bqCXA3c6pzbHbv/DDApYf9ooB7oejbevt2v+/fXiCAiIp3TczlbRCRAe8axDnbixdQprM3sOmAeUGhmY83sc8A2oDJ2xTnAKcCdzrn6Lr9gdbVfFxd3+VAiItmmx3O2iEiQYnW1cxnYx9rMbgC+m7R5iXPuXjM7H/iWma0E+gDXB/Ki8Rbrfv0COZyISLYIJWeLiASpmzorpERh7Zy7CbipjX2rgP8J/EXVYi0i0imh5GwRke4QbIN16nQF6XF1dX5dWBhuHCIiIiLSo6ybrq/L3sK6ITb6U15euHGIiIiISDjUYh2Q+ti1NL16hRuHiIiIiPSseIO1CuuAqMVaREREJDt108WL2VtYq8VaREREJKsFPdxe9hbWarEWERERyU5qsQ6YCmsRERGR7KY+1gFpbPRrFdYiIiIiWWXPcHsqrAPS1OTXkex9C0RERESykrqCBCxeWOfkhBuHiIiIiIRDLdYB6d0bhg6Ffv3CjkREREREelBuv1x6De1FpCjYUjg30KOlk8su84uIiIiIZJVR149i1PWjAj9u9rZYi4iIiIgESIW1iIiIiEgAVFiLiIiIiARAhbWIiIiISACyurCuqKgIO4TAZMq56DxSS6acB2TWuWQzfY6do/et8/TedU62vm8qrDNEppyLziO1ZMp5QGadSzbT59g5et86T+9d52Tr+5bVhbWIiIiISFDMuYCnnAmJmVUBqzr4tFJgUzeEE4ZMORedR2rJlPOA1D6XUc65srCD6EmdzNmQ2p9jKtP71nl67zonk9+3NnN2xhTWIiIiIiJhUlcQEREREZEAqLAWEREREQmACmsRERERkQBkZWFtZoVmdqeZfd/MfmFmxWHH1BFmdqKZvWpmoxO2pdU5mdlZZrbCzDab2R1mlhvbXmpmd5vZD83sNjPLCzvW/TGzY8xsqZltNbPbE7an1WcSZ2a9zGyBmR0bu59252FmN5mZiy0LYtvS7jykmT6/zmvt90Falwl/X8PS2nsX255VP39ZWVgDdwEvOueuAx4D7g45nnYzs8FAH2BW0q60OSczKwfOBM4BrgIuAa6M7X4YuNc5dy2wDLgxlCDbycz6AUcBhwOfBL5kZifEdqfNZ5LkGmB0wv20Og8zK8JfjX5ibDkntiutzkP2os+vE/bx+yBJMuHva1jaeu+y8ecv60YFMbNhwAqgxDm3K9ZSuh2Y5pxbGWpw7WRmEaAJGOOcW5lu52RmRwMvO+caY/dvAaYB3wUedc4Ni20fDiwFBjvnasKKd1/MrBCoc7FfJDObC1yNjzttPpM4MzsSmAx8C7gY/89NWp2HmcX/Sfulc64uti2tfkekJX1+ndfa74O0Ld3/voYp+b2Lbcu6n79sbLE+FtjknNsFECvuVgDHhBlURzjnokmbjiWNzsk590K8qI5ZC6wGjsPHHX/c2tjNmT0YXoc452oTiuoi4G3gP1Y6l+4AAA0PSURBVKTZZwJ74j/HOXdvwuZjSaPzMDMDPgH8GNhgZhfEdh1LGp2H7OVY9Pl12D5+H6QN6f73NUzJ7122/vxlY2E9HNictG0HMCyEWIKS7ud0KPBL0vg8Yq3wfwfygN6k57lcC9yctC2tzsN5h+O/erwTuN/MTiXNzkP2os+vE/bx+yDtp5+9TsrWn7/csAMIgQNqk7blAw0hxBKUtD0nMxsPVDnnFphZ2p4H8A7w/4BbgB8B75FG52JmJwNznHMbk3al5WfinNsK3BBrMbkSeIY0PA/ZIy1/DlNFK78Pfw85pHSin70uyrafv2xssV4LlCRt6wusCyGWoKTlOZlZDvC/wNdjm1qcR+yXsIgUPw8A51xlrAvF1fivCNPtM7kauNfMNpnZJmAk/iKdItLrPJL9FCgn/T4PaUmfXzDivw/SfvrZC05W/PxlY2H9LDAsdtEZseHcyoHnQo2qa9L1nK4GbnXO7Y7dfwaYlLB/NFAPzOnhuLpiDj4Rp9tnciFwcMKyDj9ay29Ir/NIFgXeIP0+D2lJn18w4r8P0n762QtOVvz8ZV1h7ZxbDzwJHB/bdALwSMKFcikv1pILYJCe52Rm1wHzgEIzG2tmnwO2AZVmNj32sFOAO51z9WHFuT9m1tvMZiRs+ghwR7p9Js65KufcmviCv7K7yjm3mjQ6DzMbZmYXmlkk9nvyVeD6dPs8pCV9fp3T1u9D2HGlskz4+xqW5PcuW3/+sm64PQAzK8H3hV0JDAG+4ZzbGWpQ7WRmffHjJd+FH+P5Z865qnQ6JzO7AT+0XqIlzrkpZjYKP9TbSvyYmNcnjSCSUsxsFvBXYAnwCvCac+7h2L60+UySmdlK4GLn3H/S6TzMbBq+/1418CJwu3NuSWxf2pyH7E2fX8ft6/dB9pYJf1/D0tp7BwwiC3/+srKwFhEREREJWtZ1BRERERER6Q4qrEVEREREAqDCWkREREQkACqsRUREREQCoMJaRERERCQAKqxFRERERAKgwlpEREREJAAqrEUkpZhZmZk9bmarzOzrYccjIiL7przdTBPEiEhKMbOTgWeBYvyMlhOdc1vCjUpERNqivN1MhbWIpCwzex44LpWntRcRkWbZnrfVFUREUpKZjQT+ka3JWUQk3Shvq7CWbmBmXzOzn4cdR1wQ8ZjZCWb2spkdG1BYGcXMbg34eAacB9wc5HFFpCXl664zs9Fmdo+Z3R3Q8XLM7FIzW9WOx37UzF7t7LkGmbuVtz0V1mnGzKaY2aNm9iszW2Bmzsw+GXZcSf4GPBh2EAmCiOd94KAAYklJZvYRM9tkZpVmdoSZjYz9nDkz+5mZ9Y89bpaZLTKzH5tZYWzbYGB9wCFdCPzaOdcYO75I2lG+7pR0zNcbY+v8gI4XAbYB5e147BxgSls7zWyamT0U+9m7wcyGJOwLOncrb6PCOh39BfiOc+4SYAbw2/gOMzvQzI4JI6jE13bOLXLOvRhGHK3pSDxtvYfOuaXApiDjCvLz6uqxnHN/A+4CGp1zLzvn3gcuBeqBV51z22KPew2fyK9xztXGnv4x4MkunUACM/sm8G3gn2a2BJgU1LFFepjydQelY752zu0C1gX1Ws65BmBuOx+7AV+Et7V/MVARu3urcy6xkA4sdytvN1NhnUbMbAAwEagBcM5FgW8Bzsz64ZO2hRBXaK8dpHacR7QHXyuMYz0EjPj/7Z1/zNVVHcdfb5GnRAJEZQKiUWPR1OYmSs0sI8wlDtL4UXM0aBnOLFyutRpNQzZtQbO0rXQ62krUmRn9opW2LKkBm8tF/mA9IWmtGqbgSED59Mc5F85zvd977/N4r3B53q/tjO/3fM85n88553s/z/me8zkHSecARMRO4CHg0kLWOOAf+d2rMT0innydsg8SETdGxLSImBER0yPikU6Vbcwbhe11dzkC7XWnT4IYjP6t0r6a/633e+6Y7bbdPoQH1j1EPrrmKWCDpHfnuGeATSS/pqnApyUtBZC0UtJX8jLQQT8qSTMkfUPStZL2SHpK0uclfUTSZkmXSfqRpF2SlhX5GpZXJ/tKSWtLX7O8HHqrpJsk/UbSeTm+qbwivyQtk7RV0kWSHpf0L0lL2pAxrdSnhczXtGEDTshLu7slXVnXpjfktnkgG+J22+w1siTNl7RJ0hJJ2yTd3KS8Rn3fUJ9m5JmNrcCCIvpkYI6y2wcwF1hf6HkcsKe4ny7pZkmrJf1WyaVkkaQtkuZJWifpBUnLJZ2T67hT0iWt9DOml7C9Hj72uhFVdWwmR9JUSXcpzf5+vR05DeT+XNIPJJ3fIp1td7eICIceCsAZwDOkL9S7gAnFs+3Ahfn6TKA/X59M+po+kfTFvQOYkZ99G9iYr/uA54HbSGdRfhx4rll59bJJPmbfBdbm+DeTBmvj8v184D/AW5rJa1DvU7LML+e0q4CXgdNbyKjXp6nMsg0b6LAdWAOMyvl25PgxwD1FukdJRrGtNquQNZ60vHcLcB4wu1l5dX3fUJ82368VwF/z9duBm4C9wOXF+6Ii/Vzggnw9CtgMjMr3j+R2GAO8CKwEjgMW5fuFpI/768jvoIPD0RSwvR4W9jo/v6HNdmzW15uBdxZ5os33rNaf7wGub/D8wizn2CLOtrtLwTPWPUakWcWzSIbmE8BWSY02aWwDFksaCVyQ40aTfshTgJp/7BOkHwwRsQ/YDdwfEbtIP7RTWpRXr99eBm6GmAP8L7KPLvAA6Qd5eQt59eXWyrw7p11F8v+9qIWMAfoMRmYFP4nkT7cZmFjU8cQ8o3QtqU2hzTarqO/zJAP244jYFBG/HkR5Vfq0w73A25TcQT4K3AH8ElgoaTSwO7LlzZwPbMzXc0l/vGqzIJcBN+Z2fgF4OJJf9hZgTETcF2l5/DEG1wfG9AS218PDXjegso5N5MwGxkZETZ8tg5R5LnBFRHy1zfS23V3i2MOtgGkfSSOAKRGxHficpDtJO6hvB2aWaSNir6TJwBeB79WKiIh/S9oCfJD0RT2VNJg6mLW4foXsLlRVXoWqZRnTgJGFXgckbQdObSavFRHxsqR+YBxw0iBkDFlmHa9w6PdzGvC3iLilPlE7bSZpMWkAW+MdkZaMo9R1EH1QqU8rImKbpMdI7iBjI6Jf0r2kWaQFwM8KvWvvRs1/73TSDFOtrJ1l0cX1qwzkAHZLM0cZttcD6jcc7HVJZTs26Zt3kT4ihsosYJqk0RHxUkWaA7kOtt1dxA3SW4yk8H+NiD+Rltqm1yeUNAu4OiJWRTrhoeQzwLlKxz49DnytleAW5TVjB2kGtDyGSCTfw9dLH/B0l2W0wz8Z6IeMpJmDaLP1wNlFaLi7fBDlNdSn/epwD7AU6C/0GwEsJy2b1phJ8het8RzwXknHF3Kb+vkZcxRjez2Q4WSvK+vYRM6unKePobGaNNv/zYrnL8ahTee23V3EA+ve47OSJhX3k4EN+XofabPGdNLRTmMlvUnS7Px8gqSTSL5+d5KWcTaSlxYzx1D3lS5JLcqrl62ijAeBl4CP5bLGZXm1mc8qeVVMzGkm57wb2pBR6tNKZlmPehrNXog0C3U88FNJFyttkmnVBwdlARMj4ski7C9EjCium5VX6t1QH0l9km7Pbh7NuI+0BH0/QETsBn4BPBoDTwP5EMlNpEatvdcpnYV9DVA7y7S+D+r7uadPKDCmAttrho29brcdG8oh2dI+0scXpNUJJE1qw3aLNHu8GFggaWHOe52kFaTZ+oeK9Lbd3aQbjtsO3QmkDREB/Bf4PrCWtJR0Qn5+PfAsMI+0tLONNEMwhzTTsZZkpNaTNlLsy+XtIf3QPgzsJ22QmQJ8KT//ZLPy6mQvB35H8ls7Iz+bQZrpXEPyNTwzx1fKq6h/5PzLSa4JZxXPqmS8tdSnlcyyDetkzyNt4rstt0Ut39L8/P3An3PfrMnt3G6bzWtQ1/lZ3jpgco5r1qcDyqrQ51TSf2Rwdxvv2rfq7hcBH6iLW90g3/uAvwA7gRU57hLSu3YrMIl85BhwFekP7x25rpce7t+Yg0OnArbXw8lenwb8njTrfnaLOjbr6zmklcI/kjaSPw1cneve0HaTVkX2k2aqxwMPkwb1nwKWkc7zXkNyQ6nlse3uYlBuPDNMkDQBuCoiVub7Y0izKMsiYsVhVa4FkgKYGsln0QwRSV+IiCEd5WSMeeOwvTYltt29gV1Bhh/XkNwCxhdxsyg2pR2JFEtPXnYaIpJGSLoC+OHh1sUY0xa218a2u8fwwHr48R2SL9YTkv4O/Ap4NiL+cHjVqiZv5lieb5coHftmBs9o4MGI6G+Z0hhzJGB7bcC2u6ewK4gxxhhjjDEdwDPWxhhjjDHGdAAPrI0xxhhjjOkAHlgbY4wxxhjTATywNsYYY4wxpgN4YG2MMcYYY0wH8MDaGGOMMcaYDuCBtTHGGGOMMR3AA2tjjDHGGGM6gAfWxhhjjDHGdID/A4AlFyVWGcoUAAAAAElFTkSuQmCC\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/apollo-as-202-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/apollo-as-202-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/apollo-as-202-earth.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }