{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# How to remove common systematics using basis vectors (CBVs)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Cotrending Basis Vectors are used to remove long term trends in light curves. They are built from the most common trends observed in each channel. You can read more about CBVs in [Demystifying Kepler Data](https://arxiv.org/pdf/1207.3093.pdf). They can be used to clean lightcurves of common trends experience by all targets.\n", "\n", "Below is an example workflow to remove CBVs from an light curve in Kepler; [Tabby's star](https://en.wikipedia.org/wiki/KIC_8462852)." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Downloading URL https://mast.stsci.edu/api/v0/download/file?uri=mast:Kepler/url/missions/kepler/lightcurves/0084/008462852/kplr008462852-2011073133259_llc.fits to ./mastDownload/Kepler/kplr008462852_lc_Q111111111111111111/kplr008462852-2011073133259_llc.fits ... [Done]\n" ] } ], "source": [ "from lightkurve import KeplerLightCurveFile\n", "lcf = KeplerLightCurveFile.from_archive(8462852, quarter=8)\n", "lc = lcf.SAP_FLUX \n", "\n", "#Remove nans and outliers\n", "lc = lc.remove_nans()\n", "\n", "lc.plot(linestyle='-');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The light curve above has a significant dip due to the astrophysical variability from Tabby's star. However there are other long term trends and small ramps due to heating in the spacecraft, which we would like to remove. We can use CBV's to do this. Let's import the CBV corrector from `lightkurve`." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "from lightkurve.correctors import KeplerCBVCorrector" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can use the `plot_cbvs` function to take a look at the trends we are going to remove. Let's look at just the first 8 trends." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWQAAAEWCAYAAABR8e3qAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzsnXd4HMXZwH9zRady6sVqluQu2ZZtLHdsYwPGdNM7Ni10EiAQCPBRUoAQWiC0BBJTHHq3CRhi3HEv4F5lrGYVq0vX5/tjT8eddCedpNPpZO/vefTodnfKu7sz7868886MkFKioqKiotL7aHpbABUVFRUVBVUhq6ioqIQIqkJWUVFRCRFUhayioqISIqgKWUVFRSVEUBWyioqKSohwTClkIUSOEEIKIXS9LUtPIIRIFkLsFkKE93A+UggxONBhVVRCHSHEOiHEiN7Kv1cUshCiUAhxqtvxZUKIaiHESb0hT3sIIYYKIT4UQlQKIWqFED8KIe4WQmjdPgANzr8jQoiXhRB6Z9xvhBB/8JLmHCFEWRc+HPcD/5ZSmpzpLHXmP7pV+p85z8/o4m33OkKIkc7nVymE8Oos7yw3O4UQjUKI/UKIaZ1If4oQYokQot75Xr8UQgzvRPz5Qog/+Rs+mOkJIeY53/8NnYgjhBD3CiH2CiGahRA/CyEeF0KEBUImH3k+KoR4p6fS95GnQQjxqrOuHnW+9wy3IE8DbepssOj1FrIQYh7wEnCWlHJZL8rRRjkKIQYBa4HDQL6UMha4GBgHRLsFjZNSGoF8YDJwm/P8fOBqIYRolfTVwAIppa0T8hmAeUDrArwHmOsWLhGYBFT4m3aIYgU+AK73dlEIMQv4C3AtyruYDhzwJ2EhxGRgMfA5kA4MALYCq4QQA7steS8ihIgHfg9s72TUF4AbUcpSNHAGcCrwXkAFdBKIXmwX0/gNSh0dhfLua4AX3a5/AcwUQqR1V74uIaUM+h9QiPKybwQqgXGtrk8CVqM8rK3ADLdrS4EngHVALUqlSnBeywEkoHMexwJvAKVAMfAnQOu8dg2wCngOOAr8yYuc7wCL2rkPj/yc554C/uH8HeGUcbrb9XjABIx2Hp8J7ADqnTLe4yOv6cC+VueWAg8DRW73dTvwivPcDOc5A/A8UOL8ex4wuKVzr/MZlQDXOe9psFseN7iFvQZY6XbsHtaA0sL4GTgCvApEdLOsDFaKaZvzq4Hru5jmCuBlL+f/C7zl7T7d79VZbq2ABWgAvnQr1793vs9q4N9AeFfT6+K9vQrc2vq9dRBnCGAHJrQ63x8wAyf5WRb+htJ4qQM2AtPcrj0KfIRSp+qc5dTivO8GYGsg6qwf9/oK8JTb8VnA7lZhvgXmdafcdvWvN1vItwB/BE6RUm5oOensPixCeREJwD3Ax0KIZLe4c1EURzpgQ/m6e+NN5/XBwAnAaYB7N24iSqsqBfizl/inohQivxBCpAOzgTUAUspmlFbeXLdglwC7pJRbncdvADdJKaOBkcASH8nnA7u9nC9BUQCnOY/nAm+1CvMgykduDDAamAA85JT5dJRnPAulYp5K1/kLMNSZz2AgA+WD0QYhRJYQokYIkdXZTIQQWpReSrIQYp8QokgI8XchRIQfcSOBKcCHXi5/gPIc2kVK+Q9gAUrFNkopz3G7fCVKGRiE8iwe6mZ67rIvFELc7ysdIcQElOfyakd5tuIUoEhKua6VXIdRyvJpXmO1ZT3Ku08A/gN82Gq8Yw5KfYpDKfePA+8777nF7NatOiuEmCqEqGlHxjeAE4UQ6c6ycCXKh9idnSj1JOj0pkKehfKyf2p1/irgKynlV1JKh5TyW2ADSkuyhbellNuklI3A/wGXOCupCyFEP5Ru151SykYpZTnKl/Uyt2AlUsoXpZQ2p/JsTSLKl7ojKp2FoBhoxFOJvwlc7KYs5jrPtWAFhgshYqSU1VLKTT7yiENpRXvjLWCuEGIYivnkh1bXrwT+IKUsl1JWAI+hmE1A+UD82+15PtrRzXrDaZb5FXCXlPKolLIepcJd5i28lPJnKWWclPLnLmTXD9ADFwHTUJTACfih/FCUhQbv77UUSOqCPO78XUp5WEp5FEVhXN7N9FxIKc+WUj7p7Zqz/L8M3CGldHQy6SR8l/NSINnHtdbyvSOlrHLWp2dQekzD3IL8IKX8zFmv29S3QNRZKeVKKWVcO2LuQenBFaO01PNoazOuR6lvQac3FfLNKC2I11vZWLNRFFhNyx8wFXC36Rx2+30IpXK2rkjZzvOlbum8hvJl9ZaON6pa5euLJGchiETpUn3dckFKuRLFnjvHaZ8cj9J6aOFClI/NISHEMqd90xvVeNqt3fkEOBm4A3jby/V0lOfUwiHnuZZrrZ9nV0hGuf+Nbs/7a/yszJ2kpSK+KKUslVJWAs/i+dH2RTXgwPt7TUMxoXWH1s8y3VfAAHMr8KOXj7E/VOK7nKfh53iEEOK3zkHWWuf7j8WzXnZU3wJRZzviFSAcpbEVhVJ3WreQo1HMpUGnNxVyOUpXaRrKl72Fwygt4Di3v6hWLYP+br+zUFqZrSvSYRT7V5JbOjFSSneXlo6WuvsORWH6hfOLPR+YLIRwL4hvobSMrwYWSymPuMVZL6Wcg1LoPkPpNnvjR5QPmLd8m1AK1S14V8glKIW9hSznOVBaQK2fpzuNKIq2hVQf8lWiKMoRbs87ViqDnQFFSlmNYiPv9FKFzl7ADyiDs625BPif87fHfQshWt+3r7xbP8uW59zV9PzlFOB8p/dOGYpZ5hkhxN/9iLsE6O80ebgQQvRHMXW1DLb7LAtOD5f7UJ5hvLOBUgu4N7Za32Pr40DU2Y4YDcx39uLMKAN6E1rV1zyUsaug06teFlLKEpSW3elCiOecp98BzhFCzBaKa1m4EGKGECLTLepVQojhThvQH4CPpJT2VmmXooykPyOEiBFCaIQQg0TnXOseAaYIIf7aUoGEEIOFEO8IIdp0aZyeEFcDZSit6xbeQrHN/go3c4UQIkwIcaUQIlZKaUXpQnnchxvrgLhWLjruPIAy+FLo5dq7wENC8WNOQrHrtnhrfABc4/Y8H2kVdwtwgRAiUij+xl69Hpzd5H8CzwkhUpz3lyGEmO1D3nZxumGFA2HO43Dn823h38AdQogUp2fBncBCt/jtuf3dD8wTQvxaCBEthIgXisvZZBRzDigVcoQQYoxTjkdbpXEE8OaRcZsQIlMIkYDyTt7vZnr+cg2KIhnj/NvgvJcHAYQQ1wghCr1FlFLuQbE7LxBCTHLWuxHAxyiDp985g7ZXFqJRbL8VgE4I8TAQ04HMR4AcIYTGKUcg6mxHrEcx78UKxT31VhQzSCW46nABysBe0Ol1tzfnwMHJwEVCiCecx3NQCnMFylfzXjxlfRulJVqG0v34tY/k56JU6JZR74/wzwTRItt+lEqaA2wXQtSiFNINeNpza4QQDSgFbDJwrpRSuqVTiFKwo1Dcaty5GigUQtShmHGu8iGLBeWefV0vcZpHvPEnp8w/otjsNznPIaX8L4rXxRJgH20HFZ9DGQ0/gvIxWeAjD1BaSPuANc77+Q5PG6ILoQzqNQjfg3rZKC3uFvetZjwHNf+IUrn2oAzCbMY5yOP8eDfQdnwCcJmRZgMXoPQQDqHYoKdKKfc6w+xB+dh/B+wFWj/bN1Bs/zVCiM/czv8HRakccP61POeupudCCPFfIcQDPu6pRkpZ1vKH8s7qpJS1ziD9UcxpvrgdeB3lQ90EbEN5Lue52aTbKwvfoPTS9jjjmejYvNAysFolhGgZO+lWnRVCTHPWRV/c45RtL4p+ORM43+36ucBSZ2Mx6Ag3vdEnEEIsBd6RUr7e27IEG6F4mqwATvAxCKkCCCGuQjGd/D7I+RaiuIV911HYYCOEWAz8Rkq508/wfwDOQ3HZ7BV7am8ghFiL4k65rTfyPyanGB+rOD0kcntbjlBHShnU2V99ASmlv65rLeEfFkKUo9iQv+4o/LGClHJib+avKmQVFRWvSCn9GRBUCSB9zmShoqKicqzS64N6KioqKioKnTJZJCUlyZycnB4SRcFms6HT9S1LSl+UGfqm3KrMwaMvyh2qMm/cuLFSStnhJKlOSZ6Tk8OGDRs6DtgNKisrSUrq7uzV4NIXZYa+Kbcqc/Doi3KHqsxCCL9mwKomCxUVFZUQQVXIKioqKiGCqpBVVFRUQgRVIauoqKiECKpCVlFRUQkRVIWsoqKiEiKoCllFRUUlRFAVskqnkVKyefNmmpvVBedUVAKJqpBVOk1dXR2ff/45y5cv721RVFSOKVSFrNJprFYrAJs3b8ZisfSyNCoqxw6qQlbpNA6HsoGEyWRi27ZeWcdbReWYRFXIKn5RWVmJxWJh37592O2/bPu3fv161CVcVVQCQ+gti6TS6zQ2NnLw4EFycnJYuHAhu3bt8hpu6NCh7Nmzh6KiIvr37+81jIqKiv+oCvk4RkrJ8uXLSU9PZ+fOnRiNRhISEvj000/9in/KKadQWFjI+vXrVYWsohIAVIV8HFNXV8eSJa03mfafpKQkRo8ezaZNm5g9ezZRUVEBlE5FpXfZsWMHer2eIUOGBC3PPquQ7XY7Go0GIURvi9InkVLy9ttvdxju1FNPpX///qxfv77NAJ5Wq2X8+PGsX7+eTZs2MW3atJ4SV8UPCgsLCQ8PJzU1FSklDocDrVbrEcZisWAymYiJieklKUOfmpoaTCYTH3zwgevc6NGjOf/883s875BRyHV1dR7HNpuNkpISUlJS2LNnD5988gm33XYbERERPP/889hsNgwGAwUFBYwZM4bExMQ2hU/FN2+++Wa71y+//HKGDRvmOs7Ozuaiiy6iubkZk8nk2pUhJSXFtXHBiSeeiEbT+XFiu92OlDIkd3oIZRobGzly5Ag1NTXk5+czf/78NmGuvfZasrKy2Lx5M1988UWb6wMGDKCsrIy77rqLsLCwIEgd2ixevJjVq1e3Ob9161ZOPvlkYmNjezT/Tm1yOm7cOBnoHUOklHz99desXbsWgPPOO4/PPvus0+lMnjyZ2bNnB1Q2fwnVXQpAcVHbu3cvTU1NLF++nOrqaq/hZs+ezTfffAPA8OHDueSSS/zOY8eOHXzwwQdtlLi/PPfcc9TW1jJz5kzy8vIwGo0cOXKEAQMGtAkbys/aF4GWuampia1bt7reVyBISkri+uuv54svvmDnzp1MmzaN7OxsBg8eHLA8gkF3nnVhYaHXjxrAxIkTOeOMM7oslxBio5RyXIfhelMhm0wmnnzyyYClB0oX22azMXDgQPr37x8Uk0YoK4n33nvPp5cEwL333uuy/dpsNqxWK2FhYZ3qbdjtdp5//nlSUlK4+uqr/Y5nsVh46aWXqK2tbTfctGnTmDp1KkVFRZhMJuLi4sjIyPA7n94mEOXD4XBgMploaGjg5ZdfDpBk7RMZGck999yDyWQiMjIyKHl2l64860cffdTntTPPPJMxY8Z0u/fgr0Lu1T5ioJUxwHfffQfA0qVLAbjttttITu5wb8GQYPPmzeTm5hIREeEzTGFhIWvWrCEvL4/ExEQSEhL4z3/+Q1FRkSvMGWecQV5eHps2bWpXGRcUFHgMxOl0ui6ZDbRaLQUFBSxdupSqqioSExM7jGO1Wnn88cf9Sn/FihWsWLHC49zvfvc7rFZrj3che4O6ujqMRqPL/FNUVMTrr7/ud3ytVuvhK97C0KFDGT9+PAMHDqS8vJzXXnut3XSamppYsGAB+/fv58EHH0Sv13fuRkKM0tLSDu/ZHffGSrAIegvZbDYzf/58SktLOxXv4YcfRkrpsh27s3r1ahYvXtxu/Pvvv5/w8PBOy+sPXfkqHz16lIqKCrZs2cLOnTs9rl1xxRX85z//8Tg3atQofvzxx27L2sLYsWOZMmVKwFr29fX1PPfcc0ycONEv09HChQsDumHu+eefz+jRowOWXiDxt3w0Nzfzz3/+k6NHj/qVblZWFpGRkTQ3NzN27FhGjRrVqR5hbW0tBw4cYPTo0Tz77LM0NDT4DBseHs4ZZ5xBREQEycnJxMTEhOSYjbdn3dme+AMPPBBwe3rImizWrFnD119/7XFu2rRprFixgujoaObOnUtycjJ79+6ltraWE044ASFEh4NFzc3N/OUvf2k3zMMPP9ylQaeO6IxCrqurY8eOHW2eQU9z2mmnsX37dqZMmcKIESOAwJtaPvzwQ/bv38/dd9/dboG2WCwereOCggLi4uIYNGgQBoOBF198kfDwcIYMGcJPP/3kd/7x8fHMnDkTIQRGo9Flg24p40IILBYLQoigtvbae86HDh2isbERnU7X5iPcHueeey5jx44NlIgebN68mc8//7zDcD1Vn1pwOBzs2LHD1WMUQpCdnU1zczNGoxGz2UxlZaWH+aqyspI1a9ZQWFhIZWVlp/Lr168fl19+OXFxcQG9DwhhhdzaXvPII48ghEBKiRCiW0pi27Zt/O9//+OCCy4gJSWFJ554ok2Ym266ibS0tC6l7wt/ZXY4HPzhD38IaN7+8Nvf/pbo6Og25wOtkFsGRTpSFitWrOB///uf67h1mWgpC6C4IJWUlHi4IAWCiRMnkpiYyIQJEwKarjfcn3NjYyM//vgjubm5bNy4kZUrV3Y6vTvuuMMvs1B3ac+26i5LZGQkUkoiIyOpra3FZrMRExPj+uhVVVUBtCtzy+BzVlYWZWVlNDY2snbtWg4fPtyhDPfffz/19fW8//77WK3WDsckIiMjufzyy3njjTdc5y699FLy8vI6zKurhKRCllLy2GOPAcqo/uTJk9uECbSSePzxx32uSBYoG5G7zA6HA4vFwgsvvEBTU5OrULvfuzuRkZGuAZuwsLAOV0+76aabqKysJD8/n++++46VK1dy7733UllZSUlJCRaLhcjISBYtWsT06dM5+eST/ZI7EEgpeeWVV9Bqtdx4440+u8/uFX3u3LkMHDjQ7zwqKys5ePAgS5Ysobm5mTlz5vjVmvOFRqNh9uzZ6PV6hg8f3iNmLffn7I+SA0hLS8NsNrcxXwTTrllYWEhiYiJfffVVG7OaPwwcOJDw8HB27NgBwF133cXixYvZvn07AKeffjoxMTEB/9j6Ij09ncmTJ5Ofnw8o7yUhIaFHW/kthJxCbt06bGkZtybQSsLhcPDZZ595tb8WFBRwzjnndDuPyspKoqOjefXVV326lXnDV+V0OBw0NDRQWVmJVqtl8eLFjBgxgrFjxwZUYfSEd8j69etZtGgR119/vdfp1O7miosvvthlPvEXbzJbLBaqq6tZunRplxSHO3fddVfABgpbFl4aOHAgUVFRHZrUAG6//fY291dYWEhZWRlZWVmkp6cHRDZ/cH/WJpOJtWvXsmfPHoqLi4MmQyC58847e8Qc4Q8hp5C///57li1bBuCzskLPuZCZzWY+/vhj9uzZ43G+ZSZadnZ2p9OsqKhgyZIlDB8+nI8//tjveGFhYdxxxx1ezQjBpCeetdls5plnniE3N5cLLrigzfWnn37aNXjk66PcHv7IbLfbWbhwIRkZGeTn53uYriZMmMC6des6zGfWrFmccMIJhIeH+9WCslgsNDQ0kJCQACjPwZvJzBvjxo3jrLPOCrlZp96etd1up6ysDIDi4mK++uqrHpXh2muvJTs7G5vNRlNTE++//36HHwS9Xs9NN91EYmJiyDzTkHN7a1HGQK8sRGMwGLjiiiuwWq38+c9/dp1vcZNr4corr2TAgAFoNJp2K2JTUxMvvfQSgN+tMqPR6JpteKxiMBgYM2YMGzdu9Lq+RYsyvvbaa3ussmi1WubMmeM6bqmcFosFo9FIbm4uhw4d8iiTrfn222/59ttvARgzZgzFxcXU1dUxb948ioqKXC6C+fn5mM1mXnrpJerr6zsl57XXXktTU1OP2i4DjVardQ2iZWRkEB8fz+eff+56r0ajkTPOOAOLxdJpU1JmZiZFRUWkpaWRlJRESUkJWVlZgOKSGRMTw69+9StAqX86nY5vv/2WAwcOMHnyZAYPHsy2bduYOnVqAO84uATdD7m3p8fq9XoefPBBiouLvc7KWbBgget3eno6s2bNIjY2lvj4eIQQ2Gw21q1b59PNLicnh/PPP5+YmBiklDQ1NfH0009z3nnnMWbMmJ66rZBi3LhxrFu3rt31LbrSI+kqLYO4LZ4fAwcOZODAgR4KuWUpUW9s2bLF9fsf//iHx7VPPvmkU7JoNBquueYaHA5HUJ9BTzFkyBDuuecedu7cyeHDh5k+fbrLrDZo0CC0Wi2HDh1y2YlvuukmUlJSsNvtbNiwwVWPLrnkEoYPH+53vi0TVc466yyP87m5uYG4rV4jKCaLkpISV0H+/e9/38aP2J1gznrz170HlALQ1NTk9VpaWhqlpaWMGTOG8847L5Ai9ig9+aznz59PdXU1v/nNb1w9DXd/UH8Ht1oTaJltNhtarRYhBMuWLeP7778PWNoAv/71r3nhhRcAZXryOeec02cUcSCfta/1ShwOB1artV2d0BlCddZsSJksWpRxampqwB58IDjhhBMYOXIkJpOJo0ePsmjRIsrLy72G9aWMY2JiuPDCC0OyEPQmEyZM4IMPPmDv3r2u9S1aWpqBGEgNFO4K4qSTTuKkk04CFBfKjz76qNvpJyQkkJubS3FxMbfffnu30+ur+JpEotFoQkon9DZBtR94G+TpbfR6PXq9nujoaG699Va2bdvGgQMH2LRpU4dxb775ZlJTUzvtgH48MGzYMKKjo1m3bp1LIbcMBnWma9pbjBw5kpEjRyKlpKSkhKioKL744gsyMzPJzc1l9+7dpKamsmnTJvbu3UtiYiJTp04lPz+frVu38uWXXzJp0iRAGThWP9gq/hBUhdwX1pRoqYjjx4/n4MGDLF68mJycHCZNmsSyZcsoLS3lxhtvJCUlpdft4aFM6/UtEhISXC3kvjSoKYRwDWLNnTvXdb7F/SwvL8+16WuLaaagoICCgoIgS6pyLBAUjRIREUFGRkbIuKD4Q1paGmlpaUyZMsV1rq8PGASbgoICli9fzoYNG3psmm8oEIyJBSrHB0FRyPPmzTsmV+VSaZ/o6Gjy8vLYvHkzI0eO7G1xVFRCnqB82lNTU/tUN1UlcIwfPx6TyeRas+HCCy/sZYlUVEIXta+l0qNkZ2eTkpLimjyjDm6pqPhGVcgqPYoQgvHjx7uOA73SXl+kflkRpU90PH1b5fhDdRNQ6XFGjRrFokWLeluMXkFKiWy2AdC4voyGVSVYyxoBKLpf2QVFlxRBzOwcLIfriZ6egdaobjZ6vKIqZJUex2AwcPHFF/cpL5tAIB2ShlXFmBcdpIi9PsPZKps5ukAx6TQsL0KfaST5hnykzaEq5+MMVSGrBIXOLrPZV3CY7eCQNO+sImJkElVv7cC8r6ZbaVqLGih59AcAYs8eSERuAg1rSomZlYXGoFbZYxn17aqodAHpkNhrzJQ9td51rvoD74sTtSDCdWhj9NjKm9HGhpF07UhsR01UvbXDZ5zahQeoXXgAgIaVyrKT8RcOIWp8KubCWsKyYhCa46vn0V2k3QEaEZI9tmNeIUurA2tZI/p0I0Lb8QuQNgcOsx1rWSMNq0qIGJFIVEE/13V7nYWmnyqQZjtIiBqfSmcWaGrBYbZT9vR69CmRROQno4nUETkq9GcyqsDRj/bQtOFIp+IYhsSReEUemggd9noL2mjFFKFPjSLzyWlIh0Ra7NQvL6J+SfvbFlV/vJfqjz1NIMapGRgGxRI+OB6hV8fqWyPtkuYfKzj6/m7PCwJoqb5aQeIVuVS9rZiPUm4bQ1j/4K5ZHvQ99Tqiu6s11a8sBruDxk3l2I78siBQ9IxMYk8f4BHWWtGELi4cW1Uz5sI6wofEUfbXtvenjTMQPiSeqImplP99S5vrLcSekYNxWiZIibRJNAbPBVVM+2po2ngEa0UT1qK2O/yGD4sn7vwhOJqshKUbO3vrnaarz9phtlPx6lYMg+OImZWNJix4uw/31mpetqMmj9awL8JyYki+Ph9rRRP61CgcDRaqDh4hZXSW33nVLT2M5VAdkWOSOfru7o4jeEFj1ONosKJPiyLu3EEYBnR+YlaorpzWHh7bqZlsmPfVUPVO93aRATBOzyTuzAEdB/RBSK321hNIu6T875uxVTZjGByHaWf7W6fXLy0ialI6ujhlZSmHycaRZzb6lZe9xkzj+jIa15e1G672v4XU/rfQdawx6okYnkjUhFSat1VSv7So3fim3dWUPam4Q0WMSiJqfCqGwXHYq0xo4wwIXWi0fGxVzVhLG7GWNmLaUUX8hUMwDOydrXF6GimdLSsvijFqfCpRE1PRxodjOVRH+NB41ztq+aBqYwxoMiI7lWfMjF82cIgcnYK91owmSk/DD6WED4njyPMdL3zlaLACYC1tpOK1H3ultdebFD20Emyd77n6omF5EbGzc/zqZXeHPqWQHU1Wjn68F2mxY977y8BJR8q4hbIn1xF33mCETlD9ke9R70DhaLDSuK6MxnW+FXnSr/JpXFdG89YKj/PNP1bS/KPnKnJCryHp2hFIB+jiDOiSemf2o6NRqewxs7Jp3HiEin/8RNSkNGJPz0ET3qeKVLvYqpq99pgA0h+ehCZS7zqOGN5zu0BrY5VGRPQ0ZZGj9EcnI/QahFZR/tIuKX6w/d2ry1/aQur9E1wNkmMRh8WONNlp+qnSpzI2Ts0gYlQS1R/vJXpqBpoIpby2tKKjJqaijQ6jeXsV1tJGj7jS5kD4WEY0UIS0yUJKCVJRbLX/PYi9zox5f/tbfLcQWdAP21ETloP+hQ/LicFSWAdAwuXDiBydAvziK9pC+LB4kq4dSePmcrSROsKHJVBRUk6MLZyKl7f6lZfGqEfoNKTeM87VonKYbK6RdX9JuiGf8MFdb5l2tkvqMNlo3HAER5OV+iWH6Xd3Ado4A3WLD9GwqhhtjIH4CwYTPiyhyzIFWubu0Obd5yWgSwgn9syBnWopBUtm6ZAIjcC0v4bKf/5E3DkDqfnygEeY9D9MaWNikg5J0+ZyIk9IAbtEOhxoDLqQNFlIKXHUW7BVmkArsBTWoYnU0bytEtNu3xsM9/ttAdK0DAbIAAAgAElEQVRsJyzTey9B2h1Ii8OloEEZf7JVm5BmO7qUiG55uITcJqf+4l4I6pYepu7rwg7jJP0qH02EDvP+GoyT09t07aXVTu03h1yj1O5ojHpS7y7waO2442iygk7Trp3UXebmnVXo042YdlbRuLEc6+Ff9lmLOTWLmFN97xbRvKMKW1UzhpxYyl/ybat2J/bMAURPz0TaHTiabZ3yW/VW4Zp+qsQwMBZtlPI8HCYbDatKqPv2UJv47pXbfKiO6o/3YCtvJnJsCnFnD/T5TLtDTykJaXNQ+00hUQX9qPliP+YDnh/ytIcmdtknuDcVm73BQumf1rY5n/LrE9Aa9ZQ+7n3GYMafp1JVXUVCtDJI2NueHFarlaLDh2mqasBhtfsdTxOpbzOW05OEh4eTmZmJXu9Z9vu8QjYfqKXiHz/6DKdNCCfh0mEYsmP8Ttu0t5rK+dvBLkFA2gMTXaPdgZDZG1JK7EdN6BI7b15o3llF1ZuKS1T0jEyaNpVjr7O0G8c4JZ24cwfRtLUcXWKER4vg6Pu7sTdaMe/xbEnokiOwVTS7BoIAoian0fhDabt5ZT7puV+etDmo+9/P1C87jCZST9ycwUTmB1YRBVK5SSlp/qkSfVpUu+MJcRcMxjih61O+Q6GlefS9XTRtqeg4oBfCsmNIum4Ept3VmHZXk3Dx0ABL5x0ppcs1befXG4kfkEJ8ZKxf7mr6dGPQPyJSSqqqqqivr2fAAM8BwD6tkGNMBp/eDGm/n+CyqYUKPVnhHE1W0AgP26y02jn60d42dmev6ASGgXFtlHB3SbhsGJFjUrxes5Q0UP3RHqwljUTkJxF37qCAfPggcM/a3mil9E9rfnF58kHGH6cg9N1rYYWCQgZoWFNCzWf7u51O8o35PTqIK612Kt/a4Ron0kTpqTrFwND+g9ooY2HQIs12hE6DCNMio7ToDYHvmfmLlJJdu3a12Um8z3lZ2GrNNKwoxrSpDFOTZ5ekO93Fvo63br/Qa0m8PJfmsSlU/Xt7+wnYZLeVsSYmDOwOYs8eREReAsKgbbeVEpZuJOW2MdQvL6buu0OY99cQe84gIsck96gzvsNip+Th1YSPSCThkmEeXVVpc4AQ2KtNlD3dcaMi7YEJaGNC68PfXYyT0tEaw7y6gaU/NpmSR/wbw6j4x08kXTcSaXMQnpvQrZaow2yj6s0dWEoakSab9zCNVsDgKjv61CgQuAY13bFarV2WJRB0t3wHpYVsr7dgOVxP49pSTLuriRidTOLlv+y+IaWk+PdtR4lT7xuPLj680/kFm1BoAUmrnbr/HaZ+afuTCgyD44idnYM+w0jZB9uJGZyEYWCcy8c2sqAfTRuPEJ6XQOQJKd2erGItb6L6oz1Yfq4nekZ/osb1QxMTprRoulCRfT3r+pXFrhlt7sTMyiZ8aLzfNvnwvARiTs0mLCNwfuChUD5a4zDbsZU3oU0Id40XSLsDW0UzusRwbDVmju4rw/p523EXb7TXY/KG+VAdFa/4NwgOUH2ukWFZg30OyrVgtVrb2G8DxXXXXcfChQtJSUlh27ZtPsPt3Lmzyy3koCjk1qPV/hB71gCip2V2Ol5vEEoVrsUzpXbRARpWlbjO97tnHFqj3sP04eHRYncAokf8LKVDUvbUeuw1Zo/zqfeOc9nWHSZbuy1va0UTlsI6mgfoXDJX/Gsb+pRIHCZbp2fOARhPyiT2tGyvLa1AEkrlozN4TLJotlHyWMct6ORbRhPWP7rdj239imJqF7X9eLYm4/GpWH6uw/JzPYcMlYwYN6rDd9WTCnn58uUYjUbmzp3bYwo5ZEwWLehO6Ufy5Jzj1kTRXYQQICDunEHEnTNImY2YGNFha7QnlZLQCOwNbQcjjzy/iYRLh7XpQuszjFiLf5nJqEuJxF5rVqarA+7Ta7pijokcm0LEyKQe9R0+1tBE6Mh4Yiq28iaOPKdMTIkcm0LTpnKPcC2t3qhJacSdOwihEVgrmjjyzEaiT+6P0Gup+6bQax66pAgSLs/FtLOK6JOzEBqBIScWQ04s2p31Pf7h7Ijp06dTWFjYo3kERSGHD4vHeqSJpGtHuF6mN9L/MIWjddWqMg4g+uTOzRLrKZKvz6fiNU+vGWl1eLVnuitjAFt5U5sw3tCnRZF8y2iEVoPQCizFDZS/uNl1PRADdMczQgj0/aI8vGuEQevVG6dxTSmNazzPt16jI2JEIrFnDmjjgdSRueixL7ezo6TO6zV3z4zOMDw9hkfO6f0VCYOikJOu/WWDy/SHJ2FvtKJLjMDRaKVxfRlh/aMxDI4LydWXVAKDYUAsaQ9NxLyvhqPvdW19hja4LQyT9sBEtDGeH/KwDGMb1zyVwBI/ZzDxcwZjr7NQ+nhbf2dftPjPq3gSdJOFJlLv8hzQRocRc7L/i66o9G20xjAixyjmgsb1ZdR8rrhgpT86GU24Dikl1tJGzHtraN5VhdBrsRysRVodrjR0p6dh+1ppeWU+oSrbUEEbE0b6Y1Ow15io/bqwzXIGupQI7DVmpMVB8i2jOzV/oDXttWR70oYcDELOhqxy7CN0GoyT04malObRKxJCEJZuJCzdSPRJnq0nW7UJodVQbakjdcbgYIus4gcagxZNvyiS5o1A2hxUf7oPw8BYj+VrVdonNJYPUzku6YyJShcf3sYkoRK6CJ2GhIuHHlPK+PLLL2fy5Mns3r2bzMxM3njjjYDnEZQW8qdP/YG0wcOYdMGlwchORUVFJeC8++67PZ5HUFrIBzauY9X7b3O0pP31gFVUVFSOZ4Jqsvj3XTcHMzsVFRWVPkVQFPKljz7p+l1T1v4KYioqKirHK0FRyCk5A12/3/jNr7A0++for9J9HA47az55n0Uv/BXpcCAdjo4jqaio9ApBGdQLi4jkzgWf8vyV5wPw4jWXcNMrb2JMUKeudpYDm9cTn5pOfFpGm2tVxYexNDWRkJHJrlXL+e71lzyu71q1zPV7/JyLGDZzdo/Lq6Ki4j9B80PW6vQUnHUeGxd9BsBrt8zjln8uIDKm87vhHi/UVZQz/57bmHPPg3z0p4c8rs3969+JMEbz3qP3UV9ZgcPu/y4KAOs//4j1n3/EqTfcytrPPqS+UllbOT4tA70hnJjkZObc81AHqfhPyyJW6mxMFRXfBHVQ76SrrvM4fuVXV7Lyvbdpqq2hrqJc7U67YbNY+Oft12E1NbdRxgBv3Xs7r90yj9ojZZ1Wxu589/rLLmUMUF1aTHnhfvatX8MXzz6OuamR3T+swGoydTkPgPceuY+/XX1Bt9LoK1hNJiymZtfxzqXfUl1WgpSSwzt+wm6zUldZzv6N/k81Vul9Dh8+zMyZM8nLy2PEiBH87W9/C3geQZ2pJzQafvv+Qp659GzXubWfvs/aT993HZ//yF/65FKFgWTle2+x9tMPup3O0IknMvPamzDGJ2C1mFn0t79ScNYcPnjs937F37t2NXvXrnYdX/Xk3+g3YFCXZCnZrWxF5XDY0WiOvQV+mmprqDhUSFJWNq/edDUAOn0YNquyyt26Dxd4jZc+bDj9BgxiyITJRMUnkpDe1hSlEhrodDqeeeYZxo4dS319PQUFBcyaNYvhw4cHLI9e2cJJSsmuVcv46sWnfYa5c8GnaHV9Y056V9e7Pbh5A588+ajrOK5fGjVHvHuhnHn7bxk0fhLblixm2Tv/YvCEKez5YQXRSclc9uhfCIuMJDzKiHQ4EJr2Oz4HNq0nddAQflq1HHPNUWKT+1FdWoS5qYnRp57Be4/8DrvN++4NAL9551N0nVgvwG6z8fyV5wEw9bK5TDz/Er/jtqY31xaWUmI1m1y9BUOUkc1ff8muVcsoP9j9rZHcGXXK6Uy74hoQ8NJ1lzHzmhvJzBvpMUDe04TSOs7e1hj2RjDXspgzZw633347s2bN8jgf8gvU++JoSTFVxT+z6avPKdrRdsHn4dNP5ozb7g5Yfj2FvwW3tryM5ro6PvjDA1jN/pkARs86k1NvuLW7InqlPbkddjvPXTHH6zWtXs81z7xCXL/UDvOoKj7Mf//+DEcO7AMgI3c4lz32VI/I3FOYGhsIC49g67dfseTfrwU179bc+PJ89m1YQ9bI0YRHGRFCYG5uIjohiV2rl/PNK8+7wg4smMDMub8iLrVrG7R25VnXlpdhamhAaDSuj4epsYF969eQNXIU0YnJFO34ibSheZ36qHsouf/eD2U/eQ3nkA40oguW2NR8OOPJjsM5KSwsZPr06Wzbto2YGM+FkvrsAvUJ6RkkpGcwZPxkAEr27OTd/7vXdX3H8iXs37iWSx5+gvLCA4w46ZQ+OSh0cPMGti9fwu7Vy/2OM+eeh8jIHU64sf0ta3oKjVbLb99fiMNuRwhBVdHPvHnv7QDYrVbe+PUN/Pb9hT7jSynZuvgrlr3zL3QGA+fe/QBfPPs4xbt2YLWY0YeF7n51Drudwq2b+PQvj3Up/nm/+z+qig5TX1VJ3tSTCEtI9lBs7r2YN++5jcrDh/xO+x+3XuN32AMb13Fg4zpScgZRXrifCXMuUlrdPujKwKuluQlTYwNL/v0PhkyYzNcvP+d33BbOvfsBwiIiWffFRwydOIWouASSswcQmxKa62A0NDRw4YUX8vzzz7dRxt0lJHed1tosLHv7DfZv8Bz0GDCmgMkXXUHakGE9KkNn8dWS2LF8CRU/F7Lhy098xk0ZMIjzf/dwr7gAdqUF9N4jv6N4l2IPvuvdz73agxtrqvnmlec5uGUjOaPHMvuWOzHGJ3iMHbSnzAMtsz9YzSZK9uyirqKcxa+94He87FEncPad92E1mdi5cinZ+WPoN9BzNbr2ZJYOB011tZga6jE1NiKEIK5fKus+/5CNiz53hTMmJtFQVdm1m2vF5IuuQAiBRqslfVgesSmpRCcm8exl5wBw7XOvkZCewQd/+j8sjfUkZvRnwNjxLPrbUwwaN4nh02dSunc3WSNH88kTjwREJm/owgxMuvAybBYL+oxshgwaqPjSS4nVZCIsIoLoxGR0YWGuj0hPmyysVitnn302s2fP5u67vffe+6zJwhvuhfdfd95IdWlJmzB3/edzNNrQGRhyl9liamb36hXtVurEzCymXXEN/UfkExYe4TNcT9MV5SaldFXcG1+ZT3SCZ/y963/g29dexGoyMf2qaxkz+2xXZak8fIg377kNgFvfeJeILrT+e0IhH/ppi1dPFm/c8s8FRETHdKoVGSiZD2//kS+f/wsxSSlkDMtj03+/8Lief8ps+g/PJzImjoiYGN6+79fdzjNQpA3NpXTPri7Hn3rbvWRnpLcbJjopGY1Oj8NmJdxoRCCwWsxKHZMShMBqMrkUuOSX3orNbEKrD0Oj9b2vo5SSefPmkZCQwPPPP+81DBzDCllKyYd/fJDD23/0Gvbk627mhNlne70WLEyNDdQ1NKC1Wlm+4F8c2LTeZ9jLHnuKjNzAjch2l64qij1rVvLlc09y+R//SvpQpeBZTM18P/+fbPt+MSk5gzjzjntIzOzfJu6Kd99k3Wcfkjl8JCNnzCL3xJPQ6vy3nAVKuVmamyjdt6dDRXzTq2/RWH2UhPRM9OFd2wG9p1r1dpuNxuqj6AwGtDodhsiotnm7fQR7ikHjJjH5ost55/7fcOoNtzFsyjQ+ePR+Kn4u5MaX5xOdqNx7U10thVs2Mnz6ya64K997G1NDHbtXryA+I5Oxp5/DlsVfUbxru0ce/ijkQBIZE4shKsrVINRoNewpKmXatGnk5+ejcZqcHn/8cc4880yPuMesQm6hqa4W6XDw7sP3UnukzOPamb++lyETpgB0apCgq2xc9DnL3nkD6XBgjE+gofpoh3G8tSRDga4qirJ9e1jw4N2c97v/Y1DBREr27OS/f3+WmvIyJsy5iCkXX+HTQ6aptoZXbrzKdTxo3CTOu9f/CSiBUG7bl/2vQ1vnkIlTOPOOewNSpkLJW+GLZx/3cGXsLikDBnH1k4H3x5UOB5WHD/Hlc09SXVrMtDvuIytNGUQ2REVhaWpCaARSErT5C/0GDOrQgwn68KCev7TM5rvhhddpbqjn5esvd1376oW/un6PPXMOQsCMub/qETnWf/Exyxf823XsSxmPOOkUUgYMZuSMUzhyYF9IKuPuoHUqKavZzKoPFrD2k/eJTkri0keeIDNvZLtxI2PjPI73b1jDM5eezdBJUznnrvt7TGZQfKBfnHcJNovZ47zQaPjN259gbmoEKTFEGTvVau9LnHv3A9SWlxGbkkpTbQ1HS4tJzMxCOhzUlJWQPjSPjx9/mMKtmxhx6hmc/qtfWteW5ibsdjtWUzPRCUl+KaeuIjQakrMHcN3zilfLzp07SR00xGf4sv17AYiIjcNmakZKiUarxdL8ywQdIQSdaYC2RqJs49iT9IkWcmukw8F7j9xHyZ62Oxa3kJw9gPPu/T9iklO6JIeUkn0b1rB71XJ2/7DCrzhDJk5h3NkXkD40t0t5BpuuttyOlhTx77tuJjzKiKmxgeHTZnLydTd77TJ7o7m+jpdvuMLndX14BDe/+iZhEW13zO6MzNLhwOFwIISgaOc2Pvzjg23CXPjAH8gZPdav9LpKKLWQ/cHhsOOwO6iprQ0Zubvqh2w1mdAZDC67sJRS8bGXUmlYOH3LwyIisdusWE0mGo4eRWg1WE0mopOS0enD0BsMfo9bHfMt5NYIjYbL//hXjwGm1lQcOsg/b/9lqnZETCwX3PcI/QYOdn3Zy/bv5ciBvTjsdk44/RxXeiNnzmLb99+2K8Plf/wru39YyZjZZ2HXhYVMwQ0GesMvttSz77yPYZM7t9loRHQMd7/7BTVHStm77gdW/Ge+x3WrqZkXr7mErPwxXPzQn/xOVzocVBX9zOGd28gYNpzVHy5o46njzsTzL+lxZdwX0Wi0x8xsytZ2fyGEpxlKCNeHX6vTozXqXa6mUsqgu9n2yRayO9LhoOLnwm6PKA+fNpMdK75vN4zOYGD8ORcy/pwLPF50X2sBtdAdufeuW03qoKGuAZvu0FB9lNduntthuFGnns6gE2dSfXAvBWed5zpvamygqugw7z18bzuxFTKHj8RmNnPFn58NWmU7HstHoAnFmXq+OOYH9fyhZM8uEjOz+Pu1XZ+W640RM07l9FvubDdMKBXczhBKcksp2bb0W8r27mH/xrU01lQHPI9fv/WRR+s+WITSc+4MoST38aKQ+6TJwhstdtuWCQd7167mi2cf5+Trbqa6tBiHzcbu1SswNTZ0mJYxPoFhJ57EiJNOISkzq0flVlEQQpA/8zTyZ57GLG6ncMtGPg7ApIPRp53FhDkXEZOUHAApVVR6lmNGIbdmyMQpbWaDnXrDbUgpKT+4n+TsAS4j/d51q9Hq9PQbOBh9eHivTtZQUcgZU+B6f9LhwNTYwOt33OBzt5mr//ICiZlZFO/aQdrQYSE9NVulb2IymZg+fTpmsxmbzcZFF13EY491bXq9L45ZhewLIUSbqa0tfswqoYnQaIiIjuGO+R+4utEWUzNWk4nmulqSsnJcYbNGjuo9QVWOaQwGA0uWLMFoNGK1Wpk6dSpnnHEGkyZNClgex51CVjk2CAuPICw8gqi4+N4WReU4QQiB0WgEFFu11WoN+MCwqpBVVFT6FH9Z9xd2HfW+LkZXXdVyE3K5b8J9HYaz2+0UFBSwb98+brvtNiZOnNjpvNojqFs4qaioqPRltFotW7ZsoaioiHXr1rFtW9t13LuD2kJWUTmOsZrt6A19axJIey3ZYLm9xcXFMWPGDL7++mtGjmx/uYDOoCpkFZXjEJvFzmu/Xub1Wv5JGRwta6J4dzXxaVGcMKs/qQNjqSpuJGtEAvVHTSSmG4Msce9TUVGBXq8nLi6O5uZmvvvuO+67r2MzR2dQFbKKynHAyg/3svV/hzHGG2ioNrcb9qdlxa7f1aWNLHmrrb02OiGci+4fh8MusVnsaPUaohOCP+kmmJSWljJv3jzsdjsOh4NLLrmEs88O7PK/qkJWUTnGaKw1s/rjfexZd6TNtY6Usb/UHzXx79+tbHP+xIsGM+bUY3My1ahRo9i8eXOP5qEq5FZYLXbsVgfhUe3boSwmG021FjYtPsSA8TFs+mwbZQfruOz/JhAWrj5Wld6hscbM/PtX+RX2+memucq5lBKbxcEPn+2nZE8NI0/KID5Hy2eP7/CIM/qU/mz932Gfaa76aB+rPlI2tI2MDWPqRUNoqDaz+pN9DJ3Yjz1rj5CQHsXMq3JJzopGq1P9Ctw5ZtayCBQLHllDzZEmbntV2dXA4ZAIQGgUV5qXbl7ScSICZl6Vy/ATf9nhoKnOgsMuqS5rpH9eQk+I3ml6+1l3BVVm75ibrLx+d/vLxGYMi+OcO8ag1WlwOCQaTfvuYb7kLtpdzefPbSZvSho7V5d2S+4W+g9P4KxbR/lU0OpaFscA0iFBdLyL7vL39lBV3EDJ3hrXOb8Ur8+M4fu3d7Hhq0Lqq0xtLp940WAGjE7mp++LSOpvJCYpnPQh6gQHFf9x2B1otBqkQzL//lU01Vm8hht3Zg4TzhnQpg50pIzbI3NYvKvBcvJcT8Vz5GAd37y+zWu5b4/DO47y6u1LCYvQcfbto+k3IKZbMvZVjlmF7LA7eOW2pa7jSx8aT1JmNFazHYvJhlanYemC3VQerqe2otl3Qn4w65bBfPvKvjbnfRVK926dO/0GxDB2djaZw+Kxmu1ExRloqDYRFWcI+rqsKqFHY42ZNZ/tZ9easnbDTT5/ECOmZ2CICH717jcghrl/nkLFz/VsW1GMtEt2ri5l4rkDSB8ST2OtGUuzjaULdnuNb2m28clfNwIw8dyBjJnVH52+b7nldYdjxmRxtKSRdV8eICxSh6nBysGtgdkyvTVjT89myLh+NNWZSRsUh96g9ZDZarHzDx/uRF0lrl8kw6emM3BMMsY4A0IrcNgc6MK6V1DV7n9w6I7MzQ0W/nVP28Ezb5xyTR65k9K6lI83gvGsLSYbQiN4+6EfaPbRyp84ZyCOmGpyMgeh0QoS0qLQaL2bNlSTRQiwftFB1n15sFtp5OQnkpmXwOiTlZ2SD22v4odP9jH14iEY48MxNVpJHRjrFsO7H6Y+TMttr55MbUUzG78uJGt4It/8cxvDp6Uzac5AwqP07NtYjrnRyrJ39/glW82RJlZ/vI/VH7dtVedNSWPKBYPZvrKYIwfr0GgE2fmJ7PqhjIT0KDRawZBx/YiMCVMiCDDGh7sGcYLFhq8OsvaLg9z80gy0PiqTioLNamfv+nKaGyz88Mn+DsOPPT2b0Sf3/+Ud9yFaBsCve2oqP++o4ssXtrYJs/bzA4y/SjHpOeySyqIG4lMj0YVpj7meY59qIdvtDop2VrPw71vR6jXYrf4plCsfm4RGq3yF3bn0ofEkZhi7/VK72pKwWx3UlDehN2hpqDHz6dObuiVHV5h+2VAObKlg3Jk5ZAxVCr3d5mDXD6UkZUYjNJCQHkVtRTNl+2sZOiEVvUFLdVkj9VUmzM02LM02zE025XfLf+c5i0n531ijuFtNu3Qoo2ZmdlneY7GFvGtNKVqdhrIDtcT3i/T5oR4+NZ1hE1NJHxLn9Xqg6Y1nbW62seL9PexeU0ZyVjQVP9cDMP6qeHL6D2oTPjHD6DEQGIwWst1uZ9y4cWRkZLBw4cI214/pHUOsFjsCaKgxs+DhNe3GTcmJobywjonnDmD0qVnoW3Xpzc027FZHwFsSgSq4ysaq5Qw8Idk1Em412zE3WhEawVsP/LJ9e3RieKcHTnoajUYQFqnDEKEjLEKHwe23+2i80AgmnzeI0adk+ux6+qIvK+SDP1ZSebiemMRwaitNrF/oX69u6MR+zLp2RA9L2ZbeetZ2uwNrs51wox6Lycb2FSVoU+oYOmgYdVVtx3vCInSER+kxROiw2W09rpCfffZZNmzYQF1dXcAVclBMFq/fvRxzk415T0yhodrcquvvneZ6C/+61z/bGYAxwcDF97d/v4YIHYTw2vNCCIaM7+c61mgEhgida3Dm1ldmusIBrP3yABsWFbrCj5yewbblyiyrq/44Gb1By3t/WkdznaVTPQp/mX7ZUDKGxmOIVJSuLkzjs7cx8+pcXr5F2bNQOiSrP9nH6k/ammAARp6UwfTLhvbJ7qiUkp+WFhMVG0ZKTgwHNlew8ZtCmuusXUpv7uNTjvkZcK3RajVojcqHOixcxwmzsti5cyfhRj36cC0IqCr6Zecfi7NHBhDbLxx6UB8XFRWxaNEiHnzwQZ599tmApx8UhWxuUh7Wm79XWnin3ziSQWNT2oRrrDGz9N8HOLL/x3bTy8lPpLHWwsiTMhgyvh86vW9FcCzR+h4nnjOQsadlY2q0uirtSVcM8whz3VNTvablcEh2bThERLiRsv21bPrmkNdwA0YnMe3Sob+0zgUgQR+uZeT0DJd/tj+yTz5/ED982rFNdNuyYrY5p+/Oum44KTkxrP3iAPH9IskqMNJQbcJislPxcz3L39vDnDvHYDXZyRimmFyklDRUm3HYJbUVTWQNT/RLxs5itzvQCIHQCCoO19NUZ2Hhi21toB2RPzOTn74vYvplQ9m2vJj+uQlMvWRID0jc99HqNJQ9/jimnbtASqytxkFqnWE0WkFnNIIhL5fUBx7oMNydd97JU089RX19fecE95OgKOTTbhjB4te3u46//sc21/n+eQnYrQ6a6i188Of1PtO46o+TiE4I73QX91hHb9B2abUujUaQMtBIUlISA0YlMfn8tvY5d1r8TrvD2NnZjDm1P58+s4myA3V+xfn2X54zxdYvahvmwyc6NqNd8+SJaPWaDmdgtofDIVnz6X52rCrBarIrk4Y0QvF39wO9QYvVbOeUeXkMHJOM1WInKlbZamr6pUMByJ/Rdfv68YQAEAJ9mBaHQ2K3/aKY7TYHdpvyvJ3tB1cc99+dZeHChaSkpFBQUMDSpU/RiFcAACAASURBVEu7Ib1vgqKQh4zrR8bQ+DZz392VtDdOmZdHVLyBjKHxx6WT+LGIRqvhwt8ppqUdK0tIzo4muX80NeVNVJc2Unaglk3f/BzwfN2nE0fEhDH+zBxSsmPYvqKYpjoLp9800uXvun1FMZVFDUy5cDCFWyvR6ARJmUbe+b+2Yxj+KOMZ1w5ARwTDJqZ6nA/rBT/hYwFfLdnaiiZXb7wjouIMro+hv6xatYovvviCr776CpPJRF1dHVdddRXvvPNOp9Jpj6AO6kkp2fLtYX74dB++ss0eE8ewcRkMGdfPe4AQpC8ONEHoy11zpIkFj6zxbNo4CY/SY2q0otVpCIvQkjEsnn0byj3CxKdGUl3mfVPUnmDoxH6MnJZB6sBYJGBqsBIZExbyz9kXoSR3Z6ZOV5f4N9ErIS0KKelSD3Pp0qU8/fTTfXNQrwUhBCeclsWoUzLZtqyYlR/sdV2bfMEg+uXEYEiwh0whUOld4vpFtjGVHNxVTL/0ZCJjwrBZ7R6zuKZebMZucxCT+MvIrZSShS9u5ecdRwMmV0p2NCfPy6O+0sSBLRWceNFgDJGephABfdIv+FggJTuG2vImzM3tt5aPljZ6HMf1i+z1hcGCkvvP199A46pVDNuyGU14OFqthtEn92dwQQqRMWEeg1WVlT0zw07l2CA6yeBSdK2n1HrrggohOOfXYzzONdVZWPj3rVjNdmqOeLagY5LCmXlVLmEROravLGHHihIA8k5MI29KOmmDfvEQSkw3kjNKbTyEIrEpkV7PWy126iqbvXoctS4LoHhmmZttxCRFeIw/zJgxgxkzZgRM3haCopAbVyn2u8Y1a4h2u4nO2nBUVAJBZEwYlzwwHlA8JXAoPuqtW7Qp2THMvDK3N0RU6SH0YVoS0qIwNVrRaATmZhumBt8uiS2t7LrKZsLCtT3uVBAUl4XBy5W1HYpuviUY2amo+I1Wq0Gr16jmheMIIQQRxjAMkXpiEiNIyY4hId2IMb59f+9geHgFpYWsT0kBnQ5sNmxVVegSe8YvVEVFRaUr6PQadPowIqL1WJpt6A1aLCY7pkYr+jAthm64S3aGoDn19n/5JQD2njgVhym0pvyqqKiogNJ6NkTq0WgVn/W4lEii4gzo9MFRlUFTyFEnnuj6vXvMCUhH8FYaU1FRUekLBE0hC62WrPnzXce7ho/AUlQUrOxVVFRUQp6gzkOOmjSR2PPOcx3vP3UWhVdeha2qKphiqKioqHSJnJwc8vPzGTNmDOPGdTjPo9MEfWGItCceZ/Cypa7j5o0b2XviVOwNDdgbGnxHVFFRUQkBvv/+e7Zs2UJPLEUcdIUshEDfrx/DfvRcFWvPuPHsGTee6muvQ9r8m4+uoqKicizRa/MENWFh5O3ayc5czznftr172TUyn6RbbyHp5psRYap/qIqKyi+s+GAPlYe996allF1aijepv5FplwztMJwQgtNOOw0hBDfddBM33nhjp/Nqj15fbip35w5spaVUv/supl27aVyxAoDKl1+h8uVXPMIdD2seH4tIh4PmLVvQpfQjLDOjt8XpFUw7dmDIzUVo2u+USocDy8GDaKKi0CUnYz18GG18PNrYjjd1UOl5Vq1aRXp6OuXl5cyaNYvc3FymT58esPR7XSELIdCnp5Py298CUPr992i3bKXqtdc8wu3KG442Lg5tYiKOujoiJ0wg4Zp5ROTn94bYKh0gpcS0bTt1ixZR99//YjtyhKipU8l6/Z+9LZrf2BsasFdWEpaTg3Q4MO/aRf133xE5YQJh/ftTu+grKtx2jUi86SaPchs2YACWg8o2TRVe0tdnZ2E91LWlRg15edjKykj+zW8Iy8mhadNGzHv24qirw7xvH7bycvSZmVjdPJlERASyWVkJLaKggPQ//4nGH34g7qKLEM5tj2zV1VgP/X975x0eR3Xu/8/ZvqtereqO5S5jG9sQIBAwmIABAxcwJWASCCEECBBCLgnYSQjEJNehXAIkgVwgoSfBiUOJTTH8wKbZBhe5F1mW1bWSdldbZs7vjxnJktWllbSSz+d55tnZM2fOvHt25jvv6ftxz5jR7nVjgc482f5eUy8nJweAzMxMFi1axCeffBJVQY7ZNfU0r5dIRQV7zlvYaXzHqFGE9u8n7frrSfvuDVjj218Nuj+JpWkKe0J/2B3cvZu6VavwrlpliI3dTvwppxA+dAgZDjNuVdvpCnvCQOb10dVpxxqpr72Kf+kytDovucuXEy4pIVRSQqSsnNRrr8WemTFgVYo9mX6zvwTZ5/Oh6zoJCQn4fD7mz5/Pvffey4IFC7q0NSan3+wJ1qQkrElJTCraRt2bb1Jy2w/bjRfabyw9VPWHP1D1B8P7KtjwBRZ3DC+eN8wIl5Tg/fe/qVv1b4JFRWCx4Jk7h/Trrydh/nysSUkUf/9mGtasoWjmLCZ+8flgm0zo4EGEzYYlLo6GtWupf/s/1L/1VlTSTvnW1dQ8+1yrsKRFi6hbtQoZCrWJb01LQzuq62f2Aw9Q+pOftPFyB5Lqiy9p3t932eWtjtW0mJR95J+fwX388SAlwukctlWLZWVlLFq0CIBIJMIVV1zRRoz7Ssx6yJ0RLimh4pFH8b7+epfppS5ZQtp1SwiXlhKpqiLh9NOjZWozx6KHHKmqou6NN6lbtYrAhg0AuAsLSTz3XBLPWYAtI6NV/IrHH6fykUcBo7g96rlne1Wa6cpmPRhE2O3IYLD5paz7fITLytF9DQQ2bKDsVw/0+LoAOct/jWvaNIJFRVhTUnBOmIAtNdW4RmMjh358N2nXLcFdWNilzVLXkaEQFlfvFzDVfT4q/vdx0m+6CWG1IINBNK8X+8iRoOsI65HpSbX6eoTd3ny9SE0NMhhE9/nYc+55vbahPRxjxiA1jUhlJbm//Q2uyZOxeDzUr16D5/gZ2PPyELae+YKx4CF3l754yENSkFsiIxFkKESouJiav/yV2pdf7vKc0S+9iGv69Ki9yY81Qda8XnaedjoyEMA5YYIhwud+E0dex+vBSSkpmjS5VVj8mWeQu3w5Fo8HzetFq2/AmpyMxeNG9/tp3LyZgzd9H93vJ/nSS5HBRryvryT5ssvwrlzZXB/ab9hsFHz+GcJmayVuPWWo3h8HX3iRlHFjcRcWEi4uRvN6cYwdi8XlQjgc7D7nm4QP9H65rfHvrKFx+3YscXG4Z8zg0B13kPbdG3GMzMeamNgqrhLkdohFQT6a4O7dlC9/iNDBg4R2d77CccpVVzHiv3/SZct3VwzVB66ndmu1tWCxECkvZ895C8m6715SFi/u/vkNDeyYfUJvTI06ub9bQfjgQSxJSVjcHuK/firWhATAeMmHDx/u9AXTE4b7/RHcuxffBx9gz8vj4E3fj87FrVbsebmE9x9g1F//yn6365gQ5JitQ+4tznHjyH/yiebvoeJiKh55lOD27djz82lYs6b5WM3zz1Pz/PPYMjOx5+Qw6rlnm1ubFYZXG9q1iz0XXAjtTAZlHzmyR+lZ4+OZVLStXW+5P4k79RSyly7Flp2Nf/16PLNnd1pkFjZb1MT4WMA5ZgzOMWMAmFS0rTlcRiLoDQ3ogQB1b7yJxe3CM2cOh350F41bt3aUnIGmNfdA2X/FFYT/9zECmoYlLg5H031nsQy7+uphJ8hH48jPJ/eh5a3C9i2+orneEyBSXk6kvJyiadNbxbMkJDBh/bo+e9BDgdCBA6DrBDZvoeqppwju2NHlOb3tciiEYFLRNurffRfd78c9bRq6z4dv/Xq8K1fiKphI+ve/j1ZTg2vSRITNhh4KoXu91EiJZ88e7Hn5zX2apa5DJNKmxV9GIkhdx9IiPG7evF7ZrOg5wmYzuqomJ5N23ZLm8DF/e6153/vPf+L7eB2ysRF34XTKHniw0zR1n4/Gbds6jdM0ztcxalRzqWeoMOwFuT1Gv/BXpKYR2LSJ/Vdc2WE8vb6eoslTmr8XbNqIsNtpeP99wsXFJJ1/Ptbk5IEwOWpIKal68ik0r5fqZ55pt39sR9iysogcPgzAiP/+7zb1fD3l6AZW16RJpF177ZGAFoNILA4HlowMRGVlG1EVFgu00/1K2GwML/9p+JG0cCFJC490bU295ppWx4O7d/e60bGpB1ZL7Lm5ICXWlBSAVh621PVBd76OSUEGYzpQz8yZRhFa06hYsQLvqn8TKS3t8Jztha07y7dsrbfedy+eufNwjBlNpLwCGQ7TuG0rCWec0e9/cuCrzUa1y4hMwuXl6A0N2NLSCGzcSPF3b+x1ujm/+Q0J3zjdKN7b7cOueKiIfZzjxjGpaBtbt2zBPWWKUerRNJCS4K5dRo+acMdr4h1NuKTE+Dx0qNN41sREtLq65u/CYsE1uf+r2Y5ZQW6JsFrJvPNOMu+8s1W4b906grt2E6msoOqJJzs42+Dwsp93fNDs75p+443UvvQSejBI0gXn4y4sbLcbnpQSrbaW2pdeouJ3D5P2nW9T9cc/AWDPzydcXEz+U09SfMN3e/5jO8Axbhw5Dz6APS8Pa3KyEl9FTNHk1Aibrbn+3z11KmA8L2B4u+FwGJvNhgwECO7dCz3otNCSlmIMhvesh8PU+Xx85zvfYfPmzQghePrppznxxBN7+7PaoAS5E+LmzWsuHrsKCgiXHsZZMIHI4TLcMwq7X5SKRNC9Xsp//evmoK4EviVNYgwQLi4G6JMYj/jJ3SQtWkRNKDQkW/8VipYc7TwIIRAeD+4pU9rE1erqkJEI1qQktPp69Lo6hM1GpLq6w/RtaWnYMjIQNhu33norCxYs4NVXXyUUCuH3+6P6W5Qgd5PEc85pE9bUolxZWUlaairev/+D8hUrSL/+O4QOFFPzl78MiG1Z992La8oU9l16Ge7CQvKe+L0xGY2UnfefrawcEPsUilihZbuHLTkZzDYguzlHBYDUNGQkgsXpbHVuXV0da9eu5c/mykcOhwNHlIeOK0GOEsJiIfnii0i++KLmsKyf/RQA3e9HhsNY4uObBbLy97+n4uFHsKamoh31dnZOnEjiOeeQdP5CbJmZbbr3BHfvRtjtWFNTW412a9nlSKEYrrz756co37+n3WNSlwhLz6vbMkeN5fRrjak0hdXariOzZ88eMjIyWLJkCZs2bWLWrFk8/PDDxMXF9fh6HaEEeQCweDxtwtK/9z3Sv/e95u9S00DTujVZi3PcuKjap1AouiYSifDFF1/w6KOPMnfuXG699VYefPBBfvGLX0TtGkqQYwRhtUIfhucqFMcKTZ5se/TnSL28vDzy8vKYO3cuAJdccgkPPth5v+meMvxHPCgUCkUUyMrKIj8/n+3btwOwZs0aJke5K5zykBUKhaKbPProo1x55ZWEQiHGjh3LM888E9X0lSArFApFN5kxY0a/rDbdhKqyUCgUihhBCbJCoVDECEqQFQqFIkZQgqxQKBQxQswLspQSfzi648UVCoUiFomJXhalDaUU1xezrnQd+6r3cVz6cTy+6fFun//EmU/wtdyv9aOFCoVC0f8MqiBXBao47eXT2oT/p+Q/PUrnxtWt5/x1Wp0EtSBnjjyTKelTaIw0clnBZWR4MpBSqqklFQpFj9m+fTuXXXZZ8/c9e/bw85//nNtuuy1q1xg0Qb7mjWv4ovyLTuNkxWVx5cQrmZg2kbFJYymqLmJL1RbSXGm8tvM1NF1je832NucFtSAAqw+sZvWB1QA8+WXb6S4/WvwRCY6htcSLQqEYHAoKCti4cSMAmqaRm5vLokWLonqNQRHkaf/Xei22P5z1B+ZlG/MOd7bSbaYnk1PzTgXg0oJLAQjrYQ77DvPM5md4ZccrPbLjpBdOavV9dOJo9tXt46dzf4oQgsKMQnLic5RoD0O6U1Lyh/389rPf8s2x3+SFohe4cfqNuGwu0txpuG3uVnFDWohlHy8D4N0D71Ifrm8+dvX4q7nIehE7a3eyYPSC1ssGSckrO17BY/fwXvF77K7dzZPznyTTk0lloJL1peuZNWIWae403jnwDnOz5pLsSuZg/UGswkpWXBaHfYfJ8GQQ0SNE9AjxjiMzAJb5yvCFfeyr20dRdRHvHHiHW2bewujE0WTHZ9MQauCryq9IcaaQ4cngzb1vctnEy5p/X0SPUN1YTYorBbtFLQDcxJo1axg3bhyjRo2KarpC9mBG/dmzZ8u+jlJZ/ulyntv6HAD3n3w/5487v9Xx/loyXUpJmb+MjeUb+dHaH/U6nR/O+iEXjr+QVFcqjZFGiqqLqK6txmf1keRMIqwby8lMTJ1IbnxuF6kNLkNteXpv0Mtrm19jxeYVzMmawyeHP2kT59en/JrKQCW7anfx+u7X0WXb1bLbY2bmTOZkz+GJTU90HdlkTNIYXFYX26qPzWlP4+xxJDoSKfUdWfZsQsoEdtTs4Lyx51FcX8xNhTfxbvG7TE2fysTUidgsNsYl93y2wm3btjFp0iQAav+5m9AhX7vxpNQRoud9FRw5cSQv7L5d1113HTNnzuTmm2/u1NYmhBCfSylnd5XugApyWA8z87mZACyZsoTbZ9/eJs5AisTHhz7mxtU3dvuh7Su3HH8LBakFzM2eS1gL81XlV9gtdvwRP0nOJHLjc3FYHSQ6+rZ4aHcZKoJ8+3u385/9PWtXUAwf8hPyuWvkXYwrGEeaK43K17djr5ToUscqbIBECAsC0UNBlkgkAoEjJ75DQT66NBUKhcjJyWHLli2MGDGiTfy+CPKAVllc//b1AByfeXy7YjzQnJhzIpu+tanTOGW+Mg75DlEfqudn/+9nVDd2vNRLVzyy4ZEen/PU/KeYkzUHq6Xt1JyNkUacVuewaKSUUvKnzX/i4S8e7rdrLDtpGUEtyLNbnsUiLBSkFpDiTOHlHS+3ifvjE37MeWPPw2KxkGBPQAiBLnVWfL6CheMWcvHKi1vF/9eif/HJ4U84a9RZJDmTWh2rrKzEGm/FG/TycenH/Gr9rwBYeuJSkl3JvHvgXZaetBSbxYama+yv34/b6qaouoigFmTBmAVIKdlYsZF1pes4Ne9UJqdOZl/dPs7/x/ksP3U5U9Om8tb+t1i5eyVZniweOOUBwnqY13e9Trm/nBsLbyTdnY4/4ufDkg/5vOxzpqZPZcHoBTisDrZVbWNE3Ah0qbP0o6XEO+K5IOcCCkcWEtSCPPnlk2wo38DWqq3t5q3b5iYQCUTpn2pNcX0xmtSoC9ZRF6yD0zqbM7y3U9g2UlK5pdMYBakF2Cw23njjDWbOnNmuGPeVAfOQdalT+GwhAJ9e+Skum6vdeEPBa/vs8GeMiBvB2oNrWTR+EX6vn/T0dDRdaxZOf9hPqa+UrVVbKfOX9bvQnD36bOLsPVu5IBby2hv0cvKLJ3cr7h2z7uC8rPM6tdkX9iEQOK3Odl9iHdFyocxoEwv53Bt6Y7cv7MNpdWIRFiymp+oNeimuL8Yb9LLi8xVtGuIvGHcBb+57k6AWxGaxMSNjBprUONRwiFPyTuHVHa/yu8m/I2tMVtR+W09xWB2MTx6PEILLL7+cs88+myVLlrQbd0hUWWyu3MziVYuZPWI2zyzoeMq6oXjz9sTmxkgjVY1VpLvTsQkbEonNYhRUpJQ8u/VZ/rbzbyQ4EthU0bn33hl/XvBn6kP1xNnjyInPwYKFf+z+B6nOVBaOW4jH7hn0vG6MNHLCX07oNM6J2Sfy6BmP4rQa65sNts29YSjaDLFld3siB0eqE3SpIxBEIhFj1WkkIS1EUAs2lySdNifBSBBd6tSF6mjUGnFYDW87EG7r3QsE+Yn5xNnjml8ufr+f/Px89uzZQ1JSUptzOrI15qosnt36LAAPnPLAQF0yJnHZXB029gkhuGbKNVwz5RrAGDBTHiin1FfKj97vWUPktW9e2+GxX67/ZZuwSamTuGfePYS0ELnxuViEhRWfr+D9g+/jCx9pQIm3xzM7azY3z7iZgtSCHtl0NC3FODsumxWnryA3LpdkV3Kf0lUcOzSVaCwt6o2FEAgELpsLl83VqgqpqfdIqju1V9fzeDxUVVX1weLOGTBB3lK5hSRnEllxg1fsGGpkx2eTHZ9NYUYhC0YvaHNclzphPcw+7z5W7VnFlqot7fY86A7bqrdx1b+v6jJeQ7iB94rf473i99hw9YZm776nNHURA/j7+X9nfMr4XqWjUAwnBkSQqwJVHKg/wC3H3zIQlztmsAgLTquTgtSCVt5qRI+wpWoLk1Mns6VqCyMTRxKIBEhyJLH6wGo+OvQRb+x9o8/Xf3n7y1wx6Yoen6fpGq/ueBWAF859QYmxQmEyIIL81r63AJibPXcgLnfMY7PYKMwwGlBnZM5odezC8Rdy4fgLWX7q8lZ1hFJKttds5829b/Ltad/udDBMSAtx3VvXseLzFczLnsfY5LE9su/13a83709Nn9qjcxWK4cyAzPb2QtELAExOi+6CgIroIYRgYupEbpt1W5cjEx1WBytOW4HH7uHOtXfSGGns0bXu++g+AP52/t96ba9CMRwZEEGuClQxJW1Kr+sbFbFHhieD+0++n501O3no04e6fZ436G3ePy7luP4wTaEYsgyIQj4x/wnyEvIG4lKKAeTk3JNZMmUJz2x5hrnZczlr9FldnvPgJw8CcPecu/vbPIViyDEgHvL0jOmkunrXzUQR2/zg+B8wLX0aSz9aSklDSZfx/7XnX4AxGEChGGqsWLGCKVOmMHXqVBYvXkxjY8+q67oi5lcMUcQ2dqud5acuRyK5a+1dzZMrtUdEjzTvt5yRTKEYCpSUlPDII4/w2WefsXnzZjRN48UXX4zqNZQgK/pMXkIeS09aypcVX/LYhsc6jHeg/gAAC8cuHCjTFIqoEolECAQCRCIR/H4/OTk5UU1ftbIposLZo89mXek6nt78NHOz5nJS7klt4rxz4B0Arpx05UCbpxhGvPHGGxw+fLjdY71dESgrK4tzzjmn0zi5ubnceeedjBw5ErfbzVlnncVZZ3XdbtITlIesiBp3nXAX45PH85MPf0JloLLN8ae+fAow5hFWKIYaNTU1vP766+zdu5dDhw7h8/l4/vnno3oN5SEroobb5uahUx9i8arF3P3B3Tw1/6lWcwwUpBSwvWY7HrtnEK1UDHU682TD4TB2e/+sbLJ69WrGjBlDRkYGABdddBEfffQRV13V9ZQD3UV5yIqoMj5lPHfPuZv1pet5evPTzeFSSvbV7eMbI78xiNYpFL1n5MiRrFu3Dr/fj5SSNWvWtDsDXV9QgqyIOhcddxELRi/gsQ2PsbHcWBSy3F9ObbCW6enTB9k6haJ3zJ07l0suuYSZM2cybdo0dF3nhhtuiOo1VJWFIuoIIbj3xHvZXLmZu9bexSsLX2melLzTKTsDtbB3LQRqIP04yDkeDn4GX74Eu9+FsI/0QE3fjHMmQnYh5M6E1LEw5lRIyAa7u23cprnCh8GKLIrosGzZMpYtW9Z1xF6iBFnRLyQ4Enjo6w9x9b+v5r6P7mueRGhCyoQjkSJB+PB38N6vBs6wYB3s+8DY+kp6ARw3H7QwuJOh5HOYfjlkmC+drGlKzBU9QgnyUEZK44HXzUVa9TAIC1jbadQI+aBqN9SVQEMZHFhPQkMNjJoFdYeM82YtAX8V6BHQNQhUQ8V2+PB/2qZndcBde8HZ8QCPqelTuW3Wbfzms9+woXwDGe6MIxMX/e0Gw/PtLvO+DwULqNbjSc3MBYvN+O2+CkjKN+wP+2HXaijbAvGZkD0DLFYINoDVZrwADqwzt48Mz7i+tOtrd0TldmNrya7VbaK1WnPDYjf+p/QCmHAWeEtAC8GMK2DUScbxTvJUMbxRgjxYaGHwFsOOt+DzP0NFUcdxbS5wJUND+30ve4sTYHeLeZE//WP3T9ZC8Ogs+MY9MONKQ/ja4erJV7OudB0flnxIQYrpOb7/0BExzp4Bi1+AxO51sNcrKyGhhcR5WgzJt7tg+qWdJzDh7G5dp0NCPnDEwfqn4OAn8NUrkJgHjV4I1R8R3A5/gHnsaDEv+lfreJ508LfoOmiPg7APZlwFp94BKWOU9z0MUYIcDaQ0vK/ybbDzbajZDxuj2D8x0hh1Me4SZyLMvs4oflfvgdJNMO1SmHy+4YHW7IX1T8LKH8C638NZv4DxZ7ZJxiIs3H/y/Vyy8hJjtZidq+Fdcwmpb62EsV8f2N/VVxzmQrJzbzC2i7t+iVVWVpKemmJ8qSiCV79tiO2k82HL342SSGIuuFOgbLMRL3DU6uZNy2htfL7re2v+z2HyBcaLwjp8HvHeDvoYSHqyRml7DNgip91l0BdW1CJQux8aa2HrStj4V/CV9+813SnwtVsheSR40sARbzQy1eyD8q2QO9s4FqyHhCyjusDTxWRNTdUZndDnvJYStr4Oq+8zbB13hiHMI6a0vVagEqHrpD0yG4JeuH1bt73iqNo8CPTaZl0DBEQCsOd9WHUH1B8yjuXMhENf9Cy9pHyjVJaYB/OXGfXfzsQO75NYyuu9e/eSkJBAWlpap6Lcn/2Qu0JKSVVVFfX19YwZ03rwU8ytOt1d+uUmiASNetLyrbD7HbA64ZMnjbrS/iJ3tlEnOP5Mo7eAFjZufFeSUd8ZA2/6qOV1JAif/gne/7XRaDbjSjj9HkjMbh3vL/9llCAu+kPXVQv9bfMA0u8265ohtJtfgzU/N8KSRxmORW+Z812qJ19N6sgpYBn83rHhcJiDBw92ObuapmlYre1Xnw0ELpeLvLy8Ni+FY0OQgw1GFym72yhK1h2C/9zbtj6uL6SOhfQJcPLthoCHA0ZDUuYkSBsHzoQhKRLQD0IRqIG1vzGqMqx2OOkW+Notxn8jJSwzV5O+t7rDOucBt3kAiAmbdQ1qD8AHv4UNz/U9vXMeglnXGKW1GHAumoiJvG6H7gpy7FcwBRvgzbujcxNlTISs6Ubf07gMo/iXMdHwJhJzYurGGpK4U+Ds++GE78CaZfD+g0aD3ifQSQAAE6RJREFU5TfugQxzRNMZ9/ZajBV9wGKF1DFwwWPG1h7hgPG8bXgWtv2LxoRRuLb/vf24b/zI2Nrj63fDST8AuycmvOuhROx4yFoYfJUE3lqGe9+avtXbnnA9nHJ7r+ooe0OsvpW7ot/tLv4E3rrH6I1g9xj/8Y92GsLdS4ZiXg9Fm6EduwM18P5yWPd47xJMO85ocCw4p9+cn1jN66HhIWtho87ro0eag9oZL2WQMxPO+BmMPuVIH1RFbJM/B779ttHwt+bnkHdCn8RYMci4U2DBA8Z2NMF62PMevHa90QjZHlU74cXF7R+78AlDqF1Jx/SzPTiCrGvw6hLjQW2PMafCKXfAyBPB5hxY2xTRRQiYcqGxKYYvzgSYtBB+2k73TG+J0V979X0dn/+PG9uGjfsGTFhg6EF6wTFR/THwgly2FX5/YuuwU+6EE74NiTkxW+RQKBS9JCkXTr7N2Frirz7S7/rQBqNBviW73zG2lhScCwc+NvRiyiLInDysPOqBFeSDn8Efzzjy/Z4yY3SVQqE49vCkGt4vGJ9fu/XIMSmh/rAh2JtfMwYmlW+F7auM42sfMjaA/HmQMQFGTMVdUw7peVB4efsTRsU4AyfIWuSIGM+7qf16KIVCoQDD603MNrbj5h8Jl9IQ5h1vwbaVxlD20o1QvA6AuKZ4q24HqRu9qiaea9Rvj/oaTL3IGGTliI9Jz3rgBPmdXxifUy5SYqxQKHqHEMZI0BFTjJ5UTWhh8FXgX/sYngPvGnXPm1+Fw18aGxhVHR/8pv10PWlw6XOQNdVoWBwkBkaQfZXGYIH8eXDJ013HVygUip5gtUNiDv55t+M5z5zO9UyzETFQawwzL91kTJG67Z9tz/dXwZ+/2TY8czLM+x64U2HSef1nv8nACPKr1xldYc5/NCaLCQqFYhjjTjZ6bIxrZ/mwhgr4z89g0wvtn1u+1ZhAy2KHWzZAcn6/mjowgtzohcIrjIp3hUKhiBXiM2DRE8Z2NFIapfv6Q8bYh34WYxgoQf7u+wNyGYVCoYgaQhiCHZ8xYJcc/j2tFQqFYoigBFmhUChiBCXICoVCESMoQVYoFIoYQQmyQqFQxAhKkBUKhSJGiP0VQxSKYYymSzRdEtZ0ghEdbyBMWNNx2aw0RjRCER0hIBjRqW4IsXZnBedOyyYQ1pDAYW8jDqsFm1UQjOjYLILGsI6m6+woa+DLg7V4A2H+a3Y+mi7RpWTroToCYY3sJBeJLjthTaeyIUSi284HOyuYOTKFOKeNxsYAWA+iS4nVIrBZBIfrghyo8lGQlUCcw0ZI07FaBP6QRmNYIxjWSYt3UO0LkZnowiJgY3EtVosACYGwRmFeMhKJx2EjGNFw2aw47RaEEBjrZRirS+u6sXiGzSqwmVNvOm1GPCFAAJYW+0IIAgE/cZ4aLMLotXYkrvFpabHPUedbLAIBlNUFyUl2I5EIBBLDjstm52Oz9q8PqwRZoWiB1x/GH44Yk401Rqjxh6gLhGmM6PiDEd7eWsbErASyklyENcmeigb2VPgYmxFHWNOpC0QIhDVq/CEEsOmgl4kjPAQ1wd5KH3kpbg7WBEh02bBYBLX+cI9tfPbjni9e+tBb29sNd9os2CwCf9gQxkBY42BNgLQ4BwKJ22nDIkTzi6PaFyIY0QlFdKxWgd1iIWIKfV0gTGNYJ9Ftw2oRVNQHCes6xdWtJ6wX1FLpC5ER76QxrOGyWwlGdEA2x5DmSwAgoksimo6UENR0kCCRSAm6NOSyBwsf9ZqLZ+Zh6+fVx5QgK4YUwYhOta9JJDUaGg0BbPLQ/KGW+xECIZ2d5fWENZ2IJtlf7ccqBIfrGslNdhOMaAQjOhYhCIQNj7Qr3ilqu7zYVyVe4p02Et2GgAEkue2kxztIddtxu5xIKZmel8z4zHjCms7IVA9Swv4qPwVZCeSluCmu9pPscZCX4sZmFXgcNqQEp93ClhIv3kCYWaNSSPE4sFktVDUEiXfZyExw4bRZCGk6gZBGZqITp9VKpS+IRQjS4x247EfUxNrkGXYylcFQm5tcSklFZSVpaelIKdFbCLfsSMRbfTfOKfUGiHPamkseFmF4zk5b/9fwKkFWDDoNwQgf767CZbfgDYSp9Yep9Ycorw/y3vYKRiQ6OVzX2MbT6gqLAI/DhhCQFucgPd5JWpyDEYku4pxWCvOTcdqsSClx2iy47FZq/WGSPXYyEpxkJ7lJ9thJdNlx2Y3jpd5GRiQ6cTusSAkJrqYHt+OHNVrCdnpBZo/PSfLYu440TBDCEE/Ds+79nDkZCYO3SpESZMWAE9Z0nlq7h6LD9Ww4UMPBmvaF1mGzEIroJLhs6DrMnzyC8al2nC43Lru1WVzdDituuxW3w4qnxb7DaunUA+wN+ameqKanULRECbICoLkBRQKHagN4A2Eq6oPU+EO8t72ClZsOATBrVAo7yuqxWQQXzcyjICuBghEJTBiRgNvRdQXb2h0V/Pi1Lyn1NpLosjFnTBpnT8nCYbNwynHppMc7SXbbSXTbWxWxmxhqxWiFoicoQR4G6LrEF4qw/XA9DcEIGQlOdpU3EAhp7K5oQNPB7bCwr9LPqq9KcdosZiNKz/l8fw0AHoeV59ftb05HCBiV6mHCiAQmZiUwIcv4HJ0W11ycf+7jfdy7cgtj0uL4xYVTuWx2Po4BqJdTKIYKSpBjgMawhi8YodTbSI0/RESX1PiMT38wQrUvhMUiqPaF8Ie05vrVfZU+6hojPb5eWDNENCPejs16pHj/9YIMgmGdcRlxZCa6SI93kOJxkOxxkOCyYbOIVlUAmi45UO1n++E6th9uYHtZHUWH61m9rQzT4cZhtTAuMx6HzcKm4lrOmJjJo1ccj8ehbj2F4mjUUxFlpJTUByPU+sKU1zdSH4xQ5m0kENbQJVQ1BDlQ7WfLoTqjX6fXiNNd0uOdZCY4SYt3UJifTDCiM2FEfHPPgXlj03DbrVT5QmQluohz2pg1KoUEl61NFUBfi/9Wi2BMehxj0uNYMPVIeGNYY1d5AzvK6tl+uJ6iw/XsLKvnouNzWX7J9H7vy6lQDFWUIHeCpkvqG8PUN0ao9YfZVVGPzWIU90tqAtQ3htlf7ae4qh4dC2V1Qeobw83eYUckuGyMSHSRmejkpHFpJLrtpHgcZCY6cdutOGwWUkyv1CIEiS47TrOVfyjgsluZmpvE1NzBW5tMoRiKHDOCrOuSKl8Ib8Ao9lc1hKj2GVtI06n1hzhYE6DaF6LEbNQKhDQinairRWAU6d1WxmXGccLoVFI8DjxOK+nxTtLjHSS47GQlunDaLdgtFpLcdiwWtYyVQqFoy5AVZF8wQl1jmC0ldcawSwS1gRAV9UEq6oPsKm/gQLWf8voggZBGIKx1mp5FGF2aEl12puclkRpniGlanIMkt50El40af5jpeUnEOWykxDlIcBqjrVTLv0KhiAYxKci+YIRtpXV8vLuKhmCkebCANxCmuMbfYb/VJmwWwYhEF6PTPSS57bgdhsc6MSuBeKeN9HgnKXF2PA6jSmBkqgeXPfp9VhUKhaInDJoge/1hvirxcqDaT3GNn+JqP8U1AQ5UNlATaNvIlRbnYFSah2m5SUzLNTzY4zLjCYR1PA4r4zPjSYt3kJngIllVCygUiiHIgAqy1x/mjx/uYc22craW1h0xwiLITXGTn+LhtONSGJeVgqbDyDQ3U3OSGJ0eh121zCsUimHOgAiylJJXPjvIr97YRq0/TF6Km2tPGs2Zk0YwJiOOrERX88xOqj5WoVAcqwyIIF/65Md8uq+GE0ancNPp4zltQoaqr1UoFIqjGBBBnpabzHnTc7h63ihVt6tQKBQdMCCCfO/CyQNxGYVCoRjSqJYyhUKhiBGUICsUCkWMoARZoVAoYgQlyAqFQhEjKEFWKBSKGEEJskKhUMQISpAVCoUiRlCCrFAoFDGCkLKL5S1aRhaiAtjff+YAkA5U9vM1os1QtBmGpt3K5oFjKNodqzaPklJmdBWpR4I8EAghPpNSzh5sO3rCULQZhqbdyuaBYyjaPRRtbomqslAoFIoYQQmyQqFQxAixKMhPDbYBvWAo2gxD025l88AxFO0eijY3E3N1yAqFQnGsEoseskKhUByTKEFWKBSKGKHfBFkIsU8I8ZUQYqMQ4jMzbKkQosQM2yiE+KYZPloIEWgR/oQZntAibKMQolII8TvzmFMI8ZIQYpcQYr0QYnQs2GweW2ym86UQ4k0hRLoZniqE+I8QYqf5mRJDNl9m2rtFCLG8RXjU87mndpvHpgshPjbt+0oI4TLDZ5nfdwkhHhHm2mCDnddd2Hy/EKJYCNFwVPoxm9dCCI8QYpUQosgMf7A/7Y5iXr8phNhkhj8hhLCa4VG/P6KClLJfNmAfkH5U2FLgznbijgY2dyPNz4FTzf2bgCfM/cuBl2LBZoxVWMqb0gGWA0tb7N9t7t8N/DpGbE4DDgAZ5vf/A87or3zuhd024EugsIW9VnP/E+BEQABvAOfESF53ZvM8IBtoOOqcmM1rwAOcboY5gA9a5PVgP4ud5XWi+SmA14DL++v+iMY2ZKoshBDHAZkYNwLABRjCAfAqcEaTdzTICHOLM+1JBA6Zx1ra/H/AhQNvXruMBXZIKSvM76uBi839WMjns4AvpZSbAKSUVVJKTQiRjfHAfSyNJ+tZjuTpYOd1uzab++uklKXtnBOzeS2l9Esp3zXDQsAXQF6M2N1ZXteZcWwYL5KmXgyDfX+0S38KsgTeFkJ8LoS4oUX4zWbR+OmjigljhBAbhBDvCyFOaSe9xRhv3qYMzQWKAaSUEcCL8WYcVJullGHge8BXGEI8GfiTGX9E04Nofmb20d6o2AzsAiYKo0rDhnFz5pvH+iOfe2r3BEAKId4SQnwhhLirhW0HW5x70AyDwc/rjmzujFjO62aEEMnAQmBNP9odNZuFEG9hlFrrMV4Y0D/3R9/pL9cbyDE/M4FNwKnACIzijwW4H3jajOME0sz9WRh/buJR6W0FZrX4vgXIa/F9d1Mag2kzYMe4UcdheMqPAT8149Uedb2aWMlnjAdsPfAx8Fvg7/2Vz72w+05gL8Y8BR7TxjOAE4DVLdI8BfhnjOR1uzYfld7RVRYxm9ct0rJhVA3dFkPPYnfy2oVRZTG/v+6PaGz95iFLKQ+Zn+XA34E5UsoyaRR/dOAPwBwzTlBKWWXuf27+oROa0hJCFAI281gTBzG9ONOrSwKqY8DmGWbYbmn80y8DJ5mXKDOL2Zif5X2xN4o2I6X8p5RyrpTyRGA7sNO8RNTzuad2mza8L6WslFL6gX8DM83wvBbJ5nGkemhQ87oTmzsjlvO6iaeAnVLK3/Wn3dHOayllI7ASo6oC+uH+iAb9IshCiDghRELTPkYdz+amDDBZBGw242S0aP0cCxwH7GkRdzHwwlGXWQlcY+5fArxjCuBg21wCTBZCNM3sNB/Y1o7N1wCv99beKNuMECLT/EzBaKT5Yzs29zmfe2M38BYwXRgt/Tbg68BWaRQ164UQ88w6y29xJE8HNa87srmLy8RsXpvn/xJDbG/rT7ujZbMQIr6F6NqAbwJF7djc5/sjavSH243RSLTJ3LYA95jhz2HUrX6JkSHZZvjFZrxNGI0FC49Kbw8w8agwF/AKRv3nJ8DYWLEZuBFDhL8E/smRaoI0jOqMneZnagzZ/ALGg7cVsyW6P/K5N3abx64y424GlrcIn22G7caoHmoafTqoed2FzcsxvDrd/Fwa63mNUfqQ5n290dy+EwvPYic2jwA+NeNvAR7FKGlH/f6I1qaGTisUCkWMMGS6vSkUCsVwRwmyQqFQxAhKkBUKhSJGUIKsUCgUMYISZIVCoYgRlCAr+g0hRJo4MjPXYdF6pq6P+umaxwsh/mjuXyuEkEKIM1ocX2SGXWJ+f08IMdvcb5ph7CshxFYhxC+FEE7zWIYQ4s3+sFmhaEIJsqLfkMYkLzOklDOAJ4AVTd+llCd1dX4v+W+M/qZNfIUxsKiJyzH6t3bE6VLKaRijwMZiLgkkjYmXSoUQX4uuuQrFEZQgKwYFYc4FLIQ4TRgTHb0shNghhHhQCHGlEOIT01MdZ8bLEEK8JoT41NzaCKM5umu6NGf9MvkAmCOEsAsh4oHxGIMaOkVK2YAxwOdCIUSqGfwP4Mo+/XCFohOUICtigULgVmAacDUwQUo5B2P49g/MOA9jeNgnYIw4/GM76TSN2muJxJhO9GyMeQxWdtcoaUzduBdjiDnAZxgTGCkU/YJtsA1QKIBPpTkVohBiN/C2Gf4VcLq5fybGHCFN5yQKIRKklPUt0skGKmjLi8AtGPMw3IFRrdFdWs7rWw7k9OBchaJHKEFWxALBFvt6i+86R+5RC3CilDLQSToBjHkVWiGl/EQIMRUISCl3iG7OnW5WgYwGdphBLvMaCkW/oKosFEOFt4Gbm74IIWa0E2cbRh1xe/yEHnjGZn3z48A/pJQ1ZvAE2laJKBRRQwmyYqhwCzDbXC1iK0aDWyuklEVAUtPUjUcde0OaSxAdhY3WHvq7QojNGLOWHQC+2+LY6cCqPvwGhaJT1GxvimGFEOKHQL2Usr1Gv6PjOjGmjJwqpfR2I/5a4IIWHrNCEVWUh6wYbvye1h5vu5iDQTYCj3dTjDOA/1FirOhPlIesUCgUMYLykBUKhSJGUIKsUCgUMYISZIVCoYgRlCArFApFjKAEWaFQKGKE/w/4CX/nmIg/aAAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "KeplerCBVCorrector(lcf).plot_cbvs([1,2,3,4,5,6,7,8]);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We probably don't want to correct all these trends out of Tabby's star. Let's pick just the first 5." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "corr_lc = KeplerCBVCorrector(lcf).correct(cbvs=[1,2,3,4,5])" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can plot the two light curves on the same plot window to see the difference." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEKCAYAAADjDHn2AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3Xl8VOXVwPHfmclKFpYE2YIEERBQCBBRqiKiULQVd8VawaVFbfW1trav1LYuVdu620oX3FDrKypWpZZWcUGqorIICLIrYMKaBBJCyDZz3j/unTBZSIaQycwk5/v55OPc/cxlnDPPcp9HVBVjjDGmMZ5IB2CMMSb6WbIwxhjTJEsWxhhjmmTJwhhjTJMsWRhjjGmSJQtjjDFNsmRhjDGmSZYsjDHGNMmShTHGmCbFRTqAlpKZmanZ2dmRDqOW6upq4uJi6xbHYswQm3FbzK0nFuNurZiXLl1aoKpdm9ovtu5eI7Kzs1myZEmkw6iloKCAzMzMSIdxWGIxZojNuC3m1hOLcbdWzCKyJZT9rBrKGGNMkyxZGGOMaZIlC2OMMU2yZGGMMaZJliyMMcY0yZKFMcaYJlmyMMYY0yRLFsZEieApjkOZ7ti3bx8HVq9GfT58paVUFxay8w/3U7bs85rj655HVQ/7OsZAG3ooz5hIUVV8e/eC38/uRx9l7ytzQj52dxjiKXrmmSM+R3JODp0uuZjtt/+K+F69OPqZp0k4+mj2z3qWxNyRpI0d6yQaVap37iSuWzfE4/z2VJ8P8XqPOAYTXcKWLETkaeC7wC5VPb6B7QI8BpwDlAFXqeoyd9tU4Ffurveo6rPhitOYQ1FV9s2fT3yPHlTvLiDvRz+q2ZY0ZAhVeXn4iosjGGH4HFi+nAPLlwNQlZ/PpgnfrtlW9uSTLXKNuG7dqN65k263307nK74HIjhfCyYahbNkMQt4HHjuENvPBvq7fycBfwFOEpEuwB1ALqDAUhGZq6p7whirMbXseuRRCv/2t0NuL1+9ut46b0YGXW+6kQ4nnkj56tV4O3YkrkcPkgYMQFUp++QTKrdsodMll4DHQ2FhIRkZGWhVFeLxIHFxVOblEZeRgSc5GVWlYt06ypYs5cCyZaSecQZp486gYsMGko4/HlSR+PhmvT9VrfXFXLZsGZ7kZOK6dqXkrbfY9Yf78WZmkDRoMKXvvtusazTGD1Tv3AnAznvvZee999bfSYTud91Jp0susSQSpO6/XZWvinhv8z4Hh0PCWWcpItnAm4coWfwNWKCqL7rL64CxgT9Vva6h/Q4lNzdXmzM2VKWvkkmvTyK/NJ/zjz2f20+6nb0Ve+mc1BlVJSkuqd4/TqgKCgrIyMjAfR+HfXy9/6Gryqj0VZIcn0yiN5HSylLKqstYW7SWoZlDKaksYU3RGgoPFHJyj5OZ/uF0Kn2VjOw2knG9xzGgywDe3vw2yXHJdEnqwvb929lbsZfnVj/HMZ2OYcXuFQD079yfkUeNZG/FXo7tdCzHdTmObh26sWDj22wp2Uyyz8tPTr2NA4W7yN+7lS6pXcnuNwKtrKRsyRIS+vQhvlcvttzyE3xfbaZi3To6f//7VGzcSNknn+AHjv7zDFJGj2bfO+9wYOUXVGzcQNKAAVRu/YbS999HEhJImzCBkjffRJKT0QMHOOZfb5LYr59T715WBh4P+P3OF29pKV27OmOhlX70EYnZ2cT36tXkPa72V7OtdBsb926ka3JX9v13IZ1/+XiD+8ZnZdHn+efwpKVT+dUmPKmpJB5zzGH/OwYc6dg/fvXj8/uI98ajqhSVF/HmV28ytOtQjs84niU7l7Bo2yKuH3Y9pVWlVFRX8PDSh8kvzeeOb93Byt0r2Vuxl8+2f8ZDYx/ijY1vsGzXMpK9ydyQcwPdOnTj/sX3c+nASzmuy3EIwsKNC+nXvR9ZaVkUVxTz6fZP+SDvA3qk9KB/p/7MXjebJTsb/v8wWRIZvC+NpakFjb6viz7080W20GWf0m0vXL7AjwfoctVVdP3pLXgSEg77XkXj2FCFBwrJL81naNeh+NXPPzf9k0eXPcr1Q69nYJeBXPnvKwHI7ZbLmqI17K/aX3NsWkIa8Z54isqLAOf/2VfPfbVZ3zMislRVc5vcL4LJ4k3g96r6obv8LvC/OMkiSVXvcdf/Gjigqg82dq3mJovtpduZ8OqEkPbtntKdU3udypz1jddJx3niqPZX11rXI6UH2/dvByA5LpkD1QcAyEjK4L5T76PKX8WDSx4kPSGdr4u/Zl/VvsN+L9HA61N83tA+sIO2KmuOFkat9fP5sUKXfbCzs3NsYqVy9hJlVbZw8lo/K/oKl/3Xz4D80GNJP/dcOowYTtLw4ZT06sjW/XlsLtnMluItbCnZwuaSzeTty6Naa/9bHbVHOXOFnxfHOvXuNw2/iTOPPpNnVj3Dkp1L6JHSg6lDplJcUcx7W99jZcFKpgyegkc8PLf6OXYd2NXw++0yiG9nf5tHlz1asy47PZvNJZtr7dc1uSu7DxxszTi779l8lP8RJZUlCIKinNrrVD7M/zD0m9FGPPq3agaNu4ie99xzWMe1drLYW76XzSWbSYpLYv2e9eSX5jNn/Rx2lTX82WiuUd1H8dmOz7hi0BXcNuq2Zp0jFpLFv4Df1UkWvwDGAYl1kkWZqj7UwDmmAdMAsrKyRn7++efNjrWwvJBfLfkVm0s3U+4rb/Z5osnxcjR5lTvYG18JQKcyDxl7fWzq6XwhH7/Zz4Rlyprezq+4o4ohuRIeuNBDzyIYuVHpsk958XQP+5MFUUWDfrmk71eSKw9+wQMM+8rPimNat5Pd8Zv9rMsSbn7dz9rewsLjhcwS6L4/nu2pVWzrAhUJB2NMII5eKVn0TjuarJQseqf0pnh/AZ/95wmW9Y/tDoJ9UvuwpdQZRDTJm0S5r5yzep7FO9veqbdvbmYuY3qMIX9/Pot2LmLr/q30T+9P9w7dWbN3DV7x4vP78Hq87Dyws9axWSlZdIjrwPri9Vzc92LmfO38gLp24LVc0OcCiquK2VG2g6KKInql9MLn99GjQw8q/ZWkxqeSFp9WU+IqrSql3FdOvCee+fnzeWnTS/jUx/nZ5/PM+oYb6x+eWc2Ac66kww+uDakxvbi4mI4dOx7u7QzZrgO7+PFHP6agovFSU2Mu73c5AzoO4L3891i1dxW3DriVEb1GsOvALl7f8jpXHnslnRM7U1pViqJ08HbA6znyjgRdu3aN+mQRFdVQoar0VSIIJZUlJHgTSPImsXXfVr4u/poxWWPYWbaTKl8Vi3cspn/n/nRK7ERaVRpJHZP4fNfnDOs6DI94WLRtEWuL1jKi2whyu+Vy76f3sq9yHyd2P5GTup/k/PIQOKn7STUfhL3leykqLyKrJI69c14FoMMlF7Djf26has3aI3pfng4d8JeVARDfsyfxZ08ksaKSlNEn0+Gkk/CmpjZ5Dq2spGr3bjwdOhDXuTNVvirWFa2jT8c+pCWk1dy/BO/B6oP1e9YD0Cu1FynxKZRVlfHJ9k/YuHcjp/Q8hceWPcai7YsAuPr4q3l/6/v1foE3pVdqL7ITe9BzXzzdth2g65oddP0iny7FisfjIfG4gXTIGU7yiBFs+8UvwO+nx+9+x1O9N7FhzwYGdhnIeceex10f30Xevjyy0rK4buh1LNq+iFmrZzH8qOH0TO3J0WlH41MfKfEpnNLzFPp16keFrwJwSpEe8aCqbN+/nSp/FUenHY2IUFBQQJeMLmwp2ULnxM50SurEmsI1JMcl0ye9DwCrClbhUx/b928nOz2b7I7ZrCpYxchuI6nyV5HoTcSvfjwS3iQXuEYkqnP2Ve5DEJ5e9TRPfPFEzfqMEuXat/yc+4P76HT++Y2eI1xxr9+znj8t+xML8hYccp9jOhzNdd0uZNi/1lHyxj9JyhlG8tQrKLrlFwhw1FuvkyrJFP75L5SvXUvq6adTvWM7xW/MJSE7m8rNm+l+113suOMOAJIGD6bnH35P6cKFVBcU4i8tpdsvp+NJTm7We4iFksV3gBtxekOdBPxRVUe5DdxLgRHursuAkapa1Ni1wp0smuNwP6C+0v2AcmDpUvJu+alTL98C4o46iupdu4jr0YPef/0rcRld8HbqhDQwsUo01u0eil/9gJOIHlz0IGv2rWF1wWoykjL405l/YnDG4HrH+EpKOLBiBQc+/5yyZZ9zYOXKWvf5uDVftlpjaizd64BIx+zz+7jhnRtqfkgE/PW5ZL711scNfqYhPHG/tPYl7vn0YHVY+n7l5jf8HL9FiURzfHM/u6Emi3B2nX0Rp5SQKSJ5OD2c4gFU9a/APJxEsRGn6+zV7rYiEfktsNg91d1NJYpoE+h3X/rHP+IZMJC08WfhKy6m+PU3iO92FOXr1lP8j3/U7B/oQni40iZM4Khf/IL4Ht3R8nKkQ4d21Wsk8Gs6KS6J6wddH9KXgTc9ndTTTiP1tNMA0Opqyteto/yLVaR8a3S7un+xyOvxMnPCTACKyosY99IZ+PBz/ZQD8MJwHv2/JMb957Ow/zsu3rG4JlEM+8rPLa/76VBx6P0lPp7eTz1JVV4+Jf/6F6mnnw4i+PYUUfDnvxDXvTtZjz5CQr9+bLlyChVr15L++9+RMXgwW6+9Ft/uArpcew3pEyaw+bLJAHS55hrSvz2BvXNeJe2sM8P+nsNasmhNrVmy0KoqqvfsoXzVKio2bWL3Qw+3+DU6X3klaePOIHnoUPa++ipxXbuSfvbZLX6duiL9y7G5YjFui7llzNv4Jv/70fSa5R/+28dNL3xRqy2jJeP2q59hzw0D4M8zqskscdZnv/wSSYMGNbs7c12tOFNeZEsWsc5XXMyBFSvYNv2X+AoLa9XtHw5JSkLL6zeYp55+Oll/nsH+jz5CKytJOKYficf0bfAcXaZMOezrGtNenHPsdzm733f4vw9n8Puv/sYTZ3uZ9cwwXsj9I4NyxrX49d7a/BYAfXYqx//sLjpecH6zuvPGGksWdfj27WPbbdPrPYjUWKLo+49XScjOBlX8lZVOe4DbgNnUL4PUMWNaImxj2jUR4YrTbuTs4y/i2r+OZ2NP4dIVN9PnszT+cfX7LXYdn9/HLz+4DQTuec5Hpy/azwODliyClPznP+T/5JZ663s/MZP4rCwSevXCX1GBNy3tkOfwpKSEM0RjTCO6dO7Ba9NXse6rxVz832vYkriPCx8ZyTMtlDDmb5lPtfi55i0fg+bOazeJAixZAE6D9JYrvs+BZcsASBlzGr3//OcGe1Z420Fx05hYN/CYE1l69FIufGQkW7oJd874DjPu/OyIzqmq/HLhbXQpUS4ZeOkhq43bqth+AqmF7H744ZpEkf3KKxw9c+Yhu+AZY2JDQlwCr974CQAL+x5gzuw7juh8XxV/RRU+TvlS6Xb9j5o+oI1p98mi4uuvKXziSbyZmRy35kuST6j3SIgxJkYlJqew4PRXALir4h/sLt7e7HP99J/XAfDdFXHEdzuqReKLJe0+WSRkZ9PrT3+k37x/tav6R2Pai4zs43jUexUAP/jL+Gafx7OriP75Su5Lb7ZMYDGm3ScLESF9/Hi86emRDsUYEybDJk7l7O3d+KqHsOSzNw77+B37d7Axs4pB3ygJvXuHIcLo1+6ThTGmfbh57C8B+P0Hdx32sR9v/gCAc/tNatGYYoklC2NMu9Dr5HGMWQXru1by9fY1h3XsirUfkFyhDBn9nTBFF/0sWRhj2o3/OfU2VIS7/u8Hh3Xcqq8Wkb0TUkaMDFNk0c+ShTGm3Rgw6Xv03aFsk9DnTvf5fXzdpZpu+73NHga8LbBkYYxpN0SEs7aksz1DWL32vyEd882+b6iKg5FdcsIcXXSzZGGMaVcuu/UJvD7ljUVPhbT/2s1LAei+bGs4w4p6liyMMe1Kt2OGMGC78HbZMkKZomH53x8DYORPDr8XVVtiycIY0+6c+k0KhanK17vWNbnvzsHdyChWOg88oRUii16WLIwx7c74S38OwEdLXmty3y2VO8gqiSOuc+dwhxXVwposRGSiiKwTkY0iclsD2/uIyLsislJEFohIVtC2+0VktYisEZE/io3FYYxpIQNOmkj6fmXu5y80ue9Ozz56qI3wELZkISJeYAZwNjAYuFxEBtfZ7UHgOVUdCtwN/M499lvAKcBQ4HjgROD0cMVqjGlfvCkpHLfDy560xn+DHqg+QEmCj25JXVspsugVzpLFKGCjqn6lqpXAbOC8OvsMBgJT0r0ftF2BJCABSATigZ1hjNUY084M3tOBnZ1g586vDrnP9j3fANBxXfNHq20rwpksegHfBC3nueuCrQAucl9fAKSJSIaqLsJJHtvdv7dU9fCezzfGmEaccr4zJ8Xtr15/yH3u+vfPAOj37YtbJaZoFs4Zfhoq39Xtp3Yr8LiIXAUsBPKBahE5FhgEBNow5ovIGFVdWOsCItOAaQBZWVkUFBS0YPhHrrg49KdEo0UsxgyxGbfF3Hoairt31kjiv1Z6xnU95HfHuPTRLCveTFryUa3+/RJt9zqcySIPCB7LNwvYFryDqm4DLgQQkVTgIlUtdpPAJ6pa6m77N3AyTkIJPn4mMBMgNzdXMzMzw/RWmi8aY2pKLMYMsRm3xdx66sWdmcnRb8axJSn/kO9pb8UuAAYcl0vHCLzvaLrX4ayGWgz0F5G+IpIATAbmBu8gIpkiEohhOvC0+3orcLqIxIlIPE7jtlVDGWNaVF9/Bl979x5y+47SHaSVKWm9+rRiVNEpbMlCVauBG4G3cL7oX1bV1SJyt4gEBoUfC6wTkfVAN+Bed/0cYBPwBU67xgpV/We4YjXGtE+V5WXsSfZRtK/h/jM7KwvJbOcDCAaEsxoKVZ0HzKuz7jdBr+fgJIa6x/mA68IZmzHGXDD4EhaWPsOadR9zSu4F9bbv0VK6VCdGILLoY09wG2ParYHHjAJgyfOPNLh9n1SSRlJrhhS1LFkYY9qtnv1zSKpQtieXN7h9n6eSDgWlrRxVdLJkYYxpt7ypqfTYF8fO3qn1tvn8PsqSoEvWsRGILPpYsjDGtGvpB2Cz1n+GoqTCec6hY2LH1g4pKlmyMMa0a4NT+1OUolT5qmqt31PkPBbWMbl9jzYbYMnCGNOuZXfMxu+BLdu+rLV+b6GbLFIyIhFW1LFkYYxp17K7HwfAV1uW11q/d+8O58Wb79Y9pF2yZGGMadeOyXJmwNu8c32t9XtKnAf1jvnRLa0eUzSyZGGMadcyjx5IUoWyreSbWutLyp0G7vT0oyIRVtSxZGGMadfiOnak0wEPhQdq94gq+nABACkpnSIQVfSxZGGMafc6l8Luvfm11snw4wFISrWus2DJwhhjyEjqwt50b611e1d9Tny1Et+h/gN77ZElC2NMu1e+t5CihMraK/tnk1AFkmgDCYIlC2OMYWj34RxIFA5UHahZV15dQWK1INLQpJ/tjyULY0y7l7xyEwAFQe0W5VpJos++IgPsThhj2r3scecCsGv35pp15f4KEv32FRlgd8IY0+5lduwJwM7Cg89alGsVieo91CHtTliThYhMFJF1IrJRRG5rYHsfEXlXRFaKyAIRyQradrSIvC0ia0TkSxHJDmesxpj2K6OTkyy2vvJczboKqknSsE4mGlPClixExAvMAM4GBgOXi8jgOrs9CDynqkOBu4HfBW17DnhAVQcBo4Bd4YrVGNO+dcp0koWemluzrlyqSbRkUSOcJYtRwEZV/UpVK4HZwHl19hkMBEbpej+w3U0qcao6H0BVS1W1LIyxGmPasbSMHnj8SsmBPTXrKqSaZImPYFTRJZxpsxcQPNhKHnBSnX1WABcBjwEXAGkikgEMAPaKyD+AvsA7wG2q6gs+WESmAdMAsrKyKCioP4FJJBUXF0c6hMMWizFDbMZtMbeepuJWv5+UctjjL6r5Hqnw+PGqN2LfK9F2r8OZLBrqnKx1lm8FHheRq4CFQD5Q7cZ1GjAc2Aq8BFwFPFXrZKozgZkAubm5mpmZ2XLRt5BojKkpsRgzxGbcFnPraSrulAqhpHBLzX4+UZLiEiP6fqPpXoezGioP6B20nAVsC95BVbep6oWqOhy43V1X7B77uVuFVQ28DowIY6zGmHYuVRM40DWtZrnao8R5rM0iIJzJYjHQX0T6ikgCMBmYG7yDiGSKSCCG6cDTQcd2FpGu7vI4oPY0VsYY04LSNJFSPfgEd7VHiRdLFgFNJgsRqTeYu4gMbOo4t0RwI/AWsAZ4WVVXi8jdIjLJ3W0ssE5E1gPdgHvdY304VVTvisgXOFVaT4T0jowxphlSJIlSz8F5uH0eiBN7ziIglLT5XxH5taq+DCAiPwOuxenJ1ChVnQfMq7PuN0Gv5wBzDnHsfGBoCPEZY8wRSyOJ/fEH+9BUeSDeY72hAkKphhoLXCkir4jIQpyeSqPCGpUxxrSyuPxdHPD6AfD5fagHvNZmUaPJZKGq24H/AKOBbJyH6ErDHJcxxrSqgk4eKuIVVaXK5wxXvn7P+iaOaj9CabOYj/N8xPHAOcAjIvJguAMzxpjWNCRjCCpCeXU5VZXlAOR0HRbhqKJHKNVQM1R1iqruVdVVwLeA6HpaxBhjjtA6txRxoLyEynKnV1Sc19osAkKphnq9znK1qv42fCEZY0zrG91rNAD7S/dS7U6CZA3cBzXZeiMi+zj45HUCEA+UqqrNYm6MaTMW71gC3aGsrBjiOwAQF5cQ4aiiRyglizRVTXf/knDGcpoR/tCMMab1TOx3DgAHyoqpqqoArGQR7LCf4HarpcaFIRZjjImY5KRUAPYfKKGq0tos6gqlGurCoEUPkEv9AQGNMSamdUhMB6CsbB9VyZ0BSIhLjGRIUSWUJ07ODXpdDWym/rwUxhgT01KSnUEED1Tso7KmZGFtFgFNJgtVvbo1AjHGmEhKSnZKFm+sf41pHW8CIN5KFjUOmSxE5E80Ut2kqv8TloiMMSYCkjo4bRan9ziVqirnobz4eCtZBDRWsljSalEYY0yEJbrVUFVV5fh9zuizXq+NDRXQ2J14wR1m3Bhj2ryEFCdZVFaV4/c5Awp6LFnUaKzr7GeBF26VlDHGtFmJyU41VMmyxfjd38kej81nEdBYsgieQ/uUcAdijDGRlJCUAoBnyEB8frdkYcmiRmPJ4oifpRCRiSKyTkQ2ishtDWzvIyLvishKEVkgIll1tqeLSL6IPH6ksRhjTGM8cXF4fUqVrxpnsk7wSDhnno4tjVXIHSciK3FKGP3c17jLqqqNzmInIl6cYUHGA3nAYhGZq6rBc2k/iDM/xrMiMg74HXBl0PbfAh8c1jsyxphmivMJ1VqF3+8mC2uzqNHYnRh0hOceBWxU1a8ARGQ2zsN8wcliMHCL+/p9oGaEWxEZiTMv939wnho3xpiwivdDpb8Kf6AayubgrnHIZKGqW47w3L2Ab4KW83AmUQq2AmdgwseAC4A0EckA9gAP4ZQyzjzCOIwxJiRxPqFKq/FroM3CqqECwlnGkgbW1W0HuRV4XESuAhYC+ThDivwImKeq34g0dBr3AiLTgGkAWVlZFBQUtEDYLae4OPbmiIrFmCE247aYW0+occf5hYqqCvbvd2aO3r+/LGLfK9F2r8OZLPKA3kHLWcC24B1UdRtwIYCIpAIXqWqxiIwGThORHwGpQIKIlKrqbXWOnwnMBMjNzdXMzMywvZnmisaYmhKLMUNsxm0xt55Q4o5XwedREpMSoRw6duwc0fcbTfc6nMliMdBfRPrilBgmA98L3kFEMoEiVfUD04GnAVT1iqB9rgJy6yYKY4xpafF+oQoffusNVU9jY0N9QeNjQzXaG0pVq0XkRuAtwAs8raqrReRuYImqzgXGAr8TEcWphvrx4b8FY4xpGXF4qFYf6rcnuOtq7E581/1v4Av8efe/VwBloZxcVecB8+qs+03Q6znAnCbOMQuYFcr1jDHmSMT7PVSJL6iB23pDBTTZG0pETlHV4Ce4bxORj4C7wx2cMca0pji8dZKFVUMFhHInUkTk1MCCiHwLSAlfSMYYExnx6qEa/8GH8qxkUSOUCrlrgadFpCNOG0YxcE1YozLGmAiIx0upVB0sWVibRY1QZspbCgwTkXRAVDW6Ov8aY0wLicdLFf6a3lBe6w1Vo8k7ISLdROQp4CX3GYjBInJtK8RmjDGtKk68VHn8VrJoQChpcxZO99ee7vJ64CfhCsgYYyIlXuKoFkX9zlMD4rFkERBKsshU1ZcBPzjPTwC+sEZljDEREI+Xao/id7/ivNYbqkYod2K/O7ifAojIyTiN3MYY06bESxxVHq2Z/Ei81hsqIJQy1s+AuThzWnwEdAUuCWtUxhgTAXGeOHyiqNtm4fXERzii6BFSbygROR0YiDOS7DpVrQp7ZMYY08riPF78gj2U14BQekNtAn6gqqtVdZWqVonIm60QmzHGtCoPXnwerDdUA0JJm1XAGSLyjIgkuOt6hTEmY4yJCK9bstCakoUli4BQkkWZql4GrAH+KyJ9aGQ0WmOMiVUe8eDzgq+mZGHVUAGhpE0BUNX7RWQpzjMXXcIalTHGRECcO+e2LzCfhZUsaoRyJ4KHFH9XRL4NTA1fSMYYExled+DAKq0GrM0iWGOTHx2nqmuBfBEZUWezNXAbY9ocr1uyqPJVQfzB5GEaL1n8DPgh8FAD2xQYF5aIjDEmQrxutVO1W7IQK1nUaGzyox+6/z2juScXkYnAYzjTqj6pqr+vs70PzrzbXYEi4PuqmiciOcBfgHScoUXuVdWXmhuHMcaEwhNUshBVxJ6zqNFYNdSFjR2oqv9obLuIeIEZwHggD1gsInNV9cug3R4EnlPVZ0VkHPA74EqcaVunqOoGEekJLBWRt1R1b0jvyhhjmiHOrXaq1mpEQUQiHFH0aKyMdW4j2xRoNFkAo4CNqvoVgIjMBs4DgpPFYOAW9/X7wOsAqrq+5kKq20RkF07pw5KFMSZsvBIH6jRwiz0gUEtj1VBXH+G5ewHfBC3nASfV2WcFcBFOVdUFQJqIZKhqYWAHERkFJACb6l5ARKYB0wCysrIoKCg4wpBbVnFx7I23GIsxQ2zGbTG3nlDj3v/xIhgN5b5KPH6J6HdKtN3rkFood3adAAAf7ElEQVRvROQ7wBAgKbBOVe9u6rAG1tXN1bcCj4vIVcBCIB+oDrpuD+B5YKoGHqkMPpnqTGAmQG5urmZmZjb5XlpbNMbUlFiMGWIzbou59YQSd+fTx0LlXFScJ5Yj/V4jff1gTSYLEfkr0AE4A3gSuBj4LIRz5wG9g5azgG3BO6jqNuBC9zqpwEWBaVvdaVz/BfxKVT8J4XrGGHNEvOL2hsKqoeoKpan/W6o6BdijqncBo6mdBA5lMdBfRPq6Y0pNxhnqvIaIZIrUTHI7HadnFO7+r+E0fr8S2lsxxpgj43Xnr6jCj8eSRS2hJIsD7n/L3J5JVUDfpg5yZ9S7EWd4kDXAy6q6WkTuFpFJ7m5jgXUish7oBtzrrr8UGANcJSLL3b+cUN+UMcY0x8GShc9KFnWE0mbxpoh0Ah4AluG0OzwZyslVdR4wr8664OFD5gBzGjju78DfQ7mGMca0FK83kCz8DTa6tmehTH70W/flq+48FkmBdgVjjGlL3t7yFnSGKnxWDVVHKA3cXuA7QHZgfxFBVR8Ob2jGGNO6zj32PBYWPkG1+PGolS2ChVIN9U+gHPgCqNd91Rhj2oqaaijxW5tFHaEkiyxVHRr2SIwxJsK8nngAqlFrs6gjlN5Q/xaRCWGPxBhjIiwuULLw+BGrhqollJLFJ8Br7vMQVThPZquqpoc1MmOMaWU1Q5SL4rFK91pCSRYP4TyI94WqWi2eMabNiotzq6E8SpLfShbBQqmG2gCsskRhjGnrAnNuV3vUqqHqCKVksR1YICL/BioCK63rrDGmrfG6JQufp+GRUNuzUJLF1+5fgvtnjDFtUqDNwifYcxZ1NJos3AfyUlX1560UjzHGRIzHHUjQShb1Ndpmoao+YEQrxWKMMREl7rSqfg94LF3UEko11HIRmQu8AuwPrGxqDm5jjIk1Hjn4+9me4K4tlGTRBSgExgWtC2UObmOMiSlet2QBVrKoK5RRZ490Lm5jjIkJEpQsxJJFLU0+ZyEiWSLymojsEpGdIvKqiGS1RnDGGNOaPJ6gaqgIxhGNQnko7xmc6VB7Ar1wRqF9JpSTi8hEEVknIhtF5LYGtvcRkXdFZKWILAhOQiIyVUQ2uH9TQ3s7xhjTfIGH8sC6ztYVSrLoqqrPqGq1+zcL6NrUQW632xnA2cBg4HIRGVxntwdx5tkeCtwN/M49tgtwB3ASMAq4Q0Q6h/iejDGmWby1qqFMsFCSRYGIfF9EvO7f93EavJsyCtioql+paiUwGzivzj6DgXfd1+8Hbf82MF9Vi1R1DzAfmBjCNY0xptk83qCShaWLWkJJFtcAlwI7cIb+uNhd15RewDdBy3nuumArgIvc1xcAaSKSEeKxxhjTojzBJQurhqollN5QW4FJzTh3Q3e6bs/lW4HHReQqYCGQD1SHeCwiMg2YBpCVlUVBQUEzwgyf4uLYm6o8FmOG2IzbYm49ocZdXLKv5rWoRvQ7Jdru9SGThYj8ppHjVFV/28S584DeQctZwLY6J9kGXOheLxW4SFWLRSQPGFvn2AUNBDETmAmQm5urmZmZTYTU+qIxpqbEYswQm3FbzK0nlLirDxxMDh7xRPy9Rvr6wRqrhtrfwB/AtcD/hnDuxUB/EekrIgnAZJxeVTVEJNOdVAlgOvC0+/otYIKIdHYbtie464wxJmxq9YayNotaDlmyUNWHAq9FJA24Gbgap6H6oUMdF3R8tYjciPMl7wWeVtXVInI3sERV5+KUHn4nIopTDfVj99giEfktTsIBuFtVi5rx/owxJmSBgQTBHsqrq6lRZ7sAPwWuAJ4FRri9k0KiqvOAeXXW/Sbo9RxgziGOfZqDJQ1jjAk7j9hwH4fSWJvFAzjtCTOBE1S1tNWiMsaYCLCSxaE11mbxM5yntn8FbBOREvdvn4iUtE54xhjTejw2kOAhNdZmEcozGMYY02Z4PPE1ry1V1GYJwRhjXF7vwa9EK1nUZsnCGGNcHm9wycK+HoPZ3TDGGFfwEOVWsqjNkoUxxriCBxK0VFGbJQtjjHFJrd5Q9vUYzO6GMca4xONB1BmzVMTKFsEsWRhjjEtEEHd8a2uzqM2ShTHGBAkkC3uCuzZLFsYYEyQw97aVLGqzZGGMMUE8gWooa7OoxZKFMcYEOVgNZV+PwexuGGNMEGvgbpglC2OMCRKohrKus7VZsjDGmCCBFGEli9rCmixEZKKIrBORjSJyWwPbjxaR90XkcxFZKSLnuOvjReRZEflCRNaIyPRwxmmMMQEea7NoUNjuhoh4gRnA2cBg4HIRGVxnt18BL6vqcGAy8Gd3/SVAoqqeAIwErhOR7HDFaowxAeJ2nd1Wmh/hSKJLOFPnKGCjqn6lqpXAbOC8OvsokO6+7ghsC1qfIiJxQDJQCdjsfMaYsAtUPnXfG9Ewok44k0Uv4Jug5Tx3XbA7ge+LSB4wD7jJXT8H2A9sB7YCD6pqURhjNcYY4GA1VOLRfSIbSJQ55LSqLaCh1iGts3w5MEtVHxKR0cDzInI8TqnEhzMHeGfgvyLyjqp+VesCItOAaQBZWVkUFBS09Hs4IsXFxZEO4bDFYswQm3FbzK3ncOIOdJ31+3wR/U6JtnsdzmSRB/QOWs7iYDVTwLXARABVXSQiSUAm8D3gP6paBewSkY+AXKBWslDVmcBMgNzcXM3MzAzH+zgi0RhTU2IxZojNuC3m1hNq3M6YUEpCXELE32ukrx8snNVQi4H+ItJXRBJwGrDn1tlnK3AmgIgMApKA3e76ceJIAU4G1oYxVmOMAcBrAwk2KGzJQlWrgRuBt4A1OL2eVovI3SIyyd3tZ8APRWQF8CJwlaoqTi+qVGAVTtJ5RlVXhitWY4wJCPSGsofyagtnNRSqOg+n4Tp43W+CXn8JnNLAcaU43WeNMaZV1TyUJ/acRTC7G8YYE+TgEOX29RjM7oYxxgQJlCysGqo2SxbGGBMk0LBt1VC12d0wxpgggfLE+qJ1EY0j2liyMMaYBgzKqDuUXftmycIYYxpg1VC1hbXrbKRVVVWRl5dHeXl5RK7v8/nYvXt3RK7dXLEYM0Q+7qSkJLKysoiPj49YDKZlWbKorU0ni7y8PNLS0sjOzo5Iz4aqqqqY+/KIxZghsnGrKoWFheTl5dG3b9+IxGBanliyqKVN343y8nIyMjKsC5wJKxEhIyMjYiVY07IOPpRn3xvB2nSygMPvK73lyilsuXJKmKIxbZX9IGl77KG82uxuhNm9997LkCFDGDp0KDk5OXz66ac123bv3k18fDx/+9vfah2TnZ3NCSecwLBhw5gwYQI7duyod96xY8cycOBAcnJyyMnJYc6cOQCkpqaG9w0Z005YNVRtdjfC6JNPPuHNN99k2bJlrFy5knfeeYfevQ+O2v7KK69w8skn8+KLL9Y79v3332fFihXk5uZy3333NXj+F154geXLl7N8+XIuvvjisL0PY9ojj8e+HoPZ3Qij7du3k5mZSWJiIuCMTd+zZ8+a7S+++CIPPfQQeXl55Oc3PN/vmDFj2Lhx42Ffe8GCBXz3u9+tWb7xxhuZNWsWxcXFDBw4kHXrnAeOLr/8cp544onDPr8xbZ1VQ9XWpntDBdtx331UrGl6Sozytc4+obRbJA46ju6//OUht48fP5777ruPAQMGcNZZZ3HZZZdx+umnA/DNN9+wY8cORo0axaWXXspLL73ET3/603rnePPNNznhhBMaPP8VV1xBcnIyAO+++y4ZGRlNxtyxY0cef/xxrrrqKm6++Wb27NnDD3/4wyaPM6a9sXao2ix1hlFqaipLly5l5syZdO3alcsuu4xZs2YBMHv2bC699FIAJk+eXK8q6owzziAnJ4eSkhKmT5/e4PmDq6FCSRQB48eP54QTTuDHP/4xTz75ZPPenDFtlN/9b8Unn0U0jmjTbkoWjZUAggVKFH2ef65Fruv1ehk7dixjx47lhBNO4Nlnn+Wqq67ixRdfZOfOnbzwwgsAbNu2jQ0bNtC/f3/AabM4kikV4+Li8Pv9NcvB3Tr9fj9r1qwhOTmZoqIisrKymn0dY9oadSfhTjvl1AhHEl2sZBFG69atY8OGDTXLy5cvp0+fPqxbt479+/eTn5/P5s2b2bx5M9OnT2f27Nktdu0+ffrw5ZdfUlFRQXFxMe+++27NtkceeYRBgwbx4osvcs0111BVVdVi1zUm1rmzqhLnaTe/pUMS1mQhIhNFZJ2IbBSR2xrYfrSIvC8in4vIShE5J2jbUBFZJCKrReQLEUkKZ6zhsH//fqZOncrgwYMZOnQoX375JXfeeScvvvgiF1xwQa19L7roogZ7RTVX7969ufTSSxk6dChXXHEFw4cPB2D9+vU8+eSTPPTQQ5x22mmMGTOGe+65p8Wua0ysCyQLrzf2RjIIp7ClThHx4sylPR7IAxaLyFx3KtWAX+HMzf0XERmMMwVrtojEAX8HrlTVFSKSAcTcz98RI0bw8ccf11t/55131lsXSCYAmzdvbvLcCxYsaHB9aWlpzev777+f+++/v94+a9asqXn98MMPN3ktY9oTv1sNZSWL2sJ5N0YBG1X1KwARmQ2cBwQnCwXS3dcdgW3u6wnASlVdAaCqhWGMs5aWaqswxsSmmpKFxxvROKJNOJNFL+CboOU84KQ6+9wJvC0iNwEpwFnu+gGAishbQFdgtqrW+4ksItOAaQBZWVkUFBTU2u7z+SJaH+/z+SJ27eaKxZghOuL2+Xz1PoONKS4uDmM04RGLMcPhxa1uuqisqDqsf8+WFm33OpzJoqFOylpn+XJglqo+JCKjgedF5Hg3rlOBE4Ey4F0RWaqq79Y6mepMYCZAbm6u1u09FBhOI5Iiff3miMWYIfJxe73ew+7BdiQ93iIlFmOG0ONW95srLTU94u810tcPFs4G7jygd9ByFgermQKuBV4GUNVFQBKQ6R77gaoWqGoZTlvGiDDGaowxQHADt7VZBAtnslgM9BeRviKSAEwG5tbZZytwJoCIDMJJFruBt4ChItLBbew+ndptHcYYExZ+t2Th9SZENpAoE7bUqarVInIjzhe/F3haVVeLyN3AElWdC/wMeEJEbsFJ6FepqgJ7RORhnISjwDxV/Ve4Yg129X+uBuCZic+0xuWMMVFGa3pDWQN3sLA+Z6Gq81R1gKr2U9V73XW/cRMFqvqlqp6iqsNUNUdV3w469u+qOkRVj1fVX4QzznDasWMHkydPpl+/fgwePJhzzjmH9evXs3nzZpKTk8nJyWHYsGF861vfqnlYLyMjo17j1vnnn8/LL79c7/yfffYZY8aMYeDAgRx33HH84Ac/oKysrLXeHq+//npNl9/Dcaih1Js7pPvw4cNtSHfTIuw5i4bZE9xhpKpccMEFjB07lk2bNvHll19y3333sXPnTgD69evH8uXLWbFiBVOnTuW+++4jJSWFCRMm8Prrr9ecp7i4mA8//LDWKLIAO3fu5JJLLuEPf/gD69atY82aNUycOJF9+/aFFF91dXWjy6FobrJoyKJFi5o9pPv8+fNtSHfTItyCBXFxliyCWbIIowULFhAfH8/1119fsy4nJ4fTTjut3r4lJSV07twZcIYNDx7647XXXmPixIl06NCh1jEzZsxg6tSpjB49GnBGybz44ovp1q0bRUVFnH/++QwdOpSTTz6ZlStXAs4DgdOmTWPChAlMmTKFWbNmcckll3DuuecyYcIEAB544AFOPPFEhg4dyh133FFzveeee46hQ4cybNgwrrzySj7++GPmzp3Lz3/+c3Jycti0aRObNm1i4sSJjBw5ktNOO4217ii+X3/9NaNHj+bEE0/k17/+dYP3y4Z0N9HA4yYLK1nU1m6a+//w2R9YW9T0EOWBfQJtF405rstx/O+o/z3k9tWrVzNy5MhDbt+0aRM5OTns27ePsrKymiqXiRMn8oMf/IDCwkIyMjKYPXs2N910U73jV61axdSpUxs89x133MHw4cN5/fXXee+995gyZQrLly8HYOnSpXz44YckJycza9YsFi1axMqVK+nSpQvz5s1jw4YNfPbZZ6gqkyZNYuHChWRkZHDvvffy0UcfkZmZSVFREV26dGHSpEl897vfrfmlfuaZZ/LXv/6V/v378+mnn/KjH/2I9957j5tvvpkbbriBKVOmMGPGjAZjnjBhAnfffbcN6W4iyqMCKHFx1sAdzEoWERSohtq0aROPPvoo06ZNAyAhIYFJkyYxZ84cCgoKWL58ec2v/lB9+OGHXHnllQCMGzeOwsLCmnaQSZMm1XxpgjNkeZcuXQB45513ePvttxk+fDgjRoxg7dq1bNiwgffee4+LL764pt93YP9gpaWlfPzxx1xyySXk5ORw3XXXsX37dgA++ugjLr/8coCauOo6kiHdx48fb0O6mxYR+FL02nAftbSbu9FYCSBYS/aGGjx4MK+99lpI+06aNImrrz5Ymrn88su55557UFXOO++8Bh84GzJkCEuXLuW8886rt83pVFZbYDKXlJSUWuuDl1WV6dOnc91119Xa549//GOTk8H4/X46depUU4I51PUb09wh3efPn0+PHj2aPP+h2JDuJsDjPpXntTaLWqxkEUZnnHEGFRUVteq4Fy9ezAcffFBv3w8//JB+/frVOnbDhg3MmDGj5hd5XTfeeCPPPvtsrR5Df//739mxYwdjxoyp+WJdsGABmZmZpKenN3ieYOPHj+fpp5+uGZAwPz+fXbt2ceaZZ/Lyyy9TWOgM01VUVARAWlpaTYN6eno6ffv25ZVXXgGcxLNixQoATjnllJp2mEBcddmQ7iYaBJKFx2tdZ4NZsggjEeG1115j/vz59OvXjyFDhnDnnXfWNNoG2iyGDRvGL3/5y1pVHB6Ph4suuojCwkLGjBnT4Pm7devG7NmzufXWWxk4cCCDBg3iv//9L+np6dx5550sWbKEoUOHctttt/Hss8+GFPP48eP53ve+x+jRoznhhBO4+OKL2bdvH0OGDOH222/n9NNPZ9iwYTXtBZMnT+aBBx5g+PDhbNq0iRdeeIGnnnqKYcOGMWTIEN544w0AHnvsMWbMmMGJJ554yDFvSktLbUh3E3Hplc7XYqI3McKRRBdpqLoiFuXm5uqSJUtqrVuzZg2DBg06rPO0ZDVUVVVVxMcrOlyxGDNER9yH+3krKCiIqrF/QhGLMcPhxb342kv5uGwVNzzwNgkRrHJsrXvtjruX29R+7abNIlT25LYx7VvOPY9x7IIFEU0U0ciqoYwxJkh8jx50PkQ7YXtmycIYY0yT2nyyaCttMia62efMtHVtOlkkJSVRWFho/yObsFJVCgsLSUpKinQoxoRNm27gzsrKIi8vj927d0fk+j6fD2+M9dWOxZgh8nEnJSXZw3qmTWvTySI+Pp6+fftG7Pqx2M0wFmOG2I3bmFjRpquhjDHGtAxLFsYYY5pkycIYY0yT2sxwHyKyG9gS6TjqyAQKIh3EYYrFmCE247aYW08sxt1aMfdR1a5N7dRmkkU0EpEloYy5Ek1iMWaIzbgt5tYTi3FHW8xWDWWMMaZJliyMMcY0yZJFeM2MdADNEIsxQ2zGbTG3nliMO6pitjYLY4wxTbKShTHGmCZZsmgmERkoIsuD/kpE5CdB228VERWRTHdZROSPIrJRRFaKyIhoiVlE7hSR/KD15wQdM92NeZ2IfLu1Y24sbnfbTW5sq0Xk/miJu5F7/VLQus0isjxaYm4i7hwR+cRdt0RERrn7R/PnepiILBKRL0TknyKSHnRMxO+1G8ct7md3lYi8KCJJItJXRD4VkQ3u5yXB3TfRXd7obs9u1WBV1f6O8A/wAjtw+isD9AbewnnuI9Nddw7wb0CAk4FPoyVm4E7g1gb2GQysABKBvsAmwBtFcZ8BvAMkutuOisa4634+gtY/BPwmGmNu4F6/DZztrj8HWBD0Olo/14uB09311wC/jaZ7DfQCvgaS3eWXgavc/0521/0VuMF9/SPgr+7rycBLrRmvlSxaxpnAJlUNPBT4CPALILhB6DzgOXV8AnQSkR6tHGewujE35DxgtqpWqOrXwEZgVKtEd2jBcd8A/F5VKwBUdZe7T7TFXe9ei4gAlwIvuquiLWaoHbcCgV/mHYFt7uto/lwPBBa66+cDF7mvo+lexwHJIhIHdAC2A+OAOe72Z4Hz3dfnucu42890P0etwpJFy5iM+z+9iEwC8lV1RZ19egHfBC3nuesipSZm141uNcLTItLZXRdtMUPtuAcAp7lF8g9E5ER3fbTFXfdeA5wG7FTVDe5ytMUMteP+CfCAiHwDPAhMd9dHW9zBMa8CJrmvL8Ep8UOUxKyq+Tj3citOkigGlgJ7VbW6gdhq4na3FwMZrRWvJYsj5NYnTgJeEZEOwO3AbxratYF1EemKFhyzu+ovQD8gB+dD+1Bg1wYOj1j3uQbijgM641R//Bx42f2lFTVxNxBzwOXUTiBREzM0GPcNwC2q2hu4BXgqsGsDh0fLvb4G+LGILAXSgMrArg0c3uoxuz/KzsOpCusJpABnN7BrILaIxm3J4sidDSxT1Z04X7h9gRUishnIApaJSHecXwi9g47L4mBRvrUFx4yq7lRVn6r6gSc4WCSPppihTtw48f3DrQL5DPDjjKcTTXHXjRm3yuFC4KWg/aIpZqgf91TgH+7rV4jOz0jdz/VaVZ2gqiNxEvMmd79oifks4GtV3a2qVTj391s4VXmBuYaCY6uJ293eEShqrWAtWRy5ml+IqvqFqh6lqtmqmo3zjztCVXcAc4Epbu+Rk4FiVd0e6ZgB6tQxX4BTfAcn5sluL4y+QH/gs1aLsr66v8Zfx6nfRUQGAAk4A69FU9x1YwbnS2KtquYFrYummKF+3NuA093X44BA9Vk0f66Pcv/rAX6F01gM0XOvtwIni0gHt0R8JvAl8D5wsbvPVOAN9/Vcdxl3+3vqtna3itbuAdCW/nAapAqBjofYvpmDvaEEmIHz6+YLIDdaYgaed2NaifOB7BG07XY35nW4vWGiKO4E4O84yW0ZMC6a4j7U5wOYBVzfwP4Rj7mRe30qTn36CuBTYKS7Ppo/1zcD692/3+M+hBxl9/ouYK37GX4ep4fWMTjJayNOKS7Q2y/JXd7obj+mNWO1J7iNMcY0yaqhjDHGNMmShTHGmCZZsjDGGNMkSxbGGGOaZMnCGGNMkyxZmIgTkYygEUN3SO0RcD9u5Vh+KiJfukOfvCsifepsT3fjezxo3QJ39NJAzEfVOeZicUYgznWX40XkWXc01DUiMt1dnyQin4nICnck0ruCzvGUu36liMwRkdRDxH++iDQ0ggAiUtr8O9Pg+d4JGhrGtHHWddZEFRG5EyhV1QcjdP0zcEZOLRORG4CxqnpZ0PbHgK5Akare6K5bgDNq75IGzpcG/AvnmZAbVXWJiHwPmKSqk90hYr4ExuKMUpyiqqUiEg98CNysqp+ISLqqlrjnfBjYpaq/b+B6H7vnLmhgW6mqNphkmkNEpgJZqnpvS53TRC8rWZioFvg1LCJj3cECXxaR9SLyexG5wv0l/oWI9HP36yoir4rIYvfvlMO5nqq+r6pl7uInOMMtBGIZCXTDGa47VL8F7gfKgy8DpLhDNiTjjFlUoo7Ar/9490/duAKJQtxj6v3Kc59irwgkCnHmRVjk3offBu2X6paalrn37jx3/W9F5Oag/e4Vkf8RkR4istAtNa0SkdPcXebiPDVt2gFLFiaWDMN5KvcE4EpggKqOAp4EbnL3eQx4RFVPxBmS+skjuN61OHM1BIaMeAhnwMKGPON+mf7a/UJHRIYDvVX1zTr7zgH24wzauBV4UFWL3GO84kyItAuYr6qfBg4SkWdw5mo4DvhTAzGcgvMke8BjwF/ce7EjaH05cIGqjsCZE+QhN+ancIeTcN/vZOAF4HvAW6qag/NvsBxAVfcAiSLSaiOfmsixZGFiyWJV3a7O/BWbOPgL/wsg2319FvC4+4U7F0h3q4IOi4h8H8gFHnBX/QiYp6rfNLD7Fap6As6w46cBV7pfto8AP2tg/1GAD2ek0b7Az0TkGAB1BnTMwSnRjBKR4wMHqerV7jFrgMvqnRV6ALuDlk/h4FhJzwe/PeA+EVmJM3lUL6Cbqm4GCt0kNwH4XFULcSYRutqtIjxBVfcFnWuXG5Np4yxZmFhSEfTaH7TsxxmuHJzP9GhVzXH/etX5ckNEAqWAeQ1dRETOwhk7aJKbmABG48z5sRlnDoIpIvJ7qJmXAPc6/4eTDNKA44EF7jEnA3PdRu7vAf9R1Sp1Jmz6CCcx1VDVvcACYGKd9T6c0Wovor4DOOMH1Tqkgf2uwGl3Gekmpp1Bxz2JM1vb1cDT7jUXAmOAfOB5EZkSdK4k97qmjbNkYdqat4EbAwsiklN3B1W92k0k59Td5v6q/htOotgVdMwVqnq0OqMJ34ozO9xtIhInB+dZjwe+C6xS1WJVzdSDIxB/4p5zCU7V0zhxpOAkkrVue0sn91zJuKPTuvsd664X4FycwefqWgMcG7T8EU5VEjgJIqAjTgN5ldugH9zj6zWcBHUiztTAiNMjbJeqPoFTVTUiKJbuOANmmjYuruldjIkp/wPMcKtY4nCm1bz+MI5/AEjFmcwKYKuqTmpk/0TgLTdReHGqdZ5o4hozgGdwRhoV4BlVXSkiQ4FnRcSL80PuZVV9063SelZE0t39V+BMRlTXQtz2B3W6Od4M/J/baP1q0H4vAP8UkSU47Q81iUdVK0XkfZzZ2nzu6rHAz0WkCigFAiWLkcAnenBWN9OGWddZY9oQt2vvP1X1nWYe78FpJL9ED0752ti15qrqu825loktVg1lTNtyH87cDodNRAbjzJXwblOJwrXKEkX7YSULY4wxTbKShTHGmCZZsjDGGNMkSxbGGGOaZMnCGGNMkyxZGGOMaZIlC2OMMU36f5WcQOiN/utuAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = lc.plot(color='C3', label='SAP Flux', linestyle='-')\n", "corr_lc.plot(ax=ax, color='C2', label='CBV Corrected SAP Flux', linestyle='-');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now the CBVs have been corrected only the astrophysical variability is left, the variability from spacecraft heating has been removed." ] } ], "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.6.3" } }, "nbformat": 4, "nbformat_minor": 2 }