{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 17 - Apollo-4 - 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 Apollo-4 flight test.\n", "https://en.wikipedia.org/wiki/Apollo_4" ] }, { "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('Apollo-4', 5400.0, 400.0, 0.3, 12.0, 0.0, 3.0, planet)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "# Set up entry parameters\n", "vehicle.setInitialState(120.0,0.0,0.0,10.73,0.0,-7.13,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,70.0)\n", "# bank angle = 70 deg arbitrary, actual profile will give more realistic results" ] }, { "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+17YcXAAAgAElEQVR4nOzdeXxU1f3/8dcnYZVFVlmFgIKgFVAD4g6urdoKiHtdWhW1VltaqdaKWhV/atW6FTXa9utSxa2tBZe617UqoFBbRVRAVgkgOwGSfH5/nBmYhOyZyZ2ZvJ+Px33MzJk7934mQz75cObcc8zdERERERGR+smJOgARERERkWygwlpEREREJAlUWIuIiIiIJIEKaxERERGRJFBhLSIiIiKSBCqsRURERESSoEnUASRLp06dPC8vr9r9CgsL6dy5c+oDaiB6P+lN7yf9pcN7mjFjxgp3z64fbDVqmrPLS4fPqy4Ud8NS3A0vU2OvS9xV5eysKazz8vKYPn16tfvl5+fXaL9MofeT3vR+0l86vCczWxBpABGoac4uLx0+r7pQ3A1LcTe8TI29LnFXlbM1FEREREREJAlUWIuIiIiIJEGjK6zHjRsXdQhJpfeT3vR+0l82vqdslqmfl+JuWIq74WVq7MmO29w9qQes9oRmRwE3AKe4+/xY2xjgNqAt8CjwS3cvjj3XCZgErAaaApe7+9byx83Pz/dMHNsjImJmM9w9P+o4KqKcLSJSVlU5u0F7rM2sC9AaGJbQ1gsYBYwFxgPnAT9PeNlTwJ/c/XLgc+C6BgtYRKQRU84WEamdBi2s3f0b4Nlyzb2BH7v7DHd/GPgDMBLAzIYDe7j7+7F9pwKXmFmrhopZRKQ67k5xcTHFxcVRh5JU6ZSz3aGkBLLsRywiWabBx1i7e2m5x2/Fv0KMWQx8Hbt/ODAvYd/Fsbtp+ZWpiDRO9957L02bNuWSSy6JOpSkS5ec/c470KQJjBhR3yOJiKROOl68OBS4L3a/B7Cy3PPrgO7lX1RYWEh+fv62raCgIMVhiogEubm5AJSUlNT4NQUFBdvyFdApNZE1iAbJ2Tmxv1alpVXuJiKSEjXN2Wm1QIyZ7Q4UuvusWJMDm8rt1hzY4UKYzp07Z+TE5CKS+epSWI8bN27b1ehmtiIlgaVYQ+ZsFdYiEqWa5uy06bE2s1zgAuCKhObFQPuEfQxoBSxp2OhERCpXl8I60zV0zlZhLSKZIG0Ka+Ay4FZ335zQ9jKwR8LjPGALoK5pEUkbjbGwpoFzdrywblw/YhHJNA1eWMd6MAAsoe1K4COgpZn1NbMfm9nu7j4dWGpm34nt+l3gbnff0rBRi4hULpsL63TJ2bEfsXqsRSStNegYazNrA/ww9vAcM7uH8FXi9eV2/czd/xS7fwpwtZnNJ8yn+puGiFVEpKaytbBOp5ytoSAikgkatLB293XAvbEt7obYVtlrFgDnpjg0EZE6y9bCOp1ytgprEckE6TTGWkQkI2VrYZ1ONMZaRDKBCmsRkXpSYZ16GmMtIplAhbWISD2psE49DQURkUygwlpEpJ5UWKeeCmsRyQQqrEVE6kmFdeppjLWIZAIV1iIi9aTCOvXUYy0imUCFtYhIPamwTj1dvCgimUCFtYhIPamwTj31WItIJlBhLSJSTyqsU0+FtYhkAhXWIiL1pMI69XTxoohkAhXWIiL1pMI69TTGWkQygQprEZF6UmGdehoKIiKZQIW1iEg9qbBOPRXWIpIJVFiLiNRTkyZNACguLo44kuylMdYikglUWIuI1JN6rFNPY6xFJBOosBYRqScV1qmnHmsRyQQqrEVE6kmFdeo1bRput26NNg4RkaqosBYRqScV1qnXpEnotS4tBQ1lF5F0pcJaRKSeVFg3jObNw+3mzdHGISJSGRXWIiL1pMK6YbRoEW6LiqKNQ0SkMiqsRUTqSYV1w1CPtYikuyYNfUIzOwq4ATjF3efH2loCtwDrgPbA5e6+trrnRETSQTYX1umUs1VYi0i6a9AeazPrArQGhpV7ajLwtrtfCTwL3F/D50REIpethXW65WwV1iKS7hq0sHb3bwiJdhsz6w6cDkyNNb0C/MDM8qp6rkECFhGpgWwtrNMtZ6uwFpF01+BDQdy91MwSm0YAK9x9Y+z5YjObBxwGbK3iufkNGbeISGWytbCG9MrZKqxFJN01eGFdgR7AynJt64DuQHEVz4mIpIVsLqwrEFnOVmEtIukuHQprBzaVa2tO6Pmgmue2KSwsJD8/f9vjcePGMW7cuCSGKSJSsboU1gUFBRQUFMQfdkp+VCkTWc6OT7enwlpEGlpNc3Y6FNaLCVeOJ2oDLCEk8MqeK6Nz585Mnz49JQGKiFQlJydcruLuuDvlhk5UKLGQNLMVKQ0wuSLL2a1ahdsNG2r1MhGReqtpzk6HeaxfB7rHpmjCzJoCvYB/VfOciEhaMLNtxXUjGA4SWc5u2zbcrlmTjKOJiCRfgxfWtr0rxwDcfRkwDTgi1n4k8Iy7L67quQYMWUSkWtk6zjqdcna8sF6rlQxEJE016FAQM2sD/DD28Bwzu8fdC4GLgFvMbDDQFUgcaFfVcyIiaSE3N5etW7dmVWGdbjl7553DrQprEUlXDVpYu/s64N7Yltj+LXB+Ja+p9DkRkXSRjT3W6ZazNRRERNJdOoyxFhHJeNlYWKcb9ViLSLpTYS0ikgQqrFNPY6xFJN2psBYRSQIV1qnXrl24XVl+CRoRkTShwlpEJAmaNAmXrBQXF0ccSfbq2jXcLlsWbRwiIpVRYS0ikgTqsU69bt3CrQprEUlXKqxFRJJAhXXq7bILmMHy5aAvBkQkHamwFhFJAhXWqdekCXTuDO6huBYRSTcqrEVEkkCFdcPo3j3cLtb6uyKShlRYi4gkgQrrhtG3b7j98sto4xARqYgKaxGRJFBh3TD69Qu3n38ebRwiIhVRYS0ikgQqrBtG//7hdu7caOMQEamICmsRkSRQYd0w1GMtIulMhbWISBKosG4YAwaE2//9D/SjFpF0o8JaRCQJVFg3jM6doVcvWL8e5syJOhoRkbJUWIuIJEFOTkinpaWlEUeS/YYODbcffhhtHCIi5amwFhFJgnhPdbznWlInXlj/+9/RxiEiUl5SCmszOzYZxxERyVTxnupMKKwzPWcfemi4feWVaOMQESmvSW12NrPvAlcBXdlelOcAXYCWyQ1NRCRzxAvr+JCQdJCtOXvoUGjXDr74Ar76avuiMSIiUatVYQ08Avw/4CMgPpDQgNOTGZSISKaJDwVJp8KaLM3ZTZrAEUfAM8/A88/DT38adUQiIkFtC+uv3P328o1m9lmS4hERyUhpOhQka3P2CSeEwnrKFBXWIpI+atu1cpGZja2g/QfJCEZEJFOl41AQsjhnjx4NO+0E77wD8+ZFHY2ISFDbvwAHAQ+aWUnCVgrcm4LYREQyRpoW1lmbs1u3hlGjwv0HH4w2FhGRuNoOBbkOGAdMB+KrIBhwZn0DMbOBwE+BL4B+QIG7f2xmLYFbgHVAe+Byd19b3/OJiCRTmo6xzuqc/ZOfwGOPwb33wpVXQqtWqTiLiEjN1bawfgF41t03Jzaa2QNJiOUR4AR3X2xmvYB/AgOBycCL7v5E7Ar3+4HTknA+EZGkSdMx1lmdsw88EIYPD/NZ//GPcOmlqTiLiEjN1bZrZQrwazM7NGE7DPh1EmIZCLSJ3d8E7Gxm3QlXr0+Ntb8C/MDM8pJwPhGRpEnTHuusztlmcPnl4f7118OaNak4i4hIzdW2x3oi4Su/cwCPtRlhjtSf1TOWpwljAb8H/BC4BBgBrHD3jQDuXmxm84DDgPn1PJ+ISNIUFRUB0LJlWk0PnfU5+4QT4OCD4e23YdIkuOWWVJxFRKRmatu18jtgV3fPc/c+sS0PGJOEWH4CFAEfAuvd/RmgB7Cy3H7rgO7lX1xYWEh+fv62raCgIAkhiYjUzMaNGwHYaaedarR/QUHBtnwFdEpRWFmfs83g9tiEgrffDh9+WKfDiIhUqaY5u7Y91gvdfV0F7V1reZyKtAQeJyTmO2K9HE74ijFRc2Br+Rd37tyZ6dOnJyEMEZHaq21hPW7cOMaNGweAma1IUViNImcPHQrjx8Pvfw9nngkzZ4ap+EREkqWmObu2PdY/Kd9gZicBd9TyOBV5FHjG3a8DbiIk7KWEq8oTtQGWJOF8IiJJsWXLFoqLi2nSpAlNmzaNOpxEjSZnT5oEAwfCnDnwox+Be/WvERFJttoW1h3MbCKAme1kZg8CDxK+CqwzM+sEDHb31bGmG4C1wAKge2z6JsysKdAL+Fd9zicikky17a1uQI0mZ7dsCU89BW3awJNPwm9/m+oziojsqFaFtbsfB7xsZncDHwP9gcFARSt71cYqoMjMeiS0rQRmA9OAI2JtRxJ6SBbX83wiIkmTroV1Y8vZe+0Fjz8exl3/9rfbx16LiDSUKsdYm9mhlTy1AJgL3Ar0BkYB4+sahLuXmtko4GozmwF0ASa4+1ozuwi4xcwGE8YFjqvreUREUmH9+vUAtIp4hRLlbDjuOHjgATjvPPjlL6G0NNyaNWQUItJYVXfx4i2EuUpXA6UVPP8nwtRN3ahHkgZw91nABRW0fwucX59ji4ikUmFhIQCdOqVqco8aU84Gzj0XSkrgggtgwgRYsADuuAPSa+0eEclG1RXW1wHT3X15VTuZWX2/VhQRyVjxwrpz584RR6KcHTduXBhvfc45cM898OWX8PDDEP3/fUQkm1U5xtrdn68uQcf2ezp5IYmIZJZ0KayVs8s67TR45RXo0AFeeAH22ScsJCMikipptfauiEgmSpfCWnZ0yCHw0UdwwAGwaBGMGAG//jVsKj/btohIEqiwFhGpp2XLlgGwyy67RByJVKRXL/jXv+BXvwoXM950EwweDG++GXVkIpJtVFiLiNTTV199BUCfPn0ijkQq07Qp3HwzvPMO7LknzJ0Lhx0Gp54K8+ZFHZ2IZIt6F9Zmto+ZPWpmvzSzvPqHJCKSWeKF9W677RZxJNVr7Dn7gAPCkufXXAMtWsATT8CAAWH2kOXVjk4XEalaMnqsLwSOdPfbgAPjK26JiDQGpaWlzIt1efbt2zfiaGqk0efs5s3h2mvh88/hhz+ELVvg1lshLw/Gj4clKV+AXUSyVTIK658C+8XuTwGaJuGYIiIZYeHChRQVFbHLLrvQpk2bqMOpCeXsmF13hUcegenT4fvfDxc03nEH9O0L558Ps2dHHaGIZJpaF9Zm9j0zOyN2fxjQL75crbuXuvvaJMcoIpK2Zs6cCcCQIUMijqRiytnV228/+Mc/wuwhY8eGHuwHHwwXOI4cCX/7GxQXRx2liGSCWhXWZjYJeAoYDeDuHwBnmdnpKYhNRCTtzZgxA4D99tuvmj0bnnJ27QwZAk89BZ99BpdcAq1bwxtvwJgx0KcPTJwYFpoREalMbXusDwR6Am8ltL0A3JS0iEREMsj06dMB2HfffSOOpELK2XXQvz/cdRcsXgx33gm77x7mwL7hhnB/5MiwiuOGDVFHKiLppraF9cvuvhrwhLasXxpXRKQiW7Zs4e3YUn4HHHBAxNFUSDm7Htq2hUsvhTlz4PXX4cwzoWXL0It99tnQtSuccQY8+ywUFUUdrYikg9oW1l+Z2ZXAHmZ2qpk9A1wM3Jz80ERE0tt7773Hhg0b2HPPPenRo0fU4VREOTsJcnLCio0PPwxLl0JBAQwfDuvXw2OPwahRsMsuofCeOlVFtkhjVqvC2t2nAK8BucBZwGrgaHf/QwpiExFJay+99BIARx11VMSRVEw5O/l23jnMGPLee2G89U03wb77wrp18Oij8IMfQMeOcMIJ8MADmrpPpLExd69+r6oOYJYLjHD3V5MTUt3k5+d7fKyjiEiquTv9+/fniy++4KWXXqpXcW1mM9w9P4nhVXUu5ewU+OKLcOHj00+HBWgS7bMPHH88HHMMDBsWVoEUkcxVVc6usrA2s38Aras5fidghbsfXvcQ6y/bkrSIpLf333+f4cOH07VrVxYtWkRubm6dj5Wswlo5Oz0sXgzPPw/TpsErr8DGjdufa90aDj0UjjgibHvvHYaaiEjmqCpnN6nmtbNj+8wBSoHzgY+AxGzYB+iQhDhFRDLGn//8ZwBOPfXUehXVSaacnQZ69AjDRc4/P4y3fuMNeO65UGR/9lkoup9/PuzbqVOYZeSQQ+Cgg2DQIGhS3V9mEUlb1f36/o7Qq70awMx2c/ery+9kZs+mIjgRkXS0YsUKHnroIQDOP//8iKMpQzk7zbRoAd/9btgg9Ga/9loosl99NTx+6qmwQejRHj48FNkHHwz77w+ZsaCniEA1hbW7rynX1MXMct29JN5gZiOAYSmITUQkLU2ePJmioiKOPfZY9txzz6jD2UY5O/316BFmDznzTHCHzz8PPdrvvBO2r74KRfcrr4T9c3JCL/bQoWHLz4fvfEfjtEXSVW2/cHoamGVmbwPrgIHAMcAOPSIiItlo+fLl3HrrrQBMmDAh4miqpZydxsxgjz3CdsEFoW3pUnj3XXj77VBof/QRfPxx2B54IOzTvHlYJTI/f/s2cCCkz4gkkcarVoW1u79sZkcDPwEGEaZuGuvu+lpRRBqFa665hnXr1nHssccyYsSIqMOpknJ25unWDU48MWwQVnecOROmTw/bhx/C3Lnw/vthi2vRIvRkDx4ctkGDwta+fTTvQ6SxSsZ0eznAPu4+Izkh1U02X2EuIunh7bff5tBDDyUnJ4fZs2cnbRhIA0+3p5yd4VavhhkzyhbbCxZUvG+vXqHATiy2d9tNQ0lE6qM+s4KUP9CfKmjuDJQAo+oQW0XnGA4cDHwBvA1sAG4hfI3ZHrjc3dcm41wiIjW1fv16zj77bNydX/3qV2k1troyytnZqV277dP1xX37LcyeHbZZs8L2ySfw9ddhmzZt+75Nm0L//mH4yJ57hm3gwNDWokXDvx+RbFLbMdb9gFeBxG7unsBbyQjGzH4E7ObuVyW0/Rl40d2fMLPvAvcDpyXjfCIiNeHu/OhHP+Krr75i0KBBXHPNNVGHVFPK2Y1E+/Zw2GFhiyspCcNGZs3aXnB/8kno3f7vf8OWKCcH+vYtW2wPGBAK7nbtGvb9iGSq2hbW57n7nMQGMxsAnFPfQMzsMOAM4KiEtu7A6cDFsaZXgGfMLM/d59f3nCIiNXH99dfz9NNP06ZNG6ZMmULz5s2jDqmmlLMbsdzcUBgPGACnnLK9ff16mDMH/vc/+PTT7bdffLF9+8c/yh6rU6dQYPfrF7b4/d13D1MEikhQ24sX51TQvBk4D7iinrHcDrwL3G1muwPXA70JK4RtjJ2/2MzmAYcB8+t5PhGRak2ePJlrrrkGM2PKlCkMHDgw6pBqTDlbKtK6Ney3X9gSbd4cpv+LF9v/+194PHcurFgRtnff3fF43bvvWHDvthv06aOiWxqf2o6xnkfZrxRzgS7AI/UJwsz2APYBfuTus81sAvAicAOwstzu64Du5Y9RWFhIfv72ceTjxo1j3Lhx9QlLRBq5++67j4svDp2v99xzD8cee2zSjl1QUEBBQUH8YaekHTiBcrbURvPmYYn1vfcu2+4OS5ZsL7Lnzt1+/8svw3NLlsC//rXjMTt3DsNL+vTZ8XbXXbXKpGSOmubsWs0KYmYFwGNsT9SlwDfu/nkd44wfdwxQ4O6dYo9bAoXAdcCJ7r5/wr4zgcfc/dbEY+gKcxFJFndn4sSJTJo0CYA77riDn/3sZyk7X6pmBVHOllQrKQkXR5YvuL/6CubPD73glcnNDbOW9O27Y9GdlxeKcrOGeiciNZe0WUGAR939zXIH38PMdnf3L+ocYYhjWyzuvsnM5hJ6V8rPwtkGWFKPc4mIVGrdunWcf/75PPHEE+Tk5DB58mQuiK/ekXmUsyWlcnNDIdynDxx9dNnnSkvDgjfz5oVCO34bv794cbidNy8s715eixah8O7VC3r3Lnvbq1fo8W7WrGHep0hN1baw/i7wZrm2b4DXCV8L1tVsYGcz6+Lu38TaignJuLuZtYwl7qZAL6CCL5xEROrno48+4pRTTmHu3Lm0atWKJ598MqnDPyKgnC2RyckJS7j36AEHH7zj80VFYYaSxKI7frtgQZhC8PPPw1YRM+jatWyxXb4Ab9dOvd7SsGpUWJvZTwlj59qY2eUV7PJyfYJw98/M7AVgLPAHM2sHNAP+AnwPOAKYBhwJPOPui+tzPhGRRJs3b+bmm29m0qRJbNmyhb333psnn3ySAQMGRB1anShnSyZo0WL7ku4VWbdu+zzcCxaUvf36a1i0KPSIL10K//53xcdo3Tr0bPfoAT17bt8SH3fsqOJbkqdGhbW732NmbwJnA3eVe3qLuy9NQixnAXeaWTNCD8dpsSvKLwJuMbPBQFdAV7eISNK8+eabXHjhhXz66acAXHjhhdx+++20bNky4sjqTjlbskGbNrDXXmGrSHFxuGiyoqJ7wYKwrV8fZjmJ/XpXqHnzsoV2RUV4165h2ItIdWp78WJTd99aQXsvd/86qZHVki6EEZHa+PTTT/n1r3/Ns88+C0D//v25//77GTFiRIPHksKLF5WzpdFyD8NJFi8OvduLFpW9H3+8enX1x8rNhW7dthfd3buHx/Hb+P0OHdT73RjU+eJFM7sMeN7d/xdr+pWZlf8/mwEjgRH1DVREJNU+//xzbr75Zv7v//6P0tJSdtppJyZMmMAVV1xBiwxfz1k5W2Q7s1Doduiw4xSCidav315wV1R4L1oEy5dvb3v//cqP1axZ6N1OLLgrKsA7dQpj0CX7VDcUpA+QuJDpvkArYBnbp28ywtd9IiJp68MPP+Tmm2/mr3/9K+5Obm4uF1xwAddccw3dunWLOrxkUc4WqaXWrase6w1h2sAlS0KhvXDh9rHdS5aUvb927fbhKFVp0gS6dKm6AO/WDXbZRXN9Z5oqPy53v7hc0yTgE3ffkthoZnsmOzARkfrasGEDjz/+OPfffz/xYQfNmjXjrLPOYsKECfTv3z/iCJNLOVskNZo33z6tYFU2bqy86I7fX7oUVq0KRfriai7rNQvzeXfrFnrCE2/Lt2mVy/RQ2/8H7e3uMxMbYitw7Qb8r+KXiIg0nNLSUt5++20ef/xxHnvsMdauXQtAu3btOO+88xg/fjzdu++wEGC2Us4WaUA77RSWc99tt6r3KyqCZcsqLrrjj5ctg8LCMAxl+XKYNavqY7ZuXbMCXMNQUqum0+0dQvj68KjYErmJOhCuOp+a5NhERGrE3ZkxYwZTpkzhiSeeYNGiRdueO+CAA7jgggs4+eSTM3qmj9pQzhZJby1ahNUl8/Kq3m/r1lBcx4vueDFe/nbp0jBWPL4CZlXiw1CqK8C7dg1xSu3UtMd6NfAwYfzeQeWe2wI8mMygRESqs2HDBl599VWee+45nnvuORYnfKfau3dvTj31VM444wz2ruqqpeylnC2SBZo2DeOuq/uSzR3WrNmx2K6oAI/PlFLdMBSA9u2rL8C7ddNCPIlqOo/1f8xsBHCSuyshi0iDKykp4eOPP+b111/n1Vdf5fXXX2fz5s3bnu/evTtjx47l1FNPZfjw4VgjzvLK2SKNi1kobtu1g4EDq9538+ZQaFdVgMef//bbsFU1DziEMeg1KcC7dMn+izFr/PbcfQ2V9HKY2cHu/nbSohKRRm/r1q3Mnj2bN998k9dff50333yTNWvWlNln2LBhHH/88Rx//PEMGTKkURfT5Slni0hFmjcPS7737l31fqWlsHJlxb3e5dvWrdu+KE9VzMIY7+oK8G7dMvdizOrmsf4HUN1bM2AvYJdkBSUijYu7s2DBAt5///1t28yZMykqKiqzX58+fRgxYgQjR47k6KOPpkuXLhFFnJ6Us0UkWXJywowknTtXPQ84wIYNVY//jt9fvjyMGS8shNmzqz5mfDn6xK1Xr7KPd9opee83WarrsZ5DGKu3mO1zoJZnwPHJDEpEsldRURGffvops2bNYvbs2cyePZtZs2axYsWKHfbt168fBxxwACNHjmTkyJH0rq6LRZSzRaTBtWpVs9lQiou3X4xZVQEevxizuuXoO3TYseju3Rt23x369QtDYxpadYX1zcA6d99c1U5mpqvLRaSMNWvW8PnnnzNnzpxtt5988glz5syhpKRkh/07duzI/vvvz/7778+wYcMYNmwYHTp0iCDyjKacLSJpq0mT7UM9quIelppfuLDs9vXX2+8vWhTmA1+1qvKpCDt3DgV2v37Qvz8MGgRDh4ax3qlS3QIxO3YhVWw0uspcpFEpLS1l6dKlLFiwgPnz57NgwQK+/PLLbUX08uXLK3xdTk4OAwYMYNCgQQwePJhBgwYxaNAgdt11V42RriflbBHJBmZhRpL27UMxXJHS0tD7nVhsL1wI8+dvn3YwPuzk3XfLvrZXLzjqKDjvPBg+PLmx1+vaTDPrBVwM/BQlaZGsUVxczPLly1m6dGmZLbGIXrhwIVu2bKn0GC1btqRfv37sscce9O/fnz322IOBAwey1157NZr5pNONcraIZIucnNDz3KVL6IUur7Q0LLYTL7LnzIGZM0OR/fXX8Mc/wp/+BK+8Aocfnry4al1Ym1ku8APgfOAo4BvCmD6RWnF3Nm3axMaNG9myZUuZrbi4mBYtWtCiRQtatmy57bZp06bq1ayDkpISvv32W1asWMHKlSt32AoLC8sU0MuXL6e0tLTa43bu3JnevXuTl5dH79696dOnz7ZCumfPnuRoea/IKWc3LiUlYWzqhg1hPGtJyfattLTs/dzc8NV8TTf9Oks6KikJ/97XrAlTA65aVfa2sDDM2b1oUbjdunX7a93D70ky1biwNrM+wHnAj4COwALg+8A/gQOSG5ZkuvgsD3PmzNnWyzl//nwWL17MqlWrWLlyJatWrSozD3FNmFmZQruq25rsk7hv8+bNt90m3k9sS3Wh6O5s2bKFoqKiKrcNGzawdu1a1q1bx9q1a8vcT7xds2YNK1euZPXq1bhXdi1bxTp37ky3bt3KbL1799629erVi1atWqXoJyH1pZydndasgX//G2bMgHnzwvRmCxeG9rVrQ4GRKma1K07ZPMEAACAASURBVMQTt6ZNwxa/X9PbZLwmN3f7lpNT9rau9xtb/05paShAy28lJRW3V/Tc1q1hDu3abkVF4d91/D+M5e+XmzyqWmZhzPXQoXDGGXD00cn9WVVbWJvZSYSejsOBL4DbgYeAK939xdhu71bycmkE3J2vvvqKt956iw8++IDZs2fzn//8h7Vr11b72mbNmtGmTRuaNWtWZsvNzWXz5s1s2rSJoqKibbdbt25l48aNbNy4sQHe2Y6aNm1K8+bNyc3NJScnZ4fbitpKSkooLS2lpKRk21b+cXwrP71cMrVv356OHTvSsWNHOnXqtO1+/HFiAd2lSxeaNm2aslgkdZSzs9P06XDddfDcc6HIqUqbNmGWhvJFZfnCsrJiqaJt69bQu7d1a9kev8bKrOriO154m23fEh9Xdr8uz7mHrbS06tu67BP/N5LOzMK0e23bhllCOnQIY7Pj9zt0gJ49w9ajRxhfncpp+mrSY72JsATuu8B57v55eCNWu+4vySrLli3j+eef56WXXuKtt95iyZIlO+zTpUsX9txzT/r06bNtuEDPnj3p1KkTHTp0oEOHDrRs2bJWQzuKi4srLLjrc7tp0yY2b97M5s2bKSoq2uF+YtvWrVvZmuK/Ks2aNdvWi17ZttNOO9G2bVvatm1LmzZtqrzt2LEj7du3Jzc3N6VxS9pQzs4yzz0HJ5wQegCbNIFhw8IFV/36hanFevUKxUObNmHu31R9sVabQjyxIK/pbW32rc1rKhsKU5f78WIzFUMI0l38G4D4VtlQoqqGGDVvXvutRYvw77pVq7JbvK1ly/T6BqHawtrdpwHTYhe9nG9mecCzwLbuLDMzr+33zJJxPvvsM5566immTZvGBx98UOa5jh07cvDBB3PQQQcxZMgQBg0alJLFO5o0aUKTJk0iGYIQH6axefPmbT3OibeVtZXvxY5v5R/n5uZu6w0XqSvl7Oxz/fWhqDvvPJg0CXaJaGmfnBxo1ixsjVlib25FRXh8NtF4D3D8N626+3V9zmz78JTKbuuzj8bX105tljT/GphoZk2AUUBPM7se+CNwNFCQmhAlSkuWLGHKlCn85S9/YebMmdvamzdvzhFHHMGxxx7LyJEjGTBgQNZfqGZm28Zai6Q75ezsEe8ZPfro6Ipq2S4+DCQ3N/TiiiSq9awg7l4MPA08bWa7A+MJ4/mUpLNEcXExzz33HPfddx8vvfTSttkh2rZty5gxYxg9ejRHHHGELlwTyQDK2ZnvrLPCxYo/+QkMHhwWuhCR9FSveazd/QvgZ2Y2I0nxSIQWL17Mgw8+yAMPPMDixYuBMN73uOOO44wzzuC4446jRYsWEUcpInWlnJ2ZLr4Ynn8e/vlPOPTQcDt4cNRRiUhF6lVYx7n7w8k4jpk1Az4Efubub5hZS+AWYB3QHrjc3aufakJq5YMPPuC2227jmWee2bbUdL9+/bjwwgs5++yz6dixY8QRikgyJStng/J2Q8jNhWeegVGjwmIWBx0EDz8MY8ZEHZmIlJdug2InAHkJjycDb7v7lYSLb+6PIqhsVFJSwt///ncOOeQQ9t9/f5588knMjJNOOolXX32VOXPm8Itf/EJFtYhUR3m7AbRqBVOnwg9/GObuPfFEuPrq6qfeE5GGlTaFtZkdCCwFvo097g6cDkyN7fIK8IPYFe5SR0VFRdx3330MHDiQ0aNH8/bbb7PzzjszYcIE5s2bx5NPPsnhhx+u1Q1FpFrK2w2rRYvQU33rrWGWhuuvDxc0VjDbqYhEJC0KazNrBYx19z8lNI8AVrj7Rth2Ac484LCGjzDzrV+/nttuu40+ffpw0UUXMXfuXHr37s3vf/97Fi5cyC233ELPnj2jDlNEMoTydjTM4Je/hBdegM6d4dVXw3jradOijkxEIE0Ka+By4KZybT2AleXa1gHdGySiLLFmzRomTZpEXl4el112GcuWLWPIkCFMmTKFL774gp///Oe0adMm6jBFJPMob0fo6KNh1iw46ihYsQK+//1wkWNEi9KKSEzkhbWZHQNMd/fl5Z5ywgpiiZoDFS57V1hYSH5+/ratoKBxzyS1YsUKJk6cSO/evbnqqqtYuXIlw4cPZ9q0acycOZNTTjmFJk2Scu2qiNRBQUHBtnwFdIo6ntpIRt5Wzq6/bt3gxRfhd78L8ylPngz77QcJSw6ISJLUNGdb1ItvmdnLwD4JTe2B9cDtwBnu3j9h37nANe7+WPnj5Ofn+/Tp01MdbtpbtmwZt912G/feey8bNmwAYOTIkVx11VWMHDlSY6dF0pCZzXD3/KjjqKlk5G3l7OT6+GM4/XT49NOwUt7118OECWFGERFJrqpyduQ91oQLXYYkbEuA84A/A91jUzdhZk2BXsC/IoozrS1btoxf/OIX9OnTh1tvvZUNGzZw7LHH8s477/Daa6/pgkQRSSbl7TQzZEhYROaSS8JKjb/+NYwcCQsWRB2ZSOMSeWHt7oXuvii+ASVAYWw53mnAEbFdjwSecffFUcWajpYuXcr48ePp06cPv//97ykqKmLUqFHMmDGD5557jgMPPDDqEEUkyyhvp6eWLeGuu8KFjV27wltvwaBBYUEZEWkYkRfW1bgIOMHMfgMcC4yLOJ60sXTpUn7+85/Tt29f7rjjDoqKihgzZgwff/wxf/vb39h3332jDlFEGifl7Yh997vwn//ACSfA2rVw/PHw0ENRRyXSOKTd1Wvunpdw/1vg/OiiST9Lly7l5ptv5v7776eoqAiAMWPGcPXVVzNYa9yKSASUt9NPp07w17/ClVfCzTfDOedAs2Zw2mlRRyaS3dKusJaKLVmyZFtBvXnzZgBOPPFErr76agYNGhRxdCIikm5ycuCmm6B9e7jiCjj7bOjXD/Iz5jJZkcyT7kNBGr3Fixdz6aWX0rdvX+666y42b97MiSeeyKxZs3j66adVVIuISJV+9Su46CLYuhXOO0/LoIukkgrrNBUvqHfbbTfuvvtuNm/ezNixY5k9e7YKahERqTGzMNd1z55hUZnXX486IpHspcI6zSxevJhLLrmkTEF90kknMXv2bJ566in23nvvqEMUEZEM06rV9vHVb74ZbSwi2UxjrNPEokWLuOmmm3jggQfYsmULACeddBITJ05UMS0iIvXWtm243VR+bUwRSRoV1hFbunQpN954IwUFBWzZsgUz4+STT2bixIl85zvfiTo8ERHJAqWl8Le/hftDhkQbi0g2U2EdkVWrVnHLLbdw1113sWnTJhXUIiKSMjfcADNnQrduYX5rEUkNFdYNbP369dx555387ne/Y82aNQCMHj2a6667TgW1iIgkVWkpXHNNKKzNoKAgjLcWkdRQYd1AioqKuP/++5k0aRKFhYUAHHXUUUyaNImhQ4dGHJ2IiGSb+fPh3HPhtdcgNxfuvz+swigiqaPCOsVKSkp46KGHuPbaa1m4cCEAw4cP58Ybb2TkyJERRyciItlm/Xq47bYwxd6GDWEVxkcfhWOOiToykeynwjqFXnvtNcaPH8/s2bMBGDRoEJMmTeK4447DzCKOTkREssm6dfDAA3DLLfDNN6HtpJPgnntgl12ijU2ksVBhnQJz585lwoQJPPvsswD07t2bG2+8kVNPPZWcHE0dLiIiyfPNN6F4vuceWL06tO2/f+ixPuSQaGMTaWxUWCfRmjVruO6667j77rvZunUrrVq14sorr2T8+PG0bNky6vBERCRLuIeFXu67D555JixXDnDQQXD55WEstb4YFWl4KqyTwN155plnuPTSS1m6dClmxo9//GNuuOEGunXrFnV4IiKSJVasgL/8JVyI+OmnoS0nB37wA5gwAQ4+ONr4RBo7Fdb19PXXX3PxxRczbdo0IFyY+Ic//IF999034shERCQbbN4M06bBww/D889DcXFo79oVzjsPzj8fevWKNkYRCVRY15G7c9999zFhwgQ2bNhA27Ztuemmm7jgggs0jlpEROrFHd57Dx55BJ54Ar79NrTn5sL3vgc//nFY6KVp02jjFJGyVFjXwfLlyzn33HO39VKPHTuWO++8k+7du0ccmYiIZLKvvgpT4z38MHz55fb2ffaBM8+E004LPdUikp5UWNfSSy+9xJlnnsny5ctp164d999/PyeffHLUYYmISIZavRqefDL0Tr/99vb27t3hjDNCQb333tHFJyI1p8K6htydW2+9lSuuuILS0lJGjBjBww8/zK677hp1aCIikmG2boUXXwzF9D/+EcZRA+y0E4wZA2edBYcfHoZ+iEjmUGFdA0VFRZx77rk89thjAFxzzTVMnDiRXGU8ERGphcWLw4weDzwAy5aFNjM44ohQTI8ZA61bRxujiNSdCutqrFu3jhNOOIHXX3+dVq1a8cgjjzB69OiowxIRkQzhDm+8AX/4A/z971BSEtoHDIBzzgnDPXr2jDJCEUkWFdZVWLlyJd/73vf48MMP6datGy+++CKDBg2KOiwREckAW7fC44+HFRA/+SS05eaGZcYvvhgOPVSLuIhkm7QprM1sDHAb0BZ4FPiluxebWSdgErAaaApc7u5bUx3P+vXrtxXVffr04eWXX2a33XZL9WlFRDJCuuXsdLJhAzz4INx2GyxcGNq6dYNx48KmCaREsldaFNZm1gsYBYwF9gLuBRYCtwJPAVe4+/tmdiFwHfDrVMazZcsWxowZw4cffkheXh5vvfUWPXr0SOUpRUQyRrrl7HSxaRPceWfooV61KrQNGBCWGD/9dGjWLNr4RCT10qKwBnoDP3b3YmCGmX0HGGlmbwN7uPv7sf2mAnPM7AZ335CqYMaPH8/LL79M586deemll1RUi4iUlVY5O2olJWF2j4kTYdGi0DZ8OFxxBXz/+2HJcRFpHNKisHb3t8o1LQbaAIcD8xL2W2xhQFo+8K9UxPLYY48xefJkmjVrxrRp0+jXr18qTiMikrHSKWdH7bPPwiqI770XHg8ZArfcAkceqfHTIo1Ruv4/eihwH9ADWFnuuXXADiPUCgsLyc/P37YVFBTU+qQLFy7kggsuAOCOO+5g2LBhtT6GiEhNFBQUbMtXQKeo46mnSHJ2lNzh9ttDIf3ee2EM9cMPw4wZcNRRKqpFsk1Nc3Za9FgnMrPdgUJ3n2VmDmwqt0tzYIcLYTp37sz06dPrfF535+KLL2b9+vWMGjWKCy+8sM7HEhGpzrhx4xg3bhwAZrYi4nDqLKqcHaU1a8I0eX//e3j8ox+FIrtdu0jDEpEUqmnOTqvC2sxygQuAK2JNi4H+Cc8b0ApYkuxzP//880ydOpU2bdpwzz33YOpuEBGpUpQ5Oyqffw7HHQdffAE77wwPPQQnnBB1VCKSLtJtKMhlwK3uHlvclZeBPRKezwO2AEnt5igtLeWqq64C4Nprr9XFiiIiNRNJzo7K++/DgQeGonrw4DDsQ0W1iCRKm8LazK4EPgJamllfM/sxYR7UpbErzgG+C9zt7luSee6pU6fy8ccf0717dy666KJkHlpEJCtFmbOj8OabcPjhsHIlfO978M47oKUNRKS8tBgKYmZXAdeXa/7M3f9kZqcAV5vZfKA18Jtkn3/y5MkAXHbZZbRs2TLZhxcRySpR5+yGNmtWmDZv40Y480z44x+hadOooxKRdGTuHnUMSZGfn+91uRBm/vz59O3bl2bNmrFkyRI6dOiQguhERCpnZjPcPT/qOBpSXXN2Q1u9Osz8sWABjB0LU6aEZclFpPGqKmenzVCQqDz99NO4O6NHj1ZRLSIiZVxySSiqhw6FRx9VUS0iVWv0hfXUqVMBGD16dMSRiIhIOpk+PRTTLVrAX/4CzZtHHZGIpLtGXVivWbOGd955hyZNmnDMMcdEHY6IiKSRm28Ot5dcAlqEV0RqolEX1v/+978pKSlh2LBh7LzzzlGHIyIiaWL1avjHP8IKipdeGnU0IpIpGnVh/c477wBw4IEHRhyJiIikkxdegC1b4LDDoGfPqKMRkUzRqAvrd999F1BhLSIiZf3zn+H2+OOjjUNEMkujLazdnY8++giAoUOHRhyNiIiki9JSePHFcF+X34hIbTTawnrp0qWsWrWK9u3bawlzERHZZvZs+OYb6NED9tor6mhEJJM02sL6k08+AeA73/kOZhZxNCIiki7iw0COOSZcvCgiUlONtrD+z3/+A8Dee+8dcSQiIpJOXn453B59dLRxiEjmaRJ1AFHJy8vj+OOP5+CDD446FBERSSOHHgpr1sDhh0cdiYhkGnP3qGNIivz8fJ8+fXrUYYiI1JqZzXD3/KjjaEjK2SKSqarK2Y12KIiIiIiISDKpsBYRERERSQIV1iIiIiIiSaDCWkREREQkCRpdYV1QUBB1CEml95Pe9H7SXza+p2yWqZ+X4m5YirvhZWrsyY5bhXWG0/tJb3o/6S8b31M2y9TPS3E3LMXd8DI1dhXWIiIiIiJpKGvmsTazQmBBDXbtBKxIcTgNSe8nven9pL90eE+93b1zxDE0qFrk7PLS4fOqC8XdsBR3w8vU2OsSd6U5O2sKaxERERGRKGkoiIiIiIhIEqiwFhERERFJgkZTWJtZSzO728xuNLN7zaxt1DHVh5mNMbN5ZrbSzO40syZRx5QMZtbMzGaZ2YioY0kGMxtuZpeZ2Sgz6xR1PPVhZgPN7A9mNt7MJpvZkKhjqi0zO8rM3jezvIS2rMoN2SpTP6dMz9WZmpMzMfdmUo7N1FxaSdxJ/R3NqF/wepoMvOjuT5jZd4H7gdMijqlOzKwXMAoYC+wF3AssBG6NMq4kmQDkRR1EMpjZj4Dd3P2qqGNJkkeAE9x9cezf4D+BgRHHVGNm1gVoDQwr91TW5IYsl3GfU5bk6ozLyRmcezMix2ZqLq0o7lT8jjaKixfNrDswD2jv7htj/xtZA+zl7vMjDa4OzOwQ4D13L449voXwXo6LNrL6MbMDgQHA1cA57v5GtBHVnZkdBkwEjvIs+SUzsw3Afu7+mZl1Bma5e/eo46oNM8sBSoA+7j4/23JDtsrUzynTc3Um5uRMzr2ZlGMzNZdWEHfSf0cby1CQEcAKd98IEPsBzgMOizKounL3t+L/CGIWA19HFU8ymFkrYKy7/ynqWJLkduBT4G4ze9HMDoo6oCR4GnjQzNoAPwQuiTieWnP30nJNI8ii3JDFRpCBn1Mm5+oMzsmZnHszJsdmai4tH3cqfkcbS2HdA1hZrm0dkJb/E6yDocB9UQdRT5cDN0UdRDKY2R7APsAD7v5T4FXgxdjXUJnsJ0AR8CGw3t2fiTieZMj23JAtsuVzyqRcnXE5OQtybybnWP2OxjSWwtqBTeXamgNbI4glqcxsd6DQ3WdFHUtdmdkxwHR3Xx51LEmyF7DK3WfHHt8DGHBCdCElRUvgceAx4A4zOzLieJIha3NDlsn4zymTcnUG5+RMz72ZnGP1OxrTWC5eXAy0L9fWBlgSQSxJY2a5wAXAFVHHUk+XAfuYWfxxe+BZM5vk7rdEF1adNSHhd8vdN5nZXKBjdCElxaPAqe6+2swceNzM8tx9Q9SB1UNW5oYslNGfUwbm6kzNyZmeezM5x+p3NKax9Fi/DnQ3s5YAZtYU6AX8K9Ko6u8y4FZ33xx1IPV0OjAkYVsCnEfmfGVa3mxg53JfPxYD/4sonnqLTVc12N1Xx5puANaShles11K25oZsk+mfU6bl6kzNyRmbe7Mgx+p3NKZRFNbuvgyYBhwRazoSeMbdF0cXVf2Y2ZXAR0BLM+trZj+OfY2Rcdy90N0XxTfCFbuF7r426tjqwt0/A14gTN+DmbUDmgHPRRlXPa0CisysR0LbSuDziOKpE9veBWeQnbkhG2Xy55SJuTpTc3KG596MyrGZmkvLxx1rS+rvaGMZCgJwEXCLmQ0GugLjIo6nzszsKuD6cs2fZeDV29nsLOBOM2tG+F/7aeWuPM4o7l5qZqOAq81sBtAFmJDuf2gTJVxpD3COmd3j7oVkUW7Ichn3OSlXRyIjc28m5dhMzaUVxU0Y/pHU39FGMY+1iIiIiEiqNYqhICIiIiIiqabCWkREREQkCVRYi4iIiIgkgQprEREREZEkUGEtIiIiIpIEKqxFRERERJJAhbVIJczsCDN718xGVLPfd8zsSTObmMRz/y5ZxxIRiYKZDTKzp5KZGys5zzQzOzWV50gVMzvSzN6r7u+MZI7GtECMZCgzGwhMIqxMNRQYBJzp7o+m+NQLgcE12O9roD2Qm4yTxpaC3ZqMY4mIAJjZUOD/AYcCLwHzCIt4NAXudfd/puC0XwA7k6TcWIXJwCepOHAsH48hLOf+OfBe7KnmwInAMHefX49T1PTvjGQIFdaSCf4OnOzus8wsB/hzQ5zU3T83s5U12G+tmS1N4qkPA95M4vFEpJFz9w/NbAowwN2Pj7eb2XHA38zsMnefnORzbkxybgTAzL4P/Cde0Lr788k+R5y7bwWeMLNLgDfc/aqEOJ6tyzHNbBdgpLs/4e5zzGxFksKVNKChIJLWzKwD0B/YAGHZV+BqoKGWDC1N8n41MRJ4PYnHExEB2GFpb3d/DrgW+L2ZdU3BOZOaq82sJ6GHuqFVtCz6K8Di2hzEzJoDDwMtE5qT+fdDIqbCWtKau68C5gAvmtnwWNsC4AMza2dm15vZB2b2AzNbZmafmtm+8deb2Wlm9hsze8PMbkxozzeza83sCTP7q5ntHGtvYWZ/MLOrzeyPhK8xa8XMfhWL6WQzG2pmf48d7/pYjE+YWedY+xozu6bcIZq5++bYsX5pZr8ws3fM7KZa/wBFRKr3R8KQkJMAzGw3M7vOzB40s1djxSxm1jyWy35jZm+Z2SnxA5jZT8xsopn928x+WtmJKsrJZjY2ljPPMbO5ZnaTmXU3sz/H8ulbZnZ47BCnAz2By81sVOxamFfM7OyEc5xvZpPM7GEze8jMWpnZzmZ2jZlNN7PDzWyOmX1lZnl1+YHFfiYjYj3amNkYM7vTzG4zsxfNrHXsPdxpZlNjf2+Wxn7GewInmdkvEg7ZPva3aJ2ZnV/JOXPM7GYzGxc7h5vZ3+Kfj6QJd9emLa03YC9gAeF/9X8Cdom15xKS1CrCWLd2wMvAfwn/aTwAuCq2b3tCr/dxQFtgSsLx3wF+F7v/e2Bc7H4XoISQPKuL8f8IvT67EMbi7RRrbw78E3gW6E74g7AVuAbYCdgv9ji+f3z8ePx9T47dbwNcEfVnoU2btszdgHOARZU89w1wVyyv/hXIjbU/BjwVu387MCZ2/8dAYez+acAPY/f7x3L13rHH/wdcG7tfWU7uAKwG7gCGAUcCtwJXx/adAExNiNWBvNj9lsCnwDmxx8cAjyXsOw24nzD0dXTs78WY2Pt8EZhUw5/dG7G/FTcBvyOMtz4y4fkVQK/Y/f8S/iY1B34DLAL2B86N5f034vHG9p8P3BZ77jTg60piOA+YFrvfJfZzPjrqf1faym4aYy1pz93/a2Z7AzcAPwG+b2ZHehhzXQisdfdnAMzsWuBtoC/hj0ipmf08dqh/EhL4cUDHhPZPY69tDVwUey3u/o2ZLaxFqJ2Au4Fz3X1j7BibzWwZMM/dl8TOswz4l4fxhx8REn4nwkWQxxH+CABsBM4ys8+Bewn/qRARSQUnFJvDgN7AJWYGsBLY2cKDCwidAgAPEYZCQMi108vl2m7Af8qd4xwqyMnuvsrM1gDPuvsHAGb2NbDBzLoRenhbVxi0+yYzW57Q9NPYceMeBJ6Mta8m/L34a+wc02Nx1tTrHhtjbWZPEDpz4k5y96/N7CBCQd06lv+XAF+6+/vA+7HXVnTsqbG/CR9WEdO+wKbY+/7GzFZRdkiJpAEV1pLWzCwX2NXDRSqXxoZnPA8UEHoAyvtv7LYd0At4xN0fi7XdETvm5YRC945y59qPkBDX1THc/sDBwG7ArIT28mMMS7Y94V4aS7LxYVkdPQx/wd3nmdk5wD3AeOAMYDkiIklkZu0I37b9l5A3V1aQH7sQelSbAbh7CaEzgNhrfuvu78Yel3ltggpzcoxTNlcuAi4H5gL/BvpU8RYSX9cPeC3h8ZeEYS67sGMuLqbuM5bMJPwHJG61md0MPEWIPV49l39f1Smm8trsBeBeM2tB+Cw2EXq/JY1ojLWku23j/gDcfRZwJTCgkv2bEZLYF8BSYGzik2a2f6z9ODNrWa59bezhwDrG+i5wJ/CXWOKrFTPrTsKFMGbWgzAjSn/CxYyPVfJSEZH6OJvwDdlThPx4UKynGNiWH1cAW4DvJrTvbWZtKJdrzaxJrKOivMpyckXuA+Z7mFZ1cy3ey9eEHu5tpwDWEIa6JI27l8Y6P/aK/cfkeeA6d5+ezPOUO+dUYCrhm9XTgMPcfU2qzid1o8JaMsElsaIzrgdhbFxcBzNrFrs/gjC+bjXwODDazO4xsxGx3oQNhATYCphmZsfELhQZQCjGPwKuN7NmZtaJMLZ5FzPLNbMrzOyHlcRose1qwh+f2yt4rvz+5R8fTxgPGNcPGOXua4FLKniNiEht7NATamaHAr8FLnD3QkLv8DfA82b2/dgFiofEeqinEGYPGWVmRwBnu/s6Qq691Mx+a2aHEToY4sPoEvNfZTk5LrH3eB+gs5ntBBwEtDSzvrHnthIu9ot3sCSe4z7gRDNrH3u8P/CAuxcTap4Kc7GZDY/F1bySn10O5Wqm2H8+jiN8S9meMMRwD2BXYKeECyPL94pvqSD+skFVMF4k9lm0BV4l9MrX5j8c0kA0FEQywa7Af83sOcLXZAZcmPB8MfAbM1tN6N29BMDdXzazSwlfJ44iXPz3CYCZnQD8gfCH4k/ABHd3MzuZMBXSf4EnCF/pDSEsqHAyYYhJmYVpzGwY4aKcXoSvBv8O/DYWzz+B4UCeme1G6EnpCpxsZp+xvffmHKCVuxeUe+9/trCwQzHhYiERkVqL5ZHTCcXfA0Ah4QK4nYHvuft7AO6+JZYfHyDkumcJKz20xwAAIABJREFUY6sBfkbolHiY8A1dfCaOPxJy30WEIWsXufvyWOE4HNjNzB6tLCeb2VjCuOLzzOwzd19MuJDyVkJR/TBwArA38BXwCKF3/SwLs0XtBRxrZs+7+1/NrB/wjJm9Q/jW8zexnvWTgK6x9/dfwmI53cxsr1icpxLGjf894efWlJD7BwM9YoW3E8Y2j4kd82PgX8CHhOEtLxBy+jTCBZODzWyUu8eP+wRhsZ71sUK8KzDWzOYReqKJvb78mg1fEy7OHAW0AHLM7BXCBYwNNQWtVMP0WUgms7AM7P+5e14DnW+Cu2u5cRGRLBMbltLc3dNygS4zm0CYoWV+7HELwrcNN8evzZHoaSiIZLqKhlmk5kThK9OZDXEuERFpOLGLM/dN46K6BXAVkJfQ3AnYoKI6vWgoiGSs2Pi2HxK+2jvJ3Z9K8Sm/jH1FKSIi2aWJu98bdRCVcfciC8uq3xW7sH0JYdjKVVW/UhqahoKIiIiIiCSBhoKIiIiIiCSBCmsRERERkSRQYS0iIiIikgRZc/Fip06dPC8vr9avKywspHPnzskPqJHTzzU19HNNnSh/tjNmzFjh7o3qg1XO3pHeW2bSe8tM9XlvVeXsrCms8/LymD699iuJ5ufn1+l1UjX9XFNDP9fUifJna/b/27vzODnqOv/jr0/PfWVyDblDQkIOuTUcisghLCoq6AJerAcLwft2VRZdRVZR0R8rKhoVXS/WXS8UFRUPDlFMAoQzgZDL3JN7JnPPfH5/fKtnejpzT/V09/T7+Xh8H1VdVV31rerOZz759re+ZZuzcuAsUsw+ks4tP+nc8tNozm2gmK2uICIiIiIiMVBiLSIiIiISg4JPrJcvX57tKoxLuq6ZoeuaObq2+WE8f046t/ykc8tPmTq3MX9AjJldANwAvCblefevBr4ATAC+D3zA3TuidVOB/wQOACXAh929PX2/y5Yt8/HaD0hExjczW+3uy7Jdj74oZouI9DZQzB7TFmszmwZUA6elLJsLXAJcCrwPuAp4b8rb/g+4zd0/DDwNXD9mFRYRKWCK2SIiwzOmibW77wLuSFt8NHClu6929+8CXwHOBTCzM4DF7v5gtO0vgXeZWdVY1VlEpFApZouIDM+Y97F296601/clf0KMbAO2RPPnARtTtt0Wzcbzk+knPgEnnwy//GUsuxMRGW9yJWY3PNzAquetYu1b1o52VyIiGZOL41ifCnw+mp8F7E1b3wDMTH9TfX09y5b1xO7ly5cP3jF92zZYswZ27hxNfUVEhm3FihWsWLEi+XJqNusySmMSs7uau2h8qJFEWcHfcy8iWTDUmJ1TibWZLQTq3X1NtMiB5rTNyoAjboSpq6sb/kDflZVhevjwMGsqIjI6qYmkme3JcnVGZCxjthVbOEDH2N5wLyICQ4/ZOfNffzMrAq4BPpKyeBswKWUbA6qA7bEctCrq9qfEWkRkWMY6ZiuxFpF8kDOJNfBB4CZ3b01Z9ntgccrreUAbEM8YTcnEuqkplt2JiBSQMY3ZSqxFJB+MeWIdtWAAWMqya4GHgQozO8bMrjSzhe6+CthhZsdHm74EuMXd22KpjLqCiIgMKFdithJrEckHY9rH2sxqgCuil282sy8Tfkr8VNqma939tmj+NcDHzWwTYTzVf4+tQuoKIiLSr1yK2UqsRSQfjGli7e4NwK1RSbohKv29ZzPwrxmpULLFWl1BRESOkEsxW4m1iOSDXOpjPfbUYi0ikheUWItIPlBiDUqsRURynBJrEckHhZ1YqyuIiEhe6E6s25VYi0juKuzEuro6TBsaslsPEREZkFqsRSQfFHZiPSl6jsH+/dmth4iIDEiJtYjkg8JOrKdMCdO9e7NbDxERGZCVKLEWkdxX2Il1VRWUlkJzcygiIpKT1GItIvmgsBNrM5g8Oczv25fduoiISL+sqCexdldyLSK5qbATa1B3EBGRPGAJ6/mL1ZXVqoiI9EuJtVqsRUTygrqDiEiuU2J91FFhunNndushIiIDUmItIrlOifXs2WG6dWt26yEiIgNSYi0iuU6J9Zw5YarEWkQkpyWH3OtqUydrEclNSqzVYi0ikhcS5eFPVlezEmsRyU1KrJVYi4jkhaKKIgC6WpRYi0huUmKd7AqyeXN26yEiIgNKVKjFWkRymxLrWbOgvDyMCnLwYLZrIyIi/Ugm1p3NnVmuiYhI35RYFxXB4sVhfu3a7NZFRET6pRZrEcl1SqwBli4N06eeym49RESkX919rJVYi0iOUmINSqxFRPKAWqxFJNcpsQZ4znPC9NFHs1sPERHpl/pYi0iuU2INcOqpYbpyJbie6CUikovUYi0iuW7ME2szu8DMHjSzeSnLKszsFjP7tJndamYThrIuNnPnwrRpsHcvPPts7LsXEclXuRSz1cdaRHLdmCbWZjYNqAZOS1v1VeB+d78WuAP4+hDXxVUxOP30MP/gg7HvXkQkH+VazC6qDol1Z4O6gohIbhrTxNrddxECbTczmwm8HvhltOhu4JVmNm+gdbFX7owzwvSee2LftYhIPsq1mF08sRiAjoMdcexORCR2Y94VxN3Tf8M7B9jj7k3R+g5gI3D2IOvidcEFYfq736mftYhIJJdidndifUCJtYjkply4eXEWsDdtWQMwc5B1vdTX17Ns2bLusmLFiuHV4pRTYMqU8Gjz9euH914RkWFasWJFd7wCpma7PsOQtZitxFpEsmWoMbt47KrULwea05aVAe3R/EDrutXV1bFq1aqR16KoCM4/H370o9BqfeyxI9+XiMggli9fzvLlywEwsz1Zrs5wZC1mK7EWkWwZaszOhRbrbcCktGU1wPZB1sXvwgvD9I47Bt5ORKRwZS1mK7EWkVyXC4n1n4CZZlYBYGYlwFzgnkHWxe/ii6G4GP74R6ivz8ghRETyXNZidndivV+JtYjkpmyMY23JWQB33wncCbw4Wn4+8BN33zbQuoxUbvLkcBNjZyf85CcZOYSISD7JpZidTKzb9x/Rs0REJCeM9TjWNcBbo5dvNrO6aP5twMVm9u/Ay4DlKW8baF38XvvaMP3udzN6GBGRXJdrMbt4SjEYdOztoKtdD4kRkdxjPk6Gllu2bJmP6ubFpMOHYeZMOHQIHn4YTj559PsUERmAma1292XZrsdYGmnM/sv0v9C+q50z/nEG5bPLM1AzEZGBDRSzc6GPdW6pqoI3vznMf/WrWa2KiIj0VjajDIC2HW1ZromIyJGUWPfl7W8P0+99D3buzG5dRESkW+mMUkCJtYjkJiXWfVm8GC65BFpa4POfz3ZtREQkosRaRHKZEuv+fPzjYXrrrWq1FhHJEcnEunV7a5ZrIiJyJCXW/TnllDCudXNzT5ItIiJZVX50uGGxZVNLlmsiInIkJdYDufFGKCmBb34TVq7Mdm1ERApexbEVADSvT39yuohI9imxHsiSJfC+94E7vO1t0KGnfYmIZFPFwiixfkaJtYjkHiXWg7nuOpgzB1avhk9/Otu1EREpaGUzy0iUJ2ivb6fjoBo7RCS3KLEeTE0NfPvbYf5Tn4I4HkIjIiIjYgnrabVWdxARyTFKrIfixS+G97wndAV5wxvCUxlFRCQrKpdWAtD4WGOWayIi0psS66H6zGfghBPg6afhTW+Crq5s10hEpCBVP7cagMaHlViLSG5RYj1UFRXw059CbS38/Ofqby0ikiU1z60BoPEhJdYikluUWA/HwoXwgx+AGXzsYz19r0VEZMxUnxK1WD/SiHd5lmsjItJDifVwXXQR3HxzmL/qqtB6LSIiY6a0rpSyOWV0NnZy+MnD2a6OiEg3JdYj8e53hxbrri64/HL48Y+zXSMRkYJS+6JaAA78+UCWayIi0kOJ9Uh98pPw/vdDezu85jXqFiIiMoYmnTsJgAN/UmItIrkjlsTazF4Wx37yihncdFNIsLu64Mor4ROf0GghIpLzxkPMnnjORCC0WKuftYjkimEl1mb2EjO738zWm9mGqGwCfpKZ6uU4M/j4x+FLX4JEIiTZl10GjbpTXUSybzzH7PJjyimfV07Hvg4OPahnC4hIbige5vbfAz4DPAwkm2YNeH2clco773pXGDHkda8LQ/I99RTcfjucdFK2ayYihW3cxmwzY+olU9l681b2/HQPtc+vzXaVRESGnVhvcPcvpi80s7Ux1Sd/vfSl8OCDcMklIbE+7TS48cbwxMaEurKLSFaM65g99dUhsa7/aT3HfO4YzCzbVRKRAjfcjO9tZnZpH8tfGUdl8t7ixbB6NVxzDbS1hZsbzzsP1o6Lv2Eikn/GdcyufUEtJdNKaNnQQsPKhmxXR0Rk2In1mcA3zawzpXQBt462Ima21My+YmbvM7OvmtnJ0fIKM7vFzD5tZrea2YTRHiujKivha18L41vX1cE998CJJ4bh+Zqbs107ESks4zpmW5Ex7YppAOz4xo5MHUZEZMiGm1hfDywHjgWOSSmfiKEu3wM+7e7/D7gRuD1a/lXgfne/FrgD+HoMx8q8iy8OLdVXXRWG5LvhBli6FL7/fY0cIiJjZdzH7BlXzQBg1+276GjoyOShREQGNdzE+jfAHe6+wd03R2UT8I0Y6rIUqInmm4FaM5tJuMnml9Hyu4FXmtm8GI6XeZMnwze+AffdB8cfD5s3w7/8CzzvefDb34JriCgRyahxH7OrllRR+6Jaug53sfO2nZk6jIjIkAw3sf4f4KNm9qKUcjbw0Rjq8mPCT5Y1wBXAu4BzgD3u3gTg7h3ARuDsGI43dl74QnjkEbjtNpg9O8y/5CVw+umhy4hasEUkMwoiZs95/xwAtnxuC12tiqcikj3DTaw/BrwX+C7w3ynlmhjq8nagBVgJNLr7T4BZwN607RqAmTEcb2wVFcFb3gJPPw2f/SxMnQorV8KrXgUnnBCe3Kg+2CISr4KI2VNeOYWqE6to297Gjm+rr7WIZM9wE+vPA3PcfZ67z4/KPODVMdSlgtBH74fAzWZ2PuCEnxhTlQHt6W+ur69n2bJl3WXFihUxVCkDKirg3/4tdAv50pdgzhx48snw5MZZs+CDH4T167NdSxHJsBUrVnTHK2Bqhg5TEDHbzDj6uqMB2PzJzXQcUl9rEYnXUGO2+TD6+ZrZ8939r30sv9LdbxtRTXv2cRfwWnc/YGYfA95NaGn5D3dflLLdM9GyH6a+f9myZb5q1arRVCE72tvDw2S+/OXQgp10wQXwpjeFcbGrqrJXPxHJODNb7e7LMrDfgonZ3uU8/MKHOfTXQ8z54BwWfH5BbPsWEUk1UMwebov12/vY+WXAzSOpWMo+pgInufuBaNENwCFgMzDTzCqi7UqAucA9ozleTikpgTe+Ef7+95BYX3kllJfD738PV1wB06aFGx5/9zvoUCuMiAxLwcRsSxgLv7QQDLbevJWGRzSutYiMveEm1pOjlgnMrNLMvgl8k9DHbjT2AS1mNitl2V7gUeBO4MXRsvOBn7j7tlEeLzctWwbf+hZs3w633goveAEcPhyG6LvwwtBVZPly+PWvobU127UVkdxXUDF7wrIJzHz7TLzDeeqKp+hs6RyLw4qIdBtWVxAAMzsDeANwIbATeCNw0N33j6oiZicRWldWA9OAe939HjObBHwO2ARMBz7q7o3p78/briCDefZZ+MEPQnL9zDM9y2tq4GUvCzc/XnghTJyYvTqKyKhkqitItO+CitmdTZ2sOmUVzU83M/PtM1n0lUWDv0lEZBgGitkDJtZm9qJ+Vp0GnAvcBHQBl7j7+0Zb0dEYt4l1kjs8+mgYnu9nP4M1a3rWJRJh6L6XvCQk2cuWhVFIRCQvxJVYK2YHh1Yd4uEzH8bbnEUrFjHz6vwbSEpEctdoEuu/ER4CcIAQjPvcDJjh7mWjrehojPvEOt2GDSHJvuMOeOCB3v2vJ0+G88+H886Ds8+GxYvBLHt1FZEBxZhYK2ZHdnx7B+uuXIeVGCfedSKTzpuUsWOJSGEZTWL9MmCVu+8e5ACXuvuPR1fN0Sm4xDrVoUPwpz+FpznedRds3Nh7/bRp8KIXhST77LPhOc8JrdwikhNiTKwVs1Os/8B6tn5xK4mqBCf97iRqX1Cb0eOJSGEYcWKdTwo6sU7lHsbB/t3v4J57Qtmd9jd2yhQ46yw444xQli3TkH4iWZTJPta5aixitnc5a9+yll3f3UXRhCJO/PWJ1J6p5FpERmegmF081pWRDDODY48N5R3vCIn2unU9SfY994RRR37+81AgtF6fcEJIsk8/PUwXL1artojkNUsYi7+1mK6WLur/t541569h6e1LqbukLttVE5FxSi3WhcY9jDTywAPwt7/Bgw+GGyE704alqq2FU07pXZYsgWL9X0wkbmqxzqyuji6eeccz7FixAxKw4HMLmP3+2ZjuPRGREVCLtfQwg4ULQ3njG8OypiZYvTok2clke+tW+POfQ0kqLw8t2yefHBLtk0+G446DCROycSYiIkOSKE6w6GuLKJtdxqaPb+LZDz7Lwb8cZMm3l1Bcqz+DIhKfUUcUMzsF+ADwMOFBAJtGu08ZY5WVoc/1WWf1LNu+HR5+uHfZuDE8HXJl2rMl5swJCXayHH88LF0K1dVjex4iMqhCjdlmxryPzaPqhCrWvmkte362h1VrVrHk20uY+CI9B0BE4hHHf9XfCpzv7leY2evNbJe7N8ewX8mmmTNDueiinmUHDsAjj/Qk2o89Bk89Bf/4Ryh33dV7H/Pm9STbS5aEftuLFsHUqWN6KiLSS0HH7LpL6qh+qJonLn2CxkcaeeScR5j1rlkc8+ljKKrS+P8iMjqj7mNtZiXAUe6+zcwSQLW7H4qldsOgPtZZ0tERxtR+4olQHn88TNetg/b2vt8zeXJIspNl0aIwXbgQyrI6tK5IVoxlH2vF7KCrrYvNN2xm86c3QyeUzS7jmM8fw1GvOUp9r0VkQLEOt2dmLwUmu/sPzOw0oNHdn4yhnqOS7SAtadrbw7B/jz8OTz4ZEu116+Dpp6HxiKcbB4lEaOVOJtrHHgsLFoRy9NFQWjqmpyAyVjL8SHPF7AE0PNTAuqvX0fhQiEu1L6xlwRcWMOE03TsiIn2LLbE2s/8E3gPc5e6XRstuBB519x/GUdmRypUgLYNwhx07eifayfmNG6Grn4fFJRIwd25Pop1eamrG9jxEYpSpxFoxe2i809nx7R1svHYj7fXhl7bJF01m/ifnU/M8xRYR6S3OUUFeAMwG3pSy7DfA94CsBmnJE2Y9/bfPPbf3utbW0K0kmWivXx+GBnz22dCHe9OmUP7whyP3W1d3ZLI9b15o6Z41S8MESqFSzB4CKzJmXjWToy47ii03bmHrLVvZ96t97PvVPqa8cgpzPjSH2jNr1UVERAY13Gzj9+5+wMxSm7kvjbNCUsDKysJoIkuXHrmutTUk1clEO7Vs2AD19aH87W9HvreoCGbPDkl2MtlOnZ8zR327ZbxSzB6G4tpijvnMMcx+/2z+8fl/sO3L29j7i73s/cVeapbVMPt9s6m7rI5EiR6eJSJ9G25XkNcCxwCzgPuAy4BXAe9y969kpIZDlEs/K8oY6+oKwwOmJ9ubN4dkfMeOgd9vBjNmHJlwz54dypw54YZLtVZJhmSwK4hi9ii07Wpj6y1b2f617XTs7QCgdGYpM66cwfQrp1MxvyLLNRSRbIj75sUzgDcDc4EdwO3ufvdoKzla+RCkJUtaW3u6kmze3JNwJ+e3bj3yyZPpyst7kuxkwp3+eupUJd8yIhm+eVExe5Q6mzvZ9f1dbL15K01PNnUvn3jeRGb86wymvmoqRRUaqk+kUMSaWPex8yLgHHfvo+Pr2MmnIC05pqMDtm3rnXRv2RIS7q1bQ1J+8ODg+ykr6zvpTvYpnzEDpk/X6CZyhDEebk8xe4TcnQP3HGDnt3ZS/+N6ulrCzdZFE4qYevFU6i6vY/I/TSZRqq4iIuPZiBNrM/sFMNjj86YCe9z9vJFXcfTyMUhLHmloCMl3MtFOTbqT8/v3D21fdXUhyU5NuNPnp0+HkpLMnpPkjLgSa8XssdNxsINdt+9i5207aVjZ0L28eGIxUy+ZSt1ldUw8byJF5WrJFhlvRjMqyKPRNuuALuBqwmNwU6PhfGByDPUUyV01NeHpkUuW9L9NY2NP8p2adG/fHvp5b98OO3f23Gj56KMDH7OubuDke8YMOOoo3XgpqRSzx0hxbTGz3jqLWW+dRdMzTdT/Xz27f7Sbw48eZud3drLzOztJVCaYdMEkpr5iKpNfNpmyGfq3KjLeDdZiXRttcyB6fb27f7yP7e5w94szV83B5XvrhxSIzs6QVG/f3jvhTn+9a1f/Y3qnmzQJpk0LrdzTp/fMpy876igNO5ijYmyxVszOssNrD1P/v/Xs+fkeGh/u/TCsmlNrmPzSyUx68SQmnDFBXUZE8tSIW6zdPb1j6TQzK3L37ju9zOwc4LRR11KkEBQV9SS7z31u/9t1dsLu3X0n4Dt2hJbxXbtC2b8/lLVrBz62WbjBsr/EO3U6dWp4KI/kFcXs7KtaUkXVx6uY9/F5tGxtYd+v9rHnl3s48IcDNKxsoGFlA5uv30yiMkHtWbVMOn8Sk148ieqTqrGEbn4WyXfDbb76MbDGzO4HGoClwIXAES0iIjIKRUWhq8eMGfC85/W/XVcX7N0bEuydO4+cps4nu6DU14dHzQ92/Lq6wVvBp08PLeYaDSVXKWZnUfnscmZeM5OZ18yks6mT/X/cz/7f72f/H/bT9EQT+3+7n/2/DfdmFE8pZuJZE5lw5gRqX1hLzXNr1KItkodGMtzeTODtwInAIeD/3P2ODNRtWMbrz4oisenogD17+k7A0xPxvXuHvt+Skr4T7r7ma2qUhPchw8PtKWbnoNadrRz44wH23x0S7dYtrb3WJ8oT1JxWQ+0La6k9s5YJz59AySTd0CySCzI93F4COMXdV49qRz37OwN4IbAeuB84DHyO0NoyCfiwux9Kf1+hB2mRWLW1hZbt/lq/U5cNZSjCpIqKwZPv5LSyMnPnl2PGeLg9xewc4+60bGjh4P0HOfiXgxy8/yBNTzUdsV3FsRXUnFoTyrIaak6poahKo46IjLXRjAqSvqPb+lhcB3QCl4ygbun7fwuwwN2vS1n2beAud/+Rmb0E+DrwutEeS0QGUFoKs2aFMpiWloG7oCTnd+yApibYuDGUwdTUDK0VfNo0jQ3eD8Xs/GBmVCyooGJBBdPfNB2A9r3tHHwgJNkH7z9Iw+oGmp9ppvmZZnb/cHd4YwKqjqvqSbafV0PV8VV6WI1IFg33keb3AX8AUt90JnCfu39qVBUxOxv4GHCBR5WKfsLcCExy9yYzKwYOAse5+6bU96v1QyQPNDb2n3ynv25rG/p+J0/u3do9Y0bPg3pmzQrTGTNyNgHP4CPNFbPHia72Lg4/frj7BshDKw9x+PHD4b9IqRJQuaiS6pOrqTqpiuqTqqk+uZrS6aWYumGJxCK2FmvgKndfl7bzJYTH5Y7WF4EHgFvMbCHwKeBowoMMmgDcvcPMNgJnA5tiOKaIjKXq6lAWLBh4O/fQxaS/1u/U+V27YN++UJ56qv99moWkO5lopyfeyVJREe85Z5di9jiRKElQc0ro/sHysKyzqZPGRxq7E+3GhxtpWtdE09pQ+J+e95fUlXQn2VUnVlH1nCoql1ZSVKnWbZE4xdHHej6w0t2njmIfi4GngJPd/VEz+xDhrvUbgDe4+4kp2/4V+IW7fyZ1H0cffbTX1dV1v16+fDnLly8faZVEJF+kj4yyc2cYlnDbtt4P7NmxY2hjg8+cCccc07ssWBCm06bFdvPlihUrWLFiBQCrV6/e7O7zYtnxIBSzx7fOlk6anmiicU0jjY80humaRjoPpjdtAwbl88tDkn1cJVXPqaLquCoql1Sq77ZImqHG7OF2BdlI758Ui4BpwPfc/eqRVtbMXg2sSAZ6M6sA6oHrgX9299NTtn0I+KG735S6D/2sKCID6ugIyXcy0d66tXfinXxaZkdH//uYOBGOPz6UE06Ak0+GZctG3cUkg11BFLMFd6d1S2t3on348cMcfuIwzU834x195AAG5fPKQ5L9nMowXVpJ5aJKimv1kCmROLuC/B74IT2BugvY5e5Pj6J+yXp018Xdm83sGcIfgUlp29YA20d5PBEpNMXFPTdknn5639t0doYEe8OGUJ59tmd+/frwIJ777w8lqaICzjwTXvUqeOMbQ1eX3KGYLZgZ5UeXU350OVMv7vmhoqu9i+Znmjn85GGanmji8JM9CXfLxhZaNraw987eQ2+WTCuhcnFld6lYXEHlokrK55eTKNG42yLDTay/7+73pi4ws8VmttDd14+iHo8CtWY2zd13Rcs6CMF4pplVRIG7BJgL3DOKY4mI9K2oCI4+OpRzz+29zj20eD/2WHjAzmOPwd//Dk88AXffHcqXvgQPPZRLQwUqZku/EiWJ0P3jOVVwac/yrvYumtc30/RkE4efCMl209ommp9ppn1XOwd3HeTgvb2H2bRio3xBee+ke1EFlYsrKakr0Y2TUjCGm1i/BLg3bdku4E/AKSOthLuvNbPfEP5pf8XMJgKlwA+AlwIvBu4Ezgd+4u7bRnosEZERMQujjtTVwZIl4YmYp54Kv/0t3BE9b2XdutC/e+HC7Na1h2K2DFuiJEHV0iqqllZR9889/eC9y2nd2hpukFzXRPO65u751i2tNK9rpnldM3vp3cpdPLE4tGwvrqTi2Aoqjw3TimMrKK5R1xIZX4b0jTazdxJuSqkxsw/3scnvY6jLG4H/MrNSQgvH66I7yt8GfM7MTgKm030/tIhITLq6wigkBw6E7h7J0UbSH4qzdSts3tz3UIBTpsC11+ZEUq2YLZlgCaN8bjnlc8uZfMHkXus6mztpfqa5z6S740AHDQ820PBgwxH7LJlW0ivR7k68F1boBkrJS0NaQ3DEAAAgAElEQVS+edHMTgTeBHwpbVWbu++Iu2LDpRthRKSX9ethxYowfnVTU0iak4lzsiRfHzwYunoM1bRpMH8+LFoEp5wSWq5PPz304x6BTNy8qJgtucDdad/d3pNwRw+5aXqmiZZnW+hq6X+kntIZpb2T7WTyvbBCD8GRrIrl5sVoSKWPuHt7HweY6+5bRlNJEZHYuIcROw4fHvp7JkyASZNCqas78hHr06eHJH3ePKiqyljV46KYLbnAzCidVkrptFImvmhir3XJriXJRDuZdDc/00zzhmbadrTRtqPtiP7cAGWzy45s5T62gvJjyikqV9It2TNgYm1mHwR+7e5PRov+zczSv7EGnAucE3/1RERG6Lrr4KMf7Xvdj34UxqVOJtITJoy4tTmXKGZLPkntWjLpxb0Hk/FOp2VLS68W7uR8y8YWWre20rq1lQN/OpC2UyibW9Zn95Ly+eUkSjVyiWTWYH9J5gOp/8V8LlAF7KRn+CYj9KMTEckNZvCRj4Syezfcdx985Svwpz+F9Y8+Cpdfnt06ZoZitowLVmRUzK+gYn4F/FPvdV0dXbRsaunVwp1MvFs2tdC6uZXWza3sv3t/7zcmwvjcFQvTupcsqKB8XjmJMiXdMnoDJtbu/o60Rf8JPO7uve7cMbPnxF0xEZFha23tedDLli2wdm0YDu/xx8NY1ElHH529OmaQYrYUgkRxgsqFlVQurAxj0KToauuiZWPLEV1Lmp4JI5e0bGihZUML+3+XlnQblM0po+KYCsoXlFOxoKK7lC8op2RiydidoOS14f72eYK7P5S6IHq07QLgyb7fIiIyAu7hpsODB0PZtw/q63uXPXt65rdtCyN39KemBs45B5Yvh5e/fMxOI8sUs6WgJEoT3eNop+ts6aRlQ0vvriXrm0OyvaWF1i2ttG5phT8fud/iycX9Jt1lM8uwhMbplmCow+2dRfj58ILoEbmpJhPuOv9lzHUTkXziHlqMDx+GxsYwTS3pyxobe5Lm/kpn5/DqUFQEM2fC3LkwZw4ceywcd1woS5aMi37UQ6GYLXKkovKingfipOlq66JlcwvNzzbT8myYNm/ome/Y10HDvgYaVh05ZKCVhW4rfSXd5fN0M2WhGepfmQPAdwn9985MW9cGfDPOSolIDLq6QqLb3AwtLWGaOj/cafqyvhLmrv6HzhqR8nKorQ0lOVpHepk6NUxnzgwjdxRI8jwIxWyRYUiUJqg8tpLKY49s6XZ32na2haR7Q5R0pyTg7fXtNK1tomlt05E7NiibVXZk0n1MeF0yWV1Mxpsh/QVy98fM7BzgMndXQBbp7IT29iNLW9uRy1pbw/LhTOPYtv2IUdYyr6QEqqvDcHSppb9lyaS5v1JaOvbnMA4oZovEx8wom1FG2YwyeOGR6zsaOnol3N1J94ZmWjb3jGBy8J4jhw0snljcd9K9sIKy2WV6FHweGs441gfpp5XDzF7o7vfHViuRZ58NIzi0tYVW0KGWzs6e+f4S3YGS4KGuG87DRLKprAwqKkLLb+q0r2XDnaYnylVVIbGWnKCYLTI2imuKqT6pmuqTqo9Y19XeFR73nky40xLwjgMdNK5upHF14xHvLaotouaUGqqfW03Nc2uofl41lYsq1Z87xw02jvUvgCO/KWmbAccBR8VVKSlQ27bBr34Ft98Of/5ztmszMLOQRPZVSkt7vy4rC8sGmg5lm+G+p7Q01FMKhmK2SG5JlCS6W6LTJZ9K2byhd9eS5mfDSCbt9e0c+PMBDvy5Z6zuoglF1CyrYcJpE6g5tYaa02oom6WW7VwyWIv1OkJfvW30jIGazoCCucVeYrRjBzzwAPzlLyGRfvjhnnUVFWHkhilTIJEIN6UlEsMrfSW5AyXAQ11XUhLqo0AmuUcxWyRPpD6Vsvb5tUesb93RSuNDjTQ81EDD6gYaVzeGh+L88QAH/tiTbJfOKA1J9qk13S3nZXOUbGfLYIn1Z4EGd28daCMz093l0r/OTli/Hh57LDyY49FHYc0a2LSp93aVlXDBBfCKV8Cll4Y+tiIyHIrZIuNE2Ywyyi4qY8pFU7qXtW5vpWFlA4dWHqLh7w00rGygbUcbe3+xl72/2Nu9XfHEYqpOrOpOtKtOrKLq+CqKKjRCSaYN9oCYPUPcz6vQXeaFrbExPJhj48aQRKeWjRv7vpGuuhrOOAPOPDOUs84K/XdFZEQUs0XGt7KZZZRdXMbUi6cC4F1O8/pmGlaGoQAbH22kcU0jHXs7OHjvQQ7em3LDZALK55dTuSSM8506LakrUQt3TEY1LpWZzQXeAbwTBenxxz2MJZz6QI5du0Jf6K1bQ0nOHzzybude5syBE08M5YQTQimgcYVFcoFitsj4YgmjclEllYsqmfaGaUA0POCONhrXNHL40cM0rmmk8dFGmtY20fJsCy3PtrDvV/t67ad4UjGViyu7x94un1cexuaeV07ZnDISpXrc+1ANO6sxsyLglcDVwAXALkKfvvFn+3a49Va4++4wSkVDQ+i7O3FiGFM3WaZMCWPpJsfTTc5PnQqTJ4eW2bKyse+T694z7nBjIxw40P+DOJLrkk+y2707zA91yLbycpg1C+bNg4ULe8qCBXDMMWHECBEZcwUVs0UkDA84s4yymWVMeWlPN5Ku1i6a1zfTtC6Mud09XdtEx/4ODv3tEIf+dqiPHYaxuMuODkMOlk4vPaKUTCuhuLaYouqiWFq+3Z2u5i7a97XTvqed9vq+p231bd2vu5q6KJ1eSsWiCiadP4lpV0yjtG7sh2wdcmJtZvOBq4C3AFOAzcArgN8Cz89I7bLpZz+DN74xJKTpmppC0j0cRUUhwa6pCdPq6tCnuLi4dykq6j3f2Tlw6ejoeXhHU1PPfLKMVnV174dxHHUUzJ4dyqxZPfOTJ+tmPpEcUnAxW0QGlChLUHVcFVXH9W7ocnfadrXRvK6Z5o3NtGxqoWVjS5hu6hmHu3XrgLduRAeB4gnFFNUWUVxbTHFtMYnyBFZiWLFhJUaiJAEJ8Hanq60LbwvTrpYuOg500LGvg/Z97Xjr8Ie1bV4fHlO/79f72HjtRuZ+ZC5zr50bjjlGBk2szewyQkvHecB64IvAfwPXuvtd0WYPZKyG2XD//XD55SFpffnL4R3vgJNOCjfTdXbC/v29y759PS29e/b0nt+3LyTnHR09rcNjKTmO8UAP45g4sWd+ypSQPCcTafV5FskrBRmzRWTEzIyy6WWUTS9j4tkTj1jf1d5F69ZWWja30Lazrc/SvqudjoMddDVHyfGBDloZQiI+UL1KjZLJJZTUlVAy9chpaV1pr9eJygRtO9poWN3A7v/Zzb5f7WPTJzZx4L4DnHDHCRRVjc2Nm0NpsW4mPAL3AeAqd38awMzy5AkZw9TeDsuXh0T43e+Gm28+siW2pgbmzh3eftvaQoKdLMlHQCdbnVNLcllnZ89Qc/2V4uKeB35UVvbMJ0tC/aJECkxhxWwRyahESYKK+RVUzD9yLO50Xe1ddBzsoPNgJx0HO+g41IG3Ol3tXXiH4+1R6XISpQms1HqmZQmKJxZTMrmE4knFJCoSw+5WUjKxhKqlVUy/Yjr7/7yfJ1/7JAf+cIAnLn+CE355wpg8XGfQxNrd7wTujG56udrM5gF3AN2PWDMzc8+XR9EN4qc/haeeCn2DP/e5+Lo3lJaG7hKTJ8ezPxGRPhRczBaRnJEoSVA6tRSmZrsmMOmcSZxy7yk89PyH2Pfrffzjpn8w99+G2Sg6AkNuznT3Le7+MUJ/PYDZZvapKGhfnYG6Zcett4bp+98fulGIiOShgonZIiL9qFxUydLvLQVg0yc30fKPlowfc9j9BNy9w91/7O6vIvTbex9wc+w1y4bdu+Hee0Pr8hVXZLs2IiKjNq5jtojIIKa8bAp1l9bR1dTF5us3Z/x4o+qA6+7r3f09wFtjqk923XlnGKLuvPNgwoRs10ZEJFbjLmaLiAzB/P+cDwY7v7uT1p2ju6lyMLHc2ebu341jP2ZWamZrzOyc6HWFmd1iZp82s1vNLLPZ7u9+F6aveEVGDyMikk1xxWzIgbgtIjKIykWVTL14Kt7m7PzWzoweK9eGjPgQMC/l9VeB+939WsLNN1/P2JHd4Z57wvy552bsMCIi40z24raIyBDNuHoGALt+sItM3rudM4m1mb0A2AHsj17PBF4P/DLa5G7gldGNN/Fbvx527gxjNy9ZkpFDiIiMJ1mP2yIiQzTpgkmUTC2h6akmGtf08fC/mOREYm1mVcCl7n5byuJzgD3u3gThBhxgI3B2Ripx771h+qIX6QmCIiKDyIm4LSIyRImSBFNfFcYB3PfrfZk7Tsb2PDwfBm5MWzYL2Ju2rAGY2dcO6uvrWbZsWXdZsWLF8Grwt7+F6ZlnDu99IiIjsGLFiu54RU6M+jpso4rbo47ZIiLDNPnC8CyR/b/fP+z3DjVmD+XJixllZhcCq9x9d9oqJzxBLFUZ0N7Xfurq6li1atXIK7JyZZieeurI9yEiMkTLly9n+fLlAJjZnixXZ1jiiNujjtkiIsM08byJkICDfzlIR2MHxdVDT4OHGrNzocX6g8BtZrYnqugcwg0vVcCktG1rgO2x16C5GR5/PDz++5RTYt+9iMg4k/24LSIyTCWTSqg+uRpvdxpXZ6afdS4k1q8HTk4p24GrgG8DM82sAsDMSoC5wD2x1+CRR6CzE447DqqqYt+9iMg4k/24LSIyAjWn1gDQsKohI/vPemLt7vXuvjVZgE6g3t23AHcCL442PR/4ibtvi70SyW4god+MiIgMICfitojICEw4NQytf2jloYzsP+t9rAfxNuBzZnYSMB1YnpGjJPv5qX+1iMhojU3cFhEZgZplocW68aHMdAXJucTa3eelzO8Hrs74QXXjoojIiGUlbouIjEDF4gowaN7QTFdbF4nSeDtvZL0rSNY1NsK6dVBSAieckO3aiIiIiEiGFJUXUT6vHDqhZWNL7PtXYv344+Fx5s95DpSVZbs2IiIiIpJBFYsqAGh6uin2fSuxfuyxMD3++OzWQ0REREQyrnJRJQDNT6cPuz96SqyTibW6gYiIiIiMe2VzQg+F1m2tse9bibUSaxEREZGCUTqjFIC2HW2x77uwE2t3JdYiIiIiBaRsRtRivUMt1vHatQv27oXaWpg9O9u1EREREZEM626x3q4W63g980yYLlkCZtmti4iIiIhkXOn0KLHercQ6Xhs3hun8+dmth4iIiIiMiaIJRQB0HurE3WPdtxJrUGItIiIiUiASxQkSlQlw6DzcGe++Y91bvlFiLSIiIlJwiicUA6HVOk6FnVhv2hSm8+ZlsxYiIiIiMoaKakN3kI5DHbHut7AT6927w3TGjOzWQ0RERETGTHeL9UG1WMdnz54wravLbj1EREREZMwkykMK3NXaFe9+Y91bPunqCmNYA0yenN26iIiIiMiYsdIwzLK3a1SQeOzfH5LriROhpCTbtRERERGRMZIoiVqs29RiHY99+8JUrdUiIiIiBaW7xbpNLdbxaGoK0+rq7NZDRERERMZUolQt1vFKJtYVFdmth4iIiIiMKbVYxy2ZWFdWZrceIiIiIjKmulusNSpITJRYi4iIiBSm8HwYvHOctlib2avNbKOZ7TWz/zKz4mj5VDP7upl91sy+aGbxDOHR3Bym6goiIjJsYx6zRURiZBa6ghBvXp0bibWZzQUuAS4F3gdcBbw3Wv1/wG3u/mHgaeD6WA7a1hamZWWx7E5EpFBkJWaLiMQpyqvjTqyL493diB0NXOnuHcBqMzseONfM7gcWu/uD0Xa/BNaZ2Q3ufnhUR+yMHmGZyIn/W4iI5JOxj9kiInFKNlj7OOwK4u73RQE6aRuwBTgP2Jiy3bZodtmoD5pMrIuKRr0rEZFCkpWYLSISpwy1WOdEYt2HU4GvAbOAvWnrGoCZoz5CV3QXqBJrEZHRynzMFhGJUab6WOdKV5BuZrYQqHf3NWbmQHPaJmVAe/r76uvrWbasp1Fk+fLlLF++vP8DqSuIiGTRihUrWLFiRfLl1GzWZTTGLGaLiMRpmC3WQ43ZOZVYm1kRcA3wkWjRNmBRynoDqoDt6e+tq6tj1apVQz+YuoKISBalJpJmtifL1RmRMY3ZIiJxitpVh9rHeqgxO9eaaz8I3OTurdHr3wOLU9bPA9qA0UdjdQURERmtsYvZIiIx6u4KEu/zYXInsTaza4GHgQozO8bMrgQOADuiO84BXgLc4u5toz6guoKIiIzYmMdsEZE4jefh9szsOuBTaYvXuvttZvYa4ONmtgmoBv49loOqK4iIyIhkJWaLiMQpQ8Pt5URi7e43ADf0s24z8K+xHzTZFUQt1iIiw5KVmC0iEqcCG25v7CT72IiIiIhIQRjXjzTPipib/kVEREQkT6jFOkPUYi0iIiJSWMbzI82zQi3WIiIiIoVJLdYZohZrERERkYJiGcr/lFiLiIiIiMSgcBNrdQURERERKWzqChIzdQURERERKSwZSv8KN7FWi7WIiIiIxKhwE+sktViLiIiISAwKN7GuroYZM6CmJts1EREREZExVDShiNIZpRRVFcW63+JY95ZP3vnOUERERESkoMy7bh7zrpsX+34Lt8VaRERERCRGSqxFRERERGKgxFpEREREJAZKrEVEREREYlDwifWKFSuyXYVxSdc1M3RdM0fXNj+M589J55afdG75KVPnpsR6HH9psknXNTN0XTNH1zY/jOfPSeeWn3Ru+UmJtYiIiIhIDjMfJ4/2NrN6YPMI3joV2BNzdUTXNVN0XTMnm9f2aHevy9Kxs0Ixu086t/ykc8tPozm3fmP2uEmsRURERESySV1BRERERERioMRaRERERCQGSqxFRERERGJQsIm1mVWY2S1m9mkzu9XMJmS7TvnCzF5tZhvNbK+Z/ZeZFUfLp5rZ183ss2b2RTMrSXlPv+ukNzMrNbM1ZnZO9Lrf76q+x8NjZmeY2QfN7JLoO6lrmyfG++dhZjeYmUdlTbbrMxpmdoGZPWhm81KWjYvPr69zi5bn9ec3kr/r+aK/c4vWfT/lc7sjjuMVbGINfBW4392vBe4Avp7l+uQFM5sLXAJcCrwPuAp4b7T6/4Db3P3DwNPA9SlvHWid9PYhYF7K64G+q/oeD5GZvQV4ubvf5O4/d/c96Nrmk3H7eZhZFWGEgguicml2azRyZjYNqAZOS1uV959ff+eW75/fKP6u57yBzi1at5uez+2dsRyzEEcFMbOZwEZgkrs3Rf97OQgc5+6bslq5HGdmZwF/dfeO6PXngOOATwE/dfeZ0fJZwDpgGnBCf+vc/fDYn0XuMrMXAEuAjwNvJgSyPr+rQFt/6/Q97s3MzgY+BlzgUdAbKA6ga5tTxnvMNrNkEvM1d2/JamViYGYJoBOY7+6bxtPnl35u0bK8/vxG8nc9X/5293du7n6Rmf0/4C/Az5Pr41CoLdbnAHvcvQkguqAbgbOzWal84O73pX0BtwFbgPMI1zC53bZodtkg6yQStXpc6u63pSw+h/6/qwOtk96+CDwF3GJmd5nZmeja5pNzGKefh5kZ8FrgC8AuM3tdlqs0au7elbboHMbJ55d+buPh8xvh3/W80N+5mVkl8Argf6PXL47rmIWaWM8C9qYtawBmZqEu+e5U4GsMfE11vYfmw8CNact0XUfJzBYDpwDfcPd3An8A7gLmoGubL8bt5+HBGYSuBLcAPzCzl2a5WnHT55dfhvJ3PV+dSvhlocndFxLO8VfAb8zshDgOUKiJtQPNacvKgPYs1CVvmdlCoN7d1zDwNdX1HoSZXQiscvfdaat0XUfvOGCfuz8avf4yYOja5pNx/3m4+353vw74DD39W8cLfX55Yhh/1/NO2rkB4O473P1qQl/yWPpYF2pivQ2YlLasBtiehbrkJTMrAq4BPhIt6nVNo5/HqgjXdKB1EnwQuM3M9pjZHkJr6h2E69Tfd1Xf46EpjgoA7t4MPAMUoWubLwrp8/gSMDfblYiZPr88MMy/63mlj3NLF9vnVqiJ9Z+AmWZWARANHzMXuCertcovHwRucvfW6PXvgcUp6+cRbgBbNcg6CV4PnJxSthPuXv42/X9X9T0emkeB2uiO/qQOwjXWtc0PhfR5dAEPZbsSMdPnlx+G83c936SfW7rYPreCTKzdfSdwJ5DsrH4+8JOUjvkyADO7FngYqDCzY8zsSuAAsMPMjo82ewlwi7u3ufuq/taNeeVzlLvXu/vWZCHcdV7v7lvo57uq7/HQuPta4DdEQ2CZ2USgFPgBurZ5YTx/HmY208xeb2aJqEXwA8C/Z7teoxGdB4QuV+Pq80s/t/Hy+Q3373q26jkSfZ2bmR1vZq+M1pcCbwI+F8vxCnG4PQAzm0S4iJuA6cBH3b0xq5XKA2Z2HWEInlRr3X2pmR1NGCZuE2Gsz39PGeKm33VyJDPbBLzZ3f880HdV3+OhMbOpwH8RWlrmEm5kfFLXNn+M18/DzI4j/MfvEHA/cHP0n8G8ZGY1wBWEcauvB77s7vXj4fPr69yAo8jzz2+kf9fzQX/nRuhPfTth1JO/AZ9191i6uBRsYi0iIiIiEqeC7AoiIiIiIhI3JdYiIiIiIjFQYi0iIiIiEgMl1iIiIiIiMVBiLSIiIiISAyXWIiIiIiIxUGItIiIiIhIDJdYiklPMrM7MfmFmm83sI9muj4iIDExxu4ceECMiOcXMLgT+BEwgPCFrkbvvy26tRESkP4rbPZRYi0jOMrN7gfPy6RG6IiKFrNDjtrqCiEhOMrM5wO8KNTiLiOQbxW0l1pIBZvZvZvaVbNcjKY76mNn5ZvZXMzsnpmqNK2Z2U8z7M+By4MY49ysivSlej56ZzTOzb5jZ12PaX5GZXW1mm4ew7UVm9uBIzzXO2K24HSixzjNmttTMfmpm3zSzNWbmZnZFtuuV5tfA7dmuRIo46vMP4KQY6pKTzOxlZrbHzHaY2fPNbE70PXMz+7KZTYy2O83MHjezL5hZRbRsGrAz5iq9HviWu3dE+xfJO4rXI5KP8Xp3NC2LaX8J4AAwdwjbrgKW9rfSzI4zsx9F373rzGx6yrq4Y7fiNkqs89HPgU+6+1XAKcB3kyvM7EQzOzsblUo9trs/7u73Z6MefRlOffq7hu6+DtgTZ73i/LxGuy93/zXwVaDD3f/q7v8ArgbagAfd/UC03d8JgfxD7t4cvf0VwJ2jOoEUZvYx4BPA3Wa2Flgc175Fxpji9TDlY7x29yZge1zHcvd2YPUQt91FSML7W/8EsCJ6eZO7pybSscVuxe0eSqzziJlNBhYBhwHcvQv4OOBmVksI2paFemXt2HEawnl0jeGxsrGvHwGzzex5AO6+F/gD8PKUY00EtkffvaQl7r52lMfu5u6fcvdj3X2Zuy9x93vj2rfIWFG8zqwcjNdxjwQxnPoPtm1nNE3v9xxb7Fbc7qHEOo9EQ9esA+4yszOiZZuBvxP6Nc0HlpvZWwDM7Hoz+1j0M1B3PyozW2ZmXzSz95pZk5mtM7P3m9klZrbSzF5lZj8zs0Nmdk3K+/rcX9qxrzaz76T2NYt+Dr3FzD5jZn8ys9Oi5QMeL+X9ZmbXmNkTZnaBmT1qZrvM7M1DOMaxqfUZ5JhHXMM+TIp+2m0ws6vTruknomvz0ygQD/WaHXEsM7vUzP5uZm82s2fM7MYB9tfXZ99nfQYStWw8AVyWsrgOuMiibh/AK4FfpNSzAmhKeb3EzG40s5vM7B4LXUpeY2arzOxiM7vdzA6Y2XvM7HnROe41s5cNVj+RfKJ4XTjxui/9neNAxzGz+WZ2m4XW388P5Th9HPfXZvYDMztzkO0UuzPF3VXyqADHAZsJ/0O9DTgqZd0m4Jxo/nhgQzRfR/jf9BTC/7i3AMuidV8BHojmS4F9wJcJY1G+Dtg20P7Sj03oY/Z14DvR8nJCsjYxen0pUA/UDHS8Ps57enTMa6NtbwBagKMHOUZ6fQY8Zuo17KMOm4AvAJXR+7ZEyycA/5Oy3V8IQXFI16yfY00m/Lx3M3AacP5A+0v77PuszxC/X9cBz0bzC4DPAK3Aq1O+L5ay/SuBs6L5SmAlUBm9vje6DhOAg8D1QAXwmuj15YT/3H+A6DuoojKeCorXBRGvo/WfGOJ1HOizXgksTXmPD/F7lvw8nw/8Rx/rz4mOU5yyTLE7Q0Ut1nnGQ6viCYRA80bgCTPr6yaNZ4B/MbMS4KxoWTXhH/IcINk/9inCPxjcvQ1oAH7s7ocI/9CmD7K/9Pq10vtmiIuAZo/66AI/JfyDfPUgx0vfb3KfP4y2vYHQ//eCQY7Rqz7DOWY/fumhP91KYEbKOU6JWpTeS7imMMRr1s/57iMEsDvc/e/ufvcw9tdffYbiR8AxFrqD/DPwDeC3wOVmVg00eBR5I2cCD0TzryT88Uq2grwK+FR0nQ8Af/TQL3sVMMHd/9fDz+MPM7zPQCQvKF4XRrzuQ7/nOMBxzgdq3T1Zn1XDPOapwBvc/ZND3F6xO0OKs10BGTozKwLmuPsm4N1m9i3CHdQrgNNTt3X3VjObBXwY+O/kLtx9t5mtAl5M+B/1fEIy1f3WlPkOou5C/e2vn6qm7uNYoCSlXl1mtgmYPdDxBuPuLWa2AZgITB3GMUZ8zDQd9Pz7mQtsdPeb0zcayjUzs38hJLBJiz38ZOypdR3GZ9BvfQbj7s+Y2cOE7iC17r7BzH5EaEW6DPhVSr2T341k/72jCS1MyX3tTd11ynwnvXWhbmkyzihe9zq/QojXqfq9jgN8NicS/hMxUucBx5pZtbs39rNNV3QOit0ZpAuSX0pI6f/q7msIP7UtSd/QzM4D3u7uN3gY4SHVO4BTLQz79Cjw2cEOPMj+BrKF0AKaOgyREfoejlYp8HSGjzEUO+jdDxkzO30Y1+wXwMkppc+7y4exvz7rM/TT4X+AtwAbUupXBLyH8Gq11A4AAANNSURBVLNp0umE/qJJ24AXmllVynEH7OcnMo4pXvdWSPG633Mc4DiHoveUMjI3EVr7/6uf9Qe956Zzxe4MUmKdf95lZjNTXs8C7orm2wg3aywhDO1Ua2ZlZnZ+tP4oM5tK6Ov3LcLPOA8Q/bQYSZD2v3Qzs0H2l35sS9nHz4FG4LXRviZGx0u2fPZ3vP7MiLaZFb33riEcI7U+gx0z9TzS9dV6YYRWqCrgTjO70MJNMoN9Bt3HAma4+9qU0p5yiKKU+YH2l1rvPutjZqVmtiLq5jGQ/yX8BP1jAHdvAH4D/MV7jwbyT4RuIknJ6327hbGw3wkkxzJN/wzSP+e8HqFApB+K1xRMvB7qdezzOIRYWkr4zxeEXycws5lDiN1GaD3+F+AyM7s8eu8HzOw6Qmv9H1K2V+zOpEx03FbJTCHcEOHAfuD7wHcIPyVNitb/B7AVuJjw084zhBaCiwgtHd8hBKlfEG6kaIv210T4h/ZSoJ1wg8wc4KPR+isH2l/asd8D3Efot3ZctG4ZoaXzC4S+hsdHy/s9Xj/n79H730PomnBCyrr+jjEvtT6DHTP1GqYd+2LCTXxfjq5F8n1vidafDTwefTZfiK7zUK/ZxX2c66XR8W4HZkXLBvpMe+2rn/rMJjzI4IdD+K59Ke31a4Bz05bd1Mf7XgQ8CewFrouWvYzwXbsFmEk05BjwVsIf3m9E5/rybP8bU1GJq6B4XUjxei5wP6HV/eRBznGgz/oiwi+FfyPcSP408Pbo3PuM3YRfRdoJLdWTgT8SkvqrgGsI43l/gdANJfkexe4MFosunhQIMzsKeKu7Xx+9ThBaUa5x9+uyWrlBmJkD8z30WZQRMrMPufuIhnISkbGjeC2pFLvzg7qCFJ53EroFTE5Zdh4pN6XlopSfnvSz0wiZWZGZvQH4SbbrIiJDongtit15Rol14fkaoS/WU2b2D+D3wFZ3/2t2q9W/6GaO90Qv32xh2DcZvmrg5+6+YdAtRSQXKF4LKHbnFXUFERERERGJgVqsRURERERioMRaRERERCQGSqxFRERERGKgxFpEREREJAZKrEVEREREYqDEWkREREQkBkqsRURERERioMRaRERERCQGSqxFRERERGLw/wE5EyI1K8xSvAAAAABJRU5ErkJggg==\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/apollo-4-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/apollo-4-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/apollo-4-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 }