From 203ded171bb40841428e75c93d71fe0909f2d28d Mon Sep 17 00:00:00 2001 From: RondeauG Date: Wed, 20 Nov 2024 11:52:14 -0500 Subject: [PATCH 1/6] add agent to pooch --- docs/notebooks/climate_change.ipynb | 11 +- docs/notebooks/hydrological_modelling.ipynb | 5 +- docs/notebooks/pmp.ipynb | 126 +++++--------------- 3 files changed, 37 insertions(+), 105 deletions(-) diff --git a/docs/notebooks/climate_change.ipynb b/docs/notebooks/climate_change.ipynb index 23b22587..e13b493a 100644 --- a/docs/notebooks/climate_change.ipynb +++ b/docs/notebooks/climate_change.ipynb @@ -28,17 +28,20 @@ "from xhydro.testing.helpers import deveraux\n", "\n", "D = deveraux()\n", + "downloader = pooch.HTTPDownloader(headers={\"User-Agent\": f\"xHydro-{xh.__version__}\"})\n", "\n", "# Future streamflow file (1 file - Hydrotel driven by BCC-CSM-1.1(m))\n", "streamflow_file = D.fetch(\n", - " \"cc_indicators/streamflow_BCC-CSM1.1-m_rcp45.nc\",\n", + " \"cc_indicators/streamflow_BCC-CSM1.1-m_rcp45.nc\", downloader=downloader\n", ")\n", "\n", "# Reference mean annual streamflow (QMOYAN) for 6 calibrations of Hydrotel\n", - "reference_files = D.fetch(\"cc_indicators/reference.zip\", pooch.Unzip())\n", + "reference_files = D.fetch(\n", + " \"cc_indicators/reference.zip\", pooch.Unzip(), downloader=downloader\n", + ")\n", "\n", "# Future deltas of QMOYAN (63 simulations x 6 calibrations of Hydrotel)\n", - "deltas_files = D.fetch(\"cc_indicators/deltas.zip\", pooch.Unzip())" + "deltas_files = D.fetch(\"cc_indicators/deltas.zip\", pooch.Unzip(), downloader=downloader)" ] }, { @@ -651,7 +654,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.12.3" + "version": "3.12.7" } }, "nbformat": 4, diff --git a/docs/notebooks/hydrological_modelling.ipynb b/docs/notebooks/hydrological_modelling.ipynb index 752796fc..bfb10c4d 100644 --- a/docs/notebooks/hydrological_modelling.ipynb +++ b/docs/notebooks/hydrological_modelling.ipynb @@ -111,11 +111,12 @@ "from xhydro.testing.helpers import deveraux\n", "\n", "D = deveraux()\n", + "downloader = pooch.HTTPDownloader(headers={\"User-Agent\": f\"xHydro-{xh.__version__}\"})\n", "\n", "# This notebook will use ERA5 data for a small watershed in Eastern Quebec, along with faked elevation data.\n", "\n", "# Streamflow file (1 file - Hydrotel driven by BCC-CSM-1.1(m))\n", - "meteo_file = D.fetch(\"hydro_modelling/ERA5_testdata.nc\")\n", + "meteo_file = D.fetch(\"hydro_modelling/ERA5_testdata.nc\", downloader=downloader)\n", "ds = xr.open_dataset(meteo_file)\n", "ds" ] @@ -509,7 +510,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.12.3" + "version": "3.12.7" } }, "nbformat": 4, diff --git a/docs/notebooks/pmp.ipynb b/docs/notebooks/pmp.ipynb index 6e13e7a4..5ecde174 100644 --- a/docs/notebooks/pmp.ipynb +++ b/docs/notebooks/pmp.ipynb @@ -16,7 +16,7 @@ }, { "cell_type": "code", - "execution_count": 1, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -42,36 +42,27 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ + "downloader = pooch.HTTPDownloader(headers={\"User-Agent\": f\"xHydro-{xh.__version__}\"})\n", "path_day_zip = deveraux().fetch(\n", - " \"pmp/CMIP.CCCma.CanESM5.historical.r1i1p1f1.day.gn.zarr.zip\"\n", + " \"pmp/CMIP.CCCma.CanESM5.historical.r1i1p1f1.day.gn.zarr.zip\",\n", + " pooch.Unzip(),\n", + " downloader=downloader,\n", ")\n", "path_fx_zip = deveraux().fetch(\n", - " \"pmp/CMIP.CCCma.CanESM5.historical.r1i1p1f1.fx.gn.zarr.zip\"\n", + " \"pmp/CMIP.CCCma.CanESM5.historical.r1i1p1f1.fx.gn.zarr.zip\",\n", + " pooch.Unzip(),\n", + " downloader=downloader,\n", ")\n", "\n", - "directory_to_extract_to_day = Path(\n", - " path_day_zip\n", - ").parent # To extract to the same directory as the zip file\n", - "with ZipFile(path_day_zip, \"r\") as zip_ref:\n", - " zip_ref.extractall(directory_to_extract_to_day)\n", - "\n", - "directory_to_extract_to_fx = Path(\n", - " path_fx_zip\n", - ").parent # To extract to the same directory as the zip file\n", - "with ZipFile(path_fx_zip, \"r\") as zip_ref:\n", - " zip_ref.extractall(directory_to_extract_to_fx)\n", - "\n", "path_day_zarr = (\n", - " str(directory_to_extract_to_day)\n", - " + \"/CMIP.CCCma.CanESM5.historical.r1i1p1f1.day.gn.zarr\"\n", + " Path(path_day_zip).parent / \"CMIP.CCCma.CanESM5.historical.r1i1p1f1.day.gn.zarr\"\n", ")\n", "path_fx_zarr = (\n", - " str(directory_to_extract_to_fx)\n", - " + \"/CMIP.CCCma.CanESM5.historical.r1i1p1f1.fx.gn.zarr\"\n", + " Path(path_fx_zip).parent / \"CMIP.CCCma.CanESM5.historical.r1i1p1f1.fx.gn.zarr\"\n", ")" ] }, @@ -93,7 +84,7 @@ }, { "cell_type": "code", - "execution_count": 3, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -121,7 +112,7 @@ }, { "cell_type": "code", - "execution_count": 4, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -139,7 +130,7 @@ }, { "cell_type": "code", - "execution_count": 5, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -164,7 +155,7 @@ }, { "cell_type": "code", - "execution_count": 6, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -182,7 +173,7 @@ }, { "cell_type": "code", - "execution_count": 7, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -203,7 +194,7 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -232,7 +223,7 @@ }, { "cell_type": "code", - "execution_count": 9, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -242,60 +233,18 @@ }, { "cell_type": "code", - "execution_count": 10, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "" - ] - }, - "execution_count": 10, - "metadata": {}, - "output_type": "execute_result" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAioAAAHFCAYAAADcytJ5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAABnkUlEQVR4nO3deVyN6f8/8NdpTxtlqZRspaGyhIiZouxjncUQEmMZjOyNbcQgwzB2YxsyxjRjN4x1RtliUpYY6yhiSh9Lkmg55/r94dv5OSqd0zl1TvV6zuN+POq6r/u63/eZ5O3abokQQoCIiIhIB+lpOwAiIiKiwjBRISIiIp3FRIWIiIh0FhMVIiIi0llMVIiIiEhnMVEhIiIincVEhYiIiHQWExUiIiLSWUxUiIiISGcxUXlLaGgoJBIJHj16pLE2Bw8ejNq1a5d4PPPnz8eePXuKdR91pKenY968efD19YWtrS3Mzc3h7u6Ob7/9Fq9evVK6nYiICDRp0gQmJiawt7fHuHHjkJGRUeR1mzdvhkQikR+a/H+3Z88ehbbPnz+vsbaLkvf/XpPyPqvExESNtqtNV69exahRo9C6dWuYmZlBIpEgMjJS22GVim3btmHp0qUFnpNIJAgNDZV/HxkZqfOfzbFjx9C6dWtUqlQJVatWxeDBg5GamqrUtc+fP8fYsWNRs2ZNGBsbw8XFBQsXLoRUKi3hqKmkMVEpBTNnzsTu3btL/D7aSlTu3buHpUuXolmzZli3bh327duHjz/+GKGhofjwww+hzFsafv75Z/Tr1w8tWrTAwYMHMWvWLGzevBl9+vRROo5du3YhOjoalStXVuNpFPn4+CA6OhozZszQWJvK+vzzzxEdHV3q9y1rzp8/jz179sDa2hp+fn7aDqdUvStRiY6Oxueff166AakhKioKXbp0QY0aNbB3714sW7YMx44dg5+fH7Kyst55bW5uLjp06ICtW7di2rRp2L9/P7p3746vvvoK48ePL6UnoJJioO0AKoJ69eppO4QSVadOHSQmJsLMzExe1r59e5iZmWHy5Mk4ffo02rZtW+j1UqkUkydPRseOHbF+/XoAQLt27WBhYYGAgAAcPHgQXbp0KTKOpk2bFrvnqjBVqlRBq1atcP36dY22qwwHBwc4ODiU+n3LmoEDByIwMBAAsGPHDvz+++9ajkg3tGrVStshqGTy5MlwcXHBjh07YGDw+q+mOnXqoE2bNvjxxx/xxRdfFHrtjh07cO7cOezcuVP+j5sOHTogIyMDq1atwujRo9GgQYNSeQ7SPPaoFOLhw4fo168frKysUKNGDQwZMgTPnj1TqCOEwOrVq9GkSROYmpqiSpUq+Pjjj3Hnzh2FegUN/aSlpWHo0KGwtraGubk5unXrhjt37uTrrlU2HolEghcvXiA8PFw+TOHr66upj+OdzMzMFJKUPC1btgQAJCUlvfP6s2fPIjk5GUFBQQrln3zyCczNzdXqjfL19YWbmxuio6Ph7e0NU1NT1K5dG5s2bQIAHDhwAM2aNUOlSpXg7u6OQ4cOFfteBRFCoEaNGhg9erS8TCqVokqVKtDT08PDhw/l5UuWLIGBgQHS0tIAFDz0U7t2bXz44Yc4dOgQmjVrBlNTU7i6uuLHH3/Md++zZ8+iTZs28qG0qVOnIicnJ189mUyGhQsXwtXVFcbGxqhevToGDRqE+/fvy+usWrUKenp6Ct3wixcvhkQiUXg2mUyGKlWqYOLEiap/WMWkp1cyv8YOHDiAJk2awNjYGHXq1MF3332X7/9JYmIiJBIJNm/enO/6t/8s3759G0FBQXB2dkalSpVQs2ZNdO/eHfHx8QrX5Q3R/PLLL5g+fTrs7e1haWkJf39/3LhxQ17P19cXBw4cwN27dxWGJwu7f2HOnz+PHj16wNraGiYmJmjatCl+++035T8oDXjw4AFiYmIwcOBAeZICAN7e3nBxcSnyd8Dp06chkUjy/YPmww8/hEwmK5UebSo5TFQK8dFHH8HFxQU7d+7EV199hW3btuXrQhwxYgTGjRsHf39/7NmzB6tXr8bVq1fh7e2t8BfQ22QyGbp3745t27YhJCQEu3fvhpeXFzp37lzseKKjo2FqaoquXbsiOjoa0dHRWL169TufUSqVIjc3t8hDJpMp+akp+uuvvwAAjRo1eme9K1euAAA8PDwUyg0NDeHq6io/X1wpKSkICgrC559/jr1798Ld3R1DhgzBnDlzMHXqVEyZMgU7d+6Eubk5evXqhf/++0+t+71JIpGgffv2OHbsmLzs/PnzSEtLg4mJCf788095+bFjx+Dp6Vnk0NWlS5cwceJEjB8/Hnv37oWHhweGDh2KEydOyOv8888/8PPzQ1paGjZv3owffvgBFy5cwNy5c/O198UXXyAkJAQdOnTAvn378M033+DQoUPw9vaWz/fx9/eHECJfvKampjh69Gi+Z/P393/nM8hkMqV+9rQ1v+DPP/9Ez549YWFhgYiICCxatAi//fabPMEtjv/++w82NjZYsGABDh06hFWrVsHAwABeXl4KCUieadOm4e7du9iwYQPWrVuHW7duoXv37vLPZPXq1WjTpg1sbW3lf+ZVHSo8fvw42rRpg7S0NPzwww/Yu3cvmjRpgr59+xaYfL1NU79DCvsdkFdW1O+A7Oxs6OnpwdDQUKHc2NgYAHD58uUin4V0mCAFs2bNEgDEwoULFcpHjRolTExMhEwmE0IIER0dLQCIxYsXK9RLSkoSpqamYsqUKfKywMBA4eTkJP/+wIEDAoBYs2aNwrVhYWECgJg1a5bK8QghhJmZmQgMDFT6WZ2cnASAIo8341HWpUuXhKmpqejdu3eRdefNmycAiOTk5HznOnbsKFxcXN55/aZNmwQAkZCQkO+cj4+PACDOnz8vL3v8+LHQ19cXpqam4sGDB/LyixcvCgBi+fLlhd4jJiamyOd524YNGwQAce/ePSGEEHPnzhWurq6iR48eIigoSAghRHZ2tjAzMxPTpk2TX5f3//5NTk5OwsTERNy9e1de9vLlS2FtbS1GjBghL+vbt68wNTUVKSkp8rLc3Fzh6uqq8Fldu3ZNABCjRo1SuM+5c+cEAIV4HBwcxJAhQ4QQQmRlZQkzMzMREhIiAMjjmTdvnjA0NBQZGRnv/EwCAwOV+tnz8fF5Zztv2759uwAgjh8/rtJ1b/Py8hL29vbi5cuX8rL09HRhbW2t8P8kISFBABCbNm3K10ZRf3Zyc3NFdna2cHZ2FuPHj5eXHz9+XAAQXbt2Vaj/22+/CQAiOjpaXtatWzeF3y3vun9eu29+Nq6urqJp06YiJydH4doPP/xQ2NnZCalUWmj8Qmjud8jPP/+c79nyDB8+XBgZGb3z+qVLlwoA4uTJkwrlM2fOFABEx44d33k96TbOUSlEjx49FL738PDAq1evkJqaiho1amD//v2QSCQYMGAAcnNz5fVsbW3RuHHjd86sj4qKAgB8+umnCuX9+vXD1KlTixVPcfz+++9FTlIDAHt7e5XaTUxMxIcffghHR0ds2LBB6esKW+Gi7soXOzs7eHp6yr+3trZG9erVUbt2bYVne++99wAAd+/eVet+b8vrXTh27BiCgoJw9OhRdOjQAc7Ozli4cCGA1z1iL168KLInAgCaNGmCWrVqyb83MTGBi4uLQtzHjx+Hn5+fws+Gvr4++vbti9mzZyvUA14PT76pZcuWeO+99/Dnn39i3rx5AAA/Pz95z9CZM2eQmZmJCRMmYOPGjTh69CiGDh0qX7VR0FDgm0JDQzFmzJgin9XCwqLIOpr24sULxMTEYNSoUTAxMVGIpXv37ggPDy9Wu7m5uVi4cCG2bt2K27dvKwzDXbt2LV/9gv7MA69/PjUx/+T27du4fv06vvvuO3l8ebp27Yr9+/fjxo0b8j8XBdH075Di/g4ICAjAnDlzMHz4cGzatAkNGjTAwYMHsXz5cgAlNzxIpYOJSiFsbGwUvs/rQnz58iWA13NGxP/NPyhI3bp1C2378ePHMDAwgLW1tUL5uxKOouIpjoYNGyq1IkeVP+R3795Fu3btYGBggD///DPfMxYk79keP36c7zN48uSJUm28S0HXGxkZ5Ss3MjICAJWWVCvDyckJ9erVw7Fjx9C3b19ER0dj4sSJqF+/PsaOHYsbN27Ih1G8vb2LbO/tnwXg9c/Dmz8Ljx8/hq2tbb56b5c9fvwYwOtk7m329vYKyY+/vz/Cw8Nx69YtHDt2DE2bNkX16tXlQ1v9+/fHmTNnMH369CKfoVatWkpNFNb08mxlPH36FDKZTKnPTxUTJkzAqlWrEBISAh8fH/k8pc8//7zAP8cl8Wf+TXnD05MmTcKkSZMKrFPUUn9N/Q5583fA25T5HVC1alUcOnQIgYGB8iTOxsYGS5YswdChQ1GzZs0iYyTdxUSlmKpWrQqJRIKTJ0/Kf4G8qaCyPDY2NsjNzc33BzAlJaVEYi1MvXr1lOo9mDVrllKT8u7evQtfX18IIRAZGan0ihV3d3cAQHx8PBo2bCgvz83NxfXr19GvXz+l2tFlfn5+2Lt3L6KioiCTyeDr6wsLCwvY29vj6NGjOHbsGN5///13/tyowsbGpsCfp7fL8v6CSE5Ozvf/67///kPVqlUVngF43TOU1yuUVz5jxgycOHECWVlZSvUKDRkyRKmeCR8fn1Lf96NKlSqQSCRKfX55PS5v9yoU9Bfu1q1bMWjQIMyfP1+h/NGjRxpdUq+svP+3U6dOLXQbgKJWymjqd4ibmxuA178DunbtqnAuPj5efv5dWrRogX/++QeJiYl48eIFnJ2dERsbCwD44IMPiryedBcTlWL68MMPsWDBAjx48CDfEE5RfHx8sHDhQvz6668KS+4iIiLUiuntf1UXRZPdtvfu3YOvry+kUikiIyPh5OSkdBxeXl6ws7PD5s2b0bdvX3n5jh07kJGRodJeKrrK398f69atw9KlS9GqVSv5kIafnx92796NmJiYfH+BqaNdu3bYt28fHj58KO+lkkql+PXXXxXqtW/fHsDrv0RbtGghL4+JicG1a9cUekfs7OzQsGFD7Ny5E7GxsfJ4O3TogBEjRmDJkiWwtLRUaKcwujz0Y2ZmhpYtW2LXrl1YtGiRPBl5/vx5vqXPNWrUgImJSb7Jmnv37s3XrkQiyZeIHjhwAA8ePED9+vWLFauqf+bf1KBBAzg7O+PSpUvF/tnT1O+QmjVromXLlti6dSsmTZoEfX19AK9Xrt24cQPjxo1TOqa8FZZCCCxevBj29vb45JNPlL6edA8TlWJq06YNhg8fjqCgIJw/fx4ffPABzMzMkJycjFOnTsHd3b3Qdf+dO3dGmzZtMHHiRKSnp8PT0xPR0dHYsmULgOKPp7q7uyMyMhK///477OzsYGFh8c5/EeX1ZKgrNTUV7dq1Q3JyMjZu3IjU1FSFZaxv7gdy9+5d1KtXD4GBgdi4cSOA13MnFi5ciIEDB2LEiBHo168fbt26hSlTpqBDhw7vXA2lTZs3b0ZQUBA2bdqUb47H29q3bw+JRIIjR44ozBHx9/eX7wGiTE+EsmbMmIF9+/ahffv2+Prrr1GpUiWsWrUKL168UKjXoEEDDB8+HCtWrICenh66dOmCxMREzJw5E46OjvlWuvn5+WHFihUwNTVFmzZtALze66JOnTo4cuQIevToobC8tDC1a9fW2J43mZmZ+OOPPwC8/osNeD0P7NGjRzAzM1NYsjp48GCEh4cjISHhnff/5ptv0LlzZ3To0AETJ06EVCrFt99+CzMzMzx58kReL2+e2o8//oh69eqhcePG+Pvvv7Ft27Z8bX744YfYvHkzXF1d4eHhgdjYWCxatEitvXLc3d2xa9curFmzBp6entDT00Pz5s2Vvn7t2rXo0qULOnXqhMGDB6NmzZp48uQJrl27hri4OGzfvr3I+2vKt99+iw4dOuCTTz7BqFGjkJqaiq+++gpubm4KWxcU9DsEAKZPnw53d3fY2dnh3r17+PHHH3Hu3DkcOHAApqamGouTSh8TFTWsXbsWrVq1wtq1a7F69WrIZDLY29ujTZs28j1ECqKnp4fff/8dEydOxIIFC5CdnY02bdpg69ataNWqVbG7gZctW4bRo0fjs88+Q2ZmZql1m//zzz/yvWMGDBiQ7/yb3b5CCEil0nzLTgcMGAB9fX0sWLAAmzdvhrW1NQYNGiSfyKmL8rb3L2h+x9tsbGzQpEkTXLhwQSEhyfs677ymuLm54dixY5g4cSICAwNRpUoVDBw4EB999BGGDx+uUHfNmjWoV68eNm7ciFWrVsHKygqdO3dGWFhYvnkS/v7+WLFiBdq2basw0dTf3x/r16/XaLKlrNTU1Hz/Ys77eXNyclJ4XUBGRgZMTU2L/DPWoUMH7NmzBzNmzEDfvn1ha2uLUaNG4eXLlwqJJvB6PxkAWLhwITIyMtC+fXvs378/XyK0bNkyGBoaIiwsDBkZGWjWrBl27dql1q7HwcHBuHr1KqZNm4Znz55BCKHUnJE87dq1w99//4158+Zh3LhxePr0KWxsbNCwYUOVe4rV5evriz/++ANff/01unfvjkqVKuHDDz/EokWLFHqiCvsd8vTpU4SEhCAlJQWWlpbw8fHBuXPnNJpMkXZIhCo/1VSitm3bhoCAAJw+fVqpSZX0Wl7Pxu3bt+Hk5KTUv+iVlfdLccuWLRg6dChiYmLk/2L99NNPkZCQgJiYGI3dj0qWra0tBg4ciEWLFhXr+tDQUMyePVulZICI1MMeFS355Zdf8ODBA7i7u0NPTw9nz57FokWL8MEHHzBJKaa8cf7//e9/CpNA1bF371707t07X3nehOGtW7dq5D5U8q5evYrMzEyEhIRoOxQiUgETFS3J2/Fy7ty5ePHiBezs7DB48OACdw6ld+vevbtCr4YmV1D4+voqtJ23KkkikSj9VlfSDY0aNUJ6erq2wyAiFXHoh4iIiHQWt+sjIiIincVEhYiIiHQWExUiIiLSWeV+Mq1MJsN///0HCwsLrbw3hIiIyg4hBJ4/fw57e/sSfZnhq1evkJ2drXY7RkZGCnsalUflPlH577//4OjoqO0wiIioDElKSlJr1+B3efXqFeo4mSMlVVp05SLY2toiISGhXCcr5T5RyXtXyN242rA050gXAR/3KvvvDiINYk8rvSFXmoWo68tL9D1T2dnZSEmVIiHWCZYWxf97Kf25DHU87yI7O5uJSlmWN9xjaa4HSwt9LUdDusBAXzNvKKZygokKFaA0pgpYWuiplahUFOU+USEiItJFUiGDVI2dzKRCprlgdBgTFSIiIi2QQUCG4mcq6lxblrDPiYiIiHSWVhOV0NBQSCQShcPW1lZ+XgiB0NBQ2Nvbw9TUFL6+vrh69aoWIyYiItIMmQb+qwi03qPSqFEjJCcny4/4+Hj5uYULF2LJkiVYuXIlYmJiYGtriw4dOuD58+dajJiIiEh9UiHUPioCrScqBgYGsLW1lR/VqlUD8Lo3ZenSpZg+fTr69OkDNzc3hIeHIzMzE9u2bdNy1ERERFQatJ6o3Lp1C/b29qhTpw4+++wz3LlzBwCQkJCAlJQUdOzYUV7X2NgYPj4+OHPmjLbCJSIi0oi8ybTqHBWBVlf9eHl5YcuWLXBxccHDhw8xd+5ceHt74+rVq0hJSQEA1KhRQ+GaGjVq4O7du4W2mZWVhaysLPn36enpJRM8ERGRGmQQkHLVT5G0mqh06dJF/rW7uztat26NevXqITw8HK1atQKQf9MdIcQ7N+IJCwvD7NmzSyZgIiIiKlVaH/p5k5mZGdzd3XHr1i356p+8npU8qamp+XpZ3jR16lQ8e/ZMfiQlJZVozERERMXBoR/l6FSikpWVhWvXrsHOzg516tSBra0tjh49Kj+fnZ2NqKgoeHt7F9qGsbExLC0tFQ4iIiJdw1U/ytHq0M+kSZPQvXt31KpVC6mpqZg7dy7S09MRGBgIiUSCcePGYf78+XB2doazszPmz5+PSpUqoX///toMm4iISG2y/zvUub4i0Gqicv/+ffTr1w+PHj1CtWrV0KpVK5w9exZOTk4AgClTpuDly5cYNWoUnj59Ci8vLxw5cqRE32pJREREukOriUpERMQ7z0skEoSGhiI0NLR0AiIiIiolUjVX/ahzbVnClxISERFpgVRAzbcnay4WXaZTk2mJiIiI3sQeFSIiIi3gZFrlMFEhIiLSAhkkkKLwDUyVub4i4NAPERER6Sz2qBAREWmBTLw+1Lm+ImCiQkREpAVSNYd+1Lm2LOHQDxEREeks9qgQERFpAXtUlMNEhYiISAtkQgKZUGPVjxrXliVMVIiIiLSAPSrK4RwVIiIi0lnsUSEiItICKfQgVaO/QKrBWHQZExUiIiItEGrOUREVZI4Kh36IiIhIZ7FHhYiISAs4mVY5TFSIiIi0QCr0IBVqzFGpIFvoc+iHiIiIdBZ7VIiIiLRABglkavQXyFAxulSYqBAREWkB56goh0M/REREFUBoaCgkEonCYWtrKz8vhEBoaCjs7e1hamoKX19fXL16VYsRv8ZEhYiISAvyJtOqc6iqUaNGSE5Olh/x8fHycwsXLsSSJUuwcuVKxMTEwNbWFh06dMDz5881+dgq49APERGRFryeo6LGSwmLca2BgYFCL0oeIQSWLl2K6dOno0+fPgCA8PBw1KhRA9u2bcOIESOKHae62KNCRESkBbL/20K/uEfeRNz09HSFIysrq9B73rp1C/b29qhTpw4+++wz3LlzBwCQkJCAlJQUdOzYUV7X2NgYPj4+OHPmTMl+EEVgokJERFSGOTo6wsrKSn6EhYUVWM/LywtbtmzB4cOHsX79eqSkpMDb2xuPHz9GSkoKAKBGjRoK19SoUUN+Tls49ENERKQF6m/49np5clJSEiwtLeXlxsbGBdbv0qWL/Gt3d3e0bt0a9erVQ3h4OFq1agUAkEgUh5OEEPnKSht7VIiIiLRA9n/DN+ocAGBpaalwFJaovM3MzAzu7u64deuWfN7K270nqamp+XpZShsTFSIiogooKysL165dg52dHerUqQNbW1scPXpUfj47OxtRUVHw9vbWYpQc+iEiItIKqZBAKtTY8E3FaydNmoTu3bujVq1aSE1Nxdy5c5Geno7AwEBIJBKMGzcO8+fPh7OzM5ydnTF//nxUqlQJ/fv3L3aMmsBEhYiISAvyVu8U/3rVttC/f/8++vXrh0ePHqFatWpo1aoVzp49CycnJwDAlClT8PLlS4waNQpPnz6Fl5cXjhw5AgsLi2LHqAlMVIiIiCqAiIiId56XSCQIDQ1FaGho6QSkJCYqREREWiATepCpsepHJvhSQiIiIiohpT30U1bpzKqfsLAw+WSePA8fPsTgwYNhb2+PSpUqoXPnzrh165b2giQiIqJSpROJSkxMDNatWwcPDw95mRACvXr1wp07d7B3715cuHABTk5O8Pf3x4sXL7QYLRERkfpk+P8rf4pzyLT9AKVE64lKRkYGAgICsH79elSpUkVefuvWLZw9exZr1qxBixYt0KBBA6xevRoZGRn45ZdftBgxERGR+jS14Vt5p/WnHD16NLp16wZ/f3+F8ryXKpmYmMjL9PX1YWRkhFOnTpVqjERERJqWt4W+OkdFoNXJtBEREYiLi0NMTEy+c66urnBycsLUqVOxdu1amJmZYcmSJUhJSUFycnKhbWZlZSm8OTI9Pb1EYiciIqKSp7V0LCkpCcHBwdi6datCr0keQ0ND7Ny5Ezdv3oS1tTUqVaqEyMhIdOnSBfr6+oW2GxYWpvAWSUdHx5J8DCIiomKRQaL2URFoLVGJjY1FamoqPD09YWBgAAMDA0RFRWH58uUwMDCAVCqFp6cnLl68iLS0NCQnJ+PQoUN4/Pgx6tSpU2i7U6dOxbNnz+RHUlJSKT4VERGRcjj0oxytDf34+fkhPj5eoSwoKAiurq4ICQlR6DWxsrIC8HqC7fnz5/HNN98U2q6xsbHSb44kIiIi3aa1RMXCwgJubm4KZWZmZrCxsZGXb9++HdWqVUOtWrUQHx+P4OBg9OrVCx07dtRGyERERBqj/oZv7FHRuuTkZEyYMAEPHz6EnZ0dBg0ahJkzZ2o7LCIiIrXJhAQyNd6erM61ZYlOJSqRkZEK348dOxZjx47VTjBERESkdTqVqBAREVUUMjWHfirKhm9MVIiIiLRA/bcnV4xEpWI8JREREZVJ7FEhIiLSAikkkKqxaZs615YlTFSIiIi0gEM/ymGiQkREpAVSqNcrItVcKDqtYqRjREREVCaxR4WIiEgLOPSjHCYqREREWqDuiwUryksJK8ZTEhERUZnEHhUiIiItEJBApsZkWsHlyURERFRSOPSjnIrxlERERFQmsUeFiIhIC2RCApko/vCNOteWJUxUiIiItECq5tuT1bm2LKkYT0lERERlEntUiIiItIBDP8phokJERKQFMuhBpsbAhjrXliVMVIiIiLRAKiSQqtEros61ZUnFSMeIiIioTGKPChERkRZwjopymKgQERFpgVDz7cmCO9MSERERaRd7VIiIiLRACgmkarxYUJ1ryxImKkRERFogE+rNM5EJDQajwzj0Q0RERDqLPSpERERaIFNzMq0615YlTFSIiIi0QAYJZGrMM1Hn2rKEiQoREZEWcGda5VSMfiMiIiIqk9ijQkREpAWco6IcJipERERaIIOaW+hXkDkqFSMdIyIiojKJPSpERERaINRc9SPYo1K6wsLCIJFIMG7cOHlZRkYGxowZAwcHB5iamuK9997DmjVrtBckERGRhuS9PVmdoyLQiR6VmJgYrFu3Dh4eHgrl48ePx/Hjx7F161bUrl0bR44cwahRo2Bvb4+ePXtqKVoiIiIqLVrvUcnIyEBAQADWr1+PKlWqKJyLjo5GYGAgfH19Ubt2bQwfPhyNGzfG+fPntRQtERGRZuSt+lHnqAi0/pSjR49Gt27d4O/vn+9c27ZtsW/fPjx48ABCCBw/fhw3b95Ep06dCm0vKysL6enpCgcREZGu4dCPcrQ69BMREYG4uDjExMQUeH758uUYNmwYHBwcYGBgAD09PWzYsAFt27YttM2wsDDMnj27pEImIiKiUqS1RCUpKQnBwcE4cuQITExMCqyzfPlynD17Fvv27YOTkxNOnDiBUaNGwc7OrsAeGACYOnUqJkyYIP8+PT0djo6OJfIMRERExcV3/ShHqUSlT58+Kjf8ww8/oHr16oWej42NRWpqKjw9PeVlUqkUJ06cwMqVK/Hs2TNMmzYNu3fvRrdu3QAAHh4euHjxIr777rtCExVjY2MYGxurHC8REVFpUnf4hkM/b9izZw8+/fRTmJqaKtXotm3bkJGR8c5Exc/PD/Hx8QplQUFBcHV1RUhICKRSKXJycqCnpziNRl9fHzKZTKk4iIiIdBUTFeUoPfSzfPnydyYeb9qxY0eRdSwsLODm5qZQZmZmBhsbG3m5j48PJk+eDFNTUzg5OSEqKgpbtmzBkiVLlA2biIiIyjClEpXjx4/D2tpa6UYPHjyImjVrFjuoPBEREZg6dSoCAgLw5MkTODk5Yd68eRg5cqTabRMREWkTe1SUo1Si4uPjo1Kj71qV8y6RkZEK39va2mLTpk3FaouIiEiXMVFRTrFW/chkMty+fRupqan55ot88MEHGgmMiIiISOVE5ezZs+jfvz/u3r0LIYTCOYlEAqlUqrHgiIiIyisB9ZYYi6KrlAsqJyojR45E8+bNceDAAdjZ2UEiqRhdT0RERJrEoR/lqJyo3Lp1Czt27ED9+vVLIh4iIiIiOZXf9ePl5YXbt2+XRCxEREQVBt/1oxylelQuX74s//rLL7/ExIkTkZKSAnd3dxgaGirU9fDw0GyERERE5RCHfpSjVKLSpEkTSCQShcmzQ4YMkX+dd46TaYmIiMqGsLAwTJs2DcHBwVi6dCkAYPDgwQgPD1eo5+XlhbNnz2ohwteUSlQSEhJKOg4iIqIKRZs9KjExMVi3bl2BoyCdO3dW2MPMyMio2PfRBKUSFScnJ/nXJ06cgLe3NwwMFC/Nzc3FmTNnFOoSERFRwYSQQKiRbBT32oyMDAQEBGD9+vWYO3duvvPGxsawtbUtdlyapvJk2nbt2uHJkyf5yp89e4Z27dppJCgiIqLyTgaJ2gcApKenKxxZWVnvvO/o0aPRrVs3+Pv7F3g+MjIS1atXh4uLC4YNG4bU1FSNP7sqVE5U8uaivO3x48cwMzPTSFBERESkHEdHR1hZWcmPsLCwQutGREQgLi6u0DpdunTBzz//jL/++guLFy9GTEwM2rdvX2TyU5KU3kelT58+AF5PnB08eDCMjY3l56RSKS5fvgxvb2/NR0hERFQOaWqOSlJSEiwtLeXlb/79/KakpCQEBwfjyJEjMDExKbBO37595V+7ubmhefPmcHJywoEDB+R5QGlTOlGxsrIC8LpHxcLCAqampvJzRkZGaNWqFYYNG6b5CImIiMohTc1RsbS0VEhUChMbG4vU1FR4enrKy6RSKU6cOIGVK1ciKysL+vr6CtfY2dnByckJt27dKnac6lI6Udm0aROEEBBCYMWKFbCwsCjJuIiIiEiD/Pz8EB8fr1AWFBQEV1dXhISE5EtSgNfTOpKSkmBnZ1daYeaj0hwVIQS2bduGlJSUkoqHiIioQijtnWktLCzg5uamcJiZmcHGxgZubm7IyMjApEmTEB0djcTERERGRqJ79+6oWrUqevfuXUKfQtFUSlT09PTg7OyMx48fl1Q8REREFULe0I86hybp6+sjPj4ePXv2hIuLCwIDA+Hi4oLo6GitjqKo/FLChQsXYvLkyVizZg3c3NxKIiYiIiIqBZGRkfKvTU1NcfjwYe0FUwiVE5UBAwYgMzMTjRs3hpGRkcKkWgAF7rGiC9pf7gX9SgXPhKaKpfKVG9oOgXRI72uPtB0C6ZCXGbn4s3np3EuouepH0z0qukrlRCXvfQBERERUfALAG6/QK9b1FYHKiUpgYGBJxEFERESUj8qJCvB63fWePXtw7do1SCQSNGzYED169ChwaRMRERHlJ4MEEqix4Zsa15YlKicqt2/fRteuXfHgwQM0aNAAQgjcvHkTjo6OOHDgAOrVq1cScRIREZUr2nopYVmj8rt+xo4di3r16iEpKQlxcXG4cOEC7t27hzp16mDs2LElESMREVG5U9r7qJRVKveoREVF4ezZs7C2tpaX2djYYMGCBWjTpo1GgyMiIqKKTeVExdjYGM+fP89XnpGRASMjI40ERUREVN4Joeaqnwqy7EfloZ8PP/wQw4cPx7lz5+Tv/jl79ixGjhyJHj16lESMRERE5Y6u7Uyrq1ROVJYvX4569eqhdevWMDExgYmJCdq0aYP69etj2bJlJREjERERVVAqD/1UrlwZe/fuxa1bt3D9+nUIIdCwYUPUr1+/JOIjIiIql7jqRznF2kcFAJydneHs7KzJWIiIiCoMmZBAokaywVU/hZBKpdi8eTP+/PNPpKamQiaTKZz/66+/NBYcERERVWwqJyrBwcHYvHkzunXrBjc3N0gkFSOjIyIi0iSu+lGOyolKREQEfvvtN3Tt2rUk4iEiIqoQXicq6sxR0WAwOkzlVT9GRkacOEtERESlQuVEZeLEiVi2bBmEhlO5sLAwSCQSjBs3Tl4mkUgKPBYtWqTRexMREZU27qOiHJWHfk6dOoXjx4/j4MGDaNSoEQwNDRXO79q1S+UgYmJisG7dOnh4eCiUJycnK3x/8OBBDB06FB999JHK9yAiItIl4v8Oda6vCIq1j0rv3r01FkBGRgYCAgKwfv16zJ07V+Gcra2twvd79+5Fu3btULduXY3dn4iISBu4j4pyVE5UNm3apFS906dPo3nz5jA2Nn5nvdGjR6Nbt27w9/fPl6i86eHDhzhw4ADCw8Pf2V5WVhaysrLk36enpysVLxEREekeleeoKKtLly548ODBO+tEREQgLi4OYWFhRbYXHh4OCwsL9OnT5531wsLCYGVlJT8cHR1VipuIiKhUCA0cFUCJJSpFTbZNSkpCcHAwtm7dChMTkyLb+/HHHxEQEFBk3alTp+LZs2fyIykpSaW4iYiISoW6E2k59FOyYmNjkZqaCk9PT3mZVCrFiRMnsHLlSmRlZUFfXx8AcPLkSdy4cQO//vprke0aGxsXOdxEREREZYPWEhU/Pz/Ex8crlAUFBcHV1RUhISHyJAUANm7cCE9PTzRu3Li0wyQiIioR3JlWOVpLVCwsLODm5qZQZmZmBhsbG4Xy9PR0bN++HYsXLy7tEImIiEoMV/0op8QSFU29AygiIgJCCPTr108j7REREVHJuXz5coHlEokEJiYmqFWrlkpTNEosUSnOzrWRkZH5yoYPH47hw4drICIiIiIdou6EWB3tUWnSpMk7OysMDQ3Rt29frF27VqnFNCqv+nn58iUyMzPl39+9exdLly7FkSNHFOo9f/6cG7MREREVIm+OijqHLtq9ezecnZ2xbt06XLx4ERcuXMC6devQoEEDbNu2DRs3bsRff/2FGTNmKNWeyj0qPXv2RJ8+fTBy5EikpaXBy8sLhoaGePToEZYsWYIvvvhC5YciIiKi8mHevHlYtmwZOnXqJC/z8PCAg4MDZs6cib///htmZmaYOHEivvvuuyLbU7lHJS4uDu+//z4AYMeOHahRowbu3r2LLVu2YPny5ao2R0REVDGV0w3f4uPj4eTklK/cyclJvtq3SZMm+d7nVxiVE5XMzExYWFgAAI4cOYI+ffpAT08PrVq1wt27d1VtjoiIqEIqr29PdnV1xYIFC5CdnS0vy8nJwYIFC+Dq6goAePDgAWrUqKFUeyoP/dSvXx979uxB7969cfjwYYwfPx4AkJqaCktLS1WbIyIiqrh0tFdEHatWrUKPHj3g4OAADw8PSCQSXL58GVKpFPv37wcA3LlzB6NGjVKqPZUTla+//hr9+/fH+PHj0b59e7Ru3RrA696Vpk2bqtocERERlSPe3t5ITEzE1q1bcfPmTQgh8PHHH6N///7yEZmBAwcq3Z7KicrHH3+Mtm3bIjk5WWGnWD8/P/Tu3VvV5oiIiCqk8rzhm7m5OUaOHKmRtoq1j4qtrS0yMjJw9OhRfPDBBzA1NUWLFi00tskbERFRuafuhFgdHja6efMmIiMjkZqaCplMpnDu66+/VqktlROVx48f49NPP8Xx48chkUhw69Yt1K1bF59//jkqV67Mre6JiIgqsPXr1+OLL75A1apVYWtrq9CJIZFIVE5UVF71M378eBgaGuLevXuoVKmSvLxv3744dOiQqs0RERFVUBINHLpn7ty5mDdvHlJSUuQbvuUdcXFxKrenco/KkSNHcPjwYTg4OCiUOzs7c3kyERGRssrp0M/Tp0/xySefaKw9lXtUXrx4odCTkufRo0cqvWSIiIiIyp9PPvkk32t11KFyj8oHH3yALVu24JtvvgHwerxJJpNh0aJFaNeuncYCIyIiKtfKaY9K/fr1MXPmTJw9exbu7u4wNDRUOD927FiV2lM5UVm0aBF8fX1x/vx5ZGdnY8qUKbh69SqePHmC06dPq9ocERFRxVRO3568bt06mJubIyoqClFRUQrnJBJJyScqDRs2xOXLl7F69Wro6+vjxYsX6NOnD0aPHg07OztVmyMiIqJyJCEhQaPtFXsflTlz5mg0ECIioopEiNeHOtdXBMVKVE6ePIm1a9fizp072L59O2rWrImffvoJderUQdu2bTUdIxERUflTjuaoTJgwAd988w3MzMwwYcKEd9ZdsmSJSm2rnKjs3LkTAwcOREBAAOLi4pCVlQUAeP78OebPn48//vhD1SaJiIgqnnI0R+XChQvIyckBAMTFxRW6U31xdrBXOVGZO3cufvjhBwwaNAgRERHycm9vbw4HERERVUDHjx+Xfx0ZGanRtlXeR+XGjRv44IMP8pVbWloiLS1NEzERERGVexKh/qFrcnNzYWBggCtXrmisTZUTFTs7O9y+fTtf+alTp1C3bl2NBEVERFTuCQ0cOsbAwABOTk6QSqUaa1PlRGXEiBEIDg7GuXPnIJFI8N9//+Hnn3/GpEmTMGrUKI0FRkRERGXPjBkzMHXqVDx58kQj7ak8R2XKlCl49uwZ2rVrh1evXuGDDz6AsbExJk2ahDFjxmgkKCIionKvHE2mfdPy5ctx+/Zt2Nvbw8nJCWZmZgrnVX0xoUqJilQqxalTpzBx4kRMnz4d//zzD2QyGRo2bAhzc3OVbkxERFShlaPlyW/q1auXRttTKVHR19dHp06dcO3aNVhbW6N58+YaDYaIiIjKtlmzZmm0PZWHftzd3XHnzh3UqVNHo4EQERFVKOW0RyXP+fPnce3aNUgkErz33nvw9PQsVjsqJyrz5s3DpEmT8M0338DT0zPf2JOlpWWxAiEiIqpQymmicv/+ffTr1w+nT59G5cqVAQBpaWnw9vbGL7/8AkdHR5XaU3nVT+fOnXHp0iX06NEDDg4OqFKlCqpUqYLKlSujSpUqqjZHRERE5ciQIUOQk5ODa9eu4cmTJ3jy5AmuXbsGIQSGDh2qcnsq96i8ufscERERFVM5XfVz8uRJnDlzBg0aNJCXNWjQACtWrECbNm1Ubk/lRMXHx0flmxAREZEidXeX1cWdaQGgVq1a8vf+vCk3Nxc1a9ZUuT2VE5XLly8XWC6RSGBiYoJatWrB2NhY5UCIiIgqlHI6R2XhwoX48ssvsWrVKnh6ekIikeD8+fMIDg7Gd999p3J7KicqTZo0eefbDw0NDdG3b1+sXbsWJiYmKgdEREREZdfgwYORmZkJLy8vGBi8TjPy3gE0ZMgQDBkyRF5Xmd1rVU5Udu/ejZCQEEyePBktW7aEEAIxMTFYvHgxZs2ahdzcXHz11VeYMWNGsTInIiIiKruWLl2q0faKtTx52bJl6NSpk7zMw8MDDg4OmDlzJv7++2+YmZlh4sSJKiUqYWFhmDZtGoKDgxUe8tq1awgJCUFUVBRkMhkaNWqE3377DbVq1VI1dCIiIp0hgZpzVDQWiWYFBgZqtD2VlyfHx8fDyckpX7mTkxPi4+MBvB4eSk5OVrrNmJgYrFu3Dh4eHgrl//77L9q2bQtXV1dERkbi0qVLmDlzJoeUiIiIdFRcXJw8HwCAvXv3olevXpg2bRqys7NVbk/lRMXV1RULFixQuFlOTg4WLFgAV1dXAMCDBw9Qo0YNpdrLyMhAQEAA1q9fn28flunTp6Nr165YuHAhmjZtirp166Jbt26oXr26qmETERHplrzlyeocOmjEiBG4efMmAODOnTvo27cvKlWqhO3bt2PKlCkqt6dyorJq1Srs378fDg4O8Pf3R4cOHeDg4ID9+/djzZo18sBGjRqlVHujR49Gt27d4O/vr1Auk8lw4MABuLi4oFOnTqhevTq8vLywZ8+ed7aXlZWF9PR0hYOIiEjnCA0cOujmzZto0qQJAGD79u3w8fHBtm3bsHnzZuzcuVPl9lROVLy9vZGYmIg5c+bAw8MDbm5umDNnDhISEtCqVSsAwMCBAzF58uQi24qIiEBcXBzCwsLynUtNTUVGRgYWLFiAzp0748iRI+jduzf69OmDqKioQtsMCwuDlZWV/FB1q14iIqKKICwsDBKJBOPGjZOXCSEQGhoKe3t7mJqawtfXF1evXlWpXSEEZDIZAODYsWPo2rUrAMDR0RGPHj1SOU6VJ9MCgLm5OUaOHFmcS+WSkpIQHByMI0eOFDjnJO8he/bsifHjxwN4PfflzJkz+OGHHwrdeG7q1KmYMGGC/Pv09HQmK0REpHu0uI9KYXNDFy5ciCVLlmDz5s1wcXHB3Llz0aFDB9y4cQMWFhZKtd28eXPMnTsX/v7+iIqKko+2JCQkKD0t5E0q96gAwE8//YS2bdvC3t4ed+/eBQB8//332Lt3r9JtxMbGIjU1FZ6enjAwMICBgQGioqKwfPlyGBgYwMbGBgYGBmjYsKHCde+99x7u3btXaLvGxsawtLRUOIiIiHRN3s606hzFUdjcUCEEli5diunTp6NPnz5wc3NDeHg4MjMzsW3bNqXbX7p0KeLi4jBmzBhMnz4d9evXBwDs2LED3t7eKsercqKyZs0aTJgwAV26dMHTp08hlUoBAFWqVFFp7bSfnx/i4+Nx8eJF+dG8eXMEBATg4sWLMDY2RosWLXDjxg2F627evFngqiMiIiIqWmFzQxMSEpCSkoKOHTvKy4yNjeHj44MzZ84o3b6Hhwfi4+Px7NkzzJo1S16+aNEihIeHy7//5Zdf8OLFiyLbU3noZ8WKFVi/fj169eqFBQsWyMubN2+OSZMmKd2OhYUF3NzcFMrMzMxgY2MjL588eTL69u2LDz74AO3atcOhQ4fw+++/IzIyUtWwiYiIdIuGhn7eXjRibGxc6Kts8uaGxsTE5DuXkpICAPmGZ2rUqCEfPVHH29M8RowYAS8vL9StW/ed16nco5KQkICmTZvmKzc2NlYqM1JF79698cMPP2DhwoVwd3fHhg0bsHPnTrRt21aj9yEiIip1Glr14+joqLCIpKAFKsD/nxu6devWd+5H9vZrcoQQ73x1TnEJoVyWpnKPSp06dXDx4sV8wy8HDx7MN59EVQX1lLz9XgAiIqLyQFNvT05KSlKYj1lYb8qbc0PzSKVSnDhxAitXrpRPtUhJSYGdnZ28TmpqarEmwWqKyonK5MmTMXr0aLx69QpCCPz999/45ZdfEBYWhg0bNpREjERERFQIZReO5M0NfVNQUBBcXV0REhKCunXrwtbWFkePHpWPnGRnZyMqKgrffvtticSuDJUTlaCgIOTm5mLKlCnIzMxE//79UbNmTSxbtgyfffZZScRIRERU/qi7u6yK1yozN3TcuHGYP38+nJ2d4ezsjPnz56NSpUro379/8eNUU7H2URk2bBiGDRuGR48eQSaTcUt7IiIiVWlxH5XCTJkyBS9fvsSoUaPw9OlTeHl54ciRI0rvoVISipWo5Klataqm4iAiIqJS9vbcUIlEgtDQUISGhpb4vZ2cnGBoaFhkPaUSlaZNmyo94zcuLk6pekRERBWZpibT6qrz58/j2rVrkEgkcHV1RfPmzRXOX7lyRal2lEpUevXqJf/61atXWL16NRo2bIjWrVsDAM6ePYurV68q/SJCIiKiCk8Hh3404f79++jXrx9Onz6NypUrAwDS0tLg7e2NX375ReXX2iiVqLy5s9znn3+OsWPH4ptvvslXJykpSaWbExERUfkyZMgQ5OTk4Nq1a2jQoAEA4MaNGxgyZAiGDh2KI0eOqNSeynNUtm/fjvPnz+crHzBgAJo3b44ff/xR1SaJiIgqHjWHfnS1R+XkyZM4c+aMPEkBgAYNGmDFihVo06aNyu2pvDOtqakpTp06la/81KlT79zpjoiIiN6goZ1pdU2tWrWQk5OTrzw3Nxc1a9ZUuT2Ve1TGjRuHL774ArGxsWjVqhWA13NUfvzxR3z99dcqB0BERETlx8KFC/Hll19i1apV8PT0hEQiwfnz5xEcHIzvvvtO5fZUTlS++uor1K1bF8uWLZO/9vm9997D5s2b8emnn6ocABERUYVUTifTDh48GJmZmfDy8oKBwes0Izc3FwYGBvlei/PkyZMi2yvWPiqffvopkxIiIiI1lNflyUuXLtVoe2pt+EZERET0psDAQI22p1SiYm1tjZs3byq9E22tWrVw8uTJfG9YJiIiooohNTUVqampkMlkCuUeHh4qtaNUopKWloaDBw/CyspKqUYfP34MqVSqUiBEREQVSjmdoxIbG4vAwEBcu3YNQigGKZFIVM4PlB760XRXDhERUUVWXueoBAUFwcXFBRs3bkSNGjWUfgVPYZRKVN7utiEiIiIqSEJCAnbt2oX69etrpD2VN3wjIiIiDSlnm70BgJ+fHy5duqSx9rjqh4iISBvK6RyVDRs2IDAwEFeuXIGbmxsMDQ0Vzvfo0UOl9pioEBERkcacOXMGp06dwsGDB/OdK85kWg79EBERaUHeZFp1Dl00duxYDBw4EMnJyZDJZApHcVYEs0eFiIhIG8rp0M/jx48xfvx41KhRQyPtFatH5d9//8WMGTPQr18/pKamAgAOHTqEq1evaiQoIiIiKpv69OmD48ePa6w9lXtUoqKi0KVLF7Rp0wYnTpzAvHnzUL16dVy+fBkbNmzAjh07NBYcERFReVVe91FxcXHB1KlTcerUKbi7u+ebTDt27FiV2ivW25Pnzp2LCRMmwMLCQl7erl07LFu2TNXmiIiIKqZyOvSzYcMGmJubIyoqClFRUQrnJBJJyScq8fHx2LZtW77yatWq4fHjx6o2R0REVDGV00QlISFB/nXeFvrq7E6r8hyVypUrIzk5OV/5hQsXULNmzWIHQkREROXDxo0b4ebmBhMTE5iYmMDNzQ0bNmwoVlsq96j0798fISEh2L59OyQSCWQyGU6fPo1JkyZh0KBBxQqCiIiooimvc1RmzpyJ77//Hl9++SVat24NAIiOjsb48eORmJiIuXPnqtSeyonKvHnzMHjwYNSsWRNCCDRs2BBSqRT9+/fHjBkzVG2OiIioYiqnQz9r1qzB+vXr0a9fP3lZjx494OHhgS+//LLkExVDQ0P8/PPPmDNnDi5cuACZTIamTZvC2dlZ1aaIiIionJFKpWjevHm+ck9PT+Tm5qrcXrE3fKtXrx7q1atX3MuJiIgqtnLaozJgwACsWbMGS5YsUShft24dAgICVG5PqURlwoQJSjf4dmBERESUX3mdowK8nkx75MgRtGrVCgBw9uxZJCUlYdCgQQo5hTI5g1KJyoULFxS+j42NhVQqRYMGDQAAN2/ehL6+Pjw9PZV+CCIiIip/rly5gmbNmgF4vZM98HoLk2rVquHKlSvyesouWVYqUXlzK9wlS5bAwsIC4eHhqFKlCgDg6dOnCAoKwvvvv6/cUxQgLCwM06ZNQ3BwMJYuXQoAGDx4MMLDwxXqeXl54ezZs8W+DxERkU4op0M/mtw+HyjGHJXFixfjyJEj8iQFAKpUqYK5c+eiY8eOmDhxospBxMTEYN26dfDw8Mh3rnPnzti0aZP8eyMjI5XbJyIi0jXleehHk1Te8C09PR0PHz7MV56amornz5+rHEBGRgYCAgKwfv16heQnj7GxMWxtbeWHtbW1yvcgIiKisknlRKV3794ICgrCjh07cP/+fdy/fx87duzA0KFD0adPH5UDGD16NLp16wZ/f/8Cz0dGRqJ69epwcXHBsGHD5G9rLkxWVhbS09MVDiIiIp0jNHBUACoP/fzwww+YNGkSBgwYgJycnNeNGBhg6NChWLRokUptRUREIC4uDjExMQWe79KlCz755BM4OTkhISEBM2fORPv27REbGwtjY+MCrwkLC8Ps2bNVeygiIqLSVk7nqGiayolKpUqVsHr1aixatAj//vsvhBCoX78+zMzMVGonKSkJwcHBOHLkCExMTAqs07dvX/nXbm5uaN68OZycnHDgwIFCe2+mTp2qsPQpPT0djo6OKsVGRERU0iT/d6hzfUVQ7A3fzMzMCpz8qqzY2FikpqYqLGmWSqU4ceIEVq5ciaysLOjr6ytcY2dnBycnJ9y6davQdo2NjQvtbSEiIqKyReVEpV27du9c+/zXX38p1Y6fnx/i4+MVyoKCguDq6oqQkJB8SQoAPH78GElJSbCzs1MtaCIiIl3DoR+lqJyoNGnSROH7nJwcXLx4EVeuXEFgYKDS7VhYWMDNzU2hzMzMDDY2NnBzc0NGRgZCQ0Px0Ucfwc7ODomJiZg2bRqqVq2K3r17qxo2ERGRTuHyZOWonKh8//33BZaHhoYiIyND7YDy6OvrIz4+Hlu2bEFaWhrs7OzQrl07/Prrr7CwsNDYfYiIiEh3FXuOytsGDBiAli1b4rvvvit2G5GRkfKvTU1NcfjwYQ1ERkREpIM49KMUjSUq0dHRha7eISIiogJUkGRDHSonKm8vCxZCIDk5GefPn8fMmTM1FhgRERGRyomKpaWlwqofPT09NGjQAHPmzEHHjh01GhwREVF5xcm0ylE5Udm8eXMJhEFERFTBcI6KUlR+10/dunXx+PHjfOVpaWmoW7euRoIiIiIiAorRo5KYmAipVJqvPCsrCw8ePNBIUEREROUdh36Uo3Sism/fPvnXhw8fhpWVlfx7qVSKP//8E7Vr19ZocEREROUWh36UonSi0qtXLwCARCLJtwOtoaEhateujcWLF2s0OCIiovKKPSrKUTpRkclkAIA6deogJiYGVatWLbGgiIiIiIBizFFJSEgoiThKnM0Xz2CgZ6TtMEgHGETZajsE0iFrbjTUdgikQ6SZWQBOlc7NOPSjFKUSleXLl2P48OEwMTHB8uXL31l37NixGgmMiIioXGOiohSlEpXvv/8eAQEBMDExKfSlhMDr+StMVIiIiEhTlEpU3hzuKatDP0RERLqEk2mVo/KGb3PmzEFmZma+8pcvX2LOnDkaCYqIiKjcExo4KgCVE5XZs2cjIyMjX3lmZiZmz56tkaCIiIiIgGKs+hFCKLyUMM+lS5dgbW2tkaCIiIjKO4kQkIjid4uoc21ZonSiUqVKFUgkEkgkEri4uCgkK1KpFBkZGRg5cmSJBElERFTucNWPUpROVJYuXQohBIYMGYLZs2crbKFvZGSE2rVro3Xr1iUSJBEREVVMSicqedvm16lTB97e3jA0NCyxoIiIiMo7rvpRjspzVHx8fORfv3z5Ejk5OQrnLS0t1Y+KiIiovOPQj1JUXvWTmZmJMWPGoHr16jA3N0eVKlUUDiIiIipaXo+KOkdFoHKiMnnyZPz1119YvXo1jI2NsWHDBsyePRv29vbYsmVLScRIREREFZTKQz+///47tmzZAl9fXwwZMgTvv/8+6tevDycnJ/z8888ICAgoiTiJiIjKFw79KEXlHpUnT56gTp06AF7PR3ny5AkAoG3btjhx4oRmoyMiIiqnOPSjHJUTlbp16yIxMREA0LBhQ/z2228AXve0VK5cWZOxERERUQWncqISFBSES5cuAQCmTp0qn6syfvx4TJ48WeMBEhERlUul/K6fNWvWwMPDA5aWlrC0tETr1q1x8OBB+fnBgwfLN3bNO1q1aqXmQ6pP5Tkq48ePl3/drl07XL9+HefPn0e9evXQuHFjjQZHRERUnpXm8I2DgwMWLFiA+vXrAwDCw8PRs2dPXLhwAY0aNQIAdO7cGZs2bZJfY2RkVHoBFkLlROVttWrVQq1atZCUlIQhQ4bgxx9/1ERcREREpEHdu3dX+H7evHlYs2YNzp49K09UjI2NYWtrq43wCqXy0E9hnjx5gvDwcE01R0REVL4Jof4BID09XeHIysoq8tZSqRQRERF48eKFwutvIiMjUb16dbi4uGDYsGFITU0tscdXlsYSFSIiIlKeplb9ODo6wsrKSn6EhYUVes/4+HiYm5vD2NgYI0eOxO7du9GwYUMAQJcuXfDzzz/jr7/+wuLFixETE4P27dsrlfiUJLWHfoiIiEh7kpKSFF5fY2xsXGjdBg0a4OLFi0hLS8POnTsRGBiIqKgoNGzYEH379pXXc3NzQ/PmzeHk5IQDBw6gT58+JfoM78JEhYiISBs0tOFb3ioeZRgZGckn0zZv3hwxMTFYtmwZ1q5dm6+unZ0dnJyccOvWLTWCVJ/SiUpR2VRaWpq6sRAREVUYEtnrQ53r1SWEKHRo5/Hjx0hKSoKdnZ36N1KD0nNU3hz/KuhwcnLCoEGDih1IWFgYJBIJxo0bV+D5ESNGQCKRYOnSpcW+BxERkc4o5X1Upk2bhpMnTyIxMRHx8fGYPn06IiMjERAQgIyMDEyaNAnR0dFITExEZGQkunfvjqpVq6J3796aed5iUrpH5c111ZoWExODdevWwcPDo8Dze/bswblz52Bvb19iMRAREZVnDx8+xMCBA5GcnAwrKyt4eHjg0KFD6NChA16+fIn4+Hhs2bIFaWlpsLOzQ7t27fDrr7/CwsJCq3FrfY5KRkYGAgICsH79esydOzff+QcPHmDMmDE4fPgwunXrpoUIiYiINE/d9/Woeu3GjRsLPWdqaorDhw8XP5gSpPXlyaNHj0a3bt3g7++f75xMJsPAgQMxefJk+WY0RcnKysq3ppyIiEjnaGgflfJOqz0qERERiIuLQ0xMTIHnv/32WxgYGGDs2LFKtxkWFobZs2drKkQiIiLSIq31qCQlJSE4OBhbt26FiYlJvvOxsbFYtmwZNm/eDIlEonS7U6dOxbNnz+RHUlKSJsMmIiLSCE1t+FbeaS1RiY2NRWpqKjw9PWFgYAADAwNERUVh+fLlMDAwQGRkJFJTU1GrVi35+bt372LixImoXbt2oe0aGxvL15SrsraciIioVJXyqp+ySmtDP35+foiPj1coCwoKgqurK0JCQmBnZ4dOnTopnO/UqRMGDhyIoKCg0gyViIiItERriYqFhQXc3NwUyszMzGBjYyMvt7GxUThvaGgIW1tbNGjQoNTiJCIiKgmlveqnrNL68mQiIqIKSd2VO1z1U/oiIyPfeT4xMbFU4iAiIiLdoFOJChERUUXBoR/lMFEhIiLSBg29Pbm8Y6JCRESkBexRUY7Wt9AnIiIiKgx7VIiIiLRBJl4f6lxfATBRISIi0gbOUVEKh36IiIhIZ7FHhYiISAskUHMyrcYi0W1MVIiIiLSBO9MqhUM/REREpLPYo0JERKQF3EdFOUxUiIiItIGrfpTCoR8iIiLSWexRISIi0gKJEJCoMSFWnWvLEiYqRERE2iD7v0Od6ysAJipERERawB4V5XCOChEREeks9qgQERFpA1f9KIWJChERkTZwZ1qlcOiHiIiIdBZ7VIiIiLSAO9Mqh4kKERGRNnDoRykc+iEiIiKdxR4VIiIiLZDIXh/qXF8RMFEhIiLSBg79KIVDP0RERKSz2KNCRESkDdzwTSlMVIiIiLSA7/pRDhMVIiIibeAcFaVwjgoRERHpLPaoEBERaYMAoM4S44rRocJEhYiISBs4R0U5HPohIiIinaUziUpYWBgkEgnGjRsnLwsNDYWrqyvMzMxQpUoV+Pv749y5c9oLkoiISFME/v+E2mId2n6A0qETiUpMTAzWrVsHDw8PhXIXFxesXLkS8fHxOHXqFGrXro2OHTvif//7n5YiJSIi0hC1khQ1VwyVIVpPVDIyMhAQEID169ejSpUqCuf69+8Pf39/1K1bF40aNcKSJUuQnp6Oy5cvaylaIiIiKk1aT1RGjx6Nbt26wd/f/531srOzsW7dOlhZWaFx48aF1svKykJ6errCQUREpHNkGjgqAK2u+omIiEBcXBxiYmIKrbN//3589tlnyMzMhJ2dHY4ePYqqVasWWj8sLAyzZ88uiXCJiIg0hqt+lKO1HpWkpCQEBwdj69atMDExKbReu3btcPHiRZw5cwadO3fGp59+itTU1ELrT506Fc+ePZMfSUlJJRE+ERGRejhHRSlaS1RiY2ORmpoKT09PGBgYwMDAAFFRUVi+fDkMDAwglUoBAGZmZqhfvz5atWqFjRs3wsDAABs3biy0XWNjY1haWiocREREVDZpbejHz88P8fHxCmVBQUFwdXVFSEgI9PX1C7xOCIGsrKzSCJGIiKjk8F0/StFaomJhYQE3NzeFMjMzM9jY2MDNzQ0vXrzAvHnz0KNHD9jZ2eHx48dYvXo17t+/j08++URLURMREWkIExWl6OwW+vr6+rh+/TrCw8Px6NEj2NjYoEWLFjh58iQaNWqk7fCIiIioFOhUohIZGSn/2sTEBLt27dJeMERERCVJBkCi5vUVgE4lKkRERBUFlycrR+sbvhEREREVhj0qRERE2sDJtEphokJERKQNMgFI1Eg2ZBUjUeHQDxEREeks9qgQERFpA4d+lMJEhYiISCvUfV8PExUiIiIqKexRUQrnqBAREZHOYo8KERGRNsgE1Bq+qSCrfpioEBERaYOQvT7Uub4C4NAPERER6SwmKkRERNqQN5lWnUMFa9asgYeHBywtLWFpaYnWrVvj4MGDb4QjEBoaCnt7e5iamsLX1xdXr17V9FOrjIkKERGRNsiE+ocKHBwcsGDBApw/fx7nz59H+/bt0bNnT3kysnDhQixZsgQrV65ETEwMbG1t0aFDBzx//rwknl5pTFSIiIgqgO7du6Nr165wcXGBi4sL5s2bB3Nzc5w9exZCCCxduhTTp09Hnz594ObmhvDwcGRmZmLbtm1ajZuJChERkTZoaOgnPT1d4cjKyiry1lKpFBEREXjx4gVat26NhIQEpKSkoGPHjvI6xsbG8PHxwZkzZ0rsI1AGExUiIiJtEFAzUXndjKOjI6ysrORHWFhYobeMj4+Hubk5jI2NMXLkSOzevRsNGzZESkoKAKBGjRoK9WvUqCE/py1cnkxERFSGJSUlwdLSUv69sbFxoXUbNGiAixcvIi0tDTt37kRgYCCioqLk5yUSiUJ9IUS+stLGRIWIiEgbNLSFft4qHmUYGRmhfv36AIDmzZsjJiYGy5YtQ0hICAAgJSUFdnZ28vqpqan5ellKG4d+iIiItEEmU/9QkxACWVlZqFOnDmxtbXH06FH5uezsbERFRcHb21vt+6iDPSpERETaUMovJZw2bRq6dOkCR0dHPH/+HBEREYiMjMShQ4cgkUgwbtw4zJ8/H87OznB2dsb8+fNRqVIl9O/fv/gxagATFSIiogrg4cOHGDhwIJKTk2FlZQUPDw8cOnQIHTp0AABMmTIFL1++xKhRo/D06VN4eXnhyJEjsLCw0GrcTFSIiIi0oZR7VDZu3PjO8xKJBKGhoQgNDS1+TCWAiQoREZE28O3JSuFkWiIiItJZ7FEhIiLSAiFkEKL4K3fUubYsYaJCRESkDUL1Fwvmu74C4NAPERER6Sz2qBAREWmDUHMybQXpUWGiQkREpA0yGSBRY55JBZmjwqEfIiIi0lnsUSEiItIGDv0ohYkKERGRFgiZDEKNoZ+KsjxZZ4Z+wsLC5C9FAoCcnByEhITA3d0dZmZmsLe3x6BBg/Dff/9pN1AiIiJNyNtCX52jAtCJRCUmJgbr1q2Dh4eHvCwzMxNxcXGYOXMm4uLisGvXLty8eRM9evTQYqRERERUmrQ+9JORkYGAgACsX78ec+fOlZdbWVnh6NGjCnVXrFiBli1b4t69e6hVq1Zph0pERKQ5MgFIOEelKFrvURk9ejS6desGf3//Ius+e/YMEokElStXLrROVlYW0tPTFQ4iIiKdI8TrJcbFPipGoqLVHpWIiAjExcUhJiamyLqvXr3CV199hf79+8PS0rLQemFhYZg9e7YmwyQiIiIt0VqPSlJSEoKDg7F161aYmJi8s25OTg4+++wzyGQyrF69+p11p06dimfPnsmPpKQkTYZNRESkEUIm1D4qAq31qMTGxiI1NRWenp7yMqlUihMnTmDlypXIysqCvr4+cnJy8OmnnyIhIQF//fXXO3tTAMDY2BjGxsYlHT4REZF6hAwAd6YtitYSFT8/P8THxyuUBQUFwdXVFSEhIQpJyq1bt3D8+HHY2NhoKVoiIiLSBq0lKhYWFnBzc1MoMzMzg42NDdzc3JCbm4uPP/4YcXFx2L9/P6RSKVJSUgAA1tbWMDIy0kbYREREGiFkAkKNVT+Ck2m16/79+9i3bx8AoEmTJgrnjh8/Dl9f39IPioiISFM49KMUnUpUIiMj5V/Xrl1bI9liXhu5smy126LyQbzQ13YIpEOkmVnaDoF0SN7PQ2n0VuQiR61X/eQiR3PB6DCdSlRKwvPnzwEAkY9/0nIkpDO6aDsAItJ1z58/h5WVVYm0bWRkBFtbW5xK+UPttmxtbcv9VAiJKOeDXDKZDP/99x8sLCwgkUgUzqWnp8PR0RFJSUlFribSNYxdOxi7djB27aiIsQsh8Pz5c9jb20NPr+R28Hj16hWys9Xv6TcyMipyi4+yrtz3qOjp6cHBweGddSwtLcvcH8I8jF07GLt2MHbtqGixl1RPyptMTEzKfYKhKVrfQp+IiIioMExUiIiISGdV6ETF2NgYs2bNKpM72TJ27WDs2sHYtYOxky4o95NpiYiIqOyq0D0qREREpNuYqBAREZHOYqJCREREOouJChEREemscpeohIWFoUWLFrCwsED16tXRq1cv3LhxQ6FORkYGxowZAwcHB5iamuK9997DmjVrFOpkZWXhyy+/RNWqVWFmZoYePXrg/v37ZSL2devWwdfXF5aWlpBIJEhLSyvRuDUV+5MnT/Dll1+iQYMGqFSpEmrVqoWxY8fi2bNnOh87AIwYMQL16tWDqakpqlWrhp49e+L69etlIvY8Qgh06dIFEokEe/bsKROx+/r6QiKRKByfffZZmYgdAKKjo9G+fXuYmZmhcuXK8PX1xcuXL3U69sTExHyfed6xfft2nY4dAFJSUjBw4EDY2trCzMwMzZo1w44dO0osblKTKGc6deokNm3aJK5cuSIuXrwounXrJmrVqiUyMjLkdT7//HNRr149cfz4cZGQkCDWrl0r9PX1xZ49e+R1Ro4cKWrWrCmOHj0q4uLiRLt27UTjxo1Fbm6uzsf+/fffi7CwMBEWFiYAiKdPn5ZYzJqMPT4+XvTp00fs27dP3L59W/z555/C2dlZfPTRRzofuxBCrF27VkRFRYmEhAQRGxsrunfvLhwdHcvEz0yeJUuWiC5duggAYvfu3SUWtyZj9/HxEcOGDRPJycnyIy0trUzEfubMGWFpaSnCwsLElStXxM2bN8X27dvFq1evdDr23Nxchc87OTlZzJ49W5iZmYnnz5/rdOxCCOHv7y9atGghzp07J/7991/xzTffCD09PREXF1disVPxlbtE5W2pqakCgIiKipKXNWrUSMyZM0ehXrNmzcSMGTOEEEKkpaUJQ0NDERERIT//4MEDoaenJw4dOlQ6gYvixf6m48ePl1qi8jZ1Y8/z22+/CSMjI5GTk1Nisb5NU7FfunRJABC3b98usVjfpk7sFy9eFA4ODiI5OblUEpW3FTd2Hx8fERwcXFphFqi4sXt5eb3zZ6g0aOrnvUmTJmLIkCElFmdBihu7mZmZ2LJli0Ida2trsWHDhpINmIql3A39vC1v2MDa2lpe1rZtW+zbtw8PHjyAEALHjx/HzZs30alTJwBAbGwscnJy0LFjR/k19vb2cHNzw5kzZ3Q6dl2hqdifPXsGS0tLGBiU3mupNBH7ixcvsGnTJtSpUweOjo6lEjdQ/NgzMzPRr18/rFy5Era2tqUW75vU+dx//vlnVK1aFY0aNcKkSZPkb03X5dhTU1Nx7tw5VK9eHd7e3qhRowZ8fHxw6tQpnY/9bbGxsbh48SKGDh1aKjHnKW7sbdu2xa+//oonT55AJpMhIiICWVlZ8PX1LdX4SUlaS5FKgUwmE927dxdt27ZVKM/KyhKDBg0SAISBgYEwMjJSyK5//vlnYWRklK+9Dh06iOHDh5d43EIUP/Y3aatHRROxCyHEo0ePRK1atcT06dNLOmQ5dWNftWqVMDMzEwCEq6trqfamqBP78OHDxdChQ+Xfo5R7VNSJfd26deLo0aMiPj5e/PLLL6J27drC399f52OPjo4WAIS1tbX48ccfRVxcnBg3bpwwMjISN2/e1OnY3/bFF1+I9957r6TDVaBO7GlpaaJTp07yOpaWluLIkSOlGT6poFwnKqNGjRJOTk4iKSlJoXzRokXCxcVF7Nu3T1y6dEmsWLFCmJubi6NHjwohCk9U/P39xYgRI3Q69jdpK1HRROzPnj0TXl5eonPnziI7O7u0Qlc79rS0NHHz5k0RFRUlunfvLpo1ayZevnyp07Hv3btX1K9fX2FuQWknKpr4mclz/vx5AUDExsaWdNhCiOLHfvr0aQFATJ06VeE6d3d38dVXX+l07G/KzMwUVlZW4rvvviuVmPOoE/uYMWNEy5YtxbFjx8TFixdFaGiosLKyEpcvXy7VZyDllNtEZcyYMcLBwUHcuXNHoTwzM1MYGhqK/fv3K5QPHTpUdOrUSQghxJ9//ikAiCdPnijU8fDwEF9//XXJBi7Ui/1N2khUNBF7enq6aN26tfDz8yu1v+SF0NznnicrK0tUqlRJbNu2rUTifZM6sQcHBwuJRCL09fXlBwChp6cnfHx8dDr2gshksnxzzEqKOrHfuXNHABA//fSTQp1PP/1U9O/fv2QDF5r73Lds2SIMDQ1Fampqicb7JnViv337tgAgrly5olDHz8+v1P4hSqopd3NUhBAYM2YMdu3ahb/++gt16tRROJ+Tk4OcnBzo6Sk+ur6+PmQyGQDA09MThoaGOHr0qPx8cnIyrly5Am9vb52OXVs0FXt6ejo6duwIIyMj7Nu3DyYmJmUm9sLazsrK0njMb7avbuxfffUVLl++jIsXL8oPAPj++++xadMmnY69IFevXkVOTg7s7OxKJG5AM7HXrl0b9vb2+ZbX3rx5E05OTjod+5s2btyIHj16oFq1aiUWcx5NxJ6ZmQkAOvl7lAqhrQyppHzxxRfCyspKREZGKiydy8zMlNfx8fERjRo1EsePHxd37twRmzZtEiYmJmL16tXyOiNHjhQODg7i2LFjIi4uTrRv377ElydrKvbk5GRx4cIFsX79egFAnDhxQly4cEE8fvxYp2NPT08XXl5ewt3dXdy+fVuhHV3/3P/9918xf/58cf78eXH37l1x5swZ0bNnT2FtbS0ePnyo07EXBKUw9KOJ2G/fvi1mz54tYmJiREJCgjhw4IBwdXUVTZs21fmfGSFebyVgaWkptm/fLm7duiVmzJghTExMSnRukyZ/Zm7duiUkEok4ePBgicWr6dizs7NF/fr1xfvvvy/OnTsnbt++Lb777jshkUjEgQMHSuU5SDXlLlEBUOCxadMmeZ3k5GQxePBgYW9vL0xMTESDBg3E4sWLhUwmk9d5+fKlGDNmjLC2thampqbiww8/FPfu3SsTsc+aNavIdnQx9ryhqoKOhIQEnY79wYMHokuXLqJ69erC0NBQODg4iP79+4vr16+XWNyair2wdks6UdFE7Pfu3RMffPCBsLa2FkZGRqJevXpi7NixJZqUayr2PGFhYcLBwUFUqlRJtG7dWpw8ebLMxD516lTh4OAgpFJpicas6dhv3rwp+vTpI6pXry4qVaokPDw83jlZmLRLIoQQ6vXJEBEREZWMcjdHhYiIiMoPJipERESks5ioEBERkc5iokJEREQ6i4kKERER6SwmKkRERKSzmKgQERGRzmKiQlSExMRESCQS+dbymiaRSLBnz55iXx8ZGQmJRAKJRIJevXq9s66vry/GjRtX7HvRu+X9f6hcubK2QyEqN5iokE4bPHhwkX/5ljRHR0ckJyfDzc0NwP9PDNLS0rQa19tu3LiBzZs3azuMCqGwn8vk5GQsXbq01OMhKs+YqBAVQV9fH7a2tjAwMNB2KO9UvXp1nfiXfE5OjrZD0BpbW1tYWVlpOwyicoWJCpVpUVFRaNmyJYyNjWFnZ4evvvoKubm58vO+vr4YO3YspkyZAmtra9ja2iI0NFShjevXr6Nt27YwMTFBw4YNcezYMYXhmDeHfhITE9GuXTsAQJUqVSCRSDB48GAAr9+G+/a/pps0aaJwv1u3buGDDz6Q3+vNN3TnefDgAfr27YsqVarAxsYGPXv2RGJiosqfzYsXLzBo0CCYm5vDzs4OixcvzlcnOzsbU6ZMQc2aNWFmZgYvLy9ERkYq1Fm/fj0cHR1RqVIl9O7dG0uWLFFIiEJDQ9GkSRP8+OOPqFu3LoyNjSGEwLNnzzB8+HBUr14dlpaWaN++PS5duqTQ9u+//w5PT0+YmJigbt26mD17tsL/v9DQUNSqVQvGxsawt7fH2LFjlXr2op7r8ePH6NevHxwcHFCpUiW4u7vjl19+UWhjx44dcHd3h6mpKWxsbODv748XL14gNDQU4eHh2Lt3r3yo5+3PjIg0R7f/iUj0Dg8ePEDXrl0xePBgbNmyBdevX8ewYcNgYmKikByEh4djwoQJOHfuHKKjozF48GC0adMGHTp0gEwmQ69evVCrVi2cO3cOz58/x8SJEwu9p6OjI3bu3ImPPvoIN27cgKWlJUxNTZWKVyaToU+fPqhatSrOnj2L9PT0fPNFMjMz0a5dO7z//vs4ceIEDAwMMHfuXHTu3BmXL1+GkZGR0p/P5MmTcfz4cezevRu2traYNm0aYmNj0aRJE3mdoKAgJCYmIiIiAvb29ti9ezc6d+6M+Ph4ODs74/Tp0xg5ciS+/fZb9OjRA8eOHcPMmTPz3ev27dv47bffsHPnTujr6wMAunXrBmtra/zxxx+wsrLC2rVr4efnh5s3b8La2hqHDx/GgAEDsHz5crz//vv4999/MXz4cADArFmzsGPHDnz//feIiIhAo0aNkJKSki/RKUxRz/Xq1St4enoiJCQElpaWOHDgAAYOHIi6devCy8sLycnJ6NevHxYuXIjevXvj+fPnOHnyJIQQmDRpEq5du4b09HRs2rQJAGBtba30/xciUpF234lI9G6BgYGiZ8+eBZ6bNm2aaNCggcJbUVetWiXMzc3lb3P18fERbdu2VbiuRYsWIiQkRAghxMGDB4WBgYFITk6Wnz969KjC24MTEhIEAHHhwgUhxP9/y/PTp08V2nVychLff/+9Qlnjxo3FrFmzhBBCHD58WOjr64ukpCT5+YMHDyrca+PGjfmeKSsrS5iamorDhw8X+DkUFM/z58+FkZGRiIiIkJc9fvxYmJqaiuDgYCGEELdv3xYSiUQ8ePBAoT0/Pz8xdepUIYQQffv2Fd26dVM4HxAQIKysrOTfz5o1SxgaGorU1FR52Z9//iksLS3Fq1evFK6tV6+eWLt2rRBCiPfff1/Mnz9f4fxPP/0k7OzshBBCLF68WLi4uIjs7OwCn7swyjxXQbp27SomTpwohBAiNjZWABCJiYkF1n3Xz+WmTZsUPh8iUg97VKjMunbtGlq3bg2JRCIva9OmDTIyMnD//n3UqlULAODh4aFwnZ2dHVJTUwG8noDq6OgIW1tb+fmWLVuWWLy1atWCg4ODvKx169YKdWJjY3H79m1YWFgolL969Qr//vuv0vf6999/kZ2drdC+tbU1GjRoIP8+Li4OQgi4uLgoXJuVlQUbGxsArz+f3r17K5xv2bIl9u/fr1Dm5OSEatWqKTxHRkaGvJ08L1++lD9HbGwsYmJiMG/ePPl5qVSKV69eITMzE5988gmWLl2KunXronPnzujatSu6d+9e5FwhZZ5LKpViwYIF+PXXX/HgwQNkZWUhKysLZmZmAIDGjRvDz88P7u7u6NSpEzp27IiPP/4YVapUeee9iUjzmKhQmSWEUEhS8soAKJQbGhoq1JFIJJDJZIW2UVx6enry++d5c2Lp2+fejhN4PTzk6emJn3/+OV/dNxOBohR0r7fJZDLo6+sjNjZWPlyTx9zcXN5OYZ/xm/L+gn+zbTs7uwLnbuTNb5HJZJg9ezb69OmTr46JiQkcHR1x48YNHD16FMeOHcOoUaOwaNEiREVF5ft/qupzLV68GN9//z2WLl0Kd3d3mJmZYdy4ccjOzgbwegL10aNHcebMGRw5cgQrVqzA9OnTce7cOdSpU6fQexOR5jFRoTKrYcOG2Llzp8JfpmfOnIGFhQVq1qypVBuurq64d+8eHj58iBo1agAAYmJi3nlN3jwRqVSqUF6tWjUkJyfLv09PT0dCQoJCvPfu3cN///0He3t7AEB0dLRCG82aNcOvv/4qn4BaXPXr14ehoSHOnj0r71l6+vQpbt68CR8fHwBA06ZNIZVKkZqaivfff7/AdlxdXfH3338rlJ0/f77I+zdr1gwpKSkwMDBA7dq1C61z48YN1K9fv9B2TE1N0aNHD/To0QOjR4+Gq6sr4uPj0axZs0KvUea5Tp48iZ49e2LAgAEAXic3t27dwnvvvSevI5FI0KZNG7Rp0wZff/01nJycsHv3bkyYMAFGRkb5/v8TUcngqh/Sec+ePcPFixcVjnv37mHUqFFISkrCl19+ievXr2Pv3r2YNWsWJkyYAD095X60O3TogHr16iEwMBCXL1/G6dOnMX36dAD5ezvyODk5QSKRYP/+/fjf//6HjIwMAED79u3x008/4eTJk7hy5QoCAwMV/kXv7++PBg0aYNCgQbh06RJOnjwpv1eegIAAVK1aFT179sTJkyeRkJCAqKgoBAcH4/79+0p/Zubm5hg6dCgmT56MP//8E1euXMHgwYMVPhcXFxcEBARg0KBB2LVrFxISEhATE4Nvv/0Wf/zxBwDgyy+/xB9//IElS5bg1q1bWLt2LQ4ePFhkL5S/vz9at26NXr164fDhw0hMTMSZM2cwY8YMeaLz9ddfY8uWLQgNDcXVq1dx7do1/Prrr5gxYwYAYPPmzdi4cSOuXLmCO3fu4KeffoKpqSmcnJzeeW9lnqt+/fryHpNr165hxIgRSElJkbdx7tw5zJ8/H+fPn8e9e/ewa9cu/O9//5MnMrVr18bly5dx48YNPHr0qEIvySYqcVqaG0OklMDAQAEg3xEYGCiEECIyMlK0aNFCGBkZCVtbWxESEiJycnLk1/v4+Mgnj+bp2bOn/HohhLh27Zpo06aNMDIyEq6uruL3338XAMShQ4eEEPkn0wohxJw5c4Stra2QSCTytp49eyY+/fRTYWlpKRwdHcXmzZsVJtMKIcSNGzdE27ZthZGRkXBxcRGHDh1SmEwrhBDJycli0KBBomrVqsLY2FjUrVtXDBs2TDx79qzAz6iwyb3Pnz8XAwYMEJUqVRI1atQQCxcuzPd5ZGdni6+//lrUrl1bGBoaCltbW9G7d29x+fJleZ1169aJmjVrClNTU9GrVy8xd+5cYWtrKz8/a9Ys0bhx43xxpaeniy+//FLY29sLQ0ND4ejoKAICAsS9e/fkdQ4dOiS8vb2FqampsLS0FC1bthTr1q0TQgixe/du4eXlJSwtLYWZmZlo1aqVOHbsWIGfwduKeq7Hjx+Lnj17CnNzc1G9enUxY8YMMWjQIPkE2X/++Ud06tRJVKtWTRgbGwsXFxexYsUKefupqamiQ4cOwtzcXAAQx48fl5/jZFoizZIIocRgNlEFcvr0abRt2xa3b99GvXr1tB1OkSIjI9GuXTs8ffq0VDZ8GzZsGK5fv46TJ0+W+L3Kos2bN2PcuHE6t3MxUVnFOSpU4e3evRvm5uZwdnbG7du3ERwcjDZt2pSJJOVNDg4O6N69e76Ny9T13XffoUOHDjAzM8PBgwcRHh6O1atXa/Qe5YW5uTlyc3NhYmKi7VCIyg0mKlThPX/+HFOmTEFSUhKqVq0Kf3//Andx1VVeXl64desWgP+/qkWT/v77byxcuBDPnz9H3bp1sXz5cnz++ecav4+yTp48iS5duhR6Pm/OkDbkvbjy7dVGRFR8HPohojLl5cuXePDgQaHn37WKiIjKHiYqREREpLO4PJmIiIh0FhMVIiIi0llMVIiIiEhnMVEhIiIincVEhYiIiHQWExUiIiLSWUxUiIiISGcxUSEiIiKd9f8ADxt3jrbaiK8AAAAASUVORK5CYII=", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "pmp_spring.sel(window=1).plot()" ] }, { "cell_type": "code", - "execution_count": 11, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "" - ] - }, - "execution_count": 11, - "metadata": {}, - "output_type": "execute_result" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAioAAAHFCAYAAADcytJ5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAABwJ0lEQVR4nO3de1yO9/8H8NfdOakIqchdUhplzocyFTnNme8YIYdhYyNnOcaQMdbYmNPCMOZsfB03NYeQcog5hCJWa4skUer+/P7o2/1z63Tf3Xfdt3o997ge6/5c1/W53tfdrd59TpdECCFAREREpIP0tB0AERERUWGYqBAREZHOYqJCREREOouJChEREeksJipERESks5ioEBERkc5iokJEREQ6i4kKERER6SwmKkRERKSzmKi8JSgoCBKJBP/++6/G6hw2bBgcHBxKPZ7Fixdj//79JbqOOtLS0rBo0SJ4e3vDxsYGlStXhru7O7766iu8evVK6Xp27NiBxo0bw8TEBHZ2dggICEB6enqx523atAkSiUS+afJ7t3//foW6L126pLG6i5P3vdekvPcqPj5eo/Vq040bNzB27Fi0adMGZmZmkEgkCAsL03ZYZWL79u0ICQkpcJ9EIkFQUJD8dVhYmM6/NydPnkSbNm1QqVIlVK9eHcOGDUNycrJS5z5//hzjx49HrVq1YGxsDBcXFyxduhQ5OTmlHDWVNiYqZWDOnDnYt29fqV9HW4nKw4cPERISgqZNm2LdunU4ePAg/vOf/yAoKAjdu3eHMk9p2LZtGwYOHIgWLVrgyJEjmDdvHjZt2oS+ffsqHcfevXsRERGBKlWqqHE3iry8vBAREYHZs2drrE5lffLJJ4iIiCjz675rLl26hP3798PKygodOnTQdjhlqqhEJSIiAp988knZBqSG8PBwdO3aFTVr1sSBAwfw7bff4uTJk+jQoQMyMzOLPDc7OxsdO3bE1q1bMXPmTBw6dAg9evTAjBkzMHHixDK6AyotBtoOoCJwcnLSdgilytHREfHx8TAzM5OXtW/fHmZmZpg6dSrOnj2Ltm3bFnp+Tk4Opk6dik6dOmH9+vUAAB8fH5ibm8PPzw9HjhxB165di42jSZMmJW65KkzVqlXRunVr3Lp1S6P1KqN27dqoXbt2mV/3XTNkyBD4+/sDAHbv3o1ff/1VyxHphtatW2s7BJVMnToVLi4u2L17NwwMcn81OTo6wtPTEz/++CM+++yzQs/dvXs3Lly4gD179sj/uOnYsSPS09Px/fffY9y4cahfv36Z3AdpHltUCvH3339j4MCBsLS0RM2aNTFixAg8e/ZM4RghBFavXo3GjRvD1NQUVatWxX/+8x/cv39f4biCun5SU1MxcuRIWFlZoXLlyujWrRvu37+fr7lW2XgkEglevHiBzZs3y7spvL29NfV2FMnMzEwhScnTsmVLAEBCQkKR558/fx6JiYkYPny4QvlHH32EypUrq9Ua5e3tDTc3N0RERMDDwwOmpqZwcHBAaGgoAODw4cNo2rQpKlWqBHd3dxw9erTE1yqIEAI1a9bEuHHj5GU5OTmoWrUq9PT08Pfff8vLV6xYAQMDA6SmpgIouOvHwcEB3bt3x9GjR9G0aVOYmprC1dUVP/74Y75rnz9/Hp6envKutMDAQLx+/TrfcTKZDEuXLoWrqyuMjY1hbW2NoUOH4tGjR/Jjvv/+e+jp6Sk0wy9fvhwSiUTh3mQyGapWrYrJkyer/maVkJ5e6fwYO3z4MBo3bgxjY2M4Ojri66+/zvc9iY+Ph0QiwaZNm/Kd//a/5bt372L48OFwdnZGpUqVUKtWLfTo0QMxMTEK5+V10fz888+YNWsW7OzsYGFhAV9fX9y+fVt+nLe3Nw4fPowHDx4odE8Wdv3CXLp0CT179oSVlRVMTEzQpEkT/PLLL8q/URrw+PFjREZGYsiQIfIkBQA8PDzg4uJS7M+As2fPQiKR5PuDpnv37pDJZGXSok2lh4lKIfr16wcXFxfs2bMHM2bMwPbt2/M1IY4ZMwYBAQHw9fXF/v37sXr1aty4cQMeHh4Kv4DeJpPJ0KNHD2zfvh3Tp0/Hvn370KpVK3Tp0qXE8URERMDU1BQffvghIiIiEBERgdWrVxd5jzk5OcjOzi52k8lkSr5rin7//XcAQMOGDYs87vr16wCARo0aKZQbGhrC1dVVvr+kkpKSMHz4cHzyySc4cOAA3N3dMWLECCxYsACBgYGYNm0a9uzZg8qVK6N3797466+/1LremyQSCdq3b4+TJ0/Kyy5duoTU1FSYmJjgt99+k5efPHkSzZo1K7br6urVq5g8eTImTpyIAwcOoFGjRhg5ciT++OMP+TF//vknOnTogNTUVGzatAk//PADLl++jIULF+ar77PPPsP06dPRsWNHHDx4EF9++SWOHj0KDw8P+XgfX19fCCHyxWtqaooTJ07kuzdfX98i70Emkyn12dPW+ILffvsNvXr1grm5OXbs2IFly5bhl19+kSe4JfHXX3+hWrVqWLJkCY4ePYrvv/8eBgYGaNWqlUICkmfmzJl48OABNmzYgHXr1iE2NhY9evSQvyerV6+Gp6cnbGxs5P/mVe0qPHXqFDw9PZGamooffvgBBw4cQOPGjTFgwIACk6+3aepnSGE/A/LKivsZkJWVBT09PRgaGiqUGxsbAwCuXbtW7L2QDhOkYN68eQKAWLp0qUL52LFjhYmJiZDJZEIIISIiIgQAsXz5coXjEhIShKmpqZg2bZq8zN/fX0ilUvnrw4cPCwBizZo1CucGBwcLAGLevHkqxyOEEGZmZsLf31/pe5VKpQJAsdub8Sjr6tWrwtTUVPTp06fYYxctWiQAiMTExHz7OnXqJFxcXIo8PzQ0VAAQcXFx+fZ5eXkJAOLSpUvyspSUFKGvry9MTU3F48eP5eVXrlwRAMTKlSsLvUZkZGSx9/O2DRs2CADi4cOHQgghFi5cKFxdXUXPnj3F8OHDhRBCZGVlCTMzMzFz5kz5eXnf+zdJpVJhYmIiHjx4IC97+fKlsLKyEmPGjJGXDRgwQJiamoqkpCR5WXZ2tnB1dVV4r27evCkAiLFjxypc58KFCwKAQjy1a9cWI0aMEEIIkZmZKczMzMT06dMFAHk8ixYtEoaGhiI9Pb3I98Tf31+pz56Xl1eR9bxt165dAoA4deqUSue9rVWrVsLOzk68fPlSXpaWliasrKwUvidxcXECgAgNDc1XR3H/drKzs0VWVpZwdnYWEydOlJefOnVKABAffvihwvG//PKLACAiIiLkZd26dVP42VLU9fPqffO9cXV1FU2aNBGvX79WOLd79+7C1tZW5OTkFBq/EJr7GbJt27Z895Zn9OjRwsjIqMjzQ0JCBABx+vRphfI5c+YIAKJTp05Fnk+6jWNUCtGzZ0+F140aNcKrV6+QnJyMmjVr4tChQ5BIJBg8eDCys7Plx9nY2OD9998vcmR9eHg4AKB///4K5QMHDkRgYGCJ4imJX3/9tdhBagBgZ2enUr3x8fHo3r077O3tsWHDBqXPK2yGi7ozX2xtbdGsWTP5aysrK1hbW8PBwUHh3t577z0AwIMHD9S63tvyWhdOnjyJ4cOH48SJE+jYsSOcnZ2xdOlSALktYi9evCi2JQIAGjdujDp16shfm5iYwMXFRSHuU6dOoUOHDgqfDX19fQwYMADz589XOA7I7Z58U8uWLfHee+/ht99+w6JFiwAAHTp0kLcMnTt3DhkZGZg0aRI2btyIEydOYOTIkfJZGwV1Bb4pKCgIn3/+ebH3am5uXuwxmvbixQtERkZi7NixMDExUYilR48e2Lx5c4nqzc7OxtKlS7F161bcvXtXoRvu5s2b+Y4v6N88kPv51MT4k7t37+LWrVv4+uuv5fHl+fDDD3Ho0CHcvn1b/u+iIJr+GVLSnwF+fn5YsGABRo8ejdDQUNSvXx9HjhzBypUrAZRe9yCVDSYqhahWrZrC67wmxJcvXwLIHTMi/jf+oCB169YttO6UlBQYGBjAyspKobyohKO4eEqiQYMGSs3IUeUf+YMHD+Dj4wMDAwP89ttv+e6xIHn3lpKSku89ePLkiVJ1FKWg842MjPKVGxkZAYBKU6qVIZVK4eTkhJMnT2LAgAGIiIjA5MmTUa9ePYwfPx63b9+Wd6N4eHgUW9/bnwUg9/Pw5mchJSUFNjY2+Y57uywlJQVAbjL3Njs7O4Xkx9fXF5s3b0ZsbCxOnjyJJk2awNraWt61NWjQIJw7dw6zZs0q9h7q1Kmj1EBhTU/PVsbTp08hk8mUev9UMWnSJHz//feYPn06vLy85OOUPvnkkwL/HZfGv/k35XVPT5kyBVOmTCnwmOKm+mvqZ8ibPwPepszPgOrVq+Po0aPw9/eXJ3HVqlXDihUrMHLkSNSqVavYGEl3MVEpoerVq0MikeD06dPyHyBvKqgsT7Vq1ZCdnZ3vH2BSUlKpxFoYJycnpVoP5s2bp9SgvAcPHsDb2xtCCISFhSk9Y8Xd3R0AEBMTgwYNGsjLs7OzcevWLQwcOFCpenRZhw4dcODAAYSHh0Mmk8Hb2xvm5uaws7PDiRMncPLkSXzwwQdFfm5UUa1atQI/T2+X5f2CSExMzPf9+uuvv1C9enWFewByW4byWoXyymfPno0//vgDmZmZSrUKjRgxQqmWCS8vrzJf96Nq1aqQSCRKvX95LS5vtyoU9At369atGDp0KBYvXqxQ/u+//2p0Sr2y8r63gYGBhS4DUNxMGU39DHFzcwOQ+zPgww8/VNgXExMj31+UFi1a4M8//0R8fDxevHgBZ2dnREVFAQDatWtX7Pmku5iolFD37t2xZMkSPH78OF8XTnG8vLywdOlS7Ny5U2HK3Y4dO9SK6e2/qoujyWbbhw8fwtvbGzk5OQgLC4NUKlU6jlatWsHW1habNm3CgAED5OW7d+9Genq6Smup6CpfX1+sW7cOISEhaN26tbxLo0OHDti3bx8iIyPz/QJTh4+PDw4ePIi///5b3kqVk5ODnTt3KhzXvn17ALm/RFu0aCEvj4yMxM2bNxVaR2xtbdGgQQPs2bMHUVFR8ng7duyIMWPGYMWKFbCwsFCopzC63PVjZmaGli1bYu/evVi2bJk8GXn+/Hm+qc81a9aEiYlJvsGaBw4cyFevRCLJl4gePnwYjx8/Rr169UoUq6r/5t9Uv359ODs74+rVqyX+7GnqZ0itWrXQsmVLbN26FVOmTIG+vj6A3Jlrt2/fRkBAgNIx5c2wFEJg+fLlsLOzw0cffaT0+aR7mKiUkKenJ0aPHo3hw4fj0qVLaNeuHczMzJCYmIgzZ87A3d290Hn/Xbp0gaenJyZPnoy0tDQ0a9YMERER2LJlC4CS96e6u7sjLCwMv/76K2xtbWFubl7kX0R5LRnqSk5Oho+PDxITE7Fx40YkJycrTGN9cz2QBw8ewMnJCf7+/ti4cSOA3LETS5cuxZAhQzBmzBgMHDgQsbGxmDZtGjp27FjkbCht2rRpE4YPH47Q0NB8Yzze1r59e0gkEhw/flxhjIivr698DRBlWiKUNXv2bBw8eBDt27fH3LlzUalSJXz//fd48eKFwnH169fH6NGjsWrVKujp6aFr166Ij4/HnDlzYG9vn2+mW4cOHbBq1SqYmprC09MTQO5aF46Ojjh+/Dh69uypML20MA4ODhpb8yYjIwP//e9/AeT+YgNyx4H9+++/MDMzU5iyOmzYMGzevBlxcXFFXv/LL79Ely5d0LFjR0yePBk5OTn46quvYGZmhidPnsiPyxun9uOPP8LJyQnvv/8+Ll68iO3bt+ers3v37ti0aRNcXV3RqFEjREVFYdmyZWqtlePu7o69e/dizZo1aNasGfT09NC8eXOlz1+7di26du2Kzp07Y9iwYahVqxaePHmCmzdvIjo6Grt27Sr2+pry1VdfoWPHjvjoo48wduxYJCcnY8aMGXBzc1NYuqCgnyEAMGvWLLi7u8PW1hYPHz7Ejz/+iAsXLuDw4cMwNTXVWJxU9pioqGHt2rVo3bo11q5di9WrV0Mmk8HOzg6enp7yNUQKoqenh19//RWTJ0/GkiVLkJWVBU9PT2zduhWtW7cucTPwt99+i3HjxuHjjz9GRkZGmTWb//nnn/K1YwYPHpxv/5vNvkII5OTk5Jt2OnjwYOjr62PJkiXYtGkTrKysMHToUPlATl2Ut7x/QeM73latWjU0btwYly9fVkhI8r7O268pbm5uOHnyJCZPngx/f39UrVoVQ4YMQb9+/TB69GiFY9esWQMnJyds3LgR33//PSwtLdGlSxcEBwfnGyfh6+uLVatWoW3btgoDTX19fbF+/XqNJlvKSk5OzvcXc97nTSqVKjwuID09HaampsX+G+vYsSP279+P2bNnY8CAAbCxscHYsWPx8uVLhUQTyF1PBgCWLl2K9PR0tG/fHocOHcqXCH377bcwNDREcHAw0tPT0bRpU+zdu1etVY8nTJiAGzduYObMmXj27BmEEEqNGcnj4+ODixcvYtGiRQgICMDTp09RrVo1NGjQQOWWYnV5e3vjv//9L+bOnYsePXqgUqVK6N69O5YtW6bQElXYz5CnT59i+vTpSEpKgoWFBby8vHDhwgWNJlOkHRKhyqeaStX27dvh5+eHs2fPKjWoknLltWzcvXsXUqlUqb/olZX3Q3HLli0YOXIkIiMj5X+x9u/fH3FxcYiMjNTY9ah02djYYMiQIVi2bFmJzg8KCsL8+fNVSgaISD1sUdGSn3/+GY8fP4a7uzv09PRw/vx5LFu2DO3atWOSUkJ5/fz//POPwiBQdRw4cAB9+vTJV543YHjr1q0auQ6Vvhs3biAjIwPTp0/XdihEpAImKlqSt+LlwoUL8eLFC9ja2mLYsGEFrhxKRevRo4dCq4YmZ1B4e3sr1J03K0kikSj9VFfSDQ0bNkRaWpq2wyAiFbHrh4iIiHQWl+sjIiIincVEhYiIiHQWExUiIiLSWeV+MK1MJsNff/0Fc3NzrTw3hIiI3h1CCDx//hx2dnal+jDDV69eISsrS+16jIyMFNY0Ko/KfaLy119/wd7eXtthEBHROyQhIUGtVYOL8urVKzhKKyMpOaf4g4thY2ODuLi4cp2slPtEJe9ZIdIZc6FXjr+RpDyzv7QdAemS6usuaDsE0iHZeI0z+G+pPmcqKysLSck5iIuSwsK85K02ac9lcGz2AFlZWUxU3mV53T16JiZMVAgAoG+k7QhIlxhIDLUdAumS/y3YURZDBSzM9dRKVCqKcp+oEBER6aIcIUOOGiuZ5QiZ5oLRYUxUiIiItEAGARlKnqmoc+67hG1OREREpLO0mqgEBQVBIpEobDY2NvL9QggEBQXBzs4Opqam8Pb2xo0bN7QYMRERkWbINPBfRaD1FpWGDRsiMTFRvsXExMj3LV26FCtWrMB3332HyMhI2NjYoGPHjnj+/LkWIyYiIlJfjhBqbxWB1hMVAwMD2NjYyLcaNWoAyG1NCQkJwaxZs9C3b1+4ublh8+bNyMjIwPbt27UcNREREZUFrScqsbGxsLOzg6OjIz7++GPcv38fABAXF4ekpCR06tRJfqyxsTG8vLxw7tw5bYVLRESkEXmDadXZKgKtzvpp1aoVtmzZAhcXF/z9999YuHAhPDw8cOPGDSQlJQEAatasqXBOzZo18eDBg0LrzMzMRGZmpvx1Wlpa6QRPRESkBhkEcjjrp1haTVS6du0q/9rd3R1t2rSBk5MTNm/ejNatWwPIv+iOEKLIhXiCg4Mxf/780gmYiIiIypTWu37eZGZmBnd3d8TGxspn/+S1rORJTk7O18rypsDAQDx79ky+JSQklGrMREREJcGuH+XoVKKSmZmJmzdvwtbWFo6OjrCxscGJEyfk+7OyshAeHg4PD49C6zA2NoaFhYXCRkREpGu0Mevn+fPnCAgIgFQqhampKTw8PBAZGSnfP2zYsHzLhuT1cGiLVrt+pkyZgh49eqBOnTpITk7GwoULkZaWBn9/f0gkEgQEBGDx4sVwdnaGs7MzFi9ejEqVKmHQoEHaDJuIiEhtsv9t6pyvqk8++QTXr1/HTz/9BDs7O2zduhW+vr74888/UatWLQBAly5dEBoaKj/HyEi7D0jTaqLy6NEjDBw4EP/++y9q1KiB1q1b4/z585BKpQCAadOm4eXLlxg7diyePn2KVq1a4fjx46X6VEsiIqLy6OXLl9izZw8OHDiAdu3aAchdeHX//v1Ys2YNFi5cCCC3Z+LNxVe1TauJyo4dO4rcL5FIEBQUhKCgoLIJiIiIqIzkqDnrJ+/ct2e3Ghsbw9jYON/x2dnZyMnJgYmJiUK5qakpzpw5I38dFhYGa2trVKlSBV5eXli0aBGsra1LHKe6dGqMChERUUWRI9TfAMDe3h6WlpbyLTg4uMDrmZubo02bNvjyyy/x119/IScnB1u3bsWFCxeQmJgIIHc27rZt2/D7779j+fLliIyMRPv27RWW/ShrfHoyERHROywhIUFh4khBrSl5fvrpJ4wYMQK1atWCvr4+mjZtikGDBiE6OhoAMGDAAPmxbm5uaN68OaRSKQ4fPoy+ffuW3k0UgYkKERGRFmhqMK0qM1ydnJwQHh6OFy9eIC0tDba2thgwYAAcHR0LPN7W1hZSqRSxsbFqRKoedv0QERFpgQwS5KixyVD44qfFMTMzg62tLZ4+fYpjx46hV69eBR6XkpKChIQE2Nralvha6mKLChERUQVx7NgxCCFQv3593L17F1OnTkX9+vUxfPhwpKenIygoCP369YOtrS3i4+Mxc+ZMVK9eHX369NFazExUiIiItEAmcjd1zlfVs2fPEBgYiEePHsHKygr9+vXDokWLYGhoiOzsbMTExGDLli1ITU2Fra0tfHx8sHPnTq0uC8JEhYiISAvyunDUOV9V/fv3R//+/QvcZ2pqimPHjpU4ntLCMSpERESks9iiQkREpAXaaFF5FzFRISIi0gKZkEAmSp5sqHPuu4SJChERkRawRUU5HKNCREREOostKkRERFqQAz3kqNFekKPBWHQZExUiIiItEGqOUREVZIwKu36IiIhIZ7FFhYiISAs4mFY5TFSIiIi0IEfoIUeoMUZFjeX33yXs+iEiIiKdxRYVIiIiLZBBApka7QUyVIwmFSYqREREWsAxKsph1w8RERHpLLaoEBERaYH6g2nZ9UNERESlJHeMihoPJawgXT9MVIiIiLRApuYS+hVlMC3HqBAREZHOYosKERGRFnCMinKYqBAREWmBDHpcR0UJ7PohIiIincUWFSIiIi3IERLkCDUWfFPj3HcJExUiIiItyFFz1k8Ou36IiIiItIstKkRERFogE3qQqTHrR8ZZP0RERFRa2PWjHJ3p+gkODoZEIkFAQIC87O+//8awYcNgZ2eHSpUqoUuXLoiNjdVekERERFSmdCJRiYyMxLp169CoUSN5mRACvXv3xv3793HgwAFcvnwZUqkUvr6+ePHihRajJSIiUp8M/z/zpySbTNs3UEa0nqikp6fDz88P69evR9WqVeXlsbGxOH/+PNasWYMWLVqgfv36WL16NdLT0/Hzzz9rMWIiIiL15S34ps5WEWj9LseNG4du3brB19dXoTwzMxMAYGJiIi/T19eHkZERzpw5U6YxEhERaVreEvrqbBWBVgfT7tixA9HR0YiMjMy3z9XVFVKpFIGBgVi7di3MzMywYsUKJCUlITExsdA6MzMz5UkOAKSlpZVK7ERERFT6tJaOJSQkYMKECdi6datCq0keQ0ND7NmzB3fu3IGVlRUqVaqEsLAwdO3aFfr6+oXWGxwcDEtLS/lmb29fmrdBRERUIjJI1N4qAq0lKlFRUUhOTkazZs1gYGAAAwMDhIeHY+XKlTAwMEBOTg6aNWuGK1euIDU1FYmJiTh69ChSUlLg6OhYaL2BgYF49uyZfEtISCjDuyIiIlIOu36Uo7W77NChA2JiYnDlyhX51rx5c/j5+eHKlSsKrSaWlpaoUaMGYmNjcenSJfTq1avQeo2NjWFhYaGwEREREfD8+XMEBARAKpXC1NQUHh4eCsMvhBAICgqCnZ0dTE1N4e3tjRs3bmgxYi2OUTE3N4ebm5tCmZmZGapVqyYv37VrF2rUqIE6deogJiYGEyZMQO/evdGpUydthExERKQx6i/4pvq5n3zyCa5fv46ffvoJdnZ22Lp1K3x9ffHnn3+iVq1aWLp0KVasWIFNmzbBxcUFCxcuRMeOHXH79m2Ym5uXOFZ16HS7UWJiIoYMGQJXV1eMHz8eQ4YM4dRkIiIqF2RCovamipcvX2LPnj1YunQp2rVrh3r16iEoKAiOjo5Ys2YNhBAICQnBrFmz0LdvX7i5uWHz5s3IyMjA9u3bS+ldKJ5OLaEfFham8Hr8+PEYP368doIhIiJ6B7w9u9XY2BjGxsb5jsvOzkZOTk6+CSympqY4c+YM4uLikJSUpNBrYWxsDC8vL5w7dw5jxowpnRsohk63qBAREZVXsv91/ZR0y1vwzd7eXmG2a3BwcIHXMzc3R5s2bfDll1/ir7/+Qk5ODrZu3YoLFy4gMTERSUlJAICaNWsqnFezZk35Pm3QqRYVIiKiikL9pyfnnpuQkKAwcaSg1pQ8P/30E0aMGIFatWpBX18fTZs2xaBBgxAdHS0/RiJR7FISQuQrK0tsUSEiInqHvT3TtahExcnJCeHh4UhPT0dCQgIuXryI169fw9HRETY2NgCQr/UkOTk5XytLWWKiQkREpAU5kKi9lZSZmRlsbW3x9OlTHDt2DL169ZInKydOnJAfl5WVhfDwcHh4eGjilkuEXT9ERERaoKmuH1UcO3YMQgjUr18fd+/exdSpU1G/fn0MHz4cEokEAQEBWLx4MZydneHs7IzFixejUqVKGDRoUInjVBcTFSIiIi3IAdRqFckpwTnPnj1DYGAgHj16BCsrK/Tr1w+LFi2CoaEhAGDatGl4+fIlxo4di6dPn6JVq1Y4fvy41tZQAZioEBERVRj9+/dH//79C90vkUgQFBSEoKCgsguqGExUiIiItEAbXT/vIiYqREREWqDugwX5UEIiIiIiLWOLChERkRYISCBTYzCtUOPcdwkTFSIiIi1g149yKsZdEhER0TuJLSpERERaIBMSyETJu2/UOfddwkSFiIhIC/KegqzO+RVBxbhLIiIieiexRYWIiEgL2PWjHCYqREREWiCDHmRqdGyoc+67hIkKERGRFuQICXLUaBVR59x3ScVIx4iIiOidxBYVIiIiLeAYFeUwUSEiItICoebTkwVXpiUiIiLSLraoEBERaUEOJMhR48GC6pz7LmGiQkREpAUyod44E5nQYDA6jF0/REREpLPYokJERKQFMjUH06pz7ruEiQoREZEWyCCBTI1xJuqc+y5hokJERKQFXJlWORWj3YiIiIjeSWxRISIi0gKOUVEOExUiIiItkEHNJfQryBiVipGOERER0TuJLSpERERaINSc9SPYolK2goODIZFIEBAQIC9LT0/H559/jtq1a8PU1BTvvfce1qxZo70giYiINCTv6cnqbBWBTrSoREZGYt26dWjUqJFC+cSJE3Hq1Cls3boVDg4OOH78OMaOHQs7Ozv06tVLS9ESERFRWdF6i0p6ejr8/Pywfv16VK1aVWFfREQE/P394e3tDQcHB4wePRrvv/8+Ll26pKVoiYiINCNv1o86W0Wg9bscN24cunXrBl9f33z72rZti4MHD+Lx48cQQuDUqVO4c+cOOnfuXGh9mZmZSEtLU9iIiIh0Dbt+lKPVrp8dO3YgOjoakZGRBe5fuXIlRo0ahdq1a8PAwAB6enrYsGED2rZtW2idwcHBmD9/fmmFTERERGVIay0qCQkJmDBhArZu3QoTE5MCj1m5ciXOnz+PgwcPIioqCsuXL8fYsWNx8uTJQusNDAzEs2fP5FtCQkJp3QIREVGJ5T3rR51NFdnZ2Zg9ezYcHR1hamqKunXrYsGCBZDJZPJjhg0bBolEorC1bt1a07euEqVaVPr27atyxT/88AOsra0L3R8VFYXk5GQ0a9ZMXpaTk4M//vgD3333HZ49e4aZM2di37596NatGwCgUaNGuHLlCr7++usCu4oAwNjYGMbGxirHS0REVJbU7b5R9dyvvvoKP/zwAzZv3oyGDRvi0qVLGD58OCwtLTFhwgT5cV26dEFoaKj8tZGRUYlj1ASlEpX9+/ejf//+MDU1VarS7du3Iz09vchEpUOHDoiJiVEoGz58OFxdXTF9+nTk5OTg9evX0NNTbPTR19dXyP6IiIjeRWWdqERERKBXr17yP/4dHBzw888/55ugYmxsDBsbmxLHpWlKj1FZuXJlkYnHm3bv3l3sMebm5nBzc1MoMzMzQ7Vq1eTlXl5emDp1KkxNTSGVShEeHo4tW7ZgxYoVyoZNRERUrr09aaSwnoW2bdvihx9+wJ07d+Di4oKrV6/izJkzCAkJUTguLCwM1tbWqFKlCry8vLBo0SKlf/+XBqUSlVOnTsHKykrpSo8cOYJatWqVOKg8O3bsQGBgIPz8/PDkyRNIpVIsWrQIn376qdp1ExERaZOmWlTs7e0VyufNm4egoKB8x0+fPh3Pnj2Dq6sr9PX1kZOTg0WLFmHgwIHyY7p27YqPPvoIUqkUcXFxmDNnDtq3b4+oqCitDatQKlHx8vJSqdKiZuUUJSwsTOG1jY2NQj8ZERFReaGpRCUhIQEWFhby8sISip07d2Lr1q3Yvn07GjZsiCtXriAgIAB2dnbw9/cHAAwYMEB+vJubG5o3bw6pVIrDhw+XaLyqJpRoerJMJsPdu3eRnJycb7xIu3btNBIYERERFc/CwkIhUSnM1KlTMWPGDHz88ccAAHd3dzx48ADBwcHyROVttra2kEqliI2N1WjMqlA5UTl//jwGDRqEBw8eQAihsE8ikSAnJ0djwREREZVXAlDzoYSqycjIUHmCSkpKChISEmBra1uCCDVD5UTl008/RfPmzXH48GHY2tpCIqkYK+MRERFpUlnP+unRowcWLVqEOnXqoGHDhrh8+TJWrFiBESNGAMh9pE1QUBD69esHW1tbxMfHY+bMmahevTr69OlT4jjVpXKiEhsbi927d6NevXqlEQ8RERGVglWrVmHOnDkYO3YskpOTYWdnhzFjxmDu3LkAcltXYmJisGXLFqSmpsLW1hY+Pj7YuXMnzM3NtRa3yolKq1atcPfuXSYqREREaijrFhVzc3OEhITkm46cx9TUFMeOHStxPKVFqUTl2rVr8q+/+OILTJ48GUlJSXB3d4ehoaHCsY0aNdJshEREROVQWScq7yqlEpXGjRtDIpEoDJ7N69MCIN/HwbRERESkSUolKnFxcaUdBxERUYXCFhXlKJWoSKVS+dd//PEHPDw8YGCgeGp2djbOnTuncCwREREVTAgJhBrJhjrnvkv0ij9EkY+PD548eZKv/NmzZ/Dx8dFIUEREROWdDBK1t4pA5UQlbyzK21JSUmBmZqaRoIiIiIgAFaYn563xL5FIMGzYMIVnCeTk5ODatWvw8PDQfIRERETlEMeoKEfpRMXS0hJAbouKubk5TE1N5fuMjIzQunVrjBo1SvMREhERlUMco6IcpROV0NBQCCEghMCqVau0ukodERERVQwqjVERQmD79u1ISkoqrXiIiIgqhLyuH3W2ikClREVPTw/Ozs5ISUkprXiIiIgqhLyuH3W2ikDlWT9Lly7F1KlTcf369dKIh4iIiEhO5YcSDh48GBkZGXj//fdhZGSkMKgWQIFrrOgC8weAvpG2oyBd8MQjS9shkA7pcTVb2yGQDslMz0ZYGU1gFWp231SUFhWVE5XCnrpIREREyhMA3niEXonOrwhUTlT8/f1LIw4iIiKifFROVIDcBd7279+PmzdvQiKRoEGDBujZsyf09fU1HR8REVG5JIMEEjWWwa8oS+irnKjcvXsXH374IR4/foz69etDCIE7d+7A3t4ehw8fhpOTU2nESUREVK5wwTflqDzrZ/z48XByckJCQgKio6Nx+fJlPHz4EI6Ojhg/fnxpxEhERFTucB0V5ajcohIeHo7z58/DyspKXlatWjUsWbIEnp6eGg2OiIiIKjaVExVjY2M8f/48X3l6ejqMjDj/l4iISBlCqDnrp4JM+1G566d79+4YPXo0Lly4IH/2z/nz5/Hpp5+iZ8+epREjERFRucOVaZWjcqKycuVKODk5oU2bNjAxMYGJiQk8PT1Rr149fPvtt6URIxEREVVQKnf9VKlSBQcOHEBsbCxu3boFIQQaNGiAevXqlUZ8RERE5RJn/SinROuoAICzszOcnZ01GQsREVGFIRMSSNRINjjrpxA5OTnYtGkTfvvtNyQnJ0Mmkyns//333zUWHBEREVVsKicqEyZMwKZNm9CtWze4ublBIqkYGR0REZEmcdaPclROVHbs2IFffvkFH374YWnEQ0REVCHkJirqjFHRYDA6TOVZP0ZGRhw4S0RERGVC5URl8uTJ+PbbbyE0nMoFBwdDIpEgICBAXiaRSArcli1bptFrExERlTWuo6Iclbt+zpw5g1OnTuHIkSNo2LAhDA0NFfbv3btX5SAiIyOxbt06NGrUSKE8MTFR4fWRI0cwcuRI9OvXT+VrEBER6RLxv02d83VNdnY2tm3bhs6dO8PGxkYjdZZoHZU+ffpo5OJA7tL7fn5+WL9+PRYuXKiw7+2bPHDgAHx8fFC3bl2NXZ+IiEgbyuM6KgYGBvjss89w8+ZNzdWp6gmhoaFKHXf27Fk0b94cxsbGRR43btw4dOvWDb6+vvkSlTf9/fffOHz4MDZv3lxkfZmZmcjMzJS/TktLUypeIiIiUl+rVq1w5coVSKVSjdRX4gXfitO1a1dcuXKlyNaPHTt2IDo6GpGRkcXWt3nzZpibm6Nv375FHhccHIz58+erHC8REVGZKo99PwDGjh2LSZMmISEhAc2aNYOZmZnC/reHeRRH5cG0yipusG1CQgImTJiArVu3wsTEpNj6fvzxR/j5+RV7bGBgIJ49eybfEhISVIqbiIioTKg7kFbFrp/s7GzMnj0bjo6OMDU1Rd26dbFgwQKFhVuFEAgKCoKdnR1MTU3h7e2NGzduqHSdAQMGIC4uDuPHj4enpycaN26MJk2ayP+vqlJrUSlOVFQUkpOT0axZM3lZTk4O/vjjD3z33XfIzMyEvr4+AOD06dO4ffs2du7cWWy9xsbGxXY3ERERVTRfffUVfvjhB2zevBkNGzbEpUuXMHz4cFhaWmLChAkAgKVLl2LFihXYtGkTXFxcsHDhQnTs2BG3b9+Gubm5UteJi4vTaNxaS1Q6dOiAmJgYhbLhw4fD1dUV06dPlycpALBx40Y0a9YM77//flmHSUREVCrKemXaiIgI9OrVC926dQMAODg44Oeff8alS5f+V59ASEgIZs2aJR9msXnzZtSsWRPbt2/HmDFjlLqOpsam5Cm1rp/imJubw83NTWEzMzNDtWrV4ObmJj8uLS0Nu3btwieffKKtUImIiDROU+uopKWlKWxvTih5U9u2bfHbb7/hzp07AICrV6/izJkz8pXm4+LikJSUhE6dOsnPMTY2hpeXF86dO6fSvf3000/w9PSEnZ0dHjx4AAAICQnBgQMHVH6fSi1R0dQzgHbs2AEhBAYOHKiR+oiIiMoTe3t7WFpayrfg4OACj5s+fToGDhwIV1dXGBoaokmTJggICJD/fk1KSgIA1KxZU+G8mjVryvcpY82aNZg0aRI+/PBDpKamIicnB0Du8iYhISEq31+pdf2UZOXasLCwfGWjR4/G6NGjNRARERGRDinBgNh85yN3coqFhYW8uLBxmjt37sTWrVuxfft2NGzYEFeuXEFAQADs7Ozg7+8vP+7thgYhhEqND6tWrcL69evRu3dvLFmyRF7evHlzTJkyRel68qicqLx8+RJCCFSqVAkA8ODBA+zbtw8NGjRQaC56/vy5ysEQERFVFJoao2JhYaGQqBRm6tSpmDFjBj7++GMAgLu7Ox48eIDg4GD4+/vLF1lNSkqCra2t/Lzk5OR8rSxFiYuLK3B2j7GxMV68eKF0PXlU7vrp1asXtmzZAgBITU1Fq1atsHz5cvTq1Qtr1qxROQAiIiIqfRkZGdDTU/y1r6+vL5+e7OjoCBsbG5w4cUK+PysrC+Hh4fDw8FD6Oo6Ojrhy5Uq+8iNHjqBBgwYqx61yi0p0dDS++eYbAMDu3btRs2ZNXL58GXv27MHcuXPx2WefqRwEERFRhVPGC7716NEDixYtQp06ddCwYUNcvnwZK1aswIgRIwBA/mDgxYsXw9nZGc7Ozli8eDEqVaqEQYMGKX2dqVOnYty4cXj16hWEELh48SJ+/vlnBAcHY8OGDaoFjRIkKhkZGfK51MePH0ffvn2hp6eH1q1by0f2EhERUdHK+lk/q1atwpw5czB27FgkJyfDzs4OY8aMwdy5c+XHTJs2DS9fvsTYsWPx9OlTtGrVCsePH1d6DRUgd6mR7OxsTJs2DRkZGRg0aBBq1aqFb7/9Vt7tpAqVE5V69eph//796NOnD44dO4aJEycCyO3DUqaPjIiIiP6nDJfBNzc3R0hISJEzbyQSCYKCghAUFKTWtUaNGoVRo0bh33//hUwmg7W1dYnrUnmMyty5czFlyhQ4ODigZcuWaNOmDYDc1pWSLI1LRERE5VP16tXVSlKAErSo/Oc//0Hbtm2RmJiosFJshw4d0KdPH7WCISIiqijKuuunrKSkpGDu3Lk4deoUkpOTFZ4lBABPnjxRqb4SraNiY2OD9PR0nDhxAu3atYOpqSlatGihsUXeiIiIyr1y+vTkwYMH4969exg5ciRq1qypdm6gcqKSkpKC/v3749SpU5BIJIiNjUXdunXxySefoEqVKli+fLlaAREREdG768yZMzhz5ozGns+n8hiViRMnwtDQEA8fPpQv+gbkPtb56NGjGgmKiIio/JNoYNM9rq6uePnypcbqU7lF5fjx4zh27Bhq166tUO7s7MzpyURERMoqp10/q1evxowZMzB37ly4ubnB0NBQYb+qM4RVTlRevHih0JKS599//y30+QJERERUMVSpUgXPnj1D+/btFcrznhmU95BCZamcqLRr1w5btmzBl19+CSB3zrVMJsOyZcvg4+OjanVEREQVUzltUfHz84ORkRG2b9+uncG0y5Ytg7e3Ny5duoSsrCxMmzYNN27cwJMnT3D27Fm1giEiIqowNPT0ZF1z/fp1XL58GfXr19dIfSoPpm3QoAGuXbuGFi1aoGPHjnjx4gX69u2Ly5cvw8nJSSNBERER0bupefPmSEhI0Fh9JV5HZcGCBRoLgoiIqKIRIndT53xd9MUXX2DChAmYOnUq3N3d8w2mbdSokUr1lShROX36NNauXYv79+9j165dqFWrFn766Sc4Ojqibdu2JamSiIioYimnY1QGDBgAAPKnMgO541nLbDDtnj17MGTIEPj5+SE6OhqZmZkAgOfPn2Px4sX473//q2qVREREFU85HaMSFxen0fpUTlQWLlyIH374AUOHDsWOHTvk5R4eHuwOIiIiquCkUqlG61M5Ubl9+zbatWuXr9zCwgKpqamaiImIiKjck4jcTZ3zddXjx49x9uzZAh9KOH78eJXqUjlRsbW1xd27d+Hg4KBQfubMGdStW1fV6oiIiCqmcjpGJTQ0FJ9++imMjIxQrVo1hXVUJBJJ6ScqY8aMwYQJE/Djjz9CIpHgr7/+QkREBKZMmYK5c+eqWh0RERGVI3PnzsXcuXMRGBgIPT2VV0HJR+VEZdq0aXj27Bl8fHzw6tUrtGvXDsbGxpgyZQo+//xztQMiIiKqEMrpYNqMjAx8/PHHGklSABUXfMvJyUF4eDgmT56Mf//9FxcvXsT58+fxzz//yJfUJyIiIiUIDWw6aOTIkdi1a5fG6lOpRUVfXx+dO3fGzZs3YWVlhebNm2ssECIiInr3BQcHo3v37jh69GiBC76tWLFCpfpU7vpxd3fH/fv34ejoqOqpRERElKecDqZdvHgxjh07Jn/Wz9uDaVWlcqKyaNEiTJkyBV9++SWaNWsGMzMzhf0WFhYqB0FERFThlNNEZcWKFfjxxx8xbNgwjdSncqLSpUsXAEDPnj0VMqOSLo1LRERE5YexsTE8PT01Vp/KicqpU6c0dnEiIqIKq5zO+pkwYQJWrVqFlStXaqQ+lRMVLy8vjVyYiIioIiuvK9NevHgRv//+Ow4dOoSGDRvmG0y7d+9elepTOVG5du1ageUSiQQmJiaoU6cOjI2NVa2WiIioYimnY1SqVKmCvn37aqw+lROVxo0bFzlq19DQEAMGDMDatWthYmKiVnBERET0bgkNDdVofSovG7dv3z44Oztj3bp1uHLlCi5fvox169ahfv362L59OzZu3Ijff/8ds2fP1migREREVPGUaHryt99+i86dO8vLGjVqhNq1a2POnDm4ePEizMzMMHnyZHz99ddK1xscHIyZM2diwoQJCAkJkZffvHkT06dPR3h4OGQyGRo2bIhffvkFderUUTV0IiIinSGBmmNUNBaJZjk6OhbZ83L//n2V6lM5UYmJiYFUKs1XLpVKERMTAyC3eygxMVHpOiMjI7Fu3To0atRIofzevXto27YtRo4cifnz58PS0hI3b95klxIREZGOCggIUHj9+vVrXL58GUePHsXUqVNVrk/lRMXV1RVLlizBunXrYGRkJA9iyZIlcHV1BQA8fvwYNWvWVKq+9PR0+Pn5Yf369Vi4cKHCvlmzZuHDDz/E0qVL5WV169ZVNWQiIiLdU46nJxfk+++/x6VLl1SuT+UxKt9//z0OHTqE2rVrw9fXFx07dkTt2rVx6NAhrFmzBkBus87YsWOVqm/cuHHo1q0bfH19FcplMhkOHz4MFxcXdO7cGdbW1mjVqhX2799fZH2ZmZlIS0tT2IiIiHROOX0oYWG6du2KPXv2qHyeyomKh4cH4uPjsWDBAjRq1Ahubm5YsGAB4uLi0Lp1awDAkCFDlGre2bFjB6KjoxEcHJxvX3JyMtLT07FkyRJ06dIFx48fR58+fdC3b1+Eh4cXWmdwcDAsLS3lm729vaq3SEREVO44ODhAIpHk28aNGwcAGDZsWL59eb/XNWH37t2wsrJS+TyVu34AoHLlyvj0009LcqpcQkICJkyYgOPHjxc45kQmkwEAevXqhYkTJwLIHfty7tw5/PDDD4UuPBcYGIhJkybJX6elpTFZISIi3VPG66hERkYqPObm+vXr6NixIz766CN5WZcuXRSmF+cN8VBFkyZN8j1iJykpCf/88w9Wr16tcn0lSlR++uknrF27Fvfv30dERASkUim++eYb1K1bF7169VKqjqioKCQnJ6NZs2byspycHPzxxx/47rvv8OLFCxgYGKBBgwYK57333ns4c+ZMofUaGxtzwTkiItJ5Zb0ybY0aNRReL1myBE5OTgp/+BsbG8PGxqbkQQHo3bu3wms9PT3UqFED3t7e8rGsqlA5UVmzZg3mzp2LgIAALFy4UJ6dVa1aFSEhIUonKh06dJDPEsozfPhwuLq6Yvr06TA2NkaLFi1w+/ZthWPu3LlT4KwjIiKiiujtsZjK/MGelZWFrVu3YtKkSQqtH2FhYbC2tkaVKlXg5eWFRYsWwdraWqV45s2bp9LxxVF5jMqqVauwfv16zJo1CwYG/5/nNG/ePF/iURRzc3O4ubkpbGZmZqhWrRrc3NwAAFOnTsXOnTuxfv163L17F9999x1+/fVXpQfqEhER6SwNDaa1t7dXGJtZ0LjPt+3fvx+pqakYNmyYvKxr167Ytm0bfv/9dyxfvhyRkZFo3749MjMzVbqt6OhohXzgwIED6N27N2bOnImsrCyV6gJK0KISFxeHJk2a5Cs3NjbGixcvVA6gKH369MEPP/yA4OBgjB8/HvXr18eePXvQtm1bjV6HiIiozGlojEpCQgIsLCzkxcoMf9i4cSO6du0KOzs7edmAAQPkX7u5uaF58+aQSqU4fPiwSs/uGTNmDGbMmAF3d3fcv38fAwYMQN++fbFr1y5kZGQoLOqqDJUTFUdHR1y5ciVf98uRI0fyjSdRVVhYWL6yESNGYMSIEWrVS0REpGs0NUbFwsJCIVEpzoMHD3Dy5Mlin2Jsa2sLqVSK2NhYleK6c+cOGjduDADYtWsXvLy8sH37dpw9exYff/xx6ScqU6dOxbhx4/Dq1SsIIXDx4kX8/PPPCA4OxoYNG1StjoiIiMpQaGgorK2t0a1btyKPS0lJQUJCAmxtbVWqXwghn7l78uRJdO/eHUBuF9W///6rcrwqJyrDhw9HdnY2pk2bhoyMDAwaNAi1atXCt99+i48//ljlAIiIiCokLaxMK5PJEBoaCn9/f4Vxpunp6QgKCkK/fv1ga2uL+Ph4zJw5E9WrV0efPn1Uukbz5s2xcOFC+Pr6Ijw8XL4YbFxcnNKr1r+pRNOTR40ahVGjRuHff/+FTCZTeUQwERFRhVfG66gAuS0cDx8+zDekQl9fHzExMdiyZQtSU1Nha2sLHx8f7Ny5E+bm5ipdIyQkBH5+fti/fz9mzZqFevXqAchd8M3Dw0PlmEuUqOSpXr26OqcTERFRGerUqROEyJ/hmJqa4tixYxq5RqNGjQqcBbxs2TLo6+vLX//888/o2bMnzMzMiqxPqUTl7VXmihIdHa3UcURERBVZWS/4pm1vr0I/ZswYtGrVqtiHDSuVqLy5ytyrV6+wevVqNGjQAG3atAEAnD9/Hjdu3OD6JkRERMrSQtePLimoZacgSiUqb64y98knn2D8+PH48ssv8x2TkJCgQohERERERVN5Zdpdu3Zh6NCh+coHDx5cosc3ExERVUji/7t/SrK96y0qylI5UTE1NS3woYBnzpwp8CnIREREVAANLaFf3qk86ycgIACfffYZoqKi0Lp1awC5Y1R+/PFHzJ07V+MBEhERUcWlcqIyY8YM1K1bF99++y22b98OAHjvvfewadMm9O/fX+MBEhERlUsVfDCtVCqFoaFhsceVaB2V/v37MykhIiJSQ3mfnnzp0iXcvHkTEokErq6uaN68ucL+69evK1WPWgu+EREREb3p0aNHGDhwIM6ePYsqVaoAAFJTU+Hh4YGff/4Z9vb2KtWn1GBaKysrlR4kVKdOHTx48EClQIiIiOjdN2LECLx+/Ro3b97EkydP8OTJE9y8eRNCCIwcOVLl+pRqUUlNTcWRI0dgaWmpVKUpKSnIyclRORgiIqIKo5yOUTl9+jTOnTuH+vXry8vq16+PVatWwdPTU+X6lO768ff3V7lyIiIiKlh5HaNSp04dvH79Ol95dnY2atWqpXJ9SnX9yGQylbfi1u4nIiKi8mfp0qX44osvcOnSJfky+ZcuXcKECRPw9ddfq1wfB9MSERFpi462iqhj2LBhyMjIQKtWrWBgkJtmZGdnw8DAACNGjMCIESPkxz558qTY+pioEBERaUM5HaMSEhKi0fqYqBAREZHGaHpMKxMVIiIiLSivg2nzJCcnIzk5GTKZTKG8UaNGKtXDRIWIiEgbymnXT1RUFPz9/eVrp7xJIpGovHxJiRKVe/fuITQ0FPfu3cO3334La2trHD16FPb29mjYsGFJqiQiIqJyYPjw4XBxccHGjRtRs2ZNSCQStepTOVEJDw9H165d4enpiT/++AOLFi2CtbU1rl27hg0bNmD37t1qBURERFQRlNeun7i4OOzduxf16tXTSH1KraPyphkzZmDhwoU4ceIEjIyM5OU+Pj6IiIjQSFBERETlntDApoM6dOiAq1evaqw+lVtUYmJisH379nzlNWrUQEpKikaCIiIiKvfK6RiVDRs2wN/fH9evX4ebmxsMDQ0V9vfs2VOl+lROVKpUqYLExEQ4OjoqlF++fLlES+MSERFR+XHu3DmcOXMGR44cybevJINpVe76GTRoEKZPn46kpCRIJBLIZDKcPXsWU6ZMwdChQ1WtjoiIqELKG6OizqaLxo8fjyFDhiAxMTHf43VK8sBilROVRYsWoU6dOqhVqxbS09PRoEEDtGvXDh4eHpg9e7bKARAREVVI5XSMSkpKCiZOnIiaNWtqpD6Vu34MDQ2xbds2LFiwAJcvX4ZMJkOTJk3g7OyskYCIiIjo3dW3b1+cOnUKTk5OGqmvxAu+OTk5aSwIIiKiCqecDqZ1cXFBYGAgzpw5A3d393yDacePH69SfUolKpMmTVK6whUrVqgUABERUUVUXtdR2bBhAypXrozw8HCEh4cr7JNIJKWTqFy+fFnhdVRUFHJyclC/fn0AwJ07d6Cvr49mzZqpdHEiIiIqX+Li4uRf5y2hr87qtEoNpj116pR869GjB7y9vfHo0SNER0cjOjoaCQkJ8PHxQbdu3UocSHBwMCQSCQICAuRlw4YNg0QiUdhat25d4msQERHpjHI6mBYANm7cCDc3N5iYmMDExARubm7YsGFDiepSeYzK8uXLcfz4cVStWlVeVrVqVSxcuBCdOnXC5MmTVQ4iMjIS69atK/CJil26dEFoaKj89Zur4RIREb2rymvXz5w5c/DNN9/giy++QJs2bQAAERERmDhxIuLj47Fw4UKV6lM5UUlLS8Pff/+d7+GDycnJeP78uarVIT09HX5+fli/fn2BwRsbG8PGxkbleomIiKjsrVmzBuvXr8fAgQPlZT179kSjRo3wxRdfqJyoqLyOSp8+fTB8+HDs3r0bjx49wqNHj7B7926MHDkSffv2VbU6jBs3Dt26dYOvr2+B+8PCwmBtbQ0XFxeMGjUKycnJRdaXmZmJtLQ0hY2IiEjnlNOun5ycHDRv3jxfebNmzZCdna1yfSonKj/88AO6deuGwYMHQyqVQiqVws/PD127dsXq1atVqmvHjh2Ijo5GcHBwgfu7du2Kbdu24ffff8fy5csRGRmJ9u3bIzMzs9A6g4ODYWlpKd/s7e1ViomIiKhMlHGi4uDgkG/cp0Qiwbhx43LDEQJBQUGws7ODqakpvL29cePGDZVva/DgwVizZk2+8nXr1sHPz0/l+lTu+qlUqRJWr16NZcuW4d69exBCoF69ejAzM1OpnoSEBEyYMAHHjx+HiYlJgccMGDBA/rWbmxuaN28OqVSKw4cPF9p6ExgYqDCdOi0tjckKERHpHMn/NnXOV0VkZKTCEvbXr19Hx44d8dFHHwEAli5dihUrVmDTpk1wcXHBwoUL0bFjR9y+fRvm5uYqXWvjxo04fvy4fALM+fPnkZCQgKFDhyr8jlZmSZMSL/hmZmZW4OBXZUVFRSE5OVlhSnNOTg7++OMPfPfdd8jMzIS+vr7COba2tpBKpYiNjS20XmNjYxgbG5c4LiIiovKoRo0aCq+XLFkCJycneHl5QQiBkJAQzJo1S94QsHnzZtSsWRPbt2/HmDFjlL7O9evX0bRpUwDAvXv35NeuUaMGrl+/Lj9O2SnLKicqPj4+RVb++++/K1VPhw4dEBMTo1A2fPhwuLq6Yvr06fmSFCD3+QEJCQmwtbVVLWgiIiJdo6GVad8ei6nMH+xZWVnYunUrJk2aBIlEgvv37yMpKQmdOnVSqMfLywvnzp1TKVE5deqU8vegBJUTlcaNGyu8fv36Na5cuYLr16/D399f6XrMzc3h5uamUGZmZoZq1arBzc0N6enpCAoKQr9+/WBra4v4+HjMnDkT1atXR58+fVQNm4iISKdoanry28Mb5s2bh6CgoCLP3b9/P1JTUzFs2DAAQFJSEgDke5BgzZo18eDBg5IHqQEqJyrffPNNgeVBQUFIT09XO6A8+vr6iImJwZYtW5CamgpbW1v4+Phg586dKveVERERlVcJCQmwsLCQv1Zm+MPGjRvRtWtX2NnZKZS/3WMihFBrVVlNKPEYlbcNHjwYLVu2xNdff13iOsLCwuRfm5qa4tixYxqIjIiISAdpqOvHwsJCIVEpzoMHD3Dy5Ens3btXXpa3XllSUpLC8Irk5OR8rSxlTeXpyYWJiIgodPYOERERFUALa6iEhobC2tpa4bE3jo6OsLGxwYkTJ+RlWVlZCA8Ph4eHR8kvpgEqt6i8PS1YCIHExERcunQJc+bM0VhgREREpFkymQyhoaHw9/eHgcH/pwB5z9pbvHgxnJ2d4ezsjMWLF6NSpUoYNGiQFiMuQaJiYWGh0F+lp6eH+vXrY8GCBQqjhYmIiKhw2njWz8mTJ/Hw4UOMGDEi375p06bh5cuXGDt2LJ4+fYpWrVrh+PHjWh8XqnKismnTplIIg4iIqILR0BgVVXTq1AlCFHyiRCJBUFBQsTOGyprKY1Tq1q2LlJSUfOWpqamoW7euRoIiIiIiAkrQohIfH6+wBG+ezMxMPH78WCNBERERlXfa6Pp5FymdqBw8eFD+9bFjx2BpaSl/nZOTg99++w0ODg4aDY6IiKjc0kLXz7tI6USld+/eAHL7sN5egdbQ0BAODg5Yvny5RoMjIiIqr9iiohylExWZTAYgd651ZGQkqlevXmpBEREREQElGKMSFxdXGnGUuvDADbAwz/+gQ6p4Pvh8tLZDIB1y+Gw7bYdAOiQn6xWAw2VzMXb9KEWpRGXlypUYPXo0TExMsHLlyiKPHT9+vEYCIyIiKteYqChFqUTlm2++gZ+fH0xMTAp9KCGQO36FiQoRERFpilKJypvdPe9q1w8REZEu4WBa5ai84NuCBQuQkZGRr/zly5dYsGCBRoIiIiIq99R5IKG63UbvEJUTlfnz5yM9PT1feUZGBubPn6+RoIiIiIiAEsz6EUIoPJQwz9WrV2FlZaWRoIiIiMo7iRCQFPLcHWXPrwiUTlSqVq0KiUQCiUQCFxcXhWQlJycH6enp+PTTT0slSCIionKHs36UonSiEhISAiEERowYgfnz5yssoW9kZAQHBwe0adOmVIIkIiKiiknpRCVv2XxHR0d4eHjA0NCw1IIiIiIq7zjrRzkqj1Hx8vKSf/3y5Uu8fv1aYb+FhYX6UREREZV37PpRisqzfjIyMvD555/D2toalStXRtWqVRU2IiIiKl5ei4o6W0WgcqIydepU/P7771i9ejWMjY2xYcMGzJ8/H3Z2dtiyZUtpxEhEREQVlMpdP7/++iu2bNkCb29vjBgxAh988AHq1asHqVSKbdu2wc/PrzTiJCIiKl/Y9aMUlVtUnjx5AkdHRwC541GePHkCAGjbti3++OMPzUZHRERUTrHrRzkqJyp169ZFfHw8AKBBgwb45ZdfAOS2tFSpUkWTsREREVEFp3KiMnz4cFy9ehUAEBgYKB+rMnHiREydOlXjARIREZVLfNaPUlQeozJx4kT51z4+Prh16xYuXboEJycnvP/++xoNjoiIqDyrKN036lC5ReVtderUQd++fWFlZYURI0ZoIiYiIiIiABpIVPI8efIEmzdv1lR1RERE5ZsQ6m8VgMpdP0RERKQ+LqGvHI21qBARERFpGltUiIiItIELvilF6USlb9++Re5PTU1VNxYiIqIKQyLL3dQ5vyJQuuvH0tKyyE0qlWLo0KElDiQ4OBgSiQQBAQEF7h8zZgwkEglCQkJKfA0iIiKdwXVUlKJ0i0poaGipBREZGYl169ahUaNGBe7fv38/Lly4ADs7u1KLgYiIiHSP1gfTpqenw8/PD+vXr0fVqlXz7X/8+DE+//xzbNu2DYaGhlqIkIiISPP4rB/laD1RGTduHLp16wZfX998+2QyGYYMGYKpU6eiYcOGStWXmZmJtLQ0hY2IiEjnaGEdlcePH2Pw4MGoVq0aKlWqhMaNGyMqKkq+f9iwYZBIJApb69atNXnXKtPqrJ8dO3YgOjoakZGRBe7/6quvYGBggPHjxytdZ3BwMObPn6+pEImIiMqFp0+fwtPTEz4+Pjhy5Aisra1x7969fA8U7tKli8JwDyMjozKOVJHWEpWEhARMmDABx48fh4mJSb79UVFR+PbbbxEdHQ2JRKJ0vYGBgZg0aZL8dVpaGuzt7TUSMxERkaaU9YJvX331Fezt7RWSEAcHh3zHGRsbw8bGpuSBaZjWun6ioqKQnJyMZs2awcDAAAYGBggPD8fKlSthYGCAsLAwJCcno06dOvL9Dx48wOTJkwt8Y/MYGxvDwsJCYSMiItI5Gpr18/Zwh8zMzAIvd/DgQTRv3hwfffQRrK2t0aRJE6xfvz7fcWFhYbC2toaLiwtGjRqF5ORkTd61yrTWotKhQwfExMQolA0fPhyurq6YPn06bG1t0blzZ4X9nTt3xpAhQzB8+PCyDJWIiEhnvd1rMG/ePAQFBeU77v79+1izZg0mTZqEmTNn4uLFixg/fjyMjY3ly4t07doVH330EaRSKeLi4jBnzhy0b98eUVFRMDY2LovbyUdriYq5uTnc3NwUyszMzFCtWjV5ebVq1RT2GxoawsbGBvXr1y+zOImIiEqDprp+EhISFHoPCksoZDIZmjdvjsWLFwMAmjRpghs3bmDNmjXyRGXAgAHy493c3NC8eXNIpVIcPny42IVfS4vWZ/0QERFVSBqa9fP2cIfCEhVbW1s0aNBAoey9997Dw4cPCw3R1tYWUqkUsbGxmrtvFenUs37CwsKK3B8fH18mcRAREZU3np6euH37tkLZnTt3IJVKCz0nJSUFCQkJsLW1Le3wCsUWFSIiIi0o6wXfJk6ciPPnz2Px4sW4e/cutm/fjnXr1mHcuHEAchdgnTJlCiIiIhAfH4+wsDD06NED1atXR58+fUrhHVCOTrWoEBERVRhl/PTkFi1aYN++fQgMDMSCBQvg6OiIkJAQ+Pn5AQD09fURExODLVu2IDU1Fba2tvDx8cHOnTthbm6uRqDqYaJCRESkBWW9jgoAdO/eHd27dy9wn6mpKY4dO1bygEoJu36IiIhIZ7FFhYiISBtkIndT5/wKgIkKERGRNpTxGJV3Fbt+iIiISGexRYWIiEgLJFBzMK3GItFtTFSIiIi04Y3VZUt8fgXArh8iIiLSWWxRISIi0gJtrKPyLmKiQkREpA2c9aMUdv0QERGRzmKLChERkRZIhIBEjQGx6pz7LmGiQkREpA2y/23qnF8BMFEhIiLSAraoKIdjVIiIiEhnsUWFiIhIGzjrRylMVIiIiLSBK9MqhV0/REREpLPYokJERKQFXJlWOUxUiIiItIFdP0ph1w8RERHpLLaoEBERaYFElrupc35FwESFiIhIG9j1oxR2/RAREZHOYosKERGRNnDBN6UwUSEiItICPutHOUxUiIiItIFjVJTCMSpERESks9iiQkREpA0CgDpTjCtGgwoTFSIiIm3gGBXlsOuHiIiIdJbOJCrBwcGQSCQICAiQlwUFBcHV1RVmZmaoWrUqfH19ceHCBe0FSUREpCkC/z+gtkSbtm+gbOhEohIZGYl169ahUaNGCuUuLi747rvvEBMTgzNnzsDBwQGdOnXCP//8o6VIiYiINEStJEXNGUPvEK0nKunp6fDz88P69etRtWpVhX2DBg2Cr68v6tati4YNG2LFihVIS0vDtWvXtBQtERERlSWtJyrjxo1Dt27d4OvrW+RxWVlZWLduHSwtLfH+++8XelxmZibS0tIUNiIiIp0j08CmosePH2Pw4MGoVq0aKlWqhMaNGyMqKkq+XwiBoKAg2NnZwdTUFN7e3rhx44YaN6k+rSYqO3bsQHR0NIKDgws95tChQ6hcuTJMTEzwzTff4MSJE6hevXqhxwcHB8PS0lK+2dvbl0boREREasmb9aPOpoqnT5/C09MThoaGOHLkCP78808sX74cVapUkR+zdOlSrFixAt999x0iIyNhY2ODjh074vnz5xq+e+VpLVFJSEjAhAkTsHXrVpiYmBR6nI+PD65cuYJz586hS5cu6N+/P5KTkws9PjAwEM+ePZNvCQkJpRE+ERGResp4jMpXX30Fe3t7hIaGomXLlnBwcECHDh3g5OT0v3AEQkJCMGvWLPTt2xdubm7YvHkzMjIysH379tJ4B5SitUQlKioKycnJaNasGQwMDGBgYIDw8HCsXLkSBgYGyMnJAQCYmZmhXr16aN26NTZu3AgDAwNs3Lix0HqNjY1hYWGhsBEREZVXbw93yMzMLPC4gwcPonnz5vjoo49gbW2NJk2aYP369fL9cXFxSEpKQqdOneRlxsbG8PLywrlz50r9PgqjtUSlQ4cOiImJwZUrV+Rb8+bN4efnhytXrkBfX7/A84QQhX4TiIiI3hkaalGxt7dXGPJQ2HCK+/fvY82aNXB2dsaxY8fw6aefYvz48diyZQsAICkpCQBQs2ZNhfNq1qwp36cNWluZ1tzcHG5ubgplZmZmqFatGtzc3PDixQssWrQIPXv2hK2tLVJSUrB69Wo8evQIH330kZaiJiIi0hANPZQwISFBoffA2Ni4wMNlMhmaN2+OxYsXAwCaNGmCGzduYM2aNRg6dKj8OIlE8tZlRL6ysqT1WT+F0dfXx61bt9CvXz+4uLige/fu+Oeff3D69Gk0bNhQ2+ERERHphLeHOxSWqNja2qJBgwYKZe+99x4ePnwIALCxsQGAfK0nycnJ+VpZypJOPesnLCxM/rWJiQn27t2rvWCIiIhKkwyAOg0VKk5P9vT0xO3btxXK7ty5A6lUCgBwdHSEjY0NTpw4gSZNmgDIXRokPDwcX331lRqBqkenEhUiIqKKoqwfSjhx4kR4eHhg8eLF6N+/Py5evIh169Zh3bp1ufX97zE2ixcvhrOzM5ydnbF48WJUqlQJgwYNKnGc6mKiQkREVAG0aNEC+/btQ2BgIBYsWABHR0eEhITAz89Pfsy0adPw8uVLjB07Fk+fPkWrVq1w/PhxmJubay1uJipERETaoKHBtKro3r07unfvXuh+iUSCoKAgBAUFlTwuDWOiQkREpA0yAUjUSFRkfCghERERkVaxRYWIiEgbtND18y5iokJERKQVaiYqYKJCREREpYUtKkrhGBUiIiLSWWxRISIi0gaZgFrdNxVk1g8TFSIiIm0QstxNnfMrAHb9EBERkc5iiwoREZE2cDCtUpioEBERaQPHqCiFXT9ERESks9iiQkREpA3s+lEKExUiIiJtEFAzUdFYJDqNXT9ERESks9iiQkREpA3s+lEKExUiIiJtkMkAqLFom6xiLPjGRIWIiEgb2KKiFI5RISIiIp3FFhUiIiJtYIuKUpioEBERaQNXplUKu36IiIhIZ7FFhYiISAuEkEGIks/cUefcdwkTFSIiIm0QQr3umwoyRoVdP0RERKSz2KJCRESkDULNwbQVpEWFiQoREZE2yGSARI1xJhVkjAq7foiIiEhnsUWFiIhIG9j1oxQmKkRERFogZDIINbp+Ksr0ZJ3p+gkODoZEIkFAQAAA4PXr15g+fTrc3d1hZmYGOzs7DB06FH/99Zd2AyUiItKEvCX01dkqAJ1IVCIjI7Fu3To0atRIXpaRkYHo6GjMmTMH0dHR2Lt3L+7cuYOePXtqMVIiIiIqS1pPVNLT0+Hn54f169ejatWq8nJLS0ucOHEC/fv3R/369dG6dWusWrUKUVFRePjwoRYjJiIi0gCZUH9TQVBQECQSicJmY2Mj3z9s2LB8+1u3bq3pu1aZ1seojBs3Dt26dYOvry8WLlxY5LHPnj2DRCJBlSpVCj0mMzMTmZmZ8tdpaWmaCpWIiEhzhACgzvRk1bt+GjZsiJMnT8pf6+vrK+zv0qULQkND5a+NjIxKHp+GaDVR2bFjB6KjoxEZGVnssa9evcKMGTMwaNAgWFhYFHpccHAw5s+fr8kwiYiIygUDAwOFVpS3GRsbF7lfG7TW9ZOQkIAJEyZg69atMDExKfLY169f4+OPP4ZMJsPq1auLPDYwMBDPnj2TbwkJCZoMm4iISCOETKi9Abk9B29ub/YqvC02NhZ2dnZwdHTExx9/jPv37yvsDwsLg7W1NVxcXDBq1CgkJyeX6nugDK0lKlFRUUhOTkazZs1gYGAAAwMDhIeHY+XKlTAwMEBOTg6A3CSlf//+iIuLw4kTJ4psTQFys0ELCwuFjYiISOcImfobAHt7e1haWsq34ODgAi/XqlUrbNmyBceOHcP69euRlJQEDw8PpKSkAAC6du2Kbdu24ffff8fy5csRGRmJ9u3bF5n4lAWtdf106NABMTExCmXDhw+Hq6srpk+fDn19fXmSEhsbi1OnTqFatWpaipaIiEg3JSQkKPxRbmxsXOBxXbt2lX/t7u6ONm3awMnJCZs3b8akSZMwYMAA+X43Nzc0b94cUqkUhw8fRt++fUvvBoqhtUTF3Nwcbm5uCmVmZmaoVq0a3NzckJ2djf/85z+Ijo7GoUOHkJOTg6SkJACAlZWVTgzwISIiKikhExCSkq+FIv43mLakvQdmZmZwd3dHbGxsgfttbW0hlUoL3V9WtD7rpzCPHj3CwYMHAQCNGzdW2Hfq1Cl4e3uXfVBERESaImRQb9aPeivTZmZm4ubNm/jggw8K3J+SkoKEhATY2tqqdR116VSiEhYWJv/awcFBni2qI6+OtPSKsdQwFS/79Stth0A6JCdLv/iDqMLIycr9+aCJ3z/FycZrtR71k43XKh0/ZcoU9OjRA3Xq1EFycjIWLlyItLQ0+Pv7Iz09HUFBQejXrx9sbW0RHx+PmTNnonr16ujTp0/Jg9QAnUpUSsPz588BANKm8doNhHTIHG0HQEQ67vnz57C0tCyVuo2MjGBjY4MzSf9Vuy4bGxulh0I8evQIAwcOxL///osaNWqgdevWOH/+PKRSKV6+fImYmBhs2bIFqampsLW1hY+PD3bu3Alzc3O141SHRJRF2qhFMpkMf/31F8zNzSGRSBT2paWlwd7ePt9ApHcBY9cOxq4djF07KmLsQgg8f/4cdnZ20NMrvYmxr169QlZWltr1GBkZFbvEx7uu3Leo6OnpoXbt2kUe8y5PY2bs2sHYtYOxa0dFi720WlLeZGJiUu4TDE3R+rN+iIiIiArDRIWIiIh0VoVOVIyNjTFv3rxCF8fRZYxdOxi7djB27WDspAvK/WBaIiIiendV6BYVIiIi0m1MVIiIiEhnMVEhIiIincVEhYiIiHRWuUtUgoOD0aJFC5ibm8Pa2hq9e/fG7du3FY5JT0/H559/jtq1a8PU1BTvvfce1qxZo3BMZmYmvvjiC1SvXh1mZmbo2bMnHj169E7Evm7dOnh7e8PCwgISiQSpqamlGremYn/y5Am++OIL1K9fH5UqVUKdOnUwfvx4PHv2TOdjB4AxY8bAyckJpqamqFGjBnr16oVbt269E7HnEUKga9eukEgk2L9//zsRu7e3NyQSicL28ccfvxOxA0BERATat28PMzMzVKlSBd7e3nj58qVOxx4fH5/vPc/bdu3apdOxA0BSUhKGDBkCGxsbmJmZoWnTpti9e3epxU1qEuVM586dRWhoqLh+/bq4cuWK6Natm6hTp45IT0+XH/PJJ58IJycncerUKREXFyfWrl0r9PX1xf79++XHfPrpp6JWrVrixIkTIjo6Wvj4+Ij3339fZGdn63zs33zzjQgODhbBwcECgHj69GmpxazJ2GNiYkTfvn3FwYMHxd27d8Vvv/0mnJ2dRb9+/XQ+diGEWLt2rQgPDxdxcXEiKipK9OjRQ9jb278Tn5k8K1asEF27dhUAxL59+0otbk3G7uXlJUaNGiUSExPlW2pq6jsR+7lz54SFhYUIDg4W169fF3fu3BG7du0Sr1690unYs7OzFd7vxMREMX/+fGFmZiaeP3+u07ELIYSvr69o0aKFuHDhgrh375748ssvhZ6enoiOji612Knkyl2i8rbk5GQBQISHh8vLGjZsKBYsWKBwXNOmTcXs2bOFEEKkpqYKQ0NDsWPHDvn+x48fCz09PXH06NGyCVyULPY3nTp1qswSlbepG3ueX375RRgZGYnXr1+XWqxv01TsV69eFQDE3bt3Sy3Wt6kT+5UrV0Tt2rVFYmJimSQqbytp7F5eXmLChAllFWaBShp7q1ativwMlQVNfd4bN24sRowYUWpxFqSksZuZmYktW7YoHGNlZSU2bNhQugFTiZS7rp+35XUbWFlZycvatm2LgwcP4vHjxxBC4NSpU7hz5w46d+4MAIiKisLr16/RqVMn+Tl2dnZwc3PDuXPndDp2XaGp2J89ewYLCwsYGJTdY6k0EfuLFy8QGhoKR0dH2Nvbl0ncQMljz8jIwMCBA/Hdd9/BxsamzOJ9kzrv+7Zt21C9enU0bNgQU6ZMkT81XZdjT05OxoULF2BtbQ0PDw/UrFkTXl5eOHPmjM7H/raoqChcuXIFI0eOLJOY85Q09rZt22Lnzp148uQJZDIZduzYgczMTHh7e5dp/KQkraVIZUAmk4kePXqItm3bKpRnZmaKoUOHCgDCwMBAGBkZKWTX27ZtE0ZGRvnq69ixoxg9enSpxy1EyWN/k7ZaVDQRuxBC/Pvvv6JOnTpi1qxZpR2ynLqxf//998LMzEwAEK6urmXamqJO7KNHjxYjR46Uv0YZt6ioE/u6devEiRMnRExMjPj555+Fg4OD8PX11fnYIyIiBABhZWUlfvzxRxEdHS0CAgKEkZGRuHPnjk7H/rbPPvtMvPfee6UdrgJ1Yk9NTRWdO3eWH2NhYSGOHz9eluGTCsp1ojJ27FghlUpFQkKCQvmyZcuEi4uLOHjwoLh69apYtWqVqFy5sjhx4oQQovBExdfXV4wZM0anY3+TthIVTcT+7Nkz0apVK9GlSxeRlZVVVqGrHXtqaqq4c+eOCA8PFz169BBNmzYVL1++1OnYDxw4IOrVq6cwtqCsExVNfGbyXLp0SQAQUVFRpR22EKLksZ89e1YAEIGBgQrnubu7ixkzZuh07G/KyMgQlpaW4uuvvy6TmPOoE/vnn38uWrZsKU6ePCmuXLkigoKChKWlpbh27VqZ3gMpp9wmKp9//rmoXbu2uH//vkJ5RkaGMDQ0FIcOHVIoHzlypOjcubMQQojffvtNABBPnjxROKZRo0Zi7ty5pRu4UC/2N2kjUdFE7GlpaaJNmzaiQ4cOZfZLXgjNve95MjMzRaVKlcT27dtLJd43qRP7hAkThEQiEfr6+vINgNDT0xNeXl46HXtBZDJZvjFmpUWd2O/fvy8AiJ9++knhmP79+4tBgwaVbuBCc+/7li1bhKGhoUhOTi7VeN+kTux3794VAMT169cVjunQoUOZ/SFKqil3Y1SEEPj888+xd+9e/P7773B0dFTY//r1a7x+/Rp6eoq3rq+vD5lMBgBo1qwZDA0NceLECfn+xMREXL9+HR4eHjodu7ZoKva0tDR06tQJRkZGOHjwIExMTN6Z2AurOzMzU+Mxv1m/urHPmDED165dw5UrV+QbAHzzzTcIDQ3V6dgLcuPGDbx+/Rq2tralEjegmdgdHBxgZ2eXb3rtnTt3IJVKdTr2N23cuBE9e/ZEjRo1Si3mPJqIPSMjAwB08ucoFUJbGVJp+eyzz4SlpaUICwtTmDqXkZEhP8bLy0s0bNhQnDp1Sty/f1+EhoYKExMTsXr1avkxn376qahdu7Y4efKkiI6OFu3bty/16cmaij0xMVFcvnxZrF+/XgAQf/zxh7h8+bJISUnR6djT0tJEq1athLu7u7h7965CPbr+vt+7d08sXrxYXLp0STx48ECcO3dO9OrVS1hZWYm///5bp2MvCMqg60cTsd+9e1fMnz9fREZGiri4OHH48GHh6uoqmjRpovOfGSFylxKwsLAQu3btErGxsWL27NnCxMSkVMc2afIzExsbKyQSiThy5Eipxavp2LOyskS9evXEBx98IC5cuCDu3r0rvv76ayGRSMThw4fL5D5INeUuUQFQ4BYaGio/JjExUQwbNkzY2dkJExMTUb9+fbF8+XIhk8nkx7x8+VJ8/vnnwsrKSpiamoru3buLhw8fvhOxz5s3r9h6dDH2vK6qgra4uDidjv3x48eia9euwtraWhgaGoratWuLQYMGiVu3bpVa3JqKvbB6SztR0UTsDx8+FO3atRNWVlbCyMhIODk5ifHjx5dqUq6p2PMEBweL2rVri0qVKok2bdqI06dPvzOxBwYGitq1a4ucnJxSjVnTsd+5c0f07dtXWFtbi0qVKolGjRoVOViYtEsihBDqtckQERERlY5yN0aFiIiIyg8mKkRERKSzmKgQERGRzmKiQkRERDqLiQoRERHpLCYqREREpLOYqBAREZHOYqJCVIz4+HhIJBL50vKaJpFIsH///hKfHxYWBolEAolEgt69exd5rLe3NwICAkp8LSpa3vehSpUq2g6FqNxgokI6bdiwYcX+8i1t9vb2SExMhJubG4D/TwxSU1O1Gtfbbt++jU2bNmk7jAqhsM9lYmIiQkJCyjweovKMiQpRMfT19WFjYwMDAwNth1Ika2trnfhL/vXr19oOQWtsbGxgaWmp7TCIyhUmKvROCw8PR8uWLWFsbAxbW1vMmDED2dnZ8v3e3t4YP348pk2bBisrK9jY2CAoKEihjlu3bqFt27YwMTFBgwYNcPLkSYXumDe7fuLj4+Hj4wMAqFq1KiQSCYYNGwYg92m4b/813bhxY4XrxcbGol27dvJrvfmE7jyPHz/GgAEDULVqVVSrVg29evVCfHy8yu/NixcvMHToUFSuXBm2trZYvnx5vmOysrIwbdo01KpVC2ZmZmjVqhXCwsIUjlm/fj3s7e1RqVIl9OnTBytWrFBIiIKCgtC4cWP8+OOPqFu3LoyNjSGEwLNnzzB69GhYW1vDwsIC7du3x9WrVxXq/vXXX9GsWTOYmJigbt26mD9/vsL3LygoCHXq1IGxsTHs7Owwfvx4pe69uPtKSUnBwIEDUbt2bVSqVAnu7u74+eefFerYvXs33N3dYWpqimrVqsHX1xcvXrxAUFAQNm/ejAMHDsi7et5+z4hIc3T7T0SiIjx+/Bgffvghhg0bhi1btuDWrVsYNWoUTExMFJKDzZs3Y9KkSbhw4QIiIiIwbNgweHp6omPHjpDJZOjduzfq1KmDCxcu4Pnz55g8eXKh17S3t8eePXvQr18/3L59GxYWFjA1NVUqXplMhr59+6J69eo4f/480tLS8o0XycjIgI+PDz744AP88ccfMDAwwMKFC9GlSxdcu3YNRkZGSr8/U6dOxalTp7Bv3z7Y2Nhg5syZiIqKQuPGjeXHDB8+HPHx8dixYwfs7Oywb98+dOnSBTExMXB2dsbZs2fx6aef4quvvkLPnj1x8uRJzJkzJ9+17t69i19++QV79uyBvr4+AKBbt26wsrLCf//7X1haWmLt2rXo0KED7ty5AysrKxw7dgyDBw/GypUr8cEHH+DevXsYPXo0AGDevHnYvXs3vvnmG+zYsQMNGzZEUlJSvkSnMMXd16tXr9CsWTNMnz4dFhYWOHz4MIYMGYK6deuiVatWSExMxMCBA7F06VL06dMHz58/x+nTpyGEwJQpU3Dz5k2kpaUhNDQUAGBlZaX094WIVKTdZyISFc3f31/06tWrwH0zZ84U9evXV3gq6vfffy8qV64sf5qrl5eXaNu2rcJ5LVq0ENOnTxdCCHHkyBFhYGAgEhMT5ftPnDih8PTguLg4AUBcvnxZCPH/T3l++vSpQr1SqVR88803CmXvv/++mDdvnhBCiGPHjgl9fX2RkJAg33/kyBGFa23cuDHfPWVmZgpTU1Nx7NixAt+HguJ5/vy5MDIyEjt27JCXpaSkCFNTUzFhwgQhhBB3794VEolEPH78WKG+Dh06iMDAQCGEEAMGDBDdunVT2O/n5ycsLS3lr+fNmycMDQ1FcnKyvOy3334TFhYW4tWrVwrnOjk5ibVr1wohhPjggw/E4sWLFfb/9NNPwtbWVgghxPLly4WLi4vIysoq8L4Lo8x9FeTDDz8UkydPFkIIERUVJQCI+Pj4Ao8t6nMZGhqq8P4QkXrYokLvrJs3b6JNmzaQSCTyMk9PT6Snp+PRo0eoU6cOAKBRo0YK59na2iI5ORlA7gBUe3t72NjYyPe3bNmy1OKtU6cOateuLS9r06aNwjFRUVG4e/cuzM3NFcpfvXqFe/fuKX2te/fuISsrS6F+Kysr1K9fX/46OjoaQgi4uLgonJuZmYlq1aoByH1/+vTpo7C/ZcuWOHTokEKZVCpFjRo1FO4jPT1dXk+ely9fyu8jKioKkZGRWLRokXx/Tk4OXr16hYyMDHz00UcICQlB3bp10aVLF3z44Yfo0aNHsWOFlLmvnJwcLFmyBDt37sTjx4+RmZmJzMxMmJmZAQDef/99dOjQAe7u7ujcuTM6deqE//znP6hatWqR1yYizWOiQu8sIYRCkpJXBkCh3NDQUOEYiUQCmUxWaB0lpaenJ79+njcHlr697+04gdzuoWbNmmHbtm35jn0zEShOQdd6m0wmg76+PqKiouTdNXkqV64sr6ew9/hNeb/g36zb1ta2wLEbeeNbZDIZ5s+fj759++Y7xsTEBPb29rh9+zZOnDiBkydPYuzYsVi2bBnCw8PzfU9Vva/ly5fjm2++QUhICNzd3WFmZoaAgABkZWUByB1AfeLECZw7dw7Hjx/HqlWrMGvWLFy4cAGOjo6FXpuINI+JCr2zGjRogD179ij8Mj137hzMzc1Rq1YtpepwdXXFw4cP8ffff6NmzZoAgMjIyCLPyRsnkpOTo1Beo0YNJCYmyl+npaUhLi5OId6HDx/ir7/+gp2dHQAgIiJCoY6mTZti586d8gGoJVWvXj0YGhri/Pnz8palp0+f4s6dO/Dy8gIANGnSBDk5OUhOTsYHH3xQYD2urq64ePGiQtmlS5eKvX7Tpk2RlJQEAwMDODg4FHrM7du3Ua9evULrMTU1Rc+ePdGzZ0+MGzcOrq6uiImJQdOmTQs9R5n7On36NHr16oXBgwcDyE1uYmNj8d5778mPkUgk8PT0hKenJ+bOnQupVIp9+/Zh0qRJMDIyyvf9J6LSwVk/pPOePXuGK1euKGwPHz7E2LFjkZCQgC+++AK3bt3CgQMHMG/ePEyaNAl6esp9tDt27AgnJyf4+/vj2rVrOHv2LGbNmgUgf2tHHqlUColEgkOHDuGff/5Beno6AKB9+/b46aefcPr0aVy/fh3+/v4Kf9H7+vqifv36GDp0KK5evYrTp0/Lr5XHz88P1atXR69evXD69GnExcUhPDwcEyZMwKNHj5R+zypXroyRI0di6tSp+O2333D9+nUMGzZM4X1xcXGBn58fhg4dir179yIuLg6RkZH46quv8N///hcA8MUXX+C///0vVqxYgdjYWKxduxZHjhwpthXK19cXbdq0Qe/evXHs2DHEx8fj3LlzmD17tjzRmTt3LrZs2YKgoCDcuHEDN2/exM6dOzF79mwAwKZNm7Bx40Zcv34d9+/fx08//QRTU1NIpdIir63MfdWrV0/eYnLz5k2MGTMGSUlJ8jouXLiAxYsX49KlS3j48CH27t2Lf/75R57IODg44Nq1a7h9+zb+/fffCj0lm6jUaWlsDJFS/P39BYB8m7+/vxBCiLCwMNGiRQthZGQkbGxsxPTp08Xr16/l53t5eckHj+bp1auX/HwhhLh586bw9PQURkZGwtXVVfz6668CgDh69KgQIv9gWiGEWLBggbCxsRESiURe17Nnz0T//v2FhYWFsLe3F5s2bVIYTCuEELdv3xZt27YVRkZGwsXFRRw9elRhMK0QQiQmJoqhQ4eK6tWrC2NjY1G3bl0xatQo8ezZswLfo8IG9z5//lwMHjxYVKpUSdSsWVMsXbo03/uRlZUl5s6dKxwcHIShoaGwsbERffr0EdeuXZMfs27dOlGrVi1hamoqevfuLRYuXChsbGzk++fNmyfef//9fHGlpaWJL774QtjZ2QlDQ0Nhb28v/Pz8xMOHD+XHHD16VHh4eAhTU1NhYWEhWrZsKdatWyeEEGLfvn2iVatWwsLCQpiZmYnWrVuLkydPFvgevK24+0pJSRG9evUSlStXFtbW1mL27Nli6NCh8gGyf/75p+jcubOoUaOGMDY2Fi4uLmLVqlXy+pOTk0XHjh1F5cqVBQBx6tQp+T4OpiXSLIkQSnRmE1UgZ8+eRdu2bXH37l04OTlpO5xihYWFwcfHB0+fPi2TBd9GjRqFW7du4fTp06V+rXfRpk2bEBAQoHMrFxO9qzhGhSq8ffv2oXLlynB2dsbdu3cxYcIEeHp6vhNJyptq166NHj165Fu4TF1ff/01OnbsCDMzMxw5cgSbN2/G6tWrNXqN8qJy5crIzs6GiYmJtkMhKjeYqFCF9/z5c0ybNg0JCQmoXr06fH19C1zFVVe1atUKsbGxAP5/VosmXbx4EUuXLsXz589Rt25drFy5Ep988onGr6Os06dPo2vXroXuzxszpA15D658e7YREZUcu36I6J3y8uVLPH78uND9Rc0iIqJ3DxMVIiIi0lmcnkxEREQ6i4kKERER6SwmKkRERKSzmKgQERGRzmKiQkRERDqLiQoRERHpLCYqREREpLOYqBAREZHO+j8nlEm4yc/ZawAAAABJRU5ErkJggg==", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "pmp_summer.sel(window=1).plot()" ] @@ -311,7 +260,7 @@ }, { "cell_type": "code", - "execution_count": 12, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -375,7 +324,7 @@ }, { "cell_type": "code", - "execution_count": 13, + "execution_count": null, "metadata": {}, "outputs": [], "source": [ @@ -385,30 +334,9 @@ }, { "cell_type": "code", - "execution_count": 14, + "execution_count": null, "metadata": {}, - "outputs": [ - { - "data": { - "text/plain": [ - "" - ] - }, - "execution_count": 14, - "metadata": {}, - "output_type": "execute_result" - }, - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjMAAAGwCAYAAABcnuQpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAAA++UlEQVR4nO3dfVxUZf7/8fcAyo3AeAsMSYjlTUgWZipai2mkZmrpalmmZvnNNDe3R6vpVuDm17tdy9qKXdvWm7Wy2s2y3ZV0K8xCv3hTm4qZP0OjBNm8YdASEc7vD2LWEVAGGc4ceD0fj3nkXOc653zmOMK765xzHZthGIYAAAAsys/sAgAAAC4FYQYAAFgaYQYAAFgaYQYAAFgaYQYAAFgaYQYAAFgaYQYAAFhagNkFeFt5ebkOHz6ssLAw2Ww2s8sBAAC1YBiGiouLFR0dLT+/C4+9NPowc/jwYcXExJhdBgAAqIO8vDy1b9/+gn0afZgJCwuTVHEwwsPDTa4GAADUhtPpVExMjOv3+IU0+jBTeWopPDycMAMAgMXU5hIRLgAGAACWRpgBAACWRpgBAACW1uivmQEAND1lZWUqLS01uwxcQLNmzeTv718v2yLMAAAaDcMwVFBQoBMnTphdCmqhZcuWioqKuuR54AgzAIBGozLIREREKCQkhMlSfZRhGPrhhx9UWFgoSXI4HJe0PcIMAKBRKCsrcwWZNm3amF0OLiI4OFiSVFhYqIiIiEs65cQFwACARqHyGpmQkBCTK0FtVf5dXer1TYQZAECjwqkl66ivvytOM9VVeZl0KEs6eUQKjZRi+0p+9XNVNgAAqD3CTF3krJMyZknOw/9tC4+WBi+S4oebVxcAAE0Qp5k8lbNOenO8e5CRJGd+RXvOOnPqAgDUm7JyQ1sOHNW7n3+nLQeOqqzcMLukOluxYoVatmxpdhleRZjxRHlZxYiMqvtS/9SW8XhFPwCAJWXsztcNiz7U2Je36pE1n2vsy1t1w6IPlbE736v7LSws1IMPPqjLL79cgYGBioqK0qBBg7Rly5ZL2u6dd96pr776qp6q9E2cZvLEoayqIzJuDMn5XUW/uBsbrCwAQP3I2J2vh1bvrPK/rAVFp/XQ6p1KH9dDgxMubU6UmowaNUqlpaVauXKlOnbsqCNHjuiDDz7QsWPH6rzN0tJSBQcHu26DbqwYmfHEySP12w8A4DPKyg3NfS/nQmPvmvtejldOOZ04cUKffPKJFi1apJtuukmxsbHq1auXZs+eraFDh0qquPMnPT1dQ4YMUXBwsOLi4vTWW2+5tnHw4EHZbDa9+eab6t+/v4KCgrR69eoqp5nS0tJ07bXX6i9/+Ys6dOggu92uu+66S8XFxa4+xcXFuueee9SiRQs5HA49++yz6t+/v2bMmFHvn70+EGY8ERpZv/0AAD4jO/eY8otO17jckJRfdFrZuXUfKalJaGioQkND9c4776ikpKTGfk8++aRGjRqlf//73xo3bpzGjh2rvXv3uvWZNWuWfvGLX2jv3r0aNGhQtds5cOCA3nnnHf3973/X3//+d23atEkLFy50LX/00Uf16aefat26ddq4caM2b96snTt31s+H9QLCjCdi+1bctaSa7ou3SeGXVfQDAFhKYXHNQaYu/TwREBCgFStWaOXKlWrZsqX69eunOXPm6IsvvnDrN3r0aD3wwAPq3Lmznn76afXs2VO///3v3frMmDFDI0eOVFxcnKKjo6vdX3l5uVasWKGEhATdeOONuvfee/XBBx9IqhiVWblypX73u99p4MCBSkhI0PLly1VW5rvXgxJmPOHnX3H7taSqgean94MXMt8MAFhQRFhQvfbz1KhRo3T48GGtW7dOgwYNUmZmpnr06KEVK1a4+iQlJbmtk5SUVGVkpmfPnhfdV4cOHRQWFuZ673A4XM9J+vrrr1VaWqpevXq5ltvtdnXp0qUuH6tBEGY8FT9cGrNKCj/vArDw6Ip25pkBAEvqFddaDnvQhcbe5bAHqVdca6/VEBQUpJSUFD311FPKysrSxIkTlZqaesF1zp9Ft0WLFhfdT7Nmzapso7y8XFLFQyCr225luy8izNRF/HBpxm5pwt+lUa9U/HfGLoIMAFiYv59NqcPiJdU49q7UYfHy92u4xyXEx8fr1KlTrvdbt251W75161Z17dq1Xvd5xRVXqFmzZsrOzna1OZ1O7d+/v173U5+4Nbuu/Py5/RoAGpnBCQ6lj+uhue/luF0MHGUPUuqweK/dln306FGNHj1akyZNUvfu3RUWFqbt27dr8eLFGjFihKvfW2+9pZ49e+qGG27Qq6++quzsbL3yyiv1WktYWJgmTJigX/3qV2rdurUiIiKUmpoqPz8/n33uFWEGAIBzDE5wKCU+Stm5x1RYfFoRYRWnlrw5IhMaGqrevXvr2Wef1YEDB1RaWqqYmBhNnjxZc+bMcfWbO3eu1qxZo6lTpyoqKkqvvvqq4uPj672eZ555RlOmTNFtt92m8PBwzZw5U3l5eQoK8s71QpfKZvjySbB64HQ6ZbfbVVRUpPDwcLPLAQB4yenTp5Wbm6u4uDif/aV7KWw2m9auXavbb7+9wfd96tQpXXbZZVqyZInuv//+etvuhf7OPPn9zcgMAABw89lnn+nLL79Ur169VFRUpN/85jeS5HbKy5cQZgAAQBW/+93vtG/fPjVv3lzXXXedNm/erLZt25pdVrUIMwAAWEBDXhWSmJioHTt2NNj+LhW3ZgMAAEsjzAAAAEsjzAAAAEsjzAAAAEsjzAAAAEsjzAAAAEsjzAAAcL7yMil3s7TrrxX/LS/z+i4LCwv14IMP6vLLL1dgYKCioqI0aNAgbdmyxev7tjrmmQEA4Fw566SMWZLz8H/bwqOlwYuk+OFe2+2oUaNUWlqqlStXqmPHjjpy5Ig++OADHTt2zGv7bEhnzpxR8+bNvbJtRmYAAKiUs056c7x7kJEkZ35Fe846r+z2xIkT+uSTT7Ro0SLddNNNio2NVa9evTR79mwNHTpUBw8elM1m0+eff+62js1mU2ZmpiQpMzNTNptN77//vhITExUcHKwBAwaosLBQ69ev11VXXaXw8HCNHTtWP/zwg2s7/fv31/Tp0zVjxgy1atVKkZGRWrZsmU6dOqX77rtPYWFhuuKKK7R+/Xr3Q5WTo1tvvVWhoaGKjIzUvffeq++//95tuw8//LAeffRRtW3bVikpKV45dhJhBgCACuVlFSMyqm6m3Z/aMh73yimn0NBQhYaG6p133lFJScklbSstLU0vvPCCsrKylJeXpzFjxmjp0qV67bXX9I9//EMbN27U73//e7d1Vq5cqbZt2yo7O1vTp0/XQw89pNGjR6tv377auXOnBg0apHvvvdcVgvLz85WcnKxrr71W27dvV0ZGho4cOaIxY8ZU2W5AQIA+/fRT/fGPf7ykz3UhpoeZ7777TuPGjVObNm0UEhKia6+91m0K5YkTJ8pms7m9+vTpY2LFAIBG6VBW1REZN4bk/K6iXz0LCAjQihUrtHLlSrVs2VL9+vXTnDlz9MUXX3i8rXnz5qlfv35KTEzU/fffr02bNik9PV2JiYm68cYb9fOf/1wfffSR2zrXXHONnnjiCXXq1EmzZ89WcHCw2rZtq8mTJ6tTp0566qmndPToUVc96enp6tGjh+bPn6+uXbsqMTFRf/7zn/XRRx/pq6++cm33yiuv1OLFi9WlSxd17dr10g7SBZgaZo4fP65+/fqpWbNmWr9+vXJycrRkyRK1bNnSrd/gwYOVn5/vev3zn/80p2AAQON18kj99vPQqFGjdPjwYa1bt06DBg1SZmamevTooRUrVni0ne7du7v+HBkZqZCQEHXs2NGtrbCwsMZ1/P391aZNG1199dVu60hyrbdjxw599NFHrhGl0NBQV1g5cOCAa72ePXt6VHtdmXoB8KJFixQTE6Ply5e72jp06FClX+VV3bVRUlLiNkTndDovuU4AQBMQGlm//eogKChIKSkpSklJ0VNPPaUHHnhAqamp2rx5syT3h02WlpZWu41mzZq5/myz2dzeV7aVl5fXuE5169lsNklyrVdeXq5hw4Zp0aJFVfbvcDhcf27RokXNH7YemToys27dOvXs2VOjR49WRESEEhMT9fLLL1fpl5mZqYiICHXu3FmTJ0+ukijPtWDBAtntdtcrJibGmx8BANBYxPatuGtJtho62KTwyyr6NZD4+HidOnVK7dq1k1RxrUqlcy8Gbmg9evTQnj171KFDB1155ZVur4YKMOcyNcx8/fXXSk9PV6dOnfT+++9rypQp+sUvfqFVq1a5+gwZMkSvvvqqPvzwQy1ZskTbtm3TgAEDarxAavbs2SoqKnK98vLyGurjAACszM+/4vZrSVUDzU/vBy+s6FfPjh49qgEDBmj16tX64osvlJubq7feekuLFy/WiBEjFBwcrD59+mjhwoXKycnRxx9/rCeeeKLe66itadOm6dixYxo7dqyys7P19ddfa8OGDZo0aZLKyrw/J8/5TD3NVF5erp49e2r+/PmSpMTERO3Zs0fp6ekaP368JOnOO+909U9ISFDPnj0VGxurf/zjHxo5cmSVbQYGBiowMLBhPgAAoHGJHy6NWVXDPDMLvTbPTGhoqHr37q1nn31WBw4cUGlpqWJiYjR58mTNmTNHkvTnP/9ZkyZNUs+ePdWlSxctXrxYt9xyi1fquZjo6Gh9+umnmjVrlgYNGqSSkhLFxsZq8ODB8vNr+HESm3HuCbgGFhsbq5SUFP3pT39ytaWnp2vevHn67rvvalyvU6dOeuCBBzRr1qyL7sPpdMput6uoqEjh4eH1UjcAwPecPn1aubm5iouLU1BQ0KVtrLys4q6lk0cqrpGJ7euVEZmm7kJ/Z578/jZ1ZKZfv37at2+fW9tXX32l2NjYGtc5evSo8vLy3C4wAgCgXvn5S3E3ml0FasnUa2Z++ctfauvWrZo/f77+3//7f3rttde0bNkyTZs2TZJ08uRJPfbYY9qyZYsOHjyozMxMDRs2TG3bttUdd9xhZukAAMBHmBpmrr/+eq1du1avv/66EhIS9PTTT2vp0qW65557JFXc675r1y6NGDFCnTt31oQJE9S5c2dt2bJFYWFhZpYOAAB8hOkPmrztttt02223VbssODhY77//fgNXBAAArMT0xxkAAFCfTLyvBR6qr78rwgwAoFGonLH23CdCw7dV/l2dPwOxp0w/zQQAQH3w9/dXy5YtXbPEh4SEuKbhh28xDEM//PCDCgsL1bJlS/n7X9pt74QZAECjUfkcvws99ga+o2XLlrV+9uKFEGYAoLFrQhPA2Ww2ORwORURE1PggRviGZs2aXfKITCXCDAA0Zjnrapiaf5HXpub3Bf7+/vX2ixK+jwuAAaCxylknvTnePchIkjO/oj1nnTl1AfWMMAMAjVF5WcWIjKq79fWntozHK/oBFkeYAYDG6FBW1REZN4bk/K6iH2BxhBkAaIxOHqnffoAPI8wAQGMUGlm//QAfRpgBgMYotm/FXUuqadI4mxR+WUU/wOIIMwDQGPn5V9x+LalqoPnp/eCFjXa+GTQthBkAaKzih0tjVknhDvf28OiK9kY8zwyaFibNA4DGLH641HVok5kBGE0TYQYAGjs/fynuRrOrALyG00wAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSCDMAAMDSTA8z3333ncaNG6c2bdooJCRE1157rXbs2OFabhiG0tLSFB0dreDgYPXv31979uwxsWIAAOBLTA0zx48fV79+/dSsWTOtX79eOTk5WrJkiVq2bOnqs3jxYj3zzDN64YUXtG3bNkVFRSklJUXFxcXmFQ4AAHyGzTAMw6ydP/744/r000+1efPmapcbhqHo6GjNmDFDs2bNkiSVlJQoMjJSixYt0oMPPlhlnZKSEpWUlLjeO51OxcTEqKioSOHh4d75IAAAoF45nU7Z7fZa/f42dWRm3bp16tmzp0aPHq2IiAglJibq5Zdfdi3Pzc1VQUGBbrnlFldbYGCgkpOTlZWVVe02FyxYILvd7nrFxMR4/XMAAADzmBpmvv76a6Wnp6tTp056//33NWXKFP3iF7/QqlWrJEkFBQWSpMjISLf1IiMjXcvON3v2bBUVFbleeXl53v0QAADAVAFm7ry8vFw9e/bU/PnzJUmJiYnas2eP0tPTNX78eFc/m83mtp5hGFXaKgUGBiowMNB7RQMAAJ9i6siMw+FQfHy8W9tVV12lb775RpIUFRUlSVVGYQoLC6uM1gAAgKbJ1DDTr18/7du3z63tq6++UmxsrCQpLi5OUVFR2rhxo2v5mTNntGnTJvXt27dBawUAAL7J1NNMv/zlL9W3b1/Nnz9fY8aMUXZ2tpYtW6Zly5ZJqji9NGPGDM2fP1+dOnVSp06dNH/+fIWEhOjuu+82s3QAAOAjTA0z119/vdauXavZs2frN7/5jeLi4rR06VLdc889rj4zZ87Ujz/+qKlTp+r48ePq3bu3NmzYoLCwMBMrBwAAvsLUeWYagif3qQMAAN9gmXlmAAAALhVhBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWJqpYSYtLU02m83tFRUV5Vo+ceLEKsv79OljYsUAAMDXBJhdQLdu3fSvf/3L9d7f399t+eDBg7V8+XLX++bNmzdYbQAAwPeZHmYCAgLcRmPOFxgYeMHlAACgaTP9mpn9+/crOjpacXFxuuuuu/T111+7Lc/MzFRERIQ6d+6syZMnq7Cw8ILbKykpkdPpdHsBAIDGy2YYhmHWztevX68ffvhBnTt31pEjRzRv3jx9+eWX2rNnj9q0aaM33nhDoaGhio2NVW5urp588kmdPXtWO3bsUGBgYLXbTEtL09y5c6u0FxUVKTw83NsfCQAA1AOn0ym73V6r39+mhpnznTp1SldccYVmzpypRx99tMry/Px8xcbGas2aNRo5cmS12ygpKVFJSYnrvdPpVExMDGEGAAAL8STMmH7NzLlatGihq6++Wvv37692ucPhUGxsbI3LpYprbGoatQEAAI2P6dfMnKukpER79+6Vw+GodvnRo0eVl5dX43IAAND0mBpmHnvsMW3atEm5ubn6v//7P/385z+X0+nUhAkTdPLkST322GPasmWLDh48qMzMTA0bNkxt27bVHXfcYWbZAADAh5h6munbb7/V2LFj9f3336tdu3bq06ePtm7dqtjYWP3444/atWuXVq1apRMnTsjhcOimm27SG2+8obCwMDPLBgAAPsSnLgD2Bk8uIAIAAL7Bk9/fPnXNDAAAgKcIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNIIMwAAwNICPF3h0KFD2rBhg0pLS5WcnKxu3bp5oy4AAIBa8SjMfPzxx7r11lv1ww8/VKwcEKCVK1dq7NixXikOAADgYjw6zfTkk0/qpptu0rfffqujR49q0qRJmjlzprdqAwAAuCiPwsyuXbu0YMECRUdHq1WrVlqyZIkOHz6s48eP12nnaWlpstlsbq+oqCjXcsMwlJaWpujoaAUHB6t///7as2dPnfYFAAAaJ4/CzIkTJxQREeF636JFC4WEhOjEiRN1LqBbt27Kz893vXbt2uVatnjxYj3zzDN64YUXtG3bNkVFRSklJUXFxcV13h8AAGhcPL4AOCcnRwUFBa73hmFo7969bgGje/futS8gIMBtNObc7S5dulS//vWvNXLkSEnSypUrFRkZqddee00PPvigp6UDAIBGyOMwM3DgQBmG4dZ22223yWazyTAM2Ww2lZWV1Xp7+/fvV3R0tAIDA9W7d2/Nnz9fHTt2VG5urgoKCnTLLbe4+gYGBio5OVlZWVk1hpmSkhKVlJS43judTg8/IQAAsBKPwkxubm697rx3795atWqVOnfurCNHjmjevHnq27ev9uzZ4xr9iYyMdFsnMjJShw4dqnGbCxYs0Ny5c+u1TgAA4Ls8CjOxsbH1uvMhQ4a4/nz11VcrKSlJV1xxhVauXKk+ffpIkmw2m9s6laM/NZk9e7YeffRR13un06mYmJh6rRsAAPgOj8LMN998U6t+l19+eZ2KadGiha6++mrt379ft99+uySpoKBADofD1aewsLDKaM25AgMDFRgYWKf9AwAA6/EozMTFxbn+XHndzLmjJHW5ZuZcJSUl2rt3r2688UbFxcUpKipKGzduVGJioiTpzJkz2rRpkxYtWlSn7QMAgMbHozBjs9nUvn17TZw4UcOGDVNAgMfXD7t57LHHNGzYMF1++eUqLCzUvHnz5HQ6NWHCBNlsNs2YMUPz589Xp06d1KlTJ82fP18hISG6++67L2m/AACg8fAojXz77bdauXKlVqxYoT/84Q8aN26c7r//fl111VV12vm3336rsWPH6vvvv1e7du3Up08fbd261XVtzsyZM/Xjjz9q6tSpOn78uHr37q0NGzYoLCysTvsDAACNj804/z7rWvrkk0+0fPlyvfXWW4qPj9f999+v+++/X35+vvUgbqfTKbvdrqKiIoWHh5tdDgAAqAVPfn/XOXnccMMNeuWVV7R//36FhIRoypQplzQTMAAAQF3UOcxkZWXpgQceUOfOnXXy5Em9+OKLatmyZT2WBgAAcHEeXTOTn5+vVatWafny5Tp+/LjuueceZWVlqVu3bt6qDwAA4II8njQvOjpaEyZM0PDhw9WsWTOVlZXpiy++cOvnybOZAAAALoVHFwCfe3Fv5bOYqmzwEuaZ8QYuAAYAwHo8+f1t6rOZAAAALpVHYSYiIkKPPfaY3nnnHZWWlurmm2/W888/r7Zt23qrPgAAgAvy6G6mp556SitWrNDQoUN11113aePGjXrooYe8VRsAAMBFeTQy8/bbb+uVV17RXXfdJUkaN26c+vXrp7KyMvn7+3ulQAAAgAvxaGQmLy9PN954o+t9r169FBAQoMOHD9d7YQAAALXh0chMWVmZmjdv7r6BgACdPXu2XosCgIsqL5MOZUknj0ihkVJsX8mPEWKgKfIozBiGoYkTJyowMNDVdvr0aU2ZMkUtWrRwtb399tv1VyEAnC9nnZQxS3KeMyocHi0NXiTFDzevLgCm8CjMTJgwoUrbuHHj6q0YALionHXSm+MlnTfPlTO/on3MKgIN0BB8aHS0zk/NtgomzQMakfIyaWmC+4iMG1vFCM2MXZxyArypAUZHG+Sp2QDQ4A5lXSDISJIhOb+r6AfAOypHR8//t1g5OpqzrsFLIswAsI6TR+q3HwDPlJdVjMicf5pX+m9bxuMV/RoQYQaAdYRG1m8/AJ7x0dFRwgwA64jtW3FeXrYaOtik8Msq+gGofz46OkqYAWAdfv4VFxhKqhpofno/eCEX/wLe4qOjo4QZANYSP7zi9utwh3t7eDS3ZQPe5qOjox7NMwMAPiF+uNR1qM/McQE0GZWjo2+OV0WgOfdCYPNGRwkzAKzJz1+Ku/Hi/QDUr8rR0WrnmVloyugoYQYAAHjGx0ZHCTMAAMBzPjQ6ygXAAADA0ggzAADA0ggzAADA0ggzAADA0ggzAADA0ggzAADA0rg1G95XXuYzcxEAABofwgy8K2ddDbNELuIZOgCAesFpJnhPzrqK53ecG2QkyZlf0Z6zzpy6AACNCmEG3lFeVjEi4/YQsko/tWU8XtEPAHxFeZmUu1na9deK//IzyhI4zQTvOJRVdUTGjSE5v6vo5yPTYQNo4jgtblmMzMA7Th6p334A4E2cFrc0nwkzCxYskM1m04wZM1xtEydOlM1mc3v16dPHvCJRe6GR9dsPALyF0+KW5xOnmbZt26Zly5ape/fuVZYNHjxYy5cvd71v3rx5Q5aGuortWzE868xX9T8gbBXLY/s2dGUA4I7T4pZn+sjMyZMndc899+jll19Wq1atqiwPDAxUVFSU69W6desLbq+kpEROp9PtBRP4+VecZ5Yk2c5b+NP7wQuZbwaA+Tgtbnmmh5lp06Zp6NChuvnmm6tdnpmZqYiICHXu3FmTJ09WYWHhBbe3YMEC2e121ysmJsYbZaM24odLY1ZJ4Q739vDoinYuqAPgCzgtbnmmnmZas2aNdu7cqW3btlW7fMiQIRo9erRiY2OVm5urJ598UgMGDNCOHTsUGBhY7TqzZ8/Wo48+6nrvdDoJNGaKHy51HcoMwAB8F6fFLc+0MJOXl6dHHnlEGzZsUFBQULV97rzzTtefExIS1LNnT8XGxuof//iHRo4cWe06gYGBNQYdmMTPn/PMAHxX5WnxN8er4jT4uYGG0+JWYNppph07dqiwsFDXXXedAgICFBAQoE2bNun5559XQECAysqqXjXucDgUGxur/fv3m1AxAKDR4rS4pZk2MjNw4EDt2rXLre2+++5T165dNWvWLPn7V03AR48eVV5enhwOR5VlAABcEk6LW5ZpYSYsLEwJCQlubS1atFCbNm2UkJCgkydPKi0tTaNGjZLD4dDBgwc1Z84ctW3bVnfccYdJVQMAGjVOi1uST8wzUx1/f3/t2rVLq1at0okTJ+RwOHTTTTfpjTfeUFhYmNnlAQAAH2EzDKO6S7cbDafTKbvdrqKiIoWHh5tdDgAAqAVPfn+bPs8MAADApSDMAAAASyPMAAAASyPMAAAASyPMAAAAS/PZW7MBAPWjrNxQdu4xFRafVkRYkHrFtZa/3/lPswesizADAI1Yxu58zX0vR/lFp11tDnuQUofFa3ACs6mjceA0EwA0Uhm78/XQ6p1uQUaSCopO66HVO5WxO9+kyoD6RZgBgEaorNzQ3PdyVN2sqJVtc9/LUVl5o543FU0EYQYAGqHs3GNVRmTOZUjKLzqt7NxjDVcU4CWEGQBohAqLaw4ydekH+DLCDAA0QhFhQfXaD/BlhBkAaIR6xbWWwx6kmm7AtqnirqZeca0bsizAKwgzANAI+fvZlDosXpKqBJrK96nD4plvBo0CYQYAGqnBCQ6lj+uhKLv7qaQoe5DSx/Vgnhk0GkyaBwCN2OAEh1Lio5gBGI0aYQYAGjl/P5uSrmhjdhmA13CaCQAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBphBgAAWBozAAMA8JOycoNHP1gQYQYAAEkZu/M1970c5ReddrU57EFKHRbPQzl9HKeZAABNXsbufD20eqdbkJGkgqLTemj1TmXszjepMtQGYQYA0KSVlRua+16OjGqWVbbNfS9HZeXV9YAvIMwAAJq07NxjVUZkzmVIyi86rezcYw1XFDxCmAEANGmFxTUHmbr0Q8MjzAAAmrSIsKB67YeGx91MgC8pL5MOZUknj0ihkVJsX8nP3+yqgEatV1xrOexBKig6Xe11MzZJUfaK27ThmwgzgK/IWSdlzJKch//bFh4tDV4kxQ83ry6gkfP3syl1WLweWr1TNskt0FTOMJM6LJ75ZnyYz5xmWrBggWw2m2bMmOFqMwxDaWlpio6OVnBwsPr37689e/aYVyTgLTnrpDfHuwcZSXLmV7TnrDOnLqCJGJzgUPq4Hoqyu59KirIHKX1cD+aZ8XE+MTKzbds2LVu2TN27d3drX7x4sZ555hmtWLFCnTt31rx585SSkqJ9+/YpLCzMpGqBelZeVjEiU+ONoTYp43Gp61BOOQFeNDjBoZT4KGYAtiDTR2ZOnjype+65Ry+//LJatWrlajcMQ0uXLtWvf/1rjRw5UgkJCVq5cqV++OEHvfbaayZWDNSzQ1lVR2TcGJLzu4p+ALzK38+mpCvaaMS1lynpijYEGYswPcxMmzZNQ4cO1c033+zWnpubq4KCAt1yyy2utsDAQCUnJysrq+Yf6iUlJXI6nW4vwKedPFK//QCgiTH1NNOaNWu0c+dObdu2rcqygoICSVJkZKRbe2RkpA4dOlTjNhcsWKC5c+fWb6GAN4VGXryPJ/0AoIkxbWQmLy9PjzzyiFavXq2goJrv3bfZ3If4DMOo0nau2bNnq6ioyPXKy8urt5oBr4jtW3HXkmr6Xtuk8Msq+gEAqjAtzOzYsUOFhYW67rrrFBAQoICAAG3atEnPP/+8AgICXCMylSM0lQoLC6uM1pwrMDBQ4eHhbi+Yq6zc0JYDR/Xu599py4GjPN/kfH7+FbdfS6oaaH56P3ghF/8CQA1MO800cOBA7dq1y63tvvvuU9euXTVr1ix17NhRUVFR2rhxoxITEyVJZ86c0aZNm7Ro0aLqNgkflLE7X3Pfy3F77onDHqTUYfHc6niu+OHSmFU1zDOzkHlmAOACTAszYWFhSkhIcGtr0aKF2rRp42qfMWOG5s+fr06dOqlTp06aP3++QkJCdPfdd5tRMjyUsTtfD63eWeWG44Ki03po9U7mbjhf/PCK26+ZARgAPOIT88zUZObMmfrxxx81depUHT9+XL1799aGDRuYY8YCysoNzX0v50Izp2juezlKiY/i1sdz+flLcTeaXQUAWIrNMIxGfQGD0+mU3W5XUVER1880oC0Hjmrsy1sv2u/1yX2UdEWbBqgIAGAlnvz+Nn2eGTROhcWnL97Jg34AANSEMAOviAir+Xb7uvQDAKAmhBl4Ra+41nLYgy40c4oc9ornngAAcCkIM/AKfz+bUofFS6px5hSlDovn4l8AwCUjzMBrBic4lD6uh6Ls7qeSouxB3JYNAKg3Pn1rNqxvcIJDKfFRys49psLi04oIqzi1xIgMAKC+EGbgdf5+Nm6/BgB4DaeZAACApRFmAACApRFmAACApRFmAACApRFmAACApRFmAACApRFmAACApRFmAACApTFpHgAA8FhZueEzs7sTZgAAgEcydudr7ns5yi867Wpz2IOUOizelOfucZoJAADUWsbufD20eqdbkJGkgqLTemj1TmXszm/wmggzAACgVsrKDc19L0dGNcsq2+a+l6Oy8up6eA9hBgAA1Ep27rEqIzLnMiTlF51Wdu6xhitKhBkAAFBLhcU1B5m69KsvhBkAAFArEWFB9dqvvhBmAABArfSKay2HPUg13YBtU8VdTb3iWjdkWYQZAABQO/5+NqUOi5ekKoGm8n3qsPgGn2+GMAMAAGptcIJD6eN6KMrufiopyh6k9HE9TJlnhknzAACARwYnOJQSH8UMwAAAwLr8/WxKuqKN2WVI4jQTAACwOMIMAACwNMIMAACwNMIMAACwNMIMAACwNMIMAACwNMIMAACwNMIMAACwNFPDTHp6urp3767w8HCFh4crKSlJ69evdy2fOHGibDab26tPnz4mVgwAAHyNqTMAt2/fXgsXLtSVV14pSVq5cqVGjBihzz77TN26dZMkDR48WMuXL3et07x5c1NqBQAAvsnUMDNs2DC39//7v/+r9PR0bd261RVmAgMDFRUVZUZ5AADAAnzmmpmysjKtWbNGp06dUlJSkqs9MzNTERER6ty5syZPnqzCwsILbqekpEROp9PtBQAAGi/Tw8yuXbsUGhqqwMBATZkyRWvXrlV8fLwkaciQIXr11Vf14YcfasmSJdq2bZsGDBigkpKSGre3YMEC2e121ysmJqahPgoAADCBzTAMw8wCzpw5o2+++UYnTpzQ3/72N/3pT3/Spk2bXIHmXPn5+YqNjdWaNWs0cuTIardXUlLiFnacTqdiYmJUVFSk8PBwr30OAABQf5xOp+x2e61+f5t6zYxUcUFv5QXAPXv21LZt2/Tcc8/pj3/8Y5W+DodDsbGx2r9/f43bCwwMVGBgoNfqBQAAvsX000znMwyjxtNIR48eVV5enhwORwNXBQAAfJWpIzNz5szRkCFDFBMTo+LiYq1Zs0aZmZnKyMjQyZMnlZaWplGjRsnhcOjgwYOaM2eO2rZtqzvuuMPMsgGvKSs3lJ17TIXFpxURFqReca3l72czuywA8GmmhpkjR47o3nvvVX5+vux2u7p3766MjAylpKToxx9/1K5du7Rq1SqdOHFCDodDN910k9544w2FhYWZWTbgFRm78zX3vRzlF512tTnsQUodFq/BCYxGAkBNTL8A2Ns8uYAIMEvG7nw9tHqnzv/HWDkmkz6uB4EGQJPiye9vn7tmBmhqysoNzX0vp0qQkeRqm/tejsrKG/X/dwBAnRFmAJNl5x5zO7V0PkNSftFpZecea7iiAMBCCDOAyQqLaw4ydekHAE0NYQYwWURYUL32A4CmhjADmKxXXGs57EGq6QZsmyruauoV17ohywIAyyDMACbz97MpdVjF4zvODzSV71OHxTPfDADUgDAD+IDBCQ6lj+uhKLv7qaQoexC3ZQPARZj+bCYAFQYnOJQSH8UMwADgIcIM4EP8/WxKuqKN2WUAgKVwmgkAAFgaYQYAAFgaYQYAAFgaYQYAAFgaFwDXUVm5wV0nAAD4AMJMHWTsztfc93LcHg7osAcpdVg884EAANDAOM3koYzd+Xpo9c4qTzkuKDqth1bvVMbufJMqAwCgaSLMeKCs3NDc93JkVLOssm3uezkqK6+uBwAA8AbCjAeyc49VGZE5lyEpv+i0snOPNVxRAAA0cYQZDxQW1xxk6tIPAABcOsKMByLCgi7eyYN+AADg0hFmPNArrrUc9iDVdAO2TRV3NfWKa92QZQEA0KQRZjzg72dT6rB4SaoSaCrfpw6LZ74ZoAGUlRvacuCo3v38O205cJQL74EmjHlmPDQ4waH0cT2qzDMTxTwzQINhricA57IZhtGo/3fG6XTKbrerqKhI4eHh9bZdZgAGzFE519P5P7gq//Wlj+tBoAEaAU9+fzMyU0f+fjYlXdHG7DKAJuVicz3ZVDHXU0p8FP9zATQhXDMDwDKY6wlAdQgzACyDuZ4AVIcwA8AymOsJQHUIMwAsg7meAFSHMAPAMpjrCUB1CDMALKVyrqcou/uppCh7ELdlA00Ut2YDsJzBCQ6lxEcx1xMASYQZABbFXE8AKnGaCQAAWBphBgAAWBphBgAAWJqpYSY9PV3du3dXeHi4wsPDlZSUpPXr17uWG4ahtLQ0RUdHKzg4WP3799eePXtMrBgAAPgaU8NM+/bttXDhQm3fvl3bt2/XgAEDNGLECFdgWbx4sZ555hm98MIL2rZtm6KiopSSkqLi4mIzywYAAD7EZhhGdQ+gNU3r1q3129/+VpMmTVJ0dLRmzJihWbNmSZJKSkoUGRmpRYsW6cEHH6x2/ZKSEpWUlLjeO51OxcTE1OoR4gAAwDc4nU7Z7fZa/f72mWtmysrKtGbNGp06dUpJSUnKzc1VQUGBbrnlFlefwMBAJScnKysrq8btLFiwQHa73fWKiYlpiPIBAIBJTA8zu3btUmhoqAIDAzVlyhStXbtW8fHxKigokCRFRka69Y+MjHQtq87s2bNVVFTkeuXl5Xm1fgAAYC7TJ83r0qWLPv/8c504cUJ/+9vfNGHCBG3atMm13GZzn9HTMIwqbecKDAxUYGCg1+oFAAC+xfQw07x5c1155ZWSpJ49e2rbtm167rnnXNfJFBQUyOH477NWCgsLq4zWXEjlJUFOp7MeqwYAAN5U+Xu7Npf2mh5mzmcYhkpKShQXF6eoqCht3LhRiYmJkqQzZ85o06ZNWrRoUa23V3nnE9fOAABgPcXFxbLb7RfsY2qYmTNnjoYMGaKYmBgVFxdrzZo1yszMVEZGhmw2m2bMmKH58+erU6dO6tSpk+bPn6+QkBDdfffdtd5HdHS08vLyFBYWdsHTU3VReadUXl4ed0pdBMeq9jhWtcexqj2OVe1xrGrPm8fKMAwVFxcrOjr6on1NDTNHjhzRvffeq/z8fNntdnXv3l0ZGRlKSUmRJM2cOVM//vijpk6dquPHj6t3797asGGDwsLCar0PPz8/tW/f3lsfQZJck/7h4jhWtcexqj2OVe1xrGqPY1V73jpWFxuRqWRqmHnllVcuuNxmsyktLU1paWkNUxAAALAc02/NBgAAuBSEmUsQGBio1NRUbgWvBY5V7XGsao9jVXscq9rjWNWerxwrn3ucAQAAgCcYmQEAAJZGmAEAAJZGmAEAAJZGmAEAAJZGmKnBggULdP311yssLEwRERG6/fbbtW/fvguuk5+fr7vvvltdunSRn5+fZsyY0TDFmqwux+rtt99WSkqK2rVrp/DwcCUlJen9999voIrNk56eru7du7smmEpKStL69etr7N9Uv1OS58eqqX6nJM+P1SeffKJ+/fqpTZs2Cg4OVteuXfXss882YMXm8fRYnevTTz9VQECArr32Wu8W6aMWLFjgmp2/Jmb9zCLM1GDTpk2aNm2atm7dqo0bN+rs2bO65ZZbdOrUqRrXKSkpUbt27fTrX/9a11xzTQNWa666HKuPP/5YKSkp+uc//6kdO3bopptu0rBhw/TZZ581YOUNr3379lq4cKG2b9+u7du3a8CAARoxYoT27NlTbf+m+p2SPD9WTfU7JXl+rFq0aKGHH35YH3/8sfbu3asnnnhCTzzxhJYtW9bAlTc8T49VpaKiIo0fP14DBw5soEp9y7Zt27Rs2TJ17979gv1M+5lloFYKCwsNScamTZtq1T85Odl45JFHvFuUj/L0WFWKj4835s6d66WqfFerVq2MP/3pTxft15S/U5Vqe6wqNdXvlGF4fqzuuOMOY9y4cV6syHfV5ljdeeedxhNPPGGkpqYa11xzTcMU5iOKi4uNTp06GRs3bvTo51BD/sxiZKaWioqKJEmtW7c2uRLfV5djVV5eruLi4iZ1fMvKyrRmzRqdOnVKSUlJZpfj0+pyrJrid0qq27H67LPPlJWVpeTkZC9X51tqe6yWL1+uAwcOKDU1tQGr8x3Tpk3T0KFDdfPNN5tdSo1MfTaTVRiGoUcffVQ33HCDEhISzC7Hp9X1WC1ZskSnTp3SmDFjvFidb9i1a5eSkpJ0+vRphYaGau3atYqPjze7LJ90KceqKX2npLodq/bt2+s///mPzp49q7S0ND3wwAMNVK25PDlW+/fv1+OPP67NmzcrIKDp/cpcs2aNdu7cqW3btpldygU1vb+ZOnj44Yf1xRdf6JNPPjG7FJ9Xl2P1+uuvKy0tTe+++64iIiK8WJ1v6NKliz7//HOdOHFCf/vb3zRhwgRt2rSJQFONuh6rpvadkup2rDZv3qyTJ09q69atevzxx3XllVdq7NixDVi1OWp7rMrKynT33Xdr7ty56ty5s0nVmicvL0+PPPKINmzYoKCgILPLubAGOZllYQ8//LDRvn174+uvv/ZovaZ4fUNdjtWaNWuM4OBg4+9//7sXK/NtAwcONP7nf/7nov2a4nfqfLU5VnynKtT2e1Xp6aefNjp37uzFinxXTcfq+PHjhiTD39/f9bLZbK62Dz74wIRqG87atWurfH5Jhs1mM/z9/Y2zZ89ecP2G/JnFyEwNDMPQ9OnTtXbtWmVmZiouLs7sknxWXY/V66+/rkmTJun111/X0KFDvVyl7zIMQyUlJWaXYQkXO1Z8p/7L0+9VU/4e1vTZw8PDtWvXLre2l156SR9++KH++te/NvrfCwMHDqzy+e+77z517dpVs2bNkr+/v0mVVUWYqcG0adP02muv6d1331VYWJgKCgokSXa7XcHBwZKk2bNn67vvvtOqVatc633++eeSpJMnT+o///mPPv/8czVv3rxRn0Koy7F6/fXXNX78eD333HPq06ePa53g4GDZ7XZzPkgDmDNnjoYMGaKYmBgVFxdrzZo1yszMVEZGhiS+U+fy9Fg11e+U5PmxevHFF3X55Zera9eukirmnfnd736n6dOnm/YZGoonx8rPz6/KtX8REREKCgpqEtdPhoWFVfmcLVq0UJs2bVztPvMzq0HGfyxIUrWv5cuXu/pMmDDBSE5Ovuh6sbGxDVp7Q6vLsUpOTq52nQkTJjR4/Q1p0qRJRmxsrNG8eXOjXbt2xsCBA40NGza4lvOd+i9Pj1VT/U4ZhufH6vnnnze6detmhISEGOHh4UZiYqLx0ksvGWVlZSZU37Dq8m/wXE3x1uxznX/qyFd+Ztl+2jEAAIAlMc8MAACwNMIMAACwNMIMAACwNMIMAACwNMIMAACwNMIMAACwNMIMAACwNMIMAACwNMIMAMtatmyZYmJi5Ofnp6VLlyotLU3XXnut2WVd0MGDB2Wz2VxTvgO4dMwADFhcYWGhnnzySa1fv15HjhxRq1atdM011ygtLU1JSUmSJJvNprVr1+r22283t9h65HQ61bZtWz3zzDMaNWqU7Ha7ysvLVVJSojZt2phdniRp4sSJOnHihN555x1XW1lZmf7zn/+obdu2Cgjg8XhAfeBfEmBxo0aNUmlpqVauXKmOHTvqyJEj+uCDD3Ts2LF631dpaamaNWtW79uti2+++UalpaUaOnSoHA6Hqz00NNTr+76U4+Dv76+oqKh6rgho4rz65CcAXnX8+HFDkpGZmVljn9jY2Bof+PbSSy8ZHTt2NJo1a2Z07tzZWLVqldu6koz09HRj+PDhRkhIiPHUU0+5HrT3yiuvGDExMUaLFi2MKVOmGGfPnjUWLVpkREZGGu3atTPmzZt30fpfeeUVIz4+3mjevLkRFRVlTJs2zbXs0KFDxvDhw40WLVoYYWFhxujRo42CggLDMAxj+fLlVR5kl5ubW+UhgKWlpcb06dMNu91utG7d2pg5c6Yxfvx4Y8SIEW7H59lnn3Wr65prrjFSU1MveBzOnj1rTJo0yejQoYMRFBRkdO7c2Vi6dKlrndTU1Co1fvTRR0Zubq4hyfjss89cfTMzM43rr7/edRxmzZpllJaWupYnJycb06dPN371q18ZrVq1MiIjI93qA5o6wgxgYaWlpUZoaKgxY8YM4/Tp09X2KSwsdD3FPD8/3ygsLDQMwzDefvtto1mzZsaLL75o7Nu3z1iyZInh7+9vfPjhh651JRkRERHGK6+8Yhw4cMA4ePCgkZqaaoSGhho///nPjT179hjr1q0zmjdvbgwaNMiYPn268eWXXxp//vOfDUnGli1baqz9pZdeMoKCgoylS5ca+/btM7Kzs12hory83EhMTDRuuOEGY/v27cbWrVuNHj16uJ7O+8MPPxj/+te/DElGdna2kZ+fb5w9e7ZKmJk3b57RunVr4+233zb27t1rTJkyxQgPD69TmDn/OJw5c8Z46qmnjOzsbOPrr782Vq9ebYSEhBhvvPGGYRiGUVxcbIwZM8YYPHiwkZ+fb+Tn5xslJSVVwsy3335rhISEGFOnTjX27t1rrF271mjbtq3b/pOTk43w8HAjLS3N+Oqrr4yVK1caNpvN7WnPQFNGmAEs7q9//avRqlUrIygoyOjbt68xe/Zs49///rdbH0nG2rVr3dr69u1rTJ482a1t9OjRxq233uq23owZM9z6pKamGiEhIYbT6XS1DRo0yOjQoYNRVlbmauvSpYuxYMGCGuuOjo42fv3rX1e7bMOGDYa/v7/xzTffuNr27NnjCi+GYRifffaZa0Tm3NrODTORkZHGb3/7W9f7s2fPGpdffnmdwsz5x6E6U6dONUaNGuV6P2HCBLd9GYZRJczMmTPH6NKli1FeXu7q8+KLLxqhoaGu45mcnGzccMMNbtu5/vrrjVmzZl20JqAp4G4mwOJGjRqlw4cPa926dRo0aJAyMzPVo0cPrVix4oLr7d27V/369XNr69evn/bu3evW1rNnzyrrdujQQWFhYa73kZGRio+Pl5+fn1tbYWFhtfsuLCzU4cOHNXDgwBpri4mJUUxMjKstPj5eLVu2rFJfTYqKinTkyBH16tXL1ebv76/rrruuVuufr7rj8Ic//EE9e/ZUu3btFBoaqpdfflnffPONR9vdu3evkpKSZLPZXG39+vXTyZMn9e2337raunfv7raew+Go8fgCTQ1hBmgEgoKClJKSoqeeekpZWVmaOHGiUlNTL7reub9AJckwjCptLVq0qLLe+Re/2my2atvKy8ur3W9wcPAF66qujgu1X0h1n/Fcfn5+VdpKS0urbOf84/Dmm2/ql7/8pSZNmqQNGzbo888/13333aczZ854VF91n6mynnPbPTm+QFNDmAEaofj4eJ06dcr1vlmzZiorK3Prc9VVV+mTTz5xa8vKytJVV13l9frCwsLUoUMHffDBB9Uuj4+P1zfffKO8vDxXW05OjoqKimpdn91uV2RkpLKzs11tZWVl+uyzz9z6tWvXTvn5+a73TqdTubm5F93+5s2b1bdvX02dOlWJiYm68sordeDAAbc+zZs3r3LczxcfH6+srCy3QJWVlaWwsDBddtllF60DAGEGsLSjR49qwIABWr16tb744gvl5ubqrbfe0uLFizVixAhXv8rgUFBQoOPHj0uSfvWrX2nFihX6wx/+oP379+uZZ57R22+/rccee6xBak9LS9OSJUv0/PPPa//+/dq5c6d+//vfS5Juvvlmde/eXffcc4927typ7OxsjR8/XsnJydWe7qnJ9OnTtWDBAr377rvat2+fHnnkER0/ftxtxGPAgAH6y1/+os2bN2v37t2aMGGC/P39L7rtK6+8Utu3b9f777+vr776Sk8++aS2bdvm1qdDhw764osvtG/fPn3//ffVjvhMnTpVeXl5mj59ur788ku9++67Sk1N1aOPPup22g5AzfiXAlhYaGioevfurWeffVY/+9nPlJCQoCeffFKTJ0/WCy+84Oq3ZMkSbdy4UTExMUpMTJQk3X777Xruuef029/+Vt26ddMf//hHLV++XP3792+Q2idMmKClS5fqpZdeUrdu3XTbbbdp//79kipOobzzzjtq1aqVfvazn+nmm29Wx44d9cYbb3i0j1mzZmns2LEaP368kpKSFBoaqkGDBikoKMjVZ/bs2frZz36m2267Tbfeeqtuv/12XXHFFRfd9pQpUzRy5Ejdeeed6t27t44ePaqpU6e69Zk8ebK6dOniuq7m008/rbKdyy67TP/85z+VnZ2ta665RlOmTNH999+vJ554wqPPCjRlzAAMoMkoLy/XVVddpTFjxujpp582uxwA9YQZgAE0WocOHdKGDRuUnJyskpISvfDCC8rNzdXdd99tdmkA6hGnmQA0Wn5+flqxYoWuv/569evXT7t27dK//vWvBrnIGUDD4TQTAACwNEZmAACApRFmAACApRFmAACApRFmAACApRFmAACApRFmAACApRFmAACApRFmAACApf1/DJRxiwrSbMgAAAAASUVORK5CYII=", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "sp_agg = pmp_spring_agg.sel(window=1, y=46.04, x=284.06)\n", "sm_agg = pmp_summer_agg.sel(window=1, y=46.04, x=284.06)\n", @@ -433,7 +361,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.12.3" + "version": "3.12.7" } }, "nbformat": 4, From 55514144463de754773d4081b69bfff7c639ad27 Mon Sep 17 00:00:00 2001 From: RondeauG Date: Wed, 20 Nov 2024 11:56:10 -0500 Subject: [PATCH 2/6] upd changelog --- CHANGELOG.rst | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/CHANGELOG.rst b/CHANGELOG.rst index 8176fa7f..70db10c8 100644 --- a/CHANGELOG.rst +++ b/CHANGELOG.rst @@ -8,8 +8,9 @@ Contributors to this version: Thomas-Charles Fortier Filion (:user:`TC-FF`). Internal changes ^^^^^^^^^^^^^^^^ -land_use_classification default collection has been changed to io-lulc-annual-v02 as previous one will be deprecated december 2024. (:pull:`227`). -Also added some collection, year, resolution and history attributes to xarray output of land_use_classification. (:pull:`227`). +* land_use_classification default collection has been changed to io-lulc-annual-v02 as previous one will be deprecated december 2024. (:pull:`227`). +* Also added some collection, year, resolution and history attributes to xarray output of land_use_classification. (:pull:`227`). +* Added a downloader agent to fix an issue related to ``pooch`` in recent ReadTheDocs builds. (:pull:`231`). v0.4.1 (2024-11-07) ------------------- From 4405efdbb2bdee669ec1f08e0217b55d3d035ea9 Mon Sep 17 00:00:00 2001 From: RondeauG Date: Wed, 20 Nov 2024 12:51:00 -0500 Subject: [PATCH 3/6] add missing import --- docs/notebooks/hydrological_modelling.ipynb | 1 + docs/notebooks/pmp.ipynb | 1 + 2 files changed, 2 insertions(+) diff --git a/docs/notebooks/hydrological_modelling.ipynb b/docs/notebooks/hydrological_modelling.ipynb index bfb10c4d..06c040df 100644 --- a/docs/notebooks/hydrological_modelling.ipynb +++ b/docs/notebooks/hydrological_modelling.ipynb @@ -106,6 +106,7 @@ "metadata": {}, "outputs": [], "source": [ + "import pooch\n", "import xarray as xr\n", "\n", "from xhydro.testing.helpers import deveraux\n", diff --git a/docs/notebooks/pmp.ipynb b/docs/notebooks/pmp.ipynb index 5ecde174..e825758f 100644 --- a/docs/notebooks/pmp.ipynb +++ b/docs/notebooks/pmp.ipynb @@ -25,6 +25,7 @@ "\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", + "import pooch\n", "import xarray as xr\n", "\n", "from xhydro.indicators import pmp\n", From 23ce5b016f3d5c450626738a6538859ccb9324a6 Mon Sep 17 00:00:00 2001 From: RondeauG Date: Wed, 20 Nov 2024 13:11:03 -0500 Subject: [PATCH 4/6] fix missing import --- docs/notebooks/pmp.ipynb | 2 ++ 1 file changed, 2 insertions(+) diff --git a/docs/notebooks/pmp.ipynb b/docs/notebooks/pmp.ipynb index e825758f..88057563 100644 --- a/docs/notebooks/pmp.ipynb +++ b/docs/notebooks/pmp.ipynb @@ -47,6 +47,8 @@ "metadata": {}, "outputs": [], "source": [ + "import xhydro as xh\n", + "\n", "downloader = pooch.HTTPDownloader(headers={\"User-Agent\": f\"xHydro-{xh.__version__}\"})\n", "path_day_zip = deveraux().fetch(\n", " \"pmp/CMIP.CCCma.CanESM5.historical.r1i1p1f1.day.gn.zarr.zip\",\n", From 6836d340b97a8e92202a4741238c417a44ba935d Mon Sep 17 00:00:00 2001 From: RondeauG Date: Wed, 20 Nov 2024 13:36:34 -0500 Subject: [PATCH 5/6] missed one usage --- docs/notebooks/optimal_interpolation.ipynb | 20 ++++++++++---------- 1 file changed, 10 insertions(+), 10 deletions(-) diff --git a/docs/notebooks/optimal_interpolation.ipynb b/docs/notebooks/optimal_interpolation.ipynb index 10dc3667..14480a22 100644 --- a/docs/notebooks/optimal_interpolation.ipynb +++ b/docs/notebooks/optimal_interpolation.ipynb @@ -36,8 +36,10 @@ "import xarray as xr\n", "from scipy.stats import norm\n", "\n", + "import xhydro as xh\n", "import xhydro.optimal_interpolation.optimal_interpolation_fun as opt\n", - "from xhydro.optimal_interpolation.ECF_climate_correction import general_ecf" + "from xhydro.optimal_interpolation.ECF_climate_correction import general_ecf\n", + "from xhydro.testing.helpers import deveraux" ] }, { @@ -324,16 +326,14 @@ "metadata": {}, "outputs": [], "source": [ - "# Get data with pooch\n", - "test_data_path = pooch.retrieve(\n", - " url=\"https://github.com/hydrologie/xhydro-testdata/raw/main/data/optimal_interpolation/OI_data_corrected.zip\",\n", - " known_hash=\"md5:acdf90b78b53595eb97ff0e84fc07aa8\",\n", + "# Get data\n", + "downloader = pooch.HTTPDownloader(headers={\"User-Agent\": f\"xHydro-{xh.__version__}\"})\n", + "test_data_path = deveraux().fetch(\n", + " \"optimal_interpolation/OI_data_corrected.zip\",\n", + " pooch.Unzip(),\n", + " downloader=downloader,\n", ")\n", - "\n", - "# Extract to a cache path.\n", "directory_to_extract_to = Path(test_data_path).parent\n", - "with ZipFile(test_data_path, \"r\") as zip_ref:\n", - " zip_ref.extractall(directory_to_extract_to)\n", "\n", "# Read-in all the files and set to paths that we can access later.\n", "flow_obs_info_file = directory_to_extract_to / \"A20_HYDOBS_TEST_corrected.nc\"\n", @@ -678,7 +678,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.12.3" + "version": "3.12.7" }, "vscode": { "interpreter": { From bd91a2bce6f9fda6996545bb93ff566d3120a8cc Mon Sep 17 00:00:00 2001 From: RondeauG Date: Wed, 20 Nov 2024 15:38:44 -0500 Subject: [PATCH 6/6] fix call for zipped downloads --- docs/notebooks/optimal_interpolation.ipynb | 2 +- docs/notebooks/pmp.ipynb | 6 ++++-- 2 files changed, 5 insertions(+), 3 deletions(-) diff --git a/docs/notebooks/optimal_interpolation.ipynb b/docs/notebooks/optimal_interpolation.ipynb index 14480a22..6dc66e80 100644 --- a/docs/notebooks/optimal_interpolation.ipynb +++ b/docs/notebooks/optimal_interpolation.ipynb @@ -333,7 +333,7 @@ " pooch.Unzip(),\n", " downloader=downloader,\n", ")\n", - "directory_to_extract_to = Path(test_data_path).parent\n", + "directory_to_extract_to = Path(test_data_path[0]).parent\n", "\n", "# Read-in all the files and set to paths that we can access later.\n", "flow_obs_info_file = directory_to_extract_to / \"A20_HYDOBS_TEST_corrected.nc\"\n", diff --git a/docs/notebooks/pmp.ipynb b/docs/notebooks/pmp.ipynb index 88057563..fb9e6738 100644 --- a/docs/notebooks/pmp.ipynb +++ b/docs/notebooks/pmp.ipynb @@ -62,10 +62,12 @@ ")\n", "\n", "path_day_zarr = (\n", - " Path(path_day_zip).parent / \"CMIP.CCCma.CanESM5.historical.r1i1p1f1.day.gn.zarr\"\n", + " Path(path_day_zip[0]).parents[1]\n", + " / \"CMIP.CCCma.CanESM5.historical.r1i1p1f1.day.gn.zarr\"\n", ")\n", "path_fx_zarr = (\n", - " Path(path_fx_zip).parent / \"CMIP.CCCma.CanESM5.historical.r1i1p1f1.fx.gn.zarr\"\n", + " Path(path_fx_zip[0]).parents[1]\n", + " / \"CMIP.CCCma.CanESM5.historical.r1i1p1f1.fx.gn.zarr\"\n", ")" ] },