{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Turbojet Example\n", "\n", "This example demonstrates how to use flight-mech to compute the characteristics of a turbojet." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Imports" ] }, { "cell_type": "code", "execution_count": 56, "metadata": {}, "outputs": [], "source": [ "# Python imports\n", "import os\n", "import sys\n", "sys.path.append(\"..\")\n", "\n", "# Additional imports\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "# Flight-Mech imports\n", "from flight_mech.turbine import TurbojetSingleBody" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Turbojet definition\n", "\n", "You first need to define the turbojet with its essential quantities." ] }, { "cell_type": "code", "execution_count": 57, "metadata": {}, "outputs": [], "source": [ "# Define the turbojet\n", "turbojet = TurbojetSingleBody()\n", "\n", "# Set the design inputs\n", "turbojet.compressor_efficiency = 0.86\n", "turbojet.turbine_efficiency = 0.9\n", "turbojet.OPR_design = 10\n", "turbojet.T4_max = 1700 # K" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Additionally, you need to set operating conditions." ] }, { "cell_type": "code", "execution_count": 58, "metadata": {}, "outputs": [], "source": [ "# Set the operating conditions\n", "turbojet.M0 = 0\n", "turbojet.ambient_pressure = 101325 # Pa\n", "turbojet.ambient_temperature = 285 # K" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Design mode\n", "\n", "The turbojet includes two modes, one for which the surface A4* will be considered as a variable and tuned to obtain a desired force. This mode is called the design mode. \n", "\n", "We can now tune its value to obtain the desired thrust, here 7500 N." ] }, { "cell_type": "code", "execution_count": 59, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "A4* at 7500N [m2]: 0.00813621917724611\n" ] } ], "source": [ "# Tune A4* to obtain 7500N of thrust\n", "turbojet.tune_A4_star_for_desired_thrust(7500)\n", "print(\"A4* at 7500N [m2]:\", turbojet.A4_star)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Operation mode\n", "\n", "Once the turbojet has been designed, you can switch to the operation mode. This mode allows you to compute the thrust and other quantities of the turbojet in a user-defined configuration.\n", "\n", "First let us define the operation conditions and change the mode." ] }, { "cell_type": "code", "execution_count": 60, "metadata": {}, "outputs": [], "source": [ "# Define the flight characteristics\n", "cruise_altitude = 11800 # m\n", "cruise_mach = 0.78\n", "max_mach = 0.82\n", "max_altitude = 12500 # m\n", "\n", "# Switch to operation mode\n", "turbojet.mode = \"operation\"\n", "\n", "# Set the number of points for the plots\n", "nb_points = 100" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### T4 influence" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can now for instance see the influence of a variation of T4 on the thrust. This roughly corresponds to changing the position of the throttle lever." ] }, { "cell_type": "code", "execution_count": 61, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "np.float64(3.565578125000005)" ] }, "execution_count": 61, "metadata": {}, "output_type": "execute_result" } ], "source": [ "turbojet.T4_instruction=800\n", "turbojet.tune_current_OPR()\n", "turbojet.current_OPR" ] }, { "cell_type": "code", "execution_count": 62, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAHHCAYAAABeLEexAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAUlFJREFUeJzt3Qd4k1XDxvGbblpoyyqrLFllL5XhehUEFSc4QEBkqCgOcCEOVHwFlU8RF6gg4MCBIiqIigKKguy9997QFuhu813nYPK2iEih7ZM0/991hfQZfXqSlOTumUVcLpdLAAAAfizA6QIAAAA4jUAEAAD8HoEIAAD4PQIRAADwewQiAADg9whEAADA7xGIAACA3yMQAQAAv0cgAgAAfo9ABHixWbNmqUiRIvryyy+dLorPee655+xzl9eqVq2qO++8U77iP//5j72d6bn169fP9zIB3ohABBQw8yF9JjcThrzR999/b8NGYTZnzhz7GOPj41XY7N692z62pUuX5nn4/LfbPwWzK6+80h6///7786xMQG4F5fo7AJyTjz76KMf2hx9+qOnTp/9tf506dbRmzRp5YyB6++23C3UoMoHo+eeftzVB0dHROY6tW7dOAQG+87fkTz/99LdAZB6bqelq3LhxnvyMDh06qEaNGp7tY8eO6d5779VNN91kj7mVLVv2b987adIkzZ07N0/KAZwLAhFQwLp27Zpj+88//7SB6OT9xrkGoqSkJIWHh5/TNZBTaGiofElISEi+/4yGDRvam9vBgwdtIDL7TvV77ZaSkqJHHnlEAwYM0KBBg/K9nMDp+M6fOYAfy8rK0osvvqjY2FiFhYWpdevW2rhx4yn7fyxatEiXXnqpDUJPPvmkPWaaI05Vo3Nyf5j09HRbe1CzZk37c0qVKqWLL77YBjbDnGtqh9zXdN/+zbRp03TJJZcoIiJCxYsXV/v27bVq1SrP8f/7v/+z19m2bdvfvnfgwIH2Q/3IkSOefRMnTlSzZs1UtGhRlS5d2n7o7tq167Rl2Lp1q/0Z48aN+9ux7M+PuX/sscfs19WqVfM8RvP9p3rOjM2bN+uWW25RyZIl7fPeokULTZ069ZT9wb744ot/fS1Ptnz5cvu93377rWefeZ3NvqZNm+Y49+qrr1bz5s1P2YfIlOGCCy6wX/fo0cPz2E5+TlavXq3LL7/cPpaKFSvqlVdeUX4w1zW/248++mi+XB/IDQIR4ANeeuklff311/aDwwQEU6vUpUuXv5136NAh+4FomkJef/11+6GWGyYMmEBkvu+tt97SU089pcqVK2vx4sX2+D333GP7eximic99Ox1z3ASgYsWK6eWXX9YzzzxjP3BN0HKHjFtvvdUTFk5m9rVt21YlSpSw2+bD25wfGBiooUOH6q677rLNLuZ6edHnxzTxdO7c2X49fPhwz2MsU6bMKc/ft2+fWrVqpR9//FH33XefDTum5uP666+3r9nZvpbZmaBrmu5+++03z77Zs2fbprtly5YpMTHR7jPhwjT3mUB8KqYZdvDgwfbru+++2/PYsp9vgudVV12lRo0a6dVXX1VcXJytwTGhNi9t377dPhfmd8IEW8BxLgCO6tu3r+uf/ivOnDnTHqtTp44rNTXVs3/EiBF2/4oVKzz7LrvsMrtv1KhRf7uO2f/ss8/+bX+VKlVc3bt392w3atTI1b59+7Mu78mOHj3qio6Odt1111059u/du9cVFRWVY3/Lli1dzZo1y3He/Pnz7c/68MMP7XZaWporJibGVb9+fVdycrLnvClTptjzBg0a5NlnHm/2cm7ZssVujx079l+fn2HDhtl95nv+7Tnr16+fPXf27Nk5Hne1atVcVatWdWVmZub6tTwV87pceOGFnu0OHTrYW2BgoGvatGl23+LFi+21vvnmmxy/F+bmtmDBgn98Hty/Q+7n2zBlLVeunKtjx46uM3XgwIF//J1zu/nmm12tWrXybJvzze8W4BRqiAAfYJo3svcFMc1P7qaak/u3mHPPlqmFME1ZGzZsUF4wTW2m1sbUuJh+Je6bqd0xzTozZ870nHvbbbfZZqBNmzZ59n3++ef2Md1www12e+HChdq/f7+tiTHNTW6mBsrUZJzcTFVQncwvvPBCW0PlZmrDTA2MqQEztWFn81qezJxnauqOHz9ut3///Xddc801tjbQ1BYZ5t7UtGUvS26Zsmfv92PKah7fv5UvN8zr/tVXX9laTMBbEIgAH2CarbJzNx9l71djmP4e59KJ1jSnmABTq1YtNWjQwPalMf1XzpY7WF1xxRW2ySn7zYx+MuHGzfTBMU1AJgQZptLA9BUyTYCRkZF2n7uPUe3atf/2s0wgOlUfpPxmfuapymOap9zHz+a1PFUgysjIsCOyzEg389yZfaa5K3sgqlu3ru3LdLZM36aT+4WZMv5b+c6UeQwPPvigunXr5unPBHgDRpkBPsDUqJzKiZaG/8ltX4zMzMwc2+bD1dTQfPPNNzawjB492vajGTVqlHr37p3rcps+LYbpp1KuXLm/HQ8K+t9bUIUKFewHvOkzZDqDm741pp+J6WOSF/6p8/fJz4G3vJYnO//8822tmOlHZEJVTEyMDa7mOXvnnXeUmppqA5EZ6u5E+c6UmWbCBLp3333X04fM7ejRo3afeWyMjkRBIxABfsD8hX9yh+O0tDTt2bPnb+ea2gXTrGNuZj4ZE5JMZ2t3IMrN7M/Vq1e39+YDrk2bNv96vmk2M81h5gPT1BSZD8XrrrvOc7xKlSr23hw3tU7ZmX3u46firok5+Xk4Va1Sbh6j+ZnmZ59s7dq1Ocp8rtxNVyb0mEDkbmoz9yYMffLJJ7aD9z91qHbLj9m7c8OEXDOa8aKLLjplWDI30+n8xhtvdKR88F80mQF+wAST7COUjPfee+9vtSNmlNrJ/UnMhHvmA9fNDJ03zmREV7t27Wxz15AhQ+yH4MkOHDiQY7tjx462huLTTz+1zWXXXnut5+e5a0lMuDI1VtnLZEZAmTmbTF+if2LKYYbon/w8mNqVk+XmMZp+PPPnz88xuaDp52OeXzNE3zRh5RUTfubNm2f74LgDkXlMpnnOXZPm3v9PcvPY8kOnTp1s4Dn55n4uzdfZpw0ACgo1RIAfMLU7ffr0sYHDDJs3Q7XNMHHzYZqd+fA2c9aYOX5MTZHpxGzWUcu+pII5Zph+ICbwmABjPuT+KYSMHDnS9hcx8+WY80z/IVNLYDpAm1oCM7zfzYQdM+T/tddes80npsYou+DgYPvBb2qvLrvsMttZ29SKjBgxwoaP/v37/+vzYIZ6m3sTrkw4Wr9+/d/Ocz9GM+2AKbP5uaamKns4c3viiSdsgDN9ncxzYp638ePHa8uWLbbjcF7Oam3CjhnWv2PHjhzBx9QKmSYo8xyYPkD/Fo5N53kTKs2cUOYxmQBi5lwqCKavl7mdiikDNUNwCoEI8ANmrh7zAT1mzBj98MMP9sPUjAAzkwJmZz7QzeR/pv+QqYExzT3//e9/PRMVuufpeeCBB/TZZ5/p448/tn1L/ikQGbfffrvtH2SCyLBhw+x1TedvU4ZTjYgzIejnn3+2H9amxuBkZlJE05RmrmfmxzEf6KbfjAlKJy+zcTIzG7KplTIhz/RVMiHG1C6ZIJad6ez7wgsv2NBgni/TF8o8f6cKRGY5CjP3jynLm2++aecgMjM0f/fdd6etsTobZr4jE0DN4zfzBLmZ59IEon+rHTJMuDOBzcyBZEKy6eQ8duzYAgtEgLcqYsbeO10IAAAAJ9GHCAAA+D0CEQAA8HsEIgAA4PcIRAAAwO8RiAAAgN8jEAEAAL/HPERnwMxBsnv3bjsvitPT3gMAgDNjZhYyk7yaudD+bZJUAtEZMGGoUqVKThcDAACcBTO7+7/N4k4gOgOmZsj9hJqlCAAAgPdLTEy0FRruz/HTIRCdAXczmQlDBCIAAHzLmXR3oVM1AADwewQiAADg9whEAADA7xGIAACA3yMQAQAAv0cgAgAAfo9ABAAA/B6BCAAA+D0CEQAA8HsEIgAA4PcIRAAAwO8RiAAAgN8jEAEAAEct2nZYh46lOloGAhEAAHDMgq2H1XX0fHV6708ddDAUEYgAAIAjlmw/oh5jFyg5PVPlo4uqeFiQMwUhEAEAACes3JWgOz6Yr2OpGWp5Xim927WZQoMC5RQCEQAAKFBr9yaq65h5OpqSofOrlNDo7ueraIhzYcggEAEAgAKzcf9RdXl/nuKT0tWoUrTG9rhAEaHONZW5EYgAAECB2HzgmDq/P0+HjqepXoVIfdjjQhUPC5Y3IBABAIB8t/1Qkm5/f54OHE1VXLni+qhXc0WFe0cYMghEAAAgX+08kqTO7/+pvYkpqhFTTB/3bq6SESHyJgQiAACQb/YkJNuaoV3xyTqvdIQm9G6u0sVC5W0IRAAAIF/sT0yxYWj74SRVKRWuCXe1UExkmLwRgQgAAOQ501fINJNtOXhcsSWK2jBULso7w5BBIAIAAHnKrEt2+/t/atOB46oQFaZP72qhitFF5c0IRAAAIM8cPp6mLqPnacP+YyoXGWZrhiqVDJe3IxABAIA8EZ+Upq6j52nt3qOKKR6qCXc1V9XSEfIFBCIAAHDOEpLT1W3MfK3ek6jSxUJszdB5ZYrJVxCIAADAOYehO8bM04pdCXZ+IROGzHxDvoRABAAAzlpiSrq6fzBfy3YmqER4sG0mq1W2uHwNgQgAAJyVo3+FoaU74hUdHqxPerdQXLlI+SICEQAAyLVjqRm6c+wCLdker6iiwfq4V3PVreCbYcggEAEAgFyHoR5j52vRtiOKDAuyYah+xSj5MgIRAAA4Y8f/CkMLth5R8bAgu2p9g1jfDkMGgQgAAOQiDC3whCFTM9SoUrQKAwIRAAD4V0lpGeoxboHmbz2s4qEnaoYKSxgyCEQAAOBfw1BPE4a2/BWGejdX40IUhowgpwsAAAC8vGZo7ALN+ysMfdjrwkIXhgwCEQAA+McwZIbWz88WhppULqHCiCYzAABwyg7U/hKGDAIRAAD4+2iycf4ThgyazAAAwN8mXbRD6/0kDBkEIgAA8L/lOD6Yr4Xb/CsMGQQiAAAgs1Cr6TNkluMobJMuen0foqpVq6pIkSJ/u/Xt29ceT0lJsV+XKlVKxYoVU8eOHbVv374c19i+fbvat2+v8PBwxcTE6LHHHlNGRkaOc2bNmqWmTZsqNDRUNWrU0Lhx4wr0cQIA4M0SU9J1xwf/W5vsk97+FYYcD0QLFizQnj17PLfp06fb/bfccou979+/v7777jtNnDhRv/76q3bv3q0OHTp4vj8zM9OGobS0NM2ZM0fjx4+3YWfQoEGec7Zs2WLPufzyy7V06VL169dPvXv31o8//ujAIwYAwAvD0Jj5nlXrJ9zVQg1j/SsMGUVcLpdLXsKElSlTpmjDhg1KTExUmTJlNGHCBN188832+Nq1a1WnTh3NnTtXLVq00LRp03TttdfaoFS2bFl7zqhRozRgwAAdOHBAISEh9uupU6dq5cqVnp/TqVMnxcfH64cffjijcpmyREVFKSEhQZGRkfn06AEAKFgJySYMzdOynQmKDg8uFKvWn+3nt9cMuze1PB9//LF69uxpm80WLVqk9PR0tWnTxnNOXFycKleubAORYe4bNGjgCUNGu3bt7BOwatUqzznZr+E+x30NAAD8UXxSmrqOPhGGSoQHa0LvFoUqDPlsp+rJkyfbWps777zTbu/du9fW8ERH56y2M+HHHHOfkz0MuY+7j53uHBOakpOTVbRo0b+VJTU11d7czLkAABQWR46nqcvoeVq9J1GlIkL0yV3NFVfOv1tAvKaGaMyYMbr66qtVoUIFp4uioUOH2io2961SpUpOFwkAgDxx6FiqOr//pw1DpYuF6rO7W/h9GPKaQLRt2zb9/PPPtrOzW7ly5Wwzmqk1ys6MMjPH3OecPOrMvf1v55i2xFPVDhkDBw607Y3u244dO/LokQIA4JwDR0+EobV7jyqm+IkwVLNscaeL5RW8IhCNHTvWDpk3o8HcmjVrpuDgYP3yyy+efevWrbPD7Fu2bGm3zf2KFSu0f/9+zzlmpJoJO3Xr1vWck/0a7nPc1zgVMzzfXCP7DQAAX7YvMUWd3pur9fuOqWzkiTBUI6aY08XyGo4HoqysLBuIunfvrqCg/3VpMk1VvXr10sMPP6yZM2faTtY9evSwQcaMMDPatm1rg0+3bt20bNkyO5T+6aeftnMXmVBj9OnTR5s3b9bjjz9uR6m98847+uKLL+yQfgAA/MHu+GTd9u5cbTpwXBWiwvTFPS11XhnCkFd1qjZNZabWx4wuO9nw4cMVEBBgJ2Q0nZzN6DATaNwCAwPtMP17773XBqWIiAgbrAYPHuw5p1q1anbYvQlAI0aMUGxsrEaPHm2vBQBAYbfzSJJtJttxOFmxJYrq07taqFLJcKeL5XW8ah4ib8U8RAAAX7T90IkwtCs+WVVKhdtJFytGn7r/rL9/fjteQwQAAPLe5gPHdPv787Q3MUXVSkfYmqFyUWFOF8trEYgAAChkNuw7qttHz7OjymrGFLNrk8VEEoZOh0AEAEAhsmZPop2B+tDxNMWVK66Peze38w3h9AhEAAAUEit3JajrmHmKT0pX/YqR+qhnc5WICHG6WD6BQAQAQCGwePsRdf9gvo6mZKhxpWiN73mhXb0eZ4ZABACAj5u3+ZB6jlug42mZuqBqCX1w5wUqHkYYyg0CEQAAPmz2hgO668OFSknP0kU1Sun9O85XeAgf77nFMwYAgI/6Zc0+3fvJYqVlZOny2mU0smszhQUHOl0sn0QgAgDAB01bsUcPfrZE6Zkuta1bVm/e3kShQYShs0UgAgDAx0xeskuPTFymzCyXrm1YXsNva6zgQMeXJ/VpBCIAAHzIp/O368mvV8gsvHVzs1i93LGhAgOKOF0sn0cgAgDAR3zw+xYNnrLaft2tRRU9f309BRCG8gSBCAAAH/D2zI0a9uM6+/Xdl56ngVfHqUgRwlBeIRABAODFXC6X/u+ndXp75ia7/VDrmurXpiZhKI8RiAAA8OIwZJrIxv6x1W4/cXWc+lxW3eliFUoEIgAAvJAZQfbU1yv02YIddvuFG+qpW8uqTher0CIQAQDgZdIzs/ToxGX6ZulumT7Tr9zcyI4oQ/4hEAEA4EVSMzJ1/4Qlmr56n4ICimhEpyZq37C808Uq9AhEAAB4iaS0DN3z0SLN3nBQIUEBGtW1qa6IK+t0sfwCgQgAAC+QkJxuV6xftO2IwkMCNbr7+WpVvbTTxfIbBCIAABx26Fiq7vhgvlbtTlRkWJDG9bxQTSuXcLpYfoVABACAg/YmpKjrmHnauP+YSkWE6KNezVW3QqTTxfI7BCIAAByy7dBxdRk9TzuPJKt8VJg+7t1c1csUc7pYfolABACAA9btPapuY+Zp/9FUVS0VbsNQbIlwp4vltwhEAAAUsKU74nXn2PmKT0pXXLni+rDXhYopHuZ0sfwagQgAgAI0d9Mh9R6/QMfTMtW4UrTG9bhA0eEhThfL7xGIAAAoID+v3qf7JixWWkaWWlUvpffuOF/FQvko9ga8CgAAFICvl+zUoxOX2zXKrqxbVm92bqKw4ECni4W/EIgAAMhnH87dqkHfrLJfd2haUa90bKigwACni4VsCEQAAOQTl8ult2du1P/9tN5u39mqqgZdW1cBZsVWeBUCEQAA+RSGXpy6RqN/32K3H2xdU/3b1FSRIoQhb0QgAgAgj2VkZmngpBWauGin3X7m2rrqdXE1p4uF0yAQAQCQh1LSM/Xgp0v00+p9Cgwoopc7NtTNzWKdLhb+BYEIAIA8ciw1Q3d/uFBzNh1SSFCA3urcRG3rlXO6WDgDBCIAAPLA4eNp6jF2vpbtTFBESKDe736+WlUv7XSxcIYIRAAAnKPd8cl2XbJNB46rRHiwxve8UA1jo50uFnKBQAQAwDnYuP+Y7hgzT7sTUlQhKkwf9mquGjGsWO9rCEQAAJylZX8t0nokKV3Vy0Too17NVSG6qNPFwlkgEAEAcBb+2HjQdqA2i7Q2io3S2B4XqmQEi7T6KgIRAAC5NHX5HvX/fKnSMrN0UY1Sercbi7T6Ol49AABy4eM/t+mZb1bK5ZKuaVBOw29rrNAgFmn1dQQiAADOcCmON2ds1GvTT6xL1qV5ZQ2+ob6dfBG+j0AEAMC/yMpy6fnvVmn83G12m3XJCh8CEQAAp5GakalHvlimKcv3yOSfZ6+tqzsvYl2ywoZABADAaZbi6PPRIv2+8aCCA4vo1Vsb6/pGFZwuFvIBgQgAgFM4eCxVPcYu0IpdCQoPCdS73ZrpkpplnC4W8gmBCACAk+w4nKQ7PpivLQeP27mFxvW4gKU4CjkCEQAA2azanaA7xy7QgaOpqhhdVB/1ulDnlWEpjsKOQAQAwF/mbDKzTy+yfYfiyhW3i7SWjQxzulgoAAQiAABOmn26ebWSer/7+YoMC3a6WCggBCIAgN8bP2ernvtulZ19+ur6J2afDgtm9ml/QiACAPj17NP/99M6vT1zk93u1qKKnru+HrNP+yECEQDAL6VnZmngpBX6ctFOu/3IlbV0/xU1mH3aTxGIAAB+JyktQ/d9sliz1h2wtUFDbqqv2y6o7HSx4CACEQDArxw6lqqe4xZo2c4EhQUH6J0uTXVFXFmniwWHEYgAAH5j26Hj6v7BfG09lKQS4cEac+cFalq5hNPFghcgEAEA/MLynfF2KY5Dx9MUW6KonWOoOhMu4i8EIgBAoTdz3X71/WSxktIyVbd8pF2KI4YJF5ENgQgAUKhNXLhDT0xaocwsly6pWdr2GSrOhIs4SYActmvXLnXt2lWlSpVS0aJF1aBBAy1cuDDHHBGDBg1S+fLl7fE2bdpow4YNOa5x+PBhdenSRZGRkYqOjlavXr107NixHOcsX75cl1xyicLCwlSpUiW98sorBfYYAQAFz3x+vPHLBj325XIbhm5qUlFjul9AGIL3BaIjR47ooosuUnBwsKZNm6bVq1fr1VdfVYkS/+vgZoLLG2+8oVGjRmnevHmKiIhQu3btlJKS4jnHhKFVq1Zp+vTpmjJlin777TfdfffdnuOJiYlq27atqlSpokWLFmnYsGF67rnn9N577xX4YwYA5L+MzCw9+fUKvTZ9vd3uc1l1vXZrI4UEOV4PAC9VxGUitEOeeOIJ/fHHH5o9e/Ypj5uiVahQQY888ogeffRRuy8hIUFly5bVuHHj1KlTJ61Zs0Z169bVggULdP7559tzfvjhB11zzTXauXOn/f6RI0fqqaee0t69exUSEuL52ZMnT9batWv/tZwmUEVFRdmfbWqhAADe63hqhu6fsFgz1x2QmXD6+evrqVvLqk4XCw7Izee3o1H522+/tSHmlltuUUxMjJo0aaL333/fc3zLli02xJhmMjfzwJo3b665c+fabXNvmsncYcgw5wcEBNgaJfc5l156qScMGaaWad26dbaW6mSpqan2Scx+AwB4vwNHU9X5/T9tGDJzDI3q2owwhDPiaCDavHmzrb2pWbOmfvzxR91777168MEHNX78eHvchCHD1AhlZ7bdx8y9CVPZBQUFqWTJkjnOOdU1sv+M7IYOHWqDl/tm+hwBALzbpgPH1GHkH1q+M8HOMTThrhZqW6+c08WCj3A0EGVlZalp06YaMmSIrR0y/X7uuusu21/ISQMHDrTVa+7bjh07HC0PAOD0Fm49rI4j52jH4WRVLhmur+5txYSL8J1AZEaOmf4/2dWpU0fbt2+3X5crdyLZ79u3L8c5Ztt9zNzv378/x/GMjAw78iz7Oae6RvafkV1oaKhta8x+AwB4p2kr9uj20fMUn5SuRpWiNem+VjqPCRfhS4HIjDAz/XiyW79+vR0NZlSrVs0Gll9++cVz3PTnMX2DWrZsabfNfXx8vB095jZjxgxb+2T6GrnPMSPP0tPTPeeYEWm1a9fOMaINAOBbxvy+RfdNWKy0jCy1qVNWn93VQqWLhTpdLPggRwNR//799eeff9oms40bN2rChAl2KHzfvn3t8SJFiqhfv37673//aztgr1ixQnfccYcdOXbjjTd6apSuuuoq29Q2f/58O2rt/vvvtyPQzHnG7bffbjtUm/mJzPD8zz//XCNGjNDDDz/s5MMHAJwlM6/Qc9+u0gtTVsuMlb6jZRW9262ZioYEOl00+CqXw7777jtX/fr1XaGhoa64uDjXe++9l+N4VlaW65lnnnGVLVvWntO6dWvXunXrcpxz6NAhV+fOnV3FihVzRUZGunr06OE6evRojnOWLVvmuvjii+01Klas6HrppZfOuIwJCQlmagJ7DwBwVlJqhuuu8QtcVQZMsbdRszbazwrgXD6/HZ2HyFcwDxEAeIeDx1LVe/xCLd0Rr5DAAL16ayNd1+hEawBwLp/frGUGAPAJmw8c051jF2j74SRFFQ3W+3ecrwurlXS6WCgkCEQAAK83f8th3f3RQjuSrFLJohrX40JVZyQZ8hCBCADg1b5dtluPfrFMaZlZdlj96DvOV5nijCRD3iIQAQC8kuni+s6sTRr244npWdrVK6vXb2vCSDLkCwIRAMDrpGdm6ZnJK/XZghMrBfS6uJqevKaOAs1qrUA+IBABALxKYkq6+n6yWLM3HLSr1T97XT11b8UCrchfBCIAgNfYeSRJPcct0Pp9x1Q0OFBvdm6iNnVzLs4N5AcCEQDAKyzbEa9e4xfauYZiiofqgzsvUP2KUU4XC36CQAQAcNyPq/bqoc+WKCU9S3HlitswVCG6qNPFgh8hEAEAHB1JNnr2Fg2ZtsauSXZZrTJ66/YmKh4W7HTR4GcIRAAAx0aSPfvtKk2Yt91ud21RWc9dV09BgY6uOw4/RSACADg6kqxIEenp9nXV86KqKmI2AAcQiAAABWrH4ST1Gv+/kWRvdG6iKxlJBocRiAAABWbx9iO6+0MzkizNjiQb0/0CNYhlJBmcRyACABSI75bt1iMTlyktI0t1y0dqdPfzGUkGr0EgAgDk+0iyN2ds1GvT19vtNnViNKJTE0WE8hEE78FvIwAg36RmZGrgVys0ackuu9374moayJpk8EIEIgBAvjAzTt/z0SIt2nbEBqDBN9RTl+ZVnC4WcEoEIgBAnlu/76hdk2znkWQVDwvSyC7NdHHN0k4XC/hHBCIAQJ6atW6/7p+wRMdSM1SlVLgdSVYjppjTxQJOi0AEAMizztPj52zV4CmrleWSLqxWUu92baYSESFOFw34VwQiAECeLMPx3Ler9Mlfy3Dc3CxWQ25qoJAgluGAbyAQAQDOSUJSuu6bsEh/bDxkl+F44qo43X3peSzDAZ9CIAIAnLUtB4+r17gF2nzwuMJDAu38QizDgUIbiEqWLJmri5q/ChYvXqwqVRheCQCF1e8bDuq+TxYpMSVDFaOL2pmn65SPdLpYQP4Fovj4eL3++uuKioo6o0519913nzIzM8+uRAAAr/fR3K167rvVysxyqUnlaL3X7XyVKR7qdLGA/G8y69Spk2JiYs7o3AceeODsSwQA8OrO089/t0of/3mi83SHJhU1pEMDhQUHOl00IP8DUVZWVq4uevTo0bMtDwDAS8UnpanvhMWeztOPt4tTn8voPI3CgU7VAIB/tXH/UfUev1BbDyXReRr+HYh+++23Mzrv0ksvPZfyAAC8zMx1+/XghCU6mkrnaRReZxyI/vOf//zjMXd1qbnPyMjIm5IBABxlBsmMnr1FQ6etOTHzdNWSGtm1qUoVo/M0/DgQHTly5JT7k5KSNGLECL3xxhs677zz8rJsAACHpGZk6qmvV+rLRTvtdqcLKmnwDfWZeRqF1hkHopOH3JuO1h988IGef/55BQQE6O2331b37t3zo4wAgAK0/2iK+ny0SIu3xyugiPTMtXV1Z6uqdJ5GoXZWnaonTZqkJ598UgcOHNDAgQPtMPvQUKpQAcDXrdiZoLs/Wqg9CSmKDAvSW7c31aW1yjhdLMC7AtGvv/6qAQMGaMWKFXrooYfs12cyWSMAwPt9s3SXHv9yuVIzslS9TIRGd79A1UpHOF0swLsC0TXXXKOff/5ZPXv21OTJk1WuXLn8LRkAoECY2ab/76d1Gjlrk92+vHYZjejcRJFhwU4XDSgwRVxmGMEZMP2EgoKCFBERcdp25MOHD6uwSUxMtDVhCQkJioxkqCmAwiMxJV39PluqGWv32+17Lj1Pj18Vp0DTeQjwo8/vM64hGjt2bF6UDQDgJTYfOKbeHy7U5gPHFRoUoJc7NtSNTSo6XSzAEWcciBhBBgCFx6x1+/XAp0t0NCVD5aPC9G63ZmoYG+10sQDHsHQHAPgR00vi3d826+Uf1sp0mGhWpYSdbDGmeJjTRQMcdUYzbJUsWVIHDx4844tWrlxZ27ZtO5dyAQDyWHJaph76bKlemnYiDN12fiVNuKs5YQg40xqi+Ph4TZs27YyH2B86dEiZmZnnWjYAQB7ZeSRJ93y0SKt2JyoooIieva6uuraowmSLwF/oQwQAhdyfmw/pvk8W6/DxNJWKCNE7XZqq+XmlnC4W4HuByCzTAQDwvf5C4+ds1QtT19i5hupXjNS73c63K9YDyIlO1QBQCKWkn1ic9avFJxZnvbFxBQ3t0FBFQwKdLhrglQhEAFDI7I5PVp+PF2n5zgQ7weLAq+PU6+Jq9BcCToNABACFyLzNh9R3wmIdPJamEuHBdnHWi2qUdrpYgNcjEAFAIeov9N+pa5SR5VKd8pF6r1szVSoZ7nTRAJ9AIAKAQtBf6MmvV2jS4l12+4bGFfQS/YWA/A1EgYGB2rNnj2JiYv4295DZx/xDAFBwdsUn656PFmrlrkT6CwEFGYhMteyppKamKiQk5FzKAgDIhTkbD+r+T5fY+YVMf6G3b2+qVvQXAvI3EL3xxhv23vzVMXr0aBUrVsxzzNQK/fbbb4qLizu7UgAAcvWH6fuzN9slOLJcUr0KZn6hZootQX8hIN8D0fDhwz3/EUeNGmWbztxMzVDVqlXtfgBA/jmemqHHv1quqcv32O2OTWP14k31FRZMfyGgQALRli1b7P3ll1+uSZMmqUSJEuf0gwEAubPl4HH1+WiR1u07ynpkgNN9iGbOnJlj2zSXrVixQlWqVCEkAUA++Xn1PvX/fKmOpmYopnioRnZtqmZVSjpdLKDQCMjtN/Tr109jxozxhKFLL71UTZs2VaVKlTRr1qz8KCMA+C2zBtmrP61T7w8X2jB0QdUSmvLAxYQhwOlANHHiRDVq1Mh+/d1332nr1q1au3at+vfvr6eeeiqvywcAfuvI8TT1GLdAb87YaLd7XFRVE+5qoZjIMKeLBhQ6uQ5EZr6hcuXK2a+///573XLLLapVq5Z69uxpm84AAOdu5a4EXffW7/pt/QGFBQdoRKfGeva6egoOzPXbNoAzkOv/WWXLltXq1attc9kPP/ygK6+80u5PSkrKMfIMAHB2Pl+wXR1GztHOI8mqUipcX993kW5oXNHpYgGFWq4DUY8ePXTrrbeqfv36dmRDmzZt7P558+bleh6i5557zl4j+y37NVJSUtS3b1+VKlXKznvUsWNH7du3L8c1tm/frvbt2ys8PNzOlP3YY48pIyMjxzmmb5Pp5xQaGqoaNWpo3LhxuX3YAFAgS3A88dVyDfhqhdIystSmToy+vf9iuy4ZAC8bZWZCjAlDO3bssM1lJmQYpnboiSeeyHUB6tWrp59//vl/BQr6X5FMv6SpU6fafktRUVG6//771aFDB/3xxx/2uKmlMmHINOHNmTPHLilyxx13KDg4WEOGDPFMF2DO6dOnjz755BP98ssv6t27t8qXL6927drlurwAkB92HknSvR8v1opdCTKj6B9tW1v3XlZdAQEMqQcKQhHXP63FUQBMuJo8ebKWLl36t2MJCQkqU6aMJkyYoJtvvtnuM52369Spo7lz56pFixaaNm2arr32Wu3evds25RlmcsgBAwbowIEDdsJI87UJVStXrvRcu1OnToqPj7dNfmciMTHRBjJTpshI/lIDkLdmrduvfp8vVXxSul2CY0SnJrq0VhmniwX4vNx8fue6hmjw4MGnPT5o0KBcXW/Dhg2qUKGCwsLC1LJlSw0dOlSVK1fWokWLlJ6e7mmSM0xzmjnmDkTmvkGDBp4wZJhan3vvvVerVq1SkyZN7DnZr+E+x0wf8E/Mumzmlv0JBYC8lpXl0ohfNuiNGRtk/jRtGBuld7o0ZQkOwAG5DkRff/11jm0TWkyzlGnqql69eq4CUfPmzW1/ntq1a9vmrueff16XXHKJrc3Zu3evreGJjo7O8T0m/JhjhrnPHobcx93HTneOCTnJyckqWrTo38plQpkpCwDk55B6Uyv06/oDdrtri8p65tq6Cg1icArgE4FoyZIlf9tnwsWdd96pm266KVfXuvrqqz1fN2zY0AYkM+P1F198ccqgUlAGDhyohx9+OMfjMxNPAkBeWLojXn0/Waxd8cl2SP2QmxqoQ9NYp4sF+LU8mdDCtMuZGpVnnnnmnK5jaoPMnEYbN260HaXT0tJsX5/szCgz9zxI5v7kUWfu7X87x5T5n0KX6Shujme/AcC5Ml02P5q7VbeMmmPDUNW/htQThgDn5dkMX6bDkrmdi2PHjmnTpk12BFizZs3saDEzKsxt3bp1dpi96WtkmHszGeT+/fs950yfPt0GmLp163rOyX4N9znuawBAQa1Sb5rInvlmldIzXbqqXjl9+wBD6gGfbTJ74403/vYXj+n/89FHH+VoAjsTjz76qK677jrbTGZGij377LN2+H7nzp1tr/BevXrZpquSJUvakPPAAw/YIGM6VBtt27a1wadbt2565ZVXbH+hp59+2s5d5J4OwAy3f+utt/T444/b2bRnzJhhm+TMyDMAKAgb9x+1Q+o37D+mwIAiGnh1nHpdXI1V6gFfDkTDhw/PsR0QEGCHx3fv3t32vcmNnTt32vBjlgMx17j44ov1559/2q/dP8tc30zIaEZ9mdFh77zzjuf7TXiaMmWKHVVmglJERIQtR/aRcNWqVbPhx8xpNGLECMXGxmr06NHMQQSgQHyzdJcGTlqhpLRMu0r9212a6oKqLMwKeBtH5yHyFcxDBCC3UjMy9cKU1fr4z+12u1X1UnZ+oTLFT9ReA/Cuz+9c9SEyQ+zN8PrskxwCAHLacThJt4ya6wlDD1xRQx/1ak4YAgpLk5np5GwmRjRLZgAA/m766n165IulSkzJUHR4sIbf1liX145xulgA8nqU2VNPPaUnn3xShw8fzu23AkChlZ6ZpSHfr9FdHy60YahxpWhNffASwhBQWDtVmxFbZp4gs9yGGR1mOjJnt3jx4rwsHwB4vd3xyXrg0yVatO2I3e55UTU9cXWcQoLybGYTAN4WiG688cb8KQkA+OjCrP0/X6ojSekqHhakYTc31FX1yztdLAC5xCizM8AoMwAny8jM0mvT1+udWZvsdv2KkXr79qaqUipnrTmAQrravZtZVsPMEJ2VlZVjv+l0DQCF2d6EFD346RLN33rYszDr0+3rKiyYhVkBX5XrQLR+/Xo7g/ScOXNy7DcVTWbWVUagASjMzOr0pons8PE0FQsN0tAODXRdowpOFwtAQQeiHj162LmIzAzRZs0xpp4H4C9NZMN/Xq+3Z55oIqtbPtLOOl2tNE1kgF8GoqVLl2rRokWKi4vLnxIBgJfZk5Cshz5d6mki69K8sp65liYywK8DkVlM9eDBg/lTGgDwMjPX7dfDf40io4kM8PNAZHppu7388st25fghQ4aoQYMGdvbq7BiFBaCwTLT46k/rNerXE01k9SqcGEVWlSYywH8DUXR0dI6+QqYDdevWrXOcQ6dqAIXFziNJdhTZ4u3xdvuOllX05DV1aCID/D0QzZw5M/9LAgBe4MdVe/XYxGV2+Q0z0eIrHRvq6gZMtAgUdmcUiC677DINHjxYjz76qMLDw/O/VABQwFIzMjX0+7UaN2er3TZrkb3ZuYkqleQ9D/AHZzxTdWBgoPbs2aOYGP9bqJCZqoHCbfOBY3YtslW7T/SXvOfS8/Rou9oKDmQtMsCX5ctM1azwAaAwmrR4p56evFJJaZkqGRGiV29ppMvj/O8PP8Df5WrYPZMwAigsjqdmaNA3q/TV4p12u8V5JTWiUxOVjQxzumgAvD0Q1apV619D0eHDJyYuAwBvtXJXgh78bIk2HziugCJSvza11PfyGgo0GwD8Uq4C0fPPP2/b4gDAF5mm/7F/bNVL09YqLTNL5SLDNKJTYzU/r5TTRQPgS4GoU6dOftmpGoDvO3QsVY99uVwz1u63223qlNWwmxuqRESI00UD4EuBiP5DAHzVnE0H1e+zpdp/NFUhQQF66po6drJF3tcAuDHKDEChXn7j9Z/X651Zm2TewqqXidCbnZuqbgWmzwBwloEoKyvrTE8FAMftOJxk5xZauuPE8hudLqikQdfVVXhIrte0BuAHeGcAUOh8u2y3npq0QkdTTyy/8VKHhmrfkOU3APwzAhGAQuNYaoaezTa30PlVSuj1To0VW4LlNwCcHoEIQKFgmsYe+myJth1KsnML3X9FTT14RQ0FsfwGgDNAIALg0zKzXBr16yYNn75eGVkuVYwuamuFLqha0umiAfAhBCIAPmtPQrL6f75Uf24+MUO+6Sc05KYGiioa7HTRAPgYAhEAnzRtxR49MWmFEpLTFR4SqOeuq6dbzo9lbiEAZ4VABMDnFmV9/rtV+mLhiY7TDWOj7KKs1UpHOF00AD6MQATAZyzbEa9+ny/VloPHZSqC7r2suvpfWUvBdJwGcI4IRAB8ouP0yFkb9frPG2zH6fJRYXrt1sZqWZ1FWQHkDQIRAK+fcfrhL5ZqwdYj/+s4fWMDRYXTcRpA3iEQAfBKZv3Eb5bu1jOTV9oZpyNCAjX4hvrq0LQiHacB5DkCEQCvk5CUrqe/Wanvlu22200rR+v125qocilmnAaQPwhEALzKnE0H9cgXy7QnIUWBAUX0wBU1dP/lzDgNIH8RiAB4hdSMTL3603q9P3uzXC6paqlwDb+tsZpULuF00QD4AQIRAMet3Zuofp8t1dq9R+125wsr6en2dRURylsUgILBuw0Ax2RluTTm9y0a9uM6pWVmqVREiIZ2aKC29co5XTQAfoZABMARu+KT9egXyzR38yG73TouRi91bKgyxUOdLhoAP0QgAlDgw+knL92lQZNX2eH0RYMD9cy1dW0zGcPpATiFQASgwBw5nqanJq/Q9yv22u0mlaM1/NbGqso6ZAAcRiACUCBmrtuvx79crgNHUxUUUET92tRUn8uqM5wegFcgEAHI99Xph3y/Rp/M2263a8QUs7VCDWKjnC4aAHgQiADkm0XbDuvhL5Zp26Eku93joqoacFWcwoIDnS4aAORAIAKQL5MsmpXp3/11k7JcUoWoMP3fLY3UqkZpp4sGAKdEIAKQp1bvTrSr07snWezYNFbPXl9XkWGsTg/AexGIAOSJjMwsvfvbZr3+83qlZ7pUMiJEQ25qoKvqM8kiAO9HIAJwzjYdOGYXZF26I95ut61bVkM6NFDpYkyyCMA3EIgAnNPSG2PnbNUrP6xVakaWiocF6fnr6+mmJhWZZBGATyEQATgr2w8l6dEvl2n+lsN2+5KapfXKzQ1VPqqo00UDgFwjEAHI9dIbZk4hM7dQUlqmwkMC9VT7Orr9wsrUCgHwWQQiAGdsd3yyBny1XLM3HLTbF1Yrqf+7uZEqlwp3umgAcE4IRADOqFZo4sKdemHKarsga2hQgB6/Kk49WlVVQAC1QgB8H4EIwGntTUjRE5OWa9a6A54FWc0ki9XLFHO6aACQZwhEAP6xVmjS4l16/rtVSkzJUEhQgB5tW0u9Lj5PgdQKAShkCEQA/mZfYoqenLRCv6zdb7cbxUbZWqGaZYs7XTQAyBcB8hIvvfSSHaHSr18/z76UlBT17dtXpUqVUrFixdSxY0ft27cvx/dt375d7du3V3h4uGJiYvTYY48pIyMjxzmzZs1S06ZNFRoaqho1amjcuHEF9rgAX6sV+mrRTl352q82DIUEBuixdrX11b2tCEMACjWvCEQLFizQu+++q4YNG+bY379/f3333XeaOHGifv31V+3evVsdOnTwHM/MzLRhKC0tTXPmzNH48eNt2Bk0aJDnnC1btthzLr/8ci1dutQGrt69e+vHH38s0McI+EJfod7jF+qRictsE1nD2ChNefBi9b28hoICveKtAgDyTRGX+ZPQQceOHbO1N++8847++9//qnHjxnr99deVkJCgMmXKaMKECbr55pvtuWvXrlWdOnU0d+5ctWjRQtOmTdO1115rg1LZsmXtOaNGjdKAAQN04MABhYSE2K+nTp2qlStXen5mp06dFB8frx9++OGMypiYmKioqChbpsjIyHx6JgBnmLeALxft1GAzgsz0FQoM0ENtauqeS88jCAHwabn5/Hb83c40iZkanDZt2uTYv2jRIqWnp+fYHxcXp8qVK9tAZJj7Bg0aeMKQ0a5dO/sErFq1ynPOydc257ivAfj7vEI9xi3QY18ut2HI9BWiVgiAP3K0U/Vnn32mxYsX2yazk+3du9fW8ERHR+fYb8KPOeY+J3sYch93HzvdOSY0JScnq2jRvy8zkJqaam9u5lygsNUKfb5gh16cusbOK2RGkD18ZS31vrgaQQiAX3IsEO3YsUMPPfSQpk+frrCwMHmToUOH6vnnn3e6GEC+2HE4yc4r9MfGQ3a7aeVovXJzI9WIYV4hAP7LsT8FTZPY/v37bf+hoKAgezMdp9944w37tanFMZ2lTV+f7Mwos3Llytmvzf3Jo87c2/92jmlLPFXtkDFw4EDb3ui+mfAGFIaV6cfP2ap2r/9mw1BYcICebl9HE/u0IgwB8HuO1RC1bt1aK1asyLGvR48etp+Q6QhdqVIlBQcH65dffrHD7Y1169bZYfYtW7a02+b+xRdftMHKDLk3TI2TCTt169b1nPP999/n+DnmHPc1TsUMzzc3oLDYfOCYXYNswdYjdrt5tZJ6uWNDVS0d4XTRAMC/A1Hx4sVVv379HPsiIiLsnEPu/b169dLDDz+skiVL2pDzwAMP2CBjRpgZbdu2tcGnW7dueuWVV2x/oaefftp21HYHmj59+uitt97S448/rp49e2rGjBn64osv7MgzoLDLyMzS+7O3aPjP65WWkaWIkEA9cU0ddbmwMmuQAYCvzFQ9fPhwBQQE2Boi08nZjA4zw/PdAgMDNWXKFN177702KJlA1b17dw0ePNhzTrVq1Wz4MXMajRgxQrGxsRo9erS9FlCYrd6daGuFVuxKsNuX1iqjITfVV2wJVqYHAK+bh8gXMA8RfElqRqbemrFRI2dtUkaWS1FFg/XMtXXVsWlFOxs8APiLxFx8fnt1DRGA3Fm07bAGfLVCG/cfs9tX1SunwTfWU0xx7xrJCQDehkAEFALHUzM07Md1Gj93q0ydb+lioRp8Qz1d06C800UDAJ9AIAJ83Kx1+/XU1yu1Kz7Zbt/SLFZPta+j6PAQp4sGAD6DQAT4qEPHUvXClNWavHS33Y4tUVRDOzTQJTXLOF00APA5BCLAx5hxEJOX7tLg71brSFK6zOj5nhdV08Ntayk8hP/SAHA2ePcEfGzZjacnr9Sv6w/Y7bhyxfVSx4ZqXCnnmn8AgNwhEAE+MsHiuDlb9epP65WcnqmQwAA92LqG7rmsuoJZjBUAzhmBCPByq3Yn6ImvVngmWLywWknbV6h6GdYfA4C8QiACvFRyWqZe/2W9Rs/eoswsl4qHBenJa+rotvMrsewGAOQxAhHghUwfoacnr9COwyeG0l/ToJyeu66eYiKZYBEA8gOBCPAiB4+l6r/ZhtJXiArT4Bvqq03dsk4XDQAKNQIR4CVD6Scu3Kkh09Yo/q+h9N1bVdUjbWurWCj/TQEgv/FOCzjMrDv25NcrNH/LYbtdt3yk7TTdiKH0AFBgCESAQ1LSM/XOrE0aOWuj0jNdKhocqIevrKUeF1VVEEPpAaBAEYgAB8zZeFBPTV6pLQeP2+0r4mLsYqyxJcKdLhoA+CUCEVDAnaZfnLpGXy/ZZbdjiofq2evMqvTlVKQIQ+kBwCkEIqAAZGW59PnCHXpp2lolJKfLZJ87WlTRI+1qKzIs2OniAYDfIxAB+WzNnkQ99fUKLd4e7+k0PaRDA9YfAwAvQiAC8snx1AyN+GWDxvx+Yqbp8JATnabvbEWnaQDwNgQiIB/mFJq+ep+e+3aVdiek2H1X1SunQdfVVYXook4XDwBwCgQiIA/tOJxkg9Ava/fb7dgSRe3osSvimGkaALwZgQjIA2kZWXp/9ma9OWODUtKzFBxYRHddcp4euKKmioYEOl08AMC/IBABeTCn0DPfrNSmAyfmFGp5Xim9cGM91Ygp7nTRAABniEAEnKX9iSn679Q1+nbZiYVYSxcL0dPt6+qGxhWYUwgAfAyBCMiljMwsfTh3m16bvl7HUjPsQqxdzZxCbWsrqihzCgGALyIQAbmwYOthPTN5pdbuPWq3zQKsL95YX/UrRjldNADAOSAQAWfgwNFUO8v0V4t32m1TE/TE1XG67fxKCjBVRAAAn0YgAv6leeyTedv1fz+t09GUDLuv84WV9Fi7OJWMCHG6eACAPEIgAv7BQtM89s0qu/SGUb9ipF64ob6aVC7hdNEAAHmMQASconls6LQ1mrT4xIr0kWFBeqxdbd3evIoCaR4DgEKJQAScNHps+PT1Opp6onnM9BF6/KraKlUs1OniAQDyEYEIkDR30yG75Ma6fSdGjzWoGGWX3KB5DAD8A4EIfm1PQrJenLpGU5bvsdvR4cG2eazTBZVpHgMAP0Iggl9KzcjU6Nlb9PbMjUpKy7STK3ZpbiZXrKXocEaPAYC/IRDBr7hcLs1Yu1+Dp6zWtkNJdl+zKiX0/PX1mFwRAPwYgQh+Y/OBYzYIzVp3wG7HFA/VwGvidGPjiqw9BgB+jkCEQu9oSrremrFRH/yxRemZLgUHFlGvi8/T/VfUULFQ/gsAAAhEKMSyslyatGSXXXLj4LFUu+8/tcto0LV1dV6ZYk4XDwDgRQhEKJSW7ojXs9+u0rId8Xa7WukIPXNtHV0RV9bpogEAvBCBCIXKvsQUvfzDWs8s0xEhgXqwdU31uKiaQoICnC4eAMBLEYhQKKSkZ2rM7/8bRm/c3CxWj7errZjIMKeLBwDwcgQi+Pww+h9X7dWL36/RjsPJdl+TytF67rp6alQp2uniAQB8BIEIPmvV7gS9MGW1/tx82G6XjQzVwKvr6IbGFRhGDwDIFQIRfI4ZMfbqT+v02YIdcrmk0KAA3XPpebrnsuqKYBg9AOAs8OkBn1puY+wfW+2cQsf+Wo3+2obl9cTVcYotEe508QAAPoxABJ/oJ/TDyr0aOm2tth9O8qxGP+i6urqgakmniwcAKAQIRPBqK3ed6Cc0b8thz3Ibj18Vpw5NKiqA1egBAHmEQASvnU/olR/WadKSnfQTAgDkOz5Z4FWS0zL13m+bNerXTUpOPzGf0I2NK+ixq+JUMbqo08UDABRSBCJ41bpj//fjOu1NTLH7mlUpoafb11GTyiWcLh4AoJAjEMFxczYd1ItT12jV7kS7HVuiqB051r5BeeYTAgAUCAIRHLPpwDEN/X6tfl6zz24XDw3S/VfUUPdWVRUWHOh08QAAfoRAhAJ36FiqRvyyQZ/M267MLJcCA4qoS/PKeqh1TZUqFup08QAAfohAhAJdgPWDP7bonZmbPBMrtqkTY5vHasQUd7p4AAA/RiBCgXSYnrz0RIfp3QknOkzXrxipp66pq5bVSzldPAAACETIX79vOKgh36/R6j0nOkxXiArTY1fV1g2NmFgRAOA9CETIF2v3JtoO07+uP+DpMH3f5TXU4yI6TAMAvA+BCHlqd3yyXpu+Xl8tPjHDdFBAEXVtUUUPtq6pkhEhThcPAIBTIhAhTyQkp9vZpT/4fYtSM7LsPjOP0KPtaqta6QiniwcAwGkFyEEjR45Uw4YNFRkZaW8tW7bUtGnTPMdTUlLUt29flSpVSsWKFVPHjh21b9+JOWvctm/frvbt2ys8PFwxMTF67LHHlJFxYgST26xZs9S0aVOFhoaqRo0aGjduXIE9xsIuNSNTo2dv1mXDZmrkrE02DF1YtaS+vq+V3u7SlDAEAPAJjtYQxcbG6qWXXlLNmjXlcrk0fvx43XDDDVqyZInq1aun/v37a+rUqZo4caKioqJ0//33q0OHDvrjjz/s92dmZtowVK5cOc2ZM0d79uzRHXfcoeDgYA0ZMsSes2XLFntOnz599Mknn+iXX35R7969Vb58ebVr187Jh+/zI8e+WWZGjq3Xrvhku69GTDE9cVWcWteJYYZpAIBPKeIyScSLlCxZUsOGDdPNN9+sMmXKaMKECfZrY+3atapTp47mzp2rFi1a2Nqka6+9Vrt371bZsmXtOaNGjdKAAQN04MABhYSE2K9NqFq5cqXnZ3Tq1Enx8fH64YcfzqhMiYmJNpAlJCTYmix/Zn5dfttwUC9PW+sZOVY2MlT929TSzc1iFRToaKUjAABn9fntNZ9eprbns88+0/Hjx23T2aJFi5Senq42bdp4zomLi1PlypVtIDLMfYMGDTxhyDC1PuYJWLVqleec7Ndwn+O+xqmkpqbaa2S/QVq2I163vz9P3T+Yb8OQGTn2WLvamvXo5ep0YWXCEADAZzneqXrFihU2AJn+Qqaf0Ndff626detq6dKltoYnOjo6x/km/Ozdu9d+be6zhyH3cfex051jQk5ycrKKFi36tzINHTpUzz//fJ4/Vl+15eBxO6ni1BV77HZIYIC6tayivpfXYOQYAKBQcDwQ1a5d24YfU5315Zdfqnv37vr1118dLdPAgQP18MMPe7ZNeKpUqZL8zb7EFLvm2OcLdtg1x0y3oA5NYtX/ypqKLRHudPEAACg8gcjUApmRX0azZs20YMECjRgxQrfddpvS0tJsX5/stURmlJnpRG2Y+/nz5+e4nnsUWvZzTh6ZZrZNW+KpaocMMxrN3PxVQlK6Rv66SWP/+N8Q+tZxMXaG6bhy/t2HCgBQOHldp4+srCzbh8eEIzNazIwKc1u3bp0dZm+a2Axzb5rc9u/f7zln+vTpNuyYZjf3Odmv4T7HfQ38T3Japt6ZtVGXvDLDzilkwtD5VUpoYp+WGnPnBYQhAEChFeR009TVV19tO0ofPXrUjigzcwb9+OOPtld4r169bNOVGXlmQs4DDzxgg4wZYWa0bdvWBp9u3brplVdesf2Fnn76aTt3kbuGxwy3f+utt/T444+rZ8+emjFjhr744gs78gwnpGVk6fMF2/XGjI06cDTV7qtdtrgev6q2rohjCD0AoPBzNBCZmh0zb5CZP8gEIDNJowlDV155pT0+fPhwBQQE2AkZTa2RGR32zjvveL4/MDBQU6ZM0b333muDUkREhO2DNHjwYM851apVs+HHzGlkmuLM3EejR49mDiIzsi/LpW+X7dLw6Ru0/XCS3VepZFE7hP6GxhUVyOKrAAA/4XXzEHmjwjYPkXnJf1q9T6/9tF7r9h21+0oXC9WDrWuo0wWVFRLkdS2pAADk6+e3452qUbB+33BQw35aZ+cUMiLDgnTPZdXtKvThIfw6AAD8E5+AfmLRtiN2LqG5mw/Z7aLBgep1cTXddel5iioa7HTxAABwFIGokFu1O0Gv/rReM9bu90yqeHvzynZSxTLF/XdqAQAAsiMQFVIb9x+1naXds0ubDtI3N43Vg21qqmL0qedfAgDAXxGICplth45rxM8bNHnpLmW5ZGeXvq5hBfW/spaqlY5wungAAHglAlEhsSs+WW/N2KAvFu60w+mNtnXL6uG2tZhQEQCAf0Eg8nF7E1L09syN+mzBdqVnnghC/6ldRg9fWUsNY3MujAsAAE6NQOSjzIzSI2dt0sfzttmZpo1W1Uvpkba11KxKSaeLBwCATyEQ+ZhDx1L13m+b9eHcbUpOz7T7LqhaQg9fWVstq5dyungAAPgkApGPOHI8Te/N3qzxc7YqKe1EEGpUKVqPXFlLl9QszXpjAACcAwKRl4tPStPo2Vs09o8tOv5XEGpQMUr9r6ypy2uz8CoAAHmBQOTFQWjM7yYIbdWx1Ay7r16FSLvwaus6BCEAAPISgcjLJCSla8zvm20QOvpXEKpTPlIPta6pdvXKEoQAAMgHBCIvqxEaly0IxZUrrn5tatn5hAICCEIAAOQXApEXdJYe/bvpLL3N0zR2IgjVVNu65QhCAAAUAAKRgxZtO6w7xsz3dJZ2N41RIwQAQMEiEDmoXoUoRYQGqUqpCD3UpqaurEMQAgDACQQiB4UFB+qre1sptkRROksDAOAgApHDKpUMd7oIAAD4vQCnCwAAAOA0AhEAAPB7BCIAAOD3CEQAAMDvEYgAAIDfIxABAAC/RyACAAB+j0AEAAD8HoEIAAD4PQIRAADwewQiAADg9whEAADA7xGIAACA32O1+zPgcrnsfWJiotNFAQAAZ8j9ue3+HD8dAtEZOHr0qL2vVKmS00UBAABn8TkeFRV12nOKuM4kNvm5rKws7d69W8WLF1eRIkXyPL2aoLVjxw5FRkbm6bWRe7we3oXXw/vwmngXXo/TMxHHhKEKFSooIOD0vYSoIToD5kmMjY3N159hfpH5ZfYevB7ehdfD+/CaeBdej3/2bzVDbnSqBgAAfo9ABAAA/B6ByGGhoaF69tln7T2cx+vhXXg9vA+viXfh9cg7dKoGAAB+jxoiAADg9whEAADA7xGIAACA3yMQAQAAv0cgygeZmZl65plnVK1aNRUtWlTVq1fXCy+8kGMtFfP1oEGDVL58eXtOmzZttGHDhhzXOXz4sLp06WIn24qOjlavXr107NgxBx6Rb/ntt9903XXX2ZlJzczikydPznE8r5775cuX65JLLlFYWJidKfaVV14pkMdXmF6P9PR0DRgwQA0aNFBERIQ954477rAzw2fH61Gw/0ey69Onjz3n9ddfz7Gf16RgX481a9bo+uuvt5MMmv8rF1xwgbZv3+45npKSor59+6pUqVIqVqyYOnbsqH379uW4hjm/ffv2Cg8PV0xMjB577DFlZGQUyGP0BQSifPDyyy9r5MiReuutt+wvsdk2bwRvvvmm5xyz/cYbb2jUqFGaN2+e/QVv166d/aV2M282q1at0vTp0zVlyhT7n+buu+926FH5juPHj6tRo0Z6++23T3k8L557M11+27ZtVaVKFS1atEjDhg3Tc889p/fee69AHmNheT2SkpK0ePFi+weEuZ80aZLWrVtn3/iz4/Uo2P8jbl9//bX+/PNP+0F9Ml6Tgns9Nm3apIsvvlhxcXGaNWuWDZrm/4wJmm79+/fXd999p4kTJ+rXX3+1f1R06NAhxx/qJgylpaVpzpw5Gj9+vMaNG2f/OMRfzLB75K327du7evbsmWNfhw4dXF26dLFfZ2VlucqVK+caNmyY53h8fLwrNDTU9emnn9rt1atXm+ok14IFCzznTJs2zVWkSBHXrl27Cuyx+DrzHH799dee7bx67t955x1XiRIlXKmpqZ5zBgwY4Kpdu3YBPbLC8Xqcyvz58+1527Zts9u8Hs68Jjt37nRVrFjRtXLlSleVKlVcw4cP9xzjNSnY1+O2225zde3a9R+/x7yHBQcHuyZOnOjZt2bNGnutuXPn2u3vv//eFRAQ4Nq7d6/nnJEjR7oiIyNzvEb+jBqifNCqVSv98ssvWr9+vd1etmyZfv/9d1199dV2e8uWLdq7d69tqnEz1aDNmzfX3Llz7ba5N9XQ559/vuccc75ZV83UauDs5NVzb8659NJLFRIS4jnH1DKZ2o0jR44U6GMqbBISEmyzgXkNDF4PZxa07tatm21SqVev3t+O85oU7GsxdepU1apVyz5/pqnLvF9lb1YzNXCm+Tn7+5qpTapcuXKO9zXTNF22bFnPOeZ6pibP1PSBJrN88cQTT6hTp072FzI4OFhNmjRRv379bBWzYT6Qjey/mO5t9zFzb37xswsKClLJkiU95yD38uq5N/enukb2n4HcM82Wpk9R586dPQtV8noUPNPMb57jBx988JTHeU0Kzv79+23frJdeeklXXXWVfvrpJ9100022Ocw0jbmfTxM83X9E/NP7Gq/H6bHafT744osv9Mknn2jChAn2r6ulS5faQGTa4bt37+508QCvZP7CvfXWW22nd9MHD84wtQ0jRoywfbpMTR2cryEybrjhBttPyGjcuLHtB2T6QV522WUOl7DwoIYoH5hqZnctkamiNFXP5hd56NCh9ni5cuXs/ckjAMy2+5i5N38ZZGdGA5iRHe5zkHt59dyb+1NdI/vPQO7D0LZt22wnXXftkMHrUbBmz55tn2/T3GJqfczNvC6PPPKIqlatas/hNSk4pUuXtq9B3bp1c+yvU6eOZ5SZeT5NZ+n4+PjTvq/xepwegSgfmJEzpi09u8DAQE/SN8PxzS+g6WfkZtpxTdt7y5Yt7ba5N7/c5q81txkzZthrmPZjnJ28eu7NOWZUjfkgdzMf5LVr11aJEiUK9DEVljBkpj74+eef7bDh7Hg9Cpb5A86MYjI12+6bqd02f+j9+OOP9hxek4JjmsLMEHvT9yo700fVjOAzmjVrZrtnZH9fM+ebwJT9fW3FihU5gqz7j4+Tw5bfcrpXd2HUvXt3OzpjypQpri1btrgmTZrkKl26tOvxxx/3nPPSSy+5oqOjXd98841r+fLlrhtuuMFVrVo1V3Jysuecq666ytWkSRPXvHnzXL///rurZs2ars6dOzv0qHzH0aNHXUuWLLE38yv+2muv2a/do5by4rk3ozrKli3r6tatmx2F89lnn7nCw8Nd7777riOP2Vdfj7S0NNf111/vio2NdS1dutS1Z88ezy37yBdej4L9P3Kyk0eZGbwmBfd6mM8QM4rsvffec23YsMH15ptvugIDA12zZ8/2XKNPnz6uypUru2bMmOFauHChq2XLlvbmlpGR4apfv76rbdu29v/aDz/84CpTpoxr4MCBjjxmb0QgygeJiYmuhx56yP5yhoWFuc477zzXU089leMN3gz/fuaZZ+wbhhny3bp1a9e6detyXOfQoUP2DaZYsWJ2aGSPHj3sfxyc3syZM+2bysk3E1Tz8rlftmyZ6+KLL7bXMAHYBC3k7vUwfzCc6pi5me9z4/Uo2P8jZxKIeE0K9vUYM2aMq0aNGvYzpVGjRq7JkyfnuIb5g+6+++6zUx2Y4HnTTTfZPyyy27p1q+vqq692FS1a1P6R/sgjj7jS09ML7HF6uyLmH6drqQAAAJxEHyIAAOD3CEQAAMDvEYgAAIDfIxABAAC/RyACAAB+j0AEAAD8HoEIAAD4PQIRAPxl69atdkFTczMLaObGc8895/ne119/Pd/KCCB/EIgA+Bx38Pinmwkn2R06dEixsbH22MkLYJ6KWVMt+7pQ5nonBySzCGp0dLT69etnZvzXo48+qj179tifA8D3BDldAADILRM83D7//HMNGjQox+KXxYoVy3F+r1691LBhQ+3ateuMrm8WmD15kdnspk6dqltuuUVPPPGE/dnun2luZiFnAL6HQATA55QrV87zdVRUlK35yb4vu5EjR9paIRNcpk2bds4/e8KECerRo4deffVV3X///ed8PQDegUAEoNBavXq1Bg8erHnz5mnz5s3nfL23335bDz/8sD744AN16dIlT8oIwDvQhwhAoZSamqrOnTtr2LBhqly58jlfb82aNbZGyNQ4EYaAwodABKBQGjhwoOrUqaOuXbvmyfVMZ+mmTZvagJW9DxOAwoFABKBQmjFjhiZOnKigoCB7a926td1funRpPfvss7m+XvHixe3os4iICF1++eWEIqCQoQ8RgELpq6++UnJysmd7wYIF6tmzpx0uX7169bO6ZokSJWwoatu2rf7zn/9o5syZqlChQh6WGoBTCEQACqWTQ8/BgwftvWlGM/MHnS3zvdOnT1e7du1sKJo1axahCCgEaDIDgFwyQ/1/+ukn2/x22WWXnfH8RgC8VxGXmWIVAGCX7qhWrZqWLFmS66U73KpWrWpnrzY3AL6DGiIAOEmrVq3sLTeGDBliZ6revn17vpULQP6hhggA/pKRkWFriYzQ0FBVqlTpjL/38OHD9maUKVPGNqsB8B0EIgAA4PdoMgMAAH6PQAQAAPwegQgAAPg9AhEAAPB7BCIAAOD3CEQAAMDvEYgAAIDfIxABAAC/RyACAADyd/8PlnqHIzc7yWwAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Compute thrust evolution with T4 instruction\n", "T4_array = np.linspace(800, turbojet.T4_max, nb_points)\n", "thrust_array = np.zeros(nb_points)\n", "for i, T4 in enumerate(T4_array):\n", " turbojet.T4_instruction = T4\n", " # Use try block to avoid crash if temperatures become negative\n", " try:\n", " turbojet.tune_current_OPR()\n", " thrust_array[i] = turbojet.thrust\n", " except Exception as e:\n", " print(e)\n", "\n", "# Plot thrust evolution with velocity\n", "plt.plot(T4_array, thrust_array)\n", "plt.xlabel(\"T4 [K]\")\n", "plt.ylabel(\"Thrust [N]\")\n", "plt.title(\"Thrust evolution with T4\")\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Cruise\n", "\n", "Let us now check the thrust obtained in cruise conditions." ] }, { "cell_type": "code", "execution_count": 63, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "cruise thrust [N]: 2527.981038576462\n" ] } ], "source": [ "# Switch to cruise conditions\n", "# T4 not realistic but to see what happens at full power\n", "turbojet.T4_instruction = turbojet.T4_max\n", "turbojet.altitude = cruise_altitude\n", "turbojet.M0 = cruise_mach\n", "\n", "# Compute thrust in cruise\n", "turbojet.tune_current_OPR()\n", "print(\"cruise thrust [N]:\", turbojet.thrust)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Velocity influence\n", "\n", "The velocity also has an influence on the thrust of the engine." ] }, { "cell_type": "code", "execution_count": 64, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAHHCAYAAABeLEexAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAUxJJREFUeJzt3Qd0VNXaxvE3PSSkECAJJaFKCb0KoiiKIIKK4FVEkSaoFC+giKhgF8ReQD4boF5AUbGAoEgVAVF67x2S0FJJz3zr3ThjAgEDJJn2/601zsw5ZyZ75kTOk109LBaLRQAAANyYp70LAAAAYG8EIgAA4PYIRAAAwO0RiAAAgNsjEAEAALdHIAIAAG6PQAQAANwegQgAALg9AhEAAHB7BCKghC1ZskQ8PDzk66+/tndRnM5zzz1nvruiVrVqVenTp484ixtuuMHcCnts/fr1xZlNnTrVnPe//vrL3kWBCyMQAUVA/7EuzE3DkCP66aefTNhwZStWrDCfMSEhQVzN0aNHzWdbv359sYRF/d1t3759gfs/+ugj2+83gQXOzNveBQBcweeff57v+WeffSYLFiw4b3vdunVl27Zt4oiBaOLEiS4dijQQPf/886YmKDQ0NN++HTt2iKen8/x9+Msvv5wXiPSzaXhp3Lhxkf88f39/Wbx4scTGxkpkZGS+ff/73//M/vT09CL/uUBJIhABReD+++/P93zVqlUmEJ27XV1pIDpz5owEBARc0XsgPz8/P3Emvr6+Jfrz2rRpI3/++ad8+eWX8t///te2/fDhw/Lbb7/JnXfeKd98802Jlgkoas7zJxHgYnJzc+Xll1+WypUrm7+wb7rpJtm9e3eB/T/WrFkjbdu2NUHoqaeeMvu0iaKgGp1z+8NkZWWZ2oOrrrrK/JyyZcvKtddeawKb0mO1dsj6ntbbv5k3b55cd911EhgYKEFBQdK5c2fZsmWLbf/rr79u3ufAgQPnvXb06NHmon769GnbtlmzZkmzZs2kVKlSUq5cORMmjxw5ctEy7N+/3/wM7WNyrrzfj96PHDnSPK5WrZrtM+rrC/rO1N69e+U///mPhIWFme+9VatWMnfu3AL7g3311Vf/ei7PtXHjRvPaH374wbZNz7Nua9q0ab5jO3XqJFdffXWBfYi0DC1atDCP+/bta/ts534nW7dulXbt2pnPUqlSJZkwYYIUln6mbt26yfTp0/NtnzFjhpQpU0Y6duxY4OfT77R69erm9Vqz1K9fPzl58uR5x+p57t+/v1SsWNGEUz1HjzzyiGRmZuY7LiMjQ0aMGCHly5c3v3caxI4fP17ozwFcDDVEgJ2MHz/eNNM8/vjjkpiYaC5Q9913n/zxxx/5jtMLiF4Qe/ToYUJCRETEJf0cDQPjxo2TBx98UFq2bClJSUmmr8fatWvl5ptvloceesg0uRTUxHchelzv3r3NhfDVV181tVYffPCBCVrr1q0zAePuu++WJ554woQFaxix0m0dOnQwF1OlF2+9mOuFXcsaFxcn77zzjvz+++/m/c5t4rpUejHfuXOnuYC/9dZbJnApvbAWRH/+NddcYz7Xo48+akLktGnT5Pbbbzed4fVCfDnnMi8Nuvq5li1bZt5XaW2Lvs+GDRvMeQoODjbBWZv7Bg4cWOD7aDPsCy+8IGPHjjXHaEhVWn4rDZ633HKL+R70vOhnGDVqlDRo0MD8bhVGz549zTnbs2eP1KhRw2zTgHTXXXeJj4/Pecfr75OGSj2vGoY0LH/44YfmXmtQraFbf/f091L7dmn569SpYwKSllG//7y1YUOHDjW/M88++6wJs2+//bYMGTLE1FwBV8wCoMgNHjzYcqH/vRYvXmz21a1b15KRkWHb/s4775jtmzZtsm27/vrrzbbJkyef9z66/dlnnz1ve5UqVSy9e/e2PW/UqJGlc+fOl13ecyUnJ1tCQ0MtAwYMyLc9NjbWEhISkm9769atLc2aNct33OrVq83P+uyzz8zzzMxMS3h4uKV+/fqWtLQ023Fz5swxx40dO9a2TT9v3nLu27fPPJ8yZcq/fj+vvfaa2aav+bfvbNiwYebY3377Ld/nrlatmqVq1aqWnJycSz6XBdHz0rJlS9vzbt26mZuXl5dl3rx5ZtvatWvNe33//ff5fi/0ZvXnn39e8Huw/g5Zv2+lZY2MjLR07979ouWzfjdazuzsbPOaF1980WzfunWred+lS5ean6uPtRxWZ86cOe+9ZsyYYY5btmyZbdsDDzxg8fT0zPdaq9zcXHNvff/27dvbtqnhw4eb7yohIeFfPwfwb2gyA+xE/3LO+9ev9S97/as6L21C0GMvl9ZC6F/lu3btkqKgf/nrX/P33nuvnDhxwnbz8vIyzTra+dbqnnvuMc1AWqtgpX/N62e64447zHOtrYqPj5dBgwaZphUrbYLT2oJzm6lKqpO51lpojZdV6dKlTQ2G1kxo89PlnMtz6XFaU5eammqeL1++XG699VbTMVpri5Tea21K3rJcKi173v5sWlb9fP9Wvrz0/GrtktayWTtTR0VF2T7rubTp00o7XOvviDY7Kv3MSmu/vvvuO7ntttukefPm573HuU23+v3n3aY/Oycnp8BmWeBSEYgAO4mOjs733Np8lLdfjdL+HlfSiVabUzTA1KpVyzSRaPOV9u+4XNZgdeONN5omp7w3Hf2k4cZK++BoE5C1SUMrbrSvkDbTaHOQsl7Mateufd7P0kBkj4ud/syCyqPNU9b9l3Muz6UX9OzsbFm5cqUZ6abfnW7T/mJ5A1FMTIzpy3S5tG/TueFCy/hv5Suo2UzDoDbpaXOZNuNeqL/ZqVOnTAdsbeLVcKS/H9o3SGmzotL+P9o0WNh5ki73ewYKgz5EgJ3oX9wFOdvaU/Bf2oWhfzHnpRdXraH5/vvvTWD5+OOPTT+ayZMnm35Fl0r/qrf2Izp3CLby9v7nnxXtJKsXeO0zpJ3Bte/IwYMHTb+jonChi/G534GjnMtzaa2I1oppPyK92IeHh5vgqt/ZpEmTTCdi6ygue5TvXFoDqP2Hhg0bJvv27TMB6UK0Nkn7PmkA1xovraXS3x3ty2T9HbLX5wAKQiACnJT+dXzuJIM6KufYsWPnHau1C9qso7eUlBQTkrSztTUQXcrsz9YOtXrxvtBkfXlps5k2h2kNiNYU6SgnbSKxqlKlirnX/VrrlJdus+4viLWG4NzvoaBapUv5jPoz9Wefa/v27fnKfKWsTVcaejQQWZuf9F7DkDZLaQdvPV8XUxyzd1+INpW+9NJLprbsQnMeaY3NwoULzehG7extdW6zrdYaaU3h5s2bi73cwL+hyQxwUhpMtGYhLx3Fc27tyLnDnPUv9Zo1a5oLrpUOYVaFmcVZR5bpReyVV14xQ/rPde4w6O7du5u/7LXviTaXdenSxfbzrLUkGq60xipvmXRYv87ZpH2JLkTLoSPGzv0etHblXJfyGbUfz+rVq01TlpX289HvV0fQaRNWUdHwo6PRtO+VNRDpZ9LAYa1Ju1A/ncv5bFdKQ7SO8nrjjTf+tSbn3JobHRWWlzandu3aVX788ccCZ7mm5gcliRoiwEnphenhhx82gUOHz2u/jp9//tk2pNxKL946Z43O8aM1RXrh0SHNOlzZSvcpHWKugUcvaNo/5EIhRIfY9+rVy8yXo8fpX/raFKYdoHUSv/fff992vIYdnf/mzTfflOTkZFNjlJcO2dYLv9ZeXX/99aYGwjrsXsPH8OHD//V70GHveq/hSsORDrE/l/UzPv3006bM+nO1pipvOLN68sknTYDTvk76nej3psPutZlIJyAsylmtNezoHEaHDh3KF3y0Vuj//u//zHegfYD+LRxr53kNlTonlH4mbd6y9tkpSlo79m8zmuvviJZfpx/Q0Kz94LS5Vr+/c2mw1n167rXTtAZBreXU8KydzK90ygWgsAhEgJMaMGCAucB88sknMn/+fHMx1RFgOilgXnpB18n/9KKjNTB6QdMmj7xzA+n8NDrHy8yZM+WLL74wf5lfKBAp7Tui/YM0iLz22mvmffWip2UoaESchqBff/3VXKy19uVcOoGfNqXp++n8ONZJ9zQo/dsFUZtktFZKQ572VdIQo7VLGsTy0jmOXnzxRRMa9PvSfiz6/RUUiLQjsPZ/0bK89957ZpRUw4YNTU3GxWqsLofOF6QBVD9/o0aNbNv1u9RA9G+1Q0rDnQY2nfBSQ7J21J4yZUqxBKLC0k7X+julk37q75POYaTnRX9v8tLfG60hGzNmjGki1E7Wuk3PIzOyoyR56Nj7Ev2JAAAADoY+RAAAwO0RiAAAgNsjEAEAALdHIAIAAG6PQAQAANwegQgAALg95iEqJJ2z5OjRo2YelZKcJh8AAFw+nV1IJ4XVObAuNqkqgaiQNAxFRUXZuxgAAOAy6GzwF5v1nUBUSFozZP1CdVp6AADg+HT2c63QsF7HL4RAVEjWZjINQwQiAACcy791d6FTNQAAcHsEIgAA4PYIRAAAwO0RiAAAgNsjEAEAALdHIAIAAG6PQAQAANwegQgAALg9AhEAAHB7BCIAAOD2CEQAAMDtEYgAAIDbIxABAAC7WnPglCSeybJrGQhEAADAbhLTsmTAZ2uk3RtLZHtskt3KQSACAAB28+7CXXIqNVPCAn2lRvnSdisHgQgAANjF7vgUmbZiv3k8pkuM+HjZL5YQiAAAgF28NHerZOdapH3dcLm+VnmxJwIRAAAocYu3x8uSHcfFx8tDnu4cI/ZGIAIAACUqMztXXpyz1Tzu16aaVCsXKPZGIAIAACXqs5X7Ze+JVClX2leG3FhTHAGBCAAAlJgTKRnyzq+7zOORHWtLkL+POAICEQAAKDGv/7xDkjOypX6lYLmrWZQ4CgIRAAAoERsPJ8iXfx0yj5+9rZ54eXqIo7BrIBo3bpy0aNFCgoKCJDw8XLp27So7duw477iVK1fKjTfeKIGBgRIcHCxt27aVtLQ02/5Tp07JfffdZ/aFhoZK//79JSUlJd97bNy4Ua677jrx9/eXqKgomTBhQol8RgAAIJKba5HnftgiFotI18YVpUXVMHEkdg1ES5culcGDB8uqVatkwYIFkpWVJR06dJDU1NR8YeiWW24x21evXi1//vmnDBkyRDw9/ym6hqEtW7aY95gzZ44sW7ZMBg4caNuflJRkXl+lShVZs2aNvPbaa/Lcc8/Jhx9+WOKfGQAAdzR73RFZezBBAny95MlOdcXReFgsmtUcw/Hjx01NkQYlrQVSrVq1kptvvllefPHFAl+zbds2iYmJMUGpefPmZtv8+fPl1ltvlcOHD0vFihXlgw8+kKefflpiY2PF19fXHPPkk0/Kd999J9u3by9U2TRUhYSESGJioqmJAgAAhZOcniU3vrFUjidnyKhb6sgjN9SQklLY67dD9SHSwqqwsLPVaPHx8fLHH3+YkHTNNddIRESEXH/99bJ8+fJ8NUjaTGYNQ6p9+/amBklfaz1GA5Y1DKmOHTua5rnTp0+X4CcEAMD9vL9otwlDVcsGSL9rq4ojcphAlJubK8OGDZM2bdpI/fr1zba9e/eae23eGjBggKn5adq0qdx0002ya9fZIXta66OBKS9vb28TqnSf9RgNU3lZn1uPOVdGRoZJlXlvAADg0uw5niKf/r7PPB57W4z4eXuJI3KYQKR9iTZv3iwzZ87MF5LUQw89JH379pUmTZrIW2+9JbVr15ZPP/202Dt8axWb9aYdsQEAQOFpr5wXftwqWTkWubFOuNxYJ3/lhCNxiECknaS1M/TixYulcuXKtu0VKlQw99pHKK+6devKwYMHzePIyEjTtJZXdna2GXmm+6zHxMXF5TvG+tx6zLlGjx5tmvCst0OHzg4TBAAAhbNga5ws3Xl2vTJdzd6Redo7OWoYmj17tixatEiqVauWb3/VqlVNp+hzh+Lv3LnTjBhTrVu3loSEBDN6zErfS2uXrr76atsxOvJMR7FZ6Yg0rWkqU6ZMgWXz8/Mzna/y3gAAQOGkZ+XIC3+vV/bgddUdYr0yhw1E2kz2xRdfyPTp081cRNqfR2/WOYY8PDxk5MiR8u6778rXX38tu3fvljFjxpiRYTrXkLW2SIflax8jHZb/+++/m5DVo0cPE6ZUz549TYdqfY0Oz//yyy/lnXfekREjRtjz4wMA4LI+WLJHDp9Okwoh/jLUQdYrc9hh9xp4CjJlyhTp06eP7fn48eNl4sSJphmsUaNGZlLFa6+91rZft2sI+vHHH83osu7du5sQVbp06XwTM2oA0+H55cqVk6FDh8qoUaMKXVaG3QMAUDgHT56R9m8tNavaT+zZVDo3PNsFxh4Ke/12qHmIHBmBCACAwnlw2p/y67Z4aVOzrHzR/+oLVoCUBKechwgAADi3RdvjTBjy9vSQ52+vZ9cwdCkIRAAAoMg6Uj/3w9mO1P2vrSY1w4PEWRCIAABAkfhw2V45eOqMRAT7ydCbrhJnQiACAABX7NCpMzJx8W7z+Klb60ppP29xJgQiAABwxZ77YYtkZOdK6+pl5fZGZ6e9cSYEIgAAcMUzUi/cHm9mpH6xq/N0pM6LQAQAAC5bWqZ2pN5im5HamTpS50UgAgAAl+39xbvkSEKaVAot5RQzUl8IgQgAAFyW3fEpZmSZGntbjAT4OldH6rwIRAAA4JLpQhfP/rBZsnIs0q52eekQEyHOjEAEAAAu2Y8bj8nvu0+Kr7enPOdEM1JfCIEIAABcksS0LHlxztkZqQffUFOqlA0UZ0cgAgAAl+SNX3bI8eQMqV4uUB6+obq4AgIRAAAotA2HEuTzVQfM45e61hc/by9xBQQiAABQKNk5ufLU7E1isYjc2aSSXFOznLgKAhEAACgUrRnacjRJgv29zXplroRABAAA/lVsYrq88ctO83hUpzpSPshPXAmBCAAA/CsdVZaSkS1NokPl3hbR4moIRAAA4KIWbY+TuZuOiZenh+lI7enp3HMOFYRABAAALuhMZraM+e7s4q39r60m9SqGiCsiEAEAgAt6a8FO2+Ktw9pfJa6KQAQAAAq0+UiifPr7fvNYm8qcefHWf0MgAgAA58nJtZg5h/S+c8MK0q5OuLgyAhEAADjPZyv3y8bDiRLk7y3PdokRV0cgAgAA+RxLTJPXf95hHo+6pY6EB/uLqyMQAQAAG4vFYkaVpWbmSNPoUOnZ0vXmHCoIgQgAANjM2xwrv26LEx8vDxnfvaFLzjlUEAIRAAAwEs9kydjvz8459Mj1NaRWRJC4CwIRAAAwxs3bJidSMqR6+UAZ1K6muBMCEQAAkJV7TsrMPw+Zx+O7NRR/Hy9xJwQiAADcXHpWjplzSPW8OlpaVgsTd0MgAgDAzb23aJfsO5Eq4UF+8mSnOuKOCEQAALixLUcT5f+W7jWPX7ijvgT7+4g7IhABAOCmsnNyZdQ3GyU71yK31IuUW+pHirsiEAEA4KY++m2fbD6SJMH+3vJC13rizghEAAC4ob3HU+StX3eax2O6xEh4kOsvz3ExBCIAANxMbq5Fnvxmk2Rm58p1V5WTu5pVFndHIAIAwM38b/VBWb3/lAT4eskrdzYQDw/3WJ7jYghEAAC4kaMJaTL+p23m8ciOtSUqLMDeRXIIBCIAANxoJfvR324yK9k3q1JGHmhd1d5FchgEIgAA3MTXaw7L0p3HxdfbU17t3kC83GQl+8IgEAEA4AbiktLlxTlbzeNh7a+SmuHus5J9YRCIAABwg6ayp2dvkqT0bGlYOUQGXlfd3kVyOAQiAABc3A8bjsqv2+LFx8tDXrurkXh7cfk/F98IAAAu7Hhyhjz7wxbzeOiNV0ntSJrKCkIgAgDAhZvKxn6/WRLOZElMhWB55IYa9i6SwyIQAQDgouZsPCbzNseKt6eHvPafhuJDU9kF8c0AAOCiTWVaO6QGt6sp9SqG2LtIDo1ABACACzaVPfPdJjn9d1OZBiJcHIEIAAAXHFX285Y4M6rs9f80MhMx4uL4hgAAcCHxSeky9vuzo8oevfEqiakYbO8iOQUCEQAALtRU9tTsTZKYliUNKoXIw4wqKzQCEQAALuLbtUfMBIy+Xp6mqYxRZYXHNwUAgAs4mpAmz/14tqnsv+2ZgNGpAtG4ceOkRYsWEhQUJOHh4dK1a1fZsWPHBasBO3XqJB4eHvLdd9/l23fw4EHp3LmzBAQEmPcZOXKkZGdn5ztmyZIl0rRpU/Hz85OaNWvK1KlTi/WzAQBQUvQaOeqbjZKcni1NokPlobasVeZUgWjp0qUyePBgWbVqlSxYsECysrKkQ4cOkpqaet6xb7/9tglD58rJyTFhKDMzU1asWCHTpk0zYWfs2LG2Y/bt22eOadeunaxfv16GDRsmDz74oPz888/F/hkBAChuX6w6IL/tOiH+Pp7yxn9Yq+xyeFg0VjqI48ePmxoeDUpt27a1bdcQ06VLF/nrr7+kQoUKMnv2bFObpObNm2f2HT16VCIiIsy2yZMny6hRo8z7+fr6msdz586VzZvPTlClevToIQkJCTJ//vxClS0pKUlCQkIkMTFRgoPpsQ8AcAz7T6RKp3d+k7SsHHnuthjp06aavYvkUAp7/XaoCKmFVWFhYbZtZ86ckZ49e8rEiRMlMjLyvNesXLlSGjRoYAtDqmPHjuYL2LJli+2Y9u3b53udHqPbLyQjI8O8R94bAACOJCfXIo/N2mDCUOvqZeWB1lXtXSSn5TCBKDc31zRltWnTRurXr2/bPnz4cLnmmmvkjjvuKPB1sbGx+cKQsj7XfRc7RkNOWlraBfs3aaK03qKioq74MwIAUJQ++m2vrDlwWkr7eZu1yjw9z+9agsLxFgehfYm0SWv58uW2bT/88IMsWrRI1q1bV+LlGT16tIwYMcL2XMMToQgA4Ci2xybJm7/sNI/H3hYjlcsE2LtITs0haoiGDBkic+bMkcWLF0vlypVt2zUM7dmzR0JDQ8Xb29vcVPfu3eWGG24wj7UZLS4uLt/7WZ9bm9gudIy2JZYqVarAMuloNN2f9wYAgCPIyM6RYTPXS2ZOrrSvGy7/afbPtRNOGIi0P7eGIe0kreGnWrX8HcGefPJJ2bhxo+lUbb2pt956S6ZMmWIet27dWjZt2iTx8fG21+mINQ0wMTExtmMWLlyY7731GN0OAICzeXPBTtkemyxlA31lXLeGBY7ChhM1mWkz2fTp0+X77783cxFZ+/xonx2tudGanYI6UkdHR9vCkw7T1+DTq1cvmTBhgnmPZ555xry31vKohx9+WN5//3154oknpF+/fiZ8ffXVV2bkGQAAzuSPvSflw2V7zeNx3RpI+aCz1zo4cQ3RBx98YEaWafOXDqe33r788stCv4eXl5dpbtN7rfG5//775YEHHpAXXnjBdoyGJw0/WivUqFEjeeONN+Tjjz82I80AAHAWyelZMuKrDaIT5tzdvLJ0qHd+pQFcYB4iR8Y8RAAAe3t81gb5es1hiQorJfP+29aMLoMLzkMEAAAKNn9zrAlD2l3ozbsbE4aKGIEIAAAHF5eULqO/3WgeP3x9DWlR9Z8JjFE0CEQAADiw3FyLaSo7fSZL6lUMluHta9m7SC6JQAQAgAObumK/WbjVz9tT3unRWHy9uXQXB75VAAAceDbq8fO3m8dPd64rNcOD7F0kl0UgAgDAAaVn/T0bdXautKtdXnq1qmLvIrk0AhEAAA7o9Z932GajnnBXI2ajLmYEIgAAHMzyXSfk4+X7zOMJdzVkNuoSQCACAMCBnEzJkBFfnV27876ro+WmuhH2LpJbIBABAOAgdPGIUd9slPjkDKkZXlqe6Xx2kXIUPwIRAAAO4otVB+TXbfHi6+Up7/ZoIqV8vexdJLdBIAIAwAHsiE2Wl+ZuM4+f7FRHYiqybmZJIhABAOAAQ+wfnbFOMrJz5Yba5aVvm6r2LpLbIRABAGBn437aJjvikqVcaV95jSH2dkEgAgDAjhZsjZNpKw+Yx6/9pxFD7O2EQAQAgJ0cS0yTkV9vMI8fvLaatKsdbu8iuS0CEQAAdpCTa5H/zlwvCWeypEGlEHniljr2LpJbIxABAGAH7y3aJav3nZJAXy95794mrGJvZ3z7AACUsD/2npR3F+4yj1++s4FULRdo7yK5PQIRAAAl6HRqpgz7cr3kWkTualZZujapZO8igUAEAEDJLs2hnaiPJaZL9fKB8vzt9exdJPyNQAQAQAn5ZPm+s0tzeHuafkOBft72LhL+RiACAKAErD+UIOPnbTePx3SJkXoVQ+xdJORBIAIAoJglpmXJkOlrJTvXIp0bVJD7r462d5FwDgIRAADF3G9o1Ncb5fDpNIkKKyXjujdgaQ4HRCACAKAYfbbygMzfEis+Xh4ysWdTCfb3sXeRUAACEQAAxWTj4QR5ee4283h0p7rSsHKovYuECyAQAQBQDBLPZMmg/62VzJxc6RATIX3bVLV3kXARBCIAAIqh39DjX2+w9RvSVezpN+TYCEQAABTDfEMLtsaJr5en6TcUUop+Q46OQAQAQBFac+BUnvmG6DfkLAhEAAAUkVOpmTJk+joz31CXhhXk/lZV7F0kFBKBCACAIpCTazGLtpp1ysoFyvjuDek35EQIRAAAFIH3Fu2SZTuPi7+Pp0y8r6mUZp0yp0IgAgDgCi3deVzeWbjLPH65awOpWyHY3kXCJSIQAQBwBY4kpMl/Z64Ti0Xk3pbR0r1ZZXsXCZeBQAQAwGXKyM4xky8mnMmSBpVC5NnbYuxdJFwmAhEAAJfppTnbZMOhBDPP0KT7moq/j5e9i4TLRCACAOAyfLv2sHy+6oB5/NY9jSQqLMDeRcIVIBABAHCJth5NkqdmbzKPH72xptxYJ8LeRcIVIhABAHCJi7Y+/MUaSc/KletrlZf/tq9l7yKhCBCIAAAopFwz+eI6OXjqjFm09Z0ejcXLk8kXXQGBCACAQtK5hhbvOC5+3p7ywX3NJDTA195FQhEhEAEAUAiLtsfZJl985c4GUr9SiL2LhCJEIAIA4F/sO5Eq/5253jzu1aoKky+6IAIRAAAXkZKRLQM/+0uS07OlWZUyMqYLky+6IgIRAAAXYLFY5PGvNsiu+BSJCPaTD+5rKr7eXDpdEWcVAIALmLRkj8zfEis+Xh7ywf3NJDzY395FQjEhEAEAUIDF2+Pl9V92mMcv3FFfmkaXsXeRUIwIRAAAFNiJ+uwK9j2vjjar2MO1eRfmoLCwsEt6Uw8PD1m7dq1UqVLlcssFAIBdJKdnyYDP/pKk9GxpGh3KCvZuolCBKCEhQd5++20JCQkpVAe0QYMGSU5OTlGUDwCAEp2JeviX62V3fIpEBvvL5PubiZ83K9i7g0IFItWjRw8JDw8v1LFDhw69kjIBAGAXb/26U37dFm9Gkv1fLzpRu5NC9SHKzc0tdBhSycnJUr169X89bty4cdKiRQsJCgoy79+1a1fZseNsBzZ16tQpE65q164tpUqVkujoaHn00UclMTEx3/scPHhQOnfuLAEBAeZ9Ro4cKdnZ2fmOWbJkiTRt2lT8/PykZs2aMnXq1EJ/HgCA65u78Zi8t2i3eTy+WwNpFBVq7yLBXTpVL126VAYPHiyrVq2SBQsWSFZWlnTo0EFSU1PN/qNHj5rb66+/Lps3bzYhZv78+dK/f3/be2jTnIahzMxMWbFihUybNs0cN3bsWNsx+/btM8e0a9dO1q9fL8OGDZMHH3xQfv75Z7t8bgCAY9l6NEken7XBPB5wXTXp1pSZqN2Nh0U7/RTCsmXLCvWGbdu2vezCHD9+3NTwaFC60PvMmjVL7r//fhOavL29Zd68edKlSxcTnCIiIswxkydPllGjRpn38/X1NY/nzp1rQlXeJkDtG6UBqzCSkpJMHyqtnQoODr7szwgAcCwnUjLkjvd/lyMJaXLdVeVkat+WrGDvQgp7/S50H6IbbrjhoqPKrPfnNlVdCmtT2MVGtVk/kIYhtXLlSmnQoIEtDKmOHTvKI488Ilu2bJEmTZqYY9q3b5/vffQYrSm6kIyMDHPL+4UCAFxLRnaOPPz5GhOGqpULlPfvbUoYclOFbjI7ffp0gbcjR46YPjvaN6dOnTqXXRDtp6QBpU2bNlK/fv0Cjzlx4oS8+OKLMnDgQNu22NjYfGFIWZ/rvosdoyEnLS3tgv2bNFFab1FRUZf92QAAjkcbSJ6ZvVn+OnBagvy95ePezSUkwMfexYKjB6K84UBv2hFam69atmwpM2bMkIkTJ8rGjRsvuyDal0ibtGbOnFngfg0v2g8oJiZGnnvuOSluo0ePNrVR1tuhQ4eK/WcCAErOJ8v3yaw1h0UrhCb2bCo1ype2d5FgR4VuMsvr22+/laeeesr00dHgoCPBtIbocg0ZMkTmzJlj+ilVrly5wFFrt9xyiwlhs2fPFh+ffxJ8ZGSkrF69Ot/xcXFxtn3We+u2vMdo05uOXiuIfp4r+UwAAMdeluOVn7aZx7p6fdta5e1dJDjTKDPt7NyqVSvp1auXdOvWTfbu3SuPP/74ZQcHra7UMKQhZ9GiRVKtWrUCa4Z05Jl2jv7hhx/E3z//nBCtW7eWTZs2SXx8vG2bjljTsKO1SdZjFi5cmO91eoxuBwC4l11xyfLojHWSaxG5t2WU9Lmmqr2LBGcKRLfeeqvcfPPN0rhxY9mzZ4+88sorhZq5+t+ayb744guZPn26qf3Rvj56s/brsYYhHVH2ySefmOfWY6wzYet+DT4a0jZs2GCG0j/zzDPmva1B7eGHHzbh7YknnpDt27fLpEmT5KuvvpLhw4dfUfkBAM7lZEqG9Jv2pyRnZEvLamHy/O31bQOD4N4KPeze09PTjOwKDAy86C+PTqZY6B9+gfeZMmWK9OnTx0ymqHMHFUTnFqpa9WyqP3DggBlVpsdr+Xr37i3jx4+3jURTuk8D0NatW02z3JgxY8zPKCyG3QOA848ou++jP0wn6iplA2T2oDYSFuhr72KhmBX2+l3oQKQTHhaGhhFXRCACAOell7rHvtog3647YkaUaRiqGU4naneQVNTzELlq0AEAuL5JS/aYMKRzDE26rylhCI61dAcAAMXtp03H5LWfz66T+fzt9eS6qxhRhssMRDpztE6KWFi6CKv26wEAwJ7WHTwtw79cbx7raLL7W1Wxd5HgoArVZKZrfumaYYUdVXby5EnbKDAAAOzh0KkzMuCzvyQjO1durBMuz3Sua+8iwYHRhwgA4HIS07Kk39Q/5URKpsRUCJb37m0i3l70EsEVBiJdZwwAAGeQlZMrg/63RnbFp0hEsJ980qe5BPpd1sIMcCPEZQCAyy3Y+vvukxLg6yWf9G4hFUIKXqIJyItABABwqeH1X/51yCzYqs1k9Std2YoKcB8EIgCAS/hu3RHb8Ppnb6snN9WNsHeR4EQIRAAAp7dizwkZ+fUG83jAddWkNwu24hIRiAAATr96/UOfr5GsHIvc2iBSRndieD1KIBB5eXlJfHx8gXMP6T4AAEpKfFK69JnypySnZ0uzKmXkzbsbi6d2IAKKOxBdaC3YjIwM8fVl1WAAQMlIyciWftP+lCMJaVKtXKB89EBz8ffhD3NcnkJPzPDuu++aew8PD/n444+ldOl/FsbTWamXLVsmderUucxiAABwaXMNPfLFGtl8JEnCAn1lat8W5h4o9kD01ltv2WqIJk+enK95TGuGqlatarYDAFCc9Do06puN8tuuE1LKx0s+7dNCqpQNtHex4C6BaN++fea+Xbt28u2330qZMmWKs1wAABTo9V92yLdrj4iXp4dMuq+pNI4KtXeR4I59iBYvXpwvDGlz2fr16+X06dNFXTYAAPL5fOV+mbh4j3k87s4G0q5OuL2LBHcNRMOGDZNPPvnEFobatm0rTZs2laioKFmyZElxlBEAAJm/OVbG/rDFPB5xcy25u0WUvYsEdw5Es2bNkkaNGpnHP/74o+zfv1+2b98uw4cPl6effro4yggAcHOr9p6UR2euEx3ofG/LaBl6Y017FwnuHoh0vqHIyEjz+KeffpL//Oc/UqtWLenXr59s2rSpOMoIAHBjW48myYBpf0lmdq50iImQF++oZ0Y8A3YNRBEREbJ161bTXDZ//ny5+eabzfYzZ84wMSMAoEgdOnVGek9ZLckZ2dKyapi8e28T8fZikQXYcZSZVd++feXuu++WChUqmITevn17s/2PP/5gHiIAQJE5kZIhvT75Q44nZ0idyCD5qDcTL8KBAtFzzz0n9evXl0OHDpnmMj8/P7Nda4eefPLJ4igjAMANZ6HuO+VP2X/yjFQuU0qm9WspIaV87F0suDAPy4XW4kA+SUlJEhISIomJiRIcHGzv4gCAy0rPyjFhaOXek2b26a8fbi3Vy/+zOgJQHNfvS64heuGFFy66f+zYsZf6lgAAGNk5ufLojHUmDJX285ZpfVsShlAiLjkQzZ49O9/zrKwsM4u1t7e31KhRg0AEALgs2mAx+ttN8svWOPH19jSLtTaoHGLvYsFNXHIgWrduXYHVUX369JE777yzqMoFAHCzMDRu3naZteaweHqIvHdvE2ldo6y9iwU3UiRjF7VN7vnnn5cxY8YUxdsBANzMpCV75MNle83j8d0bSsd6Z+e7A0pKkU3moJ2V9AYAwKX4bOV+ee3nHebx07fWlbubsyQHnKDJ7N133z2vmvPYsWPy+eefS6dOnYqybAAAF/ft2sMy9vuz65PpchwD2la3d5Hgpi45EL311lv5nnt6ekr58uWld+/eMnr06KIsGwDAxRdrHfn1RvO4zzVVzYKtgNMEIh1RBgDAlfht13EzvD4n1yJ3NassY7vEsD4ZnKcPkQ6x1+H1mzdvLr4SAQBc2p/7T8nAz9ZIZk6udKofKeO7NRBPHVoGOEsg8vHxkejoaLOwKwAAl2rDoQQzC3VaVo60rVVe3u7RmMVa4RAu+bfw6aeflqeeekpOnTpVPCUCALikrUeT5IFPV5t1ylpVD5P/u7+Z+HmzWCuctA/R+++/L7t375aKFStKlSpVJDAwMN/+tWvXFmX5AAAuYHd8slm5PjEtS5pGh8rHvVtIKV/CEJw4EHXt2rV4SgIAcEn7T6RKz4/+kJOpmVK/UrBM6dvSrFMGOBJWuy8kVrsHgEt36NQZuef/VsrRxHSpHREkMwe2kjKBvvYuFtxIUnGtdm+VmZkp8fHxkpubm2+7droGAOBIQprc+9EqE4aqlw+Uzx9sSRiCw7rkQLRz507p37+/rFixIt92rWjSOSQYgQYAiE1Ml54frZLDp9OkatkAmTGglYQH+du7WEDRBaK+ffuauYjmzJkjFSpUYCItAEA+8Ulnw9CBk2ckKqyUTB/QSiKCCUNwsUC0fv16WbNmjdSpU6d4SgQAcFrHkzOk58d/yN4TqVIptJRMf7CVVAwtZe9iAUU/D1FMTIycOHGieEoDAHDuMPTRKtkdnyKRwf6mmSwqLMDexQKKLhBpD23r7dVXX5UnnnhClixZIidPnsy3T28AAPcNQ7usYWhgK4kuSxiCizWZhYaG5usrpB2ob7rppnzH0KkaANzTiZTzw1C1cvkn7QVcIhAtXry4+EsCAHDqMBQR7EcYgmsHouuvv15eeOEFefzxxyUggCpQAMDZZrL7Pl4lO+POhqGZA1sThuD6naqff/55SUlJKd7SAACcZmh9jw9XEobgfsPuWeEDAJB30kUdWl8xxN/MM1SVMAR3moeISRgBwL3pchzWSRd1niFdm4yh9XC7QFSrVq1/DUWnTp260jIBABx0odaeH6+SQ6fSzAzUOs9Q5TKEIbhhINJ+RLpiLADAvew7kWpqho4lpkuVv9cmYwZquG0g6tGjh4SHhxdfaQAADmdXXLJZjkNHldUoHyj/e7CVRIawNhncdJRZcfQfGjdunLRo0UKCgoJM0Oratavs2LEj3zHp6ekyePBgKVu2rJQuXVq6d+8ucXFx+Y45ePCgdO7c2UwJoO8zcuRIyc7OzneMzqzdtGlT8fPzk5o1a8rUqVOL/PMAgKvZcjRR7vlwlQlDdSKD5MuHWhOG4N6BqDhGmS1dutSEnVWrVsmCBQskKytLOnToIKmpqbZjhg8fLj/++KPMmjXLHH/06FHp1q2bbb/OjK1hKDMzU1asWCHTpk0zYWfs2LG2Y/bt22eOadeunVmcdtiwYfLggw/Kzz//XOSfCQBcxfpDCXLvh6vkVGqmNKgUYjpQlyvtZ+9iAcXCw+JA4+mPHz9uang0+LRt21YSExOlfPnyMn36dLnrrrvMMdu3b5e6devKypUrpVWrVjJv3jzp0qWLCUoRERHmmMmTJ8uoUaPM+/n6+prHc+fOlc2bN+dr/ktISJD58+cXqmy6Tpv2n9IyBQcHF9M3AACOYdXek/LgtL8kJSNbmlUpI1P6tpBgfx97Fwu4ZIW9fl/yavfFSQurwsLCzP2aNWtMrVH79u1tx9SpU0eio6NNIFJ636BBA1sYUh07djRfwJYtW2zH5H0P6zHW9yhIRkYGC9cCcEuLt8dL709XmzDUunpZ+axfS8IQXJ7DBKLc3FzTlNWmTRupX7++2RYbG2tqeHRx2bw0/Og+6zF5w5B1v3XfxY7RkJOWlnbB/k2aKK23qKioIvy0AOCY5m48JgM++0sysnPlpjrhpmYo0O+Sxt8ATslhApH2JdImrZkzZ4ojGD16tKmxst4OHTpk7yIBQLH66s9DMnTGWsnOtchtjSrK5F7NxN/Hy97FAkqEQ8T+IUOGyJw5c2TZsmVSuXJl2/bIyEjTWVr7+uStJdJRZrrPeszq1avzvZ91FFreY84dmabPtS2xVKmC59HQ0Wh6AwB38MnyffLinK3m8b0to+Slrg3Ey5PVCeA+7FpDpP25NQzNnj1bFi1aJNWqVcu3v1mzZuLj4yMLFy60bdNh+TrMvnXr1ua53m/atEni4+Ntx+iINQ07MTExtmPyvof1GOt7AIC70n+H3/hlhy0MDbiumrxyJ2EI7seuo8wGDRpkRpB9//33Urt2bdt27bNjrbl55JFH5KeffjJD6TXkDB061GzXIfbWYfeNGzeWihUryoQJE0x/oV69eplh9a+88opt2L32S9JmuX79+pnw9eijj5qRZ9q5ujAYZQbA1eTkWuTZHzbLF6sOmucjO9aWQTfUYN1KuJTCXr/tGogu9D/dlClTpE+fPraJGR977DGZMWOGGfmlAWbSpEm25jB14MABE5x08sXAwEDp3bu3jB8/Xry9/2kR1H06p9HWrVtNs9yYMWNsP6MwCEQAXElmdq6M+Gq9zNl4TPSf4hfvqC/3t6pi72IB7hmInAmBCICrOJOZLY98sVaW7jwuPl4e8ubdjU0nasCdr98O0akaAFAydNbpvlP/lA2HEsTfx1Mm399MbqjNGpUAgQgA3MTh02fkgU9Xy97jqRIa4COf9mkhTaPL2LtYgEMgEAGAG9gRmywPfPqHxCVlSMUQf/msf0upGR5k72IBDoNABAAu7s/9p6T/1D8lKT1brgovbcJQhZCC52AD3BWBCABc2PzNx+TRmevNqDJdpPWT3s0lNMDX3sUCHA6BCABc1LQV++W5H7eIjiW+OSZC3u3RREr5shQHUBACEQC4mNxci7z683b5v6V7zfP7W0XL87fXZ/Zp4CIIRADgQrRp7ImvN8h364+a58w+DRQOgQgAXETimSx56Iu/ZNXeU+Lt6SHjuzeUu5r9s2A2gAsjEAGACzh06oyZcHF3fIqU9vOWSfc1lba1ytu7WIDTIBABgJPbeDhB+k39S06kZEhksL9M6dtC6lZgiSHgUhCIAMCJ/bo1TobOWCdpWTkmBE3p00IiQ/ztXSzA6RCIAMAJ6brcnyzfJy//tM0Mq9fmsYk9m0iQv4+9iwY4JQIRADiZ7JxcM7/QF6sOmuf3toyWF+6oJz5envYuGuC0CEQA4ESS07Nk8PR1smzncdGR9E/fWlf6X1uNYfXAFSIQAYATjSR7cNpfsiMuWUr5eMnbPRpLx3qR9i4W4BIIRADgJAu0Pvz5GjmZminhQX7ySe8W0qByiL2LBbgMAhEAOLiv1xyWp77dJJk5uVKvYrB83Ls5q9UDRYxABAAOKifXIhPyrEnWqX6kvHF3Iwnw5Z9uoKjxfxUAOGjn6eFfrpdft8Wb50NvrCnD29cSTxZoBYoFgQgAHMz+E6ky4LO/ZFd8ivh6e8prdzWUOxpXsnexAJdGIAIAB7J81wkZPH2tJKZlSUSwn/xfr+bSOCrU3sUCXB6BCAAcZObpqSv2y0tzt5m+QxqCPuzVTMKDWYYDKAkEIgCws/SsHHnmu81mNJnq3rSyvHxnffH38bJ30QC3QSACADs6mpAmD3+xRjYeThTtL/0UM08DdkEgAgA7+WPvSdNf6ERKpoQG+Mj79zaVa68qZ+9iAW6JQAQAdugv9NnKA/LinK2SnWuRuhWCTX+hqLAAexcNcFsEIgAoQWmZOfLU7E0ye90R8/z2RhXl1e4NpZQv/YUAeyIQAUAJOXAyVR76fI1sj00WL08PGd2pDv2FAAdBIAKAErBwW5wM+3K9JKdnS7nSvvJ+z6bSqnpZexcLwN8IRABQjHROobd/3SnvLdptnjerUkYm9mwqkSHMLwQ4EgIRABSTEykZ8uiMdbJiz0nzvHfrKvJ05xizHAcAx0IgAoBi8Of+UzJk+lqJS8qQAF8vGdetAeuRAQ6MQAQARTyk/qPf9sqr83eY5rKrwkvLB/c3lZrhQfYuGoCLIBABQBE5nZopj8/aIAu3x5vnXRtXlFe6NZAAX/6pBRwd/5cCQBFYc+C0DJ2+Vo4mpps+QmO7xMh9V0czpB5wEgQiALgCublnm8he+3mHmXW6WrlAeb9nE6lXMcTeRQNwCQhEAHAFo8i0iWzJjuPm+W2NKsord9aXIH8fexcNwCUiEAHAZVix+4SZaDE+OcM0kT17W4z0bEkTGeCsCEQAcAmyc3LlrV93yqQle8RiEakZXto0kdWJDLZ30QBcAQIRABTSoVNnTK2QdqBW97aMkrFd6rEwK+ACCEQAUAjfrz8iz8zeLMkZ2RLk5y3jujeQLg0r2rtYAIoIgQgALiI5PUvGfr9FZq87YluL7O17GktUWIC9iwagCBGIAOACtGls2Jfr5NCpNPH0EHn0pqtkSLua4u3FWmSAqyEQAcA5snJy5b2Fu+T9xbsl1yJSKbSUvNOjsTSvGmbvogEoJgQiAMhjz/EUGfHletlwONG2/MYLXetLMHMLAS6NQAQAfy/K+sUfB+XluVslPStXQkr5yMt31qfjNOAmCEQA3F5sYro88c1GWbbz7IzT19YsJ6//p5FEhvjbu2gASgiBCIDb0lqh79YfkWe/3yJJ6dni5+0po26pI32uqSqe2osagNsgEAFwSydTMuTp2Ztl/pZY87xR5RB54+7GZuZpAO6HQATA7fy06ZiM+W6znEzNFG9PD/nvTVfJIzfUYDg94MYIRADcqlZIJ1mcu+mYeV47IkjeuLuR1K8UYu+iAbAzAhEAtzB34zEZ8/1mOZWaKV6eHjLohhoy5Maa4ufNOmQAROxaP7xs2TK57bbbpGLFiuLh4SHfffddvv0pKSkyZMgQqVy5spQqVUpiYmJk8uTJ+Y5JT0+XwYMHS9myZaV06dLSvXt3iYuLy3fMwYMHpXPnzhIQECDh4eEycuRIyc7OLpHPCMC+4pPT5eHP18jg6WtNGNJaoe8GtZHHOtQmDAFwjECUmpoqjRo1kokTJxa4f8SIETJ//nz54osvZNu2bTJs2DATkH744QfbMcOHD5cff/xRZs2aJUuXLpWjR49Kt27dbPtzcnJMGMrMzJQVK1bItGnTZOrUqTJ27NgS+YwA7DeC7Os1h+XmN5eZjtNaK6TLbvwwtI00qEwTGYD8PCz6r4YD0Bqi2bNnS9euXW3b6tevL/fcc4+MGTPGtq1Zs2bSqVMneemllyQxMVHKly8v06dPl7vuusvs3759u9StW1dWrlwprVq1knnz5kmXLl1MUIqIiDDHaC3TqFGj5Pjx4+Lr61uo8iUlJUlISIj5mcHBwUX++QEUnSMJafLUt5tk6d/zCtWrGCwT7moo9SoShAB3k1TI67dDD6m45pprTG3QkSNHzF97ixcvlp07d0qHDh3M/jVr1khWVpa0b9/e9po6depIdHS0CURK7xs0aGALQ6pjx47mC9qyZYsdPhWA4pKTa5FPl++Tm99casKQ79/zCn0/uA1hCIDzdqp+7733ZODAgaYPkbe3t3h6espHH30kbdu2NftjY2NNDU9oaGi+12n40X3WY/KGIet+674LycjIMDcrDVAAHNfWo0ky+tuNtjXIWlQtI+O7N5Qa5ZlXCIALBKJVq1aZWqIqVaqYTtjagVo7YeetFSoO48aNk+eff75YfwaAK5eWmSPvLNwlH/2219QQBfl7y+hOdaVHiyhmmwbg/IEoLS1NnnrqKdOvSDtFq4YNG8r69evl9ddfN4EoMjLSdJZOSEjIV0uko8x0n9L71atX53tv6yg06zEFGT16tOnUnbeGKCoqqsg/J4DLt3hHvIz9frMcOpVmnt/aIFKeu62ehAezBhmAS+OwfYi0b5DetJksLy8vL8nNzbV1sPbx8ZGFCxfa9u/YscMMs2/durV5rvebNm2S+Ph42zELFiwwHat0GP+F+Pn5mWPy3gA4hrikdBn0vzXSd8qfJgxVCPGXD3s1k0n3NSMMAXC+GiKdZ2j37t225/v27TM1QGFhYaZj9PXXX2/mDNI5iLTJTIfVf/bZZ/Lmm2+a47XXeP/+/U1Njr5GQ8vQoUNNCNIRZko7YGvw6dWrl0yYMMH0G3rmmWdM05uGHgDOIzsnV75YdUBe/2WnpGRkm6H0fa+pKsNvriWBfg5b4Q3ACdh12P2SJUukXbt2523v3bu3mStIw4s2Xf3yyy9y6tQpE4q0k7XOPaTD9K0TMz722GMyY8YM0wlaR5BNmjQpX3PYgQMH5JFHHjE/LzAw0Lz/+PHjTUftwmLYPWBfaw6cNuuPbT12doBD46hQefnO+oweA1Ak12+HmYfI0RGIAPutP/bq/O3y1V+HzfOQUj4ysmNtubdltKkhAoCiuH5TxwzAIemIsel/nG0eS0zLMtvubl7ZzCtUtjTN3QCKFoEIgMP5Y+9Jee7HrbLt7+axuhWC5aWu9aRZlTB7Fw2AiyIQAXAYxxLTZNxP2+WHDUfN82B/b7MI631XR4u3l8MOigXgAghEAOwuPStHPlq2VyYt2SNpWTmiYya0j9DjHWpLWGDh1hsEgCtBIAJgNzqmY+6mY6ZWSBdkVc2qlJHnb68n9SsxegxAySEQAbCLjYcT5MU5W+XP/afN84oh/vLkrXXltoYVbNNqAEBJIRABKFFaE/Ta/O3y3fqz/YRK+XjJw9fXkIFtq0spXy97Fw+AmyIQASgRyelZ8sGSPfLJ8n2SkX12+Z1uTSrJyFtqS4WQUvYuHgA3RyACUKyycnJlxuqD8u7CXXIiJdNsu7pamDzTOUYaVKafEADHQCACUKwdpl//eYfsP3nGbKteLlCe7FRHbo6JoJ8QAIdCIAJQ5FbsOSGvztsuGw4nmuflSvvJf9tfJT1aRIkP8wkBcEAEIgBFOnLstZ93yG+7TpjnAb5eprP0gOuqsxo9AIfGv1AArtju+BR545cdMm9zrHnu4+VhJlYceuNVUj6IdccAOD4CEYDLdvDkGXln4S6Zve6w5FrEzDB9Z5NKMrx9LYkKC7B38QCg0AhEAC7Z0YQ0eW/Rbpn11yHJ1iQkYjpK61IbtSOD7F08ALhkBCIAhRabmC4fLNktM1Yfksycs3MJta1VXkbcXEsaR4Xau3gAcNkIRAAKtQq9Tqo4M08Q0rmEHu9YW1pUDbN38QDgihGIAFy0aWzy0vxBqGXVMDOE/poaZZlLCIDLIBABOM+Bk6mmRuibtYclK+dsH6GW1cJkWPurpHV1ghAA10MgAmCzKy5ZJi3ZI9+vP2JGjSkNQENvqkkQAuDSCEQAZO3B0zJ5yR75ZWucbdsNtcvLkHY1pTl9hAC4AQIR4MZrjS3bdcKMGlu195TZphVAHWMiZXC7miy8CsCtEIgAN1x9/scNR+XDZXtle2yybWZpnVBxYNsaUjO8tL2LCAAljkAEuInk9CyZsfqgTPl9vxxLTLetNaZLbPS/tppUDC1l7yICgN0QiAAXd+jUGROCvvrrkKRkZJttur5Yn2uqyn1XR0togK+9iwgAdkcgAly0f9BfB07LJ7/tk1+2xtpGjGlz2MDrqssdTSqKn7eXvYsJAA6DQAS4kPSsHNM/aNrK/bL5SJJt+3VXlTPNYm2vKi+engydB4BzEYgAF1la44tVB8waY6dSM802P29P01G637XVpFYEC64CwMUQiAAnlZtrkRV7Tsrnq/bLr9viJefvdrGKIf7Sq3VV6dEiSsoE0j8IAAqDQAQ4mYQzmfL1msPyvz8Oyr4Tqbbtuthq3zZVpX3dCPH28rRrGQHA2RCIACfqJD3jj4Myd9Mxycg+u9BqaT9v6da0ktzfqgrNYgBwBQhEgAPT/kCz1x0x8wftjk+xba8TGSS9WleRro0rSaAf/xsDwJXiX1LAwWhfoN92HTfzBi3YGmdbbb6Uj5fc3qii3Ht1tDSqHMJCqwBQhAhEgIPYezxFvll7WL5de8Q2k7RqUClE7m4RJV0bV5Qgfx+7lhEAXBWBCLCjxDNZ8uPGoyYIrTuYYNseGuBjmsPubh4lMRWD7VpGAHAHBCKghGVk58ji7cflu3VHZNH2eMnMOdtB2svTQ9peVU66N6tsRor5+zCTNACUFAIRUEL9glbvOyXfrz8iP206JknpZ9cUs3aQ7t60sllOIzzI367lBAB3RSACinGo/LpDCWYpjbkbj0l8coZtX2Swv9zRuKLc0biS1K0QRAdpALAzAhFQxCFow+FEmbfpmJkv6PDpNNu+YH9vuaV+pOkbdHX1sqaJDADgGAhEQBEsobHu0GmZtylW5m2OlSMJ/4SgAF8v6RATIV0aVpTrapVjhXkAcFAEIuAyZGbnyqq9J+XnLbHyy9Y4OZ6nOUxD0I11wuXWBhWkXe1wKeVLCAIAR0cgAgopMS1Llu48Lr9ujZPFO+IlOU/H6CA/b2n3dwi6oXZ5RogBgJMhEAEXoYun6tD4hdvizCix7L9XlFflSvvJzTER0rFehLSuUZbmMABwYgQiII/0rBz5Y98pWbw9XpbsiJf9J8/k239VeGm5qW6EtK8bLk2iy9AxGgBcBIEI4u6jwnbFp8iyncdNc5jWAllXklc+Xh7SslqY6QukkyVWLRdo1/ICAIoHgQhuJzYxXX7ffUJ+33NCVuw+KbFJ/6wbZp0jqF2d8nJD7XBpU7OclGY1eQBwefxLD5d3IiXD1Pys3HNSVuw5IXuOp+bb7+ftaeYF0mUz2tYqb5rFmCgRANwLgQguWQO0ev8p+XPfKflj30nZGZeSb79mHV1BXmt/2tQoJ82rlmFUGAC4OQIRnH6NsF3xybL2QIL8pSHowCk5dOqfiRHzrhfWqnpZaVU9zNyHBvjapbwAAMdEIIJTOZmSIRsPJ5o1wtYeOC3rDyVISsY/8wEpHfgVUzFYmlfR8BMmLauVlbBAAhAA4MIIRHDoiRC3HE2UzUcSzfpgGw4l5FsbzCrQ10saRYVKsyplpEXVMGkSHSpB/j52KTMAwDkRiOAQQ9/jkjJk67FE2XYsWbYeTZLNRxPlwDlzAFnVDC8tDSuHSNPoMuZWOzKI+YAAAFeEQIQSlZyeZTo574xLlh2xyeZ+27EkOX0mq8Djo8JKSf2KIVK/Uog0iQqV+pVDJJjaHwCAKwWiZcuWyWuvvSZr1qyRY8eOyezZs6Vr1675jtm2bZuMGjVKli5dKtnZ2RITEyPffPONREdHm/3p6eny2GOPycyZMyUjI0M6duwokyZNkoiICNt7HDx4UB555BFZvHixlC5dWnr37i3jxo0Tb2/yYHHV+OjcPrrsxZ74FDPMfXd8irmdO+ePldbwVC8XaPr+1K0Q/HcICqbzMwCgRNg1EaSmpkqjRo2kX79+0q1bt/P279mzR6699lrp37+/PP/88xIcHCxbtmwRf39/2zHDhw+XuXPnyqxZsyQkJESGDBli3uv33383+3NycqRz584SGRkpK1asMMHrgQceEB8fH3nllVdK9PO6kuycXDmWmC4HT52x3Q6cTJW9x1NNU1daVs4FX6sTH9aKDJLaEaWlVkSQ1IkMlqsiSjP0HQBgNx4W/XPeAehEeOfWEPXo0cMEl88//7zA1yQmJkr58uVl+vTpctddd5lt27dvl7p168rKlSulVatWMm/ePOnSpYscPXrUVms0efJkU+t0/Phx8fUtXA1EUlKSCVz6MzWYuTL9ldAmLJ3PJzYpTWITM+RoQpoc0dvps/da06ND3i9Ea3yiypSSGuVLmz4/el8jvLTULF9aQgJo8gIAlIzCXr8dts0oNzfX1Pw88cQTphls3bp1Uq1aNRk9erQtNGlTW1ZWlrRv3972ujp16pjmNGsg0vsGDRrka0LT99MmNK1tatKkibiDjOwcSTyTZYLOydQMOZWaaW4nUzLNTM7HkzPkeEqGxCedvc/Ms57Xhfh6e5rQEx0WIFXKBkpUWIBp9tL1viqXKSU+Xp4l8tkAALhSDhuI4uPjJSUlRcaPHy8vvfSSvPrqqzJ//nzTHKZ9ga6//nqJjY01NTyhoaH5XqvhR/cpvc8bhqz7rfsuRPsj6S1vwiwOWtuSm2sRraeziEW00iXXYjG1L1k5uZKdY5Hs3FzJyrGYldh14VFzy8oxzVJnMv++ZWRLamaOmZNHOy6npOt9tiSlZ0nCmayLNmFdSNlAX4kM8TdNXBVC/aVSaIBUKlNKKv39ODzITzwZ3QUAcAEOXUOk7rjjDtNPSDVu3Nj0A9ImLw1ExUk7XWu/peLW7rUlkpnz77UxRUGzS3ApHzNJoYYdc1/aT8oF+kr5YH8pX9pPygf5maATHuwnft706QEAuAeHDUTlypUzo8B0VFle2j9o+fLl5rF2lM7MzJSEhIR8tURxcXFmn/WY1atX53sP3W/ddyHaNDdixIh8NURRUVFS1Px9PE1/G11fy9PDQ7S+RR97e3mKt6eHaXby9vIwj7XTsS5EqkHFz8dTSvl4SSlfLwnw9ZJAX2/zWCckDPL3liA/b/M4pJSPhAb4mCCk26jRAQDAiQKRNoW1aNFCduzYkW/7zp07pUqVKuZxs2bNTKfrhQsXSvfu3c02PV6H2bdu3do81/uXX37ZNMGFh4ebbQsWLDAdq84NW3n5+fmZW3Hb+FzHYv8ZAADAgQOR9hHavXu37fm+fftk/fr1EhYWZjpGjxw5Uu655x5p27attGvXzvQh+vHHH2XJkiXmeO01rkPytSZHX6MhZ+jQoSYEaYdq1aFDBxN8evXqJRMmTDD9hp555hkZPHhwiQQeAADgBCx2tHjxYh23fd6td+/etmM++eQTS82aNS3+/v6WRo0aWb777rt875GWlmYZNGiQpUyZMpaAgADLnXfeaTl27Fi+Y/bv32/p1KmTpVSpUpZy5cpZHnvsMUtWVtYllTUxMdGUTe8BAIBzKOz122HmIXJ07jQPEQAA7nb9ZqIYAADg9ghEAADA7RGIAACA2yMQAQAAt0cgAgAAbo9ABAAA3B6BCAAAuD0CEQAAcHsEIgAA4PYIRAAAwO0RiAAAgNuz62r3zsS65JuuiQIAAJyD9br9b0u3EogKKTk52dxHRUXZuygAAOAyruO6yOuFsNp9IeXm5srRo0clKChIPDw8ijS5asg6dOjQRVfhhWPhvDknzptz4rw5pyQHOW8aczQMVaxYUTw9L9xTiBqiQtIvsXLlysX2/vrLwv/ozofz5pw4b86J8+acgh3gvF2sZsiKTtUAAMDtEYgAAIDbIxDZmZ+fnzz77LPmHs6D8+acOG/OifPmnPyc7LzRqRoAALg9aogAAIDbIxABAAC3RyACAABuj0AEAADcHoGoBEycOFGqVq0q/v7+cvXVV8vq1asvevysWbOkTp065vgGDRrITz/9VGJlxeWdt48++kiuu+46KVOmjLm1b9/+X88zHOP/N6uZM2eaWei7du1a7GXElZ+3hIQEGTx4sFSoUMGMYqpVqxb/VjrBeXv77beldu3aUqpUKTOL9fDhwyU9PV0cgo4yQ/GZOXOmxdfX1/Lpp59atmzZYhkwYIAlNDTUEhcXV+Dxv//+u8XLy8syYcIEy9atWy3PPPOMxcfHx7Jp06YSL7s7u9Tz1rNnT8vEiRMt69ats2zbts3Sp08fS0hIiOXw4cMlXnZ3dqnnzWrfvn2WSpUqWa677jrLHXfcUWLlxeWdt4yMDEvz5s0tt956q2X58uXm/C1ZssSyfv36Ei+7O5t5ieftf//7n8XPz8/c6zn7+eefLRUqVLAMHz7c4ggIRMWsZcuWlsGDB9ue5+TkWCpWrGgZN25cgcfffffdls6dO+fbdvXVV1seeuihYi8rLv+8nSs7O9sSFBRkmTZtWjGWEkVx3vRcXXPNNZaPP/7Y0rt3bwKRE5y3Dz74wFK9enVLZmZmCZYSV3re9Ngbb7wx37YRI0ZY2rRpY3EENJkVo8zMTFmzZo1pPsm7Jpo+X7lyZYGv0e15j1cdO3a84PFwjPN2rjNnzkhWVpaEhYUVY0lRFOfthRdekPDwcOnfv38JlRRXet5++OEHad26tWkyi4iIkPr168srr7wiOTk5JVhy95Z5GeftmmuuMa+xNqvt3bvXNHPeeuut4ghY3LUYnThxwvwPqv/D5qXPt2/fXuBrYmNjCzxet8Nxz9u5Ro0aZVZWPjfcwrHO2/Lly+WTTz6R9evXl1ApURTnTS+kixYtkvvuu89cUHfv3i2DBg0yf4TozMhwzPPWs2dP87prr73WrECfnZ0tDz/8sDz11FPiCKghAorY+PHjTQfd2bNnm46GcEzJycnSq1cv0yG+XLly9i4OLkFubq6p1fvwww+lWbNmcs8998jTTz8tkydPtnfRcBFLliwxNXmTJk2StWvXyrfffitz586VF198URwBNUTFSP+R9fLykri4uHzb9XlkZGSBr9Htl3I8HOO8Wb3++usmEP3666/SsGHDYi4pruS87dmzR/bv3y+33XZbvgut8vb2lh07dkiNGjVKoOTu7XL+f9ORZT4+PuZ1VnXr1jU16dqU4+vrW+zldnflLuO8jRkzxvwR8uCDD5rnOoo6NTVVBg4caAKtNrnZEzVExUj/p9S/XhYuXJjvH1x9ru3fBdHteY9XCxYsuODxcIzzpiZMmGD+0pk/f740b968hEqLyz1vOrXFpk2bTHOZ9Xb77bdLu3btzGMdEgzH/P+tTZs2ppnMGmDVzp07TVAiDDnueTtz5sx5occaah1iWVV79+p2h2GJOsxw6tSpZhj9wIEDzbDE2NhYs79Xr16WJ598Mt+we29vb8vrr79uhm8/++yzDLt3gvM2fvx4M/z066+/thw7dsx2S05OtuOncD+Xet7OxSgz5zhvBw8eNKM4hwwZYtmxY4dlzpw5lvDwcMtLL71kx0/hfmZe4nnT65metxkzZlj27t1r+eWXXyw1atQwo6sdAYGoBLz33nuW6Ohoc8HUYYqrVq2y7bv++uvNP8J5ffXVV5ZatWqZ4+vVq2eZO3euHUqNSzlvVapU0T9vzrvpPwBw7P/f8iIQOc95W7FihZmSRC/IOgT/5ZdfNlMowHHPW1ZWluW5554zIcjf398SFRVlGTRokOX06dMWR+Ch/7F3LRUAAIA90YcIAAC4PQIRAABwewQiAADg9ghEAADA7RGIAACA2yMQAQAAt0cgAgAAbo9ABAAA3B6BCIDL69Onj3h4eMjDDz983r7BgwebfXqM1cSJE6Vq1ari7+8vV199taxevbqESwygpBGIALgFXax15syZkpaWZtuWnp4u06dPl+joaNu2L7/8UkaMGCHPPvusrF27Vho1aiQdO3aU+Ph4O5UcQEkgEAFwC02bNjWh6Ntvv7Vt08cahpo0aWLb9uabb8qAAQOkb9++EhMTI5MnT5aAgAD59NNP7VRyACWBQATAbfTr10+mTJlie64hR4OPVWZmpqxZs0bat29v2+bp6Wmer1y5ssTLC6DkEIgAuI37779fli9fLgcOHDC333//3WyzOnHihOTk5EhERES+1+nz2NhYO5QYQEnxLrGfBAB2Vr58eencubNMnTpVLBaLeVyuXDl7FwuAAyAQAXC7ZrMhQ4bYRpPlpeHIy8tL4uLi8m3X55GRkSVaTgAliyYzAG7llltuMX2FsrKyzOixvHx9faVZs2aycOFC27bc3FzzvHXr1nYoLYCSQg0RALeiNUDbtm2zPT6XDrnv3bu3NG/eXFq2bClvv/22pKam5ut8DcD1EIgAuJ3g4OAL7rvnnnvk+PHjMnbsWNORunHjxjJ//vzzOloDcC0eFu1ZCAAA4MboQwQAANwegQgAALg9AhEAAHB7BCIAAOD2CEQAAMDtEYgAAIDbIxABAAC3RyACAABuj0AEAADcHoEIAAC4PQIRAABwewQiAAAg7u7/AQ0Z/IG/FrAfAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Compute thrust evolution with velocity\n", "mach_array = np.linspace(0, max_mach, nb_points)\n", "thrust_array = np.zeros(nb_points)\n", "for i, mach in enumerate(mach_array):\n", " turbojet.M0 = mach\n", " turbojet.tune_current_OPR()\n", " thrust_array[i] = turbojet.thrust\n", "\n", "# Plot thrust evolution with velocity\n", "plt.plot(mach_array, thrust_array)\n", "plt.xlabel(\"M0\")\n", "plt.ylabel(\"Thrust [N]\")\n", "plt.title(\"Thrust evolution with Mach\")\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Altitude influence\n", "\n", "The altitude also has an influence on the thrust of the engine." ] }, { "cell_type": "code", "execution_count": 65, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk0AAAHHCAYAAACiOWx7AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAXxBJREFUeJzt3Qdc1PX/B/AXewqIAwRRceLee5RpzixHmiN3Wq6yYerfzLZmU0szG47clnv1M3duUdxbFDcqyhDZ3//j/am7OEQFBb43Xs/H4zru7tPd576A9+Iz3l87TdM0EBEREdFD2T/8YSIiIiISDE1EREREWcDQRERERJQFDE1EREREWcDQRERERJQFDE1EREREWcDQRERERJQFDE1EREREWcDQRERERJQFDE1EeWjz5s2ws7PD77//rndXLM4HH3ygjl1OK1GiBPr06QNL8fTTT6tLVttWqlTJrPokx1qOeV44f/68+pmZOXNmnrweWT+GJqInJP8oZ+UigckcrVmzRgUSa7Zjxw71Hu/cuQNrc+XKFfXewsLCYIl9io+PV23N9feDKD1Hk1tElG2//fabye3Zs2dj/fr1991fvnx5HD9+HOYYmqZMmWLVwUlC04cffqhGOXx8fEweO3nyJOztLefvx//973/3BRR5bzJ6U61aNZiDh/Xpp59+QlpamklokrYiq6NVRHphaCJ6Qi+//LLJ7V27dqnQlPF+8aShST5g3N3dn+g5yJSLiwssibOzMyyZk5OT3l0gemyW8+cVkRWRv7Q//fRTFC1aFK6urmjWrBnOnDmT6XqU0NBQNGnSRIWl//u//1OPyXRfZiNDGdfnJCcnq7/iy5Qpo16nQIECaNSokQp1QtrKKJPhOQ2XR1m7di0aN24MDw8P5MuXD23btsXRo0eNj3/55ZfqeS5cuHDf/zt69Gj1wX/79m3jfYsXL0bNmjXh5uaGggULqsB5+fLlx16vkv74yPWIESPU18HBwcb3KP9/ZsdMnDt3Dp07d4avr6867vXq1cPq1aszXZ+2aNGiR34vMzp06JD6f1esWGG8T77Pcl+NGjVM2rZu3Rp169bNdP2Q9KF27drq6759+xrfW8ZjcuzYMTRt2lS9l8DAQEycOBFZMWPGDDzzzDMoXLiwCpcVKlTADz/88ND/51F9Sr+mSb4HhQoVUl/Lz6mhreF796C1Upmti5KpV7nf29tbjSb27t37gdOxJ06cwIsvvqi+v/I9q1Wrlsn3guhBONJEpIMJEyaoKaF33nkH0dHR6kOsR48e2L17t0m7W7duqQ/Nrl27qiDh5+eXrdeRD5/x48fjlVdeQZ06dRATE4N9+/Zh//79ePbZZ/Hqq6+qqZTMphMfRNrJB1LLli3x+eefq9Ev+SCVMHbgwAH1YdalSxe8++67KlAYAouB3NeiRQvkz59f3ZYPU/lwlQ9a6ev169cxadIkbN++XT1fxum07OrYsSNOnTqF+fPn45tvvlGhTBg+rDOS12/QoIF6X6+//roKmrNmzcLzzz+vFvB36NDhsb6X6UkYlve1detW9bxi27Zt6nkOHjyovk9eXl4qXMvU4sCBAzN9Hpny/eijj/D++++rNhJkhfTfQMJpq1at1HGQ74u8h5EjR6Jy5crqZ+th5PtasWJF1UdHR0esXLkSgwcPVv0aMmTIY/fJQL4H8hqDBg1Sx1X6KKpUqYLs0DQNL7zwAv7++2+89tprqg9Lly5VP6cZSbhv2LChCo+jRo1SwV9+Jtu3b48//vjjvu8vkQmNiHLUkCFDtAf9am3atEk9Vr58eS0xMdF4/6RJk9T9hw8fNt731FNPqfumTZt23/PI/ePGjbvv/uLFi2u9e/c23q5atarWtm3bx+5vRrGxsZqPj482YMAAk/uvXbumeXt7m9xfv359rWbNmibt9uzZo15r9uzZ6nZSUpJWuHBhrVKlStq9e/eM7VatWqXavf/++8b75P2m72d4eLi6PWPGjEceny+++ELdJ//Po47Z8OHDVdtt27aZvO/g4GCtRIkSWmpqara/l5mR70udOnWMtzt27KguDg4O2tq1a9V9+/fvV8+1fPlyk58LuRjs3bv3gcfB8DNkON5C+urv76916tRJe5T4+Pj77mvZsqVWsmTJ+14nq32SYy3H3ODGjRsP/HnO+LwPeo5ly5ap55g4caLxvpSUFK1x48b39aNZs2Za5cqVtYSEBON9aWlpWoMGDbQyZco89HgQcXqOSAcyspJ+bYrhr3GZFkpPpkSk7eOS0Qz5y/r06dPICTIiJVMe3bp1w82bN40XBwcHNYW0adMmY9uXXnpJTTmdPXvWeN/ChQvVe5JRASGjXpGRkWr0QqZJDGS6LyQk5L4psbxaGC+jcjJyZuDp6alGTWQ6Saa6Hud7mZG0kxG/u3fvqtsyStKmTRu1cFpGnYRcy3RV+r5kl/Q9/fo66au8v0f1T8h0qYGMosn3+qmnnlL/r9w2F/I9k5EwGbEykJ/JYcOGmbSLiorCxo0b1YhbbGys8edXRnRl5FR+Tx41LUy2jaGJSAfFihUzuW2Yqkq/zkfIFMKTLPyVaRIJOWXLllXTMTJVJutpHpchfMk6F5laSX+RXV0SgAxkTZBMN0lQEjIAJGuXZEpIpp6EYc1TuXLl7nstCU2ZrYnKbfKamfVHpnwMjz/O9zKz0JSSkoKdO3eqHXxy7OQ+Wb+WPjTJOiJZe/O4ZK1VxnVq0sdH9U/IFGnz5s3VFJYEcPk+G9bVmVNoku9JkSJFVEBML+P3Udaayc/h2LFj7/v5HTdunGqT/meYKCOuaSLSgfwVnJl/ZpYy/0s/K1JTU01uywewjPQsX75chZqff/5ZreuZNm2aWueUXYat4rKuyd/f/77H5a99g4CAABUCZL2IfNDKrsKIiAi1DionPGjBesZjYC7fy4xk8bGMrsm6Jglesthawq0cs6lTpyIxMVGFpiddY/O4/ZOfG1nULuH166+/RlBQkArwMqojP0PpywbkFvkeZ9bPx/0eG/os689kZCkzpUuXfqznJtvA0ERkgWSkIOPOoKSkJFy9evW+tjJKIVNIcomLi1NBShaIG0JTdqpslypVSl3LB7yMQDyKTNHJ1JuMpMiIk+zeateunfHx4sWLq2t5XEav0pP7DI9nxjCik/E4ZDY6lZ33KK8pr53Zjqv0fX5ShmkyCUYSmgzTenItgWnu3LlqUbp8vx4mN6qkC1n0Lf2QXWXpR9PST8HmRJ8e1la+x5lNI2b8Hsv3ZMOGDernO/1oU8bvY8mSJY1lD7Ly80uUEafniCyQhBcZoUhv+vTp9/0FLms10pMPFPlLWj4MDWTqRWSlWrb8dS5Ta5999pkqZ5DRjRs3TG536tRJjXTIzjWZmnvuueeMr2cYbZEAJiNf6fskJQ2kppWsbXoQ6YfshMt4HGSUJqPsvEdZV7Rnzx41bWYg647k+MrOQJkuyykSkGSXnQQRQ2iS9yRTgYYROcP9D5Kd9/Y4I1TpR3pkSk7KEDxKdvpkqDuWWVv5OZewmv7nSnYXyrRhxu+ZTHWmL4cgvwvfffedSTv5WZMSBj/++GOmf2Bk/PklyogjTUQWSEaJZGu1hBIpHSAfJH/++adxO72BfMDLh4TUQJIRJ1l4LVvOhw4damwjjwnZXi+hSD4spcTBg4KKfDD17NlT1ROSdrIeRKbdZNG2bOX+/vvvTT6kpD6QTO/IwlsZeUpP/uKXcCCjYLLAWBaYG0oOSEB58803H3kcZMu/XEsAkwAl5QUyMrzHMWPGqD7L68qIV/oAZyDb0CXkydorOSZy3KTkQHh4uNqSnpPVwyUQSY2nixcvmoQjGV2SD3Y5BrIm6WEkWMh6IwmeUjNL3pMsypeaVE9CykLIaJgcJylNIaM4Us1bvqeZBY7H7ZNMQcvPqYxEyvSkHG8pySCXfv36qZ8d+bns37+/Wm8kzyllEKQsg4H0UX725Hsni/Xl+ZYsWZLpuiupSyYL62WN34ABA9Tok/zMSUi+dOmS+l0ieiC9t+8R2WLJgcWLF5vcn9n2edlqXbFixUyfR7a9jxw5UitYsKDm7u6utoGfOXPmvu3zn3zyidrWLmUC3NzctJCQEO3TTz9VW/3Tb80eNmyYVqhQIc3Ozi5L5QfkfchrSpkBV1dXrVSpUlqfPn20ffv23df2p59+Us+ZL18+k7IC6S1cuFCrXr265uLiovn6+mo9evTQLl26ZNImY8kBw5b4/v37q37I83fp0kWLjIzMdAv7xx9/rAUGBmr29vYm5QcyHjNx9uxZ7cUXX1THTd6fHEMpg/C438sHiYmJUSUGpO/yfTCYM2eOeo6ePXtmaRu+lCSoUKGC5ujoaPLaD/oZyrhl/0FWrFihValSRR0DKbfw+eefa7/++ut95Ruy06fMXnvHjh2qPIWzs/N93zs5FlLiQB6rVq2a9ueff2b6HLdu3VLHy8vLS/08yNcHDhzI9Hsh399evXqp0gtOTk7q5+K5557Tfv/990ceE7JtdvKfB0cqIiIiIhJc00RERESUBQxNRERERFnA0ERERESUBQxNRERERFnA0ERERESUBQxNRERERFnA4pY5RM5pdOXKFVXILbdOa0BEREQ5SyovSfFdOV/mo4rXMjTlEAlMckJLIiIisjxSmf9RFfgZmnKIjDAZDrqcaoKIiIjMn5ySRwY9DJ/jD8PQlEMMU3ISmBiaiIiILEtWltZwITgRERFRFjA0EREREWUBQxMRERFRFjA0EREREZl7aNq6dSvatWunaiPIAqxly5YZH0tOTsbIkSNRuXJleHh4qDa9evVSW/vTi4qKQo8ePdTiax8fH/Tv3x9xcXEmbQ4dOoTGjRvD1dVVrZCfOHHifX1ZvHgxQkJCVBt5zTVr1uTiOyciIiJLo2tounv3LqpWrYopU6bc91h8fDz279+PsWPHquslS5bg5MmTeP75503aSWA6evQo1q9fj1WrVqkgNnDgQJOthC1atEDx4sURGhqKL774Ah988AGmT59ubLNjxw5069ZNBa4DBw6gffv26nLkyJFcPgJERERkKew0KYVpBmSkaenSpSqsPMjevXtRp04dXLhwAcWKFcPx48dRoUIFdX+tWrVUm3Xr1qFNmza4dOmSGp364YcfMGbMGFy7dg3Ozs6qzahRo9So1okTJ9Ttl156SQU4CV0G9erVQ7Vq1TBt2rQs9V/Cmbe3N6Kjo1lygIiIyEJk5/PbotY0yRuScCXTcGLnzp3qa0NgEs2bN1dl0Hfv3m1s06RJE2NgEi1btlSjVrdv3za2kf8vPWkj9z9IYmKiOtDpL0RERGS9LCY0JSQkqDVOMo1mSIIyelS4cGGTdo6OjvD19VWPGdr4+fmZtDHcflQbw+OZGT9+vEqmhgtPoUJERGTdLCI0yaLwLl26qJPqyXSbORg9erQa+TJc5PQpREREZL0cLSUwyTqmjRs3msw3+vv7IzIy0qR9SkqK2lEnjxnaXL9+3aSN4faj2hgez4yLi4u6EBERkW2wt4TAdPr0afz1118oUKCAyeP169fHnTt31K44AwlWaWlpqFu3rrGN7KiT5zKQnXblypVD/vz5jW02bNhg8tzSRu4nIiIi0j00ST2lsLAwdRHh4eHq64iICBVyXnzxRezbtw9z585FamqqWmMkl6SkJNW+fPnyaNWqFQYMGIA9e/Zg+/btGDp0KLp27ap2zonu3burReBSTkBKEyxcuBCTJk3CW2+9ZezHG2+8oXbdffXVV2pHnZQkkNeV5zIHu87dQkzCf6GPiIiIdKDpaNOmTVLu4L5L7969tfDw8Ewfk4v8fwa3bt3SunXrpnl6empeXl5a3759tdjYWJPXOXjwoNaoUSPNxcVFCwwM1CZMmHBfXxYtWqSVLVtWc3Z21ipWrKitXr06W+8lOjpa9U2uc9KcXee14FGrtFdm7dVSU9Ny9LmJiIhsXXQ2Pr/Npk6TpcutOk0HL95B52k7kZSahndalMXQZ8rk2HMTERHZuhhrrdNki6oG+eDj9hXV11+tP4XNJ00XvhMREVHeYGiyAC/VLoZudYpBxgTfWBCGiFvxeneJiIjI5jA0WYgPnq+AakE+iL6XjFfnhOJeUqreXSIiIrIpDE0WwsXRAT+8XAMFPZ1x/GoMRi85pIp9EhERUd5gaLIgRbzd8H33GnCwt8OysCuYsf283l0iIiKyGQxNFqZeyQL4vzbl1defrjmuajgRERFR7mNoskD9GpZA+2oBSE3TMGTufly5c0/vLhEREVk9hiYLZGdnh/Edq6BigBdu3U3Ca3NCkZDMheFERES5iaHJQrk5O2DayzWR390Jhy5F471lR7gwnIiIKBcxNFmwIF93tTDc3g74PfQSZu+8oHeXiIiIrBZDk4VrWLogRrf+Z2H4x6uOYTcXhhMREeUKhiYr8ErjYLxQLQApaRoGz92Py1wYTkRElOMYmqxkYfiEdAvDX/1tHxeGExER5TCGJitaGP5jz5rw9XDGkcsxGPUHK4YTERHlJIYmK1I0vzumpKsY/svf4Xp3iYiIyGowNFmZ+qUKYGzbfxaGf7bmOLadvqF3l4iIiKwCQ5MV6t2gBDrXLIo0DRg67wDO37yrd5eIiIgsHkOTlS4M/6RDJVQv5oPoe8l4ZfY+xCYk690tIiIii8bQZKVcHB3w48s14eflgjORcXhzYRjSZOiJiIiIHgtDkxUr7OWK6T1rwdnRHn8dj8TX60/p3SUiIiKLxdBk5aoG+WBCx8rq6+83ncGqQ1f07hIREZFFYmiyAR1rFMWAxsHq63cWH8SRy9F6d4mIiMjiMDTZiFGty6NJ2UJISE7DgNn7cCM2Ue8uERERWRSGJhshBS+/61YdJQt54Gp0gjrVSmIKT7VCRESUVQxNNsTbzQk/96oFL1dH7I+4gzFLj/BUK0RERFnE0GRjShbyxPfda8DeDvg99BJPtUJERJRFDE02SNY2vde2gvFUK5tORurdJSIiIrPH0GSj+jYsgZdqBalTrbw+7wBOX4/Vu0tERERmjaHJhk+18nH7SqgT7IvYxBT0n7UPUXeT9O4WERGR2WJosmFSKXzayzUR5OuGiKh4DJoTiqSUNL27RUREZJYYmmycr4czfuldG54ujtgdHoVxK7ijjoiIKDMMTYSyfvkwuVs12NkB8/dcxK/bz+vdJSIiIrPD0ETKMyF+GNOmvPr609XHsOkEd9QRERGlx9BERv0bBRt31A2bfwAnr3FHHRERkQFDE923o65usC/iElPQb+Ze3IzjOeqIiIgEQxNluqOuRAF3XL5zDwNn70NCMs9RR0RExNBE98kvO+r61Daeo27kH4e4o46IiGweQxNlqlQhT/zwck042NthedgVfLfxjN5dIiIi0hVDEz1Qw9IF8fELldTXX68/hRUHr+jdJSIiIt0wNNFDda9bDAMaB6uv31l8EKEXbuvdJSIiIl0wNNEjjWpdHs9W8FOnWJGF4Rej4vXuEhERUZ5jaKJHknVNk7pWQ8UAL9y6m4S+M/ci+l6y3t0iIiLKUwxNlCXuzo7qHHX+Xq44ExmHIXP3IzmVJ/clIiLbwdBEWebv7Yqfe9eCu7MD/j5zE2OX8eS+RERkOxiaKFsqBXrju27VYW8HLNh7EdO2nNO7S0RERHmCoYmyrVl5P7z/XAX19efrTmDVIZYiICIi68fQRI+lT8Ng9G1YQn391iKWIiAiIuvH0ESP7b22FdC8fGFVimDA7H24cOuu3l0iIiLKNQxN9ISlCKqjUqAXov4tRXAnPknvbhEREeUKhiZ6Ih4u/5QiCPB2xbkbdzFwdigSklP17hYREVGOY2iiJ+bn5YoZfesgn4sj9pyPwojfDyEtjaUIiIjIujA0UY4o558P03rWhKO9HVYevIIv/ndS7y4RERHlKIYmyjENSxfEhE5V1Nc/bD6Lebsj9O4SERFRjmFoohz1Ys2ieKNZGfX12OVHsOlEpN5dIiIiyhEMTZTjhjcvg041iiI1TcOQeftx+FK03l0iIiJ6YgxNlOPs7OwwvmNlNC5TEPFJqaoUwcWoeL27RURE9EQYmihXODvaY2qPGihfxAs34xLRe8Ye3L7LGk5ERGS5GJoo1+RzdcKMPv/VcJKq4azhREREloqhiXKVv7crZvarg3yujth34TbeWhTGGk5ERGSRdA1NW7duRbt27RAQEKDWwSxbtszkcU3T8P7776NIkSJwc3ND8+bNcfr0aZM2UVFR6NGjB7y8vODj44P+/fsjLi7OpM2hQ4fQuHFjuLq6IigoCBMnTryvL4sXL0ZISIhqU7lyZaxZsyaX3rXtKeuXD9N71oKTgx3WHL6Gj1YdU99bIiIiS6JraLp79y6qVq2KKVOmZPq4hJvJkydj2rRp2L17Nzw8PNCyZUskJCQY20hgOnr0KNavX49Vq1apIDZw4EDj4zExMWjRogWKFy+O0NBQfPHFF/jggw8wffp0Y5sdO3agW7duKnAdOHAA7du3V5cjR47k8hGwHfVLFcBXXaqpr2fuOI+ftp3Tu0tERETZo5kJ6crSpUuNt9PS0jR/f3/tiy++MN53584dzcXFRZs/f766fezYMfX/7d2719hm7dq1mp2dnXb58mV1e+rUqVr+/Pm1xMREY5uRI0dq5cqVM97u0qWL1rZtW5P+1K1bV3v11Vez3P/o6GjVF7mmB/tp61mt+MhV6rLswCW9u0NERDYuOhuf32a7pik8PBzXrl1TU3IG3t7eqFu3Lnbu3Kluy7VMydWqVcvYRtrb29urkSlDmyZNmsDZ2dnYRkarTp48idu3bxvbpH8dQxvD62QmMTFRjWKlv9CjvdK4JPo3ClZfv7P4ILafual3l4iIiLLEbEOTBCbh5+dncr/cNjwm14ULFzZ53NHREb6+viZtMnuO9K/xoDaGxzMzfvx4FeIMF1krRVkzpk15tK1SBMmpGl79LRRHr7D4JRERmT+zDU3mbvTo0YiOjjZeLl68qHeXLIa9vR2+7lIV9Ur6Ii4xBX1msPglERGZP7MNTf7+/ur6+vXrJvfLbcNjch0ZaXpus5SUFLWjLn2bzJ4j/Ws8qI3h8cy4uLioHXvpL5R1Lo4O+LFnLYT458ON2ET0/nUPbsUl6t0tIiIiywtNwcHBKrRs2LDBeJ+sG5K1SvXr11e35frOnTtqV5zBxo0bkZaWptY+GdrIjrrk5GRjG9lpV65cOeTPn9/YJv3rGNoYXodyh7ebE2b1q4NAHzecu3kX/WbtQ3xSit7dIiIiMr/QJPWUwsLC1MWw+Fu+joiIUHWbhg8fjk8++QQrVqzA4cOH0atXL1XTScoBiPLly6NVq1YYMGAA9uzZg+3bt2Po0KHo2rWraie6d++uFoFLOQEpTbBw4UJMmjQJb731lrEfb7zxBtatW4evvvoKJ06cUCUJ9u3bp56Lcpefl6sKTj7uTjh48Q4Gz92P5NQ0vbtFRER0P01HmzZtUtv8Ml569+5tLDswduxYzc/PT5UaaNasmXby5EmT57h165bWrVs3zdPTU/Py8tL69u2rxcbGmrQ5ePCg1qhRI/UcgYGB2oQJE+7ry6JFi7SyZctqzs7OWsWKFbXVq1dn672w5MCTCb0QpZV7b40qRfDWwjD1vSciIspt2fn8tpP/ZJKlKJtk6lB20cmicK5vejwbT1zHgNmhSE3T8OpTJTG6dXm9u0RERFYuJhuf32a7polszzMhfpjQsbL6+sct5/DTVlYNJyIi88HQRGalc60gjGodor7+dM1x/BF6Se8uERERKQxNZHZebVISr/xbNfzdPw5h0wnTshJERER6YGgisyM7J/+vTXl0qB6o1jcNmhuK0Av/nPKGiIhILwxNZLZVwye+WAVPlyuEhOQ09Ju5FyevxerdLSIismEMTWS2nBzs8UOPmqhZPD+i7yWj5y+7eboVIiLSDUMTmTU3Zwf82rs2yvnlQ2RsogpOctoVIiKivMbQRGbP290Js/vXQdH8bjh/K16dpy4m4b/T4hAREeUFhiaymNOtzOlfFwU9nXHsagxembUPCcmpeneLiIhsCEMTWYwSBT3UeeryuThiT3gUhvA8dURElIcYmsiiVAzwxi99asPF0R4bTkTincUHkZbGMwEREVHuY2gii1Mn2BfTXq4JR3s7LA+7gnErjsqJp/XuFhERWTmGJrJITUMK46suVWFnB/y26wK+Xn9K7y4REZGVY2gii/VCtUB8/EIl9fV3G8/wBL9ERJSrGJrIor1crzhGtCxnPMHvgj0ReneJiIisFEMTWbzBT5fCq0+VVF+PXnoYKw9e0btLRERkhRiayCpO8DuqVQh61C0GWQ/+5sIwbDxxXe9uERGRlWFoIqsJTrK+6YVqAUhJ0zBozn7sPHtL724REZEVYWgiq2Fvb4cvO1dF8/J+SExJwyuz9uJAxG29u0VERFaCoYmsipODPb7vXh0NShXA3aRUdZ66Y1di9O4WERFZAYYmsjquTg74qVct1Cjmg5iEFPT8ZTfORMbp3S0iIrJwDE1klTxcHDGjbx1UDPDCrbtJePnn3bgYFa93t4iIyIIxNJHV8nZzwm/966JMYU9ci0lA95934Wr0Pb27RUREFoqhiayar4cz5r5SF8ULuONi1D30+Hk3bsQm6t0tIiKyQAxNZPUKe7mq4BTg7YpzN+6qNU637ybp3S0iIrIwDE1kE4rmd8e8AfVQOJ8LTlyLRc9fdyP6XrLe3SIiIgvC0EQ2o0RBD8wbUBcFPJxx5HIM+s7Yg7jEFL27RUREFoKhiWxK6cL51OJwWSS+P+IO+s/ci3tJqXp3i4iILABDE9mcCgFemN2vDvK5OGJ3eBQG/rYPCckMTkRE9HAMTWSTqgb5YGa/2nB3dsC20zcxaE4oElMYnIiI6MEYmshm1Szui1/71Iarkz02nbyBIXMPICklTe9uERGRmWJoIptWr2QB/NyrNpwd7fHX8et4Y8EBpKQyOBER0f0YmsjmNSpTENN71oSzgz3WHrmGNxcdZHAiIqL7MDQRAXi6XGFM7VEDjvZ2WHnwCt5ZfBCpaZre3SIiIjPC0ET0r+YV/PB99xpwsLfDsrArGPE7gxMREf2HoYkonVaV/PFdt+oqOC3Zfxmj/jiENAYnIiJiaCK6X5vKRfDtS9VgbwcsDr2E/1t6mMGJiIgYmogy065qAL75Nzgt2HsR7y0/wuBERGTjGJqIHuCFaoH4qktV2NkB83ZHMDgREdk4hiaih+hQvSi+6szgREREgKPeHSAydx1rFFXXby8+qIKTHYCPX6gEe5m7IyIim8GRJqIsBqcvX/xnxGkuR5yIiGwSQxNRFnWq+V9wkhGnMcsYnIiIbAlDE1E2g5NhjdP8PREsR0BEZEMYmogeY6ru6y5VjeUIRi1hAUwiIlvA0ET0mLvqDHWcFu27hHf/OMRTrhARWTmGJqInqOM0qes/p1z5PfSSOslvSmqa3t0iIqJcwtBE9ISVwyf/G5yWHriMNxcxOBERWSuGJqIn1LZKEUzpXgNODnZYefAKhs0/gGQGJyIiq8PQRJQDWlXyxw89asLZwR5rj1zD4Ln7kZiSqne3iIgoBzE0EeWQ5hX88GOvmnB2tMf6Y9cxaM5+JCQzOBERWQuGJqIc1LRcYfzauzZcneyx8UQkBszeh3tJDE5ERNaAoYkohzUqUxAz+tSBu7MDtp2+iT4z9iAuMUXvbhER0RNiaCLKBfVLFcBv/esgn4sjdodHodcvuxGTkKx3t4iI6AkwNBHlkprFfTF3QF14uzlhf8Qd9PhpN27fTdK7W0RE9JgYmohyUZWiPpg/oB4KeDjj8OVodPtpF27EJurdLSIiegwMTUS5rEKAFxYMrIfC+Vxw4losXpq+E1ej7+ndLSIiyiaGJqI8UMYvHxa9Wh+BPm44d+MuOk/biYhb8Xp3i4iIsoGhiSiPlCjogUWv1UeJAu64dPseuvy4E2ci4/TuFhERWUNoSk1NxdixYxEcHAw3NzeUKlUKH3/8MTTtv7PJy9fvv/8+ihQpoto0b94cp0+fNnmeqKgo9OjRA15eXvDx8UH//v0RF2f6YXXo0CE0btwYrq6uCAoKwsSJE/PsfZLtkJEmGXEq6+eJazEJ6Dp9J45didG7W0REZOmh6fPPP8cPP/yA77//HsePH1e3Jcx89913xjZye/LkyZg2bRp2794NDw8PtGzZEgkJCcY2EpiOHj2K9evXY9WqVdi6dSsGDhxofDwmJgYtWrRA8eLFERoaii+++AIffPABpk+fnufvmaxfYS9XLBhYH5UCvXAzLkkFp9ALt/XuFhERPYKdln7Yxsw899xz8PPzwy+//GK8r1OnTmpEac6cOWqUKSAgAG+//Tbeeecd9Xh0dLT6f2bOnImuXbuqsFWhQgXs3bsXtWrVUm3WrVuHNm3a4NKlS+r/l2A2ZswYXLt2Dc7OzqrNqFGjsGzZMpw4cSJLfZXg5e3trV5fRrSIHkXqNvWbsRf7LtxWhTB/6lULDUsX1LtbREQ2JSYbn99mPdLUoEEDbNiwAadOnVK3Dx48iL///hutW7dWt8PDw1XQkSk5A3njdevWxc6dO9VtuZYpOUNgEtLe3t5ejUwZ2jRp0sQYmISMVp08eRK3b3MEgHKHl6sTZvevg8ZlCiI+KRV9Z+7FX8eu690tIiKyxNAkoz0yWhQSEgInJydUr14dw4cPV9NtQgKTkJGl9OS24TG5Lly4sMnjjo6O8PX1NWmT2XOkf42MEhMTVTpNfyHKLndnR/zcuxZaVvRDUkoaXp0TiuVhl/XuFhERZcIRWSABIzvs7Oywf/9+tUboSSxatAhz587FvHnzULFiRYSFhanQJFNqvXv3hp7Gjx+PDz/8UNc+kHVwcXTAlO418O7vh7DkwGUMXxiGmIQU9Kz3ZL8/RESkQ2i6c+cOvv32WzX19Siyzmjw4MFq59uTGjFihHG0SVSuXBkXLlxQgUVCk7+/v7r/+vXravecgdyuVq2a+lraREZGmjxvSkqK2lFn+P/lWv6f9Ay3DW0yGj16NN566y3jbRlpkl13RI/D0cEeX3auCk9XR8zeeQFjlx1BzL1kDH66lPojhIiILCQ0CQkuGae5HmTYsGHICfHx8WrtUXoODg5IS0tTX0spAgk1su7JEJIkvMhapUGDBqnb9evXV6FPdsXVrFlT3bdx40b1HLL2ydBGFoInJyeraUAhO+3KlSuH/PnzZ9o3FxcXdSHKKfb2dvjw+YpqrdP3m87giz9PquA0qnUIgxMRkaWsaZKAkdXAJGJjY1GyZEk8qXbt2uHTTz/F6tWrcf78eSxduhRff/01OnTooB6XDxKZrvvkk0+wYsUKHD58GL169VLTd+3bt1dtypcvj1atWmHAgAHYs2cPtm/fjqFDh6oQKO1E9+7d1SJwqd8kpQkWLlyISZMmmYwkEeUF+Zl+p2U5jGlTXt3+ces5/N/Sw0hNM9tNrkRENsOsSw5I+JLilhKWZIpNQk63bt1UMUvDTjfp/rhx41RNJRlRatSoEaZOnYqyZcsan0em4iQorVy5Uo1cSdkCqe3k6elpUtxyyJAhqjRBwYIF1WjZyJEjs9xXlhygnLZwbwRGLzkMyUttKvvjm5eqqfVPRESUc7Lz+Z3l0CQFIbNCtu7bIoYmyg1rDl/F8AVhSEpNQ6PSBfFjz5rwcMnyrDoREekRmjKuLTJ5kn/XW8i1LLK2RQxNlFv+Pn0TA3/bp2o5VQ3ywcw+tZHf47+aYkREZGbFLaXIY2aXy5cvq11usiha6ikRUc5qVKYg5g2oBx93Jxy8eAedf9yJq9H39O4WEZHNyXJokhSW/pIvXz4sXrwYderUwfz58zFlyhS1LoiIcl61IB/8/lp9+Hu54kxkHDpN3aGuiYjIzCuCL1myRJ3PTRZKv/HGG+o0J3379n3oFB4RPZnShfPh90H1UbKgB65EJ6DztB0Iu3hH724REdmMbKWcLVu2oF69eujZsyc6duyIc+fOqRPlsl4RUd4omt8di1+rj6pFvXE7Phndf9qFradu6N0tIiKbkOXQ1KZNGzz77LOqiOTZs2fx2WefZalCOBHlrAKeLmqNk+FEv/1n7eX56oiI8kC2ds/JiW49PDweWp1YaiLZIu6eo7wmJ/h9a1EYVh26qm6Pa1cBfRsG690tIiKr/fzOcsGXGTNm5ETfiCiHODvaY3LX6ijo6YKZO87jw5XHEBmbiHdbluNpV4iIbK0iuCXhSBPpRX6Fp24+q85VJzrWCMTnnarAyYEbM4iIdKnTRETmSUaVhjQtjYkvVoGDvR2W7L+MAbOlGKZtFpolIsotWQpNvr6+uHnzZpaftFixYrhw4cKT9IuIsqlLrSBM71kTrk722HzyBrpN34VbcYl6d4uIyGpkaU2TnAh37dq1Wd4td+vWLaSmpj5p34gom5qV98PcV+qpHXUHL0Wj0w87MKtfHRQv4KF314iIbGNN0+MUrTxz5gxKliwJW8E1TWROzt6IQ+9f9+DS7Xso4OGMGX1ro0pRH727RURk/Wua0tLSsn2xpcBEZG5KFfLEksENUDHAC7fuJqHr9F3YdDJS724REVk0LgQnslKF87li4av1jUUwX5m1D4v2XtS7W0REFouhiciKebo44tc+tVUZgtQ0De/+cQjfrD+lyhQQEVH2MDQRWTmp1/RV56oY2rS0uj1pw2mM+P0QklPT9O4aEZFFYWgispFaTu+0LIfPOlRWtZx+D72EfjP3IjYhWe+uERFZDIYmIhvSvW4x/NyrFtycHLDt9E10+XEXrkbf07tbRETWGZocHBwQGRmZaW0meYyIzFvTkMJY+Go9dc6641dj0GHKDnVNREQ5HJoetIA0MTERzs7O2X06ItKB1GxaOrgBShf2xLWYBHSethNbTt3Qu1tERJZfEVxMnjzZuDbi559/hqenp/Exqf69detWhISE5E4viSjHBfm644/XGuDVOfuw61yUWuP0aftK6FqnmN5dIyKy3IrgIjg4WF3LOeWKFi1qMhUnI0wlSpTARx99hLp168IWsSI4WarElFSM+uMwlh64rG4PfroU3mlRDvb2dnp3jYjIrD6/szzSFB4erq6bNm2KJUuWIH/+/E/eUyLSnYujA77uUhVB+d0weeMZTN18Fhei4lWZAlcnrlMkInrsNU2bNm0yCUwyNRcWFobbt29n96mIyEzItPtbLcrhy85V4eRgh9WHrqL7T7twKy5R764REVluaBo+fDh++eUXY2Bq0qQJatSogaCgIGzevDk3+khEeeTFmkUxq18deLk6Yn/EHXSYugNnIuP07hYRkWWGpsWLF6Nq1arq65UrV+L8+fM4ceIE3nzzTYwZMyY3+khEeahBqYJYMrghgnzdEBEVj45Tt2PHmZt6d4uIyPJCk9Rj8vf3V1+vWbMGnTt3RtmyZdGvXz8cPnw4N/pIRHlMShEsG9wQNYr5ICYhBb1+3YOFeyP07hYRkWWFJj8/Pxw7dkxNza1btw7PPvusuj8+Pp7FLYmsSAFPF8wbUA/PVw1ASpqGkX8cxvg1x5GWxpP9EpFtynZo6tu3L7p06YJKlSqpxaPNmzdX9+/evZt1moisjOyem9S1Gt5oVkbd/nHrObw2JxTxSSl6d42IyHzrNKX3+++/4+LFi2pqTmo2iVmzZsHHxwcvvPACbBHrNJG1Wx52GSMWH0JSahoqBnjh5961UMTbTe9uERHl2ef3Y4Umuh9DE9mC0AtRGDg7FLfuJqFwPhf81KsWqgb56N0tIiLzDE1S9fth3n//fdgihiayFRej4vHKrH04eT0WLo72+KpLVTxXJUDvbhERmV9oql69usnt5ORkVS3c0dERpUqVwv79+2GLGJrIlsQlpuD1+Qew8USkuj28eRm17knWORIRwdZPo2Jw4MCBTF+wT58+6NChQ3afjogskKeLo5qam7D2OH7aFo5v/zqN05Fx+PLFqnBz5i5aIrJOObamSWo0tWvXThW7tEUcaSJbJfWb3lt2BMmpGioFeqkwxQXiRGSNn9/ZLjnwIPJiciEi2/JS7WKY+0o9+Ho448jlGDz//XYciOC5KInI+mR7em7y5Mkmt2Wg6urVq/jtt9/QunXrnOwbEVmIOsG+WD6koXGB+EvTd+HzTpXRofo/JUmIiGxyei44ONjktr29PQoVKoRnnnkGo0ePRr58+WCLOD1H9M8C8eELwvDX8evq9qtNSuLdViFwsOcCcSIyT6zTpAOGJqJ/yGlWvlp/ElM2nVW3ny5XCJO6Voe3m5PeXSMiyrs1TVJeQEoLHDlyJDv/GxHZEHt7O4xoGYLvulWHq5M9Np+8gQ5Tt+PsjTi9u0ZE9ESyFZqcnJxQrFgxdbJeIqKHaVc1AL+/1gAB3q44d+Mu2k/Zjk3/1nUiIrJE2d49N2bMGPzf//0foqKicqdHRGQ1KgV6Y/nQRqhVPD9iE1LQb9ZeTN18Rm0gISKyNI9VEfzMmTNqqq548eLw8PAweZwVwbmmiSijpJQ0jFtxFPP3RKjbbasUwRcvVoG7c7Y38BIRWU5F8Pbt2z9J34jIBjk72mN8x8qq+OW45Uex+tBVnI2MU4Uwg3zd9e4eEVGWcPdcDuFIE1HW7D0fhUFzQnEzLgk+7k74vlsNNCpTUO9uEZGNismLkgNJSUmIjIxEWlqayf2yUNwWMTQRZd2VO/fw2pxQHLoUDSnhNKp1CAY0LskT/hKRdZ1G5dSpU2jcuDHc3NzUmiYpdimXEiVK3Ff4kogoMwE+blj0an28WLMo0jTgszUnMGz+AcQnpejdNSKinFvT1LdvX1WradWqVShSpAj/MiSix+Lq5KAWg1ct6o0PVx7DqkNXcSYyDj/2rIniBUw3mBARmYNsT8/JbrnQ0FCEhITkXq8sEKfniJ50ndN+3IxLhJero6og3jSksN7dIiIbEJOb03MVKlTAzZs3n6R/REQmapfwxaphjVC9mA9i/q3nNOmv0+qULERE5sI+qynMcPn888/x7rvvYvPmzbh165bJY3IhInoc/t6uWDCwHl6uVwwy/v3NX6cwYPY+RN9L1rtrRERZn56zt7c3Wbsk/0vGtUyG+2z1FCucniPKOYv3XcSYZUdUUcziBdwx7eWaKF+Ev1dEZAHFLTdt2pRTfSMieqTOtYIQ4u+lyhJcuBWvTvj7WYfK6FijqN5dIyIbluWF4B999BHeeecduLuzem9mONJElPNu303CGwvDsPXUDXVbpu7GPlcBLo4OeneNiKxErhS3dHBwwNWrV1G4MHe0ZIahiSh3pKZpmLThNCZvOK1uVwvywdQeNVStJyIis9w9x7OtEJEeHOzt8NazZTGjT214uzkh7OIdtJ28zTj6RESUV7JVcoCFLIlIL1K3ScoSyEl/b8cno/eMPfj2r1NqJIqIKC9keXpOdtDJ8NWjglNUVBRsEafniPJGQnKqqiA+f0+Eut24TEFVDNPXw1nvrhGRBcrx3XMGH374oXrivHT58mWMHDkSa9euRXx8PEqXLo0ZM2agVq1a6nHJfOPGjcNPP/2EO3fuoGHDhvjhhx9QpkwZkyA3bNgwrFy5UoW/Tp06YdKkSfD09DS2OXToEIYMGYK9e/eiUKFCqr3UoyIi8zv9yviOlVGreH6MWXYY207fVNN133evgZrF8+vdPSKyYtkKTV27ds3TheC3b99WIahp06YqNEmYOX36NPLn/+8fxokTJ2Ly5MmYNWuWOmHw2LFj0bJlSxw7dgyurq6qTY8ePdQi9vXr1yM5OVmdP2/gwIGYN2+eMWW2aNECzZs3x7Rp03D48GH069cPPj4+qh0RmZ9ONYuiYqAXBs/Zj3M37+KlH3diVOsQ9G8UzKUERJQrzHr33KhRo7B9+3Zs27Yt08el6wEBAXj77bdVOQQhw2t+fn6YOXOmCnnHjx9Xp36RESTD6NS6devQpk0bXLp0Sf3/MjI1ZswYXLt2Dc7OzsbXXrZsGU6cOJGlvnJ6jkgfsQnJGL3ksDrhr2hRwQ9fdK6qFo0TEdnM7rkVK1aooNO5c2cV1qpXr66m4QzCw8NV0JERIgN543Xr1sXOnTvVbbmWESNDYBLSXqbpdu/ebWzTpEkTY2ASMlp18uRJNdpFROYrn6sTvutWHR+/UBHODvb437HreO67bTh8KVrvrhGRlclyaEpLS8vzGk3nzp0zrk/6888/MWjQILz++utqKk5IYBIyspSe3DY8JtcZ++3o6AhfX1+TNpk9R/rXyCgxMZHn3SMyEzId17N+Cfw+qD6K5nfDxah76PTDDszacZ7lUohIn5IDeU2CWo0aNfDZZ5+pUSZZXzRgwAC17khv48ePV6NahktQUJDeXSKyeVWK+mD1643VFF1SahrGrTiKQXP286S/RGT9oalIkSJqPVJ65cuXR0TEP1uN/f391fX169dN2shtw2NyHRkZafJ4SkqK2lGXvk1mz5H+NTIaPXq0mv80XC5evPiE75aIcoKsZfqxZ02Ma1cBTg52WHf0mpquO3Tpjt5dIyILZ9ahSXbOybqi9E6dOoXixYurr2W3nISaDRs2GB+XaTJZq1S/fn11W66lFEFoaKixzcaNG9Uolqx9MrTZunWr2llnIDvtypUrZ7JTLz0XFxe1YCz9hYjMZ7qub8Ng/P5aAwT5/jdd98vf4ZyuIyLrDE1vvvkmdu3apabnzpw5o0oETJ8+XdVTMvzDOHz4cHzyySdq0biUCujVq5faEde+fXvjyFSrVq3UtN6ePXvUbryhQ4eqnXXSTnTv3l0tAu/fvz+OHj2KhQsXqjpOb731lq7vn4ieTNUgH6wa1hitKvojOVXDx6uOYcDsfepEwERE2aaZuZUrV2qVKlXSXFxctJCQEG369Okmj6elpWljx47V/Pz8VJtmzZppJ0+eNGlz69YtrVu3bpqnp6fm5eWl9e3bV4uNjTVpc/DgQa1Ro0bqOQIDA7UJEyZkq5/R0dHy56u6JiLzIv9OzN4RrpX5vzVa8ZGrtHqf/aXtCb+ld7eIyAxk5/M7y3Wa6OFYp4nI/B29Eo1h8w6oYpj2dsCbzcticNPS6qTARGSbYnKjThMRkaWrGOCNFcMaoUP1QMh5fr9afwov/7wb16IT9O4aEVkAhiYisimeLo74uktVfNm5KtydHbDz3C20nrQVG46b7qAlIsqIoYmIbI5sInmxZlGsGtYIFQO8cDs+Gf1n7cMHK44iMSVV7+4RkZliaCIim1WykCeWDG6gTvIrZu44j/ZTduBMZKzeXSMiM8TQREQ2zcXRAWOfq4AZfWqjgIczjl+NwXPf/Y15uyNY04mITDA0EREBaBpSGGvfaIzGZQoiITkN/7f0MF6bE8qaTkRkxNBERPSvwl6umNW3Dt5rW16dguXPo9fRetI27DhzU++uEZEZYGgiIkrH3t4OrzQuiaWDG6JkQQ9ci0lAj192Y/za40hKSdO7e0SkI4YmIqJMVAr0xqrXG6FbnWKQpU0/bjmHjj9sx5nIOL27RkQ6YWgiInoAd2dHjO9YGT/2rIn87k44clkWiW/Db7sucJE4kQ1iaCIieoSWFf2xbngT4yLxscuOqLpON2IT9e4aEeUhhiYioizw+3eRuJQncHa0x8YTkWj17VasP8ZK4kS2gqGJiCgbi8SlEObKoY0Q4p8Pt+4mYcDsfRj1xyHcTUzRu3tElMsYmoiIsqmcfz4sH9oQrzYpCTs7YMHei6o0QeiFKL27RkS5iKGJiOgxK4mPblMe816ph0AfN0RExaPztJ344s8TLE1AZKUYmoiInkD9UgWwdnhjdKwRiDQNmLLpLDpM3Y5T13n+OiJrw9BERPSEvFyd8HWXavihRw1VmuDolRg8N/lv/LjlLFIlSRGRVWBoIiLKIa0rF8GfbzZBs5DCSEpNw/i1J9B1+k5cuHVX764RUQ5gaCIiykGF87ni59618HmnyvBwdsDe87fVIvE5LIhJZPEYmoiIcpidnR1eql1MFcSsG+yL+KRUvLfsCHr9ugdX7tzTu3tE9JgYmoiIckmQrzvmD6inCmK6ONpj2+mbaPnNVizae5GjTkQWiKGJiCgPCmKueaMxqhfzQWxiCt7945A6Dcv1mAS9u0dE2cDQRESUB0oV8sTvrzXAqNYhcHb45zQsz369BX+EXuKoE5GFYGgiIsojDvZ2eO2pUlj1eiNUKeqNmIQUvL34IF7hqBORRWBoIiLKY2X98mHJoAYY0bKcGnXawFEnIovA0EREpANHB3sMaVoaK4eZjjrJWqer0dxhR2SOGJqIiHQ++W/6USdZ69Ti661YsCeCo05EZoahiYjITEadVr/eCNWC/tlhN2rJYfT8ZQ8uRsXr3T0i+hdDExGRmSjjlw9/DGqA99qWV3Wd/j5zEy2/3YoZ28N5DjsiM8DQRERkZjvsXmlcUlUTr/NvNfEPVx5D52k7cCYyVu/uEdk0hiYiIjMUXNADCwbUwyftK8HTxRH7I+6gzaS/8d2G00hKSdO7e0Q2iaGJiMiMq4m/XK84/vdmEzQtVwhJqWn4av0pPP/93wi7eEfv7hHZHIYmIiIzF+Djhl/71Ma3L1VDfncnnLgWiw5Tt+OjlcdwNzFF7+4R2QyGJiIiC2BnZ4f21QPx11tPoUP1QEg1gl+3h6PFN1ux+WSk3t0jsgkMTUREFqSApwu+eakaZvWrg0AfN1y+cw99ZuzF6/MP4EZsot7dI7JqDE1ERBboqbKF1Fqn/o2CYW8HrDh4Bc2/3oKFe1kUkyi3MDQREVkoDxdHjH2uApYPaYSKAV6IvpeMkX8cRtfpu3AmMk7v7hFZHYYmIiILV7moN5YPaYgxbcrDzckBu8Oj0GbSNny9/hQSklP17h6R1WBoIiKyklOxDGhS0qQ8weQNp9F60jbsOHNT7+4RWQWGJiIiKxLk667KE0ztUQOF87kg/OZddP95N95cGIabcVwoTvQkGJqIiKywPEGbykXw19tPoVf94rCzA5YeuIxnvtyMubsvII3nsSN6LHYat1nkiJiYGHh7eyM6OhpeXl56d4eIyEiqh49ZehhHr8So29WCfNTpWSoFeuvdNSKL+vzmSBMRkZWTkCQLxce1q6DOYychSk7F8sGKo4hJSNa7e0QWg6GJiMhGFor3bRiMDW8/hbZVikBm6GbuOI9mX23B8rDLrO1ElAUMTURENsTPyxVTutfA7H51EFzQQ1URf2NBGLr/tBtnImP17h6RWWNoIiKyQU3KFsK64Y3x9rNl4eJoj53nbqHVt9swfu1xngSY6AEYmoiIbJSLowOGNSujTgLcLKQwUtI0/LjlnJqyW3nwCqfsiDLg7rkcwt1zRGTpNhy/jg9XHkNEVLy63aBUAXz4fEWU8cund9eIzOLzm6EphzA0EZE1kNOuyGjT1M1nkJiSBkd7O/RuUAJvNC8DL1cnvbtHlOMYmnTA0ERE1uRiVDw+WnUM649dV7cLerpgVOsQdKweCHt7O727R5RjGJp0wNBERNZoy6kb+HDFUZy7eVfdrl7MBx+0q4iqQT56d40oRzA06YChiYisVVJKGmZsD1cnAL6blKpOy9K5ZlGMaBmCQvlc9O4e0RNhaNIBQxMRWbvrMQn4fO0JLDlwWd2W6uKvNyuNPg2C4ezIzdhkmRiadMDQRES2IvTCbXy48igOXYpWt0sW9MCYtuXxTEhhdbJgIkvC0KQDhiYisiVpaRp+D72EiX+ewM24JHVf4zIF8f5zFViigCwKQ5MOGJqIyBbFJiTj+01nMOPv80hKTYODvR1erlsMw5uXRX4PZ727R/RIDE06YGgiIlt24dZdfLr6OP73b4kCbzcnvN6sDHrWK871TmTWGJp0wNBERATsOHNT1Xc6ce2fk//KSYH/r015NC/P9U5knhiadMDQRET0j9Q0DYv2XcRX/ztpXO8kp2SRxeIVA7z17h6RCYYmHTA0ERHdv95p6uaz+OXvcFXrSQaaOtUoindalIO/t6ve3SPK9ue3RU00T5gwQQ3vDh8+3HhfQkIChgwZggIFCsDT0xOdOnXC9ev/zKkbREREoG3btnB3d0fhwoUxYsQIpKSkmLTZvHkzatSoARcXF5QuXRozZ87Ms/dFRGSN8rk6YWSrEGx46ym0qxoA+RNddtw1/XIzvl5/CncTTf8dJjJ3FhOa9u7dix9//BFVqlQxuf/NN9/EypUrsXjxYmzZsgVXrlxBx44djY+npqaqwJSUlIQdO3Zg1qxZKhC9//77xjbh4eGqTdOmTREWFqZC2SuvvII///wzT98jEZE1CvJ1x3fdqmPp4AaoWTw/7iWnquriT3+5GfN2RyAlNU3vLhJliUVMz8XFxalRoKlTp+KTTz5BtWrV8O2336qhtEKFCmHevHl48cUXVdsTJ06gfPny2LlzJ+rVq4e1a9fiueeeU2HKz89PtZk2bRpGjhyJGzduwNnZWX29evVqHDlyxPiaXbt2xZ07d7Bu3bos9ZHTc0REjyYfOWuPXMPn607gwq14dV+Zwp7qZMAsjkl6sLrpOZl+k5Gg5s2bm9wfGhqK5ORkk/tDQkJQrFgxFZqEXFeuXNkYmETLli3VQTp69KixTcbnljaG58hMYmKieo70FyIiejgJRW0qF8H6N59ShTB93J1wOjIO/WftQ7efduHgxTt6d5HIckPTggULsH//fowfP/6+x65du6ZGinx8TM+2LQFJHjO0SR+YDI8bHntYGwlC9+7dy7Rf0h9JpoZLUFDQE75TIiLbIbWb+jUKxpYRTfHqUyXV7V3novDClO0YMm+/qvtEZG7MOjRdvHgRb7zxBubOnQtXV/PaaTF69Gg1lGe4SF+JiCh7pAjm6Nblsemdp9GxRqDaYbf60FU0+2oLxi0/gptxiXp3kcgyQpNMv0VGRqr1TI6Ojuoii70nT56svpbRIFngLWuP0pPdc/7+/upruc64m85w+1FtZG7Tzc0t077JLjt5PP2FiIgeT6CPG77uUg2rhzXGU2ULISVNw6ydF/DUxE349q9TiONOOzIDZh2amjVrhsOHD6sdbYZLrVq10KNHD+PXTk5O2LBhg/H/OXnypCoxUL9+fXVbruU5JHwZrF+/XoWcChUqGNukfw5DG8NzEBFR3qgQ4IVZ/epg7it1UTnQG3eTUvHtX6dVeJq5/Z96T0R6sYjdc+k9/fTTxt1zYtCgQVizZo0qIyBBaNiwYep+KS9gKDkg7QMCAjBx4kS1fqlnz56qpMBnn31mLDlQqVIlteC8X79+2LhxI15//XW1o04WhGcFd88REeWstDQNa45cxZd/nsT5f3faBfm64a1ny+L5qoHq5MBET8rqds89zDfffKNKCkhRyyZNmqiptiVLlhgfd3BwwKpVq9S1jBy9/PLL6NWrFz766CNjm+DgYBWQZHSpatWq+Oqrr/Dzzz9nOTAREVHOs7e3w3NVArD+rafwSftKKJTPBRej7uHNhQfRZtI2/HXsuiphQJRXLG6kyVxxpImIKHfFJ6Vg5o7zmLb5LGIS/lnjVKOYD95tFYJ6JQvo3T2yUDz3nA4YmoiI8kZ0fDKmbT2LGdvDkZD8zxqnxmUKqnPaVQ0yLUFD9CgMTTpgaCIiyluRMQn4buMZLNgbgeTUfz7Knq3gh7dblEWIP/8dpqxhaNIBQxMRkT4uRsWrHXZLD1xCmiZVx4F2VQIwvHkZlCzkqXf3yMwxNOmAoYmISF9nImPxzfrTWH34qrotm+s61iiKN5qVUScNJsoMQ5MOGJqIiMzD0SvR+Gb9Kfx1/J/6fI72duhcKwhDnymtimgSpcfQpAOGJiIi8xJ28Q6++t9JbDt9U912crDDS7WDMKRpaRTxZniifzA06YChiYjIPO09H6VGnnacvaVuOzvYo1udIAx6ujT8vc3rvKaU9xiadMDQRERk3naevYVv/jqFPeFR6razoz261ymG154qxfBkw2IYmvIeQxMRkfmTjzwZcZKTAO89f1vdx/Bk22IYmvIeQxMRkYWHJwd7tebptadLccG4DYlhaMp7DE1ERJZHPgIN03aG8CQLxl+sGYTBT5diqQIbEMPQlPcYmoiILJd8FO46F4XJG05j57lbxlIF7asHqt12wQU99O4i5RKGJh0wNBERWc9uOwlPhlIFUiTzuSoBqs5TWb98enePchhDkw4YmoiIrMuBiNuYsumMsUimaFnRD0OblkHlot669o1yDkOTDhiaiIist8K4hKe1R67B8InZuExBDG1aGnWCfWEnJ7sji8XQpAOGJiIi6z+33dTNZ7E87ApS5czAAGoVz4/BTUuhabnCDE8WiqFJBwxNRES24WJUPH7cehaL9l1CUkqaui/EPx8GPV0KbSsXgaODvd5dpGxgaNIBQxMRkW2JjEnAL3+HY86uC7iblKruC/J1w8AmpdC5ZlG4Ojno3UXKAoYmHTA0ERHZpuj4ZMzeeR4zdpxH1N0kdV9BT2f0aVACPeuVgLe7k95dpIdgaNIBQxMRkW27l5SKBXsj8PO2cFy+c0/d5+HsgG51iqF/42AU8WaVcXPE0KQDhiYiIhLJqWlYfegqpm05ixPXYo2FMp+vGoABTUqifBF+RpgThiYdMDQREVF68vG6+dQNTNt8FrvDo4z3NylbCK82KYkGpQpwx50ZYGjSAUMTERE9yMGLdzB96zmsPXIV/1YrQIUiXhjYpCTaVikCJ+640w1Dkw4YmoiI6FEu3Lqrdtwt2ncRCcn/lCso4u2qFo13q1sMXq5cNJ7XGJp0wNBERERZdftuEubuvoCZOy7gZlyicdF4l9pB6NcwGEG+7np30WbEMDTlPYYmIiLKroTkVKwIu4Kf/z6HU9fjjCcIblnRH680DkaNYvm57imXMTTpgKGJiIgel3wUbzt9Ez//HY6tp24Y768a5IN+DUugTWWue8otDE06YGgiIqKccPJaLH75+xyWhV0xnqbFz8sFveqXQPc6xZDfw1nvLloVhiYdMDQREVFOkrVO83ZH4LddF3Aj9p91Ty6O9uhQPRC9G5RgvaccwtCkA4YmIiLKDYkpqapY5q/bw3Hkcozx/nolfdGnQTCereAHB1kIRY+FoUkHDE1ERJSb5OM69MJtdY67dUeuIfXfgk+BPm54uV5xdK0dxKm7x8DQpAOGJiIiyitX7txT03YL9kTgdnyyceruhWoB6iTBlQK9uOsuixiadMDQREREupQsOHgFs3acx9Er/03dyXqnLrWKon21QI4+PQJDkw4YmoiISC/yUb4/4jZm7bigpu6SUv/ZdefsYK/WPPWoVwz1S/Jcd5lhaNIBQxMREZmDO/FJWB52RZ2qJf3oUzm/fGrXXfvqAXB3dtS1j+aEoUkHDE1ERGRujl6JVmULluy/jHvJqeo+L1dHvFQ7SK19KlaAp2uJYWjKewxNRERkrqLvJWPxvouYvfMCIqLi1X0yU9e0XGH0rF8cT5UpBHsbLVsQw9CU9xiaiIjI3EmZgs0nIzFr5wWT07WUKOCOHnWLo3OtovBxt62F4zEMTXmPoYmIiCzJuRtxmLMrAotDLyI2IcVYtqBdVSlbUFyd984WxDA05T2GJiIiskTxSSlYduCKqvt0/Op/C8crB3qje91ieL5qADxcrHfhOEOTDhiaiIjI8ssW3MHcXRew6tBVY9kCTxdHteOue53iqBBgfZ9vDE06YGgiIiJrEXU3CX+EXsK8PREIv3nXeH/VIB90rxOE56pYz+gTQ5MOGJqIiMjaSETYefYW5u6JwP+OXkNy6j+RwcPZAc9XC0S3OkFqGs+Si2YyNOmAoYmIiKzZzbhENfo0f08Ezt/6p2yB4ZQtL8kpW6oHWuTOO4YmHTA0ERGRLdA0DbvORanwZHLKFkd7tKrorwpnyilbLKXuE0OTDhiaiIjIFk/ZsuzAZSzYexEnrsUa7w/0cUOnmkXRuWZRBPmad9VxhiYdMDQREZGt0jQNhy9HY+Hei1hx8Iqx7pOoV9IXnWsGoVUlf7NcPM7QpAOGJiIiIiAhORV/Hr2G30Mv4e8zN2FIGe7ODmhdqQherFkUdYN9zWb6jqFJBwxNREREpi7fuYcloZfwx/5LJovHZfquQ/VAdKgRiFKFPKEnhiYdMDQRERFlTqJG6IXbKjytOngVsYkpJrWfOlYPxHNViqCApwvyGkOTDhiaiIiIHu1eUirWH7+OpfsvYevpm+okwsLR3g5PlS2EF6oH4tnyfnBzdkBeYGjSAUMTERFR9tyITVQLx5ceuIQjl/87750Uz2xZyR8vVAtEw1IF4Ohgj9zC0KQDhiYiIqLHdyYyVp04eFnYZVy6fc94f0FPZ3XaluerBaB6kE+OVx9naNIBQxMREdGTk1iy78JtrAi7gtWHr6rz4Bk0Kl0Qc16pC70+v82vYAIRERHZLDs7O9Qu4asu77eroMoWSICSMgY1ivno2jeGJiIiIjJLTg72aFqusLrIAnLDKVv0wtBEREREZs/N2QFuyJsddQ+Se8vRiYiIiKwIQxMRERGRpYem8ePHo3bt2siXLx8KFy6M9u3b4+TJkyZtEhISMGTIEBQoUACenp7o1KkTrl+/btImIiICbdu2hbu7u3qeESNGICXlv2qkYvPmzahRowZcXFxQunRpzJw5M0/eIxEREVkGsw5NW7ZsUYFo165dWL9+PZKTk9GiRQvcvXvX2ObNN9/EypUrsXjxYtX+ypUr6Nixo/Hx1NRUFZiSkpKwY8cOzJo1SwWi999/39gmPDxctWnatCnCwsIwfPhwvPLKK/jzzz/z/D0TERGRebKoOk03btxQI0USjpo0aaJqKhQqVAjz5s3Diy++qNqcOHEC5cuXx86dO1GvXj2sXbsWzz33nApTfn5+qs20adMwcuRI9XzOzs7q69WrV+PIkSPG1+ratSvu3LmDdevWZalvrNNERERkebLz+W3WI00ZyRsSvr6+6jo0NFSNPjVv3tzYJiQkBMWKFVOhSch15cqVjYFJtGzZUh2ko0ePGtukfw5DG8NzZCYxMVE9R/oLERERWS+LCU1paWlq2qxhw4aoVKmSuu/atWtqpMjHx7TYlQQkeczQJn1gMjxueOxhbSQI3bv3Xyn3jOutJJkaLkFBQTn4bomIiMjcWExokrVNMn22YMECmIPRo0erkS/D5eLFi3p3iYiIiGy9uOXQoUOxatUqbN26FUWLFjXe7+/vrxZ4y9qj9KNNsntOHjO02bNnj8nzGXbXpW+Tcced3Ja5TTc3t0z7JLvs5EJERES2waxHmmSNugSmpUuXYuPGjQgODjZ5vGbNmnBycsKGDRuM90lJAikxUL9+fXVbrg8fPozIyEhjG9mJJ4GoQoUKxjbpn8PQxvAcRERERGa9e27w4MFqZ9zy5ctRrlw54/2yhsgwAjRo0CCsWbNGlRGQIDRs2DB1v5QXMJQcqFatGgICAjBx4kS1fqlnz56qpMBnn31mLDkg66RkCrBfv34qoL3++utqR50sCM8K7p4jIiKyPNn5/Dbr0CRnOs7MjBkz0KdPH2Nxy7fffhvz589XO9ok5EydOtU49SYuXLigwpUUsPTw8EDv3r0xYcIEODr+Nzspj0nNp2PHjqkpwLFjxxpfIysYmoiIiCyP1YQmS8LQREREZN2f3xaxENwSGLIn6zURERFZDsPndlbGkBiackhsbKy6Zr0mIiIiy/wclxGnh+H0XA4W35RTtcjJhR+0FutJUrCEMakFxam/+/H4PByPz8Px+Dwcj8/D8fhY/jGSGCSBSTaM2ds/vKgAR5pyiBzo9DWkcoP8sJnjD5y54PF5OB6fh+PxeTgen4fj8bHsY/SoESaLqNNEREREZC4YmoiIiIiygKHJAsjpWsaNG8fTtjwAj8/D8fg8HI/Pw/H4PByPj20dIy4EJyIiIsoCjjQRERERZQFDExEREVEWMDQRERERZQFDExEREVEWMDSZuSlTpqBEiRJwdXVF3bp1sWfPHlib8ePHo3bt2qqaeuHChdG+fXucPHnSpE1CQgKGDBmCAgUKwNPTE506dcL169dN2kRERKBt27Zwd3dXzzNixAikpKSYtNm8eTNq1KihdnGULl0aM2fOhKWZMGGCqjo/fPhw4308PsDly5fx8ssvq2Pg5uaGypUrY9++fcbHZc/L+++/jyJFiqjHmzdvjtOnT5s8R1RUFHr06KEK8Pn4+KB///6Ii4szaXPo0CE0btxY/U5KleOJEyfC3KWmpmLs2LEIDg5W771UqVL4+OOPTc61ZUvHZ+vWrWjXrp2qAC2/S8uWLTN5PC+PxeLFixESEqLayM/smjVrYM7HJzk5GSNHjlR99fDwUG169eqlzohhE8dHds+ReVqwYIHm7Oys/frrr9rRo0e1AQMGaD4+Ptr169c1a9KyZUttxowZ2pEjR7SwsDCtTZs2WrFixbS4uDhjm9dee00LCgrSNmzYoO3bt0+rV6+e1qBBA+PjKSkpWqVKlbTmzZtrBw4c0NasWaMVLFhQGz16tLHNuXPnNHd3d+2tt97Sjh07pn333Xeag4ODtm7dOs1S7NmzRytRooRWpUoV7Y033jDeb+vHJyoqSitevLjWp08fbffu3eq9/Pnnn9qZM2eMbSZMmKB5e3try5Yt0w4ePKg9//zzWnBwsHbv3j1jm1atWmlVq1bVdu3apW3btk0rXbq01q1bN+Pj0dHRmp+fn9ajRw/18zp//nzNzc1N+/HHHzVz9umnn2oFChTQVq1apYWHh2uLFy/WPD09tUmTJtnk8ZGf/zFjxmhLliyR1KgtXbrU5PG8Ohbbt29Xv2MTJ05Uv3Pvvfee5uTkpB0+fFgz1+Nz584d9e/IwoULtRMnTmg7d+7U6tSpo9WsWdPkOaz1+DA0mTH5QRwyZIjxdmpqqhYQEKCNHz9es2aRkZHqF3XLli3GX1L5RZF/6A2OHz+u2sgvrOGX3N7eXrt27ZqxzQ8//KB5eXlpiYmJ6va7776rVaxY0eS1XnrpJRXaLEFsbKxWpkwZbf369dpTTz1lDE08Ppo2cuRIrVGjRg98PC0tTfP399e++OIL431y3FxcXNQ/1kL+UZZjtnfvXmObtWvXanZ2dtrly5fV7alTp2r58+c3HjPDa5crV04zZ23bttX69etncl/Hjh3VB5atH5+MoSAvj0WXLl3U9ya9unXraq+++qpmLjILlZn9MSftLly4YPXHh9NzZiopKQmhoaFqWDj9+e3k9s6dO2HNoqOj1bWvr6+6luMgQ8Lpj4UM1xYrVsx4LORahm79/PyMbVq2bKlOFHn06FFjm/TPYWhjKcdTpt9kei3je+DxAVasWIFatWqhc+fOauqxevXq+Omnn4yPh4eH49q1aybvT841JVPe6Y+RTCPI8xhIe/m92717t7FNkyZN4OzsbHKMZDr59u3bMFcNGjTAhg0bcOrUKXX74MGD+Pvvv9G6dWt129aPT3p5eSws+Xcu47/ZMo0nx8Tajw9Dk5m6efOmWoeQ/kNOyG35hbZWaWlpaq1Ow4YNUalSJXWfvF/5xTL8QmZ2LOQ6s2NleOxhbSQ43Lt3D+ZswYIF2L9/v1r/lRGPD3Du3Dn88MMPKFOmDP78808MGjQIr7/+OmbNmmXyHh/2+yTXErjSc3R0VOE9O8fRHI0aNQpdu3ZVYdrJyUmFSvk9kzUnwtaPT3p5eSwe1MZSjpVhPaWscerWrZvxZLzWfHwcdXtlogeMphw5ckT9FUz/uHjxIt544w2sX79eLYakzMO2/FX72WefqdsSCuTnaNq0aejduzds3aJFizB37lzMmzcPFStWRFhYmApNsoiXx4ceV3JyMrp06aIWzssfLbaAI01mqmDBgnBwcLhvB5Tc9vf3hzUaOnQoVq1ahU2bNqFo0aLG++X9ynTlnTt3Hngs5DqzY2V47GFt5K8j2SFjrmT6LTIyUu1qk7/W5LJlyxZMnjxZfS1/edny8RGyy6lChQom95UvX17tGEz/Hh/2+yTXcpzTk92FsgsoO8fRHMlOScNok0zT9uzZE2+++aZx5NLWj096eXksHtTGEo5V8r+B6cKFC+oPOsMok7UfH4YmMyXTLTVr1lTrENL/NS2369evD2sif6VIYFq6dCk2btyotkWnJ8dBphTSHwuZ95YPRMOxkOvDhw+b/KIafpENH6bSJv1zGNqY+/Fs1qyZem8yOmC4yKiKTK0Yvrbl4yNkOjdjmQpZv1O8eHH1tfxMyT+06d+fTDvK+or0x0iCp4RUA/l5lN87Wc9iaCPbseUDI/0xKleuHPLnzw9zFR8fr9aTpCd/lMl7E7Z+fNLLy2Nhqb9zyf8GJinD8Ndff6kyH+lZ9fHRbQk6ZankgOzYmDlzptqNMHDgQFVyIP0OKGswaNAgtb138+bN2tWrV42X+Ph4ky31UoZg48aNakt9/fr11SXjlvoWLVqosgWyTb5QoUKZbqkfMWKE2l02ZcoUi9lSn1H63XPC1o+P7N5xdHRUW+tPnz6tzZ07V72XOXPmmGwjl9+f5cuXa4cOHdJeeOGFTLeRV69eXZUt+Pvvv9VuxfTbpGUXlWyT7tmzp9omLb+j8jrmtqU+o969e2uBgYHGkgOylVxKTsiOSVs8PrITVUpvyEU+Br/++mv1tWH3V14dC9lSLz+3X375pfqdGzdunO5b6h91fJKSklQJhqJFi6p/S9L/m51+J5y1Hh+GJjMntXLkw1DqNUkJAql5YW3klzKzi9RuMpB/rAYPHqy2qMovVocOHdQvaXrnz5/XWrdurWp9yAfC22+/rSUnJ5u02bRpk1atWjV1PEuWLGnyGpYcmnh8NG3lypUqGMofGiEhIdr06dNNHpet5GPHjlX/UEubZs2aaSdPnjRpc+vWLfUPu9QwknIMffv2VR8g6UndHilvIM8hQUQ+YM1dTEyM+nmRf0tcXV3V91bq8KT/kLOl4yM/55n9myPhMq+PxaJFi7SyZcuq3zkp+bF69WrNnI9PeHj4A//Nlv/P2o+PnfxHv3EuIiIiIsvANU1EREREWcDQRERERJQFDE1EREREWcDQRERERJQFDE1EREREWcDQRERERJQFDE1EREREWcDQREQWY/PmzbCzs7vvPHsZlShRAt9++22u9OHpp59WJ7t9EvIe5OLj4/PE/ZH3ani+Rx0XInoyDE1EZFZ27typzovWtm3bR7adOXNmpsFj7969GDhwoPG2BIply5bBnMyYMUOdH+9JyXv9448/cqRPRPRwDE1EZFZ++eUXDBs2TJ3M88qVK4/1HIUKFYK7uzvMmYS9woULP/HzyHv19fXNkT4R0cMxNBGR2YiLi8PChQsxaNAgNdIkI0kPm6rr27cvoqOjjdNTH3zwwX3Tc/K16NChg2pjuN2nTx+0b9/e5Dll2k2m3wzu3r2LXr16wdPTE0WKFMFXX311Xz8SExPxzjvvIDAwEB4eHuos7tK37JK+V6tWDb/++iuKFSumXnPw4MFITU3FxIkT4e/vr0LWp59+mu3nJqKcwdBERGZj0aJFCAkJQbly5fDyyy+rAPGg02M2aNBABSMvLy9cvXpVXSS8ZDZ9ZZgOkzaG21kxYsQIbNmyBcuXL8f//vc/FYb2799v0mbo0KFqSnHBggU4dOgQOnfujFatWuH06dPZfv9nz57F2rVrsW7dOsyfP1+Nukl4vHTpkurH559/jvfeew+7d+/O9nMT0ZNzzIHnICLKERISJCwJCR4yiiRhIf3oj4GzszO8vb3V6JGMwjxs+sowHfawdpmNekl/5syZg2bNmqn7Zs2ahaJFixrbREREqDAm1wEBAeo+CW4SeuT+zz77LBvvHkhLS1NBMV++fKhQoQKaNm2KkydPYs2aNbC3t1dhUoLTpk2b1IgWEeUthiYiMgsSDvbs2YOlS5eq246OjnjppZdUcMksNOU2GfVJSkoyCSeydkiCi8Hhw4fV9FnZsmXvm7IrUKBAtl9Tpg4lMBn4+fmpRfESmNLfFxkZ+RjviIieFEMTEZkFCUcpKSnGERshU3MuLi74/vvv1ahSTpIgknHqLzk5OVvPIaNREmpCQ0PVdXqyJim7nJycTG7LKFpm98mIFBHlPa5pIiLdSViaPXu2WmgdFhZmvBw8eFCFKFnfkxmZopORnkeR4JGxnUzbyRqn9OQ1DUqVKqX+v/Trh27fvm1SJqB69erqeWXkp3Tp0iaX7EwFEpFlYGgiIt2tWrVKBZL+/fujUqVKJpdOnTqpUagHTWfJaM+GDRtw8+ZNxMfHP7CdtLl27Zp6HfHMM89g3759KqzJou1x48bhyJEjJiNF0h9ZDL5x40b1mOy4Sz9VJtNyPXr0UDvslixZgvDwcDXFOH78eKxevTrHjxMR6YuhiYh0J6GoefPmmU7BSWiScCM70zLbQffaa6+ptU8yciRb8zMjI1jr169HUFCQGh0SLVu2xNixY/Huu++idu3aiI2NVeEnvS+++AKNGzdGu3btVP8aNWqEmjVrmrSRBd/y/7399ttqvZOUMZAdelI2gIisi532oP28RESUK2Rdkix4z1gn6nFJKQTZaSejaDlxahYiyhxDExGRDqHJ1dVV7bCTGkxPomLFijh37hwSEhIYmohyGXfPERHlMUPhy4w77h6H1HAy7PqTQp9ElHs40kRERESUBVwITkRERJQFDE1EREREWcDQRERERJQFDE1EREREWcDQRERERJQFDE1EREREWcDQRERERJQFDE1EREREWcDQRERERIRH+3/ABQQGwzjjkwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Compute thrust evolution with altitude\n", "turbojet.M0 = cruise_mach\n", "altitude_array = np.linspace(0, max_altitude, nb_points)\n", "thrust_array = np.zeros(nb_points)\n", "for i, altitude in enumerate(altitude_array):\n", " turbojet.altitude = altitude\n", " turbojet.tune_current_OPR()\n", " thrust_array[i] = turbojet.thrust\n", "\n", "# Plot thrust evolution with altitude\n", "plt.plot(altitude_array, thrust_array)\n", "plt.xlabel(\"Altitude [m]\")\n", "plt.ylabel(\"Thrust [N]\")\n", "plt.title(\"Thrust evolution with altitude\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "venv", "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.10.4" } }, "nbformat": 4, "nbformat_minor": 2 }