{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 15 - Apollo-AS-201 - 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 AS-201 suborbital flight test.\n", "https://en.wikipedia.org/wiki/AS-201" ] }, { "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-AS-201', 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,7.67,0.0,-9.03,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,60.0)\n", "# bank angle = 60 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+17YcXAAAgAElEQVR4nOzdeZwcZbX/8c/JZF9IQiYrSQgkbAnINgkoW9jxIsguKD8vciUCylVERAEREBRRFEEQBkRFQLksigEEwQso6yUhbIGwhEyAkIRJyDLZJzPn98dTnenpzNbTS1X3fN+vV72quqq66nRlcubM0089Ze6OiIiIiIjkplvcAYiIiIiIlAMV1iIiIiIieaDCWkREREQkD1RYi4iIiIjkgQprEREREZE8UGEtIiIiIpIH3eMOIF8qKyt93LhxrW6vra1l6NChxQuoBOkadYyuU/t0jdqXfo1mzpy5xN271AVrL2e3Rj9b7dM16hhdp47RddpcWzm7bArrcePGMWPGjFa3V1VVtblddI06StepfbpG7Uu/RmY2P+Zwiq69nN0a/Wy1T9eoY3SdOkbXaXNt5Wx1BRERERERyQMV1iIiIiIiedBlCutp06bFHULi6Rp1jK5T+3SN2qdr1Dm6bu3TNeoYXaeO0XXKjrl7cU9odihwBfAFd6+J1h0HXANsAdwBnOfuG6NtlcCVwHKgB3CBu9dnHreqqsrVB0hESpGZzXT3qrjjaIlytohIc23l7KK2WJvZcKA/MCVt3VjgGOAE4Fzgq8C30t52D3Cbu18AvA1cXrSARUS6MOVsEZHsFLWwdvfFwAMZq7cGTnf3me5+O3ADcCCAme0N7ODuL0T7TgfOMbN+xYpZRKQ9jY2N1Ndv1ihb8pKUs92hoSFMIiJJVfQ+1u7emPH636mvECMLgPej5YOAeWn7LogWE/mVqYh0PX/5y1+oqKjgpJNOijuUgkhKzn7mGejeHfbfP9cjiYgUThJvXpwM3BQtbwUszdheB4zKfFNtbS1VVVWbpurq6gKHKSICPXv2BGD9+vVZva+6unpTvgIq8x9Z0RQlZ5uFeZFvCxIRATqesxP1gBgzmwDUuvsr0SoH1mbs1gvY7DvXoUOHagBzESm6VGG9YcOGrN43bdq0TXfbm9mSvAdWBMXM2d2iZiAV1iISh47m7MS0WJtZBfA14HtpqxcAg9P2MaAf8FFxoxMRaVmvXr2A7FusS12xc3aqxbqxse39RETilJjCGvgO8HN3T//t9BiwQ9rrccAGQE3TIpIIqcI62xbrMlDUnK2uICJSCopeWEctGACWtu5CYBbQx8y2NbPTzWyCu88AFprZztGuRwDXu3uX+w0mIsnU2T7WpSIpOVtdQUSkFBS1j7WZDQBOjV6eZma/JnyV+KOMXee4+23R8heAS8yshjCe6kXFiFVEpCPKuStIknK2uoKISCkoamHt7nXAb6Ip5Ypoau0984H/KnBoIiKd0tmbF0tBknK2uoKISClIUh9rEZGSU84t1kmiriAiUgpUWIuI5KCcW6yTRF1BRKQUqLAWEcmBWqyLQy3WIlIKVFiLiORALdbFoRZrESkFKqxFRHKQ3mLtak4tGN28KCKlQIW1iEgOKioqqKiowN1paGiIO5yypa4gIlIKVFiLiOSo3B8SkwTqCiIipUCFtYhIjnQDY+GpK4iIlAIV1iIiOdINjIWnriAiUgpUWIuI5Egt1oWnriAiUgpUWIuI5Egt1oWnFmsRKQUqrEVEcqQW68JTi7WIlAIV1iIiOVKLdeHp5kURKQUqrEVEcqQW68JTVxARKQUqrEVEcpQqrNViXTjqCiIipUCFtYhIjvSAmMJTVxARKQUqrEVEcqQW68JTVxARKQUqrEVEcqQW68JTVxARKQUqrEVEcqSbFwtPXUFEpBSosBYRyZGG2yu8VFcQtViLSJKpsBYRyZFarAtPLdYiUgq6F/uEZnYocAXwBXevidb1Aa4G6oDBwAXuvrK9bSIiSVDOLdZJydm6eVFESkFRW6zNbDjQH5iSselG4Gl3vxB4ALi5g9tERGJXri3WScrZunlRREpBUQtrd19MSLSbmNko4IvA9GjV48DRZjaurW1FCVhEpAPKtcU6STlbfaxFpBQUvSuIuzdaqukhmAoscfc10faNZjYPOACob2NbTTHjFhFpTbm2WENycnaPHmFeX5/LUURECqvohXULtgKWZqyrA0YBG9vYJiKSCF3sATGx5OzoSwG6xiUWkVKVhMLagbUZ63oRWj5oZ9smtbW1VFVVbXo9bdo0pk2blscwRURa1pkHxFRXV1NdXZ16WZn/qAomlpydarHeuDHcwNi8EV1EpLA6mrOTUFgvINw5nm4A8BEhgbe2rZmhQ4cyY8aMggQoItKWzrRYpxeSZrakIIEVRiw52wy6dw+FdX19Uwu2iEgxdDRnJ2Ec6yeAUdEQTZhZD2As8FQ720REEqGLPdI8tpytftYiknRFL6yt6S4YA3D3RcCDwMHR+kOA+9x9QVvbihiyiEibyvnmxSTlbPWzFpGkK2pXEDMbAJwavTzNzH7t7rXAWcDVZrYrMAJI72jX1jYRkdiV63B7ScvZarEWkaQramHt7nXAb6Ipff0y4IxW3tPqNhGRJCjXFuuk5Wy1WItI0iWhj7WISEkr1xbrpFGLtYgknQprEZEclWuLddKoxVpEkk6FtYhIjrrYA2JioxZrEUk6FdYiIjnqYsPtxUYt1iKSdCqsRURypBbr4lBhLSJJp8JaRCRHarEujn79wnz16njjEBFpjQprEZEc6ebF4ujfP8xXrYo3DhGR1qiwFhHJkYbbKw4V1iKSdCqsRURypBbr4lBhLSJJp8JaRCRHarEuDhXWIpJ0KqxFRHKkFuviUGEtIkmnwlpEJEdqsS6OAQPCXIW1iCSVCmsRkRx1796dbt260djYyMaNG+MOp2ylWqzr6uKNQ0SkNSqsRUTyQA+JKbzBg8N86dJ44xARaY0KaxGRPNBDYgpv2LAwr62NNw4RkdaosBYRyQO1WBfe0KFh/vHH8cYhItKavBTWZvYf+TiOiEipKqUW61LN2anCWi3WIpJU3bPZ2cyOAC4GRtBUlHcDhgN98huaiEjpSOKQe+WWs4cMCfOlS6GhASoq4o1HRCRTVoU18EfgJ8AsoDFaZ8AX8xmUiEipqYiqvIaGhpgjaaascnaPHrDllvDJJ7BkCQwfHndEIiLNZVtYv+fuv8hcaWZz8hSPiEhJ6tYtNAi7e8yRNFN2OXvMmFBYz5+vwlpEkifbPtZnmdkJLaw/Oh/BiIiUKjMDoLGxsZ09i6rscva224b5vHnxxiEi0pJsC+t9gFvNrCFtagR+U4DYRERKRkJbrMsuZ2+zTZirsBaRJMq2K8jlwDRgBpDqSGjA/8s1EDPbCfgG8C6wHVDt7i+bWR/gaqAOGAxc4O4rcz2fiEg+JbTFuuxytgprEUmybAvrvwMPuHuz297N7JY8xPJH4PPuvsDMxgKPAjsBNwKPuPvd0R3uNwOn5OF8IiJ5k9AW67LL2RMmhPmcku0lLiLlLNvC+s/A983sf9PWGXAc8M0cY9kJGBAtrwUGmtkowt3rX4/WPw7cZ2bj3L0mx/OJiORNQlusyy5n77prmL/yCrhDdNlFRBIh28L6B4Sv/E4DUs0yRhgjNdckfS+hL+BngVOBc4CpwBJ3XwPg7hvNbB5wAFCT4/lERPImoS3WZZezR4wIjzb/+OMwMsi4cfk4qohIfmR78+LPgDHuPs7dt4mmcYTWj1ydDawDXgRWuft9wFbA0oz96oBRmW+ura2lqqpq01RdXZ2HkEREOibbFuvq6upN+QqoLFBYZZezzWC33cLyrFmdD15EJBsdzdnZtlh/4O51LawfkeVxWtIH+BMhMV8btXI44SvGdL2A+sw3Dx06lBkzZuQhDBGR7GXbYj1t2jSmTZsGgJktKVBYZZmzp0yBf/wDnn4ajj22U4cQEclKR3N2ti3WZ2euMLMTgWuzPE5L7gDuc/fLgasICXsh4a7ydAOAj/JwPhGRvEloH+uyzNkHHhjmTzyRz6OKiOQu28J6SzP7AYCZ9TWzW4FbCV8FdpqZVQK7uvvyaNUVwEpgPjAqGr4JM+sBjAWeyuV8IiJdRFnm7E9/Gnr2hJdfhqWZHU9ERGKUVWHt7kcCj5nZ9cDLwPbArkBLT/bKxifAOjPbKm3dUuBV4EHg4GjdIYQWkgU5nk9EJK82bNgAQK9evWKOpEm55uw+fWD//cOoINOn5/PIIiK5abOPtZnt38qm+cA7wM+BrYFjgHM7G4S7N5rZMcAlZjYTGA6c7+4rzews4Goz25XQL3BaZ88jIlIoqcK6Z8+escXQlXL2CSfA44/DPffAaacV4gwiItmztm60MbPnCWOVLgda6zhowEh3j7WZpqqqynXzoojEZcKECcydO5d33nmHCamnmHSQmc1096pcY+hKOfvjj2HkSKiogIULYciQPAYnItKGtnJ2e6OCXA7McPeP2zlBrl8rioiUtCS0WNOFcvawYXDYYfDII/C738F3vhN3RCIi7fSxdveH20vQ0X735i8kEZHSU18fRpSLs7Duajn7G98I8xtvhIaGeGMREYHsRwUREZEWJKTFuks54gjYdluYNw/uvz/uaEREVFiLiOTFunXrgGSNClLuKiqauoD88IdqtRaR+KmwFhHJUX19PWvWrKGiooK+ffvGHU6X8l//BePGwZtvwh13xB2NiHR1ORfWZra7md1hZueZ2bjcQxIRKS0rVqwAYNCgQZuewJhU5Zaze/aEyy4Ly9/9LnzySbzxiEjXlo8W6zOBQ9z9GuAzqSduiYh0FcuXhwcQDhw4MOZIOqTscvapp8J++4Uh+M4/P+5oRKQry0dh/Q1gz2j5z0CPPBxTRKRkpArrQYMGxRxJh5Rdzu7WDW6+ObRe33Yb/M//xB2RiHRVWRfWZvZZM/tStDwF2C71uFp3b3T3lXmOUUQk0ZLcYt1VcvZOO8E114Tlr34V3n033nhEpGvKqrA2syuBe4BjAdz9/4Avm9kXCxCbiEhJWLhwIQAjRoyIOZLmulrO/vrX4fjjoa4Ojj4ali2LOyIR6WqybbH+DDAa+Hfaur8DV+UtIhGREvPhhx8CMHr06Jgj2UyXytlm8NvfwqRJYZSQ44+HaHhxEZGiyLawfszdlwOetq7kH40rIpKLBBfWXS5nDxwIDz0EI0bAE0+EGxujh2KKiBRctoX1e2Z2IbCDmZ1sZvcBXwd+mv/QRERKw4IFC4BEFtZdMmdvvTU8+CBssQXccw988YsqrkWkOLIqrN39z8D/AhXAl4HlwGHufkMBYhMRKQnz5s0DYMyYMTFH0lxXztl77gmPPRZasO+9F447DlatijsqESl33bN9g7s/Dzyfem1mFWZ2sLv/M6+RiYiUgIaGBt566y0Adtxxx5ij2VxXztlTpoTi+ogjQgv2AQfA9OkwalTckYlIuWqzsDazvwH92zlGJbAEKPskLSKS6b333mP9+vWMGTOGAQMGxBqLcvbmJk+G556DI4+El16Cqiq46y6YOjXuyESkHLXXYv1qtM9bQCNwBjALmJG2zzbAlgWJTkQk4WbPng3ApEmTYo4EUM5u0fbbh+L6hBPgqafg4IPDY9C//32oqIg7OhEpJ+0V1j8DLLqrHDMb7+6XZO5kZg8UIjgRkaSbNWsWkJjCWjm7FZWV8PjjcOml8OMfww9+AA8/HIbn22mnuKMTkXLR5s2L7r4ilaAjw82s2d/3ZjYVmFKA2EREEu+ZZ54B4DOf+UzMkShnt6d7d7jiCnjkkdDP+rnnYLfd4MorNd61iORHtsPt3Qu8YmY3mdnPzOxB4DHguvyHJiKSbBs3buT558N9gfvss0/M0bRIObsFhx0Gs2eHR59v2AAXXxweKvPXv4J7++8XEWlNtsPtPQYcRrjxZQfC0E0nuPtPChCbiEiizZo1i9WrVzN+/HiGDx8edzibUc5u3aBBcMstoXvIjjvCu+/CsceG/tczZrT/fhGRlmTbYo27f+TuF7v70e5+KjDdzPYsQGwiIok2ffp0AA499NCYI2mdcnbbDj4YXn0Vrr8ettwyPK1x8mQ46igV2CKSvazGsTaz21pYPRRoAI7JR0BmtjewL/Au8DSwGrgaqAMGAxe4+8p8nEtEJBd/+9vfAPj85z8fcyQtU87umB494BvfCE9ovOoquOGGMO71gw+GYfq++13Ybz8wiztSEUm6bFustwPeB+anTb2AmfkIxsy+AnzO3X/u7n919yXAjcDT7n4h8ABwcz7OJSKSi7lz5/LKK6/Qv39/DjzwwLjDaY1ydha23BKuvhrmzYPzz4e+feGhh8KDZaqq4I47dJOjiLQt2ycvftXd30pfYWY7AqflGoiZHQB8CTg0bd0o4IvA16NVjwP3mdk4d6/J9ZwiIp31u9/9DoBjjz2WXr16xRxNq5SzO2HYsFBgf+c78Otfw29+Ex4u8//+H1xwAZxxBpx+OowdG3ekIpI02d68+FYLq9cDX81DLL8A3gSuN7NHzGwfYCqwxN3XROffCMwDDsjD+UREOmXjxo2bCuszzjgj5mhap5ydm2HD4PLL4f33w42OEyfCRx+Fh8uMGwef/Szcd59asUWkSVaFtZnNM7P30qb5hMT6l1yCMLMdgN2BW9z9G4RH7T4CjAGWZuxeB4zKPEZtbS1VVVWbpurq6lxCEhFp1f33389HH33E9ttvz7777tupY1RXV2/KV4THjOedcnZ+9OkThuZ7/fUwisjJJ4d+2Y88Ep7mOHo0nHMOPPushusTKVcdzdnmWWQBM6sG7gJSb2oEFrv72znEipkdB1S7e2X0ug9QC1wOHO/ue6Xt+xJwl7v/PP0YVVVVPkO3cItIgTU2NrLbbrvx2muvcdNNN/G1r30t52Oa2Ux3r8pDeJnHVc4ukCVLQp/rW26BN95oWr/11qHwPuUU+NSndMOjSDlqK2dne/PiHe7+pLs/FU3/Dse3CTnG2J20/t7uvhZ4B6gg3FWebgDwUY7nExHplAceeIDXXnuNUaNGcdppp8UdTnuUswukshK+9a3Qij1jBpx3Hmy1FcyfDz/9aXii48SJ8P3vw/PPQ2Nj3BGLSDFkW1gf0cK6xcA9OcbxKjDQzNKfsLCRkIxHRa0hmFkPYCzwVI7nExHJ2rp16zj//PMB+N73vpfkmxZTlLMLzAz23BN+/vPQF/upp+Css0LhPWdOGL7v058ORfe0aWGUkXXr4o5aRAqlQ4W1mX3DzJYDF5hZQ/pE6E+3OJcg3H0O8HfghOh8g4CewJ3Ag8DB0a6HAPe5+4Jczici0hm/+MUvmDt3LhMnTuTMM8+MO5xWKWfHo1s32H9/uPHGcJPj//4v/Pd/h+4hixaFbiOf+1wouo8+Oow48vbb6pctUk463MfazD4F/CdwXcamDe6+MOdAzCqBXwEzCC0ct7j7G2Y2mPCwgRpgBPB9d1+V+f5S768nIsk2e/ZsqqqqWLduHY8//jgHH3xw+2/qoEL0sVbOTg738HTHBx4I00svNd++9dZw2GFhOuigMJ62iCRXWzk725sXe7h7fQvrx7r7+znEmLOulKRFpLjWr1/P3nvvzcsvv8xXvvIVbrutpQcadl4Bb15Uzk6gBQvgscfgH/8I8yVLmrZ16wZ77BFavvffH/bdF4YMiS9WEdlcpwtrM/sO8LC7vxG9vohwc0qz3YAD3X1qfsLtnK6cpEWkcNydM888k+rqasaPH8+sWbMYMGBAXs+Rr8JaObv0NDbCyy/Do4+GQvuZZ6A+40+hnXduKrT32w9GbTZ4oYgUU1s5u70nL24DDEp7vQfQD1hE0/BNRvi6T0Sk7Fx//fVUV1fTu3dv/vSnP+W9qM4z5ewSk2qh3mOPMILI6tXw3HPwr3+F6fnnw8gjr78e+m4DjBkDe+3VNO25Z3j8uojEr83C2t2/nrHqSuB1d2/2nCkzm5jvwERE4nbPPfdw7rnnAnDbbbcxefLkmCNqm3J26evXDw45JEwA69fDiy+G0Ub+9a9QdH/wQZjuvTfsU1EBu+zSvNDeaafwEBsRKa72Wqwz7eLuzW67iJ7ANR54o+W3iIiUnr/85S+ccsopNDY2ctlll3HKKafEHVJnKGeXuF69Qj/rffeFiy6ChoYwjN8LLzRNr70WupO8/DLcfHPT+3beGXbfPbSG7757eGCNWrZFCqtDhbWZ7Uf4+vBQM5uXsXlLwl3n0/Mcm4hILO68806+8pWv0NDQwPe+9z1+8IMfxB1SVpSzy1dFBUyaFKbTTw/rVq+GmTNDkf1//wezZsHcuWHdzJlN7+3WDXbYoanQ3n330NI9dGg8n0WkHHW0xXo5cDuh/94+Gds2ALfmMygRkTi4O1dddRUXXnghAN/97nf58Y9/jJXec6mVs7uQfv2abm5MWbEitGDPmtU0vfEGvPlmmO68s2nfYcNCgb3LLqGVe5ddwlMj+/cv/mcRKXUdKqzd/TUzmwqc6O5KyCJSdlatWsW0adP405/+hJnxy1/+km9+85txh9UpytkycCAccECYUtauDTdBpgrtl18Orz/+GP75zzCl23bb5sX2LrvAdtup77ZIWzrcx9rdV9BKK4eZ7evuT+ctKhGRInr99dc58cQTmTNnDv369eP222/nuOOOizusnChnS6Y+fWDy5DClNDaGR7G/9lqYXn89zOfMgffeC9MDDzTt3707jB8PO+4YupWk5jvsoPG2RaCdwtrM/ga092WQAZOAYfkKSkSkGDZu3MjPfvYzLr30UjZs2MCkSZO499572XHHHeMOrVOUsyVb3brBuHFhOuqopvUbNoTHracX26+9BjU18NZbYcpUWdm82N5++9Dqve22obuKSFfQXov1W4S+egtoGgM1kwGfy2dQIiKFNmvWLM444wxmRnd3nXHGGfzyl7+kX2lXAMrZkhc9e4YuIDvv3Hz96tXwzjtNxfVbb4XW7bfeCk+QXLIkPOQm04gRoaV7223DPH152DAovdsYRFrWXmH9U6DO3de3tZOZ6e5yESkJixcv5qKLLuK2227D3Rk7diy33norhx56aNyh5YNythRUv36w225hSucOH33UVGS/9VYowN97D+bNg0WLwtRS0d2vX1ORve22ofV8662b5gMHFuOTieRHew+IWdLB4xyL7jIXkQRbuXIl1113HVdffTV1dXV0796dc845h8suuyzpT1PsMOVsiYsZbLVVmA4+uPm2hgZYsCAU2XPnNk2p18uWNXU1acmgQaHATi+2U/Ottw59u9XiLUmR7QNimjGzscDXgW+gJC0iCbRy5Uquv/56rrnmGpYtWwbA5z73Oa655hq23377mKMrLuVsiUNFBYwdG6apUzffvmxZU6H93nswf37oy52aL18epldeafn4/fptXmynlseOheHDQwwixZB1YW1mFcDRwBnAocBiQp8+EZHEmD9/PjfccAO33HILy5eHFLXffvtx6aWXctBBB8UcXfEoZ0vSDR4MVVVhyuQe+m2nCu3MorumBurqwhjdb7TyLNHu3UNL+ujRMGZMmNKXx4wJD8np1q1wn1G6jg4X1ma2DfBV4CvAEGA+cBTwKPDpgkQnIpIFd+fpp5/muuuu4/7776exsRGA/fffn0svvZSpU6eW4sNeOkU5W8qBWSh6hw5tPkxgintozW6p6J4/Hz74AGprm7a3pmfPUHy3VHinlisr1eVE2tduYW1mJxJaOg4C3gV+AfwBuNDdH4l2e7ZgEYqItOOjjz7i9ttv57bbbuOdd94BoHv37px88sl885vfZMqUKTFHWDzK2dKVmIUW78GDN7+hMmXdOvjww1Bkp+aZy598Em6ynDev9XP17t3Uj3zUqM2nrbaCkSM1tGBX15EW67WER+A+C3zV3d8GMLPWhnISESm4FStW8Le//Y0///nPPPLII5tap0eOHMnpp5/O2WefzahRo2KOMhbK2SJpeveGCRPC1Jo1a1ovulPLy5c33XjZloEDWy+8U8sjRkCvXvn9nJIM7RbW7v4g8GB008sZZjYOeADY9FBTMzN3V9IWkYJKFdP33HMPjz76KBs2bACgR48eHHvssZx++ukcdthhdO+e033ZJU05WyR7ffuGB9q0dT/zqlWhyP7oo7anFSvC9OabbZ+zsnLz4nvEiKZp5Mgw79s3v59VCiubR5q/D/zAzLoDxwCjzexHwG+Bw4DqwoQoIl3Z22+/zcMPP8zDDz/MU089tamYNjOmTp3KiSeeyEknnURlZWXMkSaLcrZIfvXvH54q2daDWd1Dt5LMYnvBguavFy1qeqDOq6+2fd4BA5oX3JlTqgAfOjTcqCnxyvqfwN03AvcC95rZBOBcQn8+JWkRydmqVav497//zSOPPMLDDz/Mu+++u2lbejF93HHHMWLEiBgjLQ3K2SLFYxbG1R4yBHbZpfX9GhrCTZWZRffixU0P01m0CBYuDKOe1NWFB+60d+6hQ9svwEeMCN1VdCNmYeT0t427vwt808xm5ikeEeli1qxZwzPPPMMTTzzBk08+yYsvvsjGjRs3bR88eDBHHHEE//Ef/8ERRxyhlukcKGeLJENFRVORu8cere+XGvUkvdhubaqthY8/DlN7reC9erXfCp4qxHv2zO9nL3d5+dLA3W/Px3HMrCfwIvBNd3/SzPoAVwN1wGDgAndfmY9ziUg8Fi1axHPPPcfzzz/Ps88+ywsvvEB9ff2m7d26dWPy5MkccsghHHnkkey1115dus90IeQrZ4PytkghpY96stNObe9bXx+K644U4XV17Q9BmFJZ2VRojxwZ+oK3tNy7d34+c6lL2m+r84Fxaa9vBB5x97vN7AjgZuCUOAITkeytX7+eV199leeee25TMV1TU9NsHzNjjz324MADD+TAAw9k3333ZeDAgfEELJ2hvC2SAD16NN0E2Z5VqzbvdpLe/WTRoqauKam+4K09cj5l0KDmhXZm4Z2a+vfPz+dNqsQU1mb2GWAhsCx6PQr4IuHxuwCPA/eZ2Th3r4klSBFpVV1dHS+//DKzZs3ipZdeYtasWbzxxhvNunUA9O/fnylTpvDpT3+avffem3322YfBgwfHFJrTHpwAACAASURBVLXkQnlbpDT17x+m8ePb3i/VF3zhwlBoL1zYNKW/XrSo6dHzrT0BM2XAgKaH7owe3Xw5NS/ltpVEFNZm1g84wd2/bWaXRKunAkvcfQ2EG3DMbB5wAFATS6AiQkNDA3PnzuWNN95g9uzZvPrqq8yaNWvTg1nSmRk77bQTe+2116ZCetKkSVRUVMQQueST8rZI+UvvC7777q3v19jYNBpKa8V36nVdXRiKsK3hCFPFd0tF95gxMG5cch/Ek4jCGrgAuCpj3VbA0ox1dUCXfOKDSLE1NDTw3nvvbSqgU9OcOXNYv379Zvv37NmTnXfemd1335099tiD3XffnU996lP0S2r2k1wpb4sIAN26hb7YlZXwqU+1vp87LFsWRkJJf/hO5rwjxffw4bDttrDNNmGePo0aFf4oiEPshbWZHQ7McPePMzY54Qli6XoB9bSgtraWqqqqTa+nTZvGtGnT8hmqSNlpbGxkwYIFvPPOO5tNc+fO3TRmdKYxY8YwceJEJk2atKmYnjhxIj11+3iHVVdXU129acS7khrqJB95WzlbpOsxgy23DFNrwxGmRkJpqej+8EN4/32oqQn9vxcvhuee2/wYPXvCdtvBpEmw885N8/HjO19wdzRnW9wP3zKzx4D0LxgGA6uAXwBfcvft0/Z9B/ihu9+VeZyqqiqfMWNGocMVKTlr165l/vz51NTUUFNTw7x585oVz2vXZtZBTbbaaismTZrUbJo4cSJbbLFFET9B+TOzme5e1f6eyZCPvK2cLSKd1dgYupW8917zad68MF+0qOX3DRwI++4L++8PRx/d9sN+2tJWzk5CYT2U0KKR8hzwbeAF4A1gqLuvNbMehMS9rbsvyDyOkrR0VWvWrGlWONfU1DR7vXjx4jbfP2zYMLbbbrvNpgkTJtC/3G/fTogSLKxzztvK2SJSKKtXw5w58PrrYZo9O4ztvSCjejzqKLjrruxHKmkrZ8feFcTda9Nfm1kDUOvu75vZg8DBwIPAIcB9LRXVIuXI3Vm2bBkffvghCxYsYMGCBZuW0+fLli1r8zg9evRg7NixjBs3btM0YcKETcWzhraTbClvi0iS9esHe+4ZpnQ1NVBdDT/5SXg9fTrcey+cdlr+zh17Yd2Os4CrzWxXYASgDnhS8hoaGliyZAmLFy9m0aJFzeYLFy5sVkSvW7eu3eO1VDinTyNHjtQoHFJMytsiEov6+vDkydQQgAsWhK4hc+eG+VtvwZo1TftXVsI+++Q3hsQV1u4+Lm15GXBGfNGIdMyaNWtYunQpS5YsYcmSJXz88cctFs6LFi2itraWxsbGDh13wIABjB49mq222mrTPH159OjRVFZW0q1btwJ/QpHWKW+LSD5t2BBGBqmrgxUrwlB+qWnZsuavP/kkPMBm0aIwb6+H89Zbh5bsww+HE08MT7XMp8QV1iJxcnfWrFnDJ5980qxQXrJkSZuv27oBsCVDhgxh+PDhjBgxotl8+PDhzQpo3SQoIiJJtXFjaAFOn1av3nxdaqqrg5Urm4rm9OX01y2M6Noh3bqFYfhSY2+PHNl8GL7ttoMhQ/J7DTKpsJayU19fz/Lly1m+fDnLli1rNm9v3fLly6mvb3FExzb17NmToUOHMmTIECorKxk6dOhmRXNqPmzYMHr06FGATy4iIl1NQ0MoRNetC1P6cubrzmxbu7b1wrkTvy47pKIiPCRmiy3CNGRIaFlODdXX0jRyJAwdGt/41SkqrCV2jY2NrF27ltWrV1NXV0ddXR0rV67ctJzNupUrV3aoX3JbevfuzeDBg6msrGw2pYrmltb169cPM8vTFRERkTg0NoZW2Pr6pnn6cvq6DRtan9rbnq991q0LMcWlWzfo27f51K/f5utSU//+oVBOFc0DBrS83Lt3GPO6FKmwljY1Njaybt061qxZw+rVq1mzZs2mKf11LttyLYQzVVRUMGjQoE3T4MGDN1tubd3AgQPp3bt3XuMRESkX7qGFNL3AbK0IzVzuSLGauS6f+3fkWB28/SVRzEIhmpp69crfcq9ebRfMPXuWbgFcKCqsy0RDQwPLly/nk08+4ZNPPmH58uXtFrUdKYCz7TvcWX369KFv374MGDCALbbYggEDBjSbMte19bpv375qPRaRkpLqq7p27eZTS+vXrQstlvkscDu6X7nr0QO6d295nlru2bPlqUeP1rfl6z3p+/fuHeLRr7zkUGGdcPX19dTU1PDBBx/wwQcf8OGHH26a19bWbrrJbvny5RTqYT+pordfv3707du32ZS+LnN7e/unXvfu3VujWohIWXCHpUvDU+FS06JFYeSCFSvCo5ozpxUrSqtg7dat5WKzpeWO7tdaAdvRbdnu39q2igoVqZIbFdYJ4e7Mnz+f559/nldffZU333yTN998k7lz57Kxgx2oBg0axJZbbsmQIUMYNGjQpiK2vYK3reK4T58+KnpFRDJs3AivvQYvvRTGxp0zJ0w1NZ0rkrt1gz59mk99+7a+rnfvptbLfBa4HdlPvxJEWqfCOibuzltvvcXf//53nnrqKZ5//vkWHz1tZmy99dZsvfXWjB49mjFjxjB69GhGjx7N8OHDmxXS3bvrn1NEpFDefRfuvx8efRReeCGMjtCSwYPDCAWjRoVp5MgwqsHAgTBoUNOUej1wYOjLqpZSkdKnSqyI3J0XXniBu+66i+nTp1NTU9Ns+5AhQ9hrr73YY489mDhxIjvttBM77LADffr0iSdgEZEurr4e7rkHrr0WXnyx+bZtt4UpU2DiRNhxxzCNHx9alUWka1JhXQSLFy/m5ptv5vbbb2fu3Lmb1g8ZMoTDDz+cww47jH322Yfx48frpjsRkYR47DH4+tfhnXfC6y22gKOOgs9/HvbfPzyIQkQknQrrApo9ezbXXHMNd955Jxs2bABg5MiRnHLKKZx44olMnjyZirhHMhcRkWYaG+Gii+Cqq8Lr7baD88+HU08N/ZtFRFqjwroA3nvvPX74wx9y55134u6YGccccwxnn302Bx10kIppEZGEcoczz4RbbgkjRFx2WSiqe/aMOzIRKQUqrPNo9erVXHrppfzqV7+ivr6eHj16cMYZZ3DuuecyYcKEuMMTEZF2/OpXoaju3Rv+9jc49NC4IxKRUqLCOk8efvhhzj77bObPn4+Z8eUvf5nLLruMcePGxR2aiIh0wNy58P3vh+U77lBRLSLZU2GdozVr1nDuuedSXV0NwG677UZ1dTWTJ0+OOTIREcnG5ZeHJxp+6Utw/PFxRyMipUiFdQ5mz57NF77wBWbPnk2vXr244oor+Na3vqXxpEVESsz8+XDnnaFf9Y9+FHc0IlKqVAF20kMPPcTJJ5/MqlWr2GGHHbj77rvZdddd4w5LREQ64Q9/gIYGOPlk2GabuKMRkVKlB5Nmyd259tprOfroo1m1ahWnnHIKM2fOVFEtIlKi3OH3vw/Lp58eaygiUuJUWGfB3bn44os599xzaWxs5LLLLuPOO++kX79+cYcmIiKd9MYbMG8ejBgBBx0UdzQiUsrUFaSD3J3zzjuPX/7yl1RUVPCHP/yBL33pS3GHJSIiOfrnP8P84INDH2sRkc5SYd0B7s63v/1trr32Wnr06MHdd9/NscceG3dYIiKSB08/HeZTp8YahoiUARXWHfCTn/yEa6+9lp49e3L//fdz5JFHxh2SiIjkyfPPh/mnPx1vHCJS+hLTx9rMjjOzeWa21Mx+ZWbdo/WVZnazmf3UzH5hZj2KGdett97KRRddhJlx5513qqgWESG5OTtbCxbABx/AFlvATjvFHY2IlLpEFNZmNhY4BjgBOBf4KvCtaPM9wG3ufgHwNnB5seJ64oknOPPMMwG44YYbOOGEE4p1ahGRxEpqzu6MF14I8ylToFsifiOKSClLSleQrYHT3X0jMNPMdgYONLOngR3cPUp9TAfeMrMr3H11IQOaP38+J510Eg0NDVxwwQWcddZZhTydiEgpSVzO7ix1AxGRfErE3+fu/u8oQacsAN4HDgLmpe23IFqsKmQ869at47jjjmPJkiUcfvjhXHnllYU8nYhISUlazs5FqrDee+944xCR8pCIwroFk4GbgK2ApRnb6oBRmW+ora2lqqpq01RdXd3pk1944YW89NJLbLPNNtx1111UaPwlEcmz6urqTfkKqIw7nhzFmrM7q74eZswIy3vtVfTTi0gJ6WjOTkpXkE3MbAJQ6+6vmJkDazN26QXUZ75v6NChzEhlyBz885//3DRW9Z///Ge23HLLnI8pIpJp2rRpTJs2DQAzWxJzOJ0Wd87OxWuvwdq1sN12MGRIrKGISMJ1NGcnqsXazCqArwHfi1YtAAanbTegH/BRIc6/YsUKTjvtNAAuueQSpkyZUojTiIiUhbhzdq7UDURE8i1RhTXwHeDn7r4+ev0YsEPa9nHABqAgzRyXXXYZH374IVOmTOHCCy8sxClERMpJrDk7V3/7W5jvs0+8cYhI+UhMYW1mFwKzgD5mtq2ZnQ4sBxZGd5wDHAFc7+4b8n3+2bNnc91112Fm3HTTTXTvnrheMiIiiRF3zs7V++/DP/4BPXuCRlIVkXxJRPVoZhcDP8pYPcfdbzOzLwCXmFkN0B+4qBAxfPe736WhoYGzzjqL3XffvRCnEBEpC0nI2bn6/e/BHY49Vv2rRSR/zN3jjiEvqqqqvLM3wrz00kvsueee9O3bl5qaGoYOHZrn6EREWmdmM909sUPSFUIuOTtXjY0wfjzU1MCjj8Jhh8UShoiUqLZydmK6gsTpxz/+MQBnnXWWimoRkTL35JOhqB47Fg45JO5oRKScdPnCeuHChfz1r3+le/funHfeeXGHIyIiBfbb34b5V76ix5iLSH51+ZTyhz/8gYaGBo466ihGjhwZdzgiIlJAGzfC9Olh+T//M95YRKT8dPnC+o9//CMAp59+esyRiIhIoc2YAXV1sP32sM02cUcjIuWmSxfWc+fO5Y033mDgwIEcfvjhcYcjIiIF9tRTYX7ggfHGISLlqUsX1tOj7wOPOOIIevToEXM0IiJSaC+9FOZ6sK6IFEKXLqwfeughAI466qiYIxERkWKYNSvM99gj3jhEpDx12cK6vr6eZ599FoBDNN6SiEjZW7kS3nknPG1x4sS4oxGRctRlC+tZs2axZs0att9+e4YPHx53OCIiUmCvvBLmO+8cimsRkXzrsoX1v/71LwD222+/mCMREZFieO21MN9113jjEJHy1WUL62222YYjjjiCw/QsWxGRLmH0aDjqKFB7iogUSve4A4jL8ccfz/HHHx93GCIiUiRHHx0mEZFC6bIt1iIiIiIi+aTCWkREREQkD1RYi4iIiIjkgQprEREREZE86DKFdXV1ddwhJJ6uUcfoOrVP16h9ukado+vWPl2jjtF16hhdp+yosJZNdI06RtepfbpG7dM16hxdt/bpGnWMrlPH6Dplp8sU1iIiIiIihWTuHncMeWFmtcD8NnapBJYUKZxSpWvUMbpO7dM1al/6Ndra3YfGGUyxdSBnt0Y/W+3TNeoYXaeO0XXaXKs5u2wKaxERERGROKkriIiIiIhIHqiwFhERERHJg7IvrM2sj5ldb2Y/NrPfmNkWcceUNGZ2nJnNM7OlZvYrM+sed0xJZWY9zewVM5sadyxJZmZ7m9l3zOwYM6uMO54kMbOdzOwGMzvXzG40s93ijqkUKJd3jJldYWYeTa/EHU+SmNmhZvaCmY1LW6efqzQtXaNovX6uOqgrFFA3Ao+4+91mdgRwM3BKzDElhpmNBY4BTgAmAb8BPgB+HmdcCXY+MC7uIJLMzL4CjHf3i+OOJaH+CHze3RdE//8eBXaKOaZSoFzeDjPrR7jR7NBoVWduDi1LZjYc6A9Mydikn6tIa9dIP1fZKeubF81sFDAPGOzua6KW2BXAJHeviTW4hDCz/YDn3H1j9PpqwvU5Mt7IksfMPgPsCFwCnObuT8YbUfKY2QHAD4BDvZyTSw7MbDWwp7vPMbOhwCvuPiruuJJMubxjzOxb0eJN7r4u1mASyMy6AQ3ANu5eo5+rzWVeo2idfq6yUO5dQaYCS9x9DUBUPM4DDogzqCRx93+niurIAuD9uOJJqugv9hPc/ba4Y0m4XwBvAteb2SNmtk/cASXQvcCtZjYAOBU4J+Z4SsFUlMvbZGYGnAxcAyw2sy7Z6toWd2/MWDUV/Vw1k3mN9HOVvXIvrLcClmasqwPUOtS6ycBNcQeRQBcAV8UdRJKZ2Q7A7sAt7v4N4J/AI9HXi9LkbGAd8CKwyt3vizmeUqBc3g4P9iZ8ZX89cKeZfTbmsJJOP1ft0M9V9sq9sHZgbca6XkB9DLEknplNAGrdXTcmpDGzw4EZ7v5x3LEk3CTgE3d/NXr9a8CAz8cXUiL1Af4E3AVca2aHxBxPKVAu7yB3Xxbd3/AT4Fvt7d/F6eeqg/Rz1XHlfvPiAmBwxroBwEcxxJJoZlYBfA34XtyxJNB3gN3DN2JA+Jl6wMyudPer4wsrcbqTllPcfa2ZvQMMiS+kRLoDONndl5uZA38ys3HuvjruwBJMuTx71wFPxh1EwunnKnv6uWpHubdYPwGMMrM+AGbWAxgLPBVrVMn0HeDn7r4+7kAS6IvAbmnTR8BXUZeZTK8CAzO6fmwE3ogpnsSJhh7c1d2XR6uuAFaiUUHao1yevUbgpbiDSDj9XGVPP1ftKOvC2t0XAQ8CB0erDgHuc/cF8UWVPGZ2ITAL6GNm25rZ6VG3EAHcvdbdP0xNhDuma919ZdyxJYm7zwH+Thi6ETMbBPQEHoozroT5BFhnZlulrVsKvB1TPCVBubx9ZjbKzL5oZt2iG87OAy6KO64ksaavHQ30c9WSzGukn6vslfVwewBmNhi4GqgBRgDfd/dVsQaVIGZ2MfCjjNVz3F0taK0wsxo03F6LohbZXwEzCC0/t7i7WqzTmNmuhBsYZwLDgX+5u1rI2qFc3jYzm0T4w3Yl8DRwbfTHrgBpo/DcCFwO/Nrda/Vz1aSlawQMQz9XWSn7wlpEREREpBjKuiuIiIiIiEixqLAWEREREckDFdYiIiIiInmgwlpEREREJA9UWIuIiIiI5IEKaxERERGRPFBhLdIKMzvYzJ41s6nt7Lezmf2Pmf0gj+f+Wb6OJSISBzP7lJndk8/c2Mp5HjSzkwt5jkIxs0PM7Ln2fs9I6egedwAi7TGznYArCU+tmwx8Cvh/7n5HgU/9AbBrB/Z7HxgMVOTjpNFjdevzcSwREQAzmwz8BNgf+Acwj/BAlB7Ab9z90QKc9l1gIHnKjW24EXi9EAeO8vFxwE2EJ6Q+F23qBRwPTHH3mhxO0dHfM1IiVFhLKfgrcJK7v2Jm3YDfFeOk7v62mS3twH4rzWxhHk99APCvPB5PRLo4d3/RzP4M7Ojun0utN7Mjgb+Y2Xfc/cY8n3NNnnMjAGZ2FPBaqqB194fzfY4Ud68H7jazc4An3f3itDge6MwxzWwYcKC73+3ub5nZkjyFKwmgriCSaGa2JbA9sBrA3RuBS4BiPTK0Mc/7dcSBwBN5PJ6ICMDGzBXu/hBwKfBLMxtRgHPmNVeb2WhCC3WxbXbtgMeBBdkcxMx6AbcDfdJW5/P3h8RMhbUkmrt/ArwFPGJme0fr5gP/Z2aDzOxHZvZ/Zna0mS0yszfNbI/U+83sFDO7yMyeNLMfp62vMrNLzexuM7vfzAZG63ub2Q1mdomZ/ZbwNWZWzOy7UUwnmdlkM/trdLwfRTHebWZDo/UrzOyHGYfo6e7ro2OdZ2bfNrNnzOyqrC+giEj7fkvoEnIigJmNN7PLzexWM/tnVMxiZr2iXHaRmf3bzL6QOoCZnW1mPzCz583sG62dqKWcbGYnRDnzNDN7x8yuMrNRZva7KJ/+28wOig7xRWA0cIGZHRPdC/O4mf1n2jnOMLMrzex2M/uDmfUzs4Fm9kMzm2FmB5nZW2b2npmN68wFi67J1KhFGzM7zsx+ZWbXmNkjZtY/+gy/MrPp0e+bhdE1ngicaGbfTjvk4Oh3UZ2ZndHKObuZ2U/NbFp0Djezv6T+fSQh3F2TpkRPwCRgPuGv+tuAYdH6CkKS+oTQ120Q8Bgwm/BH46eBi6N9BxNavY8EtgD+nHb8Z4CfRcu/BKZFy8OBBkLybC/G3xNafYYR+uL1jdb3Ah4FHgBGEX4h1AM/BPoCe0avU/un+o+nPveN0fIA4Htx/1to0qSpdCfgNODDVrYtBq6L8ur9QEW0/i7gnmj5F8Bx0fLpQG20fApwarS8fZSrd4le/x64NFpuLSdvCSwHrgWmAIcAPwcuifY9H5ieFqsD46LlPsCbwGnR68OBu9L2fRC4mdD19djo98Vx0ed8BLiyg9fuyeh3xVXAzwj9rQ9J274EGBstzyb8TuoFXAR8COwF/FeU959MxRvtXwNcE207BXi/lRi+CjwYLQ+PrvNhcf9caWo+qY+1JJ67zzazXYArgLOBo8zsEA99rmuBle5+H4CZXQo8DWxL+CXSaGbfig71KCGBHwkMSVv/ZvTe/sBZ0Xtx98Vm9kEWoVYC1wP/5e5romOsN7NFwDx3/yg6zyLgKQ/9D2cREn4l4SbIIwm/BADWAF82s7eB3xD+qBARKQQnFJtTgK2Bc8wMYCkw0MKLrxEaBQD+QOgKASHXzsjItSOB1zLOcRot5GR3/8TMVgAPuPv/AZjZ+8BqMxtJaOHt32LQ7mvN7OO0Vd+IjptyK/A/0frlhN8X90fnmBHF2VFPeNTH2szuJjTmpJzo7u+b2T6Egrp/lP8/Aua6+wvAC9F7Wzr29Oh3wottxLQHsDb63IvN7BOadymRBFBhLYlmZhXAGA83qfx31D3jYaCa0AKQaXY0HwSMBf7o7ndF666NjnkBodC9NuNcexISYl0nw90e2BcYD7yStj6zj2HDpg3ujVGSTXXLGuKh+wvuPs/MTgN+DZwLfAn4GBGRPDKzQYRv22YT8ubSFvLjcEKLak8Ad28gNAYQvecyd382et3svWlazMkRp3mu/BC4AHgHeB7Ypo2PkP6+7YD/TXs9l9DNZRib5+KNdH7EkpcIf4CkLDeznwL3EGJPVc+Zn6s9G2m9Nvs78Bsz6034t1hLaP2WBFEfa0m6Tf3+ANz9FeBCYMdW9u9JSGLvAguBE9I3mtle0fojzaxPxvqV0cudOhnrs8CvgDujxJcVMxtF2o0wZrYVYUSU7Qk3M97VyltFRHLxn4RvyO4h5Md9opZiYFN+XAJsAI5IW7+LmQ0gI9eaWfeooSJTazm5JTcBNR6GVV2fxWd5n9DCvekUwApCV5e8cffGqPFjUvSHycPA5e4+I5/nyTjndGA64ZvVU4AD3H1Foc4nnaPCWkrBOVHRmbIVoW9cypZm1jNankroX7cc+BNwrJn92symRq0JqwkJsB/woJkdHt0osiOhGJ8F/MjMeppZJaFv8zAzqzCz75nZqa3EaNF0CeGXzy9a2Ja5f+brzxH6A6ZsBxzj7iuBc1p4j4hINjZrCTWz/YHLgK+5ey2hdXgx8LCZHRXdoLhf1EL9Z8LoIceY2cHAf7p7HSHX/reZXWZmBxAaGFLd6NLzX2s5OSW99Xh3YKiZ9QX2AfqY2bbRtnrCzX6pBpb0c9wEHG9mg6PXewG3uPtGQs3TYi42s72juHq1cu26kVEzRX98HEn4lnIwoYvhDsAYoG/ajZGZreIbWoi/eVAt9BeJ/i22AP5JaJXP5g8OKRJ1BZFSMAaYbWYPEb4mM+DMtO0bgYvMbDmhdfccAHd/zMz+m/B14jGEm/9eBzCzzwM3EH5R3Aac7+5uZicRhkKaDdxN+EpvN8IDFU4idDFp9mAaM5tCuClnLOGrwb8Cl0XxPArsDYwzs/GElpQRwElmNoem1pvTgH7uXp3x2X9n4cEOGwk3C4mIZC3KI18kFH+3ALWEG+AGAp919+cA3H1DlB9vIeS6Bwh9qwG+SWiUuJ3wDV1qJI7fEnLfWYQua2e5+8dR4bg3MN7M7mgtJ5vZCYR+xV81sznuvoBwI+XPCUX17cDngV2A94A/ElrXv2xhtKhJwH+Y2cPufr+ZbQfcZ2bPEL71vChqWT8RGBF9vtmEh+WMNLNJUZwnE/qN/zXtuvUg5P5dga2iwtsJfZuPi475MvAU8CKhe8vfCTn9QcINk7ua2THunjru3YSH9ayKCvERwAlmNo/QEk30/sxnNrxPuDnzGKA30M3MHifcwFisIWilHaZ/CyllFh4D+3t3H1ek853v7nrcuIhImYm6pfRy90Q+oMvMzieM0FITve5N+Lbhp6l7cyR+6goipa6lbhaFOVH4yvSlYpxLRESKJ7o5c48EF9W9gYuBcWmrK4HVKqqTRV1BpGRF/dtOJXy1d6K731PgU86NvqIUEZHy0t3dfxN3EK1x93UWHqt+XXRj+0eEbisXt/1OKTZ1BRERERERyQN1BRERERERyQMV1iIiIiIieaDCWkREREQkD8rm5sXKykofN25cVu+pra1l6NChhQmoiMrhc+gzJIM+Qzxmzpy5xN1LK+gcdSZnQ2n++3aEPldp0ecqLfn+XG3l7LIprMeNG8eMGdk9SbSqqirr9yRROXwOfYZk0GeIh5nNjzuGYutMzobS/PftCH2u0qLPVVry/bnaytnqCiIiIiIikgcqrEVERERE8qBLF9bTpk2LO4S8KIfPoc+QDPoMknTl+u+rz1Va9LlKSzE/V9EfEGNmhwJXAF9Ie979ccA1wBbAHcB57r4x2lYJXAksB3oAF7h7feZxq6qqvBz7BYlI+TOzme5eFXccLVHOFhFprq2cXdQWazMbDvQHpqStGwscA5wAnAt8FfhW2tvuAW5z9wuAt4HLixawiEgXppwtIpKdohbW7r4YeCBj9dbA6e4+091vB24ADgQws72BHdz9hWjf6cA5ZtavWDGLiHRVytkiItkpeh9r7uQ4vgAAIABJREFUd2/MeP3v1FeIkQXA+9HyQcC8tH0XRIv5+cp00SLYe2846KC8HE5EpNwkJWfXzapjxp4zmPOVObkeSkSkYJI4jvVk4GfR8lbA0oztdcCozDfV1tZSVdWUu6dNm9Z+Z/VeveCFF2DgwFziFRHJWnV1NdXV1amXlXHGkqOi5OzGdY2semkV3Xp26XvuRSQmHc3ZiSqszWwCUOvur0SrHFibsVsvYLMbYYYOHZr94N9bbBHmK1ZAQwNUVGQZsYhI56QXkma2JOZwOqWYObtbr1BQN65vbGdPEZH862jOTsyf/mZWAXwN+F7a6gXA4LR9DOgHfJSXk1ZUNLVWr1yZl0OKiHQFxc7ZKqxFpBQkprAGvgP83N3Xp617DNgh7fU4YAOQvzGaBg0K82XL8nZIEZEuoKg523oZoMJaRJKt6IV11IIBYGnrLgRmAX3MbFszO93MJrj7DGChme0c7XoEcL27b8hbQIOjxpXly/N2SBGRcpGUnJ1qsfb1xX32gohINorax9rMBgCnRi9PM7NfE75K/FHGrnPc/bZo+QvAJWZWQxhP9aK8BqUWaxGRFiUpZ6sriIiUgqIW1u7/v737jpOrrvc//vps75tNNp00TAWSEAhV6XApUv2BKOgVFYId9aqooChyLQiKF0ENiB0rKkXFS5NyKZJAEkoIgZBCGpuy2Z5t398f3zO7szW7O2f2zOy8n4/HeczMOWfO+Zwpn/3sd77ne1wt8KNgirkumPp6zgbgw0kLSi3WIiK9SqWcrcJaRNJBKvWxjkasxVqFtYhIyrI89bEWkdSnwnr0aH+7s/vQqyIikiri+1g7p37WIpKaVFhPmOBvt26NNg4REemTZRmW41utXYsKaxFJTSqsJwUXBFNhLSKS0jTknoikOhXWEyf62y3hXHNGRESSQycwikiqU2Eda7FWYS0iktI6CusmFdYikppUWMd3BdEJMSIiKUsXiRGRVKfCurTUTw0NsGtX1NGIiEgf1MdaRFKdCmuAmTP97dq10cYhIiJ9yi7MBtQVRERSlwprgFmz/K0KaxGRlJVd4gvrtrq2iCMREemdCmuA2bP97auvRhuHiIj0qaOwrlVhLSKpSYU1qMVaRCQNZJeqxVpEUpsKa4A5c/ztyy9HG4eIiPRJXUFEJNWpsAaYPx+ysnxh3dgYdTQiItKLWIt1a21rxJGIiPROhTVAUREccAC0tcHKlVFHIyIivVCLtYikOhXWMYsX+9tly6KNQ0REetXRx1onL4pIilJhHXPYYf72ySejjUNERHqlFmsRSXUqrGNOOMHfPvywLm0uIpKCckpzALVYi0jqUmEdM3cuTJwI27fDSy9FHY2IiHSTXR6cvLhbJy+KSGpSYR1jBied5O8/9FC0sYiISA+5lbkAtOxsiTgSEZHeqbCOd8op/va++6KNQ0REesgbmwdAyw4V1iKSmlRYxzvzTMjJgUcegR07oo5GRETidLRYV6mwFpHUNOyFtZmdYmbPmNn0uHmFZnazmX3TzH5kZmUDWRa60aN9d5C2Nrj77qTtRkQkXaRSzs4ZlQNZ0FrdSntLe1ibFREJzbAW1mY2HigBDu+26FbgCefcl4G7gZ8McFn4zj/f3/7qV0ndjYhIqku1nG3ZRu5o32rduksnMIpI6hnWwto5tx2faDuY2STgIuDeYNaDwNlmNr2/ZUkL8t3v9ldifPRRWLMmabsREUl1qZizY91Bmquaw9qkiEhohr0riHOu++93xwM7nHMNwfJW4A3guH0sS46yMnjve/39225L2m5ERNJBquXs3HFBP+vt6mctIqknFU5enAzs7DavFpi0j2VdVFVVsXjx4o5p6dKlQ4/o8sv97c9+BvX1Q9+OiEgfli5d2pGvgMqo4xmESHN2/pR8AJo2NQ0ybBGRoRtozs4ZvpD65IDGbvPygVhzRH/LOowdO5Zly5aFE9HixXDkkfD003D77XDFFeFsV0QksGTJEpYsWQKAmaXTMESR5uyCqQUA7N24d9DPFREZqoHm7FRosd4MVHSbVwps2cey5DGDL37R37/hBmhWXz4RkUCkOTvWYr13kwprEUk9qVBYPwJMMrNCADPLBaYCj+5jWXKddRYccAC8+Sb84hdJ352ISJqINGfHWqybNqoriIiknijGsbbYXQDn3DbgPiC4njgnA3c55zb3tyzpgWZlwVe+4u9/9atQV5f0XYqIpJpUy9n5U9ViLSKpa7jHsS4FPhI8vMTMxgb3PwqcY2ZXAWcAS+Ke1t+y5LrwQjj8cNi2zXcJERHJIKmYswumBS3WbzTh2lyYmxYRSZg5NzIS0+LFi11oJy/Ge+IJOOYYKCiAF16AmTPD34eIZDQzW+6cWxx1HMMpkZz95OQnad7SzBHrjqBwRmHIkYmI9K+/nJ0KfaxT2zveAe97HzQ1wWWXwQj5R0REJF0VzS0CoGF1Q8SRiIh0pcJ6IL7/fRg7Fv71L0hkfGwREUlY0TwV1iKSmlRYD0RlJdx8s7//2c/C6tXRxiMiksE6WqxfUWEtIqlFhfVAvfvd8P73Q0ODv9/Y/RoIIiIyHIoPLAagbpVGaxKR1KLCeqDM4NZbYfZsePFF+PjH1d9aRCQCJYtKAKhbUUd7c3vE0YiIdFJhPRglJfCHP0BhIfzsZxqCT0QkArmjcimcXYhrdtS/UB91OCIiHVRYD9bChfDLX/r7V14Jf/5ztPGIiGSg0sNKAaj5d03EkYiIdFJhPRTnnw/f+pbvCnLRRfDII1FHJCKSUcoOKwOg5hkV1iKSOlRYD9WVV8JHPgJ798JZZ8HTT0cdkYhIxig/phyA6n9VM1IudCYi6U+F9VCZwS23+IvH1NfD6afDqlVRRyUikhFKDi4hZ3QOezfspWldU9ThiIgAKqwTk5XlT2I891yoroaTT4aVK6OOSkRkxLMsY9QJowDY/dDuiKMREfFUWCcqJwd+9zs47TSoqoLjj4ennoo6KhGREa/ipAoAdj+gwlpEUoMK6zDk58Nf/wrnnedbrk85BR54IOqoRERGtNGnjwZg1/27aGtqizgaEREV1uHJz/djXL///Z19rn/yk6ijEhEZsQqnF1JySAltdW1qtRaRlKDCOkw5OfDzn8MXvgBtbX7UkM98xt8XEZHQjX3XWAB2/HlHxJGIiKiwDl9WFnznO/6kxtxcuOkmOOMM3/9aRERCVfmuSgB23L1D3UFEJHKhFNZmdkYY2xlRLrkEHnwQKivhf/8XDj4YHn886qhEREZUzi6eV0zJwSW07m5l5907ow5HRDLcoAprMzvNzJ4ws9fMbF0wrQfuSk54ae7YY2HFCnjHO2DLFjjhBH/FRnUNEZFhkCk5e8KHJwCw9adbI45ERDLdYFusfwX8GbgM+GAwXRLMl95MnuwveX7llb6g/vKXfcH92mtRRyYiI19G5OzxF43H8o3dD+6mcX1j1OGISAbLGeT665xz3+s+08xeCSmekSknB779bT/G9Yc+BE8+CQsX+r7YH/uY75ctIhK+jMjZuaNzGfv/xvLWnW+x+YebmXnDzKhDEpEMNdiK7qNmdn4v888OI5gR77TT4MUX4eKLoaEBPvlJX2y/+GLUkYnIyJQxOXu/z+wHwNalW2nd0xpxNCKSqQZbWL8duN3M2uKmduBHiQZiZvPM7BYz+4yZ3WpmBwfzC83sZjP7ppn9yMzKEt1XpEaPhl//Gu66C8aO9Sc0Llrkh+irq4s6OhEZWTImZ5ctLmPUCaNoq21jy9Itw7FLEZEeBltYXwssAWYB+8dNXwshll8B33TOfR/4NvDbYP6twBPOuS8DdwMj46or73oXrFkDH/2o73v93e/CvHn+IjPORR2diIwMGZWzp3x+CgCbbtxEW71OEheR4TfYwvofwN3OuXXOuQ3BtB64LYRY5gGlwf1GoNzMJgEXAfcG8x8Ezjaz6SHsL3oVFXDrrfDMM3DoofDmm3DhhXDUURqaT0TCkFE5e/Rpoyk9vJSW7S28edObw7FLEZEuBltY/w74kpkdGzcdB3wphFj+hP/JshR4H/BJ4Hhgh3OuAcA51wq8ARwXwv5Sx2GH+eL6xz+G8eP9/WOPhXPOgdWro45ORNJXRuVsM2P/b+8PwMbrN9Kys2U4disi0mGwhfVXgE8DvwR+ETddHkIsHwOagGeBOufcXcBkoPuI/7XApBD2l1qys+Hyy/0wfF/7GhQXwz33wEEHwUUXwUsvRR2hiKSfjMvZFSdUUPEfFbTVtPHGNW8M125FRIDBF9bfBaY456Y752YE03TgXSHEUojvo3cncJOZnQw4/E+M8fKBHs0QVVVVLF68uGNaunRpCCFFoKQErrnGF9gf/agvuH/7W5g/Hy64AFatijpCEQnB0qVLO/IVUJmk3WRkzn7bDW+DbNhy6xZqltWEtl0RyVwDzdnmBnGinJkd5Zx7qpf5H3LO3TGkSDu3cT/wHudctZl9BfgUvqXlGufc7Lj11gbz7ox//uLFi92yZcsSCSE1bdrkx7u+7TZobvbzzjwTPvtZP1SfWaThiUjizGy5c25xErabsTn79c+/zqYbNlFyaAmHPnMolq1cKSLh6C9nD7bF+mO9bPwC4KahBBa3jUpgoXOuOph1HVADbAAmmVlhsF4uMBV4NJH9pZUpU+CHP4R16+CKK6CgAO67D048EQ45BH71q86CW0Skq4zN2dOumUb+lHzqltex8fqNw717EclQgy2sRwctE5hZkZndDtyO72OXiF1Ak5lNjpu3E1gF3AecFMw7GbjLObc5wf2ln8mT4aabYMMG+PrXYdw4WLEC/vM/Yfp0uO462Lo16ihFJLVkbM7OKclhzm1zAFj/1fXULq8d7hBEJAMNqisIgJkdCVwMnApsA/4T2OOc251QIGYL8a0ry4HxwGPOuUfNrAK4HlgPTAC+5JzrcSWVEdsVpC9NTXDnnfD973deuTE7G84+258EecopulS6SJpIVleQYNsZnbPXXrGWzf+zmcI5hSxevpjs4uyk71NERrb+cna/hbWZHdvHosOBE4AbgHbgXOfcZxINNBEZV1jHOAcPPQQ/+hHcfbe/2Az4VuzLLoMPfQgmTIg0RBHpX1iFtXJ2T22NbSw/bDkNLzUw9sKxHPDbAzCdmyIiCUiksH4afxGAanwy7nU1YKJzLj/RQBORsYV1vK1b4Y47/ImOGzb4ednZcNpp8P73+9bswsJoYxSRHkIsrJWze1G/up7njniOtto29v/2/ky9cuqw7FdERqZETl68FpjlnJsWN1RT92k6/mdGidrEiXDVVfD66/CPf8C55/pRQ/72N3jPe3zL9aWXwmOPQXtff3NFJI0pZ/eieF4x8349D4B1X1rHjnt3RByRiIxU/RbWzrm/O+fe2tdGnHN/Ci8kSVislfovf4HNm+EHP4DFi6GmBn76UzjuONh/f/jKV3RlR5ERRDm7b5VnVzL92ung4OULX2bPk3uiDklERiCd3TbSjRsHn/oUPPssvPwyfOlLfgi/DRv8SCIHHOAvPnPttfDKK1FHKyKSNNOunsbESyfS3tjOC2e+QP1L9VGHJCIjjArrTDJvHnzzm7B+PTz8sD+xsaLCjypyzTV++YIF8I1vwJo1UUcrIhIqM2PWj2ZReW4lrbtbWXnySupfVnEtIuFRYZ2JsrLghBN8t5Bt2+Dvf4cPfhBGjYIXXoCvfhXmzvUt2VdfDf/+t/pki8iIkJWTxbw75zHqxFE0b2tmxfErqFvVYzRAEZEhSbiwNrNFZvZrM/svM5ueeEgyrPLy4PTT/Wgi27f7Ex0/8AEoL/ct2f/933DEEf4CNUuWwL33QkND1FGLyBApZ0N2YTbz75tPxakVtFS1sOKEFdQsq4k6LBEZAcJosf4IcLJz7kbg6NilbCUN5eXBGWfAz3/ui+x//hM+/nHfJ3vbNj+M39lnQ2UlnHMO3H67PzlSRNKJcja+uD7orwcx5swxtO5qZcXxK9hxn0YLEZHEDPrKiz02YJYLjHPObTazLKDEOTfs//prHOskcg5WroR77vEt1t1f5wUL/Cgkp58ORx/tC3QRGbBkXnmxl30pZ8dpb25nzaVr2P6r7ZAFs26exeSPTd73E0UkYyUyjnVvGzvdzC4O7h+OHzN1M4Bzrj2KBC1JZgYHH+z7Xj/7LLz5Jvz4x3DmmVBUBKtWwfXX+37blZVw3nmwdCls3Bh15CIZTzm7f1l5Wcz9xVymfXUatMPaj69l7SfX0t6s80pEZPAGVVib2X8DfwTOA3DO/Rv4TzO7KAmxSaqaPBkuv9y3Xu/cCQ88AJ/9rB+6r7YW/vpXv3zaNDjwQPjc5+DBB2Hv3qgjF8koytkDY2bM+PoM5vxsDpZrbP7hZlYcv4KmN5uiDk1E0sxgW6yPBvYDHo+b9w/g26FFJOmloABOPhluvBFeeskP5feTn/irPpaU+LGzb7wRTjkFRo/23UVuvBFWrNBIIyLJp5w9CBMvmciixxeRv18+NU/VsPyQ5ex+eHfUYYlIGhlsYf2Ac64aiO+YfX6I8Ui6mzbNjx7yl7/41uxHHoEvfMH3w25ogPvv9y3YixbB+PHw7nf7biOvv+77cotImJSzB6nsiDIOfe5QKk72I4asPHklr1/5Ou171RAgIvs2qJMXzew9wP7AZHwLyAX4nxg/6Zy7JSkRDlCqnAgj/di61V+Y5sEH4aGHYNOmrsunTYOTTvIt4Cee6AtvkQyQrJMXlbOHzrU51l+7ng3XbYB2KF5YzAG/OYDiA4ujDk1EItZfzh70qCBmdiRwCTAV2Ar81jn3YKJBJirVk7R04xy89povsB980Lds79rVdZ25c+HYY+GYY/zt1KnRxCqSZMkcFUQ5OzF7ntrD6vetpmldE5bv+2Lv99n9yMrV9dVEMlWohXUvG88GjnfOPZTQhhKULkla+tDe7vtdxwrtJ57oeSGaadN8gR2bZs3yI5aIpLlhHm5POXuQWmtbee0zr7Htp9sA33o957Y5lB1WFnFkIhKFIRfWZnYPULKP7VcCO5xzJw49xMSlU5KWAWhpgeeeg8ce89Pjj8OePV3XGT++s8g+5hg46CDIzo4mXpEEhFVYK2cn165/7uLVj7xK0/omyILJn5zMjGtnkFOWE3VoIjKMEimsrwNygDVAO3AZ8DwQnw1nAKOdc58KLeIhSMckLYPQ1uYvsR4rtB97DN56q+s6JSVw+OFw1FF+OvJIGDMmmnhFBiHEwlo5O8na6ttY//X1bPreJmiD3PG57P/f+zPhkglYtn5BE8kEiRTW5cE61cHja51zX+1lvbudc+eEFfBQpGuSliFyDl591RfYjz4KTz4Jb7zRc73ZszsL7aOO8uNqq1VbUkyIhbVy9jCpfb6WtZ9YS82T/vo6JYeUMPOmmYw6ZlTEkYlIsvWXs/v9/co51+23d8abWbZzri1u48cDhyccpchgmMGcOX667DI/b9s2eOqpzmnZMl98v/oq/OIXfp3S0s5W7SOOgMMO0+gjMmIoZw+f0kWlLHpiEW/97i3WfWEddc/VseLYFVSeV8mM62ZQfIBGDxHJRIPtGPYnYKWZPQHUAvOAU4EeLSIiw27CBH859fPO84+bm2Hlyq7F9oYN/gTJh+LO25oyxRfYsenQQ2GUWp1kRFDOTiIzY/x7x1N5TiWbvruJjd/ZyI6/7GDHX3cw/n3jmf616RTuXxh1mCIyjIYy3N4k4GPAAqAG+KNz7u4kxDYo6f6zogyTrVs7i+xnn4Xly6Gurud6s2Z1LbYXLYKiouGPVzJCkofbU84eJnu37mXDf29g69KtuBaH5RgTL53I1C9OpWBaQdThiUhIkj3cXhawyDm3PKENdW7vSOAdwGvAE0A9cD2+taUCuNI5V9P9eSMxScswaGuDNWt8kR2bVqzwrd3xsrN9/+zDDoPFi/3tQQdBfn40ccuIMszD7SlnJ1njG42s//p6tv9qO7SD5Rjj3zeeKVdOoXiuuoiIpLvQCmszu6OX2WOBNufcuUOML377HwTe5py7Om7ez4D7nXO/N7PTgA84597b/bkjOUnLMGtu9iOQxBfbL73ki/B4OTm+2D7kEN+ifcghsHChH51EZBCSeOVF5ewI1b9cz4ZvbuCt377lx2gxGPv/xjL1S1MpPaQ06vBEZIjCLKwfBx4C4p/0duBx59w3EgzyOOArwCkuCCr4CfMNoMI512BmOcAe4EDn3Pr452dCkpYINTT4luxYob18uW/p7v79MfMjkcQX24sWwejR0cQtaSGJhbVydgpofL2RjddvZNvPt+Ga/Vsx6vhR7Pfp/Rhz5hgN0yeSZsIsrOc459Z0mzcXuMQ598VEgwSeBAyYCXwDmAZ81zk3OW69F4N5v4h/fiYlaUkRdXWwapW/kM3zz/vbF1+E1tae606b1rPYnjhRV44UIKmFtXJ2Ctm7ZS+bvreJrT/ZSlud/wWsYEYBkz8xmYkfnkhOuS40I5IOkt3HegbwrHOuMoFtzAFWAwc751aZ2efxZ61fB1zsnFsQt+5TwD3OuW/Fb2PatGlu7NixHY+XLFnCkiVLhhqSyNDs3eu7jcQX2ytXQmNjz3XHj/cF9sKFndPs2b6LiYx4S5cuZenSpQAsX758g3Nu+nDsVzk7eq17Wtn6s61svnkzTeuaAMgqzmLCJROY9JFJlByk7mQiqWagOXuwLdZv0PUnxWxgPPAr59xlQw3WzN4FLI0lejMrBKqAa4H/55w7Im7d54A7nXM3xG8jE1s/JE3ETpCMFdqx2+6XaAcoKPAnRcYX2wsWaPi/ES6JLdbK2SnMtTl2/m0nb/7gTaofru6YX3ZUGRMvm8i4d48ju1gXtBJJNUO+QEwvHgDupDNRtwPbnXOvJhBfLI6OWJxzjWa2Fv9HoKLbuqXAlgT3JzJ8srPhgAP8dPHFfp5z/kqRzz/vW7Rj04YN/sI23QuOadO6FtsLF8L++0NW1vAfj6QT5ewUZtlG5dmVVJ5dSd0LdWy5dQvbf7OdmqdqqHmqhtc+/RrjLxrPxMsm6mRHkTQx2BbrY51zj3WbNwd/hvlrQw7C9/lbDUxwzm0P5j0L/BC4BRgbJO5coA7Y3zm3OX4bav2QEaG62vfbji+2X3wRmpp6rltSAvPndy2258/XqCRpKIkt1srZaaatvo23/vAWW2/bSs1TnaMUlhxcwvj3jWfcRePIn6hhPkWiFGaL9WnAY93mbQceARYNITYAnHOvmNk/gPOBW8xsFJAH/AY4HTgJuA84Gbire4IWGTFGjYJjj/VTTGsrrF3btdhesaLrxW5izOBtb/NF9sEHdxbcU6boRMnMpJydZrKLs5n4wYlM/OBE6l6sY+vtW9n+y+3UraijbkUdr3/hdSpOqmD8+8dTeV4lOSU6J0MklQyoxdrMPoE/KaWv36IecM6dllAgZpXAD4BlwFTgNufcy2ZWgb/YwHpgAvAl51yPS+Wp9UMyTlVV12J75UpYvRpaWnquO2pUZ4t2bDroICjVz8upIOwWa+XskaWtqY1df9/F9l9vZ+d9O3Et/u92VlEWledWMv6i8VScUkFWnrqGiQyHUEYFMbMFwAeA/+m2qNk5tzWxEBOnJC2Cv7jN6tU9C+4dO3pff8aMrsX2/PkamSQCyegKopw9MrXsaqHqj1Vs//V29jzReQJ0dnk2ledUMvaCsYw+ZTRZ+SqyRZIlzHGsc51zPZrDzGyqc25jAjEmTElapA/O+W4jK1fCCy90TqtX97x0O/jLtM+b17PgnjRJ3UmSJIl9rJWzR7DGdY1sv3M7VX+oov6F+o752WXZVJ5dydjzx1JxagXZBRpZRCRMQy6szexzwN+dcy8Hj6/Cn/XdZTXgBOfc8eGEOzRK0iKD1NLi+26vWtW14F6/vvf1R4/uWmgvWOC7k+hkyYSFVVgrZ2euhjUNVP2piqo/VVG3orPnTXZJNqNPG82Ys8cw5owx5I7JjTBKkZEhkZMXZwDxA+geAhQD2+gcvsnw/ehEJJ3k5nYOA/ie93TOr6nxI5HEF9svvAC7dsGjj/opXnx3kgUL/O2sWepOEg3l7AxVNKeIaVdNY9pV02hY20DVXVVU/bGKuufqOgpusqD87eWMOWsMlWdXUjSnKOqwRUacwXYFOQR40TnX3G3+AbEWkqio9UMkiZyDzZt7Ftv76k6yYEHntHAhjBs3/LGngSR2BVHOznCN6xvZee9Odt67k+p/VXec+AhQOKuQMWeNYcxZYyg/ulwnP4oMUJh9rD/gnPtFt3lzgNnOuXsTCzMxStIiEWhpgVdf7Vlw99WdZMKEziI7djt3rm89z2BJLKyVs6VDa00ru/65i5337GTn33fSuqu1Y1l2STajThzF6P8YTcWpFRTNVGu2SF8SLqzN7Bj8z4dLgKXdFo8G/sc5NzXRQBOhJC2SQmLdSVat6jrV1vZcN9YlJb7YXrAgo1q3kzDcnnK29Ku9tZ2ap2o6iuyGlxu6LC/Yv4DRp45m9KmjGXXiKHJK1bVLJCaMwno+8Et8/73d3RY3A79xzl2baKCJUJIWSXHt7f6S7bEhAGNXmHz99d7Xj7VuL1oEhx4KhxziL+M+AkcmSUJhrZwtg9K0qYnd/7ubXf/cxe4Hd9O6u7M123KMsqPLqDipglEnjqLs8DJ1G5GMFtY41uXABc6528MMLixK0iJpqq7Odx+Jv5T7qlV+fnfl5b7QPuSQzmn2bMhO7+HEkjSOtXK2DIlrc9Q8W8Puf/pCu+aZGmjvXJ5VlEX528sZdcIoRp0witLFpWTlqNCWzBFaH+t+dvAO59wTCW8oAUrSIiNIe7vvp71yJTz/PCxf7qft23uuW1zsL98eX2wfcEBajUqSrD7W/exPOVsGrKW6heqHq9n98G6qH6nu0W0kuzSb8mN8oV1xQgUlB5dg2SPvlyWRmETGsb4H2NcgtQYc6JyLtEOkkrRIBti6FZ57ruu0sZfrnBQU+G5c7Z8kAAAgAElEQVQk8cX2QQf50UpSUIjjWCtnS9I1b2+m+l/V7H7EF9qNrzZ2WZ5dnk350eWUv72c8neUU3p4KdmF6f2rkki8RMaxXgNUA5vpHAO1x/aBM4cenojIAE2cCO98p59iqqp8q3Z8sf366/Dvf/spJifHF9eHHOL7bB92mG/pHlkjkihnS9Lljc9j3IXjGHeh/99s7+a9HUV29SPVNL3RxK5/7GLXP3YBYLlG6aGllL29jPJ3+II7b2xelIcgkjT7arGuBGqdc3v73YjZGOfczrCDGwy1fohIh+pqWLHCdx+JFdtr1vjxuOMVFsLixXDUUZ3T+PHDHm6ILdbK2RK5pjebqPm/GvY8sYc9T+yhbmVdj3/zCmcXdhTZZUeWUTS3CMtS9xFJD8PRx/rSqE+QUZIWkX7V1fk+2889B8uWwdNP+zG4u5sxA449Fk4+2U8Tkn+Rwgj6WCtny7Bp3dNKzdM17Pk/X2jXPF1De2N7l3Wyy7IpPayUsiPKOqa88WrVltSUtMLazKYCHwc+4ZwrHvKGQqAkLSKDtnOnL7CfegqefNJ3Hamv77rO/Pm+wD71VDjhBMgL/4/9cBXWytmSCtpb2ql7vs63aD+5h9pnatn7Zs8fWfKn5XcptEsOKVFfbUkJoRbWZpYNnA1cBpwCbA+2MznRQBOhJC0iCWtt9UP/PfwwPPggPPooNMadmDVqFJxzDlxwgS+2QzoZMpmFtXK2pIO9W/ZS80wNNc/UUPtMLTXP1tBe37VV23KM4gXFlC4upfRQPxUfVExWvob6k+EV1jjWM4BLgQ8CY4ANwKeAfwJHOeeeDCfcoVGSFpHQ7d3rW7IfeADuvddfTTJm7Fi47DL4yEdgypSEdpOkcayVsyVtuTZH/cv1nYX2MzXUv1TfZTxt8CdGFh9UTOmhpZQcWuKL7fnFZBeoZVuSJ6HC2swuwLd0nAi8BtwO/AL4snPuMyHHOmRK0iKSdK+8AnfdBb/7XWeRnZ0NH/sYfO1rMHr0kDYbZmGtnC0jVWtdK7XLaqlbXkftc7XULq/1Q/11K2Msxyg6sKijVbvkkBJKFqobiYQnkeH2ABrxl8B9ErjUOfdqsNHEz3oUEUknc+fCVVfBl7/sW7Jvvhn++Ed/++c/w9//7sfPjpZytoxIOSU5VBxfQcXxFR3zWmtbqVtRR+3yoOBeXkvDKw3Ur6ynfmU92+7Y5lfMgqI5RZQsLKF4YTElC0soWVBC3qQ8zDQaiYRnMF1BpuJbQaYDdwPHOec+GSwzF8bwIglQ64eIRGLVKrj8cn8S5NSpfqSRQfa9TlJXEOVsyUitda3Ur6yndnkttc/5grv+5Z7dSAByxuT4IjuYihcWUzxP/balf2GfvJgDnAtcDLwI/BT4D+fc0kQDTYSStIhEpr4exo2DhgY/fvbChYN6epJPXlTOlozX1thG/Uv11K+qp25lHXUr66hfWU9rdWuPdS3HKJoX17q9oITi+cXkTVDrtniJdgXpwjnXCvwJ+JOZzQQ+g28ViTRJi4gMu9ZW3/3juut8UT1mDLztbVFH1YVytghkF2ZTtriMssVlHfOcc+zdtLdLoV23so7G1xqpf6Ge+hfq4ded28ipyKH4wGKKDyqm6MAiig8qpvjAYl1FUroYdGEdzzn3GnCFmS0PKR4RkdS2aRM8/jg89BDcfbcfCxtg0iTf37qkJNr4+qGcLdLJzCiYWkDB1AIqz6rsmN9a10r9i/UdhXbdqjoaXmqgdXdrx9Uk4+WOy/UFd3zRfWAxuRW5w31IkgISKqxjnHO/DGM7ZpYHPAtc4Zz7l5kVAtcDtUAFcKVzriaMfYmI9Ms5X0SvXOmnFSv8FRs3bOi63uzZfsi9Sy+F0tJoYh2ksHI2KG/LyJNTkkP5keWUH1neMc85R/OWZt+d5MV6f/tSPQ0vNdDyVgvVb1VT/Uh1l+3kTcrraNUuPrCYogOKKJpXRO4oFdwjWSiFdYg+jz/RJuZW4H7n3O/N7DTgJ8B7owhMREaoPXtg7Vp/0mH32z17eq5fXg5vfzsccwycdRYccABkdr9L5W0Z8cyM/Mn55E/OZ/R/dA6r6dp9d5KOYju4bVjdQPOWZpq3NLP7f3d32Vbu+FyK5xVTNNcX2rHb/P3y1Yd7BEiZwtrMjga2AruDx5OAi/CX3wV4ELjLzKY759ZHEqSIpJfWVti2DTZu9C3NGzf2nKqr+35+ZaU/ETE2HXwwHHigH7talLcl41mWUTCtgIJpBYx555iO+a7N0fhGIw0vNXQpthvWNNCyvYXq7dVU/6tr7skqzvJF9tyiLoV34cxCsvI0Skm6SInC2syKgfOdc581s68Gs48HdjjnGsCfgGNmbwDHAesjCVREUkNTE2zd2ve0bZu/fest36WjP0VFMHMmzJrlu3XE344dm+mt0X1S3hbpm2UbRTOLKJpZROU5nf23XbujaWMTDa80+EL7lYaO+y1VLdQtr6NueV3XjWVD4dsKu7RwF88rpnBOobqVpKCUKKyBK4Fvd5s3GdjZbV4tMKm3DVRVVbF4cefIJ0uWLGHJkiVhxigiyeCcbzXesQOqqvwUu9/XbV3dvrcbM24cTJvmx5iO3cbfHzMmkuJ56dKlLF3aMTBHZX/rpqiE8rZytmQiyzIKpxdSOL2QMaeN6bKsZWdLR6Fdv7q+o+BueqOJxlcbaXy1kZ33dP165Y7NpXB2IUWzi7rcFr6tUFeaDNlAc/agx7EOm5mdCuQ75+4JHq8HLgEWAxc4546IW/c54E7n3A3dt6MxUUUi5Jwfbm737q7Trl095/U2v7XnWLL9ys2FCRNg4sTO2/gpNm/8eL9uikvmONbJEEbeVs4WGZi2pjYaX23s0rpdv7qexlcbaW/s5ao3AAb5U/N7FNxFs4somFaAZeuXuESEOo51EnwOWBTXYb8Cf5Ww7wX345UCW4YvNJEMECuK9+yBmho/xe7v67a6urNQbmkZegxlZb7bRWXlwG7LytRFI1rK2yLDJLsgm5IF/hLs8Vy7Y++Wvb7ofrWhy23jukb2btjL3g172f1A15MnLc8ofFthry3deeN1EZxEpUJhfREQf/3fp4DPAs8AnzezQudco5nlAlOBRyOIUSS1tLf7q/3V1vpuEXV1vd+P3fZXHNfUQFtb4jEVFEBFhZ9Gj+68P5B5ebrAQppR3haJmGUZBfsVULBfARUndv1/tr2lnaY3mnoU3A2vNtC8udn3717dwM5uPbeyS7MpnBXXpWRWIYUz/ZQ7JldF9wBEXlg756riH5tZG1DlnNtoZvcBJwH3AScDdznnNkcQpsjQtbf7FuH6+sEVw/0tb2gIN8aCAj+MXFnZ4G7LyzuL44KCcGOSlKW8LZLasnKzKJpdRNHsoh7LWutaaXytsWdL95pGWqtbqXuujrrnep7HkjMqp6PI7j7ljlPRHRN5Yb0PHwWuN7OFwARAZ7ZI+NrbobHRF73xBXD8/X097m9ZU1Ny4i4u9lf5Ky31t/3d769ALitTi7GESXlbJIXllORQenAppQd3vaCVc46WnS00rg2K7jUNvgAPptbqVmqX1VK7rLbHNrNLsn2RPatn0Z03MbO6l6RcYe2cmx53fzdwWXTRSOTa2nzR29Dgp9j9fc3rq+jtrQAOu/W3N8XFfli34uLOgncgBXFf94uLIUvjmkpqUN4WSX9mRl5lHnmVeZQfVd5lmXOOlqqWLoV242uNNK5tpGFtA2172qhbUUfdip4t3VlFWb5Pdy+Fd/7kfCxrZBXdKVdYSxpwzrfCDqXgHexzmpuH55gKCzsL39gU1uPCQp1oJyIiacvMyBuXR964PMqP7ll0t+5q7Vl0v+aL7tadrdS/UE/9C/U9t5tvnUV3t8K7YEp6jl6iwnok6d6621sx29f8wRa8w8Wss+gtKuq839+8wsLBFb5FRWr9FRERGQIzI3dMLrljcik7oqzH8pbdLTS+3rPobnytkZbtLTS83EDDyz3rCss1CvYvoHBmIUWziroU3vlT88nKSc2/2yqsh0OsD29/Re5Ai+H+1tm7d/iOKT+//8J2MIVwf/Py89XaKyIikqZyK3LJXZxL2eKeRXdrbe8t3Y2vNdK8pZnGNY00rmlkF7u6PM9yjYIZBV2K7VjxnT8t2qJbhfW+OOcvi7x+Pezc6cfr7T7V1PRfCCfr5LXuzDqL0u5TfMHa3+OBFLuFhWrhFRERkYTklOZQuqiU0kWlPZa11bd1bele29m9pHlzc8fVKLuznJ5Fd+GsQkrml5A/Ob/H+qEfU9L3kG6qq+Huu+HRR+HZZ2HdunC6PvRWzA60CB7oumrdFRERkREgu7j3C+MAtDX0XnQ3rm1k75t7/eO1jfCPrs8rP6ac6V+fTsUJ3a9jFR4V1jENDXDNNXDLLb71Od7o0TB9Oowb5+93n0pLO/vq9jYVFKjgFREREQlBdlE2JfNLKJnfS9Hd2EbTuiYa1jZ0Ft1rG6l5uoY9j+9h5YkrmfXDWUz++OSkxKbCGvwV6E45xbdQA5x0ErzznXD00TBnDowaFW18IiIiIrJP2YXZFB9YTPGBxV3mt9a2sumGTWy4dgNrr1hL+THlvbaGJ0qFNcAnP+mL6hkz4A9/gMWLo45IREREREKSU5rDjK/PoGVnC1tu2cLG6zdywK8PCH0/OgPthRfgV7/y3TXuv19FtYiIiMgINeW/pgCw4887aG9uD337Kqxvu83ffvjDMHt2tLGIiIiISNIUziikaF4R7Y3t1D7X8/LsiVJhff/9/vbii6ONQ0RERESSruxIP6Z23XM9L8GeqMwurLduhbVr/agehx0WdTQiIiIikmSFswsBaHy95zjYicrswvqll/ztwoWQo/M4RUREREa6wplBYf2aCutwxQrrA8I/K1REREREUk/+JH8FxubtzaFvO7ML6/Xr/e3MmZGGISIiIiLDI3dMLgAtO1pC33ZmF9Zbtvjbycm5+o6IiIiIpJbcSl9Yt+5sDX3bmV1Yb93qbydNijYOERERERkWOaNywKC1upX21nDHslZhDTBxYrRxiIiIiMiwsGwjuzgbgPZ6Fdbhqa72txUV0cYhIiIiIsMmq8iXwG2NbeFuN9StpZuaGn9bXh5tHCIiIiIybLIKfQnc3qAW63A0NUFzM+TlQX5+1NGIiIiIyDDJLgq6gjSO0MLazN5lZm+Y2U4z+4GZ5QTzK83sJ2b2HTP7npnlhrLDWGt1WVkomxMRySTDnrNFRELU0RWkYQR2BTGzqcC5wPnAZ4BLgU8Hi/8I3OGcuxJ4Fbg2lJ2qsBYRGZJIcraISIiyC5PTYp0q1/GeBnzIOdcKLDezg4ATzOwJYI5z7plgvXuBNWZ2nXOuPqE91gdPLy5OaDMiIhlo+HO2iEiIsgqCPtZNI7AriHPu8SBBx2wGNgInAm/Erbc5uLs44Z3u3etv1b9aRGRQIsnZIiJh8g3WuDYX6mZTorDuxWHAj4HJwM5uy2qBxK/oosJaRCQsyc/ZIiIhshwDwi+sU6UrSAczmwlUOedWmpkDGrutkg/0uLh7VVUVixd3NoosWbKEJUuW9L0jFdYiEqGlS5eydOnS2MPKKGNJxLDlbBGREFl2UFi3DqywHmjOTqnC2syygcuBLwazNgOz45YbUAxs6f7csWPHsmzZsoHvTIW1iEQovpA0sx0RhzMkw5qzRURC1NFiPcDCeqA5O9W6gnwOuME5F1S9PADMiVs+HWgGEs/GscI6Ly/hTYmIZKjhy9kiIiGKFdaEO9pe6hTWZvZl4Hmg0Mz2N7MPAdXA1uCMc4DTgJudc80J71At1iIiQzbsOVtEJESD7QoyUCnRFcTMrga+0W32K865O8zsQuCrZrYeKAGuCmWnzUGeV2EtIjIokeRsEZEQjeiTF51z1wHX9bFsA/Dh0HfaGowUlZMSL4GISNqIJGeLiIQoWS3WKdMVZNi1BwOCZ2dHG4eIiIiIDKvBnrw4UCqsszL3JRARERHJRCqsw6bCWkRERCQzBeWfa1dhHQ4V1iIiIiIZybKC4fbCratVWKuwFhEREckwQV1Ne7ibzdyqUoW1iIiISGaKdQVx6goSDhXWIiIiIhnJLGiyVot1SFRYi4iIiGSmWPmnPtYhUWEtIiIikpFiJy9qVJCwqLAWERERyUw6eTFkxcUwaRKUlUUdiYiIiIgMo5zyHPIm5pFVHG4pnBPq1tLJJz7hJxERERHJKNOumsa0q6aFvt3MbbEWEREREQmRCmsRERERkRCosBYRERERCYEKaxERERGREGR0Yb106dKoQwjFSDgOHUNq0DFIqhup76+OK73ouNLLcB6XCusRYCQch44hNegYJNWN1PdXx5VedFzpRYW1iIiIiEiaMedCvkh6RMysCtgwyKdVAjuSEM5wGwnHoWNIDTqGaExzzo2NOojhNMScDen5/g6Ejiu96LjSS9jH1WfOHjGFtYiIiIhIlNQVREREREQkBCqsRURERERCoMJaRERERCQEGVlYm1mhmd1sZt80sx+ZWVnUMQ2Emb3LzN4ws51m9gMzy4lb9mszc8F0d5Rx7ouZXRcX68pgXtq8J2b2TFz8semMYFmPY0slZnZKEP/0uHl9vvap+L70cQx9fjeC5Wnz/ZDepeJnMSy9fabT3b6+k+nKzI4zszVmttvMboo6nrCZWZ6ZrTSz46OOJSzD/Xd5RHzQh+BW4H7n3O/N7DTgJ8B7I46pX2Y2FTgXOB84EPgRsAm4IVj2FnBKsPqaSIIcADMrxp+dG4s1NipAWrwnZjYPeAr4HNASzL4deLifY0sJZjYeKAEO77aov9c+pd6X3o6hv+9G3PK0+H5Iv1LqsxiWfr6XaWtf38l0ZWblwNuBI4Gjgb+a2X3OuQejjSxUnwemRx1EWKL4u5xxo4KY2STgDaDCOdcQ/Be9BzjQObc+0uD6YWbHAE8551qDx9fjY36nmX0f+D/gr7HlqcrMPh3c/bFzrimYlzbviZlNdc5tjHs8DbjeOXdhb8eWaswsC2gDZjjn1vf32gPNfS2L8n3p5Rj6/G4Ej9Pm+yG9S6ccMRTdP9MRh5OwfX0n05WZFQJNLiiczGw58Dnn3CPRRhYOMzsamAt8FbjEOfevaCNKXBR/lzOxK8jxwA7nXANA8MV/AzguyqD2xTn3eLeiYDOw0cyKgLOAPwSPT4okwAEwMwPeA9wIbDezWGvT8aTJexJfVAfOw7da9HVsKcU5195t1vH0/dr3tywy3Y+hr+8GQDp9P6Rfx5OCn8Ww9PK9TGv9fSfTmXOuMa6oLgZWA/+KNKiQBMdzvnPujqhjCUtUf5czsbCeDOzsNq8WmBRBLIk4DP8fWINzbib+uP4G/MPM5kcbWu+cdyT+Z5mbgd+Y2emk93tyBvC3fo4t1fX32qfr+3IY8GOAdPp+SL/S9bMoXsd3ciQIWuT/AeQCRRGHE5YrgW9HHUSYovq7nImFtQMau83Lp7O/bMozs5lAlXOuoxO+c26rc+4y4I/AJyILbgCcc7udc1cD3wI+TZq+J2ZWif9ZsCY2r5djS3X9vfZp97709t2A9Pp+SK/S7rMoXl/fyTT3GvBz4CTgu9GGkjgzOxVY5px7K+pYkmG4/y5nYmG9GajoNq8U2BJBLINmZtnA5cAX+1jlf4CpwxdRQmKxput7cjZwTx/L0uV96O+1T6v3ZQDfDUif90W6SqvPongD/E6mneAf9TvwJ7GPhO5InwPuMLMdZrYDmALcbWZfiDiusA1L/s/EwvoRYFJwEgJmlot/oR+NNKqB+xxwg3Nubx/L24HnhjGeRMRiTdf3pL/COl3eh/5e+3R7X/b13YD0eV+kq3T7LIo3kO9kOluG/6cv3V0EHBw3bQEuZQR13wkMS/7PuMLaObcNuA//Ew7AycBdzrmU/3KY2ZeB54FCM9vfzD5kZgeZ2dnB8jzgA8D1UcbZFzObZGYXmVlWcFLBfwFXpeN7EpzoURD76ayvY4s0yF4EsQEY9P99SNX3pfsxBPN6+27MNLO56fL9kL6l6mcxLL19ptNdX9/JqONKhJkVmdmiuFlnAD+IKp6wOOeqnHNvxib8CDVV8d0c01FUf5czbrg9ADOrwP9xXQ9MAL7knKuLNKh9MLOrgW90m/0Kvr/ob/FnyD8NfMc5l5I/j5rZgfgTPmqAJ4CbnHOvBMvS6j0xs/OB/ZxzNwWP+zy2VGFmpcD78OMBXwv80DlX1d9rn2rvS2/HgP+pucd3wzk3LxgFJC2+H9K/VPsshqWv72W0USWmr79Xzrl5UcQTFjM7HH8S9Cv4fPJv59wfo40qfGa2nhEw3F5Uf5czsrAWEREREQlbxnUFERERERFJBhXWIiIiIiIhUGEtIiIiIhICFdYiIiIiIiFQYS0iIiIiEgIV1iIiIiIiIVBhLSIiIiISAhXWIpJSzGysmd1jZhvM7ItRxyMiIv1T3u6kC8SISEoxs1OBR4Ay/BXOZjvndkUblYiI9EV5u5MKaxFJWWb2GHCic6416lhERGTfMj1vqyuIiKQkM5sC/G+mJmcRkXSjvK3CWpLAzL5gZrdEHUdMGPGY2clm9pSZHR9SWCOKmd0Q8vYMeDfw7TC3KyJdKV8nzsymm9ltZvaTkLaXbWaXmdmGAaz7TjN7ZqjHGmbuVt72VFinGTObZ2Z/NrPbzWylmTkze1/UcXXzd+C3UQcRJ4x4NgELQ4glJZnZGWa2w8y2mtlRZjYl+Jw5M/uhmY0K1jvczF40sxvNrDCYNx7YFnJIFwE/dc61BtsXSTvK10OSjvn6reA2P6TtZQHVwNQBrLsMmNfXQjM70Mx+H3z2rjazCXHLws7dytuosE5HfwW+7py7FFgE/DK2wMwWmNlxUQQVv2/n3IvOuSeiiKM3g4mnr9fQObcG2BFmXGG+X4luyzn3d+BWoNU595RzbhNwGdAMPOOcqw7W+zc+kX/eOdcYPP0s4L6EDiCOmX0F+BrwoJm9AswJa9siw0z5epDSMV875xqALWHtyznXAiwf4Lrb8UV4X8tfApYGD29wzsUX0qHlbuXtTiqs04iZjQZmA/UAzrl24KuAM7NyfNK2COKKbN9hGsBxtA/jvqLY1u+B/czsUADn3E7gIeDMuH2NArYEn72Yuc65VxLcdwfn3Decc7Occ4udc3Odc4+FtW2R4aJ8nVwpmK/DHgliMPHva9224LZ7v+fQcrfydicV1mkkGLpmDXC/mR0ZzNsA/Bvfr2kGsMTMPghgZtea2VeCn4E6+lGZ2WIz+56ZfdrMGsxsjZl91szONbNnzew8M/uLmdWY2eVxz+t1e932fZmZ/Ty+r1nwc+jNZvYtM3vEzA4P5ve7v7jnm5ldbmYvmdkpZrbKzLab2SUD2Mes+Hj2sc8er2EvKoKfdmvN7LJur+nXgtfmz0EiHuhr1mNfZna+mf3bzC4xs7Vm9u1+ttfbe99rPP0JWjZeAi6Imz0WeKcF3T6As4F74uIsBBriHs81s2+b2Q1m9qj5LiUXmtkyMzvHzH5rZtVmdoWZHRoc404zO2Nf8YmkE+XrzMnXvenrGPvbj5nNMLM7zLf+fncg++llv383s9+Y2dv3sZ5yd7I45zSl0QQcCGzA/4d6BzAubtl64Pjg/kHAuuD+WPx/02Pw/3FvBBYHy24Bngzu5wG7gB/ix6J8L7C5v+113ze+j9lPgJ8H8wvwxdqo4PH5QBVQ2t/+ejnuCcE+vxysex3QBEzbxz66x9PvPuNfw15iWA/cCBQFz9sYzC8Dfhe33v/hk+KAXrM+9jUa//PeTcDhwMn9ba/be99rPAP8fF0NvB7cfxvwLWAv8K64z4vFrX82cExwvwh4FigKHj8WvA5lwB7gWqAQuDB4/G78P/f/RfAZ1KRpJE0oX2dEvg6Wf22Ar2N/7/WzwLy457gBfs5i7+dRwDW9LD8+2E9O3Dzl7iRNarFOM863Ks7HJ5r/BF4ys95O0lgLvN/McoFjgnkl+C/yFCDWP3Y1/guDc64ZqAX+5JyrwX/RJuxje93j20vXkyHeCTS6oI8u8Gf8F/Jd+9hf9+3GtnlnsO51+P6/p+xjH13iGcw++3Cv8/3pngUmxh3jmKBF6dP41xQG+Jr1cby78Ansbufcv51zDw5ie33FMxC/B/Y33x3k/wG3Af8E3m1mJUCtCzJv4O3Ak8H9s/F/vGKtIOcB3whe52rgYef7ZS8Dypxzf3D+5/HnGdx7IJIWlK8zI1/3os9j7Gc/JwPlzrlYPMsGuc/DgIudc18f4PrK3UmSE3UAMnBmlg1Mcc6tBz5lZj/Fn0G9FDgifl3n3F4zmwxcCfwitgnn3Ftmtgw4Cf8f9Qx8MdXx1Lj7rQTdhfraXh+hxm9jFpAbF1e7ma0H9utvf/vinGsys3XAKKByEPsY8j67aaXz+zMVeMM5d1P3lQbympnZ+/EFbMwc538ydvGxDuI96DOefXHOrTWz5/HdQcqdc+vM7Pf4VqQLgL/FxR37bMT6703DtzDFtrUzftNx99voqh11S5MRRvm6y/FlQr6O1+fr2M97swD/T8RQnQjMMrMS51xdH+u0B8eg3J1EekHSSy5x/V+dcyvxP7XN7b6imZ0IfMw5d53zIzzE+zhwmPlhn1YB39nXjvexvf5sxLeAxg9DZPi+h4nKA15N8j4GYitd+yFjZkcM4jW7Bzg4bur17PJBbK/XeAZ+OPwO+CCwLi6+bOAK/M+mMUfg+4vGbAbeYWbFcfvtt5+fyAimfN1VJuXrPo+xn/3UBM/JY2huwLf2/6CP5Xtc50nnyt1JpMI6/XzSzCbFPZ4M3B/cb8afrDEXP7RTuZnlm9nJwfJxZlaJ7+v3U/zPOE8S/LQYyKLbf+lmZvvYXvd9W9w2/grUAe8JtjUq2F+s5bOv/fVlYrDO5OC59w9gH/Hx7Guf8cfRXW+tF4ZvhSoG7jOzU82fJLOv96BjX8BE59wrcVNL3C6y4+73t734uMSjBo8AAAKMSURBVHuNx8zyzGxp0M2jP3/A/wT9JwDnXC3wD+D/XNfRQP4D300kJvZ6/9b8WNifAGJjmXZ/D7q/z2k9QoFIH5SvyZh8PdDXsdf94HNpHv6fL/C/TmBmkwaQuw3fevx+4AIze3fw3P8ys6vxrfUPxa2v3J1Myei4rSk5E/6ECAfsBn4N/Bz/U1JFsPwa4E3gHPxPO2vxLQTvxLd0/ByfpO7Bn0jRHGyvAf9FOx1owZ8gMwX4UrD8Q/1tr9u+rwAex/dbOzBYthjf0nkjvq/hQcH8PvfXx/G74PlX4LsmzI9b1tc+psfHs699xr+G3fZ9Dv4kvh8Gr0XseR8Mlh8HvBi8NzcGr/NAX7NzejnW84P9/RaYHMzr7z3tsq0+4tkPfyGDOwfwWfufbo8vBE7oNu+GXp53LPAysBO4Oph3Bv6zdjMwiWDIMeAj+D+8twXHembU3zFNmsKaUL7OpHw9FXgC3+p+8D6Osb/3+p34Xwqfxp9I/irwseDYe83d+F9FWvAt1aOBh/FF/aXA5fjxvG/Ed0OJPUe5O4mTBS+eZAgzGwd8xDl3bfA4C9+Kcrlz7upIg9sHM3PADOf7LMoQmdnnnXNDGspJRIaP8rXEU+5OD+oKknk+ge8WMDpu3onEnZSWiuJ+etLPTkNkZtlmdjFwV9SxiMiAKF+LcneaUWGdeX6M74u12sw2AQ8Abzrnnoo2rL4FJ3NcETy8xPywbzJ4JcBfnXPr9rmmiKQC5WsB5e60oq4gIiIiIiIhUIu1iIiIiEgIVFiLiIiIiIRAhbWIiIiISAhUWIuIiIiIhECFtYiIiIhICFRYi4iIiIiEQIW1iIiIiEgIVFiLiIiIiIRAhbWIiIiISAj+P/FXQJQoZluiAAAAAElFTkSuQmCC\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-as-201-earth.png',bbox_inches='tight')\n", "plt.savefig('../plots/apollo-as-201-earth.pdf', dpi=300,bbox_inches='tight')\n", "plt.savefig('../plots/apollo-as-201-earth.eps', dpi=300,bbox_inches='tight')\n", "\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }