{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# How to combine lightcurves from different Kepler quarters" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The Kepler mission took 17 quarters of continuous observations on the same targets. For applications such as finding long period planets or asteroseismology it is useful to be able to turn these quarters into one long light curve. We can do this using lightkurve's `append` function. The append function is called by doing\n", "\n", " stitched_lc = lc1.append(lc2)\n", " \n", "It works in the same way as appending a list in python.\n", "\n", "Below is an example of how to stitch together all the quarters for the exoplanet target Kepler-8b. (*Note: We'll use the jupyter notebook magic command `%%capture` to supress the output for the cell so that it doesn't print a long download message*)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "from lightkurve import KeplerLightCurveFile\n", "target = 'Kepler-8b'" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%%capture\n", "lc = KeplerLightCurveFile.from_archive(target, quarter=1).PDCSAP_FLUX\n", "for q in range(2,17):\n", " lc = lc.append(KeplerLightCurveFile.from_archive(target, quarter=q).PDCSAP_FLUX)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Above we have created a `KeplerLightCurve` object from the first quarter of the data, choosing to use the PDCSAP_FLUX. Then we have looped through all the quarters, appending each new light curve to the fist one. We can plot this new `KeplerLightCurve` using the `plot` function (and all the other functions from previous tutorials will work too!)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3XucJHV57/HP07AqYwuoiwZZjyyJQ0LEjO3qMHqMSEej4iXxFhUvoBtyOpKYRJIMx2xyko1xo5hEg2cMrqB4vOH1IF6IjhpOkmV0GUbACy2KlxUUJ8RNhiG6MM/5o6p6a2qqu6tm+lLV832/XvPa7qrq7mdrZ+vp36Wen7k7IiIiWVWGHYCIiJSLEoeIiOSixCEiIrkocYiISC5KHCIikosSh4iI5KLEISIiuShxiIhILkocIiKSy5HDDqBXtm7d6ieeeOLAP/fuu+/myCPLcRoVa++VJU5QrP1QljihfazXXnvtorsfl+e9yvE3zuDEE09k//79A//cxcVFtm7dOvDPXQ/F2ntliRMUaz+UJU5oH6uZfSfve6mrSkREclHiEBGRXJQ4REQkFyUOERHJRYlDRERyUeIQEZFclDhERCSXviUOM7vEzG43sxvb7Dcze4uZ3Wxm15tZLbH/aDP7vpld1K8YRUQkv37eAPhO4CLgsjb7nwY8PPyZBGbCPyO7gX/qY3yF0mg0AJiZmUndHknuFxEZtL4lDne/2sxO7HDIs4HL3N2Ba8zsWDM73t1vM7NHAw8GPg3s6FeMw1Cv11laWmq7f3Jysu2+tP21Wk3JREQGapglR04Avhd7fgA4wcx+CLwJeClQ7/QGZnYucC7Atm3bWFxc7FOo7R08eDDX8Z2SxnrMz89n/nvnjXWYyhJrWeIExdoPZYkTehvrMBOHpWxz4LeBT7r798zSDokd7H4xcDHAjh07fFg1Y4ZdqybP5w871jzKEmtZ4gTF2g9liRN6F+swZ1UdAB4ae74NuBWYAs4zs28DFwIvM7M9gw+vP+bm5qjVatRqNebm5pibm6NarVKtVlv7qtXqqv3RT3Jf9D4iIoM0zBbHFQQJ4v0Eg+IH3f024KzoADM7G9jh7tPDCXEwZmdnW4/TBsebzSbj4+Nr9mlsQ0SGoW+Jw8zeB5wObDWzA8CfAVsA3P1twCeBpwM3A8vAOf2KpWiiC35yJlXazKrocb1ep16vMz4+nvpeIiKD0s9ZVS/qst+BV3U55p0E03o3hXgSqNeDeQFRayTeKhERGaaRWcipjNrdszEzM9M1UbS770NEpN9UcmTIGo1G2y6rtGOixzMzM0oaIjIUShxD1m28I9oWbV9PwognHhGRjVJXVQHkmS0VTzDdypQkj4luFKzX6ywvL7OysgKwakqvWjEi0o0SR0GlTcNNJoqZmRkmJyfblimJb+9UymR+fn7VcXNzcxuOX0RGlxLHkLVrHaR98+80bVdEZFCUOIYs7Z6NPK+JVKtVgNZ9HgsLC4yNjTE+Ps7CwkKrW6qd5OtFRNpR4hiStHLpWcYvkvJ0Ky0uLq6qVaOWi4ish2ZVDUnUHRX9pI1fAKum4dbr9bazo5JTdjsdG48h/hkiIllYcAN3+e3YscP3798/8M9NfovPIs83/bytgk5lTHbu3MmWLVvW3boZpPWc12EoS5ygWPuhLHFC+1jN7Fp3z7Xukbqqhqzb7Km8F/NOU3v37NnT+sUpUpIQkXJR4hiCeOFCOFyHKv48foc4ZCuImFcRWxoiUnxKHEMQTxD1ep2pqanWDKi0rqNkt1LUQon2dboBsBMlDBFZDw2OD8Hy8jJLS0tMTU0xPj7O2NgYy8vLLCwsMDk5SbPZbCWJer1Os9ls7Zufn2d5eXnNe3aqedWp5EhaHSwRkU6UOIZgbGyMarXKvn37aDabAExMTLCyskKlUmF8fLyVJCITExPUajUqlQpjY2MArddCeuskrb5VMjkku83UChGRbjSraoP6Masi7R6P5L71zMgahRkgRVOWOEGx9kNZ4gTNqhp5nZJC2r71DJxrYFxE1kuJYwS1SwbT09Ot+ziUMERkvZQ4RkDWlkX8Po5+TO/tt04VfvOqVCrs27evZ+8nspkocYyIdhf+tKm6yRsO044bdSsrKyohL7JOmlW1ycTXM0/Ws9J0XBHJQi2OksvT5ZQ2xhHdK5JsgRRR1Dro9HdM3iAZTWmuVqutO/LTWlwikp0Sx4hI3l0eJYNIcowj/hpYfU9I0WWZWRZtq9Vqa5Jq1OISkfVR4hgB8Yt+9I06Ei3oBPCCF7yAu+66a9Vrx8bGVpU+mZiYKPQ38Xb3uHRaMTH+uiL/3UTKQomj5KIWxsLCAhAkkeXlZcbGxpiZmWFqaoqlpSUajUYraezbt2/VhbRer7fKmJSl5ZF3NpgShkjvKHGUXHThnJiYYGFhgaWlJSqVYM5DvV5fVTxxamqq9booQUSlRsoyNTVPteBeVhBOfr7IZqaSIxs0CiUHiqifpVzajZHEa4NBcK9Hcq325LYy3Q+y2f/9+6EscUJvS45oOq5sGp2Wyk0mDWBN0kjbFt0PIrKZqKtKSinr0rft1nIXkfVTi0NKJy0ZtEsIaa2MtBsda7Ua1WqVWq3G3Nxc63mlUqFarTI3N9faHh0zNzfHJz7xCd19LpuOxjg2aBT6OIuoF7H2agpuconf6H2bzSbuzuc+97kNvf+gbLZ//0EoS5ygsuoiwNrEkLzA500Y7WZQJRNG/H0XFxfzhi1Sen1LHGZ2CfAM4HZ3f0TKfgPeDDwdWAbOdvd5M5sAZoCjgXuA17n7B/oVp5RHt/GKTneEZ2l95L0XpNFocOjQIfbu3dvxdSKjpp8tjncCFwGXtdn/NODh4c8kQbKYJEgiL3P3b5jZQ4Brzewqd/9xH2OVEmh3j0aWgfJeDorHP0MtDtmM+pY43P1qMzuxwyHPBi7zYJDlGjM71syOd/fWrcvufquZ3Q4cByhxCNC55dEtQWSZgdXtvTQzSza7YY5xnAB8L/b8QLjttmiDmT0WuBfwzcGGJkXW6e7wZOHG6K75LFVx17OOu8hmNMzEYSnbWlO8zOx44N3Ay9197Z1YwTHnAucCbNu2bSjdBgcPHhz4Z67XKMQ6PT296vmePXsA2Llz55pj3Z1Dhw6xuLjIoUOH2L59O4cOHWJqaoqjjjqKyy+/PHM80edGnxf9rp1//vkceeSRre1FNgr//kVTljiht7EOM3EcAB4ae74NuBXAzI4GPgH8ibtf0+4N3P1i4GIIpuMOa1pcWabjQbljjWZNjY+Pt4o67tq1q/U4qte1srJCtVpdNU12y5YtwOpxkLQS8+1qX0WvT8Z04YUXlvqcFllZYi1LnNC7WIeZOK4AzjOz9xMMih9099vM7F7ARwnGPz44xPikYGZnZ5mammolipWVlVVl4yNRkcfJyUkqlcqqhFKv11laWlrz3t0G00epaypeej/PglbJSQhaEGvz6ud03PcBpwNbzewA8GfAFgB3fxvwSYKpuDcTzKQ6J3zpC4BfBh5oZmeH285294V+xSrl0Gg0mJiYAIJ1RqJV/aKy8FESGR8fp9lstgoSzs/PU6lUqNVqrfdqd7HLuqJitH337t29/Uv2ydTUVGrtrfn5+dy1tuLHp9X4ktGnO8c3aBTuHC2iosXaLoEULc52+lmIsR8lV8pyXssSJ+jOcZGe6zZbKt4CiVo2ZVqCtlqtpnbRiayHEocI6V1RUR9+/Jiy9ufPzs62/caZHK/IWm04rTy9bA5KHCIJ0QUxWaNqVKVNBGg0GtTr9Y6Jcz1JtN19NlIuShwiCe0uZGlVcsteq2ojJVqSqyamrZjYTXxwXuXpy0OJQ4S1U1Qj0bfs+D0k8cdluPEvrlMF4Xbl49vNMEvOqMqbNKS8NKtqg0ZhVkURDTvW+EU0fu9H9K26Wq0CsLy8zMTERCm6W6JzmqVsSpZxjl5OxU22Nob9759VWeIEzaoS6bl4iwNgaWmJSqXSukckugckOmZ5eZmxsTE+/vGPl+bCEYnGMBqNRtvB8Czl57PW7UqbaCDlphbHBo3CN44iKkusZYkzz3hMnjXc+6Us57UscYJaHCKS08zMDDt37my1MrodG+m2eFYeWdZNkXJQ4hDZJPbs2bOmsGO30vP9mnKrhFFuShwim1Sni3m3FkGn/XnXO1Hro3y6Jg4ze5C7357YdrK739S/sESk3/Jc/Lt1WeUZKE8ep4RRPllaHP/PzHa5++UAZvYa4JXAKX2NTER6anp6urWuCGS7+He6lyP+OGsiWW/CkWLJkjhOBy42s+cDDwa+Bjy2n0GJSO8lxziS8nQ7dZNMPhrrGC1dE0e4uNKngQuAFeACd1eZTZES2ch03HYti6xjFsOc1iv9kWWM4zPAbcAjCJZ3vcTMrnb38/sdnIj0xszMTGud9G6LU7W7wS9eaiX5+k7dVu3ikfLK0lX1Vnf/WPj4x2b2OILWh4iUTHyqbPLO8UiyZlUk+TxvspDRkaWr6mOJ53cD5VgvU0SA4CJ/0003cfLJJ7fWbB8bG1tVOr3ZbLZKqSwtLTE1NcXExMSa4o/x49IKInYaZI8o0ZRb15IjZvafQHTQvQjWDV9y92P6HFsuKjnSnWLtvbLECXDGGWdw5513rtpWqVRWrdMeFW2MChhWq9VVySR6TSTLSojrGc8oy3ktS5ww4JIj7n6/xIf8GppVJVIqjUaDk046KdfaIdFsqHj133YrALZrZWixptG0riKHZnaNu5/Wh3jWTS2O7hRr75UlTuhfrFlbFO3GTtKU5bxuJM54uf5BLGI10BaHmT0n9rQC7OBw15WIbFKNRqM1XjIxMdHaHk8Q8aSSljDKNC13amqq42JVc3Nzq5LBKMsyq+qZscd3A98Gnt2XaESkNJLdUtEsrWi8JOqqWlpaYnJyctV4Sryrq9lsUq/XM7VEhiVa3najx6SJtza6TZWOxM9flLz37du3rs9fjyxjHOcMIhARKZcoMSwvL7OystIaNF9YWGBlZaW1L0oYS0tLrW/j0QJZkXiLZbOKJ55uSWi9SapX2iYOM/t7OnRJufvv9iUiESmNtIHvffv2ZZp+W6Zuqrm5udSJAbt37161HG9SclJBvGsvvvTusBNBXp1aHIMfaRaR0siyvOx6X19EaeVWFhcXu66YmFzzJDp2vcmiVqutWe99EIPrcZ0Sx3vCm/1ERIT0KcZ5EmiyBRKpVqsAq9ZlL/JU5k6J44tADYJuK3f/ncGEJCJSTMkJAZ0KR7ar9TU7O7vmPpi0O/CzvOewuvs6JQ6LPX58vwMRESm6ZIujU1dV8nl81th6730ZdsKIdEoculdDRCQm72qJWZNKu+2dCksOU6XDvp83s+vN7IbY4+vN7AYzu35QAYqIFEF0r0oe7cY01qsolYg7tTh+YWBRiIgUXNYLdqfWRZbWSJ5jCzfG4e7fGWQgIiKjotv67VnWd2+32mIR1kFZV5HDTG9sdgnwDOB2d39Eyn4D3gw8HVgGznb3+XDfy4E/CQ/9S3d/V7fP62eRw05ZPW+Rs2H2UZalcByUJ9ayxAmKtReS/397HWc/rw8DLXK4Ae8ELgIua7P/acDDw59JYAaYNLMHAH/G4WKK15rZFe7+732MtaPkzTbrtVkKoIlsFtPT02zZsqVrd1OkW0KIWiNR3a+i6lvicPerzezEDoc8G7jMgybPNWZ2rJkdD5wOfMbd74DWmudPBd7Xr1j7YWpqCjhceKzRaChptBGf4hitNrd9+/aO/yFFhmG9g93JhNCpZVH0ezigQ1dVOJuqU62qR3Z98yBxXNmmq+pKYI+7/3P4fBb4Y4LEcR93/8tw+y7gLne/MOU9zgXOBdi2bdujr7vuum4hrcuZZ54JwCc+8YlV26enp7n77ru58MI1oa15TfQ87X0G5eDBgxxzzGAXboz/vXtlWOcvzTDO6Xop1o2Znp4GYM+ePa1tG40z7T2T29KOWY92sR533HE97ap6Rvjnq8I/3x3+eRbBmMRGWco277B97Ub3i4GLIRjj6Hef6JlnntmqCRPvdooujtVqlfHx8VVdW8kLZ61WG2rfbRH7jfMq2t+haPF0oljXr90d4sk487QM0t4zuS3rnelZ9Oqcdp1VZWaPd/f4nePTZvYvwF9s8LMPAA+NPd8G3BpuPz2x/Qsb/KwNqVarrSTRqTDZ0tJS1/EQdb2IlFdyZtShQ4fYsmVLa1tc2oqHncqQpO1v97p4DMMYD+k6q8rMFoDzYl1KjwP+t7t3LaDfpavqTOA8gllVk8Bb3P2x4eD4tYR1soB54NHRmEc7/V46thdlj2u14K80rOQxrJkqeeauRxYXF9m1a1emY4epqLN/0ijW3oiPb3RKHPFjZ2ZmWkliI4UMNzKu0ctZVVkSx6OBS4BjCLqMDgKviKbOdnjd+whaDluBHxLMlNoC4O5vC6fjXkQw8L0MnOPu+8PXvgL4n+Fbvc7dL+32FxlW4qhUKkxMTLRWOovMzc21XlOr1Qpx8Svyf8aktFiLUm4hruzntKiKFmu0TG60YNXY2BjLy8scddRRnHzyyauOjSZ4RAkiPvFjeXmZsbGx1rHJpNFtXY/kcXn+Lww0ccTe/Ojw+IN5PmBQ+p040qbSzs3N5foFH/aFr2j/GSNp5yWKNa25XyRFPadpFOv6Rb+HkdnZ2dYa5PElceNfICuVSmuN8qi7u1KprFroqlMrBdKTTvTn7Oxsrv8fg25xPBj4K+Ah7v40MzsFmHL3d+T5oH7rd+Jop9Mv+LATRVLR/jN26s/tVK66SIp2TjtRrL2R7Krau3fvmu7Y6IIPwap/8cQyPj6euk54o9FojZHGE008WcRFS+9OTExkusYM+gbAdwKXAq8NnzeBDwCFShyDlHXWRFESxrDlnX8elasWKaJkyzi5DdKX1IXDv/vxhBF/36hLDGit4Z58z/gxcLhrbJCyJI6t7n65mV0A4O53m9k9fY6rcPLUhylaS6Nosqw/nbwjV+dUymIjS+qWZcndLF1VXwCeS3A3d83MTgP+2t2fOID4MhtWV9XOnTvb3uFctItdEZr/Wfp2oxbHsGPNoixxgmLth7LECcOZVfUW4BHAjcBxwPPd/ct5PqjfijLG0e7CWIRB3jL8kkfnb/fu3YWPFcpxTiOKtffKEicMeIzD3a81sycCJxPc1X2Tux/K8yGbRadv00WdFTRoWW9wio9xFK3lJrLZdU0cZvZN4I3u/rbYtivd/RkdXrYpdaqtv1l1WmMgLmthNxEZviyD44eAJ5nZJPBb7v5T4IT+hlVs8YtZNIgb6XaB3GzyJoRoKuP27dtLMR1XZDPKMsYxHw6K/xHBIPkLgI+6e63jCwes32Mc7S503fo489bj76dh9Meut8UQj7XIrY5R6OMuorLEWpY4obdjHJXuhwTVat39DQRlQK4iKDwoKdotaN/rRevLIvn3Tp6f+PP44+np6bb3f4jIcGVpcTzT3T8ee/4w4OXuvtHquD01jKVjs97hvJ4if/1QhCKHcfGZZlmW5Cxiy2MUvnEWUVliLUucMKDpuGb28+7+dTNL7ZLqVuRw0AY5HTd+ActyH0fSZq9VlVYlFFa3LFRypPcUa++VJU4YXOJ4u7v/ppl9PmW3u/sZeT6o3/rd4ojXnoHV34y7lf+Ov36zjXFEonMQ6VYVdPfu3avOa1HOYdIoXDiKqCyxliVOGFJ13KIb5OB4vIvljDPO4K677moVMWs3JTeqpDnMEuvD+CVPqyoMQRG35eXlVvXQSLyiaLJWT1QRtEhG4cJRRGWJtSxxwoBuADSz53R6obt/JM8HjZLkal1R1cpms0m9Xl/TBVOv11vll4v0bXkQ4kvp1mq1VWWio3MyMRGsCRYVbqtWq7g7ZtaqEBodM4zVzkRktU73cTyzwz4HNl3iaDabqy5czWaTO++8E6BV/jheIx+C0sdZyx6PovX+vcsyHVdkM+q05vg5gwyk6OL985OTk60FWZIXuHiLI6rDrwvexuj8iRRLlvs4MLMzzeyPzOxPo59+B1Y0s7OzrQVYotZF2s19UdKYmZlhYmJiTbeViEjZZalV9TZgDHgSsBd4HvDFPsdVSGn3YSQXHIp/O9Y3ZREZRVlqVT3O3R9pZte7+5+b2ZvYhOMbESUDEdnssnRV3RX+uWxmDyEoeri9fyGJiEiRZWlxXGlmxwJvBOYJZlQV/5ZeERHpiywLOe0OH37YzK4E7uPuB/sbloiIFFWWwfEjgDOBE6PjzQx3/5v+hiYiIkWUpavq48B/ATcAK12OFRGREZclcWxz90f2PRIRESmFLLOqPmVmT+l7JCIiUgpZWhzXAB81swrBVFwjKKt+dF8jk0KLVwgG1ZMS2UyyJI43AVPADT4qNdglVaPRaFWyzWpycrLjc4C5ubkNxSUixZKlq+obwI1KGqMvb9IQkc0pS4vjNuALZvYp4CfRRk3HFRHZnLK0OG4BZoF7AfeL/ciImZubo1arUavVVj2Ob5ubm6NarVKtVlvHQLAeSbQ/+T4iMlo6tjjCm/+q7v6H63lzM3sq8GbgCGCvu+9J7H8YcAlwHHAH8BJ3PxDuewPBjYcV4DPAq9Vd1n9RBeBoXZG0pXCTy7emLYerQXKR0dUxcbj7PWZWW88bh0nnrcCTgQPAl8zsCnf/auywC4HL3P1dZnYG8HrgpWb2OODxQHT/yD8DTwS+sJ5YJJ92F/1oe3wGVacEoZlWIqMpyxjHgpldAXwQuDPamGHN8ccCN7v7twDM7P3As4F44jgF+P3w8eeBj0VvD9yHoHvMgC3ADzPEKj3U7sKffJ71OBEZDVkSxwOAfwPOiG3Lsub4CcD3Ys8PAMm5ml8GnkvQnfXrwP3M7IHuvs/MPk8wMG/ARe7+tQyxygZlbU3EJVsiye0iMlqyVMdd79rjlvZ2iefnAxeZ2dnA1cD3gbvN7OeAXwC2hcd9xsx+2d2vXvUBZucC5wJs27ZtzWp8g3DwYHkKBWeJdffuoBhylnM5PT0NwJ49e1a9Ntq3c+fO1r68ynJeyxInKNZ+KEuc0NtYs1TH3Qb8PcGYgxOMN7w6GsTu4ADw0NjzbcCt8QPc/VbgOeHnVIHnuvvBMCFc4+5L4b5PAacRJJf46y8GLgbYsWOHb926tdtfpy+G9bnrsZ5Y07qiGo0GW7Zsaduq2Lt340u2lOW8liVOUKz9UJY4oXexZumquhR4L/D88PlLwm1P7vK6LwEPN7PtBC2JFwIvjh9gZluBO9x9BbiAYIYVwHeB3zSz1xO0XJ4I/F2GWKVHGo0GzWZzzcyqaB907orSwLjI6MqSOI5z90tjz99pZr/X7UXufreZnQdcRTAd9xJ3/4qZ/QWw392vAE4HXm9mTtCaeFX48g8RjKncQNDK+bS7fzzrX0o2rtMFP22fxjdENo8siWPRzF4CvC98/iKCwfKu3P2TwCcT2/409vhDBEki+bp7gN/K8hkyGN1aEO0GyEVk9GRJHK8ALgL+luDb/7+G22QEbbSLSS0NkdGXZVbVd4FnDSAWKYBuLYqkeAsjPiaSLLsuIqOjbeIwsz9tt49gPY7dHfZLybVreSS3t0soShgio6tTi+POlG33BV4JPBBQ4hhhWVse7abqZnkvESmntonD3d8UPTaz+wGvBs4B3k+wuJOMqHatjaj7aXx8fFW3VKPRWJMsms1m61gRGS3dquM+APgD4CzgXUDN3f99EIHJ8NXrdZaXl5mYmGBmZobZ2dlV93cATE1NsbKyQqVSYWxsjOXlZcbGxlrJRURGT6cxjjcS3NV9MXBqdBe3jLYoWYyNjQEwNjZGs9lsJYhqtcrs7GzrOAjW4piYmGitILi8vMz8/DyVSpblXkSkbDq1OF5DsOLfnwCvNWuVnjKCwfGj+xybDMnKygpLS0utVkS0rVKpMD4+3lpXPL5I09TUFJVKhX379tFoNFhYWGBsbEzjGyIjqNMYh74ubkJRF1TWC340trFv377WNiULkdGW5QZA2UTyXvSVJEQ2H7UqREQkFyUOERHJRYkjxeTkZGsAWEREVtMYR6jRaLSmk0YmJydXzRwSERG1OIAgQSSTRnyfiIgcpsSRgZKHiMhhShwxlUqFubk5arUatVpNdz6LiKTQGAesGcfQvQkiIu3pKzWkVncVGZRe/v7V6/VWFWORflGLg/S1JNTqkEHo9PuW3Bc9jyZyVKvVVaXtAZaWlta8PrJ7t5bQkd4wdx92DD2xY8cO379/f0/eK146vFsCWVxcZOvWrT353H5TrL23uLjIrl272s7KW69qtboqCfRKWaaXl+nfvwxxQvtYzexad9+R573U4iD7cqgiaXqdNIC+JA2RXtEYB0GiULKQIpqbm1szu69ara7ZFs0GrFarbd/r1FNP7UuMsvmoxSGyQVH3T7zlmuVxXHy1RTjciqnX660VGNPExzG6HaMxDukVJY6YPGMbIp20+/1Jm4gBrPt3LkpG9Xo9dX33eGt6cXEx9/uLpFHiiFGykI1Izn5qtz05lpY8vlOrBGB2dnbN+7d7jWYJSj8ocbD6P2SWbgWRTrK0NrJsT7aAkwkj7T10P5IMgqbjbtAoTMcrorLE2inObq2ASJ4vJmndqVm/4JTlnEJ5Yi1LnKDpuH2nsQ5Zj2RX0npaGJ32t9un31EZNCWOFPqPKOvRqSsp0ik55EkoWbeJ9IMSB6v/w8W/NbYbjBTpJjlWFrVgk8ektWyz3JCaZZsSifSLEkeKpaUlpqamOs6fF+kk/nuTZ2pu8mbULJM1kts3+xeeRqPBwsICKysrVKvVNZNeomPg8P0yyVIsm/0cdqPB8ZipqSlWVlaAYG2O+C9eO6MwOFZEZYk1GeegvuWvZ/ptWc4pbCzWfi28llbnaxTOaeEGx83sqcCbgSOAve6+J7H/YcAlwHHAHcBL3P1AuO+/AXuBhwIOPN3dv92POKOEkUwSKk8teWWdQdXuG3DWhBO1QBqNRtvWzWbsqpqamhp2CJtC3xKHmR0BvBV4MnAA+JKZXeHuX40ddiFwmbu/y8zOAF4PvDTcdxnwOnf/jJlVgZV+xQpBCyPZski7E1ckTbv7Jzp1U8Xv+O50wW+XALq1NjZTwojs27ev1eKo1WprzuHhGLRvAAAQvUlEQVTCwgJAq7QLQLPZbBWVrFQqjI2NASo02Uk/WxyPBW52928BmNn7gWcD8cRxCvD74ePPAx8Ljz0FONLdPwPg7n39F9y3b1/q9s34H0/WJ+vvSruxjHbvlZZE4gPqne4J2YwtDgi6lOLnpduaJ+1KtWzW85dF38Y4zOx5wFPdfWf4/KXApLufFzvmvcCcu7/ZzJ4DfBjYCjwB2An8FNgOfBaYdvd7Ep9xLnAuwLZt2x593XXX9eXv0snBgwc55phjBv6566FYey8e5/T0NAB79qzqkeUFL3gBAJdffvmq7dHx3/rWt1btj7ZH+0466ST27Nmz6v3bPc4aa9H1OtbkOcp6zroZhXN63HHHFWqMw1K2JbPU+cBFZnY2cDXwfeDuMK4nAI8Cvgt8ADgbeMeqN3O/GLgYgsHxYQ1SlWVwDBRrP0Rx7t27F1j7Dfdzn/sccHjMLPqGu2XLltb+RqPBrl27Wtujb7xmxi233MKuXbtWvf8NN9xApVJh69at3HLLLQDs2rVrzbTf5LflspxT6G2s0bmLROe+F5+xGc9pPxPHAYKB7cg24Nb4Ae5+K/AcgHAc47nuftDMDgDXxbq5PgacRiJxiBRFsiz6wsJCa/yi2WyyvLzM2NhY63l0cW82m6smYYyPj7OwsMDU1BRjY2NrloKdmZmh2Wy2+uLjyShKEtEAcXz/8vIyv/iLv7jmArrZqPupN/qZOL4EPNzMthO0JF4IvDh+gJltBe5w9xXgAoIZVtFr729mx7n7j4AzgMEXohLJaHZ2lqmpKRYWFpiYmGBlZYWlpSWazSZAa8AVgot8dP9ApVJhYmKi9Xx+fr61DWi9PhrAnZycpFartbbPzs4yOTnJ/Px8K0lMTEy0klU0vbxarW64W6ZMkvdhrDdhKNGk61vicPe7zew84CqC6biXuPtXzOwvgP3ufgVwOvB6M3OCrqpXha+9x8zOB2bNzIBrgbf3K1aRjYoPxkYXdaB18U4miGjGT/S6Wq0GtB+U7XQBa7eOePI1m2k9jrTS82m63WCphJFONwBu0CjcAFREZYm1LHHC5o21n62GUTinhbsBUERk2PLW8FL3VHdKHCIyUja6GqISRndKHCIydPHZX9GAPgSTB4BVM8yq1eqqO73jx0fiNwG2u5M+a4JRC2QtjXFs0Cj0cRZRWWItS5xQ3Fj7UZQwXm6km82+qqLGOEREyL6Oe5Qgss7CkkBl2AGIyOY2NzdHtVqlUqlQq9WYm5vj1FNPbT2Ofmq1Wuu4arVKtVoFaL0u+om6qdoVnuy0JHT8dZ3eY7NTV9UGjUJTtYjKEmtZ4oRyxppnvKHdiorJY+LaHZcsD9PuuDKe0yR1VYlIKSUv6Lt37wa6r1+SFF3oO1URjt+gmUwQ0b521YcloMQhIkOX/EafvMs9Swn6SFQ3LCrzUq/XWVpaas3QgsMlYJaXlwFWlWipVCqrao3JWkocIjJ0WWcxZZ1Cu7Ky0koKQGs8JJrGGyWW5IJOUaHKqDZYtG1qaqrtuj2bkRKHiAxduzGL5L7kkrnJWlONRmPNDCnpPSUOESmkdkUe2xUh1BTawdF0XBEpjLQpsHm6sTSYPRhqcYhIYURdTtPT061V+tqNZ6S9VgZDiUNEhiqtplTaPQfdBsZlcJQ4RGSospY3b3dPh1oag6fEISKFlOfeDRksDY6LSOFMT09Tr9dTu6E6DYJrgHww1OIQkUKId0Ht2bNn1RhHMhk0m03d2T1EanGISCGk1aGq1+utelLRMTMzM8zOzq650a/ZbK6qQyX9oxaHiBRCPGFMT09zyy23AMFaGVH9qeiO8cjCwgIrKyurSoooefSfEoeIDFU8YSwsLLQeR0vCTk5OtgoUzs/Pp64YGF9WViVH+k+JQ0SGLl6dFg6vNV6pVFqFCKObAxcWFhgbG+u6Vob0jxKHiAxV2oU/fgNgvEWiJFEMShwiUmhKFsWjWVUiIpKLEoeIiOSixCEiIrkocYiISC5KHCIikosSh4iI5KLEISIiuShxiIhILkocIiKSi7n7sGPoCTP7EfCdIXz0VmBxCJ+7Hoq198oSJyjWfihLnNA+1oe5+3F53mhkEsewmNl+d98x7DiyUKy9V5Y4QbH2Q1nihN7Gqq4qERHJRYlDRERyUeLYuIuHHUAOirX3yhInKNZ+KEuc0MNYNcYhIiK5qMUhIiK5KHF0YWYPNbPPm9nXzOwrZvbqcPsDzOwzZvaN8M/7h9vNzN5iZjeb2fVmVhtwvEeY2XVmdmX4fLuZzYVxfsDM7hVuv3f4/OZw/4kDjvNYM/uQmX09PLdTBT6nvx/+299oZu8zs/sU5bya2SVmdruZ3Rjblvs8mtnLw+O/YWYvH1Ccbwz//a83s4+a2bGxfReEcd5kZr8a2/7UcNvNZjbd6zjbxRrbd76ZuZltDZ8P7Zx2itXMfic8T18xszfEtvfmvLq7fjr8AMcDtfDx/YAmcArwBmA63D4N/HX4+OnApwADTgPmBhzvHwDvBa4Mn18OvDB8/DagET7+beBt4eMXAh8YcJzvAnaGj+8FHFvEcwqcANwCHBU7n2cX5bwCvwzUgBtj23KdR+ABwLfCP+8fPr7/AOJ8CnBk+PivY3GeAnwZuDewHfgmcET4803gpPB35svAKYM4p+H2hwJXEdwvtnXY57TDeX0S8Fng3uHzB/X6vA7kP98o/QD/F3gycBNwfLjteOCm8PE/AC+KHd86bgCxbQNmgTOAK8Nf5sXYf84p4Krw8VXAVPj4yPA4G1CcRxNcjC2xvYjn9ATge+EF4MjwvP5qkc4rcGLiwpHrPAIvAv4htn3Vcf2KM7Hv14H3hI8vAC6I7bsqPMet85x2XL9jBT4E/BLwbQ4njqGe0zb//pcDv5JyXM/Oq7qqcgi7HR4FzAEPdvfbAMI/HxQeFl1oIgfCbYPwd8AfASvh8wcCP3b3u1NiacUZ7j8YHj8IJwE/Ai4Nu9X2mtl9KeA5dffvAxcC3wVuIzhP11LM8xrJex6H+TsbeQXBN3c6xDO0OM3sWcD33f3LiV2FixUYB54QdpX+k5k9Jtzes1iVODIysyrwYeD33P0/Oh2asq3vU9fM7BnA7e5+bcZYhhJn6EiC5vWMuz8KuJOgS6WdocUajg88m6Bp/xDgvsDTOsQzzPPaTbvYhhqzmb0WuBt4T7SpTTzD+r81BrwW+NO03Snbhn1OjyToHjsN+EPgcjOzDjHljlWJIwMz20KQNN7j7h8JN//QzI4P9x8P3B5uP0DQFxrZBtw6gDAfDzzLzL4NvJ+gu+rvgGPN7MiUWFpxhvuPAe4YQJzRZx9w97nw+YcIEknRzinArwC3uPuP3P0Q8BHgcRTzvEbynsehnd9w0PgZwFke9pMUMM6fJfji8OXw/9c2YN7MfqaAsRJ+9kc88EWCHoitvYxViaOLMFO/A/iau/9NbNcVQDRT4uUEYx/R9peFsy1OAw5G3Qb95O4XuPs2dz+RYFD2c+5+FvB54Hlt4ozif154/EC+Ebn7D4DvmdnJ4aY68FUKdk5D3wVOM7Ox8HchirVw5zUm73m8CniKmd0/bGE9JdzWV2b2VOCPgWe5+3Ii/hdaMENtO/Bw4IvAl4CHWzCj7V4Ev+dX9DtOd7/B3R/k7ieG/78OEEyY+QEFO6ehjxF8ccTMxgkGvBfp5Xntx2DNKP0A/52g2XY9sBD+PJ2g33oW+Eb45wPC4w14K8EshRuAHUOI+XQOz6o6KfzluBn4IIdnWtwnfH5zuP+kAcc4AewPz+vHCJrWhTynwJ8DXwduBN5NMCulEOcVeB/B2MshggvaK9dzHgnGGG4Of84ZUJw3E/StR/+v3hY7/rVhnDcBT4ttfzrBzMZvAq8d1DlN7P82hwfHh3ZOO5zXewH/J/x9nQfO6PV51Z3jIiKSi7qqREQkFyUOERHJRYlDRERyUeIQEZFclDhERCQXJQ4ZOjN7oJkthD8/MLPvx57/64Bj+QMz+2pY6XTWzB6W2H90GN9FsW1fCCuLRjE/KPGa54UVVXeEz7eY2bvM7AYLKgNfEG6/j5l90cy+HFY1/fPYe7wj3H69BVWFq23i/zUzS7vDGTNbWv+ZSX2/z4b3KMgmo+m4Uihm9r+AJXe/cEif/ySCCqfLZtYATnf334jtfzNwHHCHu58XbvsCcL677095v/sBnyCYW3+eu+83sxcT3PT2wrCcxVcJ7r35DnBfd18KqxX8M/Bqd7/GzI72sNSNmf0NQXmZPSmf96/hey+m7Fty99SEsx7hXd/b3P11vXpPKQe1OKTQom/JZnZ6WLDtcjNrmtkeMzsr/IZ+g5n9bHjccWb2YTP7Uvjz+Dyf5+6f98N3MV9DUH4hiuXRwIOBf8zxlrsJypz/V/xjgPuGJUmOAn4K/IcHolbBlvDHw7iipGHha9Z84wvvEv5JlDTCO4H3hedhd+y4atiamg/P3bPD7bstXG8mfP46M/tdMzvezK4OW1M3mtkTwkOuIKgCK5uMEoeUyS8BrwZOBV4KjLv7Y4G9wO+Ex7wZ+Ft3fwzw3HDfer2SsGKrmVWANxEUjUtzaXhh3RVe3DGzRwEPdfcrE8d+iKCw420EJU0udPc7wtccYWYLBPWlPuOH63lhZpcCPwB+Hvj7lBgeT3CncOTNBIUkHxO+LvJfwK+7e41g7YY3hTG/g7BUSfj3fSFB4cEXE5TdniD4N1gAcPd/B+5tZoOu/itDpsQhZfIld7/N3X9CUBoh+uZ/A8GaBBAUJbwovPheARwddhflYmYvAXYAbww3/TbwSXf/XsrhZ7n7qcATwp+XhhfevwVek3L8Y4F7CKrtbgdeY2YnAbj7PeEFehvwWDN7RPQidz8nfM3XgN9Y867BOhA/ij1/PEFJCghKpbT+esBfmdn1BAv+nEBQiv3bwL+FCe8pwHXu/m8EtYzOCbsRT3X3/4y91+1hTLKJKHFImfwk9ngl9nyFoJQ0BL/TU+4+Ef6ckLjQYWZR6+CTaR9iZr9CUNPnWWGSgmCxm/MsqI56IUFhuz3QWrOD8HPeS5AY7gc8AvhC+JrTgCvCAfIXA59290PufjvwLwRJqsXdfwx8AXhqYvs9wAcIWlNJdxHUylr1kpTjziIYp3l0mKR+GHvdXoIVDs8BLgk/82qClea+D7zbzF4We6/7hJ8rm4gSh4yafwTOi56Y2UTyAHc/J0wqT0/uC79t/wNB0rg99pqz3P2/eVAd9XzgMnefNrMj7fD601sISoTf6O4H3X2rH66oek34nvsJuqfOsMB9CZLK18PxmWPD9zqKoPX09fC4nwu3G/BMgqKLSV8Dfi72/F8IupsgSBaRYwgG1w+FkwHiM8c+SpCsHkNYzdWCmWW3u/vbCbqzarFYfoag6J9sIkd2P0SkVH4XeGvYDXMkcDXwP3K8/o1AFfhgOFTxXXd/Vofj7w1cFSaNIwi6ft7e5TPeClxKUL3UgEvd/XozeyTwLjM7guBL3eXufmXY7fUuMzs6PP7LQCPlfa8mHK/wYLrkq4H3hgPeH44d9x7g42a2n2C8opWE3P2nZvZ5ghUO7wk3nw78oZkdApaAqMXxaOAaP7wSomwSmo4rMkLC6cIfd/fPrvP1FYIB9ue7+zcyfNYV7j67ns+S8lJXlcho+StgbD0vNLNTCNaOmO2WNEI3KmlsTmpxiIhILmpxiIhILkocIiKSixKHiIjkosQhIiK5KHGIiEguShwiIpLL/wf0Sv1VctXMzwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.plot();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Unfortunately, as the light fell on difference pixels in each quarter and the spacecraft changed temperature during observation, these light curves all have difference offsets. We can remove them by normalizing each light curve before we append it. We can do this with the lightkurve function `normalize`." ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%%capture\n", "lc = KeplerLightCurveFile.from_archive(target, quarter=1).PDCSAP_FLUX.normalize()\n", "for q in range(2,8):\n", " lc = lc.append(KeplerLightCurveFile.from_archive(target, quarter=q).PDCSAP_FLUX.normalize())" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJztvXuYZVV54P17C0q0KC5KE0cpPumOFsaJWJYt2HEyA3Q0GBKNYj4V8TZ2SCoyMUYzUz4OMkMPsR0Vg+JXDmnbYB4vGIyk40DAtAhG2xbsFCBiF6TR0GKCLdqkKS/V1vv9sfeqXrVr733Wvp2zz6n39zz11Nlrr8u7ru+6b1FVDMMwDKMsQ70WwDAMw+hvTJEYhmEYlTBFYhiGYVTCFIlhGIZRCVMkhmEYRiVMkRiGYRiVMEViGIZhVMIUiWEYhlEJUySGYRhGJY7stQDdYM2aNXrKKaf0WozKHDp0iCOPHJwss/i0G4tPe+lWXL7+9a/vV9UTO9kbjFTtwCmnnMLtt9/eazEqs3//ftasWdNrMWrD4tNuLD7tpVtxEZHvhNizqS3DMAyjEqZIDMMwjEqYIjEMwzAqYYrEMAzDqIQpEsMwDKMSpkgMwzCMSpgiMQzDMCphisQwDMOohCkSoyNTU1NMTU31WgzDMFqKKRLDMAyjEqZIuoj17A3DGERMkRiGYRiVMEVi5DI1NcXc3FyvxTAMo8Wsitt/jfLMzs6yuLjI7Oxsr0UxDKOlmCLpIv3cGC8uLvZaBMMwWopNbRmGYRiVsBFJF+nHXn0/ymwYRndpbEQiIttE5CER+UbGexGRD4jIfSJyp4hMeu9eJyL3xn+vi81GROT/isi3RORuEdnSlOyGYRhGOE1Obf0FcE7O+xcBT4v/LgRmAETkCcAlwBnA6cAlIvL42M17VfXpwLOB54vIi5oRvVnsLIlhGINEY4pEVW8FHs6x8hLgYxrxVeB4EXkS8OvA51X1YVX9IfB54BxVnVfVm2O/fwbsBsaakt8whWcYRhi9XCM5CXjAe94Xm2WZLyEixwO/BVzRsIwDhVMMMzMzHe1u3LiR+fn5pkUyDGMA6KUikRQzzTGPHIkcCXwS+ICq7s30XORCoikzxsbG2L9/fzVpa2TPnj2l5Dlw4EClcBcWFgCCwp6fn1+x0F53GlaNT9uw+LSbQYpP2+LSS0WyDzjZex4DHozNz0yYf9F7vgq4V1X/LM9zVb0qtsv69et1zZo11SWuiUcffZSy8lSJx/DwcLAfabu1mkjDNuVLHVh82s0gxadNcenlOZLtwGvj3VvPAw6o6veAG4EXisjj40X2F8ZmiMj/Ao4D/qhXQq9WhobsyJFhGOk0NiIRkU8SjSzWiMg+op1YwwCq+mHgeuA3gPuAeeAN8buHRWQzcFvs1aWx2RjwDuBbwG4RAbhSVbc2FYdBo8rJ+p07d9YoiWEYg0RjikRVX9XhvQJvyni3DdiWMNtH+vqJYRiG0UPsZPsqwk6pG4bRBDbxbRiGYVTCFImxAjuIaBhGEUyRGCuwD1kZhlEEUyQ1sHHjRjZu3NjR3q5du7ogTTWmpqY4ePBgr8UwDKOPMEXSRWzKaLCYmpqyPO0ylubtxBRJSwipINPT0z2pRKOjo10P0zCM/sEUSU0cPHgwaHqrKnNzc4WVSZFeXNr6yI4dOwqFZxjG6sLOkVSkH9YU5ubmgm7ynZqaqv3G3yI3Dq9m2pBOriOU1XFog4xpJDtJbZNvNWAjklVA2k2+SaamppidnW3VoUWbDzeM/sBGJBXxp4Lm5+fZsGEDIyMjqb26mZkZzjjjjFLhTE1Ncffddy97dn7WwdzcXKoSadtljd3oFYeGUUWGovFo62ig24TE30+rqakp5ubmGB8fZ/PmzZXS0ferqPtBzz9TJBXxp7VcQ1xmqqtIAZubm+PgwYPBjXzbRhkQHl+3JjQoFdDFv9NZnSLplNVwuuekP1lm8/PzjIyMLJlNT08zPDy8ZG/37t2VZa1Cltwuvo4y56Dy4tC2UfHU1BQLCwts3dqe+2pNkVSgSgHbuHFjoUVsf8Tg1jHqVBBZym9iYmLpd1ZlK9KQJCt5p3n5TvSqpxcablJxuAbPTTeGKsmsm5unpqbYvXt34ZHj3NwcGzduXFI4WSPSZDwANmzYwOLiIpOTk0t++Q15Mk+byiMX96p+OPI6LX65bbpz06meFXHTLUyRlMD1gqossndym+xZpo18gKWpNFeRXUEKUXJuXcTvhfq4hsKR1dNLmvuFOtkAJuPtL+67z/tOTEwsxWN+fj63h1n2avyi0xTJBgeiuGzYsCHoin0X79nZWSYmJpby0PnlNkRkNVDOftaIxi8Tfpql5Y2/Fubeu3xw/90UbFJBDQ0NrVA4yTxK5nFW/m3cuHHJ7ujoaGZnwikmX1ml+evScHFxMXPLuu8mWf6TG038PHfvRkZGOpZJ322Zhj1kRLV3797UMtkrTJEUpI5eUJ0kG4RQ/AYlS6klK4GrTH6P09+15hpBJ8vU1NSSfFkjMP99siKn9ZKT0xqLi4tLjV2yYiWVRJW1j+RamGt8fPmSU0GuofQbtbRG2P1fXFxc6u0me/nJeLk09xt6l4YuT118fQWV3FCRHN26UVJSXhd/363fcPtlwMeVj5ARU5py98tXXqfHvRsZGcksz9PT0yum73wWFxfZvXv3ss5BWr1Kpm/ZdROHP8IrsuElzW6vRiamSArgGt8Qik5dVSFkV1bd4aXhn3FJ9trSenvuG/KwvBfrpl3S3CSniOBw4xeaN86f+fl5du/evSKv0ho0v+eclQZ79+4l/uDaMpJ20zZdJKct5+fnmZ2dXdEwJ9PVz/dkAxrSKKW9TxvFpDXOSTM3Vea7c3myuLi4YsrLd3/w4MGlUb5rzCcmJlJHtGl1a3FxcVk9cPac4vBnENymmKw1TWfu0m90dHTJbMeOHcvyLqtT6cz9suYUmF++/M5SsuPkb9xJdkoeffTR1HDzlG2TmCIJxC94Ifg9lqZ7B1lyZU2/+Q1yKBs3blwRTvIApmvk/Arry+hXwLm5OaJvm6XjNwrz8/PLpkGGhoZSpzWS8qXNZftnavyGxDVcvuLwGwNHWnq6BtRVbmff2c0rN8l3rlfs4gnLp+/yRpBpfrupqNnZ2dQ8dPL7hE5NJfHTLCm37+8ZZ5wRfFuCL29S2YakQ1p6hdRjfzTtu/fLgh+/2dnZpTLk4ytS54+vMH2ylLVTKFl2knHrZtvjaNfezpbiF6oihM6jZrlNK2ydcH7u3r07s9D5UxKhpPXADx48uGLtJmstJ81tVq8qKZvrbSaffbOsSpjMg4MHD7K4uBgc9w0bNnS0k0yb2dnZ3HxLm/5JwymVKqNNf8oqqzx0apzKrgXmpXOnEY6vUNP8TUuzZFhV1jCzws4aDfpxTdopK1dI3uVR5haMstiIpEGKnBJPDvvHx8eXCnNIY+Zoav2mm1NnaWGl9dxDcL2zTpsjXMOVHK3lNWi+nawGJoReL5T2I71ap+zUoHeSq+w5sqI4ObpVtkyRdCBtrj6UKj3+Kg1TJ9p+pUvdFFnA7EYDFTKlYRj9hCmSAKo05HkLwGmLqQ5rXOqjG6Opuu8oM4w66FY70ugaiYhsE5GHROQbGe9FRD4gIveJyJ0iMum9e52I3Bv/vc4zf46I3BW7+YCkbZOpkbobiOQ2yLm5uaWFYjf327aTtEZnerFTxjDaQtOL7X8BnJPz/kXA0+K/C4EZABF5AnAJcAZwOnCJiDw+djMT23Xu8vxvHcldHckF6zqo4zr75GFEI5+0Mx+GsVpoVJGo6q3AwzlWXgJ8TCO+ChwvIk8Cfh34vKo+rKo/BD4PnBO/O1ZVd2q0d/RjwG83JX/Z3Vo+eQfqmppysWmx7mML5sZqptdrJCcBD3jP+2KzPPN9KeaN0IvGoQ1XudsXEYtjayTGaqbXiiRtfUNLmK/0WORCoikwxsbG2L9/f2Hh8g7MFWHTpk1s2bIlyG6vlcjRRx/N2rVr2bx5M5s2bQIIlt0wjPZRpu0rSq8VyT7gZO95DHgwNj8zYf7F2Hwsxf4KVPUq4CqA9evX65o1awoLl3VgLo20y+wcd911F2XC7wWnnnrq0vTb8PAwQN/IbhjGSrpRf3t9sn078Np499bzgAOq+j3gRuCFIvL4eJH9hcCN8bt/E5Hnxbu1Xgv8Tc+k7yOGhoZa95EqwzAGg6a3/34S2AmcKiL7ROSNIvL7IvL7sZXrgb3AfcCfA38AoKoPA5uB2+K/S2MzgClga+zmn4AbmoxDCENDQ63f/jkxMbFix5hRH5a2xmqm0aktVX1Vh/cKvCnj3TZgW4r57cAv1yJgTbjL/ty2237dNTU6OjowXyLsNlU+o2wY/Y7NdVRgdHR02Q6n8fHxvj5PkJR9ZmamdsXSz+dTktODVePSz2lhGD6mSCowPj7Ojh07lhrbJhpen6OPPrqUOzfSyJNtcnIySPaq6yxpYdS5dpPmV57/oWtHbvrSn8KamZlZ0ZkoQjIt6tp2HRKfOrd4V/WrSv4PDQ0tC79o/lcJN4/JycmgdCmadnn2e7kGaoqkBMnCW5Rdu3aVcrdu3bqlwlKkAQsZJYUqwLS1gGQv3ZdrdHSUXbt2LZPX74lPTk4GrS+4NHcNf1alSa5VjY6OrvDfyeLeOTe+jMk8npiYWOo0+PaS19knZS6a137+Tk5OLqVnsrEMaciyGB0dZXx8fIWd5HNoxyWtfBVp1Pz8SYvb5ORkbn77H7nauXPnivI4MTGx5G+eX6Fy+352woXrym7RsJKMj48vuXPxcWWll+t0pkhK0MuF9Z07d7Jr1y527NixrAJl9YBCRhqhPadkY+YaumTl9XGNjPvvPjTl23eNc7IRcQ3xrl272LlzJ+Pj40ubBiYmJpZVIj88vzF2jb/vrxtJOpl27NixLE398Jz/fhr6U5ju865puO3gvjx+3HwmJyfZsWMHExMTSw19WnidwnTvfVwc3H+XJp1Gh+vWrct858s9MzOzLA+Sm0/SOhdJP11ZcnHz/ZuZmenYcDv3afF37h/3uMctvfcVtS9DUqml4dJtYmJiRSfJt+OX1/HxcXbu3Llk14XlFLqTo1NdnJubY2RkhKOPPnppbdZvB9LqUTfo9TmSvsR9ItOnV99KHh0dXeoR+98wcYQokZARi7MzMzOTGldXAd37NDnyLqNMS9Mi601+zz/tm+cTExPs3r2boaGhFWmSlXcheZkmYzJNR0ZGmJ+fZ3x8fOnmgrTGzvczLWxfTj+NJycnl31vxYWdd7bJ98+XwX0ed2RkhC1btnDuuecuvd+5c+eKDQW+nH54fll0dpxb9+lY99uXxS9naZsXnJJKxjWZfsmwZ2Zm2L9/P2vWrFkR74mJidRbLFyZdumRVR6ducsDpwz87847kp8I9t/5YWR9zsC937NnT+r7JN1SKKZICuIaorxGcXJysvbvWrjCmWz0/AY4ZOfQ6Ojosk/WZjVaSXw7oWspIyMjmXYnJiZYWFhYVtHTGra0Z99eEfnrulU5K7xnPvOZDA8PrwjLz6OdO3fmyhHaEfF3Cbow/Ot10pRlXlr5jZ8vn19estwkzbIa3OTIMQ1fPr9z4JundVSKyJH0M5k2WZ/V7pQ3aUqtCE6GqampJQXgOkBwWCnMzMws3TqR5of/PfhuYYqkAKFrEnVuBfUr38LCQmX/0nZmFSG0957Xg/MbKv+b6k2P5pr2f8uWLamniPPC7dRjzBvpJtPXjXKKxnNoaGhZT9m5379//4pRbtbaS1qYeaOXTuUnL96+Es16XwV/hJJU2Gl0Cq9oRw0O1x8X/vz8fPCUuq+QuoUpkgL4FTdkyqjqeRJ/GiarUOTJsVovXwytRJ3SNjSsIuad3hUJK6vxLRJm6NSm3xvuFJYjTen0agrYJ6ROJUeUWX7ULYtbu/LDn5+fXzFCCum0dDONTZHk4M/5hs41ugKwY8eOoCF4Hsmtpvv37+fiiy8u7Z9PLypyHQ23709VO2Xs1oVLg5CdNnNzc8tGblkkp3+SZnXSS0UAK9caqtDEqKYqfvj+GluZjks3MEWSQ96CaDeoGl6dla0ooWsXq52m06jTFFEnZmZm2LBhQ1flcoQq0KbpdfhV1166gSmSHIpkYFpvO2utpImplCRF5rLrps0Fvt+oe667zEil7Hb3pspBXhzaMHXWTdoSX1MkOTSZOf5e+zrv5krOZddNHf76i7mdrrhuS0VpI02nzfT0NMPDwz0Z2aYp0LoUajcXoavQL3JCgCIRkV9Q1YcSZqeqathG5lVO3onc5K4Mf0Gt7P7vnTt3lnKXRrKhska9GnWlXz81MEnqWiergzaV4yY2Z3STkBHJl0TkYlX9NICIvBV4I/CMRiUbEEKmBdyeb6dMhoaGun7dQduURBsamk7nWqrQ9GaANjXY/Uzb6kVbCVEkZwJXicjvAE8E7gFOb1KofiZZ4LKmBZKL0W49ZnZ2tnTYVuiNfqHunXdtZ9BH9x0Viap+T0T+Dng7sAi8XVX784MbNVKlIGSdB5iamlo25VUHWXIOWkEuymqJf6ezJ1XpRTq2fRqom2nShvhC2BrJ54HvEX1MagzYJiK3qurbmhZutdHt06j9QFsqSlPU2SGpCydT1qG3Tu7avNGjTn+aIHkwsV8Imdr6kKpeF//+kYj8CtHoxBhwmu7NdpMyjVzd6wz93BCulhFcJywd0gmZ2rou8XyI6HvqRgP0qoC2rWK0SZ42yWK0j24rl6mpKRYWFti6dWtXwgshZGrr3wCNHx8DDAMHVfW4JgUzijPoDZ71BvuTJvOtTr/T/GpqFDloZThkRHKM/ywiv43t2mo1IfOsTS/md5Mq01aG4ajr2vXVWLYKn2xX1etEZDrEroicA1wBHAFsVdUtifdPAbYBJwIPAxeo6r743bsB91Wdzap6TWy+EXgP0dcdDwKvV9X7isajKquxsJSlLmUTWtHLhlf3gcF+2eJa14J/G+LSNE3HsV/TMGRq62Xe4xCwnsNTXXnujgA+BLwA2AfcJiLbVfWbnrX3Ah9T1atF5GzgXcBrRORcYBKYAI4CbhGRG1T1EWAGeImq3iMifwD8d+D1naNanjb0uo1q9EvetaGs+VfYJM2M/qDb5ShkRPJb3u9DwLeBlwS4Ox24T1X3AojIp2J3viJ5BvCW+PfNwHWe+S3xwv4hEbkDOAf4NJESOza2dxzwYIAsPaMNDYNh9JJufqmvF/TiVnBfyafR7TQPWSN5Q0m/TwIe8J73AcmrcO8AziOa/nopcIyInBCbXyIilwMjwFkcVkCbgOtF5MfAI8DzSspntIgQhVvkG+5NUaVjYJ2KbMqmTbfT0vIwnUxFIiIfJGcKS1X/sIPfkuYs8fw24EoReT1wK/Bd4JCq3iQizwW+Anwf2Ek0GoJoBPMbqrpLRP4EuJxIuSTlvxC4EGBsbKyjBs/DfeK2jB9l3U5PR8tQW7YcXlY6cOBAkNvNmzeXCrMMIWFlpYEfn5B0yrKTlKFsmofEZc+ePZl2XHyy/KkSxyYJyZ8qrF27NtX/EBnqpFN8rrnmmo4y9CJ/0ugUF9Woqe2WnHkjktsr+r0PONl7HiMxDaWqDwIvAxCRUeA8VT0Qv7sMuCx+9wngXhE5EXiWqu6KvbgG+Lu0wFX1KuAqgPXr12uRE7pJhoeHAQqd8q3qNsud/9wvvaO8/e4uPvfff/+y56L+lLFXBpGof5QlZ578IWWhrOxVykJI/lTxPyTeVepYEar6HyJnt+plngynnnpqRzt1kqdIPh6vUZTlNuBpIrKWaKTxSuB834KIrAEeVtVFotPy22LzI4DjVfUHInIacBpwU+zsOBEZV9U5ooX8eyrIaBiGYVQkT5F8jWjnFCLyQVX9L0U8VtVDInIRcCPR9t9tqnq3iFwK3K6q24luFn6XiCjR1NabYufDRNfXQ7QOcoFTaiLyu8BnRGQR+CHwn4vIVYYqPQvf7SBeoFiX7G1Y/wihX+Q0ltPG0+BN0u22JE+R+Gsczy/juapeD1yfMHun9/ta4NoUdz8h43snqvpZ4LNl5OknBn2ny2qkHzsKddAv52lCaIOcbVSKeYqk41mR1UzTowjr+baTphuSpsuVjX7bTb/GMU+RPF1E7iQamfxi/Jv4WVX1tMalMyrTLwWz7fK1HfcJgqmpqUZO9LdZgbatjLdFjm6Sp0h+qWtSDDBtK+S9wNLAGGSsfOcoElX9TjcF6XeK9OqSdtPc9EuhbJucoRsarPLXQ7+kY8hpcKM8hS9tNPLpl4plGEa7yWpL2qgUTZGUxBSFUYW8RiLUbtvol0/t9oLQtOnXOJoi6WP6ZVtlG2QwIto6jVrXWa1+oV86B6Hk3bV1F/l3bdmurQDaWlDKFORBKvyDEIckdcdpenqa4eHhnqVVv5S3tsvXDfJGJL8Z/3enzf8y/v9qYL4xifqU5CJuiF2jfpr+ImQI/dIA1kGvPh5mRLQlPTvu2hKR56uqf7J9WkS+DFzatHD9SK8zdNBoS0VpE0XSxNKvf+mnPAtZIzlaRP6Dqv4DgIj8CnB0s2K1C6uM9eLS0123vhoZxO/L94ucRair7g9i2viEKJI3AttE5DiiNZMDdOGiRCOdtiu1tstnDDb9Uv7aLl9RQr6Q+HXgWSJyLCDueyFGd3EXtbnvIVSlTEFua+FvqvHol0apbVgvvnu0JY06KhIReSLwp8CTVfVFIvIMYIOqfqRx6YxWYA1qMfo9nVx+b9mypfYPI/V72hjphExt/QXwUeAd8fMc0ZcJTZEYmVivtDkGYX2ll52TkLBDdl8ahwlRJGtU9dMi8nZY+mDVzxuWyxhgXAVu2zUPhpGkrYq4bYQokkdF5ATiw4ki8jyiBfdVQ5sKU9u3fBY5T2MYdbNaD9j2mhBF8lZgO9E3Sb4MnAj8TqNSGSvo5kVt3apgdZ2ctrudDKO3BO3aEpH/BJxK9FGrPaq60LhkRmtY7Q1qnmJta6+2jitwynZc2pYWq4Fel8OQXVv/BLxHVT/smX1OVX8zx5lRgF4XgqQcddDruITQlnQ3ukvI9KuViWKETG0tAGeJyBnA76nqz4CTmhXLaCvdanytkY8Y1HToFJ9BjfegEqJI5lX1FSLyX4Evicj/S86twD4icg5wBXAEsFVVtyTePwXYRrTu8jBwgarui9+9Gzg3trpZVa+JzQX4X0TrND8HZlT1AyHyDDpW+foPy7Nm6Pfv0PcbIYpEAFT1f4vI14EbgSd0dCRyBPAh4AXAPuA2Edmuqt/0rL0X+JiqXi0iZwPvAl4jIucCk8AEcBRwi4jcoKqPAK8HTgaerqqLIvILgXE1AulWwa964K1bFXXQGoI60m01NZJZWBocJkSRvNP9UNUdIvLrwOsC3J0O3KeqewFE5FPASwBfkTwDeEv8+2bgOs/8FlU9BBwSkTuAc4BPA1PA+aq6GMv0UIAstWAFx0jS1rLQVrlWI91oN3qd33kftnq6qn4L+K6ITCZefy7A75OAB7znfcAZCTt3AOcRTX+9FDgmPrNyB3CJiFwOjABncVgB/SLwChF5KfB94A9V9d4U+S8ELgQYGxurZevswkK0Wa3ubbjuFtxO/h44kH98p6p8IXLUmQZ58QkJJ8/O9PQ0EI16sghN91A65U8aneJZt4x5YSbDqpo/VQjJW0daHqe5L5M/eXRKA/feUWda1R2XquSNSN4K/C7wvpR3CpzdwW/JcOfzNuBKEXk9cCvwXeCQqt4kIs8FvkKkLHYCh2I3RwE/UdX1IvIyojWWX10RkOpVwFUA69ev1zruDHIXJtZ9/1AR8sLuhnxbt26t1b8sWUPCyYtvr/LKDy+5Kyit11glPUP8T8OlzcUXX9zRXVb6NZ2+eemSvLi0SP7XKW+nNAiRswq9bIeS5H3Y6nfj/2eV9Hsf0VqGYwx4MBHGg8DLAERkFDjP3S6sqpcBl8XvPgG4Ucc+4DPx788S3QNm0PvhrVEfNo1ajTakW123PPRDWcib2npZnkNV/esOft8GPE1E1hKNNF4JnJ8IYw3wcLze8Xai0YVbqD9eVX8gIqcBpwE3xc6uIxoNbQP+E9ElkkaXKXLxnZ08HzzaflVPNxjUeJUhb2rrt3LeKZCrSOLLHS8i2uV1BLBNVe8WkUuB21V1O3Am8C4RUaKpLfd9+GGircYAjxBtC3ZTW1uAj4vIW4CDwKY8OerECs5gMugNXhK7D627FLlluF/LYN7U1huqeq6q1wPXJ8z8XWDXAtemuPsJ0c6tND9/xOHzJYaRSpsqpPXe66UNFzNafi0nZPsv8bmOfw881pmp6qVNCWUYg0AvzrcUaeCsETTqIuSurQ9zeAvuVuDlwNcalstoiLb1pNomTyi2/mN0i+RUZBvLRsiI5FdU9TQRuVNV/6eIvI8O6yPG4NPLHm+bK5TRXur6bEGb6VXdCFEkP47/z4vIk4EfAGubE6n9+JlljVq99CI925B3bZChbsrEqQ31qZsyDEq+hyiSz4nI8cB7gN1EO7bqPZW2SuhUQNtQiQwjj34ro4M89dimvAj5sNXm+OdnRORzwGPdoUHDWC20qdLm0Xb5+pV+yf9eEbLYfgTRdttTnH0RQVUvb1Y0I5R+3qlTRp42xKENMrQFa2S7Q5vTN2Rq62+BnwB3AYvNimMYxmqiLiUU4k/Vzxb0A6391C4wpqqnNS5JH+Fnlq13RNQV39WSXv1Ct8txG/K/1zL0Y9sRokhuEJEXqupNna0aeXQqGG4X2NTUVE/vsGqKpuVusndbp8z9mn9Gu2hT+QlRJF8FPisiQ0TfbxdAVfXYRiUzDKN1tKnxqos2fJa33wlRJO8DNgB3qWrQt9qN+qhSyFdbz3duzi6CzqLJshDq52orj6uJEEVyL/ANUyLGIGCNWTOUTde68sHys7eEKJLvAV8UkRuAnzpD2/6bz2or2IMU315ctthW+knWfmBQp9FCFMn98d9j4j+jJsou6vbrwm/TYYyPjzfqv2HksZpHu7mKJD6MOKqqf9IleQzDyGC1NlRNf42zm+k5qOt4uYpEVX8uIpPdEsZYSZVCXtZtvzZYvT7Z3+bbZdvQgLUxXYx6CJnamhWR7cBfAY86w4BvthsvqHqMAAAgAElEQVRG6+hXJVmVfp72q0sJZin61Vom6iREkTyB6Or4sz2zjt9sHwSaHlIbxqBg5T+MEIXej21KyO2/lb/dbqTThoLSBhkMw+hvQm7/HQM+CDyfaCTyD8CbVXVfw7L1NU32Kvqxx2LkE5KnbVjn6AUhvfh+qQv9ImdRhgLsfBTYDjwZOInoNuCPhnguIueIyB4RuU9EplPeP0VEdojInSLyxVhpuXfvFpFvxH+vSHH7QRE5GCJHWebm5hqtvO5erV6SJsPMzExjCjAZVhvSIElZmbZs2bLkvgnGx8f7eq2jF9RVvtpWTqenp1slT8gayYmq6iuOvxCRP+rkKN46/CHgBcA+4DYR2a6q3/SsvRf4mKpeLSJnA+8CXiMi5wKTwARwFHCLiNygqo/Efq8Hjg+Q3WiIZA+6X0ZJg76TLUvOELl7EcduHtDLuka+7XnaD4Qokv0icgHwyfj5VUSL7504HbhPVfcCiMingJcAviJ5BvCW+PfNwHWe+S2qegg4JCJ3AOcAn44V1HuA84GXBsjRKFYIVydVGt26G+w0//pF8Rkr6cc8C1Ek/xm4Eng/0RrJV2KzTpwEPOA97wPOSNi5AzgPuIJIKRwjIifE5peIyOXACHAWhxXQRcB2Vf2eiGQGLiIXAhcCjI2NsX///gCRl7N27VqAUm4XFhY6ug2x43PgwIFC7qano9lEN+VSVIY890l3ZeJ74MCBwmnQDdJkColvaHxC/A91lzR3vx1F0jUtf8oQUu6ywqwb3/+y8Un6k8XmzZs72qmLQ4cOdS2sEEJ2bf0z8OISfqe18smLH98GXCkirwduBb4LHFLVm0TkuURK6/vATqKRyZOB3wHODJD7KuAqgPXr12uZL6MNDw8DlPqq2tatWxvxf82aNUF+h/qfZ6fIu7JhVUnjpgiRM0vusvke4i7Lju+f++0okq5pctVZb9JGSU3nf9L/suG0rZweeeSRDA8Pt0eerBci8s4cd6qqmzv4vQ842XseAx5MePIg8LI4vFHgPFU9EL+7DLgsfvcJoluInw08FbgvHo2MiMh9qvrUDrKUog2X9zU9RdHN21fT7LTxI1dtmFpoy9RUW+QoyyDdLuznRdZ6T6/yK29E8miK2dHAG4ETgE6K5DbgaSKylmik8UqidY0lRGQN8LCqLgJvB7bF5kcAx6vqD0TkNOA04KZ4zeTfee4PNqVE+om6Dk4WKYRtqFjdpF/im9z8UMZtHmk7/KrQxnTt5nfkB4VMRaKq73O/ReQY4M3AG4BPEX3sKhdVPSQiFwE3AkcA21T1bhG5FLhdVbcTTVG9S0SUaGrrTbHzYeBL8ajjEeCCWIkYPaAtFaItclSVoW75uznSK0Kb8qutDEoadbr99wnAHwOvBq4GJlX1h6Geq+r1wPUJs3d6v68Frk1x9xOinVud/B8NlaUNdLvQ1NVbbNN+dejOrqfVQBvj3UaZ8ug3eZsib43kPUTrF1cBz1TVRg//9QNWaIpRJr3qSuM25VGZOLVJ/qJkdUD6OU5GPnkjkrcSfRHxvwPv8LbaCtFi+7ENy9Za+vFjUFaJDUeREaaVm/bQ5rXLvDWSkOtTjC4xNTXFwsJC6hbQppXNoDYmTU/ZtW1KsA7aOt1ZJ73ovOV1TvthJiTkQKJhtIKmGrE2V9C20nSaVW08u9X4Vt3CPyiK2RRJzYQUmn6mDb2jQUjH1UCVMzxtb1jbUA/ahCmSmmh7wc8jq1JUrSx1NCSDUFHbGoe2yhVK2TJiH6yrH1MkBRiEQmUVpF56OQJtQ14WlaENMg8ivU5XUyR9wszMTGsuaBsU+qEBbqKBCA3TXUI4qPS68fXpxjX6TWKKpGW0qXCn0Va5jPbSzXM0bTysuhrqjCmSmml63rXftwlCsR5vU/Htl7TqJdPT0wwPD7cujdqwHtnEukw/Y4qkJvqlgBSZy+/lN8L7ZfdOEQa9MSlCL3v6IWUry47lYTqmSLpILxvmbrNaKpxtgy7GIHYQkvSi7Pe6LJgi6WO6XWBXi3LoBqFp2Pa0brt8ebRd9n6qb6ZIusj8/Dww+IcWO+Hi2PZdaKv10kEXz02bNhVyl1Wu6174zqMbV/S3ibYoG1MkXWRkZKSyH01sE0wWxirTDkXchizm1hXfZBybrHhVGtQ0t73aVJD1Fb5e0+tGswx1HvbNu3evV5giqYmqvde29Cx8duzYUdptm+LR7wzi2lrdtyh0i6Zv4O2XdEhiiqQkZTJ8fHy8KXGMGuhFJe7XhqMITcZxNaRfP2CKpI+pcvNo2z62VOX+o6Yak9Www6jNrHYl0U/xNkVSkjLTDVWHvKHzot2sgP1S2dtwMjkkrap0DnpN0/J0q/70E1lp0u0yYookgLaeoG7qs7R1N7obN24EVq65tHUxN4+6lEG3aFoJ+uZ1l/VkZ62tCrQMWXEIXYNp245HUyQ1UXdjPkiVxsgmJH/rXFtLU+p1b3Nu27Rp07TxWvpub9BoVJGIyDnAFcARwFZV3ZJ4/xRgG3Ai8DBwgarui9+9Gzg3trpZVa+JzT8OrAcWgK8Bv6eqC03Gowr9tOMmaxtwXuFvY/zKVtoQ+03Fd5DPFjUlfxsvaOyGDG2cpmtMkYjIEcCHgBcA+4DbRGS7qn7Ts/Ze4GOqerWInA28C3iNiJwLTAITwFHALSJyg6o+AnwcuCB2/wlgE9D1nG96B1ZWg1W0sJcZQhdZXHaHLIuQdmljmW2Vbaj4ddOPcZqbmwsuM0XiF1rHpqammJubY3x8PPfeuPn5eaampirdLVdX/tQ9iul1uWlyRHI6cJ+q7gUQkU8BLwF8RfIM4C3x75uB6zzzW1T1EHBIRO4AzgE+rarXO8ci8jVgrME4ZNLtveJuXvTiiy9eUXHLHLarq+CFHLJ0DUJaY9PE6fG0BqHMaKvuytmUP02uUziaHnk23RCOj4+zZ8+ewu6aGvU0TbePGjSpSE4CHvCe9wFnJOzcAZxHNP31UuAYETkhNr9ERC4HRoCzWK6AEJFh4DXAm9MCF5ELgQsBxsbGKi1OuZ5zGT/Wrl27wu3CwsIys+npaSBafHao6gp3Bw4cYGFhAVVlYWFhyb5znydfMoykDL7Z3r17AVi3bt0KO0l/nJzuOg0/Dg6Xfs5tMj4+RdLYuU2GnZZ2yfimxT/E/7T89OMTSlZ+OMqkw/79+zvGKy9tkmH78cmKtyuPe/bsYd26dWzZsiW3PGbJl2aerHfJZ19uVyey4u3XmzJpk1duQupzMo0deemQ5S/AoUOH2L9//4o6lZVWTdOkIpEUM008vw24UkReD9wKfBc4pKo3ichzga8A3wd2AocSbv8/4FZV/VJa4Kp6FXAVwPr167Xp3UFZPZfh4WGAZbuTkmZpdk499dQVZs6uiDA8PJzrPklemE52N+fqFmRDZHdyOkJkcBx33HErzIrkU9LtxRdfDICIMD8/z8UXX7yUHyFpnszDNNk6zUuHyp8M3/nrZMjzJ0vONWvWdCwLX/jCFzJlSZZd/wqb+++/f5m/fpmZmppidnaW+++/v2MaZckXUoaz/AJW1Ik0u53siMgKGVwYrmzl1Zu8uCTz15GXDlll1ndbpf7USZOKZB9wsvc8BjzoW1DVB4GXAYjIKHCeqh6I310GXBa/+wRwr3MnIpcQLdD/XoPyl6LpO53cXG9I2CH488tJ2nDVQ9kbBLq1CaBqPLt5B1gZsqZR20DagdEi98Yl7fp1oFM88+pNlv9weIowy13Z7b+9LjdNKpLbgKeJyFqikcYrgfN9CyKyBnhYVReBtxPt4HIL9cer6g9E5DTgNOCm+N0m4NeBjbG7nuIKS6cCUpSsgpHWSFYJ21dMLi5Z/qRVnrwK6xqhrHMqmzZtqvQFviInz8uE4fufFZe6qKqQiix4h8pStz8zMzOZ6ztpbjstoPv+dvLPNb55PfYiGz7ylFeStPdZZTfNTj/QmCJR1UMichFwI9H2322qereIXArcrqrbgTOBd4mIEk1tvSl2Pgx8KR5qPkK0LdhNbX0Y+A6wM37/16p6aVPxgPYtuiYrl2vYi/TCnfuNGzeyuLjYsXcVSl4j1K2K4S86d7KTphzS8nt+fr70KKdIQ1KGoovsZcLMKxszMzNL06FlKTqKTroJiVORTwe7vHbnbepSrk0vgvdq91aj50jiHVbXJ8ze6f2+Frg2xd1PiHZupfnZqkOUs7OzQLUr4kMr/8LCwoo50TrDCLGb1UOss2Es6k9a2HVObaU1lCFTJlnmZRvNkN57md1pWezdu5cf//jHqe+S6Zs3LVQ1L9ownVZma7ArM2nKo0watXn7e6sa5X4hLQPL9DTKFoTQaYy8nptvNjQ0tGzKrMzW2DR5QnrxdU4Zlc2PIvPqjqzpxLTpvzS7VRrHMkrbj2NSniz/1q1bt7TInmUnObWbRtqaXpY/TW6z3bt379KCeieSMifz2cXXV/B1dV6a2BLfNKZICpBXcbOGwnUXgryeaRllVnRhOlSJzczMsGHDhqVDYM4Myo02qlK1clZZHE0jWV6SFOnN1o2TTVWXRiRZvW1HJ8WZ5gccHtFPTEx0lCctrCQhI7HQ8pg3y5A8G+WPWvPWGqt0ApJ+5NnpNqZISpJsUJMFK69xrlKIkjLASuVVtKHeuHFj8NA61G9nL60yurCuueaapW2VnfxJWweocxqnG/gjvvn5+WVp07TseWEn7biyXPQAn6sTrjfv4rJhw4YlO7Ozs2zcuDGzh+/s7ty5c+ld3vqLGwH6coeSV5YXFxdL3dqQRkiehihiR940d5FZhToxRRJAsqHydzrV8fncukiTsxN+76noyCRPBreIPygkK2RajzSLkM0QyXdFGjF/mqXTonfduwtDcHVkbm6uUJnI24Hl3rv4ZE2jJv0IjffQ0FBq3S7T8w8pI3Xt1uvVx/NMkVSkU8bV1TNw2xfT5uDzwuh0v1AR/LDLNHShYXTyp8qctO9HlZ1qeQ2y25CQNWLKa1DLTMH5MmQpONfT9xvHTn77ayRZ+OngGvNkA1xmV2FZ9y6ebqosLZ/KKNO0BfSkP1nT2mnbx/PWnLJGjGkyudmEXn991RRJCebn55e2zCYLT5ktuCG44ay/+OkKdVaY4+PjzM7OBsuUXHxNkjcSS3PjFvH9OCQpo+DqPmxYZP+/o8yaiS93Mm2SdlxDWFXhOaqOnLPSvOz0T1LBufi5NZOQeLvymBa3tLWcdevW5Y7afMUwMjKSugvOn7bzR0R+WEXyysXXx4Wdt56yd+/eVPNejDjBFEkhkoU/rxJ1qgA+vZjbr9qou7gnv2/h/E1WOB+XNtPT09x///1Bh858Px3JBiRPwRRZ5MyzW9cOnZGRkcwRZJpZWqNa95brpL/uklBYWZ6TowS/gc6SM60BTy6yu5Ga32nJiqfr0BUlq276HakiaVrmMGiW/bQDwg4/HZPbstPWo7qJKZKayKosaQ1O2SG1qi7bvpg15C/awIQ25D7JRjzZwObtwkkjRFYXpqtoyWkc/zmkNwvNKO+sBjVPnmSP3M/TZGOZlDmtjPk9+04k7Tj5165du6Tos+yG+Ov35oeGhpbtVMs6GxGqqJ1/nc75zM3NLV3KmAw7Lb+S75KbR5KjkTR8GbLWcNJGU66Tkbc7022E6NS56hamSArgF6aDBw+mZlpyeqjqid+s8xmhPetO4RdtSN2UXt3kTc9lkbe5oNOUi5vj9tMxpIHuJJ8/x52VtnkHOpONmr+YHCqDT9qZERfnZFhpHZPkek8IWb3jkZGR4DU7fzSSlDNtR6BTTskw/bRJ3sCbpmx886pbu/NGOJ3O3YSmt6/Yeokpkh6SN3WSVYj9Rb5kQ5rmts7NAH7PL3Qnkh/HtPnbZE8ta1toaJidyJtDzpoPTxvhhMiS9C/kTEOZg5T+c1LO5PSQL7N7V2RePe/qkCqHQEPI6pz5frltw3lnRYp08rLSyH8OWYdMKqks/PqdNQU+Nze3NLVVR52oA1MkOeQNLdMWySCsYhQ5YOXbWVhYSD33kaws/rpFUp4q8+ppO0qKNB55doss2iZ7pnWT7PH689ZlSPtOehZZ+ZN25qRuyoyiq/baHXnxTeLyosghxjzy/MvaVOKnkTtz4pRYXvol5QopE2n15nGPe1zu+mO3MUUSQFqPbWJiIrcS5VX25Hx4UydV83p9ZRVKckdJMizf/6zhvKs8mzZtWvIvqZjTGoJO8+l5vbOsKcLQHVHJiy3z0s+VjTry0PVQqzYQfllIlr8Q6piq9XvbWYq5rk0RSf/T1kgcZdYV/HKQd8rfyRfaAUvb0pss365TmeW+F5giKUGRzAptADpVHHeOpO2UbfCyzh6kkdULSzuoFiJfp4OcfqMUOoURSqeDnXlbPB1p00zJcxR5ZG0K2Lx584pr14uMzNJ6892ay08rT1mNb8h0ahnyzp4kCZnq8xWK2/5b11psVUyRFKCb85DJ3rtPVuEucpK9LGV2eIX6C81taU0bVSbPCmTJ5HrueQ1B0k0e/lRXXddwlJHHn7rz8bf/ZrlJI+v0f5qdTttUfb+SdaDqzqTkLr9kmFXLXkj96HR2qdMmE/f5Z4cdSFzF1L1VL3QB3Lfb5FC40/SXkyOvp1pWviyFF9LTT9t512kqxvcjuTMpr2EKKQNZaZRsPNKmSvMW5EManyo7mPLcVCn7WQqgKCH5E7LzLiQdO23pLSOff1q+V1NaDlMkOdQxAqlzfSLPT5+6P8YTGm63/S6SjmmV3TXQnaYF/JFImZ5flXQr4jbvhHdyd1Wessmjrp5vFWUQYrfMWS1H0Xpfh+x55atMe9HtQ86mSHLImzvtRJHhbV3XYDh6cVK+ynDep8q20LrI2unWa6p0QMq6bWo6t+k0TcvDXq0xFq3XTU0fN4kpkhyS6xRNLWj550LSpiLcOyMia8qsl8rHEbKO4s/5Vzl70YtDaKuhHNa93lDXCMr3Z8uWLR0/wdBNTJEUIGTPt6NI4+buDEqbVy8zCqpyQHEQyJtOTKPuNCnSo+zWtFc3/OkX2hzfNsuWhymSAAZpVNDL3mzbqSN/e1lGinR0jHZR92aGbpdDUyQtI6m0er0/fDUwCB2EUIrGdTWMYJOspvJQF40qEhE5B7gCOALYqqpbEu+fAmwDTgQeBi5Q1X3xu3cD58ZWN6vqNbH5WuBTwBOA3cBrVPVnTcajadxldmlU2SXU9lHHli1bVhx4s0psGP1HY4pERI4APgS8ANgH3CYi21X1m5619wIfU9WrReRs4F3Aa0TkXGASmACOAm4RkRtU9RHg3cD7VfVTIvJh4I1AI61Pt3pj/sVvq7EHWJXVonx6Ec/VkrZGNZockZwO3KeqewFE5FPASwBfkTwDeEv8+2bgOs/8FlU9BBwSkTuAc0Tkr4CzgfNje1cD/4OGFEm3aOJAXie31kAYhlEXTSqSk4AHvOd9wBkJO3cA5xFNf70UOEZETojNLxGRy4ER4CwiBXQC8KNYwTg/T0oLXEQuBC4EGBsbK7WHPHk3TxE/Nm/eXMhNSFgHDhwIDr8fsPi0G4tPe2lbXJpUJJJilrx+823AlSLyeuBW4LvAIVW9SUSeC3wF+D6wEzgU6GdkqHoVcBXA+vXrNTkXH8LWrVsLuymL+xZ7p503ZeLRZiw+7cbi017aFJehBv3eB5zsPY8BD/oWVPVBVX2Zqj4beEdsdiD+f5mqTqjqC4gUyL3AfuB4ETkyy0/DMAyjuzSpSG4DniYia0XkMcArge2+BRFZIyJOhrcT7eBCRI6Ip7gQkdOA04CbNPqgwM3Ay2M3rwP+psE4dI3QS/QMwzDaRmOKJF7HuAi4EbgH+LSq3i0il4rIi2NrZwJ7RGQOeCJwWWw+DHxJRL5JND11gbcu8t+APxaR+4jWTD7SVBwMwzCMzjR6jkRVrweuT5i90/t9LXBtirufEO3cSvNzL9GOMMMwDKMF2Mn2lmDbcQ3D6FeaXCMxDMMwVgGmSAzDMIxKmCIxDMMwKmGKxDAMw6iEKRLDMAyjEqZIDMMwjEqYIjEMwzAqYYrEMAzDqIQpEsMwDKMSEt2DONiIyPeB7/RajhpYQ3QD8qBg8Wk3Fp/20q24PEVVT+xkaVUokkFBRG5X1fW9lqMuLD7txuLTXtoWF5vaMgzDMCphisQwDMOohCmS/uKqXgtQMxafdmPxaS+tioutkRiGYRiVsBGJYRiGUQlTJC1CRE4WkZtF5B4RuVtE3hybP0FEPi8i98b/Hx+bi4h8QETuE5E7RWSytzFYjog8VkS+JiJ3xPH5n7H5WhHZFcfnGhF5TGx+VPx8X/z+lF7Kn4aIHCEi/ygin4uf+zku3xaRu0RkVkRuj836sqwBiMjxInKtiHwrrkMb+jU+InJqnC/u7xER+aO2xscUSbs4BLxVVX8JeB7wJhF5BjAN7FDVpwE74meAFwFPi/8uBNr2mcWfAmer6rOACeAcEXke8G7g/XF8fgi8Mbb/RuCHqvpU4P2xvbbxZuAe77mf4wJwlqpOeFtJ+7WsAVwB/J2qPh14FlE+9WV8VHVPnC8TwHOAeeCztDU+qmp/Lf0D/gZ4AbAHeFJs9iRgT/z7/wCv8uwv2WvbHzAC7AbOIDpIdWRsvgG4Mf59I7Ah/n1kbE96LbsXhzGiyns28DlA+jUusVzfBtYkzPqyrAHHAvcn07hf45OIwwuBL7c5PjYiaSnxVMizgV3AE1X1ewDx/1+IrZ0EPOA52xebtYZ4KmgWeAj4PPBPwI9U9VBsxZd5KT7x+wPACd2VOJc/A/4rsBg/n0D/xgVAgZtE5OsicmFs1q9lbR3wfeCj8dTjVhE5mv6Nj88rgU/Gv1sZH1MkLURERoHPAH+kqo/kWU0xa9U2PFX9uUbD8zHgdOCX0qzF/1sbHxH5TeAhVf26b5xitfVx8Xi+qk4STYu8SUT+Y47dtsfnSGASmFHVZwOPcnjaJ422xweAeM3txcBfdbKaYta1+JgiaRkiMkykRD6uqn8dG/+riDwpfv8kot49RL2Okz3nY8CD3ZK1CKr6I+CLRGs/x4vIkfErX+al+MTvjwMe7q6kmTwfeLGIfBv4FNH01p/Rn3EBQFUfjP8/RDT/fjr9W9b2AftUdVf8fC2RYunX+DheBOxW1X+Nn1sZH1MkLUJEBPgIcI+qXu692g68Lv79OqK1E2f+2njHxvOAA27Y2wZE5EQROT7+/Tjg14gWQG8GXh5bS8bHxfPlwBc0nvDtNar6dlUdU9VTiKYavqCqr6YP4wIgIkeLyDHuN9E8/Dfo07Kmqv8CPCAip8ZGG4Fv0qfx8XgVh6e1oK3x6fVCkv0tW1T7D0TD0TuB2fjvN4jm1ncA98b/nxDbF+BDROsOdwHrex2HRHxOA/4xjs83gHfG5uuArwH3EQ3Zj4rNHxs/3xe/X9frOGTE60zgc/0cl1juO+K/u4F3xOZ9WdZiGSeA2+Pydh3w+D6PzwjwA+A4z6yV8bGT7YZhGEYlbGrLMAzDqIQpEsMwDKMSpkgMwzCMSpgiMQzDMCphisQwDMOohCkSo9WIyAneDaj/IiLf9Z6/0mVZ/lhEvhnfrrpDRJ6SeH9sLN+VntkXRWSPJ/MvJNy8XERURNbHz8MicrVEt/LeIyJvj81Tb1KO330kNr9TottvRzPk/20ReWfGu4PlUybVv793N9Mag49t/zX6BhH5H8BBVX1vj8I/C9ilqvMiMgWcqaqv8N5fAZwIPKyqF8VmXwTepqq3p/h3DPB/gccAF6nq7SJyPvBiVX2liIwQHao7E/gOcLSqHoxvP/gH4M2q+lUROVbjq3RE5HKiq1y2pIT3ldjv/SnvDqpqqgIqg4i8DhhT1cvq8tNoLzYiMfoW14sWkTNF5BYR+bSIzInIFhF5ddyDv0tEfjG2d6KIfEZEbov/nl8kPFW9WVXn48evEl1D4WR5DvBE4KYCXm4G/jfwEz8Y4Oj4WpXHAT8DHtEIN2oYjv80lsspEYndrOgdisg48FOnRCT6jsrOOB02e/ZG49HW7jjtXhKbb5b4+zjx82Ui8oci8iQRuTUebX1DRH41trKd6FS2sQowRWIMCs8i+lbIM4HXAOOqejqwFfgvsZ0riL4d8lzgvPhdWd4I3AAgIkPA+4A/ybD70bihvThu7BGRZwMnq+rnEnavJbpw8HvAPwPvVdWHYzfLblLWw/dKISIfBf4FeDrwwRQZnk90jb/jCqILDp8bu3P8BHipRpc5ngW8L5b5I8RXc8TxfSXwceB8oqvzJ4jyYBZAVX8IHCUibbvx2GgAUyTGoHCbqn5PVX9KdE2EGxncBZwS//414Mq4Md4OHBtPLxVCRC4A1gPviY3+ALheVR9Isf5qVX0m8Kvx32vihvj9wFtT7J8O/Bx4MrAWeKuIrIOVNymLyC87R6r6htjNPcArVvgafbvi+97z8zl8h9Nf+tED/lRE7gT+nugq8ieq6reBH8QK8IXAP6rqD4DbgDfE047PVNV/8/x6KJbJGHBMkRiDwk+934ve8yLRFeMQlfcNGn95TlVPSjR8iIgbPVyfFoiI/BrwDqK1BhfGBuAiiW4Gfi/R5XlbAFT1u/H/fwM+QaQojgF+Gfhi7OZ5wPZ4wf18oq/8LWh0K++XiZTWEnr4JuVzEuY/B64hGm0l+THR/V/LnKTYezXROs9zYqX1r567rcDrgTcA2+IwbwX+I/Bd4C9F5LWeX4+NwzUGHFMkxmriJuAi9yAiE0kLqvqGWMn8RvJd3Bv/P0RK5CHPzatV9f/R6GbgtwEfU9VpETlSRNbEboeB3wS+oaoHVHWNqp4Su/lq7OftRNNZZ0vE0URK5luSfpPyt2J7T43NBfgt4Fspcb8HeKr3/GWi6SmIlIfjOKLF+oV4c4G/M+2zRMrruURfgESinWsPqeqfE01/TXqy/DuirzAaA86Rna0YxsDwh8CH4mmbI4Fbgd8v4P49wCjwV/FSxz+r6otz7H2hsG8AAADuSURBVB8F3BgrkSOIpor+vEMYHwI+SnRbsgAfVdU7ReQ04GoROYKoA/hpVf1cPE12tYgcG9u/A5hK8fdW4vUOjbZqvhn4RLyA/hnP3seBvxWR24nWO5aUkqr+TERuJvoq5M9j4zOBPxGRBeAg4EYkzwG+qoe/HmkMMLb91zBWCfH25L9V1b8v6X6IaMH+d1T13oCwtqvqjjJhGf2FTW0ZxurhT4m+cVEYEXkG0bdVdnRSIjHfMCWyerARiWEYhlEJG5EYhmEYlTBFYhiGYVTCFIlhGIZRCVMkhmEYRiVMkRiGYRiVMEViGIZhVOL/B1wlcBTb3wY8AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.plot();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now that the light curves are normalized they don't have significant offsets. However, there are clear long term trends in some of the early quarters of the data. We can remove this using lightkurve's `flatten` method before appending each object." ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%%capture\n", "lc = KeplerLightCurveFile.from_archive(target, quarter=1).PDCSAP_FLUX.normalize().flatten(window_length=401)\n", "for q in range(2,8):\n", " lc = lc.append(KeplerLightCurveFile.from_archive(target, quarter=q).PDCSAP_FLUX.normalize().flatten(window_length=401))" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJztvXucZkV54P99Gjpo01zUIa6hWWGijTEB23YEJ252gY6KmmgUE2+gEgnZXlmN0eyOH3dkl/kRxxU13n7tsgSD+XjBYHQnLgZML4LRcQQnDYhIww4xjJjgiA6O7aXHfvaPc6qpPnPOeevc3vect5/v5zOf6fecujxV9VQ9dTtVoqoYhmEYRllGBi2AYRiG0W3MkBiGYRiVMENiGIZhVMIMiWEYhlEJMySGYRhGJcyQGIZhGJUwQ2IYhmFUwgyJYRiGUQkzJIZhGEYlDh+0AP1gw4YNeuKJJw5ajMocPHiQww8fniKz9LQbS0976Vdavva1r+1T1eN6uRuOXO3BiSeeyC233DJoMSqzb98+NmzYMGgxasPS024sPe2lX2kRkW+FuGtsaktErhSRB0Tk6xnvRUTeJyL3iMhtIjLtvXu1iNwd/3t1/GxMRP63iHxTRO4Qke1NyW4YhmGE0+QayV8AZ+e8fy7wxPjfhcAcgIg8GrgYOB04DbhYRB4V+7lMVZ8EPBV4pog8txnRDcMwjFAaMySqehPwYI6TFwIf0YivAMeKyOOA5wCfV9UHVfX7wOeBs1V1SVVviMP+GbAbmGhKfsMwDCOMQe7aOh64z/u9N36W9XwVETkW+G1gvmEZDcMwjB4McrFdUp5pzvPIk8jhwMeB96nqnszARS4kmjJjYmKCffv2VZO2Bezfv3/QItSKpafdWHraS9vSMkhDshc4wfs9AdwfPz8j8fwL3u/LgbtV9c/yAlfVy2O3bNq0SYdlt8awpMNh6Wk3lp720qa0DHJqawfwqnj31jOA/ar6HeA64Nki8qh4kf3Z8TNE5P8DjgH+aFBCG4ZhGGtpbEQiIh8nGllsEJG9RDuxRgFU9UPAtcDzgHuAJeD8+N2DIrINuDkO6pL42QTwVuCbwG4RAfiAql7RVBoMwzCM3jRmSFT15T3eK/C6jHdXAlcmnu0lff1kKJidnQVgbm5uwJIYhmEUw87aMgzDMCphhsQwDMOohBkSwzAMoxJmSNYRs7Ozq2sxhmEYdWGGxDAMw6jEujhGviz93Ellu7UMoze2u7GdmCHJYXFxcdAiGIZhtB4zJDkcOHBg0CLUysLCwqBFGCqsd2wYEbZGYhhGELZZw8jCRiR9ZNA92JWVlYHEW5ZB59cwsB7ycD2kse3YiKQlWG/PMKrRdB2qEv6w128bkfSRqov3W7ZsYXR01HpeDdNkD7do2NbbXkuZfHB5uG3btrrFMWLMkNTAzMwMAPPzay9sLNIIpC2EV2lEQvw21Ui1qSFukjpl6JWufqS7TXlbBT8ds7OzLC4uMjk5Gewnj8XFRWZnZzufR3VjhqSPLC0trf6dVFy3ftGvBsONjpIVLC3+ojL5I6+q6VnvFdd1MKampkqHUXceFi3TrI5WUxQZ+Q/zdFM/MUMSSB0NfNpid5YiZ/WkshqFMvItLS2VqnR1xF0XbepFNzH16HTGlVMd6Q0N4/d+7/cQkUN0cPfu3YXich2ossYsT++qGlpf/5eWltZ1p6UKZkgCcD2qvOGxU0KfkCF1siFfWFhgbGys8FC8jo8ne00DZMVRJe683uri4uKaUVwa/pTgoKfU9uzZQ3zh2iFs3rwZKN/gJfPB18m0/J+cnGRubi4zDxcXF5mZmVn179yMjY2lxu/cZL3PI2kMHUWmk9LcLy4u1r4TcWFhgZmZmVV99ON0U2W7d+9m8+bN7Ny5s1A6QvHDS5azK9e2rZeaIQnAfZi4e/fuNcYirRD9CpfV4/fDSH706FeMAwcOrPb+nB/XyLvnTtFcOHmjlay0+Q1Kmp+kcelVcXxZkiT9pk33uTiXlpZYWVnpmYdp7/xK5545I53WSCQbjCTOYGX5AfjRj36UKacr17wpQzfa9PPA4fy7NPh64/IJYGRk5JDGfmVlhYWFhdWwXVp8d8kGeWZmZjU9i4uLq/H58Tpdy8sT39AnDVpoB6RXZ8Lh4nR1wp/C9dMFUcdldnY29aNjP75kPXP1YGVlZU1dy5s+nJmZYWlpabUTkcyjZGcqOUpyz9Lyz9XdQRsUMyQFSRtKu4rslG5kZGRVQbMqQV4lyup9uud+3MmK4MLNC39hYeEQg+VXeL/hmJ6ePkSuZNj+1EVWPE7hk5sK/LWhkF6va6xcg5psALN64L4svpwLCwts3rz5kFFg1hqW/3737t2MjIwcEqbfwDo580gz2n4Z+KysrBySvqRsS0tLq4YjTT9WVlZWdTQ5TeVG1r6/Xic85Bl6X7aVlRVOP/10du3atRoXwOmnnw5EuubSPTIyws6dO5mdnV0zojlw4MBqefly+Q25nz9purC0tLRap9JYWVlZbfxD8yDZaUzTHz+vfXn9jldST93ffmdhdnaWPXv28OMf//iQcAeFGZKCuMJ0w1v/WdINpE8VpPUqk6T1SrKmCHwOHDiwpiFIa8zSpgP8Hm1SDldJnJtkr9TvLaeF6afBMTMzQ3TbMqt+/Urj+0+baltZWVkjh2uMfHx5fFl8Y+Ma1SRZo7Pku15lktbwuzB9owH5vc9kmI603X5p5eFGJSFTQaEjBWeEkwbaxTM+Pp7qz6U9KYsvn+vxZ41oknUnbWTvx7F79+7VNZ+k7qSR9d7Pm+QIEjjEQKUZc5dOf5TowssqW1/PfKPi8PXfGeF+In5lHlY2bdqkt9xyS2F/aY1TGVwPrK7wejEyMtKXr9hdQxF6Jpkvl2u8Q+Scnp5eMw2QVpGyyMt7J48v165du1ZHGy7uXj2+Xbt2MTMzsyYfxsfHU+V0Izw/zCJ5MQyMj48H6Uwy73vpdVae+/57deDqIKnnvjwh+lQXfr0pi4h8TVU39XRnhiSbuhp+10PolyEZRsbHx5mfny+ch0WNeJlG3Z+SCSG0ITW6T786dVlxVx2ZhBoSOyKlD6yXnmaTpO2KC6Gon5WVlcLltXv37uAFYRi+U6WNbAZZ9930YD8wQ2J0grKLisl56KYosy3WMJqmX3cqmSExhpp+9Qh7ffdjGMOMGRLDMIwhpciUaxUaNSQicqWIPCAiX894LyLyPhG5R0RuE5Fp792rReTu+N+rvedPE5HbYz/vk6xPiQ2jj9jtk0Yb6deIvOkRyV8AZ+e8fy7wxPjfhcAcgIg8GrgYOB04DbhYRB4V+5mL3Tp/eeFXIu37AqN7tOGDLcMYZhptKVX1JuDBHCcvBD6iEV8BjhWRxwHPAT6vqg+q6veBzwNnx++OVtWdGu1b/gjwO03JbwuohmEYvRl0l/t44D7v9974Wd7zvSnPDcMwjAEx6CNS0tY3tMTzQwMWuZBoCoyJiQn27dtXWLiTTjqJ22+/vbC/NMrEbxiGUZV+tD2DNiR7gRO83xPA/fHzMxLPvxA/n0hxfwiqejlwOURftm/YsKGwcKOjo4X9ZFEmfqM7TE1N2VqM0Ur60fYMemprB/CqePfWM4D9qvod4Drg2SLyqHiR/dnAdfG7H4rIM+LdWq8C/ldTwtV9NLM7rsMYPgZ9jLdhDJJGRyQi8nGikcUGEdlLtBNrFEBVPwRcCzwPuAdYAs6P3z0oItuAm+OgLlFVt2g/S7Qb7JHA5+J/ncCu9TQMYxhp1JCo6st7vFfgdRnvrgSuTHl+C/BrtQhoGIZhVGbQU1uGYRhGxzFDYhhDjn1YazSNaVgLcBV92Cv8+Ph4JzccjIyMMD09Xap82lCmY2NjmbcVDjsjIyO1lUGVcIY9/wev5QZTU1NMT0+vuYu9DL6i+3et14GrkL3icI1u8tn4+HjqCbl1NrR5lXV6ejrTGPjPnKx+WK5cXIM8Pj6+6q4XoWWa1eCVKcdkOJOTk6VOJ65SNlm6USaeXbt2BeXD+Pj4Ie7GxsaYmpqqRc+mpqZWy74XvhunJ013KgbZSTND0kfm5uYaVSZ3pMv4+Dhzc3OlGqE0P7t27WJsbGz1H0SVIm3La1rDOTY2ttqQzc7Orjaap5xyyqoRzcLll99oJxt6R0hjmWxUxsfHVw3Erl271two54zP3Nwcc3NzzM/PMz8/z+TkJFNTU8zPz2fK7MvoGh/fWCQN7s6dO9fE7fw7nSmqN348aeXkG8QsfvVXf7WnsZyenk4NJy1O3xD7cvo6k9ZhceE5v0n/u3btYteuXczPz6+6c2TpRJkRwtzc3GrZuzjTdDfZKXSGfOfOnav5lSVHnlx5I+Okser3SNgMSQlczyetJ9trCsQ1nGkFnlb5XKUaHx/nyCOPTK24fo9/enp6tYHLMyYhPSu/QrvKMD8/v9r4unDSwvLTNzk5uSZtyUYzD9f4uHhdOrMaiLSeYFqYzki4NPlGIfk7iTMsLj4X55FHHrnaMDoZ/XQ7Q+zuDg9lampqTcOUbIySDe/OnTtXw/fjSZbV5OTkIe/9MJ2b5LNkXszPz7Nz50527dq1Gn7aVnc/r527pB65tKZ1SJI6mExfGq6sdu7cySmnnLIm/UmjlNQXZyx64XccnPH3DV9SHpenrm748vh659oTJ/PU1BRzc3OH5I1fLq6+9PucwEF/2d56pqenV79YdkoxOTnJ4uIiU1NTq8eH+w2jq0Rzc3PMzMysCc81LDMzMywtLa0qh8O/4zlpBJaXl1fdzs7OrsrlFCvNEM3Nza25rzypdA6XJj9e/3Y1P2y/cXGK79I5NzfH7Ozs6igk5EM9P5/c/Ql5aUqmxa9YY2Njq1fZ+mlKS4crJ7+8iuLHvXHjRkZHR1fzwJE0Sn5eJe+ST8v7JC58v/wXFxdZWlpaLVNXRmnltri4uFo2s7OzLCwsrJaX/27fvn1s3bp1TVpdnFmdkLQ8d/UmL3+THYO0uuMzPz+f+V1W3jtfPj8//HdOf2ZnZw+ZRfDrn8MfceelKS1trgxcvkOUX0tLS6vl4NwuLi4ekq7p6ek1+e3HuXv37r6NTMyQBFKkQPwKk9WzTatwflzJhjSpQHNzc2zevPmQ+NIYHx9frRzJ3kxaw9pL9rq+4t6+ffshxzfkjVSSFTPL2DgjDdlpCKFoOv30hDaaIVNC7plvdJJGKKlPyXDyDGjVtblepJVByDRkU2WXjNs3us6wFomnzIfG/lRvUqZkJ2BpaemQe0WOPPJITj755DWGyJejlyGuGzMkPehVIFWGkGk9dhdeSCMWOj3k9zTzyDJaoXEkfxdpiEMajVD58ox0G0hOd4WSHKEme/ppepo22srrnKTJ6ndYnJ+sxtZ3k7f21UuGPLnL4gx9MuxkA5x8Fmpos+TL6xyE1LVe8Zc1ZnVihqQAyZ6L/6woWUrXxN3faYo2yLOhnCzLy8urU0F1U9SQlaVqHFV63UldKRtWrzT4xq5KZ6OKDHX5KUOvePpdl+bm5rjgggt6ylJFt4pihiQAVyBplacrh/XV1QBkrSeEroXUIcN6ZVC61kTnxtGGDk1TlE1bV9oUHzMkBehCD7cNMjSdhl4GKTT+NuR1CGWMd52NZN7UTKh7R5WNDXWFkTd9VjQMI8IMSceoqxIZxiAJ0WO3gDwIna3D2FQlL+60jSqDxAxJSQYxL5rcjrneMaNYjTbnXxsWkPtBXdO9gy5LMyQ1McieS79I7nevSkivqg352gYZ6mKY0lKW9WCg+o0ZkiGiH41EHbuhnP9+3mM/jA1oSJqqpHfLli2N7aorQ4gBGMZy7gJmSAZA08oeEr5VuMGyXvO/TTv3+pn3vcq7rCxt0SMzJDXTdMEOWmH6RRsamjbIAIeWeYhcbdeTKvL1M21tz8e2YIakYdrSYzAGT5oumH5ks97zpEvpN0NiGD3oUoXuN74h7IpRbIN8XcmrUMyQ1ESVud+uKNUg9/PXHU5X8rxuQtLblim9JEVOTxgk61G3zJAMgKbXT/yGIFSp16PyD5I6Fp2bLqu6PnobJt0aprTUiRmSkrRZoerquQ0yjW3I1zbIYAwHobpUtM61RUfNkNRMWwrWaB9Fz6wqEk4VbLt4NSxvzJAEYYe7VaOuitaVCtsVOY1uMjs7y/LyMldccUWuG+ifDvY0JCLyi6r6QOLZyap6V4Dfs4H3AocBV6jq9sT7xwNXAscBDwLnqure+N07gOfHTrep6tXx8xngnUT3zR8AXqOq9/SSpc10ucFpQ6NZRoY0t21IS9fw86ptBxsOMqxeDJuOhYxIvigiW1X1kwAi8ibgtcCT8zyJyGHAB4FnAXuBm0Vkh6p+w3N2GfARVb1KRM4C3g6cJyLPB6aBKeAI4EYR+ZyqPgTMAS9U1TtF5D8A/wV4TXiSu4E1asagqPsImzbfVtkEbaiz/c7zEENyBnC5iPwu8FjgTuC0AH+nAfeo6h4AEfkE8ELANyRPBt4Y/30D8Bnv+Y2qehA4KCK3AmcDnwQUODp2dwxwf4As65auLdrlEVo5umyEuyx7Fk1ejFUnIfrVr7t2ukZPQ6Kq3xGRvwXeAqwAb1HVAwFhHw/c5/3eC5yecHMrcA7R9NeLgKNE5DHx84tF5N3AGHAmDxugC4BrReTHwEPAM9IiF5ELgQsBJiYmKvWulpeXgbU9tLRndZIW/v79+xuJK4tt27YdIkOSkHzIclM0PaqaGk4y/LJlU7VMXXqqhNO0XhWJ0y+fptO0ZcsWINpyXHfYjl76lqVfvQipJ3WyvLzMwYMHc+Mrm5ayhKyRfB74DvBrwARwpYjcpKpv7uU15Zkmfr8Z+ICIvAa4Cfg2cFBVrxeRpwNfBr4L7AQOxn7eCDxPVXeJyJ8A7yYyLmsjUr0cuBxg06ZNWmU//OjoKMCaPfV5C10+ZXuYaXGm/R40WXKGuimSHhFJ9ZMMP0SmonKGsmHDhmDdaEqGouTJ6+RoOk39KrM8d1n65dOGEWNImkPSUichU1sfVFU35fQDEfl1otFJL/YCJ3i/J0hMQ6nq/cCLAURkHDhHVffH7y4FLo3ffQy4W0SOA56iqrviIK4G/jZAlkoMQmm6Mq/cz29Wmp4iafu3NXm0Va6mWW/phYcvucuj39OJIVNbn0n8PghsCwj7ZuCJInIS0UjjZcArfAcisgF4UFVXiIzTlfHzw4BjVfV7InIqcCpwfeztGBGZVNVFooX8OwNk6RxdmVduA11pTNZrY98VrM6VJ2Rq64c8PCX1C8AocEBVj8nzp6oHReQi4Dqi7b9XquodInIJcIuq7iBayH+7iCjR1NbrYu+jRLvFIFoHOTc2YIjIHwCfEpEV4PvA7xdIr9FS6hyBNdlQNz1SHEYj04Y0hXx7UZecbegw9DvukBHJUf5vEfkdwnZtoarXAtcmnr3N+/sa4JoUfz8hY3uxqn4a+HRI/MZwUaVytKFydxnLPyOPwl+2q+pnRGRLE8IYxWhD5W6DDEb3Mf3pNiFTWy/2fo4Amzh095WRgVWQMLoyP90VOYtgnYFwLI/SGQlw89vev+cAPyT6sNBokLm5OVNaI5g23YVehtnZ2c7K3m9mZ2dXv7tpCyFrJOf3QxBjeBmmHm/Td403nVddLosuyx5KV9OYaUhE5P3kTGGp6usbkchYl9jRE0aThHx7UZZk4z9Mu79CyRuR3NI3KTpIVwq5aTnbnn6jHrLKuQ31oA0yrHfyDMlH3bcbhtEV7Gpho40Mu77lGZKvEh3ljoi8X1X/Y39EGn7aoFR2K1536Mqochj1ZRjT1AR5hsQ/dPGZTQtiDC9WCcOpezG/Kw1hXWe2rQeaXO8pS54hsW9FjDV0pVEysulS2XVF3+qUb9A3TpYlz5A8SURuIxqZ/HL8N/FvVdVTG5euJdiWTKMIVp79pYv5PGw6kmdIfqVvUnSQNuxiCZFhWBQ1lGG6EXLQbNmyhdHR0dy8KpqPw9aAGhGZhkRVv9VPQbpOka9yrRKFM6wNz7Cma9ioq3yGvZwLH9poDAdNLW5aA1k/XcnLrsg5TLSlvpkhMYIZtLIaEUUajy6XWZdlX2+YITEaIa+xc++2bQu5aLN71NkAtqXHaRh55J21dTv5Z22tm11bdg5Uu7HGtl6aPIW36TLqii60Xb6i5I1Ifiv+311/+5fx/68ElhqTqOM0uUW411WhTcswCNr+9f8ws337djZs2NBI2IPM2/Verk3Qc9eWiDxTVf0v27eIyJeAS5oWzugeXb8Xo18McyM2DA11V9LQFvlC1kiOFJF/o6p/DyAivw4c2axY3aNfBdoVBR8mLM8NI58QQ/Ja4EoROYZozWQ/8PuNSmV0npCP2Np2XlAbSctHM2j1YB2E+gi5IfFrwFNE5GhAVHV/82IZ64GQL6dDsIbAcJguDIaehkREHgv8KfBLqvpcEXkysFlV/7xx6YxV3ImfW7duHbQo64aQ9Z5h7NXaiNEoSsjU1l8AHwbeGv9eBK4GzJAYRkWaMkRtN3CDvHhsGDeEDLq8QwzJBlX9pIi8BUBVD4rIzxuWy8hgvVxE1aa0tEGGLtGG/KqqP21IQ5cYCXDzIxF5DPHHiSLyDKIF956IyNkicpeI3CMiW1LeP15E5kXkNhH5gohMeO/eISJfj/+91HsuInKpiCyKyJ0i8voQWdrA7Oxsz2mSYeolGYaxPggZkbwJ2EF0J8mXgOOA3+3lSUQOAz4IPAvYC9wsIjtU9Rues8uAj6jqVSJyFvB24DwReT7RNb9TwBHAjSLyOVV9CHgNcALwJFVdEZFfDEyrUQN1jhSa/OAthDaNetrAessP27JfH0G7tkTk3wEnE11qdZeqLgeEfRpwj6ruARCRTwAvBHxD8mTgjfHfNwCf8Z7fqKoHgYMicitwNvBJYBZ4haquxPI9ECBLLawHhTCKYbrQf6weto+QXVv/F3inqn7Ie/ZZVf2tHG8AxwP3eb/3Aqcn3NwKnAO8F3gRcFQ8jXYrcLGIvBsYA87kYQP0y8BLReRFwHeB16vq3SlyXwhcCDAxMVHLDpTl5ch+lg2rl/9e7/fvz59RrCpfCHXGkZeekHi2bIlmS7dv317oXZE4ipCVnjxZ3MGVdZdZmbQl/eSVT0j+ViFP/pC0pbnpVX+K0isP3HtH2TJOiyeZln7U/TxCpraWgTNF5HTgD1X1Z0RGoheS8ix5COSbgQ+IyGuAm4BvAwdV9XoReTrwZSJjsRM4GPs5AviJqm4SkRcDVwK/cUhEqpcDlwNs2rRJ65hCGR0dBSg9HdPLf0j4ee9CzuGqStU8SFIlL/LcVPVflrKy1E0ZXUiTs0r5VCFPfhe32wqfNjLJkq9OeUPrc9W4Q9LSj7qfR4ghWVLVl4rIfwK+KCK/R86pwB57idYyHBPA/b4DVb0feDGAiIwD57gPHlX1UuDS+N3HADfq2At8Kv7700Rbk40aaNOUQRtkGFb6Vc6D1CfTn/4SYkgEQFX/u4h8DbgOeHSAv5uBJ4rISUQjjZcBr1gTsMgG4MF4veMtRKMLt1B/rKp+T0ROBU4Fro+9fQY4K3b774i+a+kEvZS7C8rfBRlD6XpaqjbUbeo4GN0mxJC8zf2hqvMi8hzg1b08xd+bXERkeA4DrlTVO0TkEuAWVd0BnAG8XUSUaGrLHVk/SjT6AXgIODdeeAfYDnxURN4IHAAuCEhDKayiGVXo11buuvTU9LwZ1kO+5l1s9SRV/SbwbRGZTrz+bEjgqnotcG3imW+YrgGuSfH3E6KdW2lh/gB4fkj8VVlcXDvYWQ8KEUpI49W0IQ45GDKEQcppFGcYv0zPowv6kzcieRPwB8C7Ut4p0fSSYRjrkDY0biEy9OOolbLvfbo++5F3sdUfxP+f2T9xjDYsUHZdqYcBKwOjS+RNbb04z6Oq/nX94hiG0TT9mhoaZiNYl6Eflum5vKmt3855p4AZkg4yDD3drqShafmqGoS255+xljasS2aRN7V1fj8FaSOTk5ODFqG1DFMj1Ia0VGkAfD9dMbJtxvKwOCHbf4kPUfxV4BHumape0pRQbcEUqRpl8q8rlbgrcnYNy9duEnLW1od4+LyrK4CXAF9tWK5WY8peH5aX6bQlP5Ll07Xy8uWscrWzbRHPJ2RE8uuqeqqI3Kaq/01E3oWtjzRGGxSqDTK0ja41oEY3GBZ9CjEkP47/XxKRXwK+B5zUnEjDQT8bniJxDYPiDkMajG7TBh1Mq/dtvmr3syJyLPBOYDfRjq3BHjVpGENGGxqmprFRXTXanG8hF1tti//8lIh8FniEO6HXqIZVrO5SpcwG+bV1E7RBj+uUwepjcUIW2w8jOtvqROdeRFDVdzcrmgHtqKR1ETIUH4Z0GuVpsvyrXO3cT73sYp0Pmdr6G+AnwO3ASrPiGHXS1nWattOmuee6KFs+TRv6YdKb9UyIIZlQ1VMbl6RD2Adg7cTKwiiD6U11QgzJ50Tk2ap6fW+nhqOfSjnMN90ZhpFOm+pjiCH5CvBpERkhur9dAFXVoxuVzMjEGvbytCHv1mu5rdd0rwdCDMm7gM3A7aoacle7Ech6q1hNpzd5EZnRH9qgx22QIYSQjkxX0uITYkjuBr5uRmQwVFGqfl7800Xlz6JsWqocwdE0TRvZJkd6bRhFhtAVOZsgxJB8B/iCiHwO+Kl7aNt/q7OeFa8J7LTmbNqaN23Q/TbI0HVCDMm98b9fiP8ZHutBCddDGruAdTy6z7BOv+YakvhjxHFV/ZM+yWMYjWKNcPcIaXzNyA6WXEOiqj8Xkel+CTNMDNsxGF2gn2tCxvBgOlGdkQA3CyKyQ0TOE5EXu3+NS2YYHWRxcXFo7uE26mdycrLnetXs7GzndChkjeTRREfHn+U9szvbhwTrhdXH9u3b2bp166DFSKVf98c3QV0bBdp8sVXXCTn9t/Td7SJyNvBe4DDgClXdnnj/eOBK4DjgQeBcVd0bv3sH0WGRANtU9eqE3/cD56vqeFn5Bo0ppVGEYV2oXU8Ma50POf13Ang/8EyikcjfA29wDX6Ov8OADwLPAvYCN4vIDlX9hufsMuAjqnqViJwFvB04L74jfhqYAo4AbhSRz6nqQ3HYm4BjiyXVSMMOduwObd3CG8og1g1N5/pDyNRjUzZIAAAgAElEQVTWh4GPAb8b/z43fvasHv5OA+5R1T0AIvIJ4IWAb0ieDLwx/vsG4DPe8xtV9SBwUERuBc4GPhkbqHcCrwBeFCC/YazShgX5LjVubZC1y9Ny64UQQ3Kcqn7Y+/0XIvJHAf6OB+7zfu8FTk+4uRU4h2j660XAUSLymPj5xSLybmAMOJOHDdBFwA5V/Y6IBIgxGJpUztDK3c+vjdvQ4DRJaPqGNf1G/wjRobadohBiSPaJyLnAx+PfLydafO9FWiufPGblzcAHROQ1wE3At4GDqnq9iDwd+DLwXWAn0cjkl4hGRmf0jFzkQuBCgImJCfbt2xcgcvvYsmULEC3k7t8fXUy5vLwM0DNNIe5Cw+rlr0xc+/fvX5O+tpCWlpD0ufIpE35IPhQpqzrydf/+/aX1I406wyoT58GDB0vHHyL7tm3bSodflCppaYIQQ/L7wAeA9xAZgi/Hz3qxFzjB+z0B3O87UNX7gRcDiMg4cI67xldVLwUujd99jOjMr6cCTwDuiUcjYyJyj6o+IRm5ql4OXA6wadMmLXszWi+a7omPjo4CrN7stmHDhkOehfot6ybEX9m4ysbfiyrlUkXOkHSUDT/EjUt3XflaZ/mkhdXP+nPZZZeVTscVV1xRp1iVOfzwwxkdHa293pQlZNfWPwEvKBH2zcATReQkopHGy4jWNVYRkQ3Ag6q6AryFaAeXW6g/VlW/JyKnAqcC18drJv/K838gzYi0lWGf/llPDKIsm15QNv00ypJpSETkbTn+VFW35QWsqgdF5CLgOqLtv1eq6h0icglwi6ruIJqieruIKNHU1uti76PAF+NRx0NE24IPBqbJKEAbGsK2YYu79ZA0TINI93r7ZmRQsueNSH6U8uxI4LXAY4BcQwKgqtcC1yaevc37+xrgmhR/PyHaudUr/M5+Q1InXVb8rrBe8zaZ7hBd69pX2V2hzfU805Co6rvc3yJyFPAG4HzgE0SXXRkDIlSRmlS4ZNhl42pDpWhDBe3adxBtlWu9sH379tasj0Dv038fDfwx8ErgKmBaVb/fD8HWK22YDmg7XWvEyshbtVfflbwxhoO8NZJ3Eu2ouhw4RVUP9E2qltK1BqwsbU1nkcZ1va5zlDFWbRqJtUGmIuTJW2QasOtThXkjkjcR3Yj4X4C3eh//CdFi+9ENy9YJiih8VyrHIOlaQ9IG6sory/PuM6gyzFsjCTli3qiJKr2OtjYAZhTWN4Ms9yq6Z/panJAPEg3DGBBta9SSR+LkuanCMHVC6prianNemCEpwaBOMV1eXi79he0wVcy609DGwxnbXk5tl68uhqneNIkZEmOVkOm1YapYw5SWIrR58TZkxGO0DzMkBRgGJe/KEDpNzjbINUx0PT/LdgTq6EAUqUch04BNTxU2jRkSo1V0zViv11FNEbqYR13Rv7ZghqQldKmShTKMaWoTbWigi8rQr3td2qR7bZKlKcyQDBFtaFjaSptO6207Lq/c/Rpto67RQpvqS9dPbjZDUoKq2/jKxlnlEpu8KaMmZO7C0eZtaECMdlOXjgy6oW8aMyTrjDZMGbS9xzvshJZhk9e51tGwVj0o1NZB6sMMiRFMP3tVIZW9ijxF/Oa5bXsPc9h7wknqSm/TazghdKnszJDURFaDNwhl6JICrlfWUxmlbeHu6migrq/Uhw0zJIbRJ+puYNrQUPky9KsBbUO628ag88QMScMsLi72dFP3NMsge3tt7Y21dUTYxvxyslxwwQUDlmQtXTbEbSrfJjBDUjNtbbDyGJSShyzmtuEbCaMabTSWRr2YIVmnDLJy19Xj7deOn0F3DgZVVk1e55pMS9c7RYOiLUbaDEkAdRdWWwq/KF2Tt80UmYYMmR7tFY9RD/1sA5o6CbwJzJDURFahT05OlgqvDUeQd9XgpdGG3TbDkI9Gb5o+DLKf4YRihqRjdLExalqp22bw2jpt1jR1p6VO41+2F9823WorZkhyqGOxNUQBB62kdX7YlxZG1tSMPwffhhFDkkGXS79oay84qTfrpTzqoMp0aBnMkJQkq9LUVZmSipDWo6qyvdQqZTvp6pbhPGyqtduyh9CoIRGRs4H3AocBV6jq9sT7xwNXAscBDwLnqure+N07gOfHTrep6tXx848Cm4Bl4KvAH6rqcpPpqJthVab1QL+OZfEpu85WpwxNh9flOlEmD6p2At0Brm3Jt8YMiYgcBnwQeBawF7hZRHao6jc8Z5cBH1HVq0TkLODtwHki8nxgGpgCjgBuFJHPqepDwEeBc2P/HwMuABrNzZDpmrqvqa2z8SgrS4ic/R5Cl6ULU4xpND1dNCg56mB2dpbFxcXG60o/aVP+FqHJEclpwD2qugdARD4BvBDwDcmTgTfGf98AfMZ7fqOqHgQOisitwNnAJ1X1WudZRL4KTDSYhlbTxsrdBprsdXfFcA6SuvK/TgOxuLiIqtYWXtvpt3Ft0pAcD9zn/d4LnJ5wcytwDtH014uAo0TkMfHzi0Xk3cAYcCZrDRAiMgqcB7yhEenpTm+7KlUq/NLSUk838/PzwKGjNv8Y+UEcJ95GBtk5yPsYssk48p73wjWYMzMzTE5O5m7DX17OnwFPk6HuNdBkfpYJP2QHWr/1p0lDIinPkl2CNwMfEJHXADcB3wYOqur1IvJ04MvAd4GdwMGE3/8fuElVv5gauciFwIUAExMTpS6FOumkkwBS/brejXuXVNI0P86Ne7dlyxYg2r2UxDWyzu22bdvYv3//mnBdeO4L8bRwknH2el6URz7ykanhpKUtmUf79+/PdZ8kmSdZfupKmwvfsW/fvkPK3cdPTxppciXj8OPK8lMlfWny33XXXQBs3Lhxzbv9+/f31Ovl5WX27NnDBRdcsFoOLk0hdSDvebK8k2582VSV5eXlzDxZXl7m4MGDuXnm8iHNTZou9EpTMv408vIhL6+SaQmpP03SpCHZC5zg/Z4A7vcdqOr9wIsBRGQcOEdV98fvLgUujd99DLjb+RORi4kW6P8wK3JVvRy4HGDTpk1a5qiHe++9F4CtW7cCa638ySefDLC6fXV0dHSN37T4nJukn7Tws/DDTcbpwnHMzc1l9lqSsoSS7DGJSGo4aeEn5T3mmGNWZS4jT1b+3XvvvSwtLbF169bcr4N9P3nhOzZs2HBIuSfJkz+tLJJxJMNJy5e8nmivdKWVl0tTMm4/fjc6TytnEWF0dDRIr7PKOaT8k+n28y4pQ5LR0VH27NmTqxNZupyMK+km5NsU5yY5IsnLh2SeOL95daupI2160aQhuRl4ooicRDTSeBnwCt+BiGwAHlTVFeAtRDu43EL9sar6PRE5FTgVuD5+dwHwHGAm9tcqiuzGKEu/vm9Jxpfmx00tNDUl089vHNq45lR1O3CVNG3ZsoV77703d759bm5uddF7dna2dhnyKHLMTD9Jprds+tukh71ozJCo6kERuQi4jmj775WqeoeIXALcoqo7gDOAt4uIEk1tvS72Pgp8Me4hPES0LdhNbX0I+BawM37/16p6SVPpyCK0kH0lqkPx/Z0qyfBcr7HIQlvRI8/bvEvGpWVycrLW9S3XWKY1lIP4aj/tmWvIe+HWq6rK4yiys6+I/tRlHMvWtSIdD0eRvKirDrXFiDb6HUm8w+raxLO3eX9fA1yT4u8nRDu30sJs5UeUg+w9VGk0fb9lKkbSmIUotvNzwQUX9OzxDgNFGr42HCNTpnFaWlpa1aW6Rgp5o5yizM3NHXLadDJP8vQwaQDa2qka1Ki6lY2yUQynzFnfttS1w6QIab34fiq337AVJavHH+o3OWIMiSMrr/u5czBNFpeOmZmZ0uH6O/uy8qXKlKPTNd9tVbLyvYiRTHPj6mqRkW1aBy+5AD/oaT4zJA2RnG5KI63w65pf7SWXH97S0tLAh8hZFSvET1pF69Xw5U1X1U2IMWjLdl+fPXv28OMf/7hnWYQamrTRb5vWAfJkCRl5uDwoM7WcpiODrpNFMEOSQ7L30NRxEH6j1ovFxcXMbzey5K2KH16vHqWLc2ZmhqWlJaampnou1teVn8medCiho4M0P3VRZk4+jSz/VUc2yfwsYhw3b96c6abIdOogGtYycYboUZl3dW2QaQIzJDUTMhIpG67bVz42NlZoYc+XK0mRheleC75NUCT8PLdNjfDKNHxF5teLzsEn0+f+90cMvRpHpw8bN25c3QLfK748OcfGxnpInU+VBeqQqaK0jsSBAwcYGRnpKZM/rRYiQ52dvDat05ghKUHeTpqi/oriK01TjWKZ3lJSlpmZGVZWVnIVvcmrXKuQZlirVNbkFumqsuStZWTFEXICQRI3teXH6UhOZbmRsh9/yCK279/Hj8/F5dyU2Xk2OzvLHXfcwdjYWNDIfWRkhLGxsaA404y1L2+Tna7kiL7uXWGhmCEZIGUVrNdWz5BeosPNb/sVOauClekBpTVgLtxt27b1zINB7MFPS18yvKxGIqRzkdzdVFSWMvijgl6NjXvuvvSG3uUQMrJdWFhgZWVl9TiTNNw02M6dO3sa3jSZiuwizBqZAKysrLC0tJQ5msr7fiotzl470Pz8yCqfkAX5QY1OzJDkUKQxCukJJCtaWsWrY0osRBZf2V2FSSpqmrzOMPSS07nze3VJ+fIosxHB/51V6QZxflqZ6Yy8hdsi4STdzszMcODAAcbHxyuFm0bWTjm/x37gwIE170LKI2tXonuWds5WcpRQtIF1I5IkeXmUjKNqoz7onVhFMENSgtC1hyycH9fY+j39pL88Zapr6OwP90PpNVUyNja22miUWfjOo2zjkEW/FiR9vXHTQMmRTF2NRpWFecfGjRsZHR3NnTLL67Tk6WXWdJDvt9fo13VsxsbGWFhYYGZmJmjaq8j0T1LP8/w2vZbjP19eXk49XsemtjpGmS2zWY1v3k4sh9s7vmHDhkPiLdvDzyLETdaQ3x9iLywsFA43jzIjm6y4qy5UVl0zWVhYCJY/OS2SV/7JhiTrd1Y8vtsQ0vS2yHqMczs1NQWwqjN+nehldNJI07VHPvKRhU89cHqe7PxVJW/6GMLLIPREg6YxQ1KANIV2iu8ULu97iDK7V3z/yaOj/emjPHl74aYkyjSOLi9cT9AfLfQa6YSc5ZQlL5TLz9B1I8jvUWc1RrOzsywsLKSmPS/uZCPpN6JuDSukwUjGkdbwjYyMpG7acLj0Jr8ETyOvDJIdiTxCtgqnpW1lZe1xe86/M0wh4STr9fz8fM/dl26UlpwVyKLKB7JZBt7XDfeuLkNXFDMkASQblDQFzlp8C2nM08ILIVmJ05Q52dBnhVO0p1ZHY97rfoisOLNoav0jbe2ljNH1dwC5ckmWWZGGIG+UmWzEi/TqnZurr756dVddssFMjnD8jxKdTiTT4nQsZKtwXp0o+5X9xo0bU6fp0sILHcHnlVfazrMkacY26zQBX6Y9e/akylZ1q3VZzJAEkNXryJtSSCNZ+ZLDeZ88Be7VKJQ9mC9pHJLh+Maml8KGbo3O+iCxSH6GTFulTfH4a1RJmdNkKSpfFfJ61FnkTdfllYFLS7L8t2zZsrpGUkc6i+xSW1hYYGRkZE2Dmux1+3qaVReKyJ23Eyup9y6vyzbceWtOvabNFxcXV7dlO8pM/9WJGZIKJM8Q8imyAO4WphcXF1OniODhq0KrLqznLepV6c2ETLukNb7JKZu8Of1Q+ULWnFx4fg81b9daVuUOGQWlbQ6osg6RtRMtLcxknvnTNslGp2pvNhl/SAcpy2CmyZI3/eUbIP//nTt3psYbSkgDnZy2ytOJ5JRZ2RG925o9qKmsJGZISuDvSEoqQl6lzpp7L9KL8JU0qzFK29qbRdFF5+R6SrLS5Cl2kdNVQ/wX2UCQtficRtp0U+jU4+LiYqFpymT+FWlY6tiuWzTOZP6V/aI9+aFfMi1Fv/oPWWNJunX1LrmlPS3uEHmSncAi9TCNpM4W+XasXzsRHWZIKlDXfKQ/d5y18Okbm6wviKt8X5CUJ5SkW7+HmTUV2IYjHaDajpdeH/P1mh9PNnx5PedeepYc6czNza2GW6Ysk73rNJlCjqPp1YtvkrQOnrtmuApp36fs3r07OD29RrD+jr6ssvOntlw6q3z1XwdmSCqS1aMqQsh5PZOTkywvL685+yh02qxXY5I8XiMvXBde3u60MmTt128D/tRjlW9i0nYHJRu8rONPQkk2/HV/n1LEADTdqKWN8J1BLjLKd2VQtoOT9fGiw+9AldHrIuuug8IMSQHcvOvU1FTm1sYQZay6eOvOsAohazF7kPjbS7O2RYZMDSbnm7OmvtL8+lsnk52BpDxp38SEfNPRFFlGwR1BAvmL9b1kdPnhb/8NGWXUoWdpYSQNUpmNCK4jlvWuqFz+yDOZJ0Wm5UJGdlnfxIR0QPuFGZIAknOfaQVYZr6+CO6DxK1btxban58k76vwkO8nsiiyNuRTZutxkjwDEkLI6QRZcqbFlZx6bKpnnszzsbGxoBFDcmE6q1Hbvn07L33pS2s1jnXsLqpqsKps4Q4hTd/dxo68L/md314j2TyjmBZePzBDUpK6Cyjr47W071JCLm2qgyLfxFStlMnGrExFCPHru0kugOaNbEJ7fkWNYpkPJMs07Gl5UaZnXxdNr5NVMd5F8zcrLXnTlFXrqPuOZNDbfh1mSGqmqoLUPacdGl8T4bVtSi1Jr9FQP+TPiqPINxdpjUmWf3+XXhldq8sAJOVrOq/diD4kzpBGuaye9/pWKW8XpV/eyVMHBr2BxQxJQ9Q1tAxdPygqwyAa+bSpk61btwa5HTRlGt2sNAxqZ00ZBt1AhdA2XUlSVL7QdU3/DLZB54EZkgAGsaBVZshaRJn6PYeahv/ltCNErn70Xqv4bUpfQnQi9PbMtHDrdttVumA824YZkgLU1ZOsOkKoUpmbXmAcBGWNYnLbcwghcTSVF70WiUP1s407+fpJG9NeZudYXW7rwAxJy6m74Q+ZomljRStC2ca+a+ke5BqOMVjypoUHgRkSY93Ttcaya/J2Dcvf4jRqSETkbOC9wGHAFaq6PfH+8cCVwHHAg8C5qro3fvcO4Pmx022qenX8/CTgE8Cjgd3Aear6sybTMUhMqY06MX1qJ2XKpU1l2ZghEZHDgA8CzwL2AjeLyA5V/Ybn7DLgI6p6lYicBbwdOE9Eng9MA1PAEcCNIvI5VX0IeAfwHlX9hIh8CHgt0J4cjWlDIQ9ywTqE7du3r953YRhGd2lyRHIacI+q7gEQkU8ALwR8Q/Jk4I3x3zcAn/Ge36iqB4GDInIrcLaI/BVwFvCK2N1VwH+lYUPShkZ3vWB5bRjdo0lDcjxwn/d7L3B6ws2twDlE018vAo4SkcfEzy8WkXcDY8CZRAboMcAPYgPjwjy+sRQYnWC9GJ/1kk6jezRpSCTlWfIc5zcDHxCR1wA3Ad8GDqrq9SLydODLwHeBncDBwDCjyEUuBC4EmJiYOOSr1i6yf//+QYtQK5aedmPpaS9tS0uThmQvcIL3ewK433egqvcDLwYQkXHgHFXdH7+7FLg0fvcx4G5gH3CsiBwej0oOCdML+3LgcoBNmzbpsMzFD0s6HJaedmPpaS9tSstIg2HfDDxRRE4SkV8AXgbs8B2IyAYRcTK8hWgHFyJyWDzFhYicCpwKXK/RzTQ3AC+J/bwa+F8NpsEwDMPoQWOGJB4xXARcB9wJfFJV7xCRS0TkBbGzM4C7RGQReCzxCAQYBb4oIt8gGlWc662L/Gfgj0XkHqI1kz9vKg2GYRhGbxr9jkRVrwWuTTx7m/f3NcA1Kf5+QrRzKy3MPUQ7wgzDMIwW0OTUlmEYhrEOMENiGIZhVMIMiWEYhlEJMySGYRhGJSTaUTvciMh3gW8NWo4a2ED0Lc2wYOlpN5ae9tKvtDxeVY/r5WhdGJJhQURuUdVNg5ajLiw97cbS017alhab2jIMwzAqYYbEMAzDqIQZkm5x+aAFqBlLT7ux9LSXVqXF1kgMwzCMStiIxDAMw6iEGZIWISIniMgNInKniNwhIm+Inz9aRD4vInfH/z8qfi4i8j4RuUdEbhOR6cGmYC0i8ggR+aqI3Bqn57/Fz08SkV1xeq6OT4dGRI6If98Tvz9xkPKnEZ9M/Q8i8tn4d5fT8o8icruILIjILfGzTuoagIgcKyLXiMg34zq0uavpEZGT43Jx/x4SkT9qa3rMkLSLg8CbVPVXgGcArxORJwNbgHlVfSIwH/8GeC7wxPjfhbTv7vqfAmep6lOAKaLrkp8BvAN4T5ye7wOvjd2/Fvi+qj4BeE/srm28geg0a0eX0wJwpqpOeVtJu6prEN20+req+iTgKUTl1Mn0qOpdcblMAU8DloBP09b0qKr9a+k/ortWngXcBTwufvY44K747/8BvNxzv+qubf+IrkzeTXTd8j7g8Pj5ZuC6+O/rgM3x34fH7mTQsntpmCCqvGcBnyW6sbOTaYnl+kdgQ+JZJ3UNOBq4N5nHXU1PIg3PBr7U5vTYiKSlxFMhTwV2AY9V1e8AxP//YuzseOA+z1vr7rCPp4IWgAeAzwP/F/iBPny/jC/zanri9/uJ7pxpC38G/CdgJf79GLqbFoiuqb5eRL4m0dXU0F1d20h0LfeH46nHK0TkSLqbHp+XAR+P/25lesyQtBCJrh3+FPBHqvpQntOUZ63ahqeqP9doeD5BdI/Mr6Q5i/9vbXpE5LeAB1T1a/7jFKetT4vHM1V1mmha5HUi8m9z3LY9PYcD08Ccqj4V+BEPT/uk0fb0ABCvub0A+KteTlOe9S09ZkhahoiMEhmRj6rqX8eP/0VEHhe/fxxR7x6iXscJnvfMO+wHjar+APgC0drPsSLiLlXzZV5NT/z+GODB/kqayTOBF4jIPwKfIJre+jO6mRYAVPX++P8HiObfT6O7urYX2Kuqu+Lf1xAZlq6mx/FcYLeq/kv8u5XpMUPSIkREiK4OvlNV3+292kF0Pz2svad+B/CqeMfGM4D9btjbBkTkOBE5Nv77kcBvEi2A3gC8JHaWTI9L50uA/6PxhO+gUdW3qOqEqp5INNXwf1T1lXQwLQAicqSIHOX+JpqH/zod1TVV/WfgPhE5OX40A3yDjqbH4+U8PK0FbU3PoBeS7N+aRbV/QzQcvQ1YiP89j2hufR64O/7/0bF7AT5ItO5wO7Bp0GlIpOdU4B/i9HwdeFv8fCPwVeAeoiH7EfHzR8S/74nfbxx0GjLSdQbw2S6nJZb71vjfHcBb4+ed1LVYxinglljfPgM8quPpGQO+BxzjPWtleuzLdsMwDKMSNrVlGIZhVMIMiWEYhlEJMySGYRhGJcyQGIZhGJUwQ2IYhmFUwgyJ0WpE5DHeCaj/LCLf9n5/uc+y/LGIfCM+XXVeRB6feH90LN8HvGdfEJG7PJl/MeHnJSKiIrIp/j0qIldJdCrvnSLylvh56knK8bs/j5/fJtHpt+MZ8v+OiLwt492B8jmTGt7fuZNpjeHHtv8anUFE/itwQFUvG1D8ZwK7VHVJRGaBM1T1pd779wLHAQ+q6kXxsy8Ab1bVW1LCOwr438AvABep6i0i8grgBar6MhEZI/qo7gzgW8CRqnogPv3g74E3qOpXRORojY/SEZF3Ex3lsj0lvi/HYe9LeXdAVVMNUBlE5NXAhKpeWleYRnuxEYnRWVwvWkTOEJEbReSTIrIoIttF5JVxD/52Efnl2N1xIvIpEbk5/vfMIvGp6g2quhT//ArRMRROlqcBjwWuLxDkNuC/Az/xowGOjI9VeSTwM+AhjXCjhtH4n8ZyOSMisZ9DeociMgn81BkRie5R2RnnwzbP3Xg82tod590L4+fbJL4fJ/59qYi8XkQeJyI3xaOtr4vIb8ROdhB9lW2sA8yQGMPCU4juCjkFOA+YVNXTgCuA/xi7eS/R3SFPB86J35XltcDnAERkBHgX8CcZbj8cN7Rb48YeEXkqcIKqfjbh9hqiAwe/A/wTcJmqPhj7WXOSsj58rhQi8mHgn4EnAe9PkeGZRMf4O95LdMDh02N/jp8AL9LoMMczgXfFMv858dEccXpfBnwUeAXR0flTRGWwAKCq3weOEJG2nXhsNIAZEmNYuFlVv6OqPyU6JsKNDG4HToz//k3gA3FjvAM4Op5eKoSInAtsAt4ZP/oPwLWqel+K81eq6inAb8T/zosb4vcAb0pxfxrwc+CXgJOAN4nIRjj0JGUR+TXnSVXPj/3cCbz0kFCjuyu+6/1+Jg+f4fSXfvKAPxWR24C/IzqK/LGq+o/A92ID+GzgH1T1e8DNwPnxtOMpqvpDL6wHYpmMIccMiTEs/NT7e8X7vUJ0xDhE+r5Z45vnVPX4RMOHiLjRw7VpkYjIbwJvJVprcHFsBi6S6GTgy4gOz9sOoKrfjv//IfAxIkNxFPBrwBdiP88AdsQL7q8guuVvWaNTeb9EZLRW0YdPUj478fznwNVEo60kPyY6/2uNlxR3ryRa53labLT+xfN3BfAa4HzgyjjOm4B/C3wb+EsReZUX1iPieI0hxwyJsZ64HrjI/RCRqaQDVT0/NjLPS76Le+P/g8iIPOD5eaWq/muNTgZ+M/ARVd0iIoeLyIbY7yjwW8DXVXW/qm5Q1RNjP1+Jw7yFaDrrLIk4ksjIfFPST1L+ZuzuCfFzAX4b+GZK2u8EnuD9/hLR9BRExsNxDNFi/XK8ucDfmfZpIuP1dKIbIJFo59oDqvo/iaa/pj1Z/hXRLYzGkHN4byeGMTS8HvhgPG1zOHAT8O8L+H8nMA78VbzU8U+q+oIc90cA18VG5DCiqaL/2SOODwIfJjotWYAPq+ptInIqcJWIHEbUAfykqn42nia7SkSOjt3fCsymhHsT8XqHRls13wB8LF5A/5Tn7qPA34jILUTrHatGSVV/JiI3EN0K+fP48RnAn4jIMnAAcCOSpwFf0YdvjzSGGNv+axjrhHh78t+o6t+V9D9CtB7tFQIAAABmSURBVGD/u6p6d0BcO1R1vkxcRrewqS3DWD/8KdEdF4URkScT3a0y38uIxHzdjMj6wUYkhmEYRiVsRGIYhmFUwgyJYRiGUQkzJIZhGEYlzJAYhmEYlTBDYhiGYVTCDIlhGIZRif8HUTLBWNnVKK0AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.plot();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now all of the light curves stitch together with no visible join. We can use the `fold` and `bin` methods to form a planet transit out of all of the data for Kepler-8b." ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3X+cHXV97/HXOyFikgW5sviL9ZJQiV5iIl1XKPVaIZGKoqJSARUUlOLditJbsY33WnsLpUQFvVa4aSlCwYcoKZbeaENB9/LL8kNCGgg/zJIGfwSouKLB/KAk5nP/mJnN7Mk5Z2d3ds7uOft+Ph772HNmvjPn+z3nzPnM98d8RxGBmZnZeM2Y7AyYmVl7cyAxM7NSHEjMzKwUBxIzMyvFgcTMzEpxIDEzs1IcSMzMrBQHEjMzK8WBxMzMStlnsjPQCt3d3TFv3rzJzsaY7Nq1i332mRYfzzCXeXpwmdvHfffdNxQRB42Wrv1KNg7z5s1jzZo1k52NMRkaGqK7u3uys9FSLvP04DK3D0k/KpLOTVtmZlaKA4mZmZXiQGJmZqU4kJiZWSkOJGZmVooDiZmZleJAYmZmpTiQmJlZKQ4k1jL9/f309/dPdjbMbII5kEywsf5YTtSPa7P9TFaezNpB2eOjyPbLli2b0GNqqh2j02KKlDL6+/sZHBwEYMGCBaxYsWJ4OTDieT4dwNKlSwEYGBjYa5v8l2BwcHA47fbt25kzZw7XXXdd3dcZTbafTPbag4OD9Pf3D7/24OBg3fJkaeuVofaLWzRPjdLXy4dZM2M9HkbbV/6YHW3fzX64s+OrVna8bdq0CUkjjsH8a9X+NuSPi3r5yvJdm6/8sdusLBNNEVHdzqUrgbcBT0XEq+usF/Al4K3AduCMiFibrvsg8Ok06V9ExNXp8tcCfwfMBlYD58Yohejr64vxzLXV39/PunXrANi9ezcAXV1dAGzdunU43YwZM0akqZVfP2PGDI444ogRX4T8/rK0Cxcu5LHHHmPr1q3MmDGDOXPmDD++6667RuQt23eWt/w+s9cDWLt27fBr5POS/7IODg4OB7MFCxYMb5NtN2fOnBHrs21gT9DK3rtsX8CI18l/2fMHzFSej2gif8DypnKZq1K0zPXe86VLlw5/p/LfwSKfS+3+li5dOuKYyu87Oz7mzJmz136y72t++9pjPzu21q1bx+7du1m0aBEPPfTQXsdp7b5q1f7eZL8PwPBxnf8tqX1vyn5fJd0XEX2jpqs4kPwOsBW4pkEgeSvwMZJAchTwpYg4StILgTVAHxDAfcBrI+IXkr4PnAvcTRJI/ioibmyWj/EGkuyL1ShATKSurq66X6SJ0NvbCzAiKEy0RgdTflmzALZixQrOOussHnvssVIHQJFa03iCwngDyWjbFflRHUtNcKz5rDJ9oxpnVuZsff5HMvsBzP841q4vekyOdoKXpRnP8V1k35Mte7/yJ3hjVTSQVNq0FRG3S5rXJMmJJEEmgLslHSDppcAxwHci4mkASd8Bjpd0K7B/RNyVLr8GeCfQNJCMV+0ZeZWqCiJQbQDJ1Dugapfln2fNedkZ1NFHHz28fu3atRx99NHDabOzwnoBpvaHrbZ5MZ8u+1HL1+Sa7auReunyNTqg1MFbu8/8frPljc7SYc+Zar381ja31msiadb0km2Tz0Oj5t9GecieZ2frmd27d7N169aG39dsfVFV/shXte/xBrZ6xvp+lVFpjQQgDSTfblAj+TawPCK+lz4fAP6EJJA8PyL+Il3+p8AO4NY0/ZvS5W8A/iQi3tYsD2Watup9qbMPu6uri4GBAY466qjh5fkzpqJfiny6Ztt0dXW1rIaU5WUsZ4BV5yWTb14D9mr+y6fPn+Fm6WBkU2W9Jrz855E1ReabDbLXzfaRf4+y7wWUq5E0agbMHgN1P5tG36He3t4RTZf1mojyec0HiEzte1DvfcjW55uJOkH2/jX7cc6+L/nvU/59gda/H/fcc8+4t50SNZICVGdZjGP53juWzgbOBujp6WFoaGjMmdu5cydz585lx44dzJ49G4AdO3YAyRdm/vz5DA0NjUgzf/58ADZt2jS8n0MPPZRNmzYNpzn00ENZvnw5y5YtY9OmTcPPIRndkbWlzp07l5UrVw6nmz9//vB+MrNnz2bHjh3DgWvhwoXDr71y5UoATj75ZLZt2zac7/y2ANu2bdtr2+w1IoKFCxcO56lIf9DChQsBhreZKNm+sgM5f0DXnn3l0+YP6ux9zd6zemfAtWfJa9euZdGiRaxfv36v1230PDu5yCxZsmT488jbsmVLw/Lu3LlzOH9Znmv7uepptDxfznyea9+DJUuWDH/+9fY1d+7cuj+mEcHs2bPZtm3bhJ4J1/YxjvadatSkmlm0aBGbNm0aPiayNPnjIZ8uO243bNjAjh07mDt37l7H9MqVKzn55JMBuOCCC0Yc2wDr16/f69jI+9a3vsUJJ5ww4nUzO3bsYOHChSN+M/K/AcCI9fnjbu7cueP67Rurya6R/A1wa0R8PX2+gaQ2cgxwTER8JJ8u/bslIl6VLn9vPl0j462R1Ko9Q6sdVQHNR000OjutN7pr586dXHHFFU3z0GhUR6PRYfVeu1m+8uvy29cbndJoFFht/urV8vJ9J51y9lpPb2/vXu9xsxpJvrlvrNrlvcwPHqlX6wNG9JkMDAzsNZIw++7V1ray/oFmoxTrHZ+ZRsd2o5FUzY61s846i1mzZjVtNs2aKLP+w2YaHbMTPShkSnS2pxmZR+NAcgJwDns62/8qIo5MO9vvA3rTpGtJOtuflnQvSQf9PSSd7V+OiNXN8jBRgaQq9QJBo0DSKfIHXdZnMXv2bCSVakobzw9obbNY0f01GiDRbBTPWAJJo5E8jdQ2ydU2wY2mt7d3QvvTamuv2aCPDRs28MpXvnI4XZGTrak2VLzIyVe9AQbtZkoEEklfJ6lddAM/Bf4MmAUQEX+dDv+9FDieZPjvmRGxJt32Q8D/SHd1YURclS7vY8/w3xuBj1U1/HcytesXb6zqjdqCkcMi8+3Lmdo+jHXr1g239WeBqbbTubY22WjYcrav7Gw2/zz7cc/6Vmr7LfJ5zX7Im51ljhZI8n1Btf0jMLJ/IkuTdWJneaw37DpfpnzZa0dR5cuRf618v1K+RplP3+g6pU4/SaqnXY/nKRFIpgoHkvYw2o8qjLwwMh8IJnLYb7MLScdqtHw1KnPRM/BGF4yOtm2zazSaNa3Uvje1+2pU3nwaf7fbhwNJjgNJe+jkMjcKRs0CSb30VeZpIl6zyD46+XNupF3L3C6jtsymhbLTybTCRLzmVOi/sNZzIDGbglpRI/GPvk0Uz/5rZmaluEZiNgW5tmDtxDUSsxZYunTpXlP8m3UK10jMWqB2EsnRtKKPxGyiOJCYtYADgnUyBxKzKciBx9qJ+0jMzKwUBxIzMyvFgcTMzEpxIDEzs1IcSMymoP7+/ro3EDObihxIzMysFA//NZuCPPzX2olrJGZmVooDiVkLjLXPw30k1k4cSMzMrBT3kZi1QDvcIdFsvFwjMTOzUhxIzFrE/R7WqRxIzMysFPeRmLWI+z2sU7lGYmZmpTiQmJlZKQ4kZlOQO+atnTiQmJlZKe5sN5uC3DFv7cQ1EjMzK8WBxMzMSnEgMWsRd6Bbp6o0kEg6XtIGSRslLauz/hBJA5IekHSrpJ7cus9KejD9OyW3fKmktZLWSfqepFdUWQYzM2uuss52STOBy4DjgM3AvZJWRcTDuWQXA9dExNWSlgAXAadLOgHoBY4A9gVuk3RjRDwDrABOjIhHJP0B8GngjKrKYTZR3IFunarKGsmRwMaI2BQRzwHfAE6sSXM4MJA+viW3/nDgtojYFRHbgPuB49N1AeyfPn4B8ERF+TczswKqDCQHAz/JPd+cLsu7HzgpffwuYD9JB6bL3yJpjqRu4Fjg5Wm6s4DVkjYDpwPLK8q/mZkVUOV1JKqzLGqenwdcKukM4HbgcWBXRNws6XXAncDPgLuAXek2/x14a0TcI+mTwBdIgsvIF5fOBs4G6OnpYWhoqHyJWmjLli2TnYWWc5mnB5e581QZSDazpxYB0ENNM1REPAG8G0BSF3BSRGxJ110IXJiuuxZ4VNJBwGsi4p50F9cB/1zvxSPicuBygL6+vuju7p6gYrVOO+a5LJc5kY3u6tR+FX/OnaXKpq17gcMkzZf0POBUYFU+gaRuSVkePgVcmS6fmTZxIWkxsBi4GfgF8AJJC9JtjgMeqbAMZmY2ispqJBGxS9I5wE3ATODKiHhI0vnAmohYBRwDXCQpSJq2PppuPgu4QxLAM8BpEbELQNLvA9+UtJsksHyoqjKYTaSitYxOr41Y56l0rq2IWA2srln2mdzj64Hr62z3LMnIrXr7vAG4YWJzamZm4+VJG81apGgNwzURazeeIsXMzEpxIDEzs1IcSMzMrBQHEjMzK8WBxMzMSnEgMTOzUkYNJJJeVGfZK6vJjpmZtZsiNZI7JJ2cPZH0CXxBoJmZpYpckHgMcLmk9wAvJpnb6sgqM2XWiTz1iXWqUWskEfEkyQy7RwPzSO5ouLXifJmZWZsYtUYi6TvAk8CrSaaCv1LS7RFxXtWZM+skrolYpyrSR3JZRHwgIn4ZEQ8Cvw109l1azMyssFFrJBHxjzXPdwEXVJYjMzNrK0Watn7FnlvkPo/kXiFbI+IFVWbMzMzaQ5EayX7555LeiUdtmZlZasxXtqdNXUsqyIuZmbWhIk1b7849nQH0saepy8wmmK83sXZT5ILEt+ce7wJ+CJxYSW7MzKztFOkjObMVGTGzhGsi1m4aBhJJX6ZJE1ZEfLySHJmZWVtpViNZ07JcmJlZ22oWSL6WXnxoZi3kznZrN82G/34/e5A2c5mZme2lWY1EucevrzojZpZwTcTaTbMaia8VMZtA/f39w81WZp2kWY3kVZIeIKmZ/Eb6mPR5RMTiynNn1kEGBwcnOwtmlWgWSP5Ly3JhNg0sWLBgsrNgVomGgSQiftTKjJh1Ovd9WKca86SNZmZmeQ4kZmZWigOJmZmV0jCQSFov6YFGf0V2Lul4SRskbZS0rM76QyQNpPu8VVJPbt1nJT2Y/p2SWy5JF0oalPSIJM/5ZWY2iZqN2npb+v+j6f+vpv/fD2wfbceSZgKXAccBm4F7Ja2KiIdzyS4GromIqyUtAS4CTpd0AtALHAHsC9wm6caIeAY4A3g58KqI2C3pRQXKaWZmFWlYI4mIH6Ujt14fEX8cEevTv2XAmwvs+0hgY0RsiojngG+w931MDgcG0se35NYfDtwWEbsiYhtwP3B8uq4fOD8idqf5fKpAXszMrCJF+kjmSvqv2RNJvw3MLbDdwcBPcs83p8vy7gdOSh+/C9hP0oHp8rdImiOpGziWpBYC8BvAKZLWSLpR0mEF8mJmZhUpcofEDwNXSnoBybQpW4APFdhOdZbVTrtyHnCppDOA24HHgV0RcbOk1wF3Aj8D7iK5OyMkTV3PRkRfehvgK4E37PXi0tnA2QA9PT0MDQ0VyPLUsWXLlsnOQsu5zIlly5LuxOXLl7c6Oy3hz7nzFLlD4n3AayTtDygiir4jm9lTiwDoAZ6o2fcTwLsBJHUBJ2X7j4gLgQvTddcCj+b2+8308Q3AVQ3yfTlwOUBfX190d3cXzPbU0Y55LstlhlmzZtVd3kk6uWyNdHKZRw0kkl4M/CXwsoh4i6TDgaMj4iujbHovcJik+SQ1jVOB99Xsuxt4Ou3v+BRJ7SLrqD8gIn4uaTGwGLg53ewfgSVp2jcCnsDIOoqvgLd2U6SP5O+Am4CXpc8HgT8cbaP0pljnpNs+AqyMiIcknS/pHWmyY4ANkgaBF5PWQIBZwB2SHiapVZyWu8nWcuAkSetJRnmdVaAMZmZWkSJ9JN0RsVLSpyAJEJJ+XWTnEbEaWF2z7DO5x9cD19fZ7lmSkVv19vlL4IQir29mZtUrUiPZlo6kCgBJv0XS4W5mZlaoRvIJYBXJPUn+BTgIeE+luTIzs7ZRaNSWpDcCryQZ0rshInZWnjOzDpPdHdGd6dZpRm3akvRvwFkR8VBEPBgROyV9uwV5MzOzNlCkaWsncKyko4CPpNOd1F6hbmajcE3EOlWRzvbtEXEKyRDeOyQdwt5XqJuZ2TRVpEYigIj4nKT7SK4LeWGluTIzs7ZRJJDkr/sYkPRm4IPVZcnMzNpJw0Ai6VUR8QPgcUm9Navd2W5WEY/usnbTrEbyCeD3gUvqrAuS+a7MzGyaaxhIIuL30//Hti47ZuaaiLWbZk1b7262YUT8w8Rnx8zM2k2zpq23N1kXgAOJmZk1bdo6s5UZMTOz9lRk+C+STgAWAs/PlkXE+VVlyszM2keRubb+GjgF+BjJxYnvAQ6pOF9mHam/v394eK9ZpyhSI/ntiFgs6YGI+HNJl+D+EbNxKTIiy9eRWLspMtfWjvT/dkkvI5nEcX51WTIzs3ZSpEbybUkHAJ8H1pKM2Lqi0lyZTWOuiVi7KXJjqwvSh99M70Py/IjwrXbNzAwoEEgkzQROAOZl6SUREV+oNmtmZtYOijRtfQt4FlgP7K42O2Zm1m6KBJKeiFhceU7MzKwtFRm1daOk3608J2Zm1paK1EjuBm6QNINk6K+AiIj9K82ZWQfyNSLWiYoEkkuAo4H1EeF7tZuZ2QhFAsmjwIMOImbluSZinahIIHkSuFXSjcB/ZAs9/NfMzKBYIHks/Xte+mdmZjasaSBJL0bsiohPtig/ZtOeO+St3TQd/hsRvwZ6W5QXMzNrQ0WattZJWgX8PbAtW+h7tptVwzURazdFLkh8IfBzYAnJfdzfDrytyM4lHS9pg6SNkpbVWX+IpAFJD0i6VVJPbt1nJT2Y/p1SZ9svS9paJB9mZladIrP/juve7Wn/ymXAccBm4F5JqyLi4Vyyi4FrIuJqSUuAi4DT01v79gJHAPsCt0m6MSKeSffdBxwwnnyZmdnEKnKr3R5JN0h6StJPJX0zX3No4khgY0RsiojngG8AJ9akORwYSB/fklt/OHBbROyKiG3A/cDxaX5mktwb5Y8L5MHMzCpWpGnrKmAV8DLgYJLZgK8qsN3BwE9yzzeny/LuB05KH78L2E/Sgenyt0iaI6kbOBZ4eZruHGBVRDxZIA9mZlaxIp3tB0VEPnD8naQ/LLCd6iyrvTr+POBSSWcAtwOPA7si4mZJrwPuBH4G3AXsSm/1+x7gmFFfXDobOBugp6eHoaGhAlmeOrZsmX73DnOZpweXufMUCSRDkk4Dvp4+fy9J5/toNrOnFgHQAzyRTxARTwDvBpDUBZyU3X0xIi4ELkzXXUsyVctvAq8ANkoCmCNpY0S8ovbFI+Jy4HKAvr6+6O7uLpDlqaUd81xWp5e53jUinV7melzmzlIkkHwIuBT4IkmN4s502WjuBQ6TNJ+kpnEq8L58grTZ6umI2A18CrgyXT4TOCAifi5pMbAYuDkidgEvyW2/tV4QMTOz1ikyauvHwDvGuuOI2CXpHOAmYCZwZUQ8JOl8YE1ErCJporpIUpA0bX003XwWcEda63gGOC0NImYdzVe1WztqGEgkfabJdhERF4y284hYDayuWfaZ3OPrgevrbPcsycit0fbfNVoaMzOrVrMaybY6y+YCHwYOBEYNJGY20mg1DddErB01DCQRcUn2WNJ+wLnAmSTXg1zSaDszM5teRpv994XAHwHvB64GeiPiF63ImJmZtYdmfSSfJxmaezmwKCI8r5WZme2l2ZXtnyC5mv3TwBOSnkn/fiXpmdZkz2x66e/vHx65ZdYumvWRFJk+xczMprkiFySaWYt41Ja1I9c6zMysFAcSMzMrxYHEzMxKcSAxM7NSHEjMWsjDe60TOZCYmVkpHv5r1kIe3mudyDUSMzMrxYHEbApxH4q1IwcSMzMrxYHEzMxKcSAxM7NSPGrLbArxqC5rR66RmJlZKQ4kZi3kUVnWidy0ZdZCbrqyTuQaiZmZleJAYtZCbtqyTuRAYmZmpbiPxKyF3Edincg1EjMzK8WBxMzMSnEgMTOzUhxIzKYQj+qyduTOdrMpZHBwcLKzYDZmDiRmLZbVOOqN4FqwYEGrs2NWWqWBRNLxwJeAmcAVEbG8Zv0hwJXAQcDTwGkRsTld91nghDTpBRFxXbr8a0AfsBP4PvCRiNhZZTnMWsXDg60dVRZIJM0ELgOOAzYD90paFREP55JdDFwTEVdLWgJcBJwu6QSgFzgC2Be4TdKNEfEM8DXgtHT7a4GzAB991jYcLKzTVNnZfiSwMSI2RcRzwDeAE2vSHA4MpI9vya0/HLgtInZFxDbgfuB4gIhYHSmSGklPhWUwM7NRVNm0dTDwk9zzzcBRNWnuB04iaf56F7CfpAPT5X8m6QvAHOBYIF+TQdIs4HTg3HovLuls4GyAnp4ehoaGypanpbZs2TLZWWi56VLmZcuWAbB8+fK9ypxf16mmy+ec1+llrjKQqM6yqHl+HnCppDOA24HHgV0RcbOk1wF3Aj8D7gJ21Wz7f4DbI+KOei8eEZcDlwP09fVFd3f3eMsxadoxz2VNhzI/9thjwJ6y5ss8a9asvZZ1ok4vXz2dXOYqA8lm4OW55z3AE/kEEfEE8G4ASV3ASRGxJV13IXBhuu5a4NFsO0l/RtJB/5EK829WiUYjs5qN5jKbyqoMJPcCh0maT1LTOBV4Xz6BpG7g6YjYDXyKZARX1lF/QET8XNJiYDFwc7ruLODNwNJ0O7O24kBhnaayQBIRuySdA9xEMvz3yoh4SNL5wJqIWAUcA1wkKUiatj6abj4LuEMSwDMkw4Kzpq2/Bn4E3JWu/4eIOL+qcphNtEY1DwcYa1eVXkcSEauB1TXLPpN7fD1wfZ3tniUZuVVvn76I0sxsCvFcW2ZmVooDiZmZleJmIrMWa9QX4lFb1q4cSMxarFHA8My/1q4cSMxarFGNwzP/WrtyIDFrMddIrNO4s93MzEpxjcSsxdy0ZZ3GNRKzSeB7s1sncSAxmwI89NfamZu2zCZBbce6O9qtnTmQmE2C7du3j3ju/hFrZw4kZpNg927fAcE6h/tIzCZBb28vXV1dw7fWHRwcdPOWtS0HErNJsGLFiuHmrP7+/r2auszaiZu2zCbJ4OAgEcG2bdsAGBgYmOQcmY2PayRmkygLIoCvK7G25RqJ2STr6uryqC1ra66RmE2SrF9k69atgC9GtPblQGI2SebMmTP82CO2rJ05kJhNkgULFjB37lx6e3vd0W5tzX0kZpNkxYoVDA0N0d3dPdlZMSvFNRIzMyvFgcTMzEpxIDEzs1IcSMzMrBQHEjMzK8WBxMzMSnEgMTOzUhxIzMysFAcSMzMrRREx2XmonKSfAT+a7HyMUTcwNNmZaDGXeXpwmdvHIRFx0GiJpkUgaUeS1kRE32Tno5Vc5unBZe48btoyM7NSHEjMzKwUB5Kp6/LJzsAkcJmnB5e5w7iPxMzMSnGNxMzMSnEgmSIkvVDSdyQ9mv7/T03S7i/pcUmXtjKPE61ImSUdIekuSQ9JekDSKZOR17IkHS9pg6SNkpbVWb+vpOvS9fdImtf6XE6sAmX+I0kPp5/rgKRDJiOfE2m0MufS/Z6kkNQRI7kcSKaOZcBARBwGDKTPG7kAuK0luapWkTJvBz4QEQuB44H/LemAFuaxNEkzgcuAtwCHA++VdHhNsg8Dv4iIVwBfBD7b2lxOrIJl/legLyIWA9cDn2ttLidWwTIjaT/g48A9rc1hdRxIpo4TgavTx1cD76yXSNJrgRcDN7coX1UatcwRMRgRj6aPnwCeAka9QGqKORLYGBGbIuI54BskZc/LvxfXA0slqYV5nGijljkibomI7enTu4GeFudxohX5nCE5Efwc8GwrM1clB5Kp48UR8SRA+v9FtQkkzQAuAT7Z4rxVZdQy50k6Enge8G8tyNtEOhj4Se755nRZ3TQRsQvYAhzYktxVo0iZ8z4M3Fhpjqo3apkl/Sbw8oj4diszVrV9JjsD04mk7wIvqbPqfxbcxR8AqyPiJ+1ysjoBZc7281Lgq8AHI2L3ROSthep9WLXDJYukaSeFyyPpNKAPeGOlOape0zKnJ4JfBM5oVYZaxYGkhSLiTY3WSfqppJdGxJPpj+ZTdZIdDbxB0h8AXcDzJG2NiGb9KZNqAsqMpP2BfwI+HRF3V5TVKm0GXp573gM80SDNZkn7AC8Anm5N9ipRpMxIehPJScUbI+I/WpS3qoxW5v2AVwO3pieCLwFWSXpHRKxpWS4r4KatqWMV8MH08QeB/1ubICLeHxH/OSLmAecB10zlIFLAqGWW9DzgBpKy/n0L8zaR7gUOkzQ/Lc+pJGXPy78Xvwf8v2jvi7xGLXPazPM3wDsiou5JRJtpWuaI2BIR3RExLz2G7yYpe1sHEXAgmUqWA8dJehQ4Ln2OpD5JV0xqzqpTpMwnA78DnCFpXfp3xORkd3zSPo9zgJuAR4CVEfGQpPMlvSNN9hXgQEkbgT+i+ai9Ka9gmT9PUrP++/RzrQ2ubaVgmTuSr2w3M7NSXCMxM7NSHEjMzKwUBxIzMyvFgcTMzEpxIDEzs1IcSGxKk3Rgbtjvv6ezHmfP72xxXprOVltvVmZJt6azwWZ5flHNNiNmgZU0S9LVktZLekTSp9Llz5f0fUn3pzMh/3luH19Jlz8g6XpJXQ3y/05Jn2mwbuv435m6+/uumsxgbZ3Fw3+tbUj6X8DWiLh4kl7/WOCeiNguqR84JiJOya3/EsmEkk9HxDnpsluB8+pddJbOAvtPJPOHnRMRayS9j+QitVMlzQEeBo4BfgTMjYitkmYB3wPOjYi7Je0fEc+k+/wC8FRELK/zenem+x6qs25rRNQNQOMh6YNAT0RcOFH7tKnLNRJrW9lZtKRjJN0maaWkQUnLJb0/PYNfL+k30nQHSfqmpHvTv9eP5fWazVar8c3KXG8W2ADmptOkzAaeA56JRFZrmJX+RZqvLIgo3Wavs0NJC4D/yIJIevX1Xen7cEEuXVda21qbvncnpssvkHRuLt2Fkj4u6aWSbk9rWw9KekOaZBXw3jG8F9bGHEisU7wGOBdYBJwOLIiII4ErgI+lab4EfDHz0LkDAAADO0lEQVQiXgeclK4br+HZajX6rMxXpT+0f5r+2DebBfZ6YBvwJPBj4OKIeDrdZqakdSRzkn0nIobvZyHpKuDfgVcBX66Th9cDa3PPvwSsSN+Lf88tfxZ4V0T0AscCl6R5/grpFC5peU8Fvga8D7gpIo4g+QzWAUTEL4B9JbXzDMZWkAOJdYp7I+LJdOK/f2NPzWA9MC99/Cbg0vTHeBWwf9q8NCbaM1vt59NFw7My10n+/ohYBLwh/Ttde2aB/USd9EcCvwZeBswHPiHpUICI+HX6g90DHCnp1dlGEXFmus0jQL27SL4U+Fnu+euBr6ePv5ovHvCXkh4AvksyDfqLI+KHwM/TAPi7wL9GxM9J5pc6M212XBQRv8rt66k0T9bhHEisU+Rnjt2de76bPbNczwCOjogj0r+Da374kJTVHlbXexHtma32HbnZao8GzpH0Q+Bi4AOSlgNExOPp/18B15IEivwssD8EfotkFtg+kjP8f46InelEhv9CErSGRcQvgVtJ7hiZX/5r4DqS2latHcDza5bV6yB9P0k/z2vToPXT3HZXkEyBfiZwZfqat5PMhfY48FVJH8jt6/np61qHcyCx6eRmkkn1gOR+8LUJIuLMNMi8tXadGsxW22hWZkn7SOpOt50FvA14cJRZYH8MLFFiLkmQ+UHav3NAuq/ZJLWrH6TpXpEuF/B24Ad1yv4I8Irc838haZ6CJHhkXkDSWb8zHVyQH5l2A0nweh3JxIQoGbn2VET8LUnzV28uLy8BflgnL9ZhfD8Sm04+DlyWNtvsA9wO/LcxbJ+frRbgxxHRbFbXfYGb0iAyk6Sp6G9HeY3LgKuAB0mama6KiAckLQauVnJf8BkkM8t+O20mu1rJPVsE3A/019nv7aT9Hen09OcC16Yd6N/Mpfsa8C1Ja0j6O4aDUkQ8J+kW4Jdp7QeSEWWflLQT2ApkNZLXAnenM+Jah/PwX7NpIh2e/K2I+O44t59B0mH/noh4tMBrrYqIgfG8lrUXN22ZTR9/CcwZz4aSDgc2AgOjBZHUgw4i04drJGZmVoprJGZmVooDiZmZleJAYmZmpTiQmJlZKQ4kZmZWigOJmZmV8v8B0c2uWbTf+MAAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.fold(3.52254, phase=0.38).bin().plot();" ] } ], "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.2" } }, "nbformat": 4, "nbformat_minor": 2 }