{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "523e51a3",
   "metadata": {},
   "source": [
    "# Constrained Bayesian model: elasticity\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "4c003376",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-09-14T14:44:31.677158Z",
     "iopub.status.busy": "2026-09-14T14:44:31.675138Z",
     "iopub.status.idle": "2026-09-14T14:44:41.671097Z",
     "shell.execute_reply": "2026-09-14T14:44:41.671097Z"
    }
   },
   "outputs": [],
   "source": [
    "from pathlib import Path\n",
    "import sys\n",
    "\n",
    "PROJECT_ROOT = Path.cwd().resolve()\n",
    "while not (PROJECT_ROOT / 'data.py').exists():\n",
    "    if PROJECT_ROOT.parent == PROJECT_ROOT:\n",
    "        raise FileNotFoundError('Open this notebook from inside organized_uncertainty_analysis.')\n",
    "    PROJECT_ROOT = PROJECT_ROOT.parent\n",
    "if str(PROJECT_ROOT.parent) not in sys.path:\n",
    "    sys.path.insert(0, str(PROJECT_ROOT.parent))\n",
    "\n",
    "import os\n",
    "os.chdir(PROJECT_ROOT)\n",
    "from organized_uncertainty_analysis.bayesian.data_fit_gof import fit_full_trim\n",
    "from organized_uncertainty_analysis.bayesian import elasticity as el\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "7d369912",
   "metadata": {
    "execution": {
     "iopub.execute_input": "2026-09-14T14:44:41.679432Z",
     "iopub.status.busy": "2026-09-14T14:44:41.679432Z",
     "iopub.status.idle": "2026-09-14T14:44:49.372536Z",
     "shell.execute_reply": "2026-09-14T14:44:49.370522Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Loading SUV from: .\\suv_percentilesSLOthenUWM.mat\n",
      "Loading FLAGS from: .\\flags_combined.mat\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Dataset</th>\n",
       "      <th>Parameter</th>\n",
       "      <th>Elasticity_x50</th>\n",
       "      <th>Elasticity_s50</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>FULL</td>\n",
       "      <td>omega</td>\n",
       "      <td>-0.082300</td>\n",
       "      <td>0.710789</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>FULL</td>\n",
       "      <td>a</td>\n",
       "      <td>-1.935020</td>\n",
       "      <td>40.182491</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>FULL</td>\n",
       "      <td>b</td>\n",
       "      <td>1.973801</td>\n",
       "      <td>-42.461717</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>FULL</td>\n",
       "      <td>s</td>\n",
       "      <td>-1.944316</td>\n",
       "      <td>40.738767</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>FULL</td>\n",
       "      <td>k</td>\n",
       "      <td>2.700069</td>\n",
       "      <td>-36.685049</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>FULL</td>\n",
       "      <td>vartheta</td>\n",
       "      <td>2.944315</td>\n",
       "      <td>-41.738754</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>TRIM</td>\n",
       "      <td>omega</td>\n",
       "      <td>-0.046275</td>\n",
       "      <td>0.140809</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>TRIM</td>\n",
       "      <td>a</td>\n",
       "      <td>-0.211620</td>\n",
       "      <td>9.149163</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>TRIM</td>\n",
       "      <td>b</td>\n",
       "      <td>0.192150</td>\n",
       "      <td>-9.296605</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>TRIM</td>\n",
       "      <td>s</td>\n",
       "      <td>-0.206688</td>\n",
       "      <td>9.185517</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>TRIM</td>\n",
       "      <td>k</td>\n",
       "      <td>1.109265</td>\n",
       "      <td>-8.682383</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>TRIM</td>\n",
       "      <td>vartheta</td>\n",
       "      <td>1.206688</td>\n",
       "      <td>-10.185517</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "   Dataset Parameter  Elasticity_x50  Elasticity_s50\n",
       "0     FULL     omega       -0.082300        0.710789\n",
       "1     FULL         a       -1.935020       40.182491\n",
       "2     FULL         b        1.973801      -42.461717\n",
       "3     FULL         s       -1.944316       40.738767\n",
       "4     FULL         k        2.700069      -36.685049\n",
       "5     FULL  vartheta        2.944315      -41.738754\n",
       "6     TRIM     omega       -0.046275        0.140809\n",
       "7     TRIM         a       -0.211620        9.149163\n",
       "8     TRIM         b        0.192150       -9.296605\n",
       "9     TRIM         s       -0.206688        9.185517\n",
       "10    TRIM         k        1.109265       -8.682383\n",
       "11    TRIM  vartheta        1.206688      -10.185517"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAB6kAAALhCAYAAAANepGEAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAABuvAAAbrwFeGpEcAACPTklEQVR4nOzde5xVdb0//tcGhuswIKCiKEFeIMLMNB+Z4i0PpmjZ6djxG4mmZWWaEmj+LA1NzZDUk6jlPTWPpWUZiXnEOGqZF7oJeiQvCKKooFyG+zD794cP5sgBHC6z98yeeT4fj/1o7bU+67PeG9YY73nttVahWCwWAwAAAAAAAABl0K65CwAAAAAAAACg7RBSAwAAAAAAAFA2QmoAAAAAAAAAykZIDQAAAAAAAEDZCKkBAAAAAAAAKBshNQAAAAAAAABlI6QGAAAAAAAAoGyE1AAAAAAAAACUjZAaAAAAAAAAgLIRUgMAAAAAAABQNkJqAAAAAAAAAMpGSA0AAAAAAABA2QipAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAKBsOjR3AWyeQqHQ3CUAAAC0GcVisblLYDPomQEAAMpja/tlV1IDAAAAAAAAUDaupK5Qvs0PAABQOq7IrWx6ZgAAgNJoqn7ZldQAAAAAAAAAlI2QGgAAAAAAAICyEVIDAAAAAAAAUDZCagAAAAAAAADKRkgNAAAAAAAAQNkIqQEAAAAAAAAoGyE1AAAAAAAAAGUjpAYAAAAAAACgbDo0dwGUR7FYzJIlS7J48eKsWbOmucuhmbVv3z41NTXp3r17CoVCc5cDAAAAAABAG1IoFovF5i6CTbc2UNycv7a6urrMmjUrq1evLlVZVKiqqqoMGDAgHTr4vgoAALzblvReND9/bwAAAKXVVH1Xqw6p//KXv+S3v/1tpk2blpkzZ+bNN9/M4sWLU1NTk8GDB+fII4/M1772tfTq1WuLj/H6669n/PjxmTRpUmbPnp0uXbrkgx/8YE444YScfPLJTX6V6ub+xReLxbz88stZvnx5kqRdu3bp1KmTq2fbsGKxmJUrV6a+vj5J0rVr1/Tv3985AQAA7yLsrEz+3gAAAEpLSL0JTjvttFx99dUN7zt37pyqqqosWbKkYV2fPn1y7733Zr/99tvs+adNm5bDDz88CxYsSJJUV1dnxYoVqaurS5IMHz489957bzp16rSVn+R/be5f/LJly/Lyyy8nSXr06JG+ffumXTuPIm/r6uvr89prr2Xx4sVJkve9733p2rVrM1cFAAAth7CzMvl7AwAAKK2m6rtadVq577775rLLLstjjz2Wt99+O8uXL8/ixYuzZMmS3HLLLdl2220zf/78HHPMMVm0aNFmzb1o0aIcddRRWbBgQQYPHpwnn3wyS5YsydKlSzNx4sRUVVXlgQceyOjRo0v06TbN0qVLG5a32247ATVJ3rmifvvtt294/+7zBAAAAAAAAEqpVV9J3ZgHHngghx9+eJLk9ttvz8iRIzd53/POOy8XXXRRunTpkhkzZmTgwIHrbP/+97+fc889N+3bt88zzzyT3XffvUlq3txvJ8yZMye1tbXp2LFjdtlllyapgdbjhRdeyKpVq1JdXZ2dd965ucsBAIAWwxW5lcnfGwAAQGm5kroJfOxjH2tYfuWVVzZr31tvvTVJctxxx60XUCfJ6aefnurq6qxZsyY/+9nPtq7QrbD2ucMdOnRothpoudq3b5/kf88TAAAAAAAAKLU2HVI/8sgjDcubc5Xxc889l9mzZydJjjjiiA2Oqa6uzrBhw5K8c8U2tERrv+0CAAAAAAAA5dLmQuqVK1dm1qxZmThxYo4//vgkya677pqjjz56k+eYPn16w/LQoUM3Om7ttmeeeWYLqwUAAAAAAABoXdrMPaA7d+6clStXrrd+//33zx133JFOnTpt8lyvvvpqw3K/fv02Om7ttsWLF6e2tjbV1dWNzr2pV7bOnz9/k8atWLEia9asSX19fVavXr1J+9B21NfXZ82aNVmxYsUmn1MAAAAAAACwNdrMldR9+/bN9ttvn27dujWsO+SQQ3LllVemf//+mzXXkiVLGpa7du260XHv3vbufQAAAAAAAADaqjZzJfWsWbMalt94443cdtttufjii7PvvvvmO9/5Ti688MLmK+5disXie25fe6V1nz59Nmm+pUuXpr6+Pu3atUtVVdVW10fr0q5du7Rv3z6dO3fe5HMKAAAAAAAAtkabuZL63bbbbruMGTMm999/fwqFQr73ve9l0qRJm7x/9+7dG5aXLVu20XHv3vbufSrNvHnz8tRTT63z2tCt0yvdwQcfnEKh8J6vY445pmH82nVTp059z3lPPPHEFAqFnHjiieusnzVrVsMc7/4SxZbMBQAAQOnccsstjfaLhUIhnTp1yrbbbpu99torX/ziF/OrX/2qVfbPAAAAW6vNXEm9Ifvuu28OOOCAPPzww7nuuuty1FFHbdJ+O+64Y8Py3LlzU1NTs8Fxc+fOTZLU1NRs0vOom9OAc363eTvc/XppCtkMsy4dUZJ5d955543eAn7IkCElOSYAAACVb9WqVZk/f37mz5+fv/3tb7nlllvSv3//TJw4MUcffXRzlwcAANBitOmQOkn69euXJHn++ec3eZ+hQ4c2LE+fPj0f+MAHNjhu+vTpSQSbleakk07KuHHjmrsMAAAAWqjttttug3dMW7ZsWebPn5/Vq1c3rJs9e3Y+/elP57rrrsuXvvSlcpYJAADQYrXJ232/24svvphk827HPWjQoIYrbe+///4Njlm6dGkeeeSRJMnw4cO3skoAAACgpfjBD36Q559/fr3Xq6++2vD7gE996lMN44vFYk4//fTN+oI8AABAa9ZqQ+o1a9akWCy+55gpU6bkiSeeSPLO84g3x6hRo5Ikd9555wafJ3z11VentrY27du3z8iRIzdrbgAAAKAyVVVV5YADDshvfvObnHjiiQ3rV6xYkeuvv775CgMAAGhBWm1IPWfOnOy11175yU9+khdffHGdwHrOnDm59NJL8+lPfzrFYjG9evXK6NGj19l/3LhxKRQKKRQKGwyhx44dm759+2bZsmUZMWJEpk2bluSd509de+21Oe+885Ikp5xySnbffffSfVAAAACgRbroootSKBQa3v/xj39sxmoAAABajlb9TOq///3v+epXv5ok6dixY2pqarJ8+fIsXbq0YczAgQPzy1/+Mn379t2suXv06JFJkybl8MMPzzPPPJN99tkn3bt3z4oVKxqePTV8+PBcccUVTfeBAAAAgIrRr1+/bLvttnnjjTeSpOF/AQAA2rpWeyX1jjvumF/84hc59dRTs/fee6dPnz5ZvHhx6uvr079//xx99NG54YYbMmPGjOy1115bdIy99947M2bMyOjRo7Pbbrtl9erV6datWw444IBcf/31mTx5cjp16tTEn4xSu+CCCxquov+/LwAAANgca9asaVju1q1bM1YCAADQcrTaK6k7duyYY489Nscee+wW7T9u3LiMGzeu0XHbb799Lr/88lx++eVbdBxanp133jn9+/dv7jIAAACocP/85z+zYMGChvd77rlnM1YDAADQcrTakBq21EknnbRJX1AAAKByDTjnd81dwlabdemI5i4BeA91dXUZPXr0OutGjRrVTNUArYV/wwAArYWQGrZQu3btUl9fv86t2zakrq4uSdK+fftylAUAAEAzWbFiRebOnZs//vGPufLKK/PXv/61YduXv/zlHHrooc1YHQAAQMshpIYt1KNHj7z99tt5++2333Pc2u09e/YsQ1UAAACU2he/+MV88Ytf3KSxffv2zVlnnZUzzzyztEUBAABUkHbNXQBUqt133z1J8o9//OM9x63dPmjQoJLXBAAAQMsxaNCgfP/7389pp52Wdu38CgYAAGAtHRJsoSOOOCJJ8p//+Z9ZvXr1Bsc8+OCDeeWVV9KuXbv8y7/8SznLAwAAoES222677LLLLuu9BgwYsM5dtJ577rl88YtfzK677popU6Y0X8EAAAAtjJAattDXvva19O7dO88//3xGjhyZN954Y53tDz/8cE444YQkyciRIzNw4MDmKBMAAIAm9oMf/CDPP//8eq+XXnopb7/9dubOnZuJEydm++23T5LMmTMnn/zkJ/Pb3/62mSsHAABoGYTUsIW222673H333enZs2fuuuuu7Ljjjhk6dGg+9rGPpV+/fjnooIPy6quv5sADD8zEiRPfc66PfOQj6dOnz0ZfixcvXmf8nXfe+Z7jv/Od75TyowMAAPAedtxxx3z961/PX//617zvfe9LktTV1eXEE0/Mm2++2czVAQAANL8OzV0ALcOsS0fkxRdfzFtvvZUkad++ffbcc89Gn5n11ltv5cUXX2x4v9tuu6VHjx4lrbUlOfjggzNjxoz86Ec/yu9///u88MILWb58eXr16pXhw4fn85//fEaOHJkOHd77R+3tt99+z+319fXrvF+5cmVWrly50fG1tbWb/iEAAAAoiR122CHf//738/nPfz7JOz30tddem/PPP7+ZKwMAAGheQmqSvPON7ncHpWvWrMlf/vKXzZ5n/vz5FRtST506dYv223HHHXPppZfm0ksv3az9BgwYkGKxuFn73HLLLbnllls2ax8AAACaz9FHH5327dtnzZo1SZLf/e53QmoAAKDNc7tvkrzzbe7NDUw3ZOHChVm9enUTVAQAAACVr7q6Or179254/9xzzzVjNQAAAC2DK6lJ8s4V0Gu1a9cuVVVVm7X/2ltPF4vFLFiwIH379m3S+gAAAKBSvftL4cuXL2/GSgAAAFoGITVZtmxZli1b1vB+xx133OyQ+dlnn83SpUuTvBN4C6kBAAAgWbJkSRYsWNDwftttt23GagAA4L0NOOd3zV3CVpt16YjmLoFN4HbfrHMVdZJss802mz1Hr169GpZXrFiR2trara4LAAAAKt1vfvOb1NfXN7zfc889m7EaAACAlkFI3cbV19ev843u6urqdOrUabPneXdInawffAMAAEBbM3fu3Jx77rnrrPvsZz/bTNUAAAC0HELqNu7tt9/OmjVrGt7/37B5U1VVVaWmpqbh/VtvvbXOvAAAANBWvPrqq5k4cWI+8pGPZM6cOQ3rhw4dmlGjRjVjZQAAAC2DZ1K3ce++4rlQKGzRrb7X6tWrVxYvXpzknSu033777fTp02erawQAAICW5Fvf+lYuuuii9davWbMmCxcuzMKFC9fbttNOO+U3v/lNOnTwqxgAAACdURu2cuXKLFmypOF9TU1Nqqqqtni+bbbZJi+//HKKxWKSdwJwITUAAACtzRtvvJE33nhjk8YWCoV87nOfy3/8x39k++23L3FlAAAAlUFI3Yb93+dGb+mtvtdq3759evTo0fCN8dra2qxYsSKdO3feqnkBAACgErRr1y7V1dXp1atXPvjBD2a//fbL//t//y/vf//7m7s0AACAFkVI3Yb169cv/fr1a9I5d9111yadDwAAAJrbiSeemBNPPLG5ywAAAGg12jV3AQAAAAAAAAC0HUJqAAAAAAAAAMpGSA0AAAAAAABA2QipAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAKBshNS0aYVCYbNfAwYMaNh/wIAB623v2LFjdtxxx3zqU5/Kvffeu9Fjz5o1q2GfWbNmrbPtlltuadjWvn37vPLKK+/5OYYMGdIw/sQTT9yKPxEAAAAAAAAorQ7NXQAtxLgezV3B5hu3aKun2H///ddbt2jRokyfPn2j23fYYYf11u22227ZbrvtkiRLlizJCy+8kN/+9rf57W9/m9NOOy1XXXXVFtdYX1+fO+64I2efffYGtz/11FN59tlnt3h+AAAAAKB1GHDO75q7hK0269IRzV0CAGUgpKZNe/TRR9dbN3Xq1BxyyCEb3b4h55577jpXMC9btixjx47Ntddem4kTJ+Yzn/lMDj300M2ub9ddd82LL76Y2267baMh9W233ZYkGTx4cP7nf/5ns48BAAAAAAAA5eR231ACXbt2zY9+9KOGW4PffffdWzRPv379cuihh2b69On529/+tt72urq63HnnnampqcmnPvWpragYAAAAAAAAykNIDSXSoUOHfPjDH06SvPTSS1s8z/HHH5/kf6+Yfrf7778/b7zxRv7t3/4tXbp02eJjAAAAbIpLL700hUKh4fVeXn/99YwZMyaDBg1Kly5d0qtXrwwbNiw33HBDisVimSoGAACgJRJSQwktW7YsSdKtW7ctnuNf//Vf061bt9xxxx1Zs2bNOtvWBtdrg2wAAIBSee6553LBBRds0thp06blgx/8YC6//PLMnDkzHTp0yJIlS/Loo4/my1/+cj75yU9m5cqVJa4YAACAlkpIDSUyf/78PPHEE0mSPffcc4vnqa6uzmc+85nMmzcvDz74YMP6RYsW5d57703//v1z0EEHbXW9AAAAG1NfX5+TTz45K1asyH777feeYxctWpSjjjoqCxYsyODBg/Pkk09myZIlWbp0aSZOnJiqqqo88MADGT16dJmqBwAAoKURUkMTq62tzaOPPpqjjz46CxcuTO/evXPKKads1ZwbuuX3XXfdlRUrVuQLX/hCo7fZAwAA2BpXXXVV/vjHP2bkyJEZPnz4e46dMGFC5s2bly5duuS+++7LPvvskyTp2LFjvv71rzdcjX3ddddl5syZJa8dAACAlkdIDU3gi1/8YsMz2bp3755hw4blz3/+c/bff/889NBD2X777bdq/k984hPZcccdc88996S2tjaJW30DAADl8dJLL+Xb3/52evfunSuuuKLR8bfeemuS5LjjjsvAgQPX23766aenuro6a9asyc9+9rMmrxcAAICWr0NzFwCtwW677ZbtttsuSfLmm2/mn//8Z4rFYnbdddd84AMf2Or527dvn89//vOZMGFCfvWrX+Wggw7KI488ko9+9KMZPHjwVs8PAACwMV/+8pezdOnSXHPNNdl2223fc+xzzz2X2bNnJ0mOOOKIDY6prq7OsGHDMnny5DzwwAOb/JxraAsGnPO75i5hq826dERzlwAAQAVwJTU0gXPPPTePPvpoHn300Tz33HOZMWNGdt999/z0pz/Naaed1iTHePctv2+77bYUi0VXUQMAACV1/fXXZ8qUKTnssMMyatSoRsdPnz69YXno0KEbHbd22zPPPLP1RQIAAFBxXEkNJfCBD3wgv/zlL7PXXnvluuuuywknnJCPf/zjWzXnhz70oXzoQx/KQw89lP/5n/9Jhw4dctxxxzVRxQAAAOuaO3duzjrrrHTp0iU/+clPNmmfV199tWG5X79+Gx23dtvixYtTW1ub6urqTZq/UChs0rj58+dv0jhoaXatqW/uEraan7/Sco7QGOcIsLX8d4RycSU1lMjQoUNz4oknJnnnSuumMGrUqNTX1+eVV17JEUcc0eit9gAAALbUV77ylSxatCjjxo3L+9///k3aZ8mSJQ3LXbt23ei4d2979z4AAAC0Da6khhI655xzcvPNN+e///u/88gjj2TYsGFbNd/nP//5TJ48OUny9a9/vSlKBAAAWM/tt9+e3/3ud/nwhz+cb37zm81dToNisfie29dead2nT59ylANN7vnFlX89iZ+/0nKO0BjnCLC1/HeEcqn8Mw1asF122SWf+9znkiQXXXTRVs+3ww475MEHH8yDDz6Yww8/fKvnAwAA+L/eeOONnHnmmWnfvn2uv/76dOiw6d9v7969e8PysmXLNjru3dvevQ8AAABtg5AaSuzcc89NoVDIAw88kCeffLK5ywEAAHhP3/rWt7JgwYKccsopGTx4cGpra9d5rVq1qmHs/1234447NmybO3fuRo+xdltNTc0mP48aAACA1kNIDSU2dOjQHH300Uma5mpqAACAUnrppZeSJNdee226d+++3uv73/9+w9i1684+++wk7/Q/a02fPn2jx1i7bciQIaX4CAAAALRwnknNO8Ytau4KWoyDDz640eecrTVr1qxNGveb3/xmvXUDBgzY6HFOPPHEnHjiiZs091rjxo3LuHHjNmsfAACApjRo0KD0798/s2fPzv33359jjz12vTFLly7NI488kiQZPnx4uUsEAACgBXAlNQAAANBg6tSpKRaLG31997vfbRi7dt2VV17ZsG7UqFFJkjvvvHODX+y9+uqrU1tbm/bt22fkyJGl/jgAAAC0QEJqAAAAoMmMHTs2ffv2zbJlyzJixIhMmzYtSbJq1apce+21Oe+885Ikp5xySnbffffmLBUAAIBm4nbfAAAAQJPp0aNHJk2alMMPPzzPPPNM9tlnn3Tv3j0rVqzI6tWrk7xzm+8rrriimSsFAACgubiSGgAAAGhSe++9d2bMmJHRo0dnt912y+rVq9OtW7cccMABuf766zN58uR06tSpucsEAACgmQipAQAAgE02bty4hmdRv5ftt98+l19+eWbOnJnly5fn7bffziOPPJIvfelLadfOryMAAADaMl0hAAAAAAAAAGUjpAYAAAAAAACgbITUAAAAAAAAAJSNkBrasMaeIQcAAAAAAABNTUjdyrVr985fcV1dXTNXQku0Zs2aJP97ngAAAAAAAECpdWjuAiitzp07p7a2NqtWrUpdXV06dPBXzjvq6uqyatWqJO+cJwDQmgw453fNXcJWm3XpiOYuAQAAAABKwuWTrVy3bt0alt94443U19c3YzW0FPX19Xn99dcb3r/7PAEAAAAAAIBSclltK9elS5d07do1y5Yty6JFi7JkyZJ06tQphUKhuUujmRSLxaxcubLhCwtdu3ZNly5dmrkqAAAAAAAA2gohdStXKBTSr1+/zJo1K6tXr059fX2WL1/e3GXRQlRVVaVfv36+tAAAAAAAAEDZCKnbgA4dOmSXXXbJkiVLsmTJktTV1TV3STSzDh06pHv37unevbuAGgAAAAAAgLISUrcRhUIhNTU1qampae5SAAAAAAAAgDasXXMXAAAAAAAAAEDbIaQGAAAAAAAAoGyE1AAAAAAAAACUjZAaAAAAAAAAgLIRUgMAAAAAAABQNkJqAAAAAAAAAMpGSA0AAAAAAABA2QipAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAKBshNQAAAAAAAAAlI2QGgAAAAAAAICyEVIDAAAAAAAAUDZCagAAAAAAAADKRkgNAAAAAAAAQNkIqQEAAAAAAAAoGyE1AAAAAAAAAGUjpAYAAAAAAACgbITUAAAAAAAAAJSNkBoAAAAAAACAshFSAwAAAAAAAFA2QmoAAAAAAAAAykZIDQAAAAAAAEDZCKkBAAAAAAAAKBshNQAAAAAAAABlI6QGAAAAAAAAoGyE1AAAAAAAAACUTYfmLgAAAAAAAIDSG3DO75q7hK0269IRzV0C0ARcSQ0AAAAAAABA2bTqkHrBggW5+eab84UvfCFDhgxJt27d0qlTp+y000455phjcs8992zx3OPGjUuhUGj09fzzzzfhJwIAAAAAAACobK36dt99+/ZNXV1dw/vOnTunqqoqc+fOzdy5c/Ob3/wmRxxxRO6+++507dp1i45RVVWVXr16bXR7hw6t+o8YAAAAAAAAYLO06iup6+rqsu++++aaa67JCy+8kOXLl6e2tjYvvfRSTj755CTJ5MmT85WvfGWLj/Hxj3888+bN2+hrwIABTfRpAAAAAAAAACpfq77M96GHHsohhxyy3voBAwbkhhtuSIcOHfKTn/wkt99+ey655JLsvPPOzVAlAAAAAAAAQNvRqq+k3lBA/W5rr6ZOkqeeeqrU5QAAAAAAAAC0ea06pG5M586dG5bXrFnTjJUAAAAAAAAAtA2t+nbfjZk6dWrD8h577LFFc8yYMSNDhw7NCy+8kPbt26dfv3458MADc+qpp2avvfba7PkKhcImjZs/f/5mzw0A0FbsWlPf3CVsNf/eKy3nCAAAAEDzabNXUi9cuDDf//73kyTDhg3LoEGDtmie+fPn59lnn03Xrl2zcuXKzJw5MzfccEP23nvvfOc732nKkgEAAAAAAAAqXpu8krq+vj7HH398XnvttXTq1ClXXXXVZs+x2267Zfz48fn0pz+dgQMHpqqqKqtWrcrUqVNz7rnnZtq0abn44ouzzTbbZMyYMZs8b7FYfM/ta6+07tOnz2bXDADQVjy/uPK/i+nfe6XlHAEAAABoPpX/m5ktcMYZZ2TSpElJkmuuuSZ77rnnZs8xcuTInHXWWdl9991TVVWVJOnYsWOGDx+eRx99NB/96EeTJOPGjcuiRYuarngAAAAAAACACtbmQuqxY8dm4sSJSZIrrrgiJ510UpMfo3PnzrnkkkuSJLW1tZkyZUqTHwMAAAAAAACgErWpkPrss8/OD3/4wyTJZZddljPPPLNkx9pvv/0all988cWSHQcAAAAAAACgkrSZZ1KfddZZmTBhQpJk/PjxGTt2bDNXBAAAAAAAAND2tImQeuzYsQ1XUI8fPz5nnXVWyY/55z//uWF54MCBJT8eAAAAAAAAQCVo9SH1uwPqCRMmZMyYMVs9Z7FYTKFQ2Oj2lStX5tvf/naSpFu3bvnEJz6x1ccEAAAAAAAAaA1a9TOpv/WtbzUE1JdffvlmBdTjxo1LoVBIoVDIrFmz1tn28MMP57DDDsvtt9+eV155pWH96tWrM2XKlAwbNiyPP/54kuT8889Pz549t/qzAAAAQLn85S9/yQUXXJBPfepTGTx4cHr37p2qqqr07t07+++/fy6++OK89dZbG9z33f30e72ef/75Mn8qAAAAWopWeyX17NmzM378+CRJu3bt8oMf/CA/+MEPNjp+7Nixm/yc6mKxmClTpmTKlClJki5duqRbt25ZtGhRVq9e3XDMc845J2efffZWfhIAAAAor5tuuilXX311w/vOnTunS5cueeutt/KnP/0pf/rTn3LllVfm3nvvzX777bfBOaqqqtKrV6+NHqNDh1b7KwkAAAAa0Wo7wvr6+nWWX3/99fccX1tbu8lz77HHHpkwYUIee+yxPP3005k/f34WLlyYrl27ZsiQIRk2bFhOOeWU7LHHHltcPwAAADSXfffdNwMGDMgBBxyQwYMHN9whrLa2Nr/85S9z1lln5c0338wxxxyTmTNnpkePHuvN8fGPfzxTp04tb+EAAABUhFYbUg8YMCDFYnGL9x83blzGjRu3wW29e/dukmdbAwAAQEs0atSoDa6vrq7OCSeckB122CGHH3543njjjUyaNCkjR44sc4UAAABUslb9TGoAAACg6X3sYx9rWH7llVeasRIAAAAqkZAaAAAA2CyPPPJIw/Iuu+zSjJUAAABQiYTUAAAAQKNWrlyZWbNmZeLEiTn++OOTJLvuumuOPvroDY6fMWNGhg4dmi5duqS6ujqDBg3Kl7/85fz1r38tZ9kAAAC0QK32mdQAAADA1uvcuXNWrly53vr9998/d9xxRzp16rTB/ebPn5+33norPXv2zOLFizNz5szMnDkzN954Y84999xcdNFFm11LoVDYpHHz58/f7LmhJdi1pr65S9hqfv5KyzlCY5wjNMY5QmOcI5SLK6kBAACAjerbt2+23377dOvWrWHdIYcckiuvvDL9+/dfb/xuu+2W8ePH57nnnsuKFSuyYMGCLF26NL///e+z9957p1gs5uKLL84Pf/jDcn4MAAAAWhBXUgMAAAAbNWvWrIblN954I7fddlsuvvji7LvvvvnOd76TCy+8cJ3xI0eOXG+Ojh07Zvjw4TnwwANz4IEH5sknn8y4cePypS99KT169NjkWorF4ntuX3uldZ8+fTZ5TmhJnl9c+deT+PkrLecIjXGO0BjnCI1xjlAulX+mAQAAAGWx3XbbZcyYMbn//vtTKBTyve99L5MmTdrk/Tt37pxLLrkkSVJbW5spU6aUqlQAAABaMCE1AAAAsFn23XffHHDAAUmS6667brP23W+//RqWX3zxxSatCwAAgMogpAYAAAA2W79+/ZIkzz//fDNXAgAAQKURUgMAAACbbe1V0N27d9+s/f785z83LA8cOLBJawIAAKAyCKkBAACABmvWrEmxWHzPMVOmTMkTTzyRJDn44IMb1je238qVK/Ptb387SdKtW7d84hOf2LpiAQAAqEhCagAAAKDBnDlzstdee+UnP/lJXnzxxXWC5zlz5uTSSy/Npz/96RSLxfTq1SujR49u2P7www/nsMMOy+23355XXnmlYf3q1aszZcqUDBs2LI8//niS5Pzzz0/Pnj3L9rkAAABoOTo0dwEAAABAy/L3v/89X/3qV5MkHTt2TE1NTZYvX56lS5c2jBk4cGB++ctfpm/fvg3risVipkyZkilTpiRJunTpkm7dumXRokVZvXp1kqRdu3Y555xzcvbZZ5fxEwEAANCSCKkBAACABjvuuGN+8YtfZOrUqXn88cfz2muvZf78+Wnfvn369++fPffcM5/+9Kfz+c9/Pl26dFln3z322CMTJkzIY489lqeffjrz58/PwoUL07Vr1wwZMiTDhg3LKaeckj322KOZPh0AAAAtgZAaAAAAaNCxY8cce+yxOfbYYzd73969e2fMmDElqAoAAIDWxDOpAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAKBshNQAAAAAAAAAlI2QGgAAAAAAAICyEVIDAAAAAAAAUDZCagAAAAAAAADKRkgNAAAAAAAAQNkIqQEAAAAAAAAoGyE1AAAAAAAAAGUjpAYAAAAAAACgbITUAAAAAAAAAJSNkBoAAAAAAACAshFSAwAAAAAAAFA2QmoAAAAAAAAAykZIDQAAAAAAAEDZCKkBAAAAAAAAKBshNQAAAAAAAABlI6QGAAAAAAAAoGyE1AAAAAAAAACUjZAaAAAAAAAAgLIRUgMAAAAAAABQNkJqAAAAAAAAAMpGSA0AAAAAAABA2QipAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAKBshNQAAAAAAAAAlI2QGgAAAAAAAICyEVIDAAAAAAAAUDZCagAAAAAAAADKRkgNAAAAAAAAQNkIqQEAAAAAAAAoGyE1AAAAAAAAAGUjpAYAAAAAAACgbITUAAAAAAAAAJSNkBoAAAAAAACAshFSAwAAAAAAAFA2QmoAAAAAAAAAykZIDQAAAAAAAEDZCKkBAAAAAAAAKBshNQAAAAAAAABlI6QGAAAAAAAAoGyE1AAAAAAAAACUjZAaAAAAAAAAgLIRUgMAAAAAAABQNkJqAAAAAAAAAMpGSA0AAAAAAABA2QipAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAKBshNQAAAAAAAAAlI2QGgAAAAAAAICyEVIDAAAAAAAAUDZCagAAAAAAAADKRkgNAAAAAAAAQNkIqQEAAIB1/OUvf8kFF1yQT33qUxk8eHB69+6dqqqq9O7dO/vvv38uvvjivPXWW+85x+uvv54xY8Zk0KBB6dKlS3r16pVhw4blhhtuSLFYLNMnAQAAoCXq0NwFAAAAAC3LTTfdlKuvvrrhfefOndOlS5e89dZb+dOf/pQ//elPufLKK3Pvvfdmv/32W2//adOm5fDDD8+CBQuSJNXV1VmyZEkeffTRPProo7nrrrty7733plOnTmX7TAAAALQcrqQGAAAA1rHvvvvmsssuy2OPPZa33347y5cvz+LFi7NkyZLccsst2XbbbTN//vwcc8wxWbRo0Tr7Llq0KEcddVQWLFiQwYMH58knn8ySJUuydOnSTJw4MVVVVXnggQcyevToZvp0AAAANDchNQAAALCOUaNGZezYsfnYxz6Wnj17Nqyvrq7OCSeckNtvvz1J8sYbb2TSpEnr7DthwoTMmzcvXbp0yX333Zd99tknSdKxY8d8/etfzwUXXJAkue666zJz5szyfCAAAABaFCE1AAAAsFk+9rGPNSy/8sor62y79dZbkyTHHXdcBg4cuN6+p59+eqqrq7NmzZr87Gc/K22hAAAAtEhCagAAAGCzPPLIIw3Lu+yyS8Pyc889l9mzZydJjjjiiA3uW11dnWHDhiVJHnjggRJWCQAAQEvVobkLAAAAAFq+lStX5rXXXsukSZNy/vnnJ0l23XXXHH300Q1jpk+f3rA8dOjQjc41dOjQTJ48Oc8888xm1VAoFDZp3Pz58zdrXmgpdq2pb+4Stpqfv9JyjtAY5wiNcY7QGOcI5SKkBgAAADaqc+fOWbly5Xrr999//9xxxx3p1KlTw7pXX321Yblfv34bnXPttsWLF6e2tjbV1dVNWDEAAAAtXclC6gsvvDAnn3zyezalAAAAwLpaWj/dt2/frFixIrW1tVm6dGmS5JBDDsn48ePTv3//dcYuWbKkYblr164bnfPd25YsWbLJIXWxWHzP7WuvtO7Tp88mzQctzfOLK//JfH7+Sss5QmOcIzTGOUJjnCOUS8nOtHHjxmXgwIE55phjct999zXaSAIAAAAtr5+eNWtW5s2bl9ra2rz++uuZMGFC/va3v2XfffdtuO03AAAAbI6ShdT9+/dPXV1d7r333hx99NEZMGBAvve972Xu3LmlOiQAAABUvJbcT2+33XYZM2ZM7r///hQKhXzve9/LpEmTGrZ37969YXnZsmUbnefd2969DwAAAG1DyULql156KZMnT85nPvOZdOjQIXPmzFnn2+CTJ09u9m+DAwAAQEtTCf30vvvumwMOOCBJct111zWs33HHHRuW3ytUX7utpqbG86gBAADaoJKF1IVCIYcffnh++ctfZs6cObnkkkuyyy67NHwb/KijjmpR3wYHAACAlqBS+um1z8x+/vnnG9YNHTq0YXn69Okb3XfttiFDhpSoOgAAAFqysjz9fLvttss555yTmTNn5sEHH8znPve5VFVVtchvgwMAAEBL0ZL76RdffDHJurfrHjRoUPr3758kuf/++ze439KlS/PII48kSYYPH17iKgEAAGiJyhJSv9uhhx6aO++8M3Pnzs0Pf/jD7L777qmrq8tvf/vbhm+DX3LJJVmwYEG5SwMAAIAWq1z99Jo1axoNvKdMmZInnngiSXLwwQevs23UqFFJkjvvvDOzZs1ab9+rr746tbW1ad++fUaOHLlVtQIAAFCZyh5Sr7Vo0aK88cYbefvtt1MoFFIsFlMsFjNnzpycd955GThwYK688srmKg8AAABapFL303PmzMlee+2Vn/zkJ3nxxRfXCaznzJmTSy+9NJ/+9KdTLBbTq1evjB49ep39x44dm759+2bZsmUZMWJEpk2bliRZtWpVrr322px33nlJklNOOSW77777FtcJAABA5epQzoPV1dXlV7/6Va677rpMnTq1oZHeYYcdcvLJJ+fYY4/NAw88kB//+Md54YUXMmbMmHTu3Dlf/epXy1kmAAAAtCjl7qf//ve/N+zbsWPH1NTUZPny5Vm6dGnDmIEDB+aXv/xl+vbtu86+PXr0yKRJk3L44YfnmWeeyT777JPu3btnxYoVWb16dZJ3bvN9xRVXbOGfBgAAAJWuLFdSz5w5M2eddVb69euX//f//l8eeuihFIvFHHroobnrrrsye/bsXHjhhdljjz0yZsyYPPfcczn33HNTLBZz1VVXbfFxFyxYkJtvvjlf+MIXMmTIkHTr1i2dOnXKTjvtlGOOOSb33HPPVn+2119/PWPGjMmgQYPSpUuX9OrVK8OGDcsNN9zg+doAAABslebop3fcccf84he/yKmnnpq99947ffr0yeLFi1NfX5/+/fvn6KOPzg033JAZM2Zkr7322uAce++9d2bMmJHRo0dnt912y+rVq9OtW7cccMABuf766zN58uR06tRpa/5oAAAAqGCFYomS1FWrVuXuu+/O9ddfn4cffjhJGm4FduKJJ+arX/1qdt11143uXywW07Nnz6xcuTIrVqzYohqqqqpSV1fX8L5z585p3779Ot/8PuKII3L33Xena9eumz3/tGnTcvjhhzc876u6ujorVqxoOObw4cNz7733NmnjXSgUkkQADgDwHgac87vmLmGrzbp0RHOX0Ko5R2hMc/ZeLaGfrlR6Ziqd/3+iMc4RGuMcoTHOERrjHKExTdV3lexK6n79+uX444/Pf//3f6dYLOZjH/tYfvrTn2bu3LmZMGHCezbUyTsfcJtttmm4FdiWqKury7777ptrrrkmL7zwQpYvX57a2tq89NJLOfnkk5MkkydPzle+8pXNnnvRokU56qijsmDBggwePDhPPvlklixZkqVLl2bixImpqqrKAw88sN6zuQAAAOC9tIR+GgAAAEqpZM+kXrBgQbp3756RI0fma1/7WvbYY4/NnuOHP/xhamtrt7iGhx56KIcccsh66wcMGJAbbrghHTp0yE9+8pPcfvvtueSSS7Lzzjtv8twTJkzIvHnz0qVLl9x3330ZOHBgknee1fX1r389ixcvzrnnnpvrrrsuZ555Znbfffct/hwAAAC0HS2hnwYAAIBSKtmV1Ndee21effXVXHPNNVvUUCfJZz/72ZxwwglbXMOGAup3W3s1dZI89dRTmzX3rbfemiQ57rjjGgLqdzv99NNTXV2dNWvW5Gc/+9lmzQ0AAEDb1RL6aQAAACilkoXUX/nKV9KtW7dSTd8kOnfu3LC8Zs2aTd7vueeey+zZs5O880zrDamurs6wYcOSJA888MBWVAkAAEBbUgn9NAAAAGyNkt3uu127dtlhhx0yd+7cTRo/cODAzJkzJ3V1daUqaT1Tp05tWN6cb6dPnz69YXno0KEbHTd06NBMnjw5zzzzzCbPvfZh442ZP3/+Js8JANDW7FpT39wlbDX/3ist5wgtWSX00wAAALA1SnYldZIUi8WSjt8aCxcuzPe///0kybBhwzJo0KBN3vfVV19tWO7Xr99Gx63dtnjxYs8CAwAAYJO15H4aAAAAtlbJrqTeXKtWrUq7diXNzBvU19fn+OOPz2uvvZZOnTrlqquu2qz9lyxZ0rDctWvXjY5797YlS5akurq60bkb+8XC2iut+/Tp0+hcAABt1fOLy/PvylLy773Sco7QmpSznwYAAICm0CK62IULF+aNN95Iz549y3K8M844I5MmTUqSXHPNNdlzzz3LclwAAABoSuXupwEAAKApNNmV1P/4xz/yt7/9bZ11y5cvz6233rrRfYrFYhYuXJi777479fX1+chHPtJU5WzU2LFjM3HixCTJFVdckZNOOmmz5+jevXvD8rJly1JTU7PBccuWLdvgPgAAALBWpfTTAAAA0FSaLKS+5557cuGFF66zbvHixfniF7/Y6L7FYjGFQiHf/OY3m6qcDTr77LPzwx/+MEly2WWX5cwzz9yieXbccceG5blz5240pJ47d26SpKamZpNu9Q0AAEDbUwn9NAAAADSlJgupe/bsmf79+ze8f/nll9OuXbvstNNOG92nXbt2qampydChQ3PKKadk2LBhTVXOes4666xMmDAhSTJ+/PiMHTt2i+caOnRow/L06dPzgQ98YIPjpk+fniQZMmTIFh8LAACA1q2l99MAAADQ1JospD7jjDNyxhlnNLxv165dtt1227z00ktNdYgtNnbs2IYrqMePH5+zzjprq+YbNGhQ+vfvn9mzZ+f+++/Pscceu96YpUuX5pFHHkmSDB8+fKuOB9DWDDjnd81dwlabdemI5i4BAKgQLbmfBgAAgFJoV6qJv/vd72bMmDGlmn6TvTugnjBhwlYH1GuNGjUqSXLnnXdm1qxZ622/+uqrU1tbm/bt22fkyJFNckwAAABav5bSTwMAAECptOqQ+lvf+lZDQH355ZdvVj3jxo1LoVBIoVDYYAg9duzY9O3bN8uWLcuIESMybdq0JMmqVaty7bXX5rzzzkuSnHLKKdl99923/sMAAADQJrSEfhoAAABKqWQhdXObPXt2xo8fn+SdW6X94Ac/SN++fTf6Wvu86k3Vo0ePTJo0Kb17984zzzyTffbZJzU1Namurs6pp56aVatWZfjw4bniiitK8fEAAAAAAAAAKlKTPJP6wgsvTJL06dMnp5566jrrNtf555/fFCWlvr5+neXXX3/9PcfX1tZu9jH23nvvzJgxIz/4wQ8yadKkzJkzJ926dcvQoUNzwgkn5KSTTkq7dq32ewAAAABspZbYTwMAAECpFYrFYnFrJ2nXrl0KhUIGDRqUZ555Zp11m2vNmjVbW06rtvbPtAn+2gBarAHn/K65S9hqsy4d0dwlQJvmvyM0xjlCY8rVe+mnm5aemUrn/59ojHOExjhHaIxzhMY4R2hMU/VdTXIl9YEHHphCoZD+/fuvtw4AoKn5xzIArYV+GgAAgLaoSULqqVOnbtI6AAAA4H/ppwEAAGiLmiSkBthUreHqx8QVkAAAAAAAAFuqXXMXAAAAAAAAAEDbUbKQesmSJbn33nvzhz/8odGxDz30UO69997U1taWqhwAAACoCPppAAAAWruShdR33HFHPvOZz2Ty5MmNjr3rrrvymc98Jj//+c9LVQ4AAABUBP00AAAArV3JQup77rknSXLcccc1OvbEE09MsVjML3/5y1KVAwAAABVBPw0AAEBrV7KQ+rnnnkuSDBkypNGxH/rQh9bZBwAAANoq/TQAAACtXclC6nnz5qVnz57p3Llzo2O7dOmSbbbZJvPmzStVOQAAAFAR9NMAAAC0diULqbt27Zra2tqsWbOm0bF1dXWpra1Nx44dS1UOAAAAVAT9NAAAAK1dyULq3XbbLXV1dZkyZUqjY6dMmZLVq1dnl112KVU5AAAAUBH00wAAALR2JQupjzzyyBSLxZx11llZsmTJRsfV1tbmrLPOSqFQyJFHHlmqcgAAAKAi6KcBAABo7UoWUp922mnZZpttMn369Hz0ox/NPffck+XLlzdsX758eX71q19ln332yfTp09OjR4+cccYZpSoHAAAAKoJ+GgAAgNauQ6km7tWrV/7zP/8zxxxzTGbOnJl/+7d/S/v27dOnT58UCoW8+eabWbNmTYrFYjp37pxf/OIX6d27d6nKAQAAgIqgnwYAAKC1K9mV1EkyfPjw/PGPf8z++++fYrGYurq6zJs3L6+99lrq6upSLBZz4IEH5rHHHsthhx1WylIAAACgYuinAQAAaM1KdiX1WnvttVceeeSRPP/88/nTn/6UefPmJUl22GGHfPzjH88uu+xS6hIAAACg4uinAQAAaK1KHlKvteuuu2bXXXct1+EAAACgVdBPAwAA0NqU9HbfAAAAAAAAAPBuQmoAAAAAAAAAyqZJbvd96KGHJkne97735eabb15n3eYoFAqZMmVKU5QEAAAALZ5+GgAAgLaoSULqqVOnJkkGDx683rrNUSgUmqIcAAAAqAj6aQAAANqiJgmpv/vd7yZJ+vTps946AAAAYMP00wAAALRFTRpSN7YOAAAA+F/6aQAAANqids1dAAAAAAAAAABtR8lC6gsvvDCXX375Jo//0Y9+lAsvvLBU5QAAAEBF0E8DAADQ2pUspB43blwmTJiwyeOvuOKKXHDBBaUqBwAAACqCfhoAAIDWzu2+AQAAAAAAACibFhNSv/XWW+ncuXNzlwEAAAAVRT8NAABApWkRIfVdd92VJUuWpH///s1dCgAAAFQM/TQAAACVqENTTfQf//Ef+Y//+I911r355pt5//vfv9F9isViFi5cmMWLF6dQKGTEiBFNVQ4AAABUBP00AAAAbU2ThdQLFy7MrFmz1lm3Zs2a9dZtzCc+8Ymcf/75TVUOAAAAVAT9NAAAAG1Nk4XUxxxzTAYMGJDknW90n3TSSenRo0euvPLKje7Trl271NTUZOjQodlll12aqhQAAACoGPppAAAA2pomC6n33HPP7Lnnng3vTzrppHTp0iUnnHBCUx0CAAAAWh39NAAAAG1Nk4XU/1d9fX2ppgYAAIBWSz8NAABAa9euOQ/+9ttvZ/Hixc1ZAgAAAFQc/TQAAACVrGQh9dy5c3Prrbfm/vvvX2/bjBkzss8++6RPnz7ZZpttMmzYsMycObNUpQAAAEDF0E8DAADQ2pUspL755pvzxS9+MVOnTl1n/fLly3PkkUfmr3/9a4rFYorFYv74xz/msMMO8y1wAAAA2jz9NAAAAK1dyULqBx98MEny7//+7+us/+lPf5o5c+akV69euf7663P77bdnp512yty5c3P11VeXqhwAAACoCPppAAAAWruShdSzZs1KkgwePHid9b/61a9SKBRyySWX5OSTT87nP//5XH/99SkWi7n33ntLVQ4AAABUhJbQTy9YsCA333xzvvCFL2TIkCHp1q1bOnXqlJ122inHHHNM7rnnno3uO27cuBQKhUZfzz//fJPWDAAAQOXoUKqJ58+fn5qamnTp0qVhXX19ff70pz+lUCjk3/7t3xrW/8u//EvatWuX5557rlTlAAAAQEVoCf103759U1dX1/C+c+fOqaqqyty5czN37tz85je/yRFHHJG77747Xbt23eAcVVVV6dWr10aP0aFDyX4lAQAAQAtXsiup16xZk5UrV66z7umnn86yZcvywQ9+MNtss83/FtGuXbbZZpssXbq0VOUAAABARWgJ/XRdXV323XffXHPNNXnhhReyfPny1NbW5qWXXsrJJ5+cJJk8eXK+8pWvbHSOj3/845k3b95GXwMGDGjSmgEAAKgcJQupd9hhh6xcuTIvvfRSw7rf//73Sd5pVP+v2tra9/yGNQAAALQFLaGffuihh/L444/na1/7Wt7//vc3rB8wYEBuuOGGhnD69ttvz5w5c5r02AAAALR+JQup99tvvyTJBRdckPr6+rz55pu59tprUygUcvjhh68z9qWXXsrKlSuzww47lKocAAAAqAgtoZ8+5JBD3nP72qupk+Spp55q0mMDAADQ+pUspD7jjDOSJLfddlt69uyZnXfeOS+//HIGDhyYo446ap2x//Vf/5Uk+chHPlKqcgAAAKAiVEI/3blz54blNWvWlPXYAAAAVL4OpZp43333zU033ZRvfOMbWbJkSZJk8ODBufPOO9Ohw7qHvfXWW5M0/k1tAAAAaO0qoZ+eOnVqw/Iee+yxwTEzZszI0KFD88ILL6R9+/bp169fDjzwwJx66qnZa6+9tui4hUJhk8bNnz9/i+aH5rZrTX1zl7DV/PyVlnOExjhHaIxzhMY4RyiXkl1JnSQnnHBC5s2bl8cffzzPPfdcpk+fng996EPrjFm1alVOOeWU3HzzzRkxYkQpywEAAICK0JL76YULF+b73/9+kmTYsGEZNGjQBsfNnz8/zz77bLp27ZqVK1dm5syZueGGG7L33nvnO9/5TtnqBQAAoOUp2ZXUa3Xp0iUf/ehHN7q9Y8eOGTVqVKnLAAAAgIrSEvvp+vr6HH/88XnttdfSqVOnXHXVVeuN2W233TJ+/Ph8+tOfzsCBA1NVVZVVq1Zl6tSpOffcczNt2rRcfPHF2WabbTJmzJjNOn6xWHzP7WuvtO7Tp89mzQstxfOLS3o9SVn4+Sst5wiNcY7QGOcIjXGOUC6Vf6YBAAAAZXHGGWdk0qRJSZJrrrkme+6553pjRo4cmbPOOiu77757qqqqkrwTqA8fPjyPPvpoQ/A+bty4LFq0qHzFAwAA0GIIqQEAAIBGjR07NhMnTkySXHHFFTnppJM2e47OnTvnkksuSZLU1tZmypQpTVojAAAAlaHkIfXf//73nHLKKRkyZEhqamrSvn37jb46dCj53ccBAACgIrSkfvrss8/OD3/4wyTJZZddljPPPHOL59pvv/0all988cWtLQ0AAIAKVNIuduLEifnmN7+ZNWvWNPrcKAAAAOAdLamfPuusszJhwoQkyfjx4zN27NhmrQcAAIDKV7IrqR9//PGcccYZWbNmTU499dTcd999SZJevXrlwQcfzO23354TTzwxHTt2TJ8+fXLHHXfkoYceKlU5AAAAUBFaUj89duzYdQLqs846a6vn/POf/9ywPHDgwK2eDwAAgMpTsiupf/SjH6VYLObMM8/M5Zdf3rC+Y8eOOfTQQ5Mkn//85/ONb3wjhx9+eM4777z85S9/KVU5AAAAUBFaSj89duzYhlt8T5gwIWPGjGl0n2KxmEKhsNHtK1euzLe//e0kSbdu3fKJT3yiaYoFAACgopTsSuo//vGPKRQKOeOMM9ZZ/39vU/bhD384V111VV544YVcdtllpSoHAAAAKkJL6Ke/9a1vNQTUl19++SYF1Eny8MMP57DDDsvtt9+eV155pWH96tWrM2XKlAwbNiyPP/54kuT8889Pz549m7RuAAAAKkPJQurXX389nTp1yvve977/PVi7dlmxYsV6Yz/zmc+kqqoqv/rVr0pVDgAAAFSE5u6nZ8+enfHjxzcc9wc/+EH69u270dfa24En7wTpU6ZMyfHHH5+dd945Xbt2zbbbbptu3brlsMMOy5NPPpl27drl3HPPzdlnn91kNQMAAFBZSna7765du6aurm6ddd27d8/ixYuzcuXKdOrUqWF9VVVVunbtmpdffrlU5QAAAEBFaO5+ur6+fp3l119//T3H19bWNizvsccemTBhQh577LE8/fTTmT9/fhYuXJiuXbtmyJAhGTZsWE455ZTsscceTVYvAAAAladkIXW/fv0yY8aMLF68ODU1NUmSXXbZJX/961/z5JNP5oADDmgY++qrr2bRokXp2rVrqcoBAACAitDc/fSAAQPWu7X4purdu/cm3xocAACAtqtkt/v+0Ic+lCR57rnnGtYdfPDBKRaLufDCCxtuU7Zq1ap84xvfSBLfpAYAAKDN008DAADQ2pUspD7qqKNSLBbz85//vGHd17/+9XTq1ClTpkzJTjvtlP333z/9+vXLPffck0KhkNNOO61U5QAAAEBF0E8DAADQ2pUspD7yyCPz3e9+N7vttlvDuoEDB+aOO+5I9+7d89Zbb+Wxxx7LggULUigUcvbZZ2fkyJGlKgcAAAAqgn4aAACA1q5kz6SuqanJd7/73fXWf+Yzn8lBBx2U++67L3PmzEmPHj0yfPjw7LrrrqUqBQAAACqGfhoAAIDWrmQh9Xvp1atXvvCFLzTHoQEAAKBi6acBAABoDUp2u28AAAAAAAAA+L+E1AAAAAAAAACUTZPc7vukk05qimlSKBRy4403NslcAAAA0NLppwEAAGiLmiSkvuWWW1IoFFIsFrdo/7X7aqoBAABoS/TTAAAAtEVNElKPGjUqhUKhKaYCAACANkM/DQAAQFvUZFdSAwAAAJtHPw0AAEBb1K65CwAAAAAAAACg7RBSAwAAAAAAAFA2TXK778b84x//yO9///u8/PLLWb58eW688caGbatXr86bb76ZQqGQHXbYoRzlAAAAQEXQTwMAANAalTSkXrRoUU466aT8+te/TpIUi8UUCoX1muo999wzb7/9dmbOnJn3v//9pSwJAAAAWjz9NAAAAK1ZyW73XVdXlyOPPDK//vWv07Vr14wYMSKdO3deb1zXrl1z0kknpb6+vqH5BgAAgLZKPw0AAEBrV7KQ+sYbb8xjjz2W97///Xnuuedy7733pkePHhsc+9nPfjZJMnny5FKVAwAAABVBPw0AAEBrV7KQ+o477kihUMgVV1yRHXfc8T3H7rXXXmnXrl2eeeaZUpUDAAAAFUE/DQAAQGtXspD66aefTqFQyPDhwxsdW1VVlR49emTBggWlKgcAAAAqgn4aAACA1q5kIfXSpUvTvXv3dOzYcZPGr1q1Kh06dChVOQAAAFAR9NMAAAC0diULqfv06ZPFixdn6dKljY795z//maVLl2aHHXYoVTkAAABQEfTTAAAAtHYlC6k/+tGPJkl+97vfNTr2iiuuSJLsv//+pSoHAAAAKoJ+GgAAgNauZCH1iSeemGKxmPPPPz/z5s3b6LjLLrssP/7xj1MoFHLSSSeVqhwAAACoCPppAAAAWruSPbTqmGOOyYgRI/K73/0u++yzT77whS9kxYoVSZK77rorTz/9dO66667MnDkzSXLcccflwAMPLFU5AAAAUBH00wAAALR2JQupk+TOO+/MF77whfzmN7/JZZdd1rD+uOOOS5IUi8UkyWc+85nceOONpSwFAAAAKoZ+GgAAgNasZLf7TpJu3brlnnvuyW9/+9t86lOfSu/evVMsFlMsFlNTU5Mjjzwyv/71r/PLX/4ynTt3LmUpAAAAUDH00wAAALRmJb2Seq0RI0ZkxIgRSZK6urqsWbMmnTp1KsehAQAAoGLppwEAAGiNyhJSr3PADh3SoUPZDwsAAAAVTT8NAABAa1HS230DAAAAAAAAwLsJqQEAAAAAAAAoGyE1AAAAAAAAAGXTqkPqZcuWZfLkybnooovyr//6r3nf+96XQqGQQqGQcePGbdXc48aNa5jrvV7PP/9803wYAAAAAAAAgFagQ3MXUEpPPPFEjjzyyJIeo6qqKr169dro9g4dWvUfMQAAAAAAAMBmafUJ6jbbbJOPfOQjDa/Ro0dn3rx5TTb/xz/+8UydOrXJ5gMAAAAAAABozVp1SD1s2LC89dZb66w755xzmqkaAAAAAAAAAFr1M6nbt2/f3CUAAAAAAAAA8C5lC6mLxWLmz5+f2bNnl+uQAAAAUPH00wAAALQ2Jb/d91/+8pdcdNFFefDBB7N06dIUCoXU1dU1bH/77bdzzjnnpFAo5Ec/+lE6duxY6pKa1IwZMzJ06NC88MILad++ffr165cDDzwwp556avbaa6/Nnq9QKGzSuPnz52/23NAS7FpT39wlNAk/g6XVGs4T50hpOUdojHOExjhHqAStvZ8GAACg7SrpldS33XZb9ttvv/z6179ObW1tisViisXiOmO22WabvPTSS7n++uvzX//1X6UspyTmz5+fZ599Nl27ds3KlSszc+bM3HDDDdl7773zne98p7nLAwAAoAK1hX4aAACAtqtkV1I/++yz+fKXv5zVq1fnG9/4RkaNGpVPfvKTWbBgwXpjR40alQcffDC/+c1vMmLEiFKV1KR22223jB8/Pp/+9KczcODAVFVVZdWqVZk6dWrOPffcTJs2LRdffHG22WabjBkzZpPn/b+/dPi/1l5p3adPn62qH5rL84vL9pSBkvIzWFqt4TxxjpSWc4TGOEdojHOElqy199MAAABQspD68ssvz6pVq/L1r389V155ZZKkffv2Gxx76KGHJkkee+yxUpXT5EaOHLneuo4dO2b48OE58MADc+CBB+bJJ5/MuHHj8qUvfSk9evRohioBAACoNK29nwYAAICSXT7w0EMPpVAo5Fvf+lajY3fcccd07do1s2fPLlU5ZdW5c+dccsklSZLa2tpMmTKlmSsCAACgUrTlfhoAAIC2oWQh9auvvppu3bplp5122qTxXbp0yfLly0tVTtntt99+DcsvvvhiM1YCAABAJWnr/TQAAACtX8lC6k6dOmXVqlWNPmM5SZYvX56FCxe6JTYAAABtnn4aAACA1q5kIfWAAQOyevXq/POf/2x07H333Zc1a9ZkyJAhpSqn7P785z83LA8cOLAZKwEAAKCStPV+GgAAgNavZCH1Jz/5yRSLxfzHf/zHe45bsGBBzj777BQKhYwYMaJU5TSpxr7NvnLlynz7299OknTr1i2f+MQnylEWAAAArUBr7qcBAAAgKWFIPXr06FRXV+fHP/5xLrjggixZsmSd7cuXL88dd9yRffbZJy+99FJ69+6dr371q01ex9tvv5358+c3vOrr65Mky5YtW2d9bW3tOvuNGzcuhUIhhUIhs2bNWmfbww8/nMMOOyy33357XnnllYb1q1evzpQpUzJs2LA8/vjjSZLzzz8/PXv2bPLPBQAAQOvUUvppAAAAKJWShdTbb7997rjjjlRVVeXCCy/MtttumwULFiRJPvjBD6ZXr145/vjj8/LLL6dTp075z//8z9TU1DR5HXvttVe23XbbhtecOXOSJJdddtk660877bRNnrNYLGbKlCk5/vjjs/POO6dr167Zdttt061btxx22GF58skn065du5x77rk5++yzm/wzAQAA0Hq1lH4aAAAASqVkIXWSHHXUUXn44Yez9957Z9WqVamrq0uxWMyzzz6blStXplgsZq+99srDDz9cUbfE3mOPPTJhwoR89rOfze67754uXbpk4cKF6dKlS/bcc8+cdtpp+dvf/paLL764uUsFAACgArXWfhoAAACSpEOpD7DvvvvmiSeeyD/+8Y88+uijefXVV7NmzZr07ds3+++/f/bZZ5+SHv//3qp7U40bNy7jxo3b4LbevXtnzJgxW14UAAAANKK5+2kAAAAolZKH1Gt96EMfyoc+9KFyHQ4AAABaBf00AAAArU3Jbvd966235q677trk8b/61a9y6623lqocAAAAqAj6aQAAAFq7koXUJ554Ys4888xNHj9mzJicdNJJpSoHAAAAKoJ+GgAAgNauZCF1khSLxZKOBwAAgNZIPw0AAEBrVtKQenMsXrw4HTt2bO4yAAAAoKLopwEAAKg0LSKkfuyxx/L2229nxx13bO5SAAAAoGKUqp9esGBBbr755nzhC1/IkCFD0q1bt3Tq1Ck77bRTjjnmmNxzzz2NzvH6669nzJgxGTRoULp06ZJevXpl2LBhueGGG1z5DQAA0MZ1aKqJfvrTn+anP/3pOuveeuutHHrooRvdp1gsZuHChZkxY0YKhUIOO+ywpioHAAAAKkJL7Kf79u2burq6hvedO3dOVVVV5s6dm7lz5+Y3v/lNjjjiiNx9993p2rXrevtPmzYthx9+eBYsWJAkqa6uzpIlS/Loo4/m0UcfzV133ZV77703nTp1atK6AQAAqAxNFlLPmjUrU6dOXWfdqlWr1lu3MYMGDcq4ceOaqhwAAACoCC2xn66rq8u+++6bE088MYcffnje//73N9R60UUX5cYbb8zkyZPzla98Jbfddts6+y5atChHHXVUFixYkMGDB+e2227LPvvsk1WrVuX666/P6NGj88ADD2T06NG55pprmrRuAAAAKkOThdQHH3zwOu8vuOCCVFdXZ8yYMRvdp127dqmpqcnQoUNz8MEHp3379k1VDgAAAFSElthPP/TQQznkkEPWWz9gwIDccMMN6dChQ37yk5/k9ttvzyWXXJKdd965YcyECRMyb968dOnSJffdd18GDhyYJOnYsWO+/vWvZ/HixTn33HNz3XXX5cwzz8zuu+/epLUDAADQ8jVZSH3QQQfloIMOani/tqn+7ne/21SHAAAAgFanJfbTGwqo3+3kk0/OT37ykyTJU089tU5IfeuttyZJjjvuuIaA+t1OP/30XHLJJamtrc3PfvazXHDBBU1YOQAAAJWgXakmfumll/LEE0+UanoAAABolSqhn+7cuXPD8po1axqWn3vuucyePTtJcsQRR2xw3+rq6gwbNixJ8sADD5SwSgAAAFqqJruS+v963/veV6qpAQAAoNWqhH763c/L3mOPPRqWp0+f3rA8dOjQje4/dOjQTJ48Oc8888xmHbdQKGzSuPnz52/WvNBS7FpT39wlbDU/f6XlHKExzhEa4xyhMc4RyqVkV1IDAAAArc/ChQvz/e9/P0kybNiwDBo0qGHbq6++2rDcr1+/jc6xdtvixYtTW1tbokoBAABoqUp2JXX79u03e59CoZC6uroSVAMAAACVoSX30/X19Tn++OPz2muvpVOnTrnqqqvW2b5kyZKG5a5du250nndvW7JkSaqrqzfp+MVi8T23r73Suk+fPps0H7Q0zy+u/OtJ/PyVlnOExjhHaIxzhMY4RyiXkoXUjTWOAAAAwPpacj99xhlnZNKkSUmSa665JnvuuWczVwQAAEAlKllI/Yc//OE9ty9atCiPP/54rr/++hSLxVx99dXZfvvtS1UOAAAAVISW2k+PHTs2EydOTJJcccUVOemkk9Yb071794blZcuWpaamZoNzLVu2bIP7AAAA0DaULKQ+6KCDGh3zqU99KmeccUYOOeSQfPe7381TTz1VqnIAAACgIrTEfvrss8/OD3/4wyTJZZddljPPPHOD43bccceG5blz5240pJ47d26SpKamZpNv9Q0AAEDr0ew3lt9uu+1y9dVX57nnnsv3v//95i4HAAAAKkK5+umzzjorl112WZJk/PjxGTt27EbHDh06tGF5+vTpGx23dtuQIUOaqEoAAAAqSbOH1Mk73xLv3Llz7r777uYuBQAAACpGqfvpsWPHZsKECUneCajPOuus9xw/aNCg9O/fP0ly//33b3DM0qVL88gjjyRJhg8f3oTVAgAAUClaREhdKBTSrl27zJ49u7lLAQAAgIpRyn567NixDbf4njBhQqMB9VqjRo1Kktx5552ZNWvWetuvvvrq1NbWpn379hk5cmST1QsAAEDlaBEh9bRp07Js2bJ07dq1uUsBAACAilGqfvpb3/pWQ0B9+eWXZ8yYMZu879ixY9O3b98sW7YsI0aMyLRp05Ikq1atyrXXXpvzzjsvSXLKKadk9913b9K6AQAAqAzNHlI/8cQTOf7441MoFLL//vs3dzkAAABQEUrVT8+ePTvjx49PkrRr1y4/+MEP0rdv342+1t4OfK0ePXpk0qRJ6d27d5555pnss88+qampSXV1dU499dSsWrUqw4cPzxVXXNFkNQMAAFBZOpRq4kMPPfQ9t69YsSJz5szJq6++mmKxmI4dO+Y73/lOqcoBAACAitDc/XR9ff06y6+//vp7jq+trV1v3d57750ZM2bkBz/4QSZNmpQ5c+akW7duGTp0aE444YScdNJJadeu2b83DwAAQDMpWUg9derUTR77vve9Lz/5yU/y0Y9+tFTlAAAAQEVo7n56wIABKRaLWz3P9ttvn8svvzyXX355E1QFAABAa1KykPq73/3uex+4Q4dss8022XPPPfPxj388hUKhVKUAAABAxdBPAwAA0No1W0gNAAAArE8/DQAAQGvnAVAAAAAAAAAAlI2QGgAAAAAAAICyEVIDAAAAAAAAUDZN8kzq9u3bN8U0KRQKqaura5K5AAAAoKXTTwMAANAWNUlIXSwWm2IaAAAAaFP00wAAALRFTRJS/+EPf2iKaQAAAKBN0U8DAADQFjVJSH3QQQc1xTQAAADQpuinAQAAaIvaNXcBAAAAAAAAALQdQmoAAAAAAAAAyqbkIXWxWMyvfvWrHHvssRk4cGC6deuW6urqDBw4MJ/73Ofy61//OsVisdRlAAAAQEXRTwMAANBaNckzqTfm9ddfz7/927/lT3/6U5Ks0zy//PLLmT17dn75y19m//33zy9+8Yv07du3lOUAAABARdBPAwAA0JqVLKRetWpVDj/88Dz99NMpFovZd9998y//8i/ZaaedkiSvvPJKHnzwwTz++OP54x//mCOOOCKPP/54OnbsWKqSAAAAoMXTTwMAANDalSykvvbaa/OPf/wjNTU1uf3223PUUUetN+Z73/te7rvvvnz+85/PP/7xj/z4xz/ON77xjVKVBAAAAC2efhoAAIDWrmTPpP7FL36RQqGQq6++eoMN9VpHHnlkrr766hSLxfz85z8vVTkAAABQEfTTAAAAtHYlC6mfffbZVFVV5d///d8bHfvv//7v6dixY5599tlSlQMAAAAVQT8NAABAa1eykHr58uXp2rVrOnRo/I7iHTp0SNeuXbN8+fJSlQMAAAAVQT8NAABAa1eykHr77bfPokWLMnv27EbHzpo1KwsXLsz2229fqnIAAACgIuinAQAAaO1KFlIfeOCBKRaLGT16dIrF4kbHFYvFfPOb30yhUMhBBx1UqnIAAACgIuinAQAAaO1KFlKvbZR//etf55BDDsmUKVOyevXqhu2rV6/Ogw8+mEMOOSS//vWvUygUMnr06FKVAwAAABVBPw0AAEBr1/gDrrbQhz/84fzwhz/MN7/5zTzyyCMZPnx4OnTokD59+qRQKOTNN99MXV1dw7fCJ0yYkA9/+MOlKgcAAAAqgn4aAACA1q5kV1InyZlnnpl77703gwcPTrFYzOrVq/Paa6/l1VdfzerVq1MsFjNkyJD89re/zZlnnlnKUgAAAKBi6KcBAABozUp2JfVaRx11VI466qg8/fTTeeqpp/LGG28kSbbbbrvss88+2WOPPUpdAgAAAFQc/TQAAACtVclD6rX22GMPDTQAAABsJv00AAAArU1Jb/cNAAAAAAAAAO9WspB61apVmT17dubNm7fettra2owdOzZ77rlnPvKRj+S8887L8uXLS1UKAAAAVAz9NAAAAK1dyW73fcMNN+T000/PCSeckJtuummdbSNGjMijjz6aYrGYJPn73/+eRx55JH/4wx9SKBRKVRIAAEDrMa5Hc1ew9cYtau4KWiT9NAAAAK1dya6k/v3vf58k+fznP7/O+nvvvTePPPJICoVCRo4cmS996UupqqrKI488kttuu61U5QAAAEBF0E8DAADQ2pUspH722WeTJHvvvfc66++4444UCoV861vfym233ZbrrrsuV155ZYrFYu64445SlQMAAAAVQT8NAABAa1eykPrNN99M165ds80226yz/g9/+EOS5Etf+lLDuuOPPz7JO7cpAwAAgLZMPw0AAEBrV7KQeunSpWnXbt3pZ82alTfffDM777xzBg4c2LC+W7du6dmzZ956661SlQMAAAAVQT8NAABAa1eykLpXr16pra3NwoULG9Y99NBDSZKPf/zj642vq6tLdXV1qcoBAACAiqCfBgAAoLUrWUj9kY98JEly4403Jknq6+tz4403plAo5JBDDlln7Jtvvpna2tr07du3VOUAAABARdBPAwAA0NqVLKQ+8cQTUywWc8455+SII47Ivvvum8ceeyzV1dU59thj1xn7yCOPJEk+8IEPlKocAAAAqAj6aQAAAFq7koXUn/vc5/LFL34xa9asye9///v85S9/SefOnfPjH/84PXv2XGfsz3/+8w1+IxwAAADaGv00AAAArV2HUk5+44035qSTTsqf/vSn9OzZM4cddlgGDhy4zphVq1alR48eGTVqVI488shSlgMAAAAVQT8NAABAa1bSkDpJ9t9//+y///4b3d6xY8dcd911pS4DAAAAKop+GgAAgNaqZLf7BgAAAAAAAID/q+RXUifv3ILsv/7rv/LUU0/ljTfeSJJst9122WefffIv//Iv6dixYznKAAAAgIqinwYAAKA1KnlIPXHixFxwwQV56623Nri9V69eOf/883P66aeXuhQAAACoGPppAAAAWquShtRf+tKXcvPNN6dYLCZJdtppp/Tr1y9JMnfu3LzyyitZsGBBzjzzzPz1r3/NTTfdVMpyAAAAoCLopwEAAGjNSvZM6jvuuCM33XRTisVivvCFL2TmzJmZPXt2HnvssTz22GOZPXt2/vnPf2bUqFEpFov56U9/mjvuuKNU5QAAAEBF0E8DAADQ2pUspL722mtTKBRy+umn59Zbb82uu+663phddtklt9xyS04//fQUi8Vcc801pSoHAAAAKoJ+GgAAgNauZCH1P/7xjxQKhZx//vmNjj3//PNTKBTy9NNPl6ocAAAAqAj6aQAAAFq7koXUhUIhPXv2TO/evRsd27t37/Ts2TOFQqFU5QAAAEBF0E8DAADQ2pUspN59992zaNGi1NbWNjq2trY2ixcvzqBBg0pVDgAAAFQE/TQAAACtXclC6pNOOilr1qzJVVdd1ejYiRMnZs2aNTn55JNLVQ4AAABUBP00AAAArV2HUk381a9+Nf/93/+d8847L6tWrcqYMWNSXV29zpilS5fmhz/8Yb73ve/luOOOyymnnFKqcgAAoLKM69HcFWy9cYuauwKoSPppAAAAWrsmCalPOumkDa7v0qVLunfvngsvvDCXXXZZ9tlnn/Tr1y+FQiGvvPJKnnrqqSxfvjw9evRI586dc/LJJ+fGG29sipIAAACgxdNPAwAA0BY1SUh9yy23pFAopFgsrrP+3euWLVuWhx9+eIP7L1y4sGEOTTUAAABthX4aAACAtqhJQupRo0alUCg0xVQAAADQZuinAQAAaIua7EpqAAAAYPPopwEAAGiL2jV3AQAAAAAAAAC0HS0ipK6vr89vf/vbHHPMMc1dCgAAAFQM/TQAAACVqFlD6pkzZ+acc87JTjvtlGOOOSa//e1vm7McAAAAqAil7qeXLVuWyZMn56KLLsq//uu/5n3ve18KhUIKhULGjRv3nvuOGzeuYex7vZ5//vkmrRkAAIDK0STPpN4cy5Ytyy9+8YvcdNNN+eMf/5gkKRaLSZIPfOAD5S4HAAAAKkI5++knnngiRx555FbNUVVVlV69em10e4cOZf+VBAAAAC1E2TrCP//5z7npppvy85//PLW1tUneaaYHDx6cY489Nscee2yGDh1arnIAAACgIjRXP73NNtvkIx/5SMNr9OjRmTdv3ibv//GPfzxTp05t8roAAACofCUNqd98883cdtttufHGG/M///M/Sf73W96FQiFPPvlk9t5771KWAAAAABWnufvpYcOG5a233lpn3TnnnFOy4wEAANC2NHlIXSwWM3ny5Nx4442ZNGlS6urqUiwW06VLlxxzzDE54YQT8slPfjKJ23sDAADAWi2pn27fvn1J5wcAAKBta7KQ+oUXXshNN92Un/70p3nttddSLBZTKBRywAEHZNSoUfnc5z6X7t27N9XhAAAAoFXQTwMAANDWNFlIvdtuu6VQKKRYLOb9739/jj/++IwaNSoDBw5sqkMAAABAq9Na++kZM2Zk6NCheeGFF9K+ffv069cvBx54YE499dTstddeWzRnoVDYpHHz58/fovmhue1aU9/cJWw1P3+l5RyhMc4RGuMcoTHOEcqlyW/3/Y1vfCPjx49Px44dm3pqAAAAaLVaWz89f/78vPXWW+nZs2cWL16cmTNnZubMmbnxxhtz7rnn5qKLLmruEgEAAGgmTRZSd+zYMatWrcpVV12V22+/Pf/+7/+e448/Ph/72Mea6hCbbdmyZfnv//7vTJs2LX/5y18ybdq0zJ49O0ny3e9+N+PGjdvqY7z++usZP358Jk2alNmzZ6dLly754Ac/mBNOOCEnn3zyJn/LGwAAgLapJfbTW2O33XbL+PHj8+lPfzoDBw5MVVVVVq1alalTp+bcc8/NtGnTcvHFF2ebbbbJmDFjNmvuYrH4ntvX9uB9+vTZ4vqhOT2/uF1zl7DV/PyVlnOExjhHaIxzhMY4RyiXJjvT5s2blx/96Ef50Ic+lLfeeivXXntt9t9//wwaNCiXXHJJQzhcTk888USOPPLInHfeebnnnnuavIZp06blgx/8YC6//PLMnDkzHTp0yJIlS/Loo4/my1/+cj75yU9m5cqVTXpMAAAAWpeW2E9vjZEjR+ass87K7rvvnqqqqiTvBPHDhw/Po48+mo9+9KNJknHjxmXRokXNWSoAAADNpMlC6p49e+a0007LX//610ybNi1f+9rX0qNHj/zzn//Meeedl/e///059NBDc/PNNzfVITfJNttsk0984hM566yz8p//+Z/p27dvk8y7aNGiHHXUUVmwYEEGDx6cJ598MkuWLMnSpUszceLEVFVV5YEHHsjo0aOb5HgAAAC0Ti21ny6Fzp0755JLLkmS1NbWZsqUKc1cEQAAAM2hJNfs77XXXrn66qvz2muv5bbbbstBBx2UYrGYqVOn5ktf+lLDuAceeCB1dXWlKCFJMmzYsLz11lt58MEHM378+Bx33HHp1KlTk8w9YcKEzJs3L126dMl9992XffbZJ8k73w7/+te/ngsuuCBJct1112XmzJlNckwAAABat5bST5fSfvvt17D84osvNmMlAAAANJeS3li+U6dOGTlyZB566KE8//zzOffcc9OvX78k7zxH6rOf/Wy22267fPGLX8x9993X5A12+/btm3S+d7v11luTJMcdd1wGDhy43vbTTz891dXVWbNmTX72s5+VrA4AAABan+bupwEAAKCUyvb084EDB+aiiy7Kyy+/nPvuuy//+q//mg4dOmThwoW59dZbc/TRR2f77bcvVzlb5bnnnmt4JtgRRxyxwTHV1dUZNmxYkne+4Q4AAABbojX100ny5z//uWF5Q1/6BgAAoPUrW0i9VqFQyCc/+cncfffdmTt3biZMmJAhQ4akWCxm4cKF5S5ni0yfPr1heejQoRsdt3bbM888U/KaAAAAaN0qoZ8uFovvuX3lypX59re/nSTp1q1bPvGJT5SjLAAAAFqYDs158D59+uSb3/xmvvnNb+bPf/5zbrrppuYsZ5O9+uqrDctrb7e2IWu3LV68OLW1tamurm507kKhsEk1zJ8/f5PGQUuza019c5fQJPwMllZrOE+cI6XlHKExreIc6Ti4uUvYei34PHeOtBAt+Bxp6crRT7/99ttZs2ZNw/v6+nd+bpYtW7bO/4917ty5od99+OGH873vfS8nnnhiDj744Oy0005JktWrV+fhhx/O//f//X958sknkyTnn39+evbs2eR1AwAA0PI1a0j9bh/72MfysY99rLnL2CRLlixpWO7atetGx71725IlSzYppAYAAIDNUap+eq+99srLL7+83vrLLrssl112WcP7E044IbfcckuSd66knjJlSqZMmZIk6dKlS7p165ZFixZl9erVSZJ27drlnHPOydlnn93kNQMAAFAZWkxIzTsauzXa2iut+/TpU45yoMk9v7jsTxkoCT+DpdUazhPnSGk5R2hMqzhHOv9Pc5ew9Vrwee4caSFa8DnCltljjz0yYcKEPPbYY3n66aczf/78LFy4MF27ds2QIUMybNiwnHLKKdljjz2au1QAAACakZB6C3Tv3r1hedmyZampqdnguGXLlm1wHwAAAGjpZs2atdn79O7dO2PGjGn6YgAAAGhVKv/ygWaw4447NizPnTt3o+PWbqupqXGrbwAAAAAAAIAIqbfI0KFDG5anT5++0XFrtw0ZMqTkNQEAAAAAAABUAiH1Fhg0aFD69++fJLn//vs3OGbp0qV55JFHkiTDhw8vW20AAAAAAAAALZmQeguNGjUqSXLnnXdu8DldV199dWpra9O+ffuMHDmyzNUBAAAAAAAAtEytPqR+++23M3/+/IZXfX19kmTZsmXrrK+trV1nv3HjxqVQKKRQKGwwhB47dmz69u2bZcuWZcSIEZk2bVqSZNWqVbn22mtz3nnnJUlOOeWU7L777qX9kAAAAAAAAAAVotWH1HvttVe23XbbhtecOXOSJJdddtk660877bTNmrdHjx6ZNGlSevfunWeeeSb77LNPampqUl1dnVNPPTWrVq3K8OHDc8UVV5TiYwEAAAAAAABUpFYfUpfS3nvvnRkzZmT06NHZbbfdsnr16nTr1i0HHHBArr/++kyePDmdOnVq7jIBAAAAAAAAWowOzV1AqW3oVt2bYty4cRk3blyj47bffvtcfvnlufzyy7foOAAAAAAAAABtiSupAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAKBshNQAAAAAAAAAlI2QGgAAAAAAAICyEVIDAAAAAAAAUDZCagAAAAAAAADKRkgNAAAAAAAAQNkIqQEAAAAAAAAoGyE1AAAAAAAAAGUjpAYAAAAAAACgbITUAAAAAAAAAJSNkBoAAAAAAACAshFSAwAAAAAAAFA2QmoAAAAAAAAAykZIDQAAAAAAAEDZCKkBAAAAAAAAKBshNQAAAAAAAABlI6QGAAAAAAAAoGyE1AAAAAAAAACUjZAaAAAAAAAAgLIRUgMAAAAAAABQNkJqAAAAAAAAAMpGSA0AAAAAAABA2QipAQAAAAAAACgbITUAAAAAAAAAZSOkBgAAAAAAAPj/27v/MKvLOn/8zxl+DIwzIIOGUhpsIkqMLZG2ariShD/Q2mq7anOXVte0rdxyzd2yLQe3tt1C90ptzczNvOLKS21bbQ1ByU2xYoUwxB9Q4ZhBsAIKAwNIzPn+4Xdm8cOPGWTmfWbOPB7XNdd1mPs+97zO4eZwXvM87/ebwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAAAAAACiMkBoAAAAAAACAwgipAQAAgFdobW3N3Llz84UvfCHvec978vrXvz5VVVWpqqpKU1NTl9ZYt25dLr/88owfPz5Dhw5NQ0NDpkyZkm9+85splUo9+wAAAADo1QaWuwAAAACgd/mf//mfnHPOOa/6/kuWLMmZZ56ZDRs2JEnq6urS0tKShQsXZuHChbnzzjtzzz33pKamprtKBgAAoA9xJDUAAACwhxEjRuSMM87IFVdcke9+97s54ogjunS/TZs25dxzz82GDRty3HHH5dFHH01LS0u2bt2aG264IYMGDcr8+fNz2WWX9fAjAAAAoLdyJDUAAADwClOmTMnGjRtf8b1Pf/rTXbrv7Nmzs3bt2gwdOjQ//OEPM3bs2CTJ4MGD87GPfSybN2/OlVdemW984xv55Cc/mWOPPbbb6wcAAKB3cyQ1AAAA8AoDBgx41fe97bbbkiQf+MAHOgLq3V166aWpq6vLrl27MmfOnFf9cwAAAOi7HEkNAAAAdIsVK1bkN7/5TZLk7LPP3uucurq6TJkyJXPnzs38+fMza9asIksEAAAqXdPwclfQPZo2lbuCHiWkBgAAALrF8uXLO25PnDhxn/MmTpyYuXPn5sknnzyg9auqqro0b/369Qe0LvQWxwxrK3cJB82/v55lj9AZe4TO2CN0piL2yODjyl1C96jwve503wAAAEC3WLNmTcft1772tfuc1z62efPmbNmypcfrAgAAoHdxJDUAAADQLVpaWjpu19bW7nPe7mMtLS2pq6vr0vqlUmm/4+1HWh922GFdWg96m19t7vvHk/j317PsETpjj9AZe4TOVMQeGfJ0uUvoHhW+1/v+TgMAAAAAAACgzxBSAwAAAN2ivr6+43Zra+s+5+0+tvt9AAAA6B+E1AAAAEC3GD16dMft1atX73Ne+9iwYcO6fKpvAAAAKoeQGgAAAOgWEydO7Li9fPnyfc5rH5swYUKP1wQAAEDvI6QGAAAAusX48eNz9NFHJ0nuu+++vc7ZunVrHn744STJ9OnTC6sNAACA3kNIDQAAAHSbmTNnJkluv/32NDc37zH+ta99LVu2bMmAAQNy/vnnF1wdAAAAvYGQGgAAANjDCy+8kPXr13d8tbW1JUlaW1tf8f0tW7a84n6f+tSncsQRR6S1tTUzZszIkiVLkiQvvfRSbrzxxnzuc59Lklx88cU59thji31QAAAA9ApCagAAAGAPkyZNyuGHH97x9dxzzyVJvvKVr7zi+x//+Mdfcb/hw4fnv/7rvzJy5Mg8+eSTectb3pJhw4alrq4uH/3oR/PSSy9l+vTp+dd//ddyPCwAAAB6ASE1AAAA0K0mT56cJ554IpdddlnGjRuXnTt35pBDDsnb3va23HzzzZk7d25qamrKXSYAAABlMrDcBQAAAAC9z96uJ30gRo0alWuvvTbXXntt9xQEQGVoGl7uCrpH06ZyVwAAfZqQGgAAAAAA6B18kAGgXxBSAwAAAADdQ7gEAEAXuCY1AAAAAAAAAIURUgMAAAAAAABQGCE1AAAAAAAAAIURUgMAAAAAAABQGCE1AAAAAAAAAIURUgMAAAAAAABQmIHlLgAAAAAAAAC6pGl4uSvoHk2byl0BlJUjqQEAAAAAAAAojJAaAAAAAAAAgMIIqQEAAAAAAAAojJAaAAAAAAAAgMIIqQEAAAAAAAAojJAaAAAAAAAAgMIIqQEAAAAAAAAojJAaAAAAAAAAgMIIqQEAAAAAAAAojJAaAAAAAAAAgMIIqQEAAAAAAAAojJAaAAAAAAAAgMIIqQEAAAAAAAAoTL8IqVtaWtLU1JTGxsbU1dVl+PDhOfHEE3PNNdfkpZdeelVrNjU1paqqqtOvX/3qV938aAAAAAAAAAD6roHlLqCnPfvsszn99NPT3NycJKmtrc2OHTuyePHiLF68OHPmzMmCBQsyYsSIV7X+oEGD0tDQsM/xgQMr/ikGAAAAAAAA6LKKPpJ6165dOe+889Lc3Jwjjzwy999/f7Zu3ZrW1tbcfvvtqa+vz9KlS3P++ee/6p9xyimnZO3atfv8GjNmTPc9IAAAAAAAAIA+rqJD6ltvvTWPP/54kuR73/tepk2bliSprq7O+9///tx0001Jkrlz52bBggVlqxMAAAAAAACgv6jokPrb3/52kmTq1Kk5+eST9xj/wAc+kLFjxyZJbrvttkJrAwAAAAAAAOiPKjakbm1tzSOPPJIkOfvss/c6p6qqKmeddVaSZP78+YXVBgAAAAAAANBfVWxI/dRTT6WtrS1JMnHixH3Oax9bu3ZtNm7ceMA/54knnsjEiRMzdOjQ1NXVZfz48fnwhz+cpUuXvrrCAQAAAAAAACrYwHIX0FPWrFnTcfu1r33tPuftPrZmzZo0NDQc0M9Zv359Nm7cmEMPPTSbN2/OypUrs3Llytxyyy258sor84UvfOGA1quqquryz4W+6JhhbeUuoVv4N9izKmGf2CM9yx6hMxWxRwYfV+4SDl4v3uf2SC/Ri/cIAAAA0HMq9kjqlpaWjtu1tbX7nLf72O736cy4cePy5S9/OStWrMj27duzYcOGbN26NfPmzcvkyZNTKpXyxS9+Mddcc82rewAAAAAAAAAAFahij6Tuaeeff/4e3xs8eHCmT5+e0047LaeddloeffTRNDU15aKLLsrw4cO7tG6pVNrvePuR1ocddtiBFw29wK82V8ZnY/wb7FmVsE/skZ5lj9CZitgjQ54udwkHrxfvc3ukl+jFewQAAADoOX3/NzP7UF9f33G7tbV1n/N2H9v9PgdjyJAh+ad/+qckyZYtW7JgwYJuWRcAAAAAAACgr6vYkHr06NEdt1evXr3PebuP7X6fg3XyySd33F61alW3rQsAAAAAAADQl1VsSH388cenuvrlh7d8+fJ9zmsfO+KII9LQ0FBIbQAAAAAAAAD9VcWG1LW1tTn11FOTJPfdd99e55RKpcybNy9JMn369G79+T/72c86bo8dO7Zb1wYAAAAAAADoqyo2pE6SD33oQ0mSBx98MIsWLdpj/M477+w4FffMmTO7vG6pVNrv+I4dO/LZz342SXLIIYfkjDPO6PLaAAAAAAAAAJWs4kPqxsbGlEqlvPe9782CBQuSJG1tbbnzzjvz4Q9/OEly9tln7xEkNzU1paqqKlVVVWlubn7F2EMPPZRp06blO9/5Tn772992fH/nzp1ZsGBBpkyZ0hGKf/7zn8+hhx7acw8SAAAAAAAAoA8ZWO4CetLAgQNzzz33ZOrUqWlubs60adNSW1ubtra2bN++PUkyadKkzJkz54DWLZVKWbBgQUfoPXTo0BxyyCHZtGlTdu7cmSSprq7Opz/96fzd3/1d9z4oAAAAAAAAgD6sokPqJBkzZkyWLVuW2bNn5z/+4z/yzDPPZNCgQXnjG9+YP/uzP8ull16awYMHH9CajY2NmT17dn7605/m8ccfz/r16/Piiy+mtrY2EyZMyJQpU3LxxRensbGxhx4VAAAAAAAAQN9U8SF1ktTX12fWrFmZNWtWl+/T1NSUpqamvY6NHDkyl19+eTdVBwAAAAAAANB/VPQ1qQEAAAAAAADoXYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAB0q1tvvTVVVVWdfj3wwAPlLhUAAIAyGFjuAgAAAIDKVF1dncMPP3yf4zU1NQVWAwAAQG8hpAYAAAB6xFFHHZXm5uZylwEAAEAv43TfAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABTGNakBAACAHvH8889n8uTJWbFiRXbt2pUjjzwyp5xySi666KKcfvrpB7xeVVVVl+atX7/+gNeG3uCYYW3lLuGgrR98XLlL6B699HXEHulF7JEeY4/0LHukF7FHeow90jc4khoAAADoEa2trfn5z3+ewYMHp62tLc8880zmzJmTqVOn5sILL8zvf//7cpcIAABAGTiSGgAAAOhWo0ePzlVXXZX3vOc9GT9+fGpqarJr164sWrQoV111VR544IF861vfyiGHHJLrr7++y+uWSqX9jrcfaX3YYYcdVP1QLr/a3PePJzlsyNPlLqF79NLXEXukF7FHeow90rPskV7EHukx9kjf0Pd3GgAAANCrTJ8+PU1NTTnhhBNSU1OTJBkwYEBOOeWUzJs3L+9617uSJP/2b/+WX/7yl+UsFQAAgDIQUgMAAACFqa6uzuzZs5MkbW1t+cEPflDmigAAACiakBoAAAAo1DHHHNNxSu5Vq1aVuRoAAACKJqQGAAAAAAAAoDBCagAAAKBQv/71r7N+/fokydixY8tcDQAAAEUTUgMAAADdplQqdTp+xRVXJHn5+tTnnntuEWUBAADQiwipAQAAgG7z7LPP5qSTTspNN92UVatWdYTWbW1t+dnPfpazzz473//+95Mkl1xyScaPH1/OcgEAACiDgeUuAAAAAKgsjz76aB599NEkSU1NTerr69PS0pIdO3Z0zLngggty3XXXlatEAAAAykhIDQAAAHSbUaNG5frrr89Pf/rTPPbYY3n++efzwgsvZMiQIRk7dmxOOeWUXHjhhTn11FPLXSoAAABlIqQGAAAAus3QoUPz8Y9/PB//+MfLXQoAAAC9lGtSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFAYITUAAAAAAAAAhRFSAwAAAAAAAFCYgeUuAKBPahpe7goOXtOmclcAAAAAAAD0Q0JqAIBy8GEXAAAAAKCfcrpvAAAAAAAAAAojpAYAAAAAAACgME73DQA9wamcAQAAAABgrxxJDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhBpa7ACrLmE/fW+4SukXzP88odwkAAAAAAABQkRxJDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFEZIDQAAAAAAAEBhhNQAAAAAAAAAFGZguQsAAAAA6A/GfPrecpdw0Jr/eUa5SwAAACqAkBr2pml4uSs4eE2byl0BAAAAAAAA7MHpvgEAAAAAAAAoTL8IqVtaWtLU1JTGxsbU1dVl+PDhOfHEE3PNNdfkpZdeOqi1161bl8svvzzjx4/P0KFD09DQkClTpuSb3/xmSqVSNz0CAAAA6Ht6sh8HAACg76r4030/++yzOf3009Pc3Jwkqa2tzY4dO7J48eIsXrw4c+bMyYIFCzJixIgDXnvJkiU588wzs2HDhiRJXV1dWlpasnDhwixcuDB33nln7rnnntTU1HTnQwIAAIBeryf7cQAAAPq2ij6SeteuXTnvvPPS3NycI488Mvfff3+2bt2a1tbW3H777amvr8/SpUtz/vnnH/DamzZtyrnnnpsNGzbkuOOOy6OPPpqWlpZs3bo1N9xwQwYNGpT58+fnsssu64FHBgAAAL1XT/bjAAAA9H0VHVLfeuutefzxx5Mk3/ve9zJt2rQkSXV1dd7//vfnpptuSpLMnTs3CxYsOKC1Z8+enbVr12bo0KH54Q9/mLe85S1JksGDB+djH/tYZs2alST5xje+kZUrV3bXQwIAAIBeryf7cQAAAPq+ig6pv/3tbydJpk6dmpNPPnmP8Q984AMZO3ZskuS22247oLXb5+++xu4uvfTS1NXVZdeuXZkzZ86Blg4AAAB9Vk/24wAAAPR9FXtN6tbW1jzyyCNJkrPPPnuvc6qqqnLWWWflxhtvzPz587u89ooVK/Kb3/xmv2vX1dVlypQpmTt3bubPn99xZDUAAABUsp7sx+kFmoaXu4Lu0bSp3BUAAEC/VrFHUj/11FNpa2tLkkycOHGf89rH1q5dm40bN3Zp7eXLl+9x//2t/eSTT3ZpXQAAAOjrerIfBwAAoDJU7JHUa9as6bj92te+dp/zdh9bs2ZNGhoaun3tzZs3Z8uWLamrq+t07aqqqk7nHMg8Xp2KeHZnVcSj6LUq4tm1R3pURTy79kiPqohn1x7pURXx7NojPaoinl17pCL1ZD+uZy6/inlmvf70mIp5Zu2RHlMxz6w90mMq5pm1R3pMxTyz9kiPqZhntsL3SMUeSd3S0tJxu7a2dp/zdh/b/T7lWhsAAAD6Mj0zAAAAnanYI6n7qlKpVO4S+r32T9z7u2Bf7BE6Y4/QGXuEztgjdMYeob+y58vLaw+dsUfojD1CZ+wROmOP0Bl7pO+o2COp6+vrO263trbuc97uY7vfp1xrAwAAQF+mZwYAAKAzFRtSjx49uuP26tWr9zlv97Hd79Odaw8bNqxL16MGAACAvq4n+3EAAAAqQ8WG1Mcff3yqq19+eMuXL9/nvPaxI444Ig0NDV1ae+LEiXvcf39rT5gwoUvrAgAAQF/Xk/04AAAAlaFiQ+ra2tqceuqpSZL77rtvr3NKpVLmzZuXJJk+fXqX1x4/fnyOPvro/a69devWPPzwwwe8NgAAAPRlPdmPAwAAUBkqNqROkg996ENJkgcffDCLFi3aY/zOO+/MqlWrkiQzZ848oLXb599+++1pbm7eY/xrX/tatmzZkgEDBuT8888/wMoBAACg7+rJfhwAAIC+r+JD6sbGxpRKpbz3ve/NggULkiRtbW2588478+EPfzhJcvbZZ+eMM854xX2bmppSVVWVqqqqvYbQn/rUp3LEEUektbU1M2bMyJIlS5IkL730Um688cZ87nOfS5JcfPHFOfbYY3vwUQIAAEDvcjD9OAAAAJWvqlQqlcpdRE9qbm7O1KlTO4Lm2tratLW1Zfv27UmSSZMmZcGCBRkxYsQr7tfU1JRZs2YlSZ555pmMGTNmj7WXLFmSM888Mxs2bEiS1NfXZ/v27dm5c2eSl09Zds8996SmpqaHHh09oaqqKsnLp5+DvbFH6Iw9QmfsETpjj9AZe4S+4NX24/ReXnvojD1CZ+wROmOP0Bl7hM7YI31HRR9JnSRjxozJsmXL8vnPfz4TJ05MVVVVBg0alMmTJ2f27Nn52c9+9qob4smTJ+eJJ57IZZddlnHjxmXnzp055JBD8ra3vS0333xz5s6dK6Dug0qlkhcv9sseoTP2CJ2xR+iMPUJn7BH6gp7sxykPrz10xh6hM/YInbFH6Iw9Qmfskb6j4o+kBgAAAAAAAKD3qPgjqQEAAAAAAADoPYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAAAAABAYYTUAAAAAAAAABRGSA0AAAAAAABAYYTUANBN/vEf/zFVVVV5+9vfXu5SgD7os5/9bKqqqjJ9+vRyl0If8N3vfjenn356hg8fnvr6+px++ul5+OGHy10WAMBe6ZeBg6Vn5kDomfsGITUAdJOf//znSZI3v/nNZa4E6IvaX0MmT55c5krozV566aW8733vywc/+MH8+Mc/zs6dO7Nt27b8+Mc/zjve8Y4sXLiw3CUCAOxBvwwcLD0zXaFn7luE1ADQTTTdwMFYunRpEg03+/fJT34yd911V0aNGpW77747LS0tWb9+fWbMmJEdO3bkox/9aLlLBADYg34ZOFh6ZrpCz9y3VJVKpVK5iwCAvm7Dhg057LDDkiRPPfVUjjvuuDJXBPQlq1evzute97okya9//ev8wR/8QZkrojdatWpVxo0bl7a2tjzyyCM55ZRTOsaee+65vP71r0+pVMoTTzyRCRMmlLFSAID/o18GDpaema7QM/c9jqQGgG7Q/qnwurq6HHvssWWuBuhr2l9DRowYodlmn+bPn5+2trY0Nja+otlOkqOOOipHH310kmT58uXlKA8AYK/0y8DB0jPTFXrmvkdITb82YcKEVFVVHdDXj370o3KXTRns2rUrDz74YD7zmc9kypQpGTNmTIYMGZL6+vq8+c1vzpe+9KVs27at3GVSRu1vlv/wD/8w1dXVueuuu/KOd7wjI0eOTG1tbU488cR85zvfKXOV9AbLli3LX//1X2fChAmpq6tLTU1NXve61+XUU0/NlVdemaeffrrcJVIGnZ3+8K677sqwYcNSVVWVCy64INu3by+yPHqJZ599Nkly6KGH7nV88ODBSZKNGzcWVRJQ4fTMdIV+mc7ol+kq/TL7omemK/TMfc/AchcA5dLS0pIXX3wxw4YNy+bNm1NdXZ3DDz/8FXPa2try/PPPJ0kaGhoyZMiQvPWtby1HuZTZokWL8va3v73jzzU1NampqcnmzZuzdOnSLF26NHfffXceeuihjv/s6F/a3yy/8Y1vzHve8558//vfT3V1dYYOHZpt27Zl8eLF+Yu/+Iv8+te/zlVXXVXmaimXL3/5y7nyyiuza9euJMmQIUNSU1OT1atXZ/Xq1fnJT36SE044wenv+qF9Ndw7d+7MFVdcka9+9aupqanJ17/+9VxyySXlKJFeoLa2NkmyePHirFu3LqNGjXrF+Nq1a5Mk9fX1hdcGVB49M12lX6Yz+mW6Qr/M/uiZ6Qo9c9/jSGr6rfr6+qxZsyaXXnppkuSkk07K2rVrX/E1b968jvnLli3L6tWrc8ghh5SrZMpo9erV+chHPpJ58+Zl7dq12bZtWzZt2pSNGzfmK1/5SgYMGJBFixb55G8/1v5m+fbbb89DDz2Uf//3f8/WrVuzZcuWLFu2LJMmTUqSXH311Vm5cmU5S6VM7rvvvvz93/99Bg4cmC996Uv57W9/m23btmXz5s154YUX8sADD2TmzJk56aSTyl0qZdD+GjJ58uSO7z333HM57bTT8tWvfjVHH310Hn74Yc12Pzdt2rQkybZt2/LOd74zDz/8cFpbW5MkzzzzTFpaWpIkr3nNa8pWI1A59Mx0lX6ZzuiX6Yx+mc7omekKPXMfVIJ+7txzzy0lKV1yySV7jN1yyy2lJKWRI0eWoTL6kr/8y78sJSldeOGF5S6FMti0aVOpqqqqlKTU0NBQevrpp/eYs2rVqtLgwYNLSUpf/OIXy1Al5faud72rlKR09dVXl7sUepnnn3++lKSUpPTLX/6yVCqVSvfdd19p5MiRpSSladOmlZ5//vkyV0m53HTTTaVRo0aVGhsbS6VSqXTZZZd17Jf2r2HDhpUOPfTQjj+vWbOmzFUDlUTPzMHSL/dv+mW6Qr/M/uiZ2R89c9/mSGr6vV/84hdJXr4uzv/rscceS5K86U1vKrAi+qJhw4YlSXbs2FHmSiiHxx57LKVSKcnLp6caP378HnPGjh3b8Ylfnwzvn9qvi9N+6jJo1/6J8GHDhmXs2LG56qqrcs4552Tjxo35zGc+k3nz5uWwww4rc5WUy6JFi7Ju3bq84Q1vSJJce+21ue+++zJz5syccMIJaWhoyLZt2/Liiy8mSUaNGpUjjzyyjBUDlUbPzMHSL/dv+mW6Qr/M/uiZ2R89c9/mmtT0ay+88EKee+65JHtvqpcuXZpk7804/cvWrVszZ86c/OAHP8iyZcuyfv36jlOF7O51r3tdGaqj3NrfLB9xxBH50Ic+tM957fujra2tkLroXRobG/PYY4/lC1/4QjZu3Ji/+qu/8gtdkvzfa8jRRx+ds846Kw888ECGDx+e2267Le985zvLXB3l9vjjjydJTjzxxI7vnXnmmTnzzDNfMW/69Om5//77M3Xq1ELrAyqbnpmu0C+zP/plukK/zP7omdkfPXPfJqSmX2v/RHhVVVUaGxtfMVYqlTrGvSnq3x566KF88IMfzOrVqzu+V1tbm8MPPzzV1S+fkGLdunVJkokTJ5alRsqr/c3yn/zJn2TgwH3/19r+iT2f1uuf/uVf/iUrV67MokWLcv311+f666/P6NGjM2PGjMycOTNve9vbyl0iZdL+GrJ8+fIsX748hx9+eH7yk5/kmGOOKXNl9Aa/+93vkiRjxozZ55y2trYsXrw4ycuNN0B30TPTGf0yndEv0xX6ZfZHz8z+6Jn7Nqf7pl9rb6jf8IY3pK6u7hVjv/rVr9LS0pJEw92fPfXUUznrrLOyevXqnHbaafnP//zPrF+/Plu3bs3//u//Zu3atXn44Yc75k+aNKmM1VIu7W+W209PtjelUqnjzdDkyZMLqYve5cgjj8xPfvKTzJ8/P5dccknGjRuXNWvW5Oabb86UKVPyvve9z1ED/VT7a8j555+f6urqPP/88/nxj39c5qroLV544YUkyciRI/c5Z+HChXnhhRdSU1OTd7/73UWVBvQDemb2R79MV+iX6Qr9MvujZ2Z/9Mx9m5Cafq0r19YaNGhQjj/++AKrojeZNWtWtm3blqlTp+bBBx/Mu971rj3+w7vjjjuSJEOHDs1xxx1XjjIpo23btuXpp59Osv83Qw8++GDWr1+fQYMG7XG6GfqP6urqvOMd78jXv/71rFy5MsuWLcuMGTOSJHfddVfuvffeMldI0TZt2pRVq1YlSa688spcc801SZKPfOQjeeCBB8pZGr1E+1FoO3fu3Oecr3/960mS9773vTn00EOLKAvoJ/TM7I9+mc7olzkQ+mX2Rs9MZ/TMfZuQmn5tf6cma7+21oQJEzJ48OBC66L3eOihh5IkH/zgBzv+w9tdS0tLvvnNbyZ5eR8NGDCg0Poov1/84hfZtWtXkv87jd3e/PM//3OS5IILLsjw4cMLqY3er7GxMffcc08OP/zwJOm45iP9x9KlS1MqlTJ06NCMHz8+n/zkJ/PRj340v//97/Onf/qneeKJJ8pdImX2mte8JkmyYsWKvY4vXbo0d9xxR6qrq/MP//APRZYG9AN6ZvZHv0xn9MscDP0yiZ6ZzumZ+zYhNf3W73//+zz55JNJ9t5wL1u2LElywgknFFoXvUtra2uSZP369XuMtbW15cILL0xzc3MSpy7rr9pPOZQkd999917n3HDDDbn//vtTX1/vzVA/tGPHjv2Or1u3Lps3b07iOn39UftrSGNjY8cvbq+77rqceeaZ2bRpU2bMmLHfX+hR+f7oj/4oSfZ65MiLL76YP//zP8+uXbvyiU98wpGMQLfSM9MZ/TKd0S/TGf0yndEz0xk9c98mpKbfWrFiRbZv355k76cue+aZZ5Iko0aNKrIsepk3v/nNSZJrr72241papVIpjz76aM4666zcf//9HW+Q2ufSv7S/WW5oaMi9996bz3/+89m6dWuSZO3atbn88svzN3/zN0mSW2+9NUcddVTZaqU8vvWtb2XSpEm57rrr8stf/jKlUinJy8343XffnWnTpmXHjh05+eSTM2XKlDJXS9HaX0N2fy8yYMCA3HHHHZk4cWKeffbZnHfeeR2/BKb/mTlzZpKXT4N57bXXdlyLb8mSJfnjP/7jPPnkkzn11FM7jkAC6C56ZjqjX6Yz+mU6o1+mM3pmOqNn7uNK0E/NmTOnlKQ0YsSIvY6/8Y1vLCUpHXXUUaWrr766tGLFioIrpDdYuHBhadCgQaUkpSSl2traUm1tbSlJafTo0aV77723Y2zJkiXlLpcymDRpUilJ6YYbbiiddNJJpSSl6urqUkNDQ8feGDRoUOnmm28ud6mUycUXX9yxF5KUampqSg0NDaWqqqqO7731rW8trVu3rtylUgbHH398KUnpxhtv3GOsubm5NGrUqFKS0rvf/e7Srl27ylAhvcEFF1zQ8XoxYsSI0ujRozv+fM4555RaWlrKXSJQgfTMdEa/TGf0y3RGv0xn9Mx0hZ6573IkNf3W/q6tlSR/+7d/m5EjR+Z3v/tdbrnllgwbNqzI8uglTj311Pz3f/93zjjjjNTW1qa6ujrjxo3LrFmz8uSTT2bIkCFJkkGDBjntUD+0c+fOjmvfvOlNb8qPfvSjXH755TnqqKOyZcuWHH300bnooouyfPnyXHTRRWWulnL5xCc+kdmzZ+e8887LuHHjMmTIkLS0tGTUqFE555xzctttt+WRRx7puIYO/Udra2tWrlyZZO9HqL3+9a/PPffck6FDh+b73/9+rrjiioIrpLe45ZZb8tWvfjWNjY3Zvn17duzYkenTp+eOO+7Ivffem7q6unKXCFQgPTOd0S+zP/plukK/zP7omekqPXPfVVUq/f/n0AAAAAAAAACAHuZIagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAKI6QGAAAAAAAAoDBCagAAAAAAAAAK8/8Bw8e7bTLWSY8AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 1980x756 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from IPython.display import display\n",
    "fit_result = fit_full_trim(seed=123)\n",
    "elasticity_table = el.elasticity_full_trim(fit_result)\n",
    "elasticity_table.to_csv(PROJECT_ROOT / 'outputs' / 'bayesian_elasticity.csv', index=False)\n",
    "fig, axes = el.plot_combined_elasticity(elasticity_table)\n",
    "fig.savefig(PROJECT_ROOT / 'outputs' / 'bayesian_elasticity.png', dpi=600, bbox_inches='tight')\n",
    "fig.savefig(PROJECT_ROOT / 'outputs' / 'bayesian_elasticity.pdf', bbox_inches='tight')\n",
    "display(elasticity_table)\n",
    "import matplotlib.pyplot as plt\n",
    "display(fig)\n",
    "plt.close(fig)\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}