{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 31 - Atm. Reentry Demonstrator - 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 ESA Atmospheric Reentry Demonstrator.\n", "https://en.wikipedia.org/wiki/Atmospheric_Reentry_Demonstrator" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \n", "planet.h_skip = 120.0E3\n", "planet.loadAtmosphereModel('../atmdata/Earth/earth-gram-avg.dat', 0 , 1 , 2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('ARD', 2715.0 , 403.0, 0.00, 6.15, 0.0, 3.36, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,7.54,0.0,-2.6,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (30*60.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXxV9bX//9cKQxhljCBjZFJAATEgoBUsOHbCWai/XuuAU71WrdWq9Wq11aq3XztpjVZttbXeOlFHBKw4gDKDoFCQUVAMKIKMCVm/Pz7nJCcjCTknZ8j7+Xjsx9ln7332XtniysrnfPbnY+6OiIiIiIjUTVayAxARERERyQQqrEVERERE4kCFtYiIiIhIHKiwFhERERGJAxXWIiIiIiJxoMJaRERERCQOGic7gHjp2LGj5+bmVrm/oKCAnJyc+gsojele1YzuU83pXlVv3rx5m929Qd2g/eXsqujfUkW6JxXpnlSke1LRgd6T6nJ2xhTWubm5zJ07t8r9eXl51e6XUrpXNaP7VHO6V9Uzs7XJjqG+7S9nV0X/lirSPalI96Qi3ZOKDvSeVJez1RVERERERCQOVFiLiIiIiMRBgymsJ02alOwQ0obuVc3oPtWc7pXEi/4tVaR7UpHuSUW6JxUl4p6Yu8f9pNVe0OxE4E7gXHdfE9l2BvC/wEHAk8B17l4U2dcR+CWwFWgC3ODuheXPm5eX5+o7JCLpyMzmuXtesuOojHK2iEhZ1eXsem2xNrNOQCtgeMy2HsB44CzgGuBi4McxH/sn8Ki73wD8B/hFvQUsItKAKWeLiNROvRbW7r4JmFxuc0/gQnef5+5/Bf4InABgZiOAw9z9/cixLwJXmVnL+opZRGR/3J29e/cmO4y4S6Wc7Q5FRWEREUlV9d7H2t2Ly71/O/oVYsQGYF1k/ZvA6phjN0RWU/IrUxFpeAoLC8nKyqJly8z8ez9Vcva770KTJjB6dF3PJCKSOKn48OIw4E+R9a7AlnL7twNdyn+ooKCAvLy8kiU/Pz/BYYqIQOPGYTqAoqIiiouL93N0qfz8/JJ8BXRMTHT1ol5ytll4refHgkREgJrn7JSaIMbM+gAF7r4ossmBXeUOywYqPAiTk5Ojgc9FpN6ZGdnZ2ezZs4c9e/bQvHnzGn1u0qRJJU+km9nmRMaYKPWZs1VYi0gy1TRnp0yLtZk1Ai4FbozZvAFoF3OMAS2BjfUbnYhI1Zo2bQqQkf2sq1LfOVuFtYikg5QprIGfAPe5+56YbVOBw2Le5wJ7ATVNi0jKyM7OBmDPnj37OTKj1GvOVmEtIumg3gvrSAsGgMVsuwlYADQ3s15mdqGZ9XH3ucCnZnZE5NBTgN+7e8NpFhKRlJfJhXWq5GwV1iKSDuq1j7WZtQbOj7y9wMz+QPgq8Y5yhy5z90cj6+cCt5rZGsJ4qjfXR6wiIjWVqV1BUilnq7AWkXRQr4W1u28HHowsUXdGlqo+sxa4KMGhiYgcsExtsU6lnK3CWkTSQSr1sRYRSUuZWlinEhXWIpIOVFiLiNRRpnYFSSUqrEUkHaiwFhGpI7VYJ54KaxFJByqsRUTqSIV14qmwFpF0oMJaRKSO1BUk8VRYi0g6UGEtIlJHarFOPBXWIpIOVFiLiNSRCuvEU2EtIulAhbWISB2pK0jiqbAWkXSgwlpEpI7UYp14KqxFJB2osBYRqSMV1omnwlpE0oEKaxGROlJXkMRTYS0i6UCFtYhIHanFOvFUWItIOlBhLSJSR9HCWi3WiaPCWkTSgQprEZE6inYFUYt14qiwFpF0oMJaRKSO1BUk8VRYi0g6UGEtIlJH6gqSeCqsRSQdqLAWEakjdQVJPBXWIpIOVFiLiNSRuoIkngprEUkHKqxFROpI41gnngprEUkHKqxFROpILdaJp8JaRNKBCmsRkTpSYZ14KqxFJB00ru8LmtmJwJ3Aue6+JrKtOXAPsB1oB9zg7tv2t09EJBVkcleQVMnZKqxFJB3Ua4u1mXUCWgHDy+16AHjH3W8CJgMP1XCfiEjSZWqLdSrlbBXWIpIO6rWwdvdNhERbwsy6ABOBFyObpgHfNbPc6vbVS8AiIjWQqYV1KuVsFdYikg7qvSuIuxdbNEMGY4DN7r4zsr/IzFYDo4HCavatqc+4RUSqksldQVIlZ2dFmoH27avLWUREEqveC+tKdAW2lNu2HegCFFWzT0QkJWRqi3UVkpKzI7eYDPzbRUQySCoU1g7sKrctm9DywX72lSgoKCAvL6/k/aRJk5g0aVIcwxQRqdyBFNb5+fnk5+dH33aMf1QJk5ScHS2sG8bfLiKSamqas1OhsN5AeHI8VmtgIyGBV7WvjJycHObOnZuQAEVEqnMgXUFiC0kz25yQwBIjKTlbhbWIJFNNc3YqjGP9b6BLZIgmzKwJ0AOYsZ99IiIpoYF1BUlKzm7SJLwWFUFxcV3PJiKSGPVeWFvpUzAG4O6fAS8BYyPbxwHPuvuG6vbVY8giItXK5MI6VXK2mVqtRST11WtXEDNrDZwfeXuBmf3B3QuAy4F7zGww0BmI7WhX3T4RkaTL1FFBUi1nZ2eHonrPHmjePF5nFRGJn3otrN19O/BgZInd/iVwSRWfqXKfiEgqiG2xdnfKDU+XtlItZ6vFWkRSXSr0sRYRSWtZWVk0bhzaKQoLKwyAIXGiwlpEUp0KaxGROMjU7iCppFmz8KrCWkRSlQprEZE4yOQHGFNFq1bhdfv25MYhIlIVFdYiInGgwjrx2rYNr1u3JjcOEZGqqLAWEYkDdQVJPBXWIpLqVFiLiMSBWqwTT4W1iKQ6FdYiInHQJDI1oEYFSRwV1iKS6lRYi4jEQXTsandPciSZK1pYf/llcuMQEamKCmsRkThQYZ14OTnhddOm5MYhIlIVFdYiInGgwjrxunYNrxs2JDcOEZGqqLAWEYkDFdaJFy2sP/kkuXGIiFRFhbWISByosE68bt3Cq1qsRSRVqbAWEYkDFdaJd/DB0LgxbNkCu3cnOxoRkYpUWIuIxIEK68TLyoIePcL6qlXJjUVEpDIqrEVE4kCFdf3o3z+8fvRRcuMQEamMCmsRkThQYV0/ooX1hx8mNw4RkcrEpbA2s9PicR4RkXSVTgV1OufsAQPCqwprEUlFjWtzsJmdAtwCdKa0KM8COgHN4xuaiEj62LNnDwDZ2dlJjqRUJubsQYPC69y5yY1DRKQytSqsgSeAu4AFQHFkmwET4xmUiEi6ScXCmgzM2YMGQYsWsHIlFBSUzsYoIpIKaltYr3L335TfaGbL4hSPiEhaStHCOuNydpMmMGwYzJgBs2bBd7+b7IhERErVto/15WZ2ViXbldpEpEFL0cI6I3P2sceG1zfeSG4cIiLl1bawPhZ4xMz2xSzFwIMJiE1EJG2kaGGdkTn71FPD60svQRo9MyoiDUBtu4L8ApgEzAX2RbYZ8P/VNRAz6w/8CFgJ9AXy3X2hmTUH7gG2A+2AG9x9W12vJyISTylaWGdkzh45Ejp0gI8/hmXLSofgExFJttoW1q8Ck919T+xGM3s4DrE8AXzP3TeYWQ9gCtAfeAB4zd2fjjzh/hAwIQ7XExGJi71797J7924aNWpEixYtkh1OrIzM2Y0awXe+A48/Dv/4B9x+ezzPLiJy4GrbFeQfwM/M7PiYZTTwszjE0h9oHVnfBbQxsy6Ep9dfjGyfBnzXzHLjcD0Rkbj48ssvAWjXrl3JRDEpImNz9g9+EF4fewz27av+WBGR+lLbFuufE77yuwCI9mwzwhipV9cxlmcIfQFPBc4HrgLGAJvdfSeAuxeZ2WpgNLCmjtcTEYmLL774AoD27dsnOZIKMjZnjx4NvXrBqlUwZQqclrZT3ohIJqlti/W9QHd3z3X3QyNLLnBGHGK5AtgNzAG+dvdnga7AlnLHbQe6lP9wQUEBeXl5JUt+fn4cQhIR2b8DKazz8/NL8hXQMTGRZW7OzsqCyy4L63fdpYcYRSSxapqza9tivd7dt1eyvXMtz1OZ5sBThMR8f6SVwwlfMcbKBgrLfzgnJ4e5mopLRJJgy5ZQS7Zr167Gn5k0aRKTJk0CwMw2JySwDM/Zl10Wiup33gnjWo8ZU6fTiYhUqaY5u7Yt1leU32BmZwP31/I8lXkSeNbdfwHcTUjYnxKeKo/VGtgYh+uJiMTFhg0bAOjatWuSI6kgo3N269ZwzTVh/dpr1ddaRJKvtoV1ezP7OYCZtTCzR4BHCF8FHjAz6wgMdvetkU13AtuAtUCXyPBNmFkToAcwoy7XExGJpzVr1gCQm5ub1DgqkfE5+9proXt3WLAAHnooUVcREamZWhXW7v4tYKqZ/R5YCPQDBgOVzexVG18Au80strlnC7AYeAkYG9k2jtBCsqGO1xMRiZu1a9cC0LNnzyRHUlZDyNktW8L9kfb3G26AlSsTdSURkf2rto+1mR1fxa61wArgPqAnMB645kCDcPdiMxsP3Gpm84BOwPXuvs3MLgfuMbPBhH6Bkw70OiIiibB69Wog+S3WDTVnn346nHMO/N//wcSJoc9106aJvqqISEXm1TxKbWbvEcYq3QoUV3UYcIi7J3W6sby8PNfDiyJS34qLiznooIPYsWMHn3/+OTk5ObU+h5nNc/e8usbSkHP21q0weDCsWwcXXQQPPwypNaS4iGSK6nL2/kYF+QUw190/388F6vq1oohIWvr444/ZsWMHXbt2PaCiOs4abM5u2xaeeSaMb/3nP0OfPnDjjcmOSkQammr7WLv7K/tL0JHjnolfSCIi6WPBggUAHHXUUUmORDl72DD4299CS/XPfgZ/+EOyIxKRhqa2o4KIiEiMWbNmATB06NAkRyIQ+ltHC+qrroLf/S658YhIw6LCWkSkDqZNmwbAN7/5zSRHIlFXXAF//GNYv/rq0HpdXFWPcxGROFJhLSJygDZt2sSSJUto0aIFI0aMSHY4EuOKK0Jf60aN4O674dxz4euvkx2ViGS6OhfWZnaUmT1pZteZWW7dQxIRSQ+TJ08GYPTo0WRnJ3WQjRprSDn7wgvh1VfhoIPCg415ebBoUbKjEpFMFo8W68uAce7+v8Co6IxbIiKZ7qmnngLgnHPOSXIktdKgcvaJJ8J778HAgbB8ORxzTJhQRtOfi0gixKOw/hFwdGT9H0CTOJxTRCSlrV+/nhkzZpCdnc3pp5+e7HBqo8Hl7P79YfZsuOQS2LMHrrkGjj0WlixJdmQikmlqXVib2alm9v3I+nCgb3S6WncvdvdtcY5RRCTl/PGPf8TdGT9+PG3atEl2OFVSzg5atID8fHjhBejSBd5/H4YOhZtugu3bkx2diGSKWhXWZvZL4J/A6QDuPhv4gZlNTEBsIiIp6euvv+ahhx4C4JprDnhm8IRTzq7oe9+DDz+Eyy+HwkK46y7o1w8ee0wjh4hI3dW2xXoU0A14O2bbq8DdcYtIRCTF/e53v2Pr1q2MGjWKY445JtnhVEc5uxJt2sADD8CsWaHP9WefhQcd8/Lg9dfBPdkRiki6qm1hPdXdtwKxaSfjpsYVEalKQUEBd98d6tI77rgjydHsl3J2NUaMgJkzw2yN3bvDggVw8slw3HEwdaoKbBGpvdoW1qvM7CbgMDM7z8yeBa4Efh3/0EREUs/NN9/M9u3bOfXUU9NhUhjl7P3IyoKJE2HZMvj1r6Fjx1Bsn3RSKLBffFFdRESk5mpVWLv7P4A3gEbAD4CtwEnu/scExCYiklKmTZvGww8/TNOmTbnvvvuSHc5+KWfXXIsW8NOfwurVod91hw6hwP7ud+Hww0PXkR07kh2liKS6Wo8K4u7vuftl7n6au18E/NvMxiYgNhGRlLFlyxYuuugiAG677TYGDBiQ5IhqRjm7dlq1ghtvDAX2ffdBjx6wYgVceWXoLvKTn4TWbRGRyphX04nMzP4FtNrPOToCm909qd+J5uXl+dy5c5MZgohkqH379nHaaafx+uuvM3z4cN59910aN24ct/Ob2Tx3z4vDeZSz46yoCJ5/Hn7zmzDRTNQ3vgEXXwxnnRVau0Wk4aguZ+/vN8PiyDHLgWLgEmABEJsNDwXaxyFOEZGU9LOf/YzXX3+dnJwcnnnmmbgW1XGmnB1njRvD2WeHZfZsePhheOopePvtsPz3f8M558CECXD88dCoUbIjFpFk2t9vh3sJrdpbAcyst7vfWv4gM5uciOBERJLt3nvv5d5776Vx48b84x//oHv37skOqTrK2Qk0fHhYfvMbePrpUGRHi+2HH4ZDDoFzzw1F9rBhYJbsiEWkvlXbx9rdv4om6IhOZlbm73EzGwMMT0BsIiJJ9cADD/DTn/4UgMcffzzlRwFRzq4frVuHbiDvvx+mRb/pJjj0UPj0U7j//jA2dp8+cN118NZboTuJiDQMtX148RlgkZn9yczuNbOXgKnA7+IfmohIcrg7v/rVr7jyyisB+P3vf8/3v//9JEd1QJSzE2zgQPjlL+Hjj0Mf7Kuvhs6dYdWq0LI9enR4f8EFoa+2RhYRyWzVPrxY6QfMugBXAIOAbcA/3T3pXyumy4MwIpLaioqKuO666/jd736HmfHggw9y6aWXJvSa8Xp4sYpzK2fXs337wlB9kyeHZeXK0n3NmoW+2CefHJYBA9RlRCTdVJeza11YV3LyLOAod59XpxPVUSYnaRGpH5s3b+bcc8/ljTfeoEmTJjzxxBOce+65Cb9uIgvrSq6lnF2P3OGjj0qL7PffL7u/W7cwGc3JJ8O4cdBej5WKpLy6jApS/kSPVrI5B9gHjD+A2Cq7xgjgOGAl8A6wA7gH2A60A25w923xuJaISNS7777L97//fdauXUunTp145plnOO6445IdVp0oZyefWWiVHjAAfvYz2LQpTJc+ZQq8/jp88gk8+mhYsrJg6FA44QQYMybM/HjQQcn+CUSkNmo7ZlRfYDoQ28zdDXg7HsGY2Q+B3u5+S8y2x4DX3P1pMzsFeAiYEI/riYgUFhZy++23c9ddd1FcXMzw4cN57rnn6Nq1a7JDiwfl7BTTqROcf35Yioth8eJQZE+ZAu+8A3PnhuXee0OhffTRochWoS2SHmrVFcTMDnP35eW2HQ5c4O431ikQs9HAz4ETPRJUpG/gaqCdu+80s8bAV8BAd18T+/mG8rWiiMTP3LlzufTSS5k/fz5mxg033MDtt99O06ZN6zWORHUFUc5OLzt2hL7Zb74Zltmzy44oEm3RPvZYGDUqLN26JStakYYr0X2sDwXmuHvHOp5nHjATMKAPcAfQE7jX3bvGHLcksu0vsZ9XkhaRmvrqq6+4+eabeeCBB3B3evTowRNPPMHxxx+flHjquY+1cnaa2F+hDWGa9WiRPXIkDBkCTZokI1qRhiOefaxXU/YrxUZAJ+CJAw8vtKoARwE/dPfFZnY98BpwJ7Cl3OHbgS7lz1FQUEBeXunPOGnSJCZNmlSXsEQkwxQVFfHoo49y6623smnTJho1asSPf/xjbrvtNlq12t9M4PGVn59Pfn5+9G2dityqKGent5Yt4cQTwwKh0J41KywzZ4bX9evDZDVPPx2Oad48TE4zcmSYzGbYsNCqrZFHROqmpjm7tl1B8oG/U5qoi4FN7v6fA4wzet4zgPxoC4qZNQcKgF8AZ7r7MTHHzgf+7u73xZ5DrR8iUhV351//+hc33ngjy5YtA2DkyJH86U9/YtCgQUmOLqFdQZSzM1hxcRhxJFpkz5wJy5dXPK5TJ8jLC0V2dMnJqf94RTJF3FqsgSfd/a1yJz/MzPq4+8qqPlQDjWNjcfddZraC0LrSrtyxrYGNdbiWiDQQ7s5rr73GnXfeycyZMwHo1asXv/rVrzj77LPJyqrtHFlpRzk7g2VlhQlqBg6ESy4J2zZvDhPVvPcezJkTlk2b4OWXwxLVo0fZQvvoo6FNm+T8HCKZpLaF9SnAW+W2bQL+Tfha8EAtBtqYWSd33xTZVkRIxl3MrHkkcTcBegAz6nAtEclwxcXFTJ48mTvvvJP58+cD0LFjR2699VYuvfTSen84MYmUsxuYjh3h298OC4RxtFetCgX23Lnhdd48WLcuLM8+W/rZXr3gqKPCMmRIeD3kEHUjEamNGhXWZvYjQt+51mZ2QyWHTK1LEO6+zMxeBc4C/mhmbYGmwN+AU4GxwEvAOOBZd99Ql+uJSGbauXMnf/vb37j//vv58MMPAejUqRM/+clPuOyyy+q9H3WyKGdLlBn07h2W884L2/btg2XLyhbbixaFAnzVqrLFdk5O2UJ7yBDo2xcaNUrOzyOS6mrcx9rMBgH/Bfyu3K697v5pnQMx6wj8FphLaOF42N0/NLN2hMkG1gCdgZ+5+9flP6/+eiIN1/r163nggQfIz8/niy++AKBbt27ccMMNXHTRRTRv3jzJEVYvEX2slbOlNgoLQ7G9cCEsWBCWhQth69aKx7ZoAYMGhUJ70CA48kg44gh1JZGGI27D7ZlZE3cvrGR7D3dfV4cY60xJWqRhKS4uZvr06Tz88MM899xz7Nu3D4Bhw4Zx9dVXc/bZZ6dNl48EPryonC0HzB3Wrq1YbK9fX/nx3buXFtnR1/79ITu7fuMWSbQDfnjRzH4CvOLuH0Y2/dTMyn8BZMAJwJi6Bioisj/r1q3jscce47HHHmPt2rUANG7cmPPOO4+rr76aY445BmugnUKVsyWezCA3Nyzjx5du37w5dB1ZsAA++CAsH34YCu716+GVV0qPbdQodB0pX3D36qXuJJKZ9tfH+lCgbcz7oUBL4DNKh28ywtd9IiIJsWPHDv71r3/x+OOPM3XqVKLftOXm5nLhhRdy4YUXZsoU5HWlnC0J17EjjB0blqh9++Djj0ORvWRJ6euKFaGLybJl8M9/lh7frBkcdlho0T788PDav38owps1q/+fSSReqi2s3f3Kcpt+CSxx972xG81sQLwDE5GGbe/evbz++uv8/e9/Z/LkyezcuROA7OxszjjjDC666CJOOOGEhjBkXo0pZ0uyNGoE/fqF5cwzS7fv2hWK6vIF9yefhFbvRYvKnicrCw49tLTQji2627ZFJOXVdri9I919fuyGyAxcvYEPK/+IiEjNFBUV8c477/DUU0/xzDPPlDyICGFCl4kTJzJx4kTat2+fxCjTinK2JFXz5qVD+MX66qtQcH/0UenrRx+FVu/o8tJLZT/TqVPZgvuww8LSvbu6lUjqqOlwe98gfH14YmSK3FjtCU+dvxjn2ESkAdi9ezfTpk3j+eef51//+hebN28u2XfkkUcyYcIEzjvvPA499NAkRplelLMl1bVpA8ccE5ZYe/bAypWlhXZs8b1pU1jefLPsZ7KzQxeSaKHdr1/pervy0xWJJFhNW6y3An8l9N87tty+vcAj8QxKRDLbtm3bePnll3n++ed59dVX+frr0tHY+vbty1lnncWECRM48sgjkxhlWlPOlrSUnV06m2Ss4uLwYGRswf2f/4Qp3D/9NHQvWbKk4vk6diwtsmOL7t69IU0GDZI0U6PC2t0/MLMxwNnuroQsIrXi7qxYsYJXXnmFV155hRkzZrB3b2m336FDh3L66adz+umnM2DAgAY7qke8KGdLpsnKgp49w3LKKWX3bdsWiuxoob18een7zZvD8u67Fc936KGVF92abVLqosZ9rN39K6po5TCz49z9nbhFJSJpb/fu3bz55pslxfTHH39css/MOP744zn99NMZP348ubm5yQs0QylnS0Nx0EGQlxeWWO6wYUPZgjtadK9ZU9qXO3Z4QIBWrcp2J4kW3f36hX0i1dnfONb/Avb3z8iAgcDB8QpKRNKPu7Ny5UqmTp3Kq6++yvTp09m1a1fJ/vbt23PKKadw2mmncfLJJ9OxY8ckRpuZlLNFSplBt25h+eY3y+6L9uUu38q9fDls2QLz54elvK5dyxbdfftCnz6h9VtdSwT232K9nNBXbwOlY6CWZ8C34xmUiKSHTZs2MX36dKZNm8a0adNYX25KtqFDh3Laaadx2mmnMXz4cBrp0f1EU84WqYGq+nJDKKxjC+3o+ooVoQV8wwb497/LfibaVaVPn9JiO7r06qXZJxuS/RXWvwa2u/ue6g4yMz1dLtIAbN++nbfeeotp06Yxffp0PvjggzL7O3bsyNixYznxxBM59dRT6dKlS5IibbCUs0XqqEMHGDUqLLH27QtTvMcW3CtXhmXtWli9OixTp5b9nBn06FF50d27tybEyTT7myBmc3X7Y5yOnjIXyTjbtm3j3Xff5a233mLGjBnMmTOHoqKikv0tWrTg+OOPZ+zYsYwbN45BgwZpwpYkUs4WSZxGjULrc69eFR+g3Ls3FNXRQnvFitL1NWtC4b12LUyfXvZz0e4qVRXdLVrU248ncVLbCWLKMLMewJXAj1CSFkl7W7Zs4e233y4ppBcuXEhxcXHJ/kaNGjFixAjGjRvHuHHjGDFiBNn6jjNtKGeLJEbTpqX9rsvbuzcU1eUL7pUrQzG+fn1YyncvgdCnu3zR3bdvKLpbtkz8zyW1V+vC2swaAd8FLgFOBDYR+vSJSJrZuHFjSSH91ltvsaTcQLCNGzfmmGOO4fjjj2f06NGMGjWKNm3aJClaORDK2SLJ1bRpKIb79oVTTy27r7AQ1q2rWHSvWBGK7mif7hkzKp73kEOqLrpbt66fn00qqnFhbWaHAhcDPwQ6AGuB7wBTgJEJiU5E4qawsJDFixczc+ZMZs6cyaxZs1i7dm2ZY7KzsxkxYgSjR4/m+OOPZ8SIEbRUs0haUs4WSX1NmoRCuHdvOPnksvuKikJLdvmCe+VKWLUqTIzz6afw9tsVz9upU+VFd58+YXhCSZz9FtZmdjahpeObwErgN8BfgJvc/bXIYTMTFqGIHJAtW7Ywa9askiJ69uzZ7Ny5s8wxrVu3LimkR48ezbBhw9S1I80pZ4tkhsaNwzB+hx4KJ51Udt++faHoju1WEi26P/64dPr38hPjAOTkVF10t21bPz9bJqtJi/UuwhS4M4GL3f0/AGZW1VBOIlLP9u3bx0cffVSmkF6+fHmF4/r27cvIkSMZNWoUo0aNYsCAARoCL/MoZ4tkuEaNIDc3LOPGld1XXAyffFL5g5QrV0JBQVhmzap43g4dqi6627evj62eY30AACAASURBVJ8s/e23sHb3l4CXIg+9XGJmucBkoEn0GDMzd1fSFqkH7s7atWuZM2cOs2fPZvbs2cybN48dO3aUOa5Zs2YMHz6cUaNGMXLkSEaOHElOTk6Sopb6opwt0rBlZYXh/Xr0qDgxTnExbNxYddG9ZUtY3n+/4nnbtau66O7QQdPAR9VmSvN1wM/NrDEwHuhmZncAfwZOAvITE6JIw7Z58+aSIjr6WlBQUOG43Nxchg8fzrHHHsvIkSMZPHgwTTUVWIOlnC0i5WVllc5GOWZM2X3uoc92ZQ9SrlwJX34Jc+aEpbw2bSoOFxh9n5PTsIpuq0ujhZn1Aa4CLnH3pI62mJeX53Pnzk1mCCJ1tmPHDhYsWFDSEj179mxWr15d4bgOHTowfPjwkmXYsGFqjU5jZjbP3fPq4TrK2SJSa+6hz3ZVRff27VV/tk0b6N8/LAMGlK7n5oZCPx1Vl7PrNI61u68ErjazeXU5j0hDtG3bNhYuXMi8efOYP38+8+bNY9myZZT/Y7dFixYMHTq0TCGdm5uLNaQmAIkL5WwRORBm0LlzWI47ruw+99Bnu7KCe8UK+OoreO+9sMRq1gwOP7y00I4W3n37htFS0lWdCusod/9rPM5jZk2BOcDV7v6mmTUH7gG2A+2AG9x9WzyuJVKfvvzySxYsWFCmiF6xYkWF4xo3bszAgQPLFNEDBgygceO4/K8qAsQvZ4PytkhDZwYHHxyWY48tuy/a0v3RR2H58MPS9U8/hYULwxKraVM44gg46igYMiS8Dh4MrVrV389UF6n22/p6IDfm/QPAa+7+tJmdAjwETEhGYCI1tXnzZubPn19SQM+fP59Vq1ZVOK5p06YceeSRDB06lKOPPpqhQ4dy5JFH0qxZsyRELXLAlLdFpFKxLd0nnFB239atsGxZxaJ71SqYPz8ssefp0weGDoURI2DUqFB0p+JjRClTWJvZKOBT4MvI+y7ARML0uwDTgGfNLNfd1yQlSJEY7s4nn3zCokWLWLBgQUkxvW7dugrHNmvWjMGDB5cpogcOHKiHCyWtKW+LyIFq2zYUySNGlN2+bRssWhRashcsCMvSpaFbyYoV8PTT4bhmzeDoo2HkyFBojxoVJsZJtpQorM2sJXCWu19rZrdGNo8BNrv7TgB3LzKz1cBoYE1SApUGq7CwkGXLlrFw4cIyyxdffFHh2BYtWnDUUUcxdOjQkkK6f//+6s4hGUV5W0QS4aCD4BvfCEvU3r2hRXvu3DD+9syZobX73XfLToIzcGAY13vsWBg9OjmzTKbKb/obgLvLbesKbCm3bTvQpV4ikgbrq6++YvHixWUK6CVLlrB3794Kx7Zv354hQ4aUaY3u16+fJl2RhkB5W0TqRdOmoevHkCFw8cVh2xdfhAcio4X2rFmhZXvpUvjtb8MkOsOGwWmnwRlnhAcj6+OZ/6QX1mZ2MjDX3T8vt8sJM4jFygYKKztPQUEBeXmlI59MmjSJSZMmxTNUyTDuzvr160uK50WLFrFw4cJK+0MD9O7dmyFDhpRZunbtqtE55IDk5+eTn18ylHTHZMZSW/HI28rZIlIX7duHovm008L7PXtCoT19eljef790NJJbbw2jjZx1VijMe/Wq/fVqmrPrNI51PJjZVOComE3tgK+B3wDfd/d+MceuAP7H3f9e/jwaE1WqU1hYyEcffVShK8eXX35Z4djs7GyOOOKIMgX0oEGDOCgZ3ylJg1Bf41jHSzzytnK2iCTS9u3w5pvwwgsweXKYURJCq/W3vgW//CUMGnRg564uZ6dCYZ1DaNGImgVcC7wPfAjkuPsuM2tCSNy93H1D+fMoSUvUl19+yaJFi0paoBctWsTSpUsr7crRoUMHjjrqqJICevDgwRx22GE0SedBNCXtpGFhXee8rZwtIvWlqAjefhsefzw8/LhnT+gqcs89cO21tT9fwiaIiQd3LzM3s5ntAwrcfZ2ZvQSMBV4CxgHPVlZUS8Pk7qxevbqkgI4W0WvXrq30+D59+lToytGlSxd15RCpJeVtEUknjRuH4f5OOAHuuw/uuAP+8Ae47rowmkj5kUnqdK34nSohLgfuMbPBQGdAHfAaqN27d7N06dIyBfSiRYvYtq3ivBPNmjVj0KBBDB48uKSAPvLII2ndunUSIhdpcJS3RSRl5eTA/feHPtizZ4eRRjK6sHb33Jj1L4FLkheNJENBQUGZhwkXLlzIsmXL2LdvX4VjO3fuXKaAHjx4MH379tXQdiL1SHlbRNJBcTG89hr84hehqM7OhlNOie81VH1I0uzbt4+VK1dW6MqxcePGCsdmZWXRv3//kuI5+tq5c+ckRC4iIiLpoKgoFNHPPgvPPAPROdwOOQSeeCLM6BhPKqylXnz99dd88MEHZYroDz74gJ07d1Y4tlWrVgwePLhMS/TAgQNp0aJFEiIXERGRdLFrV5i5ceZMeOMNeOutMEJIVI8ecNVVMGlSYiaQUWEtceXubNy4sUIr9IoVK6hsBJru3btX6MrRq1cvsrKykhC9iIiIpItt28KMjIsXh77Sc+bAkiWhlTpWv35w6qlw9tlhCvRElhgqrOWAFRYWsnz58goTrGzevLnCsY0bN2bAgAFlCujBgwfToUOHJEQuIiIi6eKLL2DFijCN+dKloXheurS0W0esrKwwtfnw4aUjgXTrVn+xqrCWGtm1axeLFy9m7ty5zJ8/v9ppvtu2bVumgB4yZAj9+/cnOzu7kjOLiIhIQ+YeiueVK0MBXf61krncgDDVef/+oZAeOjRMYT50KLRqVb/xx1JhLRXs3buXDz74gLlz55YsS5Ysoaj8dytAr169KnTl6NGjh8aGFhERESAUzgUFsHYtrFkTXsuvVzJ6bomWLcOU5P36hSJ64EA44gjo3TuMUZ1KUiwcqW+FhYV8+OGHZYroxYsXV2iJzsrKYuDAgeTl5XH00UeXTPPdpk2bJEUuIiIiqWD7dvjkE9iwofR1/fqyhfOuXdWfo1WrMEJH374VXzt1ClORpwMV1g2Iu7Ny5Uree+89Zs+ezdy5c1m4cCG7d++ucOxhhx1GXl5eyTJkyBBaJfO7FREREalXxcWweXPZgrmy19hRN6rSrh307Fm65OaWXW/fPn2K5+qosM5gO3bsYM6cOcycOZNZs2bx3nvvVfpgYe/evcsU0UOHDuWgRIxBIyIiIkm3b18omD/9tOzy2Wel6xs3hsK5sHD/52vePDwg2LVr2dfY4rmhlBUqrDOEu7Nq1SpmzZpVsixevLjCbIUHH3wwI0eO5JhjjmHYsGEMHTqU9u3bJylqERERiZc9e8oWx+WL5ejy+eehuK6Jdu0qL5qjr926Qdu2mdHaHA8qrNNUcXExS5cuZcaMGbz55pu8/fbbfP7552WOadSoEUOHDmXkyJEly6GHHqoHC0VERNKEO2zdGgrk6FJZwfzZZ2FkjZpq3z7MPlh+6dw5vHbpEopnzc1WOyqs08S+fftYvHgxM2bMYMaMGbz11lt8Ue7/oJycnDJFdF5eHi1btkxSxCIiIlKVvXth06bSoji2YC6/vmdPzc7ZqFF40K+qYjm6dOoEGgE3MVRYpyh3Z8WKFUyZMoWpU6fy1ltv8dVXX5U5plu3bowePZoxY8YwevRo+vTpo9ZoERGRJHEPYy7vr1D+9NPatS63bh2K4+gSLZbLF8wdOyZ2VkHZPxXWKeSrr77ijTfeYMqUKUyZMoU1a9aU2d+zZ8+SInr06NHq1iEiIlIP9uypeetyJfOmVSorq7R1uaqCuXPncIwG5UofKqyTbNmyZbzwwgu8/PLLzJo1q8zDhu3bt+fEE0/kpJNOYuzYsfTs2TOJkYqIiGSO6Gx/NWldrmrmv8ocdFDFIrmyorlDh9B1QzKLCut6VlxczJw5c3j++ed54YUXWL58ecm+Ro0acdxxx3HyySdz0kkncfTRR9NI/9eJiIjU2O7dlT/oV379s89qNpQchNn9OnWqvmU5uuhhv4ZNhXU9cHfef/99nnzySZ577jk+/fTTkn3t2rXjO9/5Dt/73vcYO3asZjIUEREpp7i4tHV5f90xtm6t+XnbtKm+UI6ud+igvstSMyqsE2jFihX87W9/48knn+Tjjz8u2d69e3fGjx/P6aefzje+8Q0ap9pE9yIiIvUg+rDf+vWwbl2YxS92KLlo0bxpU+1al2vSsty5c5jYRCSeVNHF2Z49e3jmmWd44IEHmDlzZsn2zp07M3HiRCZMmMDRRx+thw5FRCTj7dwZiuZo4Vx+fd26cExNtG27/5blzp3D+MxqXZZkUWEdJ2vWrOGhhx7ikUceKZk2vFWrVpx55pmcf/75nHDCCeovLSIiGWXXLli9OiyrVpUu0cJ5y5b9n6NVK+jRIyzduoWJScoXzZ06QbNmif95ROpKhXUdLV26lLvuuounnnqK4uJiAAYPHswVV1zBxIkTaaUxckREJE0VF8PGjRUL5+j7mEeGKtWkSSiWe/SA7t0rvnbvHvo560tcyRQqrA/QwoULuf3223nhhRcAaNy4MRMmTODKK69kxIgR6uohIiJp48svYfly+M9/wmt0fcWKMMpGVRo3hp49oVev0uXQQ8O2Hj3g4IPVLUMalpQprM3sDOB/gYOAJ4Hr3L3IzDoCvwS2Ak2AG9y9ho8wxN+GDRu45ZZb+Mtf/oK7k52dzcUXX8z111+vcaZFpMFIl5wtpYqLQ0vzkiWlxXO0gC4oqPpzBx8ciuXyxXOvXtC1ayiuRSRIif8dzKwHMB44CxgIPAisB+4D/gnc6O7vm9llwC+An9V3jIWFhdx3333ceeed7Ny5kyZNmnDllVfy05/+lEMOOaS+wxERSZp0yNkNXUEBfPBB2WXpUtixo/LjW7SAfv3gsMPCEl3v1y9MeCIiNZMShTXQE7jQ3YuAeWZ2BHCCmb0DHObu70eOexFYbmZ3unsV6SH+Fi9ezA9/+EPmz58PwJlnnsndd99Nnz596isEEZFUktI5u6H5/HOYMwdmzw6v8+eH4ekq06ULHHEE9O9ftoDu2lX9nEXiISUKa3d/u9ymDUBr4JvA6pjjNkT6LucBM+ohLv70pz9x9dVXU1hYSM+ePXn44Yc58cQTE31pEZGUlao5uyEoKoIFC2DGDHj//VBIr11b8bhWrUIBfeSRMGhQeD3iiDDRiYgkTkoU1pUYBtwLTALKD9azHehS/gMFBQXk5eWVvJ80aRKTJk064AB2797NZZddxl/+8hcALrvsMu655x5at259wOcUEYnKz88nPz8/+rZjMmOJg6Tn7EwVLaTffDMsb78N27eXPaZlS8jLg2HDYPhwOPpoyM3VQ4Mi8VTTnJ1yhbWZ9QEK3H2RmTmwq9wh2UCFB2FycnKYO3duXGL4+uuvGT9+PNOnT6d58+Y88sgjTJw4MS7nFhGBsoWkmW1OcjgHLBVydqbZsgVeew1efjm8fvll2f19+sDo0XDssaGQPvxw0DQJIolV05ydUoW1mTUCLgVujGzaAPSL2W9AS2BjomLYsWMHJ598MjNnzqRTp0689tprDBkyJFGXExFJW6mQszPF+vXw9NPwwgswa1YYwSOqVy/45jdhzJhQUHfrlrQwRWQ/UqqwBn4C3OfueyLvpwKXxezPBfYCCWnm2LdvHxMmTGDmzJl0796d6dOn07dv30RcSkQkEyQ1Z6e7LVvg//4PnnoqdPGIatIETjgBvvWtsPTrV/U5RCS1pExhbWY3EZJvczPrBYwB3gI+NbMj3H0JcArwe3ffm4gYbr31Vl588UXatWvH1KlTVVSLiFQhFXJ2upo/H/7wh1BQRydfadYMvv1tOPdcOOkkDXEnkq5SorA2s1uAO8ptXubuj5rZucCtZrYGaAXcnIgY3nnnHe6++26ysrJ47rnnOOywwxJxGRGRtJcKOTsdvf02/PznYUSPqJNOgvPPh+99T8W0SCZIicLa3e8E7qxi31rgokRev6ioiEmTJlFcXMxNN93EmDFjEnk5EZG0luycnW4+/hiuugpefTW8b9sWLrwQLr88PIgoIpkjJQrrZHv00Uf56KOP6NWrF7feemuywxERkQzgDg88ANdfD7t2QevWcN11cM01ap0WyVQNvrDet28fd911FwC/+tWvyM7OTnJEIiKS7oqK4MorITrs7cSJcP/9kJOT3LhEJLEafGH90ksvsWbNGnr16sVZZ52V7HBERCTNuZcW1c2awWOPwXnnJTsqEakPDb6wfuSRRwC48soraaQR9kVEpI4efri0qJ4yBY4/PtkRiUh9adATnm7dupUpU6aQlZXF+eefn+xwREQkza1eDddeG9b//GcV1SINTYMurCdPnkxhYSFjxozh4IMPTnY4IiKS5q65BnbsgHPOCf2qRaRhadCF9UsvvQTAmWeemeRIREQk3U2dCpMnQ6tW4UFFEWl4GmxhXVxczBtvvAHASSedlORoREQknRUVwY9/HNZvuQUOOSS58YhIcjTYwnrhwoV88cUX9OjRg969eyc7HBERSWOPPgoffgi9epUW2CLS8DTYwvrf//43AGPHjsXMkhyNiIikq5074bbbwvpdd4GmQxBpuBpsYZ2bm8upp57KKaeckuxQREQkjbnDJZfA6NGg6RBEGrYGO471mWeeqYcWRUSkzlq2hNtvDwW2vgAVadgabIu1iIhIPKmoFhEV1iIiIiIicaDCWkREREQkDlRYi4iIiIjEQYMprPPz85MdQtrQvaoZ3aea072SeNG/pYp0TyrSPalI96SiRNwTFdZSge5Vzeg+1ZzulcSL/i1VpHtSke5JRbonFamwFhERERFJUebuyY4hLsysAFhbzSEdgc31FE66072qGd2nmtO9ql5Pd89JdhD1qQY5uyr6t1SR7klFuicV6Z5UdKD3pMqcnTGFtYiIiIhIMqkriIiIiIhIHKiwFhERERGJg4wvrM2suZn93sx+ZWYPmtlByY4pFZnZGWa22sy2mNlvzaxxsmNKdWbW1MwWmdmYZMeS6sxshJn9xMzGm1nHZMcj6Ue5vCLl7copN1dOebiUmfU3sz+a2TVm9oCZDYnXuRvC/4QPAK+5+9NmdgrwEDAhyTGlFDPrAYwHzgIGAg8C64H7khlXGrgeyE12EKnOzH4I9Hb3W5Idi6Q15fIYytvVUm4uR3m4gieA77n7hsj/S1OA/vE4cUYX1mbWBZgIXBnZNA141sxy3X1N0gJLPT2BC929CJhnZkcAJ6AEXSUzGwV8CnyZ7FhSmZmNBr4PnJjsWCR9KZdXSnm7EsrNFSkPV6o/0DqyvgtoE68TZ3pXkDHAZnffCRBJQKuB0ckMKtW4+9uRexO1AViXrHhSnZm1BM5y90eTHUsa+A3wEfB7M3vNzI5NdkCSlsagXF6G8nZFys1VUh6u6BngETNrDZwPXBWvE2d6Yd0V2FJu23agSxJiSSfDgD8lO4gUdgNwd7KDSHVmdhhwFPCwu/8ImA68ZmadkhuZpCHl8v1T3lZurkB5uEpXALuBOcDX7v5svE6c6YW1E5r4Y2UDhUmIJS2YWR+gwN0XJTuWVGRmJwNz3f3zZMeSBgYCX7j74sj7PwAGfC95IUmaUi6vhvK2cnM1lIcr1xx4Cvg7cL+ZjYvXiTO6jzXhq7F25ba1BjYmIZaUZ2aNgEuBG5MdSwr7CXCUmUXftwMmm9kv3f2e5IWVkhoTk2PcfZeZrQA6JC8kSVPK5VVQ3i6h3Fw55eHKPQmc5+5bzcyBpyLPbOyo64kzvcX630AXM2sOYGZNgB7AjKRGlbp+Atzn7nuSHUgKmwgMiVk2Ahejr2ArsxhoU+4rxyLgwyTFI+lLubxqytuBcnPllIfLiQw1ONjdt0Y23QlsI06jgmR0Ye3unwEvAWMjm8YBz7r7huRFlZrM7CZgAdDczHqZ2YWRrxclhrsXuPsn0QXYR/gKdluyY0s17r4MeJUwHBhm1hZoCryczLgk/SiXV055u5Ryc+WUhyv1BbDbzLrGbNsC/CceJ8/0riAAlwP3mNlgoDMwKcnxpBwzuwW4o9zmZXqyWuLgB8BvzawpoYVxQrmRDERqSrk8hvK21ILycAx3Lzaz8cCtZjYP6ARcH68/wszd43EeEREREZEGLaO7goiIiIiI1BcV1iIiIiIicaDCWkREREQkDlRYi4iIiIjEgQprEREREZE4UGEtIiIiIhIHKqxFqmBmY81sppmN2c9xR5jZ/5nZz+N47XvjdS4RkWQws0Fm9s945sYqrvOSmZ2XyGskipmNM7NZ+/s9I+mjIUwQI2nOzPoDvyTMljQMGAT8f+7+ZIIvvR4YXIPj1gHtgEbxuGhkuubCeJxLRATAzIYBdwHHA68DqwkT7TQBHnT3KQm47EqgDXHKjdV4AFiSiBNH8vEZhKnR/wPMiuzKBs4Ehrv7mjpcoqa/ZyRNqLCWdPACcI67LzKzLOCx+riou//HzLbU4LhtZvZpHC89GngrjucTkQbO3eeY2T+Aw93929HtZvYt4Hkz+4m7PxDna+6Mc24EwMy+A3wQLWjd/ZV4XyPK3QuBp83sKuBNd78lJo7JB3JOMzsYOMHdn3b35Wa2OU7hSgpQVxBJaWbWHugH7IAwFSlwK1BfU4YWx/m4mjgB+HcczyciAlBhGmt3fxm4Dfh/ZtY5AdeMa642s26EFur6VtkU4NOADbU5iZllA38FmsdsjufvD0kyFdaS0tz9C2A58JqZjYhsWwvMNrO2ZnaHmc02s++a2Wdm9pGZDY1+3swmmNnNZvammf0qZnuemd1mZk+b2XNm1iayvZmZ/dHMbjWzPxO+xqwVM/tpJKZzzGyYmb0QOd8dkRifNrOcyPavzOx/yp2iqbvviZzrOjO71szeNbO7a30DRUT278+ELiFnA5hZbzP7hZk9YmbTI8UsZpYdyWU3m9nbZnZu9ARmdoWZ/dzM3jOzH1V1ocpyspmdFcmZF5jZCjO728y6mNljkXz6tpl9M3KKiUA34AYzGx95Fmaamf1XzDUuMbNfmtlfzewvZtbSzNqY2f+Y2Vwz+6aZLTezVWaWeyA3LHJPxkRatDGzM8zst2b2v2b2mpm1ivwMvzWzFyO/bz6N3OMBwNlmdm3MKdtFfhdtN7NLqrhmlpn92swmRa7hZvZ89L+PpAh316IlpRdgILCW8Ff9o8DBke2NCEnqC0Jft7bAVGAp4Y/GkcAtkWPbEVq9vwUcBPwj5vzvAvdG1v8fMCmy3gnYR0ie+4vxcUKrz8GEvngtItuzgSnAZKAL4RdCIfA/QAvg6Mj76PHR/uPRn/uByHpr4MZk/7fQokVL+i7ABcAnVezbBPwuklefAxpFtv8d+Gdk/TfAGZH1C4GCyPoE4PzIer9Irj4y8v5x4LbIelU5uT2wFbgfGA6MA+4Dbo0cez3wYkysDuRG1psDHwEXRN6fDPw95tiXgIcIXV9Pj/y+OCPyc74G/LKG9+7NyO+Ku4F7Cf2tx8Xs3wz0iKwvJfxOygZuBj4BjgEuiuT9N6PxRo5fA/xvZN8EYF0VMVwMvBRZ7xS5zycl+9+VlrKL+lhLynP3pWZ2JHAncAXwHTMb56HPdQGwzd2fBTCz24B3gF6EXyLFZvbjyKmmEBL4t4AOMds/iny2FXB55LO4+yYzW1+LUDsCvwcucvedkXPsMbPPgNXuvjFync+AGR76Hy4gJPyOhIcgv0X4JQCwE/iBmf0HeJDwR4WISCI4odgcDvQErjIzgC1AGwtvLiU0CgD8hdAVAkKunVsu1x4CfFDuGhdQSU529y/M7CtgsrvPBjCzdcAOMzuE0MLbqtKg3XeZ2ecxm34UOW/UI8D/RbZvJfy+eC5yjbmROGvq3x7pY21mTxMac6LOdvd1ZnYsoaBuFcn/G4GP3f194P3IZys794uR3wlzqolpKLAr8nNvMrMvKNulRFKACmtJaWbWCOju4SGV/450z3gFyCe0AJS3NPLaFugBPOHuf49suz9yzhsIhe795a51NCEhbj/AcPsBxwG9gUUx28v3MdxXssO9OJJko92yOnjo/oK7rzazC4A/ANcA3wc+R0QkjsysLeHbtqWEvLmlkvzYidCi2hTA3fcRGgOIfOZ2d58ZeV/mszEqzckRTtlc+QlwA7ACeA84tJofIfZzfYE3Yt5/TOjmcjAVc3ERBz5iyXzCHyBRW83s18A/CbFHq+fyP9f+FFF1bfYq8KCZNSP8t9hFaP2WFKI+1pLqSvr9Abj7IuAm4PAqjm9KSGIrgU+Bs2J3mtkxke3fMrPm5bZvi7ztf4CxzgR+C/wtkvhqxcy6EPMgjJl1JYyI0o/wMOPfq/ioiEhd/BfhG7J/EvLjsZGWYqAkP24G9gKnxGw/0sxaUy7XmlnjSENFeVXl5Mr8CVjjYVjVPbX4WdYRWrhLLgF8RejqEjfuXhxp/BgY+cPkFeAX7j43ntcpd80XgRcJ36xOAEa7+1eJup4cGBXWkg6uihSdUV0JfeOi2ptZ08j6GEL/uq3AU8DpZvYHMxsTaU3YQUiALYGXzOzkyIMihxOK8QXAHWbW1Mw6Evo2H2xmjczsRjM7v4oYLbLcSvjl85tK9pU/vvz7bxP6A0b1Bca7+zbgqko+IyJSGxVaQs3seOB24FJ3LyC0Dm8CXjGz70QeUPxGpIX6H4TRQ8ab2Vjgv9x9OyHX/reZ3W5mowkNDNFudLH5r6qcHBXbenwUkGNmLYBjgeZm1iuyr5DwsF+0gSX2Gn8CzjSzdpH3xwAPu3sRoeapNBeb2YhIXNlV3LssytVMkT8+vkX4lrIdoYvhYUB3oEXMg5HlW8X3VhJ/2aAq6S8S+W9xEDCd0Cpfmz84pJ6oK4ikg+7AUjN7mfA1mQGXxewvAm42s62E1t2rANx9qpn9N+HrxPGEh/+WAJjZ94A/En5RPApc7+5uZucQhkJaCjxN+EpvCGFChXMIXUzKTExjZsMJD+X0IHw1+AJweySeKcAIr9iGiAAAIABJREFUINfMehNaUjoD55jZMkpbby4AWrp7frmf/TELEzsUER4WEhGptUgemUgo/h4GCggPwLUBTnX3WQDuvjeSHx8m5LrJhL7VAFcTGiX+SviGLjoSx58Jue9yQpe1y93980jhOALobWZPVpWTzewsQr/ii81smbtvIDxIeR+hqP4r8D3gSGAV8AShdf0HFkaLGgicZmavuPtzZtYXeNbM3iV863lzpGX9bKBz5OdbSpgs5xAzGxiJ8zxCv/EXYu5bE0LuHwx0jRTeTujbfEbknAuBGcAcQveWVwk5/SXCA5ODzWy8u0fP+zRhsp6vI4V4Z+AsM1tNaIkm8vnyczasIzycOR5oBmSZ2TTCA4z1NQSt7Ifpv4WkMwvTwD7u7rn1dL3r3V3TjYuIZJhIt5Rsd0/JCbrM7HrCCC1rIu+bEb5t+HX02RxJPnUFkXRXWTeLxFwofGU6vz6uJSIi9SfycObQFC6qmwG3ALkxmzsCO1RUpxZ1BZG0Fenfdj7hq72z3f2fCb7kx5GvKEVEJLM0dvcHkx1EVdx9t4Vp1X8XebB9I6Hbyi3Vf1Lqm7qCiIiIiIjEgbqCiIiIiIjEgQprEREREZE4UGEtIiIiIhIHGfPwYseOHT03N7dWnykoKCAnJycxATUwupfxo3sZP+lyL+fNm7fZ3VM/0Dg6kJwN6fPftD7oXpSl+1FK96KseN+P6nJ2xhTWubm5zP3/2bvz+Ljqev/jr0+WNmn2tmm60RVoS1klVDbZERDZQVBckAu9158guP3clyt49SrXn15ErxXRK/eiqKAgm7IpIGVpKRQoLRS6plva7E2aNM3n98d30kzStM0yM2eSeT8fj/OYmXPOnPPJmeTTT7/zPd/vov7NJFpZWdnv90jvdC0TR9cycYbKtTSzNVHHkGoDydkwdD7TVNC16E7Xo4uuRXeJvh77ytnqCiIiIiIikgAqrEVEREREEiCjC+v58+dHHcKwoWuZOLqWiaNrOfzoM+2ia9GdrkcXXYvuUnk9Uj5BjJmdCdwMXB433/3FwH8AxcD/AJ919/bYtrHAt4E6IBf4grvv7HncyspKV38iERmKzGyxu1dGHUdvlLNFRLrbV85OaYu1mVUAhcC8uHVTgAuBS4FPA9cAN8a97ffAHe7+BeBN4FspC1hEJIMpZ4uI9E9KC2t33wzc12P1VOBqd1/s7r8GbgNOBTCzY4FZ7v58bN8/A9ebWUGqYhYRyVTK2SIi/ZPyPtbu3tHj9dOdXyHGVAFrY89PA1bF7VsVe5rYr0wvuAAOOghqaxN6WBGRoS5dcnbj4kYWH7OY5VcvH+yhRESSJh3HsT4G+H7s+SRgW4/tjcDEnm+qrq6msrIrd8+fP7/vndVXrgzLunVQVjaQmEVE+mXBggUsWLCg8+XYKGMZpJTk7I6dHTQuaszwW+5FJCp9zdlpVVib2YFAtbu/ElvlQEuP3UYCe9wIU15ePvDBv6dMgWXLYO1aOPzwgR1DRKQf4gtJM9sacTgDksqcnVMU/rna1bhrQLGKiAxGX3N22vzf38yygX8Gvhi3ugooi9vHgAJgQ0JPfsAB4XHduoQeVkRkuEp1zs4uzAZgV5MKaxFJX2lTWAOfA25x99a4dY8Cs+JeTwPagMSO0dRZWK9du+/9RESkU0pzdnZRrLBWi7WIpLGUF9axFgwAi1v3ZWAJkG9mM8zsajM70N0XARvN7NDYrmcDt7p7W0KDmjIlPKrFWkSkm3TJ2Z0t1u2N7aR6/gURkb5KaR9rMysCPhx7eZWZ/ZjwVeJNPXZd7u53xJ5fDnzdzFYTxlP9SsIDU4u1iMge0ilnZ43IwkYY3uZ0tHaQnZediMOKiCRUSgtrd28EfhpbOt0cW/b2njXAPyU1sKlTw+OqVfveT0Qkg6Rbzs4uyqZ9Wzu7GnepsBaRtJROfayjM3UqjBgB69dDU1PU0YiISC9238CoftYikqZUWAPk5IQJYgBWrIg2FhER6dXuIfc0MoiIpCkV1p3mzAmPb7wRbRwiItIrjQwiIulOhXWn2bPDowprEZG0FD8yiIhIOlJh3emQQ8Lj0qXRxiEiIr3KKQldQdrrVFiLSHpSYd3pmGPC4wsvgMZIFRFJO7ljcwFo36bCWkTSkwrrTjNnwpgxsGULrFkTdTQiItJDzpjQYr1z286IIxER6Z0K605m8O53h+fPPRdtLCIisofcMaHFWoW1iKQrFdbxTjghPD75ZLRxiIjIHjq7guzcqsJaRNKTCut4731vePzLX9TPWkQkzajFWkTSnQrreO96F4wdG/pYv/lm1NGIiEgctViLSLpTYR0vKwvOPDM8f/jhaGMREZFudrdYq7AWkTSlwrqnCy4Ij3ffHW0cIiLSTW5FrLDevBNXdz0RSUMqrHs67zwoKAgjg7z9dtTRiIhITE5hDtkl2XTs6FA/axFJSyqsexo1Ci66KDy/885oYxERkW7yDsgDoHVda8SRiIjsSYV1bz7+8fD4s59BW1u0sYiIyG4jJ48EoHW9CmsRST8qrHtz6qlw6KGwaRP87ndRRyMiIjEjD4gV1mqxFpE0pMK6N2Zw443h+fe/Dx0d0cYjIiKAWqxFJL2psN6bK6+EyZNh6VL4/e+jjkZERIC8qaGPdcs7LRFHIiKyJxXWe5OXB1//enj+9a9De3u08YiICPmz8gFoWaHCWkTSjwrrfbnqKpg5M8zC+F//FXU0IiIZb9SsUQA0r2jGOzSWtYikFxXW+5KbC7fcEp5/5SuwYUO08YiIZLjcslxyx+XS0dKhftYiknZSXlib2Zlm9ryZTYtbl29mt5rZv5nZT82suC/bUuKCC8KkMQ0NcP31oNm+RCSDpGPOHjU71mq9vDnRhxYRGZSUFtZmVgEUAvN6bPoJ8Iy7fxm4D/hZH7clnxnceisUFcG998Ltt6f09CIiUUnXnF1wSAEATUubEn1oEZFBSWlh7e6bCYl2NzObCHwI+HNs1WPA+WY2bV/bUhJwp6lT4ac/Dc9vuAGWLUvp6UVEopCuObuosgiAxkWNiTysiMigpbwriLv3HBT6FGCruzfHtrcDq4CT97Mtta68Ej72MWhpCVOe19amPAQRkVRLx5ytwlpE0lU63Lw4CdjWY10jMHE/27qprq6msrJy97JgwYLER/rjH8Phh4dRQi6/XEPwiciALViwYHe+AsZGHU8/RJ6zR80dRVZ+Fjve3sHOmp39i15EZAD6mrNzUhfSXjnQc0DSkUBnttzXtt3Ky8tZtGhR4qOLV1gI998PxxwDjz4aZme89dbQD1tEpB/mz5/P/PnzATCzrRGH0x+R5+ysnCwKjyqk4dkGGp5rYMz7xgzoOCIifdXXnJ0OLdZVQFmPdUXAhv1si8bUqfDHP8KIEXDbbfCNb0QWiohIBNIiZ5eeXApA7RPqlici6SMdCusngYlmlg9gZrnAFODv+9kWnRNOgN/+FrKz4aab4HvfizQcEZEUSoucXXZ6qN9rH1NhLSLpI4pxrDv7TRiAu28CHgBOj60/A7jH3av2tS2FIffuoovgV78Kz7/wBfjXf9UY1yIy7KRrzi4+oZisvCy2v7Kdtuq2RB9eRGRAUj2OdRHwL7GXV5lZeez5J4ALzOwrwPuA+XFv29e2aH34w/CLX0BWFnzzm/CpT0FHzxvoRUSGpnTO2dl52ZScWAJAzcM1yTiFiEi/mQ+TVtbKykpP+s2Le3PvvfDBD0JbW2jJ/vWvw42OIiJ9YGaL3b0y6jhSKRE5u+onVbz1ybcYc/4YDrvvsARFJiKyb/vK2enQx3rou/hiePhhKCkJNzYedxy8807UUYmIDGtjLxoLBjV/qaG9QcOfikj0VFgnymmnwfPPw6xZ8NprYUi+hx6KOioRkWFr5ISRlJxYgrc6W+8fSiMWishwpcI6kWbNCsX1uedCTU14/PSnobU16shERIalcR8aB8DG2zdGHImIiArrxCspCZPIfPe7kJMDP/whHHssvPFG1JGJiAw7FR+qIKsgi/q/19O8ojnqcEQkw6mwToasrDAE3z/+ATNmwMsvw5FHws03hxscRUQkIXKKc6j4YAUAG34W3dxhIiKgwjq55s2DJUvgmmtCQf21r8HRR4fuIiIikhAT/2UiABt/vpGdNXvMni4ikjIqrJOtuBh+/nN4/HGYOTPc2HjccaHY3rw56uhERIa8oqOLKHtvGbuadrH+P9dHHY6IZDAV1qly2mnw6quhi0hOTphY5qCD4JZb1D1ERGSQpn5lKgBV/1lFe72G3hORaKiwTqX8/HBT42uvhRFDGhvh85+HQw+F++7TlOgiIgNUelIpJSeV0F7bzppvr4k6HBHJUCqso3DwwfDAA2Gc61mz4K234MIL4cQT4emno45ORGRImnnLTADW/2g9Le+0RByNiGQiFdZROuccWLo0DMk3diw8+yycdBK8//1hvYiI9FnxMcVUfLQCb3NWfnolrm8BRSTFVFhHbcQIuOGGMAX6N78JhYXw4INheL7LL4dXXok6QhGRIWPGv80guyibbfdvY8vdW6IOR0QyjArrdFFUBN/4Brz9dii0c3Phd78LBfZ558HChVFHKCKS9kZOGrm7S8hb171F2xbdHC4iqaPCOt2MGxe6hnQW2Pn5oT/28cfDqaeG5x0dUUcpIpK2Jlw7gdLTS2nf1s7yq5fjHeoSIiKpocI6XU2eHArsNWvgy18O42H/7W+h9frgg+FHP4KGhqijFBFJO2bG7F/MJqcsh5oHa1j3/XVRhyQiGUKFdborL4dvfxvWrg1jXk+bFlqzb7wRJk2CT30KVqyIOkoRkbSSNzWP2b+eDcA7X3mH2r/VRhyRiGSChBTWZva+RBxH9qGkBD77WVi5Eu69F045BZqa4NZbYfZsOPlkuPNOaNEQUyKyb5mSs8e+fyxTvjgFdsHrF7/O9uXbow5JRIa5fhXWZna2mT1jZivN7J3Yshq4JznhyR6ys+Gii+DJJ+Hll+Gf/glGjYKnnoKPfhQmTIDrrgvbRCSjKWfD9JunM+a8MbTXtvPq+16lbbNuZhSR5LH+jPNpZtXAd4AlQOcddAZ8yN3nJz68vqusrPRFixZFGUJ0Ghrg7rvh9tvhhRe61h91FFx5JVxxReg2IiJpycwWu3tlEo6rnA3s2r6Ll099mcYXGyl8VyFHPHYEuWW5KTm3iAw/+8rZ/S2sn3f3d/eyfry7bxpEjIOW0YV1vKVL4Re/CN1CamN9Cs3CiCJXXgmXXBK6lYhI2khiYa2cHdO2uY0lJy6hZWULhUcXcsSjKq5FZGD2lbP728f6E2Z2aS/rz+9/WJIUhx8eRgzZsAHuuQcuvjiMif3EE6HbSEUFXHppaOFubIw6WhFJLuXsmBEVIzjiySPIm5lH0+Imlr53KTu37Yw6LBEZZvpbWJ8A3G5mu+KWDuCngw3EzOaY2W1m9mkz+4mZHRlbn29mt5rZv5nZT82seLDnygh5eaGovuce2Lw5dBM59VRoawvrrrgijDhy/vnw3//d1botIsOJcnacvMl5HPnkkeTNyKNxUSMvnfASLat1w7eIJE5/u4LUAv8MLAJ2da4GPuLuNw0qELNFwAXuXmVmU4C/uPscM/sl8Ii7321mZwMfc/cP9ny/uoL00fr18Pvfh5FF/vEP6Pz8c3JC4X3JJaHYnjAh2jhFMkgSu4IoZ/eitaqVpecsZfur2xkxfgSHPXQYRUcVRRKLiAw9iexjfRfwcXdv7bF+0P31zGw7cLS7LzezcuAVoBJYBZS5e7OZ5QD1wFx3Xx3/fhXWA7BxI/zpT6HIfvJJ2LWra9u73gXnnhuWysowGomIJEUSC2vl7L1or2/ntYteo+7JOrJGZTH7jtmMu3xcZPGIyNCRyD7WvwW+ZGYnxS0nA18adJTwB8JXlkXAh4HrgVOAre7eDODu7YSkfXICzicTJsAnPgGPPhq6i/zyl/D+94dp1F96CW66CY49FsaPD0P53X23uoyIDC3K2XuRU5LD4Q8fTsVHKuho7mDZFctY+bmVdLR37P/NIiJ70d8W6xeBg4A6oPONBox397xBBWJWANwHTAb+w91/bmafJ3xleXjcfguB+939O/Hvj7r1Y1hpaQnTpz/4YFhWr+7alp0N8+bBGWeE5dhjYcSIqCIVGRaS2GKtnL0f7k7VbVW8/em38Xan5KQS5vzPHPIOGNTlEZFhbF85O6efx/o+8LC7dxtOIkGzeOUDvwEmAT80s1WEfwh63lkyEtjjVu7q6moqK7t+xvnz5zN/fqTDtA5d+flwzjlhufVWWL68q8h+5hlYuDAsN90UJqc56aSuQvuwwyArIRN6igxrCxYsYMGCBZ0vxybpNMrZ+2FmTL5uMoVHFrLssmXUP1XPi4e9yME/OZiKD1WkPB4RSU99zdn9bbE+zt0X9rL+ane/o99Rdj/GI8AV7l5nZl8DPgXcCHzD3Q+O2++t2Lq74t+fLq0fw15DQ5jl8bHHwvL66923l5fDaaeFYvukk+CQQ1Roi+xHEluslbP7oW1LGyuuXcG2+7cBMO6KcRx464GMGKtv5USkSyJvXrzT3T/SY91lwC/cfcBDKpnZWOBVd58Qe23ASuBjwCNAubu3mFku0ATMcPeq+GOkY5LOCBs3hjGyOwvt9eu7bx89Gk48savQPuqoMAKJiOyWxMJaObuf3J2Nt29k5Y0r6WjuIGdMDgf+x4FUfLSC8GOKSKZL5M2Lo2MtE5jZKDO7HbgdeHGQMdYAO8wsft7tbcBS4AHg9Ni6M4B7eiZoidCECWFGx1/+EtauDd1GfvIT+OAHwzTqNTVw//3wuc+FvtmlpfDe98LNN4d+3Nu3R/0TiAxnytn9ZGZMvHYila9UUnpqKe3b2ll+1XJeOf0Vmlc0Rx2eiKS5frVYA5jZscCVwFnAJuCjQL27D2q4CDM7Avg/wGKgAnjK3f9uZmXA94DVwHjgS+7e1PP96dr6kdHcw42PTz0FTz8dHt96q/s+2dlhtshjj4XjjguPBx4YpmEXyRDJarGOHVs5e4Dcnc13bmblZ1bSvq0dyzEmXTeJqV+bSu5oTYcukqkG3BXEzE7ay6Z5wKnALUAHcKG7f3qwgQ7GUEjSAmza1FVkP/ssvPJK9/GzAcaM6V5oz5sHRZq8QYavRBXWytnJ0ba1jXe++A6b7tgEDjllOUz7xjQmfmIiWSN0D4lIphlMYf0cMIcwVNPeBvc0YIK7jxxsoIMxlJK0xNm+HRYvDqOMPPdceNy8ufs+ZnDooXDMMWGprAyjj4yM9FdOJGESWFgrZydR48uNvP3Zt6l7og6AvBl5TP3qVCo+UkFWjgpskUwxmML6fcAid9+ynxNc6u5/GFyYgzMUk7T0wh3WrOleaC9ZAu3t3ffLzQ1dSCoru5a5c8N6kSEmgYW1cnaSuTvbHtjGO//3HZqXhz7XeTPzmPa1aYy7cpwKbJEMkLBRQdLZUE3S0gctLaG4XrQoLC++CCtWhCI8Xl4eHHlk92J79mxNxy5pL5l9rNPVUM/ZHe0dbPnNFtZ8aw0tK8PQ3fkH5XPA5w+g4iMVZOcp74gMVyqsZfhpaOhebC9aBCtX7rnfqFFhiL+jjgpF95FHhm4l6kYiaUSF9dDV0d7Blru2sPpbq9nx9g4AcsflMun6SUz6xCRyx+hbNJHhRoW1ZIba2tBfO77YXrNmz/1ycmDOnFBkxxfcZWWpj1kEFdbDQUd7B9W/q2bd99fR9HIYBCUrP4vxV49n8vWTGTVrVMQRikiiqLCWzFVdHVq2lyyBl18OS2/dSACmTu0qsjsL7ilTNPSfJJ0K6+HD3al7so51319HzSM1u9eXnlbKxE9MZOwFY8nKVT9skaEsqYW1mR0FfBZYQpgIYPWgDjhAwzVJSxJs3w6vvhqK7M6Ce+lS2LFjz32Li0PXkcMPDyORHHZYeK3WbUmgVBbWytmp0/RaE1U/qmLzXZvpaA6DtIyYMIIJ10xgwrUTyDsgL+IIRWQgkl1Y/wy4wN3Hm9mHgD+6e8ugDjoAmZCkJYna2+HNN7tatTuL7q1be99/8uSuQrtzmT1bfbdlQFJcWCtnp9jOup1svnMzG366geY3YrM3GpSdUcb4j41n7EVjyR6lmx1FhopkF9a5wDh3rzKzLKDQ3RsGddAByKQkLSniHsbUfvXV7svrr/feup2dDbNm7VlwT50KWfrqV/YuxYW1cnZE3J36p+qp+mkVW/+4FW8L//5mF2VTflk54z82npITS7AsdT8TSWcJLazN7BxgtLv/r5nNA5rcfVkC4hyUTEzSEpFdu+Dtt/csuFeu7L3vdmFh6D7Ss+AeMyb1sUtaSvKU5srZaWhnzU623L2FTf+9icbnG3evz5uex7grxlH+gXIKjyjEdI+HSNpJWGFtZt8GbgAecfdLY+u+Cyx197sSEexAZXqSljTQ3AzLlu1ZcPecSbLT+PFhUptDD+1aDjkk9OuWjJKswlo5e2hoXtHMpl9vYvOdm2ld17p7ff7B+Yz7QCiyCw4tUJEtkiYSWVg/CVwEfMzdfxRbdzJwp7tPSUSwA6UkLWmrurr37iTbt/e+/5QpXYV2Z+E9Zw7k56c2bkmZJBbWytlDiHc4dU/VUf27aqr/UM3O6p27t42aPYryD5RTfkk5BYepyBaJ0r5ydk4/j/Wou9eZWXw1funAQxPJAOXlcNppYenU0RHG2H79dXjtta5l+XJYuzYsDz3Utb8ZzJzZvXV77lw4+GAYMSL1P5MMFcrZQ4hlGWWnlFF2ShkH/ueB1P+9ni2/20L1PdU0L29mzbfWsOZba8iblseY88cw9oKxlLynRMP3iaSR/rZYXwHMACYBTwOXEVpDrnf325ISYR+p9UOGhfb20H+7s9DuLLzffDP07e4pJyfcMBnfun3ooTBjhqZyH0KS2GKtnD0MdOzsoO7JOrb8bgvb/ryNnVu6WrJzSnMYfe5oxl4wltFnjyanqL/tZSLSX4m+efFY4CpgCrAR+I27PzbYIAdLSVqGtdbWMLFNfAv366/DO+/0fsNkXl7oPtKzhVsT3qSlJN+8qJw9jPgup+GFBrbet5Vt922jeXnz7m02wig9uZTR54xm9NmjGTV7lLqMiCRBsofbywZOcffHB3WgQVKSloy0fTu88Ub31u3XXoP163vfv6goFNjxrdtz54YbKfUPcGRSPNyecvYw0vxmM1vvD0V2/bP10NG1beSUkYw+ezSjzxlN2Wll5BSrNVskEQZcWJvZ/UDhfo4/Ftjq7qftZ7+kUpIWiVNXF0Yoie+//frrsGVL7/uPHr1nd5K5czUkYIokqrBWzs5sbdVt1D5aS80jNdQ8UtPt5kfLMYpPKA6F9lmjw1B+Gi9bZEAGc/Pi0tg+Kwj/D76WMA1ufDacDoxOQJwikiilpXD88WGJt2VLV8t2fAt3TQ089VRY4o0fv2fBrSEB05lydgYbUT6Cig9VUPGhCrzDaVrSRM0jNWx7eBsNCxuo/3s99X+vZ9WXVpEzOofSU0spO72MstPKyD84X91GRBJgfy3WJbF96mKvv+XuX+9lv/vc/YLkhbl/av0QGSB32LBhz2J72bL9DwkYX3BrSMABS2CLtXK29Gpn7U5qH6+l5uEaah+rpXVta7ftIyaNoOy0MkpPK6XstDLypuRFFKlI+htwi7W71/dYVWFm2e6+e3gCMzsFmDfoKEUkGmYwaVJYzjqra33nkIA9C+6+DAnYWXDPnRtGLdGQgCmhnC17k1uWy7hLxzHu0nG4Ozve2UHtE7XUPl5L3RN1tFW1sfnOzWy+M0xolX9gPqWnlVJ6cikl7ykh7wAV2iJ90d87Gf4AvGJmzwCNwBzgLGCPFhERGeKysmD69LCcd17X+vghAeML7jffDNO6r1wJf/pT1/45OWG87Z79t2fODNskmZSzZQ9mRv7MfPJn5jPx2om4O9tf207dE3XUPlFL3d/qaFnZQsvKFjYu2AhA3rQ8St5TQslJJZS+p1RdR0T2YiDD7U0E/g9wONAA/N7d70tCbP2irxVFItbaGorrnkMCvv1270MCjhwZuo/0LLinTg1FfQZJ8nB7ytnSLx3tHTS91ETtE7XUP11P/T/q2VXffRz93HG5lLynhNKTQot24eGFWLYKbckMyR5uLws4yt0XD+pAXcc7FjgRWAk8A2wHvkdobSkDvuDuDT3fpyQtkqaam8OQgD0L7rVre9+/oKCr2D7yyLAcccSwvmEyxcPtKWdLv/gup+nVplBkP1VP3dN17Ny8s9s+2cXZlJxQQvHxxZQcV0LRvCJNViPDVsIKazO7o5fV5cAud79wgPHFH//jwEx3/2rcul8Cj7j73WZ2NvAxd/9gz/cqSYsMMfX14QbJngX3pk2973/ggXDUUaHQPuqosIwfn9qYkySJMy8qZ0vCuTstK1t2F9n1T9WzY9WO7jtlQcFhBZQcV0LxccUUH1dM/oHqPiLDQyIL66eBx4H4N50APO3uNw0yyJOBrwFneiyo2FeYq4Ayd282sxygHpjr7qvj368kLTJMbN0aCuxXX4UlS+Dll0PR3da2574VFV1Fducyc+aQm+wmiYW1crakRGtVK3VP19GwsIGGhQ00LWnC27vXF7ljc3cX2cXHFVN8TDHZBdkRRSwycIksrGe5+4oe62YDV7n7FwcbJPAsYMCBwE3AVOD77j4pbr/XYuv+O/79StIiw1hbW+hO0llodz427NHDIEx2M29eWN797vA4dmzqY+6HJBbWytkSiV0tu2hc3EjDs6HQrl9Yv0f3EbKh8PDCUGQfW0zxvGLyD8rXxDWS9pLdx3o68KK7D/hfLjObBbwBHOnuS83s84S71m8GrnT3w+P2XQjc7+7fiT/G1KlTvby8fPfr+fPnM3/+/IGGJCLprqMDVq8ORXZnob1oEWzevOe+M2Z0Fdnvfje8613h5skILViwgAULFgCwePHiNe4+LRXnVc6Z6yPiAAAgAElEQVSWKLg7O1bvoOHZUGQ3LGyg6ZUm6H5PJNkl2RQfU0zRvCKK54XHkROi/VsVgb7n7P62WK+i+1eK2UAFcKe7XzvQYM3sYmBBZ6I3s3ygGvgWcIm7vztu35eAu9z9lvhjqPVDRHCHdevg+efhhRfC4+LF4QbKeHl5cNxxcPLJYXn3uyOd3CaJLdbK2ZK2dm3fReOiRuoX1tP4QiMNLzTQVrVnl68Rk0bsLrKL5xVTVFlETrFujJToDGZK854eBe6iK1F3AJvd/c1BxNcZx+5Y3L3FzN4i/CNQ1mPfImDDIM8nIsORWZgVcsoUuOyysK69PfTZ7iy0Fy4MN00++WRYIExgM28enHoqnHNOeJ49LPp+KmdL2souyKb05DAJTafWqlYaXmzYXWg3vthIW1UbW/+4la1/3Bp2Mhg1e1QotGOt24WHF5I1MrOG6ZT01N8W65Pc/ake62YR7jBfOeAgQp+/N4Dx7r45tu5F4MfAbUB5LHHnAk3ADHevij+GWj9EpM+qq+Hpp+Hvfw/L0qXdx9oePRrOPjsU2eecA2PGJDWcJLZYK2fLkOYdTstbLaHIjhXbTS834W3daxcbYRQeWUjRMUUUVYZl1OxRZOWo2JbES2SL9dnAUz3WbQaeBI4aQGwAuPtyM3sYuBS4zcxKgRHA/wLnAKcDDwBnAPf0TNAiIv1SXg4XXxwWgNpaeOYZePRRePBBeOcduOuusGRnw5lnwhVXwIUXQklJtLH3j3K2DGmWZYyaNYpRs0Yx/iNheM2O1g6aljbR+GLj7oK7eXkzjS800vhC4+73ZuVnUXhU4e5Cu6iyiFEHj9JENpJUfWqxNrPrCDelFO1ll0fd/exBBWI2FvgRsAiYAvzc3ZeZWRlhsoHVwHjgS+7e1PP9av0QkYRwh7fegoceCssTT8Cu2B1WI0eGYvxTnwr9shM0rF+iW6yVsyXTtNe307i4kcZFXcseY2sDWQVZFL2rq9AuOrpII5FIvyVkVBAzOxz4GPCfPTa1ufvGwYU4eErSIpIUW7fCH/4Av/0tPPVUV5eRykr42tfgvPMGXWAnoyuIcrZkup3bdtL4Uvdiu3Vt6x77ZRdlU3R0KLQLjw4t3PkzNZmN7F0ix7HOdfedvayf4u57mZ84NZSkRSTp1qyBn/4Ufv5zqKkJ6048Ef7rv8I07AOUxD7Wytkicdqq27q1bDctbqJ1/Z7Fdk5pTiiyj+5q3c6blqdiW4BBFNZm9jngIXdfFnv9FcJd3912A05191MSE+7AKEmLSMq0tMCCBXDTTbBtWxi+7yc/gY9/fECHS1RhrZwt0n+tm1ppWtwUiu3FjWEkkk17DvuXMzqnW6FdVFnEyANGqtjOQIO5eXE6UBr3+l1AAbCJruGbjNCPTkQkM+Tnww03wFVXwWc+A3fcAVdfDbm58OEPRxmZcrZIP40cP5KR545kzLldo/+0bmjtKrQXhWJ7Z/VOah+tpfbR2t375Y7N7VZoFx5dyMhJKrYzWX+7grwLeM3d23qsP6SzhSQqav0Qkcj8v/8XCuzSUqiqglGj+vX2JHYFUc4WSQB3p3V9j2J7USPt29r32De3onuxXVRZxMjxmj1yOEnkcHuHuftLPQ4+C5gJRJqkRUQic+ON8IMfwPr1sGIFHDXgkewSTTlbJAHMjLwD8sg7II/yi8qB2DTta3Z0dSOJLTs376TmwRpqHqzZ/f6RB4wMY2wfE5s98ugicko0e+Rw1KdP1czeQ/j68MzYFLnxRhPuOv9zgmMTEUl/u3bBv/97KKpHjYKZM6OOSDlbJAXMjPxp+eRPy6f8krhie9WOrkL7xdDC3bquldZ1rWy9d+vu9+fPyg9FdqzgLjyykOy8YTHja0br63+X6oBfE/rvndBjWxtweyKDEhFJe+7w8MPwjW9AZ5eG//xPKC6ONq5AOVskAmZG/ox88mfkM+4D44Awe2TziuauCW1ebKTp5SZaVrTQsqKFzXduDu/NMQoOLwit2rGp2gsOKdCENkNMnwprd3/VzE4BLnN3JWQRyVx1dWFGxgUL4JVXwrpJk8IwfOedF21sMcrZIunDsoyCOQUUzClg/Edjs0e2dc0e2VlwNy9rpumlJppeamLjz8JQ81mjYhPazCvaXXDnzdCwf+mszx183L2evbRymNmJ7v5MwqISEUknLS3wl7/A738Pf/xjeA1QUQGf+xz8y79AYWG0MfagnC2SvrJGZFFcWUxxZTF8Iqxrb2yn6aUmGl4MrdqNLzSyY/UO6p+pp/6Z+t3vzRmdQ9ExRZQcV0LxccUUv7tY/bXTyD4/CTO7H9jfvxYGzAXGJSooEZHIbdkCjz0Gf/4zPPAANMXNyn366XDNNXDRRWGa8zShnC0ydOUU5VB6cimlJ3eNmNlW3Rb6ar/QuLvg3rllJ7V/qaX2L7Fh/wxGHTKqq9A+rphRs0ZpmvaI7O+/OCsIffWq6BoDtScD3p/IoEREUq6+Hp57Dp58Ev76V1iypPv2ykq49FK47DKYMSOaGPdPOVtkGBlRPoIx54xhzDlhjG13p3VtKw3PN9CwsIH6hfU0vdRE8+vNNL/ezMbbQxeSnLIcio8t3l1oF88rJqdYrdqpsL+r/O9Ao7vvOd9nHDPT3eUiMnR0dMDKlaGQfvbZsLz2WrghsVNeHpx0Epx1VmiZnj49unj7TjlbZBgzM/Km5pE3NW/3zZG7duwKXUhihXbDwgbaNrRR83ANNQ/HhvwzKDi0gOLjiik5voSS95SQN119tZNhn4W1u2/d1/Y4F6G7zEUkHbW2wuuvhxbozuWVV2D79u775ebC0UfDiSfCe98bHvPzo4l5gJSzRTJPdl52KJaPL+EADgit2utauwrtZxtoWtLE9le3s/3V7WxcEFq1R0waQel7Sik5qYTSk0oZNUfdRxJhUN8LmNkU4JPAdShJi0iUWlvhzTdh2bKwvP56eHzrLWjfc3Y0Jk8OhfQJJ8Dxx4fneXmpjzuFlLNFhj8zI29KHnlT8hh3eaxVu2UXjYsbQ7H9j3AzZFtVG1t+u4Utv90CQM6YnN2Fdsl7Sig8spCsnKwof5Qhqd+FtZllA+cD1wJnApsJffpERJLLHbZuDcXyW291L6RXrgxdPHoyg1mzwmyIncuRR0J5eerjj4Bytohk52dTemIppSeWwudjY2u/0Uzd03XUP1VP3d/raNvQxtY/bWXrn8IXX9mF2RSfUEzpyaWUnVFG0buKNKZ2H/S5sDaz6cA1wMeBMcAa4DzgL8BxSYlORDJTTU1X8dxzqa/v/T1ZWXDQQXDIIWGZOzc8zpoVZkTMMMrZIrI3lmUUzC2gYG4Bk/5l0u4ZI+ueqqP+6Xrqn6qnZWXL7tFHVn15FTmlOZSeVkrZ6WWUnVFG/kH56qPdi/0W1mZ2GaGl4zRgJfAD4L+BL7v7I7Hdnk1ahCIy/OzcCevWwapV8M474bFzeeutUFjvTXFxKKA7l85C+uCDh31Xjr5QzhaR/oqfMXLCVRMAaN3QSt1TddQ9WUftY7XseGcHW+/dunta9pEHjNxdZJeeXsrI8ekz9GiU+tJi3UKYAvdZ4Bp3fxPAzPY2lJOIZDp32LSpe8EcX0SvW9d7t41OBQXdi+f4pbw8dO+QvVHOFpFBGzlxJBVXVFBxRQUALataqH28ltrHaql7vI7Wda1s+tUmNv1qEwCFRxYy+tzRjDl3DMXzijO228h+C2t3fwB4IHbTy7VmNg24D8jt3MfMzN2VtEUyxc6dUFUFa9fuuaxeHZbO2Ql7YxZuHpwxIwxjF78ceCCMH6/ieYCUs0UkGfKn55N/TT4Tr5mIdzhNS5t2F9l1f6+j6eUmml5uYu2315I7NpfRZ49m9LmjGX3WaHLLcvd/gmHC+ptbzSwHuBC4EngN+AXwXndfkPjw+q6ystIXLVoUZQgiw4M71Nb2XjR3Lhs2dB/zuTdjxnQvmOOL6ClT0mrGwqiZ2WJ3r0zSsZWzRSSpdrXsou5vddQ8VMO2B7exY9WOro3ZUHJ8CWMvHEv5JeXkTR36Xfb2lbP7XVj3OPCBwPXAte4e6d1BStIifeAe+i9XVXUt69eHx3XrugrnnmM895SVBRMnhgK5t2X69NAXWvokmYV1j/MoZ4tIUrk7zcub2fbgNmoerKH+mXq8vavWLDqmiPJLyhl7yVhGHTg0byzfV84e1DjW7r4SuMHMFg/mOCKSADt3hn7NnYVyz8K5c9mxY//HKiqCqVP3XjhPnBgmVJEhRTlbRJLNzCiYU0DBnAKmfG4K7fXt1DxSQ/U91Wx7cBuNLzbS+GIj73zxHQqOKGDcZeMY96Fx5E8fWhNy7U1CJo53918n4jhmNgJ4EbjB3f9mZvnA94BGoAz4grs3JOJcIkNKQ0Moijds2LNY7ny+efP+u2dAaEmeNCkskyd3PY8vnEtKkv8zSWQSlbNBeVtE9i2nJIdxl49j3OXj2NW8i5q/xIrsP29j+yvbWfXKKlZ9dRUlJ5VQ8ZEKxl02jpyShJSnkUi3yD8PTIt7/RPgEXe/28zOBn4GfDCKwESSoq0NNm7sKpg7H+Ofb9gATU37P5ZZuOmvZ8Ec/3zSpNAaLZI4ytsi0ifZo7Ipv6ic8ovK6WjtoPaxWjb/ZjNb791K/VNh/OyV169kzAVjmHD1BMrOKBty06ynTWFtZscDG4Ha2OuJwIcI0+8CPAbcY2bT3H11JEGK9FXnDIH7K5i3bOnb8fLyuhfHvRXP48ere4aklPK2iAxU1sgsxpw7hjHnjqG9sZ3qe6rZ/OvN1P2tjuq7q6m+u5r8A/OZ+ImJjP/4+CEzskhaFNZmVgBc6u6fMbOvx1afAmx192YAd283s1XAycDqSAIVgXBj3/4K5g0bQmv0/mRlQUVFV4E8cWL3x87npaUafk7SivK2iCRKTlEOE66awISrJrBj7Q42/XoTGxdspGVlC29/9m1WfXUV4z44jgM+cwAFcwuiDnef0qKwBr4AfLfHuknAth7rGoGJvR2gurqaysquGzTnz5/P/PnzExmjDHft7aGf8r4K5qqqvU+p3VNp6Z6Fcs+CuaICctLlz1BSacGCBSxYsHvEu7FRxjJAg8rbytki0pu8KXlM++o0pnxxCjUP1lB1WxW1j9ay6Y5NbLpjE2MvHMuUL02heF5qR57qa84e1HB7iWBmZwEj3f3+2OvVwFVAJXCZu787bt+XgLvc/Zaex9HQTbJX7lBX13uRHP988+Z9zwbYacSIvRfKnY8TJ4bZA0X6IFXD7SVKIvK2craI9FXzm82s/9F6Nv5iI94a6tbS00uZfvN0So5N/c32SRtuL0E+BxxlXV9zlxFmCftB7Hm8ImBD6kKTtLdjR1fXi721MG/YsO9ZAOONG7fvgnnSpDDxibplSGZT3haRlBl18CgOvu1gpn5tKut/uJ4NP9lA3eN1LHl8CeWXljP9O9PTZkzsdGixLgfip2BbCHwGeB5YBpS7e4uZ5QJNwAx3r+p5HLV+DDMdHeHGvn21MFdVhclO+qKwcP8F8/jxoTVaJMWGYIv1oPO2craIDNTO2p2su2Ud63+wno4dHViOMfnGyUz75jSyC7KTfv60brF29+r412a2C6h297Vm9gBwOvAAcAZwT29FtQwxDQ37v/lv48bQ53l/cnJgwoR9F8waYk4koZS3RSRKuWW5zPj2DCZ+YiKrv76aTb/axLpb1lF9TzUH/+xgRp85OrLYIi+s9+MTwPfM7AhgPKA7W9JZW1uY+W9/N//1ZUxmCF0ueiuU4wvm8vIwsoaIpAvlbRFJibzJecy+YzYTPzGRFdesYPvS7Sx971ImXTeJmbfMJGtk6uuDtCus3X1a3PNa4NroopHddu0KhfGaNbB6dXhcu7Z7wdzXMZnz8/c9UsakSaEVOi8vqT+SiCSG8raIRKn4mGKOXnQ0625Zx+pvrKbqx1XUL6xn7u/npnyq9LQrrCUibW1hauz4wjn++bp1+++akZUV+invb8QMjcksIiIiCZSVm8XUL02l7Iwyln1gGU2Lm3jp2Jc4/KHDKTo6dd1BVVhnio6O0LL89tuwciW88073ArqqKgxLty8VFTB1KkybFh6nTOma/U9jMouIiEjEio8p5uglR7PsA8uofbSWJScv4bA/H0bZqT0HLEoOVUHDSXt76J6xcmVXAd35/O23w9B0e2MWiuTOojm+gO4sovNT+3WKiIiISH/lluZy2AOHsfzq5Wz53y28dv5rHPHkERRXJn9SGRXWQ407VFfDG2+EZflyWLEiFNCrV++7u0Z5OcycCQceCDNmwPTpXYXz5Mkaak5ERESGhawRWcz59Rxw2HLXFl4951WOXnw0eVOSe/+WCut01dERumh0FtCdRfQbb+x77ObJk0Ph3FlAxz8Wp3b6TxEREZGoWJYx+1ez2bl1J7V/rWXZB5dx5N+OJCs3eaOFqLBOB83N8Npr8PLLYXnllbBs3977/sXFMGdOWGbPDstBB4UWaHXXEBEREQHCTY1z/ncOi45cRMOzDay7ZR1TvzQ1aedTYZ1qzc2weDE891x4fOUVePPN0ELd0/jxXQV0/DJhgkbVEBEREemDEWNHMPuXs1n63qWsuXkNFR+pIG9ycrqEqLBOJnd4661QRHcuS5eGMaHj5eTA3LlwxBFw5JHh8YgjQp9oERERERmU0WeOZuzFY9l671bWfmctB992cFLOo8I60dasgSeegMcfD48bN3bfnp0NRx0Fxx4LlZXh+SGHwMiR0cQrIiIikgGm/es0tt67lU2/3MT0b00nd0xuws+hwnqwdu6Ep5+GP/0JHnooDGsXb9w4OOGEUEgfeywcfTQUFEQTq4iIiEiGKjy0kLL3llH711q2/HYLkz45KeHnUGE9EDt3wsMPwx/+AA88ALW1XdtKSuCUU+C00+D000NrtPpDi4iIiESu4sMVobD+vQrr6L3+OtxxB/zP/8CWLV3r58yBCy+E888P3Ts0+6CIiIhI2hl7/lgsx6h/up72hnZyihNbs6kC3B93ePJJ+Pd/h7/+tWv9IYfARz4CF10Es2ZFF5+IiIiI9ElOSQ6F7yqk8YVGGp5vYPSZoxN7/IQebbh54QX49Kfh2WfD64ICuPJKuPpqmDdPXTxEREREhpiS40tCYf2cCuvUqK8PBfUvfxlejxkDN9wAn/wkjE7sByAiIiIiqVNwWBhEonlFc8KPrcK6pxdfhMsvh1WrYMSIUGB/+cuaDlxERERkGMg/KMxS3fJWS8KPrcI63kMPwaWXQktLGF/6rrvCdOEiIiIiMizkzwiF9Y53diT82FkJP+JQ9be/hZE9Wlrg4x+HhQtVVIuIiIgMM7njwsQwO2t24h2e0GOrsAZYvx4uuSSMT33ddfCLX2gmRBEREZFhKCs3i+ySbOiA9rr2xB47oUcbqq6/Hmpq4Kyz4Ic/1GgfIiIiIsNY7thYq/W2nQk9rgrrhQvDdOSFhXD77ZCdHXVEIiIiIpJEOSXhNsP2+mHaYm1mF5vZKjPbZmY/MrOc2PqxZvYzM/t3M/uBmeUm9MQ/+EF4vP56mDw5oYcWERmuIsvZIiIJkJUXSmBvHYZ9rM1sCnAhcCnwaeAa4MbY5t8Dd7j7F4A3gW8l7MSNjfDnP4euH5/8ZMIOKyIynEWWs0VEEqSzsO7Y0ZHQ46bLcHtTgavdvR1YbGaHAqea2TPALHd/Prbfn4EVZnazu28f9FkfewxaW+H442HSpEEfTkQkQ0STs0VEEiRZhXVatFi7+9OxBN2pClgLnAasituvKva0MiEnfu658Hj66Qk5nIhIJogsZ4uIJMiwLqx7cQzwX8AkYFuPbY3AxIScZfHi8FipnC8iMgipydkiIgmSNXJ4dwXZzcwOBKrd/RUzc6DnfJMjgT3GRqmurqYyrkCeP38+8+fP3/fJVq4Mj3PmDCpmEZH+WrBgAQsWLOh8OTbKWAYjpTlbRCRBLDsMrdzXCWL6mrPTqrA2s2zgn4EvxlZVAQfHbTegANjQ873l5eUsWrSo7yfbtQuqYt9SHnDAQEMWERmQ+ELSzLZGHM6ApDRni4gkUmefjT42WPc1Z6dbV5DPAbe4e2vs9aPArLjt04A2YPDZeMsWaG+HsWMhL2/QhxMRyUCpy9kiIglkWf1rse6rtCmszezLwBIg38xmmNnVQB2wMXbHOcDZwK3u3jboE9bWhscxYwZ9KBGRTJPynC0ikkj9bLHuq7ToCmJmXwVu6rF6ubvfYWaXA183s9VAIfCVhJy0sTE8Fhcn5HAiIpkikpwtIpJAyWqxTovC2t1vBm7ey7Y1wD8l/KQNDeGxqCjhhxYRGc4iydkiIomUpBbrtOkKknJNTeFRhbWIiIhIRhn2faxTri3W5W/kyGjjEBEREZHUUot1gu2MDauakxa9YUREREQkRdRinWjtsdl4VViLiIiIZBa1WCdYZ2GdmxttHCIiIiKSUmqxTjR1BRERERHJTBZ7VIt1gqgriIiIiEhmihXW7mqxTgx1BRERERHJTJ0t1omtqzO4sN61KzxmZe4lEBEREclEZp1N1ok9buZWlZ2FdXZ2tHGIiIiISGqpxTrBCgpg4kQoKYk6EhERERFJoZySHEZMGEF2YWIbWDP3zr3rrguLiIiIiGSUqV+ZytSvTE34cTO3xVpEREREJIFUWIuIiIiIJIAKaxERERGRBFBhLSIiIiKSABldWC9YsCDqEIYNXcvE0bVMHF3L4UefaRddi+50PbroWnSXyuuhwloSQtcycXQtE0fXcvjRZ9pF16I7XY8uuhbdqbAWERERERlizD3BU85ExMyqgTX9fNtYYGsSwslEupaJo2uZOEPlWk519/Kog0ilAeZsGDqfaSroWnSn69FF16K7RF+PvebsYVNYi4iIiIhESV1BREREREQSQIW1iIiIiEgCqLAWEREREUmAnKgDiIKZ5QPfAxqBMuAL7t4QbVRDh5ldDPwHUAz8D/BZd283s7HAt4E6IJdwXXdGF+nQYGYjgBeBG9z9b/r9HBwzOxY4EVgJPANsR9dzyNPfRZe95eBoo4pez1wacTiR6pkH3T0jb2Q0sznAdYTrcBCwwN1fTuY5M7KwBn4CPOLud5vZ2cDPgA9GHNOQYGZTgAuBS4G5wE+BdcAtwO+BL7r782b2L8C3gC9FFesQ8nlgWtxr/X4OkJl9HJjp7l+NW/dLdD2HA/1dsN8cnOl65tKM1FsezGB3Ahe4e1Xsb+cvwJxknjDjRgUxs4nAKqDM3ZvNLAeoB+a6++pIgxsCzOw9wMLO1hEz+x4hud8E3OvuE2PrJwErgAp33x5VvOnOzI4HZgNfB64C3kS/nwNiZicDXwPO9Fhi09/78KDPscvecrC7nxttZNHqmUsztcW6tzyYycxsO3C0uy83s3Lglc46JVkysY/1KcBWd28GiCWnVcDJUQY1VLj70z2+cqwC1gKnEa5j535VsaeVKQxvSDGzAuBSd78jbvUp6PdzoH4AvAHcamaPmNkJ6HoOF6egzxHYZw7OWHvJpZmqtzyYyf4A3G5mRcCHgeuTfcJMLKwnAdt6rGsEkvo/mGHsGOC/0HUdiC8A3+2xTtdxAMxsFnAU8HN3vw54HHgEOABdz+FAfxd715mDM1lvuTTj7C0PmllFtJFF6v8AOwh975vc/Z5knzATC2sHWnqsGwnoJrt+MrMDgWp3fwVd134xs7OARe6+pccmXceBmQvUuPvS2OsfA4au53Chz7EXPXJwRtpHLs1Ee8uDF0QXUuTygd8AdwE/NLMzkn3CTLx5sYpwR3m8ImBDBLEMWWaWDfwz8MXYqirg4LjtBhSg67o3nwOOCpcJCL+T9xG+xtPvZ//lEJfP3L3FzN4CstH1HA6Ut3voJQdnql5zqZl9292/F11YkdhbHhwTXUiR+x/gCnevMzMHfmNm05J571cmtlg/CUyMDd2EmeUCU4C/RxrV0PM54BZ3b429fhSYFbd9GtAGLEpxXEPFh4Aj45YNwDXAL9Hv50AsBUp6fOXZTriuup5Dn/L2nnrm4Ey1t1yaid1j9pYHl0UUT6RiQwAf4e51sVU3Aw0keVSQjCus3X0T8ABwemzVGcA9cTfbyX6Y2ZeBJUC+mc0ws6sJY1dvNLNDY7udDdzq7m1RxZnO3L3a3dd3LsAuwle6a9HvZ7+5+3LgYcIQZJhZKTAC+F90PYc85e3uesvBsW4hGWcfuTTjxjjfRx58MMq4IlQD7IiNUtZpG2H0raTJuOH2AMysjDDRwGpgPPAld2+KNKghwsy+ShhaL95yd59jZlMJQx2tBgqBr2jSgr4xs9XEhojS7+fAxFonfkT4lmQK4QaeZbqew4M+x2BfOTiKeNJNfC6NOJRI7C0PRhtVdMzsCMINjIuBCuApd0/qN10ZWViLiIiIiCRaxnUFERERERFJBhXWIiIiIiIJoMJaRERERCQBVFiLiIiIiCSACmsRERERkQRQYS0iIiIikgAqrEVEREREEkCFtYikFTMrN7P7zWyNmX0x6nhERGTflLe7aIIYEUkrZnYW8CRQDCwHDnb3mmijEhGRvVHe7qLCWkTSlpk9BZzm7u1RxyIiIvuX6XlbXUFEJC2Z2QHAXzM1OYuIDDXK2yqsJQnM7P+a2W1Rx9EpEfGY2RlmttDMTklQWMOKmd2S4OMZ8AHgu4k8roh0p3w9eGY2zcx+bmY/S9Dxss3sWjNb04d9zzWz5wf6syYydytvByqshxgzm2Nm95rZ7Wb2ipm5mX046rh6eAj4TdRBxElEPOuAIxIQS1oys/eZ2VYz22hmx5nZAbHfMzezH5tZaWy/eWb2mpn9h5nlx9ZVAJsSHNKHgF+4e3vs+CJDjvL1gAzFfL0l9jgyQcfLAuqAKX3YdxEwZ28bzWyumd0d+937qpmNj9uW6NytvI0K66HoT8C/uvs1wFHArzs3mNnhZnZyFEHFn9vdX3P3Z6KIozf9iWdv19DdV8jK4ZAAAAjTSURBVABbExlXIj+vwR7L3R8CfgK0u/tCd18HXAu0Ac+7e11svxcIifzz7t4Se/t5wAOD+gHimNnXgG8Cj5nZcmBWoo4tkmLK1/00FPO1uzcDGxJ1LnffCSzu476bCUX43ra/DiyIvbzF3eML6YTlbuXtLiqshxAzGw0cDGwHcPcO4OuAm1kJIWlbBHFFdu5E6sPP0ZHCc0VxrLuByWZ2NIC7bwMeB94fd65SYEPsd6/TbHdfPshz7+buN7n7Qe5e6e6z3f2pRB1bJFWUr5MrDfN1okeC6E/8+9t3V+yxZ7/nhOVu5e0uKqyHkNjQNSuAR8zs2Ni6NcALhH5N04H5ZvZxADP7lpl9LfY10O5+VGZWaWY/MLMbzazZzFaY2WfM7EIze9HMLjKzP5pZg5n9c9z7ej1ej3Nfa2a/iu9rFvs69FYz+46ZPWlm82Lr93m+uPebmf2zmb1uZmea2VIz22xmV/XhHAfFx7Ofc+5xDXtRFvtqt9HMru1xTb8Zuzb3xhJxX6/ZHucys0vN7AUzu8rM3jKz7+7jeL199r3Gsy+xlo3XgcviVpcD51qs2wdwPnB/XJz5QHPc69lm9l0zu8XM/m6hS8nlZrbIzC4ws9+YWZ2Z3WBmR8d+xm1m9r79xScylChfZ06+7s3efsZ9ncfMppvZHRZaf7/fl/P0ct6HzOx/zeyE/eyn3J0s7q5lCC3AXGAN4X+odwDj4ratBk6JPT8UeCf2vJzwv+kxhP9xrwUqY9tuA56NPR8B1AA/JoxF+UGgal/H63luQh+znwG/iq3PIxRrpbHXlwLVQNG+ztfLzz0+ds4vx/a9GdgBTN3POXrGs89zxl/DXmJYDfwHMCr2vrWx9cXAb+P2+wchKfbpmu3lXKMJX+/9EJgHnLGv4/X47HuNp4+/X18F3o49nwl8B2gFLo77fbG4/c8H3hN7Pgp4ERgVe/1U7DoUA/XAt4B84PLY6w8Q/nP/WWK/g1q0DKcF5euMyNex7d/s43Xc12f9IjAn7j3ex9+zzs/zOOAbvWw/JXaenLh1yt1JWtRiPcR4aFU8jJBoPgq8bma93aTxFvARM8sF3hNbV0j4Qz4A6Owf+wbhDwZ3bwMagT+4ewPhD238fo7XM75Wut8McS7Q4rE+usC9hD/Ii/dzvp7H7TzmXbF9byb0/z1zP+foFk9/zrkXf/bQn+5FYELczzgm1qJ0I+GaQh+v2V5+3hpCArvP3V9w98f6cby9xdMXdwMzLHQHuQT4OfAX4ANmVgg0eizzxpwAPBt7fj7hH6/OVpCLgJti17kOeMJDv+xFQLG7/87D1+NL6N9nIDIkKF9nRr7uxV5/xn2c5wygxN0741nUz3MeA1zp7v/ax/2Vu5MkJ+oApO/MLBs4wN1XA58ys18Q7qBeALw7fl///+3dT4gcRRTH8e/LkvWQQGIIgezGaA7BgH/IQclFFFdRNMJejAoSMCIkxMgeRERYEMSLkAX/HUQRIohRUQhBMSB68e9BEEVUElhFI4IQRA2i2bjPw6tOaprp6Zmxx81kfh9omJnurqqu3qnUVL2uuP9tZpPAw8BLRRLu/ouZfQbcQPyi3kR0ps6cmr0+TQoXqkqvoqh5GpuB5Vm5Fs3se2BDp/zquPtfZjYPrAbW9pBH33mWnObs92cj8J27P1k+qJs6M7OdRAe2cKnHlLHnZe3hHlSWp467HzOzz4lwkFXuPm9mrxGjSDuAt7NyF38bRfzexcQIU5HWiTzp7PU/tFpEYWlynlF73XJ9o9Be5yrrscO9uZL4EdGvKWCzma1095MVxyyma1DbPUCqkOGynCz+1d2/IKbatpQPNLMpYK+7P+6xwkPufuBqi2WfvgSeqMu4Jr1OfiBGQPNliIyIPfyvxoGjA86jGz/TGoeMmW3roc4OA1uzre3T5T2k17Y83V8OrwK7gPmsfGPADDFtWthGxIsWfgKuMbMVWb4d4/xEzmNqr1uNUntdeY0d8vk9nTNOf/YTo/1PVez/zc8+dK62e4DUsR4+D5jZRPZ+EjiSXp8iHtbYQizttMrMLjCzG9P+dWa2loj1e5GYxvmYNLWYLKP0K93MrCa9ct6WpXEIOAncldJanfIrRj6r8quyPh0zmc490kUeeXnq8syvo6zd6IURo1ArgLfM7GaLh2Tq7sGZvID17v5tti1kWYxlrzull5e7bXnMbNzMnk9hHp28TkxBvwHg7n8A7wAfeetqIDcRYSKFor4PWqyFvQ8o1jIt34PyfR7qFQpEKqi9ZmTa627rsW0+RFs6Tvz4gpidwMwmumi7jRg93gnsMLM70rkPmtksMVr/Xna82u5BGkTgtrbBbMQDEQ78CrwMHCCmki5M+x8FjgPTxNTOMWKEYDsx0nGAaKQOEw9SnErp/Ul80W4BFogHZC4CHkn77+2UXinvGeADIm7tsrTvKmKkc46INbw8fV6ZX8X1ezp/hghNuCLbV5XHJXl56vLM67CU9zTxEN+zqS6K83al/dcBX6V7M5fquds6m25zrben/A4Ck+mzTve0Ja2K8mwg/iODV7r4W3u69P5O4PrSZ/vbnHct8DVwAphNn91K/K09A0yQlhwD9hD/8L6QrvW2pf6OadPW1Iba61FqrzcCHxKj7ltrrrHTvd5OzBR+SjxIfhTYm669bdtNzIosECPVa4D3iU79fcBuYj3vOSIMpThHbfcAN0uVJyPCzNYBe9z9sfR+GTGKstvdZ5e0cDXMzIFNHjGL0icze8jd+1rKSUT+P2qvJae2ezgoFGT07CPCAtZkn02RPZR2LsqmnjTt1CczGzOzu4E3l7osItIVtdeitnvIqGM9ep4jYrG+MbMfgXeB4+7+ydIWq1p6mGMmvb3HYtk36d1K4JC7z9ceKSLnArXXAmq7h4pCQUREREREGqARaxERERGRBqhjLSIiIiLSAHWsRUREREQaoI61iIiIiEgD1LEWEREREWmAOtYiIiIiIg1Qx1pEREREpAHqWIuIiIiINEAdaxERERGRBvwLwbmUvFj2iLAAAAAASUVORK5CYII=\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/ard-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/ard-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/ard-earth.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }