{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "2c72937c-8213-4523-a0b6-ae70c5b77eb7",
   "metadata": {
    "editable": true,
    "slideshow": {
     "slide_type": ""
    },
    "tags": []
   },
   "source": [
    "# Using `pymf6`\n",
    "\n",
    "We can achieve the same effect with `pymf6`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "49045458-1e11-4b6c-aac5-899ccf846a79",
   "metadata": {},
   "outputs": [],
   "source": [
    "from modflowapi import Callbacks as states\n",
    "from xmipy.errors import XMIError\n",
    "\n",
    "from pymf6.mf6 import MF6\n",
    "from pymf6.api import create_mutable_bc\n",
    "\n",
    "from helpers import plot"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4ef1ed48-177a-4324-b71f-51dc765c7d44",
   "metadata": {},
   "source": [
    "We set up some initial data:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "e25c7887-3a7a-418e-8e90-522fff495087",
   "metadata": {},
   "outputs": [],
   "source": [
    "h_mean = 320.0\n",
    "name = 'rivercond'\n",
    "flux = []\n",
    "chd = []"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "01406375-ef8b-4f4c-b75b-8c34f6d4d25e",
   "metadata": {},
   "source": [
    "Now, we create an instance of `MF6`, get a reference to the flow model and the model loop:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "1949d16a-37d0-4769-9472-ed60e1228a29",
   "metadata": {},
   "outputs": [],
   "source": [
    "mf6 = MF6(name)\n",
    "gwf = mf6.models['gwf6']['rivercond']\n",
    "loop = mf6.model_loop()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3bd90dc0-c558-4c46-ab53-ad9b48afa1e0",
   "metadata": {},
   "source": [
    "We create a mutable river boundary condition and get the initial value of the river conductance:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "297c1e2a-b1fa-44cc-b1f2-8e6812ff06dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "riv = create_mutable_bc(gwf.riv)\n",
    "condref = riv.cond[0]"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2ec6ea40-907c-46fa-b205-57ae9ad74b16",
   "metadata": {},
   "source": [
    "Now, we can loop over all time steps, modified the conductance and collect the flux values in both cells:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "09171870-00cd-469e-8928-d9e9ac7dc055",
   "metadata": {},
   "outputs": [],
   "source": [
    "for model in loop:\n",
    "    if model.state == states.iteration_start:\n",
    "        if gwf.X[0, 0, 0] > h_mean:\n",
    "            riv.cond = condref\n",
    "        else:\n",
    "            riv.cond = condref * 0.10\n",
    "    elif model.state == states.timestep_end:\n",
    "        flux.append(mf6.sim_values.river[0])\n",
    "        chd.append(mf6.sim_values.chd_0[0])"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8a160da3-5f60-4282-8126-ddd4dfb3d800",
   "metadata": {},
   "source": [
    "The results are the same:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "4cad53e2-b6c4-4ad3-81e8-ec1df1fb40d7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc0AAAH/CAYAAAA194a7AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAXblJREFUeJzt3XtcjYfjB/DP6R7psIpcojCJXKZcYu6USsnYyqUm1tbMsJgxzOXLwr6rzebOmBjNJSulacj9mnKZXEaUKREVUqme3x++zm9nFY9O9Zxz+rxfr/N66bl+zoPz6XnOc5EJgiCAiIiIXklH6gBERESagqVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSTXehg0bIJPJynxNnTpVMZ21tTXGjBlT7fnmzp0LmUxWLeu6efMmZDIZNmzYUC3rk0KfPn3Qp08fpWEymQxz586VJA9pFj2pAxCpi/Xr16N169ZKwxo1aiRRmv/3wQcfYNCgQdWyroYNG+L48eNo0aJFtayPSNOwNIn+x97eHo6OjlLHKKVJkyZo0qRJla6juLgYRUVFMDQ0RLdu3ap0XWV5+vQpjIyMqm2PmqiieHiWqIJeHNa9efOm0vD4+HjIZDLEx8cDAK5duwZTU1O8++67StPt378furq6mD179kvXU9bhWWtrawwePBixsbHo1KkTjI2N0bp1a/z000+vzP3iEOySJUuwYMEC2NjYwNDQEAcOHCh1eHbXrl2QyWTYt29fqeWsWLECMpkM58+fVww7c+YMPD098cYbb8DIyAhvvfUWfv31V6X5Xmy3vXv3YuzYsbCwsECtWrVQUFBQbubs7GxMmTIFzZs3h6GhIerXrw83NzdcvnxZMU1hYSEWLFiA1q1bw9DQEBYWFvD398e9e/deuU2IxGJpEv3Pi72tf74qw5tvvok1a9Zg+/btWLp0KQAgIyMDI0eORM+ePSv8Xdq5c+cwZcoUfPbZZ/jtt9/Qvn17jBs3DocOHRI1/9KlS7F//37897//xZ49e0odmgaAwYMHo379+li/fn2pcRs2bECnTp3Qvn17AMCBAwfQo0cPZGdnY+XKlfjtt9/QsWNHeHt7l/kd6dixY6Gvr4+wsDBs374d+vr6ZeZ89OgR3n77baxatQr+/v6IiorCypUr0apVK6SnpwMASkpKMGTIECxatAgjR45EdHQ0Fi1ahLi4OPTp0wdPnz4VtU2IXoWHZ4n+p6zDks+ePYOenur/Tby9vXHw4EF8/vnn6NKlC2bOnAlBELBlyxbo6upWaJn379/H0aNH0bRpUwBAr169sG/fPvzyyy/o1avXK+c3MjLC77//rlRW/95r1tPTw+jRo7FixQrk5ORALpcDAJKTk3Hq1Cn88MMPimnHjx+Ptm3bYv/+/Ypt5uLigvv37+PLL7+En58fdHT+//f0/v37Y9WqVa/M+d133+HPP/9EXFwcBgwYoBj+zjvvKP7866+/IjY2Fjt27FAa3qFDB3Tu3BkbNmzAxx9//Mp1Eb0K9zSJ/mfjxo04ffq00qsyCvOF0NBQtG3bFn379kV8fDw2bdqEhg0bVnh5HTt2VBQm8LwEW7VqhVu3boma39PTs9y9u38aO3Ysnj59ivDwcMWw9evXw9DQECNHjgQA/PXXX7h8+TJGjRoFAEp7625ubkhPT8eVK1eUljts2DBROffs2YNWrVopFea/7d69G3Xr1oWHh4fSujt27AhLS0vFoXIiVbE0if7Hzs4Ojo6OSq/K9KJk8vPz0bFjRwwcOFCl5ZmZmZW5DrGHIsUWdtu2bdG5c2fFIdri4mJs2rQJQ4YMwRtvvAEAuHv3LgBg6tSp0NfXV3qNHz8ewPM944qs/969e688Eeru3bvIzs6GgYFBqfVnZGSUWjdRRfHwLFEFGRkZAUCpE1jK+4C+ePEivvrqK3Tu3BmnT59GSEgIgoKCqjxneV7nTFV/f3+MHz8eycnJuHHjBtLT0+Hv768Yb25uDgCYMWOG0uHRf7K1ta3Q+i0sLHD79u2XTmNubg4zMzPExsaWOb5OnTqi1kX0KixNogqytrYGAJw/f16pECIjI0tN++TJE7z77ruwtrbGgQMHMH36dEyfPh09evRA165dqytyhY0YMQJBQUHYsGEDbty4gcaNG8PZ2Vkx3tbWFm+++SbOnTuHr7/+ulLX7erqiq+++gr79+9Hv379ypxm8ODB2Lp1K4qLizVie5LmYmkSVVDnzp1ha2uLqVOnoqioCPXq1UNERASOHDlSatrAwECkpqbi1KlTqF27Nr799lscP34cPj4+SExMRN26dav/DbyGunXrYujQodiwYQOys7MxdepUpZN6AGDVqlVwdXWFi4sLxowZg8aNG+PBgwdITk7G2bNnsW3btgqte/LkyQgPD8eQIUMwffp0dOnSBU+fPsXBgwcxePBg9O3bFz4+Pti8eTPc3NwwadIkdOnSBfr6+rh9+zYOHDiAIUOGYOjQoZWxKaiG43eaRBWkq6uLqKgotG7dGoGBgfDz84OhoSF+/PFHpenWrl2LTZs2YdmyZWjbti0AwMDAAOHh4Xjw4IHSYU515u/vj8zMTBQWFpZ5O8G+ffvi1KlTqFu3LiZPnowBAwbg448/xh9//PHSk3hepU6dOjhy5AjGjRuH1atXw93dHQEBAbhy5Yrijk26urqIjIzEl19+iZ07d2Lo0KHw8vLCokWLYGRkhHbt2lV4/UT/JBMEQZA6BBERkSbgniYREZFILE0iIiKRWJpEREQiaVVpLl++HDY2NjAyMoKDgwMOHz780ukLCgowc+ZMNGvWDIaGhmjRooWoG14TEVHNpDWXnISHh2Py5MlYvnw5evTooTj9/dKlS0q3Gvun9957D3fv3sW6devQsmVLZGZmVtpNuomISPtozdmzXbt2RadOnbBixQrFMDs7O3h5eSE4OLjU9LGxsfDx8cGNGzcUtwIjIiJ6Ga3Y0ywsLERCQgKmT5+uNNzZ2RnHjh0rc57IyEg4OjpiyZIlCAsLQ+3ateHp6Yn//Oc/MDY2LnOegoICpVumlZSU4MGDBzAzM+PDc4mI1IAgCHj06BEaNWpU6gYclUErSvP+/fsoLi5GgwYNlIY3aNAAGRkZZc5z48YNHDlyBEZGRoiIiMD9+/cxfvx4PHjwoNzvNYODgzFv3rxKz09ERJUrLS3tlTf6rwitKM0X/r23JwhCuXuAJSUlkMlk2Lx5s+IZgSEhIRg+fDiWLVtW5t7mjBkzlG6wnZOTg6ZNmyItLQ2mpqaV+E6IiKgicnNzYWVlVWU36deK0jQ3N4eurm6pvcrMzMxSe58vNGzYEI0bN1YUJvD8O1BBEHD79m28+eabpeYxNDSEoaFhqeGmpqYsTSIiNVJVX5lpxSUnBgYGcHBwQFxcnNLwuLg4dO/evcx5evTogTt37uDx48eKYVevXoWOjk6V7NITEZHm04rSBICgoCCsXbsWP/30E5KTk/HZZ58hNTUVgYGBAJ4fWvXz81NMP3LkSJiZmcHf3x+XLl3CoUOH8Pnnn2Ps2LHlnghEREQ1m1YcngUAb29vZGVlYf78+UhPT4e9vT1iYmLQrFkzAEB6ejpSU1MV05uYmCAuLg6ffvopHB0dYWZmhvfeew8LFiyQ6i0QEZGa05rrNKWQm5sLuVyOnJwcfqdJpAEEQUBRURGKi4uljkIq0NfXh66ubpnjqvpzWWv2NImIXqawsBDp6enIy8uTOgqpSCaToUmTJjAxMan2dbM0iUjrlZSUICUlBbq6umjUqBEMDAx4QxINJQgC7t27p7jKobw9zqrC0iQirVdYWIiSkhJYWVmhVq1aUschFVlYWODmzZt49uxZtZem1pw9S0T0KlVxWzWqflIeJeC/ICIiIpFYmkRERCKxNImINJxMJsOuXbskWffRo0fRrl076Ovrw8vLC/Hx8ZDJZMjOzpYkT1VjaRIRqakxY8ZAJpNBJpNBT08PTZs2xccff4yHDx8qTZeeng5XV1dJMgYFBaFjx45ISUnBhg0bJMlQnViaRERqbNCgQUhPT8fNmzexdu1aREVFYfz48UrTWFpalvkwicpSXFyMkpKSMsddv34d/fr1Q5MmTVC3bt0qy6AuWJpERGrM0NAQlpaWaNKkCZydneHt7Y29e/cqTfPPw7NOTk6YPn260vh79+5BX18fBw4cAPD8Epxp06ahcePGqF27Nrp27Yr4+HjF9Bs2bEDdunWxe/dutGnTBoaGhrh165bSMm/evAmZTIasrCyMHTsWMpmszD3NuXPnomPHjkrDvvvuO1hbWwMA8vPz0bZtW3z44YeK8SkpKZDL5VizZs1rbKnqwes0iajGysvLw+XLl6t1na1bt67wtaI3btxAbGws9PX1y51m1KhR+OabbxAcHKy4NCM8PBwNGjRA7969AQD+/v64efMmtm7dikaNGiEiIgKDBg3ChQsXFI9FzMvLQ3BwMNauXQszMzPUr19faT1WVlZIT0+Hra0t5s+fD29vb8jlcpw8efK13pORkRE2b96Mrl27ws3NDR4eHvD19UXfvn0REBDwWsuqDixNIqqxLl++DAcHh2pdZ0JCAjp16iR6+t27d8PExATFxcXIz88HAISEhJQ7vbe3Nz777DMcOXIEPXv2BAD88ssvGDlyJHR0dHD9+nVs2bIFt2/fRqNGjQAAU6dORWxsLNavX4+vv/4aAPDs2TMsX74cHTp0KHM9urq6sLS0hEwmg1wuh6Wlpej39G8dO3bEggULEBAQgBEjRuD69euSndj0KixNIqqxWrdujYSEhGpf5+vo27cvVqxYgby8PKxduxZXr17Fp59+Wu70FhYWGDhwIDZv3oyePXsiJSUFx48fx4oVKwAAZ8+ehSAIaNWqldJ8BQUFMDMzU/xsYGCA9u3bv1ZWVUyZMgW//fYbfvjhB+zZswfm5ubVtu7XwdIkohqrVq1ar7XXJ4XatWujZcuWAIClS5eib9++mDdvHv7zn/+UO8+oUaMwadIk/PDDD/jll1/Qtm1bxR5jSUkJdHV1kZCQUOoWdP+8AbqxsXGl3HlHR0cH/36Y1rNnz0pNl5mZiStXrkBXVxfXrl3DoEGDVF53VeCJQEREGmTOnDn473//izt37pQ7jZeXF/Lz8xEbG4tffvkFo0ePVox76623UFxcjMzMTLRs2VLppcoh1vJYWFggIyNDqTiTkpJKTTd27FjY29tj48aNmDZtGi5dulTpWSoDS5OISIP06dMHbdu2VXz3WJbatWtjyJAhmD17NpKTkzFy5EjFuFatWmHUqFHw8/PDzp07kZKSgtOnT2Px4sWIiYmpkrz37t3DkiVLcP36dSxbtgx79uxRmmbZsmU4fvw4Nm7ciJEjR2L48OEYNWoUCgsLKz2PqliaREQaJigoCGvWrEFaWlq504waNQrnzp1Dz5490bRpU6Vx69evh5+fH6ZMmQJbW1t4enri5MmTsLKyqvSsdnZ2WL58OZYtW4YOHTrg1KlTmDp1qmL85cuX8fnnn2P58uWK9S9btgzZ2dmYPXt2pedRlUz498FmEq2qnxBORJUjPz8fKSkpsLGxgZGRkdRxSEUv+/us6s9l7mkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmEdUYPO9RO0j598jSJCKt9+IG53l5eRInocrw4vrNf9/RqDrwNnpEpPV0dXVRt25dZGZmAnh++7zKuEUcVb+SkhLcu3cPtWrVgp5e9VcYS5OIaoQXt4h7UZykuXR0dNC0aVNJfvFhaRJRjSCTydCwYUPUr1+/zBuGk+YwMDCAjo403y6yNImoRtHV1ZXkuzDSDjwRiIiISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCSSVpXm8uXLYWNjAyMjIzg4OODw4cOi5jt69Cj09PTQsWPHqg1IREQaTWtKMzw8HJMnT8bMmTORmJiInj17wtXVFampqS+dLycnB35+fujfv381JSUiIk0lEwRBkDpEZejatSs6deqEFStWKIbZ2dnBy8sLwcHB5c7n4+ODN998E7q6uti1axeSkpJErzM3NxdyuRw5OTkwNTVVJT4REVWCqv5c1oo9zcLCQiQkJMDZ2VlpuLOzM44dO1bufOvXr8f169cxZ84cUespKChAbm6u0ouIiGoOrSjN+/fvo7i4GA0aNFAa3qBBA2RkZJQ5z7Vr1zB9+nRs3rwZenp6otYTHBwMuVyueFlZWamcnYiINIdWlOYLMplM6WdBEEoNA4Di4mKMHDkS8+bNQ6tWrUQvf8aMGcjJyVG80tLSVM5MRESaQ9wulpozNzeHrq5uqb3KzMzMUnufAPDo0SOcOXMGiYmJmDBhAgCgpKQEgiBAT08Pe/fuRb9+/UrNZ2hoCENDw6p5E0REpPa0Yk/TwMAADg4OiIuLUxoeFxeH7t27l5re1NQUFy5cQFJSkuIVGBgIW1tbJCUloWvXrtUVnYiINIhW7GkCQFBQEHx9feHo6AgnJyesXr0aqampCAwMBPD80Orff/+NjRs3QkdHB/b29krz169fH0ZGRqWGExERvaA1pent7Y2srCzMnz8f6enpsLe3R0xMDJo1awYASE9Pf+U1m0RERC+jNddpSoHXaRIRqRdep0lERKQmWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEkmrSnP58uWwsbGBkZERHBwccPjw4XKn3blzJwYOHAgLCwuYmprCyckJv//+ezWmJSIiTaM1pRkeHo7Jkydj5syZSExMRM+ePeHq6orU1NQypz906BAGDhyImJgYJCQkoG/fvvDw8EBiYmI1JyciIk0hEwRBkDpEZejatSs6deqEFStWKIbZ2dnBy8sLwcHBopbRtm1beHt746uvvipzfEFBAQoKChQ/5+bmwsrKCjk5OTA1NVXtDRARkcpyc3Mhl8ur7HNZK/Y0CwsLkZCQAGdnZ6Xhzs7OOHbsmKhllJSU4NGjR3jjjTfKnSY4OBhyuVzxsrKyUik3ERFpFq0ozfv376O4uBgNGjRQGt6gQQNkZGSIWsa3336LJ0+e4L333it3mhkzZiAnJ0fxSktLUyk3ERFpFj2pA1QmmUym9LMgCKWGlWXLli2YO3cufvvtN9SvX7/c6QwNDWFoaKhyTiIi0kxaUZrm5ubQ1dUttVeZmZlZau/z38LDwzFu3Dhs27YNAwYMqMqYRESk4bTi8KyBgQEcHBwQFxenNDwuLg7du3cvd74tW7ZgzJgx+OWXX+Du7l7VMYmISMNpxZ4mAAQFBcHX1xeOjo5wcnLC6tWrkZqaisDAQADPv4/8+++/sXHjRgDPC9PPzw/ff/89unXrpthLNTY2hlwul+x9EBGR+tKa0vT29kZWVhbmz5+P9PR02NvbIyYmBs2aNQMApKenK12zuWrVKhQVFeGTTz7BJ598ohj+/vvvY8OGDdUdn4iINIDWXKcphaq+HoiIiF4Pr9MkIiJSEyxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYlUoUeDBQUFiZ42JCSkIqsgIiJSOxUqzcTERKWfExISUFxcDFtbWwDA1atXoaurCwcHB9UTEhERqYkKleaBAwcUfw4JCUGdOnXw888/o169egCAhw8fwt/fHz179qyclERERGpA5YdQN27cGHv37kXbtm2Vhl+8eBHOzs64c+eOSgHVGR9CTUSkXtT+IdS5ubm4e/duqeGZmZl49OiRqosnIiJSGyqX5tChQ+Hv74/t27fj9u3buH37NrZv345x48bhnXfeqYyMREREaqFC32n+08qVKzF16lSMHj0az549e75QPT2MGzcO33zzjcoBiYiI1EWFv9P88ssv4eXlhS5dugAAnjx5guvXr0MQBLRs2RK1a9eu1KDqiN9pEhGpl6r+XK7wnmZ6ejoGDx4MXV1deHh4wMvLC/3794ehoWFl5iMiIlIbFf5Oc/369bh79y5+/fVX1K1bF0FBQTA3N8c777yDDRs24P79+5WZk4iISHIqX3LyT8nJyYiKisJvv/2GM2fOoGvXrvD09MSIESPQuHHjylqN2uDhWSIi9VLVn8uVWpr/dO/ePURGRiIyMhI9e/bE1KlTq2I1kmJpEhGpF40tzZqApUlEpF7U7uYGDx8+xIMHDwA835vcsWMHLl68WOnBiIiI1M1rnT27du1aBAcHo6SkBNOmTcPmzZvRvn17zJkzBxMnTsSHH35YVTmJSAOdOXMGoaGhKCgokDoKqSldXV2MHTsWLi4uUkcR5bUOz3bo0AEnT55EXl4emjZtipSUFFhYWCA3Nxe9evVCUlJSFUZVPzw8S1S+06dPo0ePHmjevDmsra2ljkNqKjMzE4mJidixY0el3EVOra7T1NXVhZGREYyMjNCyZUtYWFgAAExNTSGTySo9HBFppuLiYowfPx5t27bFyZMnYWBgIHUkUlMlJSXw8vLCxIkT4ezsDBMTE6kjvdRrfaepp6eH/Px8AMDBgwcVw3ljdiL6p7Vr1+LMmTNYtmwZC5NeSkdHB0uXLkVWVhbmz58vdZxXeq3S3L9/v+KOP3K5XDH86dOnWLduXeUmIyKNdO/ePcyYMQP+/v7o3r271HFIA1hbW2PmzJkIDQ3FpUuXpI7zUrzkRAX8TpOotA8++AA7duzA1atXFV/hEL1KQUEB7O3tYWVlhX379lX4Kz+1+k6zPPn5+Th//jwyMzNRUlKiNM7T07MyVkFEGuD48eNYt24dli9fzsKk12JoaIgffvgBrq6u2Lp1K0aMGCF1pDKpvKcZGxsLPz+/Mu81K5PJUFxcrMri1Rr3NIn+X3FxMRwdHaGrq4uTJ09CV1dX6kikgYYNG4bjx4/j8uXLFfpcVbubG/zbhAkT8O677yI9PR0lJSVKL20uTCJStnLlSpw7dw7Lly9nYVKFhYaGIicnB3PnzpU6SplULs3MzEwEBQWhQYMGlZGHiDTQ3bt3MXPmTHzwwQeKZ+wSVUTTpk0xe/ZsLF26FBcuXJA6Tikql+bw4cMRHx9fCVGISFN98cUX0NXVRXBwsNRRSAsEBQWhZcuW+OSTT6Bu56qq/J1mXl4e3n33XVhYWKBdu3bQ19dXGj9x4kSVAqozfqdJBBw5cgQ9e/bE6tWrERAQIHUc0hJ//PEHBg4ciI0bN8LX11f0fGr/lJO1a9ciMDAQxsbGMDMzUzpNWCaT4caNGyqHVFcsTarpioqK0KlTJxgbG+P48ePQ0VH54BWRgre3N+Lj43HlyhXUrVtX1DxqfyLQrFmzMH/+fOTk5ODmzZtISUlRvLS5MIkIWL58OS5evIjly5ezMKnShYSEIC8vD3PmzJE6ioLK/8oLCwvh7e3N/zBENcy9e/cwZ84cBAQEwMHBQeo4pIUaN26MWbNmYdmyZfjzzz+ljgOgEkrz/fffR3h4eGVkISINMnv2bAiCgAULFkgdhbTY5MmTYWNjg0mTJqnFSUEq3xGouLgYS5Yswe+//4727duXOhEoJCRE1VUQkZpJSkrC6tWrERISwjv/UJUyNDREaGgoPDw8sGvXLgwdOlTSPCqfCNS3b9/yFy6TYf/+/aosXq3xRCCqiQRBQJ8+fXDv3j2cO3eu1C/KRJVNEAS4u7sjOTkZly5dgrGxcbnTqv29Zw8cOFAZOYhIQ2zbtg2HDh3C77//zsKkaiGTyRAaGgp7e3t8++23mDVrlmRZePYOEYmWl5eHqVOnwtPTE87OzlLHoRrE1tYWkydPxtdff420tDTJcrA0iUi0JUuW4O7duzxXgSQxe/ZsmJqaYtq0aZJlYGkSkSi3bt3C4sWLERQUhBYtWkgdh2ogU1NTLFq0CFu3bsWhQ4ckycCHUKuAJwJRTfLee+/hyJEjuHLlCurUqSN1HKqhSkpK4OTkhIKCAiQkJJR6oo7a3xGIiLRffHw8tm3bhsWLF7MwSVI6OjpYunQpzp07hzVr1lT/+qt04To66NevHxISEqpyNURUhYqKijBp0iR069YNo0aNkjoOEbp27Yr3338fs2bNwoMHD6p13VVamj/99BN69+6t1U86IdJ2a9aswfnz5/H999/zdpmkNoKDg1FYWFjt96Xld5oq4HeapO0ePHiAVq1awcPDA+vXr5c6DpGSb775BjNmzEBiYiLatWsHgN9pEpGE5s6di8LCQj5cmtTSpEmT0KJFi2q9L22llObhw4cxevRoODk54e+//wYAhIWF4ciRI5WxeCKSwItHfs2ePRuWlpZSxyEqxcDAAKGhoThw4AB27NhRLetUuTR37NgBFxcXGBsbIzExEQUFBQCAR48e4euvv1Y5IBFVP0EQMHnyZDRv3hyTJk2SOg5Rudzc3ODm5obPP/8cT58+rfL1qVyaCxYswMqVK7FmzRql+1B2794dZ8+eVXXxRCSBXbt2Yd++fQgNDYWBgYHUcYheKiQkBLdv366WO1WpXJpXrlxBr169Sg03NTVFdna2qosnomqWn5+PKVOmwNXVFe7u7lLHIXolW1tbTJw4EcHBwbhz506Vrkvl0mzYsCH++uuvUsOPHDmC5s2bq7p4Iqpm3377LdLS0hAaGip1FCLRZs+ejVq1amHu3LlVuh6VS/Ojjz7CpEmTcPLkSchkMty5cwebN2/G1KlTMX78+MrIKNry5cthY2MDIyMjODg44PDhwy+d/uDBg3BwcICRkRGaN2+OlStXVlNSIvV0+/ZtfP3115g0aRJsbW2ljkMkWt26dbFw4UKEh4dX7YqESvDll18KxsbGgkwmE2QymWBkZCTMmjWrMhYt2tatWwV9fX1hzZo1wqVLl4RJkyYJtWvXFm7dulXm9Ddu3BBq1aolTJo0Sbh06ZKwZs0aQV9fX9i+fbvodebk5AgAhJycnMp6G0SSGjVqlFC/fn0hOztb6ihEr62oqEiwt7ev0s9llW9ukJqaiiZNmiA/Px+XLl1CSUkJ2rRpg9q1ayMtLQ1NmzatnHZ/ha5du6JTp05YsWKFYpidnR28vLzKvMbsiy++QGRkJJKTkxXDAgMDce7cORw/flzUOnlzA9Imx44dQ48ePbB27VqMGzdO6jhEFRITEwN3d3f1vbmBjY0N7t+/j1q1asHR0RFdunSBiYkJHjx4ABsbm8rI+EqFhYVISEgo9VBcZ2dnHDt2rMx5jh8/Xmp6FxcXnDlzBs+ePStznoKCAuTm5iq9gOcnQxFpuhkzZsDBwQH+/v5SRyGqsLfffrtKl6+n6gLK21F9/PgxjIyMVF28KPfv30dxcTEaNGigNLxBgwbIyMgoc56MjIwypy8qKsL9+/fRsGHDUvMEBwdj3rx5pYaPHDkSdnZ2AABjY2OEh4djwYIFGDt2LAIDA5WmnTlzJhITEzF8+HAsXrxYqXB79OiBIUOG4MKFC6hVqxZWrVqlNO+OHTsUzzP08fFRGjdhwgRkZWWhS5cu2LJlC06dOqUYZ29vjwkTJmDPnj2wtbXFkiVLlOb96aefsHLlSsyaNQvvvPMOioqKFON8fX1Rq1YtNGzYEEePHkVcXJxinJWVFb755husWLECzs7OmDlzptJyQ0NDsWXLFsyaNQv+/v7IyspSjBsyZAhsbW1RVFSE1NRUbN++XTHO1NQUmzZtwoIFC+Dn54cJEyYoLXfu3Lk4evQo/Pz8MGfOHNy4cUMxrnfv3nB2dsb169cV7+2fIiMjsWDBAnz66afw9fVVGhcUFITU1FT06dMH69atQ2JiomLcW2+9hXHjxiE+Ph5NmzYtdWp7WFgYfvjhB8yaNQuenp5K48aOHQsAaNGiBfbu3YuDBw8qxjVv3hzz5s3Dxo0b0aNHj1InMfz444/YuHEjZs2ahdGjRyt+UQOA4cOHo2nTptDT08OVK1fw22+/KcaZmZlh/fr1WLBgAUaMGIHPPvtMabkLFy7E3r178fHHH+Pzzz/HX3/9hUOHDqFjx45YtmwZevTogfT0dOTl5SEsLEwxn56eHnbu3IkFCxYgMDBQ8d5emDZtGq5cuQJXV1f8+OOPuHjxomJcly5dMGLECJw6dQpmZmb48ccflebdunUrQkJC8MUXX2DYsGFK4z766CPk5eWhXbt2+O2333D06FHFOFtbW3zxxRfYvn073nrrLSxcuFBp3pUrV+Knn37CrFmz4O3trXQt34gRI2BmZga5XI7ExETExMQoxjVs2BA//vgjQkJCMGTIkFIPPl6yZAl+++03BAUFYcKECUhPT1eMc3Nzw1tvvYWcnBxkZWVhy5YtinH8jPh/VfEZUd5OT2Wp8OHZoKAgAMD333+PgIAA1KpVSzGuuLgYJ0+ehK6urtI/7qpy584dNG7cGMeOHYOTk5Ni+MKFCxEWFobLly+XmqdVq1bw9/fHjBkzFMOOHj2Kt99+G+np6WXeAaWgoEBx8wbg+eFZKysrdOnSBSdPnqzkd0VUfTZt2gRfX99y/+0TaYqq/tqswnuaL34LFwQBFy5cULoA2sDAAB06dMDUqVNVTyiCubk5dHV1S+1VZmZmltqbfMHS0rLM6fX09GBmZlbmPIaGhjA0NCw1/PTp08jKyip3PiJ1FxMTAwcHBxYm0StUuDQPHDgAAPD398f3338v6YkwBgYGcHBwQFxcHIYOHaoYHhcXhyFDhpQ5j5OTE6KiopSG7d27F46Ojkp3NhJDEAT8/vvvGDly5OuHJ5JYcXExYmNjSx0GJ6LSVD4RaP369Wpx5mhQUBDWrl2Ln376CcnJyfjss8+Qmpqq+L5gxowZ8PPzU0wfGBiIW7duISgoCMnJyfjpp5+wbt26Cu0dt2/fXum7ECJNcvLkSTx8+BBubm5SRyFSeyqfCPTCpUuXkJqaisLCQqXh/z4poqp4e3sjKysL8+fPR3p6Ouzt7RETE4NmzZoBANLT05GamqqY3sbGBjExMfjss8+wbNkyNGrUCEuXLi11AoIYzs7OWL9+PYqLi6Grq1tp74moOsTExMDMzAydO3eWOgqR2lP5Os0bN25g6NChuHDhAmQymeJsWplMBuD5oR9t9eIL57i4OAwcOLDUiUhEmqBTp05o06YNNm3aJHUUIpWp/UOoJ02aBBsbG9y9exe1atXCn3/+iUOHDsHR0RHx8fGVEFH9OTg44I033uAhWtI4d+7cQWJiIlxdXaWOQqQRVC7N48ePY/78+bCwsICOjg50dHTw9ttvIzg4GBMnTqyMjGpPV1cXgwYNQnR0tNRRiF5LbGwsZDIZXFxcpI5CpBFULs3i4mKYmJgAeH7px4vHsjRr1qxG3SnH3d0diYmJVf5YGqLKtGfPHnTt2hXm5uZSRyHSCCqXpr29Pc6fPw/g+f1flyxZgqNHj2L+/Pk16tFgLi4ukMlk2LNnj9RRiER59uwZ9u7dy7NmiV6DyqU5a9YslJSUAAAWLFiAW7duoWfPnoiJicHSpUtVDqgpzMzM0K1bN36vSRrj2LFjyM3NZWkSvQaVLzmxs7NDkyZNADy/j+alS5fw4MED1KtXD2lpaSoH1CTu7u5YvHgxCgsLle6QRKSOYmJi0KBBA7z11ltSRyHSGJX2lJN/euONN6r1KSfqws3NDY8ePcKRI0ekjkL0SjExMXB1dYWOjsofA0Q1hsr/W9ThKSfqomPHjmjYsCEP0ZLaS01NxcWLF3lolug1Vfjw7IunnMhkMnz11VdlPuWkY8eOKgfUJDKZDG5uboiJicF///tfqeMQlWvPnj3Q1dXFwIEDpY5CpFG04ikn6sTNzQ3r1q1DSkpKjTs8TZojJiYGPXr0QN26daWOQqRRtOIpJ+pkwIAB0NfXR0xMDD755BOp4xCVUlBQgD/++ANfffWV1FGINI7WPOVEXZiamiouuSFSR4cOHUJeXh6/zySqAJ42VwXc3Nywf/9+5OXlSR2FqJSYmBg0adIE9vb2Ukch0jgszSrg5uaG/Pz8GnPDetIsMTExcHNzUzyJiIjEY2lWgdatWyue10mkTv766y9cvXqVh2aJKoilWQVeXHoSHR1d7nWsRFLYs2cP9PX10b9/f6mjEGkklUtz1KhRWL16Na5evVoZebSGu7s7bt68icuXL0sdhUhhz5496N27t+LJRET0elQuTRMTE4SEhKB169Zo1KgRRowYgZUrV9b4sujTpw+MjIx4iJbURl5eHg4cOMAHThOpQOXSXLVqFS5fvow7d+4gJCQEcrkc33//Pdq2bYuGDRtWRkaNZGxsjH79+vHB1KQ29u/fj/z8fLi7u0sdhUhjVdp3mnXq1EG9evVQr1491K1bF3p6erC0tKysxWskd3d3HD58GLm5uVJHIUJ0dDRatmyJVq1aSR2FSGOpXJpffPEFunXrBnNzc8yaNQuFhYWYMWMG7t69q7jVXk3l6uqKoqIixMXFSR2FajhBELB79264u7vzUhMiFaj8PM1vvvkGFhYWmDNnDoYMGQI7O7vKyKUVbGxsYGdnh5iYGAwbNkzqOFSDXbhwAbdv3+ahWSIVqbynmZiYiJkzZ+LUqVPo1asXLC0t4e3tjRUrViA5ObkyMmo0d3d3xMTE8NITklR0dDRq166NXr16SR2FSKPJhEr+ND937hy+++47bNq0CSUlJSguLq7MxauV3NxcyOVy5OTklHv/3QMHDqBfv35ISEhAp06dqjkh0XM9evRA/fr1ERERIXUUoiol5nNZFSofngWe723Gx8cjPj5eceJLx44d0bdv38pYvEbr0aMH6tSpg5iYGJYmSSIrKwsnTpzA6tWrpY5CpPFULs169erh8ePH6NChA/r06YOAgAD06tWLTz75HwMDAzg7OyMmJgazZs2SOg7VQLGxsSgpKeGt84gqgcqlGRYWxpJ8BTc3N3zwwQe4f/8+zM3NpY5DNUx0dDQ6depUo6+bJqosKp8INHjwYBbmK7i6ukIQBPz+++9SR6EapqioCHv27OFZs0SVpFK+08zOzsa6deuQnJwMmUwGOzs7jBs3DnK5vDIWr/EaNmyITp06ISYmBqNGjZI6DtUgx48fR3Z2NgYPHix1FCKtoPKe5pkzZ9CiRQuEhobiwYMHuH//PkJDQ9GiRQucPXu2MjJqBTc3N8TGxmr12cSkfqKjo1G/fn04OjpKHYVIK6hcmp999hk8PT1x8+ZN7Ny5ExEREUhJScHgwYMxefLkSoioHdzc3PDgwQOcPHlS6ihUg0RHR8PV1RU6OnwKIFFlqJQ9zS+++AJ6ev9/pFdPTw/Tpk3DmTNnVF281ujSpQvMzMz41BOqNrdu3cLFixf5fSZRJVK5NE1NTZGamlpqeFpaGurUqaPq4rWGrq4uBg0axKeeULWJjo6Gnp4enJ2dpY5CpDVULk1vb2+MGzcO4eHhSEtLw+3bt7F161Z88MEHGDFiRGVk1Bpubm5ISkrC33//LXUUqgGio6PRs2dPnpBHVIlUPnv2v//9L2QyGfz8/FBUVAQA0NfXx8cff4xFixapHFCbuLi4QEdHB7GxsRg3bpzUcUiL5eXlYf/+/ViwYIHUUYi0SqXdezYvLw/Xr1+HIAho2bIlatWqVRmLVWsVucdhjx490KBBA+zcubOK01FNtnv3bnh4eCA5ORmtW7eWOg5Rtanqe89W2il1tWrVQrt27dC+ffsaUZgV5e7ujr1796KgoEDqKKTFoqOj0aJFC9ja2kodhUirVOjwbFBQkOhpQ0JCKrIKreXh4YGZM2ciPj4eLi4uUschLSQIAqKjozF06FA+cJqoklWoNBMTE0VNx/+wpdnb26NZs2aIiopiaVKVuHDhAtLS0nipCVEVqFBpHjhwADdu3IC1tTUvmn5NMpkMHh4eiIyMxA8//MBfLKjSvXjgdO/evaWOQqR1Ktx4b775Ju7fv6/42dvbG3fv3q2UUNrOw8MDqampuHDhgtRRSAtFR0djwIABMDQ0lDoKkdapcGn++6TbmJgYPHnyROVANUHv3r1hYmKCqKgoqaOQlsnKysLx48d5g3aiKsJjqxIwNDSEs7Mzdu/eLXUU0jJ84DRR1apwacpkslLfx/H7OfE8PDxw8uRJZGZmSh2FtEh0dDTeeustNGrUSOooRFqpwncEEgQBY8aMUXxvkp+fj8DAQNSuXVtpOl7EX7YXewLR0dHw9/eXOA1pg6KiIsTGxuKTTz6ROgqR1qpwab7//vtKP48ePVrlMDVJ/fr10a1bN0RFRbE0qVIcO3YMDx8+5PeZRFWowqW5fv36ysxRI3l4eGDhwoXIz8+HkZGR1HFIw0VGRsLS0hKdO3eWOgqR1uKJQBIaPHgwnjx5gvj4eKmjkBaIjIzE4MGDee00URXi/y4Jvbg7EM+iJVVduXIF165dg6enp9RRiLQaS1NCL+4OFBUVVeq6V6LXERkZCWNjY/Tv31/qKERajaUpMd4diCpDZGQkBg4cyCcMEVUxlqbEeHcgUtW9e/dw7NgxHpolqgYsTYkZGhrCxcWFpUkVFhMTg5KSEj7VhKgasDTVgIeHB06dOsUb3lOFREVFoWvXrrC0tJQ6CpHWY2mqAVdXVwDP9xiIXkd+fj5iY2N5aJaomrA01cA/7w5E9Dri4+Px5MkTliZRNWFpqgkPDw/s3bsX+fn5UkchDRIZGQkbGxu0bdtW6ihENQJLU014eHjw7kD0WgRBQFRUFDw9PfmEIaJqwtJUE23btoW1tTUP0ZJoSUlJuH37Njw8PKSOQlRjsDTVBO8ORK8rMjIScrkcvXr1kjoKUY2hFaX58OFD+Pr6Qi6XQy6Xw9fXF9nZ2eVO/+zZM3zxxRdo164dateujUaNGsHPzw937typvtBl8PT0RFpaGpKSkiTNQZohMjISrq6u0NfXlzoKUY2hFaU5cuRIJCUlITY2FrGxsUhKSoKvr2+50+fl5eHs2bOYPXs2zp49i507d+Lq1auSn4HYu3dvyOVy7Nq1S9IcpP5u376Ns2fPSv5vlqimkQkafiwwOTkZbdq0wYkTJ9C1a1cAwIkTJ+Dk5ITLly/D1tZW1HJOnz6NLl264NatW2jatKmoeXJzcyGXy5GTkwNTU9MKv4d/Gj16NC5evMi9TXqpFStWYOLEicjMzES9evWkjkOkNqric/mfNH5P8/jx45DL5YrCBIBu3bpBLpfj2LFjopeTk5MDmUyGunXrljtNQUEBcnNzlV6VzcvLC+fOnUNKSkqlL5u0R2RkJHr27MnCJKpmGl+aGRkZqF+/fqnh9evXR0ZGhqhl5OfnY/r06Rg5cuRLfzMJDg5WfG8ql8thZWVV4dzlcXFxgaGhIX777bdKXzZph8ePH2P//v08NEskAbUtzblz50Imk730debMGQAo8xo1QRBEXbv27Nkz+Pj4oKSkBMuXL3/ptDNmzEBOTo7ilZaWVrE39xJ16tRB//79+b0mlev3339HYWEhLzUhkoCe1AHKM2HCBPj4+Lx0Gmtra5w/f77MG53fu3cPDRo0eOn8z549w3vvvYeUlBTs37//lce/DQ0NYWho+OrwKvLy8kJgYCDu378Pc3PzKl8faZaIiAi0a9cOLVq0kDoKUY2jtqVpbm4uqjCcnJyQk5ODU6dOoUuXLgCAkydPIicnB927dy93vheFee3aNRw4cABmZmaVll1VHh4e+Oijj7B7926MGTNG6jikRgoLC7F7925MnDhR6ihENZLaHp4Vy87ODoMGDUJAQABOnDiBEydOICAgAIMHD1Y6c7Z169aIiIgAABQVFWH48OE4c+YMNm/ejOLiYmRkZCAjIwOFhYVSvRUFS0tLODk58RAtlRIfH4+cnBwMHTpU6ihENZLGlyYAbN68Ge3atYOzszOcnZ3Rvn17hIWFKU1z5coV5OTkAHh+jVtkZCRu376Njh07omHDhorX65xxW5W8vLywd+9e5OXlSR2F1EhERASaNWuGjh07Sh2FqEbS+Os0pVSV1wNdu3YNrVq1wq5duzBkyJBKXTZpppKSEjRp0gTe3t4IDQ2VOg6RWuJ1mjXUm2++CTs7Ox6iJYWTJ08iPT2dh2aJJMTSVGNeXl6IiopCUVGR1FFIDezatQsWFhbo0aOH1FGIaiyWphrz8vJCVlYWjh49KnUUkpggCIiIiICnpyd0dXWljkNUY7E01ZijoyMaNWrEQ7SES5cu4dq1azw0SyQxlqYa09HRwZAhQ7Br1y4+Y7OGi4iIgImJCfr37y91FKIajaWp5ry8vHDz5k2cP39e6igkoYiICLi5ucHIyEjqKEQ1GktTzfXp0wdyuRw7d+6UOgpJ5NatWzh79iwPzRKpAZammjMwMICnpyd27NghdRSSyK5du2BgYAA3NzepoxDVeCxNDTBs2DD8+eefuHz5stRRSAIRERHo379/lVyoTUSvh6WpAZydnWFiYsK9zRro3r17OHz4MA/NEqkJlqYGMDY2hru7O7Zv3y51FKpmUVFREASBD5wmUhMsTQ0xfPhwJCUl4fr161JHoWoUERGB7t27v/LZsERUPViaGsLV1RXGxsY8RFuD5OTkYO/evRg+fLjUUYjof1iaGqJ27dpwdXXlIdoaJCoqCoWFhRg2bJjUUYjof1iaGmT48OE4ffo0bt26JXUUqgbbtm1Dt27dYGVlJXUUIvoflqYGcXd3h4GBAW90UAPk5ubi999/x7vvvit1FCL6B5amBjE1NYWLiwu/16wBdu/ejYKCAh6aJVIzLE0NM2zYMBw9ehR37tyROgpVoe3bt6NLly5o1qyZ1FGI6B9YmhrG09MTenp6iIiIkDoKVZHHjx9jz549PGuWSA2xNDVMvXr10L9/f55Fq8Wio6ORn5/P0iRSQyxNDTR8+HAcOnQImZmZUkehKrBt2zY4OjrCxsZG6ihE9C8sTQ00ZMgQAM+ffkHa5cmTJ4iJieFeJpGaYmlqIAsLC/Tp04eHaLVQTEwMnj59ytIkUlMsTQ313nvvYf/+/TxEq2W2bduGt956Cy1atJA6ChGVgaWpoV5cv8drNrVHXl4eoqOjeUMDIjXG0tRQ5ubmGDhwILZu3Sp1FKok0dHRyMvL46FZIjXG0tRgPj4+OHz4MG7fvi11FKoEW7ZsgYODA958802poxBROViaGszLywv6+vrYtm2b1FFIRTk5OYiJicGIESOkjkJEL8HS1GByuRyurq4IDw+XOgqpaNeuXSgsLIS3t7fUUYjoJViaGs7HxwcnT55ESkqK1FFIBVu2bEHPnj3RpEkTqaMQ0UuwNDWch4cHatWqxb1NDXbv3j388ccf8PHxkToKEb0CS1PD1a5dG4MHD+ZZtBrsxU0qeNYskfpjaWoBHx8fnDt3DsnJyVJHoQrYsmULBg4cCAsLC6mjENErsDS1gKurK+rUqcNDtBooLS0Nhw8f5qFZIg3B0tQCRkZGGDp0KMLDwyEIgtRx6DX8+uuvMDQ0xNChQ6WOQkQisDS1hI+PDy5fvoykpCSpo9Br2LJlC9zd3WFqaip1FCISgaWpJQYMGAALCwts2rRJ6igk0rVr15CQkMAbGhBpEJamltDX14ePjw9++eUXFBcXSx2HRNi6dStMTEzg7u4udRQiEomlqUV8fX2RkZGBffv2SR2FXkEQBGzatAlDhw6FsbGx1HGISCSWphZxdHREq1atEBYWJnUUeoVTp07h6tWr8PPzkzoKEb0GlqYWkclk8PX1xc6dO/H48WOp49BLhIWFoXHjxujbt6/UUYjoNbA0tczo0aORl5eHiIgIqaNQOQoLC7FlyxaMHj0aurq6UschotfA0tQy1tbW6NmzJw/RqrGYmBg8ePAAvr6+UkchotfE0tRCvr6+2LdvH+7cuSN1FCrDxo0b0alTJ7Rt21bqKET0mliaWmj48OHQ09PDli1bpI5C/5KVlYXdu3fzBCAiDcXS1EL16tWDh4cHD9GqofDwcJSUlPCGBkQaiqWppXx9fXHu3DlcuHBB6ij0Dxs3boSrqyvq168vdRQiqgCWppZydXWFmZkZ9zbVyJUrV3Dy5EkemiXSYCxNLWVgYAAfHx+EhYWhqKhI6jgEYNOmTZDL5fDw8JA6ChFVEEtTi40bNw4ZGRmIjY2VOkqNV1JSgrCwMHh7e8PIyEjqOERUQSxNLfbWW2+hY8eOWLdundRRaryDBw/i1q1bPDRLpOFYmlpu7Nix2L17N+7evSt1lBpt7dq1sLW1Rffu3aWOQkQqYGlquVGjRkFHR4cnBEno4cOH2LFjB8aNGweZTCZ1HCJSAUtTy73xxhsYOnQofvrpJwiCIHWcGmnz5s0oLi7moVkiLcDSrAHGjh2L5ORknDx5UuooNY4gCFizZg08PT3RoEEDqeMQkYpYmjVA//790bRpU54QJIGEhAScP38e48aNkzoKEVUClmYNoKurizFjxmDr1q148uSJ1HFqlLVr16Jx48ZwcXGROgoRVQKWZg3h7++PJ0+eIDw8XOooNcaTJ0+wZcsWjB07ls/NJNISLM0awtraGi4uLli5cqXUUWqM7du3Izc3F/7+/lJHIaJKwtKsQT7++GOcPn0aCQkJUkepEdasWYP+/fvDxsZG6ihEVElYmjWIm5sbmjRpwr3NanDu3DkcPXoUH3/8sdRRiKgSaUVpPnz4EL6+vpDL5ZDL5fD19UV2drbo+T/66CPIZDJ89913VZZRHejp6SEgIAC//PILcnJypI6j1VasWIFGjRrB09NT6ihEVIm0ojRHjhyJpKQkxMbGIjY2FklJSfD19RU1765du3Dy5Ek0atSoilOqhw8++AAFBQXYtGmT1FG0Vk5ODjZt2oQPP/wQ+vr6Uschokqk8aWZnJyM2NhYrF27Fk5OTnBycsKaNWuwe/duXLly5aXz/v3335gwYQI2b94s6sOtoKAAubm5Si9N06hRIwwZMgQrVqzgHYKqSFhYGPLz8xEQECB1FCKqZBpfmsePH4dcLkfXrl0Vw7p16wa5XI5jx46VO19JSQl8fX3x+eefo23btqLWFRwcrDgELJfLYWVlpXJ+KQQGBuLPP//E0aNHpY6idQRBwPLlyzF06NAac/SCqCbR+NLMyMhA/fr1Sw2vX78+MjIyyp1v8eLF0NPTw8SJE0Wva8aMGcjJyVG80tLSKpRZav3790fLli2xYsUKqaNonYMHDyI5ORnjx4+XOgoRVQG1Lc25c+dCJpO99HXmzBkAKPPJEYIglPtEiYSEBHz//ffYsGHDaz11wtDQEKampkovTaSjo4PAwEBs27btpb9Y0Otbvnw5WrdujT59+kgdhYiqgNqW5oQJE5CcnPzSl729PSwtLct8VuS9e/fKvUH24cOHkZmZiaZNm0JPTw96enq4desWpkyZAmtr6yp+Z+ph3LhxMDAw4N5mJUpPT0dERATGjx/PR4ARaSk9qQOUx9zcHObm5q+czsnJCTk5OTh16hS6dOkCADh58iRycnLKfeCvr68vBgwYoDTMxcUFvr6+NebuLXXr1sWYMWOwcuVKfPnllzA0NJQ6ksZbsWIFDA0N+QgwIi2mtnuaYtnZ2WHQoEEICAjAiRMncOLECQQEBGDw4MGwtbVVTNe6dWtEREQAAMzMzGBvb6/00tfXh6WlpdI82u7TTz9FZmYmtm7dKnUUjff06VOsWLECY8eOhVwulzoOEVURjS9N4PlDftu1awdnZ2c4Ozujffv2CAsLU5rmypUrvKD/X2xtbeHm5obvvvuOl5+oaNOmTcjKysKkSZOkjkJEVUgm8NOywnJzcyGXy5GTk6OxJwXt3bsXLi4uOHjwIHr16iV1HI1UUlICe3t72NraKo5mEJE0qvpzWSv2NKniBg4cCDs7O62/hWBV+v3335GcnIygoCCpoxBRFWNp1nAymQyTJk3Cb7/9hpSUFKnjaKTQ0FA4Ojri7bffljoKEVUxlibB19cXdevW5d5mBVy4cAFxcXH47LPPeJkJUQ3A0iTUqlULn376KdauXYv79+9LHUejhIaGonHjxnj33XeljkJE1YClSQCe30wCAH744QeJk2iOv//+G5s3b8bEiRP5NBOiGoKlSQCe30wiICAAP/zwAx4/fix1HI3w7bffwtjYGIGBgVJHIaJqwtIkhSlTpuDRo0dYvXq11FHU3v3797Fq1Sp8+umnGnu5ERG9PpYmKVhZWWH06NH49ttvUVBQIHUctfbipCnezICoZmFpkpJp06YhPT0dmzdvljqK2srJycGPP/6IwMBAUfdHJiLtwdIkJXZ2dvDy8sLixYtRXFwsdRy1tGzZMjx9+hRTpkyROgoRVTOWJpXy5Zdf4urVq9iyZYvUUdTOkydPEBoairFjx6JRo0ZSxyGiasbSpFIcHR3h6emJefPmoaioSOo4amX16tV4+PAhpk2bJnUUIpIAS5PKNH/+fPz111/YuHGj1FHUxpMnT7Bo0SL4+fnBxsZG6jhEJAGWJpWpQ4cOGD58OObPn4/CwkKp46iFpUuX4uHDh/jqq6+kjkJEEmFpUrnmzp2L1NRUrFu3TuooksvOzsaSJUvw4YcfwtraWuo4RCQRliaVq23bthg5ciQWLlyI/Px8qeNIKiQkBPn5+Zg5c6bUUYhIQixNeqk5c+YgIyMDK1eulDqKZO7du4fQ0FBMmDABDRs2lDoOEUmIpUkv9eabb8Lf3x8LFixAdna21HEkMW/ePOjq6uKLL76QOgoRSYylSa80f/585OfnY+HChVJHqXaXL1/GypUrMXPmTN79h4hYmvRqDRs2xLRp07B06VKkpKRIHadaTZs2DVZWVvj000+ljkJEaoClSaJMmTIF5ubmmD59utRRqs2BAwcQFRWF4OBgGBkZSR2HiNQAS5NEqV27NhYuXIhff/0VR48elTpOlSsuLkZQUBC6du0Kb29vqeMQkZpgaZJofn5+cHR0xPjx47X+9norVqxAUlISvvvuO8hkMqnjEJGaYGmSaDo6OlixYgUuXLiAZcuWSR2nyty9exezZs1CQEAAunXrJnUcIlIjLE16LY6OjggMDMTs2bORnp4udZwq8fnnn0NPTw/BwcFSRyEiNcPSpNe2cOFCGBkZaeXzJA8ePIiwsDAsXrwYZmZmUschIjXD0qTXVq9ePXzzzTfYsmULoqOjpY5TaZ4+fYoPP/wQTk5O8Pf3lzoOEakhliZViJ+fHwYNGoQPP/wQDx8+lDpOpZg9ezZu3bqFn376CTo6/K9BRKXxk4EqRCaTYfXq1Xj8+DGCgoKkjqOyY8eOISQkBP/5z3/QunVrqeMQkZpiaVKFWVlZISQkBBs2bNDow7RPnz6Fv78/unTpohW/ABBR1WFpkkrGjh2LQYMGYdy4ccjIyJA6ToVMmTIFt27dwvr166Grqyt1HCJSYyxNUolMJsP69esBAL6+vigpKZE40evZtm0bVqxYge+++w52dnZSxyEiNcfSJJVZWlpi06ZN2LdvHxYtWiR1HNFSUlLwwQcf4N1338VHH30kdRwi0gAsTaoUAwYMwJdffomvvvoKhw4dkjrOK+Xn58PHxwdmZmZYs2YNb5VHRKKwNKnSzJ07F7169cKwYcPU+hFigiDgww8/xPnz5xEeHg65XC51JCLSECxNqjR6enrYtm0b5HI5PD098ejRI6kjlWnJkiUICwvDTz/9hM6dO0sdh4g0CEuTKpWZmRmioqKQmpqKkSNHqt3TUHbt2oUZM2Zg1qxZGDFihNRxiEjDsDSp0tnZ2SE8PBx79uzBuHHj1OaM2n379sHHxwfDhg3DvHnzpI5DRBqIpUlVYtCgQdi0aRPCwsLwySefQBAESfMcPXoUnp6e6Nu3LzZt2sTb5BFRhehJHYC0l4+PD/Ly8jBu3DgYGhoiJCREkrI6ceIE3Nzc0LlzZ+zYsQOGhobVnoGItANLk6rU2LFjUVBQgE8++QSZmZlYv359tZZWVFQUvL294eDggKioKNSqVava1k1E2ofHqKjKffzxx/j111+xc+dOuLq6Ijs7u1rWu2rVKnh5ecHNzQ1xcXGoU6dOtayXiLQXS5OqxfDhwxEXF4fExEQ4OjoiMTGxytb14rmYgYGB+OSTTxAeHg4jI6MqWx8R1RwsTao2PXv2xJkzZ2BqagonJyd88803lX5JyunTp9GlSxeEhYVh3bp1WLp0KW/CTkSVhqVJ1apFixY4duwYPvnkE0yfPh1du3bFwYMHVV7uvXv3MHHiRHTr1g36+vo4deoUxo4dWwmJiYj+H0uTqp2RkRG+/fZbHDt2DDo6OujTpw9cXFywd+/e1740JS0tDdOnT0eLFi3w888/Y9GiRTh16hTatWtXRemJqCaTCVJfQKfBcnNzIZfLkZOTA1NTU6njaKSSkhJs27YNixYtQlJSEqytrfHOO+9gwIAB6Ny5M8zNzZWmLywsxMWLF3H06FHs2rUL8fHxMDExwYcffogvvvii1PREVLNU9ecyS1MFLM3KIwgCjhw5gl9++QW7du1SPNC6Tp06qF+/PnR0dPDo0SPcvXsXgiBAT08P/fr1w/Dhw+Hj48MzY4kIAEtTrbE0q4YgCLhx4wbOnDmDtLQ03Lt3D4IgoFatWmjcuDHatGmDTp06wdjYWOqoRKRmqvpzmTc3ILUjk8nQokULtGjRQuooRERKeCIQERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpG0ojQfPnwIX19fyOVyyOVy+Pr6Ijs7+5XzJScnw9PTE3K5HHXq1EG3bt2Qmppa9YGJiEgjaUVpjhw5EklJSYiNjUVsbCySkpLg6+v70nmuX7+Ot99+G61bt0Z8fDzOnTuH2bNnw8jIqJpSExGRptH4p5wkJyejTZs2OHHiBLp27QoAOHHiBJycnHD58mXY2tqWOZ+Pjw/09fURFhZW4XXzKSdEROqFTzl5hePHj0MulysKEwC6desGuVyOY8eOlVmaJSUliI6OxrRp0+Di4oLExETY2NhgxowZ8PLyKnddBQUFKCgoUPyck5MD4PlfEhERSe/F53FV7Q9qfGlmZGSgfv36pYbXr19f8SDjf8vMzMTjx4+xaNEiLFiwAIsXL0ZsbCzeeecdHDhwAL179y5zvuDgYMybN6/UcCsrK9XeBBERVaqsrCzI5fJKX67alubcuXPLLKh/On36NIDnz1/8N0EQyhwOPN/TBIAhQ4bgs88+AwB07NgRx44dw8qVK8stzRkzZiAoKEjxc3Z2Npo1a4bU1NQq+cupbLm5ubCyskJaWprGHE5m5urBzNWDmateTk4OmjZtijfeeKNKlq+2pTlhwgT4+Pi8dBpra2ucP38ed+/eLTXu3r17aNCgQZnzmZubQ09PD23atFEabmdnhyNHjpS7PkNDQxgaGpYaLpfLNeIf0wumpqYalRdg5urCzNWDmauejk7VnOeqtqVpbm4Oc3PzV07n5OSEnJwcnDp1Cl26dAEAnDx5Ejk5OejevXuZ8xgYGKBz5864cuWK0vCrV6+iWbNmqocnIiKtpPGXnNjZ2WHQoEEICAjAiRMncOLECQQEBGDw4MFKJwG1bt0aERERip8///xzhIeHY82aNfjrr7/w448/IioqCuPHj5fibRARkQbQ+NIEgM2bN6Ndu3ZwdnaGs7Mz2rdvX+pSkitXrijOdgWAoUOHYuXKlViyZAnatWuHtWvXYseOHXj77bdFr9fQ0BBz5swp85CtOtK0vAAzVxdmrh7MXPWqOq/GX6dJRERUXbRiT5OIiKg6sDSJiIhEYmkSERGJxNIkIiISiaVZQcuXL4eNjQ2MjIzg4OCAw4cPSx1JYe7cuZDJZEovS0tLxXhBEDB37lw0atQIxsbG6NOnD/78889qzXjo0CF4eHigUaNGkMlk2LVrl9J4MRkLCgrw6aefwtzcHLVr14anpydu374tSd4xY8aU2ubdunWTLC/w/LaPnTt3Rp06dVC/fn14eXmVujZZ3bazmMzqtK1XrFiB9u3bKy78d3Jywp49exTj1W37ismsTtu3PMHBwZDJZJg8ebJiWLVta4Fe29atWwV9fX1hzZo1wqVLl4RJkyYJtWvXFm7duiV1NEEQBGHOnDlC27ZthfT0dMUrMzNTMX7RokVCnTp1hB07dggXLlwQvL29hYYNGwq5ubnVljEmJkaYOXOmsGPHDgGAEBERoTReTMbAwEChcePGQlxcnHD27Fmhb9++QocOHYSioqJqz/v+++8LgwYNUtrmWVlZStNUZ15BEAQXFxdh/fr1wsWLF4WkpCTB3d1daNq0qfD48WPFNOq2ncVkVqdtHRkZKURHRwtXrlwRrly5Inz55ZeCvr6+cPHiRUEQ1G/7ismsTtu3LKdOnRKsra2F9u3bC5MmTVIMr65tzdKsgC5dugiBgYFKw1q3bi1Mnz5dokTK5syZI3To0KHMcSUlJYKlpaWwaNEixbD8/HxBLpcLK1eurKaEyv5dQmIyZmdnC/r6+sLWrVsV0/z999+Cjo6OEBsbW615BeH5B82QIUPKnUfKvC9kZmYKAISDBw8KgqD+27mszIKg/tu6Xr16wtq1azVi+/47syCo9/Z99OiR8OabbwpxcXFC7969FaVZnduah2dfU2FhIRISEuDs7Kw03NnZGceOHZMoVWnXrl1Do0aNYGNjAx8fH9y4cQMAkJKSgoyMDKX8hoaG6N27t9rkF5MxISEBz549U5qmUaNGsLe3l+x9xMfHo379+mjVqhUCAgKQmZmpGKcOeV/c3OPFjaw1YTv/O/ML6riti4uLsXXrVjx58gROTk4asX3/nfkFddy+APDJJ5/A3d0dAwYMUBpendtabe89q67u37+P4uLiUjeDb9CgQbmPIqtuXbt2xcaNG9GqVSvcvXsXCxYsQPfu3fHnn38qMpaV/9atW1LELUVMxoyMDBgYGKBevXqlppHi78HV1RXvvvsumjVrhpSUFMyePRv9+vVDQkICDA0NJc8rCAKCgoLw9ttvw97eHoD6b+eyMgPqt60vXLgAJycn5Ofnw8TEBBEREWjTpo3ig1gdt295mQH1274vbN26FWfPnlU83eqfqvPfMkuzgv792DHhJY8iq26urq6KP7dr1w5OTk5o0aIFfv75Z8UX+uqc/4WKZJTqfXh7eyv+bG9vD0dHRzRr1gzR0dF45513yp2vuvJOmDAB58+fL/MpPuq6ncvLrG7b2tbWFklJScjOzsaOHTvw/vvv4+DBg4rx6rh9y8vcpk0btdu+AJCWloZJkyZh7969MDIyKne66tjWPDz7mszNzaGrq1vqN5PMzMxyH0Umtdq1a6Ndu3a4du2a4ixadc4vJqOlpSUKCwvx8OHDcqeRUsOGDdGsWTNcu3YNgLR5P/30U0RGRuLAgQNo0qSJYrg6b+fyMpdF6m1tYGCAli1bwtHREcHBwejQoQO+//57td6+5WUui9TbF3h+aDUzMxMODg7Q09ODnp4eDh48iKVLl0JPT0+x3urY1izN12RgYAAHBwfExcUpDY+Liyv3UWRSKygoQHJyMho2bAgbGxtYWloq5S8sLMTBgwfVJr+YjA4ODtDX11eaJj09HRcvXlSL95GVlYW0tDQ0bNgQgDR5BUHAhAkTsHPnTuzfvx82NjZK49VxO78qc1nUYVv/kyAIKCgoUMvt+6rMZVGH7du/f39cuHABSUlJipejoyNGjRqFpKQkNG/evPq2dQVOYKrxXlxysm7dOuHSpUvC5MmThdq1aws3b96UOpogCIIwZcoUIT4+Xrhx44Zw4sQJYfDgwUKdOnUU+RYtWiTI5XJh586dwoULF4QRI0ZU+yUnjx49EhITE4XExEQBgBASEiIkJiYqLtsRkzEwMFBo0qSJ8Mcffwhnz54V+vXrV2Wnvb8s76NHj4QpU6YIx44dE1JSUoQDBw4ITk5OQuPGjSXLKwiC8PHHHwtyuVyIj49XunwgLy9PMY26bedXZVa3bT1jxgzh0KFDQkpKinD+/Hnhyy+/FHR0dIS9e/cKgqB+2/dVmdVt+77MP8+eFYTq29YszQpatmyZ0KxZM8HAwEDo1KmT0inxUntxfZK+vr7QqFEj4Z133hH+/PNPxfiSkhJhzpw5gqWlpWBoaCj06tVLuHDhQrVmPHDggACg1Ov9998XnfHp06fChAkThDfeeEMwNjYWBg8eLKSmplZ73ry8PMHZ2VmwsLAQ9PX1haZNmwrvv/9+qSzVmVcQhDLzAhDWr1+vmEbdtvOrMqvbth47dqzic8DCwkLo37+/ojAFQf2276syq9v2fZl/l2Z1bWs+GoyIiEgkfqdJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSTWeTCbDrl27qnw91tbW+O6779RmOZq27leJj4+HTCZDdna21FFIi7E0SatlZmbio48+QtOmTWFoaAhLS0u4uLjg+PHjimnS09OVHqemLjZs2IC6deuWGn769Gl8+OGHWrtuInXG52mSVhs2bBiePXuGn3/+Gc2bN8fdu3exb98+PHjwQDHNi0c4aQoLC4sauW4idcA9TdJa2dnZOHLkCBYvXoy+ffuiWbNm6NKlC2bMmAF3d3fFdP88PHvz5k3IZDL8+uuv6NmzJ4yNjdG5c2dcvXoVp0+fhqOjI0xMTDBo0CDcu3dPsYw+ffpg8uTJSuv38vLCmDFjys0XEhKCdu3aoXbt2rCyssL48ePx+PFjAM8PNfr7+yMnJwcymQwymQxz584FUPoQaWpqKoYMGQITExOYmprivffew927dxXj586di44dOyIsLAzW1taQy+Xw8fHBo0ePysz1OuuWyWRYtWoVBg8ejFq1asHOzg7Hjx/HX3/9hT59+qB27dpwcnLC9evXldYRFRUFBwcHGBkZoXnz5pg3bx6KiorK3VZliYmJQatWrWBsbIy+ffvi5s2bSuOzsrIwYsQINGnSBLVq1UK7du2wZcsWxfiNGzfCzMys1COxhg0bBj8/PwDAuXPn0LdvX9SpUwempqZwcHDAmTNnXisnaReWJmktExMTmJiYYNeuXeU+K7A8c+bMwaxZs3D27Fno6elhxIgRmDZtGr7//nscPnwY169fx1dffaVSPh0dHSxduhQXL17Ezz//jP3792PatGkAgO7du+O7776Dqakp0tPTkZ6ejqlTp5ZahiAI8PLywoMHD3Dw4EHExcXh+vXr8Pb2Vpru+vXr2LVrF3bv3o3du3fj4MGDWLRoUZm5xK77hf/85z/w8/NDUlISWrdujZEjR+Kjjz7CjBkzFAUzYcIExfS///47Ro8ejYkTJ+LSpUtYtWoVNmzYgIULF4redmlpaXjnnXfg5uaGpKQkfPDBB5g+fbrSNPn5+XBwcMDu3btx8eJFfPjhh/D19cXJkycBAO+++y6Ki4sRGRmpmOf+/fvYvXs3/P39AQCjRo1CkyZNcPr0aSQkJGD69OnQ19cXnZO0kMrPZyFSY9u3bxfq1asnGBkZCd27dxdmzJghnDt3TmkaAEJERIQgCIKQkpIiABDWrl2rGL9lyxYBgLBv3z7FsODgYMHW1lbx878fUyQIgjBkyBDFo84EQRCaNWsmhIaGlpv1119/FczMzBQ/r1+/XpDL5aWm++dy9u7dK+jq6io93ujPP/8UAAinTp0SBEEQ5syZI9SqVUvpuYKff/650LVr13KziFm3IDzfdrNmzVL8fPz4cQGAsG7dOsWwLVu2CEZGRoqfe/bsKXz99ddKyw0LCxMaNmxYbp5/mzFjhmBnZyeUlJQohn3xxRcCAOHhw4flzufm5iZMmTJF8fPHH38suLq6Kn7+7rvvhObNmyuWW6dOHWHDhg2ic5H2454mabVhw4bhzp07iIyMhIuLC+Lj49GpUyds2LDhpfO1b99e8ecGDRoAANq1a6c0LDMzU6VsBw4cwMCBA9G4cWPUqVMHfn5+yMrKwpMnT0QvIzk5GVZWVrCyslIMa9OmDerWrYvk5GTFMGtra9SpU0fxc8OGDVXO/4KYbZWfn4/c3FwAQEJCAubPn684EmBiYoKAgACkp6cjLy9P1DqTk5PRrVs3yGQyxTAnJyelaYqLi7Fw4UK0b98eZmZmMDExwd69e5GamqqYJiAgAHv37sXff/8NAFi/fj3GjBmjWG5QUBA++OADDBgwAIsWLSp1mJlqHpYmaT0jIyMMHDgQX331FY4dO4YxY8Zgzpw5L53nn4fgXnyA/ntYSUmJ4mcdHR0I/3o07bNnz8pd/q1bt+Dm5gZ7e3vs2LEDCQkJWLZs2Svn+zdBEJSKo7zh/z6k+O/8qhCzrQAo1ldSUoJ58+YhKSlJ8bpw4QKuXbsGIyMjUev897Yuy7fffovQ0FBMmzYN+/fvR1JSElxcXFBYWKiY5q233kKHDh2wceNGnD17FhcuXFD6Hnru3Ln4888/4e7ujv3796NNmzaIiIgQlZG0E8+epRqnTZs2lX5dpoWFBdLT0xU/FxcX4+LFi+jbt2+Z0585cwZFRUX49ttvoaPz/HfXX3/9VWkaAwMDFBcXv3S9bdq0QWpqKtLS0hR7m5cuXUJOTg7s7Owq/H7ErLuiOnXqhCtXrqBly5YVXkZZf4cnTpxQ+vnw4cMYMmQIRo8eDeB5WV+7dq3Udvnggw8QGhqKv//+GwMGDFDaaweAVq1aoVWrVvjss88wYsQIrF+/HkOHDq1wdtJs3NMkrZWVlYV+/fph06ZNOH/+PFJSUrBt2zYsWbIEQ4YMqdR19evXD9HR0YiOjsbly5cxfvz4l15k36JFCxQVFeGHH37AjRs3EBYWhpUrVypNY21tjcePH2Pfvn24f/9+mYcuBwwYgPbt22PUqFE4e/YsTp06BT8/P/Tu3RuOjo4Vfj9i1l1RX331FTZu3KjYi0tOTkZ4eDhmzZolehmBgYG4fv06goKCcOXKFfzyyy+lDrm3bNkScXFxOHbsGJKTk/HRRx8hIyOj1LJGjRqFv//+G2vWrMHYsWMVw58+fYoJEyYgPj4et27dwtGjR3H69GmVfhkhzcfSJK1lYmKCrl27IjQ0FL169YK9vT1mz56NgIAA/Pjjj5W6rrFjx+L9999XFJaNjU25e5kA0LFjR4SEhGDx4sWwt7fH5s2bERwcrDRN9+7dERgYCG9vb1hYWGDJkiWllvPicpl69eqhV69eGDBgAJo3b47w8HCV3o+YdVeUi4sLdu/ejbi4OHTu3BndunVDSEgImjVrpphmzJgx6NOnT7nLaNq0KXbs2IGoqCh06NABK1euxNdff600zezZs9GpUye4uLigT58+sLS0hJeXV6llmZqaYtiwYTAxMVEar6uri6ysLPj5+aFVq1Z477334Orqinnz5qm6CUiDyQQxXw4QEVWjPn36oE+fPorrQ6vawIEDYWdnh6VLl1bL+khzsTSJSK08evQIbdq0QXJyMkxMTKp0XQ8ePMDevXsxatQoXLp0Cba2tlW6PtJ8PBGIiNRKnTp1kJaWVi3r6tSpEx4+fIjFixezMEkU7mkSERGJxBOBiIiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCL9H7hJf86hwyz8AAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 450x500 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot(flux);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "f5aa3a2d-f85f-465a-97c9-69bcf83a8529",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAc0AAAH/CAYAAAA194a7AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAXjZJREFUeJzt3XlcTQnjBvDnti90UUJUGCSyVZYYO6Ga7LKVZZjB27zIMtaxjBHesYyxjiWjQYwlW1IvlX1LWRNDFEq2ylqq8/tjXvc3V8VRt869t+f7+dzPR2d97sF9OueeRSYIggAiIiL6JB2pAxAREWkKliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmlXqbNm2CTCbL9zVx4kTFdNWrV8fQoUNLPN/s2bMhk8lKZF13796FTCbDpk2bSmR9UmjXrh3atWunNEwmk2H27NmS5CHNoid1ACJ1ERAQgLp16yoNs7KykijN/xsxYgS6du1aIuuqUqUKTp8+jS+++KJE1kekaViaRP/j4OAAZ2dnqWPkUa1aNVSrVq1Y15GTk4Ps7GwYGhqiRYsWxbqu/Lx58wZGRkYltkdNVFg8PEtUSO8P6969e1dpeGRkJGQyGSIjIwEAt27dgpmZGfr27as03dGjR6Grq4uZM2d+dD35HZ6tXr06PDw8EBoaCkdHRxgbG6Nu3brYuHHjJ3O/PwS7aNEizJs3DzVq1IChoSEiIiLyHJ4NDg6GTCbDkSNH8ixn9erVkMlkuHz5smLYhQsX4OnpiQoVKsDIyAhNmjTBjh07lOZ7v93CwsIwfPhwVKxYESYmJsjMzCwwc1paGiZMmICaNWvC0NAQlpaWcHNzw40bNxTTZGVlYd68eahbty4MDQ1RsWJFDBs2DI8fP/7kNiESi6VJ9D/v97b++VKF2rVrY926ddi5cyeWL18OAEhJScHAgQPRunXrQn+XdunSJUyYMAHjx4/H3r170bBhQ3z99dc4duyYqPmXL1+Oo0eP4ueff8ahQ4fyHJoGAA8PD1haWiIgICDPuE2bNsHR0RENGzYEAERERKBVq1ZIS0vDmjVrsHfvXjRu3BheXl75fkc6fPhw6OvrIzAwEDt37oS+vn6+OV+8eIEvv/wSa9euxbBhw7B//36sWbMGderUQXJyMgAgNzcX3bt3x4IFCzBw4EAcPHgQCxYsQHh4ONq1a4c3b96I2iZEn8LDs0T/k99hyXfv3kFPr+j/Tby8vBAVFYVJkyahWbNmmD59OgRBwLZt26Crq1uoZT558gQnT56EjY0NAKBNmzY4cuQItm7dijZt2nxyfiMjIxw+fFiprD7ca9bT08PgwYOxevVqpKenQy6XAwDi4uJw7tw5/Prrr4ppx4wZg/r16+Po0aOKbdalSxc8efIE06ZNg4+PD3R0/v/39I4dO2Lt2rWfzLls2TJcu3YN4eHh6NSpk2J4r169FH/esWMHQkNDsWvXLqXhjRo1QtOmTbFp0yaMHj36k+si+hTuaRL9z+bNm3H+/HmllyoK872lS5eifv36aN++PSIjI/HHH3+gSpUqhV5e48aNFYUJ/F2CderUwb1790TN7+npWeDe3T8NHz4cb968wfbt2xXDAgICYGhoiIEDBwIA/vrrL9y4cQODBg0CAKW9dTc3NyQnJyM+Pl5pub179xaV89ChQ6hTp45SYX7owIEDKFeuHL766iuldTdu3BiVK1dWHConKiqWJtH/2Nvbw9nZWemlSu9L5u3bt2jcuDE6d+5cpOWZm5vnuw6xhyLFFnb9+vXRtGlTxSHanJwc/PHHH+jevTsqVKgAAHj06BEAYOLEidDX11d6jRkzBsDfe8aFWf/jx48/eSLUo0ePkJaWBgMDgzzrT0lJybNuosLi4VmiQjIyMgKAPCewFPQBffXqVfzwww9o2rQpzp8/jyVLlsDPz6/Ycxbkc85UHTZsGMaMGYO4uDjcuXMHycnJGDZsmGK8hYUFAGDq1KlKh0f/yc7OrlDrr1ixIu7fv//RaSwsLGBubo7Q0NB8x5ctW1bUuog+haVJVEjVq1cHAFy+fFmpEPbt25dn2levXqFv376oXr06IiIiMGXKFEyZMgWtWrVC8+bNSypyoQ0YMAB+fn7YtGkT7ty5g6pVq8LV1VUx3s7ODrVr18alS5cwf/58la67W7du+OGHH3D06FF06NAh32k8PDwQFBSEnJwcjdiepLlYmkSF1LRpU9jZ2WHixInIzs5G+fLlsWfPHpw4cSLPtKNGjUJiYiLOnTsHU1NTLF68GKdPn0b//v0RExODcuXKlfwb+AzlypVDz549sWnTJqSlpWHixIlKJ/UAwNq1a9GtWzd06dIFQ4cORdWqVfHs2TPExcXh4sWL+PPPPwu17nHjxmH79u3o3r07pkyZgmbNmuHNmzeIioqCh4cH2rdvj/79+2PLli1wc3PD2LFj0axZM+jr6+P+/fuIiIhA9+7d0bNnT1VsCirl+J0mUSHp6upi//79qFu3LkaNGgUfHx8YGhpixYoVStOtX78ef/zxB1auXIn69esDAAwMDLB9+3Y8e/ZM6TCnOhs2bBhSU1ORlZWV7+0E27dvj3PnzqFcuXIYN24cOnXqhNGjR+O///3vR0/i+ZSyZcvixIkT+Prrr/Hbb7/B3d0dI0eORHx8vOKOTbq6uti3bx+mTZuG3bt3o2fPnujRowcWLFgAIyMjNGjQoNDrJ/onmSAIgtQhiIiINAH3NImIiERiaRIREYnE0iQiIhJJq0pz1apVqFGjBoyMjODk5ITjx49/dPrMzExMnz4dtra2MDQ0xBdffCHqhtdERFQ6ac0lJ9u3b8e4ceOwatUqtGrVSnH6+/Xr15VuNfZP/fr1w6NHj7BhwwbUqlULqampKrtJNxERaR+tOXu2efPmcHR0xOrVqxXD7O3t0aNHD/j7++eZPjQ0FP3798edO3cUtwIjIiL6GK3Y08zKykJ0dDSmTJmiNNzV1RWnTp3Kd559+/bB2dkZixYtQmBgIExNTeHp6Ykff/wRxsbG+c6TmZmpdMu03NxcPHv2DObm5nx4LhGRGhAEAS9evICVlVWeG3CoglaU5pMnT5CTk4NKlSopDa9UqRJSUlLynefOnTs4ceIEjIyMsGfPHjx58gRjxozBs2fPCvxe09/fH3PmzFF5fiIiUq2kpKRP3ui/MLSiNN/7cG9PEIQC9wBzc3Mhk8mwZcsWxTMClyxZgj59+mDlypX57m1OnTpV6Qbb6enpsLGxQVJSEszMzFT4ToiIqDAyMjJgbW1dbDfp14rStLCwgK6ubp69ytTU1Dx7n+9VqVIFVatWVRQm8Pd3oIIg4P79+6hdu3aeeQwNDWFoaJhnuJmZGUuTiEiNFNdXZlpxyYmBgQGcnJwQHh6uNDw8PBwtW7bMd55WrVrh4cOHePnypWLYzZs3oaOjUyy79EREpPm0ojQBwM/PD+vXr8fGjRsRFxeH8ePHIzExEaNGjQLw96FVHx8fxfQDBw6Eubk5hg0bhuvXr+PYsWOYNGkShg8fXuCJQEREVLppxeFZAPDy8sLTp08xd+5cJCcnw8HBASEhIbC1tQUAJCcnIzExUTF9mTJlEB4eju+++w7Ozs4wNzdHv379MG/ePKneAhERqTmtuU5TChkZGZDL5UhPT+d3mkQaQBAEZGdnIycnR+ooVAT6+vrQ1dXNd1xxfy5rzZ4mEdHHZGVlITk5Ga9fv5Y6ChWRTCZDtWrVUKZMmRJfN0uTiLRebm4uEhISoKurCysrKxgYGPCGJBpKEAQ8fvxYcZVDQXucxYWlSURaLysrC7m5ubC2toaJiYnUcaiIKlasiLt37+Ldu3clXppac/YsEdGnFMdt1ajkSXmUgP+CiIiIRGJpEhERicTSJCLScDKZDMHBwZKs++TJk2jQoAH09fXRo0cPREZGQiaTIS0tTZI8xY2lSWpLEAQ8f/4cN2/exM2bN5GUlMSHhFOpMnToUMhkMshkMujp6cHGxgajR4/G8+fPlaZLTk5Gt27dJMno5+eHxo0bIyEhAZs2bZIkQ0ni2bOkNgRBQHR0NIKDg3HixAlER0cr3RsY+PtEjjp16qBFixbo1q0b3NzcJLlWi6ikdO3aFQEBAcjOzsb169cxfPhwpKWlYdu2bYppKleuXKwZcnJyIJPJ8j2R6vbt2xg1alSpuWc39zRJcq9evcKKFStQt25dNG3aFGvWrIGFhQVmzpyJP//8E5GRkTh27BgOHjyI1atXo0OHDrh48SK8vLxgaWmJb7/9FtevX5f6bRAVC0NDQ1SuXBnVqlWDq6srvLy8EBYWpjTNPw/Puri4YMqUKUrjHz9+DH19fURERAD4+xKcyZMno2rVqjA1NUXz5s0RGRmpmH7Tpk0oV64cDhw4gHr16sHQ0BD37t1TWubdu3chk8nw9OlTDB8+HDKZLN89zdmzZ6Nx48ZKw5YtW4bq1asDAN6+fYv69evjm2++UYxPSEiAXC7HunXrPmNLlQzuaZJksrOzsW7dOsyZMwdPnjxB7969sXLlSrRr1w56ep/+p3nnzh1s3boVq1atwrp16+Dt7Y0ff/wRNjY2JZCetMHr169x48aNEl1n3bp1C32t6J07dxAaGgp9ff0Cpxk0aBD+85//wN/fX3Fpxvbt21GpUiW0bdsWADBs2DDcvXsXQUFBsLKywp49e9C1a1dcuXJF8VjE169fw9/fH+vXr4e5uTksLS2V1mNtbY3k5GTY2dlh7ty58PLyglwux9mzZz/rPRkZGWHLli1o3rw53Nzc8NVXX8Hb2xvt27fHyJEjP2tZJUKgQktPTxcACOnp6VJH0TixsbGCo6OjIJPJBG9vbyEhIaHQy8rMzBRWr14tWFpaCqampsIvv/wiZGdnqy4sabw3b94I169fF968eaM0PDo6WgBQoq/o6GjRuYcMGSLo6uoKpqamgpGRkWIZS5YsUZoOgLBnzx5BEAQhNTVV0NPTE44dO6YY7+LiIkyaNEkQBEH466+/BJlMJjx48EBpGR07dhSmTp0qCIIgBAQECACE2NjYT2aUy+VCQECA4ueIiAgBgPD8+XNBEARh1qxZQqNGjZTmWbp0qWBra6s0bNGiRYKFhYXw3XffCZUrVxYeP35c4DoL+vsUhOL/XOaeJpUoQRCwbNkyTJ48GXXr1sWZM2fQrFmzIi3TwMAAo0aNwsCBAzFt2jSMHTsWwcHB2Lp1a7F/10OarW7duoiOji7xdX6O9u3bY/Xq1Xj9+jXWr1+Pmzdv4rvvvitw+ooVK6Jz587YsmULWrdujYSEBJw+fRqrV68GAFy8eBGCIKBOnTpK82VmZsLc3Fzxs4GBARo2bPhZWYtiwoQJ2Lt3L3799VccOnQIFhYWJbbuz8HSpBLz6tUr+Pj4YPfu3ZgwYQLmz58PAwMDlS3fzMwMK1asQO/evTFw4EA0adIEu3fvhouLi8rWQdrFxMQEjo6OUsf4KFNTU9SqVQsAsHz5crRv3x5z5szBjz/+WOA8gwYNwtixY/Hrr79i69atqF+/Pho1agTg7/vw6urqIjo6Os8t6P55Up2xsbFK7ryjo6MD4YOHab179y7PdKmpqYiPj4euri5u3bqFrl27FnndxYEnAlGJePz4MTp27IjDhw9jz549+Pnnn1VamP/Uvn17xMbGolatWujQoQP27t1bLOshksKsWbPw888/4+HDhwVO06NHD7x9+xahoaHYunUrBg8erBjXpEkT5OTkIDU1FbVq1VJ6FceRmYoVKyIlJUWpOGNjY/NMN3z4cDg4OGDz5s2YPHmy2p7cx9KkYnfnzh20atUKCQkJiIyMRI8ePYp9nZUqVUJ4eDg8PDzQq1cvrFmzptjXSVQS2rVrh/r162P+/PkFTmNqaoru3btj5syZiIuLw8CBAxXj6tSpg0GDBimO+iQkJOD8+fNYuHAhQkJCiiXv48ePsWjRIty+fRsrV67EoUOHlKZZuXIlTp8+jc2bN2PgwIHo06cPBg0ahKysLJXnKSqWJhWre/fuoV27dhAEAadPn4azs3OJrdvIyAhBQUHw9fXF6NGjsXLlyhJbN1Fx8vPzw7p165CUlFTgNIMGDcKlS5fQunXrPGeUBwQEwMfHBxMmTICdnR08PT1x9uxZWFtbqzyrvb09Vq1ahZUrV6JRo0Y4d+4cJk6cqBh/48YNTJo0CatWrVKsf+XKlUhLS8PMmTNVnqeoZMKHB5tJtOJ+Qrime/jwIVq3bg0AOHbsGKpWrSpJDkEQMGHCBCxduhQbN27EsGHDJMlB0nn79i0SEhJQo0YNGBkZSR2Hiuhjf5/F/bnME4GoWDx58gSdOnXCu3fvJC1M4O8LvxcvXozXr19jxIgRMDU1Rb9+/STLQ0Sai6VJKpeVlYVevXrhyZMnOHHihOLOH1KSyWRYtWoVXrx4AR8fH1hbW/OsWiL6bPxOk1RKEASMHj0aZ8+eRXBwcJ5rwaSko6ODjRs3wtnZGT169EBiYqLUkYhIw7A0SaXef2+4bt06tGzZUuo4eRgaGmL37t0wNjaGp6dnnhvCExF9DEuTVObYsWOYNGkSJk2aBB8fH6njFMjS0hL79+9XPJ2B58KVHvy71g5S/j2yNEklHj9+jAEDBqB169bw9/eXOs4nNWjQAGvXrsWWLVtKxTMAS7v3Nzh//fq1xElIFd5fv/nhHY1KAk8EoiLLzc3F0KFDkZWVhS1btkjyD7kwBg4ciCNHjsDX1xctWrSAvb291JGomOjq6qJcuXJITU0F8Pft81Rxizgqebm5uXj8+DFMTExEPQ1J1ViaVGTLli1DSEgIDh48KOmlJYWxfPlynD59Gv369cP58+d5DZ8We3+LuPfFSZpLR0cHNjY2kvziw5sbFAFvbgDExcWhSZMmGDNmDJYsWSJ1nEK5cuUKnJ2dMX78eCxYsEDqOFTMcnJy8r1hOGkOAwMD6Ojk/+1icX8uszSLoLSXZk5ODlq1aoXnz58jNjYWxsbGUkcqtAULFmD69Ok4deoUmjdvLnUcIiqk4v5c5olAVGhLlizBuXPnEBAQoNGFCQATJ06Eo6Mjhg0bhrdv30odh4jUFEuTCiU+Ph4zZ86En5+fWl6P+bn09PSwadMm3L59G3PmzJE6DhGpKZYmfbb3d/2xtrb+6INwNU39+vXxww8/4Oeff8bVq1eljkNEaoilSZ9t27ZtiIiIwIoVKzT+sOyHJk6ciJo1a+Jf//oXL4QnojxYmvRZ0tPTMWHCBPTp0wddunSROo7KGRoaYuXKlTh27Bj++OMPqeMQkZphadJnmTVrFl68eIGlS5dKHaXYdOrUCV5eXpg4cSLS0tKkjkNEaoSlSaJdu3YNv/76K2bNmoVq1apJHadYvX/+5qxZs6SOQkRqhKVJor3/vm/s2LFSRyl2VatWxbRp07Bq1SrcvHlT6jhEpCZYmiRKWFgYQkNDsXDhQhgYGEgdp0SMGzcOVlZW+P7776WOQkRqgqVJn5STk4OJEyeidevW6Nmzp9RxSoyxsTHmz5+P4OBgHDt2TOo4RKQGWJr0SQEBAbhy5QoWL15c6p4MMWDAADg7O2PChAnIzc2VOg4RSYylSR/15s0bzJo1CwMHDkTTpk2ljlPidHR0sGTJEly4cAF//vmn1HGISGIsTfqoNWvW4NGjR6X61nKtW7dGt27d8MMPPyA7O1vqOEQkIZYmFejVq1dYsGABhg4dilq1akkdR1Lz5s3DzZs3ecMDolKOpUkFWrFiBZ4/f46ZM2dKHUVyjo6O6N27N+bMmYOsrCyp4xCRRFialK+MjAwsWrQII0eOhK2trdRx1MKcOXNw7949rF+/XuooRCQRlibla9myZXj16hWmTZsmdRS1Ub9+fQwaNAjz5s3jMzeJSimWJuXx/PlzLFmyBKNHj0bVqlWljqNWZsyYgZSUFPz+++9SRyEiCbA0KY/ly5cjKysLU6ZMkTqK2rGzs0OfPn2wcOFCnklLVAqxNEnJy5cvsXz5cowcORKVKlWSOo5amjZtGhISErBt2zapoxBRCWNpkpL169cjIyMDfn5+UkdRW40bN4a7uzv8/f15lyCiUoalSQpZWVlYvHgxBg4cyDNmP2H69OmIi4vDnj17pI5CRCWIpUkKW7Zswf379/lUDxFcXFzQvn17zJ8/H4IgSB2HiEoIS5MAALm5uVi4cCG6d++OevXqSR1HI0yfPh0XL17E4cOHpY5CRCWEpUkAgL179yI+Pp5nzH6GDh06oGnTpvj555+ljkJEJYSlSRAEAQsWLEC7du3QokULqeNoDJlMBj8/Pxw5cgSXL1+WOg4RlQCWJuHUqVM4d+4cJk2aJHUUjdO7d29YW1tj6dKlUkchohLA0iQsW7YMdnZ26Nq1q9RRNI6+vj6+++47bN26FSkpKVLHIaJixtIs5e7du4fdu3fj3//+N3R0+M+hMEaOHAl9fX2sWrVK6ihEVMz4KVnKrVy5EmZmZvDx8ZE6isYqV64chg8fjlWrVuHNmzdSxyGiYsTSLMVevXqFdevWYcSIEShTpozUcTTa2LFj8ezZMwQGBkodhYiKEUuzFNu8eTMyMjLg6+srdRSN98UXX6BHjx5YunQpb61HpMW0qjRXrVqFGjVqwMjICE5OTjh+/Lio+U6ePAk9PT00bty4eAOqkdzcXPzyyy/o2bMnb5mnIuPGjcONGzdw9OhRqaMQUTHRmtLcvn07xo0bh+nTpyMmJgatW7dGt27dkJiY+NH50tPT4ePjg44dO5ZQUvUQFhaG+Ph4jB07VuooWqN169ZwcHDgCUFEWkwmaMmNM5s3bw5HR0esXr1aMcze3h49evSAv79/gfP1798ftWvXhq6uLoKDgxEbGyt6nRkZGZDL5UhPT4eZmVlR4pc4Dw8PPHjwABcvXoRMJpM6jtZYvXo1fH19ce/ePVSrVk3qOESlTnF/LmvFnmZWVhaio6Ph6uqqNNzV1RWnTp0qcL6AgADcvn0bs2bNErWezMxMZGRkKL000b179xASEoIxY8awMFVs8ODBMDExwW+//SZ1FCIqBlpRmk+ePEFOTk6ehyZXqlSpwAvOb926hSlTpmDLli3Q09MTtR5/f3/I5XLFy9rausjZpbBu3TqULVsWAwYMkDqK1ilbtix8fHywbt06ZGVlSR2HiFRMK0rzvQ/3mgRByHdPKicnBwMHDsScOXNQp04d0cufOnUq0tPTFa+kpKQiZy5pWVlZWL9+Pby9vXmZSTEZPXo0UlJSEBwcLHUUIlIxcbtYas7CwgK6urp59ipTU1Pz7H0CwIsXL3DhwgXExMQoLrfIzc2FIAjQ09NDWFgYOnTokGc+Q0NDGBoaFs+bKCF79+7Fo0ePMGrUKKmjaC0HBwe0adMGq1atQr9+/aSOQ0QqpBV7mgYGBnByckJ4eLjS8PDwcLRs2TLP9GZmZrhy5QpiY2MVr1GjRsHOzg6xsbFo3rx5SUUvcWvWrMGXX34JBwcHqaNotTFjxiAqKgrXrl2TOgoRqZBW7GkCgJ+fH7y9veHs7AwXFxf89ttvSExMVOxRTZ06FQ8ePMDmzZuho6OTpzQsLS1hZGSk1WUSHx+Po0eP4o8//pA6itbr2bMnKlWqhNWrV2PFihVSxyEiFdGa0vTy8sLTp08xd+5cJCcnw8HBASEhIYoL95OTkz95zaa2W7t2LczNzdG7d2+po2g9AwMDxf1o//Of/8DY2FjqSESkAlpznaYUNOk6zTdv3qBq1aoYMWIEFi1aJHWcUuGvv/5C7dq1ERgYiMGDB0sdh6hU4HWapBJ79uzB8+fPMXLkSKmjlBq1atVCu3btsGHDBqmjEJGKsDRLiQ0bNqBNmzaoXbu21FFKlREjRiAyMhK3bt2SOgoRqQBLsxRISEjA0aNHMXz4cKmjlDq9evVCuXLlsHHjRqmjEJEKsDRLgYCAAJQtWxZ9+vSROkqpY2xsjMGDB2PTpk3Izs6WOg4RFRFLU8vl5ORg06ZN6N+/P0xNTaWOUyp9/fXXSElJQUhIiNRRiKiIWJpa7siRI0hKSuKhWQk1btwYTk5OWL9+vdRRiKiIWJpabsOGDahXr55W3+VIE4wYMQIhISF4+PCh1FGIqAhYmlrs6dOnCA4OxvDhw/kIMIkNGDAA+vr62Lp1q9RRiKgIWJpabOvWrcjNzYW3t7fUUUo9uVyO7t274/fffwfvJ0KkuViaWiwgIAAeHh6wtLSUOgoB8PHxwdWrV3Hp0iWpoxBRIbE0tdS1a9cQExMDHx8fqaPQ/7i6usLS0hKbN2+WOgoRFRJLU0sFBgaifPnycHNzkzoK/Y+enh4GDRqELVu28JpNIg3F0tRCubm52LJlC7y8vDT+odnaxtvbG6mpqQgLC5M6ChEVAktTC0VFReH+/ft8soYaaty4MRwcHHiIlkhDsTS1UGBgIGrWrImWLVtKHYU+IJPJ4OPjg+DgYKSlpUkdh4g+E0tTy7x+/Ro7d+7E4MGDeW2mmho0aBDevXuHnTt3Sh2FiD4TS1PL7Nu3Dy9evOChWTVmZWWFTp068RAtkQZiaWqZwMBAtGjRgs/NVHM+Pj44fvw47ty5I3UUIvoMLE0tkpqaisOHD3MvUwP06NEDZcqUQWBgoNRRiOgzsDS1SFBQEGQyGby8vKSOQp9gamqK3r17Y+vWrbytHpEGYWlqkcDAQLi5ucHCwkLqKCTCgAEDcPPmTcTExEgdhYhEYmlqifj4eFy4cIE3Z9cgHTp0gIWFBbZt2yZ1FCISiaWpJYKCglC2bFl4eHhIHYVE0tfXR9++fREUFITc3Fyp4xCRCCxNLSAIAoKCgtCjRw8YGRlJHYc+w4ABA3D//n2cOnVK6ihEJAJLUwtcvnwZN27cQP/+/aWOQp+pVatWqFatGg/REmkIlqYWCAoKQoUKFdCpUyepo9Bn0tHRgZeXF/78808++YRIA7A0NZwgCNi+fTt69+4NAwMDqeNQIQwYMACPHz/G0aNHpY5CRJ/A0tRw58+fR0JCAg/NajBHR0fUrl2bh2iJNABLU8MFBQWhUqVKaNu2rdRRqJBkMhkGDBiA3bt34+3bt1LHIaKPYGlqsNzcXOzYsQN9+/aFrq6u1HGoCPr374+MjAwcOnRI6ihE9BEsTQ128uRJPHjwgIdmtYC9vT0aNWqEoKAgqaMQ0UewNDVYUFAQrK2t4eLiInUUUoEBAwZg//79ePHihdRRiKgALE0NlZ2djT///BNeXl7Q0eFfozbw8vLCmzdvcPDgQamjEFEB+GmroSIjI/H48WMemtUi1atXR9OmTfHnn39KHYWICsDS1FBBQUH44osv4OjoKHUUUqE+ffogJCQEL1++lDoKEeWDpamB3r17hz179qBfv36QyWRSxyEV6tOnD96+fYuQkBCpoxBRPliaGigqKgrPnj1Dnz59pI5CKlazZk04OjryEC2RmmJpaqBdu3ahevXqaNKkidRRqBj07dsXISEhePXqldRRiOgDLE0Nk5OTg927d6NPnz48NKul+vTpg9evX/NGB0RqiKWpYU6ePInU1FT07t1b6ihUTGrVqoXGjRtj586dUkchog+wNDXMzp07Ua1aNTRr1kzqKFSM+vTpgwMHDuDNmzdSRyGif2BpapDc3Fzs3r0bvXr14g0NtFzfvn3x6tUrHqIlUjP85NUgZ8+exYMHD3jWbClQp04dNGjQgIdoidQMS1OD7Nq1C5UqVULLli2ljkIloG/fvti/fz8P0RKpEZamhhAEATt37kSvXr34GLBSok+fPnj58iXCwsKkjkJE/8PS1BAXL17EvXv3eNZsKWJvb4/69evzRgdEaoSlqSF27twJc3NztG3bVuooVIJ69eqFAwcOICsrS+ooRASWpkZ4f2i2R48e0NPTkzoOlaCePXsiPT0dkZGRUkchIrA0NcKVK1fw119/8dBsKdS4cWPY2tpiz549UkchIrA0NcLOnTshl8vRsWNHqaNQCZPJZOjZsyf27t2L3NxcqeMQlXosTQ2wa9cueHp6wsDAQOooJIGePXsiOTkZZ8+elToKUanH0lRz8fHxuH79Og/NlmKtWrVCxYoVeYiWSA2wNNVccHAwTExM4OrqKnUUkoiuri48PT2xZ88eCIIgdRyiUo2lqeaCg4PRpUsXGBsbSx2FJNSzZ0/89ddfuHbtmtRRiEo1lqYaS05OxpkzZ9CjRw+po5DEOnbsiDJlyvAQLZHEWJpqbP/+/dDV1YW7u7vUUUhiRkZGcHNzQ3BwsNRRiEo1lqYaCw4ORps2bWBubi51FFIDPXv2VNxOkYikwdJUUxkZGThy5Ai6d+8udRRSE25ubjAwMODeJpGEWJpqKjQ0FFlZWSxNUjAzM0PHjh35vSaRhFiaaio4OBiNGzdG9erVpY5CaqRnz544fvw4Hj9+LHUUolKJpamGsrKycPDgQZ41S3l4enpCEATs379f6ihEpZJWleaqVatQo0YNGBkZwcnJCcePHy9w2t27d6Nz586oWLEizMzM4OLigsOHD5dg2oJFRkYiIyODpUl5VKpUCa1ateIhWiKJaE1pbt++HePGjcP06dMRExOD1q1bo1u3bkhMTMx3+mPHjqFz584ICQlBdHQ02rdvj6+++goxMTElnDyvvXv3onr16mjYsKHUUUgN9ezZE+Hh4Xj58qXUUYhKHZmgJfflat68ORwdHbF69WrFMHt7e/To0QP+/v6illG/fn14eXnhhx9+yHd8ZmYmMjMzFT9nZGTA2toa6enpMDMzK9ob+J/c3FzY2Nigb9++WLp0qUqWSdrl9u3bqFWrFnbt2oVevXpJHYdIrWRkZEAul6v0c/mftGJPMysrC9HR0Xnuz+rq6opTp06JWkZubi5evHiBChUqFDiNv78/5HK54mVtbV2k3PmJjo7GgwcPeNYsFeiLL75AvXr1+L0mkQS0ojSfPHmCnJwcVKpUSWl4pUqVkJKSImoZixcvxqtXr9CvX78Cp5k6dSrS09MVr6SkpCLlzk9wcDAqVKiAL7/8UuXLJu3h6emJAwcOICcnR+ooRKWKVpTmezKZTOlnQRDyDMvPtm3bMHv2bGzfvh2WlpYFTmdoaAgzMzOll6oFBwfjq6++gp6ensqXTdrD09MTT548wZkzZ6SOQlSqaEVpWlhYQFdXN89eZWpqap69zw9t374dX3/9NXbs2IFOnToVZ8xPunXrFq5fv86zZumTmjVrBktLS+zbt0/qKESlilaUpoGBAZycnBAeHq40PDw8HC1btixwvm3btmHo0KHYunWrWtwUfe/evTAyMkLnzp2ljkJq7v2N/FmaRCVLK0oTAPz8/LB+/Xps3LgRcXFxGD9+PBITEzFq1CgAf38f6ePjo5h+27Zt8PHxweLFi9GiRQukpKQgJSUF6enpUr0F7N+/H506dYKpqalkGUhzeHp64saNG7h165bUUYhKDa0pTS8vLyxbtgxz585F48aNcezYMYSEhMDW1hbA38+m/Oc1m2vXrkV2djb+9a9/oUqVKorX2LFjJcn/7NkznDx5El999ZUk6yfN07lzZxgaGvIsWqISpDXXaUpBldcDbdmyBYMHD8aDBw9gZWWlooSk7Tw8PPDy5UtERkZKHYVILfA6zVJi//79cHJyYmHSZ/H09MSJEyfw9OlTqaMQlQosTTXw7t07hIaG8tAsfTYPDw/k5OTg0KFDUkchKhVYmmrgxIkTSE9Ph4eHh9RRSMNYWVnB2dmZ32sSlRCWphrYv38/rKys4OjoKHUU0kCenp44dOgQsrKypI5CpPVYmhJ7/2xEDw8PUXcvIvqQp6cnXrx4gaioKKmjEGk9lqbE4uPj8ddff/H7TCq0hg0bwsbGhjc6ICoBLE2J7d+/H8bGxujYsaPUUUhDyWQyfPXVV9i3bx94BRlR8WJpSuz9XYCMjY2ljkIazNPTE4mJibhy5YrUUYi0GktTQk+fPuVdgEgl2rZti7Jly/IQLVExY2lKKDQ0FLm5uWpxs3jSbIaGhujatStLk6iYsTQlxLsAkSp5enri/PnzePjwodRRiLQWS1MivAsQqZqbmxt0dHRw8OBBqaMQaS2WpkSOHz+O9PR0liapTIUKFdCyZUuWJlEx0ivMTH5+fqKnXbJkSWFWofX279+PqlWrokmTJlJHIS3i7u6OH3/8EW/fvoWRkZHUcYi0TqFKMyYmRunn6Oho5OTkwM7ODgBw8+ZN6OrqwsnJqegJtRDvAkTFxd3dHVOnTkVUVBS6dOkidRwirVOo0oyIiFD8ecmSJShbtix+//13lC9fHgDw/PlzDBs2DK1bt1ZNSi0THx+P27dv45dffpE6CmkZBwcH2NjY4ODBgyxNomJQ5O80Fy9eDH9/f0VhAkD58uUxb948LF68uKiL10rv7wLUoUMHqaOQlpHJZHB3d8eBAwd4dyCiYlDk0szIyMCjR4/yDE9NTcWLFy+KunitxLsAUXFyd3dHQkICbty4IXUUIq1T5NLs2bMnhg0bhp07d+L+/fu4f/8+du7cia+//hq9evVSRUatwrsAUXFr3749jIyMeBYtUTEocmmuWbMG7u7uGDx4MGxtbWFra4tBgwahW7duWLVqlSoyapVDhw4hNzeXD5ymYmNiYoKOHTuyNImKQaFLc9q0aTh37hxMTEywatUqPH36FDExMbh48SKePXuGVatWwdTUVJVZtcKBAwfg5OSEKlWqSB2FtJi7uztOnDiBtLQ0qaMQaZVCl2ZycjI8PDxQpUoVfPPNN4iKioKdnR0aNWrEsixAdnY2Dh8+zL1MKnbu7u7Izs5GWFiY1FGItEqhSzMgIACPHj3Cjh07UK5cOfj5+cHCwgK9evXCpk2b8OTJE1Xm1AqnT59GWloa3NzcpI5CWs7GxgYODg48REukYkX6TlMmk6F169ZYtGgRbty4gXPnzqFFixZYt24dqlatijZt2uDnn3/GgwcPVJVXo4WEhKBixYpwdnaWOgqVAh4eHorv0IlINVR671l7e3tMnjwZJ0+exP379zFkyBAcP34c27ZtU+VqNFZISAi6desGHR3e8peKn7u7Ox4/fozz589LHYVIa8gEXgFdaBkZGZDL5UhPT4eZmdlHp01KSoKNjQ2CgoLg5eVVQgmpNMvOzoalpSV8fX0xd+5cqeMQlYjP+VwujM/e5Xn+/DmePXsGAHj8+DF27dqFq1evqjyYtjl06BB0dXXh6uoqdRQqJfT09NC1a1ccOHBA6ihEWuOzSnP9+vVwdnaGk5MTVq9ejZ49e+LIkSPo378/fvvtt+LKqBVCQkLQsmVLpdsNEhU3Dw8PxMTE8MHURCryWTds//XXX3Ht2jW8fv0aNjY2SEhIQMWKFZGRkYE2bdrgm2++Ka6cGi0zMxP//e9/MWPGDKmjUCnTtWtX6OjoICQkBCNGjJA6DpHG+6w9TV1dXRgZGaFChQqoVasWKlasCAAwMzPjI64+4tixY3j16hXc3d2ljkKlTIUKFeDi4sJLT4hU5LNKU09PD2/fvgUAREVFKYbzxuwfFxISgmrVqsHBwUHqKFQKeXh4IDw8HJmZmVJHIdJ4n1WaR48ehaGhIQBALpcrhr958wYbNmxQbTItcvDgQbi5uXFvnCTh7u6OV69eKf2iS0SF81mlWaZMmXw/+C0tLeHo6KiyUNrk1q1buHXrFg/NkmQcHBxgbW3NQ7REKvBZJwIV5O3bt7h8+TJSU1Pz3H3E09NTFavQWIcOHYKBgQEfOE2S+eeDqZctW8YjHkRFUOTSDA0NhY+PT773mpXJZMjJySnqKjTawYMH0bZtW5QpU0bqKFSKeXh4YM2aNYiPj0fdunWljkOksYp8PzdfX1/07dsXycnJyM3NVXqV9sJ89eoVIiMjeWiWJMcHUxOpRpFLMzU1FX5+fqhUqZIq8miVo0ePIisri081IcmZmJigQ4cOLE2iIipyafbp0weRkZEqiKJ9Dh48iFq1aqF27dpSRyGCu7s7jh8/jvT0dKmjEGmsIt+w/fXr1+jbty8qVqyIBg0aQF9fX2n8v//97yIFVGcfuzGwIAiwtbVFr169sGzZMmkCEv3DvXv3UL16dezYsQN9+/aVOg5RsSjuG7YX+USgrVu34vDhwzA2NkZkZKTSmXkymUyrS/Njrl27hqSkJB6aJbVha2sLBwcHHDp0iKVJVEhFLs0ZM2Zg7ty5mDJlCp8T+Q8HDx6EiYkJ2rRpI3UUIgU3Nzds3rwZubm5/P9KVAhF/l+TlZUFLy8v/gf8QEhICDp16gQjIyOpoxApuLm5ISUlBbGxsVJHIdJIRW66IUOGYPv27arIojWeP3+OkydP8tAsqZ2WLVvCzMwMISEhUkch0khFPjybk5ODRYsW4fDhw2jYsGGeE4GWLFlS1FVonPDwcOTk5LA0Se3o6+ujc+fOCAkJ4aPqiAqhyKV55coVNGnSBABw9epVpXGl9XZdISEhaNCgAaytraWOQpSHm5sbRowYgSdPnsDCwkLqOEQapcilGRERoYocWiM3NxeHDh3CsGHDpI5ClK+uXbtCEASEhYVh4MCBUsch0ig8e0fFoqOjkZqaylvnkdqysrJCkyZN+L0mUSGwNFUsJCQEcrkcLi4uUkchKpCbmxtCQ0NL/f2hiT4XS1PFDh48iC5dukBPTyVPXSMqFm5ubnj69CnOnz8vdRQijcLSVKFHjx7h/PnzPDRLaq958+aoUKECD9ESfSaWpgodPnwYwN8nWhCpM11dXXTp0oWlSfSZirU0dXR00KFDB0RHRxfnatTGwYMH0bRpU1haWkodheiT3NzcEB0djZSUFKmjEGmMYi3NjRs3om3btqXipu3Z2dk4fPgwD82SxujSpQtkMhlCQ0OljkKkMYr8aLDS7J+PoLl06RLatGmDc+fOoWnTplJHIxKlRYsWsLGxwY4dO6SOQqQSxf1oMH6nqSIHDx6EpaUlnJycpI5CJJqbmxvCwsLw7t07qaMQaQSVlObx48cxePBguLi44MGDBwCAwMBAnDhxQhWL1wghISHo1q0bn/ZCGqVbt25IT0/H6dOnpY5CpBGK/Am/a9cudOnSBcbGxoiJiUFmZiYA4MWLF5g/f36RA2qC+/fv48qVK7xBO2kcJycnVKxYkWfREolU5NKcN28e1qxZg3Xr1ik94aRly5a4ePFiURevEcLCwqCrqwtXV1epoxB9Fh0dHXTr1o2lSSRSkUszPj4ebdq0yTPczMwMaWlpRV28RggLC0OrVq1Qrlw5qaMQfTY3NzdcuXIFSUlJUkchUntFLs0qVargr7/+yjP8xIkTqFmzZlEXrxEiIyN5aJY0lqurK3R0dHDo0CGpoxCpvSKX5rfffouxY8fi7NmzkMlkePjwIbZs2YKJEydizJgxqsgo2qpVq1CjRg0YGRnByckJx48f/+j0UVFRcHJygpGREWrWrIk1a9YUar1v3rzh9ZmkscqXL4+WLVvyEC2RGIIKTJs2TTA2NhZkMpkgk8kEIyMjYcaMGapYtGhBQUGCvr6+sG7dOuH69evC2LFjBVNTU+HevXv5Tn/nzh3BxMREGDt2rHD9+nVh3bp1gr6+vrBz507R60xPTxcACFWrVhVyc3NV9VaIStz8+fMFU1NT4e3bt1JHISqS95/L6enpxbL8It/cIDExEdWqVcPbt29x/fp15Obmol69ejA1NUVSUhJsbGxU0+6f0Lx5czg6OmL16tWKYfb29ujRowf8/f3zTP/9999j3759iIuLUwwbNWoULl26JPr0+/cX0Q4bNgwbN24s+psgksilS5fQuHFjhIeHo1OnTlLHISo0tb+5QY0aNfDkyROYmJjA2dkZzZo1Q5kyZfDs2TPUqFFDFRk/KSsrC9HR0XnOXnV1dcWpU6fynef06dN5pu/SpQsuXLhQ4IXemZmZyMjIUHq9n49IkzVs2BBVq1bF9u3bpY5CpNaK/NDHgnZUX758CSMjo6IuXpQnT54gJycHlSpVUhpeqVKlAm9GnZKSku/02dnZePLkCapUqZJnHn9/f8yZMyfP8ICAAGzZsgUAYGxsjO3bt2PevHkYPnw4Ro0apTTt9OnTERMTgz59+mDhwoWIj49XjGvVqhW6d++OK1euwMTEBGvXrlWad9euXVi4cCH8/PzQv39/pXG+vr54+vQpmjVrhm3btuHcuXOKcQ4ODvD19cWhQ4dgZ2eHRYsWKc27ceNGrFmzBjNmzECvXr2QnZ2tGOft7Q0TExNUqVIFJ0+eRHh4uGKctbU1/vOf/2D16tVwdXXF9OnTlZa7dOlSbNu2DTNmzMCwYcPw9OlTxbju3bvDzs4O2dnZSExMxM6dOxXjzMzM8Mcff2DevHnw8fGBr6+v0nJnz56NkydPwsfHB7NmzcKdO3cU49q2bQtXV1fcvn1b8d7+ad++fZg3bx6+++47eHt7K43z8/NDYmIi2rVrhw0bNiAmJkYxrkmTJvj6668RGRkJGxsbLFmyRGnewMBA/Prrr5gxYwY8PT2Vxg0fPhwA8MUXXyAsLAxRUVGKcTVr1sScOXOwefNmtGrVCrNnz1aad8WKFdi8eTNmzJiBwYMHK35RA4A+ffrAxsYGenp6iI+Px969exXjzM3NERAQgHnz5mHAgAEYP3680nJ/+uknhIWFYfTo0Zg0aRKSkpIgl8uxfv16WFhYoG/fvkhOTsbr168RGBiomE9PTw+7d+/GvHnzMGrUKMV7e2/y5MmIj49Ht27dsGLFCly9elUxrlmzZhgwYADOnTsHc3NzrFixQmneoKAgLFmyBN9//z169+6tNO7bb7/F69ev0aBBA+zduxcnT55UjLOzs8P333+PnTt3okmTJvjpp5+U5l2zZg02btyIGTNmwMvLC2/evFGMGzBgAMzNzSGXyxETE6P0vW6VKlWwYsUKLFmyBN27d8fkyZOVlrto0SLs3bsXfn5+8PX1RXJysmKcm5sbmjRpgvT0dDx9+hTbtm1TjONnxP8rjs+IV69eoTgV+vCsn58fAOCXX37ByJEjYWJiohiXk5ODs2fPQldXV+kfd3F5+PAhqlatilOnTsHFxUUx/KeffkJgYCBu3LiRZ546depg2LBhmDp1qmLYyZMn8eWXXyI5ORmVK1fOM09mZqbi5g3A34cBrK2ti+0wAFFJysrKQsOGDWFlZYUjR45AJpNJHYnosw0fPhwBAQHF9rlc6D3N97+FC4KAK1euwMDAQDHOwMAAjRo1wsSJE4ueUAQLCwvo6urm2atMTU3Nszf5XuXKlfOdXk9PD+bm5vnOY2hoCENDQ9WEJlIzBgYGWLp0Kdzc3LB79+48e3tE6i42NhYBAQHFuo5Cl2ZERAQAYNiwYfjll18k3dMyMDCAk5MTwsPD0bNnT8Xw8PBwdO/ePd95XFxcsH//fqVhYWFhcHZ2VrqzEVFp0q1bN7i7u2PChAlwc3ODsbGx1JGIRBEEAePGjUPt2rVx69atYl2RVnh/ycmGDRuE69evC+PGjRNMTU2Fu3fvCoIgCFOmTBG8vb0V07+/5GT8+PHC9evXhQ0bNhT6kpPiOrWZSArx8fGCvr6+8OOPP0odhUi0P//8UwAg7Ny5s1g/l4t8ItB7169fR2JiIrKyspSGf3hSRHHx8vLC06dPMXfuXCQnJ8PBwQEhISGwtbUFACQnJyMxMVExfY0aNRASEoLx48dj5cqVsLKywvLly3lIikq9OnXqYOzYsfD398fQoUNRrVo1qSMRfdSbN28wadIkuLu7o3PnzsW6riJfp3nnzh307NkTV65cgUwmU5xN+/4kgpycnKKnVFPFfT0QkVQyMjJQp04ddOjQAVu3bpU6DtFH/fTTT5g9ezauXbuGypUrq/d1mmPHjkWNGjXw6NEjmJiY4Nq1azh27BicnZ0RGRmpgohEVNLMzMzg7++Pbdu2larn4pLmuX//PubPn49///vfqFOnTrGvr8h7mhYWFjh69CgaNmwIuVyOc+fOwc7ODkePHsWECROUrnXTNtzTJG2Wm5uLFi1aIDs7G+fPn4eurq7UkYjyGDRoEP773//i5s2bkMvl6n9HoJycHJQpUwbA3wX68OFDAICtra3SRblEpFl0dHSwfPlyxMTE8DaRpJZOnjyJrVu3wt/fH3K5vETWWeTSdHBwwOXLlwH8ff/XRYsW4eTJk5g7d26peTQYkbZq0aIFvL29MW3atFLzfFzSDDk5Ofjuu+/g7OyMoUOHlth6i1yaM2bMQG5uLgBg3rx5uHfvHlq3bo2QkBAsX768yAGJSFoLFizAmzdv8r2FJJFUNm7ciJiYGCxfvhw6OkWuMtFU9pSTf4Z+9uwZypcvX6JPOZECv9Ok0sLf3x8//PADLl++DHt7e6njUCn3/Plz1KlTB25ubvj999+Vxqn9d5rvn3LyTxUqVCjRp5wQUfEaP348bGxsMH78+AIf0kBUUubMmYO3b99iwYIFJb7uIpdmQf+BSvIpJ0RUvIyMjLBkyRIcPnwYBw4ckDoOlWLXrl3DihUrMHPmzHyfRlXctOIpJ1Lh4VkqTQRBQJcuXXDnzh1cu3aNDy+gEicIAjp37ox79+7h6tWr+f4bLO7PZa14ygkRFT+ZTIZly5ahYcOGWLZsGb7//nupI1EpExwcjCNHjuDAgQOS/dJW5BOB1OEpJ1LhniaVRuPGjcOGDRtw8+ZNSQ6PUen05s0b1KtXD/b29koPC/+Q2p8IFBAQwMIgKkVmzZoFQ0NDTJkyReooVIosXrwYDx48wNKlSyXNUXIXtxCRVihfvjx++uknbN68GWfPnpU6DpUCSUlJ8Pf3x9ixY2FnZydpliIfni3NeHiWSqucnBw4OzvDwMAAp0+fLtGLy6n0GTBgACIiInDz5s1Pftaq/eFZIip9dHV1sXz5cpw7dw6BgYFSxyEtdvz4cQQFBWHBggVqsXPCPc0i4J4mlXb9+/dHVFQU4uPj+X+AVC47OxvOzs4wNDQUfURD7fc0Bw0ahN9++w03b95URR4i0iD/+c9/kJGRgdmzZ0sdhbTQmjVrcPnyZfz6669q8xVAkVOUKVMGS5YsQd26dWFlZYUBAwZgzZo1uHHjhiryEZEas7a2xg8//IDly5fjypUrUschLfLo0SPMmDEDI0eORLNmzaSOo6Cyw7MpKSmIjIxEZGQkoqKicPPmTVhaWiI5OVkVi1dLPDxLBGRlZaFRo0awsLDAsWPHIJPJpI5EWmDIkCE4ePAg4uPjYW5uLno+tT88+17ZsmVRvnx5lC9fHuXKlYOenh4qV66sqsUTkZoyMDDAihUrcOLECZ4URCpx/PhxbN68GQsWLPiswiwJRd7T/P777xEVFYVLly7BwcEBbdq0Qdu2bdGmTRuUK1dORTHVE/c0if5f//79ERERgfj4eK3/v0/F5927d3B0dISpqSlOnTr12d9lFvfncpFLU0dHBxUrVsT48ePRvXv3UvWsPZYm0f978OAB6tati6FDh+LXX3+VOg5pqKVLl2LChAm4cOECHB0dP3t+tT88GxMTg+nTp+PcuXNo06YNKleuDC8vL6xevRpxcXGqyEhEGqBq1aqYPXs2Vq1apXigA9HnePjwIWbNmoXRo0cXqjBLgsqv07x06RKWLVuGP/74A7m5ucjJyVHl4tUK9zSJlL179w5NmjSBmZkZTpw4oTaXCZBmGDRoEMLDwxEfH4/y5csXahlq+2iwf4qJiVGcOXv8+HFkZGSgcePGaN++vSoWT0QaQl9fHytXrkS7du2wadMmDB8+XOpIpCEiIiKwdetWBAQEFLowS0KR9zTLly+Ply9folGjRmjXrh3atWuHNm3alIo9L+5pEuVv8ODBOHz4MOLj41GhQgWp45Cae/fuHRo1aoTy5cvj+PHjRTpCofZ7moGBgaWmJIlInJ9//hn79+/H9OnTsXr1aqnjkJpbtmwZ4uPjcfHiRbU/pF/k0vTw8FBFDiLSIpUrV8bcuXMxfvx4WFtbo06dOlJHIjWVkpKC2bNnw9fXF40aNZI6ziep5DvNtLQ0bNiwAXFxcZDJZLC3t8fXX38NuVyuisUTkQb617/+hYSEBMyYMQN8LgR9TL9+/fDTTz9JHUOUIn+neeHCBXTp0gXGxsZo1qwZBEHAhQsX8ObNG4SFhantacOqwO80iT7t1atXyM7OljoGqSkdHR2ULVtWZctT+5sbtG7dGrVq1cK6deugp/f3jmt2djZGjBiBO3fu4NixYyoJqo5YmkRE6kXtS9PY2BgxMTGoW7eu0vDr16/D2dkZr1+/LlJAdcbSJCJSL2p/RyAzMzMkJibmGZ6UlKTSXW4iIiKpFbk0vby88PXXX2P79u1ISkrC/fv3ERQUhBEjRmDAgAGqyEhERKQWinz27M8//wyZTAYfHx/Fl/36+voYPXo0FixYUOSARERE6kJl9559/fo1bt++DUEQUKtWLZiYmKhisWqN32kSEakXtb8j0HsmJiZo0KCBqhZHRESkdgpVmn5+fqKnXbJkSWFWQUREpHYKVZpin5Unk8kKs3giIiK1VKjSjIiIwJ07d1C9enW1v7kuERGRqhS68WrXro0nT54ofvby8sKjR49UEoqIiEgdFbo0PzzpNiQkBK9evSpyICIiInXFY6tEREQiFbo0ZTJZnhN9eOIPERFps0JfpykIAoYOHQpDQ0MAwNu3bzFq1CiYmpoqTbd79+6iJSQiIlIThS7NIUOGKP08ePDgIochIiJSZ4UuzYCAAFXmICIiUns8EYiIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRNKK0nz+/Dm8vb0hl8shl8vh7e2NtLS0Aqd/9+4dvv/+ezRo0ACmpqawsrKCj48PHj58WHKhiYhI42hFaQ4cOBCxsbEIDQ1FaGgoYmNj4e3tXeD0r1+/xsWLFzFz5kxcvHgRu3fvxs2bN+Hp6VmCqYmISNPIBEEQpA5RFHFxcahXrx7OnDmD5s2bAwDOnDkDFxcX3LhxA3Z2dqKWc/78eTRr1gz37t2DjY2NqHkyMjIgl8uRnp4OMzOzQr8HIiJSjeL+XNb4Pc3Tp09DLpcrChMAWrRoAblcjlOnToleTnp6OmQyGcqVK1fgNJmZmcjIyFB6ERFR6aHxpZmSkgJLS8s8wy0tLZGSkiJqGW/fvsWUKVMwcODAj/5m4u/vr/jeVC6Xw9rautC5iYhI86htac6ePRsymeyjrwsXLgAAZDJZnvkFQch3+IfevXuH/v37Izc3F6tWrfrotFOnTkV6errilZSUVLg3R0REGklP6gAF8fX1Rf/+/T86TfXq1XH58mU8evQoz7jHjx+jUqVKH53/3bt36NevHxISEnD06NFPHv82NDSEoaHhp8MTEZFWUtvStLCwgIWFxSenc3FxQXp6Os6dO4dmzZoBAM6ePYv09HS0bNmywPneF+atW7cQEREBc3NzlWUnIiLtpLaHZ8Wyt7dH165dMXLkSJw5cwZnzpzByJEj4eHhoXTmbN26dbFnzx4AQHZ2Nvr06YMLFy5gy5YtyMnJQUpKClJSUpCVlSXVWyEiIjWn8aUJAFu2bEGDBg3g6uoKV1dXNGzYEIGBgUrTxMfHIz09HQBw//597Nu3D/fv30fjxo1RpUoVxetzzrglIqLSReOv05QSr9MkIlIvvE6TiIhITbA0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCSSVpTm8+fP4e3tDblcDrlcDm9vb6SlpYme/9tvv4VMJsOyZcuKLSMREWk+rSjNgQMHIjY2FqGhoQgNDUVsbCy8vb1FzRscHIyzZ8/CysqqmFMSEZGm05M6QFHFxcUhNDQUZ86cQfPmzQEA69atg4uLC+Lj42FnZ1fgvA8ePICvry8OHz4Md3f3T64rMzMTmZmZip8zMjKK/gaIiEhjaPye5unTpyGXyxWFCQAtWrSAXC7HqVOnCpwvNzcX3t7emDRpEurXry9qXf7+/opDwHK5HNbW1kXOT0REmkPjSzMlJQWWlpZ5hltaWiIlJaXA+RYuXAg9PT38+9//Fr2uqVOnIj09XfFKSkoqVGYiItJMaluas2fPhkwm++jrwoULAACZTJZnfkEQ8h0OANHR0fjll1+wadOmAqfJj6GhIczMzJReRERUeqjtd5q+vr7o37//R6epXr06Ll++jEePHuUZ9/jxY1SqVCnf+Y4fP47U1FTY2NgohuXk5GDChAlYtmwZ7t69W6TsRESkndS2NC0sLGBhYfHJ6VxcXJCeno5z586hWbNmAICzZ88iPT0dLVu2zHceb29vdOrUSWlYly5d4O3tjWHDhhU9PBERaSW1LU2x7O3t0bVrV4wcORJr164FAHzzzTfw8PBQOnO2bt268Pf3R8+ePWFubg5zc3Ol5ejr66Ny5cofPduWiIhKN7X9TvNzbNmyBQ0aNICrqytcXV3RsGFDBAYGKk0THx+P9PR0iRISEZE2kAmCIEgdQlNlZGRALpcjPT2dJwUREamB4v5c1oo9TSIiopLA0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpEhERicTSJCIiEomlSUREJBJLk4iISCSWJhERkUgsTSIiIpFYmkRERCKxNImIiERiaRIREYnE0iQiIhKJpUlERCQSS5OIiEgkliYREZFILE0iIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJK0ozefPn8Pb2xtyuRxyuRze3t5IS0v75HxxcXHw9PSEXC5H2bJl0aJFCyQmJhZ/YCIi0khaUZoDBw5EbGwsQkNDERoaitjYWHh7e390ntu3b+PLL79E3bp1ERkZiUuXLmHmzJkwMjIqodRERKRpZIIgCFKHKIq4uDjUq1cPZ86cQfPmzQEAZ86cgYuLC27cuAE7O7t85+vfvz/09fURGBhY6HVnZGRALpcjPT0dZmZmhV4OERGpRnF/LuupfIkl7PTp05DL5YrCBIAWLVpALpfj1KlT+ZZmbm4uDh48iMmTJ6NLly6IiYlBjRo1MHXqVPTo0aPAdWVmZiIzM1Pxc3p6OoC//5KIiEh67z+Pi2t/UONLMyUlBZaWlnmGW1paIiUlJd95UlNT8fLlSyxYsADz5s3DwoULERoail69eiEiIgJt27bNdz5/f3/MmTMnz3Bra+uivQkiIlKpp0+fQi6Xq3y5aluas2fPzreg/un8+fMAAJlMlmecIAj5Dgf+3tMEgO7du2P8+PEAgMaNG+PUqVNYs2ZNgaU5depU+Pn5KX5OS0uDra0tEhMTi+UvR9UyMjJgbW2NpKQkjTmczMwlg5lLBjMXv/T0dNjY2KBChQrFsny1LU1fX1/079//o9NUr14dly9fxqNHj/KMe/z4MSpVqpTvfBYWFtDT00O9evWUhtvb2+PEiRMFrs/Q0BCGhoZ5hsvlco34x/SemZmZRuUFmLmkMHPJYObip6NTPOe5qm1pWlhYwMLC4pPTubi4ID09HefOnUOzZs0AAGfPnkV6ejpatmyZ7zwGBgZo2rQp4uPjlYbfvHkTtra2RQ9PRERaSeMvObG3t0fXrl0xcuRInDlzBmfOnMHIkSPh4eGhdBJQ3bp1sWfPHsXPkyZNwvbt27Fu3Tr89ddfWLFiBfbv348xY8ZI8TaIiEgDaHxpAsCWLVvQoEEDuLq6wtXVFQ0bNsxzKUl8fLzibFcA6NmzJ9asWYNFixahQYMGWL9+PXbt2oUvv/xS9HoNDQ0xa9asfA/ZqiNNywswc0lh5pLBzMWvuPNq/HWaREREJUUr9jSJiIhKAkuTiIhIJJYmERGRSCxNIiIikViahbRq1SrUqFEDRkZGcHJywvHjx6WOpDB79mzIZDKlV+XKlRXjBUHA7NmzYWVlBWNjY7Rr1w7Xrl0r0YzHjh3DV199BSsrK8hkMgQHByuNF5MxMzMT3333HSwsLGBqagpPT0/cv39fkrxDhw7Ns81btGghWV7g79s+Nm3aFGXLloWlpSV69OiR59pkddvOYjKr07ZevXo1GjZsqLjw38XFBYcOHVKMV7ftKyazOm3fgvj7+0Mmk2HcuHGKYSW2rQX6bEFBQYK+vr6wbt064fr168LYsWMFU1NT4d69e1JHEwRBEGbNmiXUr19fSE5OVrxSU1MV4xcsWCCULVtW2LVrl3DlyhXBy8tLqFKlipCRkVFiGUNCQoTp06cLu3btEgAIe/bsURovJuOoUaOEqlWrCuHh4cLFixeF9u3bC40aNRKys7NLPO+QIUOErl27Km3zp0+fKk1TknkFQRC6dOkiBAQECFevXhViY2MFd3d3wcbGRnj58qViGnXbzmIyq9O23rdvn3Dw4EEhPj5eiI+PF6ZNmybo6+sLV69eFQRB/bavmMzqtH3zc+7cOaF69epCw4YNhbFjxyqGl9S2ZmkWQrNmzYRRo0YpDatbt64wZcoUiRIpmzVrltCoUaN8x+Xm5gqVK1cWFixYoBj29u1bQS6XC2vWrCmhhMo+LCExGdPS0gR9fX0hKChIMc2DBw8EHR0dITQ0tETzCsLfHzTdu3cvcB4p876XmpoqABCioqIEQVD/7ZxfZkFQ/21dvnx5Yf369RqxfT/MLAjqvX1fvHgh1K5dWwgPDxfatm2rKM2S3NY8PPuZsrKyEB0dDVdXV6Xhrq6uOHXqlESp8rp16xasrKxQo0YN9O/fH3fu3AEAJCQkICUlRSm/oaEh2rZtqzb5xWSMjo7Gu3fvlKaxsrKCg4ODZO8jMjISlpaWqFOnDkaOHInU1FTFOHXI+/7mHu9vZK0J2/nDzO+p47bOyclBUFAQXr16BRcXF43Yvh9mfk8dty8A/Otf/4K7uzs6deqkNLwkt7Xa3ntWXT158gQ5OTl5bgZfqVKlAh9FVtKaN2+OzZs3o06dOnj06BHmzZuHli1b4tq1a4qM+eW/d++eFHHzEJMxJSUFBgYGKF++fJ5ppPh76NatG/r27QtbW1skJCRg5syZ6NChA6Kjo2FoaCh5XkEQ4Ofnhy+//BIODg4A1H8755cZUL9tfeXKFbi4uODt27coU6YM9uzZg3r16ik+iNVx+xaUGVC/7fteUFAQLl68qHi61T+V5L9llmYhffjYMeEjjyIrad26dVP8uUGDBnBxccEXX3yB33//XfGFvjrnf68wGaV6H15eXoo/Ozg4wNnZGba2tjh48CB69epV4HwlldfX1xeXL1/O9yk+6rqdC8qsbtvazs4OsbGxSEtLw65duzBkyBBERUUpxqvj9i0oc7169dRu+wJAUlISxo4di7CwMBgZGRU4XUlsax6e/UwWFhbQ1dXN85tJampqgY8ik5qpqSkaNGiAW7duKc6iVef8YjJWrlwZWVlZeP78eYHTSKlKlSqwtbXFrVu3AEib97vvvsO+ffsQERGBatWqKYar83YuKHN+pN7WBgYGqFWrFpydneHv749GjRrhl19+UevtW1Dm/Ei9fYG/D62mpqbCyckJenp60NPTQ1RUFJYvXw49PT3FektiW7M0P5OBgQGcnJwQHh6uNDw8PLzAR5FJLTMzE3FxcahSpQpq1KiBypUrK+XPyspCVFSU2uQXk9HJyQn6+vpK0yQnJ+Pq1atq8T6ePn2KpKQkVKlSBYA0eQVBgK+vL3bv3o2jR4+iRo0aSuPVcTt/KnN+1GFb/5MgCMjMzFTL7fupzPlRh+3bsWNHXLlyBbGxsYqXs7MzBg0ahNjYWNSsWbPktnUhTmAq9d5fcrJhwwbh+vXrwrhx4wRTU1Ph7t27UkcTBEEQJkyYIERGRgp37twRzpw5I3h4eAhly5ZV5FuwYIEgl8uF3bt3C1euXBEGDBhQ4pecvHjxQoiJiRFiYmIEAMKSJUuEmJgYxWU7YjKOGjVKqFatmvDf//5XuHjxotChQ4diO+39Y3lfvHghTJgwQTh16pSQkJAgRERECC4uLkLVqlUlyysIgjB69GhBLpcLkZGRSpcPvH79WjGNum3nT2VWt209depU4dixY0JCQoJw+fJlYdq0aYKOjo4QFhYmCIL6bd9PZVa37fsx/zx7VhBKbluzNAtp5cqVgq2trWBgYCA4OjoqnRIvtffXJ+nr6wtWVlZCr169hGvXrinG5+bmCrNmzRIqV64sGBoaCm3atBGuXLlSohkjIiIEAHleQ4YMEZ3xzZs3gq+vr1ChQgXB2NhY8PDwEBITE0s87+vXrwVXV1ehYsWKgr6+vmBjYyMMGTIkT5aSzCsIQr55AQgBAQGKadRtO38qs7pt6+HDhys+BypWrCh07NhRUZiCoH7b91OZ1W37fsyHpVlS25qPBiMiIhKJ32kSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSCxNIiIikViaREREIrE0iYiIRGJpUqknk8kQHBxc7OupXr06li1bpjbL0bR1f0pkZCRkMhnS0tKkjkJajKVJWi01NRXffvstbGxsYGhoiMqVK6NLly44ffq0Yprk5GSlx6mpi02bNqFcuXJ5hp8/fx7ffPON1q6bSJ3xeZqk1Xr37o13797h999/R82aNfHo0SMcOXIEz549U0zz/hFOmqJixYqlct1E6oB7mqS10tLScOLECSxcuBDt27eHra0tmjVrhqlTp8Ld3V0x3T8Pz969excymQw7duxA69atYWxsjKZNm+LmzZs4f/48nJ2dUaZMGXTt2hWPHz9WLKNdu3YYN26c0vp79OiBoUOHFphvyZIlaNCgAUxNTWFtbY0xY8bg5cuXAP4+1Dhs2DCkp6dDJpNBJpNh9uzZAPIeIk1MTET37t1RpkwZmJmZoV+/fnj06JFi/OzZs9G4cWMEBgaievXqkMvl6N+/P168eJFvrs9Zt0wmw9q1a+Hh4QETExPY29vj9OnT+Ouvv9CuXTuYmprCxcUFt2/fVlrH/v374eTkBCMjI9SsWRNz5sxBdnZ2gdsqPyEhIahTpw6MjY3Rvn173L17V2n806dPMWDAAFSrVg0mJiZo0KABtm3bphi/efNmmJub53kkVu/eveHj4wMAuHTpEtq3b4+yZcvCzMwMTk5OuHDhwmflJO3C0iStVaZMGZQpUwbBwcEFPiuwILNmzcKMGTNw8eJF6OnpYcCAAZg8eTJ++eUXHD9+HLdv38YPP/xQpHw6OjpYvnw5rl69it9//x1Hjx7F5MmTAQAtW7bEsmXLYGZmhuTkZCQnJ2PixIl5liEIAnr06IFnz54hKioK4eHhuH37Nry8vJSmu337NoKDg3HgwAEcOHAAUVFRWLBgQb65xK77vR9//BE+Pj6IjY1F3bp1MXDgQHz77beYOnWqomB8fX0V0x8+fBiDBw/Gv//9b1y/fh1r167Fpk2b8NNPP4nedklJSejVqxfc3NwQGxuLESNGYMqUKUrTvH37Fk5OTjhw4ACuXr2Kb775Bt7e3jh79iwAoG/fvsjJycG+ffsU8zx58gQHDhzAsGHDAACDBg1CtWrVcP78eURHR2PKlCnQ19cXnZO0UJGfz0Kkxnbu3CmUL19eMDIyElq2bClMnTpVuHTpktI0AIQ9e/YIgiAICQkJAgBh/fr1ivHbtm0TAAhHjhxRDPP39xfs7OwUP3/4mCJBEITu3bsrHnUmCIJga2srLF26tMCsO3bsEMzNzRU/BwQECHK5PM90/1xOWFiYoKurq/R4o2vXrgkAhHPnzgmCIAizZs0STExMlJ4rOGnSJKF58+YFZhGzbkH4e9vNmDFD8fPp06cFAMKGDRsUw7Zt2yYYGRkpfm7durUwf/58peUGBgYKVapUKTDPh6ZOnSrY29sLubm5imHff/+9AEB4/vx5gfO5ubkJEyZMUPw8evRooVu3boqfly1bJtSsWVOx3LJlywqbNm0SnYu0H/c0Sav17t0bDx8+xL59+9ClSxdERkbC0dERmzZt+uh8DRs2VPy5UqVKAIAGDRooDUtNTS1StoiICHTu3BlVq1ZF2bJl4ePjg6dPn+LVq1eilxEXFwdra2tYW1srhtWrVw/lypVDXFycYlj16tVRtmxZxc9VqlQpcv73xGyrt2/fIiMjAwAQHR2NuXPnKo4ElClTBiNHjkRycjJev34tap1xcXFo0aIFZDKZYpiLi4vSNDk5Ofjpp5/QsGFDmJubo0yZMggLC0NiYqJimpEjRyIsLAwPHjwAAAQEBGDo0KGK5fr5+WHEiBHo1KkTFixYkOcwM5U+LE3SekZGRujcuTN++OEHnDp1CkOHDsWsWbM+Os8/D8G9/wD9cFhubq7iZx0dHQgfPJr23bt3BS7/3r17cHNzg4ODA3bt2oXo6GisXLnyk/N9SBAEpeIoaPiHhxQ/zF8UYrYVAMX6cnNzMWfOHMTGxipeV65cwa1bt2BkZCRqnR9u6/wsXrwYS5cuxeTJk3H06FHExsaiS5cuyMrKUkzTpEkTNGrUCJs3b8bFixdx5coVpe+hZ8+ejWvXrsHd3R1Hjx5FvXr1sGfPHlEZSTvx7FkqderVq6fy6zIrVqyI5ORkxc85OTm4evUq2rdvn+/0Fy5cQHZ2NhYvXgwdnb9/d92xY4fSNAYGBsjJyfnoeuvVq4fExEQkJSUp9javX7+O9PR02NvbF/r9iFl3YTk6OiI+Ph61atUq9DLy+zs8c+aM0s/Hjx9H9+7dMXjwYAB/l/WtW7fybJcRI0Zg6dKlePDgATp16qS01w4AderUQZ06dTB+/HgMGDAAAQEB6NmzZ6Gzk2bjniZpradPn6JDhw74448/cPnyZSQkJODPP//EokWL0L17d5Wuq0OHDjh48CAOHjyIGzduYMyYMR+9yP6LL75AdnY2fv31V9y5cweBgYFYs2aN0jTVq1fHy5cvceTIETx58iTfQ5edOnVCw4YNMWjQIFy8eBHnzp2Dj48P2rZtC2dn50K/HzHrLqwffvgBmzdvVuzFxcXFYfv27ZgxY4boZYwaNQq3b9+Gn58f4uPjsXXr1jyH3GvVqoXw8HCcOnUKcXFx+Pbbb5GSkpJnWYMGDcKDBw+wbt06DB8+XDH8zZs38PX1RWRkJO7du4eTJ0/i/PnzRfplhDQfS5O0VpkyZdC8eXMsXboUbdq0gYODA2bOnImRI0dixYoVKl3X8OHDMWTIEEVh1ahRo8C9TABo3LgxlixZgoULF8LBwQFbtmyBv7+/0jQtW7bEqFGj4OXlhYoVK2LRokV5lvP+cpny5cujTZs26NSpE2rWrInt27cX6f2IWXdhdenSBQcOHEB4eDiaNm2KFi1aYMmSJbC1tVVMM3ToULRr167AZdjY2GDXrl3Yv38/GjVqhDVr1mD+/PlK08ycOROOjo7o0qUL2rVrh8qVK6NHjx55lmVmZobevXujTJkySuN1dXXx9OlT+Pj4oE6dOujXrx+6deuGOXPmFHUTkAaTCWK+HCAiKkHt2rVDu3btFNeHFrfOnTvD3t4ey5cvL5H1keZiaRKRWnnx4gXq1auHuLg4lClTpljX9ezZM4SFhWHQoEG4fv067OzsinV9pPl4IhARqZWyZcsiKSmpRNbl6OiI58+fY+HChSxMEoV7mkRERCLxRCAiIiKRWJpEREQisTSJiIhEYmkSERGJxNIkIiISiaVJREQkEkuTiIhIJJYmERGRSP8HQDDc7lMrU8gAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 450x500 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot(chd);"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.13.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}