{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 27 - OREX - 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 OREX entry probe.\n", "http://adsabs.harvard.edu/full/1995ESASP.367..271I" ] }, { "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('OREX', 761 , 56, 0.00, 9.08, 0.0, 1.35, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,7.43,0.0,-3.17,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": "iVBORw0KGgoAAAANSUhEUgAAAtgAAAHtCAYAAAAjjAmGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOzdeZhU1bX38e9qZpkEm0EmmxmUUQERRVFRo4aIxiEao0lEcmOi3kzXJEaNUa8ZiEluYogkMQ5xikOiCMrrhII4ACIgMsokg9CoyChN0+v9Y1d1V090N13Vp6rr93me/VTVOafOWdXi6tW79tnb3B0REREREUmOnKgDEBERERGpT1Rgi4iIiIgkkQpsEREREZEkUoEtIiIiIpJEKrBFRERERJJIBbaIiIiISBI1jDqAZMnNzfW8vLxqH5+fn0+7du1SF1A9oZ9T9elnVX36WZU2f/78be6eVT+QmubsuGz9t6PPnV30udNbdXJ2vSmw8/LymDdvXrWPHzZsWI2Oz1b6OVWfflbVp59VaWa2LuoY6lpNc3Zctv7b0efOLvrc6a06OVtDREREREREkkgFtoiIiIhIEmVtgT1x4sSoQ8gI+jlVn35W1aeflRyqbP23o8+dXfS5M5+5e91e0OwM4HbgEndfG9t2AfBboBXwT+AH7l4Y25cL3AFsBxoBN7j7/rLnHTZsmGfCuB0RkbLMbL67D4s6joooZ4uIlFadnF2nPdhm1gFoAYxI2NYNGA9cCHwPmAD8d8LbHgfudfcbgBXAL+osYBGRLKacLSJyaOq0wHb3LcDTZTYfBXzT3ee7+wPA3cCpAGY2Eujr7m/Fjp0KXGtmzesqZhGRqrg7hYWFHDhwIOpQkiqdcrY7HDgAhYW1PZOISOrV+Rhsdy8q83pW/KvFmI3A+tjz04A1CcdujD1Ny69SRSQ7TZgwgUaNGnH//fdHHUrSpUvOfv11aNgQxoyp7ZlERFIvHW9yHA78Jfa8M/Bxmf07gU5l35Sfn8+wYcOK25QpU1IcpohI0KhRIwAKCgqq/Z4pU6YU5ysgNzWR1Yk6ydk5sd9WRUUHPUxEJCVqmrPTaqEZM+sF5Lv7wtgmB/aWOawJUO6GmXbt2mXE5OQiUv80btwYgH379lX7PRMnTiy+Y97MtqUksBSry5ytAltEolTTnJ02Pdhm1gD4FvDjhM0bgTYJxxjQHNhUt9GJiFSuSZMmQM16sDNdXedsFdgikknSpsAGfghMcvfELqAXgL4Jr/OAAkBd1SKSNg6lB7seqNOcHS+w69l9pCJST9V5gR3r0QCwhG0/BRYAzcysh5l908x6ufs8YLOZDYgd+gXgj+6ePd1EIpL24gV2fezBTpecrR5sEckkdToG28xaApfHXn7dzP5E+IrxtjKHLnP3e2PPLwFuNrO1hPlYb6yLWEVEqqu+DhFJp5zdoEF4VIEtIpmgTgtsd98JTI61uNtjrbL3rAOuSnFoIiKHrL4OEUmnnK0ebBHJJOk0BltEJCPV1x7sdKICW0QyiQpsEZFaqq892OlENzmKSCZRgS0iUkv1+SbHdKEebBHJJCqwRURqSUNEUk83OYpIJlGBLSJSSxoiknrqwRaRTKICW0SkljREJPVUYItIJlGBLSJSS/EhIurBTh0V2CKSSVRgi4jUknqwU0+ziIhIJlGBLSJSS7rJMfV0k6OIZBIV2CIitaSbHFNPQ0REJJOowBYRqSUNEUk9FdgikklUYIuI1JKGiKSeCmwRySQqsEVEaklDRFKvUaPwqL9hRCQTqMAWEakl9WCnXrNm4XHvXnCPNhYRkaqowBYRqSX1YKdeo0ZhJpEDB2D//qijERE5OBXYIiK1pJsc60ZiL7aISDpTgS0iUksaIlI3DjssPKrAFpF0pwJbRKSWGjZsCEBhYSFFmuYiZdSDLSKZQgW2iEgtmZl6seuACmwRyRQqsEVEkkA3OqZefIjInj3RxiEiUhUV2CIiSaAbHVPv8MPD46efRhuHiEhVGtb1Bc3sDOB24BJ3Xxvb1gz4NbATaAPc4O47qtonIpIu6usQkXTK2bm54fHjj5NxNhGR1KnTHmwz6wC0AEaU2fVnYLa7/xR4GrinmvtERNJCfRwikm45+4gjwuO2bck6o4hIatRpge3uWwgJt5iZdQIuA6bGNr0IfMnM8g62r04CFhGppvrYg51uOVs92CKSKep8iIi7F5lZ4qYxwDZ33xPbX2hma4BTgP0H2be2LuMWETmY+tiDDemVs+M92Pn5tT2TiEhq1XmBXYHOQNn+iJ1AJ6DwIPtERNJGTk74QtDdI44k5SLL2V26hMcPP0zG2UREUicdCmwHys5q2oTQE0IV+4rl5+czbNiw4tcTJ05k4sSJSQxTRCR5pkyZwpQpU+Ivc6OMpYYiy9l5eeFxzZoaRCsikgQ1zdnpUGBvJNxpnqglsImQyCvbV0q7du2YN29eSgIUEamu6vZgJxaUZpZJt+1FlrO7dw+Pa9aAO5QeuSIikjo1zdnpMA/2K0Cn2NROmFkjoBvwahX7RETShmVPtRdZzm7TBlq1gl27dKOjiKS3Oi+wreS3kAG4+0fAs8Dpse1jgSfdfePB9tVhyCIi1VbfxmCnU842g759w/P33kvGGUVEUqOu58FuCfxX7OXXzaxd7Pm3gfPM7EbgHCBxIN7B9omIpIX62IOdjjn72GPD4zvvJPOsIiLJVadjsN19JzA51hK3fwpcXcl7Kt0nIpJu6lMPdjrmbBXYIpIJ0mEMtohIxquPPdjpKD7xyOzZ4UZHEZF0pAJbRCSJ6lMPdjoaPDgsOLNuHaxaFXU0IiIVU4EtIpIE+/eHqZ4bNWoUcST1W4MGMHZseD5jRrSxiIhURgW2iEgS7N0b1ldp1qxZxJHUf+ecEx4ffzzaOEREKqMCW0QkCT7//HMAmjZtGnEk9d/550OzZvDaa7B2bdTRiIiUpwJbRCQJ1INdd1q2hPHjw/O//z3aWEREKqICW0QkCVRg161rrgmPd98dVnYUEUknKrBFRGrJ3TVEpI6ddBKceCJ8+mkoskVE0okKbBGRWtqxYwdFRUW0aNGChg3rdP2urHbzzeHxjjvgo4+ijUVEJJEKbBGRWsrPzwegXbt2VRwpyXTmmXDuubBzJ/z4x1FHIyJSQgW2iEgtxQvs3NzciCPJPr/7HTRuDPffD88+G3U0IiKBCmwRkVpSD3Z0evcOQ0QArroKtmyJNh4REVCBLSJSa5s3bwagQ4cOEUeSnb7/fRgzBrZuhQsugNj9piIikVGBLSJSSx988AEAPXr0iDiS7JSTA488Al27wpw5oSf7wIGooxKRbKYCW0SkluIFdq9evSKOJHt17AjPPAPNm8PDD8OECSqyRSQ6KrBFRGpp1apVAPTs2TPiSLLbkCHhRsfDDoP77oMrroB9+6KOSkSyUVIKbDM7JxnnERHJNAUFBSxbtgyAPn36RBxN9dTnnD1mDEyfXtKTffrpELsHVUSkztSowDazL5jZbDNbZWarY20t8GRqwhMRSW+LFi2ioKCAvn370rp166jDKSVbc/Ypp8Ds2dC5M7z+OgwbFsZmi4jUlZouOfYgcCewACiKbTPgsmQGJSKSKebOnQvA8OHDI46kQlmbs4cMgbffhvPPD48nnwy33AI/+QlosU0RSbWappnV7n5X2Y1mtixJ8YiIZJTXXnsNgOOPPz7iSCqU1Tm7UyeYNQtuugl+/euwtPqTT8I990B6/ucSkfqipmOwv21mF1aw/UvJCEZEJJMUFhYyY8YMAM4666yIo6lQ1ufsxo3hV7+CF16AvDxYuBBOOAG+/W0tSiMiqVPTAvtE4G9mdiChFQGTUxCbiEhae+utt/j000/p2bMnvXv3jjqciihnx4wdC0uWwA03QIMG8Je/QM+e8POfw86dUUcnIvVNTQvsXwATgd5Aj4T289oGYmb9zexuM/uemf3ZzIbEtjczsz+a2f+a2WQza1Xba4mIJMMjjzwCwLhx4yKOpFLK2QkOOwx++Ut4910YNw5274Zbb4UePeC22+CTT6KOUETqC3P36h9s9jDwDXffV2Z7R3f/qFaBmM0DznP3jWbWDZjh7v3N7B/A8+7+mJl9AbjS3S8t+/5hw4b5vHnzahOCiEi1ff7553Tq1IlPP/2Ud999l8GDBx/yucxsvrsPS2J48fMqZx/ErFmhR/uNN8Lr5s1h4kS4/no46qhIQxORNFadnF3THuxHgZ+Y2ckJ7RTgJ4ccZYn+QMvY871AazPrRLjbfWps+4vAl8wsLwnXExE5ZI8//jiffvopQ4YMqVVxnWLK2QcxenSYxu+ll+DMM0OP9u9+B927wxe/CFOnajVIETk0NZ1F5CbCV41fB+Jd3wZ0BK6vZSxPEMYKng1cDlwLjAG2ufseAHcvNLM1wCnA2lpeT0TkkBQVFXHnnXcCcN1110UczUEpZ1fBDE47LbQFC2DSJHjiCZg2LbSuXeHKK+Gyy6B//6ijFZFMUdMe7N8AXd09z927x1oecEESYrkG+ByYC+xy9yeBzsDHZY7bCXQq++b8/HyGDRtW3KZMmZKEkEREynv88cdZunQp3bp14/LLLz+kc0yZMqU4XwG5SQ2whHJ2DQwdCg89BBs3wm9+A716wYcfwu23w9FHh/2/+U3YJiLZpaY5u6ZjsE9w9zcq2P5Nd7+3RpGWP0cucB4hQd8Qez4EuMjdj0847h3gYXeflPj+dBjPJyL13+7du+nfvz8ffvgh99xzDxMnTqz1OVM4Bls5uxaKimDmzLDk+hNPwGeflewbNgy+9CU47zwYODD0hItIdkjFGOxrKrjIRcDva3ieivwTeNLdfwH8EngE2Ay0KXNcS2BTEq4nIlJjt956Kx9++CHHHnssV111VdThVEU5uxZycsLQkb/9LcyZ/e9/w4UXQrNmMG9eWLhm8OAwC8n118Pzz8OePVFHLSLpoKYFdlszuwnAzA4zs78BfyN8RXjIYj0hg919e2zT7cAOYB3QycyaxY5rBHQDXq3N9UREDsXLL7/MpEmTyMnJ4e6776ZBgwZRh1QV5ewkadIExo+Hxx+HbdvgmWdgwgTo0AHWroX/+z84+2xo0wZOPRX+93/DEu26SVIkO9XoJkd3P9fMRprZH4GzgI+AwcBnB39nlT4BPjezzu6+MbbtY2AR8CxweuxxLKHHZGPFpxERSY1Nmzbxta99DXfnpptuYuTIkVGHVCXl7NQ47LAwj/a4cWEYydtvh4L7hRdg/vwwrGTmTLjxRjj8cDjllDBjyUknwbHHQqNGUX8CEUm1g47BNrOTK9k1AjgVmAQUAePd/Xu1CsRsMOHrzPlAB+A1d3/VzNoAvybcgd4R+Im77yr7/kwYzycimWnXrl2cfPLJLFiwgNGjR/Pyyy/TsGFNJ2GqXLLGYCtnR++TT+Dll0Ox/eKLsHp16f3NmsHIkSUF98iR0LJlxecSkfRUnZxdVYH9JmGu0+2EpFzhYcCR7t7kUANNhvqarEUkWnv27GH8+PG88MIL9OrVizfeeIPc3ORO+pHEAls5O82sXh0WtJk1C2bPhuXLS+/PyQnT/w0fXtIGDQpDUkQkPVUnZ1fVBfMLYJ67b63iQhfWNDgRkXS3e/duxo0bxyuvvEK7du2YPn160ovrJFPOTjM9eoR25ZXh9datYXGb2bNDe+cdWLIktPvuC8c0bhyK7OHDYcSIMKykf38NLRHJJDWapi+dZUtviIjUjY0bNzJ+/HjmzZvHkUceyUsvvUT/FK00kqpp+tKZcnbw+efw7rswd25JW74cyv5qbtQIjjkmzFoyZEh4HDwY2raNJm6RbJaMHmwRkazz9ttvM378eDZv3kz37t2ZMWMGvXv3jjosqYeaNg3jsBPvmd2xI9ws+fbboeB+91344IPw+O67cP/9Jcd27VpScA8aFIrwXr1CL7iIREcFtohITFFREZMmTeJnP/sZ+/fv55RTTuGJJ55I92EhUs+0ahWm+jv11JJtO3fC4sWhwF64MDwuXhxWlfzwQ5g6teTYhg2hT5+w+uQxx5Q89u6twlukrqjAFhEB1qxZw4QJE3j55ZcBuPbaa5k0aRKNVZFIGmjZEkaNCi3uwAFYtaqk6H7vvTCWe80aeP/90J54ouT4Bg1CkR0vuvv2DYV4795hOkERSZ5aF9hmNhT4AbCAMN/p2tqeU0SkrhQUFPDb3/6W2267jb1799K+fXv+8Y9/cM4550QdWkooZ9cfDRqEIrlvX7jkkpLte/bAsmWhwF6ypORx9eqwfdkyePLJ0udq1y4U24mtd+8w3KRZs7r9XCL1QTJ6sP8LGOvul5vZZWa2xd33JuG8IiIp4+5MnTqVH//4xyxduhSAr371q9x11120b98+4uhSSjm7njvssDDzyLHHlt6+d29J4f3++7BiBaxcGR7z80N7/fXS7zEL47zjBXePHtCzZ8nsKJrDW6RitZ5FJLYUbnt332hmOUALd9+RlOhqQHeki0h1vf7669xwww28HqsmevfuzeTJkzn99NMjiacuZxFRzpayiopg06ZQaMdbvPBevRoKCyt/b25u+aI73jp3Dr3sIvVNSmYRMbOzgbbu/pCZjQB2ufv7AO5eBNR5ohYRqYq7M3PmTO68805eeOEFAI444ghuuukm/uu//osm9XRlD+VsqUpODnTpEtppp5Xet38/rF1bUnSvWROK7g8+CM+3bQvt7bfLn7dxY8jLK1949+gB3buHmzlF6qsaFdhmdgdwPfA88JC7v21mvzSzRe7+cEoiFBGphQMHDjB16lTuvPNO3o5VAc2bN+d73/seP/zhD2ndunXEEaaOcrbUVqNGYWhIRbNUFhXBRx+FgjtedMefr14d9sV7xCvSpk0owCtrKsAlk9W0B3sU0AW4MmHbc8CDgJK1iKSNrVu3cu+99/KXv/yFdevWAaHH+vrrr+c73/kObbNjhQ7lbEmZnBzo1Cm0k04qv3/37pIe78T2wQewbh18+mloCxZUfP7EAvyoo8oX4PX4b2OpB2paYL/g7tvNLHHgtpbcFZG0UFRUxOzZs5kyZQqPP/44BQUFAPTs2ZPrrruOq666iubNm0ccZZ1SzpbING8OAwaEVpZ7uKly7drKW1UF+OGHV977fdRRoQA3S/rHEqmWmhbYq83sp0BnM/sKcBFwPnBt0iMTEammVatW8eCDD/Lggw+yZs0aAMyMcePGcc0113DmmWeSk5MTcZSRUM6WtGQG7duHNmJE+f2JBfi6dRUX4Nu3l6xuWZEWLaBbt9C6di15Hn/dpQvU01svJA3UqMB290fNbCTwdeAKYDNwpru/mILYREQq9dFHH/Hvf/+bf/7zn8yZM6d4e5cuXfja177GxIkTycvLiy7ANKCcLZmqOgX4tm2V936vWwe7dpVMSViZjh1LF99ln7dvH4bCiNRUjWcRcfc3gTfjr82sgZmd7u4vJTUyEZEyNmzYwFNPPcWTTz7JrFmziE8z2rx5c7785S9zxRVXMGbMGBpobrBiytlSH5mFxXHatYPhw8vvdw893OvXh6Xk168v/3zjxnAj5kcfwdy5FV+nceOSoruiQrxrV80FLhU7aIFtZs8ALao4Ry6wDVCyFpGkcneWLVvGtGnTeOqpp3jjjTeK9zVu3JgzzzyTiy++mPPPP58WLapKVfWfcrZIYBZukmzTBgYPrviYwsJQXFdUfMdff/xxuCnzgw8qv9bhh5cU24mFd/x5586hUJfsUlUP9qLYMcuBIuBqwvK6iasDdAey4nZ8EUm9zz//nJkzZzJt2jSmTZtWPKYaoGnTppx99tlceOGFfPGLX6SV5vEqSzlbpJoaNiyZ/3vUqIqP2b27pPCurCd8+/bQFi2q+BxmJUNRKirAu3aFDh00FKW+qarA/g1htcftAGbW091vLnuQmT2diuBEpP5zd1auXMmLL77I888/z0svvcSePXuK9+fm5nL22Wczbtw4zj77bPVUH5xytkgSNW8O/fqFVpH4zZgffli+EI9v27QJNm8OraIFeSDMN96lS+UFeNeuoadcs6JkjoMW2O7+WZlNHcysgbsfiG8wszFABbcgiIhUbNOmTbz00kvFbcOGDaX2Dx06lHPPPZdzzz2X4cOHa0x1NSlni9StxJsxjzuu4mMKC0ORXbbwTny+bVuYMzzhC7tyWrQoX3gnFuBdu0KzZqn5nFJzNb3J8QlgoZnNBnYC/YGzgHI9JCIicVu3bmX27Nm88sorvPTSSyxdurTU/nbt2nHaaacxduxYzjnnHDp16hRRpPWOcrZIxBo2LLk58sQTKz5mzx7YsKHyXvD168OsKEuXhlaZ3NzKe8G7dYMjjwzxSOrVdJq+F8zsTOAaYBCwHbjQ3fV1o4gAYcjH2rVrmTVrVnFbvnx5qWOaN2/OySefzNixYzn99NMZOHBgts5TnVLK2SKZ4bDDoE+f0CoSnxXlYL3gGzaEnvBt2ypfnCe++ubBbsrMzdVQlGQ4lGn6NgE/i782sxwzO87d5yc1MhHJCAcOHGDJkiXFxfTs2bPZuHFjqWOaNWvGyJEji4vqESNG0Fi31dcJ5WyRzJc4K8qgQRUfU1QEW7ZUXoCvXx9mTdmwIbTKNG0axoNXNAxFUxNWX40KbDO7t4LN7YADwPhkBBRbFOEkYBUwG9gN/Jrw9WYb4AZ335GMa4lIzW3evJm33nqLt956izfffJN58+axa9euUse0bduWk046idGjRzN69GiGDh2qgjoCytki2SMnJwwBOfJIOP74io8pKAjzfx+sCN++HVatCq0yrVtXXHjn5UGPHiGGbP9SsqY92L0Jc6d6wrYuwKxkBGNm3wB6untib8s/gOfd/TEz+wJwD3BpMq4nIge3d+9eFixYwJtvvllcUK9fv77ccXl5eZx44omMHj2ak046if79+2vIR3pQzhaRYo0bQ/fuoVVm166DjwX/8EP47DNYvDi0ijRtGgrtnj1LHuMtLy87lqivaYE9wd1LDaY0s36EZXhrxcxOAb4KnJGwrRNwGfCd2KYXgSfNLM/d19b2miJSoqCggPfee4933nmH+fPnM3fuXBYuXEhhYWGp41q2bMmIESM4/vjji1uHDh0iilqqoJwtIjXSogX07x9aRdzDAjwVFd5r1oRFebZtq3yZerMwBCWx6O7dO0yF2KtXKM7rg5re5Li8gs37gAnAj2sZy13AHOCPZtYLuA04Ctjm7nti1y80szXAKcDaWl5PJGsVFBSwePFi5s+fX9wWL15MQUFBqeNycnIYOHAgI0eO5Pjjj2fkyJH069dP0+ZlCOVsEUk2s3AjZG4uHHtsxcfs2AGrV5esghlvq1fDunUlhfnMmaXfl5MTetf79QsFfnwO8n794IgjUv7RkqqmY7DXUPqrxgZAB+DB2gRhZn2BocA33H2Rmf0IeB64Hfi4zOE7gXJzeOXn5zNs2LDi1xMnTmTixIm1CUukXti3b1+FxfT+/fvLHdunTx+OO+64Uq2l7mZJiSlTpjBlypT4y9xUXEM5W0Si0KoVDBkSWln794ciO7HwXrkyTD+YWJRPm1b6fbm5MHAgDB5ccu7+/etuGfqa5mxz96qOKTnYbArwMCUJuwjY4u4rahhn2fNeAExx99zY62ZAPvAL4MvufnzCse8AD7v7pMRzDBs2zOfNS1wNWCT7bNu2jYULF/Luu+8WPy5durTcMA8zK1dMDx06VEuPR8TM5rv7sKqPrPF5lbNFJGPs2xdurly2LBTcy5aVtN27yx/fqBEcfXQotkeMCEveDxwIqf6StTo5u6ZjsP/p7q+VuUhfM+vl7ge537RKDRNjcfe9ZraS0NvSpsyxLYFNtbiWSMY7cOAAq1atKlVML1y4sNz0eBCK6f79+5cqpocMGaKe6eygnC0iGaNJEzjmmNASuYepBRctgnffDW3hwtDzvXBhaPffH45t0QJGjgyL+pxxRphRJYrFdWp6yS8Ar5XZtgV4hfB14aFaBLQ2sw7uviW2rZCQlDuZWbNYAm8EdANercW1RDLKrl27WLx4cale6cWLF7Nnz55yxzZv3pzBgwczePBghgwZwuDBgxkwYADNmzePIHJJA8rZIpLxzEqmBDz33JLtO3eGmUwWLIA33oA5c8KNli++GNqtt8Lhh8NZZ8F558H48XW3nHy1Cmwz+y5hbF1LM7uhgkNeqE0Q7r7MzJ4DLgTuNrPDgcbAQ8DZwOnAs8BY4El3L99NJ5LhDhw4wOrVq1m8eHFxW7RoEatWraKioVxdu3YtVUgPGTKEHj16aHo8Uc4WkazQsmUYFjJqFHwnNnfR5s2h0H7tNXjuudDL/dhjobVuDV/5CvzgB2HmklSq9hhsMxsEXAn8X5ldBe6+udaBmOUCfwDmEXo8/uru75tZG8KiBWuBjsBP3H1X2fdrPJ9kkq1bt5YrpJcsWcLevXvLHduoUSOOPvroUoX0oEGDOCLTbqmWSqViDLZytohIGNM9bRo89BDMnRu2NWgA11wDv/71oU0LWJ2cXdObHBu5e7mpB8ysm7uXX32iDilZSzras2cP77//frlieuvWrRUe37lzZwYOHFjcBg0aRP/+/bUKYj2XwpsclbNFRGLeew9+9zu4776wtPxpp8GMGTUfo13rmxzN7IfAdHePTxX+P2ZW9t5MA04FxtQsPJH6o6bDO1q2bMmAAQNKFdIDBgygbdu2EUQv9YVytohI5QYMgL//PQwnOftsePlleOYZuOCC5F+rqpq9O3B4wutjgebAR5RM+2SErwFFskJNhnc0aNCAvn37liqkBw4cyFFHHYWZRRC91HPK2SIiVRgyJMxUsnVraKlw0ALb3b9TZtMdwHvuXmq5NzM7OtmBiURt7969xcM7Fi1aVKPhHfFCul+/fjRp0qSOI5dspZwtInJw69bBtdfCK69A8+YwblxqrlPTafoGuvs7iRtiK3r1BCpYcV4k/bk7a9euLS6g48X0ihUrKCoqKnd8fHhHvIgeOHCghndIulLOFpGsd+BAKKjvvx8efRQKC8Nqk088AZ07p+aa1Z2mbzTha8UzYkvvJmpLuEt9apJjE0m6zz77rFwhvXjxYnbu3Fnu2AYNGtC/f//iQlrDOyRTKGeLSLbbujXMhTxnChYAACAASURBVD1jRmhbYjP25+TAZZfBr34FXbqk7vrV7cHeDjxAGN93Ypl9BcDfkhmUSG0VFhayYsWKUoX0okWLWL++4okTOnTowKBBg0oV0/3796fpoczfIxI95WwRyRr794cFZ958s6StXFn6mF694KtfhSuugB49Uh9TtQpsd19sZmOAi9xdiVnShruzZcuWUkX0okWLeP/99ykoKCh3fNOmTTnmmGOKi+h4Qd2+ffsIohdJDeVsEamvtm8vv2T6e+9B2V/5zZrBySfDmWeGJdMHDAgrQtaVao/BdvfPqKTXw8xOcvfZSYtKpAJ79+5lyZIl5YZ45OfnV3h89+7dyxXSvXr1omFNJ7wUyUDK2SKSyQoKYMUKWLIktEWLQjG9dm3Fx/fqBSecACNHhjZwIDRqVKchl1LVPNjPAC2qOIcBxwDqApSk+eSTT3j33XdZsGBBcVu2bFmFNx22atWq3PCOAQMG0KpVqwgiF4mOcraIZJrCwrDa4nvvlRTTS5aE4rqwsPzxTZuG3ughQ2Dw4PA4aFC4aTGdVNWVt5wwlm8jJXOolmXAF5MZlGQPd2fjxo2lCukFCxawbt26csc2aNCAo48+ulSv9KBBg+jatatuOhQJlLNFJC0VFMDq1bBsWelievny8sM7IAzn6NkzzFd9zDElRXWfPjVfeTEKVYX4K2Cnu+872EFmprvRpUpFRUWsXLmyXDG9bdu2csc2a9aMQYMGMXTo0OI2YMAAmjVrFkHkIhlDOVtEIlNUBBs3ht7neFu+PDyuWRP2V+Soo0qK6HhB3b8/HHZY3cafTFUtNFO+8qnY+eiudElw4MABli9fzty5c5k3bx7vvPMOCxcuZPfu3eWObdOmTalCeujQofTp00djpUVqSDlbRFKtoCAs1rJmTeiRjj+uWBFm7qhgUWMg9Eh37x56oBN7pfv3h5Yt6/Yz1IVaVTBm1g34DvBdlKyzVnyhlrlz5xa3+fPns2vXrnLHdunSpVwx3a1bNw3xEKkDytkiUpWiIvjoo9IFdOLzDRvAKxuABnToEIrosq1HjzB+OlvUuMA2swbAl4CrgTOALYQxf5IlPvroo1LF9Lx58yoc5tG1a1eGDx/O8OHDOe644xgyZAjt2rWLIGKR7KWcLSKJCgtDkbx+feiJjrfE159/Xvn7c3KgW7fQG92jR3iM90z37g2HH153nyWdVbvANrPuwATgG8ARwDpgHDADOCEl0Unk9u3bx4IFC5gzZw5z5szhrbfeYsOGDeWOy83NLS6m461Dhw4RRCwioJwtkq127z548bxxY+VjoeNyc8sX0PHn3bpFO/1dpqiywDaziwg9H6cBq4C7gPuBn7r787HD5qQsQqlTW7ZsKS6m58yZw/z589m3r/T9Ui1btuS4444rVUxr+XCR9KCcLVJ/ucPHH1dePK9bF/YfjBl06hRuLOzWLTyWbfVxTHRdq04P9l7C0rpzgAnuvgLAzA4yAkcyQVFREe+99x6vv/56cUG9evXqcscdc8wxnHDCCYwaNYqRI0fSt29fcnJyIohYRKpBOVskQx04AJs2lS6YyxbTe/Yc/ByNG0PXruWL5ngx3aULNGlSN58nm1VZYLv7s8CzsZtjrjazPOBpoPgLAjMz94MNeZd0UFRUxOLFi5k5cyYzZ87ktdde45NPPil1TPPmzTn++OMZNWpUcUHdpk2biCIWkZpSzhZJX/v2le9xTmwbNlS8uEqiVq0qLpzjrUOHME5aolWTpdLXAzeZWUNgPNDFzG4D/g6cCUxJTYhyqMoW1K+++iqffvppqWO6dOnC6NGjOfHEExk1ahQDBw7U9Hgi9YBytkjd27mz8uJ57dowO0dV2reveNhGvOkmwsxQ40rK3QuBJ4AnzKwX8D3CeD8l6zSwdu1aZsyYwf/7f/+PV155pVxB3bVrV0499VROOeUUxowZQ/fu3TV2WqQeU84WSQ532Lat8gJ63Too8yu3nAYNoHPnyovnbt1A66nVD7XqqnT3VcD1ZjY/SfFIDe3evZuZM2cyY8YMZsyYwYoVK0rt79atG2PGjClueXl5KqhFspRytkjlkjH+uWnTym8cPOqoUFzrS+LskJT/zO7+QDLOY2aNgbnA9e4+08yaAb8GdgJtgBvcfUcyrpXJVqxYwTPPPMNzzz3H7NmzKSgoKN7XqlUrxo4dy5lnnsmZZ55J9+7dI4xURNJRsnI2KG9L5igoKD3+ee3amo9/bt364MM32rcPs3SIpNvfUT8C8hJe/xl43t0fM7MvAPcAl0YRWJSKioqYO3cu//nPf3j66adZunRp8T4zY8SIEZx11lmcddZZHH/88RpDLSJ1SXlb0kK8Bzq+8mC8rV0bHqsz/3OHDgcvoFu3rpOPIvVA2lRiZjYK2Ax8GnvdCbiMsKwvwIvAk2aW5+5rIwmyDhUWFvLyyy/z1FNP8cwzz7B58+bifYcffjjnnnsu48aNY+zYsRxxxBERRioi2Up5W+qSO2zdWr5wjrf162H//srfH1+BMC+v8vHP2bSUt6RWWhTYZtYcuNDdv29mN8c2jwG2ufseCDfqmNka4BRgbSSBplhRURFz5szhkUce4V//+lep5ce7devG+PHjOe+88xg9ejSNtIySiERIeVtSYd++UCyvXBna6tUlhfTatVWPge7YMRTQ8dUH4y0vL8wN3bhx6j+DCKRJgQ3cAPyyzLbOQNn1iHYCneokojq0ePFiHnzwQR599FE+/PDD4u39+vXj4osv5vzzz2fw4MG6OVFE0klW5205dAUFoXBetaqkkI4/X7/+4MM42rQpXzjHnx91FBx2WJ19DJGDirzANrOzgHnuvrXMLiesSJaoCVDhF0D5+fkMGzas+PXEiROZOHFiMkNNql27dvHoo4/y17/+lbfffrt4e9euXfnKV77CZZddpqJapB6bMmUKU6YUz5SXG2UsNZWMvJ1pOVtq5sCB0Ou8bFn5QnrdusqL6Jwc6NEDevWC3r2hZ8/SxbTGQEtUapqzLerFvMzsBWBowqY2wC7gLuCr7t4n4diVwC3u/nDZ8wwbNsznzZuX6nBrbcGCBUyePJlHHnmEXbt2AWHmj0svvZTLL7+cUaNGaRlykSxjZvPdfVjVR6aHZOTtTMnZcnB79sCKFbB0aSimly4NbeXKMNyjIjk5obe5d+/Q4sV0796hiNYwDkl31cnZkfdgE26IaZLw+g3g+8BbwI/MrJm77zWzRkA34NUIYqyVoqIipk2bxl133cXMmTOLt5944olcffXVXHTRRRym77VEJHPU+7wtpeXnlxTQiY/r1oWbDyvSpQv06wd9+pQupLt3VxEt9V/kBba75ye+NrMDQL67rzezZ4HTgWeBscCT7r4xgjAPyb59+7jvvvu46667iheAadmyJVdddRVXX301Rx99dMQRiojUXH3O29luzx54/31YtAgWLy55zM+v+PiGDUPR3K8f9O9f8ti3L7RsWbexi6STyAvsKnwb+LWZDQY6AhkxQG///v3cd9993HbbbcU3LXbr1o3rr7+eCRMm0KpVq4gjFBFJmYzM29mmqCiMkS5bSK9cWXGPdMuW5Yvofv3CGGlNaiVSXtoV2O6el/D8U+Dq6KKpGXfnkUce4aabbmL16tUADBgwgBtvvJELL7xQC8CISL2UyXk7GxQVwQcfwLx5oc2fD++8Azt3lj+2QYNQOA8cCIMGlTx27aoVCkVqQhVfkixYsIBrr72W119/HYC+ffvy85//nIsvvlg3LYqISJ1wD+Oi33orFNLxgnpHBYvVH3lkKJ4TC+l+/aBJk/LHikjNqMCupV27dnHDDTcwefJk3J327dtz5513csUVV6jHWkREUqqwEBYuhNdfh9mzw+OmTeWP69QJjjsOhg0L7bjjwrLgIpIaqgBrYdasWVx55ZWsWbOGBg0acN1113HLLbfQWhN1iohICuzfD2+/DS++CLNmwZtvwu7dpY9p2xZGjixdTHfSUj8idUoF9iEoKirilltu4Y477sDdGTx4MA888ACDBg2KOjQREalH3MN0eC++CC+8ADNnlh873asXnHhiSevXL8w1LSLRUYFdQ9u3b+fyyy9n2rRp5OTk8JOf/IRbbrmFxprUU0REkqCwMPRO//vf8PTTYfnwRH37whlnwJgxoaDu2DGSMEXkIFRg18DGjRs544wzWLp0KW3btuXRRx/ljDPOiDosERHJcAcOhB7qRx+FqVPhk09K9rVrB2PHhqJ67Ngwo4eIpDcV2NW0Zs0aTj/9dNasWcOAAQN45pln6N69e9RhiYhIBnv/fbj/fnjwQdi8uWR7nz5w/vmhDR+uIR8imUYFdjVs2rSJMWPGsH79ekaMGMFzzz1H27Ztow5LREQyUFERTJsGv/0tvJqwiHyvXvC1r8GFF4aFXDTvtEjmUoFdhR07dnDOOeewfv16Ro4cyYwZM7QSo4iI1FhRETz0ENxxByxfHra1bAmXXAJf/zqMGqWiWqS+UIF9EO7ON77xDRYuXEifPn2YOnWqimsREamxF16AH/0ozFkN0K0bXH89TJgA+rUiUv+owD6IyZMn89RTT9GqVSumTZtGbm5u1CGJiEgG2b4d/vu/wzhrCDco/uIXcPnloLXIROov/e9didWrV/ODH/wAgClTptCrV6+IIxIRkUyyYAGcdx58+CE0bQo33xyK7WbNoo5MRFJNBXYF3J1rr72Wzz//nK9+9atccsklUYckIiIZZNo0uPhi2LMHRoyABx4I81eLSHbQxD8VmD59OtOnT6dVq1ZMmjQp6nBERCSDzJgBF1wQiusrr4TXXlNxLZJtVGCX4e7cfPPNANx888101BJZIiJSTQsXhrmrCwrguuvgH/+AJk2ijkpE6poK7DKeffZZ3nnnHTp27Mg111wTdTgiIpIhPvkkFNd794b5rH//e027J5KtVGCX8Zvf/AaA//mf/6GZ7kQREZFq+ta3YM0aOO44uOceFdci2UwFdoLFixcza9YsWrZsyYQJE6IOR0REMsS//gVPPAEtWoRH9c+IZDcV2AkmT54MwBVXXEHLli0jjkZERDLB9u3w3e+G55MmQV5epOGISBpQgR3z+eef89BDDwHw7W9/O+JoREQkU/ziF5CfD6NHw8SJUUcjIulABXbM9OnT2bFjB8cddxzHHHNM1OGIiEgGWL4c/vjHMN76D3/QuGsRCVRgxzz88MMAXHbZZRFHIiIimeL734fCQpgwAYYOjToaEUkXaVNgm9kFZrbGzD42sz+YWcPY9lwzu8fMfmVmd5lZo2Rf+7PPPuPZZ5/FzLRqo4hINUSZs9PF9OmhtWoFt98edTQikk7SosA2s27AeOBC4HvABOC/Y7sfB+519xuAFcAvkn39559/nn379jF69Gg6d+6c7NOLiNQrUefsdFBQEHqvAW6+Gdq3jzYeEUkvDaMOIOYo4JvuXgjMN7MBwKlmNhvo6+5vxY6bCiw3s9vdfXeyLj59+nQAvvjFLybrlCIi9VmkOTsd3H13GH/duzdce23U0YhIukmLHmx3nxVL1HEbgfXAacCahOM2xp4OS9a1i4qKeP755wE4++yzk3VaEZF6K8qcnQ7y8+HWW8Pz3/0OGjeONh4RST9pUWBXYDjwF6Az8HGZfTuBTmXfkJ+fz7Bhw4rblClTqnWhd999l61bt9KlSxfNHiIidWbKlCnF+QrIjTqeWqqznJ0Ofv97+OwzOOssOOecqKMRkbpQ05ydLkNEiplZLyDf3ReamQN7yxzSBNhf9n3t2rVj3rx5Nb7erFmzADj99NMxza8kInVk4sSJTIxNmmxm2yIO55DVdc6O2q5dEFuTjJtv1rR8Itmipjk7rXqwzawB8C3gx7FNG4E2CfsNaA5sStY133jjDQBGjRqVrFOKiGSFKHJ21O69Fz79FEaNCk1EpCJpVWADPwQmufu+2OsXgL4J+/OAAiBp3R7xAvuEE05I1ilFRLJFnefsKLnDn/4Unv/wh9HGIiLpLW0KbDP7KbAAaGZmPczsm8B2YHPsDnWALwB/dPeCZFwzPz+f9evX06JFC44++uhknFJEJCtEkbOjNm8erFwJHTrAl74UdTQiks7SYgy2mf0MuK3M5mXufq+ZXQLcbGZrgRbAjcm67rJlywA4+uijadCgQbJOKyJSr0WVs6P2+OPh8ZJLQL8yRORg0qLAdvfbgQrXwXL3dcBVqbhuvMDu169fKk4vIlIvRZWzo/bqq+Hx3HOjjUNE0l/aDBGJwtKlSwHo379/xJGIiEg627sX3nkHcnJAt+yISFWyusBesyash9CzZ8+IIxERkXT2wQdQWAi9ekHLllFHIyLpLqsL7A0bNgDQtWvXiCMREZF0tnJleOzVK9o4RCQzqMAGunTpEnEkIiKSztatC4/du0cbh4hkhqwtsAsKCtiyZQs5OTl07Ngx6nBERCSN5eeHxw4doo1DRDJD1hbYmzdvxt058sgjadgwLSZTERGRNLV1a3hs1y7aOEQkM2R1gQ1w5JFHRhyJiIiku3gPdvv20cYhIpkhawvsbdu2AdBO3REiIlIF9WCLSE1kbYGdH+uOUIEtIiJV+eST8Ni2bbRxiEhmyNoCO96DnZubG3EkIiKS7nbuDI+tW0cbh4hkhqy9uy8vL48vfelLDBkyJOpQREQkzZ15JmzZAocfHnUkIpIJsrbAvuiii7jooouiDkNERDLAP/4RdQQikkmydoiIiIiIiEgqqMAWEREREUkiFdgiIiIiIkmkAltEREREJImytsCeMmVK1CFkBP2cqk8/q+rTz0oOVbb+29Hnzi763JlPBbYclH5O1aefVfXpZyWHKlv/7ehzZxd97syXtQW2iIiIiEgqmLtHHUNSmFk+sK4Gb8kFtqUonPpEP6fq08+q+vSzKu0od28XdRB16RBydly2/tvR584u+tzprcqcXW8KbBERERGRdKAhIiIiIiIiSaQCW0REREQkibKuwDazZmb2RzP7XzObbGatoo4pHZnZBWa2xsw+NrM/mFnDqGNKd2bW2MwWmtmYqGNJd2Y20sx+aGbjzSw36ngk/WVr7s7mXJzNOTXbcqSZ9Tezu83se2b2ZzMbEnVMtZU1/6Mm+DPwvLs/ZmZfAO4BLo04prRiZt2A8cCFwDHAZOBDYFKUcWWAHwF5UQeR7szsG0BPd/9Z1LFIRsm63K1cnJ05NUtz5IPAee6+MfbvfgbQP+KYaiWrCmwz6wRcBnwntulF4Ekzy3P3tZEFln6OAr7p7oXAfDMbAJxK9iT1GjOzUcBm4NOoY0lnZnYK8FXgjKhjkcyRxbk7a3NxtubULM6R/YGWsed7gdYRxpIU2TZEZAywzd33AMSS1hrglCiDSjfuPiv2s4nbCKyPKp50Z2bNgQvd/d6oY8kAdwFLgT+a2fNmdmLUAUlGGEMW5u5szcVZnlOzNUc+AfzNzFoClwPXRhxPrWVbgd0Z+LjMtp1ApwhiySTDgb9EHUQauwH4ZdRBpDsz6wsMBf7q7t8FXgKeN7MO0UYmGUC5O8iWXJyVOTXLc+Q1wOfAXGCXuz8ZcTy1lm0FthO+ekjUBNgfQSwZwcx6AfnuvjDqWNKRmZ0FzHP3rVHHkgGOAT5x90Wx138CDDgvupAkQ2R97s6WXJzlOTWbc2Qz4BHgYeD3ZjY24nhqLavGYBO+XmtTZltLYFMEsaQ9M2sAfAv4cdSxpLEfAkPNLP66DfC0md3h7r+OLqy01JCEnOPue81sJXBEdCFJhsjq3J1luTibc2o258h/Al9x9+1m5sAjsXssdkcd2KHKth7sV4BOZtYMwMwaAd2AVyONKn39EJjk7vuiDiSNXQYMSWibgAlkx9e4NbUIaF3m685C4P2I4pHMke25O5tycTbn1KzMkbFpCAe7+/bYptuBHWT4LCJZVWC7+0fAs8DpsU1jgSfdfWN0UaUnM/spsABoZmY9zOybsa8oJYG757v7hngDDhC+xt0RdWzpxt2XAc8RphzDzA4HGgPTooxL0l825+5sy8XZnFOzOEd+AnxuZp0Ttn0MrIgonqTItiEiAN8Gfm1mg4GOwMSI40k7ZvYz4LYym5dl6R3dklxXAH8ws8aEHshLy8ySIFKZrMvdysVZKetypLsXmdl44GYzmw90AH6U6X9UmbtHHYOIiIiISL2RVUNERERERERSTQW2iIiIiEgSqcAWEREREUkiFdgiIiIiIkmkAltEREREJIlUYIuIiIiIJJEKbJFKmNnpZjbHzMZUcdwAM/uXmd2UxGv/JlnnEhGJgpkNMrPHk5kbK7nOs2b2lVReI1XMbKyZvVHV7xnJPNm40IxkGDPrD9xBWO1pODAI+Jq7/zPFl/4QGFyN49YDbYAGybhobBno/ck4l4gIgJkNB+4ETgb+H7CGsGBPI2Cyu89IwWVXAa1JUm48iD8D76XixLF8fAFhqfYVwBuxXU2ALwMj3H1tLS5R3d8zkmFUYEsm+A9wsbsvNLMc4B91cVF3X2FmH1fjuB1mtjmJlz4FeC2J5xORLOfuc83sUaCfu38xvt3MzgX+bWY/dPc/J/mae5KcGwEws3HA4nhh6+7Tk32NOHffDzxmZtcCM939ZwlxPH0o5zSz9sCp7v6Yuy83s21JClfSiIaISFozs7ZAH2A3hCVVgZuBulqCtCjJx1XHqcArSTyfiAhAuSW33X0a8HPgd2bWMQXXTGquNrMuhB7rulbRcuUvAhtrchIzawI8ADRL2JzM3x+SJlRgS1pz90+A5cDzZjYytm0d8LaZHW5mt5nZ22b2JTP7yMyWmtmx8feb2aVmdqOZzTSz/03YPszMfm5mj5nZU2bWOra9qZndbWY3m9nfCV9v1oiZ/U8spovNbLiZ/Sd2vttiMT5mZu1i2z8zs1vKnKKxu++LnesHZvZ9M3vdzH5Z4x+giEjV/k4YKnIRgJn1NLNfmNnfzOylWFGLmTWJ5bIbzWyWmV0SP4GZXWNmN5nZm2b23couVFFONrMLYznz62a20sx+aWadzOwfsXw6y8xOi53iMqALcIOZjY/dK/OimV2ZcI2rzewOM3vAzO43s+Zm1trMbjGzeWZ2mpktN7PVZpZ3KD+w2M9kTKyHGzO7wMz+YGa/NbPnzaxF7DP8wcymxn7fbI79jI8GLjKz7yecsk3sd9FOM7u6kmvmmNmvzGxi7BpuZv+O//eRNOPuampp3YBjgHWEv/LvBdrHtjcgJKtPCGPhDgdeAJYQ/ng8AfhZ7Ng2hF7wc4FWwKMJ538d+E3s+e+AibHnHYADhCRaVYz3EXqB2hPG6h0W294EmAE8DXQi/GLYD9wCHAYcF3sdPz4+vjz+uf8ce94S+HHU/y3U1NQytwFfBzZUsm8L8H+xvPoU0CC2/WHg8djzu4ALYs+/CeTHnl8KXB573ieWqwfGXt8H/Dz2vLKc3BbYDvweGAGMBSYBN8eO/REwNSFWB/Jiz5sBS4Gvx16fBTyccOyzwD2EIbHnx35fXBD7nM8Dd1TzZzcz9rvil8BvCOOxxybs3wZ0iz1fQvid1AS4EdgAHA9cFcv7M+Pxxo5fC/w2tu9SYH0lMUwAno097xD7OZ8Z9b8rtYqbxmBL2nP3JWY2ELgduAYYZ2ZjPYzJzgd2uPuTAGb2c2A20IPwy6TIzP47dqoZhER+LnBEwvalsfe2AL4dey/uvsXMPqxBqLnAH4Gr3H1P7Bz7zOwjYI27b4pd5yPgVQ/jExcQEn8u4WbJcwm/DAD2AFeY2QpgMuGPCxGRVHBC0TkCOAq41swAPgZaW3jxLULnAMD9hCESEHLtvDK59khgcZlrfJ0KcrK7f2JmnwFPu/vbAGa2HthtZkcSenxbVBi0+14z25qw6bux88b9DfhXbPt2wu+Lp2LXmBeLs7pe8dgYbDN7jNCpE3eRu683sxMJhXWLWP7fBHzg7m8Bb8XeW9G5p8Z+J8w9SEzHAntjn3uLmX1C6aEmkkZUYEtaM7MGQFcPN7NcFxu2MR2YQugRKGtJ7PFwoBvwoLs/HNv2+9g5byAUvL8vc63jCIlx5yGG2wc4CegJLEzYXnYM4oHiHe5FsWQbH651hIdhMbj7GjP7OvAn4HvAV4GtiIgkkZkdTvj2bQkhb35cQX7sQOhhbQzg7gcInQLE3nOru8+JvS713gQV5uQYp3Su3ADcAKwE3gS6H+QjJL6vN/BywusPCMNf2lM+Fxdy6DOcvEP4QyRuu5n9CnicEHu8ii77uapSSOW12XPAZDNrSvhvsZfQGy5pSGOwJd0VjwsEcPeFwE+BfpUc35iQzFYBm4ELE3ea2fGx7eeaWbMy23fEXvY/xFjnAH8AHoolwBoxs04k3DBjZp0JM6j0Idz0+HAlbxURqY0rCd+YPU7IjyfGeo6B4vy4DSgAvpCwfaCZtaRMrjWzhrEOi7Iqy8kV+Quw1sN0rPtq8FnWE3q8iy8BfEYYApM07l4U6wQ5JvYHynTgF+4+L5nXKXPNqcBUwjetlwKnuPtnqbqe1I4KbMkE18aKz7jOhLFzcW3NrHHs+RjC+LvtwCPA+Wb2JzMbE+td2E1IhM2BZ83srNgNJf0IRfkC4DYza2xmuYSxz+3NrIGZ/djMLq8kRou1mwm/hO6qYF/Z48u+/iJhvGBcb2C8u+8Arq3gPSIiNVGuZ9TMTgZuBb7l7vmE3uItwHQzGxe7kXF0rMf6UcJsI+PN7HTgSnffSci115nZrWZ2CqGjIT68LjH/VZaT4xJ7k4cC7czsMOBEoJmZ9Yjt20+4KTDe0ZJ4jb8AXzazNrHXxwN/dfdCQs1TYS42s5GxuJpU8rPLoUzNFPsj5FzCt5ZtCEMP+wJdgcMSbqAs20teUEH8pYOqYBxJ7L9FK+AlQi99Tf7wkDqmISKSCboCS8xsGuHrMwP+K2F/IXCjmW0n9PZeC+DuL5jZdYSvGccTbhJ8D8DMzgPuJvzCuBf4kbu7mV1MmEJpCfAY4au+Tpau3QAAIABJREFUIYSFGS4mDD0ptcCNmY0g3LzTjfCV4X+AW2PxzABGAnlm1pPQs9IRuNjMllHSm/N1oLm7Tynz2f9hYYGIQsJNRSIiNRbLI5cRisC/AvmEG+VaA2e7+xsA7l4Qy49/JeS6pwljrwGuJ3ROPED4xi4+c8ffCbnv24ShbN92962xAnIk0NPM/llZTjazCwnjjieY2TJ330i44XISobh+ADgPGAisBh4k9LZfYWF2qWOAc8xsurs/ZWa9gSfN7HXCt6A3xnraLwI6xj7fEsKiO0ea2TGxOL9CGFf+n4SfWyNC7h8MdI4V4E4Y+3xB7JzvAq8CcwnDXp4j5PRnCTdWDjaz8e4eP+9jhEV/dsUK8o7AhWa2htAzTez9Zdd8WE+4iXM80BTIMbMXCTc61tXUtVJNpv8mksksLC97n7vn1dH1fuTuWsZcRKSeiQ1XaeLuabnQl5n9iDCjy9rY66aEbx9+Fb93R9KHhohIpqto+EVqLhS+Sn2nLq4lIiJ1J3YT57FpXFw3BX4G5CVszgV2q7hOTxoiIhkrNv7tcsJXfhe5++MpvuQHsa8uRUSkfmno7pOjDqIy7v65heXa/y92A/wmwnCWnx38nRIVDREREREREUkiDREREREREUkiFdgiIiIiIkmkAltEREREJInqzU2Oubm5npeXV+P35efn065du+QHlGb0OeufbPms2fA5/z97dx4fV13vf/z1ydo0SZOmSdOWrtCVtmAh7CBFFhH8CRdRryiKXCzuyFUuXkW9Cnpd0HtxtypuV7xeQUUQUERk31qwLUtpS/c9TZsmaZOmaT6/P74zzWSaPbMl834+HucxM2fOnPM5k+Qzn3zn+/2epUuX7nL34X2ScZSzk0PvT8/0/vRO71HPamtr2bhxY685e9gU2FOnTmXJkv5fobSmpmZArxtqdJ7DT7acazacp5ltSHcMqaacnRx6f3qm96d3eo96VlNTw8aNG3vN2eoiIiIiIiKSQCqwRUREREQSKOsL7EWLFqU7hJTQeQ4/2XKu2XKe0jf6feiZ3p+e6f3pnd6jnvX1/Un5hWbM7HzgFuAd7r4+su4y4BvAKOB/gE+4e1vkuUrgS0A9kA/c6O4H4/dbU1Pj6jMkIkORmS1195p0x9EV5WwRkc76krNT2oJtZtVACXByzLrJwKXA5cD1wDXAx2Ne9lvgdne/EVgFfDFlAYuIZDHlbBGRgUlpge3uO4C741ZPAa5296Xu/gvgu8A5AGZ2KjDL3Z+JbHsP8FEzK05yoGEREcliQyZni4hkmJT3wXb39rjHj0W/WozYAmyM3H8DsC5m2y2Ru8n7KvUDH4BJk2Djxt63FREZ5jIlZzc808DSU5ey8uqVg92ViEjSZeIgx5OAH0TuHwXUxT3fCEyIf1FtbS01NTWHl8WLFw/s6Js3w5Yt8MQTA3u9iEgfLF68+HC+AirTHc8gpCZn50LjM400Lm1MQMgiIv3T35ydUReaMbPpQK27L4uscqA5brNC4IgBM1VVVYmZGP2MM+BPfwoF9hVXDH5/IiJdWLRo0eHR6Ga2K83hDEgqc/aIKSMAaNnQMqBYRUQGo785O2NasM0sF7gW+FTM6i3A6JhtDCgGtiYtkDPPDLePP560Q4iIDHWpztn5lfnkFOVwaO8h2va29f4CEZE0ypgCG/gkcKu7H4hZ9yAwK+bxVKAVSN7cTjU1kJ8PK1bA3r1JO4yIyBCX0pxtZmrFFpEhI+UFdqRFA8Bi1n0aeAEoMrOjzexqM5vu7kuAbWY2L7LphcC33b01aQEWFYUi2x0eeyxphxERGQoyKWcXTikEVGCLSOZLaR9sMysF3h15eJWZfYfwFePNcZuudPfbI/ffAXzOzNYT5mP9TNIDveACeOopuP9+ePObk344EZFMlGk5O9qCfWDjgV62FBFJr5QW2O7eCHw/skTdElm6e80G4F+SHFpnF10EX/gC3HdfaMk+3IAjIpI9Mi1nq4uIiAwVmdQHO3PU1EBVFaxfD6+8ku5oREQEGDEtFNjNr8VPVCIikllUYHclJye0YgPcdVd6YxEREQBGzhoJwP6V+9MciYhIz1Rgd+ed7wy3v/ylLpsuIpIBogV28+pm2tvae9laRCR9VGB359xzoboaVq+G555LdzQiIlkvtziXwsmF+EGnZa36YYtI5lKB3Z28vI4rOf74x+mNRUREABg5J9JN5BV1ExGRzKUCuyeRS2Lyy19CbW16YxEREUbOVj9sEcl8KrB7Mns2XHwxtLTAd7+b7mhERLJe8dxiAJqWNaU5EhGR7qnA7s0NN4Tb//ovqKtLbywiIlmutKYUgMaljWmORESkeyqwe3P22eHKjg0N8OUvpzsaEZGsVjy3GCswmlc107a3Ld3hiIh0SQV2X3zlK+H2O9+BNWvSG4uISBbLKcih5PgSABqfVyu2iGQmFdh9sWABvPe90NoK116rebFFRNLocDeR51Rgi0hmUoHdV7feCpWV8Le/wU9/mu5oRESy1qhTRwGw97G9aY5ERKRrKrD7qrIyDHQEuO46ePXV9MYjIpKlys8pB6D+kXpd0VFEMpIK7P5417vgHe+ApiZ429tgv+ZhFRFJtRGTRlA0vYhDjYdoel7T9YlI5lGB3R9msHgxzJgBK1bAlVfCoUPpjkpEJOuUvyG0Yu95aE+aIxEROZIK7P4aNQp+/3soK4Pf/S50F9GgRxGRlBp93mgAdt+3O82RiIgcSQX2QMydC3ffDQUF4QqPN96oIltEJIUqLqzACo29T+yldUdrusMREelEBfZAnX02/PrXkJcHX/86XH+9imwRkRTJK80LrdgOu+7ele5wREQ6UYE9GJddBnfeCfn5cNtt8C//EubKFhGRpKu6rAqA2t/WpjkSEZHOVGAP1iWXhO4iRUVhfuwLL4Q9GnQjIpJslZdWYoXGnof20LKhJd3hiIgcpgI7Ed70JnjkERg3Dh5+GE47DV5+Od1RiYgMa/kV+VT9UxU4bPvptnSHIyJymArsRDnpJHjmGZg/P1yEpqYGbr9d/bJFRJJo/DXjAdh++3ZddEZEMkbKC2wzO9/MnjGzqTHriszs22b2ZTP7vpmN6stzGWfyZHjySXjPe6C5OfTJfuc7oVb9A0VkaMr0nF1+TjlF04s4sOkAu+7SYEcRyQwpLbDNrBooAU6Oe+p7wOPu/mngbuCHfXwu85SUwM9/Dj/7GYwcCb/5DcyZA3fcodZsERlShkLOthxj0icnAbDxKxtx5VkRyQApLbDdfQch4R5mZhOAK4B7Iqv+CrzFzKb29FxKAh6M9743XO3x3HOhri5cZv3Nb4Y1a9IdmYhInwyVnF393mryq/Np+kcTu/+sC8+ISPqlvIuIu8d3klsI7HL3/ZHn24B1wNm9PJf5jj4aHnwQfvKTcOXH++6DY48NF6ZpbEx3dCIivRoKOTt3RC6TPhFasdfeuBY/pFZsEUmvTBjkeBRQF7euEZjQy3Od1NbWUlNTc3hZvHhxUoLtNzO4+mpYuRKuugoOHoSvfQ1mzgzdSA4dSneEIpIGixcvPpyvgMp0x9MPGZmzj/rIURROLmTf8n1s/8X2Qe1LRCRef3O2paO/mpk5MM3d15vZJ4G3ufspMc8/D9wRedjlc+5+a+w+a2pqfMmSJSmIfpCefRY+9rEw4wiEFu0vfjFctMYsvbGJSFqY2VJ3r0l3HN0ZKjl7x6928Mq7X6FgfAEnvXwS+eX5Cd2/iAj0LWdnQgv2FmB03LpSYGsvzw1NJ58cZhr5xS9gypQwX/bll4dp/e67TwMhRSTTZWzOHvvOsYw6bRSt21pZ+29rU3FIEZEuZUKB/TAwwcyKAMwsH5gMPNLLc0NXTg5ceSWsWgXf/S6MHw/PPw8XXwwLFsCvfhW6koiIZJ6MzdmWY8z68SyswNj2o23s+Zuuqisi6ZGOebCj/SAMwN23A/cC50bWnwfc5e5benouhSEnT0EBfOhDYWaRr38dqqth2TJ497th+nS47TYNhhSRtBpqObv42GKmfHYKAK+85xVad7Wm6tAiIoeleh7sUuADkYdXmVlV5P4HgUvM7DPARcCimJf19NzwMHIkfPKTsH49/OhHYQDkxo3w8Y/DUUfBhz8ML76Y7ihFJMsM1Zw9+cbJjDp9FK1bWll55Uq8XV3vRCS10jLIMRmGzCDHvmhvh3vugVtvhccf71h/5pmhxfuyy6CwMH3xiUhCZfogx2RIds5u2dTCkgVLaKtrY8rnpzDtP6Yl7Vgikl2GyiBHiZeTA5dcAo89BsuXh6K6pCQU21dcEfpsf+hD8PTTGhQpItKFEZNGMOeXcyAHNnxhg6buE5GUUoGd6ebPDwMht26F738fXvc62LMn3D/tNJg9G770JdiwId2RiohklDFvGsP026YD8Oo1r7LnIQ16FJHUUIE9VJSWwgc+AC+8EAZCfuITMG5cmInkpptg6tTQheS222DL8BgDKiIyWBM/MpGJ10/EDzor3rKC+kfq0x2SiGQBFdhD0XHHhf7ZmzaFubP/+Z9hxAh44okwMHLiRDjjjFBsb96c7mhFRNLqmFuPYdxV42jf387yi5ZT/6iKbBFJLhXYQ1leHrzpTfDrX8POneH2sstCsf3kk6HYnjQpdCX58pdhxQr12RaRrBOdH7v6PdWhyL5wObvu2ZXusERkGFOBPVyUloaW7LvugtrazsX200/DZz4TWr6nTg3T/j3wALS0pDtqEZGUsFxj9u2zGXf1ONqb23nx0hfZ+sOhe1FgEclsKrCHo5KSzsX2738PV18dLmSzcSN873uh5XvMGLj00jBg8rXX0h21iEhSWW5oyZ7y+SnQDqs+sIo1n1hDe1t7ukMTkWFG82Bnk/Z2WLIE7r03zLP9j390fv7oo+GCC+CNb4RzzoGysvTEKZJlNA926m398VZWf3A13uaUnV3Gsf97LIXjdH0BEeldX3K2CuxstmkT3H8/PPgg/PWvUB8z8Cc3F049NRTcF1wANTWhz7eIJJwK7PSof7yel9/2Mq3bWykYX8CcO+YweuHotMYkIplPBbb03aFDoXX7L38Jy1NPhXVR5eWwcGFo2V64EObNCxfEEZFBU4GdPge2H+Dld7zM3kf3gsHE6ycy7UvTyB2Rm+7QRCRDqcCWgWtogIcf7ii416zp/HxlJZx9dkfRfeyxYJaWUEWGOhXY6dV+sJ0NX9rAhls2wCEYeexI5vxiDqUnlqY7NBHJQCqwJXHWrQsFd3SJv5hNVVUothcuDBe8mTs3dDMRkV6pwM4MDc828MqVr9C8qhlyYOLHJjL1i1PJK1X3OBHpoAJbksM9zDry8MPw97+H223bOm8zalTow33GGXD66XDKKWEqQRE5ggrszHFo/yHWfWYdm7+1GdqhYEIB02+bTtVbqzB9SyciqMCWVHGH1atDof3II+EiNxs2dN4mJweOPz4U22ecES5+M2WKupWIoAI7EzW+0MiqD66i8ZlGAMrPKeeYrx+jbiMiogJb0mjLllBoP/FEuH3hBWhr67zN2LFw8slw0kkdt2PGpCdekTRSgZ2ZvN3Z9qNtrP33tbTtCflr7BVjmfalaRRNLUpzdCKSLiqwJXPs3w/PPddRdD/1FOzefeR2Rx/dUXCffDIsWADFxamPVySFVGBntoN7DrLxyxvZ/K3NeKtjBcaED05g8r9NpnCC5s4WyTYpK7DN7CJ3v2/QOxqEoZSshdCtZN06ePbZsDz3HCxdCs3NnbfLyQlTAp50EpxwQii4jztORbcMK6kusJWzB6Z5fTPrPrOOnXfsBMAKjQnvn8CkGycxYuKINEcnIqmS8ALbzC4EbgLG0XGZ9Ryg2t3T+n3ZUEzWEqetDV5+uXPRvWJF5/m4IRTdM2eGYnvBgo7Cu6IiPXGLDFKyCmzl7ORoWtbEhls2UHtnLQBWYIx73zgmfWISI2eMTHN0IpJsySiwa4H/BF4A2qOrgSvcfdFAA02EoZyspQf794dLuj/3XOjH/cILoQiP788NMHlyR9EdXSZO1EBKyXhJLLCVs5Oo6cUmNn5pIzt/sxMcMBjz5jFMvH4i5QvLNeuIyDCVjAL7GXc/pYv149x9+wBiTJjhkKylj1pa4KWXOgruF16AZctCMR5v9GiYP7/zMm9emEZQJEMkscBWzk6Bfa/sY9M3NrHjf3bgB8JnasnrSph4/USq3l6lq0KKDDPJKLBPAI529zvj1i9y98UDCzMxhlOylgE4dAhWrepcdP/jH1BX1/X2U6aEvtyxhffMmZCfn9q4RUhqga2cnUKtO1vZ+v2tbPneFg7uPAhAXkUe4947jvGLxlM8W2NHRIaDZBTYHwVuBmInAjXA3X1Q/6Kb2RzgI8AaYAaw2N3/YWZFwNeARmA0cKO7N8S/fjgmaxkk93ABnBUrOi8vvwwHDhy5fUEBzJ7dUXjPnQtz5sDUqaHft0iSJLHAVs5Og0Mth9j5651s+fYWml5oOry+7PVlTFg0gcq3VqpVW2QIS0aBvQe4FlgCREeeGXClu9880EAj+14CXOLuW8xsMvBnd59jZj8FHnD330QG7LzX3d8Z//rhnKwlwdraYM0aWL68c+G9dm3X2xcVhcJ7zhw49tiwzJkDxxyjFm9JiCQW2MrZadawpIFti7ex444dtO8L3eDzyvOoensV1e+upuyMMixHfbVFhpJkFNh3AO9z9wNx6wfdn8/M9gEnuvtKM6sClgE1wDpgtLvvN7M8YC8w193Xx74+W5K1JFFTU+jbHS28X3kltHZv3dr19vn5oVtJfOE9cyaM0JRd0ndJLLCVszNEW2MbO3+9k62Lt9K0tKNVe8TUEYx911iq312tLiQiQ0QyCuy3ACcAf4tdDVzm7tcNKMqOff8cOAZ4E3ANsBEoBL7u7kfFbPdiZN3PY1+fbclaUmjv3o5i++WXO+6vX9/19jk5oXU7WmzPmhWWmTPD1Ss1s4DESWKBrZydgZpebGLnr3ay41c7OLCp43+fkhNLGPv2sVS9tYqiY3SlSJFMlYwC+zlCX7t6wqREEJL1OHcfVJOdmRUDdwMTgW+4+4/M7AbCV5nHxWz3FPBHd//P2Ndnc7KWNNm3D159taPwjhbfa9ZAe3vXrykr6yi6Y29nzNDFc7JYEgts5ewM5u1O/aP17PjlDmrvrOVQQ8ec/yULSqh6axVVl1cxcpbm1hbJJH3J2Xn93OfXgfvdvTHuQBf1N7guFAG/Bo4C/tvM1hE+EOIu7UchcDD+xbW1tdTUdJzrokWLWLQordO8ynBXXBwucnPCCZ3XHzgQZjR55ZVw++qrHbd794Y5vZ977sj9TZzYubU7ejtlCuRqQNRws3jxYhYvPjyRR2WSDqOcncEsxxi9cDSjF45mxndmsPv+3dTeVUvdPXU0vdBE0wtNrLtpHcXziql8ayWVb6mkZEGJ5tcWSYP+5uz+tmCf5u5PdbH+ane/vc876nrfDwD/7O71ZvZZ4GPAx4HPu/vMmO1WR9bdEft6tYZIxnOH2trOBXf0/muvwcEjapCgoACmTw+t3NOnd74/caKK72EgiS3YytlD0KGWQ+x5cE8otu+uo62+48JaBRMKGHPxGMa8eQyjzx1NbrH+/kVSLRldRH7p7lfGrXsb8BN3H/CVO8ysEljh7uMjj40w9dN7gQeAKndvNrN8oIkwr+uW2H0oWcuQ1tYW+nTHFt/R2+4GWUIovo8+unPxHS3AJ02CvP5+SSXpkMQCWzl7iGtvbWfP3/aw6/e7qLu3jtatrYefs0Jj9BtGM+biMVRcXEHRVPXbFkmFZHQRqTCzz7r7zWY2EvgW8Dagi++7+2U30GJmR8Uk4TpgOXAvcG7k9jzgrvhELTLk5eV1FMcXX9z5ucZGWL069O2OXVavhu3bYeXKsMTLz4dp045s9Z4+PXQ70RSD2UA5e4jLKchhzIVjGHPhGNydpmVN1N1bR929dTQ+28ju+3ez+/7d8BEomlHE6PNHh+Wc0eSV6R9skXTpVws2gJmdCrwLeCOwHXgPsNfd9wwqELPjgQ8BS4Fq4FF3f8TMRhMuWrAeGAf8u7s3xb9erSGSlZqaQveS2KI7en9LDzVNXl64gE58q/f06WF9QUGqzkBIXgt2ZN/K2cNU645W6u6vo+6eOvb8dU+nQZLkwqhTRjH6/NFUnF9B6cml5OTrglkiiTDoLiJm9vpunjoZOAe4FWgHLnX36wcaaCIoWYvE2b+/c/EdW4Rv2tT963JyQgt3tOCOvZ02TcV3EiSqwFbOzl7tbe00PtvIngf3sPvB3TQ83dBxaSEgtzSX8oXllJ9dTtnZZZS8roScPBXcIgORiAL7aWAOYYqnbuYdw4Dx7l440EATQclapB+am2HduiNbvdesgQ0bwoDMrnRVfEfvq/gesAQW2MrZAkBbQxv1f68PBfdfdtO8qvPkLrmluZSdWUb52eWULyyn5IQStXCL9FEiCuyLgCXuvrOXA13u7ncOLMzEULIWSZADBzoX39ECfPVq2Lix+zm+c3I6up3Et36r+O5RAgts5WzpUsvGFur/Xk/9I2Fpea2l0/M5xTmUnREK7rLXl1F6Yim5RZqhRKQrCZ9FJJMpWYukQLT4ji26o7e9tXx3VXzPmKE+3yS3D3amUs5Or5bNLex9ZO/hgju+hdvyjZITSig7vYxRp4+i7PQyCiek9UsPkYyhAltEUmegxXdubke3kxkzYPbsjmXChKy4tLwKbEm3A9sOUP9IPXsf2cveJ/ay78V9Hdf+jCicUthRcJ9WRvHxxerHLVlJBbaIZIauiu/o/Z6K75KSzgV3dJk+HQqHT2uaCmzJNG1722h4toGGJxvY++ReGp5u6DxLCZAzModRJ4+i9OTSw7eFEwt1pUkZ9pIxD7aISP8VFnYUx/Fii+/o1S2jc3vv2gVLloQlVk5OuMDOnDlHFt8VFak5J5FhLK8sj4rzK6g4P/w9+SFn38v7OgrupxpoXt0c+nX/vf7w6/Kr80OxfVLp4dv8Cs25L9ln0C3YZrYA+ATwAuGCAusTEFe/qTVEZBjatatzwR1d1q7tfrBlVVUovOfP71jmzYNRA75wYdKlsgVbOVsSpXVnKw1PN9D4XCMNzzXQ+GwjbXvajthuxDEjGHVSR0t3yYISckdqAKUMXSnpImJmPwQucfdxZnYF8Ht3b+7tdYmmZC2SRQ4cCC3e8YX3ypWwb1/Xr5kypXPRPX8+zJqVEVe0THGBrZwtSeHuNL/WTONzjTQ+G4rupuebaG+O+2c4F4rnFFNyQgmlJ5RSckIJJceXkDdKX6rL0JCqAjsfGOvuW8wsByhx94ZB7XQAlKxFBHfYvBleeglWrOhYXn4ZWluP3L6gAF73OjjppI5l1qww8DKFUlxgK2dLyrS3tbP/pf00PNtwuOje9+K+ThfBiSqaURSK7gWRontBCQWV2T3DkGSmpBTYZvYmoMLdf2VmJwNN7v7yIOJMCCVrEelWW1to8Y4tulesCF1N4pWUQE0NnHUWnHMOnHoqFBUlNbwkXypdOVsyyqH9h9i3Yh+NzzfS9HwTjS80sm/FPrz1yHqkcHIhJQs6WrpLF5RSMKFAAyklrRJeYJvZl4DrgAfc/fLIuq8Ay939jsEEO1hK1iLSb/X1sHQpPPdcxxJ/GfnCwlBkn3suvOUtcNxxCZ86MFkFtnK2DBXtre3se3kfTc830fRCUyi+/9FE+/4jx1rkj82n5HWhW0nJ8SUUH1/MyFkjdSVKSZlkFNgPA/8EvNfdb4usOxv4pbtPHkywg6VkLSIJsWMHPP00/P3v8PDDsGxZ5+enTYNLL4Urrggt3QmQxAJbOVuGLD/k7F+9P7RyP99I0wtNND3fRFv9kQMprcAonlt8uOCOFt+awUSSIRnT9D3o7vVmFluVX97/0EREMlR1NVxySVgA6urg0Ufh/vvh7rvDlIL/9V9hOfFE+NCH4MorM2KwZBeUs2XIslyjeHYxxbOLqb6iGggDKVvWt9C0rIl9y/bRtKyJpmVNtKxtCQX4C02d9lE4sbBTwV1yfAlF04uwXHUxkeTqbwv2PwNHA0cBjwFvI7SOfNTdv5uUCPtIrSEiknSHDoXW7TvvhF/8AnbvDutnzoTvfhfOO29Au01iC7ZytmSFtoY29q3YR9M/mg4X3ftW7DtyBhPCBXKK5xV36mJScpxmMZG+S9Ygx1OBq4DJwDbg1+7+14EGmShK1iKSUi0t8Nvfws03hwGUOTnw4x/D+97X710leZCjcrZkJT/kNK9pPlxwR1u9D2w+0OX2I6aNCAX3cR3F94hpI7ActXZLZ6mapi8XWOjuDw1qR4OkZC0iaXHwINxyC3zxizBiBGzdCqNH92sXKZ6mTzlbstrB3Qc7FdxNy5rY91LXs5jkluRSPL+4U+FdPL+YvFK1dmezQffBNrM/AiW9HKcS2AWkNVmLiKRFfj7ccAN861thVpING/pdYCeKcrZI7/Ir8hl9zmhGn9Pxd9p+sJ39r+5n3/JIwR25bd3WSsNTDTQ81Xmq+BFHd9HaPVWt3dKht3/Blke2eRVoB95PuLxubLPDNKAiKdGJiGQyd3joIbjuulBcz5gBc+emMyLlbJEByMnPoWReCSXzSg4PqARorW09XGxHC+99L++jZW0LLWtb2PX7XYe3zS3JDQX3cR0zmRTPU2t3turtp/51QjeSegAzO8bdPxe/kZndnYzgREQy0rZtcNddoc91dBq/GTPgT39K92wiytkiCVRQVUDBuQWMPreL1u7oLCbLQ1eT1u2tNDzZQMOTca3dx4yg5LjOAyrV2j389Vhgu/veuFXVZpbr7ocvcmpmC4GTkxCbiEhmaGsLF6F56CH485/hiSdC6zXA2LGhBfv665N+xcfeKGeLJF+n1u53xbR272w9XGw3LQ8t3vtf3k/Lay20vBbX2l3a0dod7WpSPL+YvBK1dg8X/f1J3gksM7PHgUZgDvCITUeeAAAgAElEQVRG4IgWEhGRIWvnTnj22Y7lqaegIaZVqqAALrwQ3vEOeOtbw9UeM5NytkiKFIwtoOK8CirO6+iB1d4aWrtj+3U3LWvi4I6DNDzRQMMTMXnFoOiYok79uouPKw6t3bo0/JDTrwLb3R80swuADwHHAfXA5e6urxtFZOhpawtT7K1Y0bEsWwbr1x+57cyZ4XLp554b5rsuK0t5uP2lnC2SXjkFOZTML6Fkfuexx607Iq3dMUX3/lf207ymmeY1zez6XUxr96iOmUwOD6ycX0JucW6qT0f6od/fRbj7VuCm6GMzyzGzE919aSICiszZeiawBngc2Ad8jdD6Mhq40d0but+DiEicpqZQSK9a1bG89BK8/DIc6GJO3OLicBn0k0+GU04Jy8SJqY87AZSzRTJPQXUBFedXUHF+XGv3yv0d3Uwi/bt7bO2OXqUyMrByxBS1dmeKfhXYZnZ7F6urgEPApYMNxszeBxzj7rEfBj8FHnD335jZhcAPgXcO9lgiMoy4w549oeV53TpYu7ZzQb1tW/evnTIF5s/vvMyZA7lDv3VIOVtk6MgpyAl9so8rgXd3rG/d0dppMGXT8tC3+3Br912dW7tLXldCaU0ppSeGpWhGkQZUpkF/W7BnEOZOjZ2NfSLhEryDYmZnA+8Czo9ZNwG4AvhwZNVfgbvMbKq7rx/sMUVkCGlsDMVztIiOv23ooZG0sBCmTw/dPKLL7Nkwbx6MGpWiE0gL5WyRIa6guoCKCyqouCCutfuV/YcHU0a7mhzceZC9j+5l76Md451zS3MpOSGm6K4ppegYFd3J1t8C+xp3fzV2hZnNJlyGd7C+CTwJfNvMpgM3A1OAXe6+H8Dd28xsHXA2sD4BxxSRTLBvH2zeDJs2dX9bX9/zPkpKYNo0mDo13M6cGabOmzkTJk0aFi3SA6CcLTIM5RTkHO6TzZUd6w9sP0DT8000Lm2kcUkjjUsbad3Syt5H9rL3kZiie1QupSeEYrvkxJKOolvdSxKmv4McX+1i9QHgGuBTAw3CzGYBC4D3uftyM7sBeAC4BaiL27wRmBC/j9raWmpqOq5auWjRIhYtWjTQkEQkEaJdN7Zu7Vi6KqD37Ol9XyNGdBTP8bfTpkFFBQyhD4fFixezePHi6MPKZBxDOVskuxSOK6TwokLGXDTm8LoD2w/QtLTpcMHduLSR1q2t1P+9nvq/dzRc5JXnUXpyKaNOG8WoU0cx6pRR5I9O67z+GaW/OdvcvbdtOjYOLRGxL8gFqoFfuvv7+xFn/H4vAxa7e2XkcRFQC3wReKu7nxKz7fPAHe5+a+w+ampqfMmS2IuViUjSuMPu3aFg3rat+9tt27oeRBivsDAMIpw4MbQ2R29j748ZM6QK6P4ws6XuXtP7lv3er3K2iBzhwLYDh1u5o8V36/bWI7YbOXtkKLYjRXfx3GIsd3jm4f7oS87ubxeRB4E76EjY7cAOd181gPji4zgci7s3m9lqwofB6LhtS4GtgzyeiHSlrS3MAb1jB2zfHpZooRxtgY4+bj0yGXeprAzGj4cJE8JtbNEcva2sHLbFc5opZ4vIEQrHF1L45kIq39zRENuyuYWGpxvC8lQDjUsb2b9yP/tX7mf7z7YD4XLw0Vbu8teXM+r0Ubo4Tjf6+678j7s/GrvCzGaZ2XR3XzOIOJYDZWZW7e47IuvaCEl5gpkVRRJ4PjAZeGQQxxLJLocOwa5dHUVzbPEcv66uruMKhb0pK+somru7HT8+THkn6aKcLSJ9MmLiCEZcPoKxl48FwkDKpn80HS64G55uoGV9C/V/q6f+b/Vs/NJGLM8oObGE8rPLKX99OWVnlpFXpoIb+l9gXwg8GrduB/AwoT/egLj7SjO7H7gc+K6ZlQMFwK+ANwHnAvcC5wF3ufuWgR5LZFhobw9dNLoqkuNva2vD9n1hFi79XV0N48aF2+6K55Ejk3uOkgjK2SIyIDkFOYw6eRSjTh4FHwvrDmw/QMPTDex9PMxU0vh8I43PhGXT1zZBDpS8LlJwn1NO+cJy8kqzs+Du01mb2UcIg1dKzezGLjZ5MAGxvAe4zcwKCC0e74yMQP8g8DUzOx4YB2gUjAxP7mGmjJ6K5ej9nTtDd46+GjOmo2Du6bayEvKyMxkOJ8rZIpIMheMKqbq0iqpLqwBoa2xj7xOh2K5/pJ7G5xpper6Jpueb2Pxfm7F8Y9Tpo8I0g2+soGRBSdZMD9jnQY5mdhzwXuBbcU+1unsPV3FIDQ2YkYzkHq4i2F2XjPh1fe3XDFBe3reieexYyNdI8EyWjEGOytkikmqH9h+i4akG6h+pZ89f99DwTEMY+RGRX5nP6PNHU/HGCiourqCgsiB9wQ5CX3J2f2cRyXf3g12sn+zuGwcQY8IoWUtK7d/ftz7N27dDc3Pf91ta2nWRHL+uujrMvCHDQhJnEVHOFpG0OVh/kPqH6tn9l93s/vNuDmyImVUqB8rOKqPy0koqL6mkaFpR+gLtp0HPImJmnwTuc/eXI6v+zczir9ZgwDnAwoEGKpIRDhwIXS/ii+WuCufGxr7vt6ioo0juqXiurla/ZhkU5WwRyST55flUvbWKqrdW4e40r2pm9192U3dvHfUP1x++AM5r179G8fHFVF5aSfU7qxk5a+h/FvbW2XIaUB7z+ASgGNhOx7RPRuhnJ5J5upp2rruW575c7CSqoKDnFubYdSUlmoJOUkU5W0QykpkxctZIRs4aycSPTqRtbxt199ex6w+72H3fbvYt28e+ZfvY8IUNlJ5USvW7qhn7z2MpqB6i3Uj62UXkBOBFd2+NW39sTItJWujrxizT2toxR3PsVQJjlx07wvR0ff0dz8sL/ZX7UjiXlaloloRJYhcR5WwRyXjtB9rZ87c91P5fLbV31XKo8VB4Ihcqzq9g/LXjGfPmMeTk5aQ30IhkXGhmvrs/H3eQWcAxQFqTtQwT0Rbn7orm6FJb27f9mXV0v+itpbmiAnIy449XJEGUs0Uk4+UU5jDmTWMY86YxzPjeDOruqWPH/+xg9/272f1AWAonFjL+2vGMv2Y8heMyfwxSX6fpO4vwteL5kUvvxqogjFK/J8GxyXDiHrpgbNoEW7b03Orclzmbc3JCYTxhQtdL9CInlZWQG98FVWR4U84WkaEqtyiXsW8fy9i3j6V1Vys7frmDrd/fSvPqZtZ/dj0bvrCB6ndXM/nfJzNyZub21e5rC3Y98AtC/74z4p5rBX6cyKBkCGpqCsXzpk2wcWPH/dhl//7e9xNtce6qYI59PHasCmeR7ilni8iQV1BZwKTrJzHxuons+dsetn5vK7vu3sX2n21n+8+3U/X2KqbcNIWSeSXpDvUIfSqw3X2FmS0E3ubuSszZ5tAh2LwZ1q3rvniur+99P6WlMGkSTJwIRx3VdctzdbXmbBYZJOVsERlOLMeoOK+CivMqaH6tmY1f3cj2n22n9je11P62lvH/Mp5pN0/LqAGR/Rrk2O1OzM5098cTEM+AacDMILiHwYDr1sHateE2uqxdG4rq3q4aOGJEKJwnTYLJk8Nt/FJWlprzERlikjXIsYfjKWeLyJDWsqmFjV/dyLYfbsPbnNySXKZ8fgqTrp+E5SZ3EoJEzIP9R6C3dncD5gJj+xeepFR7e2hpXrUqLKtXdy6o9+3r+fXjx8O0aTBlStfFc2WlZtUQSTPlbBHJFiMmjWDmd2Yy8aMTee2Tr1F3bx1rb1jLrrt2MeunsyieXZzW+HrrIvIqoS/fFjrmUI1nwJsTGZQMkHuYXSO2iI69f+BA968tK4Ojjw5FdOxy9NGhqC4aOldYEsliytkiklVGzhrJ/HvmU3dfHa8uepWGpxtYumApM38wk3HvTd+U/70V2F8FGt29h8oMzEyj0VPJPcy48dJL8OKLYXnpJXj1Vdi7t/vXjR8PM2eGZfp0OOaYjkJ69OjUxS8iyaKcLSJZacxFYzjpxZNYc90advxiByuvWknj0kaO+eYxaZk/u8cC29139XE//4RGpSdHXR0sX95RTEdvuxtUWFbWUUTHLjNmhEGGIjJsKWeLSDbLL89nzs/nUHZmGas/vJot395C67ZW5twxh5z81BbZ/b3QTCdmNhn4MPARlKwHr7YWnn8eli7tWDZs6HrbigqYNy8sc+eG5dhj1RdaRLqlnC0i2WDC+ycwcs5IVly8gto7a/E2Z+6dc5M++DFWvwtsM8sF3gK8Hzgf2EHo8yf9ceAALFkCTzwBTz0ViulNm47crqgI5s8PS2xBPW6cCmkR6ZVytohko/Izyzn+oeNZfsFydv1hF6/d8BrTvzk9Zcfvc4FtZtOAa4D3AWOADcD/A/4MnJaU6IaT+np47LFQUD/+eCiu4wcdFhfDggVwwglw4olhmTUL8gb1RYOIZCHlbBHJdqNqRjHvD/NYdt4yNv/XZkprSqm+ojolx+61cjOztxFaPt4ArAG+Cfwc+LS7PxDZ7MmkRThUtbXBc8/Bn/8Mf/kLPPPMkZcAnzcPzjgjLCedFPpJ6+qEIjIIytkiIh3KX1/OjG/PYNUHVrH6I6spP6ecwvGFST9uX5pGmwmX1n0SuMbdVwGY2eCvUDPcNDXBn/4Ev/tdKKpjByLm5cHpp8NZZ8GZZ8Jpp2nmDhFJBuVsEZEY4xeNZ9fdu9h9/25eu+E1jv2fY5N+zF4LbHe/F7g3Mjjm/WY2FbgbOHw9azMzT8QlIYeiffvgD3+AO++EBx6AlpaO52bMgDe+ES64ABYu1CweIpJ0ytkiIp2ZGTO+N4NnZz7Lzjt2MvnGyZTM7+2aXIPT58697r4R+KyZ5QGXAhPN7GbgJ8AFwOLkhJihnn8efvQjuOMOaGjoWH/GGXD55XDJJWF+aRGRNFDOFhHpUDS1iAnXTmDLd7aw6dZNzPn5nKQer9+TArp7m7vf6e7/ROjXdz3w3wmPLBO1t8Pvfw+nnBIGIP7gB6G4PuUU+M53YMuWMIDx4x9XcS0iGSGrc7aISIyJH58IQO3/1XJwz8GkHmtQ01O4+xrgOjNbmqB4MpM7/O//ws03wyuvhHWjR8OVV8I114Qp9EREMlzW5GwRkS4UHVNE+bnl1D9Uz67f72L81eOTdqyEXNbG3X+RiP2YWYGZLTOzhZHHRWb2bTP7spl938xGJeI4/fLss2Fw4hVXhOJ68mT41rdg82a47TYV1yIy5CQqZ0OG5m0RkW5UXVYFwO77dyf1OKm/OHvPbgCmxjz+HvC4u3+aMEjnhymLpLUVbrwRTj0Vnn4aqqtDn+s1a+CjH4WRI1MWiohIBsucvC0i0ouKiyoA2P3gbrw9eWO9M6bANrPTgW3AnsjjCcAVwD2RTf4KvCUyIj65Nm8O0+h97Wvhaok33ACrVoXuIPn5vb9eRCQLZFTeFhHpg6KpRRROLOTQ3kM0r25O2nEyosA2s2Lgcne/PWb1QmCXu++HMFAHWAecndRgVqwIrdbPPx8GKj72WCi0R+lbThGRqIzK2yIi/VByQpiir/H5xqQdIyMKbOBG4Ctx644C6uLWNQITutpBbW0tNTU1h5fFiwcwA9WaNXDeeWE2kLPOCpczP/30/u9HRKQXixcvPpyvgMp0xzMAg8rbCcnZIiIDUHJ8KLD3vbivz6/pb84e1CwiiWBmbwSWuPvOuKeccEWyWIVAl/OqVFVVsWTJkoEH0tAAF14IO3fC+efDH/8II0YMfH8iIj1YtGgRixYtAsDMdqU5nH5JRN4edM4WERmgEdNCfdeyoaWXLTv0N2envcAGPgksMLPo49GEgTHfjNyPVQpsTUoUH/kIvPYaHH883HWXimsRke5lRt4WERmAEVNDjXdgw4GkHSMTCuwrCC0cUU8B/wo8A9xgZkXu3mxm+cBk4JGER/DII/DLX0JRUZjvWpc0FxHpSfrztojIAI2YHGnB3tT3Fuz+SnuB7e61sY/N7BBQ6+4bzexe4FzgXuA84C5335LgAOBTnwr3P/UpmD07obsXERlu0p63RUQGIa8ilL9te9qSd4yk7TkxPgh8zcyOB8YBixJ+hCVLwjzXY8bAv/5rwncvIpJlkp+3RUQGIa8sDwwONRzCDzmWa72/qL/HSPgeB8ndp8bc3wO8P6kH/MlPwu1VV0FJSVIPJSIyHKU8b4uIDILlGHllebTVt9FW30b+mMRf4yRTpulLD3e4995w/13vSm8sIiIiIpISeeWRbiL1yekmkt0F9iuvhDmvq6vD7CEiIiIiMuzlFIUSuL2lPTn7T8peh4roHKxnnQU52f1WiIiIiGSLnMJIgX1ABXbiLV8ebtV6LSIiIpI1rDAMbFSBnQwrV4bbuXPTG4eIiIiIpEy0BdtbPTn7T8peh4rNm8Pt5MnpjUNEREREUkZdRJJpS+TaBxMmpDcOEREREUkZFdjJcvAg7NoVBjeOHZvuaEREREQkRSw/9MH2g+oiklhNTeG2tBRyc9Mbi4iIiIikTrQCTk59ncUFdmNjuC0tTW8cIiIiIpJSlhNpwW5XC3ZiqcAWERERyU7RCjg5XbBVYKvAFhEREckuasFOlpaWcFtUlN44RERERCS1LHKrFuwEO3gw3ObnpzcOEREREUkptWAniwpsERERkeykPthJEi2w8/LSG4eIiIiIpJRasJOlrS3cqgVbREREJLuoBTtJ1EVEREREJCupBTtZ1EVEREREJDtFZxHRlRwTTC3YIiIiItkpUgGrBTvRDh0Kt7m56Y1DRERERFIq2kVk2PfBNrPLzGydmdWZ2W1mlhdZX2lmPzSzr5rZN80sMU3O7tEDJ2R3IiLZJOU5W0QkkSLl37BuwTazycClwOXA9cA1wMcjT/8WuN3dbwRWAV9MyEFVYIuIDEhacraISCIluQ92pozwmwJc7e5twFIzmwecY2aPA7Pc/ZnIdvcAr5rZLe6+b1BHVIEtIjJQqc/ZIiIJZEmu/zKiBdvdH4sk6qgtwEbgDcC6mO22RO7WJOCg4VYFtohIv6QlZ4uIJEOWzSJyEvAD4CigLu65RmDCoI+gAltEJFGSn7NFRBIpyeVfpnQROczMpgO17r7MzBxojtukEDgY/7ra2lpqajoaSRYtWsSiRYu6P1B7ZNhoTqb+jyEiw9nixYtZvHhx9GFlOmMZjJTlbBGRZOhjC3Z/c3ZGFdhmlgtcC3wqsmoLMDPmeQOKga3xr62qqmLJkiV9P5hasEUkjWILSjPbleZwBiSlOVtEJJH6Wf71N2dnWvPtJ4Fb3f1A5PGDwKyY56cCrcDgs7IKbBGRwUpdzhYRSQL3YTxNH4CZfRp4ASgys6PN7GqgHtgWGaEOcCHwbXdvHfQBVWCLiAxYynO2iEgiZUMfbDO7Cbg5bvVKd7/dzN4BfM7M1gMlwGcSclAV2CIiA5KWnC0ikgzDeR5sd78FuKWb5zYA/5KEg4ZbFdgiIv2SlpwtIpJAWTEPdlqowBYRERHJblk2D3byqcAWERERyU5JvlS6CmwV2CIiIiLZJcnlnwpsFdgiIiIiWWnYT9OXciqwRURERLKTWrCTRAW2iIiISHZTH+wkUYEtIiIiklU0TZ+IiIiISDKoBVtEREREJAHUB1tEREREJAnUgp1gJSUwfjyUlqY7EhERERFJobyyPArGF5BTnJxSOC8pex0KPvzhsIiIiIhIVpnymSlM+cyUpO0/e1uwRURERESSQAW2iIiIiEgCqcAWEREREUkgFdgiIiIiIgmU9QX24sWL0x1CSug8h59sOddsOU/pG/0+9EzvT8/0/vRO71HP+vr+qMDOkl8knefwky3nmi3nKX2j34ee6f3pmd6f3uk96pkKbBERERGRNDD3JF3CJsXMrBbYMICXVgK7EhxOJtJ5Dj/Zcq7ZcJ5T3L0q3UGkknJ20uj96Znen97pPepZJVDcW84eNgW2iIiIiEgmUBcREREREZEEUoEtIiIiIpJAKrBFRERERBIoawtsMysys2+b2ZfN7PtmNirdMSWSmZ1vZs+Y2dSYdcPqnM3sMjNbZ2Z1ZnabmeVF1lea2Q/N7Ktm9k0zy093rINlZmeb2atmtsfM/jtm/bD6mUaZWYGZLTOzhZHHw/I8pe/0O9C7rvK+BN19XkiH7j5npLP4z6fuZG2BDXwPeNzdPw3cDfwwzfEkjJlVAyXAyXFPDZtzNrPJwKXA5cD1wDXAxyNP/xa43d1vBFYBX0xLkAliZmXAGcCpwLuBD5vZeZGnh83PNM4NwNSYx8P1PKXv9DvQgx7yftbr5fNC6PVzRjqL/3zqUlbOImJmE4B1wGh33x/5T3YvMNfd16c1uAQxsxzgEDDN3dcPt3M2s7OAp9y9LfL4a8Bc4Gbgd+4+IbL+KOBVoNrd96Ur3sEwsyKgxSN/rGa2FPgk4byGzc80ysxOB2YDnwOuIvyTNOzOU/puuOWvZInP+2kOJ2N093nh7henN7LM0d3njLs/nN7IMkv855O7/727bbO1BXshsMvd9wNE/ujWAWenM6hEcvf2uFULGUbn7O6PRZNlxBZgI/AGwnlFt9sSuVuTwvASyt2bY5JeMfAK8HeG2c8UDp/f5e5+e8zqhQyz85R+W4h+B3rVRd4Xevy8kIgePmckopvPp25la4F9FFAXt64RmJCGWFJluJ/zScAPGMbnGWmFuR/IB0YyPM/1RuArceuG43lK/+h3QBIp+nkhcbr4nJEOXX0+dStbC2wHmuPWFQIH0xBLqgzbczaz6UCtuy9jGJ8nsAb4GXAu8HWG2bma2RuBJe6+M+6pYXWeMiD6HZCEiPu8kCPFf84IPX4+dStbR9FuAUbHrSsFtqYhllQZludsZrnAtcCnIqu2ADNjnjegmCF+ngDuvg243czaCYMsHmN4/Uw/CSwIPzIgnNvdwDcZXucp/Tcs85ekVhefFxKni88ZCbr8fDKzL7n717p6Qba2YD8MTIh06icyjdtk4JG0RpVcw/WcPwnc6u4HIo8fBGbFPD8VaAWWpDiuZFpCKDiG28/0CuB1MctWwmj/nzK8zlP6b7j9rkt6xH9eSPeinzMSdPf51G1Xo6wssN19O3Av4SsQgPOAu2IGxA151vFvlsHwPGcz+zTwAlBkZkeb2dVAPbDNzOZFNrsQ+La7t6YrzsEys5FmtiBm1UXAbcPtZ+rute6+OboQZkOodfeNDKPzlP4bbr/ryRKf96VDV58Xke4iQvefM+mKJ9P08PnU0N1rsrWLCMAHga+Z2fHAOGBRmuNJGDMrJcxjCXCVmX3H3WsZRudsZjcRpuSLtdLdbzezdwCfM7P1hHlhP5Pq+BJsHvAnM1sJPA086+5/ijw3bH6mvciW85Tu6XegBz3k/azX0+dFOuLJUD19zsgAZOU82CIiIiIiyZKVXURERERERJJFBbaIiIiISAKpwBYRERERSSAV2CIiIiIiCaQCW0REREQkgVRgi4iIiIgkkApsEREREZEEUoEtIhnFzKrM7I9mtsHMPpXueEREpGfK20fShWZEJKOY2RuBh4FRwEpgprvvTm9UIiLSHeXtI6nAFpGMZWaPAm9w97Z0xyIiIr1T3g7URUREMpKZTQL+ku1JWkRkqFDe7qACWxLOzP7NzL6b7jiiEhGPmZ1nZk+Z2cIEhTWsmNmtCd6fAW8HvpLI/YpIZ8rXg2dmU83sR2b2wwTtL9fM3m9mG/qw7cVm9sxAzzWRuVt5uzMV2EOMmc0xs9+Z2Y/NbJmZuZm9O91xxbkP+HW6g4iRiHg2AccnIJaMZGYXmdkuM9tmZqeZ2aTI75mb2XfMrDyy3clm9qKZfcPMiiLrqoHtCQ7pCuAn7t4W2b/IkKN8PSBDMV/vjNwWJmh/OUA9MLkP2y4B5nT3pJnNNbPfRH73bjKzcTHPJTp3K2/HUIE99PwB+IK7XwMsAH4RfcLMjjOzs9MRVOyx3f1Fd388HXF0pT/xdPceuvurwK5ExpXIn9dg9+Xu9wHfA9rc/Sl33wS8H2gFnnH3+sh2zxIS+g3u3hx5+f8D7h3UCcQws88C/wH81cxWArMStW+RFFO+7qehmK/dfT+wNVHHcveDwNI+bruDUIx39/xLwOLIw1vdPbagTljuVt4+kgrsIcTMKoCZwD4Ad28HPge4mZURkrelIa60HTuR+nAe7Sk8Vjr29RtgopmdCODudcBDwJtjjlUObI387kXNdveVgzz2Ye5+s7vPcPcad5/t7o8mat8iqaJ8nVwZmK8TPWNEf+LvbdtDkdv4ftEJy93K20dSgT2ERKa8eRV4wMxOjazbADxL6Pc0DVhkZu8DMLMvmtlnI18PHe5nZWY1ZvZNM/u4me03s1fN7F/N7FIze87M/snMfm9mDWZ2bczrutxf3LHfb2Y/i+2LFvma9Ntm9p9m9rCZnRxZ3+PxYl5vZnatmb1kZueb2XIz22FmV/XhGDNi4+nlmEe8h10YHfnKt9HM3h/3nv5H5L35XSQh9/U9O+JYZna5mT1rZleZ2Woz+0oP++vqZ99lPD2JtHS8BLwtZnUVcLFFuoMAbwH+GBNnEbA/5vFsM/uKmd1qZo9Y6GryDjNbYmaXmNmvzazezK4zsxMj51hnZhf1Fp/IUKJ8nT35uivdnWNPxzGzaWZ2u4XW4K/35ThdHPc+M/uVmZ3Ry3bK3cnm7lqG0ALMBTYQ/mO9HRgb89x6YGHk/jxgbeR+FeG/6zGE/8A3AjWR574LPBm5XwDsBr5DmMvyncCWnvYXf2xCH7QfAj+LrB9BKNrKI48vB2qB0p6O18V5j4sc89ORbW8BWoApvRwjPp4ejxn7HnYRw3rgG8DIyOs2RtaPAv43ZrsnCMmxT+9ZN8eqIHzt99/AycB5Pe0v7mffZTx9/P26CXgtcv8Y4D+BA8BlMb8vFrP9W4CzIvdHAs8BIyOPH428D6OAvcAXgSLgHZHHbyf8k/8JIr+DWrQMpwXl66zI15Hn/6OP72NPP+vngDkxr/E+/p5Ff56nAZ/v4vmFkePkxaxT7i7jTqEAAAV8SURBVE7yohbsIcZDK+N8QsJ5D/CSmXU1mGM1cKWZ5QNnRdaVEP6gJwHR/rOvEP5wcPdWoBG4090bCH9w43rZX3x8B+g8aOJioNkjfXiB3xH+MC/r5Xjx+43u847ItrcQ+gef38sxOsXTn2N24x4P/e2eA8bHnOOYSAvTxwnvKfTxPevmfHcTEtnd7v6su/+1H/vrLp6++A1wtIVuIm8FfgT8GXi7mZUAjR7JwBFnAE9G7r+F8CEWbRX5J+DmyPtcD/zNQ7/tJcAod/8/D1+bv0D/fgYiQ4LydXbk6y50e449HOc8oMzdo/Es6ecxTwLe5e5f6OP2yt1JlpfuAKTvzCwXmOTu64GPmdlPCCOuFwOnxG7r7gfM7CjgRuDn0V24+04zWwKcS/gPexqhqDr80pj7bUS6EXW3v25Cjd3HDCA/Jq52M1sPTOzpeL1x9xYzWwuUA5X9OMaAjxmnjY6/n8nAOnf/7/iN+vKemdmVhEI2apaHr5I9NtZ+/Ay6jac37r7azF4gdBMpc/e1ZvYbQqvS24A/xcQd/d2I9u+bQmhxiu6rLnbXMfcP0Vk76q4mw4zydafzy4Z8Havb97GHn81xhH8mBuoNwAwzK3H3pm62aY+cg3J3CuiNGVryiekf6+7LCF/BzY7f0MzeAHzI3W/xMCNErA8DJ1mYLmo58NXeDtzL/nqykdAiGjt9kRH6Jg5WAbAqycfoi2107qeMmZ3Sj/fsj8DrYpYuR6P3Y39dxtP30+F/gfcBa2PiywWuI3ydGnUKoT9p1BbgTDMrjjluj/0ARYYx5evOsilfd3uOPRynIfKaAgbmVkLr/23dPL/XOwanK3engArsoeejZjYh5vFRwAOR+62EQR2zCVNClZlZoZmdF3l+rJlVEvoC/oTw9c6TRL5yjMgh7r92M7Ne9hd/bIvZxx+AJuCfI/sqjxwv2hLa3fG6Mz6yzVGR1z7Qh2PExtPbMWPPI15XrRlGaJUqBu41szdaGEzT28/g8LGA8e6+MmY5GHOI3Jj7Pe0vNu4u4zGzAjNbHOn+0ZP/I3w1fSeAuzcC9wNPeOfZQy4gdB+Jir7fv7Ywl/ZHgOhcqPE/g/if85Ce0UCkG8rXZE2+7uv72OVxCLm0gPBPGIRvKzCzCX3I3UZoTb4SeJuZvT3y2k+Y2U2E1vuHYrZX7k6FZHTs1pKchTBwwoE9wP8APyN8xTQ68vzngc3AJYSvfFYTWgwuJrR8/IyQrP5IGHDRGtnffsIf3JuAg4SBNJOAf488f3VP+4s79nXAY4R+bXMjz9UQWj6/QeiLOC+yvtvjdXP+Hnn9dYQuC/NjnuvuGFNj4+ntmLHvYdyxLyEM9vtO5L2Ivu59kefPBl6M/Gy+EXmf+/qeXdLFuV4eOd6vgaMi63r6mXbaVzfxTCRcEOGOPvyufSvu8TuAc+LW3drF614PvAzUATdF1l1E+F37NjCByFRlwAcIH8A/ipzrm9P9N6ZFS6IWlK+zKV9PBh4ntMK/rpdz7OlnfTHhm8OnCQPOVwEfipx7l7mb8C3JQULLdQXwN0Jxfw1wLWE+8G8QuqdEX6PcnYLFIm+iZAkzGwt8wN2/GHmcQ2hVudbdb0prcL0wMwemeejTKANkZje4+4CmgBKR1FG+lljK3UOLuohkn48QugtUxKx7AzGD1zJRzFdS+jpqgMws18zeBdyV7lhEpE+Ur0W5e4hSgZ19fkDoq/WKmW0CHgQ2u/tT6Q2re5FBH9dFHl5lYbo46b8S4A/uvrbXLUUkEyhfCyh3D0nqIiLy/9uxYxoAAAAAQf1bW8ITUjgBAEYONgAAjAQ2AACMBDYAAIwENgAAjAQ2AACMBDYAAIwENgAAjAQ2AACMBDYAAIwCFCelCZvr8UIAAAAASUVORK5CYII=\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/orex-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/orex-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/orex-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 }