{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# One Fluid Dispersion Relation\n", "author: Louis Richard\n", "\n", "Solves one fluid dispersion relation" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Load IGRF coefficients ...\n" ] } ], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "\n", "from pyrfu.dispersion import one_fluid_dispersion" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Local conditions" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Magnetic field" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "b_0 = 10e-9" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Angle of propagation" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "theta = 5.0" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Particles" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "ions = {\"n\": 10e6, \"t\": 10, \"gamma\": 1}\n", "electrons = {\"n\": 10e6, \"t\": 10, \"gamma\": 1}" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Solve one fluid dispersion relation" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "wc_1, wc_2, wc_3 = one_fluid_dispersion(10e-9, 5, ions, electrons)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAqEAAAImCAYAAABjF1M+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAAxOAAAMTgF/d4wjAACUjklEQVR4nOzdd1hV9QMG8Pfey94IIoqgiHsLam7cONK0zFGusuFOU0sty1X2KzOzYW5FTXMPVBwIuPdWhlsciAiy4a7z+wMlCVRALt873s/z3Ac459xz3ntReDnje2SSJEkgIiIiIipBctEBiIiIiMj0sIQSERERUYljCSUiIiKiEscSSkREREQljiWUiIiIiEocSygRERERlTiWUCIiIiIqcSyhRERERFTiWEKJiIiIqMSxhBLpgeXLl0Mmk+U8bG1tUbFiRfTs2RPr1q3Df29sNnXqVMhkMkFpX58h5n+dzFu2bMGcOXOKdZ2i3L9/H4MGDYKrqyvs7e3Rp08fPHnyRHSsVwoLC8v1f+zZw8nJqUDPj4mJQa9eveDo6AgHBwe8/fbbuHPnzmsvS2TKWEKJ9Mj69etx9OhR7Ny5EzNmzIClpSX69euHDh06ICMjI2e5jz76CEePHhWY9PUYev7CelEJNbT34ebNm2jcuDGSk5OxevVqzJ8/H7t378bIkSNFRyuwefPm4ejRozmPffv2vfI56enpaNu2LSIjI7FixQqsXLkSV69eRZs2bZCWllbkZYlMnZnoAET0r/r166Ny5co5Xw8YMADvvvsu3n33XXzxxRf47bffAADly5dH+fLlRcV8oaysLFhaWr5yOX3JX9C8uqIv70NBSJKEfv36oX79+ti0aVPOHtzo6Gj873//w+LFi2FlZSU45avVqFEDTZo0KdRzFi1ahBs3biAqKirn/2fdunVRpUoVLFiwAJ9//nmRliUyddwTSqTn3nnnHbz11ltYtGgR0tPTAeQ9jBsdHY2ePXvCzc0NVlZW8PLywrvvvgu1Wp2zzLPnXLx4EW3atIGNjQ3Kli2Lb775BlqtNtc2z58/j+7du8PZ2RnW1tZo3rw5Dh48mGuZZ+u7dOkSAgICYGdnh969excqy/OCg4PRtGlTWFtbw9HRET169EBUVFS+27x69Sq6du0KOzs7VKhQAdOnT8/zGv7rRXkL85r/69q1axgwYAC8vb1hbW2NSpUqYdiwYUhMTMxZZvDgwVixYgXu3buXcwi4YsWKJfI+FOR7UVCbN2/G8ePHMWfOnFyZvby8oFQqcf/+/UKv01Bs27YNTZo0yfUHore3N5o3b46tW7cWeVkiU8cSSmQAunTpgqysLJw6dSrf+V27dsW9e/dyDo/+8MMPsLS0zLeY9ejRA+3bt8eWLVvw3nvvYcaMGZg+fXrO/DNnzqBZs2ZISEjAokWLsHHjRri4uKB9+/Y4ffp0nvW99dZb8Pf3x7Zt2zB27NhCZXkmODg4p0z9888/mD9/Pi5duoQWLVrg3r17eZbv2bMn2rZtiy1btqBHjx749ttvsWLFioK8lXnyFuU1P3P//n14enpi7ty52L17N7755huEhISgS5cuOctMmTIFXbp0QenSpXMOAW/evLlE3oeCfC9u3boFmUyGqVOnvvR9W7p0KZo2bYpKlSpBrVbnPFJTUwEAZma6ObAmSVKu7b3oodFoCrS+999/HwqFAi4uLnjvvfcKdK7m5cuXUbt27TzTa9WqhStXrhR5WSKTJxGRcMuWLZMASFevXs13fnBwsARAWrt2rSRJkvTtt99Kz/77Pnr0SAIgbd269aXbePacWbNm5Zr+0UcfSXZ2dlJiYqIkSZLUtm1bqXr16lJWVlbOMmq1Wqpevbr01ltv5Vnf3Llzc6YVNsszfn5+UuXKlSWVSpUz7caNG5KZmZk0duzYPM9bunRprvXVrl1b6tChQ4G2+XzeZwrymv+bOT8qlUo6ePCgBEA6c+ZMzvRBgwZJHh4eL8z0THG+DwX9Xty6dUtSKBTStGnTXrhMVlaWZGNjIwHI92Fubp4rc3EKDQ194Xaff/j7+790PWfOnJHGjRsnbdu2TQoLC5N++eUXqXTp0lK5cuWkhw8fvvS55ubm0pdffpln+ldffSUpFIoiL0tk6nhOKJEBkJ5eHZ/fldQuLi6oVKkSJk6ciIcPH6J169aoUqXKC9f1/CFoAOjbty8WL16MS5cuwc/PD+Hh4Zg8eTLkcnmuw7bt27fH6tWr86yvZ8+eRc4CAGlpaThz5gwmT56ca2/as0OY4eHheZ7TtWvXXF/Xrl0bZ8+efel28ssLABkZGYV+zc8olUrMnj0bgYGBuH37NjIzM3PmRUVFoUGDBgXKBBT/+1DQ70WFChVeeXj+ypUrSE9Pxx9//IHGjRvnmte3b184OzvnZE5JSYGHhwcSExOhUChyljtz5gxatGiBqlWrQqvVIikpCZ9++ikmT5780m37+fnh5MmTL10GAOzt7V86v0GDBrm+H/7+/mjVqhUaN26MefPmYebMma/cBhEVLx6OJzIAMTExAICyZcvmmSeTybB37140bNgQkyZNQtWqVVGpUiXMnz8/33WVKVMm36/v3buHhIQEaDQazJgxA+bm5rkev//+OxITE/McVn8+U2GzAEBiYiIkScr3tbm7uyMhISHP9FKlSuX62tLSMlcBfJn/bqcor/mZSZMmYerUqejfvz927NiBEydOYNOmTQBQ4DzPFPf7UJTvxYvcunULANCiRQs0bNgw5+Hp6YmbN2+iU6dOOcuePn0aDRo0yFVAn03v1KkTzp07hwsXLuDvv//GtGnT8gw/9l92dnaoX7/+Kx/Pn4NZUL6+vqhateorS66zs3Ou83yfSUhIgLOzc5GXJTJ13BNKZAB27NgBKysr+Pn55Tu/UqVKCAwMhCRJOH/+PH7//XcMHz4cFStWROfOnXMt+/DhQ1SqVCnX1wDg4eEBJycnyOVyjBgxAgMHDsx3W3J57r9d/7t3tjBZgOxf2jKZDLGxsXnmxcbG5ilar+u/eYvymp9Zu3YtBg4ciK+//jpn2rNzJAtLF+9DYb8XL/JsT+l/i+XKlSshk8kwePDgnGknTpzAG2+8AQA4fPgwPvvsM0ybNg2nTp3K2Yuq1Wpx4MABtG3b9pXjpIaHh6NNmzavzOjv74+wsLACv6bnvSpDrVq1cPny5TzTr1y5gpo1axZ5WSJTxz2hRHpu48aN2LZtG4YOHQobG5uXLiuTyVC/fv2cMSkvXbqUZ5l169bl+nrt2rWws7NDnTp1YGtri5YtW+L8+fPw9fXNtdfr2aOgCpIFAGxtbeHn54f169fnurjk9u3bOHLkCFq3bl3gbRbF67zm9PR0mJub55q2bNmyPMtZWlrmGuf1RTl09T4U9HvxIs+u5n++XMXGxuJ///sfPvnkE/j4+ORMP3nyJBo2bIjp06djypQp2Lp1K7p27YrTp09j8eLFqF+/PkqXLo2NGzdi/fr1r9z2s8Pxr3osWLCgUK8JAE6dOoWoqKg8pxj8V/fu3XHs2DHcuHEjZ9qtW7dw+PBhdO/evcjLEpk67gkl0iPnzp1DfHw8lEol7ty5g6CgIKxfvx4dOnTArFmz8n3OhQsX8Nlnn6FPnz6oXLkyNBoNli9fDjMzM7Rt2zbP8osWLYJWq0WjRo2we/duLF68GFOnToWjoyMAYM6cOWjVqhUCAgIwZMgQlC1bFvHx8Thz5gw0Gg1++OGHF+YvbJZnZsyYga5du+LNN9/E8OHDkZqaim+//RaOjo4YN25cId/Fwivqa+7UqRNWrFiBOnXqoHLlyti0aROOHDmSZ7maNWsiISEB8+fPR8OGDWFlZYU6derkWa4434eCfi9u374NHx8ffPPNN/jmm2/yXZefnx9q1KiBSZMmwcrKCllZWZgyZQoqV66M2bNn51r2xIkTuHjxIt5++23s3bsXCoUCSqUSFy9exMOHD+Hk5ISsrCwEBARg0aJFOSMUvIi9vX2h/vh5kffffx/e3t7w9fWFk5MTzp49i1mzZsHDwwOjR4/OWS48PBzt2rXD0qVLc/aMf/zxx/j999/x1ltvYebMmZDJZJgyZQo8PT3x6aef5tpOYZYlMnUsoUR65N133wUAWFlZwc3NDb6+vli7di169er1wkOG7u7u8PLywpw5c3D37t2cghMUFJTv4futW7di1KhRmDFjBhwdHfH1119jypQpOfN9fX1x8uRJTJs2DaNHj0ZSUhJKly4NX19fDB069KX5C5vlmU6dOmHHjh2YNm0aevfuDQsLC7Ru3Ro//vgjypUrV5C37rUU9TX/9ttvkCQJX331FYDsobTWrFmTZ8/aRx99hGPHjmHy5Ml48uQJKlSokHOe5fOK830o6PdCkiRoNJqXDqElk8mwefNmfPzxx+jduzdcXV3Rv39/TJkyJdfe+bi4OKjVajRt2hQpKSk5h+8vXrwINze3nFtkWlpaolatWvkOO6UrtWvXxpo1a/Dbb78hPT0d7u7uePvttzFt2jS4urrmLJff+2Fra4v9+/dj7NixGDBgACRJQrt27TB37lzY2dnl2k5hliUydTLpVWeF69jRo0fx1Vdf4eTJk1AoFKhZsyYOHTr0wnOwiKhopk6dimnTpkGlUulsTEcybUFBQViwYAE2b96Mbt26oXnz5vj666+xYMECbN26FTt37gSQfaFd8+bNERgYqPPTLYhIfwn9TXT06FF07twZ8+bNQ1BQECwsLHD69OlXniRORET65+TJk/D19YWZmRk2bNiANm3awN3dHadPn8bx48dRv359mJmZwdzcHD///DMLKJGJE7ontGXLlmjcuDF+/vlnURGITAb3hBIRkT4Rdsw7PT0dR44cgUKhQOPGjeHi4gI/Pz9s3LhRVCQiozZ16lRIksQCSkREekHYntC7d+/C09MTbm5uCAoKQoMGDbBt2zb07dsX4eHhaNq0ac6yWq0W9+/fh729PQ/VExEREekZSZKQkpKCcuXKFfi6HmG7RJ7dYm3w4MFo1KgRAODtt99GmzZtsGXLllwl9P79+/D09BSSk4iIiIgKJiYmBuXLly/QssJKqKOjI3x8fAq0Z/NZYY2JiYGDg4Ouo5mkyZMn4/vvvxcdw6jp6j1WKpXw9PREeHg4qlevXuzrNxT8N6x7fI91i+9v8bvw6ALGh4/H21Xexu21t1843jK9vuTkZHh6euZ0toIQenLYqFGjMGvWLPTt2xd169ZFUFAQwsPDMX369FzLPSuqDg4OLKE6YmFhwfdWx3T1Hh87dgy2trZo1KiRSZ+uwn/Dusf3WLf4/havxxmPMeH4BIxtNhZ9qvfB55s/5/tbAgrze0hoCf3ss8+Qnp6Obt264cmTJ6hSpQr++eefnPsOE9GrHTx4EC1atDDpAkpE9F8u1i7Y0H0DPOw8REehFxB+meykSZMwadIk0TFMXkBAgOgIRk9X7/Gz2wyaOv4b1j2+x7rF97f4PV9A+f7qH+F3TCqI5ORkODo6IikpibvSiZ6j0WhQqlQphIaGwtfXV3QcIiIyUUXparw3JpEBO3fuHGQyGerVqyc6ChGRMAawP43yIfxwPBEVXXh4OFq2bAmFQiE6ChGZsMzMTCiVSiHbfpL5BF8d+grD6w9HLddaQjKYCgsLC1hZWRXb+lhCiQxYWFgY/P39RccgIhOWmZkJb29vxMbGCs2xCquEbt8UuLu74+bNm8VWRFlCiQyURqPBwYMHMWXKFNFRiMiEKZVKxMbGcixvI/dsHFClUskSSmTqLl68CI1GgwYNGoiOQkTEsbyp0HhhEpGBCg8PR4sWLWBmxr8liYjI8LCEEhkong9KRESGjCWUyABptVocOHCAJZSIiAwWSyiRAbp8+TKysrLg5+cnOgoREVGRsIQSGaDw8HA0b94c5ubmoqMQEREVCUsokQHi+aBERGToWEKJDIwkSTwflIjIRPTu3RsVKlR44fxr167BwsICQ4cOLdD6IiIiIJPJsGfPnuKKWGQsoUQGJiIiAqmpqWjUqJHoKEREpGO1atVCTEwMUlNT850/adIkWFtbY/r06QVa3+nTpwEADRs2LLaMRcUSSmRgwsPD0axZM1hYWIiOQkREOla7dm1IkoSIiIg8844dO4YNGzZg0qRJcHNzK9D6Tp8+DW9vb5QqVaq4oxYaSyiRgeH5oEREJefcuXMwNzfHkiVL8p0/bdo0yGQyxMXF6WT7tWvXBgBcuXIlz7wJEyagQoUKGDNmTIHXd+bMmVx7QZOTk9GrVy+4u7sjPDz8tfMWBksokQGRJAnh4eEsoUREJWTEiBGoVasWPvzww3zn16hRAwBw4sSJPPMkSYJarX7lQ6PRvHD7lStXhqWlJS5fvpxr+pYtW3Do0CHMmjWrwPdylyQJZ8+ezSmhFy9eRMOGDXHv3j2cPn26xH+3sIQSGZDo6Gg8efIEjRs3Fh2FiKjIlBolUpWp+T402vwLmUareeFzlBqlTnLu3r0bR44cweTJkyGTyfJdxsvLCwBw+/btPPPCw8Nhbm7+yke7du1emEGhUKB69eq59oSq1WpMnDgRb7zxBvr27Vvg1xMdHY2UlBQ0bNgQgYGBaNKkCdq3b4/w8HB4eHgUeD3FhTedJjIg4eHhaNq0aYH/6iUi0keLLy7G/PPz8523uftmVHaunGf6zaSb6LmtZ77PGVZvGIbXH16sGQFg8eLFcHR0RM+e/243MDAQ1atXz9kZ8OyCIa1Wm+f5fn5+OHny5Cu3Y29v/9L5tWvXxpEjR3K+XrRoEaKionD48OE85fjSpUto0KABvvrqK0ydOjXXvDNnzgAAFixYgO3bt+Ovv/7CwIEDX5lPV1hCiQwIzwclImPwUZ2PMLBm/uXH2sw63+nejt442u9ovvMsFLq5UDMsLAzNmjXLuTFIZmYmPvnkEyxatCinhMbGxmbn8/bO83w7OzvUr1//ldt50V7WZ2rVqoW///4b6enp0Gq1mDZtGt599100a9Ysz7KjR4+Gv78/zp07l2fe6dOn4eTkhI0bN2LUqFFCCyjAw/FEBoPngxKRsbBQWMDOwi7fh0KuyPc5Crnihc/RRQlNSkpCfHw8fHx8cqaFhoYiKysL5cuXzzXNzMwMzZs3z7OO4jgcD+S+Qv7HH3/EkydP8L///S/Pcv/88w80Gg2+++47nD17Ns/806dPIyAgAPPmzcO8efOwadOmwrwlxY57QokMxPXr1xEfH48mTZqIjkJEZPTUajWA7L2fz/z111+wsrKCUpl9Dmp8fDzWrVuHbt26wdnZOc86ivNwPACEhIRgzpw5GDVqVJ49r2lpaZg4cSK2bt0KHx8f3L17F4mJiTm5nl2UNHHiRAwfPhxRUVEYMGAAPD09hY07zRJKZCDCw8PRuHFjWFvnf6iKiIiKj4uLCzw9PbFx40Z07NgR58+fx9WrV9G5c2f88ccfyMrKwnfffQcAmDVrVr7rsLe3L5ZB4StWrAhbW1t8++23sLOzw9dff51nme+++w5du3ZF3bp1AQA+Pj44e/Ys2rZtCyB7R0ZSUlLO6QFz5szB9evX0b17dxw/fjznAquSxMPxRAYiLCwMrVu3Fh2DiMhkLF++HK6urhg4cCAOHjyIrVu34osvvsClS5fQp08fWFlZISwsDNWqVdNpDplMhpo1ayIzMxNTp06Fo6NjrvnXrl3DvHnzsGXLFlSsWBEVK1bEvXv3cp0X+uxOSc9KqEKhwNq1a1GmTBl07doVycnJOn0N+ZFJkiSV+FYLKTk5GY6OjkhKSoKDg4PoOEQlTpIkVKhQAUuXLkX79u1FxyEiysHf0eK9+eab6Ny5M0aMGJEzbebMmYiKisLKlSuLZRuv+j4X5d8B94QSGYAbN24gNjY23yshiYjIdAUFBeH69ev49NNPc02vVatWvlfI6xOWUCIDEBoaiiZNmsDGxkZ0FCIi0iNvvvkmIiIiYGaW+zKfnj174uLFi4JSFQxLKJEBCA0NRZs2bUTHICIiKjYsoUR6TpIkllAiIjI6LKFEei4qKgoJCQkcH5SIiIwKSyiRngsNDUWzZs14v3giIjIqLKFEei40NDRnsGEiIn1z5uEZ0RHIQLGEEukxSZIQFhbG80GJSC/tvb0X48LGiY5BBoollEiPXb58GWlpacLu60tE9DKH7h3C9ObTRccgA8V7xxPpsdDQULRo0QIWFhaioxAR5TGt2TQht3sk48A9oUR6bP/+/TwUT0RERokllEhPabVahIeH86IkIiIySiyhRHrq/PnzUKvV8PX1FR2FiExcQmYCHqU/Eh2DjAxLKJGeCg0NRatWrfLcD5iIqCTFpMRgwM4BWBO5RnQUKoKIiAjIZDLs2bNHdJQ8WEKJ9BRv1UlEol15fAUDdg5Aq/KtMLLBSNFxqAhOnz4NAGjYsKHgJHmxhBLpIbVajfDwcJZQIhLm6P2jGLJ7CAbVGoQvGn0BuYyVwRCdPn0a3t7eKFWqlOgoefBfFJEeOnPmDMzMzFCvXj3RUYjIBO24sQOfhX6GyW9Mxge1P4BMJhMdSZhz587B3NwcS5YsyXf+tGnTIJPJEBcXp9Mc69evh7+/PxwcHGBnZ4dGjRphx44dr3zemTNncu0FTU5ORq9eveDu7o7w8HBdRn4lllAiPRQaGgp/f38oFArRUYjIBD1Ie4C5reeim0830VGEGzFiBGrVqoUPP/ww3/k1atQAAJw4cSLPPEmSoFarX/nQaDQvzTBu3Dj07dsXtWrVwqpVq7BmzRo0a9bslWO0SpKEs2fP5pTQixcvomHDhrh37x5Onz4Nf3//grwFOsMrHoj0UGhoKLp06SI6BhGZqI/qfKTT9WtUWmjU2nznmVkqIJfn3fOq1UpQZ+Vf1hRmcijMi3+/2u7du3HkyBH8888/L9wb7OXlBQC4fft2nnkFPa3K398fYWFh+c5bu3Yt5syZg/Xr16NXr14507t1e/UfCNHR0UhJSUHDhg0RGBiIYcOGYdCgQZg7d65e3ASFJZRIz6hUKhw6dAg//fST6ChERDpxOvgWTu64le+8vt80hks5uzzTE2PTsHZ63r2NANCoa0U07lapOCMCABYvXgxHR0f07NkzZ1pgYCCqV6+Oxo0bAwBSU1MBZI/t/F9+fn44efLkK7djb2//wnnffvstunXrlquAFtSZM2cAAAsWLMD27dvx119/YeDAgYVej66whBLpmZMnT8La2hq1atUSHYWISCf8OlVE/fZe+c4zs8z/NCRnd1t8/EurfOcpzHRzdmFYWBiaNWsGc3NzAEBmZiY++eQTLFq0KKeExsbGAgC8vb3zPN/Ozg7169d/5XZetJf1xo0biI6Oxtdff12k/KdPn4aTkxM2btyIUaNG6VUBBXhOKJHeeXarTrmc/z2JSLcepT9CdGJ0iW9XYS6HhbVZvo/8DsUDgFwue+FzdHEoPikpCfHx8fDx8cmZFhoaiqysLJQvXz7XNDMzMzRv3jzPOsLDw2Fubv7KR7t27fLNcP/+fQBAuXLlXpgzNDQUzZs3h6+vL6pWrYrffvstZ97p06cREBCAefPmYd68edi0aVOh3wdd4p5QIj0TGhpapMMuRESFcTPpJobtG4a2Xm3xRaMvRMfRO2q1GkD23s9n/vrrL1hZWUGpVAIA4uPjsW7dOnTr1g3Ozs551vG6h+Oflc/Lly+/sKj269cPp06dQvny5SFJEp48eQLg34uSJk6ciOHDhyMqKgoDBgyAp6cnGjVq9MpMJYEllEiPZGVl4ciRI/jjjz9ERyEiI3bh0QWMCBmBd6q8g898PxMdRy+5uLjA09MTGzduRMeOHXH+/HlcvXoVnTt3xh9//IGsrCx89913AIBZs2bluw57e/vXGiTe29sbLVq0wLfffgsAqFu3Lh4+fIidO3fiq6++QtWqVVG+fHmMHDkSffv2xZtvvplThq9fv46kpKSc0wHmzJmD69evo3v37jh+/HjOBVUi8XgfkR45duwYnJycUK1aNdFRiMhIHbh7AB/v+RhD6w3FGL8xJj0G6KssX74crq6uGDhwIA4ePIitW7fiiy++wKVLl9CnTx9YWVkhLCxMZz+zZTIZNm7ciF69emH27NkICAjAhAkToNVqc85BPXbsGEaPHo2DBw+iSpUqSEpKAvDvnZKelVCFQoG1a9eiTJky6Nq16yuHdyoJMkmSJNEhXiU5ORmOjo5ISkqCg4OD6DhEOvPNN9/g+vXrWL16tegoRGSENl/djFknZmF68+noVLFTsayTv6PFiYyMRNWqVSGXy3Hv3j3Uq1cPt2/fhq2tbbFv61Xf56L8O+DheCI9EhISgiFDhoiOQURGytrMGn+0+wON3PXjnEB6Pb/88gtCQ0Nha2sLKysrrFq1SicFVFdYQon0REpKCk6cOIE1a9aIjkJERqqTd/Hs/ST9sGDBAtERXgvPCSXSEwcOHEDFihX14mRxIiIiXWMJJdIT+/bte+EQHERERMaGJZRIT4SEhKB9+/aiYxCREYhNi8WR+0dExyB6KZZQIj3w8OFDXL58GW3atBEdhYgM3LXEa3h/5/sIjwkXHYXopVhCifTA/v37Ua9ePbi4uIiOQkQG7PTD0xgYPBC9qvTCxMYTRccheileHU+kB0JCQng+KBG9ln2392HyocmY0GgC3q36rug4RK/EPaFEeoDngxLR69h0dRMmH5qMH1r+wAJKBoN7QokEu3HjBu7du4cWLVqIjkJEBqqSYyUs7LAQ9d3qC8ugD7eBJN3RxfeXJZRIsH379qFp06YGdZcLItIvIsunhYUF3N3d4enpKSwDlQx3d3dYWFgU2/pYQokE4/mgRGTIrKyscPPmTSiVStFRSMcsLCxgZWVVbOtjCSUSSKvVYv/+/fjss89ERyEiKjIrK6tiLSdkGnhhEpFAFy5cQGZmJho1aiQ6ChEZgLspd7Hp6ibRMYiKBUsokUAhISHw9/eHubm56ChEpOciEyIxYNcARCdGQ5Ik0XGIXhsPxxMJFBISgg4dOoiOQUR67viD4xgTOgYf1fkIH9b+EDKZTHQkotfGPaFEgiiVShw4cIDjgxLRS+26uQuj9o/CxMYTMaTOEBZQMhrcE0okyPHjx2Fra4vatWuLjkJEemrLtS2YdXwW5rSegxYeHEuYjIuwPaFTp06FQqGAnZ1dzqNfv36i4hCVuJCQELRt25Z7NYjohRq5N8KSgCUsoGSUhO4Jbdq0KQ4dOiQyApEwISEhGDRokOgYRKTHPOw84GHnIToGkU7wnFAiAVJTU3Hs2DGeD0pERCZLaAk9e/YsSpcujQoVKuC9997DzZs3RcYhKjEHDhyAl5cXKlasKDoKERGREMJKaK9evXDlyhXExcXhyJEjkMlkaN++PVJTU0VFIioxvFUnET3vdvJtzDszj+N/kkkRdk7o81cEe3h4YOnSpXB0dMSRI0fQsWPHfJ8zefJkWFhYAAACAgIQEBBQIlmJitu+ffswefJk0TGISA9cfHQRI0JGoEflHpAgQQZerEiGYffu3di9ezeA7GEHC0sm6cmfXUqlEo6OjtiyZUuecpmcnAxHR0ckJSXBwcFBUEKi4hEbGwsPDw/ExcXBxcVFdBwiEujg3YMYHz4eIxuMxICaA0THISqyonQ1YYfj161bh/j4eADAw4cP8dFHH6FMmTJo1qyZqEhEJWLfvn1o0KABCyiRidtybQvGhY/D1GZTWUDJJAkroatWrUKNGjVgY2MDX19fqFQq7Nu3D/b29qIiEZWIvXv3vvCUEyIyDbtu7sL/TvwPv7X9DZ29O4uOQySE3hyOfxkejidjIUkSPDw8sHr1arRp00Z0HCISJFWZigdpD1DFuYroKETFoihdjbftJCpBly9fRlJSEk87ITJxdhZ2qGLBAkqmjYPVE5WgPXv2wN/fH5aWlqKjEBERCcUSSlSCeD4oERFRNpZQohKSmZmJ8PBwdOjQQXQUIioh159cx8SDE6HSqERHIdI7LKFEJeTIkSNwdnZGzZo1RUchohJwNu4sBu4aiPJ25WEm5yUYRP/F/xVEJWTPnj3o0KEDZDLeDYXI2IXcCcGkg5MwvuF49K7WW3QcIr3EPaFEJYTngxKZhnVR6zDp4CTMajGLBZToJbgnlKgEPHr0COfOnUP79u1FRyEiHQqPCcevZ37FX+3/gm8ZX9FxiPQaSyhRCQgJCUHdunXh5uYmOgoR6VBzj+b4581/UN6+vOgoRHqPh+OJSsDevXt5VTyRCTCTm7GAEhUQSyiRjkmShD179vB8UCIiouewhBLpWFRUFOLj49GiRQvRUYioGEmSJDoCkUFjCSXSsT179qBVq1awsrISHYWIiklUQhQ+2P0BUpWpoqMQGSyWUCId4/mgRMblxIMT+CD4AzQv1xy25rai4xAZLJZQIh1SKpUIDQ1lCSUyEsG3gjFy/0hMaDQBH9f9mDefIHoNHKKJSIeOHTsGOzs71KlTR3QUInpNq66swm9nf8Ns/9loVb6V6DhEBo8llEiH9uzZg/bt20Mu50EHIkN2Nu4sFl5YiMUdF6NOaf5RSVQc+JuRSId4q04i49DArQG29tjKAkpUjFhCiXQkISEBp06d4q06iYyEs5Wz6AhERoUllEhHQkJCULNmTZQrV050FCIiIr3DEkqkI8HBwejcubPoGERUSBqtRnQEIpPAEkqkA5IkITg4GAEBAaKjEFEhXI6/jB5beyA+I150FCKjxxJKpAOXLl3CkydPeKtOIgNy+N5hDNkzBL2q9oKrtavoOERGj0M0EelAcHAw2rZtC0tLS9FRiKgAtl/fjhnHZuDbpt+ia6WuouMQmQSWUCId2L17N3r27Ck6BhG9giRJWHZ5GRZdWIRf2/yKpuWaio5EZDJ4OJ6omKWmpuLgwYPo1KmT6ChE9Ao3km5gdcRqLA1YygJKVMK4J5SomIWFhcHT0xM+Pj6ioxDRK/g4+WBHzx2wMrMSHYXI5HBPKFExCw4O5l5QIgPCAkokBksoUTFjCSUiIno1llCiYnTt2jXExMSgdevWoqMQ0X+kq9JFRyCi57CEEhWj4OBgtGzZEnZ2dqKjENFzzsWdQ5dNXXA7+bboKET0FEsoUTHavXs3D8UT6ZnQO6H4ZO8nGFpvKCo4VBAdh4ieYgklKiZZWVnYv38/b9VJpEc2RG/Alwe/xPctvkff6n1FxyGi53CIJqJicujQITg5OaF27dqioxCZPEmS8NeFv7DyykrMbz8ffmX8REciov9gCSUqJs+uipfJZKKjEJm8h+kPsefWHqzotAJVnKuIjkNE+eDheKJiwqGZiPSHu607NnbfyAJKpMdYQomKwd27dxEREYH27duLjkJET8ll/BVHpM/4P5SoGOzZsweNGzeGs7Oz6ChEREQGgSWUqBjwUDyROPEZ8aIjEJk8lUpT6OewhBK9JrVajb1797KEEglwMvYkum/ujkvxl0RHITJZwaG38MvXhwr9PF4dT/SaTpw4AYVCAT8/DgFDVJL23NqDrw9/jYmNJ6K2K4dGIyppF6IeY+uKy7BPVMGlsmOhn88SSvSagoOD0bFjRygUCtFRiEzG3xF/49czv+KnVj/B39NfdBwik/IkXYm5e6Oh2BMLa3cb9BrZAA52Ej4fX7j1sIQSvaZdu3Zh1KhRomMQmQRJkjDv7DxsiN6AhR0Xol7peqIjEZkMlUaLVcduY+6+q6jv6YSJXzZEDS8nAEBycnKh18cSSvQaYmNjcebMGZ4PSlRCUlQpuBR/CYGdA+Ht6C06DpFJkCQJoVFxmLkjAnKZDHP71kebam6vvV6WUKLXEBwcjIYNG8LN7fX/MxLRqzlYOGBRx0WiYxCZjHMRj/B76DWcik/B2PZV8d4bXjBXFM917SyhRK9h586d6Nq1q+gYRERExSo2Ph0rllyA1c00VKntiJ/Ht4GjjXmxboMllKiIVCoVdu/ejS+++EJ0FCIiomKRmaXGitWXkXIqHpKdGVoNr4MGdXVztI8llKiIjhw5Amtra/j6+oqOQmSUbiXdgqe9JxRyjjxBpGuSJGHbvpu4tP0WzCSgVreK6NqpEmQymc62ycHqiYpo586d6Ny5M+Ry/jciKm5H7h9B3x19cSL2hOgoREYv4kEy3l98HAvDrqNMPReMmd0Kb3b20WkBBbgnlKjIdu7ciW+++UZ0DCKjE3QjCNOPTsc3Tb9B03JNRcchMlrxqVn4eU80Np25i0HNKmJEfz84WhfveZ8vwxJKVAR37txBZGQkOnToIDoKkVFZcXkF5p+fj7mt56KZRzPRcYiMUpZag+WHb+H3/dfQrLILdo9phYqutiWegyWUqAh27tyJ5s2bw8nJSXQUIqOglbSYfWo2dtzYgSUBS1DLpZboSERGR5IkbNpxDUfCYnClvBkWDPRDMx9XYXlYQomKYOfOnejSpYvoGERGQ61VIzkrGas6r4Kng6foOERG5/i5WOxZHQnbVA3qvFEaP/SvDXMzsdc0yCRJkoQmKIDk5GQ4OjoiKSkJDg4OouOQicvMzISLiwuOHz+O2rVri45DRET0QjEPUrBqyQXY3M2EtpIdBgypA1cXm2LfTlG6GveEEhVSeHg4XF1dUasWDxcSEZF+ylRpsCTkGjK33oXMyRwdxtRHzeouomPlwhJKVEjPDsXreugKIiKiwpIkCUEXHuCHXZEoZWuBcR9UR+vGHqJj5YsllKiQdu7ciTlz5oiOQWSwrjy+ggoOFWBrXvJX4xIZs3MxTzAj6ApiEtLxRafqeLuBB+Ry/d1hwlG2iQohOjoaMTExaNu2regoRAYpPCYcHwR/gOMPjouOQmQ0bt1Lxtg1Z9Fv4TE093FB6PjW6OVXXq8LKMA9oUSFsnPnTrRu3Rq2ttyDQ1RYm69uxqwTszCz+Uy09eIfckSvKyVNiWXLL0J96Qmsa9tj3zh/eDhZi45VYCyhRIWwc+dOvPnmm6JjEBkUSZKw8MJCrLi8An+0+wON3BuJjkRk0DQaLf7ZGoU7++8DZnK80bcK2vh7iY5VaCyhRAWUmpqK8PBw/Pnnn6KjEBkMjVaDWSdmIfROKJZ3Xo6qzlVFRyIyaAeO30X4P1dhnamFZwt39O1TAwqFYZ5dyRJKVEAhISGoWLEiKleuLDoKkcGQy+QoZVUKq7qsQlm7sqLjEBmsu4np+F9wFNJOxqOBpxMGflgXjo6WomO9FpZQogLauXMnunbtKjoGkUGRyWQYXn+46BhEBistS435Ydex+NANdKlTFt9/2xLujlaiYxULllCiApAkCTt37sSyZctERyEiIhOg1UrYcOYuftodhYouNlj3aVPULe8kOlaxYgklKoDz58/jyZMnaNmypegoRERk5PYfjsHa/Tdw2VyDqd1qoUsdd6O8QQpLKFEBbN++HQEBAbC0NOzzb4h06czDMyhrW5bnfhIVUfTNJ1i37CLs4pSoV8cJ8z6pDytzhehYOsMSSlQA27dvx4gRI0THINJbe2/vxVeHvsL05tNZQokKKTE5E0uXXoA8MgUWZSzx1mQ/eHs5iY6lcyyhRK/w4MEDnDlzBl26dBEdhUgvrY1cizmn5+DHVj+itWdr0XGIDIZGK+HvPdfwYPsdaCzl8P+gOpq9oZ/3edcFllCiV9ixYwfeeOMNlC5dWnQUIr0iSRJ+O/sb1kWvw8IOC1Hfrb7oSEQG48i1eEwPugJllhrD25fH292rQG6g430WFUso0Sts374d3bp1Ex2DSK9IkoRvj3yLYw+OIbBzICo5VhIdicgg3IxPw/c7I3Ds+mOMblcFA5tVgKWZ8Z73+TJ6U7l79uwJmUyGffv2iY5ClCMjIwN79+7lrTqJ/kMmk8G3jC9WdVnFAkpUAIlpSswMuoJOcw/Azd4SYRNa4+NWlUy2gAJ6sic0MDAQ6enpomMQ5bF//36UKVMGtWrVEh2FSO/0qNxDdAQivadSaRC49goeHYvDzZo22DqyOaq7O4iOpReEl9C7d+/i66+/xqFDh1ChQgXRcYhyeXYo3hjHZyMiIt3aGXoLZ7bcgIUGqNbBE190rwy5XG8OQgsntIRKkoQPP/wQX3/9Nby8vERGIcpDkiQEBQXxLklERFQoF6PisWXFFdgnqlC6bikMGFQbNjbmomPpHaEldP78+ZAkCZ988onIGET5Onv2LJKTk+Hv7y86CpFQR+8fha25LeqWris6CpFeS0xT4td90dDsjYVjGRv0GtkA5cvZi46lt4SV0OvXr2PGjBk4duxYgZ8zefJkWFhYAAACAgIQEBCgq3hEOXdJevZvjsgU7bixA9OOTsO0ZtNYQoleQKXRYuXR2/g15CoaeDlh4sRGqF7eUXQsndu9ezd2794NAFAqlYV+vkySJKm4QxXE8uXL8cknn8DB4d+Tcx8/fgwHBwf06dMHCxcuzJmenJwMR0dHJCUl5VqeSJcaNmyI0aNHY+DAgaKjEAmx4vIKzD8/Hz/7/4zmHs1FxyHSO5IkYX9kHL7bGQG5TIavutZAm2puomMJUZSuJqyEpqenIyEhIdc0T09PrFmzBh07dkSpUqVyprOEUkm7f/8+vLy8EBsbC1dXV9FxiEqUVtJizqk52H5jO/5s9ydquXJ0CKL/OnvlEf7cfw0nE1Iwtn1VvPeGF8xNbLD55xWlqwk7HG9jYwMbG5s8011dXXMVUCIRgoKC0KRJExZQMknTjk7DydiTWNl5JbwceNEo0fMePEpD4JILsLqVjio1HTB7fBs48qKjIhE+RNPzBO2UJcqDd0kiU9apYieMbjAaLtYuoqMQ6Y3MLDWWr7qE1NOPIdmZwX94HdSva5qH3ouLsMPxhcHD8VSS0tPT4eLigtOnT6NmzZqi4xARkUCSJGHb/pu4vPUWzCSgVpcK6NKpEseP/g+DOhxPpK9CQkJQrlw51KhRQ3QUIiIS6Mr9ZMwIuoLHMSnoU88N/d+vBUsrVqfiwneS6D94lyQiItP2KCULP++JwpZz9zCoaUWMGOgHByue91ncWEKJnqPVahEUFISVK1eKjkKkcwfuHkC6Oh2dKnYSHYVIL2SqNFh2+Bb+CL2GFpVdsXtMK1RwsRUdy2ixhBI958yZM0hLS0PLli1FRyHSqc1XN2PWiVmY3ny66ChEwkmShI07ruHk/hhc8jLD4kEN0aQSL8zTNZZQouds374dnTp14l2SyGhJkoTFFxdj2aVl+L3t72hctrHoSERCHTsbi71/R8I2TYOajd0w8/2aMDdXiI5lElhCiZ6zZcsWTJw4UXQMIp3QaDWYdWIW9t/Zj2WdlqFaqWqiIxEJc+dBClYtuQjbuxmw97HDgCF14VLKWnQsk8ISSvTUjRs3EBERgS5duoiOQqQT3x//HidiT2BVl1UoZ1dOdBwiITJVGizafw2ZW+5C7myODmPqo2Z1HnoXgSWU6KmtW7eiTZs2cHR0FB2FSCf6Vu+LUQ1GwcnKSXQUohInSRK2nb+PH4Oj4GpngQkf10QL37KiY5k0llCip7Zu3YrevXuLjkGkM1Wcq4iOQCTE2TuJmBF0BfeeZODLTtXRo74H5HIOwycaSygRgPj4eBw6dAirVq0SHYWIiIrJzbvJmBd+DcGRcfikVSV86l8JNhasPvqC3wkiAEFBQfD19UX58uVFRyEioteUkqbE0mUXobn8BLY17BAyzh/lnHjRkb5hCSVC9lXxPXr0EB2DqFiE3AnBzaSb+KjOR6KjEJUojUaLf7ZE4U7ofcBMjib9qqB1Ky/RsegFWELJ5KWnp2PPnj34/vvvRUchem3rotZh9qnZmNVylugoRCUq/NhdHFh3FdaZWni1LIs+vatDoZCLjkUvwRJKJm/v3r3w8PBAjRo1REchKjJJkvDHuT+wJnINFnRYgAZuDURHIioRdxPT8cOuSCSejkdTT2cMHFIXDg6WomNRAbCEksl7diheJuOVkmSY1Fo1ZhybgSP3jyCwcyB8nHxERyLSudQsNeaHXcOSQzfRtU45zPm2Fco4WImORYXAEkomTa1WY/v27di6davoKERF9svpX3Ax/iJWdV6FMrZlRMch0imNVsLG03fx054oeLvYYv2nzVCnPMd3NkQsoWTSDh8+DIVCgSZNmoiOQlRkg2sNxqf1PoWDhYPoKEQ6te9wDDbsvY5L1lpM614LnWu78yiWAWMJJZO2ZcsWdO/eHQqFQnQUoiIrbVNadAQinYq6mYj1yy7BLk6JurWdMPfT+rAy589tQ8cSSiZLkiRs3boV8+bNEx2FiIjykZCciWVLLkAelQLLMpboMbkhKnrx0LuxYAklk3Xx4kXExcWhXbt2oqMQEdFz1Bot/g65jgdb70CyVMD/gxpo+kY50bGomLGEksnasmULOnXqBGtr3kWDDEPwzWAcfXAU05pNEx2FSGcOXY3HzB1XkKVUY1RHT/R4szLkHO/TKLGEksnasmULPv/8c9ExiApk5ZWV+P3s75jtP1t0FCKduPEoFd/vjMDxmwn4rF0VDGxaERZmLJ/GjCWUTNLt27dx8eJFdO3aVXQUopfSSlrMPT0XW65tweKOi1GndB3RkYiKVWJKFn4Lu47Vx2+jd0NPhI1vDRc7DjZvClhCySRt27YNrVq1grOzs+goRC+k0qgw5cgUnIs7h5VdVqKCQwXRkYiKjVKlQeDaK3h8NA53a9hg28gWqOZuLzoWlSCWUDJJz+6SRKSvJEnC4ouLcePJDazqsgqu1q6iIxEVmx37b+Ls1puw0ADVArwwoZsP5HIeejc1MkmSJNEhXiU5ORmOjo5ISkqCgwMHY6bX8/jxY7i7u+P69evw8vISHYcoX5IkIVOTCa2kha25reg4RMXifGQ8tgZegX2iCrb1SmHAoNqwtjYXHYuKQVG6GveEksnZunUrGjRowAJKek0mk8HajCM3kHFITFNi7r5oZO6LhVsZG/Qd7Yuy7naiY5FgLKFkcjZu3Ih33nlHdAwiIqOnVGux8tht/LovGo0qlsLEyY1RpSyPaFI2llAyKUlJSdi3bx9+/fVX0VGIctFKWshlPCeOjIMkSQiJiMP3OyNgppDh9/d80aoqby9LubGEkkkJCgpC9erVUblyZdFRiHJsu74N265vw6IOiyCTyUTHIXotZ67EYcHeaziZnIaxHaqiXyNPmHGwecoHSyiZFB6KJ30iSRKWXFqCJReX4Jc2v7CAkkG7/ygNgUsuwPpWOqpUt8eP41vDkRcd0UuwhJLJSEtLQ3BwMGbMmCE6ChE0Wg1+PPkj9tzeg2WdlqF6qeqiIxEVSUaWGstXXkLa6ceAgxlaj6iDenXcRMciA8ASSiYjODgYXl5eqFmzpugoZOKyNFmYfHAyohOjsbLzSpS3Ly86ElGhSZKEbaE3cWnzLZgBqPNWRXQOqMQ9+lRgLKFkMjZu3Ii3336bPyBJuM1XNyM2LRaBnQPhbMW7dpHhuXQvCTOCruDR/VT0b+CG/v1rw8JCIToWGRgOVk8mISsrC6VLl0ZoaCj8/PxExyETp5W0UGqUsDKzEh2FqFDikjMxe08Utp67jw+ae2NEGx/YW/G8T+Jg9UQvtHfvXri4uMDX11d0FCLIZXIWUDIomSoNlhy6iT9Dr8G/WmnsHesPLxcb0bHIwLGEkkngoXgiosKTJAnrg67ifMhdXKhohqWDG+GNSi6iY5GRYAklo6dSqbBt2zZs375ddBQyQVmaLFgqLEXHICq0o2ceYO+aKNilaVD1DTdMfa8mzM153icVH5ZQMnrh4eGwtLREkyZNREchE7M+ej1WX1mNDd03wEzOH7dkGG7fT8HqJRdgey8Tjj52GDCkLkqVshYdi4wQfyqS0du4cSN69uwJuZx37KCSIUkS5p+fj1URq/Bb299YQMkgZCg1WBh2HRlbYqBwskDHsfVRoxoPvZPu8CcjGTWNRoPNmzfj77//Fh2FTIRaq8bMYzNx6N4hBHYKRGVn3iKW9JtWK2Hb+fv4X3Ak3Bys8OXQWmhW1110LDIBLKFk1I4cOQK1Wo1WrVqJjkImIEOdgS8OfIG7KXexqssquNvyFznptzN3EjF9+xXEJmViYufq6F6vHORyXsBJJYMllIzaxo0b0aNHD5iZ8Z866d6Re0eQnJWM5Z2Ww9HSUXQcohe6eTcZ88KuYXd0HD5pVQmftKoEGwv+nKSSxX9xZLQkScKmTZuwYMEC0VHIRLSr0A6tPVtDIecVxKSfktOUWLbsAjSXk+BQ1RYh4/xR1pEXHZEYLKFktE6dOoXk5GS0a9dOdBQyISygpI80Gi3WbI7C3bD7gJkcTd6ritYtPUXHIhPHEkpGa926dejevTssLCxERyEiEib8+D0cWBcN6wwJFVqWRe/e1aFQcLQQEo8llIySJElYt24d/vzzT9FRyEglZibC2cpZdAyiF4pJSMcPuyLx8Fw82niWwoAP68DBgTdOIP3BP4XIKJ04cQJJSUno0KGD6ChkhFZHrMY7295BhjpDdBSiPFIyVfhfcCQ6/BIOGwsF/pzSCiPGNGQBJb3DPaFklNatW4eePXvyUDwVK0mSMPfMXGy6ugl/tPsD1ma8oIP0h0YrYf2pGMzeE41KpW2xYWgz1PbgKA2kv1hCyehotVqsX7+eV8VTsVJpVfj28Lc4E3cGgZ0D4e3oLToSUY59h+9gy+4bOG+nxcwetRBQyx0yGcf7JP3GEkpG5/jx40hNTeVV8VRs0lXp+DzsczzOfIyVnVeitE1p0ZGIAACRNxKxftkl2D9SomYtJ8weWg9W5vzVToaB/1LJ6Pzzzz88FE/FKjoxGnKZHMsClsHOwk50HCIkJGdi2ZILkEelwMrdCj0mN0RFLx56J8MikyRJEh3iVZKTk+Ho6IikpCQ4ODiIjkN6TKvVwtPTE8uWLUPHjh1FxyEjIkkSD2+ScGqNFn+H3cSDTbcgWSnQuk9VNGlcTnQsoiJ1Ne4JJaNy5MgRZGVloU2bNqKjkJFhASXRDl59hBlBV6BSazG6SwW81akS5BzvkwwYSygZlXXr1uHtt9+Gubm56ChERMXiWlwqvt8ZgZO3EvBZuyoY2LQiLMxYPsnwsYSS0dBoNNiwYQMCAwNFRyEDdif5DrwcvETHIMLjpEz8Fn4da07cQZ9Gnpj9bhuUsuW57mQ8WELJaBw+fBgqlQqtW7cWHYUMkCRJWH55ORZfXIydb++EoyUv8iAxlCoNVqy5jMRjj/Comg2CRrVAlTL2omMRFbtCl9B169ahd+/eePLkCZycnHQQiaho1q1bh3feeQdmZvzbigpHK2nx08mfsOvmLizuuJgFlISQJAk799/C2W03YakBqnbywoRulXk+MhmtQv+2fnbF05tvvonk5GTUqFEDfn5+8PPz47iMJMyzQ/F///236ChkYJQaJSYfmozIhEis7LISnvaeoiORCTofGY+tgVfgkKiCW11nDBhcB9bWPLedjFuhS2inTp0AAIcOHYJWq0VERAROnz6N7du3s4SSMAcOHIAkSfD39xcdhQxIijIFY0LHIF2VjsDOgShlVUp0JDIxCWlK/LI3Gin7H8DTzQ59RvmibFmORUumocjHLdu1a4c9e/agVq1aqFatGtq2bVucuYgKZd26dejVqxcUCoXoKGRAEjMT4WbjhilNpsDG3EZ0HDIhSrUWgUdv4deQq3jDuxS+nNwEVdx53ieZliKXUEmScn7hm5mZoX///ggLCyuuXEQFplarsXHjRqxfv150FDIwXg5emNVylugYZEIkScK+iDh8vzMCFgo5/nzfFy2r8DawZJqKXELNzMzw+PFjuLi4AMj+j0UkQnh4OBQKBVq0aCE6ChHRC52+HIdFu6/iRFo6Pu9QFX0becKMg82TCStyCZ06dSratm2LoUOHQqVSsYSSMP/88w8PxROR3roXl4aVSy/A6lY6qla1x48TWsPBihcdERW5hDZr1gzbt2/HmjVrkJmZiTVr1hRnLqICUSqV2LhxI7Zu3So6Cum5c3HnUK90PQ53QyUmI0uNZasuIf3UY8DBDG1H1EHdOm6iYxHpjdcaUNHLywtffvllcWUhKrQ9e/bAzs4OzZo1Ex2F9JQkSVhwYQECrwRiU/dNcLd1Fx2JjJwkSdgadguXN92EOYA6PbzRqaM3/wAi+o8Cl9DQ0FDs2rULaWlpqF69Ot577z24uLhg+PDh+PPPP3WZkeiF1qxZg379+kEu53lVlJdGq8F3x79D+N1wLO+0nAWUdO7SvSRMD7qCB7FpGOLrhvferwULC54qRJQfmVSAkzkHDx6MlStXwtvbG87OzoiKioKVlRU2bNiQM2i9LiUnJ8PR0RFJSUk5g+UTpaWlwc3NDUeOHEG9evVExyE9k6nOxJcHvsSt5Fv4q/1fKGtXVnQkMmJxyZn4aXcUtl+4jw+be2N4m8qws+Td28h0FKWrvXL30V9//YWwsDAcP34c165dw8mTJxEbG4tRo0ahe/fuyMrKKlLYadOmwcfHB46OjnB1dUVAQADOnTtXpHWRadq2bRsqVqyIunXrio5CeiYpKwmf7P0ECZkJCOwcyAJKOpOp0uD3/VfRZnYY0lUa7B3rjy86VWcBJSqAV/4vWbhwIRYsWICGDRvmTLOxscGUKVNgY2ODCRMmFGnDffv2xejRo+Hs7AylUonffvsNAQEBuH//Pq9ypgJ5diie51lRfmq71sboBqNhZWYlOgoZIUmSsG77VVwJuYsL3uZY/mFjNKrIO24RFcYrD8fb29vj8ePHsLCwyDNPkiQcPHgQrVq1eq0QWVlZmD9/PsaOHYu4uDiULp174F4ejqf/SkhIgLu7OyIiIuDj4yM6DhGZkCNnHmDfmijYpWng+oYb+r1XE+bm3HlCpq0oXe2Ve0LlcjkyMzPzLaEqlQq3b98ufNKnduzYgffffx9JSUmQyWQYO3ZsngJKlJ+NGzfC19eXBZSISsyt+8lYveQi7O5lwtHHHgM/qgNnZ2vRsYgM1ivPCW3UqBH+/vvvfOd9/vnnGDx4cJE33rVrVzx58gSPHz/Gzz//jKZNmxZ5XWRa/v77b/Tr1090DCIyAelKNX7ZE4VV352EWYYWAZ83wGcTGrOAEr2mV+4J/fLLL9GzZ09kZmaid+/ecHZ2xvnz5zFr1ixERkbC0tLytUOUKlUKn332GZydnVG1atUXXuk8efLknD2yAQEBCAgIeO1tk+G5d+8eDh069MI/jsi07L+zHy3Lt4S5nHegoeKl1UrYcu4efgyOQlknK3w5vA6a1OJg80TP7N69G7t37waQffOYwirQEE3Lli3DmDFjkJqaCiD7XNBWrVph9erVqF69OlJSUgq94f9Sq9VwcHBAYGAgevXqlWsezwml582ZMwc7d+7Evn37REchgSRJwryz87A+ej1Wdl4Jb0dv0ZHIiJy+nYDp26/gUUoWvuxcHd3rleNFkEQvoZNzQjds2IA+ffqgV69eOHToEJKTk1G9evWcvZWTJ08uUthff/0Vffv2RZkyZfDo0SN89dVXsLCwQPPmzYu0PjIdf//9N4YPHy46Bgmk0qow9chUnIo9hcDOgSygVGyuxyRhXshV7L0Rj2H+PvioZSVYc7B5Ip145Z7QFi1a4Pz582jfvj3eeecddO/evVj2Rr755ps4efIkUlNT4eDggEaNGuGbb77JNRTUM9wTSs9ER0ejTp06ePjwIZycnETHIQHSVekYFz4Oj9If4c/2f8LNhodH6fUlpSmxdOkFSFeSkFLZFkM/ro8yDhzei6igdLIn9NChQ3jw4AE2b96M5cuXY+jQoWjZsiXeeecd9OjRA66urkUKGxQUVKTnkWlbs2YNunTpwgJqopKykjB071DYmNtgWadlsLewFx2JDJxGo8Xfm6NwN/Q+5OZyNH2vKlq19BQdi8gkFOic0OfFx8dj27Zt2LRpE0JDQ9G4cWO888476N27N9zcdLNHgntCCcg+B7BGjRqYPn06evfuLToOCaDSqLAyYiX61+gPC0XeYeOICmP/sbs4tO4qbDMllPd3R69e1aFQvHLQGCLKR1G6WqFL6H83uH37dmzevBlNmzbFuHHjirqqV26HJZTOnDkDf39/PHz4EDY2NqLjEJGBuvM4Hd/vjMD9S48RUK4UBnxQBw4Orz/SC5Ep02kJjY2NxcmTJ5Geno6yZcuiefPmJXZ7TZZQAoDx48cjLi4OgYGBoqMQkQFKyVTh99BrWHHkFt6q54FxAVXhZs/zPomKg07OCX2mXLns4SmedVZnZ2d89dVXGDt2LIetIJ3TarVYu3YtFi1aJDoKERkYjVbCPydj8POeKFQtY4+Nw5qhVjlH0bGITF6BS+jp06dRo0YNWFhYIDY2FkFBQfjuu++QmZlZ5GGaiArqwIEDyMrKQvv27UVHoRKy+epmtPVqC0dLlgUqur2H7mDHrus46wTMersOOtQswx0nRHritc4JvXbtGtq3b49bt24VY6S8eDiePvzwQ9jZ2WHevHmio5COaSUtZp+ajR03dmBRx0Wo6lxVdCQyQBE3ErF+2SU4xCthWdMJgz+tB2uLAu93IaJC0unh+GeSkpJw9+5dlC1bFnfu3EFycnKhgxIVRnp6OjZs2ID9+/eLjkI6ptQo8fXhr3Ep/hJWdV4FTwcOlUOF8/hJJpYuvQCz6BTYlLXC25MbwcuTOy+I9FGhS+iNGzfQrFkzKJVKVK9eHXPnztVBLKJ/bdu2DR4eHvDz8xMdhXQoVZmKMaFjkKJKwcrOK+Fi7SI6EhkQtUaL1QduInbDLcBKgWYf1sAbjcuJjkVEL1HoEtqgQQMkJiZix44d+P7772FnZ6eLXEQ5AgMDMWDAAJ7HZcTSVGn4YPcHcLFywbKAZbAx5xBcVHDh0Y8wM+gKNJKEsW9VRNd23pBzvE8ivVfgc0KHDRuGWbNm5bpTzcOHD9G0aVPcuHFDV/kA8JxQU/bw4UN4enri2rVr8PLyEh2HdESSJATdCEIn704wl5uLjkMG4lpcCmbuiMDZO08wpn0V9G9SAeYsn0RC6PSc0Dt37sDb2xsffvghOnXqhAoVKuDkyZNISUkpcmCiV1mzZg2aN2/OAmrkZDIZuvl0Ex2DDMSjhAz8fuAa1p66i36NvTC3T3042fAOWkSGpsAldMeOHdiwYQNmzZqFX375JefQ6MSJE3UWjigwMBCjRo0SHYOI9ECWUoMVay7jyfFHeFLFGjtGt0BlN3vRsYioiIo0RNPDhw8RExOD0qVLo0KFCrrIlQsPx5umy5cvo1GjRoiNjeX3nciESZKEoJCbOLf9Fqw0QLUAT3R7szLPEyfSIyUyRBMAlClTBmXKlCnKU4kKbOXKlXjrrbdYQI2IJEkIvBKINp5t4OXAUyzo1c5GPMK2wCtweKJG2fql0H9gbVhZ87xhImPAkXtJL2k0GqxatYq36TQiGq0Gs07Mwv47+9G0XFPRcUjPxadmYc7eaCSGxaKKmx36ftYQZd1tRcciomLEEkp6KSwsDGq1Gh06dBAdhYpBpjoTkw5OwvWk61jVZRXK2XH8RspfllqD5Ydv4ff919DExwWTvmqKSm4cCpDIGLGEkl4KDAxEv379YGbGf6KGLikrCaP3j4ZG0iCwUyCcrJxERyI9JEkSdl+Oxfc7I2FjocCCAX5oVtlVdCwi0iH+hie9k5aWho0bN+LgwYOio9BrUmlU+GD3B/Cw88CPrX6EtZm16Eikh05ejMOSXdE4pczAuI7V0LuhJxRyXnREZOxYQknvbNmyBRUrVkT9+vVFR6HXZK4wx5eNvoRfGT+YyfnjhnK79zAVK5ZegO3tDFStbI+fxr8BeytedERkKvhbgfROYGAgBg4cyOFXjMQbZd8QHYH0THqmGstWXkT6mQQoHMzQdmRd1KldWnQsIiphLKGkVx48eID9+/djyZIloqMQUTGTJAmbw28hYuNNmMuAej280bGjN//gJDJRLKGkV/7++2/4+/ujfPnyoqMQUTE6H/MEM4Ku4EF8Oj7xK4P33qsJcwuF6FhEJBBLKOkNSZKwbNkyfPnll6KjUCFJkoT55+ejkXsjNHJvJDoO6ZHYpEz8GByJnZce4OOWlTD0w8awteSvHiJiCSU9curUKdy5cwfvvPOO6ChUCGqtGtOPTsexB8fQqWIn0XFIT6RnqbHo4E0sOHAd7WuUQci41vBw4ugIRPQvllDSG0uXLkXfvn1hY2MjOgoVULoqHRMOTMCDtAdY2Xklytjydr6mTquVsG57NKJD7uJSJUusHNIYfhVKiY5FRHqIJZT0QkZGBtasWYNdu3aJjkIFlJiZiJEhI2GhsMDyTsvhYOEgOhIJdvD0fexfEwX7dC0qNS2DSX1rwNyc530SUf5YQkkvbNmyBWXLlkWTJk1ER6EC0EpafLr3U5S3L49ZLWfBUmEpOhIJdONeMv5echH29zNRqrI9BgypCydnK9GxiEjPsYSSXli6dCk++OADDtViIOQyOb5v8T28Hb2hkHNPl6lKy1Ljr7BrSN5+F7aOlugy3hdVKjuLjkVEBoIllIS7c+cOwsPDERgYKDoKFUJl58qiI5AgWq2EjWfu4qfdUfAqZYOJo+qjYTXe552ICocllIRbsWIFAgICULZsWdFRiOgVjt94jBk7riAxTYVvutVE1zpleQSDiIqEJZSE0mq1WLZsGWbPni06ChG9xLU7SZi3Lxohtx5jeJvKGNLCG1a86IiIXgNLKAl14MABpKSk4M033xQdhfKhlbT45fQvqOpcFd18uomOQwI8Sc3CkqUXIYtIglslW4ROaA03e150RESvjyWUhFq6dCn69+8PCwsL0VHoP1QaFb4+/DUuPLqAd6u+KzoOlTC1RovVGyNxP/wB5OZyNHu/Klq28BQdi4iMCEsoCZOcnIwNGzbg2LFjoqPQf6Sp0jAmdAySspKwsstKuFrzohNTEnL0Lg6tvwq7LAk+rcrhnV7VoFDIRcciIiPDEkrCrFu3DjVq1EDdunVFR6HnxGfEY/i+4XC2csayTstga24rOhKVkJvxafh+ZwTuRCSgR8VS6P9BHdjbcwxYItINllASZunSpfjwww9Fx6D/mBA+AT5OPpjebDrMFeai41AJSMpQ4beQq1h57Dbe8SuPWV/5w9WO5ZOIdEsmSZIkOsSrJCcnw9HREUlJSXBw4K0BjUFkZCTq16+P+/fvo1Qp3ldan8Slx8HV2hVyGQ+/Gju1Ros1J+5gzt5o1CrniK/frIHq7vwZS0SFV5Suxj2hJMSyZcvw1ltvsYDqITcbN9ERqAQEH7iNvTtu4KwrMPvdemhb3Y3jfRJRiWIJpRKnVquxcuVKLF26VHQUIpNz6VoCNi6/BMd4FarUdsLMj+rB2oq/Coio5PEnD5W4HTt2wNzcHB06dBAdxaQ9OxOHe79MQ/yTDCxdcgHmV1NhV9YKvb6qj/KePPROROKwhFKJW7hwIYYMGQKFgndbEUWj1eCHEz/A1doVn9b7VHQc0iGVRouVB28ibv0tyKwUaPFRDTRqWE50LCIillAqWXfu3MHevXuxYMEC0VFMVpYmC5MOTsLVxKv4q8NfouOQjkiShNCoOMzcEQG5TIbPe/mgs38FyOTc801E+oEllErU0qVL0alTJ5QvX150FJOUrEzG6P2jodKqENg5EM5WzqIjkQ5ExaZg5o4ruHQvCWM7VMV7jb1gxsHmiUjPsIRSiVGr1ViyZAnmz58vOopJik2LxbB9w+Bh54Gf/H+CtZm16EhUzGIfp+O3sGvYcPYe+jepgN/7+cLRhmO9EpF+YgmlEhMcHAxJktCpUyfRUUzSDyd+QL3S9fB1k69hJud/fWOSqVRj+d+XkXwiHlk+Ntj1WUtUKm0nOhYR0UvxNxGVmGcXJJmZ8Z+dCN+1+A42Zja8Gt6ISJKELXtv4tKOW7DSADW7VMCErj78HhORQWAboBJx7949BAcH4/fffxcdxWTxHvDG5dTlOAStjIBTkhrlGrig/8DasOR4n0RkQPgTi0rE0qVL0aFDB3h5eYmOQmTQ4lIy8fPuaDw6HIs6pe3Rd2xDuJfhHxhEZHhYQknnNBoNFi9ejHnz5omOYhIkSYJaq4a5ghekGJNMlQZLD9/En6HX0bKKK779qjkquLJ8EpHhYgklndu7dy9UKhW6du0qOorRU2vVmHlsJuQyOb5p+o3oOFQMJEnCzouxmLUrAo7W5lg8qCGaVHIRHYuI6LWxhJLOLVy4EB9++CEvSNKxDHUGvgj/AndT72J+ew6DZQyOX4jFip1XcVKThS8CquEd3/KQc7B5IjISbAWkUw8ePEBQUBB+/vln0VGM2pPMJxixfwTMZGZY3mk5HC0dRUei1xATm4rApRdgdycD1Xzs8NPYprC15I9rIjIu/KlGOrVs2TK0bdsW3t7eoqMYrfup9/Hp3k9R2akyfmj1AywVlqIjURGlZaixdOVFZJ5NgLmjOdqProtaNUuLjkVEpBMsoaQzWq0Wixcvxk8//SQ6ilFbemkp3ij7BiY1ngSFXCE6DhWBJEnYdOAWIjbcgIVMhgZvV0L79hU53icRGTWWUNKZkJAQpKeno3v37qKjGLWJjSdCIVOwsBioM3cSMX37FcQ9ycTwN8qgT+8aMLfgHxNEZPxYQkln5s+fjw8//BDm5hwqSJd4C07DdO9JBv63KxJ7rzzEp/6V8EmrSrCx4PeSiEwHf+KRTsTExGDHjh345ZdfREch0iupmSosPHADiw7eROfa7tg/3h9lHa1FxyIiKnEsoaQTCxcuRKdOnVChQgXRUYyGVtIiVZUKBwsH0VGoCLRaCWu2ReNGyF1E+VhizSdNUN/TSXQsIiJhWEKp2CmVSixatAiBgYGioxgNlUaFb498izRVGn5t+6voOFRI4SfuIXRdNBzSJfg0K4NJfWvAzIznfRKRaWMJpWK3adMm2Nvbo3379qKjGIU0VRo+D/sciZmJ+LP9n6LjUCFci0nC30svwvFBFkpXcUD/IXXg6GQlOhYRkV5gCaVi9+eff2LYsGGQy+Wioxi8+Ix4DN83HI6WjlgasBR2FnaiI1EBpGSq8EfodSTvvAtnB0u8Od4XPpWdRcciItIrLKFUrC5evIhTp05h69atoqMYvDvJd/Dp3k9Rt3RdzGw+E+YKjjKg7zRaCetPxWD2nij4lLbDpDENUN+H93knIsoPSygVqz///BP9+vWDszP3+ryuPbf3oJ1XO3ze8HPIZdyrrO+OXI/HjKAIpGWpMbNHbQTUcufYrUREL8ESSsUmKSkJK1euRHh4uOgoRmFI7SEsMQYg6tYT/LY7CuH3n2Bk28oY3LwiLHnRERHRK7GEUrFZuXIlateuDT8/P9FRjAILqH5LSMnC0qUXII9Mhoe3LUIntIarnaXoWEREBoMllIqFJEn4888/MXHiRNFRiHRKpdZg5cZIPDwQC4WFAs0GVkOzpuVFxyIiMjgsoVQswsLCEBcXh969e4uOYnC0khYP0x6irF1Z0VHoFfYcjcHRdVdhpwSqtvZAj7erQqHg+bpEREXBEkrF4s8//8SQIUNgZcUxEAtDqVFi8qHJeJzxGEsDlvIQvJ66FpeK73dG4NbVBPSu5Ir+g+vAzt5CdCwiIoPGEkqv7d69e9i2bRsiIyNFRzEoKcoUfBb6GTLVmfi93e8soHooMU2JX0OuYs2JO+jbyBOzJ7dBKVuWTyKi4iDsONLEiRNRp04dODg4oGzZsujXrx9iYmJExaHXsGjRInTo0AHe3t6ioxiMh2kPMSh4EKwUVljccTFKWZUSHYmeo9JosfTQTbSeHYZbj9MQNKoFpr1VmwWUiKgYCSuhMpkMy5cvR3x8PCIiIiCTydCtWzdRcaiIlEolFixYgOHDh4uOYjBuPLmBAbsGoJZLLfza9lfYmNuIjkRPSZKEoLBbmDwxDGuO38bcvvWx/IPGqFLGXnQ0IiKjI5MkSRIdAgDOnTuHBg0aICEhIc9A58nJyXB0dERSUhIcHBwEJaT8rFq1CtOnT0dkZCRv01lAu27uwtXEqxjVYBQPweuR89GPsXnFZTgnqGBTxxn9P6wDayvepYqIqCCK0tX05pzQPXv2oEKFCrzTjgGRJAm//vorRo8ezQJaCJ29O6Ozd2fRMeiphwnpWLrkAqyup8GhnDXe/boBynlwzycRka7pRQndt28fpk2bho0bN4qOQoVw5MgRXL16FYMHDxYdhajQstQaLD94E4833IaZtRlafVwTfn4cJouIqKQIL6FBQUHo378/Vq1ahU6dOr102cmTJ8PCIvvCgICAAAQEBJRERHqBX3/9FUOGDIGdnZ3oKEQFJkkSdl9+iFm7ImBtrsC496qifdPykMl5agQRUWHs3r0bu3fvBpB9jUhhCT0ndPXq1Rg+fDjWrVv30kLJc0L1z507d1ClShVERkbyqvgXUGvVuJp4FTVcaoiOQk9dvp+EGUFXcPVhKsZ1rIY+jTyhYPkkInptBnVO6O+//44pU6YgKCgILVu2FBWDiuiPP/5A165dWUBfIFOdiS8OfIGH6Q+xpusayGU8Z1ak+/FpmBdyDVsu3cegZhWxcGBDOPCiIyIioYTtCZXJZDAzM4OlpWWu6bt27cpTSrknVL+kpaWhfPny2Lp1K1q1aiU6jt5JykrCyJCRkMvkmNd2HhwtHUVHMlkZWWosXX0JaaceI6mSDT79oB68XDgkFhFRcTOoPaF6MjIUFcHKlStRsWJF7sHOx/3U+xi6bygqOVbCDy1/gJUZb2MqgiRJ2LT7Bi7vvA0bCajTtQI6dfHhkFhERHpE+IVJZFi0Wi1+/fVXTJw4kb/Q/yM6MRrD9g5DG682mNR4EhRyhehIJunExYfYsToCzkkaePq64v0BtWBhxR91RET6hj+ZqVD27t2LhIQE9O3bV3QUvZOuSke/Gv0wpPYQFnQBYpMy8ePuSMQej8Mbrg547/M6KO1mKzoWERG9AEsoFcrcuXMxbNiwPOfyElDfrT7qu9UXHcPkZCg1WHjgBhYcuI72Ncrgx69boLwzz/skItJ3LKFUYJGRkQgNDcWyZctERyGCVith2/n7+F9wJNwcrLByyBvwq8A7rhERGQqWUCqwefPmoU+fPnB3dxcdhUzc4bOxWB0UjbNmKnzZqTq61ysHOcf7JCIyKCyhVCCJiYlYsWIFDh48KDqKcCqtCscfHEcLjxaio5icmw+SsXLJRTjezUQNH3vMHt0MNpb8MUZEZIj405sKZMGCBfDz84Ovr6/oKEKlq9LxefjnSMhIwBvub8BcwQHPS0JKugpLAi9CdT4R1k7m6DSmPqpVdxEdi4iIXgNLKL1SVlYWfv31VyxcuFB0FKEeZzzGiJARsDO3w9KApSygJUCrlfDPwZu4vuEmLGVy+L3jgzbtKnD0ASIiI8ASSq+0evVqODk5oWvXrqKjCBOTHIOh+4ailkstzGwxExYKC9GRjN7xG48xY8cVPElVYWSLsuj1djWYmXPsVSIiY8ESSi+l1Woxe/ZsjB8/HnK5ad7/PDIhEp/u/RRdvLtgQqMJvA+8jt15nI5ZuyJwIPoRhrepjCEtvGHF8klEZHRYQumldu7cicTERPTv3190FGHsLewxtN5Q9KveT3QUo5acrsQf4dex4sgtvFXPA6ETWsPNnrc9JSIyViyh9FI//fQTRo8ebdKD03vYebCA6pBGK2HV5kjEhN7HLR8rbBjaDLU9HEXHIiIiHWMJpRc6ceIEzpw5gy1btoiOQkZq37G7OLj+KhwzJFRpURYTe1eDmRkPvRMRmQKWUHqhn376CR999BGcnXkXGipekbcSsXbZJTg/VMK9mgPe+7AOHB156J2IyJSwhFK+rl+/jm3btiE6Olp0lBKj1CgRdCMIPSv35BBAOpKUrsKvIVeRtPceyjlYo/sXvvCuxD9yiIhMEUso5WvOnDno1asXKlSoIDpKiUhRpmBM6BikqdIQUDEAtua2oiMZFZVGi7+P38Ev+6JRx8MRk8Y3Qk0vJ9GxiIhIIJZQyiM+Ph7Lli3DkSNHREcpEY/SH2HYvmFwtXHF0rZLYWNuIzqSUQmNisN3OyKglSTM6V0Pbaq5cU8zERGxhFJef/zxB5o3b4769euLjqJzN5NuYujeoWjo3hBTm02FuZx3QSoul64n4PddUTj6OAVj2ldB/yYVYK7gGKtERJSNJZRySU9Px++//47Vq1eLjqJzUQlR+GjPR+hVtRdGNxjNvXPF5FFSJpYsPQ+L6FT4VLTFDxNaw8mGd5giIqLcWEIplxUrVsDDwwMdOnQQHUXnPOw88GXjL/FmpTdFRzEKWSoNlq+7gsdH4mBuqUDzQTXwRpNyomMREZGeYgmlHGq1GrNnz8a0adNMYq+gnYUdC2gxkCQJu47E4MT6a7BXATXblkf3nlUhlxv/vyEiIio6llDKsXbtWgBA3759BSchQxHxIBkzd1zB3ZgUDKjigvcG1YaNHQ+9ExHRq7GEEgBAq9Vi1qxZ+PLLL2Fmxn8W9HLxqVn4eU80Np25i4FNK+DP9/3gaM2LuoiIqOB4qSoBALZs2YKkpCQMGjRIdJRil6nOxPzz86HUKEVHMXhZag3+Cr+ONj+FISEtC7vHtMJXXWuygBIRUaFxlxdBkiR8//33GD9+PCwtLUXHKVZJWUkYtX8UJEnC+zXeh4WCh4qLQpIkbN1/E8d23sLl8mZYOLAhmvq4iI5FREQGjCWUsGfPHty+fRsff/yx6CjFKjYtFkP3DoWXgxd+bPUjrMx4b/KiOHXlEbauugKXRDWq1iuFqYPrwMqKPzqIiOj18DcJ4bvvvsOYMWNga2s8t6q8mngVQ/cNRavyrfDVG1/BTM5/6oV1Lz4Ny5ZcgO3NdLiWt0GfUX5wL2snOhYRERkJ/mY2cQcPHsT58+exbds20VGKzY2kGxgUPAgDag7A0LpDTWK4qeKUqdJgUfh1PNl8B1Y25mg7tDbq1S8jOhYRERkZllAT991332HkyJFwcnISHaXYeNlnH35v4dFCdBSDIkkStp2/jx+Do+BiZ4HxH9RES9+ykHG8TyIi0gGWUBN26tQpHDx4ECtXrhQdpViZyc1YQAvpzJ1EzAi6gvtPMvBFQHX0bODBweaJiEinWEJN2KxZs/DJJ5+gdOnSoqOQIHcepuKXvdEIjorDJ60q4VP/SrCx4I8FIiLSPf62MVFXrlzBzp07ce3aNdFRSIC0TBUWrbwI5dlEOHnbYP94f5R1tBYdi4iITAhLqImaNWsWBg4cCA8PD9FRiixdlY45p+dgeP3hKGVVSnQcg6DVSli76xqu7r4DW0mGBt290T7AmxdvERFRiWMJNUHXr1/H+vXrceXKFdFRiiwhMwEjQ0bCyswK5nLeracgDp59gN1rouCaooVPo9Lo835NmFvyRwAREYnB30AmaMaMGejXrx8qVaokOkqRxKTEYNi+Yaheqjq+b/E974L0Cncep+OH4AjEnXuMNq6O6PdFHZRytREdi4iITBxLqImJjo7G2rVrcfnyZdFRiiTicQSG7RuGTt6d8EWjLyCXyUVH0lspmSr8HnoNK47cQvd65TB1Siu4OfCuUUREpB9YQk3MjBkz0L9/f/j4+IiOUmixabH4cPeH+Ljux/ig1gc8j/EFNFoJ/5yMwZy9UfApbYcNQ5uhtoej6FhERES5sISakMjISKxfvx6RkZGioxSJu607FgcsRi2XWqKj6K39p+5h3bariLCXMLNHHQTUKsOyTkREeokl1IRMnz4dAwcORMWKFUVHKTIW0PxF303CqiUX4fIgC/Uq22PuaF9YcbxPIiLSY/wtZSKuXLmCzZs3G+xeUMpfYooSiwIvQH4pCU7OFuj6eQNUrsrhqoiISP+xhJqIadOmYfDgwahQoYLoKFQMVBotVoffRMymW7BWyNGkTxU09/fkoXciIjIYLKEm4NKlS9i2bRuio6NFRymQVGUqvj3yLUY2GAlvR2/RcfROaFQcvtsRAa1Wi9HtPfFmVx+YmStExyIiIioUllATMG3aNAwZMgSenp6io7xSfEY8hu0bhlJWpVDGpozoOHrl6sMUzNgRgfMxTzC2fRW836QCzBUcooqIiAwTS6iRO3/+PIKCggziHvG3km5h6L6haODWANObTYe5gndCAoD45Cz8uv8q1p2KwXtveGFe3/pwsuEA/UREZNhYQo3ctGnT8PHHH+v9PeIvPLqAkSEj0bNKT4zxHcNzGwEo1Vos23AFjw4+REpVG+z8rCV8StuJjkVERFQsWEKN2NmzZxEcHKz3e0GTspIwdN9QDK83HP1r9hcdRzhJkhB06A6Ob7oOlyyglr8HuveqCgUPvRMRkRFhCTViU6ZMwdChQ1GuXDnRUV7K0dIRG7ptQDk7/c5ZEs5de4z1Ky7D7ZEKFWo6od/g2rBzsBQdi4iIqNixhBqpgwcP4sCBA1i2bJnoKAVi6gU0LiUTc/ZEI+3AQ/g42eDtSfXgWYG32iQiIuPFEmqEJEnCxIkTMX78eJQuXVp0HHqJTJUGSw7dxPyw62hZxRVfTGoKb3ee90lERMaPJdQIBQUF4erVqwgODhYdJQ9JknjREbLfhx0XH+CHXZFwtDbHkkEN8UYlF9GxiIiISgxLqJHRaDSYPHkypkyZAnt7e9FxcknKSsK48HEYXm84fMv4io4jzMnIR5gfFIWLWZn4IqAa3vEtD7mcxZyIiEwLS6iRWb16NVJTU/HJJ5+IjpJLbFoshu0bhvL25VHDpYboOELcfZyGJUsvwP5GOupVtMNv45vA1pL/BYmIyDTxN6ARycrKwjfffIOZM2fC0lJ/rqi+lngNQ/cNRQuPFvi6ydcwk5vWP7u0LBUWrbmC9BPxsLM2Q9uPaqG+n7voWEREREKZVhswcgsWLICDgwP69esnOkqOMw/PYNT+UXi/xvsYVm+YSZ0PqtVKWB92E5e23oSzRoYGnbwQ8GZlHnonIiICS6jRSElJwcyZM7F06VIoFArRcQAAmepMTDgwAZ/5fobe1XqLjlOiTt5KwIygK0h9koXBtVzQp38tWNrwNqRERETPsIQaiTlz5qBatWro2rWr6Cg5rMyssPmtzXCwcBAdpcTEJKTjh12RCI2KwzB/H3zcqhKszPXjjwIiIiJ9whJqBB49eoTZs2cjODhY7w53m0oBTclU4c+w61h2+Ca61S2H0PGtUcbBSnQsIiIivcUSagS+++47tGnTBs2bNxcdxeRotBL+Dr6KS3ticLeKNTYMbYbaHrzTERER0auwhBq469evY+HChTh+/LjQHCqtCuZy0zrnMfTsA+xeE4myKRLq+rliRv9asLDifykiIqKC4G9MAzdx4kS8//77qFOnjrAMiZmJGLl/JAbVHISOFTsKy1FSrt5Lwopll+B6LxNeXrboM74OXNxsRcciIiIyKCyhBuzQoUMIDg7G1atXhWW4l3oPQ/cORRXnKvD39BeWoyQkpavw675oqPc8gIutBTqPqItqtUuLjkVERGSQWEINlFarxeeff46JEyfC3V3MwOeRCZEYtm8Y2nu1x8TGE6GQG+dV4CqNFquP3cbckKuoW94Jn4+oh3rVXCHjeJ9ERERFxhJqoNasWYPY2Fh8/vnnQrZ//MFxjAkdgyF1hmBI7SF6d1V+cZAkCaFRcfhuRwQA4Jfe9dGmupvgVERERMaBJdQApaenY+LEifjhhx9gbW1d4ttXa9X48eSPmNh4It6q/FaJb78kXL79BLN3ReJsXDLGtq+K997wgrlCLjoWERGR0WAJNUC//PILypYtK+z2nGZyM6ztuhbmCuO7Gv5RciYWBF6E+ZVk1Kxsj7nj28CRdzoiIiIqdiyhBiY2NhY//PADdu3aBblc3J45YyugWWoNlm6JRGzYAzjJFWjcqwqatfE0ytMMiIiI9AFLqIGZMmUKOnXqhBYtWoiOYhQkScK2I3dwZNN1lMsA6jUrh+69q8HMwjgvsiIiItIXLKEG5Pz581i1ahUuXbpUYttMV6XDUmFplFe+X7j7BDODIqCJSUNnD0f0/qAOHJx5q00iIqKSwBJqICRJwrhx4zBy5Ej4+PiUyDbjM+IxfN9w9KzSE/2qizn/VBdikzLx4+5I7Lz4AENaeGPYB41gZ8n/CkRERCWJv3kNxI4dO3Du3Dls2LChRLZ3O/k2hu4dinpu9dCrSq8S2aaupSvVWHjgBhYeuIF2Ncpg3+f+KO9sIzoWERGRSRJWQteuXYs//vgD58+fR0pKClQqFczM2Inzk5mZiTFjxmDmzJlwcnLS+fYuxV/C8H3D0aNyD4zxGwO5zLCHJtJqJawPv4m9wTcR72WFlUPegF8FZ9GxiIiITJqw1ufs7Izhw4cjIyMDQ4YMERXDIPz8889wdHTExx9/rPNtHbp3COPCxmFE/REYWGugzrena0cjH2HjqsvwfKxBi+pO6P9JfZiZG9/5rURERIZGWAkNCAgAAISFhYmKYBBu376NWbNmYe/evVAodFueJEnCqohVmNpsKjp7d9bptnTtVlwaFgWeh9P1DHiXtkKPLxrAw9tRdCwiIiJ6ise/9dy4cePw7rvvomnTpjrflkwmw/x28w16bMzkTBX+2B0N5Z4HKG1hBv+BNVC/SVmDfk1ERETGiCVUj+3duxf79u1DVFRUiW3TUMuaWqPF2pMx+GVvNKq52+PTd6qgRStPKMwM+3xWIiIiY2VQJXTy5MmwsLAAkH04/9khfWOkVCoxevRoTJ8+HWXKlBEdR6+FRz/CdzuuQKWR8MM7ddG+hpvBlmkiIiIRJEkCJAnQarM/StK/0557PD9t97592LNvH4Ds3lJYBlVCv//+ezg4OIiOUSJ+/fVXmJubY/jw4TpZf7IyGWYyM9iYG+4QRVH3kzArOBJnY5Iwul0VDGhSARbc80lE9EqSVgtoNNmFQqOBpNEC0tNpWi2g1eae9my5p/Og1f673PPTnj3nuc9f/BzpBctKL/5cq4Wk1bx6GenfzwEpe/3aZ9uQ/v1ckvJ/zn/mAdILn5ddzLQvWPa5eU8LngQp99e5it/TfM9Ny7P8K77OUx61WkhA/uXyuWlF4Q3g06efp2o0mF/I5wsroRqNBiqVKqc5Z2VlQa1Ww8LCQug90fXBvXv3MH36dOzYsUMnw1Y9THuIYSHD0KFCBwyrN6zY169rj1Oz8Oeai8C5J6hWzxlzJ7SGk42F6FhEJJgkSYBaDUmthqTRZH+u0UBSawC16t/PNep8Pldnlyy1BpJGnf2L+7lp0P533rPnarOn5fr4dHmNNu8yWg2g1jz38bnip3k6Pddy/52fz3Ka54rZf5d9Vviem55dlApBJgPkckAuh0wuBxSK7KNNz75+fp5cDijkkMmeW06hAOSy7GnPls1nes7zZTLIFHLgldNlOZ9DBsjkiqfrf25b+UyXmSmy1/X88//7HIUCgCz3dnK2/+x5sqfrkWcvq3j6up5fj0wGyJ/f3r/rytmOTJbzkMmfW78MT79+bvl85/+73pxtPMv0bJn/fv38c2TPfY+zv8i93NNOJpPJ8iz372uSITklBXBzK9Q/LWEldOXKlfjggw9yvrazswMAhIaGonXr1oJS6YcvvvgC3bt3R6tWrYp93Tee3MCn+z5F07JN8XEd3Q/5VJyy1Bos2xmNm/vuoZxGjjqdKsK/SyUoFKb9RwvR68opbypVAR5Pl1Nnfw2VKrv0PZun/nc+/jtdrX5u+n+mqf4tg/l+rdbklMv/fg2V6mnxe0W5MjeHTKHI/sWvUEBmZpb9+dOPOdNfNO3pc6CQQ6Z49lwFZHJFdrGRK56bJwcUZpDJ5ZBZWECusM6el7PM0/U995zc85+uQ67I+ZjzPPlLPsqfe87TgpP/MvK8054rjf8WSgVPb6ICkRfhcLxMkoq4D7YEJScnw9HREUlJSUZ/OP7AgQPo2rUroqKiUK5cuWJd99m4sxgZMhL9qvfDiPojDOYHiyRJCDp5F/s2XIVPClCmgSveer8GrGzNRUcjKhJJq4WkVELKyoI2KwuSUgVJmQUpK/uhVSohZSkhqbKXkZTK55ZTZi+rVGZPf/rx33lPHypV7s+ffa1SQatSAsrc5TI/MnPznAfMzbNLm4VF9sdn88zMnk43zy5u5hb/TjMzA8yfLmv23LLmZtnL5jdNkf25zCyfrxWKp89R5Po8uww+/TxnObPcn5v4ETYiXStKVzOoc0KNnVKpxIgRI/DNN98UewHdf2c/Jh6ciPENx6N3td7Fum5dOhfzBDODrsD9egbqOtvi7TG14VrOTnQsMjKSVgspMxPazExIGRnQZmZCm5EJKTMD2sys7I8ZmZCyMv/9OjPrua8zoc3KhJSZ9W+xzMyEVpmVe9qzx3/3GMjlkFlaQm5hkV3yrKyyP1paQG7+dJqlZfZHCwvILMwht7TMLnwWFpA7OP47z9wcMgvznM/luaa/4OPzj6dFE2ZmBvOHKhEZJpZQPTJ79mzIZDKMGTOm2Nd95P4RzGoxC+0qtCv2devCvScZ+Ck4ErsvP8RHLb3x8SBvOPC8T5MmSVL2nr/0dGjT0qFNT4OUng5NWhq06enZ5TE9Hdr0px8z0rOnp2dAm/Hs8dzXzxVOKTPz3w3JZJBZW0NuZQWZlSXkVs8+t3ruoyVkVtbZHy2tIHd1gbnl0+Utn06zsswujpZWkFs+K5aW2Z9bZs+TP/0o4y2LicgE8XC8nrh27Rrq1auHkJAQNGnSRHQcYVKz1Pgr7DoWH7qBzrXLYkJANZRzshYdi16DpFRmF8XUVGhTUqBJTc3+PDU1+/O0NGhTn85PS4M2LfujJi3t6dfpTz+mZV9UAQByOeS2tpDb2GQ/rK0ht7GBzNYGcuvnpz2dbmWd83V2wbSG3Noq+3Nrm+zPrayyl7Ww4B5AIqJC4uF4AyVJEoYOHYoPPvjAZAuoRith7cGbWLv3BqzK2+CfT5qinqeT6FgEQFKpoElOhiYpGdrkJGhSUqBJToY2ORma5BRoU7I/alKSoU1OgSY1JddHKSsre0VmZlDY2UFuZwe5vT3ktjZQ2D792s4OcjtbmLuUgsLWNvvrZyXT1jb3w8Yme+8hiyIRkUFjCdUDK1euREREBDZt2iQ6ihBhVx5i3ZoIVHmsRa+qzhj4qS8Lhg5IGk12mUxMhObJk6cfk7I/T/rPx+RkaJKeQJuUnL0HEoDM0hIKBwfIHR2gsHeA3MEeCgdHKOztoXB2hkWFCpDb22UvY2cPhYN99kf77NLJ4khERM9jCRUsPj4en3/+ORYuXFgspxo8yXyCTE0m3G3diyGdbkXFJmP+qktwu5mB2g4W6DyyNrxruoqOZTAktRrqhARoEhKgfvwYmmefJyRmf0xMgCYhMbtsJiZCk5QESFJ2mXR2zn44OkLh5ASFk2N2kfT2fjrNMfujgwPkjtmfyy0tRb9kIiIyIiyhgo0fPx4tWrRAz549X3td91PvY+i+oWjl0QrjG40vhnS68SglC79tj4DqSDx85GZo8nYV+LXxhFzOvWSSVgvNkydQP4qH+tEjqOMfQRMf/+/Xjx9Dk/AY6vjH0Dx5AkgS5A4OMHN2hsLFBYpSzjBzLgVFqVKwrOwDRalSUDiXgsLZKXsZZ2fIrXmOLRERiccSKlBISAg2bdqEy5cvv/ZhyqiEKAzbNwxtvdpirN/YYkpYvDJVGiw5dBPzw66jeYVSeLtlebTpXhkW1qbxz1CblgbVwzioH8ZCFfsQ6ri4XA9VXBzUjx4BajXkdnYwK10aZq6uMCvtCrPSpWFZvTpsXV1h5lIKCpdnH10gt+CoAUREZHhM47e/HsrIyMDQoUMxc+ZMeHp6vta6TsaexGf7P8Pg2oPxcZ2P9e68O61Wwtbz9/BTcBRc7CyxeFBDNKnkIjpWsdIqlVDHxkJ1/wFUDx5A9eA+1A9ioXrwIKd0alNSILO0hJl7GZi7lYFZmTIwc3ODTUM/mLm5/ftwdeXeSiIiMnosoYJ89913cHZ2xogRI15rPcG3gvHN4W8wqfEk9Kzy+of0i9uxG4/x3Y4IPE7NwoRO1fBWPQ+DPOwuKZVQ3b8P5b17UN27B9Xdpx+fPtSPHgFmZjB3d4d52bIwL1cWZmXLwqpWTZi5u8Pc3R1mZcpA4eSkd38kEBERicASKsDFixfx888/4+jRo1AoFK+1rodpDzHbfzZalS/++8y/juuPUvHrusuwiUpBp4DyGNKxCqzMX++16pomNQ3K27egunMHyjsxUN2Nyf4YEwNVbCxkCgXMy5WDefnyMPfwgGXVqrBr0wbmHuVgXq4czFxds28fSERERK/EElrCVCoVBg0ahLFjx6J+/fqvvb5BtQa9fqhiFJ+ahXlBEUg4+gjV1ArUbuOF5u0qwUJPCqikVEJ55w6ybtyA8tZtKG/dgvL2bShv34YmPh5yBwdYVKgAC09PmHt6wrF+fZh7esLC0xNmZcrw/tNERETFhCW0hP3www9QKpX49ttvRUcpVhlKDZaEXcPJ3Xfgl6FAvVql0bFfddiXshKSR5OaBuW1q8i6fj27cN64CeWNG1DevQuZuTksvL1hUbECLCpWhE3jRtnFs2JFHi4nIiIqISyhJejChQuYNWsWDh48CEsjGXNRo5Ww6cxd/LwnGl3j5WjjaI+A4TVQ1sexRLavzcpC1rVryIq+iqxrV7M/v3oV6vsPsocp8vGBhU8l2DZtCuf334dlJW+YubtzjyYREZFgLKEl5Nlh+M8//xx+fn6Ffv7jjMeIS49DDZcaOkhXeJIkISzqEX7YFYk0pRqTulRHm/KlYF/KCjIdXXikfvQImZGRyIyMRFZkFDKjIqG8eQtyW1tYVq0Cy8qVYdeyFVw+HALLKpVhVqqUTnIQERHR62MJLSHff/891Go1pkyZUujnxiTH4NN9n6JZuWb42uVrHaQrnPMxTzBrVwSiYlMwsm0V9G/iBUuz4jvnU5IkqOPikHn5MjIvXUbm5cvIuHIZmvjHsPDygmX16rCqXg0OXbvCqno1mJUty0PoREREBoYltAScO3cO//vf/3Do0KFCH4a/HH8Zw0OG481Kb2Jcw3E6Slgwt+LTMHtXJI5FPELvVhWxcGBDOFiZv/Z6NcnJyLh4EZkXLiDj3HlkXL4MTUICLH0qwapmTdg2bwaXTz6GVfXqkNvaFsMrISIiItFYQnVMqVRi8ODBGD9+PHx9fQv13MP3DmNc+DgMqzdM6FXwccmZ+C3kKi4evo82Skt0b1IRHTtVL9K6JK0WWdeuIePMmezCeeEClDdvwrxcOVjXqwfbZk3h8umnsKpeDXIbm2J+JURERKQvWEJ17Pvvv4ckSfj668IdRg+6EYTpR6fj26bfomulrjpK93JP0pX4K/wG9oTdRme1FbqZ26JV3yqo7OdW4HVolUpkXrqE9NOnkXH6DNLPnoWkVMK6bl1Y168Pt/HjYV23DsxcXXX4SoiIiEjfyCRJkkSHeJXk5GQ4OjoiKSkJDg4OouMU2NmzZ9G8eXMcPnwYDRo0KNRzQ+6EwMbMBk3LNdVRuhdLy1Jj2eGb+Hv/DXTSWMEtTULDzhVRv50nzCxefu6npFIh4+IlpB07ivRjx5Fx7hzkdnaw8fOFta8fbPx8YVWjBmTmr38Yn4iIiPRDUboaS6iOZGRkoFGjRujVqxemTp0qOk6BZKo0+Pv4HfwZdg2epWwwtIYHzGMy8MZblWDrmP+5rJIkISsyEmlHjyHt+DFknDwFmYUFbJo0gW2TN2DTuDEsvL154RAREZERK0pX4+F4Hfnyyy9hZ2eHr776SnSUV8pQarD6+G38FX4DZRwsMevtumhfw+2FxVGTmoq0I0eQeuAA0g4chCY1FbaNGsG2aVO4jR0Ly6pVOQ4nERERvRRLqA7s2rULy5cvx9mzZ2Gux4ed05VqrDp2GwsP3EBZR2v88HYdtHtB+VTevo2UfSFIPXAA6adPw6JCBdi1aoVy//sB1n5+kFtYCHgFREREZKhYQotZXFwcPvjgA8ybNw8+Pj6vXj49DpEJkWhVvlUJpMuWkqnCqmN3sOjgDVS1tcZ49zLo/UFtyBW5915mXb+OlD17kLx7D5TXr8OmaRPYB3RE2e9mwqJ8+RLLS0RERMaHJbQYSZKEIUOGoFWrVhg06NVDKt1IuoFhe4ehmUezEimhD5IysOzwLaw5fgc1S9vhC/cyeHIhAWWaWECjkSBXAFnXriF55y4k79kNVcxd2LZsAZchH8KudWso7O11npGIiIhMA0toMfrrr79w7tw5nD9//pUX4pyLO4eR+0eid9XeGNVglE5zRTxIxqIDNxB08QHaVHXF/+pXwoPDsbAur0H7LxvC2VFC8qYNeLJpE7KiomDXujVKDx8O21b+UNhxcHgiIiIqfiyhxSQiIgITJkzA9u3bUeoV9ywPiwnDlwe+xBi/MehXvZ9O8kiShINX47H40E2cuPkY7/iWx84RLXBqSQQSlPFo078a3NKvIenX6bi6Zw8sq1WD09s94dClCxQGMgIBERERGS6W0GKQlZWF9957DyNHjkSbNm1euuz269sx49gMzGw+Ex0rdiz2LGlZamw6ew/LD9/Ek3QV3m9SAXN614OrXfYQS2adysHm3F4kj5uC+xkZcOzeHRXXr4NV1arFnoWIiIjoRVhCi8HXX38NuVyO6dOnv3LZSk6V8Ee7P9DIvVGxZrjzOB2BR2/hn1MxqOBig+GtK6Nr3bKwMs8eXF559x4SV65E1oYNQNWqcJswAfbt2nLQeCIiIhKCJfQ17dq1C3/99RdOnjwJiwIMU1TLpVaxbVuSJBy5/hjLDt/CgehH6FCrDJYNboT65R0hl8kgV8iRce4cHi9fgdSQENi1bwevJYthXb9+sWUgIiIiKgqW0NcQExODAQMG4M8//0T16tVLbLvpSjU2n72HFUdu4VFKFvo19sKMHrVQ1tEaty89xj8zT6KWjxrOuxcgMyICTu++i0q7dsGivEeJZSQiIiJ6GZbQIlKpVOjTpw969uyJAQMGlMg2YxLSserYbaw9GYNyTtb4qGUldK9XDlbmCjy+n4ptK87h4bUEVE45CsudQbD9YCDKz/+TQysRERGR3mEJLaJJkyYhLS0N8+bNy3f+g9QHOHjvIHpX6/3a2zp9OwGLD95ESEQc2td0w6KBDdGoojNkMhkyUpUIX38NEYfvwTPjMppGboL7gN5wnreH5ZOIiIj0FktoEWzduhULFy7EqVOnYG1tnWf+1cSrGLpvKFqXbw1Jkl45Zmh+1Botdl9+iMWHbuDaw1T0e8MLoRNaw8Mp9/ZOrTqBuFPX0Dh6I7ze6wrneduhsLMr8msjIiIiKgksoYV08+ZNDB48GIsWLULVfIY1Ohl7Ep+FfoZBNQfhk7qfFLqApmSqsO7UXSw7fBMA8EFzb/Rp5Ak7y9zfqsyICMTNnYvSp86g+geDUWruGpZPIiIiMhgsoYWQlZWF3r174/3330efPn3yzN9zaw++Pvw1vmz0Jd6p+k6h1h2fmoVlh28i8OhtVHazw+QuNdCxZhmY/ed+7srbt/Fo3m9I2bcPzv36odwPP8DM2fm1XhcRERFRSWMJLYTx48dDkiT8/PPPeebtvLET045Ow0+tfoK/p3+B13k3MR2LDtzAP6di0KSSC5YMaoTG3rnvuJSerMT9c3dgu38VkjZthkO3N+GzayfMy5V77ddEREREJAJLaAH9888/WLVqFc6cOQNLS8s88/3K+GFRx0WoW7pugdZ39WEK5odfR9CFB+hQsww2DG2G2h6OuZbRqLQ4G3wdp3fdgkvcebzhngDvzZtg6eNTLK+JiIiISBSW0AI4f/48PvroI/z999/w9vbOd5kytmVQxrbMK9d19WEK5oZcxd4rD/F2Aw/sHtMK3q62uZaRJAnXT8Xi8OqLQOIj+ElnUfOr92HToEGxvB4iIiIi0VhCXyE+Ph49evTAl19+iW7duhV5PdcfpWJeyFXsuhSLd/3KI2x8a5RzyntlfUpiJoJnH8aTuHRUeXIIDT7tCId2PxXpCnsiIiIifcUS+hJqtRp9+vSBr68vJk+eXKR13IpPw7yQqwi6+ADv+Hpg/zh/lHe2yXfZ9DNnEffjbDhkVYB/rwYo3ec7yMz4LSIiIiLjw4bzEuPHj0dcXBy2bt0KuTz7KvWYlBhsiN6Az3w/g1wmf+FzHyRlYO7eq9h87h7eqlcO+8b6w8sl//KpvH0bcT/PQeqhQ3D58EN0+GAw5La2+S5LREREZAxYQl9gxYoVCAwMxMmTJ2H3dPzNK4+vYPi+4ejs3fmFz0vOVOGvsOtYdvgW2tVwy/ecz2fUiYmInz8fT/5ZB8e33oJP8C6Yu7np5PUQERER6ROW0HycOHECI0aMwJYtW+Dz9Er0o/eP4vOwz/Fp3U8xqNagPOdoZqk1WH3sDn7bfxXV3O2x9pMmqOfplO/6Y6PjcWDRcfgcm4/SNT3gvWE9LKtU0fXLIiIiItIbLKH/ERsbi7fffhvTp09H+/btAQBBN4Iw/eh0TGkyBd18cl+cJEkStl94gJ92R8LaXIE5veujdbXS+V5IlJKQgYN/HsSdO1pUTL8En+8mw6lVsxJ5XURERET6hCX0OZmZmXjnnXfQpk0bjB07FgAQcjsEM4/NxNzWc9HMI3dhvHQvCVO3XUZMYjrGdayGd3zLQyHPWz5VSg1OrDiOi6dSUDo5Et26VUS5vhMgk7/4nFIiIiIiYyaTJEkSHeJVkpOT4ejoiKSkJDg4OOhkG1qtFu+//z5u3bqF/fv3w9o6e/ikdFU6YlJiUK1UtZxlE9KUmL0nCpvO3MWHzb0xok1l2Frm3+eVd+5g26xDyEhRomE9LaqN7Ae5dd6hmYiIiIgMVVG6GveEPjVlyhQcP34cx44dyymgAGBjbpNTQNUaLVYfv4Of90ShsXcpBH/WChVfcNGRJikJ8fP/QuKaNfB78214fPspLNxfPZg9ERERkSlgCQWwdOlSzJ8/H0eOHIHbC65OP3EzAd9svQSlWotf+zVAm2r5LycplUhcuxbxf/wJq1q1UHHdP7CqVi3fZYmIiIhMlcmX0H379mHUqFHYsWMHqlevnmd+UroKPwRHYOu5+xjdrgo+bO4NC7O853IqM1SI33MAGX/+BJiZodxPP8K2ZUve6YiIiIgoHyZdQi9fvoxevXph/vz5qNigIr4I/wIzWsyApcISkiRhx8UHmLrtCmqWc8DuMa3gWSrvYPOSVsLFTadxYk8sSidcRKsPBsOpVy/e6YiIiIjoJUy2KcXGxqJLly4YM2YM6neujwG7BqBX1V6wkFvg3pMMfLPlEs7FPME33Wqie71y+e7RvHPiBg6uuIDMdDUaVEhGvZ9HwszeXsCrISIiIjIsJllC09LS0K1bN7Rq1QpthrTBx3s+xmjf0ehX7T0sO3wLP++JQte6ZREyzh9ONhZ5np90NxFh88LxINECVa0eoOm3b8K6oqeAV0JERERkmEyuhCqVSvTq1Qu2trboNrkbJhyYgBnNZ6CmQ0v0XXQMsUmZWDSoIZr5uOZ5rqTVImnLVlyfvwYo3wbvfvQGXJp3EfAqiIiIiAybSZVQrVaLwYMH4+HDh/jxnx/xzclv8Hvb33E9xh2dlx1EjwblsGxwo3zH/Ew7fgIP//cDtE+SUGX8OPh27syLjoiIiIiKyGRKqCRJGD16NE6dOoVDhw7BtbQrXO28MTf4MSIfXMWf7/uiVdXSeZ6nvHULD2fPRvrRY3D59FOUGjQQcktLAa+AiIiIyHiYTAmdNm0atmzZgkOHDqF06dLYdv4+vtl6He2qu2H32FZwtDbPtXzC9Yc48VcIvPbOhnOP7vDZsxtmLi6C0hMREREZF5Moob/99ht+++03HDx4EK7/b+/eg2u6GzWOf7dcGiHJ3pJoXeNanNLm0lBM8qIqotUX6asi9KUc8ko1ZRzHoEQZoeNSHZWmWrRu5URdyihSxqhEL5FQ6tCp+4kTTSIXIkKS80fHbnNSBMle6fJ8ZvaMrP1ba579m5h5stZe69eoKRM2ZHDg52zmRzxN345PVBhbnH+dg0t2c+qiG025StPV6/DsVPn5oSIiIiLy4ExfQteuXcvUqVNJTk6m1LMJL71/gKY2d3a+FUJDDzf7uNJbpaSv2M/h74vwvJnLi4Na0/zlMQYmFxERETEvU5fQHTt2EDMzhp6JPfmh2J2liamMf74N0aGtqVPn95uKio7/xBdLjnPzZjldA515KnoYdVxc7nJkEREREXkYpi2hu3fv5rWpr9Hm7baUlDzD2tTLrBndmSC/BvYxNy9f5tf336dg+w78B43hyfFDcbV5GZhaRERE5NFgyhKanJzMsFnDaD7Bj/IrL+PnHc67bz6Dl/tvZzfLiovJXbmSnOUfUy8khFZfbsO1mR42LyIiIuIopiuhe/fuZfjC4TT696bcuBTFlB4RDH/OD4vFwq2btyjYsZOcJYtx9vam2fKPcA8KMjqyiIiIyCPHVCV0//79/OOtKJpOaIZr9hhWRkUQ2NwGwKkt35Gy4wJNs78j+K1YPPv3x1KnjsGJRURERB5NpimhBw4c4OWho2g/KoEmJfVIHBuKr8dj/G/6aQ4s/47ckvp0anaD4MUzcKlfz+i4IiIiIo80U5TQgwcPMiB6Ko3+uZh+wa2Y2q8DN7Lz+WrOFs5csdKybjF9/yMEj5ZNjI4qIiIiIpighO7fv5/BcZ/Q4OXJxP8jgAFPNyIvKYnv1/1Aoe9TDPjnEzQKCTM6poiIiIj8wV+6hG7Zvp3oFQdp1D2CteNC8TtzjDMDYygrLqbLpEl4vtAbi8Vy7wOJiIiIiEP9ZUvo/LXv82n+BtoEjSXp7+24Oes/+Z/Dh/H517+wDYuijqur0RFFRERE5A4MvT28vLycmTNn0rhxY+rVq0doaCjHjh27535vJb7D6uJPaXW9OwuuXiUvagguTZrQevcuvF8fqQIqIiIiUssZWkIXLFjAihUr2LVrF9nZ2XTv3p2wsDCuXr16x32GLpnEXqcvGXx8MC/9dw8yc91pkfRfPDHjbZxtNgemN5ddu3YZHcH0NMc1S/Nb8zTHNUvzW7M0v7WPoSV02bJlTJo0iU6dOlG3bl1mz55NSUkJmzdv/tPxA5a8SVHJZaK/nUyLPD96h3kwKHEEbm3bOji5+eg/Z83THNcszW/N0xzXLM1vzdL81j6GfSc0Pz+fs2fP0rlz59/DODsTEBBAeno6w4cPr7RPcFYgLW74EfxvEBDTDydXF0dGFhEREZFqYlgJLSgoAMBqtVbYbrPZ7O/dVl5eDkDbogL6Tw+gro+Va8XXofi6Q7I+CkpKSirNu1QvzXHN0vzWPM1xzdL81izNb826Pbe3O1tVWMrvZ3Q1ys/Px2q1kpKSQteuXe3b+/TpQ8eOHVm0aJF928WLF2nWrJkRMUVERESkii5cuEDTpk2rNNawM6FeXl60aNGC77//3l5Cb926RUZGRqVL8Y0bN+bChQt4eHjouZ8iIiIitUx5eTmFhYU0bty4yvsY+pzQcePGsWDBAnr16kXr1q2ZM2cOLi4uDBw4sMK4OnXqVLlVi4iIiIjjeXl53dd4Q0vopEmTKCwspHfv3hQUFPDss8/y1VdfUb9+fSNjiYiIiEgNM+w7oSIiIiLy6DL0OaEiIiIi8miq1SX0QZf1lKr5/PPPCQkJwdPTE4vFwq1bt4yOZCpTpkyhU6dOeHp60qhRIyIjI7lw4YLRsUxl1qxZtG7dGi8vL3x8fAgLCyMjI8PoWKY1cOBALBYLycnJRkcxjbi4OJycnKhfv779FRkZaXQs00lNTaVXr154eHhgtVrp1q0bZWVlRscyhaeeeqrC76+7uzsWi+WOCw/9Ua0uoQ+yrKdUnc1mY9y4cbz33ntGRzEli8XCqlWryM7O5sSJE1gsFvr37290LFMZMmQIP/zwA/n5+WRmZtKnTx/CwsIoLS01OprpfPbZZxQVFRkdw5S6du3K1atX7a/169cbHclUUlNTCQ8PZ8SIEWRlZZGdnc3ixYv1tJ1qcvz48Qq/v/PmzcPb25vw8PB77lurS+j9Lusp9ycsLIzIyEhatWpldBRTio+PJygoCFdXV6xWK5MnT+bIkSNcuXLF6Gim0a5dO2w2G/DblRMnJycuX75Mbm6uwcnM5eLFi0yfPp3ly5cbHUXkvk2ePJlRo0bx2muv4e7ujrOzM126dFEJrSEJCQmMGjUKNze3e46ttSX0Xst6ivzV7N69Gz8/P3tpkuqxY8cOrFYrbm5uTJw4kQkTJuDr62t0LNMoLy/n9ddfZ/r06TRv3tzoOKaUnp6Or68vfn5+DB06lDNnzhgdyTSKiopISUnBycmJzp074+3tTVBQEJs2bTI6mint3buXU6dOER0dXaXxtbaE3s+yniK1XXJyMrNmzeLDDz80OorpvPjii+Tl5ZGTk8PChQsrrMAmDy8hIYHy8nLGjBljdBRTeuWVV/jpp5+4fPkyKSkpWCwWevfura+dVZPc3FzKysr49NNP+eCDD8jKymLatGlERkaSmppqdDzTWbZsGX379qVly5ZVGm/oc0LvxtPTE4C8vLwK269cuUKTJk0MSCTyYLZv386wYcNYs2YNffv2NTqOaTVo0IDY2FhsNhtPPvkkzzzzjNGR/vJ++eUXZs+ezaFDh4yOYlodO3a0/7tJkyasWLECLy8vUlJS6NOnj4HJzMHDwwOAESNGEBwcDMCgQYPo2bMnW7Zs0R+t1SgzM5OtW7eyZcuWKu9Ta8+E/nFZz9tuL+sZEBBgYDKRqlu7di1RUVFs2LCh0kpgUv3Kysq4efMmP//8s9FRTOHAgQPk5OQQFBSEj48PPj4+AEREROjMaA2xWCxYLBb0CO/q4eXlRevWrfX9Twf46KOPaNasWZVuSLqt1pZQ+H1Zz2PHjnH9+nVmzpz5p8t6yoMpLS2luLiYkpISAG7cuEFxcbEeW1FNli5dyhtvvMH27dsJCwszOo4pLVmyhKysLAB+/fVXxo0bh6urK927dzc4mTkMHjyY06dPk5GRYX8BJCYmMm/ePGPDmcTGjRvJzs4GICsri9GjR/P444/TrVs3g5OZx/jx41m1ahUZGRmUlZWxbds29u/fz6BBg4yOZhq3bt1i+fLljB07ljp1ql4ta+3leNCynjVt9erVjBw50v7z7Xndt28fPXr0MCiVeYwfPx5nZ+dKfxXu3LmTkJAQg1KZy549e5g7dy5Xr17F09OT4OBgkpOTadSokdHRTMHd3R13d/dK2318fGjQoIEBicxnzZo1xMTEcO3aNWw2G6GhoSQnJ9svI8vDi42NpaioiP79+5OXl0fbtm3ZsGEDXbp0MTqaaWzdupWcnBxGjRp1X/tp2U4RERERcbhafTleRERERMxJJVREREREHE4lVEREREQcTiVURERERBxOJVREREREHE4lVEREREQcTiVURERERBxOJVREREREHE4lVEREREQcTiVUROQeoqKimDBhgtExRERMRSVUROQe0tLSCAwMrLS9TZs2LFmypNL20tJSAgMDGT16dJWO3717dzZu3Gj/ec+ePTz//PN4enpSr149/P39SUhIQKssi4iZqISKiNxFYWEhp06d+tMS6u/vz9GjRyttX7p0KRcvXuTdd9+95/EzMzPJyMigX79+AKxcuZIBAwYwaNAgzp8/T05ODvHx8cybN48RI0Y89OcREaktVEJFRO4iPT2dunXr0r59ewC+/fZb2rVrR1xcHAEBARw5cqTC+MzMTN5++20WLVpEgwYN7nn8zZs307t3b+rXr09BQQGxsbG88847xMTEYLVacXNzIzw8nKSkJFavXs3XX39dI59TRMTRVEJFRO4iLS2Np59+GoA5c+YQGRnJxx9/TFxcHP7+/hw/fpzS0lL7+AkTJtClSxeGDRtWpeN/8cUXREREAPDNN99QWFjIq6++WmlccHAwLVu2ZOfOndXwqUREjOdsdAARkdrs8OHDeHt706NHD3x9fUlLS8NmswEQEBBAcXExp06dokOHDuzZs4dt27bx448/VjhGWVkZgYGBdO3alYSEBPv27OxsUlNT2bRpEwAFBQUA9uP/fzabjWvXrt0x686dO8nPz2fIkCEP9ZlFRBxBJVRE5C7S0tLIy8vDw8ODvXv34uLiYn+vcePGNGzYkKNHj9KqVStiYmKYNm0abdq0qXCMTz75hD59+nDw4MEK27du3UpISAhWqxXAfsn/yJEjdOvWrcLYa9eucfLkSYYPH37HrOHh4Q/zUUVEHEqX40VE7uB28du6dSsNGzZk5MiRle5Q9/f358iRI8ybNw9nZ2cmT55c4f3CwkI2btzI3LlzyczMrLD/Hy/F3z5WSEgI8fHxlbIsWrSIxx57jKFDhwKQkpLCCy+8QM+ePenZsyfw2132ubm51fb5RURqkkqoiMgdpKen4+bmRmBgIF9++SVHjx5l4sSJFcb4+/uzY8cO5s+fT2JiIq6urhXej4+PZ+rUqTg7O9OyZUvOnj0LQH5+Pvv27WPAgAEVxq9bt46TJ08SExNDcXExZWVlLFy4kAULFpCUlISvry/nzp0jOjqatWvXsm/fPnbv3k1paSn5+flVuhlKRKQ2UAkVEbmDtLQ0/P39cXJywmq1smvXLrZs2cLcuXPtYwICAjh69ChDhw4lJCSkwv7nzp1j/fr1zJ8/n759+3L69Gn73fTbt28nODiYhg0b2sfPnTuX9u3bc/78eZYtW8aGDRs4cOAAkyZNoqioiJdeeono6GjWr1/P2LFj7fu6uLhw4sQJOnTo4IBZERGpHpZyPf1YRKRGREVFMWPGDNq1awdAYmIiWVlZzJgxg4iICP72t7/x5ptv3vdx4+LisNlsxMbGAnDz5k3WrVvHpUuXmDJlSrV+BhGRmqIzoSIiNeDQoUMUFBTYCyhAq1at7A+3f+655xg8ePADHXv06NEkJSURGhpKaGgoly5dIi0tjaCgoGrJLiLiCDoTKiIiIiIOpzOhIiIiIuJwKqEiIiIi4nAqoSIiIiLicCqhIiIiIuJwKqEiIiIi4nAqoSIiIiLicCqhIiIiIuJwKqEiIiIi4nAqoSIiIiLicP8HQQMtGfn7A/AAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "v_a, c_s, wc_p = [wc_1.attrs[k] for k in [\"v_a\", \"c_s\", \"wc_p\"]]\n", "\n", "k = wc_1.k.data\n", "\n", "f, ax = plt.subplots(1, figsize=(9, 7))\n", "f.subplots_adjust(hspace=0, left=0.07, right=0.93, bottom=0.10, top=0.95)\n", "ax.plot(k * v_a / wc_p, wc_1 / wc_p, color=\"k\")\n", "ax.plot(k * v_a / wc_p, wc_2 / wc_p, color=\"tab:blue\")\n", "ax.plot(k * v_a / wc_p, wc_3 / wc_p, color=\"tab:red\")\n", "ax.plot(\n", " k * v_a / wc_p,\n", " v_a * k / wc_p,\n", " color=\"tab:green\",\n", " linestyle=\"--\",\n", " label=\"$\\\\omega = V_A k$\",\n", ")\n", "ax.plot(\n", " k * v_a / wc_p,\n", " c_s * k / wc_p,\n", " color=\"tab:purple\",\n", " linestyle=\"--\",\n", " label=\"$\\\\omega = c_s k$\",\n", ")\n", "\n", "ax.set_xlim([0, 7])\n", "ax.set_ylim([0, 6])\n", "ax.legend()\n", "ax.set_xlabel(\"$k V_A / \\\\Omega_{ci}$\")\n", "ax.set_ylabel(\"$\\\\omega / \\\\Omega_{ci}$\")\n", "ax.set_title(f\"Dispersion relations: $\\\\theta_{{kB}}$ = {theta:3.2f}\")\n", "\n", "f.savefig(\"../../_static/example_dispersion_one_fluid_nb_thumbnail.png\", dpi=100)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.10.12" } }, "nbformat": 4, "nbformat_minor": 4 }