{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from lbmpy.session import *\n", "from lbmpy.parameterization import ScalingWidget\n", "from lbmpy.parameterization import Scaling" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# Tutorial 05: Nondimensionalization and Scaling\n", "\n", "## Short version\n", "\n", "Play around with the widget below and make sure that.\n", "If you can't see the widget start this tutorial interactively in a Jupyter notebook.\n", "\n", "- max lattice velocity $u_l < 0.1$, ( smaller means smaller errors)\n", "- relaxation rate should not be too close to 0 or 2, otherwise simulation might get unstable. Keep in mind that advanced models (MRT, Cumulant, Entropic) can simulate more \"extreme\" relaxation rates." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", " \n", " " ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "8013073325cc4f7da83800d88fc944bf", "version_major": 2, "version_minor": 0 }, "text/plain": [ "VBox(children=(HBox(children=(Label(value='Scaling', layout=Layout(width='200px')), Select(options=('diffusive…" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "p = ScalingWidget()\n", "p.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Long version\n", "\n", "There are three different scaling/setup strategies available. In all strategies you first have to first fix the parameters of your physical setup. These are in this example:\n", "\n", "- typical physical length\n", "- maximal physical velocity\n", "- kinematic viscosity\n", "\n", "Additionally a first guess for the spatial resolution $\\Delta x$ has to be made, either by directly specifying it, or by selecting how much cells should be used to resolve the typical length. \n", "\n", "In the **diffusive scaling** approach, the relaxation has to be chosen/guessed. Valid values are $\\omega \\in (0,2)$, where values too close to the interval boundaries might lead to instable simulations. From this information the time step length and the maximum velocity in lattice units is determined. Compressibility errors and Ma-number errors get smaller with smaller lattice velocity. When using diffusive scaling, lattice velocities are decreased by a factor of 2, when doubling the spatial resolution. However also four times as many time steps have to be done in this case.\n", "\n", "The **acoustic scaling** approach lets you choose a fixed value for the maximum lattice velocity (choose something smaller than 1, or better smaller than $0.1$ there!). Doubling the spatial resolution in this approach has the following effects: time step length decreases only by a factor of 2 instead of 4 as in diffusive scaling, the relaxation rate is chosen automatically and may tend to instable values, and since the lattice velocity is fixed, the compressibility error and Ma error remain constant.\n", "\n", "The third scaling approach lets you choose $\\Delta t$ freely, while relaxation rate and lattice velocity are determined automatically.|" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Example: Channel Flow\n", "\n", "To illustrate these concepts, we create a scenario which is defined in physical units.\n", "We use the `Scaling` class, that provides the basis of the widget shown above.\n", "\n", "![](../img/channel_with_dimensions.svg)\n", "\n", "\n", "| Quantity | Value |\n", "| ------------------------------- | -----------------------------|\n", "| Viscosity | $10^{-6} \\, \\frac{m^2}{s}$ |\n", "| Time to simulate | $3\\, s$ |\n", "| Max. inflow velocity | $2 \\frac{cm}{s}$ |\n", "| Resolution of cylinder diameter | $30$ cells |\n", "\n", "\n", "The typical length of the scenario is chosen to be the cylinder diameter, that is resolved with 30 cells." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "DiffusiveScalingResult(dt=0.0002436647173489281, lattice_velocity=0.02923976608187137)\n" ] } ], "source": [ "cm = 1 / 100\n", "sc = Scaling(physical_length=0.5 * cm, physical_velocity=2*cm, kinematic_viscosity=1e-6, \n", " cells_per_length=30)\n", "scaling_result = sc.diffusive_scaling(1.9)\n", "print(scaling_result)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In the diffusive scaling approach we fix the relaxation rate and compute the time step length and the lattice velocity. Here we have to check that the lattice velocity is not too high (smaller than 0.1 as a rule of thumb).\n", "From this information we can determine the number of time steps corresponding to the desired 3 seconds simulation time." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAADIAAAAOCAYAAABth09nAAAACXBIWXMAAA7EAAAOxAGVKw4bAAACMElEQVRIDa2W7VHDMAyG2x4D8LEBbMDHBrABrNBuAMev9i9sAEzAtRvACMAGdIQeG5Tn8Vm5JjRN0qI7n2RJ1ivZspP+crnsBU0mk1PkKeMM+Sf0wdEdI9/l+Tl84Rz9V9Yllv1GWbcPd90D+vesKzH0G3F1zjFrsfvj8VigF4ZJmZxBD1hYKiQHeoJfYU+E/IBwy7hCTknCjWfSUUgP+RqdG3SDPIOra4Wbfd2IjdgDAv4wBBD41YU1ZNJFcvqwxh2yYJMMGiIMsZl8UJzEfSiwt8V1SSP2IAK34Jf4fJOAO7lKJrmP3l2TbDOLcyQy6ZC35I3Yex0Cm/DphqRSgdj1O1iNiy5O52lV30FuxG5dCMnc1AB7p3rYSxc+fNG7m6k1kJ9D34WzrhG7dSHrgAGwCFsqXpPCLdss4oJhkR+F8R+EKvZOhZCPl3xG0MdqbuhMPp0Ssq31CS9erar/FvMSdpfLXsIiKft9bnIlw5oJPj65XvgpcrpLa9xaq9Zhb1UIgXxiD+HFNyWyQOeDkO5N6DKP1rLdtqY67M6tRSDb5ARenARyenrhc2yfZon856OqfhciZi12pxMhkDt9Aa9ebgEWOUlbyHsjX6XzPPEp7UxN2NUTOcoIh/BSIgRy171g78jV78Elurjw1SJ72CzUuzFCLsVFJ9XiamRNI3bfn0Yc4xfD/hXQ18Y2ecOW3n74N/PUQvAqfWE/CyWycYrWQ3bdn59G/BpxjYlfI/Yvg8gTUZYF1vwAAAAASUVORK5CYII=\n", "text/latex": [ "$\\displaystyle 12312$" ], "text/plain": [ "12312" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "total_time_steps = round(3 / scaling_result.dt)\n", "total_time_steps" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Alternatively we could have used one of the other two scaling strategies:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "AcousticScalingResult(relaxation_rate=1.9577133907595927, lattice_velocity=0.012000000000000002)\n", "FixedLatticeVelocityScalingResult(relaxation_rate=1.6949152542372883, dt=0.0008333333333333333)\n" ] } ], "source": [ "print(sc.acoustic_scaling(dt=1e-4))\n", "print(sc.fixed_lattice_velocity_scaling(0.1))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The scaling class also computes the Reynolds number and grid spacing for us." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARUAAAAVCAYAAABomXpWAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAGOElEQVR4Ae2c7XEcNwyGzx4XoDgdyB1IVgWROpDlCmx1YE1+Sf88dgexK/BHB3EHHquDqIMo6kB5n82S4q3IW5wJzuguxAyXJPYFQHC5IHbvY3F7e7tIy/n5+W7a7+3l+enz0eejr4G7NZCLF48XCV1cXLxRdy9h9WafgT4D/8MZUCzYVXmtsjPjPjjiRqRHRF1IJ45VHag+GxiTg/gEmy8q+2rfTE4PXfHfjfy/VT9TeSfe1ciLlRUXBSaNWvmJupXdWltWeW9ccEp6Z68b2MR+EP0k3mXopLUV641jDN46rfpa2G6h08sf6SEecL+X6EaYXzipmqBC/wP9JxzUIRr9rnqffqCR/1H9a5XnKrvh3LQW9od4b1V/5dwo+0P1kUoMLGqbcFP9oV8rH/RY6lpbVvkGOK6n9bpxTVk8ZxrHN+ZFNX0KG0Mk8U1YbxwD8NZp1dfCdgudDfw50DhZD/HeZdwjHaqOyYdsv1fhXv+schMef8gw/ggSoQag8kLlVLxPgT+tdf61eDuqh4DCebXJZuhHvVYc8jmqlc/pLPFqbVnlvXH4I52m6zb6TvAgKxkCysgjKOUWkxXrjWNY3jqt+lrYbqHT3R+tCRKC07Ro4MSKb+LFex1nRNznw5PKkKmoc4IgZ36SXkgulyp/F/+NdBNwCDJWXGkYtfIlvTl+rS2rvDcu50uWp2tCisvj0W8pQPyjtE/bivXGbZvtTfFH4+TezRGvNFizSyTeB5V/VM4e68DCyu1KS0IzHdIhHpGmFPRyHrLi/kPfP9bK39dY5tTassp748oe3T/DRkJWQ8CfIyvWG8e4vHVa9bWw3UKnuz9aE9NMZCEemchbHCgQ9/vJEx3YldLUt4DPs2WIVHmOnlpxJUW18iW9OX6tLau8Ny7nywyP92RXGgfZykuV8IL9i3jTNWHFeuNwwVunVV8L2y10tvCHcUbSeuB92p7q+C4lnrxrsGaOCCoMKL73uDtvbj0dkat2OwKPFVcyXCtf0pvj19qyynvjcr6s4oUN4Xm6WNQmjX2lku5WVqw3jvF767Tqa2G7hc4W/jDOlMhShncmKXPS/kv9Q17UMqDco8sEX9X91ShtxZXU1cqX9Ob4tbas8t64wRcFjLAQ2X2GjwITJz+r/TFgQi3eSqw3jvF467Tqa2G7hc4W/jDOlGSDLOVQ9TR7TWG0iSO7BBV2yxuVn6VVASnsxKTVVlxpHLXyJb05fq0tq7w3LufLHC+890pxfOxP0CGLTcmK9cYxBm+dVn0tbLfQ2cIfxgnxzianfziZHMDsEFSqSNErBKSw86X6Ao/ndhMuFU7btfKprrl2rS2rvDduzq/0fGI7XJf0dGizQy2sWG/cttneFH8Y54SO1V+1AQb4kKAQVACHmz+cXLcmLRoW4EQwZCohbbLiJmpit1Y+KjI0am1Z5b1xBtciBNurrn26O1mx3jgG663Tqq+F7RY6W/jDONlQWB/c26s2H6AQ2GuCCgsnFxDENhNfvJmmygjvq1wmO5gVh2yOzPLjZOR0WHlmWwWFVnlvXGE4WTYv6HPXnuvGR80s1kBWrDcO+946rfpa2G6hs4U/jBMK97U1U7kiqFyqHCA9Q+GFYcg+IlyLjxd916pJkwYab+oTdV6NrMUaOL4sd6vCs32kdeQlxCcYS/JRkaGxji1ha8bqOncZ11ZdNz7d4duR8a2+2uw2S9cNneKbsN64bbO9Kf4wzpFYD5AlUxmSiEf66fKeBPhewtLvPAY1OojPTgrxJS0MEITIbv7UufipgdqcY3FinBezBCp+CwQ+0ho4Pp5iwfOSKNKa8sgVfwAZlRYaa9qqGavr3OGOxm66biOW68YYIDaNe9dtOKOD9Jqw3jjse+u06mthu4XORv6QybI5T79igAtLJPvcs6fDr5THDr/xWQoASxIb2JE/ZE7c7AS6Tn0G+gw0mgHdYwQfEo1nPP5A7D5LGcHA3fwDf+XQA8rmX8fuwcOfAb5pSxxZDEFFNx6PMfzZCtFmK0i+kM7zGNapz0CfgYYzMMYN4sfwOiRkKpjkl4e8Rd4W4l+r3m+LM92PPgMPeAaIG/FJJ/7zGwMeI85xvxkf8OXrQ+sz8IBmQLGCf337qjp+p+lfyqmDJd848coAAAAASUVORK5CYII=\n", "text/latex": [ "$\\displaystyle \\left( 100.0, \\ 0.000166666666666667\\right)$" ], "text/plain": [ "(100.0, 0.00016666666666666666)" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sc.reynolds_number, sc.dx" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "All *lbmpy* functions get their inputs in lattice coordinates. For geometry information this means, that all quantities have to be passed in cells. So physical quantities in meters have to be divided by dx first." ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAC6UAAAPhCAYAAAC7K+h5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAB7CAAAewgFu0HU+AACvV0lEQVR4nOzde7ReVWHv/d9Mdu4hhEAgISAKyFVRCEQLKK9i9aDl9cI4irWvRaHi0GqtmqJVe9S3HlRoreCgFm/QDj3g8chri9IWEC9QbGJApWCCXBQQhASIkHs2e75/8AQeN/uSTJLs7OzPZ4xn7LnWmmuu+QTCP3yzUmqtAQAAAAAAAAAAAACAFuNGegMAAAAAAAAAAAAAAIxeonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmPSO9gbGklDIpyXM7h8uTPDaC2wEAAAAAAAAAAGDsGJ9kdmd8U611/UhuBoCdiyh9+3puksUjvQkAAAAAAAAAAADGtGOS/HikNwHAzmPcSG8AAAAAAAAAAAAAAIDRy5vSt6/lmwaLFi3K3LlzR3IvAAAAAAAAAAAAjBH33XdfFixYsOlw+VBzAWBLidK3r8c2DebOnZt99tlnJPcCAAAAAAAAAADA2PTY8FMAYPONG+kNAAAAAAAAAAAAAAAweonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKCZKB0AAAAAAAAAAAAAgGaidAAAAAAAAAAAAAAAmonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKCZKB0AAAAAAAAAAAAAgGaidAAAAAAAAAAAAAAAmonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKCZKB0AAAAAAAAAAAAAgGaidAAAAAAAAAAAAAAAmonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKCZKB0AAAAAAAAAAAAAgGaidAAAAAAAAAAAAAAAmonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKCZKB0AAAAAAAAAAAAAgGaidAAAAAAAAAAAAAAAmonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKCZKB0AAAAAAAAAAAAAgGaidAAAAAAAAAAAAAAAmonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKCZKB0AAAAAAAAAAAAAgGaidAAAAAAAAAAAAAAAmonSAQAAAAAAAAAAAABoJkoHAAAAAAAAAAAAAKBZz0hvYKz66le/mt12222ktwEAAAAAAAAAAMAY8PDDD4/0FgDYiW3TN6WXUvYspfxBKeXjpZQrSikrSim187loM9eYWkp5XSnl70spi0spD5dSNpZSHiylXF9K+WgpZc4W7GlqKeUvOms9VEpZXUpZWkr5m1LKfs1fFgAAAAAAAAAAAABgDNrWb0q//+ncXEo5Isl1SaYPcHlWkhd2Pn9eSnlbrfXSYdY7MMl3kjy736WDO58zSilvqrVe/nT2DQAAAAAAAAAAAAAwVmzTN6X3c1eSf9/Ce2bkySD9uiQfTPL7SY5K8ook/5CkrzPvq6WUkwZbqJSyS5Jv58kg/QtJTkxybJIPJVnVWefSUsrzt3CfAAAAAAAAAAAAAABj0rZ+U/rHkyxOsrjWen8p5ZlJ7tyC+/uSfD3Jx2qttwxw/d9LKVckuSzJ+CTnl1KeXWutA8xdmOSgzvgvaq3ndF27vpTyvSTfTzI1yd8l+b+2YJ8AAAAAAAAAAAAAAGPSNn1Teq31f9RaL6+13t94/3/UWt8wSJC+ac63knyzc3hAkiP7zymlTEjy7s7hz5P8zUDPSvKlzuEJpZRjWvYMAAAAAAAAAAAAADCWbNMofTu6pmt8wADXX5Jk18744lpr3yDrXNQ1fu1W2BcAAAAAAAAAAAAAwE6tZ6Q3sJVM6ho/NsD147vG3x9inR8nWZNkapLjtsK+BnXM4g9kzrSnnn/l3Mu35WMBAAAAAAAAAADYye1y6Aeecm7jxI0jsBMAxoqd5U3pJ3SNfz7A9cO6xksHW6TW2pvkts7hoVthXwAAAAAAAAAAAAAAO7VR/6b0Usrzkryqc3hTrXWgKH2fzs/VtdaVwyx5d5IjkswupUyqta7fgr3sM8yUOZu7FgAAAAAAAAAAAADAaDCqo/RSyqQkX0wyvnPqQ4NM3aXzc9VmLLu6azw9yWZH6Xk8aAcAAAAAAAAAAAAAGDPGjfQGnqbPJTm6M7641vovg8yb3Pm5YTPW7I7Qp7RuDAAAAAAAAAAAAABgLBi1b0ovpXwwyRmdw8VJ3jnE9HWdnxM3Y+lJXeO1W7itfYe5PieP7xUAAAAAAAAAAAAAYKcwKqP0UsqZSf5n53BpklfWWlcPccujnZ/TN2P5aV3jVVuyr1rrPUNdL6VsyXIAAAAAAAAAAAAAADu8cSO9gS1VSnljkgs6h79K8vu11hXD3LYpFp9WSpk5zNxNbztfXmtd37ZLAAAAAAAAAAAAAICxYVRF6aWU/zvJP+bxfd+X5MTh3k7ecUvX+JAh1u9JckDn8Oet+wQAAAAAAAAAAAAAGCtGTZReSjkxydeT9CR5MI+/If32zbz92q7xCUPMOzrJtM74ui3eJAAAAAAAAAAAAADAGDMqovRSyrFJvpVkUpLfJnlFrfXmLVjie537kuSPSyllkHmndY0v28JtAgAAAAAAAAAAAACMOTt8lF5KeX6Sb+fxN5ivTvKqWuuSLVmj1rohyXmdw0OTvH+A5/xektM7h9+vtS5u3TMAAAAAAAAAAAAAwFjRsy0XL6Ucn+TArlN7dI0PLKWc1j2/1npRv/sPSPJvSWZ2Tn04yW9LKc8Z4rEP1FofGOD8OUnekOSgJJ8upRyY5JIka5O8JMlf5vFfj7VJ3jPE+gAAAAAAAAAAAAAAdGzTKD3JGUn+eJBrx3U+3S7qd/yiJHt2HX9mM575sSQf7X+y1vpoKeVVSb6T5NlJ3tb5dHskyZtqrT/ZjOcAAAAAAAAAAAAAAIx540Z6A9tTrfW2JEcmOSvJj5OsTLImybI8HrwfUWu9fMQ2CAAAAAAAAAAAAAAwymzTN6XXWk9LctrTuP+iPPXt6U9LrXV1kk93PgAAAAAAAAAAAAAAPA1j6k3pAAAAAAAAAAAAAABsXaJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGgmSgcAAAAAAAAAAAAAoJkoHQAAAAAAAAAAAACAZqJ0AAAAAAAAAAAAAACaidIBAAAAAAAAAAAAAGi2TaP0UsqepZQ/KKV8vJRyRSllRSmldj4XNax3UinlslLKPaWU9Z2fl5VSTtqCNXpKKW8vpfywlLK8lLK2lHJ7KeUfSimHb+meAAAAAAAAAAAAAADGsp5tvP79W2ORUsq4JBcmOb3fpXmdz2tKKV9McmattW+IdfZI8p0kx/S7tH+StyX541LKn9Zav7g19g0AAAAAAAAAAAAAsLPbpm9K7+euJP/eeO8n8mSQfmOSNyZZ0Pl5Y+f8GUn+erAFSinjk1yWJ4P0byY5KckLkrw7yQNJJiX5hy158zoAAAAAAAAAAAAAwFi2rd+U/vEki5MsrrXeX0p5ZpI7t2SBUspBSd7fOfxxkhfXWtd2jheXUv45yfeTHJ1kYSnly7XW2wZY6o+THN8ZX1BrfWfXtUWllCuSLEkyI8l5pZRDa629W7JXAAAAAAAAAAAAAICxZpu+Kb3W+j9qrZfXWu9/Gsu8J0/G8+/qCtI3PWNNknd1DnuS/Pkg62wK2x9KsnCAvd6W5OzO4YFJXtu+ZQAAAAAAAAAAAACAsWGbRulPVymlJHl153BprfVHA83rnF/WOXx1577udQ5Kcmjn8OudkH0gF3WNRekAAAAAAAAAAAAAAMPYoaP0JM9Ksndn/P1h5m66Pi/JM/tdO36AeU9Ra/1Nkls7h8dt3hYBAAAAAAAAAAAAAMauHT1KP6xrvHSYud3XD+13rWWdfUsp04aZCwAAAAAAAAAAAAAwpvWM9AaGsU/X+J5h5t7dNd53K6xTOvctG2b+E0op+wwzZc7mrgUAAAAAAAAAAAAAMBrs6FH6Ll3jVcPMXd01nr6N1hnO3cNPAQAAAAAAAAAAAADYeYwb6Q0MY3LXeMMwc9d3jadso3UAAAAAAAAAAAAAAOiyo78pfV3XeOIwcyd1jdcOs866DG6odYaz7zDX5yRZvIVrAgAAAAAAAAAAAADssHb0KP3RrvH0YeZO6xqvGmadoaL0odYZUq31nqGul1K2ZDkAAAAAAAAAAAAAgB3euJHewDC6I+99hpnb/Zbyu7fCOrXffQAAAAAAAAAAAAAA9LOjR+m3dI0PGWZu9/Wfb4V17q61rh5mLgAAAAAAAAAAAADAmLajR+l3Jrm3Mz5hmLkv7vz8dZJf9rt2bdd40HVKKXOSHNQ5vG7ztggAAAAAAAAAAAAAMHbt0FF6rbUm+Vbn8JBSygsHmtc5v+kN59/q3Ne9zq158u3pry+lTB3kkad1jS9r2jQAAAAAAAAAAAAAwBiyQ0fpHX+X5LHO+PxSypTui53j8zuHvZ35Azm383NWkk/3v1hKOSDJBzuHt0WUDgAAAAAAAAAAAAAwrJ5tuXgp5fgkB3ad2qNrfGAp5bTu+bXWi/qvUWu9tZRyTpIPJDk6yXWllE8luT3JAUnOSnJkZ/o5tdZfDLKdi5O8NclxSd5ZSpmT5AtJHk6yIMlHksxI0pfk3bXW3s3/pgAAAAAAAAAAAAAAY9M2jdKTnJHkjwe5dlzn0+2iQeZ+KMmeeTwqPzLJJQPM+VKSDw+2kVrrY6WU1yT5TpJjkpzS+XRbn+RPa61XDLYOAAAAAAAAAAAAAABPGjfSG9gctda+WuvpSV6V5FtJ7k2yofPzW0leWWs9o9baN8w6K5Icm+QdSa5N8mCSdUnuyONvTZ9fa/3iNvsiAAAAAAAAAAAAAAA7mW36pvRa62lJTtuK630nj7/p/Oms0Zvk7zsfAAAAAAAAAAAAAACehlHxpnQAAAAAAAAAAAAAAHZMonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaDaqovRSysRSyhmllH8rpdxXSllfSllVSllWSvlKKeXYzVznpFLKZaWUezpr3NM5PmlbfwcAAAAAAAAAAAAAgJ1Jz0hvYHOVUvZL8u0kh/e7NDHJQZ3PaaWU85P8Wa21DrDGuCQXJjm936V5nc9rSilfTHJmrbVvK38FAAAAAAAAAAAAAICdzqh4U3opZUJ+N0j/WZLTkvxekpcn+XiS1Z1r70py1iBLfSJPBuk3JnljkgWdnzd2zp+R5K+33u4BAAAAAAAAAAAAAHZeo+VN6a/Ok0H69UleVGt9rOv6laWUf+5cm5DkrFLKubXW3k0TSikHJXl/5/DHSV5ca13bOV7cuf/7SY5OsrCU8uVa623b7isBAAAAAAAAAAAAAIx+o+JN6UmO7Rqf3S9IT5LUWpckubxzODPJof2mvCdPRvjv6grSN92/Jo+/ZT2deX/+9LYMAAAAAAAAAAAAALDzGy1R+sSu8R1DzLt9oHtKKSWPv209SZbWWn800M2d88s6h6/u3AcAAAAAAAAAAAAAwCBGS5S+rGu8/xDzDuj8rEl+0XX+WUn27oy/P8yzNl2fl+SZm7k/AAAAAAAAAAAAAIAxabRE6f8rySOd8VmllPH9J5RSjkzyqs7h12qtj3RdPqxrvHSYZ3VfP3RLNwoAAAAAAAAAAAAAMJb0jPQGNketdUUp5f/J43H6cUkWl1L+LsmtSaZ3zr0vycQkN3TG3fbpGt8zzOPu7hrvuyX7LKXsM8yUOVuyHgAAAAAAAAAAAADAjm5UROlJUmv951LK/DwenJ+e5OJ+U+5P8pEkX6i1rul3bZeu8aphHrW6azx9C7d59/BTAAAAAAAAAAAAAAB2HuNGegObq5QyMcmbk7w6SRlgyl5J/ijJywa4NrlrvGGYR63vGk/Zkj0CAAAAAAAAAAAAAIw1oyJKL6VMS3JVkg8mmZXk00kOTTIpya5JXp7k2iRHJ/n/Sinv7bfEuq7xxGEeN6lrvHYLt7rvMJ9jtnA9AAAAAAAAAAAAAIAdWs9Ib2AzfTTJizrj02utF3dd25DkylLKNUn+PclLkpxTSrm61vrTzpxHu+ZPH+ZZ07rGq7Zkk7XWe4a6XspAL3gHAAAAAAAAAAAAABi9dvg3pZfHS+63dg5v7RekP6HW2pvkI53DcUlO67rcHYvvM8wj9+0a3735OwUAAAAAAAAAAAAAGHt2+Cg9yV5JZnXGNw4zd0nX+JCu8S2DnB9I9/WfDzMXAAAAAAAAAAAAAGBMGw1Rem/XuGeYuRMGue/OJPd2xicMs8aLOz9/neSXw20OAAAAAAAAAAAAAGAsGw1R+kNJHumMf6+UMlSY3h2c37lpUGutSb7VOTyklPLCgW7unN/0pvRvde4DAAAAAAAAAAAAAGAQO3yUXmvtS/LtzuHeST400LxSym5JPtV16vJ+U/4uyWOd8fmllCn97p+S5PzOYW9nPgAAAAAAAAAAAAAAQ9jho/SOjydZ0xl/tJTyz6WUU0opR5ZSfq+U8udJfpLksM6cq2ut/969QK311iTndA6PTnJdKeUNpZSjSylvSHJd53ySnFNr/cW2/EIAAAAAAAAAAAAAADuDnpHewOaotS4tpbw6yf9KskeSkzufgXw3yX8f5NqHkuyZ5K1JjkxyyQBzvpTkw09rwwAAAAAAAAAAAAAAY8RoeVN6aq1XJTkkyVlJvpdkeZKNSdYmuTPJ15O8JsnLaq0PD7JGX6319CSvSvKtJPcm2dD5+a0kr6y1nlFr7dumXwYAAAAAAAAAAAAAYCcxKt6Uvkmt9cEkn+58ns4630nyna2yKQAAAAAAAAAAAACAMWzUvCkdAAAAAAAAAAAAAIAdjygdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaCZKBwAAAAAAAAAAAACgmSgdAAAAAAAAAAAAAIBmonQAAAAAAAAAAAAAAJqJ0gEAAAAAAAAAAAAAaNYz0hsAAABg21u4cspIb2GbOGfm2pHeAgAAAAAAAACMed6UDgAAAAAAAAAAAABAM1E6AAAAAAAAAAAAAADNROkAAAAAAAAAAAAAADQTpQMAAAAAAAAAAAAA0EyUDgAAAAAAAAAAAABAM1E6AAAAAAAAAAAAAADNekZ6AwAAAPyuhSunjPQWRo1t8Wt1zsy1W31NAAAAAAAAANiZeVM6AAAAAAAAAAAAAADNROkAAAAAAAAAAAAAADQTpQMAAAAAAAAAAAAA0EyUDgAAAAAAAAAAAABAM1E6AAAAAAAAAAAAAADNROkAAAAAAAAAAAAAADTrGekNAAAAjEULV04Z6S0wiKH+2Zwzc+123AkAAAAAAAAAjA6j8k3ppZRnlFI+Vkr5cSlleSllXSnl7lLKD0spHy+lPGeY+08qpVxWSrmnlLK+8/OyUspJ2+s7AAAAAAAAAAAAAADsDEbdm9JLKe9KcnaSaf0u7dP5HJ9kRpL3DHDvuCQXJjm936V5nc9rSilfTHJmrbVv6+4cAAAAAAAAAAAAAGDnM6qi9FLKh5P8v53DW5N8IcniJL9NsnuSI5O8NslgQfkn8mSQfmOSTye5PckBSf6ic/8ZSZYn+cut/w0AAAAAAAAAAAAAAHYuoyZKL6WcmCeD9H9MckatdWO/aVcnObeUMnGA+w9K8v7O4Y+TvLjWurZzvLiU8s9Jvp/k6CQLSylfrrXetrW/BwAAAAAAAAAAAADAzmTcSG9gc5RSxiX5+87hT5OcPkCQ/oRa64YBTr8nT0b47+oK0jfdsybJuzqHPUn+/OnsGQAAAAAAAAAAAABgLBgVUXqSlyd5dmf8qVpr75bcXEopSV7dOVxaa/3RQPM655d1Dl/duQ8AAAAAAAAAAAAAgEH0DD9lh/DfOz9rkss3nSylzEqye5IHa60PDXH/s5Ls3Rl/f5hnfT/JwUnmJXlmkjsb9gsAAIwRC1dOGektsB21/PM+Z+ba4ScBAAAAAACwxZYsWbJvkmOTvDDJvkmmjuyOAHZoa5LcneRHSf5j/vz5d2/NxUdLlP7Czs9f1lofLaX8YZIPJnnOpgmllFuTfCHJ+bXW9f3uP6xrvHSYZ3VfPzSidAAAAAAAAAAAANhhLFmypCQ5I8mZScq4ceOmjBs3bnKScSO7M4AdWl9fX99hfX19L05SlyxZ8vkkX5o/f37dGovv8FF6KWVckkM6hytKKZ9N8u4Bph6U5Jwkry2lvKrWurLr2j5d43uGeWR39b/vFu51n2GmzNmS9QAAAAAAAAAAAICn+B9J/qCnp2ePnp6e3ZOMGzduXN+4ceP6kmyVuBJgJ1P6+vrG9fX1jUvS19vb+2Bvb+/b83hj/dGt8YAdPkpPsmue/NNLz01yTJL7kixM8p0k6zrnPpXH36h+bJIvJ3ld1xq7dI1XDfO81V3j6Vu41636GnsAAAAAAAAAAADgSUuWLDkwyR9MmDBhzvjx43ebNWvWg7vuuutvp0yZsr6UMtLbA9hh1Vqzdu3aSb/97W93feihh2aXUno2btz4B0uWLPmn+fPn3/501x8Nf1XFtK7x5CRrkryk1vrVWuvDtda1tdYfJHlpkp925r22lPKCfvdtsmGY563vGk9p3TQAAAAAAAAAAACw1b0uyfjx48fPnD179gNz5859YOrUqYJ0gGGUUjJ16tT1c+fOfWD27NkPjB8/fmaS8fndF4E3Gw1R+rp+x1+stS7rP6nWujbJh7pOvWGQNSYO87xJXeO1m7XDJ+07zOeYLVwPAAAAAAAAAAAAeNKR48eP3yVJmTVr1sMjvRmA0ajz38/S+e/pkVtjzZ6tscg29mi/438fYu7VSXrz+PfqDsC715g+zPO638y+atjddam13jPUdX8SCwAAAAAAAAAAAJ6WGaWUnvHjx/f29PT0jfRmAEajnp6evvHjx/f29vb2JNl1q6y5NRbZlmqt60spy5PM7py6e4i560opK5LM6ZqfJN2x+D7DPHLfrvGgzwIAAMaWhSunjPQWGKWG+nfnnJlb+hd0AQAAAAAAjHklSUopdaQ3AjCadf13dKu8dXvc1lhkO7i5azx+mLmbrvd2nbula3zIMPd3X//5MHMBAAAAAAAAAAAAAMa00RKl/6BrvP9gk0opM5Ls0Tn8ddelO5Pc2xmfMMyzXtx1/y83f4sAAAAAAAAAAAAAAGPPaInS/0/X+LVDzHttnnyF/A83nay11iTf6hweUkp54UA3d85velP6tzr3AQAAAAAAAAAAAAAwiFERpddaf5bkis7hG0spJ/afU0qZk+SvO4cbknyl35S/S/JYZ3x+KWVKv/unJDm/c9jbmQ8AAAAAAAAAAAAAwBBGRZTe8Z4kK/P4ni8vpZxdSnlRKeXoUso7kixOsk9n7kdqrb/uvrnWemuSczqHRye5rpTyhs79b0hyXed8kpxTa/3Ftv06AAAAAAAAAAAAADuH8847b/dSyvxSyvxly5ZNHOn9jLQFCxYcXEqZv2DBgoNHei+wPfSM9AY2V6311lLKyUm+kWSvJB/ofH5nWpJP1Fo/PcgyH0qyZ5K3JjkyySUDzPlSkg9vlU0DAAAAAAAAAAAAI+L3//b7hz60ZsOEkd7HtjZr6sSNV773hJ+P9D6AsW3UROlJUmu9tpRyeJJ3JXlNkmclmZjkviTfS3J+rfXGIe7vS3J6KeX/JHlbkmOS7JFkRR5/0/o/1Fqv2JbfAQAA2HEtXDllpLfAGDPUv3PnzFy7HXcCAAAAAACw83lozYYJD67a+aN0gB3BqIrSk6TW+mCSj3Y+rWt8J8l3ttKWAAAAAAAAAAAAAOAJixYtWjbSe4DtadxIbwAAAAAAAAAAAAAAgNFLlA4AAAAAAAAAAAAAQDNROgAAAAAAAAAAAABDWr58+fh3vOMd8571rGcdPnny5KNmzZr1vGOPPfagL3/5y7tt7hrLli2bePrpp+974IEHHj5t2rQjp0yZcuR+++33nD/8wz/cb9GiRVOGureUMr+UMv+9733v3knyL//yL7u87GUvO2DPPfc8YvLkyUftv//+hy9cuHDuI4888jtt7KWXXrrrCSeccOCmeQcccMDhH/zgB+esW7euDPasdevWla997Wu7vvnNb37Gc57znENnzJjx/J6enqNmzpz5/COOOOKQ9773vXvfd999PUPtd8GCBQeXUuYvWLDg4IF+HTZ9n/POO2/3JLnssstmvPSlLz1wjz32eN7EiROPmjdv3nPf9KY3PeP222+fMNRzYEcx5G8IAAAAAAAAAAAAAMa2G264YfJ/+2//7aDly5c/EUivX7++5/rrr9/l+uuv3+WKK6548EUvetGjQ63xuc99bvf3ve99+23YsOF3YvC77rpr0l133TXp61//+h4LFy789dlnn/2b4fbzl3/5l3M++clPzqu1PnHuzjvvnHzuuefuffXVV+/6ve9979bp06f3nX766ftedNFFe3bfe8cdd0z+5Cc/Oe+6667b5bvf/e4venqemtK+6U1v2u+b3/zm7v3P//a3vx1/0003TbvpppumfeUrX5l96aWX3vbyl7989XD7Hc473/nOeRdccMGc7nP33nvvxK997Wuzr7jiit2uuuqqZUcdddS6p/sc2JZE6QAAAAAAAAAAAAAM6KGHHhr3qle96tmbgvRXvepVD7/5zW9eMXfu3N5bbrll8vnnn7/XN77xjd2XLl066JvOL7nkkl3f/e53P7PWmqlTp/adeeaZ97/iFa94pKenp/7whz+c/tnPfnbOypUrez75yU/Omzlz5mNnnXXW8sHWuuqqq2bcdNNN057//Oevfvvb3/7AYYcdtu6BBx7oOe+88/b8wQ9+sOuNN9447cMf/vCcWbNmPXbRRRft+eIXv/i3b33rW1cccMABG371q19NOOecc+b+9Kc/nfbDH/5wxt/+7d/O/ou/+IunPKu3t7fss88+60866aSVCxYsWP2sZz1rw4QJE+odd9wx8corr5zxv//3/95j5cqVPaeeeuqBN910083z5s3rbf31veiii2bfeOON04455phVp59++vLDDjts3UMPPTT+4osv3v2yyy7b/eGHH+5561vf+syf/OQnS1ufAdtD6f5TImxbpZR9ktydJFefksyZ9tQ5r5x7+XbeFQAAsMnClUP+bXCwXZ0zc+1IbwEAAAAAABildjn0A085t/GhjVn23mWbDvettd6zXTe1lSxZsuQ7PT09h02ePHnmwQcffNtQc+f/9ZVHPLhqw4Sh5uwMdp8+ceOSD//+z7bV+meeeeY+F1544V5J8oEPfOApbzJfv359OfHEEw+87rrrZmw6t3Tp0psOPvjgDZuu77vvvs9dvnz5hKlTp/ZdeeWVS4899tjf+Z9ht95668Tjjz/+kOXLl0+YPHly3x133HHT3Llzfyf0LqXM3zR+xSte8fDll19+R/dbznt7e3P00Ucf8tOf/nTatGnT+np7e8ub3vSm5V/60pfu7l7n0UcfHXfIIYccfu+990486KCD1i5btuyW/t/55ptvnnTooYeuHzdu3IC/JosWLZrykpe85JA1a9aMe/e7333fZz/72Xv7z1mwYMHBixcvnn7MMcesWrRo0bLua8uWLZt4yCGHPHfT8amnnrriq1/96q/6P+/UU0/d79JLL90jSa699tpbjjvuOP8Tka1m2bJlB65bt25lb2/vLfPnz3/l011v4N8tAAAAAAAAAAAAAIxp69atK5dccskeSXLQQQet/cQnPvGb/nMmTZpUL7744l/29PQM+Jbkf/qnf5q56S3r73nPe+7rH6R31t7wsY997J7OM8ddcMEFuw+2p8mTJ/ddfPHFv+oO0pOkp6cnb3nLW5YnyerVq8fttttuGy+44IKn/OGLXXbZpe/1r3/9g0ly6623TnnwwQfH959z+OGHDxqkJ8mCBQvWnnrqqSuS5Iorrpg56MTNMHv27I1f/vKX7xroeR/84Aef+PW+5pprdnk6z4FtTZQOAAAAAAAAAAAAwFNce+21Ux955JHxSfLGN77xwcFC7QMOOGDj8ccf/8hA166++uoZSVJKyTvf+c4Vgz3rLW95y8PTp09/LEmuueaaGYPNO+644x7Za6+9Hhvo2vz5858I3k866aSVkyZNGjCUf/7zn79m03jZsmUTB3vWJsuXLx9/8803T/rxj388efHixZMXL148eebMmb1Jcvvtt09Zv359GW6Nwbzyla98eMqUKQPu83nPe976qVOn9iXJHXfcMan1GbA99Aw/BQAAAAAAAAAAAICx5ic/+cmUTeMXvvCFq4eaO3/+/NXf+973du1/funSpVOSZN68eev33nvv3sHunzx5cj3ssMPWLFq0aJdbb711ymDzDjzwwPWDXZs1a9YT6x900EHrBpu32267PRG1b4ru+1u0aNGUc845Z6/vfe97M1asWDFhsLX6+vqyYsWK8fPmzRv0uw3lkEMOGXSfSTJjxozeNWvWTFy1apUXUbNDE6UDAAAAAAAAAAAA8BQPPfTQE53p3LlzNw41d6+99hrw+sqVK8cnye677z5stL3nnntuTJJHHnlk0L5105vDBzJ+/PjNmtf9xvfe3t6nvOX8M5/5zB4LFy58xmOPPbZZb0BfvXp1czA+1D6TJ/e6uXuBkeJPTQAAAAAAAAAAAAAwpFKeXhP9dO/fXm688cbJm4L0WbNm9X7kIx+554c//OHPf/Ob3/xk3bp1N9Ral9Ral3zmM5/55aZ7aq0juGPYMXhTOgAAMKYsXDno3/IGO5Sh/l09Z+ba7bgTAAAAAAAAxqrddtvtibeb33vvvROOOOKI9YPNvf/++ycMdH7mzJmPJcmKFSuGbVYfeOCBCUkyY8aMYd+qvq184Qtf2P2xxx4r48ePz1VXXbXsyCOPXDfQvO63yAPelA4AAAAAAAAAAADAAJ7//Oc/8bakH/3oR9OGmnvDDTcMeP2QQw5ZmyS//vWvJ917772Dhtzr168vt9xyy9QkOeigg0bsLU1Lly6dkiQHH3zwmsGC9GTw7wtjlSgdAAAAAAAAAAAAgKc4/vjj18yYMeOxJLn00kt37+vrG3DenXfeOeHaa6+dMdC1E0888ZEkqbXmggsu2H2wZ1100UW7rVq1anySvOQlL3nkaW++UW9vb0mSNWvWDNrY/upXv5pw9dVX77r9dgU7PlE6AAAAAAAAAAAAAE8xZcqU+vrXv35F8vgbxP/qr/5qr/5zNm7cmNNOO22/jRs3loHW+KM/+qOVs2fP3pgkn/3sZ+cuWrRoSv85t91224SPfOQj+yTJ5MmT+97xjnc8uHW/yebbf//91yXJXXfdNfnKK698ytvQH3300XGvf/3rn7Vu3ToNLnTxGwIAAAAAAAAAAACAAX3yk5+8b6+99tqYJJ/4xCf2Ofnkk5/1jW98Y8a111479cILL9ztqKOOOvQHP/jBrocffviage6fPHlyPe+8835VSsmqVavGv/SlLz1k4cKFc6+88spp3/3ud6d97GMf2/MFL3jBYcuXL5+QJB/96EfvmTt3bu/2/I7d3vKWtzyYJH19fTnllFOe/YEPfGDOFVdcMf2aa66Z+qlPfWr2c5/73MMWLVq0y1FHHbVqpPYIO6Kekd4AAAAAAAAAAAAAADum3Xff/bHLL7/81pNOOumgFStWTLj88stnXX755bO655xyyikPvvjFL370z/7sz5450Bqnnnrqb5cvX/7L97///futXr163Lnnnrv3ueeeu3f3nPHjx2fhwoW/Puuss5Zvw68zrBNOOGHN+973vnv/5m/+Zu9HH310/Kc+9al5n/rUp35nzp/8yZ/c/5znPGftDTfcMH2Etgk7HG9KBwAAAAAAAAAAAGBQRx999Lr/+q//uvntb3/7b/bbb7/1EydOrDNnzux9wQte8OjnP//5O7/xjW/8crg13vWudz34s5/97L/e8pa3PLD//vuvmzJlSt/kyZP79t133/Wnnnrqiv/4j/+45eyzz/7Ndvg6wzr33HPvu+SSS2477rjjHpkxY8ZjEyZMqHvttdfGl7/85Su/+c1v/uLCCy+8Z6T3CDuaUmsd6T2MGaWUfZLcnSRXn5LMmfbUOa+ce/l23hUAAOx8Fq6cMtJbgBFzzsy1I70FAAAAAABghO1y6Aeecm7jQxuz7L3LNh3uW2sdlVHtkiVLvtPT03PY5MmTZx588MG3DTV3/l9fecSDqzZM2F57Gym7T5+4ccmHf/9nI70PYHRZtmzZgevWrVvZ29t7y/z581/5dNfr2RqbAgAAAAAAAAAAANiRzJo6ceNI72F7GCvfE9ixidIBAAAAAAAAAACAnc6V7z3h5yO9B4CxYtxIbwAAAAAAAAAAAAAAgNFLlA4AAAAAAAAAAAAAQDNROgAAAAAAAAAAAAAAzUTpAAAAAAAAAAAAAAA0E6UDAAAAAAAAAAAAANBMlA4AAAAAAAAAAAAAQDNROgAAAAAAAAAAAAAAzUTpAAAAAAAAAAAAAAA0E6UDAAAAAAAAAAAAANBMlA4AAAAAAAAAAAAAQDNROgAAAAAAAAAAAAAAzUTpAAAAAAAAAAAAAAA0E6UDAAAAAAAAAAAAANBMlA4AAAAAAAAAAAAAQDNROgAAAAAAAAAAAAAAzUTpAAAAAAAAAAAAAAA0E6UDAAAAAAAAAAAAANBMlA4AAAAAAAAAAAAAQLOekd4AAABAi4Urp4z0FmCHNNTvjXNmrt2OOwEAAAAAABhZt//lcw7tfXT5hJHex7bWs8vsjQf8z//6+bZY+5RTTnnmN7/5zd2TZOnSpTcdfPDBG4a7Z968ec+99957J+69994bfv3rX9/Ufa2UMn+w+yZNmlRnzpzZe9hhh6153ete9/CZZ5754IQJg//j2/Scga5NmDCh7rrrrr3Pfvaz15588skr//RP//TBXXbZpW+4vQPtROkAAAAAAAAAAADATqf30eUTHnvkgZ0+St9ZrF+/vtx///0T7r///l2vueaaXf/+7/9+z3/913/9xb777tu7pWtt3LixrFixYsKKFSsmXH/99TMuuOCCvS6//PJfPO95z1u/LfYOiNIBAAAAAAAAAAAA2I4OP/zwNV/5ylfu7D73yCOPjP/pT3865cILL9xz2bJlU2655Zapr3nNaw5YsmTJsqHWmj179sZvf/vbt3afW7Nmzbibb7558le+8pXZP/7xj6ffc889k04++eRnL1u27OYpU6bUbfGdYKwTpQMAAAAAAAAAAACw3UydOrXvmGOOWdf//Iknnrj6bW9724NHHHHEYXfeeefkG264YfpVV1017WUve9nqwdaaMGFCHWitE044Yc2ZZ5750LHHHnvQokWLdrn77rsn/eM//uNuZ5555kNb+/sAybiR3gAAAAAAAAAAAAAAJMn06dPrGWec8cCm4x/96EfTWtcaP3583vOe99y/6Xjx4sVTn+7+gIGJ0gEAAAAAAAAAAADYYey///4bNo3Xr19fns5az372s9d3raWbhW3Eby4AAAAAAAAAAAAAdhi//OUvJ24aP+MZz9gw1Nzh3H777ZO21lrA4ETpAAAAAAAAAAAAAOwQVq1aVb74xS/umSRTpkzpO/nkkx9pXauvry+f/exn90qSUkpe97rXrdxK2wT66RnpDQAAAAAAAAAAAAAwdqxZs2bc4sWLJ3efW7Vq1bif/OQnU7/whS/Mvv322yeXUvJXf/VX98yZM+exodbauHFj6b/W2rVrx91yyy2TL7rooj3+8z//c5ckefvb3/6b+fPnr9v63wZIROkAAAAAAAAAAAAAbEc333zz1AULFhw+2PXjjjvukbPOOus3J5988qPDrbV8+fIJQ6313Oc+d/X73//+37z5zW9e2bhdYDOMG+kNAAAAAAAAAAAAAMAm//mf/7nL+eefv+dtt9024emudfPNN0/7/Oc/P3vRokVTtsbegIGJ0gEAAAAAAAAAAADYbo455phVtdYl3Z9169bdsHTp0pvOPvvsu6ZPn/7YlVdeOfPYY4899MYbb5w81Fp77733hv5rbdiwYckdd9zxs8997nN3zpkzZ8P1118/48QTTzz43/7t36Zvr+8IY40oHf5/9u48Tquy7AP4dWZhFmbYkX1RlEVcQc19ScrsdSnsNZdM3NKs3LdKrSxzKU0xc09zRStLowxNeV1AAsFcQMAVFBCGZYSBGZjlvH84QyPMDPAIDMv3+/k8n+c+577Ofa7zzPD0R7+5BQAAAAAAAAAAAFaTJMnKcZqma3XN2tatKi8vL+3Xr9+Kyy67rOSZZ56ZlpOTk5aUlOSefvrpvdZ1rdzc3Nh2220rv/e97y0cO3bs1NatW1eXlZVln3rqqdtWVlZm1B/QNKF0AAAAAAAAAAAAAFaTn59fUzdeunTpWmVOy8vLsyIiCgsLa9ZU25g99tij4qCDDvokImLSpElFr7/+el6ma/Xq1aty6NChCyIiZs+e3eJvf/tbq0zXAhonlA4AAAAAAAAAAADAatq1a1dVN541a1bumurLy8uTJUuW5EREtG7dumpN9U3p27dvRd140qRJBZ9nrf79+69c67XXXvtcawENE0oHAAAAAAAAAAAAYDW77LJLed14woQJhWuqHzduXEF1dXVEROy4447layhvUlVVVVI3rqysTJqqXYu1GlwXWH+E0gEAAAAAAAAAAABYzWGHHbYkOzs7jYh4/PHH29XU1DRZ/4c//KF93XjIkCGLP8+9//Of/7SsG/fu3XvF51lr4sSJK9fq0aPH51oLaJhQOgAAAAAAAAAAAACr6dGjR9Xhhx++KCJiypQphT/+8Y87N1b75JNPFj/00EMdIyK6du264oQTTijN9L4jRoxoPWHChKKIiDZt2lQdfPDBSzNd66WXXiocOXJku4iI3Nzc9Mgjj/xcYXmgYTnN3QAAAAAAAAAAAAAAm6bf/va3H40dO7bVwoULc6699tpuL774YvHxxx+/YMCAActzc3PTGTNmtPjb3/7W+s9//nP76urqJCsrK26//fYPcnIaj6guW7Ysa8KECfn1z61YsSKZOXNmi5EjR7Z+9NFHO9Sdv+KKK2bl5uY2ulZlZWWy6lpVVVXJnDlzcp9++ulWf/jDHzquWLEiiYg466yzPu7WrVtVpp8F0DihdAAAAAAAAAAAAAAa1KtXr8rnnntu6tChQ7d/77338seMGdNqzJgxrRqqLS4urr7zzjvfP/LII5c0tebkyZML99prr4FN1eTk5KSXXnrprAsuuGB+U3UlJSW5a1orSZIYNmzYvJtuuml2U3VA5oTSAQAAAAAAAAAAAGjUrrvuunzq1KmT77777nZPPvlkm9dff73lwoULc6qrq5PWrVtX7bDDDhVDhgz55Nxzz53fvn376kzukZ2dHUVFRdU9e/as2H///ZecffbZ83faaaflmayVlZUVhYWF1T169Fi+5557lp155pnz99133/JM1gLWjlA6AAAAAAAAAAAAAE3Kzc2N7373uwu/+93vLsx0jTRNJ66vfmbNmvXG+loL+PyE0gEAAAAAAAAAAIAtTk5xx8rm7mFj2FqeE9i0CaUDAAAAAAAAAAAAW5w+v3zzrebuAWBrkdXcDQAAAAAAAAAAAAAAsPkSSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAABgvRk+fHj7JEkGJ0kyeNq0aS2au59MTZs2rUXdcwwfPrx9c/ezJfHZbnlymrsBAAAAAAAAAAAAgPXta098bcCiikW5zd3HhtY2v23lX4/+61sbYu2RI0cWH3nkkX3Xtv7mm2/+4JxzzlmwIXrZXF1wwQVdf/Ob33RZ9Xxubm5aVFRUXVxcXN27d++K3Xfffdlhhx22+PDDDy9rjj5Zs/nz52e/+OKLLceNG1f4yiuvtHzjjTdalpSU5EZE7LnnnmXjx4+ftr7v+dhjj7W66667Or7++ustFy1alNO2bduqXXbZZekZZ5xRcuyxxy5e3/f7PITSAQAAAAAAAAAAgC3OoopFuQsrFm7xoXQ2T5WVlcmiRYtyFi1alDNz5sy8F154ofXNN9/cZbvttqv48Y9/PPvUU09d1Nw9Nqfhw4e3P/fcc3tHREydOvWNfv36rWjmlmLXXXfdcfbs2Rtl5//q6uo44YQTej322GMd6p+fN29e7r/+9a82//rXv9o8/vjj8x966KEZ2dnZG6OlNRJKBwAAAAAAAAAAAKBJJ554Ysm55547r6mabbfdtjIi4pxzzllgx/TV3XTTTR/su+++SyMi0jRNFi1alD137tycCRMmtHz66adbT58+veC9997LP+2007YbNWrU/IcffniTCRyvb/369VuRpunE5u4jU+3bt6/aZZddlo4ePbr1hlj/nHPO6VYXSB8wYMCy88477+N+/fotnzZtWt5NN93U+a233ip89NFHO3To0KHqt7/97awN0cO6EkoHAAAAAAAAAAAAoEnbbLNN1Z577lnR3H1szvr06bOioc/w29/+dmlEzHr44Ydbf+973+tdWlqa89hjj3Vo27Zt9e233/7Rxu+Uhpx22mnztttuu+X777//0u23374yIiJJksHr+z6vv/563h133NEpImLgwIHLxo0bN7WoqCiNiDjooIOWHX/88Z/ss88+/SZPnlx4++23dzrrrLPm77TTTsvXdx/rKqu5GwAAAAAAAAAAAACArd0JJ5zwyQsvvDC1ZcuWNRERd9xxR6eXXnqpsLn74lNXXXXV3GHDhpXWBdI3lOuvv75TdXV1EhFx8803z6wLpNcpLi6uufnmm2dGRFRXVyfXXXddpw3Zz9oSSgcAAAAAAAAAAABgvRk+fHj7JEkGJ0kyeNq0aS3qzx1//PG96uZ+//vft21sjYceeqh1Xd3QoUN7N1Tz5ptv5p122mk9+vbtu2NxcfFu+fn5g7p3777zMccc0/uFF15YY5i7qqoqrr322o677LJL/6Kiot2Li4t323HHHQdceeWVncrLy5N1fvD1YOedd15++eWXr9wd/Ze//GXnpupfeumlwhNOOKFn7969dyosLNy9oKBg9969e+904okn9nz99dfzGrtu1Z9RdXV1/PrXv+6w++6792/VqtVuBQUFu/ft23fHSy+9tPOSJUuazBu/+OKLhccee2yv3r1771RQULB7Xl7eoM6dO+8ycODAASeddFLPhx56qHVNTc1nrpk2bVqLuvsPHz68fd35kSNHFidJMvjcc8/tXXeuf//+O9fV1r1GjhxZPGfOnJwWLVoMSpJk8AknnNCzqR4jIh5++OGVv1N33313o797zammpiaefvrpNhER2267bcWhhx66tKG6Qw89dGnv3r0rIiKefvrpNqt+vs1BKB0AAAAAAAAAAACAjeKOO+74sFevXssjIs4///xe7777bu6qNTNnzsz5wQ9+0DsiomvXrivuvvvumavWXHnllZ123333gb///e+3efvttwvKysqyly9fnsyaNavF448/3v7ggw8ecN5553VtrI9PPvkka++99+73wx/+sOcbb7zRcunSpVllZWXZb731VuHPf/7z7nvssUf/+fPnZ6/HR19r3/ve9xYUFxdXR0SMHj269fLly1cLyFdXV8fpp5/e/cADDxzwyCOPdJwxY0ZeeXl5VkVFRdaMGTPyHn744Y6DBg0a+Otf/7rDmu5XVlaWdcABB/S9+OKLe/3nP/9puWTJkuyKioqst99+u+D666/vtt9++/VdvHhxg5njn/3sZ9scfPDBA/74xz92mDFjRl5FRUXWihUrkrlz5+ZOmTKl8MEHH+z4rW99a/s1Bdsz0aVLl6ohQ4aURkSMHDmy3bJly5r8Q4L77ruvQ0RE69atq0888cTS9d3P+jB16tQWJSUluRERe++995KmavfZZ5+yiIh58+blTp8+vUVTtRuDUDoAAAAAAAAAAAAAG0WrVq1q/vCHP7yXnZ2dLl68OPtb3/rWtqvu8vytb31r20WLFuVkZ2fHvffe+167du0+U3DFFVd0+vnPf969qqoq6du3b/l11103869//ev0F1544a3bbrvt/d12221pmqZx8803d/nFL36xTUN9HHPMMdtOnDixKCJi5513XnrHHXe898ILL7w1YsSIdw4//PBFU6ZMKTzrrLN6b6jPoSnFxcU1gwYNKouIqKioyBozZsxqu74PGzas5z333NMpTdPYY489ym666aYPRo4cOe3//u//3rrhhhtmbL/99hXV1dXJxRdf3Ouhhx5q3dT9zjjjjN7//ve/i4cOHbpgxIgR77z44otv3X///e/utttuSyMi3njjjZY//OEPu6x63b///e+Cq666qkdNTU1069ZtxU9+8pOPnnjiieljxoyZ8tRTT0274YYbZhx55JELCwoK1nob7wMPPHDp+PHjJ19yySWz6s79+c9/fnv8+PGT678OPPDApRERp59++vyIiCVLlmQ/8MADje5+Pnv27JzRo0e3joj42te+tqCgoCBd2542ptdee62gbty/f/+Kpmr79etX3tB1zSWnuRsAAAAAAAAAAAAAYNM2b968nAkTJuQ3Nt+1a9eqbt26Va3NWgcddNCyCy+8cPb111/fbdy4ccU//elPO1111VVzIyKuvvrqbV588cVWERFnn332nC9/+ctL6187ceLE/GuvvbZbRMT5558/59e//vXsrKz/7tF8wAEHLDvjjDMWHnPMMds+8cQT7a655ppuZ5555oKOHTtW19WMGDGi9bPPPtumtpdPnnnmmXdyc/+7Yfs3v/nNTy666KIuN9xwQ6M7rW9ou+6667Lnn3++dUTEW2+9lffFL35x5efwl7/8pdWDDz7YMSLixhtvnHH++efPr3/tQQcdtOyss85acOihh+4wbty44ksuuaTnscce+0b9Z6zv1VdfbXnrrbe+f/bZZy+sO7f//vsv+8Y3vvHJrrvuOuDtt98uePjhhzvcdNNNs+qv8cgjj7StqamJgoKCmpdffvmtHj16fObn/5WvfKXsggsumL9gwYLsoqKitQqmt2rVqmbPPfesePnllyvrzg0cOLCiX79+KxqqP/rooxd37dp1xezZs1vcf//97c8888yFDdXddddd7aqqqpKIiO985zvzG6rZFHz44Ycrdzzv0aNHZVO1PXv2XDk/c+bMhn+4G5Gd0gEAAAAAAAAAAABo0kMPPdRxr732GtjY64YbbmhwR/LGXH311R/X7QZ+7bXXdnv55ZcLJk6cmP/zn/+8e0TETjvttOyGG26Yvep11113Xaeqqqpk4MCBy1YNpNfJzs6Ou+++e2aLFi3SZcuWZd1///2f2UH79ttv7xgR0aJFi/S+++6b0VBY+7rrrpuzww47lK82sZG0b99+ZcB70aJFn9mE+vrrr+8cEXHYYYctWjWQXqewsDC99dZbZ0ZEzJ49u8XIkSNbNXavL3/5y6X1A+l1CgoK0jPOOGNeRERpaWnOpEmTPrMb99y5c3MjInr37l2xaiB9lWepzs7Obmz6c8nOzo7jjz9+fkTEuHHjWr3zzjsNhrMfeuihDhERAwYMWLbvvvs22891TZYsWbLyF7q4uLi6qdqWLVuunC8rK9swH/A6EEoHAAAAAAAAAAAAYKPKycmJhx9++P2ioqLqysrK5KSTTtruxBNP3G758uVJfn5+zUMPPfReQ2Hxf/3rX20iIo466qhFDQXS63To0KG6LlQ+bty4orrzVVVV8e9//7s4ImK//fZb3Lt37wZ3o87Ozo7jjjtuwed6yM+h/s7i9YPKCxcuzBo/fnxxRMTQoUMXNbXGoEGDKtq0aVMVETFmzJiWjdWdcMIJjT7nF77whWV14+nTp+fVn+vcuXNlRMS7775bMHr06MKmetmQvvvd787PysqKmpqauPPOOzusOv/iiy8Wvv322wUREd/61rc22V3SIyIqKipW/qxbtGiRNlWbn5+/cr68vLzZM+HN3gAAAAAAAAAAAAAAm7bzzz9/TpqmExt73Xjjjavtar4m/fr1W3HdddfNjIh4991386dNm1YQEfHzn//8w1122WX5qvXTp09vUbdr+DXXXNMtSZLBTb0mT55cGBExb968lTuNT5kyJa8u+Dt48OClTfW3zz77NDm/IdUPordq1WplQP3ll18urKn59PDMM8/cbk2fQWlpaU5ExMcff9zgDuIRETvttFNFY3MdOnRYuQP64sWLP5M7/va3v70wJycnXbFiRfKlL32p/xe/+MXtr7/++o4TJkzIr+txY+jTp0/l/vvv/0lExIgRI9qvOn/XXXd1iPg05H3GGWestiP8piQ/P3/lB7dixYqkqdqKioqV8wUFBRvvA2+EUDoAAAAAAAAAAAAAzeLss89euOeee5bVHe+7776LL7roogZ3s54zZ05OQ+fXpLy8PLtuXFJSsnKNbbbZpsFd0ut07dq1yfkNaf78+StD5O3atVsZDJ87d26j4fKmNLWTdv1d2VeVnb3yo4vq6urPhKR33333irvuuuu9Vq1aVVdXVyejR49ufemll/bca6+9BrZv337Xr33ta9v+85//LFpt0Q3g1FNPnR8R8eGHH+Y99dRTK+9ZXl6ePPHEE+0iIr70pS+VduzYsXpj9JOp4uLi+jvkZzdVu3Tp0pXzRUVFzf5cGf3jBAAAaG6/alPe6NzFpQUbsRPYtDT1bwMAAAAAAAA2NS+++GLhq6++2rLueOrUqYWzZs3K6datW9WqtVVVVStD0eedd96cE044Ya12va4f9K0vSZrciLpZvf766yv/j++BAweu3Mm8quq/H8uvfvWrGQcddFBZrIUOHTpskNDysGHDSo888sg37r333rbPPPNM6wkTJhQtWrQop7S0NOeJJ55o98QTT7QbOnTogscee+yD+gH39e2444775MILL6xasGBBzj333NPh8MMPL4uIePDBB9ssXrw4O+K/wfVNWY8ePVbUjT/88MMm/wBh5syZK+d79uzZbH9AUUcoHQAAAAAAAAAAAICNrqysLDn55JO3raqqSgoLC2vKy8uzFi5cmHPSSSf1fu65595ZtX6bbbZZmcjOzc1N99xzz4pVa9akQ4cOa73r+OzZszPalfzzWrJkSdakSZOKIiIKCwtr9tlnn2V1c/V3+i4sLKzJ5DNY39q3b1990UUXza/b4X7SpEn5f/rTn9rcfffd25SUlOQ+/vjj7X/5y18uu+KKK+ZtqB7y8vLSb3zjGwvuuOOOTv/4xz/afvLJJzNbt25dc//993eIiOjSpcuKo48+evGGuv/6suuuu67chWzq1Kn5TdVOmzZt5R8u1L+uuTS6FT8AAAAAAAAAAAAAbChnnXVWj/fffz8/IuLXv/71jNNPP31uRMTo0aNbX3/99R1XrR8wYMDyoqKi6oiIcePGFWVyzx133HF5fn5+TUTExIkTWzZVO27cuCbnN5Tf/va37cvKyrIjIr74xS+W5ub+Nxv/hS98YVndDu9jxozJ6DPY0AYNGlTxy1/+8uOxY8e+VVBQUBMR8Ze//KXduqyRyS72Z599dklERHl5edZ9993X9t13380dO3Zsq4iIY489dsGG3Kl9fenfv/+Kjh07VkZEjBs3rrip2rp/A9tss01l3759VzRVuzEIpQMAAAAAAAAAAACwUT366KOtH3rooY4REUccccTC7373uwuHDx8+q1+/fuURET/5yU+6v/7663n1r8nJyYmDDz74k4iIMWPGtJo0aVKTO0k3JDc3N/baa68ldWvMmDGjwd3Qq6ur45FHHmm/rut/Xm+88Ube1Vdf3b3u+Ec/+tHH9ee7du1ateuuuy6NiHjyySfbzZ49O2dj97i2tt9++8revXtXREQsXLhwnfqs+8OBiIiKioq1SqjvsssuywcPHlwWEfHAAw90uPPOO9vX1NREkiRx1llnzV+X+zeXrKys+PKXv1waEfH+++/nP/vssw3+YcSzzz7bsu4POr785S+XZmU1fyS8+TsAAAAAAAAAAAAAYKsxe/bsnO9973u9IyI6d+684t57750ZEZGfn58++OCD7+Xl5aUVFRVZJ5544naVlZWfufaKK674ODs7O2pqauLYY4/t8+677zYYKo+IqKqqittuu63dqjVnnnlmSUTEihUrklNOOaVXVVXVatf+6Ec/6vz2228XfO6HXQePPPJI6wMPPLD/0qVLsyIivve97328zz77lK9ad9lll82JiCgrK8s++uij+8yfP7/RLcDLy8uTa665puOyZcvWfevxNXjggQfaNHXvd955J7cuON2jR4/l67J29+7dV/7gp06dutZ/fDBs2LD5ERETJ04suueeezpFROy1115L+vfv3+w7iUdEjBw5sjhJksFJkgw+5phjejdUc8kll8yt29X93HPP7VlWVvaZn11ZWVly7rnn9oyIyM7OTi+55JK5G7rvtbHJ/nUEAAAAAAAAAAAAAFuek046qdeCBQtysrKy4q677vqgQ4cO1XVze+yxR8WPf/zjj6688soeb775ZuHFF1/c9aabbppdN7/XXnuVX3nllR/+5Cc/6fHuu+/m77bbbgNPPPHE+UOGDFncpUuXyvLy8qx333037+WXX275j3/8o21JSUnu+PHjJ/fp02dlyPmEE0745O677/5k9OjRrUePHt168ODB/b///e/PHTBgwPKPP/4457777uvw97//ve3AgQOXTZ48uXB9Pfe7777bYsKECfkREWmaJqWlpVlz587NnTBhQuGoUaPaTJ8+fWUI/vjjj58/fPjwWQ2t881vfvOTUaNGzbv33nu3eeWVV4oGDBgw8Nvf/nbJgQceWNaxY8eqsrKyrGnTpuW/9NJLRf/85z/bLl68OPvss89eEBHp+nqWiIhbbrml03e+851tDznkkE8OPvjgJQMHDixv27Zt9YIFC3LGjx9fePfdd29TUVGRFRHxne98p2Rd1t57772X5eXlpcuXL0+uuuqqri1atEj79OmzvG5H8N69e68oKipa7XmGDRu28Ic//GGPsrKy7AULFuRERHz7299eL7ukjx07tuCVV15p8PehpKQkZ/jw4Z/ZWf/kk09e1Lp165qG6puyyy67LD/zzDM//t3vftd58uTJhXvttVf/888//+N+/fotnzZtWt5vfvObzm+99VZhRMRZZ501d+edd16nwP+GIpQOAAAAAAAAAAAAwEbx61//usNzzz3XJiLijDPOmHvEEUcsWbXmxz/+8bxRo0a1HjNmTKvf/va3XY444ohPhgwZsrRu/sorr5zXsmXLmssvv7xHWVlZ9h133NHpjjvu6NTQ/XJzc9OCgoLVwst//vOf3xsyZMgOkyZNKnr99ddbfuc739mu/vyAAQOW3X777TMOOOCAAZ/7oWudd955vddU06dPn4orrrhi1sknn1zaVN3dd9/9Ybt27apuuummLvPnz8+98cYbu954440N1hYUFNRkZ2ev10B6nYqKiqynnnqq7VNPPdW2ofmsrKy48MILZ5900kml67Ju27Zta0455ZS5t99+e+cpU6YUDh06dIf683/729+mN/S7U1RUlB599NELH3rooY61x9Xf/va3F63LvRvzpz/9qe1vfvObLg3NffDBB/nnnntu7/rnDjvssCWtW7fOaIf24cOHzyopKcn54x//2OGtt94qXPX3MyLi2GOPnX/zzTc3+IcLzSGruRsAAAAAAAAAAAAAYMv3xhtv5F155ZU9IiL69+9fftNNNzUYqM3KyooHH3zwgzZt2lRVV1fHKaecsu0nn3zymczrhRdeOH/69OlvXHjhhbMHDRpU1qZNm6rs7Oy0oKCgplevXssPO+ywRddff/2MGTNmvL7TTjuttpN027Zta8aNGzft6quvnrnTTjstKywsrGnZsmVN//79y3/4wx/OeuWVV6Z27NixasN8EhE5OTlpmzZtqrp37778wAMP/OTcc8+d889//nPaO++8M3lNgfSITz+jX//613PefPPNN88666yPBw4cuKx169bV2dnZ0bJly5o+ffpUHHXUUQtvueWWD2bNmvVaQ7uKf15//OMf3/v1r38948gjj1zYv3//8g4dOlRmZ2enhYWFNdtvv33FiSeeWDJ27Ngp119//ZxM1r/11ltn3XDDDTMGDx5cVvdsa2PYsGEL6sZHHnnkwg3x7BtadnZ2PPbYYzNGjBjxzqGHHlrasWPHytzc3LRjx46Vhx56aOmjjz769qOPPjpjbT+TjSFJ083uc95sJUnSPSI+jIh49piIzi1Xr/lql5EbuSsAANjyXFxasOYi2EL9qk15c7cAAAAAAAA0s+IBl612rnJhZUy7YFrdYY80TT/aqE2tJxMnTvxHTk7Ojvn5+W369ev3TlO1X3viawMWVSzK3Vi9NZe2+W0r/3r0X99q7j5gU3HDDTd0uOiii3pFRDz33HNvHXLIIcuau6dN0bRp07avqKgoraqqmjJ48OCvft71ctZHUwAAAAAAAAAAAACbEkFt2Do98MADHSIidthhh3KB9I0na80lAAAAAAAAAAAAAACbtqeeeqrotddeaxkRccopp5Q0dz9bEzulAwAAAAAAAAAAAACbpenTp7eoqKhI/vOf/xT88Ic/7BER0aFDh8of/OAH85u7t62JUDoAAAAAAAAAAAAAsFk65JBD+s2ePbtF/XPXX3/9h0VFRWlz9bQ1ymruBgAAAAAAAAAAAAAAPo+WLVvW7LbbbksfeeSRd0455ZRFzd3P1sZO6QAAAAAAAAAAAADAZmnWrFlvNHcP2CkdAAAAAAAAAAAAAIDPQSgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMZymrsBAACA9e1Xbcobnbu4tGAjdgIbRlO/4wAAAAAAAACwsW3WO6UnSXJdkiRpvdfBa3HN4UmS/CVJko+SJFle+/6XJEkO3/AdAwAAAAAAAAAAAABsWTbbndKTJNktIi5Yh/qsiLgzIk5bZapb7etrSZLcHRFnpmlas776BAAAAAAAAAAAAADYkm2WO6XXC5jnRMS8tbzs6vhvIP3ViDg+IvaqfX+19vzpEfGL9dcpAAAAAAAAAAAAAMCWbbMMpUfEORGxZ0RMjYh71lScJEnfiLio9vCViNgvTdMRaZpOSNN0RETsX3s+IuLiJEm23wA9AwAAAAAAAAAAAABscTa7UHqSJD0j4ue1h2dFxIq1uOy8+HRX9YiIH6RpWl5/Mk3TZRHxg9rDnIg4//N3CgAAAAAAAAAAAACw5dvsQukRcWtEFEXEH9I0fX5NxUmSJBFxdO3h1DRNxzVUV3t+Wu3h0bXXAQAAAAAAAAAAAADQhJw1l2w6kiQ5NiKOiIiFEXHRWl62bUR0rR2vKcT+fET0i4huEdE7It5f9y4BAIBN2a/alDc6d3FpwUbsBJrW1O8qAAAAAAAAAGxKNpud0pMkaRMRN9ceXpqm6fy1vHTHeuOpa6itPz9gLdcHAAAAAAAAAAAAYDMybdq0FkmSDE6SZPDw4cPbrzo/fPjw9nXz06ZNa9EcPcLmZHPaKf36iOgcEWMi4p51uK57vfFHa6j9sN64xzrcIyIikiTpvoaSzuu6JgAAAAAAAAAAALDu/vjHPw6oqKjIbe4+NrT8/PzK//3f/31rQ60/cuTI4iOPPLJv3fFXv/rVRX//+9/fa+qaY445pvfjjz/ePiIiTdOJG6q3119/Pe/WW2/tOHbs2OKZM2fmLV26NCs3Nzdt165dVc+ePZfvuuuuyw4++OAl//M//7OkuLi4ZkP1AWwmofQkSQ6IiNMjoioizkrTNF2Hy4vrjcvWULu03rhoHe5R58M1lwAAAAAAAAAAAAAbWkVFRW55efkWH0rf2J566qm248ePL9hrr73Km7OPCy+8sMvNN9/cpbq6Oql/vrq6Opk9e3aL2bNntxg3blzxHXfc0em73/3ux7/73e9mNVevsDXY5EPpSZK0iIg7IyKJiN+kafrmOi6RX2+8Yg21y+uNC9bxPgAAAAAAAAAAAABbtDRN4/LLL+/69NNPv9tcPfzoRz/qfOONN3aNiCgqKqr+1re+VXLwwQeXde/efcXy5cuzZs6cmTtu3Liip59+uvX777+fv6b1GnLOOecsOOeccxas385hy7XJh9Ij4kcR0T8iZkbEzzK4vqLeuMUaavPqjTP5C54ea5jvHBETMlgXAAAAAAAAAAAAoFm1adOmqrS0NOeZZ55pM2bMmIL99ttvo++WPmfOnJy6QHqnTp0qX3rppbe23377ylXrhg0bVhoRH40ePbqwpKRkc8jLwmYtq7kbaEqSJP0j4oe1hz9I03RpBsssqTcuWkNty3rjsnW9UZqmHzX1ioiP13VNAAAAAAAAAAAAgE3B6aefPq9FixZpRMQVV1zRrTl6eOKJJ1otX748iYg499xz5zQUSK/vkEMOWXbssccu3jjdwdZrkw6lR8T58enu5u9FRGGSJMet+oqInerVf7HeXF3A/KN6893XcL/6O51/+Lm7BwAAAAAAAAAAANhC9OjRY8Xxxx9fEhExevTo1qNHjy7MdK3Zs2fnnHPOOV0HDBiwY3Fx8W55eXmDunXrtvPXvva1bUeNGtXoJsQzZ85sUTfu27fv8kzvvybDhw9vnyTJ4CRJBk+bNq3FqvN77bVXvyRJBu+11179IiJee+21vOOPP75Xt27dds7LyxvUsWPHXb761a9u9+yzz7ZcfXXY8mzqofS82vftIuKRRl7H1Ku/ot75jrXnptSb77+G+9WffyuzlgEAAAAAAAAAAAC2TD/72c8+zs/Pr4nIfLf0xx9/vFW/fv12uuWWW7pMnTq1oKysLHvFihXJ7NmzWzzxxBPtvvKVr/T79re/3bO6unq1a+t2ao+ImDx5cn7GD7IePfbYY6322WefHUeMGNFh9uzZLVasWJHMnz8/96mnnmr75S9/uf9VV121TXP3CBvaph5KXx/ej4jZteOD1lB7YO37rIj4YEM1BAAAAAAAAAAAALA56tWrV+W3vvWtkoiIMWPGtGpqV/OGjB07tuC4447bvqysLDsnJyc97bTT5v7tb3+b/n//939v/epXv5rRrVu3FRERDzzwQMfvf//7q4Xe99hjj2V145tuuqnLyy+/XPB5n+nzmDdvXu7pp5++XXZ2dnrZZZfNeuaZZ6Y+88wzUy+99NJZRUVF1TU1NfGTn/ykxwMPPNCmOfuEDW2TDqWnaTosTdOkqVdE/KzeJYfUm/ugdo00Ip6one+fJMneDd2r9nzdTulP1F4HAAAAAAAAAAAAQD0//elPPy4oKKiJiPjJT37SdV2uPfPMM3tXVlYm2dnZ8eijj75z9913f3TEEUcsOeigg5ZddNFF81999dUpffr0qYiIuPPOOzu/8sorn9kN/fDDD1/St2/f8oiIRYsW5ey33347fuELX+h72WWXdX7yySeLFyxYkL2+nnNtzJgxIy9N03j++eenXnPNNR8PGTJk6ZAhQ5Zee+21H48ePXpqUVFRdUTERRdd1HP58uXJxuwNNqZNOpS+Ht0UEXX/DYdbkiT5zF/F1B7fUntYVVsPAAAAAAAAAAAAwCq6detWNWzYsHkREf/+97+L//a3vxWvzXWjR48ufPPNNwsjIr75zW+WDB06dPGqNR07dqz+3e9+90FERE1NTdx8883b1J/Pzs6OP//5z+/27NlzeUREmqYxfvz44uuuu67b0Ucf3bdjx4679e3bd8fvfve73V5//fW8z/moa+WCCy6YM2jQoIpVz++xxx4V55577pyIT3dUf/jhh9tsjH6gOWwVofQ0TadHxK9qD/eIiDFJknwzSZI9kiT5ZkSMqT0fEfGrNE3fbo4+AQAAAAAAAAAAADYHV1555cctW7asiYj46U9/ula7pY8aNapV3fjMM8+c31jdl7/85aXbbbddRUTEiy++2GrV+Z122mn5m2++OeXyyy//qFevXsvrz6VpGm+//XbB7bff3nnQoEE7nXHGGd0rKyvX9rHWWZIkcdZZZy1obP7ss89ekCSfbpD+zDPPrFV4HzZHW0UovdaPI+L3tePdI2JEREyofd+99vw9EXH5xm8NAAAAAAAAAAAAYPPRuXPn6tNPP31uRMSkSZOK/vznP68WHl/V5MmTCyIicnNz03322WdZU7W777770oiImTNn5lVUVCSrzhcXF9f8/Oc/n/vBBx+8+eabb77529/+9v1TTjll3s4777w0K+vTeGx1dXXcfffdnb75zW/2XvcnXDvdunVb3qVLl6rG5rt27VrVtWvXFRERU6dOLdxQfUBzy2nuBjaWNE1rIuK0JEn+HBHfiYg9I6JDRMyPT8Ppd6Rp+lQztggAADSzX7Upb3Tu4tKCjdgJW4umfucAAAAAAABgU3f55ZfP/f3vf7/NkiVLsq+66qquxxxzzOKm6ktLS3MiIlq3bl2Vm5vb5NqdOnWqjPh05/OSkpLsHj16NBr8Hjhw4PKBAwcuj4iFERHvv/9+7g9/+MOujz76aIeIiL/85S/tR40aNf+www4rW7cnXLP27ds32ledDh06VM6aNatFaWlp9vq+P2wqNvud0tM0/Wmapknt6//Wov4faZp+LU3Tbmma5tW+f00gHQAAAAAAAAAAAGDtdejQofrMM8+cGxHx+uuvt3zkkUdar811SbLaxufr1bbbbls5YsSIGUccccTCunMjRoxouyHutaGfBTYXm30oHQAAAAAAAAAAAIDm8aMf/WhumzZtqiIifvGLX3RtqraurrS0NKeysrLJdefOnZsb8Wnou2PHjtWZ9HbGGWfMrxu///77eZmssSbz58/PWYua3IiINm3aZPQcsDkQSgcAAAAAAAAAAAAgI23btq05++yz50ZETJkypfD+++9v01jtwIEDyyMiKisrk5dffrmwqXX/85//tIyI6Nmz5/L8/Pw0k9569OixMvm+oXY0nzVrVt7HH3+c3dj87Nmzc2bPnt0iIqJ///7lG6QJ2AQIpQMAAAAAAAAAAACQsUsvvXReu3btVu6WnqYNZ8gPO+ywxXXjO++8s0Nj6/3rX/9q+e677+ZHRBxwwAGL68/V1NSsdV/1g++9evVavtYXroM0TeOOO+5o39j8bbfd1r7u8xgyZMjixupgcyeUDgAAAAAAAAAAAEDGWrVqVfODH/zg44iIt99+u2D06NGtG6o75JBDlg0cOHBZRMSIESM6PPHEE8Wr1ixYsCD77LPP7hURkZWVFeecc868+vM33nhjh+OPP77Xm2++mddUT9OnT2/xi1/8olvd8de//vXSdX6wtXTjjTd2fe2111brZ9KkSfk33XRTl4iIjh07Vp5wwgkbrAdobjnN3QAAAAAAAAAAAAAAm7eLLrpo3u9+97tOJSUluaWlpY3mU++8884PDj744AGVlZXJ//7v/+4wbNiweUcffXRpcXFxzYQJEwpvuummzh999FFeRMR3vvOdj/fcc8+K+tevWLEiGTFiRIcRI0Z0GDRoUNkhhxyyeI899ljWtWvXyqysrJg5c2buc889VzxixIiOS5cuzYqIOPTQQ0uPPvroJRviuXv27Ll84cKFOQcccMCA73//+x8feuihSyIinn322eJbbrmlc1lZWXZExPXXXz8zPz+/4S3kYQsglA4AALAWftWmvMHzF5cWbORO2Nw09rsDAAAAAAAAW5KioqL0vPPOm/PjH/+4Z1N1++67b/mIESPeOfnkk7crKyvLvuuuuzrdddddnVatO+mkk0p++9vfzlr1/DbbbFOVm5ubVlZWJpMmTSqaNGlSUVP3O/LIIxc+/PDDH6zzA62lTp06VV533XUfnnLKKdtdc8013a655prPzGdlZcXll1/+0bBhw0o3VA+wKchq7gYAAAAAAAAAAAAA2Pydf/758zt37rxiTXVDhw5dPG3atDe///3vf9y/f//yoqKi6hYtWqRdunRZcdRRRy385z//Oe3++++fmZ2dvdq1p59++qI5c+a8ds8997x30kknley2225L27dvX5Wbm5vm5uambdq0qdp1112XnnbaaXNfeOGFt5588sn3i4qKNugO5ccdd9wnY8aMeesb3/jGgq5du67Izc1N27VrV3XYYYctGjVq1NSf/exnczfk/WFTYKd0AAAAAAAAAAAAYIuTn59f2dw9bAwb+jmPOOKIJWmaTlyb2oKCgnTOnDlvrE1t165dq2655ZZZEbHabuhr0r59++pTTz110amnnrpoXa+t069fvxVNPdc555yz4JxzzlmwtusNGjSo4o9//OMHmfYDmzuhdAAAAAAAAAAAAGCL87//+79vNXcPAFuLrOZuAAAAAAAAAAAAAACAzZdQOgAAAAAAAAAAAAAAGctp7gYAAAA2Z79qU57RdReXFqznTtgYMv15AwAAAAAAAMCWzE7pAAAAAAAAAAAAAABkzE7pAAAAAAAAAAAAALAOxo8fP625e4BNiZ3SAQAAAAAAAAAAAADImFA6AAAAAAAAAAAAAAAZE0oHAAAAAAAAAAAAACBjQukAAAAAAAAAAAAAAGQsp7kbAAAA2Br9qk15o3MXlxZsxE5YVVM/GwAAAAAAAABgdXZKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMZymrsBAAAAPutXbcrX+5oXlxas9zU3BRviswIAAAAAAAAA1o2d0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAbxDHHHNM7SZLB3bp127m5e9nUjRw5sjhJksFJkgweOXJkcXP3A+sip7kbAAAAAAAAAAAAAFjfHvnZvweUl63Ibe4+NrSCohaVx//kC281dx/A1k0oHQAAAAAAAAAAANjilJetyC1fUrnFh9Jp2vDhw9ufe+65vSMipk6d+ka/fv1WNHNLsEXKau4GAAAAAAAAAAAAAADYfAmlAwAAAAAAAAAAAACQsZzmbgAAAIAN71dtypu7BQAAAAAAAABgC2WndAAAAAAAAAAAAACaVFFRkVx77bUdv/CFL/Rt27btrrm5uYM6dOiw60EHHbT97373u3bV1dVrtc7777+fe/rpp3fv3bv3TgUFBbu3bdt214MPPnj7P/3pT63WdO3999/fZsiQIX06deq0S4sWLQa1bNly9+7du+88ePDgfueee27X0aNHF9bVjhw5sjhJksHnnntu77pz/fv33zlJksH1XyNHjiyuf49nn3225TnnnNN1r7326tehQ4ddc3NzBxUVFe3ep0+fgSeeeGLPiRMn5q/tZ/boo4+2Pvroo7ft3r37zgUFBbvn5eUN6tat286HHXZYn+HDh7dfsmRJxjneBx54oM3hhx++XZcuXXbOy8sbVFxcvNtOO+004MILL+xSUlKSnem6kCk7pQMAAAAAAAAAAADQqGnTprU4/PDDd3j//fc/E8hesGBBzgsvvND6hRdeaH3PPfd0/Mc//vFOp06dGk2nv/DCC4Vf//rXd1i4cOHK/GpFRUXW888/3/r5559vffrpp8+96667Plr1uqqqqjjqqKO2e+qpp9rWP19ZWZksW7asxaxZs1pMmjSp6Nlnn2395ptvvpXpcw4fPrx9/RB7vfsn7733Xv57772X/+ijj3b8xS9+MfOyyy4raWydjz/+OHvo0KF9Xn755eJV52bPnt1i9uzZLZ5++uk2ERHnnHPOgnXpsaSkJPuoo47qM27cuM+svWLFiuzJkycXTp48ufC+++7b5rHHHnvn0EMPXboua8PnIZQOAAAAAAAAAAAAQIM++eSTrCFDhvT96KOP8iIihgwZUnrKKafM79GjR+U777yTd9ttt20zYcKEokmTJhUdfvjhO4wfP35qTs7q8dTy8vKs448/vk9ZWVn22Wef/fGRRx75SV5eXs3YsWOLbr755s4lJSW5d999d6eePXuuuOKKK+bVv/b666/vWBdIHzRoUNmwYcPm9+3bd3lRUVF1SUlJzn/+85/CZ555ptWSJUtW7hB+4IEHLh0/fvzkP/3pT22uv/76bhERf/7zn9/u0aPHivpr9+vXb+VxVVVV0qpVq+ovfelLpQcccMCS/v37Ly8qKqr+8MMPW0ycOLHw7rvv3qa0tDTnRz/6Uc8dd9yx4qijjlqy6nMuWbIk68ADD+z39ttvF0REDBw4cNmwYcNKdt111/L8/Pz0gw8+aPH8888XPfnkk+3W9WdRXl6eHHzwwX2nTJlSmJ2dHUceeeSCww8//JPtt99+RWVlZTJ69OiiO+64o9PChQtzhg4dusOECROm9O3bd8WaV4bPTygdAAAAAAAAAAAAgAZdcsklXesC6eecc86cm2++eXbd3AEHHLDs5JNPXvT1r3992yeffLLdq6++2vKGG27oeOmll662i/iiRYtylixZkj755JPTDz/88LK684cccsiyE088cdHee+89YO7cubnXXHNNt9NOO21h165dq+pqHn/88XYREbvsssvScePGTcvNzf3M2l/72teW/PSnP507d+7claH0Vq1a1ey5554VL7/8cmXduYEDB1bUD6Gv6utf//onZ5xxxsLi4uKaVabKjzvuuE8uu+yyefvuu2+/6dOnF/zsZz/retRRR01bdY3zzz+/a10g/aSTTiq57777ZmZlZa2cP+CAA5addNJJpRUVFbNKSkqyV72+KZdccknXKVOmFBYXF1f//e9/n37AAQcsqz9/2GGHlZ166qkL999///4lJSW5F110Ubcnn3zy/XW5B2Qqa80lAAAAAAAAAAAAAGxtysvLk0ceeaRDRMT2229fccMNN8xetSYrKyvuvffeGW3atKmKiLjrrru2aWy9E044oaR+IL1O7969K3/xi198WHvPrNtvv719/fn58+fnRkTsueeeZasG0uvr1KlT9Vo+WoO23XbbygYC6Su1b9+++oorrpgdETFp0qSijz/++DOh8vnz52c/9NBDHSM+3SH997///WcC6fXl5+enPXr0qGpwsgGffPJJ1n333dcxIuKyyy6btWogvU7fvn1XXHjhhXMiIp566qm2ixcvlhVmo/CLBgAAAAAAAAAAAMBqXnrppcIlS5ZkR0Qcf/zx83Nychqsa9euXc3//M//LIqIePfdd/NnzJjRYHL8jDPOWNDYvU466aTS4uLi6oiI0aNHt6o/17Fjx8qIiGeeeabNnDlzGm5iA1i8eHHWtGnTWrzyyiv5EyZMyJ8wYUJ+bm5uWjc/fvz4wvr1f//734srKiqyIiLOOuuseY19Xpl46qmnisvKyrIjIr71rW8taqr20EMPXRIRUVVVlYwZM6awqVpYXzbaP0wAAAAAAAAAAAAANh+vvfZaQd143333XdpU7Re+8IWldbuET5o0qaBXr16V9edzc3PTvffeu8HdvSMi8vLy0h133HHZv//97+Jp06YV1J878cQT57/yyitFM2fOzOvbt+9OX/nKV0qHDBmyeMiQIUv69OlT2diamZgzZ07O1Vdf3WnkyJFtZ86cmZemaaO18+bN+0wOd9KkSSsD4EOGDFmyPvuaMGHCyrV79eq169peN2vWrMa3lof1SCgdAAAAAAAAAAAAgNUsXLhwZc60S5cuTYa/u3btunJ+/vz52avOt27dumpNO4fX7Yi+ePHizxSed955C95999382267rVNZWVn2n/70p/Z/+tOf2kdE9OjRY/lhhx1Wev7558/bcccdV6zVgzXixRdfLDzqqKN2KC0tXat8bXl5eVb94/nz56+8rmfPnus1LD9v3ryMwuXLli3LWnMVfH5C6QAAAAAAAAAAAAA0KUmSZr3+lltumfWDH/yg5Pe//327559/vtV//vOflhUVFVkffvhh3t13393pvvvu2+bqq6/+8JJLLinJZP2KiorkhBNO2K60tDQnJycnHTZs2LyhQ4eWDhw4sKJjx47VBQUFaUTElClTWgwcOHDniIimdlFf36qrq1eOX3rppSktWrRYq5tvu+226zUcD40RSgcAAAAAAAAAAABgNe3atauqG8+ePTt3l112Wd5Y7ezZs1fu5N2hQ4fqVedLS0tzqqqqoqnd0ktKSnIjIlq1alXV0Hzfvn1XXHvttR9HxMfLly9PXnjhhcJHHnmk3cMPP9xx+fLlyWWXXdZzv/32K9tvv/3K1+4J/2vkyJHFH330UV5ExHXXXTfzggsumN9QXf3d0FfVoUOHlX3PnDkzt3///p9r5/b62rdvv/Iz7dy5c1WfPn2Ezdmk2JIfAAAAAAAAAAAAgNXsuuuuK8PdY8eObdlU7fjx41fODxo0aLVQeGVlZTJu3LjCxq6vrKyMt956qzAiom/fvmsMlefl5aVf+tKXlv7+97//8M4773wv4tOdy0eMGNGuft3a7tD+xhtvFNSNTznllIWN1Y0bN67Rz2Hw4MHL6sbPPPNM8VrdeC0NGjRo5drPPfdc0fpcG9YHoXQAAAAAAAAAAAAAVrP//vsvKy4uro6IGDFiRPvq6tU2QI+IiEWLFmWNHDmybUREnz59Knr16tXgLt733HNP+8bu9cADD7RdvHhxdkTEIYccsnhd+jziiCNW1i9YsOAzO5nn5+fX1I0rKioaTahXVVWtnFuyZEmD+drq6ur4wx/+0KGxNb761a8uKSgoqImIuOOOO7apqmpww/eMHHnkkYvrnuW2227rVFNTs6ZLYKMSSgcAAAAAAAAAAABgNQUFBenxxx8/PyLi7bffLrjkkku6rFpTU1MTp512Ws/S0tKciIgzzjhjXmPrPfjggx1HjRq12i7fM2fOzLn88su7R3waIj/rrLMW1J//3e9+166yssGce0REPPHEE63rxr17915ef6579+4rL5w6dWp+Y2v07du3om58++23Nxg8//73v99typQpje723qFDh+oTTzyxJCJi8uTJhaeddlqPxsLjy5cvT2bNmpXT4GQjaw8bNmxeRMSrr77a8vTTT+/R2B8JRER8+OGHOTfeeGOjAXpY39b6lxkAAAAAAAAAAACArcu11147+x//+Eebjz76KO/GG2/sOnny5IJhw4Yt6N69e+U777zT4rbbbttm/PjxxRERu+2229ILL7ywpKF12rZtW1VQUFBz9NFH9z399NPnHnnkkZ/k5+fXjB07tuVvfvObLiUlJbkREZdeeunsbt26fWaL8e9973vbXnnllT2+8pWvLNp3333L+vbtu7ygoKBmzpw5uaNGjWr14IMPdoyIKCwsrDn11FM/E2jfe++9l+Xl5aXLly9Prrrqqq4tWrRI+/Tpszwr69N9nXv37r2iqKgo/frXv774ggsuqFq4cGHO9ddf3+2DDz5occwxx5Rus802VVOnTs279957O7788svFgwYNKps0adJqwfo6N9544+znn3++1dtvv11w//33bzNx4sSiU045pWS33XYrz8vLS2fMmJH7wgsvFP/1r39t9+Mf/3jWOeecs6CxtRpae+zYscWvv/56y3vvvXebsWPHFp988sklgwcPXlZUVFSzYMGCnDfeeCP/ueeea/XCCy+07tu3b/kFF1wwf23Xh89DKB0AAAAAAAAAAACABrVt27bmX//61/TDDz98h/fffz9/1KhRbUeNGtV21bpBgwaV/eMf/3gnJ6fhaGpBQUHNI4888u7RRx+9w6233tr51ltv7bxqzbBhw+b99Kc/ndvQ9QsWLMh56KGHOj700EMdG5ovKiqqvvfee9/bfvvtP7Oletu2bWtOOeWUubfffnvnKVOmFA4dOnSH+vN/+9vfph9xxBFLWrVqVXPnnXe+f+KJJ26/fPny5JFHHun4yCOPfOZee+2115Jbb7115p577jmwwYeMiOLi4prnn39++lFHHdXnlVdeKZo8eXLhRRdd1Kux+nVRUFCQ/t///d/04447btunn366zbRp0wp+9KMf9WysvqioqPGt1GE9E0oHAAAAAAAAAAAAoFH9+vVbMWXKlCm/+c1vOvzlL39p9/bbb+eXlZVlt27dunrHHXdcdtxxxy0488wzF2ZnZze5zoEHHrhswoQJU66++urOzz77bOt58+a1KCgoqN55552Xff/735977LHHLm7ougkTJkz+61//2nrs2LHFM2fOzFuwYEHOkiVLcgoLC6u32267ii9+8YuLzz///Hk9evSoauj6W2+9ddYOO+yw/OGHH27/zjvvFJSVlWVXV6+e1z7mmGMW9+rVa8ovfvGLLmPHji1etGhRTnFxcfX2229f8c1vfnPBeeedN/+dd95psabPq0uXLlUTJkyYdv/997d55JFH2r366qtFixYtykmSJDp27Lhi5513XnbUUUeVnnzyyYvWtNaq2rZtWzNq1Kh3R40aVXTvvfe2//e//11UUlLSoqKiIikqKqrp0aPH8t13333pEUccUTp06NAGP0/YEJI0TZu7h61GkiTdI+LDiIhnj4no3HL1mq92GbmRuwIAAAAAAAAAAGBLUjzgstXOVS6sjGkXTKs77JGm6Ucbtan1ZOLEif/IycnZMT8/v02/fv3eaar2kZ/9e0B52YrcjdVbcykoalF5/E++8FZz9wFsXqZNm7Z9RUVFaVVV1ZTBgwd/9fOuZ6d0AAAAAAAAAAAAYIsjqA2w8WQ1dwMAAAAAAAAAAAAAAGy+hNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAACAzUUaEZGmadLcjQBszup9j6brYz2hdAAAAAAAAAAAAGBzsThN06rq6uqcqqoqGUiADFRVVWVVV1fnpGlaFRGfrI81fSEDAAAAAAAAAAAAm4tXq6url0REunDhwrbN3QzA5qj2+zOt/T59dX2sKZQOAAAAAAAAAAAAbC4ej4jq6urq0pKSkm3mzJmzzbJly/LSNG3uvgA2aWmaxrJly/LmzJmzTUlJyTbV1dWlEVEdn36vfm4562MRAAAAAAAAAAAAgA1t8ODB70ycOHFkZWXlEWmaVi1cuLD9woUL22dlZdUkSVKTJIl0OsAq0jRN0jTNqqmpyYqImqqqqpKqqqr5ETFy8ODB766PewilAwAAAAAAAAAAAJuTn0XER1VVVWdVVVUtyMrKKsjKysqPiKzmbgxgE1ZTU1NTUVNTUx4RaUTcFhG/X1+LC6UDAAAAAAAAAAAAm43BgwenEXH3xIkTR0XEPjU1NfvU1NR0j4iWzdwawKZsaUR8FBEvR8TLgwcP/nB9Li6UDgAAAAAAAAAAAGx2agOVH0bEY83dC8DWzn+qAgAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGObRSg9SZI9kiS5MkmSp5Mk+ShJkuVJkpQlSTI9SZJ7kyTZfx3XOzxJkr/UW+uj2uPDN9QzAAAAAAAAAAAAAABsiXKau4E1SZLkhYg4oIGpFhGxQ+1rWJIk90fEGWmarmhirayIuDMiTltlqlvt62tJktwdEWemaVqzPvoHAAAAAAAAAAAAANiSbQ47pXetfZ8dETdHxDciYq+I2CciLoiIWbXz346I+9aw1tXx30D6qxFxfO1ax9ceR0ScHhG/WA99AwAAAAAAAAAAAABs8Tb5ndIjYmpE/Cgi/pymafUqc+OSJHkgIsZERN+IOD5JktvTNH1h1UWSJOkbERfVHr4SEQemaVpeezwhSZInI+L5iNgjIi5OkuT3aZq+swGeBwAAAAAAAAAAAABgi7HJ75SepukRaZo+1kAgvW5+fkRcWO/UNxpZ6rz4bwj/B/UC6XXrLIuIH9Qe5kTE+Rk3DQAAAAAAAAAAAACwldjkQ+lraXS9cZ9VJ5MkSSLi6NrDqWmajmtokdrz02oPj669DgAAAAAAAAAAAACARmwpofS8euOGdlTfNiK61o6fX8NadfPdIqL352sLAAAAAAAAAAAAAGDLtqWE0g+qN36rgfkd642nrmGt+vMDMu4IAAAAAAAAAAAAAGArkNPcDXxeSZJkRcRl9U491kBZ93rjj9aw5If1xj3WsZfuayjpvC7rAQAAAAAAAAAAAABs6jb7UHpEnB8Re9WOH0/TdGIDNcX1xmVrWG9pvXHROvby4ZpLAAAAAAAAAAAAAAC2HFnN3cDnkSTJQRFxbe3hvIj4biOl+fXGK9aw7PJ644IMWwMAAAAAAAAAAAAA2CpstjulJ0kyMCL+Ep8+Q0VE/G+apvMaKa+oN26xhqXz6o3L17GtHmuY7xwRE9ZxTQAAAAAAAAAAAACATdZmGUpPkmTbiHg6ItpGRHVEHJem6QtNXLKk3rhoDcu3rDcuW5e+0jT9qKn5JEnWZTkAAAAAAAAAAAAAgE1eVnM3sK6SJOkaEf+KiK4RkUbEqWmaPrGGy+qHxbuvobb+bucfrnuHAAAAAAAAAAAAAABbj80qlJ4kSYeIeCYitqs99YM0Te9fi0un1Bv3X0Nt/fm31qE9AAAAAAAAAAAAAICtzmYTSk+SpHVEjIqIHWtPXZam6a1refn7ETG7dnzQGmoPrH2fFREfrEuPAAAAAAAAAAAAAABbm80ilJ4kSWFE/D0iBtWeujpN0+vW9vo0TdOIeKL2sH+SJHs3cp+94787pT9Rex0AAAAAAAAAAAAAAI3Y5EPpSZK0iIi/RMR+taduTtP08gyWuikiqmvHtyRJUrDKfQoi4pbaw6raegAAAAAAAAAAAAAAmpDT3A2shUci4su14+ci4p4kSXZqon5FmqbTVz2Zpun0JEl+FRGXRcQeETEmSZLrIuLdiOgTEZdGxO615b9K0/Tt9fUAAAAAAAAAAAAAAABbqs0hlD603viLEfH6GupnRETvRuZ+HBHbRMSp8WkAfUQDNfdERCY7sQMAAAAAAAAAAAAAbHWymruBjSlN05o0TU+LiP+JiCciYnZErKh9fyIivpqm6elpmtY0Y5sAAAAAAAAAAAAAAJuNTX6n9DRNkw2w5j8i4h/re10AAAAAAAAAAAAAgK3NVrVTOgAAAAAAAAAAAAAA65dQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY0LpAAAAAAAAAAAAAABkTCgdAAAAAAAAAAAAAICMCaUDAAAAAAAAAAAAAJAxoXQAAAAAAAAAAAAAADImlA4AAAAAAAAAAAAAQMaE0gEAAAAAAAAAAAAAyJhQOgAAAAAAAAAAAAAAGRNKBwAAAAAAAAAAAAAgY1tlKD1Jkl5JktyQJMnUJEmWJkmyMEmSCUmSXJwkSWFz9wcAAAAAAAAAAAAAsLnIae4GNrYkSY6MiAcjolW904URsUft6/QkSf4nTdN3mqM/AAAAAAAAAAAAAIDNyVa1U3qSJLtHxKPxaSC9LCJ+HBH7RsShEXFXbVnfiPh7kiTFzdIkAAAAAAAAAAAAAMBmZGvbKf3miCiIiKqI+HKapi/Xm3suSZK3I+L6+DSYfmFE/HSjdwgAAAAAAAAAAAAAsBnZanZKT5Jkr4g4oPbwnlUC6XVuiIi3asfnJkmSu1GaAwAAAAAAAAAAAADYTG01ofSI+Fq98b0NFaRpWhMR99cetomIQzZsSwAAAAAAAAAAAAAAm7ec5m5gI9q/9n1pRExsou75euP9IuLpDdHMhD2vjbZt2652/kcxa0PcDgAAAAAAAAAAgK3F8h+sdmrRikVxWVzWDM0AsDXYmnZKH1D7/k6aplVN1E1t4BoAAAAAAAAAAAAAABqwVeyUniRJfkR0qD38qKnaNE0XJUmyNCJaRkSPdbxP9zWUdF6X9QAAAAAAAAAAAAAANnVbRSg9IorrjcvWor4ulF60jvf5cB3rAQAAAAAAAAAAAAA2a1nN3cBGkl9vvGIt6pfXvhdsgF4AAAAAAAAAAAAAALYYW8tO6RX1xi3Woj6v9r18He/TYw3znSNiwjquCQAAAAAAAAAAAACwydpaQulL6o2L1qK+Ze172brcJE3Tj5qaT5JkXZYDAAAAAAAAAAAAANjkJWmaNncPG0WSJPMjon1EvJam6W5N1LWNiIW1h39M0/TY9dhD94j4MCLiww8/jO7du6+vpQEAAAAAAAAAAKBRH330UfTo0aPusMeaNmEFgHWR1dwNbERTat+3T5KkqR3i+9cbv7UB+wEAAAAAAAAAAAAA2OxtTaH0l2rfW0bE4CbqDqo3HrPh2gEAAAAAAAAAAAAA2PxtTaH0v9Ybn9JQQZIkWRHx7drD0ogYvWFbAgAAAAAAAAAAAADYvG01ofQ0TcdHxIu1h6clSbJPA2UXRsSA2vHNaZpWbpTmAAAAAAAAAAAAAAA2UznN3cBGdm5EjImIgoh4OkmSX8anu6EXRMRxEfGd2rrpEXFDs3QIAAAAAAAAAAAAALAZ2apC6WmavpokyTcj4sGIaBURv2ygbHpE/E+apks2anMAAAAAAAAAAAAAAJuhrOZuYGNL0/RvEbFLRPwmPg2gL4uI0oh4JSIujYjd0zR9p9kaBAAAAAAAAAAAAADYjGxVO6XXSdN0RkRcUPsCAAAAAAAAAAAAACBDW91O6QAAAAAAAAAAAAAArD9C6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGNC6QAAAAAAAAAAAAAAZEwoHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoAAAAAAAAAAAAAABkTSgcAAAAAAAAAAAAAIGM5zd3AVia7bjBnzpzm7AMAAAAAAAAAAICtyCqZtezG6gAgE0maps3dw1YjSZI9ImJCc/cBAAAAAAAAAADAVm3PNE1fae4mANhyZDV3AwAAAAAAAAAAAAAAbL7slL4RJUmSFxE71x5mR8S42vGeEfFxszQFsGXpHP/9L1L4bgX4/HyvAqx/vlsB1i/fqwDrn+9WgPXL9yrA+ue7NXPZEdGxdvxGmqbLm7MZALYsOc3dwNak9n/EX4mISJKke72pj9M0/ah5ugLYciRJUv/QdyvA5+R7FWD9890KsH75XgVY/3y3AqxfvlcB1j/frZ/bjOZuAIAtU1ZzNwAAAAAAAAAAAAAAwOZLKB0AAAAAAAAAAAAAgIwJpQMAAAAAAAAAAAAAkDGhdAAAAAAAAAAAAAAAMiaUDgAAAAAAAAAAAABAxoTSAQAAAAAAAAAAAADImFA6AAAAAAAAAAAAAAAZS9I0be4eAAAAAAAAAAAAAADYTNkpHQAAAAAAAAAAAACAjAmlAwAAAAAAAAAAAACQMaF0AAAAAAAAAAAAAAAyJpQOAAAAAAAAAAAAAEDGhNIBAAAAAAAAAAAAAMiYUDoA/H979xotZ1UecPz/QEBC5KJErQTlFlB0uVpqQLkooWJbBSoiWMRWI6hYlQpVEfEGWNRqXeJdEQWsq1wUghplKSBEQFwGKtYqGINCQS6Bcg+QGHn64d2z5s0wM2dmzpw55+T8f2tlzd7z7nfPni9PnrPfPXtLkiRJkiRJkiRJkiRJkqSBuShdkiRJkiRJkiRJkiRJkiRJkjQwF6VLkiRJkiRJkiRJkiRJkiRJkgbmonRJkiRJkiRJkiRJkiRJkiRJ0sBclD4JImLbiPhkRNwQEasi4p6IWBYR746ITSd7fJI0FURE9vjv8h76ellELI6IWyNidXldHBEvG8FXkaSRiIinRsQBEXFyRFwUEXfXYuWZA/Q37tgZEbMi4i0RcUVE3BURj0TEjRHx5Yh4br9jkqRRGkZcjYhFfeS1i3rob9OIOK7MIdxT5hRuKHMM2473O0vSRIuIBRHxwYj4YS3PfCgilkfEGRGxd5/9mbNKmtGGEVfNWSWpKSI2j4jDSsxaGhErIuL+iFgTESsj4vIS47bqsb89I+IbEXFzRDwaEXdExA8i4jV9jus1JdbfUfq5ufS7x2DfVJJGYxhxNSIW9pGvntjDmJwHkCRpAkVmTvYYZpSIOBD4BrB5hybLgf0zc8XoRiVJU09E9Pof1NLMXNihjw2A04Aju9x/OnBUZj7W3wglaWoZI26elZmLeuxnKLEzIuYC3wd269BkNfD2zDy9l3FJ0qgNI66WRTtn9PiRb8jMM7v0NZ8qru7UockDwGszc0mPnydJIxURPwZe1EPTrwNvysw1XfoyZ5U04w0rrpqzSlJTROwHXNxD07uBf8jMH3Tp60TgA3TeKPB7wCGZ+WiXPmYD3wJe3qHJY8DJmXlSD2OWpJEbRlyNiIXAZT1+5EmZeWKX8TgPIEnSBJs12QOYSSJiV+BcYDbwEPBRqsRpNnAY8CZgZ+B7EbEgMx+crLFK0hTyReALXa6v6nLtFJoPqH8OfBy4EdgROA7YFXgjcBdwwrhHKklTx/8CNwB/PcC9446dEbEhsJjmpN4FwFeAe4AXAO8Hngp8OSL+kJkXDTBOSRql8cTVhr8Bbuty/dZOFyJiM6qH1Y3FPV8BzgEeAfYF3kv14/dzI2KvzLxuHOOUpImydXm9DfgmcAVVfN0Q2AN4JzAPeB2wEXB4l77MWSVpuHG1wZxVkuAWqmf415by7VQLy7cBDgEOBuYC34mI3TPzF60dRMRRwIdK9UbgI8AvqWL3O6ji4v7A1+gen79Gc0H6ZcCnqeL086jy3B2BEyPi9sw8bcDvK0kTbdxxteYIYFmX6ys7XXAeQJKk0XCn9BGq7VqxFnhxZl7dcv3dVA9QYIxf70nS+q62M+VA8TAidgZ+RfUDrGuo4u4jteubAkuBBVRxeRdPqZA0nUXESVQTccsy886I2A74fbnc646+Q4mdEXEE8NVS/UJmvq3l+nyqycfNgRWln7W9f1tJmnhDiquLaO46uX1m3jTgWE6m2l0N4LjM/ETL9T2p4vMsupwkJEmTKSKWUO3We35m/qnN9bnAVVSbdgDsk5k/btPOnFWSGGpcXYQ5qyQB1YLFdjG1pc1BVIsaARZn5sEt158M/A7YgurHQs/PzLvrn1HuP7C8tW9mXt7mc/4KuLRUvwu8sj62EuevBZ4J3AfskJn39vRFJWlEhhRXF9LcKb1tzOxxLM4DSJI0Ap2OitKQRcTuNI9R/GrrgvTik8D1pfyOiNhoJIOTpPXTMTRPBDm6/oAaIDMfBo4u1VnAsaMbmiQNX2Z+KDOXZOad4+jmGIYTO99VXu8B3t1mrCuoTg0CmA+8cvAhS9LEGFJcHbcyN/DPpXo91dzBOjLzJzQfqOwTEZ2On5WkSZOZB2TmeZ0eRpeFOu+svXVIh66OwZxVkoYZV8fNnFXS+mKshZOlzYXAb0r1RW2avJFqQTrAe+oL0muf8Vag8VmPy0WLRr66Fnhr69hKv+8p1S3L50rSlDKkuDoszgNIkjQCLkofnYNq5TPaNcjMx6h2tYDqD8d9J3ZIkrR+iogAXlGqN2TmT9u1K+83/sB9RblPkmakYcXOsnPlLqV6XlkU1M6ZtbITe5LU2b40H2afVeYO2jmzVjauSpquLquVd2y9aM4qSX3rGleHyJxV0kzzYHndpM21g8rrA8AF7W7OzFuBS0r1JRGxWf16qb+kVC8p7du5oHwOGFclTW/d4uq4OQ8gSdLouCh9dPYur6uojnvpZGmtvNfEDUeS1mvbA1uX8tJuDWvX5wHbTdSAJGkaGFbs3LtNu8fJzDuA5aVq3itJnfUUV4FrgMbDFOOqpOnqCbVyu93UzFklqT9jxdVhMWeVNGNExLOAvyjVG1qubQzsXqpXZ+aaLl014uUTgAUt13YDNm5p9zil/8YPNXfzJHZJ01G3uDpEzgNIkjQiLkofncYv7lZk5tou7eoJ1i4dW0nSzHFoRPw6Ih6OiAcj4rcRcVZEdDtN4jm18lh/uBp3JakyrNg5SD/PiIg5Y7SVpOnujIi4LSLWRMTdEfHTiPjXiJg3xn09xdUy17CiVM1rJU1X+9TK17e5bs4qSf0ZK662MmeVpDYiYtOI2Cki/oVqMeOscunUlqY7AxuW8qjz1VnATmO0laQpoY+42uqUiLg5IlZHxL0R8fOI+FTZCb0b5wEkSRoRF6WPQERsAswt1U5HawGQmfdS7aYO8IyJHJckTRPPoZqMmw08EZgPvA74UUQsjogt2tyzTa3cNe4Ct9TKxl1JM9mwYucg/UTLfZK0PloIPB3YCNgKeAHwPmBFRBzV5b5GfFyVmfeN8RmNuPqUiHhC15aSNMVExAbA8bW3zmvTzJxVknrUY1xttRBzVkkCICIWRURGRFI9v18OfBJ4WmnyMeA/W26bzHy1XT+SNGUMGFdb7Qk8k+o0iS2pdlg/Brg+Ik6MiOhwn/MAkiSNyKyxm2gINquVH+qh/SpgDtXiS0maqR4GvgNcSvVr5IeAp1Dt7vMWqociBwHfjoiXZuYfa/f2E3dX1crGXUkz2bBipzFYktb1O+AC4GqaDzR2AF4FHAJsAnwpIjIzT2tzfyOu9jqf0PBEYPVAI5akyXEssHspX5CZ17ZpY84qSb3rJa42mLNKUu+uA96cmcvaXDNflaT+XUfnuNpwO1W+eiVV7rqWanH6AVQb2m0EfIhqsfoJbe43rkqSNCIuSh+NTWrlNT20b0zAzZ6AsUjSdDGvw646F0fEZ4GLgF2pFqn/E/CZWpt+4m79oYdxV9JMNqzYaQyWpKbFwFmZmS3vLwPOjYgDqB6mbAR8KiK+k5l3tLRtxNV+5hPAuCppGomIfah2RANYSfV3fjvmrJLUgz7iKpizSlInFwLXlPJsYEfg1cArgbMj4pjMXNJyj/mqJHV2If3HVajy0m1bNqkD+C/gwog4DfghsAVwfEScm5m/aGlrXJUkaUQ2mOwBzBCP1sob99C+cVzhIxMwFkmaFrod85qZd1Lt0NP4w/Polib9xN36EbHGXUkz2bBipzFYkorMvL/N4p769SXAyaW6KXBkm2aNuNrPfAIYVyVNExHxXKoFkbOoYt6hmbmyQ3NzVkkaQ59x1ZxVkjrIzPsy83/Kv2WZeU5mHky1I+8OVCf5Lmq5zXxVkjoYMK6SmavaLEivX/8Z8PZSjVq5zrgqSdKIuCh9NB6slXs52mVOee3lmENJmpEy83fAxaU6PyK2rl3uJ+7OqZWNu5JmsmHFTmOwJPXnNKCxCGifNtcbcbWf+QQwrkqaBiJie6rdzJ4E/Ak4LDN/3OUWc1ZJ6mKAuNorc1ZJKjLzP4BvUq21+FxEPLl22XxVkvo0Rlzt1TnAA6XcLV8F46okSRPKRekjkJmPAv9Xqtt0axsRT6KZ4NwykeOSpPXAr2vlebXyrbVy17gLPKNWNu5KmsmGFTsH6Sdb7pOkGaPsWtmYM5jXpkkjPs6JiC3H6K4RV+/KzNVdW0rSJCs/Lr8E2JoqHzwiM789xm3mrJLUwYBxtSfmrJL0OI34Ogf429r7k5mvtutHkqaLTnG1J5m5Flheqt3yVXAeQJKkCeWi9NFpLJycHxGzurR7dq18/QSOR5LWB52Ola0vVn92hzbtrht3Jc1kw4qdg/RzS2auGqOtJK3POuW10GNcLXMNO5aqea2kKS0i5lKdfrZDeevozPx6D7eas0pSG+OIq/0wZ5Wkprtq5W1r5eVUJ1XA6PPVtcBvx2grSVNVp7jaj3Hnqy3XnQeQJGkALkofnSvL6xzg+V3a1Y+RuWrihiNJ64Xn1Mq31cq/r9XbHc9V9+Ly+gfgpuEMS5KmpWHFzitr5Y79RMSfATuXqnmvpBkrIp4CzC3V29o06SmuAgtonrxmXJU0ZUXEFsAPaP5Nf3xmfr7H281ZJanFOONqr59hzipJ66rvwvtQo5CZa4CfleoeEbFxlz4a8XI1cE3LtWXAmpZ2j1P6f2Hjnsz84xjjlqSpqm1c7VX58WPj7/eB81XnASRJGj8XpY/OhbXyG9o1iIgNgNeV6n3AZRM7JEmaviJie+ClpXpjZv6hcS0zk+YRX8+OiBe23l/6eCHNXzp/u9wnSTPSsGJnZi6nubPPqyNi0w4fuahWXjzQoCVp/fBmIEp5aZvrlwP3l/LrIyLatAHjqqRpoOSG3wP+srx1Smb+W6/3m7NK0rrGG1f7YM4qSes6tFb+Zcu1C8vr5sDB7W6OiG2A/Ur10sx8sH691C8t1f1K+3YOLp8DxlVJ01u3uNqLvwe2KOXH5avOA0iSNDouSh+RzPwZcEWpHhkRe7Rp9k5gl1L+tL9kljRTRcSB5dfMna4/DTgfaOww8YU2zU6leUTiZyNidksfs4HPlura0l6SZrpTGU7s/Pfy+mTg460XI2JH4L2lugIn9iSthyJiu4jYdYw2BwAfLNVHgDNa25Rd1j5TqrsA72rTzx7AkaW6NDOXDTpuSZooZRfHxcBe5a1PZ+b7B+jqVMxZJWkocdWcVZLWFRGLImKTMdocC7y8VH9Pcw1Aw+k0f6jzsYjYquX+Damea21Y3vpEh49q5KuzgM+X++r9zAUaP0S6r3yuJE0p442rEfGkiFg4xv27A58r1QS+2KGp8wCSJI1AuCns6JSJvauA2VTHzXyEajf02cBhVDtNACwHFrT+IlqSZoqIuAnYiGrh+dVUx2w/QnVE7ELgKJrHxV4J7JeZq9v081Hg+FL9OdXk3I3AjsB7gMYDl49m5gnD/yaSNDoRsTcwv/bWXJoPNK6i5aFEZp7ZoZ9xx87ygGQpzQfj5wNfAe4Fdgc+ADwVeAw4IDMv6uU7StIojTeulocll1Hls98FfgGsLJd3AA4p/xq7SL4tM9v92JKI2IzqKO/G0bGnAedQ5cj7AicATyz1PTPzuh6/piSNTEScT3OnyB8Bx1A9LO5kTdnJrF1f5qySZrxhxFVzVklaV3k+tRlVbnglVY75UHnvecBraeaPa4D9M/OSNv0cBXypVG8ETqHa+Xdrqni9b7l2dmYe3mU8Z1OtI4AqXp8K3FbG8j6q/BfgqMw8rZ/vKkmjMN64GhHbUS1U/2+qkyiuBW6n+rH6M4EDgH+kuZndJzLzuA5jcR5AkqQRcFH6iEXEgcA3aB6j1Wo5VZK1YnSjkqSppfxxum0PTc8H3piZ93XoZwOqPySP6NLHV4E3Z+ZjfQ5TkqaUiDgTeH2v7TOz7VHaw4qdZaee7wO7dWiyGnh7ZrqDj6QpabxxtbbAZywPA8eO9fA4IuZTxdWdOjR5AHhtZi7p4TMlaeQiot+J6Jszc7sOfZmzSprxhhFXzVklaV19PJ+6FTgiMy/u0tdJVAsc287DUsXLV2Xmo136mA18i+YOwq0eAz6cmSf2MGZJGrnxxtXaovSx/An4MHBydlkI5zyAJEkTz0XpkyAitgXeAewPbEP1a78VwDeBz2Xmw5M4PEmadBGxD7APsAfVjjxzqX7M8xBwC/AT4KzMvLrH/l5OdRrFbqWvu4FlwJf9hbOk9cWwFqXX+ht37IyIWcCbgMOpju+eQ7WTz6VUx4r/qtfxStKoDWFR+mbA31HltAuAp1PF01lUu+/8iioenp6ZK+lBRMwB3gYcSrWL+8ZU+fH3qeLqzb2OV5JGbZiL0mt9mrNKmrGGtCjdnFWSaiLiWVTP8PeiimFPA7aiOuVhJXAdsAQ4r5dn+hGxJ1VMfFHp6z6qUynOyMyz+xjX4cAi4M+BLYE7gSuo1hb09KxMkibDeONqRGxMM1/dHZhHla9uAtwP/Aa4nCpfvanHMTkPIEnSBHJRuiRJkiRJkiRJkiRJkiRJkiRpYBtM9gAkSZIkSZIkSZIkSZIkSZIkSdOXi9IlSZIkSZIkSZIkSZIkSZIkSQNzUbokSZIkSZIkSZIkSZIkSZIkaWAuSpckSZIkSZIkSZIkSZIkSZIkDcxF6ZIkSZIkSZIkSZIkSZIkSZKkgbkoXZIkSZIkSZIkSZIkSZIkSZI0MBelS5IkSZIkSZIkSZIkSZIkSZIG5qJ0SZIkSZIkSZIkSZIkSZIkSdLAXJQuSZIkSZIkSZIkSZIkSZIkSRqYi9IlSZIkSZIkSZIkSZIkSZIkSQNzUbokSZIkSZIkSZIkSZIkSZIkaWAuSpckSZIkSZIkSZIkSZIkSZIkDcxF6ZIkSZIkSZIkSZIkSZIkSZKkgbkoXZIkSZIkSZIkSZIkSZIkSZI0MBelS5IkSZIkSZIkSZIkSZIkSZIG5qJ0SZIkSZIkSZIkSZIkSZIkSdLAXJQuSZIkSZIkSZIkSZIkSZIkSRqYi9IlSZIkSZIkSZIkSZIkSZIkSQNzUbokSZIkSZIkSZIkSZIkSZIkaWAuSpckSZIkSZIkSZIkSZIkSZIkDcxF6ZIkSZIkSZIkSZIkSZIkSZKkgbkoXZIkSZIkSZIkSZIkSZIkSZI0MBelS5IkSZIkSZIkSZIkSZIkSZIG5qJ0SZIkSZIkSZIkSZIkSZIkSdLAXJQuSZIkSZIkSZIkSZIkSZIkSRrY/wPgpLDrWXQ1PwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "domain_size_in_cells = (round(6*cm / sc.dx), round(2*cm / sc.dx))\n", "domain_size_in_cells\n", "\n", "scenario1 = create_channel(domain_size_in_cells, u_max=scaling_result.lattice_velocity,\n", " relaxation_rate=1.9, optimization={'openmp': 4})\n", "\n", "obstacle_midpoint = (round(2 * cm / sc.dx), \n", " round(0.8*cm / sc.dx))\n", "obstacle_radius = round(0.5 * cm / 2 / sc.dx)\n", "add_sphere(scenario1.boundary_handling, obstacle_midpoint, obstacle_radius, NoSlip(\"obstacle\"))\n", "plt.figure(dpi=200)\n", "plt.boundary_handling(scenario1.boundary_handling)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAFkAAAAVCAYAAAAtkUK4AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAEiklEQVRYCe2Z7VEUQRCGF4oATs3gyAAwAiEDwAiADLT8Bf8szQCJACEDMAKFDCAD8TLA5xm2x9lj97g74IAqu2que3p6e3re6fnYver6+roqy+7ubr+s/5eb+NyFRxt+81VBe3t7H6guFar/4uQI9Gsc85NzzoxEwzrsLfxjUtQ/1PuIO3W1B7f+Bf1prWsw9F8aiqo6RHde6gqb3+gXKfq7LG0eUsa3iXNEWUYeDPtG55hi3CvIV9bRN+L2OXQxvs7YsTFZB/BvPpNApiJ4P+DLKoNqvQAEyBWyk2HAG8jHha2Bqje4NAFw60twgUyEfIbwGZ6ehdu3ujXkBwO69nuAXwETOIF+hb4BMnXj3oevwRMhC6RAGVNOJuSxY69t38EH8zduK53u13LJtqlsYyiwQdHpp1DUXEDN2mhXLYAZONr014OXk+Ogrbf1j3o6og8zyUQwQQ5HeHHsOYm04xmz2rgcUyJ0k8buePRdBcibOEmpnTz++3G52JklEXZZLnROgpnS8IGtmZAzhPYNyq0liO4nZRXbHnzWtEqHFy19mywmhJkuTRQ7z4mFuPYW+BGgnG16C6LNjl5FXV7bK5aZZyaYObcmQMOCHFBjIuq26N/2nOXFc48pOka3tK7YY+Knid1xbS7wY6aVS5xqOxGIHaXlhVyC5Z53ic5sfk+JQ+EIXfINj2Bp7qTXnS2P1EBcZmgbOZaK9vN7xO7Y1wRZgMqspNokOrFDAX5Lcbn/opQUAK5gG6d0hfyHskUxOwPArozRX/gpfc9cJl7H6zYRY5k29gt8rAqyA/ME7iQ6Fdi0lyK7vZzB0+0CHsC45Mrs1t93ygH6sVYKtm986BmQB94xcX+dIJa22MW178HnLI3KrkY/dGxWau9WEABrE/uqcpBXHm1cLaMmMjLFbeZJiTG5qt36ym1k2tjFpCfInURHZqdLZ5hiu/BGEBMUfNjWum9B0S7owxS6tokatn20OjF6TXsNL29E1T1iTwksyM5SDHJ4AGaiW0NXe9i7HYyyCfC0c68bpsjkcbeV4efvXWeMboOL8JzByCZHxDtN7GJyJcgCEI4QGzSg5t4kL8nlLwUoLrE2H75BerULO/e6eBYxk3ae4rkfZAOcCdGXq/XWJwV0Ah9bxdixF0GbPJeCfE7x1tBGnq4nZQMB2bEA7AQocPfpU3h6w9EeWZtNypZ1CZ0H4xVcH4k67HzWm4kr6SEoDqVYMdknfZgcAuiLh6/XuaArxzhW7NnxjZCSZ45Pc86ih1j+vlAaovfqlpcQskH5PSOyM5ujE+TIQAfkNwonMRN127UbUDzonOA2O68/UutHnZum0b/0JXiSY7BfY3HlntCWbkJw+2lbhajT6hKoROPGXtjreyc+EFnxStYAJIyfihOPGe8KcUJeFBGzE+dkLs7Xkae3uGc4CvfJFwdwjaNbbdo+E8gMxKVTnqRPjjcxubyf/N48DRDEbhaLZ9qSIpP15b478vVaoxnSNkFO8sY1w9Du7Eoc8+fT/M+Ij9UzsP6CB3fn6B/bAOz82O+1N94Nqr/9zvLPVE3uNgAAAABJRU5ErkJggg==\n", "text/latex": [ "$\\displaystyle \\left( 360, \\ 120\\right)$" ], "text/plain": [ "(360, 120)" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "domain_size_in_cells" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAABMAAAAPCAYAAAAGRPQsAAAACXBIWXMAAA7EAAAOxAGVKw4bAAABIUlEQVQ4Ea2T0W3CMBRFTcUAEd3AbAAdgY7QFToCfCZ/SGwA3QC6QdsRyAgdAXWD9BxDEIQgJVWvdOPnm+fr5D07VFUVujLP80kzFy2DUX0YLlAUxYTpDk6Jfy5e1eEXesakPAnGYupjeHr5RnyATzDCezBHuOk3fIfLemPN/IIXGIjnDCbeQ0lOym1LeGgT/6r9q9lVA7p8Eb/5Sp6Ff4TW15qlhvT9Mk22LF7BBbHcE88YQy8zFj1DG5ZAbEc/4Vqhl5kLWqBhxDh2NiP5A+5bzGop62zGCg+0NWtipMBGZR+zDQvGTSfmFt+63dTMdou02zE8P9eYpULXCnNvjEi3YuBtR/RyC3fxVzw3FtY6bRgTiCOBx0G4oXd1gZ46/AvpNYSy4bXeMQAAAABJRU5ErkJggg==\n", "text/latex": [ "$\\displaystyle 15$" ], "text/plain": [ "15" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "obstacle_radius" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAE8AAAAVCAYAAAAHIbMXAAAACXBIWXMAAA7EAAAOxAGVKw4bAAAESklEQVRYCe2Y0VEbMRCGDZMCHNKB6QBCBYEOgFQQ6CAZnuCNIR0QKmCggyQVMKED6ABwB+T7lFtFPgtj48vED+yMbqWVtLv3a7XSXe/x8bFXlsPDw0HZfq3/waeGy3KvoKOjo8801wrRa/UvAoMGnyxZMrIkOrZhG/AvSdB6IBfUC8o69WGru4dsgCzmvqd+bxv5dWXsSSO7g69SThh32x7XdRsbvsM+fL/UTbtP+6CQ2b5A/qOQ9WgbXEP4N+UJPBoO/glfVxjUyM9oC4SAaPwt8hHwaAvcKXwLnoi6AGlsi3p2gvovZMfwSwfCta3Mcf8UQPTfYOcavgPPRFvf24AaKOfIk58xmLa+foAPY9v6oqcxILgDKDsUFZ+HvMKd3zZuFAqyTiRCzx6VPjw7RN0xtsfsp0kdPbDjQo5R45OAtOkTgjIao18/fd9egLeLkhSKMWJGvsn4G3QYRSUZcYJlZEqu+Ng2RnZF2azMRzw/odcd4yJZ2mTayDum3dluo0ucxKu/zMNcN+92EaRbdNWc036AKsimgDaFffv/BX1sXrqm24Xbpv87Jfx0XHU3Ngr0d/cND1HPOanpnIlhdCSHFJNd8R795pnSsWLISHVlpNVBA7tu1ydTAv2XFsYYRA/UTTdGo2Dm9EK7JPHactt6EJhIOyUMC5zbNU7gAOap6NT+NAA7birCB+2btyOyq/Pod/EjbRlx7oBaeon54jUQPB2ubaUY+FLuQeGqfp1BwbsZxk4z1GtJgPLkeMYYdS6qEWdUCbo5XHmNxCuBZ0RMioba5IkyjLpNzIHldp60QBGVdxMVz9CJbU/2J7drqGrGeU3yTqrPprHw+4x2bTcYyX0jr1PCmE6vNE5k3bRjgWrOhGzi9srKnqlgy8jxlJ9Gn9s0UkvSzDxznVGoX7VDLAWcB4YREc5TfTlh1DBfhcfK9aj7InJfJLaEopIi8uY6uAqF2vRrKd8xm741uJ9ZyvXnmCLIQ/gI6S9FEMO3sl+87gVPJekFy95Z6xjSsdrnnYBG3tHpdMFs6ffLxhM5vwT16ku15lWbzHURxhYC+YNyeLm4HigCWotSQRrTg0xAb5d5XFM2KM9RJPOxldA4kwXGF/ZTJxdkJu0EClwQ7+E5EVPXwV2KN/pEjcxrQ+3mH8NewrVlKUkg/Y4dkdM2/XhdqYGaFnuJXy1GjJPd42OEPELfva8BwVahilNEwdPRjaxGRpTGElFXh9E3pHhAuHB+66o3U6PTdvVHRB44RQVdHhwucOQvt+MV8nQTgIvBAcUUFjT2YyA6Gt/20788ALwRxEX7d4dP25T+IvmFPwPKjT65bSUjYeTDPkn//8McaoQuEnkyp7ydwMNBt59J09BeCMIXt3dn974uXqrBR5xSuorIU7eJ89lLZRdOTKljDydn+TqZUu1cw8Qn79D8J1mVDbL+YVg0p+d64y4mg4k/GPzczKfvb8cEaioADf0lAAAAAElFTkSuQmCC\n", "text/latex": [ "$\\displaystyle \\left( 120, \\ 48\\right)$" ], "text/plain": [ "(120, 48)" ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "obstacle_midpoint" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Finally we run the simulation and create a video of the velocity magnitude plot." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "if 'is_test_run' not in globals():\n", " scenario1.run(30000) # initial steps\n", "\n", " def run():\n", " scenario1.run(100)\n", " return scenario1.velocity[:, :]\n", "\n", " animation = plt.vector_field_magnitude_animation(run, frames=600, rescale=True)\n", " set_display_mode('video')\n", " res = display_animation(animation)\n", "else:\n", " scenario1.run(10)\n", " res = None\n", "res" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.2" } }, "nbformat": 4, "nbformat_minor": 2 }