{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Removing noise from K2 light curves using Self Flat Fielding (`SFFCorrector`)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "You can use `lightkurve` to remove the spacecraft motion noise from K2 data. Targets in K2 data move over multiple pixels during the exposure due to thruster firings. This can be corrected using the Self Flat Fielding method (SFF), which you can read more about [here](https://docs.lightkurve.org/tutorials/04-replicate-vanderburg-2014-k2sff).\n", "\n", "This tutorial demonstrates how you can apply the method on your light curves Using `lightkurve`.\n", "\n", "Let's start by downloading a K2 light curve of an exoplanet host star." ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:23:09.582635Z", "iopub.status.busy": "2023-11-03T14:23:09.580391Z", "iopub.status.idle": "2023-11-03T14:23:11.611065Z", "shell.execute_reply": "2023-11-03T14:23:11.610751Z" } }, "outputs": [], "source": [ "from lightkurve import search_lightcurve\n", "lc = search_lightcurve(\"EPIC 247887989\", author=\"K2\").download() # returns a LightCurve\n", "\n", "# Remove nans and outliers\n", "lc = lc.remove_nans().remove_outliers() \n", "\n", "# Remove long term trends\n", "lc = lc.flatten(window_length=401)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:23:11.613081Z", "iopub.status.busy": "2023-11-03T14:23:11.612966Z", "iopub.status.idle": "2023-11-03T14:23:11.859673Z", "shell.execute_reply": "2023-11-03T14:23:11.859336Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAvoAAAF4CAYAAAAlj12bAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAADJYklEQVR4nOydeXxVxfn/P/dCLiFAhMgSQBC0lU0bgyzGYBFqACWgXwxFFGlQaWtsLUUFmrSK1bSgVqNt08WFiKBSUle0shhEwVTQYmQLiyAiCkFCCCFkIbm/P/g9hzlz55x7zrn35i4879crL+XcOXNmn2eeeeYZl9fr9YJhGIZhGIZhmJjCHe4EMAzDMAzDMAwTfFjQZxiGYRiGYZgYhAV9hmEYhmEYholBWNBnGIZhGIZhmBiEBX2GYRiGYRiGiUFY0GcYhmEYhmGYGIQFfYZhGIZhGIaJQVjQZxiGYRiGYZgYpHW4ExDrNDc3o7KyEm3btoXL5Qp3chiGYRiGYZgowev14tSpU0hKSoLbbV8/z4J+iKmsrMSMGTPCnQyGYRiGYRgmSlm0aBE6d+5s+z0W9ENM27ZtAZypoISEhDCnhmEYhmEYhokWamtrMWPGDE2etAsL+iGGzHV++9vf+my5jB8/HuPHjw9HshiGYRiGYZgowan5Nwv6LcQTTzzBGn2GYRiGYRimxWCvOwzDMAzDMAwTg7CgzzAMwzAMwzAxSMSa7pw6dQqvvvoqdu3ahV27dqGmpga/+tWvcO2111p6v7GxEUuXLsXatWtRU1ODPn36YNq0aUhNTQ0orMiyZcuwZMkS9O7dG3/9618d5VOktrYWx44dQ3Nzc8BxMUxL4na70alTJzZPYxiGYZgIImIF/erqarzyyivo0qUL+vbtiy1btth6v6CgABs2bMDEiRPRo0cPvPfee3jooYeQn5+PQYMGOQ5LfPfdd1i+fDni4+Md55Fobm7G8uXL8d///jfguBgmnFx55ZWYPHmyI1+/DMMwDMMEl4gV9JOSkrB48WJ06tQJu3fvxuzZsy2/u2vXLnzwwQeYMWMGJk2aBAAYPXo0fvGLX6CoqAiPPfaYo7Aizz//PPr164fm5mZUV1cHkFNoQv6ECRNw0UUXoVWrVgHFxzAtTVNTE/bu3Yu33noLADBlypQwp4hhGIZhmIgV9OPi4tCpUydH727YsAFutxvjxo3Tnnk8HmRkZGDx4sU4cuQIunTpYjsssXXrVmzYsAFPPfUU/vGPfzhKI1FbW6sJ+aNHjw4oLoYJJ3369AEAvPXWW5gwYQKb8TAMwzBMmInJ/fW9e/eiZ8+ePoLGJZdcAgDYt2+fo7DAGc3lP/7xD4wZM0YTbALh2LFjAICLLroo4LgYJtxQO6Z2zTAMwzBM+IhYjX4gVFZWKncD6NnRo0cdhQWAd999F0eOHMEjjzxiK021tbWGz5ubm+FyudDU1ASXy8X2zUzUQmZnfKCcYRiGYcJPTAr6DQ0NiIuL83nu8Xi0352Era6uxtKlSzFlyhScd955ttI0Y8YMw98SExNx4MABNDQ0ICkpCeeff76tuBmGYRiGYRhGJiYFfY/Hg8bGRp/nJLSTEG837JIlS9C+fXtkZmbaTtOiRYuUNssHDx7EX/7yF/Tq1QsXXHCB4yuOGYZhGIZh/FFcXIyioiJkZ2cjKysr3MlhQkxMCvpJSUk+JjfAWbthUWNuNew333yDlStX4s4770RlZaUWrrGxEU1NTTh8+DASEhLQoUMHZZoSEhKUgn58fDzcbjfcbnfUe9t56aWX8PLLLyt/mz59OiZPngwAuOOOO1BRUQHgjP/1Ll26YPDgwbj11lu1nZInn3wSe/bs8bmf4PDhw1i+fDn+97//4dixY4iPj8cll1yCMWPGID093TBt//nPf/DRRx/hyy+/RH19PXr37o2srCxceeWVhu+88cYbePbZZzFkyBA8+OCD2vMnn3wSJSUlfvPZ1NSE119/HatXr8aRI0fQqVMnXHXVVZg6dSratm2rvVNRUYEXXngBW7ZsQV1dHXr27ImsrCxdfn7zm99g69atym/ef//9+OEPfwjA906ICy+8ED/5yU+QkpKie8fKNwFg48aNWLp0Kb766it07NgR1157LW6++eaob6sMwzDnKkVFRTh06BCKiopY0D8HiElBv2/fvvj8889RW1urE6537typ/W437NGjR9Hc3Ix//vOf+Oc//+nzzTvvvBMTJ07EzJkzQ5KnaMHj8SA/P9/nuey5KD09HTfeeCOampqwc+dOvPTSS/jyyy+xYMECwzMK5eXlmD9/PhITE5GVlYVevXrh1KlT+OSTT/D444+jR48euroV+de//oXBgwfj+uuvR3x8PDZs2ID8/HzMmjULP/rRj3zCHzt2DC+//DI6duzo89vNN9+M6667Tvfsww8/xJtvvokhQ4bovrls2TLceuut6NevH/bv348XX3wRlZWVuO+++wCcEcxpETFz5ky0b98ea9euxcKFCzF//nwMHjwYAHDXXXf5nPN488038dFHH+mE+GeeeQZr167Fbbfdhp49e2LNmjWYP38+HnvsMXzve9+z9c3y8nLk5+fjhz/8IX7yk5/gq6++wosvvoi6ujrccccdynJmGIZhIpvs7GxNo8/EPlEv6NfV1eHIkSNITEzUtMHp6el47bXX8O6772q+8RsbG7FmzRr069dPJ3RaDdu7d2/k5ub6fH/JkiU4deoUZs6cie7du4c6uxGP2+1G//79/Ybr2LGjFm7QoEFoaGjA0qVL8cUXX+D73/++T/iGhgYsXLgQ559/Ph577DHdomzYsGG47rrr0K5dO8PvFRQU6M5VpKam4vDhw3jttdeUgv6iRYswfPhwbedBpHv37j51/cILL6BXr166hca6deswcuRITcP/gx/8ANXV1fj3v/+NpqYmtGrVCl988QW+/vpr/OEPf8Bll10GAEhJScG2bduwfv16Teju3bu3Tzoef/xxpKamavk6evSotus0YcIEAMDgwYNxzz334JVXXsFvf/tbALD8zZdffhl9+/bFvffeq8Xl9XqxePFiTJo0ybH7W4ZhGCZ8ZGVlsSb/HCKiBf0VK1bg5MmTmmnNxo0btf/PzMxEu3btsHv3buTm5mLq1Km45ZZbAAD9+vVDeno6Fi9ejOPHj6N79+4oKSlBRUUF7rnnHt03rIY977zzkJaW5pPGN998EwCUvzHWIW3z4cOHlYL++vXr8d133+F3v/ud0gTKSJNPqA5PX3zxxVi1apXP823btuG///0v/v73vxtemCZy9OhRbN++HbfeeqvueVNTk8/io127dvB6vbowAHR5crvdSEhI0IWT2bFjBw4fPoxp06Zpz/bt24fm5makpqZqz1wuFy6//HK8/fbbaGxsRFxcnOVv7t27F9dee63uu4MHD8bzzz+PzZs3870PDMMwDBPhRLSg/9prr+k0qqWlpSgtLQUAXHPNNaYa3NmzZ2PJkiWarXKfPn3wwAMP4NJLLw0obDSwceNGnD59GldddVVYvk+CpIg/m+7Dhw8DOHNmQsXWrVvhdrt9bM0DYfv27bjgggt0z+iehB//+MeGaZFZt24dmpubNTt5YsyYMXj11VcxfPhwXHLJJThw4ADeeustjBs3TiuP/v37o3fv3njxxRfx85//XDOjOXjwIO6++27Tb8bHx+vOGNChctmLVFxcHBobG3H48GFccMEFlr+p8khF/z5w4IClsmEYhmEYJnxEtKD/3HPP+Q1z2WWX4a233vJ57vF4cPvtt+P222/3G4edsDJ//OMfbb/jBK/Xi/r6esTHx5uG27x5M3JycgAATz31lOkB1VBQV1eHG2+80ef5ggULMGjQIO3fXq8XTU1NOH36NHbt2oV//etfSE5OxsUXX6yM9+jRozjvvPPQpk2boKTz/fffx44dO3zMsd555x3U1dXhhhtusBzXunXr0L9/fyQnJ+ueT548GY2Njfjd736nacqvueYa3TmOVq1aIT8/H4888oj23OPx4P777zc0gWpqasL69esxbNgwXXvo0aMHAGDXrl3o1q2b9pzOm9TU1Nj6Zo8ePbBr1y7dt8vLywEAJ06csFo8DMMwDMOEiYgW9JkzeL1e3H333di0aRPy8vKUgjQhuuf0p0VvaGhASUkJ+vfvH5RbfoEzAuOCBQt8nsua83feeQfvvPOO9u/vf//7+MUvfmEqyAfL9ei+fftQWFiIa6+9VmdyVVVVhaVLl2L27NnKuxVUHDhwAHv37sXPfvYzn99WrFiBN998E3feeScuuugifPXVV1iyZAn+8Y9/4K677gIA1NfXY8GCBfB6vcjNzUVCQgI2bNiAxx9/HPPnz9ds6EU2b96M48ePY+TIkbrnF154IQYNGoSioiJ06dIFPXr0wJo1a3y89Vj95vXXX4+nn34ab775JkaNGqUdxnW73ewGlmEYhmGiABb0o4CGhgZ88skn8Hq9+Oijj0wF/csvvxzPPvssTp8+rfMAo+Jvf/sbXnzxRbRt2xarV6/2u1tgBbfbrbSxlxkxYgQmTZqE1q1bo3PnzoZuSYnzzz8fZWVlaGho0N1tYJeKigrMnz8fl1xyiY9pzNKlS9GnTx8MHDhQ0343NTWhubkZNTU1aNu2rc/iad26dWjVqhWuvvpq3fPq6mo8//zzmDFjhnYw9tJLL0VCQgL+9Kc/YeLEiejZsydWr16NXbt2YdGiRdo5gpSUFHz77bdYvHix8ozAunXr0KFDB+3QrMisWbOwcOFC3H///QCArl274uabb8ZLL72kmSJZ/eaPfvQjfPnll3j++efxzDPPoHXr1pg6dSrefPNNy2ZNDMMwDMOEDxb0W4jZs2f7uI0cP348fvCDH/h9t02bNsjLy8NHH31kyX3n5ZdfbilNpJUNh3b2vPPOs7QgIC677DKsXr0aZWVlGDp0qKNvHj9+HA888AA6duyI3NxctG6tb/5ff/01tm3bhqlTp/q8O3XqVMyfPx9XXHGF7vkHH3yAlJQUn8O+hw4dQmNjIy666CLdc/r3t99+i549e+Krr77C+eef7/P+RRddpPTVX19fj48//hjXXHONT/oBIDk5GU8++SQOHz6M+vp69OzZE2+88QaSkpLQtWtXALD8TbfbjZkzZ+KWW25BRUUFunTpgqamJrz44ovo16+fz7cZhmEYhoksWNBvIZ544gmltxirhxpvuOEGW3bjVrjrrrswcOBAXHLJJUHR5ocS8oy0ePFiDBo0yKcsv/zyS7Rr187HXz9x6tQpzJ8/H6dPn8Yf/vAHZV3MnDkTJ0+e1D175pln4PF48JOf/MTHvGnnzp349ttvcfPNN/vERen44osvdGcT9uzZAwCaDX3Xrl3x3Xff4fjx4zrBe8+ePZpgLrJx40acOnXKx2xHhuKvr6/HqlWrkJGRof1m95vt2rXTvBotWbIE3bp1C+qhaIZhGIZhQgML+ucwcXFxPu4TA6W5uVk7sCnSsWNHn8OqdvB4PJg7dy7mz5+P2bNn44YbbkCvXr1QW1uLzZs3Y+XKlfjTn/5kKOj/4Q9/wL59+3DPPfegoqJC582JDqDK2nfgjJAbHx+vtJVft24dPB6P0rVqp06dcOWVV2LJkiVoamrCxRdfjK+++govvfQSLr/8cvTq1QsAMHLkSCxfvhzz589HVlYW2rZtiw0bNuDzzz/H7Nmzld/s0qULBg4cqMznihUrkJCQgC5duuDw4cN444034PF4dD6TrX5z165d2Lp1K/r27YuGhgZ8/PHHWLt2LebPn8834zIMwzBMFMCCPhNUGhoaNPtwkYyMDJ87DOzSv39/FBQUoLi4GP/6179QVVWF+Ph4XHLJJbj//vtNfel/9tlnAIAnn3zS5zeV1yZ/iJ5v2rZtqwzz61//Gq+88gr+85//4OjRo+jUqRNGjhyp87ffpUsX5OfnY8mSJfj73/+O+vp69OjRA7Nnz8aoUaN08dXU1ODTTz/FDTfcYGhu1djYiJdffhnfffcdEhMTkZaWhmnTpul2bKx+s3Xr1vjoo4/wyiuvAAAuueQS/PGPf7R0IRrDMAzDMOHH5TW7lYcJmNraWkyZMgXLli0zNN154oknMHv2bE3LyzDRCrdnhmEYhgke/uRIf7j9B2EYhmEYhmEYJtpgQZ9hGIZhGIZhYhAW9BmGYRiGYRgmBmFBn2EYhmEYhmFiEPa600IYXZhFl1s1NTWFIVUME1yoHcttnWEYhmGYlocF/RbC6MKs2tpaAMDevXt9LmRimGhj7969AM7cI8AwDMMwTHhhQT/MJCQk4Morr9R8uV900UV8GRETdTQ1NWHv3r146623cOWVVzpyAcYwDMMwTHBhQT8CmDx5MgBnFzcxTCRx5ZVXau2ZYRiGYZjwwoJ+BOB2uzFlyhRMmDABx44dQ3Nzc7iTxDC2cLvd6NSpE2vyGYZhGCaCYEE/gkhISGBBiWEYhmEYhgkK7BqDYRiGYRiGYWIQFvQZhmEYhmEYJgZhQZ9hGIZhGIZhYhAW9BmGYRiGYRgmBuHDuC2E0c2448ePD1OKGIZhGIZhmFiGBf0WwuhmXIZhGIZhGIYJBWy6wzAMwzAMwzAxCAv6DMMwDMMwTExSXFyMzMxMFBcXhzspYYEFfYZhGIZhmCjjXBdgrVJUVIRDhw6hqKgo3EkJCyzoMwzDMAzDRBnnugBrlezsbCQnJyM7OzvcSQkLfBiXYRiGYRgmysjOzkZRUdE5K8BaJSsrC1lZWeFORthgQZ9hGIZhGCbKONcFWMYabLrDMAzDMAzDMDEIa/RbCL4wi2EYhmEYhmlJWNBvIfjCLIZhGIZhGKYlYdMdhmEYhmEYholBWNBnGIZhGIZhmBgkYk13Tp06hVdffRW7du3Crl27UFNTg1/96le49tprLb3f2NiIpUuXYu3ataipqUGfPn0wbdo0pKamOg67f/9+vPzyy9izZw+OHTuGNm3aoHfv3pg0aRKGDRsWlHwzDMMwDMMwTDCIWI1+dXU1XnnlFRw4cAB9+/a1/X5BQQFef/11jBw5EjNnzoTb7cZDDz2Ebdu2OQ575MgRnDp1Cj/60Y/w05/+FDfffDMA4OGHH8a7777rLKMMwzAMwzAMEwIiVtBPSkrC4sWL8fzzz2PGjBm23t21axc++OADTJ8+HbfffjvGjRuH/Px8dO3a1ecGOTthhwwZgoceeghTp07F2LFjMXHiROTn56Nv3754/fXXA8swwzAxC19VzzAMw4SDiBX04+Li0KlTJ0fvbtiwAW63G+PGjdOeeTweZGRkoLy8HEeOHHEUVkWrVq3QuXNnnDx50lFaGYaJffiqeoZhGCYcRKygHwh79+5Fz549fdxZXnLJJQCAffv2OQpL1NXV4fjx4/j222/x+uuv49NPP0VKSkqws8EwTIyQnZ2N5ORkvqqeYRiGaVEi9jBuIFRWVip3A+jZ0aNHHYUlnnvuOc0m3+12Iy0tDT//+c9N01RbW2sp7XFxcYiLi7MUlmGY6ICvqmcYhmHCQUwK+g0NDUph2ePxaL87CUtMnDgR6enpqKysxIcffojm5mY0NjaapsnqOYOpU6filltusRSWYRiGYRiGYYyISUHf4/EoBW8S2kmItxuW6NWrF3r16gUAGD16NH73u9/h4Ycfxp/+9Ce4XC5lmhYtWmTpZlzW5jMMwzAMwzDBICZt9JOSknDs2DGf5/Ts/PPPdxTWiPT0dOzevRsHDx40DJOQkGDpjwV9piVgLzAMwzAME/vEpKDft29fHDx40McufufOndrvTsIaQdp/q3b4DBNu2AsMwzAMw8Q+US/o19XV4cCBAzh+/Lj2LD09Hc3NzbpLrBobG7FmzRr069cPXbp0cRS2qqrK5/unT59GSUkJPB6PZs7DMJEOe4FhGIZhmNgnom30V6xYgZMnT2qebzZu3Kj9f2ZmJtq1a4fdu3cjNzdXd4i1X79+SE9Px+LFi3H8+HF0794dJSUlqKiowD333KP7hp2wf/3rX1FbW4tLL70USUlJqKqqwvvvv4+vv/4ad9xxB9q2bdsCpcIwgcNeYBiGYRgm9oloQf+1115DRUWF9u/S0lKUlpYCAK655hq0a9fO8N3Zs2djyZIlWLt2LWpqatCnTx888MADuPTSSx2Hvfrqq7F69Wq88847OHHiBNq2bYvvfe97yM7OxvDhw4OUa4ZhGIZhGIYJHJfX6/WGOxGxTG1tLaZMmYJly5ZZ8rrDMAzDMAzDMEDgcmTU2+gzDMMwDMMwDONLRJvuxBKzZ8+G261fV40fPx7jx48PU4oYhmEYhmGYWIYF/RbiiSeeYNMdhmEYhmEYpsVwZLpTU1NjOezu3budfIJhGIZhGIZhmABwJOj/4he/QFlZmWkYr9eLZcuWYc6cOY4SxjAMwzAMwzCMcxwJ+seOHcMDDzyA5557DqdPn/b5vaKiAr/5zW+wdOlSUxeYDMMwDMMwDMOEBkeC/h/+8Ad07twZb775JmbPno39+/drv61duxb33HMPtm/fjtTUVPz5z38OWmIZhmEYhmEYhrGGo8O4gwYNwp///GcUFhbigw8+wOzZs3Hrrbfiiy++wIcffoi4uDjMnDkTEyZMCHZ6GYZhGIZhGIaxgGOvOwkJCbjvvvswZMgQFBYW4oUXXgAAXHzxxZg9ezZ69eoVtEQyDMMwDMMwDGOPgNxrNjU14csvv0R9fT3ogt3Tp0+jubk5KIljGIZhGIZhGMYZjm/GPXjwIO677z689tprSEpKwm9/+1v88Ic/xP79+zF79my88cYbwUwnwzAMwzAMwzA2cKTR/89//oPnn38e9fX1SE9Px91334327dtj2LBhGDp0KP72t7/h+eefx6effopZs2YhKSkp2OmOOvhmXIZhGIZhGKYlcSTo/+1vf0Pbtm0xa9YsjB49WvfbyJEjMWjQIDzxxBP47LPP8Mtf/hJLly4NSmKjGb4Zl2EYhmEYhmlJHJnu9O/fH0899ZSPkE907twZ+fn5+MlPfoJTp04FlECGYRiGYRiGYezjSKO/YMECHzMUGZfLhZtuugmDBw92lDCGYRiGYRiGYZzjSKPvT8gX6du3r5NPMAzDMAzDMAwTAI697jAMwzAMwzAME7k4Mt254YYbbIVnV5sMwzAMwzAM07I4EvQ7d+6sfO71enHs2DE0NTUBALp06eI8ZQzDMAzDMAzDOMaRoP/cc88Z/tbc3IzPP/8c//znP9GnTx/cf//9jhPHMAzDMAzDMIwzHAn6Zrjdblx++eV44IEHcM8996C4uBiTJ08O9meiDr4wi2EYhmEYhmlJgi7oE8nJyejXrx/WrFnDgj74wiyGYRiGYRimZQmp1502bdrgu+++C+UnGIZhGIZhGIZREDJBv6KiAlu3bkXHjh1D9QmGYRiGYRiGYQxwZLpTUlJi+NupU6fwzTffYO3atTh16hSuv/56x4ljGIZhGIZhGMYZjgT9goICuFwuw9+9Xi8AYNSoUbj11ludpYxhGIZhGIZhGMc4EvRvvvlmQ0G/devW6NSpEy699FIkJycHlDiGYRiGYRiGYZzhSNC/5ZZbgp0OhmEYhmEYhmGCSEi97jAMwzAMwzAMEx5Y0GcYhmEYhmGYGMSS6U5eXl5AH8nPzw/o/ViAb8ZlGIZhGIZhWhJLgv6WLVscf8DMO8+5BN+MyzAMwzAMw7QklgT9Z599NtTpYBiGYRiGYRgmiFgS9Lt27RrqdDAMwzAMwzAME0QcuddsCU6dOoVXX30Vu3btwq5du1BTU4Nf/epXuPbaay2939jYiKVLl2Lt2rWoqalBnz59MG3aNKSmpjoOu2vXLpSUlODzzz9HRUUFOnTogH79+uG2225Dz549g5JvhmEYhmEYhgkGlrzuPPXUU1i9enWo06Kjuroar7zyCg4cOIC+ffvafr+goACvv/46Ro4ciZkzZ8LtduOhhx7Ctm3bHIf997//jY8++ggpKSmYOXMmxo0bh23btmHWrFnYv3+/47wyDMMwDMMwTLCxJOi/99572L59u/K3G264AU8//XRQEwUASUlJWLx4MZ5//nnMmDHD1ru7du3CBx98gOnTp+P222/HuHHjkJ+fj65du6KoqMhx2BtvvBHPPfccfvazn2Hs2LGYMmUKFi5ciKamJhQXFweYY4ZhGIZhGIYJHgH70fd6vfB6vcFIi464uDh06tTJ0bsbNmyA2+3GuHHjtGcejwcZGRkoLy/HkSNHHIUdMGAA4uLidN/q0aMHevfujQMHDjhKK8MwDMMwDMOEgoi10Q+EvXv3omfPnj7uLC+55BIAwL59+9ClSxfbYVV4vV5UVVWhd+/epmmqra21lPa4uDifxQTDMAzDMAzD2CUmBf3KykrlbgA9O3r0qKOwKt5//30cPXoUt956q2k4q+ZHU6dOxS233GIpLMMwDMMwDMMYEZOCfkNDg1Ir7vF4tN+dhJU5cOAA/v73v6N///4YPXq0aZoWLVpk6cIs1uYzDMMwDMMwwSAmBX2Px4PGxkaf5yS0kxBvN6zIsWPH8Pvf/x4JCQmYN28eWrVqZZqmhIQEvhmXYRiGYRiGaTEsC/obNmzAli1bfJ67XC7D3+j3Z555xnkKHZCUlKQ0uTl27BgA4Pzzz3cUljh58iTmz5+PkydPYsGCBcowDMMwDMMwDBNOLAv6dXV1qKurs/2by+VylrIA6Nu3Lz7//HPU1tbqtOg7d+7UfncSFjij6X/44Ydx8OBBPPLII34P4TIMwzAMwzBMOLAk6D/77LOhTodj6urqcOTIESQmJuK8884DAKSnp+O1117Du+++i0mTJgE4c/vtmjVr0K9fP50XHTthm5qa8Oijj6K8vBy//e1v0b9//xbMKcMwDMMwDMNYx5Kg37Vr11CnQ8mKFStw8uRJzbRm48aN2v9nZmaiXbt22L17N3Jzc3Xeavr164f09HQsXrwYx48fR/fu3VFSUoKKigrcc889um/YCfv888/j448/xrBhw3DixAmsXbtW9/uoUaNCVRQMwzAMwzAMY4uIPoz72muvoaKiQvt3aWkpSktLAQDXXHMN2rVrZ/ju7NmzsWTJEqxduxY1NTXo06cPHnjgAVx66aWOw+7duxfAmQXHxo0bfeJhQZ9hGIZhGIaJFFzeUFxry2jU1tZiypQpWLZsGXvdYRiGYRiGYSwTqBzpDkGaGIZhGIZhGIYJMxFtuhNLzJ49G263fl01fvx4jB8/PkwpYhiGYRiGYWIZFvRbiCeeeIJNdxiGYRiGYZgWg013GIZhGIZhGCYGYUGfYRiGYRiGYWIQFvQZhmEYhmEYJgaxZKNfUlIS0EdGjx4d0PsMwzAMwzAMw9jDkqBfUFAAl8tlO3Kv1wuXy8WCPsMwDMMwDMO0MJYE/ZtvvtlH0P/222/x/vvvo02bNkhNTUXXrl0BABUVFfjss89QV1eHa665Bt27dw9+qhmGYRiGYRiGMcWSoH/LLbfo/v3NN99g9uzZuOaaazBz5kx06NBB93tNTQ2eeeYZbNy4EY8//njwUsswDMMwDMMwjCUcHcZ94YUX0L59e/zqV7/yEfIBoH379rjnnnvQrl07vPDCCwEnkmEYhmEYhmEYezi6MGvLli1ITU1Fq1atDMO0atUK/fr1w+bNmx0nLpbgm3EZJvopLi5GUVERsrOzkZWVFe7kMAzDMIwpjgT9hoYGHDt2zG+4qqoqNDY2OvlEzME34zJM9FNUVIRDhw6hsLCQBX6GYRgm4nFkutOnTx9s27YNn332mWGYsrIybN26FX369HGYNIZhmMgiOzsbycnJAIBDhw6hqKgovAliGIZhGBMcafQnT56M/Px8/P73v8fIkSNx9dVXo0uXLgCAI0eOYP369Xj//fcBgLVdDMPEDFlZWcjKytKZ8DAMwzBMpOJI0B8+fDjuuusuPPvss3jvvfd8LtTyer2Ii4vDT3/6UwwfPjwoCWUYhokUSOBnGIZhmEjGkaAPANdddx2uuOIKrF69Gtu3b0dlZSUAICkpCQMHDsS1116Lbt26BS2hDMMwDMMwDMNYx7GgDwBdu3bFrbfeGqy0MEzUwt5YGIZhGIaJNBwdxmUYRg95Y+HDmQzDMAzDRAoBafS/+uorrFy5Ert370Z1dTWGDx+OGTNmAAB27NiB3bt3Y9SoUcpLtRgmlsjOzubDmQzDMAzDRBSOBf3XX38dL7zwApqamgAALpcL1dXVujDPPfcc4uLicN111wWWyhiAL8yKbfhwJsMwDMMwkYYjQX/Tpk14/vnn0a1bN9x+++0YOHAgbrvtNl2YAQMGIDExER9//DEL+uALsxiGYRiGYZiWxZGN/uuvv474+Hg8/PDDSEtLw3nnnacMd9FFF+HgwYMBJZBhGIZhGIaJfoqLi5GZmYni4uJwJ+WcwZGg/8UXX6Bfv37aDZFGJCYm4tixY44SxjAMwzAMw8QO7Lii5XEk6Dc2NqJt27Z+w1VVVaFVq1ZOPsEwDMMwDMPEENnZ2UhOTmbHFS2II0G/W7du2Ldvn2mYxsZGfPnll+jRo4ejhDEMYw/eEmUYhmEimaysLKxYsYKdV7QgjgT94cOHo6KiAq+//rphmFdffRXV1dW46qqrnKaNYRgb8JYowzAMwzAijgT9m266CUlJSVi0aBEWLlyIDz74AMAZU53S0lI88cQTeOmll9CtWzd2H8kwLUS4t0R5R4FhGIYxgueI8ODyer1eJy8ePHgQCxYswP79++FyueD1euFyuQAAXq8XvXr1Ql5e3jlvulNbW4spU6Zg2bJl7F6TiWkyMzNx6NAhJCcnY8WKFeFODsMwMU5xcbF2USGbgkQ+PEc4I1A50vGFWT179sRTTz2FjRs3YvPmzaioqEBzczM6d+6Myy+/HFdddRUfxGWYcwi+HZhhnMNCq31Ec0Uus8iH54jw4Fijz1iDVmI9e/bkm3EZhmEYJazttA8vjphzgbBo9Ldu3Qq3242BAweahjt48CCOHTuGSy+91MlnYgq+GZdhGIYxgrWd9snKymIBn2H84Ogwbm5uLn7zm9/gkUceQV1dnWG45cuXIy8vz3HiGIZhIgGzQ2R8wIwJBux2kGGYUOBI0AcAl8uFjRs34v7770dFRUUw08QwTBhhwdUXM9el7NaUYRiGiVQcC/ojR47EhAkTsH//ftx7773YunVrMNOFU6dOYenSpXjwwQcxdepUTJgwAWvWrLH8fmNjI4qKivCTn/wEN910E+69915s3rw5oLCBpolhogEWXH0xc10abremDMMwLUleXh6GDRvGFhtRgmNB3+12Y+bMmbj77rtRU1OD3/3ud3j33XeDlrDq6mq88sorOHDgAPr27Wv7/YKCArz++usYOXIkZs6cCbfbjYceegjbtm1zHDbQNDFMNMCCqy9mZhVscsEwzLnE6tWr0dzcjNWrV4c7KYwFHAv6xNixY/HII4+gXbt2+Nvf/oZ//OMfaG5uDjhhSUlJWLx4MZ5//nnMmDHD1ru7du3CBx98gOnTp+P222/HuHHjkJ+fj65du/poKe2EDSRNDBMtsODKMAzDGJGRkQG3242MjIyQfofNSINDwII+AAwaNAhPPPEEevfujXfeeQcPPvggampqAoozLi4OnTp1cvTuhg0b4Ha7MW7cOO2Zx+NBRkYGysvLceTIEUdhA0kTwzAMwzBMtJOfn4+NGzciPz/fcRxWhHjRjDQShP5ISIMTgiLoA0DXrl3x2GOPYfjw4SgrK8Ps2bNx6NChYEVvi71796Jnz54+7iwvueQSAMC+ffschQ2E2tpaS3+NjY1B+Z4VorXRMgzDnGuw5ycmlrByFkw0I42Es2ORkAYnBE3QB4D4+Hjk5ubixz/+MQ4dOoQdO3YEM3rLVFZWKjXv9Ozo0aOOwgbCjBkzMGXKFL9/y5cvD8r3rBCtjZYJLSw0+IfLiAkmdrWbdn5jmEjEylkw0Yw0Es6ORUIanODowqzRo0ebXpY1bdo0XHjhhVi8eDHCcfFuQ0MD4uLifJ57PB7tdydhA2HRokWWLsxSpSVU8AUtjAq+Vt4/52oZ8U2kzvBXboWFhaiurkZhYaFhuZqN1zyWM9GG3cvOIuFytEhIgxMcafRnzZrl9xDG1VdfjWeeeQbPPvuso4QFgsfjUZrAkNBOQrzdsIGQkJBg6a8lBX0+dMmoiFatRUtyrpZRrGmOg7EzYyWOwsJCHDp0CIWFhY6/YbZQ4LGcYaxzru3IBtV0J1JISkrCsWPHfJ7Ts/PPP99RWIY5F2ChwT/nahnF2gInGAuXYMSRk5OD5ORk5OTk+P1GSwop55pAxJwbxJrCwh+WBP2tW7di69atmpab/m31r6Xp27cvDh48iNraWt3znTt3ar87CcswDHMuE2sLnGAsXKzE4U+QNyvX4uJi1NbWIj4+HrW1tdruQCiFFBLwCwoKAtqJYBirtOSi0m6/j/YFryUb/dzcXLhcLhQWFqJnz57av63yxhtvOE6gP+rq6nDkyBEkJibivPPOAwCkp6fjtddew7vvvotJkyYBOHP77Zo1a9CvXz906dJFe99OWIaJRtiumgkl0dy+gmFzS3GQMKAqh0C+Q/b7LpcLdXV1cLlcSExMRHZ2dsjKnjSeduZ5hgmEgoIC1NXVoaCgIOTjiN3+GO1nsiwJ+qNGjYLL5dIOk9K/Q82KFStw8uRJzfPNxo0btf/PzMxEu3btsHv3buTm5mLq1Km45ZZbAAD9+vVDeno6Fi9ejOPHj6N79+4oKSlBRUUF7rnnHt037IS1mqZoIJonZ8Ye0T5IMZENt68zhLoc2rRpg/r6es3BRVZWFkaPHo3q6moUFBRoWvecnJyAv0+He1NSUlBWVhYVplqxMKfFQh6sIue1vr4eALT/2nk31ET7YXeXNxxucSxyxx13oKKiQvnbs88+i27dumHLli0+gj5w5jDtkiVL8P7776OmpgZ9+vTBtGnTMHjwYJ+47IS1kiaR2tpaTJkyBcuWLbPkdaelyMzMxKFDh5CcnIwVK1aEOzlMCDmXJg+m5TkX25eYZwA+QnEwy0H8Fmn3ExMTUVJSogn6LpdLWwCEekyP1PqOhTktFvJgFTmveXl5WL16NTIyMvxexHUulRMQBDnSy4SUkydPejMzM70nT54Md1J05ObmeocOHerNzc0Nd1LOGZYvX+4dP368d/ny5eFOCsPEFC3dt8aPH++94oorvOPHj9f+f+jQoSH/vpxP+vdtt93mHTJkiDc9PT3kaRDzHgjBrrNYGF8DnZejqQycpnX58uXeUaNGeUeNGhUV+QwGgcqRMel1hzmL0SGSsrIyNDc3o6ysLEwpO/c41076M0xL0dJ9SzzMl52dDbfbjebm5pB/Xz60S/+urKyE1+tFx44dQ65lD5bnpWDXWTQfFKd5urS0NKB5OVhl2hKHT53WV1FREaqrq5GQkBCVdR0OLNnol5SUBPSR0aNHB/R+LDB79my43fp11fjx4zF+/PiQftfIbjTabc6iES5zhgkNLd23VIf5wtm3g51/Ms9RmSMF69IgozTbMQ0yCxuJJkZGaaJ5Oj4+Hm63GykpKY7iT0lJQUVFheP35fQE+7xJMOok2G3djslQtGLJRn/ixImODt96vV64XK6Qet2JdMJto++vY0XiYMhEJ9yWGMYeweozwe57ZANNOxUtaQttx/7aLGwk2nEbpYnqr7a2FtXV1Y7THKw8h2osj8Q6GTZsGJqbm+F2u7Fx48ZwJ0dJoHKkJdOdm2++2dHf1KlTcfPNN9tOFBM8VNtjtC2Xl5eHRx99lM1JoohI9ufLpkmRQzDbSSS3uWjHSZ9R1UewTTZSUlKQnJyMjIyMFr8gzY5pkFnYQE2MxHIOVh8wShPN02lpaQFp9INlVhUqMygxfZEyrmRkZMDtdiMjIyOs6QgpQTwvwCiIxMO44uGxljpExgSHYB2ECwXRdBAsGERyfoPZTiK5zdklWE4IgnUg0EobksOo6iNYbTGUcUcbqgPXRn0gWO0hlvqaP0aNGuW94oorvKNGjTIMI7a9c7Uder18GJdxAK2qSVszZ84cNrWIEoKlsQkF0XwYzgmkRS0sLIwIzZRIMNtJS7S5ltLurV69Gs3NzVi9enVA6SA3l9XV1QFp0a30GbmdkbZd1IoCUMZjt1yprlNSUrT3InmnLpTtJiUlRdOu++sDdEDUaXuQd1Kc9LVI0ZBbpaGhQfdfFWLbi+R2GOmwoH+OUltbi9LSUranjjLCLUxH22RiBad5oskfQMRNQMFsJ8Fuc6E0PfGHv216u+lwuVyWhTKzdmb2m9zOysrKtPrwt9i0mx+q69LSUi1ef0Kunf4T7PHDTv7sftuOZ7rs7GwkJiZqNxbbhfIh1q3TOCJpHDLD4/EAOHNBllGdyN6tIlXJFekEJOhXVFTg7bffxjPPPIOnn34aTz31lM/f008/Hay0MkFC1D48+uijMSW0MaEl2iYTKzjNEwlFOTk5YZmAonXRpSrvlprEU1NT0bVrV6Smpip/t5oOqvO5c+ciKyvLUl2YtTOz32T77aSkJO1b/habwShXfws9O/2HwtK8k5eXh2HDhiEvL08LY6dd28mf3X4uavT9vZuVlYWSkhKUlJQ4EtKDUU/RJAhT3dJFb0VFRcp6F9teuJVc0YxjQf/ll1/GT3/6U/zzn//EihUr8N5772l/1ODp30xkQdoHl8sVct/P0SqMRAPhKFsrk0m01XmgE2S4JqBgmg61ZJ2pynvz5s2oqKjA5s2bg/otOV9WBDYrdSmHsyJEqvJtxWRD9rFeXl6uc31otth02jZzcnKQmJiofd8Mu4dnxTsHVKZUdgRyMX+ikwlVW7bbz6m8aefbyq6G0bft5MMpTuII11hNysYOHTpo5SovAp1gNT/RNkcFiiNB/8MPP8TLL7+Mzp074+6778bll18OAHjooYdw11134dJLL4XX68UNN9yARx55JJjpZYIAaR/GjBmjO+EfisYfjM5rhEobdC5A9VRYWNji2nWjyURsO3Y1Z+EedIMxyYZz0QUEbjoUyE6N3byrytuf7bxT5HwFS+uZl5eHoUOHYsSIESguLtZpf41Q5duKyQaFAYDk5GR07doVAJCUlGQatxOoLjdv3oyamhrdrq9RPdv5dlZWFubMmaPVgcqUSq4jq+2Lymn16tXKtmzmgS7Qfqv6drjHNSvQHFJYWBjUeI08FskL25ycHK1O5EWgHI8VrI5jdvMdDXVphiNB/+2330br1q3xhz/8AWPGjEGnTp0AnNkWve6665Cfn4/bb78db731ls8lUecqs2fPRk5Oju7v7bffDmuaZBvEUJhlqDpvsLAjHASqcYkk5InfqdASzMFLbDt2hSl/7S4aBtlwmDQF03QoEAE4GHkPlYs7OV/+3A1b1QauXLkSXq8XdXV1KCoqcnzTuLhAMGrnlAcSiioqKgAA5eXltr5lBRKAVq1ahebmZgDQxm6jerZ79kCsg/z8fGzcuFF3UZGT3RLA18lEoOY81K9ycnL8Coa0S+7xeDQ7/WCOCdEwBoqIeVf9f2lpqc878iJQjscKTheJdvITjTiSwr/88ksMGDBA0yzQZVpe4e6tG2+8ET179sS//vWvICQz+nniiSdQWFio+wv1rbiEvwmEOkUobPzkzhvMAcuOcEAddeXKlSHXYIQaeeJ3qsUL5uAlth272kV/7S4aBtlw2scGQ1MZiEY4GHlXCXzBwI5XGxobCgoKTOMU22F8fDyys7MtafRViOYhRu1czkMo/X6TB5RWrVohOTkZY8eO9TkMmZSUpNtJNeufBQUFlsrUDKvti8opPz/fcls2i5u0zEVFRaaeYShsQkIC6urqkJCQoL0brDHByhjoZA4SFzPBhPKekpKC2tpabfHjbxeS6hCAbS9EtJtMcxBgXG52TNPE/IRjfA8KTnxyTpo0yfvoo49q//7rX//qnThxovfEiRO6cI899ph36tSpjvx+xgqR4Ec/knzzhiMty5cv96anp3uvuOIK3V+gPo9FIqWM7fgatuuXONThQx3PuURLtken9dNS9Sp/x+jfNC4MGTLEdrqt+FlX5ZV8iaenpwfFD3ugqHyby2kX719R/S4yZMgQS2UabILVtsT6MWtD8rNgt20r8Tnp86Hug0ZpMvouPadyN7rjR/W+lbZrJW2RSFj86CclJeH48ePav88//3wAwP79+3XhvvvuO237jwkfodrOsoJsR+9U8xUIRUVFqKur83lu1eexlfKirdva2lrTcIGWver94uJijB49GqNHj7Zlt29Xi2tXsx4sTbzddIpl5PQcRyD1FMr+ZTXultRAOa1n8T2xDcttO9CylNMn/5vaF2kaXS6X6fdU7dHprhRpVD0eD6qrqzVtsBmhbF8qDa+c9oyMDLhcLng8HhQXF5v2TzoHNmbMGEfpcZrXYO8CejwepTmRVbemTvNB7wH6exJU8Tnp86HaKSCMTNNkzT2ZzdH8BcDU5FeV7lOnTun+CxjPHXl5eTh8+DDi4+N1d0YEmt9IxZGgf+GFF+LgwYPavy+77DJ4vV689NJLmkD14YcfYvv27ejdu3dwUso4xsjmccGCBSE/yCrb0Vu1ZQ1WZysuLkZVVZXh74HacYrQ4TWzcIFOQKr3i4rOukttaGgImYBndyIJ13anWEZG5zj8tS+jerLiDjCUpkZGAoRsXx7I4sguduuZhPqqqiqdPbPqwqFgHBYUzQjMzAHI9t3qeSKV4GJU3kZlJLvOVClA/LWvYAomVhYx+fn56Natm3Y+wYxATbKc9qVgjT1UL2lpacr4gbMmKGJa5e8HYzHs77kT8zsrirdAxjNxvjeau8RDzNXV1QCAhIQEH3t9EVX9nj59GgDQ2Njotz+sXr0aXq8XDQ0NKCsrM1UEiETrIsCRoD9s2DAcPXpUE9YGDhyIyy67DFu2bMHUqVNxyy234PHHH4fL5cKUKVOCmmAmcMROHWwvFwR1iP79++vsSa0OwE4GF1UnLDLQ5gNnNHeqQVH2EiDaGJqlt7m5GW632zRcoBOQLLQUF5/xpU3nZEKJ3YmkJd1OinUmlrGRPbO/9iXGIcZNC4eVK1caCl6hXOAYCRBGnkasEshkbueMQHFxMR599FFUV1f72DMHcuGQlfTRJVCrVq3S2fESpKkmu3sz8vLysGDBAr+LEFEjSwsa1Y7FunXrNFt9GX/tK1SeUwhV/aoExFAIQU77kmhfH8gOq5FiSnUQXkyrXGZO82EkiBuNT3axongLZDwT31XFQ89onCYOHTpkOn+ofhswYAAAIC4uzu9YJs4LNPbQbrxZfkOpxAkpTux9GhoavAcOHPBWV1frbIj+8pe/eKdNm+a98cYbvXfffbd37dq1juyJYolw2+ib2ZNeccUV3ttuuy3oNnrLly/X7Did2r85sR1U2ejl5uZ6hwwZoqWH7P6GDh3qzc3N9RuPVTs+K+kNpj2knC6Km84iiGXQkvizhw7V96jO7NqCWkmX2B5yc3O1diSXvZVvhMp2Nzc3NyD7YH/v2I3TqN/Q8yFDhliyR8/NzTXtq3agenRiSy8jjifU1/zZ7qvKRCwP+q/cd2677TZL45VsR26HYNSv6pkcr5P6tNs2/Y0JqnzQWB/IOGGWLqfxWJl/ArE1b4nzNVbHweXLl2t9avjw4bbS4/X6r0uzdAVzng8FYbHRj4uLwwUXXIAOHTpozxISEnD33XfjxRdfxGuvvYa//OUvuOaaa4K1HmEcolqB0op13rx5qKysDPoK1ap224xgaYNLS0vh9Xp1HqG8Xi/mzJmDsrKyoNk62/HuEcgFR7IfYkoXfZ+uFQ8XlEfyvR1qDQjFD6hdjRp932n7ys/Px7x58xzvSqmeBWK/S2YCoqcRJxpes/IgLbzKP7gVj16q3Za5c+daukXUrttKs7IkTxtmuwZW26uo/c/JyTEsI7Gfmmk0yZbd+/9vCRXTUl5ebloGpFUG4EizL6edMDvfYpYX8ZlcnqIpndG5DBl/dSL/Lo8JKSkpGD16NEaMGOHzLUozYM0LjJ0+KvZDp+OglfmHNNJVVVWmZanqG07HQTv5MRsH5bmQ+ud9991nKz2AfW90YrqCOc9HIuzkPsZRbf2JjdXJFqC/cBTnnDlzTC9WCuQbKsxcZomC/pgxY0wFrrS0NG3AsLoFbAUakE+cOGF5kDSyzzW6YMeOuzSzMnYqfGb/fzOi5uZmFBYWhtSMBTjbvtPS0pTl4fT7Yv7lMqX+s3nzZgwbNkxzISjXp6pvqezD7UyaYrqCuYgyq2954S4KMPT/Cxcu1B1mo8kzKytLl05VfzL7tp36MxJYCboo0GyBYfV7+fn52LRpE9avX6/lkRxPpKSkKPupSkigZ6mpqfB4PHC5XNpYLZo1WEkT2Sj7cwUpQmUm1i/VB/nSV5l3muXFzL5fNJkoEs5lmC1O/NWJ/Lss8JWVlWmmYnQGRD7k6u8uikD7mtNxSB5rpk+f7nMmR3TtKV5wFuw8OMmPkekrvd/Q0KBbDNGBdAA+i0B/c5I/IVx+X5SNolWAtwoL+jGOP42Y2MCtDgT+wpl1GqN3g3GgkQY88UBfTk6OzvYvMTERpaWluoOActpkIdpKWqwIxpQ+r9cLl8vl10OP6tv+Blg7A5ZZvswWQv7yKO70EZs3bw7I64SsPab/Jz/kVtq3UbwqZOFUFQcJQnV1dT5tqbi4GIWFhaitrdXFp1qg2RECRE1YVVWVTjAk0tLS4HK50NDQYGkR509ANlu4k1Cp0kQbtVkjDWwgt2erBFazsEaeq2hxanfCz87O1saZsrIy2wqUov9/lsjr9Wpt2apfeCq/pqYmADDc1TM6w0QLlP79++t2hFq1ahWwv36576SmpqJr165ITU3VlAJ24zD7XV5kAmcEOpfLhdatWxteZuX0MLUZooLA7s6AvJuycuVKNDc3Y/v27T5nckiYjo+PB3DmMDndWyDWub+Dt3aUO1bnGVF4Fz0GUR2JbVUsY3ERaDSuWIXyJXukc3rJXTTiWNDfv38/CgoKcOedd+Kmm27CDTfcYPjHhO9mXDsDlNWwTnYBCHmwMeqETt1wynnIyjpzYRdtCTY0NPgcBPSXfyvlIgsrRsIEaTc6dOhg6KFHZeog5icQzYM4gZjlS/Wbv4FWtSMi3rSpMv2Q82u04Hv00Ud1bYSeA85uBxbzoqorK3Xepk0bnzjFdIsTlVl8duqU4gGgCYYrV67UmSWQuZrsFYWE2YULF/oI2mYCsiykiPVLE7XL5dJpVMXDbaqDifLhNyNXelYn98LCQh+BVUS1EyKaDagECzvQOKM6jGklD1QmiYmJPu7+/EFj5YABA0x381TpEOtbvmk3ISEh6JeYyQvouXPn+qSZ6mr69Om2XeOq8lhWVgav14uEhARNW2xXcHfixUpecBilT4XZre/Jycmak4uUlBQUFBRo3moIckAhfs9MsC0uLtbGBaPdFbtzvZE2X9wRVC2GjA7nZ2dnIz4+HocPHzZtE0ZzCJUdxUf9JikpKSo96djB5RVtGiyyZcsWzJ8/H42NjXC5XGjfvr22mlTx3HPPBZTIaKa2thZTpkzBsmXLtEEmUqDBKCUlRdNC2RUgMzMzcejQISQnJ2vCgBh3dnY2Nm/ejNWrVyMjI0NzZUXhR48ejerqarRu3RrNzc3IyMhAfn6+9jwxMRElJSVBy+/ChQs1M564uDjce++9QdmuEzWKos2n2+3GnDlzAED3O2ktVGVulHejyUNOh1mYYcOGaULdxo0bDeOgwZ4GYLPnBLUFynNWVpaWl/j4eHTs2FGnTaM2ILYhKjf6TSzXxMREJCQk6DTDTtqsXE5yeuzEQeXR0NCAuro6LZ8pKSmaBxVVWQVKcXExFixY4PM8OTkZtbW1qK6uhsvlwty5c7VvU10A0NWRlXYFACNGjNDyuH79ei0dqneNxgWj3/Py8rQxggRLKt+GhgZ4PB7DcpT7NQB88sknht+jOqdyomdmbTsQrJavKq1G7dFJ+1Wlo7i4GAUFBaivr8eYMWOQn5+vhQPO9MWBAwdi8eLFuriM6stfPlXvGeWfEMcq+RvivwEo2wuFEevbaK4C1OOKWd5Uv6nq0N/4aVZO8r/F+A8fPqztErdp00bXR63kTS5zsX+L+TMqPyOM2jGNQ/L4ZIWhQ4dqed20aZOl7/obo0jJ4G9eDieBypGONPpFRUVobGzEj3/8Y7z88stYunQpnnvuOcM/Jnz4s70VtwGdbJ0baUZETYKonaDwpLWirf+mpibdIa1QUFRUpBMGGhsb/WpWrGgxaAD3eDy6q75FLSVpTQFofv3taubNtEFGOyMyRu4m5e+otJsq0ygRlWaWNDazZs1Snguh94wOLIqaUvGQlbxd7+TCk9raWhQWFmp28ykpKbjyyisxZMgQTJ8+XVm+eXl5mu0oAM3mm6irq8OhQ4dQWlpqag9uVzsmk5WVhbFjx8LlcmnmD3FxccjOztbKXJ5ERdttccuc6qSgoABDhgzBiBEjlOmqr6/X/Vd8t0jYzQD87/rJO3YqbSO1Q9G2WoXcrwHfMzriLsLmzZsBQNuZSEpKwsKFC9HQ0IC0tDSfvASKXW2w0Ziq2pUgocTpOZii/28y1K1bN03wpvSS4Ld9+3af98hsbdWqVbq4/Gmr7bhzHDhwoO5iLtU3xH9Te6FdW9qxAczt8CkOcddQngvN8qb6TWxvFE9hYaGWPqM2RnWcmpqq202R7yMQ61y8lGzWrFnaeCsLuKoxk2z9k5KStDScPn1a11+pTOze0WLULnNycuByueD1em0fHKddVHk31ey7/saojIwMra7sXDgZTTgS9L/88kv069cP06ZNizgtNaPHbIASG7rZLXRmGE1iYmcThUsKT5p9j8eD5ORkzQcupcHOoVKryGZAVvxlW5m8ZIHkT3/6EwDotvJJCAbg96IZo7zLA5hq4gdgGAbQ28gaQZMULVrEScHsTgHZfIGeybap4sCbl5dnOBmJ8VqxUVZN/iJyedFihuzmy8rKtAONonAj2rCvXr1ad4BwxIgRGDJkiOFdDUbIaTQy9TIzcSorK8PcuXM1Qb+pqUkrQ9XEJgr3dXV1OmGGBD76TdU+RYGCzIBGjBjhYwoE6OtM5QVIFvhUQoHYZ8zs7kXBUCxfscyI6upqrFy5EocOHcLKlSuRnZ2N8vJyzdTJzn0E/hZrTheghHyuxUi4t2MvrRJI5XKndNIOG5Ur1fno0aO1Nte6dWvTuOQ4U1JS/N4gTvlZvHixz8VcZkoCEla7du3q44VIFPhVcxXNf4D6RlazvKl+UylFaKFNygCriwaxDEUlA+VFXASYmYyJ4zj1WTKp3LFjh/adpqYm5WJHvh3YH0btUjzH1dDQ4HO42KyfiAsZs+/K459Zuaampmp1BZydPwNVxkQSjgT9Dh06oGvXrsFOCxMCzAYo6hBlZWWWPTs4QXU7IqWLbhzcu3cvgLO2vvIgYdTprEy2NDiSKYXb7ca8efM0jxlmWDkrQIIxodopkM8KmJWz2cRNWgfaRaCJjMpTdi1mJgQblV1Wlt47Cb1DQq58tsFK2o3SQYKVmctRf3Ust3GjOpPLS2UDSoJLXFycbpIQ7b/pPQA6Ad/tdmPs2LGWFqhymlUTkVGZyZpHGotbtWql7TiohG85TaIwI5aV0QJY7Mfi4pa06UaXJ8lX08u2u3l5eXj00Uc17xeE2Gfat29vWJbU5iZOnOhTp1RmhYWFOHHihM+7hYWFut2t/v372z4YbbQosLoANXpPXnQ4Ee7F8U/l8Um1EBeVBp988olmtiMujqnO27Zt6xOXKk0UZ1lZmU4Attu3CVoE0e4McPaMQUVFhU6DD6hdZ4rpJmVUWlqaTlkhLhTtXN6kSjsttOPj4w0XV3IdGSknzJQYRt+nOli5cqVWf23atEFycrJOQ96tWzcA0DynGd0OrELeKTCqV6obj8ejW3D4c0XqdFGrakeqhbM4f1qZK6MFR4L+kCFDsHPnTu2UPxO5mJ32FzWVq1at0jyEBKNR+5vUZM0+CUxt2rRRmmOo4hPTb9QZxcHxxIkTSExMVHoPMYK8uqxatcpUA1VSUqJpvuLi4pCSkuJzuInC0SFVlV9nM8wGequTjfhvM3/rYlnSO3YXg2IcovAtCnp0qIzcrKnaiz+/8HLeZW0xpUN2O5iQkKCzlc3KysJ///tfJCcn6xZromZ57969qKmpQVpaGnJycrSzSa1btzYVRv2lWTURqeqOXJdSX21ubtaEssbGRqxatUoTxGTPTllZWboFqSgQU1klJyf7LIBF7T21V1ooiRpd0exO7K+0S0L/pXZcU1OjvWfmxtHMXExMI924Ky5Eqd01NDT4mPcA0J1ZAIDKykrLXprMlCji70lJSRg2bBiSkpJsOTsQ+5vRYshfGsUxw8glrzyOqvJFfTY+Ph6JiYkYM2aM1pb8CXXywk6MX1y8qrS5cj+RF0G0O1NUVKTdbUCmPvSukcmO+C1xzDDbjbIzN8ppJ/fPHo9HW6yJ9UQLISNFDXkPUi3EVfOjaoyhdu5yuZCYmKiZVM6aNUv7TTwfAdjzTmO0SJXLTq4bcaHhr19ZQY5DNT8afUe16JLbabThSNC/7bbb4HK58PTTT+PkyZPBThMTRKjRqmzPRE2l1+vVJlOrmiezb5qZeKiESBKYSOthdUVu5N9bTEd8fLxmE2imjTZDdB9oxMSJE5GcnIx7771X8/RATJ8+XbP9FrWhdrx8yAO9ysTH30RpRSMiT/z0jmjuo9LcyN8WfayLbjAp/wkJCaisrERzc7NmwqXSZBG0zetvoCUhlIRcUUPpdrtxwQUXmHqYUE0SpPGrq6tDc3MzVq5cic2bN2P9+vX45JNP0LlzZ51ZiN2+o6qXzZs3o6KiQtNailveRkqWNm3aaFpt2bNTcXGxz2JHFAKMTCpU7ZUWrZ07d9bC0Q6B3P9lk5/a2lptwUICmtm5ESu7auJYIHqtISGlvr4eiYmJSqcRtNAAoLNVluP3J0jJfYLKdfv27WhubkZ5ebmlvkfCttjfzBZDRmmUhXMjgYY8FhmZR4mLqI4dO6KkpAT5+fnaAsyfuZPY32XzPBI+xbNMZnGJiyDRrIscSnTo0MHH9MzKLqOR0Ef9hVzVmrmhVZWt2Jfo+9XV1bp7J8R8yd8Xn9Oc0rFjR78KHSPat2+vLdREs2vaPZMXAjk5OT5KGrMxmMJ27drVp8+a1e3IkSOV3necIPYhs0UxLZ6Kiop08hH9f2lpqWE7jTYcCfrnnXce/vSnP2HPnj2444478Lvf/Q5PPvkknnrqKZ+/p59+OthpZmwgCjnydmRSUhLcbre2vQmcOSiq2j5U3WBo1OnJ3VdDQ4OPZlBedFCnJts7umVSXiioOj8NKmT7r8o7TU7kxs2uNjotLc2SuQ19j4RHcYLNycnRbL63b9+uCVX+JmASGOi/5DaRJh95t4a0flYFTXmhoJr4RaFJddOjOMkbDeS04KEDr2LdGl14JcZFW8cALGlVZC0wTYLAGWF0x44dSu2u+L7c1khgFDXYK1eu1NIhasvkcvOHUT9SCXdUZ2PGjNEtkIEzGvqLLrpIt+Mgb90bnSUQy0wsX1FYpD9xMUCCcWJiouaqTtasiyY/dCCRTAZIoG3fvj1KS0uV5UVlb3RQX0yjx+PRuXIlQcPr9aKmpgYjR47U3avhcrnQv39/7d+yi0nCiiBl1CcIf/7oZUFS7ANWDtHLCyJx/FMdChfnBtE8Su7HRoso1c6Dk7Ij4VPW9qugvpmfn68rE9oZpl1b1fuyb3rxW6o+n5eXp/WX+vp6wwWRqv/SWPn4448rXYR6vV4sWLBAm1NV35cFVrOy8bdzn5mZqXPDSeUlL4jIbEm8uVrU6PtbiFHYiooKNDc3o7S0VLmrS1B8RrtNTqA4zeYJCkOKGQCaokN14ZxYNoHsNIQLR4J+bW0tFi5ciK+//hq1tbUoKyvD2rVr8d577yn/mPBBg4PK9oyuVq+oqNDC19XVabbswFnXdSqTEaNOr/LOIYY3G4wB694JSBCiLUX5ZlxRSwnY83JDA/XKlSuRlpZmepsmIQqU8gRLZj0DBw7UtKGzZs1SHmaXBQZxMCKtgnyQksI3NzdbvozLaFscODvxl5aW+tTxiRMntEWiaNMsT0Rk2uJyuTRBnm6pJCHQaFtYnBRoN4BMP0StipGQrJrEiVatWmmLN5UtvZmZhlxf4kE/mgjmzJnjM5HKi2XZ9lY1MZkJd6mpqT5b7gkJCVqfXrlyJQoKCrS2T2UqXhykOuwta61IWGxoaNAOwFVXV2sX8pBgXF1drS2iSDtu1gbr6upQVVWl3a5r5sPenzZNTCOdGSBhTGxX1M6oXOPj4+H1elFZWYmxY8eaCtJGgpjYTlSCLz2bN28eUlNTfdqVKHyKArXcn1TnnGSMDjgbLTrFuUG1MCalkGhmR2MRLfZpoWaG6nCkXHdkxmVHmyvml9qI2a6tylOQGeICW1RWyKafZsLv6dOntW+S4kjEbEeX5qCCggLdDlFhYaHPjbEjRozA0KFDlbd0Uzy0aKmtrdXGgqqqKr/zhDgW+1uIyX0AOHs+gsZxUb4wuldHTpMdkyl5vDAaz0XFjNgHaHdZHh8D3WkIJ4786P/lL3/BqlWrcOGFFyo1SzKXXXZZQImMZsj/ac+ePXWaJAAYP348xo8f3+JpEretysrKkJKSglWrVvn44nW5XNr2PwAfv7dGW2RGvpJJY0U+a0XhS/R9W1VVpdM8UljxezS4in7qi4qKtHfJBz3FK/pgt9JRRZ/jZj57VWVrxVfy9OnTNS2/nDaxfkh4p3xu3rwZK1eu1N4rKSnxqU/ZP7hcR1bqjQTV+Ph4NDQ0ICMjA6mpqbrbR1X1KMctl798Z4NRWkT/+/X19VrbJCGO0qjy7WwUp+wH3iic6FN67NixSE1N1cIB0Mq6tLQUJ06cgNfr1X0/Ly8Pq1atQqtWrbRzAOIiSpwAaWGqisesPMU+IbY3sX0Qcl+QvyG3X/F7cp8V7Xc/+eQTrc30798flZWV+O6777QFGeVVzg/FJ0LaePneBrncVf1XbP/y3QVUPqJvdbEuiLFjx9q+GMrufReq8hfvtZgzZ46uH8t9VqxnVXnI4wCVm9GYqypDszSTlhY4U191dXW6tmwWv1HbM5sT/KVP1VbN2gjdO9G6dWt07txZGzuM7pSg+1/69++Pr7/+2rDsze4USEpKQnl5OTwej1Ze9F0AugW3XHe0AKZxTyxrALo2Q89cLhe6deumKwNxLiNEL0PiPKEqE6M7AYzKWkQsm9LSUp/+Is8PRr76rYxdRu2FytFfP7Wap3AQFj/6H3/8MTp37ozHHnsMEyZMwNChQ3HZZZcZ/jHAE088oWmv6C8cQj6g3/6k/5Jpy9y5c7VwtAYk7afsl9tohWvmwtHj8cDlcvncRidqCkQ3gCoPGqLWiQYjWWsBnOnI5LO+urralm2lqHnxer1+NQmi2YuRZxpRsyD7phbTJtaPqCUmDThBNusAdPV5wQUXADhjViFrm2iHRmWfLmvIRHOX0tJSna06afLpoCFtTcvfEw9DqrZoVXbO4sREQj7ZeFNZ0E6DyrezeDZA1J6LAqhZOYhxrV69WpcnsW4SEhJ02mPgzMRGXi1Onz6tueHMztZ7+RHLt7q6Gh06dFBqysRvU3t2uVw+5hnU3lJTUw1t0c20cVTu8i6AWOcpKSlavPHx8ZoNfEZGBiorK5Gdna2bhMR0iqgOLXfs2FETwuWdqsLCQhQUFODw4cM6LytiGqlO5LsLyCRJvAlb1OYRpHG1ojH0h5GGV1X+Zq6HxfepnlXnqKh9A/B5v6jI/NZjQC3wyOYWSUlJunGZ7gypra3Vdvf8eSaTTUfJZMWKSYSqTKkuScNN+ZfzAEDr4y6XCwkJCdrOsnhWRP4O7aJMnDgRNTU1PmVPu3ni+SM5bZWVlZgzZ45mmpqUlIQTJ04gPj4e8+bNU3o2o+/LZnrZ/9/UkVybUrrFch8zZowmtFPeKZ6BAwdqO1m080vjEX2fxl1x11jVbuXdBqN+I84pqjNl4jgomxgbtR8zxL5A+VKZ4lAd2T2/Fq04EvTr6+vRr18/Uy0+EzlY2fYSG7k4oJPZihXTFcJooisqKtK21svLy7XBcvTo0ZowlJWVhZycHJ05gnhg0MymkiBTEfoeYTbZyWlfvXq1tp1Pz628B6gHKjHMo48+qts2BPQLGtHMY/PmzTqXmqLASK7JxMk+Ly9PW0SUl5frtmjJht9oE0/cRpXNXWpra5GZmYl169ahubkZX3/9NVasWKEzFREnBTIXoEmQ8qoSBmTzBVH7NGDAAE24z8/P12meCKN2QBeykEBPgr7H40FR0dlLlsRDvqQ5bd26tWa/bXSwPDs7WzP/IAFUNAsQ7fnJXIv6EZWvOPmpTEPEb4v9Z+XKlUqXibTQ9Hg8mDdvnm5iFfuMPCZQuasOiNNkXVpaCo/Ho2nfabKXhRNq016v18cki75DbVg8pyKPG6IQQPk2OoxqNMZNnz4dhw4d0i4To3KYM2eObv5q06aNqSJANr0CzpqmnThxQmeDbSSUqMYs0SSH8qByg0n9ns5IiGFkJwTi9+k9M49QcrmL/6a637Fjh7YbJpv70Big8soiL0JEoVZWbKjapVGZyooVq2alALTFqNfr9Vlgq+pOXixZFU7FeqE80p0N8vk1itPfeQfRyQMtpKncExMTUVZW5uN4gw719+rVCxs3btQWHaIdvmrxSwsJs7m2vr7e0uFpozmbvt3Q0KAz8RQVP/7cm4qI9U7/b2SKQ9gxC4pWHAn6F110EY4dOxbstDBBRGUD7E9YJURttplLLdXkB/hqb0Q/zjRZkUYIgM9ALWp5SGAUL1YCfA8cyQJGaWmpJogBZzSQVtxqZmefdWG4atUqvwNvXl4ehg4diu+++07bhiWhRfRIk5eXh8OHDwM4M4C2b98eycnJaN++vc8OgKi9I3eJ5KmhsLAQOTk5mptOebIXBU3aLqXLgMhGFTi7GBJR2VASp0+fxqFDh3zOXYg2zTQppKSkaLa8dB7D4/FowoB8KE6+OZk0VgC0BYV8joC8ztTV1fm0P/FsgOxWMTExEWlpaaitrdWEeQC6yaG6uhqnT59Gt27dUFlZabhDk5WVpcVPAii5iYuPj8d9992nOzMit5thw4Zh8+bNPho4sTxF+2W5bMxcJlL6/Pk1J62k2DflRY24KyPePArARzihvmt0gJ0ET6on8ZyKLGRR+sV2qrKhz8vLw4IFC3T5obKkRW9jY6OPgEEekz755BNN+DHq56JASZpOVf0bmZj4cz8pOilQ1SktFDt27Gir3mlHw8wWXF6ci4sI+o3atXhpkhXhVG5nRhpiObycVrPzRNTXRSWC/A2aH7xer85sS15giwoOeVdDdvpANyurFulGGJ27kXfZjbT8olcoWkj7W3iI46uReQr1W1JOJCcn+9SRqEwCzoylAwYMUCpw7AjoogKD5A0aH1VnDsyg8klKStKNNXJ+5e8bHd6V56poxZGgf/PNN6O8vByffvppsNOjcerUKSxduhQPPvggpk6digkTJmDNmjWW3yc/2D/5yU9w00034d5771Vu+4YybDgRB0yzgVUFdTZ/GnBRIBU7ojgoi2HKyspQUlKi2SsCZw/QykKBvI1NgmNBQYEub6LANGfOHJ0glJWVhY4dOwKA0iWZCA1MYl16vV6sWrXKdJBYvXq1ZqZB31QN0BTO5XLpJgeVj2caoBITE3X+hckFquiXXx5EyZvJvHnzfGyORYGXFkMqSENJAzVw9pIX0X92cXEx8vPzNc2xOLGI3yEvMLTVTwsOCkdeT/r374+srCzdba8y1JbpUhfA90Ab1bvX69UWlBRfQkKCdiC4ubkZXq9Xc70oakD9eQGh9kL3ANDELd7cKAtZosAnTr5GvrpJwBA1a3PnztUtXuV0Wb1RWhYOqG+uX7/eZ/eOtLqkGSPBjtqwKJyQIEGLUVFLS4t1qrPCwkKlG1d5V4PaCWks5clY1vKLgiXlUbwATa7DvLw8zVOWFURTD1lwUwmqKmFXRBRa/d0aq/Jc4u9CI1mQl+OWx1pxEUECPfUjsV3Jwingq4CRTfeono0Wt9T//DkToDzl5ORofV02nZHbJLVbMmEZO3asX+GvSNjVIA8yZgsysY/LziEA9eWRqnfFPIoLKdErFJnA0hxvtDsotlGjhRSVkWh6J+dNVCaR8uPrr79WKofsKBflMRc4uzipr6+3JbvQDb87duywtMil7xsd9l+5cqW2Yx3NOBL0W7dujfHjx+Phhx/Gk08+iffffx9btmzB1q1blX9OqK6uxiuvvIIDBw6gb9++tt8vKCjA66+/jpEjR2LmzJlwu9146KGHsG3bthYLG07EAdPKQkSc9MhcgA65ylp70tJXVVVpGgCj0/uUjtatW+Pw4cM+K2PZpEH1Pk3WwJmtQlH4EgUmEoREQcfqIocGJvHWQMC//3zRNR+Zt5BGjAbolJQU7WzCgAEDfPKv0njQIU4SHMeOHav9prI5pEGeLkCh+Og5TThko0nxiBOLqPWiRQpt11O8ov9so3KhiYUmVNLk03tt2rTRCUeVlZUAzmjvaXEh1yOg15aKnqKMNMfJyWduXk5ISNDZuYoTqOypQ2VioxJKqL1UVlbqTC/EXafRo0fju+++AwB89913Oq2zmUcd0Y+zvLMFwHDxKh7GVAlRFEbUtBldJqQqS1mwB3wFO3GClzXV1KZolweAoZco8V3ahVKdsykuLtb61tixYzWhlwRL4Ix2srGx0UdjJ/Z5M68/wNl+RDcfU3nJ55HE8YaUEHRRFqC+oVUsX+ojCxcuVF6mp/JS5e9CIzPbf1U65LZAyhp/d5DIdT969GhtUS/fk2GmuVeNL7IQLI6b/sZ4+hYJ5osXL/YRtsX4ZZM5fxpz8X3ZdbQVQVNVfiLkYSsr6+x5GTIZJBNYo4UyoF9cmNWxmRlScXGxdrM03dVRVVXlc8jXzPzMCJUMQPkcM2aMLk/+NOykGKP/Wpn/5XMi1HZHjBjhN+3RgiOvOxMnTtROggMw1L4Rb7zxhu2ENTY2oqamBp06dcLu3bsxe/Zs/OpXv8K1117r991du3bh3nvvxYwZMzBp0iQAZ4SaX/ziFzjvvPPw2GOPhTwsEehp6UAQvYcAvh4wiOLiYkNvKuLBSDqJT9vY9Ew8vS96xyEvGOT2ToyDBCJxe9voWUpKCt577z2cPn0acXFxuPfee3UdX+XhR86fHLcMHaIE4NO25UPIYrwLFy7UhfV6vT4eZkRPAsCZyUL2FiSmTeVZpajorEchl8uFMWPG+HjmMINszwHgggsu8PEEQcKs6ImJysPtdmPjxo2G5WmUXnFSp/9XeSSi8P486NBkRIeBy8vLlfUuviN6rlDVo1m7k/MnekeS601sCxRO9uxCkNcfQm7D1Odat26N5uZmnRcksa7kuqd0Up9UpZfKmDxqpaWladq4tLQ0W23KzBuH6DFK7A/0jBYvq1ev1rVDcQFBZQ5AJ1DMmzfPJz/0rqquRE9EclrF/ksaa6O8q9qKmTeQoUOHart4mzZtsjQOiWOxmF6xXFUeeazES33Pbj1biV8OJwqOqr5H6RG9IZnlx4nXFVV+AbX3JjF+QO1FyKwc6H3RS5m8+FflH9B7iBLDGnl1EtOiag+qsigoKEB9fT0GDBigHZxXKQmM4hHHFtnTj+ihR0wzHa4Xx0jVN6zM4QR5qVJ5F1LlQ5ZFxPI2QuWlyIlXrmASFq87o0aNwqhRozRNL/3b6M8JcXFx6NSpk6N3N2zYALfbjXHjxmnPPB4PMjIyUF5ejiNHjoQ8bLghrQQtwszMcIqKzh42MrO3bG5uRmFhobayly+bEQ98kgaPDvgRdNBX1jyotBn0bPXq1ZppDJlO0coegOa33OgqdiNNiYjKjMXtdhsK+cAZoVUU8klrDECnwSTBivwQi9uEstkGHR4EoPOsQvbxtJiQb6P0Zwcsapfopk7RZIW+sWrVKtTW1iI1NVXLS9euXQHo7crF+pO1JhSX6GvbyF5YHJhJa5qSkqKNLaKGjOwvm5ubsXfvXp1Wzt+ZFKOdGZVNrvy+rNGTtWcknFFbqK2t9blkSjyYKx52B3w1slQOTU1NurML/g6i08Frt9uNxMRETastXnAGnF2Q0q2m4nkQsS0SRm1LZUYiao9pTKFzK3SYl9qcbBKkuoyJtNwyYn7kw5S0OKCdDerbZmMgCVNWNNbizoCZxrBVq1a6/xrVm0z79u19zkrQt0krDZxtr2bxippm0sobafaN6tlfulX22NnZZ8+TdOjQweddGgvEw99G4wugb2v+xjoqL9Fk1OggMADt/BSZwhnVqTxWiPMNpQ+AzsxHLBOj9FFajMr5xIkThnUims0Z1UtBQYFmA799+3afHTd/bQg4285JoKd+Rjv/qvfksjaah63c+kyQpt/o8LycD1GG8LdrpyIxMVFpBhttOBL0f/3rX2PWrFmW/1qavXv3omfPnj4TxCWXXAIA2LdvX8jDhhsaSMirANnnGZnXUCeWXWLKrrno8BkJwWR2AZw5GCle2kQeOujP5XJp15Mb2SOKgys969+/v+6in+zsbN3lJ9SZyY5f7shmk7EK0Ve7UZmJkMaKtkYBaBfMiHGWlZX5bBOKV6wD0LndXLVqFUaPHq0JuF6vV7PBJ9tweTEgC/7iZESLDDFNostBUQAkARuAZiYjDsh0CPnKK6/0ccmoEozF30QbXDHdNEivW7dOG5SBs4KcaJ8qHwoWF0ziN6zebCymQ04jpcvIzIUWykRTU5N2YdqsWbPQsWNHbcwgwVdEbp/0vTFjxmj9CNALDCphhw5enz59GgkJCboFBS0CGhoatG1tl8uFjIwM7RCu0eau0QRtZjIiCwYkaNbV1ekOF5N5FaVftZCi+iDBkeqIFgK0UKW6osUkLYYB34uOxMUZxWvWz4uLz7rqFW15ZcFCNC2gOrejgaNxW77N1mgxbSS0yKYkVAZmgqwVhYhRmuUzCCozSjltdAA8MTFR20Eke2i5bdu5nZXKS+z3omJIXpzS+SnyhONv4SQqsmQ7ftpxFc8kWEmfWC4qU0o7dSLXvThWDhw40MeBg7ywB/ybyBiZ3NIYSQ4PxPxRnquqqnRmaWZmjHKZkBmSv8PzBLV1Gues3HBPeZg3b56Py99oxZGg/6tf/Uq7fCISqaysVO4G0LOjR4+GPKwMafn8/TU2NlrMpTlmk6gZqkG0trZWu/FSPJBHWlZRCKcdgQ4dOqCurg4ejwcdO3aEx+PR+RyXv6MaXMkt2N69e+H1etG5c2dtYCFhpVWrVj6H6Iw0MWaaNNFrDw06NJAbHaKjCVw0XaOJOiEhQVsEif6KZbMWGoRJKKMbdIGzh29FARc4I/DJh6CMBGhxMqJFhtlhTuLUqVPa/2dkZGi20JRWumDt9OnTOq2cuFVKLtNSUlJ02jqjWzgJcWISD5aJE8GAAQN0CxlqA7RgIr/X5L/Zn7Alp0P0m02otP/iu2S/LXqioMm0oaFB6ebQbGu/rKxM60eixpMOtRr1Z5fLhezsbN0hTXLLR30yOTlZW5zm5ORo5aYSeuW2ZVRmImZ9jnx0UxyyQK7qvyUlJZrgSGVFwop49kFMF41FqkOqtDij/PrT9BUVFenMD42EM3ExbPVgtIhRmcrl6U95IbY7+YCmUd0YxWmkQZftsQH1fSCAvr+IZ1BIWJQXi2YKAlogUHu0ouGnejl06JDP4lQWNM3io3Yj736LShRyAuDvALQsKMsLJgB+7xhQKXTIU01tbS2Sk8/64Z83bx4WL16s1b2sZFKVF7lNNhpvjM5OyDeg02/yLg5gfkBZtYgUv2N1l0l1x4Yq/WLcAHTOAqIZR4L+wYMHte3ISKShoQFxcXE+z0lQERt2qMLKzJgxA1OmTPH7t3z5cou5NEfscGaTKKDXhspXY5PgSkI6TRgUL/kFTkhIQEJCgnZwiCY5ALoDeGSXbuSbXEQ+eS8OmrSiJ6ERgHYgz4rdqaq8xAFfND2gPIidvbi4WCtXMmkCoCs/mgBFf8XigFlUVKSVa1paGjIzMzFx4kR88sknmDdvns4VKSGbfYiXZgFnvZnImgxRSCOBRfTlTGkhyFQqPj4epaWlWLBggaatqqur03kDonDk758WGKLLNKpLMg2iupcH7OLiYrRp08anLouLi1FaWqptpVZWVuoWMuRuk3zki36vjXY7zAZ3s0uGjBaptI1eWVmpLcLIzM3j8SgP5hlNoPQcONtvZTekcp+gPkfmZvIFaKTRIm0b3c0gllebNm10np3IFI9caqoWPlb6G13iRoimayRw23FTCJx1eUmQiYNcHzRWqUxuOnTo4HPpmQqx/MSdAbn+ZA9S/nZhjNi8ebNfUzwrZVVXV6ctgu18X8RfG125ciVSUlI0LzPyglB+XyVgUvmKO8Ci5lfcqREP9tM4KtcvzVsURnQQIC8SysrKtHs6zPIrmvhQeHHxJLqVledJI2WRWC80X9TW1uLQoTMe5mQTQaPD76JChxQlTU1NpuY9ouJGRnabbLTYNyor1fhEz1XuUI2geAD41LEKeTfDn4CuSr8YB42LRgu2aMGRoN+jRw9tAotEPB6PUjNOA4ssOIUirMyiRYuwbNkyv3+TJ0+2mEv/kJ2+OPCYTQwnTpzAunXrNM2f2MHj4+N12khZawb4anNo4iSTENqyp0GY/t+o88gn71W2pWlpaZpWrk2bNkr3e1bKSbzcafXq1dqgO2fOHJ3Jhmj3qELUnJJAJpa5rJmib1DeCgoKMGzYMLz55pua+UVdXR0SExN12kHS1jY1NWkTCbXD2tpare5kLzli2fXv39/HJpq00qQF8ng8OmGqVatWSE4+492H3GqOHTvWx7c6mVvRroFoW6ny4DF9+nQMGTJEW1CQ+0z5Ii15J4DaH6WXTDlke1Ijswd58hQ1ve3bt1cKXP60qeLkRIuntLQ05eRnFBc9l3c0rNyiqkqnqEEkbRsJ9KKXLXkMs+rmzp8QKe5KyZ5rgDN23ICvFx8xfjPtGo0xqjolbavK5Ib6tr87NoxMFeT6o1080aSRMBKKVGGMzktYRdyhVN2mC/jWmZH21KyNEuT1TLWYlc9yqARMKl8a7+rr6zXNL42J5L2INOj0/xSvWL+kcKHxh7TGixcvVi4S6OJG1SV0VE60i1lfX6+ctzZv3qy5EpZ3XgC9oCqWPaWBzhKQkkU82wboFXLFxWc934k76dnZ2RgzZow2b/prIzTGqExkRLfJlB9RCSC3DXHHUR6fFixYoD1XuUMVoXimT5+ujfXUnlUuMEVEBYloJmR0fk81JstKFnF+9rdwiFQcCfoZGRnYunUrDhw4EOz0BIWkpCTlhV707Pzzzw95WBnSePv7U+0YOIWEI39u0URbQHFwoYGOBE754FB2dramLVRpc0SBni6hevTRR322eo20HTQ4p6amaq4KyXcwTfpk3gD4v6XPrJzo0Cxw9jAVaYMB31sdqZxat26tE8BJkBevVid7dnLXRRONeACY3qurq0Nzc7N2aAqAj8AHnLVZpUlBhC63Mropk7Zxv/76axw6dAgLFy7U8kjaH/qv6Aua0i6Xi3zwkkyLaAub0iteBS8LDeLZBHK3KPrbFzWqNPGo0rt582atjZEWjCYLldmDLPzQbx06dNDZeYttyp82VRQiaRFaWlqqFM7JPE3lAlfUugPqLW5x4ShrscxMNKgsAejGCPm2V9JQDxgwwDTP/oTYjIwMTfAS3QWKwoZZHLQAEzXvtNgdO3aszqyH4iEBke7YsGISYxf5fdUBZcLfAlEMI++YOUmXfHGZ2YJX/DagNsGRtcpZWVna7eGkAVYJTvJZDn9+/4GzLhwTExNRX1+v3c4LnBGoxXsPxHtfxBu5vV6v8v4U1XgInL24sbS01GfH69ChQ5pbYDoIKmuXVYdKxbFA5YRh4cKFPgsLMq1s3bq1qQkRnXlJSEjQ3RqemppqaAojojJRkc9Z+DOPEcdWVf7Fti7fnm7UDygemgMLCgpQWFioO8hvpFigsY2gxaJqHBfLQFyQ0IKS5jLAmgVCJONI0J8wYQJ+9KMfITc3F6+//jq++eaboNmWB4O+ffvi4MGDPrbbO3fu1H4PddhIQOxQZho38YAoDTKJiYk6wRbwNf0RJwrVip80aXTVOAm/69atA3BmO1++TIUgjcXo0aNRUFCg3VZKvoPp8GNSUpJWH2QLn5KSghEjRmDo0KGWbrQz6rwqDReVKZVTQkKCbjBU2RPTYS9xl0Re6IjejICzh6bkw4YECRTdunXTJk2jHSVxQJeFXuCM1lnU0siTeUlJiabdkQUycdeIzLIyMzN12/LyoS9ZgC0uLta80sTFxWlCtuhv30ijKiOaCIlbuEYTmDzpbN68GYcPH9byJJo+GWmFjMjKytIWoSLiN408TtBELt7EavQN0XzNygE1sSxJ6CLBlGx3ybPT119/DeCMIwKzvPubvPPz89GtWzetD4hpseIPnX4TNe+0A0e3QcvxkIBI2uZABHqrmB1QtpIGCkMLYjv2/aq45DshzBYl8i6HXA+qXRV58alazMrfMSsj8T4QSjtpqFu1aqXtdog7UmLb8GeqKgunYp7lxa+8ezhr1iztIKhKu2x2qFScY8U00cWFYr2IZqliWYuLYlnxAfjOV1bHK1pEA9DNxeIcLJoGUpuUF4qq/NNikKC5zuxyOoqHlJ4krNOiRjaBlctZHA+pjPzdcC8+Ly8v9/GoZkVhGsk4EvRvuOEGrFq1CsePH8eiRYtw1113ISsrCzfccIPyL5TU1dXhwIEDOH78uPYsPT0dzc3NePfdd7VnjY2NWLNmDfr164cuXbqEPGy4KS4u1nUoo44hH1oZOXKkTttCg5zK9EeemEm4q66u1m6bFG9UlDXvKjePhCjoiIczyWadNKXbt2/XtNrNzc1ISEhAaWmpZh9OWgQzRM2ULGDLGi4rE7E8mYjaTDGPF1xwAVwul5ZP8gCSmJiIiRMnauVNwmpBQYGPlqyiokKbNOVLfeS0yQK3+DvVnZGJhEogS0lJwcKFC7VBVV4YUjsSbXLFeqV0FBUV4fTp00hOTkZpaalusqfJ26pwLbtfEy9nMtM4k4mOeMiYNGWyFxejOFSLadWBTLEsqW14PB7lAVhxZ8joO2Z9lDDyokHtaPXq1ZpJAfUvccIWd8tkAYDypFqQivhbDJgJwma/yV5JKCwJiCrBywmqfMt1IZ9xUoVRxSte1kTjcXZ2tmZSYmbbbcdUUXzHSOA2KmtRIJR3LIzaZHZ2tk8bM2sHqt0D4IzLURrn6QyP6h1/fUGlvKE4aGEh39YulwcJ7fHx8Th06BBGjBiB4uJi3Q60Pzel4m6fURmYXWanUnzI85W/8YoghRwA3cVm4lhNc5/Yx8V5YMSIEVi1ahX69+/vM16TGVBiYiJOnDjh4+ZSZTa0ceNG3HvvvTqlGjkZEPNqphSgyxJzcnJM3ZBSecq7U4B6JzEacSTod+7cGZ07d0aXLl20/zf7c8qKFSuwbNkyTdu1ceNGzZb95MmTAIDdu3cjJycHb7/9tvZev379kJ6ejsWLF2PRokV49913kZeXh4qKCp/KClXYcCN3UqOOIQt/8uCvGvTlyYg6vqhRFn3oix5PEhMTtQmYEN08EqLGQvQaQDacpO0VIUFaPhTtb6ADzg4uEydO1DTOZE5BLhrFA1xUpoD57aAU96ZNm7B+/XqdMEwHmUVImy2mmfJTV1enxUv1KWoqZHMWVdoIygMNolR3dgY22h4X00jp8ng8WjsSbXKzs8/616aDxPLWtdjmjBaoRtBNpbRgFdNHB8BU9rE0KZIQIS7K5LIW/fz7O5gm7qbQN0WhW9Z0i31L1GaKgpX8HdU3ZIx2Dqg+6EB5YWGhZpaWmpqqCRvUB1NSUjQNl9nBYidCqBFGcZGZkdHBPjNvHnbJy8vDggULlItUsS7EMzr0zMyNoRyH/P9mQpHdviF/i9z20n/9IQqE8uKA4qXdJ1n4FtsYPaNFjL/7T+jmYuCM8mDWrFmGrjtlVAsxo11kCl9UVIScnBzT3UPa+QLgs0sllrFqgU3fEM8VqdIqz71m9S3OSdRnxfFKjJe+Q+UOnPXuI+4gi+5P5V0DEvxp90blp1/8FnBm91tUZqnilMtYVKqJd9pYUQqQp8GFCxdaWmQbmUZSWZopMSIdRzfjthR33HGH7qp7kWeffRbdunXDli1bkJubi6lTp+KWW27Rfm9oaMCSJUvw/vvvo6amBn369MG0adMwePBgn7hCFRYI3824pNEHzG+Ck2/fk7XRslBbVKS/jRI4e4tgdna2dgNft27dNM3u2LFjfW72VN3UR/HTt83Sa3Q7bHGxtdsuVflKSkrS2YoTcj7p/xMTE1FTU6O7SViOW5UP+p1u66utrUVTUxPGjBmD1NRUn3cpz/Hx8ejYsaNp+RjVbU5Ojq490OAq2jOKW7JyfYj/T95O4uPjUV9frw3e4kCsKl+5TVE7Em9WFMnLy8OqVavQpk0bjBw5UnkDpHyrYmZmpq6did9R1aU4qVm5YVIVn5V6F9NFuzB047DYHsSbY1esWKG8adLolkfxG5QmYvr06di+fTsGDhyIxYsX68YHurlavoVWjod2e+iwcocOHQxvNKX2JcZhlj65jI3GG7ms5d/EsnRyg6xROLqRk6DbjVXjpTzumt1wapRn+n85LlX7tjoeyN8S61J1y6jVsqHnFJ+qfuUyoHwA+pvIxRtI5dvK5famQm5jRv82uuVaHNv8jbXi+DRr1ixlOxPnIvkMiTzuGfUPuU3L/zXqb/7KhMyPrKZDVf/FxcXahVzAGZNT8eZdsbwJeZ7ZvHmz6c24VvutjHjDrVF/8Vfmcj7Mxq5QEqgcGdGCfiwQLkFfhDpjfX09xowZo+tM/iaZEydOaO7n2rdvrw3WdKW9LHDIgwjg28loUKLfSMARPcyI8YvXpAPQBk/VxCl+38yThthxRcGNoMmH4pHLxoqgarVujCZOGkBbt26tbUGaLVjkiaaqqgp1dXXatezioC4OruIiTBzExDIC4FNeJHyaTfD+8mj2LglXbrcbXbt29Zkw0tLSNEFAFprl8nj88cdx+vRpDBw4EBMnTvQRBqlOzQRDccFFuxRGgofZgvnNN9/UCd1yeauulCdBR7XAtrvQEBcQAHyEHiMFgZg+ajeygEt9lXz32xW0VW1OHA9kwUYWIEQhisYfEh6tTtaqcCTU0XRJgr4ZZosOf5j1GSdCj+o9WiTTuBvIOEbfMGo7RgIi3fYtLqLmzZvnM5dQuvwpsOQyT0pKQnl5uSZEqspBFEZpvqPxP1DBbsSIEToHFzQXqsY9o7r1J6j7W/TJczyVX1pamo9iwSwdgLFihfKmWhiIcoRqAQZAOWcbLeatIo9JAAwVBkZKKFGBYEVxGioClSMdme6oqKmp0dwnMpEBbUvRdjwdsBMRt8Dooo2CggIUFZ31nw9AJ7SLPoLlQyq0Ndq/f38fO0dReypuucsu/GhblQ5diRdsZGWdvSyroaHB0DbU31YbhaPBT0T04kEDj5h+AKau+ezY5Bpt6xcVnb2gh+zFVYMLCTaULvEZnX0gMxqqx6KiIp2JllFZiIfoxDCi3TnZq5rZXqvSCPgeUJO3lOn234yMDO37wNmDcmJbNrPDzsrK0vJfXl6u23JWmUwQRnXj8XhMDwbTrpL4rtjPjFwwUh7F9iR6xALMbzc12842a+/ixEH2wwkJCT6uRcX0yb7pabygvkquAsW0mJkBqmy7qf2lpaXptuLFA5iUVtFEMTExUWfWR+1EjNusj9J3k5KSdNv6mzZt0hYWZq6UCWpTqrKw+q6qzzg5VKyKT7xl1Mx1odktqWI5ZmWp3Wuq0p2VddbNoqxrJJMech0p31ugMhsTTUTEvk2OG8xMUcWxhZQ7ZFbjr63IZSBDZieyS9n27dv7nEuTTfBoHExKStKZ74lznGy2aWbKQqaQNGfn5+crx26zuGQTQOpvdI+KXAbUJuS7KsS+mJ2td38ryi2HDp259M3Ik5UZWVlnzjF07NjRsrMCWQlJYc3adjQQkKD/ySef4IEHHsDkyZNx66234tZbb8XkyZPx4IMP4pNPPglWGmOC2bNnIycnR/cnnisIBWKjJcHaTCCiQ6/19fU6u834+HjtcIvoMYU6q2gDSJ4wvv76a58bScWJT/SlK7rwy87O1gTcVq1aafbc4kEc0eZbnsBIICebayOPJeJgJtr6Dhw4UJsUxMGOtM8q22h5QKQ0GX2bfhc97dBlPzTI0WKIyt9ocFIJcDRw0t0FaWlpSE5OxsCBA5WCu7iIEdMrLgbEvKombdlGXJz8VGmUJ2aahER7+crKSk2opHDiQTma6ObNm+f3shuVRwh5svE3CVi96bSoqMjn7IJRncmCklyOWVln/eeTwGtXyFNp6MQ8HDp0SFdmcj2o2juNDzQ5ixO+VTd0/voUjSWiG1pCDCe35YSEBJ3bWapzWeAxsnemdl9eXu4TxoprSMJKmwr2u0ZCpxyf3O6ojamEKfmWVBFxLFOdt7GSRxKqSRgGoC0aZeUhtTOxjcntqLi4WOc1y196qF3Q2NK+fXudC1hVWxHLTy4DsYwobvEwqChsA/DxECf3P6PFSmpqqmn5Fhf7Hh5XeVpSjf1GyOMoCdPkfpsUG/L4T+VKygNAP7+KCx/KPwn4ALT8q/LkD1EpII+fct2qwopzstP+HG4cm+4888wzWLFihTahJSQkwOVyaYdkXS4XJkyYgDvvvDN4qY1Cwmm6Y2ebVzbvocEdOKPhTk1N9bF9J7MZ0aSBbLZpu13ccjQyE6LOJm7jAWf9CVOaUlNTNXOW+Ph4zJo1SxePvBVLmJmUiNpaMsNRmZJYNQmieM1MUlS23pROyr9oTuHUJEhOg8okwWhrVUynVdMDI1MfVdkbpUX8nnhmwt/WciBmEmZlKH7HX38S00BmHqJtPbV71VkJM1t2ue0B8GsupSprue2KNtCyeYS/shRt1ufNmwfA//kao3QBettswt/5IbHM5Wfilr1Vkzf5N9nsQ0xzS9nqmrVBwHjss2qeROMMmfmpzCHFdmJky2zFvMLqfCSPSaqx0GjsnjNnjrJPWcVsXDIqZ5rDRJNGwl/bE5VxqjM7KvMjs7Sq8gLAdA4IRpsWTQznzp3rUwcqcyOj78ttSixPwPgMkQiZ2rVq1crQ7JXCGZ0PCNR0KJiExXTnww8/xFtvvYXExET89Kc/xcsvv4xXXnlF++/PfvYznHfeeXjrrbfw4YcfOvkEEyB2bTmL/r+pSLdu3ZCfn69p2YEzGi5RSwmc9ZhSXFys8z5AW3Sie0qalIy0whQuJSUFVVVVAM660aQ4y8rKdB5BGhoalNpl4OyqXHUDp1g+ovcQAD7boSLkNs/flj2Vu5HHA7EcZHeYomYZMHcHKX7PyJWbvC2r0hKSdonqVtToqUxczNIhXiriTyMp1rmcftKmibepivHIJjWiRk91G7HVMpOR25ZV7yllZWU6ryCi1ojK28wbluwNRb75s6hI7RlJlS+qF5Wbzvz8fJ1LOXmHxswlHWn2xo4da0nrqdLEZWcbe3EBfHdQVH1d9V15y55+U+2emJk6VFZW6jSpchu3gtW2ZoScP7ktyXkn7anKjEJMj3xpoejCWEa+JVVVjmauII3yYpQ2AJo7ZpVrXKOxm84eHT582GcX1Eo9GNWv2RxDY41YBnRLuJW2l52d7bPbLvc/uR2q0qCC2gLdLWOklRbjcaIxB/Qmhqpxjf6tml+N5ASx/6t845u1Nbq75vTp06bmNuKutdxGxEVYtGryCUeC/ttvv424uDgsWLAAmZmZaNeunfZbQkICxo8fjz/+8Y9o3bo13nnnnaAllrGOPKj6G+jkzibaD5Oph7g17/F4NPMUonXr1jqhwao9qShQkiBPmjgj93n9+/c3zI+VrU0SlsjVF23VycINlRu5kpRdqaniFe1yzYRklUAlbyOrJmzVtrG/BYH4PbE+aAFDqMwk/E0olG/xvIa/ulctIuS8iMIkYOw6kuoSOCMM+3OnJtvPi78FIpSpJn/ZtESc2On2XrmcaIFD/xVv/hQXUfLOkqotiAs5VT8yu7zIDLLvpQlSdSOs2QIHUC9GCSu7J/620+V2a6evAL5mDnIbt4L4TSfty0hgSklJMRRKzWyJ5fFJdp9qZJbmz+zJigmILFCauSUWy1peRBiN3VlZWZqAR4ogVdxGWKlflcmhXEakDBLNTY2+v3nzZjQ0NPg4yTAqNxErYyydJaJx1t/ZGVU/tYJZXxa/oVIeyPlQlbF4XsfKxYl0P0nr1q1Nb9MVy1auI/rN7H6SaMGRoP/ll1/iBz/4AXr27GkYpmfPnvjBD36Affv2OU4c4xx5cDDTRqom1ezss/7OS0tLkZWVhbS0NNTX16OhoQEXXHABAOguvEpISHAkNIhpli9CEW35acCfN28eKisrDQdu6tCUZ1UYKp+5c+ciLS0Njz76qPKwGXV+unjM7Xab3pIqTsR27fpkrQ+gv2lYTpNKc2IUzgjy+210eQtgPKFQevPy8lBVVWXoz9yIvLw8nfZNzovo11hsv7Kml94jUy/5AiW5TIysFf2Vlz8bfSsL26ysLHg8Hu1GTLEMSZMmHkIW80cHzAFfIUMMJ2vLqG6pH/l7R8RMOBXLS7VgMFvg+CszO7snVg4hWsmrjJwn1fv+NKDieKC6iM4fch7EBbKRUGqUT9JYx8fHo6qqCqNHjwYATQDzJ6gTqkUdYL6okecYVV+T24u8aAagG1Pnzp3rI/AZ3U4rl4k/oc+o3VO66cCwqt5pnBDdDRvVidH9FiJWFWYicvqttv3sbN8zEFYR64nautmt3kbIbaO4uBirVq1Cc3Mz1q1bZ2mxTIfn//vf/2rzOzkbkReoKoWMkUlgIIqgcOJI0G9sbNQmVTPi4+PR2Njo5BNMgNgZHIw0NHSTLXCmka9cuVLTau/YsUP7jS68srJ96y/NdEGQfDOvOODLnVKGhASy7zPTMJAWyGiwpe80NTUBOLOwqa6uxsqVK029YphpUIyQ60G8UVaVJrlMRKyaGVBc/jR6ZuldvXq1zsTKDHGwJDt2MhsoLCxEVVWVdtjbKlQGF110EQAgLi7OcLKmHYz4+HifvPoTBpxMuIDvgUWxPqmtkoep6upq7N27F+3bt9cWAtTe/d10qUqf2I9FDaPZO4RoEmW2WFYt0lRpExftgUICyXfffac7yEjpVglhdutPzpOR6ZB8iZYIvUMHi1V14ASz8jbKJ6WV7k2Q0yy3U9kEjJ7RAqi0tFRXzuICQB7LjLSlYvr91Q/FsWrVKk3xImN0SZoYt9iuRUHU366FmG4AtjTfRnkjE9lWrVqZ7kKqBEyz51YdR6jSaUVjbgbtuADQLkmzIyCrdtJIOSNeGmkV2auf2RhFi10al8X2YUVxFqk4EvS7d++OrVu36vzDytTV1WHr1q3o3r2748QxzpEHATNtpNGkIb4jN+7WrVsDOOuRx0zotJNO2dOGlW09M6xss9NgK55LIOj7tAARDwaZTdp2tYcqwVz0LiRr+/1p32hS9+fy1oqZE6VP1tCRFr9r164+t8n604gVFRVpblLbtGmjpVd0pSoit0XVgEumLk1NTT6TtWyC1bFjR2RlZelcB8ptzarg46/cqJzIbSxdMy/2RbqRF4BWBmI50ORpdGOxmWZZpWG0kgd/3xTLS9UmrQh3ch4oParxSu4D5F1HdhssCt+BTMpWxrNsaRdSTidBi7s2bdoExQTA3y6b2e4CmTUY/Q74ng8S61Il7BYWFmoCFV0CKI5lVhZNYvrF3VjRTIvswOVFih1tqyyI+lvEilC6ybRS5VZSbvdm7knJRLaxsdGvlzZRw222Yy3mr6qqyvICQcZOmYphqexEL0p2tPuqnTRqr2Q6bHQGRQXt8owZM8Zvf6ayBuDjctbunB5JOBL0R4wYgePHjyM/Px/ffPONz+/ffvst/vjHP6K6uhpXX311wIlk7CMLOfI2lIjRoCtr0BMTE7UO17ZtWyQnJ2PWrFk6TYndgzxWtdhyeLNBQ3TBaNYpaQDeu3cvAF+/5iKkKSJ3mP6Qt5v9lQcJJ+LCxEi49adpBc4edCSNilm9WIlPpaEjLX5FRYVOeFaFF9NFgyX5mJ41a5bWvozcM1rZzZG37cVw4gAuvktbwirXgVYFHzPEchIFfMoL1fGsWbO0xWRcXJz2vuw728jjk5lwa6bh9deGzL6pSoOZYG9VY1tUVGSoPRf7fna270FG+q5T8wNAb5Imuz6UycrKwty5czUFQHFxse4+EkJctDtBtdD2t5CW36c5ID8/X3dYXMyL0cFaqkvS1mdnZ+vGWQDaohCAz1hmVQmk6qf0rLS0FO3bt1e+58/US0QWRP0tYlVlTZpv+cC3GD/FSwugVatWKc2FaF5VLTpERYGo4T50yNi/PI39wBmlgZnppxn+wonlIvbNzZs3AzjrVIEO6hotqmSoTEiYpz5G46Rdf/ZGuzyqvJBZWE5Ojs8dL2YyVKTjyL1mfX095s2bhy+++AJutxsXX3wxunbtCgA4cuQI9uzZg+bmZnzve9/DH//4R01zdy4SLveaYsOkTpic7Mw9FMUl3gQr34wJ+Lrzom+bdQ65A/m7Lt7M1Zscr7+b7IYOHapzF+rvllPR3SPhr0ytuuX0N5Co6tNqnKIbT1V6raRRTp9YvmlpaT43JIvhAV83gFbqxyhPTt9Tla94c6VRXdK7/m7SNfqu0a3UMqLLw5qaGr/uWQMtz0AmL6M0BDIZym1KdbMuuRUN9BZXM2RXgATd2Gr2TnJyMg4fPqyde9m0aZOW9kDKxsjlqtw+jL4TqAtFeSwx+q7qNnM77cOsr4nuFmV3m/7mDav5U6VNLDu64NHspl0Z+VZlf32a8isrKOQbos3cPhYXF+tukJdd/Fodu4zGfFH51Lp1a929FXSrsFgXdtu/WXu1M1dabXtWXZU67T+BEqgc6diP/qlTp7B48WKsXr3aR/vq8XiQkZGB6dOno23btk6ijxnC6UefCNYko/JnKzZ8o4GAwljpoGZCkZ2JngZ+wFiAE4U8swlC9ssPnPH7bcVfu9WFiR2s1qc8+QLWBXG7mA2EZr6bAdgS2lS+oY0wKic572YCstMFs5P3rNaJlUlHJTAFUygO1cRn5m8bOJOvBQsWaP82E76dMn36dGzfvh2JiYm6ey7i4+Oxfv16v4LE5s2bDX1zq7AjAItCl50+a9TmVQsqFbLffX8CE9Wf2LfttBmr/uyNysdfOchjothHs7OzdWMCcLYvivPPxo0bTfOgyg+gvzNCtUCXffTL9wXYmU9VYxD928rYoIoD8N1tEaG6D0TQD0RmkdsqLc5UcoOoxCsvL0f//v1RWVlpOmeEQ6MfNkGfqKurwxdffKGZPSQlJeHiiy+2dFj3XIAqqGfPnrpr2QFg/PjxGD9+fJhS5h+VsEAd28qFTvLk4u+SHyMhTrWj4E+TTdoMusDDaCAUtR5Ggj59/+jRo9rhcqsChriY8PeOE+HbqibKivAdCFYmYfHil9TUVMOLuszi97doEQdwo7pVXdBiJOyLl8GR2YUVjbnqQh0rk4QToc/s+6RdM+oHToWnYKVThi6w6dq1KyoqKpTCsrgTFwoNm3gZmAhp6O0stKwIDKHUFpp9D4DhgspqfozCicIVxW3lciKaa1T93E6ajBa6qt0auY/KAq1YLmZ5MENUUhkpPGhONbogy0l/MpvDrdS90aJDrBNaGIu4XC6dss5OG7fTd8zeF+UFsT+ryl8OQ4v6SCEsF2aJxMfHY9CgQbj66qtx9dVXY9CgQSzkK3jiiSdQWFio+4tkIR9Qu7DLzjb2LSvbMop2jkVF5gf7AGPbWnHgpau0ReTvFhUVaRd3GQn5BF3KYuZthvLdtm1b7YyC1YNF5FHG5XJZtk8tKlJfhiNDk4cVryhmvxUXm5+tMLIHFjGzwaXfduzYodmqkn06lb8qjeJ3yQa3tLRU8woBnBmoySZa9iCi0mEUF/seei4qUtu3FxcXa2cy6uvrbflQF8tXVTaqMjWqT7N+5e/7ZDapcjnqr/2YtUEraRDft9KGgLMH8SoqKgxd9Y4ZM8bn8Lc/rH4fgI9rxuTkZO0wH2DtUJ5R2ameW4nPKUbfo7MN/fv3D8m3U1NTffyqm93ZQOkkz1JlZWV+bbHN2pccH8VBZZ2RkaGba8T2bDQPFRef8Tg0Z84cTYA1O2grkpWV5WPznZmZiaSkJK0tA2fOApWXlyvLyeo5B0qreGBXnMPpUDPVverSQkIex1QeeRYvXoxPPvlEu+8E8HXla6eN2+k7ZqSlpWl1bXQmQy4LwuwCuWgkYEGfiV1UndNMaLHit95ME280kIiLC9XgLw8AVr5F71nxTCOG7dChg62DRXTAkoQEI/z5lDdKE2khVP7rrQjftAgz81Jid4A1QvSyQ/GS8Az4Hlr29136XTUoi8IMuWqVv0m7AeRbXJ7YCwsLNW04uR+1KhAZlb1RXxE1foC+PimfdnxS0/dnzZpleDDV3+LbyF+6GaLAJbZhq22IhCyPx2O4AExNTUW3bt0wa9YsQ02vvHC1c1iTboKlHdiKigqdcGdnoSWn39+YalVwpHz6W7zI3yNtJ93lUFlZqfRqJsdptf4oXGFhoY/21Uy5IAtcqjFQTpuqfZGHLDpUKd/ESoI83ctCZjtifo3mIbkMyCzTnx98MV7ZfWd5ebnmyIAOQRvdzk6I5WD0/9TeGxoafOKiBRfVvZVLHf0pKYCz/cbOfRkq7PQdFVS2tLihS7pUl3XJZUE3hQ8YMEBZrtGKJdOdl19+OaCPTJ06NaD3o5lIsNEPNVbtN/3h1BzAqf2c0Xaq2XetbCk7gcwF7Nh+2km/v3jMbHadblerviPWkz87dn+mTKptaTGcvzj92Q0HesBPhVFfEc+UAL6ma3br2YoZmL9+Y/fAp9N35HRbNfEz+l0sS6o7q3Upl5sV075gYmccsGoOoepn8sFOf/3BqdmEVTNNAH7zYuWApvxdf6ZS9F0r54Tk8YYOXQNnLBuMFp5W4rIzd4nplw8oq57L7dZsHLZjjmp0ZiVY8oBTzEwxVWFVeRfHCyrXcOUHCIIc6bXAhAkTvBMnTvROmDDB8t/EiRO1v3OZkydPejMzM70nT54Md1JCxvLly73jx4/3Ll++PKD3xo8f773iiiu848ePD0UyLaXBX9ihQ4eGJI25ubneIUOGeNPT022Vo1H67eaLwqrqoCXqxWkbkhHTaiVOszDBSpNIbm6ud+jQod7c3Fzd8/T0dO8VV1zhveKKK7xDhgzx+d1uWsRycFp/Rmk1i09Mp5Pyo7iHDh1q+J6/eEeNGqWV5ahRoyy9Y5Q3iovicZIeOxiVudXvqp6JeZLjl/PrdNyQn48aNco7atQo0zoUw1C6brvtNsv90Ur+VW1VbqPimG4l/xRnenq6d+jQoVrftfp+oIjfENunWd9zWq+q94YMGaL1r/Hjx2tpkNuQ/DyQfNrFX737QyxXaptiPbc0gcqRljT6drcsjh49qnnjcblceOONN+yvQGKEWNPoB/P0uepQZLDitoqdb9Iqn0w5gunNRPTOE4yDQFbLVtagktYvGDsmZoSqrv3FK+5OpKamtnh7M9JKmmn0gcC8VgDOvCrZ1e7Kz53siFD9GHm/sIIdjZ7qXVHTamUHL9yu91RpEbWPgO8uF7lIDmR30ijfTjS6/rwt2fm+iN1dYtUuoL84xTGFzGCsegsKBDvt3KnHI9XOAOXX4/Ggrq7OsJ4DzXcg/SrQQ+9mu2Axq9G3SmVlpfcf//iHd9KkSd4JEyZ4p0yZ4l26dGkwPxF1xJpGP5ja3ZbQfvjDTn7EVX6wtdykVSKNbqBY3S2xokENJv60PaFuE1TOQ4cOtVSHwU6PPy3o8OHDlRr9QNqb3V0fep6bm2vrPVn7bUUbLiO2x3C0DzkdVnZEImEck9Ni1L9yc3O1cSbQ3UmjfJMG9Morr7Rc//7anOqboSh3q3Ea7YL621EJ9LsidsYEpxp90mKnp6f7hLvtttu8V1xxhfe2224L6vcDeScU8UVC/w5UjgzKYdzjx4/jueeew8yZM/H222+jdevW+PGPf4xnn30Wt9xySzA+wUQIVg/EWMHOAZ1QYSc/4k21Zu85Obwjevrwd3DXCnLZGqWXnofq8iGRvLw8LFiwwMcXs4h84M1ffFYPLhLi7blW6t5OegKBDgB27tzZx2MFEFi/M8qDv+crV65EVVWV5fhkxP5iFdWhTCffDhSxvOWyl9udE28ooTrYR4dNASi9xtChUeoDchn7S5/4u1G+6YCjeImSlXQbHZgUCXX9UzoAX+cAYt7FdIjtQ1UmoRpn7IwJZm3UzKuX6kZnCr9jxw4AQHl5ud/vO8lfsOUDK/GZtf8333zT9nwTMQSyyjh+/Lj3+eef906ePNk7YcIE7+TJk70vvPCCt7q6OpBoY4pY0+gz/gnHWYNgYceu3Z/2TUbctXBqf62Kb+jQoZa+74Rga3PCoR12otGnerKjVVfZlwdytsBOXlqSQNpdqLSw/r4h2qIPGTLE0H7eX/r82fSLtvd2xwcrGGnSncZhJ24jzX0w2mQgNvWBftPMll6VDgo/fPhw3bkKs/oOJD8t1efFPkK79mJbEHfCWpqwaPRramqwePFi3HnnnXjttdcAAJMmTcKzzz6L6dOno0OHDkFdjDCRTyy4oAoWwdz1aAmMNFVGiBpfO1oa0qaPHTvW1MWqVS0Oxde/f/+Qtb1ga5X8tY1Q7HIZxWn2fODAgQCAuLg4ZGdn29LkirsRdjV5VvpObW0tCgoKDO99CBbFxcUYMWIEhg4dqtPiibtCdhHzR2Wal5cXkEtLs29QPhYsWKC5bu3QoYOhK11/5S//LqexSHBfS9p5wFc7LuNk/nA6zlrZ4VLFbaS5d1pPootRALo+JccZivmVvgGod1cB8/EoISEBGzduRGVlpfLOAqvxWE1nqHfxigR3wwB82sLAgQMd9/uwY2dVUFNT433xxRe9U6ZM8U6YMMGblZXlfe6557xVVVWOVhnnArQS+9nPfua96667dH8rVqwId/KCRjRrse0SCVrFYOLPxlRG1Hy0lG2/GUZtL9B6ioV6dmobK9evUzvjUO2IGO02BBPxW6HQ4oXqPILR+RzRo5NVzz52vyVq9OXvm9WV1fkjGPOMkTbezm6mURu3U2divfire7tjnFEdB2Mnwmj3IRQ7OKrvhYpA2kWoCVSjb0nQP3nypHfp0qWagH/TTTd5n332We+xY8ccffRcIppNd+w08GBOHpFOrC1qwnFQKpiDp1FcgdZTqOvZahmohCerOMmDSphuycnO7Ft0AJCE1lCkRxRcxMOIoRJgnApI/to9LdKWLz9zqJKE/GAtjI2EIX/mPXby4jScVez2DwovmnY4jc+OS2W7Y5yReVkszF0tLXhHQpm1iOnO7bffjmXLluH06dOYMGECnnnmGdxxxx3o2LFjiPcbmEAIdLvPzpaZ2fXmkYzTLePExETU1tbGhKmSk23VQE1Mgrkda5SWQE2oQm2CZfW2VjKHMDK3MMNJHuhA7MCBA00PGYYKs3KhA4AAMHfu3JCkh9pmWVkZ1q9fj+TkZNTV1dkqe7NxhX7bvHkzgDM3/TopW6M+lJ2dDZfLhebmZhQWFiIrKwvr16/Hpk2bkJ+fr7WJlJQUXRrt3B4sf19l9kLxA/CbP6P2ZXZQNBjY7R8UHoDStIZu5LUSX1lZmXYrrr/82B3jjMzLVOGjzey2pUx5CCc3hEcalgT92tpaAEBjYyPeeecd3H777fi///s/y39MeAi0Q9gZBKPNLp1w6g0gISHBkeAVjYRiImiJ9hKoUBAOr1CqsqaFpexFxQpO8iBfCx9I/oPddtq0aQPgzF0ToaqXlJQUuFwuVFVVobi42FFbJaG5oKDAJ/9Oz7iIFBcXo7a2VtkmsrKyfM7Jqc5YkM/3YMwPonBPwn6g8QMt52XHbltKS0tTnlUoKyvzic+oDwRjDDRKf35+PjZu3Ohzkzl5ZSoqKtLSY6eMI2FR0NKyRrQqMUUsH8b1njHzQVNTk+0/JjwE2iGsalnMwkYyZpOlP6J1YeOEUEy20dhego3K/aSqrMn1ZklJSYuUVzC1fk7ajplbzlmzZiE5ORmzZs2ylQ47kKa1rq5O04g7bav19fU++Rc1g26329EYIh56VaVLLkNVPcj1bNcdqlguqsVDMMbISBtnjQR6MZ1yXzHqA3ZcXjpF5YJYTo+dMm5pbbqKlp47Iq0NOiKIZkSMgmi20TciEmzWgkGs5CPURMJhJCtESzoJq4dXIyFfwXZlGMksX77cO2TIEM0W2+jwoRWbc5X9PZWlmZtLK2lUndsw+q5RmgOpH1UaIrX9BgsreXFyRsFfHE5R2eoHWudW342leg83LXIYl3FOLAr6sdKBYyUfzBlCsXALZRsxS6/Z4cZwcK71FdG5gHzDbyCCnHgzrRXPQWZxq9qFeBA31F5szLzG+PtGLLenYOQtWOVj5/baYBMJ45ZINLe5iLgZlzm3iBWzi1jJB3OGQLZYjbbKQ7lVbZZe1eHGcG4dt2RfiQQ7YDO7XLk+7LQRupmW8Ge6Yxa3ql3QM7r1Vj5wayUOq1g9O6L6RiSYgIQKua84ac/B6m+VlZW6/7YkkTBuicRym/NLkBcejEQsavSZyCKaNRWRQrBcDgZKqH1SRzrLly+3rI0OdTqs+tF2qtEfPny40nRHjM+u22I5LcEwIwkFkZIOFeF24RlMIrmcW5poLotA5UiX1+v1hnuxEcvU1tZiypQp6Nmzp3bjGjF+/HiMHz8+TCmLLehW1+z/7wow3LRkejIzM3Ho0CEkJydrN1FGWnlEOpFSXqq6PJeg/LvdbsyZMydq265Ze6I8EvHx8Vi/fr3P78nJyaitrUV1dTUSExNRUlLiN365/ADowp3r7csKwS6jSBlbmOiF5Mhly5YhISHB9vtsutNCPPHEEygsLNT9sZAfPCJtW64l03OubY2Hgkgx44q07e6WhvIfzUI+YN7/yC83UV9fr/vdShswij87Oxtut1vzoS8LmOd6+zLDiS98K0TK2BJsIsHEjrEGC/qMbVq6g1v5XqRNYC2ZHtVEEmnlwVgjVoUCq4Qy/y05blH/U9nIk/0/ALhcLowZM0b3rlgGRi4vjfp3VlYW5syZo7zUSY47lAS7rFui7sx84TO+sDIpemBBn7GNqoOHciC2MqBEmoAU7vSE+/uMfyJZIxbJaXNKSwkmoqlGaWmpz22zpHUHgG7duvlcaiRi1I/N7jgpLCxEbW2tz6VOLYm/srbbvlqi7iJJORIN/S+SyosxJ6Jt9BsbG7F06VKsXbsWNTU16NOnD6ZNm4bU1FS/7+7ZswcvvviidmV6v379MGPGDFx00UUtFhYI3LYqElHZHIbS9pNtHJlYJJLtpSM5bU5pqXEkEBv7YH0bQFjrzl/+7Lavc20OiMX+xzgnpm30CwoK8Prrr2PkyJGYOXMm3G43HnroIWzbts30vT179mDu3Lk4dOgQpk6diptvvhnffPMNfvOb3+Drr79ukbCxTEubirB2uuWIBk1SrBDJGrFITptTQj2OqGy8ZdObvLw8DB06FAUFBSERWq26vAw28rjhr6zttq9Q1F0kj3WR0v8iuYwYGwTRA1BQ2blzpzczM9P773//W3tWX1/vnTlzpve+++4zfXf+/Pnem2++2Xv8+HHt2dGjR72TJ0/25ufnt0hYgt1rRi7R7G7LH07zFmmXnFgllusyHASrPM+lerHSd8h1KF1qFSvlEo3jRjSmuaXhMooMYvbCrA0bNsDtdmPcuHHaM4/Hg4yMDJSXl+PIkSOG727btg0pKSlITEzUniUlJWHQoEHYtGkTTp06FfKwTOQTy4eJnOYtUjRJdonlugwEpxq5YJVnKOsl0rSNVvpORkYGXC4XAKC5uVlZLpGWLytE47gRjWluabiMYoOIFfT37t2Lnj17+tgjXXLJJQCAffv2Gb7b2NiINm3a+Dxv06YNTp8+jf3794c8LBP5xPIg5jRv0WomFct1GQjhXvA5iceqoBtpizsrfSc/Px+bNm3CvHnzLN2KbEYkLQiiddxgjDnXzkXEMhEr6FdWVqJTp04+z+nZ0aNHDd+94IILsHPnTjQ1NWnPGhsbsWvXLp93QxVWpra21tJfY2OjYRxMcInlySmW86Yi1vIbLCEuGhd8VgXdaF7cmZWv1XxF2kIn2oi08oukhRugL59ISxtjj4gV9BsaGhAXF+fz3OPxaL8bcf311+PgwYN4+umn8dVXX2H//v148skncezYMZ93QxVWZsaMGZgyZYrfv+XLl9soJYaJvAmCCZxgCSHhXgA5yYdVQTfceXOKv/5qNV/RvNCJBCKt/CJt4SGWT6SlrSWIpXm1dbgTYITH41Fqt0mYJoFfxXXXXYcjR47gtdde01yafe9738OkSZPwr3/9C/Hx8SEPK7No0SJLbpFUixuGMUMchKNF6Im0beFISw9NroEKIeHOl5N8ZGVlRUQdBIJZuQerv8ZCOYWTSCu/YPX5YCGXTySlrSWIxnnVkCAfDg4av/3tb7133XWXz/PPPvvMm5mZ6f3444/9xnHixAnv1q1bvfv27fN6vV7vCy+84M3MzPTu37+/xcKy153IJha8gljJQ6TlM9K8OURaeoJFrOYr0jEr90jriwzD+BJJ/TRmve707dsXBw8eRG1tre75zp07td/90b59ewwaNAh9+vQBAHz22Wfo3LkzLrjgghYLy4SeQLbYYmFL0spWf6TlM9K2zcORnpbYGg53OTvJYyxsmZuVe7SaHDHMuUQs9dOIFfTT09PR3NyMd999V3vW2NiINWvWoF+/fujSpQsAoK6uDgcOHMDx48dN4/vwww+xe/duTJw4Ubt+vKXDxhKRNBkHIsSGWxBqKSItn5E2iIYjPS2x+Ap3ORcWFuLQoUMoLCy0/I6Vcomk8UeFv3KP9PSfC3AdMOcKEWuj369fP6Snp2Px4sU4fvw4unfvjpKSElRUVOCee+7Rwu3evRu5ubmYOnUqbrnlFgDA1q1b8corryA1NRUdOnTAzp07sWbNGgwePBgTJ07UfSdUYWOdSLJfC8S2MVR2muG2jZaJNHtUJvJsciMFK+USSeOPE6I9/bFApNVBpM0ZkZYexjkRK+gDwOzZs7FkyRKsXbsWNTU16NOnDx544AFceumlpu+df/75cLvdePXVV3Hq1Cl069YN06ZNw4033ohWrVq1SNhYJ5KElEgUYiNtEmEij0hst8EmJycnJIdxI2n8cUK0p9+MaBEQI60OIm3OiLT0MM5xeb1eb7gTEcvU1tZiypQpWLZsmSWvO0xsEC2TXSwQrWUdrelmGDMyMzNx6NAhJCcnY8WKFeFOTtQQaeNBpKXnXCZQOfLcMSpnmBYk3LbR55L9aTBs3cNRXpF2QJphgkGknQeKFsI9Z8hEWnoY57Cg30LMnj0bOTk5ur+333473MliYpRzSYgMhmARjvJigcg559JCNtoIpoDI9Rw+uOxjBzbdCTFsumMMbw2GDi5be3B5RRdsHhJ9OOljXM/hg8s+cmDTHSZqOZe0zi0Nb7vag8sruuDdkOjDyXjP9Rw+uOxjBxb0mbDBAwnDxB4tseXPCzNzItHswsl435L1HIllFk64j8UObLoTYth0h2GYcwne8g8eTk3KuA7sw2XGRCpsusMwDMNEDLxTFzycmjdyHdiHy4yJVVijH2JYo88wDMM4gQ+JBwaXHxMLsEafiVrYJpJhGMYYtpMODHb4wDAs6DNhhAdhhmEYJlSwOQ7DsKDPhBEehBkm9uCdOiZS4B0RhmEb/ZBDtlU9e/aE261fV40fPx7jx48PU8oYhhFhe97gwN5LmEiB+zQTCwRqo986BGliFDzxxBN8GJdhIhjRlIyFAuekpKSgoqICKSkp4U5KTMNCrH8KCwtRXV2NwsJCLiPmnIVNdxiGYcCmZMGirKwMzc3NKCsrC3dSYho+48QwjBVY0GcYhgHb8wYLXjC1DFzO/snJyUFycjJycnLCnRSGCRtsusMwDMMwUUZWVhYvSv3AZcQwrNFnmHMC9oTiHy6j4MAmJQzDMJEDC/oMcw5QWFiIQ4cOobCwMNxJiVhYQA0ObFLCqOCFNMOEBxb0GYZhwAJqsOCzDowKXkgzTHhgQZ9hzgH4UBrDMOGEF9IMEx74wqwQwxdmMUx0wBc9MQzDMJEGX5gVJfCFWb7whS9MJJGdna21R4ZhGIaJBdh0hwkbbLPJRBJsW84wDMPEGizoM2GDbTaDA3uzYBiGYRhGBQv6TNhgDWpw4J2R6IQXaMy5BLd3hgkPLOgzTJTDOyPBoaUFEV6gMaEk0gRrbu8MEx5Y0GcYhkHLCyK8QGNCSaQJ1tzeGSY8sKDPhI1I0zhFK5E2oUcrLS2IsOkaE0oiTbDm9s4w4YEFfSZssIAaHCJtQo9WWBBhYgluzwzDACzoM2GEBdTgwBM6wzAyvGPKMAzAN+OGHL4Zl2EYhmlp+KZnhokN+GbcKIFvxmUYhmFaCr7pmWEYgE13GCbq4S16hmFk2KSPYRiABX2GiXr4UDPDMAzDMCoi2nSnsbERS5cuxdq1a1FTU4M+ffpg2rRpSE1N9fvunj178OKLL2LHjh0AgH79+mHGjBm46KKLAgr7zTffYMmSJdi+fTtOnDiBLl26YOTIkfi///s/xMfHB5hjhrEPb9EzDMMwoaK4uFibY3iHKPqI6MO4jz32GDZs2ICJEyeiR48eeO+997B7927k5+dj0KBBhu/t2bMHc+fORefOnTFu3Dh4vV68/fbbqKmpwZ/+9CdccMEFjsIeOXIEv/zlL9GuXTuMGzcOHTp0QHl5Od577z0MHz4cv/3tb33SEughCoZhGIZhmHDBB7vDS8wext21axc++OADzJgxA5MmTQIAjB49Gr/4xS9QVFSExx57zPDdpUuXwuPx4LHHHkNiYiIA4JprrsHPf/5zLF68GLm5uY7Crl27FidPnsTChQtx4YUXAoC2OCgpKUFNTQ3at28f9LJgGIZhGIYJB7xrHN1ErI3+hg0b4Ha7MW7cOO2Zx+NBRkYGysvLceTIEcN3t23bhpSUFE1wB4CkpCQMGjQImzZtwqlTpxyFra2tBQB07NhR971OnTrB7XajdeuIXTcxDMMwDMPYhg92RzcRK+jv3bsXPXv29NmmuOSSSwAA+/btM3y3sbERbdq08Xnepk0bnD59Gvv373cU9rLLLgMA/PnPf8bevXtx5MgRfPjhh/jPf/6DzMxMUxv92tpaS3+NjY2GcYSDxsZGvPTSSxGXrkiDy8kaXE7W4HKyDpeVNbicrMHlZA0uJ2tEQjlFrAq6srISnTp18nlOz44ePWr47gUXXICdO3eiqakJrVq1AnCmsHft2uXzrp2wV1xxBaZNm4Z//etf+Pjjj7XnP/7xj3HbbbeZ5mfGjBmmvxNTp07FLbfcYilsS9DY2IiXX34ZN954I+Li4sKdnIiFy8kaXE7W4HKyDpeVNbicrMHlZA0uJ2tEQjlFrKDf0NCgLBSPx6P9bsT111+PwsJCPP3007jpppvg9XqxbNkyHDt2zOddO2EBoGvXrrj00ktx1VVXoUOHDvjkk0+wfPlydOrUCZmZmYZpWrRokaVDFNxhGIZhGIZhmGAQsYK+x+NRbnWQ4E0Cv4rrrrsOR44cwWuvvYaSkhIAwPe+9z1MmjQJ//rXv3QmNnbCfvDBB/jLX/6Cf/zjH+jcuTMA4KqrrkJzczOKiorwwx/+UGfrL5KQkMBedxiGYRiGYZgWI2IF/aSkJKV5Dmnazz//fNP3p0+fjkmTJmH//v1o164d+vTpg8WLFwMAevbs6SjsO++8g4svvlgT8onhw4fjvffew969e3H55ZfbzivDMAzDMAzDBJuIPYzbt29fHDx4UPN0Q+zcuVP73R/t27fHoEGD0KdPHwDAZ599hs6dO+t849sJW1VVhebmZp93T58+DQBoamrym6a3337bb5hIezcQnH6Xyyn073I5Rf674fhmNL4bjm9G47uBEG3jTLSVU7je5XJqmXfD+c2IFfTT09PR3NyMd999V3vW2NiINWvWoF+/fujSpQsAoK6uDgcOHMDx48dN4/vwww+xe/duTJw4EW63ebaNwvbo0QNffPEFDh48qAv/wQcfwO12a4sEM86lBhrId7mcQv8ul1PkvxuOb0bju+H4ZjS+GwjRNs5EWzmF610up5Z5N5zfjFjTnX79+iE9PR2LFy/G8ePH0b17d5SUlKCiogL33HOPFm737t3Izc3VeavZunUrXnnlFaSmpqJDhw7YuXMn1qxZg8GDB2PixIm679gJO2nSJHz66aeYN28exo8fjw4dOmDTpk349NNPMWbMGL/mRAzDMAzDMAzTUkSsoA8As2fPxpIlS7B27VrU1NSgT58+eOCBB3DppZeavnf++efD7Xbj1VdfxalTp9CtWzdMmzYNN954o+ZC00nYSy+9FI899hheeuklvPPOOzhx4gS6deuG2267DTfddFPQ888wDMMwDMMwToloQd/j8eD222/H7bffbhjmsssuw1tvvaV71r17d/z+97+39A07YYEzF3bNnz/fcniGYRiGYRiGCQcRLejHAl6vFwC0Q8XNzc0+B4ytEo536Z2W/i6XU+jf5XKK3He5nKwTSFlFW14DeZfHcmuca32Pyym07waznEietIvL6/RNxhLfffed5VtxGYZhGIZhGEZm0aJFPu7drcCCfohpbm5GZWUl2rZtC5fLFe7kMAzDMAzDMFGC1+vFqVOnkJSU5NdrpAoW9BmGYRiGYRgmBolYP/oMwzAMwzAMwziHBX2GYRiGYRiGiUFY0GcYhmEYhmGYGITda8Ygu3btQklJCT7//HNUVFSgQ4cO6NevH2677Tb07NlTF3bPnj148cUXsWPHDgBnbiSeMWMGLrroIkfhAKCxsRFLly7VXXQ2bdo0pKamhijHzti/fz9efvll7NmzB8eOHUObNm3Qu3dvTJo0CcOGDdOFtZMnq2GjpZyA0JTVqVOn8Oqrr2LXrl3YtWsXampq8Ktf/QrXXnutz/ejpaxCUU52+vO5XE6h6s/hpCXytGzZMixZsgS9e/fGX//616DE2dKEopy2bNmC3Nxc5fcee+wx9O/f33ackUAo29SePXvw8ssvY/v27WhoaEBycjLGjh2LiRMnOo4zXISinJ588kmUlJQYfrOoqAjnn3++rTitwIdxY5A//vGP2LFjB9LT09GnTx9UVVVhxYoVqKurw+OPP44LL7wQwJlOOXfuXHTu3Bnjxo2D1+vF22+/jZqaGvzpT3/CBRdcYCsc8dhjj2HDhg2YOHEievTogffeew+7d+9Gfn4+Bg0a1OLlYcQnn3yCt956C/3790dSUhLq6+vx0UcfYdu2bbj77rsxbtw4LaydPFkNGy3lBISmrA4fPow777wTXbp0QXJyMrZs2WIo6EdLWYWinKz2ZztxhptQlFOo+nM4CXWevvvuO/z85z+Hy+VC165dfQT9c7mcSNCfMGECvv/97+u+N3jwYJx33nm244wEQtWm/ve//+Hhhx/GxRdfjBEjRqBt27b49ttv4fV6dS7Go6WsQlFO5eXl+Pbbb3Xf8Xq9KCwsRNeuXVFYWGg7Tkt4mZhj+/bt3oaGBt2zgwcPev/v//7P+/jjj2vP5s+f77355pu9x48f154dPXrUO3nyZG9+fr7tcF6v17tz505vZmam99///rf2rL6+3jtz5kzvfffdF7Q8horTp097f/nLX3p/9rOfac/s5Mlq2GgvJ6838LJqaGjwVlZWer1er3fXrl3ezMxM7+rVq32+E+1lFWg5We3P53o5tVSc4SaYeVq4cKE3NzfXO2/ePG9OTo7ut3O9nD7//HNvZmamd/369abfifZy8noDL6uTJ096p02b5s3Pz/c2NTUZfifayyoU48nWrVu9mZmZ3mXLlgUtThm20Y9BBgwYgLi4ON2zHj16oHfv3jhw4ID2bNu2bUhJSUFiYqL2LCkpCYMGDcKmTZtw6tQpW+EAYMOGDXC73brVrsfjQUZGBsrLy3HkyJGg5zeYtGrVCp07d8bJkye1Z3byZDVstJcTEHhZxcXFoVOnTn6/E+1lFWg5We3P53o5tVSc4SZYedq6dSs2bNiAmTNnKr/D5XSW2tpaNDU1KX+L9nICAi+rdevWoaqqCrfddhvcbjfq6urQ3Nzs851oL6tQjCfr1q2Dy+XCyJEjgxanDNvonyN4vV5UVVWhd+/e2rPGxka0adPGJ2ybNm1w+vRp7N+/H/3797ccDgD27t2Lnj17IiEhQRf2kksuAQDs27cPXbp0CWbWAqaurg719fWora3Fxx9/jE8//RRXX3219rudPFkNG43lBAS3rKwSjWUV6nJS9Wcup9DFGW6Cnaempib84x//wJgxY9CnTx/lN7mczvDUU0/h1KlTcLvdGDRoEGbMmKEz5YnGcgKCW1afffYZEhIScPToUeTn5+PgwYOIj4/HqFGjcOedd8Lj8diOM1II5Xhy+vRprF+/Hv3790e3bt2CEqcKFvTPEd5//30cPXoUt956q/bsggsuwM6dO9HU1IRWrVoBOCP879q1CwBw9OhRW+EAoLKyUqmlpWdi2Ejhueeew7vvvgsAcLvdSEtLw89//nPtdzt5sho2GssJCG5ZWSUayyrU5aTqz1xOoYsz3AQ7T++++y6OHDmCRx55xPCb53o5tW7dGldddRWGDBmCxMREfPXVV3jttdcwb948PProo7j44ottxxlJBLOsvvnmGzQ1NeGRRx5BRkYGpk+fji1btmDFihU4efIk7r//fttxRgqhHE/+97//4cSJE7jmmmt0z4NdTizonwMcOHAAf//739G/f3+MHj1ae3799dejsLAQTz/9NG666SZ4vV4sW7YMx44dAwA0NDTYCkf/L5sZANBW9GLYSGHixIlIT09HZWUlPvzwQzQ3N6OxsVH73U6erIaNxnICgltWVonGsgplORn1Zy6n0MUZboKZp+rqaixduhRTpkzRHSiVOdfLacCAARgwYID27+HDhyM9PR2//OUvsXjxYjz00EO244wkgllWpPW+7rrr8LOf/QwAcNVVV+H06dN49913ceutt6JHjx5RWVahHE/WrVuH1q1bY8SIEbrnwS4nttGPcY4dO4bf//73SEhIwLx58zSNPABcd911mDx5Mj744APcfffd+MUvfoFvv/0WkyZNAgDEx8fbCgecaYhiJyCoYVJDjSR69eqFyy+/HKNHj8aDDz6IU6dO4eGHH4b3/zukspMnq2GjsZyA4JaVVaKxrEJVTmb9mcspdHGGm2DmacmSJWjfvj0yMzNNv3mul5OKHj164Morr8Tnn3+u2exHYzkBwZ/3AOCHP/yhLizZnZeXl9uOM1IIVZs6deoUPv74Y6SmpurOPwYSpxEs6McwJ0+exPz583Hy5Ek89NBDmn9WkenTp+PFF1/EggUL8Oc//xlPPvmk1oBFH91WwyUlJWmafhF6pkpDpJGeno7du3fj4MGDAOzlyWrYWCgnILCyskoslFUwyslff+Zyark4w43TPH3zzTdYuXIlJkyYgMrKShw+fBiHDx9GY2MjmpqacPjwYZw4ccJWnJFMKOq+c+fOOH36NOrr64MWZyQQ6LwHAB07dtSFpR2jmpoa23FGKsFqU//9739RX1/vY7YTSJxGsKAfozQ0NODhhx/GwYMH8cADD+gO7cm0b98egwYN0g5lffbZZ+jcubOPf3wr4fr27YuDBw+itrZW9+7OnTu13yMdWjVTHuzkyWrYWCgnILCyskoslFWg5WSlP3M5tVyc4cZpno4ePYrm5mb885//xJ133qn97dy5EwcPHsSdd96JV155xVackUwo6v7QoUPweDzaTnYslBMQWFl973vfA+BrO15ZWQkAmsY6FsoqWG3q/fffR9u2bX0u3wokTiNY0I9Bmpqa8Oijj6K8vBzz5s3T3eDnjw8//BC7d+/GxIkT4XYbNw+jcOnp6WhubtYOrwBnDu6uWbMG/fr1i6gT9VVVVT7PTp8+jZKSEng8HvTq1QuAvTxZDRtN5QSEpqysEk1lFYpystqfz/VyCmcbDRXBzlPv3r2Rm5vr89e7d2906dIFubm5yMjIsBVnJBCKuj9+/LhPnPv27cPGjRuRmpqqzXvRVE5AaMqKbMxXr16ti3fVqlVo1aoVLrvsMttxhptQjifHjx9HWVkZrrzySp3pMxHscuLDuDHI888/j48//hjDhg3DiRMnsHbtWt3vo0aNAnDGj/Irr7yC1NRUdOjQATt37sSaNWswePBg3ZXVVsMBQL9+/ZCeno7Fixfj+PHj6N69O0pKSlBRUYF77rkn9Jm3wV//+lfU1tbi0ksvRVJSEqqqqvD+++/j66+/xh133IG2bdsCsJcnq2GjqZyA0JQVAM0rA2mCNm7cqP1/ZmYm2rVrF1VlFYpystqfz/VyClUbDSfBztN5552HtLQ0n++8+eabAKD77VwuJwB49NFH4fF40L9/f3Ts2BFfffUVVq5ciTZt2uAnP/mJFi6aygkITVldfPHFyMjIwOrVq9HU1IRLL70UW7ZswYYNGzB58mTN1CSayiqU48mHH36IpqYmpdmO0zjNcHnJ0JqJGX7zm99g69athr+/9dZbAIBvv/0Wf/vb3/DFF1/g1KlT6NatG0aPHo0bb7xRd+LbajiioaEBS5Yswfvvv4+amhr06dMH06ZNw+DBg4Of2QD44IMPsHr1anz55Zc4ceIE2rZti+9973vIzMzE8OHDdWHt5Mlq2GgpJyB0ZXXHHXegoqJC+c1nn31W8y0cLWUVinKy2p/txBluQlFOoWqj4aSl8vSb3/wG1dXV+Otf/xq0OFuSUJTTm2++iXXr1uHbb79FbW0tzjvvPPzgBz/A1KlT0aNHD0dxRgKhalOnT5/G8uXLsWbNGlRWVqJLly4YP348brjhBsdxhpNQ9r377rsPhw8fRlFRkc6hQiBxmsGCPsMwDMMwDMPEIGyjzzAMwzAMwzAxCAv6DMMwDMMwDBODsKDPMAzDMAzDMDEIC/oMwzAMwzAME4OwoM8wDMMwDMMwMQgL+gzDMAzDMAwTg7CgzzAMwzAMwzAxCAv6DMMwDMMwDBODtA53AhiGiU4mTJhgK3zXrl3x3HPPaTe9ijffxip1dXXYvHkzNm3ahO3bt6OiogJutxvdu3fHVVddhRtvvFG7St0fv/3tb1FWVgYAWLRoETp37qz7fcuWLcjNzTV8v1+/fnj88cf9fmfr1q3Izc2F1+tFRkaG4ZXru3btwquvvoodO3bg+PHjiI+Px4UXXoiMjAz86Ec/gsvl0oVfs2YN/ve//2Hfvn2oqqpCXV0dEhMT0b9/f9xwww0YOHCgMi1r167Fnj17UFlZiZqaGsTHx6Nv37649tprMWrUKJ/vmCG32Xnz5iE9PV3795NPPomSkhKf9xISEtCrVy+MHDkS119/vc9tlmvWrMFTTz2F0aNH49e//rXP+19//TXy8vJQWVmJcePGIScnBy6XS+sLf/jDH3DZZZdp4eXbkN1uN9q2bYvExET06dMHKSkpuOaaa9CuXTvLebcC5V9OT6Rxzz33YN++fdq/p06diltuuSWMKWKYyIUFfYZhHDF69GifZzt27MC3336Lvn37om/fvrrfEhMTWyppEcO6devwl7/8BQDQq1cvDBs2DKdOncKOHTvw0ksv4YMPPsAf//hHdOzY0TSeNWvWoKysDC6XC/4uM+/evTsGDBigfO6PxsZGLb1mbNiwAY8++iiam5tx8cUXY+DAgaiursa2bduwfft2fPbZZ7jvvvt077z99tvYt28fLrzwQgwcOBBxcXE4ePAgPvroI5SWluKuu+7Cddddp3vn448/xqpVq9CzZ09cdNFFaN++PY4ePYpt27Zhy5Yt+PTTT3H//ff7Ta9IfHw8rrrqKgBnFp8qBgwYoJVXc3MzKioqUF5ejp07d+LTTz/Fgw8+aHmBceDAAeTl5eHYsWO4/vrr8fOf/9zyu4MHD9baxqlTp/Ddd99h06ZNKC0txQsvvICf/vSnuPbaay3FFUsMGzYMffv2xbfffosdO3aEOzkME9GwoM8wjCNUmssnn3wS3377La688kpDDduvf/1r1NfX4/zzzw91EsNO69atMXbsWNxwww3o1auX9ryyshIPPfQQ9u7di2eeecZUWD1+/Dief/55pKam4uDBg6ioqDD95oABA5R1Y4Vly5bhm2++QUZGBlatWqUM09TUhL///e9obm7Gvffei2uuuUb77cCBA5g7dy7WrVuHMWPG4Ac/+IH2289//nP06tULCQkJuvg+/vhj/OEPf8Czzz6Lq666Cuedd572W0ZGBm688UaftvLNN9/gN7/5DT744AOMHDkSw4YNs5zHxMREv+UzZswYHwF69+7dmDdvHj799FOUlpZqiwUzvvrqK+Tl5aGqqgoTJkzAT3/6U8vpBICsrCwfzfrJkyfxxhtvYNmyZXjqqafQ1NSEsWPH2oo32pk2bRqAMwtgFvQZxhy20WcYpkXp2rUrevXqhdatY1/P8KMf/Qi/+MUvdEI+ACQlJeGuu+4CAJSWlqKxsdEwjmeeeQb19fVa+FCxf/9+/Pvf/0ZGRoZyR4D4+uuvUVVVhZ49e+qEfODMrgU92717t+63fv36+Qj5ADB8+HBcdtllaGhoQHl5ue633r17KxeEPXr0wPXXXw8A+Pzzz61kL2C+//3va2Y+27Zt8xt+//79mpB/44032hbyjWjXrh1uueUWzJo1CwDwz3/+E8eOHQtK3AzDxB6xP9MyDBNRGNnoT5gwAV27dsU///lPLF++HCUlJTh69Ci6du2Km266SdOwlpWVYdmyZdizZw/cbjeGDRuGO++8U2ka1NTUhJUrV6KkpARfffUVmpqa0LNnT/zoRz9CZmamj611S9KnTx8AZ8xlTpw4gaSkJJ8wn376KdatW4dp06ZZMr1xitfrxV//+le0a9cO2dnZ+Pjjjw3DxsXFWYqzQ4cOlr9P9WBn8efknUCh3YampibTcCTkHz9+HJMmTcKMGTOCnpZRo0bh3Xffxfbt27Fq1SpMmTLF8rurV6/GW2+9hYMHDyIhIQGDBw9Gdna2Yfht27bhww8/xLZt23DkyBE0NjaiS5cuuPLKK5GVlYX27dtrYTds2IAFCxbghz/8oeFO1V/+8hes/H/t3XlQ1PX/wPEnh8cueKAIhEcqKqOIIqhkYyoeoCUemOWR5Y0zOaiVOY7pMGphoRZjk9mIeOtMgUAqg3igeIGM1XKKCAoKsSSCcqgg/P5gdmPdBRdF88f39fiT/Vy77M683u/36/V6R0ezbNky7e+6pKSEI0eOkJCQQGFhIaamprRv3x5HR0cmTpxInz59jH5/Qoh/SaAvhHitfPvtt6hUKpydnbGzsyM5OZmgoCAAFAoFgYGBODo64urqSnp6OmfOnKGgoIBNmzbp5D4/evSI9evXo1KpaNOmDY6OjrRs2ZKMjAx27typLV41Nf1vFjYLCgqA2kDVUFD88OFDfvrpJ7p06YKPj4/R183Pz2fPnj08ePCAtm3b0q9fP1xdXRt8n8ePHyctLY0VK1Y8M0C3tbXljTfe4M6dO8TGxuql7sTGxmJpacmwYcOMet6//voLlUqFpaUljo6ORp1TWFhIVFQUAIMHDzbqnKaQmZkJoLdCU1d2djZfffUV9+/fZ/r06Xz88ccv7XlGjBhBamoqKpXK6EB/9+7dhIaGYm5ujrOzM0qlkqtXr5KUlKRXV6MREhJCdna2thD48ePHZGVlERoaypUrV9i8ebO2qNzd3R0rKysuXbrE/fv39QbgFRUVnDt3DqVSyfDhwwEoLy/n888/p6CgAGtra1xcXDAzM6OwsJC4uDjs7Owk0BfiOUmgL4R4bajVahQKBTt27NDOnqpUKtasWcO+ffuorKxkzZo1DBkyBKgNEFauXElqaipJSUk6OeG7du1CpVLxzjvv8Omnn2o7lJSXlxMYGEh8fDzR0dF6BaCvSmRkJFBbcGlolvzAgQOo1Wq++eYbo2fRobYg+um85e7du7N69Wrs7e31jr979y579+5lwIABBgusn2ZmZsby5ctZv349W7ZsITw8HHt7e0pKSkhJSaFr164sX7683gHDyZMnSUpKorKykvz8fDIzM7GwsGDlypU6M8N1paenExUVRXV1NUVFRaSmplJdXc1HH31E//79jfhUnt+TJ08oLCzk6NGjJCcnY21tjYeHh8Fj79y5w5o1a3jw4AEffvihNpf8ZdEE5rdv3zbq+PT0dMLCwrCwsODrr7/GwcEBqA2+N27cSEJCgsHzZsyYQd++fXW6/FRWVrJjxw6io6MJDw9n5syZQO3AdezYsfz666+cOXOGyZMn61zr3LlzVFRUMGHCBFq3bg3AxYsXKSgowN3dXW/wXVJSQnFxsXEfiBBCj+ToCyFeK4sWLdIpyBwwYAA9e/akqKgINzc3bZAPtW0PNYWIddsRFhcXc+LECaytrVm2bJlOgKJUKvHz88Pc3Jzjx4+/gnekLzExkZiYGMzNzQ0Gg5mZmURGRjJ69Gij2xwqlUp8fHzYvHkzBw8e5ODBg2zcuBFHR0du3rzJ2rVrKSsr0zvv559/5vHjx42qAejXrx8BAQHY2dlx48YN4uLiUKlUmJiY4OLigp2dXb3npqamcvr0aeLi4sjMzKRNmzb4+fnh6upa7zn5+fmcPn2a2NhYVCoV1dXVzJo1q1ErHY0RFBSEt7c33t7eTJkyhUWLFhEREcHIkSPZvHmzwVoDgGvXrvHgwQP69Onz0oN8+LeTVWlpqVHHR0VFUVNTg7e3tzbIh9qVMl9f33q7AQ0ePFivlWeLFi1YtGgRZmZmeqle48ePx9TUlOjoaL1rxcTEAOgUEJeUlAC1v/WnV57atWvHm2++adT7E0Lokxl9IcRrw9zc3OAMrZ2dHVlZWQwaNMjga1DbyUYjKSmJqqoq3NzcaNWqld45VlZW2Nvbc+vWLR49emTwmJclNzeXLVu2UFNTw7x58/TSJZ48ecK2bduwsLBg/vz5Rl/XwcFBJ3gDGDhwIP3792fNmjWkpKRw/Phxpk+frn394sWLXL58mRkzZtClSxej73X27FmCgoJwdHRk5cqVdOvWjaKiIsLCwggPDycpKYnAwECDKxF+fn74+flRUVHBnTt3CA0NJSAgAC8vL5YuXWrwfh4eHnh4eFBZWYlareb06dMcPnyYK1eu4O/vX+9KwPOq214T4N69e2RmZnL+/HksLS21Ae7TevbsSV5eHhkZGezevbvBvPem8KxWq0/TFBGPGDFC77Vu3brRo0cPsrKyDJ579+5dEhISuH37NuXl5VRXVwO1v9m8vDydY21sbHB1dSUxMZG0tDRtcffNmze5du0avXr10vmu9urVC4CwsDDat2/P4MGD6x1MCSEaRwJ9IcRro3379gYDKM0Sv6EOLJrc4LqdazQtKKOjow3OKtZVWlraYKCvaW/5NE9PT5ycnBq89tPu3r2Lv78/paWlTJkyhUmTJukdExkZSVZWFn5+fjorG8/LzMyMadOmkZKSwtWrV7WBfnl5OTt27MDe3p4PPvjA6Ovl5eXxww8/0K5dO9atW6f9/O3t7Vm6dClFRUVcuXKFmJgYbWccQxQKBb169WLVqlVUVlYSHR3NoEGDdDawelqLFi3o3Lkzc+bMoU2bNgQHB3PgwAF8fX2Nfn5jGGqvqUn5OnbsGJaWlgZn7Lt37868efNYv349oaGhKBSKRhXJNtb9+/cB4wufNYPhTp06GXzdxsbGYKAfHh7Onj17qKqqMvrZJkyYQGJiItHR0dpAX/NbfLod6MCBA5k8eTKRkZEEBgZiZmaGg4MDLi4ujBs3rsEVIiFEwyTQF0K8Np5VGGts4axmtrFnz57a7jb1eVbXlocPHxrcLdXZ2blRgf6DBw9Yu3YtarWasWPH1jtbn5CQgImJCadOndK7r6aN4qZNm2jRogXvv/8+bm5uz7y3Jje/bhvGGzduUFRUhI2NDevWrTN4n8TERFavXo2VlRVffvklUJtjXVVVhaurq8FdfYcPH86VK1dISUlpMNCva9SoUcTHxxMfH99goF+Xh4cHwcHBXL58uckDfUOUSiVz584lMTGRo0eP1pua4+LiwqpVqwgICGD//v0oFAqDA7qmoAnKGyoOflHp6ekEBwdjYWHB0qVLcXZ2xsrKSrta88knn+ispmm4ublhbW3N+fPnWbx4Mebm5sTGxqJQKAyuKCxcuJDx48cTHx/Pn3/+SVpamnb35S+++MLo74UQQpcE+kKIZsfa2hqozSV/0SDQ1taW33///YWuUVFRgb+/P7m5uQwbNoylS5c2uDtqTU1Ng73ar127BtT26TeGJofb0MqFWq2udxOue/fuce/ePZ0dZP/55x8AvZxtDc3fjc0bh39zzTW52sawtLTE1NRUO6v9KmjawZaVlVFSUlLviou7uzsrVqxg69at7Ny5E6VS+VJ2sI2LiwPQKUJviJWVFWq1msLCQoODA0Pfg8uXLwMwZ84cve/bo0eP6u3hb2ZmhpeXFwcOHCA2NhalUklpaSmenp71puV06dKFLl26MG3aNB4/fszRo0cJCQlh+/btEugL8Zwk0BdCNDuaor6EhAQWLFjwn27OVVlZycaNG8nIyMDV1ZWVK1c22L8/ICCg3tcWLFiAWq0mJCREO5gxxsWLFwF08qKdnZ3rHcCcPHmSoKAgxo0bh5+fn85rVlZWgP6GWBqav9cdHDyLppC6MXsFaDrvvMq0Dk1LVBMTk2fWdYwcOZKHDx/y448/sm3bNlq3bq1tJ9kUTp8+TVpaGq1atcLT09Ooc5ycnFCr1Zw/f17bJUcjNzeX7OxsvXM0AzZDaXMXLlxosE7A09OTw4cPEx0drQ3ujd3Ft2XLlvj4+BAREUFRURHFxcW0b9/eqHOFEP+SrjtCiGanY8eOjBs3DrVaTWBgoMFZx7y8PC5cuPBSn+PJkycEBgaiUqlwcnJi9erVjWqV2RgREREUFhbq/K2mpoaoqCgiIiIwMTExOpWmIe7u7gDa4t660tPTiYiIANCZgc3NzSUuLk5vB+CamhrOnTtHaGgoJiYmeu09w8LCDK4MZGRksG3bNoCXMlNuSHl5OSEhIQD0799fWzfSEC8vLxYuXEh1dTVbtmwhMTHxhZ+jrKyMQ4cOafeWWLJkidG1HJpWshERETpB/cOHD/nll18MBu2atK+YmBidHP2cnBx2797d4P06dOjA0KFDycrKIjk5me7duxvsh3/p0iW9XZGhtvtUcXExCoWi3hUkIUTDZEZfCNEsLVq0iIKCAi5evMjVq1fp0aMHnTp14tGjR+Tk5JCfn4+7u/tLTQk4duwYly5dAmrTU7Zv327wuPnz579w4W1kZCS7du3CwcEBW1tbKisruXnzJgUFBZiamrJ48WJtd5MX0atXL6ZOncqRI0fYvn07x44do2vXrhQVFXHt2jWqq6vx8vLCxcVFe05xcTHfffcdFhYWODg4YGVlRVlZGTk5OajVakxNTVmwYIFeEBgSEsK+fftwcHDAxsaGqqoq/v77b22QOnz48JeS/37ixAmSkpJ0nv/69evaTciWLFli9LUmT55MRUUFBw4cICAgAH9/f72WqfWlcf3222+cPHkSqA3G7969y40bN6iqqkKpVOLr62vU3gcaffv21f7vPvvsMwYMGIBSqSQ5OZkWLVowdOhQvV76Y8eOJTw8nISEBJYsWULv3r0pLS0lOTkZd3d3rl+/Xm/qF9QOLjS/gfHjxxs8Jjk5mcjISDp27EjPnj1RKpUUFRWRkpKibaX6sgbIQjR3EugLIZqlVq1a4e/vz9mzZzl16hTZ2dlcv36dtm3bYmNjg4eHh8GiwKZUdzZaE+wYMmvWrBcO9KdMmcIff/xBTk4Oubm5VFVV0aFDB0aNGoW3t3eT7iw6f/58+vbtS1RUFJmZmdy5cweFQkH//v3x9PRk5MiROsd369aN2bNnk5SURF5eHmlpaZiYmGBtbc24ceN49913DQ5CfH19UalUZGdnc+vWLaqqqmjXrh3u7u6MGTPG6N13G+vpTcdatmyJra0tY8aMwcfHR5u+ZKwZM2ZQUVFBWFgYGzZsYMOGDTg6OmpXOOpbHbh69SpQW4SuUCho06YNQ4YMYeDAgXh4eDxXC8r58+fTuXNnjh49SlJSEhYWFri4uDB37lz27t2rd3zbtm3ZunUru3fvJjk5mYSEBGxtbZk9ezZTp05l8eLFDd7PyckJc3NzTE1N9b4XGmPGjMHU1JSUlBSuX79OWVkZVlZWDB48mEmTJjFw4MBGv08hRC2TmsY24hVCCCH+H/L29sbGxobg4OD/+lGA2gLX4uJi9u/f3yStVF9HZ8+eZfPmzYwePZoVK1Y06bU1tSQzZ85k1qxZTXptIZoLmdEXQgjxP+P+/ft8//33AEycOJHevXv/J88RFxdHcXExXbt2bbZBflVVFaGhoQC89957TXbd/fv3U1hYSH5+fpNdU4jmSgJ9IYQQ/zPq7oswdOjQVx7oBwcHk5mZSWpqKoBe95vmID4+nsuXL5ORkUFOTg5vvfVWk6aOJSQkGOwQJITQJ6k7QgghxCuyYMECSkpK6N69Oz4+Prz99tv/9SM1uYMHD3Lo0CEsLS1xc3PD19fX6N17hRBNSwJ9IYQQQgghmiHpoy+EEEIIIUQzJIG+EEIIIYQQzZAE+kIIIYQQQjRDEugLIYQQQgjRDEmgL4QQQgghRDMkgb4QQgghhBDNkAT6QgghhBBCNEMS6AshhBBCCNEMSaAvhBBCCCFEM/R/GBEjmngjYGsAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.scatter();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This light curve of the object [K2-133](https://exoplanetarchive.ipac.caltech.edu/cgi-bin/DisplayOverview/nph-DisplayOverview?objname=K2-133+b&type=CONFIRMED_PLANET), which is known to host an exoplanet with a period of 3.0712 days. The light curve shows a lot of motion noise on K2's typical 6-hour motion timescale.\n", "\n", "Let's plot the folded version of it to see what the signal of the known transiting exoplanet looks like." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:23:11.861749Z", "iopub.status.busy": "2023-11-03T14:23:11.861514Z", "iopub.status.idle": "2023-11-03T14:23:11.883909Z", "shell.execute_reply": "2023-11-03T14:23:11.883603Z" } }, "outputs": [], "source": [ "lc.fold(period=3.0712).scatter();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can see the hint of an exoplanet transit close to the center, but the motion of the spacecraft has made it difficult to make out above the noise.\n", "\n", "We can use the `SFFCorrector` class to remove this motion. An in-depth look into how the algorithm works can be found [here](https://docs.lightkurve.org/tutorials/04-replicate-vanderburg-2014-k2sff.html). You can tune the algorithm using a number of optional keywords, including:\n", "\n", "* `degree` : *int*\n", " The degree of polynomials in the splines in time and arclength. \n", "* `niters` : *int*\n", " Number of iterations\n", "* `bins` : *int*\n", " Number of bins to be used to create the piece-wise interpolation of arclength vs flux correction.\n", "* `windows` : *int*\n", " Number of windows to subdivide the data. The SFF algorithm is run independently in each window.\n", "\n", "This [tutorial](https://docs.lightkurve.org/tutorials/04-replicate-vanderburg-2014-lightkurve.html) will teach you more about how to tune these parameters. For this problem, we will use the defaults, but increase the number of windows to 20." ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:23:11.885823Z", "iopub.status.busy": "2023-11-03T14:23:11.885682Z", "iopub.status.idle": "2023-11-03T14:23:12.299812Z", "shell.execute_reply": "2023-11-03T14:23:12.297217Z" } }, "outputs": [], "source": [ "corr_lc = lc.to_corrector(\"sff\").correct(windows=20)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Note: this is identical to the following command:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:23:12.305072Z", "iopub.status.busy": "2023-11-03T14:23:12.304677Z", "iopub.status.idle": "2023-11-03T14:23:12.859149Z", "shell.execute_reply": "2023-11-03T14:23:12.858258Z" } }, "outputs": [], "source": [ "from lightkurve.correctors import SFFCorrector\n", "corr_lc = SFFCorrector(lc).correct(windows=20)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now when we compare the two light curves we can see the clear signal from the exoplanet." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:23:12.871095Z", "iopub.status.busy": "2023-11-03T14:23:12.868311Z", "iopub.status.idle": "2023-11-03T14:23:12.953811Z", "shell.execute_reply": "2023-11-03T14:23:12.952856Z" } }, "outputs": [], "source": [ "ax = lc.fold(period=3.0712).scatter(color='red', alpha=0.5, label='Original light curve')\n", "ax = corr_lc.fold(period=3.0712).scatter(ax=ax, color='blue', alpha=0.5, label='Motion noise removed using SFF');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "VoilĂ ! Correcting motion systematics can vastly improve the signal-to-noise ratio in your lightcurve." ] } ], "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.9.13" } }, "nbformat": 4, "nbformat_minor": 4 }