{ "cells": [ { "cell_type": "markdown", "id": "3d12d1e1", "metadata": {}, "source": [ "# Hydraulic utilities\n", "\n", "This tutorial covers two utility modules:\n", "\n", "1. **`compare_series`** — model evaluation metrics (NSE, KGE, RMSE, …)\n", "2. **`friction_law`** — pipe/channel friction coefficients (Colebrook, Barr-Bathurst)" ] }, { "cell_type": "markdown", "id": "ee4b83af", "metadata": {}, "source": [ "## Part 1 — compare_series: model evaluation metrics\n", "\n", "Compare observed vs. simulated time series using standard hydrological metrics." ] }, { "cell_type": "code", "execution_count": 1, "id": "937215d8", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from wolfhece.compare_series import (\n", " Nash_Sutcliffe_efficiency,\n", " Kling_Gupta_efficiency,\n", " Root_Mean_Square_Error,\n", " Mean_Absolute_Error,\n", " Pearson_Correlation_Coefficient,\n", ")" ] }, { "cell_type": "code", "execution_count": 2, "id": "ab8f7edb", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "NSE = 0.9671 (1 = perfect)\n", "KGE = 0.9609 (1 = perfect)\n", "RMSE = 0.3644\n", "MAE = 0.2981\n", "r = 0.9843\n" ] } ], "source": [ "# Synthetic example: observed vs. simulated discharge\n", "np.random.seed(42)\n", "t = np.linspace(0, 10, 200)\n", "observed = 5.0 + 3.0 * np.sin(t) + np.random.normal(0, 0.3, len(t))\n", "simulated = 5.1 + 2.9 * np.sin(t + 0.1)\n", "\n", "nse = Nash_Sutcliffe_efficiency(observed, simulated)\n", "kge = Kling_Gupta_efficiency(observed, simulated)\n", "rmse = Root_Mean_Square_Error(observed, simulated)\n", "mae = Mean_Absolute_Error(observed, simulated)\n", "r = Pearson_Correlation_Coefficient(observed, simulated)\n", "\n", "print(f\"NSE = {nse:.4f} (1 = perfect)\")\n", "print(f\"KGE = {kge:.4f} (1 = perfect)\")\n", "print(f\"RMSE = {rmse:.4f}\")\n", "print(f\"MAE = {mae:.4f}\")\n", "print(f\"r = {r:.4f}\")" ] }, { "cell_type": "code", "execution_count": 3, "id": "498150b6", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxYAAAEiCAYAAABkykQ1AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAXfBJREFUeJzt3Qd4VEUXBuAvCZCE3nvvvUlHAelFigIiggIKCoKA4g+iKCpIVTrSBETpXQRBEBEU6UgTpAuh94QWSjL/c2bYsAlJ2GQ32fa9z7OkbMndwr33zJw5x0cppUBERERERGQHX3vuTERERERExMCCiIiIiIgcgjMWRERERERkNwYWRERERERkNwYWRERERERkNwYWRERERERkNwYWRERERERkNwYWRERERERkNwYWRERERERkNwYWRERERERkNwYWROQVvvvuO/j4+CAgIABnz5594vpatWqhZMmSkX53//59jB07FuXKlUPq1KmRNm1alChRAm+99Rb+/fffJx47psvWrVsd8hymT5+OYsWK6edQqFAhjB8/3ub77tq1Cw0bNtTPI1WqVKhfvz727NkT7W3leQ8ZMgRFixbVfytLlixo0qQJzpw5E3Gbjh07xvqco3uNn+b333/X9128ePET2/PCCy/A19cXM2bMiHTdyZMn0aNHDxQuXBjJkyfXl+LFi6N79+7Yt29fpNt+9tlnsW7zhQsXYK/w8HCMGDEC+fLl069d6dKlMW/ePJvvv27dOjz77LP6eaRLlw6tWrXCf//9F+1tb968ib59++q/5e/vjxw5cujb37lzJ+I258+fx4cffojnn39ev+/yPOV1JiJKCEkS5FGJiFzUvXv3MGzYMJtOylu2bInVq1ejbdu26NKlCx48eKADipUrV6JatWr6xNvaF198oU/yoipYsKDd2z1lyhR07dpVb9P777+PP/74Az179tQnkf369Yv1vrt379Ynq7ly5cLAgQP1ye8333yDmjVrYvv27ShSpEjEbeU5ShDx119/6ecsJ8bXr1/Htm3bEBwcjJw5c+rbvf3226hbt26kv6OU0tuYN29efZLrCLI9crL8888/Y9q0aXjjjTcirpP3oU2bNkiSJAnatWuHMmXK6OBD3qOlS5di0qRJOvDIkydPpMeU36dMmfKJvyWBo70+/vhj/fmS165ixYr48ccf8eqrr+oT+ldeeSXW+8rzad68OcqXL68fIyQkRAe28t79/fffyJQpU8Rt5b2Q90+CPQl05TN2+fJl/bmQz7gEJuLw4cMYPny4DkRLlSqFLVu22P0ciYhipIiIvMDMmTOV7PLKli2r/P391dmzZyNdX7NmTVWiRImIn7dv365v/+WXXz7xWA8fPlRXrlx54rF37NiRINt+584dlSFDBtWkSZNIv2/Xrp1KkSKFunbtWqz3b9y4sUqXLl2kbT537pxKmTKleumllyLddvjw4Spp0qRq27Ztcd7OP/74I8bXzBYbNmzQ91+0aJH++f79+6pFixbKx8dHTZ06NdJtjx07pp97sWLF9HOJ6sGDB2rs2LHq9OnTEb8bOHCgfvzLly+rhHDmzBn92nXv3j3id+Hh4eq5555TOXPm1J+b2BQvXlwVLFhQ3bt3L+J3e/bsUb6+vur999+PdNtu3bqptGnTqhMnTsT6mCEhIerq1av6e3ld5fnL60xElBCYCkVEXuWjjz5CWFiYHhGOzfHjx/XX6tWrP3Gdn58fMmTIgMSyYcMGXL16Fe+8806k30u6z+3bt7Fq1apY7y+j2DK7YL3N2bJl0yPeMkp+69Yt/TuZyZAR8hdffBGVKlXCw4cPI6XVPM3cuXP1yLyM0NtL/raM8MuIv8wwyAyANUk3kuc+c+ZM/VyiklkMmdGRWZrEItsqMyzW75O8Ht26ddMzC7HNFly7dg0HDx7Ur32yZMkifi+zMJL+Nn/+/Ijf3bhxQz9vmamQGTJJFZNZiuhI+lP69Okd9hyJiGLDwIKIvIqciL3++us6rebcuXMx3s6SPjNnzhx9kmsLSU+5cuVKpIsEBNYkrSjqbaK7WJ/QSxqMqFChQqTHeuaZZ3Tqj+X6mMhJZ2Bg4BO/l3QZOSk9cOCA/llObOU1kfQnOWlNkSKFvsjPEtzERk6oFy5cqFPEJBXKHvJ6S/rZsmXLMHHiRJ12FZUERJL+U7ly5Tg/vpzER3295WTdmi3vkVysT+jlfZDXSwIBaxKkWa6PieVxYnqf5H2xrAH5888/ERoaqp+/pInJ9XI/CYJjWjdDRJQYGFgQkdeRPHg5eZXc85hUqVJFj+hLACLrCmQUXtYlnD59Osb7yKyA5MFbX6KuNZCF4FFvE91FRuStF+DKLEnmzJkjPZaMbMssRGwBkpA1FLKAXGZqLCSgkHUTwrLQ+ujRo/rr6NGj9QJfWdchI+NyEisLv6Muhrb2yy+/6CBK1jrYSxYbywLuCRMm6NH+qGTtgTznqIvthQQI1if+d+/ejfb1iPp6y/ttzZb3SC7WC7PlfZKF7jJLYc0yoxLb+yT3kzUemzdvjvR7eU0l4Ivuferfvz+CgoLw/fff6wBMZtlq166tt4OIyBm4eJuIvE7+/Pnx2muvYerUqfokNrpUGjk5lJPlr776CrNnz9YnkHKR9KOXX35Zn3RHXewrJ3dSnciaBATWZAYkupPd6LbRQm5vnR5jTSoPPe3xJDVHTtDffPNNXUVIUp4GDx4ccQJqub8lJUqqDcnouiWNSE5WZXRcgh15LWJKg0qaNKl+bex18eJFncoU3UJ4S2AholuALdW99u7dG/HzyJEj8cEHH0S6zZIlS3R1LGsy0xC1OpMtpEqYhbyOUp0puvfIcn1MZOZJZmYk2JWAQRapy/OU90uCwOjeJ/mMrl+/PuJ1kKC1atWq+nMo7y8RUWJjYEFEXmnAgAH44Ycf9FoLWVcQHTlJlNkNuchJ+MaNG/VtJeVHTqKjnmRLykvUdKWooluz8TSS5mI5uYxKZhOiS5+xJpWaZGRbTrJnzZqlfyfbKSetX375ZcSJqeVxZBut1ybkzp1bVyaSSlHRkRNdWV/QoEEDh6w9kQBmzJgxOs1n7dq1T7xmsm7A8nejkoBPAiMJTtq3bx/t49eoUQMZM2aMdRuiVryyhbx+0a11kPfIcn1spKqYzLLI87esAZKywBIQTp48+Yn3qWnTppGCK5l1kWAspveJiCihMRWKiLySzAjIiafMWtiSOiKzGrKYeNOmTbp0pwQXtq69sCYlQSVX/mkX65Nm+duSxnTp0qVIjyXBhqTKZM+e/al/VwIIOdmWhdyS0rRjxw49cyEssyyWx5G0nKgkDUvWh0Rn+fLlek2II9KgLM9XZgzSpEmjS99az0AI+b3cxrI2xJqsuZCgID4BnDVb3iO5WM9CyDbJ76TsrjXL5+tp75PMSn377bc6ZUo+Z1IqVmbNZO2OzGhYyhbH930iIkpoDCyIyKtnLZ621iIqmamQxcyyWFlGl+NKehvICejTLpKCZVG2bFn9defOnZEeS36W4MBy/dNIwzWZeZB+BuLXX3/V60cs/Tjk9/L8omtuJye71n0UoqZ3ych5s2bN4MjAT06q5YRaZkIs6wosJOA4duyY7sOREGx5j+SyYMGCiPvI+yAB1qFDhyI9lmUti63vkwQMzz33nA74JKCU9S4SMFlmJ2TRvojr+0RElNCYCkVEXqtAgQJ61kLSZ6QKlOT1W8iJrKRCSRpQ1MXBUjZUTtLjcwIXnzUWssZBSoZK2dXGjRtH/F5+lopAcpJtYVm0LNttaZIWHTkhllkLCWDk5N2SYiSPLxWXpMmcJeCQE2VJr4muOpPMwEiAIlWcYvt78SGBjpTSrVevnr7IwmbLYnhJ45J1HbIWQdYZRB29jzprEFfxWWMhze3ee+89vchfFp5btkPSmGS7pWKW9SyGzETIZ1CCuZjI+yO3tW7oKIvPpQytpJ/Je21J65K0MUl5e/fdd+P1nImI7MXAgoi8mqyfkLUWknZifZIo6TdSCapRo0Z69FhO7GWEWNYoyKiwrAGIujBbunTLCXlUckJpCRTiu8Zi0KBBeuF469at9Qi+pDTJGg9JcbLuUyAntJ9//rkuDysLmYWk1Uj+vuTryxoIqRAl1Z6k0lOvXr0i/a0hQ4boE3UJZqQPhBg3bpz+G9IDJLoARWZ9YkuD+uyzz57YJlvJYmTpoi3rCSS4kOctz0HS0SSwkIBGTrQtnbflRF66bct1EjBZOoVbk4pT0S38lse3BCjxWWMhf6t37956LYvMaMnslKSJyTZLQGn9eZEF2vJZkm21lOeV91MWlssaENk+Cdgk5a5z586647o1qdwl2yszUBLwSZAyatQoPcsRtZKWZSH3P//8o7/K511K1lpm7YiIHCZB2u4REbmY2Lpjd+jQQV9n3Xn74sWLatiwYbojd7Zs2VSSJEl09+ratWurxYsXR/vYMV3kekeQ7tNFihRRyZIlUwUKFFCjR4/WnZ2tWbpLW3dXli7V9evXVxkzZtRdx4sWLaqGDh0aqcOztV27dqm6devqztapUqVSzZs3V0eOHIn2tlWqVFGZM2eOtat0nz59dPfsQ4cOxanztrUFCxboDtQVK1bU3aStn5t0oZaO1QEBASowMFA/v65du+qu1dG9NjFdHNGROiwsTA0ZMkTlyZNHv0/ymZo9e3aMn7mTJ09G/E66ndeoUUN/zuS5lClTRk2ePPmJ99hi3bp1+vWX26ZPn1699tpr6vz580/cLrbnTETkSD7yj+PCFCIiIjxRLUtSzRYtWsSXhojIgzGwICKiBCO9GGQtinSEjtqRmoiIPAsDCyIiIiIishvLzRIRERERkd0YWBARERERkd0YWBARERERkd0YWBARERERkXc3yAsPD9eNqqRbrI+Pj7M3h4iIiIjIo0hnips3byJ79uy68ajHBhYSVOTKlcvZm0FERERE5NGCgoKQM2dOzw0sZKbC8kRTp07t7M0hIiIiIvK4fkQykG857/bYwMKS/iRBBQMLIiIiIqKEYcuyAy7eJiIiIiIiuzGwICIiIiIiuzGwICIiIiIiuzGwICIiIiIiuzGwICIiIiIiuzGwICIiIiIiuzGwII82efJk5M2bV38lIiLicYgo4fgo6dPtxg070qRJg+DgYPaxoGhJUHHq1CnkyZMH//33H18lIiJKVDwOkbuLy/k2ZyzIo3344Yc6qJCvREREPA4RJRzOWBARERERUbQ4Y0FERERERImKqVBERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZEbKRHRERugo1fyZUxsCACMGzYMN1IT74SERElttgCBuvreLwiV+bUwCIsLAyffPIJ8uXLh8DAQBQoUACDBg2CGzcDJzfFRnpEROTM4CG2gMH6OkccrzjrQR7ZIG/IkCEYNWoUZs2ahRIlSmDnzp3o1KkTvvzyS/Ts2dOhDTuIiIiInEmCCgkQJDD477//Il1nmY2QgKFr1642X+fo7SBy2wZ5f/31F5o3b44mTZroD3mrVq1Qv359bN++3ZmbRS6KIyxEROTOx6fYZhskYJCT/OgCh9iuc/R2ELn1jMXUqVOxdu1aFC5cGHv37tWBhcxitGvX7onb37t3T1+sI6hcuXJxxsJLFM2dG2FBQXgmSxbMHzUKuHYNuH0buHULCA0FfH2BJEnMRSLq9OnNJWtWIF8+IFMmwMfH2U+DiIjc1d27wNmz+vLrrFnY+OOPaFKjBrZv2IB7wcFIlTo1unbvHvNxKHPmBD0OcSaCnD1j4dTAIjw8HB999BFGjBgBPz8/veZC0qD69+8f7e0/++wzfP7550/8nqlQrsFhU7XykTx+HNi6Fdi1Czh0CDh4EAgKsmv7HiRLhuNKwa9sWRRq3RooWxaoWBFImzbW+zl6CpqIiFycHIdOnMD6L7/EicWL0SBXLuSWgaxTp+x73MBAIH9+oHRpoEwZcxyqVAlIl84hm83jFXl1YDF//nz873//w8iRI/Uaiz179qB37956xqJDhw5P3J4zFq4t3iMllkDil1+AtWuBzZuBq1ejv22KFECOHOaSIQOQMqW5BARIpCoVAYAHDyTaNDMa8jjnzpkRpug+6jJyJDv3GjWA2rWBunXN37DzeXHnTkTkZk6c0Mehk1OmINX+/cgox5ToJE+uj0FnfXyw98wZ5C9TBkWfeQbw9398HHr40ByHrl9/fBw6cybm41CpUkDNmsDzzwP16pnjGpGLiNOaZuVEOXPmVBMmTIj0u0GDBqkiRYrYdP/g4GD5H6q/kvNNmjRJ5cmTR399qvBwpbZuVeqDD5TKn192tZEvyZIpVbWqUj17KjV1qlJ//KHUlSvx/7uhoWru55+rjpkyqa3NminVurVSBQs++Xf9/ZVq1EgeVKlLl+L+vB6R28tnU74SEZFr0fv13LnV4g8/VKpvX3UjU6YnjgehgDqfL5+akSqV+v3VV81x6PJlc/yKj3v3lDp2TKlVq5QaMkSpNm2UKlQo+uNfgwZKyfnRhQuOe75xPI4Rxed826mBRfr06dU333wT6XdDhgxRheQ/mg0YWLihI0eU6t9fqdy5I+9IkyZVqlYttbV5c9U8a1Y1Zdy4xDmJP3dOqQULlOrWTal8+SJvk5+fCTLmzFHq9u04PSx34kRELur4cTUudWp1PMoJ/X1AbfX314NPLeJ4HLKLBA+LFinVo4dSBQpEPg75+ipVv75Ss2YpdfNmvP8EB7vIKwKLDh06qBw5cqiVK1eqkydPqqVLl6qMGTOqvn372nR/BhZuQkZpZs9WqkaNyDvMFCmUeuUVpRYvjthhxmfnF/UkPt4n9TIKdeCAUkOHKlWhQuRtTZNGqV69lDp0yObtICIiF3H/vlLz5ilVu3bkYEJmB1q3Vr907qxK5Mrl/P23HIfkODNihFKVK0c+DqVOrVT37uY4ZQPrYxKPT+QVgUVISIjq1auXyp07twoICFD58+dXH3/8sbonJ6I2YGDh4mTKePBgpbJlizz6IrMAMjpz584Td7F15xfb7Rw2MvPvv0p9+umTMxl16ii1Zs0T0+EcESIicjFXr5rBohw5Hu/DfXxMqpHMVsdxNtpWDjuRP3pUqc8/fzJ1t2ZNpVaujDUti8ck8rrAwl4MLFzU2bNmdD8w8PFOUIKLL75QKijIIX8ith2mw0dmwsJMICFrMyQwsjynsmWVmjtXqYcPE+bvEhFRnFj2w98NG6bU++8rlTz543121qxKDRyo1KlTCf6qOvykXo5D69Yp9dJLJk3X8pxKllTq+++VevDgibvwmESOwsCCnBdQSI6oLIC27PTKl1fqhx/UlPHjHXrS7bQdphyQ3nvPpHFZnmPRombkS3b8rrStRERepkLOnGoUoO7KrIRlH12mjFLffaeLeNjLEbPqdpMBOil8kirV4+coa1NloCuG4xCRPRhYUOK6cUOpjz6KPEPx7LNmdOXRNK0rTsnateOX6XWZgUmfPvLB6+efn7ipKz53IiKPEhKiU1f1mgnLPrlKFaVWr45/FadouNT+/Pp1U10qY8bIMxgrVjj0ORMFxyEVyjexauCSB5I63RMmAAULSht105G0WjVg/Xpg0ybTE+JRh1FpLid9IORrdD0fpFeEfE1M0vRO+lPI1ziTTqqffGLqnn/2GZAqFbB3L9C4MdCwIfDPPxE3je25ExGRHaRnhBw7ChUCvvgCSe/fBypXNn2R/vrL7I8d2Ona0fvzmI5/Nh0XpbmrNBSW49DgwUCaNMCBA0CzZqYXxr59DtlGojhx5ziMayycaONGpUqVejxKIr1Hli+P1yiJs0aAHDpVLT02+vQxZXMtpWqleoeMKBERkeNt3mzWulmOQ7LAWaoMutFofUzHv3gdF69dU6pfP9MHw1Is5e23zQw7kR04Y0EJ5+JFoH170yF0/34zcj9pkhklad48XiNDMY0AJfRMRteuXXUnbflq99+VLuBffQUcPAi0aGFG0SZOBIoWBebNi77bKhERxSjGffHly0CnTkD16sCePWbkftw4s/9t2dKhMxQJLabjn60zI5Feo3TpZCoe+PdfoHVr0wV8yhSgSBHg++95HKLEodwYZywSkYwAzZypVLp0j8v1yUhIDN2wHTEj4KyZDIf83fXrzSyOZSStbl3dlElwMTcRUTz2xXIckoal1msK3nxTzRwxwmsLZMR6vJLMghIlIpeoPXzYGZtJbo6Lt8mxTp5Uql69xzuncuWU2rEjwU/QnXUC7rC/KxVIBg16XCVLyh6OH6/y5s7tOov/iIhcVKR98enTSjVp8vg4VLq0Ulu2uN6C6kT21OOV9AWT0ruW4ioBAUqNGhVRJp3IFgwsyDFkdGj6dKVSpny8Q5IdVDT1sqPDkXmrBkcyUvTogHimcGH1bPbsXjm6RkQU5+PQDz+YrtOyD5X1A9J4VTppP8JjjY0DhDJzbgnMqldX6sgRfhjJJlxjQfa7dAl48UXgzTeBW7eAZ581FSb69QOSJLFrDYOnsHkthlTN+u03U0EreXLkOHIEf4SEoGvy5Jg8aVK81pE4q5IWEVGiuXoVaNMGeO01ICTEVHuSNRUff4zJ06dH7AM9/VjjEHnzAmvXmjUXKVMCmzcDZcsC334b7+MQUXR8JBKBmwoJCUGaNGkQHByM1KlTO3tzPIfsfGRHLsFF0qSmjF2fPoCfn7O3zKXIjljK1coCOzmo2UTKAnbsCPzxh/5xRfLkeP3OHaSNy2PE928TEbkLGYyRQiHnz5vBrIEDZUVzxMAW94F2OHXKLH7fsEH/uDp5crS/cwepeDwhB5xvs48FPfbggamJ3aCBCSpKlAB27MDk1KmRt0ABjmY4op55/vxmZy7Bmp8fmt25gwN+fhglo3JxwN4YROSRpKKeBBHSB0mCCqmst3UrMGBApNly7gPtkCcPsG4dMHy4Hjxs9Og4NKZVK0e8g+TtlBtjVSgHCgpSqlq1x/mXXbsqdedOgiyMc+d8WIdu+9atShUoYF5v6X8xdqxb1V8nInKos2cjrUdTnTsrdfs2X+SEtHOnUoULP+6/NGKEUmFhfM0pEq6xoLiRtJxnnjFdSmWKa+FC05siMDBBRobs6njtZA7ddskX3r0bkFEimS3q1cvUHg8OdsSmEhG5jy1bzHFo40azBmDOHGDaNL0uzVG4Ni0a8prv2gW8+qqZLerb1/Riun7dYa87eRemQnkzGaOQAKJ2bZP6VKYM8Pff5uTWiqMXxrnzFLbDt90SyElzJ1nPsmQJULUqcOyYYx6fiMjVffutabp64YJJwZUBFznRdXAw4c6DWglKArnZs+XFAvz9gZ9+MgNf0miPKI64eNtb3bsHdO8OTJ9ufpYcf/k+RQpnb5n32r4deOkl4OxZ00F10SKgTh1nbxURUcK4fx947z3gm2/Mz7L/mzXLnOg6iPUibxkQkqBCvrKCVPQWf/QRqowYgZwye5EmDTBvHtCokcPeD3JPXLxNsTt3DqhVywQSvr5mAZfsPBhUOFelSnqxvB4pkmloWUQ/fryZWSIi8iQXL5oF2hJU+PiYghaLFzs0qIg6y8yytE/3wdy5KB8Whh0ycyFpuS+8AHz9NY9DZDOmQnlrHqtU2UibFvj5Z5NTKTt2cr5s2YDffwdef93ku/bsiYM1apiRvRgwb5iI3MrOnUCFCmZ9n6SDrlhhelNMmeLwfgoMJuJGArDkefLg76++Mn2swsOBDz4w5WlDQ2O8H49DFEG5MVaFiqN580zXUnnbS5QwHaG9rKqT2wgPV4PTplVhj6qjnC1YUJXLmTPa19zRVbuIiBLM4sXqgVTBA9S1LFmU+vffiKu4L3MxUqVw3DhTLUqORVWqKHXuXLTnAawe6dmCg4P1eYZ8fRrOWHgD2SXI6EPbtmbkWyo+yIyFdISOB+sFcBylSCA+PsgwdCjeyJwZ9wICkP3YMSw9cwaLBg3yqMXwRORFpEhF69ZI8uABfgJQI1kyoEiRiKu5L3MxPj6YnDQp2mXIgFCpziXnDZKyu3//EwvhWT2SIig3xhkLGzx8qNS77z6uC96rl/mdHaxHKjjClAgOHlQ3MmXS79/d5MmV2rw5Mf4qEZFjSF+EPn0ijkMHatRQ+XLn5qy3G7Ac42tmz65U0aLmPUyTRi1//32bMhfim+HAzAj3Pd9mYOHJpMHdiy8+DipGjXL4n+B//kR6PS9fVqpyZfM+BgYqtWKFg/8yEVECCA1Vqk2bx8ehYcPYCNRNj0kzvvpKbfP3N++jfF2yxK7BR54/eGZgwXKznurKFaBZM7NYW6abf/gBePllZ28VxaE0ovQOieT2bfMeyoJ7Pz9gyhSzuI6IyBVJdbsXXzRN76RPz4wZQPv2zt4qsuP4dPHUKfwYGIj6d++aoi8TJwLdusWrxG+sxztyKSw36+1OnACqVzdBhVR+WrfOpqCC6yWcL9Y8VSkHvHy5qc4hFaM6dzYlGlmOlohcQKRjyOnTwHPPmaBCKj+tXs2gws3JcSlLnjw4OXIk8NZb5tjzzjvAJ59EOg7ZWuKXa2o8E2csPM0//5ja4NLBNHduszMvXtymu7KRkJuQHbjsyL/80vwso0XS70JmMaKwdJtlQygiSmiWY8jz2bPjNxnNlmaf2bOb41Dp0nwDPO04JMVEBg40P7/xhplFT5LE2VtGTp6xYGDhSXbvBurXB65eBUqVAtasMTt1G1mfhFoqPnCK0oVNmKD7XOgdfMuWwNy5Ju3NCqeaiSixyDFk+RdfYOmtW0h+8yZQrJg5DskgF3mmadOkWYjpdyHp1wsXAtJcjzwKU6G80V9/Ac8/b4KKihVNk7VHQYWtKU7WU5aconQDPXqYnbgEE0uWAC+99EQDI76PRJRYuj7zDNaEhpqgomxZkwbloKCCqbouqksXYNkyICDANDps3hy4c8fZW0VOxBkLT/Dbb2akQBb3Sk7rypUmp/URjlp7uLVrzc5cgop69cw6DKk5TkSUWKSLdpMmgAQVVaqY9CdZ4+cgPI65wXlI06YmqKhVC/jpJyBlSmdvFTkIZyy8yapVQOPGJqiQNCiZdo6S/8ZRaw8n77scxGVxtyzUl8/DrVvO3ioi8hay32nQwAQVclIpgx0ODCoEj2MurnZt4JdfgFSpTMaEfB6Cg529VeQEnLFwZ4sXA6++Cjx4YEasFyxgbqM327wZaNTIHNyrVTNladOkcfZWEZEnk/SX1q2B+/fN/kfSMgMDnb1V5Czbt5ug4sYNk5Ytg53p00dczYIi7omLt73BnDnA66+bBVNt2wKzZpk64eTdduwwMxiyU69QwYwgWe3UiYgcZtEiM7j18KFZ4yUFJLhwl/7+26TlyprPMmXMjFamTPp1YUqbe2IqlKebN+9xUCEN0qT5HYMKEjJCtGEDkCEDsHOnmZ6+fJmvDRE5lsxMyKCWBBXt2sU6Y86F116mXDmTDpUlC7B3r0mPkxL4TGnzCr7O3gCKI6kCJJ1LJaiQBmlTp0bbv4C8mFRjsd6pS18TGTmKhfWBnycBRBQrKRDxyiumUedrr5kZ81j6F1jKl8vX2PY95EFKljRVwaQ65cGDQJ06epArtoZ55Bm4xsLdRojatDE7844dgenTAV/GhhSDw4cfjxSVL48Z7drhi3Hjom2WZz09LdjDhIiiJdV+pG+OrO2TNKjvv3/q4FZsefVMjfFwx48DNWuaZonSJFGqR8mMOrkVpkJ5wwjRt99i8tSpHOmhmBUpAqxfb3Jbd+9Guf79cS2GUUPriiusvkJE0ZKCEK1amaBCBrlkpsKGGfPYRqm5v/FwBQqYYCJrVmDfvsdrAMljccbCjUeIYhrpsbXqAqszeAnZmUvzxGvXsN3fH/uGDUPn3r2dvVVE5E6kEIRUH7x3zwQXstYvlvQnokgkHUpm0GXNX6VKZkF3lNL45Lo4Y+ElI0QxjfTEls8an9uRm5PpZ9mJp02LSvfuobOUh7SxMyrzn4kIv/4KtGhhgooXXzTVnxhUUFwUL24+R1KlUErSst+Sx2KCvouJdCInaSxSwk/qg0twMXt2pJ15TNPLtk4tcwra88QYCJQv/7h5kVSNkpMD6dT9FAw+ibycdNRu1szsL6Sz8vz5rEJIdg9y6b5Llk7d5FmUGwsODlbyFOSrp8iTJ49+Ti2yZlUqRQql5C1q3lyp+/edvWnkRp8f+RqtP/+M+Fz9GhiopowfH+vjTZo0ST+WfCUiL7Nzp1KpUpnjUOPGSoWG2rxP4L6DYrRt2+PPVd26St29y8+MB51vM7BwMbIzbpAtmwoNDDT/6erV0ztzW+7HE0Cy6XOwYYO64+OjP18rkidX6uFDvnBEFNk//yiVIYM5DtWsqdSdO7YNXjxi6+3IS23e/HjwtEULpR48iPEzw/Mb9wosmArlYrrWqYM14eHwv3sXqFoV3zZpgrxFijy1xretKSvMmfdsNtUIr1UL63v0wH0ATWUaunt32bUn5mYSkSs7efJx5+QKFQBZlxUYqK9iqi05RLVqpjCNNFWUqpdduuDDvn3tWjdKLkK5MY9LhTp9WqncuU0EX7q0Uteu2TzqY2tEz1EkirBokVK+vubz1r8/XxgiUurcOaUKFDD7heLFlbp8ma8KJZzly5Xy8zOft969lQoPf+ImnLFwPs5YuCMpwSYjRKdPA4UKAWvXAunS2Tw6ZGs3Sy7YpghSEGDKFPP90KHAyJF8cYi82bVrps+ANDXLl88chzJmdPZWkSeTEsYzZpjvx4wBvvzyiZuwW7d7cXofi7Nnz6Jfv35YvXo17ty5g4IFC2LmzJmoINOvDqyr69KkWUzt2sDffwO5cgF//gnkzu3srSJvMWIE0K+f+X7aNKBzZ2dvEREltps3gbp1TSnQbNnMcSh/fr4PlDjGjgUs/ZUmTDApuuQy3KaPxfXr11G9enUkTZpUBxYHDx7E119/jXTp0sFryFoKKbkmQYV0SJZSbAwqKDH17fs4sHj7bWDxYr7+RN7E0p9CggrpMyAzFQwqKDH16gUMHGi+79EDmDOHr7+bcmrbzOHDhyNXrlx6hsIin0y/eouHD4G2bc3IkESA0megSBFnbxV5I0mFkjQImbGQ7u5p0pjUPCLybOHhwOuvm75JKVMCq1cDJUs6e6vIG0lgIceh8eOBDh3MceiFF5y9VRRHNs1YpE+fPk6XDBky6BX8T7NixQqd8tS6dWtkzpwZ5cqVwzQ5sfEGkoH2zjvAjz+aqghSdaNcuVjvwopOlGB8fIBJk4DWrU2Xdxm93LGDLziRpx+HZKR44ULT9G7ZMqBSJZvvzmMSOfw4JOss2rcHwsKAl18Gtmzhi+yJayx8fX0xZswYnV/1NPJw77zzDg4cOID8T5lKDQgI0F/ff/99HVzs2LEDvXr10jurDhKtRnHv3j19sc75khkPt1xj8emnwKBB8uICixaZDttPIR2VJWCTxdyyUJvI4aTLu4wQSUqepObJTr1AAb7QRJ5oyBDg44/N9/PmAa+8Eqe785hECUIGt1q0AH7+2aTmSZfuokX5YjtRnNY021JmysfHR128eNHmslQpU6ZUx48ff+rtkiZNqqpWrRrpd++++66qUqVKtLcfOHCgLr0a9eJ25WYnTjSl1eQSh47GLLlGiSIkRKly5cznU8pOxuH/PhG5ienTHx+HxoyJ10PwmEQJ5tYtpSpVMp9PKbcvZZDJc8rNhoeH61QlW928efOpsxUiW7ZsKF68eKTfFStWDKel5Go0+vfvr6MlyyUoKAjuwjJlvPatt8zCJPHZZ1JHzebHYMk1ShSpUpmRIlnvJGUnmzQBbt3ii0/kKaQxmRyLhJQyl3SoeKQ48ZhECSZFCmDlSlN+X1LrGzWSYXO+4G7AqeVmX331VR0c/PHHHxG/e++997Bt2zb89ddfHlVuVnbSeU+dwi8A/C3VdySnXXIKiVzRkSOmO6p035WduqwHkjxsInJfcmytUwcIDQU6djQ9BKIch5jiRC7jxAmgalXg0iXzuZVBr2TJnL1VXickIcvNzpo1C6tWrYr4uW/fvkibNi2qVatm04JtaxJEbN26FUOGDMGxY8cwd+5cTJ06Fd09sH7xyPbtscLHxwQVsp5i4kQGFeTaChc2I0aBgaZSjATDNo5DcFEnkQs6eNCsoZKgonFjYOrUaI9DbKRKLkOyXySYkIplUrmsUydTySwaPO646YxFkSJFMGnSJNSuXRtbtmxB3bp1MXr0aKxcuRJJkiTB0qVL47QBcj9JcTp69KguNSsLubt06WLTfd1mxkJStqpUAc6dA2rUMGVlHy1cJ3KLtAlZSCc78wEDTNGBp+CIJ5GLOX/eHIck1Vi+/vqrSTchcgfSW0XScqVM/wcfACNHPnETHnfcdMZCUpekO7ZYvnw5WrZsibfeegtDhw6NlNJkqxdeeAH79+9HaGgoDh06ZHNQ4VZdtSWNRIKKEiVMOgmDCnIn0sDRkms9ePDj72MZKeKIJ5ELkTVSMlMhQYXkrMtgAYMKcif165u0PfHVV6YsbRQ87riGOAcWKVOmxFXJudYB5FrUe9RES0rH3pUu0hS5dKekPf3zD5A9u5nOS5uWrxC5Hwn4H3VFDevWDWu6do0UTAwbNkynQspXwUWdRC5CRniljOzu3UDGjCatUb7GAVNMyCW89ppp5ipFhd57D+uiDETzuOOmgYUEEp07d9aXI0eOoLHkaULOnf/RJxn0iGSYvfkmsGGDyQ2UdSm5c0f78nCnTW5h4EDMS5kSfgBqTZ2KtZ99FhFMcKSIyIUb4MnxR2bKZaYiHn1pog4cEDlNv374LlUqffJa49tvzTkWuXdgMXHiRFStWhWXL1/GkiVLdJdtsWvXLrRt2zYhttE9ffIJMHs24OcHLF4MlC0b40250ya34OODkGHD8GtgIAKUwtxbt1A7e3YdVHCkiMh5Yhyc+vpr4JtvzALtOXPM2op44MABuQwfH9wbNgyrkic3xXBk/d/+/c7eKrJma3OM6dOnq8uXLyt3bdiRqKZMedx4aMaMp96cTYbI7RoXVa78uHHR+fPO3iIir5YnTx59LJSvERYufHwcGjXKmZtH5Hh37yr13HPm850jh1JBQbHenOdZiXe+bXNVKKkCJb0lypcvj+bNm+tLUSe3WHfJqlCyjqJZMyAszOSkSxM8Ik9z+TJQvTpw9CjwzDPAxo1cDErkJJZ1TpbZQ2zebGr+37sHvPsuMHYsy5uT57l+HXj2WVNGuXRpQAoIxXAuyIpRLlgV6rfffsP58+fxzjvv6LSnSpUqoVChQujTpw82bdqku3N7vV27gJdfNkGFNB56tNg1OlxXQW4tUyYTREsqpHzuX33VfO6JKNFFSkWUYL95cxNUyCDX6NEMKsgzpUtnjkNZswL79pnzrwcPor0p0/ncoPP2/fv3dbCxYsUK/PTTT7oilCzkbtasGRo1aoQUiVDKzqVmLP77z+SvXrwoK9z1YrnJ06dHHkWywuiZPKaLb+3a5iSmZ08zMkpEznHliulSfOwYUKEC8PvvnEkkjyaDtKs+/xxLr11DUqnEKZWipkxhMO1OfSwskiVLhoYNG+Kbb77RvS3WrFmjT5YHDRqEUaNGwatcu2Z6VUhQIdNxslg7adJYF2UzeiaPUK0a8MMP5vtx4xhYEDmLlHuXGQoJKqRC48qVTwQVUWfKOXNO7k7Or1ZeuIBuadKYYGLaNGDECH0dP99uNGMhEYtELjF58OABkiZNCq+YsZCRWmncsmkTkDMnsHUrkCNH9HmvRJ5KuqD27Wt27EuXmkodRJQ4JBVZ0kCWLDG9kmQmsVixJ24WdaacM+fk7iKdZ0kalMyci/nzkbdfv0ifd3KxGYsBAwboN2f48OHImTMnPv300xhvmxhBhUuQmEzWUkhQIS+05Po9CioES3CS1/jgA/nAm/8Tst5i+3ZnbxGR95CgXoIKOfYuXx5tUBHdTDlnzsndRTrPkkIFvXubKzp0wJhWrSJ93i04k+EiMxaS9iQN8TZv3gx/f3/UqFEDBw4cgDO5xIyF5JT362eCCsk1J/LmDr+SiiGdfTNnNrN3+fI5e6uIPNvEiUCPHuZ76Z3Urp2zt4jIeaSISKtWJsBOn94chwoVinQTztS5yIxFkiRJULFiRWTLlg3p06dHxowZ47FpHki6mh4/zqCCKEkSYMEC0wzy0iWgcWNTDpCIEoZ00rakfgwezKCCSJoSSzPIihXN+lc5DklRAyucqUtYNgcWzzzzDBbIScOjilByoUes0p+IvFqqVGbRqKw3+vdf4KWXzDokInKsnTuBV14x6ys6dwY++siuh2N6CHmM5MlN0C1FDKSYgZRfDg2NuJpp6i5abvbhw4d6FgPengpFRE+SmuLSuOjmTRyuXBkNzp/Hh/37s4gBkaPLmzdoYE6i7FzbyPQQ8jiHDpnKhTdumOIG8+YBvvEuhurVQhKj3KwEFbdu3dJ/zPpCRBRRdtnPD0W2bUPH06ejLbtMRHEk6YWS3mEpb75wod1BhWB6CHkcKWKwbJn5/yH/T+yc1SPbxDmwOHnyJJo0aaIb4En0ki5dOn1Jmzat/uptOH1MFAMpw/yoXv5nAGbUrMn/L0T2kLRCSS+UkVhJwV21ylQkdACmh5BHqlULmD7dfD98uGmeR66VClW9enXIXXr16oUsWbLAR+rWW6lZsya8KRUq6vQxe1cQRSGjREOHyjQnXk2fHvMuXWJdcaK4kkP166+byk+ylumPP4AyZZ56Nx6TiAB88QUwcKCeRdepg9LUmBLkfDvOgUXKlCmxa9cuFClSBM7mCoFF1J0281SJopDFpVICc/583AsMxAtp06Llp59yvQVRXEjvqEGDzImRzFTI2gob8JhE9Cgw79QJmDVLTmRNYC4VDMn5ayyk5GxQUFBc7+axok4fM0+VKApZLDdzpl7M7X/3LtYlTYqucezMzZRD8mozZpigQkgqh41BheAxiUiG0X2AqVNNa4Bbt4AmTYAzZ/RLw+OLY8V5xuL48eP6JLp9+/YoWbLkE122S8tiMi+asSAiG129ClStChw9KvWrgY0bgRQpor0pZwKJHlm3zizWlgaUH39s+lUwxYkofqRCVPXqwMGDpvjBH38gb+nSkVLaKZFTobZu3YpXX3010osv6yzkYeRrmHQ9TCQMLIjcjNQUl+BCGhY1bWoqdkhqRxRcu0QELBwwAI2GDEEqOUy/+qpZX/FoXSNTnIji6dQpU675wgU9+ze1aVMMGTkyIqWdEjmwKF68OIoVK4a+fftGu3hbIr7EwsCCyA1t2QI8/7ypcNOjBzBuXMTJkgUXnJLXO3sW5/LkQfawMGz190eV4GDA35//R4gc1WBSig3duWMaTEqaVJTjECVSYCFlZvfu3YuCBQvC2RhYELkp6XHRurX5fvRooHdvZ28Rkeu4eRN47jlg714cS5IEm4YNwxt9+jh7q4g8y4oVgKz3k9PgIUOA/v2dvUXeuXi7du3aOrAgIoq3Vq2AkSPN9++/DyxdyheTSDx4YLoEy3E2c2YUPHIkzkEFF6MS2aBZM2Ds2Mdl0efO5cvmAHGesZg6dSoGDx6MN954A6VKlXpi8XYzeaMSCWcsiNyY7HokFeqbb4CAAOC338z6CyJv/j/x9tvAtGlAYKApcFCxYpwfhusviOJABrdk5jxZMmDtWpMiRYmXCuUrpSNjejAu3iaiuJBKN9JJWBoWZcxo1l+4QJolkVNIOoZUfpJcbyls0Lx5vB6Ga5SI4thrSWYJlywB0qYF/voLKFaML2FipUKFh4fHeEnMilBE5P4mf/stiu3Zg0tS9EEqRUk3VPlK5G2k4pMEFUIKGsQzqIiuvxIRxUIGzH/4wcyYSzlaOQ5JxahYMN3QgTMWroSpUETuzZKyUSFnTuxIkgSQMtayc1+/3qSCEHkDSQNs2NCsr/jgg8frj4go8ciglhx/pCz6U3oteVu6YYijZyzGjRuH0NBQmzdAIrmbUtWCiMiGrsBvykjtzz8D6dKZdKj27QHOgJI3OHAAePHFx4u2hw939hYReaeMGTH3tddwVWYwdu0CXnnFpOtGgx3t7Zyx8PPzw4ULF5ApUybYQqKZPXv2IH/+/EhInLEg8jCbNgH16gH37wPvvQeMGhXrzZlLTm7t7FnTqOvMGVNeVhaOSiEDInIKmYnIduoUNvj4IEBOj7t1AyZO9PoeFyGOXrwtC7ZLliyJJJKqYIP9+/fj8OHDDCyIKO7mzwfatjXfSynAnj1jvKm3TUeTBwkJMcHEvn1A0aLA5s1A+vTO3ioir2YZrJrSoAEaSHU2OUUeMQL43//gzUIcHVh8/vnncd6IXr16Ia2srk9AnLEg8hyRZh+ky/CHH5pRIulxIU2MnnYfLlQldyFpT02aAOvWAVmyAFu3SpTs7K0iImtjxpiZc8uAV5s2Xvv6hCRkuVlXwsCCyHNEmn04eRJ45x2JHExqyIYNJmWEyN3JIbdTJ2DWLLMwVBaIykJRInI9vXubmXPpcfHrr2aW0QuFJGS5WSKihBBpMZzMVIwfb0Z1pXBE06bA8eN84cn9SQaABBV+fsDChTqoYOlKIhf19demuML9+wht2BB1cuTQ/18pZpyxICLXdesWUKuWqdAhjfM2b8bkpUuZ/kTuacYM4M03zfdTpgBvvaW/5VohIhd25w5QuzawbRtkJd/LOXJguxRc8CIhnLEgIo+QMiWwcqXJP5fa4i+8gHFDhuiUKVlbQeQ2fvklIpDARx89/p6lK4lcW/LkwE8/4UbmzJCVUD/L76SR3iOccYyMMxZE5PqOHAGqVQOuXsXp4sVR+9YtfNC/Pxdsk3vYvRuoWdPMwEmPlu+/9/rylURuR9b+yXFIunLXqGEGCwICvGLGMSQxZiyuXLmi/xARUYIrXBhYtUqPHOU+eBDHatVC17ff5gtPrk9m2ho1MkGFpFNMn47JU6bokxHmahO5ridmIvLlA1avlmZtpudSu3a6kSub5dkRWNy4cQPdu3dHxowZkSVLFqRLlw5Zs2ZF//79cUdy0IiIEkrlymaxqyx6lRFf6dYdBaekyaXIyGaDBsClS7icKxdKHj2KyTNm6DQ+pvMRubZo/5+WLQv8+KOpEiWl0Hv00INcMlPBkudxDCyuXbuGypUrY9asWWjZsiW+/vprfWnWrBnGjx+PGjVqIDQ0FNu3b8e4ceNsfVgiIttJlShpWiSGDjWVo6zwhI1chszoy0zFiRNA/vxopBT+CQqKKDwQUQEtDhg4EyWeGP+fSkGROXNMOqPMZgwaxLfFmrJRr169VMmSJdWFCxeeuO78+fOqVKlSqlWrVip16tTqu+++U4khODhYenDor0TkRQYPlm4ASvn4KLVwYcSvJ02apPLkyaO/xsbW2xHFS2ioUrVrm89o5sxKHT3qkM+c3F+OefKViJxs4kTzf1wuU6YoTxYch/NtmxdvS57ZlClT0ECmdaOxZs0aNG7cGAMHDtSXxMAGeUReSnZbPXoA33xjpqTXrAGef97mu3vDYjtykvBwoG1bk7YnVc2kAV758g55aHaaJ3Ixn35qZix8fYElS4AWLeCJEmTx9vnz51GiRIkYry9ZsiR8fX3jHVTI9LCPjw96S5dDIqLYyBS0pFy+9JJuXIRmzYCdO21+zbjYjhIs4O3VywQVSZMCy5Y5LKgQksPNXG4iF2t42aWLGVB45RWzqNvL2RxYyILt2Eb2Tp48icyZM8drI3bs2KFnQ0qXLh2v+xORF5JF3JLnKjMVUnGnYUPg0CGb7soTNEoQQ4YAEyaYwPeHH4C6deP8EFxHQeRG5P+6zJw3bw7cu2cGufbv9+r/zzYHFpIC9fHHH+O+jA5Gce/ePXzyySdoKAf2OLp16xbatWuHadOm6SpTREQ2CwgwFToqVtQ9LlC/PnDqFF9ASnwTJwIDBpjvx4wB2rSJ18OwAAGRm0mSBJg3D3j2WSA42AxyPToOeeP/Z5sDiy+++AKHDx9GoUKFMGLECKxYsQI//vijfrHkd4cOHcJnn30W5w2Q8rVNmjRB3XiM7BARIVUq4OefgWLFgDNngHr1dHlPIkd56qijzE7Imh9LznXPnvH+W0zTI3JDgYHAihWALBk4d84chy5ciPT/2WtmL+KyKvzEiROqYcOGytfXV/n4+OiLfN+gQQN19OjROK8ynzdvnq40dffuXf1zzZo1dfWpmISGhuoV6ZZLUFAQq0IRkREUJGVzTIWOcuWUunEj4pVhFSiyR6zVmJYuVcrPz3zu5PgVHs4Xm8hbnTnz+DhUqpRSV696RFW3uFSFilODvHz58mH16tW66/bWrVv15fLly7oiVMGCBeMU0AQFBaFXr16YM2cOAiSdwQZDhw7Vq9Itl1y5csXpbxKRB8uZE1i3DpC1Xn//DTRtCty9G6fpaK8ZUSLHzCL8+qtZsBkWBnTsCIwaZXKuicg75cgBrF8PZMtm1lpIL5ubN71qNtLmcrOOtnz5crz44ovwkwWYj4SFhenKUFJdStZtWF8n5HdysS5/JcGFLeWviMhL7NkD1KxpGpRJQ72lSyO6HcsOPbbuqCxDSzbbssUszr5zB2jZEpg/3+RaExH98485DsnaP/kq6brJk3tFuVmnBRY3b97UI4jWOnXqhKJFi6Jfv366fO3TsI8FEUXrzz/NQm6ZsYjDSR/7BJBN9u413Xdv3DCfM8mt9vfni0dEj+3aBdSubQa5ZOZi+XLTd8kNuUVgEZ1atWqhbNmyGCMVNWzAwIKIYrKyZ0/UGz8e+nRPGpbJAtsos6BEcXbkCPDcc6ZAQPXqwC+/AClS8IUkotgHuVq1MtWj3HBmM0Ea5BERuZMeK1agFYAH8oPszDt3Nk2MiOLr6FHTN0WCirJlgZUrGVQQUcxr9Z599vFMxeLFXnEccqkZi7jijAURPS2taUr9+mgwY4ZZYPvWW3IFF9hS3B0/bnKlz541JSU3bAAyZeIrSURPX6u3fLmZsbAchyZNAnzdZ2yfMxZE5PUs3bUbTJ1q0qBkJy7f9+olhQC9/vWhODh50sxUSFAh/VKk6guDCiKKxofRVX9q0QL4/vvHx6F33vHYmQv3CZeIiOJL1ljIrIUYPx7o25fBBdlGRhwlqAgKAooWBX77DciSJdqbslwxEXV9NKj1RAXCV18FZs0yM+ZTpkQKLjxp38HAgoi8Q4cOZmcuvvoK+PhjBhcUu9OnTVUXqWBYuLAJKrJmjXQT6xMCW/ulEJGXat/ezFxYgovu3XVw4Un7DgYWROQ9JLdVZizE0KF65mLypEkeM1JEDnTmjJmpkDQoaQArQYU0vYrC+oTAWxpgEZGdwcWsRzMXctzp3h0f9u3rMfsOLt4mIu8jwUXPnvrb71Klwhs3byK39UI78m6S9iQzFceOAfnzAxs3ms7u0WDvEyKKl++/Bzp2NDPn3boBEya47IJut+1jEVesCkVE8SYL6CQHVinMTZkSIcOHo6vkvJJ3k2BCOmpL+lPevCaoyJ3b2VtFRJ4eXHTtCkyc6JLBBatCERHZkhY1c6beib966xa6bt9uSgHGwJMW11EMDh4EatR4vKZi0yYGFUSUICbLMeXTT/Hba689Tovq1AlTJk5062MNZyyIyLtJ8zzZsUtQ8corZgQpaVLbapOT59i923TIvXoVKFUKWLcuxupPtmKaFBHFJNIx5csvTYGRsDD8EhiIZnfvIpsLHWs4Y0FEFJdStAsXmmBi/nwTXNy//8TNuDDXg/31l1moLUFFxYrA77/bHVQIT6r0QkSO9aF1sYd27YAlS3SH7gZ37+LXgAAM6N3bLV9yzlgQEYmVK4GWLU1QISPXspNPmZKvjaeTZnfNmgF37pg0qJ9+Ap6yONFWnLEgojjvj5o3B27fBqpWBVatAtKlg7NxxoKIKK5eeMGcVCZPDqxda0awL13i6+jJ5P1u0sQEFQ0aAKtX2xRU2LreJsZGWURE0alTxwQXEkxs2QLUqgVcvAh3whkLIiJrsohbTjavXDH9C375xZQcJc8ybZop8Shra1580ay18fe36a5cb0NECWrfPjNzLkFFoULAr786tZAEZyyIiOKrUiVg82ZTalRKj1arBvz9N19PTyFlHT/5xFQFk6BCSj3KGhsbgwrB9TZElKBKlwb+/BPIkwc4ehR47jmTHuUGOGNBRBSd8+eBhg3NyFGqVMDy5aZpGrkvWT/TpYup/CUGDjQXKfVIRORqzpwB6tUDpMfSu+86bTPYIC+KsLAwPHjwwDnvBiWIZMmSwdcFm8iQhwkOBlq0MFWCpGrUDz8Abdo4e6soPkJCzOJ8SSnw8wOmTAHefDPGm3PhNRG5hLt3gcBAp24CA4tHpKn4hQsXcOPGDae+IeR4ElTky5dPBxhECSo01PS5WLzY/Dx4MPDRR5g8ZYouIyppMVyc6+LOngUaNzazTylSAIsWAY0axXoXrqMgIjIYWDxy/vx5HVRkzpwZyZMnhw+nuz1CeHg4zp07h6RJkyJ37tx8XynhSS7+++8D48aZn195BUU2b8aRoCA2zHMRMc4w7N9vFuMHBZneFFK+8Zln4v94REReJiQkBGnSpEFwcDBSP6VynseusZD0pyNHjuigIkOGDE7bRkoY8p5LcFGwYEEdYBAlCkmf6dEDePgQl/LkQdOHD9FpwACeeLqAaGcYZFF2p06mnGyRIsCaNWZRPhER2YxVoYCINRUyU0Gex5ICJQEkUaJ5+21g3TogQwZkPnUK28LD0bV8eb4BLiBSpSbZL8hXWQ8jQYUsfpTu2k8JKmztT0FERNHz+NWvTH/yTHxfyWmkYZH0uihRwlSOkm7Nc+bwDXGyiGZ0L79sUp+GDzdX/O9/wM8/A+nTP/UxJPVJZj3kKxGRK5jsZgMeHh9YeDL5oI0ZMwaewtOeD3kwaZgnI+BNmwL37gHt2wN9++oUKXIiWU9RsaJpaihVVKTp3YgRQJIkNt2d/SmIyNUMc7MBDwYWLiooKAhvvPEGsmfPrtN+ZIq/V69euHr1qrM3jYiErOtatsyk3IiRI4E6dfD98OHxGl1yt1EplyNVu6pWBU6cMClPW7boRfbxmvXgYm0ichEfWqd5ugEGFi7oxIkTqFChAo4ePYp58+bh2LFj+mRj/fr1qFq1Kq5du+aU7ZL1DFKRiYgekX4IQ4diXefOuCVV5zZtQqOPPkL+eIwuuduolMuQGaOePYHWrU1n2jp1gJ07gTJlnL1lRER2c7cBDwYWLqh79+56lmLt2rWoWbOmLqnaqFEj/Prrrzh79iw+/vjjiNvevHkTbdu2RYoUKZAjRw5MnDgx4jop+PXZZ5/p+/v7++vZj55yAH7k3r17+OCDD/T95P6VK1fG79II7JHvvvsOadOmxYoVK1C8eHH9GN9++y0CAgKe6A0isym1rboS//nnn3juuecQGBiIXLly6b9726od/aVLl9C0aVN9vfSjmMMcdXJjXdatQ3mlcChpUmQKD8c6AEuKF5cqEh47KuUSjh41sxTjxz9eTyGVn1gJkIjIKRhYuBiZjfjll1/wzjvv6JNua1mzZkW7du2wYMECHTSIkSNHokyZMvj777/1CYmc4K+TqjVyYrNkCUaPHo0pU6bo2Y/ly5ejVKlSEY/Xo0cPbNmyBfPnz8e+ffvQunVrNGzYUN/W4s6dOxg+fLgOKP755x/99yXYkMe2nsmQbZLrxPHjx/XjtGzZUj+uXCeBhvw9i44dO+p0rw0bNmDx4sX45ptvdLBB5I7k/979PHmw+auvgDfegB+AZ1avNgu7JTXHA0elnEr2f999B0hFrr//NoHEypXRrqeImmLGlDMiogSk3FhwcLCcXeuvUd29e1cdPHhQf3WESZMmqTx58uivCWnr1q36OS1btiza60eNGqWvv3jxot6ehg0bRrq+TZs2qlGjRvr7r7/+WhUuXFjdv3//icc5deqU8vPzU2fPno30+zp16qj+/fvr72fOnKn/1p49eyLdplevXqp27doRP//yyy/K399fXb9+Xf/85ptvqrfeeivSff744w/l6+ur34/Dhw/rx92+fXvE9YcOHdK/Gz16tE2vk6PfXyKHmj9fqTRp5PRXqVSplJoxQ6nwcL7IjnDpklIvvmheW7nUrKnUmTMx3lz2k7Jvka/R/UxERPE/346KMxYumv9sa99CWXMR9edDhw7p72UG4u7du8ifPz+6dOmCZcuW4eGjqjX79+/XMw2FCxdGypQpIy4bN27UMw4WkpJVunTpSH9DZiYkZUoa1AlJY2rSpImeyRB79+7VaVTWj9ugQQO9PuPkyZN6+5IkSYJnrLrfFi1aNOL+RG5P+ifs3Qs8+6zkK+pZDLzwAnD2rLO3zL39+CNQsqRZNC+NMYcMAdavB3LksDnFjClnREQJh4GFjRLrYCSdpKVHgyU4iEp+ny5dOmTKlOmpjyVrGw4fPqzTjCStStKratSooZsH3rp1C35+fti1axf27NkTcZHHHzt2bMRjyP2i9oyoWLEiChQooFOoJHCRgMWSBiXksd9+++1IjyvBhqRYyf2IvEKePICsWZJ+CtLQUXopyEmxpPAoxZScuLh40QRrLVrIAi3zOkovkf79zQL6OKSYMeWMiCjhMLCwUWIdjDJkyIB69erpYEBO2q1duHBBzw60adMm4mR/69atkW4jPxcrVixSYCCLpMeNG6dnGWRNhcxWlCtXTs9YyLoGCWasL7KW42kkkJBt+emnn+Dr66tnLCzKly+PgwcPPvG4cpEZEJmdkJkTCWosJACKuiCcyO3JSa/0t5B1ABUqAPIZ79RJVy5aOGgQq0DZupZC9mkLF5rXs18/YMcOoGzZGO/GdRRERM7BwMIFTZgwQVdskvShTZs26UXOa9as0QGHVHD68ssvI267efNmjBgxAkeOHNEVoRYtWqQXcAtJR5o+fToOHDigS9jOnj1bBxoy8yIpUBIcvP7661i6dKlOUdq+fTuGDh2KVatWPXUb5b67d+/W29KqVStdMcqiX79++Ouvv/RibZmtkJmKH3/8MWLxdpEiRfTibpnV2LZtmw4wOnfu/MRidSKPIRWitmzB1hYtcFcGBTZswLpLlzAmTRoMeP99Z2+da9qzB3juOROIXb8OlCtnZikkHTUgINa7snQvEZFzMLBwQYUKFcLOnTv12oiXX35Zpw+99dZbeP755/WMQ/r06SNu26dPH31bmYEYPHgwRo0apQMSIWsWpk2bhurVq+t1ElKuVmYYZFZEzJw5UwcW8hhyst+iRQvs2LFDl6d9Gpl9qFSpkq76ZJ0GJeRvyVoNCXak5Kxs26effqrL3VrI35afpZzuSy+9pJ9f5syZHfgqErmYJEnwyt9/o7hS2BAQAL+HD9ErOBidpZLU/PlmdJ6kNJ6UrANkDdbmzUDy5CadTIIKqQJlw6wE11EQETmHj6zghpsKCQlBmjRpEBwcjNTSBddKaGioHoWXHgnSd4E8C99fSgxy0iqj33Ki6og0yIjH69cPXWWdVJ8+wOnT5spq1QAJMqIUZPAakvo5bpxuOIjgYPM7WVchr0nOnE/cXIIKKaghM7CSpkpERIl/vh0VZyyIiBIppSZirVa3bkCrVsC//wKDBplR+b/+MsFF06YmDchbSBPB6dOBwoVlqsEEFVKJ7rffzExONEGF4KwEEZHrYWBBRBSDBD95lXVFAwYAR44Ab75pFidLozdZT9C6NWBV4MDjFinfuwdMmSK5n0DnzsCZM4CkYX7/vVns/vzzsd6d1Z2IiFwPU6HILTEVijySdL0fOBCYN+/x7+rVM5WQatcGopR+dst0IFlD8e23wPjxJpgQWbKY6lnvvPPUhdlERJS4mApFROSOZPR+7lzpYAm0b29mMNatA+rWBUqUkJJxj9cfuFs60L59gKSASWqTBEoSVEhBB+mbc/IkINWxGFQQEbk1pkIREbkaaQD3ww/AsWOmQlKKFNIdE3j3XdNlWkqwrl+Prl26JEp/nXiTMrGTJpkeHmXKSO6WWaQtayhmzgSOHwd69jQpYU7g1qlkREQuiIEFEZGrypvXpAydO2e+SqO427dN0ziZxZDu3jLSLx2+Hz50jZPrq1eBGTOAxo1NipOkN8lakaRJgZYtdQ8PvTi9Y8c4z1A4elvZ74KIyLG4xoLcEtdYkFeS6uDS22H2bNOJWmYELNKlAxo2NIuea9WSZjNPrMmIjt3rNGQRtvSYkJSttWtNV+zw8MfXlyplZlgktUtK7NrB0WtKHF1OmIjI29dYMLAgt8TAgjyZTSe8ckK/ejWwfDmwahVw5Urk67NlAypVMhWm5CLpVblymZmDuP4tC1nfISlZUiZ3925g2zZTwUlKxlqTVCepaiUzFDLL4iAMBIiIEh8DC554ejwGFuTJ4jwyHxaGljlyoOzFi6jv74/KMrNx//4TNwv39YWvlHSVFKuMGYEMGcxFUpKSJDGXsDCTbiWXkBDg/HlzOXsWuHQp+r8vj1WnDlC/vqliJQEMERF5XWCRJNG2ihzGx8cHy5YtQ4sWLRL85KZ37976kti+++47/Xdv3LiR6H+byJHiM8out7XcxyZ+fqj32Wf6Ppk+/BCVO3Qw6Ukym7BnDw7NnYu8Dx4gUFKUJFCxJ41IKjkVLWpmQCpXNpf8+W1KuyIiIs/m1MBi6NChWLp0Kf79918EBgaiWrVqGD58OIoUKQJvdvnyZXz66adYtWoVLl68iHTp0qFMmTL6d9WrV8f58+f171wNgwGi2BcI2xpYyO1ium1MgcoT96lZ01wAbKxSBY2HDsXn77yD1599Fjh1yiyylov0lZC0Kln8LRdfX1OFSrqBp0plUqokmJCvBQoATxmtshfTnYiI3JhyogYNGqiZM2eqAwcOqD179qjGjRur3Llzq1u3btl0/+DgYCVPQb5GdffuXXXw4EH91d0899xzqnLlyuq3335T//33n9q2bZsaMmSI+vHHHxN1O/LkyaNGjx5t8+3lvUyTJo1D/vbTHsud31/yLpMmTdL/l+SrI+4v38t+T766gqjbZ+/zdbXnR0Tk7YJjOd+OyqmBRVSXLl3SG75x40avDSyuX7+un9Pvv/8e423k+mXLlunvT548qX9esGCBevbZZ1VAQICqUKGCOnz4sNq+fbt65plnVIoUKVTDhg3162tRs2ZN1atXr0iP27x5c9WhQ4cYA4uvv/5alSxZUiVPnlzlzJlTdevWTd28eVNft2HDBr0d1peBAwfq60JDQ1WfPn1U9uzZ9X0rVaqkbx81kMiVK5cKDAxULVq0UF999RUDC6JoTrQdHag4evusf7b1b1nfztGBChEReWlgcfToUb3h+/fvT5jAIjxcKZkNSeyL/F0bPXjwQKVMmVL17t1bn5DbGlgULVpUrVmzRj/nKlWq6ICiVq1a6s8//1S7d+9WBQsWVF27drUrsJDvZRZF/ub69etVkSJFdHAh7t27p8aMGaNSp06tzp8/ry+WoKNz586qWrVqatOmTerYsWNq5MiRyt/fXx05ckRfv3XrVuXr66uGDx+uA6KxY8eqtGnTMrAgSoRAwF6xBQK2BhmxbRNnMIiInMstA4uwsDDVpEkTVb169RhvIyfa8qQsl6CgoLgFFnKSbyrBJ+7FxtQui8WLF6t06dLp2Qc5Ie/fv7/au3dvrIHFt99+G3H9vHnz9O/k5N9i6NChOhCwJ7CIatGiRSpDhgyxpi+dOnVK+fn5qbNnz0b6fZ06dfTzEm3bttVpcNbatGnDwIIoASTmDEBMQUZctokzFkRE7hNYuEzn7e7du+PAgQOYP39+rIu9pdyV5ZLLQ0satmzZEufOncOKFSvQsGFD/P777yhfvrxeHB2T0lI3/pEs0u1W96UqFel3l2IqFWmjX3/9FXXq1EGOHDmQKlUqvPbaa7h69Sru3LkT433279+PsLAwFC5cGClTpoy4bNy4EcePH9e3OXToECpLZRkrVatWtWtbibxdTF2qZYG3lLBNjIZw1n9LFptL+dzoKl3Ftk2Jub1ERGQflwgsevTogZUrV2LDhg3ImTNnjLfr37+/rqFruQQFBcXtD0mVk1u3Ev8ifzeOAgICUK9ePXzyySf466+/0LFjRwwcODDG2ye1anol5Wij+124VTdcX19fma2K9BgPoja5siIH9hdeeEEHMEuWLMGuXbswceJEfd39aOrlW9y6dQt+fn769nv27Im4SDAxduzYp74ORGR/NSpHBibxxQCBiMjzOTWwkBNbCSqkJ8Nvv/2GfPnyxXp7f39/3ZjD+hIncsItZRQT++KA+u7FixfHbWlY5SCZMmXSZWstZFZBZoxiIoGBBCZff/01qlSpomcgZFbFWrJkyfTjWCtXrpz+ncyWFCxYMNIla9as+jbFihXDNunga2Xr1q0OeqZE3im2GYL4BCaODjSIiMjz+Do7/Wn27NmYO3euTq25cOGCvty9exfeSlKLateurV+Xffv24eTJk1i0aBFGjBiB5s2bO+zvyN+QPhlykT4i3bp1i7UZnQQCMqMxfvx4nDhxAj/88MMTJxhy0iEzFOvXr8eVK1d0ipQEIO3atcPrr7+ue5bI89m+fbtOa5O/LXr27Ik1a9bgq6++wtGjRzFhwgT9MxE5b4YgamBi7wwIERF5PqcGFpMmTdIpTbVq1UK2bNkiLgsWLIC3kvUHst5g9OjRqFGjBkqWLKnTobp06aJPuB3ljTfeQIcOHfQJf82aNZE/f348//zzMd5eGvSNGjVKNzCUbZozZ44ODqxJg0M5iWnTpo2eEZFgSMycOVP/nT59+ujmh9IxfMeOHcidO7e+XmZApk2bplOj5O+sXbsWAwYMcNhzJSLbWM9KRA1M7J0BISIiz+cjK7jhpkJCQvQibglOoqZFhYaG6tFxSa+S9QrkWfj+EjmeBBUyKyEBhAQVREREIbGcb7vk4m0iInK+xJyV4JoNIiLPwxkLckucsSByb5wdISJyD5yxICIil8Y1G0REnieJszeAiIi8jywKZ9M7IiLPwjUWREQUJ1wfQUREXhlYuHHRK4oF31ci52FPCyIi8qrAImnSpPqrNGkjz3P//n391c/Pz9mbQuR1uD6CiIi8qiqUOH/+vO4mnTlzZiRPnhw+Pj5O2U5yrPDwcJw7d04Hj9Jkj+8rERERkfOrQnn04u2sWbPqr5cuXXL2ppCD+fr6MqggIiIiciEeHVjISHa2bNn0jMWDBw+cvTnkQMmSJdPBBRERERG5Bo8OLCwkD5+5+ERERERECYdDvkREREREZDcGFkREREREZDcGFkRERERE5N1rLCyVcqUMFhEREREROZblPNuWDhVuHVjcvHlTf82VK5ezN4WIiIiIyGPJebf0s/DYBnmWRmmpUqVyWpM0ieIksAkKCnpq0xDyDvxMED8XxH0F8fhBnnJOIaGCBBXZs2d/aql/t56xkCeXM2dOuAJ5sxlYED8TxH0F8fhBPKcgTzvPfNpMhQUXbxMRERERkd0YWBARERERkd0YWNjJ398fAwcO1F+J+Jkg7iuIxw/iOQV563mmWy/eJiIiIiIi18AZCyIiIiIishsDCyIiIiIishsDCyIiIiIishsDCztNnDgRefPmRUBAACpXrozt27fb/66QWxo6dCgqVqyoGzZmzpwZLVq0wOHDh529WeRChg0bppt59u7d29mbQk509uxZtG/fHhkyZEBgYCBKlSqFnTt38j3xUmFhYfjkk0+QL18+/XkoUKAABg0apJuSkXfYtGkTmjZtqhvQyTFi+fLlka6Xz8Knn36KbNmy6c9I3bp1cfToUbgiBhZ2WLBgAd5//329Wn/37t0oU6YMGjRogEuXLjnuHSK3sXHjRnTv3h1bt27FunXr8ODBA9SvXx+3b9929qaRC9ixYwemTJmC0qVLO3tTyImuX7+O6tWrI2nSpFi9ejUOHjyIr7/+GunSpeP74qWGDx+OSZMmYcKECTh06JD+ecSIERg/fryzN40Sye3bt/U5pAxWR0c+D+PGjcPkyZOxbds2pEiRQp9vhoaGutx7xKpQdpAZChmhlp2BCA8P123X3333XXz44YeOeo/ITV2+fFnPXEjAUaNGDWdvDjnRrVu3UL58eXzzzTcYPHgwypYtizFjxvA98UJybNi8eTP++OMPZ28KuYgXXngBWbJkwfTp0yN+17JlSz0yPXv2bKduGyU+Hx8fLFu2TGc9WGYrZCajT58++OCDD/TvgoOD9Wfmu+++wyuvvOJSbxNnLOLp/v372LVrl56OingxfX31z1u2bHHU+0NuTP7ji/Tp0zt7U8jJZCarSZMmkfYX5J1WrFiBChUqoHXr1nrgoVy5cpg2bZqzN4ucqFq1ali/fj2OHDmif967dy/+/PNPNGrUiO8L4eTJk7hw4UKk40eaNGn04LYrnm8mcfYGuKsrV67ovEiJGK3Jz//++6/Ttotcg8xeSR69pDyULFnS2ZtDTjR//nydKimpUEQnTpzQaS+SRvvRRx/pz0XPnj2RLFkydOjQgS+Ql85ihYSEoGjRovDz89PnFl9++SXatWvn7E0jF3DhwgX9NbrzTct1roSBBVECjVAfOHBAjzqR9woKCkKvXr30mhsp8EAkgw4yYzFkyBD9YsiMhewrJHeagYV3WrhwIebMmYO5c+eiRIkS2LNnjx6YkvQXfibI3TAVKp4yZsyoRxYuXrwY6ffyc9asWR3x3pCb6tGjB1auXIkNGzYgZ86czt4cciJJl5RiDrK+IkmSJPoia25kEZ58LyOT5F2kqkvx4sUj/a5YsWI4ffq007aJnOt///ufnrWQXHmpEPbaa6/hvffe05UGibI+Oqd0l/NNBhbxJNPWzzzzjM6LtB6Jkp+rVq3qqPeH3IgssJKgQhZd/fbbb7p0IHm3OnXqYP/+/XoE0nKR0WpJcZDvZXCCvIukR0YtQy259Xny5HHaNpFz3blzR6/RtCb7BjmnIMqXL58OIKzPNyV1TqpDueL5JlOh7CA5sjJNKScKlSpV0lVepGRYp06dHPcOkVulP8lU9o8//qh7WVhyH2WRlVT3IO8jn4Ooa2ykTKD0L+DaG+8kI9GyWFdSoV5++WXd+2jq1Kn6Qt5J+hfImorcuXPrVKi///4bo0aNwhtvvOHsTaNErBx47NixSAu2ZfBJir/I50JS46SiYKFChXSgIX1PJFXOUjnKpSiyy/jx41Xu3LlVsmTJVKVKldTWrVv5inop+e8U3WXmzJnO3jRyITVr1lS9evVy9maQE/3000+qZMmSyt/fXxUtWlRNnTqV74cXCwkJ0fsEOZcICAhQ+fPnVx9//LG6d++eszeNEsmGDRuiPX/o0KGDvj48PFx98sknKkuWLHq/UadOHXX48GGXfH/Yx4KIiIiIiOzGNRZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERORQHTt2RIsWLfiqEhF5mSTO3gAiInIfPj4+sV4/cOBAjB07FkqpRNsmIiJyDQwsiIjIZufPn4/4fsGCBfj0009x+PDhiN+lTJlSX4iIyPswFYqIiGyWNWvWiEuaNGn0DIb17ySoiJoKVatWLbz77rvo3bs30qVLhyxZsmDatGm4ffs2OnXqhFSpUqFgwYJYvXp1pL914MABNGrUSD+m3Oe1117DlStX+G4REbkoBhZERJTgZs2ahYwZM2L79u06yOjWrRtat26NatWqYffu3ahfv74OHO7cuaNvf+PGDdSuXRvlypXDzp07sWbNGly8eBEvv/wy3y0iIhfFwIKIiBJcmTJlMGDAABQqVAj9+/dHQECADjS6dOmifycpVVevXsW+ffv07SdMmKCDiiFDhqBo0aL6+xkzZmDDhg04cuQI3zEiIhfENRZERJTgSpcuHfG9n58fMmTIgFKlSkX8TlKdxKVLl/TXvXv36iAiuvUax48fR+HChfmuERG5GAYWRESU4JImTRrpZ1mbYf07S7Wp8PBw/fXWrVto2rQphg8f/sRjZcuWLcG3l4iI4o6BBRERuZzy5ctjyZIlyJs3L5Ik4aGKiMgdcI0FERG5nO7du+PatWto27YtduzYodOffvnlF11FKiwszNmbR0RE0WBgQURELid79uzYvHmzDiKkYpSsx5BytWnTpoWvLw9dRESuyEexPSoREREREdmJwz5ERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERGQ3BhZERERERAR7/R/qAlWSh+64JAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "fig, ax = plt.subplots(figsize=(8, 3))\n", "ax.plot(t, observed, 'k.', ms=2, label='Observed')\n", "ax.plot(t, simulated, 'r-', label='Simulated')\n", "ax.set_xlabel('Time'); ax.set_ylabel('Q [m³/s]')\n", "ax.set_title(f'NSE={nse:.3f}, KGE={kge:.3f}')\n", "ax.legend()\n", "plt.tight_layout(); plt.show()" ] }, { "cell_type": "markdown", "id": "b7e0c82b", "metadata": {}, "source": [ "### Available metrics\n", "\n", "| Function | Metric |\n", "|----------|--------|\n", "| `Nash_Sutcliffe_efficiency` | NSE |\n", "| `Kling_Gupta_efficiency` | KGE |\n", "| `modified_Kling_Gupta_efficiency` | KGE' (CV-based) |\n", "| `Root_Mean_Square_Error` | RMSE |\n", "| `Mean_Absolute_Error` | MAE |\n", "| `Mean_Absolute_Percentage_Error` | MAPE |\n", "| `Pearson_Correlation_Coefficient` | Pearson r |\n", "| `Spearman_Rank_Correlation_Coefficient` | Spearman ρ |\n", "| `Dynamic_Time_Warping_distance` | DTW distance |\n", "| `normalize_series` | Normalize to [0,1] |" ] }, { "cell_type": "markdown", "id": "e0caaa97", "metadata": {}, "source": [ "---\n", "\n", "## Part 2 — friction_law: hydraulic friction coefficients\n", "\n", "Numba-accelerated friction factor computations for pipe/channel flows." ] }, { "cell_type": "code", "execution_count": 4, "id": "8cf79419", "metadata": {}, "outputs": [], "source": [ "from wolfhece.friction_law import f_colebrook, f_barr_bathurst" ] }, { "cell_type": "code", "execution_count": 5, "id": "a77e775d", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Colebrook : f = 0.022175\n", "Barr-Bathurst : f = 0.022184\n" ] } ], "source": [ "# Friction factor for a smooth pipe\n", "k_sur_D = 0.001 # relative roughness k/D\n", "Re = 1e5 # Reynolds number\n", "\n", "f_cole = f_colebrook(k_sur_D, Re)\n", "f_barr = f_barr_bathurst(k_sur_D, Re)\n", "\n", "print(f\"Colebrook : f = {f_cole:.6f}\")\n", "print(f\"Barr-Bathurst : f = {f_barr:.6f}\")" ] }, { "cell_type": "code", "execution_count": 6, "id": "d3fd4cfb", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxYAAAHpCAYAAAAf5apCAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAmv9JREFUeJzt3Qd8k9X6B/Bf0r1Ld0uhLVD23giyZQmI4wqKiriuXlQU0evmjwuvV3FccQ/kiorrIoKgsjfI3rtAC5QO6N5N/p/npCltaSFpkqZJf18/xyRvzrsP8D45S6PX6/UgIiIiIiKygNaSlYmIiIiIiBhYEBERERGRVbDGgoiIiIiILMbAgoiIiIiILMbAgoiIiIiILMbAgoiIiIiILMbAgoiIiIiILMbAgoiIiIiILMbAgoiIiIiILMbAgojIjmJjY3H33XfbdB+rV6+GRqNRr0ayT9l3Q5aTk4OwsDDMnz+/3t7nkydPqnv35ptvoj6R8xk9evQV86Snp8PHxwe//fZbnR0XEdkXAwsiclpz585VD2WS1q9ff9n3er0eTZo0Ud9f7SGJnM+7774LPz8/TJgw4bLvdu3ahTvuuEOVDw8PDwQFBWHo0KH48ssvUVpaapfjdTTBwcG477778MILL9j7UIiojjCwICKn5+npiW+++eay5WvWrEFSUpJ6cGxoPv30Uxw+fBgNVXFxsQos5MHXxcWl0nefffYZunfvjlWrVmHixIn44IMP8OKLL8LLywv33nsv/vWvf9ntuB3Ngw8+iB07dmDlypX2PhQiqgOudbETIiJ7GjVqFH744Qe89957cHW99NeeBBvdunVDWloaGho3Nze77Dc3N1c1j7G3xYsXIzU1Fbfeemul5Zs3b1YPw3369FFNeKRGw+ixxx7Dtm3bsG/fPjgKe1/vNm3aoH379qr2cPDgwXY7DiKqG6yxICKnd9ttt6n23n/++Wf5sqKiIvz444+4/fbba3wge+KJJ8qbwrRq1Uq1c5fmUxWVlJTg5ZdfRvPmzVU+aXv+7LPPorCwsFI+We+VV15BdHQ0vL29MWjQIOzfv79SnhMnTqhmWW+//fZlx7Nx40b13bfffnvFc5UamHHjxqmHSek/8Pjjj192LDX1sZDzu+aaa1QTFvl1XoIuuUZV5efn49FHH0VISIh68B47dizOnDmjju///u//yvPJe1l24MABdZ0bNWqEfv36qe/27NmjjqFZs2aqRikiIgL33HOPuk8VGbdx5MgR1TQpICAAoaGhqnmNXNPExETccMMN8Pf3V9t46623YIqFCxeq85f7VtHMmTPV/qTfRcWgwkhqMir2lTC1nFQnIyNDBSvGdVu0aKFqQ3Q6XbX5pVzExMSoezNgwIDLAhw5Ll9fXxw/flwF03L8UuNii/Jcna+++koF7k8++WSl5ddddx1+/fVXk64JETk4PRGRk/ryyy/lSUb/119/6a+55hr9nXfeWf7dwoUL9VqtVn/mzBl9TEyM/vrrry//TqfT6QcPHqzXaDT6++67T//+++/rx4wZo7b12GOPVdrHpEmT1PJbbrlFP2fOHP1dd92lPo8bN65Svueff14tHzVqlNrePffco4+KitKHhISobRj17dtX361bt8vO5R//+Ifez89Pn5ubW+P55uXl6Vu2bKn39PTUP/XUU/p33nlHbatjx45q36tWrap03HLeFUVHR6v9yPHNnj1b37NnT7Xe4sWLK+W79dZb1XK5nnLO8rlTp05q2YwZM8rzyXtZ1rZtW/0NN9yg/+CDD1R+8eabb+qvvfZa/UsvvaT/5JNP9FOnTtV7eXmpfcr1r7qNzp0762+77Ta1DblXskyOsVWrVvqHHnpILZdrJ8vXrFmjv5oWLVrob7rppkrL5Nq6ubmpe28Kc8qJXOuK91n2JfclODhY/+yzz+o/+ugjVXZkW3ItjBISEtT2OnTooI+NjdX/61//0s+cOVMfFBSkDw0N1ScnJ5fnle17eHjomzdvrt7LNufNm2eT8lz1z8zHH3+stv/cc89ddp2+/vprtY29e/eadF2JyHExsCCiBhFYyMOUPJjLw7f429/+ph80aFC1D0kSdMh6r7zySqXtycOWPDwdO3ZMfd61a5fKJw9rFU2fPl0tX7lypfqckpKid3d3V/uo+NAsD5SSr+IDpzygybKDBw+WLysqKrosAKmOBBKy7vfff1/pAVYeok0JLIzXpuJ+27dvX+lBe/v27dU+kN599901BhYSEFRVdV/i22+/VfnXrl172TYeeOCB8mUlJSUqCJJ78frrr5cvv3jxogpOrnadiouL1bpPPPFEpeW7d+9W+6r4YH8lppaT6gKLl19+We/j46M/cuRIpXWffvppvYuLi/706dOVAgs5r6SkpPJ8W7ZsUcsff/zxy4IC2UZtjtPU8lz1z8y7776rtiPnVJ2NGzeq9RcsWHCFq0lEzoBNoYioQZC29NKER9rWZ2dnq9eamkFJ23rp0CvNfSqSpiTyg8zSpUvL84lp06Zdlk8sWbJEvS5fvlw1vXrkkUdUMxsjaQZT3XFK06CKQ6D+/vvvqh+INAW6EjmeyMhI3HLLLeXLpNnVAw88AFNIExujixcvIjMzE9dee63qfGu0bNky9fqPf/yj0rpybjWRPgtX2ldBQYE6v969e6vPFfdnJJ2sjeTeSJMkuRfSmdooMDBQNfGRJmVXcuHCBbWuNM2qKCsrS71W1wTKknJSHenzI9dWjkHO3Zhk5CkZdWrt2rWV8kvztsaNG5d/7tmzJ3r16lXtUK4PPfSQTctzRW+88QamTp2qmnA9//zz1Z6r8To3xL5MRA0NO28TUYMg7fLloU06bOfl5amHt4oP4BWdOnUKUVFRlz1gSkdU4/fGV61Wq9rGVyRt/eUht2I+ER8ff9kxVX24lfXGjBmjjlPaugsJMuSh8mqdX2U/ciwVgxchD9umkGBL+oHIUKsV29RX3J7xnOPi4iqtW/UaVFQ1r/HhXvozfPfdd0hJSan0nQQ0VTVt2rTSZ+lrIQGY9POourxqP42aVG3zL/00hASepjC1nFTn6NGjqp+JlIHqVL0mVcuOaNmyJb7//vtKy6SPg/TjsWV5rjiqmgQb//znPy/rV1Hdda5aLonI+TCwIKIGQ2oo7r//fiQnJ2PkyJHqYckarP3AdNddd6lftKXDdocOHbBo0SJVQyAPfbaybt061Qm7f//+anhVqfmQkaNk3obqhuo1R8XaiYo1M3J+8kDauXNn1elYOi2PGDGi2s7LVYeErWmZuFonYZmTQu6Z1MpUJA/U8mC+d+9e2Jqco3Rqfuqpp6r9XoKG2pAO15aWE1PLc7t27VQH9P/+97/4+9//Xm0AKYzXuWoQSETOh4EFETUYN954o3oAkiFFFyxYUGM+GXlHmi/JL9cVf+U9dOhQ+ffGV3lAlF+fjb/+ivPnz6sHror5hOSTUZCMZLjTqg+3Qh6u5ZdsqamQ5i5Sw3LnnXde9fxkPzJSkDxYV3w4NGW+ip9++knVAEizq4rzekhgUXUfcs4JCQmVfkU/duwYTCXnvGLFClVjIfNDGMn1qQsSPMioR3IOFUmzMakVkjkXZLQpGUHpSkwtJ9WR/cvM31KLZorqro2MlGXK7OnWLs9GEijIqGEy0teQIUPUJJRSM1KV8TpX3CYROSf2sSCiBkN+Ff/www/VEKbS3KgmMlSnNJV6//33LxvuUx7YpbbDmE+88847lfLNnj1bvV5//fXqVR4e5df///znP5V+Ta+6XsUHXxkiV5q5yPj/UmvRsWPHq56fHM/Zs2crDRErQcknn3xy1XXl1385t4qzSp88eVINy1rR8OHD1avUalQk52YqY01D1ZqFmq6HLcg8FTInRVUzZsxQxyWBnDz4V7V9+3Y1rKo55aQ6UmOzadMmFchVJQ/xMuxrRXIfZEhfo61bt2LLli1X3IetynNF0uxKghbpvyQ1MNU1Q5NrJk3UpIaDiJwbayyIqEGZNGnSVfNI0CHzTDz33HPq4bpTp074448/8Msvv6gO18a5D2S5bE8e3OVhUOYWkAc+efCUzrayDSG1D9OnT8esWbMwevRo9QC3c+dO1Wm2puYh0hxKJvST2Z9NnelZmnnJw6OsKw9z0pxJmqnIL/FXIw+N8gAptSXSZEza+M+ZM0c1D5K+AEYyt8XNN9+sHj7lIVI6XEtbe/n13NRmNNKXQZpcScdfmQFb+o/I9a1ag2BLMveFXBs57orNjmQeDzlvaXrWunVrFWBIzYz82r969WrVLE36oZhTTqojTcBkW1IeZP4Jua4y14Q0w5LAULZXsWzIfZCaAemYLf1f5PrLfCM1NaWyZXmuSo5Ntjdw4EAVeEqNj7G/ipD5Y+QY2MeCqAGw97BURER1MdzslVQdblZkZ2eroTxlrgmZ2yA+Pl7/73//u9JwscahS2Vegbi4OJWvSZMm+meeeUZfUFBQKV9paanKFxkZqYYOHThwoH7fvn2XDUNaUbt27dRcGxWHGb2aU6dO6ceOHav39vZWQ9TK0KnLli0zabjZzz//XJ2nzIXQunVrdf2Mw71WJEPYTpkyRc2l4Ovrq+Y4OHz4sMpXcfhX47qpqamXHaec04033qgPDAzUBwQEqOF/z549W+OQtVW3Iccvw7VWNWDAAHXdrqawsFBdn5qGSJVhdW+//fby+9+oUSP9kCFD9F999ZW6l+aWk+rus6wrZUWGA5bhiOV4ZL4VmeNDhvqtONysbPOtt95S5Uvuj8wBIsPjmnJNbFGeq/szI0PgypDO/fv3Lx9OWIZNluNfvnx5jfeCiJyHRv5n7+CGiIgu16VLF9XRWPoj1HcykpQc79dff10+23N9J6NuSR8S6VNQU0dwsozUiMjQuVKDxhoLIufHPhZERPWQtP+Xh3Vp1lTfSHv6qqRpjoxGJE2cHMXjjz+u+lHIkLdkfdJU7rPPPlNNxxhUEDUMrLEgIqpHZFQn+XX3rbfeUhOKyWRvMlpTfSKjOckxSpt76WgufUUkyUR8H3/8sb0Pj4iI7IQ1FkRE9Yh03J08ebLq1Pztt9/Wu6DC2MFZJriTpkQyK7N0gJaRtqTTMxERNVyssSAiIiIiIouxxoKIiIiIiCzGeSyuQGYglcmmZKZSdjwjIiIiooZGr9eruXyioqLUIB1XwsDiCiSoaNKkibXvDxERERGRQ0lMTER0dPQV8zCwuAKpqTBeyIqziNqyhiQ1NVXN0nu1iNBe27ZkO+aua05+U/JaK4+jsve5sXzX7npYK5+977+t2fv8nLl812Y9lm/nKdu2PgZHK98N8dkkKytL/dBufC6+EgYWV2Bs/iRBRV0FFgUFBWpftviDa41tW7Idc9c1J78pea2Vx1HZ+9xYvmt3PayVz97339bsfX7OXL5rsx7Lt/OUbVsfg6OV74b8bKIpey6+Euf714WIiIiIiOocAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrIYAwsiIiIiIrKYq+WbIFvR6/XQ5+dbZVs6nU5tS5eXB2i1dtmOueuak9+UvNbK46jsfW623L8zl29r5bP3/bc1e5+fM5fv2qzH8u08ZdvWx+Bo5duezyZ6vR71nUbvCEdpJ1lZWQgICEBmZib8/f1tvj8pOCkpKQgLC4NWq1WF7HDXbjbfLxERERHVb/Hb/oKrr2+9fh52vp+tiIiIiIiozrEpVD2m8fJCqx3brVYbkpqaitDQUFUbYo/tmLuuOflNyWutPI7K3udmy/07c/m2Vj57339bs/f5OXP5rs16LN/OU7ZtfQyOVr7t+Wyi8fJCfcfAoh7TaDTQeHtbZ2M6nSqQWm9vy/5SsGQ75q5rTn5T8lorj6Oy97nZcv/OXL6tlc/e99/W7H1+zly+a7Mey7fzlG1bH4OjlW87PptoNBrUd074rwsREREREdU1BhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGQxBhZERERERGSxBhFY3HjjjWjUqBFuueUWex8KEREREZFTahCBxdSpUzFv3jx7HwYRERERkdNqEIHFwIED4efnZ+/DICIiIiJyWnYPLNauXYsxY8YgKioKGo0GCxcuvCzPnDlzEBsbC09PT/Tq1Qtbt261y7ESEREREVH1XGFnubm56NSpE+655x7cdNNNl32/YMECTJs2DR999JEKKt555x0MHz4chw8fRlhYmMrTuXNnlJSUXLbuH3/8oQIWUxUWFqpklJWVpV51Op1Ktib70Ov1NtmXtbZtyXbMXdec/KbktVYeR2Xvc2P5rt31sFY+e99/W7P3+Tlz+a7NeizfzlO2bX0Mjla+G+Kzic6M/do9sBg5cqRKNZk9ezbuv/9+TJ48WX2WAGPJkiX44osv8PTTT6tlu3btssqxzJo1CzNnzrxseWpqKgoKClAXNy4zM1MVHq1WWy+3bcl2zF3XnPym5LVWHkdl73Nj+a7d9bBWPnvff1uz9/k5c/muzXos385Ttm19DI5Wvhvis0l2drbjBBZXUlRUhO3bt+OZZ54pXyYXdOjQodi0aZPV9yf7kdqRijUWTZo0QWhoKPz9/WFrUnCkOZjszxZ/cK2xbUu2Y+665uQ3Ja+18jgqe58by3ftroe18tn7/tuavc/Pmct3bdZj+Xaesm3rY3C08t0Qn008PT2dI7BIS0tDaWkpwsPDKy2Xz4cOHTJ5OxKI7N69WzW7io6Oxg8//IA+ffpcls/Dw0OlquQm1tWNlIJjq/1Za9uWbMfcdc3Jb0pea+VxVPY+N5bv2l0Pa+Wz9/23NXufnzOX79qsx/LtPGXb1sfgaOW7oT2baM3YZ70OLKxl+fLl9j4EIiIiIiKnVq9/tgoJCYGLiwvOnz9fabl8joiIsNtxERERERGRAwUW7u7u6NatG1asWFGpnZl8rq4pExERERER2Yfdm0Ll5OTg2LFj5Z8TEhLUKE9BQUFo2rSp6kw9adIkdO/eHT179lTDzUpfCeMoUUREREREZH92Dyy2bduGQYMGlX82jsokwcTcuXMxfvx4Ndzriy++iOTkZDVnxbJlyy7r0E1ERERERA04sBg4cKAal/dKHn74YZWIiIiIiKh+qtd9LIiIiIiIyDEwsCAiIiIiIosxsCAiIiIiIosxsCAiIiIiIosxsCAiIiIiIosxsCAiIiIiIscfbtYRyGzfkupiPzL0ri32Za1tW7Idc9c1J78pea2Vx1HZ+9xYvmt3PayVz97339bsfX7OXL5rsx7Lt/OUbVsfg6OV74b4bKIzY78MLKoxZ84clUpLS9VnmaCvoKAAdXHjMjMzVeHRarX1ctuWbMfcdc3Jb0pea+VxVPY+N5bv2l0Pa+Wz9/23NXufnzOX79qsx/LtPGXb1sfgaOW7IT6bZGdnm5yXgUU1pkyZolJWVhYCAgIQGhoKf39/2JoUHI1Go/Zniz+41ti2Jdsxd11z8puS11p5HJW9z43lu3bXw1r57H3/bc3e5+fM5bs267F8O0/ZtvUxOFr5bojPJp6enibnZWBhArmJdXUjpeDYan/W2rYl2zF3XXPym5LXWnkclb3PjeW7dtfDWvnsff9tzd7n58zluzbrsXw7T9m29TE4WvluaM8mWjP26Zz/uhARERERUZ1iYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERAwsiIiIiIjI/lhjQUREREREFmNgQUREREREFmNgQUREREREFmNgQUREREREFmNgQUREREREFnO1fBPOT6fTqVQX+9Hr9TbZl7W2bcl2zF3XnPym5LVWHkdl73Nj+a7d9bBWPnvff1uz9/k5c/muzXos385Ttm19DI5Wvhvis4nOjP0ysKjGnDlzVCotLVWfU1NTUVBQgLq4cZmZmarwaLXaerltS7Zj7rrm5Dclr7XyOCp7nxvLd+2uh7Xy2fv+25q9z8+Zy3dt1mP5dp6ybetjcLTy3RCfTbKzs03Oy8CiGlOmTFEpKysLAQEBCA0Nhb+/P2xNCo5Go1H7s8UfXGts25LtmLuuOflNyWutPI7K3ufG8l2762GtfPa+/7Zm7/Nz5vJdm/VYvp2nbNv6GBytfDfEZxNPT0+T8zKwMIHcxLq6kVJwbLU/a23bku2Yu645+U3Ja608jsre58byXbvrYa189r7/tmbv83Pm8l2b9Vi+nads2/oYHK18N7RnE60Z+3TOf12IiIiIiKhOMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLMbAgIiIiIiKLuVq+Ceen0+lUqov96PV6m+zLWtu2ZDvmrmtOflPyWiuPo7L3ubF81+56WCufve+/rdn7/Jy5fNdmPZZv5ynbtj4GRyvfDfHZRGfGfhlYVGPOnDkqlZaWqs+pqakoKChAXdy4zMxMVXi0Wm293LYl2zF3XXPym5LXWnkclb3PjeW7dtfDWvnsff9tzd7n58zluzbrsXw7T9m29TE4WvluiM8m2dnZJudlYFGNKVOmqJSVlYWAgACEhobC398ftiYFR6PRqP3Z4g+uNbZtyXbMXdec/KbktVYeR2Xvc2P5rt31sFY+e99/W7P3+Tlz+a7NeizfzlO2bX0Mjla+G+Kziaenp8l5GViYQG5iXd1IKTi22p+1tm3Jdsxd15z8puS1Vh5HZe9zY/mu3fWwVj57339bs/f5OXP5rs16LN/OU7ZtfQyOVr4b2rOJ1ox9Oue/LkREREREVKcYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcUYWBARERERkcVcLd+E89PpdCrVxX70er1N9mWtbVuyHXPXNSe/KXmtlcdR2fvcWL5rdz2slc/e99/W7H1+zly+a7Mey7f12Lts2/oYHK18N8RnE50Z+2VgUY05c+aoVFpaqj6npqaioKAAdXHjMjMzVeHRarX1ctuWbMfcdc3Jb0pea+VxVPY+N5bv2l0Pa+Wz9/23NXufnzOX79qsx/LtPGXb1sfgaOW7IT6bZGdnm5yXgUU1pkyZolJWVhYCAgIQGhoKf39/2JoUHI1Go/Zniz+41ti2Jdsxd11z8puS11p5HJW9z43lu3bXw1r57H3/bc3e5+fM5bs267F8O0/ZtvUxOFr5bojPJp6enibnZWBhArmJdXUjpeDYan/W2rYl2zF3XXPym5LXWnkclb3PjeW7dtfDWvnsff9tzd7n58zluzbrsXw7T9m29TE4WvluaM8mWjP26Zz/uhARERERUZ1iYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERBZjYEFERERERHUTWHTt2hUXL15U71966SXk5eVZvmciIiIiImpYgcXBgweRm5ur3s+cORM5OTm2Pi4iIiIiInIgrqZk6ty5MyZPnox+/fpBr9fjzTffhK+vb7V5X3zxRWsfIxEREREROUNgMXfuXMyYMQOLFy+GRqPB0qVL4ep6+aryHQMLIiIiIqKGx6TAolWrVvjuu+/Ue61WixUrViAsLMzWx0ZERERERM4UWFSk0+lscyREREREROSwONwsERERERFZjIEFERERERFZjIFFPabT63Ch4IK9D4OIiIiIyLqBRWlpKdauXYuMjAxzVqNa2nJuC4b8MASPrXoMqxNXo1hXzGtJRERERI7fedvFxQXDhg1TE+YFBgba7qhI+Sv5L5ToSrDi9AqVgj2DMbrZaIxrMQ4tGrXgVSIiIiIix20K1b59e5w4ccI2R0OVPNr1Ufw89mdMajsJQZ5BSC9Ix1cHvsKNi27EbYtvw3eHvkNmYSavGhERERE5XmDxyiuvYPr06WqyvHPnziErK6tSIuuKbxSP6T2mY/nfluO9Qe9hSNMhcNW4Yl/6Pry65VUM/n4wnlzzJNafWa9qN4iIiIiIHGIei1GjRqnXsWPHqpm2jfR6vfos/TCcjczdURfzd8g+5DpWty8XuGBA9ACVpEP3bwm/4Zfjv+DIxSNYdnKZStJUakTsCIyKG4V2we0q3Z8rbdtax2jtdc3Jb0pea+VxVPY+N1vu35nLt7Xy2fv+25q9z8+Zy3dt1mP5dp6ybetjcLTy3RCfTXRm7NfswGLVqlVwdnPmzFHJGCSlpqaioKCgTm5cZmamKjwyw/mVDAseptKxrGP4/czvWHVulWoqNf/QfJWivaMxOHIwhkQNQZR3lFnbttYxWrquOflNyWutPI7K3udmy/07c/m2Vj57339bs/f5OXP5rs16LN/OU7ZtfQyOVr4b4rNJdna27QKLAQMGwNlNmTJFJWnaFRAQgNDQUPj7+9t8v1JwpJZB9mdqwQkLC8M1La7B87rnsensJlWTsSpxFZLykjDv+DyVOoZ0xMjYkejm203lt/QPrrnHWNt1zclvSl5r5XFU9j43W+7fWtuuj+XbWvnsff9tzd7n58zluzbrsXw7T9m29TE4WvluiM8mnp6etgsshAw3+/nnn6vRoUS7du1wzz33qIdwZyQ3sa5upBSc2uzPQ+uBgU0HqpRbnIuVp1diyYkl2HRuE/ak7VFJq9GiT2QfjIgbgUFNBiHAI6BOj7E265qT35S81srjqOx9brbcv7W2XR/Lt7Xy2fv+25q9z8+Zy3dt1mP5dp6ybetjcLTy3dCeTbRm7NPswGLbtm0YPnw4vLy80LNnT7Vs9uzZePXVV/HHH3+ga9eu5m6SrMzHzQdjmo9RKS0/DcsSlmHxicXYn74fG85uUMlV64proq7B8NjhKsjwc/fjfSAiIiKiWjM7sHj88cdVx+1PP/0Urq6G1UtKSnDffffhscceUxPoUf0R4hWCO9regdtb344dCTvwV/Zf+OPUHziWcQxrk9aq5KZ1Q9+ovhgWO0wFGb7uvvY+bCIiIiJyMLWqsagYVKiNuLriqaeeQvfu3a19fGRF0T7R6BrXFQ91fgjHM47jj5N/qNGkTmSewOqk1Sq5a93Rt3FfXBdzHfpH9691cykiIiIialjMDiykE/Pp06fRunXrSssTExPh58fmNI6ieWBzFWBIOnbxGH4/9btqMnUy66Tq/C1J5svoHtFdzZ0xuOlghHmH2fuwiYiIiMhZAovx48fj3nvvxZtvvolrrrlGLduwYQOefPJJ3HbbbbY4RrKxFo1aqPSPTv/A0Yyj+P3k76rztzSX2nxus0oyGZ+MLiVNpTr5dEIYGGQQERERkQWBhQQU0jP9rrvuUn0rhJubGx566CG8/vrr5m6O6hG5ry0btVTpkS6P4FTWKRVgrDi9ArtTd5ePLiWa7W2majIGNBmA9sHt4aJ1sffhExEREZEjBRbu7u549913MWvWLBw/flwta968Oby9vW1xfGRHMf4xmNx+skqpeamqedTyU8uxNXmr6pdxYu8JfLr3UwR5BqFf435qVnAZaYqdv4mIiIgaHrMDC5mvQgIL6U/RoUOH8uW5ubl45JFH8MUXX1j7GKkeCPUOxa2tbsUt8bfgxJkTOFR0SHX23nBmAy4UXMCi44tUkn4Z3SK6oX/j/qo2Q4ITIiIiInJ+ZgcWX331lWryVLWjdn5+PubNm8fAogHwdfPFqMajMLr5aBTrirErZRfWJK7BmqQ1qvP3lnNbVPr3tn8j1j8W10Zfq4az7RLaxd6HTkRERET2DiyysrKg1+tVys7OrjS9d2lpKX777TeEhbFDb0Mjc2D0iOih0vQe01W/DJkbQ4KM7cnbVaBx8sBJ/PfAf9VQtu0btcfA2IGq6VSLwBaqXwcRERERNaDAIjAwUD0Eqg6+LVte9r0snzlzprWPjxyMNH26s+2dKuUU5WDj2Y2G2b7PbMD5vPPYkb5DpdnbZyPMKwx9ovqoeTP6RPZBoGegvQ+fiIiIiGwdWKxatUrVVgwePBg//fQTgoKCKnXojomJQVRUVG2Pg5yQdOKW2bwlSdmR+TL+OPIH9mbvxbbz25CSn4Jfjv+ikgYatAlug14RvdAzsie6hnWFp8ulWjEiIiIicpLAYsCAAeo1ISEBTZs2ZRMWMovUaMmkfDfH3oyHwh5Csb4Y289vx8YzhhoNmTPjQPoBlb7c/6XqBN4+pD3a+bXDIN0gdA7vDA8XD151IiIiImfpvL1y5Ur4+vrib3/7W6XlP/zwA/Ly8jBp0iRrHh85KQkSZGhaSdMxHSl5KWoY263ntqqO32dzz2JX6i6V5p+Yr/J3DuusajS6hXdTQYe7i7u9T4OIiIiIahtYyPwVH3/88WXLpeP2Aw88wMCCaiXMOwyjm41WSSRlJ2Hz2c1Yd2od9mTsQVp+WvloU0J1BA9pr4IMSR2CLw19TEREREQOEFicPn0acXFxly2XPhbyHZE1RPtF46b4m9AvoB9CQ0NxKueUqs2QWg1pQiVzZ+xI2aGSTNKn1WjR3K85ekb1VPNodAnrghCvEN4MIiIiovoaWEjNxJ49exAbG1tp+e7duxEcHGzNYyNdKbD/f0DbcYCL2bfKqfpnNAtoptKE1hNUR3AZ1laCCgkydpzfgaScJBzNOqrS/EPz1XrRvtHoGNpRpU6hnRAfEG/vUyEiIiJyWmY/rd5222149NFH1QR5/fv3V8vWrFmDqVOnYsKECbY4xoZrzwJg4UPA2n8D170ExA+Tp2w0dBJoxAbEqiS1GuJczjmsProaxwuOY3vKdhzPOK6CDUm/JfxW3nyqhX8LdIvshk5hnVSwEe4dzoEIiIiIiOwRWLz88ss4efIkhgwZAldXw+o6nQ533XUXXnvtNWscE5XTAF6NgNRDwDe3ArHXAsNeAaI68xpVIQHCoMhBGB82HlqtFllFWdiXtg97UveotDdtLzIKM3Ag44BK/z34X7WezKVhrNWQPhutAlvx2hIRERHVRWAhc1YsWLBABRjS/MnLywsdOnRQfSzIuko6jMcPWe3wt/wf4Lr1Y+DkOuCTAUDH8cDgF4DAJrzkNfB39y8fdUqo5lOZp7D+xHqcLDyJPWl7cOTiETWXxvLTy1UyauzdGO1D26NtSFu0DW6LNkFtEOARwGtNREREdAW1brgvs29XNwM3Wc83W0/jxd8SMafRILx83Q0YeOZjaPZ+b2gitX8h0Psh4NppgCcfek1pPtXUvymGRg1V/YSkViO/JB8H0w+qGo3dqbuxP30/zuWew5m8Mzhz6gx+P/V7+fqNfRurAEMm8TMGG8Fe7FNEREREZFFgkZSUhEWLFqlRoIqKiip9N3v27NpskqoR7u+JqABPJF3Mx+SF+ejT7F68dtNdiNvxL0PtxYZ3gJ3/BQb8E+h+D+DixutoBi9XL3QN76qSUXpeOjad2IRzunM4dOEQDl44iMTsRJzJOaNSxZoNGSK3dVBrxAfGI76RIcX5x8GN94GIiIgaILMDixUrVmDs2LFo1qwZDh06hPbt26s+F9LUpGvXSw9oZLnh7SLQPz4UH645jo/XHMemE+kYkgDc3vM1PH3jKfiuewlIOwIsfQrY8jFw3Uyg9Wh28LZAI89G6B7SvbxWQ2QWZuLwhcMqyJBaDanlkFGpZFI/SWuT1pavLzOGS6dyCTJaNmqJFoEt1Psonyh2EiciIiKnZnZg8cwzz2D69OmYOXOmGhnqp59+Ug9hEydOxIgRI2xzlA2Yl7sLpl3XErd2j8as3w5hyd5z+HpLIn7d44Fpg7/BHR5r4LLmdeDCcWDBHUCT3sDwV4Ho7vY+dKch/St6RvZUySi3OFfVaEjAcSzjGI5ePIqjGUfVcvksaWnC0vL8Pm4+5UGGvErNRlxAHMJ9wtUcHEREREQNLrA4ePAgvv32W8PKrq7Iz8+Hr68vXnrpJdxwww146KGHbHGcDV50I2/MmdgVd55Ix8xfD+DguSzMWHIE88PjMfP65eiT/DWw8X0gcTPw2RCg3Y3AkBlA0OWTGZLlJFAwzvptJLV20kfDGGRI53B5fzLzpAo4pB+HpKrNsWL9DUPnSqChkn8cYvxj4OnqyVtFREREzhtY+Pj4lPeriIyMxPHjx9GuXTv1OS0tzfpHSJX0bhaMxY/0w7dbT+OtPw7jyPkc3DZvP65rOwIz7pyA6J1vA7vmGybWO7gY6PkA0H864BnIK1kHHcSjfKNUGtBkQPny4tJinMw6WR5wnMg4gYSsBCRmJRo6kF84qFKlbcGwLRVw+Mep4KOJXxOVIn0j4aptuBMmEhERUf1k9tNJ7969sX79erRp0wajRo3CE088gb179+Lnn39W35HtuWg1uKN3DMZ0jMI7K45g3qZT+PPAeaw5nIr7rn0ED9/7ALxXzwSOrwA2zwF2fg30fRSIu5m3xw6kM7exc3dFxbpiJGUnqRoNCTQSMg3pROYJZBdll3cY33BmQ6X1pB+HBBfGQKNiivaLVrUgRERERPU+sJBRn3JyctR76Wch72Vei/j4eI4IVccCvN0wY0w73N6zKV5afADrjqbhg9XH8eN2Dzw98j2M630Y2uUvAuf3QbvyZYR6fwQMegboehdHkKoH3LRu5c2fBmFQpSZVFwouqCBDajok0JDaDRmdSlKRrqj8fXVk0j8JMCRF+kSW16JEeEVAq2N/DiIiIrJjYPHee+/hgQcegKenp+pXIRPiGZtFffTRRzY6NDJVfLgf5t3TE8sPpuCVJQdwKj0P077fjXlNAvF/Y35F54t/Qr/qFbhknAaWTAM2zQGGvAC0HccRpOppkyqZI0NS94jKnfB1ep0aiUqCCqntOJ19ujzIkOAjuzhbTfonaUfKjmq3H+oVqmo8ZKQq46sEHhKERPhEwNfNlyNYERERkW0Ci2nTpmHChAkqsIiLi8O5c+fUSFBUvx5Gr2sbjv4tQ/DF+pN4f+VR7ErMwLgPNuGmrq0wfeJa+O76FH47P4ZGRpD64W4gqgsw9P+AZgPtffhkIhlBSh7+JfWI6FHpO6npkKFxJciQgONszlmczT2Lcznnyl8LSguQmp+qkkwMWB1pShXuHa6SzNUhI1fJq6QI7wj1GuQZBBetC+8bERERmRdYREVFqWFlpU+FPLzIBHkFBQXV5m3atKkpmyQb8XB1wUMDm+Omro3xr2WH8POOMyot3ZuMid2G4rG/3wffHR8Dm94Hzu4E5t0ANBtkCDCiOvO+OHhwGegZqFKHUEOtYkWlpaU4euYoir2KkZyXrAIPGcWq4mtWUZbqUC5NsCTVxEXjglDvUBVkSA1IsKehhiXEK6T8vUqewfB287bxmRMREVF9YFJg8fzzz+ORRx7Bww8/rB5eevSo/EupkIBDvpOHF2ej0+lUqov9yHW0xr5Cfd3x5i0dcWevpnjlt0PYfuoiPtt8Dr8euIAnh92BcVMmw2XDW8C2L6E5sQr4ZBX07W6EftDzQFAzmxyjueuak9+UvNbK46jkvALcAhDaKBTtgg0juVWVV5ynajPO550vnwBQUvnn/BSk5aehVF+K5Nxkla7G29W7PMjw1fgiMiASId6GAEQCkUYejVQwJK9+7n61mtfDWvetPpZva+Vz5rJdH87Plvu3d/muzXos385Ttm19DI5Wvhvis4nOjP2aFFhI/4rbbrsNp06dQseOHbF8+XIEBwfDWc2ZM0clY5CUmppaYw2NtW9cZmamKjzGWZ8tFekBvD8uDiuOBOA/65JwPqsQ03/cg8/CvDG1/wPoPuFW+P71LjyPLoZGDVH7K/Jb34KcblOg8wmz6jGau645+U3Ja608jsrUc/OCF2K1sYj1jQV8L/++VFeKi0UXkVaQhrTCNFwsvKg+Xyi8gIyiDPUqn2V5oa4QeSV5yMvOu9TZPKXmY9RCCz83P/i7+6sgKMA9QH2WV/ksy/3d/NVn46uPqw+gh1XuW30s39bK58xluz6cny33b61t13Y7tVmP5dt5yratj8HRyndDfDbJzs62/qhQMst2+/bt8eWXX6Jv377w8PCAs5oyZYpKWVlZCAgIQGhoKPz9/W2+Xyk4Uusj+7N2wZkQGop+zQKx5GguPlxzAodS8vDQj0cwvF04nh7xIWIGPQmsfAmaY8vhfeA7eB35BehxH/R9HwO8g6xyjOaua05+U/JaK4+jsua5RSLyqnnkL0AJKqSGI70gHal5qTidehoFrgW4WFAWmOSnIaMwQyWZRFAHHTKLM1VKRPWjXlUlc374uvvCx8VH1Xz4e/irmg8JQkx59XDxqNfl21r5nLls14fzs+X+rbXt2m6nNuuxfDtP2bb1MTha+W6Izyaenp62G2520qRJaGjkJtbVjZSCY6v9ebpJ/4sWGN8zBrP/PILvtp7G7/vPY9WhVNzdNxZTbv4OAee3AMtnQpO0Fdj0H2i2zwX6TDEkT3+Lj9Hcdc3Jb0pea+VxVHV9bn4ufvDz8EMc4tRfjCneKWrgh+r2LxMJSoAhNR0ZBYZgQ30uuHjpfZXvJBjRQ6/m/ZD/kvOv3jSrKnetu5pJ3Zjc4Y5A70A1Opb0D5HXit/7uPuoWhIVzJQtkzyeLp42K9/WyufMZbs+nJ8t92+tbdd2O7VZj+Xbecq2rY/B0cp3Q3s20ZqxT07f2wCF+HrgtRs74K4+MXh1yUE1/8Una0/gx+1JeHxoPG67exlcTywHVr4MJO8F1rwObP0YkNqLHvfZ+/DJiScSlA7hkkxVVFqkOpxnFmTi1PlTcPVxVUPuSqAhyyu9Fmap98aUU5SjghKZF6SosEgFLeUu1vIctG6qT4mXm5caXUuCDXlVn10My9RyV0/1XWlBKUIuhKjgxLi8Yh71qvVUNT/SBM1ZAwIiInIODCwasNYR/mr+i9WHU9X8F8dTc/HCL/vx1aZTeO76Lhj0wFrg4C/AqteAtCPA8hnQbP4A3p0fAPpPAdw5wzPZl7uLu+oAHuQRBJ9CnxprQ6ojc4LkFOeoAENqPiRJAHI27SxcvV3Vw7xxuSTJK53bK74avyssLSyfTT2zKFMlW3DVuqqmW5Lk3I2vEqRIzQtKAT8vP3i4Vs6j3mvdUZxfjKD0IBW0VPxOAiIJ7NRrxeRS/Wc5jtp0siciIufGwKKBk6q1Qa3D0C8+BN9uPY23/zyCYyk5mPzlX+jfMhTPXz8ULR8aA+z9Hlg9C5qM0/Bf/zL0e78EBvwT6HQ74MJiRI5HHoxVZ3D3S/2nVHMt15qba9VEAoqcwhycSj4FnwAf1Wldhu01poKSgsveS3ByMeci4AbD8tIC5BfnX3pfYX2jEl2JShLM2JurxlUFGhJkmBKMVF0m68mwxZLK32sN72Xb8r7id+pVvq/wnfGzCnSgRXZmNoJ1wWr7xnUqbbvKdqV/jrxXzQvkP41WfW98L0n+jiQiItPwiZAUNxct7uoTixs6Ncb7q45i7saTWHskFSOOpuK2nk3x2NCbEdr+Fuh2fAX96jfgkpkELHoEWP8OMOhZoN1N0giPV5MaJHlYDvAIQIRXBMICTQtKVBCTcvUgRkanS0xOhF8jP5ToS1TtiCRpBlbxvQQsqRdT4eHjYWjeVfZ9eb6SQmTkZkDrplXfl69XWqD6t0hwJEGLvBo/G5Pkk2ZjFcmxlJSUwNmVBx9lgUbFJMsl8DB+r9fp4eZqCGqulteU5cb3cumluZ6Xpxc0Wo0KfGQd9Z+0u5YAqOy9/GesTTIuN36W9/n5+fBN9FXbMX5Xdd3LtqUHcnNz4ZfqV368Ffdt3HbFdWWdnJwc+Gf6q2CtuuOU/2SQh+ysbPjn+as/B7LMuD1jUCfvpYYxKzMLgQWB5cdenldj2E5mRiYCiwLL/zxV/L7idsveXraPivkr3v/LtlHhfcX81W3jsv2ZsI2qxybndiHnAnI9citdo0rltJoAuLp8NS6vdpGm0t9X6XnpKM4pvuzvq2qPx8RjFPLnRgbz0Ofqy++tKetWzSfHKKMSavNN7PtQtr6sJ01hXfJd4OLiYtJ+pDmte6F7pZrb6o5Rr9Or5rmehZ7q3KrLJ9sz1phXPO6q17/ij0xOFVjIP3Jz587FihUr1D+KVce2XblypTWPj+pYgLcbnru+LSb2isHrSw9h2f5kzN9yGgt3nsHfBzTHPX0nISfqOoSdWgTthncAmcX7p3uBdbMNAUbr6+VPDe8bkZXIP0LS1KmRZ6OrjgqV4lVzoGJqIFMT6eNRMdioGnxUXaaClKp5yj5LYCPbk+BEvepK1NwoxldZZvxszFP1exXY6Cp/V1hUCLgYmrnJfqpuq3zdsu2aQgIqlbdyXEVEVKfCPMPw59/+dL7AYurUqSqwuP7669Xws6wmdk6xIT746M5u2HIiHa/9dhC7kzLVSFJfbz6Fe3tF4J6BU6DtPhnY/CGw8T9Ayn5gwUQgoiMw8Bmg1UgGGERORDUh0rrAE6YPO1iXahM4yS/BEmwYXyUgkVTxvfFz1TxV80rAkpaehoBGAerXX/ms8sh/urJtQG9Yjitv97LlOp0aR97b17v813nZhryq/+SzXlf+Xv6r+tm4vbzcPHh6G+6h5DG+Vs1bcfuy/7z8PDXkZFmO8n1IwFXtsej1yC/Mh7u7u7oelfZRtn21rgSBxcVwdbv0OGLMV/VzUXER3NzcypcZgz3jsZQUl5Rvp7ptGO95+fuqecreG8+v4rKq+au+r7iNituvusxw2OZtQ33W6WWin2oD3Kq1iTUtq3o+V8tbdT3jRMhX3beJ+6i0TN5qLDu/6vZbU16qR4HFd999h++//x6jRo2yzRFRvdKrWTD+94++WLz3HP79+yEkXsjHrOWn8NOedDwzqg0G9n8SGhkpatP7wJaPgeQ9wHe3AZGdDQFGy+EMMIioXpKHJOl3Idyks4ulgU1pCsJCa1cjdNVtW1DbZOl2arOeqetYK5+1rlF9VB/OzZbHUF/Kd3VzRFQXlEjLnZTUyn/Wawp+dDodUlNSERpWtu1qYhz5seFq2zNuxxGYfQfl14cWLVrY5mioXtJqNRjbKQrLpw3Ac6Naw9/DBUekg/fcvzDxsy3Yd9EFGPIiMHUP0O9xwM0HOLcL+HY88Olg4Mgf8lOCvU+DiIiIqFqq30uVVG3fqgqDRxiTcaAI1wrJOFhFpQEuXC5PMkKfjNqnXsuSccS+ikmWO2Vg8cQTT+Ddd9+tscqJnJeHqwvu7ReHHye3x/3XxsHdRYuNx9Mx+j/r8fiCXUgq8gKG/h/w2B6g71TAzRs4uwP45m/AZ0OAo38ywCAiIiJyUmY3hVq/fj1WrVqFpUuXol27duXtHY1+/vlnax4f1UP+nq54ZmRrNYrUm38cxi+7zuJ/O89gyd5zmHxNLP4xqAUCrnsJ6PMIsPFdYOtnwJnt0H57K4LCOwNDngfih7KJFBEREZETMbvGIjAwEDfeeCMGDBiAkJAQBAQEVErUcDQJ8sa7E7pg0cN90btZEIpKdPh47QkM+PcqfL4+AUWewcCwVww1GH0eht7VC+7nd0H7zS3AF8OB4ytZg0FERETUUGssvvzyS9scCTmsjtGB+Pb+3lh1OAWzfjuEoyk5eHnxAXy18SSeGtEK13eIhGb4q9D3eRh5f86C98EF0CRuAf57I9CkF9D/SaAFazCIiIiIHFmtu9+npqaqZlGS5D01bNLJaXDrcCydei1m3dQBoX4eOH0hDw9/sxM3frARm0+kA77hyO77LPSP7AR6PQS4eAASYMy/Bfh0EHBoCWswiIiIiBpKYCGzb95zzz2IjIxE//79VYqKisK9996LvLw82xwlOQxXF62aqXv19IF4bGg8vN1dsCsxAxM+2YzJc7fhaGoe4BcBjHy9vImUoZP3TuC724GP+gH7/ydjq9n7VIiIiIjIloHFtGnTsGbNGvz666/IyMhQ6ZdfflHLZMQoIuHj4YrHhrbE6icH4s7eMXDVarDmSCrumn8Q077fjcQLZQHG8FeBx/YC/aYB7n7A+X3AD3cDH/QG9iwAdKbNjktEREREDhZY/PTTT/j8888xcuRI+Pv7qyST5X366af48ccfbXOU5LDC/Dzx8rj2ag6M0R0i1ZQvC3edxeC3VuP/Fu1Hek4h4BMCDJ1hqMEY8DTgGQCkHYZ24YMI+W4UsPNroLTY3qdCRERERNYMLKS5U3h4+GXLZaZDNoWimsSG+OC92zpj7m2t0a9FMIpL9Zi78ST6v7EK7y4/itzCEsA7CBj0jKEGY/AL0HsFwTXrFLS/PgK81xX463OgpJAXmYiIiMgZAos+ffpgxowZKCgoKF+Wn5+PmTNnqu+IrqR1uA/m3dMTX9/bCx0aByC3qBRvLz+ihqiVUaRkyFpVY9F/OvRTdyOrzz+h9wkDMk8DS6YB73YCNs0BCnN4oYmIiIgcebhZmXV7+PDhiI6ORqdOndSy3bt3w9PTE7///rstjpGcUL/4EFzTvC9+23cOb/5+GCfT8zBj0X41/8UTw1piTMcowN0XeZ3uge/AqdDs+hpY/w6QfRb4/Vlg7b+Bnn8Hev3dUNNBRERERI4VWLRv3x5Hjx7F/PnzcejQIbXstttuw8SJE+Hl5WWLYyQnpdVqMLpjFIa3i8CCvxLx7oqjaojaqd/twidrT+DJYS3RKkAPuHkZAohudwO7vwU2vAtcOAGseR3Y+B/D8t4PAag8CzwRERER1ePAQnh7e+P++++3/tFQg+TmosUdvWNwU9fG+GJ9Aj5ecwL7z2bh7rnb0C3aD8+P8UCXmCDA1cMQRHS5EzjwC7B+NpC8F9g8B5qtn8C/5Vhg8D+BsFb2PiUiIiKiBsekwGLRokVqFCg3Nzf1/krGjh1rrWOjBsbb3RUPD47H7b1i8MGqY5i36SS2J2Xjxg83YXi7cDwxrBVahvsBWheg/U1AuxuBYyuA9W9Dc2o9vA/9BP2hn4E2Y4BrpwFRXex9SkREREQNhkmBxbhx45CcnKxGfpL3V5p9ubS01JrHRw1QkI87nh/dFpOuicHri/di6cF0/L7/PP44cB43dm6s5sdoGuwtBQ6IH6qS7tRmFK18HZ6nVgEHFxlSs0GGACP2WkNeIiIiIrLvqFA6nU4FFcb3NSUGFWRNjQO98MKwWCx9tB9Gto+AXg/8vPOMmgPjuf/tRXLmpZHJ0KQnMkZ+BN2DG4CO4wGNC3BiFfDVGOCzIYamUzoGvURERET1ZrjZefPmobDw8rkEioqK1HdE1hYf7ocP7+iGXx/uhwEtQ1Gi02P+ltNqiNpXlxzAhdyiS5nD2gI3fQI8uhPocT/g6gmc2Q58fxc0c3rAa998oDiPN4mIiIjI3oHF5MmTkZmZedny7Oxs9R2RrXSIDsBX9/TEggd6o0dsIxSW6PDpugQ1yd7by48ip7BCjUSjGOD6Nw2T7fV/EvBqBM3FBASsfwmadzoAq14DctN4s4iIiIjsFVjo9XrVl6KqpKQkBAQEWOu4iGrUq1kwvv97H8yd3APtovyRU1iC/6w8hpu+3KuGqc0vqhBg+IYBg58HHt8P3Yh/ocQvGpr8C8CafwFvtwN+fQxIP86rTURERFRXw8126dJFBRSShgwZAlfXS6tK34qEhASMGDHC0uMhMomUw4GtwtA/PhTL9ifjrT8O43hqLl5fdhhfbDiJRwa3wPgeTeHuWhY7u/sAPR9AWtPRCLuwFVqZ/+LsDmD7l8D2uUDr64FrHgWie/AOEBEREdkysDCOBrVr1y4187avr2/5d+7u7oiNjcXNN99cm2MgsmiSvVEdIjG0dSj+u+4wvth6Hmcy8vHCL/vx8doTeHxoS4zr0hgu2rJaNq0r0HacYajaUxuBje8BR5YBhxarpInuAY92k4CQCbJx3hkiIiIiawcWM2bMUK8SQEyYMAEeHh6mrkpkc64uWlzfNhi392uJH7afUU2jki7m44kfduPDNcdVgDG8rWFkM0Wa88X2NaTUw4YZvPcsgCbpLzRK+gv6v942zObd6TbA41IQTURERETVM/sn2bZt26pai6q2bNmCbdu2mbs5IqvycHXBXX1isfbJQXh6ZGsEeLnhWEoOpnyzA9e/vwGrj11U/YQqCW0F3PC+6uit7zcNOo8AaC4cB36bDrzdFvjjBSAjkXeKiIiIyJqBxZQpU5CYePlD1pkzZ9R3RPWBl7sLHhzQHOv+OQhTh8TDz8MVh5Oz8fTiExg7ZyNWHDx/eYDhFwH94BeQescq6Ea+AQQ1BwoyDc2l3u0EfD8JOL1FRjCw12kREREROU9gceDAAXTt2rXazt3yHVF94u/phseva6kCjCkDm8PbTYv9Z7Nw71fbMO6DjVhzJPWyAEPv5mOYA+PhbcBtC4C4AYC+FDiwEPhiGPDpYGDPD0Bpsd3Oi4iIiMjhAwvpW3H+/PnLlp87d67SSFFE9UmgtzueGNYSP9/TAQ/0j4OXmwt2J2Zg0hdb8bePNmHjsbTLazCk83arEcCkRcBDG4EudwAuHobRpH6+D5D5MNa9BeRdsNdpERERETluYDFs2DA888wzlSbJy8jIwLPPPovrrrvO2sdHZFWBXq54ekRrrH1qEO7tFwcPVy22nbqI2z/bggmfbMbWhBqChPB2wA1z1HwYGPQc4BsOZJ8DVrwEzG4L/DoVSDnEu0VEREQNltmBxZtvvqn6WMTExGDQoEEqxcXFITk5GW+99ZZtjpLIykL9PPDC6LYqwJjUJwbuLlpsSbiACZ9uwSM/H8GO0xerX9E3FBjwlGFG7xs/BiI6AiX5hrkwPugFfDUWOLgYKC3hPSMiIqIGxezAonHjxtizZw/eeOMNNUJUt27d8O6772Lv3r1o0qSJbY6SyEbC/T0x84b2WP3kQEzs1RRuLhr8dTobt3y0WTWT2pWYUf2Krh5ApwnA39cCk5cCrUcDGi2QsAZYMBF4rzOwbjaQm8Z7R0RERA1CrTpF+Pj44IEHHrD+0RDZSVSgF169sQP+3j8O//5tH347eEF17JY0tE0YHhvaEu0bB1y+osyHEXONIWWcBrZ9AWz/CshMBFbMBFa/DrS/Ceh5P9C4mz1OjYiIiKj+BBaLFi3CyJEj4ebmpt5fydixY611bER1LrqRN567LhbTRrTD+6tO4H87k7D8YIpKQ1qHYerQeHSMDqx+5cCmwND/AwY8Dez/Gdj6CXB2J7D7W0OSwKLnA4ZZv6XGg4iIiKihBRbjxo1TfSjCwsLU+5poNBqUlpZa8/iI7CIm2Adv3doJ/xjUHP9ZcRSLdp/FikMpKg1qFYqpQ1uic5MaAgw3T6Dz7YaUtN0QYEigcWY78L+/A78/B3SbBHS/BwiIrutTIyIiIrJfYKHT6ap931DIOdfFecs+ZMhTW+zLWtu2ZDvmrmtOflPy1iZPXLA3Zt/aCQ8Pao45q4/jl11nsepwqkoDWobg0cEt0KVpo5oPLKoLMO5D4LqXgB3zoNn+BTRZZ9Uwtfr1bwOtRkHfbTLQbKChj4aDli9779+Zy7e18tn7/tuavc/Pmct3bdZj+Xaesm3rY3C08l0fnk3qmjn7NauPRXFxMUaMGIGPPvoI8fHxcFZz5sxRyVj7kpqaioKCgjq5cTKMrxQercyhUA+3bcl2zF3XnPym5LUkjy+Afw6IxG0dG2Hu1mT8figda46kqdSrqT/u7R2JjlGS6wpa3QnE3waPkyvhve9reJzdAhxaDM2hxSjxb4L8Nrciv/XN0HkFw9HKl73378zl21r57H3/bc3e5+fM5bs267F8O0/ZtvUxOFr5rk/PJnUlOzvbNoGF9LGQEaGc3ZQpU1TKyspCQEAAQkND4e/vb/P9SsGR5mSyP1v8wbXGti3ZjrnrmpPflLzWyBMWBnRv1RQn03Pxwerj+N/Os9hyOkulvi2CVQ1Gj9igK59YxB1A7zugSzkIzfYvgT0L4JqVCL8tb8H3r/eANqOh7zoZiO1n6BzuAOXL3vt35vJtrXz2vv+2Zu/zc+byXZv1WL6dp2zb+hgcrXzXx2cTW/P09LTdqFB33HEHPv/8c7z++utoKOQm1tWNlIJjq/1Za9uWbMfcdc3Jb0pea+VpFuqHN//WGY8Obok5q47hpx1J2HAsXaVrmgdj6pB49Gp2lZqHiHbA9W8C180E9v0MbP8SGumHsf9/0Oz/HxDcApBmUtJXw/sqwUo9KF/23r8zl29r5bP3/bc1e5+fM5fv2qzH8u08ZdvWx+Bo5bu+PpvYijn7NDuwKCkpwRdffIHly5erOSxk6NmKZs+ebe4miRxW02Bv/OuWjnh4cAt8sPoYftiWhI3H01Xq3SwIU4e0RJ/mVwkw3H2Arnca0rndwLYvgb0/AOnHgD+eM8zu3fYGQ2fvpr2tWotBREREZC1mBxb79u1D165d1fsjR45cFk0RNURNgrwx66aOmDJIAozj+GFbIjafuIDNJzajZ1wQHhsSrwKMq/4ZiewEjHkHGPYysPdHw7wYyXuAvd8bUmhrQy1Gx1utVotBREREVKeBxYkTJxAXF4dVq1ZZZcdEzjoPxms3dlABxoerj+H7v5KwNeECbv9sC7o0DcTDg1pgcOuwqwcYHn5A98lAt7uBszsMtRj7fgJSDwHL/gn8+YJhtm+p5YgbKPWUdXWKRERERNUy+WlERoGS0ZGMxo8fj/Pnz5u6OlGD0jjQC6+M64A1Tw3EpD4x8HDVYufpDNz71TZc/956LNlzDqU6/dU3JAGITKx3w/vAE4eAUW8CER2A0iLD3Bj/vRF4tyOwapZh5m8iIiKi+h5YyBBXFf3222/Izc21xTEROY3IAC/MvKE91v9zMP4+oBl83F1w4FwWpnyzA8PeXoOftiehuNTE8aE9A4Ce9wMPrgceWAP0uN+wLDMRWPM68E5HYN44QxOqYtsPj0xERERUEdtPENWBUD8PPDOyDTY8PViNGOXv6Yrjqbl44ofdGPzWaszfcgqFJWbMWh/V2TCi1BOHgZs+A+IGSPgPnFgF/HQv8FYr4LengOS9tjwtIiIiIvMDC2kTXrVdODtrE5kn0Nsdj1/XUgUY/xzRGiG+7ki8kI/n/rcP/d9Yhc/XJyCvqMT0Dbp5AR3/BkxaBEzdDfR/CvBvDBRkAFs/Bj7qB3zcH9j6KZB3gbeLiIiI7N95W5pC3X333fDw8FCfZSbqBx988LLhZn/++WfrHyWRk/HzdMNDA5vj7mti8d1fp/HJ2hM4l1mAlxcfUPNi3NsvDnf2iYG/p5vpG20UCwx+Dhj4NHB8FbDzv8ChJYYhbM/thmbZMwiMGQj0vBuIHwa4utvyFImIiKiBMTmwmDRp0mUT5RGRZbzcXTC5bxwm9orBzzuS8OGa4ziVnod//34YH605jsnXxKrvG/mYEQRoXYD4oYaUm65m9sbub6FJ3gPPhD8BSV5BQIdbgE4TgKiunBuDiIiI6i6w+PLLLy3fGxFVy91Viwk9m+KWbtFYvOecqrU4mpKD91Yew2frEzCxV1Pcf20zhPl7mncFfYKBPv9QSXduL/I3fwHv40ugyTkPbP3EkEJaGgKMjuOBgGjeISIiIqoVdt4mqkdcXbQY16Uxfn+sPz66oyvaN/ZHXlEpPl2XgH5vrMILC/ch6WJe7TYe3g7Zff4J/WP7gDt+Ajr8DXD1AtKOGGb3frs98NUYYNc3QGGOtU+NiIiInJzZM28Tke1ptRqMaB+J4e0isPpIKuasPIZtpy7iv5tP4dutpzG2UxT+PqA5WkX41WLjrkCLoYZUkAUcXATs/g44uQ5IWGtIS54wTMAnwUfzQYCLGX09iIiIqEFiYEFUj8nIa4NahWFgy1BsSbigmkitO5qGn3eeUWlI6zDVCbx7bFDtduDpD3S5w5Bkgj3pj7HrW+DCcWDv94Yk/THa3mAIMpr24SzfREREVC0GFkQOEmD0bhas0t6kTNWx+7d957DiUIpK3WMaqQBDghCp7aiVwKZA/yeBa6cDSduAvT8A+/8H5KYA2780JBnKtt2NhiAjshM7fRMREVE5BhZEDqZDdADmTOyKhLRcfLL2OH7afkY1k7r3q21oFe6nZvge0ykKbi617EIl89U06WFIw18zNJGS2bwP/gpknQE2vW9IwS2A9rcYRpcKibf2aRIREZGDYedtIgcVF+KDWTd1xPp/DlLBhK+HKw6fz8a073dj4L9X48sNZk62Vx0XV0Mfi3FzgCePAuPnG2osXD2B9GPAmteB97sbJuHb8B6QmWSt0yMiIiIHwxoLIgcnQ9A+M7IN/jGwBeZvOYUv1p/EmYx8zPz1AN5bcRSTronFpD6xCPCy8I+7qwfQZrQhFWYDh34D9v0IHFtRPgkf/nwBiO5h6JMhSZpXERERUYPAwILISQR4uang4p6+cfhpRxI+XnMCpy/k4Z3lR9X78T2iMa61P8LCrLAzDz+g03hDkkn4DiwE9v0EnNoIJP1lSH88b5h8rzzIiLHCjomIiKi+YmBB5GQ83VzUTN7juzfB0n3JqqP3/rNZmLvxFP67GRjX+QIeHNAc8eG1GKq2pkn4etxrSNnJhr4YB34BTm0Azu4wpOUzoInsBJ+mQ4EetwMhLayzbyIiIqo3GFgQOfFke9KJe3THSDVE7Yerj2PTiXT8tOOMSkPbGIaq7RZTy6Fqq+MXAfS835ByUi4FGSfXQXNuN/ykudSWt4DwDkA7qckYx47fREREToKBBVEDGKq2f8tQ9GsRjFV7EvD93ov448B5LD+YopIMVXvftc1wXdtwuNR2qNrq+IZdqsnITYPu4K8o3vUj3M9sgub8XkDSyleA0DZA61FA6+uByC6cJ4OIiMhBMbAgakDaRfjgw45xSEjPwydrTuDnnUlqqNptp7YjNtgb917bDLd0jYaXu4t1d+wTAnSdhIvRIxHm6wrNkaWGfhknVgOpBw1p3VuAXyTQaqQhyIjtD7i6W/c4iIiIyGYYWBA1QM1DffGvWzpi2rCW+GrjSXy9+RROpufhhYX7MPuPw7izTyzu6hODEF8P6+/cOwjoeqch5V8Ejv4JHFoCHFsOZJ8Dtn1hSO5+QPx1hiBDXj0DrH8sREREZDUMLIgasHB/Tzw1ojWmDGqB77cl4vP1CUi6mK+GqZVO3zd3jcZ918apQMQmvBoBHW81pOICw2R8hxYDh5cCOeeB/T8bktYViL3WEGRIjUZAtG2Oh4iIiGqNgQURwcfDFZP7xuHO3jH4ff95NaP37qRMfLv1tErS0fv+a5uhZ1yQ6rNhE26ehpoJSde/bRhNSoIMmS8j7TBwYpUh/TYdiOgIxA8DWg4HGneTniS8i0RERHbGwIKILv2F4KLF9R0jMapDBP46eRGfrjuB5QcvdfTuFB2gOnqPbB+h8tqMVgtEdzekof8HpB0DDi8xNJlK3Aok7zGkdW8CXkHQtBgCz7BegO+NgG8I7ygREZEdMLAgostIrYTUTkg6npqjmkj9tD1J1WI88u1ONA70wr394jC+RxNV22FzMu9FyFSg71QgJ9XQH+Po78CxlUD+BWj2/oBA/AD9yqcMM3+rmo/hQEQHORneYSIiojrAwIKIrkj6V7x2YwdMu64l/rtJJtk7hTMZ+Xhp8QG8s/wIJvaOwd3XxKr+GnXCNxTofJshlZYAiVugP/I7Sg4thduFI+qzSjKUrYwypYKMYUDcAMDTn3ebiIjIRhhYEJFJZISox69rqSbV+2lHEj5bl4CEtFw18d5n605gbKfGuL9/HFpH1OHDu4t06u4LfdM+SO/4D4R5FEJ7fAVw5A8gYY1hlKkd8wxJ4wI06Qk0GwQ0HwREdTWsT0RERFbBf1WJyCyebi6Y2CsGt/VoqvpfSD8M6Y8hwYakvi2CcU/fOAxqFQatNSfcM0VAE6D7PYYko0ydWm8YzvboH8CFE8DpTYa0+jXAIwCIu9YQZDQfDAQ1q9tjJSIicjIMLIioViRoGNYuQqWdpy+qGoyl+85hw7F0lWTCPRlp6pZu0XXTD6O6UaZaDDWkkf8CLp4EjpeNLCUT8xVklo06tdiQPzAGmuaD4BHcFfAbA/gE1f0xExEROTAGFkRksS5NG2HOxEZIupin+mHIELUy4d6MRfvx5h+HMaFHE9zZuynsOo92o1ig+2RD0pUCZ3cBJ1Yagg3pk5FxCprtc9EIc6H/8zFDU6m4/oZajSa9AHcfex49ERFRvcfAgoisJrqRN54Z1QaPDonHzzuS8OWGkziRlotP1yWokaUGNA/Eg4Nd0TMu2HbzYZhC6wJEdzOk/k8ChdnAyQ3QH1+J0qPL4XrxOHBmmyGtn22YoE/my5BJ+mL7lQUa3vY7fiIionqIgQURWZ00fbqzT6zqi7HmSCq+2JCAdUfTsOpYBlYd24IOjQNwT79YXN8hCu6uNpwPw1QefkCrEdDHD0Na1xSEeRZDm7DWMBN4wjogK+nSaFMyd4bWzRBoxJUFGtE9Adc6GhWLiIionmJgQUQ27YcxqHWYSgfPZuKjlYew7NAF7D2TiccX7MZrvx3CXb1jcHuvpgj29ag/d8K/MdBloiHp9aqZlAowTq43BBtZZ4DEzYa09t+Aizs0jbvBN7QL0Hoo0LSnIVghIiJqQBhYEFGdaBXhh2eGxuCFGzphwbYkfLXxJFKyC/HWn0fwn1XHcGPnxpjcL7Zuh6s1hTTZkv4ZkrreaQg0LiaUBRnrDQFH9lloTm+Cr4w4tf0DQKMFwtsDTXsbmk017QMENLb3mRAREdkUAwsiqlNBPu6YMqgF7r+2mRpFSvpe7EnKxIJtiSrZdbhaUwMNGZpWUte7DIHGhRPQJaxD4eHl8EzdDU3GaSB5jyFt/eTSULjlgUZvIKS1vc+EiIjIqhhYEJFdSN+KGzo3xthOUdhx+iK+WH+y0nC1cSE+uLN3DG7pHg1/T7f6e5ck0AhuDjSKQ2b0CHiEhUGTkwyclqZSWwzzZiTvBTITgb2SfjCs5uGHRmGdgOb9gSY9DKNQeQXa+2yIiIhqjYEFEdmVjA7VLSZIJeNwtd9sPa1m9X5p8QE1XO1NXRvjrj6xaBnuIP0W/KOA9jcZkijMMYwwJcGGpKS/oCnMhkfiekCSUXA8EN3d0DFckjSnkhGpiIiIHAD/xSKiejlc7f92nsG8TSdx5HwOvt58WqVrmgerAGNomzC4utSD0aRM5eELNBtoSEJXCl3yPuQc+BN+GQegObPD0G8j/agh7f7WkE86hUd0hF+jNkC8jD7Vw9AEy55D9RIREdWAgQUR1cvhau/oHYOJvZpi04l0zNt4Cn8cSMbG4+kqNQ70wsTeTTGhR1PVZ8PhyDwaER2Qpw2HrzSd0mqB3HTg7A4gSebP2G5I+RegObMNPlLbse+/hnW9GhmaTUV1BiI7AREdDR3LiYiI7IyBBRHV62ZS1zQPUelMRj7mbz6F7/5KVO/fWHYY7yw/ijEdo3BXn6YIr8fdMEziEwzEX2dIomz0KV3SNuQfXQfviwehObcHyL8IHF9hSEaeAYaaDf8WQPPeQFQXILiFIYAhIiKqIwwsiMghSC3FUyNaq2ZSS/acw1ebTqrRpH7akaRS+wgf3NO/GNd3jIKHqxM8UBtHnwqMRXZYf3hJzYauBEjZb6jNkCDj3G4g5QBQkAnNyXXwwTpgz5eG9d28DX00pFYjsiMQ3gHQB9n7rIiIyIkxsCAih+Lp5oKbu0WrDt27EjMwb9MpLN5zFvuSczHt+z1q0r3bejZVs35HBDjZbNiu7obaCElGpcVA6iHozu5C/okt8M48Ak3yPqA4F0jaakjS+gpAuMbFMIJVWFsgvF3Za1sVvECaYxEREVmAgQUROWwzqS5NG6n09IhW+Hz1ISzan47krEL8Z+UxfLD6OEa0i8BdfWLQMy5I5XdKLm6qvwbC2iE76jpDzQYMc2uoGo1zu1Tthj55DzTSjCrtiCEdWHhpG24+QFjrskCjvSHYCGtnaJ5FRERkIgYWROTwQv08cE+vSDwxqgNWHErF3I0nsTXhApbsPadS6wg/3NknBuM6N1Ydw52e9K0IiTekDreoRfrSUqSe3IeQ0hRoUw8amlCd3w+kHjbUbhg7jFfkEwaEtgJCWpa9yjZbGYbTddZAjYiIaq0B/AtLRA2Fm4sWozpEqnTwXJZqJvW/nUk4lJyN5/63D7N+O4QbuzRWI061inCQOTGsRaOBziccCOsAtCzrIC5KSwy1G+f3lQUbBwz9OC6eBHJTDOnkusrbcve9FGSEtjQEHvI+KM5Qg0JERA0SAwsickptIv0x66YOeHpEa/ywPRHfbDmNE2m5+O/mUyp1j2mkAoyRHSKco7N3bbm4GoIDSSib0M84qZ/UZqimU/J61PBZgpCiHODsTkOqSCbzaxRn6McRJMn4vhkQ0ISjVBEROTkGFkTk1AK83XDftc1wb784NQfG15tlTozz2HbqokovLXbH37pHY2LPGEQ3crLO3pZO6hfdzZAqKikyTOZXHnQcKXt/1NCkyjjJX1VaN8N8G8ZAo2IKbMqgg4jICTCwIKIGQTpv920RotL5rAIs+CsR3249jXOZBfh4zQmV+seHYHTrAIwLDoE7R0mqeWQq6W8hqSKZdyPrjCHAkFqNSikBKC28QtDhCgREGwKMwBhDahRz6bNvOEetIiJyAAwsiKjBCff3VPNh/GNgc6w8lIKvt5zG2iOpWHs0TaW3155Rs3pP6NlE5SUTSGduCQ4kNR9U+TtdKZB1Frhw3BBopMtrQtnnsqBD+nRIqo6LBxDYpELgIa9NAf/GQEBjwC8SkKF0iYjIrhhYEFGD5eqixbB2ESqdSs9VM3sv+MtQi/H28iN4b+VRDGsbrvpiXNM82HmHrK2LUapUYNAEaDaw8nc6HZB9Dsg4DWScAi6euvReUuaZstqOY4ZULQ00vuEI8gqFJjjWEGxI0CGjV0mgI68SfLBjORGRTTGwICICEBPsg6dHtsbEToHYnlKKb7Ykqj4YS/clqxQX4oPxPZrglm7RCPH14DWzFmlyJoGApJg+l38vEwBKEysJNioGHZlJhiQ1IbpiaHKS4Z6TDKTurWFHGkOTKhV0RBkCD98wwDfCsNwv3PDqHcz+HkREtcTAgoioAg9XLcZ1jsBNXZvgUHIW5m8+jf/tPIOEtFy8vvQQ3vrjMK5rG65m9+7bPARaLWsxbMqlrNO3pLhqvpcaj7w06DKSkJl4AAGaHGizzxoCDqntkKCkLPiABB6Sqs7XUZE0qfIJNQQdfhFlwYcEHRGXL3P3seWZExE5HAYWREQ1aB3hj5fHtVc1Gb/uPotv/0rE7sQM/LY3WaUmQV6qL8bfVC2GO6+jvWo85EHfOwSFrlFAWNjlHb3Lgo/yGg5jsJGTUhZsyOt5IDdNZhK8FIAk77nyvl29AJ8QVcuh8Q5GgNYHmqBowNewTI7J8L2kIMAzkJ3QicipMbAgIroKma17Qk/pzN0UB85m4bu/DLUYiRfy8e/fD2P2n0cwpHUYRsT7YUxI6GXPtVRPgg9JjbvWnE+aXUlwUTHYyD5veDUuy5bX80BJAVCSD2QmqiT1Vl5XOw6pDVEBR7Ah4PAKBLwaGZKn8X01y9xYM0JEjoGBBRGRGdpG+eOlG9rjmZFtsGTvOXy39bTqiyFzY0h6c00Sxndvilt7RCMy4KqPmlTfml35RxrSlcjQuoXZhlqQvAsqGNHlpiIn5RT8NAXQ5BuWIS/dkCc3HSjKNtSGGGczTzX9sDQaF4R5+EMjtR7GoEMCEM8AwMMP8PA3vF722f/Se3ZcJ6I6wMCCiKgWvNxdVEduSUfOZ6s5MX7aloizGYYRpd5dcQSDWoWpvhgDW4WqEajIScjoYPLQLkkm+BM6HfJSUuAbFgZNdVVWJYWGQKM84EgH8i8C+RmG14KMyz9L0FJaCI2+FJoCWXax9sfs6gWNpz9CXLyg8QmqEnhI8jX0GXH3Bdy8L7139wZcveGSXQB46S7lk5G+iIiqYGBBRGShluF+eOH6Nri7SyM1otSCv5KwJeECVhxKUSnc3wO3dm+Cv3VrgqbB3rzeDZGrR9loVFHmrVecD13eBVw4cxxBXlpoCzMvBR+FWUBBlqH2RJbLa/nnstfiPMN2SvKhyck3/KOfWcN8ITWQMCn0svPxMgQd5QGIT1lAYghGNK5e8CvWQ+PfyPCdqyfg5lX26g24eRq24eIB15x8AFFlQYyXIZ8kF3dDEEdEDoOBhQl0Op1KNSktLUVxcbHFN0P2UVRUhLy8PGit3EjbWtu2ZDvmrmtOflPyWiuPJdzc3ODiYp9f+uTc9Hr9Fcuyo+7fWtu2ZDuyjruLBmM7RmJc58Y4kZqDBduS8NP2JJzPKsR/Vh5TqU+zINXZ+7o2YSbty9Rjulo+e99/W7P3+dls/y4e0PmEo7iRFrrQUPM7f0u/kaIcFXDoCjKRmXwagd4u0BRdCkI0EoRInqI8oCi37H2uISgp+6wvyoWmOA8afdn5Sf8SSVLzUg0JB0zpGSJnE1LDd3rZSlkwonHzQojGDRoPH+glSJOgo+Kr1h3+JTrAJwB6t7LlLh7V5634KpMvuhpfK37nBmjdyl5dy/K6ARptgyvbzv73t7nrmZNfZ0Jea+WxJXP2y8CiGnPmzFFJAgaRmpqKgoKCy/LJTc7JyVHfWWviLLl5WVlZVtmWrbZtyXbMXdec/KbktVae2pIy4+npCV9f3zqfbE3OKzMzUx2DLYIme+7fWtu2ZDtV1/UFcG+3INzZKRBrT2Rg8f50bDmVhU0nLqjk467FgFgf3NQlB+0iai4Pph7T1fLZ+/7bmr3Pr/6Xby/oNB7I9PZAYUCAWdsx7j/A3x9afbEKMLQSZJTkGYKN4nxoissCD6kZkYCkOB9FuRnwdAW00pxLmoLJd+p9QXlSneCL86HVFUFTWra8LHjRQG8IbiTlX/mBRf70VFcXaO2/ZfUSWGhdoZdgQ+tW9nr5Z8MyQ1BSU57qP7sYgheNfC57hQbaomLkevlA4+IKvQwCIEnrYniv1nEpPzZZX6+tkqf8fYU8GtmfYR/Vr6MtSxqU6jXIysqGvrQI2rL1rVWbZO+/v81dz5z8OhPyWiuPLWVnZ5ucV6OXo6RqycNlQEAALl68CH9//8u+T05OVjc6NDQU3t7eVnlQlJoP+VXbFqy1bUu2Y+665uQ3Ja+18tSG/FGTmhAJVKVcRUREoC7JX0yybymv9nrwstX+rbVtS7ZjyrpnMvLx844z+HF7EhIvSvMPg5bhvqoWY1znKARXmXzP1GO6Wj57339bs/f5OXP5rs16tS638kgic44U55cHHao5WFEeMtKS0cjHExpdEVAqSYIVw6u+uBC5WRfg4+GmghTjcpVPgpoKeS//rprX0iJodCUmX6OGSgUqlZIxINFUWV4xaSoEUFI2NCjV6eHiKoFYbbZnSHIs0uLA3cMDmrLtGoIf47rG96j0nQRuBYVF6kc/Q/+o6tYr61tlCHeRn18IL3nuq5ofZe8r5M3LLyh7RqxwHhXWkTBang28fXyr2f+l7WQXaeA74GG7/P0mz8ONGjVSz7zVPQ9XxBoLE8hNrHojpTZDLnBYWBiCg4NhrQdPV1dXlaz9a7a1tm3Jdsxd15z8puS1Vh5LGAPQlJQUhIeH13mzKNl3deXZGfZvrW1bsp2rrdskyAdTh7bEI4Pjsel4GuZtOIbVxzJw5HwOXv3tEN74/TCGtA5XM3xfGx9S3uHb1GO6Wj57339bs/f5OXP5rs16tS638vei9MGoSKdDiXsKNGFhNQbOV+w8XxsqyCkxNCeTYKe0pOy1us8V8xk/F13huxq2oSs1jB4mr/JZXwp9aQkK8nPh6e6mOvIbvjPmM66ju/S+bD01f0v5ezPyGZu7mUDVLpmR35YPo/Ivtmct1zOn55upzfyMeaX2+krkKcAPV+fnGwXNoEft8vebOftkYFFLxj4V8qBIZA5jmZEyZK/+FmRfMlt3n+bBaO5XCk+/Rli8Nxk/bEvE7qRMLNufrJJ0+L65azRu7trY5H/EiMiK5Mcl6Vdh56F69dIMJiUFHtYMmq4aUElT8LJXvfQzLUFqSgpCQ4KhleWSRwUVZcHKZcn4vQQtpVf8XldajIsXL6CRNNHT6C/PU+36FbdRdoylJcjOyoKfny+08kSvGuSUHavhQlZZZnjV6XXIyc6Gr493NevJq1q5fH29Xofc3Bz4yA+Fhg1fOhZU3rahpUIuvL28DPUO1exftpefnwcvqTGpdr+GPAV6j6vPlVMPMLCwUF23kyfHxzJDFfl7ueGO3jEqHUrOwvd/JeF/Ow0dvj9YfVylLo19cVvvIlzfMUpN1kdEZMN/pACXsr9njEGVdB52zzPMlWLt4EanQ7FXChAWZtm2dTrkp6TAz9ztVKjtMmU9CfRyUlLgbUKgp9fpkJ2SAq8r5JU8WSkp8LxKHrUd1H/OWR9OROSAWkf448UxbbHl2aH4cGJXDGoVqn5B23kmB0/9tBfdX1mOaQt2YcOxNOh07B5HRET1C3/6IiKqZ9xdtRjZIVKlsxfzMG/dYSw7nIGT6Xn4eecZlaICPDGuS2Pc3C0azUOv1oqXiIjI9hhYNEADBw5E586d8c4779j7UIjoKiICPHF3z0g8eX1H7ErKwk87krB491mczSwobyrVuUkgbuoShV5RbgjjFSUiIjthUyiq0cyZM3HHHXeo97GxseUjd8hQY3Fxcbj11luxcuVKm19B6fz04osvIjIyEl5eXhg6dCiOHj161fVkLhI5bhlCrlevXti6dWul72X+kSlTpqhRvWReiZtvvhnnz5+vlOf06dO4/vrrVYdrGQHsySefREnJpSEIz507h9tvvx0tW7ZU1+axxx6z4pkTXSJ//rrFNMJrN3bA1ueGYs7tXTG4dRhctBrsSszAi4sOYPSne/CP+Tuw/MB5FJc650R4RERUfzGwoBr98ssvGDt2bPnnl156CWfPnsW+ffvw1VdfITAwUD3kv/rqqza9im+88Qbee+89fPTRR9iyZQt8fHwwfPjwaictNFqwYAGmTZuGGTNmYMeOHejUqRNGjBihhnk1evzxx/Hrr7/ihx9+wJo1a9S5SXBRcUhhCSpkXOyNGzeqc547d64KcowKCwvV+OvPP/+82gdRXfB0c8H1HSPxxd09sPmZIXj++jZoE+mH4lI9lu0/j/vmbUPv11Zg5q/7se+MYVIlIiIiW2NgYUVqWLGiErskSx4clixZoiZsmz9/fvmyxMRE7N+/Xz2MG/n5+alJ3Zo2bYr+/fvjk08+wQsvvKAetA8fPgxbkPN699131YP7DTfcgI4dO2LevHkqCFi4cGGN682ePRv3338/Jk+ejLZt26qgRGodJDAQMgfJ559/rvINHjwY3bp1w5dffqkCCAlexB9//IEDBw7g66+/Vk3HRo4ciZdfflnVhEiwIaRGRI7vrrvuUteQqK6F+nngvmubYckj/fDfiW1wb79YhPh6ID23CF9uOInR/1mPke+uw6frTiA1x1BuiYiIbIF9LKwov7gUbV/8HfZw4KXh8HY3/3Z+8803ePDBB9Xr6NGjy5cvWrRI9cW42gyLU6dOVQ/bUrvx1FNPVZtHti8P51eSk5NT7fKEhAQ1w7nUjBjJA7w0bdq0aRMmTJhw2Try0L99+3Y888wz5cukmZJsY/PmzeqzfC/zSFTcbuvWrVXQJHn69u2rtt+hQwc1kZ2R1JQ89NBDKujq0qXLFc+JqK7Fh3qjb7tYPDOyDdYeTcVP28/gzwPncSg5G7OWHlZjrvdpflZ1+h7RPgL+nvYdn5+IiJwLA4sGTH55f+6551RzoAEDBlT6TgIFqSG4mqCgINX34OTJkzXmkSZU06dPV7UP0j/BnFmtjX0eKj7cGz9LwFGdtLQ01Yyp6jpynAcPHlTvZV13d3fVnKum7cprdfs1fkdUX8mM3YNbh6uUmVeMxXvP4ucdZ7D91EVsPJ6u0vML92FomzCM69wYA1uFqZGoiIiILMHAwoq83FxUzUFt1ebBu+K+zfHjjz+q/gYbNmxAjx49Kn2XlZWl+hxIUyFTj/tKxysP9JIsOT8iqp0AbzdM7BWD23o0wa6jidiQVIhfdp/DsZQc/LY3WaUALzeM6hCJcZ2j0CM2SM0MTkREZC4GFlYkD8u1aY5kpB68taiTB29pxiOdmr/44gt079690v6WLl2q+iU0adLkqttJT09HamqqGiWqJpY0hTLWEEjNhYwKZSSfpd9DdUJCQuDi4nLZCE8SSBm3J31FpMlURkZGpVoLWUe+M+apOpKUcZvGPESOJCrAA1Pim+DhwfHYfzYLv+w6g192nUVKdiG+3XpapcaBXhjbOUrVZLSK8LP3IRMRkQNh3XcD1bx5c6xatUo1eXrkkUdq1QxKSMdl6b8wbty4KzaF2rVrF3bu3Im//vpLvcrniqkmErDIQ/yKFSsq1ahIB+s+ffpUu440cZLO2BXX0el06nPv3r3VZ/nezc2tUh7pgC7DyxrzyPb37t1baSSpP//8U/U7kcCLyFHJDwntGwfguevbYtMzQzD/vl74W7do+Hm44kxGPj5cfRzD31mLEe+sxUdrjuNsRr69D5mIiBwAaywaMJl7QYIL6aQttSQyYZ40VZIaC+kTUVV2drbqW5Cfn69GjZJRpD777DPMmjULLVq0sElTKMknHcRfeeUVxMfHq0BDRqKKioqqFMxIJ2wZGvfRRx9Vn2Wo2UmTJqnamJ49e6pzy83NVcuMHcDvvfdelU/6iUiwIAGWBBPSMVwMGzZMBRB33nmnGvJWzl1Gp5K5Lzw8PMr3bQyMpNZFam/kswQ3DD7IEcg8GH1bhKj08rj2WHEwBQt3ncHqwymq0/frSw/hX8sOoVdcEMZ2MnT6DvJxt/dhExFRPcTAooFr1aqVmuROggtpPjRq1Cg1WVzXrl0vyyvDykqSh2apRZBf9uUX/0GDBtn0GGW0qby8PDzwwAOq6VK/fv2wbNkyNfGd0fHjx1WnbaPx48erh3w5XgkIpNmUBEwVO2O//fbbqrZF5q6Q+ShkxCfp0G4k12Px4sVqFCgJOGT+DAlMpAamooqjQ8loUzLCVkxMzBU7tBPV5/kxJGXkFWHJ3nP4ZedZbD15AZtPGNILv+xDvxYhGN0xEsPaRaj+GURERIKBRQO0evXqSp/btGlT3ndAfvEfM2bMZesYH5Lt0QFb9iMP81Uf6KsOS1txRmzx8MMPq2RkPHYjCUwkkKgYTFTNIwHCb7/9dsXj4+Rj5IwCvd1Vp29JSRfz8Ovuc1i856zqm7HmSKpKz/1vH/q3DMWYTpEY2iYcPh78J4WIqCHjvwJUSfv27Wvsu0BEDVN0I288NLC5SidSc7B4zzn8uvssjqbkYPnB8yp5umkxpHW4qskY1DpM1X4QEVHDwsCCKpHmRkRENWkW6otHh8SrdDg5WwUYUpNxMj1PNZ2S5OPuguvahmNMpyhcGx/KOTKIiBoIBhZERFQrMhxtq4hWeGJYS+w7k6UCDKnNkJGlFu46q5K/p6vq8D26YxSuaR6sJu8jIiLnxMCCiIgs7gfVITpApX+OaI2diRmqJuO3vefUHBnfb0tSqZG3G4a1jcCIDhHo2zyENRlERE6GgQUREVmNzNrdLaaRSi+MboutCRdUTcbSfcm4kFuEBdsSVfLzdMV1bcIxskMkro0PYZ8MIiInwMCCiIhsNkdGn+bBKs0c204FGRJgLNufjNTsQvy884xK0idjcJtwjGofgQGtQuHtzn+aiIgcEf/2JiIi2/9j46LFNS1CVPq/se2w4/RF1VRq2b5knMssUE2nJMnoUoNahal+GUPahMOXQ9gSETkMBhZERFTnNRk9YoNUeuH6ttidlKECjN/2nUPihXxVqyHJ3VWL/vGhGNk+AkPbhnMyPiKieo6BBRER2bVPRpemjVR6emRrNQHf0n3nsHRvMk6k5ZbPk+HmIs2qQjCsbbgayjbc35N3jYionmFg0QANHDgQnTt3xjvvvGPvQyEiqjS6VPvGASpNH9YKR87nqOZSEmjI+7VHUlV6fuE+dGoSqIIMSS3CfHkViYjqAQ4oTjWaOXMm7rjjDvU+NjZW/aOv1Wrh7++PuLg43HrrrVi5cqXNr6Ber8eLL76IyMhIeHl5YejQoTh69OhV15szZ446bk9PT/Tq1Qtbt26t9H1BQQGmTJmC4OBg+Pr64uabb8b58+cr5Xn00UfRrVs3eHh4qGCMiOqG/H0j82Q8fl1L/PH4AKx4YoAayrZr00D1/e7EDPz798O47u21GPzWGsxaegi7z+agVKfnLSIishMGFlSjX375BWPHji3//NJLL+Hs2bPYt28fvvrqKwQGBqqH/FdffdWmV/GNN97Ae++9h48++ghbtmyBj48Phg8frgKDmixYsADTpk3DjBkzsGPHDnTq1AkjRoxASkpKeZ7HH38cv/76K3744QesWbNGnZsEF1Xdc889GD9+vM3Oj4iurnmoLx4a2Bw//6Mvtj47BK/d2AGDWoXC3UWLhLRcfLouAX///jB6z1qJf/64BysOnkdBcSkvLRFRHWJTKGvS64HiPMvWLykBdK7yc51567p5m79OmSVLluD222/HBx98gIkTJ6pliYmJ2L9/v3oYN/Lz80NERARKSkrQrFkzDBgwQNUiSG3CLbfcglatWsEWtRXvvvsunn/+edxwww1q2bx58xAeHo6FCxdiwoQJ1a43e/Zs3H///Zg8ebL6LEGJnOfcuXPx7LPPIjMzE59//jm++eYbDB48WOX58ssv0aZNGxW89O3bVy2TgEakpqZiz549Vj8/IjJfmL8nbu/VVKWcwhKsOZyKPw4kY8WB80ivMFeGl5sLBrQMxbB24RjcOgyB3u683ERENsTAwpokqHgtqtarS1jgVtuVnz0LuPuYvZo8WD/44IPqdfTo0eXLFy1apPpiSLOnK5k6dSpefvllVbvx1FNPVZtHtv/1119fcTs5OTnVLk9ISEBycrKqGTEKCAhQTZs2bdpUbWBRVFSE7du345lnnilfJk24ZBubN29Wn+X74uLiSttt3bo1mjZtqvIYAwsiqt9kONrrO0ZiZPtwnD13Hgm5Llh+MAV/HDivhrGVOTMkGUaiaqQCDElSAyLNrYiIyHoYWDRg0gfhueeeU82BpPahIgkUjDUEVxIUFISwsDCcPHmyxjzShGr69Omq9kFqO1xdXU3+B93Y50FqKCqSzxJwVCctLQ2lpaWXrSPHefDgQfVe1nV3d1fNuUzdLhHVb64uGvRtEYJrW4apuTJkhKk/9ierIONQcjY2n7ig0mu/HULTIO/yIKNXsyB4uLrY+/CJiBweAwtrkuZIUnNQS7V58K60bzP8+OOPqr/Bhg0b0KNHj0rfZWVlqT4H0lTI1OO+0vHKA70ki86PiKiWI0xNG9YKp9PzsOLQeaw8lIItJy7g9IU8zN14UiVvdxcVkAxpHYZBrcM4lC0RUS0xsLAmeViuRXOkSn0stCWAay36WJipS5cuqlPzF198ge7du1d60F+6dCnatm2LJk2aXHU76enpqv+BjBJVE0uaQhlrHaTmQvpzGMnnmkZpCgkJgYuLy2UjPEkgZdye9BWRJlMZGRmVai1kHfmOiJxL02BvTO4bp1JuYQnWH0vDqkMpKtBIyS7EnwfOqyTaRfmXBxmdogPVXBtERHR1DCwaqObNm+Ott95S/SjkIfz99983uxmUkI7V0n9h3LhxNmkKJQGLPOivWLGiPJCQGhXpYP3QQw9Vu440cZIhYmUd43HpdDr12biOfO/m5qaWGUeCOnz4ME6fPo3evXubdGxE5Jh8PFwxvF2ESvL3kjSZkgBDkswCLp8lvbfyGIJ93DGwlaHJ1LUtQ+DvWeuecERETo+BRQPWsmVLrFq1SgUX8rAvE+bJg7/UWEggUFV2drbqf5Cfn69GjZo/fz4+++wzzJo1Cy1atLBJUyjJJx3EX3nlFcTHx6tA44UXXkBUVFSlYEY6YcvQuDLvhJChZidNmqRqY3r27KnOLTc3Vy0zdgC/9957VT7pJyKd1B955BH06dNHdQw3OnbsmKpNMZ73rl271HKp0ZEAhoicp8nUo0PikZZTiNWHU1VthkzGJ6NM/bQjSSXpAC7zaPSPD0X/lqHo0DiAtRlERBUwsGjgZIhYmeTOWHMxatQoNVlc165dL8srw8pKkgdqqUWQX/blF/9BgwbZ9BhltKm8vDw88MADqulSv379sGzZMjXxndHx48dVp20jmXdCmmjJ8UpQILUdEjBV7ND99ttvq9oWqbEoLCxUc2NIh/aK7rvvPtXfpGITMuNoVTL5HhE5lxBfD9zSLVql4lId/jp5obzJ1PHUXPx18qJKb/15BI283dBPgoz4EPRrEaxG9iMiasgYWDRAq1evrvRZ5m4w9keQX/zHjBlz2TrGUZ/s0QFb9iPNqSTVRB705bgqevjhh1UyMh67kQQmEkhUDCaq5ql6rYio4XBz0eKa5iEqPXd9WyReyMPao6mqJmPjsXRczCvGr7vPqiSah3hhcJsLqjajR2wQPN040hQRNSwMLKiS9u3bq+ZARERUWZMgb0zsFaOS1GbsSsxQQYakPWcycTwtH8fXJahZwD3dtOgVF6yCjAEtQzhvBhE1CAwsqBJpbkRERFevzZBaCUlPDGuFtOwCLN1xArvOF6kRp85nFWLNkVSVXgYQFeCpgoxrWkgNSLBqckVE5GwYWBAREVkoyMcd17UKwsRrw1TzzcPns7HuSJpqOrUl4QLOZhbgu78SVRKtI/zKmlkFqwn6/DjaFBE5AQYWREREViSBResIf5Xu798M+UWl2JKQjnVH07DxeDoOnstSM4FL+mJDghptqmN0gAoy+jYPQdeYRuyfQUQOiYEFERGRDXm5u6i5MCSJ9JxCbDqRroKMjcfScDI9DztPZ6g0Z9VxeLhq0T22UXmNRttIP94fInIIDCyIiIjqULCvB0Z3jFJJnMnIx4Zjadh0PF29ykzgG47J+3T1vbe7CzpE+KBf62z0aRaMjtGBcHfV8p4RUb3DwIKIiMiOGgd64dbuTVSSIa+Pp+aooGLj8TRsPnEBmfnF2HI6SyUhI051bdpIjTol/TM6Nwlk0ykiqhcYWBAREdWj/hktwvxUmnRNLHQ6PQ6ey8TyPadxIK1ITc53IbfI0IzquKFGQ2ovJLjoHReEXs2CVdAhza+IiOoaAwsiIqJ6SqvVoE2kP4JdwhAWZhhx6mhKDracSMfmhAvYcuIC0nIKsTXhgkpYeQyuWg3aNw5A95hGqq9Gt5gghPpxeFsisj0GFg3QwIED0blzZ7zzzjv2PhQiIjKDBBYtw/1UurNPrGo6dSItVwUYMvKUvCZnFajJ+yR9tj5BrRcb7I1uMY0Q38gFg9p7IT7cXwUtRETWxN5fVKOZM2fijjvuUO9jY2PVP2harRb+/v6Ii4vDrbfeipUrV9r8Cso/nC+++CIiIyPh5eWFoUOH4ujRo1ddb86cOeq4PT090atXL2zdurXS95988okKsuR85NwyMjJseBZERNYnf3c1D/XF7b2a4t0JXbDpmcFY99QgvD2+Eyb2aopW4X7QaKBGnvppxxm8vuI0hr+7Hl1f+RP3zv0LH6w+pmo6CopLeXuIyGKssaAa/fLLL3j66afLP7/00ku47777kJeXh6SkJMyfP1895L/88st47rnnbHYl33jjDbz33nv46quvVEDzwgsvYPjw4Thw4IAKGqqzYMECTJs2DR999JEKKqR2ZsSIEdi3bx+iogwjsch5yDJJzzzzjM2On4ioLgONJkHeKt3YJVotk87fO05fxLaEC9h49DwOpuQhI68YKw6lqCTcXAxNrqSvRpemgejcpJGq5ZDtERGZioGFlX9Zzy/Jt2j90tJSuOhdzP7L3MvVq9b/ACxZsgS33347PvjgA0ycOFEtS0xMxP79+9VDt5Gfnx8iIiJQUlKCZs2aYcCAAaoWQWoTbrnlFrRq1QrWJtfk3XffxfPPP48bbrhBLZs3bx7Cw8OxcOFCTJgwodr1Zs+ejfvvvx+TJ09WnyXAkPOcO3cunn32WbXsscceU6+rV6+2+nETEdUXAV5uGNQqDAPiQ3Bn50A0Cg7BweQcbDt5AdtOXsS2UxdVP409SZkqzdt0Sq0X6O2mAo1O0QGI9QMG+AYiyLf6H3OIiAQDCyuSoKLXN73sUrK23L4F3m7eZq/3zTff4MEHH1Svo0ePLl++aNGi8mZCVzJ16lRVYyG1G0899VS1eWT7X3/99RW3k5OTU+3yhIQEJCcnq5oRo4CAAFULsWnTpmoDi6KiImzfvr1SLYQ04ZJtbN68+YrHQUTk7NxcDKNISbrvWsMPOEkX87FT+mXIRH2JF7H/bJaq1Vh9OFUlZeExNAvxMayrajUC0SrCDx6uHIGKiAwYWDRg0gdBmjD9+uuvqvahIgkUjDUEVxIUFKRGKjl58mSNeaQJ1fTp09U/XlLb4erqanLtyvnz59Wr1FBUJJ8l4KhOWlqaqvmpuo4c58GDB03aLxFRQ2w+NbaToaloUYkOB89lYefpiyrg2H4yHUkZhaqjuKSfd54pb0IlwUWHxgFqJKqOjQPRMsIXbuwYTtQgMbCwImmOJDUHFjeFcqldUyhz/Pjjj0hJScGGDRvQo0ePSt9lZWVhzZo1+Pzzz00+7isdrzzQS6pNYEFERHVP5sboJM2gmgTiLp1O/Xvh6hOIvWeyVKAhAcfeM5mqVmPfmSyVgMRKwUaLRu7oGV+oZgqXUaw4WziR82NgYUXysFyb5khG6sFbUzcP3l26dMGOHTvwxRdfoHv37pX2t3TpUrRt2xZNmjS56nbS09ORmpqqOlXXxJKmUMZaB6m5kP4cRvJZhsytTkhIiArOjLUdRvIPY9VaDCIiMk2QjzsGtQ5TSRibUEmAoVKS4VU6ixuCDWDhvjSV191Fi9aRfmgXFYC2kX5oG+WP1hH+8HLj4JREzoSBRQPVvHlzvPXWW6ofhTyEv//++2Y3gxLSsVr6L4wbN84mTaEkYJEO4ytWrCgPJKRGZcuWLXjooYeqXcfd3R3dunVT6xiPS6fTqc81rUNERLVvQjWqg+GHH/l7PvFCPvYkXcSWo+dw/GIJ9p3JRFZBSXnn8EvrA02DvNE8yB2dY7LQrnGAGpkqMsCTtdpEDoqBRQPWsmVLrFq1SgUX8rAvQ7LKg7/UWEggUFV2drbq15Cfn69GjZLhZj/77DPMmjULLVq0sElTKMknHcRfeeUVxMfHlw83K0PGVgxmpGP22LFj8eijj6rPMtTspEmTVG1Mz5491bnl5uaqZUZyLpKOHTumPu/duxe+vr5q23K8RERkHvk7u2mwN6IbeaJ7uEv5bOEq2DiTgQNns3DgXJbqv3E+qxCn0vNUWnn00jxCMhpVmwh/VashgUabSD80C659awAiqjsMLBo4GSJWJrkz1lyMGjVKPVx37dr1srwyrKwkqRGQWoTevXurWoBBgwbZ9BhltCmZc+KBBx5Qk9j169cPy5YtqzSHxfHjx1WnbaPx48erJlpyvBI8SG2HBEwVm0LJELQyCaBR//791asES/fcc49Nz4mIqKEFG5JGdzR0DhfpOYXYfzYTfx09i8RsPQ6ey8ax1BzVb2PTiXSVjKQveHSgB9pEBapJ/1pG+KnX2BAfNcoVEdUPDCwaoKrzNrRp06a8P4L84j9mzJjL1jGO+mSPDtiyH2lOJakmMiytHFdFDz/8sEpGxmM3+r//+z+VKqqah4iIbCPY1wP9WoSgpb9O1WxIs1qZAfxYSo6q1TDWbhw5n62CjdMXC3H64nn8vv9S/znpKN4sxBfx4b6I8tGga3O96rshzbNcODIVUZ1jYEGVtG/fHn369OFVISKiOufp5qKGrZVU8QeflKwCbDmciJRCFxxLycXh89k4kpyN3KJS9V6SsumsepERqOKCfdA8zAfNQ33LU7NQH/h48NGHyFb4p4sqkeZGRERE9YXUWof6eaBnU//ymg1jwHE2s0AFGIeSs7D3VCpOZ5bgaEoOCkt0lQOOCqRzuAQYFQMOCUAi/NlpnMhSDCyIiIjIIQOOxoFeKg1oGYKUFF/DQCHQ4GxGvuqvcTwlB8dTc3E8NQcnUnOQllOEc5kFKm04dqkPh/B00yI22Acxwd7qVfpvGN9L0KFl0yqiq2JgQURERE5D+lYYh8Ed1KryCH8ZeUXlgYYh2DC8l5GpCop1OKRqPy6v5fBw1aogI0YCjvJXHzVcbmSgJzuQE5VhYEFEREQNQqC3O7rFSGpUaXlRiQ5nMvJxMj0Xp9JycTI9z/A+PQ+JF/JU06oj53NUqkoqMiIDvNC4kReiG3mhSSMZblfeS3DjpWo7XDlyFTUQDCyIiIioQVOdvUN8VEKryt+VlBqDDplzIxcn0wyvCem5auZxY1AiaWtC9TUo0q+jasAhr1GBngj3Z40HOQ8GFkREREQ1PSi5SDMo6W/hAyC00nc6nR5pOYVIvJiPpIt5KtAwvkpNhwQbxaX6suX51W5fRm4P9fVQwUdEgKeq/aj4PtzPHdpSHe8POQQGFkRERES1IB26w/w9VaravMoYeKRkS+AhwUYeki7kl703vCZnFqjAQ/JI2p2UWeO+QnzdVaAhAUeUCjy8EO7vgTA/TzVqVpifh5q1vK7mmCKqDgMLIiIiIhsFHhIISOoRG1Rt4JGeW6QCjHOZ+UjOMoxYdS4jX70aP0tzKxnRStLeMzUHHzJhoNR+SKBhSJeCDuMyeR/i66HmDCGyNgYWDdDAgQPRuXNnvPPOO/Y+FCIiogYdeBgf+DtEX5oUsKLS0lIcPX0OJW6+SM4qxLmsAiRn5uNcRkFZTUcBUrMLcTGvWNV+yNwekq4mwMtN1YIE+3ggyMcdQeq9u+G9z6Xlwb7uCPDk4yKZhiWFajRz5kwcPXoUX3/9NWJjY3Hq1Cm13NPTE+Hh4ejZsycefPBBDB482KZXUSZBmjFjBj799FNkZGSgb9+++PDDDxEfH3/F9ebMmYN///vfSE5ORqdOnfDee++ha9eulQKsNWvWXDZB4Pvvv2+zcyEiIjKHNG0K9HJFWJg/2kcbJgesjqFWo1AFGRJwpFYIOiouk1RUqkNmfrFKMvyuKXzdXRDs51EWfJS9+roj0MtNBSnSDMtfXr3c1XtZ5u3uwqZZDQwDC6rRL7/8gqeffrr880svvYT77rsPeXl5SEpKwvz58zF06FC8/PLLeO6552x2Jd944w0VFHz11VeIi4vDCy+8gOHDh+PAgQMqyKnOggULMG3aNHz00Ufo1auXqp0ZMWIE9u3bh6ioqPJ8999/vzovIy8vL5udBxERkS1HtooK9FLpaj/WZeWXGIKOnEJcyC1SKT2n6NL73EvLJen0QE5RKXLUyFh5Jh+TNM0K8HJHgJerGurXGIQEeBsCEFnu7+UKFOYhpkCWe8DP01UlLzcGJY6IgYUVyR9WfX6+RevrSkqgc3U1O8LXeHnV+leBJUuW4Pbbb8cHH3yAiRMnqmWJiYnYv3+/ehg38vPzQ0REBEpKStCsWTMMGDAAkZGRePHFF3HLLbegVasqY/RZgVyTd999F88//zxuuOEGtWzevHmqxmThwoWYMGFCtevNnj1bBQ2TJ09WnyXAkPOcO3cunn322fJ83t7e6pwq7k/Oj4iIyBnJs4I82EuKD/e7an7pB3IxrxBHTycDHr7IyC9W/UIu5EgAUoSsspoPWS4TEGbmlyAzv0g1y5IktSiSgKvVjBy7bJheXw9XlYzBhp+nW/lnX/nscWmZr4cLSvJz0KTYA/5e7up7H3dXNaM6O7TXHQYWViRBxeGu3WAPrXZsh8bb2+z1vvnmG9WcSV5Hjx5dvnzRokWqqZC/v/8V1586daqqsZDajaeeeqraPLJ9aU51JTk5l086JBISElRTJqkZMQoICFC1EJs2bao2sCgqKsL27dvxzDPPlC/TarVqG5s3b66UV2pd5NgkuBgzZowKYNzd3a94rERERA2pH0gjb3fEBnkiLCxI/Xt6NfIjXV5RqSHgyJOgwxCAGN6XBSJ58lqkXtOy8pBbDGQXFCOnsETVkJTq9OXNtcxzpNIn+c3V280FXu6u8PFwUTUhPh6uqpmWJFUzUlqEkMAL8Dbmkdey72VZxVcvNy1yCksRWKKDh5uGQUsVDCwaMOmDIE2Yfv31V1X7UJEECsYagisJCgpCWFgYTp48WWMeaWo0ffr08toAVzNqZM6fP69epYaiIvksAUd10tLSVGe3quvIcR48eLD8s9TSxMTEqKZRe/bswT//+U8cPnxYNaMiIiKi2pF/4+XhXdLVmmbpdDqkpKSof6MlaDEGJRJgSKCRXVBS9r4EOQUlyCoLPuS98Tu1rKAEF3MLUFCiVw/+slzo9UBuUalKadX/hlkmzezzlFnXZXQtlVy16tVDfdbC01Xea6EpLUaA7zl4ubvAw7Usr3xfto6HqxaF+TkIS9HBy8NVfZZ8hlfDe1etBF0lCEP9x8DCiqQ5ktQc1FZtHrwr7tscP/74o/qDvGHDBvTo0aPSd1lZWapT8+eff27ycV/peOUvC0mWnJ8tSEdtow4dOqhmXUOGDMHx48dt0qyLiIiITA9KZFZyU1UNUKQJV0FJKXILS5GvAosS5KlUtqy4RL3mFZYg5WImNG6eyC+Wz6WGPEUlZeuVlq8nefOKS1Wwovaph2F5UelVji7DhDOo+QdaEerrhi3PXuojWl8xsLDyH4baNEcykgdvbUkJtHXw4N2lSxfs2LEDX3zxBbp3715pf0uXLkXbtm3RpEmTq24nPT0dqampqlN1TSxpCmWsdZCaC3nwN5LPMmRudUJCQuDi4lJe22Ekf+FUrcWoSJpXCQYWREREjt+Ey9B8ydWsgORqpEVE0rnz8G8UjKJSPQqKS1FQrCt7LUVByaX3+UUlSLuYCTdPbxSWVMhbYvi+sFinApys3ALoNS5q3cISnRrhq7CktOzV8Nndxf4/yJqCgUUD1bx5c7z11luqH4U8hFccYtXUZlBCOlbLH8Rx48bZpCmUBCzS/2HFihXlgYTUqGzZsgUPPfRQtetIH4lu3bqpdYzHJX9xyOea1hG7du1SrxU7cxMREREZyfOLNFGS0a2uFojoVNDiccWgxZTARvKcS678Y2l9xcCiAWvZsiVWrVqlggt52JchWeXBX2osJBCoKjs7W/VryM/PV6NGScfnzz77DLNmzUKLFi1s0hRK8kkH8VdeeUXNW2Ecblb6RVQMZqRj9tixY/Hoo4+qzzLU7KRJk1RtjMy3IeeWm5urlhlrJaTD+qhRoxAcHKz6WDz++OPo378/OnbsWIurSURERGQbMkqWI2Bg0cBJX4KVK1eW11zIg7avr2+lieSMZFhZSVIjIL/q9+7dW9UCDBo0yKbHKKNNydwZ0idCJsjr168fli1bVmkOCwkUpNO20fjx41UTLTleCYaktkMCJmNTKDmH5cuXlwcc0uzr5ptvtul8HERERETOjIFFA7R69epKn9u0aVPeH0F+8ZdhV6syjvpkjw7Ysh9pTlVxIrvqhqWtOv/Eww8/rFJ1c1RIIFF11u2qeYiIiIjIdAwsqJL27dujT58+vCpEREREZBYGFlTjEKxERERERKa6+rhaREREREREV8HAwkLSJp+IZYaIiIgaOgYWteTm5qZeZbQiInMYy4yxDBERERE5A/axqCUZmjUwMFBNaiK8vb0tHiXJliMuWWvblmzH3HXNyW9KXmvlqS3ZtgQVUmak7EgZIiIiInIWDCwsYJyh2RhcWOPBU2ZXlJkXbfFQa41tW7Idc9c1J78pea2Vx1ISVHB2byIiInI2DCwsIA+ekZGRalbp4uJii2+GPNCmp6ermaCvNk28vbZtyXbMXdec/KbktVYeS0jzJ9ZUEBERkTNiYGEF8qBojYdFeaiVB0+ZUdoWgYU1tm3Jdsxd15z8puS1Vh4iIiIiuhyfnIiIiIiIyGIMLIiIiIiIyGIMLIiIiIiIyGLsY2HC5HdZWVmoC9K+Pzs722Z9LKyxbUu2Y+665uQ3Ja+18jgqe58by3ftroe18tn7/tuavc/Pmct3bdZj+Xaesm3rY3C08t0Qn02yyp6DTZkUmoHFFchNFE2aNLHWvSEiIiIicsjn4oCAgCvm0ehNCT8aKIkQz549Cz8/P5vNaVBVjx498Ndff9XrbVuyHXPXNSe/KXmvlkeicgkkExMT4e/vD2djy/Jl7/07c/m2Rj5nL9uC5dt216g267F8O0/Zdva/v/lscmUSKkhQERUVddUaE9ZYXIFcvOjoaNQlGbbWVv/oW2vblmzH3HXNyW9KXlO3J3mc8eHLluXL3vt35vJtzXzOWrYFy7ftrlFt1mP5dp6y7ex/f/PZ5OquVlNh5HwNbR3clClT6v22LdmOueuak9+UvLa8vo7A3ufP8l2762HtfM7K3ufvzOW7NuuxfDtP2bb1MTha+eazSc3YFIqoSnMRicozMzPt/usQkTWxbJMzY/kmZ5blQM8mrLEgqsDDwwMzZsxQr0TOhGWbnBnLNzkzDwd6NmGNBRERERERWYw1FkREREREZDEGFkREREREZDEGFkREREREZDEGFkREREREZDEGFkREREREZDEGFkRmyMvLQ0xMDKZPn87rRk4lNjYWHTt2ROfOnTFo0CB7Hw6R1SQkJKgy3bZtW3To0AG5ubm8uuQUDh8+rP7ONiYvLy8sXLjQrsfkate9EzmYV199Fb1797b3YRDZxMaNG+Hr68urS07l7rvvxiuvvIJrr70WFy5ccIi5AIhM0apVK+zatUu9z8nJUT8QXXfddbAn1lgQmejo0aM4dOgQRo4cyWtGROQA9u/fDzc3NxVUiKCgILi68jdVcj6LFi3CkCFD4OPjY9fjYGBBDcLatWsxZswYREVFQaPRVFtVOGfOHBXte3p6olevXti6dWul76X506xZs+rwqInqrnzLegMGDECPHj0wf/58XnpyirItPwhJLZxso2vXrnjttdfq+AyIbPt3t9H333+P8ePHw94YWFCDIG1qO3XqpP6AVmfBggWYNm0aZsyYgR07dqi8w4cPR0pKivr+l19+QcuWLVUicrbyLdavX4/t27erX73k4WvPnj11eAZEtinbJSUlWLduHT744ANs2rQJf/75p0pEzvJ3t8jKylJNWUeNGgW70xM1MFLs//e//1Va1rNnT/2UKVPKP5eWluqjoqL0s2bNUp+ffvppfXR0tD4mJkYfHBys9/f318+cObPOj53IFuW7qunTp+u//PJLXmxy+LK9ceNG/bBhw8q/f+ONN1Qicqa/u+fNm6efOHGivj5gjQU1eEVFReqX2qFDh5ZfC61Wqz7LL1xCmkAlJibi5MmTePPNN3H//ffjxRdfbPDXjpyjfMuvZtnZ2eUdAFeuXIl27drZ7ZiJrFW2pWmf/Lp78eJF6HQ61fSkTZs2vMDkFOW7vjWDEuzBRA1eWloaSktLER4eXulayGfprE3k7OX7/PnzuPHGG9V7ySuBszyQETl62ZaO2tK0r3///tJCA8OGDcPo0aPtdMRE1n82yczMVP0ufvrpJ9QHDCyIajF0IZEzadasGXbv3m3vwyCyCRnJj6P5kbMKCAhQPw7VF2wKRQ1eSEgIXFxcLvuDKZ8jIiIa/PUhx8byTc6KZZucWYiDPpswsKAGz93dHd26dcOKFSvKr4W0xZXPffr0afDXhxwbyzc5K5ZtcmbuDvpswqZQ1CBIh9Rjx46Vf05ISFCzVcpkSU2bNlXDuU2aNAndu3dHz5498c4776gOrZMnT7brcROZguWbnBXLNjmzHGd8NrH3sFREdWHVqlVqKLeqadKkSeV5/vOf/+ibNm2qd3d3V0O8bd68mTeHHALLNzkrlm1yZquc8NlEI/+zd3BDRERERESOjX0siIiIiIjIYgwsiIiIiIjIYgwsiIiIiIjIYgwsiIiIiIjIYgwsiIiIiIjIYgwsiIiIiIjIYgwsiIiIiIjIYgwsiIiIiIjIYgwsiIiIiIjIYgwsiIjIpgYOHIjHHnvsinliY2PxzjvvONyd+L//+z907tzZ3odBRFQvMLAgInJAd999NzQajUpubm6Ii4vDU089hYKCAnsfGhERNVCu9j4AIiKqnREjRuDLL79EcXExtm/fjkmTJqlA41//+hcvqYOTeyoBIxGRI2GNBRGRg/Lw8EBERASaNGmCcePGYejQofjzzz/Lv9fpdJg1a5aqzfDy8kKnTp3w448/qu/0ej1atGiBN998s9I2d+3apYKTY8eOqc/y/rPPPsONN94Ib29vxMfHY9GiRZXWWbNmDXr27KmOJzIyEk8//TRKSkpqPO6UlBSMGTNGHZMc2/z58yt9L8cmTYyaNm2qthkVFYVHH330qs2R/vvf/6omVQEBAZgwYQKys7Ov2NRK1pF1jeRcP/74Y4wePVqda5s2bbBp0yZ1LaQ5l4+PD6655hocP378smOQ9eQ+yHq33norMjMzK30v11C25+npidatW+ODDz4o/+7kyZNq3wsWLMCAAQNUnqrXhIjIETCwICJyAvv27cPGjRvh7u5evkyCinnz5uGjjz7C/v378fjjj+OOO+5QgYA8yN5zzz2qxqMi+dy/f38VdBjNnDlTPSzv2bMHo0aNwsSJE3HhwgX13ZkzZ9SyHj16YPfu3fjwww/x+eef45VXXrliM67ExESsWrVKBTrykC3BhtFPP/2Et99+Wz2sHz16FAsXLkSHDh2ueP7ysC/5Fi9erJKc4+uvv272dXz55Zdx1113qQBLAoDbb78df//73/HMM89g27ZtKuh5+OGHK60jgcf333+PX3/9FcuWLcPOnTvxj3/8o/x7CRJefPFFvPrqqzh48CBee+01vPDCC/jqq68qbUcCsqlTp6o8w4cPN/vYiYjsTk9ERA5n0qRJehcXF72Pj4/ew8NDL3+da7Va/Y8//qi+Lygo0Ht7e+s3btxYab17771Xf9ttt6n3Z86cUdvYsmWL+lxUVKQPCQnRz507tzy/bPf5558v/5yTk6OWLV26VH1+9tln9a1atdLrdLryPHPmzNH7+vrqS0tL1ecBAwbop06dqt4fPnxYrb9169by/AcPHlTL3n77bfX5rbfe0rds2VIdjylmzJihzjUrK6t82ZNPPqnv1atX+eeYmJjy7Rt16tRJrVvTuW7atEkt+/zzz8uXffvtt3pPT89K+5ZrmJSUVL5Mro3ci3PnzqnPzZs313/zzTeV9v3yyy/r+/Tpo94nJCSo/bzzzjsmnS8RUX3FPhZERA5q0KBBqoYgNzdX/cLv6uqKm2++ufxX9Ly8PFx33XWV1ikqKkKXLl3Ue2lidP311+OLL75QTZnkF/fCwkL87W9/q7ROx44dy99LcyB/f//yGgb5db1Pnz6qBsSob9++yMnJQVJSkmrOVJHkl+Ps1q1b+TKpGQgMDCz/LPuXZkvNmjVT/UikRkSaTsl6NZGmTn5+fuWfpUlWxVoQU1U81/DwcPVasbZElkkH+aysLHUdhJxj48aNy/PI9ZBmaIcPH1bHJLUp9957L+6///7yPNJUTJpsVdS9e3ezj5eIqD5hYEFE5KDkId/YZEmCA+lDIc2Q5CFWHuzFkiVLKj30Cum3YHTffffhzjvvVIGJNIMaP3686idQUdVOxBJEyIOzrUhfBXkoX758ueozIs2K/v3vf6vmTTV1aL7aMWq1WtWMqWoH6SttxxgsVbfM1PM33odPP/0UvXr1qvSdi4vLZfeTiMiRMbAgInIC8uD87LPPYtq0aapfQNu2bVUAcfr0adUhuCZSGyAPtFLzIf0D1q5da9Z+pUOy9ImQh3bjQ/eGDRvUL/XR0dGX5ZfaCfm1Xkaxkn4ZQoKIjIyMSvmkY7fUUkiaMmWKWm/v3r3o2rUraiM0NBTnzp0r/yw1DgkJCbAGucZnz55VNUBi8+bN6n60atVK1XDI8hMnTqi+KUREzoydt4mInIQ0IZJfwefMmaMe7KdPn646bEsnYWmOs2PHDvznP/+p1GlY8ktnaumcLCM+STMec0htgnTEfuSRR3Do0CH88ssvmDFjhgpw5OG6KnnYluZN0iF6y5YtKsCQWhMJJIzmzp2ral6kQ7o8kH/99dfq+5iYmFpfm8GDB6tRo9atW6cCFBmat2qNQW3JKE6yPem8LtuXEayks7uM2GXs/C4d6d977z0cOXJE7V9qh2bPnm2V/RMR1RcMLIiInIT0QZARi9544w3V70JGOJLRh+ShVmoW5IFemkbJEK8VSdMp6XsxefJks/cpzax+++03bN26VTXFevDBB9X2nn/++RrXkYdq+RVfalJuuukmPPDAAwgLCyv/XvpbSNMh6ashfR6kSZT0/wgODkZtSeAk+5OhZKVfiQzP27x5c1iDNEeT85Dan2HDhqljrjicrAROMtysnLf015DjkOCp6n0gInJ0GunBbe+DICIi+5Ff2YcMGaJqHowdlomIiMzFwIKIqIGSEaBSU1NVMx5ptsNJ2YiIyBJsCkVE1EB9++23qt+CdJyW5lNERESWYI0FERERERFZjDUWRERERERkMQYWRERERERkMQYWRERERERkMQYWRERERERkMQYWRERERERkMQYWRERERERkMQYWRERERERkMQYWREREREQES/0/bQjwdvtBt7QAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Moody diagram: f vs. Re for various roughness values\n", "Re_values = np.logspace(3.5, 7, 200)\n", "roughness = [0.0001, 0.001, 0.01, 0.05]\n", "\n", "fig, ax = plt.subplots(figsize=(8, 5))\n", "for k_D in roughness:\n", " f_vals = [f_colebrook(k_D, Re) for Re in Re_values]\n", " ax.loglog(Re_values, f_vals, label=f'k/D = {k_D}')\n", "\n", "ax.set_xlabel('Reynolds number')\n", "ax.set_ylabel('Friction factor f')\n", "ax.set_title('Moody diagram (Colebrook)')\n", "ax.legend()\n", "ax.grid(True, which='both', alpha=0.3)\n", "plt.tight_layout(); plt.show()" ] }, { "cell_type": "markdown", "id": "a081d613", "metadata": {}, "source": [ "### Available functions\n", "\n", "| Function | Description |\n", "|----------|-------------|\n", "| `f_colebrook(k/D, Re)` | Colebrook with laminar/turbulent transitions |\n", "| `f_colebrook_pure(k/D, Re)` | Pure iterative Colebrook (turbulent only) |\n", "| `f_barr_bathurst(k/D, Re)` | Barr-Bathurst approximation |\n", "\n", "> **Note:** For open-channel flow, use $Re' = 4 \\cdot Re$ (hydraulic radius convention)." ] } ], "metadata": { "kernelspec": { "display_name": "python311", "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.9" } }, "nbformat": 4, "nbformat_minor": 5 }