{ "cells": [ { "cell_type": "code", "execution_count": null, "id": "7a43bd85", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "TRITON_tseries.zarr already written. Not overwriting.\n" ] } ], "source": [ "import TRITON_SWMM_toolkit.examples as examples\n", "import TRITON_SWMM_toolkit.analysis as exp\n", "from TRITON_SWMM_toolkit.post_processing import TRITONSWMM_sim_post_processing\n", "import TRITON_SWMM_toolkit.post_processing as pp\n", "from TRITON_SWMM_toolkit.post_processing import *\n", "from importlib import reload\n", "reload(examples)\n", "reload(exp)\n", "reload(pp)\n", "norfolk_irene = examples.retrieve_norfolk_singlesim_testcase(download_if_exists=False)\n", "run = norfolk_irene.system.analysis._retreive_sim_run_object(0)\n", "\n", "# working on run methods\n", "self = TRITONSWMM_sim_post_processing(run)\n", "self._run.latest_sim_status()\n", "\n", "self.export_TRITON_outputs(verbose = True)" ] }, { "cell_type": "code", "execution_count": 4, "id": "f5957c12", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "SWMM_node_tseries.zarr and SWMM_link_tseries.zarr already written. Not overwriting.\n" ] } ], "source": [ "self.export_SWMM_outputs(verbose = True)" ] }, { "cell_type": "code", "execution_count": null, "id": "eedfae2e", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "5920cca0", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import TRITON_SWMM_toolkit.post_processing as pst\n", "from importlib import reload\n", "reload(pst)\n", "\n", "\n" ] }, { "cell_type": "markdown", "id": "08cf2481", "metadata": {}, "source": [ "# plotting" ] }, { "cell_type": "code", "execution_count": null, "id": "14e909ff", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAe4AAAGOCAYAAABR4+AMAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAb/hJREFUeJztnQeUE3XXxu8uvUqR3pUiTRDpTUBAQFCKNBEUEAEpIgIvRQSRoqh0QRCkdxBQRKqFaqF3BUVAemfpsMl3nqvhy2aTzSSZPvfHmbMkm0wmk+zc/23PjXK73W4SBEEQBMESRBt9AIIgCIIgKEcMtyAIgiBYCDHcgiAIgmAhxHALgiAIgoUQwy0IgiAIFkIMtyAIgiBYCDHcgiAIgmAhxHALgiAIgoVIrORBLpeLTp8+TWnSpKGoqCjtj0oQBEHQHehxxcTEUPbs2Sk6Wju/7s6dO3Tv3j1V9pU0aVJKnjw5OQlFhhtGO1euXNofjSAIgmA4J0+epJw5c2pmtPPlSU1nz8eqsr+sWbPSsWPHHGW8FRlueNqeDzNt2rRaH5MgCIJgANevX2cnzXPN1wJ42jDax3bkobRpIvPqr8e4KN/Tx3mfYrh98ITHYbTFcAtqgz+6H3/8kY4fP87biRMneEPIrnz58lStWjWqXr06Pfroo3LyBUEH9EiJwmhHaridiiKPWxC0ZNq0afTmm29yTg25tdy5c1OePHkoWbJktG7dOpo4cSKlTp2afvnlFypSpIh8GIJgA2LdLop1R74PJyLLHcFwUPyYJEkSzn0hHbNlyxaaN28eTZ8+nf7444+H+bY2bdrQ/fv3jT5cQRBUwEVuVbZQgSOQL18+Dq0//fTTtGnTpgQfP3fuXCpRogSlTJmSsmXLRm3btqVLly6RkYjhFgwnRYoUbJADVbHCaM+ePZv27NlD/fr148pXQRCEUFm4cCH16NGDBgwYQLt27aIqVapQ3bp1OTXnj82bN7PD0L59ezpw4AAtXryYfvvtN3r99dfJSMRwC4aDlSy4fft2wMeULl2aRo4cSZ9++im1bt2arl69quMRCoKgNi6V/oXCqFGj2AjD8BYuXJjGjBnDxXiTJk3y+/iff/6Z8ubNS927d2cvvXLlytSxY0favn07GYkYbsFwkMsOZrjB22+/TQsWLKCvvvqKsmTJQi+99BKdO3dOp6MUBEFNYt1uVTZPNbz3dvfuXfJXBLtjxw6qXbt2nPtxe+vWreSPihUr0j///EOrVq3iSB+uN0uWLKHnn3+ejEQMt2A4KDpDxXj69OmDPrZ58+Z09OhR+vDDDzkX3q5dO12OURAE85IrVy565JFHHm4jRoyI95iLFy9SbGwsL/q9we2zZ88GNNzIceO6A6EX9IynS5eOxo8fT0YiVeWCoWAVu2jRImrUqBElTqzs64jKc3jf58+fZw9cEATrEW5xmTee5/tqjHiieEpa3XANCtT+dvDgQQ6Tv/fee/Tcc8/RmTNnqHfv3tSpUyfuhjEKMdyCoezevZtVj5o2bRryc7FSRugqoT88QRDMCYxurEqGO60CjRFE9RIlShTPu4YD4OuFe4DnXqlSJTbW4Mknn6RUqVJxUdvQoUO5ytwIJFQuGIqnyAy926GCPzbkxW/cuBHRMdy6dYtq1qxJBQsWpFdffZVX0oFCZ4IgWLMdLGnSpNz+BW0Ib3AbIfFA1wbfbhcYf2Bkd4sYbsFQEPYGf//9d8jP9aySIy1QQ459w4YNVKhQIdq/fz+98cYblCNHDqpRowZNnjyZc2OCIFifnj170tSpU+nLL7+kQ4cOccoNrWAIfQO0m6L9y0ODBg24GBZV53/99RfX1SB0XrZs2YfXLiMQwy0YyrJlyzgfhZVwqKBQRA3D/euvv7Iy2/Lly7nqFKEzGGysrKHohtdBfgt/7FeuXInotQRBUL+qXCkoMkML2JAhQ6hkyZK0ceNGrhiHUiNADtu7p/u1117jFrIJEyZQsWLFOKWHBT6MuZFEuRX4+yivR6XetWvXRKtcUJVSpUrRE088wUppoYI/Mqx6ly5dSo0bNw77GPBcGOQffvgh3u9gxLF/CDfgjxwFdDDiuAC88MIL8vcg2Ao9rvWe1zh8KAuliVCrPCbGRU8UPuc42yQet2AoyDtBOS1UHjx4wEIIEG/BSjhSj7tcuXJ+f5c5c2bq3LkzD0FBP+cnn3xCly9fZhEY/A5GH0b95s2bER2DIAiCUsRwC4aCanBMAQsVGNDvvvuOxRBQVBYuyK2fOnWKypQpE/Sx8O6R30KeC1PMUFWKNpQWLVqwEYcXjta2SIvlBMEJoKJcjc2JiOEWDAOhZ0gKvvjiiyE9D4b+448/Zo8bOsORMHr0aBZ+8VVTCgaq4Hv16sW6xRCEeffdd/knjHemTJnYE4dwA0J4giDEB5PB1NiciBhuwRC+//57zhFXqFCBWrZsqfh5qASFzjC82r59+0aca5syZQq99dZblCZNmrD38/jjj3M1Kgrb/vzzT/rggw/o9OnT9Morr7AnXr9+fZoxYwaH2AVBECJFBFgEzYA28DfffMPV2jCSuI0pYPgJT7tatWoc6g40FcxfLrpq1apc2DJu3DieGhYJyJNjlCiqRBOiVnRTWudarGifjz32GHvi2BBGR/Up3iOkWVGljhYzaKw3bNiQPXNBcCoYDxLpNG0XORMx3IImwMjC07xw4QI99dRT3BeNQjLM3UZBWq1atdhjxm2lYG4utIaRX8Ys3UhBiBxGG8VnWFDUqVOHi8zgOSP3/eyzz3LLSCT6yfDmsaECHq1vMOLoGcX2zDPPsBFHeD1jxowRvx9BsBIuiqJYiop4H05E2sEE1XG5XCxQACOLPG+RIkVUC6/DmHbt2pVz00q1zRMC4Wv0aqMy3Bv0liMygFC62rN3sZhBFAJGHO8JEQfkxPE61atXVxyBEAQrt4PtPJiFUkfYDnYjxkWlikg7mCBEzJw5c9hrxQQdtYw2QJgZCkYTJ05UbTpPhgwZaP78+bR3717u10ZoH/lp5NChoAbP+MCBA6QmCJF36NCB1qxZw68FPWRotkN2tUCBAjR8+HC+XxDsjMutzuZExOMWVAcqaAiNf/3116rvG3lpCP3DwK1YsYK0BDN90SOeL18+NrJaDjKBDhJmAn/xxRfcUgZvH/rJhQsX5slpqHoXT1ywk8f9y4Gsqnjc5YqeFQEWQYiE33//nXbu3MlSgWqDHHfx4sX5NaAxrJWxxqzv999/nyVOYTwxhGDt2rWkJVgUYAoRqs+RD4fEIvrGf/rpJ255w0Jl5MiRHGYXBDuA/LYamxOR4jRBVRB2xkq9Xr16qp/ZPn360OHDhzkMD6lULXLzAwYMYC1jtHFB7hR5egDVNL2AN+IpYIMnjgr8zz//nGcC4/hg4CG7Ci8c50FGmgqCs5AqGBtw6dIlQ0fMecAxQHO8SZMmqlR9++Lx4uEJw8iqCarIUen+6aefcs4ZOWaEqzEZ7I8//qC2bduSEcAoo9d95syZrPA2duxYXhghD166dGkqX748y7GqfT4EQWtc7ihVNicihtvCoC0K1cgYEI8NXi5CvFAk0xvkxCAFeuTIEXr55Zc1eQ0opSEHjOK09u3bc1hbKWjz6tatG/Xu3ZsVznzz5jCOuB9hcTwGIKeMNi2Eqc2QX8axoAIetQNYrEHyFceOSvRs2bLxOMIFCxbwgkMQzI6EysPH+KuREBbwslC4hFnSyIdCQxseLzwy9AfDeOql1AUvG5XSVapU4RnZMCRagZYptJjNnj2bXwteuMeAJQRm8CLcPGvWLC5uQ1W6x0tFRfnZs2c5v4zKbiuAXnj0nUNyFd8FRAT27NnDKnRFixbl9yKV6YJgT6Sq3IJALATGB142Zsl6y3XCeCPP3KVLF+5FxiB4qHlBpQyh1VAET5SCyVrwSPFaGNEZyYhNpSDXjfcJDxOh7MGDB9OgQYMCPr5///5s8NHaBeGXzz77jBc4CLtDiAULD4TIPd62Vdm/fz+/h9WrV/NtfB7ofYfQCxZWUG8TBDNUlX+/P5cqVeU1ip10XFW5GG4LglApBnOgehuqZP5APnTYsGG0bds21s/GYI5UqVJxYROMOMLqMP4JFTZBHAQGERd7GABM4cI+8BwYavzEBkMIowBBESVAQhQolRFNCHiZUDeD1w0PNBCTJ0/mxQwkTiHcgvcGGdJz586xAUe7F+Zx49zYAbwveOJ4nzDiJ06c4HB6s2bNeJoZFltS1CYYabg37MtNqSI03DdjXPRs8RNiuBM60U5b1ZgVtEJBPhMFVUpAGBlGHhdybGirgsAIcrcYZwmD7L3Bgz948CB76DCK0ARHCxby1zB8/grhKleuzPvVGywWmjZtykVkCcmGQm0NHjpC4p6oAxYzyJdDNQ3va9++fbb0SPF5IaWC6ATeK85B3rx52YBjC7aAE5yDGG5rIB63BYGXXaJECc5jhgOKlzZs2MCeO0LHCDXDQ/OQNWtWNmp58uThHCo0xv0ZA+8NHriSAq5QBnYoza+3atUqwUWlx1Bh9CY2J4P2NhQvwohj0YM6CERTPEY82MAVwd7oabjX7sujisddu/hxxzmVYrgtCMK5+KJu375dNQ/x6tWr7FHDs4YhR+EWKpghAmJmEElA2BeRAxgjzMn2BZO6UI2OKvx06dIZcpxmrZVYv349nzdEcLBYw6hVLG6gfmeGSnrBvob7u735VDHcdZ88JoY7oRPttFWNWUHeGope8DZDmWVtVyBQgglbntA35mMj9I0NeWsUnQUrXnM6OFfwwCHwgnw4FmyYJ966dWsOoyMNgUiMYG/EcFsDUU6zIJi8hYpxKHsJRAMrfkr53CXpkRco3kIGRmfo0KHUr18/OVUJAMEcGGqEyzdv3szTyzDzHDKrAJEdKLlh8SNzxAU1wEhOV4QdyS4yXnjKCMRwWxDkKCE+EsmsaDvhyZkj146ebPxMkSIFGyNEiiQ8rhxU3CMVgw0Ss0gvABQeoksBE9QQxUBeXBAiQQ2t8ViHapVLjttiwChB5QtFRqgUljykvqhdXGclUOQHWVhU8KM7QQrZ7IeeofJlewpQqjSR1ejcjImlRiWOOC6NK9UnFgKynQgFw2AjhKmV0fb0WQv/fz48m1ONNkCOG90IUG3DBDVBEIxBDLeFQJ72m2++ocWLF2smK1orUXNN9mtlnGysfcHUNEjtohBQECLPcUe+ORHJcVuIX3/9lWrUqMHylWrhdC9SKXKO/h9UmP/1118GfhqCHUBhWqwUp4WFeNwWAoNDIF8JOVOtQuPrYheKkRISBCp50BCYMmUKtyYiZykIgn6I4bYQmISFaulJkyap6klik7y2oJT69etzuBwCPdATQKERBtkg7w0lPojioHhSEBIi1h2tyuZEnPmuLQqqJjG+EQMzIJihJhIKFpSSK1cubruD3v3u3btpzpw5PIEMPd7FihVj1TW0i8Ejh2QuNPVv374tJ1iIFypXY3MiznzXFqZbt2506dIlHpghCEaCPnlo5kMrHpKyx44d435vTCTD4JKOHTuyWFC+fPl4dCpmqWN4jSAIkSF93BakZs2a3JKDWdxmQ82RnYK1gbLf6dOn6cKFC+x5YwY6btetW5d69uzJRh2RI8yVFz0C5/Vxz95VnFJG2Md9KyaWWj+1T/q4BfNTvHhxnrFtRjw5c+/2Mqu3mEn+P/zWMaj7QbSlf//+7JHPnDmTiytxHy7eWbJkoaJFi9KKFSv8josV7AsqytXYnIgz37XFQSEQ8oaY4GV2uEo9dqGljaBED9QBUaI2bdpwXvynn37iqWSLFi2iHDlyUMOGDalKlSr05ZdfsocuCEJgxHBbEOQMMVMbYUerIUZQwOCXqlWr8kS3pk2b0rp167jNEfcjDw6FNhjxTz/9lOs5BHvickersjkRZ75rGxhugNCjIFgdGOznnnuOC9vOnDnD1egoZsOIUXzX0WZmheiSYI1Q+cSJE/l7heJKdEDge5dQCy6+n74b0jtGIobbgiAvCKzocQtCsO92+/bt6euvv6aTJ0/y/5Efx8jRBw8eyMmzEViKxbqjItpcIb7mwoULqUePHrwo3LVrF0d2UCyJGfT+GDt2LC8mPRu+kxkyZOBIkZGI4bYYqMhEnhDjF1GVK/z/EBDBXmDu9+jRo2nu3Lnc/ogCNkGIhFGjRvFiECkZiAiNGTOGdQkCiVqhgBKpG88GxcArV66wnoaRiFa5xUBLDeYhoxXMEzJ3EoHazSR3bk/gZc+YMYP/nyZNGqMPR1ARNQRUXP8931d2N1myZLx5g7qgHTt2UN++fePcX7t2bdq6daui15s2bRq34+bJk4eMRDxui5EkSRKWk4RqleT9BLsDMZe1a9fSypUr+QIr2Ac1JU9z5crF3rFnGzFiRLzXwxx5XDs9qUYPuI1Z88FAqPy7775jb91oxOO2oHIaLmaNGzemnDlzsmoV8n+QmnQCdvOsZTpb8Jxk/vz5qV69ejp9IoIVOXnyZBzBGF9v2xsUl3kD/QDf+/yByE+6dOm4ddFoxOO2GPhCfvvttxzaadCgAUtNQpAFuRunooUxl5y58UBVbdmyZdSsWTNFF1bBufO406ZNG2fzZ7ih0JcoUaJ43jUU/ny9cF9g3KEx0Lp1a9YjMBox3BYEF7EKFSpwWwPCN3369KF33nmH89+COtjNs7caGEqCufO3bt3ii6VgP/SeDpY0aVJu/4JugDe4jSl3CQHBoKNHj3JhmxmQULnFwZcRfa7379+nrl278tjPdu3aGX1YggLEq/cP1NSgZQ7xFbSFYdKYIKhBz549eSFYunRpdn6gGYBWsE6dOvHv+/Xrx5K8s2bNileUVq5cOdOkJMXjtokHPnDgQL7AffLJJ0YfjiUNkrSUmQMUEKFm46mnnqJ9+/ZJQZqNMUKApXnz5twCNmTIENbR37hxI3foeKrEEcH07enGwJWlS5eaxtsG4nFbHOReECpHmBytM5jVbTQSZg4NKVD7f9CzDaZPn845ScG+uCCg4o6sdsEVxvPffPNN3vzhaT30BlXqSNmYCTHcFgfFO/CyEeLp3r07iwQIoeM70cx7MIpW0QhZ4MQFrTooAELRpRhtQQiMGG6LMnXqVK4oh5LP888/T8OHDzf6kGyDVkab9y1Fb375559/WBEQ4XEnd0g4CZcKYzldDs32iuG2IMuXL6cOHTpQo0aNqGPHjobr5gphet5R0ZouEqwCitCQ00YLDzQKqlWrZvQhCTqgxnQvl0Ong4nhtiBfffUVPfnkk/xTUL+YTg+vWDzv/weaBChK++OPP6hAgQKan3vBHMRSFG+R7sOJiOG2GJA5XbNmjeEi93ZEjKk+7N27l2bOnEnHjx9nxatDhw6x7j4U0gRBCI4z4wwWZs+ePaz0U6dOHaMPxTaIwdaHq1evcgElwuKQMkWbDfpi0RXx888/izqaQ0PlkW5ORDxui4EcIERXgin9GIGV25qsetxWMtqQ5sXPkSNHsua+GaQjBeOIVSHUHUvORAy3xYDIimc8HZR/BMEKfPTRR5zHPnjwoCPH0QqCmjgzzmBh6tatS48//jiNHz+ezIZ4rYI/oDoFWd4BAwaI0RYeIqHy8BHDbTGio6Opc+fOtHjxYs4RCoKZ2bVrF2tDQ2oShlsQjBoyYiec+a4tDqpwM2fOzANFBMHMoBgNLV6QMJXRnIKgDpLjtmA72Jw5c6h///5S3COYGgxw2Lx5M61YsUIWmUI83F7ztCPZhxMRw23BUDkUppInT270oQhCggwbNoyFgurXry9nSoiHGqHuWIeGysVwWwCIVGCQCJSlLly4QNevX+epYIJgVsaNG0dr167lfm0sNgVBUA8x3AZy8+ZNmj17Nve23r9/nzdUjDdr1ixOaLFJkyb0999/U6VKlahUqVJ8G7lDQTAjEyZMoLfeeouFVURHXzDbWE87IIbbQNlHj0HGvNckSZJQ4sSJ6dSpU1zIA28FAhWXL1+m3377jaZNm0bt2rUz6nAFQRETJ05kcZVevXpxC5gUpAmBiFVhOlisQ+urxXAbRN++fTmEeODAASpYsODD+zdt2kQ1a9ak3Llzs0FHPhtAcEUMt2Bmbty4wUa7U6dOrI4mRltICPG4w0cMtwHExMTQhg0b+OLmbbRBlSpVaN26dbR+/Xq6ffs2b2XKlGHhFUEwM4gQoesB31cx2oKgHWK4dQYDQiBEAdnSF1980e9jqlatypsgWM1ww2ijbkOiQ0IwXBTNWyS4JFQuaM3hw4epdOnSFBsby4U7efPmlZMu2Ip+/fpR48aN6YcffqDq1asbfTiCiYl1R/EW6T6ciHjcOvLuu+9SxowZaefOnfxTEOwGokhI97z88sv8Pc+WLZvRhyQItsOZJXkG8fXXX1PXrl1tZ7QxzlMQAAouoaOPFsf58+fLSRGCFqdFujkR8bh1AuFx9GlnyJCB7IRnBrev8Q53UpiVZ3oL/5IlSxbKlSsXnTlzRk6JEBC3O5onhEW6DycihlsnUIwGPO1ddkMNY1srUfO4t8WImxrPYs3fZ581a1Yx3IKgEWK4dcLuhlsV3K64RiDKmatpO6RHYLjPnj2r2/EI1iOWoniLdB9ORK6MOnH37t2HLTN2IZLctr/n+npu62IXhr1/IbzPwN+C89atW2HluiHIIgiBcLnVyHOTIxHDraMuOUiZMqVeL2kJwyGFbeZNb2CQDSJEzz//fEj7gKH/9ttvqV69eqocpyAIcZFQuU4cOXKEf2KIiN3zm0rwfl6k+7I63osXrc9BKK8F9TO0d/mLEiX03G+++Ya97VatWkV4tIKdcalQnOaS4jRBa/EVeC958uSRE20QUuwWOsuXLw/5OXPmzKHy5cvbZpEqaIOLoniLdB9ORDxuHQ03dMkTJUpEVkZP7zASfL14vq1xsZuSvL0/9DyPerwWdPb/97//af46grUR5bTwkRy3job7iSeeIKuj1YVflXYyr5w59ue9T74dYbFbQjn5UHL1/h5rp3z/M888Q6tXr35YkCkIgrqIx62j4X799dfJypjdsERi/CN9b0pz9p7fBRKu8fdc9LdbqcIe6oANGjSg7Nmz01NPPUX58+fnUbUvvfSS0YcmmAjJcYePeNw6cO3aNRajsIPH/RAVw85qefG+Ai6hHkOg4/D8LqHHPDwGL0McyCgHuh0Qt4usRP369enQoUPUpUsXSp8+PW3evJmaNm1KK1euNPrQBLPluCNtByPJcQsaV5QXKlTIPufYhMYEXql3AVq41epKHq91Jbz3+zBzLUEgsEgdMmTIw7YyVKe/+uqrtGvXLsqdO7fRhycIlkY8bh2IiYnhn+nSpSOrYsmK7KhoTY225//em7/H2il/HQ5oK5sxYwalSJGCBg8ebPThCCbB/V9VeSSbOwyPe+LEiZQvXz5Knjw5Pf3007Rp06YEH49ajQEDBnBHEDqD0C3x5ZdfkpFIjlsHHjx48O/JTiynW0viDToJMS+sZFCKXgbYcoukICBkjgE7GLQjCECN6V6uEJ+/cOFC6tGjBxvvSpUq0eTJk6lu3bp08ODBgJGgZs2a0blz52jatGlcr3H+/PmH13SjEI9bBzzSj1ZVTfMuqDKzUVGSgw5EIG852GMC5v8V1gAEOl67eeg7duygffv2UYsWLYw+FMHBjBo1itq3b8+FwoULF6YxY8bwJLtJkyb5fTy6I3766SdatWoVF1jmzZuXypYtSxUrViQjEcOtAyhMS5IkifXncHsZJLMa74RIKJyd0HMCEsg4I//vpwYgnha7n6pyxa8dAUYsCjCL/tFHH6U6dero/tqCuavKI92UAt19LCBr164d537c3rp1a8DvbenSpWnkyJGUI0cO1uLo1asX3b59m4xEYrc6cPr0acqWLRsPXrAiVjTSwfLSvgVsWhfoeReaBWoF44WADkV/3segZ51HpkyZLC9AJJgzVH79+vU49yMX7TuJ8eLFixQbG8vz4r3B7UCT7P766y/uikA+fNmyZbyPN998ky5fvmxontualsSihlswhoS87FCMV0KPDeV3AR+rY6W+3osxo3OCgr3JlSsXPfLIIw+3ESNGJFgs6Q26Hnzv8+Byufh3c+fO5RA5Bucg3I5iSyO9bjHcOoXKIUYhJFCZHUEPdqREanQTatsKtu9wFwNmI1jlPC566O3euHGjrsclmJdIK8pdXlrnJ0+eZL0Mz9avX794r4dUDSI+vt41is18vXAPcLgQIsdiwANy4zD2//zzDxmFGG4dEI+bdJu9rVQvPBINcV8pVaW/CwUrGW0lixAUpVWrVo2aNGlCx48f1/XYBHMS+SzuqIeh8rRp08bZfMPkAFPu0P61bt26OPfjdqBiM1Se4/rtPVv+jz/+4LRnzpw5ySjEcOsA8ipWzW/rgZpGKhxj67kvFIU0JUproRyzWgbfrKA4c8mSJRx2nDBhgtGHI9jMcCulZ8+eNHXqVM5PIwL09ttv04kTJ6hTp078e3jqbdq0efj4l19+mYuK27Ztyy1jiBj17t2b2rVrx7oERiHFaTqAhv2jR4/q8VKWIJz8shpTycxuEHF8tZPYt10KF0BU8P74449GH4rgUJo3b06XLl1iVT+kMIsVK8atXp5xy7gPhtxD6tSp2SPv1q0bV5fjO4y+7qFDhxr4LsRw6wKa9tFWIHgpmoUYGg/WOmU1aid9mdbemxfv/rX3F5CdqVq1Ks2fP58Le+Cx3Llzh2d+HzhwgPOGIFWqVHxBLVGiBBcdBSocEqyNEQIsAFXh2PyBojN/8r2+4XWjEY9bBwoUKEDHjh1j1SiEDB2PV/W0JaVUVZjD7c9oWwV/n5nSKAoKfVCpe/XqVfZs4L3s3buX7/f8beB32DwywU8++ST17duXFa4E+2CU4bYDYrh1AFJ6yHNDNs/IggazEEmo21/1stGG327RAC0V8uBNg1KlSnE1Lxa1u3fvZu/ag6did8+ePbxBvaphw4b0559/yt+PIEhxmvagGnHQoEEPqxqFuIQ7BMRoY+2L2Y4nFPwthLR6P6jqxbhPFPdMmTKFtm/fHsdoA4TGESLHeFAMd/j222+5uPPzzz+XfnAbgcRI5ENGnIl43BqDfN7OnTtpzZo1lDlzZq1fzhGY1Uia9bjUDpFH8j7TpEkTclU52ntat25Nw4YN4+dCchIGXXLf1kZC5eEjPUoagnFwmEKD3JyvPq6ZMXvY1+zHp9VxPpe8FWmBUkMciZhMpGAIxLZt27g9Z+DAgdyaIwhORQy3RkC9BwYbvYIorLESZvYcnWq0gTs2lpwKFK+eeuqph7Pt/QlsCNbCiD5uuyCGWyP69OnDIXK0EVSuXFmrl3EkZs1zKyESeVe1W8WssgjyMG7cOJozZw7NmzePWywFayOGO3zEcGvkbeMCg1xclSpVtHgJwaYV3Xqi9sJH68/BM9Rhy5Ytho9VFAQjkeI0DZg1axbPfsXAdkE9AvV8m9HzViK9aias0E+P3DaUq7AgxlxlyE+KLoJ1keK08BGPWwP+/vtv/vn9999rsXvBpnh61I2clBaMhI5NlRnnCYAqcrSS4e8KbWTvv/++Jq8j6IPbHaXK5kTEcGvAhx9+SK+88gpvw4cPfyjlKESGP7UuW4bJfeZy6/E+FVeWB5Gq1ePzqFChAmsjYOYywuaCNVFzrKfTkFC5BiB8h+kz+fLl435TiNp/8sknju47jTQU6+/5Zg/thoq/92PLhYkKYIoTNM779+9PP/30k9GHIwi6IoZbI2Ck33vvPc7Jde3aldWhvMfFOQ27GVm9UCv8rKWMqRGfM9rDGjduTB988AFHtJy8KLYqkuMOHzHcGtO5c2eeAXv9+nWtX8r2WKGASiu837evEQ5k1O0+4xta51AjFKNtTdTIUbsdmuMWw60xZ8+e5QrzvHnzav1StkXCxXEJljJwyvl69tlnaezYsXT48GEevSgITkGK03SqMDeT4fa9sGt1oVe7qMqOXqPqFekq6IlbhZo1a1LKlClp8WL7v1c7IgIs4SMet06GO0+ePGQW9CzyinTfTvEejRai4TYvt8tSBj9FihRcN4LhIxhEAg8cQ0xgzDE+FL+XMLp5kVB5+Ijh1hAMRUCvKWZw44IihIZ3LlcMeHCsZHTVYsyYMXT16lUWZXnw4EGc36GrY/To0fTCCy+IARdshYTKNZwMhklGqVOnppUrV2r1MrbF21A7uShNTRJc/Pj0jgfbj1kWUhg2gtG5586d457utWvX0rJly1hyuFChQtSwYUOe63306FGjD1Xw43FHGi53S3GaoLbs6YkTJ2jVqlVUuHBhR57cSIytvypqIXyUnkMliyQzLqIyZMhAFStWjHMfFs4rVqygt956i4oXL067du2SIjYTAVmqSLWp3ORMxOPWiKVLl1L16tUda7TVxIyGwmoEK1jz/N5O5xr5bXjc+/bto6RJk9KiRYuMPiRBUAUx3Bpw8+ZN+vHHH+n555/XYveCIIQACteyZcvGETCXS3lKQNAWkTwNHzHcGoAhCMhxi+EWBHOQO3dumjZtGgu2LFmyxOjDEWTISESI4daABQsWUIECBahgwYJa7F4QgiJ1AXH59ttveQwoBpS0bduWjhw5It8ig5E+7vARw61B3/bChQt5/KAQHmJ0BC0G/1SpUoXmzZvHYfMGDRpw6FwQrIgYbpUZNWoUPfLII/T666+rvWvHoPVcZ0F7zPrZQU8BnR5IZaEK/eDBgzyk5Nq1axQTE2P04TkKVJSrsTkRMdwqcvHiRZo6dSp169aNlZsEwSjMUB3uLcNqJvLnz88932ghw9Q+eOPp0qXj25g4BiMu6KecFunmREQ5TUU+++wz/okxnk6DL85R0arKZprB+Aj2JHv27JzznjlzJou4wGhDxKVHjx7cNtahQwejD1EQAiKGWyVu3bpF48eP5xD5o48+Sk7DI0sqxtZ+RDLL28zfB3jZEGfxBoprGzZsEMOtA6JVHj4SKleJTZs20aVLl+jNN98kJ2K2cKig72dil88fE8dguDGKV9AWqSoPHzHcKvHzzz/zVCKntoCpobpll4u/3Qjns7WqChtkUi9fvvww7SUIZkQMtwqcP3+eq8lfffVVio6WUyqYAz0Mp90WW0WLFqU33niDp/qh2FTQDqkqDx+xMirQv39/SpQoEX3wwQdq7M622O0iL9iTIUOGcM0KhJQErQ13pFXlZClu3LhB169fj7OFgxhuFQRXvvzyS16hZ8yYMdLdOdYD9Bh1M7YPCf/i73PxfKZWDIsHIlOmTNzjDbU1QTuc0g527Ngxlr9GizA0PtKnT88biiPxMxzEcEcI9I8h6tCuXbtIdyUIpsZOxjkY+HtevXo1jRkzxuhDEVRm4sSJlC9fPkqePDk9/fTTXFgcCAyLwpQ53+3w4cOKX69Vq1Z05coVdvBQ+IhZFth++OEH/hkO0g4WAQ8ePOAPAwUtWgmu1ErUnNbFLiSntJN5/i8IahBui2KbNm3owIED9Pbbb1PWrFmpRYsW8oFoMY9bhX2EAuSo0asP412pUiWaPHky1a1blxX0MIgmEL///jtPmfOOyihl7969tGPHDipUqBCphXjcEQDpxNOnT2va86mV0ZZwtGAHvFMrCYXyw2HEiBH0yiuvsBEP1zMSzBUqHzVqFLVv3571NgoXLswRlVy5ctGkSZMSfB6mymEB59lQ06SUMmXK0MmTJ0lNxOOOgC+++IJDLaVKlSIrwV68Sb1asx6XlbGKME6w4/Q1zN4tZ97GW633ig4RpMLQNdKwYUNWWitZsqQq+xb05969e+z59u3bN879tWvXpq1btyb43Keeeoru3LlDRYoUoXfffZeqV6+u+HUhg92pUyc6deoUFStWjCV2vXnyySdDfCdiuMPmn3/+YY8bIRerXWjNGnq3gnER1MXXGHu++/7SJkZ8P5ImTcrzu3GhRkgVecknnnhC9+OwJSrGyq/7VGdDxhabN2jvi42NpSxZssS5H7fPnj3rd/eYJDdlyhR20DCYZvbs2fTss89y7rtq1aqKDvHChQv0559/8jhZD8iTY7gNfuKYQkU87jCZPn06Fze0bNmS1IYvXA7JbQv2XxAFql0IlK4JZ+Gq5XtE8SkqzGG8EfZEXhR1LUKEqFEV7v73+Qh3ezNo0CAaPHiw36fAWMbZxX8G1B/IS3vnpjHPHWHvTz75RLHhRqEjPPb58+fzIiHQa4WCGO4wwAoJITQUrHgXLKiJU4y2VcK4gvp4f+6eITVm/d7jgvvLL79Q586duUoYHtfYsWMpRYoURh+aQMTG1Pta7OttA8yQQG7a17tGKsTXC0+I8uXLs6a9Uo4fP05ff/01T6VTCylOCwMMEzlx4oQjp4CpDiaKCbZGifwpP8akRtvb80aoFDlL/MQFfPv27XT//n2+HqxcuZLmzp3LBauCvsppadOmjbP5M9xIeyDkvW7dujj34zb69pWya9cuDqErpUaNGrRnzx5SE/G4QwR/oChO6NKlC4c/hPCpnaQFX6ylDcze2CmqgjAnqpLLli1LzZo149C5J1/p3Sq0bds2evzxxw09VrNjxHSwnj17UuvWral06dIc9kb+Gtd0FI+Bfv36cRHZrFmz+DaqzvPmzctSuChug6e9dOlS3pTSoEEDbivct28fFS9ePF5x2gsvvEChIoY7BPDHCYMNxZthw4aR04n0grz2/v9LStrlwu70z9lORjohcAHeuXMnG+ijR4/yfG9UB8OrK1euHP3vf//jojbBXDRv3pynOELW9syZM1zljSLjPHny8O9xHwy5BxjrXr16sTFHWgQGHPUO9erVU/yankUBXtOXcIvTotzeS8UAoGIPUm3Xrl3TLKdrBfCH2LRpU1q2bBm3hzgZ8ZKFhL4b8fLXDlqczZw5k1577TWeGAgjbiX0uNZ7XiPvtIEUnTJ5RPty3bpDf7f/wHG2SRKMCsH6Bv1/CGs43WgLghAYiLbAk4PXrcAvciwyHSx+FEepUIsYboXcvn2be/HEaFt73rKgPb7fC6d9V1C5/OGHH9JPP/1EAwYMYJ1qIYE+7kg3Gw2sQqGjEsRwKyRlypRcVXru3LlIPhvhP0RyVbAzyIH26dOHZVPRaoQJUYKgFlKcFoJSWkxMTFiqSd5Gykmehz+JSu/7nZb7VAunFIBZGXjba9as4f9j3G+44xvtjBFV5XZBPG6FHDp0KGxdWU+o0GkXWyVylU47J2aOVHgGdiQ0E91OkRKt3gu8bMhiIkL33XffcUUyOlEEP0iYPCzEcCvE5XLxT7R7CMpx4oJFTRIarGGUUEpCht3M+DuXar+P/fv3U//+/bmFCOpqderU4WElgqAmEipXwM2bN2n58uX/nrDEcsoEbTDLYA0tDJqZDbSaqQdIW6ZOnZp7dkMZ/ehEJFQePmKFgoBmfcjkQd926NChPItVUAcneeKhjqw0M6b/3CCj63bFOVY9jhke9rhx49jL9ie5KWg3HcwOYHiNUs10MdxBgEoOROKR47bSOD8zTxezukcXCkrfp2/hntGY3jgngN7fe4h/oP3ro48+omrVqvHwEcHZjBs3TvFju3fvzj9DmTgnhjsIEKCHJrmVjDbzn8chWAsrG0wt0gWBOhPMALSnMd5x8eLF3H87fPhw6t27t4TIFYOK8EirwqPIjIwePVrR4yB56jHcoSCGOwibNm2ixo0bk9Uw0wXOX8jYrMenNk6KLmiBmb8nUFKEzjXy2Zi5nCNHDqMPyVrYOFR+TOO+fSl3DNK7jTB55cqVNf0QBHvjpIWKGuh5riJ5LUwJhMeEqNz06dNp8+bNInEqBAQDS37//Xd68OABRYoY7gQ4cuQI/4TusBA5D4VXEjWX0+lQgvWJG1VFHw4YC4n52xhugZB5lSpVeE63oBCHSJ7eunWLR8FCfRPTxTzTxxAiR21EOIjhTgDPykgqRAVB8EfLli1p5cqV3H2CSWAInQsKgeqZGpvJwYzvPXv2cNdB8uT/Pw2tZs2atHBheIWUYrgDgKk+v/zyC/8fK6VAOCF/qYbghnc7lFmr3QXnpQ7UEpNBz/YzzzzDNTEyEUwZTpkOtnz5cpowYQKnXJFa8VCkSBEeXBUOYrj9cP78eapevToNHDiQevToQZkyZXr4O39/5E4w3pFe7MxyoRb0w2l/FwiVnz59mv766y+jD0UwERcuXKDMmTP7FfbyNuShIIbbD59++int3r2bw16esn5vY+U7NMTORkmti68MFFHPu7MK/v4uzPS3ovZnUalSJb4QIyQqKMAhOe4yZcqwHogHj7H+4osvuE4iHKQdzA+o/MMfYd26dR/e591XaqaLjxHtTKG+fycZK8EaeC8k1fp+YgIYhotMmjSJ28PC9aYcgxo5arf5zzGGzkBN7+DBg1w3BYGeAwcO0LZt23iKXDiIx+2H2NhY0SQXBBvjHSlTcyEOAZYdO3bQDz/8oNo+BWtTsWJF2rJlC1eXP/7447R27VqWNoXhhpx2OIjHHcBwJ0mSxO8Jc5K37e89K404iJctOJFatWpRiRIlaOTIkVSjRg2jD8fURLn/3SLdhxUoXrw4zZw5U7X9icftB/TZyTCRyHDiAicYdq+HEP7NX/bp04fWrFlDCxYskArzhHBIjrt69eo0bdo01rRXCzHcPty5c4cOHz7Mq2ZBHUMlxkpwUlSmWbNm1KBBA+7xzpgxI1WtWpVV1tRQzBKsB7xtfP5wBps0acLtYVBRiwQx3H5K9F0uF0VHy6nxhxhhQS3s+l1KnDgxz+VGdfk777xD2bJl4wKlYcOGGX1o5sIhAizjxo2jU6dO0YoVKyhNmjT06quvshF/4403pDhNLbBChqLNrFmzyIqSkYIgmAMIsgwYMIDVsd577z0eRhJuFbEtcUioHMARrF27Ns2YMYPOnTvHs7d//fXXsOsgpDjND23atKHWrVuzmEL27NnDOrGCM7wrQVACQqXff/895ztRZdywYUN6/fXXKV26dAk+b+/evfT+++/TyZMnKW/evNymevToUW5ZxfxmbEmTJpUPwSKcPXuWax/mzJnDny16vMNB4sF+uH37NksYBvujMoKHOeMo+ejshERQ7A2uJxDh+PLLLzmqh7axjh07Bn1e165daf369TzoCNMKEXqHMBRypG3btqV8+fKxVrolcYjHff36dZ4eh46DXLlycZ8/aiD++OOPh7LaoSIetx+2b9/OU1wS0igXBCWoIdgjqnP2IHXq1PTaa6/xBnEnDJ7AKNBAY4OvXr3Kvb8IucPgg/v37z9sVT106BC9/fbb1KpVK44OpkqViiyFjedxe4OebYjzoGhx+PDhYXvZ3ojb5offfvuNSpcuTWbG7IM6rOpBqn3c3spc3vUJ/tTofH/vbbDVVPgSjKdw4cJ0/Phx1jfv1q3bw4pzeGcIo8Ijg771zp07qWDBgg+f560vgX189NFH/Bw8znI4pDhtxYoVHC0ZM2aMKkYbiOH2AXmoXbt2cSGBXCgdSFS0JsY7ITxGWqv2Oe/3I4WNwT8HPRg1ahTduHGDw6bYUKSEvDeMNeprrly5wjO+ccHHz0AgMpgiRQoudBLMCWwJOpWQ8kBRWkxMDN+PKAm+A+EgoXIvfv75Z2rcuDGHppo2bUpTWy5xnDa5YN5hHEq/h4GMjyxEg38Oep4j5L07derEMpgQbUEYHG1juPbkzJlTcesZZDMRJbQaTlFOO378OGuVQ9jr7t27nOtGWxjU9aAb8vnnn4e8T/G4vejSpQuHpaJ+SkfPJW7u2AterUTN494O8/1b8ry5XZot1Lw9ai0Xg/727/tZ6O15yzjcBM5NrVoc5YN2NXLWSo2291znI0eOkOVwSHHaW2+9xalXRFEQHfHQqFEj2rBhQ1j7FMPtBYrRsCr6kw7QNffleBdASxqiCPPn/kaZBsOq58mqxx0IM6nWmeU47HKc3sB7g3CUYE5QgIh2QN+2vTx58rAwSziI4fYCQgkvvvgixWa/Rb/R9zyXO5AXKtjr4hfpcZs5dxznffm0ESo5ZiXvLdz3L4JCkYP2MhgAVKAL5gP5bQyu8gX1C1h0hYMYbi8gtoLiAYgdAGjKPrzwaRhCNTVeF/pQL85Gn6+EqriNOmbDjbvbFdbTgp2XSM+b0d8Vs5PQ96Z58+bscaNX2EpEeeW5w97IGqkQVJR7D6JBUdqgQYOoXr16Ye1TDLcfjh07xj979epFTsc7bO65uMYxiImax4lGeP5vlguxHjnlUAoYzXJeQjkerd6bDKEJ7VwF4rHHHmO1R4izYNazkDATJ05k4ZrkyZNzYd+mTZtICZipjWLAkiVLUiiMHj2apW5Ri4BiNKjdQQUPURK084WDVJX7AWX7qPiEPKFUlSurcjbckzRBNbLZjHIwrHa8VkWP7wYihRcvXuRUHxTaLDELXI0+bHdUyOnQHj16sPGGfCzOG8RwDh48SLlz5w74PIzkxOLo2WefZa3xUCO5u3fvpvnz53O/PULn7du3Z+Ec72K1UIhyu91B6/LQ4P/II4/wwadNm5bsTufOnXl1BS1ZK16Q1cTf+w9kpH0NudnPm/cxJqRO5vs7K7y3QIgKm32BNwfDjWKo1atXs7hLqOhxrfe8Rp4Rwyg6efKI9uW6c4eO9xug+HjLlStHpUqV4t55byEb9NCjFS8QLVq0oAIFCrBDhxQqDLGRiMftB+QfIFEnBDdWvspgnsdY1bD5w1vVzPu2Fd+jFY/Z7Jjlu4DQL4xK/fr1OXcKkRfkv53gbCkB+u47duygvn37xhNI2bp1a8DnoXbgzz//ZEW7oUOHKnotjHVVygsvvEChIobbB4QxEDZBL6WTw78JeZq+Fynf3Lfv/WZ9L4EeFygF4H2/md6bYAxmvD4g9AqjgZnPGGKCHmIYBuRykWPFhjYkjJm0k1b59evX49ydLFky3rxBKgHV3dAO9wa3MbXLH+iPh6FHHhz5baXAg1cCCtX8VZwHwwSfnvkGjCAPAUUjM/XBmuFiYbZzoec5dsp7F5Rj1u8Eho0sWbKENSkwB/zvv//mWeDwxFHIhhYkCILgd1DyMoqIK8rd/6+chqlbCL97toTC3jCW3iBb7HsfgEFFIRnGqnrrxSt1AJVs4RhtIB63DwiZIOSEIgThX8zSShUpCR1foN/JcA8h3O+U0SBqCG8RG4wT+oYPHDjAEUXU73z44YccWl+2bBnLrlrZ4z558mSclICvtw0effRRzlH7etfnz5+P54UDaIrDkYOqHcarAhhbnEt436jgV1IEiP56tSdNiuH2Y7gxwcWsw+n1ulDEMdZR0RQVrV7HpL/JWGrtNyEDHO5rmvniLAhKgEcJrxQbdLMB5FVRUQ3vdOrUqZY+kWnTpg2ay8c1HSmDdevWsdyoB9xGUZ+/fe7bty/OfahGxyAqRDTQUqaEdOnScYSjWrVqPAcDY1wjHcEqhtsLrKRQTY7pPE4vgrGqsfKtFLfq+xAErSlRogQbMMN6vw2Yx92zZ0++vsOQVqhQgaZMmcIpBaRGQb9+/bi/etasWVwHUKxYsTjPx/Q2RGR9708I9HBj+/HHH2nChAlc/Y/Kdo8hx+IpVMRwe4EPEKPWKlasGPKJFJQTzpSrcERAtBiPKQsBwU7AcMCLxHUP/cZ2nw7WvHlzunTpEuf8z5w5wwZ41apVXLAHcB/sgJpggYANKQvktDHJDRPBIKn98ccfh5XnFsP9HyhI80gG4iQL5sfXkIpRFbT4jtn5++W51sGY+AsX25E333yTN3/MmDEjwecOHjyYt1A5fPgwe9wezxsjXBs0aMALp3AQw03E4RK0TiRJkoTncWfKlInMiB0vHJG+Xz0rwH0L1Zz2eQj2AwVYIEOGDI5QTjOCrFmzsqFGIRvC4/3796fixYtHtE/Ht4NhpQmjjapBCPUvXbpUtQ9MsB92FJgRnItnLgO0s3XHIfO4s2bNyqJeCMFjQ3U/bkeCYw03GvZRqICQSVJKTmPHjmWPW7DGqMqEjkvrUZFitJ2D3T9r9Hjjuqd3fttJ7N69m/XNBwwYQA8ePKCBAwdyVBfyq74qbkpxbKgc1YNffPEFrzRzU35zKAmZCLMa7GAXVN82NkHQ6rtmF48bwzXQ36w3RhSnGQVawqBgh1YwDDdZsWIFzZs3j/vE0U8fKo403KgcXLBgAU+J2fXJMWsMdfUqxtK6utnsRtupF1lBUJv9+/cr7ke2QzuYEUDgBgVp2CCAkzFjRh4Ag3GfmEAZDo50SdCzh148hMqtZrjMfGyC/eH56xqlIeS7rf/44jVr1nCLlKAdqKFCb3iHDh0ehs0h4IK6qqJFi4a1T8d53DWjXqLNtI6yUm4uGrASeoqvyEVU8PvdiF1oie9GOONLAyn6aaX0ZyTw/GBI4PFhNrQhqBAqJwt43JBUVQJC5nAqEVYPhuM87qnHPqG7dJtSUmq+HWiQiBUuTloi1dPB0boIzuj3lNDc9WDGy+jzEo5x9fe+/L2HcN+X0ecEQHykS5curJoGve1p06b5HbChCw6pKlfK8OHD6fLly4oe6zjDDa1e0PWzdkYfiumxg2ehJWacHufUBZvvYgMh/XDfgxbv3WiD7QFyn5gr/dFHH7HXbVh+G4jhjie5rRTHhco9c1sxKcYpoWK+iMUuJLtidCjTakYuWEhZT+Pt+9kFCnEHC33HGyIT4ffd9xyEWxCakMceyfkN93iQX23VqhW98847Yb+2YDyOM9wXLlx4KBZv14txPNwuxQ8VPe7QUOMibEf85Yd9OyKULo7NsIg20+drhvOhBk5qB1ObaKcabrPKmmrxBxvqUA+rXRiMDFl7XtsM+UszLjpCGaeakDdtlbSE7/dAaYRALzDaEqIrgrVxnOE+cuQI/7RaRbkWWM3QCBqgg0hNIKMbKDxvBQOdkMEOtlgx8v01a9aMvvvuO9bOFqyL40LlqKJECwSa4K2E1kbWX67RahdQI4/XaufKLMdshmNQA+/wv7+/HTWFkyLZB0ZIJkuWzBzyzg4RYFEKRFlSpEih6LGO8rj37dtHmzdvps6dO5PVsKonIghOIqG/U6P/flG1jIry2rVrk5ly3JFuVnAW/QHdckhve8Bc8GzZsinap6MM96RJkzhE3rBhQwkTB8Hoi4wgOB21o2zbtm2jPXv2cB+3oB+o4G/SpEmcHm3M5y5btiwtWrQorH06xnC7XC7WJ2/bti2HiZQILdgd3xyjVXt4BcGOqP13+Nlnn1H+/PmpVq1aZBocIL6ya9cubsPDDO5169bx51CqVCkqVqwYS6CGg2Ny3Hv37qUrV67Qc8895/f3YqzkPAjOAdoGaJN0yt89ZDcXL17MspqmmYTokBx3vnz5aOPGjfT2229TnTp1eBLbrFmzqEWLFmHv0ySfoPZgMguKMjADVYiPFduZBCFcINDiFKMNpk6dygbjtddeM/pQHMnKlStp/vz5VLFiRdYix0jp06dPh70/RxnuChUq8FQwQRAEp4AiqM8//5xatmxJGTJkILPglOK0jh07chtenz592PNG9BdOJELnkuMOkt/GCatWrZp+n5bVPOyoaMt6IHpHC0J9PYlkCEby7bff0smTJ81XlOYQrfItW7bQL7/8Qr169eKBLiiQRgX5kCFDqF278GZmJHZKGxjy28888ww5HTEiBp13m+vFC+YEg0QwKrJy5cr09NNPk5lwiuTpjh072MP2BQupmjVrhrXPxE7Kb5cvXz7e75yqNR3v/YagZ25G9Pz8Qn0ts3y3nPpddyowGCjGzZkzJy1dutTow3EsyfwYbQ+FChUKa5+OMdww2k7Pb/OFGyFx8fwcid2m3gmBWb9+PfcOFy5cmCVO06dPb77T5ZCqcrBkyRLOZ584cYLu3btH3uzcuZNCxfbFaVALCpbftqsH4nfggdsV8OJtxcpyqx2vngT6POWc2bvtq27dutyrXbp0ae4bNqXRdlCOe9y4cawfgomU6OmG8Aokt//66y/+rMLBEYYb+e1AUnJ2NdqhGmQrTWDyhxij+IQ6PlNQHz0Xw7jWwUAgRA7vDkY7TZo0ury2EJiJEyfSlClTaMKECTydDdXl+Gy6d+9O165do3CwfagcYgNp06YN+wTZBW9VNDteyK264NDzvFhxcIyg3GgPHTqUq5W/+eYbql+/vulPnVOK0xAeR/82wBCRmJgY/n/r1q05hQuDHiq2N9wADe9GGm41i4KUXHyDvZ7dLt56vh8rGz9fI+57nxA+/mZwe9+n5fcGRhuDkyZPnkyDBw+2hNF2Uo47a9asdOnSJcqTJw9vP//8M5UoUYKOHTvGn1042D5UDlKnTm0Kj1sNTzeUP347etZGG227YOW0iNnRu1YEeVMY7fHjx9OgQYN0e11BGTVq1OAoCGjfvj1Ln6IGoXnz5tSoUSMKB9t73H/++ScdPHiQT5bR6OVxR/J4K3s4WmPH8xgK0otuzu8GRnXCOXnjjTfIUjjE454yZQqLgAH01KMwbdOmTdSgQYOwR0wntrvUX48ePVjm7+WXXzb6cAQNcLox1RNpIzQf33//PY0ZM4bz2yh8shJOyXFHR0dzCxjavlD1j75uj/DK6tWr2YCHim0NN3IHqNpDDyPCFCgKsAPe+TOn5LC98Q1B2jWiIAjBgBFo1aoVh2L/97//yQkzKTDOKERDntsXSKDGxsaGvM9oO5+sSZMmsbh+uL1yaqGVYbFi37XauVkx2oITQegVk75w0Z89ezZP/rIcDunj7tq1Kw8ZOXPmDH9u3ls4RtvWhhteNuTkUAygtqEMx2BqYWCUFBgFO07P7622AOD3HRWt23E7cZGkJXIuIwPhcUQTMdc5kEaF2XHKdLDz589Tz549KUuWLKrt07aG+/Lly5QjRw4ORTidhC6S/lqErILkXK2J1b5nZuO3336jvn378rSpOnXqkGVxiMf90ksvsey2mtg2x3316lXu31YTj3iJGcKzSo8hVI/cLO/Pl0B9x2YeLmI2zPDZitGOvHanQ4cOVLJkSRo2bJhKn4rzlMw+/vhjDl0XLVqUoxdVqlTx+9jNmzdz/cDhw4fp1q1b3IeN+dqhdClBYKVp06ZcSY4Z3EmSJInze9RihYptDTfCE6VKlVK9DcboC5/TMUo4xIqfu6+htuJ7EOKye/du2rNnD8/YtloVuRnawRYuXMidRjDelSpV4v531EChZTh37tzxHp8qVSrOUT/55JP8fxhyGG78X2n73bx582jNmjVcIA3P2zsKjP+HY7htGSpH9R6+3P7GeArWRwyQdZHPLjJmzpzJudLatWuT1YlSaQuFUaNGcd3T66+/zpPT4G3nypWLC5n98dRTT1HLli3ZM8+bNy+98sorPCoV3rNS3n33XRoyZAiLgP3999+smObZMGgkHKLtWlGOir169eqpsj8751L9XUjNGs4Uta/Qz5dgH+7fv8/eG1rAEie2bbA0LK5fvx5nu3v3brzHoJcaA1h8Fz24vXXrVsUqdXjsM888o/jY8LpQSUM/t1rY0nDv27ePV1HZs2cP+BizGicjUJIHN/p8iRGyJkgzCeowd+5cunDhArVp08Yep1TF4rRcuXLRI4888nAbMWJEvJe7ePEit1/5Vnfj9tmzZxM81Jw5c7JwCkaldunShT12pbz66qscolcTWy7boJSGVZdZC3TMiJJz4tTBFHq/b6VdAFZ4//6iVXb8/mgtB4uitP79+1OLFi14QIUdUFM57eTJkzwF0gOMbMDn+HQa4dwG6z5CaPzGjRs8IAQV/fnz5+cQuhKwWBg5ciTnuZEr9y1OQ/g+VGxpuJE7wESWhLDjxUNtjPayHU9UNJH7X41jI1Hzb0UWzuGB3CiqoG3jbatM2rRp4xhufzz66KMsVOPrXaOQOViPdb58+fgnqsLPnTvHU9iUGm5EgJErB/v374/zu3DblW1puLds2RKwvN9MLV2Rejtae4K++3WqIbdDBbtZvvNmOAYt0LoOBrlV4DEAtkDnqvKkSZPS008/TevWrYszlQu3X3zxReUv6Xb7zaEH4ocffiC1saXhRv92oFC5HS4eRkzG0vu1BEGIOwEMSpDBIomWQ2cBlZ49e7JuOHLVFSpU4MldJ06c4KldoF+/fnTq1ClWpAOfffYZt4k98cQTfBvtYJ988gl169aNjMSWhhurKs/8U72MolY5Ln/HoHboUov9CuZCPlvrcuDAAVq2bBlNmzbN6EOxPM2bN+d2YbRnIfVQrFgxWrVqFQurANwHQ+4B3Ukw5ki/opL/8ccfpw8//JB7uY0kyg2/PwjwXlGphz60YHkEs1Rfot8Osqfp06cP+flOKsISwy1YmfXr13OFLy68ffr04RZQu8kcv/nmm7R8+XLOc2stuqLHtd7zGsXeGE6JkiaPaF+x9+7Q/in9LWOb1MK2HjfA/NNnn3025Od78uBWINLcpb88dkL707p6VrBPPtsfuMBC9AI5W/S3oigI4hZFihShggULcsWtZ3ISfApU5MJgocAH6lYQwcCkJU/octu2bexBwXDXr1+fChQoQLVq1eK/e9yXPHlyrjDG66RJk4asxu3bt7l3Gy1IlldKM4Fyml2wpeHGBSB16tTcbB+O4QZmvfDVTtKC1t5fwP/nxQUqjwXBAt9dAAOEsG/VqlXZoCIMvHjxYrp582aCz0uZMiUrXS1ZsoQGDRrE1b1QRoQnioItSElu2LCBli5dSmvXrmVJS2/gjeH3yG1aLbeNxU67du3IbqjZDuY0bGm4oVCDXARW6XYLh3uMtpajQhP8vXjblgIRkqjoqDjfG6OAl400FgqCMCjDA7xr9OEePXr0YU8t/obxExvEL9COg/vu3LnD/bCLFi3ihTmmY6FHFo+rWbMmbwD7g+AGHo8NeUoYeoTSEVrHT7OrjyHVh95ttB3heiYIts5xr1ixgho2bMirVcgDBgofWtV4C4IVgZFF2Bu9rHobTXj0s2fPpqlTp7LBR3U2FK2gW43wuhlBzh5DMDCZSq+Z23rmuIu3VyfHvW+a83Lctouz4kvRuXNnzne9/PLLD+8X4ywIxoEwNTxlSFEa4elimhNafrZv386eP2YkwygirQbdabT/YGyj2hw6dIh7hLFQgKevdKgEcvtwPF577TXdjLZRofJINydiO8MNRRusvpDj8lSXerS2/RWciUEXBG1BKBwzjRGqRiTMaDDLevz48XT69GkO3WMhAe8bBhKLfkwWVINffvmFKleuzB4zwt1oUcVCAV4+1LeCiXagNckTMRQE2xpu5LTHjRtHAwcOZNF5QV1kYITK5xOLSQsN4YBWM3LTyA9jQ5uSR9ErIVB8hvA0+l/N1KqF+ciIyiEa8Oeff7KoBowrDHuDBg3o119/DXvfKJBDYSwK6qBvPXr0aO4F/vjjjzmVB48fXnUgsKCAHnaZMmXItqg4ZMRp2CLHjRU9Cl7ef/99Ps69e/fGaZ2QXmVBiIwHDx5wyHfjxo08BhGXDXiU8FphXDB4oXHjxn6fi6EYOXLkYKELK7zPBQsW0LBhw9hThvFFiBvvHYY+EGht+/zzzx+O30QxHM4TiuhQEe87sKJatWpUsWJFDuEj/462L3j78MZxXitVqsRphd69e5Oe6JnjfvI1dXLce2dIjtty4I/k7bff5tU/Kkqx0vXtd5Q5zoIQPjDS+PvC39ZXX33FLVf4efz4cW7HwkW4SZMm7I375omhQoWFNELRVgBhc4g3oYAOoxhhUBHmRhgdalnoG/fn68TExNBbb73FW69evXiKF9refI02wByFL774gltWsaFiPlOmTHz+kGuHJCe0KHr06KHTuxashqU9bqz+IcaAfBFyVl27do3X8+yOjU1YUEQqywWd4LC422W5ugp4n++++y7NmDHDrwHGJWT69On89wfBE/RYI1UF0RMMcECu9o8//ghLxdAM4NhhULGhzQzvCyHw0z9epe3XNz4UdoEONkLiiC4MGDCADTAMs9KoIcLnMPbw7Hfv3s35cL3R0+Mu8ao6Hveemc7zuC1ruHEsngsEqkOhOWslBSnBmVhtoehprfzggw/YeCfE77//zl4qKqkxqAHjE2GskcaywwxpGFcUjUHsBQVsCP0j1A1xFOTH4TnjdxMmTOBIBIw2rk3wvpWAnDdqATB+8rHHHiMj0NVwt1HJcM8Sw20Zw41cEIpeEIaDDKIW+FahW+ViK5gb3wWlWT1xpKEgRYoiKRgpMxWWmQF435MmTeKFCcRSUNAGrxntZQDRiatXryoq4DMLYritgSWryleuXElffvklh6W0MtrA+0JqtouqYF18v0tQozPj9wvTqFBt/dFHH4nR9gOifcOHD2cDDuP9zz//UJs2bVgjvVy5clyEBi9c8E8UVPJU2JyILqFyrNxR2IEcEYYIYEUKrWJUYv72229cwIIw1AsvvMBhp4SUjLCyxVAC6BN/++23cS4oUj0uCOqASmd42ij4hOKYoAzIrKL2BkIvMOCowbFSbl9Pj7vkK8NUCZXvnjPAVNFgPdBFwggGFsUraJHImDEjT7vBhCBQqlQpbo1ANSryQehzRHUlBgmg0jJRokRx9oUVLvpJUZXpL3QXnTIluYIoIHmHKq2WcxQEPRg7duzDucWCcpCfRltcoNY44f+RISMmD5WjAhNGFqGld955hwswzp8/z9vq1au5IhxhOfSEYmAAClzQfgEjjz5GhKIAQlEo/MA+ULnpDxjtoIMyJAQuCH65cuUKh31RjIYWMAn1CoIDQ+WovkTfJy4EWI1euHBB0cUD4/tQnYpcNsJ2ECRA2xdCUci7eR+HkrYvQRtkPrc9wN8WCqu2bt3KaSu0MsGAo71LcA56hsqfelmdUPmueRIqVxV4yjC46EvET0gJKgE5Ic/YP6ihIcy+ZcsWypAhA7333nvxvlDRTzxOUReuqHvwQjz8ab177rfTosmJixFEuSDxCWMNtTAtiz4FAUio3KQ5bij/wMhCOjDc1RvEDbp3785bIFbvGxbBUQrhYidj7Q2Mtt0WI8GAp122bFlusxQEwaGGG32fkEWcP3++o6r9nIBdDZrTjLU3Bw8epAoVKhh9GIKTUGNIiJsciSaG+86dO6wkhJBb8+baTD+S1i/9cZJRc5q3jcl6SEMJgl5IqNxkhhuaxn///Xe8Pms1cdKFVbAPZvPq0dmB0ZYYaiEtTILg0HYw6O1+8sknPC3oiSeeICOLpgTBTN8ls31P0a0BY40IGdJa0dGWFFIUrIrM4zaPx71582Zu14LXLQhWQWshHrMZbcgGo+ULYysxGCN37txGH5Lg0HC5EDqqL7E9FwBIk+qJ2S6MgrXQcmY7DxHxeh2jwWAeTPyC9DC0EsqXL2/0IQmCYKTHDaUl9Fuj2AXa47oQFc3TlQTBjJipJxy1J5AehoohputhdoAgGAK0vyIdEuJ2psuuSXEahopAdEUP6uZ6i6KTJiXXvXu6vJ4gWAXo/2OkZJkyZShp0qS0aNEiLkRDwejMmTPFaAuGIlXlJjLc0BPHEBDIm2rFc6naUFTSpBT1aAZyX7tO0enSUtT9+5q9niBEihHDbAYOHEijRo1iHYXMmTPT0aNHqVWrVjxDGsJGgmAo0sdtDsON4SHo20aoHDk0rYw2RUVRVOLERDduktvlIrp5i9zicQsmxojc9pw5c9jDLlSoEE/6Qi67adOmlBh/O4IgWBZV/oIxpwQTvnr16kVPPvkkrVmzRjO1NNft2/xzzY2ZVDdbF/a8Y69eNUXRj0CGFiTKdyAuKVOm5GJREVYRzEiU698t0n04kYgNN6rH27VrRytWrGBt8g8//JCSJUtGWuK5QH935jP+KRXlziSQoTabyIlRYALT2bNnjT4MQfCPhMqNMdwofoHiEipVYbj1qCL3d0FO6CJtRG5RMBZ/n7XVvwehHv8PP/zAI3X79u2r8ZEJgmApw71z50422q+++qpurV9KvCkJoQq+WNVgh3P8SF298847nNPWalaAIESKVJUbIMAChbSaNWtS5cqVWeJUD5SGxK1+kRYEfyIuSr//UC5EGxj6tbWaFSAIqvVxR7o5kLAN9/vvv09FihSh9evXa9r65U10suSUKEP6kJSwEsqDxtm81K0EwYwiLkoXpBs3bmSDDWU0QRDiMnHiRBYKS548Oad6N23aRIGAhn+tWrUoU6ZMXHCN0bcovrak4YYqGgx27969NS9E8yYqeTKKUmi4Q0aU1wSbAMNdsmRJLk4TBLOHyiPdQmHhwoVcRD1gwACOSlWpUoXq1q1LJ06cCPi3BMO9atUqbneuXr06NWjQgJ9rJFFuJMSCcP36db4IXLt2jVcdr7/+Oq1evZqOHTumq/pS3dw96Hq5XJR2818PK8pDRUm40duzUfp4qWQWzMCDBw/oscce46lfY8aMMfpwBIvhe63X8jXK1f+AEidJHtG+Hty/Q7+sHKj4eMuVK0elSpViESIPhQsXZt2RESNGKHrNokWLcu2IkW2W0eHM74WwQ9euXXWXTPzuxBiKTUoUmzfrwxB3qAQLN/r+PmDIHfroDkfa8MzHuHHj6NSpU9SmTRujD0UQdOP69etxtrt378Z7zL1799hrrl27dpz7cXvr1q2KXsflclFMTAyLjFmqqnzy5Mk8t/eNN94gI7iTIZrupU5Nj25PRO7Y2LA83WB570C/V1rNHshj99wfqOo9UMuPPwPp8fKDHZeRkQDUDfgO2JDIhHacPHmSvYAuXbqwVyEITqkqz5UrV5z7Bw0aRIMHD45z38WLFyk2NpayZMkS537cVqp38Omnn/Ic+2bNmpGlDDfeYI4cOQxZcdR8ZjhFFUtB7sRE0SlTkuvWrbCNtzdKn+95XELtZsF+7+92sPvDfZyR/dKCvpw7d45bMhEu/OCDD+T0C46aDnby5Mk4ofKEaq98Oy2QLVbSfTF//nxeDECzBNr/luvjRjWe3jw2bzg9mjc5pTzvotNVibJmz0KJYm6SO+YGxcbEqGKAlIbRg4nAhGtYAy0g/C0GlDwvkmMJl1BTEULkQEsBBTTwBNauXStFaYLjPO60adMGzXGj+ylRokTxvGukf329cH9Fbe3bt+dRuGiDNpqQDffjjz9OM2bMULxKUYtkKe6TO1Fyupcar+kmd7IkFHX5PnvdWhgDtUK6/sLF4RCq6Ewoz9UbCZerx/79++m5557jMbpbtmzhNhdBEOKD0bZo/1q3bh01atTo4f24/eKLL1JCnjZkvfHz+eefJzMQsuEuWLAgS52ifD5PnjykF2lS3KWbmdJSbAqivEVPUUyB7JQ6URRFXb+hiXFQy9iFYrQVhev/WwhY2fhZ9bj1BDPtMcUr2OK4ZcuWlDFjRr74BPMaBMHpWuU9e/ak1q1bU+nSpbkne8qUKWzLOnXqxL+HaBGKO2fNmsW3YaxR6Dl27FhWIvR461goG9luGbLh9hS9/Prrr7oZ7urfv0NJEiWmB6n/vX3qcjrKHvOA1ux4n5yGZyEgxs++IJpVtmxZTklBLyFVqlQBH3fkyBEaOXKkGG3Bchghedq8eXMecTtkyBA6c+YMFStWjHu0PbYM93n3dKMYGy2WKPjE5gEy34g8W6qPu0yZMlxZ/ssvv2h+gPnGf0pp812lGzeTU9TxlORK5qbsRc/S5pojNXk9K3uygj3AgJAaNWrw/2fPnk2vvPJKQCEkjNH95ptvqH79+jofpWBH9OzjrvjcEFX6uLeueU/T4zUjYTUjd+7cmX777TcuhtEad9oHFBOT4t//577NoZFLMf49EOH/8ZZzFazDlStXqG3bthyW81zkAjFq1CjKmTMn57gFwXK43OpsDiQsw41VPhz1Pn360NGjR0lLCuQ5S+Xy/U3u2ChKmfIuJb4VRW63DE4IhkQNgmO2RQ1CcjDaMNYLFiygEiVK0M8//+z3sadPn6a5c+fSW2+9pbsQkiComuOOdHMgYRlu5Ln79+/PifsCBQqw8gzCdVqEyQukvUA/H81HrthoSpXsHiW9QnTv5H/JbpMbPLMZBsG8KRGEveFl4+8IuTPk3OrUqcMDDaDW5AsEJpAD79ChgyHHKwiCxfq4kd8eNmwYvfvuu9zXBt1XiD989NFH7IWrRfT9KPrheAFKciI53c38gBJFuyjRPaJEt+N73FevXqXly5fThg0bqG/fvqwnCxDOP378OF24cIH79fDT+/+oDsSUM2woVMibN6+MQnQAZjLakFusVq0aL4Lxf+gpAxhu/E3t3r07jhIa8t5Tp07lilgZJCJYFVzFIy5OI2cSVnGaPyC1CMUmGHCU2qPkHkZQKfAq0PriaX/JN3cEuW4npuRnEtODVMhlEGUreZZuLM5MV6PO0uRGtVgtCgb48OHDNG/ePLpz5w6lT5+e22heeuklvr106VI+bg/4HUa0QfkGP2/cuEEHDx58mEuEwYcX071795AMuD8PTokkqROIF3mIilalt90uoBDt8uXLHBb3FjeCtjJavbAYHjhwIN+H7zS0FKpWrcoRL0GwanFapWcHU+LEERanPbhDWzYMdlxxWlged6D53JBBXblyJX3yyScsDYfxZx07dqTs2bOzwfS3ofx++vTp7EXA6OP/8JKjfr1J0Vdv0K3bNyj29k2KvXGDDv90nmLWHaJ7d2OozuRP/3sDSejxQo+xl42pZVDGwbjRn376idtoUJHbokULNtTY8IXxJ3mH3j0I0ONiiLFvODbv8n87hWX1Jp7ym8EjVL0/DyM/G0R8hg4dylXky5Yti6dICMEI/E1gUbx582YWjcB3FL2kaGcRBMGZqOZxewPDiwli48ePpwMHDgR9PIQj4EEg7B6IRKmSUYoMySll8iKUNXNJSpwuA7kyp6VkdxNR9H0iVxKidAeu/atde/Isrb78BYULvG2EIadNm0b58+enrFmz8qbn7HE74q31boeFTCQRlUWLFrGEItJOWHRi8xfhuX37Nj8W/aTbtm2j1KlTc30JhCIEwcoed+Ua6njcm793nsetieH2gF0fOnSIDTkqZv1tyDE/88wz7OF+/vnnNPjrtUTpUlHSJKkpWXRqcmVOTa5sKah2lYN0MzYp7Z9bhFJcdlPaP26QO0k0uRNH0/V8ycmViGj7tJ78upEaBoQjcUwQmfEG9yHn6MlBCqFjlfSBv6luar0XhMDRmYFWLmggIxyuBEw2AogqCYLlDXd1lQz3D84z3KqFyv0BDwJFX0p4fmN3OpMnLaV+tQzF3klEUTcSEy5TSa5H03PPbKeR2bZSimzHiBTYzEiNAkKW8G4wBg5hSYTzMcgBKQB4Rght2p3aSVrQ2vsLNNFcN3u1vdLjC+d7hsUsUjrHjh2jmTNnKjbaQAy2YCeiMO8iwulgUZFOF7MomhruUPi26jgyEwhhevLi8I4AetaRi3QCvkY7Klq9+k2ze9ta5r+x6EM9B1JJErkRBEG3Pm6ngmHt//zzD3tNwOyeo2q4XXEMeSRY6ZxpscBA9TgQiVLB8bhU2hyIaTxuK4Cc5N27d6la9IuUNOrfQjW7FFolhJrvz+7nSkl+21MxLghORkLl4SMed4geN7hDtyzrRQraEuy7IIZbEIRIEcMdAuhHB/foTrzfifEWlEQUYLhRPyGFZoLjEa3ysJFQeQh4+rhdfpTtnR4CFpQBwy1hckHgFot/t0hwO7OqXDzusHA7zmirEVFAdTo2O+I9QjWhc3X//n3WLBAEQQgXuYKEQJo0abg3vfsX7WlBh9W2MtjBxETUeK+ePnCriLCEgvd7sdP7EgStwICRiIeMuMmRiOEOMVSeO3du7ueWi3P4OPncIb/tb0ynIDgOCZWHjYTKwyCUqWFWIpi0pxD5ecF3R4HKsCAIQkDE4w4BXHBPnz79sLrc7kY7kGFSK9Rt5R74cI8dhls8bkEginL9u0VClEODV2K4Q8Az5MHJqOl9W9VoR3Ls4nELwn9IqDxsJFQeAqgGLlSoEO3fvz/kE221cLO/47XaezBzjlvC5YLjkT7usBGPO0SKFy9O+/bts7V3qUrrl4KxmE7k8ccf53G2+A55htcIgiCEgnjcYRhueNxm9pgiNbwwtHGMbVR0yAbYsw8x2nGpVq0atxX269ePB9YIgtO1yiPdnIgY7jAMNwbBnzx5UvFz9AwxKxEBCdn4+um/FoMcfkvhtGnTaPv27VSwYEGeyS0Ijs5xR7o5EDHcYRhuEE643C45aMl1R0bTpk3pyJEj1LJlS3rttdfYkAuCIChFctwhAgGWtGnTcrj8+eefJ6MIx+NVq/1KD2/bd3FgNw8f36GpU6eyB/76669T3rx56dlnnzX6sARBP+AsR9rO5SZHIh53GO08yFHeuhV3tKcVZDd97zOrbrj3AsPOYXl8lz777DOqVKkSvfPOO9LfLTgKyXGHjxjuMPu5QxnLqKfhCcXQeXLXWuw7IsIohrOy8f74449pz549NHfuXKMPRxAECyCGOwyMnvBk9xxzqAsKq1OhQgVq0qQJDRgwgG7fvm304QiCjn3ckRankSMRwx0iaOG5dOkSPfHEE2QUZvZEw1lUeI/EdCrDhw+nM2fO0Pjx440+FEHQB6kqDxsx3CGybds2/lm5cmXSA6sZtHAWFXbOYysFrWEdO3akYcOG0e+//2704QiCYGLEcIfIhQsXKEmSJJQ5c2bSE6sZcCF0hg4dSjlz5qTnnnuOvW9BsDUulTYHIoZbAfCAunXrRg0bNuQq4FSpUulaXa2mqIpecPjbpFXrZiVdunS0evVqlkStW7cuXbt2zehDEgTNkKry8BHDHYSLFy9S2bJladmyZXTv3j3Wmu7RowfpgW9LlJXCyd6Ka4JycuXKxcb777//psaNG9ONGzfk9An2xKAc98SJEylfvnyUPHlyevrpp2nTpk0BH4vI18svv8zDpTAgSK9rfzBEgCUIULW6e/cuHT16lDJlyqTPpyI4mmLFitHXX39N9erVo9KlS9OiRYtkIIkgqMDChQvZ+MJ4Qz9h8uTJHN06ePAgi2v5gms/rvvo+Bg9erRpPgPxuBMAoUpU+UKaUoy2oCdVq1alHTt2sLJaiRIl2DuoUqWKDCYR7IMBHveoUaOoffv2rFZYuHBhGjNmDEe5Jk2a5PfxUDQcO3YstWnThh555BEyC+JxJyCy0qFDB4qJiaH33ntP30/Fk8+OiqYoj9CL+98qjLX3F5AZiJNvj4qWsLgGIDz3888/0+LFi+nKlSs0YsQIGjhwIE2fPl2LlxMEfVFjSIj73+dj8JM3WPBi8wapTiyG+/btG+f+2rVr09atW8lKJHaqUUYOERfDkiVLcjHQxo0bOZ948+ZNDpvMnj2bTp06xRdN5EOMwDtHXDtJi4eKYmbASvl2K5MiRQpe7QNMpFuyZInRhyQIpiNXrlxxbg8aNIgGDx4cr14J1/4sWbLEuR+3z549S1bCcYYbxrhmzZp0+PBhvp0+fXoOQ3q336C6t0WLFhxSQY7RDIYxEk8bRt8snroQPvAqJGUj2AYEEaNU2Af9u6jF4B4Pvt62r8ywN27M9fa5z+w4ynBjtVWrVi32qr/99tuH7TeXL1/mnAd6aDFAJGnSpJb7IBPC7XKoLqDNgFeQNWtWow9DEFRtB4t0HwBG29tw++PRRx/lGRO+3vX58+fjeeFmx1GG+6+//qJDhw7RqlWruJIQVKxYkWzPf/lxwdrggoNCNUEQQgcOGdq/1q1bR40aNXp4P26/+OKLZCUcZbhPnDjBPwsUKEBOQvLR9jHcUFUTBFugYnGaUnr27EmtW7fmFCiG+0yZMoXtQqdOnfj3/fr143TqrFmzHj5n9+7d/BM1UFDOxG0sAooUKUJG4SjDnTFjRv6J0LgdQG5mwoQJNGfOHBaGwZfx+eef52pkwV7gs4bhzpYtm9GHIgjqgBReVISG2xXa85s3b85DooYMGcJ1TdBMQAQ2T548/Hvc53HwPDz11FMP/4+q9Hnz5vHjUeBsFI4y3J6LHlZUdgDV7927d2ehDgw/mT9/PvXu3ZumTp1Kbdu2DWlfqLg8cOAA7xMb/o8vZsqUKTlCAfWgcuXKcfWmnfL/VgEdEBgnKzluQYiMN998kzd/zJgxw++i2Ww4ynCjIhdG6N133+XCtGeeeYZl7KzK999/zyEbtKyhivLWrVv02muvcSsEwkFKZoajUA867J7eYDwHk6oQBipTpgzduXOHFwXNmjXj30OEoHjx4qzkhZ+eLVhhiBAZWEgBj2cgCJbHgFC5XXCU4YaRXrFiBQ8LqVGjBsvede7cmaxsuJGngVcMUBEPsRj0pqMPPZDXjep6LF727dvHBgFVlePGjaNnn32W8ufPz4sBXxA+2rt3Lz8HP3/66SeWC8S+PAbFY9Ch8AVRAysviswGuiCw8PQO2wmCtVHBcJMzDbfjrqyQudu8eTP/P0OGDGRVdu3aRWvXruWedG9QdYyFSbt27Vi8Axd75HO8Wb58OX344YccAsICZvv27ex1w8v2Z7QBdHzr16/PxRsIye/fv5+9dRRqYJEAjxxCNgg1oWIfuSN48dD6FdQx3EiJyGJIsA0GDRmxA47yuD3s3LmTf6I1wEzAkLpcLu419AcMIzxkhK6hClS0aFG/02qgroUL/fHjx9loe8KsAB7y8OHD2WDjMZHg0dH2bVHasmULjRw5khcP8OwRyndE251GoNYACyWkQARBEBxnuFGZ+8EHH3CFOSqxzcKff/7JYxwRkn7hhRfYU0YIHB4r2hDQf/7rr7+yl4s8NMLRmHSTOnXqePuC4cc+POF0TJeCshC8YBgAeOs//PCDZu8FU3eQksAxw/tGa4UY7vD4559/6KWXXuLvAsSDBME2cEW4vlXldsFxhhvTYDC4YcOGDaapjl6zZg1LrELZBwpuMKoYbgIjDTnWVKlS0WOPPcb5a+S0ESnw5LWDMXfuXFq6dCl99dVXHNZG+BzvHYV5eqQlkJM1sm3C6mAy0bFjx+jHH3801XQiQVBFGCpScSi3M8WlHGe40UCPamk9DFew9h6Eu3FBRsFXnTp1uD8Q1e5qAgOPCnNsRoCxeJCVFcIDqk4NGjTg76wgCILjitNQCQ0PNEeOHIYdA0LV8Jzhjc6cOZPKly9P06ZNo2+++UZ1o20GPEIFZuyFNDuo9t+zZ0+8AkRBsAVSnBY2jvK40caEnDFCxxCyQC4Y+WJsCENC9u6NN97Q5LVRWIZZyqjmhoFGPhu59uzZs5OdgeFGLzgiHZkzZzb6cCwFUhpADLdgSyTHHTaOMtzo2UZbE4qmUF0NY+rZEK6G4UZYUi1ZSShdoTgM+WUUa8F4ffTRR1wJHqjtym54hFmQrxdCY/369VxQKGppgiA41nDDs8YUGH+TYKBfjtFuy5YtCyiHF+rsZFQBoxIchWWvvPIKtWrVynECGp6+Y49Qi6AMqOChaLFp06ZyygR7IsppYeMow50QEGOBchh6oCM13LjownP//fffOa+O1i2zVLDrjacnHf3pgnIgXXv16lVq06aNnDbBnnA3WKSSp+RIHFWcFgz0y8LQIqQdLvfu3eP9QI0MU2eqVq3qWKMNxOMOj6NHjzoyQiMIQnDEcHuBEDqMLMKTmNsKpaq//vqLlIJcOS62KCpCTltER8TjDhcUMF67di3s5wuC6ZGq8rARw+0F1MoQKkdYF73HkO0cMWKEIi/766+/5lnYyJFDqUwqgeOGyiXHHRqlSpWi7777LqLojyCYGqTP1NgciBhuHzCgY+PGjTyTun379uw9ozrcl9u3b9OmTZs4H44qdHjrkFOF4fdX/Ob0ULnkuEMD/f0ocETIXBBsiXjcYSPFaQnQsmVLmjJlChti9HejZcwz1vLIkSNsjHLlykUdOnTgEDnGWgpxEY87dGCwe/fuzVEbGHBBEARvxHAHGZaBCVoQS0HYEpXnmDf93HPPUa9evXgqFnTDZdRiYMTjDp3u3bvTpUuXWAPAyYWNgs2RdrCwEcMdBPRiY8gD8tgIicuFVD+PG8/BkA14oGXLluU+e7ONYlUbjECFFC62fPnyGX04gqAdopwWNmK4FYCpXYL+Hnf//v25QNCbiRMnUseOHW0Z5cAIT7w3dDUYNRRGEATzY7+rn2AqPONHISSCnmTUDWDmuBKOHz9ONWrU4Elqp06doi5dunAxYPr06VnjHTroBQoUoKJFi1LJkiXZK69cuTLXJCh9DTOARc2MGTOodOnSPML1888/l8iOYHvcbpcqmxMRj1vQfKznsGHD6LfffuMhI5iCBiOOmePB0g5QDkO4HD3N2MaPH0/NmjWjbdu28bAYVPsjheH7E6Mw0R2wefNmxXPLtQST0bCQgPztzp07ue2wUaNGHArH/c2bN+cZ8VjUIMKAWgpBcESOm8PlEe7DgYjhFjQHIW8PdevWZaN1+PBhHm0aCExwg1Y3PFEPMPRQosOWEBiFifnV8FwhpGMka9eupW7dutEff/zBtxEpQAQBuXtoBaDQEXUAUOwL9r4EQRCAhMoFXalWrRr/3LFjR8DHnD59mtvvmjRpEpZWd8GCBdlTT5MmDRkBwvpoD8yZMycbZhhrKOmdOXOGf4dIADztIkWKUOrUqfm2GG3BcUgfd9iIxy3oCkLemJaGkDmMm2+4HOHxFi1a8NjTcHO9mHsOL1YvY4iwPQzx+fPnadeuXfTuu+9yiB6DQipUqED169eP8z6givb++++zFsCoUaM4dC4IjgMFq1ER5qjdkuMWBF2AYW3bti0rzaFXHkYaHjKqqhEyRigZhj2can6o3n366adsEAsVKkR6tQxCRc/Dq6++SqNHj+YiukC89957uhybIAj2QzxuQXfgiaJ6etq0aTR37lw22vBIc+TIwYNZYNQQ7g4H5I2zZs1KPXr0IL1Af7+HX375havbBUFQUlgmxWnhIIZbMAT0KmNTG4Thz507xyNVUQinR7/3ggULqE+fPvTyyy9zWHzlypVivAUhCG6Xi9wRhsrdDg2VS3GaYCuwGEA/NAwoZGn1ANECKLpt2bKF8ufPT9WrV6e+ffty+5cMVxGEAEhxWtiI4RZsxSOPPMLh6iFDhnDfN3Ld6O3WA+Tk169fT+3ateM0QLly5Sh37tzcDoa+dTHigiCogRhuwXbAA0boGi1lmLKFwTDIfesBqsmxYEDrF4w1WtqWL1/OCnBQfRME4T8gvqLG5kDEcAu2JFmyZPTZZ5+xUhmK1ZDvrlKlCmuBwxufPn26prKoaA1DpXny5MlZxQ38/vvvmr2eIFgzVO6KcHOTE5HiNMHWYPQqPN8lS5bQ7Nmzafv27fT1119zARukSKFxDm1z5KbhLWNDXzVmq6NNLVyg+DZ48GAulsM+X3rpJRoxYoSq700QBGcihltwROjct4odo0IRPkcYe8yYMXTt2rV4HjuGoqAXG5Xp2CDq4vk/trRp07K+OIrRfIVioJBWrFgxll8VBCE+bpeb3FGRecxu8bgFwTnA6GJgCTb88UP97NatW7xBCAYFbpBlvXnzJveZo7AMQ0zw03N79+7dNHXqVHriiSfYeCdJkoQSJ07MRh1a61BNEwQhANzKpb9y2sSJE+njjz/mOhRMFsTCHWm0QEAUCjMPDhw4wPLFqJ/p1KkTGYl43ILjgbeMXDQ2TOaCxnj58uWDnhcYfCi1QZoV071g2B88eMA/cSGAgpsgCOZh4cKFLM4E4w3VxsmTJ3P9y8GDB7kDxJdjx45RvXr1qEOHDjRnzhxu+USRKdJpKDw1iii3glgDwopos0E4EZ6KIAiCYD/0uNZ7XqNaVCNKHJUkon09cN+nH93LFB8vWjQxK2DSpEkP78OUQhSQ+qtB+d///sc1MYcOHXp4H7xtpMAwXtjUHrfHtuOEC4IgCPbEc43XI3f8wH034iEhD+i+4sdCzwHpL4gjeVO7dm3aunWr3+fAOOP33mDiH3QaEFlDesy0hjsmJoZ/5sqVS+vjEQRBEAwG13x4xVqAbg20aG4+u0qV/WXNmpXu3LkTr7gUmzcXL17k+pQsWbLEuR+3MfDIH7jf3+OREsP+vOcUmM5we+Yanzx5UkLlgiAINva44aBpOcsetSTIHaulaDhy5Mh4xnXQoEHcjukP3w4QRBcSGh/s7/H+7jed4UbuADkEfJhGHqwgCIKgHbjGe671WuIpBlVrTHCfPn3i3OfrbXskidHS6etdQyzJ1/B7e/P+Ho/ukYwZM5KpldNgrGG8xWgLgiDYFyte65MlS8bH7L35M9wI0WMY0Lp16+Lcj9sYJ+wPtHT6Pn7t2rU8yMio/DYQyVNBEATBEfTs2ZO1F7788kuuFH/77bdZ+tjTl92vXz9q06bNw8fj/uPHj/Pz8Hg8D4Vpek0eDIT0cQuCIAiOoHnz5nTp0iWeVwABFqgbrlq1ivLkycO/x33eMwzy5cvHv4eBx+wDCLCMGzfO0B5uxX3cgiAIgiCYAwmVC4IgCIKFEMMtCIIgCBZCDLcgCIIgWAgx3IIgCIJgIcRwC4IgCIKFEMMtCIIgCBZCDLcgCIIgWAgx3IIgCIJgIcRwC4IgCIKFEMMtCIIgCBZCDLcgCIIgWAgx3IIgCIJA1uH/AIjX2yDmx0HgAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import TRITON_SWMM_toolkit.examples as examples\n", "norfolk_multisim = examples.GetTS_TestCases.retreive_norfolk_multi_sim_test_case()\n", "\n", "event_iloc = 0\n", "analysis = norfolk_multisim.system.analysis\n", "\n", "self = analysis.plot.max_wlevel(event_iloc)" ] } ], "metadata": { "kernelspec": { "display_name": "triton_swmm_toolkit", "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.11.14" } }, "nbformat": 4, "nbformat_minor": 5 }