{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "9e47daa3", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import cantera as ct\n", "import pde\n", "\n", "import ipywidgets as widgets\n", "\n", "from matplotlib import pyplot as plt\n", "%matplotlib widget" ] }, { "cell_type": "code", "execution_count": 2, "id": "74655ac0", "metadata": {}, "outputs": [], "source": [ "grid = pde.CartesianGrid([[0, 0.14]], 32, periodic=False)\n", "field = pde.ScalarField(grid, 0)" ] }, { "cell_type": "code", "execution_count": 3, "id": "d447beea", "metadata": {}, "outputs": [], "source": [ "eq = pde.PDE({\"c\": \"0.0001 * laplace(c)\"}, bc=[{\"derivative\": 1}, {\"value\": 0}])" ] }, { "cell_type": "code", "execution_count": 4, "id": "9d3fc1f0", "metadata": {}, "outputs": [ { "data": { "application/json": { "ascii": false, "bar_format": null, "colour": null, "elapsed": 0.01196742057800293, "initial": 0, "n": 0, "ncols": null, "nrows": 29, "postfix": null, "prefix": "", "rate": null, "total": 1, "unit": "it", "unit_divisor": 1000, "unit_scale": false }, "application/vnd.jupyter.widget-view+json": { "model_id": "f882e9d7d6044daf90c18abcf20eb901", "version_major": 2, "version_minor": 0 }, "text/plain": [ " 0%| | 0/1.0 [00:00, ?it/s]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "solution = eq.solve(field, t_range=1, dt=0.001)" ] }, { "cell_type": "code", "execution_count": 5, "id": "9564e6fc", "metadata": {}, "outputs": [ { "data": { "application/json": { "ascii": false, "bar_format": null, "colour": null, "elapsed": 0.010971546173095703, "initial": 0, "n": 0, "ncols": null, "nrows": 29, "postfix": null, "prefix": "", "rate": null, "total": 1, "unit": "it", "unit_divisor": 1000, "unit_scale": false }, "application/vnd.jupyter.widget-view+json": { "model_id": "372e9e77cf9f415197ff441271fd1c8a", "version_major": 2, "version_minor": 0 }, "text/plain": [ " 0%| | 0/1.0 [00:00, ?it/s]" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "application/vnd.jupyter.widget-view+json": { "model_id": "906c8eb3a30b4a5e897e8e930d9e2c6b", "version_major": 2, "version_minor": 0 }, "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAABDeElEQVR4nO3de1xU953/8fcwDIwXQIQIXoiCGi81RoHEYiRu2ixe0jZm/a1mm5p0t03L1qwim9YYk6ZNmxAbt+svjZc10fySplE3pTa2Na3YRKqRpJECSZQkWi8QhCAmMhoUGDi/P3BGEUSEGc4M5/V8dB6Bw3fOfM73MZ15+z3n+z02wzAMAQAAwDJCzC4AAAAAPYsACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxYSaXUAwa25u1vHjxxURESGbzWZ2OQAAoBMMw9Dp06c1ZMgQhYRYcyyMANgNx48fV0JCgtllAACALigvL9ewYcPMLsMUBMBuiIiIkNTyBoqMjDS5GgAA0Bkul0sJCQne73ErIgB2g+e0b2RkJAEQAIAgY+XLt6x54hsAAMDCCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEBEQDXrFmjxMREOZ1OpaSkaPfu3R22z8/PV0pKipxOp5KSkrRu3bpWf9+/f7/mzp2rESNGyGazadWqVR3uLycnRzabTVlZWd08EgAAgMBnegDcsmWLsrKytHz5chUVFSk9PV2zZs1SWVlZu+2PHDmi2bNnKz09XUVFRXrooYe0aNEi5ebmetvU1dUpKSlJTz75pOLj4zt8/XfeeUfr16/XxIkTfXpcAAAAgcr0APjzn/9c3/rWt/Ttb39b48aN06pVq5SQkKC1a9e2237dunW69tprtWrVKo0bN07f/va39W//9m9auXKlt82NN96op556SnfddZfCw8Mv+9pnzpzR3XffrWeffVbR0dE+PzYAAIBAZGoAbGhoUGFhoTIyMlptz8jI0N69e9t9TkFBQZv2M2bM0L59+9TY2HhVr79w4ULdfvvtuu222zrVvr6+Xi6Xq9UDAAAg2JgaAGtqatTU1KS4uLhW2+Pi4lRVVdXuc6qqqtpt73a7VVNT0+nX3rx5swoLC5WTk9Pp5+Tk5CgqKsr7SEhI6PRzAQAAAoXpp4CltvfiMwyjw/vztde+ve2XU15ersWLF+tXv/qVnE5np+tctmyZamtrvY/y8vJOPxcAACBQhJr54rGxsbLb7W1G+6qrq9uM8nnEx8e32z40NFQxMTGdet3CwkJVV1crJSXFu62pqUl/+ctf9Mwzz6i+vl52u73N88LDwzu8ptBX/vh+pba/V6Vbx16jOycP8/vrAQAAazF1BDAsLEwpKSnKy8trtT0vL09Tp05t9zlpaWlt2u/YsUOpqalyOBydet0vf/nLeu+991RcXOx9pKam6u6771ZxcXG74a8nfVB1WttKjmvvoZOm1gEAAHonU0cAJSk7O1sLFixQamqq0tLStH79epWVlSkzM1NSy2nXiooKvfjii5KkzMxMPfPMM8rOztZ9992ngoICbdiwQZs2bfLus6GhQQcOHPD+XFFRoeLiYvXv31+jRo1SRESEJkyY0KqOfv36KSYmps12M4weFCFJOlh9xuRKAABAb2R6AJw/f75Onjypxx57TJWVlZowYYK2b9+u4cOHS5IqKytbrQmYmJio7du3a8mSJVq9erWGDBmip59+WnPnzvW2OX78uCZPnuz9feXKlVq5cqWmT5+uXbt29dixddV1cf0lSYeqz1zxekgAAICrZTM8Myhw1Vwul6KiolRbW6vIyEif7bfB3azxP/yj3M2G9j74JQ0Z0Mdn+wYAwOr89f0dTAJiFjBaCwsNUWJsP0nSR5+cNrkaAADQ2xAAA9Toi04DAwAA+BIBMEB5JoIwAggAAHyNABigPCOAzAQGAAC+RgAMUNfFtYwAHvrkjJinAwAAfIkAGKBGxPRTaIhNp+vdqnKdM7scAADQixAAA1RYaIhGeGcCcxoYAAD4DgEwgI0edP46QCaCAAAAHyIABrDR568DPMgIIAAA8CECYADzjgBWMwIIAAB8hwAYwK67aASQmcAAAMBXCIABLDG2n+znZwJ/4qo3uxwAANBLEAADWFhoiEbE9JXEHUEAAIDvEAADnPc0MHcEAQAAPkIADHAsBQMAAHyNABjgRjMCCAAAfIwAGOBGx7WMAH70yWlmAgMAAJ8gAAY470zgc25Vn2YmMAAA6D4CYIALD7VrODOBAQCADxEAg8B1g7glHAAA8B0CYBDwXAfILeEAAIAvEACDwOg4RgABAIDvEACDgGctQGYCAwAAXyAABoGka/opxCa5zrl1gpnAAACgmwiAQSA81K4RMf0kSR9xGhgAAHQTATBIMBEEAAD4CgEwSIw+vxQMI4AAAKC7CIBBwjMCeIgRQAAA0E0EwCBx8QggM4EBAEB3EACDhGcmcO3ZRp04w0xgAADQdQTAIOF02DX8/ExgFoQGAADdQQAMIp4FoQ9+wnWAAACg6wiAQcQzEeSjakYAAQBA1xEAg8h15+8JfIhTwAAAoBsIgEFklOeewNXcExgAAHQdATCIjLymv0Js0qm6RtWcaTC7HAAAEKQIgEHE6bDr2oF9JTERBAAAdF1ABMA1a9YoMTFRTqdTKSkp2r17d4ft8/PzlZKSIqfTqaSkJK1bt67V3/fv36+5c+dqxIgRstlsWrVqVZt95OTk6MYbb1RERIQGDRqkOXPm6MMPP/TlYfnF6PPXAR5kIggAAOgi0wPgli1blJWVpeXLl6uoqEjp6emaNWuWysrK2m1/5MgRzZ49W+np6SoqKtJDDz2kRYsWKTc319umrq5OSUlJevLJJxUfH9/ufvLz87Vw4UK99dZbysvLk9vtVkZGhj7//HO/HKeveJaC+YgRQAAA0EU2w+TZBFOmTFFycrLWrl3r3TZu3DjNmTNHOTk5bdovXbpU27ZtU2lpqXdbZmamSkpKVFBQ0Kb9iBEjlJWVpaysrA7rOHHihAYNGqT8/Hzdcsstnard5XIpKipKtbW1ioyM7NRzuuu3RRXK2lKsm0YM1P9mpvXIawIA0JuY8f0daEwdAWxoaFBhYaEyMjJabc/IyNDevXvbfU5BQUGb9jNmzNC+ffvU2NjY5Vpqa2slSQMHDuzyPnoCM4EBAEB3mRoAa2pq1NTUpLi4uFbb4+LiVFVV1e5zqqqq2m3vdrtVU1PTpToMw1B2dramTZumCRMmXLZdfX29XC5Xq0dPGzWov2zMBAYAAN1g+jWAkmSz2Vr9bhhGm21Xat/e9s66//779e6772rTpk0dtsvJyVFUVJT3kZCQ0KXX645WM4GruQ4QAABcPVMDYGxsrOx2e5vRvurq6jajfB7x8fHttg8NDVVMTMxV1/Af//Ef2rZtm9544w0NGzasw7bLli1TbW2t91FeXn7Vr+cLowednwnMHUEAAEAXmBoAw8LClJKSory8vFbb8/LyNHXq1Hafk5aW1qb9jh07lJqaKofD0enXNgxD999/v37zm9/o9ddfV2Ji4hWfEx4ersjIyFYPM3juCcwIIAAA6IpQswvIzs7WggULlJqaqrS0NK1fv15lZWXKzMyU1DLqVlFRoRdffFFSy4zfZ555RtnZ2brvvvtUUFCgDRs2tDp929DQoAMHDnh/rqioUHFxsfr3769Ro0ZJkhYuXKiXX35Zr776qiIiIryjilFRUerTp09PdsFVuy7OsxQMI4AAAODqmb4MjNSyEPTPfvYzVVZWasKECfrv//5v71Is3/zmN3X06FHt2rXL2z4/P19LlizR/v37NWTIEC1dutQbGCXp6NGj7Y7oTZ8+3bufy10v+Pzzz+ub3/xmp+o2axr5+xW1+sov9mhgvzD97ZF/7LHXBQCgN2AZmAAJgMHKrDfQ2YYmjX/0jzIMad/Dtym2f3iPvTYAAMGOABggs4BxdfqE2ZUQ7bknMKeBAQDA1SEABqnrmAgCAAC6iAAYpEaxFAwAAOgiAmCQujATmBFAAABwdQiAQcqzGPShakYAAQDA1SEABinPPYFPft6gk2fqzS4HAAAEEQJgkOoTZtew6JYFqw8yCggAAK4CATCIXeedCMJ1gAAAoPMIgEFslHcpGEYAAQBA5xEAg5hnBJCZwAAA4GoQAIPY6PMjgMwEBgAAV4MAGMRGDWoJgDVnGvTp5w0mVwMAAIIFATCI9Q0LvTATmNPAAACgkwiAQe66uPPXAXIaGAAAdBIBMMiNPn8a+BAjgAAAoJMIgEFutGcE8BNGAAEAQOcQAIOcZwSQtQABAEBnEQCD3IWZwPX6jJnAAACgEwiAQa5feKiGDuCewAAAoPMIgL3AdecXhOaOIAAAoDMIgL2AZyIIdwQBAACdQQDsBTwTQRgBBAAAnUEA7AU8I4BcAwgAADqDANgLeGYCnzhdr1N1zAQGAAAdIwD2Av2ZCQwAAK4CAbCXGM1MYAAA0EkEwF7Ce0cQbgkHAACugADYS1yYCMIIIAAA6BgBsJdgBBAAAHQWAbCX8IwAVp+uV21do8nVAACAQEYA7CX6h4dqSJRTEqeBAQBAxwiAvYhnFPAjTgMDAIAOEAB7Ee91gIwAAgCADhAAe5HrPDOBGQEEAAAdIAD2IqPiGAEEAABXRgDsRTyngD9xMRMYAABcHgGwF4lwOjQsuuWewPuP15pcDQAACFQEwF7mhmEDJEklHxMAAQBA+wIiAK5Zs0aJiYlyOp1KSUnR7t27O2yfn5+vlJQUOZ1OJSUlad26da3+vn//fs2dO1cjRoyQzWbTqlWrfPK6wWDisChJ0rsfnzK3EAAAELBMD4BbtmxRVlaWli9frqKiIqWnp2vWrFkqKytrt/2RI0c0e/Zspaenq6ioSA899JAWLVqk3Nxcb5u6ujolJSXpySefVHx8vE9eN1hMPD8C+C4jgAAA4DJshmEYZhYwZcoUJScna+3atd5t48aN05w5c5STk9Om/dKlS7Vt2zaVlpZ6t2VmZqqkpEQFBQVt2o8YMUJZWVnKysrq1uu2x+VyKSoqSrW1tYqMjOzUc/ztTL1b1//oTzIMad/Dtym2f7jZJQEAEFAC8fu7p5k6AtjQ0KDCwkJlZGS02p6RkaG9e/e2+5yCgoI27WfMmKF9+/apsbFzM1+78rqSVF9fL5fL1eoRaPqHh2rkNS2zgTkNDAAA2mNqAKypqVFTU5Pi4uJabY+Li1NVVVW7z6mqqmq3vdvtVk1Njd9eV5JycnIUFRXlfSQkJHTq9Xqa5zrAknJOAwMAgLZMvwZQkmw2W6vfDcNos+1K7dvb7uvXXbZsmWpra72P8vLyq3q9nnJhJvApU+sAAACBKdTMF4+NjZXdbm8z6lZdXd1mdM4jPj6+3fahoaGKiYnx2+tKUnh4uMLDA/+augszgWuvGGoBAID1mDoCGBYWppSUFOXl5bXanpeXp6lTp7b7nLS0tDbtd+zYodTUVDkcDr+9bjAZNzhSDrtNn37eoI8/O2t2OQAAIMCYfgo4Oztbzz33nDZu3KjS0lItWbJEZWVlyszMlNRy2vWee+7xts/MzNSxY8eUnZ2t0tJSbdy4URs2bNADDzzgbdPQ0KDi4mIVFxeroaFBFRUVKi4u1qFDhzr9usHM6bBrbHzLrCaWgwEAAJcy9RSwJM2fP18nT57UY489psrKSk2YMEHbt2/X8OHDJUmVlZWt1uZLTEzU9u3btWTJEq1evVpDhgzR008/rblz53rbHD9+XJMnT/b+vnLlSq1cuVLTp0/Xrl27OvW6wW7isCi9V1Grdz8+pdsnDja7HAAAEEBMXwcwmAXyOkL/+065fpD7rr6YNFCbv5NmdjkAAASMQP7+7immnwKGf0xMaJkI8n6FS83NZHwAAHABAbCXGnVNf/Vx2HWm3q3DNWfMLgcAAAQQAmAvFWoP0YShLcPaxSwIDQAALkIA7MUmnl8QmlvCAQCAixEAe7EbEgZIkkpYCgYAAFyEANiL3XD+jiClx11qcDebXA0AAAgUBMBe7NqBfTWgr0MNTc36sOq02eUAAIAAQQDsxWw2m64f2jIKWMJ1gAAA4DwCYC93AxNBAADAJQiAvdzE89cBck9gAADgQQDs5TwzgT/65LTqGtzmFgMAAAICAbCXi4t0Ki4yXM1Gy23hAAAACIAWwHWAAADgYgRAC2BBaAAAcDECoAVcmAhyytxCAABAQCAAWsDEoQMkScdO1ulUXYO5xQAAANMRAC0gqq9DI2L6SmI5GAAAQAC0jIlMBAEAAOcRAC3Ccx1gcTkjgAAAWB0B0CI8M4EZAQQAAARAi/jCkEiF2KTq0/Wqqj1ndjkAAMBEBECL6BsWquviIiRJJYwCAgBgaQRAC+GOIAAAQCIAWsrEBM+C0EwEAQDAygiAFuIZASwpPyXDMMwtBgAAmIYAaCFj4iMUFhoi1zm3jp6sM7scAABgEgKghTjsIRo/OFIS1wECAGBlBECLueH8gtAlLAgNAIBlEQAthlvCAQAAAqDFeO4I8v7xWrmbms0tBgAAmIIAaDFJsf0UER6qc43NOlh9xuxyAACACQiAFhMSYtOEoZ71AE+ZWwwAADAFAdCCPAtCl7AgNAAAlkQAtKCLF4QGAADWQwC0oInnl4L5sOq0zjU2mVwNAADoaQRACxo6oI9i+oXJ3WzoQKXL7HIAAEAPIwBakM1m844CvstpYAAALCcgAuCaNWuUmJgop9OplJQU7d69u8P2+fn5SklJkdPpVFJSktatW9emTW5ursaPH6/w8HCNHz9eW7dubfV3t9uthx9+WImJierTp4+SkpL02GOPqbnZGmvjedYDfJeJIAAAWI7pAXDLli3KysrS8uXLVVRUpPT0dM2aNUtlZWXttj9y5Ihmz56t9PR0FRUV6aGHHtKiRYuUm5vrbVNQUKD58+drwYIFKikp0YIFCzRv3jy9/fbb3jYrVqzQunXr9Mwzz6i0tFQ/+9nP9NRTT+kXv/iF3485EHgngrAUDAAAlmMzDMMws4ApU6YoOTlZa9eu9W4bN26c5syZo5ycnDbtly5dqm3btqm0tNS7LTMzUyUlJSooKJAkzZ8/Xy6XS6+99pq3zcyZMxUdHa1NmzZJkr7yla8oLi5OGzZs8LaZO3eu+vbtq1/+8pedqt3lcikqKkq1tbWKjIy8ugM32ckz9Ur56U7ZbNK7j2YowukwuyQAAHpEMH9/+4qpI4ANDQ0qLCxURkZGq+0ZGRnau3dvu88pKCho037GjBnat2+fGhsbO2xz8T6nTZumP//5z/roo48kSSUlJdqzZ49mz57d7eMKBjH9wzV0QB8ZhvReBaeBAQCwklAzX7ympkZNTU2Ki4trtT0uLk5VVVXtPqeqqqrd9m63WzU1NRo8ePBl21y8z6VLl6q2tlZjx46V3W5XU1OTHn/8cf3Lv/zLZeutr69XfX2993eXK7hn0N6QEKWKU2dVUl6rqSNjzS4HAAD0ENOvAZRaZqVezDCMNtuu1P7S7Vfa55YtW/TSSy/p5Zdf1t/+9je98MILWrlypV544YXLvm5OTo6ioqK8j4SEhCsfXACbeP46QG4JBwCAtZgaAGNjY2W329uM9lVXV7cZwfOIj49vt31oaKhiYmI6bHPxPr///e/rwQcf1F133aXrr79eCxYs0JIlS9q97tBj2bJlqq2t9T7Ky8uv6ngDjXcpGGYCAwBgKaYGwLCwMKWkpCgvL6/V9ry8PE2dOrXd56SlpbVpv2PHDqWmpsrhcHTY5uJ91tXVKSSk9eHb7fYOl4EJDw9XZGRkq0cwu35olGw2qeLUWdWcqb/yEwAAQK9g+ing7OxsPffcc9q4caNKS0u1ZMkSlZWVKTMzU1LLqNs999zjbZ+Zmaljx44pOztbpaWl2rhxozZs2KAHHnjA22bx4sXasWOHVqxYoQ8++EArVqzQzp07lZWV5W3z1a9+VY8//rj+8Ic/6OjRo9q6dat+/vOf68477+yxYzdbhNOhkdf0l8RpYAAArMTUSSBSy5ItJ0+e1GOPPabKykpNmDBB27dv1/DhwyVJlZWVrdYETExM1Pbt27VkyRKtXr1aQ4YM0dNPP625c+d620ydOlWbN2/Www8/rEceeUQjR47Uli1bNGXKFG+bX/ziF3rkkUf0ve99T9XV1RoyZIi++93v6oc//GHPHXwAmDgsSoeqz6ikvFZfGtv+aXcAANC7mL4OYDDrDesIvbD3qB7dtl+3jrlGz//rTWaXAwCA3/WG7+/uMv0UMMx18UQQ/i0AAIA1EAAtbtzgSIWG2HTy8wZ9/NlZs8sBAAA9gABocU6HXWMHR0hiORgAAKyCAAgWhAYAwGIIgNAN568DLCEAAgBgCQRAeEcA369wqbmZiSAAAPR2BEBo9KD+6uOw60y9W4drzphdDgAA8DMCIBRqD9GEoS3rIJWUMxEEAIDejgAISRdOA3MdIAAAvR8BEJKk5GujJUlvH/7U5EoAAIC/EQAhSUobGSObTfrwk9OqPn3O7HIAAIAfEQAhSRrYL0xfGNJyHeCbh2pMrgYAAPgTARBe00ZdI0nafZAACABAb0YAhFf66FhJ0p6DNTIM1gMEAKC3IgDCK2V4tMJDQ1R9ul6HqlkPEACA3ooACC+nw66bEgdK4jQwAAC9GQEQrUwbdf40MBNBAADotQiAaOXm8wHwrcMn1djUbHI1AADAHwiAaGX84EjF9AtTXUOTispOmV0OAADwAwIgWgkJsWmq5zTwwRMmVwMAAPyBAIg20s8HwN1cBwgAQK9EAEQbN59fD7Ck/JRc5xpNrgYAAPgaARBtDB3QR0mx/dRsSAV/P2l2OQAAwMcIgGjXtIvuCgIAAHoXAiDa5VkP8E2uAwQAoNchAKJdXxwZI3uITYdrPlfFqbNmlwMAAHyIAIh2RTodumFYlCSWgwEAoLchAOKypo2+RpK05xATQQAA6E0IgLisi68DbG42TK4GAAD4CgEQlzX52gHqF2bXp5836ECly+xyAACAjxAAcVkOe4i+mBQjSdrDbGAAAHoNAiA65FkPkOVgAADoPQiA6JDnOsC/HvlU5xqbTK4GAAD4AgEQHRo1qL/iIsNV727WvqOfmV0OAADwAQIgOmSz2TRtlGc5GE4DAwDQGxAAcUXTRnsmgrAgNAAAvQEBEFd08/nrAPcfd+nTzxtMrgYAAHQXARBXNCjCqbHxETIMZgMDANAbBEQAXLNmjRITE+V0OpWSkqLdu3d32D4/P18pKSlyOp1KSkrSunXr2rTJzc3V+PHjFR4ervHjx2vr1q1t2lRUVOgb3/iGYmJi1LdvX02aNEmFhYU+O67e5OK7ggAAgOBmegDcsmWLsrKytHz5chUVFSk9PV2zZs1SWVlZu+2PHDmi2bNnKz09XUVFRXrooYe0aNEi5ebmetsUFBRo/vz5WrBggUpKSrRgwQLNmzdPb7/9trfNZ599pptvvlkOh0OvvfaaDhw4oP/6r//SgAED/H3IQenm8+sB7j5YI8PgtnAAAAQzm2Hyt/mUKVOUnJystWvXereNGzdOc+bMUU5OTpv2S5cu1bZt21RaWurdlpmZqZKSEhUUFEiS5s+fL5fLpddee83bZubMmYqOjtamTZskSQ8++KDefPPNK442dsTlcikqKkq1tbWKjIzs8n6CQV2DW5N+nKeGpma98cA/KDG2n9klAQDQJVb6/r4cU0cAGxoaVFhYqIyMjFbbMzIytHfv3nafU1BQ0Kb9jBkztG/fPjU2NnbY5uJ9btu2Tampqfrnf/5nDRo0SJMnT9azzz7bYb319fVyuVytHlbRNyxUycMHSGI5GAAAgp2pAbCmpkZNTU2Ki4trtT0uLk5VVVXtPqeqqqrd9m63WzU1NR22uXifhw8f1tq1azV69Gj96U9/UmZmphYtWqQXX3zxsvXm5OQoKirK+0hISLiq4w126aPPrwd4kOVgAAAIZqZfAyi1LDZ8McMw2my7UvtLt19pn83NzUpOTtYTTzyhyZMn67vf/a7uu+++VqeiL7Vs2TLV1tZ6H+Xl5Vc+uF7EsxzM3r+flLup2eRqAABAV5kaAGNjY2W329uM9lVXV7cZwfOIj49vt31oaKhiYmI6bHPxPgcPHqzx48e3ajNu3LjLTj6RpPDwcEVGRrZ6WMn1Q6MU1ceh0+fcerei1uxyAABAF5kaAMPCwpSSkqK8vLxW2/Py8jR16tR2n5OWltam/Y4dO5SamiqHw9Fhm4v3efPNN+vDDz9s1eajjz7S8OHDu3w8vZ09xKapI1tC9psHuQ4QAIBgZfop4OzsbD333HPauHGjSktLtWTJEpWVlSkzM1NSy2nXe+65x9s+MzNTx44dU3Z2tkpLS7Vx40Zt2LBBDzzwgLfN4sWLtWPHDq1YsUIffPCBVqxYoZ07dyorK8vbZsmSJXrrrbf0xBNP6NChQ3r55Ze1fv16LVy4sMeOPRh5TgPvZiIIAADBywgAq1evNoYPH26EhYUZycnJRn5+vvdv9957rzF9+vRW7Xft2mVMnjzZCAsLM0aMGGGsXbu2zT5feeUVY8yYMYbD4TDGjh1r5Obmtmnzu9/9zpgwYYIRHh5ujB071li/fv1V1V1bW2tIMmpra6/qecHsaM0ZY/jS3xujHvqDceZco9nlAABw1az4/X0p09cBDGZWXUco/Wevq/zTs3r+mzfq1rGDzC4HAICrYtXv74uZfgoYwWfaqJblYHZzHSAAAEGJAIir5rkv8J5DrAcIAEAwIgDiqk0dGSObTfrokzP6xHXO7HIAAMBVIgDiqkX3C9P1Q6MkSW8yGxgAgKBDAESXeJaD2cN1gAAABB0CILok3XsdYI2YSA4AQHAhAKJLUkZEy+kIUfXpeh2sPmN2OQAA4CoQANEl4aF23ZTYcls4loMBACC4EADRZdNGtQTAPQdZDgYAgGBCAESXeRaEfvvIp2pwN5tcDQAA6CwCILpsbHyEYvuHqa6hSUVln5ldDgAA6CQCILosJMR2YTkY1gMEACBoEADRLZ4AyEQQAACCBwEQ3ZI+uiUAvvvxKdXWNZpcDQAA6AwCILplcFQfjR7UX82G9OcPPjG7HAAA0AkEQHTbV28YIknaWlRhciUAAKAzCIDotjmThkqS3jxUo2rXOZOrAQAAV0IARLddG9NXKcOj1WxI20qOm10OAAC4AgIgfGLO5JZRQE4DAwAQ+AiA8ImvXD9YoSE27T/u0kefnDa7HAAA0AECIHwiul+Y/mHMIEnSbxkFBAAgoBEA4TN3nj8N/GrxcTU3GyZXAwAALocACJ/58rhBiggPVcWps3rn6KdmlwMAAC6DAAifcTrsmn39YEnSb4s5DQwAQKAiAMKnPLOBf/9upc41NplcDQAAaA8BED41JXGgBkc5dfqcW7s+rDa7HAAA0A4CIHwqJMSmOyaxJiAAAIGMAAif88wGfuODEzpV12ByNQAA4FIEQPjcmPgIjRscqYamZv3hvUqzywEAAJcgAMIv7pw8RBKLQgMAEIgIgPCLr90wVDab9M7Rz1T+aZ3Z5QAAgIsQAOEX8VFOTR0ZI0l6lTUBAQAIKARA+M2ci2YDGwa3hgMAIFAQAOE3MyfEKzw0RH8/8bner3CZXQ4AADiPAAi/iXA69I/j4ySxJiAAAIGEAAi/8qwJuK3kuNxNzSZXAwAAJAIg/OyW667RwH5hqjlTrzf/ftLscgAAgAIkAK5Zs0aJiYlyOp1KSUnR7t27O2yfn5+vlJQUOZ1OJSUlad26dW3a5Obmavz48QoPD9f48eO1devWy+4vJydHNptNWVlZ3T0UXMJhD9FXJw6WxJqAAAAECtMD4JYtW5SVlaXly5erqKhI6enpmjVrlsrKytptf+TIEc2ePVvp6ekqKirSQw89pEWLFik3N9fbpqCgQPPnz9eCBQtUUlKiBQsWaN68eXr77bfb7O+dd97R+vXrNXHiRL8do9XNOX8a+I/vV+nzerfJ1QAAAJth8vocU6ZMUXJystauXevdNm7cOM2ZM0c5OTlt2i9dulTbtm1TaWmpd1tmZqZKSkpUUFAgSZo/f75cLpdee+01b5uZM2cqOjpamzZt8m47c+aMkpOTtWbNGv30pz/VpEmTtGrVqk7X7nK5FBUVpdraWkVGRl7NYVuKYRi6deUuHT1Zp/+ef4PunDzM7JIAABbG97fJI4ANDQ0qLCxURkZGq+0ZGRnau3dvu88pKCho037GjBnat2+fGhsbO2xz6T4XLlyo22+/Xbfddlt3DwUdsNls3lHArUXHTa4GAACYGgBramrU1NSkuLi4Vtvj4uJUVVXV7nOqqqrabe92u1VTU9Nhm4v3uXnzZhUWFrY7yng59fX1crlcrR7oHM+i0HsOnlD16XMmVwMAgLWZfg2g1DJCdDHDMNpsu1L7S7d3tM/y8nItXrxYv/rVr+R0OjtdZ05OjqKioryPhISETj/X6kbE9tPkaweo2ZB+V1JpdjkAAFiaqQEwNjZWdru9zWhfdXV1mxE8j/j4+Hbbh4aGKiYmpsM2nn0WFhaqurpaKSkpCg0NVWhoqPLz8/X0008rNDRUTU1N7b72smXLVFtb632Ul5d36bityrMmILOBAQAwl6kBMCwsTCkpKcrLy2u1PS8vT1OnTm33OWlpaW3a79ixQ6mpqXI4HB228ezzy1/+st577z0VFxd7H6mpqbr77rtVXFwsu93e7muHh4crMjKy1QOdd/v1gxUaYtN7FbU6VH3a7HIAALCsULMLyM7O1oIFC5Samqq0tDStX79eZWVlyszMlNQy6lZRUaEXX3xRUsuM32eeeUbZ2dm67777VFBQoA0bNrSa3bt48WLdcsstWrFihe644w69+uqr2rlzp/bs2SNJioiI0IQJE1rV0a9fP8XExLTZDt+J6R+u6dddoz9/UK3fFh3XAzPGmF0SAACWZPo1gPPnz9eqVav02GOPadKkSfrLX/6i7du3a/jw4ZKkysrKVmsCJiYmavv27dq1a5cmTZqkn/zkJ3r66ac1d+5cb5upU6dq8+bNev755zVx4kT9v//3/7RlyxZNmTKlx48PrXlmA/+2uELNzaauQAQAgGWZvg5gMGMdoat3tqFJNz6+U2fq3XolM003jhhodkkAAIvh+zsARgBhLX3C7Jo5IV6S9Ju/MRkEAAAzEADR4/7p/GngP7x7XPXu9mdcAwAA/yEAosdNSYpRfKRTrnNuvfHBCbPLAQDAcgiA6HH2EJvumDREEmsCAgBgBgIgTOGZDfz6B9WqrWs0uRoAAKyFAAhTjBscqbHxEWpoatb297k1HAAAPYkACNN4RgE3/7VMrEYEAEDPIQDCNHOTh8npCFHJx7XK/4jJIAAA9BQCIExzTUS4Fnyx5Y4vq3YeZBQQAIAeQgCEqb5zy0g5HSEqLj/FKCAAAD2EAAhTMQoIAEDPIwDCdIwCAgDQswiAMB2jgAAA9CwCIAICo4AAAPQcAiACAqOAAAD0HAIgAgajgAAA9AwCIAIGo4AAAPQMAiACCqOAAAD4HwEQAYVRQAAA/I8AiIDDKCAAAP5FAETAuXgU8L8ZBQQAwOcIgAhInlHAkvJT2sUoIAAAPkUAREDiWkAAAPyHAIiAxSggAAD+QQBEwGIUEAAA/yAAIqAxCggAgO8RABHQGAUEAMD3CIAIeIwCAgDgWwRABDxGAQEA8C0CIIICo4AAAPgOARBBgVFAAAB8hwCIoMEoIAAAvkEARNBgFBAAAN8gACKoMAoIAED3EQARVBgFBACg+wiACDoXjwJuKzludjkAAAQdAiCCzjUR4cqcPlKS9PDW91X+aZ3JFQEAEFwCIgCuWbNGiYmJcjqdSklJ0e7duztsn5+fr5SUFDmdTiUlJWndunVt2uTm5mr8+PEKDw/X+PHjtXXr1lZ/z8nJ0Y033qiIiAgNGjRIc+bM0YcffujT44L/3H/rKCVfO0Cn693K2lIsd1Oz2SUBABA0TA+AW7ZsUVZWlpYvX66ioiKlp6dr1qxZKisra7f9kSNHNHv2bKWnp6uoqEgPPfSQFi1apNzcXG+bgoICzZ8/XwsWLFBJSYkWLFigefPm6e233/a2yc/P18KFC/XWW28pLy9PbrdbGRkZ+vzzz/1+zOi+UHuI/u9dkxURHqrCY5/p6dcPmV0SAABBw2aYfBX9lClTlJycrLVr13q3jRs3TnPmzFFOTk6b9kuXLtW2bdtUWlrq3ZaZmamSkhIVFBRIkubPny+Xy6XXXnvN22bmzJmKjo7Wpk2b2q3jxIkTGjRokPLz83XLLbd0qnaXy6WoqCjV1tYqMjKyU8+Bb71aXKHFm4sVYpM2fydNNyUONLskAECA4/vb5BHAhoYGFRYWKiMjo9X2jIwM7d27t93nFBQUtGk/Y8YM7du3T42NjR22udw+Jam2tlaSNHAgASKY3DFpqOYmD1OzIWVtLlJtXaPZJQEAEPBMDYA1NTVqampSXFxcq+1xcXGqqqpq9zlVVVXttne73aqpqemwzeX2aRiGsrOzNW3aNE2YMOGy9dbX18vlcrV6wHw/vuMLGhHTV8drz2nZ1ndZGgYAgCsw/RpASbLZbK1+NwyjzbYrtb90+9Xs8/7779e777572dPDHjk5OYqKivI+EhISOmyPntE/PFT/967JCg2xaft7VdryTrnZJQEAENBMDYCxsbGy2+1tRuaqq6vbjOB5xMfHt9s+NDRUMTExHbZpb5//8R//oW3btumNN97QsGHDOqx32bJlqq2t9T7KywkageKGhAF6YMYYSdKPf3dAh6rPmFwRAACBy9QAGBYWppSUFOXl5bXanpeXp6lTp7b7nLS0tDbtd+zYodTUVDkcjg7bXLxPwzB0//336ze/+Y1ef/11JSYmXrHe8PBwRUZGtnogcHwnPUnTRsXqbGOTFm0qUr27yeySAAAISKafAs7OztZzzz2njRs3qrS0VEuWLFFZWZkyMzMltYy63XPPPd72mZmZOnbsmLKzs1VaWqqNGzdqw4YNeuCBB7xtFi9erB07dmjFihX64IMPtGLFCu3cuVNZWVneNgsXLtRLL72kl19+WREREaqqqlJVVZXOnj3bY8cO3woJsenn827QwH5hOlDp0s/+yLqOAAC0ywgAq1evNoYPH26EhYUZycnJRn5+vvdv9957rzF9+vRW7Xft2mVMnjzZCAsLM0aMGGGsXbu2zT5feeUVY8yYMYbD4TDGjh1r5Obmtvq7pHYfzz//fKfrrq2tNSQZtbW1V3W88K+8/VXG8KW/N4Yv/b3xxgefmF0OACDA8P1tGKavAxjMWEcocD366vt6oeCYYvuH6bXFt+iaiHCzSwIABAi+vwPgFDDgD8tmj9PY+AjVnGnQf75SouZm/p0DAIAHARC9ktNh19P/MlnhoSH6y0cntPHNI2aXBABAwCAAote6Li5CD39lvCRpxR8/0PsVtSZXBABAYCAAolf7xpRrlTE+To1NhhZtLlJdg9vskgAAMB0BEL2azWbTirkTFR/p1OETn+vH2w6YXRIAAKYjAKLXi+4Xpp/Pv0E2m7RlX7n+8G6l2SUBAGAqAiAsYerIWP379JGSpGW/eVcVp1jwGwBgXQRAWMaSf7xOkxIGyHXOre/+cp8++7zB7JIAADAFARCW4bCH6Om7Jiu6r0PvV7j0z/9ToMpaRgIBANZDAISlXBvTV69kpik+0qlD1Wf0f9YW6PCJM2aXBQBAjyIAwnJGDYrQr/89TUmx/VRx6qz+eV0BawQCACyFAAhLGhbdV/+bmaYJQyN18vMG3bX+Lb11+KTZZQEA0CMIgLCs2P7h2nTfFzUlcaDO1Lt1z8a/Ku/AJ2aXBQCA3xEAYWkRTode+LebdNu4ODW4m5X5UqFyCz82uywAAPyKAAjLczrsWveNZM1NHqamZkP/+UqJntt92OyyAADwGwIgICnUHqKn/s9EfWtaoiTpp38o1co/fSjDMEyuDAAA3yMAAueFhNj08O3j9P0ZYyRJz7xxSA//9n01NRMCAQC9CwEQuIjNZtPCW0fp8TsnyGaTfvV2mRZvLlKDu9ns0gAA8BkCINCOu6cM1y/+ZbIcdpt+/26lvvXCO6prcJtdFgAAPkEABC7jKxOH6Ll7b1Qfh127D9bo7ufe1qk67h8MAAh+BECgA9Ovu0a/um+Kovo4VFR2SvP+p0CHqk+bXRYAAN1CAASuIPnaaP3vd9M0KCJcH31yRjNX7dbjfzig0+cazS4NAIAuIQACnTAmPkK/XXizbhsXJ3ezoWd3H9GX/itfW4s+ZqkYAEDQsRl8e3WZy+VSVFSUamtrFRkZaXY56CFvfFitx353QEdqPpck3TgiWj/62hf0hSFRJlcGAOgMvr8JgN3CG8i66t1N2rDniH7x50M629ikEFvLzOH/zLhOA/qGmV0eAKADfH9zChjokvBQu773D6P05/+crq9MHKxmQ/rlW8d068pd2vTXMhaPBgAENEYAu4F/QcBj799r9KNt+/XRJ2ckSROHRenHX/uCJl8bbXJlAIBL8f1NAOwW3kC4WGNTs14sOKZVeR/pdH3LotHzUofpBzPHKrZ/uMnVAQA8+P7mFDDgMw57iL41LVF/fmC65iYPkyT9776PdevKXXr+zSNyN3E7OQBAYGAEsBv4FwQ6UnjsMz267X29X+GSJCUM7KN/mjxM/5Q8VMNj+plcHQBYF9/fBMBu4Q2EK2lqNrT5nTI99acPdaruwsLRqcOjdWfyUH3l+iGK6uswsUIAsB6+vwmA3cIbCJ1V1+DWjv2fKPdvH+vNQzXyTBIOs4fotvGD9E+Th2n6mGvksHNVBgD4G9/fBMBu4Q2ErvjEdU6vFlcot7BCH35y4b7CMf3C9NUbhmhu8jBNGBopm81mYpUA0Hvx/U0A7BbeQOgOwzB0oNKl3/ytQq8WV6jmTIP3b6MH9dc/JQ/TnMlDNDiqj4lVAkDvw/c3AbBbeAPBV9xNzdp9sEa/KarQjv1Vqne3zBi22aS0pBjdPCpWkxIGaOKwKEU4uWYQALqD728CYLfwBoI/uM416rX3KpX7twr99cinrf5ms0nXDYrQ5GsHaFLCAE2+NlqjBvWXPYTTxQDQWXx/EwC7hTcQ/K380zrlHfhEReWnVFT2mT7+7GybNv3DQ3VDQpQmJ0RrUsIATbp2AAtPA0AH+P4mAHYLbyD0tOrT51RcdsobCN/9uFZ1DU1t2l07sK8mXztA4wZHKiG6r4ZF99Gw6D4a2C+MySUALI/v7wAJgGvWrNFTTz2lyspKfeELX9CqVauUnp5+2fb5+fnKzs7W/v37NWTIEP3gBz9QZmZmqza5ubl65JFH9Pe//10jR47U448/rjvvvLNbr3sp3kAwm7upWQerz6iorCUQFpWf0qHqM5dt38dh17DoPhp6PhAO84bDlv/GEBABWADf31Ko2QVs2bJFWVlZWrNmjW6++Wb9z//8j2bNmqUDBw7o2muvbdP+yJEjmj17tu677z699NJLevPNN/W9731P11xzjebOnStJKigo0Pz58/WTn/xEd955p7Zu3ap58+Zpz549mjJlSpdeFwhEofYQjRscqXGDI/X1KS3v29qzjXr341MqKjulv584o48/O6uPP6vTJ656nW1s0sHqMzp4mZDodIRoWHRfDR3QMloY1cehAX0d3v8O6BOmSO/PDkX2cbB2IQAEIdNHAKdMmaLk5GStXbvWu23cuHGaM2eOcnJy2rRfunSptm3bptLSUu+2zMxMlZSUqKCgQJI0f/58uVwuvfbaa942M2fOVHR0tDZt2tSl120P/4JAMKl3N+n4qXP6+LM6byis+Ozs+Z/P6pPT59SVT4P+4aGK6nMhJPYPD1WfMLucoXY5HSFyen+2q48jRE6H/fzj4p9bfnfYQ+QICVGo3dbyOP+zd1uIjRFKAN3G97fJI4ANDQ0qLCzUgw8+2Gp7RkaG9u7d2+5zCgoKlJGR0WrbjBkztGHDBjU2NsrhcKigoEBLlixp02bVqlVdfl1Jqq+vV319vfd3l8t1xWMEAkV4qF2Jsf2UGNv+fYjr3U2qPHVOH392VhWn6nSqrlGnzjaq9myjauta/nvqbINOnf/59Dm3JOlMvVtn6t2qONV2goo/2ENagqDDHnIhJIbYFGKTQkJsCrG187Pt/M8hkt3WEiJDbC37ssmm8/+T7Xxbm02yqeW/0sXbJJvN5m3bskXedp5oeuH3S/5+yfZWOrepwwDcUTQ2IzcT1eELMyfEa+aEwWaX0euYGgBramrU1NSkuLi4Vtvj4uJUVVXV7nOqqqrabe92u1VTU6PBgwdfto1nn115XUnKycnRj3/8404fHxBMwkPtGhHbTyMuExAv5W5q1ulzbm9IPFXX4A2G5xqbVO9u1rnGJp1taNI5d5PONjTrnLtJ9Y1NOtvYpHONzd6/1Tc262xjkxqbmuVuMuRublZjU/vDkU3NhpqaDe9aiQB6txGx/QiAfmD6NYBS23/RGobR8b9y22l/6fbO7PNqX3fZsmXKzs72/u5yuZSQkHDZ9kBvFmoPUXS/MEX3C/PbazQ1Gy2hsNlQU5OhxuYLAfHioNjUbKjZMNRstDzHMDzbWv5/3XT+b82XtGs2DBmGZMjz3wufJ4ahi/7est2QpIvan//V277ld+OS31s3uDTWXnra/eKrcq7UtrO6c51PAMwThMVNvjba7BJ6JVMDYGxsrOx2e5tRt+rq6jajcx7x8fHttg8NDVVMTEyHbTz77MrrSlJ4eLjCw1lfDegp9hCb7CF2s8sAgF7H1Ol7YWFhSklJUV5eXqvteXl5mjp1arvPSUtLa9N+x44dSk1NlcPh6LCNZ59deV0AAIDewvRTwNnZ2VqwYIFSU1OVlpam9evXq6yszLuu37Jly1RRUaEXX3xRUsuM32eeeUbZ2dm67777VFBQoA0bNnhn90rS4sWLdcstt2jFihW644479Oqrr2rnzp3as2dPp18XAACgtzI9AM6fP18nT57UY489psrKSk2YMEHbt2/X8OHDJUmVlZUqKyvztk9MTNT27du1ZMkSrV69WkOGDNHTTz/tXQNQkqZOnarNmzfr4Ycf1iOPPKKRI0dqy5Yt3jUAO/O6AAAAvZXp6wAGM9YRAgAg+PD9bfI1gAAAAOh5BEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxpt8KLph5bqLicrlMrgQAAHSW53vbyjdDIwB2w+nTpyVJCQkJJlcCAACu1unTpxUVFWV2GabgXsDd0NzcrOPHjysiIkI2m+2qnutyuZSQkKDy8nLL3ofQg75oQT9cQF+0oB8uoC9a0A8XdKcvDMPQ6dOnNWTIEIWEWPNqOEYAuyEkJETDhg3r1j4iIyMt/39iD/qiBf1wAX3Rgn64gL5oQT9c0NW+sOrIn4c1Yy8AAICFEQABAAAshgBokvDwcD366KMKDw83uxTT0Rct6IcL6IsW9MMF9EUL+uEC+qJ7mAQCAABgMYwAAgAAWAwBEAAAwGIIgAAAABZDAAQAALAYAqAPrVmzRomJiXI6nUpJSdHu3bs7bJ+fn6+UlBQ5nU4lJSVp3bp1bdrk5uZq/PjxCg8P1/jx47V161Z/le8zvu6HZ599Vunp6YqOjlZ0dLRuu+02/fWvf/XnIfiMP94THps3b5bNZtOcOXN8XLXv+aMfTp06pYULF2rw4MFyOp0aN26ctm/f7q9D8Bl/9MWqVas0ZswY9enTRwkJCVqyZInOnTvnr0Pwiavph8rKSn3961/XmDFjFBISoqysrHbbBePnpeT7vgjWz0x/vCc8gunzsscY8InNmzcbDofDePbZZ40DBw4YixcvNvr162ccO3as3faHDx82+vbtayxevNg4cOCA8eyzzxoOh8P49a9/7W2zd+9ew263G0888YRRWlpqPPHEE0ZoaKjx1ltv9dRhXTV/9MPXv/51Y/Xq1UZRUZFRWlpq/Ou//qsRFRVlfPzxxz11WF3ij77wOHr0qDF06FAjPT3duOOOO/x8JN3jj36or683UlNTjdmzZxt79uwxjh49auzevdsoLi7uqcPqEn/0xUsvvWSEh4cbv/rVr4wjR44Yf/rTn4zBgwcbWVlZPXVYV+1q++HIkSPGokWLjBdeeMGYNGmSsXjx4jZtgvHz0jD80xfB+Jnpj37wCKbPy55EAPSRm266ycjMzGy1bezYscaDDz7Ybvsf/OAHxtixY1tt++53v2t88Ytf9P4+b948Y+bMma3azJgxw7jrrrt8VLXv+aMfLuV2u42IiAjjhRde6H7BfuSvvnC73cbNN99sPPfcc8a9994b8B9o/uiHtWvXGklJSUZDQ4PvC/Yjf/TFwoULjS996Uut2mRnZxvTpk3zUdW+d7X9cLHp06e3+2UfjJ+XhuGfvrhUMHxm+qsfgu3zsidxCtgHGhoaVFhYqIyMjFbbMzIytHfv3nafU1BQ0Kb9jBkztG/fPjU2NnbY5nL7NJu/+uFSdXV1amxs1MCBA31TuB/4sy8ee+wxXXPNNfrWt77l+8J9zF/9sG3bNqWlpWnhwoWKi4vThAkT9MQTT6ipqck/B+ID/uqLadOmqbCw0HuK7/Dhw9q+fbtuv/12PxxF93WlHzoj2D4vJf/1xaUC/TPTn/0QTJ+XPS3U7AJ6g5qaGjU1NSkuLq7V9ri4OFVVVbX7nKqqqnbbu91u1dTUaPDgwZdtc7l9ms1f/XCpBx98UEOHDtVtt93mu+J9zF998eabb2rDhg0qLi72V+k+5a9+OHz4sF5//XXdfffd2r59uw4ePKiFCxfK7Xbrhz/8od+Opzv81Rd33XWXTpw4oWnTpskwDLndbv37v/+7HnzwQb8dS3d0pR86I9g+LyX/9cWlAv0z01/9EGyflz2NAOhDNput1e+GYbTZdqX2l26/2n0GAn/0g8fPfvYzbdq0Sbt27ZLT6fRBtf7ly744ffq0vvGNb+jZZ59VbGys74v1I1+/J5qbmzVo0CCtX79edrtdKSkpOn78uJ566qmADYAevu6LXbt26fHHH9eaNWs0ZcoUHTp0SIsXL9bgwYP1yCOP+Lh63/HHZ1swfl5K/q07mD4zfdkPwfx52VMIgD4QGxsru93e5l8q1dXVbf5F4xEfH99u+9DQUMXExHTY5nL7NJu/+sFj5cqVeuKJJ7Rz505NnDjRt8X7mD/6Yv/+/Tp69Ki++tWvev/e3NwsSQoNDdWHH36okSNH+vhIusdf74nBgwfL4XDIbrd724wbN05VVVVqaGhQWFiYj4+k+/zVF4888ogWLFigb3/725Kk66+/Xp9//rm+853vaPny5QoJCawrfbrSD50RbJ+Xkv/6wiNYPjP90Q9///vfg+7zsqcF1idDkAoLC1NKSory8vJabc/Ly9PUqVPbfU5aWlqb9jt27FBqaqocDkeHbS63T7P5qx8k6amnntJPfvIT/fGPf1Rqaqrvi/cxf/TF2LFj9d5776m4uNj7+NrXvqZbb71VxcXFSkhI8NvxdJW/3hM333yzDh065P1Al6SPPvpIgwcPDsjwJ/mvL+rq6tqEPLvdLqNlkp8Pj8A3utIPnRFsn5eS//pCCq7PTH/0QzB+Xva4np510lt5prBv2LDBOHDggJGVlWX069fPOHr0qGEYhvHggw8aCxYs8Lb3LO+wZMkS48CBA8aGDRvaLO/w5ptvGna73XjyySeN0tJS48knnwz4ZQ380Q8rVqwwwsLCjF//+tdGZWWl93H69OkeP76r4Y++uFQwzGrzRz+UlZUZ/fv3N+6//37jww8/NH7/+98bgwYNMn7605/2+PFdDX/0xaOPPmpEREQYmzZtMg4fPmzs2LHDGDlypDFv3rweP77Outp+MAzDKCoqMoqKioyUlBTj61//ulFUVGTs37/f+/dg/Lw0DP/0RTB+ZvqjHy4VDJ+XPYkA6EOrV682hg8fboSFhRnJyclGfn6+92/33nuvMX369Fbtd+3aZUyePNkICwszRowYYaxdu7bNPl955RVjzJgxhsPhMMaOHWvk5ub6+zC6zdf9MHz4cENSm8ejjz7aA0fTPf54T1wsWD7Q/NEPe/fuNaZMmWKEh4cbSUlJxuOPP2643W5/H0q3+bovGhsbjR/96EfGyJEjDafTaSQkJBjf+973jM8++6wHjqbrrrYf2vsMGD58eKs2wfh5aRi+74tg/cz0x3viYsHyedlTbIYRgOcIAAAA4DdcAwgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACwHknTpxQfHy8nnjiCe+2t99+W2FhYdqxY4eJlQGAb3EvYAC4yPbt2zVnzhzt3btXY8eO1eTJk3X77bdr1apVZpcGAD5DAASASyxcuFA7d+7UjTfeqJKSEr3zzjtyOp1mlwUAPkMABIBLnD17VhMmTFB5ebn27duniRMnml0SAPgU1wACwCUOHz6s48ePq7m5WceOHTO7HADwOUYAAeAiDQ0NuummmzRp0iSNHTtWP//5z/Xee+8pLi7O7NIAwGcIgABwke9///v69a9/rZKSEvXv31+33nqrIiIi9Pvf/97s0gDAZzgFDADn7dq1S6tWrdIvf/lLRUZGKiQkRL/85S+1Z88erV271uzyAMBnGAEEAACwGEYAAQAALIYACAAAYDEEQAAAAIshAAIAAFgMARAAAMBiCIAAAAAWQwAEAACwGAIgAACAxRAAAQAALIYACAAAYDEEQAAAAIshAAIAAFjM/wcAyvXLwKwKLQAAAABJRU5ErkJggg==", "text/html": [ "\n", "