From 737b45fe62e9c659b3d35c72a4e12ff518f24cf5 Mon Sep 17 00:00:00 2001 From: Sauravroy34 Date: Mon, 17 Mar 2025 17:20:23 +0530 Subject: [PATCH] Added notebook for TESS light curve --- Lightcurve/tess_light_curve.ipynb | 289 ++++++++++++++++++++++++++++++ 1 file changed, 289 insertions(+) create mode 100644 Lightcurve/tess_light_curve.ipynb diff --git a/Lightcurve/tess_light_curve.ipynb b/Lightcurve/tess_light_curve.ipynb new file mode 100644 index 0000000..ba4fcc2 --- /dev/null +++ b/Lightcurve/tess_light_curve.ipynb @@ -0,0 +1,289 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## TESS lightcurve using Stingray" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "##### TESS light curves are FITS format files that contain the output of the photometric extraction and subsequent systematics removal (cotrending) performed by the SPOC algorithms. The flux and respective uncertainties are provided at each cadence, with NaNs filling in any missing data values" + ] + }, + { + "cell_type": "code", + "execution_count": 168, + "metadata": {}, + "outputs": [], + "source": [ + "from astropy.io import fits\n", + "import numpy as np\n", + "from stingray import Lightcurve\n", + "import matplotlib.pyplot as plt" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "#### Dataset source : https://exo.mast.stsci.edu/exomast_planet.html?planet=HD100546b" + ] + }, + { + "cell_type": "code", + "execution_count": 153, + "metadata": {}, + "outputs": [], + "source": [ + "file = \"/home/saurav/Desktop/workspace/notebooks/Lightcurve/hlsp_tess-spoc_tess_phot_0000000280773136-s0037_tess_v1_lc.fits\"\n" + ] + }, + { + "cell_type": "code", + "execution_count": 182, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Filename: /home/saurav/Desktop/workspace/notebooks/Lightcurve/hlsp_tess-spoc_tess_phot_0000000280773136-s0037_tess_v1_lc.fits\n", + "No. Name Ver Type Cards Dimensions Format\n", + " 0 PRIMARY 1 PrimaryHDU 44 () \n", + " 1 LIGHTCURVE 1 BinTableHDU 167 3649R x 20C [D, E, J, E, E, E, E, E, E, J, D, E, D, E, D, E, D, E, E, E] \n", + " 2 APERTURE 1 ImageHDU 49 (13, 12) int32 \n", + "None\n", + "TIC 280773136\n" + ] + } + ], + "source": [ + "with fits.open(file) as hdul:\n", + " print(hdul.info())\n", + " img = hdul[2].data\n", + "Target = hdul[0].header[\"OBJECT\"]\n", + "print(Target)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Lets plot the ImageHDU (aperture pixels) " + ] + }, + { + "cell_type": "code", + "execution_count": 183, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAecAAAG1CAYAAADDSzyoAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIiJJREFUeJzt3XlwVfX9xvHnQswCyQ2LkICENYEQMGxWCVphCIugiAqCI4h7wDKytCoypYP+KiUdLa4UiUVEJaKGgtYqMVJZVKCsNToIWFBTSQiK5BLIAsn5/dHhTq9ZSEhIvt/D+zVzZ8y5J+d8bwZ48nzOuVeP4ziOAACAMZo09gIAAEAgwhkAAMMQzgAAGIZwBgDAMIQzAACGIZwBADAM4QwAgGEIZwAADEM4AwBgGMIZAADDEM64aHg8nho9NmzYoG+++UYej0dPPfVUheMcOXJEDz30kOLj49WsWTM1b95cAwYM0BNPPKHjx49Xu4b169frnnvuUffu3dWsWTN17dpV9913n3JzcyvsW15erhdffFF9+/ZVeHi4oqKiNGrUKH322WcV9i0pKdGcOXPUvn17hYWF6aqrrlJWVlbAPmdfU1WP+++/37/vXXfdVe2+33//vX/fP/zhDxo4cKDatGmj0NBQxcXFadasWTp69GiFdS5YsEA33nijoqKi5PF49Nhjj1X6c1qzZo1Gjhyp9u3bKyQkRB06dND48eP1xRdfVLr/iRMn9Mgjj6hLly4KCQnRZZddpvHjx+vUqVOV7g+YLqixFwA0lNdeey3g61dffVVZWVkVtvfs2VNFRUWVHmP79u0aPXq0CgsLNXnyZA0YMECStGPHDqWmpmrTpk368MMPq1zDnDlzdOzYMd16662Ki4vTwYMH9cILL+i9997Tnj17FB0d7d/34Ycf1qJFizR58mT96le/0vHjx7V06VINHjxYn376qa688kr/vnfddZcyMjI0a9YsxcXF6ZVXXtHo0aP18ccf65prrpEktWnTpsJrlaR169Zp5cqVGjFihH/b1KlTNWzYsID9HMfRtGnT1LlzZ1122WX+7Tt37lTfvn112223KSIiQnv37tVLL72kv//979qzZ4+aN2/u33fevHmKjo5Wv379lJmZWeXPKTs7Wy1bttTMmTN16aWXKi8vTy+//LKuvPJKbdmyRX369PHvW1BQoMGDB+s///mPUlJSFBsbq6NHj2rz5s0qKSlRs2bNqjwPYCwHuEhNnz7dqeqvwKFDhxxJzpNPPunf9tNPPzmXXXaZExUV5ezdu7fC9+Tl5Tm///3vqz3nxo0bnbKysgrbJDm//e1v/dtOnz7thIWFOePHjw/Y9+DBg44kZ8aMGf5t27Ztq7DWoqIip1u3bk5SUlK163Ecx0lOTna8Xq9TVFRU7X6bN292JDkLFiw45zEzMjIcSc4bb7wRsP3QoUOO4zjO0aNHHUnO/Pnzz3mss/Ly8pygoCBn6tSpAdsfeOABp0WLFs7BgwdrfCzAdIy1gRpaunSpvv/+ey1atEjx8fEVno+KitK8efOqPca1116rJk2aVNjWqlUr7d2717/t9OnTKioqUlRUVMC+bdu2VZMmTRQWFubflpGRoaZNmyolJcW/LTQ0VPfee6+2bNminJycKteTm5urjz/+WLfccotCQ0OrXXt6ero8Ho9uv/32aveTpM6dO0tShTH/2e3no23btmrWrFnAMY8fP67ly5crJSVFXbp0UWlpqUpKSs77HIApCGeght59912FhYVp/Pjx9XrcwsJCFRYW6tJLL/VvO3vd+JVXXtHKlSv13Xff6fPPP9ddd92lli1bBgTx7t271b17d3m93oDjnh1779mzp8pzr1q1SuXl5Zo0aVK1azx9+rTeeustDRo0qNKAdRxHP/zwg/Ly8rR582bNmDFDTZs21ZAhQ879A6jG8ePHdfToUWVnZ+u+++6Tz+dTcnKy//lPPvlExcXFio2N1fjx49WsWTOFhYXp6quvrvZ1A6bjmjNQQ3v37lX37t0VHBxcr8d95plnVFpaqokTJwZsf/311zVx4kRNnjzZv61r16769NNP1bVrV/+23NxctWvXrsJxz247fPhwledeuXKl2rVrp6FDh1a7xszMTP34449VhviRI0cC1tChQwelp6dXOmGojYEDB2rfvn2SpPDwcM2bN0/33nuv//kDBw5IkubOnatu3brp1VdfVUFBgR5//HENHTpUX375ZaU/G8B0hDNQQz6fTxEREfV6zE2bNunxxx/XhAkTKgRkRESEevXqpaSkJCUnJysvL0+pqam66aabtHnzZn/TLioqUkhISIVjnx1TV3Vz2/79+7Vz507Nnj27wqj959LT03XJJZdowoQJlT7fqlUrZWVlqbi4WLt379Zf//pXFRYWnvP1n8vy5cvl8/l08OBBLV++XEVFRSorK/Ov9+w5PB6P1q9fr/DwcElSv379lJSUpMWLF+uJJ56o8zqAhkY4AzXk9Xp14sSJejveV199pZtvvlm9e/fWX/7yl4Dnzpw5o2HDhmnIkCF6/vnn/duHDRumXr166cknn9Qf//hHSf8dgVd2nbW4uNj/fGVWrlwpSeccaRcWFuqdd97RyJEj1bp160r3CQ4O9t/dfcMNNyg5OVlXX3212rZtqxtuuKHa41cnKSnJ/9+33XabevbsKUn+t7idfW1jxozxB7P038bdpUuXSt92BtiAa85ADcXHx2v//v0qLS2t87FycnI0YsQIRUZG6v3336/QyDdt2qQvvvhCN954Y8D2uLg49ezZU59++ql/W7t27Sp9n/TZbe3bt690Denp6erRo4f/7WBVWbt2rU6dOnXOEP9fgwYNUrt27fy/ANSHli1baujQoQHHPPvafn7jnPTfG8h++umnejs/0JAIZ6CGxowZo6KiIq1evbpOx/nxxx81YsQIlZSUKDMzs9JrokeOHJEklZWVVXju9OnTOnPmjP/rvn37av/+/fL5fAH7bdu2zf/8z23btk1ff/11jQJ35cqVCg8Pr/CLwrkUFxeroKCgVt9zLkVFRQHHPPuLxf9+KMpZhw8fVps2ber1/EBDIZyBGpo2bZratWun3/zmN9q/f3+F5/Pz8895ffPkyZMaPXq0vv/+e73//vuKi4urdL/u3btL+u/d1P9r165d2rdvn/r16+ffNn78eJWVlSktLc2/raSkRMuXL9dVV12lmJiYCsdPT0+XpHO+Lero0aP66KOPdPPNN1f6YR4nT56s9FO4Vq9erZ9++klXXHFFtcevSn5+foVt33zzjdavXx9wzB49eqhPnz5655139MMPP/i3f/jhh8rJydHw4cPP6/xAY+OaM1BDLVu21Jo1azR69Gj17ds34BPCdu3apTfeeCPgGmllJk2apH/+85+65557tHfv3oD3NoeHh+umm26S9N9GOHz4cK1YsUI+n08jRoxQbm6unn/+eYWFhWnWrFn+77vqqqt06623au7cucrPz1dsbKxWrFihb775RsuWLauwhrKyMr355psaOHCgunXrVu1633zzTZ05c6bKhn3gwAENGzZMEydOVHx8vJo0aaIdO3bo9ddfV+fOnTVz5syA/V977TV9++23/kDftGmT/xeaO+64Q506dZIkXX755UpOTlbfvn3VsmVLHThwQMuWLdPp06eVmpoacMynn35aw4cP1zXXXKOpU6eqoKBAixYtUvfu3fXAAw9U+/oAYzX2p6AAjaW2nxB21uHDh53Zs2c73bt3d0JDQ51mzZo5AwYMcBYsWOAUFBRUe85OnTo5kip9dOrUKWDfU6dOOf/3f//nJCQkOGFhYU5kZKRzww03OLt3765w3KKiIuehhx5yoqOjnZCQEOcXv/iFs27dukrXsG7dOkeS89xzz1W7VsdxnIEDBzpt27Z1zpw5U+nzR48edVJSUpz4+HinefPmTnBwsBMXF+fMmjXLOXr0aIX9Bw8eXOXr//jjj/37zZ8/37niiiucli1bOkFBQU779u2d2267zfn8888rXUdWVpYzcOBAJzQ01GnVqpVzxx13OLm5ued8fYCpPI7jOI3wOwEAAKgC15wBADAM4QwAgGEIZwAADEM4AwBgGMIZAADDEM4AABiGcAYAwDCEMwAAhiGcAQAwDOEMAIBhCGcAAAxDOAMAYBjCGQAAwxDOAAAYhnAGAMAwhDMAAIYhnAEAMAzhDACAYQhnAAAMQzgDAGAYwhkAAMMQzgAAGIZwBgDAMIQzAACGIZwBADAM4QwAgGEIZwAADEM4AwBgGMIZAADDEM4AABiGcAYAwDBBjb0AAADOR3FxsUpLS+t8nODgYIWGhtbDiuoP4QwAsE5xcbG6dOmivLy8Oh8rOjpahw4dMiqgCWcAgHVKS0uVl5ennJwceb3e8z6Oz+dTTEyMSktL7Qxnj8dzIdeBn0lJSWnsJQCukpaW1thLuKg4jtMg54mIiFBERMR5f39DrbO2aM4AAGs5jlOngDU1nLlbGwAAw9CcAQDWcmtzJpwBANZyazgz1gYAwDA0ZwCAtdzanAlnAIC13BrOjLUBADAMzRkAYC23NmfCGQBgLcIZAADDuDWcueYMAIBhaM4AAGu5tTkTzgAAa7k1nBlrAwBgGJozAMBabm3OhDMAwFpuDWfG2gAAGIbmDACwllubM+EMALCWW8OZsTYAAIahOQMArOXW5kw4AwCsRTgDAGAYt4Yz15wBADAMzRkAYDVT229dEM4AAGsx1gYAAA2C5gwAsJZbmzPhXI2UlJTGXgIaQFpaWmMvAQ3gYvz7fDH82XZrODPWBgDgPKSmpsrj8WjWrFmSpGPHjunBBx9Ujx49FBYWpo4dO2rGjBkqKCio9bFpzgAAazVWc96+fbuWLl2qxMRE/7bDhw/r8OHDeuqpp5SQkKBvv/1W06ZN0+HDh5WRkVGr4xPOAABrNUY4FxYWatKkSXrppZf0xBNP+Lf37t1bq1ev9n/drVs3LViwQJMnT9aZM2cUFFTzyGWsDQCw1tlwrstDknw+X8CjpKSkynNOnz5d119/vYYNG3bO9RUUFMjr9dYqmCXCGQAAxcTEKDIy0v9YuHBhpfutWrVKu3btqvL5//XDDz/o97///XndjMhYGwBgrfoaa+fk5Mjr9fq3h4SEVNg3JydHM2fOVFZWlkJDQ6s9rs/n0/XXX6+EhAQ99thjtV4X4QwAsFZ9hbPX6w0I58rs3LlT+fn56t+/v39bWVmZNm3apBdeeEElJSVq2rSpTpw4oeuuu04RERFas2aNLrnkklqvi3AGAKAGkpOTlZ2dHbDt7rvvVnx8vObMmaOmTZvK5/Np5MiRCgkJ0bvvvnvOhl0VwhkAYK2GvFs7IiJCvXv3DtjWvHlztW7dWr1795bP59OIESN06tQpvf766/6byySpTZs2atq0aY3PRTgDAKxl0ieE7dq1S9u2bZMkxcbGBjx36NAhde7cucbHIpwBADhPGzZs8P/3kCFD6i3sCWcAgLVMas71iXAGAFjLreHMh5AAAGAYmjMAwFpubc6EMwDAWoQzAACGcWs4c80ZAADD0JwBANZya3MmnAEA1nJrODPWBgDAMDRnAIDVTG2/dUE4AwCsxVgbAAA0CJozAMBabm3OhDMAwFpuDWfG2gAAGIbmDACwllubM+EMALCWW8PZ49RwZR6P50KvBQBcKSUlpbGX0OCWLl16QY/v8/kUGRmprVu3Kjw8/LyPU1hYqIEDB6qgoEBer7ceV1g3XHMGAMAwjLUBANZy61ibcAYAWMut4cxYGwAAw9CcAQDWcmtzJpwBANZyazgz1gYAwDA0ZwCAtdzanAlnAIC13BrOjLUBADAMzRkAYC23NmfCGQBgLcIZAADDuDWcueYMAIBhaM4AAGu5tTkTzgAAa7k1nBlrAwBgGJozAMBqprbfuiCcAQDWYqwNAAAaBM0ZAGAttzZnwhkAYC23hjNjbQAADENzBgBYy63NmXAGAFiLcAYAwDBuDWeuOQMAYBiaMwDAWm5tzoQzAMBabg1nxtoAABimxs05JSXlQq6jWmlpaY12brhfY/7m7PF4Gu3caDiN9W9YY/673VDc2pwZawMArOXWcGasDQCAYWjOAABrubU5E84AAGu5NZwZawMAYBiaMwDAWm5tzoQzAMBahDMAAIZxazhzzRkAAMPQnAEA1nJrcyacAQDWcms4M9YGAMAwNGcAgNVMbb91QXMGAFjr7Fi7Lo/zlZqaKo/Ho1mzZvm3FRcXa/r06WrdurXCw8M1btw4HTlypNbHJpwBAKil7du3a+nSpUpMTAzYPnv2bP3tb3/T22+/rY0bN+rw4cO65ZZban18whkAYK3GaM6FhYWaNGmSXnrpJbVs2dK/vaCgQMuWLdOiRYs0dOhQDRgwQMuXL9dnn32mrVu31uochDMAwFr1Fc4+ny/gUVJSUuU5p0+fruuvv17Dhg0L2L5z506dPn06YHt8fLw6duyoLVu21Op1Ec4AgIteTEyMIiMj/Y+FCxdWut+qVau0a9euSp/Py8tTcHCwWrRoEbA9KipKeXl5tVoPd2sDAKxVX+9zzsnJkdfr9W8PCQmpsG9OTo5mzpyprKwshYaGnvc5a4LmDACwVn2Ntb1eb8CjsnDeuXOn8vPz1b9/fwUFBSkoKEgbN27Uc889p6CgIEVFRam0tFTHjx8P+L4jR44oOjq6Vq+L5gwAsFZDfkJYcnKysrOzA7bdfffdio+P15w5cxQTE6NLLrlE69ev17hx4yRJ+/bt03fffaekpKRarYtwBgCgBiIiItS7d++Abc2bN1fr1q392++99179+te/VqtWreT1evXggw8qKSlJAwcOrNW5CGcAgLVM+2ztp59+Wk2aNNG4ceNUUlKikSNH6s9//nOtj0M4AwCs1djhvGHDhoCvQ0NDtXjxYi1evLhOx+WGMAAADENzBgBYq7Gb84VCOAMArOXWcGasDQCAYWjOAABrubU51zic09LSLuQ6jJSSktJo574Yf96m/iW5kC7G1+zxeBp7CXARt4YzY20AAAzDWBsAYC23NmfCGQBgLcIZAADDuDWcueYMAIBhaM4AAGu5tTkTzgAAa7k1nBlrAwBgGJozAMBabm3OhDMAwGqmBmxdMNYGAMAwNGcAgLUYawMAYBi3hjNjbQAADENzBgBYy63NmXAGAFiLcAYAwDBuDWeuOQMAYBiaMwDAWm5tzoQzAMBabg1nxtoAABiG5gwAsJZbmzPhDACwllvDmbE2AACGoTkDAKzl1uZMOAMArOXWcGasDQCAYWjOAABrubU5E84AAGsRzgAAGIZwbkQpKSmNct60tLRGOW9jMvUPKtyjsf6MeTyeRjkvcD6sCGcAACpDcwYAwDBuDWfeSgUAgGFozgAAa7m1ORPOAACrmRqwdcFYGwAAw9CcAQDWYqwNAIBh3BrOjLUBADAMzRkAYC23NmfCGQBgLcIZAADDuDWcueYMAIBhaM4AAGu5tTkTzgAAa7k1nBlrAwBgGJozAMBabm3OhDMAwFpuDWfG2gAAGIbmDACwllubM+EMALCWW8OZsTYAAIahOQMArOXW5kw4AwCsRTgDAGAYt4Yz15wBADAM4QwAsNbZ5lyXR20sWbJEiYmJ8nq98nq9SkpK0gcffOB/Pi8vT3fccYeio6PVvHlz9e/fX6tXr6716yKcAQDWauhw7tChg1JTU7Vz507t2LFDQ4cO1dixY/Xll19KkqZMmaJ9+/bp3XffVXZ2tm655RZNmDBBu3fvrtV5CGcAAGpozJgxGj16tOLi4tS9e3ctWLBA4eHh2rp1qyTps88+04MPPqgrr7xSXbt21bx589SiRQvt3LmzVuchnAEA1qqv5uzz+QIeJSUl5zx3WVmZVq1apZMnTyopKUmSNGjQIL355ps6duyYysvLtWrVKhUXF2vIkCG1el01vls7JSWlVgdG3Zh6ByEAmKY+/r2MiYkJ+Hr+/Pl67LHHKt03OztbSUlJKi4uVnh4uNasWaOEhARJ0ltvvaWJEyeqdevWCgoKUrNmzbRmzRrFxsbWaj28lQoAcNHLycmR1+v1fx0SElLlvj169NCePXtUUFCgjIwM3Xnnndq4caMSEhL0u9/9TsePH9dHH32kSy+9VGvXrtWECRO0efNmXX755TVeD+EMALBWfb3P+ezd1zURHBzsb8IDBgzQ9u3b9eyzz+qRRx7RCy+8oC+++EK9evWSJPXp00ebN2/W4sWL9eKLL9Z4XYQzAMBaJnwISXl5uUpKSnTq1ClJUpMmgbdzNW3aVOXl5bU6JuEMAEANzZ07V6NGjVLHjh114sQJpaena8OGDcrMzFR8fLxiY2M1depUPfXUU2rdurXWrl2rrKwsvffee7U6D+EMALBWQzfn/Px8TZkyRbm5uYqMjFRiYqIyMzM1fPhwSdL777+vRx99VGPGjFFhYaFiY2O1YsUKjR49ulbnIZwBANZq6HBetmxZtc/HxcWd1yeC/RzhDACwlgnXnC8EPoQEAADD0JwBANZya3MmnAEA1nJrODPWBgDAMDRnAIC13NqcCWcAgLXcGs6MtQEAMAzNGQBgLbc2Z8IZAGAtt4YzY20AAAxDcwYAWMutzZlwBgBYi3AGAMAwbg1nrjkDAGAYmjMAwFpubc6EMwDAWm4NZ8baAAAYhuYMALCWW5sz4QwAsJqpAVsXjLUBADAMzRkAYC3G2gAAGOaiD+e0tLQLuQ78jMfjaZTzmvoHFe7RWH+2AZvQnAEA1rromzMAAKYhnAEAMIxbw5m3UgEAYBiaMwDAWm5tzoQzAMBabg1nxtoAABiG5gwAsJZbmzPhDACwllvDmbE2AACGoTkDAKzl1uZMOAMArOXWcGasDQCAYWjOAABrubU5E84AAGsRzgAAGMat4cw1ZwAADENzBgBYy63NmXAGAFjLreHMWBsAAMPQnAEA1nJrcyacAQBWMzVg64KxNgAAhqE5AwCsxVgbAADDuDWcGWsDAGAYmjMAwFpubc6EMwDAWoQzAACGcWs4c80ZAADDWNGcU1JSGuW8aWlpjXJe4ELzeDyNvQSgXri1OVsRzgAAVMat4cxYGwAAw9CcAQDWcmtzJpwBANZyazgz1gYAwDA0ZwCAtdzanAlnAIC13BrOjLUBADAM4QwAsNbZ5lyXR20sWbJEiYmJ8nq98nq9SkpK0gcffBCwz5YtWzR06FA1b95cXq9X1157rYqKimp1HsbaAABrNfRYu0OHDkpNTVVcXJwcx9GKFSs0duxY7d69W7169dKWLVt03XXXae7cuXr++ecVFBSkf/3rX2rSpHZdmHAGAFirocN5zJgxAV8vWLBAS5Ys0datW9WrVy/Nnj1bM2bM0KOPPurfp0ePHrVeF2NtAMBFz+fzBTxKSkrO+T1lZWVatWqVTp48qaSkJOXn52vbtm1q27atBg0apKioKA0ePFiffPJJrddDOAMArFVf15xjYmIUGRnpfyxcuLDKc2ZnZys8PFwhISGaNm2a1qxZo4SEBB08eFCS9Nhjj+n+++/XunXr1L9/fyUnJ+vAgQO1el2MtQEA1qqvsXZOTo68Xq9/e0hISJXf06NHD+3Zs0cFBQXKyMjQnXfeqY0bN6q8vFySNHXqVN19992SpH79+mn9+vV6+eWXqw38nyOcAQAXvbN3X9dEcHCwYmNjJUkDBgzQ9u3b9eyzz/qvMyckJATs37NnT3333Xe1Wg9jbQCAtRr6rVSVKS8vV0lJiTp37qz27dtr3759Ac/v379fnTp1qtUxac4AAGs19N3ac+fO1ahRo9SxY0edOHFC6enp2rBhgzIzM+XxePTwww9r/vz56tOnj/r27asVK1boq6++UkZGRq3OQzgDAFBD+fn5mjJlinJzcxUZGanExERlZmZq+PDhkqRZs2apuLhYs2fP1rFjx9SnTx9lZWWpW7dutToP4QwAsFpDfj72smXLzrnPo48+GvA+5/NBOAMArOXW//EF4QwAsJZbw5m7tQEAMAzNGQBgLbc2Z8IZAGAtt4YzY20AAAxDcwYAWMutzZlwBgBYy63hzFgbAADD0JwBANZya3OucTinpKRcyHXAEB6Pp7GXAAA15tZwZqwNAIBhGGsDAKzl1uZMOAMArEU4AwBgGLeGM9ecAQAwDM0ZAGAttzZnwhkAYC23hjNjbQAADENzBgBYy63NmXAGAFjLreHMWBsAAMPQnAEA1nJrcyacAQDWcms4M9YGAMAwNGcAgNVMbb91QTgDAKzl1rE24QwAsJZbw5lrzgAAGIbmDACwllubM+EMALCWW8OZsTYAAIahOQMArOXW5kw4AwCs5dZwZqwNAIBhaM4AAGu5tTkTzgAAa7k1nBlrAwBgGJozAMBabm3OhDMAwFqEMwAAhnFrOHucGq7M4/Fc6LXAACkpKY29hAaXlpbW2EtAA+DPdsO60KHn8/kUGRmp+Ph4NW3a9LyPU1ZWpq+++koFBQXyer31uMK6oTkDAKzl1uZMOAMArOXWcOatVAAAGIbmDACwllubM+EMALCWW8OZsTYAAIahOQMArOXW5kw4AwCs5dZwZqwNAIBhaM4AAKuZ2n7rgnAGAFjLrWNtwhkAYC23hjPXnAEAMAzNGQBgLbc2Z8IZAGAtt4YzY20AAAxDcwYAWMutzZlwBgBYy63hzFgbAADD0JwBANZya3MmnAEA1nJrODPWBgCghpYsWaLExER5vV55vV4lJSXpgw8+qLCf4zgaNWqUPB6P1q5dW+vz0JwBANZq6ObcoUMHpaamKi4uTo7jaMWKFRo7dqx2796tXr16+fd75pln5PF4zntdhDMAwFoNHc5jxowJ+HrBggVasmSJtm7d6g/nPXv26E9/+pN27Nihdu3ande6CGcAgLXqK5x9Pl/A9pCQEIWEhFT7vWVlZXr77bd18uRJJSUlSZJOnTql22+/XYsXL1Z0dPR5r4trzgCAi15MTIwiIyP9j4ULF1a5b3Z2tsLDwxUSEqJp06ZpzZo1SkhIkCTNnj1bgwYN0tixY+u0HpozAMBa9dWcc3Jy5PV6/dura809evTQnj17VFBQoIyMDN15553auHGjvv76a/3jH//Q7t27z3s9ZxHOAABr1Vc4n737uiaCg4MVGxsrSRowYIC2b9+uZ599VmFhYfr3v/+tFi1aBOw/btw4/fKXv9SGDRtqvC7CGQCAOigvL1dJSYkef/xx3XfffQHPXX755Xr66acr3Eh2LoSzoVJSUhp7CRcNfta40NLS0hp7Ca7V0Hdrz507V6NGjVLHjh114sQJpaena8OGDcrMzFR0dHSlN4F17NhRXbp0qdV5CGcAgLUaOpzz8/M1ZcoU5ebmKjIyUomJicrMzNTw4cPPew2VIZwBAKihZcuW1Wr/8/3FgXAGAFjLrZ+tTTgDAKzl1nDmQ0gAADAMzRkAYDVT229dEM4AAGvVNZhNDXbCGQBgLbeGM9ecAQAwDM0ZAGAttzZnwhkAYC23hjNjbQAADENzBgBYy63NmXAGAFjLreHMWBsAAMPQnAEA1nJrcyacAQDWcms4M9YGAMAwNGcAgLXc2pwJZwCAtQhnAAAM49Zw5pozAACGoTkDAKzl1uZMOAMArOXWcGasDQCAYWjOAABrubU5E84AAGu5NZwZawMAYBiaMwDAWm5tzoQzAMBabg1nj2PqygAAqILP51NkZKQkyePxnPdxzkZgQUGBvF5vvaytPtCcAQBWc2PH5IYwAIB1goODFR0dXS/Hio6OVnBwcL0cq74w1gYAWKm4uFilpaV1Pk5wcLBCQ0PrYUX1h3AGAMAwjLUBADAM4QwAgGEIZwAADEM4AwBgGMIZAADDEM4AABiGcAYAwDD/D2dAnv4JkVBAAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "fig, ax = plt.subplots()\n", + "cax = ax.imshow(img, cmap=plt.cm.grey, origin=\"lower\")\n", + "cbar = fig.colorbar(cax)\n", + "fig.suptitle(f\"{Target}\")\n", + "plt.axis(\"off\")\n", + "plt.show()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Since the data has NAN in its better to clean the values " + ] + }, + { + "cell_type": "code", + "execution_count": 184, + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "True\n" + ] + } + ], + "source": [ + "with fits.open(file) as hdulist:\n", + " tess_bjds = hdulist[1].data['TIME']\n", + " sap_fluxes = hdulist[1].data['SAP_FLUX']\n", + " \n", + " \n", + " sap_fluxes_error = hdulist[1].data[\"SAP_FLUX_ERR\"]\n", + " pdcsap_fluxes = hdulist[1].data['PDCSAP_FLUX']\n", + " pdcsap_fluxes_error = hdulist[1].data[\"PDCSAP_FLUX_ERR\"]\n", + "\n", + "diff = np.abs(tess_bjds[0] - tess_bjds[1])\n", + "nan_mask = np.isnan(tess_bjds)\n", + "check = True in nan_mask\n", + "print(check)\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Over here time values has NAN too, to maintain sequnces we cant just replace it with some values is better we add diff to previous not NAN value as shown below" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [ + "tess_bjds_clean = tess_bjds.copy()\n", + "for i in range(0, len(tess_bjds_clean)):\n", + " if nan_mask[i]: \n", + " tess_bjds_clean[i] = tess_bjds_clean[i - 1] + diff " + ] + }, + { + "cell_type": "code", + "execution_count": 146, + "metadata": {}, + "outputs": [], + "source": [ + "sap_fluxes_corrected = np.nan_to_num(sap_fluxes,nan = np.nanmedian(sap_fluxes))\n", + "sap_fluxes_corrected[0:10]\n", + "t = np.linspace(0,100,3649)\n", + "sap_fluxes_error_corrected = np.nan_to_num(sap_fluxes_error , nan = np.nanmedian(sap_fluxes_error))\n" + ] + }, + { + "cell_type": "code", + "execution_count": 149, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 149, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlUAAAGwCAYAAACAZ5AeAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAATedJREFUeJzt3XlcVPX+P/DXgMwAwrAvIohbgrhgotCUmhoxGter5S2zMnNNw1Lpqtd7Da1bYfZtsdzq9jXs981SKyt3CRVLcUNRRCU1FBIG3JhBZJ35/P7wcmQEFPDAsLyej8c8Ys55n895nxM6L882CiGEABERERHdFytLN0BERETUEjBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhk0MbSDbQmJpMJ2dnZcHR0hEKhsHQ7REREVAtCCBQUFMDHxwdWVjUfj2KoakTZ2dnw8/OzdBtERERUD1lZWfD19a1xPkNVI3J0dARw63+KWq22cDdERERUGwaDAX5+ftLneE0YqhpRxSk/tVrNUEVERNTM3OvSHV6oTkRERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkwFBFREREJAOGKiIiIiIZWDRUrVy5Er1795a+tkWj0WDbtm3S/MGDB0OhUJi9pk2bZjZGZmYmIiMjYW9vD09PT8yZMwfl5eVmNXv27EHfvn2hUqnQtWtXxMXFVell+fLl6NixI2xtbREWFoZDhw6ZzS8uLkZUVBTc3Nzg4OCA0aNHIzc3V76dQURERM2aRUOVr68vFi9ejOTkZBw5cgRDhw7FyJEjkZaWJtVMmTIFOTk50mvJkiXSPKPRiMjISJSWlmL//v1Ys2YN4uLiEBMTI9VkZGQgMjISQ4YMQUpKCmbNmoXJkydjx44dUs26desQHR2NhQsX4ujRowgODoZWq0VeXp5UM3v2bGzatAkbNmxAYmIisrOz8dRTTzXwHiIiIqJmQzQxLi4u4osvvhBCCPHoo4+KmTNn1li7detWYWVlJXQ6nTRt5cqVQq1Wi5KSEiGEEHPnzhU9evQwW27MmDFCq9VK70NDQ0VUVJT03mg0Ch8fHxEbGyuEECI/P1/Y2NiIDRs2SDWnT58WAERSUlKN/RUXFwu9Xi+9srKyBACh1+trsSeIWo+i0nJRUma0dBtERNXS6/W1+vxuMtdUGY1GfPvttygsLIRGo5Gmf/3113B3d0fPnj0xf/583Lx5U5qXlJSEXr16wcvLS5qm1WphMBiko11JSUkIDw83W5dWq0VSUhIAoLS0FMnJyWY1VlZWCA8Pl2qSk5NRVlZmVhMYGIgOHTpINdWJjY2Fk5OT9PLz86vPriFq0W6UlCPwje3otmAbEn+/bOl2iIjqrY2lG0hNTYVGo0FxcTEcHBywceNGBAUFAQCee+45+Pv7w8fHBydOnMC8efOQnp6OH374AQCg0+nMAhUA6b1Op7trjcFgQFFREa5fvw6j0VhtzZkzZ6QxlEolnJ2dq9RUrKc68+fPR3R0tPTeYDA0WLA6eUmP9Uey8PKjXdDe2a5B1kHUEDYcyZJ+Hr/6EHxd7GBtpcCGaRp4OtpasDMiorqxeKgKCAhASkoK9Ho9vvvuO4wfPx6JiYkICgrC1KlTpbpevXqhXbt2eOyxx3D+/Hl06dLFgl3XjkqlgkqlapR1PbVyP0rLTfgq6SLCOrli3cuaey9E1MhKyo3YfSYP5SYBb7Ut/N3aIm7/BbOaP68XAQBC30nA6beGwU5pLc37+4bjOJ1jwJqJoXB3aJw/W3IrM5qgAPD4R3uRnV+Ehzq7Yc3EUAC39s+N4nK4NdNtI2rtLB6qlEolunbtCgAICQnB4cOHsXTpUnz22WdVasPCwgAA586dQ5cuXeDt7V3lLr2KO/K8vb2l/955l15ubi7UajXs7OxgbW0Na2vramsqj1FaWor8/Hyzo1WVayyttNwk/Xww4xqivj6KT8Y+CGsrhQW7osrO5d3An9dv4tFuHlAomt7/l6xrN+HSVgkHVdW/FowmgSMXrmHBjycBANtmDkQbayss/eUsDmZcRbCfM+ZqA6BQKCCEqHb78m+Wos9b8XXqqXvMdlxYHAkAmPvdcXyX/CcAoN/bv1QJXM3B1wcv4l8bT5pNS/z9MlL/1MPbyRb93/kFAKDt4YXPxvWzRItEdB8sHqruZDKZUFJSUu28lJQUAEC7du0AABqNBu+88w7y8vLg6ekJAIiPj4darZZOIWo0GmzdutVsnPj4eOm6LaVSiZCQECQkJGDUqFFSDwkJCZgxYwaAW2HPxsYGCQkJGD16NAAgPT0dmZmZZtd/WVKwrxOO/6mX3m9JzcGW1BwAQPrbw6BqU7sPH5NJYNOJbCitreDtZIvevs6yBzMhBM7oCuDhqKrX0YbiMiP2n7+Cvh1c4GyvBABcuVECK4UCrm2VsvZaXxuP/QnXtio82s0DHf+xpcr8iqBQFyaTwIlLenTxaAtHW5v77rHMaMLXBy6ig5s9JsYdkab/Ev0ouno6AAB2pOlwLu8G3t+RbrZs139tw/+bFIqPfvkdALD//FWs3HPerKZiG3X6Yvzl019x5UbpPXuyV1ojakhXs/V1/McWXFgcifVH/jSr7R6zHf38XfDd9IcB3DoF/u7W09h//ircHVRY+mwfPNzFrckE2Op+DyqMWPab2fsdabkY9vFebJ81qKHbIiIZKYQQwlIrnz9/PoYPH44OHTqgoKAAa9euxXvvvYcdO3agc+fOWLt2LZ544gm4ubnhxIkTmD17Nnx9fZGYmAjg1sXtffr0gY+PD5YsWQKdTodx48Zh8uTJePfddwHceqRCz549ERUVhYkTJ2LXrl147bXXsGXLFmi1WgC3Hqkwfvx4fPbZZwgNDcXHH3+M9evX48yZM9K1VtOnT8fWrVsRFxcHtVqNV199FQCwf//+Wm+vwWCAk5MT9Ho91Gq1nLsSI5f9huN/6jHwAXf8evZKlfn9O94KIP95seZ//Qoh0Gn+1irTHVVt4GjbBu/9rTce7uIOnaEYQgj4utjDUFyGD3f+ju+T/0RBye3ng/3P08H4W4gvACDXUIyfUi7hXN4N7DyVi/ybZVLdkw+2x4fPBJt98M377gR+Pp6NT8c+iPCg29e6FZcZcS7vBv7y6e0PoIzYJxD4xnaUVDpSVzmwlJQb8cyqJNgr22DtlLAqH7BCCJy8ZICyjRUCvB1r3DeVXSssRd9/3zri8lTf9hAC6OLRFtMe7QJDcTkO/nEV078+CgD43/H9MGnNkSpjJLz+KLp43AouHf+xBVYKYOqgLnh5UGe4VBMML14txKPv75HeW1spsHZyGMI6u9Wq58oyrhRi49E/8cmuczXWLH22D2Z+m1LnsSuzUgDbZg6C9uO9tapPiXlcCsnA3UNIZbPDu2Fm+AM11qcuisDsdSk4lHENPs52NQaVS/lFKC4zSv9f5GAyCXT+Z9U/U3VhZ2ONhzq7Ynf67Yv4v5+uwd7fr+DkJT1c2yrx/tPB99tqi5ajL4K1QgFPdd2u0TOaBMpNplr/o5Rartp+fls0VE2aNAkJCQnIycmBk5MTevfujXnz5uHxxx9HVlYWXnjhBZw8eRKFhYXw8/PDk08+iQULFpht0MWLFzF9+nTs2bMHbdu2xfjx47F48WK0aXP7INyePXswe/ZsnDp1Cr6+vnjjjTfw0ksvmfWybNkyvP/++9DpdOjTpw8++eQT6XQjcOvhn6+//jq++eYblJSUQKvVYsWKFXU6/dcYoWr1S/1QUFyOtQczcTDjWpW6bl4O+OGVR8xO8RSXGaH9eC8uXr1Zpf5+nPn3MER8tBeZ1+497um3hmFrag5SL+mrXGNTVxXbOOvbFPxy2vy07ufjQhDRwxs3S8sRFLPDbF6wnzN+fOVhZF67iWc+S0J7ZzvEjOiBbak5ULWxQr+OrhjQ1b3eH5JRQ7pg+e7bR3MmD+iEL37LqNdYd2rvbIdL+UXwdbHD2NAOmPhIJ+nUmL6oDG2sFGiraoO3N5+q1zoDvByRnluAp/q2xw9HL5nNG9DVHU/0aodFP6fB19UOf1wurHGcqYM6459PdAcAXL1x64h0TdcPXblRgn5v/1Jl+oZpGhzKuFbl6FltbX51AExCYEtqDn46lo03/hIEfVEZ/rkxFQCw9bWBCPKR58/n3zfcPmVZ4fy7T0hHf40mgS6Vfp9eG9oVA7t54OlVNd9VXJP6HP0sKjXik11nEezrhGE929V5+aYu11CMsHcTzKadfWc4bKxv3/heXGZEuUlUOe1dORBX/r1taiI/+RW/5xagzCjwosYfb43saemWWqRmEapam4YKVZU/KFe/1A9DA28f3UnXFeBY5nX844fUKst5OKoQ1E5d7W3s8bMH4fGPaneEobnZOXsQIiywbRcWR9b66EtNvNW20BmKZenHyc4GAd6OUNvaIP9mKTZM0yAoZgfKTSYIAbRztkV7ZzvETQiFrc3tf6lXPqI5pp8f3vtbb7Nx95+7gulfH4W+qMxs+rl3hqONdd2e4nLn0dMHPB0QH/0oAGDmt8fwU0p2lWUOzH8MJ/7Mx9T/l1yndVXmYn//p1cB4Hqlo7KOtm1w+F/hZvsSAFL/1ONkth6RvdtB/d/TutcKS/FzyiWkZOXjx2q28V7c2iohACTOGQwHVRuYBLD9pA7O9jY48aceK3afMzuyDABTBnbCP5/oftfTpdtP6lBuMkEBBcqMJuiLyuDtZIsdaTqkZOaju48ay5/re9fe0rL1iD+Vixce8q/V6X8hBErKb63r6VVJKC03YdOrA7Dgx1Q82MEF0x69ddPSn9dvIi3bAP3NMgzv5Q1HW5sa/7xdWByJm6XlWH84C4s2nQIArHqhLwxF5fggPh0vajrC1sYa/958ymy5+vwOy6mguAxn825AbdsGBzOuYU/6ZcSfqvrNHu4OKvTwUSOyVzs805+P8pEDQ1UT1FChKvKTX5GWbQBQNVRV2JOeh5e+PHzPsdZMDMWj3TwA3DpNdCzzOv7S2wdndAb8eb0I6boC9PFzxk8pl/BjSjYUCsDPxR575w6RxqjuL7L+HV2w5G/BMAmBKwUl6NPBGao21hi9cj+SL16vUh/o7QhvJ1vs+e8pD6W1FUL8XbDorz1gbQUkX7yOed/fCooxfwnCxAGdcDa3oF5B8IWHOuD/DmTWebnkBeH4MSUb7ZxscSjjGlzbKqG2bYMPdv6OsM5uGNTNHTE/3Xpe2rxhgZg+uAtKy03otmCb2ThPh/giZkQQNp/IwdbUHLi1VWL/+avIK7h9beG4h/zx5l97wOq/RzhOXtKbnQbt4tEWbg4qHKrm6GR19s4Zgg5u9nXe5roQQiA9twAKKNDNy6He1zadzjHg3a2nMXVQZwx8wMNs3gtfHMRv526f7q58lCnhdK506nX7rIEY9vGv9dyS+7f4qV54NrRDvZat+Cv6ckEJ3BxU6PLPrfB0VCGvoAQ7ZtX+9Gp9VD76VZd/EHg4qhA/e5DZ6VwhBCbEHZb+TAOA2rYNHG1toO3hjeG9vNHJvW2VoHWv9Y7s44OiUiN2VhMuauJsb2N2GUIFdwflPa/9C+3kipul5VgyOvi+j2jWdFNHZfk3S7Hx2CW8uenUXetq8s8nAjF1UNO/W76pY6hqghoqVK0/nIW5358AUHOoqrD/3BV8fSgTB/+4his3bn9ov/NkTzwf5i9LPyaTwOp9GXh7y2l0dLPH5y/2Qzevmq9XGvHpb7hWWIqSchPWvfwQnOxs7ut2+ez8Ijzy3i60U9vCyV6JbTMHAgCGfrCnyqmpjNgnoFAocCm/CI8s3mU2r2d7Nbp6OOCZ/n54uIs7tpzIwcZjl2ASAjMfewDBfs737KXMaIJOXww/19sBxmgSKCguw5/Xi+DtZCvrowGMJoHE3/Mw97sTCPRWo8xogtrOBjbWCpzOKYCnowqrX+qPttXc4ddclZQbUVxmgpNd1aNLlT+0vtyXIX0wnXtnOKytFCgzCnywMx1XC0vh5qDEmH5+MMn8V6Kdsk2DPjuu4ho/Z3sb2FhbYY42AHO/O1Hr5TWd3ZD0x9Ua5599ZzjGrz6E/edrrrmbAC9HtLFWSP/wu5exoX4YHOCJbak5CPJR492tZ+q13srS3tSisLQcoe8k3Lv4Dg91dsWBP2r3j5X+HV0wNNALfwvxhbuDUrobNi3bAH83ezja2mD1bxkwFJfh41/OVlm+vbMdPnwmGB/G/45T2QYYhcDNUuM91/vB08EI8HbEv348ieNZ+VXm1+fUMJljqGqCGipU7U7Pw4T/HoW6V6i6042ScrRVWjeZO6QaWuV/9aYuiqhyF50QAjdKymW5u46ansKScti3gt/3cqMJPxy9BG8nW/xx+QaKykwY3bc94k/n4nddAcaGdUCg9+2/g+pyFCp5QTic7ZXINRTDW22LojIjyowmONsr8cupXEz+quqNGY1lTD8/zAx/ANqP96Kg+NbpzemDu2DesEAAQLd/bUOp8fZNLVMGdsIz/fxqPML90sMdseivPQDc2kc21reCuCVtfnUALl69icKScuz5PQ9P9/PDkABPs5q9v1/GR7/8jmOZ+QBu/+OR6o+hqglqiqGqtdHfLMO5yzfQt4Mz/5Ih+i8hBC5cvQmTEPjzehE0nd2qnKYGgJ9nPILevs73HC/p/FVsPpGNrw9mws7GGsXlRggBvP+33ni6nx/0RWVIy9YjrJMbrK0UKCwpx5bUnBqPsHk6qvDN1Ifg72oPK4UCf1y5gXZOdvifnenwUtsiO78IIf4uGNmnvbRMabkJNtaKKn/O/7h8Ayoba7Ojh5cLSqRnhAHAphkD0MvXqcbt0+mL8e/Np2BtpcD+81fNjvrX1uAAD7NTofdS16NN+qIyBL+5E8Ctu5Af687PhfvBUNUEMVQRUXOSkpWPHWk6ZFwuxMfP9qlykb3cMq/exNHM6+jk3hbO9jb44tcMhHV2xV96+zToeoFbwdIkcN/P5RNC4PKNEvzj+1QYTQIPdnBGP39XvPC/B6WaP959Qro+8k5XbpRA1cYKDqo2UCgUMJlEjbX3UnEEsoOr+XWvVHe1/fxuORdWEBGRrPr4OaNPLa4dlEsHN3uzGyj+ParxHg+gUChgLcPBa4VCAU9HW6x+qb/Z9LQ3tUjLNqB/R5e7HiW/8xrL+gaqW70AQqBWj7UheVju3lAiIqJWoq2qDUI7uTbqZQezHuvWaOtqCt7fcQb/+P4Eki/W7saChsBQRURE1AI93c9X+lmnl+f5dk3VZ4nnsXz3eXx7OEu6QN8SGKqIiIhaIJ9KF+P/lHLpLpXyKzeakH/z3t/3WR/FZUazhwv/c2MqYrfdfvTG0/0s98BTXlNFRETUwv169gpefrThHwIavT4FSeevIqfSkbEVz/dFeHcvKNvU7ziOobgM9jbWMAqBv3zyG87m3aixdv7wwGqfWddYGKqIiIhaqM7ubfHHlUJcuFrzd3LK5bn/HKj2IbGv/PcL5r+eHIZHurrXerzK34pQG2snh+HhOozfEHj6j4iIqIV6uKsbAODP60UNvq7Kgaq6Ozef/+Igxq8+VKuxdp/Ju2ugWvxUL/T+77PElNZWeOnhjhYPVACPVBEREbVYtx7WWvfvNq2tGyXlUFpbIfCN2w+LXfpsH4zs0x7jHrr11Wcz1h7F5hM5AIDE3y+j4z+23PVhpv934CIW/HhSeu/naodAbzX+GuyDsE6u8FTbAkC9v0+zITFUERERtVA9fWp+Mvz9+p8d6Vi2+1yV6ZWfbA8Ay57ri4EPZGLe96nStIoHky4aEYRnQztID5bNv1lqFqg+HxeCiB7eDdF+g+DpPyIiohaqk3tb6ees+3gI6LHM6/hwZzqy82+dRgxYsK3aQJUS83i1y4/p3wEn39RWmb5o0ykEvrEdvRbtQEm5EX3eipfmvfRwx2YVqAAeqSIiImqx7JS3v1poR5oOkwd2rvMYQgg8uWI/AOCTXVWDVIWe7dVwtlfWON9B1QarXuiLN35Kw+UC8+9LLCguR8CC7dL7QG9H6cusmxOGKiIiolZg04mcOocqk0mg8z+31jg/I/YJALe+Dqc2X6kzrGc7DOvZTnp/55dZV9g+a1Cd+mwqePqPiIioBfNS3/o+weNZ+XVetqZANayHN7bNHAiFQgGFQlHv7yj0cFThyzu+J7EiqDVHPFJFRETUgj35oC9WJZ6v83IVF5NXuNsde/djSKAnMmKfgL6o7K6nD5sDHqkiIiJqwR7t5nHfYzRUoKqgUCiafaACGKqIiIhatGC/249VuFjLJ6v/cipX+jnh9Udl76mlYqgiIiJqweyVt6/0ia8Ulu5m8le3n2bexcNB9p5aKoYqIiKiVqK67+Yj+fBCdSIiohauo5s9Lly9ib2/X66x5o/LN7AjLRdndAZp2lcTQxujvRaDoYqIiKiF03Rxw4WrN1FuEtXOT8nKx6jl+6pMHyTDRe6tCU//ERERtXCPBXpJP5eWm6rMry5QRQ3p0qA9tUQ8UkVERNTCDQn0lH5edzgT4zQdAQAXrhTidM7t032qNlZY/7IGHVzt4dK2+T/ioLExVBEREbVw1pWeeL542xmM03TE06v24/CF62Z1u/4+GO2d7Rq7vRaDp/+IiIhaARd7GwBAYakRe9LzqgSqvh2cGajuE49UERERtQKvRwRgwY8nAQAvfXlYmr5mYqgsT10nHqkiIiJqFZ7p51dlWif3tgxUMmKoIiIiagWUbcw/8l94qAN2/32wZZppoXj6j4iIqJVIiXkc//tbBnr7OuPxIK97L0B1wlBFRETUSjjbK/F6RICl22ixePqPiIiISAYMVUREREQyYKgiIiIikoFFQ9XKlSvRu3dvqNVqqNVqaDQabNu2TZpfXFyMqKgouLm5wcHBAaNHj0Zubq7ZGJmZmYiMjIS9vT08PT0xZ84clJeXm9Xs2bMHffv2hUqlQteuXREXF1ell+XLl6Njx46wtbVFWFgYDh06ZDa/Nr0QERFR62XRUOXr64vFixcjOTkZR44cwdChQzFy5EikpaUBAGbPno1NmzZhw4YNSExMRHZ2Np566ilpeaPRiMjISJSWlmL//v1Ys2YN4uLiEBMTI9VkZGQgMjISQ4YMQUpKCmbNmoXJkydjx44dUs26desQHR2NhQsX4ujRowgODoZWq0VeXp5Uc69eiIiIqJUTTYyLi4v44osvRH5+vrCxsREbNmyQ5p0+fVoAEElJSUIIIbZu3SqsrKyETqeTalauXCnUarUoKSkRQggxd+5c0aNHD7N1jBkzRmi1Wul9aGioiIqKkt4bjUbh4+MjYmNjhRCiVr3Uhl6vFwCEXq+v9TK1setMrvCft1n4z9ssEk7r7r0AERER1VptP7+bzDVVRqMR3377LQoLC6HRaJCcnIyysjKEh4dLNYGBgejQoQOSkpIAAElJSejVqxe8vG4/a0Or1cJgMEhHu5KSkszGqKipGKO0tBTJyclmNVZWVggPD5dqatNLdUpKSmAwGMxeRERE1DJZPFSlpqbCwcEBKpUK06ZNw8aNGxEUFASdTgelUglnZ2ezei8vL+h0OgCATqczC1QV8yvm3a3GYDCgqKgIV65cgdForLam8hj36qU6sbGxcHJykl5+flW/IoCIiIhaBouHqoCAAKSkpODgwYOYPn06xo8fj1OnTlm6LVnMnz8fer1eemVlZVm6JSIiImogFn+iulKpRNeuXQEAISEhOHz4MJYuXYoxY8agtLQU+fn5ZkeIcnNz4e3tDQDw9vaucpdexR15lWvuvEsvNzcXarUadnZ2sLa2hrW1dbU1lce4Vy/VUalUUKlUddgbRERE1FxZ/EjVnUwmE0pKShASEgIbGxskJCRI89LT05GZmQmNRgMA0Gg0SE1NNbtLLz4+Hmq1GkFBQVJN5TEqairGUCqVCAkJMasxmUxISEiQamrTCxEREbVuFj1SNX/+fAwfPhwdOnRAQUEB1q5diz179mDHjh1wcnLCpEmTEB0dDVdXV6jVarz66qvQaDR46KGHAAAREREICgrCuHHjsGTJEuh0OixYsABRUVHSEaJp06Zh2bJlmDt3LiZOnIhdu3Zh/fr12LJli9RHdHQ0xo8fj379+iE0NBQff/wxCgsLMWHCBACoVS9ERETUulk0VOXl5eHFF19ETk4OnJyc0Lt3b+zYsQOPP/44AOCjjz6ClZUVRo8ejZKSEmi1WqxYsUJa3traGps3b8b06dOh0WjQtm1bjB8/Hm+99ZZU06lTJ2zZsgWzZ8/G0qVL4evriy+++AJarVaqGTNmDC5fvoyYmBjodDr06dMH27dvN7t4/V69EBERUeumEEIISzfRWhgMBjg5OUGv10OtVss27u70PEz48jAAYPVL/TA00OseSxAREVFt1fbzu8ldU0VERETUHDFUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkwFBFREREJAOGKiIiIiIZMFQRERERyYChioiIiEgGDFVEREREMmCoIiIiIpIBQxURERGRDBiqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkwFBFREREJAOGKiIiIiIZMFQRERERyYChioiIiEgGDFVEREREMmCoIiIiIpKBRUNVbGws+vfvD0dHR3h6emLUqFFIT083qxk8eDAUCoXZa9q0aWY1mZmZiIyMhL29PTw9PTFnzhyUl5eb1ezZswd9+/aFSqVC165dERcXV6Wf5cuXo2PHjrC1tUVYWBgOHTpkNr+4uBhRUVFwc3ODg4MDRo8ejdzcXHl2BhERETVrFg1ViYmJiIqKwoEDBxAfH4+ysjJERESgsLDQrG7KlCnIycmRXkuWLJHmGY1GREZGorS0FPv378eaNWsQFxeHmJgYqSYjIwORkZEYMmQIUlJSMGvWLEyePBk7duyQatatW4fo6GgsXLgQR48eRXBwMLRaLfLy8qSa2bNnY9OmTdiwYQMSExORnZ2Np556qgH3UN2ZTJbugIiIqJUSTUheXp4AIBITE6Vpjz76qJg5c2aNy2zdulVYWVkJnU4nTVu5cqVQq9WipKRECCHE3LlzRY8ePcyWGzNmjNBqtdL70NBQERUVJb03Go3Cx8dHxMbGCiGEyM/PFzY2NmLDhg1SzenTpwUAkZSUVKvt0+v1AoDQ6/W1qq+tXWdyhf+8zcJ/3mZxJscg69hEREStXW0/v5vUNVV6vR4A4Orqajb966+/hru7O3r27In58+fj5s2b0rykpCT06tULXl5e0jStVguDwYC0tDSpJjw83GxMrVaLpKQkAEBpaSmSk5PNaqysrBAeHi7VJCcno6yszKwmMDAQHTp0kGruVFJSAoPBYPZqaApFg6+CiIiIqtHG0g1UMJlMmDVrFh555BH07NlTmv7cc8/B398fPj4+OHHiBObNm4f09HT88MMPAACdTmcWqABI73U63V1rDAYDioqKcP36dRiNxmprzpw5I42hVCrh7OxcpaZiPXeKjY3Fm2++Wcc9QURERM1RkwlVUVFROHnyJH777Tez6VOnTpV+7tWrF9q1a4fHHnsM58+fR5cuXRq7zTqZP38+oqOjpfcGgwF+fn4W7IiIiIgaSpM4/Tdjxgxs3rwZu3fvhq+v711rw8LCAADnzp0DAHh7e1e5A6/ivbe3911r1Go17Ozs4O7uDmtr62prKo9RWlqK/Pz8GmvupFKpoFarzV5ERETUMlk0VAkhMGPGDGzcuBG7du1Cp06d7rlMSkoKAKBdu3YAAI1Gg9TUVLO79OLj46FWqxEUFCTVJCQkmI0THx8PjUYDAFAqlQgJCTGrMZlMSEhIkGpCQkJgY2NjVpOeno7MzEyphoiIiFovi57+i4qKwtq1a/HTTz/B0dFRujbJyckJdnZ2OH/+PNauXYsnnngCbm5uOHHiBGbPno1Bgwahd+/eAICIiAgEBQVh3LhxWLJkCXQ6HRYsWICoqCioVCoAwLRp07Bs2TLMnTsXEydOxK5du7B+/Xps2bJF6iU6Ohrjx49Hv379EBoaio8//hiFhYWYMGGC1NOkSZMQHR0NV1dXqNVqvPrqq9BoNHjooYcaec8RERFRk9M4NyNWD0C1ry+//FIIIURmZqYYNGiQcHV1FSqVSnTt2lXMmTOnyi2NFy5cEMOHDxd2dnbC3d1dvP7666KsrMysZvfu3aJPnz5CqVSKzp07S+uo7NNPPxUdOnQQSqVShIaGigMHDpjNLyoqEq+88opwcXER9vb24sknnxQ5OTm13t7GeKRCuo6PVCAiIpJTbT+/FUIIYblI17oYDAY4OTlBr9fLen3V7vQ8TPjyMABg5+xB6OblKNvYRERErV1tP7+bxIXqRERERM0dQxURERGRDBiqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkwFBFREREJAOGKiIiIiIZMFQRERERyYChioiIiEgGDFVEREREMmCoIiIiIpIBQxURERGRDBiqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGVg0VMXGxqJ///5wdHSEp6cnRo0ahfT0dLOa4uJiREVFwc3NDQ4ODhg9ejRyc3PNajIzMxEZGQl7e3t4enpizpw5KC8vN6vZs2cP+vbtC5VKha5duyIuLq5KP8uXL0fHjh1ha2uLsLAwHDp0qM69EBERUetk0VCVmJiIqKgoHDhwAPHx8SgrK0NERAQKCwulmtmzZ2PTpk3YsGEDEhMTkZ2djaeeekqabzQaERkZidLSUuzfvx9r1qxBXFwcYmJipJqMjAxERkZiyJAhSElJwaxZszB58mTs2LFDqlm3bh2io6OxcOFCHD16FMHBwdBqtcjLy6t1L0RERNSKiSYkLy9PABCJiYlCCCHy8/OFjY2N2LBhg1Rz+vRpAUAkJSUJIYTYunWrsLKyEjqdTqpZuXKlUKvVoqSkRAghxNy5c0WPHj3M1jVmzBih1Wql96GhoSIqKkp6bzQahY+Pj4iNja11L/ei1+sFAKHX62tVX1u7zuQK/3mbhf+8zSJdZ5B1bCIiotautp/fTeqaKr1eDwBwdXUFACQnJ6OsrAzh4eFSTWBgIDp06ICkpCQAQFJSEnr16gUvLy+pRqvVwmAwIC0tTaqpPEZFTcUYpaWlSE5ONquxsrJCeHi4VFObXu5UUlICg8Fg9iIiIqKWqcmEKpPJhFmzZuGRRx5Bz549AQA6nQ5KpRLOzs5mtV5eXtDpdFJN5UBVMb9i3t1qDAYDioqKcOXKFRiNxmprKo9xr17uFBsbCycnJ+nl5+dXy71BREREzU2TCVVRUVE4efIkvv32W0u3Ipv58+dDr9dLr6ysLEu3RERERA2kjaUbAIAZM2Zg8+bN2Lt3L3x9faXp3t7eKC0tRX5+vtkRotzcXHh7e0s1d96lV3FHXuWaO+/Sy83NhVqthp2dHaytrWFtbV1tTeUx7tXLnVQqFVQqVR32RP3YtrFu8HUQERHR3Vn0SJUQAjNmzMDGjRuxa9cudOrUyWx+SEgIbGxskJCQIE1LT09HZmYmNBoNAECj0SA1NdXsLr34+Hio1WoEBQVJNZXHqKipGEOpVCIkJMSsxmQyISEhQaqpTS+W4uGotOj6iYiIqJ5HqrKysqBQKKSjSocOHcLatWsRFBSEqVOn1nqcqKgorF27Fj/99BMcHR2la5OcnJxgZ2cHJycnTJo0CdHR0XB1dYVarcarr74KjUaDhx56CAAQERGBoKAgjBs3DkuWLIFOp8OCBQsQFRUlHSWaNm0ali1bhrlz52LixInYtWsX1q9fjy1btki9REdHY/z48ejXrx9CQ0Px8ccfo7CwEBMmTJB6ulcvRERE1IrV59bCAQMGiK+++koIIUROTo5Qq9VCo9EId3d38eabb9Z6HADVvr788kuppqioSLzyyivCxcVF2NvbiyeffFLk5OSYjXPhwgUxfPhwYWdnJ9zd3cXrr78uysrKzGp2794t+vTpI5RKpejcubPZOip8+umnokOHDkKpVIrQ0FBx4MABs/m16eVuGuqRCmdzDXykAhERUQOp7ee3Qggh6hrEXFxccODAAQQEBOCTTz7BunXrsG/fPuzcuRPTpk3DH3/8IWfuazEMBgOcnJyg1+uhVqtlG/dcXgHCP9wLANg5exC6eTnKNjYREVFrV9vP73pdU1VWViadWvvll1/w17/+FcCt5zbl5OTUZ0giIiKiZq1eoapHjx5YtWoVfv31V8THx2PYsGEAgOzsbLi5ucnaIBEREVFzUK9Q9d577+Gzzz7D4MGDMXbsWAQHBwMAfv75Z4SGhsraIBEREVFzUK+7/wYPHowrV67AYDDAxcVFmj516lS0bdtWtuaIiIiImot6HakaOnQoCgoKzAIVcOs7+8aMGSNLY0RERETNSb1C1Z49e1BaWlplenFxMX799df7boqIiIiouanT6b8TJ05IP586dcrsi4SNRiO2b9+O9u3by9cdERERUTNRp1DVp08fKBQKKBQKDB06tMp8Ozs7fPrpp7I1R0RERNRc1ClUZWRkQAiBzp0749ChQ/Dw8JDmKZVKeHp6wtqaX+5LRERErU+dQpW/vz+AW182TERERES31euRCgBw9uxZ7N69G3l5eVVCVkxMzH03RkRERNSc1CtU/ec//8H06dPh7u4Ob29vKBQKaZ5CoWCoIiIiolanXqHq7bffxjvvvIN58+bJ3Q8RERFRs1Sv51Rdv34dTz/9tNy9EBERETVb9QpVTz/9NHbu3Cl3L0RERETNVr1O/3Xt2hVvvPEGDhw4gF69esHGxsZs/muvvSZLc0RERETNRb1C1eeffw4HBwckJiYiMTHRbJ5CoWCoIiIiolanXqEqIyND7j6IiIiImrV6XVNFRERERObqdaRq4sSJd52/evXqejVDRERE1FzVK1Rdv37d7H1ZWRlOnjyJ/Pz8ar9omYiIiKilq1eo2rhxY5VpJpMJ06dPR5cuXe67KSIiIqLmRrZrqqysrBAdHY2PPvpIriGJiIiImg1ZL1Q/f/48ysvL5RySiIiIqFmo1+m/6Ohos/dCCOTk5GDLli0YP368LI0RERERNSf1ClXHjh0ze29lZQUPDw988MEH97wzkIiIiKglqleo2r17t9x9EBERETVr9QpVFS5fvoz09HQAQEBAADw8PGRpioiIiKi5qdeF6oWFhZg4cSLatWuHQYMGYdCgQfDx8cGkSZNw8+ZNuXskIiIiavLqFaqio6ORmJiITZs2IT8/H/n5+fjpp5+QmJiI119/Xe4eiYiIiJq8ep3++/777/Hdd99h8ODB0rQnnngCdnZ2eOaZZ7By5Uq5+iMiIiJqFup1pOrmzZvw8vKqMt3T05On/4iIiKhVqleo0mg0WLhwIYqLi6VpRUVFePPNN6HRaGRrjoiIiKi5qNfpv48//hjDhg2Dr68vgoODAQDHjx+HSqXCzp07ZW2QiIiIqDmoV6jq1asXzp49i6+//hpnzpwBAIwdOxbPP/887OzsZG2QiIiIqDmoV6iKjY2Fl5cXpkyZYjZ99erVuHz5MubNmydLc0RERETNRb2uqfrss88QGBhYZXqPHj2watWq+26KiIiIqLmpV6jS6XRo165dlekeHh7Iycm576aIiIiImpt6hSo/Pz/s27evyvR9+/bBx8en1uPs3bsXI0aMgI+PDxQKBX788Uez+S+99BIUCoXZa9iwYWY1165dw/PPPw+1Wg1nZ2dMmjQJN27cMKs5ceIEBg4cCFtbW/j5+WHJkiVVetmwYQMCAwNha2uLXr16YevWrWbzhRCIiYlBu3btYGdnh/DwcJw9e7bW20pEREQtW71C1ZQpUzBr1ix8+eWXuHjxIi5evIjVq1dj9uzZVa6zupvCwkIEBwdj+fLlNdYMGzYMOTk50uubb74xm//8888jLS0N8fHx2Lx5M/bu3YupU6dK8w0GAyIiIuDv74/k5GS8//77WLRoET7//HOpZv/+/Rg7diwmTZqEY8eOYdSoURg1ahROnjwp1SxZsgSffPIJVq1ahYMHD6Jt27bQarVmj5UgIiKiVkzUg8lkEnPnzhW2trbCyspKWFlZCXt7e/Hmm2/WZzghhBAAxMaNG82mjR8/XowcObLGZU6dOiUAiMOHD0vTtm3bJhQKhbh06ZIQQogVK1YIFxcXUVJSItXMmzdPBAQESO+feeYZERkZaTZ2WFiYePnll4UQt7bX29tbvP/++9L8/Px8oVKpxDfffFNjf8XFxUKv10uvrKwsAUDo9fqad0Q9nM01CP95m4X/vM0iXWeQdWwiIqLWTq/X1+rzu15HqhQKBd577z1cvnwZBw4cwPHjx3Ht2jXExMTImfcAAHv27IGnpycCAgIwffp0XL16VZqXlJQEZ2dn9OvXT5oWHh4OKysrHDx4UKoZNGgQlEqlVKPVapGeno7r169LNeHh4Wbr1Wq1SEpKAgBkZGRAp9OZ1Tg5OSEsLEyqqU5sbCycnJykl5+f333sCSIiImrK6hWqKjg4OKB///7o2bMnVCqVXD1Jhg0bhq+++goJCQl47733kJiYiOHDh8NoNAK4dcG8p6en2TJt2rSBq6srdDqdVHPnV+pUvL9XTeX5lZerrqY68+fPh16vl15ZWVl12n4iIiJqPur1nKrG8uyzz0o/9+rVC71790aXLl2wZ88ePPbYYxbsrHZUKlWDhE0iIiJqeu7rSFVj69y5M9zd3XHu3DkAgLe3N/Ly8sxqysvLce3aNXh7e0s1ubm5ZjUV7+9VU3l+5eWqqyEiIqLWrVmFqj///BNXr16VnpGl0WiQn5+P5ORkqWbXrl0wmUwICwuTavbu3YuysjKpJj4+HgEBAXBxcZFqEhISzNYVHx8vfTl0p06d4O3tbVZjMBhw8OBBfoE0ERERAbBwqLpx4wZSUlKQkpIC4NYF4SkpKcjMzMSNGzcwZ84cHDhwABcuXEBCQgJGjhyJrl27QqvVAgC6d++OYcOGYcqUKTh06BD27duHGTNm4Nlnn5Wel/Xcc89BqVRi0qRJSEtLw7p167B06VJER0dLfcycORPbt2/HBx98gDNnzmDRokU4cuQIZsyYAeDWhfmzZs3C22+/jZ9//hmpqal48cUX4ePjg1GjRjXqPiMiIqImqpHuRqzW7t27BYAqr/Hjx4ubN2+KiIgI4eHhIWxsbIS/v7+YMmWK0Ol0ZmNcvXpVjB07Vjg4OAi1Wi0mTJggCgoKzGqOHz8uBgwYIFQqlWjfvr1YvHhxlV7Wr18vunXrJpRKpejRo4fYsmWL2XyTySTeeOMN4eXlJVQqlXjsscdEenp6nba3trdk1hUfqUBERNRwavv5rRBCCAtmulbFYDDAyckJer0earVatnHP5RUg/MO9AICdswehm5ejbGMTERG1drX9/G5W11QRERERNVUMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkwFBFREREJAOGKiIiIiIZMFQRERERyYChioiIiEgGDFVEREREMmCoIiIiIpIBQxURERGRDBiqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkYNFQtXfvXowYMQI+Pj5QKBT48ccfzeYLIRATE4N27drBzs4O4eHhOHv2rFnNtWvX8Pzzz0OtVsPZ2RmTJk3CjRs3zGpOnDiBgQMHwtbWFn5+fliyZEmVXjZs2IDAwEDY2tqiV69e2Lp1a517ISIiotbLoqGqsLAQwcHBWL58ebXzlyxZgk8++QSrVq3CwYMH0bZtW2i1WhQXF0s1zz//PNLS0hAfH4/Nmzdj7969mDp1qjTfYDAgIiIC/v7+SE5Oxvvvv49Fixbh888/l2r279+PsWPHYtKkSTh27BhGjRqFUaNG4eTJk3XqhYiIiFox0UQAEBs3bpTem0wm4e3tLd5//31pWn5+vlCpVOKbb74RQghx6tQpAUAcPnxYqtm2bZtQKBTi0qVLQgghVqxYIVxcXERJSYlUM2/ePBEQECC9f+aZZ0RkZKRZP2FhYeLll1+udS/VKS4uFnq9XnplZWUJAEKv19dl19zT2VyD8J+3WfjP2yzSdQZZxyYiImrt9Hp9rT6/m+w1VRkZGdDpdAgPD5emOTk5ISwsDElJSQCApKQkODs7o1+/flJNeHg4rKyscPDgQalm0KBBUCqVUo1Wq0V6ejquX78u1VReT0VNxXpq00t1YmNj4eTkJL38/PzquzuIiIioiWuyoUqn0wEAvLy8zKZ7eXlJ83Q6HTw9Pc3mt2nTBq6urmY11Y1ReR011VSef69eqjN//nzo9XrplZWVdY+tJiIiouaqjaUbaMlUKhVUKpWl2yAiIqJG0GSPVHl7ewMAcnNzzabn5uZK87y9vZGXl2c2v7y8HNeuXTOrqW6Myuuoqaby/Hv1QkRERK1bkw1VnTp1gre3NxISEqRpBoMBBw8ehEajAQBoNBrk5+cjOTlZqtm1axdMJhPCwsKkmr1796KsrEyqiY+PR0BAAFxcXKSayuupqKlYT216ISIiotbNoqHqxo0bSElJQUpKCoBbF4SnpKQgMzMTCoUCs2bNwttvv42ff/4ZqampePHFF+Hj44NRo0YBALp3745hw4ZhypQpOHToEPbt24cZM2bg2WefhY+PDwDgueeeg1KpxKRJk5CWloZ169Zh6dKliI6OlvqYOXMmtm/fjg8++ABnzpzBokWLcOTIEcyYMQMAatULERERtXKNdDditXbv3i0AVHmNHz9eCHHrUQZvvPGG8PLyEiqVSjz22GMiPT3dbIyrV6+KsWPHCgcHB6FWq8WECRNEQUGBWc3x48fFgAEDhEqlEu3btxeLFy+u0sv69etFt27dhFKpFD169BBbtmwxm1+bXu6ltrdk1hUfqUBERNRwavv5rRBCCAtmulbFYDDAyckJer0earVatnHP5RUg/MO9AICdswehm5ejbGMTERG1drX9/G6y11QRERERNScMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkwFBFREREJAOGKiIiIiIZMFQRERERyYChioiIiEgGDFVEREREMmCoIiIiIpIBQxURERGRDBiqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhk0KRD1aJFi6BQKMxegYGB0vzi4mJERUXBzc0NDg4OGD16NHJzc83GyMzMRGRkJOzt7eHp6Yk5c+agvLzcrGbPnj3o27cvVCoVunbtiri4uCq9LF++HB07doStrS3CwsJw6NChBtlmIiIiap6adKgCgB49eiAnJ0d6/fbbb9K82bNnY9OmTdiwYQMSExORnZ2Np556SppvNBoRGRmJ0tJS7N+/H2vWrEFcXBxiYmKkmoyMDERGRmLIkCFISUnBrFmzMHnyZOzYsUOqWbduHaKjo7Fw4UIcPXoUwcHB0Gq1yMvLa5ydQERERE2faMIWLlwogoODq52Xn58vbGxsxIYNG6Rpp0+fFgBEUlKSEEKIrVu3CisrK6HT6aSalStXCrVaLUpKSoQQQsydO1f06NHDbOwxY8YIrVYrvQ8NDRVRUVHSe6PRKHx8fERsbGydtkev1wsAQq/X12m5ezmbaxD+8zYL/3mbRbrOIOvYRERErV1tP7+b/JGqs2fPwsfHB507d8bzzz+PzMxMAEBycjLKysoQHh4u1QYGBqJDhw5ISkoCACQlJaFXr17w8vKSarRaLQwGA9LS0qSaymNU1FSMUVpaiuTkZLMaKysrhIeHSzU1KSkpgcFgMHsRERFRy9SkQ1VYWBji4uKwfft2rFy5EhkZGRg4cCAKCgqg0+mgVCrh7OxstoyXlxd0Oh0AQKfTmQWqivkV8+5WYzAYUFRUhCtXrsBoNFZbUzFGTWJjY+Hk5CS9/Pz86rwPiIiIqHloY+kG7mb48OHSz71790ZYWBj8/f2xfv162NnZWbCz2pk/fz6io6Ol9waDgcGKiIiohWrSR6ru5OzsjG7duuHcuXPw9vZGaWkp8vPzzWpyc3Ph7e0NAPD29q5yN2DF+3vVqNVq2NnZwd3dHdbW1tXWVIxRE5VKBbVabfYiIiKilqlZhaobN27g/PnzaNeuHUJCQmBjY4OEhARpfnp6OjIzM6HRaAAAGo0GqampZnfpxcfHQ61WIygoSKqpPEZFTcUYSqUSISEhZjUmkwkJCQlSDREREVGTDlV///vfkZiYiAsXLmD//v148sknYW1tjbFjx8LJyQmTJk1CdHQ0du/ejeTkZEyYMAEajQYPPfQQACAiIgJBQUEYN24cjh8/jh07dmDBggWIioqCSqUCAEybNg1//PEH5s6dizNnzmDFihVYv349Zs+eLfURHR2N//znP1izZg1Onz6N6dOno7CwEBMmTLDIfiEiIqKmp0lfU/Xnn39i7NixuHr1Kjw8PDBgwAAcOHAAHh4eAICPPvoIVlZWGD16NEpKSqDVarFixQppeWtra2zevBnTp0+HRqNB27ZtMX78eLz11ltSTadOnbBlyxbMnj0bS5cuha+vL7744gtotVqpZsyYMbh8+TJiYmKg0+nQp08fbN++vcrF60RERNR6KYQQwtJNtBYGgwFOTk7Q6/WyXl91Lq8A4R/uBQDsnD0I3bwcZRubiIiotavt53eTPv1HRERE1FwwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUMVERERkQwYqoiIiIhkwFBFREREJAOGKiIiIiIZMFQRERERyYChioiIiEgGDFVEREREMmCoIiIiIpIBQxURERGRDBiqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZMBQRURERCQDhioiIiIiGTBUEREREcmAoYqIiIhIBgxVRERERDJgqCIiIiKSAUNVHS1fvhwdO3aEra0twsLCcOjQIUu3RERERE0AQ1UdrFu3DtHR0Vi4cCGOHj2K4OBgaLVa5OXlWbo1IiIisjCGqjr48MMPMWXKFEyYMAFBQUFYtWoV7O3tsXr1aku3RkRERBbWxtINNBelpaVITk7G/PnzpWlWVlYIDw9HUlJStcuUlJSgpKREem8wGBqkN3+3tujbwRkmAXR0a9sg6yAiIqK7Y6iqpStXrsBoNMLLy8tsupeXF86cOVPtMrGxsXjzzTcbvDcbayv88MojDb4eIiIiqhlP/zWg+fPnQ6/XS6+srCxLt0REREQNhEeqasnd3R3W1tbIzc01m56bmwtvb+9ql1GpVFCpVI3RHhEREVkYj1TVklKpREhICBISEqRpJpMJCQkJ0Gg0FuyMiIiImgIeqaqD6OhojB8/Hv369UNoaCg+/vhjFBYWYsKECZZujYiIiCyMoaoOxowZg8uXLyMmJgY6nQ59+vTB9u3bq1y8TkRERK2PQgghLN1Ea2EwGODk5AS9Xg+1Wm3pdoiIiKgWavv5zWuqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEM+ET1RlTxnFWDwWDhToiIiKi2Kj637/W8dIaqRlRQUAAA8PPzs3AnREREVFcFBQVwcnKqcT6/pqYRmUwmZGdnw9HREQqFwtLtNDkGgwF+fn7Iysri1/g0EO7jhsd93PC4jxse97E5IQQKCgrg4+MDK6uar5zikapGZGVlBV9fX0u30eSp1Wr+IW5g3McNj/u44XEfNzzu49vudoSqAi9UJyIiIpIBQxURERGRDBiqqMlQqVRYuHAhVCqVpVtpsbiPGx73ccPjPm543Mf1wwvViYiIiGTAI1VEREREMmCoIiIiIpIBQxURERGRDBiqiIiIiGTAUEWyiY2NRf/+/eHo6AhPT0+MGjUK6enpZjUvv/wyunTpAjs7O3h4eGDkyJE4c+aMWc1rr72GkJAQqFQq9OnTp9p1nThxAgMHDoStrS38/PywZMmShtqsJqWx9vGFCxegUCiqvA4cONCQm9ckyLGPjx8/jrFjx8LPzw92dnbo3r07li5dWmVde/bsQd++faFSqdC1a1fExcU19OY1CY21j/fs2VPt77FOp2uU7bQkOfbx1atXMWzYMPj4+EClUsHPzw8zZsyo8v21rfX3uDoMVSSbxMREREVF4cCBA4iPj0dZWRkiIiJQWFgo1YSEhODLL7/E6dOnsWPHDgghEBERAaPRaDbWxIkTMWbMmGrXYzAYEBERAX9/fyQnJ+P999/HokWL8Pnnnzfo9jUFjbWPK/zyyy/IycmRXiEhIQ2yXU2JHPs4OTkZnp6e+L//+z+kpaXhX//6F+bPn49ly5ZJY2RkZCAyMhJDhgxBSkoKZs2ahcmTJ2PHjh2Nvs2NrbH2cYX09HSz32NPT89G21ZLkWMfW1lZYeTIkfj555/x+++/Iy4uDr/88gumTZsmjdGaf4+rJYgaSF5engAgEhMTa6w5fvy4ACDOnTtXZd7ChQtFcHBwlekrVqwQLi4uoqSkRJo2b948ERAQIEvfzUlD7eOMjAwBQBw7dkzGbpun+93HFV555RUxZMgQ6f3cuXNFjx49zGrGjBkjtFrt/TfdzDTUPt69e7cAIK5fvy5nu82SXPt46dKlwtfXV3rP32NzPFJFDUav1wMAXF1dq51fWFiIL7/8Ep06dYKfn1+tx01KSsKgQYOgVCqlaVqtFunp6bh+/fr9Nd3MNNQ+rvDXv/4Vnp6eGDBgAH7++ef76rW5kmsf6/V6szGSkpIQHh5uVqPVapGUlCRD181LQ+3jCn369EG7du3w+OOPY9++ffI03czIsY+zs7Pxww8/4NFHH5Wm8ffYHEMVNQiTyYRZs2bhkUceQc+ePc3mrVixAg4ODnBwcMC2bdsQHx9vFpDuRafTwcvLy2xaxfvWcK1EhYbcxw4ODvjggw+wYcMGbNmyBQMGDMCoUaNaXbCSax/v378f69atw9SpU6VpNf0eGwwGFBUVyb8xTVRD7uN27dph1apV+P777/H999/Dz88PgwcPxtGjRxt0m5qa+93HY8eOhb29Pdq3bw+1Wo0vvvhCmsff4ztY+lAZtUzTpk0T/v7+Iisrq8q8/Px88fvvv4vExEQxYsQI0bdvX1FUVFSlrqZTU48//riYOnWq2bS0tDQBQJw6dUq2bWjqGnIfV2fcuHFiwIAB99t2syLHPk5NTRXu7u7i3//+t9n0Bx54QLz77rtm07Zs2SIAiJs3b8q7IU1YQ+7j6gwaNEi88MILsvTeXNzvPs7JyRGnT58WP/30kwgKChLTp0+X5vH32BxDFckuKipK+Pr6ij/++OOetSUlJcLe3l6sXbu2yryaPvDHjRsnRo4caTZt165dAoC4du1afdtuVhp6H1dn2bJlwtvbu66tNlty7OO0tDTh6ekp/vnPf1ZZZuDAgWLmzJlm01avXi3UavV99d2cNPQ+rs7f//538dBDD9Wr3+ZIrr8rKvz6668CgMjOzhZC8Pf4Tjz9R7IRQmDGjBnYuHEjdu3ahU6dOtVqGSEESkpKar0ejUaDvXv3oqysTJoWHx+PgIAAuLi41Kv35qKx9nF1UlJS0K5du/saozmQax+npaVhyJAhGD9+PN55550qy2g0GiQkJJhNi4+Ph0ajuf+NaOIaax9Xh7/Hd1/mXn9XmEwmAJBqWvPvcbUsEOSohZo+fbpwcnISe/bsETk5OdKr4hDw+fPnxbvvviuOHDkiLl68KPbt2ydGjBghXF1dRW5urjTO2bNnxbFjx8TLL78sunXrJo4dOyaOHTsm3e2Xn58vvLy8xLhx48TJkyfFt99+K+zt7cVnn31mke1uTI21j+Pi4sTatWvF6dOnxenTp8U777wjrKysxOrVqy2y3Y1Jjn2cmpoqPDw8xAsvvGA2Rl5enrSeP/74Q9jb24s5c+aI06dPi+XLlwtra2uxfft2i2x3Y2qsffzRRx+JH3/8UZw9e1akpqaKmTNnCisrK/HLL79YZLsbkxz7eMuWLWL16tUiNTVVZGRkiM2bN4vu3buLRx55RFpPa/49rg5DFckGQLWvL7/8UgghxKVLl8Tw4cOFp6ensLGxEb6+vuK5554TZ86cMRvn0UcfrXacjIwMqeb48eNiwIABQqVSifbt24vFixc34pZaTmPt47i4ONG9e3dhb28v1Gq1CA0NFRs2bGjkrbUMOfbxwoULqx3D39/fbF27d+8Wffr0EUqlUnTu3FlaR0vXWPv4vffeE126dBG2trbC1dVVDB48WOzatauRt9Yy5NjHu3btEhqNRjg5OQlbW1vxwAMPiHnz5lV5REVr/T2ujkIIIeQ++kVERETU2vCaKiIiIiIZMFQRERERyYChioiIiEgGDFVEREREMmCoIiIiIpIBQxURERGRDBiqiIiIiGTAUEVEREQkA4YqImrVXnrpJYwaNcpi6x83bhzefffdWtU+++yz+OCDDxq4IyKqLz5RnYhaLIVCcdf5CxcuxOzZsyGEgLOzc+M0Vcnx48cxdOhQXLx4EQ4ODvesP3nyJAYNGoSMjAw4OTk1QodEVBcMVUTUYul0OunndevWISYmBunp6dI0BweHWoWZhjJ58mS0adMGq1atqvUy/fv3x0svvYSoqKgG7IyI6oOn/4ioxfL29pZeTk5OUCgUZtMcHByqnP4bPHgwXn31VcyaNQsuLi7w8vLCf/7zHxQWFmLChAlwdHRE165dsW3bNrN1nTx5EsOHD4eDgwO8vLwwbtw4XLlypcbejEYjvvvuO4wYMcJs+ooVK/DAAw/A1tYWXl5e+Nvf/mY2f8SIEfj222/vf+cQkewYqoiI7rBmzRq4u7vj0KFDePXVVzF9+nQ8/fTTePjhh3H06FFERERg3LhxuHnzJgAgPz8fQ4cOxYMPPogjR45g+/btyM3NxTPPPFPjOk6cOAG9Xo9+/fpJ044cOYLXXnsNb731FtLT07F9+3YMGjTIbLnQ0FAcOnQIJSUlDbPxRFRvDFVERHcIDg7GggUL8MADD2D+/PmwtbWFu7s7pkyZggceeAAxMTG4evUqTpw4AQBYtmwZHnzwQbz77rsIDAzEgw8+iNWrV2P37t34/fffq13HxYsXYW1tDU9PT2laZmYm2rZti7/85S/w9/fHgw8+iNdee81sOR8fH5SWlpqd2iSipoGhiojoDr1795Z+tra2hpubG3r16iVN8/LyAgDk5eUBuHXB+e7du6VrtBwcHBAYGAgAOH/+fLXrKCoqgkqlMruY/vHHH4e/vz86d+6McePG4euvv5aOhlWws7MDgCrTicjyGKqIiO5gY2Nj9l6hUJhNqwhCJpMJAHDjxg2MGDECKSkpZq+zZ89WOX1Xwd3dHTdv3kRpaak0zdHREUePHsU333yDdu3aISYmBsHBwcjPz5dqrl27BgDw8PCQZVuJSD4MVURE96lv375IS0tDx44d0bVrV7NX27Ztq12mT58+AIBTp06ZTW/Tpg3Cw8OxZMkSnDhxAhcuXMCuXbuk+SdPnoSvry/c3d0bbHuIqH4YqoiI7lNUVBSuXbuGsWPH4vDhwzh//jx27NiBCRMmwGg0VruMh4cH+vbti99++02atnnzZnzyySdISUnBxYsX8dVXX8FkMiEgIECq+fXXXxEREdHg20REdcdQRUR0n3x8fLBv3z4YjUZERESgV69emDVrFpydnWFlVfNfs5MnT8bXX38tvXd2dsYPP/yAoUOHonv37li1ahW++eYb9OjRAwBQXFyMH3/8EVOmTGnwbSKiuuPDP4mILKSoqAgBAQFYt24dNBrNPetXrlyJjRs3YufOnY3QHRHVFY9UERFZiJ2dHb766qu7PiS0MhsbG3z66acN3BUR1RePVBERERHJgEeqiIiIiGTAUEVEREQkA4YqIiIiIhkwVBERERHJgKGKiIiISAYMVUREREQyYKgiIiIikgFDFREREZEMGKqIiIiIZPD/AVr1YBidjf8OAAAAAElFTkSuQmCC", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "\n", + "lc = Lightcurve(tess_bjds_clean,sap_fluxes_corrected,err=sap_fluxes_error_corrected)\n", + "lc.plot()" + ] + }, + { + "cell_type": "code", + "execution_count": 150, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([353359.56, 353359.56, 353359.56, ..., 367221.88, 367348.25,\n", + " 367650.34], dtype='>f4')" + ] + }, + "execution_count": 150, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "pdcsap_fluxes_corrected = np.nan_to_num(pdcsap_fluxes,nan = np.nanmedian(pdcsap_fluxes))\n", + "pdcsap_fluxes_error_corrected = np.nan_to_num(pdcsap_fluxes_error , nan = np.nanmedian(pdcsap_fluxes_error))\n", + "\n", + "pdcsap_fluxes_error_corrected\n", + "pdcsap_fluxes_corrected\n", + "\n" + ] + }, + { + "cell_type": "code", + "execution_count": 151, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 151, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlUAAAGwCAYAAACAZ5AeAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAiINJREFUeJzt3Xd4VFX6B/DvTEgmjUkvhIQkELrUIBALAgJBWRdWVCwrLCIKC7qQVSCroosFFlcRlbKKgPtbFUQXV0GIkSYltEAgVOlFUmjJQHpm7u+Pyb25d0oySSaZku/nefIw5czMmcsk951z3vMelSAIAoiIiIioQdSO7gARERGRO2BQRURERGQHDKqIiIiI7IBBFREREZEdMKgiIiIisgMGVURERER2wKCKiIiIyA5aOLoDzYnBYMCVK1fQsmVLqFQqR3eHiIiIbCAIAm7duoWoqCio1dbHoxhUNaErV64gJibG0d0gIiKierh06RKio6Ot3s+gqgm1bNkSgPE/RavVOrg3REREZAudToeYmBjpPG4Ng6omJE75abVaBlVEREQuprbUHSaqExEREdkBgyoiIiIiO2BQRURERGQHDKqIiIiI7IBBFREREZEdMKgiIiIisgMGVURERER2wKCKiIiIyA4YVBERERHZAYMqIiIiIjtgUEVERERkBwyqiIiIiOyAQRURERG5hdzCUtwuq3TY6zOoIiIiIpc398fj6D93E/6X9ZvD+sCgioiIiFzeocsFAICtJ686rA8MqoiIiMjlHb5cCAAY2TPKYX1gUEVEREQur7hcDwAorTA4rA8MqoiIiMhttA3zc9hrM6giIiIitxEd5OOw12ZQRURERC7t+u0y6bKvVwuH9YNBFREREbm00srqPCp/DYMqIiIiogbRtHBsWMOgioiIiFxavq4UAFBpEBzaDwZVRERE5NJO5t4CAOgZVBERERHV37EcnaO7AIBBFREREbm4c9eKAACdIls6tB8MqoiIiMiltQrwVvzrKAyqiIiIyKXtPH0dANAnLtih/WBQRURERC4tQqsBABQUlzu0HwyqiIiIyKUduFgAAOgSpXVoPxhUERERkVtw5BY1AIMqIiIicmFXb1Xv++fj6eHAnjg4qFqyZAm6d+8OrVYLrVaLpKQkbNiwQdEmIyMDgwcPhp+fH7RaLQYMGICSkhLp/ri4OKhUKsXPvHnzFM9x+PBh3HvvvfD29kZMTAzmz59v1pc1a9agU6dO8Pb2Rrdu3fDjjz8q7hcEAbNnz0arVq3g4+ODIUOG4NSpU3Y8GkRERFRXt0orpMtJ7UIc2BMHB1XR0dGYN28eMjMzsX//fgwePBgjR47E0aNHARgDquHDh2PYsGHYu3cv9u3bh6lTp0KtVnZ7zpw5yMnJkX5eeOEF6T6dTodhw4YhNjYWmZmZePfdd/HGG2/gk08+kdrs2rULTzzxBCZMmICDBw9i1KhRGDVqFI4cOSK1mT9/Pj788EMsXboUe/bsgZ+fH5KTk1FaWtrIR4mIiIiskY9UeXo4eAJOcDJBQUHCsmXLBEEQhH79+gmvvvpqje1jY2OFBQsWWL1/8eLFQlBQkFBWVibdNnPmTKFjx47S9ccee0wYMWKE4nH9+vUTnn/+eUEQBMFgMAiRkZHCu+++K91fUFAgaDQa4auvvrL5vRUWFgoAhMLCQpsfQ0RERNZ9teeCEDtznZD45k+N9hq2nr+dJqdKr9dj1apVKCoqQlJSEvLz87Fnzx6Eh4fjrrvuQkREBO677z7s2LHD7LHz5s1DSEgIevXqhXfffReVlZXSfRkZGRgwYAC8vLyk25KTk3Hy5EncvHlTajNkyBDFcyYnJyMjIwMAcO7cOeTm5iraBAQEoF+/flIbS8rKyqDT6RQ/REREZD/nrhurqetKK2tp2fgcHlRlZ2fD398fGo0GkyZNwtq1a9GlSxecPXsWAPDGG29g4sSJ2LhxI3r37o37779fkcv04osvYtWqVdiyZQuef/55vPPOO5gxY4Z0f25uLiIiIhSvKV7Pzc2tsY38fvnjLLWxZO7cuQgICJB+YmJi6nRsiIiIqGbZlwsBAMG+XrW0bHyOXXsIoGPHjsjKykJhYSG++eYbjBs3Dtu2bYPBYAAAPP/88xg/fjwAoFevXti0aROWL1+OuXPnAgBSUlKk5+revTu8vLzw/PPPY+7cudBoNE3/hmRSU1MV/dPpdAysiIiI7Ojs1ap9/1o5dt8/wAlGqry8vJCQkIDExETMnTsXPXr0wMKFC9GqVSsAQJcuXRTtO3fujIsXL1p9vn79+qGyshLnz58HAERGRiIvL0/RRrweGRlZYxv5/fLHWWpjiUajkVY2ij9ERERkP21CfAEAXR1c+BNwgqDKlMFgQFlZGeLi4hAVFYWTJ08q7v/1118RGxtr9fFZWVlQq9UIDw8HACQlJeGXX35BRUX1ksv09HR07NgRQUFBUptNmzYpnic9PR1JSUkAgPj4eERGRira6HQ67NmzR2pDRERETW/vuRsAgPhQfwf3xMHTf6mpqXjggQfQpk0b3Lp1C19++SW2bt2KtLQ0qFQqvPzyy3j99dfRo0cP9OzZE59//jlOnDiBb775BoAxwXzPnj0YNGgQWrZsiYyMDEyfPh1//OMfpYDpySefxN///ndMmDABM2fOxJEjR7Bw4UIsWLBA6sdf/vIX3HfffXjvvfcwYsQIrFq1Cvv375fKLqhUKkybNg1vvfUW2rdvj/j4eLz22muIiorCqFGjmvy4ERERkZKfl2MLfwJwbEmFZ555RoiNjRW8vLyEsLAw4f777xd++km5JHLu3LlCdHS04OvrKyQlJQnbt2+X7svMzBT69esnBAQECN7e3kLnzp2Fd955RygtLVU8x6FDh4R77rlH0Gg0QuvWrYV58+aZ9eXrr78WOnToIHh5eQldu3YV1q9fr7jfYDAIr732mhARESFoNBrh/vvvF06ePFmn98uSCkRERPYVO3OdEDtznbD//I1Gew1bz98qQRAERwd2zYVOp0NAQAAKCwuZX0VERGQHcbPWAwB2zRqMqECfRnkNW8/fTpdTRURERGSL0gq9dFmlcmBHqjCoIiIiIpd0vahcuhyp9XZgT4wYVBEREZFLyi0skS6rnGCoikEVERERuSQxkAr09XRwT4wYVBEREbmozAs3EDdrPZLmbkJzXnem9WZQRURERA0wekkGACCnsBRpR/Nqae1+cgpKHd0FBQZVREREbmD32euO7kKTu1FsTFS/eKPYwT0xYlBFRETkoiK0GunybwUlNbR0TzlV77lDhOO3qAEYVBEREbmsPF2ZdLmorNKBPXGMllW5VIG+Xg7uiRGDKiIiIhdUqTcormc0w+m/f2w8AQBoE+zr4J4YMagiIiJyQdm/FSquRwc1zhYtruDC9SJHdwEAgyoiIiKXVFKuV1y/dKN55VTJt6iZMbyTA3tSjUEVERGRC7rgJCveHOWxf2VIl/vEBjmwJ9UYVBEREbkgH08PR3fBoQ5frp7+dIYtagAGVURERC5JzKnq1jpAui1f51zFMBuLfKXj6N7RDuyJEoMqIiIiFxRUtd9dYUmFdNv5681jSnDMJ9VTf+8+0t2BPVFiUEVEROTC7k4IkS4/9q8MPLJklwN70/j2nruBI7/ppOtqtXNM/QFAC0d3gIiIiGz3ytpsfLHnonTdoCxXhf0XbjZxj5rWqr3V733Fn+50YE/McaSKiIjIBew9dwND39+mCKgA4HpRGXrGBCpuk5cbcDe60up8qkGdwh3YE3MMqoiIiFzAY//KwKn822a3x4f6ISrQW3Fbp9c24oGF25uqa03q5+N5AJynirocgyoiIiInd+6a9YrhkQE+eGFwe7Pbj+foUFbpviNWD3SLdHQXzDCoIiIicnJpR3Ot3pcQ7m91ixpBaKweOYZ8WtOZSimIGFQRERE5uV9zb1m9z7uFGi29PfHhE73g7ak8rZdVGqw8yjVdkJWMaBfm78CeWMagioiIyMn9eCTH6n2tAoyjVL/vEYWn+sUq7rt2u6xR++VIHk5USkHEoIqIiMjJhfhprN4nn/oL9Ve2K6twr5GqCr3x/YT6ezm4J5YxqCIiInJiJeV6/FZQAgCYPqQDWnorS0zKi1+aBhu5upLG72AT2n32OgDg2u1yB/fEMgZVRERETuqzHefQefZG6Xrv2EBkv5EMazNffeKCFddL3WykSu0kGydbw6CKiIjISb257pji+l3tQgEAP00fgMfvjMFP0wco7o8P9cPch7tJ1zcdz2/8TjahFbvOAQCGdI5wcE8sY1BFRETkpExX84nJ2QnhLTFvdHd0iGhp9pgn+raBp4ex3fnr1utbuaLOkVoAwNErhQ7uiWUMqoiIiJyUfPru9Ye62Py4YV2MhTGjAi3Xr3JVYkX5yQPbObgnljGoIiIicnJP9WuD8XfH29w+MTYIAODcGUh1J1RVM71dVllLS8dwaFC1ZMkSdO/eHVqtFlqtFklJSdiwYYOiTUZGBgYPHgw/Pz9otVoMGDAAJSXVqxlu3LiBp556ClqtFoGBgZgwYQJu31bujXT48GHce++98Pb2RkxMDObPn2/WlzVr1qBTp07w9vZGt27d8OOPPyruFwQBs2fPRqtWreDj44MhQ4bg1KlTdjwaRERE1Sr11aNUE+6xPaCSu3ijuPZGLuR8VfHP+BA/B/fEMocGVdHR0Zg3bx4yMzOxf/9+DB48GCNHjsTRo0cBGAOq4cOHY9iwYdi7dy/27duHqVOnQq2u7vZTTz2Fo0ePIj09HevWrcMvv/yC5557Trpfp9Nh2LBhiI2NRWZmJt5991288cYb+OSTT6Q2u3btwhNPPIEJEybg4MGDGDVqFEaNGoUjR45IbebPn48PP/wQS5cuxZ49e+Dn54fk5GSUlpY2wZGqnSAIyC10jr4QEVHNzly9jZTVWSgoNi8NsOPUNcTNWo9n/71fui3E33qdKkuKqkZybhQ5Z+mBhnLGwp8AoBIE59oZKDg4GO+++y4mTJiA/v37Y+jQoXjzzTcttj1+/Di6dOmCffv2oU+fPgCAjRs34sEHH8Tly5cRFRWFJUuW4JVXXkFubi68vIz1O2bNmoXvvvsOJ06cAACMGTMGRUVFWLdunfTc/fv3R8+ePbF06VIIgoCoqCj89a9/xUsvvQQAKCwsREREBFauXInHH3/cpvem0+kQEBCAwsJCaLXaeh8jS9776SQ+2nwao3tH473Hetj1uYmIyL7iZq2XLp+fN8LqfaITbw6Ht6eHzc+/Zv8lvPzNYfRuE4j//vnu+nfUiVToDWj/inE26+eUAUgIN0/Sbyy2nr+dJqdKr9dj1apVKCoqQlJSEvLz87Fnzx6Eh4fjrrvuQkREBO677z7s2LFDekxGRgYCAwOlgAoAhgwZArVajT179khtBgwYIAVUAJCcnIyTJ0/i5s2bUpshQ4Yo+pOcnIyMjAwAwLlz55Cbm6toExAQgH79+kltLCkrK4NOp1P8NJaPNp8GAHx74HKjvQYREdnfyI934LGlGTXmCdUloAIArY9nQ7vldCr11WNA4tY8zsbhQVV2djb8/f2h0WgwadIkrF27Fl26dMHZs2cBAG+88QYmTpyIjRs3onfv3rj//vulXKbc3FyEh4crnq9FixYIDg5Gbm6u1CYiQlnPQrxeWxv5/fLHWWpjydy5cxEQECD9xMTE2H5gGuDDTcz1IiJyVqYTRIcuF2Lv+Ru44/U0lLvZBsj2dEM2Veqs038OD6o6duyIrKws7NmzB5MnT8a4ceNw7NgxGAzGD9bzzz+P8ePHo1evXliwYAE6duyI5cuXO7jXtklNTUVhYaH0c+nSpSZ53ffTf22S1yEiororqyFw+tvabLPbNC3qf6o+cLGg3o91NvL8s7qO3DWVFrU3aVxeXl5ISEgAACQmJmLfvn1YuHAhZs2aBQDo0kVZl6Nz5864ePEiACAyMhL5+cpqsZWVlbhx4wYiIyOlNnl5eYo24vXa2sjvF29r1aqVok3Pnj2tvjeNRgONpm7JhfaiK62A1tv9hn+JiNzZN5nmKRyt61FrSm8wjoYF+brPeeD7rCuO7kKtHD5SZcpgMKCsrAxxcXGIiorCyZMnFff/+uuviI2NBQAkJSWhoKAAmZmZ0v2bN2+GwWBAv379pDa//PILKioqpDbp6eno2LEjgoKCpDabNm1SvE56ejqSkpIAAPHx8YiMjFS00el02LNnj9TG0QJM5s+v3SpzUE+IiKgmdV2pffZa3auitwn2BQB4NWCUy9nctLBS0tk49Ginpqbil19+wfnz55GdnY3U1FRs3boVTz31FFQqFV5++WV8+OGH+Oabb3D69Gm89tprOHHiBCZMmADAOGo1fPhwTJw4EXv37sXOnTsxdepUPP7444iKigIAPPnkk/Dy8sKECRNw9OhRrF69GgsXLkRKSorUj7/85S/YuHEj3nvvPZw4cQJvvPEG9u/fj6lTpwIAVCoVpk2bhrfeegvff/89srOzMXbsWERFRWHUqFFNftwsEX+BRIUlFVZaEhGRLfafv4H0Y3m1N6yDDdk5GPjPrWa317RPcKsA7zq/jpPvO1wv4nnt/k7htbR0HIdO/+Xn52Ps2LHIyclBQEAAunfvjrS0NAwdOhQAMG3aNJSWlmL69Om4ceMGevTogfT0dLRrV12e/osvvsDUqVNx//33Q61WY/To0fjwww+l+wMCAvDTTz9hypQpSExMRGhoKGbPnq2oZXXXXXfhyy+/xKuvvoq//e1vaN++Pb777jvccccdUpsZM2agqKgIzz33HAoKCnDPPfdg48aN8Pau+4e9MZjm7N0qdc5qs0RErkBvEPDIUuPq7i0vDUR8aMOLTQqCgMlfHDC7/Y7WWswf3QMPfrhdcfvaP9+F1fsuYdYDnRr82u5AzEVz5pWNDg2qPvvss1rbzJo1S8qvsiQ4OBhffvlljc/RvXt3bN++vcY2jz76KB599FGr96tUKsyZMwdz5sypucNO4uzV2xjQIczR3SAickliThIAnMy9ZZegyjRB3UOtQsrQDnj8zhiz4p7/N6EverUJQq82QQ1+XXdxKs+4W4o9/i8ai8MT1alxaJx0ZQQRkau5VWqfdIrT+dVbqM19uBsevzMGKtk83dyHuyH1v9kIa6nBve3t86U4T+c++bW/FRi3qHOukuVKDKrcTAu1CpUGAb/m3XJ0V4iIXNbV29XByGc7zuHRPg2vMygPqp7o28bs/if6tsHQLhEItMP0VqBvdcHrCr0Bnh6unbAur+3VLtx5R6pc+yiTmcqqIWt3WvFBRNTU5EU4T+TeQp6u4Xur7jl3vdY2of4atLBDAOTv5V5jJuWyzaX7xAY7sCc1c6+jTogO8sHlmyXwcMelH0REDnI8R4cIbd0XJuXpSjFu+V4M6BCGED9j3lR0kHNuseLMjl2p3uatpbfzhi7O2zOqFzGW2nmm9m9ERERkG3niel288OVBnMi9hRO5t6S/z0GyqbmmUKkX4Opptn9YvEu67Kdx3tCFc0RupqDYmFAZH+JbS0siIrLGdH++c/UowAkAe8/fkD2n8d+4Jli9Jh/NybXD1KWjPdgt0tFdsAmDKjfTPToAAKBjnSoionq7dLNEcT2njlXQAaCsUm/x9uFdGz9AUKtV8PUyDk+5QzLIj9m5AIBXR3R2cE9qxqDKzYjJlYcvFzi2I0RELsy0oPJnO86htMJykGTNuOV7Ld4eE9w0OVXumFtbVFa3/4OmxqDKzfhWrfiIC3HeJadEpias3IeHPtph9Zs9kTPIvHDT7Da9QcCG7BxcKSjBfw9cxl9WHYShKv9q99kbZu0BoFUAE9XrorC4uk7Yw71bO7AntXPebC+ql8iq1SmHOFJFLqKwpAKbTuQDAA5fLsSdcc67XJqat33nb8Db0wPnrhVhdO/WUKlUeP7/9uPn4/mKdv/LulLj83g1Uc0osYL7lcKSJsnjaiyXC4qlyzHBzp0vzKDKzYi1PCr0AqZ8cQALxvRkzSpyatdlRRbFHBAiR7O02u+Dn0/hg59PAQCu3S7D/Z3CzQIqWwT4Ns3edeL5wGCopaGTu3673NFdsBnPtm5GvifS+uwcZF0qcFxniOro4vXi2hsRNQFxn7k4Kyup5204gaELfrHpub5+Pslu/aqLTpEtHfK69nbhhvHvQqwLrGpnUOVm2oYph3jPXL1tpSWR8/k1j59Xcg6+GuOo6a3SSoxLim3Qc/WNd+yU9vUi197/78jlQgDABRf40sWgys2YTp+UV7r4uC+5vRtF1UP7fhpO/5FzOHixAADQLToAI3vVnhydGBuElt4tsOmv9+HQ7GGK2+W6tNLatZ81OXvVWFtLV2KfDaEd5dJNYzDVq02gYztiAwZVbiispUa6LC/tT+SMSmTL1C2triJyhKu3jKM7N4srbMr1mzOyK7LfSEa7MH8E+HpKZRNM6yody2m6v8k9YwIBwOU3UxZnXMRtfpwZE9Xd0CdPJ0ol/ZvyF5iooaICudScnIN/1VYoES01tX4un70nHl2jAhS3/fLyIJRVGuBtsj9MgE/TJKkDgNbHPU7xeTpjgNsx0t/BPamda4evZFGvNkF4qEcUACD7t0IH94bIdttPXXV0F4gAAHvOGWtMJcYGwc+r5uDk1d91MbtNpVIpAqrH+kQDABY+3tN+nWwGTudX51maBq7OiEGVm7o3IVS6bLqHFZGzinfhWjrkvA5fLsDJ3FsW77t0oxi/5pnfd62q1Mftskp4mJZXl5l4b7xNfZj/SA/8+tYDGNgx3Kb29lRU7rpFdad+eUC6HKF1/uk/BlVuamSvKOnylXrsWUXUVG7L9qlk/E/2dulGMX7/8U4kf/ALisuVe6KWVxpw/3vbMGzBL2YbJrcKMBZS7t3GmGj+7D3VwdPZdx7EgjE98NrvumDG8E4296WpawbeLjO+X/loj6s5IQuGE8Kcv0SEe0y4khlNi+ph50s3itGauSrkpORB/0/H8rD24GVEB/mysjrZxRFZCsSF68XoLFt9V1hSIRXIvHijWDFSKm6gLC78efV3XfBIn2hEBfpArVbhD72im6L7DaKq2ko51N/LwT2xj6YqmtoQHKlyY+Jemmo33FST3Ifpyqrpqw/h0aUZ0v5pRA2x4/Q16fL01VmK+/afr96br0Q2RZYrC/S9PatPk50itdB6O/+JXdTRDYp/ijW+bJ1mdTQGVW5M3F+qyGTIm8gVVLj63hq12PbrVXzyyxnmPDYycT9UwDiVVFBcrrguyi0skS73n7tJuhzq7/x5PLU5dNl1FyztrVow0KOqPISzY1DlxsTNNN//6VcH94QAYMoXBxA3az0uXC+qvTFh/eEcR3ehUY1bvhfv/HhC2kuOGsd76cq/f/Kq3C29qzNgfKpGTOXFaAHXrvF0+aYxUGzKMg72Jk5dusrAtet+WshmLKvgeKUVeqzPNgYJi7acdnBvXIMrbaJaV/ITd2mF667McnY3i8w/Q/KEdHEUBIC0T2qObMQKAPw0rpt63DfemGTfoobVi87uWtXfAVfJC2ZQRdTIBEHAgPlbpOsXbzj//lXO4N2fTjq6C43GIJvyO25lqT81XLGFgFW+Aq91UPWJWsypyilwn9XS7pBPK74FVyinAHD1H1Gj+/bAb8i/Vb2haVwIazHZoltr5y/0Vx9Lt52RqnUDwLErHEluLCdkO0rEh/rh3LUixUiVfG9UsfxArs59giqRwYXz9sSuy1e0OzMGVUSNbN6GE4rru85cd1BPnNM1WcApd7PY8tTNydxbGNYlAmoXnNLYfCLP7PPgW0u1bqq/LSfzpcvilKt89GaPbPpPXNUnr+n0cO/aN1J2ZmIwtflEfi0tnZN8mlzj6RoTa67RSyIXJlZmFkUHuUZuQFO5ettyUHX2qnlC/+D3tmLSfzIVJ0tXYilPjNPBjadVgPF3LdDXE4M6hgFQ5rDJ61JtrvpMVVTVrQr288J7j/Zoqq42isqq7O5OLlpaIf9W9aihq5SyYFBFZGcVegMe/yQDL3510OL9HKlSkq/AkpMvhQeMuWniVMDX+y81drcaBUelmtbBizcBAMldIqUin8dlU4Lpx/Kkyx3CjYHHF3suAgD6tw2GysVzkuJdPNUg1wV3A2FQRWRnJ3JuYffZG/j+0BWzbTEA7m9nK9PcFvmS6oRw59+t3pJz18y3C6lpX7mmIAiC2xZaFauhX7pZDH3VexRrU4kjUqK9VYVAxTjKHepTuboPN7leuREGVWQ3lXoDnv5sD8Yu34tKvfXCjVdvleHpz/bg3xnnm65zTei8rA6VrqQ6qHok0bithekeY2Sdtc+Rq65qsjRSpXdgQKMrrUB86o9o+7cfUVhc4bB+NJaDFwsAAP3bhiA6yBeAcUrJYBDQ/pUNZu3lo6HJXSObqptkxYGq/z9X4tCgasmSJejevTu0Wi20Wi2SkpKwYUP1B33gwIFQqVSKn0mTJimew/R+lUqFVatWKdps3boVvXv3hkajQUJCAlauXGnWl0WLFiEuLg7e3t7o168f9u7dq7i/tLQUU6ZMQUhICPz9/TF69Gjk5eWZPY+zEle2NKa9529g+6lr+OXXqzh/3XqeyLOf78P2U9cw+39HXbZGT03f7r09q1epyPNl/tg/ttH75W50pZY/t676uSmx0m/TUZOmcv9726TLhy4X1Pt5BEGo8YuUowT5GgtHXrtdhpiqXMbWgT64bqF+lUiM18UNlclxgv2M/39/uivOsR2pA4cGVdHR0Zg3bx4yMzOxf/9+DB48GCNHjsTRo0elNhMnTkROTo70M3/+fLPnWbFihaLNqFGjpPvOnTuHESNGYNCgQcjKysK0adPw7LPPIi0tTWqzevVqpKSk4PXXX8eBAwfQo0cPJCcnIz+/Ohl2+vTp+OGHH7BmzRps27YNV65cwcMPP9w4B6YRXKwhyLEX+d5ZRTUEcUF+1Zt7uuLu6QaDgJGLdlr9di8vtCffgiRGlqCuK3W/UYHGIF+1pCupPmYnXLC2U3mlAe+mWa69JS5mWLP/Eib/J7PJRq+uylZemlYSr4sXV2Uh4ZUNir30nIFnVU2qbq0D0CHCmDN15moRdp25ZrF9rq5UGqnyt5Lr54r2nb/p6C7Ui1jy4sFurRzcE9s5NKh66KGH8OCDD6J9+/bo0KED3n77bfj7+2P37t1SG19fX0RGRko/Wq3W7HkCAwMVbby9q79hLF26FPHx8XjvvffQuXNnTJ06FY888ggWLFggtXn//fcxceJEjB8/Hl26dMHSpUvh6+uL5cuXAwAKCwvx2Wef4f3338fgwYORmJiIFStWYNeuXYq+miorK4NOp1P8OMr3h640+mtsP1X9h+qqlWXygHFneFF9ZnHOXSvCa98dwak8x5xYbxaX43DVXlrLdpytsa08LyhYFkzKA9Dm5Ny1Ijy2NAMjPtyOdYevYPHW01i05YyijTzH6KU1h6TL8kC0XZjr5VRN+Hyf1fsq9caRnpe/OYwNR3LxY3bTbNHTNqw6v68htYx+qPr78sjSjAb3yZ7EOlUeahXk7+7Yleq/xfd1CJMu7z5bvYjEVVab1cTVR93EGRZXqp7iNDlVer0eq1atQlFREZKSkqTbv/jiC4SGhuKOO+5AamoqiovNR1ymTJmC0NBQ9O3bF8uXL1eMDmRkZGDIkCGK9snJycjIMP7yl5eXIzMzU9FGrVZjyJAhUpvMzExUVFQo2nTq1Alt2rSR2lgyd+5cBAQESD8xMTF1PCr289Xei436/JtP5GHlrvPS9V/zzQOe0go90o7mKqbNLtVjOfnEf+/H/+2+gD9+tqdefW2oK7KKy7VNq8pHCOUriX4rKLHU3O099elu7D1/A0ev6DD1y4OYv9F85MbaKM0WF621I5J/6TD1W0GJovp3TV9K7KmlLHCw5QvOoi2n8eGmU7hwvQgLfz6FuFnr8fDinY3Yw4YRSwqUVOilbU78vDwUX+zmje4mXb54vfr3Uj6N76paBxrzyFwoJlEQv2BFulBQ6PDxzezsbCQlJaG0tBT+/v5Yu3YtunTpAgB48sknERsbi6ioKBw+fBgzZ87EyZMn8d///ld6/Jw5czB48GD4+vrip59+wp///Gfcvn0bL774IgAgNzcXERERiteMiIiATqdDSUkJbt68Cb1eb7HNiRMnpOfw8vJCYGCgWZvc3Fyr7y01NRUpKSnSdZ1O57DASv5HxJ70BgFPfrpbUUQPAPIsLIX99JezZpub1mf6T3xMns544jEYBPxnzwWE+WvwQBMME+fJRp8EATh0qQAdI1tKf4TPyhLRxfjA9JuiCxc4bpArNi6RHn93HFbsPK+4Tf4ZM/28OTuhlv/wX/NuKUpI1DR9bk+HZXlUJ3Nr/l28fLNYmr58X/Z77MzJxIG+nrhRVI7YYD8paCwq10sbDP+ueyu0CvCBn5cHisr1+OXUVQBAbIivo7pMFrjSptYOD6o6duyIrKwsFBYW4ptvvsG4ceOwbds2dOnSBc8995zUrlu3bmjVqhXuv/9+nDlzBu3atQMAvPbaa1KbXr16oaioCO+++64UVDmSRqOBRuPey3LPXy+yeIIzHYkpq9SbBVQAcO5aw3O9Dl66idn/M+bhZb8xTPHtuzHIk4pX7jqPlbvO48FukVj8VCIAwMuj+nvhmavGE1VZVW5AbIgvLjRBfltjEQShQbV7vDzUUr0gUwE+nigsqcDCx3vCy0ONFTvPI0oWjMrre0W50DdXQJmg/kTfGHy1V1lna92hHAxoXz0NVdNCD3uSx3q+XjWPzIhfYmqTf6sU4S2d4/9HfH9eLdSKbU7E4rFRVaNXRVXT8ZkXjLlHHi66utSdCILg0JWx9eXw8M/LywsJCQlITEzE3Llz0aNHDyxcuNBi2379+gEATp8+bfX5+vXrh8uXL6OszPgHIDIy0myVXl5eHrRaLXx8fBAaGgoPDw+LbSIjI6XnKC8vR0FBgdU2zdVvNy1PY1XKfhlW77uIjq9utNju2wOXzW47lXcLGXUokCkP6n441Pi5KDcsbJ/yY7ZsxFL2B1ksr2CaBFwge46nP9uD7m+kKYoSOqNf827hzrd/xnsmGx1bqsVljbWACgBeGJyA7DeGYWTP1tJGt/IArqMNVaF/KyjBgPlb8J/dF2zuU2PTGwRcu1X9/z1reGezNnGhvoqpZGsFUe3JdPRsXy1J5ou3WP+7K/fF7sZNNagLsXyJIAhSzSqgeno138o+f2WVzreSsbmR75cqjiy6AocHVaYMBoMUEJnKysoCALRqZX2KJysrC0FBQdIIUVJSEjZt2qRok56eLuVteXl5ITExUdHGYDBg06ZNUpvExER4enoq2pw8eRIXL15U5H81R79aSRbfevKqdHnmt9k1PofpH/ehC37BE5/utpqIHuir/AU7Kks6feP7o4r7BEHA6n0XsdWO25pY2moEUE6liMST451xQQAgjVIVyFYNbj91DbrSSkz6T2atr33pRrEiILOna7fLFBvMmvrkl7O4drscH22uPrmOWrQTXWanYdn2mhP2RWJeizWmo4zWkqc3ncjHGgtV1f/232xcvFGMV787YlN/mkL/uZsw4N0t0nXPFuajIJkXbioC7x2nredf2YvpF5fa/m822ZjTdulm442ylVcaaixU+sWeCxi/Yq8ULHlVrf4LMPmbcbPq908sIjvrgU6K+90t59HWaXdnIv/Vd6X8NocGVampqfjll19w/vx5ZGdnIzU1FVu3bsVTTz2FM2fO4M0330RmZibOnz+P77//HmPHjsWAAQPQvXt3AMAPP/yAZcuW4ciRIzh9+jSWLFmCd955By+88IL0GpMmTcLZs2cxY8YMnDhxAosXL8bXX3+N6dOnS21SUlLw6aef4vPPP8fx48cxefJkFBUVYfz48QCAgIAATJgwASkpKdiyZQsyMzMxfvx4JCUloX///k170Brg7NX6lS/Ye+4G1h2+YjEv5K31xxXXH60qcOnpYfvw+aUb1X/Arsv2gfvzFwcstpcHJGWVehy8UL1c2HQk5MhvOsz8Nht/WrHPbrWArM3vv/PjcTz00Q68vf6YdJtpnkq31gEAgBZVx6essnpaqLiWFYFzfjiGe+dvQc856XYfFt915hr6vPUzOrxqXhBRdOCC+bLsrEsFAMw/B9bYerISP2o5haXS58709V/+5rDZ4+SFV51Bpd5glnTu6aFG2rQBeLp/LEb2jAJgrPwtL2jatpGr7guCgCeXKRd6XKv6slBWqcfX+y7V+++F+JmwRm8QcDr/dq15ZqYuXC9Crzk/4Q+Ld1p97Ctrj2DLyatYZDKqZi0lQAyq+sYH16kvrkI+4ulqVfPFDdVbuNLSPzg4pyo/Px9jx45FTk4OAgIC0L17d6SlpWHo0KG4dOkSfv75Z3zwwQcoKipCTEwMRo8ejVdffVV6vKenJxYtWoTp06dDEAQkJCRI5RFE8fHxWL9+PaZPn46FCxciOjoay5YtQ3JystRmzJgxuHr1KmbPno3c3Fz07NkTGzduVCSvL1iwAGq1GqNHj0ZZWRmSk5OxePHipjlQdvLV3ot4ZUSXOj2mvNKAx/5lXOEYMckbPaIDce/8zfD29MAPL9yjaLtz1mDcLCrHmszLqNDb/gv8XdZvePH+9gCUw+6nbEhiv3a7HBrZt5j2JtuXFMmmpvJ0pVJVZbmfjubCs4UagzqG29TfAxct13zZfdZ8+kSsPyQGjuIfOfH4GGRxXtco83Ihcst3npMul1bo4aex36/vP2Sr8Eor9Ba/GXaIaKlIwrc3+WKKEH8vs/tbereQRhhE5ZUGaTQCACJaejtVztolC9Pjnh5qdIxsiTdH3YG1By/jf1lX4OmhVmy8bbpFjz3N+vYwVu0zH+X7+XgeNmTnYHLVl5kgX08cnD3M4nN89EQvrDt8BU/1i8W8DSdwLEeHYV0i8NOxPIsbYcv1eSsdN4srMGVQO7yc3KnGtnLbT11DUbkehy4XorTCgA82/YoWapXF5xBHr+Ujr5bOzUM6G//GR9cySueqXGnazJT4O1DpYsGgQ4Oqzz77zOp9MTEx2LZtm9X7AWD48OEYPnx4ra8zcOBAHDxoeXNb0dSpUzF16lSr93t7e2PRokVYtGhRra/nrOSlAGpz5LdC/HD4CrpGBUi3nc6/DU0LtZSwmnakOo/ozVF3oHWgj+KPWFFZJU7aUEtKHqTUVivHUrFN+bYvpoHYZdlJraC4AtFBysdeuF6E5/7POO229aWBiLNhhKC2aRJLxHwgMWFZrJ8jyKrnWJtWtOTijWJ0blVzEFYXvrIg6kZRuZTAK1q97yI2Hq3+/958Ig+DOylXzNamthWo8qDS0nYuvl7mQVVppV4RVAX5VZ9E5AnTgiDg6BUdOkW2RIsmXElUW+V3VdVi9+2nrilKLhy+XNjgRQHWmAZUfeODsbcqL3GybHRYfqxvyqYm9786BKH+GjzUwzjKNqCqztOy7WfxU9UGxbrSCqt1nsTnXbTlTJ2CKnkfdp25hn9tM045T7inLYL9vBTH+mTeLcV1j6rdNsJbaqRcnUBfT+mz4GvHLyhkH5eryu04em/MunK6nCpqPHWpffP8/2XiX9vO4sWvqoPR1P9mY+HP1RtcyguKitNascHVI0F5ulI8vHiX1dcY1NH4xzi7qpAmoEx89/Y0/3jeKqtbaQj5c1hapv7tgd+ky9byw0ztv1D35fxi4U+x6KeY7C4PdGvLRfGxsv2NPcjP3ZYSz03z4p79fL9ZG3lBRUs2n6h5W6eisuqToEYWKF29bf1ze+S3QsV1cVsSAOj79iZpunH8yn343Uc7MLqJi1PmmYw4meYD1jSN+8Nh+y+6sPR6v68KjmpSKpumtrbRcGJs9TeWIe/V/IW4Pjxln4m1B6t/bysNBpRV6tHpterFMP3iQ6TpIwCI0Jr3+fE720iX/RlUOR1xZbirFS5lUNWM6OuQw2At90WerJotO6GJQZVarZKmuHJMkiODTE4oPWOMf4TloyJq2beS0gqDIndi99nr+P3HykKDVyz009Z9Dk/n31bsgq6xMRkyrAG715uekOSrjwpq2NBWbxAUy/Lz7VwcMlOWr2T6/2aJpVhg6bYz5jfKTF9dXR39D71am90vP2nKpzbza1jK/+Sn1XlB8zacMBuFef7/jMGfuHDi0KWCOufyNITp6LBpYF9TPaSjVwqt3ldflvYe7FdDPpE42pNrw2einWzq3drnU/7+LQU6NZEPWKyTBZwGA7D5uDKJ/ufjeXhXNqUtjvjJ+yX+zSLnJM5ayL8ouQIGVc2I6bfmhhKDAB9PD8UQ7a2qTXBNRy52zhosfSMM8vVEl6rpnuzfrJ885BufPv7JbrPSBJa2e8mRBVryoEtevkFvEDDkfeW3abGo6PlrRYibtR4D5m+BJQ0JaO6o+kN+9DfjsSkzSZ63dvKav/GE4voXdigZIAgCvj90BYcuFShy2UxP/NaCENOVgnUZCbWU62EtObu2ArFiPywFdSct7BFYU1kHezOdvft8fF+bH1vTSkx7igm2HtiJeV5iEn2ohVw3kel0n6VRsZGLqr8U2Vr3qja/FZRYnCb9b9VolrVRqOgg5RT3iTdrTyVxNS1kC4ZuNtKq4cYi/u02/X9ydgyqmpHaVpeJ6rqyzPTbrzhca3q7r1cLfDflLswY3hFbXhpoMX/quslUj9hklZVtdr7L+s3sNvGk+Z/dF/DOj9XByNf7L0vlCCytBBSnCsWqyhdvFJv1B6jeWqdnTKB0W20JoeLJSQxYQvy9UFxeifErlPvB9Z+7yWL5h70mNYTsESDvOXcDL351UHGiA4BjOcpAxDSPSWQ6TZhx9joWpP+K6auzcKGWVXjyPeeSu0bglQc7W03UF6dlrY2e1hQk3ZMQCsC2P8zbfr2K/1qom9YQplPKwTUEJaYq67DYw1bygPmtUXdg3Qv31LhcfWdVaQfxBCcvoGnJqueqV0NvOGI+fSkPkH08PVBQXF7ryGFphR75t0qtBtdqVc2rjeVTrpv/eh8AID7UD+0jlItavD09FPl57kC+6lE+ve4KxNFq0ylzZ+denyCySAxybN2FvqE1ncThWktTDQnhLfHngQkI9PVCx4jqYo5ikCOOcokeWboLpRV6zPqv5VpX/5XlRIVXFfcTl4ZbqlV0qCp/y1LS9IELBQCUBSd3nL6GknI9Mi/clIJNXVUfJ9wTj/mPdMe3k5Ow52/348SbwxWbJj8sm+IS31f3aONIVZ6uFF/usRwo/mnFPpw3WWXnVZVQK548bhZXNHga630LFe4B5UhfQXE5PrOyabSlIH3hplNYe/A3/PXrQxYeUS1BNlU0/I5ITBzQ1my0Ia5qamynSc2muxNCFCfR3MLayzTcrOWzX1qhx7jle5Hy9SHsr6UIZl2YBiztw5UFTGv6AnPwkuVVpnKHLxcoVg1aU1qhR8rXWfjLquocyT/2j5VGTq25UWT8PRG7WVtJjP5tQ6TL3x1UfuEx/byWVOjRc0464lN/rPFvzuOf7Ea/dzbhfA27L/yaZ300U75YpW2YP86+8yC2vDTQ4mIId1RbpXxnJS6eCHOS6vy2ah6fKjf31rpjUrBgSV2SMHMLSzHBQhJyTeSjDkD1MPPFWpa2yzfJFMsdtDD5xnnherHV4EOufbi/tPKvphVXG4/k4L4OYWbBG1D9jeiCLKBZvOUM/rIqC4CxgOeaSXdJ9wX4eEoroEQdIvyl0grvj+kpTUGIxFGVKwWl2H/e+klz4D+34vupd6N7dCAAQFs1EtY60EfawuSFrw7i4yd7W32O2uy1sn9eiCzva8j7v1g9addUE2q/SU0p+Qn1y4n9FInopoGGKCbYF+evF5tNK374eC+E+GsQN2s9AONJM8HKc2ypyqUqkgWAZZUGsxEXebDwyNIMu307lk9PLxjTw2wlk+kqS7nacknGr9grvb/z80bU2DbtaK7iC0hN03hyx6pWqYqfAVvqOfWNC8be8zfMAsqavtTN/fEEBlopaSLWvTIdrRVdvFFsVoG+daCP9H9q+l7VNawmc611Zu4trKUGvxWUWCyF4cw4UuUGxN3tvTzUaBtaPQLwwB2RaKFWYfrQDtJttY1upB+veYUWYFyZ1Um2ZYhpXRrxj+mp/Oqpj1dHmG/NIf+ju/XkVRy8eNPit9E5646Z3WZqw1/ulabjrFViB6rLIcgDBfFEJwYC8hOqvCTEvvM3Uak3oGVVkGppSun9x3oiqW0I3hp1B4DqSurD7zBuZySOGoZrNYgPq7l8w77zN5F1qQDv/Hgc6VVL1Uf2rB79Wnc4p1H2xtpxuroafk2jINk1BPKm7ny7ejeCnjGB6CAbpTQNykViOY9btSw8qOkIdIxoabZwwVLu1yWT1ZQFxRV2+RFz1eY93A1/6BVt9rqmXyIAoEfVaGamhWKrcltkuxbUVrrB9AvONRvLd4hB4Q6TacCaDOhgnHI1HWGsqaSGLaVXrFGrVIppzS0vDVQEWWItKlv8sX8sAOPfTnIsMZiu6YuHM+JIlRuYNqQ9nrk7HqH+XgiUfbtd9GRvFJZUKIKXG0XlipEIU55WvhacnzcCZZV6bDySi37xIXjvp5M4UZUEHN5S+Xxh/hqcu1YkDcmrVcCz97at8T38L+s3m7YW6RcfjD3nbiA+1E9Rn6qFh1qqAF3TSVg8UckDwcf6xOCrvRellXhFNeSe3Tt/i/T8lpJjowJ98JUsr+TLif0VRUfFKYcL14ulJOo/3RWHtKO5Zqvu3rQQTJoWxdxx+hru6xBm1q42NQXX8VWBuaV90eTHfe6GE2b3W7L91FVFcObloUYLLzU2TrsXKqisTsOIU4SWRhUB4+qt7N8KUVhcoUjqvichVAoCKgwGDDNZkGApAJGPVP2ccp9N78tWvl4eNp8Yds0ajO2nruLQ5exacyC7RwfgcFVge/VWWY0J56bbtNg6UiUGK+LvuC01g8R+W8vFs8S0aG/1c9m2klcs9nlv+1DEh/pJf5sA4M4426ulvzqiM/50V1y9atE5O11p3crROAtXK2DKoMoNhLf0hqXZD7VahSA/rzrl3hyxsIxbLDKpaeEhjZQ83jcGazKNSb3yPCIA8NMoh/1rmjJIahuCjLPXsc9kKqxrlFaxp59o9fNJEAQBo5fswjmT+7pHB2LHaWNRwNQHzEfGACBCaxwpql7VBLSpOhnlFJZiy4l8i68rkgc+piUiLPH0UCuquIvHystDLQ1rX71VhkGdwm2a5hybFIen+sWi3d9+BGA+wmIPm6pGKzeb7PXWLz4Yq59PkqbdRD6eHhbz50SHTUa0xGC0U2TNxUtjZCOBlgI8MRC6eqtMMQry72f6YvfZ63hy2R6L1b3zdWXoGmUsIPnWuuPoExekKKSaYOUE3xi8TaYhowJ90CrA+L5rm4KUH9fayoiY/h+Ir1EbcQRJDKZG2lDTSl6vSl4EtKbiw0F+loM8W/50HcvRSTXzxJHbAB9P6TNRl5EqlUpVY3DqisQgV1dL8V1nFehiQRWn/5qZ2kr+i3WU5HlYp/PNh+YTY4PxfxP6YlxSLNZMUm4qbTqCU1NulWBl8ibAxxNPVw3Fm1KpVFKFcqC6zo58RZ+15ehi15ZsNS6/79UmCP6yqYL/HvwNwTbm09SnYKC4oqVcb8DPVbV1urbW4o2HumLqoIQaH/vuI8Y9L+WjBep6VNyu0BsUhVtF4uiFeGx9ZAmuJ94cjtXPG/+fx/SJUTwuJtgHZ9950Oz5xFEG0+k2Wyskh8lGQC1tWSTWePKW9VOlMn6Z0FgoHCt//dIKPZ78dA+O5ejw74wL2FFVzXxIZ9u2KrIXeTAhTlmJU8SmdctKyvVIfDMdcbPW4+L1YsUIsWnbSpMVkd9kKlc1Wio1YYmYB7X7bNXmyzZ83uS7MHR/4yfp8rrD5p850RkbtqSyRvxdBqr/fml9qn835ZebI+lLgovlJolas6QCObPaCjseqprHlk8p+VhZcn1v+zD8feQdZpuVDu2i/GbYJcr6CiP5t1o5gyDgFQt5WKKbRdUnkXvbG3M45DvNm27GLJ7HxW/s4ojRrdIKaGVBVVSAt8UEbvlScVF9tjuxNL2qaWFcyv1ScscaRydCZSdR8eRv6/SI3Oe7zkvJ93KPVQVLBy8WAKg+Ud+TEKqYQg5tqRxVGN41Emq1Cm+O7Kq4XVx1JR9F2j5jkM39bBtWPWJkaXPfiKpVQSUWjkGYv/UVQ1/vv6Sovg1U7zN22cJefU1FnOaUb9gt/3IwctEOqW7b6v0XFfXScnXV/X5pzSEkvLLBLJCSMw1sP3k6UZGHZDoKK+6naMtKQ9Mpamla3so0rnhf/q1Ssy9Ddc0YFP/2/GdCP3SKbImZwzs1ylY/rsTVEr0B5WpdV5v+Y1DVDKhU1VXOa/sFE1ccyQvFPXBHqzq9XlxI7fvniR5JjLF4e4VegLenB/a9MgTP3hMPAPjH6G7S/fIkcTHHSV43avfZ64rdzU2XhIvVnJ+9py3ayU7eubpSRWVokXypeEOEa81P9oWyY/3fyXfhlQctB5PyPA+x5syZq7dx97zNiJu1Hu+m2Zbj9Nb64xZv9zQJEsVimrvOKBOOTYs8iiNKT/RtoxhdFHOXxFGRRxKj6z21knXJfFq6WHr+22YJ+5b+EIv/5+ssbP8i5tp1qWVT68YkLoCIkH1GOry6AYu3ngagLBtg+kVno2wfTjGYemmNsazFxH+br+Y1XS03rGskst9IxtaXBuKfj/aQRiUB43SwGIRZy32SM/0ciVsaiYstRnQz/3tSrjeg79ub8OhS5bZWOTYkxsuJdaZiQ/ywcdoATB7Yrk6PJ+cg34bL9Eu7s2NQ1UyIJ8Laql6Ly43lQYS6jp+SpHYhGNmzOvfiWA3bbQRbWTYullsIa6nBq7/rgvPzRmCMbK+urq2rT35iTpRKpcLEe40B2O2ySsVUp/jNWxx9Eau4B/h64o7WAdJoVZtgX8XGtnIdIhon16ZDpHwlnD8mDmiL76ferWgzY3hHxQlN/PL91d5LUqC4aEvN28TUpmebQOmy3iBIo5qmW+uYLnSIrQqiW3io8eaoO6Tpq0s3jP0S/39syUEzJW4lIq+GLwb+Yq6guPTaeJuxjb9J0DB9SAcE+9W+LUp9plPtRfxyYLrn5fyNJ3E8R5nnZ7p6L+2oMWApq1Tmtp29elsKZgDggzE9MaxLBJb/6U6LfYgL9cMjidFIkH3R+MfGE1JejvwLSE3kXwx+KyjB94euSOVETJPm5UxLw9R1cSv38HMPlvZpdRUMqpoJ8aRTWlHz1hdicrJ8qfe4u+Lq/HoLH+8lXTadDpSzlu/wz0d61Pj84gkbAEbJygxYG/l6dUQXxfUWVZGimIcl7kWXcea6FASYkgeazw2oeTVjTUz3vbNUj6h7dCAeSTQuwZ8zsiv+PDBBMY1hbQXn/e9trXFKsKZaQX1kU7GXZZs7v3h/e0W7YD/lSdE0EBGDMTF5XSzT0aYOI5gi0+lQTQu1FPiLxWPPXL0tBZniCI7piKxKZTk30NSAeqyktBcxaLY0XbVq70VFwVMpx0mmtEKPVJMiuYNNNjYe1as1Phnbp9ain2q1Shrtk29ybuuWIRMHtJXKrpRVGhQbs8t/Xy251YBVaq62/L4pWdrSy1ldlir4u16I4no9pnoR6yXV9kVcPOmWlutx8q3h2PO3+2tdpVWbmqZUVCoV/GSJxkue6o3z80YokqQt+V134xRCqL9GsfluOws1j3bNGoyBHatPlrfLKqVgLjbY2P5qVa7IydxbiqFnwFjMEAD+PDABgzuFY2iXCMwc3gn1tWBMT0WdL2tTKv98tAfOvvMgxibFmd0XZWXn9jNXi/Bu2kmL9wHGIpByYq5YctcIxTD7fe9ulS7f1U459Wm6K0xUoLIv4ns7d804XdW5lfG6afK0LeR5coByMYKY31OpN+Bc1So/cYRKpVJheNfqWkPRQT4YZKW4pJxXDdudNJbNf70PY5NisWxcH+m2E28OV0yTZV0uRIVs2xpLqy03HlEW+DQlLnSw1eN9jV9Q5CVGrH3hsGT83XEWb+8bH4ylf0zEkM7h2PfKEPzzUeUXKPmq00pD3T4zpvlcVL1XaW3V8J2JmMMZYmVVqDNjUNVMiN9+a/sWKOZ1dInSQtPCQ5HfUVeb/nof3n2kOx7ubV70UE7+R9tSzpElbcP8seWlgfhp+gDF7ZaSx8NaahQjOxuyq3NqxKlNMXC0VONKPFFHBnhj+Z/uxKdj+9i8gs2alVUb68aG+Nb4nq1Vf46zsvkwAKzYed7qffItK+Y/0h3924bg/LwR+NfTfaw+pq3JlE9MsHI0wLTOlHhsTPOaauqzNWEmNdDk00HiqiY/TQtptZ98els+Cmrpc7znb/dLK0erNX1Q1TbMH3NG3iFNowLGk8qip3pLgZW4gER0wcKK2mmrs2p8nWGyINMWvWLMF5HUJel7VC/zEakjf08GYCyGu2zcnQhrqZG+IIlmfntYuiwfkab6EVMsatrj0dmI23TF1mN029EYVDUT4sbAv9WyukkMKuyx31K7MH882ifGLHG1JvIRnNrEh/qZ1ciyxSvfHUGeTplbZmkz3+7RAVCpgL//vqvZfQ0VGeCN8/NGYNvLtq+Gk3usj+VpTqDmvb6uV+XihPh54dHEmoNdAHj/MfNpWJVJ4GFa/FWclsuvOsZ1HGxQCKuhUK2YPyOuVgSA/m2rg6Rusiku0ynWlKEdEKH1Nvu81eXz1xTiQmseGWpRh+C+rquotN7mU/N1+TJhafNlSzlP4oIU0f2yulJeLZSvN21Ie3z0RC+zz5yoJXOqzNTny4yjiZ+zUzZM2TsbBlXNhPgH1TSBV06e5Go6pdNU/OzwR/EvJjlApieeO2QBVEuN8bjIt00BjLk830+9B2ffedApiwGanty8ZLkHponlchuOGEfpQvy9LI46bHt5oHR5xZ/utDjKKH9prXcLs9G0G1WrGc9fL4IgCNL+cXVeH4+aR0YKS4xfAH4rKMEpCxvqypOefb08cFZWgV+ssm36f+tsy7dNP5emIrTeZoHKt5PvMgusu9WSQ2VJH5NK5NOGtLfS0ron+1UvLpEHTqbCWmqkOm3iohqDQZCmPLu00mLHzEGYNqQDHuoRhb2vDMEnTydKj+8RHYA5I7s2+/IJ7kLcoube9o7LcawvhvXNRFyIHw5cLMCBCwVW28inTiKacGdwsSK3veqpyOs5RWg10h/aCffE47Md5xTbf4grkUxPruIqQWf+I/1w79ZSDs3xOcPx9f5LSP1vNi7eKIYgCBb7nhDuj33nb1pdVRUb4ocTbw6HpoXa6ntvJUsGHtHdfHm8OIrw8/F8pHx9CEG+nrhZXFFjQF8f7cKN38C9PdVYuOkUAGVRyz6xQVJ5gZhgX9zbPhSnq4pMiisUTQNuR67+s6S2gL5tmJ/ZatXE2CDFZ3zHzEGKqv628jJJEu7Squ65lZPvayftFGA6lWtKnpT84aZTeD/9V+l6pcFg9h6GdI7AZ+P6oG2YP+JdcDSGrBM/K7m11FV0RhypaiYuVa3mqmn/J3n+Qk07udvbtpcH4k93xeHHv9xrl+dLbFOdCyKf5hNPpNZGw1ZYWWburF4b0QWT7muH76feDQ+1SjF1NeFz89pEgLEEAwA81N36diPenh41BpPyj4alFUXyFaZrD/4m7QFXn5IKANBWdsL835S7Zbf7m72ePAh5rE8MPh3bB7tmDbY6bfWUSdX+mpb7O0JtOY2P9olRBIaWisc2ZB87eWmE0FqCIktign2x95X7cXzOcJsf88uvVxUBFVBdl01OrVbh/s4RDKhqIZYeOXy5wLEdqYOMM8bVrfaqD9iUGFQ1E+LWI7vOmC/FFp29Vv+tIhoiXOuNN37ftcGrDEVtLawABABt1dSOWOgRgGLl4aBO1avDLK0idDZBfl6Y9UAndI8OBKCc9tt8Il9abVdaoZcui98ASyvrv7xanphuuj8gAKkUhKkQG+pEWXJZtmqph6zAq6WtaOT79qnVKgztEiEtsy+TVesWS4Z4eqgdWkahNq0DfZA+fQC+nXwXnrk7Hn97ULkaMi7EF61kU/XiCOvK8cYvCFEB3g0abZ1wTzziQ/3QNz4YvWTHvi7CW3rXupoXqK68b2mVWs96vjZVF9K0pU6bsxC/HN0oqr2Cv7Ph9F8zUVOejUic+nDF2iBy1la5WLrddLXg4qd6Y/fZ62a1mVyB6SjF9aJy+Hh5oN/bm1BSocfZdx6Upnjva2AgER3kg8s3S/CGhST+uxNC0TbUT5HDBNQ/X87aPo7xFlYGzUi2XuriwvXq/kTKRoDee7QH/rLqYL3qsTWF9lV5VYmxQVVFQKsr5wf6eGFYl0h8tsO4vbi4BH1gx3Ccnzeiwa+tVquw5aWBDX4eW9SUUF3X0gpU7a52oThjYWNxZyaWY+lW9YXRlbj22ZNs9jvZdI9gZet3cXi4PkmtzubbyUkY1iUCu1Pvl26ztmJI7sFurTBn5B02BaHOxjRh+eqtMqSsPiTVNDonCyoaUioDANa/cC9WP9ffYiFHD7UKmy2ciE1zdGxlbdTQdIq6V5tAqRK/JXdY2YMyrKUGX07sj+Q6lhxwBNPE9ZbeLRS1mVx5Kqymau01Fa0l96OrWoRiqR6bs2NQ1UyI+54B1jeNFZem33bhLQJEibHB+GRsnxpPsu5GpVJhqyyY+a2gBFdkUymr9l6ULjckzwYw5h71axtSY+5dU+R8y4sD3t0utMa27vC5Ns0NC/LzwgDZCqn2jbSVUlOoKeiuT6I9uS4xTcHWCv7OhEFVMyGvli3mLpgSc05iQ9zzD5jpSq/XftfFSkvXJZ9CyThzHSfzqlfDfbr9nHS5vqNGdfHGQ9VTgx8+0auGljWraeri7T9Ub7Jd276WQ2T1jywViXVV8lV1rjxL1j3a+gh5jxruI/cjrtI13TjcFTCnqhnae+6Gxb2/xNyVO+NMq0y7B9NRq7sTXG9liS38vDxQVK7Hyl3nHdqPsUmxaBPiiyBfrwYlGoulMCzpHVv9vLUlQw/sGIZPx/ZRJLO7ou0zBuHfGefx7L3G/SdD/Lzg7alGaYUBf+hd8756zszTQ41XR3TGW+uPm93X0Olqci1eLdQorzQwqCLXMD/tBJ65J97s9ptVRRtb2rmekLNqH+5c1bPtJcDHU7H1jyn5iE1jUqlUNu23V5u/DGmPVgHeFvsdLqun9sAdNedEqVSqGjf3dhUxwb54RbZBeAsPNTJm3Y/SSj1aBbjedIncs/e2xYR74qFSqRA3a710+0A7fI6aO1s2FHcW4qSCK+7l6D5j4FQrcdsJK3nqUsFAd81fMN0up6H79zmrmvbxA+pfL8pRtN6eePbetlZXh51550HsfeV+9HPBmjb2EuTn5fIBlUgsAbH8T8bP8aOJ0TaVZCDLcqoKaJpuL+XM5LXnXA2Dqmbkj1WFDmury9MU+TaO0FymELpY2MdQTpw2chceapVixIrcw+BOETg/bwTefdR8/0mynbivqatsqFxYXJ3zK88FdhXuefYki8SCavkWEnov3aje9b6hK8NcQU1Jsa7OQ63CO3/ohrahfnhzZFese+Ee6b672oVIhWCJyP3VtMG6M6qQrbZwtr04bdE8kmcIgLGyNgCczTevnP7Ep7uly7Xt0eUODss223VHT/Zro9jM9tOxfRDi74Xesi18iIgay64z17Dj1DX8dVjHOqVa6KysTncVDh2pWrJkCbp37w6tVgutVoukpCRs2LBBun/gwIFQqVSKn0mTJime4+LFixgxYgR8fX0RHh6Ol19+GZWVyno0W7duRe/evaHRaJCQkICVK1ea9WXRokWIi4uDt7c3+vXrh7179yruLy0txZQpUxASEgJ/f3+MHj0aeXl59jsYTSCkqqDlLQv1euS1q0xzj9zRve1rrmnkboZ2iWBARURNwmAQ8OSne7B46xlM/LflfUitOX/dtaq/m3Lo2TM6Ohrz5s1DZmYm9u/fj8GDB2PkyJE4evSo1GbixInIycmRfubPny/dp9frMWLECJSXl2PXrl34/PPPsXLlSsyePVtqc+7cOYwYMQKDBg1CVlYWpk2bhmeffRZpaWlSm9WrVyMlJQWvv/46Dhw4gB49eiA5ORn5+dX7mk2fPh0//PAD1qxZg23btuHKlSt4+OGHG/kI2VeoDSsp/uKC27PUh7sm4xMROZq8Evq5a3ULknafvWHv7jQph07/PfTQQ4rrb7/9NpYsWYLdu3eja1dj4UBfX19ERlpeKv3TTz/h2LFj+PnnnxEREYGePXvizTffxMyZM/HGG2/Ay8sLS5cuRXx8PN577z0AQOfOnbFjxw4sWLAAycnJAID3338fEydOxPjx4wEAS5cuxfr167F8+XLMmjULhYWF+Oyzz/Dll19i8ODBAIAVK1agc+fO2L17N/r372+xf2VlZSgrq85f0ul0DThaDRcdWB1IlFXqoWlhnGsXN9sFgBHdWzV5v5rSe4/2QNalAkwf2sHRXSEickt7z1UHRhdl+bq2uHbbeM5s5aK7YTjNPI9er8eqVatQVFSEpKQk6fYvvvgCoaGhuOOOO5Camori4ur/oIyMDHTr1g0REdW1Z5KTk6HT6aTRroyMDAwZMkTxWsnJycjIyAAAlJeXIzMzU9FGrVZjyJAhUpvMzExUVFQo2nTq1Alt2rSR2lgyd+5cBAQESD8xMTH1OTR2ExFQnSslX7IqT1yvaf8tdzA6MRpvjroDwX6uV/+EiMgVyFNM9AYrNXysOJVnzPl1xS1qgHoGVZcuXcLly5el63v37sW0adPwySef1Pm5srOz4e/vD41Gg0mTJmHt2rXo0sVY2O7JJ5/Ef/7zH2zZsgWpqan4v//7P/zxj3+UHpubm6sIqABI13Nzc2tso9PpUFJSgmvXrkGv11tsI38OLy8vBAYGWm1jSWpqKgoLC6WfS5cu1eHI2F8LdfV/95d7qveByymszqdy19pNRETNWWFJ02xK/f2hK3jxq4OK2+JmrcfGI9bPlXLZvxkXEcWFuObm4PUKqp588kls2bIFgDHgGDp0KPbu3YtXXnkFc+bMqdNzdezYEVlZWdizZw8mT56McePG4dixYwCA5557DsnJyejWrRueeuop/Pvf/8batWtx5syZ+nS7yWk0GikJX/xxJHnA9I+NJ6TL568ZR/+8mkGCOhFRcyLu6XrsStOkn5gGVKJJ/8m06fHtq7aR6tzKsefL+qrXWfTIkSPo27cvAODrr7/GHXfcgV27duGLL76wuLKuJl5eXkhISEBiYiLmzp2LHj16YOHChRbb9uvXDwBw+vRpAEBkZKTZCjzxupiHZa2NVquFj48PQkND4eHhYbGN/DnKy8tRUFBgtY0ra+FhDLbK9a5bxZaIiMyJdarEOoWOZLBhKvBUVcmfTi5aT69eQVVFRQU0GmN+zs8//4zf//73AIx5Rjk5OQ3qkMFgUCR3y2VlZQEAWrUyJlMnJSUhOztbsUovPT0dWq1WmkJMSkrCpk2bFM+Tnp4u5W15eXkhMTFR0cZgMGDTpk1Sm8TERHh6eiranDx5EhcvXlTkf7kaMUH9SNVwa2Isl9wTEbkT043km9KT/dooVp2fq0O5BFfdmqheQVXXrl2xdOlSbN++Henp6Rg+fDgA4MqVKwgJsX3/rdTUVPzyyy84f/48srOzkZqaiq1bt+Kpp57CmTNn8OabbyIzMxPnz5/H999/j7Fjx2LAgAHo3r07AGDYsGHo0qULnn76aRw6dAhpaWl49dVXMWXKFCnomzRpEs6ePYsZM2bgxIkTWLx4Mb7++mtMnz5d6kdKSgo+/fRTfP755zh+/DgmT56MoqIiaTVgQEAAJkyYgJSUFGzZsgWZmZkYP348kpKSrK78c1ZbXxooXRb3hBJXZ4gbKhMRETXEivF34p0/dMP+V4dKt2WcuV7jY3Sl1YU/XTWnql4lFf7xj3/gD3/4A959912MGzcOPXoY92b6/vvvpWlBW+Tn52Ps2LHIyclBQEAAunfvjrS0NAwdOhSXLl3Czz//jA8++ABFRUWIiYnB6NGj8eqrr0qP9/DwwLp16zB58mQkJSXBz88P48aNU+R1xcfHY/369Zg+fToWLlyI6OhoLFu2TCqnAABjxozB1atXMXv2bOTm5qJnz57YuHGjInl9wYIFUKvVGD16NMrKypCcnIzFixfX5/A5lHxT2os3itEqwBtpR41Tn2evunbRNSIicpyr8pXkoeYryU/m3qrx8Xp99fSg1gW3qAHqGVQNHDgQ165dg06nQ1BQ9ZTRc889Bz8/26PLzz77zOp9MTEx2LZtW63PERsbix9//LHW/h48aDl5TjR16lRMnTrV6v3e3t5YtGgRFi1aVGufXEVRWaWikvqwLhE1tCYiIrJOHjS1CanO4Wob5oezV4sQVEspmxtuMFtSr+m/wYMH49atW4qACgCCg4MxZswYu3SMGk/HCGMC4K95t6SkQAB47zHuBk9ERPVjbYuZexKM24IdvlxQ4+NLyqsrsbtqeZ96BVVbt25Febl5RFlaWort27c3uFPUuMR56/PXixV5VC29XXO4lYiIamYQ6laEsz5aWAmExBmRzAs3bXoeV62mDtRx+u/w4cPS5WPHjikKX+r1emzcuBGtW7e2X++oUQT6eiGnsBRHr+gQqTV+eDtEuHcldSKi5kiMpfadty2gaYhDVSNRgzuFK25vG+qHzQBulVaaP0imwg3K+tQpqOrZsydUKhVUKpW0B56cj48PPvroI7t1jhpH1ygtjufocDxHhz/0igIA+Hi65vJVIiKyTtPCOCHVFNu+hPobV91fKShR3H5P+1As23EOAFBYUoEAK0no4kr060Wum1tVp6Dq3LlzEAQBbdu2xd69exEWFibd5+XlhfDwcHh48OTs7O6I0uKbquK24nY17r7nHxFRcxSubbqptIMXCwAA/dsqSysF+VYnqPf4+084P2+ExceLAaD4ryuqU1AVGxsLwFgck1zX433b4I0fjFsBnb9u/GZw6WbddhInIiKSEzeqv2xyPmkbpqwKcP12GUKqRrXkKqpKKrhqNXWgniUVAODUqVPYsmUL8vPzzYKs2bNnN7hj1Hi8LUz1JbW1vWgrERGRqeM5xv0F+8YHK243XQR17lqRxaBKLMlQUqE3u89V1Cuo+vTTTzF58mSEhoYiMjISKlV1xr9KpWJQ5QKeH9AW//rlrHRd3wQrQ4iIyH15VU3b6UrME9J3zRqMu+ZtBgCcvVaEPnHBZm3EbWyKy5pZUPXWW2/h7bffxsyZM+3dH2oiwSZF2MoqOKVLRET15+dlDCnaW1hNHhVYnSh/IsdyZXVxw+VAX9ct71OvbLCbN2/i0UcftXdfqAmN7KksfREb4vgdzImIyHWdyjcGS14elkOLNsHG84y4ys+UuGqwe3Sg/TvXROoVVD366KP46aef7N0XakIRWuV89iOJMQ7qCRERuYPbZcZpv3Ir9aZigo2jVdZW94kLp0Jq2c7GmdVr+i8hIQGvvfYadu/ejW7dusHTUzlU9+KLL9qlc9R4VCoVWmpa4FbVL4GPF0thEBFR/YX5a3ClsBTRQZZnPhJjg7Hz9HWsz86BuItuSbkemRduom98MEL8vVBYUuHS56N6BVWffPIJ/P39sW3bNrNNj1UqFYMqF7HuxXvw6fazmHRfO0d3hYiIGpG4VUxjqkqJsjr951G1qE2lAh5bmgG9IEhb10wb0h5nrxoT1dtHNLOSCufOnbN3P8gBYkP88Naobo7uBhERNRK/Jhr1EQQBubpS42VYXk1+V0IIFvxs3Dpn7/kbivs++PmUdNnTwzU3UwbqmVNFREREzi+oifKT5FV5TFeXiwKtbE9jKj7Ur/ZGTqpeI1XPPPNMjfcvX768Xp0hIiIi1yMfm/L1shxaWNvzz1SIn3lhUFdRr6Dq5k3lbtcVFRU4cuQICgoKLG60TERERO5Lvomyn8bylGNYS/NgacX4OzF+xT7FbV7NZe8/0dq1a81uMxgMmDx5Mtq1Y9IzERFRc2KQzf9pWlgOquS7rwDAquf6o2dMYGN2q8nZLRxUq9VISUnBggUL7PWURERE5ALqs7owMTbIbC/aUH/XrVEF2DlR/cyZM6isNN/zh4iIiNyX3lD3/WM9LZRe8NfUawLNadSr9ykpKYrrgiAgJycH69evx7hx4+zSMSIiInINJ3ONW9S0lu3xZ0nK0A54P/1XRdV0b081Sqv2n/X3boZB1cGDBxXX1Wo1wsLC8N5779W6MpCIiIiant4gwEPdODWgsi4VAAB0JRU1tnthcAL6xAYpCnx2itRKj48OdO19aOsVVG3ZssXe/SAiIiI7k5cxyCkssbqFTEMVVgVTkQHeNbZTqVS4KyFUcVvrQB8pqDp46aaFR7mOBuVUXb16FTt27MCOHTtw9epVe/WJiIiI7MBS3lJjOHqlEABwZ3xwnR+bMqyDdPmbSXfZrU+OUK+RqqKiIrzwwgv497//DYPBOA/q4eGBsWPH4qOPPoKvr2sP3xEREbkLTQs1yioNjfoaN4uNI1WttDWPVFnSLswf5+eNsHeXHKJeIWxKSgq2bduGH374AQUFBSgoKMD//vc/bNu2DX/961/t3UciIiJyYhFaY2HP6OCaE9XdXb1Gqr799lt88803GDhwoHTbgw8+CB8fHzz22GNYsmSJvfpHRERETi5PVwYAiAtx3X377KFeI1XFxcWIiIgwuz08PBzFxcUN7hQRERG5HtOq6c1NvYKqpKQkvP766ygtLZVuKykpwd///nckJSXZrXNERERkH0Ld63PWWVRg3XOq3Em9pv8++OADDB8+HNHR0ejRowcA4NChQ9BoNPjpp5/s2kEiIiKqPzFJ/beCEsQE228hmcEg4J8/nVRsNaMCR6rqrFu3bjh16hTmzp2Lnj17omfPnpg3bx5Onz6Nrl272vw8S5YsQffu3aHVaqHVapGUlIQNGzaYtRMEAQ888ABUKhW+++47xX0qlcrsZ9WqVYo2W7duRe/evaHRaJCQkICVK1eavcaiRYsQFxcHb29v9OvXD3v37lXcX1paiilTpiAkJAT+/v4YPXo08vLybH6vREREjmTvcOfAxZtYvPUM3k//VbqtpYtXRG+oer37uXPnIiIiAhMnTlTcvnz5cly9ehUzZ8606Xmio6Mxb948tG/fHoIg4PPPP8fIkSNx8OBBRXD2wQcf1DhPu2LFCgwfPly6HhgYKF0+d+4cRowYgUmTJuGLL77Apk2b8Oyzz6JVq1ZITk4GAKxevRopKSlYunQp+vXrhw8++ADJyck4efIkwsPDAQDTp0/H+vXrsWbNGgQEBGDq1Kl4+OGHsXPnTpveKxERkSO0C/PDmatFdn9eS/v9mW6Q3OwI9RAbGyvs3LnT7Pbdu3cLcXFx9XlKSVBQkLBs2TLp+sGDB4XWrVsLOTk5AgBh7dq1ivaWbpObMWOG0LVrV8VtY8aMEZKTk6Xrffv2FaZMmSJd1+v1QlRUlDB37lxBEAShoKBA8PT0FNasWSO1OX78uABAyMjIsPm9FRYWCgCEwsJCmx9DRETUEIP/uUWInblO2H3mml2fd/W+i0LszHXST4dXfrTr8zsTW8/f9Zr+y83NRatWrcxuDwsLQ05OTr2CO71ej1WrVqGoqEhKdi8uLsaTTz6JRYsWITIy0upjp0yZgtDQUPTt2xfLly+HIMvGy8jIwJAhQxTtk5OTkZGRAQAoLy9HZmamoo1arcaQIUOkNpmZmaioqFC06dSpE9q0aSO1saSsrAw6nU7xQ0RE5A5KK/SK64amyIR3cvWa/ouJicHOnTsRHx+vuH3nzp2Iioqq03NlZ2cjKSkJpaWl8Pf3x9q1a9GlSxcAxim3u+66CyNHjrT6+Dlz5mDw4MHw9fXFTz/9hD//+c+4ffs2XnzxRQDGANC0/ENERAR0Oh1KSkpw8+ZN6PV6i21OnDghPYeXl5diWlFsk5uba7Vvc+fOxd///nebjwUREZGryC0sVVyv0DOoqldQNXHiREybNg0VFRUYPHgwAGDTpk2YMWNGnSuqd+zYEVlZWSgsLMQ333yDcePGYdu2bTh9+jQ2b96MgwcP1vj41157Tbrcq1cvFBUV4d1335WCKkdKTU1FSkqKdF2n0yEmJsaBPSIiouYqV1dae6M68DdJSm8X1rwLfwL1DKpefvllXL9+HX/+859RXl4OAPD29sbMmTORmppap+fy8vJCQkICACAxMRH79u3DwoUL4ePjgzNnzpiNDo0ePRr33nsvtm7davH5+vXrhzfffBNlZWXQaDSIjIw0W6WXl5cHrVYLHx8feHh4wMPDw2IbccoxMjIS5eXlKCgoUPRH3sYSjUYDjUZj45EgIiKyPzFJvdzO+/9tOp5v8XWas3rlVKlUKvzjH//A1atXsXv3bhw6dAg3btzA7NmzG9whg8GAsrIyzJo1C4cPH0ZWVpb0AwALFizAihUrrD4+KysLQUFBUjCTlJSETZs2Kdqkp6dLeVteXl5ITExUtDEYDNi0aZPUJjExEZ6enoo2J0+exMWLF1nslIiInNq97UMBAGo7VzvPvHBTcf2ehFC7Pr8ralBBCX9/f9x55531fnxqaioeeOABtGnTBrdu3cKXX36JrVu3Ii0tDZGRkRZHgdq0aSPlcv3www/Iy8tD//794e3tjfT0dLzzzjt46aWXpPaTJk3Cxx9/jBkzZuCZZ57B5s2b8fXXX2P9+vVSm5SUFIwbNw59+vRB37598cEHH6CoqAjjx48HAAQEBGDChAlISUlBcHAwtFotXnjhBSQlJaF///71fv9ERESNzd7BFADk3zKfSjx6pdDur+NqHFqlKz8/H2PHjkVOTg4CAgLQvXt3pKWlYejQoTY93tPTE4sWLcL06dMhCAISEhLw/vvvK+pnxcfHY/369Zg+fToWLlyI6OhoLFu2TKpRBQBjxozB1atXMXv2bOTm5qJnz57YuHGjInl9wYIFUKvVGD16NMrKypCcnIzFixfb72AQERG5iHHL95nd9mgf5gyrBIFrIJuKTqdDQEAACgsLodVqHd0dIiJqBsYt34ttv17Fe4/2wOjEaLs859QvD2DdYWMJpX2vDEFhSQUSwv3t8tzOyNbzd71yqoiIiMg1iPWjsi4V2O05N58wJqk/1a8Nwlpq3DqgqgsGVURERG7sSkEJACDUv/6r0csq9Tidf0u6Ls5xVbI2lQKDKiIiIjfWr21Ig59j2qosDHn/F6w7fAUAUFJVTX1Qp7AGP7c7YVBFRERENdpwxLh7yNQvlQW5g/1Yi1GOQRURERHZTL7nX0ywjwN74nwYVBERETUDAuyT/3SzuFy6HNHS2y7P6S4YVBEREbkxg8EYTGWcuW6X50s/Vr2tWyPUFXVpDKqIiIjcWGFJBQAgNsS3Xo8vq9Qrrm8/dU26rGJUpcCgioiIyI11iw5o0ONNS4SfyNU16PncGYMqIiIistmlG8a6Vz1jAh3bESfEoIqIiIisqjRYTnDP15lvqtzcMagiIiIiq85dLbJ4e0tvzybuifNjUEVERNQM1Hfvv2u3yyze3icuqAG9cU8MqoiIiNzYjdvGulL1HVmyFIw93Ls13hp1R0O65ZYYVBEREbmxrq21AABNi/qd8r89cBkA0FLTQrrt77/vynIKFrSovQkRERG5Kg91w8ZPLt80rvYL9PPEf57tBw+1ivlUVjCoIiIiolr9dWhH9GAZhRpx+o+IiIhq1T7C39FdcHoMqoiIiJoB08rotqjUG6TL3p4eduyNe2JQRURE5MaEqmgq42zdN1SWx2Ghfho79ch9MagiIiJyYy2qEtXbBNd9Q+XbpZX27o5bY1BFRETkxiIDvOv92KXbzkiXtT5c21YbBlVERERk0e2y6pEq1qWqHYMqIiIisuiLPRcd3QWXwqCKiIiIatQ21M/RXXAJDKqIiIioRs/f19bRXXAJDKqIiIjIjMFQXVAh1J/lFGzBoIqIiKgZuHijuN6P7d0myI49cV8MqoiIiNxYQFUpBHUdF++Vy6qpk20YVBEREbkxrbcngLqXRMjTlUqXA3w87dond8WgioiIiMxcuF49Xaiu6zBXM+XQoGrJkiXo3r07tFottFotkpKSsGHDBrN2giDggQcegEqlwnfffae47+LFixgxYgR8fX0RHh6Ol19+GZWVyrL6W7duRe/evaHRaJCQkICVK1eavcaiRYsQFxcHb29v9OvXD3v37lXcX1paiilTpiAkJAT+/v4YPXo08vLyGnwMiIiInFHmhZuO7oLLcWhQFR0djXnz5iEzMxP79+/H4MGDMXLkSBw9elTR7oMPPrA4bKnX6zFixAiUl5dj165d+Pzzz7Fy5UrMnj1banPu3DmMGDECgwYNQlZWFqZNm4Znn30WaWlpUpvVq1cjJSUFr7/+Og4cOIAePXogOTkZ+fn5Upvp06fjhx9+wJo1a7Bt2zZcuXIFDz/8cCMcFSIiIsfTlVY4uguuR3AyQUFBwrJly6TrBw8eFFq3bi3k5OQIAIS1a9dK9/3444+CWq0WcnNzpduWLFkiaLVaoaysTBAEQZgxY4bQtWtXxWuMGTNGSE5Olq737dtXmDJlinRdr9cLUVFRwty5cwVBEISCggLB09NTWLNmjdTm+PHjAgAhIyPD5vdWWFgoABAKCwttfgwREVFD5BWWCLEz1wltU9fX6XHjV+wVYmeuE8Yt39NIPXMdtp6/nSanSq/XY9WqVSgqKkJSUhIAoLi4GE8++SQWLVqEyMhIs8dkZGSgW7duiIiIkG5LTk6GTqeTRrsyMjIwZMgQxeOSk5ORkZEBACgvL0dmZqaijVqtxpAhQ6Q2mZmZqKioULTp1KkT2rRpI7WxpKysDDqdTvFDRETkCHpZ3SlbtPQ2rhpkNpXtHB5UZWdnw9/fHxqNBpMmTcLatWvRpUsXAMYpt7vuugsjR460+Njc3FxFQAVAup6bm1tjG51Oh5KSEly7dg16vd5iG/lzeHl5ITAw0GobS+bOnYuAgADpJyYmppajQURE1HjKK20vkyDGYHfGBzdSb9xPC0d3oGPHjsjKykJhYSG++eYbjBs3Dtu2bcPp06exefNmHDx40NFdrLfU1FSkpKRI13U6HQMrIiJqUvJq6NeLytAqwMemx+06fQ0AINRtgKtZc3hQ5eXlhYSEBABAYmIi9u3bh4ULF8LHxwdnzpwxGx0aPXo07r33XmzduhWRkZFmq/TEFXnidGFkZKTZKr28vDxotVr4+PjAw8MDHh4eFtvIn6O8vBwFBQWK/sjbWKLRaKDRsLQ/ERE5jlqtgpeHus7FPK8XlQMA6ljeqllz+PSfKYPBgLKyMsyaNQuHDx9GVlaW9AMACxYswIoVKwAASUlJyM7OVqzSS09Ph1arlaYQk5KSsGnTJsVrpKenS3lbXl5eSExMVLQxGAzYtGmT1CYxMRGenp6KNidPnsTFixelNkRERO6kdaBxRCs22M/BPXEdDh2pSk1NxQMPPIA2bdrg1q1b+PLLL7F161akpaUhMjLS4ihQmzZtEB8fDwAYNmwYunTpgqeffhrz589Hbm4uXn31VUyZMkUaIZo0aRI+/vhjzJgxA8888ww2b96Mr7/+GuvXr5eeMyUlBePGjUOfPn3Qt29ffPDBBygqKsL48eMBAAEBAZgwYQJSUlIQHBwMrVaLF154AUlJSejfv38THCkiIqKmJY5QtQ6ybbqQHBxU5efnY+zYscjJyUFAQAC6d++OtLQ0DB061KbHe3h4YN26dZg8eTKSkpLg5+eHcePGYc6cOVKb+Ph4rF+/HtOnT8fChQsRHR2NZcuWITk5WWozZswYXL16FbNnz0Zubi569uyJjRs3KpLXFyxYALVajdGjR6OsrAzJyclYvHix/Q4GERERuTSVIDAFranodDoEBASgsLAQWq3W0d0hIqJmosMrG1CuNyAjdbDNiepxs4wzOt9NuRs9YwIbsXfOz9bzt9PlVBEREVHjqKjkOEpjYlBFRETk5sSVf1cKS+r8WDFhnWrHoIqIiMjNaauqo3uobauPIM8MYkkF2zGoIiIicnMh/nWrmZirK5Uu+2scXtLSZTCoIiIiIgX5PoHenh4O7IlrYVBFREREFmlaMEyoCx4tIiIiUrhVWunoLrgkBlVERETNxJUC21b/FZZUAADKKuu2X2Bzx6CKiIjIzZ27VgQAMNhY71tc8Nc2lPv+1QWDKiIiIjd3T0IoAEAF2+ojXLtdbmzPcgp1wqCKiIjIzdU1OCoqM+ZUnbla1Ai9cV8MqoiIiEhBLBLaNz7YwT1xLQyqiIiIyCIf1qiqEwZVREREzcSZq7cd3QW3xqCKiIjIzZ3MvQUAaKG27bQvllSgumFQRURE5Obu6xAGAPCw8az/W1U9q5vF5Y3VJbfEoIqIiMjNiYnnttJ6eyr+JdswqCIiIiKFq7dLAQBtQnwd3BPXwqCKiIiIFHILjUEVc6vqhkEVERERKQT7eQEAQqr+JdswqCIiImomSir0NrUTtwiMDPBuxN64HwZVREREbq6s0gAAyP5NZ1P7LSfzAVQHV2QbBlVERERuTm8wRketA20beeoQ0RIAUGbjyBYZMagiIiJycx0i/Ov1uISq4Ipsw6CKiIiIFE7lczub+mBQRURERArFZZUAgFJO/9UJgyoiIiJSqNAbc7Biglj8sy4YVBEREbk5cRXf1pNXbWpfrjeuFvT0qNv2Ns0dgyoiIiI3d7vcOJ2XEG5bwrq4V2DrIJ9G65M7YlBFRETk5jpHauvUXizB4KHiSFVdMKgiIiIiiSI5nTFVnTg0qFqyZAm6d+8OrVYLrVaLpKQkbNiwQbr/+eefR7t27eDj44OwsDCMHDkSJ06cUDyHSqUy+1m1apWizdatW9G7d29oNBokJCRg5cqVZn1ZtGgR4uLi4O3tjX79+mHv3r2K+0tLSzFlyhSEhITA398fo0ePRl5env0OBhERkRO4WVwuXQ7z1ziwJ67HoUFVdHQ05s2bh8zMTOzfvx+DBw/GyJEjcfToUQBAYmIiVqxYgePHjyMtLQ2CIGDYsGHQ65VLPFesWIGcnBzpZ9SoUdJ9586dw4gRIzBo0CBkZWVh2rRpePbZZ5GWlia1Wb16NVJSUvD666/jwIED6NGjB5KTk5Gfny+1mT59On744QesWbMG27Ztw5UrV/Dwww837gEiIiJyEE8P40AF1YHgZIKCgoRly5ZZvO/QoUMCAOH06dPSbQCEtWvXWn2+GTNmCF27dlXcNmbMGCE5OVm63rdvX2HKlCnSdb1eL0RFRQlz584VBEEQCgoKBE9PT2HNmjVSm+PHjwsAhIyMDKuvXVpaKhQWFko/ly5dEgAIhYWFVh9DRERkb2sPXBZiZ64TRi/eWWvbKwXFQuzMdULC39Y3Qc9cQ2FhoU3nb6fJqdLr9Vi1ahWKioqQlJRkdn9RURFWrFiB+Ph4xMTEKO6bMmUKQkND0bdvXyxfvhyCbAfIjIwMDBkyRNE+OTkZGRkZAIDy8nJkZmYq2qjVagwZMkRqk5mZiYqKCkWbTp06oU2bNlIbS+bOnYuAgADpx7TfRERETaGyKvF8/4Wbtba9frtc8RiyncODquzsbPj7+0Oj0WDSpElYu3YtunTpIt2/ePFi+Pv7w9/fHxs2bEB6ejq8vLyk++fMmYOvv/4a6enpGD16NP785z/jo48+ku7Pzc1FRESE4jUjIiKg0+lQUlKCa9euQa/XW2yTm5srPYeXlxcCAwOttrEkNTUVhYWF0s+lS5fqfHyIiIgayl/jAQCID/Wrte3tqmrqAmOqOmvh6A507NgRWVlZKCwsxDfffINx48Zh27ZtUmD11FNPYejQocjJycE///lPPPbYY9i5cye8vY07bb/22mvSc/Xq1QtFRUV499138eKLLzrk/chpNBpoNEzyIyIixwqpQ8L5uWtFAICwljx/1ZXDR6q8vLyQkJCAxMREzJ07Fz169MDChQul+wMCAtC+fXsMGDAA33zzDU6cOIG1a9dafb5+/frh8uXLKCsrAwBERkaardLLy8uDVquFj48PQkND4eHhYbFNZGSk9Bzl5eUoKCiw2oaIiMgd+HgaR7VKy7nvX105PKgyZTAYpIDIlCAIEATB6v0AkJWVhaCgIGmEKCkpCZs2bVK0SU9Pl/K2vLy8kJiYqGhjMBiwadMmqU1iYiI8PT0VbU6ePImLFy9azP8iIiJyVYcuFwAAuscEOLYjLsih03+pqal44IEH0KZNG9y6dQtffvkltm7dirS0NJw9exarV6/GsGHDEBYWhsuXL2PevHnw8fHBgw8+CAD44YcfkJeXh/79+8Pb2xvp6el455138NJLL0mvMWnSJHz88ceYMWMGnnnmGWzevBlff/011q9fL7VJSUnBuHHj0KdPH/Tt2xcffPABioqKMH78eADG0bIJEyYgJSUFwcHB0Gq1eOGFF5CUlIT+/fs37UEjIiKqJ3Fqz5L/yziP1/53FC01xtDg8s2SpuqW23BoUJWfn4+xY8ciJycHAQEB6N69O9LS0jB06FBcuXIF27dvxwcffICbN28iIiICAwYMwK5duxAeHg4A8PT0xKJFizB9+nQIgoCEhAS8//77mDhxovQa8fHxWL9+PaZPn46FCxciOjoay5YtQ3JystRmzJgxuHr1KmbPno3c3Fz07NkTGzduVCSvL1iwAGq1GqNHj0ZZWRmSk5OxePHipjtYRERE9WSwYSXfa/8z1oi8VZWo3jGiZaP2yR2pBIH5/U1Fp9MhICAAhYWF0Grrtg8TERFRff1WUIK7522Gl4cav779gMU2cbPWK66PTYrFnJF3NEX3nJ6t52+ny6kiIiIi+1LXUhjd0viKmtXU64xBFRERUTOXf8t8AVgPJqrXGYMqIiKiZqLSYLB4+/7z5pXWOVJVdwyqiIiImgmDAFTozQOr0/m3zW67My64KbrkVhhUERERubnwlt7S5ZtF5Wb3+1VtYwMAnh4qDOkcgahAnybpmztx+DY1RERE1Lg81Cp4qFXQ11Ja4Q+9WmPBmJ5N0yk3xJEqIiKiZq6wpMLRXXALDKqIiIiaAbFsgqUA6mxVpfVbpQyuGoJBFRERUTMgzvzpSivN7gvzN+6X6+3pYXYf2Y5BFRERUTPQJtjX6n0nc28BAOJC/JqqO26JQRUREREBAK5bWBlItmNQRURE1Mz5aYzFABLC/R3cE9fGoIqIiKgZuV1mnlMl8vNiTlVDMKgiIiJqBi7eKAYAXLWwz9/Px/OaujtuiUEVERFRMyBO7Xm1MD/1d4psCQCoqKU4KNWMQRUREVEzEN5SU2ubeK7+axAGVURERM1ISbl5TtWJqpIK1DAMqoiIiJoBMUH97NUiq20qDYam6o5bYlBFRETUDKhUKgBAgK+n4vZi2chVx6rcKqofBlVERETNQMcIyzWoKvTVyemh/rXnXZF1DKqIiIiakUq9coVfsYUcK6ofBlVERETNgDgitf/CTcXtl2+WSJdbqFVN2id3w6CKiIioGRBHpKKDfKTbzl8rwk9Hc6XrYt4V1U8LR3eAiIiIGl+XVgFIO1pdOb1Sb8DAf26Vrof4eTmgV+6FI1VERETN0LEcneJ6EIOqBmNQRURE1IwculQAANh3XplbxXSqhmNQRURE1Axcu23cSFnrbaxTVaFXFvq8q11ok/fJ3TCoIiIiaga6tQ4AUL2h8q95yq1pBnUKb/I+uRsmqhMRETUDHibzewE+xhGrIF9PzH6oCwa050hVQzGoIiIiaob2nb8BABjZszX+0Cvawb1xDw6d/luyZAm6d+8OrVYLrVaLpKQkbNiwQbr/+eefR7t27eDj44OwsDCMHDkSJ06cUDzHxYsXMWLECPj6+iI8PBwvv/wyKiuV1WG3bt2K3r17Q6PRICEhAStXrjTry6JFixAXFwdvb2/069cPe/fuVdxfWlqKKVOmICQkBP7+/hg9ejTy8vLMnoeIiMiZiZsmt1CrFdep4RwaVEVHR2PevHnIzMzE/v37MXjwYIwcORJHjx4FACQmJmLFihU4fvw40tLSIAgChg0bBr1eDwDQ6/UYMWIEysvLsWvXLnz++edYuXIlZs+eLb3GuXPnMGLECAwaNAhZWVmYNm0ann32WaSlpUltVq9ejZSUFLz++us4cOAAevTogeTkZOTn50ttpk+fjh9++AFr1qzBtm3bcOXKFTz88MNNdKSIiIgaxiAYK6rvPH0dAJCnKwUAtA70dVif3I7gZIKCgoRly5ZZvO/QoUMCAOH06dOCIAjCjz/+KKjVaiE3N1dqs2TJEkGr1QplZWWCIAjCjBkzhK5duyqeZ8yYMUJycrJ0vW/fvsKUKVOk63q9XoiKihLmzp0rCIIgFBQUCJ6ensKaNWukNsePHxcACBkZGTa/t8LCQgGAUFhYaPNjiIiI7OG7g5eF2JnrhMH/3CIIgiDEzlwnxM5cJ3yz/5JjO+YCbD1/O83qP71ej1WrVqGoqAhJSUlm9xcVFWHFihWIj49HTEwMACAjIwPdunVDRESE1C45ORk6nU4a7crIyMCQIUMUz5WcnIyMjAwAQHl5OTIzMxVt1Go1hgwZIrXJzMxERUWFok2nTp3Qpk0bqY0lZWVl0Ol0ih8iIiJHiNR6K64HVxX7DPLzdER33JLDg6rs7Gz4+/tDo9Fg0qRJWLt2Lbp06SLdv3jxYvj7+8Pf3x8bNmxAeno6vLyMH4Tc3FxFQAVAup6bm1tjG51Oh5KSEly7dg16vd5iG/lzeHl5ITAw0GobS+bOnYuAgADpRwwGiYiIHE1dtc9fVKBPLS3JVg4Pqjp27IisrCzs2bMHkydPxrhx43Ds2DHp/qeeegoHDx7Etm3b0KFDBzz22GMoLS11YI9tl5qaisLCQunn0qVLju4SERE1c7pS42IusRgo2Y/DSyp4eXkhISEBgDExfd++fVi4cCH+9a9/AYA0ytO+fXv0798fQUFBWLt2LZ544glERkaardITV+RFRkZK/5qu0svLy4NWq4WPjw88PDzg4eFhsY38OcrLy1FQUKAYrZK3sUSj0UCj0dTjqBAREdmXvipR/eqtMpRW6KXbxXpV1HAOH6kyZTAYUFZmOXoWBAGCIEj3JyUlITs7W7FKLz09HVqtVppCTEpKwqZNmxTPk56eLuVteXl5ITExUdHGYDBg06ZNUpvExER4enoq2pw8eRIXL160mP9FRETkbCKqcqq8PdUoLq8Oqlp6M6iyF4eOVKWmpuKBBx5AmzZtcOvWLXz55ZfYunUr0tLScPbsWaxevRrDhg1DWFgYLl++jHnz5sHHxwcPPvggAGDYsGHo0qULnn76acyfPx+5ubl49dVXMWXKFGmEaNKkSfj4448xY8YMPPPMM9i8eTO+/vprrF+/XupHSkoKxo0bhz59+qBv37744IMPUFRUhPHjxwMwjpZNmDABKSkpCA4OhlarxQsvvICkpCT079+/6Q8cERFRHWlaVI+jXCkokS77axw+aeU2HHok8/PzMXbsWOTk5CAgIADdu3dHWloahg4diitXrmD79u344IMPcPPmTURERGDAgAHYtWsXwsON+xN5eHhg3bp1mDx5MpKSkuDn54dx48Zhzpw50mvEx8dj/fr1mD59OhYuXIjo6GgsW7YMycnJUpsxY8bg6tWrmD17NnJzc9GzZ09s3LhRkby+YMECqNVqjB49GmVlZUhOTsbixYub7mARERHZgSAAv/toh6O74ZZUglA1yUqNTqfTISAgAIWFhdBqtY7uDhERNSOXbxbjnn9sMbv9/LwRDuiNa7H1/O10OVVERERkfyqVqvZG1CAMqoiIiJqBqADv2htRgzCoIiIiagYsjVQt/WNvB/TEfTGoIiIiaqaG39HK0V1wKwyqiIiImiEvD4YA9sYjSkRE1AwN6BDm6C64HQZVREREzURr2ebJZ67edmBP3BODKiIiomYi/1apdLlvXLADe+KeGFQRERE1E8F+XtLlovJKB/bEPTGoIiIiaibiQvyky8ypsj8GVURERM2E3lC9M11CuL8De+KeGFQRERE1E7m66pyqFmpuW2NvDKqIiIiaifCWGulylGwlINkHgyoiIqJmoo9sxZ8HN1i2OwZVREREzcRvBSXS5QAfTwf2xD21cHQHiIiIqGlMH9IBhy8XYEjnCKiZU2V3DKqIiIiaiYRwf2yfMdjR3XBbnP4jIiIisgMGVURERER2wKCKiIiIyA4YVBERERHZAYMqIiIiIjtgUEVERERkBwyqiIiIiOyAQRURERGRHTCoIiIiIrIDBlVEREREdsCgioiIiMgOGFQRERER2QGDKiIiIiI7YFBFREREZActHN2B5kQQBACATqdzcE+IiIjIVuJ5WzyPW8OgqgndunULABATE+PgnhAREVFd3bp1CwEBAVbvVwm1hV1kNwaDAVeuXEHLli2hUqkc3R2no9PpEBMTg0uXLkGr1Tq6O26Jx7jx8Rg3Ph7jxsdjrCQIAm7duoWoqCio1dYzpzhS1YTUajWio6Md3Q2np9Vq+UvcyHiMGx+PcePjMW58PMbVahqhEjFRnYiIiMgOGFQRERER2QGDKnIaGo0Gr7/+OjQajaO74rZ4jBsfj3Hj4zFufDzG9cNEdSIiIiI74EgVERERkR0wqCIiIiKyAwZVRERERHbAoIqIiIjIDhhUkd3MnTsXd955J1q2bInw8HCMGjUKJ0+eVLR5/vnn0a5dO/j4+CAsLAwjR47EiRMnFG1efPFFJCYmQqPRoGfPnhZf6/Dhw7j33nvh7e2NmJgYzJ8/v7HellNpqmN8/vx5qFQqs5/du3c35ttzCvY4xocOHcITTzyBmJgY+Pj4oHPnzli4cKHZa23duhW9e/eGRqNBQkICVq5c2dhvzyk01THeunWrxc9xbm5uk7xPR7LHMb5+/TqGDx+OqKgoaDQaxMTEYOrUqWb71zbXz7ElDKrIbrZt24YpU6Zg9+7dSE9PR0VFBYYNG4aioiKpTWJiIlasWIHjx48jLS0NgiBg2LBh0Ov1iud65plnMGbMGIuvo9PpMGzYMMTGxiIzMxPvvvsu3njjDXzyySeN+v6cQVMdY9HPP/+MnJwc6ScxMbFR3pczsccxzszMRHh4OP7zn//g6NGjeOWVV5CamoqPP/5Yeo5z585hxIgRGDRoELKysjBt2jQ8++yzSEtLa/L33NSa6hiLTp48qfgch4eHN9l7dRR7HGO1Wo2RI0fi+++/x6+//oqVK1fi559/xqRJk6TnaM6fY4sEokaSn58vABC2bdtmtc2hQ4cEAMLp06fN7nv99deFHj16mN2+ePFiISgoSCgrK5NumzlzptCxY0e79NuVNNYxPnfunABAOHjwoB1765oaeoxFf/7zn4VBgwZJ12fMmCF07dpV0WbMmDFCcnJywzvtYhrrGG/ZskUAINy8edOe3XVJ9jrGCxcuFKKjo6Xr/BwrcaSKGk1hYSEAIDg42OL9RUVFWLFiBeLj4xETE2Pz82ZkZGDAgAHw8vKSbktOTsbJkydx8+bNhnXaxTTWMRb9/ve/R3h4OO655x58//33Deqrq7LXMS4sLFQ8R0ZGBoYMGaJok5ycjIyMDDv02rU01jEW9ezZE61atcLQoUOxc+dO+3TaxdjjGF+5cgX//e9/cd9990m38XOsxKCKGoXBYMC0adNw991344477lDct3jxYvj7+8Pf3x8bNmxAenq6IkCqTW5uLiIiIhS3idebQ66EqDGPsb+/P9577z2sWbMG69evxz333INRo0Y1u8DKXsd4165dWL16NZ577jnpNmufY51Oh5KSEvu/GSfVmMe4VatWWLp0Kb799lt8++23iImJwcCBA3HgwIFGfU/OpqHH+IknnoCvry9at24NrVaLZcuWSffxc2zC0UNl5J4mTZokxMbGCpcuXTK7r6CgQPj111+Fbdu2CQ899JDQu3dvoaSkxKydtampoUOHCs8995zitqNHjwoAhGPHjtntPTi7xjzGljz99NPCPffc09BuuxR7HOPs7GwhNDRUePPNNxW3t2/fXnjnnXcUt61fv14AIBQXF9v3jTixxjzGlgwYMED44x//aJe+u4qGHuOcnBzh+PHjwv/+9z+hS5cuwuTJk6X7+DlWYlBFdjdlyhQhOjpaOHv2bK1ty8rKBF9fX+HLL780u8/aCf/pp58WRo4cqbht8+bNAgDhxo0b9e22S2nsY2zJxx9/LERGRta1qy7LHsf46NGjQnh4uPC3v/3N7DH33nuv8Je//EVx2/LlywWtVtugfruSxj7Glrz00ktC//7969VfV2SvvxWi7du3CwCEK1euCILAz7EpTv+R3QiCgKlTp2Lt2rXYvHkz4uPjbXqMIAgoKyuz+XWSkpLwyy+/oKKiQrotPT0dHTt2RFBQUL367iqa6hhbkpWVhVatWjXoOVyBvY7x0aNHMWjQIIwbNw5vv/222WOSkpKwadMmxW3p6elISkpq+Jtwck11jC3h57jmx9T2t8JgMACA1KY5f44tckAgR25q8uTJQkBAgLB161YhJydH+hGHgM+cOSO88847wv79+4ULFy4IO3fuFB566CEhODhYyMvLk57n1KlTwsGDB4Xnn39e6NChg3Dw4EHh4MGD0mq/goICISIiQnj66aeFI0eOCKtWrRJ8fX2Ff/3rXw55302pqY7xypUrhS+//FI4fvy4cPz4ceHtt98W1Gq1sHz5coe876Zkj2OcnZ0thIWFCX/84x8Vz5Gfny+9ztmzZwVfX1/h5ZdfFo4fPy4sWrRI8PDwEDZu3OiQ992UmuoYL1iwQPjuu++EU6dOCdnZ2cJf/vIXQa1WCz///LND3ndTsscxXr9+vbB8+XIhOztbOHfunLBu3Tqhc+fOwt133y29TnP+HFvCoIrsBoDFnxUrVgiCIAi//fab8MADDwjh4eGCp6enEB0dLTz55JPCiRMnFM9z3333WXyec+fOSW0OHTok3HPPPYJGoxFat24tzJs3rwnfqeM01TFeuXKl0LlzZ8HX11fQarVC3759hTVr1jTxu3UMexzj119/3eJzxMbGKl5ry5YtQs+ePQUvLy+hbdu20mu4u6Y6xv/4xz+Edu3aCd7e3kJwcLAwcOBAYfPmzU38bh3DHsd48+bNQlJSkhAQECB4e3sL7du3F2bOnGlWoqK5fo4tUQmCINh79IuIiIiouWFOFREREZEdMKgiIiIisgMGVURERER2wKCKiIiIyA4YVBERERHZAYMqIiIiIjtgUEVERERkBwyqiIiIiOyAQRURNWt/+tOfMGrUKIe9/tNPP4133nnHpraPP/443nvvvUbuERHVFyuqE5HbUqlUNd7/+uuvY/r06RAEAYGBgU3TKZlDhw5h8ODBuHDhAvz9/Wttf+TIEQwYMADnzp1DQEBAE/SQiOqCQRURua3c3Fzp8urVqzF79mycPHlSus3f39+mYKaxPPvss2jRogWWLl1q82PuvPNO/OlPf8KUKVMasWdEVB+c/iMitxUZGSn9BAQEQKVSKW7z9/c3m/4bOHAgXnjhBUybNg1BQUGIiIjAp59+iqKiIowfPx4tW7ZEQkICNmzYoHitI0eO4IEHHoC/vz8iIiLw9NNP49q1a1b7ptfr8c033+Chhx5S3L548WK0b98e3t7eiIiIwCOPPKK4/6GHHsKqVasafnCIyO4YVBERmfj8888RGhqKvXv34oUXXsDkyZPx6KOP4q677sKBAwcwbNgwPP300yguLgYAFBQUYPDgwejVqxf279+PjRs3Ii8vD4899pjV1zh8+DAKCwvRp08f6bb9+/fjxRdfxJw5c3Dy5Els3LgRAwYMUDyub9++2Lt3L8rKyhrnzRNRvTGoIiIy0aNHD7z66qto3749UlNT4e3tjdDQUEycOBHt27fH7Nmzcf36dRw+fBgA8PHHH6NXr15455130KlTJ/Tq1QvLly/Hli1b8Ouvv1p8jQsXLsDDwwPh4eHSbRcvXoSfnx9+97vfITY2Fr169cKLL76oeFxUVBTKy8sVU5tE5BwYVBERmejevbt02cPDAyEhIejWrZt0W0REBAAgPz8fgDHhfMuWLVKOlr+/Pzp16gQAOHPmjMXXKCkpgUajUSTTDx06FLGxsWjbti2efvppfPHFF9JomMjHxwcAzG4nIsdjUEVEZMLT01NxXaVSKW4TAyGDwQAAuH37Nh566CFkZWUpfk6dOmU2fScKDQ1FcXExysvLpdtatmyJAwcO4KuvvkKrVq0we/Zs9OjRAwUFBVKbGzduAADCwsLs8l6JyH4YVBERNVDv3r1x9OhRxMXFISEhQfHj5+dn8TE9e/YEABw7dkxxe4sWLTBkyBDMnz8fhw8fxvnz57F582bp/iNHjiA6OhqhoaGN9n6IqH4YVBERNdCUKVNw48YNPPHEE9i3bx/OnDmDtLQ0jB8/Hnq93uJjwsLC0Lt3b+zYsUO6bd26dfjwww+RlZWFCxcu4N///jcMBgM6duwotdm+fTuGDRvW6O+JiOqOQRURUQNFRUVh586d0Ov1GDZsGLp164Zp06YhMDAQarX1P7PPPvssvvjiC+l6YGAg/vvf/2Lw4MHo3Lkzli5diq+++gpdu3YFAJSWluK7777DxIkTG/09EVHdsfgnEZGDlJSUoGPHjli9ejWSkpJqbb9kyRKsXbsWP/30UxP0jojqiiNVREQO4uPjg3//+981FgmV8/T0xEcffdTIvSKi+uJIFREREZEdcKSKiIiIyA4YVBERERHZAYMqIiIiIjtgUEVERERkBwyqiIiIiOyAQRURERGRHTCoIiIiIrIDBlVEREREdsCgioiIiMgO/h85ahs3P+AJpAAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "lc2 = Lightcurve(tess_bjds_clean,pdcsap_fluxes_corrected,err = pdcsap_fluxes_error_corrected)\n", + "lc2.plot()" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.12.3" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +}