{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Advanced Quality Analysis Example" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from IPython.display import display, HTML\n", "display(HTML(''))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Example of running an advanced quality analysis of Cl2 in Net2 using EPANET-MSX -- more information on how to run advanced quality analysis using EPANET-MSX can be found in the [documentation](https://epyt-flow.readthedocs.io/en/stable/tut.quality.html#advanced-quality-analysis)." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Note: you may need to restart the kernel to use updated packages.\n" ] } ], "source": [ "%pip install epyt-flow --quiet" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "from epyt_flow.simulation import ScenarioSimulator\n", "from epyt_flow.visualization import ScenarioVisualizer, black_colors\n", "from epyt_flow.utils import to_seconds" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Create a new scenario for chlorine analysis in the [Net2 network](https://epyt-flow.readthedocs.io/en/stable/epyt_flow.data.html#epyt_flow.data.networks.load_net2):" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "scenario = ScenarioSimulator(f_inp_in=\"net2-cl2.inp\", f_msx_in=\"net2-cl2.msx\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Set simulation duration to five days:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "scenario.set_general_parameters(simulation_duration=to_seconds(days=5))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Monitor the \"CL2\" bulk species at every node by calling [set_bulk_species_node_sensors()](https://epyt-flow.readthedocs.io/en/stable/epyt_flow.simulation.html#epyt_flow.simulation.scenario_simulator.ScenarioSimulator.set_bulk_species_node_sensors):" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "scenario.set_bulk_species_node_sensors(sensor_info={\"CL2\": scenario.sensor_config.nodes})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Run the entire simulation:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "scada_data = scenario.run_simulation()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Show \"CL2\" concentration at the 11th node over the entire simulation:" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjkAAAGwCAYAAABLvHTgAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/TGe4hAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC3V0lEQVR4nO2deZwU1dX3f9Xb7AswwzAsw6YIyCariLtEzOISoiHGBXmUJCpG5dEIRiVqImiMmiivqE/QLBqJj2geo5LoICpIAFkEFFEQmAGGZRhmYLZe6/2jp6qre7p7qqrv1t33+/mQ2D09Pbdv17117jm/c46iqqoKiUQikUgkkgzDwXsAEolEIpFIJDSQRo5EIpFIJJKMRBo5EolEIpFIMhJp5EgkEolEIslIpJEjkUgkEokkI5FGjkQikUgkkoxEGjkSiUQikUgyEhfvAfAiFArh4MGDKCoqgqIovIcjkUgkEonEBKqq4uTJk+jduzccjuS+mqw1cg4ePIh+/frxHoZEIpFIJBIb1NbWom/fvklfk7VGTlFREYDwJBUXF3MejUQikUgkEjOcOHEC/fr10+/jychaI0cLURUXF0sjRyKRSCSSNMOM1EQKjyUSiUQikWQk0siRSCQSiUSSkUgjRyKRSCQSSUYijRyJRCKRSCQZiTRyJBKJRCKRZCTSyJFIJBKJRJKRSCNHIpFIJBJJRiKNHIlEIpFIJBmJNHIkEolEIpFkJNLIkUgkEolEkpFII0cikUgkEklGIo0cCXHqmtrwye561DW18R5KViHnPb2R3x8f5LxnNlnboFNCjrqmNuypb8HAsgJ89NVRzF++DSEVcCjAwukjMWNCFe8hZiRy3tMf7TvcdqAJj777pfz+GCHnPXtQVFVVeQ+CBydOnEBJSQmamppkF/IUWLahRr+5xsOpKFg97wJUluSxHViGI+c9/Un2Hcrvjx7ZMO/GA1C6f5Z4WLl/S0+OxDZ1TW1Jb7QAEFRV7K1vzciFxgs57+lPV9+h/P7oUNfUhnnLtyHR0T4T5n3Zhhr9M0rvlNTkSFJgT31L0hstED4ZDSjLZzOgLEHOe/rT1Xcovz867KlvSWjgAOk/77FGXEgF7l2+Pav1RtLIkdhmYFlBp+cUw387FOCR6SPS+lQkIl3NuzNL5j1WMGp8nOxnIlDVPfmNNBu+Px7065Z4ThWk/7zvOdrZiNO8U9mKDFcJRmws1fgYgFBx1p5FuSjMcaHZGwAQPgU9Mn0E/lD9NQ40tuPpq8/Ad0f15jzKzKNXcS6Kc1040R6Z94evOB33vbkdIRV449YpGNW3lO8gKZBMaP39M/rgjc0HEFIjBp+Kzj8TxX1f2xBtbDkVBddP7o8XP9mL0jw39/FlKodPeKMeOxUFl42uxBtbDmJEn5K0n/ddR5vjPn+gsRWf7Fbj3keS3XPiPU43pJEjEEZBXFcbtwgb9Ya9DWj2BlCU68Kz14zF4J6FqCzJwz+31uFAYztavEGu48tUdtSdxIn2AHJdDjx3/TgMqShCZUkeFn+wGwca2+APZkYugXFzXbXzKO59I76WIqQCr286oD9Wk/xMc9+fO6Sc64b9zrY6AMB3R/bCtWcOwICyfBTkuPDiJ3vR2ObHyXY/inLd3MaXqazaeRQAMHVYBW48eyAGlOXj8Akv3thyEEdPerv4bbHxBoJ44eNvAITvF8Z1cNdrW/XngcQHgK4ei3DfsYo0cgQhVojY1cYt0kZ9yem9cPap5frzA3oU4OOv67H3WAuvoWU0q746AgA465QynDekp/58n9I8HGhsw4HGNozr343X8CyRyFNp9NSQhre4NBhS8e72QwCAK8f3w+TBPfSfdS/woKHFh33HWjGiTwmX8WUy2tq5ZEQvfd5djrBq4/DJdvgCIXhc4qs44nlb/t8Hu1Hb0Ibyohy8OvtMHOkw2q5+4T/673V1AOjqsQj3HatII0cQzIhJjYi0UX9nVGXUz/r3COsN9jWkTxw4nVyyH3acRs8/rTzq+b7d8rB+L7D/uLjznijkZDxhxp5CScNbXPrp3gbUN3tRnOvClMFlUT+r6p6PhhYfahrSw8hJp3Vz5GQ7th84AQA4b0hk7ZQVepDjcsAbCKGuqQ39e3TWvIlAoto+Rm8LAJxzahkG9yzE4J6F+GR3PdEx8L7v2EEaOYIQT0yaDBE26qMnE2/UAFBzTNybrZHYMOHC6SNx7pByITfvk+1+bNx3HABwvsGLAwB9OkSVB46LIa6NxTjPsYZMov9OhAJAUcKnS6ei4IozeuPNzQcRVNWwwaQAquFnyzcd0A0o3uJSzQP6reG9OnkNBvTIx5baRuxLg7WTTusGAD76KnzDH9GnGOVFOfrziqKgT7c8fHO0BQeOi2nkJKrtE+ttAYA3Nx/A3dNOQ2VJHgaWFcDRsU5IwPu+Ywdp5AhCr+JcdMt343irH0CCjRuRG8Cvv893o9a8OHE36g6Dbe+xFqiqCkVROv2+KMQLE857fZt+A3UowD2XDMXIviVCbNxrdtUjEFIxqKwAVT2iN5s+pR1GTqN4Rk7sPFvdc2ONmkemj8C5Q8qxt74VA8ryUVmSh7umnaY/BhD1s8E9C/HYip2YMKAbV01ByOAB/e6oXp1+XtVxg90neKjXzLoRzehZtTMcqoo9HADhtfPN0RbsT4O10xUhFbq3pbIkDwunj8S9y7cnPABo95hEjzUjKl2zZaWRIwj7jrXieKsfLgfwwvUTMLSyqNPGHQiG8K0nPkJ7IITRHLNnwht1h3Ay3kbd4ck52R5AY6sf3Qo8TMdnhe0HTnTaPFQgqs7Ewne/BMBeeBcvFKAJJ8+LCVUBYntyrIZjzRg1AKI2XG1TNz7WOPfUcjy2Yid2Hm7manhvrDmOIye9KMpxYcopZZ1+PkAL9QruyYn3fcaum3mvb9NvqLzXTiAYwsdfhz05sWFeIBzqBTJj7cR6W2ZMqIpaOwASHg7iPW7xBrHi80O47sz+aSc6BqSRIwwffx2+eY3r3x0XDI2cNGI37omDeuCjr45i3Z5jGN6bTzuKf39xCIdPeFHgccbdqHPdTvQqzsWhE+3Y19AqpJFT19SG7QdO4NEVX5r+HZbCu3ihgB+O76cbOeefFv80CoQ9OaJ50AaWFXQKURkNmdgTphmjxgpDKorgcTrQ1OZHbUNbJy8YK/6+oRYAMOWUMuS4nJ1+runZagTWs9U1tWHDnoYuX6fq/8N/7ZzSswhNbX4U57owpl9pp98R2QuaTMoQz/sSz9uS7AAQ72fGxxcMLceKzw/h84MnSHwc5kgjhzGJhHraKePcIZ1PGUYmDeweNnK+acCsKQOpjlXDOOa3t9bh12/vAAC0+IJ4c/OBuNZ9VY/8sJFzrCXupsKT2Pi2x6kgEFI7CWDjwUJ4Fy8UcO/y7ch1O3HoRDs8TgWTBnbv9Hu9OzbqVl8Qx1v96C6QcVlZkofTexdje8dGGc+QAUDMqInF43JgWO9ifFbbiM/2NzIzcrS1M6BHPl77dD9e27gfAPCvLw5h2YaaTmunqnv4hnawqQ3eQDCuIcQTY8sAIw50eHKS/C7PtfODsX0AAOMHdIfL2Tl7qm+38PUgomi/JM8Np6Ig2DHpTkXBLy45DaP6lib0vpBk0sBwFtpn+xvR5gsizyPWNdkV0shhSLwTxowJVQgEQ1i7+xgA4Ow4nhEjZw4K39zW721gclrvqhFkotNZ/+75WL+nQTi3e7z4diCk4o1bzkKrL4QBZfn46Kujegw7FhbCu3ju6aCq4vZXtwAAfEEV/9jS2bjMdTtRXpSDoye9OHC8TSgjJxAM6dl2D18xAlOH9UwYcqLFqD4l+Ky2EdsONOHS0fSLVCYyCICwxyre2ikr9KDA40SLL4jahjac0rOQ+jjNoq0d4+dRADzz4zMwtn+3qHUTz+jhuXb+3mFcfvDlkbjGpR7qFdCTs+6bBgRVFb2Kc/HkjNEYEEffFOt9IUn/HvmoKM7B4RNebK49jrMGJ79HiYb4BQEyhEQnjLqmNny2vxEnvQGU5ru7TBsd2acUuW4HGlp82HUkfnVLWmOOR6KS4Zr4WDQjJ94mGFKBVl8Ikwf3QGVJHmZMqMLqeRfgb7PPxPzvDIXRjmQhvNMyIpKRqB9NxO0u1rxvqW3EyfbwNf7jiVVcxIsj+4bX1me1jdT/VjyDIJZ4a0dRFF18XNMglvg4kQ6ne0FOp3WzZv6FWPSDkfrrFEai1a7Wjor4a0dbN3WN7QjSKM6UAh9+FQ5RXzC0HJMHlzFfO4qi6N6cdd90HaYUDWnkMCLRCWNvfaue2jjllDI4u7i7eVwOvdDbf0zExVMhlUaQmvhYtCwRTRtiJN5nqCzJw+TBPfDTcwfjyR+O6fjdfCbCOy0jIhmJjEvtRLpfMAHlRx0b9dkmrnFaaGL97QeaEKJ8I0tl7fTvHgndicSAOCG+2M+grRvN6LlsdLiG1n9NGcBs7Txw6fCkr4m3diqKc+FyhMPWR0620xyiZT7q0Guee2pyKQNNJnVEENbtOcZtDHaRRg4j4p0wFAADyvJ10fE5XYSqNCJWNd0LLtGYtecSidwAcQsCVpbkRdXISPYZNDRN0aEmL9RkR3OCzJhQhbFV4b/73xcP6fQ9JLpB9hXU7f6hSc0ZTQaXFyDPHQ4FfVNP1wua0topE1N8HNumxczaOa1XODlCK43BgstG99H/+42bJ5taO06HgsrSXABiZVjtP96Kb462wOlQcJbJ+wMNtHvO5ppGeAPp1a5HanIYoZ3O73l9m/6cQwHqmtrx2f4mAMDZp5o1cjSrmq4uJ3bMigIs6qh90ZXIrX+HgPLoSS9afQHke8S41HYdacaRk164HMCS68bj9N7FXbp/e5fmwaEAbf4g6pt9UUYSTRwd3+upPQujal0ku7n0LRXPk3O8xYet+xsB8D2NupwOjOhTjA17j2Pr/iac0rOI2t+qLMnDQ5ePwH1vbgcQXTfG7NoRzQu6oqO+z1mDe+C2C081JXLt1+GV2t/A7nr0B0MAwobLGf27m147fUrzUNvQhv3H2zB+ALPhJkXz8o/pV4qSPH69zAaXF6CsMAf1zV58VtuEiXESH0RFjDtPljBjQhWe/+gb7D7agm75Hhxv9eGGpesRDKmo6panK/y7YnS/UnhcDhw96cWe+hYMKqcnTpwxoQrLNx3Auj0NuO87w3SXc1ebW0m+G6X5bjS2+rHvWCuGVfJJd4/lX5+HN+opp5Rj6rAKU7/jcTlQWRLuC1XT0MrMyNE2a7fT0anWRaL5F7FWzupd9VBV4LSKIvQqyeU6lpF9SnUjZ/rYvlT/1uVjeutGzsr/Pl/XqXW1dvoLWitnRcfauWJMn6h+W8no13E9svRK+QLaugkfEkyvndJ8AA1CeUG1MC/PwwGg6XK64+1tdVj3zbG0MnJkuIoxWrTjgUuHweNUcKI9AACoPd6GZRtqTL1Hrtuph1DWUdblANCFtz2Lrd2g+ncXb7PWjJxLRnQuYpiMft01DwnDzbqjm7i7I+XVqHdIRHijFitcpW/UQ/hnZYzuFxYfa54lmhi7wWsaNTNoRk7t8VZhRLC1Da34/OAJOBRg6nBzhwMg8rkPn2xHu59NmMN4ONAwtXYE07MFgiGs6eg9Fa/4J2smGiII6YQ0chjj7ThlFOe6ozbBRKr/RJzZccH987ODpn/HLv6Ym61Z+gtWov5AYxu27m+CogDfsrBRA0C/buz7ccXbrLtC26ib2vw42c5OB5EIVVUjwkmOehyNkR3Zi58fPIFAx/zSQvMouBwKHBbE1pUleXA7FfiDKvW1bRbtcDBpYA9LpQm6F3iQ73FCVdkZ3tp+5bG4X/UVrCCgMSNxpADNWjXx8fo9DcJl/iVDGjmM0URb9c3eToWzEmXMxKOt41S0ZvcxTFm00rQXyA7aZp3jsmrkiCU+/leHpmDCgO4oK7QWctKbjjL8LNpN2OMyf4MszHGhND8cuxdhs16z6xgOn/Aix6VgwgD+Lu4BPQpQlOOCNxDCV4fpio8jYRNr68bpUHSjWhQvqNZvy6oHVFEU5mvHzuEAiIj2RSkI+PbWcOuccVXduGUkGtlc0wgA8AVDOO+3q6jec0gijRzGeP3hBRgv+8Jssay6pjb8cfUe/bGx5g4N7G7WonUj1zQF37a4UQMRAWUtww1QO5G6HNbmXa+Vw9ntvmxDDa774zoAgDcQLmDIG4dD0evlvLl5P1VPiU83Uq1vsyLpco6caNc73198ujUPKGAUH7M1clxOa4aB5gU92NEWhSfLNtTgxU/2AgBWdhQw5EldUxt++UYkaUalfM8hiTBGzuLFizFgwADk5uZi0qRJWL9+fdLXNzY24tZbb0VlZSVycnIwZMgQvPPOO4xGax8tXNW3Wz4WTh8JZ4fgxUw6pkaymjs08NvcrI3dyHmz/UAT1nfEkqednoKRwzBLxGfzRMqzD09dUxs+2V2PT/cdw7zXt0V5K0XZFDVB6vMf76HqBdUOB/aMnI5QrwBhAa1a8PDKrjMR46GHepkZOfbCVZUleVAUoN0fwrEWH42hJUVbOx/uPBJubtqBVSkDDVjfc0giRHbVsmXLMHfuXCxZsgSTJk3CU089hWnTpmHnzp3o2bNzI0Kfz4dvfetb6NmzJ/73f/8Xffr0wb59+1BaWsp+8BYIhVT9xpXjMp8xE4vmBTJedDRLpnttbtaa8PhgYxt8gZCtzZ4EyzbURG0aH3991HJhMk14XNfE7rP4bYSrAH4ZVl21AGHRu6gr6pra9LRcgG7jSN2TY/FmCxgadXL25CzbUIPH/7UTALCj7kTclghdUdWxdlgdEOyGqzwuB3oWhdsXHDjeZjmknQqirx3W9xySCOHJeeKJJzB79mzMmjULw4cPx5IlS5Cfn4+lS5fGff3SpUvR0NCAN998E1OmTMGAAQNw3nnnYfTo0YxHbg2fQeiY4w43OTOj+o8ltiIu7ZLpkU3D2s22vCgHuS4HQiqwpfY4jaF1iV5e3/CcnVNReWEOct3hz3KQlYDSZphQ8+TsZ+jJMdMCRIRNcU99S0paOCuk5skJz9PnB09wO8HXNbVh3vLUPQr9GGtydA+oxcMBwMcLmg5rR7vnsG5xQwLuRo7P58PGjRsxdepU/TmHw4GpU6di7dq1cX/n//7v/zB58mTceuutqKiowIgRI/DII48gGEycouj1enHixImof6zR9DiAdRFvLDMmVOH7Z4SbDN4wmW7JdKP3yQp//7QW7R0b/Yzn/8Mlrrz7CBk3Kx8Bpb2stkhHZXYbdaI2BmYq/LIkFS2cVfwpeHI+PxDen2oaWqknFiRiw56GTr237KydKj3U28pE6xKwuW4APt3I02XtzJhQhT/8aAwAoH/3PCZtOkjAPVxVX1+PYDCIiopoQVtFRQW+/PLLuL/zzTffYOXKlbjmmmvwzjvvYNeuXbjlllvg9/uxYMGCuL+zcOFCPPjgg8THbwUts8qhhNNKU6WqozKql1EqrJVNQzudaCTqukyTuqY2/M/H33R63u5NrV+3fHx1uJmJ+FhVVfhDqWWJsAxXJXJnL79lst7dnfcmDYRPpL/5/kj92nRQ9ILa9eTUNbXhyfe/0h/TDKkl+vvbD5zAYys677921o5mOJz0BtDU5kdpvvkUdDvonmeLgn2AT6g3XdYOEGnT0dRR3y0d4O7JsUMoFELPnj3x/PPPY9y4cZgxYwZ++ctfYsmSJQl/Z/78+WhqatL/1dbWMhxxGK+eiu0k0oqhrKPybv1Jb8rvlQw7wmPeQrVlG2pw1sKVWNVRiE6b7VRORSzd7sGQqp+irYYJNZd7fbOXWQE23Z3d8VgLoY7u181yOJY2V0+sQnmH3uJ/Zo6ndiL1BuyFeXmunWUbajBl0UrM/vOn2N/YjlyXI2WPQp7HqVcJZ7F2/GkWroqVHzgEXjuaTqmx1a/Ps+hw9+SUlZXB6XTi8OHDUc8fPnwYvXrFz4KprKyE2+2G0+nUnxs2bBgOHToEn88Hj6fzSSEnJwc5OeyEZPHQPDk5bjK2ZXlh+HPWN9MzckIh1Va2Ak+hmqYlMN4nFADP/PgMjO3fzfamwbIPj7FQpFVPTmm+G7luB9r9IXxW24hJg8yV4E+VGROq8O72Q1i18yj+++LThHZnV5bm4mizFzT3abtZibzWTjxtiC8Ywhu3nJWyR6Gqez6OnvSitqENozq6wdPCbskLIOLJ2XnoJOqa2pgZGFeN66f3CHz75+cI0wYnltI8N5wOBcGQimPNPu5tWszA3ZPj8Xgwbtw4VFdX68+FQiFUV1dj8uTJcX9nypQp2LVrF0KhyA711VdfobKyMq6BIwrtfnvalkRoVnV9M710R79hjt0Wxq2dToxRuYevOJ3JprGnvqWTliAEoHtBTkp/n6UmxyhSt7pZ//3TWv1a+9ELbLVQWhhWM8BFpWeHZ+HIyXZqfyMSrnJ28cpoOnnFwEbkGc+DFFKBVl8oZY8Cyx5WdrVsALC9o91H7fE2ploo4z6rhZtFxOFQ0KOA/uGaJNyNHACYO3cuXnjhBfzpT3/Cjh07cPPNN6OlpQWzZs0CAFx//fWYP3++/vqbb74ZDQ0NuP322/HVV1/h7bffxiOPPIJbb72V10cwhTFcRYKIkUPvYtM2asC6gHLGhCqsvucC5HV4rsb060Z0bIkYWFaAWEc1iZOwlkbO0uUOWAt3JNJCscrQie23JSpaH7YjJyiuHV14bD1sMmNCFZ6cMQZA+KbHwitGU5TN8oBgV/Ad1kJ9rT+mXWTVSCqeW9Zo952j0sgxz4wZM/D444/jgQcewJgxY7BlyxasWLFCFyPX1NSgrq5Of32/fv3wr3/9Cxs2bMCoUaPw85//HLfffjvmzZvH6yOYQg9XkfLkdJxGW31BtProCMFSMXIAoHdpPkZ3NBPdfrCJ1LCSUlmSh8mDIi0ESGUnaEXNmtr8aGqj2xcqkiGiWNJv8dZC2U17Z03Ek0Nvo7YbrtLQen3VHm9j0ocskTaEhAdJD/UyEO3bLXnBc+34A/Y9t6xhpQUlBXdNjsacOXMwZ86cuD9btWpVp+cmT56M//znP5RHRRbdk0NIk1Pgcerai/qTPlT1IP91RloLWGsyaGRknxL855sGbD/QhB+O70dyeAnpX1aIT75pwFXj+mLuxUOIbNQFOS70KPDgWIsPtQ2tKKHYNE8vTW8xQ4R30S67hdhY07Mo7Mk5yiJcZXMuuhd40Kc0Dwca2/D5wRM4k4G2asaEKjz8zy/Q7A3ildlnEvubLEX7+p5lcd55rh1t3TgdihB9qpJRpmtB2VeFtoPYO1GG4fWTDVcpikLddZhKQTONER3GwPYDbDw5QGTTGFReSFTLwOpE6rN5GuVdtMvuKZo1LDw5diuFGzm9d1iAynLtaJFSLdOIBFq46sDxNgSTVb0jgF1Dm3WRVSM+/VAj9roBoGcmSk2OpBOkw1UAfV2O3f5JRjQj54u6E3pnbdrQutmy0hakEuqYMaEKi6aPAgAM61XENMspfTQ54XVz+AR9T04qczGS4wGB5HdYUZwLt1NBIKRS17j4U9RCXdLR2+6W8wYzWzt2+23xgIUWlCTiz2gGEREep5GRQ+A0OrBHAQo8TrT7Q/imnk3DQU3TQrrHVD9GfXj8gdSMhaG9igAAjZS1Q7GkS7iqokN4XN/so+ZZSFWTAxi8oAfZVGhXVRWBUEQPRgqnQ9GLAtJeO6kezLSWGu0BNgcywFjbR+x1AwBlRTK7SpIA0tlVgCE+epJOfDSVJoMaDoeC03uHN+tt+9mcSH02NS1dwcqTk+pGrXkqjp70IkQ5PGAkoN/YxXa79yjwQFHCRRcbKHWcJhnq3X20GS1e+lVmjVk+VjUtXdHP0N6BJrpo3+a8lzMIZcaSLmFeAOhRoAmPpSZHEoPXT7YYIEDfk0PiNAoAp/fp0BYwyrCitWloGVZfdRQLo0UgxfGXFeZAUYBASMXxVnabUSo1Sljicjr0zZpWrRwSB4TyohxUFOdAVcNdwGljLF1AOnSi1cr5ZHc91bUTaetgb+3oRg7FUGYs6bJuABmukiSBTrgq7Mk51kI3XJWqscBaW0DKOIvls45iYXUn2qkWC0t103M7Heje0SOI5YmUhIaLFbTFx6k06DSirZ1tDNaO3fpMZjje4TF7c8tBymsnRS+onnnH3pOTFpqcjnBVQyu9UC9JxJ/RDIJKuKqIruvQR8hY0Nzunx88wSR8QuNkVNfUht/+a6f+mGaxMBLaFr5ud/G3Fi2kR+vETiK7CkAk1MvEyAmvG0UB0VTmuqY2vLv9kP6Y5trxBVILV+nXBct1kyb1pQCge3441KuqoBbqJYn4M5pBpGV2VYq1PjQGlxci1+1Aqy/IRHxM42bLsliY3RRyI1zc7oSuFxZUFNGtekwiuwqIeHI+P8AuXOV2Oog0EdbYU9+C2KMNrbWTuicnvG6avQFqRVZj0de74Fo2IBzq1bzE6RCyEn8nyiD0OjkUNDm06+SkulE7HQqGdzSd+5yBLkcvpkfQ5U6z7H0sJIw0ze3O1pOjnaLF36xpn9hJhUw1L+jXR06izUe3q3yqepZEsFw7gZD9FHIAKMxxIc8d9rbTbPthJJ00OUB66XLSY0YzBBrhKq0w08n2ANr95DdAktoWXVvAIMNKS8Em6VGgWfY+FhIp8MYMKxaoqpqmmhxKwmNC4aqK4hyUFeYgpAI7DtH15tBKZa4sycP8bw/THzsprh0tXGU3O0xRFOYhq3QK8wLplUaeHjOaIdAIVxXnufQb+TEK8VFS4SoAOL3DyFlDObsCiHT1Jb1pzJhQhbNPCZe6/8W0odSKhZGogMqi07aRgCGWlw6bdTllT5f2Haa63hVFwciO7MT/23KQcmYSPY/CdZP76//97h3nUls7ZLygbA8I6SQ8BgyenDRII0+PGc0QaGRXKYqCHnqtHPILkqQn5+Dx8Oa8o+4k1ewKgG7diYri8Okz1EllQA6i4SpmLnd66cc0qNCFx2KHegHo+piXPtnLJDOJxveX63bqYSCSe2AsJNZ+JNTL5oBAKouVFTJcJYlLRJNDLlwF0L3gSGWI1DW14Q8rv9Yf08yuAFKvGJyMbvluAEBTK71qwnq2BYFwFTOXe8DoyRF/s+5ZHEkVVlXyBquPUKn+uqY2fPDlEf0xi6w+klo2I9raaaS5dgjVJwJYhqvSU5NDSwtKkvSY0QyBRrgKQMSTQ+GCI6WxYJmZBNCNcZd2bNQ0i+yR6GVjdLnTuInHol0rpNOPaaHp2XzBEJWbLilNDsvMJB/FwwEAlHRk5bBYOyl5QSl7+WJJp7YOgPGeI8NVEgM0hMeA0ZND/oLTBbwpLj6W2RUAmRTsRJR2bNRUT6Mhci73Nn8QzUxaAtBJP6aFx+XQPQs0Tuy+jkNNqgZDumX1JUP3glLsqUbCYNAMYFbhKlpZbbTQO5EzzNy0izRyGEJDkwMYXIc0NupgeKNO1eWuZSZpS1gBvewKwNC/hkq4ioGRQ+BEnedxoijHBYCN2z3dNmqArvbCT6hJLNOsvhTTr7tC94JSLCJH4jo0hjJZkE5ZiYDU5EgSQKN3FWBs7UDBk0Owm/eMCVWYOXkAAOD7Z/Shll0BZEK4isz4yxm63dPN5Q5EwhKHKcyPj+ChZsaEKlR1Dxs1v//RGHpZfSmmX3dFqR6uounJSa3iMUC/5Ucs/hSrNLNGb+3Q4mPaANgO6TGjGYKPUriqvIie65BkCjkADOpZAABooVhJNBRS9XRmOuGqDvEkC5d7iuNnmUZOW89BA5qeHNKnc81zQ/OeQrsbNotwFYmsNm3dNLT49PejSbqlkGvNbQMhlep3SYL0mNEMgXa4imZ2FamNmsUJSdOzAHRORhFNjo+aoJfUDZJls8FIqCN9thWaAlM/IeGxBosQCm1NTmkeC+Fx6oZat3yPXqOKRUiGhAaPJR6XAyV5YYNV9JBV+uxGGYCeXUU8XEXPyCHdzbucQe0WTY8D0LnhaqdRf1BFK6Uy+8TCVQzd7rS9ADSgWfTNS3jtsDggBAilvSeilEEKuebFTeUzOByKvnZYHBBolryghSaTED2NPH1mNAPQ6+QQz66KxLmNBdlIQLpIlX5TaaaX1mycAxqbRp7bqW+gtE6kJNo6AIYbI4MmnTJcFUFVVeKhXhbfJW0BrNELSgvNg5aqroilLifd2joAdLN6SZI+M5oB0ApXdcv36HVJGgiLj7XFR2rM2unIFwjhRBsdXQ7tei2KolA/kZJo6wCwLQiYjht1BSXhsZ+CN7G8iH4BNmbFAGlqcgh5FMsZVj0m7TFnQRlFLShJ0mdG05xQSCXWyyYWh0NB94IO1yHhC45UQTONXLcTxblaWjOdzYNF9VDaaeSkPgPLTuRpmV1luJGR9CxGtbggFq6iH+qlLYBlmUKe6mdgWRCQZl0vWpSnSRp5+uxGaY7PsOnlEm7rANDT5dBwX2sCSlo33gCDTIUS/URKZ7Mm0dYBYBuuitxc0mej1m5k7f4QThIsmGjMyCEnPGYg2qd8QNDCVSfaAwhSSBMLhlQ9+yxlPVshSy9o+oV6yyhW2idJ+sxomqPpcQA6zenKKJXZJu3JAeh3+GUhgO2m18qh5ckhYzBop/8T7QG0++mIpDV8abhR57qdKNI8iwRP7NrhwEEwZKqtm6Y2P7XvMuKNo7N2tIwcgE4aeZQeL9UDQrG2TzE4IBDOYmVBD6nJkRjRMqucDoVKoa108uSUU67dQrugGRBJhW2k5Hb3h8gYDMV5Lt1ApZ0lko4bNUCnlhCNw0FJnlv3TtI6PeuaHAed79DtdOhVuGmI9o1GTsp6Ng6h3nQqv5AuVY/TZ0bTHFqiYw3dk0N4QdIQxEVCKHQ3apobRmkBXQElKYNBURRmWSLpKDwGgIqO8Oknu+qJdfb2UbgGFUWhXhKAZIXzROhrh4IX1Cj4Tl3PxkGTQ8mDRgNa9xzSpNdulMbQ6kCuoVnVpFs76FWaSWpyKJ+QWISraBc1I/kZIuFBum73iEGcPhs1ALR2VN9+5oPdmLJoJZZtqEn5PSOeHDrVzWndeEmXjIiH7gWl6MlxOpSUw4RauKq+2Uu9dUE6HhB0T04LvaKoJEifGU1z2inVyNHQLrhdR04SO40Chk2PpCenmLYmh0V2VUd5esqaHBKfgZXbPR01OXVNbfistkl/HFKBe5dvT3kN0RJh0zZYAyH6N1ua5RdIGmllhTlQlHBxQZoVmoH0FB4by4F8faSZ82gSkz4zmubo4SrC1Y41Pq87AQDYduAEsdMoYHBfk9TkFNLV5ERqfdDcqOl6ckgaDKxSYWnrOWiwp74FsWfQoKpib31rSu9LQ5MDgHoVXhaVd2munQAhLZv2Ht07xsou1Js+XtB/bDmg//clT31E7J5DmvTZjdIcmuGquqY2vLhmj/6Y1GkUIN+7CqCfCssilZl2k06ynhw2TTojvZrSZ6MeWFaA2NE6FQUDyvJTel9aRk4mhHppNukkrcdj1RaFRFNRltQ1tWH+8m36Y5L3HNKkx4xmABHhMflw1Z76FsSGREmcRgE6wmOtkuhJSmnNLOLbtIsBBgjqW1ht1OmoK6gsycNN5wzUHzsVBY9MH6F3/LaLl9Jc0D4g0G7rAACleVr5BfKeHJ/e0oFQ2r5W04tynal0Wzt76lsQK1Midc8hTXrMaAYQ6VtFfsoHlhUgVmNH4jQKGITHBMddnOvS34+G251FfDuiK/BRESVqn4FE6IdFpVyAXNo7a2ZMqAIA5LkdWD3vAv1xKpDuQK5B2ysXYLJ2Ir32SEPaWGBVEJBkmI0FNO85pEmPGc0AaHUgB8Kn0YXTR+qPHQqInEYBOic7RVEMJ1LymzWL9gJaUbOQCqKVcjVIzjszT06audw1tPlp84d0D12q0EghBxhocliUX6Ao2ietIaSdJKHhJ9zMlTbaPUezcxSQu+eQJj1mNAOgGa4CwqfRob2KAACLpo8ichoNhlS99DpxAWUhPTGsbuRQaM6pket2Iq+jPQfNVFgS4SpjKmxtQ0vK75eIdGzrANDxLNLW5NQ3+6i0RdAbw1LV5FAUHhM+lGmesx11TVT1JnqiQRrp2WZMqMLN5w8GAEw7vReRew4NpJHDCNrFAAGgb7ewFR0gtPlFlUgnngob3qxpdFRmlcrcjWIqLEmvSPWOI/p/n/fbVdSyINIxhRyI9SySuR5peUTKCj1QlPABpIFCtW0W2pASminkhIvq7epIjV635zjRrNVY0k2To3FKz0IAQDMFbzYp0mtG0xivn24xQCAi6D1MSCRnbCpK/ERKMa2ZVHPLrqCZCktK31LX1IZfvsEmCyIdu5BraJ5FUvVnaHlyXE4HehRoac00Qr0sakzRLAZITstW19SGV9ZHjBoWayddwlUaPQnfc2iQXjOaxtAOVwEgXr7f2EnZTbj2CU0BZaSgGV3XbymlVFhVVYmd7FhmQaTraRQgn5rtpWTkAJHCnzR0IiSz+hKheUBbfMGoPYYEJI0FmlmrsaTr2qmgnO1HgvSa0TSGdjFAINKDh1S6o7F6qIOwvoWmGFY/jVIuSqdrCwiHDYIhVd9cUzXUWGZBpKsmByAv6KXpEdHTmimsHR9BT0giinLdUDoukcY2smuHZBd1VmsnfKjRrpf0Wjva4aCpzU+lHAgJpJHDCBbhKtKeHJouVF2TQ2OjDpCNyyeihFJBQJJNBllmQfgYVMulBelmjLTCVYCxtQNF0T7F79DpUPTsRNIZViSNy8qSPDzyfTpZq0ai1nuahXqL81z6NU47A80u6TWjaQyLcJXuySGsK6Cx8Oh6cti4fmkJj31Rgu/UP8OMCVWYMaEfAODHZ1ZRy4JIV5c7AOIlDXzB8KGGzgGBvpFDu2p1pCAgaSOH7DX4o4lVul7rjzMnUFk7xgSPdNPkKIpiCFmJqctJrxlNY2h3IQeiazqQSC+lVesDiIz1WDOZsRrRsstobxi0OpEHKGS1VfUIu9i1opQ0SGvhsWY4EMr2o+nJKaepZ2OUIUdLtE+jLYW2V3Xq/0EIP+FDDWsi4mPpyclq9IrHFDU5PQrC6aUhFTjWkvoFR7OfSo+CHDi0sRJOI2fVB4ZWN2VjfF5RyOysZQURo5IW6azJIV0VmkZjWw2aFaxZtHUA6BUEpLH2exRqa4dWM97wmB1KOJSXbuieHEEzrKSRwwgW4SqX06FnXpDYALWNmob3yelQ9M2DdMgq0oWc7oZBKxWWRjfvHoXhsR6jUFtFg4VolRZaCKiekGeRZnYVzf5VrLph0yoISEPwXdaRsk/rgKCnvaehFweg3zQ2VdJzVtMQFuEqgGxqNk2XO0BPW8BKG0KrE7mPwo2mu75R0zNyWNUnokF3gxeURJE9mt7EckMKuRqb45wirFpzaMJj8qJ98uPX1w6lA0K6tXSIRQufynBVlsMihRwgmyWiiSdpbXi0auUEKIYKjJRSSiGn0fld8/AdayF/Y9Rg5QWgQbjIHrnrkcZ3qKF5ctr8QeKVZlk1WaXlBQ1QCJlqHud6ap6c9F03APmEF9JII4cRkS7k9MJVgPGCI6HJodO3SoOWtoCGJyQemifnRHuAqHjaTyEVWwtXtftDaPXRqWeRrlVbNUhm/NH0guZ7XCjMcQEgGyIgWYSyK7oV0MpMJB/60dYOjTYaADsdFC1oZvuRID1nNQ1hHa4iUWabtrFAOqNFI6LJoZ1d5db/m2TVY3+I/KaX73HpDUVphaxYtASgCcnNWls7ObS9oAQPCCSLUHZFiZ5CTiu7iqQXlG6ol1VGGy00z6KorR3Sc1bTEBbCYwAoJ+jJ0WPFlMZMq38VzcwWIy6nA0UdJ2qSm7U/QMe41E6k9QQy7+KRznVyALJGDsnKu/Eoo3BAIFmEsisi4So6mhyi4SrKmYk0Q5ssqOjwyB9v9RNv00GC9JzVNISVJqeCpMudcvhBu6nsPnqSaNM72jcYI6UU3O60PCK0U2FZFZKjBcnWDnp2lZPSAaFjrOu/OUZs7ZAuQpkMeuUXaKSQa4cDHxU9G6vwOi1K8936PYK0V54E0shhBIu2DoChrw2JcFWA7k3rs/1NAICvj7RgyqKVWLahpovfMAdLjwINASWt8WvdqxuoeXLS3O2eRpmJmnHw13U1xNYOjSKUiaCeQk5w3jVPji8QIi70BtJ/3SiKYsiwEi9klZ6zmoawClcZXe6hFMWwNIWkdU1teO7D3frjkArcu3w7kVOpn2G9Fj0VluCJlNbJTjNy6ikXNUvXzbqcoBDeRynkCITXzppd9fpjUmsnsm7IFaFMhNb3zRsIEW3sSOOAkOdxIt9DT8/GKm2fJrSkByQQalYXL16MAQMGIDc3F5MmTcL69esTvvall16CoihR/3JzcxmO1hoRI4fulGsWdSCkpnxK8lJcfHvqWxBrgwVVFXvrW1N+b5ZhExonUlpCRJrhKpaZObTQ26IQcLnT1FnsqW9B7PGFxNph+f0V5bj06r5E9Wy0DggUi2mme1YiYDxcS09OQpYtW4a5c+diwYIF2LRpE0aPHo1p06bhyJEjCX+nuLgYdXV1+r99+/YxHHFi6pra8Mnu+qiTlZ5dRVmT43Y69BN7qsWZaG7UA8sKEFvB3KkoGFCWn/J7s2rrABjK05PMrqI073qWCIVwlTEzJ103a2PGUqraCz27itLaiXW0kFg7LLUhiqLo2YlEvaAUyi8AdMXHPoYaQlpopUtELAgozG70xBNPYPbs2Zg1axaGDx+OJUuWID8/H0uXLk34O4qioFevXvq/iooKhiOOz0uf7MFZi1bixy+s02PlwZBqaJFAN1wFkGvgR1NXUFmSh4XTR+qPHQrwyPQRqCzJS/m9A4wKmgF0Gg1qm56LcB+bHhRTYaMyc9J0s9bWDYkiezQN7cqSPPz03EH6Y6eiEFk7ehFNRlk+2gGBjieH1gGB3tpJVw8oQK+wKwmEmFWfz4eNGzdi6tSp+nMOhwNTp07F2rVrE/5ec3Mz+vfvj379+uHyyy/H559/nvC1Xq8XJ06ciPpHmrqmNjz41hf6iVaLldcci7iRaYerAENBwBStah/lcuMzJlRhQEd37N//aAxmTKgi8r4s3e7aaXTXkWZiWS60xt9dO41S2KhZZubQwlhkL9UMK9ppwdPH9gUAFOY4sXreBUTWDutwY4EnPNd761uIvWcgRMcb1Z1i/6p0D/MCYnciF2JW6+vrEQwGO3liKioqcOjQobi/c9ppp2Hp0qX4xz/+gb/+9a8IhUI466yzsH///rivX7hwIUpKSvR//fr1I/459tS3INbLHVRV7Dp6Un/MwsghZVWzqN+gnZ6dBEXCtOrMxGPn4bCx/J9vGohludDqAdWDwUYNkPdAsaQnoRIMXsoHBM24bvUF9TolqaJ7EBmsm2UbarD1QDi78pdvbCeXWRmgUyMr0tpBanLiQbNpbKqk7axOnjwZ119/PcaMGYPzzjsPy5cvR3l5OZ577rm4r58/fz6ampr0f7W1tcTHlEhnonlWnA6FSadZUhcci2yZkrzwjZekpsXHyP1b19SGv38aMapJZ7mQ3vS0/lUNLb6UM+9iMW7UtDNzaFJGyMihnUJe3GHkhFTgJKG0ZlZZPnVNbZi/fJv+WAW5zMqIoUbpgEDDC8rwUEYLzZMjhccJKCsrg9PpxOHDh6OeP3z4MHr16mXqPdxuN8444wzs2rUr7s9zcnJQXFwc9Y80lSV5+M33O+tMSjtu5Cy8OIBRBJaqJod+jJ5GUbAAhbYI8UjkuUs5y4Wyyz0QUnGinXARtoDWMyh9N2qAXNVj2qfzXLcTuR1JDCcIHRA0LRttjwKLzErSaydyQKDhBSXfb4s12sG6vtkX5dUVASFm1ePxYNy4caiurtafC4VCqK6uxuTJk029RzAYxLZt21BZWUlrmKa4emKV7kr+039NxIwJVcz6VmmQcrnTrngMRNzujW3k2yLQ3qxpZYjRaNAJhI3V4tywDoK02z3da+RokBDtB4Ih/SZO84BAukYTq++QZmYlLeOSrmg//ddO93yPHqam1a3dLsRn9euvv8Z5551n+ffmzp2LF154AX/605+wY8cO3HzzzWhpacGsWbMAANdffz3mz5+vv/6hhx7Cv//9b3zzzTfYtGkTrr32Wuzbtw833XQTsc9iF11n0uG2Z1UIMPL3SQmPw8YZaW2IET0Fm0JbBNpehcqSPNzz7aH6YyehDDGam54xZEWSTNioAYPbPYW1YxRhU/WCEg71aocDFutm4fSRehq8AoKZlRQqHgORFHK6mpz09YI6HJGqx6IVBHSRfkOfz4fVq1db/r0ZM2bg6NGjeOCBB3Do0CGMGTMGK1as0MXINTU1cBjEqcePH8fs2bNx6NAhdOvWDePGjcMnn3yC4cOHE/ssdukeE79lVSNHo6I44nJXVdW2RkJPe6epySHcqE9VVaZehR9PrMLCd74EAFT/9/kYUFaQ8nvScrkD4Wvzm/oW4uLjTNioAUO4KoX5MTYppKpnyyfrBWWZyjxjQhUaW/1Y+O6XmDy4B7HMSlprX/PkNLSEK8k7CIrrM8UL2rM4F3VN7cK1diBu5KTCnDlzMGfOnLg/W7VqVdTjJ598Ek8++SSDUVkn4toMb5ReP5tqxxqaRe0LhtDY6ke3DqPLKrogjmLdkxLC4aqAIdjPIltBS4MFgMJcMsuJ5qZnbDZIEho9g3hA4jSqfX+KQjfTjHS4inWWT99u4fBUIEhOBK99BtLzrlU2D6lAY5tfP8iSgJb3iTWkZBKksTyrP/vZz/DCCy/g008/hc9HpwdOutNdb4SoeXLYhqtyXE49DJTKBRfR5NAbN+mqpwHGRekcDgUFHX1tWghludBq6wAYWzvQ8eSk/2mUnCeHdqaZtnaIhasYd8PWDgWkssMAQ2YiYYPB43LoRiVp8XHGrB1BjRzLR89t27bh5ZdfRktLC9xuN4YPH46xY8di3LhxGDt2bFRIKVvpEVN0jbXwGAAqinLR2OrH4RPtOK1Xka33YJHaqBljpDJEjHoIFg06gfBm3eIL4mQ7oVReivWJygroCCgzZaMuN2iWfIGQre+Advq4BumWIqwr72qFF5u9BPV4FNPgexR60NTmR32zD6f0JPe+mRLqjRShTfNw1Zo1a6CqKnbu3IlNmzbp/9544w00NjYCQFrXySBBJH4b48lhpMkBwifSnYdPpubJYbBZa+LJRsKnUYDdibQgxwXAS8yTQ6utAxDx5NASHqf7Rt0t3wOnAgRV4IuDTRhT1c3ye9CqcxRLJFxFSpPD1lAt6vDkNBM6HAB0+2+VFeTgm6MtxA8ItPptsSZjPDlA2IgZOnQohg4dih//+Mf689988w02btyIzZs3ExtgOtI95rQc0eSwCVcBEW3Bhr0NmHJKD1uZCywqHmviyVZfEN5AMOU5MrrcWRnbRfqJlJQnh2a4qkOTQzhclSkb9Wsba6FFPL//7CdYNH2kZVEsK0+OJtpP23BVx7pp8QaJvSfNWj89KDW4zRQvqBbq/aY+3OKGRLYcCSzP6gMPPICNGzfG/dmgQYNw1VVX4ZFHHkl5YOlMJLuqQ3jMIVylndSXbai13W6ARZ2cohyXXjODxGZNU8+SCE1bQMzIodTWAeic+UeKTNioO1XitVnB2hcMr3fq4SriwmPG4aqOdeMLhvQ9MhWCIRVBis15tbVDOo1cXztpLjzeUtMIANhb30qsxQ0JLM/q/v378e1vfxt9+/bFzTffjHfffVcKkGOIrUUSER6zuYjrmtrw4c6j+mO77Qb8DE6kDoeil6gnUSuHRzqmlmFFzpNDL/RTRlt4nMYbNalKvKy8WiW0hMeMvkNjZiKJkFVU/zQKa4e2aD+dQ711TW145oNItwFSLW5IYPlqXrp0KQ4dOoS//e1vKCoqwh133IGysjL84Ac/wJ///Gc0NDTQGGdaoVn8jW1+BIIh5tlVe+pbEJuUaWuzZmQwRKoep75Zs3a5AwZPDinhMcXTqNaD53hr+NokRSZs1KQq8bLwgAI0hMcda4dRg1WnQ0F+R2YiiQNCtB6PRiHNaK0lKVj12qMJzVYdqWJrVh0OB8455xw89thj2LlzJ9atW4dJkybhueeeQ+/evXHuuefi8ccfx4EDB0iPNy3olu+BooTd3cdb/fD62RYDHFhWgFg5ip3N2stYW0DC7U6rJUIyCklrcihmiJTme/QbeQMhwSqQGRu1Xom347HdSrzMsqvyyBbS5OEF1dYOicxEv7F8BJUDgubJIRyuYtQYlSY0W3WkCpFZHTZsGH7xi19gzZo1qK2txcyZM/Hxxx/jb3/7G4m3TzucDkX3TjS0+JiHqypL8nDLeYMj41EUW5s1q+JgJOt9+Bk15zRC3Mih6I1yOpROdZxIkAkbNRCuxHvdmf0BAFeO62urEq+PUe80LVzV5g8S0bToBwSGIUeSejbNM+l0KHBSyUzUCmlK4XEs2gFBQ2tOLYL4mPislpeX48Ybb8Q//vEP3HXXXaTfPm3Q47ctXubhKgC4cnw/AECu24HV8y5IbbNmVO+DRCqsn0Ftn1iIh6sob3qx2X8k0CvNpnG4SqN3t/DGHIxtMW8SFlmJQDgFWyEp2udwQNAzEwmsHZrp40Ak1EuvxlR6r50ZE6owsMNz8/sfjSHWqiNVLF/Nr7/+OsaMGaM/njdvHpYuXYqNGzfC6xUrP54nxhsJj+wqzTvS7g/pYlMrBENqpJNyOnlyeGRXaamwPlJ1cuh+hkizQXLrlXVLAJoU52oFKu19n6wOBw6HEhEfkwj1MtbkAGQ9Ofrap1QEVDu4NrX5o/qTpUomhHo1tL2FVSFWM1iuk/Piiy/ihhtu0B8vXrwYwWAQ7e3tcDqdGDZsGD766COUlpYSHGb60cMQEtDr5DAsBqhlLAHhRWnV0IlqMphOmpw01xUAEbc7NSOnkPyJNJM26uK88Pd5ot3e9ehlaPCV5LnR2OonckDw8QhX5ZBr7UA7O6w0zw2HEs4cOt7q0yv8pkomhKs0NKPb7tqhgeVZ/fzzz3HxxRdHPbdt2zZ88803WL58OdxuN5YsWUJsgOmKsR4Jj3CV06HoFUXtbIBGI4dZ5VYSGzUH128BJU2Oh1LvrTJDKJUUmbRRRzw59q5HmnWOYiFZK4fPASE8fiLhKsqhaodDQXcK4mPa650l2uHarheUBpav5rq6OpSUlOiPnU4nFEXBgAEDcOmll+Luu+/GW2+9RXSQ6UikfL6XS7gKSK1LsS8qHZPu4islWJ6eRzHAIr1yK9mKx7Rcvj1oCo8zaKO265ljlUIORMZK4oCgaXJYlgHQWzsQ6F8VoFh6QaOMQtXjTBHtA0BxbmpeUBpYntWysjLs3btXf3zo0CH0799ffzxmzBh88cUXRAaXzvSI0uSwD1cBqTW/NG7UtNsjkKz3wUr0aYS08Jh2Kq9mgJOs3EqznD5r9I3a5vXISpMDhEsCAGTWDo/WHIUEhccsdGE0Qr0062KxhnSBShJYntULL7wQf/zjH/XHubm5cDojYRiHwwG/X5wPyIuocBWH3lWAsfml9QXJotqxBkkjh2Zzy0QUENQVAPTd19q1ubejxwwJeNRYoYXmHWn2BRCKrXBmgkgKOf1rUBftk8hM1DPk2B8QiGhyAvQz/LRw1ad7G4itnYwK9ealFuqlgeVZvfvuu/HKK6/g97//fdyfr1mzBoMGDUp5YOlOlPA4jcNVLLQtJQSLmvHYMIzhKtVm2rER2iG3z2qPAwC+PtJCrMdMJrnctRCKqtq7+bL0JpLUs3GpFk4w1MvC0NZaOvx1XQ3xtZMZXlBNeJzGmpyRI0fir3/9K+6++25MnToVr7/+OmpqanDw4EH8/e9/x/z583HNNdfQGGta0cPQv0oLV+W62XpySlLwkLB1uUcU+UEbJ2cjuoHAIVwVUsOF2VKF5mZd19SG5z76Rn9MqsdMptT6AMIe19yO0LKdEymrSuEAWS+otnZY3myLaKSQUxp/XVMb1u4+pj8mt3a0PSv9146WmZjW4SoAuPLKK7Fu3ToEAgFcddVVGDhwIPr164cf/ehHmDRpEubOnUt6nGlHd71HkA+tPj6enFQyL3wcTqOqCpxMUbDGo15LntuplzRPdbNWVZVqYT1aPWa0jZqlFoomRbn2U2FZhu5S8dbGwrOtAwlNToDy2ifVE9CIqqoZGeoVKVxluU6OxhlnnIFVq1ahpqYG27Ztw8mTJ3H66adj5MiRXf9yFtAtP3LjPnyiHQB7TU4qIjAfw/CD2+lAgceJFl8Qja1+XUxpBx6aHEVRUJDjwsn2AJrbA+hZZP+9giEVWsSLxmat9ZgxGjokesxk0kYNhMXHR096baXCstSz0QhXsaxaTbJOjn4NUvKIDCwrgAJEGTqprp2AYSFmwtopTuFwQAvLs/rrX/8a77zzDg4fPgwAqKqqwne/+1386Ec/kgaOAZfTobuSeWdX2TFyWHtENMMm1c2aR/8dgFz/KtpNBitL8nDfd4frj52EesxkkngSMJxIU/DksFg72rohcXLm4QUlmZlIO1xVWZKHy8f01h/b7QloxNg5PRM0OSJmV1n25DzwwAN6SnGvXr0wduxYjBs3Tv//Pn36EB9kutKjwBPlRuYnPLaeecFSkwOEx3qgsS3lxcGrvQApIye6PhGdz3Dl+L546J/hMg8r7zof/XsUpPyemaTJAVIrCOhj2JCXaN83Hno2goU0dU8UxZYCZw7qgTe3HMTYfqVYfO3Y1A8HAeOhJv3XTrGhnZA3EGQevYiHZSNnwoQJqKurw6xZs1BWVoZNmzZh+fLleOSRRxAMBlFeXo6xY8finXfeoTHetKJHQQ52H23RH7MPV9mvocGqk7IGqc060oWc7YZB6kQaiDJy6HwGY2+f8iLrfc3ioW3WmXAaBYyeHPvZVSw1OU1tfoRCKhwphGl5tkRp9QURDKkpdQ9nUTlYS68vzHMT6bCtHWoUBVQ6p7OmKCfcNDasrwwgpzANjZx169bhpZdewr333osJEybgiSeewODBg+H1erFlyxZs2rQJmzdvpjHWtEMTH2swFx6nkl3FuKgeqSwR7WbLstYHQK5JZ8TlrlArwmjUXBjDY6mQiZocQPzsKs3ICanhuj6aB8oOfg56Nu1wAIS9OSV59sfPQkeoHTyMh5FUMBqWtIuussDhUFDYoU+00zORypjs/NINN9yAr776CqeddhrGjh2L+fPnIxgMYtKkSbj55pvx/PPPkx5nWtK9MMbIYazJMWZeWK3fwlJ4DJDLEuGlDSGVJcLC5W68iZHerFmKVmmSkiaHoZGT646ku6faiZxHhlyOy6l7/1INWbFo66CtywChwwGPzu+0SbX3G2lsXw2FhYV47LHH8Omnn2L79u045ZRT8Oc//5nk2NKeshhPDmtXvuYdCYRUPY3dLKw3PFIFASOaHLabBqmqxyyKMCqKorvGAynWJdLgpYWihbZR2+lfxdqrRUrsyaugI6lQL4vxa0a8FhZPFdqd03lQkkKolwYpzWwgEIDX68XVV1+Nvn37YtasWWhoaCA1trTHGK5yORTmIZQ8d+SUZDVryddRpZm5JsdGCwojtDMsEkGqciurarmaN8dPzJPDJ6uNFlpRMzunUdb900pJHRB46dl08TEpLy698UfCVYTCvBz6hdFGtIKAljU5ixYtwrZt27Bt2zZ8+eWXyM3NxahRozBx4kT89Kc/jepQnu10N8QjWetxgPCJvTjPjfpmLxpbfehTal4ox1yTo/fgSc9wVREx4TGbTc/tdMAbCKVcYVqDdXiTNqnU+9Czq1h5cggfEHhlJtrt+q7hY7B2tHAVucNBZnlAAfHCVZaNnHvvvRcDBgzAzJkzcfXVV2PIkCE0xpUR9DB4cnIYt3TQKM0PGzlWrWrWGx4x4TGnVGby4SrKnhzN7U7oRBrg5AWgRZEuPLYRrgqwDUGQCFcFQ6pu8DIX7RNq7cDigKOtG1Jh3kxbN4AxXCWGkWP5ajjnnHNw7NgxPPjggxg3bhymTJmC2267DS+++CI+++wzBIOp9+7JFHoYhMc8PDmAfQ+JV9+o2Sw+XZNDzMhJ03AVg07KQCRcFSCmLcjUFHI74SrGBwQCon0/g9IFiSgitHYCDPR42r5CSrCfmeEqsQoCWvbkfPjhhwCAr7/+Ghs3bsSmTZuwadMmvPzyy2hsbEROTg5GjhyJ9evXEx9sumHU5PAycuye8iJuVDYeqEidHDIZItzCVYQqHtO+QRLPEsnUcJXgKeQAGU8Oz/YCmicnPcJVZD2gmVYpHDCuHTGEx7Z7V5166qk49dRT8aMf/Uh/bs+ePfj0009lnZwOuuUbjRw+4apIvN6q8Ji1J0fbqH1QVdV2zQhe2QoFHs3ISc2TyWrTI+1292VYlogmnjzpDVgusqeL9lnXmErFkxOgX2k7EeRaomheUJp1cjoOBzK7KiElmmhfkHCVbSMnHgMHDsTAgQNx1VVXkXzbtMXtdKAkz42mNj/zGjkaduvPsBZPahu1PxhOd9c0LlbhVXcikgab2sJmkUIefn+ybvdMbeug2iiyxzpcVaL3fbMvPNa+PweHyrvEUsgZXIMuwtlVvEpe0ES0TuSZYz4KiqbL4afJsdfagbUbNZV0dyMsXNbxIHUaZZVdRdLtHgypelfzTNHk5LqduifG6mbNOjORRCFNnhWriwh7cmjOO+nsKl77FU1Ey67KnJkVFC3Dile4KpK1ZO2Ux7pBp6IokdBaCv2r/IwzWzQiwmMy4Sra8+4kKDyOFq1mzpZiR1tgzFJiLTxORZPDUzheSCgzkYUez6OHq6SWLRGp9H2jQebMrKBo4uMWbwB1TW3M/75dUSLr0yhAZrPmlZJpTIMNpbABavNOu39QJFyV+mbNonM6D3RdjoUQpM+obUmjvm8BjtqQwg5jkly4ikEKuRQeJ4RUBW5SZM7MCsrxDhfy5tpGTFm0Ess21DD9+yU2s5Z4FHcjIqDkHK4CUmvSmY7C42jRagZpC3Ktn0iNBh8zTQ6BmworLVg8CnPCXm5ydXLoa3JIt3Wg2TmdNcZq4VZ7JtIgJeFxdXU1qqurceTIEYRivvSlS5emNLBMoK6pDRv2RNpchFTg3uXbce6QclSWmK8+nAp2a2jw8OTo2oJUNmtO7t8clwMuh4JASEWLN4gim92gWYXb3HoKOYlwFf3O6TywI6D0cTD4NN1dqy8IbyBoKzSufYc0G8MmojCHkCeHQc0Zbd2oajg0mapIO5M1OYGQijZ/EPkeovlNlrE9sw8++CAuvvhiVFdXo76+HsePH4/6JwH21Lcg1o4Nqir21rcyG0OJTaU7j3LjJJp08sryURRFzwhLpQdPgJGeQ9uc/SQ8OQw6p/OgONd6KqzxcMDK4CvKdUH7U3a9Oaz7bRkhVfGYhXjaWKSThPg4E9dOvseph9tFqJVj28RasmQJXnrpJVx33XUkx5NRDCwrgEMBjPcRp6JgQFk+szGUdqSXnvQG4A+GTG8ArIXHAJkmnayMhHgU5rjQ1OZPqVYOq7BBRFuQ+kbNM9RBkyIbwmPNE8fy+nM4FBTnhktVnGjzo2dRruX34FkCINK7KjUNBws9nnH/JBnqzaRwldYzsaHFh6Y2P3qVWL8eSWJ7Jfp8Ppx11lkkx5JxVJbkYeH0kXB2HLOcioJHpo9gFqoCIqdRwKLbnUO2hRZa++rQSdsibZ7ZCiSadGoud9r9g9wEs0R4egFoUmyjqBmPMC+QesVwXlo2ILpaeCoaDm3t0NyzjAkBZEK9mSc8Bux5QWlhe2ZvuukmvPLKKyTHkpHMmFCF1fMuwN9mn4nV8y7AjAlVTP++y+nQ61BYcWWzbjIIALuPNgMAPth51LZIW89O4nAiJRGuYhUm1HtXEcgSYaGF4IGdeh8+Dp4cAMjvaAD8zdEWW78f6ZnGL4U8pAJtfvteUBbVg40aHBI1pjJRkwPYl0nQwHa4qr29Hc8//zzef/99jBo1Cm53tNDyiSeeSHlwmUJlSR5T700sJflunPQGLAl69dL0jBZfXVMb/rHloP7Yrkibh5ZII1IQkMBGzariMYEsEZ6F5Ghip0lnpL0FOyN72YYa7Dh0EgBwz+tboUK1fJjiWXk33+OEonRUl24P2Baqsii/oCgK3E4F/qBKtMZUpq4dEdLIbRs5W7duxZgxYwAA27dvj/pZJmVYZAIleW7sP95m6YLTi4Mx2qyTibTNGjnGyrs8Ng0SrR1YhQ2cBCseBzJUk6O73C1oclh7cuqa2jB/+Tb9sQp7hwOehqqiKCjMceFkewDN3gB62nwfZuUXHA74g0EiXlAWndN5IFLVY9tGzgcffEByHBKK2Kk/E9ms2VRqHlhWoJ/mNKyKtKMq7/LIEvGkniXCvE4O0RTyzDyNWvLk6EJSNutmT30LYmVVVg8HALt2IokoMhg5dokczBisHT+Z7KpMDVdF9Gz8s6sya2YlcSnVU7PNZy2xFpNWluThxikD9cd2RNrGTYd2xeB4RFJhUw9X0Z53vU6OFB4nRDuNnrRSDDDA9mSuZXAasZPBybvBKokmnawOCDRE+5nUhRxI43DV3Llz8fDDD6OgoABz585N+lqpyRGHyAVnfbNmuelNHV6B/1m9B5UluVh+y1mWdUzG0AuPkxEJ4TGrtg4ky9NnqianxEZ2FY/DwcLpIzHv9W1QASiArQxO3toQEv2rWBlqkea2UpOTiLQNV23evBl+v1//70RITY5Y2Kk/4+VwOtc2J7fTYUuorW0YToeSciVSOxQRaNLJKvTjotCgM/M0OZGNWlVVU/sajxTyGROqsGlfI5Z9Wotrz+xvK4OTd9iERP8qVmuHZN83noJvmtgJ9dLCkpFj1OFITU76oPe2ManJUVWVS5aSVvUzaNMNzPtmq7ncrYQ3YmHV1kFLFSYhPM7Y02jHugmpQIsvGNWfLBFeTinkpQXhsebYvG4CnL/DopzU9GzG7u/s+r4R0ORkaPkFkZp0ZtbMSuJitbt3IKTqAmCWJ1K9+Z1NNzBvASyZtg5sTnZEhccMirDxIMfl0A1ms253XgZfqhor3o0iC1M0cqKSDpiFq+QBIRF2MhNpkVkzK4lLJFxlbaMGWIeryGzUvE+jqYSrWIUNSAqPM1WToyiKoRO5ubXDox0KEJl7n02j1cexQSeQuhfUeB2nY7gqU4XHIoSrMmtmJXGxqnSP7qTMMlyVmieHh1jaSEGKp1GAXeVZki73TN2oAWMncnPfKTcjp8MD4w/Y9YLyNVRT9YL6Ge5ZuseZ4NrJOE1OrgxXxWXx4sUYMGAAcnNzMWnSJKxfv97U77366qtQFAVXXHEF3QGmKaUWu3trG7WisE3FTvWEFGAUk09Eqi53gN2mR7StQ4YKjwGj292aF5R16C5Vz5yuyeEUrtI1OTY9Odq8OxRQTzrQvF1kMhMzW5PT7A0gRMBbnArCzOyyZcswd+5cLFiwAJs2bcLo0aMxbdo0HDlyJOnv7d27F3fddRfOOeccRiNNP/RigG0+Uw3wjOEHlplyqXoXeLZ0AAg16GRWDJCk8Lhjo+YU6qCJVbc7v3BVeO3YDVfphfS4ZVeldkBgGTJ1p6gdNMKzZxhNtL1QVVMrC0AC2xWPAaC6uhrV1dU4cuQIQjE3pqVLl1p6ryeeeAKzZ8/GrFmzAABLlizB22+/jaVLl2LevHlxfycYDOKaa67Bgw8+iI8//hiNjY22Pkemo1nV/qCKNn+wy94w2kadw3jhGVsNmE3ZNRLZMPiGq9r8QQSCIVsbF+sU8iCRDBG+XgCaWK334eXlyekwquyGqyL1mTjXybHtyWFnpGlzRLZOTmatnVy3EzkuB7yBEE60+fV7EA9sXxEPPvggLr74YlRXV6O+vh7Hjx+P+mcFn8+HjRs3YurUqZGBORyYOnUq1q5dm/D3HnroIfTs2RM33nhjl3/D6/XixIkTUf+yhXyPU19EZkJWrMqjx2L0BNhJI+ctgC3IiZTyb/HZEx+zrtrqJ1m1NcNOo0CkPL3Zm6+PUQmAWFLOruJsqKbqyQkw1IWRLKTJ2/tME1HSyG17cpYsWYKXXnoJ1113XcqDqK+vRzAYREVFRdTzFRUV+PLLL+P+zurVq/HHP/4RW7ZsMfU3Fi5ciAcffDDVoaYliqKgJM+N+mYfGlv96F2avNBeRMDLduEZPTCBkAqr7X9499/JcTnhcTngC4TQ7A3YOr2wSuV16pockuLJzNuorWZXcdPkuFILofD+DiOZiamGq+gbaZEsUNn3LRnFeW4cOenlnmFle2Z9Ph/OOusskmMxzcmTJ3HdddfhhRdeQFlZmanfmT9/PpqamvR/tbW1lEcpFgWesMXwzdHmLl/rC4a9ELzSYAF7mzXvjRowiI9TdLvTDhu4iZ5GM3ujBtIgu0pLIbebXcVbtJ+iJ4flNUiyTg5v7zNNrIr2aWHbk3PTTTfhlVdewf3335/yIMrKyuB0OnH48OGo5w8fPoxevXp1ev3u3buxd+9eXHrppfpzmibI5XJh586dGDx4cNTv5OTkICcnJ+WxpiPLNtRgX0MbAOC2v21Giy+QtPR7pAonW9e1MZMrlXAVL00OEDZyGlp8wgsodV2BDFclRd+oTZ5GtariXr/9Wkl2cBEKV/FaO6lrclgKj7XsqtQ9OQFGnlselFg8INDCtpHT3t6O559/Hu+//z5GjRoFtzvaNW+lQafH48G4ceNQXV2tp4GHQiFUV1djzpw5nV4/dOhQbNu2Leq5++67DydPnsTvf/979OvXz/oHylDqmtowf3lkrlQA9y7fjnOHlCfsDxXpv2MxXpQixtRPO6ckETwKqdbKYRWu0m5mJITHkTo5mbdRF1kIVy3bUIN/fRE+qD29chf6dMuz1UfKDh5C4Sp+hTTD8+wNhOALhCx7wvwMa2RFskClFzQZohQEtG3kbN26FWPGjAEAbN++PepndtKO586di5kzZ2L8+PGYOHEinnrqKbS0tOjZVtdffz369OmDhQsXIjc3FyNGjIj6/dLSUgDo9Hy2s6e+BbFrMaiq2FvfmtDI0TYMnz+IuqY2W80y7aAoCtxOBf6gaivezXujBlLXFrDSFenCY4LhqozU5GidyLs4jdo5TJAk5XAV5+8wSrTvDcDj8lj6fZbhtkh2lQxXJUOUgoC2jRzSDTpnzJiBo0eP4oEHHsChQ4cwZswYrFixQhcj19TUwJGBdThoM7CsAA4FUYaOU1EwoCw/4e98sDNcm2h3fQumLFqJhdNHMjuRuhwO+INBW1oREVy/hSnWymG16ZEUHmfDRt3VadTOYYIkKYerOH+HLqcDuS4H2gMhfFPfjHEF3S39vp9hsoSbUN83YyPkjFw7eWJocoSa2Tlz5mDfvn3wer1Yt24dJk2apP9s1apVeOmllxL+7ksvvYQ333yT/iDTjMqSPCycPhLabV8B8Mj0EQk33rqmNryyrkZ/HFLDJ9K6pjb6g0VqTTpFqB6qhavsFMAybnq0tRFuki53Ttl4LNBc7seavUnXgHaYMNLVYYIk5MJVfA4IyzbUoL3jOrpyyVos21DTxW9EwzLpINLWIbW1EzQ2Qs7AtaNpcnYfbWF2/4hHSjPb2NiI3/3ud7jppptw00034YknnkBTUxOpsUkIMWNCFaaP7QMAuOGsAUm9MnvqWxC7dLUTKQtSadKpGwgcPX6agHJH3QnLC5vlpkfS5c77BkmTVTuPAgCavUFMWbQy4c1XO0xoOJTkhwnS6OHHFMNVPAzVTqE+GwcrlkkHkbYOqXlyjGsvE/VsOw+dBACs3lWfdO3QxvYV/emnn2Lw4MF48skn0dDQgIaGBjz55JMYPHgwNm3aRHKMEgKU5odj3Dnu5GLigWUFiF1uLE+kqTTp9Af4h6v2Hw8bg/+7cb/lhR216TFqMkhGeMyneCRt6prasOjdHfrjrryaMyZUYWivIgDAY1eOYhbiBVLPluMZNkkW6jMLSyONlBfU2IIj07ygdU1tWL7pgP6YdUTAiO2ZvfPOO3HZZZdh7969WL58OZYvX449e/bge9/7Hu644w6CQ5SQwGUyjlxZkodLR/fWHzsVhemJNJXGkbxrfdQ1tWH11/X6Y6sL29jVmFUKOZkmg5kZrkrl5ltRnEtpVPFJ53AViVAfSyMt0vctVU9O5PdZNkJmAe+IgJGUPDn33HMPXK6IdtnlcuEXv/gFPv30UyKDk5DDY2FhjupbAgA4a3APrJ53AdsTaQrVRHmL+FJd2MZQA+2bTURXQLL/TmYZOXZuvtwqHjtTM1p1TwgHb1xsqE+xEepjmXTgTuEgZsR4rbBshMwC3ho1I7av6OLiYtTUdHbF19bWoqioKKVBSchjxZ2tbXiVJXnMPDgaEeGxDU8O54JmqYb6Ii53hfqm5yboyWEllmaNHZ2N7tVibCxohwP7Xcg7xs1JzzZjQhUmDugGALj/u8MsH6xYJh2Qauvg7yi6mmnrBoisHWfHPsY6ImDEdgr5jBkzcOONN+Lxxx/X2zusWbMGd999N66++mpiA5SQQTceTAgTAxxd19oma6fiMe+2DpUlebh8TG+8ueUgAOsLm6VwmmhBs0Dm1smZMaEKj63YiWMtPrx4wwScd1rPpK/nNRfuFLISjb/HUwCrFV7UMhStwHbtkBHtZ2qYV2PGhCqcO6Qce+tbMaAsn4uBA6Rg5Dz++ONQFAXXX389AoEAVFWFx+PBzTffjEWLFpEco4QAHgtZSzzDD+meQn7W4DK8ueUgxvYrxeJrx1pa2GybDJJv0Jmpm3Wexwm0RNLJkxGpFs6nC7mqhg8ITosaD15NeY1oc+ZLwYvLJFxFaO1k+roB+EQDYrFt5Hg8Hvz+97/HwoULsXv3bgDA4MGDkZ/PPuYm6RorxoMW0uLhRnWloC0ICLBpaCfhglyX5cXtZ3iDdFIRHmee2x0w3HxNeEF51Qwyhsf8wRCcDmstWbTDD09vXCpVm5kKj7UM0BS9oBHPc2auG1GwZOTMnTsXDz/8MAoKCjB37tykr7XSu0pCH7cF4XGAY9hHF/WlJDzmt2lE6s9YHz+rlg6AcaOWnpyu8FjQu3Dz5BiueV8whNwuSkXEIsJ3aGWPioVpWwdCDTojIcLMXDeiYMnI2bx5M/x+v/7ficg0pXgmEHGxmhce8/HkpCA8FiBclUpPKJYx+lSzcYwEMrROjkaOSU+OqqrcvFpGwbDV7zRcaZu/CFa7fuwUNOTT1iFVTw7//SobsGTkGPtVke5dJaGL28JplOepLpXMBRGEfKnUK2GZHWYUHquqmtLBRIR5p4nZcBXPMv0OhwKnQ0EwpFq+9lgWoUyGFraxkyHGMvSTauFFDRG8Z9mA7dmtqamBqsb/kuOllkv4YkXrwlV47LDvyeGZFaaRmq6AnS7CWHzMTiabERHChDSJCGKTf6fGa5aHV0ubf6vXnvFAIYQmJ4WkAxeTcBVZ4bHU5NDF9hUxcOBAHD16tNPzx44dw8CBA1MalIQ8bgvtEjRDiEcVzlQEsSK0FyDRe4ulrgBIPY2cpXHGA+1zebswHnxRxRx56NnsXXta2jvA+YCgh6vEXjukQr2+gAxXscD27CZycTc3NyM3l21Jc0nX6FoRExsgz/BDpC+M/XAVzwadKYknGXpEjAZsKuXpgyFV9wRl6mZtNlxl9EDwOCDoRoLF79M4bqup5ySxUpU9FpZeXFKifRmuYoPlFHItq0pRFNx///1RKePBYBDr1q3DmDFjiA1QQgYrLtaAofIua1IptCVC2EQvypZCuIrlaRRI7URqvCFlapaIxxXOVDJr5HhcfMr02w1XidJewEqqfiwsvbikPDkyu4oNlo0cLatKVVVs27YNHo9H/5nH48Ho0aNx1113kRuhhAhWTklchccO+/Fu3TgTIFxlq6AZ0zo5ChQlXDwulXBVlJGTodoCsynkekE6Tidzl91wlQCHA+Pft1cIlGWoN7Xq0hpSk8MGy0aOllU1a9Ys/P73v0dxcTHxQUnIY0l4HGIn4osllXYDvNs6AKmFq7SN+mS7H3VNbdQrhbocCvxBNaUePFGZORzDhDQx62FgaaTGw2MzXMWzOaeRVITHJ9vDpU2aO/6fJnaNyVhEqNCeDdiuePziiy8CAL744gvU1NTA5/NF/fyyyy5LbWQSorgtpGfyzFJKpeJxRJOTnrqCtbvrAQBbapswZdFKLJw+kmoHeJfDAX8wSCRc5XIocHCcd5qYrZPjDfD1iNgNlbLs+5QMu0basg01+M83DQCAR975EsV5bqrrhlhbBwFaaWQDto2cPXv24IorrsC2bdugKIqeTq7FdIPBIJkRSohgJY4sRLgqlYrHPMNVLnuFwuqa2vDm5oP645AK3Lt8O84dUk7No+NyKoA/Nbe7CD2PaGM+hZzvXNit3yJK2MRO+YW6pjbMX75Nf6yCxbqxrxs0ohuXMlxFFdur8ec//zkGDhyII0eOID8/H59//jk++ugjjB8/HqtWrSI4RAkJLPWu4uhGTUl4LEA3bKPLPVEdqXjsqW9B7KuDqoq99a0ERxeNNtZU6uRkw0btMXnz5V3CwG2zYrAIhwPA6AW1tm5iL1/a68aVwkHMiAjh9WzAtidn7dq1WLlyJcrKyuBwOOBwOHD22Wdj4cKF+PnPf5607YOEPVa0IjxvXKkU2tI2HRHaOgDhmL3Z0MXAsgIoQJSh41QUDCij1/DWmULhRY1Mr5EDRIwWs3VyeM2Fx6YgVpT2AnY0OQPLCuBQEGXo0F43xOrkCDLvmY7t2Q0GgygqKgIAlJWV4eDBsKu9f//+2LlzJ5nRSYhhpVCYnqXEIUZvt6AZYAyd8E8hB6zdbCpL8nDRsJ76Y6ei4JHpI6iKj1MJDWrwDtGwIF2Ex6mGq3h/h3Y0OZUleVg4faT+2KGA+rohnV3Fe94zHduzO2LECHz22WcAgEmTJuGxxx7DmjVr8NBDD2HQoEHEBighg9tCTyWeKaVOC5WZYxHhRGr821Yrt57WK3xouGREBVbPu4CqeBIgoy3QU3ddWRCu6uKa9HLWJ6UcruKuybFX52fGhCr0KQ0XoF1y7Vjq6yaVg5gRXXicwWtHBGyHq+677z60tLQAAB566CF873vfwznnnIMePXpg2bJlxAYoIYN+ygt23ZBRq+TJI4U8lQ6/IpyMjJldVlNhtc29f/cC6unjAJkePNmQIRLx5CRPpuBtLNgNV4nSXiCVzETN3qgsoRem0nClsEcZkZocNtg2cqZNm6b/9ymnnIIvv/wSDQ0N6NatG9eqmZL4eCxoRXgKeHXvgsUQiqqq+smK54lUURR4nA74giHhtRHaiTQV4bE255m8UVsPVzmpjykedsNVES0bZ0+Oy75n0c/Qo6iHq1IVHmd4OxRRsDW7fr8fF110Eb7++uuo57t37y4NHEExioi7OoEEdE8OB+Gxw94pKaooHfeiZvY+g4+xpkMPDaZg5Ghj9gZCqGtqIzIu0cgxmUIeER5zqpOTcrhKDE+OnbYOXoaib+1woKqpHRCaWsO15dr8stwKTWxdEW63G1u3biU9FglFjAZL1/U++J0w7KY1R7UX4FzUzG3yphgL69APiaJmH+w4AiCcyjtl0Uos21BDZGwiYT6FnLMmx6aezS9IuCqVise6gcnggOCymVxgZNmGGry97RAAYMmq3Rm5bkTB9hVx7bXX4o9//CPJsUgoYrzxd3VT46ktsJu5IFIPpYj+yaI2gvG8pyo8rmtqw1/+s09/rBUwzDSPjtlwlZfhjTYeemqzxQMC6+suER4LyRGxsNS3RCUX2BhrogKGmbZuRMG2JicQCGDp0qV4//33MW7cOBQUFET9/Iknnkh5cBJyOBwKnA4FwZDa5SYY4OnJsd1kUKu4HQnD8MJ+vRK2N0ktNGjX5Z6sgCEL4TQrzNbJ4Z3dp+lR7HYhF8WTYzXcFgiGdOExE0+Ow3zoPx7JChhm0roRBdtGzvbt2zF27FgAwFdffRX1M6nLERO3M2zkJNsEVVWN9IBKQ0+O2+ngfv25bfbg8TEWfEeaodpzufMoYMgD013IBamTY/X7DAhS0NFjU3hs/F5YzL3xEGVHfMyjgGE2Y9vI0bqRS9IHt8OBdoSSekmMp3qe2VXWhccdRo4ATSIjPXjELspmLCtgh8qSPFw+pjfe3BIuBMqigCEPzIaruFc8TtFIEMWTo7VEMXtYMdajYvEZFEWB26nAH1RteXK0Aob3vB4OWbEoYJjN2L4iampqEvbmqamRIioRMeNhMG6QXOrk2KzCK0r/HcCojbBXJ4fVZyAhPB7ZtxQAcNbgHkwKGPLAago57y7kdsNVvPuPRWtdzBsP3o5m0IoSHUqiie41s3lAmDGhCsMriwEAi6aPysh1Iwq2d9OBAwfi6NGjnZ4/duwYBg4cmNKgJHRwmci+MLpfWW0YRuyKYUUpaAYQ0OSwEh7brKtiRLuh9i7Ny9iTqOkU8jQNV4miyfHYFPQaPWisQtUkauVoq65XSS6BEUkSYfuqTuRObG5uRm6u/NJExG3CgDCK/rh0IbfpydFez1tXAKQermJWJ6djow6m4MlhmbrLC48zXNzPbLiKl7Ggh6ssXne6Jod3F3KXPSOHR5NYEk06vR0VtHMyeO2IgGVNzty5cwGE45L3338/8vMjYqlgMIh169ZhzJgxxAYoIYeZ8ISm13FwylKyWzJd2xQDoXBROp5eBbviadZdiSOhwVS6kItjXNIiXRp0um1ed8c7itJ5ORelczoUXZBrJeTGw9A24xXvCm3cOW4+FbKzBctGzubNmwGEPTnbtm2Dx+PRf+bxeDB69Gjcdddd5EYoIYaZUBBv17XdGjP/2n4YAHD4hBdTFq3EwukjucW5Ix4zq5qcYNTv04Zkg07eXgCaaJ8tEFIRCqlwJDD+eQuP7YQfl22owd8/3Q8A+PPafRjeu5irPsTtdMAbCFkqCMjDyLFbk8gIyyrN2YxlI0fLqpo1axb+8Ic/oKioiPigJHQwc/PlXuvDad27UNfUhhdWf6M/1orSnTuknItHx26jQea9qwgIj3nf2FlgvHn6giHkOuKfvLmvHYttHRIVpeO1boDwXHsDIUuGN4/sMBJrR/Oc5bgzd+2IgO0U8hdffBHV1dWorq7GkSNHEIrRUCxdujTlwUnI4jZRFyXAOUPETgr5nvoWxCb68SyuFUmFtanJYZ1CTqB3FW/RKk2M34c3EEJugvACb6+WVcG7iEXp7BwQuISrSHpBM3jtiIBtI+ehhx7Cgw8+iPHjx6OyspJ7ATZJ15gRxEYKAfJyuVsXHotWlE47UVs95TEXHusVj6XwOBlGgz+ZVoS38Niq0SpiUbrIHmXByOFgLNhNkNBQVVUPV0lPDl1sGznPPvssXnrpJVx33XUkxyOhiJmFybv6qZ2shcqSPMyY0A+vbqgFwL8onV0BaOQmyebAYLdbupFsMHIURYHH5YCvC60I9zo5FsNVsUXpFAGK0umtKaxkVzGuLwWknl3lD6q69znHJYXHNLF9Vfh8Ppx11lkkxyKhjDlNDt/CYHYzk8ZWdQMAnNGvlHtROrfNSsKsQz9EXO5ZYOQAQI4JD4OeLZMm4SogXJRu4oDw2lnwveHci9Lp4Sobnpwclp4cm/uUhtGIkynkdLE9uzfddBNeeeUVkmORUCbiYUiWXRX+GY9CgID9Bp1azYmK4lzuRensN0pkW6/EbnVpI6wLGPLCTBq5MJmJljVW4e+uZzH/+mZmannFEqkUzu4atLtPaRjT9aUmhy62w1Xt7e14/vnn8f7772PUqFFwu91RP5ddyMXDjKiX+0ZtM4QiUnzbTgp5MKTqfcNYbXrOFDdqgL/YlhVmjBwf56J6VsNVGiIVpdPnOWi+Zg8XTU6Knhxj+niikgQSMtg2crZu3aoX/du+fXvUz6QIWUzMZC5op3p+p1F75dJFylSwkyFifC0rbYGLRBqsvllntq7AzM2XdZ2jWOxqwfQDggDaEDvVwnlmV9nV5GRLmFcEZBfyLMLM6SNS64NvCrmqhr0bZqsui7Rp2HK5G42cNBIe8xbbskIzXL1Jw1W8a0zZDfWK4wVNJYWcaZ2cFEO9Xs76rWxCznAW4TIhiI0Ij/mGq4xjMYNXSCPHeoYIEIn304Zkg04R5p0mVjQ5vG5cdtKvja8X4Yarh9xseEHZenK61jcmQ6Q5z3RSmuGPP/4Y1157LSZPnowDBw4AAP7yl79g9erVRAYnIYvHZaJ3Fe8UcsMN3sqJ1CeSy91lPWxg9KCxitETrXic4Zu1KU0O9zo5dsNV4TCbCN+hJmC307uKpcEQCVfZ9eSIM+eZju0Zfv311zFt2jTk5eVh8+bN8Hq9AICmpiY88sgjxAYoIYeZvlCRYoB8U8gBIGgjw0KETcNOCjkPwTdJ4XGmn0g1oz9Z/RbeImxjjy0reP0CHRDshKt4tHVIsbmtSDqoTMf2VfHrX/8aS5YswQsvvBCVWTVlyhRs2rSJyOAkZNE3kCQLM8BZV2BMXbciPhbJ/at5SKwUNPNy8AKQEB77OXsvWGHNk8Op/IKNGjOAWPqQiMDbup6NbXZVajWmfALpoDId2zO8c+dOnHvuuZ2eLykpQWNjYypjklDCTHiCt5BUUZRIZWbBN7pE2Enl5eHJsdMMNRbe3gtW5FjQ5PDrQm7duA6FVKG8cbY8OVy6kKd2QNDDVQLsV5mO7Rnu1asXdu3a1en51atXY9CgQSkNSkIHM1k/uvCYkfg1HnZqUOi1PgQ4GdkSHnO40ZgJX3aFsd5HJhPxMMSfq2BI1XtApVO4KqryboLGoyyx1buKhxc0RdG+SBltmY7tGZ49ezZuv/12rFu3Doqi4ODBg3j55Zdx11134eabbyY5RgkhrKWQ81t8dqqJ+gS62XpspPLy8KDZ8ZjFwrt4JCs8Xdx8jc/zTiE3FpbsCmNKvBhrx45on192lX1PjtTksML2VTFv3jz8+Mc/xkUXXYTm5mace+65uOmmm/DTn/4Ut912m633XLx4MQYMGIDc3FxMmjQJ69evT/ja5cuXY/z48SgtLUVBQQHGjBmDv/zlL3Y/TlZgxsMQ4ByuAuxtICKmkFs5jfLR5BAQHguk56CJdl0lqpMTXecofcovaB5QRRGj1lFXHrN48LgG7dYk0sgWD6gI2C4GqCgKfvnLX+Luu+/Grl270NzcjOHDh6OwsNDW+y1btgxz587FkiVLMGnSJDz11FOYNm0adu7ciZ49e3Z6fffu3fHLX/4SQ4cOhcfjwT//+U/MmjULPXv2xLRp0+x+rIzGTPE3zf3K82TutJGdJFJ2lZ1wGw8Pmm5M2ixoFgiGuIdoWNGV8DiqYjUnY8F4wzR7841kVjmEqFQfEU9bD7mxDVel2NbBL054PdNJeYY9Hg+GDx+OiRMn2jZwgHCvq9mzZ2PWrFkYPnw4lixZgvz8fCxdujTu688//3x8//vfx7BhwzB48GDcfvvtGDVqlKzRkwQtjpzslMS7CzlgFMTa8OQIcDLy2Mi88PMQT2phQbsZIobrKOONnI62FYnWjjFcystYMN7kzYreRQub6F5QK72rAux7hqXc1kEgsXemY3uGFy5cGNcAWbp0KR599FFL7+Xz+bBx40ZMnTo1MjCHA1OnTsXatWu7/H1VVVFdXZ0w4wsAvF4vTpw4EfUv29CyfpItTC1ExNNYsFNNNJKSyX+zTkV4nE5NBo2n7YzX5Jj05PA09pwOBZp9Zbb8gmjhRo/LvieH5dpJua2Dn//1ki3YnuHnnnsOQ4cO7fT86aefjiVLllh6r/r6egSDQVRUVEQ9X1FRgUOHDiX8vaamJhQWFsLj8eC73/0unn76aXzrW9+K+9qFCxeipKRE/9evXz9LY8wEzCxMzbDg6snRPQzWC4KJ4MmxUydHd7m72AuPzYpUY/EGI3oOV4Z3Uu7KyOFdI0fDat80kbISAbu9qzoao3Lw5NitkyOaBy2TsX1VHDp0CJWVlZ2eLy8vR11dXUqDMktRURG2bNmCDRs24De/+Q3mzp2LVatWxX3t/Pnz0dTUpP+rra1lMkaRiLiCu04h53ky1wwsKzdfkTQ5dvrvcEmDJVTQjGeIhhVd1cnhoQuJh24kWAxXiXA4ACJGoteGno2pJyfF7CrRPGiZjG3hcb9+/bBmzRoMHDgw6vk1a9agd+/elt6rrKwMTqcThw8fjnr+8OHD6NWrV8LfczgcOOWUUwAAY8aMwY4dO7Bw4UKcf/75nV6bk5ODnJwcS+PKNMxkLQlh5NioQaGfSAXYNDw24vU8NmpXii53HmPmRVdtHfS54Hz9WRWTi+ZRsFNIk0vvKk14bDu7SvauYkVKdXLuuOMOvPjii9i3bx/27duHpUuX4s4778Ts2bMtvZfH48G4ceNQXV2tPxcKhVBdXY3Jkyebfp9QKKT30JJ0xlwKeUe4imP4wc4pSaSTkb3sKi1cxSENNlVPjgBzThuz4SreBl+kfIHZ7CqxwlWpVDzmUn5B1skRHtuenLvvvhvHjh3DLbfcAp/PB1VVkZeXh3vuuQfz5s2z/H5z587FzJkzMX78eEycOBFPPfUUWlpaMGvWLADA9ddfjz59+mDhwoUAwhqb8ePHY/DgwfB6vXjnnXfwl7/8Bc8++6zdj5TxmKp4LEAKuZ0wikg33FSqtrK8STpTTIMVac5p01WdHBGEx4B1TYtIfauAyDjsNLfl09ZBHhBEJ6U6OY8++ijuv/9+7NixA3l5eTj11FNth4RmzJiBo0eP4oEHHsChQ4cwZswYrFixQhcj19TUwGFoNdDS0oJbbrkF+/fvR15eHoYOHYq//vWvmDFjht2PlPGYqngsgIDSThhFpB5KdlLII5oOdvPutqF9MqKl+Yow57TpKlzFw5sQD6vhKp9gHoVUCmkyTSE3UY4jGSKF1zMd20YOAFRXV6O6uhpHjhxBKGZRJapvk4w5c+Zgzpw5cX8WKyj+9a9/jV//+teW/0Y2Y0Yrom2OXNs6WAyjhEKqUPqQlFLIeWSI2DVyAvy9fqyIhKvi12/hYaTGw3K4SjCPgrsLYzIeXFqipOjJEc2DlsnYNnIefPBBPPTQQxg/fjwqKyszPrsiEzBTpdOnp5Dzz64yaySI12Qw0t1bVVVTa4NLF3JHihkiAqXt0yYd6uQA1g1s0TwKHjuZiRwK60XaOshQr+jYNnKWLFmCl156Cddddx3J8UgoomcuJFmYQvSuckSMBDOI1mTQKB72B1V4TNS+4aHJ0QzZkBr2hjksis2zaaPuqqeSKOEqq9XCRRPAum2I9iNrh91niBwYZZ0c0bG9In0+H8466yySY5FQxkwZf5FSyM16GKI7QPP3KBoNFbObNY/eVU6DUWO2Qq4RHlWaeZHThVZElLmwnl2lVQoX4zvsqtt7PCKZiSz1bGQ8OaLMeyZje4ZvuukmvPLKKyTHIqGMKeExh5ttLFbbOhhFxyKETY3p91ZDbjwyRAB74uOs9OQkLAYoSJ0cixlzmnhclHCV1YrNvPR4qWtyOuY9Cw4IvLEdrmpvb8fzzz+P999/H6NGjYLb7Y76+RNPPJHy4CRkMZVCLkSDTmunJL3WhyAbtdZDSFXNCyj9PGp9OIweJ2nkJMN8WwfOKeQuq2tHrO+wq7BgLLyaxOoFS1OtkyM9OdSxbeRs3boVY8aMAQBs37496mcinKYlnTET7w4IkKVkVZMjWkdfRVHgdjrgC4Qse6OYZogYPE52xMdeQUI0LOjq5iuc8NhidpUo2hCrwmnj65iK9p3W9qhYREvdz2RsGzkffPAByXFIGGAmNVvTZvCseOwyMU4jolSbNeLRjByT2gIeN0mHQ4FDCQuP7WzWuvdJEOOSJtq1lagYoCieHG3dmvWECJddZVGT4+OUdJBqtXDRUvczGTnDWYQWgvIFQ1DV+IszEq7iWSfHWmqziGETq1kimlCUtaGmG5R2NDnZ6MnpUnjMuU6OFq5K04rHmnjYqpbN5VAsZwemgp3WLUZEC7FnMikVA2xsbMQf//hH7NixAwAwfPhw3HjjjSgpKSEyOAlZjDejYEiNq7sRI1xlrUidiKciqwJKXlltbocCH+yFq0Q0LmlhDFfFq30kSsVtq9W2I9lVYoRNjOM3U2NKC8uxnnd3CocDwBhiF2PeMxnbV8ann36KwYMH48knn0RDQwMaGhrw5JNPYvDgwdi0aRPJMUoIYfTOJFqcYgiP7XlyRNowrGoLfJxCP3b6hGmI4r1gQU5HDRY1QWgvXcNV+s1WEG9cbI2pruDVWsRqFpuRoDEjLAsOCLyx7cm58847cdlll+GFF16AyxV+m0AggJtuugl33HEHPvroI2KDlJDBqLPxBUPIjXN6S8cUcjE9OdY2QV4Gg9XicUay0ZMDhD937PoQRnish6vMrh2xupAbPci+YKjL+fRyMi5T0eQYQ54yXEUf20bOp59+GmXgAIDL5cIvfvELjB8/nsjgJGQxbgSJFiePPjCxOLVigGabDAqoDbHag4fXTdKpt3awvlmLaFzSItbIKYjpQyyKJ8dyF3ItXCXId2icP38gBHTR75lXzzqrjVCNSCOHLbZnuLi4GDU1NZ2er62tRVFRUUqDktDB2ZFNAyTeBAMCeHLcFm+8mohPpJutVU2Oj9O869WlUxIeixMmpIXToegGYTzDVZQGsVbDKKKlkBvn2cxn8HESTkfq5KgJkzgSoXnPHArfBI9swfYMz5gxAzfeeCOWLVuG2tpa1NbW4tVXX8VNN92Eq6++muQYJQRJphVRVTWSrcCzd5VlA0Gs0yhg6BNmOhU2vPGxd7vbb9IZSSHPfE0OkDy9WRThsX7dWQ1XibR2nOZ1Rbw8aKlUCxfNsMx0bIerHn/8cSiKguuvvx6BQAAA4Ha7cfPNN2PRokXEBighi9vpgDdBkTrjYuV5IrWqExFRG+Kx+Bl4aaFSER6LGCakicflQJs/GLdWjijhKruCd5HWjtvpQLs/ZKpWDq8wb2wShxV7RVY7ZottI8fj8eD3v/89Fi5ciN27dwMABg8ejPz8fGKDk5An2cndGLLg6Ua1XPFY0I0aiIShusLPyRsVmWv72gKRvAA0SVYrRwQtGxAJ9aZruAqwlgbPSxcW258uXhJHIjTvWbYcDnhjeZZXrlyJ4cOH48SJEwCA/Px8jBw5EiNHjoTf78fpp5+Ojz/+mPhAJWRIdnI3uofFqHicngXNAMM8W6x4zN6TY788vQgd61mih6vianLEMLSth6vE8yp4XOa9UTzaoYT/XtdJHIkQcc4zGcuz/NRTT2H27NkoLi7u9LOSkhL89Kc/lc05BSbZSc+4WLkKjy12+BWxTo7HYgp5JBWW7WatC49ldlWX5CTx5IjSWsRquErEyrtWMhP9+jXIdu1rTXiBSCscs4i4X2Uylq/szz77DJdccknCn1988cXYuHFjSoOS0EOvoxFnYWobo0OJpBbzwGrFY1FEn0bsNhpMJ+GxiGFCmiQLV/kEKe5mtT6TiOEqXXhswgvKUxfmtnlA8ApiEGcLlmf58OHDcLvdCX/ucrlw9OjRlAYloYdeETVOl2JRwg8uu72rBNo0rGtywq/jlgore1d1SaS1Q7DTz8QTHnf9fQZDqh6m5G2cGbFyQIgY2uwPZS6LHmcN3Xsmw1VMsDzLffr0wfbt2xP+fOvWraisrExpUBJ6RHquJK71IcpGbb5qq3geBSsbdTCk6plt3DQ5dlLINaNYoHmnSbIUclEOCHYMBECscFWOBU2On6OhrdckshquyrLDAW8sz/J3vvMd3H///Whvb+/0s7a2NixYsADf+973iAxOQp5km2CAk4gvFqfFzUPEWh/aydKM8WD8Lpj3rrKYyWZE1xZkyWatGdHxUsiFER5bKF2grRtArLWje0HjeJtj4XnAsdvaQbSmqJmO5RTy++67D8uXL8eQIUMwZ84cnHbaaQCAL7/8EosXL0YwGMQvf/lL4gOVkMGdpC9UpBCgIBu1ReEx7xuMES0MZCZcZRRYMhceE+jBI9K80ySpJkeQkKl+iLFgIDgdCvc1b8SK8JhnmNBlUf+kIWLx0kzGspFTUVGBTz75BDfffDPmz5+vl7RWFAXTpk3D4sWLUVFRQXygEjIku6kFhClNnxkFzQCTLnfDTZP13KfSoFOU8CYr0iGF3EoXctH6VmlYqRbOc97ttkQRsQ1NJmOrGGD//v3xzjvv4Pjx49i1axdUVcWpp56Kbt26kR6fhDDJsi+0Gx3Plg7Gv2+6GKCAMW6t1YG5jVozFhQoCtu5dxp68FhFRC0UTcx4cniHepNlT8YiYpgXsNZklOcBx25mooh1vTIZ2xWPAaBbt26YMGECqbFIGJDMw6DFwHkWAgQiYzTbE0Y7kYp0s7WyUfMUrUaaodpJIc+uE2nyFHIxhMceG+Eq0b4/Tc9mpRggF+GxzZYosk4OW8S6uiXUcSU5uWunP94btdVOypEYtzibhpUUcq8AuoLUupBnxzaSrBggrxIAsVhZOyLWyAEi6yCewDsWEbKrrIZ6pSeHLXKWsww96ydJMUDeRo7VrAUReyhFPoPguoIUhMd+QQrgsSKRJodnCYBYdD1LGoerrNT6ScvsKkHnPVORs5xl6Fk/SU6jvHUFrgzoQm6l8izP06jb5mnUeGPPFk9OohRyniUAYrETrhKtKJ2l3lUCeEEtZ1cJeCjLZOQsZxmRYoCJKx7zTidNFlKLh1dAbYiV0yjPDtZ2hcdGI1mkeadJIk2O0bPD2+CzcuONZFeJFa6yo2fj4smxm10l4KEsk5GznGUkywgQJYXcataCiCcjK7U+eGpytLkO2qzaGn4PceadJh5n2BiI9eQYjR7eXlBLFY8F1VRZ6l3F0WCw68kRVQuVqYh1dUuooy3MeILYSDFA3uEqTVeQzg06rYSr+GlbXEmKQyZDpBs7KxJ5cozhRtYlAGLxWNGzCNpDKdIjTGxDza6ezSdomDBTkbOcZSQTxAYEKe7mspjWLGJXXyu6Aj9PXYHucrfnyfG4+N/YWZHo5itKjRzAYrhKQA8oYLWQJr8Dgl09mx5eF2i/ymTkLGcZSXtXhcTYrDUjJ6QCIRPeHBFj3FbK6wuRBmv3NJpFG3XEkxPdhVykRqVGzZ1WjT4RooZNLLWm4OrJsecFFVXwnanIWc4ykvau6lh82umeF0bhc1epsKqqCllcS9+oTZzy9EJyLvbGpd2CZiLd2FmhGXSdhMcBMbRsQLQ3sKvvVNRU5mTtM2LRvWgcyy9ITY7YiHV1S6iTrC9UQJRaHwZPUlceBuNGLpInx0rYgGdzR9vCYwFDhLRJGK4SpL5UeAyRtdPVtSdipXDAmiZHiPILNj052bR2eCJnOcvQNpB4C9MviLbA6EnqKj3TuBGKdCK1Uq+EZ6NLvUKuTIPtkq6ExyJcf8ZrqKubr4iVwgFjuMpKdhVHL6hVTY6ggu9MRc5ylhG5qcUpBiikJyf5BhJVr0Wgk5El8STH0E8kQ8TaRs2ztg8vPAnCVTyF47EY+8515QnR6+QIdrO1kpkY8SiyN9QipS4s6tkENS4zFbGubgl1kmkwIsUA+d64FEWBU89cMKcrcDsVODg3FjWi1/oQPFyVqvDYk0UbdaKKx16OmqpYFEUxbSQIq8mxE67i2tbBqicn+7ygPJGznGV4khYDFCdWbLbRoKjakHQpyma1JpGGiK00aJPo5usX7Bo021NJVAGslVAvz/R9Xd9oM9QrmnGZqchZzjKSZQRo3h3enhzA/EYt6s02mfYpFj9HT0CqFY+zM4VcXOExYL7atqi6KkvVwrl6cqzV89LwCdiGJpORs5xlJOupJEoXcsB8k85M2Kh5zrvVPmEaPA0zXiTU5AhWcdtsuMonaLjKncCYjEVVVa5zn3KdHMHmPVORs5xlJNsARal4DBjDVeZc7qLcYDTMhtsAviE3p8Xq0hrZmAabkzBcJU6dHMBOuEqMcWuYNdLCBQ/D/81Hz2a9WriqqlJ4zBixrm4JdfSFmUx4LICAN9k4jYhYCBAwtnUwn0LO1eUuNTldkihc5RXIAwoY+9OZza4Sa+3k6GvHQmYl13CV+bXjDxoMsyxaOzyRs5xl6KckgVPIAYMruItTkridlMPjCYZUBE3W+uESrrLZZFCk0CYrEtbJEczgMyt6FzW7KllI3Yjx83EV7VtYO15DSxDR5j1TkbOcZSTbAEUpBhgeQ3oLj61UnuVZZ8Vuk0FR550m2o00EFKjeqqJZvB5LIarRPsOtXmMTdWPRbsGFSUSdmWJy8baMRrI0shhg5zlLCOZ8RBp0Mn/sjC7gXgFzVSI7iFkzhvFJQ3WpifHJ6iegybGa8wYCuJZdTceZluKiKvJMeuJinhxFYVHZqL1tcN7zNmIWFe3hDrJ4vU+PYWc/2Vh9uYr6s3WSnl9nm0BnEkqYCdDtLRpFhiNHKOXgWf/pHiYNRJE1bOZ1eTwzmqz0p9OQ1TDMpORM51lJPXkCFSq320yhVxUI8fpUKB50Lu+2fDTQqVaml6UGzsLjJ/VGHYQTXhsVtMiviZHbD2eW8+uMr929P1KsFYamYyc6SwjWXqmSNoCsynkPs6nuWSYrZUjQp0cmV3VNYqiRGrlGL5TPYVckLkw39Yh/PNcwW64bkP9mVCS65L3NWjPk9MRXhdgj80W5ExnGcmKv/Hshh2L1XCViJuGx+SJWi9Nn1ZVW7PPyAHiZ1iJdDgALGhatBRywcJVxnWQLIzKP1xlX5MjWtp+JiPGqpQwQxNHxgsDac+J0dbBWsVj0TZqILJZm9YWcC0GKFPIzRDPyBHN4LMarhJl3BrGdZDsM3g5ZiUC9jITRQ2vZzJyprMM3ZMTJz1TpMqtTpPtBkRNgwUMnci7SIWNnEj5ZYhYFh5n6WYdr7WDeMLjrsMogWAIWiRItO8wKjMxydrh7cVNpU6OiPtVpiLUTC9evBgDBgxAbm4uJk2ahPXr1yd87QsvvIBzzjkH3bp1Q7du3TB16tSkr5eE0b0LcWLd2o1OhIrHbpPtBkQ7RRvRRd5dFgPkFyZ0pSo8FnDeaRLpRB4p6uYVSLAPmAtXGbPDRPOCOh2K7mFMpmfjWSkcMN9fz0gkRJhd64Ynwsz0smXLMHfuXCxYsACbNm3C6NGjMW3aNBw5ciTu61etWoWrr74aH3zwAdauXYt+/frh4osvxoEDBxiPPL1wJ+mppLd1EOBEGql4nJ7FAAGjJqcrQy18w+QtPFZV61ki2Rquikoh56ipioeZcJVx/CKuHTNeUN6eHLfJ1jNGZN8q9ghzdT/xxBOYPXs2Zs2aheHDh2PJkiXIz8/H0qVL477+5Zdfxi233IIxY8Zg6NCh+J//+R+EQiFUV1fHfb3X68WJEyei/mUj2gaoqujUbkBbrCK43SOivvQsTQ8YDLUuw1X8e1cBna+HZPgEulZYkinhKq9uWCtcqgV3hZkDAn/hsY3sKr+4h7JMRYiZ9vl82LhxI6ZOnao/53A4MHXqVKxdu9bUe7S2tsLv96N79+5xf75w4UKUlJTo//r160dk7OmGK0m7gYgnh/+m5zYpiBXZk2M1hZyn8BiwWu8jO7UFcYXHgoXu3CYOCPrNVhDDLBYzDW55r307zW1FPpRlKkLMdH19PYLBICoqKqKer6iowKFDh0y9xz333IPevXtHGUpG5s+fj6amJv1fbW1tyuNOR5K1GxAyhdxkc0sRN2urjQb5FANMfD0kg/cNhhcRTU6cOjmCXIMR47rrcJWoqczuOB6zWHhroVw2wlWy4jF7XLwHQIJFixbh1VdfxapVq5Cbmxv3NTk5OcjJyWE8MvFI1m7AL5CA0mz9FpGzfMxqciKpsPyaDALWwlX+LA1X5STx5IhwOADMhVFEXjdAfGMylkj3dz6GWiptHbLtcMATIYycsrIyOJ1OHD58OOr5w4cPo1evXkl/9/HHH8eiRYvw/vvvY9SoUTSHmRE4HQoUJazJiV2cAZE8OVoKeRc3XpE3DbfL3CbIU1tgDFdZSYXNWk9OHA8Dz2KO8fCYCVcJHjYxkyHG24trNnvSiMh1vTIVIa5wj8eDcePGRYmGNRHx5MmTE/7eY489hocffhgrVqzA+PHjWQw1I4jURoksTlVVIynkAnhyXJY9OeJtGubDVfy8IoqimO74bkQ0HQor4oarBAuZWgpXCbhuAJNGDufu764kmaqJEN2DlokI4ckBgLlz52LmzJkYP348Jk6ciKeeegotLS2YNWsWAOD6669Hnz59sHDhQgDAo48+igceeACvvPIKBgwYoGt3CgsLUVhYyO1zpANuhwIforN+giEVWgaxCJu16RozAnsUIi00Em+CwZCqh4l4edBcTgWBkGotFTbLU8jjC4/5Hw4Ac2EU0YvSeUykkPM2LpM1O06E6POeiQhj5MyYMQNHjx7FAw88gEOHDmHMmDFYsWKFLkauqamBwxG5MJ599ln4fD5ceeWVUe+zYMEC/OpXv2I59LTD7XIAvmDUyd1oTAhRJ8fkKckr2CnaiMdEuMr4M17hDrfDgXaErJ1IBZ53mmifN16dHI9TDK+IqXCV4EXpItlVZjw56VMMUGTPc6YijJEDAHPmzMGcOXPi/mzVqlVRj/fu3Ut/QBmK5mHwBSKGjdH9LkLFY5fJFHKvX9yTkZkMEeO88ytPH55rS3VyBPag0SS+J6fDEyeKJ0c/ICRJv9aK0gnWgVzDSsiNmwfU0HpGVVUoStffv8gawkxFznQWEq/5pdGYECEEEUkh70KTExT3RGpGk2MMGfLKajPbJ8xItmoL4mlyeFasjofbRGaSqB3INfS1YyZcxblODmD+gCBTyNkjZzoLiXfz1TYMhwIhKqBGdAXpq8kxI56M1CZSTJ0EaWC247sR2YXcKDwWK53eVDHANMmuSmao8daFGcP6ZjOsNINYVA9aJiJnOguJJ0wU7aYV6QtjLrtKRCPHYyJDjPdGDZg3KDVCIVXf1EWcd5rkxEshFyzTzG3i+xQ9bJJjRpPD2Yvriiq/YO6AoM+7IPtsNiBnOgvxxMkKEKnaMWAU9ZmreCziidRlQlcgQiE5zaA063KP0hEJOO80ifXkiJAdF4sZL4joYRMzDTr5h6sSF1ZNhB4mFLTSdCYi5hUuoUo8T05Av9nyD1UBxgadXQmPxdUWmAtX8T9Rm61JpOEVQEfEi1hNjl9Ag89cuErcdQOY07Px9oJqhVUB6DXGukLkQ1mmImc6C4l389U2ExHSxwFDg06TwmNRbjBGPCbqlfCu9QEYhMcmPTlRN3ZBrhdWxFY8NnpLRDH4zIWrBNfkxNE+xSJC6MdtsX+VrJPDHjnTWYg7TjaNCDdbIy4TJzljqECUcRuxUrWV5w3SbJ8wDZ/h5sJLLM0LrU+SdoM1Zv+Icg2aue4iYRMxxhyLmb5vInlBTWtyBK9PlInImc5C4hWxCgjU0gEwl/FjPOWJeDKKnEbF1uSYqatiRGSxN206h6v4Z8fFYsrIEayAYSxWigFy1bOZbN2iEQlXiTnvmUj27VKSuAtTuxGLUAgQiC60lQjN9QuIeTKykkLO9zRqTXgsUrd61kSEx8GO/+d/o43FUrhKUE+OLjwWOLsKsF5+QXpy2CNnOgtxxxMeh8TarJ0O8+nXDkUcLZERj4kN0C/ATdJqg07R049pkkiTI9JcWBMeizNuI2aqhfsDAhwQbGpyRJ33TETOdBYSbxMUrk6OiRRy0W+2egp5knCVCFooM/onIyLe2Fmh3ZwaWn2oa2oT1JNjonSB4NlVpsJVAuxZVjU5os97JiLOypQwI179FqO2QATM3HhFbxJpSnisbdQc+x65TXjNjIh4Y2fFJ7vrAQB761sxZdFKvLOtDoBY12A8T20sontyPGbWvwCHHP3AaLGtQzYeEHghZzoLiZdNo22IooR9zNx4tfi2R9BTkZmbjU+ANFj9NGpRkyPSjZ0FdU1t+J/Ve/THIRX4f6t2ARDrpmUqXOUXXZPTdUFDEQ45EdF+1weEoKFSuKjGZSYiZzoLiaSQG4sBipWKbUYMK4LwMBnm0mD5V8vV59qyy13MeafFnvoWqDGXo3Z5iuIBBcxl/IhQYyYZZjQ5EU8Ov7k3W7QUiP4sohqXmYic6SwkXq+iiCdHjM064l3oepMT9WZrRhuha6G4iifNtdDQECFMwIOBZQWITT7UHos0Fy5DZpIaa5V1oIerBG0vYCWFnGcavJXsKmM2qKjGZSYiZzoLicSRO1c8FkVnYaaSqOg3W81wMZMhxtflbk94LMq1worKkjwsnD4yytD5/hl9AIg1F8ZrKZEn1Cd4lo+ZUK8QxQAt1JjSDEunQxFGFpANyJnOQuLV0YikkAvmyTFR60NYI8dEvF4EXYHtiseCzjtNZkyowup7LsDAsgIAwM7DJwGIZeQYx5Jo/aSN8DhBZmLIoG/huWdpNuTRk94uXyu65zlTkbOdhcTL+tEWoHaq542ViseibhpuV9ceEr8A2VVmO75riGCY8aR3aT5uveAUAMD2AycAiHUNGkPOiYS76dKg05tg/MbPxcvYXrahBltqGwEA97+5Hcs21CR9veiHskxFznYW4opj5ERORWJcEmaKbIler8VUQTMBQj/6XJvtpKylkAs67yz43qhKdC/w6I9FWTdAJNQLJPbOadlVoq4dXZOTYO3wNnLqmtowf/k2/bEK4N7l21HX1Jbwd0T3nmUqcrazEL0Sr1F4LECjSCNOE6Ee8TNE0iSF3NH5ekiG9nlyBJ13FuS6nfjRhH7645Pt/qQ3OJY4HIph/YS/07qmNnyyu14fo+g33K5qTBkPDm4O3uc99S2IdXwGVRV761sT/o7o3rNMRcwrXEKVSDFAg/BYME+OmSJbomtDrKSQi9C7SjbotMa1Z/aH1pNzw97jmLJoZZchC1YYDexlG2owZdFK/PiFdZiyaCVeXV9jyK4S8zvU0sK7MnLcTgUODv32EmXaDSjLT/g7kbpeYs55piJnOwuJd3IPCJpCHgypXafBCnoyMlOvRIRMJatNBqWRE0ZREFU3J6R2HbJghebJ2bq/EfOWb9O9DiEVuPeNSJhF1LWjpYUnCvXyLkipZdoZt8urJ1ahsiQv4e+IXtcrU5GznYVoN6foFHL+YRMjRhd0IiNB9Jut20StDzEadForTe8TrNwAL/bUt3R6rquQBQuWbahBizesubn1lc0JCxgC4t5wNSF+ohpTIujCZkyowup5F+Ky0b0BAMeafUlfL7oOKlORs52FaDc1Y+NIzZAQzZMDJPYwCG/kWEgh56mFcskUclvEC1k4FSVpyII2sYLYeBjHLMqhJpauNDmi6PEqS/Jw8/mDAQArdx7BiXZ/wteKroPKVORsZyHxwhMiZPkYiTZyEpzmgh0nI0HGHIs2lyE1cVE2vwAubKvCY9HnnRWRkEV4/pyKgkemj0gasqBNPEEsABhtsV9cMhRA+PvjoWcxgydOZqJRPC1CIUCNob2KcErPQvgCIbz3+eGEr9M+S6svKERIM1tw8R6AhD3xTkkBwUIQ0WmwCTQ5frHFk0ZXuj8YgtPRWf+gedNE6F1ltkGn9OREmDGhCucOKcfe+lYMKMvnauAAEe+S8at0Kgqeu24sbvrzRgDAiN4lAMT2KGjXljcQNgg++uoo5ndoixwKMPucQeHXCbBfKYqCS0f1xpPvf4W3th7ED8b11X9W19SGPfUtGFhWgDW7wh3st+5vwpRFK7Fw+kjMmFDFa9hZgzRyspB4gli/AGETIw6Hom/WicIoPsFTmd0xRdly4/QJEsGDpo3zcFMb6praurxR+wVr5sqbypI87saNhuZdunf5dgRVVfcuTR3eC2cO6o7/fNOA1zbWAhD3cAAAb2+tAxBe/5MXroz6WUgFXvj4GwDiGNqXjq7Ek+9/hdVf16OhxYfuBR4s21CjG2axaCL1c4eUC3PtZCrSyMlCXHHqt2ineFEqHgNhD4MvEEroYRDdoxAlnk5U1EwAAeXmmuMAgPUdadBdnTBFn/dsJ5F36ftn9MF/vmnAu9sPARA3s6quqQ2/fvuLpK/RtgRRrsFB5YUY0acY2w+cwOIPvsZlo3snNHA0NJG6NHLoIsYVImGKOw2KAQIR4W5CT47gN1uHQ+mywzfvrLa6pja8ufmg/thMGrRX8HmXhD06kwf3iLqBXjKiEp6OgwMg7veXSFdkRJMSiRJeB4CqbmHB+R9X78Xliz/p8jPwFqlnC+JcIRJmxNXkhPiHTWKJrdoaiygZFsnQvGZd1vvg1LtqT30LYme3qzRoEUJsEuuU5Llx0bCe+mNRNTnxstYUQC+8CAA/HB+uNi3K2q9rasOKzw8lfY2CiHEmgkg9W5DhqixEC0lFhav0FHIxNg3AWPW4iyaDcbQuouB2OtDuDyWu3MpZ8J1IqJrshCm6J0CSmCvO6KOHq/yBkCkNFmsS6YomDeyB7/z+I7T6Q2j1hTP8ROmflsj7pK0t7TOIJFLPFqSRk4VoXgNjCEU04TFgrN+SKIVcfE+OJ47I24ivozMxLyNHu6HMW75NLxrX1QkzHeZdEp/zTytHrsuB9kAIu+tbhM3ySaQr+uGEKrz0yV68uz0sTBblGkx0WFh+y2S0+kJRn0EaN2wR4wqRMEX35BhCKKKlkAPxPU5GNANBZI9CV0XNROhdNWNCFd6+7WxoX/3IPqVJXx/x5IhjEEvM0dDiQ7th3YvUiiKWeLqiH00Mh6m0dSNKyC1RzaTR/bp1+gwStkhPThbijlMXRYQeSrFEihZ2ockRZKOLR6Q8vZg9eDSG9y7BxcN74d3th/Daxlos6H26/jNjrY/KkjzDmMUNE0rik6wVRTrciIf2KsbYqlJsqmkEADS2+oQJuYlWM0kSRho5WYg7Tgq5pnsRpa0DENEHJQxXpUGZdN2g7EJ4LIJx+cMJ/fDu9kN4c/MBzPv2UOS4nFG1PhwKsHD6SKnJSWPsaLDiEQwG4fcnbmFAk+sn9sbh4ycBAHuPNOKqxR9i7reG4Dsje3MZj5FuOQq69SkAALS3t3MeTXrjdrvhJHCQkkZOFhLPePAHxCvwFkm/Ts8UciAyn5o3KtYr4hUodf/cU8vRqzgXh060o3rHEZxRVRpV60MLbZQX5QAQe94l8Ukk6rXidWhubsb+/fuhxnb+ZMTg3BB+dUHPqOcUfyN27W7XMzIl6Y+iKOjbty8KCwtTeh9p5GQh2g3VF1UMsMOTI9Am4e7KkyNA36eucBnmetmGGl3gq3lFROrB43Qo+MG4Plj8wW4s+XA3pg7r2SljJKiqaPdrYmlxrhWJeVIJqwSDQezfvx/5+fkoLy+HorC/Blq8foSOd9YQVXbLR0GOvKVlAqqq4ujRo9i/fz9OPfXUlDw68orIQiLGQ+cGnSKlkMerzGxE9+QIrA3R5vpQYxvufWO7XpNG84qoEMuDlu8Jbwlb9zdh6/6mTj8PhzrEEn1KrGO3FYXf74eqqigvL0deHh/NicPlgeNkMKq+kwIFRQX5QhwWJGQoLy/H3r174ff7UzJy5BWRhcTrjh0QsB+R1hahK+GxyD14NI/+oyt2xi26p300ETQ5dU1t+N2/d3Z63ujc+8G4vvr3IbJxKaELDw+OhsflQJ9ueVA6eqsrUNCnW640cDIMUteY9ORkIUZxsdYdO1IMUJwQhHlPjpib27INNdhS2wgAaGzrLNI0CkBFKGqWqKDZH350BtbsqsffNtRi477jhn5b4lwrkuyie0EOCnPc8AVC8Lgc0sCRJEReGVmIsXGkdioXKctHw+lIXgzQ26ENaWz1MRuTWeqa2jB/+bao54xl3QFgVN8S/b9FMNTildN3KgrGDeiGe787DMW5LnxztMXgyeE/Zkn24nE5UJjrkgaOJCny6shCjIJRLbU5IGDF42RtHV5dX6PX+Zn+7CdYtqGG6di6Ip5XREXYK3Lfd4YBALYdOKH/TIR5T1TQrLIkD0W5btwwZWDU6xtaxDMuJRKJNW644QZcccUVvIdBDWnkZCHGNEstq8ovZMXj6AaddU1t+GR3PTbta4jykohYtTWZV+TGcwZieGWxrodyOxWuGgcjMyZUYfW8C/C32Wdi9bwLosr9zzprADwGY2zaUx8JZ1xKJIk4dOgQbrvtNgwaNAg5OTno168fLr30UlRXV+uvGTBgAJ566qm4v//ZZ5/h6quvRr9+/ZCXl4dhw4bh97//vem/X1NTg7vuugujR49GWVkZBg0ahCuvvBIrVqxI9aNxZ/ny5bj44ovRo0cPKIqCLVu2dHrN888/j/PPPx/FxcVQFAWNjY1MxibOHU3CDEVRonoqqaoaSSEXwKOgYcwCW7ahBlMWrcSPX1iH6c+utdw5mzXJvCKKomD2uRGviENRhDLQ4pXTB4D2QDCqB5eIxqUkfdAOLSyun71792LcuHFYuXIlfvvb32Lbtm1YsWIFLrjgAtx6662m3mPjxo3o2bMn/vrXv+Lzzz/HL3/5S8yfPx/PPPNMl7/7l7/8BSNGjMCBAwfwq1/9CtXV1fjb3/6GM888Ez/5yU9w/fXXIxgMpvoxudHS0oKzzz4bjz76aMLXtLa24pJLLsG9997LcGRSeJy1OBwAgsDBxlZUFOXoWUAi6Sw0g6uhxYdn/rkrrihWw07VVtokq0fyvVG9seAfn+NEewDeQEjYRolG9tS3JDQuZQn77ERVVbT5rd+cX9+4Hwv+73O9kvaDl52OH4zra+k98txO0x7QW265BYqiYP369SgoKNCfP/300/Ff//Vfpt4j9nWDBg3C2rVrsXz5csyZMyfh77311lu4++678e9//xtnnnlm1M8mTZqEm2++GT/4wQ9wxx134Omnn074PgMGDMBPfvIT7Nq1C6+99hq6deuG++67Dz/5yU/012zbtg2333471q5di/z8fPzgBz/AE088oRfUCwaDuPvuu7F06VI4nU7ceOONnYo6hkIhPProo3j++edx6NAhDBkyBPfffz+uvPLKhGO77rrrAISNyUTccccdAIBVq1YlfA0NpJGThSzbUIN2f9hzM+O5/+Chy0foPxOqTk6HQPposzeugaNlJ9mp2sqKRPVI6pu9ONke0B9rXpFzh5QL+TkAci0BJJlDmz+I4Q/8K6X3CKnA/f/4HPf/43NLv/fFQ9P0uk7JaGhowIoVK/Cb3/wmysDRKC0ttfR3jTQ1NaF79+4Jf+7z+TBnzhy89NJLOPPMM7F69WrccccdqK2txfe//320trZi2rRpePnllzFkyBDccccdGDx4cML3+93vfoeHH34Y9957L/73f/8XN998M8477zycdtppaGlpwbRp0zB58mRs2LABR44cwU033aT/fe33X3rpJSxduhTDhg3D7373O7zxxhu48MIL9b+xcOFC/PWvf8WSJUtw6qmn4qOPPsK1116L8vJynHfeebbnihfi3NEkTIjN+gmpwAP/2K4/Fqnisb/Dfbt297FOP3MqCt645ay42pF0IJlXRFSSheAkElHZtWsXVFXF0KFDib7vJ598gmXLlkV5UmL58MMPUV5ejksuuQSNjY24/PLL8d3vfhf/+te/UFZWhldeeQV+vx89evTAd77zHbz33ntJ/+Z3vvMd3HLLLTjllFNwzz33oKysDB988AEA4JVXXkF7ezv+/Oc/Y8SIEbjwwgvxzDPP4C9/+QsOHz4MAHjqqacwf/58TJ8+HcOGDcOSJUtQUhLJ8vR6vXjkkUewdOlSTJs2DYMGDcINN9yAa6+9Fs899xyBWWOP9ORkGfGyfoyP65u96NuN/8l82YYa/N9ndQCAvcfCN34F4Qwl7eY6ul83fgNMkXT1ishOyxIjeW4nvnhomqXfOdTUjqlPfBh17TsU4P2556FXSa6lv20GGj22tm/fjssvvxwLFizAxRdfnPB127Ztw1lnnQUgbBT16NEDDz74IABgzJgxWLZsmf7ayspKHD9+POnfHTVqlP7fiqKgV69eOHLkCABgx44dGD16dJS3asqUKQiFQti5cydyc3NRV1eHSZMm6T93uVwYP368Pke7du1Ca2srvvWtb0X9XZ/PhzPOOCPp2ERFGjlZRrybq5FzH/uAuzYkUY2ZN26ZjDa/mhE3VxKNEnlhtyWAJPNQFMVUyMjIoPLCuNf+oPLUGjEm4tRTT4WiKPjyyy+JvN8XX3yBiy66CD/5yU9w3333JX1tIBDQ21/4fL5O4TJj88lNmzbhpz/9adL3c7vdUY8VRUEoQQNjOzQ3NwMA3n77bfTp0yfqZzk5OcT+DkukkZNlxN5cYxFBG5KoxkybX8XkwT24jIkG0isiyVZYXvvdu3fHtGnTsHjxYvz85z/vZGg0Njaa1uV8/vnnuPDCCzFz5kz85je/6fL1p5xyii60nTBhAr788kv84x//wKWXXoq33noLn332Gdra2vDb3/4WtbW1uOyyy6x+PJ1hw4bhpZdeQktLi/4Z16xZA4fDgdNOOw0lJSWorKzEunXrcO655wIIG2EbN27E2LFjAQDDhw9HTk4Oampq0lJ/Ew9p5GQhxg3mWIsXc17ZHPVz3hkz6RrKsYP0ikiyFZbX/uLFizFlyhRMnDgRDz30EEaNGoVAIID33nsPzz77LHbs2KG/9sCBA53qvPTv3x8HDhzAhRdeiGnTpmHu3Lk4dOgQAMDpdKK8vDzu3506dSpmz56Nr776CkOGDMHixYtx9dVXw+fzYcKECZg2bRpuv/12fPvb30Z1dXVK3pJrrrkGCxYswMyZM/GrX/0KR48exW233YbrrrsOFRUVAIDbb78dixYtwqmnnoqhQ4fiiSeeiKpXU1RUhLvuugt33nknQqEQzj77bDQ1NWHNmjUoLi7GzJkz4/7thoYG1NTU4ODBgwCAnTvDPfB69eqFXr16AQjXKTp06BB27doFIBzKKyoqQlVVVVLxdsqoWUpTU5MKQG1qauI9FK4cbGxVB877p9r/nsi/QfPeVg82tnId16vr96mD5r2tj+fV9fu4jkcikahqW1ub+sUXX6htbW28h2KZgwcPqrfeeqvav39/1ePxqH369FEvu+wy9YMPPtBf079/fxVhx3HUv7/85S/qggUL4v6sf//+Sf/uo48+qo4ePVqtr69XVVVVvV6vevDgQVVVVbW+vl5tbTW31/bv31998skno54bPXq0umDBAv3x1q1b1QsuuEDNzc1Vu3fvrs6ePVs9efKk/nO/36/efvvtanFxsVpaWqrOnTtXvf7669XLL79cf00oFFKfeuop9bTTTlPdbrdaXl6uTps2Tf3www8Tju3FF1+MOzfGsSWavxdffDHueya71qzcv4Uxcp555hm1f//+ak5Ojjpx4kR13bp1CV+7fft2dfr06foFGfvFm0EaORFENSgONraqn+yq525wSSSSMOls5PAiFAqpP/vZz9S+ffuqzz//vHrkyBFVVVW1ublZ/d///V919OjR6oYNGziPUjxIGTlChKuWLVuGuXPnYsmSJZg0aRKeeuopTJs2DTt37kTPnj07vb61tRWDBg3CVVddhTvvvJPDiDMLUbUhMpQjkUjSHUVR8Oyzz+Lb3/42HnvsMfzsZz+Dy+VCIBDA+PHjcd9992H8+PG8h5mxKKpKIb/OIpMmTcKECRP08tihUAj9+vXDbbfdhnnz5iX93QEDBuCOO+7Qqykmwuv1wuv16o9PnDiBfv36oampCcXFxSl/BolEIsl02tvbsWfPHgwcOBC5uebTvSUR2traUF9fj9LSUhQVFfEejrAku9ZOnDiBkpISU/dv7sUAfT4fNm7ciKlTp+rPORwOTJ06FWvXriX2dxYuXIiSkhL9X79+/Yi9t0QikUgkZsjLy0O/fv2kgcMI7kZOfX09gsGgrv7WqKio0NXrJJg/fz6ampr0f7W1tcTeWyKRSCQSiXgIoclhQU5OTtoWM5JIJBKREEDlIMlwSF1j3D05ZWVlcDqdem8NjcOHD+v59RKJRCLhj9MZbqXg8/k4j0SS6WjXmHbN2YW7J8fj8WDcuHGorq7GFVdcASAsPK6urk7avl4ikUgkbHG5XMjPz8fRo0fhdrvhcHA/J0sykFAohKNHjyI/Px8uV2pmCncjBwDmzp2LmTNnYvz48Zg4cSKeeuoptLS0YNasWQCA66+/Hn369MHChQsBhC28L774Qv9vrUJlYWEhTjnlFG6fQyKRSDIZRVFQWVmJPXv2YN++fbyHI8lgHA4HqqqqoChKSu8jhJEzY8YMHD16FA888AAOHTqEMWPGYMWKFboYuaamJurEcPDgwaiOqI8//jgef/xxnHfeeXqfEIlEIpGQx+Px4NRTT5UhKwlVPB4PEU+hEHVyeGAlz14ikUgkEokYpFWdHIlEIpFIJBIaSCNHIpFIJBJJRiKNHIlEIpFIJBmJEMJjHmhSpBMnTnAeiUQikUgkErNo920zkuKsNXJOnjwJALKHlUQikUgkacjJkydRUlKS9DVZm10VCoVw8OBBFBUVpZyHH4vW4by2tlZmbnWBnCvzyLkyj5wr88i5soacL/PQmitVVXHy5En07t27yzTzrPXkOBwO9O3bl+rfKC4ulovAJHKuzCPnyjxyrswj58oacr7MQ2OuuvLgaEjhsUQikUgkkoxEGjkSiUQikUgyEmnkUCAnJwcLFixATk4O76EIj5wr88i5Mo+cK/PIubKGnC/ziDBXWSs8lkgkEolEktlIT45EIpFIJJKMRBo5EolEIpFIMhJp5EgkEolEIslIpJEjkUgkEokkI5FGDmEWL16MAQMGIDc3F5MmTcL69et5D4k7CxcuxIQJE1BUVISePXviiiuuwM6dO6Ne097ejltvvRU9evRAYWEhfvCDH+Dw4cOcRiwOixYtgqIouOOOO/Tn5FxFOHDgAK699lr06NEDeXl5GDlyJD799FP956qq4oEHHkBlZSXy8vIwdepUfP311xxHzI9gMIj7778fAwcORF5eHgYPHoyHH344qv9Pts7XRx99hEsvvRS9e/eGoih48803o35uZl4aGhpwzTXXoLi4GKWlpbjxxhvR3NzM8FOwIdlc+f1+3HPPPRg5ciQKCgrQu3dvXH/99Th48GDUe7CcK2nkEGTZsmWYO3cuFixYgE2bNmH06NGYNm0ajhw5wntoXPnwww9x66234j//+Q/ee+89+P1+XHzxxWhpadFfc+edd+Ktt97Ca6+9hg8//BAHDx7E9OnTOY6aPxs2bMBzzz2HUaNGRT0v5yrM8ePHMWXKFLjdbrz77rv44osv8Lvf/Q7dunXTX/PYY4/hD3/4A5YsWYJ169ahoKAA06ZNQ3t7O8eR8+HRRx/Fs88+i2eeeQY7duzAo48+isceewxPP/20/ppsna+WlhaMHj0aixcvjvtzM/NyzTXX4PPPP8d7772Hf/7zn/joo4/wk5/8hNVHYEayuWptbcWmTZtw//33Y9OmTVi+fDl27tyJyy67LOp1TOdKlRBj4sSJ6q233qo/DgaDau/evdWFCxdyHJV4HDlyRAWgfvjhh6qqqmpjY6PqdrvV1157TX/Njh07VADq2rVreQ2TKydPnlRPPfVU9b333lPPO+889fbbb1dVVc6VkXvuuUc9++yzE/48FAqpvXr1Un/729/qzzU2Nqo5OTnq3/72NxZDFIrvfve76n/9139FPTd9+nT1mmuuUVVVzpcGAPWNN97QH5uZly+++EIFoG7YsEF/zbvvvqsqiqIeOHCA2dhZEztX8Vi/fr0KQN23b5+qquznSnpyCOHz+bBx40ZMnTpVf87hcGDq1KlYu3Ytx5GJR1NTEwCge/fuAICNGzfC7/dHzd3QoUNRVVWVtXN366234rvf/W7UnAByroz83//9H8aPH4+rrroKPXv2xBlnnIEXXnhB//mePXtw6NChqLkqKSnBpEmTsm6uAOCss85CdXU1vvrqKwDAZ599htWrV+Pb3/42ADlfiTAzL2vXrkVpaSnGjx+vv2bq1KlwOBxYt24d8zGLRFNTExRFQWlpKQD2c5W1DTpJU19fj2AwiIqKiqjnKyoq8OWXX3IalXiEQiHccccdmDJlCkaMGAEAOHToEDwej74INCoqKnDo0CEOo+TLq6++ik2bNmHDhg2dfibnKsI333yDZ599FnPnzsW9996LDRs24Oc//zk8Hg9mzpypz0e8NZltcwUA8+bNw4kTJzB06FA4nU4Eg0H85je/wTXXXAMAcr4SYGZeDh06hJ49e0b93OVyoXv37lk9d+3t7bjnnntw9dVX6w06Wc+VNHIkTLn11luxfft2rF69mvdQhKS2tha333473nvvPeTm5vIejtCEQiGMHz8ejzzyCADgjDPOwPbt27FkyRLMnDmT8+jE4+9//ztefvllvPLKKzj99NOxZcsW3HHHHejdu7ecLwlx/H4/fvjDH0JVVTz77LPcxiHDVYQoKyuD0+nslOVy+PBh9OrVi9OoxGLOnDn45z//iQ8++AB9+/bVn+/Vqxd8Ph8aGxujXp+Nc7dx40YcOXIEY8eOhcvlgsvlwocffog//OEPcLlcqKiokHPVQWVlJYYPHx713LBhw1BTUwMA+nzINRnm7rvvxrx58/CjH/0II0eOxHXXXYc777wTCxcuBCDnKxFm5qVXr16dEkwCgQAaGhqycu40A2ffvn147733dC8OwH6upJFDCI/Hg3HjxqG6ulp/LhQKobq6GpMnT+Y4Mv6oqoo5c+bgjTfewMqVKzFw4MCon48bNw5utztq7nbu3Imampqsm7uLLroI27Ztw5YtW/R/48ePxzXXXKP/t5yrMFOmTOlUiuCrr75C//79AQADBw5Er169oubqxIkTWLduXdbNFRDOfHE4ord8p9OJUCgEQM5XIszMy+TJk9HY2IiNGzfqr1m5ciVCoRAmTZrEfMw80Qycr7/+Gu+//z569OgR9XPmc0VcypzFvPrqq2pOTo760ksvqV988YX6k5/8RC0tLVUPHTrEe2hcufnmm9WSkhJ11apVal1dnf6vtbVVf83PfvYztaqqSl25cqX66aefqpMnT1YnT57McdTiYMyuUlU5Vxrr169XXS6X+pvf/Eb9+uuv1ZdfflnNz89X//rXv+qvWbRokVpaWqr+4x//ULdu3apefvnl6sCBA9W2tjaOI+fDzJkz1T59+qj//Oc/1T179qjLly9Xy8rK1F/84hf6a7J1vk6ePKlu3rxZ3bx5swpAfeKJJ9TNmzfrGUFm5uWSSy5RzzjjDHXdunXq6tWr1VNPPVW9+uqreX0kaiSbK5/Pp1522WVq37591S1btkTt916vV38PlnMljRzCPP3002pVVZXq8XjUiRMnqv/5z394D4k7AOL+e/HFF/XXtLW1qbfccovarVs3NT8/X/3+97+v1tXV8Ru0QMQaOXKuIrz11lvqiBEj1JycHHXo0KHq888/H/XzUCik3n///WpFRYWak5OjXnTRRerOnTs5jZYvJ06cUG+//Xa1qqpKzc3NVQcNGqT+8pe/jLr5ZOt8ffDBB3H3qJkzZ6qqam5ejh07pl599dVqYWGhWlxcrM6aNUs9efIkh09Dl2RztWfPnoT7/QcffKC/B8u5UlTVUO5SIpFIJBKJJEOQmhyJRCKRSCQZiTRyJBKJRCKRZCTSyJFIJBKJRJKRSCNHIpFIJBJJRiKNHIlEIpFIJBmJNHIkEolEIpFkJNLIkUgkEolEkpFII0cikUgkEklGIo0ciSTLuOGGG3DFFVdw+/vXXXed3jmcN+effz7uuOMO3sMgzrx583DbbbfxHoZEwh1Z8VgiySAURUn68wULFuDOO++EqqooLS1lMygDn332GS688ELs27cPhYWF+vM7duzAPffcgw8//BCBQADDhw/H66+/jqqqKgBAe3s7/vu//xuvvvoqvF4vpk2bhv/3//4fKioqUhpPQ0MD3G43ioqKUnqfRNxwww1obGzEm2++SeX9E1FfX49BgwZhy5YtGDRoENO/LZGIhPTkSCQZRF1dnf7vqaeeQnFxcdRzd911F0pKSrgYOADw9NNP46qrrooycHbv3o2zzz4bQ4cOxapVq7B161bcf//9yM3N1V9z55134q233sJrr72GDz/8EAcPHsT06dNTHk/37t2pGTg8KSsrw7Rp0/Dss8/yHopEwhcqHbEkEgl3XnzxRbWkpKTT8zNnzlQvv/xy/fF5552nzpkzR7399tvV0tJStWfPnurzzz+vNjc3qzfccINaWFioDh48WH3nnXei3mfbtm3qJZdcohYUFKg9e/ZUr732WvXo0aMJxxMIBNSSkhL1n//8Z9TzM2bMUK+99tqEv9fY2Ki63W71tdde05/bsWOHCkBdu3atqqqRpoErVqxQx4wZo+bm5qoXXHCBevjwYfWdd95Rhw4dqhYVFalXX3212tLSEvXZjc1P+/fvr/7mN79RZ82apRYWFqr9+vVTn3vuuYRjU1VVfe2119QRI0aoubm5avfu3dWLLrpIbW5uVhcsWJCwSWFNTY161VVXqSUlJWq3bt3Uyy67TN2zZ4/+ntp39Ktf/UotKytTi4qK1J/+9KdRzTQT/V2NP/3pT2rfvn2Tjl0iyXSkJ0cikeBPf/oTysrKsH79etx22224+eabcdVVV+Gss87Cpk2bcPHFF+O6665Da2srAKCxsREXXnghzjjjDHz66adYsWIFDh8+jB/+8IcJ/8bWrVvR1NSE8ePH68+FQiG8/fbbGDJkCKZNm4aePXti0qRJUeGdjRs3wu/3Y+rUqfpzQ4cORVVVFdauXRv1N371q1/hmWeewSeffILa2lr88Ic/xFNPPYVXXnkFb7/9Nv7973/j6aefTjoXv/vd7zB+/Hhs3rwZt9xyC26++Wbs3Lkz7mvr6upw9dVX47/+67+wY8cOrFq1CtOnT4eqqrjrrrvwwx/+EJdcconuSTvrrLPg9/sxbdo0FBUV4eOPP8aaNWtQWFiISy65BD6fT3/v6upq/T3/9re/Yfny5XjwwQe7/LsaEydOxP79+7F3796kn1ciyWh4W1kSiYQOVjw5Z599tv44EAioBQUF6nXXXac/V1dXF+U5efjhh9WLL7446n1ra2tVAOrOnTvjjueNN95QnU6nGgqFOr1vfn6++sQTT6ibN29WFy5cqCqKoq5atUpVVVV9+eWXVY/H0+n9JkyYoP7iF79QVTXiyXn//ff1ny9cuFAFoO7evVt/7qc//ak6bdq0qM8e68kxepVCoZDas2dP9dlnn437mTZu3KgCUPfu3Rv357Fzraqq+pe//EU97bTToubB6/WqeXl56r/+9S/997p37x7ldXr22WfVwsJCNRgMdvl3VVVVm5qaVAD6PEok2YiLn3klkUhEYdSoUfp/O51O9OjRAyNHjtSf0wS+R44cARAWEH/wwQdR2hqN3bt3Y8iQIZ2eb2trQ05OTpQ4OhQKAQAuv/xy3HnnnQCAMWPG4JNPPsGSJUtw3nnn2f4cFRUVyM/PjxLeVlRUYP369abfQ1EU9OrVS//csYwePRoXXXQRRo4ciWnTpuHiiy/GlVdeiW7duiV8/88++wy7du3qpAVqb2/H7t27o947Pz9ffzx58mQ0NzejtrbW1N/Ny8sDAN37JpFkIzJcJZFI4Ha7ox4rihL1nGaYaEZJc3MzLr30UmzZsiXq39dff41zzz037t8oKytDa2trVEimrKwMLpcLw4cPj3rtsGHDUFNTAwDo1asXfD4fGhsbo15z+PBh9OrVK+HniP0M2nPaZ0iEld9xOp1477338O6772L48OF4+umncdppp2HPnj0J37+5uRnjxo3rNHdfffUVfvzjHycdm5W/29DQAAAoLy839Z4SSSYijRyJRGKZsWPH4vPPP8eAAQNwyimnRP0rKCiI+ztjxowBAHzxxRf6cx6PBxMmTOikefnqq6/Qv39/AMC4cePgdrtRXV2t/3znzp2oqanB5MmTCX8y6yiKgilTpuDBBx/E5s2b4fF48MYbbwAIf75gMBj1+rFjx+Lrr79Gz549O81dSUmJ/rrPPvsMbW1t+uP//Oc/KCwsRL9+/br8uwCwfft2uN1unH766TQ/vkQiNNLIkUgklrn11lvR0NCAq6++Ghs2bMDu3bvxr3/9C7Nmzep0U9coLy/H2LFjsXr16qjn7777bixbtgwvvPACdu3ahWeeeQZvvfUWbrnlFgBASUkJbrzxRsydOxcffPABNm7ciFmzZmHy5Mk488wzqX/WZKxbtw6PPPIIPv30U9TU1GD58uU4evQohg0bBgAYMGAAtm7dip07d6K+vh5+vx/XXHMNysrKcPnll+Pjjz/Gnj17sGrVKvz85z/H/v379ff2+Xy48cYb8cUXX+Cdd97BggULMGfOHDgcji7/LgB8/PHHOOecc/SwlUSSjUgjRyKRWKZ3795Ys2YNgsEgLr74YowcORJ33HEHSktL4XAk3lZuuukmvPzyy1HPff/738eSJUvw2GOPYeTIkfif//kfvP766zj77LP11zz55JP43ve+hx/84Ac499xz0atXLyxfvpza5zNLcXExPvroI3znO9/BkCFDcN999+F3v/sdvv3tbwMAZs+ejdNOOw3jx49HeXk51qxZg/z8fHz00UeoqqrC9OnTMWzYMNx4441ob29HcXGx/t4XXXQRTj31VJx77rmYMWMGLrvsMvzqV78y9XcB4NVXX8Xs2bOZzodEIhqy4rFEImFGW1sbTjvtNCxbtkyIUJOopFop+d1338V///d/Y+vWrXC5ZH6JJHuRnhyJRMKMvLw8/PnPf0Z9fT3voWQ0LS0tePHFF6WBI8l65AqQSCRMOf/883kPIeO58soreQ9BIhECGa6SSCQSiUSSkchwlUQikUgkkoxEGjkSiUQikUgyEmnkSCQSiUQiyUikkSORSCQSiSQjkUaORCKRSCSSjEQaORKJRCKRSDISaeRIJBKJRCLJSKSRI5FIJBKJJCP5/wgrzhaPUlZkAAAAAElFTkSuQmCC", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "scada_data.plot_bulk_species_node_concentration({\"CL2\": [scenario.sensor_config.nodes[10]]})" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "\n", "\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " Once\n", " \n", " Loop\n", " \n", " Reflect\n", " \n", " \n", "\n", "\n", "\n", "\n" ], "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from IPython.display import HTML\n", "\n", "scenario.place_bulk_species_node_sensors_everywhere(['CL2'])\n", "vis = ScenarioVisualizer(scenario, color_scheme=black_colors)\n", "vis.color_nodes(data=scada_data, parameter='bulk_species_concentration', \n", " species='CL2', use_sensor_data=True, statistic='time_step', \n", " pit=(0, 60))\n", "anim = vis.show_animation(return_animation=True, duration=10)\n", "HTML(anim.to_jshtml(fps=15))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Do not forget to close the simulation!" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "scenario.close()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.12.4" } }, "nbformat": 4, "nbformat_minor": 4 }