{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 33 - Stardust - 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 Stardust sample return aeroshell.\n", "https://en.wikipedia.org/wiki/Stardust_(spacecraft)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "# Set up the planet and atmosphere model.\n", "planet=Planet(\"EARTH\") \n", "planet.h_skip = 125.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('Stardust', 45.8, 60.0, 0.00, 0.52, 0.0, 0.23, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(125.0,0.0,0.0,12.6,0.0,-8.2,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+17YcXAAAgAElEQVR4nOzde5xVdb3/8deHYYCR+2XkJjAiInjHRsHqKCresswIS6tfdTTpZ8eO2S+PpmmWdUrrmB5vxTHPUdPsFJZpibfUTAEFAQUEQRnu4nAHQUDm8/vju/bMns1cmbX32nv2+/l4fB9r77X3Xuuz18z+zGe++7u+y9wdERERERFpmw5JByAiIiIi0h6osBYRERERiYEKaxERERGRGKiwFhERERGJgQprEREREZEYqLAWEREREYlBx6QDiEu/fv28oqKixc+vrq6mvLw8ewHFTPFmX6HFrHizK5fxzp49e727F87BiUFrc3ZKof0e5Tsdz3jpeMYrX49nUzm73RTWFRUVzJo1q8XPr6ysbNXzk6Z4s6/QYla82ZXLeM1seU52lEdam7NTCu33KN/peMZLxzNe+Xo8m8rZGgoiIiIiIhIDFdYiIiIiIjEo2sJ68uTJSYfQKoo3+wotZsWbXYUWb7HQzyVeOp7x0vGMVyEeT3P3pGOIRWVlpefjOBwRkeaY2Wx3r0w6jlxSzhaRQtVUzi7aHmsRERERkTi1m1lBRESS4u7U1NTg7nTsqLSaDe5QUxNul5QkG4uISGPUYy0i0kbLli2jY8eOjBw5MulQ2q2XXoKOHeGkk5KORESkcSqsRUTaqEOHkEprUl2qEjuzsGwnpwWJSDulwlpEpI1UWGdfdIhVWItIXlNhLSLSRiqssy/VY61DLCL5LOeFtZmdbmYzzawibd1EM1tmZhvM7DYz65j2WD8z+5WZ3WRmt5hZaa5jFhFpSnsurPMlZ2soiIgUgpwW1mbWH+gGnJC2bihwHjAJuAL4GvCttJf9HrjX3a8C3gJ+mLOARURaoL0W1vmUszUUREQKQU4La3dfBzyasXoYcJG7z3b3+4E7gVMAzGwccJi7z4ye+xjwTTPrmquYRUSa014L63zK2RoKIiKFIOdDQdy9JuP+i+7+Ydqq1cCK6PapwLK0566ObhbVFcpEJL+118Ia8idnayiIiBSCfDx58Xjgl9HtwcCGjMe3AYNyGpGISBPac2HdAjnJ2RoKIiKFIK8uEWZmI4Bqd58XrXJgZ8bTOgN7Ml9bXV1NZWVdp8jkyZOZPHlytkIVEam1P4X1lClTmDJlSupuv/ijyr5c5mwNBRGRJLU0Z+dNYW1mJcDXgavTVq8GRqY9x4CuwJrM15eXlzNr1qxshykiso/9KazTC0kzW5+VwLIo1zlbQ0FEJEktzdn5NBTkO8DP3X1X2rqngcPS7lcAuwFV0CKSN4p0KEhOc7aGgohIIUhiHmtL3Uxbdw0wBygzs+FmdpGZjXD3WcBaMzsyeupZwO3uvju3UYuINK49F9b5krM1FERECkFOh4KYWXfgS9Hdr5rZHYSvEm/MeOoid783uv154HozqyLMp3ptLmIVEWmp9lpY51PO1lAQESkEOS2s3X0bcHfUUn4UtcZesxy4OMuhiYjst1RhvXfv3oQjiVc+5WwNBRGRQpBPY6xFRApSe+2xzicaCiIihUCFtYhIG5WUlAAqrLNJQ0FEpBCosBYRaaO68/vAVfllhYaCiEghUGEtIhIDDQfJLg0FEZFCoMJaRCQGKqyzS0NBRKQQqLAWEYmBCuvs0lAQESkEKqxFRGKgwjq7NBRERAqBCmsRkRiosM4uDQURkUKgwlpEJAYqrLNLQ0FEpBCosBYRiYEK6+zSUBARKQQqrEVEYqDCOrs0FERECoEKaxGRGKiwzi4NBRGRQqDCWkQkBiqssytVWO/dm2wcIiJNUWEtIhIDFdbZ1blzWO7enWwcIiJNUWEtIhIDFdbZ1aVLWH7wQbJxiIg0RYW1iEgMUoX1Xo1VyIpOncJy1y6NsxaR/KXCWkQkBiUlJYB6rLOlQ4e64lrDQUQkX6mwFhGJgYaCZJ+Gg4hIvlNhLSISAxXW2afCWkTynQprEZEYqLDOPhXWIpLvVFiLiMRAhXX2pQrrnTuTjUNEpDEdc71DMzsd+BHweXevitaVATcD24DewFXuvrW5x0RE8kV7LazzKWf36hWWmzfHsTURkfjltMfazPoD3YATMh66C/iHu18DPAr8qoWPiYjkhfZYWOdbzu7TJyw3boxriyIi8cppYe3u6wiJtpaZDQK+ADwWrXoGONfMKpp6LCcBi4i0UHssrPMtZ/fuHZYqrEUkX+V8KIi715hZ+qrxwHp33xE9/qGZLQNOBvY08VhVLuMWEWlKeyysIb9ydqrHetOmtm5JRCQ7cl5YN2AwsCFj3TZgEPBhE4/VU11dTWVlZe39yZMnM3ny5HgjFRFpRGsL6ylTpjBlypTU3X7ZiSorEsvZqcJ6Q+YeRESyrKU5Ox8Kawcyz/HuTOj5oJnHapWXlzNr1qz4oxMRaYHWFtbphaSZrc9aYPFLLGf37x+Wa9e26mUiIm3W0pydD4X1asKZ4+m6A2sICbyxx0RE8kZ7HQrSgMRy9rBhYbl8eRxbExGJXz7MY/0cMCiaogkzKwWGAi8085iISN4oosI6sZytwlpE8l3OC2urOwvGANz9XeBx4LRo/QRgqruvbuqxHIYsItKs9lpY51POHjo0LJcvB/c4tigiEq+cDgUxs+7Al6K7XzWzO9y9GrgUuNnMjgEGAOlnsDT1mIhIXvF2VPHlW87u2RPKy6G6GlasqOvBFhHJFzktrN19G3B31NLXbwIuaeQ1jT4mIpIvPvzwQwBKS0sTjiQ++ZizjzkGnnkG5s1TYS0i+ScfxliLiBS8VGHdsWM+nBPefh1zTFjOm5dsHCIiDVFhLSISAxXWuTFmTFjOnJlsHCIiDVFhLSISAxXWuXHyyWH5wguwZ5/ZsUVEkqXCWkQkBnv37gVUWGfbQQfByJGwfTvommAikm9UWIuIxEA91rlz5plh+cgjycYhIpJJhbWISAxShXVJSUnCkbR/n/98WD78MLSzacNFpMCpsBYRiYF6rHPnxBPDVHurVsHf/pZ0NCIidVRYi4jEYPfu3UD7msc6X3XoAJdEM2X/4hfJxiIikk6FtYhIDLZt2wZA9+7dE46kOHz961BWBn/9KyxYkHQ0IiJBLIW1mX0iju2IiBSiXbt2sWfPHkpLS+ncuXPS4TSrPeTsfv3goovC7WuuSTYWEZGUVhXWZnaWmf3DzJaa2TtRqwKmZic8EZH8l6+91e09Z3/ve9C1K/z5z2FeaxGRpLX2LJsHgJ8Ac4DUudgGfCHOoERECkm+Fta085w9YABceSXccEMYGjJnThgeIiKSlNYW1u+4+y2ZK81sUUzxiIgUnO3btwN5WVi3+5x91VXwv/8LCxfCd78Lt96adEQiUsxaO8b6UjOb1MD6c+MIRkSkEG3atAmAnj17JhzJPtp9zu7SBe67D0pK4Lbb4Le/TToiESlmrS2sPwbcY2Z701oNcHcWYhMRKQhr164FYMCAAQlHso+iyNmVlXXT7l10EcycmWw8IlK8WjsU5IfAZGAWsDdaZ8D/iTMoEZFCkiqsBw4cmHAk+yianH3ZZfD663DPPXDWWeHCMWPGJB2ViBSb1hbWTwCPuvuu9JVm9l/xhSQiUljeffddIC8L66LJ2WZw112wYQP88Y8wYQI8/ni4SqOISK60dijIw8B3zeyktHYy8N0sxCYiUhDWrFkD5OVQkKLK2aWl8PDDcO65sHEjnHoqTG0XEwuKSKFobY/1dcChwFcBj9YZMAC4PL6wREQKx9KlSwEYPnx4wpHso+hydqdOoZi+7DL41a9g0qQwJd+PfxwKbxGRbGptj/XPgCHuXuHuB0etApgYf2giIvnP3Vm0KMxeN2rUqISj2UdR5uyOHeHuu+FnPwuzhfzsZ3DSSbBkSdKRiUh719rCeqW7b2tgfZu//zSz0WZ2p5ldYWZ3mdmx0foyM7vdzP7dzO42sx5t3ZeISFyqq6vZtGkTPXv2pH///kmHk6loc7YZfOc74YqMBx0EM2bAUUfBT34Ce/YkEZGIFIPWFtbfyFxhZucDcUzJ/wDw7+7+C+CnQGo20ruAf7j7NcCjwK9i2JeISCzeeOMNIPRWm1nC0eyj6HP2xz4Gc+fCl78Mu3bBNdeE2UL++ldwb/71IiKt0drCuo+ZXQdgZgeY2T3APcCrMcQyGkhdtmwn0NPMBhEuvftYtP4Z4Fwzq4hhfyIibTZ9+nQAxo4dm3AkDVLOBvr2DReRefppGD4cFiyAc84JJze+8kqSkYlIe9OqwtrdzwGeNrPbgbnASOAYoKEre7XWHwgXMugOfAn4JjAeWO/uO6L9fwgsA06OYX8iIm328ssvA/DRj3404Uj2pZxd34QJoaj+j/+A3r3h+edh7Niw/umn1YMtIm3XZGGdMUXTSWZ2EtAJWA4sAb4PDAOujyGWbwAfEHpStrv7VGAwsCHjeduAQZkvrq6uprKysrZNmTIlhpBERBq3Z88eXnrpJQBObOWEyVOmTKnNV0C/OOJRzm5ely7w7W/D22/DVVdBt27w7LNwxhnhCo733Qc7d+YkFBEpIC3N2eZN/ItuZjMIX/dtBmoaexow0N0773+4YGb9gE8TEvNV0e1jgfPdfWza814DHnL3n6e/vrKy0mfNmtWWEEREWuWZZ57h9NNP5/DDD2fBggX7vR0zm+3ulW2NRzm79TZtCjOI3HYbvPdeWNerVxiTPXkyHHFEsvGJSP5pKmc3N4/1D4FZ7v5eMzuI42vF3wAXuPtmM3PCiTDfAnpnPK87sCaG/YmItMmjjz4KwKc//emEI6mlnN1KvXuHExqvuAIeegimTAnjrv/zP0OrrIQLLoDPfQ6GDEk6WhHJd00OBXH3vzaXoKPn/aEtQUQ9H8e4++Zo1Y+ArYSvLweZWVn0vFJgKPBCW/YnItJWu3bt4uGHHwbgM5/5TMLRBMrZ+6+sDC6+GGbOhDlz4NJLoUcPmDUrTNs3dGiYC/vOO2HlyqSjFZF81dpZQbJlI/CBmQ1OW7cBeB14HDgtWjcBmOruq3Mcn4hIPX/84x9Zv349Rx99dGrMXTFp1zn72GPhrrvg3XfDVRwnTQpjs198MVzRcejQMGXf9dfDq69CTWODbkSk6LT2kuZZ4e41ZnYecL2ZzQb6A1e6+1YzuxS42cyOIVzUYHKSsYqIuDu/+MUvAJg8eXI+zl+dVcWSs8vKYOLE0LZtg0cfhUcegaeeCnNjz50LN94IAweGmUUmTIDTToPBg5vftoi0T02evFhI8uVEGBFp//7yl7/wyU9+kvLycpYtW0bXrl3btL24Tl4sJIWcsz/4IEzV99hjoWUODTnssLoi++Mfh/LyRMIUkSxpy8mLLdn4GOD/AXMIX/lVtXWbIiL5as+ePVx99dUAXH311W0uqnNNObvtunSBs84K7Y47wtzYzz4LzzwTCu7Fi0O7887w/JEjQ4H9sY+F5aGHhkuui0j7E8dQkP8LTHD3L5nZF8xsnbtrFlARaZd+/vOfM3/+fA455BAuvfTSpMPZH8rZMTKDI48M7fLLYc+eMO76mWfguefCyZBvvRXavfeG15SXhyJ73Dg4/nj4yEegZ89k34eIxKPNQ0Gis74PdPfVZtYB6ObuW2OJrhUK+WtFESkMc+fOZdy4cezatYunn36aCRMmxLLdXA4FUc7OrT17wljsf/wjtJdegnXr9n3eyJGhyD7++DDF35gxcMABuY9XRJoX61AQMzsb6OPuD5rZCYQrbi2EcEILYcolEZF2ZfPmzXz2s59l165dXHLJJbEV1dmmnJ2s0tK6gvmKK8Jl099+OxTYr74a2ty5db3aDz4YXldSAqNHw9FHwzHH1C0HDNAwEpF81qoeazP7MXA5MM3dJ0Xrfgq87u4PZSfElimW3g8Ryb0PPviAT3ziEzz33HMcd9xxvPTSS3Tp0iW27Werx1o5uzDs3g1vvBGK7FmzwnLBAti7d9/n9utXv9A+6qhwsmSBDfUXKWhx9lh/FDgI+ErauieAB4BEk7SISDbs2bOHCy64gOeee46BAwcyderUWIvqLFPOLgCdOoVx1h/5SN26HTtg/nyYNw9efz20efNg/fpwouSzz9bfxrBhoYd79GgYNarudr9+uX0vIsWutYX102mXr02J49K4IiJ5Z/v27Zx//vlMmzaN3r1789RTT1FRUZF0WK2hnF2gDjgATjghtBT3MLVfqtieNy8U30uWwPLloU2bVn87/frVL7hHjgyzkhx8cBimIiLxam1h/Y6ZXQMMNrMLgPOBzwDfjD0yEZEErVmzhs985jO88sor9OvXjyeeeIIjjzwy6bBaSzm7HTELV30cOhQ+9am69Xv2hHHbixbBm2/WtUWLQg/3iy+Glq6kBCoqYMSIUGint2HDVHSL7K9WzwpiZuOArwJDgbXAb939mfhDax2N1xORuDz33HNccMEFvPfee1RUVPDkk08ycuTIrO0vm7OCKGcXL3dYtap+ob1kSWgrV4bHG9KxYyi604vtESNg+PBQdHfunNO3IZJ3Yp0VxN1nADPSNl5iZqe5+7NNvExEJO/t2rWLG2+8kZ/85CfU1NRw6qmn8tBDD9G/f/+kQ9tvytnFywyGDAntjDPqP/bBB/DOO3WF9tKl9YvupUtDe+KJfbc5eHAYSjJ8eFim3x44EDp0yN17FMk3TRbWZvZnoFsz2+gHrAeUpEWkYM2YMYOLL76YhQsXYmZce+21/OAHP6CkpCTp0FpMOVtaqksXOPzw0DLt3Nlw0b1sWSi6V60KLXN4CYTe7IqKuoI7swDv3Tvrb00kUc31WL8ePWcxUANcQrgMbvr3dwcDfbISnYhIlq1Zs4Zrr72W++67D3dn5MiR/PrXv+bjH/940qHtD+VsabOyMjjiiNAy7dkTiutly0J75536t6ur6y7p3pAePerGiQ8bVnc71QYNCkNRRApVc7++PyOMw94MYGaHuPv1mU8ys0ezEZyISLZs2bKFW2+9lZtvvpkdO3ZQWlrKt7/9bb7//e9TVlaWdHj7Szlbsqq0NPRADx/e8OPbt0NV1b4Fd+r21q1hJpP58xt+fYcOYahJQ0V3quny75LPmiys3X1Lxqr+Zlbi7rXT1pvZeOAEREQKwMaNG7ntttu47bbb2LIlpLiJEydy0003MWLEiISjaxvlbElat25w5JGhZXKHjRthxYrG29q1oUd85crG99GzZxg3ftBBoQgfPLjudmrZp4+uUCnJaO0XLn8A5pnZP4BtwGjgTGCfHhERkXzy9ttvc+edd3LPPfewbds2AE4++WR++MMfctJJJyUcXdYoZ0veMIO+fUMbM6bh5+zeHcZvN1Z4L18OW7aE1livN4Qx5I0V3anbAwZo2InEr1W/Uu7+tJmdAXwDOBrYDExyd32tKCJ5p6amhqeeeorbb7+dJ554gtT0oqeffjrXXXcd//RP/5RwhNmlnC2FplOnpoeauMOmTaHIXr06FOENLbdsCXN7v/124/vq0CEU15kF+ODBYXaTgQPDmO9evdT7LS23P9PtrQG+l7pvZh3M7CPuPjvWyERE9tOSJUt44IEHeOCBB6iqqgKgc+fOXHjhhVx22WV8JP3a0e2ccra0J2ZhmEefPnDssY0/b/v2pgvvVavgvfdgzZrQXn218W117lxXZKcX3Jnr+vZVAS6tLKzN7N4GVpcDe4HzYolIRGQ/VFdX88gjj3D//ffz8ssv164fOnQol156KRdffDHl5eUJRph7ytlSrLp1g8MOC60xu3eHMd2rV9cV26tWhUJ77dq6tm1bOCEz+h+9UZ06hR7wxorw1O1+/TTXd3vW2h7rQwlzn6Zfr+kgoIHZLEVEsmv16tX88Y9/ZOrUqfz973+npqYGgK5duzJp0iS+/OUvM378eDoU718x5WyRRnTqFGYfGTas6edt315XZKcX3Zm3t2ypGwvelI4doX//xovvVDvwQI0BL0St/ZF9zd3rzU5pZqMIl8sVEcmqmpoa5s6dy7Rp03j88ceZPn167WOlpaWcccYZXHjhhUycOJFu3Zq7TkpRUM4WaaNu3eou7d6UnTtbVoBv3FjXS96UDh1CcZ1ZcGcW4QMGhH8SJD+09uTFhqZ83wV8Dbg6joDMbBzwcWAp8A/gfeBmwhntvYGr3H1rHPsSkfy3ceNGnnrqKZ544gmefPJJ1q1bV/tYly5dOPvss5k4cSKf/OQn6dWrV4KR5h/lbJHcKStr+sTLlF274N136xfdmcNP1q4NF9t5993Q5sxpept9+zbe852+vnCn6C8crR1jvYz6XymWAP2BB+IIxsz+GTjE3dNPtPlvYJq7/87MzgJ+BVwYx/5EJP9s2bKFF198keeff57nn3+eOXPm1A7xADjooIM4++yzOfvssznjjDPo2rVrgtHmN+VskfzTuXPLhqDs2QPr1u1bcGcW4evWwYYNoTU1BSGEOcCb6v1Ote7ddSLm/mrtUJCngYeoS9Q1wDp3f6utgZjZycAXgdPT1g0CvgD8S7TqGWCqmVW4e1Vb9ykiydu0aRMvvfRSo4V0aWkp48eP5+yzz+ass87iiCOOwJTxW0o5W6RAlZaGKQAPOqjp5+3dG3q3Mwvwhorw1BzgixY1vc2uXZvu+U613r1VgGdqbWH9G3f/e/oKMzvMzEa4+9I2xnIL8DJwu5mNAG4EhgHr3X0HgLt/GPXAnAxUtXF/IpJje/fuZcGCBcyYMYPp06czY8YMFmVk+I4dO3LiiScyfvx4xo8fz4knnqhe6f2nnC3SzpWUhHHWAwY0fuEdqLvyZUPDTjLXvf8+LF0aWlM6dw77TRXcgweHq2Kmrow5ZEh4rLQ03vecz1pbWJ8F/D1j3TrgOaCJH2fTzOyw6PX/7O6vm9mVwDTgR8CGjKdvAwbt775EJDfcnbVr1zJ79mxmzpzJ9OnTeeWVV9i+fXu953Xu3JnKykoV0tmhnC0iQP0rXx51VOPPc4etWxsuwDOL8G3bwtUwly9vfHupC/GkCu70ojvVBgwI/yC0By0qrM3sMkLC7G5mVzXwlKfbGMcRwEZ3fz26fwfwfcLXlzszntsZ2JO5gerqaiorK2vvT548mcmTJ7cxLBFpCXdn+fLlvPbaa/Va+omGKRUVFYwbN44TTzyRcePGceyxx9KpCE9pnzJlClOmTEnd7RfntpWzRWR/mYWx2D17wujRTT/3/ffrF9urVsHKlXXLlSvDyZepC/HMnNnwdkpKQs/2sGF1J4AefHDd7QEDkp/7u6U521KX+G2OmR0NfAX4z4yHdrv72v0JMm3bnwOmuHuvtHVzgP8l9IiMTFu/BPi+uz+Uvo3KykqfNWtWW8IQkRbYvXs3ixcv5o033mDu3Lm1RfSmTZv2eW7Pnj0ZM2YMY8eOZdy4cYwbN44BAwYkEHV+M7PZ7l7Z/DNbtU3lbBFJXOpCPKlCO9XSi+/33mt6G126QEVFXaE9fDiMGhUK/6FDc190N5WzWzwUJPq672p336fnwcyGunszU6I36XWgp5n1d/dUF9eHwBpgkJmVuftOMysFhgIvtGFfItICNTU1VFVVMX/+fN544w3eeOMN5s+fz+LFi/nwww/3eX55eTnHHXdcvXbwwQfrRMOEKGeLSD5oyYV4PvggzOu9fDksWwbvvFO/rV8fTrhs6KTLAw4IRfbhh4dCe/RoOOaY0OOdxJ+fJgtrM/sO8Fd3Xxit+jczyxwFY8ApwPj9DcLdF5nZE8Ak4E4z6wV0Ah4EzgZOAx4HJgBT3b2ZadVFpKVqampYtWoVixYt4s0336wtoOfPn8/777+/z/PNjBEjRnDUUUdx9NFH1xbRgwcPVhGdMOVsESlEXbrAIYeE1pBt2+oX3EuXhiJ74cIw3eBrr4WWrm9fqKwMbexYOOmkMLwl25rrsT4YSL/iwnFAV+Bd6qZvMiCO73a/DNxmZp0IPRwXRmeUXwrcbGbHRPvRIDyR/fDBBx+wZMkSFi1aVFtEL1q0iMWLF7Njx44GXzNw4ECOPPJIjjrqqNrl4YcfzgEHHJDj6KWFlLNFpN3p3h2OPjq0TBs3wptvhrZwYWivvRamIHzyydAgjOMeOxbOPBO++MXGi/i2avEYawAzOw6Y7+67M9YfntZDkgiN1xOpm4lj6dKlvPXWW/WK6GXLltHY571///6MGjWKUaNGceSRR9a2fv1iPadOGpGNMdbRdpWzRaTouIex26++GtpLL8GMGZA+ivGss+CWW5o/QbMhsYyxjhzl7vU626Nplw4BEk3SIsWipqaGtWvXsmTJEpYuXcrSpUvr3W6s97mkpIRDDjmktoAeNWoUo0eP5rDDDqN37945fheSI8rZIlJ0zMJJjUOHwmc/G9Zt3QrPPw9/+ENo06bB9Omhd7u5y9C3Rkun2/snwteHp0eT/afrQzjr/LH4whIpbjU1NaxevbrBwnnp0qXs3Jk5o1mdfv36MWLECA499NB6BfQhhxxSlNPaFSPlbBGR+nr0gHPPDe273w0nO27ZAk8/DV//enz7aWmP9WbgfsL4vY9lPLYbuCe+kESKw+bNm1m2bBnLli3jnXfeqbesqqpi165djb62vLy8tngeMWJEvdu9evVq9HVSNJSzRUTSrFwJzz0Hjz4Kjz8e1nXvDhMmxLufFhXW7v6GmY0Hznd3JWSRFti1axfLly9vtHhuaN7ndOXl5bXFcvrykEMOUfEsTVLOFpFiVVMTiuj582HuXJgzJwz3WJb23Z1Z6Lm+6ab4T2JszTzWW2ikl8PMPu7u/4gtKpECsHfvXtasWUNVVRVVVVW1BXOqeF69enWjJwsClJWVcfDBBzN8+PB6y1Tr3r17Dt+NtDfK2SLSXr3/ft3FZZYvD9PvLVkCb70Vbn/wwb6v6dEDTj4ZTjsNJk2CwYOzE1tz81j/GejWzDaMcHnbA+MKSiQf7Nmzh1WrVrF8+XKqqqpql6nbK1eubPBCKSkdOnRg2LBhtYVyeuE8fPhwDjzwQM37LLFSzhaRQlZTAxs2hCs1vvtuaOmXSl+xIs7SMTMAACAASURBVLQNG5rezoAB4aIxY8bAsceG5ejR0LG1U3bsh+Z2sZgwVm81dXOgZjLgk3EGJZILu3btYuXKlfWK5vTiefXq1dTU1DS5jQEDBlBRUcGwYcP2KZyHDBlCaWlpjt6NCKCcLSJ5pqYGNm0K80qnWnrRnL5ct67+lHiN6dQJhgypayNGwKGHwsiRYdmjR/bfV2OaK6xvAra5e+NnUQFmprPLJe/s3LmTFStW1OtlTl+uXbu2yaEaZsbgwYOpqKioLZ7Tl0OHDqVLly45fEcizVLOFpGs2rs3XJQlvVCurob33tt3XXV1uBz53r0t337v3jBwYOh1Ti0HDQqXRB86NBTSBx4IHTpk7z22RZOFtbuvb+F2PoPOMpcc2759O8uXL99niEZquW7duiZf36FDB4YOHVpbKGcWz0OGDNH0dFJQlLNFpDV27w5F8oYNoaXfTm8bN4YC+b33wu1mvszdR8+eUF5e19KL5szbnTtn573mSptGm5jZUOBfgMtQkpaYbdmypdHxzVVVVWxoZpBVx44dGTp0aIO9zRUVFQwePJiOuRhwJZInlLNF2if3MCdzQwVxY4Xyhg2wbdv+7a9Pn/qFcmY78MC62/36haEbxaLVVYWZlQDnApcApwPrCGP6RFpl586dVFVV1c6kkTmP8+bNTf9ade7cubZwbqh4HjhwICUlJTl6NyL5STlbpDC4w/btYTxyY23jxsYfa81wi5SSklAk9+0bWvrthlp5eVjq9KHGtbiwNrODga8B/wz0BZYDnwKeBE7MSnRS0Pbu3cuqVav2KZpTbe3atU2+vqysrNHe5mHDhtG/f3865OsgK5GEKWeL5F5LiuPGCuXNm1t24l5jundvWXGc/niPHvk7VrlQNVtYm9n5hJ6OU4GlwC3AfcA17j4tetrLWYtQ8pa7s379+kYL5+XLlzc5HV3Hjh3rTUeXPqNGRUUF5eXlmo5OpJWUs0XabufOUPSmCt/MZVOtLcXxAQeEk/da2vr0qbtdTMMt8llLeqx3Ei6B+zLwNXd/C8DMGp9OQdqNVPG8ZMmSfdrSpUvZ1swArYEDBzY4j/PBBx+sMc4i2aGcLUIocDdvbrgwbqpo3rgRdjU5r07Tysr2LXpb2gr9xD1pQWHt7o8Dj0cnvVxiZhXAo0DtCBszM29q3jLJexs3bmyweF6yZAlbtmxp9HU9e/ZssGg++OCDqaiooKysLIfvQkSUs6U9cQ9X2WuqCG7ssSb+dDWrU6dQGKeK4/QiubmCWcVxcWvNJc1XANeZWUfgPOAgM7sR+DVwBjAlOyFKXNydNWvWsHDhQhYsWMDChQtZuHAhb775Jhs3bmz0dT169ODQQw9tsPXt2zeH70BEWko5W/JNauaK9evrt9Rcx+ktvWDes2f/9mcGvXrtWxw3tMxcV1YWXi/SWq3+Ht7dPwT+APzBzEYAVxDG8ylJ5wl3Z8WKFbWFc3rbunVrg6/p2rVro8WzxjqLFC7lbMmWHTsaLoobWrd+fZjebX/GH6fGHTdVCDe07NlTJ+ZJ7rVpgKu7LwUuN7PZMcUjrbRnzx4WLlzI3LlzmTNnDnPnzmXu3LmNDt/o06cPRxxxBIcffnjtcvTo0QwcOFDFs0g7p5wtzdm5M1xWet26uktMp5aZ67Zvb/32e/QI8xo31FJzHqfPXqGhFVJoYjlzzN3vj2M70rTdu3czd+5cXnnlFebMmcOcOXNYsGABu3fv3ue5/fr148gjj+Twww+vV0Sr91lElLOLz/vvw6pVsHp13XLNmn2L50a+1GxQp07hQiBNFcnprW9fFcnS/mlKhjy2atUqZsyYwfTp05kxYwazZ89mVwOnKo8YMYJjjz2WMWPG1C4HDBigAlpEpJ1zD+OQV62qXzhnFtHNXG+rVmkp9O9f1wYMqL9Mv92rl8Yhi2RSYZ1Hli1bxt/+9jf+9re/8cILL7B69ep9njNq1CjGjh3LRz7yEcaMGcPRRx9Njx49EohWRERyYetWWLYMqqrCMv12VVXLLkvdqRMcdBAMHly3HDQoFMnpxXPv3iqWRdoirwprM+sEvApc7u7Pm1kZcDOwDegNXOXurfiiKr+tX7+eJ598sraYrqqqqvd4z549GTduXG0bO3YsvXv3TiZYEZEGFFvezpaNG2Hx4tDeeiu0VBG9aVPTr+3eHYYM2bdwTl/27auCWSQX8qqwBq4EKtLu3wVMc/ffmdlZwK+AC5MILC6LFy/mz3/+M3/+8595+eWXqampqX2sd+/ejB8/nlNPPZVTTjmF0aNH65LdIpLv2n3ejos7rFwJr78OCxfWFdKLF4dZMxpTVgYVFXDwwaGlbqeW6mUWyR95U1ib2UeBtcCm6P4g4AvAv0RPeQaYamYV7l6VSJD7acmSJTz44IM8/PDDLF68uHZ9aWkpp512GmeccQannnoqxxxzDCUlJQlGKiLScu05b7fVrl2wYAHMm1e/Ndb73LUrjBwJhx0W2siRMHx4KJwPPFCFs0ihyIvC2sy6ApPc/dtmdn20ejyw3t13QJiL1cyWAScDVYkE2gqbN2/mwQcf5IEHHmDmzJm163v37s0555zDueeey5lnnqnx0SJSkNpj3m6LVatg+vS69tpr0MCETfTtC8ccA0ceWVdEH3ZYGLKh4lmk8OVFYQ1cBfw0Y91gYEPGum3AoIY2UF1dTWVlZe39yZMnM3ny5DhjbJEFCxZwxx13cP/997Njxw4AunXrxsSJE/nSl77EKaecQseO+XLYRSQpU6ZMYcqU2mu09Esylv3UprydLzl7f61aBc8+C08/Dc8/H2beyHTYYXDssaGQTrVBg1RAixSilubsxCs8MzsTmOXu72U85MDOjHWdgQYvblpeXs6sWbOyEGHLTJ8+nRtuuIGnnnqqdt1pp53GxRdfzLnnnkvXrl0Ti01E8k96IWlmTYywzT9x5O2kc3Zr7dkDL7wAf/5zKKYXLar/eK9eMHYsnHhiaCecENaJSPvQ0pydeGENfAcYkzbncm/gUeCW6Ha67sCa3IXWvFmzZnHdddcxbdo0IFwa/Mtf/jKXXXYZhx9+eMLRiYhkRUHn7ZbatQuefBIeeSQU1Onjo7t1g/Hj4fTT4bTTYPRoXT5bRPKjsP4CoUcjZTrwbWAmcKWZlbn7TjMrBYYCLyQQ4z7Wr1/PNddcwz333IO7061bNy6//HK+/e1v06dPn6TDExHJpoLM2y31+uvw61/Db34TpsFLGT0aPvMZOPvs0DtdWppcjCKSnxIvrN29Ov2+me0Fqt19hZk9DpwGPA5MAKa6ewMj2XLrt7/9LZdddhkbN26ktLSUyy+/nKuvvpq+ffsmHZqISNYVYt5uTk0NPPoo3HwzzJhRt/7oo+Hznw8F9ejRycUnIoUh8cK6GZcCN5vZMcAAINEzW7Zv385ll13GfffdB4Qx1HfccQejRo1KMiwRkXySV3m7OTU18OCD8OMfh/mkAXr2hC9+ES6+GI47Ltn4RKSw5F1h7e4Vabc3AZckF02dFStWcM455zB//nzKysq49dZbueSSSzCd3i0iRS5f83ZzXnwRrrgCZs8O94cOhSuvhIsuggMOSDY2ESlMeVdY56N58+Zx9tlns3btWkaNGsXUqVN1YqKISIF6//1QQN99d7g/eDD86Eehl1rjpkWkLVRYN2PhwoVMmDCB9evXM378eB555BF698486V1ERArB/PkwcSIsWRKK6O9+F/7t38KVD0VE2kqFdRNWrVpVW1SfddZZ/OlPf6Jz587Nv1BERPLOU0/B+efD1q1w1FHwwAPhoi0iInFRYd2IXbt2MWnSJNauXctJJ53EI488oqJaRKRAPfkkfOpT4UIv558P990HZWVJRyUi7Y2ms2/ENddcw8yZMxk6dChTp06lTBlYRKQgvfZaGP6xZw9885vw8MMqqkUkO1RYN+DVV1/l1ltvpaSkhN///vf069foJeFFRCSPbd0Kn/sc7NgBX/kK3HqrrpAoItmj9JLB3bn00kupqanhiiuu4IQTTkg6JBER2U/XXQdvvx3GUv/ylyqqRSS7lGIy/OlPf2L27NkMGDCAG264IelwRERkPy1ZAnfdFYrp++6DLl2SjkhE2jsV1mncnR/84AdAGGPdVfMviYgUrNtugw8/DENANPuHiOSCCus0L7/8MvPmzePAAw/kkksK4sJhIiLSgJ07w3R6EK6uKCKSCyqs0/zyl78E4KKLLqKLvjMUESlYzz8fTlwcMybMWS0ikgsqrCM7duxg6tSpAOqtFhEpcE8+GZaf+ESycYhIcVFhHZk2bRo7d+7khBNOYPjw4UmHIyIibfDKK2F58snJxiEixUWFdeSRRx4BYOLEiQlHIiIibbF3L7z+erg9ZkyysYhIcVFhDdTU1PDEE08AcN555yUcjYiItMXbb8P778NBB4Gu7yUiuaTCGnj99dfZuHEjQ4YMYeTIkUmHIyIibbBwYVjqpEURyTUV1sBzzz0HwKmnnoqZJRyNiIi0xZtvhuWoUcnGISLFR4U18Le//Q2AU045JeFIRESkrRYtCsvRo5ONQ0SKT9EX1u7Oyy+/DMDJOn1cRKTgpQpr9ViLSK4VfWG9bNkyNm7cSHl5OcOGDUs6HBERaQN3FdYikpy8KazNbKKZLTOzDWZ2m5l1jNb3M7NfmdlNZnaLmZXGud9XX30VgBNOOEHjq0VEWiipnN2ctWvDFRf79NGMICKSe3lRWJvZUOA8YBJwBfA14FvRw78H7nX3q4C3gB/Gue9XoqsIHH/88XFuVkSk3UoyZzcnfXy1+kpEJNc6Jh1AZBhwkbt/CMw2syOBU8zsH8Bh7j4zet5jwGIz+5G7vx/HjlM91iqsRURaLLGc3RwNAxGRJOVFj7W7vxgl6JTVwArgVGBZ2vNWRzcrY9ovb7zxBgBjdHkuEZEWSSpnt4Sm2hORJOVFYd2A44FfAoOBDRmPbQMGZb6gurqaysrK2jZlypRmd/Lee++xefNmevbsyYABA+KIW0SkRaZMmVKbr4BCHw2ck5zdEuqxFpFsaGnOzpehILXMbARQ7e7zzMyBnRlP6QzsyXxdeXk5s2bNatW+FkUZeNSoUTpxUURyavLkyUyePBkAM1ufcDj7LZc5uyWWRf3lhxwS+6ZFpIi1NGfnVY+1mZUAXweujlatBnqnPW5AV2BNHPtLFdajdRUBEZFWy3XObk5NDaxcGW5r9lQRSUJeFdbAd4Cfu/uu6P7TwGFpj1cAu4FYujmWLFkCwMiRI+PYnIhIsclpzm7Ou+/C7t1QXg4HHJCLPYqI1Jc3hbWZXQPMAcrMbLiZXQRsBtZGZ5wDnAXc7u6749jn6tXhvJohQ4bEsTkRkaKRRM5uzvLlYaneahFJSl6MsTaz7wE3Zqxe5O73mtnngevNrAroBlwb137XrAnfTg4ePDiuTYqItHtJ5ezmpIaBDB2aqz2KiNSXF4W1u/8I+FEjjy0HLs7GflM91oMG7XPCuoiINCKpnN2cTZvCsm/fJPYuIpJHQ0GSkOqxHjhwYMKRiIhIW23dGpY9eiQbh4gUr6ItrHfu3MnOnTvp1KkT3bt3TzocERFpIxXWIpK0oi2sN2/eDEDv3r01h7WISDuwbVtYqq9ERJJStIX1pmgwXq9evRKORERE4rArmvSvS5dk4xCR4lW0hXV6j7WIiBS+3dGkfp06JRuHiBSvoi+s1WMtItI+qLAWkaQVbWGdGgqiHmsRkfZBhbWIJK1oC+stW7YA0EOnj4uItAsqrEUkaUVbWG/fvh1AU+2JiLQTKqxFJGl5ceXFJAwfPpxzzz2Xo446KulQREQkBmPHQkkJ6JpfIpKUoi2sJ02axKRJk5IOQ0REYnL99UlHICLFrmiHgoiIiIiIxEmFtYiIiIhIDFRYi4iIiIjEQIW1iIiIiEgMirawnjJlStIhtIrizb5Ci1nxZlehxVss9HOJl45nvHQ841WIx1OFdYFQvNlXaDEr3uwqtHiLhX4u8dLxjJeOZ7wK8XgWbWEtIiIiIhInc/ekY4iFmVUDy1vxkn7A+iyFkw2KN/sKLWbFm125jHeYu5fnaF95YT9ydkqh/R7lOx3PeOl4xitfj2ejObvdFNYiIiIiIknSUBARERERkRiosBYRERERiUHRFdZmVmZmt5vZv5vZ3WbWI+mYmmJmE81smZltMLPbzKxj0jG1hJl1MrN5ZjY+6VhayszGmdl3zOw8M+uXdDyNMbPRZnanmV1hZneZ2bFJx5TJzE43s5lmVpG2Lm8/e43EW5CfvfYqn39/CkmhfTbzWWM5wsz6mdmvzOwmM7vFzEqTjrUQmNnJZrbYzDaZ2a1p6wvq97PoCmvgLuAf7n4N8Cjwq4TjaZSZDQXOAyYBVwBfA76VaFAtdyVQkXQQLWVm/wx80t1/7u5/cvd8PFki5QHg3939F8BPgd8mHE89ZtYf6AackPFQXn72Goq3wD977VVe/v4UkkL7bOazZnLE74F73f0q4C3gh4kEWUDMrCfwMWAc8CXgX8xsQvRwQf1+FtXJi2Y2CFgG9Hb3HdF/l1uAI9y9KtHgGmBm/wRMd/cPo/s3E2I9J9nImmZmHwVGAdcDX3X355ONqGlmdjJwHXC6F8AHwszeBz7i7ovMrByY5+6Dko4rnZl1APYCB7t7Vb5/9hqItyA/e+1Vvv/+FJJC+2zmq8ZyBHAj8EgqJ5vZYGAx0N/d308q3nxnZmXAB6m/wWY2G/gO4dgV1O9nsfVYjwfWu/sOgOgDsQw4OcmgGuPuL6Y+tJHVwIqk4mkJM+sKTHL3e5OOpRVuAd4EbjezaWb2saQDasYfgHvMrDvhP/tvJhzPPty9JmPVePL4s5cZbyF+9tq58eTx708hKbTPZr5qIkecSjh+qeetjm5W5jC8guPuO9OK6q6Ev8nPU4C/n8VWWA8GNmSs2wbkVW9fE44Hfpl0EM24ijA8oSCY2WHAGOC/3P0y4FlgWvSVab76BvAB8Cqw3d2nJhxPS+izJ21R6L8/+UzHNh6pHKHj2QbRNwFPAKXAARTg8Sy2wtqBnRnrOgN7EoilVcxsBFDt7vOSjqUxZnYmMMvd30s6llY4Atjo7q9H9+8ADPh0ciE1q4wwrvoh4Na0cWj5TJ89aYuC/f0pADq2bZSRI3Q822Yp8D/AacDPKMDjWWxnua8Gemes6w6sSSCWFjOzEuDrwNVJx9KM7wBjzCx1vzfwqJn92N1vTi6sJnUk7XPg7jvNbAnQN7mQmvUb4AJ332xmDvzWzCryfPyePnvSFgX5+1MgdGzboIEcsRoYmfa4AV3R8WwRd18L3GtmNYRJEF6kwH4/i63H+jlgUDRInmgKnKHAC4lG1bzvAD93911JB9KMLwDHprU1hDOl8/kr9NeBnhlDPz4EFiYUT5OiaQCPcffN0aofAVuB0clF1SL67ElbFOrvTyHQsW2bzBzxNHBY2uMVwG5gVo7jKnSzCP+kFNzvZ1EV1u7+LvA44SsGgAnA1LSTC/KOmV0DzAHKzGy4mV0Ufe2Ud9y92t1XpRrhzPNqd9+adGyNcfdFhPFckwDMrBfQCfhLknE1YSPwQXSmecoGwpROecPqvrYwyP/PXma80bqC+ey1d/n++1NICu2zmc8ayhHAZmCtmR0ZPe0s4HZ3351UnIXAzA4wszFpqz4B3FaIv5/FNhQE4FLgZjM7BhgATE44nkaZ2fcIU/ekW1RgM24Ugi8Dt5lZJ8J/whdmnO2dN9y9xszOA66PpiPqD1yZT/+8pM1WAvBVM7vD3avJ089eQ/ESvtrVZy+/5OXvTyEptM9mPmvq77OZfZ6Qo6sI84Zfm+v4CtCRwF/MbBEwA3jF3VMdXAX1+1lU81iLiIiIiGRLUQ0FERERERHJFhXWIiIiIiIxUGEtIiIiIhIDFdYiIiIiIjFQYS0iIiIiEgMV1iIiIiIiMVBhLdIIMzvNzF42s/HNPO9IM/tfM7suxn3/LK5tiYgkwcyONrPfx5kbG9nP42Z2QTb3kS1mNsHMpjf3d0YKRzFeIEYKjJmNBn5MuOrg8cDRwP9x999kedcrgWNa8LwVQG+gJI6dRpds3RPHtkREAMzseOAnwEnAU8AywsU2SoG73f3JLOx2KdCTmHJjE+4C5mdjw1E+ngj8knCF2+nRQ52BzwInuHtVG3bR0r8zUiBUWEsh+BPwOXefZ2YdgP/OxU7d/S0z29CC5201s7Ux7vpk4O8xbk9Eipy7v2pmDwOj3P2TqfVmdg7wRzP7jrvfFfM+d8ScGwEws08Bb6QKWnf/a9z7SHH3PcDvzOybwPPu/r20OB7dn22a2YHAKe7+O3dfbGbrYwpX8oCGgkheM7M+wEjgfQiX9AauB3J1ydCamJ/XEqcAz8W4PRERgA8zV0SXjb4B+IWZDcjCPmPN1WZ2EKGHOtf2OXbAM8Dq1mzEzDoD9wNlaavj/PshCVNhLXnN3TcCi4FpZjYuWrcceMXMepnZjWb2ipmda2bvmtmbZnZc6vVmdqGZXWtmz5vZv6etrzSzG8zsd2b2iJn1jNZ3MbM7zex6M/s14WvMVjGzf4ti+pyZHW9mf4q2d2MU4+/MrDxav8XMvp+xiU7uviva1v8zs2+b2Utm9tNWH0ARkeb9mjAk5HwAMzvEzH5oZveY2bNRMYuZdY5y2bVm9qKZfT61ATP7hpldZ2YzzOyyxnbUUE42s0lRzvyqmS0xs5+a2SAz++8on75oZqdGm/gCcBBwlZmdF50L84yZfSVtH5eY2Y/N7H4zu8/MuppZTzP7vpnNMrNTzWyxmb1jZhX7c8CiYzI+6tHGzCaa2W1m9h9mNs3MukXv4TYzeyz6e7M2OsaHA+eb2bfTNtk7+lu0zcwuaWSfHczsJjObHO3DzeyPqZ+P5Al3V1PL6wYcASwn/Fd/L3BgtL6EkKQ2Esa69QKeBhYQ/mk8Efhe9NzehF7vc4AewMNp238J+Fl0+xfA5Oh2f2AvIXk2F+P/EHp9DiSMxTsgWt8ZeBJ4FBhE+IOwB/g+cADwkeh+6vmp8eOp931XdLs7cHXSPws1NbXCbcBXgVWNPLYO+M8orz4ClETrHwJ+H92+BZgY3b4IqI5uXwh8Kbo9MsrVR0X3/we4IbrdWE7uA2wGbgVOACYAPweuj557JfBYWqwOVES3y4A3ga9G988EHkp77uPArwhDXz8T/b2YGL3PacCPW3jsno/+VvwU+BlhvPWEtMfXA0Oj2wsIf5M6A9cCq4CxwMVR3n8+FW/0/CrgP6LHLgRWNBLD14DHo9v9o+N8RtK/V2r1m8ZYS95z9wVmdhTwI+AbwKfMbIKHMdfVwFZ3nwpgZjcA/wCGE/6I1JjZt6JNPUlI4OcAfdPWvxm9thtwafRa3H2dma1sRaj9gNuBi919R7SNXWb2LrDM3ddE+3kXeMHD+MM5hITfj3AS5DmEPwIAO4Avm9lbwN2EfypERLLBCcXmCcAw4JtmBrAB6GnhztcJnQIA9xGGQkDItbMycu1A4I2MfXyVBnKyu280sy3Ao+7+CoCZrQDeN7OBhB7ebg0G7b7TzN5LW3VZtN2Ue4D/jdZvJvy9eCTax6wozpZ6zqMx1mb2O0JnTsr57r7CzD5GKKi7Rfl/DfC2u88EZkavbWjbj0V/E15tIqbjgJ3R+15nZhupP6RE8oAKa8lrZlYCDPFwksq/RsMz/gpMIfQAZFoQLXsBQ4EH3P2haN2t0TavIhS6t2bs6yOEhLhtP8MdCXwcOASYl7Y+c4zh3toH3GuiJJsaltXXw/AX3H2ZmX0VuAO4Avgi8B4iIjEys16Eb9sWEPLmhgbyY39Cj2onAHffS+gMIHrND9z95eh+vdemaTAnR5z6uXIVcBWwBJgBHNzEW0h/3aHA39Luv00Y5nIg++biD9n/GUteI/wDkrLZzG4Cfk+IPVU9Z76v5nxI47XZE8DdZtaF8LPYSej9ljyiMdaS72rH/QG4+zzgGmBUI8/vREhiS4G1wKT0B81sbLT+HDMry1i/Nbo7ej9jfRm4DXgwSnytYmaDSDsRxswGE2ZEGUk4mfGhRl4qItIWXyF8Q/Z7Qn78WNRTDNTmx/XAbuCstPVHmVl3MnKtmXWMOioyNZaTG/JLoMrDtKq7WvFeVhB6uGt3AWwhDHWJjbvXRJ0fR0T/mPwV+KG7z4pzPxn7fAx4jPDN6oXAye6+JVv7k/2jwloKwTejojNlMGFsXEofM+sU3R5PGF+3Gfgt8Bkzu8PMxke9Ce8TEmBX4HEzOzM6UWQUoRifA9xoZp3MrB9hbPOBZlZiZleb2ZcaidGidj3hj88tDTyW+fzM+58kjAdMORQ4z923At9s4DUiIq2xT0+omZ0E/AD4urtXE3qH1wF/NbNPRSco/lPUQ/0wYfaQ88zsNOAr7r6NkGv/1cx+YGYnEzoYUsPo0vNfYzk5Jb33eAxQbmYHAB8DysxsePTYHsLJfqkOlvR9/BL4rJn1ju6PBf7L3T8k1DwN5mIzGxfF1bmRY9eBjJop+ufjHMK3lL0JQwwPA4YAB6SdGJnZK767gfjrB9XAeJHoZ9EDeJbQK9+afzgkRzQURArBEGCBmf2F8DWZAf837fEPgWvNbDOhd/ebAO7+tJn9K+HrxPMIJ//NBzCzTwN3Ev5Q3Atc6e5uZp8jTIW0APgd4Su9YwkXVPgcYYhJvQvTmNkJhJNyhhK+GvwT8IMonieBcUCFmR1C6EkZAHzOzBZR13vzVaCru0/JeO//beHCDh8SThYSEWm1KI98gVD8/RdQTTgBridwtrtPB3D33VF+/C9CrnuUMLYa4HJCp8T9hG/oUjNx/JqQ+y4lDFm71N3fiwrHccAhZvabxnKymU0ijCv+mpktcvfVhBMpf04oqu8HPg0cBbwDPEDoXf+yhdmijgA+YWZ/dfdHzOxQYKqZvUT41vPaqGf9fGBA9P4WEC6WM9DMjojivIAwbvxPacetlJD7+EaytgAAIABJREFUjwEGR4W3E8Y2T4y2ORd4AXiVMLzlCUJOf5xwwuQxZnaeu6e2+zvCxXq2R4X4AGCSmS0j9EQTvT7zmg0rCCdnngd0ATqY2TOEExhzNQWtNMP0s5BCZuEysP/j7hU52t+V7q7LjYuItDPRsJTO7p6XF+gysysJM7RURfe7EL5tuCl1bo4kT0NBpNA1NMwiOzsKX5m+lot9iYhI7kQnZx6Xx0V1F+B7QEXa6n7A+yqq84uGgkjBisa3fYnw1d757v77LO/y7egrShERaV86uvvdSQfRGHf/wMJl1f8zOrF9DWHYyveafqXkmoaCiIiIiIjEQENBRERERERioMJaRERERCQGKqxFRERERGLQbk5e7Nevn1dUVLT6ddXV1ZSXl8cfUAHSsaijY1FHx6JOto7F7Nmz17t7UR1k5ezmFdN7heJ6v8X0XqH9vd+mcna7KawrKiqYNav1VxKtrKzcr9e1RzoWdXQs6uhY1MnWsTCz5bFvNM8pZzevmN4rFNf7Lab3Cu3v/TaVszUUREREREQkBiqsRURERERiUPSF9eTJk5MOIW/oWNTRsaijY1FHxyJ5xfQzKKb3CsX1fovpvUJxvd92c4GYyspKb0/jd0SkeJjZbHevTDqOXFLOFpFC1VTOLvoeaxERERGROKiwFhERERGJgQprgLvvhmOPhV//OulIRESkAdte28bs42ez6J8XJR2KiEij2s081m2yaRPMmweLFycdiYiINMD3ONtmbVN3kIjkNaUogCFDwnLlymTjEBGRBpV0KwFg7/a9CUciItI4FdagwlpEJM+psBaRQqDCGuCgg8JShbWISF5SYS0ihSDnhbWZnW5mM82sIm3dRDNbZmYbzOw2M+uY9lg/M/uVmd1kZreYWWnsQaUK69WrYa+StohISr7kbBXWIlIIclpYm1l/oBtwQtq6ocB5wCTgCuBrwLfSXvZ74F53vwp4C/j/7d13nGRVmf/xz9M5h5npyZFJTBIGhowyBBVQEVkQxdUFlWHFgKsirrDiCqvIgoIYW2VVUAxgROCnBBFEcQYGBpjA5NCTekJP53x+f5zb3dU1nftWV3XV9/16nVfduvfWvadu1Tzz9Klzz/lS6BXLyYGyMp9U79sX+uFFREajRIrZlmVYhuGaHe3N7WEcUkQkdCOaWDvn9gG/i1o9A/igc+4F59xPgG8BZwOY2anAfOfc88G+fwA+bmb5oVdu+nT/uHVr6IcWERmNEilmm1lXq3WdWq1FJDGNeFcQ51x71PNnnHOtEasqgB3B8jnA1oh9K4LF8Kf+nT/fP2rIPRGRTokUs9UdREQSXSLevHgS8N1geQpwMGp7DTA5+kWVlZUsW7ass5SXlw/urEqsRWQElZeXd8YrYFy86zMMIxazlViLSLwMNGYn1AQxZjYHqHTOvRysckBD1G7ZQEv0a8vKyli1atXQT37ssf5RibWIjIAVK1awYsUKAMzsQJyrMyQjHbOVWItIvAw0ZidMi7WZpQPXAJ+LWF0BlEbsY0A+sDv0CnS0WK/XdLkiIv2JR8xWYi0iiS5hEmvgM8AdzrmmiHV/BuZHPJ8JNAPDaJruxdy5YAabN0NjY+iHFxFJMiMeszsT6xol1iKSmOIxjrV1LEas+zywGsg1s2PM7INmNsc5twrYY2aLg13PB+5xzjWHXrG8PN8dpLUVXnkl9MOLiIxGiRSzM0p878XWqtZ+9hQRiY8R7WNtZoXAvwZPrzSzb+J/Srwlatf1zrl7g+XLgS+Y2Tb8eKo3xqyCy5bBunXwwgtw0kkxO42IyGiQaDE7Y4z/L6vl4FFdtkVEEsKIJtbOuRrgO0HpcGtQenvNduBDMa6ad+KJcN99MJybIEVEkkSixezMsX4Sx9ZDarEWkcSUSH2s429ZMNTqypXxrYeIiBxFLdYikuiUWEc64QTIzoY1a+Bg9FCsIiISTx0t1i2HlFiLSGJSYh0pNxdOP90vP/VUfOsiIiLdZI4JuoIcVFcQEUlMSqyjnXuuf3ziifjWQ0REuskYq64gIpLYlFhH60isH30UnItvXUREpFNHi7W6gohIolJiHe3kk2HSJNi+XaODiIgkkMyyILHe34JTw4eIJCAl1tHS0uCyy/zyL38Z37qIiEinjMIM0gvSaW9o1yQxIpKQlFj35PLL/eP990OLfnIUEUkUWVOyAGiqaOpnTxGRkafEuiennQYLF8LevfDb38a7NiIiEsiekg1Ac0Uos6SLiIRKiXVPzODaa/3yt78d37qIiEin7Kk+sVaLtYgkIiXWvXn/+6GgAP7yF83EKCKSIDparJVYi0giUmLdm6Ii+OhH/fLNN8e3LiIiAkQk1juVWItI4lFi3ZfPfMa3Wj/6KDz7bLxrIyKS8nKOyQGgYXNDnGsiInI0JdZ9GTcOPvUpv/zRj0KrhncSEYmnvHl5ADS8rsRaRBKPEuv+3HADzJoFa9bAXXfFuzYiIikte0Y2lmk07Wqira4t3tUREelGiXV/8vLgW9/yyzfdBC+/HN/6iIiksLSMNHJn5wJQv7E+zrUREelOifVAXHABrFgBTU3wnvdAXV28ayQikrJy5/vEumGDuoOISGJRYj1QX/+6nzRm/Xq45hpwLt41EhFJSfmL8wGofak2zjUREelOifVA5eXBL34B+fnw05/Cl78c7xqJiKSkwhMLAahZVRPnmoiIdKfEejAWL4YHHvAzM950E9x/f7xrJCKScgqXBYn1CzU4/XooIglkxBNrM3uzmT1vZjMj1uWa2T1m9mUz+46ZFQ1kW1y84x1wxx1++d/+zbdii4gkqUSM2dlTs8ksy6T1cCuNWxvDPryIyJCNaGJtZhOAAuDkqE3fBp51zn0e+B3wvQFui49PfcrPxtjeDu97H9x3X7xrJCISukSN2WZG0ak+Xz/yzJGwDy8iMmQjmlg75/bhA20nM5sMXAH8IVj1OHCRmc3sa9uIVLgvN9/su4O0tcEHPgBf+YpuaBSRpJLIMbvknBIADj95OOxDi4gM2Yh3BXHOtUetWg4ccM7VB9tbga3AWf1siy8zuOUWP2mMGXz+8/CRj0Bzc7xrJiISmkSN2aXnlAJQ9WSV+lmLSMJIhJsXpwAHo9bVAJP72dZNZWUly5Yt6yzl5eUxqexRrrsOfvlLyM6G730Pzj4bdu8emXOLyKhVXl7eGa+AcfGuzyAkRMzOX5xP5rhMmnY1Ub9eE8WISGwNNGZnjFyVeuWA6FH+s4GWYLmvbZ3KyspYtWpV+LUbiEsvhenT4V/+BZ57Dk48EX71KzjzzPjUR0QS3ooVK1ixYgUAZnYgztUZjISI2ZZmjHnbGPb9eB8Hfn2A/Bvzh3wsEZH+DDRmJ0KLdQVQGrWuENjdz7bEcvLJ8MILcNZZsHevf7zpJnUNEZFkkzAxu+xfygCofKgyFocXERm0REisnwImm1kugJllAtOBp/vZlnjGj4fHH4cbbvA3Mv7P/8Cpp8Jrr8W7ZiIiYUmYmF365lLSC9OpXV1Lw2ZNby4i8RePcaytYxHAObcXeBg4N1h/HvCQc66ir20jWOXByciA226Dv/wFZs6E1ath6VK48UaoVz9AERldEjlmp+ekM+5i39Vxz//ticUpREQGZaTHsS4E/j14eqWZlQXLHwHeaWY3AhcCKyJe1te2xPWmN8HLL8PVV0NLi58CfeFC+P3v410zEZEBGQ0xe9LVkwDY+8O9tLdED2AiIjKyLFmGKVq2bJmL282L/XnuObj2Wp9oA5x7rm/V9neWikiKM7MXnHMpFRDCitnOOVYuWkn9unoWPbios9+1iEis9BWzE6GPdfI7/XRYtQruvhuKi+GJJ+Ckk+Dyy2HjxnjXTkRk1DIzJl/rR/Pb8dUdGtNaROJKifVIyciAT3wCtmyB66/3417/8pe+e8hVV8Hrr8e7hiIio9KkD04ic3wmNStrOPwnzcQoIvGjxHqkjRkDt9/uW6qvugra2+FHP4Jjj/Ut2B3dRUREZEDS89KZ9ulpAGz74ja1WotI3Cixjpdp0+Dee2HDBn+DY0aGb8E+/ni44AJ49FGfdIuISL8mf2QymeMzqf5HNft/vj/e1RGRFKXEOt7mzIHyct9F5LrrIDcXHnsMLrzQt2J/4xtQXR3vWoqIJLSMwgyO+fIxAGz57Bba6triXCMRSUVKrBPF1Klw112wcyd89at+ivSNG32yPWUKXHMNPP+8n3hGRESOMvHKiRScUEDTria23Lgl3tURkRSkxDrRjB0Ln/0sbN4MDz0Ey5dDba1v1T71VFiyBL72NdivnzpFRCJZujG/fD6kQ8XdFVQ9XRXvKolIilFinagyMuCSS+Cpp+DVV+HTn4ayMj89+qc/7VuxL7kEfvMbaGyMd21FRBJC4YmFzLhxBgDrr1xPy6GWONdIRFKJEuvRYNEiuOMO2LXLJ9Jvf7u/sfE3v/HJ9cSJ8KEPwZNPQpv6FYpIaptx4wwKlxXSuK2Rde9bh2tTFzoRGRlKrEeTrCy4+GL4wx98X+zbb/ejiBw54kcYOfdcP9rIpz7lJ6RRf2wRSUFpWWksenARmeMyOfTYIbZ+YWu8qyQiKSKUxNrMLgzjODIIkyf7iWZWr/bdQ266CY45Bvbsga9/3c/sOG8e3HAD/OMfGrpPRDqlQszOmZHDwl8shDTY8eUdVHy3It5VEpEUMKjE2szON7NnzWyTmW0JyjbgodhUTwZk4UK45RbYtMkn0R//OIwf75/ffjucdppvyb72Wnj8cWhRn0ORVJDqMbv0nFLmfXceABuv3cj+X+mmbxGJrYxB7n8f8BVgNdDRBGrAFWFWSobIDE45xZevfQ2efdb3w/7Nb3zXke98x5eSEnjLW/xENOef7/toi0gySvmYPfnqybRUtrD1xq2se986LMMoe1dZvKslIknKBjP1q5k975w7pYf1E51ze0Ot2SAtW7bMrVq1Kp5VSFzOwYsvdiXZa9d23750qU+yL7zQJ+UZg/17S0SGw8xecM4ti8FxFbMB5xybr9/Mrjt3QTos+PECJrxvwoicW0SST18xe7B9rD9iZpf2sP6iwVdLRowZnHgi3Hqr74+9aRPcc49PpHNyfD/tL38ZzjzTD+l32WXw3e/C66/rBkiR0U0xGzAzZv/vbKbfOB3aYN3717Hzrp0MpmFJRGQgBtti/XHgFqAwcjXgnHPpIddtUNRiPUQNDfD00/Doo75s3Nh9+9SpcM45XWXatPjUUySJxbDFWjE7yvbbtrP1P/0oIZNWTGLuN+eSlqkBskRk4PqK2YNNrA8D1wCrgI4Bkw14v3PuluFWdDiUWIdk82Z/g+OTT/py4ED37XPm+AT7TW/yLdzTp/sWcREZshgm1orZPdj3wD7WX7Ue1+QoObuEBT9bQPbE7LjURURGnzAT658BVznnmqLWp1R/vZTR3u67jjzxhE+yn34aqqu77zNlik+wO8qSJZAe14YwkVEnhom1YnYvqp+v5tWLX6V5bzOZ4zNZcP8Cxrx5TNzqIyKjR5iJ9UXACcCTkauBS5xz1w2rlsMU7yCdElpb/U2QTz0Ff/ubH3Xk8OHu+xQWwumnwxlnwKmn+vG0S0riU1+RUSKGibVidh+adjex7l/XUfVUFRhMu34aM/97Juk5ahwQkd6FmVivBOYCVUDHCw2Y6JzLGW5FhyMRgnTKaW+H9et9gt1RtvYww9n8+X60kZNP9uW44/wskiICxDSxVszuh2tzbP/ydrZ9cRu0Q96xecy/dz7FpxXHu2oikqDCTKzfDTzqnKuJWn+hc+6RYVZyAfAxYBP+P4Jy59xLZpYL3A7UAKXADc656ujXJ0qQTnkVFb41+29/g3/+04840tTUfZ+sLD/E3ymn+BbtE0/0s0SqC4mkqBgm1orZA3TkuSOs/+B6GjY0gMGUT0xh1n/PIqNYw4+KSHd9xezBRoyd0QE6EMYMI/cB73TOVZjZdOD/AQuAbwOPOed+YWbnA98D3hvC+SQWpkyBd7/bF4DmZlizxifZzz/vH9ev98vPP9/1urw835K9dKkvJ5wAixZBtm4oEhkGxewBKj69mGUvLWP7l7az4/YdVNxdwf4H9nPMl49h4pUTsXTdpC0i/Rtsi/V9zrn3R627DPihc65oWBUxqwNOdM6tN7My4GVgGbAVKHXO1ZtZBnAEWOSc2xb5+kRr/ZA+HDkCq1b5xHrlSt+qvX370ftlZvrkuiPZPv54f3Ok+mxLkolhi7Vi9hDUrK5h48c3Uv0339BecEIBs++cTeny0jjXTEQSQZhdQf4I/MM5d4uZ5QHfAC4DVjnnzh1mJX8MzAYuAD4M7ACygf91zk2J2O/VYN2PI1+fyEFaBuDgQXjpJZ9kv/iif9ywoecJaqZO9Qn24sX+cckSOPZYP9mNyCgUw8RaMXuInHPs//l+tnx2C027fHe20vNKmXnLTIpPVf9rkVQWWlcQ59zbzOxUM7sHeCuwFzgO3yIxXNcCvwNWAnc65x4ys+uBg1H71QCTQzifJJKxY+Hcc33pUFvru5F0JNtr1vjh/3bt8uXRR7v2TU+HuXO7J9uLF8Mxx6jvtqQsxeyhMzMmvHcC4y4ax847d7Lzzp0cfvwwhx8/zJi3jWHmzTMpOmlYjf4ikoT6TKzN7E29bNoObATuAGYAFwP/Mcy65AIPAFOAu8xsK/4u9oao/bKBlugXV1ZWsmxZ1x8PK1asYMWKFcOsksRVQYEfuu/007vWtbXBli3wyivw6qv+8ZVX/IyR69f78uCDXfvn5PgbI4891pcFC/zjvHm+X7dInJSXl1NeXt7xdFwYx1TMDl96fjozvzCTKR+bws47d7Lr7l0c+uMhDv3xEMVnFTP9+umMuWAMlqY+2CLJbKAxu8+uIGb2D/zNKFVAe2+7AZOcc8O6y8zMHgPe45yrMrP/Aj4BfBK42Tk3L2K/jcG6n0W+PtF/VpQYa2yEdeu6J9uvvOJHKenNjBldiXZk0l1WptkkZUSF1RVEMTv2mvc3s/POnez+7m7aqv1klnkL85j26WmMf+940nP1C5lIshtyH2szuxDfF29/Pye41Dn3YF/79PP6ccArzrlJwXPDD+H0b8BjQJlzrsHMMoFa4BjnXLeMabQGaYmx6mrfV3v9ep94d7Rqb9zoJ7zpSWmpH3t77lw/hfvcuV3LunFSYiDExFoxe4S0Hmll9/d3s+uuXTRXNAOQUZrBxCsnMvmayeTN1y9iIskqtJsXY8XM0oDNwJnB0E0GPA+cB5QD9zvnHjazC4D3O+euiD7GaA/SMsJaWnyXko5EuyPpXrfu6GnbI40b1z3R7lieOxeK1N9ShiZWNy/GimJ2l/bmdvb/fD8V36ygZmXXyIYl55QwecVkxl40Vq3YIkkm4RNrADM7Dn8zzAvABOCvzrmnzawUP9nANvzYq//pnKuNfn2yBGmJM+dg717foh1dNm2ChujuoxHKyron3bNn+5snjznGJ+TqXiK9GG2JNShm96TmhRp2f3c3+362j/Z63xMnvSidssvKmPiBiRSfWay+2CJJYFQk1sOVjEFaEkx7O+zZ03PCvWmT7+fdm4ICn2DPmtWVbHeUmTM1VGCKG42J9XAlc8xuqWph3/372PfjfdSs6mrFzp6RzYR/ncD4d48nf0k+pj+2RUYlJdYisdbe7m+U7Ei0N26ErVt9d5PNm/vuXgJ+xsqeku5jjoGJE9XaneSUWCevunV17LtvH/vu30fTzqbO9blzcyn7lzLKLi2j4IQCJdkio0hME2szWwp8GlgNPBQ9u9ZISZUgLaOQc3D4cFeiHVm2bvWzTvZ2IyVAbq5v1Y5s4Y4spaVKvEe5kUysFbPjw7U7qp6uYv8D+znw2wO0VHaNQJgzM4dxl4xj3MXjKDqtiLSMtDjWVET6E+vE+nvAO51zE83sCuA3zrk+OqLGRqoFaUkira1+wpvIZDsy+T5woO/XFxb6BHvGjKOT7pkzYcwYJd4JboQTa8XsOGtvbefIM0eofKiSA78+QPOe5s5tGSUZlL61lLFvG8uY88eQVZYVx5qKSE9inVhnAuODO8PTgALnXD+/e4cvlYO0JLmamu7J9vbtsG2bL1u3+u19KSjoSrJ7Sr7HjlXiHWcjnFgrZicQ1+6o/ns1lb+u5ODDB2l4PeJvHIPCkwo7k+zCEwuxdP1bFYm3UBPrYPikMc65n5rZyUCtc25tCPUcFgVpSUnOQVVVV6LdU+mvf3d+fvdEOzr51ogmMRfLxFoxe3Sp31TPoUcOcfCRg1Q9VYVr7vo/Or04nZLlJZSeW0rpuaXkLchT32yROAgtsTaz/wGuAx5zzl0arLsNWBM9q9ZIU5AW6UVfiffWrf0n3nl5fSfemqly2GKVWCtmj25tdW0cfvIwB/94kMOPH6Zxc/eRh7ImZVFyTgml55RSsryEnFk5SrRFRkBfMTtjkMc6HZiKn12rw6PAfUBcg7SI9KKkBI4/3peedCTekV1MIktVFaxd60tPOm6u7C35Hj9eiXf8KGaPYun56Yx7xzjGvWMcAI3bGzn8xOHO0rynmf0/3c/+n/qJNrMmZ1F8ZnFnKXhDgbqOiIywwSbWf3bOVZlZZDP3pWFWSERG2EAS7+3bu0YwiU7CDx/2M1auW9fz63Nzj27ljnw+YYIS79hRzE4iOTNymPTBSUz64CScc9SvrfdJ9pOHOfLsEZp3N1P5y0oqf1kJQHphOkWnFflE+43FFJ1URHq+ZoEUiaXBJtZbzOzzwBQzew9wGfAu4OOh10xEEkNJiS/HHdfz9iNHem/t3r4dDh3qmjq+Jzk5Ryfes2Z1zV5ZUhL+e0oditlJyszIX5RP/qJ8pn5iKq7dUb+hniPPHuksjVsaOfynwxz+02H/onTIX5xP0clFFJ1SROHJheQvzFertkiIhnLz4qnAlcB0YA/wgHPu8fCrNjjqryeSoKqre0+8t23ziXdfxo7tSrKjH5Okm0mMb15UzE5RTbubuiXatWtqoa37Pmn5aRSeWNiZaBedXET2tGz11RbpQ6yH20sHljvnnhjWgYZJQVpklKqpOXoIwS1b/AyWmzdDQx9DLBcU9J50T50KaaNjoo0RHm5PMTtFtdW1UbO6hpp/1lD9z2pqnq+hcVvjUftlTcyi4IQCCpYWULi0kIKlBboxUiTCkG9eNLPfAwX9HH8ccACIa5AWkVGqsBAWL/YlmnOwZ49PsDsS7cjHqip46SVfomVnd3Up6Ui2586FBQtg2rSkaOmOppgtfUnPT6fkzBJKzuzqXtW8v5malT7Rrn6+mpp/1tC8t5lDjxzi0CNdvyZllGRQcLxPtjtK3rF5miVSJEp/fazXBPtsANqBq/HT4EY2M8wCxsSkdiKS2sxg8mRf3vjGo7cfOtRzwr1pE+zb13vf7vx8n2AvWAALF/qyYIGfMj59VN/cpZgtg5I1PouxbxvL2LeNBcA5R8PmBmpX11K7upaaF2uoXV1Ly/4Wqv5SRdVfqjpfm5aTRv6SfJ9oH19A/pJ88hfnk1mSGa+3IxJ3/SXW/4vvLlIFYGaznXNfiN7JzH4Xi8qJiPRpzBg4+WRfotXW+kQ7MunesMGPXrJ/P6xa5Uuk7GyYP9/fqLl0KZxwgh8tpbh4ZN7P8Clmy7CYGXlz8sibk8f4y8YDPtlu3tPsE+3VNZ1Jd+PWRmpW1lCzsvvsr9nTsslfnO8T7SX5FCwJWrez1botya/PxNo5dyRq1QQzS3fOdd7+YGbLgR7+VxMRiaOCAp8g9zSaycGDXUMErl3b9bhzJ6xZ48t993XtP3u2T7JPPx3e9CZ/zARs2VbMllgwM7InZ5M9ObuzZRugpaqF2pd8kl23po7aV2qpf62epp1NNO1s4tCjETcmp0PevLzOVu2OhDtnVg6WlnzdsiR1DXa4vQeBl83sWaAGWAC8FTiqRUREJGGNHQtnnulLpJoan2CvXt1V1qzpavn+1a/8foWFcMYZ8Pa3w8UXw5QpI/8eBkYxW2ImsyST0uWllC4v7Vzn2nxXkrpX6qh71Sfbda/U0bCpgfp19dSvq+8cZxsgLS+NvPl55C3MI39BPnkL8shbkEfunFzSMtXCLaPPUIbbmwxcC7wBqAZ+5ZyL+8+KusNcRGKipcX30165Ep55xpfNm7vvc9ZZcP31cOGFQ7opMsbD7SlmS9y1NbRRv66eule6ku26V+to3t3c4/6WYeTOzSVvQR75C7sS7rz5eaTnJd6vRZJaYj3cXhqw1Dn3wrAONEwK0iIyYnbvhieegN/8Bh57rGtIwA99CH7wg0EfboSH21PMloTRcrilsyW7bl0d9Wv9ck/DAAJgkDMzpzPR7ky65+eROUY3TcrICC2xNrN7e1hdBrQ55y4eYv1CoSAtInFRUwNf+Yov4G+QnDdvUIeIVWKtmC2jVVt9G/Ub6qlfGyTc6/xyw6YGXGvPeUvmuExy5+f6riXz88id55dzZ+eSlqVuJRKeIY9j3YO5+LFPI7/VU4Fnhlg3EZHR7Xe/g+98xy+PHQsTJsS3Pt0pZsuolJ6XTuHSQgqXFnZb397S3tlfu25tXWdrd/3r9bQcaKHlQAvVf6uOOhjkzsrtTLTz5ud1JuBZE7M08Y2EarCJ9YedcxsiV5jZsfjpckMRTL97JrAJeBaoA27H33hTCtzgnKvu/QgiIjHW1gbPPw/f/z786Ed+3Xnnwbe/nWhD8ylmS1JJy0wjf0E++QvyKbukrHO9c46miiYaNjRQ/3o99Rvq/fIG362kYVMDDZsauk16A5BemN494e5YnpdHer76csvgDSqxjg7QgSbgw8DnhlsZM7sKmO2cuyli3f8BjznnfmFm5wPfA9473HOJiAzY/v1+dsfVq/1NjE895SenAT/s3q23wg03JNxsjorZkirMjJypOeRMzaH03NJu29oa22jc3Oi7lgSl4XWfdLceaqX2hVpqX6g96pjZU7N9y/a8iFbueXmbrsz8AAAgAElEQVTkzMjB0hPr37okjkEl1ma2le4/KaYDE4D7en7FoI59FvA+4M0R6yYDVwAfDVY9DjxkZjOdc9uGe04REQDa2+HAAT+O9datsHGjn1Rm40Zf9u49+jUzZsBll8HVVw+6T/VIUcwWgfScdPIX5ZO/KP+obc0Hmnts5W7Y1EDTriaadjVR9URVt9dYtpE7p6tlu7Nf97w8MsfqBspUN9iuIH8GfkZXoG4H9jnnXg+hLl8DngPuMbM5wC3ADOCAc64ewDnXGvxHcRawLYRzikgya2vzLcv79/uyb58vFRWwa1dXqaiA5p6H/QL8ZDPHH+/L0qV+eL1jjkm4FuoeKGaL9CFrXBZZ47IoPqN7F6721naatjf12MrdvLuZ+tfqqX+t/qjjZYzN6DHhzp2Tq5knU8RgE+v7nXN/jVxhZvPNbI5zbtNQK2Fm84GlwFXOuTVmdj3wGHArcDBq9xpg8lDPJSKjXEsL7NnTlSTv29eVNEc+7t/vW6Hb2wd23DFjYOpU3xI9d64vc+b4x2nTIG1U/qeomC0yBGkZaeTOziV3di5jLxzbbVtrTSsNGxuOSrgbXm+g9WAr1c9VU/1c1G0FacEwgdEJ9/xcsqdk6wbKJDLYxPp84K9R6/YBT+GD7FAtAg4559YEz78J3IxvZWmI2jcbaIk+QGVlJcuWdY18smLFClasWDGMKonIiHMODh/2rcg7dnSV7du7lnfvHniyDD5hHj++e5kyxSfR06b5xylTIC8vdu+rB+Xl5ZSXl3c8HRej0yhmi4QsozCDwhMKKTyh+4glzjma9zQflXDXb6incWsjjVt84bHux0vLS+sx4c6bl0dG0WDTNImVgcbsAY1jbWYfw7dEFPayy5+dc+cPtpIRx383UO6cK4lYtxr4Jb5FZF7E+o3Azc65n0UeQ2OiiiSwpiaorPR9lffu9S3Oe/Z0LUc+NjX1fSwzmDQJJk70Q9uNH9/12LHc8XzcOMhM/D6PYY9jrZgtkljam9tp2NxwVMLd8HoDLZVH/d3ZKWti1tEJ9/w8cmblkJYxKn9FSwrDHsfaOfdNM/sr8G/AN6I2Nzvn9gyzjmuAYjOb4JzbF6xrBXYDk80s1znXYGaZwHTg6WGeT0SGqr0dqqp8N4vKSv8YudzTupqagR+/sNC3Ik+f7rtlTJ/evUyZAllZsXt/SUAxWySxpGV1DRMYreVQC/WvRyXcGxqo31hP895mmvc2c+TpI91eY1nmRytZkEfewjzyF+STt9An3+rLHV8D/o0h6Ef3OefcUX9amdl059yOoVbCObfezB4FLgW+ZWYlQBbwU+AC4FzgYeA84CHnXMVQzyUiEdrb4cgROHiwKxk+cKDv5wcPDq4rBkBGhm89njixq7U58rFjeeJEyD/6Px4ZPMVskdEhc0wmxacWU3xq9xsoXbujcUdjj6OWNO1sou7VOuperet+sDTInZ3bPdlekEfesXlkFKhbyUjo8yqb2WeAR5xza4NVnzWz6BHTDTgbWD7MunwAuNvMsvAtHO8N7ij/CHC7mR0HTATUCU+kN/X1vpU48ua9vpLlgwf9yBmDVVQEZWW+jBvX9djbcnHxaBhBY9RTzBZJHpZm5M7MJXdmLmPeOqbbttbaVurX+2neI2ehbNjcQMNGXw7+rvt9xNnTs8lfmN/Vyr0on/zF+WQUKuEOU39XcxZQEvH8BCAf2EvX8E2GD57D4pw7gB8TNXr9YeDq4R5fZFRqafGJcmSy3FF6WldX1/8xoxUVdSXC48b5abl7ez52rC+joN9yilLMFkkBGQUZFC0romhZUbf1bY1tfsSStV3Jdt3aOhpeb6BpRxNNO5o49Fj32SdzZuaQvySf/DfkU7CkgPwl+eTOy1Uf7iHqM7F2zn00atX/AK8657oN+GpmC8OumEhSamjo3mocWXpKmA8fHtzxs7K6buIrK/OPfSXKY8aov3ISUcwWSW3pOekULCmgYElBt/Xtre00bm6kbl1dZ9Jd96pPvBu3NdK4rZGDf+hq4bYsI29BXmei3VE0NGD/Btv+v8Q592LkimA809nA2p5fIpKkmpqO7l7RU4ncp/7oCQX6lJbWlSBHJsuRJXJdYaG6XEgkxWwRIS0jzY8sMj8PLu5a397STsPGBmrX1FL3Sl1nadzWSN3LddS93P1X0IzSDPIX+yS74PgCCk8oJH9xvm6YjDCgxNrM3oj/+fDNwSxakcbg7zr/Q8h1Exk5ra09J8l9Jc6DGemiQ1ZW9z7Ika3IPSXMY8aM1olJJI4Us0VkINIy08hfmE/+wnx4T9f61upW6l7rSrRrX6mlbk0drYdbOfLMEY480zVKiWUY+YvzKTjBJ9oFJxRQ8IYC0vOjb+9IDQNtsa4CfoLvv3dG1LZm4AdhVkpkWCKHg4se/q23UlU1+POkpx+dIPdX8vPVoiwjQTFbRIYsoyiD4tOKKT6ta6QS5xzNu5s7k+zal2qpXV1L/YZ6v/xSLXvv3et3ToO8Y/M6E+3CEwopXFaYEsn2QMexfsXMlgOXOecUkGXkOOdbhvtLjKNbmQc7HJzZ0X2Re7qRL7JopAtJUIrZIhI2MyN7SjbZU7IZe37XNO+tta3UvVxHzYs11L5YS82LNdS95vty16+tZ9/9wVD36VCwpICi04ooOrWIotOKyJ2Tm3R9tgczjvURemnlMLMznXPPhlYrSV593bwXWSJbmVt6n5WqV8XF3ZPg3rpfdKwvKfEt0CJJQjFbREZCRkEGxWcUU3xGV+t2W0Mbda/WdSbaNStrqF1T29myvfs7u/1rx2ZQdGqRbx0/s5iiU4tGfX/t/sax/j1Q0Nc++H58i4DxYVVKRomWlq4+yAPpbjGUm/fAd58YTHcLDQcnKUoxW0QSQXpuOkUnFVF0UtdwgG11bdSsqqH6H9Uc+fsRqv9eTcv+Fg798RCH/uiHAEzLSaPo9CJKzi6h9OxSCk8qJC1rdCXa/bVYb8D31augawzUaAa8PcxKSRw1NsK+fb7s3du1HF327vUz9g1WVlb/LcnRSXJubvjvUyQ5KWaLSEJKz0+n5KwSSs7yQ+0752jc1kj1P6qpfq6aqqerqHuljqonq6h6soptbCMtL43iM4sZe+FYxr59LLmzEz8fMOd6i71gZuOAGudcU58HMRvrnDvY1z6xtmzZMrdq1ap4ViFxtbb6RHj3btizp+ckuWO5unrgx01P77sPck+loED9kkWimNkLzrllIRxHMVtERq3mymaqnq6i6ilf6td1/5U7d34uY9/uk+ySN5Zg6fHJJ/qK2f1NEHNggOd4F7rLPD6qq6Giou+yb9/Ab+bLzPRDvU2YABMn+sfo0rG+tFRDwYkkEMVsERnNssqyGH/peMZf6nuqNe9r5tCfDnHwjwc59NghGjY0sGvDLnbduYusiVmUXV7GhCsmUHhSYcLcBDmsCeLNbDrwUeBjKEiHr63NJ8bbtnUv27d3Jc21tf0fx8wnwlOmwKRJvSfKHclygnw5RSRcitkiMppkTchi4vsnMvH9E2lvaaf6uWoOPnyQyl9X0rilkYq7K6i4u4Kc2TlMvnoyEz80kaxx8Z1NeNCJtZmlAxcBVwNvBvbh+/TJUBw8CBs2wKZNRyfQO3f6bhx9yc31CXNfZdIk3cwnkqIUs0UkGaRlpnX20T7m9mOoWVnD/gf2s//n+2nc3MiWz21h681bmfDeCUz99FQKFvd3H3dsDDixNrNZwIeBq4CxwHbgHcD/A06LSe2SRXMzbNzoE+jXX/ePHcsH++nmOGkSzJwJM2b4x47lqVN90lxSohZmETmKYraIJCszo+jkIopOLmL2HbM59KdDVHyrgkOPHGLvj/ay98d7GX/5eGZ+caafxn0E9ZtYm9ll+JaOc4BNwNeAHwOfd849Fuz2XMxqONocOAAvv+zLSy/5x7Vre295LiiAefN8mTWrK3meOROmT4ecnBGsvIiMdorZIpJKLN0Ye8FYxl4wlobNDey6axe7y3ez/+f72f+r/Uy/fjozvjCD9NyRmatiIC3WDfgpcJ8DPuycex3AzHofTiRVtLXBmjXwt7/58txzsGPH0fuZwezZMH9+V5k3zz9OmqQWZxEJk2K2iKSk3Nm5zL1nLtOun8a2L21j77172XHbDip/XcmihxaNSPeQfhNr59zDwMPBTS9Xm9lM4HdAZ6ddMzPX17h9yWTrVnjkEV+eecZPtx0pLw/e8AY47jg4/nhflizxk5yIiMSYYraIpLqc6Tkc+4NjmfShSWz48Abq19bz4qkvsvDnCxn39nExPfdgpjTfAfyXmWUAFwNTzewW4IfAW4Dy2FQxAWzdCj/5CfziF7BuXfdts2bBGWd0lYULNTW2iMRdSsdsERGg+LRiTlx5Ihuu3sD+n+3ntUteY/HvFzP2/LExO+egRwVxzrUCDwIPmtkc4D/w/fmSK0i3t8Mf/gBf/zo8/XTX+qIieMtb4MIL/eOUKfGro4hIP1ImZouI9CA9L50F9y8gsyyTirsrWHvZWk5cfSJ5c2JzU+OwxrF2zm0CrjOzF0KqT2L44x/hhhvgtdf889xcuOQS+MAH4OyzNXSdiIxKSRuzRUT6YGbM+focmiuaqXywkvUfWM/Svy2NyaQyw0qsOzjnfhLGceLu0CG45hp48EH/fMoUuP56uOoq31ItIpIEkiZmi4gMkJkx7/vzqHqmiuq/V3Pgtwcoe1dZ6OfRfNQdNm6EU07xSXV+Ptx5J2zZAtddp6RaREREZJTLLMlkxo0zANh1166YnCOhEmszyzKzl81sefA818zuMbMvm9l3zCw2GW5FBZx7rp/98Pjj4dVX4VOfgqz4TospIpLo4ha3RUSGYOK/TSQtN40jfz1C447G0I+fUIk1cD0wM+L5t4FnnXOfxw8X9b3Qz9jeDldc4acPP+00ePZZPzmLiIgMxMjHbRGRIcooymDMW8cAcPjJw6EfP2ESazM7HdgDHA6eTwauAP4Q7PI4cFEwJmt4fvpT+OtfYfx4+P3vNd60iMgAxS1ui4gMQ/EbiwE48uyR0I+dEIm1meUDlzrn7o1YvRw44Jyrh84ho7YCZ4V2Yufgq1/1y7fdBuNiO2i4iEiyiFvcFhEZpoKlfgbG+vX1oR87lFFBQnADcFvUuinAwah1NcDkng5QWVnJsmXLOp+vWLGCFStW9H3WNWv8kHrjx8P73jfYOouIDFl5eTnl5Z1DSY/Gv+qHFbeHFLNFREKQOzsXgMbNA+9jPdCYHffE2szeCqxyzu2P2uSAhqh12UBLT8cpKytj1apVgzv5I4/4x4su0o2KIjKiIhNJMzsQ5+oMShhxe0gxW0QkBNlTsyEdmvc2097cTlpW/x04BhqzE6EryGeAe83sQFDRafgbXvKB0qh9C4HdoZ25I6ifpV8pRUQGIX5xW0RkmCzNyCjxbcutVa2hHjvuLdb4G12yI57/HfgU8DxwvZnlOucazCwTmA483cMxhmbNGv/4hjeEdkgRkRQQv7gtIhKCzNJMWg+20nq4lazx4fVaiHti7ZyrjHxuZm1ApXNuh5k9DJwLPAycBzzknKsI6cSwY4dfnjMnlEOKiKSCuMVtEZGQZJT6FLjlUI89jId+3FCPFr6PALeb2XHARCC8O1sOH4bmZj+rYl5eaIcVEUlxsYvbIiIhScvzvaHbG9tDPW7CJdbOuZkRy4eBq2Nyor17/ePEiTE5vIhIqhixuC0iEpKOGxZdswv3uKEebTSprfWPxcXxrYeIiIiIjCjLMgDam8NtsU7dxLquzj+qG4iIiIhISknLVot1uOqD2XaUWIuIiIiklI6uIGqxDosSaxEREZGUZJm+K4harMPS3Owfs7P73k9EREREkkuQAbt2JdbhaGvzj2mpewlEREREUpGZb7Em3LxaiTXp6fGth4iIiIiMrI4MWIl1SNqDzupqsRYRERFJLR0N1uoKEhK1WIuIiIikJHUFCVtHi7USaxEREZHUEuTVSqzDopsXRURERFJTR4O1U1eQcKgriIiIiEhKsjR1BQmXbl4UERERSU0dXUHCnXgxhRPrjqZ/JdYiIiIiqUVdQUREREREQqCbF0VEREREhk/D7YUt5KZ/ERERERklNPNijHT8xSIiIiIiqUEzL4qIiIiIDJ+6goiIiIiIhCFGHRYSJrE2s0vMbKuZHTSzu80sI1g/zsy+Z2ZfNbOvmVlmKCdUH2sRkSEb8ZgtIjIKJERibWbTgYuBS4H/AD4MfDLY/CvgXufcDcDrwJdCPnmohxMRSXZxjdkiIgksI94VCMwAPuicawVeMLPFwNlm9iww3zn3fLDfH4ANZnarc64uXpUVEUlxitkiIj1IiBZr59wzQYDuUAHsAM4BtkbsVxEsLhvB6omISATFbBFJGily8+JJwHeBKcDBqG01wOToF1RWVrJs2bLOUl5e3vcZ1MdaROKovLy8M14B4+Jdn2GKfcwWEQnTIHsCDzRmJ0pXkE5mNgeodM69bGYOaIjaJRtoiX5dWVkZq1atGsoJh1JNEZFhWbFiBStWrADAzA7EuTpDNuIxW0QkDgYasxOqxdrM0oFrgM8FqyqA0ojtBuQDu0e+diIiEkkxW0Sku4RKrIHPAHc455qC538G5kdsnwk0A2rmEBGJP8VsERnVXMhdgxMmsTazzwOrgVwzO8bMPghUAXuCO84Bzgfucc41D/uE6mMtIjJkIx6zRUTCFKOewAnRx9rMbgJuiVq93jl3r5ldDnzBzLYBBcCNIZ881MOJiCS7uMZsEZEElhCJtXPuVuDWXrZtBz40sjUSEZHeKGaLiPQsYbqCiIiIiIiMqBQZxzr21MdaREREREKUuol1B/WxFhEREUkpFqP8T4m1iIiIiEgIlFiLiIiIiIQgdRNr9bEWERERSW26eTFk6mMtIiIiklpilP4psRYRERERCYESaxERERGREKRuYl1QAJMmQWFhvGsiIiIiIiMovSidrElZpOenh3rchJjSPC4++lFfRERERCSlzLxpJjNvmhn6cVO3xVpEREREJERKrEVEREREQqDEWkREREQkBEqsRURERERCkPKJdXl5ebyrkDB0LbroWnTRteiiaxF/qfQZpNJ7hdR6v6n0XiG13q8S6xT6sPuja9FF16KLrkUXXYv4S6XPIJXeK6TW+02l9wqp9X5TPrEWEREREQmDOefiXYdQmFklsH0ILx0HHAi5OqOVrkUXXYsuuhZdYnUtZjjnymJw3ISlmD0gqfReIbXebyq9V0i+99trzE6axFpEREREJJ7UFUREREREJARKrEVEREREQqDEWkREREQkBCmbWJtZrpndY2ZfNrPvmFlRvOsUK2Z2iZltNbODZna3mWUE68eZ2ffM7Ktm9jUzy4x4Ta/bkoGZZZnZy2a2PHje6/chFb4rZnaqmX3GzC4OPvuUux5mtsDMvmVm/2Fm3zaz44P1KXctElGqXWszu9XMXFBejnd9wmZmbzaz581sZsS6pP2Me3q/wfqk+pyHkm8kHedcShbg/4DLg+XzgQfiXacYvc/pwE+AE4EPAHXAZ4JtTwGnBMv/Dnwl4nW9bkuGAtwIHAGW9/d9SPbvCnAVcGvUupS7HsAqYEqwPB1Yl6rXIhFLKl1rIB/4LnBeUObGu04hv78JwLsAB8xM9s+4j/ebVJ/zUPONZCspOSqImU0GtgKlzrn64C+qI8Ai59y2uFYuZGb2RuDvzrnW4PntwCLgFuDXzrnJwfopwAZ8AFjS2zbnXN3Iv4twmdnpwLHAF4Argdfp5fsANPe2LRm+K2Z2FvBfwJtdEAz6+vdBEl8PM6sDTnTOrTezMuBlYBkpeC0STSrFbAAz+2Sw+F3nXGNcKxMjZpYGtAGznHPbkv0zjn6/wbqk+pyHkm8kQ04RLVW7giwHDjjn6gGCL8FW4Kx4VioWnHPPdHzJAxXADuAc/Hvu2K8iWFzWz7ZRzczygUudc/dGrF5O79+HvrYlg68B64B7zOwxMzuD1L0eDwI/MLNC4F+Bj5O61yLRLCdFrrWZGfAe4E5gn5m9N85VignnXHvUquUk8Wcc/X6T8XMeYr6RdFI1sZ4CHIxaVwNMjkNdRtpJ+J+e+roGyXx9bgBui1qXktfCzOYDS4HvO+c+BjwBPAZMIwWvB3At0AisBGqdcw+Rot+NBJQy19p5p+In1LgH+KmZXRDnao2ElPmMIWU+54HkG0knI94ViBMHNEStywZa4lCXEWNmc4BK59zLZtbXNUjK62NmbwVWOef2R23q7/0m3bUILAIOOefWBM+/CdxM6l6PXOAB/H8Cd5nZVlL3WiSapIxJfXHOHQZuClo2Pwk8GucqxVrKfcaQvJ/zIPKNpJOqiXUFUBq1rhDYHYe6jAgzSweuAT4XrKoA5kVsN/yNFLv72TaafQZY6t8O4L8Dv8N3h+jt++D62DbaZRARA5xzDWa2EUgnNa/H/cB7nHNVwX8ED+D/o0vFa5FoUi5mR/gG8Jd4V2IEpPJnDEn0OQ8y30g6qdoV5ClgspnlAgTDvkwHno5rrWLrM8Adzrmm4PmfgfkR22fib8Za1c+20ewK4PiIshv4MP5O9N6+D8n8XVkDFJvZhIh1rfjrklLXw8zGAcc556qCVbcC1cB2UuxaJKhUvtbtwIvxrsQISOXPGJLrcx5MvpF0UjKxds7tBR4Gzg1WnQc8FNGhPqmY2eeB1UCumR1jZh8EqoA9ZrY42O184B7nXLNzblVv20a88iFyzlU653Z1FPwd2pXOuR308n1I5u+Kc249/mfHSwHMrATIAn5K6l2PQ0BjcLd6h4P4Pz5S7VoknFS61mY22cyuMLO0oGXv0/jhQZOKdf10aJD8n3H0+03Wz3mw+Ua86hlLKTncHoCZlQK3A9uAicB/Oudq41qpGDCzm/BD3URa75xbYGYz8EPObQMKgBsjhsnpdVuyMLNtwJXOub/09X1I5u9K0FJ7N77lYDr+Rsa1qXg9zOw4/A2ML+CHnfyrc+7pVLwWiShVrrWZLcL/wVsNPAvcFfwRnDQiRt75NvAl4JvOucpk/Yx7er/AeJLscx5qvpFsUjaxFhEREREJU0p2BRERERERCZsSaxERERGRECixFhEREREJgRJrEREREZEQKLEWEREREQmBEmsRERERkRAosRYRERERCYESaxFJKGZWZma/N7PtZva5eNdHRET6prjdRRPEiEhCMbO3Ak8BRcB6YJ5z7lB8ayUiIr1R3O6ixFpEEpaZ/RU4J1mnvhURSTapHrfVFUREEpKZTQP+lKrBWURktFHcVmItMWBmnzWzb8W7Hh3CqI+ZnWdmfzez5SFVK6mY2R0hH8+AdwO3hXlcEelO8Xr4zGymmX3fzL4X0vHSzexqM9s+gH3fZmbPD/W9hhm7Fbc9JdajjJktMLNfm9kPzOxlM3Nm9q/xrleUR4AH4l2JCGHUZydwXAh1SUhmdqGZHTCzPWZ2mplNC75nzsy+aWYlwX4nm9mrZnanmeUG6yYAe0Ou0hXAD51zrcHxRUYdxeshGY3xen/wmB3S8dKAKmD6APZdBSzobaOZLTKzXwTfvZvMbGLEtrBjt+I2SqxHo98C/+2c+zCwFPhJxwYze4OZnRWPSkWe2zn3qnPu2XjUoyeDqU9v19A5twE4EGa9wvy8hnss59wjwLeBVufc351zO4GrgWbgeedcVbDfP/GB/HrnXEPw8ncADw/rDUQws/8Cvgg8bmbrgflhHVtkhCleD9JojNfOuXpgd1jncs61AC8McN99+CS8t+2vAeXB0zucc5GJdGixW3G7ixLrUcTMxgDzgDoA51w78AXAmVkxPmhbHOoVt3OHaQDvo30EzxWPY/0CmGpmJwI45w4CTwBvjzhXCbA7+O51ONY5t36Y5+7knLvFOTfXObfMOXesc+6vYR1bZKQoXsdWAsbrsEeCGEz9+9u3LXiM7vccWuxW3O6ixHoUCYau2QA8ZmanBuu2A//E92uaBawws6sAzOxLZvZfwc9Anf2ozGyZmX3NzD5pZvVmtsHMPmVmF5vZSjN7l5n9xsyqzeyaiNf1eLyoc19tZj+K7GsW/Bx6j5l9xcyeMrOTg/V9ni/i9WZm15jZa2b2ZjNbY2b7zOzKAZxjbmR9+jnnUdewB6XBT7s1ZnZ11DX9YnBtfh0E4oFes6POZWaXmtk/zexKM9toZrf1cbyePvse69OXoGXjNeCyiNVlwNss6PYBXAT8PqKeuUB9xPNjzew2M7vDzJ4236XkcjNbZWbvNLMHzKzKzK4zsxOD93jQzC7sr34io4niderE65709h77Oo+ZzTKze823/v7vQM7Tw3kfMbOfmtkZ/eyn2B0rzjmVUVSARcB2/F+o9wLjI7ZtA5YHy4uBLcFyGf6v6bH4v7h3AMuCbd8CnguWs4BDwDfxY1G+F6jo63jR58b3Mfse8KNgfQ4+WSsJnl8KVAKFfZ2vh/c9MTjn54N9bwUagRn9nCO6Pn2eM/Ia9lCHbcCdQF7wuh3B+iLg5xH7/Q0fFAd0zXo51xj8z3t3AScD5/V1vKjPvsf6DPD7dROwOVieDXwFaAIuifi+WMT+FwFvDJbzgJVAXvD8r8F1KAKOAF8CcoHLg+fvxv9x/2mC76CKSjIVFK9TIl4H2784wOvY12e9ElgQ8Ro3wO9Zx+d5GnBzD9uXB+fJiFin2B2johbrUcb5VsUl+EDzAeA1M+vpJo2NwPvNLBN4Y7CuAP8PeRrQ0T92Hf4fDM65ZqAGeNA5V43/hzaxn+NF16+J7jdDvA1ocEEfXeDX+H+Ql/RzvujjdhzzZ8G+t+L7/765n3N0q89gztmLPzjfn24lMCniPY4NWpQ+ib+mMMBr1sv7PYQPYL9zzv3TOff4II7XW30G4hfAMea7g/wL8H3g/wHvNrMCoMYFkTdwBvBcsHwR/j+vjlaQdwG3BNe5CnjS+X7Zq4Ai59wvnf95fDWD+wxERgXF69SI10/Fj7sAAAS6SURBVD3o9T32cZ7zgGLnXEd9Vg3ynCcB73PO/fcA91fsjpGMeFdABs7M0oFpzrltwCfM7If4O6jLgVMi93XONZnZFOAG4Mcdh3DO7TezVcC5+L+oZ+GTqc6XRiy3EnQX6u14vVQ18hhzgcyIerWb2TZgal/n649zrtHMtgAlwLhBnGPI54zSSte/n+nAVufcXdE7DeSamdn78Qlsh/nO/2TsIus6iM+g1/r0xzm30cxW47uDFDvntpjZL/CtSJcBf4yod8d3o6P/3gx8C1PHsQ5GHjpiuY3u2lG3NEkyitfd3l8qxOtIvV7HPj6bN+D/iBiqc4C5ZlbgnKvtZZ/24D0odseQLsjokklE/1fn3Mv4n9qOjd7RzM4BrnXO3er8CA+RPgqcZH7YpzXAV/s7cT/H68sOfAto5DBEhu97OFxZwOsxPsdA7KF7P2TM7JRBXLPfA8dHlB7vLh/E8Xqsz8DfDj8HrgK2RNQvHbgO/7Nph1Pw/UU7VABnmll+xHn77OcnksQUr7tLpXjd63vs4zzVwWuyGJo78K39d/ey/YjruulcsTuGlFiPPh83s8kRz6cAjwXLzfibNY7FD+1UbGbZZnZesH28mY3D9/X7If5nnOcIfloMpBH1V7qZWT/Hiz63RRzjt0At8J7gWCXB+TpaPns7X28mBftMCV772ADOEVmf/s4Z+T6i9dR6YfhWqHzgYTN7q/mbZPr7DDrPBUxyzq2PKC0Rp0iPWO7reJH17rE+ZpZlZuVBN4++/BL/E/SDAM65GuBR4G+u+2ggb8F3E+nQcb0fMD8W9seAjrFMoz+D6M95VI9QINILxWtSJl4P9Dr2eB58LM3C//EF/tcJzGzyAGK34VuP3w9cZmbvDl77aTO7Cd9a/0TE/ordsRSLjtsqsSn4GyIccBi4H/gR/qek0mD7zcAu4J34n3Y24lsI3oZv6fgRPkj9Hn8jRXNwvHr8P7QLgBb8DTLTgP8Mtn+wr+NFnfs64Bl8v7VFwbZl+JbOO/F9DRcH63s9Xy/v3wWvvw7fNWFJxLbezjEzsj79nTPyGkad+534m/i+GVyLjtddFWw/C3g1+GzuDK7zQK/ZO3t4r5cG53sAmBKs6+sz7XasXuozFT+Rwc8G8F37RtTzy4Gzo9bd0cPr3gSsBQ4CNwXrLsR/1+4BJhMMOQb8O/4/3u8H7/Xt8f43pqISVkHxOpXi9XTgWXyr+/H9vMe+Puu34X8p/Af+RvLXgWuD995j7Mb/KtKCb6keAzyJT+o/DFyDH8/7Tnw3lI7XKHbHsFhw8SRFmNl44N+dc18KnqfhW1Gucc7dFNfK9cPMHDDL+T6LMkRmdr1zbkhDOYnIyFG8lkiK3aODuoKkno/huwWMiVh3DhE3pSWiiJ+e9LPTEJlZupm9D3go3nURkQFRvBbF7lFGiXXq+S6+L9Y6M9sJ/BnY5Zz7e3yr1bvgZo7rgqdXmh/2TQavAPitc25Lv3uKSCJQvBZQ7B5V1BVERERERCQEarEWEREREQmBEmsRERERkRAosRYRERERCYESaxERERGRECixFhEREREJgRJrEREREZEQKLEWEREREQmBEmsRERERkRAosRYRERERCcH/B9pzC98/S4NiAAAAAElFTkSuQmCC\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/stardust-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/stardust-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/stardust-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 }