{ "cells": [ { "cell_type": "raw", "id": "beba14b5", "metadata": {}, "source": [ "---\n", "aliases:\n", "- /2026/07/pysm-3.4.6-release\n", "categories:\n", "- python\n", "- pysm\n", "- cosmology\n", "date: '2026-07-20'\n", "description: A practical introduction to the new FLAMINGO kinetic SZ presets and documented c1 cosmology in PySM 3.4.6.\n", "title: 'PySM 3.4.6 released: FLAMINGO kSZ models and c1 documentation'\n", "---" ] }, { "cell_type": "markdown", "id": "16c0d5e9", "metadata": {}, "source": [ "[PySM 3.4.6](https://github.com/galsci/pysm/releases/tag/3.4.6) is out! Install it from [PyPI](https://pypi.org/project/pysm3/3.4.6/) or conda-forge:\n", "\n", " pip install --upgrade pysm3==3.4.6\n", "\n", "This release adds two lensed kinetic Sunyaev–Zeldovich (kSZ) presets based on the [FLAMINGO simulations](https://doi.org/10.1093/mnras/stag625), and documents the cosmological parameters used by the long-standing `c1` CMB model. This notebook shows how to use the new presets with the standard PySM API." ] }, { "cell_type": "code", "execution_count": 1, "id": "fa1f9a40", "metadata": { "execution": { "iopub.execute_input": "2026-07-21T02:44:16.115143Z", "iopub.status.busy": "2026-07-21T02:44:16.114789Z", "iopub.status.idle": "2026-07-21T02:44:16.121475Z", "shell.execute_reply": "2026-07-21T02:44:16.120041Z" } }, "outputs": [], "source": [ "import healpy as hp\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pysm3\n", "from IPython.display import Markdown, display\n", "from pysm3 import units as u" ] }, { "cell_type": "markdown", "id": "350605f6", "metadata": {}, "source": [ "## New FLAMINGO kSZ presets\n", "\n", "PySM 3.4.6 adds two models that can be selected through `preset_strings`:\n", "\n", "- `ksz5`: the fiducial FLAMINGO `L1_m9` feedback model\n", "- `ksz6`: the strongest-feedback `fgas-8sigma` model\n", "\n", "Both are lensed, full-sky HEALPix templates at native NSIDE 4096. The following summary identifies the model and template used by each preset." ] }, { "cell_type": "code", "execution_count": 2, "id": "e7451ad2", "metadata": { "execution": { "iopub.execute_input": "2026-07-21T02:44:16.125021Z", "iopub.status.busy": "2026-07-21T02:44:16.124738Z", "iopub.status.idle": "2026-07-21T02:44:16.185971Z", "shell.execute_reply": "2026-07-21T02:44:16.184789Z" } }, "outputs": [ { "data": { "text/markdown": [ "| preset | FLAMINGO model | native NSIDE | template |\n", "|---|---|---:|---|\n", "| `ksz5` | Fiducial L1_m9 | 4096 | `flamingo/ksz/L1_m9_lensed_kSZ_uKCMB.fits` |\n", "| `ksz6` | Strong feedback (fgas-8sigma) | 4096 | `flamingo/ksz/fgas-8sigma_lensed_kSZ_uKCMB.fits` |" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "models = {\n", " \"ksz5\": (\"Fiducial L1_m9\", \"flamingo/ksz/L1_m9_lensed_kSZ_uKCMB.fits\"),\n", " \"ksz6\": (\"Strong feedback (fgas-8sigma)\",\n", " \"flamingo/ksz/fgas-8sigma_lensed_kSZ_uKCMB.fits\"),\n", "}\n", "rows = [\"| preset | FLAMINGO model | native NSIDE | template |\",\n", " \"|---|---|---:|---|\"]\n", "for name, (label, template) in models.items():\n", " rows.append(f\"| `{name}` | {label} | 4096 | `{template}` |\")\n", "display(Markdown(\"\\n\".join(rows)))" ] }, { "cell_type": "markdown", "id": "013800ec", "metadata": {}, "source": [ "## Generate and compare the kSZ maps with PySM\n", "\n", "The workflow is the same as for every other PySM preset: create a `Sky`, call `get_emission`, and work with the returned map. The example below compares the new FLAMINGO models with the existing WebSky and AGORA kSZ presets at 143 GHz. The embedded result was generated at NSIDE 2048 with this code in the [release validation notebook](https://github.com/galsci/pysm/blob/3.4.6/docs/preprocess-templates/verify_templates/compare_flamingo_ksz.ipynb).\n", "\n", "The full comparison downloads 9.75 GiB into the Astropy cache. Run it on a compute node with at least 15 GiB of free persistent shared storage and 32 GiB of RAM; PySM temporarily uses additional memory while downgrading the native NSIDE 8192 AGORA maps in double precision. To try only `ksz5` and `ksz6`, set `presets = [\"ksz5\", \"ksz6\"]`: those templates occupy 3.00 GiB, so allow 4 GiB of free cache space and 16 GiB of RAM." ] }, { "cell_type": "code", "execution_count": 3, "id": "5ccfd716", "metadata": { "execution": { "iopub.execute_input": "2026-07-21T02:44:16.189331Z", "iopub.status.busy": "2026-07-21T02:44:16.189064Z", "iopub.status.idle": "2026-07-21T02:44:16.248229Z", "shell.execute_reply": "2026-07-21T02:44:16.246818Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlIAAAG7CAYAAAAIfHukAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvctoD+AAAAAlwSFlzAAAPYQAAD2EBqD+naQAA4KRJREFUeJzsnXd4VEXbh+9t2U3vCaGlQOi9ShMQQaUIiIUivOBnQbGi2BugoljfV1FAAiKigIA0QTooiHSkdwIppPey/Xx/LLvJZnfTSKHMfV17kTMzZ2bOIdnzO8888zwySZIkBAKBQCAQCAQVRl7bExAIBAKBQCC4WRFCSiAQCAQCgaCSCCElEAgEAoFAUEmEkBIIBAKBQCCoJEJICQQCgUAgEFQSIaQEAoFAIBAIKokQUgKBQCAQCASVRFnbE7iVMZvNJCYm4u3tjUwmq+3pCAQCgUAgKAeSJJGbm0vdunWRy0u3OQkhVY0kJibSoEGD2p6GQCAQCASCShAXF0f9+vVLbSOEVDXi7e0NWP4jfHx8ank2AoFAIBAIykNOTg4NGjSwPcdLQwipasS6nOfj4yOElEAgEAgENxnlccsRzuYCgUAgEAgElUQIKYFAIBAIBIJKIoSUQCAQCAQCQSURQkogEAgEAoGgkgghJRAIBAKBQFBJxK49gaAWMBgMmEym2p6GQCC4hkqlQqFQ1PY0BDchQkgJBDVITk4OaWlp6HS62p6KQCAohkwmw9fXlzp16ohMFIIKIYSUQFBD5OTkkJCQgJeXF0FBQahUKvGFLRDcAEiSRH5+Pqmpqbi7u+Pn51fbUxLcRAghJRDUEGlpaXh5eVG/fn0hoASCGwx3d3d0Oh0pKSn4+vqKv1FBuRHO5gJBDWAwGNDpdOILWiC4gfHx8cFkMgn/RUGFEEJKIKgBrF/MKpWqlmciEAhcoVRaFmmMRmMtz0RwMyGEVBm88847+Pn54efnx7Fjx2p7OoKbHGGNEghuXMTfp6AyCCFVBm+99RaxsbFER0cLc69AIBAIBAI7bgshde7cOf7991+X9UajkZMnT3L58mWHOo1Gg5+fn4gvIhDcwCQkJLBu3TrS0tJuin4FAsGtwy0tpH7++We6detGx44d6devn9M2W7dupUGDBgwYMIBWrVrRs2dPUlJSanimAsGtQ3JyMuvWrePcuXMOdQaDgXXr1rF+/Xqn5544cYJ169aRmZlZoTG3b9/OkCFDOHLkSLnPMZvNXLx4kR07dnDw4EGysrKqpF+BQHB7cUsLqQMHDvDZZ58xbdo0p/UZGRk8+OCD/N///R/x8fEkJSWh0+l4/PHHa3imAsGtQ25uLkOGDOHdd991qNu7dy9Dhgxh0KBBnDlzxqH+hRdeYOjQoUiSVK1zXLZsGU2bNqVRo0b07duXTp06ERgYyIABA9iyZUu1ji0QCG4tbuk4Ul988QUA+/fvd1r/66+/otPpeOONNwDw9PTk1Vdf5ZFHHiEpKYk6depUaDydTmcXsTonJ6eSMxcIbl4aN25M/fr12blzp0Pdjh078PHxQavVsmPHDpo2bWqrMxgM7Nmzh7Zt2xIQEFBt81uzZg2PPPIIAKGhoTRt2hS9Xs/Zs2fZvHkzISEh3H333dU2vkAguLW4pS1SZXHo0CGaN2+Op6enraxr165IkmQz5S9YsAA/Pz8OHDhAr1696Nu3r8v+ZsyYga+vr+3ToEGD6r4EgeCGpG/fvly9epXTp0/ble/YsYM777yTzp07s337dru6vXv3UlBQ4PRv7NKlS2zZsoW///6bwsLCMsc/ffo0W7Zscer3OGfOHACmT59OQkICO3fuZM+ePaSlpbFlyxb69+9fZv96vZ7169fz+++/s3btWg4dOuS0XUpKCuvWrePKlStl9ikQCG5ObmmLVFmkp6cTGBhoVxYUFGSrAxgzZgzDhw+31VvjjDjjjTfeYPLkybbjnJwcIaYEpSOZQZte27MoHU0gyCr2ztWnTx8WLVrEjh07aNasGVBkcZo6dSqZmZnMnz/f7pwdO3bYzrVy+vRpHnvsMfbs2WMr8/Hx4dNPP+XJJ590GDcvL49BgwbZ+WCNGTOGmJgY1Go1ANnZ2chkMqZMmWK3iUQmk7n0pSxOTk4Ow4cPZ8+ePSxdupQPPviAK1euEBcX5/D9MHPmTL744guny5gCgeDW4LYWUiqVyiF5rPVt1xo40c3NDTc3t3L1p1arbV/WAkG50KbDopDankXpjE0B9+AKnWK1Ku3YsYOJEycCRRan3r17k5mZyUcffcTp06dtQmvHjh0oFAruvPNOABITE7nzzjtJTU2lbt26NG3alLy8PA4dOsRTTz1FaGgoQ4cOtRv3zTff5Pz583Tq1Am5XM6RI0dYvHgxPj4+fPvttwB06dKF3bt3s2zZMsaOHVuh60pMTGTgwIHEx8ezdetWunXrRnp6OhMmTGDdunUMGzbM1lav1/Pjjz/Sr18/oqOjKzSOQCC4ebitl/YaNmxIYmKiXZn1uGHDhrUxJYHgliAyMpLw8HCblQksQsnb25sOHTrQo0cPlEqlbXlPr9ezZ88e2rdvj6+vL2BZKk9NTeW///0vcXFxbNu2jX379nHs2DFCQkL44IMPHMZNTk7mwIED7N+/n71793Lo0CFCQkL4/vvvSU5OBuDdd9+lTZs2jBs3jsjISEaNGsXHH3/Mrl27MJvNLq/p9OnTdO/enaysLHbt2kW3bt0AeOSRR/D39+f777+3a79q1SpSU1N5+umnr+teCgSCG5vbWkj169eP8+fPc/bsWVvZ2rVr8fPzo0OHDrU4M4Hg5qdv374kJydz6tQpwCKkevTogUKhwNPTk44dO9qE1r59+xz8o37//XfCwsKIiopi/fr1rFu3jnXr1nHhwgU6derEoUOH0Gq1dmO+9NJLtGnTxnbcsmVLpkyZgtFoZNeuXQD4+flx8OBBli1bxp133smJEyd499136dWrF82aNWPfvn0O1/L333/To0cP/Pz82LNnj82KBpZkt+PHj2fjxo3Ex8fbyufNm0fdunW5//77r/9mCgSCG5ZbWkidP3+eAwcOEBcXh9Fo5MCBAxw4cMD25du/f3/69OnDww8/zLp165g7dy4ffPAB77//frmX8wQCgXOsomj79u02i1Nx/6c+ffrYhJT13+JCKj4+nqtXrzJkyBCHz/r16zGbzQ7xpoqLKCutW7cGICkpyVamVCp56KGHWLhwIUePHiUrK4tffvmFlJQUhg0b5hBTavr06WRkZDBnzhzCwsIcxpg4cSJms9nm9xUbG8uWLVt4/PHHS/WrFAgENz+39F/4rFmz+OuvvwDLlmyrr8bSpUtp1KgRMpmMtWvX8vHHHzNz5kw8PDyYM2dOhf0mBIJKowm0+CDdyGgCy27jhOJ+Um3atLH5R1np3bs3n3zyCSdPnmTHjh0olUp69uxpq3d3d8fHx4c77rjD5RglX3ism0SclXl4eLjsx8PDg5EjR3LixAk++OADdu7caed/9fXXX/Phhx9y//33s3XrVlq1amV3fpMmTejbty/z58/n7bffZt68ecjlcp544gmXYwoEgluDW1pIffnll2W28fLycuprcT3MmjWLWbNmidx8grKRySvsyH2z0KBBAxo1amQTUtblPCs9e/ZEoVCwceNG9uzZQ6dOnfD29rbVd+7cmQMHDjBv3jynMd0SExMddt0uWLCA8ePH2yWfXbBgAVBkrbp8+TLh4eFO52wN15CRkWFX3rhxY3bu3Mldd91F37592bJlC23btrVr8/TTT/PQQw+xYcMGfvjhBwYPHkz9+vXLvE8CgeDm5pYWUrXFpEmTmDRpEjk5OTbHWYHgdqRv377MmzeP77//nu7du9t2wwJ4e3vTvn17vvrqK6fxo15++WUGDhxIx44dGT9+PC1atECpVHLx4kXWr1+PRqNh8+bNdufs27ePu+66i7FjxyKTyfj5559toscq4saMGUNaWhr3338/jRo1Ijg4mNTUVFatWsUff/yBXC6nR48eDtcSFRVlE1N33XUXmzdvtvOlHDZsGGFhYUycOJGEhARiYmKq8lYKBIIbFUlQbWRnZ0uAlJ2dXdtTEdQyhYWF0smTJ6XCwsLankqNsnjxYgmQAOmDDz5wqH/55Zdt9Zs2bXKo//zzzyWFQmFrU/wzbtw4W7tFixbZxvDw8LBr5+fnJx05csTW9tFHH3XaHyApFArps88+c+h38+bNtrK4uDipcePGkp+fn7Rv3z67+b7zzjsSIEVFRUlms/m67p2g5rld/04FjlTk+S0sUgKBoNq46667GDRoEAADBw50qB82bBinT59GJpM5tQJNnjyZIUOGsHjxYk6dOoVSqSQqKoqBAwfawg8A1K9fn0GDBjF48GCGDRvGnDlzSEhIoEmTJjzzzDN2gXEXLVrEe++9x5o1azh79izJycl4eXnRokULHnzwQbuYT9Z+g4OD7cp27tzJCy+8wFdffcU333yDv78/APfddx/Tp0/nqaeeslteFAgEty4ySarm7KC3MdalvezsbHx8fGp7OoJaRKvVcunSJSIjI9FoNLU9HUE18dxzzzF79mwSEhIICbnBA60KHBB/pwIrFXl+C4uUQCAQXAfWAJ1Hjx7lu+++Y9CgQUJECQS3EUJICQQCwXVw+vRphgwZAoBCoeDtt9+u5RkJBIKaRAgpgUAguA78/f0ZNGgQAQEBPPbYY3Tu3Lm2pyQQCGoQIaSqARFHSiC4fWjatCnr1q2r7WkIBIJa4pZOEVNbTJo0iZMnT7J///7anopAIBAIBIJqRAgpgUAgEAgEgkoihJRAIBAIBAJBJRFCSiAQCAQCgaCSCCElEAgEAoFAUEmEkBIIBAKBQCCoJEJICQQCgUAgEFQSIaQEAoFAIBAIKokQUtXArFmzaNGihYhwLBAInJKYmEidOnVYvnx5tY+1detW/Pz8WLNmTbWPVZP06tWLXr162Y7Xr19PcHAwly5dqsVZCW5HhJCqBkRAToGgCL1eT0REBH5+fixcuLDUtjk5OXz++efcfffdREZGUqdOHdq1a8fgwYP55ptvSE1NdXpeSkoKb731Fp07dyYsLIzw8HDuueceYmJiMBqNDu0/++wz/Pz8bJ+AgAAaNWrEww8/zJ49e0qdY35+PnXr1sXPz6/SQuiNN97A19eX4cOHV+r8imAwGMjOzkav11f7WDVJbm4uubm5tuP77ruP+vXr89prr9XirAS3JZKg2sjOzpYAKTs7u7anIqhlCgsLpZMnT0qFhYW1PZUaZ+nSpRIgubm5Sd27d3fZbvfu3VKdOnUkwOUnMjLS4bxt27ZJfn5+Ls/p2rWrlJqaanfO9OnTXbZXKBTSpk2bXM4zJibGdj0DBgyo8P04d+6cJJPJpG+//bbC51aGDRs2SID066+/1sh4NUXbtm2ltm3b2pUtXLhQAqTjx49Xqs/b+e9UYE9Fnt/CIiUQCKqVmJgYQkJCmDJlCn///TenT592aHPp0iUGDRpEUlISvXv35tdff+XChQskJydz9OhRfv/9d5577jmCg4Ptzrt48SLDhg0jKyuLgQMHsnXrVq5evUpsbCwLFy4kMjKSvXv38vDDDyNJksO4S5YsITMzk7S0NA4dOsSDDz6IyWRi2rRppV5PeHg4zz77LFu2bOHKlSsVuh+zZ89GpVLxyCOPVOg8Qdk88MADeHp6Mnv27NqeiuB2ovp13e2LsEgJrNTUm27qiRPSudWrpdQTJ6p1nPJy5coVSS6XS5MnT5ZiY2MlmUwmvfLKKw7txo4dKwHS448/LpnNZpf9layznvfYY485bZ+amipFRkZKgLR69WpbudUitXbtWrv2Op1O0mg0kre3t9P+Tp8+LQHSu+++K504cUICpPfff9/lfJ1Rv359qU+fPg7lX3/9teTr6yudPXvWoa5bt25256xYsULy9fWVtmzZIq1du1bq1q2bFBoaKnXp0kVas2aN3bmuLFJGo1H67rvvpB49ekihoaFSRESENH78eOnixYt27Soylslkkr7++mupW7duUr169aTmzZtLo0ePlg4fPlypsSXJ8n/41FNPSQ0bNpTCw8Olxx9/XMrIyHBqkZIkSRo8eLAUFBRU6u+RK4RFSmClIs9vIaSqESGkBFbK8wWdn5JS6U/2lSvSsrvvlj4F22fZ3XdL+amplk8l+qwKpk6darfU0q9fPyk0NFQyGAy2NgaDQfLy8pI8PT0r9LdiPc/Dw0PKzMx02W7x4sUSID366KO2MldCymQySV5eXlJwcLDTvqZMmSLJZDLbA79r165SREREuR/a58+flwDp1Vdfdaj79NNPJUA6deqUQ13Lli2ljh072o5/+eUXCZCeeuopSSaT2S1NyuVyadeuXba2zoSUyWSShg0b5nRpMygoSLpw4UKlxnrjjTec9qlUKis1dl5entS8eXOHdp07d5ZatmzpVEhZ/2+PHTtWjv8Re4SQElgRS3sCwU3ItyEhlf4s79+fy1u2IFcqCWnfHrlSyeUtW/h99GgWNG9eqT6vF0mSWLBgAV26dKFly5YAPPbYYyQnJ7Nu3Tpbu/Pnz5OXl0fXrl3x8fGx6+O7776zcwr38/PjwIEDDuf5+fm5nMe9994LwL///lvqfLOyspg8eTJ5eXn06dPHod5oNPLjjz/St29fIiMjbdcTGxvL1q1by7wfAMePHwegSZMm5WpfFt9//z3vvvsup06dIjY2ltdffx2z2cysWbNKPW/evHmsWrWKHj16sHHjRhISEjh79izvv/8+6enpvPrqq5Uaa+nSpQQHB7NmzRoSEhI4f/48S5cupUOHDpUa+4svvuDUqVN07tyZnTt3Eh8fz9q1a0lLS+PEiRNOr61Zs2YAHD16tEL3UiCoLMranoBAILh+Ms6cQa5UMvbQIYJbtyb16FEWdezI5c2bUZciMqqTrVu3Ehsba7eL6oEHHsDPz4+YmBiGDRsGQHZ2NoCD/xOATqez1Vux7sKzloeUIfoCAgJQKBRkZWU51I0cORKlUonJZCIvLw+A5s2b8/nnnzu0XbduHcnJyXz22Wd257/44ovExMRw9913lzoPsOwutM6pKpg4cSLvv/++7XjGjBksW7asTNG4cOFCAgMD2bx5M+7u7rby9957j/Pnz7NixQqMRiNKZdEjojxjeXt7Ex0dzZAhQ2xl1t2QlRl75cqVeHl5sXbtWkJDQwGoV68e9erVsxNnxbHeW1c7PAWCqkZYpASCW4SgVq0Ibt0agOA2bQhq1QoAyWSqlfnExMSg0WgYNWqUrUyj0TBy5Eg2bNjA1atXAfDy8gJwEEwATz/9NJmZmWRmZvLiiy/a1Xl7ewOQmZlZ6jxyc3MxmUy29sXJz88nOzvbJqLc3d358ssvadCggdPr8fHxYcSIEbYy6/Fvv/1W5jwAm8O7XF41X739+vVzKGvUqBE5OTmlnnfixAmysrIICwtzsPgtX76cwsJCByFSnrE+/fRT9u7dS69evZg+fTrr168nPz+/0mOfO3eOjh072kSUlfbt2xMWFub02qz31mw2l3oPBIKqQggpgeAWIe34cVKPHQMg9ehR0q4tI8kUihqfS2ZmJqtWrUKv1xMeHm73sFy4cCEmk8kWU6pRo0ao1Wr27t2LTqez60etVtvOU6vVdnVRUVGo1Wr27dtHYWGhy7ns2LEDgBYtWjjUWXftXb16lTVr1uDv78+wYcMcloWuXr3Khg0byM/PdxAAK1asQKfTsXjx4jLvi9XqlpGR4VAnk8kAnO4uLHlfrHh4eDiUyeVyp30Ux2g0YjKZyM7OdvhotVqnY5ZnrP79+3Pp0iWeeeYZUlNTmT59OiEhIbz99tuVGluSJDurWHHc3NyclqenpwPOLZwCQXUglvaqgVmzZjFr1ixMtWQJENycPHNt2acyrB05krht21jUoQNBrVqRdvw4ZqOR8P79GfTzz1DGg7WqWbx4se2h6MzSBDB//nxef/11NBoN9957L6tXr+bjjz/mvffeK9cYJc+bOnWqQxudTmfrz7qUWBxPT0+bf9WQIUMICwvjjjvu4LnnnmPnzp22dlbxV9r1xMTE8Oyzz5Y65+bNmwOWsA0lsT744+Pjbe0AkpOTuXTpEu3atSu174pgtSQdPnzYZZuS/mrlxc/Pj1GjRtkskb/99hsPPPAAHTp04IEHHqjQ2OHh4Rw+fJjc3Fw7i+KFCxe4cuWKU984670tfg8FgupEWKSqARHZXFAZPIKDK/0ZsnQp4f37YzYaSTlyxE5EeQQFVarP6yEmJgaFQsGZM2dsS3PFP6NGjeLcuXP8+eefAEydOhW1Ws3777/PhAkTOHjwoE24aLVa/v77b3bv3u0wznvvvYebmxvTpk3jueees6UHMZvN7Nq1i759+3L48GHatm1brrhNnTp14sknn+TPP/9k9erVtvL58+ejVquJjY11ej1DhgzhyJEjHDp0qNT+mzVrRmhoKPv27XOoi46OBmDatGlcvXoVSZI4deoUDz30UJW/lI0aNYrY2FieffZZEhMT8fHxwdfXl9zcXH799VfefvvtCi8/mkwm7rvvPtauXWuzCuXn53Pw4EEAW/ywiow9ePBgMjIyGD16NHFxcUiSxPHjx3nkkUdcWt3++ecffHx8qlR4CgSlUn2bBwUi/IHAyu0UR+rgwYMSIN1///0u2+zYsUMCpHHjxtnKVq5cKXl4eNhtrff09LTb9j5u3DhJq9Xa9fXzzz9Lbm5utjYajUZSKBS246ioKLst9ZLkOvyBJElScnKy5O3tLTVt2lQyGAzSzp07JUAaPXq0y+tZt26dBEjPPPNMmffnmWeekTw9PaW8vDyHug4dOtjmbb2m1q1bS40bN3Ya/mDDhg0Ofdxzzz1SvXr1bMfOwh9otVqpV69edve6+D3s3bt3hccyGAx2/1fu7u62cAlKpVI6dOhQhcdOSUmRwsLCbHUqlUoCpEaNGklNmjRxCH+g0+kkf39/6T//+U+p/weuEOEPBFZE+AOB4DYlqEULGt9/P0FO/IFqipiYGACeeOIJl2169+5N06ZNWb58uc1Zefjw4Rw7doxJkyYRFRWFJEnk5+fj4+PDkCFD2Lx5MwsXLnTwlRo1ahQHDx5k9OjRBAYGotVqMZvNREdH89Zbb3Hw4EGioqLKPf+QkBBeffVVzpw5w5w5c8p1Pffddx8NGjTg559/ti1puuLpp58mPz+flStXOtQtWbKEO++8E5lMhtlsZvDgwWzatMnhmq8XtVrN5s2bmT59OtHR0ZjNZkwmE+Hh4UycOJFvv/22wn0qlUr++OMPhg0bRlBQEIWFhXh5eXHPPfewbds22rdvX+Gxg4OD2blzJwMGDEChUGAymejbt6/Djj8r69atIzMzk4kTJ1b+5ggEFUQmSTXsPHEbkZOTg6+vL9nZ2ZX2NxDcGmi1Wi5dukRkZCQajaa2p1Ot5OXlYTQa8fX1tTlQO6OgoAC9Xo+3tzcKJw7xJpMJvV7v9IFZGlqtFpVK5bRPKzqdzvagd+bMbDKZyM3NRaVS2R70pcWqKn49Pj4+ZS6LDRs2jISEBJfL/3q9HoVCYbuG3NxcZDKZbYejwWAgPz/f6fzz8/MxmUy27xyj0UheXh6enp6oVCqn4xkMBgCn9RUZqzharbZcv+uljV0ca9gL6xysCYuL+0716dMHpVLJli1byhzXGbfT36mgdCry/BbO5gKBoEqxPuzLwsPDw+lOMCsKhaLCIgoo1wNQrVaXauVRKBRlCqeSlHU9xZk5cyatWrVizZo13H///Q71JXeklQzdoFKpXM7P09PT7lipVJZ5LaWJmIqMVZzyCpGyBJSVkiKu5D3ZunUrf/31V5l+agJBVSOElEAgENQwTZo0IT093eXWfkHF6dGjB5mZmcL6L6hxxF+xQCAQ1ALOAoQKKo9GoxHLcYJaQTibCwQCgUAgEFQSIaQEAoFAIBAIKokQUgKBQCAQCASVRAgpgUAgEAgEgkoihFQ1MGvWLFq0aEHnzp1reyoCgUAgEAiqESGkqgGRa08gEAgEgtsDIaQEAoFAIBAIKokQUgKBQCAQCASVRAgpgUAgEAgEgkoihJRAIBAIBAJBJREpYgQCQbXz9ttvk5WVxeOPP067du1KbRsfH8+WLVu4cOEChYWF1KlTh/r163P33XcTFBTk8rxTp06xadMm4uPjcXNzo1mzZgwcOJDAwECHths3bmTt2rW2Y7lcjr+/Px07duS+++4rM5Huq6++SkFBAZMmTaJ58+alX7xAILilEUJKIBBUKwcPHuTDDz8EICcnhx9//NFpu4KCAl588UUWLFiA0Wh0qFcqlTzyyCP89NNPduU5OTk88cQTLFu2zOEcLy8vpk+fzosvvmhXvn//fmbNmuV0Hi1btmTLli3UqVPHaf3OnTv59NNPATCbzXz77bdO2wkEgtsDIaSqmf+7O4ZF3xzgmTf6IpPJans6AkGNM3/+fJRKJXfffTcrVqzgm2++wcfHx66NwWDgnnvuYdeuXahUKoYOHUqHDh3w8vIiOTmZuLg4tmzZwq5du+zO0+v1DBgwgL179+Lp6cnDDz9MixYt0Ol07Ny5k82bN/PSSy9RUFDAm2++6TC3iRMn0rJlS0wmE1euXOGnn37ixIkTvPbaayxcuNDl9bi7u9O9e3d++eUXPv/8c9zd3avuhgkEgpsLSVBtZGdnS988tVX68vFN0tED8bU9HUEtUlhYKJ08eVIqLCys7anUKIWFhZKfn580dOhQadu2bRIgzZ4926HdJ598IgFSs2bNpPPnzzvty2AwSLt27bIr+/jjj23nXblyxeGc1atXSyqVSlIoFNLZs2dt5dOnT5cAae3atXbtz507J8lkMqlOnTpO55CdnS15eHhIY8aMkdasWSMB0k8//VTmfRDcHNyuf6cCR7KzsyVAys7OLrOtsEjVACqFgi1LT9G6Y73ankqFKSzQE/O/v7kan82wMe3o3C28tqd0a2E2Q3p6bc+idAIDQV65fSkrVqwgKyuLxx57jD59+hAVFUVMTAxPPfWUXbvvvvsOmUzG0qVLadSokdO+lEolPXr0sCubPXs2AIsXL6ZBgwYO59x///1MmTKFjz76iLlz59qW5FzRuHFjAgICyMvLc1r/yy+/UFBQwGOPPUbv3r0JCwsjJiaGMWPGlNqvQCC4dRFCqoaQtFJtT6FSfDfjL0g0UVflzfLPDhDxbQDBod61Pa1bh/R0CAmp7VmUTkoKBAdX6tSYmBhCQ0MZOHAgMpmM8ePH8+6773Ls2DFat24NWJzLY2Njad++PW3atLE7f8eOHSxfvtyu7JVXXiEiIoKEhARiY2Np06YNHTp0cDmH8ePH89FHH7F79+4y57t7927S09NdpneKiYkhIiKCvn0tS/Vjx47l008/5eLFi0RFRZXZv0AguPUQ4Q9qCLXHzalZM87kolZZ5t4w2J/l8w7V8owENwuXLl1ix44djB07FqXS8jv0n//8B7lczvz5823tkpKSAJxaoo4cOcKsWbPsPtb2V69eBaBJkyalzqNRo0bI5XJb++LMnj2bZ599lqeffprBgwfTt29fFAoF77zzjkPb48ePs3//fsaPH2/zd5wwYQKSJLFgwYLy3BKBQHALcnM+3W9CjB43p0UqLMDeKTgvQVtLMxHcbMyfPx9JkpgwYYKtrGHDhvTr149FixbxySef4ObmhkKhAMBkMjn00bdvX77++msA1qxZw+bNm2118mvLjc7OK44kSUiSZBunOL///rvdsbu7OytXruTee+91aBsTE2Ozqllp1qwZ3bp144cffmDq1Km2OQkEgtsHIaRqgJSsPNwC3Wp7GhXGoHfcgi6Xi52HgrIxm8388MMPeHp6OoQHyM/PJz09ndWrV/PQQw9Rv359AE6ePOnQT9u2bWnbti1gWQIsLqSsPlHHjh0rdS7Hjx9HkiSnPlTWXXtarZb9+/ezbNkyXnvtNXr16oWnp6etnV6v56effsLX19fBz8poNBIfH8/GjRu57777Sp2LQCC49RBCqhqwLkGYTCae7/sdOQVaPPJre1YVJy3ZcdImj1qYyK1MYKDFB+lGxklAy7LYuHEj8fHxAC7jNcXExPDQQw8RHBxMq1atOH78ONu2beOuu+4q1xjFz9u0aRMDBgxw2s46ft++fR3qBg0axODBg23HnTp14tVXX+Wdd97hiy++sJWvXr2atLS0Mq9HCCmB4PZDCKlqYNKkSUyaNImcnBwWvXqA7AIt5N98Jv/UxFyHMpn7zXcdNzRyeaUduW9kYmJiAHjnnXcIceJMP2fOHDZv3kxcXBwNGjRg8uTJPPbYY4wePZrFixfTr18/h3Oys7Mdyl566SX+7//+j3HjxrFq1SruuOMOW53ZbOarr77i+++/x9PTkyeffLLMeb/yyiusWLGCWbNm8eyzz9ocyK3X88EHH+Dr6+tw3v/+9z/WrFlDWlpaqdHXBQLBrYcQUjVATr4WZV7lb3V+no6F3/5DdkYhj/xfR6Kia+bBm55kb5EymsxodYYaGVtw85KWlsbatWtp3rw506ZNc9rG3d2dxx9/nB9++IF33nmHCRMmsHXrVhYvXszdd99Ns2bN6NChA97e3qSlpXH48GEuXrxIvXr17JboJkyYwIYNG1i+fDndu3enW7dutoCcu3bt4tKlS8jlcubOnesyUnlxZDIZn3/+OT179uT1119n2bJltiXFzp0789Zbbzk9z2Qy8eKLL7Jo0SJeeumlyt04gUBwUyLMCzVAdoGW/Hy9XVlWZiEfvrSB18f9xqY1p0o9/39vb8dwzoBPlhuz3/gTnc7Rd6k6yEwtsDvW6g0UlLgOgaAkixYtQq/X88QTT7hsM3LkSLy9vVmwYAGSJNnO++yzzwgKCuL06dP8/PPPzJkzhxUrVnD58mXGjRvHwYMHqVevKB6bTCZjyZIlvPPOO7i7u/P3338zb948Fi1axKVLl4iKimLNmjWMHj263PPv0aMHw4cP59dff2XPnj0sWLAAs9lc6vWMHTsWjUZjs1wJBILbB2GRqgFyCnR45ersyr6f/hd+hWr8PNQcXhFLq/Zh1G3g53CuJEkos0Gptmje8CB/tq89w70Ptqz+eWfY79DTGowOglAgKEnjxo35+uuvSw1S6enpyc8//0xsbCzZ2dn4+fkhk8l4+eWXef755zlw4ADnz59Hq9VSp04devbsib+/v9O+FAoF06ZN49VXX2XXrl12SYs7derkdCfdvffei5+fH61atXLa55dffsldd91FYWEhLVu25Ouvv2bUqFEurycgIICff/6ZhIQE8vPz7RzVBQLBrY0QUjXAmL7t0Rvtt2h7FBTdei93Ncu+PciLMxz9QpKv5uKhtt/xl3wh57rmI0kSa1ccY/+eK3TrFcl9Q1s4zQNYkGMv/rR6AwUmIaQEpTNkyJBytSvu5F0clUpFt27d6NatW4XG9fLychq2wBmdOnWiU6dOLuvDw8N59tlnKzT+8OHDK9ReIBDcGgghVQMo5HKU8tLjSOUkFvL+k+uQJInRL3ahactQAC6eSHVo6+FTeigFSZIoyNej1qhQKh3fxrdvOsu01zYAsGHVCfz83ene2zEqszbfgJqi2DtavZFCsxBSAoFAIBBYET5SNYRSIcdgcB04MNDdg2C5ByEKT5Z+vN/mNxJ/IcuhbWn9aAsNTHlsJS898CvPPbSUxHjHnU5zvtpld/zJe5uc9mXQ2o/TuG4QbYLLdtgVCAQCgeB2QQipGkImk5GTWb6o4IGeHpw4lAhAupMQBHqta2fzDctPEK7yo21UXVqFhrL02wMObTIS82kVXoduzcLp3z6auu4+TnoCSe9oRZM7WQIUCAQCgeB2RSzt1SDZGQUEhniiK0UIWcnPsiyh5Wbo8ERlV2fQubZI/bP6Ag2DipxyTVcdxwr186Zzk6It5GqV80z3mB2LlDKhvQUCgUAgsCKeijVIbqbFeTuvhBO3M6yGH2OBo2gqTUjVD/SzO3ZXqxza5Gntx/fSqDGbHa1PcWlZZObZh0BQKhSlLi0KBAKBQHA7IYRUDZKbbVnay8suf+Jfudl+Kc1sNqN3Ziq6RnkW3gr09g7j7moVKYn2OwFNJjN7T15h65Hz9vORy2zXIRAIBALB7Y4QUjWINZxAeSxSOq2RnGwtHm72FqX1B05TgOvo4ievJDuUWR3XrRjMjkLs8tl0u+MvP9xmaWtytD6V19dLIBAIBIJbHeEjVYPk51osQRISRy4m4qFWIZfJCPX3xttdbddWrzWREJfFit3H8HZX4+2hxsddQ3a+Fm2hayF1LjGNluFFO+usoRA8vSz9GwwmdFojWr0RjVvRf3/CpaLdfUcOxLNk4UFLe6Oj6MrNLKzE1QsEAoFAcOshhFQNUphvEUAGs5nDFxJs5XUDfLinY1O7tlcSM5n17V/ka/Xka/UkZRbt3issRUhpDfbO5TKZjNSreXhGW4RUdnoBY/p2cDgv/WqRw/npE0VWLZPZjFmS7Hbr5WWXbVETCAQCgeB2QCzt1SDaa0KqsMBeCCkV9v8NOQVaZn35J1cuZTjvp8C5kDLoTej0jrv00pOKRFJ6cr5DPUBuetFyXcldhYYSUdnLszQpEAgEAsHtgLBI1SD6QotAKWlRUpTIBWZy4sNUHK3WuZDKyizALEnojSbclEURyTNTinbeZaYXODsVXV6ReNLlG/Dz1GA0S07nUpDn2iImEAgEAsHthBBSNYg1bEFJi5RCbr/Xzt/Lg14tIzGZzZjMEhm5BZxLTLPVlzzfStY13yWd3mAnpLLTi3yaSiYitiIZihzSDakGhndv7fI6tHkiTYxAIBAIBCCEVI1iFVIlncXPX00nNjkThUKGQi6nbVRdmtUPsdXna/WE+nujlMtRKuQo1M5XZK3WJq3BhHex8tysIvGU62LHnVquwGgwoVQpMBrNKEsJpKAtKDugqOD2RqvVcvz4cZf1kZGRBAYG2toFBwcTHh5e7v7T0tKIjY3F39+fRo0aOW1z5swZcnNzUSqVtGvXzmmblJQUrly5AkDTpk3x9vZ2OL9jx462pN7WMj8/Pxo3buzQX3x8PElJSbRs2RJ3d3enY2ZnZ5OYmIhcLqdu3bp2Y5ZGYmIiaWlp+Pr60rBhQ6eJxkvD2fWUxfnz58nKyrL9f90I6HQ6Ll68iFwuJyoqCpXKMVaeQFCTCB+pamDWrFm0aNGCzp0725WbDJZlssICPUq5HG93NX6e7gT5eOLjqUFnMFGgMzj4OalVSqLrBhFZJ4AGwX74qjVOx40/nckjd7Yj2NfTrrwgp8iCVJDr3JqkcVNx4vBVAMzG0hMs60pxdhcIwPIA7ty5s8vP2rVr7dpNnz69Qv1PnjyZzp07c+edd2JyEqID4KmnnqJz5860b9+eEydOOG3zzDPP2OZ08OBBp+frdDqHsk6dOpGR4ejD+NVXX9G5c2fOnDljVy5JEgsWLKB9+/b4+/vTokULmjVrho+PDyEhIUycOJEjR4449Gcymfjyyy+JiIigXr16tG3bloiICOrUqcObb75Jfr5zn8fS7kfx63HGmTNneOONN2jatCnR0dF2/1+1SWFhIc8//7zd/QsNDeXzzz+v7akJbnOERaoamDRpEpMmTSInJ4dXHlyDwWRCbzTRsEkAAPlxWsb262h3TmZeAav2WL7sSzqfGyUzymKaV+HibTI3Q4uHk0jmuvwiYVaY61oEHd+bSNsu9TEbzZSmsUsmMxYIXBEUFOTU0hQUFFTpPnNyclixYgV169YlMTGRP/74g0GDBjltq1Kp8Pb2Zu7cufz3v/+1q0tOTmbNmjXUqVOHpKSkCs0hOzubjz76iM8++6zMtiaTiYcffpiVK1cCoFariYiIQKFQkJCQQGpqKnPmzMFoNDJv3jzbeQaDgaFDh7JhwwYAAgMDqVevHqmpqVy9epUZM2bwxx9/sG3bNvz8/Co0/9JYsGABn3zyCQB+fn5kZWVVWd/Xw6OPPsrKlSttliiZTMalS5d45ZVXkMlkTJ48ubanKLhNERapaubQhQSOxSZxJj6V9HzL0puzFC/+Xh48elcH/tOvk10cKAAz9haiks7pVgqc7Ka7nJJJnr6oXOfCvwog4XymZTxT6RYpIaQE5WXo0KEcOHDA4TN48OBK9/nLL79QUFDAvHnz8PX1JSYmxmVbuVzOuHHjWLRoEYWF9vHPFixYgEwmY8yYMRUaX61WM2DAAGbNmkVcXFyZ7d9//31WrlxJaGgov/zyCzk5OZw+fZoTJ06QlZVFWloas2fPpn379g7nbdiwgbCwMNavX09qair//vsviYmJ/PPPP7Ro0YLDhw8zceLECs2/LJo1a8aMGTM4ffo0X3755XX3d+zYMZuFzmQycfbsWeLj4x3aJScnc+LECfR6R6v5iRMnWLlyJQ0bNuTEiRNcuHCB8+fPc/z4cRo0aMB7771HTk6Ow3kCQU0gLFI1SP41J22jwYybEw2rUigcygAHuauQyzEaTSiV9u21+Ua7fi9cTefP4xfp5N7QVqbXmtDgfBxtpkVkSSXy7mXmFXA5JROD0YzBZKJ5xzrOThdUArNZIvsGD3Dq6++OXF4xf5zqJCYmhqioKO69915GjhzJ/PnzSUlJISQkxGn7J598kq+++opff/2VcePGAZaltu+//57hw4dXyjr28ccf07FjR9577z3mz5/vsl1OTg5ffvklbm5ubNu2jRYtWji0CQwM5KmnnnI477///S8qlYpNmzbRqlUru/quXbuybds2WrVqxbJly5g6dSpNm9rHoqss48ePt/28Z8+e6+5vyJAhRERE8P777zN27FibiBo8eDDLli1Dq9Uybtw41q1bB1islUuWLKFfv362Po4dOwbAxIkTadasma28efPmTJw4kbfeeov169czcuTI656vQFBRhJCqQQquCSmzofSls5LIVXIoZgSSy2UU5Onx8StyZt2/5zJXL2UTHupvK9PqLcIoq9iD2qgzQQkhFZ+WRVpOPslZeeTn6ZBKNLmSlsXhC4m24/qtA8o9d0HpZGcW0r/L17U9jVLZvO85/AM9KnVuWloaBw4csCvz8PBwKijKw/Hjx9m/fz/Tpk1DJpPx2GOPMWfOHH788UdeeeUVp+c0b96cXr16MXfuXJuQ2rJlCxcvXmTevHns3bu3wvNo3749I0eOtI3r6np27txJfn4+o0ePdmhjdeS2olAobFYp63kPP/ywg4iyEhoaypNPPslHH33Ehg0bqkxIVQdJSUkMGzYMs9lM8+bNuXTpEuvWreOdd97h8OHD7Ny5k+joaLKzs0lJSWH8+PFcvHjR5khuddy3bgwojrXs8OHDQkgJagWxtFeDFORbhJSpDGfukqjcHS1IxYNi/rX5HOu/PmonoqAoynlxIVXSkfzAuXg2Hz7H4QuJJKbncOpYkkNuPjeNvd4uzBfhDwTlY/Xq1Q6O5qNHj650fzExMcjlcpvVpEuXLrRs2bJUqxBYrFK7d+/m5MmTAMydO5fo6Gj69OlT6blMnz4duVzOW2+95bLN+fPnbfMsyYsvvmh3X3r16mWrO3fuHADdunUrdQ533HGH3Tg3KmfOnGHs2LGkpaVx8uRJzp07R2hoKF9++SWJiYmcPn2as2fPkpiYyLhx44iPj2fLli2283v06IFGo2Hu3LlMmTKFrVu3snXrVl555RUWLlwIUGE/N4GgqhBCqgbJy7OIH6kMH6SSqD0cDYf5xXbi7Vl7kSAfT4c21t1/WZkFReKoxNh6o/0OwasJOZRwyULtbj9+gQs/q/jLmVw4ner8IgS3JUFBQXTs2NHu07Jly0r1pdfr+emnn+jXrx8NGjSwlU+YMIFTp06Vugz14IMPEhAQwNy5c0lJSWH16tU88cQTFQ4hUJxGjRrx5JNPsmrVKpdjW/2yPD0d/z6jo6Nt98TDw8PpeV5eXqXOwVpfUOA80O6NQnBwsG2JE6B+/foMHz4cs9nMJ598YgsloVAoeOaZZwC4dOmS7fygoCA+++wzJEnis88+4+677+buu+/miy++sDmZO/OtEghqArG0V4PotEaMRrPFB6kCElbjqSIPe8FTPIyBNtOAp6fjbj2rRcpoMJOfp8fLW428RHwovcHecVynNyIrIaTcPez7Lihw/MKaM/0v9HE6zJKEW10VT0/tU9ZlCW4Dhg4darcT7XpYvXo1aWlp3HHHHXbLhc2bNwcs1ipXFhyNRsO4ceP48ccf8ff3RyaT2fkCVZZ33nmHhQsX8vrrr7Nz506H+oAAyzJ4QkKCQ11xR+5OnTpx+vRp27E1ZpMzp+ziWOuvZxdkTdC0aVOUSvvHjfUaSwpra3lmZqZd+aRJk2jXrh0xMTGcPXuWgIAA2+5oKLrXAkFNI4RUNXN/1xaolArLR6EgO6MQWekZYBxw0ygwmyU7h9+CvKKlPU+1m9PzrD5SYFne8/JWo5SVSEeD/WT0OiOUEFvuJURayVx/h/6+gi5Oh0IuRw5oEwxkpRfgV0m/mtsJX393Nu97rranUSq+/s4DS9Y01t1506dPdxp3aunSpfz3v/91av2BIqfzDz74gBEjRhAcHHzdcwoNDeWll15i+vTprF+/3qG+U6dOAKxdu5b33nuv3P1aY9CtXr2a999/32W71atX241zo6JwtZGmlLqSLgZgWeLr0aOHXZlVELdp06byExQIrgMhpKqZwBJLbjmZhZQUKmWhUisxmc3I5UVfOD/HHMCkhJ59G7Hl37N0impAHX/7CMmtI8KINhjRqJRcjc0iJNQLVYmdfu5ebpBcdKzXmxwtUp72Qq340p4kSayae4RQTdEShJtSwaXTabTv0RBB6cjlsko7ct9OxMXFsXnzZgICAoiMjHSoz8jI4NKlSyxbtowJEyY47aN58+Y8+uijnDp1imeffbbK5jZlyhRmz57NG2+8YbfTDLAtZR48eJBPPvmE1157rVx9dujQgRYtWnDkyBFmzpzJq6++6tBm+fLlrFy5koCAAJdxtG4lCgsLHaLFHzp0iJ9//hmlUsn9999fSzMT3O4IIVXDbPjtBJgr5iPlplFgkiSK24UuX8zgxceX88GXQ0jNyGdrzjlG9GiNxq2oVYNgP9vPKfG5XPBPZeexi7ipFKiVStyUCgJCPeBCUb96nRF5Cb8RD5WKYd1aolQoUCnkmIq9KV48m4a/wjHSena281Q0AoEznO3us9KpUyd++OEHzGYzH374odO4SZcvXyYqKoqYmBiXQgpg0aJFVTZnK97e3rz11lu8+OKLDsErZTIZs2bN4u677+b1119n48aNjBo1iiZNmqBUKklKSmLr1q2cOnXKTiTIZDK+/vpr+vfvz2uvvcaePXsYM2YM9evXJzU1lVWrVrFw4UIkSeLTTz91aYVzxqFDh2y+SsUJDw8nODiYnJwczp49C0BsbKztX+v/T5MmTfDx8angXbp+PvjgA+Li4hg6dCghISH8888/zJgxA4PBwEsvvURYWFiNz0kgACGkqh2T2WwXQFOKNeGlUVeoD7VG6RCUU3mtz+lvWKIe640mNh06S7fm4dRv6I8u2375LSulgBPHkricUuR30DDCn86+DejfPhq5TI5CIUOfYnAQUu6ebvh7FVlNjCYTkiQhk8k4fzQFN1XRr5EkSazYfYwGd9/YPhuCG4vVq1fblqlKYjabWbBgAZ6eni53/IWHh9O/f382btzI6dOn7WIN1QRPP/00X331lU14FKd3796sWrWKcePGsX37drZv3+7QplmzZvzyyy92ZXfddRdLlixhwoQJrFq1ilWrVtnVq1QqPvnkEx577LEKzbXk0piVr7/+mmeffZa///6b++67z65u6tSpTJ06FYANGzZw7733VmjMqkCtVrNo0SIHMfzAAw/w8ccf1/h8BAIrQkhVM3qjEXcnb38VQe2uRCoppK75Fei0RU7o6bkFbD95gZ2/vMSnT23EQ1ZkncrN1HIhKd2ujxZtwnAzKqkf5GcrM2rNKEoEX/TyVZNHUQgFpUJBYaEBDw838jLto6nnFerILdSRmpxXuYsV3BK4u7vTsWNHIiIiytWuNI4fP05AQACjRo0q1RLywgsvkJaWxp9//kmzZs1o2rQpBkPZeSHDwsLo2LGjQ/Lgpk2bkpeXh7zYi5CrPt3c3Jg5c6YttUrJXXiDBg0iNjaWxYsXs3PnThISEpDJZERERDBw4EAefPBBp75CDz30EL1792b+/Pns3r2b9PR0fHx86NSpE+PHj3eaONkV1utxhTWgqa+vb6n/J76+vuUeEyy+S86sRfXq1aNjx46o1fYvlhqNho4dO1K3bl278nfeeYdmzZqxYsUKEhISCA8P5+GHH2bo0KEVmo9AUNXIJGcefYIqIScnh2+f+wtvF1ngyyIzrwCFXE7fcc34c8lZPJVFguzwhQRikzORyyAjr0jkNIwMYOWWJ/jy+c246Yu+mPUeZvaevsL5M0XhCSa/fRcpR7NRFRY9KHQaE1v+PItCIUMhl6OQy3n+7T4cX2u/e+ihdzsTUtebhZ/tIe980fip2fms23eSh8Z24LX3+5d6fYUFehZ+tYeCPAP9H2yOr787DSL8r2tL+o2KVqvl0qVLREZGotE4TzotEAhqF/F3KrCSk5ODr68v2dnZZS5lC4tUNeMqM31ZZOYV2pIYP/RWZwwyM8lZuYT6Wd6a2zeqR/tG9cjILWD1P0WZ7YNDLE7fak8Vkr5oR15hrp6L59LsxmjZJoy0E7l2ZQa9iZRs+7fW4Hr2b+oAuVlaQup6U5hn/3ZujUuVVg6L1JevbMXXrEaNnD/nnuHAuXgi2wfy7me3vuOsQHAzc+HCBYfwBCUJCAggKiqqhmYkENQeQkhVM9kF2fiVEVTPGcWX19RqJbkaPVu2naFeoC8DOjSx1QV4e+Dn6U6hXo/OYCI41DKWp5+avGIRzfOzdJiLObkrFDKatghll/IcxSNwmgyOsRmc7SrLy7Us6enyDSiL7ULUXYtdlZqS63BOSXzN9ib9TtH1yU3RcuFMKo2aXv/WdIFAUD28/PLLLn3arDzyyCMsWbKkhmYkENQeQkhVM4cubiU8ZJzTupSsPFRKuZ0jt5XiDupqjRJ3d4u/k1WoFGd4d0suLpPZDErLeT4BGvIuFQmpuoG+NAz2IzU7j0K9kUZNgtG4q5Cr5JiLJfJzJqQ8vN0wS5KdE3rBtRQ1Bq0JZbFfI921AJ9pKflOr9mKqxVlbw8N21edodFrQkgJBDcqjRs3LtO3rVGjRjU0G4GgdhFCqprZ+u9CPvz4LY5sTsSYbxEZJrOZ9JwCft9/CrlMxn/utg+md+JyEmk5RULETa1E467CQ62iZ8sIl2Mp5HLUXhbBFVDHk0TsTe/92kUDFodw87UlX6VSjr64kHKSvkatUWIymZEXi0GVf21Jz6Q324XFat4gBD9PDWqVErPZbOeoWxy91vWSp1kr3PYEghuZzz77rLanIBDcMAghVQM07x5KcIQ/Ex/5mXydwc6qpFA4Co0Tl5PJ1xWlYdFoLELK213t1HpVHJ9Ai2N7aH0fjrto4+WuRrrm/65UySme8MVZHkA3N6VDHCttnuUsyShBiew0YQEWlZaelEdwXedOegajEaPJjNLJ9VvJyday6Pu9mM0SY5/oit8NEmFbIBAIBAIrQkjVEHUjfO1211lRyh13qBnN9strbmol7h4qPDVlh1HwD7EIrbCGpe8y8PS1+Ccp3ey3XMskcFerMJnMmMwScqUMuVzmEMeq0BrdvJR0Nwmx2S6FlNFgJl+rx9fT9c6YVyau5NC+OAD277nCj785XyIVCAQCgaC2EEKqhtBoVPj4asgpEfFb4WTpq1N0fXQGi8VGklnSiGg05RNSIdd22PkEuDsEAy2Ob7BFcJUUUsG+Xoy8s53t2GiyKCWpRN4YXYHFqqYoJVRB7Jl02nVv4LTOYDBhNrtWYRlp+TYRBXDy6FXir2RRv6Gfy3MEAoFAIKhphJCqQYJDvRyEVL5Oz5I/j6C8FrOpjr833ZqH2+qN18InKI0yOkU7FyXFCQv3AywpJgwmk0shFRRmSSnhpnadTBTALFnETolcx+gKLUJK6aJ/gMunM1zW6fUmjC5S5cgUUJCvdyjPyxFpZwQCgUBwYyGEVA0SEurNhbNFsZwUchlKhdzmq601GMjKt1/+s+a1k5UdoBmAwDpFObdMctdO23UaWqITl7RIlcRs3V1XYglSrzViNJhYt+8UapUStVLBna2i7NLF5KU6LmVaMehNll2GzpCwC9Vg4xYM1CkQCASCmxshpGoCQw5kpxAc6kl0vSDaR9XDYDLh52nvPH38chKJ6dl2ZWZJgsyTuF1dDdhnli+J3mRCoSyyEMmVMnChpepF+AGg1pQupKynyxUyim3uw6gzkZOtJTu/yEq05/QVercuCsCnlpQY9EZUbo6/Znq966U9s1nCZDLTvlFd5DI5crkMuUyGodAx9INAIBAIBLWJEFLVTMdw8FrXDnTpBGdPJC/obpe+Tq3C69AqvI5dmYQEfz2Fm7YeZQmpkrvnVJ5KclK0eGnUyItZlIxmMx4+ljm4aUr/FbDm+FOo5PZCSm92WKZMy7GPZu6uVnH8wFXaO/GTKsjTEezrPFCp2Sih15toF1XPrtzqlyUQCAQCwY2CawcXQZXw9kCQ6yzJgj1l7oSH+FfofEmSIGkXhUZ1mW3d/ewFmm+0Jyt2H2PHsQt25QazyZbPTl2mkLIg95Jz4nIS/15M5MC5OPIkvYOQ0ktmDCVS4pw5kuS03+zkQlRK59Ywk9FMYZ6jj5RBV7l0OwKBwB6TyUR8fDxXrlyp7akIBDc9wiJVzQxrV/SzvyaJXNduQy6wiIcQ7wJSS09tRVA9ewtP4yaW6OAe6hKmqmL6xcNXTXxaFiazhMlsJqpOoH3ba4YsTYgb+84W7aILauLjIKR8fNXoJROqYgNoc507d+lLEUUmg5nsFMcbZdALIXWzYDAYSEhIcFkfFBSE17XUSVevXkWv1xMeHu6yfUny8/NJTU3Fw8ODkJAQp23y8vJIS7P4JIaGhuLuJHm4JElcvnwZAC8vL4KCgmx1V69eRafTERERYSvLyckhIyMDmUzmdL5arZakpCT8/f3x9fV1OX+DwUB6ejoeHh5lJkQteU25ubkEBASgVpf9cuVsfk899RS//PILBoMBX19fsrKyKtxPbZGWloaHhwceHqXH0xMIahJhkaoGZs2aRYsWLejcubNdeYin6weLKyTJIkQ61dvuMq2KlYZNAuyOO3ZtgMpNgYfa3lKlKLZTzz/Ug82Hz7Ht3/PsPHaRvWdKvKFe+w3x8LDvo7BQT26JXXTePhqkEv7gJqNzPyhX1qVD5+PRqozodY7LeEYhpG4azpw5Q2RkpMvP8uXLbW2HDBlS4XQib7zxBpGRkbRv395lYvDly5fbxvvoo4+cttm6dautzSuvvGJXN3z4cCIjI+3K5s6dS2RkJBEREWzevNmhv3/++YfIyEi+/vprhzqj0cicOXPo0qUL7u7uhIWF4evri5eXF127duXLL7+0Cb/iFBYWMmPGDJo0aYK3tzd169bFw8ODbt268euvv7q8R8749ttv+fHHH5HL5TRs2JCGDRtW6Pza4NKlSzz66KN4eHgQHByMp6cnERER/Pe//y01hIpAUFMIIVUNTJo0iZMnT7J//3678ro+lyvcl8xsEStebum0DVlWatvo1vZv5u4ebrTvXN/BIqXxKjp2U9sbJRUldufJrh27lxBSBfkGsjMK8dS4oVYpUcrl+PhqHMIkmIzOxZ/BiVACiK4XjNZkdCq0hEXq5sPLy4vw8HCHj1clEnlb0el0LF68mICAABITE9mwYUOp7QMCApg/f75TwTV37lwCAwOdnFU2b7zxRpkvN1Zyc3O56667mDhxou17ITQ0FH9/f/Lz89m3bx+TJ0+2E5gAmZmZ9O7dmzfffJNz587h4eFBWFgYCoWCf/75h4cffpinn3663HPesWMH3t7exMfHc/nyZY4ePVr+C64FcnJy6NmzJ4sXL0an0xEaGoqXlxeXL1/mxRdf5M0336ztKQoEQkjVJHV8Km6RkklFVp9eDebyWJvhDLhzG+cSUu3amcxm/EIczd3d74zCvYSQ8vIvWhIoGUeqZNwpq5Dy8LTvo7BAT25cIQ/3asvoPu0Z268jjX0Cbb5Xtnm5sEgZnSRHBjCbzej1Jqe5+ISQckSSJAw6U6mf8j7sq4NHHnmE2NhYh8+DDz5Y6T5XrVpFRkYGs2bNws/Pj/nz55fafty4cSQmJrJu3Tq78pSUFFavXs24cRWPmN+vXz8OHjxYbovQ//3f//HXX38RFRXF8uXLycnJISkpiYyMDPLz89m/fz+TJ08mONg+WfcTTzzB/v37adKkCZs2bSI3N5fExESysrKYO3cunp6ezJ49mzlz5pQ6fn5+PrGxsVy5coXg4GDy8vKIjY11WNYzmUwkJSVhMFgs4TqdjtjYWPLy8pz0aqlPSkpyaRUsTkZGhst+XPHHH3+QmJhIr169iI+PJykpiZycHFatWoVKpSrz/14gqAmEkKpBVPKK7zqTYX+OuzKHBv6xJGXm2pUXGo020VOcbndGOiztefkWCSm1g0WqhJBSXLNIudsLqYICg4OwkStkDhYps6ulPReiyGSW0OtMTusNBiGkSmLUm5n7ws5SP0b9rbX8ERMTQ0hICCNGjGDkyJGsW7eOlJQUl+27du1K27ZtHcTGggULkCSJ8ePHV3gOL7zwAvXq1ePtt9/GaCz97/rff//l119/pU6dOvz999+MGDHCzsfHw8ODTp068fnnnzNixAi781asWEFISAh//vkn/fv3tyUB9/Dw4IknnuC3334DYNq0aaXOY+3atURGRvLvv/9y8eJF23LmV199ZWvz6aefEhwcTFhYGN7e3kyePJk9e/YQGRnJkiVLbO0MBgNz5syhXbt2eHp6EhYWhr+/P88++yxarf1yv9Fo5NVXX8Xf35/AwEC8vb2JjIzkv//9b5n3DbBd70MPPURYWBhgCTY8dOhQIiIinCZF1+v1zJw5kw4dOuDr64tKpcLNzQ21Wo1arUaj0bBq1aoyxxYIykulnc2Tk5Od+ghUltDQUPr3719l/d0qyGWOu9c0pHP+ajqxKZn4emrw9XAnKtr58kRk40BUKvsvm0YtihxqSy7tlUwiLL8mpDQlhJROa8SkN9spcZlCZvFNL/bcNjtJggy4fLibzGZ0OqNTIWUy1J5lpSYozHX8vy6L8ljpCvP0GPWlxwtzhrt32SmJysJq+SiOm5sbdevWrVR/ly9fZuvWrbz44ouoVComTJjA7Nmz+fHHHx18nIrz5JNP8txzz3HlyhUaNmyIJEnMmzeP4cOHu3RWLw13d3fee+89nnzySebNm8fEiRNdtl2zZg1gEV+hoaF2dXFxcXbWHI1GQ506lhAoa9euBeC5555zOM9K//796d27Nzt37uTgwYN07drVaTtPT0/Cw8O5evUqZrOZevUsoUX8/PwA+Pzzz3n11VdtbQG+/PJLDh8+7NDX5cuXbdfr6emJRqMhPT2dWbNmkZmZyeLFi21tp0+fzqeffgpAYGAgJpOJ2NhYXnzxRUaMGEH9+vVd3TYA7rnnHurXr89nn31GnTp1aNOmDfn5+fz888+cO3eOt99+2669TqfjvvvuY/v27YwcOZLnnnuO1NRUZs6cSXp6Ot27dycyMtLBf1UguB4qLaROnTrF2LFjq2wivXv3vi2EVLBqLamGIeVuL3cW0lxrCadgNJlJzykgPaeA/iObOz9fLqduhB85V6/tgpNDs05htnqVyvKAVchlBHh74O9lv7NJobQIKUOekUGdm6NUWFLZmCQzRqMZt2JSSqGUYZRLnLucilmSMJklGrW1X6qwYnRhXbJYpIxOl/5Mt7hFav6UXdd1/oSZPVFdW6o16EwseNXS37KP9qPLr7g1dNLsu65rPgBLly5l6dKldmVt27blyJEjlepvwYIFmM1mHnvsMQC6dOlCq1atiImJKVVIPfroo0yZMoV58+Yxbdo0tm3bxvnz55k9e3al5gHw2GOP8fnnnzNt2jTGjRvncifZmTNnAOjRo4dDXefOnUlOTrYd9+7dmx07dgBw+vRpAHr16lXqPHr16sXOnTs5c+aMSyE1ZMgQhgwZQqdOncjLy7P1DZZlv6lTp+Lh4cHixYsZOnQoAKtXr2bMmDEOfbm5ufHCCy/w/PPPExVlCcB76dIlJkyYwC+//MJHH31k29G4fft2AgIC2LlzJ61atQIs4vHLL79EpVI59F0Sb29v9uzZw5QpU3j44Ydt5QEBAXz77bcOAvaDDz5g+/btfPHFF7z00ku28hYtWjBkyBDatWvHrFmzyhxXIKgI1x3+IDQ0lHvvvbfS5yclJbFx48brncZNw4Ot/sdWBpB8PgF5wXkyTQNKba9wYpFCm84Do9qx8pcjgMVa9NCY9i77GDqpHX+vvICuwECn+yLsgnDK5TJ6torE39OdIB9Ph3OtkdLlcjkhfkUOwnqjEbPBTPHVYYVKgd7TxK6Tsbay+h3tdxJaceUjJcNirXK2Q8/VOQILKrXCJqRuFLy8vBycuStrjZIkiR9++IEuXbrQsmVLW/mECRN4+eWX2bNnD926dXN6ro+PD4888gjz58/nvffeY+7cuTRq1Ii77rrLTshUBIVCwYcffsiDDz7IV1995dLxOT8/H7CIgpI0aNAAjUYDYAvDUPK80sIoALbwCRX1P7Ly999/k5ubyxtvvMGwYcNs5cOGDeOll17iww8/tGvfsGFD25JgTk4O2dnZyGQynn32WXbu3MmePXtsQioyMtJmBbTSoEEDvvjiC9txRkYGOTk5dmPUq1fPJrSWLFnC2rVrkclkBAUFkZ+fT0ZGBt999x19+/alWbNmgOV+/e9//6NLly68+OKLdv3162cJZnz+/PlK3SOBoDSuW0g1a9aMH374odLn79ix47YSUnKZmf69k0D9CttPdSMzuxJCSpfBS2/2xcdXQ1JiDg+P7YBfgOu4Kj5B7tz7ZCuX9ZEhAQ5LejascaRKLO0p5HKMRnuxo1DJHJYRS/o15eZokctlLp3QQ/298dC7YXJmkTIJIXWz8cgjjzBv3rwq6WvLli1cvnyZ8ePH2y0X9uzZE7lcTkxMjEshBfDUU0+xYMECfvjhB1atWsX06dMdNkdUlBEjRtC1a1dmzpzpcnnP398ShNeZYCu+s7fkXKznXb16lXbt2rmcQ1JSkl37ihIXZ4kP5+zede/e3aFMkiRmzpzJ7NmzHZZtwf46P/74Y8aMGUP9+vW588476dixI3fffTc9e/a0Xe+7777rYCU6deoUzZo1Y/Xq1UyZMoWBAwfy/fff20T477//zsiRI7n//vs5ceIEKpWKbdu2kZOTw5gxYxzupU6nA4qWLQWCqkQ4m9cGSbtBn41aUXZ0ToXcSRt9Fu4aOc9O6c0HXw6hTYd6jm0qgFlyLlB2nbhEWkEBAO6ejk7pJcWOyk1hWyq0UtyHZ8HXe5jy8Apefmg55066dg5Gcm59Ehap0im5W+9WIyYmBoCpU6faxaTq2rUrZrOZpUuXlmqVsTqdP/fcc5V2MnfGxx9/THZ2NjNmzHBa3769xVpcVpgGV+etX7/eZRtJkmz9dujQoUL9W1EqLX/bhYWO3zXOyj7++GNef/11YmNj8fb2pl69eoSHh9tEjnXHH0BYWBjbtm3j6NGjjBw50iZ07rjjDjIyMgCL71TJ8BhWa5Q1HMRXX31lZ8kcNGgQEyZM4Ny5cxw7dgzA5s9lvW/FOXjwIAAdO3as4N0RCMqm0hapTp06cfjw4euKB1OV/dxUXLFsw9aoShdSB8/H0y/yT+eV+izQBDmvqyBmJz7c6bkFnEtMY/QLXQBw93R0PJZMkt1vkMpNgeRm/yZotUhlZxWSuDeD1hFhlIUMi8WqpK3AGpNKW2jg9IlkIqICSrXE3Ww89mnPCp9j0JtY9NYeAJtPVEkefrMzKrcba8mvomRkZLBq1Src3d2dOocXFBSQmprKsmXLbP5Tznj++eeZNm0aAwYMqJSTuTP69OnDvffeyzfffEObNm0c6keMGMGUKVOYM2cOI0eO5I477ihXvyNGjOCVV17h+++/t4mPknz22WecOnWK9u3b07Rp00rNv0mTJgCsW7fOzg8Jihzli/Prr78SEhLCzp07bctqYBGKAwcOtGtrNBpRKpVEREQQERHBo48+yuTJk23Lg9OmTWPq1KlMnTrV6dwKrr3IpaenEx0dbVeXmmoJAWNdArWGcnDmq2Z1gC++K1IgqCoqLaS8vLxKNTfXdD83FYUW07dGWVBqM4VchlrpPMUK2vQqE1ISjkpKbzASGuZNr7ssEaedCSl5Canj5q7ErLK3hBj0ZtatPM4303cwpGuLcs1HhiX+lLJE/2aTmayMAsYOW8jVhBx8/d2Zs3gUjZs6d2i/2ajMLjllOSxP7l5uN5zvlDOcLROBJZ3MTz/9hE6nY8aMGbz++usObS5fvkxUVBQxMTGlCqnHHnus1PrKMmPGDDp06MAHH3zgUFe3bl3efPNN3nvvPfr06cOTTz7J4MGDCQ8Px2g0Eh8fb9uhZ7UOWc97/fXXmTZtGv369WPy5MkMGTIEf39/rly5wg8//MBPP/2EQqGw8zmqKJ07dyY8PJyffvqJoKAgRo0aBVh8k3766SeH9h4eHhQUFHDgwAFkMhlarZZ//vnHqRjq1KkT99xzD3379iUiIoK8vDwWLFgA4DSKe0m6du3KypUrGTNmDB9++CFt27YlLy+PZcuWsWTJEtzd3W3PD+tOxD///NPO8rR9+3Z++OEHRo4caSf8BIKqQuTaq0V81KUsb2FZPlMrXey20mVU2TycCSmD2cxPq8fbLD4eXo4PeaXc/uHsplFiKvG8jr+SyW/fHiq3iAKQy2TXwiaUEFKSmd+W/svVBItjanZmIcsXH+b1aaX7md3KKN3kPPnf3mW2udExmUwO6VisLFiwgPnz56NUKpkwYYLTNuHh4fTv35+NGzdy+vTpGn9gtmvXjlGjRvHzzz87rX/nnXcoLCzkk08+4euvv3aaQqZDhw7MnTvXruy9996zhRb44IMPHISaWq1m/vz59OnTp9JzVygUfPfdd9x///18+eWXfPnll7a6J598krlz59rl9Rs/fjxPPPGEw67tl156ye5cgOzsbGbOnMnMmTPtyjUaTbkE7dNPP82iRYs4fvy4TeAVZ8aMGTYn/kceeYR3332XqVOn4u7uTpMmTdi9ezeffPIJrVu35rvvviv7ZggElUAIqWpm3i549Om30RxzfFP1cssjJSvPbjecVm/kVFwyZrNESnYe6uYuLFK5sXB5DeRdgRaToI6jU2h5cRadqW5DX/wDiwUN9HTcqqxW2qsmtUZJQaae/u2jkctkyOVyjCYT9YP8KjQfuUyGuURqmVNxyTRrVJdZn9kvdd7uQkomk92Q1iY3NzfCw8PtkgC7om7duqVaJ7Kzs8nKymLcuHEu4ymBJTXT6dOn+f3332nWrJktPU1ZDsYKhcLpXMPCwhwSE/v6+hIeHu40AfL06dP5559/MJlMtvhMVmQyGTNmzGDs2LEsWLCAPXv2kJqaikajoVmzZowYMYIHH3zQIcCkXC7nm2++YfTo0cyfP5/Dhw+Tm5tLcHAwvXr14plnnqlQvry6devalsKKY4299Mknn3D27Fnq1avH888/T05ODnPnziUgoGj37eOPP467uzvz588nISGBBg0a8Mwzz9CmTRtWrlxpt8vw4MGDLFiwgC1bthAbG4u/vz8dOnTgueeeK9dSpLe3NwcOHGDOnDls3LiR2NhYZDIZrVq14sknn+Suu4pCdNSvX59169bxzDPP2NLmBAYGMnHiRN5///3by31EUKPIpArmj1i9ejVr167FaDTSr18/Ro0aZWeOtraJi4vj2WefrdLJ3mzk5OTg6+tLbtwBvNZ3cqgv1Kv4dM9KAr2LBMvOYxe5mJRuO/7sgR/pE30S5CowFxNVMgVI15Z1FGoYkwCayuUMm/n4H3gq7S1O8iAFT39QZOmQJIlZE7eVusupcd9Qkq/kkHuhbCd6gH8vJtI2ynErvNFkQu9mxsNUJN6OX06iSa86LI7Z79D+wIXXyjVebaLVarl06RKRkZG27e4CwY2CyWRCoXAU5Pfccw+bNm3i8uXLN0WCYysZGRkYDAZCQkIqtDNT/J0KrFif39nZ2bYQI66okEVq2rRpvPfee7i7u2M0Glm4cCH/+9//+PXXX4mIiLC1i4uLswv4drtjdnfuYK1SmhySBJtKZDPXKA0gU0JIV0gq5kwsFfONMeng6p8QObxyE3TyPaPxsLdAyWQyTGYJpcL1l5LGQ2WLO1UW+Vo9V1KznAophVxOmimfA7tPI5fLUMjlaPUG6uc5j0klEAiujwMHDjBjxgwmTJhAdHQ0SUlJzJ07l02bNtG1a9ebSkQBdhY0gaC6KbeQSkxMZMaMGcTExDB+/HgKCwuZM2cOb7/9Nt27d2fjxo20bt26Oud686LyAjc/y067YijlZpQlTPktGobi66HBaDZbRJXcA3wbg3ud0sfIuVD5+TnJ0efp4+gTZQmT4FoouXuoUKjKJ6QKdHoH0WhFJpORk6clu8A+b9fa5ccI9vUkp0CHWqXAU+2Gp8aNnGwtPr7i7VEgqCwKhYLVq1ezevVqu3JPT0+n/lwCgaCIcgup/fv3c+edd9ocBD09PZk8eTIDBgxg0KBB9O7dm99//73UgHi3NZ71HYQUwMZDZ1Ep5CiupV65t2NT6vgXRUA24wcBoaAuI9ieSedYZiwAoxY0pb+dyeQ4OEp5+akd2jkLk1Acdy83hzhSrjCazC6FFEB+btH1NK0XTKcmDXBTOu874VImPu3KDqsgEAic06lTJ9atW8cPP/xgS2nTqVMnpkyZQvPmztNPCaqJq3/Bqdng0xjav2lx3RDc0JRbSGVmZtqybxenVatW/P3339x77730799fZNV2RWA7yDzuUNyyqY49h4sUitksIS9mIdIoC6FuX4tTeWloSzjrxq6B7WPAkA9tp0DXT1yeKpPLoMQuem9/RwtPyd19p+JSSEzPRqGQo5TLuT+4Xbl3iBlMZkylKLPCfIs/mJtSQZemDV1HXgcSL2fTXAgpgeC6GDRoEIMGDartadw+mHSWl103P7D6cRVchXV9wBok2aSDrh+77kOfC1mnIW49+DW3PGMOTbfUNRwCfRZYfGfzrkBeHAR3EsKsGii3kGrWrBlvvfUWOp3ObissWOJ3/PXXXwwZMoTBgwdz11132ZJZCq7R7nVI3g25l4rKFGpenX4/n3zyL1kpaXgUHEMut89K7q4sgLp3wZW1pfevTS362WyCXRPBcC3K878zIWIYhDqxFkpmMOSA3N5q5RfkGNROktkLn5SsXK6kZtmOvf3dUZbbImUqNeWLtsBAsK8nXZqULqIA0hMql2NMIBAIqgXJDJIEchffh4k7LYIJwCMM+iwE7whYf2+RiAL49xNo9n+QehDy4wEJ9r5avjlcWQs/Otk16xttmVtIV8g5Dyl7odFI6PGNRXC5+YJPFORcAn02BLYtEnoCp5RbSHXp0gU3NzdefvllvvnmG4d6Pz8/Nm3axKhRo1i9ejWTJk2q0one9AS0hIdOQvLfkPIPGHKh8WgaBDTjmx+agUnH/KdfpuR+N7W3l+UXX12G82Rxi1TKXsubTXGOfeVcSP07E3NhEHja9x8Q4glph2H/25Ydgl0+Qirxt6Qo4d+lcVeWeyu+h1qFl7vrNyNtoYHwIH+70BCuyEor3y5BgUAgqDSFqXBqjsXC0+wJkCstggTshcaV3+GPwZafIx+Efkss511eDQnbQCoRG7DgKqwvJYTL0iZVex3Z5yz/5hRL4HxhieXjCoU7tHwWOr1veS4kboewOy0fQfmFlFwu57PPPmPz5s0cP36cVq0ck+C6u7uzYsUK3n33XWGRcoZSA/XusnxKolBjVtVzWGLzqN/O8keqLiO0QXGLlDPrVewqi9gqGQ393GIkpjg0DwxWw5phRUuKOecoMM8kISkbk8mM0WwmM89ewKjVynKnIqnj71NqkE4/d3f8PB1j9TijMMdJYmeBQCCoCjJOwD8vQ/zGYmXHILgL7LwWIFamsLyoRj0Ef79Q1O7Scph3C4RrNBXC0U8tn5J0+wqaPwmbR1iWDzu+a7kPtxEV+h8eMWJEmbmKFAoFH374Ybn6KywsdBrY7rbFzY+SJin38K6WH8pwGLezSF12zI+FWQ9nfoC2rxQrM0HmcbLz0nBzq4NSIUcul5OeW4Bad87eLyvrNHnmFHafdB44UeOusgSHrKIo2n3aNCp3W6P+5knQW8GwbQKBoAZx+Ps0m2C5o9GAk98BxSKlSyZLeJriIWpuF/a8aPlY2fIw9F5gWaI05kPTx0DhBmYjKG/N5/11SeULFy7g6+tbrujFJbl48SLvv/8+P/744/VM4ZZCUno7lLlHXTOdBnVwGkLBhjbNYmb+91PIPOm8TfxGeyGVHwdAXMolNh21j6D+bq5jahq1ORWnQaewLOsBuLld39vXuYRUEjNy6N48ApWLXXolMRmKvvyyMgo4cfQq0c1CCKnjeD9rC2uwQ4PBIF4eBILawlAAuRfApAePUMtu6mIYDQbQpqNcfAeYcmppkrcAO4ulcvr7+aKfA9vD8GtBldMOWVZbgh2DVd9sXNdTLy4ujocffpjt27eXGfmzOBcvXqRPnz5i+a8Eknso5NqXqQKu5R9TesCAVXBomkVQtXq+yFkRLLs//nwczsx3PUBxR3eA7LMAaA2O6V9s6+jFUBsSAMsXj1wmQ6mQYzCZkSQJjcbSh0rjKH4uJWWgUiqoH+TrUFeSE1eSycwrpEvThqgon5AyGy3OmYnx2fxn+I9kZhTg5a3m+yWjiW4WUq4+qhuVSoVarSY7Oxtvb+8KRVsWCARVRPbZIh+lgiTQZVmsSdeyRuRkSShyz6Ew5bruQ1B50g87X+r8P63l/0HpuMnpZuC6F28PHTrE4MGD2bhxY7netK0iKi4ujjvuuON6h7+l8PDxJjulwK7M7oFbtzfU3Wr5WZfl2EFpIgosS3VmU9FOkiyLkBrTeRcf/FG0ZNs80gDZcQ6nqw3xWIVU/SBf+rWLBsBsNpNvuBauQO34KxWXloWvh6ZMIXU5OZPMvELqBvjg7uZE3Lni2sregm/3kJlhuX95uTp+/H4f0z8fXP5+qpmgoCASEhKIj4/H19cXlUolBJVAUN0Y8qEw2WKFKunobbAE/ZUkyNdDTnYWYVcXInOagVRQbcQUC7fT+iVoOQl8yu/eUdtcl5Bq3749bdq04a+//mLEiBGsXr0alcr1A/DChQv07duXuLg4unbtyvfff389w99y3DemFb+8v9f2cNXhuLxmw83Xkjqm5BdDcdzrQGFS0bHZYPGfOvqpxYKV/i8A/Zv9y88HenAxrQ4ebjpeGvwv5GQ6dKdRFgXJLB6SQC6X22JfqTWOv1J3tirb8ng8Pon9pyzirUvTiqWjsLo1/Lb0X7vyDatO3FBCymq1TUtLIyEhoZZnIxDcZJiNlu87udu1KMJOkMz2/qLG8u7olZCZCvBLXY9v9p9lNxdUH8e+tHyCu1gEVeMxliVAyWTJOXsDcl1CytfXl02bNtGrVy82bNjAo48+yi+//OKQwRwsIqpPnz7Ex8fTtWtXNm7caJclXACBYV6EdwziyqF0zEjc85+WrhvLZJYdeMWFkl1n7Sxr0Qu8wVQs1crmBxyaeqr1/DhuFqeT61LXN5OQ+uGQ7Sik1Moi0VY3oMRS7jXDijOLVHkwKCzLc3KZDH+vivkQmY1mpr62vlLj1jQ+Pj74+PhgMBgwmW4eJ3mBoFa5sgH2vWT52acR9FsGbp6gz4ED71r8Qjt/ANvHVXIACZUhHYU5v8qmLLhOUvfBjn2w4z9FZcFdYOhuS+iJG4jrnk1oaChbtmyhZ8+eLFu2DB8fHwdLkxBR5WfIk20x6EyWZL1l5a0rTUh1/5/ll807whL5tgw0KgPt6l+2HGQcc9pG4daE/u2bOF+iuyak1J4q4tOyqB/kV+aYxWnSKoQjxxPK7WBeHD9PdzZtP2N5aZFAcS3RsUwG2zaewc/fgw5dGlS43+pEpVKVar0VCATXyDwFO4q9AKZeht87wUOn4MfQovI/etf83CpJSq4POVp3GgUlO411aTTL2XOxCT6aAlrXiyMj3xNvjRa10siOcy34Nz6cHo3O0KnhRaf9F/p0ZcdhL3xkiXQfPhSZb7RlNUIyQ6OHocF9kLwHNg6Bev0t0c6vrKvmq64CUvfBPBX0/BZaPF3bs7FRJbKuYcOGbNmyhV69ejFv3jz8/Pz49FNLvAkhoipOeYNa4h4MjoYj8I6EOj2u/RxRLiFVHuQKH5d+TjKF5dvAy8eNzYfPMaRrC4J8PMvVr7uPEkWwRVSoyxkZvSQDOjQFLNuXrUujVzNy+OrtbegMRkY+1YkJTzsGJE2Iy+LU8STadqhHcOiNs8tPILjtyb1ssUZc3emkLhbm35y7X7ecbsW76x9Fb4CejU7Rr+lxlh3qxh2RZxnd9TCH4xrw6oqHy+xn0b7evDm5Lvc1P4i7soCTew8jl0msTZrI0t/Sbe1C/vJi9GOt2LhGwbBH2jKod0vycnV4hdxD7tAc0lPziW4aiOL0LEtYh0YjLcE3E7bAttH2g4Z2twSVrm12PWP5KD3ggUPg17RWpyOTqjCwzb///kufPn3Iyspi+vTpjBo16rYWUTk5Ofj6+pKdnV2hXY3lZsvDcPFXx/L2b0Pna/mW/nrakgCzCvjx8EvkSs59jgpkBqZ8dw96nZHuLT6nb5tGRISWEfsKuJJ0lndHruTXy++RfDKn1Gjn18OJxCS+XWP/pXDqWBKTxi5BiQKzwswPK8dRP7yM5NACgaDqMemhMAXcQywxh/TZ8INfbc/qukjL8yIpx58d+jc5erKQ/nWW0TIsjnE/PlfbUyuTrxc8RNeekWgLDSxZeJBLF9J5cHQ72nasb/G1PfwRJGwh9dy/uLvp8VLryu60Omn3Jqi8LD50Sg/LpqrQbhDWq9JdVuT5XaULjW3btmX9+vX079+fd955hy+++ILMzMzbUkTVCCWjlFsJLyZ2vMOrbDiVXOcQed1Kfq4lurg1snl5lujOxKeSl30Uz5zNqOWvVpuIAgj393dICP3T7H3c264Znho3cgq0LJi1h3dmDqy2OQgEtz0J2+D3fkXHdfvCHZ9bLE8uXApqkzydml0XmlHPN4PW9YrtZG5wL8T9wZ6L0Xy5bTBKhYmHO+zhSN69+GgPsedSEy6lW5cdswE4dGxYjc+/sjw3wfEFfcOqEwCMGN2OgcP+j1W7O7N2+TE8PVV8+EE7eg68oyh8wZUNsP8N24amaufIR87L710HDas/EXeVe2x169aN1atXM2jQICGiqhtNsGOZ2h+CigU484pwfX5YH0CyKPkrv5c5nJtC61JI+WnMUJiCzD0EpUqOSlG2kFLIZbi7Wd5kPNV5ZFB+q93eM1eoG+BDg2C/crX3cleTm6PF169oOSDzYh51w+sA4OOh4fBhsZNOIKg2Uvbaiyiw5Gxb2aF25lMGOqOSR394jvgsywvr3U2P0rRFCP3+Mw6Duj5zf9/N1g1nbO2LQshU3gpyM7Di5yOs+PmI7Tg/38CLL+1nyK5C7rm/BfNn7SGsni8vvL6bgKBr7h1mgyUMReZJcPOG5W1qZrJ/DIbH8i2+Yaqy87YWR3F1S7nbXpeQOnbsGG+88YbTuuDgYOLj43Fzc2PMmDFO27Ru3ZoZM2ZczxRub5xZpOrdbZ9x3DvC9fnNn4TGoyw/LwwCXbpjm3r9IWEzcE1IuSBUsxsOboYeX6NSKXArh0VKoZDjrrLEn/L2clRoWr0RjZNI6WcTUjl5JZmG5RRRAAlp2bwycSVf//CwLXhoy2siykr7RvXK3Z9AIChBfgL8+YQlGW70OIufjdqvqH7VjRU3UJLgaEJD3JQmwgNS2XiyLXsuNSEt35tHO//Fq6vG2rXfcqYNW87ArN821dKMb2zWrjjO2hXHrx3F8ftvxxn+SFvuvLsxXXpE4Obmi6xOd0v1fRssWTiUHhaBXTxXbFUzv5ivbs/Z0OKpss85MQvP3c+We4jrElLp6en8/nvploy//vrLZV1eXt71DC9QO/FBqn+P/XFpQqq4g54m0LmQin7UJqTUChcxWSQt/RsvhZNXIG49KvM43MrhNK6Qy3C7ZpHy9cwHivyTTGYzGw+dYUiX5g7hNKzWrgDv8kfBrRfkizpHyYaVJxg+ul25zxMIBOVk/zsQt8Hy84F3LB+Avoss3yM3AJIEB+Oi2HyqDSuOuBZ2ryZE1NykbmF+W/qvXXy/16cNoNddjQgI7Y9q8L1FDQuSIe0g7Hu9epd4d020fKw8rrfEpipIAqWnxVoGUAERBdcppKKioq7LotSwYcUCLwpK4OZkKaz+APtj91BQuFuyd5fEt0nRz5pA61J+ETI5+DUrGk7paJHKztzCsz1j8FGnWApyL9G+USQearcyp++hdsPtWtBRPw/7vFYKuZycfC0Ltx6kdUQYnaKLcmIpFXLc1SrUqor9+gb5eLJh0fFShdQvCw4wasLNn/tJIKhxzi5wXr59rOVTG+ihcBmcO+dDanRdtvu05Y9T7at8GDdTAaGF58hVBZOlrlvl/d8qfPzuJj5+t+g4snEgXy94GA9PX3TqXgQ+8G+RH6tktvhYVefS7zwnz6l6d1e4m+sSUg0bNuT111+/ni4E10NoD0uUX7PF0ZuQO8CrRLwkmQwihsGFX+zL5W72a8bOrFvuoXZ+WJoSFimdwYiv8lCRiLpG3aBQSqKUZWOU7P3kgnw8Mel6AisJ9HC0hrmpFBh1ZowlAlcqFXL8PCu39blteOlfct9/vZuHx3VAoSgjhpdAcDtgyLdsiZfJLXF7lO6WpZjzSyCwLTT5j+U7Jj+xtmdqw2BSkLzXmyM7g4lLUZBZeAk3KYfcNBV7Q4NRKQsxKFx/f3gZ0mmRuY0ctxDO+XTDJHd82KpMBSgkI0rJQIAujq4pv6LAhBk5u0NHk+DZqjov8Zbh0vl0Bvf6zq7sied74Oam4KFHO+AV1B4eOQs7xkPeZUtAztjfqndSCeX3jbJyY4UHFVQMTQDctdiSyFgTBD1mOW93x6eOQqpOT/tjdaDjee51LJYq63Aq+zyASoWcAE/H5VmT2T6iRmp2HhM7PsMOzRrSj9snA1XKLVYuL2WqXQwoALVKSYHOgNFkdhjXv5JCKiWr9OXknGwt2kIDnl7Vt4NQILhp+GMwXN1h+TlhK/T4H6zuYUnXAZZYPpjBVLPb3/NzQKkEtQpQAHIwp8mYtv5B1iW2p3fSfMIKLYnXrTLI25jOPQnfAFCg8MHDZLGCX/Fszf6g4Rjk7rTM2krrzKIHabuMDWwLe5xUTaQl/54k0S1lCQ0KTjidlxwzd6QsY0WkEFKV5fv/7QZg3td/883Ch9m64TLNW89h4OiWRc+HrDOwrFkpvdQsVSqksrOzycvLo1494bRbY0Q9aPmUhmc96LcUtj5if15xNE6ElEcdUHnbcvppVPYWKYVcjo9bgcNpJYXPsdgkgu5MISjcmwNbLhNZp8j6pZJb+lSYMjGYzbgV2+3n46Ghx4BGnPz7ql1/SoWCDm3rQyW+u49fTkKSJExGs9N6GaDXmfCs2AYPgeDWI2FrkYgCiFsPS0qkYnLmMlBN6LPgn/kK9p2xiDgFctQKL/JNeaBQcdlvAP/4diRUe56wwvMAyJVKglq3Ju3YMczGohRXVhEF0DD/GA3zXfvl3HV1XoXmqZL0NM/aQXjuYdTmAo759yfJowkqs5ZsVShOQ5k74blX+/D1zB0VGvtWQqcz8sTIn23HRqOZoQ9d2+3n1xSevPbCvudlOPZFLcywiCoVUocPH6Zv377UqVOHzp072z6dOnUiKKhoh9ny5cvZvHkzr7zyCtHR0VU5BYErGj1sCXR3cQUEd4JmT9jXu7JIyWQWy1dhChqlo2jSuDlGAzeXsEjJZTKUCjMqNwUqpf2SmZt1uVCXiUwjA0NRXbCPJ3f0iqRri4Zc3Fu0q8Pfyx1ffw+ykxznUxaFOgMfvbWRpyc736LsplKi17tOBB1/OZOjhxNp076uCN4puHUx5MPvFfcVuR70uZC7D1S+kOemYtGBO/nnUihuoT682n8127+4SvH4KybMFJhyLNmpTDoi0teiKThNnWuWKLlSydhDhwhu3ZrUo0dZ1LGjnZiqTtpm/GH7uUvaStvPOW4hbKr7DGpTPvlKfz5IWco9+f8i8/ODc+dINanRaY2275ZHxnVg144LhNbxpnV7i4EiLjaTv/+8SNMWocjlMib9ZymFBQYef7Y79w5twaVz6Ux5pmj5a87Po5j39d/s32NJAdYg3J+4y85SYhTR954mNIzwZ+GcvVV1S66b6a9vYPrrG+jSI5x3Px5InbrXfIS7fQ6dP6zVSPdVKqQCAgJo1KgRFy5cYO3ataxdu9ZWFxkZSefOnWnfvj179uxhzZo1yGQyZs+umqjb1UV+fj4vvfQSu3btokOHDnzzzTf4+fnV9rQqR8Qwy8cZrixSYBFZhSn4qHMdmgR5O/4KmUoEy1dccx5Uu8lRlogv5aa4JoZ0GYRG+JBxrmjprXVkGFKhmeadwuyEFECXgRFsnn8SAKPJRKHeiHc5Anre0awhO34/S4euDVi45QD/udvesTy6bhCFeQan5547ncL4EYvQaY24e6iY/+ujRDcLKXNMgeCGw2wEZPahUoqzdVSVDbU3tjF/nW9O67qXuafFUVt53lnYNhvydRAcBUfOFj/LgIKt9ABIge3l3MhlFVEAQa1aEdy6NQDBbdoQ1KoVKUeOXO/lXBc++hQejH3fdnz82scvKwufunVpNGwYBe7uuE2eTEjbtmjcVdx9n/0SVoMIfx6J6Gg7/uvYZLv6iKhANu59liMH4mnUJIiIqEDadqzPgX8u4+fvQbOW9j6s+3bH8sy4pbbj2YtH0ekOy0aw517tQ3paPpJZwstHzcpfjuDj607LNnV46J6YqrkpFWTf7ss2v6qmLUOZ+ukgGjcNtliodFmWmFWaIFjf32JVrQGqVEi1adOG8+fPk5GRwf79++0+ly5d4tKlSyxbtszWPiUlpZTebgzeeecd8vLy+PXXX5k1axYvvfQSCxa42J1yM+PM2dwqpK6JLC+1o39RA7/zDmUOFqlrQkol16Mq4cRts3LpMml9Rz12njtjV6+UyWnYMhC/Oh5kXbNAdRoYQZMuddB4qUi5nMtrr6/m3g7ly7UU6ONJ//ZNePfldXi7a1iz9wT3d21pq+/cpAF5Oc7jZf3v4x3odUbUKiXaQgM/z9/PezOrP2quQFClFE8b5Szyc9xGuLLW8bxyYDTAX4vg6kVo3AE8utZj0tLHAVhysAeyXBhw6igcgj90EHvtvISzLrusNGnHj5N67JjNIpV2/HjZJ9USWUCWwcCVXy0Rxff++CP3L1lCk0ceKfU8VwQGedLv3qLvRKVSzh09I5227dIjgq8XPMT+PVfo0j3cJqKK92Vl9ITOtp8PXHiN35YcYdG8/TRpHkz33lEs/+kwJ48lVWrOleHMiWRGDpxPVHQQcxaPwj/Qr6hy0DVft/wEWFzf6flVRZXm2iuNixcvsmzZMmbMmEFhYSEjRoxg6tSpNGnSpOyTr5PExEQKCgpo3Lix03pJkoiLi8PT05PAQHvLTKNGjfjjjz+Ijo4mKyuLyMhIMjNLN4taqfZce1VJwlZHU36/pZYlweVtIcPyJrn23Ayu5HYBoI7ncUY0fcGhq3e2/kYd/6Lr3XPqMj+NeYzjrQ6ydnYSPh4aW12U52zua/oraIIpHB7P/Cm77Pq6d2IrGrULQVdo5Ny+JNx93IhqF2znlP6fYQvpUqfEbsUySM7KJdTPeZLibo82okNPx9Q6PZt9Tv8OTQjy8SQtJ5/Nh86y6/TLFRpXIKhRzAbLcr5cBZHDLRkMNt5v32bEETg4FXIvWVKfHPm4YmMUQtIKy271RGD77qKqPM8G5Jg9qFtoeUHSyj3RmPOv65IqglypJKhVK9KOH7/uZT1X/lbVSTOg/siRNH/qKdQ9eoDKEkwYsxmMRpKPHePI7Nn4RkTQ6ZVXUKprf5OM0WjmlYkr2bX9Qo2P3euuRox7sisymYy2HevZPScAyy/pmjshebfzDkoge4qaz7VXGlFRUbz++us88MADdOzYEZlMVu3+UatXr+abb75h165deHp6kpaW5tBm9+7djB07lqysLPLz8+nbty+//PIL/v6WNeqkpCTq17eoWT8/P8xmM1qtFo1G49DXTY1TH6lrJmBreAXg3qipHEsdhlmS0zpktdOuSkpzq0VKacqxE1EA7qpry4W6DNy9VAQ28CI9zmL5kitlNGhmsZSp3ZW06u38reKhB9pz+W/H/9vScCWiALSFRV+S506ncOZkCp27NaRJvWCCfCxvZ0E+njSu6yLXoUBwo7BpRNnWpRXtin5OP1J620TQ/QnKuqDoBchgx1Q44OLd0is/juL7NmpSRAGYjcYqWc6rCn+rygix08DpJUvYsmQJbYE7sVivtgNJQPFetJmZ9PnsswrNqTpQKuV8Ne9BTCaz7bs/LSWPfw8m8Ppzzp8ZVcVf2y7w17YiAbf//KtIEmgLDag1SktYm6HXXtYNeXDpN7i8BnwbO7xAGIN7AOUTXDUe/qBJkyY899xzzJgxgwEDBjB+/PhqG2v9+vW88sorDBw4kA8//NChPjs7m6FDh/Kf//yHTz/9lOzsbPr06cOTTz7Jr9dMrP7+/mRmZuLu7o5er8dsNt96IgqKIro6K4seC/vfAkCl0NKhzpJSuwoLsFfvXZtaTMW65DjAfkenRnnN2VwygSGPARNa8teys+gKjHQeHImbpuxf0cZNgisspEpDr7V8Pe3bHctzj/2KyWjGx1fDiC6t7dp1blIxK5hAcN1IEpxZAEl/QoP7oJGTpZ+LK+DCEnAPrvQSndOhs2H9NDglgTcw/Cr43ulaRN1KXK+/VVUIsX+vfVxx4PPPbwghZaV4LL7gUG/uHtiMH+r68M9fl0iIz0Ymg4BAT36Y/U+1zaFz45m2n5u3qsMXcx8gOPTac03lBU3GWj4AXeyDixfk5ADlyxFcK3GkOne2rLOuW7euWoXUnDlzADh16pTT+uXLl5OXl8f777+PXC7H39+f119/nUcffZTU1FSCg4Pp1asX8+fP5+2332b+/Pn06NHD5Xg6nQ6drmhPfk5Ojsu2Nxxe4ZZdeoXX1rc1QeB/LRZK9Dg4Pc9i+r8OZCl/ASPtylTyYjEMdOkE1I1g6IvXIg9r0y1LDpIJWr1o2T3ohOhOoRzfkUByrOV+n09MIytfS4uGIeWKsF4S3TUhNXPqZtQKBb4+nqTlVHyHoEBQ5VxaDn/+n+Xnswstf6f1iiUCTj0IW8oIh1JJrvxiEVEAucCP22DwtmoZ6objev2tasrxXX/xIm6RkeUOsVDTtGpXl1bt7IMiKxQyYmbtqfaxTx1P4pP3NvP+p4Pw8q7aJdAqFVIFBQVIkoSnp2fZjQG9Xl92o2pk//79tGjRAm/vImtM9+7dMZvNHD58mAEDBvDhhx9y//33M3PmTIKDg1mzZo3L/mbMmMHUqVNrYupVj1wJveZY8hBJEvScZQmXAOBVHx48Cis7QrYTr9BWL1icV8sIyqfK+IeSQsrDrVgcmmNfQfevIOM4/PMKxG8sqovbCMP+cfoFoVDKGfJ8W+JOZbJg7j/8dcIi+I7FXqXfnU1oqC7fW4UVg86yxTovVcfw7q1xUyrIK6zZgIMCgVO2jrY//vNJGFXMF+Uf+x1cVYIezP/Cr4cdq9ZVQfe14XtUUcxGI4s6dKi4v5XSsnqRdvocVw8dJahlS9KOHyft9DlLnVFriRovOY9r5wyZTxiywEaYk09CQYZd3f8aNWJMv34Ynn0W/+RkvHv1ghYtyt13bfD05Dvp3C2c/Hw9Pfo0QqGQ8e/BBF57dhXpqVW7FLxj8zn6tPvKdhzRKIB2Hevz9Mt32jnVXzqfxmvPlj+CepUKqX379nH33XfTokULunTpYosj1bp1a1RWJzmwLZsNGlS7O57S09Pt4lsBtmOrP1VUVBTHjx8nJyenTIezN954g8mTi77IcnJyaNDgJlr+ibjf8nGGygvCh8LRT+3Lrdm0fZvAnhdBMuOmyUOvLfKM8FVavoEb+F2FePvTG/glFB0c/y/4RsPZHyF1n33D1H2QdwW8HZ3AAdQeKhp3DKHnwEZs3GKxQKrVSp55tzfrPjlS1pXbYV3a6xRdHzelZXu4l5PQCuaa2achEFjCFaQeAKnEAzz3ov3x1T8rPURaIqRfBY0HhEWCmwaIhbyvYXmhNxY7VNVSXbGe9gQ/wj25B2iuvYBZ4cM83+52sZ1KMuCee2g5dSqrR47kYmws3vXqMbRuXWL37ycRyAHSXPhbhQMPbduGtH0HCfNWcxQJ/dWjFCg1ZP7f77Z2K+emATsBUP/H8YXceG4rhj3fgaHo5VKuVODX+V5ypbqYU88j5SWjHmz5DpYKs9H99gxSrv0uucVbt8LWom3/YyZPJuzzz8u+abVIp2723+vtOtVn4z+WxMFarYHTx5N5/JHFVT5u7IUMYi9ksGpZUVgOdw8VhQXOQ+C4okqFlKenJ+7u7hw7doxjx44RE2OJM6HRaGjXrh3R0dGcO3eO/fv38+qrr/Lkk09W5fAVRqFQOFjFrEtzSqX9rSnPrju1Wo36Btg1UW14OolYr74WlLLlM9DIYm269zKs+eoISAB67ov+HwBebmkEak6Rrm0OQJD3ZbzUGdfaXaO0rNs5F1wKKSv3DW2Jl7eak8eSuOueJoRHOl8OLA2j3vJ2WNLXy6GdufxvkQJBpTGbYG0f1zuN5l7fMk5+NqydDfElNNkALM7MltW7qhdRUPElr0te7Yn16kCKexSSTEHd/JOEaC/S3WxiV4sh5HgHMnxkO5Y+2Bzle+/AqVPwxBPUbXUna/67kYYbfmGY+2VOjXucrQs2oIi+Cx+znnofPISiY3OGX7xIQUoKKm9vLv6bhfFUEt19sgj0k3HhiVdIVgXi76MmtlU7zv/1J/76q3R/7WU2XwrhbEpvuL93pe+FMrofyuh+DuWFOH9Qy9x90YxeTGHMQDDaW8wVzQej6v4M5qvH+HnBt0hffIFHSAij//4bj5AQVF5ejjvablA0GhXtOtXnwIXXMJsl/vfxdn6K2V9t41VUREE1hD+QJInz589z+PBh2+fIkSMkJyfbtXNzc6NJkya0bNmSFi1a0LJlS1q2bEnjxo0dRMz18tVXX/HBBx847NqbMmUKa9as4cyZothFJ0+epGXLluzevZvu3btf17g3VfiD8nBxOWx5yL5s4Cao39+haVp8HimbvqZu3jf4aYoSmuqMnvybMgKAtkO7o1YUlH854s550Oz/HMsvr4O/XwBtqiVzd6vnoW4fW/WsiRVz5EiW8ujzQFNO/ZZQart8nZ5XF9xbob4FggpzcUXV+T0dBv4A/ICRkGuA798Fs1Q7S2ylWaTO+dzBwcD7aZ++jtDCC8R6d+CMby8kmX0sunmhe2n38hgYMaLM8c7tT2ZTjPM8edXJhJk9Uakt1m2DzsSCV3eVcUbVov31CaTCTMsSolGHm18dRu/+G99w+4DCSjf5TSGwUpJyGd5vrs2ftbo4ePH12gl/YA1rEB0dzcMPP2wrv3r1qk1UWQXWiRMnOF7Caa9v375s21YzHox9+vThs88+4+LFi0RFRQHw+++/4+PjQ4cOHWpkDjcVzixSbn5OmwbV9yKoRToctc8Kr1bm06Xuj5aD4DFQfwAcfB8M5XDMz7noWGbUws7HLCIKLJnBY3+DPguhyTgAGrQMIO5EhuO5Ljh3MpX1W04xoX/nUtsZjCa++HAroWE+PDKuI8oS6W8EgkqhTbe8GGSdsiynH3zv+vs0A+eA4okkjsA/16pqK52KK9+j0749UbQdzOaVz+F5dRiqjFTunLgXSWexAg8Y3JyP/uvohmDUm0g8n4W7txu56Vq8AzQEN7T4wJ7dn8TmmJPVfk3OUKkVNiFVXqpSfGke+t6hbOkXpwD7jVjdRzSmbb8GtrAFNyohdbxZueUJ9u2+TFZmAfO++Zv8vNrzua6xXXthYWGEhYUxcOBAW1lubq6dsDpy5EiVWqPi4+PJy8sjJSUFk8nE6dOnAYvfk5ubG/fddx/dunVj1KhRfPLJJyQmJjJt2jTeeuutWzPEwfXi7iQdiqKUXXGqMqxwPo0s/wZ3hMTtZY9f0h8EIO1AkYgqzj+v2IRU9+GNWZ90lNx05xHLS+LrqaFZ/dJTv+gMRtbtPYnhb8sXe0pSLi+9eVe5+hcISuXQNDh/zR8k7VC5TzMb4Z/vIPEKRDSD9n6giAUU2J6XWiAfWAUUj1pQrbvKlBpkfg0tvjw6+xemJPdo6hSes42ll2vYVHcSRu86zPt2LP6BHhDYHGjOXyfvJD0tH7PJXLSFvRh6rZEFr+6yLc1XJbVhUaqM+Lpe/l5xnr9XnOeeEXIC23TBw8eNnNQCrl7Ipl4TP3yCPIAbw3IVGubDkActv7Njn+jK1YRsvv38TzasrnmxXCvhD6x4e3vTq1cvevUqSiBblSuNH374Idu3Wx7QoaGhDBs2DIC1a9cSHR2NXC7n999/Z+rUqUyePBkPDw9mzpzJ008/XWVzuKXwigDPBpAfZznWBIFfc9ftNaUErFR5g9e1VASa4PKN78wiVXDVedti4iqovhdjP+jG8k8OkBJbtq9H64gwp+UpWXlsPHSGQG8PlAoFBlPRF/bimP3Ua+DHiNHt7OKnGPSWHYAqt5r9QhTcQEiS5SVAHVDkU+gMk84S7uP4/yo1zLFF8Pc1A3/sPthRgXOrLZ2Kuz/q+79C7lcfvTYf85rnkTJjATgYOIRzvj3w0SfTu4MH2TJf+o69lw4qBU2ahZB4NJN125OIbBuE2kPFkS1XSL5kEWJ3jmxC8qUczuy1OFqr1Arbbtvq4HpFTfG5Vec8yyv4ymq3cYUZVriO7/T4F71Qe6hc1tcGYfV8mf7FEKZ/MYSrCdlcOJvGJ+9t4mpC9YchqrSQMpvN6PV65HI5bm4Vj9Xjqp+qVLnfffddmW38/f356quvqmxMgFmzZjFr1ixMpur7g6kV5AroPc+Sq0syQff/lW6RanAfIMPem/waHd6xpK0A55YuZ+RetKSIT/nHskvQOxxySolttf4eaD0ZGtyDTCajefe6pMSecd2+DAp0eowmM8lZjjkHAWa+v5ncHC3/N8niWzf/2z3M/vIvvH3UTPt8MD36NKr02IKbFMkMm4Zboie7+UL/FfZxn4qzshNkVl7AbL6OuIaV3t5fBpe7z6CZnyUjgZvGk4wer1Ow+UOa+rtTf8CD9Gwdzn0DWpBwOpM6Ub5Ed7JkU9jz23n+3WJ5Ybt8PN2h3z+X2IdhqU5xUhVUxoJVGfFVXsF3vcJw3uS/APCv44G7jxt3PhJts1ZZqU2rVVg9X8Lq+dKzr8UocuLoVf4z/MdqG6/SQurPP/+kb9++9O7dmx07dlR6AlXVz43EpEmTmDRpks3Z/Jai/gD7uDWl4RMJgzbD6RjwbwntXoesk6DQWMIcWHFlkfJtCtnFhI82DZY2KQoa2v1/pT944jdB4g54NAE0QTTr4M75nbkkJBQtC8gVMnQYUZnK/lJRyC2WprAAb1qFh1GoN3DwXBxagxF/Lw8KdHq+++IvmrUMRSaTMfvLvzCbJbKztHw1Y7sQUrcjCdssIgpAnw17X4UHDlqsVOlH4PJaUHpA49GVElGZKbDxO4hPLLttWVxPOhUTChTYP+wT3ZvSrLF9GrCAetG8nH6evDQtPdO0yGSw4dtjABzdFs+meVXnCF6Vy3E1ZVEqTk07pDujrHuYmVRAZlIBS6Y77qJ78r+9a3xp0hUt24Rx4MJrAKSn5TPz/c1s3VD5l+qS1OrSnuA2oF4/+zfwgNaObVxZpHp9B+v6YWfRKiwWM+Xv58se36yHk99B+7dQ/t6doSFnOavsx9G8ieQbQ+k8OJJNy09BOb4bFXI5apWCfu2iUSksXxBuSgUaNyWhft4YjCa2HjnHC/+3HIBgX0/aRtZFbzRx8Fw8Op0RtVr8yd1WnCxhFU87ZAk4u7zE38HeKRXuOu0y/PDRdcytCkh0b0qcZyvivNrgr0ug39W5ALQfMpyCJuPBifE27ngGW3446dRQXZVUpY/RjSBqykN5BV95213PPZz7wk6CGnhxz+Ot8Av1KPuEGiIwyJNPvhlGTraWrRtO4+Wt5o3nXQfaLg/X/a0eFxfHZ9eR3+fChZrPEC24wXAmpGQKCO1h2SmYH+9YXxFyLlgeYNlnkMmgaeAWmgZugZEXwKceG7acguSyu1EoZETVCbSJKIDwkCKfF5VSQcfo+qzbdwqFXEb/9k1Qqyx/YmqVgrSUPOo18Lu+axHcXJQMoAmOIqqC7L7QhK3HW+GxfuV19VNezvncwcGgYbZjmWTCT3+VAqUfOkVR4F19SAuCJ+2meZs6tL8nHP38k1w84rgRZMuC2tk5dztQXsFXlcKwNKtVWlwei9+zrDlHdw6l33+ao7hBdjf7+GoYPrIdAHfd25TjRxLZvvEsCqWcB0a1w9tXjq/v6+Xq67qF1MWLF5kypeJvUwKBDWdLe77RFv8rn6jrF1ImHeTGOpafmgtdP0apLt8ftlIuJyK09ACfwb5eeGncqBPgYxNRAPWD/Ei5miuE1O2GueLB/QD0WbDjS0hMhywDhMuhnRl0DeryjjSIgfFfVu08XZClCrWJqEAfD+oF+pKSlUdSZtHLxOOTujF6QmcWvbGHvDgt++Ni2f97bI3MrzRqejnurnHNANj2o2V3eL/xzWnatQ7ZKYXI5DJ8g90xm8wc3R5PfpaO1n3q4+7jRn6WDm9/DXlZWkwGiYC6RalKdIVGLhxKwTtAQ93Gflw+kc6xHfHEn66ZTNHluYfltVqd25/Muf2WN1Z3H0smijuGNcLtBrDSKxRy2nasT9uO9W1lFcmVW+kraNiwIa+99lplT3fAGsdJcBvizCLlHVX073WkvQDgwhLLpyRpBwDQao2oyvGnoFIqCPQpO4/kQ73aOi1PisuGLjdRyiBBxcmLg9T9ENwZChIhbkOlutn3PRwttop9wQwXgLiCznTJW141cy0HW+tZnHUDvD0Y3Lk58mt+gpsPnyXQ24MOjetjOmtk0RvVn3S2olSV1eWhNzoR1MCbKyfSiT+TiV+IB5FtgvD0c57Fonl3+6S8xZe15Ao57e5uaF8fYqn3DXZc/lK7K2nRo6i/qHbBRLULtu1uv3QkjU3zj2MySLi5KzDoTUgmaHpHKB3vi8DD2424oyls/KFy/kDVtaRZmGPg2PYEjm1PoF4TPxp3DKFOI1/863jeMBarilBpIRUVFcXHH39clXMR3K44s0hZQyf4VKPA1lne6vJydHiX40/Bx+P6YoulXnW+209wEyOZ4cp6iy9e0i44VnlLkSEREr8F3zD457zzNg3SV1e6/4pyxqcHBrnld75TdH2biALo375JlYxRWxG/2/VvSMKZTFKvFIVD6To0ila96pEcm8OB9bF4+LjR8+FovAMs9yCidRARrUsJ6VKDWHfDRbUPZvzHPTHqzXj5qynI0WPUm/AJcre1bdS1LsHbEm3X6mYqQC8pQKlGMpuQya/Pj+x6LX8JZ7NIOJtlO27cMQSfIHfa9W+A+/+3d9/hUVXpA8e/UzLJpFcCIZAeekmA0KQpVVBARVGKWNaGuvZ1/dnW1XVXV2XXxd7FXkAFAUF6R0LvLSGBEJKQXmaSmfn9MWmTmUlmkkmD9/M8eZg599xzzwCTeefcc97j3fiMAC2p9cfUhHD3ty7zqhxi9e5qfcxVdHkAlBXpQd3wSFNT5WWWOH3O0UOZ7N99joFDuxIZHdQMvRJNsuFPcPSjRp9uMsFvh/ty4mgI6hXrKMQAWUrM+cZblskrFJVPCMYLRzkcNA5lvxvpj5JDZzLpHNQ8q4+bM+lkTGIHMk7kUlJgfXu1/1Vd6D60I+u/OEpZcTlJ10QTO8A8Mh7RO4iI3u3nvebh5QaVv748fa0DD4VCwXWPJXJqbxZabw1dOqsxvvgiq/77LvsVoIoYjmb0Y9X19Zv/hzHnFO4TX0ShaXiSuCPBrzMB84ldFwDYuzaNe/47usG22wIJpJrBJZtHqrkolBCUADm7qwqg++3mhz7NOCKlzzP/UWKAFtgKcdv6FNz/s4lb7xqMh7bhZHbJO9K4Z9ZXGI0mPLRufL5kLlGxbeMb8WXv3HrIO+JUEFV0Dja8DWU6GHoDdEqCH/YM5p+/TSe2YBsDq5eOujiIUrmByUSJwhNPQ828j3y3DmwOnYV3+UVCeo4iMd68IXhOQTHhRhMd/M0TyRNibGwN1cZ0ivFj6iMJFslwwTzHqCi3jPTDuRzdfh6tjxvjbu+Fh5cbXv7uXPf4gFbqcctSa1TED+pY/Vz5yitMePRRJowfT8XB1WwOjuVU5wQK0ndhOLQUjOXovrsDVYd4Ep5+Fb8uYaTvSyP1mK6eq9jXmIDZoDdW75ManRDCmFnd8fBuW0lAq0gg1Qwu6TxSzWXIq7D6JnNwk/hMzfYxHQaDtqNl2gNX0eeByciIETFc2Jvv0ClGT/Dz8KDwomPbzdTm6e7G+//dTGZGAc/+8+oG6//776sxGs1zIcpKy1ny7T7ZhqYlGHSw559QeBp63gcdkiyP73sdtj3qdLMrXoOUyru7Zz+Ege4xvLJiCv7l5xiYvaTp/bZBPeh23BJnYSrNY/2h8+RkXaBn0Wao0HHEbyQBEZ157qlJbFtSs2uAI/MAXaUpt4V8Q7S4aZT0vKIzvUaEWQVRYJ5j5K71JijMm35XyfxEC6GhsHs36p07GXXTTYw6+Dlb8/LYXHnYVHSBwNgwRt02BID+Y7vyWUKCOdeY2gPN2KdRRQxtka6e2p3Fqd3mFaBTH04gvFs9OwS0AoXJlXuyCAtVgZQju0cLwGgwr3JS15mLdPEgfN+7/nM9O0H4RDj2sXPXnJfPhsUZ7F971qJ42c7DTB5kvf2Nwl/J3S+NJHlFKjt+qSerug2H0zLZduQMQHVyuPoMjPmXVZkj54lGyNwK6+8wb55trIDSOvkwtKHQa745KWwDW7iUFsGyN8xJMmNCYeLtsKM8nq2vWGbj1ivc0Zga9w3fEQrvDnjM+sqi7GJhCYE+npzPLWDQqAjS9rXM6i9n3fCXAYRGmb+EFufpWPPZYc4cqtl4vEOEDzP+Wv+m4sJ5pTk5/Dh5MhnbtwMObGbt4Y8qtBszvnsXz/BY9i7bx8G1p1H4dLRot75be42ZJ3fjUwPxDfG0ubmyqzKqO/P5LSNSou1Qqsw/dQX2MueUytxsfaw2j/pTE9iky6X/6GAOrD2DCfO1/ziejlppe+WIyl2JSqVk0OQoNv50HHelY2+h1Mzc6iCqSnGRjr//dTk7t6Qy+Ioonv7HRDy9zHMcbH2/CQp2fKTAZDLx41d7Wf3rEXr06cjdD10hyUDrs+k+yDts/3hpJvzxrN3DxWdgz9eg9YIstKSklwJwNAOOvgTLO09iEpaBVHMGUTnu4Wh6z6Tu8ohAH/Ocl44Bvm02iAIIDKvJT+Xl7841D/anXGdg39o0jAYTfUaH13O2aCxtUBA3b9qELj+fX6ZO5czmzfVvHVSWR1CAjvBB5pXKY+4cStrC+8jO1eI+8cXqavUFR4257fftP/6we6w1MqrLb1bRPtQdpbKiqJmg7gx9Hr7qfKbHP8yRnPEEalOhqJjVxdNtd6P20lwnVunmFBZbPF/wjzUs+rBmW4Xflh4mcXAXbrglAYDsC0X0juiIu5uagpIy8kvKCAxxPJDau+ssLz+zEoCdW1MJDPJk9p1JDZx1mTLozdu1NJKxAr54GQqqpzaVWtWZdPY/jW7fIe6+KNw8MBVd4GzMzRSEjyEuom0GG9c9nkhwuA+5GcXsWpmKWqNk0JQovHxr0gmoNdZvLjd3FQMmRrZgTy9PSrUabVAQU5YsYdnMmaT+/nu9WwfZ2vDaWFFB2dfz8Jj5SYv1uzVJICXaB4UD3zB8Ip1vV5cLulw6eR+kk7d5n68jaUPtjkipNSrzUiuFArWbChzIt1hUqmPv6QyLstpBVJV/PvNbdSBVkF9GTKeg6hEEgNRCx0cQ3nxlncXzBS+vlUDKnjPLmnT6zp8CKTBebLhic1BpcBv9BOrYMQDoyiuIdXP+17or0xDEDQrFP9STPavO0Gd0Z3SlBg5uMN86HzApgk4x/gB0iPRl0t1Ny/Iumo9ncDAzVq8m+9Ah8rZupfTOO1lpo569Da9N+WmUfjjZou5tR45QrvDjx38nY6r84lHfPDlHbwnqyyr45C/mOxbbfjrJ4GujrW7vNecmyg6943Jzc1m2bBkjR46ka9f6l6OfO3eOb775hocfftglHRQCaDiQUijAO9L2Ma9w+9nRL+6DCssRhGDvQtzsJIULKFoLH46CzmPx8XuSkmzrLUBK9eVoNTWrS7y17sSGBVNRYcBgMqEvryAzr/6cUrqyCnw9LRP+5RaZJ7jv3ZXOuws2oVAquO3eoSQmdbGaK5B6upU+2NuD9FWwdg6UF8GQ12DTPU1q7ovjw+jKUhd1zkFqd0CB+9T/oAyOrS52b0QQBU1PQxDcxRuNh5reozoTm9gBhVJB0pSo6uOjb+nW6LZF6wru2ZPgnj3hjjsIfOstzs6fT90UyXY3vK6wXJTzcWwkc3fvZs7fh/HZ/20xlzXytl/tY7UDsH1r0tm3xvr3/a0vD8Pd03rVX1WAZTKZqNDXrJgt1xvQNHgnpLINRyqlpqYyZ84cAgMD+eyzz5g8ebLduhcuXODTTz+VQEq4VkOBVOR02yNSYWPAszOcWGT7vC1/NqdfqCXYqwA/L63N6tEeW8wT4tOWE+Q2gxIirOoUlRYSEhnG5k2nKCrTEeTjxYheNR8qBSVl/LB5f70vpzC3DLXK8jVn5xVRrjfwyF0/kJ9n/gW1fVMKA4d05b8fzbDYakGjaRu7rrdJ6+bVTCZ3IojamRrDi8uvo6zCjUeuXMqEnvswmhScLg2hSdnOVG7gpoWyhrekUHbqh2bC31C4+zTlii7V78ouXHFjXGt3Q7SAzvfdR+f77qPnu++y/p57yASc/cr2WUICM7fsaI7u1evTv26xWV41p6pCb+S9P6+3OPb6HY6NVjv89cXPz4+ysjKuueYannjiCV566SVUKvllLVqIqp5vBip36Pe47cSeuovQcWT9bZss8/a4ufnRO6KjzaqB2tTqxxGqpaQx3/qSumJufTyWwddF8+jdP1KeZ6R7l5ptcAzGhvMElRZa3zNMS8/j8w92VAdRVf7Ydob1q48zbnLNKkM3N3lvWjGZYPsT5q1bGnHqyyuncTbfnKjxHyuvY1TcIYoOaVGaGrefHoAiMArNxJdQ+oRiSN2K/rfnzZOuqqg0qPvdiMEjEMXZP3AbMKdZgqimpCEYMMn6y4S4tHnffTeTJ06E+HgO6/U4e3P862FJoPbgtkOH8K1zl6v2Lbz6/l9aHNPXPK59y6+kUM+ip+vfvqiqnabsx+hwIBUZGcmiRYuYMWMG//rXv9i6dStfffUVYWFhDZ98mZGEnM2g32NwutYeY+EToOc9kLMXom6wnwFdl+v0NjN+nrYDnSDtCTp5Hah+HuW3hU3p1oFUTPAZWHEN3ROeYs5dSfz2meVu9wajicjQAHp0CcWj8lbMyfM57DudAYuToO9j6HWW6R4qDOY+vfWa7X0H339zi0UgFdMhiCuiIjEYjBiMJtKy8xp+4ZeaijJz5nF7o5ENKCjT8tRPN7MzNYbYoPOcya3ZyqhY78EVr5tXJfUvq/+2nlKtJrhPH7L377da9aROnI3SJxQAVcRQPOZ+j6kgA8Px3zGcTUYz+gmUIXG4AfSZ2qjX4YimzInS+rSPbTyEi0VEwPnzhAU2YrU0QEUZH8dHc+Wbb5Iwf77N+Uv1/b+0d6z2LT83Xc0XSnsBlq12atct1RVbHa/LqRvqvXv3ZufOndx111189dVXJCQk8OWXX3LVVVc508wlTxJyNoOQJIi/FY59as7pk/gsdBwGkdMs63UeC2dX1zxPfK4muaeDvDXWg9VxAb8zqut/qP1e93U/T1RMMadPWq6m6x56GrJ2wG/TKE35lQFxlqungn29GNM31qJsQGw45RUGnv04ik0nD+LhnsXkpJrASKEwz3/RlVfQLTyEpPiuVBgMbDhwmrM5+ejKyjm0L4PiIj0DhnSlQmfAM7DmAy6nsISKCqPlqsNLRfE58yhTUII5fYbJBAcXwpYHmtTsL/sHsC3FvKfc0Wzr7N79c5bSPb/+AMRuHh7PYNTxY1HHjLaor3D3QRHigzIknraZw9nS2Nt6tnYXRGsKCMBPryesSxfOZWbWW9XeF4o1DzzAmgce4JHycpTq5lv/Zi/AaqhuqQNZSpzutbe3N19++SUjRozg4YcfZvz48Tz77LM888wzFptaCuFSCgWM/gSGvgEqrf10CINeggs7zIkV/btD1HSryeQN8dWcx0eTQaG+EwA+mvOMjfwnSoX1SNXYkYf4RT+Y82nlBGm2E+SVT7hvzQiUR/khSkhw6LrlFQZ+PZgIgJfW8tuZSqnkltEJfLE2maT4rqhVStQqJYO7deHHLfmcTctn7vTPzH8FQyNQ6C37ajAaKSrU4R9ge+5Xu5SzD5aOMd++BXOusWvWQ+4Bh4Iokwne2jiej7eas8U/dtXPBHkV8un2UfhrS6qDqJoTjPTN/Y3Y/K0N54Cq/P8Z1Lcv/vE9KdcZ8O/Wi6C+iWQdOorHDe+i0Po7+4pdLqJ3EJF9gvD006ArqeCP5SkoFAqGTo+hXGfEw0tN8spUMk6YM//HDexAaVE550/mE9k3mOgEGxuOi8uLmxvXHTlC8pw58OuvJBqNbAD21arSYGJP4HU3N675/nvir7uOu/4zyualqiaE52eV8OOryTbr2LvlZ0/V6FNTVqs2Ovy79957GTx4MDNmzOD5559n8+bNfPHFF41tTgjHuDewNUCHJLjpGBSegsC+4OYFGn/zqFTBSXMdNx/zvKqybJtNKBQmxke9yJb0ezChYFjnd20GUQCavU9yfQhg5/PEQ38EHAykfL1qgkNb210A3DiyH+pax/y8tKiUCgzGmgSeO7emMrSH5bwVg8F+ILV/91nSUvMYPjoaP/92EmgVn4Uf+lmWZW7GcPBd1n/yFgpFL0bFHUKpsL1xQ16pJ+uO9eDzzcNRKAyYFCr+/fu11cfdjGV4GvMoUfkBJmaefsqhbql6TMFt0DwUWvP/0yKwnMA66GW0rZiQO3FCV4ZOj7V7vMcw66kakX2CyT1fjKHCRHC4t42zxOXOw9+fYb/8Uv18fEoKAVFRVP3PD+7dm5A+5lQXIX37Ety7t81Vfr/ccAM+XbsyZ9cuPINt7yuq8TBvznztg/05vDWD4zstR8LsBUP2AixbKwOdnS/VpHG0xMREkpOTmTdvHkuWLCEhIYEnn3yyKU0K0XSeoeafKgoFjHwfNtwFhlLzkvc9/7QbSAF09DrCdd0eanJX3NWl4GAC635RYRw/m010pyA6+NlOvulmY4GHu5uaEp3lhGdVndFhg9FEbk4xv3y/j2NHspg8rRdjr+7O3/7yK798b15BGNrJh+9/uxOtZxuf85J7CL7rZfPQc09tYMWhOdXPJ/VK5qpuBxgdZx4lNJng9TVT+GbnEIZlfskNJc9R4BbMxtBbKdSYo+EOpScZnrkId6N5JPOkT0ORjwJVzGgUfuG4DZrX5JfXXGb/fSh+IY0LlAM6ttz+e+ISEBnJoORkPBMTWY7tpJ211b7tV3jmDJ/07s24d94hdupUu7mfuvQMpEvPQMbf0Qt9WQXvP2R7/mgVZ0abated88rABus3+Yakn58fixcv5vXXX+fJJ5/kgQceoF+/fg2fKERLChsDM4/XPD/8XotcVuNmnWeqPjdc0dfpa0SGBnLojOW3MlWdvFIGo5G339jIjs3mVYcbfz/BI5mF1UEUQGZGIcuWHOSGWxLYv/ssRw9f4IrRMXQMc26fSJPJRFlpOW5qUKvtbPtjz6nvzfPgAvtC7wcg9RfzYoHOV5lv2W68y7zAwIaLxV6sOGQ5+rf8YCLLDyby6dz/0atTOnvPRvDVH1fQueQg4SXm4Mq3PJtu+Zv4I8SczT4hZ2l1EAUQU2idPLU2t2H3ou5zvc1jrkx02RTxSaGNDqKEaJSEBHqtWsWhceNItZO0E2zf9ivJzOSn6eb3473nz+MVGmrvKoB5VKnqdmBFuYFzJ/LISS9i59KUBrvZYqv2PD096dXL9jfAKo888ghDhw7l5ptvRt2ME8aEcIlmynBbl4d786/cHNytK8fSs6iolVahdkJQgH7RYXy9frdF2esvrrFqa8PqEyhQVG8v4+Wt4aUF1/LQneYVk2PGx/P8q1fj5e1udS6Yh8yffOAn1q8+TmRQFq/f9B1dp78O0bYDDQsXD8DqGebHZ5bCnn/UepGvwvbH6z39fIG/3WP/WXs1793yHssPmgOtQVk/WhyPLdxeHUgF6DOszq+PvSAKmp7osimGTo/BzV2Fyk1J9yG203kI0azGjmXy55+zbM4cUu0k7Wzott/bHTsy5Ztv6H7jjXYvo1AoaiaSu6uI6d+BmP4dGDApkm2LT7JndZrdc+v7cvPIh5OZ80r9k+jBwUAqPj7eoflPQ4cOJSUlxZEmhWhdRhsBTmAfuFh/okxnaVWFLm3PnptG9uOLdbvrreOmVlGqtxwhUykV9I0Kw8/Tg6PpF9iy/hRb1p+qPl5cpK8OogDW/naMpOERTJ/ZH0wm8zY5tWxce4L1q80jfyk5IXyxJZG/dnzYPOlf0cBilB1/rX54MqsDeoMbcSEZ6CrUeNkJorIKfcgu9qVb6DnKDfYDloMZ5pWTP+wZgkdFAR5G6yXN15z5J+6GOuVKNXj4QVmejf8zChRetudxtCaFUkHvUZ3pOybcvKWREK3Ic/ZsZsyeTXaXLmSlp3MCOFrreEO3/QB+u/POegMpe1QqJcNviGP4DeaEsRV6A2eP5VFSoGPNZ0caPF9fJzO7PU0eOiooKOC///0vhw4domvXrjz44IOSW0q0AzYmjw9dAL/fVO/cKWdpla5rqz5qdcMfmBq1mroTtvpHd6ZvlHl1YmRoAN9u3Gs136qufz23irdf30hxkY77Hx/FnD8Nrj72n5fXWtT9Yc8Q/jphCegLLBOmlhdB6s/m7XuMFbDxbig4AcDHW0ezcMNEi3Ym9NjD36/5hkMZ4Xy3ewhuSgNRwRd4a8N4dBUaErucYt6QdTb7qzAZ6Jh3iOuemYnCy0BswTab9bwq8iwLPANxv/pfKIOiMWafQPfrk1CaC9oAFGoPPG5pOD9VUxJd1mfETXFs/Oa4VfmcF4fiGyy38ETbE5yWRvCiRcTNn49XQQFVa+7s7dVXm76wkKKMDLw7dWpSH9QaFRG9gzCZTMQkdODItgxSD+Sg9Xbj6PaGR57sUZhMJtvLWur44IMPeOCBB/jiiy+47rrrANDr9QwYMIADBw6gUqkwGAyEh4ezb98+AgIaWF11GajKI5Wfn4+vr3PzTEQz+2U0ZFhuB8BdJvNs5L3/shgdaYrTud359fRCl7QFsOVwCsN6RNo89tnvuzAYjXQN8eeq/tZbdqzYdZSMi5bbkNw2znIi9cHU8+w4Zn8Y3JbV/z2Fvz4ZlO5MeOEacvIsj79+/aeMfPI7OPQ2lBdCvyfI+G4eby2L5mhmGCqFiYIyLdP77WBI1HHmfW6d5BQgwLOI3BLHV435687ibighrmBr9XyodM+e1Y8boh4wB7eB8xy+Xkv50xsj0WjVZKUVoisup/BiGTlni4kbFEpopPyeEe1A5dSKs8BXTp73UGkpanfbUwuaoqyonM0/HKcoV0fCuK74d3Fz+PPb4RGpN998k7/97W/VQRTAp59+yuHDh/nxxx+59tpr2b9/P1OnTuXtt9/mqaccWy58KZLM5u1A/79aBlIRlUvfFQpwD3LZZTzd6s+KO6jTJ+zMmOdQW2cu5HLyXA5qlZKkeOtM7reM7o9CobBasVfFy73hFI++no5t0mlx3aeC+HX+HwCY9KMByy1MHvnhVqLXfchnt76Hh1s5HP+M5766i+Q0y4zz72wazzubxtu9jq0gSmEy4G4oRq/yxKio+XUWm7+VgTk/WdWvG0Sp+92Iuvd1GAvOUb7uFUyGctzHP48ytG0mmrzxqUFotObXGdKl7ey3J4RTysrgH/+g8wsv8Bjwb0fPM5lY4OHBXampVlvLNJWHtxtX3Vrzvi8oaHjvyyoOjUgVFhbSpUsXLly4gEZTszR62rRpaLVavvqqJqZ87733+Omnn1i2zNnddy49MiLVhhkNsGIypK8Ej2CY8DOEDjUfO/mt+RafCxTpA/n0wHd2j8/pNYvPDzY8/7D2aJJSoWBgXDi97OwHWJ+lOw4xqk8MPlp3SnR6PN0tUx2kZF5k7b6TRHUMJMjHk5TMXHKLSugY4EtxmZ68YtvJTe8buYJbBm7m2nee4GKJ7Q94lcLA1b12E+hVxKfbRzvd97rcDUWMyviIQP058jQdWd/xNuLytxDnVYK7oQRjzkkUWn9MxVnW85uUKtxGPYY63jJwqzi0FHXPKU3umzP8Qz3JyyyxeWzw1Gjik0JRKBR4+WlQ2skvJkS7pVBQArzl5GnRkycz/eefUTRTInBnPr8dGpE6d+4coaGhFkEUwM6dO3nppZcsynr37s0HH3zgZJeFaGFKFUxaDkWp4B4ImlpvFI3rtvXRqOxPVjSZTGjdch1qp0d4BwJ9tHQP74DW3Y3jZxs392pKUs03rrpBFECXEH9mj0nErXLOVZ/ITpRXGHBTqzCZTGw4cJqiMh0jekXhplbxx/F0TpzL5q0NEzlyvjMKFAzvGUl0xyAuFpawbt9JinV6AAwmFb8caDgni6N6564mUG/egNhff56rzr1LwOBZqHtPs6hnzDyM7te/gL5mdFDVbaJVEAW0eBA19LoYEsdHYDAYObU7i98+OGhxfMDECLt5dIS4JCxZgue0aTwI/NfGYXtby5xatoxNTz/NiH/8w8ZZLcuhUC4sLIzU1FRyc2t+6R86dIhz586RlJRkUTc7O5uQENk2QLQDCgX4RFoGUeDSQEqltM7GufHAKZZsPcAvOw7hZuO4LRGhASTFd8XX0wM3lYqeXevPq9JYKqWyOoiqUvVcoVAwKD6cId3M/dBq3BjaPQK3ylGSNcf6EOIfRHznENQqJR38vekd2TzL7tVGHXF1Jo17K/VWQRSAMrQHqtgrawrcfVB3v7pZ+uWsxPHmDPQqlZK4gaH0GVW5p58CrpgRJ0GUuPRNnQoXL6JZvJjHgJG1DlXlmJqbnMycXbus9uLb/vLL7FqwoCV7a5NDgZSPjw/9+/dn3rx5HDlyhOPHj/Pggw8SExNDz56Wcwk++eSTBnNOCdGm1Q2smkClMGAwWq4Q7BsVRkJMZzSVWcp7BP1afcxbk0l6Rv0ZeluTp7uGIN+aLNdqlZIuITULS3pGWm4Q7eqAT6lQ0DuiI6MjPVB4WX5hU/jaXy3sNsw8gV3ZdQges75C2aG7S/tV29DrYpj+aAKDr43m2gf7o3Kz/Wt2zktDrcpG3tyNmc8mMfuFIfS7qkuz9VGINiUgAKZNA72eJOCKymJbOabqWvvww3zUvfnez45weLL5a6+9xrhx4/j555/NJ6rV/PijZVK77777jsWLF7N7d/35bIRo01R2lo8PXQB+8aDWwsGFcO536DIJTnxZb3MVBqPFBHA/Lw/8vDxIy8oDYFSX/xCsPUmZwYfewb+QGuzFfV9DoU5LdEwAH0x5iI/2foGbm+tGylzJy8MNdeXrc1NbT2gP8fMir6iUcoNlQOnnZR5dM5lMlBsMFJTUjM5V7R/YKcAHPy8PUi/kUaovZ2iPCOI7hwBdMIX/j7Kv54KhHExGFCr7k+mrjrlPesluHVcYMDGiepQpLM4cYF7/xAAObjyHrricnLNFlBWXM/yGOHyDbP8/CwqT/ezEZcrNDUwmYufMYdOiRQ7lmAK4ePQo/1YoeMyxJAQu53D6A4CUlBQWL16M0Whk/Pjx9KmMFAEqKipYuHAhKpWK+++/v1k6297IZPN2ymSCb7tBfmWeHt9YuOmo/YSS79W5/eIbA+XFUHoegOeWz2fzyXjGJ8bj5VEzL2nr4VQWzbrdZpN5JZ6cyw8gevi1eJx6m4U7vwBV+81OrSuvYOmOwxSUmOeM9Y7oyKD4mhEXk8nEtiOpnDp/kav6x9IxwPL9UqYv56dtB7lpZH+rtk1FWeg3vI4yMAq3IXfZ7YN+/WtoRj3qmhdkQ5/R4Qy7PsYqSakQwnnfRUSQeuaMeY5UPTmmaus0ZAiztm51yfWd+fx2KpACyMvLQ6vV4l5PHoc//viDgQNdN6m0vZJAqh3L2Qvb/wKYIOmfEJxgv27dQCqwD0xZC2nLYe0c0nMD+duvMwjpMIHOQTWjSjuOnuHTm29zqDvfHv4vWaWNu2U+Le4hlhxf0KhzXen4uWw2HTwNWOevas/cPFRcdWsPYhI6tHZXhLhklGRns6x3b1Iz7SfKtDcR3RUjU858fju9bjAlJYU5c+ZgNNrIDA18++233Hvvvc42K0TbEtQPrl4BV6+sP4iyxVgOHkEQNxvCxhAecJH3Z72LwWj55lY6sWzX061xW83oykvxUDueD8VZt71yBXf9ZxR3/WcUt71yRb1148LM26mM7B1db732os+YcO58fQR3LRhFdH9ZYCOEK3kGBzPj3DnmzZvHNGBSneP1TUT/t0LByaVLW6yvTgdSnTp1Ys2aNfz5z3+2Ovb9998za9Ysxo+3n1RPiEuOr+UEa7r/qebxoJo5OZo6q+F8tNbpB+xxU9nO39SQ8vIKPFRFjTrXEVWb8jq6Oe9t4wYR08l1CU9bgq14N6JPECNvisfd0zz3SlbXCdEMlEqCP/6Y2EWL6AX8Gaj67dHQRPTF11zDN1deycWjR2luTgdSoaGhLFu2jI8++ogXX3yxuvyHH37g5ptv5v7777fKLSXEJW34QlBWBkU+UdCt1ryn0KEw6mPoNJqOAZaJKruFO34rSKvOa1TX1MoSNOrmC6QudQoFDL0ulvjBlqsPk6ZEtVKPhLgMzZoFFRW4AddUFlVNRAfsTkRPW7uWj7p3Z+3DDzc4v6opGrVp8eDBg/nuu++YNm0aHTt2JDAwkJkzZ3LffffxxhtvuLqPQrRtXSbADfuh4Dh0HGGdPqHbPPPPsjXW53a/E440nMA20m8b+7OmO921Tr6ZqBU6Aj1Oc7HM9R/+zbUpb1tw68vD8fBSo9aoyMssITejhLwLJfQZFU5IV9meRYgWpVKByUTw558TMXcuqQ5sdlxl14IF7FqwgPuysvAMDnZ515yebF7bp59+yp133olCoeCuu+7if//7X731Dx48yC+//EJAQACJiYkkJiaiUlnfDsjLy8PPz6/dDpfX3mvv2LFjMtlcALDwHutAav6zKlg6usFzjSYlS46/RkZR3+qyKL/NXBX5Lz7Y+7Pd82L91zIh+kUyinqy7szD6Ay+FJe7/hfJpeiehaNRyZYsQrQ5JXPmsGzRIlLtHLc3Cd2ZVX0uX7WXn5/PqlWrbB576623uHDhAs8//3x1mb+/P2PHjrWqO2XKFCoqKigrK+OPP/5Aq9Uya9Ys7r33Xrp161Zd7/fff+ett97i+++/b7fBFMiqPWHJZiD1dz9Y7NgK19IKX45dvAoFJroHrURTOW9qYfLvds/pGbyUMV0tR4nrq29LgHsquboIp86pMBhRt6MgJCjcm07RfhzYcBaApGuiGDRZbt8J0SYVFYGPD3uBupFJ1ST0qtxTnw8YYDVa9aeUFPwi6v+d5vK99k6fPs2MGTPqrVP7eL9+/dizZ49VnSuvvJJ3332Xhx9+mI8//pjDhw/z3XffkZCQwNSpU3nhhReIi4vjqquuYvHixaxevZpx48Y50kUh2jytjxulheWWhW6OJ1/Uqgvo12GxVXmE7zZSC4bYPMddWWxVNin6WZafeqHeaw0Ne5+zRf04U5DkdBAFtHoQ1XdMOPvWpjtUt/eozgy/PhaVm5IewzuhUisJ6ixJMYVos7y9wWik3z/+QaennyYbqNofwtYk9At14pH3IyMBmL50KTGTJ9u8RO6RI/S0ecSaQyNSer2ec+fOOdgkaDQawsKst2swmUx8+eWXvPPOO2zdupXY2FiSkpIICgri119/JSUlhT//+c/07t2bF198kc8++4whQ2x/QLQHMiIlaju9L5tf39pX/XzkzHj6DFLCF52db8wvHvKPAZBXFs66Mw9ToophgP//6Ba0GpNJgd7oiQITGlWJ1en1jUp18trPtPhHKNB15LfTz5CvC0NvbJ7AwpiXRsXOj1H1no6qU5+GT3DQ4Guj2P7z6XrrdIj0ZcaTku9OiHbNZAKlknPAlzg2IlXbI+XlFqkTSrKzWXbzzaSuXg3AXc2RkNNVcnJyWLt2LVu2bGHPnj2kp6dz/vx53N3dyc7OZvbs2Xz++eet0TWXkUBK1GY0mkhemUrKvmw6xvgxZGo0alMxfOLk1i8BvWDiMsjeBauuN5d5hsHM41B8DrJ2QKdRoMuD7633pgLrQCrQ4zSTop9Hb9QSrD2JUmHOE1f12+H9vT9TbvTC1Uo/vhb05lEz9+vfRRkc65J2YxJDOJmcZVU+++9D2bnsNEqlgqRrovEOsJ9YWAjRTqSkQFQU3wGp4FQ2dIDwUaO4+rPP8O3ale/GjSN19erqeVbT1q1ru4FUfZYvX87dd9/Nli1bCA8Pb+3uNJoEUqJBRgN84ODi2fE/gX938O4Kag9zWXYyXDwIkVOtVwvq8uDTAKtm6DyWhb/81aIo0m8zk2Oeta47ZhFETefdh7dTUW75q2JS9LOsT3uQEjuT143ZxzEVZaGKHGbzuElfTNnH11Y/d5/2JspQRwfTwXgxBWVgpM1jQ6ZFs23JKYuykTPj6TO6/f4+EULUY+pUSn7+mWXg9CR0W/WqRrUKCgpcn9m8JUyaNIkDBw4QGBjY2l0RonkpbSSxvOJt23X94sA/viaIAghOhPg51kEUgLu/dVlQf5i8isg+lkkxB3S0s/GyeyCoPbG1kUG0/2bm9b6JEI9tNk9VBsdRcXQFJr31rUUAhcZyhMuZICr+1HvEswKvldbB34ib4uk1ojN+ITWbAod3D6DnCOvpBkKIS8TixXi++SYzgHlA3d+I9WVCr6v2PCtHNCqPlD3bt29n1apVpKWloVQqCQ0NpV+/flx11VVOj8jICI64bHS7HY5+ZH7sEwU97gLvLrBiimU9LxcEAm7m/EcDro7kwplCSvL1dA9aTqjnEeu6Gj8IGwVAaKQvGSfzLY+P+RyFRwgdVJ5k/VFufT5gTNlM2ddz0c793uZxhW8YpgLz/Etj0QWU3pZJSsu3vI3h9EY8ZlkGen3eeZeO0X5w4gTnrrmd3/vfRYFfGIMGedJ3jHnU6cb/G0RmSgF+wVp8g7UIIS5hSiXcfz/cey/BSUnclZxMPnACWEv9k9CrRqqy9u3DZDBUJ/t0NJhySSCVk5PDzTffbDdFglqtZubMmTz99NMWaQ6EEJhHoAL7gi4Xet4LCiWEXQXaDlB6wVwnKAHcGvHlovtdcOS9muf9HgegY5Qfc18aRkXKGtzX/NvynLArzf3pNR/UnuSsWYPmtw8g5s6aLk9QmvcSBFS+xwDrFXLle74xPzDobXbNmJuKqSCj+nnF1nfQjDOPMJnKy9B9Mw9TcRaorecxVY82xcYStu475uzeDX27Q61FLhoPNV26y6i2EJcVlQo2bYKICPyysqhac1w7OMrat686K3rdyemf9u+PsVayz2nr1zd4ySbPkSovL2fIkCEkJycDEBERQZ8+fVCpVGRnZ7N//34KCsybpmo0GhYsWHDZbGosc6REk2RsgJ3/ByotDH0dAm1PHK9XwWlYNR1y9kLsLTD6E1C61RzP3Ao/1ZnDNGk5dJkIgC4jg3fCwigHVD2mEJA4lW43jGHQ5EjUbubbkqse+pBjdbKmm8oKKPv6VtAVgFKF9k+/WXWt7Id7MGUftyhTdk5EGRSDIWUzAxLiKDlyhNzMTNSjHiS7+1QAeo0IY/Ss7s7/XQghLi9LlsD06fVOQu/Qvz9zd++uPuWTPn0stptxZNVek0ek3n33XZKTk/Hw8OCjjz7i5ptvtjhuMBhYv349L774ImvXruW+++6juLiYxx57rKmXFuLS1mkkXLuxaW34RsF15i85KGxMiVRbr8TTZ+SwYUwsF9LTuajTUXXTznB4KdmHlzJvUa3vXmVlpG9aDwMtA6my7+40B1FgnlBfR9mXszAVnrcqN55Nxng2mfgrrmDEqlWo3GqCvvOn8jEZTXSMcXKVoxDi8jRtGkyYwOSVK82T0CsqrHJK1R2punikZprDx8BdDlymySNSAwYMIDk5mVdffbXe4MhkMvHiiy/y7LPPolKp2LBhA8OG2V7Nc6mQESnR5uWfgG/iLIq2fdKPTVv32j3lsdRUCA9n19VXs3HlSoxh/VGGD4AKPRh0GPPSMKZabsOg6n0dbsPuRaFQUnF0BeXrXrVq1we4c/FimDzZIoASQogmeestmD+fbCAP8K/8c0nlYXvpEh4H120RY09JSQk+PubJq9nZ2QQE2FhqXcfdd9/Ne++9x4QJE1ixYkVjL90uSCAl2rySDFhUM6+oogwW/Ln+Ux4FiqOieOd0/Qkv61L4dwG1R/XtvBtjYghasoTjr72GLieH3q+9hldcXAOtCCFEI6xYAZMmWRQdAOqLQlokkDpw4AB9+vQhLi6OY8eOOXROTk4OkZGRFBUVcfr0aSIrU7VfiiSQEm2evqA6IajJCF88BOd19Z8yBTgKVM1ucjQ3S13XRUcTffJkY3othBDO27kTkpIsikyYA6ozwOE61R0NpJqUR6pqEnlISIjD5wQFBTF1qnnS6MaNTZz/IYRomlpzpM783nAQBbAUyyDK0dwsdel0DlxMCCFcZdAg83YNy5dXFymAPsBk4B4gphHNNimQ0uvNy5rdnJzPMHz4cAD++OOPply+zVq4cCE9e/Zk0KBBrd0VISyZTKQ89xwH5sxBn5ZmTgjacQQAaXucb85WbhZHndRonL+gEEI01cSJ5oBqr+VcUG9gOvBY5Y+jWiWzedWGxhcuXGiNyze7+fPnc+jQIXbu3NnaXRHCwo7x4/n+hRdYsWgRX8TEYNTpYMznEH0TKJxfDVe14gUga98+i2XDtXUG6m4/nvTGG05fTwghXKZvX/jNOjWLs1ySkDM3N5fjx48TGxuLQqFosL6fn/kXdl5enisuL4Rw0IbKHc0BcsrL2TNmDEWFhegyMtibk1/PmbbVTlxX3waho4Dg6dPRZWeTdfw4vebOJeTaa23WFUKIFjNuHOTnw/z5sGhRo5pwSSC1b98+4uPj8ff3JzExkQEDBjBw4EAGDhxIdHS03fNKS0tdcXkhRCOt2bq14UoNMNrIzQIQEBjI8PnzCUtLwzc2Fh5+mKs8PZt8PSGEcClfX/j8c/MPQEkJxevWweTJDp3epECqS5cuTJ8+neTkZFJTU8nLy2PNmjWsWbOmuk5AQEB1YFX1pxDi0qZxd2fi4sV0HjmytbsihBDO8fTEcMUVDldvUiAVExPDjz/+CMDFixdJTk62+Dlx4gS5ubmsXr2a1bVuKWi1soGoEC3NUFLSLO128PXlQuUK3iq3nz6Nd6dOzXI9IYRoS1xyaw8gMDCQsWPHMnbs2OqywsJCdu/ebRFcHTlyRG7pCdFCsn7/nW3z5+Pm4cGg//zH5e2PHDOGTnPm8M3tt1eXDZ8yRYIoIcRlw2WBlC0+Pj6MHDmSkbWG90tLS9m7dy/JycnV6ROEEK5n0On4ceJECisngJ+dMMGl7V8zbx7x770HKhUjt2zh6OLFhHbvzoBGTtgUQoj2qMl77Qn7JLO5aG6G0lKSb7qJ/JQU+j33HCHXXw+Y97bc/eCDrPnf/5rluvFRUVx76lSztC2EEK3Nmc/vZh2REkI0r01XXsnObdsAOHDDDdxz+jTuXbvyQ48epDi4bZMjugYHcyY7u/r5oH/+02VtCyFEe9YqCTmFEK5RFUQBVAB777mHk6+95tIgCmDoa6/hoVIB0DMhgY4zZri0fSGEaK9kREqItk6nM29nsGABrFsHU6aYk8fZSH7rX1jIlpdecnkXusydy13XX095URFeoaEub18IIdorCaSEaINMBQWcfeYZ0n/4gYNnz6IExgHhACtXQkQEXHMNaswjUVV8vb3JzXc+Q3lto4D1tZ73iooCQOPlhcbLy+Y5QghxuZJASoi2xmRiWWQkR3JzLYpXArcDWcD5a6+ly7Zt1N2QRZ+dTXkTLx8IXA3sAnyBkV991cQWhRDi0iWBlBBtTN7ixVZBFEAusB9YDRgBhtTdBhh+Sk5u1DW9AQNQDrjNnUvMwIH0HDTI5jWEEELUkEBKiDYm/ddf7R5raJ9yRzKzdQAu1Cm7p/aTq6+Gm25yoCUhhBCyak+I1lZQAGvWwLlzLm86CRgARAMBwEBgDlB7uviwuicNGODyfgghxKVKRqSEaE0XLsDAgZCWZt6B/Pff0V+86LLmo6mcoF7HTOD4jTeiHTmSyOxseP5584EZMyA21mXXF0KIS50EUkK0gsylS8nfu5fInBw0aWnmwoICePllSi7UvfHWeMUAPXrA4cMwcSJERcGBA7jNnUvPO++sqThpEhQXw6hRLru2EEJcDiSQagYLFy5k4cKFGAyG1u6KaCtMJjh6FPz8OPTqq/z6xhvVh/wwT/b2AGJ//JHS7t1dcslOQLfJk+GXX8y5qDw87FdOSnLJNYUQ4nIje+01I9lrT1S75Rb46isMajVvVNRNWlBjAFDQqRPHMzKafMk4YOq2bTB4cJPbEkKIy4kzn98y2VyI5nbgAHz1FaXAonqCKAAtUFZS4pLLalUqGWkSQohmJrf2hGhm6a+8wi7guAN1tYBOp2uw3kTgZANteiQl2dxGRgghhOtIICVEM8vNzXUoiALzPKkr/f0pOn8eHXAAsHWTL/SKK+gF5Pj7kxkayvIPP7Sq4yl74gkhRLOTQEqIZlZWVuZw3fXAHXl5qCqf9wP2Aqvq1FM/9xyKsWMJBoIMBi4uWcL2nByLOp7hthIfCCGEcCWZIyVEM0s9dMjhugXAG2VlfAicryxzt1HPzd+/+rFCpWLE11+jrVPHKyLCuY4KIYRwmgRSQjSj43/5CymNyFieC6ytfGwraYHKx8fiualjR7R16np26eL0dYUQQjhHbu0J0UxMRiO//fvfjT4/A/PmxJ2vvRb3n3+magq6j0KBR2SkRV1FfDy3BwVBTg4GoMTLC88JExp9bSGEEI6RESkhmkl5VhalRmOjzzcC+YDbjz8ydsgQtIAXMPa221C417nhp9HAG2+ApycqjQaf119HVev2nxBCiOYhI1JCuIihvJyMFSvoMGoUGl9f9GfPNqk9NXARCFCp6LF1Kz1OngSVCuqMRlWbM8e8V55CAXUDLSGEEM1CAikhXECfk8On4eHkl5WhVigYefPNHFqypNHtqYBeQNCsWTWFMTENn1jfNjBCCCFcTgIpIVxg7z33kF+Z5qDCZGLNl1/arOcPzAPeBmqn3ew3aBDB/fvjFhDAildewQCcVKsZ/MQTzdpvIYQQTSOBlBAukPzzz/Uevx0oBsI+/BDuuIO6ucsTnn6a4GuvBSB4zBiyNmwgYu5cfF20gbEQQojmIYGUEE4w6vUUbNuGz8CBqDw9q8t15eX1npcPRMXHw+23oz98mLh//5sSoLTyR1vrtl3HiRPpOHFis/RfCCGEa0kgJYSDSg4f5qt+/cgtLyfI3Z0Z27bh3b8/+pwc9CZTvefmAFGVO4hr/P2ZWrdCt27N0mchhBDNS9IfCGGLyQTr1pl/ioshK4tts2aRWznylKPTceDxxwEo2b+/weYKACoDKQIDrSuo5TuNEEK0RxJICWHL/PkwZoz5x9sbY2goybt3W1Q5vH49APrMzAabywfw8zM/uflmc96nKtdf76JOCyGEaGkSSAlRV24uvP02xYChsuiCjVt3FUYjpoICTj/9dL3NRVK5X17ViJS/P/zyC4weDTNnwn/+46qeCyGEaGFyP0GIulJTAdgOnAUmYV5xV5eHRsOmxES2nzxZb3MGoANwQa+nQ1Xh+PHmHyGEEO2ajEgJUUfGH3/wGZAMZAKfAOcxJ8isLbO0lF0NBFEAacA6ICskxKX9FEII0fokkBKXL70e/vY3uPJKeO01MJhv5O3973+5UKfqVszZxuuqqPNcBVwLTAf61yrrN3w43f75T5d1XQghRNugMJkaWLctGq2goAA/Pz/y8/PxrZofI9qODz+EO++0KDoyYgRLN25sdJNaYH7tgmHDYPPmRrcnhBCi5Tnz+S1zpMTlq04QpYcmBVFgTq5pnDAB5YwZUFoKt93WpPaEEEK0bRJICVGp4dlOjinT6fC84w4XtSaEEKItk0BKXLYygFxgI+CFeUK5K5QUFODZcDUhhBCXAAmkxGXJWFbGl0DVBMFCJ89XAZ52zispKmpK14QQQrQjsmqvGSxcuJCePXsyaNCg1u6KsKNozx5srbIIc/D8UcBdwBgbx7IKChrdLyGEEO2LBFLNYP78+Rw6dIidO3e2dleEHUV799os9wAGOnC+2623ojhwgAEGg9WxjDBHwzEhhBDtnQRS4rJUcPiwzXItMPSmmxp8Y1QUF0OvXqBUMmzWLItjg2XLFyGEuGxIICUuOfWmRjt3Dj78kMJ162weVnXrhvvXX3P/vn141XMN7x49qh8Peu89ku68k6jERKZ88AEhV1zRuI4LIYRod2SyubhkGAsL+S0hgROnTtE1LIyr9+1DHRhYUyEz0zyKlJeH3engCgUAmj59uPWLLzj56qsEdO/Onm3bOJKSAoBWqSTqL3+pPsXN05OR77/fPC9KCCFEmyaBlLhkHHrwQQ5U7n137OxZ0kJC8NFq8fTyYsxDD+H14Yeo8vJQArvstKGoDKQAPG+5hT633AJASE4OfrfcQmlODgMWLEDtVd94lRBCiMuFBFKiXTIZDJz6+9/J2raN2PvvJ3jKFPb98INFnVKjkdLiYiguZtVTTxEE2J5iXkMTGmqz3D0oiBErV7qm80IIIS4ZEkiJdkd34QJLevUiLTsbgE0rVxLu78+5QvvZoNIBnQNt93vtNdd0UgghxGVBJpuLdmf3vHnkZ2fTC+hfWZael9fgeVkNHB95zz34JyY2rXNCCCEuKzIiJdqdC/v3UwCUAZHAHgfOGQWsr+f4/BMn0MbENL1zQgghLisSSIl2p7gyc/hJHN9ouL4gKsjDA4/o6KZ2SwghxGVIAinR7hSXlrqknU6Ayt2dK7/91mK1nhBCCOEoCaREu1NcXt7kNjoBs5YuhcmTm94hIYQQly2ZbC7aFf3589QNoxr6T6yxURYQHS1BlBBCiCaTQEq0K8UHDliVBdqoV1uwt7dVmX9EhIt6JIQQ4nImt/ZEi6koK+PI3/6G2seHbk8+iULpRBxfUADffkvxmjVWh7IbODUkMJBzRZabwniHhzt+bSGEEMIOCaREi/kxMpIzmZkAnPn+e8YnJ1Oans62m2+m9OJFBi1YQMi4cdYnlpfD8OFw4AD2U27aF9SxI5w5Y1FWe9NhIYQQorEkkBItIvOnn6qDKIB9u3dzVVkZq664gmOpqQCcnjiRu/PyUPv4WJ68fDkcOEAKsMyBa/UFxmHOZF4CeD7xBPtvvZWs4mLAvOlwl/nzm/6ihBBCXPYkkBItIuPbb63Kyk6frg6iwLw3XsprrxE8diwb5syhuLCQ0J49GdipE9uAfQ5eqwxQAB6Ax8CBMH06V/v4sObWW6koL2fEv/6Fm6+vC16VEEKIy50EUqJFFJ87Z1V24aefrMoKDx5k67//TWbl6NHZjRs5AyQ4ca0ygEWLICAAxo4FpZKQ8eO5KSOjUX0XQggh7JFASrSI4lq39apkrFplVZa+Y0d1EFUlGwgBooFTDlxLBzBuHHTo0IieCiGEEI6T9AeiReSeP29VlmEjlcHROpPCq5wHrgF6Vj4PUKls1gsEfAGCghrTTSGEEMIpEkiJFpFbuT9ebRlZWQ6fnwG4AVe//TaP5ORwh43s5h2B24GpAHYCLSGEEMKV5NaeaHblFy9SZDBYlZeZTA63kQIweDDcfTdKO/viuVc9eP55J3sohBBCNI4EUqLZFW7dijuVc5fq0RWwfWMPSoEzQ4fStZ7NhTv6+8MXX8CkSY3qpxBCCOEsubUnml2gWs0DwJVAH8zznIbbqDe0nja0QMeHHrIo69e3b/VjFTBgxw64+mqoJ9gSQgghXElGpETzq1yxl1ir6LCNat/YKAvTaFAAI555Bk2d/fFGrVuHZupU8tPTSXjqKTzj4lzVYyGEEMIhEkiJ5lcn9cEfwLpaz/tiO9mmVqHglpISUCptjjJpAgIYtWGDCzsqhBBCOEdu7YnmVyuQMgBb6hweDETaOK2rp6d59Z3cqhNCCNFGyYiUaH61ckgVAvo6h09hngNV2yBgVPfuzdsvIYQQoolkREo0v1ojUvk2DvfHOpAqBfMWL0IIIUQbJiNSonmYTByZNYtjy5eTm5fHAMAT6y1eQjBvMKwNC4Na+/GVAvj7t1BnhRBCiMaRQEq43KmXX2bfggWcuHChumyFnbq+lX9qe/WSQEoIIUS7I4GUcKn9t9/Oyo8/dri+L8CcOWhDQqDWJsYSSAkhhGgPZI6UcKnkr792qr5f//7w6adow8MtyktA5kgJIYRo8ySQEi6VVVrqVH3fsWNBoUBbJ9mmDkirtdpPCCGEaIskkBJNdvipp/jIy4s3GpHvybMygPKMibE6diE9vcl9E0IIIZqTzJESTVJy4gS/vvwypkaeHzh2LAAesbFWx4zl5U3omRBCCNH8JJBqQGlpKfn55uxH/v7+eHh4tHKPWpc+J4e9f/oTJoOBfm+9RdbixQ4HUVqVCq1azUWdDoB+3bvjWZl0U+XlZVU/fMoUV3VbCCGEaBYSSDVg8eLFPPLII+Tn5/Pxxx8zc+bM1u5Sq/qlZ09OV6Y1OLluHVP+9S8mACtt1I0CpmKeOO6WnY3KwwM3Dw/SXn0VFAq6PP64Rf2RY8awYe1aALp4e9PxT39qzpcihBBCNJkEUg245ZZbuOWWW5g3b15rd6XVlRw7Vh1EAZwtKKBswwb6AFlAcp36/pj/g/k+/zwEBVWXd33ySZvtJ/3+O+GvvkpZRgYRzz2HQilT+IQQQrRt7T6Q0uv1rFy5Ep1Oxw033GCzzrlz59i2bRteXl6MHDkSrbZmQ5LCwkKKi4utzvHw8MBf8hhZyNu40ars/OrVhACxWAdSWoBOneDOOx27gEJB2BNPNK2TQgghRAtq14HUCy+8wHvvvYeHhwd5eXk2A6mPPvqIBx54gCFDhpCZmUl+fj6rVq2ie+XcnJdffpmPPvrI6rxrrrmG999/v9lfQ1tXkpKCQqVC26ULebt2WR2/kJUFQJiNcztdfz188IEk1hRCCHHJatf3ToKDg9m1axf333+/zeNpaWnce++9LFiwgN9//519+/bRp08f7rjjjuo6//jHPzh//rzVjwRRsGPaNN6OiuLtrl3Zc+ut5B05YnHcCxhd+VgNJNY65qNU0vXDDyWIEkIIcUlr14HUfffdR2hoqN3j33//PVqttnp+k1Kp5MEHH2TLli2cPn3aoWvo9XrOnz9fvXrvQq05QnXpdDoKCgosftqCrCVL+NLXl4/c3Tny6KMAGLOy0C1fjuniRZvnlKWmsvmnnzABRmDDZ5+Rc+KERZ14QFXr+ShgvLs7w4cPZ9bmzaj8/Jrj5QghhBBtRrsOpBqyf/9+unXrhpubW3VZnz59ADhw4IBDbfzxxx/079+f9evX89xzzzF06FC7dV9++WX8/Pyqf7p06dK0F+Aiv8+bx7nCQi7q9ax4/XUu/vADX4aF8ebVV/NVx46U7t1rdc75b77BUOu5HjiVlmZRxx8sUh+ogL6vvsrQTZvwHjLE9S9ECCGEaGPa9RyphuTn5xNQZ7+2wMDA6mOOGDZsGOcd3Krkr3/9K4888kj184KCgjYRTKXXeq0VwHczZlBoModA58rL2f/QQyStXQsVFfDdd2AykZdcd+o41E2P6Q98BhQC7oAHMPLcOSIQQgghLg+XdCCl1WrJzMy0KCssLKw+5mru7u64u7u7vN2mqKh8vbVVBVFVNqxbRxLA7NnwzTcA5DjQtj9QVusnHzA0w9+rEEII0VZd0rf2oqOjSU1NtSireh4dHd0aXWpx+owMxyoWFlYHUeBYIOUH1N2i2KNzZ0e7JoQQQrR7l3QgNXnyZNLT09m6dWt12ddff014eDj9+vVrxZ61HL0DtyXdgfzkZN4GPgK+AM40cI4PYMB8q7A2CaSEEEJcTtr1rb1169aRnp7Orl270Ol0LFq0CDDngPLz82Pw4MHMnTuXGTNm8NBDD3Hu3DnefPNNvv76a5SXSdZsRwKpG4Cy/fspBqxTk5p1BGq3VAi8aaOeR2Skkz0UQggh2q92HUjt2rWL3bt3AzB16lRWrFgBwOjRo/GrXHr/8ccf8+WXX7JhwwY8PT3ZtGkTgwcPbrU+tzR9PekaqvgAFyv/Hm1RYg6cHOEeIVPNhRBCXD7adSD1aGVOpPoolUpmz57N7NmzW6BHZgsXLmThwoUYDIaGK7uYyWRiw8SJHFm/ns6xscROmlRvfQXgCZw/dMhuHT9gAObbfWlYz4uqogFUMtlcCCHEZURhMtVZwiVcpqCgAD8/P/Lz8/H19W16gyYT7NoFfn4QF2ezypGnnmLpyy9XPw/w8CC3rMxuk97APcDB4GCWZ2dbHNNiDpr6AeOqugBk+/qyTKcjW6ezqO+jUnF3Rd1ZU0IIIUT74szn9+UxUehSMWsWDBoE3brBwoU2q6x99VWL5/UFUWAOpAB0dYIogA6Yt4GpfSNUAYRMmEBAZT6u2rQaTb3XEkIIIS41Eki1Fzt2wFdfmR+bTFBr1KmKPjubYidHhHwq/9TZOHYGmPTYY/guXQpVQZJCAQ89RMTw4dZt1Ul+KoQQQlzq2vUcqcvK//5n+fzsWSgvh1rb35x+4w2nmz0OrAKsN4mBmx9/nLBXXjE/2bULfv8dRo2C/v3p26sXBUeOsOvAAQyYI/J+Tzzh9PWFEEKI9kzmSDUjV86R2hMWRmBGBmHUin4vXICQkOo6v8bHc+j48SZdp0piv35cuWdPg/WKTpzg/G+/EThkCIGJiS65thBCCNGanPn8lhGpdkCXmcnvGRmYMP+DhWOe/O138WJ1IGWsqODUiRM2z48CTjt5TY2DgZ93bCyxsbFOti6EEEJcGmSOVFul18OBA5CXR/r771M1bFgBpAAlgC49vbr6hcWLKbMzuGhre+YrMac1sMfd378RnRZCCCEuLxJINYOFCxfSs2dPBg0a1LgGCgth5Ejo0wdiYkh9912rKl8A6Rs3Vj9Pr7VPXl0FNsoigbmAm41jAO4ycVwIIYRokARSzWD+/PkcOnSInTt32q6QkgI//2wOmGw4NHs2+7ZvxwBw8SJnao081VaWmWlewXfgAAlpadwMjACCKo8HAT2A4ZjTFlicC7ir1XTp0MFm2+7BwfZfoBBCCCEAmSPV8lauhEmTzAFQXBzs3g1eXtWHyzMyWP/zzxQD24EhwEU7TemysmD2bPjyS1RA58qfJMy3AGuPNv2B5T56OoDu3fFTq82T1utwtxNgCSGEEKKGjEi1tBdfNAdRAMePw4cfWhzee/vt1QFPPrASMNppqmD/foq//JK6M6MUgFudyeIedeqUAfTpg5+dTYbdQ0PreRFCCCGEAAmkWpRRp6No0ybLwv/7v+qH5dnZ7Fi50uKwfz3tHTp+nLeBt4Cvgc1VB8LD4corAVgG/Ark1DlXB9C3L77du1u1Gw94RUfX+1qEEEIIIbf2WtTFX3/lE8zbroRW/gwtKkJZWgpaLSeeeYaSOivv+gPr7LRXVbcUSKfWrbzRoyEkBBNwBKxGrKDWiFRJiUW5ApgMqGwEWEIIIYSwJCNSLSjz118B81ylU8ABKv8BKlffHfz2W4v6Xak/RUFdVZPMmTsXoqKowHYQBbAJKO7YEb8hQyzKTUAhmDdGFkIIIUS9ZESqBWXu2GHxvHoW0vLlFAUEkHrRclp5b8z5ohwVDOb9+MaNA50OfQP1Pfv3R6FS4Y7lXnv5gL+7uxNXFkIIIS5PEki1oAunLfOLlwFrgZz336frkSNWo0cXgW1OtB+0aBHMnAmALiSE9Q3UV6hUAPgCWbXKbeWdEkIIIYQ1CaSawcKFC1m4cCEGg6G6zKTTkVknb1R65Q/FxaStWGHVjjNBFEDQ+PHVj5XR0Ryqp+7gMWOqH/thGUjZyoQuhBBCCGsyR6oZ2ErIefHXXymv5xxDnefuiropNOvnq1ajqbWBsVutx7X5uLsTFx/PwM8/ry7zq3MbTwIpIYQQwjESSLWQqonmjlAC8fHxTrUfHBTUcCXg7rIyph49irZz5+oyz0mTqq/rDXg6eW0hhBDiciWBVAupO9G8Pv4aDd5OZhYPiopytkvV+r75Jh4qFUagCOg8d26j2xJCCCEuJzJHqoVkp6VZPNcqlZQabecs9w8IwLNTJ4syBXAt8JOd9oP79bMq81QqKbFzDYu+hIdzW0oKJ99+m+Dhwwm7+uoGzxFCCCGEjEi1mOI6iS+j6hlBCujcGa/wcIsyXyAOsJeUoION4Gfso49aPO9bT5JNr/Bw+r70kgRRQgghhBMkkGohJXrLrE6dhw2zW9c/Ph7POnvgVYVhY9zcqBuCeSgUBE+ebNVO7AsvEFc5FypAo2Hgf//rbLeFEEIIUQ+5tdcCjGVllNbZ+iV00iRUn39utVpPDQQkJOAVF2dRXg7ogd49enBu3z6LY+EdO1bnhKpN6eHB1PR09NnZqP39Uarln1sIIYRwJflkbQH6U6eskm169+qFn0bDxVojVROBXoDpxhttpkooATQhIebcU7V0Tkqq9/qa4GDnOy2EEEKIBsmtvRZQdvy4VZk2Lo6AwECLslxA4eGBsmtXNF26WEW5JYC+pITSOuXhN9zgyu4KIYQQwkEyItUM6mY2L0tJsTjuoVCg0mrxDw+H8+ery3MBYmJAqUQBVNRptwgIO3yY+4AczFnRzwGhM2Y00ysRQgghRH1kRKoZ1M5s3gEI9PLiWmAsMAxIqLzVFtCtm8V5eQB15kbV9jPAPfegwLxBcX/g6qAglLLBsBBCCNEqJJBqZtcD3jod8ZgDn2HA8N69AQhITLSomwuYoqPtthUAMGMG1J5Y/vHHLu2vEEIIIRwngVQziwUUWVmWhZVZywNGjLAoLgeKa210PKrO8VH33AOJibBxIzz9NKxcCVOmNEe3hRBCCOEAmSPVzHywH0j51BqR8gT8gfI+farL+r35JrkjRnCusJDYTp2Ifukl84GhQ80/QgghhGhVEkg1sxKAOtvDVAVSCpWK7vHxHDl2jBLApFLhfeON1dU0/foxPjMTLlyAzp1B8kAJIYQQbYp8MjezPwGa9estC2ttSDxu2zb85s6lNCuLAW+8gZuPj2VdrRYiIpq/o0IIIYRwmgRSzcwTKMO8R56iqrBWIOUeEMCIX35p+Y4JIYQQoskkkGoB72GeSO5Z+XN1YSEhrdslIYQQQriABFLNrGqPPIDiyh9VaGjrdUgIIYQQLiPpD5pZ3e1cADzrJOIUQgghRPskI1LNyGQyWQVSSkDn54e+oKA1uiSEEEKIBhRUfkabTKYG60og1Qyq9trT6/WctHH80YCAFu+TEEIIIZxTWFiIn59fvXUUJkfCLdEoRqOR+Ph4du3ahUKhaPiEBgwaNIidO3e2+PnOnudo/YKCArp06UJaWhq+vr5O9+ty0NR/85bSWv1szuu6sm15715e5H3betd1Vdsmk4kBAwZw7NgxlMr6Z0HJiFQzUiqVaDSaBqNZR6lUqib90mrs+c6e52x9X19f+WVsR1P/zVtKa/WzOa/ryrblvXt5kfdt613XlW1rNJoGgyiQyebNbv78+W2mrcae7+x5rnzNl7v28nfZWv1szuvKe1c0Vnv5e5T3rWvaklt7otUUFBTg5+dHfn5+u/j2JoQwk/euEDVkREq0Gnd3d5577jnc3d1buytCCCfIe1eIGjIiJYQQQgjRSDIiJYQQQgjRSBJICSGEEEI0kgRSQgghhBCNJIGUaJM2bNhAYmIiWq2WUaNGkZKS0tpdEkI44K233iIsLAx/f38eeeQRh7bYEKI9k0BKtEnffvstH330ERcvXqRHjx4899xzrd0lIUQD8vPzOX36NHv37uXIkSOsXLmStWvXtna3hGhWktlcNJvt27dTWFjI2LFjbR4vLi4mOTkZrVZLQkICKpWq+tj//ve/6sf9+/fn2LFjzd5fIQTk5uayY8cOYmJiiI2NtVnn5MmTnDlzhri4OMLDw6vL/fz8ePXVV6moqCAjIwOVSkVISEhLdV2IViEjUsLl3nvvPXr27Mm0adOYNm2azTpLly4lPDyc+fPnM336dHr27MnJk9ZbPO/du5fPP/+c//u//2vmXgtxeUtLS+P222+nV69eTJ8+nUWLFlnVqaioYObMmfTv358nn3ySuLg4Hn30UYs6Op0ODw8PunbtypAhQ+jdu3dLvQQhWoUEUsLlcnJy+P777+3ejsvJyWHWrFk8+uij7Nu3j9OnT9O1a1fmzZtnUW/dunX8+c9/5scffyQoKKgFei7E5SsjI4MrrriCEydOEBYWZrPOG2+8werVq9m/fz/bt29n48aNvPnmm/zwww/Vddzd3amoqCA7O5uzZ8/y2WeftdRLEKJVSCAlXO6vf/0rPXv2tHt88eLF6PV6HnroIQDUajWPP/44mzZtqh6V+uqrr3j++ef54YcfCAoKwmAwtETXhbhsJSUlcfvtt+Pp6Wm3zqeffsrMmTOJjIwEYODAgYwdO5ZPPvkEMH/5effdd8nNzaWkpASdToder2+B3gvReiSQEi1uz549xMbG4u3tXV2WkJAAmG/lATzzzDNs2rSJ0NBQPDw8GD9+fKv0VQhhptfrOXz4cPV7tUpCQkL1+3bYsGEcO3aM+Ph4hgwZQq9evaxGmoW41Mhkc9Hi8vLyCAwMtCirunWXm5sLwIkTJ1q8X0II+woKCjAajTbfu1XvW41Gw2uvvcZrr73WGl0UolXIiJRocRqNhtLSUouykpKS6mNCiLan6r1p670r71txOZNASrS4yMhI0tPTLcqqnlfNvRBCtC2+vr4EBgaSlpZmUZ6eni7vW3FZk0BKtLiJEyeSkZHB9u3bq8t++OEH/P39GTx4cCv2TAhRnwkTJvDTTz9VP9fr9SxdupSJEye2Yq+EaF0yR0q43J49ezh//jyHDh3CYDCwYsUKwDwR1dfXl6SkJG688UZmzpzJM888Q3Z2Nn//+99544035BaBEK2ktLSU9evXA+bbdSdOnGDFihUEBgaSlJQEwHPPPUdSUhJz585l0qRJfP755wDVK3CFuBwpTLIRknCxl156iY0bN1qVv/nmm8TFxQFQXl7OO++8w9q1a9FqtcycOZNrrrmmpbsqhKiUmZnJrbfealXet29fXnnllernR48e5c033+TMmTN069aNhx9+2G7eKSEuBxJICSGEEEI0ksyREkIIIYRoJAmkhBBCCCEaSQIpIYQQQohGkkBKCCGEEKKRJJASQgghhGgkCaSEEEIIIRpJAikhhBBCiEaSQEoIIRrJaDRy6tQpdDpda3dFCNFKJJASQohG2LZtGxEREcTExDB+/PjW7o4QopVIZnMhhHBSeXk5sbGxdOjQgYEDB/Lee+9RXl6OUinfTYW43Mi7XgghnPTLL79w5swZnn32WfR6PT4+PhJECXGZkne+EEI46dtvvyUgIIBJkyaxZ88eoqKiWrtLQohWIoGUEEI4wWg08ttvvzFu3DhKSkrYu3cvI0eObO1uCSFaiQRSQgjhhOTkZHJzcxk1ahSrVq3CYDAwduzY1u6WEKKVSCAlhBBO2LRpEwBDhgzh66+/Jjg4mAkTJrRyr4QQrUUCKSGEcMK2bdtwc3MjJCSEn3/+mdmzZ6PRaFq7W0KIViKBlBBCOGHPnj3ExMTw4YcfYjAYuP/++1u7S0KIViR5pIQQwkEGgwGtVsuIESPYvXs3U6ZM4bPPPmvtbgkhWpEEUkII4aDc3FwCAwPRarWUl5dz4MABunXr1trdEkK0Irm1J4QQDqraU6+0tJT77rtPgighhIxICSGEo0pKSvDy8iIoKIjjx48TEBDQ2l0SQrQyGZESQggHbd++HYBJkyZJECWEACSQEkIIh5SVlXHPPfcAsGrVKsrKyigsLOTKK6/knXfeaeXeCSFaiwRSQgjhgCeeeIKzZ8/y9ddfU1hYyOjRoxkyZAheXl7cdtttrd09IUQrkTlSQgjhgFOnTqHVaunUqROnT59m5cqVREREMGnSpNbumhCiFUkgJYQQQgjRSHJrTwghhBCikSSQEkIIIYRoJAmkhBBCCCEaSQIpIYQQQohGkkBKCCGEEKKRJJASQgghhGgkCaSEEEIIIRpJAikhhBBCiEaSQEoIIYQQopEkkBJCCCGEaCQJpIQQQgghGkkCKSGEEEKIRpJASgghhBCikSSQEkIIIYRopP8Hcp8ck1zegrIAAAAASUVORK5CYII=", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "presets = [\"ksz1\", \"ksz2\", \"ksz3\", \"ksz5\", \"ksz6\"]\n", "titles = {\n", " \"ksz1\": \"WebSky\",\n", " \"ksz2\": \"AGORA (unlensed)\",\n", " \"ksz3\": \"AGORA\",\n", " \"ksz5\": \"FLAMINGO L1_m9\",\n", " \"ksz6\": r\"FLAMINGO fgas-$8\\sigma$\",\n", "}\n", "colors = {\"ksz1\": \"#ff9900\", \"ksz2\": \"#8b0000\",\n", " \"ksz3\": \"#ff0000\", \"ksz5\": \"#332288\",\n", " \"ksz6\": \"#9467bd\"}\n", "styles = {\"ksz1\": \"-\", \"ksz2\": \"--\", \"ksz3\": \"-\",\n", " \"ksz5\": \"-\", \"ksz6\": \"--\"}\n", "markers = {\"ksz2\": \"o\", \"ksz6\": \"s\"}\n", "\n", "lmax = 6000\n", "spectra = {}\n", "for preset in presets:\n", " sky = pysm3.Sky(\n", " nside=2048, preset_strings=[preset], output_unit=u.uK_CMB\n", " )\n", " intensity = sky.get_emission(143 * u.GHz)[0].copy()\n", " spectra[preset] = hp.anafast(intensity.value, lmax=lmax)\n", " del intensity, sky\n", "\n", "for preset, cl in spectra.items():\n", " ell = np.arange(len(cl))\n", " dl = ell * (ell + 1) * cl / (2 * np.pi)\n", " plt.loglog(ell, dl, color=colors[preset],\n", " linestyle=styles[preset], linewidth=3,\n", " marker=markers.get(preset), markevery=300,\n", " markersize=4, markerfacecolor=\"white\",\n", " label=titles[preset])\n", "\n", "plt.xlim(50, lmax)\n", "plt.ylim(0.05, None)\n", "plt.xlabel(r\"$\\ell$\")\n", "plt.ylabel(r\"$D_\\ell^{\\rm kSZ}\\;[\\mu\\mathrm{K}^2]$\")\n", "plt.legend()" ] }, { "cell_type": "markdown", "id": "6dddb960", "metadata": {}, "source": [ "The two purple curves are the new FLAMINGO models. The `fgas-8sigma` run has stronger baryonic feedback, which distributes gas more broadly and reduces small-scale kSZ power. Users can now evaluate this modeling uncertainty by changing the PySM preset from `ksz5` to `ksz6`." ] }, { "cell_type": "markdown", "id": "4e31dfa6", "metadata": {}, "source": [ "## Cosmological parameters for the `c1` model\n", "\n", "The `c1` preset generates a lensed CMB realization from unlensed CAMB spectra. Version 3.4.6 documents the original CAMB [`params.ini`](https://github.com/galsci/pysm/blob/3.4.6/docs/cmb_c1_params.ini) and compares it with the Planck 2015 TT,TE,EE+lowP best fit. The principal input values are summarized below." ] }, { "cell_type": "code", "execution_count": 4, "id": "c3ca9caf", "metadata": { "execution": { "iopub.execute_input": "2026-07-21T02:44:16.252410Z", "iopub.status.busy": "2026-07-21T02:44:16.252086Z", "iopub.status.idle": "2026-07-21T02:44:16.279611Z", "shell.execute_reply": "2026-07-21T02:44:16.278406Z" } }, "outputs": [ { "data": { "text/markdown": [ "| CAMB key | meaning | c1 value |\n", "|---|---|---:|\n", "| `ombh2` | $\\Omega_b h^2$ | 0.02225 |\n", "| `omch2` | $\\Omega_c h^2$ | 0.1198 |\n", "| `hubble` | $H_0$ [km/s/Mpc] | 70 |\n", "| `re_optical_depth` | $\\tau$ | 0.06 |\n", "| `scalar_spectral_index(1)` | $n_s$ | 0.9645 |\n", "| `scalar_amp(1)` | $A_s$ | 2.1e-9 |\n", "| `l_max_scalar` | maximum $\\ell$ | 2200 |" ], "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "c1_parameters = {\n", " \"ombh2\": (r\"$\\Omega_b h^2$\", \"0.02225\"),\n", " \"omch2\": (r\"$\\Omega_c h^2$\", \"0.1198\"),\n", " \"hubble\": (r\"$H_0$ [km/s/Mpc]\", \"70\"),\n", " \"re_optical_depth\": (r\"$\\tau$\", \"0.06\"),\n", " \"scalar_spectral_index(1)\": (r\"$n_s$\", \"0.9645\"),\n", " \"scalar_amp(1)\": (r\"$A_s$\", \"2.1e-9\"),\n", " \"l_max_scalar\": (r\"maximum $\\ell$\", \"2200\"),\n", "}\n", "rows = [\"| CAMB key | meaning | c1 value |\", \"|---|---|---:|\"]\n", "rows += [f\"| `{key}` | {meaning} | {value} |\"\n", " for key, (meaning, value) in c1_parameters.items()]\n", "display(Markdown(\"\\n\".join(rows)))" ] }, { "cell_type": "markdown", "id": "618b1a7e", "metadata": {}, "source": [ "The primary density and spectral-index parameters match the documented Planck 2015 values. Some secondary values differ slightly, including `hubble = 70`, `tau = 0.06`, and `A_s = 2.1e-9`; these values are now explicit and can be recorded alongside simulations that use `c1`.\n", "\n", "See the [PySM model documentation](https://pysm3.readthedocs.io/en/latest/models.html) for details, and report problems on the [issue tracker](https://github.com/galsci/pysm/issues)." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.13.6" } }, "nbformat": 4, "nbformat_minor": 5 }