{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# IP3 demo" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:55.849353Z", "iopub.status.busy": "2025-08-18T03:36:55.849210Z", "iopub.status.idle": "2025-08-18T03:36:56.263997Z", "shell.execute_reply": "2025-08-18T03:36:56.263636Z" } }, "outputs": [ { "data": { "text/plain": [ "1.0" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from neuron import n, rxd\n", "\n", "n.load_file(\"stdrun.hoc\")" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:56.356579Z", "iopub.status.busy": "2025-08-18T03:36:56.354334Z", "iopub.status.idle": "2025-08-18T03:36:56.361568Z", "shell.execute_reply": "2025-08-18T03:36:56.359307Z" } }, "outputs": [], "source": [ "dend = n.Section(name=\"dend\")\n", "dend.L = 100\n", "dend.nseg = 101" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:56.363212Z", "iopub.status.busy": "2025-08-18T03:36:56.363070Z", "iopub.status.idle": "2025-08-18T03:36:56.368291Z", "shell.execute_reply": "2025-08-18T03:36:56.367926Z" } }, "outputs": [], "source": [ "def my_initial(node):\n", " if 0.4 < node.segment.x < 0.6:\n", " return 1\n", " else:\n", " return 0\n", "\n", "\n", "cyt = rxd.Region(n.allsec(), name=\"cyt\", nrn_region=\"i\")\n", "ip3 = rxd.Species(cyt, name=\"ip3\", d=1, initial=my_initial)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:56.370333Z", "iopub.status.busy": "2025-08-18T03:36:56.369757Z", "iopub.status.idle": "2025-08-18T03:36:57.038640Z", "shell.execute_reply": "2025-08-18T03:36:57.038197Z" } }, "outputs": [], "source": [ "from matplotlib import pyplot\n", "\n", "\n", "def plot_it():\n", " xs = [seg.x * dend.L for seg in dend]\n", " ys = [seg.ip3i for seg in dend]\n", " pyplot.plot(xs, ys)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:57.040719Z", "iopub.status.busy": "2025-08-18T03:36:57.040496Z", "iopub.status.idle": "2025-08-18T03:36:57.354730Z", "shell.execute_reply": "2025-08-18T03:36:57.354326Z" } }, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, '[IP3]')" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAir1JREFUeJzs3Xd4FNUax/HvtlRSCUmABEIv0ntRAUEBFTsKKiAKCoqCWAARUFGxwcUCojQbCIiCiIqFKtIDofcACSUJ6b3tzv1jyJKYACmzuynv53n2ZrM75WSubH45855zdIqiKAghhBBCVBJ6RzdACCGEEEJLEm6EEEIIUalIuBFCCCFEpSLhRgghhBCVioQbIYQQQlQqEm6EEEIIUalIuBFCCCFEpWJ0dAPszWKxcPHiRTw8PNDpdI5ujhBCCCGKQVEUUlJSqFWrFnr99ftmqly4uXjxIsHBwY5uhhBCCCFKITIykqCgoOtuU+XCjYeHB6BeHE9PTwe3RgghhBDFkZycTHBwsPX3+PVUuXCTdyvK09NTwo0QQghRwRSnpEQKioUQQghRqUi4EUIIIUSlIuFGCCGEEJWKhBshhBBCVCoSboQQQghRqUi4EUIIIUSlIuFGCCGEEJWKhBshhBBCVCoSboQQQghRqUi4EUIIIUSl4tBws2XLFgYMGECtWrXQ6XSsXr36hvts2rSJdu3a4ezsTMOGDfnqq69s3k4hhBBCVBwODTdpaWm0bt2aOXPmFGv7M2fOcNddd9GrVy/CwsIYN24cI0aM4I8//rBxS4UQQghRUTh04cz+/fvTv3//Ym8/b9486tWrx8yZMwFo1qwZW7du5X//+x99+/a1VTOFqJIUReFiUiaKoji6KZWeh4sJL1eTo5shRKVRoVYF3759O3369CnwWt++fRk3btw198nKyiIrK8v6fXJysq2aJ0SlMnZZGGv2X3R0M6oEo17H8me60L6ur6ObIkSlUKHCTVRUFAEBAQVeCwgIIDk5mYyMDFxdXQvtM2PGDN588017NVGISmNfZAIATgY9Op2DG1OJ5Zgt5FoUDl9MlnAjhEYqVLgpjUmTJjF+/Hjr98nJyQQHBzuwRUJUDBaL+vWHUV1pHezt0LZUZmOW7mXtgUuYLXL7TwitVKhwExgYSHR0dIHXoqOj8fT0LLLXBsDZ2RlnZ2d7NE+ISiX3Srox6KXbxpbyrq+EGyG0U6HmuenatSvr168v8Npff/1F165dHdQiISov85WeGwk3tiXhRgjtOTTcpKamEhYWRlhYGKAO9Q4LCyMiIgJQbykNHTrUuv2oUaMIDw/n1Vdf5dixY8ydO5cVK1bw4osvOqL5QlRqZum5sQvDlYKmXAk3QmjGoeFmz549tG3blrZt2wIwfvx42rZty9SpUwG4dOmSNegA1KtXj19//ZW//vqL1q1bM3PmTBYsWCDDwIWwgbyeBAk3tmU0qNfXIuFGCM04tOamZ8+e151Do6jZh3v27Mm+ffts2CohBEDe71qDDJWyKf2V62uW+YSE0EyFqrkRQtiPFBTbh1FqboTQnIQbIUSRLFJQbBd6CTdCaE7CjRCiSHk9N0YJNzYlPTdCaE/CjRCiEEVRrDU3egk3NiU9N0JoT8KNEKKQ/L9opaDYtmQouBDak3AjhCgk/8gdg0HCjS3l3ZayyGgpITQj4UYIUUheMTFIz42tyW0pIbQn4UYIUUhuvnQjo6VsSwqKhdCehBshRCEFem4k3NiU9NwIoT0JN0KIQgr03MhtKZuSnhshtCfhRghRSF5BsU4nQ8FtTZZfEEJ7Em6EEIVYF82UXhuby7vtJ0PBhdCOhBshRCGyIrj9WIeCS7gRQjMSboQQhUi4sR+99NwIoTkJN0KIQiTc2I/03AihPQk3QohC8mbLlXBje1JQLIT2JNwIIQrJu0UiK4LbntEgQ8GF0JqEGyFEIXm/aPUyWsrmrD03Em6E0IyEGyFEIWbpubEbo179GJaCYiG0I+FGCFGItedGwo3NGa58CktBsRDakXAjhChERkvZT95tKem5EUI7Em6EEIVIuLGfvIJii4yWEkIzEm6EEIXkDUuW5RdsTwqKhdCehBshRCHSc2M/eQXFEm6E0I6EGyFEIRJu7OdKtpFwI4SGJNwIIQqRoeD2Iz03QmhPwo0QohAZCm4/eUPBZfkFIbQj4UYIUYj1tpQUFNucdSi4WcKNEFqRcCOEKMQsC2faTd5tKRkKLoR2JNwIIQqRgmL7ySsolkn8hNCOhBshRCESbuzH2nMj4UYIzUi4EUIUIuHGfqSgWAjtSbgRQhQiQ8Htx5A3FFwKioXQjIQbIUQheb0IehktZXN5I9Kk50YI7Ui4EUIUYu25MUi4sTWDQVYFF0JrEm6EEIVYJ/GTnhuby+u5kYJiIbQj4UYIUYgUFNuPDAUXQnsSboQQhUi4sZ+8oeAgvTdCaEXCjRCiEOsMxXJbyubyX2MpKhZCGxJuhBCF5A1LloJi2zPku8ayMrgQ2pBwI4QoRIaC20+BnhsJN0JoQsKNEKIQmcTPfvLXNUlRsRDakHAjhCjEOhRcwo3N5Q83UlAshDYk3AghCrGOlpLbUjaXPz9Kz40Q2pBwI4QoxBpupKDY5nQ6nbX3xiKjpYTQhIQbIUQhMhTcvqzrS0nPjRCakHAjhChECortK6/nRsKNENqQcCOEKEQKiu1Lwo0Q2pJwI4QoRHpu7Csv3EhBsRDakHAjhChEem7sSwqKhdCWhBshRCEyFNy+8maCzjVLuBFCCxJuhBCFWEdLSc+NXRil50YITUm4EUIUklf7IeHGPqTmRghtSbgRQhRikYJiu5LRUkJoS8KNEKIQKSi2LykoFkJbEm6EEIXIUHD7st6WkoJiITQh4UYIUUheQbFeRkvZRd6oNOm5EUIbDg83c+bMISQkBBcXFzp37syuXbuuu/3s2bNp0qQJrq6uBAcH8+KLL5KZmWmn1gpRNVh7bmThTLuQgmIhtOXQcLN8+XLGjx/PtGnT2Lt3L61bt6Zv377ExMQUuf3SpUuZOHEi06ZN4+jRoyxcuJDly5fz2muv2bnlQlRu1pob6bmxC2vNjYQbITTh0HAza9YsRo4cyfDhw2nevDnz5s3Dzc2NRYsWFbn9tm3b6N69O48++ighISHccccdDB48+Lq9PVlZWSQnJxd4CCGuT4aC25deem6E0JTDwk12djahoaH06dPnamP0evr06cP27duL3Kdbt26EhoZaw0x4eDi//fYbd9555zXPM2PGDLy8vKyP4OBgbX8QISohGQpuX0YZCi6EpoyOOnFsbCxms5mAgIACrwcEBHDs2LEi93n00UeJjY3l5ptvRlEUcnNzGTVq1HVvS02aNInx48dbv09OTpaAI8QNSEGxfUlBsRDacnhBcUls2rSJd999l7lz57J3715++uknfv31V6ZPn37NfZydnfH09CzwEEJcnxQU25cUFAuhLYf13Pj5+WEwGIiOji7wenR0NIGBgUXuM2XKFIYMGcKIESMAaNmyJWlpaTz99NNMnjwZvb5CZTUhyi0pKLYvKSgWQlsOSwNOTk60b9+e9evXW1+zWCysX7+erl27FrlPenp6oQBjMBgAUKQ7VwjNXJ3ET/5gsAfpuRFCWw7ruQEYP348w4YNo0OHDnTq1InZs2eTlpbG8OHDARg6dCi1a9dmxowZAAwYMIBZs2bRtm1bOnfuzKlTp5gyZQoDBgywhhwhRNldXX7BwQ2pIqTnRghtOTTcPPLII1y+fJmpU6cSFRVFmzZtWLdunbXIOCIiokBPzeuvv45Op+P111/nwoUL1KhRgwEDBvDOO+846kcQolLKCzcGuS1lF3m3/6TnRghtODTcAIwZM4YxY8YU+d6mTZsKfG80Gpk2bRrTpk2zQ8uEqLryRktJQbF9WIeCy+11ITQhnc5CiELyFnCUgmL7yLstZTZbHNwSISoHCTdCiELy5luRgmL7sIYb6bgRQhPyySWEKEQKiu1LCoqF0JZ8dAkhCpGh4PYlQ8GF0JZ8cgkhCskrbDXIJ4RdyPILQmhLPrqEEIWYzXnhRj4i7MFwZVRarhTdCKEJ+eQSQhRi7bmR0VJ2kXedZSi4ENqQcCOEKCRXCortyjpayiJDwYXQgnx0CSEKsUhBsV1dDTcObogQlYR8cgkhCsm7PSLZxj6sQ8HltpQQmpCPLiFEARaLQt7vWOm5sQ/rUHApKBZCE/LJJYQoIH9RqxQU24cMBRdCWxJuhBAFmPNNJGeQhTPt4uokflJ0I4QWJNwIIQooEG6k58YupKBYCG1JuBFCFJB/CQApubEPGQouhLbko0sIUUD+xRuloNg+pOdGCG3JJ5cQooD8BcV6uStlF1JQLIS2JNwIIQrIq7kx6HXopObGLmRVcCG0JeFGCFGANdxIsLEb6yR+Em6E0ISEGyFEAfl7boR9yFBwIbQl4UYIUYCEG/uTgmIhtCXhRghRgHVFcMk2dpN3C1CGgguhDQk3QogC8kbsGA3y8WAv1p4bKbkRQhPy6SWEKCBv8UZ9RS0oroC9HzKJnxDaMjq6AUKI8sXac1MR7kvFHIULoRB9BKIPQcwRyEiA6g3BvzkE3KQ+6nYDFy9Ht/aaroYb6boRQgsSboQQBVSIguJLB2DDdDj5Z9HvXz6mPg7/pH7v4gXdx0HnUeDkZrdmFtfVoeAObogQlYSEGyFEAbnlOdzEnYaN78ChH9XvdQa1VyagxZVemubg5gexJ9SenOgjcH4XJJyF9W/Cznlw6yvQbhgYnRz6o+QnQ8GF0JaEGyFEAXm3pcpVuLFY1FDz72yw5KqvtXgQer4Gfg0Lb+9TFxrdfmVfMxz8Qd0/MQJ+exm2z4GBi6FWW7v9CNdjHS0ld6WE0IQUFAshCsgrKC434SYnA1YOh38+UoNNozvgmX/goUVFB5v/0hug9SAYEwp3fgTu/pBwBhbfCcfX2b79xWAwSEGxEFqScCOEKMDac1MeRkulxcE398KR1aA3wf1fwGM/QM1WJT+W0Qk6jYTn90D9XpCTDssGw675mje7pK7Oc+PghghRSUi4EUIUYJ3Ez9E9N3GnYWEfiNypFgQPWaX2wJSVi5cakNo+DopFvU31x2SHVvPKUHAhtCXhRghRQN7ijQ4dCh5zFBb0gfhw8K4DT/0F9W7R7vgGE9zzGdz2uvr99s9g9ShQHFP0IkPBhdCWhBshRAFmR/fcpMXC0kcgI14t+B2xHmo00f48Op06cuqB+aA3woHlsOk97c9TDNah4JJthNCEhBshRAG5juy5yc2CZY9B4jnwqQeP/wTV/G17zlYPw92z1eeb34ODK217viLIUHAhtCXhRghRgMMKihUFfhkLkTvA2QseXQFuvvY5d7sh0O159fnqZ+H8Hvuc94q8ay3ZRghtSLgRQhTgsEn8/p0N+79XJ+Z7+Cuo0di+5+/zJjS5E8xZ8P1gSIy026ml50YIbUm4EUIUYHFEuDm6Fv5+U33e/31ocJv9zp1Hb1DrbwJaQFoMfD8IslLtcuqrBcV2OZ0QlZ6EGyFEAXYfCp58EX5+FlCg09PqXDSO4lwNBi9TJ/qLPgR/TbXLaWUouBDaknAjhCjArkPB8+psMpPUkVF9Z9j+nDfiHQwPLlCf71kI4ZtsfkoZCi6EtiTcCCEKsPbc2KOgOGypurK3wQnumweGcrLcXf0e0HGE+vzn5yErxaanuzpDsYQbIbQg4UYIUYBZsVPPTdIFWDdJfd7rNfBvatvzlVSfN8G7LiRFwJ9TbHoqa8+NgyYRFKKykXAjhCjALgXFigK/vABZSVC7PXR93nbnKi3nanDvHPV56GI4vcFmp7JO4iclN0JoQsKNEKIAuwwF3/cdnPobDM5w3+fl53bUf9W7RS1yBvX2VGayTU5jlKHgQmhKwo0QogCb99wkX4Q/XlOf93rNNksraKnPG+ATAsnn4S/b3J7S51t+QZFbU0KUmYQbIUQBNu+5+ftNyEpWb0d1K4e3o/7LyR3unas+D/0aLoZpfor89U1SVCxE2Um4EUIUYNPlFy6EwoFl6vM7P1QnzqsIQrpDy4GAAn9M1nz18PxzCklRsRBlJ+FGCFFArtlGk/gpV4IBQKtBas9NRdJ7Ghhd4NxWOLZW00PnD5LScyNE2Um4EUIUYLOh4Ed+hojtYHSF3vaZ+VdT3sHQdYz6/M8p6grmGjHIbSkhNCXhRghRgE0KinMyry5l0P0F8Kqt3bHt6eYXoVoAJJyBXfM1O2z+ay0DpoQoOwk3QogCbFJQvOsLSDwHHjWh+1jtjmtvztXgtisjpjZ/AGlxmhw2/20pGQ4uRNlJuBFCFGAtKNYq3KRehi0fqc97T1VHH1VkbR6FwJbqBISbtFkLS6/XkZdvpKBYiLKTcCOEKCCvoFizcLNphjr0u2ZrtZC4otMb4I531Od7FsHl45oc1iiLZwqhGQk3QogCNB0KnhgBe79Wn/d9F/SV5COnfg9o3B8UM2x+X5ND6mXxTCE0U0k+aYQQWsmr+dBkKPg/M8GSC/V6QMjNZT9eedLryizLh37SpPfGID03QmhGwo0QogDzlXrWMg8FT4yEfUvU5z0nlu1Y5VHNVtD0bkBRi4vLSMKNENqRcCOEKMB8peemzDU3W2eBJQfq3Qp1u2nQsnKox6vq10M/wuUTZTqUhBshtOPwcDNnzhxCQkJwcXGhc+fO7Nq167rbJyYm8txzz1GzZk2cnZ1p3Lgxv/32m51aK0Tll9dzU6Zwk3Qe9n6rPu9RCXtt8tRsDU3uBBTY8mGZDmUtKJbRUkKUmUPDzfLlyxk/fjzTpk1j7969tG7dmr59+xITE1Pk9tnZ2dx+++2cPXuWlStXcvz4cebPn0/t2hV0QjAhyiFNCor/udJrE3KLui5TZWbtvVkJsadKfRgpKBZCOw4NN7NmzWLkyJEMHz6c5s2bM2/ePNzc3Fi0aFGR2y9atIj4+HhWr15N9+7dCQkJoUePHrRu3fqa58jKyiI5ObnAQwhxbWWexC/pAuzL67WZoFGryrFaba+MnLKUqfdGhoILoR2HhZvs7GxCQ0Pp06fP1cbo9fTp04ft27cXuc+aNWvo2rUrzz33HAEBAbRo0YJ3330Xs9l8zfPMmDEDLy8v6yM4OFjzn0WIyqTMyy9s/R+Ys6HuzVDvFg1bVo71vBLiDq6AuNOlOoRewo0QmnFYuImNjcVsNhMQEFDg9YCAAKKioorcJzw8nJUrV2I2m/ntt9+YMmUKM2fO5O23377meSZNmkRSUpL1ERkZqenPIURlU6ah4MkXr85r07MK9NrkqdUWGvUtU++NFBQLoR2HFxSXhMViwd/fny+//JL27dvzyCOPMHnyZObNm3fNfZydnfH09CzwEEJcW5mGgu+Yq/ba1Omq1ttUJXlh7sAKtaC6hCTcCKEdh4UbPz8/DAYD0dHRBV6Pjo4mMDCwyH1q1qxJ48aNMRgM1teaNWtGVFQU2dnZNm2vEFWFdSh4SQuKM5Mh9EqvTfdxoMUMxxVJ7fZqoFPMsPPaf3Bdi0EKioXQjMPCjZOTE+3bt2f9+vXW1ywWC+vXr6dr165F7tO9e3dOnTqFJd+quSdOnKBmzZo4OTnZvM1CVAVXlpYqec3Nvm/VNaSqN4JGd2jfsIqg6xj1a+jXatgrAYMMBRdCMw69LTV+/Hjmz5/P119/zdGjRxk9ejRpaWkMHz4cgKFDhzJp0iTr9qNHjyY+Pp6xY8dy4sQJfv31V959912ee+45R/0IQlQ6pSooNufCjs/V593GVJ41pEqq0R3g11gNeXkjxopJbksJoR2jI0/+yCOPcPnyZaZOnUpUVBRt2rRh3bp11iLjiIgI9Pk+JIODg/njjz948cUXadWqFbVr12bs2LFMmFCFCheFsLHc0sxQfGQ1JEWCm1/lWPm7tPR66Poc/DJWDXudngFD8T5mZSi4ENpxaLgBGDNmDGPGjCnyvU2bNhV6rWvXruzYscPGrRKi6rKUdIZiRYHtn6nPOz0NJhfbNKyiaDUINrythr0jq6HlQ8XaTYaCC6GdKtp3LIS4lhL33JzbBhf3gdEFOj5lw5ZVECYX6DhSfb7tUzX8FYP03AihHQk3QogCrAXFxR3ttO1T9WvrweDuZ5tGVTQdn1LD3qUwOPdvsXaxLr8gBcVClJmEGyFEASVaFTz2JJz4HdCptSZC5e4HbR5Vn2/7rFi7SEGxENqRcCOEKKBEq4Jvn6N+bdIf/BrZrlEVUZfnAJ0a/mJP3nBzCTdCaEfCjRCigGIPBc9IgP3L1Oddix4UUKX5NVRDH8CuL2+4uYQbIbQj4UYIUUCxC4r3LYHcDAhoAXW72aFlFVCnp9WvYd9DVsp1N5WCYiG0I+FGCFGApTgzFFsssHuB+rzjiKq31EJx1e+pzticnXK1l+sapKBYCO1IuBFCFFCsnpvTGyDhDDh7QauH7dSyCkinU8MfqGHwOsHFaJCeGyG0IuFGCFGAdRK/6/XG7J6vfm3zKDi5275RFVmbwWByh8vH4OzWa26ml4UzhdCMhBshRAE37LlJOAsn/lCf5/VKiGtzyde7lRcKiyAFxUJoR8KNEKKAGw4F37MIUKDBbeqIIHFjna7MWHx0LSRfLHITCTdCaEfCjRCigOtO4peTAXu/UZ/nLTEgbizgJqjbHRQz7Flc5CZ5twFzJdwIUWYOXzhTCFG+mK83z82hn9T5bbzqQOO+dm7ZVek56YRdDiMsJoyotCgSMhOIz4onITOBrNwsvFy88HX2xdvFG18XXxp6N6RDQAfqedVD56iRXR1HqEsxhH4Ft74CRqcCb+cVFFsk3AhRZhJuhBAFWK63tlRezUiH4aA32K1NiqKw//J+NkRuIDQ6lCOxR8hVcq+5fUxGTJGv+7r40j6gPZ0CO9EvpB/eLt42anERmg2AaoGQGgVH1xRaLVyGgguhHQk3QogCrllQfCFUXf3b4ATthtqlLanZqawNX8sPJ37gRMKJAu/VdK9J+4D21POqh4+LD77Ovvi4+OBsdCYpM4n4rHgSMxO5nHGZQ7GH2H95P/GZ8fx17i/+OvcXH+7+kDtC7uDhJg/TpkYb2/foGEzQ/gnY/J46LPw/4UYm8RNCO8UON76+viU6sE6nY+/evdStW7fEjRJCOI7lWgXFebUiN91v89W/o9Ki+PLAl6wNX0tGbgYAzgZn+tTtQ7da3Wgf0J7a1WqX6JjZ5mwOxx1mT9Qe/jz3J8fij7E2fC1rw9fS0LshQ5sP5Z4G92CwZY9U+ydgy4cQsR1ijoF/U+tbegk3Qmim2OEmMTGR2bNn4+XldcNtFUXh2WefxWw2l6lxQgj7y+u5MeYPN5nJcOhH9Xn74TY7d1pOGosOLeLrw1+TZc4CoJ5XPR5u/DADGgzAy/nGnz/X4mRwoq1/W9r6t2VEyxEcij3EihMrWHdmHacSTzF121SWHF3CSx1eomutrlr9SAV51lTXmzq2FvZ+Df1mWN+SnhshtFOi21KDBg3C39+/WNs+//zzpWqQEMJxFEWx1tzo84ebQyshJx38mkCdLpqf12wxs+rUKj7b9xlxmXEAtPNvx5i2Y+gQ0EHzW0Y6nY6WNVrSskZLXun4CitPrGTBgQUcTzjO0389za1Bt/JS+5eo711f0/MC0G6YGm72fw+9p4HJBZCeGyG0VOxwY8nrqy6mlJTrLxInhCh/8v9iLVBQHPqV+rX9MM3XkTqXfI6JWyZyKO4QAHU86jC+/Xhuq3ObXUY2eTp58mSLJ3mg4QPMOzCP5ceWs+X8Fv698C/PtH6GkS1HYtRrWJ7YsDd4BkHyeTj6C7QaCMhQcCG0JPPcCCGs8o/UMVwZmszFfXBpv1pI3GqQZudSFIUfT/zIwF8GcijuEB5OHkzoOIHV966md93edh+y7e3izcROE1l17yp6BvfErJiZGzaX4euGE5kSqd2J9AZoN0R9nhcauXpbyiKjpYQosxKHm7i4ODZu3Eh8fDwAsbGxvP/++7z11lscPXpU8wYKIewnfwettecm9Gv1a7N7wL26JudJyExg3MZxvLH9DTJyM+gc2Jmf7vmJx5s/jslg0uQcpRXiFcKnt33KjFtmUM1UjbDLYQz8ZSBrTq9B0Sp4tH0cdHo4txViTwFyW0oILZUo3OzatYsGDRrQu3dvGjZsSGhoKJ06dWLhwoV88803tG/fnr1799qqrUIIG8vNl24Meh1kpcLBH9QX2j+hyTn2X97Pg2seZEPkBox6Iy+1f4kv7/iSQPdATY6vlbvr383Ke1bSzr8daTlpTN46mde2vmYtdC4TryBodIf6fO9XgBQUC6GlEoWbyZMnM3DgQJKSknjttde477776N27NydOnODUqVMMGjSI6dOn26qtQggbK9Bzo9epI6SyU8G3AYTcXObj/37md55c9ySXMy5T36s+S+9cyhMtnkCvK593yGtXq82ivot4vu3zGHQG1oav5ak/niIuI67sB283TP0athRys6TnRggNlegTJTQ0lPHjx+Ph4cHYsWO5ePEiI0deXV9mzJgx7N69W/NGCiHso0DPjU6nDleGMhcSK4rC5/s/59Utr5JtyaZnUE++v+t7mlVvVtYm25xBb+DpVk/zxe1f4OHkwf7L+3nst8c4lXCqbAdudAd41IT0ODj2q/TcCKGhEoWb7OxsXF1dATCZTLi5ueHnd3UyLz8/P+LiNPiLRgjhEHkFxTod6GMOqbMS603Q+tFSHzPLnMWkrZOYGzYXgKHNhzK712zcTG6atNleOtfszJI7lxDsEcyF1AsM+X0I/174t/QHNBih7dXCYll+QQjtlCjcBAcHEx4ebv1+2bJl1KxZ0/r9pUuXCoQdIUTFYl00U6e7Wkjc9C6oVqNUx0vPSWf036P5NfxXDDoDU7tO5ZWOr9h2FmAbqudVj6V3LqV9QHtSc1J5bv1z/HL6l9IfsN0QQAdnNuOdeR6QoeBCaKFE4WbQoEHExFxdkO6uu+6y9uQArFmzhk6dOmnXOiGEXeWFGzd9DhxYob7YflipjpWancqov0exO2o37iZ3Pu/zOQMbD9SqqQ7j7eLNl7d/yT0N7sGsmJm8dTKrTq4q5cHqqPPeAE0urgZkVXAhtFCimammTZt23fcnT56MwVAx/yITQlwNN331uyErCbzqQL2eJT5OcnYyo/8azYHYA3g4efBFny9oWaOlto11ICeDE9O7T8fV6Mry48uZum0quUpu6cJb2yFw6m8aXPwFPb2k50YIDZR42s3k5GR27txJdnY2nTp1okaNq93Vbm4V6x66EKKgvHDzoH6j+kLbx0BfspFMSVlJPPPXMxyOO4yXsxdf3v4lzas317qpDqfX6ZnceTJGvZElR5fw1va3yLXkMrjp4JIdqEl/cPXFLSOGW/X7sVjK15B4ISqiEn1qhYWF0bRpU/r168eAAQNo2LAhf/zxh63aJoSwM4uiEKSLoQuHAB20KVkhcVJWEiP+HMHhuMP4OPuw8I6FlTLY5NHpdEzoOIFhzdVbd+/ufJfvjnxXsoMYnaHVIwA8bNgsBcVCaKBE4WbChAnUq1ePrVu3EhoaSu/evRkzZoyt2iaEsLNci8JAwxb1m/o91JqQYsrIzeC59c9xLP4Y1V2qs6jvIpr4NrFRS8sPnU7HSx1eYkTLEQC8v/t91pxeU7KDtH0MgD76UNxyErRuohBVToluS4WGhvLnn3/Srl07ABYtWoSvry/Jycl4enrapIFCCPsx5+bykGGz+k3eMOViyLHk8NKml9h/eT+eTp7Mv2M+DX0a2qiV5Y9Op+OFti+QY87h6yNfM/XfqXg7e3Nr0K3FO0BgSxK9b8I78TBd0tYDfW3aXiEquxL13MTHxxMUFGT93tvbG3d3d5nbRohKwvX8v9TWxZGMuzoEvBgsioVp/07jnwv/4GJwYU7vOTTyaWTjlpY/Op2O8R3GM6D+AMyKmZc2vURYTFix9z9f90EAeqT9AXJrSogyKfGc50eOHOHAgQPWh6IoHD16tMBrQoiKyef4cgD+NtwCJtcbbK2atWcWv4T/gkFnYGbPmbTxb2PDFpZvep2eN7u/yS21byHTnMlz658r9kzGUXXvIksxUTf3rLoSuxCi1Eo8Wqp3796FVsa9++670el0KIqCTqfDbDZr1kAhhJ2kx+N1Th0g8LupDw8UY5evDn3F10fUyf7e7PZm8W/DVGImvYmZPWcy8s+R7L+8n2f+foYldy654cKgFmdv1lk6cq9hG+z7Dmq3s1OLhah8ShRuzpw5Y6t2CCEc7dCP6C3ZHLXU4bTxxvUyGyI2MCt0FgAvtX+Jexvea+sWVhiuRlfm9J7DsN+HcTrpNM9veJ6v+3193SUnjAYdK8w91HBzcCX0fafYvWdCiIJKFG7q1q1rq3YIIRxt37cALDf3xHCDuW2Oxx9n4j8TUVB4pMkjPNHiCTs0sGLxcvbi8z6fM+jXQRyLP8bkrZOZ2XPmNVdA1+t0bLPcRLTen4CsGDi6FlpV/BmdhXCEEoWb4tbTtGrVqlSNEUI4yKUDcGk/Fr2J1ebuBOqvvQJ4XEYcz294nozcDDrX7MyEThPs2NCKpWa1mnzc62Oe/ONJ/o74mzlhc3i+7fNFbmvU61HQ86dTH4ZkLoV930i4EaKUShRu2rRpY62tuRapuRGiAgpbCkBs7d4knvSg9jXCTbY5mxc3vciltEvU8ajDzB4zMelN9mxphdPGvw1vdHuDyVsn8+WBL2ng1YA7699ZaLu8zrI/jL0YwlI4swUSI0o015AQQiU1N0JUdbnZcFBdJPNiyANwEoxFhBtFUXhr+1vsi9mHh8mDT3t/ipezl71bWyHd0+AeTiWeYvGhxUzdNpVgj+BCa20Zr6Sbi/hDvVvVcBP2PfSUnjEhSqrY4ebAgQO0aNECfTHXmTl8+DBNmjTBaCzxgCwhhD2d/APS46BaAJf9uwP70RcRbpYeW8rPp3/GoDPwUY+PqO9V3/5trcDGth1LeGI4m89vZtzGcawYsILqrtWt7xuufLSaFQXaPHYl3CyBW18p8fpeQlR1xU4ebdu2JSoqqsBCmdfTtWtXwsLCqF9fPgCFKNeu3JKi1SPkYgDAoCsYbsJiwvho90cAvNThJbrV7mbXJoLac5QbFUX22bNknztH9pmzZEdEYE5OxpKWhiU9HUtaGkpODno3N/Xh7o7ezQ1jgD9OISE41a2LU0gIziEh6N3d7dp+g97A+7e+z6O/Pkp4UjgTtkzgi9u/wKBXr7n+yjXPNSvQbAD8+jIknoOI7RDS3a5tFaKiK3a4URSFKVOmFHvl7+zs7FI3SghhJ6kxcOLK4rdtHsMcpdbTGfL13MRnxvPy5pfJVXLpG9KXx5s9bpemKWYzWcePk75nD+l7QknfswdzfHyx9rUkJ19/A70e56ZNcGvfAbcOHXDr0B5j9erX30cD7iZ3/tfzfwz6dRA7o3YyJ2wOL7R7Abh6W8qiKODkDjfdp45gC1si4UaIEip2uLn11ls5fvx4sQ/ctWtXXF1ljgYhyrWDP4Bihtrtwb8p5osXgKvhxmwxM2HLBKLTo6nnVY83u72JTnftkVRlpeTkkLZjB8m/ryNl/XosSUkFNzAacQoKytcLUxeDb3VrD43e3R2dyYglPQNLeprao5OaRs7Fi1d7fM6dwxwXR9aRo2QdOUrCt+oQeJfmzfHo3w/Pfv1wCg622c9Y37s+b3V7i1e2vML8g/NpXaM1PYJ7WO88mS1XBmy0eUwNN4dXQ/8PwLmazdokRGVT7HCzadMmGzZDCGF3igL7lqjP2zwKXP3Fmhdu5u6fy45LO3A1ujKrxyzcTdrfylEUhYzQUBJXryb1r78x5ws0end3XNu3w61DR9w6tMelRQv0Tk5lPmdOTAwZoaGk795D+p49ZJ04QeaRI2QeOcLlmbNwadECz/798br/Poy+vmU+33/1q9ePfTH7WHpsKZO2TmL53csx6r2BfOGmThfwrQ/x4XB0jfX/IyHEjUm1rxBV1aX9EHMYDM7QQl20MX+42XJ+C18e+BKAaV2nab7KtyUzk+S1a4n/bglZx45ZXzf4+eF5x+149OuHW7t26GwwKMHk74+pf388+/cHIDc+npS//yZl3TrSduwk89AhMg8d4vLHH+N51134PP4YrjfdpGkbXu7wMofiDnHg8gFe2vQSb3SYC1wpKAbQ6dRAs+FtNYRKuBGi2CTcCFFV5RUSN70LXH2Aq+HGok/kta1vATCoySDuql+8FcKLIzcujvivviJxxQ/WXhqdiwteA+7G8+4BuHVoj85g0Ox8xWH09cXn4YfxefhhNej8+ReJP/5I5sGDJK1aRdKqVbi2bYvv8CfwuP12TW7NmQwmZvaYycO/PMzR+KN8fXwO0AGzOd88Yq0GwYZ34NxWiD8DvvXKfF4hqgIJN0JURblZ1rltaPOY9WW118DCaeaTnJVE8+rNeaXjK5qc0pyYSNzCRcQvWYKSng6AqXZtfB59FO8HH8Dg7a3JecrK6OuLz6BH8Bn0CBn79xP/3RKS160jY98+Luzbh3PzZtR4/nmq9exZ5pAT6B7IjFtmMOrvUfx2biXGaq6Ys/PN8O4dDPV7QPgm2L8Mek0q2w8nRBUhkycIURWdWAcZCeBRExr0sr5stig4Vd9EMsdwNbrywa0f4GQoW42LOSWFy59+xqk+txM3fz5KejouLVoQNOczGvz5B9WferLcBJv/cm3dmtoffkCjDeupPuoZ9G5uZB05yvnRz3J20CDStm0r8zm61+7OEzc9AYBLzR/J1ScU3KDNldFpYUvBYinz+YSoCiTcCFEV5d2Saj0I9FdvAUWmHcGpxt8AvN7ldep6ln6xXMViIXHVak7360/snDlYUlNxbtqUoLlzCPlhBR69e9v99lNpGWvUwH/cOBqs/5vqI55C5+JC5v4DRDz5FJGjRpMdEVGm47/Q9gUaeTdDZ0zHGLAMsyXfEjZN7wJnT0iKgLP/lPEnEaJqkHAjRFWTEg0n/1Kft75apJqcncyv0TPR6SwE6LsyoP6AUp8i88gRzj32OJcmTcIcF4dTSAi1Z/+Pej/9iMdtt9l0OLktGX188H/5ZRr+9Sc+Q4aA0Ujqpk2E3z2Ay598giUjo1THNRlMvN7pbRSLEwa3cBYcXHD1TSc3uOl+9fn+7zX4KYSo/CTcCFHVHFyhzm0T1BFqNAbU4djTt08nJTcGS7YvzZ2GlyqAmFNTiXprOmceGkjGvn3o3Nzwf/kl6q/5Gc9+/dBVkmUEjDVqEDj5Neqv+Rn3bl1RsrOJnfs54XfdTcrGjaU6ZohnCJlR9wLw+f7PCYsJu/pmXl3UkZ8hK6WMrRei8qscnzRCiOJRlKu3pPINLf759M+sO7sOPQYyLgzCWV+8mcjzS9u+nfB77iFhqVob4nlnfxr89ivVR4xAp8HcNOWRc/36BC9cSO2PP8ZYsyY5Fy9yfvSzXJwwscB8PcVh0OnITWpHTlIbzIo6eWJqdqr6ZnAn8G0AOelwZI0NfhIhKhcJN0JUJZfCIOaIOrfNTQ8AcCH1Au/teg+A9p6DsGTWwWgofq+NJS2NS2++ScTwJ8m9eAlTcDB1vlpM7VmzMAUG2uKnKFd0Oh2efe+gwa9r8X3qSdDrSfr5Z8IH3EPqli3FPo7BoAN0ZEbdRy332lxMu8j7u9/PO8nVMJoXToUQ1yThRoiqJOxKzUazu8HVG7PFzOStk0nLSaOtf1taVFNvi+iLeUsqfc8ewu+7n8TvlwHg8+hg6q9ehXuXLjZpfnmmd3Mj4JVXqLvkO5zq1iU3JobIp5/h4uuvY05Nu+H+1sVKLS5M7fwWOnSsPrWa9RHr1ddbDwJ0V+e8EUJck4QbIaqKAnPbqL0A3x75ltDoUNyMbrxz8zsoivqRYNRfP9woFgux8+ZxbugwciIjMdaqSZ3FiwicOtXuq22XN25t21Jv9Sp8hw0FnY6klT9y9qGHyLzB2nz5Fytt4deW4S2GA/DmtjeJzYgFryCo31PdYP8yWzVfiEqhXISbOXPmEBISgouLC507d2bXrl3F2m/ZsmXodDruu+8+2zZQiMrgxB9X57ap34vj8cf5ZN8nALza8VWCPYKtMxTrrxNucmNjiRwxksuzPwaLBa9776H+mjW4d+1qlx+jItC7uhIwaRJ1vv4KY0AA2WfPcnbgwyQsX4GiKEXukz/cWCwKz7V5jsY+jUnISuCNbW+o++UVFu+XOW+EuB6Hh5vly5czfvx4pk2bxt69e2ndujV9+/YlJibmuvudPXuWl19+mVtuucVOLRWigss3t022Yua1ra+RY8mhZ1BPHmik1t9Y15a6xm2ptB07Cb//ftK2bUPn4kLNd9+l1vvvY6gmK1YXxb1TJ+qtXoV7j1tRsrOJmjaNiy+9jDk1tdC2+fNkrkXByeDEjFtmYNKb2Hx+Mz+d/OnqnDeJEXDuXzv+JEJULA4PN7NmzWLkyJEMHz6c5s2bM2/ePNzc3Fi0aNE19zGbzTz22GO8+eab1K9f346tFaKCSo2Bk3+qz1s/ypywOZxIOIGPsw/Tuk2zDvu2hpv/FBQrikLcwkVEPPkk5suxODVsQL0fVuD9wP12/TEqIqOPD8Gff47/Ky+DwUDyb79xduDDZIUXrJvR6XTW3hvLld6dxj6NGdtuLADv736fyKy4q3PeSGGxENfk0HCTnZ1NaGgoffr0sb6m1+vp06cP27dvv+Z+b731Fv7+/jz11FM3PEdWVhbJyckFHkJUOQeuzm0TpqTz1eGvAJjWbRp+rn7WzfJWpM7fc2PJzOTiqxOI+fBD9TbUffdRb8UKnBs1suuPUJHp9HqqP/UUdb/9FmNgINlnznD2kUcKjaYy/CdkAgxpPoQOAR3IyM1gyr9TsLQerL5x5GfIKtwDJIRwcLiJjY3FbDYTEBBQ4PWAgACioqKK3Gfr1q0sXLiQ+fPnF+scM2bMwMvLy/oIDg4uc7uFqFAUBcKWAJDZ6mH1F6RiYUD9AfSu07vApnm/VPMKinOiojj3+BCSf/kFDAYCprxOzRnvoncr+Tw4AtzataXeyh9wbdcOS0oKkc+MIm7BAmsdTl7PTf5wo9fpmd59Oq5GV0KjQ/k+7fSVOW/S4KjMeSNEURx+W6okUlJSGDJkCPPnz8fPz+/GOwCTJk0iKSnJ+oiMjLRxK4UoZy7tt85t85n5MmeTz1LDtQYTOk0otGn+guL0ffs489BAMg8dwuDtTZ2FC/F97LEKu3RCeWH086PuV4vxHjgQFIWYj2Zy8eVXsGRmFhluAII8gnip/UsAzN77MRE33a2+sW+JXdsuREVhdOTJ/fz8MBgMREdHF3g9OjqawCIm/zp9+jRnz55lwICra95YrowYMBqNHD9+nAYNGhTYx9nZGWdnZxu0XogK4kqvzb5Gt/LNieUAvNHtDbycvQptmvdLtVbYv0R8MxslOxvnxo0JmjsHp6Ag+7W5ktM5ORH41ps4N2tK9LszSP71V3IuXMC7/kBScSbXUnhE1cAmA/kr4i92XtrJlIwTLEaHPm/OG996DvgphCi/HNpz4+TkRPv27Vm/fr31NYvFwvr16+laxLDSpk2bcvDgQcLCwqyPe+65h169ehEWFia3nIT4r9wsOPgDGTodU3TxKCjc2+Bebg26tcjNzWYLD5zcRMsFH6BkZ1PtttsI+X6pBBsb0Ol0+D76KHUWLEDv6UlGWBjT/5pNzdRYa0Fxfnqdnre6vYWb0Y29cYdZUq+N+oYspilEIQ6/LTV+/Hjmz5/P119/zdGjRxk9ejRpaWkMH65OYDV06FAmTZoEgIuLCy1atCjw8Pb2xsPDgxYtWuBUSdevEaLUjv8GGQl8ElCbc5mX8Xfz59VOrxa5qWI203ntYkYeXguAz2OPEfTpJ1V+Uj5bc+/SmZDvl2KqVYvAlMv8b8un5B44UOS2tarV4qUO6u2pj3VJnDUa1VmnZc4bIQpweLh55JFH+Oijj5g6dSpt2rQhLCyMdevWWYuMIyIiuHTpkoNbKUQFFbaUvc7OLHFV/6m/0fUNPJ08C21mycjg/JjnuWmXOlz8zMARBLw+GZ3BYNfmVlXODRoQsnwZZ6oH45WdBi89R/Kffxa57cDGA+lSswtZSi5TAvwxJ0XA2X/s3GIhyjedcq3pMiup5ORkvLy8SEpKwtOz8Ie8EJVG8iUyZzdnYM0AzjqZuL/h/bzV/a1Cm5mTk4kcNZqMvXvJNZh4v90geo8azPDuUsdhb7dN/42hfy+gc/RR0OsJfGMaPg8/XGi7S6mXuH/N/aTlpDExLp7HQu6CB750QIuFsJ+S/P52eM+NEMJGDizjcy8PzjqZ8Hf15+WOLxfaJPfyZc4NHUbG3r3oPT1Z9dhEttZufcO1pYRt5Dq78lbnJ8jtfw9YLERNnUZsEdNe1KxWk/HtxwPwsY8354+vhcwkezdXiHJLwo0QlZGicHj/t3zlpf5183qX1wvdjso+f56zjz1O1rFjGPz8qPvtN5wPagxcf20pYTsGvQ6L3kDmC69SfeRIAC7PnEXMRx8VWpPqocYPqZP76fW84eOOcugnRzRZiHJJwo0QlVBOxHamOKVh0enoX6cPver0KvB+1smTnHv0MXIiIjAFBRGydAkuTZoUmsRP2FfePDe5FvB/aby6ZAMQt2AhUVOnopjN1m31Oj1vdnsTF52Bna4u/LR/gUPaLER5JOFGiEpowc73OOnkhI/OyMSuUwq8l3nkCOeGDCU3JgbnRg2pu2QJTnXqAFeXX9DLRH0Okbf8Qt5Q8OpPPUXg9LdAryfxh5VcfOVVlJwc6/Z1POswpoW6DM1H+kSiz++wf6OFKIck3AhRyZyMOciX6eEATGoyBF8XX+t7GQcOcO6J4ZgTE3Fp2ZK6336LKcDf+r6158Yg4cYRrvbcXL0F5TNwILVnzQKjkeTffuPCSy+jZGdb33+8zbO0xIVUvZ63/51a6PaVEFWRhBshKhGzxcy0zS+Tq9PRMxv6dRhrfS997z4ihj+JJTkZ17ZtqbN4EQZv7//sLz03jmRdFfw/MxR79utL0CefoDOZSPnzT86PHYflSsAx6A281WIkRkVhU+Ylfg//1e7tFqK8kXAjRCWy5OgSDqZfxMNs4fWQ+6zz1KTt2kXEiBFY0tJw69iROgvmY6hWrdD+eT0GBqm5cQh9ET03eTxu60XQ3DnonJ1J3biR88+NwZKZCUDD1k/wTKoadt7b8TYJmQn2a7QQ5ZCEGyEqifMp5/ls3ycAjE9IIKCDWouRtmMHkU8/g5Kejnu3rgR/+cU1Zx22SEGxQxmvsXBmnmq33ELwvM/RubiQ9s8/nH/2WTXgGJ14qv69NMzOJiE3jQ92f2DPZgtR7ki4EaISUBSFt7a/RYY5i44ZmTzo1x58QkjbtYvIUaNRMjNxv/UWgj7/HL2r6zWPIwXFjvXfguKiuHftSp35X6JzcyNt23bOP/sclqwsTO2G8dblePSKwtrwtWy9sNVezRai3JFwI0QlsOb0GrZf2o6zojAtNh5du2Gk79lD5DOj1GBzyy0Effopemfn6x5HCoodq6iC4qK4dexInS+/uBJwtnF+zPNYfBrR0u8mHktOAeCt7W+RnpNu8zYLUR5JuBGigovNiLXehhidkEhdYzXSM4KIePoZlIwM3Lt3J+izGwcbkIJiR7tWQXFR3Dp0UG9Rubqqt6heeAFLi8GMSUiitkXHpbRLfHLlNqUQVY2EGyEquPd2vUdydjLNdK4MS0oh3fN2Ikc9Z62xCZrzWbGCDeTrudHLR4MjFLfnJo97p04Ef36lBmfzFi58tRtXXJgaEw3A0qNLCYsJs1VzhSi35BNMiApsY8RG/jj7BwadnjcuRpATbyJy4V4s6em4delC0Jw56F1cin08a8+NfDI4REl6bvK4d+lM8Odz1VFUW/7lwoFGdE3L5B6nQBQU3tj2Btnm7BsfSIhKRD7ChKigUrNTeXvn2wAM9W5F/agcIrf4Y0nLwLVDe4Lnzrlu8XBR8sKNQW5LOUTe7cDi9tzkce/aVR0mbjKRcjiWizu9eSX8GL7OPpxOOs3CQwtt0Vwhyi0JN0JUULP3ziYmPYZgj2BGHjxGxKbqmDMVXFq1InjeF+jd3Ep8zLzRUlJQ7BjWoeClmGW4Wvfu1P74Y3Um43NuZPyrY5L/LQB8eeBLTiee1rStQpRnEm6EqID2xexj+fHlALzpeSdRK5MwZxpwbtyIOvO/xFCt6HlsbiTXLAXFjpR3W8pstpRqf4/belH7ww9ApyMx3J02C/+kR+1bybXk8sa2N7AopTuuEBWNhBshKphsczbTtk0D4FG//vhMXURuhgGnGq7UWbwYg5dXqY+dN7+KFBQ7hjXclGF5KM/+/ak55RUAEkJTeWmrK25GN8Iuh7Hi+AotmilEuSefYEJUMPMPzudM0hnqWHx4aM4BcuIzMLnnUmfWNIzVq5fp2FJQ7FilKSguivejwwnoXwuAzGW/8u65joB6KzMqLapsjRSiApCPMCEqkJMJJ1lwcAGuWQrv/uRCbvg5jK5m6tzniqn9gDIfX4aCO1ZJh4Jfj+/T4/BvnQxArW/W88SpWqTlpPH2jrdl5XBR6cknmBAVhNli5o3tb6DLyuHdXzxxOhmJwVVPnZ5xOPUYqkl3S14hq0E+GRyiOMsvFFvjflTv4Eb1ZuqMxXeujOTmYzo2n9/MH2f/KPvxhSjH5CNMiApi2fFlHI7az8trdNQ+mYDezZXgW6Jx9gHaDtHkHGYpKHYo66rgZSm6yWMwQdvHqdEqBe923qAojFljpu0pCzN2zSAxM7Hs5xCinJJwI0QFcDH1Ip+EzubZXy20PZGLztmZ4OFtcfXNgSZ3gkeAJucxS0GxQ5VlKHiR2g9Dp4PARkfwvL0nerOFl1cp+J+M46M9H2lzDiHKIfkEE6KcUxSF6dvf4pHf07jlsAJGI7Vnvo9b6l/qBu2f0OxcuVJQ7FDW0VIWjYZs+4RAg9vQ6aDW3f5U69kTU67ChB/MhG1bzfaL27U5jxDljHyECVHO/XrmVwJW/EP/UDV41JoxAw+/OMhKAu+6UL+XZueySEGxQ10NNxoetP1wAHQHllJ75ge4dmiPexZMXm7ms7Wvy8rholKSTzAhyrH4zHh2zHmDR/5Rf9sFTJ6M14C7Yc9idYP2wzTtZsm7HSLZxjGsQ8G1HM3UpD9UC4C0GPTnNhA8dy6mJo3xToMRiy6ycJPcnhKVj3yECVGOrZg7jsG/pgHgM/oZfIc8DtGH4fwu0BuhzeOanctiUcj7nSo9N45h0LKg2HpQtbAYgNCvMHh6ErJgAeZaNQhMhEbTl3H4zC7tzidEOSCfYEKUUzt//oLuC3ajByz330HAC2PVN0K/Ur9qWEgMBYtYZeFMx9B0KHh+7YYCOji9AeLPYKxRg8ZfLyHd05m6MQqRo0eRlZas7TmFcCAJN0KUQ/F7d+I85WOMFrjQOYTm7/wPnU4H2emwX11TSstCYrg6gR+AQRbOdIirk/hpvAbUlcJiAPZ+A4BTcDBBC74kzUVH3bMZ7Bk5GCU3V9vzCuEgEm6EKGeyws8Q8fQzOGcrHG/oQtfPl6HLu010eJVNConhP+FGem4cwiYFxXnywvC+78CcA0BAq04kTh9NthF894ZzYuJ4mb1YVAoSboQoR3Kiozk1fAjOqVmcqgkBs2fi7pZvIcxQ2xQSQ8Ep/6XkxjE0HwqeX77CYo79an25991j+OOplph1YFn7FzEzZ2p/biHsTD7ChCgnzImJnHvqKfTRcVz0hQMT7qFrw9uubnBpP5zfrXkhcZ78izVKQbFj2LTnxmC6OpP17gXWl3U6HUOf/oTFd7sAEL9gIXGLv7JBA4SwH/kEE6IcsGRkEDn6WXJOnSa+Gsx9ogZjer1WcKNd89Wvze/VtJA4T/6CYr3clXIImxUU5+kwHHR6OPsPxByzvhzoHkinpyaypKf6KyHm/fdJ+vln27RBCDuQcCOEgym5uVx4cTwZ+/aR6gLvDDLwbN9peDnnux2VHg8HV6rPO460STvyam4Mep1avCzsTstVwYvkFaSOsgPYPb/AWwObDCTynvb80kltw8XJr5O6ZYtt2iGEjUm4EcKBFEXh0tRppG7aRI5Rx/sPGWja4Q561+1dcMOwJZCbAQEtoU4Xm7TFGm4k2DiMdRI/W4UbgE5XwvH+ZZB5dfi3Xqfnze5vsaKPC1tu0kFuLufHjiMjLMx2bRHCRiTcCOFAl2f9j6SffkLR65h1n45LDbx4rfN/bkdZLLB7ofq80wiwUfjI33MjHMNmQ8Hzq9cD/BpDdiocWF7grRCvEEa1fZbP79JzsKERJSODyGdGkXX6tO3aI4QNSLgRwkHiv/6auPnqrYH5/Y2ENtIzsdNE/Fz9Cm54ej0knAFnL2g50GbtkXDjeDYtKM6j0129tblrPvynvueJm56gqX8LPrhXITrEE3NSEhEjRpJz6ZINGyWEtiTcCOEASb+sJXrGewCsv7Mmf7dS6BHUg7vr311447xC4raPg5O7zdpkXRFcso3D5N0StMlQ8PxaDwKnahB7HM4UrKsx6o1M7z4ds4uJ1+5NIzuoBrmXLhExciTmxETbtksIjUi4EcLOUv/ZysVJkwC4dFd7vmgVg4fJgyldphQu5I0/Ayf/VJ93fMqm7coboWM0yMeCo1h7bmw9j56LpxpwoFBhMUAjn0Y80+oZUtx0THkoB71/DbJPnSZy9LNYMjJs3Dghyk4+xYSwo4wDBzg/dizk5qLv25NXWx8BnY5XOr5CgHsRw7v3LAQUaNgHqjewadvyFmvUS0Gxw9h0Er//6jhC/XrsV0g6X+jtp1o+RVPfppxxTWXFM03Qe3qSsW8fF8a9iJKTY/v2CVEGEm6EsJOs8HAin34GJT0dt5u7807vJLKUHLrX6s59De8rvEN2Ouz9Vn1uo+Hf+Vl7buS+lMNcDTd2WALBvxmE3AKKBfYsLvS2SW9ievfpGHVGVuTu4OK0J9A5O5O6eTOXpkyVZRpEuSbhRgg7yImKImLECMyJibi0asWWUZ3Zl3AQd5M707pOK3pemUMrITMRvOtAo9tt3kYpKHa8q0PB7XTCvN6b0K8gJ7PQ2019m/JUS/V26Bupy/D84C0wGEhavZrLskyDKMck3AhhY+bERCJHjiT34iWc6tXD8uFEZh+dB8D49uOpWa1m4Z0UBbbPVZ93HAl6g83bmSvhxuHsMhQ8v6Z3gWcQpMfCwR+K3OSZVs/QyKcR8ZnxfOS8icDpbwEQt2AhcYsK9/gIUR5IuBHChizp6eo8ISdPYfT3p9aX85hy6EOyLdl0q9WNgY2vMbT79Aa4fFQd0dJ+mH3aqki4cTTraCl73fExmKDz0+rzHXMLDQsHMBlMvNP9HYw6I3+d+4vtbVzwf/klAGI++IDEVavt1Fghik/CjRA2omRnqzO87t+P3suLOgsXsCTxTw7GHsTD5MGb3d689jIHO6702rQdAi5eRW+jsbyCYgk3jmMw2LGgOE+7YWByh5gjEL6xyE2aVW/G063VEPTOzncwD74H3+HDAbj0+uukbCh6PyEcRcKNEDagWCxcnPQaaf/8g87VleB5n3O2uoU5YXMAmNBpAoHugUXvHHMUTv2tLnDY+Rm7tdnacyOjpRzm6jw3djypqze0u7Ja+PY519xsRMsRNK/enOTsZN7c8SY1XnkZr/vvB7OZCy++SPru3fZprxDFIOFGCI0pikL0uzNI/vVXMBoJ+uRjTK1a8Pq/r5NryaVncE/uaXDPtQ+Q12vT9C7wrWefRpNvEj/puXEYuw4Fz6/zKECnhup8q4XnZ9Krt6dMehNbzm/h59M/U3P6W1S77TaUrCwiRz9L5rGi9xXC3iTcCKGx2LlzSfjuO9DpqPXee1S75Ra+PPglx+KP4eXsde3RUQBpsbD/yno/XcfYr9FcXaxRhoI7jl2HgufnW08N03A1XBehoU9DnmvzHAAf7P6AqMzL1J41E7cOHbCkphIxYiTZ587Zo8VCXJeEGyE0FP/td8R++hkAAZMn43X3XRy4fID5B9RZYF/v/HrhtaPy270QzFlQqx0Ed7ZHk63M0nPjcNah4I6YQiYvTO9fpobsa3jipidoVaMVqTmpvP7v6+DsRNDnc3Fu2hRzbCwRTz5FTnS0nRotRNEk3AihkaQ1a4h+5x0A/J4fg+/jj5Gek85rW1/DrJjpH9KffvX6XfsAOZlXp8Lv+pzNVv++llzpuXE4uw8Fz69OFzVUm7Ngz6JrbmbQG3j35ndxNbqyK2oX3x35DoOHB3UWzMdUtw45Fy4Q8dRT5CYk2LHxQhQk4UYIDaRs2MjFSa8B4DN0CH7PPgvAzD0zOZd8jgC3ACZ3mXz9gxz8AdIuq/OONL/X1k0uRAqKHS/v2jsi26DTqaEa1MVac7OuuWldz7q83OFlAD7e+zEnE05i9POjzsJFGAMC1HWonhmFOTXNHi0XohAJN0KUUdrOXVwYNw7MZrzuvZeAiRPR6XRsOb+FFSdWAPD2zW/j5XydId0Wy9Vah85Pq/OP2JlM4ud4Du25ATVUe9aGtBg4sOK6mw5sPJAeQT3ItmQz8Z+JZJuzcQqqTZ2FCzB4e5N54ADnx4zBknXtkCSErUi4EaIMMg4e5Pyzz6JkZ1Ptttuo+c7b6PR64jPjmfrvVACGNB9Cl5pdrn+gk3+o84w4eajzjjiARcKNw10tKHZUA0xXRk4B/34MFvM1N9XpdLzR7Q18nH04kXCCz/aptWbODRsSPP9L9G5upO/YwYUXx8tCm8LuJNwIUUqZx08QMWIklrQ03Dp3pvb/ZqEzGlEUhTe2vUFcZhwNvRsytt3Y6x9IUWDLR+rzjk+p8444gAwFdzyHDQXPr8NwcPGGuJNw9Jfrburn6scb3d4A4KvDX7E7Sp3rxrVlS4LmzkXn5ETqhg1cnPQaiiN/JlHlSLgRohSyz54l4qmnsCQl4dq6NcFz56B3dgbgp5M/sTFyI0a9kfdueQ9ng/P1D3b2H7iwBwzO0OVZO7S+aDIU3PEcNhQ8P2ePq5NH/jOzyCUZ8rutzm080OgBFBQmb51MUlYSAO5dOlP749lgNJK8di1Rb74lK4kLu5FwI0QJ5Vy8yLknn8QcG4tz06YEf/kFend3AMITw3l/9/sAPN/2eZr4NrnxAf+ZpX5tNwQ8AmzV7Buy9txIQbHDXJ2h2MEhoPModUmGqANwev0NN5/QcQLBHsFcSrvEW9uvhhiPXr2o/cH7oNeTuHw5MR98KAFH2IWEGyFKIDc2lojhT1pX+K6zcAEGL7VQOMucxatbXiUjN4MuNbvwxE1P3PiAF0LV9Xx0Buj2gm0bfwNmRXpuHM3ac+PoAODmC+2fUJ/nhe/rbW5y44NbP8CoM/LnuT/56eRP1vc877yTmldWEo9fvJjYudeeJFAIrZSLcDNnzhxCQkJwcXGhc+fO7Nq165rbzp8/n1tuuQUfHx98fHzo06fPdbcXQiu5CQlEDB9O9rlzmGrVos6ihRirV7e+/7/Q/3E84Tg+zj68e/O76HXF+OeV94uj1cPgU9dGLS8eKSh2POskfuWhPKXbGNCb4Ny/cG77DTdv4deC59s9D8B7u94jPDHc+p73gw8SMGkiALGffkbcwoW2abMQVzg83Cxfvpzx48czbdo09u7dS+vWrenbty8xMTFFbr9p0yYGDx7Mxo0b2b59O8HBwdxxxx1cuHDBzi0XVYk5KYmIJ58i6+QpjP7+1Fm8CFPNmtb3N0duZsnRJYA67LuGW40bHzTmGBxbqz7vPs4GrS4ZGQrueEZHDwXPz7MWtHlUfb71xr03oM5e3KVmFzLNmby65VWyzFeHgfsOG0aNcWpxfcyHHxH/7XeaN1mIPA4PN7NmzWLkyJEMHz6c5s2bM2/ePNzc3Fi0qOgZMpcsWcKzzz5LmzZtaNq0KQsWLMBisbB+/Y3vCwtRGuYra+ZkHT2Kwc+POl8txqnu1V6WmPQYpvw7BYDHmz3OrUG3Fu/A/85Wvza9G/ybatzqkpOeG8fT51t+oVzUpnQfq65Of/JPuHTghpvrdXrevfldfF18OZ5wnP+F/q/A+36jRlF9tDrUPPqdd0hYfv25dIQoLYeGm+zsbEJDQ+nTp4/1Nb1eT58+fdi+/cbdoADp6enk5OTg6+tb5PtZWVkkJycXeAhRXJa0NCJHPk3mwYMYvL2ps2ghzvXrW983W8y8tvU1ErISaOrblBfbv1i8AyecuzpJ2i3jbdDykpOCYsfLPzu0w4uKAao3gJvuV58Xs/emhlsNpnefDsCSo0vYFLmp4PsvvIDvU08CEPXGGySuWq1RY4W4yqHhJjY2FrPZTEBAwREiAQEBREVFFesYEyZMoFatWgUCUn4zZszAy8vL+ggODi5zu0XVYElPJ3L0s2Ts24fe05M6ixbi0rhxgW2+OPAFOy/txNXoyge3foCTwal4B/9nJihmqN8TarfXvvGlYJGCYoczGPKFm/LQcwNw85XwfXi1eiu1GG4NupUhzYcAMHnrZC6kXi0b0Ol0+L/8Mj6PPw6KwqXJk0n6Za3WrRZVnMNvS5XFe++9x7Jly1i1ahUuLi5FbjNp0iSSkpKsj8jISDu3UlRElowMIkc/S/quXejd3amzYD4uzZsX2GbbxW3M2z8PgCldplDPq17xDh53GvZdqTfoOUnLZpdJrlkm8XO0ctdzAxDYApoNABTY9G6xd3ux3Yu09GtJcnYyL296mRzz1VmKdTodAZNfw/vhh8Fi4eKECSSt/dUGjRdVlUPDjZ+fHwaDgejo6AKvR0dHExgYeN19P/roI9577z3+/PNPWrVqdc3tnJ2d8fT0LPAQ4nosGRlEjhpN+s6d6N3dCV4wH9f//DcWnRbNpH8moaDwUOOHGNBgQPFPsPkDtdem4e3qSszlhAwFd7z89U7lJtwA9HwN0MGRn+HS/mLtYjKY+KjHR3g6eXIo7hAzQ2cWeF+n0xH4xjS8HnpQDTivvkrSrxJwhDYcGm6cnJxo3759gWLgvOLgrl27XnO/Dz74gOnTp7Nu3To6dOhgj6aKKsLaY5Mv2Li1bVtgmxxLDq9ueZX4zHia+jZlYqeJxT/B5eNwYLn6/LYbrBJuZ1JQ7Hj5r315GDBlFdAcWj6kPt9Y/N6bWtVq8e7N6vZLji7hj7N/FHhfp9dT8623rgacVyTgCG04/LbU+PHjmT9/Pl9//TVHjx5l9OjRpKWlMXz4cACGDh3KpElXu+7ff/99pkyZwqJFiwgJCSEqKoqoqChSU1Md9SOISsKSkUHks8+SvmMHejc3gucXDjYAn+77lL0xe3E3uTOzx8wbL6+Q36YZgKKOkKpV+NiOJEPBHS//balyMRw8vx4T1ckmT6yD83uKv1twD55soRYQT9s2jXPJ5wq8bw04Dz5gDTjJv/2madNF1ePwcPPII4/w0UcfMXXqVNq0aUNYWBjr1q2zFhlHRERw6dIl6/aff/452dnZPPTQQ9SsWdP6+Oijjxz1I4hKwJKWpt6K2n4l2CxYgFu7wuFjY8RGFh9aDMD07tOp41mn+CeJOgiHVwE66PWaRi3XTl5BsYQbx9HrdeTlm3JTUJzHryG0Gaw+3/B2iXZ9vu3ztPNvR1pOGuM3jScjN6PA+zq9nprTp+P1gBpwLrz8Cklr1mjVclEFGR3dAIAxY8YwZsyYIt/btGlTge/Pnj1r+waJKsWcmkrk08+QsXeveitq/vwig014UjiTtqq9iI83e5zb695eshPldee3eAACbiprszWXV1As4caxjHodOWalfNXc5Ln1Vdi/XF0y5OxWCLm5WLsZ9UY+7PEhA38ZyImEE0zbNo33b3kfXb6eKp1eT823p4NeR9LKH7k4YSJKTg7eDz5oq59GVGIO77kRwpHyZh7O2LtXHe69eFGRwSY1O5WxG8aSlpNGO/92jO9QwrlpzofC8d/UCdHK0Qip/Kw9NzLPjUPpy8vimUXxqQvthqrPN7xzwxXD8/N382dmj5kYdUZ+P/M73xz5ptA2ebeovAcPujJM/HUSvv9eq9aLKkTCjaiychMSODd8OJkHDmDw9qbuV4sLjYoCsCgWXtv6GmeTzxLgFsDMnjMx6U0lO9kGdVIzWg8Gv0YatF57eTUeMhTcsayLZ5bHcANw68tgcIaIbcVaMTy/DoEdeLnjywDMCp3Fzks7C22j0+sJnDoV32FqiIp68y3ivvqqzM0WVYuEG1El5cTEEDF0GFlHjmKoXp06X39daB6bPF8c+IKNkRtx0jsxu9ds/Fz9SnayU3+r3fh6E/R4VYPW24b5Sv2qDAV3rHIfbjxrQaeR6vO/poHFXKLdH236KPc0uAeLYuHlzS8XmOAvj06nw3/iRKqPVM8T8977xM6bVz6WpBAVgoQbUeVknz/PucceJ+vkSYw1alD3229wadK4yG03R25mbthcAF7v8jot/FqU7GTmXPjjdfV552fAJ6QMLbctGQpePlhXBi/Pv8hveQlcvCH60NUJKYtJp9MxpcsUmldvTmJWIuM2jitUYJy3XY3xL+L33HMAXJ79MTEffiQBRxSLhBtRpWSdPMm5wY+SExmJKTiYukuXFFgrKr9TCaeY+I86h82gJoO4v9H9JT/hvm/g8lFw9VG788sxGQpePlxdGbwc/xJ384UeE9TnG96GrJQS7e5idOHjXh/j6+LLsfhjTPl3Chal8NB3nU5HjefH4D9BPVf8okVETZ2KYi5Zb5GoeiTciCoj4+BBzj0+hNzLl3Fu1Ii6S77D6RprjcVlxDFmwxhSc1JpH9CeVzuW4nZSZrJadAlqEbGrTxlab3tSUFw+lOuC4vw6jgDfBpAWA1tnl3j3QPdAPurxEUa9kT/O/sHn+z+/5rbVhz9BzXfeBr2exB9WcuGll1Gys8vQeFHZSbgRVULajh1EDHsCc1ISLq1bUffbbzD5+xe5bZY5i3Ebx3Eh9QLBHsHM7jkbk6GEBcSgLo6ZHgvVG0GHJ8v4E9ie9NyUD8byXnOTx+gEd1wplN/+GSSWfN2+joEdmdplKgDz9s9jbfi1F9D0fvBBas+aBSYTKevWEfncGCxpaaVquqj8JNyISi/p11+JGPk0lvR03Lp2oe6iRRi8vYvcVlEUpv47lbDLYXg4efBZ78/wdil62+tKOAs71Fod7ngbShOO7ExqbsoHfUUJNwBN7oSQWyA3E9a/WapD3N/ofoa3UGekn/rvVMJiwq65rWe/vgTPnYvO1ZW0f/7h3BPDyY2PL9V5ReUm4UZUanFffcXFl16GnBw8+vUj+Isv0Lu7X3P7Lw58wW9nfsOgMzCr5yzqexVdj3NDf78J5myo1wMa9y1l6+1LhoKXD+V+tFR+Oh30fQfQwcEfSrQsQ37j2o3jtuDbyLHkMHbjWM6nnL/mttVuuZm6i9U/UDIPHuTs4MFkR5a810hUbhJuRKWkWCxEv/8BMe+9D4DPkCHUnjUTvZPTNff5Lfw35oTNAWByl8l0qVnKFbsjdsLhn4ArH/wVpIZFhoKXDxUq3ADUbA1tHlOfr5tUqhU/9To9M26ZQTPfZsRnxjNm/RiSspKuub1rmzbUXboUU+3a5JyL4OzgR8k4fLi0P4GohCTciErHkp3NxVcnEL9YXQPK/+WXCHhtEjr9tf9z335xO5P/VVfpHtp8KAMbDyzdyc05sPZF9XnbxyGwZemO4wDmK7+UpKDYsQwVpaA4v9teB5M7nN8FYSUbGp7HzeTGp7d9ir+rP6eTTvPChhfIzM285vbO9etR9/ulODdtijk2loghQ0n9Z2tpfwJRyUi4EZVKbkICEU8+SfLatWA0UvO9GVQfMaLAGjb/dTTuKOM2jiPXkkvfkL681OGl0jdgx1yIOQyuvtCndDUIjnJlaSmpuXEwa89NRZrPxbPm1cVg/5wCqZdLdZgA9wDm9pmLh8mDvTF7mfjPRMzXmSTQ5O9P3e++xa1LFyzp6USOGkXCsuWlOreoXCTciEoj68wZzg4aRMaeUPTVqhH8xTy877vvuvtEpkQy+u/RpOem0ymwE+/e/C56XSn/WSScg40z1Od3vA3u1Ut3HAeRguLyocLdlsrTeZTaU5mZCH++XurDNPFtwse3fYxJb2J9xHre3fnudSfuM1SrRp0vv8Dr3nvBbCbqjTeIfv8DlFLcHhOVh4QbUSmk7drFuUGDyTkXgalWLUK+X0q17t2vu09cRhyj/hpFXGYcTXyaMLvXbJwM167JuS5Fgd9ehtwMqHsztHm0dMdxoLyCYgk3jlVhhoL/l8EId38M6ODAMgjfVOpDdQzsyHu3vIcOHStOrOCLA19cd3udkxM135tBjbEvABC/eDHnX3gBS3p6qdsgKjYJN6LCS/zxJyKeGmGdwyZkxXKcG11/ccqU7BSeXf8sESkR1K5Wm8/7fI6Hk0fpG3HkZzj5p7p+1N3/qzBFxPnl/aEr4caxKtRQ8P8Kan913am14yHn2jUzN3JHyB1M6jwJgDlhc1hxfMV1t9fpdPiNHk2tjz5C5+RE6t/rOTdkKDlRUaVug6i4JNyICkvJzSXq3Xe5NHmydah33a+/xuh3/YUt03PSefbvZzkSdwQfZx/m9ZlHDbcapW9IZhL8fmUq+ptfhBpFr1NV3knPTflQYXtu8tz2OlQLhPjTsHVWmQ41uOlgRrZUw9L0HdP5+dTPN9zH6+67qPPVYgw+PmQePsyZhwaSvndfmdohKh4JN6JCyk1IIGLESBK++RYAvzFj1KHeLi7X3S8zN5PnNzxvnaTvi9u/IMQrpGyNWT8dUqPAt766oGAFZS0oroC9TpWJdfmFilRQnJ+LF/RXp2Dgn1lw+USZDvd82+d5tKl6m3fqtqmsO7Puhvu4tWtHyA8rcG7cGHNsLOeGDSNx5coytUNULBJuRIWTeeIEZx9+hPQdO9C5uVH700+oMea56w71Bsg2ZzNu0zh2Re3C3eTOvD7zaFa9WdkaE74Zds9Xn981C0zXD1flmVl6bsqFCltQnF/ze6FRX7DkwOpRYM4t9aF0Oh0TOk3gwUYPYlEsTPxnIhsiNtxwP6egIEK+X4rHHXdATg6XXp9C1NvvoOTklLotouKQcCMqlKS1v3J20GB1Ve+gIEK+/x7P22+/4X45lhxe2fwK/174F1ejK3N6z6FVjVZla0xGIqx+Vn3ebhg06FW24zmYWWpuyoVKEW50OrX2zMULLoSq66yVgV6nZ0qXKdxd/27MipmXN7/M1gs3ntNG7+5O7dn/w++F5wFI+O47zg0fTk5MTJnaI8o/CTeiQrBkZxP11nQuvvwyypU1okJ+WIFLkxvXt+SYc3h186tsiNyAk96Jj3t9TPuA9mVv1G+vQPJ58KkHfd8t+/EcTIaClw+VItwAeNVWezMBNr8P50PLdDiD3sD07tO5ve7t6jING8ay5fyWG+6n0+up8eyzBM35DL27Oxl7QjnzwIOk7dpVpvaI8k3CjSj3ci5e5NzjQ0hYuhSA6qOeoc6CBRh9fG64b5Y5i7Ebx/J3xN+Y9Cb+1+t/dK3VteyNOvQjHFwBOj088CU4Vyv7MR1MCorLhwpfUJxfy4egxYOgmOGnkZBdtlW8jXoj79/yPr2Ce5FtyWbsxrGsP7e+WPt69O5NyMofrHU4EU8MJ3b+/OvOoSMqLgk3olxL3byZMw88SOaBA+i9vAia9zn+48ahMxhuuG96TjrPrX+Ofy78g4vBhc96f8atQbeWvVHJF9VhrqAWEAd3KvsxywGLzFBcLlT4guL/uvMj8Kiljp76a2qZD2cymJjZcyZ31L2DXEsuL21+id/P/F6sfZ3r1SNk+TJ1wj+LhcszZ3H+uTHkJiSUuV2ifJFwI8olS3Y2Ue++S+QzozAnJuJy003U+/FHPHr2LNb+qdmpjP57NDsv7cTN6MbcPnPpVqubBg2zqHU2mYlQsw30mFD2Y5YT0nNTPhgNlajnBsDNF+5TF6Rl9wI4+VeZD2nSm3j/1ve5p8E9mBUzE7ZMYPWp1cXaV+/qSs33ZhD45pvoTCZSN2zgzH33y22qSkbCjSh3ssLPcPaRQdZh3j5DhlB36RKcgmoXa//4zHhG/jmSvTF78TB58OUdX9IxsKM2jdsxB8I3gtEFHpgPBpM2xy0HrJP4yVBwh9JXxIUzb6TBberyDKD+cZBS9on1jHoj07tPZ2DjgSgoTPl3CkuOLinWvjqdDp9HHiZk+TKcQkLIjY4mYtgTXP7kE5Tc0o/sEuWHhBtRbiiKQuLKlZx58EGyjh7F4OND0OdzCZz8Gnpn52IdIzI5kiG/DeFQ3CG8nb1Z0HcBrWu01qaBZ/6Bv6apz+94u8JO1nct0nNTPlSaguL/6vMG1GgGaTHwwxNgLvuQ7LxRVI83exyA93a9x6w9s7AoxVtXyqV5c+r9uBKvBx8ARSF27uecGzKU7PMXytw24VgSbkS5kHv5MueffY5Lr09BycjArWsX6q1ejUev4g+vPhx7mMd/f9y6pMI3/b+hefXm2jQw6QKsHK4WRrZ6BDqO0Oa45YgMBS8fKm24MbnCI9+BsydEbC/T4pr56XQ6Xu34KmPbjQVg8eHFTPpnEjnFDE96d3dqvfMOtWZ+hL5aNTL27ePMPfeQ8MMPUmxcgUm4EQ6X/PvvhN89gNSNG9GZTPi//BJ1Fi7EFOBf7GP8c/4fhv8xnPjMeJr5NuO7O7+jnlc9bRqYmwU/DIO0yxDQEu6eXSHXjroRmcSvfMi7LZhb2cINgF9DuP/KIpg758GB668XVVw6nY4RLUfwzs3vYNQZ+e3Mb4xeP5qU7JRiH8Prrruot3oVru3bY0lPJ2rKVCKfeYacaJkTpyKScCMcJjchgQvjx3PhxfGYk5Jwbt6MkB9XUn3EiBvONpzfiuMreH7D82TkZtCtVjcW91uMn+v115cqkXWT4PxudUKyR74BJzftjl2OmGWem3Ihr6DYUhnDDUDTO+HWV9Tna16AqIOaHfqeBvcwp/cc3Ixu7Ly0k2HrhnEx9WKx93cKCqLuN1/j/+qr6JycSNvyD+H33EPSL79IL04FI+FG2J2iKCStWUP4nXeR/NvvYDDg9+yz1Fu2DJfGxa9jybHk8PaOt5m+YzpmxcyA+gP47LbPcDe5a9fYfUtgz0JABw8sUNePqqQssrZUuVDphoIXpeckaNAbcjNg+eOQod1Q7G61r/6BczLhJIPWDmJP1J5i768zGKj+5HDq/fQjLjfdhCUpiYuvvKr24lyQWpyKQsKNsKvsyEgiR4zk4qsTMCck4NyoESHLvqfGC8+jc3Iq9nHiM+N5+s+nWX58OTp0jG03lndufgeTlqOXzvwDa8epz3tOhMZ3aHfsckgKisuHSjWJ37XoDfDgAvCuAwlnYfkQ9favRppXb873d31PM99mJGQlMPLPkaw4XrJbYM4NG6qfTWNfQGcykbblH07fPYC4r76SEVUVgIQbYRdKdjZxCxcSPuAe0v79F52TEzXGjaPeTz/i2rJliY51PP44g9cOZk/0HtxN7nxy2yeMaDkCnZY9DtGHYdljYM6GZgPg1le1O3Y5ZZGC4nJBXxXCDajz3wxaCk4ecPYfWPXM1f8INRDoHsjX/b+mf0h/cpVcpu+Yzts73ibHUvxRWjqTCb/Ro6n382pcO7RHycgg5r33OfvIIDIOHdasrUJ7Em6EzaX++y/h991PzIcfoWRm4ta5M/XX/IzfqGfQmYrf06IoCqtOruLx3x7nYtpF6njUYcmdS+gZ3FPbBidGwncPQlYS1Omm3o4qQQ1QRZXXc2OUcONQVaLnJk9gSxj0HehNcHgV/PEaaHg7ztXoyvu3vs/YdmPRoWP58eUMXzecS6mXSnQc5/r1qfvNNwS+9SZ6Dw8yDx/m7MCBXJo6TWY3Lqcq/ye2cJjs8xc4//zzRD41guzwcAy+vtR85x3qfLUYp5CQEh0rPSed17a+xtRtU8k0Z9K9VneW3rWUBt4NtG10erwabFIuqXNyDF4KJhdtz1EOKYpirbnRS7hxqCrTc5Onfk+4f576fOfnsO0TTQ+fN5Lq09s+xcPkwf7L+3nol4fYGLGxZMfR6/F5+GHq/7oWz7vvBkUhccUKTvfrT/ySJXKrqpyRcCM0Z05JIeZ/swm/6y5S/vobDAZ8hw2lwbrf8X7wgRLfPjoef5xH1j7C2vC16HV6Xmj7AnP7zMXL2UvbhudkwPeDIPY4eNaGx1eC640X56wM8v8ilYJix6rUQ8GvpeVDcMc76vO/psL+ZZqfokdwD1YMWEGL6i1Izk7mhY0v8MHuD4o9H04ek78/tT/6kLrffYtzkyZYkpKInv42Zx58iNSt/2reblE6Em6EZpTsbOK//Y7Td/Ql7osvULKycOvcmXqrfiJg0iQMnp4lOp5FsbDk6BIe++0xziafxd/Nn0V9FzGy1Uj0Oo3/081OU4NN5E51yPfjP4JXkLbnKMfyj8wxGCTcOFLebSlLZR4tVZRuY6DrGPX56mfh4ErNTxHkEcQ3/b9hSPMhAHx75FuG/D6E8KTwEh/LrUMH6v24koApr6P38iLr+HEiR4wg4smnyDxyROumixKScCPKTLFYSP79d07fPYDod97BnJCAU716BH32KXW+Wlyi4d15LqReYOSfI3lv13tkmbO4ufbNrBywkvYB7bX/AbJSYMlACN8ETtVg8HLwb6b9ecqx/HWc0nPjWFXutlR+t0+H1o+qM4H/NBLClmp+CpPBxKsdX+WTXp/g6eTJ4bjDPPzLw3x9+GvMFnOJjqUzGvF97DEarPsd32HDwGQibds2zjzwIBdefZXs8+c1b78oHgk3otQUi4XkdX9w5t77uPDieHIiIjD4+RH4xhvU/2UNHn36lPgWlKIorDyxkgd+foBdUbtwNbryWufXmNN7Dj4uNrhFlJEI394P5/5Vp4UfsgrqdtX+POVcbr50I6OlHKtKFRT/l14P986BdsNAsag9OHsW2+RUver04sd7fqR7re5kmbP4aM9HPPnHk0QmR5b4WEYfHwImTaTB77/hedddACSv+YXT/fpzacoUWavKASTciBIrEGrGjSPr5En0Hh74PT+Ghn+sw2fQI+iMxhIfNzI5ktF/j+bN7W+SnptOO/92rBywksFNB2t/GwrU4uFv7r0y+7A3DFsDwZ20P08FUKDnRsKNQ1XpnhtQA86Aj6HTM4CizjW180ubnCrQPZDP+3zOtK7TcDO6sTdmLw/+8iDfHvmWXEvJC4SdgoKoPfMjQlauxL17d8jNJfGHlZzu149LU6ZKyLGjkv8GElWWJTub5DVriFu0mOxw9R61vlo1fIcNw3fY0BLX1OTJNmez6NAiFhxcQJY5Cye9Ey+0e4HHmz2OQW/Q8ke4KjEClj4CMUfAzQ+GrlaHpVZRBXpu5LaUQ1Xpnps8Oh30fx+MTrDtU/j9FciIhx4TNF/XTafT8VDjh+haqytT/p3C7qjdfLD7A34+9TOvd3mdNv5tSnxM1xY3UWfhAtL37iP2s89I27aNxB9+IPGnn/C8806qP/UkLk2bavpziIIk3IgbMicnk7BsOfHffoP5cixwJdQMHYLvsGEYvEo/amnbxW28u/NdziWfA6Bzzc683vl1QrxCtGh60SJ2wvLH1IUwqwXC0J/Bv2p/0OQVFOt0MhTc0arE8gvFodOpNThGV9jyAWyaAZePw31z1RXGNVa7Wm0W3LGAn07+xP9C/8fxhOMM+X0IDzZ6kHHtxuHt4l3iY7q1a0udRQtJ37uX2M/mkLZtG8m//ELyL7/g3r071Z96EreuXbWdgFQAEm7EdWQeP0HC0qXqonHp6QAYAwLwHTYM74cHYqhWrdTHPpN0htmhs9kQuQEAP1c/Xu34Kv1C+tn2H3rY9/DLC+rMw4EtYfCyKjUq6lqsi2bKh6zD5d0WrFJDwa9Fp4PbJoN3MKx9EQ7/BAlnYND34FlT89PpdXoeavwQt9W5jVl7ZvHz6Z/58eSP/B3xN8+0eoZBTQaVaokXt3btqLNoIRmHDhO/aCHJ6/4g7d9/Sfv3X5wbN8bn0cF4DRiA3l3DdfGqOJ1SxZY6TU5OxsvLi6SkJDxLeRulMlOys0lZv56EJUtJ33N1sTnnRo3wfepJvO68s0RrQP1XXEYcn+//nJUnVmJWzOh1egY1GcSYtmPwcPLQ4kcomsUC69+Ef2er3zcbAPd/AU7yYQJwPiGdm9/fiLNRz/G3+zu6OVXat9vPMuXnw/RvEcjnj9tgdGBFdXarugZVRjx41FSXbqjdzqanDI0O5e0db3Mq8RQAQdWCGNt+LH3r9i3TH2HZkZHEf/U1iT/+iJKZCai94V733YfPo4Nxrl95F+gti5L8/pZwIwDIPHqUxFWrSP5lLea86cQNBjx698bn0Udx69ypTP+YU7JTWHJ0CYsPLSY9V+0F6hnUk3Htx2k/y/B/JV2A1aPgzBb1+1tehl6Tq8SSCsV1Li6NHh9uws3JwJG3+jm6OVXakp3nmLzqELc3D2D+0A6Obk75En9GnY/q8jEwOMPtb6qFxzb8t5xryWX1qdXMCZtDbIZ6W76lX0teaPcCnQM7l+lz0ZyUROKqVSR8/z055yKsr7u2aYPXA/fj2b8/Bg8b/tFXwUi4uQ4JN1flREeTsm4diat/JuvoUevrxho18B74EN4PP4wpMLBM50jKSmLJ0SV8d/Q7UrJTALip+k281OElOgZ2LNOxi+XwKvhlHGQmgskNBnwCrQba/rwVTPjlVG6buRkPFyMH3+jr6OZUact3RzDhx4P0burPwifs8G+koslMhp+ehhO/q9836K3W4XiU7bPqRtJz0vn6yNcsPrSYjNwMANrUaMMzrZ+he63uZQo5isVC2r/bSFi6lNTNm63DF3UuLnjcfjteA+7GvWvXEq3FVxlJuLmOqh5uci9fJvmPP0n+/XcyQkOtr+tMJqr17o33/ffh3r17qYZy5xeXEceSo0tYemwpaTlpANT3qs+o1qPoG9LXNkO788tKgd8nQNgS9fta7eDBBVDdxr1EFdSpmBT6zNqCt5uJsKl3OLo5VdoPeyJ5ZeUBejapwVfDq+bUBDekKLB7Afz5OuRmglt1uOczaHqnzU8dmxHL/APz+fHkj2SZswBoUb0FI1uNpGdwzzJ/tuXExJD8yy8k/rSK7NOnra8bvLyodnsfPPv3x71z5zJ/RldEEm6uo6qFG0VRyDp5ktRNm0nduJGMsLACq+66tmuH51134nnnnRh9yj5J3rH4Y3x35Dt+O/MbORZ1zZbGPo15utXT3F73dtuHGkWB47/B7xMhKQLQwS0vQc+JUIpCwKriWFQy/Wb/g181J/a8frujm1Olrdp3nheX7+eWRn58+1RnRzenfIs5Bj+OgOiD6vctH4Y7ptu8Fwfgcvplvjr8FSuOryDTrNbN1PGow6PNHuW+hvfhbipbPZ+iKGQePEjS6p9J/vNPzLGx1vcM3t5U63Er1Xr2xP3mm6vMrSsJN9dRFcKNOTmZ9F27SNu2ndTNm8m5UHDiKNfWrfHo3w/Pfv3KfNsJIMeSw+bIzSw9tpTdUbutr7f0a8lTLZ+iV3Av24cagLjTam/Nqb/U773qqKsNh3S3/bkruMMXk7jrk634eziza3IfRzenSvs57AJjl4XRrUF1lo7s4ujmlH+5WbD+Ldg+B1DAyQN6TYJOT9vlD5q4jDi+OfINP5z4wXrrvZqpGg80eoCHmzxMXc+6ZT6HYjaTvnsPyet+J+WPP6/WRQIYjbh16EC1m7vj1rUrLs2aoauk9YQSbq6jMoYbc1ISGfv3k757D2k7dpB5+HCBKWd1Tk64de2CR69eVOvRA1NNbYZQhieGs+rUKtacXkN8ZjwABp2BO+rewWPNH6N1jdaanOeGslJh6/9g2yfqEG+9Cbo9D7e+LKOhiunA+UTu+exfanq5sH1Sb0c3p0pbe+AiY5buo1M9X1Y8U/WWAim1C3vht5fhwpXb7TWaqRMB1u9hl9On56Sz5vQalhxdwtnks9bX2/m344FGD3B73dtxM7mV+TxKbi4Z+/aRsmkTqRs3WSdUzWPw8sKtc2fcOnfCrV07nBs3Rmew0WSodibh5joqerixZGeTdfIkWUePqoFm3z6yT50utJ1T/fq4d+mC+83dce/SBb1b2f9RgXq/+e9zf7M2fC37L++3vl7dpTr3NbyPQU0HEehu+y5hQK2r2TVfncE0Qw1XNOgN/T8Av4b2aUMlsS8igfvnbiPIx5WtE25zdHOqtHWHLjHqu710qOvDytHdHN2cisVigX3fwt9vXP1MCLlFvS0dcrN9mqBY2HZxG0uPLuXfi/9iUdQ/NN1N7vQL6Ue/ev3oENABo16bmpnsc+dI3byZtO07SN+1C0taWoH39W5uuLZpjWubtri0aIFL82YYAwIq5MSBEm6uo6KEG8ViIefCBbLDw8k6HU7WqVNkHj1K1qlTkJNTaHunkBBc27bFrXMn3Lt2xRQQoFlb4jPjWR+xnj/O/MHu6N3Wf6wGnYFbg27l/ob3c3PQzZj0dqppKSrU+DaAPm+o89dUwH+0jrbnbDwPzdtO3epubH6ll6ObU6X9eTiKp78NpU2wN6ufk1uqpZIeDxvfhdCv4Ertn71DDkB0WjRrTq9h1alVRKZcXZDT18WX2+veTt+QvrTzb6fZMjNKbi4ZBw+SvmMH6XtCyQgLKxR2AAw+Prg0a4Zzs6Y412+Ac4P6ONWvX+oldOxFws11lJdwY8nKwhwbS25sLDkxMeRcuEDOhYvknD+vhppz51CysorcV+/lhUuzZri2uAnXdu1wbdMGo6+vZm1TFIUTCSfYcn4LW85vYf/l/Shc/c+kpV9L+ob05a76d+Hn6qfZeW/o8gnYs1CdZTgrSX3NtwH0eBVaPASGqjd6QCs7w+N45Msd1K/hzoaXejq6OVXahmPRPPnVHloFebFmjP1+EVdKiZHqLeu931wNObXaQscRcNMD4KRNj/aNWBQLodGh/HbmN/4+9zeJWYnW97ydvbm59s3cGnQr3Wp1w8u59MvZ/JdiNpN16jQZ+/aRsW8fmUeOkBUeDmZzkdsbavjhVKcuptq1MNWujVPt2hhr1sRYowbGGjUweHk5tJ5Hws112Crc5M0ZY8nKRsnKQsnOwpKVhZKRgTklFUtyMuaUFCzJyeQmJGBJTr7hMXUmE04hITg1aIBz/fo4N2uKa/PmGGvV0rxLMTotml1Ru9gdtZvtl7YTlRZV4P1mvs3oG9KXviF9CfKw43IFuVlw/Hc11ORNwgcSajS27XQsj87fSSP/avw13j41CqJom47H8MTi3dxUy5NfX7jF0c2pHPJCzr5v1bo8ABcvaPM4tH8CajS2W1NyLDnsurSLdWfXsT5ivbUIGdTe8NY1WtO5Zmc6BXaiVY1WOBlKPyN8USyZmWSdOEHmkaNknThB1plwsk+HkxsTc+OdjUaMvr4YvL3Re3hg8PBA7+mBoVo1dE7O6Jyd0Tk7oXd2xhgYiNddd2nadgk312GrcJO+dx/nHn20RPvoTCYMNfww1qiBU+3amPI9nOrUwRQUZJNCMEVROJt8lv2X97P/8n72RO0pUAAH4GxwpkvNLtwadCu3Bt1qvzoagNxsOLMZDv0Ex3692kuj00PjftDxKah/m8wwrKF/Tl5myMJdNA30YN24Wx3dnCpt68lYHl+4U/6/sIW0WDXg7FkEiVdnBCagJbS4H266H3ztt/RBriWXsJgway/56aSC9ZMuBhfa+LehnX87WtdoTYsaLfB0ss0dB3NqKtlnzly9e3DhAjnnL5AbdYnc2LiCI7SKwbV1a0KWL9O0jSX5/S1/8mrEWMMPz7vvtqbW/CnW4OGJwdMDvYcnBo9qGHx9Mfr5off0tHlRl6IonE89z7H4YxyLP8aRuCMcuHyA5OyCPUc6dDSr3ozOgZ3pGNiRDoEdcDVqv/LuNaXFQvgmOLVenacmM/Hqe9UCoe2Vv7C8g+3XpirEunCmrAjucHmZ3SwLZ2rP3Q9ufhG6vaB+1uxZCKf+VufJiT6oDimv2Ub9I6rBbVC7vU17ho16Ix0CO9AhsAPjO4wnMiWSHZd2sOvSLnZF7SI+M54dl3aw49IOQP2cru9Vn5Y1WtLUtylNfZvSxKcJ1ZxKv4hxHkO1ari2bIlry5ZFvq/k5JAbH0/u5VgsyUmYk1MwpyRjSU7BkpaKkp199c5FVhamYMd+Vku40YhTcDC1P/rQYec3W8xEpUcRnhjOmaQzhCepX08knCA1J7XQ9s4GZ26qfhOtarSirX9b2ge01/Re7w2lxcGFPRCxHU5vgEv7C77v7g/N74UWD0BwF+mlsbG8X6RGCTcOZ7zy37qEGxvSG6DxHeojPR6O/qIu1XJmC1wKUx+b3wNnL6h/K9TrAUEdIaCFTcNOsEcwwR7BDGw8EEVROJ14mt3Ru9Ve9pj9nE89z+mk04V6eII9gmng3YD6XvWtjxCvEE0XI9aZTJgCAjQdrGJLEm4qiGxzNjHpMUSnRxOdFk1UehQXUy9yPuU851PPcyH1ArmW3CL3NelNNPRuSLPqzWji04RWNVrRxKcJJnvN2JseD9GHIeaIOhfF+d0QX3j4OgEtoUEvaHQ71O2ufgAJu8j7RaqXcONwhryem6pVMeA4br7Qfpj6SItVe45Pb4DTG9Ue5KO/qA9Q16er1RaCOkBgK/BvDn6NbDJZoE6no6FPQxr6NGRw08GAOmHggcsHOBx3mOPxxzmWcIyotCgiUyKJTIlkU+SmAsfwcvYiqFoQQR5BBFULoqZ7TQLcAwhwCyDAPQAfZ58KOSS8OMpFuJkzZw4ffvghUVFRtG7dmk8//ZROna69psoPP/zAlClTOHv2LI0aNeL999/nzjttv6ZIWSiKQrYlmyxzFuk56aTnppOek05aThqpOakkZyWTnJ1MUlYSydnJJGQmEJcZR3xmPHEZcYVuIxXFpDdR17Mu9bzqWdN7A+8G1Peub/th2tlpkHAOEs6oK/cmnFFnDI45AqnRRe/j11j9a6heD6jfEzwqxl8ElZH1tlQl/aCrSPRX/j/INUu4sTt3P2g3VH1YzGoPzqkNag/z+T1q/d+5f9VHHr1JDTg1moJvPfCpd/WrR6Cmf6RVd61Orzq96FXn6nQNiZmJHE84TnhSeIGe+8sZl0nKSiIpK4nDcYeLPJ5Rb8TXxZfqLtWp7lodXxdfvJ298XTyxMvZC08nTzycPHA3ueNucsfN6IabyQ1XoytOBifN5uqxBYe3bPny5YwfP5558+bRuXNnZs+eTd++fTl+/Dj+/v6Ftt+2bRuDBw9mxowZ3H333SxdupT77ruPvXv30qJFCwf8BKqjcUd5dcur5FpyMStmzBYzuUouuZZcssxZ1gXWysJJ71Qgdddyr0WwR7A1lfu7+Ws2XwIWM2QmqX+5ZCZBepx6Kyk9Vv3rJi0Gki9ByiVIvqBucz3edSHgJghsCUGdoHY79S8mUS7k9RJIzY3j5d2WskjPjWPpDWrNTe326vcWC8SdVHuez+9R/3CLOQpZyVeeHyl8DJ1BDTietdRHtUA1QLn7gduVry7e4Oqtjt4yuZV4ni5vF2861+xM55oF1yFLz0nnfOp5zqecJzIlkvMp59We/yu9/3GZceRacolJjyEmvRgjpYpg0BlwNjhbg45BZ7B+bV69OR/2cFyphsNHS3Xu3JmOHTvy2WefAWCxWAgODub5559n4sSJhbZ/5JFHSEtLY+3atdbXunTpQps2bZg3b94Nz2er0VKh53fyxPoRxdrWoDPgZnDBRe+Mm8EFV4MLHkY3PAzueBjdqGZww9tYDW+TOz5GD3yM7vgY3fHSu6JDAcUCihmdxQJKLlhy0SkWMOegs+SCRf2qs+Sor5mz0Jmz0ZmzwZyFLjcTXW7Gla+Z6HPS0eWkostOR5+Thi4nFX124TqdG7E4eZLrXZdcr5CrD7+m5Pg2RtGg4E3Yzt9HonnjlyOynlE5kLfOV3V3J34eI5P4lWuKgiHlPKbYYxgTT2NMPIsx6SzGpHMYkiPRKUXPJ3PNw+mNWJw8UEzuKE7VsJjc1ecmVxSjC4rBRX1ucEYxOKEYnMHgdOW5CfRGFL0JDCb1q86AojeAzgh6PYrOoAYunZ5sLMTnpJJgTiMxN5X43FQSc1NJyU0jJTedFHM6KbnppOamk2HJItOcSbo5k0xL8f5Qb+HVlO/v+6E0V/WaKsxoqezsbEJDQ5k0aZL1Nb1eT58+fdi+fXuR+2zfvp3x48cXeK1v376sXr26yO2zsrLIyjcZXnIx5pcpDd2FWL66GI0BBZMCBhQMV766KApOFvWrs6JgBCrK38dpijPJuJOouBOneBKPJ3GK+ojGhyjFl0uKL9GKLymZrpD8358sCdhd1KFFOSQ9N46X13MTl5bNze9vdHBrRPEYgSZXHioDZvxIIlAXT6Aunpq6eGroEvElBT9dMr66ZHxJxlOXjhdpGHUWdJZcDJkJkFmyYdelVbsU+1iALJ2ObJ2OrHwPM2DWgRkdZh1ERblo3NqScWi4iY2NxWw2E/Cf6uuAgACOHTtW5D5RUVFFbh8VFVXk9jNmzODNN9/UpsHX4W50JzDz6uXM3x2mXIkyOejIufJe3muK2heDgt76ugKY0QM6LOiwXHnPgh4zepQrX68+DOSix6LoycFILkZyMJCLgRyMZGMiO9/XLMWJTJzIxJlMnMjAmTRcSFdcSMeFNFxIxp1U3Mi51n8iuv98BZy1uJDCYUwGPf1a2HE+I1Gken7utA7y4lhUyo03FuWYnkT8SMSPQr/NFAr+kkDBjSw8SMOTdFzJxE2XiTtZuJOBM9m4koUzObjosnEhGxO5OJGD05WvRswYMWMiFyO5mDCj11kwYMFw5T31t4f6G8WABT0WdFd+++T9ptFd+e2U/zeSDq58vfKbS1H/xwSYAI8rP0z+P410zo7trXd4zY2tTZo0qUBPT3JyMsE2GH/fpH1PaH9R8+MKIaoWJ6Oen2XZBVHBtXbw+R0abvz8/DAYDERHFxxNEx0dTWBg0X9BBgYGlmh7Z2dnnJ2lT0EIIYSoKhw6M5qTkxPt27dn/fr11tcsFgvr16+na9euRe7TtWvXAtsD/PXXX9fcXgghhBBVi8NvS40fP55hw4bRoUMHOnXqxOzZs0lLS2P48OEADB06lNq1azNjxgwAxo4dS48ePZg5cyZ33XUXy5YtY8+ePXz55ZeO/DGEEEIIUU44PNw88sgjXL58malTpxIVFUWbNm1Yt26dtWg4IiICfb6p97t168bSpUt5/fXXee2112jUqBGrV6926Bw3QgghhCg/HD7Pjb3Zap4bIYQQQthOSX5/y2qEQgghhKhUJNwIIYQQolKRcCOEEEKISkXCjRBCCCEqFQk3QgghhKhUJNwIIYQQolKRcCOEEEKISkXCjRBCCCEqFQk3QgghhKhUHL78gr3lTcicnJzs4JYIIYQQorjyfm8XZ2GFKhduUlJSAAgODnZwS4QQQghRUikpKXh5eV13myq3tpTFYuHixYt4eHig0+k0OWZycjLBwcFERkbKelU2JtfafuRa249ca/uRa20/Wl9rRVFISUmhVq1aBRbULkqV67nR6/UEBQXZ5Nienp7yj8VO5Frbj1xr+5FrbT9yre1Hy2t9ox6bPFJQLIQQQohKRcKNEEIIISoVCTcacHZ2Ztq0aTg7Ozu6KZWeXGv7kWttP3Kt7Ueutf048lpXuYJiIYQQQlRu0nMjhBBCiEpFwo0QQgghKhUJN0IIIYSoVCTcCCGEEKJSkXBTRnPmzCEkJAQXFxc6d+7Mrl27HN2kCm/GjBl07NgRDw8P/P39ue+++zh+/HiBbTIzM3nuueeoXr061apV48EHHyQ6OtpBLa483nvvPXQ6HePGjbO+JtdaOxcuXODxxx+nevXquLq60rJlS/bs2WN9X1EUpk6dSs2aNXF1daVPnz6cPHnSgS2umMxmM1OmTKFevXq4urrSoEEDpk+fXmBNIrnWpbNlyxYGDBhArVq10Ol0rF69usD7xbmu8fHxPPbYY3h6euLt7c1TTz1Famqqtg1VRKktW7ZMcXJyUhYtWqQcPnxYGTlypOLt7a1ER0c7umkVWt++fZXFixcrhw4dUsLCwpQ777xTqVOnjpKammrdZtSoUUpwcLCyfv16Zc+ePUqXLl2Ubt26ObDVFd+uXbuUkJAQpVWrVsrYsWOtr8u11kZ8fLxSt25d5YknnlB27typhIeHK3/88Ydy6tQp6zbvvfee4uXlpaxevVrZv3+/cs899yj16tVTMjIyHNjyiuedd95Rqlevrqxdu1Y5c+aM8sMPPyjVqlVTPv74Y+s2cq1L57ffflMmT56s/PTTTwqgrFq1qsD7xbmu/fr1U1q3bq3s2LFD+eeff5SGDRsqgwcP1rSdEm7KoFOnTspzzz1n/d5sNiu1atVSZsyY4cBWVT4xMTEKoGzevFlRFEVJTExUTCaT8sMPP1i3OXr0qAIo27dvd1QzK7SUlBSlUaNGyl9//aX06NHDGm7kWmtnwoQJys0333zN9y0WixIYGKh8+OGH1tcSExMVZ2dn5fvvv7dHEyuNu+66S3nyyScLvPbAAw8ojz32mKIocq218t9wU5zreuTIEQVQdu/ebd3m999/V3Q6nXLhwgXN2ia3pUopOzub0NBQ+vTpY31Nr9fTp08ftm/f7sCWVT5JSUkA+Pr6AhAaGkpOTk6Ba9+0aVPq1Kkj176UnnvuOe66664C1xTkWmtpzZo1dOjQgYEDB+Lv70/btm2ZP3++9f0zZ84QFRVV4Fp7eXnRuXNnudYl1K1bN9avX8+JEycA2L9/P1u3bqV///6AXGtbKc513b59O97e3nTo0MG6TZ8+fdDr9ezcuVOztlS5hTO1Ehsbi9lsJiAgoMDrAQEBHDt2zEGtqnwsFgvjxo2je/futGjRAoCoqCicnJzw9vYusG1AQABRUVEOaGXFtmzZMvbu3cvu3bsLvSfXWjvh4eF8/vnnjB8/ntdee43du3fzwgsv4OTkxLBhw6zXs6jPFLnWJTNx4kSSk5Np2rQpBoMBs9nMO++8w2OPPQYg19pGinNdo6Ki8Pf3L/C+0WjE19dX02sv4UaUa8899xyHDh1i69atjm5KpRQZGcnYsWP566+/cHFxcXRzKjWLxUKHDh149913AWjbti2HDh1i3rx5DBs2zMGtq1xWrFjBkiVLWLp0KTfddBNhYWGMGzeOWrVqybWuIuS2VCn5+flhMBgKjRqJjo4mMDDQQa2qXMaMGcPatWvZuHEjQUFB1tcDAwPJzs4mMTGxwPZy7UsuNDSUmJgY2rVrh9FoxGg0snnzZj755BOMRiMBAQFyrTVSs2ZNmjdvXuC1Zs2aERERAWC9nvKZUnavvPIKEydOZNCgQbRs2ZIhQ4bw4osvMmPGDECuta0U57oGBgYSExNT4P3c3Fzi4+M1vfYSbkrJycmJ9u3bs379eutrFouF9evX07VrVwe2rOJTFIUxY8awatUqNmzYQL169Qq83759e0wmU4Frf/z4cSIiIuTal1Dv3r05ePAgYWFh1keHDh147LHHrM/lWmuje/fuhaY0OHHiBHXr1gWgXr16BAYGFrjWycnJ7Ny5U651CaWnp6PXF/z1ZjAYsFgsgFxrWynOde3atSuJiYmEhoZat9mwYQMWi4XOnTtr1xjNSpOroGXLlinOzs7KV199pRw5ckR5+umnFW9vbyUqKsrRTavQRo8erXh5eSmbNm1SLl26ZH2kp6dbtxk1apRSp04dZcOGDcqePXuUrl27Kl27dnVgqyuP/KOlFEWutVZ27dqlGI1G5Z133lFOnjypLFmyRHFzc1O+++476zbvvfee4u3trfz888/KgQMHlHvvvVeGJ5fCsGHDlNq1a1uHgv/000+Kn5+f8uqrr1q3kWtdOikpKcq+ffuUffv2KYAya9YsZd++fcq5c+cURSnede3Xr5/Stm1bZefOncrWrVuVRo0ayVDw8ubTTz9V6tSpozg5OSmdOnVSduzY4egmVXhAkY/Fixdbt8nIyFCeffZZxcfHR3Fzc1Puv/9+5dKlS45rdCXy33Aj11o7v/zyi9KiRQvF2dlZadq0qfLll18WeN9isShTpkxRAgICFGdnZ6V3797K8ePHHdTaiis5OVkZO3asUqdOHcXFxUWpX7++MnnyZCUrK8u6jVzr0tm4cWORn8/Dhg1TFKV41zUuLk4ZPHiwUq1aNcXT01MZPny4kpKSomk7dYqSb8pGIYQQQogKTmpuhBBCCFGpSLgRQgghRKUi4UYIIYQQlYqEGyGEEEJUKhJuhBBCCFGpSLgRQgghRKUi4UYIIYQQlYqEGyGEEEJUKhJuhBAVysKFC7njjjtsfp5169bRpk0b63pEQoiKQ8KNEKLCyMzMZMqUKUybNs3m5+rXrx8mk4klS5bY/FxCCG1JuBFCVBgrV67E09OT7t272+V8TzzxBJ988oldziWE0I6EGyGE3V2+fJnAwEDeffdd62vbtm3DycmJ9evXX3O/ZcuWMWDAgAKv9ezZk3HjxhV47b777uOJJ56wfh8SEsLbb7/N0KFDqVatGnXr1mXNmjVcvnyZe++9l2rVqtGqVSv27NlT4DgDBgxgz549nD59uvQ/rBDC7iTcCCHsrkaNGixatIg33niDPXv2kJKSwpAhQxgzZgy9e/e+5n5bt26lQ4cOpTrn//73P7p3786+ff9v7+5VWgnCMI4/K2j8wkYkkEIQMUJAUcHCKxBMHSwSU/iB3kEQwUZEUItgES9gEQshYGMtSQoR8XSxWkkKBUXUtVAR9ZzigKghCknMmuX/q5Z3Z+ad8mGW3f2jYDCoiYkJRaNRRSIRHR8fq7u7W9FoVO//JdzZ2Smv16t0Ol1STwDOINwAcMTY2JhmZmYUDoc1NzenlpYWraysFB1/e3sr27bl8/lK7jc7O6uenh4tLi7q7u5Ow8PDCoVC8vv9isViOjk50cXFxYd5Pp9P+Xy+pJ4AnEG4AeCY9fV1PT8/a2dnR1tbW/J4PEXHPjw8SJIaGxtL6tXf3/927fV6JUl9fX0FtcvLyw/zmpqadH9/X1JPAM4g3ABwjGVZOj8/1+vrq3K53Jdj29vbZRiGbm5uvl335eWloFZfX/92bRhG0drnV7+vr6/V0dHxbU8AvwfhBoAjnp6eFIlEND4+rqWlJU1PTxecmrzX0NCgQCCgbDZbcO/zo6TT09OK7PHx8VGWZWlwcLAi6wGoDsINAEcsLCzItm1tbGwoFovJ7/drcnLyyzmjo6PKZDIF9d3dXSWTSVmWpeXlZWWzWeXzeZ2dnZW1x4ODA3k8Ho2MjJS1DoDqItwAqLr9/X3F43GZpqm2tjbV1dXJNE2l02ltbm4WnTc1NaW9vT3Ztv2hHgwGtbq6qkAgoFQqpUQiocPDQ5mmWdY+t7e3FQ6H1dzcXNY6AKrL+Pv+vUcA+OVCoZCGhoY0Pz8v6f93bgYGBhSPxyva5+rqSr29vTo6OlJXV1dF1wbwszi5AVBT1tbW1Nra+uN9crmcEokEwQaoQZzcAKhpP3VyA6B2EW4AAICr8FgKAAC4CuEGAAC4CuEGAAC4CuEGAAC4CuEGAAC4CuEGAAC4CuEGAAC4CuEGAAC4yj8olOzNKUanWAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "n.finitialize(-65)\n", "for tstop in [0, 50, 100, 150]:\n", " n.continuerun(tstop)\n", " plot_it()\n", "pyplot.xlabel(\"x (µm)\")\n", "pyplot.ylabel(\"[IP3]\")" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:57.358504Z", "iopub.status.busy": "2025-08-18T03:36:57.356505Z", "iopub.status.idle": "2025-08-18T03:36:57.373473Z", "shell.execute_reply": "2025-08-18T03:36:57.373080Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "crossed 100 nM at t = 2.8250000000000064 ms\n" ] } ], "source": [ "nM = 1e-6 # relative to NEURON's default mM\n", "\n", "n.finitialize(-65)\n", "while dend(0.7).ip3i < 100 * nM:\n", " n.fadvance()\n", "\n", "print(\"crossed 100 nM at t = {} ms\".format(n.t))" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "execution": { "iopub.execute_input": "2025-08-18T03:36:57.376864Z", "iopub.status.busy": "2025-08-18T03:36:57.376704Z", "iopub.status.idle": "2025-08-18T03:36:57.754259Z", "shell.execute_reply": "2025-08-18T03:36:57.753889Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "peak ip3 = 0.25437182627374194\n", "final ip3 = 0.2103198308254167\n" ] } ], "source": [ "ip3_vec = n.Vector()\n", "ip3_vec.record(dend(0.7)._ref_ip3i)\n", "t_vec = n.Vector()\n", "t_vec.record(n._ref_t)\n", "\n", "n.finitialize(-65)\n", "n.continuerun(1000)\n", "max_ip3 = max(ip3_vec)\n", "print(\"peak ip3 =\", max_ip3)\n", "print(\"final ip3 =\", dend(0.7).ip3i)" ] } ], "metadata": { "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.12.10" } }, "nbformat": 4, "nbformat_minor": 2 }