{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 14 - Fire-II - 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 Fire II technology demonstrator for the Apollo heating re-entry environment. See: https://en.wikipedia.org/wiki/Project_FIRE" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \n", "planet.loadAtmosphereModel('../atmdata/Earth/earth-gram-avg.dat', 0 , 1 ,2, 3)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "# Set up the vehicle\n", "vehicle=Vehicle('Fire-II', 86.5, 180.0, 0.0, 0.354, 0.0, 0.805, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,11.35,0.0,-14.70,0.0,0.0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Set up solver\n", "vehicle.setSolverParams(1E-6)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Propogate vehicle entry trajectory\n", "vehicle.propogateEntry (2400.0,0.1,0.0)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "# import rcParams to set figure font type\n", "from matplotlib import rcParams" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtYAAAHtCAYAAAA9RTk1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeXxcddn//9eVNi2lLW2h6ZJu6a4tWyFl33cFFbmLCPgDVAzgLbeCICibIKig+FVR0KDojSCb5RYpyC4gIEsLFNoChe5Nt7SlO9A2uX5/fM4kk2m2SWZyZibv5+NxHjNz5sw510yaa65+cp3PMXdHRERERETapyjuAERERERECoEKaxERERGRDFBhLSIiIiKSASqsRUREREQyQIW1iIiIiEgGqLAWEREREcmArnEHkCn9+/f3srKyVm1bXV1NSUlJdgPKMr2H3KD3kDvy+X3MmDFjtbvnZ/BtlE7OTpbPP+eOos+odfQ5tUyfUeOay9kFU1iXlZUxffr0Vm1bXl7e6m1zld5DbtB7yB35/D7MbFHcMXS0dHJ2snz+OXcUfUato8+pZfqMGtdczlYriIiIiIhIBqiwFhERERHJgE5ZWFdUVMQdQrvpPeQGvYfcUSjvQ5qnn3PL9Bm1jj6nlukzSp+5e8ce0OxY4HrgNHdfGK07BbgZ2AW4C/iuu2+PnusP3ACsA4qBy9x9W+p+y8vLXX1AIpKPzGyGu5fHHUdjlLNFRBpqLmd36Ii1mQ0EegH7Ja0bDpwMTAEuAs4FvpP0sgeAO9z9MmAucF2HBSwi0okpZ4uIpKdDC2t3Xwk8lLJ6BPA1d5/h7ncCvwWOBDCzA4Dx7v5KtO3DwIVm1rOjYhYRaYm7U1NTQ01NTdyhZFQu5Wx32L49LCIiuarDe6zdvTbl8b8Tf0KMVAGLo/tHAQuStq2K7ubkn0xFpHN64IEH6Nq1K6effnrcoWRcruTsl16C4mI4/PD27klEJHty8eTFycDvovtDgDUpz28ESlNfVF1dTXl5ed1SWVmZ5TBFRIKiopBK0xmxrqysrMtXQP/sRNYhOjRnd/BpQSIiQOtzdk5dIMbMxgDV7j4zWuXARymbdQd2OBGmpKREk5iLSCy6dOkCpFdYV1RU1J1xb2arsxJYlnVkzjaLDqDCWkRi0NqcnTMj1mbWBTgPuDxpdRXQL2kbA3oCyzo2OhGRpiVGrGtra1vYsnB0dM5WYS0i+SBnCmvgEuDn7v5J0rongfFJj8uArYCGpkUkZ7RlxLoAdGjOVmEtIvmgwwvraAQDwJLW/QB4A+hhZqPM7GtmNsbdpwPLzWz3aNMTgFvcfWvHRi0i0rRCLqxzJWersBaRfNChPdZm1hv4SvTwHDP7DeFPiT9K2fRdd78jun8acLWZLSTMp3pFR8QqItJahdoKkks5u668FxHJYR1aWLv7RuC2aEm4Plqaes0i4OtZDk1EpM0KdcQ6F3O2RqxFJJflUo+1iEheKtTCOpeoFURE8oEKaxGRdirUVpBcosJaRPKBCmsRkXbSiHX2qbAWkXygwlpEpJ0ShbVGrLNHhbWI5AMV1iIi7dSWS5pLejQriIjkAxXWIiLtpFaQjqMRaxHJZSqsRUTaSa0g2adWEBHJByqsRUTaSa0g2afCWkTygQprEZF2UitI9qmwFpF8oMJaRKSd1AqSfSqsRSQfqLAWEWkntYJkn2YFEZF8oMJaRKSd1ArScTRiLSK5TIW1iEg7qRUk+9QKIiL5QIW1iEg7qRUk+1RYi0g+UGEtItJOagXJPhXWIpIPVFiLiLSTWkGyT4W1iOQDFdYiIu2kVpDs06wgIpIPVFiLiLSTWkE6jkasRSSXqbAWEWkntYJkn1pBRCQfqLAWEWkntYJknwprEckHKqxFRNpJrSDZp8JaRPJB144+oJkdC1wPnObuC6N1PYCbgI1AP+Ayd9/Q0nMiIrmgkFtBciVnq7AWkXzQoSPWZjYQ6AXsl/LUrcAL7v4D4CHg9618TkQkdoXaCpJLOVuzgohIPujQwtrdVxISbR0zKwXOAB6OVj0FfN7Mypp7rkMCFhFphUJtBcnFnK0RaxHJZR3eCuLutdZw6OEIYLW7b4me325mC4DDgW3NPLewI+MWEWlKIbeC5ErOViuIiOSDDi+sGzEEWJOybiNQCmxv5jkRkZyQXHjW1tbWtYYUqFhytgprEckHuVBYO/BRyrruhJEPWniuTnV1NeXl5XWPKyoqqKioyGCYIiJN69KlCzU1NdTU1LSqsK6srKSysjLxsH9Wg8usWHK2CmsRiVNrc3YuFNZVhDPHk/UGlhESeFPPNVBSUsL06dOzEqCISEsShXVr20GSC0kzW53N2DIslpytwlpE4tTanJ0Lf6/8F1AaTdGEmRUDw4HnWnhORCRnFOrMII2IJWdrVhARyQcdXlhbfTOiAbj7CmAacHS0/hhgqrtXNfdcB4YsItKiQp0ZJFdydqK7Zvv29u5JRCR7OrQVxMx6A1+JHp5jZr9x92rgAuAmM9sLGAQkN9o195yISE4oxJlBcilnd+8ebrft0K0tIpI7OrSwdveNwG3Rkrz+Q+AbTbymyedERHJFIbaC5FLO7tYt3H7ySab3LCKSObnQYy0ikvcKtRUkVyRGrLdujTcOEZHmqLAWEcmAQmwFySXFxeF261bNDCIiuUuFtYhIBhRiK0guKSqCrlHzovqsRSRXqbAWEckAtYJkn9pBRCTXqbAWEckAtYJkn05gFJFcp8JaRCQD1AqSfYnCWiPWIpKrVFiLiGSAWkGyL9EK8vHH8cYhItIUFdYiIhmgVpDs69Ur3G7aFG8cIiJNUWEtIpIBagXJvl12CbcbNsQbh4hIU1RYi4hkgFpBsq9Pn3CrwlpEcpUKaxGRDFArSPZpxFpEcp0KaxGRDFArSPapsBaRXKfCWkQkA9QKkn2Jwnr9+njjEBFpigprEZEMUCtI9vXtG27XrYs3DhGRpqiwFhHJALWCZN+AAeF25cp44xARaYoKaxGRDHB3AMws5kgK1+DB4Xb58njjEBFpSkYKazP7bCb2IyKSrxItIImWkFyWrzl70KBwu2JFvHGIiDSlazobm9kJwJXAIOqL8iJgINAjs6GJiOSPRAtIoiUkFxRazk6MWKuwFpFclVZhDfwF+AnwBpA4Q8eAMzIZlIhIvsnREeuCytkDB4JZKKy3bYPi4rgjEhFpKN3Cer67/yJ1pZm9m6F4RETyUi6OWFNgObt7dxg6FJYsgUWLYMyYuCMSEWko3W+AC8xsSiPrP5+JYERE8lWisM6xEeuCy9ljx4bb99+PNw4RkcakW1gfDPzBzGqSllrgtizEJiKSNxKtIDk2Yl1wOTsxSv3BB/HGISLSmHRbQa4DKoDpQGKyVgP+v/YGYmafBr4FfACMBSrd/U0z6wHcBGwE+gGXubsuaCsiOSVHR6wLLmePGxdu33knU3sUEcmcdAvrfwIPufsnySvN7PYMxPIX4AvuXmVmw4HHgU8DtwKPuft90RnuvwdOz8DxREQyJkdPXiy4nD1pUridMSNTexQRyZx0C+t7ge+b2TNJ6ww4Bfh2O2P5NNA7uv8R0MfMSglnr/93tP4pYKqZlbn7wnYeT0QkY3L05MWCy9n77BNuZ87UzCAiknvSLayvIvzJ7xzAo3VGmCO1vUn6b4RewM8AXwEuBI4AVrv7FgB3325mC4DDgYXtPJ6ISMbkaCtIweXsvn3DCYzvvw+zZtWPYIuI5IJ0h1Z+Bgxz9zJ3HxktZYTRj/b6JvAx8Bqwyd2nAkOANSnbbQRKU19cXV1NeXl53VJZWZmBkEREWmfr1q0AFLdyCLWysrIuXwH9sxRWQebsAw8Mt//6V9sCFxFJV2tzdroj1kvcfWMj6weluZ/G9ADuISTmX0ajHE74E2Oy7sC21BeXlJQwffr0DIQhIpK+LVu2ANCzZ89WbV9RUUFFRQUAZrY6S2EVZM4+7ji48054/HG4+OI27UJEJC2tzdnpjlh/M3WFmZ0K/DLN/TTmLmCqu18H/JSQsJcTzipP1htYloHjiYhkhLvXFdY777xzzNE0UJA5+9hjw+3zz8NHqWW8iEiM0i2sdzWzqwDMbGcz+wPwB8KfAtvMzPoDe7n7umjV9cAGYBFQGk3fhJkVA8OB59pzPBGRTNq2bRs1NTUUFxe3uhWkgxRkzh4wACZPho8/hkceyeSeRUTaJ63C2t1PBJ40s1uAN4FxwF5AY1f2Ssda4GMzG5K0bg3wFjANODpadwxhhKSqnccTEcmYzZs3Azk3Wl3QOfuMM8LtXXdles8iIm3XbI+1mR3WxFOLgPeBnwMjgJOBi9oahLvXmtnJwNVmNgMYCFzq7hvM7ALgJjPbi9AXWNHW44iIZEOutIF0ppz95S/Dd78Ljz4KK1fCwIHZOIqISHpaOnnxJsJcpeuA2kaev4MwddNg2pGkAdx9JnBeI+s/BL7Rnn2LiGRTuicuZlGnydmDBsHnPgcPPQS/+Q386EfZPqKISMtaKqyvA6a7+6rmNjKz9v5ZUUQkb61bF1qNe/fu3cKWWdepcvYll4TC+tZb4fLLIf7/14hIZ9dsj7W7P9pSgo62+1vmQhIRyS+rV4eZl0pKSmKNo7Pl7IMPhv33h7Vr4de/jjsaEZH0ZwUREZEU1dXVQPyFdWdjBtdfH+7/+Meh11pEJE4qrEVE2ilRWPfvn60LKEpTjjkGTjoJNm2CK66IOxoR6exUWIuItJNGrOP1s59BcTH88Y/haowiInFpd2FtZpPM7C4z+66ZlbU/JBGR/LJixQoABgwYEHMkLSvEnP2pT8F114X7X/86rFkTbzwi0nllYsT6fOAYd78ZOChxxS0Rkc5iwYIFAJSVlcUbSOsUZM6+5BI48ECoqoLTToPt2+OOSEQ6o0wU1t8C9o3u3wvk1PV8RUSyLVFYjxo1KuZIWqUgc3bXrnDffeFCMU8/DRdfDO5xRyUinU3ahbWZfcbMzozu7weMTVyu1t1r3X1DhmMUEclZW7duZenSpRQVFTF8+PC4w9lBZ8rZw4bB1Kmh3/qWW+AnP4k7IhHpbNIqrM3sBuAB4IsA7v4qcJaZnZGF2EREct7ChQtxd4YOHUpxcW4N/nbGnH3wwfCXv4Sp+K64QvNbi0jHSnfE+iBgKPDvpHX/BH6asYhERPLIW2+9BcDEiRNjjqRRnTJnn3Ya3H57uP/tb8ONN6otREQ6RrqF9ZPuvg5ITlEFcWlcEZG2eOONNwCYNGlSzJE0qtPm7K9/PVzq3Cxc7vzSS6G2Nu6oRKTQpVtYzzezHwDjzezLZjYV+G/gxsyHJiKS+3K8sO7UOfuCC+Cvfw091zffDFOmwIaC6SgXkVyUVmHt7vcCzwBdgLOAdcBx7v7bLMQmIpLT3J0ZM2YAuVlYK2fDl78M06bBLrvA//0f7L8/vPtu3FGJSKHqmu4L3P1l4OXEYzPrYmZHu/vTGY1MRCTHzZkzh1WrVjF48OCcnWpPORuOOw6mT4dTToFZs2DyZPjlL+FrXwutIiIimdJsYW1m/wB6tbCP/sBqoNMkaRERgKefDmnvqKOOwnKgQlPObtrYsfDyy/CNb8A998C558Lf/w6VlTB4cNzRiUihaGnE+q1om/eAWuAbwBvA9KRtRgK7ZiU6EZEcliisjz766JgjqaOc3YyePeHuu+HEE+Fb3wotIp/6FFx/fejH7pr233BFRBpqKY38DLDorHLMbLS7X526kZk9lI3gRERy1caNG3nyyScBOPbYY2OOpo5ydgvM4Mwz4fDD4fzz4ZFH4H/+B+64I7SHHH543BGKSD5r9uRFd1+fSNCRgWbWJXkbMzsC2C8LsYmI5KyHH36Yjz76iIMPPpihQ4fGHQ6gnJ2OoUPh4YdDO8iIEfDmm3DEEWE0O5qaXEQkbelOt/c3YKaZ/c7MfmZm04AnAV3bSkQ6lXvvvReA0047LeZImqWc3Qwz+MIXYM4cuPZa6NULHn0U9t47zCaiAltE0pXudHtPAscRTnwZT5i6aYq7/yQLsYmI5KQlS5bw6KOP0rVrV0499dS4w2mScnbr7LwzXH01zJ8frtTYtSvcdx/stRd87nPw0ktxRygi+SLdEWvcfZm7X+nun3f3rwAPm9m+WYhNRCQn3XbbbdTU1DBlyhQGDRoUdzjNUs5uvZKS0Gc9b14osHv0CCc4Hnww7Lcf3HknfPxx3FGKSC4zd295q8TGZnc0sroEqHH3kzMSkNkBwCHAB8ALwGbgJmAj0A+4zN13uHZWeXm5T58+PXW1iEhGbdiwgZEjR7J27VpefPFFDjrooHbv08xmuHt5BsJL3a9ydjtUV4dC+7bb4MMPw7r+/cNUfeefH3qzRaTzaS5npztiPRZYDCxKWroDM9oVYcTMvgqc5O4/d/e/u/tq4FbgBXf/AfAQ8PtMHEtEpC1uvvlm1q5dyyGHHMKBBx4YdzgtUc5uh5ISuOEGqKoKs4bssw+sXg0//SmMHBkuPHPXXbB5c9yRikiuSHfEery7v5ey7lPAOe5+ebsCMTscuAo41qOgzKwUWAD0c/ctZtYVWA9MdPeFya/Ph9EPEclvq1atYvTo0WzatInnn3+eQw89NCP7zeKItXJ2BrnDK6/Ab38L998PW7eG9T17wpQpcNZZYWaRorSbLEUkn2RsxDo1QUc+Ac5tS2ApfgG8A9xiZo+Z2cHAEcBqd98SHX87IWlrplER6XAXX3wxmzZt4rOf/WzGiupsUs7OLDM44AD4y19gxQr43e/goIPCiPX//i8cfTQMHx7mxX7uOaipiTtiEeloaRXWZrbAzOYnLYsIifX/2hOEmY0HJgG3u/u3CJfafQwYBqxJ2XwjUJq6j+rqasrLy+uWysrK9oQkItLA448/zt13302PHj349a/bP1tdZWVlXb4iXGY845Szs6dfPzjvPHjxRZg7F666CsrKQtvILbeEkevS0rDNE0/Atm1xRywi7dHanJ1uK0gl8Fcg8aJaYKW7z21HrJjZKUClu/ePHvcAqoHrgP9y9/2Ttn0d+Ku7/zx5H/n4Z0URyQ9r1qxh0qRJLFmyhBtvvJHvfe97Gd1/FltBlLM7UG0tvPYaPPggTJ0aZhdJ6NcPTjgBPvtZOP740L8tIvmpuZzd0iXNU93l7s+n7Hy8mY1x9w/aHGGIoy4Wd//IzN4HuhDOKk/WG1jWjmOJiLRabW0tZ511FkuWLGH//ffnoosuijukdChnd6CiIth//7D89KfhAjOJInv2bLjnnrCYweTJ8JnPhEK7vFx92SKFIt1f5RMaWbcSeKCdcbwF9DGzgUnrthOScWk0GoKZFQPDgefaeTwRkVa55pprePTRR9l11125//77KS4ujjukdChnx8QsXGDm2mth1qzQLvKrX4XR6m7d4NVXw3P77w8DB8IZZ8Dtt4dR7jT+kCwiOaZVhbWZfcvM1gGXmVlN8kLop1vZniDc/V3gn8CU6Hh9gW7A3cA04Oho02OAqe5e1Z7jiYi0xu233871119PUVERd999N8OHD487pFZRzs49Y8eGkxofewzWrIGHH4ZvfjP0Za9eHUayKypgzJiw7qtfDRekWbo07shFJB2t7rE2sz2Bs4HUs3a2uvvydgdi1h/4FTCdMMJxu7vPMbN+hIsNLAQGAd93902pry/Efj0Ric99993HmWeeSU1NDb///e+pqKjI2rGy0WOtnJ0f3OG99+Dpp+Ff/wrL2rUNtxk7Fg47LFwB8qCDYNy4MCIuIvFoLmene/JisbvvcG6zmQ1398XtiLHdlKRFJFPuuusuzj77bGpra7n66qu59tprs3q8LJ68qJydZ2prQ2/2v/4FzzwTpu3buLHhNrvtFgrsRKFdXh4uvy4iHaPNJy+a2SXAo+4+J1r1PTPrkroZcCRh/lIRkbzl7txyyy185zvfwd255ppruOaaa+IOq9WUs/NfURHsvXdYLroItm+H11+HF16Al14K0/utWBFaSR5+OLymuDhcFXL//UORve++MH48dEn9yYtI1rU0K8hIoG/S432AnsAK6qdvMsKf+0RE8tb27dv59re/za233grAj3/8Y77//e/HHFXalLMLTNeusN9+Ybn44tA6snBhfZH90kthhPuVV8KS0LMnTJpUX2iXl4cWEs0+IpJdzRbW7v7fKatuAGa5+9bklWY2IdOBiYh0lGXLlnHmmWfy7LPP0r17d/70pz9x+umnxx1W2pSzC58ZjBwZljPPDOs2bAhF9fTpYZkxAxYtCqPcL7xQ/9pevcLI9qRJsOeeYZk4UW0kIpmU7jzWe7j768kroitwjQbmNP4SEZHc9cgjj3DOOeewevVqBg4cyIMPPshBBx0Ud1iZopzdCeyyCxx7bFgSqqtDgT1jRn2xvWQJPP98WBKKisLJkXvtVV9s77lnuDS7TpAUSV+rCmszO5Tw58NjzWxBytO7Es46fzjDsYmIZM26dev43ve+x+233w7Acccdx5133snAgQNbeGXuU86WkpJwpccTkmYyX7UqFNgzZ4blrbfCjCSJ5f7767ft0wf22CMsEyaEke0JE2DAABXcIs1p7Yj1OuBOQv/ewSnPbQX+kMmgRESyxd2ZOnUqF154IStWrKC4uJjrr7+eSy65hKLCaUBVzpYdDBgQrvb4mc/Ur/v4Y3jnnVBkJ4rtmTPD3NqprSQQZiRJLrQTtwMHquAWgVYW1u7+tpkdAZzq7krIIpKX3nrrLS699FKeeOIJAA466CBuv/12JkworJZj5WxprZ12Cj3XkybVr3OHlStDgT17NsyZU3+7Zg38+99hSbbrrg0L7QkTwsmSQ4fqhEnpXFrdY+3u62lilMPMDnH3Fxp7TkQkbsuWLePqq6/mjjvuwN3p06cPP/nJTzjvvPMKaZS6AeVsaSszGDQoLMcfX7/eHaqqGhbas2eHZe3axke4d9op9HCPG7fjsttuGuWWwtPSPNb/AHq1sA8DJgIDMhWUiEgmVFVVcdNNN1FZWcnHH39M165d+eY3v8lVV11F//794w4v45SzJZvMwgj00KFw3HH1691h2bKGhfZ778HcuWHk++23w5KqX7+GhfbYsWG2k1GjVHRL/mppxPo9Qq9eFfVzoKYy4KRMBiUi0h4LFy7kxhtv5I477mDr1jDT3CmnnMJPfvITxo0bF3N0WaWcLR3ODIYMCUvyzCQA69fD+++HIjt1+fDDHeffTujVKxTYo0bVF9uJ+2VlmiJQcldLhfWNwEZ3/6S5jcxMZ5eLSKzcneeee45f//rXPPTQQ9TW1mJmnHrqqVx55ZXsueeecYfYEZSzJaf06RMuTlOecvHnRB/33Lmh8H7vPZg3D+bPD8uGDeFEyrfeany/gwfvWGyPGBGmCRw2DLp3z/pbE2lUSxeIWd3K/XwRnWUuIjHYuHEj99xzD7fccguzZs0CoLi4mDPOOIPvf//7BXdiYnOUsyVfJPdxH3ZYw+fcw2j2/PmwYEF9sZ24v2gRLF8elhdfbHrfI0Y0vfTu3THvUzqfdC8Q04CZDQf+G/gWStIi0kFqa2t59tln+fOf/8zUqVPZsmULAIMGDeL888+noqKCwYMHxxxl7lHOlnxgFmYZ2XXXHUe6AbZvDydRJhfcixbVL0uX1hfeL7/c+DH69Wu+8O7fXz3e0jZpF9Zm1gX4PPAN4FhgJaGnT0Qkq2bPns29997LnXfeyeLFi+vWH3rooZx//vlMmTKFbt26xRhh7lHOlkLTtWt9AXzkkTs+nyi8k4vt5GXx4jAi/uGH8OabjR+jR4/QVpJcbA8dGtpMhg0L99XnLY1pdWFtZiOBc4GvArsBi4DPAY8DB2YlOhHp1Nyd2bNn88ADD/DAAw/wzjvv1D03YsQIzj77bM466yxGjx4dY5S5STlbOqvkwrsx7uEqlE0V3osWhZMuE1ekbEr//g2L7eSie9iwcDKner07nxYLazM7lTDScRTwAfAL4H+BH7j7Y9FmL2UtQhHpVGpqanjllVeYNm0af//73xsU07vuuitf/OIXOf300znyyCMLdg7q9lDOFmmeWbhS5MCBsN9+jW+zfn0Y2U4e5V6yJCxLl4Zl9eqwNDXqDeEYycV26jJ4MBQXZ+d9SjxaM2L9EeESuC8B57r7XAAza2oqJxGRtKxbt47HH3+cadOm8c9//pM1a9bUPZcopk899VSOOuooivUt1BLlbJF26tMH9tgjLI2prQ2zmiQX24n7iWXZsrDNypUwfXrj+ykqCidaJo92py6DB4M63PJHi4W1u08DpkUnvXzDzMqAh4C6bzczM3dX0haRVtm6dSuvvPIKTz31FE8//TQvv/wyNTU1dc+PGjWKk046iZNOOokjjjhCxXQalLNFsq+oKBS8gwc3Peq9fTusWNF88b1iRSjAly1rfD5vqB9hTxTaQ4bsWHwPGaKe71yRziXNFwNXmVlX4GRgqJn9CPgjcBxQmZ0QRSTf1dTU8Pbbb/PMM8/w1FNP8fzzz7N58+a657t06cLhhx9eV0yPHz8e0yn57aKcLRKvrl3rC98DmzirYdu2UFQnCu2qqvpWk8SyfHkowFesaHrkG8LVKhsb8U4uvjXNYPalPSuIu28H/gb8zczGABcR+vmUpEUEgC1btvDKK6/w4osv8sILL/DSSy+xcePGBttMmDCBY445hqOPPprDDz+cPn36xBRtYVPOFsldxcXNn2gJ9SPfycV2agFeVQVr1oRl5sym97XLLs0X30OHQt++mmqwPdo1j7W7fwB828xmZCgeEclDVVVVDQrp119/ne3btzfYZtSoURx22GEcffTRHH300ZpnOgbK2SL5J3nkuym1tVBdveNod+qyYQPMmROWpuy8c+Oj3cmP+/cP7TCyo3YV1gnufmcm9mNm3YDXgG+7+7Nm1gO4CdgI9AMuc/cNmTiWiLTN2rVree211xosy5cvb7BNUVER++yzD4cccgiHHHIIBx98MKWlpTFFLKkylbNBeVskFxQV1c90su++jW/jDmvXNt5ukliWLIFNm8Kl5ufObfp43bo13uudvAwcCF26ZOf95rKMFNYZdClQlvT4VuAxd7/PzE4Afg+cHkdgIp3R2rVrmTlzJm+88Qavvvoqr732GvPnz99hu+NajQkAACAASURBVD59+lBeXs7BBx/MIYccwgEHHEBvNfN1FsrbInnALPRh77Yb7Lln09tt2LBjwZ1ajK9dG654uWBB0/vp0gVKS5svwAtxxpOcKazN7CBgOfBh9LgUOINw+V2Ap4CpZlbm7gtjCVKkQNXW1jJ//nxmzpzJm2++WXe7ZMmSHbbt0aMHkyZNYvLkyXXLmDFjNKd0J6S8LVJ4dtkFJkwIS1O2bGl+5Hvp0nARnsRJmU1JnvFk5EgYPToso0aF26FD82/UOycKazPrCUxx94vN7Opo9RHAanffAuEEHDNbABwOLIwlUJECUF1dzezZs5kzZw6zZs1i5syZvPXWW2zatGmHbXv06MGee+7J3nvvTXl5OZMnT2bixIl07ZoTqUNipLwt0nntvDOMHRuWpnzySZjRpLniu6UZT7p1g7Ky+oI7eRk5MjenGMyVb8fLgJ+mrBsCrElZtxFQo6ZIKyQK6EQRnbitrq5udPvS0lL22msv9t5777rbMWPG0CXfhgukoyhvi0iTuncPRXFZWdPbJM/1PX8+zJvXcFmxovl+7xEjwsj6pz9df/vpT0O/ftl4R60Te2FtZscD0919VcpTTriCWLLuwLbG9lNdXU15eXnd44qKCioqKjIZqkjOqampYdGiRcydO5f33nuP9957r66IXr16daOv6d27NxMmTGDixIlMmDCBvfbai7322ouSkpIOjr5zq6yspLKybsa7/nHGkq5M5G3lbBFpaa7vzZsbL7jnzau/3PyiRfDPfzZ83aBBMGkSTJ5cvwwc2L5YW5uzLe6Lb5nZk8CkpFX9gE3AL4Az3X1c0rbvA9e4+19T91NeXu7Tm5s5XSRPuTurVq1i7ty5dQV04v68efPYunVro69LLaAnTpzIxIkTGTp0qC6+kmPMbIa7l7e8ZW7IRN5WzhaR9ti2LRTYc+bAO++EZc4cePdd+Cj1v/eES8YfeiicdBIcfzzsumvbj91czs6FwrqEMKKR8B/gYuAVYA5Q4u4fmVkxIXGPcveq1P0oSUs+c3eqq6uZP38+8+bNY968eQ2K6A0bmp6tbOjQoYwbN45x48YxduzYuiJaBXT+yMPCut15WzlbRLKhtjbMVjJjBrz2WlhmzAjTCCZ06QInnwzXXdf8SZpNaS5nx94K4u4NGj7NrAaodvfFZjYNOBqYBhwDTG2sqBbJB1u3bmXRokV1xXPy7fz58xs9eTChb9++jB8/vq6ATixjxoyhV69eHfguRJS3RSR3FRXVn+D4pS+FdTU1YST7iSdg2jR4/nmYOhUeewxefhl23z1zx4+9sG7BBcBNZrYXMAhQA57krMSo88KFC1m4cOEOxfOSJUuora1t8vX9+vVj1KhRjB49mlGjRjUopHfbbTeNPku+UN4WkZzSpQtMnBiW886Dykq46KLQw/388wVeWLt7WdL9D4FvxBeNSL3a2lqWL1/OokWLWLhwIYsWLdrh/keNNXZFioqKKCsra1A8J25HjRpFvzhPYxZpB+VtEclFNTWhLWTWLHjlFXjhhdAa8skn4fmxY+tHtTMl5wprkbhs376dpUuX7lAwJ26XLFnS5ImCCf369aOsrIyysrIGRfPo0aMZPnw43QrtElMiIiIx2rYNFi+uvxJkYpk7N5zQmDreZQb77w9nngnnnpv5ubBVWEun4O6sW7eOJUuWsHjxYpYsWVK3LF68mIULF1JVVUVNTU2z+xkwYAAjRoygrKys0VtdxltERKT9amthzZpwEZlly+pvU+9XVYVtmzJ0aGgBmTQJDjkEDjoou/Ncq7CWgrB58+YGxXJy0Zy4v3nz5mb3YWYMGTKkQaGcfH/48OHsvPPOHfSORERECse2baFQrq6G1avD0tz9VavCa1piFornUaPC1RgTy+jRoaDu2zf77y2ZCmvJeVu3bqWqqqrJgnnJkiWsXbu2xf306tWL4cOHM2zYsLol8bisrIxhw4apVUNERKQJtbWwYQN8+OGOy7p1ja//8MNQKK9fn/7x+vWDwYOhtLT+Nvn+4MEwZEi4ymOuUGEtsaqtrWXFihVNFsxLlixhxYoVtDTferdu3RoUzKmF87Bhw+jTp49m1hARkU6ptjbMgrFhQyhyN2yoX5IfJ+6vX79jgbx+PbT18idFRbDbbtC/P5SUNLxt7H5JSeb7nzuCCmvJGnfnww8/bLRoTtyvqqpiWwt/6ykqKqK0tLTJ0eZhw4ZRUlJCUVFRB70zERGRjlFTEwrijRvD0lwx3Nz9DRvaXhQn6907jCQnL337Nr+upCTcdoavaRXW0mZbtmxptmhuTV8zQP/+/RstlhP3S0tL6dpV/1RFRCT31dSEq/wlCuHk+409bmmbLVsyF9vOO0OfPrDLLvW3zd1PLZb79AF9HTdPH480atu2bSxbtqzRGTQS99esWdPiflL7mlPvDx06lB75+LceEREpCNu3p1cIt/S4mcsZtFnPnmGkuHfvHQvg1hTIffqE16oozj59xJ1QbW0t1dXVTY4yL168mBUrVjR7lUAIfc1Dhw5tcqRZfc0iIpJpqSPCLY0Gt/Tcxx9nPsZeveoL4eT7bXncs2fnaKEoFCqsC5C7s3r1ahYsWNDosnjx4hYvdGJmLfY1DxgwQH3NIiLSrNQe4dYWvE2tz/SIsFlmC+Gdd1Yh3JmpsM5TGzZsaLJwXrhwYYu9zbvttluTBfPw4cMpLS2luLi4g96NiIjkCvdQvCZOfGvvqHAme4QTUovZxgrcltYnnlMhLJmkwjqHbdy4kffff5+5c+cyd+7cuvsffPBBi/M29+nTh5EjRza6jBgxgp49e3bQuxARkY6SXBQ3t6xb1/zz27dnNq7kHuHWFrxNrVdrhOQyFdYxq6mpYd68ecyZM6dBET137lxWrFjR5Ot22mmnukK5rKxsh+K5Xzav1ykiIlnjHkZ+164NV6pbu7bxJTGvcGqRnImiuHv3+pPf2lMIq0dYOhsV1h3E3Vm2bBmzZs3i7bffrrudM2cOHzdx5kT37t0ZM2YMY8eOZdy4cXXLmDFjGDRokE4KFBHJce6h2F21ClauDLerVoWCubmiuT3FcaIoTnfp27f+fi5dyU4kn6iwzgJ3Z+nSpbz66qu8+uqrvPbaa7z55pt8+OGHjW4/dOhQJk6cyPjx4+uK57FjxzJs2DC6dOnSwdGLiEhLtmyBpUvDUlVVXzAniufkIrqFc8Ub1bMn7Lpr48tuu4Xbvn0bFsMqikXip8I6A7Zv387rr7/OM888w0svvcSrr77KypUrd9iuX79+7LHHHuy+++51t7vvvjt9+/aNIWoREWmMOyxfDu+/DwsXhuJ5yZKGty2c5tJAr14wYAAMHBhuBwwIV6JrqnDu1w922ilrb09EskiFdRstWLCAf/zjHzzzzDM8++yzbNiwocHz/fr1Y/Lkyey3337st99+7LPPPpSWlqp9Q0QkR3zyCcyaBW+9FYroxPLBB2F6uOYUF8PQoWEZMqS+aE69HTAgzDohIp2DCus0zJs3j3vuuYcHH3yQN954o8FzY8aM4aijjuKwww5j//33Z/To0SqiRURyRG1tKKCfew5mzIA334R33mm6l3m33WDsWBg1CoYNqy+iE/dLSnRCnojsSIV1C2pra/nHP/7BbbfdxhNPPFG3vlevXpx44omccMIJHHXUUQwfPjzGKEVEJNXatfDgg/DII6GgTj3NxQw+9SnYe28YPz4U0olFEyuJSFuosG6CuzNt2jR+8IMfMGvWLCBMcXfqqafypS99iWOOOYad1AQnIpJT3OGJJ+DXvw63ySPSw4fDEUfAgQeGYnqPPcJJgiIimaLCuhGrVq2ioqKChx56CIBhw4Zx0UUXcfbZZ7PrrrvGHJ2IiDTmP/+B//5vSHTqdekCxx8P//VfcMwxMHJkvPGJSOFTYZ3irbfe4jOf+QzLli1jl1124dprr+WCCy6gu+YvEhHJSe5wzTXwox+Fx4MGwbe/DeeeC/37xxubiHQuOVNYm9kpwM3ALsBdwHfdfbuZ9QduANYBxcBl7r4tGzHMnj2bww47jPXr13PIIYdw9913q3daRKQRuZCzEy66CH71q3Ay4eWXw5VXQo8e2TyiiEjjcuKcZjMbDpwMTAEuAs4FvhM9/QBwh7tfBswFrstGDBs2bOBzn/sc69ev5wtf+AJPPvmkimoRkUbkQs5OePDBUFQXF8PDD8MNN6ioFpH45MqI9Qjga+6+HZhhZrsDR5rZC8B4d38l2u5h4D0zu97dW5hlND0//OEPWbBgAZMmTeKee+7RiYkiIk2LPWdDaAG5/PJw/+ab4bOfzfQRRETSkxMj1u7+7yhBJ1QBi4GjgAVJ21VFd8szefwVK1bwm9/8BjPjD3/4Az003CEi0qS4c3bCq6+GC7oMGQIXXJCNI4iIpCcnCutGTAZ+BwwB1qQ8txEoTX1BdXU15eXldUtlZWWrD1ZZWcm2bds4+eST2WeffdoTt4hIq1RWVtblKyDfT7Hr0Jyd8MIL4fbEE6Frrvz9VUQKUmtzds6lIjMbA1S7+0wzc+CjlE26AzucCFNSUsL06dPbdMz7778fgPPPP79NrxcRSVdFRQUVFRUAmNnqmMNpszhydsK774bbPfZo125ERFrU2pydUyPWZtYFOA+IuuaoAvolPW9AT2BZpo75wQcfMHv2bPr06cORRx6Zqd2KiBS8OHJ2stXRV1vpDuPhIiLxyKnCGrgE+Lm7fxI9fhIYn/R8GbAVaN8wR5JnnnkGgOOPP57i4uJM7VZEpDPo8JydbO3acKvLj4tIrsiZwtrMfgC8AfQws1Fm9jXCPKjLozPOAU4AbnH3rZk67gtRk96hhx6aqV2KiBS8uHJ2svXrw22fPtnYu4hI+nKix9rMrgR+lLL6XXe/w8xOA642s4VAL+CKTB47UVgfcsghmdytiEjBijNnJ9scTeDXu3e2jiAikp6cKKzd/Xrg+iaeWwR8PRvHXbt2LQsWLGDnnXdmD539IiLSKnHl7FSbNoXbnj074mgiIi3LmVaQOMyaNQuAiRMn0qVLl5ijERGRdCQK61694o1DRCShUxfWb7/9NoBGq0VE8ox7fSuIRqxFJFeosEaFtYhIvvnoo1Bc77QT6A+OIpIrOnVh/d577wEwYcKEmCMREZF0JEar1QYiIrmkUxfWCxYsAGDkyJExRyIiIunQiYsikos6bWG9fft2li5dipkxfPjwuMMREZE0aMRaRHJRpy2slyxZQk1NDaWlpXTv3j3ucEREJA0asRaRXNRpC+uFCxcCUFZWFmscIiKSPo1Yi0gu6rSFdVVVFQDDhg2LORIREUmXRqxFJBd12sJ6xYoVAAwaNCjmSEREJF26OIyI5KJOW1iXlZXx+c9/nn333TfuUEREJE2lpXDyyTB5ctyRiIjU6xp3AHGZMmUKU6ZMiTsMERFpgyOPDIuISC7ptCPWIiIiIiKZpMJaRERERCQDVFiLiIiIiGSACmsRERERkQzolIV1ZWVl3CG0m95DbtB7yB2F8j6kefo5t0yfUevoc2qZPqP0qbDOU3oPuUHvIXcUyvuQ5unn3DJ9Rq2jz6ll+ozS1ykLaxERERGRTDN3jzuGjDCzamBRKzfvD6zOYjgdQe8hN+g95I58fh8j3L0k7iA6Upo5O1k+/5w7ij6j1tHn1DJ9Ro1rMmcXTGEtIiIiIhIntYKIiIiIiGSACmsRERERkQzoVIW1mfUws1vM7MdmdpuZ7RJ3TOkys1PMbIGZrTGzX5lZ17hjaisz62ZmM83siLhjaQ8zO8DMLjGzk82sf9zxpMPMPm1mvzWzi8zsVjPbO+6YWsPMjjWzV8ysLGld3v9+S9P0821cU98JZtbfzH5vZjea2S/MrDjuWOOW+p2jf1ONS/1O0+eUnrwtytroVuAxd7/PzE4Afg+cHnNMrWZmw4GTgSnAROA2YAnw8zjjaodLgbK4g2gPM/sqMNrdr4w7ljb6C/AFd6+K/n09Dnw65piaZWYDgV7AfilP5fXvt7RIP98ULXwnPABc7u6vmNn5wHXA9+OKNUekfufo31SKxr7TzOxP6HNqtU5z8qKZlQILgH7uviX6X/16YKK7L4w1uFYys0OB/7j79ujxTYT4T4w3svSZ2UHAp4CrgXPc/dl4I0qfmR0OXAUc63n6i2Rmm4F93f1dMysBZrp7adxxtcTMioAaYKS7LyyE329pmn6+jWvqOwH4EfBg4nfZzIYA7wED3X1zXPHGKfU7B5iL/k010Nh3mn730teZWkGOAFa7+xaAKBEtAA6PM6h0uPu/Ewk0UgUsjiuetjKznsAUd78j7lja6RfAO8AtZvaYmR0cd0Bt8DfgD2bWG/gKcGHM8bSKu9emrDqCPP/9lmYdgX6+O2jmO+EowueT2K4qulvegeHljCa+c45A/6ZSNfaddgT6nNLSmQrrIcCalHUbgZwfnWvGZOB3cQfRBpcBP407iPYws/HAJOB2d/8W8DTwWNSmkE++CXwMvAZscvepMcfTVoX4+y319PNtncR3gj6vhhr7ztFnlKSp7zRgGPqc0tKZCmsHPkpZ1x3YFkMs7WZmY4Bqd58ZdyzpMLPjgenuviruWNppIrDW3d+KHv8GMOAL8YXUJj2Ae4C/Ar80s2NijqetCur3W3agn28LUr4T9HlFmvnO0WfUUFPfafqc0tSZTl6sAvqlrOsNLIshlnYxsy7AecDlccfSBpcAk8ws8bgf8JCZ3eDuN8UXVtq6kvT74+4fmdn7wG7xhdQmdwFfdvd1ZubAPWZWlod9mAXz+y2N0s+3GY18J1QB45KeN6AnnfPzavQ7h9D2oH9T9Zr6TuuCPqe0dKYR638BpWbWAyCaemg48FysUbXNJcDP3f2TuANpgzOAvZOWZcC55F9Ly1tAn5TWj+3AnJjiSVs0NeBe7r4uWnU9sIEcnxWkCYX0+y070s+3eanfCU8C45OeLwO2AtM7OK5c0NR3zp/Qv6lkTX2nLUOfU1o6TWHt7iuAacDR0apjgKlJJ3XkBTP7AfAG0MPMRpnZ16I/AeYFd69296WJhTCzQ7W7b4g7tnS4+7vAPwnTXGFmfYFuwCNxxpWmtcDH0YwBCWsIZ8vnNKsffjIonN9vaZx+vk1r7DsBWAcsN7Pdo81OAG5x961xxRmXZr5zFqN/U3Wa+U67G31Oaek00+0BmFk/4CZgITAI+L67b4o1qDSY2ZWEaZSSvevu+TjCCICZLSR/p9vrD/yKMAo0nHDSR96MWAOY2V6EExhnAAOB5909p0cikmYwuZUwN+9v3L0633+/pXn6+e6oue8EMxtBmFpuIWHe9ytSZhDplJK/c/RvqqGmvtP0OaWnUxXWIiIiIiLZ0mlaQUREREREskmFtYiIiIhIBqiwFhERERHJABXWIiIiIiIZoMJaRERERCQDVFiLiIiIiGSACmuRJpjZ0Wb2kpkd0cJ2u5vZ/WZ2VQaP/bNM7UtEJA5mtqeZPZDJ3NjEcaaZ2ZezeYxsMbNjzOw/LX3PSP7o2vImIvEys08DNxCuFDgZ2BP4/9z9riwfegmwVyu2Wwz0A7pk4qDRJWO3ZWJfIiIAZjYZ+AlwGPAEsIBwsY9i4DZ3fzwLh/0A6EOGcmMzbgVmZWPHUT4+Bfgd4aq0/4me6g78F7Cfuy9sxyFa+z0jeUKFteSDvwNfcveZZlYE/KkjDuruc81sTSu222BmyzN46MOB5zO4PxHp5Nz9NTO7F/iUu5+UWG9mJwL/Z2aXuPutGT7mlgznRgDM7HPA24mC1t0fzfQxEtx9G3CfmV0IPOvuVybF8VBb9mlmA4Aj3f0+d3/PzFZnKFzJAWoFkZxmZrsC44DNAO5eS7hMb0ddMrQ2w9u1xpHAvzK4PxERgB0uae7ujwA/BP6fmQ3KwjEzmqvNbChhhLqjNXY5+KeAqnR2YmbdgTuBHkmrM/n9ITFTYS05zd3XAu8Bj5nZAdG6RcCrZtbXzH5kZq+a2efNbIWZvWNm+yReb2anm9kVZvasmf04aX25mf3QzO4zswfNrE+0ficz+62ZXW1mfyT8GTMtZva9KKYvmdlkM/t7tL8fRTHeZ2Yl0fr1ZnZNyi66ufsn0b6+a2YXm9mLZvbTtD9AEZGW/ZHQEnIqgJmNNrPrzOwPZvZ0VMxiZt2jXHaFmf3bzE5L7MDMvmlmV5nZy2b2raYO1FhONrMpUc48x8zeN7Ofmlmpmf0pyqf/NrOjol2cAQwFLjOzk6NzYZ4ys7OTjvENM7vBzO40s/81s55m1sfMrjGz6WZ2lJm9Z2bzzaysLR9Y9JkcEY1oY2anmNmvzOxmM3vMzHpF7+FXZvZw9H2zPPqMJwCnmtnFSbvsF30XbTSzbzRxzCIzu9HMKqJjuJn9X+LnIznC3bVoyekFmAgsIvyv/g5gQLS+CyFJrSX0uvUFngRmE/7TeCBwZbRtP8Ko94nALsC9Sft/EfhZdP//ARXR/YFADSF5thTjnwmjPgMIvXg7R+u7A48DDwGlhC+EbcA1wM7AvtHjxPaJ/vHE+741ut8buDzun4UWLVrydwHOAZY28dxK4NdRXn0Q6BKt/yvwQHT/F8Ap0f2vAdXR/dOBr0T3x0W5eo/o8Z+BH0b3m8rJuwLrgF8C+wHHAD8Hro62vRR4OClWB8qi+z2Ad4BzosfHA39N2nYa8HtC6+sXo++LU6L3+RhwQys/u2ej74qfAj8j9Fsfk/T8amB4dH824TupO3AFsBTYH/h6lPefTcQbbb8QuDl67nRgcRMxnAtMi+4PjD7n4+L+d6Wl4aIea8l57j7bzPYArge+CXzOzI7x0HNdDWxw96kAZvZD4AVgFOFLpNbMvhPt6nFCAj8R2C1p/TvRa3sBF0Svxd1XmtmSNELtD9wCfN3dt0T7+MTMVgAL3H1ZdJwVwHMe+g/fICT8/oSTIE8kfAkAbAHOMrO5wG2E/1SIiGSDE4rN/YARwIVmBrAG6GPhwXmEQQGA/yW0QkDItdNTcu1g4O2UY5xDIznZ3dea2XrgIXd/FcDMFgObzWwwYYS3V6NBu39kZquSVn0r2m/CH4D7o/XrCN8XD0bHmB7F2Vr/8qjH2szuIwzmJJzq7ovN7GBCQd0ryv/LgHnu/grwSvTaxvb9cPSd8FozMe0DfBS975VmtpaGLSWSA1RYS04zsy7AMA8nqfxP1J7xKFBJGAFINTu67QsMB/7i7n+N1v0y2udlhEL3lynH2peQEDe2MdxxwCHAaGBm0vrUHsOauifca6Mkm2jL2s1D+wvuvsDMzgF+A1wEnAmsQkQkg8ysL+GvbbMJeXNNI/lxIGFEtRuAu9cQBgOIXnOtu78UPW7w2iSN5uSI0zBXLgUuA94HXgZGNvMWkl83Fngm6fE8QpvLAHbMxdtp+4wlrxP+A5KwzsxuBB4gxJ6onlPfV0u203Rt9k/gNjPbifCz+Igw+i05RD3Wkuvq+v4A3H0m8APgU01s342QxD4AlgNTkp80s/2j9SeaWY+U9Ruih59uY6wvAb8C7o4SX1rMrJSkE2HMbAhhRpRxhJMZ/9rES0VE2uNswl/IHiDkx4OjkWKgLj+uBrYCJySt38PMepOSa82sazRQkaqpnNyY3wELPUyr+kka72UxYYS77hDAekKrS8a4e200+DEx+o/Jo8B17j49k8dJOebDwMOEv6yeDhzu7uuzdTxpGxXWkg8ujIrOhCGE3riEXc2sW3T/CEJ/3TrgHuCLZvYbMzsiGk3YTEiAPYFpZnZ8dKLIpwjF+BvAj8ysm5n1J/Q2DzCzLmZ2uZl9pYkYLVquJnz5/KKR51K3T318EqEfMGEscLK7bwAubOQ1IiLp2GEk1MwOA64FznP3asLo8ErgUTP7XHSC4qHRCPW9hNlDTjazo4Gz3X0jIdf+j5lda2aHEwYYEm10yfmvqZyckDx6PAkoMbOdgYOBHmY2KnpuG+Fkv8QAS/Ixfgf8l5n1ix7vD9zu7tsJNU+judjMDoji6t7EZ1dESs0U/efjRMJfKfsRWgzHA8OAnZNOjEwdFd/aSPwNg2qkXyT6WewCPE0YlU/nPxzSQdQKIvlgGDDbzB4h/JnMgPOTnt8OXGFm6wijuxcCuPuTZvY/hD8nnkw4+W8WgJl9Afgt4YviDuBSd3cz+xJhKqTZwH2EP+ntTbigwpcILSYNLkxjZvsRTsoZTvjT4N+Ba6N4HgcOAMrMbDRhJGUQ8CUze5f60ZtzgJ7uXpny3v9k4cIO2wknC4mIpC3KI2cQir/bgWrCCXB9gM+4+38A3H1rlB9vJ+S6hwi91QDfJgxK3En4C11iJo4/EnLfBYSWtQvcfVVUOB4AjDazu5rKyWY2hdBXfK6ZvevuVYQTKX9OKKrvBL4A7AHMB/5CGF0/y8JsUROBz5rZo+7+oJmNBaaa2YuEv3peEY2snwoMit7fbMLFcgab2cQozi8T+sb/nvS5FRNy/17AkKjwdkJv8ynRPt8EngNeI7S3/JOQ06cRTpjcy8xOdvfEfu8jXKxnU1SIDwKmmNkCwkg00etTr9mwmHBy5snATkCRmT1FOIGxo6aglRaYfhaSzyxcBvbP7l7WQce71N11uXERkQITtaV0d/ecvECXmV1KmKFlYfR4J8JfG25MnJsj8VMriOS7xtossnOg8CfT1zviWCIi0nGikzP3yeGieifgSqAsaXV/YLOK6tyiVhDJW1F/21cIf9o71d0fyPIh50V/ohQRkcLS1d1vizuIprj7xxYuq/7r6MT2ZYS2lSubf6V0NLWCiIiIiIhkgFpBREREREQyQIW1iIiIiEgGqLAWEREREcmAgjl5sX///l5WVpb266qrqykpKcl8QB2oEN4DFMb7KIT3AIXxPvLpPcyYMWO1u+dHsBnSWXJ2vsULirmj5FvM+RYvZC/m5nJ2wRTWZWVlTJ+e/pVEy8vL2/S6XFII7wEK430UwnuAwngf+fQezGxR3DF0tM6SpcuYfQAAIABJREFUs/MtXlDMHSXfYs63eCF7MTeXs9UKIiIiIiKSASqsRUREREQyoNMX1hUVFXGH0G6F8B6gMN5HIbwHKIz3UQjvQXaUbz/XfIsXFHNHybeY8y1eiCfmDr9AjJkdC1wPnJZ0vftTgJuBXYC7gO+6+/bouf7ADcA6oBi4zN23pe63vLzc8633R0QEwMxmuHt53HE0RjlbRKSh5nJ2h45Ym9lAoBewX9K64cDJwBTgIuBc4DtJL3sAuMPdLwPmAtd1WMAiIp2YcraISHo6tLB295XAQymrRwBfc/cZ7n4n8FvgSAAzOwAY7+6vRNs+DFxoZj07KmYRkc5KOVtEJD0d3mPt7rUpj/+d+BNipApYHN0/CliQtG1VdDczfzK98UaYNAnuuy8juxMRKTS5krM3vr6R6ftO592vvtveXYmIZE0uzmM9GfhZdH8IsCbl+Y1AaeqLqqurKS+vz90VFRUtN60vXw5vvglVVc1vJyKSYZWVlVRWViYe9o8zlnbqkJxd+3Etm17fRFG3Tn/OvYjEoLU5O6cKazMbA1S7+8xolQMfpWzWHdjhRJiSkpL0JwHvGf11cvPmNCMVEWmf5ELSzFbHHE6bdGTOLuoeCurarbUtbCkiknmtzdk5819/M+sCnAdcnrS6CuiXtI0BPYFlGTloorDesiUjuxMR6Sw6OmdbNwPAt3bsTFYiIunImcIauAT4ubt/krTuSWB80uMyYCuQmTmaNGItItJWHZqzEy0gtZ/8/+zdeZwcVbn/8c8z+74kM9l3EkgIgQBhkR0CCoioCOjFDbwQXC8ueJXFDbkKiv5QXIMiohcuXFDZlFUuBMFAYiCQBQJJJhCyTCazZWaS2c7vj1M909OZnrV7unv6+3696tXdVdVVT1VXn3nm9KlzVGMtIslrxBProAYDwMLmXQ2sAvLNbJaZfcrMZjvnVgDbzOyQYNUzgVucc60xCUaJtYhIn5KlzA7VWKspiIgksxFtY21mxcDHgpcXm9nP8D8lfjdi1fXOuduC5x8Gvmlmm/H9qV4Ts4CUWIuIRJVMZXaoxlpNQUQkmY1oYu2cawR+GUwh1wdTtPdUAf8el4CUWIuIRJVMZbblqsZaRJJfMrWxHnlKrEVEUoJqrEUkFSixBiXWIiJJrquNtW5eFJEklt6JdUGBf1RiLSKS1DKyu2usnVOttYgkp/ROrEtL/WN9fWLjEBGRPlmmQaZ/7tqVWItIckrvxHrMGP+4e3di4xARkX6FRl9UO2sRSVbpnVgXFUF2th95ce/eREcjIiJ96BokRj2DiEiSSu/E2gzGjvXPa2oSG4uIiPRJNzCKSLJL78Qa1BxERCRFqMs9EUl2SqxVYy0ikhI0rLmIJDsl1qHEeteuxMYhIiJ9ysgP2li3KLEWkeSkxHrSJP+4dWti4xARkT5lFvr+9jqaOhIciYhI75RYT5niH99+O7FxiIhInzKLgsR6jxJrEUlOSqyVWIuIpIRQjXVnk5qCiEhyUmKtxFpEJCWoxlpEkp0S61Bi/dZbiY1DRET6pDbWIpLslFhPneoHinn7bWhtTXQ0IiIShWqsRSTZKbHOy4OZM6GjA954I9HRiIhIFF2JtWqsRSRJKbEGmDvXP65bl9g4REQkqoxC/ydLNdYikqyUWAPMm+cf169PbBwiIhKVmoKISLJTYg3difWaNYmNQ0REosoqzgKgo0GJtYgkJyXWAAsX+seVKxMbh4iIRJU11ifWbTVtCY5ERKR3SqwBFiyAnBx4/XWoq0t0NCIi0ovssdmAEmsRSV4jnlib2RlmttzMZoTNyzezW8zse2b2SzMrGciymMnJgcMP989XrIj55kVEUlUyldlKrEUk2Y1oYm1m44Ei4OiIRb8AnnXOXQ3cD/x6gMti56ij/OPy5XHZvIhIqkm2Mrsrsd6lxFpEktOIJtbOuR34graLmU0CLgIeDGY9AZxrZjP6Whbz4E44wT8+9VTMNy0ikoqSrczOKsuCDOio76CzvTMWmxQRiakRbwrinIssDU8BdjnnmoPl7cAm4OR+lsXWaaf5x2efhZaWmG9eRCQVJVOZbZlGVrm/gbF9d3ssNikiElPJcPPiZKAmYl4jMKmfZT1UV1ezaNGirmnp0qWDi6Ky0rez3rcPli0b3HtFRAZp6dKlXeUVUJHoeAYhoWV2dkXQHKRazUFEZOQMtMzOGrmQonJAZBVxLhAqNfta1qWyspIVw73x8IwzYNUqeOwxePe7h7ctEZE+LFmyhCVLlgBgZrsSHM5gJLTMzp2US8trLezbuo/C+YWDfr+IyFAMtMxOhhrrrUB5xLxi4J1+lsXe2Wf7xz//GZyLyy5ERFJcQsvs3Cm5AOx7e1+sNikiEjPJkFg/BUwys3wAM8sGpgFP97Ms9k44AcaPh40bfc21iIhESmiZnTs1SKzfUmItIsknEf1YW+gpgHNuO/AQsDiYfzpwn3Nua1/L4hJcZiacf75/fs89cdmFiEgqSbYyuyuxVo21iCShke7Huhj4dPDyYjOrDJ5/Bni/mV0DnA0sCXtbX8ti74IL/ONdd0FHR1x3JSKSzJKxzFZTEBFJZiN686JzrhH4ZTCFz68FLovynqjL4uLEE2HmTNi0CR59tLvdtYhImknGMjtvWh4Ae6v2xmsXIiJDlgxtrJNLRgZcfrl//stf9r2uiIiMqPwD8gFoebNFg8SISNJRYt2bT30KcnLg4YfhzTcTHY2IiAQyCzPJnZqLa3Xs3axaaxFJLkqse1NZCRdd5Lvcu/HGREcjIiJhCg4qAKDlNY2SKyLJRYl1NFddBWZw++2wZUuioxERkUDBXJ9YN69vTnAkIiI9KbGO5sAD4cMfhrY2+P73Ex2NiIgEQol105qmBEciItKTEuu+fOMb/mbGW2+FNWsSHY2IiABFhxcB0PivxgRHIiLSkxLrvhx8sO8hpKMDvvQlDXMuIpIEihYWQQY0vdpER4vGGxCR5KHEuj/XXQdlZfD443DffYmORkQk7WUWZFI4vxA6YM/LexIdjohIFyXW/amogO99zz///Odh9+7ExiMiIhQvKgag8UU1BxGR5KHEeiAuvxxOOgl27PBNQkREJKFKjikBoH5ZfYIjERHppsR6IDIy4De/gbw8uOMOuP/+REckIpLWyk4rA6D277W4Tt3/IiLJQYn1QM2Z090k5JJLoKoqsfGIiKSx/Nn55E7Jpb2mnaZX1O2eiCQHJdaDccUV8N73Qm0tfOQjvo9rEREZcWZG2eKg1vrJ2gRHIyLiKbEejIwM+P3vYcoU+Oc/4XOfUxd8IiIJMuaMMQDUPFST4EhERDwl1oM1dizce69vb33rrfCDHyQ6IhGRtDTmvWOwbKPu6Tpaq1sTHY6IiBLrITnmGPjjH8EMvv51uOuuREckIpJ2ssuyKV9cDp1Q84BqrUUk8WKSWJvZ2bHYTkr50Ie6a6s//nENHiMiKWM0ldkV51UAsPOenQmORERkkIm1mZ1pZs+a2RtmtjGYNgPpmVV+5Stw1VV+yPOPfAT+8pdERyQi0iUdyuzK8yqxHKP28Vr2btmb6HBEJM1lDXL9PwDfB1YBncE8Ay6KZVApwwz+67+gvR1++EM4/3zf7vqSSxIdmYgIpEGZnT02m8rzKtn5PzvZ/rvtzPjWjESHJCJpbLCJ9Ubn3I8jZ5rZ+hjFk3rM4MYbISfHJ9mf+hRs3+7bXpslOjoRSW9pUWZPvHQiO/9nJ9tu28b0a6djmSp7RSQxBtvG+jNmdn4v88+NRTApywyuvx5uucU/v/pqn2Dv1c+SIpJQaVFml51aRt4Beezbso9df9mV6HBEJI0NNrE+HviNmXWETZ3AL4cbiJnNM7Ofm9mXzOwXZrYwmJ9vZreY2ffM7JdmVjLcfcXN5z8P99wD+flw++1w4omwZUuioxKR9JUWZbZlGFO/PBWALTdswWl8ARFJkMEm1tcBS4A5wKyw6dsxiOUPwPecc/8PuAEI9WH3C+BZ59zVwP3Ar2Owr/g5/3x4/nmYORNWrIAjj4T77090VCKSntKmzJ5wyQSyK7NpXNGokRhFJGEGm1j/DbjfObfROVcVTJuBW2MQyzygOHjeApSa2ST8TTYPBvOfAM41sxkx2F/8HHaYT6rPPBN27YIPfAAuvRQaGxMdmYikl7QpszPzM5nyxSkAbP7mZtVai0hCDDax/h/gKjM7KWw6GbgqBrHci//Jshj4GPAF4BRgl3OuGcA51w5sAk6Owf7ia8wYePhhuPlmyM2F3/4WFi6Exx9PdGQikj7Sqsye/IXJZI/LpuH5Bqrvqx6JXYqI9DDYxPobwBeBO4Dfh02XxyCWzwJ7gReBPc65+4DJQORwWo3ApBjsL/4yMuCKK2DlSp9Ub9wI73637/P6nXcSHZ2IjH5pVWZnFWcx4zszANj4tY107O0Yid2KiHQZbGL9Q2Cqc26Gc25mMM0AzotBLPn4Nnp3Ajeb2emAw//EGC4XaIt8c3V1NYsWLeqali5dGoOQYmT+fFi+HL7/fX9j4913w9y58KMfqecQkTS1dOnSrvIKqIjTbtKuzJ546UQKDi5g78a9VF1fFZNtiogMtMy2wbRDM7N3Oeee72X+p5xztw0p0u5tPAJ8xDlXZ2bfAP4DX9PyLefcgWHrbQjm3Rn+/kWLFrkVK1YMJ4SRUVXla7FDNzROmwbf+Y4fFj0zM7GxiUhCmNlK59yiOGw3Lcvs+n/Us+rEVVimceTKIyk6tCgu+xGR9NRXmT3YGuvP9rLxC4CbhxJY2DYqgMOcc3XBrOuBBqAKmGRm+cF62cA04Onh7C+hpk/3Q58//DAccojvju+SS+DQQ+FPf4LOzv63ISIyMGlZZpceX8qkz07CtTvWX7xeTUJEZMQMNrEeE9RMYGYFZvYb4Df4NnbDsRvYa2aTw+bVAKuBh4DFwbzTgfucc1uHub/EO/tseOkl+MMfYMYMWLsWPvQhWLAA7rgD2vb75VREZLDStsye9b1Z5M3KY8+qPbx55ZsjvXsRSVODagoCYGbHAh8F3gNsBz4B1DvnhtVxqJkdhq9dWQmMB55xzj1tZuXAD4DNwATgKufcnsj3p0xTkN60tsKtt8INN8Dbb/t506bBF7/oa7PLyhIbn4jEVbyaggTbTtsyu2FFA6uOW4Vrcxx8z8GMu2BcXPcnIumhrzK7z8TazE6Ksuho4FTgJqAT+IBz7kvDDXQ4UjqxDmlthTvvhBtvhPXr/bz8fPjoR+Fzn/M9i4jIqBOrxFpl9v7e/unbvHHFG2QUZHD4M4dTfGRx/28SEenDcBLrf+IHAajDF8a9rgZMdM7lDjfQ4RgViXVIZyc8+CDccgs8+WT3/OOOg89+Fs47zyfcIjIqxDCxVpkdwTnfznrHHTvImZjDEcuPIG9qXtz3KyKjV19ldlY/770OWOGc29nPDs4fanDSi4wMeP/7/bR+Pfzyl3D77fDcc34qKYELLoBPfhJOOAHMEh2xiCQHldkRzIyDbj2IfVv2Ufd/dbx8xssc/vTh5IzPSXRoIjIK9XnzonPur/0V0MF698YuJOlh7lz4yU9g61b41a/gqKOgocGP5HjSSXDAAfDtb8MbbyQ6UhFJMJXZvcvIyWD+n+ZTeGghLa+18NLil2itbk10WCIyCg22VxBJlKIiuPxyeOEF34PIVVfBlCmwaZPvB3vOHDjySH8DpJJsEZEessuzOeyJwyiYX0DzmmZeOuUl9r6lAbpEJLaUWKeiefPge9+DzZvhiSf84DJFRfCvf/mEe84cOPxw+K//gldfhUH2/CIiMhrlVOaw8MmFPrle28yq41bRtLYp0WGJyCiixDqVZWbC4sW+3+vqaj+a48c+5ttgv/QSXHut7xd75kz4/OfhkUc0hLqIpLWc8TkcvuxwSk8oZd/b+1h1/Cp2P7Y70WGJyCgx7MTazA43sz+a2VfMbMbwQ5IhycuDc8/1A87s3Ol7Fbn4Yhg3zg+j/vOfw1lnwdix/qbIn//c3xip2myRtKIy2zcLOfSxQ6n4YAXtde2sPms1VTdUMdhxHUREIsWixvrTwOnOuR8Bx4WGspUEys2Fc86B3/0Otm2D5cvhG9/wzUOam+GBB3wN9rx5MHWqT8D/+Ee/roiMdiqzgcz8TObfO5/p35oOnbDpqk2sOW8NbTUa9VZEhm7QIy/utwGzbGCcc26rmWUARc65hphENwijqh/reNq61TcJefJJ3z67urrn8oMPhtNPh1NOgRNPhIqKhIQpkk7iOfJiL/tSmR1h14O7WPexdXQ0dJAzMYeDbjuIsWeOTXRYIpKkhjxATJSNnQWMcc79t5kdDexxzq2NQZzDkkyFdMro7PQ3N4aS7KefhqaIG3kOPth363fiif5xypTExCoyisV5SHOV2QPQsqmFdR9fR8M//P8Ykz4ziVk3zCKrpL/hHkQk3cQssTaz/wKuAB5xzp0fzLsBWO2cuzMWwQ5VshXSKam11Xfn98QT8Mwz8M9/QktLz3VmzuxOtI8/Hg480A9oIyJDFq/EWmX24LgOx5YfbmHzNzfj2hw5E3OYffNsKi+oxDQQl4gEhjPyYqTjgCnAJ8Pm/Q34A5DQQlpiICfHj+R4wgn+dWsrrFzpk+xnnoFnn/X9Zm/aBL//vV+nrAyOOQaOPdY/HnMMjBmTuGMQkXAqswfBMo3pX5/O2LPG8trlr9G4vJG1H15L+W3lzPnpHAoOLEh0iCKS5AabWD/unKszs/Bq7rQZGjft5OTAu97lp699DTo6YPVqWLbMJ9rPPw/vvAOPPuqnkAMP9Il2KNlesACysxN3HCLpS2X2EBQdVsQRzx3Btlu3sfHrG6l9tJYX57/IxCUTmfHNGRoOXUSiGmxTkI8As4DJwDLgAuCDwBeccz+PS4QDlIw/K6aFt9/2TUb++U/f+8iKFfv3lZ2TA4ceCkcc0T0tWOC7CBSReDYFUZk9TK07W9l0zSa23bYNOiGzKJOpV05lypenkFWs9tci6SjWNy8eC1wMTAO2AXc5554YbpDDlSqF9KjX1uZrtcOT7Q0b9l8vMxPmz++ZbB92mB9BUiTNxPnmRZXZMdC0tomNX99IzYM1AGSVZzHli1OY/IXJZJfrFzmRdBLTxLqXjWcCpzjnnhzWhoYp1QrptNLQ4EeC/Ne/uqd163yvJOHMYPZsX5u9YAEccoh/nD3bJ+Iio9QId7enMnsY6p6pY+PVG7t6D8kszmTy5yYz5UtTyBmnJiIi6WDIibWZPQD0V4VYAexyzp029BCHL1UL6bTV3OxrtsOT7Vdf9TXekfLyfLd/oYQ7NE2Y4JNxkRQXq8RaZfbIcM5R/0w9VddXUftELQCWa4y/aDyTvzCZ4sOLExyhiMTTcHoFWR2s8xrQCVwGrALCS8OZgLqBkMEpKOi+wTGktdUPs/7KKz2nt97qTr7DjR3rm5PMm9c9zZ3rR5NUwi3pSWX2CDAzyk4uo+zkMhqWN1D1vSpqHqhh+++2s/132yk9oZTJ/zGZig9UkJGt7khF0kl/ifUP8bXadQBmdoBz7puRK5nZ/fEITtJM6CbHQw/tOb+uztdmRybcNTXdXQGGKyz0CXYo0Q4l3Qcc4PchMnqpzB5hJceUsOD+BTRvaGbrz7ey/XfbqX+2nvpn68mZlMOET0xgwiUT1FWfSJroM7F2ztVHzBpvZpnOuY7QDDM7BTg6DrGJeGVlPfvXBnDO90iydq2v5V63rnuqrvb9b69c2XM7WVk+uZ43D+bM6TlNmqRabkl5KrMTp2BOAXNunsPM785kxx072PqzrTSvb2bLDVvYcsMWSo4vYeIlE6m8sFK9iYiMYoP9dt8LvGxmzwKNwDzgPcB+NSIicWXmm3xMnQrveU/PZTU13cl2eNK9eTO89pqfIhUU+JskIxPuOXNg/Hgl3ZKqVGaPsKziLCZ/bjKTPjuJ+mfr2f677ey8ZycN/2ig4R8NbPiPDVR8oIJxHx7HmPeMISNXTUVERpOhdLc3CfgscCjQAPyvcy7hPyum8o0wMkJaWuD1132yvWFDz2nXrujvKy7uPemePRsqK5V0y7DFubs9ldkJ1r6nnep7q9l+23bql3X/qJBZmtmVZJcvLicjR0m2SCqId3d7GcDhzrmV/a48sO0dC5wAvAE8CzQBP8DXtpQDX3PONUS+L50KaYmDurr9k+3QVFsb/X1FRTBrlm9iMmtWz+fTp6tNtwzICHe3pzI7gVo2trDz7p1U31PNnpf2dM3PKs+i4gMVjD13LGPOGENmoboYFUlWMUuszey2XmZXAh3OuQ8MMb7w7V8CHOCcuzZs3u+AR5xzd5vZmcAnnXP/FvnedC2kZQTU1PSecL/xBtRHNmkNk5Hhm6qEJ9vhj+XlI3cMktTiOPKiyuwk1vxaMzvv2cnOu3fSvKa5a77lGuWLy6k4t4Kx54wld3JuAqMUkUixTKyXAU8C4W86HljmnPvuMIM8GfgGcIYLggp+wtwElDvnms0sC6gH5jvnNoe/X4W0JMTu3bBxI7z5Zs/HjRt9N4GRg+CEKyvrPemeNcsn5Fm6wSldxDGxVpmdIprWNLHr/l3semAXjS809vjEio4sYux7xzLmPWMoPrqYjCw1GRFJpOH0Yx3pUudcjzu/zGwufrjc4fox8Bxwi5nNBr4LTMcPZNAM4JxrN7NNwMnA5hjsU2R4xozx06Jevl+trVBV1TPZDk/A6+p6758bfFI9fbpPsmfO3H+qqFDbbhkIldkponB+IYXzC5l+9XT2bd/H7od3s+vBXdQ+VsuelXvYs3IPVddVkVmSSfnicsrfXc6Yd48hf1Z+okMXkTCxaGM9E3jROVcxjG0cBKwDFjrnVpvZV/F3rV8PfNQ5d2jYus8DDzjnvh++jenTp7vKysqu10uWLGHJkiVDDUkkvpzz3QJGS7rfeafv9xcVwYwZ3Yl2ZAJe1N/ge5JoS5cuZenSpQCsXLmyyjk3YyT2qzI7tXS0dFD39zp2P7Kb3Y/tpuX1lh7L8w7IY8wZYyhbXEbZSWUaVl0kTgZaZg+2Kcgmev6kmAmMB/7gnLtsqMGa2XnA0lBBb2b5QDVwHfAh59wxYev+C7jTOXdT+Db0s6KMKi0tvnvAN9+ETZv2nxr2uxesp4qK3mu6Z87UTZVJKI5NQVRmjzItm1uofbzWT0/U0l7b3mN5wcEFflTIU/zIkDnj9V0XibVYNgV5HLiT7oK6E9jhnHt9GPGF4uiKxTnXYmYb8H8EIu/wKgb6qc4TSXH5+d0jRkZyzvdU0lvCHZp27fLTiy/u/34zmDy596R71iw/WE6G2nCOEiqzR5n8GfnkX5bPpMsm4TocjSsb2f3Ybuqfrqf+H/U0r22meW0z7/zSn/KCuQWUnlxK2UlllBxXQt70PEzNyETiZrCJ9R+dcz3Gjzazg8xstnPujWHEsRooNbPxzrkdwbx2fGE8yczyg4I7G5gGPD2MfYmkNrPutt1HHrn/8s5O2LYtetL99tvd07Jl+78/O9sn3qEBeHqbxo5VG+/UoDJ7FLNMo+ToEkqOLoFrobO1k8YVjdT9Xx11T9dR/2w9zeubaV7fzLZfbwMgZ2IOJceVUPquUkqOK6H4iGINUiMSQ4NNrM8EnomYtwN4Cjh8qEE459ab2d+A84Gfm1kZkAP8N3AWsBh4CDgduM85t3Wo+xIZ9TIyfGI8eXLPYeBDWlt9jyXREu+dO30zlM2bo+8jL88n2FOmRE++S0uVfCeeyuw0kpGTQelxpZQeV8r0q6fT2RYk2kGS3fB8A63bWtl13y523ecHxbIco/jI4q5ku/ioYnKn5qpWW2SIBtTG2sw+j78ppTjKKo87584cViBmFcBPgBX4Go5bnXNrzawcP9jAZmACcJVzbk/k+9VeTyRGmpt9bfZbb+0/heb31X93SFFR9KQ7lJDrJksg9m2sVWZLb5xztLzeQv1z9TQ810D98/U9+s8Oya7MpvioYooXBdNRxeROUF/aIiEx6cfazA4FPgn8NGJRq3Nu2/BCHD4V0iIjqLGx98Q7fGre/w/2fsrK+m5yMmWKrx0f5eJx86LKbBmItto2GpY30PBcAw3/bKBxZSPtu9v3Wy9nck53or2omKKFRUq2JW3FcoCYbOdcWy/zpznntgwjxmFTIS2SRJzz/XT3lXi//Tbs29f/tior96/pDp8mT/btwlNYHHsFUZktg+KcY++mvTSuaPTTi400rmyko7Fjv3Wzx2dTdFgRRQuD6bAi8g/M1wA2MuoNuVcQM7sS+Ktzbm0w6z/NLDNyNeBU4JThBioio4SZH7K9vBwOPbT3dUJ9eYc3MYmctm7161RX9z6QTmhfEyb0XuM9bZrv7aSyMi3ae6vMluEyM/Jn5ZM/K59xF44DwHU6Wja00PBiA40rGv2ANS/voW1HG7WP1VL7WG3X+zPyMig8pJCihUUUHuYfiw4tIqtEI8lKeujvSp8JlIW9PgIoBLbT3X2T4dvRiYgMnBmMG+en3no3AejogB07em/nHZq2beueXnih9+0UFPQcUGfmzJ6vy8p6f1/qUZktMWcZRsFBBRQcVMCEj/lLx3U69m7ey56X97DnpT1dj/uq9nXVdofLm5VH4SF+dMnQY8HcAvVIIqPOYJuCHAG86pxrjZh/cFgNSULoZ0WRNNXW5pPq3mq8q6p87ya1tX1vo6ysO8k+8ECYO7d7Ki2N+yHEsSmIymwZUW21bTStbuqRcDe92oRr7SXXyIT82fndCXeQdOfPyScjWwm3JK9YDhCzwDnX4/fYYGjbA4CEFtIikqays32Tj2nToq9TV+cT7FCXguHPN23yy1et8lOkCRO6k+z5833t+sKFfhCf5KcyW0ZUdnm2H/nx5O4fTjrbOml+rZnmNc00rWny06tNtLzRQstrfgp1/wdg2UEN+fwCn2wf7Gu382fnq4Zbkt6AEmszOxH/8+EZwRC54cbg7zp/MMaxiYhzPipDAAAgAElEQVTERlmZT4YXLtx/Waitd2gI+ddeg/Xr/fTaa7B9u5/+7/+635OZ6ZPs446D974XFi9OqkRbZbYkk4zsDIoOKaLokJ7da3bs7fAD2KxppunV7oR776a9/vWrTVRT3f2GTMiflU/B3AIK5hX4x+B5dllq38Aso8dAa6zrgDvw7feOj1jWCvwmlkGJiIyY8LbeRx/dc1lnJ2zZ0p1or14NK1bAmjX++erV8Ktf+aT6wgvhu9/1N00mnspsSXqZeZkULyymeGHP7tY7mjpoWucT6+Y1fuTIpnU+4W7Z0ELLhhZqHqzp8Z7s8dkUzivskWwXzC3QYDcy4gbTj3UpcIFzLikLZLXXE5ER09Tkm408+SQ89JBPtgEmTfK13IMc+CZO/VirzJZRpWNvBy0bWrqGaW9e19z1vLOls9f3ZBRm+GYlYcl24bxCNSuRYYlJG2vnXD1RajnM7ATn3LNDjE9EJLUUFsLxx/t+tQ84AG6+GVauhHfe8e21k2BESZXZMtpk5mVStKCIogU9v1+u07HvrX1dNdvhSXfbzjb2/GsPe/4VMfhnRtCsZF5Ys5J5PunOKlXXgDJ0/fVj/QDQ318IA+YD42IVlIhIUgj1OLJ1q+/mr6rKNwNZswbWrYM9YX+sMzLgmmt8sp0gKrMlHVmGkTc9j7zpeYx5z5gey9p2t+1Xw93VrOSNFlre2L9ZSc7EnO4a7rBa7pxJOWpWIv3q79+y1/Bt9bbS3QdqJAPOiWVQIiJx0d4Ou3dDTU30aefO7kR6505/c2M048bBokX+5sULL0xoUh1QmS0SJntMNqXHlVJ6XM9uM7ualazrTrqb1jXR8loLrdtaad3WSt1TdT3ek1mc2WvCnXdAnkablC79JdY3Ao3OuT7HHTYz3V0uIvHX2elvIqyvh4aGnlNj4/7zwpfV1vr3DYYZTJzoh02fPNnfmHjwwb5HkIMPhoqK+Bzn0KnMFhmAvpqV7K3a2yPhDtVyt9e0+yHeX+w5+I1lG/mzu5uVFB5SSNGCYHh39ceddvpMrJ1zu/paHuaD6C5zEYmntjY48URYvnzo28jI8MOsjx0bfaqo8DXPkyf7PqyzUqe9pcpskeGxDCN/Zj75M/MZe/bYHstaq1u7k+2whHtf1b6ueT22lW2+VntBMNpk8Jg3PU9NSkaxYf3FMLNpwOeAz6NCWkTiobERnnnG3yAYSqqPOQZKSvafioujzy8r81NG+tYgqcwWGbqcyhxyKnMoO7Gsx/yOpg6aXw+ak6z13QQ2vdLE3o17aXrFPw+XWZzpE+2wZLtwQSE5FTkjeTgSJ4NOrM0sEzgXuAw4A9iBb9MnIjJ0zsGOHd19Rq9dC88957u16wy60ioogHvvhbPOSmysKURltkh8ZRZmUnx4McWH9+yPu31PO81rm31y/WoTe17xw7u37Wij4fkGGp5v6LF+zoSc/ZLtwoMLySzMHMnDkWEacGJtZjOBS4FLgLFAFfA+4FHgXXGJTkRGB+d8zXPoxsC33vLT2293P27c2Hsb6Kwsf4PgOefAv/+77yta+qUyWySxsoqyKDm6hJKjS3rMb61u7Uq2m17pTrhbt7fSur2V2idqu1c2yJuVR9GCou5ke0Eh+XPydcNkkuo3sTazC/A1HacBbwA/Bn4PXO2ceyRY7bm4RSgiyaWz0yfAtbU9p7o63+NGdbVPoHfu7Pm8tbX/bZeVwdy5fjroIDjqKDj2WN9vtAyIymyR5JZTmUPOaTmUn1beNS9002Qo2Q49Nq9vZu+be9n75l52/aX7FgrLsR43SoZquTXSZOINpMa6BT8E7nPApc651wHMbGBDNopIcmhrg+ZmX3Mc6kFjoI91dX4K9awxwBFbeygshMpKX+M8daq/QXDq1O7nM2b45fqjMFwqs0VSTPhNkxXv6+5tqLO1k+bXm3sk202vNLF3816aXm6i6eUmdrKza/2ssVkUH1FM8ZHFFB1ZRPGRxeTN0M2SI6nfxNo59xDwUHDTy2VmNgO4H8gOrWNm5gY6NrqI7K+9Hfbu9YlvaGpq6vk43HltbbGLt6TE965RXu5rmUPPy8t9386Vlf4x9LyyUrXOI0RltsjokZGTQdEhRRQd0rNbwPbGdprWNPWo4d7z8h7aa9qpfbyW2se7m5NklWdRdIRPsouPKqb0hFJyJ+SO9KGkjcEMab4F+IaZZQEfAKaY2XeB3wLvBpbGJ0SRJNHRAU88AZs3+yR4377YPYZuzounzEx/81+o54zi4p7Poz2GetQIJdGlpSnVBV26UpktMnplFWdRemwppcd2D3zjnB/avXFlI40rG9nzrz00rmykbWcbdU/WUfdk9z3L+bPzKT2plNITSyk7sYy8WarVjhUbTqWFmc0GvgBc5pwriFlUQ7Bo0SK3YsWKRIYgo9WGDXDnnfCb3/ib7OIhIwNyc33iW1jY++Nwl2Vnq5lFkjKzlc65RSOwH5XZImnEOce+rfu6kuyG5xuof66ezqaelTm5U3MpX1xO+enllJ1WRu5E1Wj3pa8ye1jVTs65N4ArzGzlcLYjklTa2uCFF+DRR+HPf4ZXX+1edsABcOqpkJ/vE+G8vL4fB7JOXp5qgGVEqMwWSS9mRt6UPPKm5FFxrm+73dneyZ5Ve6hfVk/9snrqltWx7619bL99O9tv3w5AwfyC7kT75DKySvQ3aqBicqacc3fEYjtmlgO8CFzhnPs/M8sHfgA0AuXA15xzDX1tQ2TQ2trg5Zfhqafg73+HZct82+SQ0lJ43/vg4ot9Up3GA4zI6BCrMhtUboukmoysDEqOKqHkqBKmfnkqrtOxZ/Ueap+o9U1GnqmjeU0zzWua2frTrZAJJUeV+CR7cRml7yolI1d/B6NJtn9BvgrMCHv9C+AR59zdZnYm8Gvg3xIRmIwSzkFVlR/Bb/lyXzO9cqVv6xxu3jw47TQ491w45RTI0YhYIlGo3BZJYZZhFC8spnhhMdOunEZnaycN/2yg9slaap+opWF5Aw3/9FPV9VVk5GdQemIp5aeXU764nKKFRViGmhmGJE1ibWbHAduA2uD1JOAi/PC7AE8A95nZDOfc5oQEKall3z5Yt87XRq9e7R9ffhl27dp/3dmz4eSTfTJ96qkwceLIxyuSYlRui4w+GTkZlJ1URtlJZcz8zkzaG9qpe8bf/Fj7RC1NrzZR+1gttY/5nkeyxmZRfqqvzS4/vZz8A/LT+kbIpEiszawQON8592Uz+2Yw+xRgl3OuGcA5125mm4CTgc0JCVSSU3u776lj3To/vfKKT6DXrfPLIo0dC0cfDccc4x+PPtrPE5EBU7ktkh6ySrKoOKeCinN8G+192/dR9/e6rhrtfVv2UX1vNdX3VgOQOy23qza7fHE5OePT6xffpEisga8BN0TMmwzURMxrBHodz7i6uppFi7pv0FyyZAlLliyJZYySaE1N8PrrPmFev7778fXXex/VzwzmzIHDDoNDD/WPhx0G06apdwxJuKVLl7J0aVePdxV9rZukhlVuq8wWSU25E3IZf9F4xl80HuccLW+2dNVm1/7dJ9rbb9vO9tv8jZCFhxRSfkY5Y84aQ9lJZSnbPnugZfawutuLBTN7D5DrnHsgeL0ZuBhYBFzgnDsmbN1/AXc6526K3I66bhol6urgjTfgzTf9Y/i0fXv0902Z4ttFz50L8+f7BPqQQ6CoKPp7RJLESHW3FyuxKLdVZouMPq7TseelPV212fXL6uls6e7aL6Mwg/LTyhlz9hjGnjWWvOl5CYx26OLW3V6MXAkcHtYepxw/StiPg+fhioF3Ri40ibn2dt8X9ObNftq0qWcSXRNZ2RUmJ8d3dzdvXncSPW8eHHigH8REREaKym0R2Y9lmB9S/Yhipn11Gp37Oql/vp7aR2up+WsNTaubqHmwhpoHa9jABgoOLmDs2WMZc9YYSk8oJSMnNWuzwyVDjXUlEN4T+fPAl4HlwFqg0jnXYmbZwB5glnNua+R2VPuRJNrbYevW7sQ5PIHevNkn1R0d0d9fUOBvJJw92yfRoeezZ8PkyX70QJFRJgVrrIddbqvMFkk/e9/ey+5HdrP7b7upfbyWjsbufCCzKJPyM8qp+GAFY88ZS3Z5dgIj7VtS11g756rDX5tZB1DtnNtiZg8Bi4GHgNOB+3pLqmUEhRLnqqreE+e33uo7cTbzCfKMGd1TeAI9YYLaP4skOZXbIjIUeVPymHTpJCZdOonO1k7qn6tn9199ot30ahO7/ryLXX/ehWUZZaeWUXFeBRXvr0ipkSATnlj34zPAD8zsMGACoDtb4q211SfHocQ58rG/GmeASZN6Js4zZ3Y/nzrVjzYoIqOVym0R6VdGTgblp5RTfko5B/zgAPZu2UvNgzVU/6mauqfrqH28ltrHa9nw2Q2UvKuEygsqGfeRceROSO4cIuFNQWJFPysO0N69sGVL70nz5s3wzjt+EJW+TJoE06f7KTxpnjnTJ855qXkzgkiipFpTkFhQmS0i0bTVtLHrQV97vfvR3bh9QV6SAWPePYbxHxtPxQcqyCxMTPPQpG4KIjHW1BS9trmqqu+eNcAP1z1lik+Up0/v+agaZxEREYmz7LHZTLx4IhMvnkj7nnZ2/3U3O+7c4ZuNPOKnjMIMKs+rZNJnJlFybEnSDEqjxDrVNDREr22uqup9VMFwWVk+OQ5PlsMT6ClTIDt5bxgQERGR9JFVlMW4C8cx7sJxtNW0sfN/d7Ljjzto+EcDO/6wgx1/2EHR4UVM+uwkxl80nsyCxHZyoMQ62TQ3d/fj3FsCXVfX9/tzcvwAKL0lzTNm+GYc6llDREREUkz22Gwmf3oykz89mZaNLWy7dRvbfrONPav28Pplr7PxPzcy5YopTL5iMtlliakkVGKdCK2tvheNDRv8qIHhj2+91fd78/P3T5bDHydM8M05REREREap/Fn5zPr+LKZ/azrV/1vN1p9tpfGFRjZ/ezNv/fgtplwxhalXTiWrZGRTXSXW8dTS4ofdfuUVePVVWLvWJ9CbNkXvWSMrC2bN8l3PzZy5f/JcWanu6ERERESAzLxMJnx8AhM+PoG6p+vY/J3N1D1VR9V3q9h26zZm/WAW4z82fsTaYCuxjgXnfFONlSt9Av3KK3564w3o7Nx/fTOfKB94IMyZ4x9Dz2fM8Mm1iIiIiAxY2cllLPz7QuqW1fHmV9+kcXkj6z+xnh137GDuHXNHpD9sZXBDUV8PL74Iy5d3Tzt37r9eZqYfdnvBAj/Nnw8HHeQHRFGXdCIiIiIxV3ZiGUc8dwQ7/rCDN698k9onallx6AoOvudgyk8tj+u+lVgPxL598Nxz8Pjjflq5cv++nseOhUWL4LDD4JBDfCI9d64SaBEREZERZhnGhE9OoPzd5az/xHpqn6hl9ZmrmffHeYy7YFzc9qvEOpqWFvjrX+Huu+Hhh31vHSHZ2XD44XDMMd3TAQeo7bOIiIhIEsmdmMuhjx7KG196g60/3cq6i9aRPTab8tPiU3OtxDrSpk3wk5/AbbdBY2P3/AUL4Iwz/HTSSVBQkLgYRURERGRALMOYffNsLNt4+0dvs+b8NRy19qi4DI+uxDqkpgauuQZuvbX7hsNFi+AjH4ELLvB9Q4uIiIhIyjEzDrjxAJpebaL20VrevPJNDv7jwTHfjzo8Bli2DObNg1//2vcB/fGPw6pV/gbFr3xFSbWIiIhIirNM48BfHYjlGDvv3EnLppaY70OJ9d//DosXQ3U1nHwyrF4Nd9wBCxcmOjIRERERiaH8GfmM+/A4cLD999tjvv30Tqxranwzj7Y2+PSn4YknfM21iIiIiIxKlRdWAlD7aG3Mt53eifVNN8Hu3XDqqfCzn2lgFhEREZFRruzEMgAaVzXiOlw/aw9O+ibWzsFdd/nn3/mOH8xFREREREa1rNIscibm4PY59r61N6bbTt/Eevt2Pwx5aSkcf3yioxERERGREZI71Xe117qtNabbTd/E+vXX/eO8eb4nEBERERFJC1mlvvlve317TLebvhnl7t3+cVz8hrUUERERkeQTSqw76jtiut30Tazr6vxjaWli4xARERGREWW5BkBna2dMt5u+ifW+ff5RQ5OLiIiIpBXL9In1qO0VxMzOM7NNZlZjZj8xs6xgfoWZ/drMbjSzH5tZdkx22BFU/at9tYjIoI14mS0iEkOhxJrYtgRJjsTazKYBHwDOB74EXAp8MVj8v8BtzrmvAa8D18Vkp51B1b+62RMRGZSElNkiIjFkWUGNdXtsa6yTZUSU6cCnnHPtwEozOwQ41cyeBQ5yzi0P1nsQeM3MrnfONQ1rj6qxFhEZqpEvs0VEYmhUNwVxzi0LCuiQrcAW4DRgU9h6W4Oni4a9U9VYi4gMSULKbBGRWArSv1GZWPfiKOBXwGSgJmJZIzBp2HsIJdZmw96UiEiai3+ZLSISQ5YR5H+xzauTpilIFzObDVQ75142Mwe0RKySC7RFvq+6uppFi7orRZYsWcKSJUui78gFZ1JNQUQkAZYuXcrSpUtDLysSGctwjFiZLSISDwNMrAdaZidVYm1mmcDlwNeDWVuBA8OWG1AIvBP53srKSlasWDESYYqIDFt4ImlmuxIczpCozBaRlDXIBgsDLbOTrbr2SuAm51zQyTSPAweFLZ8BtAIqjUVEEk9ltoikNOdGaRtrM7saWAXkm9ksM/sUUAdsC+44BzgTuMU515qoOEVERGW2iKS4ON1ilxRNQczsWuC7EbPXO+duM7MPA980s81AEXDNSMcnIiLdVGaLyKgxGm9edM5dD1wfZVkV8O9x2GnMNykikg4SUmaLiMSQxalXuKRpCpIw6m5PREREJD3FuJ5VibWIiIiIpJdQvaoSaxERERGRYYhTgwUl1iIiIiKSlkZtd3siIiIiIiNCNdYxpl5BRERERNKb2ljHmHoFEREREUkr6m5PRERERCSWVGMtIiIiIjIMamMtIiIiIhJDqrGOkaIimDgRiosTHYmIiIiIjKCs0ixyJuaQURjbVDgrpltLJZ/7nJ9EREREJK1Mv2Y606+ZHvPtpm+NtYiIiIhIDCmxFhERERGJASXWIiIiIiIxoMRaRERERCQG0j6xXrp0aaJDGLbRcAwwOo5jNBwDjI7jGA3HIPtLtc811eIFxTxSUi3mVIsXEhOzEusUvFAijYZjgNFxHKPhGGB0HMdoOAbZX6p9rqkWLyjmkZJqMadavKDEWkREREQkZZlzMR5yJkHMrBqoGsJbK4BdMQ5npI2GY4DRcRyj4RhgdBxHKh3DdOdcZaKDGElpVGanWrygmEdKqsWcavFC/GKOWmaPmsRaRERERCSR1BRERERERCQGlFiLiIiIiMSAEmsRERERkRhI28TazPLN7BYz+56Z/dLMShIdU2/M7Dwz22RmNWb2EzPLClv2RzNzwXR/2PwKM/u1md1oZj82s+zERN+TmV0fFu/Lwbyon0OyfUZmtjws/tB0drBsv2ML5ifFMZjZGUH8MwYSW7J+LlGOI+p3JFieUt8TGfg1ZmZnmdnPzOwXZnb+SMcZFkef12DYenPMrCPsejx8pGONiKfXcitinaQ4x0EsUcvgiPXKzKwhbJ0PjnCcgypve3n/iJ/zoZStEe8f0Wu7t3iD+f1e08F68TvHzrm0nIDfAR8Onp8J3JXomHqJcRpwB3Ak8AmgCbgybNmPgdODaWrY+54Cjgmefxr4fhIcSyHwq7B45/T3OSTTZwTMA24GTgSODaZXgbxox5YsxwCMBz4IOGDGQGJLxs+lt+Po6zsStjxlvieaBn6NAQcC/wIygtfPAIckINY+r8GIdW8E3h1ciycn+BxHLbeS7RwH+45aBvey7leDsiJ0bJkjGOegy9tEn/OhlK29bGPEru0+znG/1/RInOMRudCSbQImAfuAguB1VnDRzEhkXL3EeSKQFfb6B8DDwfP/B5wfvjyYfyzwTtjrycAeoDDBx/LFYMoLmxf1c0i2zwiYFvF6OnB3tGNLtusM/+tUeKE5pHOf6GPq5TiifkeC1yn1PdE08O8N/g/o98JeXwP8LgHx9nkNhs0fDzwETEn0OQ7i6bXcSsZzHOw7ahkcMT8XeJwoSdUIxTrg8jZZzvlgy9aI9474tR0ZbzCv32t6JM5xujYFOQXY5ZxrBnDOtQObgJMTGVQk59yyILaQrcAWMysA3gfcE7xeHLbOafhjCW1ja/B0UbzjjcbMDPgI8CNgh5n9W7DoFKJ/Dn0tG3HOuS0Rsz4I/KWPY4MkOgbnXGfErFMY2rnva1ncRR5HtO8IQKp9T6TLKQzsGuvxGQKvBe8dUX1dgxH+DTgVeMvMfmdmhSMSYC/6KbfCJcU5huhlcC+rngMcA7xuZg+bWUXcg4swyPI2UkLO+WDK1l6M+LUdGe8grmmI8zlO18R6MlATMa8R/19lMjsK+JVzrtk5Nxt/HA8DfzOzBcE6SXdszjsW31H7LcB/m9lZ9B1r0h1HhLPx/71HOzZI7mMY6rlP5mOC4DsCkGrfE+ky0M8mcr1k+fy6rsFwzrmbgRLgXPzP1L8d4bjCY+mr3AqXrOcYgjI4cqZz7j6gFJ+0zgAeMLNE5zqDKW+S9Zz3el1Dclzbg7imIc7nOGpD9FHOAS0R83KBtgTEMiBmNhuods51NcZ3zm0DLgtq5j4PXE4SH5tzrha4NvjP8ov4n+v6ijUpjyOoAdnrnGsIzevl2P5GEn8W9B9bKn4u+31HIPW+JzLgzyZyvYR/ftGuwRDnXAfwoJm9CbxkZhOD6zMhopRbPVYhyc4x9F4Gh3P+9/1nzOw0YB2+Bvv5EQxxv5AYeHmTdOe8v+sakufaHsA1DXE+x4n+Ly5RtgLlEfOKgXcSEEu/zCwTnwx8PcoqP8XfaAARxxZcXIUk17GF4u3rc0jmz+hc4IEoy6J+FoFkOYahnvukPKYBfEcg9b4n6Wqg11jkegm9Dgd4DQLgnFsLPEn39Zho4d+NcEl1jsP0VQZ3cc7tAO4m8ed5MOVmUp3zwVzXkFTXdrRrGuJ8jtM1sX4KmGRm+QDmu9maBjyd0KiiuxK4yTm3L8ryTvwdruBrgQ8KWzYDaAVWxC26wQvF29fnkMyfUV+FevhnkczHMNRzn6zH1N93BFLve5KuBnqNRX6GBwbvTZSBXIPhmoD1cYxnMMK/G+GS7RyHDCixDnQAL8UxloEYTLmZbOd8sNc1JMe1He2ahjif47RMrJ1z2/F3sIZuZjoduC/sBqakYWZXA6uAfDObZWafMrNDzOzcYHkO8En8Hbs451YA28zskGATZwK3OOdaExA+AGY2ycwuMrOMoGbwK8A1fX0OyfoZBTdl5Dnndgavez02SK7rLIgNwPqLLZk/l8jjCOb19h2ZbWZzU+l7Il5f15iZXWlmoT+IvwLeA13XxWn4WqoRF+UanGNm15nZhGCdi81sXPD8ROBp51x9guKNWm4l6zkOiSyDg3ldMZvZ+WY2K3g+F9jtnHtthGMccHkbGT8JOueDLFst0dd2ZLx9XdPB8pE7x0PtTiTVJ/zPALcGJ/4WoCjRMfUS47X4tkDh0zr8l3MnsBz4CTAp4n3T8TcPfAPft2RWIuIPi2c+/m7iV4MLeu5APodk/IzwXbd9cSDHlizHgP+Z6zPB9fMdoHI45z5Rx9TbcUT7jgTrp9T3RFOPz6bXawxYCZwftt7H8P2U/xh4f4JijVZO5+N7Hjg2WO9+YBdwF/CJBJ/fvsrkpDvHEbH3KIMjYwZ+BtQBfwL+g6Cv4hGMbyjlbULP+RDK1oRe21Hi7e9v8YidYwt2ICIiIiIiw5CWTUFERERERGJNibWIiIiISAwosRYRERERiQEl1iIiIiIiMaDEWkREREQkBpRYi4iIiIjEgBJrEREREZEYUGItIknFzCrN7AEzqzKzryc6HhER6ZvK7W4aIEZEkoqZvQd4CigB1gMHOud2JzYqERGJRuV2NyXWIpK0zOwZ4DTnXHuiYxERkf6le7mtpiAikpTMbCrwWLoWziIiqUblthJriQMz+08z+3mi4wiJRTxmdrqZPW9mp8QorFHFzG6K8fYMuBC4IZbbFZGeVF4Pn5nNMLNbzezXMdpeppldZmZVA1j3vWa2fKjHGsuyW+W2p8Q6xZjZPDP7k5n9xsxeNjNnZh9LdFwR/grcleggwsQinreAw2IQS1Iys7PNbJeZbTOzd5nZ1OA6c2b2MzMrC9Y72sxeNbMfmVl+MG88sD3GIV0E/NY51x5sXyTlqLweklQsr3cGj7kx2l4GUAdMG8C6K4B50Raa2Xwzuzu49q41swlhy2JddqvcRol1KvoL8B3n3KXA4cAdoQVmdqiZnZyIoML37Zx71Tn3bCLi6M1g4ol2Dp1zrwG7YhlXLD+v4W7LOfdX4BdAu3PueefcW8BlQCuw3DlXF6z3Ar4g/6pzriV4+/uAh4Z1AGHM7BvAt4EnzGw9cFCsti0ywlReD1IqltfOuWbgnVjtyznXBqwc4Lo78El4tOVrgKXBy5ucc+GJdMzKbpXb3ZRYpxAzGwMcCDQBOOc6gW8CzsxK8YW2JSCuhO07lgZwHJ0juK9EbOtuYIqZHQngnKsBngTOCdtXGfBOcO2FzHXOrR/mvrs4577rnJvjnFvknJvrnHsmVtsWGSkqr+MrCcvrWPcEMZj4+1u3I3iMbPccs7Jb5XY3JdYpJOi65jXgETM7NphXBbyAb9c0E1hiZpcAmNl1ZvaN4GegrnZUZrbIzH5sZl80s2Yze83MvmxmHzCzF83sg2b2ZzNrMLPLw97X6/Yi9n2Zmd0e3tYs+Dn0FjP7vpk9ZWZHB/P73F/Y+83MLjezNWZ2hpmtNrMdZnbxAPYxJzyefva53znsRXnw026jmV0WcU6/HZybPwUF8UDP2X77MrPzzewFM7vYzDaY2Q19bK+3z77XePoS1GysAS4Im10JvNeCZh/AucADYXHmA81hr3IzxOMAAAbWSURBVOea2Q1mdpOZPW2+ScmHzWyFmb3fzO4yszozu8LMjgyOscbMzu4vPpFUovI6fcrr3kQ7xr72Y2Yzzew287W/PxzIfnrZ71/N7L/N7Ph+1lPZHS/OOU0pNAHzgSr8f6i3AePClm0GTgmeHwJsDJ5X4v+bHov/j3sLsChY9nPgueB5DrAb+Bm+L8p/A7b2tb3IfePbmP0auD2Yn4dP1sqC1+cD1UBxX/vr5bgnBPu8Olj3emAvML2ffUTG0+c+w89hLzFsBn4EFATv2xLMLwH+J2y9f+ALxQGdsyj7GoP/ee9m4Gjg9L62F/HZ9xrPAK+va4E3g+cHAN8H9gHnhV0vFrb+ucCJwfMC4EWgIHj9THAeSoB64DogH/hw8PpC/D/3XyG4BjVpGk0TKq/TorwOln97gOexr8/6RWBe2HvcAK+z0Of5LuBbvSw/JdhPVtg8ld1xmlRjnWKcr1VcgC9oPgGsMbPebtLYAHzczLKBE4N5Rfgv8lQg1D52Hf4Lg3OuFWgE7nXONeC/aBP62V5kfPvoeTPEe4EWF7TRBf6E/0Ke18/+Ircb2uadwbrX49v/ntHPPnrEM5h9RvGg8+3pXgQmhh3j2KBG6Yv4cwoDPGdRjnc3vgC73zn3gnPuiUFsL1o8A3E3MMt8c5APAbcCjwIXmlkR0OiCkjdwPPBc8Pxc/B+vUC3IB4HvBue5Dvi78+2yVwAlzrl7nP95fBWD+wxEUoLK6/Qor3sR9Rj72M/pQKlzLhTPikHu8yjgo8657wxwfZXdcZKV6ABk4MwsE5jqnNsM/IeZ/RZ/B/VS4JjwdZ1z+8xsMvA14PehTTjndprZCmAx/j/qmfhkquutYc/bCZoLRdtelFDDtzEHyA6Lq9PMNgNT+tpff5xze81sI1AGVAxiH0PeZ4R2ur8/04BNzrmbI1cayDkzs4/jE9iQg5z/ydiFxzqIzyBqPP1xzm0ws1X45iClzrmNZnY3vhbpAuDhsLhD10ao/d50fA1TaFs14ZsOe95BT52oWZqMMiqvexxfOpTX4aKexz4+m0Px/0QM1WnAHDMrcs7tibJOZ3AMKrvjSCcktWQT1v7VOfcy/qe2uZErmtlpwGedc9c738NDuM8BR5nv9mk1cGN/O+5ne33Zgq8BDe+GyPBtD4crB3g9zvsYiG30bIeMmR0ziHP2ALAwbOr17vJBbK/XeAZ+OPwPcAmwMSy+TOAK/M+mIcfg24uGbAVOMLPCsP322c5PZBRTed1TOpXXUY+xj/00BO/JYWhuwtf2/yTK8nrXfdO5yu44UmKder5gZpPCXk8GHgmet+Jv1piL79qp1Mxyzez0YPk4M6vAt/X7Lf5nnOcIfloMZBDxX7qZWT/bi9y3hW3jL8Ae4CPBtsqC/YVqPqPtL5qJwTqTg/c+MoB9hMfT3z7DjyNSb7UXhq+FKgQeMrP3mL9Jpr/PoGtfwETn3PqwqS1sF5lhz/vaXnjcvcZjZjlmtjRo5tGXe/A/Qd8L4JxrBP4G/MP17A3k3fhmIiGh832X+b6wPw+E+jKN/AwiP+eU7qFAJAqV16RNeT3Q89jrfvBlaQ7+ny/wv05gZpMGUHYbvvb448AFZnZh8N6vmNm1+Nr6J8PWV9kdT/FouK0pPhP+hggH1AJ/BG7H/5RUHiz/FvA28H78Tzsb8DUE78XXdNyOL6QewN9I0Rpsrxn/RTsLaMPfIDMVuCpY/qm+thex7yuAZfh2a/ODZYvwNZ0/wrc1PCSYH3V/UY7fBe+/At80YUHYsmj7mBEeT3/7DD+HEft+P/4mvp8F5yL0vkuC5ScDrwafzY+C8zzQc/b+Xo71/GB/dwGTg3l9faY9thUlnin4gQzuHMC19tOI1x8GTo2Yd1Mv7zsJWAvUANcG887GX2u3AJMIuhwDPo3/w3trcKznJPo7pklTrCZUXqdTeT0NeBZf676wn2Ps67N+L/6Xwn/ibyR/HfhscOy9lt34X0Xa8DXVY4C/45P6S4HL8f15/wjfDCX0HpXdcZwsOHmSJsxsHPBp59x1wesMfC3K5c65axMaXD/MzAEznW+zKENkZl91zg2pKycRGTkqryWcyu7UoKYg6efz+GYBY8LmnUbYTWnJKOynJ/3sNERmlmlmHwXuS3QsIjIgKq9FZXeKUWKdfn6Fb4u1zszeAh4H3nbOPZ/YsKILbua4Inh5sflu32TwioC/OOc29rumiCQDldcCKrtTipqCiIiIiIjEgGqsRURERERiQIm1iIiIiEgMKLEWEREREfn/7daxAAAAAMAgf+tJ7CyKBmINAAADsQYAgIFYAwDAQKwBAGAg1gAAMBBrAAAYBLDUxMzK3i7IAAAAAElFTkSuQmCC\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/fire-II-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/fire-II-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/fire-II-earth.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }