From 22f7053fb60d40dcfd8450409aab7cbe1c08d83a Mon Sep 17 00:00:00 2001 From: lgrcia <20612771+lgrcia@users.noreply.github.com> Date: Tue, 28 May 2024 14:28:43 +0200 Subject: [PATCH] docs: more context --- docs/notebooks/tutorials/tess_search.ipynb | 46 +++++++++++++--------- 1 file changed, 28 insertions(+), 18 deletions(-) diff --git a/docs/notebooks/tutorials/tess_search.ipynb b/docs/notebooks/tutorials/tess_search.ipynb index c1a84b4..529de4e 100644 --- a/docs/notebooks/tutorials/tess_search.ipynb +++ b/docs/notebooks/tutorials/tess_search.ipynb @@ -38,20 +38,9 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": 12, "metadata": {}, - "outputs": [ - { - "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxYAAAEiCAYAAABkykQ1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB6f0lEQVR4nO3de3gU5d0//vfuZpNsgIQQA2EFFsoq4RCRiIfY4lcaNAXUL4+21dZ61u9Da1XEA0XbeniqtFq1th5bD6jVUvso1gKaC1MRbakHYmwIxhqFDRICEeKGQ0g2u/v7I797uGd2ZndmdzbZJO/XdXFpNrOTe0738XPf44hGo1EQERERERGlwNnfCSAiIiIiooGPDQsiIiIiIkoZGxZERERERJQyNiyIiIiIiChlbFgQEREREVHK2LAgIiIiIqKUsWFBREREREQpY8OCiIiIiIhSltXfCchEkUgELS0tGDFiBBwOR38nh4iIiIioX0SjUezfvx9erxdOZ/wxCTYsdLS0tGD8+PH9nQwiIiIiooywY8cOjBs3Lu42bFjoGDFiBIDeE5ifn9/PqSEiIiIi6h8dHR0YP368Uj+Ohw0LHSL8KT8/nw0LIiIiIhryzEwP4ORtIiIiIiJKGRsWRERERESUMjYsiIiIiIgoZZxjQURERDRIhcNhhEKh/k4GZTC32w2Xy2XLvtiwICIiIhpkotEoWltb8dVXX/V3UmgAGDlyJEpKSlJ+fxsbFkRERESDjGhUjB49Gnl5eXzhL+mKRqM4dOgQ9uzZAwAYO3ZsSvtjw4KIiIgGvUAggKamJvj9fvh8vv5OTlqFw2GlUVFUVNTfyaEM5/F4AAB79uzB6NGjUwqL4uRtIiIiGvSamprQ2dmJpqam/k5K2ok5FXl5ef2cEhooxL2S6nwcNiyIiIho0PP7/fB4PPD7/f2dlD7D8Ccyy657haFQRERENOj5fL5BHwJF1N84YkFERERElCGi0Sj+3//7fxg1ahQcDgfq6upw+umnY8mSJf2dtITYsCAiIiKijHDppZfC4XDA4XAgOzsbfr8fd955J3p6egAAGzZsUH7vdDpRUFCAWbNm4eabb8auXbti9tfR0YFbb70VpaWlyM3NRUlJCebNm4eXX34Z0WjUMB0fffQRzjnnHIwePRq5ubmYOHEizj//fGX1JNmKFSvgcrlw7733xvxu5cqVqvSOGzcOl112me5+hNdffx0rV67EmjVrsGvXLsyYMcPMqcsIbFgQERERUcb41re+hV27duHTTz/FDTfcgNtvvz2m0v7JJ5+gpaUF77//PpYtW4Y33ngDM2bMQH19vbLNV199hVNPPRXPPvssli9fjtraWmzcuBHnn38+br75ZgSDQd2/39bWhsrKSowaNQrV1dX4+OOP8fTTT8Pr9eLgwYMx2z/11FO4+eab8dRTT+nuLz8/H7t27cIXX3yBP/zhD3jttddw0UUXGR7/Z599hrFjx+LUU09FSUkJsrIGzswFNiyIiIaYQCCAmpoaBAKB/k4KEVGMnJwclJSUwOfz4Yc//CHmzZuHV199VbXN6NGjUVJSgmOPPRYXXHAB/vGPf6C4uBg//OEPlW1uueUWbN++He+++y4uueQSTJs2Dcceeyyuuuoq1NXVYfjw4bp//x//+AeCwSCeeOIJzJo1C5MmTcLcuXPxwAMPYNKkSapt33rrLXR2duLOO+9ER0cH/vnPf8bsz+FwoKSkBF6vF/Pnz8e1116LN954A52dnTHbXnrppbjmmmvQ3NwMh8OBiRMn6qbR4XDglVdeUX02cuRIrFy5EgDw7LPPYvjw4fj000+V3//oRz9CaWkpDh06pLtPO7BhQUQ0xAylZTeJaODzeDzo7u5OuM3ixYvxj3/8A3v27EEkEsGqVatw4YUXwuv1xmw/fPhww5GAkpIS9PT0YPXq1XHDpQDgySefxPe+9z243W5873vfw5NPPmnqeCKRiBLeJXvwwQdx5513Yty4cdi1axfef//9hPvTc/HFF2PBggW48MIL0dPTg7Vr1+KJJ57A888/n9ZliNmwICIaYobisptElLz+GuWMRqN44403UF1djW9+85sJty8tLQUAbN++HV9++SXa29uVz6w45ZRTcMstt+D73/8+jjrqKMyfPx/33nsvdu/erdquo6MD//u//4sf/OAHAIAf/OAHePHFF3HgwAHDfX/66ad47LHHMHv2bIwYMSLm9wUFBRgxYgRcLhdKSkpQXFxsOf3C448/jl27duHaa6/FFVdcgdtvvx0nnHBC0vszgw0LIqIhxufzobKykktvEpEpfT3KuWbNGgwfPhy5ubmYP38+zj//fNx+++0JvydGFxwOR8KRBuHuu+/G8OHDlX/Nzc0AgLvuugutra147LHHMH36dDz22GMoLS1VzeH405/+hMmTJ2PmzJkAgOOPPx4+nw9//vOfVX8jGAxi+PDhyMvLw5QpUzBmzBg8//zzptKXisLCQjz55JN49NFHMXnyZPzkJz9J+99kw4KIiIiIDPX1KOfcuXNRV1eHTz/9FJ2dnXjmmWcwbNiwhN/7+OOPAQATJ05EcXExRo4cicbGxrjfWbx4Merq6pR/cthUUVERvvOd7+DXv/41Pv74Y3i9Xvz6179Wfv/kk0+ioaEBWVlZyr+tW7fGTOIeMWIE6urqsGXLFhw8eBAbN27Esccea+WUxNBrPOm9NXvjxo1wuVzYtWuX7sRzuw2caeZERERE1Of6+uWCw4YNs9yI6ezsxO9//3ucdtppSvjQBRdcgOeeew633XZbzDyLAwcOIDc3F6NGjcKoUaMS7j87OxuTJ09WKuf19fX44IMPsGHDBtX39+3bh9NPPx2NjY1KGJbT6bS9UVZcXKxaXvfTTz+NmZT9z3/+E7/61a/wt7/9DcuWLcOPf/xjPPPMM7amQ4sNCyIiIiIaUPbs2YPDhw9j//792Lx5M+655x58+eWXePnll5Vt7rrrLmzYsAEnn3wy7rrrLsyePRtutxtvv/02VqxYgffffx8jR46M2feaNWuwatUqXHDBBTj22GMRjUbxt7/9DevWrcPTTz8NoHe04qSTTsJpp50W8/0TTzwRTz75pO57LezyzW9+Ew899BAqKioQDoexbNkyuN1u5ff79+/HRRddhGuvvRbz58/HuHHjcOKJJ+Lss8/Gt7/97bSliw0LIiIiIhpQpkyZAofDgeHDh+NrX/sazjzzTCxduhQlJSXKNqNGjcK//vUv/PKXv8QvfvELBAIBFBYWoqysDPfeey8KCgp09z1t2jTk5eXhhhtuwI4dO5CTk4NjjjkGTzzxBC666CJ0d3fjj3/8I5YtW6b7/fPOOw/33Xcf7r777rQcOwDcd999uOyyyzBnzhx4vV48+OCD2Lx5s/L76667DsOGDVPSUFZWhrvvvhv//d//jYqKChx99NFpSZcjanZ2yxDS0dGBgoICBINB5Ofn93dyiIiIiEw7fPgwtm3bhkmTJiE3N7e/k0MDQLx7xkq9mJO3iYiIiIgoZWxYEBERERFRytiwICIiIiKilLFhQUREREREKWPDgoiIiGgQ4vo8ZJZd9wobFkRERESDiHifgfaFaURGxL0ivwsjGf36HouNGzfi3nvvxebNm7Fr1y6sXr0aixYtivudDRs2YOnSpWhoaMD48ePx05/+FJdeeqnutr/85S+xfPlyXHfddfjNb35je/qJiIiIMo3L5cLIkSOxZ88eAEBeXh4cDkc/p4oyUTQaxaFDh7Bnzx6MHDkSLpcrpf31a8Pi4MGDmDlzJi6//HKce+65Cbfftm0bFi5ciMWLF+P5559HTU0NrrzySowdOxZVVVWqbd9//308/vjjOO6449KVfCIiIqKMJF4UJxoXRPGMHDlS9XLBZPVrw2L+/PmYP3++6e0fe+wxTJo0Cffddx8AYOrUqXjnnXfwwAMPqBoWBw4cwIUXXog//OEP+MUvfmF7uomIiIgymcPhwNixYzF69GiEQqH+Tg5lMLfbnfJIhdCvDQurNm3ahHnz5qk+q6qqwpIlS1SfXX311Vi4cCHmzZtnqmHR1dWFrq4u5eeOjg5b0ktERETUn1wul22VRqJEBlTDorW1FWPGjFF9NmbMGHR0dKCzsxMejwerVq1CbW0t3n//fdP7XbFiBe644w67k0tERERENGQMqlWhduzYgeuuuw7PP/88cnNzTX9v+fLlCAaDyr8dO3akMZVERERERIPPgBqxKCkpwe7du1Wf7d69G/n5+fB4PNi8eTP27NmD8vJy5ffhcBgbN27EQw89hK6uLt3hwJycHOTk5KQ9/UREREREg9WAalhUVFRg3bp1qs/Wr1+PiooKAEBlZSXq6+tVv7/ssstQWlqKZcuWMcaQiIiIiChN+rVhceDAATQ1NSk/b9u2DXV1dRg1ahQmTJiA5cuXY+fOnXj22WcBAIsXL8ZDDz2Em2++GZdffjn+/ve/48UXX8TatWsBACNGjMCMGTNUf2PYsGEoKiqK+ZyIiIiIiOzTr3MsPvjgA8yaNQuzZs0CACxduhSzZs3Cz3/+cwDArl270NzcrGw/adIkrF27FuvXr8fMmTNx33334Yknnoh5hwUREREREfUtRzQajfZ3IjJNR0cHCgoKEAwGkZ+f39/JISIiIiLqF1bqxYNqVSgiIiIiIuofbFgQEREREVHK2LAgIiIiIqKUsWFBREREREQpY8OCiIiIiIhSxoYFERERERGljA0LIiIiIiJKGRsWRERERESUMjYsiIiIiIgoZWxYEBERERFRytiwICIiIiKilLFhQUREREREKWPDgoiIiIiIUsaGBRERERERpYwNCyIiIiIiShkbFkRERERElDI2LIiIiIiIKGVsWBARERERUcpsbVhEo1E7d0dERGkQCARQU1ODQCDQ30khIqJBxHLD4tJLL8XBgwdjPt++fTtOO+00WxJFRETp09TUhM7OTjQ1NfV3UoiIaBCx3LD46KOPcNxxx2HTpk3KZ8888wxmzpyJo446ytbEERGR/fx+PzweD/x+f38nhYiIBpEsq1947733cMstt+D000/HDTfcgKamJrz22mu4//77cdVVV6UjjUREZCOfzwefz9ffySAiokHGcsPC7Xbj3nvvRV5eHv7nf/4HWVlZeOutt1BRUZGO9BERERER0QBgORQqFArhhhtuwK9+9SssX74cFRUVOPfcc7Fu3bp0pI+IiIiIiAYAyyMWs2fPxqFDh7BhwwaccsopiEajuOeee3Duuefi8ssvxyOPPJKOdBIRERERUQazPGIxe/Zs1NXV4ZRTTgEAOBwOLFu2DJs2bcLGjRst7Wvjxo04++yz4fV64XA48MorryT8zoYNG1BeXo6cnBz4/X6sXLlS9ftHH30Uxx13HPLz85Gfn4+Kigq89tprltJFRERERETWWG5YPPnkkxg2bFjM57NmzcLmzZst7evgwYOYOXMmHn74YVPbb9u2DQsXLsTcuXNRV1eHJUuW4Morr0R1dbWyzbhx4/DLX/4SmzdvxgcffIBvfvOb+L//9/+ioaHBUtqIiIiIiMg8R9TiW+2am5vj/n7ChAnJJcThwOrVq7Fo0SLDbZYtW4a1a9diy5YtymcXXHABvvrqK7z++uuG3xs1ahTuvfdeXHHFFabS0tHRgYKCAgSDQeTn55s+BiIiIiKiwcRKvdjyHIuJEyfC4XAY/j4cDlvdpWmbNm3CvHnzVJ9VVVVhyZIlhmn5y1/+goMHD8ZdtaqrqwtdXV3Kzx0dHbakl4iIiIhoqLDcsPjwww9VP4dCIXz44Ye4//77cdddd9mWMD2tra0YM2aM6rMxY8ago6MDnZ2d8Hg8AID6+npUVFTg8OHDGD58OFavXo1p06YZ7nfFihW444470pp2IiIiIqLBzHLDYubMmTGfzZ49G16vF/feey/OPfdcWxKWiilTpqCurg7BYBD/+7//i0suuQRvvfWWYeNi+fLlWLp0qfJzR0cHxo8f31fJJSIiIiIa8Cw3LIxMmTIF77//vl2701VSUoLdu3erPtu9ezfy8/OV0QoAyM7Oht/vBwCccMIJeP/99/Hggw/i8ccf191vTk4OcnJy0pdwIiIiIqJBznLDQjv/IBqNYteuXbj99ttxzDHH2JYwPRUVFTEv4lu/fn3Ct35HIhHVHAoiIiIiIrKX5YbFyJEjYyZvR6NRjB8/HqtWrbK0rwMHDqCpqUn5edu2bairq8OoUaMwYcIELF++HDt37sSzzz4LAFi8eDEeeugh3Hzzzbj88svx97//HS+++CLWrl2r7GP58uWYP38+JkyYgP379+OFF17Ahg0bVEvSEhERERGRvSw3LN58803Vz06nE8XFxfD7/cjKsra7Dz74AHPnzlV+FvMcLrnkEqxcuRK7du1SLW87adIkrF27Ftdffz0efPBBjBs3Dk888QSqqqqUbfbs2YOLL74Yu3btQkFBAY477jhUV1fjjDPOsHqoRERERERkkuX3WAwFfI8FEREREVEa3mPx6quvmv7j55xzjultiYiIiIhocDDVsIj3NmyZw+FI6wvyiIiIiIgoM5lqWEQikXSng4iIiIiIBjCnmY1GjRqFvXv3AgAuv/xy7N+/P62JIiIiIiKigcVUw6K7uxvBYBAA8Mwzz+Dw4cNpTRQREREREQ0spkKhKioqsGjRIpxwwgmIRqO49tprVW+6lj311FO2JpCIiIiIiDKfqYbFH//4RzzwwAP47LPP4HA4EAwGOWpBREREREQKy++xmDRpEj744AMUFRWlK039ju+xICIiIiJKw3ssZNu2bUs6YURERERENDiZmrxNREREREQUDxsWRERERESUMjYsiIiIiIgoZWxYEBENMYFAADU1NQgEAv2dFCIiGkRMTd7u6OgwvUOuokRElNmamprQ2dmJpqYm+Hy+/k4OERENEqYaFiNHjoTD4TC1w3A4nFKCiIgovfx+P5qamuD3+/s7KURENIiYali8+eabyv9v374dP/nJT3DppZeioqICALBp0yY888wzWLFiRXpSSUREtvH5fBypICIi21l+QV5lZSWuvPJKfO9731N9/sILL+D3v/89NmzYYGf6+gVfkEdEREREZK1ebHny9qZNmzB79uyYz2fPno333nvP6u6IiIiIiGgQsNywGD9+PP7whz/EfP7EE09g/PjxtiSKiIiIiIgGFlNzLGQPPPAAzjvvPLz22ms4+eSTAQDvvfcePv30U7z00ku2J5CIiOwVCASUyduca0FERHaxPGKxYMEC/Oc//8HZZ5+Nffv2Yd++fTj77LPxn//8BwsWLEhHGomIyEbycrNERER2sTx5eyjg5G0iGsw4YkFERGaldfI2ALz99tv4wQ9+gFNPPRU7d+4EADz33HN45513ktkdERH1IZ/Ph8rKSjYqiIjIVpYbFi+99BKqqqrg8XhQW1uLrq4uAEAwGMTdd99tewKJiIiIiCjzWW5Y/OIXv8Bjjz2GP/zhD3C73crnX//611FbW2tr4oiIiIiIaGCw3LD45JNPcNppp8V8XlBQgK+++srSvjZu3Iizzz4bXq8XDocDr7zySsLvbNiwAeXl5cjJyYHf78fKlStVv1+xYgVOPPFEjBgxAqNHj8aiRYvwySefWEoXERERERFZY7lhUVJSoruSyDvvvIOvfe1rlvZ18OBBzJw5Ew8//LCp7bdt24aFCxdi7ty5qKurw5IlS3DllVeiurpa2eatt97C1VdfjX/9619Yv349QqEQzjzzTBw8eNBS2oiIiIiIyDzL77G46qqrcN111+Gpp56Cw+FAS0sLNm3ahBtvvBE/+9nPLO1r/vz5mD9/vuntH3vsMUyaNAn33XcfAGDq1Kl455138MADD6CqqgoA8Prrr6u+s3LlSowePRqbN2/WHWkhIiIiIqLUWW5Y/OQnP0EkEkFlZSUOHTqE0047DTk5ObjxxhtxzTXXpCONik2bNmHevHmqz6qqqrBkyRLD7wSDQQDAqFGjDLfp6upSJqEDvctqERER0eDD5ZaJ0sdyw8LhcODWW2/FTTfdhKamJhw4cADTpk3D8OHD05E+ldbWVowZM0b12ZgxY9DR0YHOzk54PB7V7yKRCJYsWYKvf/3rmDFjhuF+V6xYgTvuuCMtaSYiIqL+JxoUPT09CIVCaGpqYsOCyGaW51g8++yz+Pjjj5GdnY1p06bhpJNOwvDhw3H48GE8++yz6Uhj0q6++mps2bIFq1atirvd8uXLEQwGlX87duzooxQSERFRXxBvnAcAj8cDv9/fzykiGnwsNywuvfRSnHTSSXjppZdUnweDQVx22WW2JUxPSUkJdu/erfps9+7dyM/Pjxmt+PGPf4w1a9bgzTffxLhx4+LuNycnB/n5+ap/RERENHj4/X54PB6UlpbyBZFEaWI5FAoA7rjjDlx00UWor6/H7bffbnOSjFVUVGDdunWqz9avX4+Kigrl52g0imuuuQarV6/Ghg0bMGnSpD5LHxEREWUmn8/HxgRRmlkesQCAH/zgB/j73/+Oxx9/HN/+9reVoUWrDhw4gLq6OtTV1QHoXU62rq4Ozc3NAHpDlC6++GJl+8WLF+Pzzz/HzTffjMbGRjzyyCN48cUXcf311yvbXH311fjjH/+IF154ASNGjEBraytaW1uTTiMRERERESVmuWHhcDgAAKeccgreffddNDU14dRTT8X27dst//EPPvgAs2bNwqxZswAAS5cuxaxZs/Dzn/8cALBr1y6lkQEAkyZNwtq1a7F+/XrMnDkT9913H5544gllqVkAePTRRxEMBnH66adj7Nixyr8///nPltNHRERERETmOKLRaNTKF5xOJ1pbWzF69GgAwKFDh3DhhReipqYGBw8eRDgcTktC+1JHRwcKCgoQDAY534KIiIiIhiwr9WLLIxa33XabamnZvLw8rF69Gtdffz1fQEdERERENERZHrEYCjhiQURERERkrV5salWoV199FfPnz4fb7carr75quJ3D4cDZZ59tLbVERERERDTgmRqxkOdVOJ3G0VMOh4NzLIiIiIiIBgnbRywikYju/xMREREREQFJvseCiIiIaCAJBAKoqalBIBDo76QQDVqmRix++9vfmt7htddem3RiiIiIiNKhqakJnZ2daGpq4hu4idLE1ByLSZMmmduZw4HPP/885UT1N86xICIiGlwCgQCamprg9/vZsCCywPY5Ftu2bbMlYURERNS/hmoF2+fzDanjJeoPnGNBREQ0hMghQUREdjI1YqH1xRdf4NVXX0VzczO6u7tVv7v//vttSRgRERHZz+/3KyMWRER2stywqKmpwTnnnIOvfe1raGxsxIwZM7B9+3ZEo1GUl5enI41ERERkE4YEEVG6WA6FWr58OW688UbU19cjNzcXL730Enbs2IH/83/+D77zne+kI41ERGQjLrtJRETpYLlh8fHHH+Piiy8GAGRlZaGzsxPDhw/HnXfeiV/96le2J5CIiOzFGHsiIkoHyw2LYcOGKfMqxo4di88++0z53ZdffmlfyoiIKC38fj88Hg9j7ImIyFaW51iccsopeOeddzB16lQsWLAAN9xwA+rr6/Hyyy/jlFNOSUcaiYjIRoyxJyKidLDcsLj//vtx4MABAMAdd9yBAwcO4M9//jOOOeYYrghFRERERDREmXrz9lDDN28TEREREaXhzdtGDhw4gEgkovqMFXEiIiIioqHH8uTtbdu2YeHChRg2bBgKCgpQWFiIwsJCjBw5EoWFhelIIxER0YDGJX6JaCiwPGLxgx/8ANFoFE899RTGjBkDh8ORjnQRERENGvISv5w4T0SDleWGxUcffYTNmzdjypQp6UgPGQgEAmhqaoLf72ehlEF4XYjIDL/fr+QVRESDleVQqBNPPBE7duxIR1ooDr7QKjPxuhCRGT6fD5WVleyAIKJBzfKIxRNPPIHFixdj586dmDFjBtxut+r3xx13nG2JoyPS3ds11Hre7Tpe9kISEWU+Oc8HMKTKO6K+ZLlh0dbWhs8++wyXXXaZ8pnD4UA0GoXD4UA4HLY1gdQr3S+0Gmrxv3YdL180RkSU+bSjy0OpvCPqS5ZDoS6//HLMmjULmzZtwueff45t27ap/mvFxo0bcfbZZ8Pr9cLhcOCVV15J+J0NGzagvLwcOTk58Pv9WLlyZcr7pN6ed4/HM2R63ofa8VrFFWyIaDDx+/1wu93o6elBYWEh838bDYbyYjAcQ6awPGIRCATw6quv2vJAHjx4EDNnzsTll1+Oc889N+H2YqnbxYsX4/nnn0dNTQ2uvPJKjB07FlVVVUntk3oNtZ73oXa8Vg21ESwiGvxCoRAAoL29HZWVlf2cmsFjMJQXg+EYMoXlhsU3v/lNfPTRR7Y0LObPn4/58+eb3v6xxx7DpEmTcN999wEApk6dinfeeQcPPPCA0rCwuk/qNdTmWFB8nDtCRIOJvMAG8zV7DYbyYjAcQ6aw3LA4++yzcf3116O+vh5lZWUxk7fPOecc2xKntWnTJsybN0/1WVVVFZYsWZLSfru6utDV1aX83NHRkdL+BiK21knGEZ3BjR0JNNTIFUfe8/YaDOXFYDiGTGG5YbF48WIAwJ133hnzu3RP3m5tbcWYMWNUn40ZMwYdHR3o7OyEx+NJar8rVqzAHXfcYUcSByy21omGjsbGRoRCITQ2NrIwpSGBFUeivmF58nYkEjH8N1BXhFq+fDmCwaDybyi+p4NrrBMRJTbQJ3kO9PQTUWaz1LAIhULIysrCli1b0pWeuEpKSrB7927VZ7t370Z+fn7SoxUAkJOTg/z8fNU/IiODpWAeLMdB1pWWlsLj8aC0tLS/kzLgiLDR+vp61NbW9ndyLONLPSnT6ZVNg628GmzHI7PUsHC73ZgwYUK/jUxUVFSgpqZG9dn69etRUVHRL+npb4P5xsxkg6VgbmxsRGdnJxobG/s7KdTHOEKZPDlctKWlpR9TkpyhutQ2y8uBQ6+MTabczbRrLqdnsNQj9FgOhbr11ltxyy23YN++fSn/8QMHDqCurg51dXUAepeTraurQ3NzM4DeEKWLL75Y2X7x4sX4/PPPcfPNN6OxsRGPPPIIXnzxRVx//fWm9zmYDOYbM5NoM6eBXjCL4xmooYtE/cnn88Hr9QKA8t9EMqmCM1QblSwvBw69MjaZcjfTrrmcnoFej4jHEY1Go1a+MGvWLDQ1NSEUCsHn82HYsGGq31sZGt6wYQPmzp0b8/kll1yClStX4tJLL8X27duxYcMG1Xeuv/56bN26FePGjcPPfvYzXHrppab3aUZHRwcKCgoQDAYzOiwqk1Z2yaS02K2mpkZZHGAwrH0ujsftdiMrK2tQXjOKbzA/r5nI7jxkMF+/dB1bOs/ZYL4eA1ltbS1aWlrg9XpRXl5u236Tvd4D+T6xUi+2vCrUokWLkk1XjNNPPx3x2jV6DYHTTz8dH374YdL7HEzsXOXC6g2v3V5vudqB/BDJBtuKWVx2cXCy8rzJ8wQA8D4wYFceZnceksry4JmeL6dr6fN0rgrF5dozU3t7u+q/gD33P693fJYbFrfddls60kH9zOqDot1er+AcLA/fYFumcLAdD/Wy8rz5/X6lUTHQn890sisPs/uZS6Whkun5sl2NML0KZLp6sM2kOVMbdHanqy+O0+zfsKteov17ydyjgUBAyXO3bt2aUfeA3Sw3LITNmzfj448/BgBMnz4ds2bNsi1RZF2qD3NhYSE6OztRWFhoanvtg6VXcA60nv5MzfgzGc9Z5rDyvIlrNZCez/6QqXlYKg0Vq3l9InbkAdp92JGX6FUgxWT7lpYW0w2LQCCgLHBRWlqqmzYzac7UBp3d6eqL4zT7N+yql2j/XjL3qDzPY7DPb7TcsNizZw8uuOACbNiwASNHjgQAfPXVV5g7dy5WrVqF4uJiu9NIcYgMuaenB6FQKOmHua2tTfXfROSKifxzOgqIvpKpGX8m4znLDOK5KywsjHkm9cg9t7xu+gZiT64ZeuEhqbAjD0hHPqJXgfR6vcp9H498rcScUpHOZNNnZyPVznvJ7nT19PTA7XantTFuJc16dRKr58yOc+T3+9HQ0IBIJGJ60YeByvKqUNdccw3279+PhoYG7Nu3D/v27cOWLVvQ0dGBa6+9Nh1ppDhEhgzA9hUGEq1kordcqdVVGNK9WoqV/ZtZpSGTVnfRYzV9qR6P0TnL9PM0kOmdW/HctbS0mHr+5J5bo/0O9Wto94oymbJCjd/vh9vtRk9Pj+61tXrdE+3PDDF6YtcoChC7+lUgEEB7ezvKysoSjlaIa9XY2Iienh44nc64lWUz50ykB0DKz5Wd95Kdq4SJRlhWVm+fdTLHaeVcWg3bTpYd58jn82HBggU466yzbA3Dy0SWRyxef/11vPHGG5g6dary2bRp0/Dwww/jzDPPtDVxlJhdE3GLi4vR0tKiGnGy2otktbdCjjnU/o1UVl2Qh63lAiLR/uL1ZIgeXpfLhXA4nFJ8ZjqZuWba3rhUegqNzhlHMtJHW1CKkQqgt2LW3t6e8PlzOp2IRCJwOo/0LWmvWWNjI0KhEBobG4fENdSLoxYVy0Ag0O89w3blI0aLbYi/YXXuTbz9mSWPmBsdp5Xj19t269atCIfDpuLbxbUSkQBiBT0jZo8/mfOr/X5jYyPC4bAtowJmzqmV6yHf41bzD230RSoLS8j1AFGn6e9wRu35Std8n0xgecQiEonA7XbHfO52uxGJRGxJFMUyehOlHP5g1MKXvyv+v7a2VrU/vVCoRD344oEV/5V7K8z2JAjanqpkexlEGkRYmDgGACn1Woie3XA4rDonZnpX9I4lEAiguroa1dXVKfcIy2kwM+oip8eutbTtetfHUO8lN0M+t2LUsK2tDZWVlSgvLzfVszZ9+nS43W64XC4lX+iLEIZMpn1Om5ubVXkJ0L/3p1E+opefJ2I0QiDv28p9IO7JwsLCpPK1np4e5b9GL+60UibobSvi2s3Et4seavGGesC4/DDz7Ijr1NDQoHyW7MT7UCiESCRiWM5auUfNvCTV6LxrP0+14Sv2J18fo/NtpswVz257e7thntiXz7N8vgKBgGrUeLCVeZYbFt/85jdx3XXXqYbQd+7cieuvv35QrPGfqfQebr3wB70HRc48tN+pr6/HmjVrlIw9FAop3403/Cc/GCJW12plUt5OO7cj2Yqp3+9XemELCwtjCohkK00FBQXKf8U5Eb1PiQo7vVABbQPIKvk6a3vLElUsRXq6urrQ2Nhoy0iK9v6U02El886UcJHBTlxv0aOo1ykgnpnS0tL+TGqf0eY5wWBQ+V1hYSFqamqwdevWpN9WH69hYDZUUzy3ouJuNQROMJpnIc5BWVmZpTxBPO/t7e1J5WtiiXjtUvFWO020x6GdXwFAaUyboS0/xH2gDUPUPjva6yrKYNH56na7k8pzxT1g1IgxWyYJIj2RSMTwXjQ673JjsqamRlXPqKmpQXFxsaX8Q+zP5XIpn+mF15kpI8R5cjqdcUP0+rK8kc+j/PdcLpdSFxssjQvLDYuHHnoIHR0dmDhxIiZPnozJkydj0qRJ6OjowO9+97t0pJEQ+xDX1tYqvSRer1d1w8Z7UMR+tJOH5AxdZAyJegQEkSYAMTGtZgtN7WiX1XhU8bcAICcnB0BvoSl/bjZGUowmrFu3DuvWrUN1dTUOHToEAOju7la2k89BdnY21qxZo/uCSJ/Ph6ysLFVhq20AWZXKqINITyQSSao3Vi8Wv6urSzkW7e/NZt7sNTdHLsDlyr/V3je5B1ev8Wtn7PVAJHcmtLe3x/SmWqX3nFqp2Og9t3J+bjYPiPecJdMhoK38m6nQaclvMpfvaSuhrEbHIZSXl8Pj8SAcDpsq4/T2J+4D7ciO2evqdDpNVbaNzr3P50NpaalhWJbc4NVrBGn/hrifnU6nYZqN8gHtORH7ikQi6OzsjDtSoEfsLy8vDwDgcDh0G6hmyjufz4eqqirk5OQoYVVyxIbVkfVkRwb1jk+EWYoG/LRp05RtBkuHmuU5FuPHj0dtbS3eeOMN5SaeOnUq5s2bZ3vi6AgRyy7e4CombANAUVGRKkZPG8crMmjtSk1FRUXKKgVaRvGicviVGOVoa2tT9XyKOO9EK1Vph6nXrVsHoLcFL5b1Mxu7Krarr6+H1+tFT08Purq6LMWzamM8BRGLrg2BEpP6XC6X0rtptIyhNr5aPjYzq7PI5729vV1pjFhdfUvsJzs7G52dnXC5XPD7/arY30SxrdrCvqenR7mH2tvbVYWvyETNxJaLnj+PxzNkK7OJBAIB1b0pX/vq6mpLsckOh0P571CcU6GlPf45c+Yov5Ofv7a2NoTDYVRXVxsuP6olYr57enpQX1+PvXv3ory83PI8DrG9+P9kVrgRzxkA7N27V9URJJcTZu4Hbb5RVlaGqqoqpZwyO4+gvLxclW/K3xF5jNX5fvL1AnrLQbHUrqj8Wm2wmF2qV5vnactgvfTqzQ3UO17t7+TvyiMi2nxYS6/XXByjFfJ8CuBIh0V3d3dSc5NEWSo6O7XpMbrn9UKxxPUCjhyv1WVj5Xtc3leq85zk505MLxCNwb6Yj5lOSb3HwuFw4IwzzsAZZ5xhd3pIh5xx5OXlqRoVgPom13tQ4j08LpdL1bBwOp2qTFBLZGqCyEzE8KXc6HG73apRFu3DIj/0gHpY1mql1O8/8sIvOUxP/n2iGFBxbGKoORwOK2lyuVyqUD+5cBYND71l5AKBgGqJOVEQiEJPpC0RkTZxvjo7O3XDFcweo9iPGDUR506IV5mQJzZ2dnYq+xCNFPF9EbJhdglqsd/Bkrmmg7YxblS5NVPwGcWcy+GQqcRMD1ThcBg1NTXIzs5GMBhUJleKc/Daa68hEomo8qlE5PwCONIBYXXi8969exEKheByubB3717VBFUxad9M40TOK0WF3ko6RD4jRioFsd9UJqpr8zA5zzTTsyx3Dsnli/gc6K20ejweVYNFbKN3DkUZLK6hGA0X5yASiajyTDndNTU1KCwsVDq66uvrUVBQgO7ubuV45PxXpMHoeOXGqJgALLjdbmRnZyvfM9qH6BgTDh8+rPy/1WWIRSNUS4wMmbmXxIR0l8uF3NxcVXrMLoEv37979+6NqQfolU1mOgbkPNfr9SZcHMPMXBO5DiCnU3TSptpw6W9JNSxqampQU1ODPXv2xPR2P/XUU7YkjI6QCyU57tfr9aKtrc1Ub5dej5veSIXovUxUqZR7zoDeTGTMmDGqhyQSicSEZ8n7FRmYaK2LCo5cQfX5fNi7d6+ql0+P2K+2gizOk1zYGD208rHp9eIBwNtvv41gMKgMRcsFhN53tm7dqpznlpYWFBUVKWns6elBVVWVUvjEy4jE+Zavm15vm15vlryNtjEXCoVUEwqFeCEf2sJejFhkZ2er0iEmGYp7Il5mKadTFOCprAoyWInrJ1Yni0QiaGhoUCqcDocDWVlZSVXoSktLVSN8gPHI5WAkOlREZVM8Jy0tLWhpaYHb7UZpaanq2TDbu6t97kRjvLa2VjVymIi8iISc12qfsXiVG5/Ph+bmZgSDQSXcS06jOKZ4HUxyJVtr69atmD9/Pny+3pVvxCiy2ZVv9J5/bQ+v0bFpO74EEfYlyqxoNIrOzk5V/q1XTmkbFIK2lx440iDXjrZoIwyAI+W43r5FPhrvmRPf0Vaee3p6lAp9vA4duU4hFnEQlfns7GzD7+nRjqLKxLOS6JrJnXTafYVCIaxbt04VyaBHLr+19QD5fT1y2WSmjBH7zc7O1n33j/bY9O4j0QAU91t3dzfC4XDMfSHOl53LLvcHyw2LO+64A3feeSdmz56NsWPHKsPpZD9t2JH2gRMteZEJi8xEr+fKKIMDjiw9CfTe2PF6K+VKpbby2dbWhrKyMuVhDYfDSqECxPbMyxm9UYYh98jEe1uqOFfiBUhCWVmZqrBxu92GPeJmRntEgSB6e7q7u2MWLZAzFvkcFRQUqHo/REEkRjQaGhoMKwWi4ggcuV7yJEmxnbanS5vJ6fVG6TUw5Ql0RsTfbGxsjInXljN3h8OBaDQaN7OUKxNynjJUKrVmyAsmyPeV3HiLRqMoLi5W7rN4566goEBVuZQL3kQ9noOReNa1vcCCXNkVWltbUV1dDcA4DwOOdKKIZ7ekpET1d0RlUKTDiJy/FRQUKHO/5HxfDt0wen7EXLFgMKhUhkUa29ralPwxmQVZ5HvT7Juu5bJOLuf0GkrxRlaMKpXaUGFRHgSDQWXVJ23eKZ9HQTQuAehWprVpkkdhd+3apYT3iA4APXLnoR65DMnNzYXD4VDqB9FoFNFoNGGHjnyetOmQ/742REuvQi2PfIjOTnFuRF4eL6zO7z/y4jjRYaIlRghFWvTSJJfV2hEL8f/l5eVKGLBg1LjQRorI4c5FRUWq8k/OG/RGm8Tf17u2TqdTuW7i/tCO0pgZBckklidvP/bYY1i5ciXeffddvPLKK1i9erXqH9lHZIBtbW3IysqKCbMJhUKqhzoUChmuDmI0advtdqOkpES1nr1opMSbVNjY2BhTIZULR5nRJC7xs8hwBHmClZw5xHtbpThX4gVIIkNrbm5WPeBZWVm6E/ASEWnKzc0F0JuhG036kifay7q7u3WXd5RDwORjaWo6ssqXfB60K6eIIfHq6uqYVZ60k9PE33e5XKoKvNvtRllZGcrKypQRpHgT1ES6RKaqXfrQ5/MpFVaR3nhD7HJhIh/fUKnUmmG0EpHX61X1PJtdIWjYsGGq/wLqRRMGUkFmJyuhINoGvhG/36/qwBGjIDJROYk3QbS8vFx5PoPBIIqLi1FVVaUsMwyoR22Nnh/5czGxVeQVQOKlucUEa6/Xqyo7AHU+LU/KjkeUN9pzIuL05Twx3oRbn8+HsrIyVZpaWlqUcwr0LuIhnheHw6E7H0z8TZnD4UBVVZXSK62lN+oknqfy8nKl7ADUeZx4+Z5IU6JzJf8NMcpQXFwccx0KCgpizpO8mInI63t6elBYWKiUB3JeInr39e5vvdGhlpYWVYMrHA6jtrbWcEQD6D1HYsEVbaNCr+NalDnxnjnt+RZpCwQCupV7vf3Ix679TmNjo3Iu5TQbzdmJd00jkUjMqmihUEi1ZLN8/w8Elkcsuru7ceqpp6YjLaQhhqbFzS0mC8sVem0l0+l0IicnRzeDEze6nHmLxoiWWDVEL8zKaOhT9CpoadOoRzxIpaWlSu9FfX29Em9ZUFAQt8dLO4wv/mYwGFTCqVpaWlBYWIiioqK4PbHyBGcRYy16YdxuN8466yzVtiKGVh4p8vl8SgYudHV1KedaXMfs7GwlRregoEDZF3Ck50ObeWvPp/Yayj1U2hEXUWnSi6vfu3evqrcpXnysdhRIey71MvB451t734gCcahVaq0qKysDoD+vKN4IkXYkEIBy/8q9fUMlHE3uEdQLczEj3v2t1xGjzcuFeD3N2rxXHgnQ9rAbLWkqjlWMWAHAli1blJ5rMfoRb66TnK+sWbNG+Vw7OlBUVIT29nYUFRXpnhu9Y5IlGvXWO67s7OyY8yrfz/I5ysrKUoUOinytvr5eqfCLtEWjUWV0qri4WDdfjjdSqB1NEfv1eDyorKxEbW0tgsEgdu/eHTe0WQ5lA45MKM7JyVGlKRgMxszDE8fX2NioXG9RfohjDQaDqK2tVe4BUT5p8xN5NEYv/xHk3xmthqU30gSoyzoRCiV3sMR75uR5GkDvvSl/1+FwYOzYsYZzJuSRFK1QKKS8cFHsS4wU6YXTtbe3q+4lMxKNgmQyyw2LK6+8Ei+88AJ+9rOfpSM9JNH2nMkxkdqfRYzkiBEjVKuZ6NEWaHIBI0SjUWRlZemunGG11ax9wGUul0s5DjFyIadNfFde5lWPOFctLS26k73kYX4xx6SxsRF79+6NCR3Tho21tLQoPSfaCXMiQ9FbLUKbWcrHJf5fnHePx4Pu7m5lP5WVlTET7IyIBqDIhOPFyBpl4OI4ZT09PYYVC3k/olES7x7RFnDyhD1txi3iUElNFKrhcFgZNpcrGFrxet6111quzMZbHGKwknsEKysrLb+nwuFwqJ4BsU8xz0ovvKOkpEQ15wrozQ/HjBmjW9kxCtESvZra/RhV4kSjSb7O0WhUFbcvh1wluv7ytvKkZrnCGa+BaqY8EWmKN1cuXrivEW14j5yv6XWUyedI5P2CqFRu3bpVSY+cd3766ae6aejs7FRWdAN6G1NbtmyJG5qjvZ+ys7Nx6NChmDRp96NdeEPv2MTxCaLBos1PjDrQ4jGaK+nz+Qwr8YJ4a/qYMWMSTujWe37lUDSg955va2tT5hKJdMgSpUfeF9Bbt5InzwPG837MEqutDaT3xFluWBw+fBi///3v8cYbb+C4446LeQv3/fffb1vihjqjSmAoFFLNZQCOVMDjxWeKgsnhcKgaFyJOV+b1epWefe3KGYkqu3LMoPy39UYctI2OcDgckzkCR94TYTQJUD5Xcvrcbrdq3Xa5IqsXh2rU0yHHsOuFJWlX+QBie5aMiLCptrY2Ve+/GIoV50O7WoZMvCcD0L8H5N48s0SDqaGhAfX19XC5XJg2bZqqt1K7Epe4R7SrxWjJE/ZkTqdTSf9QqNBaIRps8nmLd2/Fu9baeyk3NxfhcFh13RwOB3JzcwdML1kq5B7BeD3oRuR8QZ74btSIB448M9p5adrlw+VQID16FXPxXgIgtrJk5l0cch6XaJW2gwcPKv8vH7PZBmq8zg5t+BhwZLVBozl7RtcuNzcXXV1dqrKloaFBeR+MXr4WTzQaVa2IJzcM5JEPoPcaxOtgk/N5sW+teMcmd1BVVlYqo0ja/chzJM0co5m5VvGunx7tKJvYt3aFSkHu6ReLFojPzIzqCUaRE0Zzdqx0oMqRImbm/Zghd/gOtLLQcsPi3//+N44//ngAva1hGSdy20sO4dGrbIvP5IdIhNOIG1yvZ1hb6dd7mFtbW9Ha2gqXy6UaGjf7gGjT2tLSEjMyoPcyOaMeAnnilPivaGTIPTnav11cXKykWcQPy39DVLDEMG9DQ0NM2ltbW3XPv+zQoUMxk+CM4jm1enp6lBAk+R0O2oJSWzCJNIXDYZSUlCjnRoQN6K20pF2ZJl6PjMgsxXdEgSlncNrhXfndITJtT5VRZiunZ6CvjJEOZsIKBb1JmPJoh0yv0iMm3A+kAi1ZIq/V3pOJnpGsrCz09PTohlWYySv15loY9cQbEeGd8gRYADEVW8GoAudwOJCfn49gMKhMCpZHHYwqN0Z5nFias7CwEK2trejs7MRrr72mdE4IPp/PsNKsV6coLi7W7VzSa3jL9O5xvWWDxX8T9aCL7wO9vfqRSERZMlXO3+rr65VRjHiysrKUScIOh8PyOyDkTintwgyCmQ4mcc6j0Si2bNmCGTNmxPSWayd1WyHPN5BDs4wavHrXU3wmyghtWScTz6ienp4e5Xk5fPiw6pybbRC43W7lusnvhRF8Pl9SDYtDhw4p+c9AKwstNyzefPPNdKSDDBQVFani3gWjG1Vexs7M8CKg34Ml9+qLyddmhzwjkYhuRVxkIOLBT9QrKBfa8lCjdqWReIWJPGTa2dkJr9eL1tZW5fhEYdPa2qqKFdYej2DUwBDzYOSC3KjHQ7sPOaRN2zsY79jk3i157oYodOWVOOIdl17lSax+Iq+SI79Mz6giEO9e064Mk6jSZnb98qEkXq+nHlFQxruP4km0ms9gYBRiBMS/n4Ejo6PRaDRmpEFvLX0z5EqfmXDI1tZWVVrD4XDcCaXyPDZZNBpVyo/Dhw/D7Xar0p+dnR0zcmG0wENubq7SISbnt3qdE4B6qWNtmrT5pQjJ0VspJ5n5MfJbwuXOuETXXtbd3W0YwgaYGyUS518cr/Y8iZCdRHMa4o0umQkX04YLNTQ0xJxnOT/RW648EZFGOTTLyvkWxApN8SY3x3t+otGokqdqz3m8Bq9MvnftzC/lvyuHGA6EuYeWV4WivrV161bTmaS8RKhYNjaZh1VLZJRiZRMj8u+MelblVVTi7QvozRC8Xi+ysrLgdDqVt1xrV8+Ily7tuWtpaVGdE6/XC4/HY/o8mekxbmhowJo1awy3jbcP7YpVZsJQtKNP69atQ3V1tZKhhsNhZU1zvR5Avd4rMQclEAigqKgIHo8H06ZNA9DbqE2mkqrtqUqmwTvUaVc6SSTefaR3L2j3L3pOB7NkKv96xOpKYkGHeBOW45FHAOQKnFEeJ4d06hEjlXKly0x+p33Gg8GgMgoiVlky6rSQK2vavyU6J8TKV/JS4VoulysmvxQrRYmVpOQ0+Hw+VFVVxYRoG5GX7hbnJ5lGeDgcVjrg9Ea+jMij3NqX9+ktwOL3+w2vtZjoHm/1oGR6vuW3lItrJu/Hah1DrFonzlFlZaXpl6jqqa+vR2FhoTLXMC8vL6n9yB1nNTU1CVezkv++7LXXXovJM80s3x6PdpnlTJfUC/Ko75itXIkeZnGTR6PRpIbf9Mgt+HjD8mYyGBGXarSClJZ2knQkEolZ7UL0LFjN4OT5GnKvZby5DGZoR0PMCoVCSgYrMrimpibL6dGeW5fLpfTyaQtpvYn78j7keHHRs5QsUdESvbCJwsvsaBQPJrW1tXHvA73zKb9hVrzTRUx+zcrKipk8r92/Xs8pGdM+L4mIzhJ5nXw5DEY7mTgVhYWFqknCevTuIb2wKHm+g1XySI58vsQL6+T8qLi4WDXiARypQIuGsfidPGncbLkp71dektvqZFt5qVkrDVWj/FRUJAGoRofilenysYjGlzjXoqyz+lZtmZllla0SIbLJNO7le1Ucl97SsGaI1fXk/DHZCdd6o3LalwebIeZtFhYWKvc1MDCWYOeIRYZLNG8lKytLmZQsD0nKhYPeWuNmuVwupfVtdpWiRPuz0nqX4z1lInMTvQvJpEuEEADq9eGTbVS43e6kClpZS0uLKka7s7MzpUYO0DtUa5QZdXR0xP1uYWGh6r0cyYxUCKJnUfQIWpkvQIkrLHrnU+6NFWvqix498fbZRAZafK8VevO8UiFWhDF6z4KWiMuXK0OiU0jkuyJfSfR+g3jKysrQ3t6e8PnVu4dE/lNUVKRUwOJtLxi9AE6+j8U7FARtpVDbqND726JsE51e2lFpM8TblLXvTjJLbsgkW9bKRN4vRoeqq6stdRSKFZTkhqB4V4nb7TasV5gZERXXLNWyTkh2xHDs2LHKM5FqHiXCBUVjQu+dX1aI0DoASd1TDocDc+bMQWVlJdra2pT3mem9DywTsWGR4YwyZ0G0rrW91Lm5uUqB1NbWlnRvl2h9A8kNEWvl5eVZ2ofehHWn04muri6sW7dOyYCTqaRGIhHVw5/KsTkcDpSWlqY0pCuITC4VouBwOp3KSiB6GWWi87Z79+6EQ/BmicJXVLy0Ewv1DPYwHCvMnC89cgw5kPjNvlqiUjIY2REGJVfSgsEgmpub0dnZib1795qqnBjlzVu3blW9DExe9toKp9OJxsZGSyvCycRL5ETFVm5cxOvQMdPZEwqF4s6lSlRuyW/CTqUnV/R4W11mWFZfX4/q6mqMGDEi6X3o0b50ziyxwqIgrl9paalhvm+mE8vo3VfJSCVESF5aXvzX7XZbDhcFeo+psLAQHo8HeXl5qgZGMuS5jcncU9FodEDnuWxYZDirKy4Ihw8fVi1BmmwvisPhiHlrcyqSGaaURaNRZVWTRIWOmUJYzCNINWwsGo2isbERu3btSmk/wJGYaJfLlfR1Ew1SkXGn0hOnXXUrUWPXiEiL6DnXW+ZYayDEk/YVM+dLS7zjRiws8Pbbbyf1t1OpcGWyVHolgd5nQVtJk1ewS2UBAr1wnmQ7UJINEdH+TbH4hh0rQIo3Q+uxcl3s6MG1K7wklfOsx+VypXSPRqNR1chCU1OTrXlqqqMW4XA4qYaANg2i8Z2VlZX0CH9raysqKytVz28y4VB2rY5aX1+PNWvWKOVmXl6eMscl07FhkYHkiW3JEjGxHo8HxcXFSY9YiJdw1dTUqEKHUn14Uh0ulodzjXo9zBbCdlWa5JfU2UHv5XFmlZaWwuPxKA3TVAoT7WTtZMPh5BEisxPjBkI8aV9JZjJ7T0+P6nrZWekZDJKdYC3oPQtyRSnVEd7i4mJbQmushqDGY1cYY3d3t+HEXbOdIPL8imQ7h8TLDaurq21bMMKuymU4HEZ5eXlK90BxcbESvSBWYrKD1egDPQUFBSmF+hYXF6s6usT7MJIhyqdUOxscDoeq7E01ikF+Rxknb1PS5JfrJNsjEI1GlaXvUu2N16uMpFKwJFpmNJG8vDzMmTMHNTU16OzsTKkwsKsx4HK5bF/FKJV9imve3Nwcd43wviRWmgLMVRzcbveAiCfNZHoTcZMpyO2qlGaaZJbKTMTOt8a3tbWp3iWTrHA4DK/Xq7t0uVWp5t+CWFo7lfMVDofjLhdshvalp9p38yS7T5nH40n6GtbU1KRURskT7c0sPx+P3kT+ZLndbkvXXu99FKKOI5Zh3bt3b0plXVNTEyorK1M6NpfLpXrvx+7du5Pel0z7zq1M1q8jFhs3bsTZZ58Nr9cLh8OBV155JeF3NmzYgPLycuUthytXrozZ5uGHH8bEiRORm5uLk08+Ge+99579iU8juUch2d5h0TthR+u2oKAAbrfblp4zILVGCdDb0FmzZg2ys7NTHoq1a4RBzsxSHdrV22eygsGg6Um6fcXsPZlsGOBglWzlXv5esr2D4uVRg006ngu792lXD3NbW1vKz1QoFLKtkSmW1k41n7Mr3l8oLi6OmaRuVUFBgVJeOp3OlBqGyc4h1Er2XREyUbG145yHw2FLlWS9upCYyyQmNaeaLrEaYyqKi4sRCARQXV2NtWvX2tapJ/LugfB+p35tWBw8eBAzZ87Eww8/bGr7bdu2YeHChZg7dy7q6uqwZMkSXHnllaiurla2+fOf/4ylS5fitttuQ21tLWbOnImqqirs2bMnXYdhO5/Ph7KyMrjd7qQzlOLiYmWJMrMNAqMK8bBhw5CVlYVIJGJL48KuynwwGEw63h9I35viU13FKZOI8+twOOD1euOuKGKGHEYxWHvB0yWZ9dmdTifGjBljy98fCEPwVqUrD7BLJBJBc3OzLfvq6emxZQnyUCiU8nnzer3KIg6Zlg+0trbaMsovJnFnSqeOXeWuXQ25SCRiSyVZTk+yC1wI9fX1KV97cf/YvfKhXR27fcERzZA1Hx0OB1avXo1FixYZbrNs2TKsXbsWW7ZsUT674IIL8NVXX+H1118HAJx88sk48cQT8dBDDwHovXnHjx+Pa665Bj/5yU9MpaWjo0NZTzs/Pz/5g0qRCPXJBF6vV1mVJENuGRQUFGDChAm2va+D9LndbmWd/VTf8WGFx+NRDSkPdUZvho/HrrAVQP3el8HitddeszVMMJPyx0wlRgNSDcuRDdbz3pf5bX+w67qJ5Z3jvWdrIHM4HJgxY0a/vnnbSr144DSBAGzatAnz5s1TfVZVVYVNmzYB6I1t3bx5s2obp9OJefPmKdvo6erqQkdHh+pfJkh2ecB0EHGCmZR5i6UdKb3klU7sKOTMhopl0r2WCZKZVGhnb+lAGIK3yu7e5Gg0mnKv6WBXX1+PrVu32nruB2teIS8bPhjZdd3Em8G7urps2V+mEe+2sXP+VjoNqLu1tbU1Zlh/zJgx6OjoQGdnJ7788kuEw2HdbeQVjbRWrFiBgoIC5d/48ePTkn6rMmkFl0yY/KvHznM0GDNvs5VRp9PZZ8dv9l4arHH9yRpsowWZwGzFxkroTybl25kqU8sTIVNC5MQk8kwJp8pkmTaXMB3k94plssFXk0rC8uXLEQwGlX87duzo7yQBUGduyWR0LpcrYzLITDdYz1W8BrUsEon02fFbWXFlIGSifaU/G1niBZBD1WDMG/pSWVmZbsdFpp7XTBkBSXbxFhq8BsIS7ANqudmSkpKYpbt2796N/Px8ZRKYy+XS3aakpMRwvzk5OcjJyUlLmpNVW1urytySyegyvVcIsDcGPBV2natMOR7BSloy8X4ZCJloX0nH0qhmRaPRQbn0r94Slnoy6ZkeiOrr65X4bFl/VOAzLY8mMsvr9Q6IfHhAjVhUVFSgpqZG9dn69etRUVEBoHdOwgknnKDaJhKJoKamRtlmoLB7Cb1MNdAzeG089UA9HofDoVoekTKP0b2Vm5ub8kudhir2CPedTAkRG6h59FDCMiiW0+kcMOGw/Xr1Dhw4gLq6OtTV1QHoXU62rq5OmZC7fPlyXHzxxcr2ixcvxueff46bb74ZjY2NeOSRR/Diiy/i+uuvV7ZZunQp/vCHP+CZZ57Bxx9/jB/+8Ic4ePAgLrvssj49tlTxwbJfOobdM6WwTJbL5YLX61WWlU1lieN0YCjUEUaNh8OHD6e9I8Kud7PQ4JaOPDZd4VKZGoZFbPzpiUQiqK2t7e9kmNKvoVAffPAB5s6dq/y8dOlSAMAll1yClStXYteuXapVfyZNmoS1a9fi+uuvx4MPPohx48bhiSeeQFVVlbLN+eefj7a2Nvz85z9Ha2srjj/+eLz++uu2reXeV0aMGDHgK62ZxOl0ZlSFOVNkZ2er3sbrcrky6jwNhLeM9pWioiK0trbqFrrpDu/IxDC5dHI6nXC5XLrzgTJhaVO9sKJMkJWVlfJbq7XSda77+xqmC0O9Bq+WlhYUFRVlfDhUxrzHIpNkwnsstGvWizdM251py4wKK6sZVW5uLsLhMHp6ekxl3m63O63HRcZcLhcikYjqpXVWK5Fut9v0tbaK77I4orq6Wvc5EZXgcDictgqFw+HAwoUL07Lv/vT222+brqC7XC6MGTNGWXY3mTxLr1Fidp6HrKysDI2Njbbkm6k0lERHhLjv0pkXZILB/l4Jii8TOhX6q0wctO+xGMpCoVDaK98HDhzQ/dxqZSUcDqO0tNT0A1hcXAyPx6M0ntKhL5dTHUjC4bByndxud1LnKBKJYMaMGXYnDQBHLGTa51CEJzkcjrQvtdjfhWm6WFkX3ul0oqWlJaV8WO88Wm1U5ObmKm/2TZXb7VbCIJMhGrMirMjutw1nApEnikYFQ6j6Ryr3qZbT6bRc3ygrK8OMGTP6tR7hcrkGxIImrGkNcHbGPofDYXg8HtMZZ0FBge4yguFwGI2Njab/bktLC/x+f9qWs3S73SgpKYHL5UrL/gcyt9sNr9cLj8eD0tLSpK5BOBxGfX29rRm/0N7ebvs+ByrtcyYqlqIi53A4Uq4omv3bg4WZhmtBQQE8Ho9qJC8SifTbhHk7e8zD4bAtDRQ7GxPprLgncx+L9Ijznsyx2vn8FBQU2P48Op1OeDyejH65o9UGuCjXysrKYu6pSCRi6b53Op1oamrC3r17YzpwzjrrLEvXw+12K+e5oKDAUgMnEolkfBgUwIZFRrKyXr2dhYzX60VhYaHpjLO7uxs+ny+mwu5wOCz3ntbX12Pv3r2WviOL92BnZWWp5hHYQS+zMiKvuKT93GohKmdKqXK5XCgtLUV7ezv8fj98Ph98Ph/KysqSGj3SZvypVBDcbjc8Hs+A6J3pK8XFxaqfw+Gw6jpFo9GkwmrMGDFihO37zASJ3ibu8XgwYcIEAOr7ORwOo6ioCGeddZalv2fHs2tHxVvsIxNj8dM54uFyuSznbeI5S6UTT5SRogKfigMHDth+3aZPn47CwsK0zNtJZiTcjuekqKhICRmKd0+Zua4ulwudnZ0xi2QUFBSgpqYm7usMZE6nE1lZWdi/fz+A3vAiK/WSgTIayIZFBoq3Eo4d4ULagkn0yBUVFVlaXcbv9yMQCMQ8GOFwOKnJnqmsbKPNaOWMrLOz09bKlsvlgs/nM90znJWVhenTp8f0cEajUcsZRSQSwYQJE5Ku/MtEL0xnZyeampoQCARQXV2NxsZGlJaWJiwMEv0+Go2irKwMZWVlcbcTPe1yA7W4uBiVlZUDonemr+iN3mjvwcLCwrSMLmTiRGE7JKqg+f1+5RnRhlOKfNrsyIXL5cKcOXNQVlam6h222lAw2t7KCIpevpNKhdfOUQav16t7Dydb2ZTzydLS0qRG9LKysmLKtGTy30gkkrCzJNG51HYoJJIoPxD7SnVlOaN0FxcXY8GCBZbur+7ubuV6J5uf1dfXo7a2NuHKgmY6Z7XXvqysDF6vF8FgEJ2dndi1a5epZyAajaKzs1PJd+TnUHx/MKzAx4ZFBvL7/boPk9frRVVVVcKKWjxOpxNjx45VfnY4HMrD0djYaKlw8vl8qodWL+NP57wJI263O6bCYLUCHy/d06ZNAxDbg2wkFAqhqalJt2JoNZwiHA4r5zyVlXrE/VVYWKiMDDQ1NSlzeZqamjB9+nTD77tcLkyfPh1nnXVWzP2orXz5fL64lQJxbeTvJepJHoq0FRK32w2/3680yNxuN9rb2y31ZmZy6ENfiJcvlJWVwefzwe/3K6GC8otURRhVeXm5qfOoXZlQzGnTjmjGq1iIF8HqSaZiKKc7Ozsbbrc7qUaCeEmtzGqF0OVyoaysDOXl5bp5j9y4Nbtvt9ut6iTZu3dvUvO2tM+e2J/Z/Fvu2GpqajL8nsvlwowZM+KOhGmvT6I0JMoPSktLdSvfZstut9uNsrIyZGVlKXMX5DS1tLQgEAiYGn0W3yssLMShQ4eU9Cd7X7a0tCTsVDSz30gkopprU19fr3rexGhxInJZpz2/4nfxGjoDZX4PGxYZyOfz6RZ44uUoPp9PeQBdLpelm83lcqkquNq/U15ebqkxIGcW3d3dSk+6eHDM9hCJURM7elvtCHky6rF3u93Yu3cvampqYt7wrkdUGnp6epRKvBxf2d7ebrlxIRoBcoFhpbdDjA6EQiG0tbUpw8U9PT3KdSssLIw7T0bMq6iursbevXuVcyUaHKJXVtwfZnq8zTbUhioRpiZkZWXB5/Mp5624uFipBJuZT5Sbm4s5c+aYet4HSoFmlVHDwuFwKKN4onHR0NCAzs5OZZv29nYEAgHU1NSYur9FvitGQETHQCQSUT3LiXpQxWIXqfZsut1uJcwL6H1Gq6qqklr9KxgMKveeSJe8KISZ3mq5o8Tn8xnelwUFBaYazyLU0+fzKdu3tLRYnrclGpilpaXweDzwer1KjH5bW5uS18Ur56LRqJJHZmdno7W1VXc7p9OJhoYGrF27Nu6+xLnyer0oKipKugNPPNci7xc/a/OPeOWyKAtCoRBycnJQVVWFoqIi1TaNjY1oamqKe8+K8hCIHZ1NdlEA7ZLRenNJzK5cKRg9n2brHR6PB9OnT9e9X8R7pYwMlLluAyOVQ5D2ZhcPvCjIRI9uOBw2XeiLir5ej434HQBVBmpEPJhyTH5XV5cSRrNgwQJUVVWpMvV4RMGUTAys3USok166Q6EQWlpaVBUDI263GwsXLkROTg5CoRDa29tRWVmp9MSIkSKzBZ3oGRIVHXHdxRJ4Ho8HxxxzTMJ9lJaWKmkX/xWjFaJgMDsnRZyPSCQCj8eDadOmKb1fIpQpEAgkvEfFXA95vwPlZUB9STxvcqNNnLe2tjY0NTXB7/cro2rxiAKytLRU1Rmg99wPtYaFCFkQ97K2IQ+ow6TMENcrOzs76fSK/KKyshLz5s2LW2EXlRSja6ftqc7NzcWaNWtQW1urugfMVGa8Xi98Ph8qKytVk5w9Hg+Ki4tNh6LK6ZE7GuSKoNmwPKfTqYRSimNwOp2mes5Fh5BcFonj0+bX4vOpU6call0ul0uprAeDQd2yRVSC5eW/jYjvt7e3K3m3/LfNPq/RaFT5vvx35XBmr9cbNz3ysYhzK19Hka7Ozs64jeZgMKh0vok6ivbes1I38Hq9mDZtmipvmz59OiorKzFnzhzT+3E4HKiqqtKdS2pVQUGBUi76/f6YkZjs7GwUFRUZPnNsWJCtsrKyEAgEUF9fH1OQmV3tSFQUtRmjCLESlUBROdH2OsiCwaAyyVzMNxC9OPX19aoJ6GYaFg6HQ+n1yMrKUsUfG8Xc6hFx0CIjETH+8nK2iTIEUSkTacjNzdVtaCUKf+jp6VGGgPUmIot0mK1oiN43oPeci96frKwsVThTPKFQCHv37lXuGfFfOY1682aMyKtKyRUtOR1NTU0JC0o5oxXEEDqpiYqMeF5Fb2MoFFLOvXZ0Q4/oMZOfX0A/pGawrqgWr/NEfmbl8FQRsiOHSSUitgdiK8ZGizsYkTuGxN8XMd8ifUBvvnLw4MG4z574vtfrVSp9LS0tqpFro/xbzmPFaLr890VazS7Rq80jReeZ2+1WzU1JlH9rf19bW6scQ0lJialnQ7ybQ4SFAkc69QoLC+F2u1WdccCRzhk906ZNUxrwcvrEvSRWV5RHDRIRFXDx/MuNN71rrjda4PV6VfeAXAkX17G9vT3uyJ58HOIel/eZlZWljLIlIhrN7e3tSkeXXA5biUZoaWlR0hOJRJQwYlGmmH3mxLHL19rj8WDGjBnKMZq9ZvLy1j6fD1VVVarrUlhYiPr6esNnbqCM6vfrm7fJHI/Hg+zsbNTX1yufFRcXY9++fTh8+LDhG2L1NDU1obCwEJ2dncpLteS3OYqKoZmXL4kKDNCbkcjfkX/n9Xrjxv9qX5AnV5Cj0Sja29uxYMECAMC6desMHzqHw6FspyU3mnp6ehAKheDxeNDV1aXaX0FBAerr67F161Yl43E4HEq4UFFRkXIOE8U0i94geRKyyNTkFx4a9cBpz1tjY6Oyn9raWuUaigaHqGQm0tLSorxgS6RJrAgF9L6ITcvoxVDyW+8FcX5qamrg9/vh9/tV966WPPqlLZzl+2iokxv94pzoVWYKCwuVaxjvhU7yvSKeOaP7J11LQfe38vJylJeXx7yQtKCgIKYSoHcfis+NXrTncDgwY8YM1Xe1z7XVRRxECJa4F0Te5PP5UF5ejkAggMbGRhw+fDjuaEp9fb3yckU5NMdsIzISiSAnJyfmvEybNk153q3M+zB66Vc4HFbu5+Li4oT7dDgcyM3NVRop8vaiU02kuaGhwbA8EXmeaMiJ8q2npwdVVVUIBAJoaGhAfX29UkHfsmWL7rVsbGxEVVWV6nmVX3RWU1OjlEkADK+bSFNBQQHmzJmjfE+PqDeI+1Kbf7vdblV4tZyOnJwcJU8oLCyMKScFec6AfB+I56KmpkY5FlHviEduyGvzupqaGlXaE9VPtB1/kUhE1fEyffr0uOWSIOoBcp6rzYMTvXFePGd6o2XysSbqGDQKocs0HLHIUHLPrd/vjym0WlpalIxCNC7M8Pv9Suaqt8KJ2R444EiPPHCk9a0N0wAQd+QjNzdX6cnRhvqIz/x+P2pra7FmzZq4ox9OpzNhD3dPT49Seeru7o7ZnzjPYhhYeyxGw+F6HA6H8l3R2yUKp6ysLFXImV6PR7wCVPwuHA7HrWTqEWELIjNMlJl5PB7d/er1gInz09bWpjRQEzUM5DhvuaeLy82q6Y0Gac9XWVmZ0tuXKC5Z7q0X101uQIh7Utwvg5U2z/B6veju7o4512LbmpqamO+IHnUtMQ9GVl5erpqLJoegxetFlZdg1rsXRNrE3zXTWBG98nI+KEZszTQmu7u7Y86FPAdQrNNvJoRFux+RP4qOMxF2mYjL5VJ15sgjOWJEVpyneOWJKF9FXi+2Ff+Vw+NaWlqwd+/emHOuvZ7ack3+XFzbeJPLw+EwysrKlFAe+fmX/5Y41nhhY8XFxaipqUFtba1yT8vpkMu6RFEH4XBY97mQ95fo2sn5jDwqK+9LnLu8vLy4+wKOjA7I5axcpsj3qRExZ1DQpks7X8pIJBLRzQu0+5RHII32MxBwxCIDBQIB1Y0qT6KVeyDF/4uWeaLYU4fDEdPq1v6/6Gmora1FS0uLMqqhty9RKdXrqRDEfowcPnxYtzdQ+5mZngXxYj6jSpC24q09LofDgbFjxyrHLYceacm9DEa9Xrm5uTEZkFw50Dtubc8pcOTtuHJBJHo95YxR7mES59zpdKrS5vF4lF4q+RhkpaWlaGxsVN6qa9TLZOUdKvFGreR7yKhnmPSvl9H5SvS8OJ1OVfiKvD+g95p0d3cjHA6jtbVVGX0ajNdGr2Hd09MTU/kT28q9noIIU9XSyzvFaIN4zkKhEFwul+7IiUybH2nvBTlt8fKlgoICHDp0SDV5XJDDWcwQ4SXa74jKeHd3tzKqaTSqIxjl3cXFxUpYVHFxMVpbW5W3fes1niKRiDIKC/R2bMnv6qmurkYoFEJjYyMKCgoM0yT2L8f7h8NhpQIvn+OCgoK45ZwIYdF7XuXRJyB+h5K2zNXuT+xLhNTocblcyly4zs5OJX/Xjq4LYqRBW5YAR0b2enp6YuYkacu4vXv3oqWlxfCcx1sJUByXOAdm5tlo6zTyvsQIWGlpKYqKigzPVaIGgzZSw0x64pHrXwMZGxYZSDtRUL5pZ8yYoWwjF/RygaSXAQC9w5basBcAupm5KBiMHiwxMS/Rw2Kml8KMROFUZshDsQ6HI2b4UlS25AqXXvgJcCQDCAQCysQ87TC9fG7E3y4uLtat0MU7Tr3vaNMppwnoLUz1MjwxyiS2NepBMTP64XK5lB4q7TkSE0PFOSgvL9e9fnqVN9In33PxKvo+ny9hYRdvhFP8nerqalXjcrCGpWlD9cR96vF4Yo7XqDFuNOqnd55FhU6EIsmM8rmsrCzlbxg9u9rwEXGvyMfm8XhiJq4a5XHxRjK9Xq9SEYwX3iGHRMqhZWZoQ48EEV5jNCIjVikCjuRlRvevnCZtQ0X8vygLtXmafA3kMB2gt1PJ4XCgp6cHkUgk7gi33ugTAMMKuDyiITckROOpsrIyJj0y7Si3/F09Iu165/vQoUNKw0/k5UbnW5RZeqG2QG/ZtGbNGni93piyTdwLYlUnEeIV74WgZjoYGxsbE65aGa+z0ufzoaGhIe73E6VHj1Fdx+oKkv2FoVAZyOgBF71JesOE4obzer1x3wIZbwlRM2mQmXmBmUiXXqiPPOnPKMRAEOED8UIFxFrzRuQekdzc3JjhVL19G2X68u9FaJO8f+3Qv9Eyelrl5eUx58zKOx0CgQDWrVunTISXiQm+Zu4BMyFxordS7xzJ96gceqAlJpOTeYnuSUB/uWR5pRvtc6L3/BmFEAw2Pt+Rle20ixHobauX7xmdG73JlvK2LpdLlW8ZdTpoe4SB2Gumlzbt6jxWj0kvzAborfgUFxcbVspEWIdYQU9UyLWThOVQsHh5t3ysZt9DoQ3vFccu7uvS0lLV+dCrOMvnzOg8yX9D3DujRo0CcGRp4HjPjpw+ORRmzpw5uuGmchki8gL5PIt9CkYTi8XxlJeXxy3LRZr0zo8IT5PnVsrvRtKKtzCI2H+8DkQR5jZnzhycddZZho2KeBVwbbrkhRkEK28LN7v0sZWFSIxW5osXVp5JOGKRgXw+X0wvU6IwBLkHO95wupU0ADDs+TT70Il0aYfBtXHb8qRxvd4zkSaj5R1FhmOWPKlNCIVCeOONNzBv3jzlM6MeSr3fyxV27ZB1ov3IREMimXW79ZbFBKyv6mPU4yloJ0gaHZvR92ViFRqxvbYXbTD2kqfCzL0kzpl87qPRKHJzc3WfE71exqESliYmOwOx4UZmGeWX8sIY2m2N8rl4xJyGeD3x8giEPLlWrNi3d+/euKOmcjrl/WrLFVEBNOrRlRs1hYWFqnME9K7Q1NbWZjgqLo8QmGlMC9qJsnoj9PLPW7du1U1DonBYmTzKJ1ZnAxAz2hLvu7LGxkbdUDbtBGB5ZKi1tRWdnZ2ora1Vrq9eeWl1mVSRvrVr18aUSW63G8XFxWhvb1cav4DxRHxtp5aYjO5yuTB8+HAEg0HdSrV2tEjQG+FLFM6n9/z5fD7V/S2W7xWT2OOR02AULWIUMmhEHt0R+YlYDGYg5MlsWGQoMQwqVn8wS69V7Ha7EQ6HlYzSLKOQmES9S3q0lXhtL7y8Ko22sJQLyuzsbN2GhZlKeGlpqVLRktMjDzlr5w0kqlzJhWU4HFZltGZi4QW5Qq0931bOtTa2Wgzvi7RlZ2db6nnWVgbEsn1ZWVmqwiPesHMickGiF/c7EDLSvmIUtqJH79xnZ2frhlEZNVas/L2BShsakezx+nw+bNmyRXf/ep0kyZxPuYISLyxL5KEiP5MrO3JD3gxxD8QLO9Ejh5jJLwcU51quEEYikYTz9bTzAkXng1weWJ0jAhyZO6GVaM6elpVluuOJF4Y6f/581c/yORIVY3F9xe+0nXpmXlirFQgElO/Jc17khpd2rogZ4v4Mh8MIBoOG10/veQkEAkoIVjQaRU9Pj2peYyI9PT1Kuebz+WIaBZFIxFQ+UF5erqwW6ff7lbkkWsmM+JaWlirludGqUpmIDYsMNWHCBDQ2NuLQoUOqiWha2oJfrzIhekzkCVBmHj6j5UuNVjeIJ94cCbnHUK9SLheUer0v0WjUVI+8Nq5UZI7Dhg1TMt5kMl25IIhEIigqKtItuOP1jGor1ILVFXm0GbDc45Fqb6yZeFwtUfnRGyEC1JPJxfZyL/tAyUj7Srx4cS3tuQSONKiNGu7afVr5ewOVdrQxlePVdnCYCSHTnn+9vFJ+h4bcEw8gJk8Xx5NMmIgebb4rFpNIlBdoe4ZF2vRGMBPNtdIbdQgEArYsvyl3OKVCr+xNZplm+frJ95PD4YhbF9Bb0AM4Mo9EbxEQs+SRuKKiIuUFqnLDy0xjWYw8iHtHlMN6y9QL8eYAaRtgZkdj5O+Kv+lyuVQNCyujDPKx64UaJ7uy3kAdNWbDIkPp3fiC/KBpC369jFv7/gazD4v2wRUterMxrjKxModeYST/HfGCHJk83CvW0xarVYmM18zyc0DshDsRnyok8zZN7coQ8sRBWbxrqt2H2+1OOITeV+Rwps7OTku9gtqMUawSJibPawu5VMJEhgKzIXXimddOSBWruBg13LWNDfGsD5UGnl7HhhXy+dabhKpHe/5FD6icH4g322vzRqMQNvG53AtrNj1a2jLFSgeFXoNAG+LjdDqT6tBpbGyMCTtJpkHo8/nQ3Nys2/GR7GixyN+SIc6Z/A4IIHEojN6CHsCRhUNCoRCKi4tTzlPNjEIb0d4PgUAA7e3tcZ87o84N+X1cYqK82fBhuTNBnnuj7XRMJh/QRgsAwK5duyzvB4gtL5PtHOxrnLydofx+/fWugdhlBbVrM8uTnL1er6rHycokTO12otJt5h0OWiINouGgrVQbHas4JvFuhEgkArfbray1LphZfk7LrsqSz9f7Dg95wrZe5hvvOH0+n6ogsuuFZPJ+tG9EN0t7LMkULGLypWjERaNR5W3vWuJ6D4QMtK+ZPTfieRNvZRf3ZjAY1A2D0uYL8Z7XwUY0+EOhEFpbW1M63rFjxwLozXeLioriLkgh6J1/oHf0VEx0Nsobjb4rGoRyxTuZfBtQr/fvcrmwd+9e3eNKtACHII8uizc8m50/YaSgoCClBQa05Yf2TdJmiF5voDd/M/OOoHi0k92TWUEvEAioOs5aWlqSKgPEIgSiAi8mpadSTokOq0T5jNE9Lu7n7OxsZXTO7BLooszWK4OKiooMf2eGqHvprTBmlVxepno/9SWOWGSoeENgcq+l0XZi0pHP51Pe0pxoqVO9NABHeuDFw5HMiIWcZu3QptXhPjFEKb8NOpll2Hy+I8tyOp1OUxO14hHr0gP6jZZEx6m9rkDyMe7y98rKypRraCVmWJsuq2FQMm04hdfrHRLx+/1Fe83kpZAbGhpMhVFZjZceqOQeeSsvoNK7f+XV30T8f6JedL18QTwvPT09ysie2dhzOR1Cqss6l5eXK8cj7iPtcZkJm9NWakOhELxeb1L5ijyhN9X8Q4zsOJ1OLFiwIOn9yOFnqa6kpte7byWUGdDvBDIaUY9HXnJWrASVzHwWQbuoh5UwOPk7Io87fPiw6p1eZv6+9t5JNexT7FPMQZFHUPRW9zJDGxYpL96QydiwyGCJ3qFg9B3xwGonziXTYyU/dGbW5I63H20IktUMTn6hlHYuQjJD/PKLCEtKSpLahyzVeEj5+8mGrwkikxTrdBu96NBqupIlV97EvApxHwzm+P3+Iq6ZmGMjV3j0wkf0wmnE7+SfByM5PMfsandA7HkTc9LkSnyyjTM5NDKZ518O9TCaYGv1morQE71wOvE3Ex2vXly8XvirGXbGn0+fPt2WhrTcMy32lezLJfXmT1rNL+XKd6K5DIn2ow2xTiXfluchpDL/AFCvfDdhwgRT39U7l6l2psidZ9p3ZCUTYg2oG/SA9dWl+gtDoTKYfPNb+Y6svr4+7rrSZojwC7H+d6qZrzbkwko69OJWxaQ2q5qajizNmmyYQLqIay+WvbN6zsXwMQAl1jrZoWuzIQ7xiBA9+f4xGuKm9BCFm7aQM7oOyeQ/A9X06dPh8Xgwffp009/RnjdRaRaLW6Qa0idXTKw+I0ahHsleUxHrDfT2muodV6LjlRte8j2YyvNvR94E2Bd+qV3hLpVnqLGxUekcAqznl3LDpL29XSlLkslz5fxbhGaJl62mKpWyV3teGxoakg4/TPUekPep1+hOlt/v1128IZNxxCKDJdOClnsoREGQbI+Qll09RHJvaDIVZu3xJbu+szz5K9MeVrmHKJlVuMT2jY2NyvLAyV67VIeI5TSZHXkje8jhIqIHXNs4TxRqkGnPRjokcy9qv2Pn+ZIrTKmEnGglm0Y5HCOVHt1QKASPx6MsvOF0OlM6NrvyJlkqozrapUdFGu24J6zeo9q5mKmGjelFHSR73uVlVJMJrRbkfA3oHY1NNvwwVUYhjamG6smjVdnZ2QOizGTDIoOlWthpM7hk2fECqXjpTPZ7YuUUILnCTp78lWkPayqNL0EU5G63G42NjWhsbEzq+onMW/RQpRJ/yvkUfUv7nFm5n9jws8bO82VHRVBPsmmUlzJNpUdXbuSK1f1SkY7Gr52NlVTuCaOXwpklOs6ys7Ntv5dSPe9yZTmVEQttOJQdcyXtYmd+IN7flZ2dbcv+0s0RTXa6+iDW0dGhvDQtPz+/v5PT7+Rl7/SWPByoBntlVztPA0j++ol7oL++T0SDh/ZlapmUD6daLmRKXqddrjZeevqjLLTrb8rHaefonh1ECGGySz0L8lvBzzrrLDuSZpmVejHnWFBCfn/85WCtsCsm1o792RVXaze744ZLS0tTvn6pzofgfIqBye7nlSjepORMYGesfX8S6fB6vQnT0x/XwK7yVxxnpjUqgCMhhEYvBzZLzN1JZvXL/tDvDYuHH34YEydORG5uLk4++WS89957htuGQiHceeedmDx5MnJzczFz5ky8/vrrqm3279+PJUuWwOfzwePx4NRTT8X777+f7sOwXSYV6EYTAZNhdwaWyv4y6RzL7DhH8rHZcf1SLQQytRFH8WVapY/MydS8DYi9pzKlIm6XTMnrRDrKy8sTpmcgX4NMOd967GoQFBUVwePxoKioyI5kpV2/Niz+/Oc/Y+nSpbjttttQW1uLmTNnoqqqCnv27NHd/qc//Skef/xx/O53v8PWrVuxePFi/Nd//Rc+/PBDZZsrr7wS69evx3PPPYf6+nqceeaZmDdvHnbu3NlXh2WLwVqg252BiYlfyUwAy9RzbMc5kpebzdQKBvWdZCuaA7nCMZRlat4GxN5TmVwx7C993TDkNUiP8vJynHXWWSkvZZ/Jz7Oefp1jcfLJJ+PEE0/EQw89BKB3Rv/48eNxzTXX4Cc/+UnM9l6vF7feeiuuvvpq5bPzzjsPHo8Hf/zjH9HZ2YkRI0bgr3/9KxYuXKhsc8IJJ2D+/Pn4xS9+YSpdmTDHYrDH/9sllXjWTD3HdqRLO7+iv+N9qX9lStw39Y1MzdvIHLufV94PA1smXD8r9eJ+WxWqu7sbmzdvxvLly5XPnE4n5s2bh02bNul+p6urK+YNhh6PB++88w6A3rdehsPhuNsMFFyRxZxUVqfI1HMsls9L5i3Zgna5WfY4D21DaelYyty8jcyx+3lNx7K81HcG2vPcbw2LL7/8EuFwGGPGjFF9PmbMGNVbGWVVVVW4//77cdppp2Hy5MmoqanByy+/rLxReMSIEaioqMD//M//YOrUqRgzZgz+9Kc/YdOmTXEf0K6uLnR1dSk/d3R02HCE1BfseOAyoTcgHeR14+04rsF6noaCgVYwEQ1ldj+v7FigvtTvk7etePDBB3HMMcegtLQU2dnZ+PGPf4zLLrtMeSshADz33HOIRqM4+uijkZOTg9/+9rf43ve+p9pGa8WKFSgoKFD+jR8/vi8OhzJEpsUvijecJ/OWbJnd8fGZdp6IiCgxzqGgvtRvDYujjjoKLpcLu3fvVn2+e/dulJSU6H6nuLgYr7zyCg4ePKi8tG348OH42te+pmwzefJkvPXWWzhw4AB27NiB9957D6FQSLWN1vLlyxEMBpV/O3bssOcgaUDItAmqdhUCdhcmmXaeiIiIKLP0WyhUdnY2TjjhBNTU1GDRokUAeidv19TU4Mc//nHc7+bm5uLoo49GKBTCSy+9hO9+97sx2wwbNgzDhg1De3s7qqurcc899xjuLycnBzk5OSkdDw1cDBMxh+eJiIiI4um3hgUALF26FJdccglmz56Nk046Cb/5zW9w8OBBXHbZZQCAiy++GEcffTRWrFgBAHj33Xexc+dOHH/88di5cyduv/12RCIR3Hzzzco+q6urEY1GMWXKFDQ1NeGmm25CaWmpsk8iIiIiIrJfvzYszj//fLS1teHnP/85Wltbcfzxx+P1119XJnQ3Nzer5kYcPnwYP/3pT/H5559j+PDhWLBgAZ577jmMHDlS2SYYDGL58uX44osvMGrUKJx33nm466674Ha7+/rwiIiIiIiGjH59j0WmyoT3WBARERER9Tcr9eIBtSoUERERERFlJjYsiIiIiIgoZWxYEBERERFRytiwICIiIiKilPXrqlCZSsxn7+jo6OeUEBERERH1H1EfNrPeExsWOvbv3w8AGD9+fD+nhIiIiIio/+3fvx8FBQVxt+FyszoikQhaWlowYsQIOBwO3W06Ojowfvx47Nixg0vSDjG89kMXr/3QxWs/dPHaD0287kdEo1Hs378fXq9X9X45PRyx0OF0OjFu3DhT2+bn5w/5G26o4rUfunjthy5e+6GL135o4nXvlWikQuDkbSIiIiIiShkbFkRERERElDI2LJKUk5OD2267DTk5Of2dFOpjvPZDF6/90MVrP3Tx2g9NvO7J4eRtIiIiIiJKGUcsiIiIiIgoZWxYEBERERFRytiwICIiIiKilA3phsXGjRtx9tlnw+v1wuFw4JVXXjHcdvHixXA4HPjNb36j+vyuu+7Cqaeeiry8PIwcOVL3u83NzVi4cCHy8vIwevRo3HTTTejp6bHvQMiyvrr2Docj5t+qVavsOxCyLNVrv337dlxxxRWYNGkSPB4PJk+ejNtuuw3d3d2q7/773//GnDlzkJubi/Hjx+Oee+5J0xGRWX1x7bdv36773P/rX/9K45FRInbk+eeccw4mTJiA3NxcjB07FhdddBFaWlpU2/C5zyx9cd35zKsN6YbFwYMHMXPmTDz88MNxt1u9ejX+9a9/wev1xvyuu7sb3/nOd/DDH/5Q97vhcBgLFy5Ed3c3/vnPf+KZZ57BypUr8fOf/9yWY6Dk9MW1F55++mns2rVL+bdo0aJUkk4pSvXaNzY2IhKJ4PHHH0dDQwMeeOABPPbYY7jllluUbTo6OnDmmWfC5/Nh8+bNuPfee3H77bfj97//fVqOiczpi2svvPHGG6rn/oQTTrD1WMgaO/L8uXPn4sUXX8Qnn3yCl156CZ999hm+/e1vK7/nc595+uK6C3zm/39Rikaj0SiA6OrVq2M+/+KLL6JHH310dMuWLVGfzxd94IEHdL//9NNPRwsKCmI+X7duXdTpdEZbW1uVzx599NFofn5+tKury6bUUyrSde3j7ZsyQ6rXXrjnnnuikyZNUn5+5JFHooWFhapnfNmyZdEpU6bYlXRKUbqu/bZt26IAoh9++KG9CSbb2HXt//rXv0YdDke0u7s7Go3yuc906brufObVhvSIRSKRSAQXXXQRbrrpJkyfPj2pfWzatAllZWUYM2aM8llVVRU6OjrQ0NBgV1LJZnZce+Hqq6/GUUcdhZNOOglPPfUUolzhOaMlc+2DwSBGjRql/Lxp0yacdtppyM7OVj6rqqrCJ598gvb2dtvTTPaw49oL55xzDkaPHo1vfOMbePXVV+1OKtnM6rXft28fnn/+eZx66qlwu90A+NwPRHZcd4HPfC82LOL41a9+haysLFx77bVJ76O1tVXVqACg/Nza2ppS+ih97Lj2AHDnnXfixRdfxPr163HeeefhRz/6EX73u9/ZlEpKB6vXvqmpCb/73e/w3//938pnfO4HJjuu/fDhw3HffffhL3/5C9auXYtvfOMbWLRo0ZCuaAwEZq/9smXLMGzYMBQVFaG5uRl//etfld/xuR947LjufObVsvo7AZlq8+bNePDBB1FbWwuHw9HfyaE+ZOe1/9nPfqb8/6xZs3Dw4EHce++9KTdYKD2sXvudO3fiW9/6Fr7zne/gqquu6oMUUrrYde2POuooLF26VPn5xBNPREtLC+69916cc845aUk7pcbKtb/ppptwxRVXIBAI4I477sDFF1+MNWvWsJ4wANl13fnMq3HEwsDbb7+NPXv2YMKECcjKykJWVhYCgQBuuOEGTJw40fR+SkpKsHv3btVn4ueSkhI7k0w2seva6zn55JPxxRdfoKury57Ekq2sXPuWlhbMnTsXp556aszkTD73A49d117PySefjKampjSlnFJl5dofddRROPbYY3HGGWdg1apVWLdunbL6D5/7gcWu665nKD/zHLEwcNFFF2HevHmqz6qqqnDRRRfhsssuM72fiooK3HXXXdizZw9Gjx4NAFi/fj3y8/Mxbdo0W9NM9rDr2uupq6tDYWEhcnJyUtoPpYfZa79z507MnTsXJ5xwAp5++mk4neo+moqKCtx6660IhUJKHO769esxZcoUFBYWpv9AyDK7rr2euro6jB071vY0kz2SzfMjkQgAKB1FfO4HFruuu56h/MwP6YbFgQMHVC3Kbdu2oa6uDqNGjcKECRNQVFSk2t7tdqOkpARTpkxRPmtubsa+ffvQ3NyMcDiMuro6AIDf78fw4cNx5plnYtq0abjoootwzz33oLW1FT/96U9x9dVXs3LZj/ri2v/tb3/D7t27ccoppyA3Nxfr16/H3XffjRtvvLFPjpH0pXrtd+7cidNPPx0+nw+//vWv0dbWpmwreiW///3v44477sAVV1yBZcuWYcuWLXjwwQfxwAMP9MERkpG+uPbPPPMMsrOzMWvWLADAyy+/jKeeegpPPPFEug+P4kj12r/77rt4//338Y1vfAOFhYX47LPP8LOf/QyTJ09GRUUFAD73magvrjufeY3+XpaqP7355ptRADH/LrnkEt3t9ZYhu+SSS3T38eabbyrbbN++PTp//vyox+OJHnXUUdEbbrghGgqF0ndglFBfXPvXXnstevzxx0eHDx8eHTZsWHTmzJnRxx57LBoOh9N7cBRXqtf+6aef1v2+Njv96KOPot/4xjeiOTk50aOPPjr6y1/+Mo1HRWb0xbVfuXJldOrUqdG8vLxofn5+9KSTTor+5S9/SfORUSKpXvt///vf0blz50ZHjRoVzcnJiU6cODG6ePHi6BdffKH6Hp/7zNIX153PvJojGuXal0RERERElBpO3iYiIiIiopSxYUFERERERCljw4KIiIiIiFLGhgUREREREaWMDQsiIiIiIkoZGxZERERERJQyNiyIiIiIiChlbFgQEREREVHK2LAgIqK4NmzYAIfDga+++qpf/n5NTQ2mTp2KcDgMALj99ttx/PHHJ7Uvh8OBV155xfD327dvh8PhQF1dHQBg69atGDduHA4ePJjU3yMiGkrYsCAiIsXpp5+OJUuWqD479dRTsWvXLhQUFPRLmm6++Wb89Kc/hcvlMv2dSy+9FIsWLYr5fNeuXZg/f77p/UybNg2nnHIK7r//ftPfISIaqtiwICKiuLKzs1FSUgKHw9Hnf/udd97BZ599hvPOO8+W/ZWUlCAnJ8fSdy677DI8+uij6OnpsSUNRESDFRsWREQEoLeX/6233sKDDz4Ih8MBh8OB7du3x4RCrVy5EiNHjsSaNWswZcoU5OXl4dvf/jYOHTqEZ555BhMnTkRhYSGuvfZaJXwJALq6unDjjTfi6KOPxrBhw3DyySdjw4YNcdO0atUqnHHGGcjNzY353eOPP47x48cjLy8P3/3udxEMBgH0hko988wz+Otf/6och/g72lCo9957D7NmzUJubi5mz56NDz/8MObvnHHGGdi3bx/eeustayeUiGiIyervBBARUWZ48MEH8Z///AczZszAnXfeCQAoLi7G9u3bY7Y9dOgQfvvb32LVqlXYv38/zj33XPzXf/0XRo4ciXXr1uHzzz/Heeedh69//es4//zzAQA//vGPsXXrVqxatQperxerV6/Gt771LdTX1+OYY47RTdPbb7+N73//+zGfNzU14cUXX8Tf/vY3dHR04IorrsCPfvQjPP/887jxxhvx8ccfo6OjA08//TQAYNSoUTH7OHDgAM466yycccYZ+OMf/4ht27bhuuuui9kuOzsbxx9/PN5++21UVlaaPp9EREMNGxZERAQAKCgoQHZ2NvLy8lBSUhJ321AohEcffRSTJ08GAHz729/Gc889h927d2P48OGYNm0a5s6dizfffBPnn38+mpub8fTTT6O5uRlerxcAcOONN+L111/H008/jbvvvlv37wQCAWV72eHDh/Hss8/i6KOPBgD87ne/w8KFC3HfffehpKQEHo8HXV1dcY/jhRdeQCQSwZNPPonc3FxMnz4dX3zxBX74wx/GbOv1ehEIBOKeEyKioY4NCyIisiwvL09pVADAmDFjMHHiRAwfPlz12Z49ewAA9fX1CIfDOPbYY1X76erqQlFRkeHf6ezs1A2DmjBhgtKoAICKigpEIhF88sknCRtFwscff4zjjjtOtf+KigrdbT0eDw4dOmRqv0REQxUbFkREZJnb7Vb97HA4dD+LRCIAesOOXC4XNm/eHLO6k9wY0TrqqKPQ3t5uU6qTt2/fPlVDioiIYnHyNhERKbKzs1UTru0ya9YshMNh7NmzB36/X/Uv3gjDrFmzsHXr1pjPm5ub0dLSovz8r3/9C06nE1OmTDF9HFOnTsW///1vHD58WLUfPVu2bMGsWbPi7o+IaKhjw4KIiBQTJ07Eu+++i+3bt+PLL79URhxSdeyxx+LCCy/ExRdfjJdffhnbtm3De++9hxUrVmDt2rWG36uqqsI777wT83lubi4uueQSfPTRR3j77bdx7bXX4rvf/a7SSJk4cSL+/e9/45NPPsGXX36JUCgUs4/vf//7cDgcuOqqq7B161asW7cOv/71r2O22759O3bu3Il58+alcAaIiAY/NiyIiEhx4403wuVyYdq0aSguLkZzc7Nt+3766adx8cUX44YbbsCUKVOwaNEivP/++5gwYYLhdy688EI0NDTgk08+UX3u9/tx7rnnYsGCBTjzzDNx3HHH4ZFHHlF+f9VVV2HKlCmYPXs2iouL8Y9//CNm38OHD8ff/vY31NfXY9asWbj11lvxq1/9Kma7P/3pTzjzzDPh8/lSOHoiosHPEY1Go/2dCCIiIiM33XQTOjo68Pjjj6e0n66uLuTm5mL9+vWmRx+6u7txzDHH4IUXXsDXv/71lP4+EdFgxxELIiLKaLfeeit8Pl9KYVkdHR3405/+BKfTidLSUtPfa25uxi233MJGBRGRCRyxICKiQe/666/HCy+8gJtuugk33nhjfyeHiGhQYsOCiIiIiIhSxlAoIiIiIiJKGRsWRERERESUMjYsiIiIiIgoZWxYEBERERFRytiwICIiIiKilLFhQUREREREKWPDgoiIiIiIUsaGBRERERERpYwNCyIiIiIiStn/BwjhZtNevA21AAAAAElFTkSuQmCC", - "text/plain": [ - "
" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import lightkurve as lk\n", @@ -72,9 +61,26 @@ "# normalize\n", "flux_median = np.median(flux)\n", "flux /= flux_median\n", - "error /= flux_median\n", - "\n", - "# plot\n", + "error /= flux_median" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Let's plot the light curve." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": { + "tags": [ + "hide-input" + ] + }, + "outputs": [], + "source": [ "plt.figure(figsize=(8, 3))\n", "plt.plot(time, flux, \".\", c=\"0.7\", ms=2, label=\"PDC-SAP flux\")\n", "plt.xlabel(\"time (btjd)\")\n", @@ -87,7 +93,11 @@ "cell_type": "markdown", "metadata": {}, "source": [ - "For this tutorial we decided to focus on a single sector, but searching and combining different TESS sectors can be done following the approach presented in the [ground-based search tutorial](./ground_based.ipynb)." + "We see that TOI-540 display some stellar variability, likely due to the presence of starspots rotating with the star. A traditional approach consists in modeling and removing this signal before searching for transits, which could affect the transit signal themselves. Luckily, nuance is equipped to search for transits while modeling stellar variability.\n", + "\n", + "```{note}\n", + "For this tutorial we decided to focus on a single TESS sector, but searching and combining different TESS sectors can be done following the approach presented in the [ground-based search tutorial](./ground_based.ipynb).\n", + "```" ] }, { @@ -96,7 +106,7 @@ "source": [ "## GP kernel optimization\n", "\n", - "The Next step is to define and optimize the GP kernel that will help model the covariance of the data (mostly stellar variability). Here, we will use a mixture of two SHO kernels, implemented in `nuance.kernels` and representative of a wide range of stellar variability signals." + "The next step is to define and optimize the GP kernel that will help model the covariance of the data (mostly stellar variability). Here, we will use a mixture of two SHO kernels, implemented in `nuance.kernels` and representative of a wide range of stellar variability signals." ] }, {