{ "cells": [ { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "6eUXGz8G94ZQ" }, "source": [ "# How to estimate a star's mass and radius using asteroseismology" ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "CqMHYBsI96w6" }, "source": [ "## Learning Goals\n", "\n", "In this tutorial you will learn:\n", "\n", "- What the periodogram of a star with solar-like oscillations looks like.\n", "- How the oscillations can be characterized by two key metrics ($\\nu_{\\rm max}$ and $\\Delta\\nu$). \n", "- How these metrics can be used to estimate a star's mass, radius, and surface gravity." ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "xsLsB9PV-IdJ" }, "source": [ "## Introduction\n", "\n", "Stars of all shapes and sizes oscillate. In many cases, these oscillations are created by standing waves formed inside of the star, which disturb the stellar surface. The study of these types of internal oscillations is called *asteroseismology*.\n", "\n", "This tutorial focuses on a very common type of oscillating star: *solar-like oscillators*. This type includes main sequence stars like the Sun as well as many red giant stars. Stars of this type experience convection in the outer layers of their atmospheres. As material moves up and down in these layers, the turbulent motion forms damped standing waves throughout the stellar interior. Because these waves probe the full stellar interior, we can estimate fundamental properties of a star by studying how these waves disturb the surface of the star. Properties which can be inferred in this way include mass, radius, surface gravity, and even age for the most advanced studies.\n", "\n", "In this tutorial we will explain two key metrics asteroseismologists use to characterize the oscillations of a Sun-like star: the average frequency spacing $\\Delta\\nu$, and the frequency of maximum oscillation amplitude $\\nu_{\\rm max}$. We will also demonstrate how these metrics can be used to estimate the mass, radius, and surface gravity of a star." ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "_-ldXN0k-DaL" }, "source": [ "## Imports\n", "\n", "This tutorial requires **[Lightkurve](https://docs.lightkurve.org)** (for its features to handle *Kepler* data) and **[NumPy](https://numpy.org/)** (to help plot some additional details on figures). We will use **[Matplotlib](https://matplotlib.org/)** for plotting." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "colab": {}, "colab_type": "code", "execution": { "iopub.execute_input": "2023-11-03T14:23:44.665634Z", "iopub.status.busy": "2023-11-03T14:23:44.665128Z", "iopub.status.idle": "2023-11-03T14:23:45.913960Z", "shell.execute_reply": "2023-11-03T14:23:45.913649Z" }, "id": "9vxrKDrS-FjI" }, "outputs": [], "source": [ "import lightkurve as lk\n", "import numpy as np\n", "%matplotlib inline" ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "Ahkq0CCwrokE" }, "source": [ "---" ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "szMCbZPa-YI5" }, "source": [ "## 1. Plotting the Frequency Spectrum of a Sun-Like Star" ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "Zq1yZeY2-aJi" }, "source": [ "The frequency spectrum (also known as periodogram) of a solar-like oscillator like the Sun shows a rich spectrum of modes of oscillation, equally spaced in frequency, colloquially referred to as a \"comb.\"\n", "\n", "We will investigate this spectrum using the Sun-like star KIC 10963065 as an example. This star is also colloquially known as \"Rudy\" (cf. [Lund et al. 2016](https://arxiv.org/pdf/1612.00436.pdf)). In order to get a good resolution in the frequency domain, we start by downloading and combining four quarters of Short Cadence *Kepler* data as follows:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "colab": {}, "colab_type": "code", "execution": { "iopub.execute_input": "2023-11-03T14:23:45.916102Z", "iopub.status.busy": "2023-11-03T14:23:45.915947Z", "iopub.status.idle": "2023-11-03T14:23:47.787837Z", "shell.execute_reply": "2023-11-03T14:23:47.787471Z" }, "id": "_MqdHCoz-bSg" }, "outputs": [], "source": [ "search_result = lk.search_lightcurve('KIC 10963065',\n", " cadence='short', \n", " author='Kepler',\n", " quarter=(2,5,6,7))\n", "lc = search_result.download_all().stitch()" ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "83w84X7O-cFt" }, "source": [ "In one of the companion tutorials, we already identified the location of the modes of oscillation, so we'll zoom in on that region right away when creating the periodogram. We also set the [`normalization='psd'`](https://docs.lightkurve.org/reference/api/lightkurve.periodogram.LombScarglePeriodogram.html#lightkurve.periodogram.LombScarglePeriodogram.from_lightcurve) keyword argument, to match the common conventions in asteroseismology of solar-like oscillators." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 411 }, "colab_type": "code", "execution": { "iopub.execute_input": "2023-11-03T14:23:47.789883Z", "iopub.status.busy": "2023-11-03T14:23:47.789757Z", "iopub.status.idle": "2023-11-03T14:23:48.671621Z", "shell.execute_reply": "2023-11-03T14:23:48.671234Z" }, "executionInfo": { "elapsed": 3856, "status": "ok", "timestamp": 1600230552390, "user": { "displayName": "Geert Barentsen", "photoUrl": "https://lh3.googleusercontent.com/a-/AOh14Gj8sjdnDeqdejfe7OoouYPIclAQV0KSTpsU469Jyeo=s64", "userId": "05704237875861987058" }, "user_tz": 420 }, "id": "YdQYN0b-iL1e", "outputId": "eb3d43a1-499c-421e-808f-43713a602df4" }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuEAAAGQCAYAAAAA1r5rAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAABqxklEQVR4nO3dd3gU1f4/8PcGEiDSexORIkHCtdFREQTEGwQLijS5NAuiYCxfROkgCEpRQPGCSBcIvUhN6F2khiTUAKEkpJCyqZv9/cFv9m7fmd0pm93363l8JLOz55ydMzv7mTNnPqMzGo1GEBERERGRagK0bgARERERkb9hEE5EREREpDIG4UREREREKmMQTkRERESkMgbhREREREQqYxBORERERKQyBuFERERERCpjEE5EREREpLLiWjeA5Ldt2zbs2LED8fHxeOedd9C7d29RrxERERGROjgS7oMqVKiA3r17o02bNpJeIyIiIiJ1cCTcB7Vu3RoAcPLkSUmvEREREZE6GIQrIDs7G+vWrUNcXBzi4uKQmZmJ4cOHo2PHjnbXz8/Px/LlyxEVFYXMzEzUrVsXffv2xTPPPKNyy4mIiIhIDZyOooD09HT8+eefuHnzJh5//HGX68+aNQsbNmxAu3btMGTIEAQEBGD8+PG4cOGCCq0lIiIiIrVxJFyEzMxMnD171uE86n379qFly5YoWbIkAKBixYpYsmQJKlSogEuXLiE8PNxh2XFxcdi/fz8GDBiAN998EwDQoUMHDBs2DH/88QemT59uWverr77CxYsX7ZbzzjvvoF+/fu5+RCIiIiJSEYNwEbZu3YoVK1Zg5MiRpjnVgp07d+Lnn3/Ghx9+iLCwMABAYGAgKlSoIKrsQ4cOISAgAF26dDEtCwoKQqdOnbBkyRIkJSWhSpUqAIBp06bJ9ImIiIiISEucjiJCjx490Lx5c0yfPh1nzpwxLT98+DDmzp2Ll156Cf/+97/dKvvq1auoVasWgoODLZY/8cQTAIBr165JLtNgMCAvLw+FhYWmfxsMBpevEREREZE6GISLUKxYMXz11Vdo1KgRJk+ejNjYWJw5cwY//PADnnvuOYwYMQI6nc6tslNSUuyOmgvLkpOTJZe5atUqvPXWW9i5cydWr16Nt956C1FRUS5fIyIiIiJ1cDqKSEFBQRg9ejRGjRqF8ePHo6CgAE888QT+7//+D8WKFXO73Ly8PAQGBtqtT3hdqt69ezt8CI+z14iIiIhIHRwJlyA4OBgDBw5ERkYGsrOzMWTIEJQoUcKjMoOCgpCfn2+zXAi+hWCciIiIiHwHg3AJ7t69ix9//BG1a9dG1apVMWXKFKSkpHhUZsWKFZGammqzXFhWqVIlj8onIiIiIu/DIFyklJQUjB49GsWKFcPEiRMxceJE5ObmYsyYMcjIyHC73McffxwJCQnQ6/UWy2NjY02vExEREZFvYRAuQmZmJsaOHYucnBxMmjQJlStXRs2aNTF+/Hjcv38f48ePR05Ojltlt23bFoWFhdi+fbtpWX5+Pnbv3o1GjRqZ0hMSERERke/gjZkibN26FUlJSZgyZQpq1qxpWl6vXj2MHj0aY8aMwZ49e0x5wgFgy5YtyMrKMmU3OX78uOnfXbt2xSOPPAIAaNSoEdq2bYslS5bgwYMHqFGjBiIjI5GYmIhPP/1UxU9JRERERGrRGY1Go9aN8HYGgwG3bt3CY489Zvf1a9euoW7duhZpCgcNGoTExES76y9YsADVqlUz/Z2Xl4dly5Zh7969yMzMRN26ddG3b188++yz8n4QIiIiIvIKDMKJiIiIiFTGOeFERERERCpjEE5EREREpDIG4UREREREKmN2FDsKCwuRkpKCUqVKWdxsSURERETkiNFoRHZ2NipWrIiAAOdj3QzC7UhJScGAAQO0bgYRERERFUGLFi1C5cqVna7DINyOUqVKAXi4AYODgzVuDREREREVBXq9HgMGDDDFks4wCLdDmIISHBzMIJyIiIiIJBEznZlBuBPh4eE283nCwsIsnoxJRERERCQVg3AnZsyYwZFwIiIiIpIdUxQSEREREamMI+FERERUpOj1eqSmpqKwsFDrppAfCQgIQIUKFWSbJcEgnIiIiIqEwsJCrFmzBkePHtW6KeTHWrVqhbfffttlHnBXGIQTERFRkSAE4K+99hrq1auHYsWKad0k8iMGgwFXr17F5s2bAQA9e/b0qDwG4UREROT19Hq9KQDv0KGD1s0hP1W3bl0AwObNm/Haa695NDWFQbgTTFFIRETkHVJTUwEA9erV07gl5O+EfTA1NZVBuFKYopCIiMg7CDdhcgoKaU3YBz29MZgpComIiIiIVMYgnIj8wpdffon09HStm0FERASA01GIyE9ERUXhk08+QdmyZbVuChH5uRUrVmD9+vVYs2aNxfKFCxdi48aNGDZsGDp37uxwvQcPHiAiIgLHjx9HUlISihcvjnr16qFDhw54+eWXHU7ZuX37NtavX4/Y2FjEx8ejdu3amDt3rt11d+7cibVr1yIpKQm1atVCv3790KJFC4t1bt68iYULF+L8+fMoXrw4mjdvjsGDB6NcuXI25e3ZswebNm3CzZs3UbJkSTRs2BCjRo1CiRIlAADr1q3Dvn37cO/ePRQUFKB69ero0qULwsLCoNPpTOUYjUZERERg27ZtSE9Px+OPP47BgwcjJCTEps4TJ05g9erVuHbtmmkbhYeHo3LlyqZ+WLlypc37hg4dildffdXudpETg3AiIiIijf3xxx/YuHEjhg4dis6dOztc7/bt2/jmm29QWFiI7t27o0GDBsjPz8fZs2exYMEClC1bFq1atbL73hs3buDkyZN44oknYDQaHc5p3r9/P+bMmYN33nkH//rXv3DgwAF89913mDp1qinY1ev1+Oabb1CpUiV88cUXyM3NxZIlSzBhwgRMnz7dIrHFqlWrsHbtWrz99tsICQlBeno6zpw5Y1F/ZmYmnn/+eTz22GMICgrCmTNn8Ntvv0Gv1+Odd94xrRcREYEVK1agf//+ePzxx7F161aMGTMGP/30E6pXr25aLyoqCj/99BPeeOMN9OvXD9nZ2bhw4QLy8vIsPmtQUBAmT55sscy8HCUxCCciIiLS0NKlS7F27Vp89NFH6NKli9N1f/zxRxgMBsycOROVKlUyLX/uuecQFhYGvV7v8L0tWrQwBegzZ87E5cuX7a63YsUKvPDCC+jbty8A4F//+heuX7+OP//8E+PGjQMAbN26FVlZWZg9ezYqVKgAAKhZsybCw8Nx7NgxtG7dGgBw69YtrFy5Et9++y2aNWtmqqNt27YWdb733nsWfz/99NNISkrCnj17TEF4Xl4eIiIi8MYbb+D1118HADz55JP48MMPsW7dOgwdOhQAkJGRgV9//RVDhgzBv//9b1OZLVu2tPmsAQEBdkfR1cA54UREREQaWbFiBVavXo0PPvjAImC058KFC4iLi8Pbb79tEYALqlataspjbY+YJzzevXsXCQkJeP755y2Wv/DCCzhz5gzy8/MBAFevXsXjjz9uCsABoGHDhihTpgyOHz9uWrZ7925Uq1bNIgAXq0yZMigoKDD9ffHiRej1eou2BQYGonXr1vj7779Nyw4ePIjCwkJ06tRJcp1q4ki4E8wTTkREREpZtWoVVq5cicGDB6Nr164u1z9//jyAh6PeSrl16xYAoHbt2hbLH330URQUFODu3bt49NFHkZ+fj8DAQJv3BwYG4ubNm6a/Y2Nj8dhjj2HVqlXYvHkzsrKy0LBhQwwaNAiNGjWyeb/BYEBubi4uXLiAqKgovPvuu6LatmnTJuTm5qJEiRKIiYlB7dq1ERkZiVWrViE5ORmPPfYY3nvvPZuTgby8PPTp0weZmZmoVasWunfvjldeeUXiVnMPg3AnmCeciIiIlJCTk4Nly5ahc+fO6N69u6j3JCcnAwCqVKmiWLsyMzMBAKVLl7ZYLvwtvF6jRg3s2bPHFPgCQGJiIlJTU1GqVCnT+1JTU3HlyhXEx8fjo48+QokSJbB69WqMGTMG8+fPR/ny5U3r3r59Gx988IHp7549e5qmnQh1BwYGIigoyKZtRqMRmZmZKFGiBNLS0pCQkIBly5ZhwIABqFChArZu3YpJkyZh9uzZeOyxx0yfoX///qhfvz7y8vKwb98+zJkzB1lZWXjzzTc93JKuMQgnIiIinzNlyhQkJibKWmbVqlXx9ddfy1JWUFAQGjZsiH379uHll1/Gk08+KUu5annllVewefNmzJ07F/3790dubi7mzp1rkckEeJjNJDs7G99//z0ef/xxAECjRo0waNAgbNmyxTTvHHh4cjFjxgzk5OTgwoULiIiIgE6nQ58+fSS1rbCwENnZ2fj8889N88CbNm2KDz74AGvXrkV4eDgAoH379hbva968OQoKCrB69Wp069YNxYsrGyYzCCciIiKfI1ewrJSAgACMHj0ao0aNwoQJEzB16lSn87kBmOaBJyUloWbNmoq0SxjxzsrKspjvbT1CXrt2bXz66af47bffEBUVBQBo3bo1mjVrhuzsbIvyypQpYwrAgYdzvevVq4cbN25Y1B0YGIiGDRsCeBg0BwcHY+HChfj3v/+NChUqoHTp0sjPz0deXp7FaHhmZiZ0Op2pbcL///Wvf5nWKV68OEJDQxEfH+/08z///PM4dOgQ7ty5g0cffVTMJnMbb8wkIiIi0sAjjzyC8ePHo2zZshg7dizu3r3rdP3Q0FAAwD///KNYm4T51sL8a8GtW7dQvHhxi/R9HTp0wNKlSzFnzhz88ccfGDVqFO7cuWMx17tOnToO6xJu8nSkfv36KCwsxL179yzalpCQYNO2KlWqmKbFOKvTOkWhlhiEExEREWmkfPnymDhxIgBgzJgxSE1NdbhukyZN8MQTT2D16tVISUmxeT0pKQnXr1/3qD3Vq1dHrVq1cOjQIYvlBw4cwFNPPWVzM2ZgYCAee+wxVKpUCWfOnMHt27fx8ssvm15v3rw5MjIycPXqVdOy9PR0XL16FfXr13falujoaOh0OlSrVg0A0LhxYwQHB+PgwYOmdQoKCnD48GGLm1WbN28OADh9+rRpWX5+Ps6fP48GDRo4rXP//v145JFHUKNGDafryYHTUYiIiIg0VK1aNUyYMAEjR47E2LFjMWXKFDzyyCN21/38888xatQofPbZZ3j99ddND+s5f/48tm7dis8++8zhtJacnBxTKr+kpCTo9XpTsB0aGmp60mWvXr3w448/onr16qaH9cTFxWHq1KkWZa1YsQKhoaEIDAxEbGws1qxZg169ellkL2nVqhUaNmyIqVOnol+/fggKCsKaNWtQvHhxU7a5rKwsjB8/Hi+99BJq1KgBg8GAc+fOYdOmTejSpYtpWkxQUBB69OiBlStXoly5cnjsscewbds2ZGRkWNxI2aBBA7Rp0wZz585FZmam6cbMtLQ0i/VGjBiBDh06oHbt2sjLy8PevXtx5MgRDBkyRPH54ACDcKeYopCIiIjU8Nhjj2Hs2LH49ttvMWHCBEyYMMHuejVr1sSsWbOwdu1a/PXXX0hKSkJgYCDq1auHIUOGmEaB7Xnw4IFFIA3A9Pd3332Hpk2bAgDatWuH3NxcREREICIiArVr18aoUaMsHmqj0+kQHx+P3bt3IycnB7Vr18ZHH32Ejh07WpQfEBCAsWPHYsGCBZg7dy7y8/PRpEkTTJ061SK4rlmzJjZs2IDk5GSUKFECNWrUwMcff2xz82SPHj0AAOvXr8eDBw9Qr149TJgwweYpl5999hkWL16MxYsXQ6/Xo0GDBpg4caLFCUqNGjWwceNGpKWlAQDq1q2Lzz//HC+99JLDbSgnndFoNKpSUxGi1+vRs2dPrFq1iikKiXxEs2bNsG7dOqdzBYnIe928eRMzZsxAeHi44jfMETnjbF+UEkP6/Ej4zJkzERkZ6fD1P/74w+5Tp4iIiIiIlOLzQfirr76Kp59+2mKZ0WjEvHnzULVqVQbgRH6EF/6IiMhb+HwQHhISYjGHCQAuXLiA3Nxc1eb8EBERERGZ88sUhfv27YNOp0O7du20bgoRERER+SG/C8ILCgpw8OBBhISEmPJOEhERkXcTspUZDAaNW0L+TtgHrTPoSeXz01GsnTp1ChkZGaKmouj1elFlBgYG2iSvJyIiIvkI6eyuXr3q8vHuREoSHjwk7JPu8rsgfN++fShevDief/55l+sOGDBAVJm9evVC7969PW0aERERORAcHIxWrVph8+bNAIB69eqhWLFiGreK/InBYMDVq1exefNmtGrVyuM01n4VhGdnZ+PYsWN45plnULZsWZfrL1q0SNQG5ig4ERGR8t5++20AMAXiRFpo1aqVaV/0hF8F4UePHpWUFSU4OJgP6yEiIvISAQEB6NmzJ1577TWkpqaisLBQ6yaRHwkICECFChVkiw39Kgjfu3cvSpUqhRYtWmjdFCIiInITB8nIF/hNdpQHDx7gzJkzaNWqFUqWLKl1c4iIiIjIj/lNEH7gwAEYDAY+oIeIiIiINOc3QfjevXtRvnx5PPXUU1o3hYiIiIj8nN/MCf/hhx8kvyc8PNwmEXtYWBjCwsLkahYRERER+SG/CcLdMWPGDN74QURERESy85vpKERERERE3oJBOBERERGRyhiEE5HfMBqNWjeBiIgIAINwIiIiIiLVMQgnIiIiIlIZg3AiIiIiIpUxRaETzBNOREREREpgEO4E84QTERERkRI4HYWIiIiISGUMwomIiIiIVMYgnIiIiIhIZQzCiYiIiIhUxiCciIiIiEhlzI7iBFMUEhEREZESGIQ7wRSFRERERKQETkchIiIiIlIZg3Ai8hs6nU7rJhAREQFgEE5EfsRoNGrdBCIiIgAMwomIiIiIVMcgnIiIiIhIZQzCiYiIiIhUxhSFTjBPOBEREREpgUG4E8wTTkRERERK8Isg/PLly1i5ciWio6ORl5eH6tWr45VXXkG3bt20bhoRERER+SGfD8JPnTqFiRMnon79+ujZsydKlSqFO3fuIDk5WeumEREREZGf8ukgXK/XY+bMmWjevDlGjhxpM7+biIiIiEgLPh2V7tu3D2lpaejXrx8CAgKQk5ODwsJCrZtFRERERH7Op0fCT58+jeDgYCQnJ2Py5MlISEhAyZIl0b59ewwePBhBQUFO36/X60XVExgYiMDAQDmaTERERER+wKeD8Nu3b8NgMGDSpEno1KkT3nvvPZw7dw5btmxBVlYWvvzyS6fvHzBggKh6evXqhd69e8vRZCIiIiLyAz4dhOfk5CA3NxevvvoqPvjgAwBAmzZtUFBQgO3bt6NPnz6oWbOmw/cvWrRIVIpCjoITERERkRQ+PSdcmG7y4osvWixv164dACAmJsbp+4ODg0X9xyCciIiIiKTw6SC8YsWKAIDy5ctbLC9XrhwAIDMzU+0mERERERH5dhDeoEEDALDJCZ6SkgIAKFu2rOptIiIiIiLy6SD8+eefBwDs2rXLYvnOnTtRrFgxNG3aVItmEREREZGf8+kbM+vXr49OnTph165dMBgMCA0Nxblz53Do0CG8/fbbqFSpktZNJCIiIiI/5NNBOAAMHToUVapUwe7du3H06FFUqVIFgwcPRvfu3V2+Nzw83OYpm2FhYQgLC1OquURERETkB3w+CC9evDh69eqFXr16SX7vjBkzRKUoJCIiIiKSwqfnhBMREREReSMG4UREREREKmMQTkRERESkMgbhREREREQqYxBORERERKQyBuFERERERCrz+RSFnmCecCIiIiJSAoNwJ5gnnIiIiIiUwOkoREREREQqYxBORERERKQyBuFERERERCpjEE5EfsNoNGrdBCIiIgAMwomIiIiIVMfsKE4wRSERERERKYFBuBNMUUhERERESuB0FCIiIiIilTEIJyIiIiJSGYNwIiIiIiKVeTwnPDIyUo52AAA6dOggW1lERERERN7K4yB81qxZ0Ol0crSFQTgRERER+QVZsqPUrVsXrVq1cvv9R48exfXr1+VoChERERGR15MlCK9Xrx569erl9vvv3bvnlUE484QTERERkRI8DsKDg4NRokQJj8oICgpCqVKlPG2K7JgnnMi3yDV1joiIyFMeB+F//vmnx40YOnQohg4d6nE59pw7dw6jRo2y+9r06dMREhKiSL1ERERERI74zRMzX3vtNTRs2NBiWY0aNTRqDRERERH5M78Jwps0aYK2bdtq3QwiIrIyd+5cPPfccx7d4E9EVNSo+rCe7t27q1mdDb1eD4PBoGkbiEg7RqNR6yaQHZs2bcL58+e1bgYRkapUHQnX8gdw9uzZyM7ORkBAAJo0aYIBAwbYTE+xptfrRZUdGBiIwMBAOZpJRERERH5A1SBci8wExYsXR5s2bdCsWTOULVsWN27cwPr16zFy5EhMmzYN9evXd/jeAQMGiKqjV69e6N27t1xNJiIiIiIfp+mc8O7du6NGjRoYNWoU6tSpY/FabGwsvvrqK2zcuNGjOho3bozGjRub/m7ZsiXatm2LTz75BEuWLMH48eMdvnfRokWiUhRyFJyIiIiIpFB1Trg1o9GIwMBAfPHFF9i7d69q9dasWROtWrXC2bNnnc4RDw4OFvUfg3AiIiIikkLTIBwAhg0bhnfeeQezZs3C/PnzVbtxsnLlyigoKEBubq4q9RERERERCTRPUajT6dCjRw+EhIRg+vTpuHz5MkaOHKl4vXfv3kVQUBBKliypeF1EREREROY0HwkXhIaGYtasWQgMDMTw4cNx+vRpWcp98OCBzbJr167h+PHjeOaZZxAQ4DWbgIiIiIj8hOYj4eYqVKiASZMmYenSpVixYoUsZU6bNg1BQUEICQlB+fLlcePGDezYsQMlSpRA//79ZamDiIiIiEgKTYPwyZMno3bt2hbLAgIC0L9/fzRr1gx37971uI6WLVti37592LhxI/R6PcqVK4fWrVujV69eqFmzpsflExERERFJpWkQ3rRpU4evNWnSBE2aNPG4jm7duqFbt25uvTc8PNxmukpYWBjCwsI8bhcRERER+S/Vn5j5559/in5ypk6nw7vvvqtwqxybMWOGqDzhRERERERSqD4SvmPHDptlKSkpKFu2LIoXt22OlkE4EREREZESFAnC//77b/zyyy9YsGCBxXKdTodFixZZLDMYDHjjjTcwduxYNGjQQInmEBERERF5FUXy8+Xk5CApKUmJoomI/J7BYFDtwWZERKQMSSPhy5YtE7XerVu33GoMERG5NmPGDOTk5GD06NFaN4WIiNwkKQhfvXo1HnnkEZc3K+bl5XnUKCIiJYi9KdzbJSYmIjc3V+tmEBGRByQF4dWrV0eTJk0wfPhwp+sdOnQI06ZNk9QQnU4naX01MEUhERERESlBUhAeEhKCixcvilrX3oiT0WjEyJEj7a7/008/oVSpUjbLp06dKqWJsmKKQiIiIiJSgqQgvFu3boiOjna5XmhoKCZPnmz3terVq4taRkRERETkqyQF4Q0aNBCVRrBcuXJ2n4ap0+kwYsQIKVUSEREREfkcRVIUyiEyMhIffvih1s0gIiIiIpKd6k/MNNetWzenN2RWrlxZxdYQEREREalD0yD8448/tgnCc3NzcfnyZRw+fJhTV4iIiIjIJ2kahL/yyisOX6tVqxYWLFiA2bNnq9giIiIiIiLlaRqEO9OoUSOsWrVK0zYwTzgRERERKcFrg/BTp07h2Wef1bQNzBNORERS3Lx5E5988gk2bNigdVOIyMtpGoTPmDHDZllhYSESEhJw7do1tGzZ0mKd8PBwNZtHREQkSVpaGm7duqV1M4ioCNA0CL9//77d5aVKlcKTTz6JjIwMZGRkqNwqIiIiIiJlqRqEWz/K/rvvvlOzeiIiIiIir6Dqw3o2bdqkZnVERERERF7J4yA8IyMDOTk5HpWRk5PDaSdEpDjrq3GkvcLCQiQnJ2vdDCIi1Xk8HaVv3754+eWX8emnn7pdxq+//oqoqChs3LjR0+bIiikKiYiUxRMjIvJXHgfhRqPRZw+iTFFIRERS+OrvIRHJT5YbM6Ojoz16smV0dLQczXBp1apVWLZsGerUqYO5c+eqUicRERERkTVZgvA7d+7gzp07HpWh0+nkaIpD9+/fx5o1a1CyZElF6yEiIiIicsXjIHzy5MlytENxv//+Oxo1aoTCwkKkp6dr3RwiIiIi8mMeB+FNmzaVox2KOn/+PA4dOoTZs2dj/vz5WjeHiIiIiPycqnnCtWAwGDB//nx07twZdevW1bo5RERERETaPrZeDdu3b0dSUhImTZok+b16vV7UeoGBgQgMDJRcPhER+RZmRyEisXw6CE9PT8fy5cvRs2dPlCtXTvL7BwwYIGq9Xr16oXfv3pLLJyJ1KX0DOBGDcCISy6eD8GXLlqF06dLo2rWrW+9ftGiRqDzhHAUnIvJe58+fx9y5c/HLL79o3RQiIhOfDcJv376NHTt2YPDgwUhJSTEtz8/Ph8FgwL179xAcHIwyZco4LCM4OJgP6yEiKuLu3r2LEydOaN0MIiILPhuEJycno7CwEL/99ht+++03m9cHDx6Mbt26YciQIRq0jojIfZcuXUKpUqW0bgYREXlA9iB8wYIFKF26NN599125i5akTp06GDVqlM3yZcuWITs7G0OGDEGNGjU0aBkRkWdu3bqldROIiMhDsgfhW7duRcuWLeUuVrJy5cqhdevWNss3bdoEAHZfIyIidfFmWSLyV7LnCa9UqRIKCwvlLpaIyGPMXEFERN5C9pHwVq1aITIyEnq93itvapwyZYrodcPDwxEQYHmeEhYWhrCwMLmbRURERER+RPYgvHfv3jh37hzGjx+P999/H/Xr15e7CtXMmDHDK08kiIiIiKhokz0Inzx5MgIDA3Hx4kWEh4ejQoUKqFKlCoKCghyuT0RERETkT2QPws+dO2f6t9FoREpKikWebnO8IYeIiJTG3xoi8kaKpCgkIiIS4969e1o3gYhIE7IH4VWrVpW7SCIi8lGOrpTKiSPhROSNZE9RSERE5K+YBpOIxFLssfXp6enYu3cv4uLikJ6ejqeeegpvvfUWACA+Ph53797FU089hZIlSyrVBI8xRSEREUnBIJyIxFIkCD948CB+/vln5OTkwGg0QqfToVKlSqbXk5OT8d1332HEiBFo3769Ek2QBVMUEhEREZESZJ+OEhMTgx9++AHFihXDwIED8eOPP9qMDDz11FMIDg7GkSNH5K6eiIiKEI4cE5G/kn0kfPXq1dDpdJgwYQIaNGhgd51ixYqhfv36iI+Pl7t6IiIiIiKvp8hIeEhIiMMAXFChQgWkpqbKXT0REZEFZkchIm8kexCem5uLcuXKuVwvMzNT7qqJiIg0xek1RCSW7EF4pUqVcOPGDafrGI1GxMfHo1q1anJXT0RERDKJjo7Gzp07tW4GkU+SPQh/9tlnkZCQgP379ztcZ+fOnbh//z6aNWsmd/VERFQEKTllxLrs69evK1aXr9m9ezd++eUXrZtB5JNkvzGzR48e2LdvH2bOnImrV6+iVatWAICcnBxcuXIFR44cwbp161CuXDl0795d7uplxTzhRL6FUwW8l5p906NHD5w8eVK1+oiI7JE9CK9cuTLGjh2LKVOmYN26dVi/fj10Oh0OHz6Mw4cPw2g0oly5cvj2229Rvnx5uauXFfOEExEREZESFHlYT0hICH755Rfs2rULp0+fxr1792A0GlG5cmU8/fTT6NKlCx555BElqiYiGRUWFtpcDSJSQn5+vtZNICJSlWKPrQ8ODkb37t29fsoJETnWokULn7psf+XKFZfpU0kbCxcuxEcffaR1M4iIVCP7EFd6errcRRIRyeLo0aOq1LN//358//33qtRF3oX3HRCRWLIH4f369cMnn3yC+fPn49ChQwzKicjvxMXFITIyUutmEBGRF5N9Okq9evVw9epVxMfHY9u2bQCARx99FKGhoWjatClCQ0NFPcyHiIiIiMhXyR6Ez5w5E3q9HufPn8e5c+dw7tw5XLt2DTdu3MBff/0FAKhVqxaaNm2Kpk2b4vnnn5e7CbJhikIicgenJJAvUWt/FupRMmc8kTdR5MbM4OBgtGjRAi1atAAAi6D8/PnzuHr1KhISErB9+3avDsKZopCI3MVAQhyesHg3NffjqVOnolSpUhgxYoRqdRJpSbHsKBaVFC+OEiVKoESJEggKCkJAQAAMBoMqX+74+HisXLkSly9fRmpqKkqUKIE6dergzTffNJ0kEJF/UDPgY3DpPXhCVDQkJCQwfTH5FUWC8Ly8PFy8eNE08n3p0iUUFBTAaDSiSpUqaNeunWk6itKSkpKQnZ2Nl19+GRUrVkRubi4OHz6MiRMn4uOPP0aXLl0UbwMR+RcGfURE5IrsQfjIkSNtgu4XXnjBFHRXrVpV7iqdatasGZo1a2axLCwsDJ999hk2bNjAIJzIhX/++QfPPPOM1s0oUjgK7r/Y90QkluxBeHR0NHQ6HerUqYN3330Xbdq08bpRoWLFiqFy5cq4dOmS1k0h8nrJycmq1HPw4EE8/fTTKF26tOxlX758WfYyXfG2454/Y18QkTeSPQhv0aIFoqOjER8fj2nTpqFUqVJo0qSJKT1h/fr1NTkg5uTkIDc3F3q9HseOHcPff/+NF154QfV2EJF9I0aMwE8//YQ2bdrIXvbixYtlL5PkwQCZiPyV7EH4t99+C6PRiGvXruHs2bM4f/48oqOjceLECeh0OpugXK1HSC9cuBDbt28HAAQEBKB169b48MMPnb5Hr9eLKjswMBCBgYEet5GIlMWpAt7H1/okISFB6yYQURGhyI2ZOp0O9erVQ7169fD666+bgnIhb/g///yDkydPAgA2btyoRBNsdOvWDW3btkVKSgoOHDiAwsJC5OfnO33PgAEDRJXdq1cv9O7dW45mEnmdtLQ0rZsgG18L+JSSmJiIt99+G/v27dO6KUXOpEmTtG6C7Pi9IVKGKikKk5KScP36ddN/QvCr5mXIRx99FI8++igAoEOHDhg9ejQmTpyIH3/80WE7Fi1aJCpPOEfByZfdu3dP6yaQyvR6PbKysrRuBnkBThciUo4iQXhiYqLFEzOTkpIAPDybLl68OJ588knTY+y10rZtW8ydOxcJCQmoXbu23XWCg4P5sB4iFfnKDz5HDomIyBXZg/DBgwdbBN2BgYFo3LixaQ5448aNERQUJHe1kuXl5QEQP++biIouLYJiXzmhIO+VmZmJRx55hPsaURElexCekpKCJ5980pQXvFGjRpoG3WlpaShfvrzFsoKCAkRGRiIoKMg0RYWIfB+fmOl91NhOvhqkvvTSS/j1119tnoUhJ+7HRMqRPQhftWqVV82Rnjt3LvR6PUJDQ1GxYkWkpaVh7969uHXrFgYNGoRSpUpp3UQi8jG+GvQVVb7cH64SDMhBre3HgJ/8jexBuDcF4ADwwgsvYNeuXdi2bRsyMjJQqlQpNGjQAP/5z3/QsmVLrZtH5PXUDGB8JVgq6sEErxiQVnzlGEAkhqLZUWJiYnDhwgXTE/cqVaqEJk2aICQkRMlqLbz44ot48cUX3XpveHg4AgICLJaFhYUhLCxMjqYRkcp4QkFKioqK0roJRFSEKBKEJyQkYMaMGaZHRQsjHcKPUoMGDfD555+jZs2aSlQvmxkzZjA7ChERiSIMOKlpzpw5aN26ter1EpHnFLkx8+uvv0ZaWhoqVqyItm3bolq1agAepi48dOgQLl26hK+//hozZ85ExYoV5W4CEZEFTnkgXxUbG6t4Hfz+EClDkRsz09LS0L17d7z33ns2c8T/85//YPHixdi4cSPWrFmDDz74QO4mEBHZxWDC+3Dajndj/xApJ8D1KtKcPHkStWrVwqBBg+zepFm8eHEMHDgQtWrVwvHjx+WunoiKMF/6wReel0BERGSP7EF4amoq6tev73QdnU6H+vXrIzU1Ve7qiagI+eWXX5CYmKh1M8gKs6O4x5c+i1a4DcmfyB6EBwcH4/79+y7Xu3//Pm96JCoClBydXrhwIa5fv65Y+VphIEEk3bFjxxAdHa11M4hUI3sQHhISgosXL+LEiRMO1zl58iQuXryoaqpCIiLyT740zcnX3blzR+smEKlG9hsze/TogZMnT+K7777D888/j3bt2llkR9m/fz/2798PnU6HHj16yF29rJgnnMg3cmurPTLNoM97bd26VesmEBEBUCAIDwkJwfDhwzF37lzs27cP+/fvt3jdaDQiKCgIH3/8sdePhDNPOJFvBJRCEK5WMO5r01EKCgrwyiuvYM+ePVo3xWNjx47VuglERAAUelhP+/bt0bRpU+zYsQPR0dFISUkBAFSsWBFNmjRBp06dUKVKFSWqJiKZKR1QmpfvCwG/LyosLMSDBw+0bobX87WTLwDsdyIFKfbY+sqVK6NPnz5KFU9ERAphdhQSrF+/HsWLKxYqEPk12b5ZJ0+exNGjR5GUlITAwEDUrVsXHTt2RPXq1eWqgog04AtzwkkaBsbuycvL07oJiigsLNS6CUQ+SZYg/IcffsCBAwcA/O/gfeLECaxfvx5fffUVWrZsKUc1REREXkuNR8gTke/wOAjfuXMn9u/fj2LFiqF9+/aoV68esrOzceLECcTExGDmzJlYuHAhHnnkETnaS0Q+hlMfiIjIH3kchEdGRkKn02HcuHF46qmnTMvffvttzJo1C1FRUThy5Ag6duzoaVWqY4pCIiJlqTEFydemOc2bNw9Dhw5VrT6evLonIyMDZcqU0boZ5MU8DsKvX7+ORo0aWQTggnfeeQeRkZFF9ol4TFFIRKSsQYMGad2EIuf3339XNQgn97Rv3x7r1q1DnTp1tG4KeSmPn5iZnZ2NGjVq2H1NuClTr9d7Wg0REZFbfG0knIoOg8GgdRPIi3kchBuNRpspG6bC//9yXsoiIm/AYMw9PIYTEcnP4yCciHybL6QoVPuJmb6ksLAQbdu21boZREQ+R5YUhZGRkYiMjLT7mk6nc/r6xo0b5WgCEfkAjlR7H564EPcBImXIEoTzC0qknaysLAQGBiIoKEjrppAP4okREZEyPA7CN23aJEc7iMhNH3/8MVq2bImPPvpIkfIZhJEvunXrFmrXrq11M4jIj8n22HpfxDzhVBRkZGQgNzdX62aQD/G1q5v2TiTXrl2L4cOHK14P+Tdf+y6RvBiEO8E84eSuli1b4tixY1o3QxZq/ogoHcTwB5GIiLyFzwfhcXFxiIyMxNmzZ5GYmIgyZcqgUaNG6NevH2rVqqV188hHqZkbloGlaykpKQC4rUhZauxfmZmZitdBROrw+RSFa9euxeHDh/HUU09hyJAh6NKlCy5cuIARI0YgPj5e6+YRua2wsBD79u0DoOwIsi9cYlf7gRkM9kkp8+fP17oJJIHSx88RI0YoWj4py+dHwl9//XV88cUXCAwMNC174YUXMGzYMERERODzzz/XsHVE7svOzsbnn3+u+COR1ZwionSecCIlqXHCyn2ZzB08eFDrJpAHfH4kvHHjxhYBOADUrFkTderUwc2bNzVqFRH5Mk+CsZiYGKxatUrG1niGQR+p4cGDB1o3gUh1Ph+E22M0GpGWloayZcs6XU+v14v6Lz8/X6WWE5En1JpP60ngevToUfz2228OX4+NjXW7bDEKCgoQHR2taB1E1tavX691E4hU5/PTUezZu3cvkpOT0adPH6frDRgwQFR5vXr1Qu/eveVoGvkAtW+cMhqNPjFvWw1HjhzRugke27lzp6Llx8bGYuLEiYrW4at41YCscZ8gZ/wuCL958yZ+/fVXhISEoEOHDk7XXbRokagUhdbTXci/xcTEqFofg3Ai53zp++FLn8WcLwSr9+/fR/ny5VG8uN+FVuQmv5qOkpqaigkTJiA4OBgjR45EsWLFnK4fHBws6j8G4UTuM//xXbRokYYtKRp8IVjxVb4aIJM4r776Knbt2qV1M6gI8ZvTtaysLIwbNw5ZWVmYOnUqKlWqpHWTiGShdFB26dIl1UbbDx06pHgd5F/y8/Nx8uRJrZtRJM2dO1e1unzh5NJoNKKwsFDrZlAR4hcj4Xl5eZg4cSISEhIwZswYxVO6EalBjdR+ALBnzx6kpaUpVj5Jo/TJkK+N5sbHx/MKi5vU3G6+EIT7kmbNmiEjI0PrZvg8nw/CDQYDpk2bhpiYGIwcORIhISFaN4mIVOZLgSWDFe/lS/uZ2nxlv7b+HEV5nygoKNC6CT7P56ej/P777zh27BhatGiBjIwMREVFWbzevn17jVpG5JmifHA3p+SP75UrVzTJP+wrAYUvcNQXSvQR+52IpPD5IPzq1asAgOPHj+P48eM2rzMIJ1+g5lMti5KIiAicOXMGJUuW1Lopkmi5va33paLa91Q07N27FxEREXj66ae1boos+P0hKXw+CJ8yZYrb7w0PD0dAgOWMnbCwMISFhXnaLCL6/3xlRF9QUFAAg8Hg9vt9bXv4E/addHfu3MHZs2d9JggnksLng3BPzJgxQ1SecH+ydetWnoRoaNy4cRg3bpzFMo60OKbFtpkxYwZWr16ter2kPa2+iytXrkTLli1Rr149Ter3lF6v17oJsuHxmKTw+RszSV5jx47VuglFRkFBgeyp0bZs2WL6d0JCAgA+rMcVnU6n6vZJSUlRra6cnBzZy+S+5L1yc3PtBnk//vgjPvvsMw1aJB9/D16HDx+udRNIAwzCiWQmBDFJSUn48MMPFavnvffeU6xsa0r+QPr7j69U5kHyihUrNGwJWVP6BKZPnz7YvHmzonWoTdhm/n4c8MZnJPh7n6iBQTiRzJQYnbRHOEDyQOmYr28jPhjENV/q+8TERGRlZWndDEWo0U+FhYU+tT9Q0ccgnEhmWlxW5BQCx3xt25gHEQwo3KdF6kr6n71796JPnz44c+aManWOGjUK//3vf1Wrj8gVBuGkury8PD4EQAZqBpc3btxQra6iztcC4/z8fK2boIgTJ07IXqbSfe9LJ5RxcXGIjY3Frl27AAB37941vabUjZoJCQm4f/++ImWbt19pw4cPx4YNG1Srj5TD7ChOeEOKwsLCQpw6dQrNmjVTrU6lff7552jatCnef/99rZviE9QI+oYMGSL7TaYCzje35azdS5YsUbEl4I+9BL4UJKvNfGDmxRdfVOx4o5Tk5GTV6rp8+TKaNGmiWn2kHAbhTnhDikK9Xo8PP/ywyB2QnLl//z7S09O1boZP4Y+/Y0Vt22jd3rNnz1r87UnOc2+g5omY1n2npJycnCL30KuiymAwIDo6Gk2bNtW6KaQwTkchG7dv38bcuXM9LseXRu/dofQPsnn5QrpCsqVGEHbjxg2kpqYqXo+1tLQ02cucNm2axd9F9WqC4ODBg1o3wSdkZ2fLWl5sbKys5WnN0fdk6dKluHz5sqSy7ty5gwEDBsjRLPJyDMKLOKlfbrFlLlq0SPZySTl//fWXallZiiKlA8kRI0aoPk0EgGLzW33JL7/8okm9vvZ9PHbsmKzl7du3T9byxNDihHLz5s2q3nxKRQuD8CLu3Xff1boJRVZSUpLWTfCIdfo9X74U7gm1tsvSpUsBeP5DL+X9RX2UWktKbDvz9IFanJQpSYkBH3NqfU+VrmfixImKli/YvHkz7ty5o2gdPL4oj0E4+aXMzEy8+uqrNssNBoPiB7aiKioqSusmSFZUT1CKWnuLilu3bilW9unTpy2+I+5mlRk/frxmU2hOnDiBGTNm2H3NFwKy69evY+3atYrWYTAYPM7+JWZb37lzB6dOnfKoHm/y5Zdfyj7lqShgEE4uDRo0SJX5e1rMqbV25swZvPbaa1o3w6kePXrg2rVrpr/V+nE8d+6cKvXITafTIS8vz/T3wIEDNWwNaUnJAGzLli0Wf7s7pW/v3r24fv26DC2S7uLFiz7zRM7MzEybZUpOEfriiy9kK8vVMV2tY76agwFRUVGKpab0ZsyO4oQ3pCi01qNHD0RERKha55kzZ5CUlIRGjRopWk+nTp00zwJTFEZ7rl+/jrS0NJsDpNwHzKKwLcSKiYkx/ds6+0dRYj266kt9RP/jqF95hUS8l156SZXfE6PRCJ1O5/T+DKX6Tc394e7duyhRogQqVKigWB3+eDxjEO6EN6QotKbVCImc/PGLpgT+IPsuR9+R3NxcRevNyMgQtayoKgrfGa3b6Gjf43HbvubNm8se7Hvjtv7iiy/QpEkTfP3111o3xadwOoqXioiIwLfffqt1M0zMfxiUSIvmi9Sa3ybcYKr1jzc9ZP4DumPHDsnv17IfrdMTApb7sS896t0bAx2BmJHwkSNHyl6vlvsej1//I3bfnDdvnsItIaUxCPdSCQkJuHjxoiZ1uzoYms+vVaJ8X6HkqKW9g7TS01G8OWjxVsuXL9e6CZJY93Hfvn0t/i7qN04VhWOP2Dbu3r1b4ZZY4vffewj7yN27dzVuiWsbNmzQLJYpChiEkw1vONgWFha6nV1AcPDgQYd3+hd169ev9/hkyJvYu4lKDt6wL8tJ7ROjmJgYn9uGArkfM75hwwZJ68+ePRvnz5+3+5o721w4Hpw/fx63b9+W/H6tbd26VbW6xGxfXzq+ysHd48Avv/yCQ4cOydwa38Eg3A+0bt3a4zJOnTplOrArMZp0+vRpi783btyI/v37e1RmdHQ0du3aJek9cn42JUfd/vrrL8XK1oKSIzpqjH7euHFD8ToA2Mw9LQoju97EfHupcXKxd+9ehylPV6xYITnrlLP+FoLG0aNHIyIiAunp6RZ5ywXHjx9HYmKipHp9jZibs9u0aSNbfVK/p/b2zU6dOsnVHLfbofTxxldP+J1hEO7l5NjpPR1RBoDFixdj27ZtAJT5Ig4ePNji7/T0dFkCM+FLLdfNJJcvX/abA4X159y2bZvH+W+14Ev9ZX2yqvZnK+pB/x9//GH6txrbbtSoUW7l/F64cKHHx7+vvvoKM2fOtFkeHh6OPXv2OHxfUfq+iJ2Sc/PmTYu/vf1ZEPb6wBtS+BalfaOoYBDuxbTa4b3lh9bTz6/T6UxlWI+IC4+ylhpUvvvuu0hJSfGoXUVVcnKy3ZE1TxmNRs0P7p07d1ak3MLCQkXKFaix3czr8JZjgxzU3nZS6PV6j59QmZ+fL/n4VtT6V+zNqW+88YbF35GRkUo0RzbekifcnBr7hiefqyjMj7eHQbgT4eHhGDp0qMV/Ss9by8nJwYABAxStwxVvSFElxxfeWRlHjhwB8DCXrFTubIdmzZpJfo+58ePHWzygxx65D5L79u2TtTytiNkuwolVQUEBVq9e7VF95vuH3CNu1vueGn0kJQg/e/as6ODxl19+KZJXVqQQckhL4en32Ly/7B2rxBy/7AU0chz/fb2/fU12djb+/PNPAJ71v6P90FGZ//73vyXX0bVr1yKZuY1BuBMzZszAvHnzLP5T+kE9ubm5XvtkQrlGD8R+mdUI+q2foObsgQtSXb161WbZ5cuX3TpQbN68GfHx8U7XkTsIt7ctlBoN0XokXJCdnW03TZ+7AgICcPLkSfz999+ylKf1dnLV/zNnzsSKFStElbVw4UKHU+V27Nih6NMN1aJlfznrK+vXrO9p6Nq1q6TyxGrVqpXHZVizF9i7utFbjX65d++e3eVxcXEuB7o8aZ87VysnTJhgd3l8fDx++OEH0986nQ7//POP3XXT09Ml1blw4UJ88803FsuEz+1P9yz4fBCenZ2N5cuXY+zYsejVqxdee+011VM7SSHshN5yWXDv3r2mfwtfMqUyWZgzn0oi1oQJE3DmzBno9Xq358+NGjXKrffZc/jwYZtlQ4cOVezR2XI/8leNNIjO6lKrXFcjh54ICAjA8uXLJacqlNqONWvWKJYSU2xbcnJyLAKAL7/80u06v/nmG594MBkAm6cuuyJ8x6R+/+z9dojtuzfffBMAMGvWLNMy65Mgd74bRqNR1oENe8aPH2+zTOl89hs3bsTvv//u8PVVq1bh119/tfta7969FU1f265dO8nv2bRpk93lw4cPN/07NjYWERERGDJkiN11O3ToIKmv7927h1u3bklrqA/y+SA8PT0df/75J27evInHH39c6+ZIpkUwbl7nF198YfP60aNHba4ItG3bFgaDwWZdd0ezXH1ug8Fg86OwadMmXLt2Db/99hveeust/PrrrxbrrF69ukifYbsKsqVOoxg4cKDT19UKwqX8uLszrcdVm48dOya5zJMnT2LSpEkAxF09kfIZ3dnG33//vSzBjif9O2/ePIurP1FRUaZ/3717127QYm+7fP/99263wdu4Mx3FXdYjkZ7W+/zzz1v87U4QfuvWLXTp0sXtNly/fh0xMTFO1/F0v8/NzUVERISodS9cuIAbN27gn3/+MU1ptGfmzJmyT99wJTw83PRvreZH27sHRtgPO3bsqHZzigSfD8IrVqyIJUuW4Pfff9d8rrUUzkZDlCamTutLbbm5uZLu6BbmmblryJAhDs/eDQaD3Utj06ZNQ3R0tEf1Ap73ibvvHzNmjNPXpT49zVWaLqUfkW5+0qbUfn7hwgWHuZgFwoni4MGDRc/h/vDDD015oa3zQ7uaQ3337l1RI+N5eXk2feDN01GEttpbJy4uTvT+uWbNGoty0tLSnJ7MT5482fRvOTJBycmd/lJ6Tri7UlJSsGDBAknvsTcwI8XixYtdTg+z/oxz5swRXX5MTAySkpIwdepUUetPmTIFK1eulCVf/4IFC2yuKrvbXwUFBdi/f7/pbzVzrovlaBqm0Wj0+J6poszng/DAwEBUqFBB62Z4hevXr6NPnz6K1zNr1iynP7jOfmTMz+adHZDi4+Pdmlttb2RfKk9+2LxlmpErBoPB5QlLbm6uR/N2W7Zs6fZ7xXJ0k6DBYLB5kNPp06dlf4CLTqez6fOLFy/aTR1nbfr06Rg7dqys7fGUu/uvebCRnp6O2bNnA3D+XRLqGjx4MBYvXuxwvfXr15v+3bp1a6+6OUvNkXBrctd79epVh1Ms1GoDAPTq1cvp6+ZpKB0R9jvrJ8IKnA1AWJ8kAvZHgLds2WL6t73t8OuvvzqcN+7oexEdHY09e/bYlCfH1CGp5JguFh0djebNm5v+lvuk8cMPP5StPKX4fBDuCb1eL+o/OUdfzEc7xR7AHR0wUlJSTA/YSU9Pxy+//ILY2FhTYGJ+kDDn7oFT+AJdvXrVIpOH0Wi0e5OiPcLZvJg22FvHei55SkqKx6Mxaho1apTpIC+3nJwcUT9QwMMA2dV0lJkzZ2L06NF239+sWTNJozFKpD50Jj8/33QDodiRw40bN9oss/7xc5SS0J0fl/T0dGRkZHhcjqfMP7dOp3PrnhDz49rdu3exdOlSu+vZO5bl5eVJyqpRlL7vzhiNRuzevdt0TBR7cqH2A4mk8jRt56VLlyz+VuIzWk+dSE9PR0pKiukkXent6qj8HTt24P/+7/8cBu9i/fe///Xo/QDQo0cPj96fnZ3tcRusmU9dSktLs3m4mTcqrnUDvJnY6Su9evVC7969ZanT3s18ruzYsQPdunWzWW6e+7iwsND043Tp0iU0aNAA48aNw4IFCywupycmJlpc3jXn6ovv7MD0zjvv4OTJk5g8ebLNQdSe+/fvOz25kXIQlHtus5IH4JiYGFSrVk2RslNTUzFnzhz85z//cbsM8+2WmZnpNCDbu3ev6GxC8+fPd7tN7nAUeDvr24kTJ6J79+4Wy86fP2/xHuvviHAi7U6ubXsn4fbaJ1zKlWO/tFeG9eO7X3rpJSxfvhxPPPGEw8/iblumT5/u1vvMSflunzp1Cs8++6zHdTrjybFm3bp1KFu2LF588UXJc2o9HUxxJicnByVLlnSr/Pz8fJQoUcLlesJ+p8Tx1lWZ1gHijBkzcPv2bcUfZS+0y9E8d3evPN64cQMVK1ZE6dKlATw83jq6wdIeT6aL3L9/X3TGJClOnz4NnU6Hp556yrRs2LBhpsC7qFx15ki4E4sWLcKqVatc/vf2229r2k6xKdXs7ZTWdycnJCSYDgBSR5TEHCzNLx07s3TpUpdBuCcBgKt0f82aNfN4dFZo3/Hjxy2W375922ZUa8uWLab6lMomoNWo2K1bt5zeKBQeHq766KWwLbZt2ybrHP+kpCSPyrIu9+jRoxb3T7jTVqmZWZwR9uk+ffrgzJkzDtdxd5uaf+fc/RGV8uAlR1dyBEePHnWrDeY8DQY8eeS5o/zMQpmZmZmiA0vz77D1TZuubN682XTCLvbz3Lx5E5mZmS5T9tr7jFOmTLFZZp4xRer+aZ3TWuljqaOrBe5m1urfv7/NvStycDY1yDxwF67oSbnpev/+/Rbz3K39/vvvWLZsmcPXGYT7gODgYFH/BQYGyl63lBszlcqn+9prr9ldbp43NDk52ZTGUOmDlPmXWo25lgkJCXaXi/1s27ZtA/AwLaG5TZs24f3337cob9y4cabUUn/99Zc7zXVJaPfSpUtNwYXUqQXm29zViYzgu+++w08//SS6XDUIP3IxMTEW/Snm4TfmV3KsA86BAwda/O3OpXfr/cud7C3mxMw/d8fgwYM9LuPUqVMWf0u5YiDHcc/Vd3nYsGE2y+ydUDobKXT3YT3mbfvss88klyGm3mHDhmHhwoUWyxxlYXIn7Z/QhvHjx5uyMUnZHtb5y4Xvp6u+t3fy9PLLL4uu1xHztiuZKUrq76f1iYrYuMHdVL4CMVe1zZlPtbS3/czbvXbtWrtTAKVatWqVx2UoiUG4F3M3L6u9NHyuRnwAWNwgAYhLmL9s2TKbmx09GQlTwldffaV1E2wIc1wvXbqEF198UZE6Zs6ciYsXL5r+Fvpk9uzZ2LVrF4CHUwukBIorV640jVpfvHjRrX6+fPmyzbQOb7l5Tbh3wTo4NGd+pcJVu1NSUjyajiJmmT0FBQUe/8DaY9128ytp7sxFtp4bLmXqmNTRWE+cPn3a9O+uXbuisLBQ9Gdct26daUTvwYMHkq76CIMNBw4ckNReo9GIkydPOrxaIcjMzFRkbq49wndL7P5vb73PP/8cABxmM5Gyz7vzAC1HQbjSv3dXrlxx+vonn3xi8bezY4eQOnT58uUeTU1UmtFoFJVj//Llyzb3jGRkZJiuNssxxU1JDMJ9xMCBAxEZGYm7d+/afeTrkSNHXB78PB1dtvewCKU4O+hZv+bskpYzV65cMd28+sUXX5gecS/nj1ZGRoZiP4LLly9HbGys3dfMp8g4SlVobxvPmzfPrRuRd+7cafp3amqqw6sMwMMrAfZGsoQRDbGXz80DJ2vCjUmOThjFPvhIp9O5PFnV6XSS7vWwHgkV+wNvPlK5e/duu089NLdlyxabbBeuvrvWr+/YscPues7aLFwhAmyvNjl6n6OUpsKVHHeDILHvsx71N09R6cq5c+dMQfTLL79skT/dGZ1Oh8LCQo8ee+8qX7Qcx2pXebztEXPir9PpLIIw86sN5sdM8z40TwjgjNFoxAcffGCqx5FmzZpZHAs8SR08evRoi4ffifHpp59i165dkp8BYU2n0yEmJgbXr183ZS4rKCjQNKWnve2ekpKCF154AYD9bVxQUIBPP/3UYtm7775r8zv3008/mfrX2zEIL0KcPV3q4sWLSEpKEnVwU+qsfeTIkZLfc/fuXbfb4+4PiDAKLHjnnXfsrnfgwAGMGzcOa9aswd69e00/+O4G9cD/2qz1fDXzfNiOMk+Iyaku9Ql+5ssdWblyJYYNG2azLwsjGocOHTItS01NdZgdyN68UIEwxUPsVRvhhsu5c+daLB8yZIjoS6bCVQnz7XH9+nWHI3JSvxf9+/c3/Ts/P99lnvczZ85gwYIF+PbbbyXVY87eFLSYmBgcPHjQ4Xus75Ew52jqjaMTVeHE2HpbmU/3kot5Srb4+HikpKS4VY7UwMeTOeFS1nc3IHOU5s8RnU6HFi1aiFrP0Wffs2eP3SsKYvvd3pUFR/cAOfoeSe2XI0eOmPYhYZu7ms5x/vx5tx68Y28f2LFjBxYtWmRxIqH075CjbTps2DC76WP1er3pu26vffn5+aYBDVcDcUUlS5JfBOFbtmzBqlWrTMHX8ePHTTdVqp0aTSwhODDf0V5//XWH61v/6Np7sIvcXzjrEUkhODKvJyoqypQT2FpERAS6du3q9qNrdTodFi1ahJ07d5oCGTFBlflIHACH6ROFwMt6FNP6x2rlypW4d++eTbnWMjIyTFMZhG20efNmp+9xl/k0FODh3EqpKaUcZcNxdnCzDujT09NN22v58uXYu3evy/4RXnd0cmD+OPTevXtj3LhxTstzVoezf5sTsg+JuZHTUTBj7zu5fv16hycL9kbClTiBjoyMFL2umFSBMTExovPxW9+T4Op4fODAAVGjicJ0IlftTU1NtbuOvX2/R48ebk/xMRqNWLdunaT37N+/X5EBCkfBlxIP53L2BEVPnThxAsDD41HPnj0dngy4emKn0B4xUzMcjYRbpxJ19l5zju7JEcpPT0+XNVhWah67o/3UXkpcg8GAo0ePIi4uzuY16+eE6HQ6i6kl1vU6GtTSepBLCr8IwtevX49ly5aZbng7cuQIli1bhmXLlrmV81YtUnek1NRU03vsBYRy75ht2rRx+JpwA43wcAF7hHl91l9g8x/FnJwci1EwIWDOysrCsWPHMHfuXCxdutR0WfzOnTuiL0mae/Dggc38SeFv6/ZZ/1j9+OOPeP311zFmzBinNwzZe3Kjo9FmvV5v96mfrgiXhvv16wfgf21PSUmRlGvZGfPRG6PRiIKCAixYsACxsbFo27atxbr9+/c3BURbt27F8ePH7R6w//nnH9O/zU8gXKXGSkpKsrmyIYXYkXApoyrW01OE6QfC9896tFfsScns2bNFB2Wu9lmDwSA6U5E568wj7sx1N+cs6LZX3saNG+1OA3G0XVq1auW0/oKCAixfvhzp6enYuXMnZs+ejdzcXISFhdkNrKxPLvR6vahBBKPRiO+++86irQaDAdOnT8dnn31mkVVH+NzC53Rnu1rfHCxM75o+fbpqUxCSkpLsnvhL+Tz2AjWB8BnPnj2LK1euOJwW4yrblHDzp/BMDWv22nvgwAGL5Z06dXJah8CdDCXuBuHmN6KaB6vCFB+5TujXr1+Pjz76yO5r9j6v9QCROfPf+suXL2PPnj344YcfTFfZrefhF6Vg2xG/yBNuffe3WOHh4TY3BoSFhYnOfeyJO3fuICUlRdQZtljCDvv999/bzBtXav6UmCDH+otkfrl68uTJKCgoQIMGDQD8L984YDnyLtSxYMECtw4uUVFRmDRpkt3XrLMCmAc0QvAvBGnOUpqZn82bX3KzZ+bMmbhx4wZeeOEFlyM55vr27ev0cr89nh6MhUDnzp07FsGqTqezCLLi4uLQqFEj0eUKd/wXFBSgeHHnh6q7d++ievXqSE5ORqVKlQCI+1zO0lzayyEt5qDv6iYqYW57Tk6OxfZKSkpCfHy8zQj44cOHkZmZiaVLl+LNN990Wb85YSSybdu2Fg+usO4nd3n64BVnHGVPcNRed592+vPPP2Pr1q2mE3zhuCLcBGjO+vOOGDHCrToFBQUFOHDgABYvXmz36oHUOeGDBw+2WT87O9uUP1kYWBAGKpQMYjwdWRdzD5McxA7E/fPPPyhWrBgA8feMCAoLC5Gfn286YXM0gmvv6tevv/7q9Cq4I46y2ZiP5ou58dGV+/fvmzJlCe0WTnzsXTmKiIgQVa7wfTYYDA6PMwzCfdyMGTMQHBysSd3bt28H8L/0YmIuwS5dulRUzvKsrCzTXErg4UHInTvFxQgICDB9Mc2nEThjPgcuOzsb+fn5+OWXX5y+x3z+sauDs735gI4CcMBx+qfjx4/jm2++sanfEfPgTLjcbT0VRhiVEW6amTVrluQHRNi7g99ZsCT1alBCQgKefvppm+XWc6ONRqNN2Zs3bxadx9n8xNBVerSuXbti7NixGD9+PLZv3+5yiorQF472F4PBgPfff9/miWvuzM8UWPf1sGHDcPr0adSpUwfZ2dl47733LEZEzdslfH6pgYej1JBa/ng5G910xVFQajQabQYspPSVeVAl9JO9eejmxwLhmNS4cWOX5bu6YmA0GjFnzhyULVvW4rsqdbTv9OnTqF27tsu616xZY3OSIcdV4Rs3bmD+/Plo3769wxF3sZ/HVWo56xv03CWMsoo5aXD35HXFihUW21fok2nTpiEoKMhhOmDA/bn69phPf5RzJNl8HzMYDE4Hjty5x8EXgm1HGIQXEeajJLdu3ULNmjXtfjkdZboAHk7DEZgfEMwDcjly/1o/cEP4UXGUFcA6Q4Jw2VbgKvAICAiQnMLLE0J7zPN/iwnC3b3BSribXaz33nvPpgxn2/DLL7+U9HjfsWPH4pVXXgFg/zOZB1n2TiDcuXHM+jHLaWlpKF++vMUy4XJySkoKjh49inr16rks2zoIF9M2uUbgzLfTsGHDbB70Y16PMKVHTJaV3Nxc002D1lNjNmzY4HBU7d69e5JvOlZyJNwRe/tcYWEhJkyYYLHMVXYYV2UKXD0gRcw2MA/ejEYjvv76a5t1Tp8+bZPR58yZM6hZs6bL8q2JSZ9nPT1Nyvayp1u3bqbv4KVLl1CrVi2PynN2c69YUk8sIiIicOHCBY/rtWY9tVA4ATIYDLh06ZLpd9toNGLPnj02o8We3pzrbBDDWRYpRxw9VGvJkiUuT2bsze92Ji8vz+bzT58+HcnJyShTpoyksryRX8wJL8rsHUBff/11REdH230SnrNMB2rMf8/Pz7fIsbxlyxaXTxF0NtJr7+BjPcdPp9OZLl2pERQ4q8P8RMcTnuRatzfnTkxZwnaVchnZXrmObsQVSL257d69ezYp6iZPnmzTD+YZX8RatGiRRfovNYNKYfQ1NTXV7o+kecAt3ITm6GZZcxEREZgzZw4Ay/7R6/V2r/gI37EDBw5IvlFdqUxLwMN23b171+Ik3dEl9MjISKcDEGLqsiYElPZunDWf+uLOPQX27mNwNMVB6sO7bt26hcWLF9u07+DBg7KOKFoH8Y7mVLtLjvLMb/QTY+rUqTY3yzt7MmzLli1FlWsdhFtPGRQeZAQ8/A0Rvu/mdYkhZeqtMMJsnaVEzEPYhAfLmZcFPLwS4upKl9QBs6ioKJvPv2rVKly+fNmr7+kTi0G4l3M0z3H37t2SR0iVYP2wjdatW0suw9WPmPXon/UTKNW+VOWsvc5GzaS0MyUlRdYAR8xId5cuXXD69Gmb3NH2OPssrm78GTt2rMvyzdkbub1z547NfFwhp7tAbLBvfsOsEqNgrmRkZMi2D3/66acWwZF5H1hfTRB4Uvfvv/9umgYiZn+1l5bMkbNnz1p89zdt2oTc3Fy7J4lK5NoX+4RIqfce2MuSA3g2TccZ4cRyxIgRFgFWTk4Orl696vZx5ocffnB547Q95qkelebsoVtiOZqqqNfrRd+w7eymeDHbX+zTYYUpUmJOpoV56tbfybfeektUXQJ70x/lptPp7LZLmDpmfdIpdQqnlhiEe5Hvv//eZpmjs9Jly5Yp3RzVODvA2JvCYj1aqXY+0DVr1jjMiCDXQUiOS7HmFixY4PR14cf01q1bom46Ej6nvXm3cl0NcCYmJsbh1Axh/3Anu4zA3nQBgaOn9XnC1X7jKv2lwHqbmH83lBrldzZfvVmzZhb1vvvuu6LLnThxosWo94QJE2wCc7l48r2VYwqTklcUzLe/eXDWsWNHANKfeyDcr+TqCqcj3vyURnsKCgrsfnc8ufLiiKP5z1If375y5UqX60RGRiIxMdHujZKe3PuixL588eJFu7GQsK3MM7AdPXoUW7dulb0NSmEQ7kXE/tD6GuHmRrGsD4hajF668yjcJUuWSFpfjs81efJklw+EMLdnzx5RwZow5cTdHO9KEoIET34MnN0M5c60FzX9/PPPpn87mrtpfbPe5s2bPT7xc5QKrkOHDh6Va84bR7ikjO5rwdWN0FIv6QtTHoU0tFIVtSkEGRkZkqcESaHkIFJubq5NYgGBs2O3u/cHOHrYlqeEq23W28re1DRvffaLI7wx0wm1UxTa23nEXhIldZk/tdGcsxE160uS1jegKqVXr16i1z1w4ACee+45l+vZux/Bmrs/0p4Svkda3DToiqPL444eGOUpR9vA+sRs/PjxbtchBPSORlQ9Cboc3Z9gfYLlKAuMWNY3h/sTT64YOVLUAiEtWT+S3lHQ7A7r5zYowXzw0GAwKHpVx3ouui9kTWEQ7oSWKQrJMW8MrgRSHogj58FWTkqlq1SLOw+hUYvao6aO7gUwDzql5jy2x50HZIlhfcVJ+IG3nu8v5/MUlOBvQem9e/fwxBNPyFaeO3PPi6ILFy54/ZUVazdu3DA9mwFQdmqVNblvBtYCp6NQkaPml5xITr4YjPXv31/Vm+0A109B9DauMkJYj4YWNfYeUubt07a8kdQHrcmpWbNmbp/sSE3zqiRXzxTxNgzCqcjxlUvHUuZqk2+YO3euZnULl27lHlVUYjqDK71791a9TnJs2LBhNsuK2oguyUOpLD9iiUmx6E0YhBMROTB69OgiN+rqiPlDMuQerdJ69IvIFxTVK2XCw8FIOs4JJyJy4K+//kJQUJDWzZCF+UNIzHOjy8H6YUpEJJ2Y1ILkWzgSTkTkhNQcvUWB3D/2Yp7kSURElhiEExGRRxylEiQiIsc4HcUJtfOEExEVRb4yb56ISE0Mwp1gnnAiIiIiUgKnoxARERERqYxBOBERERGRyhiEExERERGpjEE4EREREZHKGIQTEREREamM2VGcYIpCIiIiIlICg3AnmKKQiIiIiJTg80F4fn4+li9fjqioKGRmZqJu3bro27cvnnnmGa2bRkRERER+yufnhM+aNQsbNmxAu3btMGTIEAQEBGD8+PG4cOGC1k0jIiIiIj/l00F4XFwc9u/fj/feew8DBw5Ely5dMHnyZFStWhV//PGH1s0jIiIiIj/l00H4oUOHEBAQgC5dupiWBQUFoVOnToiJiUFSUpKGrSMiIiIif+XTQfjVq1dRq1Ytm5srn3jiCQDAtWvXtGgWEREREfk5nw7CU1JSUKFCBZvlwrLk5GSn79fr9aL+y8/PV6T9RL6uatWqNsuKF1fmfvFJkyYpUi4AhIaGKla2tVKlSmHnzp2q1SeHMmXKOHytUaNGKFmyJKpVq6Zae/bs2aNY2T179lSsbAAoXbq06d9Vq1ZFgwYNFKnnpZdeMv172rRpitThax555BGLvytVqqRRS6io8OnsKHl5eQgMDLRZHhQUZHrdmQEDBoiqp1evXujdu7f0Blo5efKkx2UQkWPmU9OKOh4vPKPk9vvyyy8VK1sr3N+I5OfTQXhQUJDdUWoh+BaCcUcWLVokKk+4vUCfiIiIiMgRnw7CK1asaHfKSWpqKgDXl4qCg4P5sB4iIiIikp1Pzwl//PHHkZCQAL1eb7E8NjbW9DrJKz8/HytWrOA8eS/B/vAu7A/vwb7wLuwP78L+UIdPB+Ft27ZFYWEhtm/fblqWn5+P3bt3o1GjRqhSpYqGrfNN+fn5WLlyJb+4XoL94V3YH96DfeFd2B/ehf2hDp+ejtKoUSO0bdsWS5YswYMHD1CjRg1ERkYiMTERn376qdbNIyIiIiI/5dNBOACEh4dj2bJliIqKQmZmJurWrYsxY8aomlKMiIiIiMiczwfhQUFBGDhwIAYOHCj5veHh4QgIsJyxExYWhrCwMLmaR0RERER+yKfnhHtqxowZmDdvnuk/NQLwrVu3Funy1aLG5/CVOtTgK9uK/eEd5atVhxp8ZVuxP/yvDqX5ynbypA4G4RJ4e2d6Q/lq8YW+UKsONfjKtmJ/eEf5atWhBl/ZVuwP/6tDab6ynRiEExEREREVIQzCiYiIiIhUxiCciIiIiEhlPp8dxR1GoxEAbJ60WVhYaLNMbkrXoXT5QtlFfTv5Sh3sD++qw1f6g30hni9sKzXqYH94Vx1q9IcvbCd7dQj/FmJJZ3RGMWv5mfv372PAgAFaN4OIiIiIiqBFixahcuXKTtdhEG5HYWEhUlJSUKpUKeh0Oq2bQ0RERERFgNFoRHZ2NipWrGjzrBlrDMKJiIiIiFTGGzOJiIiIiFTGIJyIiIiISGUMwomIiIiIVMYUhQQAyM7Oxrp16xAXF4e4uDhkZmZi+PDh6Nixo931L1++jJUrVyI6Ohp5eXmoXr06XnnlFXTr1s1ivfz8fCxfvhxRUVHIzMxE3bp10bdvXzzzzDNurecPpPTF7du3sWzZMkRHRyMjIwNVqlRBu3bt8MYbb6BkyZIW67Iv3BMXF4fIyEicPXsWiYmJKFOmDBo1aoR+/fqhVq1aFutK2XbsD/eI7Q8p/QawP9whdRsLVq1ahWXLlqFOnTqYO3euzevsC/dI7Q/+jmuPI+EEAEhPT8eff/6Jmzdv4vHHH3e67qlTp/Dll1/iwYMH6NmzJ95//300b94cycnJNuvOmjULGzZsQLt27TBkyBAEBARg/PjxuHDhglvr+QOxfZGUlITw8HDExsYiLCwMQ4YMQUhICFasWIEffvjBZn32hXvWrl2Lw4cP46mnnsKQIUPQpUsXXLhwASNGjEB8fLzFulK2HfvDPWL7Q0q/AewPd0jdxsDDFMBr1qyxGSQwx75wj5T+4O+4lzASGY3GvLw8Y0pKitFoNBrj4uKMXbt2Ne7atctmvaysLGPfvn2NkydPNhoMBqdlxsbGGrt27Wpcu3ataVlubq5xyJAhxi+++ELyev5CbF+sWrXK2LVrV+P169ctls+YMcPYtWtXY0ZGhmkZ+8J90dHRxry8PItlCQkJxjfeeMP4ww8/mJZJ2XbsD/eJ7Q+x6xmN7A93SdnGgu+//944atQo48iRI41Dhw61eZ194T6x/cHfce/BkXACAAQGBqJChQou19u3bx/S0tLQr18/BAQEICcnB4WFhXbXPXToEAICAtClSxfTsqCgIHTq1AkxMTFISkqStJ6/ENsXwlO5ypcvb7G8QoUKCAgIQPHi/5ttxr5wX+PGjREYGGixrGbNmqhTpw5u3rxpWiZl27E/3Ce2P8SuB7A/3CVlGwPA+fPncejQIQwZMsRhmewL94ntD/6Oew8G4STJ6dOnERwcjOTkZHz44Yd4++230bNnT8ybNw95eXkW6169ehW1atVCcHCwxfInnngCAHDt2jVJ65Glpk2bAgB+/vlnXL16FUlJSThw4AD++usvdO3a1eJyL/tCXkajEWlpaShbtqxpmZRtx/6Ql73+kLIe+0M+jraxwWDA/Pnz0blzZ9StW9fh+9kX8rLXH/wd9x68MZMkuX37NgwGAyZNmoROnTrhvffew7lz57BlyxZkZWXhyy+/NK2bkpJid0RXWCbMPRO7Hll67rnn0LdvX6xevRrHjh0zLX/nnXfQr18/i3XZF/Lau3cvkpOT0adPH9MyKduO/SEve/0hZT32h3wcbePt27cjKSkJkyZNcvp+9oW87PUHf8e9B4NwkiQnJwe5ubl49dVX8cEHHwAA2rRpg4KCAmzfvh19+vRBzZo1AQB5eXk2l8aAh5eohNelrEe2qlatitDQULRp0wZlypTByZMnsWbNGlSoUAFdu3Y1rce+kM/Nmzfx66+/IiQkBB06dDAtl7Lt2B/ycdQfUtZjf8jD0TZOT0/H8uXL0bNnT5QrV85pGewL+TjqD/6Oew9ORyFJhC/Uiy++aLG8Xbt2AICYmBiLdfPz823KEL6MQlli1yNL+/fvx5w5c/DJJ5/glVdeQZs2bfDpp5+iQ4cO+OOPP5Cenm5al30hj9TUVEyYMAHBwcEYOXIkihUrZnpNyrZjf8jDWX9IWY/94Tln23jZsmUoXbq0xcCAI+wLebg6VgH8HfcGDMJJkooVKwKwvRlQGN3IzMy0WDc1NdWmDGFZpUqVJK1HlrZt24b69eujcuXKFstbtmyJ3NxcXL161bSMfeG5rKwsjBs3DllZWRg/frzNtpCy7dgfnnPVH1LWY394xtk2vn37Nnbs2IHXXnsNKSkpuHfvHu7du4f8/HwYDAbcu3cPGRkZpvXZF54Tc6wC+DvuDRiEkyQNGjQAYDu/KyUlBQAsbv54/PHHkZCQYMriIYiNjTW9LmU9spSWlmb3jvaCggIAD2+EErAvPJOXl4eJEyciISEBY8aMQZ06dWzWkbLt2B+eEdMfUtZjf7jP1TZOTk5GYWEhfvvtNwwePNj0X2xsLBISEjB48GD8+eefpvXZF54Rs8/zd9x7MAgnSZ5//nkAwK5duyyW79y5E8WKFTNl7ACAtm3borCwENu3bzcty8/Px+7du9GoUSNUqVJF0npkqWbNmrhy5QoSEhIslu/fvx8BAQEWGQjYF+4zGAyYNm0aYmJiMHLkSISEhNhdT8q2Y3+4T2x/iF0PYH+4S8w2rlOnDkaNGmXzX506dVClShWMGjUKnTp1Mq3PvnCf2H2ev+PegzdmkolwZ7Rwdnz8+HHTv7t27YpHHnkE9evXR6dOnbBr1y4YDAaEhobi3LlzOHToEN5++22LS06NGjVC27ZtsWTJEjx48AA1atRAZGQkEhMT8emnn0pez5+I6Ys333wTf//9N0aOHImwsDCUKVMGJ06cwN9//43OnTuzL2Ty+++/49ixY2jRogUyMjIQFRVl8Xr79u0BSNt27A/3ie0PsesB7A93idnG5cqVQ+vWrW3eu2nTJgCweY194T6x+zx/x72Hzmg0GrVuBHmHQYMGITEx0e5rCxYsQLVq1QA8nO6wZs0a7N69GykpKahSpQrCwsLQvXt3m/fl5eVh2bJl2Lt3LzIzM1G3bl307dsXzz77rFvr+QuxfREXF4cVK1bg6tWryMjIQLVq1dChQwe89dZbNjefsS/c8/XXX+P8+fMOX9+8ebPp31K2HfvDPWL7Q0q/AewPd0jdxtbvTU9Px9y5c21eY1+4R0p/8HfcOzAIJyIiIiJSGeeEExERERGpjEE4EREREZHKGIQTEREREamMQTgRERERkcoYhBMRERERqYxBOBERERGRyhiEExERERGpjEE4EREREZHKGIQTEREREamMQTgRERERkcoYhBMRERERqay41g0gIlLDa6+95vT10NBQTJkyRaXWkCfOnTuHUaNGWSxbtmwZypUr51G5wj6yefNmh+vs3r0bs2fPRocOHfDZZ5+5Xde7776LrKws09/Dhw9Hx44d3S6PiIoeBuFE5Fc6dOhgd3nt2rVVbgl5qkaNGmjcuDEAICgoSOPWSPPiiy8iNzcX165dw7Vr17RuDhFpgEE4EfkVT0Yvybs0bty4yPbn0KFDAQArVqxgEE7kpzgnnIiIiIhIZRwJJyKycu/ePQwePBihoaEYPXo0VqxYgSNHjiA5ORlhYWEYMmQIACApKQlr167F33//jeTkZJQoUQIhISF45513TNMkrB07dgxr1qzBtWvXUKJECYSGhqJ///7Yt28fVq5caTE3WJj77Gj+8cyZMxEZGYnvvvsOTZs2tXhNStvMP++4ceOwcuVK7N+/H6mpqahSpQo6d+6Mt956CzqdzqYNSUlJWLduHU6dOoX79+8jKCgI1atXR4sWLdC9e3cEBwfj0qVLCA8PR0hICKZPn253u6xevRpLly7Fu+++iz59+rjuJCc82W7uGjRoEBITE52u42yuORH5HwbhREQO5OXl4euvv0ZiYiJCQ0NRv359lC5dGgAQExOD8ePHIzMzE7Vq1UKzZs2Qnp6Of/75B6dOncIXX3yBF154waK8v/76C/PmzYNOp8OTTz6JihUrIjY2Fp9//jmaN28uW7vdaRsAFBQUYMyYMbh58yZCQ0ORm5uL8+fPY/HixcjOzka/fv0s1r9w4QImTpyIrKwsVK1aFc2bN0deXh5u3bqFFStWoGXLlqhXrx4aNmyI+vXrIyYmBvHx8XjssccsyjEajdi1axcCAgLQqVMn2baDmtq0aYP09HSb5WlpaTh16hQCAnjhmYgsMQgnInIgLi4OISEh+O9//2sKvgFAr9djypQp0Ov1+Pzzz/HSSy+ZXrt06RLGjBmDn3/+Gf/6179MGTsSExOxYMECFC9eHKNHj8azzz4L4GHgO3v2bOzdu1eWNrvTNkFMTAxCQ0OxYMECBAcHm97zxRdfYOPGjejRowdKlSoFAMjIyMCUKVOQlZWFAQMG4PXXX7cINGNiYlCxYkXT36+++irmzJmDnTt3mq4kCM6cOYO7d+/iueeeQ9WqVWXZDmobNGiQzTLhJA4A+vfvr3aTiMjLMQgnIr/iKFXhggULUK1aNZvl77//vkUADgC7du1CSkoKXn/9dYsgFwAaNmyInj17YuHChYiKisLrr79uek9eXh46dOhgCsABoHjx4hgyZAiOHDmC3Nxczz6cm20TBAQE4OOPPzYF4MJ7nnvuOZw4cQKXL182Td/YuXMnHjx4gGeffRZvvvmmTTtCQkIs/m7Xrh1+//13REVF4T//+Q8CAwNNr+3cuRMA8Morr3jy0WXjKp2lWD///DPi4uLQvn17u9uIiPwbg3Ai8iuOUhSWLFnSZlnFihXRsGFDm+X//PMPgIdTEOxp0qQJgIejyILo6GgAsDsNpGzZsnjmmWdw9OhRF613zZ22CapUqWI3VWOtWrVw4sQJpKSkmJadPn0aANClSxdR7SpZsiReeuklbNu2DYcOHTKdIDx48ABHjhxBhQoV0KJFC1FlKc3RPgIAd+7cwcWLF12WERERgb1796JRo0b45JNP5GweEfkIBuFE5FekpLSrUqWK3eXCDXhfffWV0/ebzxFOTk4GAIfTLeSahuFO2wSVK1e2u64wBSU/P9+07P79+wAe5uoW69VXX8W2bduwc+dOUxAeFRWFgoICdOzYEcWKFRNdlpKc7SO7d+92GYQfP34cS5cuReXKlfHNN99YjPoTEQkYhBMROeAoeCosLAQAtG3bFiVKlHD4fqUfAGQ0Gm2WedI2e9lP5FS3bl00btwY586dw+3bt1GzZk3s3LkTOp0OnTt3VrRutdy4cQM//PADihcvjm+++QYVKlTQuklE5KUYhBMRSVS5cmUkJCSgR48eaNCggaj3VKxYEQkJCUhMTESdOnVsXk9KSrJZVrz4w0N0Tk6O3TKF0WhP2+aOypUr49atW7hz5w7q1q0r+n1dunTBxYsXsXPnTrRo0QI3b97E008/jerVq8veRvORe3P2rgLIIT09HRMnTkR2dja++uorRbc/ERV9zJlERCTR008/DQA4cuSI6Pc8+eSTAICDBw/avJaRkWGay21OGEW9ffu23fdcuXJFlra5Q6hnx44dkt73/PPPo0yZMtizZw+2bdsGAIqNgtvbbgaDAfHx8bLXZTAY8P333+Pu3bvo2bOn3bn/RETmGIQTEUnUpUsXlC9fHuvWrcP27dtNU0AEBoMBp06dsgj2OnbsiMDAQOzbt890UyPwMEXhggUL7I52V69eHVWqVMH169ctbtrMycnBnDlzoNfrZWmbOzp37oyyZcvi77//xsaNG22mxsTExCAtLc3mfUFBQXj55ZeRlpaGffv2oVy5cmjVqpVHbXHkypUrOHTokMWyiIgI01UHR1cY3DF//nycPXsWrVq18vhhQ0TkHzgdhYhIotKlS+Obb77BxIkTMXfuXKxatQqPPfYYSpcujdTUVFy5cgVZWVkYNWqU6cE01atXx6BBg/Drr79i7NixaNKkCSpUqICYmBhkZWXhpZdespsrvFevXvjpp58wZcoUhIaGomTJkoiLi0NwcDBatmyJY8eOedw2d5QpUwYjR47ExIkTsWDBAmzevBkNGzZEXl4ebt68iTt37mD27NkoX768zXu7dOliCtw7dOig2I2LVapUwffff4+QkBBUrFgRN27cwM2bN9GgQQNcvnwZCxcuxKlTp/DBBx94VE9SUhL++usvAA/TPM6aNcvuelJuCiYi38cgnIjIDSEhIZgzZw42btyIEydO4Pz58wAezv0ODQ1F69atTVM2BGFhYahYsSLWrl2L2NhYBAUFoUmTJujfvz8OHDhgt55OnTpBp9Nhw4YNiI6ORunSpdGiRQv0798fCxculK1t7mjatCl++uknrFu3Dn///TeOHj2KUqVKoVq1aujTp4/Ded61atVCpUqVcP/+fUVvyGzatCmaNGmCtWvX4tKlS6hSpQo+/vhjtGzZEmPGjMG9e/fspqaUyvxqw+HDhx2uxyCciMzpjPZuryciIlWtWLECK1euxPDhw9GxY0etm6OomJgYfPnllwgNDcWUKVMkv//cuXMYNWoUOnToYDewdfW6N/GnficiSxwJJyIiVa1atQoA0LVrV4/KuXjxImbOnAkA+PDDD035zIuCefPmITc3F9euXdO6KUSkEQbhRESkuIsXL2LXrl2Ij49HXFwc6tevj9atW3tU5p07d3Dnzh0AwMCBA4tUEL5//35kZWVp3Qwi0hCDcCIiUlxCQgJ27dqFUqVKoVmzZvjoo48QEOBegq6mTZti8+bNMrdQXX/++afWTSAijXFOOBERERGRypgnnIiIiIhIZQzCiYiIiIhUxiCciIiIiEhlDMKJiIiIiFTGIJyIiIiISGUMwomIiIiIVMYgnIiIiIhIZQzCiYiIiIhUxiCciIiIiEhlDMKJiIiIiFT2/wD36YJoPu29vgAAAABJRU5ErkJggg==", "text/plain": [ "
\n", "\n" ], "text/plain": [ "