{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 47 - Venus Flagship - Decadal Study - 2010 (Concept)" ] }, { "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 Venus Flagship Mission, 2010 Decadal Mission Concept Study (page 159). https://vfm.jpl.nasa.gov/files/Venus+Flagship+Mission+Study+090501-compressed.pdf" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"VENUS\") \n", "planet.h_skip = 220.0E3\n", "planet.h_trap = 10.0E3\n", "planet.loadAtmosphereModel('../atmdata/Venus/venus-gram-avg.dat', 0, 1, 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle1=Vehicle('venus-flagship', 2020.0, 349.0, 0.0, np.pi*2.65**2.0*0.25, 0.0, 0.30, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle1.setInitialState(220.0,0.0,0.0,11.1,0.0,-13.8,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle1.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle1.propogateEntry (120.0*60.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAt0AAAHtCAYAAADFpcLCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzde5gcdZX/8feZyUyY3BMy5AZhwISAEVhhQBSFiCAoiuCCKOsqKoyiuOhPdlFu6wIqF2VRXNFR0UVlRYQVuSqKgqwCDgZEDCGBXCA3JglJJteZzJzfH1WddDpz6Ut1V1f15/U89XR3VXf1+aaT0yffPlVl7o6IiIiIiJRPXdwBiIiIiIiknYpuEREREZEyU9EtIiIiIlJmKrpFRERERMpMRbeIiIiISJmp6BYRERERKbNhcQcQlYkTJ3pLS0tBr+ns7KS5ubk8AVVQGsaRhjFAOsaRhjFAssbx5JNPrnb3ZAQbkWJyNiTrcy1WLYwRamOcGmN6ZI+z2JydmqK7paWFjo6Ogl7T2tpa8GuqURrGkYYxQDrGkYYxQLLGYWZL4o6h0orJ2ZCsz7VYtTBGqI1xaozpkT3OYnO22ktERERERMpMRbeIiIiISJnVdNHd1tYWdwiRSMM40jAGSMc40jAGSM84ZFe18LnWwhihNsapMaZHFOM0d48glPi1trZ6LfQUiUj6mNmT7t4adxyVpJwtIklVbM6u6ZluEREREZFKUNEtIlICd6enp4e+vr64Q0k1d+jtDRYRkSRS0S0iUoLDDz+cxsZGnnrqqbhDSbU//hGGDYNjjok7EhGR4qjoFhEpwbBhweUOenp6Yo5ERESqmYpuEZESNDQ0ACq6RURkcCq6RURKoKK7slJywi0RqUEqukVEStDY2AhAd3d3zJGkm1ncEYiIlEZFt4hICTTTLSIi+VDRLSJSAhXdlaX2EhFJKhXdIiIlUNFdGWovEZGkU9EtIlIC9XSLiEg+VHSLiJRAM92VpfYSEUkqFd0iIiVQ0V0Zai8RkaRT0S0iUgIV3SIikg8V3SIiJVBPt4iI5ENFt4hICTTTXVnq6RaRpKpo0W1m7zWzRWa2xsy+bmbDwvUTzew7ZnaNmV1vZg1Zrxlwm4hI3NJcdFdTzlZPt4gkXcWKbjObDpwKnA58FjgH+Ey4+XbgZne/CHgeuCLrpYNtExGJVVqLbuVsEZFoDavge+0LfNTdtwNPmtnrgLea2aPALHd/PHze3cB8M7sKOHigbe6+qYKxi4j0K9PTnbaimyrN2WovEZGkqthMt7v/IUzeGcuApcBxwKKs5y0L77YOsU1EJHaZme60HUhZbTlb7SUiknRxHkh5BPBtYBqwJmdbFzB1iG276OzspLW1dcfS3t5ehpBFRHZVbHtJe3v7jnwFTIw+ssgpZ4tIzYoiZ1eyvWQHM5sBdLr702bmwJacpwwHeoDBtu2iubmZjo6OcoQrIjKgYovutrY22traADCz1ZEHFqFqytlqLxGROESRsys+021m9cDHgc+Hq5YB47O2GzASWD7ENhGR2KW4pxuonpyt9hIRSbo42ksuBL7q7tvCxw8Cs7K2twDdQMcQ20REYpfWnu4sytkiIhGo9Hm6LwbmAk1mtr+ZfRRYB6wIj4wHOAm40d273b1joG2VjFtEZCBpPWUgKGeLiESpYj3dZnYpcGXO6ufc/WYzOxO43MwWA6OAS7KeM9g2EZFYpbXortacrZ5uEUmqihXd7n4VcNUA25YAHyt0m4hI3NLa011tOVs93SKSdHGeMlBEJPFqoKdbREQioKJbRKQEaW0vqVZqLxGRpFLRLSJSAhXdlaH2EhFJOhXdIiIlUNEtIiL5UNEtIlKCzIGU6umuDLWXiEhSqegWESmBZrorQ+0lIpJ0KrpFREqgoltERPKholtEpAQquitL7SUiklQqukVESqCe7spQe4mIJJ2KbhGREmimW0RE8qGiW0SkBCq6RUQkHyq6RURKoKK7stTTLSJJpaJbRKQE6umuDPV0i0jSqegWESmBZrpFRCQfKrpFREqgoruy1F4iIkmloltEpAT19fWYGe5Ob29v3OGkltpLRCTpVHSLiJQo09et2W4RERmIim4RkRJlWkx0MGX5qb1ERJJKRbeISInU111+ai8RkaRT0S0iUiIV3SIiMhQV3SIiJVJPt4iIDEVFt4hIidTTXTnq6RaRpFLRLSJSIrWXlJ96ukUk6VR0i4iUSEW3iIgMRUW3iEiJ1NNdOWovEZGkUtEtIlIi9XSXn9pLRCTpVHSLiJRI7SUiIjIUFd0iIiVS0V05ai8RkaRS0S0iUiL1dJef2ktEJOlUdIuIlEg93SIiMhQV3SIiJVJ7iYiIDCWSotvM3hnFfkREkihpRXeSc7Z6ukUkqQoqus3sJDN71MwWmtmL4bIYuKOAfZxgZo+bWUvO+qvMzMPl6az1TWZ2o5l92cxuMrMxhcQsIlJu1drTnaacrZ5uEUm6YQU+/0fAV4C5QF+4zoCz8nmxmU0CRgFH5qwfCUwETghXLcna/C3gAXe/zcxOAr4DfKDAuEVEyqaKe7qVs0VEqkShRfeL7n597kozey6fF7v7KjO7q59N5wLPAY+6+9as/U4l+HL4VLjqN8AdZtbi7osLjF1EpCyquL0kdTlb7SUiklSF9nSfZ2an97P+lHx34O592Y/NzID3A18DVplZ9ozIHGC1u28OX7sdWAQcW2DcIiJlU8VFd2pyttpLRCTpCi26jwa+Z2a9WUsfcFOxAXjgKIKfKm8EfmJm7wg3TwPW5LykC5iau5/Ozk5aW1t3LO3t7cWGJCJSkGJ6utvb23fkK4L8Vw7K2SIiEYgiZxfaXnIF0AZ0AL3hOgP+uZg3z+burwKXhrMonwHuBxzYkvPU4cBu32zNzc10dHSUGoaISMGKmelua2ujra0NADNbXZbAUpiz1V4iInGIImcXWnTfD9zl7tuyV5rZd4t58wF8A/h9eH8ZMD5n+2hgeYTvJyJSkio+kDI1OVvtJSKSdIUW3T8FvmBmD2WtM+C9wAURxdQH/CW8/zvgu2bW5O5bzKwBmA48HNF7iYiUrIp7upWzRUSqRKFF92XATOBsgp8RIUjgk8kzgYc/RWZelznafQ7Bl4MDnwMuAXD3lWZ2D/A24B7geOAOd19WYNwiImVTxUV36nK22ktEJKkKLbqvA+53967slfle3czMRgMfDB+ebWbfJPgp8mrgYuBR4IacU0udB1xrZocSfFG0FRiziEhZVevFcUhRzlZ7iYgkXaFF90u5yTs0OZ8Xh6+9iV2PnO8k+PlxoNe8SnBOWBGRqlTFPd3K2SIiVaLQUwZ+MneFmZ0B3BBNOCIiyVPF7SXK2SIiVaLQonuCmV0GYGYjzOx7wPeAP0cemYhIQlRx0Z26nK2ebhFJqoLaS9z9ZDM7ysxuBE4EVgKHAuvLEZyISBJUa093mnK2erpFJOkGLbrN7JgBNi0BFgBfBfYFTgU+G21oIiLJUC093crZIiLVa6iZ7muBg4B1BOdizXUzwWmkpqAELiI1qoraS1Kfs9VeIiJJNVTRfQXQ4e6vDPYkMzs9upBERJKlioru1OZstZeISNINWnS7+3357MTdfx5NOCIiyVMtPd3K2SIi1avQs5eIiEiOaunprgVqLxGRpFLRLSJSoipqL0kttZeISNKp6BYRKZGKbhERGUrJRbeZvd7MfmxmnzOzltJDEhFJlmrp6c6HcraISDyimOn+BHC8u38NeJOZNUWwTxGRxEhYT3eic7Z6ukUkqaIous8HDg/v/xRoiGCfIiKJkbD2kkTmbPV0i0jSFVx0m9k7zOyfwvtHAjPdfRmAu/e5+4aIYxQRqWrVXHQrZ4uIVIeCim4z+xJwO3AagLs/AXzIzM4qQ2wiIolQrT3daczZai8RkaQqdKb7TcDewB+y1t0PXB1ZRCIiCVPFPd2pydlqLxGRpCu06H7Q3dcB2XMNibucsIhIlKq4vUQ5W0SkSgx6Gfh+vGhmFwPTzOz9wBkEP1t+OvLIREQSooqL7tTlbLWXiEhSFTTT7e4/BR4C6oEPAeuAt7v7f5UhNhGRRKjWnu405Wy1l4hI0hU60427PwY8lnlsZvVm9jZ3/22kkYmIJEQVz3QrZ4uIVIlBi24z+yUwaoh9TARWA0rgIlKTsg+kdHcspmlZ5WwRkeo11Ez3X8PnzAf6gHOBuUBH1nP2AyaUJToRkQSoq6ujrq6Ovr4+ent7GTas4B8Ro5L6nK2ebhFJqqG+Ga4DLDz6HTN7jbtfnvskM7urHMGJiCRFY2MjW7dupaenJ86iO7U5Wz3dIpJ0gx5I6e7rM8k7NMnM6rOfY2ZzgCPLEJuISGJUQ1+3craISPUqdDrm58DTZvYo0AUcBJwI7DaTIiJSS6r0AjmpydmZHw+q8FhVEZG8FHrKwAeBtxMchDOL4PRTp7v7V8oQm4hIYlTDTHeuNOXs8eOD27Vr441DRKRYxZwycDlwaeaxmdWZ2eHu/mSkkYmIJEgVn6s7FTl7xAhobIStW2Hz5uCxiEiSFFR0m9nN/axuBnqBUyOJSEQkgapxpjtNOdsMJk6E5cuhsxP23TfuiEREClNQewkwE1gKLMlahgOJmjEREYlalfZ0pypn779/cPv88/HGISJSjELbS85x9/nZK8zsQODsyCISEUmgapzpJmU5e/ZsePRReOYZOOGEuKMRESlMoQdSzu9n9TbgnGjCERFJpmosutOWs9/whuD2wQfjjUNEpBiF9nQvArKvB1YPTAJ+VMA+TgCuAs5098XhuibgWoJTWo0HLnL3DUNtExGpFtu2bQNgjz32iDmSndKWs08+Gerq4Le/hZUrYfLkKPYqIlIZhbaXPAjcys4k3gescve8OuzMbBIwit0vzPAt4AF3v83MTgK+A3wgj20iIlVhw4agrhwzZkzMkewiVTl7r73gPe+B//1f+NrX4LrrotiriEhlmLsP/azMk82OcfdHctbNAnrdfWGe+6gjOHJ+P3dfbGZTgUXAeHffbGbDgPXAbKB7oG2ZGZeM1tZW7+joyHssIiJRGjt2LBs2bGDdunWMHTu2oNea2ZPu3hp1TGnM2X/+c9BmUlcHDz8MRx9d8C5EREpSbM4u9OwlJ/WzbhVwe747cPe+nFVzgNXuvjncvp0gaR87xDYRkarQ1dXFhg0bGD58OKNHj447nGypy9lHHAGf/Sz09sK73w2PPRbVnkVEyiuvotvMzjezdcBFZtabvQBrCJJ4saaF+8jWBUwdYpuISFVYsGABADNmzKCurtC5jOilPWdfcw2ceiq8+ioceyx8+csQttSLiFStvL4d3P2bwDHADcD+Ocve7t7fbEq+HNiSs2440DPEtl10dnbS2tq6Y2lvby8hJBGR/D3xxBMAzJ49O+/XtLe378hXwMQo40l7zh42DH72Mzj/fOjuhksugQMPhG98AzZuLGwwIiL5iCJn530gpbv/1cw+7+67JU8zm+7uS4sJAFhGcIR7ttHAcoIEPtC2XTQ3N6OebhGJw9133w3ACQWcPLqtrY22tjYAzGx11DGlPWc3NMCNNwYz3uefD889BxdcAF/4ArzrXfC+9wXn8q6u41pFJKmiyNmDHkhpZhcC97n738PHlxCccmqXpwFvdfc5eb+pmbPzoJzJwEKg2d23mFkDsJFgRqZ3oG3uvix7nzqQUkTiMG/ePGbPnk1DQwNLly5l0qRJBe8jqgMpazVn9/bCL38J118fXDwno74eWlvhrW+FY46Bww6DIj4eEZFdFJuzh5rp3g8Yl/X4MGAksJKdp6AyIO+zpZqZZb0Od19pZvcAbwPuAY4H7sgk6MG2iYjEadu2bZxzzjm4Ox/5yEeKKrgjVpM5u74eTjstWJYsgdtvhzvuCM508vjjwXL11cFzp06F178+WGbPhgMOCJZRo8oZoYhI4acMPAz4m7t356x/bWZmZYjXjwY+SHAe1yuAb7p7p5mNJ7iYwmKCL4MvuPvG8DUDbsummW4RqaRNmzZx1lln8ctf/pIpU6bwt7/9jQkTJhS1rzKeMrCmc3ZXF/zhD/C73wWF91NPBev6M3UqzJoFM2dCSwvssw9Mnx7cTpsGjY1lDVVEEqTYnF1o0f1hd//vnHWzgAPc/e5C3zxKKrpFpFIeeeQRzjvvPP7+978zbtw4Hn74YQ455JCi91fGols5O0tfH7zwAsydGxTgzz0Hzz8PCxYEB2QOxCy4+uX06bD33sH9SZOC2+xl0iQV5yK1oFztJZmdv4Xgp8UTwssKZ5sAfAOINYGLiJSTu/PII49w3XXXce+99wJw4IEHcuedd3LQQQfFHN2ulLP7V1cXzGTPnBkcaJnR2wtLl+4swJcu3bm89BIsXw4rVgTL448P/h4TJgTF9557BsuECTvv5z7O3N9jj/KOW0SqQ75nL1kH3ELQL5h7/a9u4HtRBiUiUg3cnWeffZbbb7+dW2+9lYULg4s4NjU18W//9m9cdNFFNDU1xRxlv5SzC1BfD/vtFywnnrj79p6eoPBeujS4XbkSVq0KbjPLqlXBsnZtsBSiqQnGjoVx44LbQu+PGRP8h0JEqlteRbe7P2Nmc4Az3F3JWkRSyd158cUX+dOf/sRDDz3Er3/9a5Yt23kM4OTJk/nEJz7BJz/5SZqbm2OMdHDK2dFqaIB99w2WwfT1wZo1QfG9Zk2wrF27835/j9esgS1bgmXlyuLiMwsOBM0so0eX/niPPYL9ikh0CjlP93oGmB0xsze7+6P9bRMRqUa9vb0sXLiQZ555hmeeeYa5c+fy2GOP0dnZucvzJk2axDve8Q7OOussjjvuOOrrc8/AV52Usyuvrg6am4MlX+6waROsXx8s69btejvQ/ex1XV07lyjHMnJksIwYsfvS3/p812XWNzaqsJfaMmjRbWa/BIY6kZIBs4G9ogpKRCQKfX19rFixghdeeIGFCxeycOFCXnjhBRYsWMC8efPYunXrbq9pbm7mjW98I0cffTQnnngihxxyCJaQykA5O3myZ6mnTStuH729QcG9cePuS3/r81m3dWv0hXyuurqgtWaPPXbeZt8v162KfYnLUDPd8wl6A5ex8xyvuQx4V5RBiYjkY8OGDbz88su89NJL/d4uWrSILVtyr0q+0z777MPBBx/MIYccwiGHHMIb3vAG9ttvv8QU2f1Qzq5B9fVBf/e4cUM/N1/btwcz8Js3775EsX7TpuCMMZs2BUslmfVfjBdTwBe6NDSo4K9lQxXd1wBd7r5tsCeZWc0dBS8i0XN3urq6WLVqFatWreKVV17Zcb+/x115TMNNnDiRGTNmMGPGDF7zmtfsuD3ooIMYF2WVUh2UsyUSw4btPGCzXLZvD2bUt24NetorddvTs7OPvtIyBf9Ay/DhxRXz+S6a5Y/XoEW3u+d7bfnT0NHwIpKlu7ubtWvXsmbNGtauXbvL/f7WrV27ls7Ozn5bPgayxx57sM8++7DPPvuw995773bb0tLC2HJWDVVGOVuSZNiwna01ldTbG02xv2ULbNu2c1+DLZnnxVnwZ5S7sB9s/8OH1/aZdvI+kLI/ZjYd+BRwPkrgIqnR29tLV1cX69evZ926daxfv77fZaBt69atY1ORvxmPHDmSvfbai0mTJjFp0qRB748fPz7JrSAVp5wtErTjZA4SrbTe3vwL9UKXfPbb3R08b9u24CDcODQ29l+QZz6TUaOKvx01qrqL+oKLbjOrB04BzgVOAFYR9BCKSEzcna1bt7Jx40a6urp23GbfH2zdxo0bdyma82nbGEp9fT177rknEyZMYMKECTvuD7Zuzz33ZFSlp71STjlbpHrU1+88g0sc+vqKL/qj+M/Ctm1B4d/dDRs2RD8+s+C89ZljHPpbxo8PTv+5//5wwAGVvThV3kW3me0HnAN8BNgTWAK8G/gV8MayRCeSMpnieNOmTWzatInNmzfvuJ+95K4fqFDOLqJ7e3sjjXXMmDGMHTu232XcuHFDbhs9erRmoWOknC0iuTJnjInrml59fUHBnVuMb9kSHGC7cWNwYG3ubX/r+rvduHHnqTSXLBk6nvHj4aGH4B/+ofxjhzyKbjM7g2CG5DhgIXA98N/Axe7+QPi0P5YtQpEK6+npGbQALmX95s2b6evrK0vcw4cPZ9SoUYwePZrRo0fvuJ97O9C27CJ79OjRiTkftexKOVtEqlVd3c52kozsA2rzmVnPnP2mv4J8/XpYsQKWLQsuPDWUV1+Fl16qoqIb2EJw2eA/Aue4+/MAZjbQ6ahEyipTFGcXt8Xezy2QN2/eTE9PT1njHz58OCNHjmTkyJGMGDFix/3cJXtbPsVzQ0NDWeOWxFDOFhEguPjS9u07WzriWDI95AMV0xH/SNuv0aNhwoRgmToVWlpgv/3gmGPgiCPK//4ZQxbd7n4PcE94AM65ZtYC3AXs+IY3M3N3JXQBdp0pLqTozfd+uYviurq6vArhfLdlrx8xYgTDhpV0/LLIoJSzRYrnHhSB27cHZxrJZ8n3ueXY51DP7e6O+090aNmz3/mcEaWpaeiDKjP399wzaCFpbIx7lIFCLgO/FLjMzIYBpwJ7m9mVwPeBtwPt5QlRolZIUVzM/e3bt5c1/tyiOLeozb0/1Pbc+42NjepFlsRTzpZK6+vbveDLvi30fj7PLWTJt5BNm2HDgqIzrqWhYfCCupbmoQoeqrtvB34O/NzMZgCfJegfVAIvk+3bt9PV1cWGDRsKWrq6uvotjMtdFNfX15dc+A72OhXFIvlTzk62zCnmMj/JZ+7nPs7dVsnCN3O/Em0ClWAWFIpDLcOG5fe8OJ/b2Fjdp9CrNSX9/8LdFwIXmNmTEcWTWj09PbtcDCT7IiHZ91999dXdiufNmzdHGkt2UVyOmeKGhgYVxSJVSDm7ND09u1/WPHN/sHXZF1HJp2DOflzmOZLIZRd7+d4v5Lm598tRyKpIlXKJZFLf3W+JYj9JtGXLFhYvXszy5cv7XVasWMGaNWvYUMIJKc2MMWPG5LWMHTt2x/2Beosbq6W5SURiUas5u7cXXnklOLvBqlWwbl1wtoN163a/n3m8fv3O4jmu1oNMT2umr7W/+9mPc3/ar1QRPGyYLjEuMpga6qQpzebNm/nLX/5CR0cH8+fPZ8GCBTz//PO89NJLeb2+rq5ulwuCZF8MJPv++PHjdymcM8WzZo5FRIb28sswbx7Mnx8sL7wQFNkrVkBnZ9B3XKy6uuAArREjdr8dbF1T0+BF8mD3GxpUyIqkhYruAbg7Tz/9NHfeeSf33nsvTz/9dL8XH6mvr6elpYW9996bqVOn7rZMmTKF5uZmxowZQ51+sxIRidT27XDvvXDHHfDww7B06eDPb26GyZNh0qTg9GFjx+56tbrsx2PHBsuoUUEB3dioAlhEiqeiux8PPvggl19+OY899tiOdXV1dRxyyCEceeSRvO51r2PmzJnMnDmTlpYWnR9ZRCQGv/41fPzjsHjxznXjxsEhh8CsWcEyYwZMm7az0Fa6FpG4qOjO4u5ceOGFXH/99QBMmDCBM844g9NOO42jjz6aUaNGxRyhiIhAUHC/851Bn/YBB8DHPgYnnQSve50OhBOR6qSiO8uVV17J9ddfT0NDA1/84he54IILGDlyZNxhiYhIjgsvDAruz34Wrr22ts71KyLJpDQVWrBgAVdddRVmxi9+8Qve+c53xh2SiIj0o7MTnnkmOEjx6qtVcItIMuhHuNB//ud/0tPTw4c//GEV3CIiVayzM7jdZ5/qubyziMhQVHQTnA7wlluC09ZeeOGFMUcjIiKD2bQpuB0xIt44REQKoaIbuO+++9i0aRNHHnkks2fPjjscEREZRKbo1rHtIpIkKrqBn/3sZwCceeaZMUciIiJD6eoKbkePjjcOEZFC1HzR3dPTw3333QfA6aefHnM0IiIylEzRrZluEUmSmi+6n3rqKTZt2sTMmTOZPn163OGIiMgQNmwIbseOjTcOEZFC1HzR/eijjwLwlre8JeZIREQkH+vXB7cqukUkSWq+6H722WcBOOyww2KORERE8pGZ6R4zJt44REQKUVVFt5ldZWYeLk+H65rM7EYz+7KZ3WRmkabZ+fPnAzBr1qwodysiknpx5GyAzZuDW50yUESSpGqu42VmI4GJwAnhqiXh7beAB9z9NjM7CfgO8IGo3nfBggUAHHDAAVHtUkQk9eLK2RBc/h2gvj7KvYqIlFfVFN3AucBzwKPuvhXAzKYCZwGfCp/zG+AOM2tx98WlvmF3dzerVq2ivr6eadOmlbo7EZFaUvGcnaGiW0SSqCraS8zMgPcDXwNWmVlmVmQOsNrdNwO4+3ZgEXBsFO+7cuVKACZNmkS9sreISF7iytkZKrpFJImqYqbb3R04yszGA58DfmJm64BpwJqcp3cBU3P30dnZSWtr647HbW1ttLW1Dfq+y5cvB2Dq1N12JyJSVu3t7bS3t2ceTowzlkLFlbMzVHSLSKVFkbOroujOcPdXgUvDWZTPAA8CW3KeNhzoyX1tc3MzHR0dBb2fim4RiUt2kWlmq2MOpyiVztkZKrpFpNKiyNlV0V7Sj28A04FlwPicbaOB5VG8yYoVKwCYMmVKFLsTEalVFcnZGX19wa2KbhFJkmotuvuAvwC/A6aaWROAmTUQJPaHo3iTNWuCX0EnTkzUL7siItWmIjl7x5uFRbdZlHsVESmvqii6zWyqmZ1lZnXhz5SfAy5x95XAPcDbwqceD9zh7suieN9XX30VgPHjcydmRERkIHHl7N3jKMdeRUTKo1p6uscDVwMXA48CN2SdXuo84FozOxSYDOR3pE0e1q5dG7y5im4RkULEkrNFRJKsKopud3+W4CfI/ra9SnA+2MhlZronTJhQjt2LiKRSXDlbRCTJqqK9JC5qLxERSR73uCMQESmcim5UdIuIJJF6ukUkSWq66O7q6gJg9OjRMUciIiIiImlW00X3li3BNRxGjBgRcyQiIiIikmY1XXRv3rwZgKamppgjERGRfKmnW0SSqGaLbnffMdOtoltEJHnU0y0iSVKzRXd3dzd9fX00NDTQ0NAQdzgiIiIikmI1W3RrlltEREREKqVmi+5MP7cOohQRERGRcqvZov6MPbcAACAASURBVFsz3SIiyaQDKUUkiVR0q+gWEUkkHUgpIklSs0V3d3c3AI2NjTFHIiIiIiJpV7NF9/bt2wEYNmxYzJGIiIiISNrVbNHd09MDoNMFiogkjHq6RSSJVHSr6BYRSST1dItIkqjoVtEtIiIiImVWs0W3erpFREREpFJqtujWTLeIiIiIVIqKbhXdIiIiIlJmNVt0q71ERERERCqlZotuzXSLiIiISKWo6FbRLSKSSDpft4gkiYpuFd0iIomi83OLSBLVbNGtnm4RERERqRQV3Sq6RURERKTMarbo7u3tBaC+vj7mSEREpBjq6RaRJFHRraJbRCRR1NMtIkmkoltFt4iIiIiUWc02NM+YMYNTTjmF2bNnxx2KiIgU4LDDoKsL9t477khERPJnnpKmuNbWVu/o6Ig7DBGRgpnZk+7eGncclaScLSJJVWzOrtn2EhERERGRSlHRLSIiIiJSZono6TazJuBaoAsYD1zk7hvijUpERPqjnC0isrukzHR/C3jU3S8G7gK+E3M8IiIyMOVsEZEcVV90m9lU4Czg7nDVb4BTzKyl1H23t7eXuouqkIZxpGEMkI5xpGEMkJ5xJE05czbUxudaC2OE2hinxpgeUYyz6otuYA6w2t03A7j7dmARcGypO07LX5Q0jCMNY4B0jCMNY4D0jCOB5lCmnA218bnWwhihNsapMaZHFONMQk/3NGBNzrouYGr2is7OTlpbd569pa2tjba2tvJHJyJShPb29uwkPjHOWCKmnC0iqZOTs4rK2VV/nm4zuxA4w93fkLXuL8Ct7v7VrHWdwJICdz8RWB1JoPFKwzjSMAZIxzjSMAZI1jj2dffmuIOIQplzNiTrcy1WLYwRamOcGmN6ZI+zqJydhJnuZQRHv2cbDSzPXpGWLywRkYRTzhYR6UcSZronAwuBZnffYmYNwEZgf3dfFm90IiKSTTlbRKR/VX8gpbuvBO4B3hauOh64Q8lbRKT6KGeLiPSv6me6AcxsPMGFFhYDk4EvuPvGEvaXigs3mNl7ga8BY4AfA58LzxSQOGbWCPwZuMDdfx9zOEUzs6OANxPM9D3q7onqczOzg4DzCeKfCbS7+1PxRpUfMzsBuAo4090Xh+tS8W89aaLO2eE+U/lZDpTHzWwi8CVgHdBAMN6e+CItXW6eT/Fnusv3ALCJlIxzoO+INHyWhX6HFDVmd6+5BfhB+IcKcBLwP3HHVMQYpgO3AIcDHyL4R31h3HGVMJ5LgPXAnLhjKWEMHwGuijuOEsfQAUwL708H5sUdU55xTwJOAxxoyVqf+H/rWtL7WQ6Wx4HfAW8I738C+Erc8UYw3l3yfEo/092+B9I0zoG+I5I+xmK+Q4oZcyJmuqMUXrhhETDe3Teb2TCCJDDbw//ZJIGZvQX4k4cz22Z2LcEYTo43ssKZ2ZuAA4HLgbM9gTPdZnYscBlwgif4H5WZbQIOd/fnzKwZeNrdpw71umpgZnVAL7Cfuy9Oy791SU/ezjVQHgeuBO7M/Nszs2nAfGCSu2+KK95S5OZ54HlS9pn29z2Qtr+7/X1HAK2kYIyFfIcA3QNtG2zMVd/TXQZzKOOFGyrF3f/gu7aSLAOWxhVPscxsJHC6u98cdywluh6YB9xoZg+Y2dFxB1SknwPfM7PRwAeBT8ccT97cvS9n1RxS8G9dgJR+loPk8eMIxpd5XqYfvpUEGiDPzyF9n2l/3wNzSNc4+/uOmEMKxljgd8hg2wZUi0V3XhduSKAjgG/HHUQRLgKujjuIUpjZLOD1wHfd/Xzgt8ADZjYp3siK8klgK0Hf5UZ3vyPmeEqR1n/rtahWPstMHk/bePvL86ka40DfA8A+pGic9P8dkarPMstg4ypqzLVYdDuwJWfdcCCxB6iY2Qyg092fjjuWQpjZiUCHu78Sdywlmg2sdfe/ho+/CRjwnvhCKloT8D/ArcANZnZ8zPGUInX/1mtY6j/LnDyemvEOkudTM8bQQN8DaRtnf98RaRtjxmDjKmrMSbg4TtTyunBDUphZPfBx4PNxx1KEC4HXm1nm8XjgLjP7krtfG19YBRtG1r8lD85NvADYM76QivZj4P3uvs7MHPgfM2tJaB9pqv6t17hUf5b95PFlwAFZ2w0YSTLH22+eJ2jFSNNnOtD3QD3pGudu3xHAZ0jXGDMGyzs+yLYB1eJM9++AqeGpXggv3DAdeDjWqIp3IfBVd98WdyBFOAv4h6xlOXAOyWuT+SswNqedZDvw95jiKUp4irJD3X1duOoqYANwUHxRlSRt/9ZrWdo/y9w8/iAwK2t7C8GBWx0VjisKA+X5H5Cuz3Sg74HlpGScg3xHLCElY8wxWN4pKifVXNHtKbpwg5ldDMwFmsxsfzP7aPgTZSK4e6e7v5xZCI4a7vSEndvT3Z8D7gdOBzCzcUAjcG+ccRVhLbA1PFNCxhqCswxUPds5lWaQrn/rtS7Nn2V/eZzg3NwrzOx14dNOAm509+644izWIHl+KSn6TAf5HvgJ6RnnQN8RfyUFYyzkO6TYnFRzpwyE8ly4odLM7FKC00ple87dkzoriZktJrmnDJwIfJ1gJmo6wcE0iZrpBjCzQwkOlHmS4Lylj7h71c9WZB1J/y3gCuCb7t6Zhn/rEkjjZzlYHjezfQlOr7cYGAVc4gm9+Fm27Dyfts90oO+BNI1zoO+IpI+xmO+QYsZck0W3iIiIiEgl1Vx7iYiIiIhIpanoFhEREREpMxXdIiIiIiJlpqJbRERERKTMVHSLiIiIiJSZim4RERERkTJT0S0yADN7m5n90czmDPG815nZz8zssgjf+7qo9iUiEgczO8TMbo8yNw7wPveY2fvL+R7lYmbHm9mfhvqekXQYFncAIkMxs4OALxFcDesI4BDgn939x2V+65eAQ/N43lJgPFAfxZuGl5PtiWJfIiIAZnYE8BXgGODXwCKCC3o0ADe5+6/K8LYLgbFElBsH8S3gb+XYcZiP3wt8m+DqvH8KNw0H/hE40t0Xl/AW+X7PSAqo6JYk+AXwPnd/2szqgB9U4k3d/XkzW5PH8zaY2YoI3/pY4JEI9yciNc7d/2xmPwUOdPd3Zdab2cnA/5rZhe7+rYjfc3PEuREAM3s38Eym2HX3+6J+jwx37wFuM7NPA79390uz4rirmH2a2V7AW939Nnefb2arIwpXqpzaS6SqmdkE4ABgE4C79xFcGrlSl1Lti/h5+Xgr8LsI9yciArDbZeTd/V7gi8B/mtnkMrxnpLnazPYmmNmutN3+7IDfAMsK2YmZDQduAZqyVkf5/SFVTEW3VDV3XwvMBx4ws6PCdUuAJ8xsnJldaWZPmNkpZrbSzOaZ2WGZ15vZB8zsEjP7vZl9OWt9q5l90cxuM7M7zWxsuH4PM/svM7vczL5P8NNoQczs38KY3mdmR5jZL8L9XRnGeJuZNYfr15vZv+fsotHdt4X7+pyZ/T8z+z8zu7rgP0ARkaF9n6DN5AwAM3uNmV1hZt8zs9+GhS5mNjzMZZeY2R/M7MzMDszsk2Z2mZk9ZmbnD/RG/eVkMzs9zJlnm9kCM7vazKaa2Q/CfPoHMzsu3MVZwN7ARWZ2anjszW/M7MNZ73GumX3JzG4xs/82s5FmNtbM/t3MOszsODObb2YvmllLMX9g4Z/JnHAmHDN7r5l93cy+ZmYPmNmocAxfN7O7w++bFeGf8WuBM8zs/2Xtcnz4XdRlZucO8J51ZnaNmbWF7+Fm9r+Zz0cSwN21aKnqBZgNLCGYDbgZ2CtcX0+QwNYS9NaNAx4EniX4D+UbgUvD544nmC0/GRgD/DRr//8HXBfe/0+gLbw/CeglSKxDxfhDgtmivQh6/0aE64cDvwLuAqYSfFn0AP8OjAAODx9nnp/pV8+M+1vh/dHA5+P+LLRo0ZLcBTgbeHmAbauAb4R59U6gPlx/K3B7eP964L3h/Y8CneH9DwAfDO8fEObqg8PHPwS+GN4fKCdPANYBNwBHAscDXwUuD5/7r8DdWbE60BLebwLmAWeHj08Ebs167j3AdwjaaU8Lvy/eG47zAeBLef7Z/T78rrgauI6gv/v4rO2rgenh/WcJvpOGA5cALwNvAD4W5v3fZ+INn78Y+Fq47QPA0gFiOAe4J7w/Kfxzfnvcf6+05L+op1uqnrs/a2YHA1cBnwTebWbHe9Dj3QlscPc7AMzsi8CjwP4EXzB9ZvaZcFe/IkjuJwN7Zq2fF752FHBe+FrcfZWZvVRAqBOBG4GPufvmcB/bzGwlsMjdl4fvsxJ42IN+x7kEXwYTCQ7IPJngCwJgM/AhM3seuIngPxwiIuXgBIXokcC+wKfNDGANMNaCBx8nmDAA+G+C9goIcm1HTq6dAjyT8x5n009Odve1ZrYeuMvdnwAws6XAJjObQjAzPKrfoN23mNkrWavOD/eb8T3gZ+H6dQTfF3eG79ERxpmv33nY021mtxFM9GSc4e5LzexogmJ7VJj/lwMvuPvjwOPha/vb993hd8KfB4npMGBLOO5VZraWXdtUpMqp6JaqZmb1wD4eHDDzL2HLx31AO8HMQa5nw9txwHTgR+5+a7juhnCfFxEUwTfkvNfhBMmyq8hwDwDeDLwGeDprfW5PY++ODe59YQLOtHrt6UFLDe6+yMzOBr4JfBb4J+AVREQiZGbjCH6le5Ygb67pJz9OIpiJbQRw916CiQLC1/yHu/8xfLzLa7P0m5NDzq658mXgImAB8Biw3yBDyH7dTOChrMcvELTO7MXuuXg7xZ9Z5S8E/znJWGdm1wC3E8SeqaxzxzWU7Qxcm90P3GRmexB8FlsIZs0lIdTTLdVuR58hgLs/DVwMHDjA8xsJEtxCYAVwevZGM3tDuP5kM2vKWb8hfHhQkbH+Efg68JMwKRbEzKaSdVCOmU0jOHPLAQQHVt46wEtFRErxYYJf1m4nyI9HhzPMwI78uBroBk7KWn+wmY0mJ9ea2bBwEiPXQDm5P98GFntwathtBYxlKcHM+I63ANYTtM9Ext37womR2eF/Wu4DrnD3jijfJ+c97wbuJvhF9gPAse6+vlzvJ9FT0S1J8OmwIM2YRtCLlzHBzBrD+3MI+vnWAf8DnGZm3zSzOeEsxCaC5DgSuMfMTgwPWjmQoFCfC1xpZo1mNpGgl3ovM6s3s8+b2QcHiNHC5XKCL6br+9mW+/zcx+8i6D/MmAmc6u4bgE/38xoRkULsNoNqZscA/wF83N07CWaVVwH3mdm7w4Ml3xLObP+U4Cwnp5rZ24APu3sXQa79FzP7DzM7lmDyIdOal53/BsrJGdmzzq8Hms1sBHA00GRm+4fbeggOPMxMvmS/x7eBfzSz8eHjNwDfdfftBDVPv7nYzI4K4xo+wJ9dHTk1U/gfk5MJft0cT9C2OAvYBxiRdZBm7mx6dz/x7xpUPz0o4WcxBvgtwWx+If8ZkSqg9hJJgn2AZ83sXoKf3gz4RNb27cAlZraOYFb40wDu/qCZ/QvBT5SnEhyI+DcAM3sP8F8EXyI3A//q7m5m7yM4ndOzwG0EPxP+A8HFJN5H0Layy0V5zOxIggOEphP83PgL4D/CeH4FHAW0mNlrCGZgJgPvM7Pn2DnrczYw0t3bc8b+AwsuarGd4MAlEZGChXnkLILC8LtAJ8HBeGOBd7j7nwDcvTvMj98lyHV3EfRyA1xAMGFxC8Eve5kzhnyfIPedR9AGd567vxIWlUcBrzGzHw+Uk83sdII+5nPM7Dl3X0ZwUOdXCQruW4D3AAcDLwI/IpiV/5AFZ7WaDbzTzO5z9zvNbCZwh5n9H8GvpZeEM/JnAJPD8T1LcKGgKWY2O4zz/QR96r/I+nNrIMj9hwLTwqLcCXqp3xvu8yngYeDPBC0z9xPk9HsIDt481MxOdffMfm8juFDRxrBInwycbmaLCGawCV+fe02KpQQHip4K7AHUmdlvCA6mrNRpdKUEps9JksyCS+f+0N1bKvR+/+ruukS7iEjKhK0uw929Ki9OZmb/SnAmmcXh4z0IfqW4JnMskFQ3tZdI0vXXulGeNwp+hv1LJd5LREQqJzxQ9LAqLrj3AC4FWrJWTwQ2qeBODrWXSGKF/XQfJPi58Ax3v73Mb/lC+LOniIikyzB3vynuIAbi7lstuBT9N8KD7JcTtMJcOvgrpZqovUREREREpMzUXiIiIiIiUmYqukVEREREykxFt4iIiIhImaXmQMqJEyd6S0tLwa/r7Oykubk5+oDKSDFXThLjTmLMkMy4o4r5ySefXO3uyRp8iWopZ/dH46guaRkHpGcs1TyOYnN2aorulpYWOjoKv/pqa2trUa+Lk2KunCTGncSYIZlxRxWzmS2JIJxEqaWc3R+No7qkZRyQnrFU8ziKzdlqLxERERERKTMV3SIiIiIiZVbzRXdbW1vcIRRMMVdOEuNOYsyQzLiTGHPSpeXPXOOoLmkZB6RnLGkZR7bUXByntbXVq7X3R0RkMGb2pLu3xh1HJSlni0hSFZuza36mW0RERESk3FR0i4iIiIiUWe0W3T/8IRx2GHz963FHIiIiQ+h6qouOwzp47iPPxR2KiEhRarfoXrUK5s6F5cvjjkRERIbQt6mPjXM3svn5zXGHIiJSlNotujNSciCpiEhNUMoWkYSq3aLbLO4IREQkX0rZIpJwtVt0i4iIiIhUiIputZeIiCSHUraIJFTtFt1qLxERSQ6lbBFJuNotujM00y0ikhhpuYqyiNSe2i26NdMtIpIcStkiknC1W3RnaNZERERERMqsdotuzXSLiCSP5klEJKFqt+gWEZHEME2UiEjCqehWe4mISHIoZYtIQtVu0a1ZExGR5FDKFpGEq2jRbWbvNbNFZrbGzL5uZsPC9RPN7Dtmdo2ZXW9mDVmvGXBbJDTTLSLSr6rM2SIiCVWxotvMpgOnAqcDnwXOAT4Tbr4duNndLwKeB67Ieulg20oJKJLdiIikUdXl7AzNk4hIQg2r4HvtC3zU3bcDT5rZ64C3mtmjwCx3fzx83t3AfDO7Cjh4oG3uvqmCsYuI1JrqytmaJxGRhKvYTLe7/yFM3hnLgKXAccCirOctC++2DrEtqsAi25WISFpUa87WFSlFJKniPJDyCODbwDRgTc62LmDqENt20dnZSWtr646lvb198HdXe4mIxKi9vX1HvgImxh1PHmLO2UVGLSISgShydiXbS3YwsxlAp7s/bWYObMl5ynCgh6B7b6Btu2hubqajo6PwYDRrIiIxaGtro62tDQAzWx1zOIOqqpwtIhKDKHJ2xWe6zawe+Djw+XDVMmB81nYDRgLLh9hWaiAl70JEJO2qJmdnaJ5ERBIqjvaSC4Gvuvu28PGDwKys7S1AN9AxxLZoaKZbRGQwVZGzdUVKEUm6Sp+n+2JgLtBkZvub2UeBdcCK8Mh4gJOAG9292907BtoWQTAl70JEJM2qKmdnaJ5ERBKqYj3dZnYpcGXO6ufc/WYzOxO43MwWA6OAS7KeM9g2EREpg6rL2ZonEZGEq1jR7e5XAVcNsG0J8LFCt0UUWNl2LSKSVFWbs0VEEirOUwbGS+0lIiLJo3kSEUmo2i26MzTTLSJS/TRPIiIJV7tFt2a6RUQSR1ekFJGkqt2iW0REkkPzJCKScCq6NWsiIiIiImVWu0W32ktERJJH8yQiklC1W3RnaKZbRKTq6YqUIpJ0tVt0K4GLiCSP5klEJKFqt+jO0Ey3iEj10zyJiCRc7RbdmukWERERkQqp3aJbRESSRz9OikhCqehWe4mISPXTj5MiknC1W3SrvUREJHF0RUoRSaraLbozlMBFREREpMxqt+jWTLeISHIoZYtIwtVu0S0iIsmjHydFJKFUdKu9RESk6umKlCKSdLVbdCuBi4iIiEiF1G7RnaGZbhGR5FDKFpGEqt2iWzPdIiLJoZQtIglXu0V3hma6RUSSQylbRBKqdotuzXSLiCSHUraIJFztFt0iIiIiIhWiolvtJSIiiaHLwItIUtVu0a32EhGR5FDKFpGEq92iO0OzJiIiyaGULSIJVbtFt2a6RUQSQ1ekFJGkq92iW0RERESkQlR0q71ERCQ5lLJFJKFqt+jWT5UiIsmhlC0iCVe7RXeGZrpFRJJDKVtEEqp2i27NdIuIJIdStogkXCRFt5m9s4DnnmBmj5tZS876q8zMw+XprPVNZnajmX3ZzG4yszFRxLyDZrpFpMYkOmeLiCTUsEKebGYnAZcCk9lZsNcBk4CmPF4/CRgFHJmzfiQwETghXLUka/O3gAfc/bbw/b8DfKCQuAcIpuRdiIhUs1Tl7JCuSCkiSVVQ0Q38CPgKMBfoC9cZcFY+L3b3VWZ2Vz+bzgWeAx51962ZlWY2Ndz3p8JVvwHuMLMWd19cYOwiIrUmPTlb8yQiknCFFt0vuvv1uSvN7Ll8d+DufdkXObDgwfuBI4D/MLNPuPv/hJvnAKvdfXP42u1mtgg4FlhcYOwDBRTJbkREqlAKc3YkexERqbhCe7rPM7PT+1l/SrEBeOAogp8qbwR+YmbvCDdPA9bkvKQLmJq7n87OTlpbW3cs7e3tg7+x2ktEJEbt7e078hVB/iuH1ORsXZFSROIURc4udKb7aOBKM7sta50RzD0MUeUOzt1fBS4NZ1E+A9wf7ndLzlOHAz25r29ubqajo6OYNy78NSIiJWpra6OtrQ0AM1tdprdJX84WEYlBFDm70KL7CqAN6AB6w3UG/HMxbz6AbwC/D+8vA8bnbB8NLC/5XTRrIiLpl56cnaF5EhFJqEKL7vuBu9x9W/ZKM/tudCHRB/wlvP874Ltm1uTuW8ysAZgOPBzh+4mIpFV6crbmSUQk4Qrt6f4p8AUzOyZrORb4Qr47sJ2NeRY+nmpmZ5lZXbjtc8AlAO6+ErgHeFv4muOBO9x9WYFxD0ztJSKSXinM2ZHtSUSkogqd6b4MmAmczc7UZwTngL1gqBeb2Wjgg+HDs83smwQ/RV4NXAw8CtyQc2qp84BrzezQ8H3aCox5oGAi2Y2ISBVLUc6OZC8iIrEptOi+Drjf3buyV+Z7dbPwdTeFS0Ynwc+PA73mVYJzwpaHZrpFJL3Sl7NFRBKq0PaSl3KTd2hyFMFUlGa6RST90pOzQ7oipYgkVaFF9ydzV5jZGcAN0YQTAyVwEUmv9ORszZOISMIVWnRPMLPLAMxshJl9D/ge8OfIIyu3hobgtme308eKiKRFenJ2huZJRCShCurpdveTzewoM7sROBFYCRwKrC9HcGXV1BTcbsm9joOISDqkKWfripQiknSDFt1mdswAm5YAC4CvAvsCpwKfjTa0Mttjj+BWRbeIpESqc7aISMINNdN9LXAQsI7gAgi5bibotJtC0hJ4ZqZ769Z44xARiU56c3aG2ktEJKGGKrqvADrc/ZXBnmRmp0cXUoWovURE0ie9OVvdJSKScIMW3e5+Xz47cfefRxNOBam9RERSJtU5O0Mz3SKSUIWevSQ9RowIbjdtijcOEREZkjUGU9193f11zYiIVL/aLbonTgxuV6/WubpFRKpcfVM9AH1bVXSLSDLVbtE9alTQYrJli2a7RUSqXN0ewddV3xYV3SKSTCUX3Wb2ejP7sZl9zsxaSg+pQsxg0qTg/iuDHnMkIpIaSc3Z1mhg4D2O9+rXSRFJnihmuj8BHO/uXwPeZGZNEeyzMqZMCW5ffjneOEREKieROdvMds52q8VERBIoiqL7fODw8P5PgYYI9lkZs2YFt/PmxRuHiEjlJDZn1zWp6BaR5Cq46Dazd5jZP4X3jwRmuvsyAHfvc/cNEcdYPq99bXD797/HG4eISJmkKWdniu7eTb0xRyIiUriCim4z+xJwO3AagLs/AXzIzM4qQ2zld/DBwe2f/xxvHCIiZZC2nN0wMZiU71ndE3MkIiKFK3Sm+03A3sAfstbdD1wdWUSV9OY3Q309PPEErF8fdzQiIlFLVc5u3KsRgO5XumOORESkcIUW3Q+6+zp2vSZY8i4nnDF6NBx1FPT2wq9+FXc0IiJRS1XObtgrnOl+RTPdIpI8hRbdL5rZxcAsM3u/md0BfAq4JvrQKuTMM4Pb738/3jhERKKXqpzdODmY6d62fFvMkYiIFK6gotvdfwo8BNQDHwLWAW939/8qQ2yV8U//BMOHw69/DU89FXc0IiKRSVvObpoRnN1wy/wtMUciIlK4YYW+wN0fAx7LPDazejN7m7v/NtLIKmXCBDjvPLjhBvjc5+DBB6Gudi/UKSLpkqacPfK1IwHY9HddRVhEkmfQotvMfgmMGmIfE4HVQOIS+A6XXAK33AIPPQQ33ggXXBB3RCIiBUt7zh45Oyy6n9lEX3cfdY2aIBGR5Bhqpvuv4XPmA33AucBcoCPrOfsBE8oSXaVMnBj0dJ92WjDbPWMGnHxy3FGJiBQq1Tm7Yc8GRsweweZnN7PhsQ2MO2Zc3CGJiORtqKL7OsDCo98xs9e4++W5TzKzu8oRXEWdeipcfDF8+ctw+unws5/Bu98dd1QiIoVIfc4e/7bxbH52M2vuW6OiW0QSZdDf5tx9fSZ5hyaZWX32c8xsDnBkGWKrvKuugo9/HLZuDYrwa64JTicoIpIAtZCzm/+xGYBVP16F9/oQzxYRqR6FNsT9HHjazL5tZteZ2T3Ag8A3og8tBmZw001wxRXQ1wef/zwcdxzMmxd3ZCIixUhdzh77lrE0zWiie1k3r9z+StzhiIjkrdBTBj4IvJ3gIJxZBKefOt3dv1KG2OJhBpddBvfcA5MnwyOPBJeLP+88WLUq7uhERPKWxpxtZuzzb/sAsPiyxfRt64s5IhGR/BR86Le7L3f3S939FHf/IHC3mR1ehtjidfLJ8Le/BcU2wLe/HRxgedFFsGJFvLGJiOQpjTl78tmTaZrVxJaFW1h02aK4wxERyYu5598TZ2Y397O6Geh191Mji6oIra2t3tHRMfQTizFvXlBs33138LixMbiS5cc+BsccE8yOi4gUDH6l3QAAIABJREFUycyedPfWMuw3tTl7/Z/WM/fNc8HhtT99LXu9b68IoxMRGVixObvQme6ZwFJgSdYyHHiy0DdOlIMOgl/+Eh5/HP7xH6GnB370I5gzBw44AC6/HObOhQL+AyMiUgGpzdlj3ziW/a/eHxzm/fM81ty3Ju6QREQGVehM9yx3n5+z7kDgbHf/fNTBFaKsM925XnwRbr4ZfvhDWLZs5/rp0+GUU+Ckk+DYY2HUUNeoEBEp60x3qnO2u7PwgoUsu3EZ1MMB3zqAKedOwfTro4iUUUVmunOTd2gbcE6hb5xo++8fnF5wyRL49a+Dvu8pU2DpUvjmN+Fd7wouLz9nTnDe744O2L497qhFpMakPWebGTO+PoPpX5gOvfD8x5/n7x/4Oz2v9sQdmojIbgqd6V4EZL+gHpgE/Mjdz81zHycAVwFnuvvicF0TcC3QBYwHLnL3DUNty1bRme7+9PXBE0/A/fcHhfgTTwTrMkaOhCOPhKOOgje+Mbhtbo4vXhGpGmWc6a6ZnL3ylpUs+NQCejf20jCxgf2+vB9TPjoFq9est4hEq9icXWjR3Q7cys4k3gescvfn83z9JOBNwJ3AflkJ/AfAA+5+m5mdBHzY3T8w1LZssRfduV59FR56KCjAf/OboCUl1957w6H/v707j5PsrOs9/vlVV/XePdM9S2Y6mclMSDKEAIk6kHAxZoREkUgwXBIUF8RLBgW85HUBkRgVgSugCRcBIxlEkReLoHBZwuIFCWIgmARiImCSIcvMOJPMvnRP712/+8dzavp0Te1de33fr9fzqlPnnDrnd05VP/Xrp57znIuWlnPPha6u09cVkbZVw6S7o+rsyR9P8tD/eIjj3zoOQP+F/Zz9lrNZ87I1JJJlD9YlIpJTvZLun3H3b2XN20K4Ev7HJW4jASwQVeBmNgY8Boy4+6SZJYHjwIXAbL5lmco/o+mS7mwHDsB3vwt33RUe77kHTp48fb2+PjjvPNiyJZTzz1+cXrGi/nGLSM3VMOnuuDrb3Tn46YM88qZHmNkzA0Dv5l7GfmeMda9YR/fa7qrvU0Q6S6V1drLM9V8AfCtr3n7gDuAnStmAu6ezLnLZBhxy98lo+Xz0k+jlwFyBZY+XGXtjrV0bLrK8+urwPJ2GRx6Bf/93uP/+xbJnDzzwQCjZzjgDNm+Gs88OZdOmpY8DA/U8IhFpfh1XZ5sZa1+2ltXXrGb/x/az+127mdo5xaO/9yiP3fgYq65exdrr1jL6wlGSQ+V+BYqIVK6kGsfMXkfo0zdkZm/OscrXlhHDmUD2WE/jwBgwX2DZEgcPHmTr1sV/OrZv38727duXEVaNJRKhRfu88+DaaxfnHzsGDz8cykMPLZadO8MdMffvDy3luaxaFZLvsbHFsn790um1a9V9RaQJ7Nixgx07dmSerq7mtlVnQ6I7wfrfWs+6V6zj8JcP88SHnuDwlw5z6LOHOPTZQ1iPMfrzo6x64SpGrhyh75y+qu1bRNpPNerskruXmNkzgVcA78taNOvuZd2i0cycxZ8q3whc6+6XxJZ/n9APkXzL3P3m+DabvnvJcqXTYXjCxx8Po6ZkP+7aBbOzxbeTSIQW8zVrFsvq1Uufx8voqJJ0kRqrRfcS1dmnm9k7w4F/OMDBfzzIie+cWHKJae85vYxcOcKKn17BiuesoPecXg09KCI51bx7ibs/YGa/7+6njcVkZhvdfXe5O4/sJVzhHjcE7CNUifmWdZZEAjZsCOWyy05fnk6HVvDdu8Nt6vftCyV7+uDB8FjqrezNQgv66tUwMlJe6e/X3TpFGkR19ul6zuxhww0b2HDDBmb2zXD4i4c58rUjHPvnY0w/Os0Ttz3BE7eFujG1NsXwpcMMXzLM4MWDDF40SPdYtxJxEalYwZbuqEXjy+7+o+j5HxCGnFqyGvCz7r6t5J0ubTVZB/wYWOPuU2aWAiaAcwgX7+Rc5u5749ts+5buapmdDcn5oUMhAS9Wjh6tfF+p1NIkfHh4sQwNLX2eXeLLU6nqHb9IE6pWS7fq7Mr4gjN+7zhHv3GUE3ed4MRdJ5g7dPpY38nRJIMXDTLwzAEGnj5A/5Z++s/vJ7U2pWRcpIPUqqV7M7Ay9vwngQHgSRZ/mDNgXak7tMWayQDc/Ukzux14PnA7cAXwmUwFXWiZVKC7e7HFvBRzc3DkSEjAjx0LSXipZXo6jNpy4MDyYu7tPT0hHxoKd/wcGAiPmRJ/Xmg6qQuopC2pzq6AdRnDl4RWbQgjoEw/Os3xu44zfs84Jx84ycT9E8wfmefYHcc4dsexJa/vGu6i//x++s7vC4/n9dF7di+9m3rpXt+NJZSQi0j5Qwb+JPADd5/Nmv+0TMtKkdcPAb8G3Aq8DfiAux80sxHCzRQeJ3wZvMXdJ6LX5F0W10ytJhKZnl5MwI8dgxMnQhkfX5yOl3zzFxaqH1tPT/7EfGAgdI1ZTunuVtcaKVkNhwxUnV0l7s7M3plTCfjkjyaZfHiSqYenmD+W/47DljJ6NvSEJPzsXnrO7qF3Uy89Z/XQM9ZD91g3yRVJtZSLtJB6jdP9Cnf/u6x5W4Dz3f2L5e68mlqtApcSucPU1NKE/PhxmJgI45xPTFQ2Hb9baC0kEstP3LNLX9/pz9X1pi3UMOlWnV1j7s7coTmmHp4KSfjOKaZ2TjG9a5rpXdPMHSh+S/pEX4LusW561ockPJOM94z10L2um9SaFKm1KVKrU7rJj0gTqOmFlGZ2GeGnxSujMVfjRglXxze0Apc2ZbaYZJ5xRnW26Q4zM7mT8YkJmJxcfpmdXdxeLSWTS5PxXNPVmNfXF/6RkJagOrt+zIzuNd10r+lmxXNPv4HZwtQCM7tnmH58+lQiPr1rmtm9s8zsm2Fm7wzpk2mmH5lm+pHpovtLjiRJrU3RvSZKxtek6F6bezo1miLRo79bkWZRasfWY8BHCf0Fn5u1bBb462oGJVJTZqGfeG9vGJWlFubnQwt9NRL4eDl5cul25+fDrwDj47U5jrienuUn8aW8pqdHXXOWT3V2k+jq6woXXG7pz7vO/Pg8s0+EJHx2X9bjk7PMHZwL5dAc80fnmT86z9RDUyXtPzGQIDWaCsn6aIrkaGmPif6EuryIVFlJSbe7/4eZbSOMv6rKWqSYZHLxgs9acQ8Xuk5OLibi8YS83HmFlk9NhV8HZmaWN6JNKeK/bmT628cvni30mD1veBhWruy4C2dVZ7eW5FCS5FCS/vPzJ+YQRlmZOxIS8NkDi8l4vun5o/OkT6aZOTnDzJ6ZsmKybiM5kiS5Iiorw2PXiq5T06fNjz8f7lJXGJEs5YzTfZw8rSNm9tPufmfVohKR4szCBZvd3SGxrKV0OiTclSb05ST5c3OhRf/kyTBqTjUMD58+jvzo6OL0+vVhRJ+NG+EpT2mLrjSqs9uPdS12ZRl42kDR9d2dhYkF5o/MM3dkrvjj0cXn6ak0c/vnmNtfvE96Pl2DXacl46cl7sMhQe8a6grTmcdoXtdAl1rcpW0UTLrN7AvAYJFtGHAhsLZaQYlIk0kkQtePvr5ws6RaynTNySTeme4zExPlPx4/vnTknF27iu9/9Wq47jp42tNCYp4pq1aFx+Hhpk3KVWdLnJmdakXvPbu3rNcuTC2ErizHo3JsnoXjC6emc86PLVs4scDCRCize0u4W3I+iZC8ZxLyruEukkMhKc81L1/ynhxOkuhTlxlprGIt3Q8R+gbuZckNc5cw4BerGZSIdLBqd81Jp0PCffRoGHM+Ppb8kSOh7NsHe/bAww/Dk0/Crbfm314isdhSvmpVuMB382bYsgVe9CI488zqxF0Z1dlSFV19XXT1ddEz1lPR6z3tLIznSMaPLywm7cfmwzonlj4unFhgfjwk7umpdEjgT1Rh6NgERRP1rsEyykCXxmCXshS7I+VqYNzdC3YGM7NV7n642sGVo12GnxKRBnKHO+6Ab3wDDh9eTMqPHAnPi7WUDw7Czp2wruR7zwBVvSOl6mxpK+n5dEjEM0n5iYWiiXq+5enp6g8Vm+hPlJagl5rQD3SRSDXnL2myqCZDBrr7oRK3cw26Gl5EGmFuLtyIaXo6dEvJns41L9fyTJeWTNeUTIk/L2blyjD6SoOozpZ2k0gmSIwkSI0s/54E6bmsBD4rKc90hylW5sfnWZhYIH0yTXoylFLGYy+V9Vh5Le5FSnIoqaEjm8SyLuk3s43Aa4HXoQpcpDm5hy4W8/OFy9xc8XWWs345r5mbWxwtZWYmjHsef56ZNz1dmzuW5tPfD2vWhLJ6dShnnAHnngsXXgiXXtrUNyxSnS2dLJFKkBgNQyhWg6edhcnSk/VSi8848zPzzB/Of6fTclnKlnajGeo67XlyKKvrTZ71dHFr5cpOus2sC7gauB64EthP6EMo0jrcQ3I3OxsSvLm5xelc8xqdsC43wW1nmYs8e3tPfyxlOj4vPkRh9vTAQFMn1PmozhapDUsYycEkycHqDUnq7qRn0tVL4qOWfZ9z5o/MM3+kCt8HFo1ME0vOkyuTYYjJkSSpkTDe+6npkWQY/z2a7hrq3KS95E+KmW0GXgW8ElgF7AJeBPwT8JyaRCetIZPAZrdE5mqZLJTYFkt8S51XymvaPRHNlkiECxRLLalUbdcv5TU9PYuluzv3897elkyE60F1tkjrMTO6ervo6u2CKt67LT2TXuxKMx51kYn3hc/qdrNkvRNLX5OeXOymM0sFI9N0QXJluAlTanUq3GF1bfepx+4zupfMS42msK72SNKLJt1mdi2hheR5wI+B9wB/B9zo7l+NVvtOzSKUwtxDEpk93nH24/R08aS4lHXyvabABblNq6srJG+p1OJjfDo+r5ykstYJayUJbpMOcSfVpzpbRLIlehJ0r+mGNcvfVno+vaQVff5ENDrN0cWSGfd9/ug8c0fD2O+Z6fTJNPOHQ/eZqZ0l3Fk1Aak1KXrO7KF3Yy89G3pOld4NvfRs7KF7fXdL3IyplJbuKcJtg78DvMrdHwYwsxbMshrEfel4w+PjYQiz+PPMxVqFEud8NxpJV/+K7LIlEostj/EWyezWyVKT3ELzKnlN9jwlotK+VGeLSM0kkgkSKxOkVlb2K2N6Nn0qAZ87NMfcgTlm98+GxwOLj5l580fnT92oaeL7uS9ot5TR95Q++rb00b+lP5QL+hm8eJCuvq7lHG5VFU263f124PboApzrzWwT8Hng1Nk2M/NCYw+2A/eQGB88GIYOO3wYDh1anM48P3Ikd0Jdy9OTSoV+qf39uR/jfVbzJcS5Sjnrd9httkWalepsEWlmie4E3WeEbiSlSM+mmT0wy8x/zTCzZ7FM75k+NT37xCyTD04y+eAkh1kcDdWSxsAzBhi+ZJiV21Yy+gujJIcbl6+Ucxv43cAfmlkS+CXgLDN7O/Bh4OeAHbUJsY7SaXj0UXjgAfjRj2D37lD27AllfLzybff3L97wI1OGh5c+HxwM6+VLnvM9KuEVkSwdUWeLSNtLdCfoPauX3rN64dLc6yxMLjC1cyok3g+FcvI/TnLyhyeZuG+Cifsm2PfBfVjKGH3BKBvetIGVl62s74FQ5OY4RV9sdi7wu8D17t5ftagqUPGNFqam4BOfgM99Dr75zcJj8fb3w9q14S50mbJ69dLnq1aFZDqeUA8OKjEWkbyqdXOcEvbT+nW2iEiJ5ifmmfjeBCe+e4LDXz7M8TuPQ9Qjd90r13H+B88n0V1+V9Oa3BynGHf/MfB6M/vecrbTMHffDS95Cezduzhv/Xq46CJ4xjNg0ybYsAE2bgyPIyPQocPciEjra/k6W0SkDMnBJCsvX8nKy1ey8c0bmd0/y95b97Lnz/fw5N8+SXIkybm3nFu/eKqxEXf/aDW2U1dTU3DNNbBvH1x8MbzmNXDVVTA21ujIRERqqiXrbBGRZeo+o5vNf7KZ0ReMct9/u499t+5j01s3kRyqT2+Ezu3z8O1vh4T7ggtCi7fG+hURERFpeyues4LBiweZ+PcJTv7gJCues6Iu++3cMdN27QqPz362Em4RERGRDtJzVg8Acwfn6rbPzk26T54Mj0NDjY1DREREROoq0R9S4IXJhfrts257ajZz0X823aWNEykiIiIi7SHRE1Jgn6nfLQs6N+kWERERkc4UDUZXz/uEdW7SrZuxiYiIiHQkS2Sy7vrts3OT7gyNuy0iIiLSWTLpX7p+u+zcpFst3SIiIiKdSd1LGkAt3SIiIiIdRd1L6kkt3SIiIiKdSd1LGkAt3SIiIiKdRd1L6kgt3SIiIiId6VT3ErV015FaukVEREQ6Syb969Q+3Wb2DjPzqNwfzeszs/eb2Z+a2V+Z2XCj4xQREdXZItLCogzY0/XLupN121MRZjYArAaujGbtih5vBb7q7p8ysxcAtwG/0oAQRUQkojpbRFqZNaCnQ9Mk3cD1wIPAne4+DWBmY8DLgddG63wd+IyZbXL3xxsSpYiIgOpsEZGyNEX3Egv/bvwycAuw38wyrSLbgEPuPgng7vPAY8DljYhTRERUZ4uIVKIpWro9jNdyqZmNAG8APm5mx4AzgcNZq48DY3UOUUREIqqzRaRt1PFCyqZIujPc/ShwU9SKcgPwNWAqa7UeYC77tQcPHmTr1q2nnm/fvp3t27fXMFoRkcrt2LGDHTt2ZJ6ubmQslVKdLSItq8wu3dWos5sq6Y55H/BNYC8wkrVsCNiX/YI1a9Zw77331j4yEZEqiCeZZnaoweEsl+psEWlr1aizm6JPdw5p4PvAHcCYmfUBmFkK2Aj8SwNjExGRpVRni4gU0RRJt5mNmdnLzSwR/Uz5BuAP3P1J4Hbg+dGqVwCfcfe9jYpVRKTTqc4WkbbRgX26R4B3ATcCdwLvjQ0v9TvAn5nZRcA6QJ3+REQaS3W2iLS2BtyQvCmSbnf/IeEnyFzLjhLGgxURkSagOltEpHxN0b1ERERERKSdKekWERERkY4UbjtQH0q6RURERKSzNKBPt5JuEREREZEaU9ItIiIiIlJjSrpFREREpDPVcZxuJd0iIiIi0lHCfb3qS0m3iIiIiEiNKekWEREREakxJd0iIiIi0pnUp1tEREREpEY0TreIiIiISPtR0i0iIiIiUmNKukVERESkM6lPt4iIiIhIjahPt4iIiIhI+1HSLSIiIiJSY0q6RURERKQjudevU7eSbhERERHpLOrTLSIiIiLSfpR0i4iIiIjUmJJuEREREelMGqdbRERERKQ2zOrfqVtJt4iIiIhIjSnpFhEREZHOpO4lIiIiIiI1oiEDRURERETaj5JuEREREZEa69yke3AQ1q+HoaFGRyIiIiIiddQ13EX3+m66Brrqts9k3fbUbF772lBEREREpKNsumkTm27aVNd9dm5Lt4iIiIhInSjpFhERERGpsZboXmJmfcCfAePACPBmdz/R2KhERCQX1dkiIqdrlZbuW4E73f1G4PPAbQ2OR0RE8lOdLSKSpemTbjMbA14OfDGa9XXgajPbVI3t79ixoxqbqSvFXD+tGHcrxgytGXcrxlxrqrNLo+NoLu1yHNA+x9IuxxHX9Ek3sA045O6TAO4+DzwGXF6Njbfim6qY66cV427FmKE1427FmOtgG6qzi9JxNJd2OQ5on2Npl+OIa4U+3WcCh7PmjQNj8RkHDx5k69atp55v376d7du31z46EZEK7NixI/6lsrqRsVSZ6mwRaTvVqLPN3asXUQ2Y2RuBa939kti87wOfcPebY/MOArsq2MVq4NCyA60vxVw/rRh3K8YMrRl3tWI+293XVGE7Dac6u2Q6jubSLscB7XMszXwcFdXZrdDSvZdw9XvcELAvPqNdvrBERFqc6mwRkRxaoaV7HfBjYI27T5lZCpgAznH3vY2NTkRE4lRni4jk1vQXUrr7k8DtwPOjWVcAn1HlLSLSfFRni4jk1vQt3QBmNkK40cLjwDrgLe4+0dCgREQkJ9XZIiKna/qWbgB3P+ru17v7/3b3361G5W1mfWb2fjP7UzP7KzMbrkas1WBm7zAzj8r90by88TbqWMzsSjP7t/j4u5XGWa9jyBVzNP+0c95EMb/EzB4zs8Nm9hdmlozmrzaz28zs3Wb2nuhnfJazrNYxR8s+FjvXn2+WmKP9XG5mD5nZUTN7b2x+U3+um02n1dnlylfftIJy6/1mVe53QTOq5LuhGVXyfdGy3L0jC/C3wMui6RcAn2x0TFEsA8AHCT/JXgGcVyzeRhwLcAZwDeDAplJiafQxFIg55zlvkpg3Ah8Ffgr4DeAk8MZo2R3AJdH0bwPvjL2uomV1iHkj8J7Yud7QDDFH210B3Ei4CPAqYA64otk/151S2uVcFqpvmr0UqENb6r2p5Lug2Uql3w3NVir9vmjV0vAAGvQmjwEzQH/0PBm90ZuaILYbotJbSryNPBbCLyWnKq1K46znMWTHnO+cN8t5By4DkrHnfwZ8CbgU2BebfybhYrWBSpfVOuZo+v8AL40vj+Y3NOZou31EXe6i598DfrYVPtftXtrpXOarb1qlZNehrfrelPNd0Iylku+GRsdcznFE0zm/L1q5tET3khrYRg3vmFYpMzPgl4FbgP1m9ivRom3kj7fQsppy93TWrErjrNsxZMdc4Jwv53iqGe+/RtvO2AvsBp4X7S+zXuYita3LWFbTmM2sH3gR8Ono+fNj6zQ05mi7Ux7V9GY2APwn8E1a4HPdAbbRBueySH3TEsqs95tWmd8FTafC74amU+H3Rcvq1KS7pDum1ZsHlxIGhH8/8HEz+wUKx9tMx1JpnA07hgLnnCJxNSrmZxF+/mylc/0s4IPuPunu50YxfAn4ipk9I1qnaWI2s8uArwApoH8ZsTXT32ara4tzWaS+aVV6b5pDKd8NraCU74uW1Qo3x6kFB6ay5vUQ+nA2nLsfBW6K/vO+AfgaheNtlmMpdl6b9hhynPOvsLzjqTozOxc46O73m1mh2CpdVnXxmDPz3P0J4PqoJeN1wKubKWbCGNMfIfzM+efAo0X23zSfkTbW1HV2ufLUN61K702DlfHd0NTK+L5oWZ3a0l3SHdOawPsIFxIUireZjqXSOJvpGDLnHJooZjPrIlQ2v58rtugLYiBXbGUsq3XM2fKe60bFDKGSd/e/Ad5I+Im8HT7Xra5dz2X8b6BV6b1poDK/G5pWmd8XLatTk+47gDEz6wOIhtPZCPxLQ6M6XRr4PoXjbaZjqTTOZjqGzDmnSFz1jvmNwM3uPhM9/xqwJbZ8EzAL3LuMZbWOOVv8XDdLzHH3Er7A2uFz3era9VzG/wZald6bxirnu6GZlfN90bI6Mun2Jr1jmpmNmdnLzSwR/Xf6BuAPCsXbyGOJYgQwKHxem+UYsmPOd86Xczw1iPlG4D6gz8zOMbPfAo4BT5jZ06PVXgC8391n3f3eSpbVOmYze7qZXR0t7wZeQejCQZPE3G9mPxGb9ULgL1rhc93u2uVcFqpvWkU59X4DwitZOd8Fzarc74ZGxVlMud8Xrawl7khZC9aEd0wzswsJ/cdOAHcC73X3B6NleeNtxLGY2RDwa8CtwNuAD7j7wUrjrMcx5IoZWEuec94kMd8EvD1r9oPufoGZnQ38UbT/QcI/aPPR6ypaVsuYCf3xPkm4sv67wLvdfV/sdQ2LOdrHswkX7DwYxXe3u/9DtKxpP9edoh3OZaE6vhVUUu83o0q+C5pNpd8NzabS74tW1bFJt4iIiIhIvXRk9xIRERERkXpS0i0iIiIiUmNKukVEREREakxJt4iIiIhIjSnpFhERERGpMSXdIiIiIiI1pqRbRERERKTGlHSLSFMxszVm9gUz22Vmv9/oeEREpDDV26XRzXFEpKmY2c8DdwDDhDuTne/uRxoblYiI5KN6uzRKukWkaZnZt4DnNestjEVEZCnV2/mpe4mINCUz2wD8P1XcIiKtQfV2YUq6perM7PfM7C8bHUdGNeIxsyvM7C4z21alsNqKmd1c5e0ZcB3wrmpuV0SWUn29fGa2ycw+ZGa3VWl7XWZ2vZntKmHdq8zs3yo91mrW3aq3i1PS3WLM7AIz+6yZ/bWZ3W9mbma/1ui4snwZ+GSjg4ipRjx7gIuqEEtTMrMXmtkhM3vCzJ5jZhuiz5mb2QfMbGW03rPN7AdmdouZ9UXzzgCerHJILwc+7O7z0fZFWo7q64q0Yn19IHrsqdL2EsAxYGMJ694LXJBvoZldaGafij57N5nZutiyatfdqreLUNLdej4H/Im7vwr4CeCjmQVm9kwzu7wRQcX37e4/cPc7GxFHLuXEk+8cuvtDwKFqxlXN92u523L3LwO3AvPufpe77wGuB2aBf3P3Y9F6dxMq+Te5+1T08hcBty/rAGLM7A+BtwJfN7MHgS3V2rZInam+LlMr1tfuPgnsq9a+3H0O+F6J6+4nJOj5lv8Q2BE9vdnd40l21epu1dulUdLdQsxsFDgfOAng7mngjwA3sxWECt0aEFfD9l1NJRxHuo77asS2PgWcZWY/BeDuh4F/Bn4xtq+VwL7os5fxVHd/cJn7PsXd3+7u57n7Vnd/qrt/q1rbFqkX1de11YT1dbVHpSgn/mLrLkSP2f2sq1Z3q94ujZLuFhINv/MQ8FUzuzSatwu4m9CPajOw3cxeCWBmbzOzP4x+WjrVb8vMtprZe8zsBjObNLOHzOx/mdkvmdk9ZnaNmf1fMzthZq+OvS7n9rL2fb2ZfSTety36ifX9ZvZOM7vDzJ4dzS+4v9jrzcxebWY/NLMrzewBM9tvZr9Zwj7Oi8dTZJ+nncMcRqKfi8fN7Pqsc/rW6Nx8NqqkSz1np+3LzF5qZneb2W+a2U4ze1eB7eV673PGU0jUIvJD4NrY7DXAVRZ1JQGuBr4Qi7MPmIw9f6qZvcvMbjazf7HQTeVlZnavmb3YzD5pZsfM7PVm9lPRMR42sxcWi0+klai+7pz6Opd8x1hoP2a22czAhPglAAAGnElEQVT+xkKr8Z+Xsp8c+/2ymX3czJ5bZD3V3Y3g7iotVIALgV2E/2z/BlgbW/Y4sC2afjrwaDS9hvBf+CrCf+q7ga3Rsr8EvhNNdwNHgA8Qxtr8FWBvoe1l75vQp+024CPR/F5CIrcyev5S4CAwVGh/OY57XbTPG6N13wFMA2cX2Ud2PAX3GT+HOWJ4HLgF6I9etzuaPwz8fWy9bxMqzJLOWZ59jRJ+Mnwv8GzgikLby3rvc8ZT4ufrJuCRaPopwDuBGeAlsc+Lxda/Grgsmu4H7gH6o+ffis7DMHAceBvQB7wsen4d4R//NxB9BlVU2qmg+roj6uto+VtLPI+F3ut7gAtir/ESP2eZ9/M5wB/nWL4t2k8yNk91dwOKWrpbjIfWyGcQKqHfAH5oZrkuGNkJ/LqZpYDLonmDhD/yDUCmP+5/Ev6YcPdZYBz4R3c/QfgjXFdke9nxzbD0woyrgCmP+gQDnyX8sb6kyP6yt5vZ5ieidd9B6G98ZZF9LImnnH3m8UUP/ffuAdbHjnFV1BJ1A+GcQonnLM/xHiFUbp9397vd/etlbC9fPKX4FHCOhS4m/x34EPBPwHVmNgiMe1QrR54LfCeavprwxZZpPbkGeHt0no8B3/DQD/xeYNjdP+3hJ/f7KO89EGkJqq87o77OIe8xFtjPFcAKd8/Ec2+Z+3wW8Kvu/iclrq+6uwGSjQ5ASmdmXcAGd38c+J9m9mHCld47gEvi67r7jJmdCbwZ+LvMJtz9gJndCzyf8J/4ZkKideqlsel5oi5I+baXJ9T4Ns4DUrG40mb2OHBWof0V4+7TZvYosBJYXcY+Kt5nlnkW/342Ao+5+3uzVyrlnJnZrxOS24wtHn6G9nisZbwHeeMpxt13mtl9hC4mK9z9UTP7FKH16VrgS7G4M5+NTH/BswktU5ltHY5vOja9wFJp1NVN2ozq6yXH1wn1dVze81jgvXkm4R+MSj0POM/MBt19Is866egYVHc3iE5Wa0kR62/r7vcTfr57avaKZvY84DXu/g4PI1HEvRZ4loWhqx4A3l1sx0W2V8huQstpfCglI/R1XK5u4OEa76MUT7C03zNmdkkZ5+wLwMWxkvMq+DK2lzOe0g+HvwdeCTwai68LeD3hp9iMSwj9UzP2Aj9tZgOx/RbsVyjSxlRfL9VJ9XXeYyywnxPRa7qpzM2EXwn+Is/y4754Abzq7gZR0t16ftfMxmLPzwS+Gk3PEi4ceSpheKoVZtZjZldEy9ea2WpC38IPE34a+g7Rz5WRBFn/3ZuZFdle9r4tto3PARPAL0fbWhntL9Nimm9/+ayP1jkzeu1XS9hHPJ5i+4wfR7ZcrR5GaL0aAG43s5+3cMFOsffg1L6A9e7+YKzMxXbRFZsutL143DnjMbNuM9sRdR0p5NOEn7X/EcDdx4GvAN/2paOW/Byh60lG5nx/0sJY368DMmO1Zr8H2e9zS4+kIJKH6ms6pr4u9Tzm3A+hLu0m/GMG4VcNzGyshLrbCK3Ovw5ca2bXRa99g5ndRGjl/+fY+qq7G6UWHcVValMIF2c4cBT4GPARws9TI9HyPwb+C3gx4eeinYSWhasILSQfIVRgXyBc1DEbbW+S8Ef4C8Ac4WKdDcBbouW/VWh7Wft+PfCvhH5yF0bLthJaSG8h9G18ejQ/7/7yHL9Hr389obvDM2LL8u1jUzyeYvuMn8Osfb+YcEHhB6JzkXndK6PllwM/iN6bW6LzXOo5e3GOY31ptL9PAmdG8wq9p0u2lSeeswg3cfhECZ+192U9fxnws1nzbs7xup8BfgQcBm6K5r2Q8Fl7PzBGNGwa8NuEL+UPRcf6i43+G1NRqVZB9XUn1dcbgTsJrfUXFznGQu/1VYRfGL9LuKj9YeA10bHnrLsJv6bMEVq4R4FvEBL+VwGvJoxXfguha0vmNaq7G1QsOrHSIcxsLfDb7v626HmC0Pryane/qaHBFWFmDmz20EdSKmRmb3L3ioajEpH6UX0tcaq7W5+6l3Se1xG6GozG5j2P2AVyzSj2c5Z+yqqQmXWZ2a8Cn2l0LCJSEtXXorq7jSjp7jwfJPT9+k8z2wN8Dfgvd7+rsWHlF11Y8vro6W9aGLpOyjcIfM7dHy26pog0A9XXAqq724a6l4iIiIiI1JhaukVEREREakxJt4iIiIhIjSnpFhERERGpMSXdIiIiIiI1pqRbRERERKTGlHSLiIiIiNSYkm4RERERkRpT0i0iIiIiUmNKukVEREREauz/A9F2cf9Werm3AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig = plt.figure(figsize=(12,8))\n", "plt.rc('font',family='Times New Roman')\n", "params = {'mathtext.default': 'regular' } \n", "plt.rcParams.update(params)\n", "\n", "plt.subplot(2, 2, 1)\n", "plt.plot(vehicle1.v_kmsc, vehicle1.h_kmc, 'k-', linewidth=2.0)\n", "\n", "plt.xlabel('Speed, km/s',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 2)\n", "plt.plot(vehicle1.acc_net_g, vehicle1.h_kmc, 'b-', linewidth=2.0)\n", "plt.xlabel('Deceleration, Earth g',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.subplot(2, 2, 3)\n", "plt.plot(vehicle1.q_stag_total, vehicle1.h_kmc,'r-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-rate, '+r'$W/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "\n", "plt.subplot(2, 2, 4)\n", "plt.plot(vehicle1.heatload/1.0E3, vehicle1.h_kmc, 'm-', linewidth=2.0)\n", "plt.xlabel('Stagnation point heat-load, '+r'$kJ/cm^2$',fontsize=14)\n", "plt.ylabel('Altitude, km', fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.savefig('../plots/decadal-2010-venus-flagship.png',bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-venus-flagship.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-venus-flagship.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAF4CAYAAACYdLRaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZzV8/7A8de7VdpEJdGGEkU3xi5FIrJfO5FtcC2Xnz0url1kufa5dMkaImS5pCxZUqJLJGnTpqm0UKqZef/+eH+POY2Zmjnne873LO/n4/F9fL/n+z3n+33P1Mx7PruoKs4551yYakUdgHPOudzjycU551zoPLk455wLnScX55xzofPk4pxzLnSeXJxzzoWuTtQBZILmzZtr+/btE/pscXExLVq0CDcgt0H+fU8//56nX6Z/z7/44otFqlppgJ5cgPbt2zNhwoSEPltQUJDwZ13i/Puefv49T79M/56LyKyqrnm1mHPOudB5cnHOORc6Ty5JKiwsjDqEvOTf9/Tz73n6ZfP3XHxuMSgoKNBMrtd0zrlMJCJfqGpBZde85OKccy50nlycc86FzpOLc8650Hlycc45FzpPLs4550LnycU551zoPLk455wLnc8t5pzLbKqwZo0d16oFtWvb3mU0Ty7OuWisWgUzZsCsWTBzZvl+3jxYvty2ZctsX1Ly58/XqgWNG8Omm667tWoFHTpA+/bl+yZN0vu1ucxKLiLSEygCWgJPqurFwfkGwCBgBdAMuFJVlydzzTmXRsuWwVdfwZdfwsSJtk2ZAqWl1ft83bogYu8vK7PSTFmZ3XfZMktS69OyJey0E3TrVr7ffnuoVy/5r81VKmOSi4g0BfYG9gD2AkaIyEhVHQU8BLytqsNEpC/wKHBi8NFErznnUmX5cvjoIxg9Gt57DyZN+vN7ateGbbe1kkX79tCune232go22cRKG02aWOmkfv11P6tqiWbFCliyxLbFi22bN89KQDNm2H7mTFi4EEaNsi2mfn3YdVfo0QP22Qf22sue60KRMXOLBaWM3zUISES+AC4DvgdmAM1UdaWI1AGWAV2ANYlcU9WZ8c/2ucWcS5KqlUpGjIB334Xx49ctldSrZyWGnXeG7t1tv+OO0KBBemKbPdsS3P/+Z/tJk+CHH9Z9n4jFeMghcNhhsNtulgBdldY3t1jGlFxUdVXsWEQaAt8B72MljUWqujJ4X4mIzAB6AmsTvDYzXV+XczmrtBQ+/hheftmSyqy4daNq17aSwP7727bnnrDRRtHEKWKlonbt4PDDy8//8gt88omVsD76yBJiLPHcdhu0aAH9+sERR8DBB/+59OTWK2OSS4yI9ABuAeYDGwNbAosrvG0F0BooSfDaOoqLiykoKE++hYWFWT3VtXMpo2q/hJ94Al56CYqLy69tsQUceaT9Qt53X6vOymTNmlms/frZ61WrYOxYeP1122bOtK/ziSeso8DJJ8Ppp1vJK48VFRVRVFQUe9m8qvdlTLVYjIhsARyMNcS/AEwHjlXV3ePeMxF4NnhZ42uqelf8M71azLkNmDsXnnoKnnzSGuJjtt0WjjrKtt13z50uwqowebIlmWHD1m0z6tYNzjzTEk2jRtHFmAGyolosRlXnA0NEpAy4HPgI6+kVrzEwD9AErznnNqSkxH65PvqotaOUldn5li3hlFOgf3/7RSsSbZypIAJdu9p29dXWnjRkCDzzjCWaiy6CG26w/YUXWsnGrSPjSi4xItIVuBs4FZgGtFDVVSJSF/gV2BooTeSaqs6Nf5aXXJyLs3AhPPYYPPII/PSTnatXz9orTjsNDjrIugbno9Wr4dVX4Z574LPP7FzDhnDOOXDppdD6T7XuOS0rFgsTkY1FJL4y8xDgPlVdAIwEegfnDwCGq+rcRK+l+mtxLiuNH2+lkTZt4JprLLF07Gi/SOfPhxdfhEMPzd/EAtaof9xx1hFgzBg48ED47Te4+277Xt18s7XduMwpuYjIbsAbwBTgM+BzVX0xuNYMa4OZCbQCrlbVX5O5Fs9LLi5vqcLbb8Mdd8AHH9g5EUsi558PffrkTjtKqkycaEnllVfsdfv2MHiwtUPlYpVhnPWVXDImuUTJk4vLO2vXWkP1oEHw9dd2rkkTKCyEv/3Npk1xNTNmDPz97+Xfz/33t6rFjh2jjSuFsqJazDmXBqtXw8MPWy+v/v3tF+EWW1iSmT0b7rzTE0ui9tvPSjEPPmgN/KNHW7flIUOshJhnPLk4lw/WrLFeXx07Wslk9mzo3Bkef9ymSbn8cmjaNOoos1+dOvb9nToVjj/e2mPOPBOOPdamqMkjnlycy2W//w7//jd06gTnnmuN9F27WuP85Mlwxhk+8jwVNtsMnnsOhg61waTDh9vUMh9+GHVkaZNx41ycc3HWrLFSxsKFNhp+4ULbFi+2dpOSEttKS+310qX2vkWLbL88bhLwHXaA66+HY47xRvp0ELGqx332sXFBn3wCBxxgpcX+/aOOLuU8uTiXCVRh2jSYMAG+/Ra++872P/xQ+Vom1VW7NnTpYgMBjz3WJ2KMQocO1hPv8svh3nvh1FPhxx8t0edwbzJPLs5FQdV+wbz/fvk2t5IhWCLQtq0tgNWypU2m2LIlNG9u1Vm1a1s9f506dty0qb0ntjVt6qWUTFCnjo0X2mYb61H2z3/av/9jj+VstaQnF+fSRdV6E73wgrV5VFzgqnlzm0m4a1erwtphB9huO9h442jideG74AIryRx/PDz9tFVdjhgR3YzRKeTJxblUmzQJnn/eksr06eXnN9sMeva0Lqy9elky8VJG7uvXz6b4P+gg+O9/4a9/tWULcqwE48nFuVRYvdpKJw8+WD4HFVj11l//alOI7L23t4Hkq+7dbRzMfvvBm29ae9hLL+XUssueXJwL06xZNir7scesxxZYu8eJJ1pVSI8enlCc6drVll3ef3+bffr44610myNzt3kZ3LkwTJ1qMwZvsw3cfrsllm7doKjIGuofftiqvjyxuHjdutlyBptsYm0v55+fM6P5Pbk4l4xvv7UVCrff3gbMAZx0ki3/++WXcPbZNiW7c1XZeWd46y1r1P/3v23SyxzgycW5RPz4I5xwglVtPPuslUjOPttKMM88Y72+cngMgwvZHnuU/3FyxRXlMyxnMU8uztXE0qVw2WVWUhk2zOrHzzvPBjsWFcHWW0cdoctWxx4Lt95q1WInn2wDarOYJxfnqmPtWnjgAZtNePBgGzU/YICNqn/oIWjXLuoIXS646io4/XRbcOyoo8o7hWQhTy7ObciHH1rD64UX2pxePXvaX5X/+Y+t2uhcWESst+Eee8CcOTZVTFlZ1FElxJOLc1VZsgTOOsuSyXffWanllVdsUaidd446Oper6tWzKtdNN7WG/kGDoo4oIZ5cnKtI1Rrpt9/eZrCtV88mGfzmGzjySG+od6nXti089ZQdX3ttVk7V78nFuXjz59v68SefbFPb77uvTd9yww05Nz2Hy3CHHGJtMKWlNgg3yxYb8+TiXMzw4bDjjjYdxyab2Cj7MWNsxUbnonDTTdatfd48uOiiqKOpEU8uzi1bZqPrjznGGuwPPNCqwM480yeSdNGqUweefNJmxn7mGfsDKEv4T47Lb598YsvPDh0KDRpYd+O334Ytt4w6MufMttuWN+qfe65V12YBTy4uP6nC3XdbT7DZs2HXXW26lvPP9wZ7l3nOOw9697ZxL+eckxXzj3lycfln6VI4+mi49FIbDHnppTYX2HbbRR2Zc5WrVQuGDIHGjW2CyyyYHsaTi8svX34Ju+xiP6BNm9oP6V135cw05y6HtW0Lt91mx3//O6xYEW08G+DJxeWPF1+0BbqmT7dBkBMn2rgV57LFuedaFe6cOdY9PoN5cnG5TxX++U9b/XHVKjjjDKsG80kmXbapXdumh6lVC+67D776KuqIquTJxeW2lSttavwbbrAfyLvvtvErG20UdWTOJWbnnW2eu9JSa+jP0LnHPLm43PXzz9Yb7IUXrCH09dfhkku8N5jLfjfeCFtsAZ99Bs89F3U0lfLk4nLT9OnWvjJhAnToAJ9+atNpOJcLmjSxtV/ApohZuTLaeCrhycXlnq++ssTy44/Qvbslli5doo7KuXCdeqpVkc2ZYz0eM4wnF5db3n/fqsIWLID997fXm28edVTOha9WLbjnHju+4w6YOzfaeCrw5OJyx8iRcNBBsHy5LRn75ptWfeBcrtp3XxsQvHIl/OMfUUezDk8uLje8+qr9kK1ZA3/7mzVy+hT5Lh/ccYd1UR461NoaM4QnF5f9XnnFZjReu9Z6gz3wgP2wOZcPtt0WTjnFuibHRvBnAE8uLrsNH26DI0tK4LLLYPBg72rs8s/AgdYG8+STMGtW1NEAnlxcNnv5ZTj+eEssV15p05J7YnH5qFMnW61y7VqrJssAnlxcdnrnHRt5X1oKV19t1QGeWFw+u+Ya+xl4/PGM6DmWUclFRI4WkRkislhE7hOROnHXnhYRDbZX4843F5FHReQOEblbROpW55rLYp98AkcdZX+lXXQR3HKLJxbntt8ejjjCOrVkwJT8GZNcRKQtcCRwDHAJcBZwcdy1hUCfYLsg7qMvAkNU9UpgKnBjNa+5bDRpEvTrZ10vTzvN+vl7YnHOHHaY7d99N9o4ANEMWdFMRHoAn6pqSfB6ENBFVfuJyD3Ax8CI2PXgPXsAL6tq6+D1lsD3wObAjlVdU9Xf4p9dUFCgEyZMSPnX6JI0bZqNvF+40EouL7xga4w758xPP9m6L40bw+LFKV+nSES+UNWCyq5lTMlFVT+KTxzAXGC2iGwMHAa8ELzuHfee/YEZcfeIVTQWbOCayzbFxdC3ryWWPn1sHIsnFufW1aYNdO5sC4mNGxdpKBmTXCqxK/CIqq5U1W2BLYE3gLdEZMfgPVsCiyt8bgXQegPX1lFcXExBQcEfW1FRUZhfh0vWqlVw+OHlc4UNH+4DJJ2rSp8+tn/nnZTcvqio6I/flUDzqt5XrT/9RGQTYLmqlgWvWwKdVHVsGMFW8rxtgWJVnRQ7p6rzgbODkswFwDmAAqsqfLw+sHYD19bRokULvFosQ5WWwskn29TibdvaFC+NG0cdlXOZq08fuP9+a3e5Mfxm5sLCQgoLCwEQkUVVvW+DJRcRuRL4ARgtIgNFpJaqLgRS8ue9iNTGEsdVVbzlX0Db4Hgu0CzuswI0BOZt4JrLFpddZj1fmja1ucJa/6ng6ZyL16uXVRl//jksXRpZGNWpFjsX6KyqvYChwN9FpCnw23o/lbjLgLtUdXUV18uAicHxu8B2cdfaA2uACRu45rLBQw/Bvfdao+Qrr/i0+c5VR+PGsOeetkLl6NGRhVGd5PI/YBmAqs5R1XuAfkDo082KyEDgS6CBiGwtImeISFcROTy4Xg84DRgUxDMBmC8iXYNb9AXuV9U167sWdtwuBcaMsTEsYMsS77dftPE4l00OPND2EXZJrk6by5nAflhJAABVfVZEtg4zEBG5FripwukpWPvKYyIyA/gMuFVVl8W953jgOhGZCTQCrqnmNZeppk+3KfNLS+GKK2xRJOdc9R1wgE3BP2ZMZCFkzDiXKPk4lwyyYgXstRd8840tS/zaaz7DsXM1tXYtNGsGv/0G8+bBFluk5DGRjHMREV+w3NVMWRn072+JpXNnePZZTyzOJaJuXRtwDPDBB5GEkHRyEZG+IjJWRKaJyPRgmwkMTz48l1duvdUW/dpkEyuxNG0adUTOZa9evWz//vuRPD6MIc5PAbdhDfFlwTkBTgrh3i5fvPMOXHedzRP23HPQsWPUETmX3WLJJaKSSxjJZbqq3l3xpIhMCeHeLh/Mng0nnQSqcMMNNs2Lcy45BQWw8cYwZQosWACtWqX18WG0uZwnIsdUcv7wEO7tct3q1bZE8eLFcPDB1sPFOZe8iNtdwkgue2NdhUvjtjLg4RDu7XLdJZfA+PHQrh089ZQt1eqcC0eE7S5hVIvdCBRiI99Lg3MC9A/h3i6XvfgiPPww1Ktnk1FutlnUETmXWyJsdwkjubwFvFpxuhYR+XcI93a5avp0OOssO777bthll2jjcS4XFRTYDOLffQdLlsCmm6bt0WHUQTwPXC0i+8ZtPYGrQ7i3y0Vr1sAJJ8Dy5bbo19/+FnVEzuWmevXK/3BL8/ouYZRc/gF0BAZg09yDVYu1Av4ewv1drhk4sLyd5fHHfZli51Jpzz3hk0/g00+t00yahJFc7gTeUtUV8SdFJH1fhcseb70FgwfbyPvnnrMpKpxzqbPnnrb/9NO0PjaMarFPgGtFpBGAiOwlIkeq6lsh3NvlkuJiOP10O7755vL/9M651In9nI0bZ5PBpkkYyeVpYHegLoCqfgJsGUyf75xRtQb8n3+2HiyXXx51RM7lh9atrQp6xQr49tu0PTaM5DJeVXup6i9x574ALgnh3i5XPPZY+XxhTz7pE1I6l04RVI2FkVxKRKR+7IWINMB6iv0cwr1dLpg6FS6+2I4ffhjatl3/+51z4YoguYTRoP8E8JGIzALqAftgvcV8+hcHJSU2jf7KlXDyyXDiiVFH5Fz+2WMP23/+edoemXRyUdXvg3EthwPtgBeAN1R1abL3djngzjvtP3SbNvDAA1FH41x+2mknq4qeMsX+0Nt445Q/ssbJRUTuBOYDI1V1KoCqrgKGhRyby3Zffw3XX2/Hjz9u67Q459Jvo41ghx3sZ/Lrr2H33VP+yBq3uajq5VgPsR4icqeI3CAie4v4SDgXZ+1aOO00259zDvTpE3VEzuW37t1tP3FiWh6XULWYqi4EHgcQkY2APsBtIqLA58B/VXVlaFG67HPrrfDll9C+vVWNOeei1b07DB1qP5dpkFByEZH6sYkqVfV34HXg9aD0sjtwmYg0UtUrwgvVZY3//c8GSQIMGQKNG0cbj3OuvOSSyckFuBn40yg4VVXgs2Bz+ai01AZLlpTAeefBfvtFHZFzDuAvf7H9119bdXXduil9XKLjXPqKyBARuVlEdgs1Ipfd7rvPJqXcaiu4/faoo3HOxTRtCttsY6u/fvddyh+XaHI5Q1XPwCat3E5EHhCRW0Vk7xBjc9lm+nS49lo7fuQRaNIk2nicc+uKVY199VXKH5VQclHV8cF+mao+paoXALcAHUVkoojcH4x9cflC1XqFrVplAyX79Ys6IudcRV262D4NJZdEG/SbqOry4Lg+cDBwHHAYUBuYAbQMK0iXBZ5+GkaNsqWK77sv6micc5XZfnvbZ2pyAW4WkXHAoUA/bEbkd4BzgNdU9deQ4nPZYMkSuPRSO77rLmjRItp4nHOV69zZ9lOmpPxRiSaXC4CzgdHAFGAUUKSqM0OKy2WTgQNtrZaePW3gpHMuM3XqZCu/Tptmy43Xq5eyRyWaXF4BzlbVJQAisjVwrIi0B6YDL6nqrFAidJnt00/h0UetW+PDD/uSxc5lsgYNoEMH63wzbZpNCZMiifYW+3sssQCo6nSs9LIcG//yo4jcG0J8LpPFxrIAXHZZeX2ucy5zpandJdHksiWAiHQPuiD/AEzApoEZDGyjqheHFKPLVI88ApMm2RQvsS7IzrnMlqZ2l4SrxURkLdAGmIyt6TJMVaeFFZjLcIsXw3XX2fE996RlCm/nXAjSVHJJNLnUBf6NJZT0LcrsMsf118Mvv0Dv3nDEEVFH45yrrm22sf3MmSl9TKLJ5VRVfSvUSFz2+Ppra7yvVQvuvdcb8Z3LJu3a2T7FySXREfqeWPKVKlxyCZSVWWN+165RR+Scq4mttrI/DOfNs+7IKZJog77LV6+9Bu+9B82awT//GXU0zrmaqlsXWre2PxTnzEnZYzy5uOpbvRr+7//s+MYbbaoX51z2ad/e9rNSNxwxo5KLiBwtIjNEZLGI3CcidYLzzUXkURG5Q0TuFpG6cZ9J6JpLwL332uCrLl3g3HOjjsY5l6hYu0u2JZdg/MvTInJpMGq/Op9pCxwJHANcApwFxMbKvAgMUdUrganAjXEfTfSaq4n588tXl7znHqiTaF8Q51zksjW5AOcCB6jqYGAvEWlQjc+0w9aJ+UJVhwIPAvuJyB7Adqo6Lnjf68CFItIw0WshfY35ZeBA+PVXOPxw6NMn6micc8nI4uRyAbBLcPw8Ni5mvVT1I1UtiTs1F5gN7I9N4R9739zgsCCJa64mvvoKnnjCGgIHD446Gudcslq1sv3PP6fsEUknFxFpHLRp3BG83g04IfbLXFXLYmu/1NCuwCPYVDOLK1xbAbRO4pqriauvtv3558O220Ybi3MuebFlMYqLU/aIMEouTwB7A80AVPVzoL6I3JnoDUVkW6BYVScBCqyq8Jb6wNokrq2juLiYgoKCP7aioqJEQ889778Pb78NjRvDNddEHY1zLgwtg7UcE0guRUVFf/yuBJpX9b4wWmXXquo+IvL3uHOzgTuwGZJrRERqY4uOXRWcmgt0irsuQENgXhLX1tGiRQsmTJhQ01BznypcFfwzXH45NK/y/5FzLpskUXIpLCyksLAQABFZVNX7wii5xGY/0+Bh9YErsF/uibgMuEtVVwev3wW2i7veHliDzcKc6DVXHSNGwLhx9lfOJZdEHY1zLiyNG9tCYb/9BqsqVvCEI4zkMlJEngP6icgjwI9YY/75Nb2RiAwEvgQaiMjWInIGsBSYLyKxeUb6Aver6hpVnZDItUS/0LxSUmI9xAD+8Q9o1CjaeJxz4RFJebtL0tViqvqFiBQC/YC2wNXAm6pasTF9vUTkWuCmCqenqOoQETkeuE5EZgKNgPjK/0SvufV58klb72HrrSEoAjvnckiLFjB3LixcCG3bhn77pJOLiJymqk9iXY5j57YTkb1U9fXq3kdVbwZuruLaLODMMK+59Vi1Cm64wY5vuiml62w75yISa0NdXKNyQLUlnFxEpAcgQB8RmVHh8qbAv7CBiy7bPPigTWj3l7/ACSdEHY1zLhUaN7b9r7+m5PbJlFyWAkOBDlhX5HhrgMeSuLeLytKlcOutdnzbbTY1t3Mu98TaUTMtuajq1yLSEzhOVddJJCKyWU3bXFyGGDzYVpjs2RMOOijqaJxzqZLi5JLUn6WqulxVHxORZiLSRkTaBhNQdhCRISHF6NJlyRK47z47vvVWX2HSuVzWMJhm8bffUnL7MBr07wQuxca5xH4bKfB1svd2aXbPPbBiBRx4IOy1V9TROOdSKVOrxeK0AVoCPYDxqjpHRA6gGpNVugwSX2q5/vpoY3HOpV4mV4sFvlHVRVjPsFODcx8C94Rwb5cuXmpxLr/EqsUyuORSV0SmAf2BH0RkDNCEDFvl0q2Hl1qcyz/169t+TWomLQljhP71IvIC8IOqfhpMZNYdeCnp6Fx6eKnFufxTu7btS0tTcvswGvRvBUap6mQAVR0DjEn2vi5NvNTiXH5KcXIJo+rqGOBPi4GJyBYh3NulmpdanMtPmV5yAW4BTgrWptfgXG3gLODkEO7vUsVLLc7lr9jsG2VlKbl9GMnlb9iiXEdTnlwAWoVwb5dK997rpRbn8lUWlFwGAaNV9Zf4kyLy1xDu7VJl+XK4/347vu66aGNxzqVfpre5qOrwWGIJpoCpHTuf7L1dChUV2SSVPXrA3hXnHXXO5bzY9E6q639fgsIei/I9sE3I93RhW70a7r7bjq+6KtpYnHPRiI1viY13CVnYycVnOswGTz4J8+dDt25w8MFRR+Oci0IsuaRoMUAfRZ9vSkpg0CA7vuoqn/nYuXyV6clFRE6Le3k2MD9Y5viwZO/tUuCll+DHH2GbbeCYY6KOxjkXlRQnl7CXOZ6NTf3iyxxnIlW4/XY7vvxyqBNGZ0HnXFbK1OSCL3Ocfd56CyZNglat4LTTNvx+51zuSnGDfkqWOXYZ6rbbbP9//wcbbRRtLM65aK1aZftM7C0WW+a44nkR2TSZ+7oU+PRTGDsWNtkEzjkn6micc1FbHkwJ2aRJSm4fSm8xEWkWDKBsKyJtga1FZEgY93YhiY1rOffclP1ncs5lkWXLbN+0aUpuH8aU+3cCl2LzisX6tSrwdbL3diGZMQNeftka8C+4IOponHOZIMUllzC6C7UBWgI9gPGqOkdEDgDqhnBvF4Z//ctmPj35ZNhyy6ijcc5lgljJJYOrxb5R1UVYt+NTg3MfAveEcG+XrKVL4bGgWeySS6KNxTmXOWIll0ytFgPqisg0oD/wg4iMAZrgo/8zw2OPwa+/wn77QffuUUfjnMsUKS65JJ1cVPV6EXkBmKqqn4rIImwg5UtJR+eSs3atVYkBXHpptLE45zLL4sW23zQ1nXvDaNBvTHl12JXAb0Cxqs5O9t4uScOHw08/wXbb+QSVzrlyZWVQXGzHLVum5BFhVF09gY3Qbwagqp8D9YNeZC4qqjB4sB1fckn5kqbOOffLLzaJbdOmmTmIMrBWVfcBJsedmw2cEcK9XaI+/hgmTIDNNoNTT93w+51z+ePnn22/+eYpe0QYyeW7YK8AIlIfuAKYG8K9XaJiSxifcw40aBBtLM65zLJwoe1TmFzC6C02UkSeAzYVkR2AQ4GGwOEh3NslYu5cGzRZu7aNyHfOuXhpKLmE0VvsCxEpBPoBbYGrgTdVdXGy93YJKiqy+tS//hXatIk6GudcplmwwPaZnFwAVHUF8Hz8ORE5TVWfDOP+rgbWrIFHH7XjCy+MNhbnXGaaNcv2bdum7BFhrER5pojMEpHVIlIabGWAT1wZheHDrcjbtSvsu2/U0TjnMlEsubRrl7JHhFFyuQcoBMYDJcE5AU4I4d6uph54wPYXXAAi63+vcy4/pSG5hNFb7F1ghKr+qKqzgm0m8EwiNxORPiIyTkTaVzh/s4hosE2KO99ARO4XkVtF5GERaVKdazlp4kT45BPru37yyVFH45zLVJlYchGRinUtQ4GbROT1CuePAS6q4b03BxoBu1U43xBoDvQJTs2Ku/wQ8LaqDhORvsCjwInVuJZ7HnzQ9qefDo0aRRuLcy4z/fYbLFoE9eqltEFfVLVmHxD5BOgGLALKqnobsLmq1niAhYjUAkqBDkEJCBG5OLj8iKr+Hvfe1sAMoJmqrhSROsAyoAuwpqprsfvGFBQU6BGMekoAACAASURBVIQJE2oaamZZssSm0//9d5g6FTp2jDoi51wm+vZb6NIFtt0WfvghqVuJyBeqWlDZtUTaXP4B/KSqUzfw0ITGuahqmcS1FYi9OAHYFfiniJyrqs8Fl3sBi1R1ZfDZEhGZAfQE1q7n2sxEYstoTz5pieXAAz2xOOeqNjX41Z3i3xM1Ti6q+l413/dazcOp9D4K7CEizbAVL58RkaWq+hawJVBxPM0KoDXWuaCqa+soLi6moKA8+RYWFlJYWBhG+OmhamNbwAdNOufWb8oU22+3XUIfLyoqoij2+8aaKyoVyjiXdFDVX4Brg5LMxcBb2JQzqyq8tT5WamED1/7QokULsrpa7KOP7D9Mq1Zw6KFRR+Ocy2Tff2/7zp0T+nj8H9/BEiuVysapcv+FzQQANn9ZswrXGwPzNnAtt8QGTZ55JtT11aWdc+uRZMmlusIYRHlaJee2E5HDkr13FcqAicHxGKC1iDQInlsXSzwfbOBa7li8GF56yca0nH121NE45zKZannJJcXJJeFqMRHpgfUK6xM0lMfbFCthVOyeXJ37xlrzJXjdGmu4fx6rBrsUuAZAVReIyEigNzASOAAYrqpzg89WeS1nPPmkTfly8MEp7bPunMsBxcW2lkuTJlaNnkLJtLksxca4dMAWC4u3BnispjcMVrU8JXg5QEQewKq2bgcGAmOBeyt0JT4PGCQi3YBW2GwB1bmW/eIb8rOpA4JzLhqTgvHnO+6Y8hk8Ek4uqvq1iPQEjlPVGieSKu65Ang42GKKKW9jqewzvwCV1get71pO+PBDK+K2bu0N+c65DfvqK9v/5S8pf1RSbS6quryqxCIiA5K5t6uGWKnlzDOhTtZ0/HPORSVbkguAiJwtInNEZG2FWZEfDyE+V5WlS21BMBE4w1eUds5VQxqTSxh/7t4FnIPPipxezz9vI/J794b27aOOxjmX6Vatsm7ItWvb9C8pFkZyeQ8Yqaq/xp8UkadCuLeryhNP2H7AgCijcM5li6++grIySywNajztY42FkVyGAreKyEtx5wSbFdmXQkyF776DceOsO+HRR0cdjXMuG3zyie333DMtjwsjuQwEtgEOx8ahgLXltMKTS2r85z+2P/542HjjaGNxzmWHTz+1/V57peVxYSSXW4A3VXWdObuC9VNc2EpK4KmgxtGrxJxz1aEKH39sx9mSXFT11SouzarivEvGf/8LCxbY1A1pKt4657LcrFn2e2PTTaFTp7Q8MpGVKO8HnlPVT4LXD2OzDa/zNmw1ydR3Scg3sSqxAQNSPsLWOZcjYqWWPfdM2++NREous7GpX2I2AlYDCyhvcxGsHcaFadEieO01qFUL+vePOhrnXLZ4L1iGa9+Kq9SnTiKLhd1Z4dQtwGxVXRN/UkT+k0xgrhLPPQdr10LfvraksXPObYgqvPuuHffpk7bHhtHmMk1EGorI8cBWwPfAG6rqbS5hi1WJnX56tHE457LH1KkwZw60aAHduqXtsUknFxHZFXgD6348E6sSu0lEjlXVb5O9vwtMmgRffgnNmsHhh0cdjXMuW8RKLb17W5V6moTxpDuwqe1bqGqBqu6CTcHvi7mHaehQ2594Imy0UbSxOOeyRyy5HHBAWh8bxjiXz1R1ePwJVV0qIqtDuLcDKC21ucQATjll/e91zrmYVatg1Cg7PvDAtD46jJJLo2CRL8CWExaR/sAeIdzbAXzwAcybBx06wB7+bXXOVdO778LKlbDLLtCmTVofHUbJ5SHgUxFZCdQFtgZ+BQ4L4d4O4NlnbX/SST62xTlXfSNG2P6oo9L+6DB6i00RkV2AfkBHbBzMCFVdley9HTat/kvBnKAnnxxtLM657FFSYuPiAI48Mu2PD6vrQEtgE2AtsCLYuzC89RYsW2aL+2y/fdTROOeyxccfw+LF0LEj7LBD2h8fRlfk04BHgDXAdKAeMCjoijw52fvnvWeesb2XWpxzNfHcc7Y/+uhIqtPD6op8F9BSVburahegD3BBCPfOb8uWwciR9h/jBF/Y0zlXTatXw7BhdhxRD9MwksvPwB2q+kfXY1WdC3wXey0ivuhIIl5+2f6T9OwJW20VdTTOuWwxciQsXWrV6V27RhJCGMnlUuACEWkbt+0ObCsibUSkHXB5CM/JP7EqsZNOijYO51x2ia35FOEEt6KqG37X+m4g8g2wodYiVdXaST0ohQoKCnTChAlRh7Gu+fNtcso6deDnn23aF+ec25BFi6B1axt8PWcObLFFyh4lIl+oakFl18IoufwL2EJVa1W2AbWB/wvhOfll2DCbzfSQQzyxOOeqb8gQmz394INTmlg2JIxxLkUbuK4i8lCyz8k7L7xg+xNPjDYO51z2KC2Fhx+247/9LdJQ0jJFpqr6uJeamDMHPv3UJqjs1y/qaJxz2eLtt2HmTJsqqm/fSENJ3/zLrvqGB/OAHnIINGoUbSzOuezx4IO2P++8tE6vXxlPLpnoxRdtf8wx0cbhnMseU6ZYyaV+fTjjjKijSV1yEZEBqbp3Tps716ZtqF8fDj006micc9li0CDrBDRgAGy2WdTRJJ9cRORsEZkjImtFpDTYyoDHQ4gv/7z8su0PPhgaN17/e51zDuCnn2xsS61acHlmDCsMY8r9u4BzgPFASXBOAJ+vJBGxKrFjj402Dudc9hg82GZBPvFE2GabqKMBwkku7wEjVfXX+JMi8lQI984v8+fD2LFWJXaYL4fjnKuG4mL497/t+Moro40lThjJZShwq4i8FHdOgGOAC0O4f/4YPtzqTPv29Sox51z13HqrrTZ56KHQrVvU0fwhjOQyENgGOByIzSVTC2iFJ5eaiS0K5r3EnHPVMWsWPPSQzZx+yy1RR7OOMJLLLcCbFQdKiki0I3iyzcKF8OGHUK+eV4k556rnhhtgzRpb72mnnaKOZh1hTP/yqog0FJETgK2AKViyeTvp6PLJ669blVjv3tC0adTROOcy3eTJMHSoTW57441RR/MnYaxEuSvwBlYVNhNrb7k5WIny22TvnzdefdX2Eax17ZzLMqpw8cVQVgbnngtbbx11RH8SRrXYHcB5wMsazN8vIpsANwIX1fRmItIHuBk4XlVnBucaAIOAFUAz4EpVXZ7MtYzy22/w7rt27FVizrkNeeUVGDXKZkz/5z+jjqZSYYzQ/0xVh2vcwjCquhRYvZ7PVEpENgcaAbtVuPQQMFZVBwKvAo+GcC1zvPMO/P477LFHpFNkO+eywMqVcMkldnzLLdC8ebTxVCGM5NJIRP7oNysidUWkP7BHTW+kqj9jSeAPItIaOAl4PTg1CjhcRNoneq2mcaXciBG2P+KIaONwzmW+22+H2bNtCePCwqijqVIY1WIPAZ+KyEqgLrA18CuQUP2OqpaJSPypXsAiVV0ZXC8RkRlAT2BtgtdmJhJbSpSU2HrX4O0tzrn1mzzZkgvA/fdD7Yxd4DeU3mJTRGQXoB/QEZgNjFDVVcneO7AlsLjCuRVAa2y6mUSuraO4uJiCgvKVOgsLCylM118EY8fCkiXQqRN07pyeZzrnsk9pqc12vHatlVj22SeSMIqKiigq+mONyCrr5MIouaCqq4GX48+JyPaq+l0YtwcqJqr6WMmEJK79oUWLFkyYMCHJMBMUqxLzUotzbn3uvRc+/xy22spmQI5I/B/fIrKoqvfVOLmIyP3Ac6r6SfD6YeyX9jpvwxrlu9T0/pWYi/X0itcYmIclnkSuZQbV8i7I3t7inKvKDz/Atdfa8aOPZsVYuEQa9GcDS+NebwSsCc7PitsqVkklagzQOuhWjIjUBdoCHyRxLTN8+60tSdqiBey+e9TROOcy0dq1NgL/99+hf39boTYL1Ljkoqp3Vjj1MfBk/PQvIrIdCSYXKW/Nl+B5C0RkJNAbGAkcAAxX1bnB+xO6lhHefNP2hxyS0Q1zzrkIXX89jB8P7drBv/4VdTTVlnCbi4j0wBJAL+B7WbeLVzPgcuD+Gt6zMXBK8HKAiDygqsXYIM1BItINmxAzvrU90WvRi08uzjlX0ejR1jusVi14+mnYZJOoI6o2iRv7WLMPiuyITbffAfilwuU1wDOqmnkT3lSioKBA096gv2yZDX5StfUYmlVsHnLO5bXiYhvLMm+elV5uuCHqiP5ERL5Q1YLKriVcclHVr0WkJ3Ccqj6WcHT5atQoG+PSo4cnFufcukpLbVXJefNgr73KG/OzSFIj9FV1eWWJRUQ2Tea+ecGrxJxzVfnHP+C996BlSxg2zGY+zjJhTP+CiDQTkTYi0lZE2gJbi8iQMO6dk1Q9uTjnKjdiBNx2m3XyGTbMxrVkoTCm3L8TuBQbVxJr1Ffg62TvnbO++goWLIAtt4Qdd4w6Gudcppg8GU47zY5vuw169Yo0nGSEUdZqA7QEegDjVXWOiByAzTPmKhNfall3HjXnXL4qLoZDD4Xly+HYY+Gyy6KOKClhVIt9o6qLsNmHTw3OfQjcE8K9c5NXiTnn4v3+u00BNXMm7LorPPFE1v/hGUbJpa6ITAP6Az+IyBigCSG15+ScxYvhs8+gbl044ICoo3HORa2sDM46Cz75xNpXXn0VNt446qiSFsasyNeLyAvAD6r6qYgsBv4CvJR0dLno3XftP9N++0GjRlFH45yL2lVXwTPPQMOG8PrrObNgYCilC1WdrKprguPRqno3sH8Y9845seWMDzoo2jicc9EbPBjuvNO6Gg8fboMmc0TSyUVEzhaROSKyVkRKYxvweAjx5RbV8uTiVWLO5benny5vtH/iiZz7gzOMNpe7gHOA8dgCXWBdkk8I4d65Zdo0+Oknm/alW7eoo3HOReWVV2DAADsePNhmPc4xYSSX94CRqvpr/EkReSqEe+eWUaNs37u3TUTnnMs/b74Jxx9vU7xcfTX83/9FHVFKhJFchgK3ikh8A74AxwAXhnD/3BFLLl4l5lx+GjUKjj7a1mi55BK45ZaoI0qZMJLLQGAb4HBsZD5YW04rPLmUKy216bPBk4tz+WjUKDj8cFi9Gs47z6rDsnwsy/qEkVxuAd6MXywMQET6hnDv3PHFF7B0KWyzDbRvH3U0zrl0evNNK7GsXg1nngkPPJDTiQXC6YrcoWJiCZSIyI0i4uv3gleJOZevRoyw0ferV8Pf/gZFRXnR5hrGV3iuiIwVkS9E5AqAYN36EcCjwHYisl8Iz8lunlycyz9PPw3HHFPexvLAA3mRWCCc5NIG+AoYArQUkf7AlsDGwXr1T2HtMvlr1Sr4+GMrBu/vY0udywv33AP9+1t768CBOd/GUlEoDfqqel/shYicBNTDljpGVVVEOoXwnOw1bhysWWOjbzf1ddScy2mqlkxuv91e3323lVryTBjJpa2IdAd+B3oDuwJzgGVgC4kB9UN4Tvb68EPb77tvtHE451JrzRo4+2wYOtQW+/rPf6z0kofCSC7PYu0rbbBR+hdgVWSfici1QGfg3RCek708uTiX+375xXqEvf++zWr8wgvQr1/UUUUmjFmRvwDaiUgLVS0GEJFeqrpYRFpgyebhZJ+Ttdassam0AXr0iDYW51xqTJ9u6zN9/z20agUjR8Iuu0QdVaRC6bYQVH1tJCJtRaQt0EFEhqhqsaper6oLwnhOVpo40Rr0O3eGli2jjsY5F7ZRo2yBr++/t2XLx43L+8QC4cyKfCewGJgZt40Ddk723jnBq8Scy02qcNddNpvxkiVWchk7Ftq2jTqyjBBGm0sboCXQAxivqnNE5ACgbgj3zn6x5NKzZ7RxOOfC89tvNtJ+2DB7/Y9/wA035M0YluoI4zvxjaouAl4HTg3OfQjcE8K9s1tpKXz0kR17e4tzuWH6dNhzT0ssjRrByy/DjTd6YqkgjJJLXRGZBvQHfhCRMUATQmrPyWr/+x8sXw4dOkCbNlFH45xL1n//CyeeaD3DOnWyqV223z7qqDJSGL3FrheRF4AfVPVTEVkEdAde2sBHc5+3tziXG0pKbHr8G2+EsjI47DB46ilo2jTqyDJWGA36p6nqZFWNjcgfA7wB+FKLY8fa3qvEnMtes2ZBr17WpqIK119vJRZPLOuVcMlFRHpgi4L1EZEZFS43A+7H2mHy12ef2X7PPaONwzmXmGHD4JxzYNkyaN3aSis+P2C1JFMtthRbhbIDsHeFa2uAx5K4d/abOxfmzIEmTWyMi3Mue6xYARddBE88Ya+PPBIeeww22yzSsLJJwslFVb8WkZ7Acaqa34mkMuPG2X733b0XiXPZZPx4OOkkmDYNGjSw2Y0LC/NqRuMwJJRcgpmPwZY1HhOc2xEowuYSews4T1WXhRFkVopVie2xR7RxOOeqZ/Vqa7S/7TZrwN9pJ3juOdhhh6gjy0qJ/kn9NNAHGKuqPwZziL2HDaj8O/A1cG84IWYpTy7OZY9x42DnneGmmyyx/P3vds4TS8ISrRb7n6qeHvf6dmATYBdV/RpARP6TbHBZa+1amDDBjnf3VZ6dy1grV9ro+nvvtS7GHTvC4497D88QJFpy+Tx2ICK7AwOAh2KJJVA7ibiy29df22SVHTt6A6Bzmer9963q6+677fUVV8CkSZ5YQpJoyWUzEekKrMSqyGYD18YuisgWQP4uZBBrzPcqMecyzy+/wNVXw6OP2uuddrLSSkFBtHHlmERLLtdgI/CnASVAP1X9VUTaichlwGdYNVmoRORmEdFgmxScayAi94vIrSLysIg0iXt/lddSavx42++2W1oe55yrhrIy607cqZMllrp1bcT9+PGeWFIgoZKLqk4BOotIM1X9Je7ST9jgyfvDCC6eiDQEmmMdCQBmBfuHgLdVdZiI9AUeBU6sxrXUmTjR9r6mg3OZ4fPP4YILyv/w69ULHngAunSJNKxcJqoadQzVIiIXB4ePqOrvwbnWwAygmaquFJE6wDKgCzaQs9Jrqjoz/t4FBQU6IdYAn6zff4fGjW1G5BUroGHDcO7rnKu5hQth4ECr9gLYcksYPBiOO87HrYRARL5Q1UqLfVkxuk9EBDgBGAz8LCKx0kcvYJGqrgRQ1RIsofTcwLXU+eYb68rYubMnFueisnYt3H+/VYE9/rhVgV19NUyZAscf74klDcKYcj/l1IpXewTLKV8KPCMiS4EtsVUw460AWmNtQVVdS50vv7T9zr4Qp3NppwrDh1tp5Ycf7FzfvnDffZZoXNpkRXKJCdp3rg1KMhcD7wKrKrytPrA2OF7ftT8UFxdTENegV1hYSGFhYWJBxtpbundP7PPOucR88IF1J/48GCnRqRMMGgSHH+4llRAVFRVRVFQUe9m8qvdlVXKJ8y/gfWAuNgNzvMbAPGxqmqquraNFixaE1uYSSy5ecnEuPb75Bq66Ct54w15vvrlNj3/mmVYd5kIV/8d3sH5XpbKizaUSZcBEbF6z1iLSAEBE6gJtgQ82cC01Skps9UmAv/wlZY9xzgEzZ8IZZ0C3bpZYGjWyrsXTpsG553piiVhWJBcRaS0iJ4lIraBK7FLgGlVdAIwEegdvPQAYrqpz13ctZYFOmWK9xTp0gGYVC03OuVDMmAFnn20zYPznPzbr+AUXwI8/2lQujRpFHaEje6rFmmHzlw0ExgL3xnUnPg8YJCLdgFZAfGPJ+q6FL1Zq6eaLcDoXuunTbdbioUOtlqBWLejfH667DrbdNuroXAVZkVxUdTJWpVXZtV+As2t6LSUmT7Z9165pe6RzOe/HH8uTSmmpJZVTT4Vrr7XSi8tIWZFcskYsufioX+eSN2kS3HknPP+8JZXateG00+CaazypZAFPLmH65hvbe3JxLjGqMHq0dSF+5x07V7s2DBhgScWrv7KGJ5ewrFxpdcJ16sB220UdjXPZpaQEXnrJkkpsIHLDhnDWWXDJJdCuXbTxuRrz5BKWKVPsr66OHaFevaijcS47LF1qPb7+9S/rWgzQsiVcdBGcdx5summk4bnEeXIJi7e3OFd9kyfbrMRDh1qpH+wPs8sus8b6jTaKNj6XNE8uYfHk4tz6lZbC66/bhJKjR5ef790bLrwQDj3U2ldcTvDkEpYpU2y/ww7RxuFcppk/H554whbomhUsw9SwoZVQLrjAf2ZylCeXsEybZnvvIumclVLeeQf+/W947TV7Ddbb6/zz4fTToWnTaGN0KeXJJQxlZTbQC2CbbaKNxbkozZ0LQ4bYcsKzZ9u52rXhqKOgsBAOPNAGQbqc58klDHPn2pxiLVtCkyZRR+Nceq1aBa++ao3z//2v/bEFsPXW1pV4wADYYotIQ3Tp58klDLEqMR/g5fJFWRl89JEllJdeguXL7XzdunDMMTax5P77eyklj3lyCYMnF5cvvv8ennoKnn66vHEeYLfdrIH++OOheZXrR7k84sklDN6Y73LZokXwwgtWShk3rvx827Zwyik2M3HnztHF5zKSJ5cweMnF5ZrFi+GVVyypjB5d3turUSM49lgrpey7r1d7uSp5cgmD9xRzuWDxYhgxwhLKe++VJ5TateGggyyhHHkkbLxxtHG6rODJJQw//WT7tpUuOeNc5lqyZN2EUlJi52vXtm7Dxx1nCWWzzaKN02UdTy7JWrnSfkDr1oUWLaKOxrkNmz3bpmF59VUYM2bdhNKnT3lC8YZ5lwRPLsmaO9f2W23l9c8uM6nCV19ZMnnttfIp7cESygEHWEI56ihPKC40nlySFasS22qraONwLt6aNfDBB+UJJfb/FGxer4MOgsMPh379PKG4lPDkkqw5c2zvycVF7aef4O234a23YNQoWLGi/FqrVpZMjjjCBjf6lPYuxTy5JCuWXNq0iTYOl39Wr4axYy2ZvP12+bIPMV26WDI54ggoKPBqW5dWnlyS5dViLl1UbYT8e+/ZHF6jR8Nvv5Vfb9TI2k/69rXNlwZ2EfLkkqz5823funW0cbjcNHu2JZPRo22bN2/d6zvtZInk4INhr718iW2XMTy5JKu42PYtW0Ybh8sNCxeWJ5LRo8sH6Ma0bGltJrESypZbRhOncxvgySVZixbZ3nvcuEQsWAAffwwffmjJ5Jtv1r3etCn06mUJZf/9rR1FJJJQnasJTy7JipVcfACl2xBVmDrVGuFjW2xeupgGDaBHj/JksvPOvq68y0qeXJKxdi388ov9JdmsWdTRuEyzZo0NWIxPJrGSbkzDhrDnnrDPPrDffrD77lC/fjTxOhciTy7JWLLE9ptt5n9d5jtVmD4dPv/cpqX//HOYONG6C8dr1cpKJnvvbQmlWzeo4z+GLvf4/+pkxKrEvL0l/yxaBOPHlyeSzz+3WYUr6tzZksk++9jWoYO3mbi84MklGfElF5ebVK27+Zdf2jZxou1nzvzze1u2tGqt3XazfUGBV5e6vOXJJRmx6TUaN442DheOsjLr+htLJLFt4cI/v7dBA9hll3WTSdu2XipxLuDJJRm//mp7Ty7Zp7jYpkv55pvy/aRJ687HFdO0KXTvvu7WubO3lTi3Hv7TkQwvuWS+pUv/nEQmT668NAKwxRblCWTnnW3fvr2XSJyrIU8uyYgll0aNoo0j35WV2TQpU6fa9v33tn37bfl6OxU1amQDErt2Ld/vuKP15nLOJc2TSzK8Wiy9Fi0qTyCxJDJ1Kvzww5+7/MZstBHssMO6SaRLF28fcS7FPLkkw0su4Vqzxkog06fDjBm2xY6nTy/vnVeZVq1gu+2gU6fybYcdrOuvj0FyLu08uSTj999t36BBtHFki9WrrZrqp59si08eM2bY2jhlZVV/vlGjdRNI7LhjR2jSJH1fh3Nugzy5JGPtWtvXrRttHJlgzRpLHHPmlCePisdVNaLH1Kpl1VVbb20ljg4dyo+33ho239yrspzLEjmfXESkATAIWAE0A65U1eWh3LykxPa52iW1rMyqohYssO3nn8uPK24V58yqTO3aNkX8VlvZyp2xBBJLHm3a+HokzuWIHP2tuI6HgLdVdZiI9AUeBU4M5c6x5JINJRdVayNavNgSweLFf97izy9caMkk9jVuSK1atmBaLHG0afPn41atvP3DuTyR08lFRFoDJwHnB6dGAcNFpL2qzkz6AekouZSUwKpVsHJl+bZiBSxbBsuX2z7+uLL90qVWAolV49XEJptYUth8c9tXtTVvnrslOOdcjeX6b4NewCJVXQmgqiUiMgPoCcxM+u7VTS6//QannAIjRiT9yKQ0bGjzoFW2NW++7uuWLS2hbLRRtDE757JSrieXLYGKU9WuANZZ8L64uJiCgoI/XhcWFlJYWLjhu8dKAhtKLh9+mN7E8sILNmVJ06bWi6ppU9h0U08UzrmkFRUVUVRUFHtZ5ZTwuZ5cFFhV4Vx9YJ36oRYtWjBhwoSa3/3kk22KkK5d1/++3r3h8svh+eetwbpePWunqVfPFoZq2NC6M2+88bpbgwb2npr0kGrdGo49tuZfi3POVUP8H98iUmVPnlxPLnOxHmLxGgPzQrn7UUfZtiH16sGgQbY551weqBV1ACk2BmgddEdGROoCbYEPIo3KOedyXE4nF1VdAIwEegenDgCGq2oVsxk655wLQ65XiwGcBwwSkW5AK6AaLfXOOeeSkfPJRVV/Ac6OOg7nnMsnOV0t5pxzLhqeXJxzzoXOk4tzzrnQeXJxzjkXOk8uzjnnQufJJUlxc+y4NPLve/r59zz9svl77sklSdn8j5/N/Puefv49T79s/p57cnHOORc6UdWoY4iciBQDsxL8eHOgGmv8upD59z39/Huefpn+PW+nqi0qu+DJxTnnXOi8Wsw551zoPLk455wLnScX55xzofPkkiARaSAi94vIrSLysIg0iTqmfCEifURknIi0jzqWfCAiR4vIDBFZLCL3iUjOz6aeCUSkp4h8LyK/iMi9UcdTU55cEvcQMFZVBwKvAo9GHE9eEJHNgUbAblHHkg9EpC1wJHAMcAlwFnBxpEHlARFpCuwN7AGcApwvIgdEG1XNeG+xBIhIa2AG0ExV/DIObQAACIBJREFUVwZ/yS0DuqjqzEiDywMiUgsoBTr49zu1RKQH8KmqlgSvB2H/z/tFG1luC5Zm/12DX9Ai8gVwmaqOiTay6vOSS2J6AYtUdSVA8IM3A+gZZVD5QlXLoo4hX6jqR7HEEpgLzI4qnnyhqqviEktD4Dvg/UiDqiFPLonZElhc4dwKoHUEsTiXTrsCj0QdRL4ISo5vAXWBjSMOp0Y8uSRGgVUVztUH1kYQi3NpISLbAsWqOinqWPLINOAJoDdwZ7Sh1Iz3+kjMXKBZhXONgXkRxOJcyolIbeAc4KqoY8knqjofGCIiZcDlUcdTE15yScwYoHXQ6IaI1AXaAh9EGpVzqXMZcJeqro46kDw1AfujNmt4ckmAqi4ARmJFVYADgOGqmlX/+NlKRCR2GGkgeUJEBgJfAg1EZGsROSOoInMpIiIbi0j3uFOHAPdFFU8ivFoscecBg0SkG9AKKIw4nrwgIo2xfv8AA0TkAVUtjjKmXCYi1wI3VTg9RVWHRBFPHukKvCEiU4DPgM9V9Y2IY6oRH+finHMudF4t5pxzLnSeXJxzzoXOk4tzzrnQeXJxzjkXOk8uzjnnQufJxTnnXOg8uTiXABG5QkQejDoO5zKVJxeXUURkexF5WUQeE5FJIqIicsqGP5l2bwLPRR2Ec5nKk4vLNCOAf6rqWUB3YGj8RRG5MIqgRKSliBwfe62q36jq2ChiySVR/Xu61PPk4jKGiGwKdAJ+gz8WBbsOW+IAEdkJuDmCuOpjSa5Bup+dy0RkAHB01HG41PDk4jKGqi4BvgfeFpE9gnOzgM+DX/ADgCYicruI9BSR1iLyn6D94yMR2T92LxHpKyK3ish1IlIqIl+KyCnBpItTRGRfERktIr+IyKHBZ+qIyAMicpWIjBSR/wtudxiwA3CsiPyfiHQUkSdE5NG45+0jIveKyF0i8q6IdAzOV/m8eCKyq4gMF5HrRaRIRFaIyPMispGIdAq+zkdF5HER+TE4XyAiN4jIsKAqsWlwr7+KyKUi8oyIfBacayQiN4rI1SLyU/D92V1EporIDcF7jhaR1SLSq6bPXJ/gvveJyGAReTuIpQNwDLBN3L/nk8Hzzg9iHCMimwbf6yUi8njN/ke5SKmqb75lzAZ0AWYBZcAQoGXctfb2X/aP13cB1wXHlwOvB8fNsJVBmwev3wCeDY5bYCWhq7GSyNXAx8G1Q4HRwfGuwIq4Z70PDAiO6wOPAk8ErzcH/gfUCV5fBnwL1F7f8yp83XWwmbbHAB2BXbAS3FXARkBRcM8dgbOBpsDzcZ//GLgzOP4g7vwtwf584LjguAfQNzh+Ergh7v2zsGW8a/TMDfybLgLaBseTgb8GxwOA9yt8T8dipdcmwM/Aw8CmwBbAGqBT1P9Hfave5iUXl1FUdTL2y+wB4FRgcjDzdGWKgMdFZAusZNEoON8xOI6tFvod5VVai4L9c6q6ChiPzWoN8BFwmYhsDOwZd7+KMa4GFsSdOgWYquVrzT8GbA/ss4Hnxd+zJHjvB6r6g6p+ATwLHKqqvwP/397ZhcZVRHH89zdaS4tgKQZsEmsN2oCgCEoLRVNswIKIolTffFDyaAlUKfbFUqRYRFqwaKAvLQZRREGRooII+iBqDfj50Jg0xkK0YPqBCC2a48M5l72NN9mNLuwmnB8M3Jm9c8/M7HLPzDmzc6aBr8zsOzM7gh/BvlbSkKSh6GNBj6RhSZ3Ujmn/Az/F+wH8BT5aiJ7blGjPYmUuxA4zm5K0BVci/xrX0pj+ZGYnzewCcBL4wsxmzINm/Vo1dkl7kkfuJ22DPNphj5lNAjvDDHIcVyKbKqqcBnYDY/ix5Bui/HtgCo+3816Uvwm+7JEuCwPzF77CwMzOS9oE9AOf1mlu+aV8Mx7jnHjOOUkzQPdC8hrgB1zJFfLKMm8ATpnZoYp6D+GhcSfwFd2rwGvATcBbeOCpHQ3IX4zMhTgn6UDIPs38cXjmKrq/5+RnSVP+kiG/qKSduIrSS888VvseoG+e+4eBSTMbAS6W6v0J7AIGJD0OjJjZG/WES3oC2GxmB4HfF9HuKXzldNnjcP/R/2EFPnuvYhq4XxENFSAUI8AF4E7cBHdYUi9ulnou2nk1cCDunaVxZbeQzEokXYtPEPaZ2YkG5STLgFQuSbvxlKR1pXwX8EFcXwKQtEZSH75V+bowY22hFilxFb6r7HXga+BbSSuibvGbL8+ei+s7cLPPlUSUUUm9klaG7EJuUaeodwzoktQfdW4BxszsRB15VVxfut4KvFLKl5XAcWA18L6k+yQNUlPCT5vZrJm9jEeQBHhMUreZTQD7Sm34DbgrnPUP476OzlhFNixTvgmi6v9IvbgPbK2kjUAPsErSjdTGdGXky2NaUC+ftCutdvpkylQk3IlswFlgBDftHAPWxOcCPgK+wZ2+g8B54F3gEeAM8CA+M/8cmMHNUIavRG6POoZvce4GDuGrnu34bH8afyHfi5twng/ZT8bzB/GNBZ/hPodb4/PtIfMF3F/UHeXzyqvo/1HcH7ILj/64M8o34s7zX4D+0v39uAnwLPAScEWUz8a47QaGomwvMI6vBA8C66O8CzefjUWfR4H9uBlyMTJHgYmKPnUAH+JK7NkYmy+BTmAdMIn/t6kv+v4jcBvurzoDvA2sBx6NcRsGrmn1bzVT/ZSRKJNlR6wuBszscOQ78F1o28xNXm2JpKO4mW9vi5vyn5D0jJm92Op2JO1BmsWS5cgeYIOk1ZHvAAZw5347U2UWWhJIuofaDrQkSeWSLEv24zu4fpZ0CngH+MTMxlvbrPmRdDewGdgmae7mgKXAuJl93OpGJO1DmsWSJEmSppMrlyRJkqTppHJJkiRJmk4qlyRJkqTppHJJkiRJmk4qlyRJkqTppHJJkiRJms4/rKjqJm5g0vUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(6,6))\n", "plt.plot(vehicle1.stag_pres_atm[0:1700], vehicle1.q_stag_total[0:1700], 'r-', linewidth=2.0)\n", "plt.xlabel('Stagnation pressure, atm',fontsize=14)\n", "plt.ylabel('Stagnation point heat-rate, '+r'$W/cm^2$',fontsize=14)\n", "ax=plt.gca()\n", "ax.tick_params(direction='in')\n", "ax.yaxis.set_ticks_position('both')\n", "ax.xaxis.set_ticks_position('both')\n", "ax.tick_params(axis='x',labelsize=14)\n", "ax.tick_params(axis='y',labelsize=14)\n", "\n", "plt.savefig('../plots/decadal-2010-venus-flagship-thermal.png',bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-venus-flagship-thermal.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/decadal-2010-venus-flagship-thermal.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }