{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# How to remove K2 motion systematics with SFF" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "You can use `lightkurve` to remove the motion from K2 data. Targets in K2 data move over multiple pixels during the exposure due to thruster firings. This can be corrected using the self flat fielding method (SFF), which you can read more about [here](LINK?). There is a full tutorial on how SFF works [here](http://lightkurve.keplerscience.org/tutorials/motion-correction/replicate-vanderburg-2014-k2sff.html). Here we will show a quick tutorial on how to use the correct function for light curves.\n", "\n", "Let's download a K2 light curve of an exoplanet host star. Remember that for K2 data, we can use the same Kepler functions `KeplerTargetPixelFile` and `KeplerLightCurveFile` as they are the same instrument." ] }, { "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:K2/url/missions/k2/lightcurves/c13/247800000/87000/ktwo247887989-c13_llc.fits to ./mastDownload/K2/ktwo247887989-c13_lc/ktwo247887989-c13_llc.fits ... [Done]\n" ] } ], "source": [ "from lightkurve import KeplerLightCurveFile\n", "lcf = KeplerLightCurveFile.from_archive(247887989) #KeplerLightCurveFile\n", "lc = lcf.PDCSAP_FLUX #KeplerLightCurve\n", "\n", "#Remove nans and outliers\n", "lc = lc.remove_nans().remove_outliers() \n", "\n", "#Remove long term trends\n", "lc = lc.flatten(window_length=401)" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEKCAYAAADXdbjqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJztnX2YJkdV6H9nkoHMZMgHTOBClkuyyiwSDcOyJiz4ETIaAohBxQCiiQGMDkFR0Wt8EKKsH3vlQ0B4Xm5cNxAfBCKIBgyGMBAjuIkJ65APIZOwAVlAwxJYmMyKu8y5f3TXbE1NdXf122/P+87M+T3PPPP2V/Xp6qo6p06dqhZVxTAMwzCaMNRvAQzDMIy1jykTwzAMozGmTAzDMIzGmDIxDMMwGmPKxDAMw2iMKRPDMAyjMaZMDMMwjMaYMjEMwzAaY8rEMAzDaMyx/RZgNRgfH9fTTjuttfSPHDnCsccOXlaaXOkMokxgctXF5EonRaZPf/rTB1T1lJT0BuvpWuK0007jtttuay39AwcOMD4+3lr63WJypTOIMoHJVReTK50UmUTki6npmZvLMAzDaIwpE8MwDKMxrSkTEdktIveLyJ0Fx0VE3ioi94rI7SKy1Tt2sYjck/9d7O1/iIhcKSJzIvI5EfmZtuQ3DMMw0mmzZ/JO4PyS488CHp//XQp0AETk4cAVwNnAWcAVInJyfs2rgftVdQJ4IvBPrUhuGIZh1KK1AXhVvUlETis55QLgas0+qHKziJwkIo8GzgFuUNUHAETkBjKl9B7gJcAT8vQXgQNtyW8YhmGk089orlOBL3nb+/N90f0iclK+vUNEzgE+D7xCVf8rlriIXErW42HTpk0cONCe3jl48GBraTfB5EpnEGUCk6suJlc6vZapn8pEIvu0ZP+xwCbgU6r6myLym8AbgF+IJa6qVwJXAmzbtk3bDssbtLA/h8mVziDKBCZXXUyudHopUz+jufYDj/W2NwFfKdn/dWAB+GC+/2+ArRiGYRh9p5/K5Frgojyq66nAQVX9KnA9cJ6InJwPvJ8HXJ+PrXyIbEwFYAr49z7IbaxRpqenmZ6e7rcYhrEuac3NJSLvIWv4x0VkP1mE1jCAqr4DuA54NnAvWY/jkvzYAyKyA7g1T+p1bjAe+B3gr0TkzcDX3DWGYRhGf2kzmutFFccVuKzg2G5gd2T/F4Ef6YmAhmEYRs+wGfCGYRhGY0yZGIZhGI0xZWIYhmE0ZkMsQW8YAJ1Op98iGMa6xXomhmEYRmNMmRiGYRiNMWViGIZhNMaUiWEYhtEYUyaGYRhGY0yZGIZhGI0xZWIYhmE0xpSJYRiG0RhTJoZhGEZjTJkYhmEYjTFlYhiGYTTGlIlhGIbRGFMmhmEYRmNMmRiGYRiNMWViGIZhNMaUiWEYhtEYUyaGYRhGY0yZGIZhGI0xZWIYhmE0xpSJYRiG0RhTJoZhGEZjTJkYhmEYjTFlYhiGYTTGlIlhGIbRGFMmhmEYRmNMmRiGYRiNMWViGIZhNMaUiWEYhtEYUyaGYRhGY0yZGIZhGI0xZWIYhmE0xpSJYRiG0RhTJoZhGEZjTJkYhmEYjTFlYhiGYTSmNWUiIrtF5H4RubPguIjIW0XkXhG5XUS2escuFpF78r+LI9deW5SuYRiGsfq02TN5J3B+yfFnAY/P/y4FOgAi8nDgCuBs4CzgChE52V0kIj8NzLcjsmEYhtENrSkTVb0JeKDklAuAqzXjZuAkEXk08EzgBlV9QFW/AdxArpREZAz4TeAP25LbMAzDqM+xfbz3qcCXvO39+b6i/QA7gDcCC1WJi8ilZD0eNm3axIEDB3ogcpyDBw+2lnYTTK50BlEmMLnqYnKl02uZ+qlMJLJPi/aLyCTwvar6GyJyWlXiqnolcCXAtm3bdHx8vIGo1bSdfreYXOkMokxgctXF5EqnlzL1M5prP/BYb3sT8JWS/duBp4jIF4BPAhMicuOqSGoYhmGU0k9lci1wUR7V9VTgoKp+FbgeOE9ETs4H3s8DrlfVjqo+RlVPA34ImFPVc/olvGEYhnGU1txcIvIe4BxgXET2k0VoDQOo6juA64BnA/eSjYFckh97QER2ALfmSb1OVcsG8g3DMIw+05oyUdUXVRxX4LKCY7uB3SXXfgH4/ibyGYZhGL3DZsAbhmEYjTFlYhiGYTTGlIlhGIbRGFMmhmEYRmNMmRiGYRiNMWViGIZhNMaUiWEYhtEYUyaGYRhGY0yZGIZhGI2pVCYi8sjIvi3tiGMYhmGsRVJ6Jv8sIhe6DRF5FfDB9kQyDMMw1hopa3OdA1wpIj8LPAr4LNnndA3DMAwDSOiZ5MvC/yPZ90ROI/vUrn2D3TAMw1iismciIjcAXyVbpXcTsFtEblLV32pbOMPoJdPT0wB0Op0+S2IY64+UMZO3q+pFqvpNVb0TeBoweB80NgzDMPpGZc9EVf8u2D4C7GhNIsMwjDVCv3u7/b6/T4qb69uA5psPIfta4ryqntimYIZhGMbaIaVn8jB/W0Seh0VzGYZhGB61Z8Dnbq9zW5DFWANMT08vda0NwzAcKW6un/Y2h4BtHHV7GcaaYRD8yobRS+bm5votwhIpkxaf6/0+AnwBuKAVaQzDMIw1ScqYySWrIYhhGMZao9+93YmJib7e36dQmYjIn1PizlLVX2tFIsMwDGPNUdYzuW3VpDAMwzDWNGXK5N35BEXDWKLf3XrDMAaTMmXyr8BWyFxeqvqrqyOSYRiGkcIgGXdl80zE+/30tgUxDMMw1i5lysTmkhiGYRhJlLm5niAit5P1UL4n/02+rap6ZuvSGYZhGGuCMmXyfasmhWEYhrGmKVQmqvrF1RTEMAzDSGeQlp+HLhZ6NAzDMIwQUyaGYRhGY0yZGIZhGI0pW5vrDsrX5rJoLsMw1j2DNjYxqJT1TH6CbPn5f8z/Xpz/XQe8v33RDMMwjCKcchuUj9VVRnOJyNNV1Z8Bf7mIfAp4XdvCGYZhGGuDlDGT40Xkh9yGiDwNOL49kQzDMIy1RsqXFl8K7BaRE8nGUA4CL2lVKsMwDGNNUdkzUdVPq+qTgDOBSVWdVNW9VdeJyG4RuV9E7iw4LiLyVhG5V0RuF5Gt3rGLReSe/O/ifN+oiPyDiHxORO4SkZ3pj2kYhmG0SaUyEZFHichfAu9T1YMi8kQReWlC2u8Ezi85/izg8fnfpUAnv9/DgSuAs4GzgCtE5OT8mjeo6hOAJwNPF5FnJchhGIbRNZ1OZyAjuaanp5mbm+u3GEukjJm8E7geeEy+PQf8etVFqnoT8EDJKRcAV2vGzcBJIvJo4JnADar6gKp+A7gBOF9VF1T1E3na/wPsBTYlyN8K09PTAxNFYRjGxsIpkomJiYFRdCljJuOqeo2I/C6Aqh4Rke/24N6nAl/ytvfn+4r2LyEiJ5GFLb+lKHERuZSsx8OmTZs4cOBAD0Q+yuHDhwE4cOAABw8e7GnavcLkSmcQZQKTqy4bRa7Dhw+jqhw+fLjrtq3XMqUokwdF5BHkExhF5Klkg/BNkcg+LdlPfv9jgfcAb1XVfUWJq+qVwJUA27Zt0/Hx8WbSBgwPDwPg0u11+r3C5EpnEGUCk6suG0Gu4eFhRITh4eFG6fZSphRl8irgWrJvmnwKOAX42R7cez/wWG97E/CVfP85wf4bve0rgXtU9c09kMEwDMPoAUnRXMCPAk8Dfhk4Q1U/04N7XwtclEd1PRU4qKpfJRufOU9ETs4H3s/L9yEifwicSMKYjWEYxnpmkMZLIKFnIiKfB16vqu/w9n1YVX+i4rr3kPUwxkVkP1mE1jBAntZ1wLOBe4EF4JL82AMisgO4NU/qdfm+TcCrgc8Be0UE4G2quiv9cY31jq2jZBj9IcXNdRh4hoicDfxyHkl1asU1qOqLKo4rcFnBsd3A7mDffuLjKX3BGivDMIyjpCiTBVV9gYj8H+CfReRCSlYTXs+Y1WsYxiAwiG1QijIRAFX9UxH5NNn4xcNblcpozPT0NIcPH2bXLvMCGobRPinK5LXuh6rOiMgzgYvbE8kwDMNYa5R9HOsJqvo54Mv+ulk5H25XLMMw+o25dY06lPVMXgX8EvDGyDEFzm1FIsNogDV8axtTYGuXso9j/VL+/xmrJ45hGIaxFilzc/102YWq+re9F2ew8a0ls6AMwzCOUubmem7JMQU2nDJZS3Q6nZ4vbrnWMIVvGKtHmZvrktUUxDAMw1i7pIQGIyLPAc4AjnP7VPV1bQllGN1gPZHeYvk42AxaeU9Zm+sdwCjwDGAX8HzgX1uWyzCMDUjYMA5ag2kUk/Klxaep6kXAN1T1D4DtLF863jCMNcggfi10enqaqampgZMrlUHM09Uixc11KP+/ICKPAb4OnN6eSIOLbyWZpWSsd6xX0B/War6nKJMP55/JfT3Zd9eVzN1lrDPWaiEuYr08Rz/wvzHezbUQz//1Vsaasp7yo1KZqOqO/OcHROTDwHGqOpgfWjYMo+/Mzc01unZ+fr5RGt3Qz0Z9bm6O6enppXuv9rP3ipQB+GOA5wCnufNFBFV9U7uibWz6bbH0+/7dsJZk7TdNeh5tsrCw0G8RWqWX9WrQynuKm+tDwH8DdwCL7YozeBRVurYb27VonQySAuqnLKt177bvs7Cw0PdyOAhlqW16nc/9KvspymSTqp7ZuiRrENc93bFjR/XJa4CmBbrfDY9Rj0H7hjjA6OgoCwsLPe8xrVYD2036i4uLzM/PMzU11TM5QtfZapCiTD4iIuep6kdbl2bASX0xvSi4/ejuu3sOqgvE6B2+4i8qr6Ojo6VloJtybgZHlgdTU1Mr8tbVv9HR0eS0yt7BaiuUFGVyM/BBERki+x68kH3C/YRWJRsw/ALQ7cspevHOIpmZmWkmZEP8QjwILo4yBsmlVsRakLENyhTQWjZQUt9n1Xm+odjpdDj77LOBrIfSCxlnZ2eBekqpF6QokzeSTVS8Q1U35Lffi6hqcMNC1VbjvH37dgD27NmTLEsVftc7Vcmt98HTfnD55ZczPDy8phRSE1mdsllLz9uN0bC4uFjYHrg88NOtcw+nlFa7PqYok3uAO02R1MMVlBRLbH5+ful3ndmzqzHTtk6B7IVl1Q2xirYeBjTbJKVcdtM7dS5Sd4+6edZGHjtrvcxS7/V9i1xMnU5nyfgru3ZiYoLZ2VkWFxdr1fN+ehNSlMlXgRtF5CPAd9zOjRIa7KyC8CV1Op3CATMXK982Vb2iw4cPs2vX6s0vHRpKWZ2nXVzF66c7xTUgvVb2oVJbjTG8hYWFWn73ubk5FhYWartYeqWwn/Oc5wBwyy23NEqnLinjE2F5cHkbEhpwReOXgza2maJM7sv/HpL/bSjcC4tZ6G1FnEA9K79XPYLY8/hpt2mhh2mX3WtqaoqFhQUmJyd7LkcTBrEHU7fBiRlOrqF01Hm+qjxJVbhlLp+UfK8KJmiC60E4Wdy2n29+xGfd+hrrIfrv1c8Tn4EaM8knLI6p6m+vkjwDS/hiYpXUt4r37t0LFEduxHAFxldcRdbO9PR0qQU4NzdHqmcyVlGHhoYKC33sfMgqSVu9k6IGo8iF4fKy2wZ+tRRDm8EXMcUQu1/VM7rGLCzDRdctLi4WukfDa2K966Ly5Y75svjPWKaY5ufnlwamw3t0Q5FrdWFhoXDMIjRKi+ru6OjoMneYn17YHjiFFcvH1Q79LlUmqvpdEdm6WsIMKq4AOEu4ypryC627ttvuaLc+0FhlruqK+xV1dHR0xVhO0TPErCZY7iIMC3ZVcEJZA+GezW9sFhYWlnosTSwyl56r+EUVNYW2KnJTBeTyKtbYOBetMwq6saJT8N8bVFvR/jhMylhObHmSqmepUjB+T8iV6enp6aV0Y67tUM5Qhl6NcRTV99Ukxc01KyLXAn8DPOh2bpRvwM/NzS2zuP1C7SrjgQMHChtbNwkL0hoX37JJoWhsJkyjqCvs43fPIeud+JXcWXe+e8ldMzY2tuz80DVSRdgbC5WCr6BjMg8NDa3IuyKLOlWetgIK/HzxG+4YO3fuZHx8HDha9tw7b2pdVxkXvkFRNzLIlaO6eR/2NIsaRN94iPWIZ2dno+UvZWKg70at40orYn5+ngsvvBARWdo3NDTExMTEsucbGxtb0XPx63dRBFjYRvnpryYpPomHky07fy7Zd+GfC/xEm0INIikNSy/mZtSxquta4HNzc8sq2dTUVLRy+dZiUU8kfE5Xqd0ifUV+Xt/y948vLi4uRa64yhEu+Le4uMjU1FTSu3C9q27x7+EMhxD/eeowOztbqBx7hT/WV6YIQqXvzi1zVZXd0+WbU+7uWVPzyb37GGX1K6VMhAonNgA+OzvL/Pw8i4uLK96R71JyspS9R7/8HTp0aMVx32hzsvl1zu95OqVSFtjj50FZ6HFbpKwabN+Cz0l9Oa434v6nWrh1XTT+2EwVztoPC1mR9e6fW+Tq8RuPsEfgGqOi4IW6hHkY3jfWI+k3Ze6oWGNVZf26dxS+89h9/EHhoaGhrixqVxaLGrDQnZlaFn2cjO7Z/HHDUMk5A8PRTcSkbygtLi6yd+/epfv4eebL5uTyy3MoS9X9QhYXF1c8M7BsQL2oDIcRi0Xn1Z0n1pSUVYM3AX8OPJ3sWyafBF6pqvtblm3gcN3SmNsqLIghYaGKNRxN3CopA99Vk5n8htlZS76ScW4sd32qde0azsnJSebm5rj88suX5AgtQ99aq6LI2u3lpK0ypRQLF4ejjXuqCxLS5kIUUfWcVeXSH5+KGQVwdFA4ZHZ2Nhpg4q4bGxtbtr8snL5IPnf/qsY77MWGafj7U8tG2bP5FAWsFMldVr7rjtOFrrAyedokxc11FXAt8BjgVLJVhK9qU6i1gu/mCC0sZ7mMjo4uNcJFbiW/YIUFt8jvGVtHKUbsnmEhiw3EOreTc0f51lhVIfWVUliR5ufnueuuuwp7bPPz8yvu4SutOtao/wzuPRW5pi6//PKlSuxX5Pn5eRYWFpifn+96DKgK9+5iriUnl8PPi7m5Ofbu3RvNx6J3NDMzkxxSHbr6HNPT00sGgXvHRb2Sovfl8tRP393D1Zcyl1eMsK64d+nS8J+7jrfA/U91/Q0NDTE2NsbMzMyK832Xln8Pt7ilXw+L7uWXT//8mIJazVnwKcrkFFW9SlWP5H/vBE5pWa6BxL1Ef3AwZkH43Urnu3dRKHVebpV7wm9YipSOu2eVhV3UUJbJXNfy8WVwSrYMv5Gtiy+za+y2b9++TIYixeJf61t4Lp9C94u/PTU1FVW4ZeMrMcvScccddyw11GW936o8ct9VdzK4MhpeV/ROehl1FMuforrkX9ft/SDN1RczyJwxUnfe1/z8PNu3b4+mOTo6uqQs/WcOx2SasNpzTCBNmRwQkZ8XkWPyv58nG5Df8DgrGpZXwrpLSfsv3i/wVRZk1TwQd06RW8UdD2Pw4ah7wl3vIra6JbTq/J5JzF2WOl/Fl8m/pshajz0rwF133bVs8DWWjj/o6uef7yaKWeLbt29ndna2UGkXPat/7vT09IoGok4vzc/7UA5/+fOiMuXfa25ujn379iXfO/XLiWH5aEpsDDKW11WylSn7qv1Vz+SMTWdwNlEkzj09NjbWl1nxKTX2JcCFwH+SLa3y/HzfhsIvhEWWfsz6cQ1UWYHt5sWH6/tUWV2xhmdiYiI6nyClIS+zfFxXPkwn1kiXWdupA4fuXt30lPw5CCnuO6cEQ9dblZ/bj05LncMRWq0ppEaYddNoTU9PL4XNlhG+9zphqn4Unh9Z5o8F1MV3M9cZ12yKL/MZZ5yx7FgYIVjmPXBh9z6xnpwfmt/pdFZ9eaPKu6nqf6jqT6rqKar6SFV9nqp+cTWEWw8U9Qr27t0btZBjM35T03WkDIyHLrsQp2iKFEMq3V7nnq/K+naN++jo6FKlq7qnC0FuQsxNkzLg7vd86uRNOMZQRJG7KBz/CuVyYbwpcrjrymbAx+6dkq6Tz1nrqbi8DHunLr99L0KIa9TbcA25NF1PLhwL9Mcjw+smJyejvYwixeyMHed6rZuHTSkszSLy2pK/16yahAOEe1Gx+RfdFP6UWb91LFJHioXuGhS/gvmx7r6f2FmG4YBpUWPou8ZChoaGGBkZSXiiNHdhGJ2T0sNw1B3D8u8Hzddm89NqOls+RTH5wQ3hc69W5E+V67bsfYyNjRVe76LNYmWrKm9CJVvn+pjxEnoxFhcXOXToUKM8dmOMQ0NDS6H8bga+fz+/HSqaH9UWZbn8YOQP4KXA77Qs15rDf7mxQulm1PqVIRYaC1mlKRufKLOq61jcRUs9hIUwrDBjY2N0Op2uKsfo6CibN29ecc/VxHcz9eP+Mfw10dx2iq8+pMwCr6Kqx+RHjm3evLm0rMXSKJuLEnu+bmbPp6bt4yL2/OdvOj7o4/ea/MmcfpSZfz8/wMdFeDn3X8wQCI1S/5swq7k2V6EyUdU3uj/gSmAEuAR4L7C56Lr1iOtFuMbHjxBxBcWf1ZwSquk3ZP4L37t3LxMTE0thhb1YcydWsVPDHGFlr8stZZFy39igcZ3B2zKrMJxXkNqIujGPJu678P5NPjng5t842bqhrjuj7D6pMlSNEXVDWMcc4XJGIW7MITbTvK5iCMtS2bOEdWNsbIytW7dGxznC82NjSWV1fmFhYUWdcmOyRfdbTUrvLiIPF5E/BG4nm+C4VVV/R1XvXxXpBoCyZTnCngawbIA1JFYRYktN1BkUTrlHTP6yAeu9e/cusyDDnkqqTEUK68EHH4ycvZJYCGnTChP2Spr0TNqK4e9mPKeuLHXzsWgOQ5FCDveF61AV3cP9OVero2hpH/99xiIrIa7ow3WsquQq21+W992Et4cBCP7+VKOh7tp4vaBszOT1wK3At4EfUNXfV9VvpCYsIrtF5H4RubPguIjIW0XkXhG53V+dWEQuFpF78r+Lvf1PEZE78mveKv7KaS1TFPk0NzfHhRdeWGhR+cSiutygY+hnLSsIZfcoCj8Maaqk6vQCeoHfyJSd0408RXldZdHWjYAromjxvjo0zZei64os3qrepQtTdczOzibN1A/H59z9Q7er+x/OVxoZGWFycrKnVnrMoBkbG4vmuXvGlLXhXA/Zvesit55zefoTKGG5R6PpWnS9oCzHX0U26/33gK+IyLfyv2+LyLcS0n4ncH7J8WcBj8//LgU6kPWGgCuAs4GzgCtE5OT8mk5+rruuLP2e4F50nW+nl1EUJlx2XTiLPTzXb/Ri6fhROv3sCnd7b1dxnXWaYhHWIXWeS5n8YeNZh/A5UhreEDde50iZc1OFv7ZcmeJ0x2LK1zcCUg0L//2G68nFIpuqGtEiN6//u67R4+Qom5cUMjo6uiLc3R+nLJocGX7Souh+vgehaD5Vm5SNmQyp6oiqPkxVT/D+HqaqJ1QlrKo3AQ+UnHIBcLVm3AycJCKPBp4J3KCqD+Q9oRuA8/NjJ6jqnvx79FcDz6vxrLUJJ435VFmOZVZtakOR6tMv+ghOrCfVTU+hzpyEMvdBHcvJzz83qOysVWeN9qJnUGd+Sso8FJ+tW49+CqjqXcbWj2pCVRptGBUTExNs3bo1qmyKerJF5/uEPZOynsehQ4fodDorypoLaunlc/vh876s4TlVPaXUtdX8AIG5ubllBmadyM+2SPmeSVucCnzJ296f7yvbvz+yP4qIXErWi2HTpk0cOHCgtoB33333kn//7rvvXnZsfn6ew4cPo6rRQb+yAhJWFvdNC3//4cOHV4S9Fj3D4cOHVxwbGhriuOOOa2S9uzBe98XGWMN7/PHHL4U9ut+OkZGRZeMjZbKEaRd9JXJ2djYa/ulXtDqNcSijv//000/njjvu6CpdyN5LKF8qZefHjAT/XqnpxxRy1Ttw57j3PDIysvT+Dx8+zL59+5alUZS/jrKVhv3enrOyR0ZGlqIBR0ZG2LFjB8997nOXrjnuuOM4cODAsjI7MjJSWvb8MlzF8ccfD6zMb5cvTiYft46Zqi7J5uebS8/JrKps3ryZffv2Fb7X8N247bDMnXvuuVxzzTXRNA4ePFj5vHXo5/B/bLxDu9gfRVWvVNVtqrrtUY96FOPj47X//CEZEVkx2H7fffcRG7bpxrL0G8PR0VGGh4eXHVtcXOQ1r3kNL3jBC5b2O2treHh4SWbH6OhoVLY6jI6OsmXLFkQk6uoYGxtjy5YtS/t8RTI0NLTi/lUugLJr/eVdwmNOlq1bty57R0XWoG+hxgwBt99/B3XnEA0NDS27vtf4zzk2NsZ9993HXXfd1SjN8P1u2bJl2WRQ/xxXHg4dOrTUUxweHkZEllnrRfnrE4bihvfz/7Zs2cLw8DCTk5N8/OMfZ3x8fNm7ERHGx8f5+Mc/ztjY2DIZi1xtrgwPDQ0t6ynFyo+IICIMDw9z3333LaXn8uDQoUO85jWvWdaOuHLgZHNp+LIMDw+zZcuWpTrl8tIvQ36eHDp0CBFZagNEhC1btrBr165l17h7xv5OPPHEyjawDv1UJvuBx3rbm4CvVOzfFNm/KpQtd9CrQWaXlr96akwOH9/F5bviJiYmkqOmHH6F3rp1KzMzM6XuM4dfmV0jV+TaSB0kX1hYWCZPag/Lf0cx90LoMilSEr14p6luh25dj45uJqeVvQd/cmBszpFPGHQSjiEU5W84phOWlaI8KXPnDA0NLc1/8ed7hYZK+NvNn3FjM/6gfugai5XDlDD7kZGR0nlo7rn8b5SE4x51xnZ6Pa6YdM9Vu9NKrgUuyqO6ngocVNWvAtcD54nIyfnA+3nA9fmxb4vIU/MorouAv++b9JTPoPYLZEoD6p9TNNmoaLnxbqjy4ce+Y+/OD0MUJyYmotEtVQOf4XnhOFNR3qZWkE6nw549e5ZZm/5grpO916RE9sWuScGPaNu6dSu33HJLaXBIrPEsi4pzyjZssMtWeXayu/WgXPop5SFF7hh+6Ksbu1hcXFwWZeZkdgqxarKfwx/kdgP8YQRXOI4Ye163lpnLy1ARuvIehj1XfdPHf39+9JhPLw3cVFpTJiLyHmAPsEVE9ovIS0XkV0TkV/JTrgP2AfcCfwG8HEBVHwB2kIUl3wq8Lt8HMA3syq/5PPDp2hfyAAAd30lEQVSRtuSHeKMVWlQprpvUyBlHanx42HNItYRdYQxn3hZZ7GFXHo4q0pjicw1K1XO7LnrZMt0xQneMO68sCmpycnJpMllY8YqCJeoslpeqOB3+xLY6jWyZu63ofrExlomJCfbs2VP47K7HG65iXTRnqWgWdtjglr1zR1Xou99bSvmQVNiIuzx3RpC/31ceo6OjS2XKRbb5IbruGidXKHdR4Eqsp1R34c1wNeJ+D75DiwPwqvqiiuMKXFZwbDewO7L/NuD7eyJgAqOjo0uD4+G+Nmka0lc2VuAosshjSmBsbGypEjkrqhcyOgvQVdxYZXIyu3v6is1/jnAGuZ8H7rqxsTEWFhZWDPpOTEwsC7aIyTozM8P27dujA9f+fVPpdqnxomXKi1axdvnkl2M3v6HIaPGX1G+Kny9uSSFH1eB70X7fHeTShaPv3JXNIoXlylrZ/B5Xzp2MscUY3XnOtVekmBcWFpaCRsJl8cP7dDqdws8BuPLrrpuZmSlUpt0EjDSlv/Pv1wBDQ0NLFWB+fr7WUtpliyHGzncsLi5Gv1+Rilv0ropwpdu5ublSa9W5bpxfOVW+KvdIqCB8/JBgNwbiFJGbs9PpdJiZmVlRSUOKGoROp8PmzZsrewllrh73nHV7GSlfPUxJs8gCDsuV//yuEYq9b99fXzXXypULv8zGGrjJyUk6nc5ST6Eov6oawZSZ4W4cJnRX+ddUfVDO4doAX9bY6gxlxsGhQ4dWuM/c/1i+lrnPY2XbGRnhMkcDsWqwEW986ljksXj3snPrEFZCf2HAKou3qAC7dIqu990TrjFPtV7DOP9YRaxS1Kn5GQ6yV7F9+3b27du3VFGLGjmXL93MsC47vyoMOOWZmzQaobvHX4E5VArhuUW4Btl/bjfO0c1qzRBXeuGS+ilBJ76brEgO/9lj+eCnFW77davKFV532ZNu8241MGWSQJOFFev0TnqFv5Cc++++Se0X/rIoHUeoXPxZyZC5KqamppZVUH88JrUx9ClSBk5hxipouMpuuAik/9/HucD8EFbf+vN/x+4bNnC+onH57s7pthwU9Tr8KL6ilarn5uYKPzLmL8tRFEQRUuWecz3EskY0FmZeRCySqmgQ3OEbBK4XG36/x237vd2y3op/rU+n01kWtp6y/ljso1hOBv9dhiHTIX4AQpjeuhqAX2/UeTllEUtlNBmIjRXikZGRJd+p33A4v2wYduhb3q4B8gtmbOAxJKzgYbe+Tk/NKWK35H1VlIujalFAn5iiCcNh/c/aOtmKnsM/5vfiXPpVnxfoxi0aS2Pr1q1LYwMxy9d/p9PT00uuQn8Aulti7z9s1CFt4UJ/4qL/jjqdTqHbMrzelW33Xv0yGpZvVyeqnj/Wew0VRTi2Gj6vk78oyszPOz8P/GvC4AJ3bLUNWDBlkkwdZRCG+RZFwYQ0cVXELLTNmzeXumT8ili1cF1IkayuMQqVilMORTirLtaTmZ+fX/pMccx/7X8noxeVyCncMpefayhiqym7ZUX8HoGbu+EatKK5IeEYysjISNflwh8nKFtKpI6rJdaA++OKIb1s1MK0wsH98FjKfJcYrtEPB/bde3X1uSrNUF7Xq5+dnV0W5r+4uFianjvmeznC0HyXvjvm3sdGmWeypkgJ83Whj90OTHc752F6erp0klfsWFUhcxO/fOvIFerY8/kWku/XDa21Os8Yk91FIsUGkiHeK6nqDTjKPixVlF+x8GE/L2ZmZpY1tE5RxUKq/bW84OgkPH8bVvrlw3yKNezhgHFYTsN5G/4cDXd90Tdsytw7rk5MTk4uDYqHVn0YKlukNMJebux6yPK8yLBz+RYLUimaG1NkZPlzPPz7hAuzwlGjoI5S8+fItDEfqteYMikh7O5XWVlFL9z3Sac0akNDQ8kT0WCli8vJsXPnzqXfVe4Av8AWzTLupmC7hmZ6enrZctvuOfw/WK4MXHc9lu9+w1qlGF3F78bCL3L5+HKFjXevKv/k5CQ7d+5cup/7iFbdeRWh3L6cEFeUvhJ0jWVRr8TN3fEb6pAwcCXWaPurQxfJ7lyePkXuLl8ROHde2AOLzZJ37y8mi8tb1ztx6aUoiti9QlyP3vUuwmeNPaeTyS+r4RjRamDKpATfTeDHklcphFi0lO+TDhtIl16n01m2vlS3DVJKZFjRAKi7NuY2KGvEwgY0ZekHVwF9N6DLO9crKnNldIsbUygLOljtDwtVkTK463DKsyg8NzZxE8pDgIveZ9mcDV/hhy6n2IztIku/KALQT6+KxcVF9u7dG53ZDyvLlu+W9ct2GAmW+k4efPDB5HFX14NzRpg/gx6Wr1Lgy+23V/2I+DJlUkJYUMosZDhaAMt80G6Jjz179nRlNTjrrUzRxO5d1HDGuvEpA5vuWtcw+dZw6gBukQUfTsQK56E4HzMU98p83BcxXSPlXBwx90qqj7/MEm2yHHisIXe92qJ3Htvvyog7Fj5X0cS+uqGqsHJtLV8ZxBpcZzA4+apmiTtXmUuvDjG3Zd2eY8r7DCMeffxylRLEElIWyh/K2c+wYVMmCcTCe4u2qwaaw3RjjYSr0DGfeq99p1WKyZFa8J3MYXRU0b1cJfX3u+9ah6TIWVXpfQvYr3hnnHFG9PyJiYmuFb+jrEz47zjWm02hjuIK506Ek/l6gQs6Semd+oyOji4tERTmV9mqE370UlU9nZ+fXzZ24uq2kzNcoqhM8fg9iHC/b5D5+eArXj+E2Y3XhYP+4Tpy4X5nJIX5UTS+1Cb9/J7JwBNa3P5+YKkyhv7g1Modayyci6dKLv/alKU1iiq1izjy75nSAFT1XvwK2pSYwvVdIuFM/ircCgM+O3fuXLa8P5Q35kWV1EWANWmgyxR3XYs6dk2YV355jl0/MzPD2WefnXxfR51eXiiTC18vCpuFo8v8+GmMjIyUBqQAS4PjruyHDa9fNlLrcjhu5kKaY3L4vXlf9vD6+fn5qEvdV0iuDXLp9WN+icN6Jl0QVtI64Zt1/fGuMPlzLYr8vjFiDWL4BckUecJ4+HDuSehXLsKXPbRAY935Itei70pLbWD9yKtYxSuzQp1VWDVeVtSAu2NFLsgwuqmOMvLvFQtHdYOzIWHotlvVoBvC+zo3YmweRdH4SphnrtcUegbCPFxcXIx6D5oSjpt0Oh1mZ2eZnZ2tvWJFEf7zNZE/5uKqCuTpNdYzSaRM419zzTVLH5LphWUaEqtoPhMTE6WL5qVEkIX3i0XtOFKWlPFDil2BdgslOtzYSNmzOattdnY2av25Gd5h2j5+hfLnD/SSohDSbijLj1jjEIZGuzRijfYgLMUR9oSKZPXPreoFux590ZcJ/TJdlAcxZe96vf4YU1mYtE+RzOH4Svj+yhS6K7/hBM6pqakVQTWrGckF1jNJJjWst9eEPu26g4dlE++KiA2g+9ZyOKch5QNadY+7c/zGxO9ZhJZoN+tl+Vx++eVLvb6i78kUkbJmWJ2xNFjZsKUOwrryEev1hbg89V2rYejwxMTEMjn8QAV/mZimDVfV81VZ7mVL0Y+NjbFnz57SMPOqstPEfVS0ppg/8TWc2e7k8omVobK6vZrjJWA9k1LC8YtU/2nqoHaqDD5VMpRF/ZQRDvS57bAidDqdZb2AsnDScGZ0aoX0rbNOp7PMX++nF+az265qdGOuFygf5HX3Lkrbt7ZjE+KK8OcpFBGm5a6JuQSL7uv7551lG2tAq2R3cvrfIC/LF6doyvI2Na/C913VI3f3h/Ll9FMjF31iocRAtMdSx5jzZSlb28xnamqq9c9ipGA9kwbUtWDbTLfMugorYdmExKL9/vGqyldlqfuWbKzx92cw1yFFkczMzDReeHHQcFGBsZnv/rt3Pb1uexH+jOzNmzdX9pL9iY9hb6Yu7jn8McNw9YGJiYmlSZ5l+NFV4dyNbqjqAYbvJZb/zvMR5k9dmVw+1zFoesX6qE0tUrdhL7PU28Jf/yvWoKa6SPz02uwiu+9hh/css25TuvyppCpFnzDwoUjeum5I6E1PNvxeuJMRVhoPYRSTOyfMD/8Zm7qy6gYVhBSV4Vi6biXfUGZ3bpXLOgwd9qnrMi4yWFJmw5fNDSuiaT43wdxcNfALVlGUVNXge1MlUzZfITVKzLkHfLdIeG24v8iNVRb9VHSdv95UEUV5WJTvVQP5Pr4Lz3cP9KKn0gsjos67rKKoRxp7VveZgirZAF72spct2y47v6hsVb2vWNmMRXzFiD13UQg0pK8wMDk5uWJR0ao8SFnt2pfNBQvEjB3nZi6i18E/dTBlkki3ft3VpE5DlmppFikRfxZ66nXOGty5c+dS9FsRfh4WrUhbpAhTcW6TXkQ4pc5QLpKjiRKqamB6cZ+wXKe8wypSZel1L79o3KWo3objkFV1p+y4y0fXi3TvJLUMV82h6SemTNY5bbjb2gqv9XHhjlNTU4XWcq/Xz+qXRdcL6soenl9lADVVeKtNN4ZGirLvdmJgzOXmy1lEVQBOrC72y5g1ZbLG6PVkqdR0Y+6uug2MUwoHDhyoI2oSdWTxo9X831Uunjaoir5q+z5tX9tvfNmLFEtVTyNc3shfaDH1fnffffeKAILYdc5lm+Kuii342E+DyJTJBqRO4xCrgGWTzNYr/bL2YvSicS+bod9r5dHUHRlLq5u0Y8/lesDdkDq+ND09zebNm5e5s3pRf/xlYPq9yCOYMukp/YjkqqJsAD2F2PIgvRwgbkLZ0iW9phfvtM3KXtSwFckda5QH4Z32g7KB+di5vVAEZYootTz7M//reBnawpRJCXUb4kGyXtukqYLqNXXlSHF9tMFaHpPpBYNSXkJ62XMqSr8N1274GYZ+tz+mTIyBpI01zsL04ahFXrSmUy/pd2VvShMjYtAMkJA64eWrEYxQttJC+NnqQcGUiTGwlC3TAc0a53ACYhuWo6Mbi3etK562aFMpxZTEarv/Uuft+BNUB6WsmDLpIet5UHpQLcq1QDduCJffbSq5fjHovZQQX86mMrcVjTkImDIxorgKP2gF1metKO9+R9mkslYa943OoL4nUyY9ZBAb3kGxpNpgEPO7jEHOS8NoiimTEqzyx7F8qccgDpb2k7XSo1wrtB2NloopEyOKq/D98m+3fb+11qsZBDbqLPq1Qr/z2JTJKrPWBh8HmSYW7mrm/3r5bsogYPVmcDFl0kOarBw7aKyVQeO1wEafrBjSRq9wPUe/rRVMmRhRwgZwEHtUa8VVtVbkrGIQy4AxOJgyMZIYxB7VIMpkVGPvbX1iysQwWma9WPK9UgLrpadmLMeUiREl/Oxut8t0t0mTcR1z2axN2npvVg6aY8pklTGrrHfYwLYxqPTDWOm3gWTKJJFevai1YgGtFTkNwxgMWg2AF5HzReRuEblXRC6PHH+ciMyIyO0icqOIbPKO/V8RuTP/e4G3f0pE9orIrIh8UkS+t81nMDImJiasV7XBsTJglNFaz0REjgHeDvw4sB+4VUSuVdV/9057A3C1qr5LRM4F/gT4BRF5DrAVmAQeCvyTiHxEVb8FdIALVPWzIvJy4PeAX2zrOYyMQeypWMO2NmlSliwSbHBp0811FnCvqu4DEJH3AhcAvjJ5IvAb+e9PAH/n7f8nVT0CHBGRzwDnA9cACpyQn3ci8JUWn6Hn9NuvaRgxUsqllVmjjDaVyanAl7zt/cDZwTmfAX4GeAvwU8DDROQR+f4rRORNwCjwDI4qoZcB14nIIeBbwFNjNxeRS4FLATZt2tR4Zqz7El8snYMHD/YknV5TR67VpFdy7dixA+guL8P3sN7zqoq65bJf+XX66acDxXIOynvsR/nq9ztsU5lIZJ8G278FvE1EfhG4CfgycERVPyoiPwj8C/A1YA9wJL/mN4Bnq+otIvLbwJvIFMzyG6leCVwJsG3bNh0fH2/0MLt27So9npr+8PBwrfObslr3qUu/5Yq9z37LVMRqyNVNuexHfqXIOQjvsR/lq6qNitFLmdocgN8PPNbb3kTgklLVr6jqT6vqk4FX5/sO5v//SFUnVfXHyRTTPSJyCvAkVb0lT+J9wNNafAbDMAwjgTZ7JrcCjxeR08l6HC8Efs4/QUTGgQdUdRH4XWB3vv8Y4CRV/bqInAmcCXw0v+xEEZlQ1Tmywf3PtvgMPccGEA1j/bMRx0ZbUyaqekREXgFcDxwD7FbVu0TkdcBtqnotcA7wJyKiZG6uy/LLh4F/FhHIxkV+Ph+MR0R+CfiAiCwC3wBe0tYzGIYxWGykxnmt0eqkRVW9Drgu2Pda7/f7gfdHrvtvsoiuWJofBD7YW0lXDwtnNQYRa6SNpthXewzDMIzGmDIxDMMwGmPKxDAMw2iMLfS4yphv2jDWPxuxnlvPxDAMw2iMKRPDMAyjMaZMDMMwjMaYMjEMwzAaY8rEMAzDaIwpE8MwDKMxpkwMwzCMxpgyMQzDMBpjysQwDMNojKiGHz9cf4jI14AvtniLcaD97/DWx+RKZxBlApOrLiZXOikyPU5VT0lJbEMok7YRkdtUdVu/5QgxudIZRJnA5KqLyZVOr2UyN5dhGIbRGFMmhmEYRmNMmfSGK/stQAEmVzqDKBOYXHUxudLpqUw2ZmIYhmE0xnomhmEYRmNMmUQQkceKyCdE5LMicpeIvDLf/yQR2SMid4jIh0TkhHz/Q0Tkqnz/Z0TkHC+tp+T77xWRt4qINJDrOBH51/wed4nIH+T7TxeRW0TkHhF5n4g8JN//0Hz73vz4aV5av5vvv1tEntmtTF3K9SMisldEjojI84O0Ls7Pv0dELl5luX5TRP5dRG4XkRkReVyv5epCpl/Jy8+siHxSRJ7opdW3d+hd93wRURHZ5u3rZ9n6RRH5Wp5fsyLyMi+tvpWt/NiFefm6S0T+utdydZFXf+bl05yIfLORTKpqf8Ef8Ghga/77YcAc8ETgVuBH8/0vAXbkvy8Drsp/PxL4NDCUb/8rsB0Q4CPAsxrIJcBY/nsYuAV4KnAN8MJ8/zuA6fz3y4F35L9fCLwv//1E4DPAQ4HTgc8Dx6yiXKcBZwJXA8/30nk4sC//f3L+++RVlOsZwGj+e9rLr57J1YVMJ3jX/iTwj4PwDr26cRNwM7BtEOQCfhF4WySdfpetxwP/5u4JPLLfZSu49leB3U1ksp5JBFX9qqruzX9/G/gscCqwhazyANwA/Ez++4nATH7+/cA3gW0i8miyxmCPZm/pauB5DeRSVZ3PN4fzPwXOBd6f73+Xd48L8m3y41MiIvn+96rqd1T1PuBe4KzVkktVv6CqtwOLQVLPBG5Q1QdU9RtkeXz+Ksr1CVVdyPffDGzqtVxdyPQt7/Lj83Ohz+8wZwfwp8B/e/sGQa4YfS1bwC8Bb8/v7dqJnsrVMK9eBLyniUymTCqQzDX0ZDItfyeZdQjws8Bj89+fAS4QkWNF5HTgKfmxU4H9XnL7831N5DlGRGaB+8le8ueBb6rqkcg9TgW+BJAfPwg8wt/fJ7mKGCS5XkrWk+y5XHVlEpHLROTzZA33r7UhU125ROTJwGNV9cNBMoPwDn8md1W+X0RcHe23XBPAhIh8SkRuFhHXOPe1bOXXPI6sF/nxJjKZMilBRMaADwC/nluILwEuE5FPk3Xx/yc/dTdZht8GvBn4F+AIWbczpFH4nKp+V1Unyazms4DvK7lH0f37LVcRAyGXiPw8sA14fRty1ZVJVd+uqt8D/A7we23IVEcuERkC/gx4VeR4v9/hh4DTVPVM4GMc7Zn3W65jyVxd55D1AnaJyEm9lqvLevhC4P2q+t18uyuZTJkUICLDZIrk3ar6twCq+jlVPU9Vn0LWJfx8vv+Iqv6Gqk6q6gXAScA9ZApmk5fsJuArvZBPVb8J3EjmEz1JRI6N3GM/ee8pP34i8IC/v09yFdF3uUTkx4BXAz+pqt9pU64u8uq9HHVR9DOvHgZ8P3CjiHwhP+/afBC+r+9QVb/uvbe/IPMS0G+58vv/vaoezt1/d5Mpl0EoWy/kqIvLyVpfpthAykb/I9PMVwNvDva7QbOh/PhL8u1R4Pj8948DN3nX3Jq/UDcA/+wGcp0CnJT/HgH+GfgJ4G9YPsD28vz3ZSwfgL8m/30GywdJ99FskLSWXN5172TlAPx9ZIN+J+e/H76K+fVkMgPh8UE6PZOrC5ke7137XOC2QXqH+f4bOToA31e5gEd71/4UcPOAlK3zgXflv8fJ3EiP6GfZyre3AF8gn3PYJK+6ysj1/gf8EFm37nZgNv97NvBKssiuOWCnewFk0Ul3kw3Uf4xspU2X1jaysZbPA2/zX1oXcp1JFhFye57ma/P9m8mixu7NC85D8/3H5dv35sc3e2m9OpfpbhpEmHUp1w+SWT8PAl8H7vLSekl+/r3AJass18eA//Le+bW9lqsLmd4C3JXL8wngjEF4h8G1N5Irk37LBfxJnl+fyfPrCQNStgR4E/DvwB3kjXs/y1Z+7PeBnZG0astkM+ANwzCMxtiYiWEYhtEYUyaGYRhGY0yZGIZhGI0xZWIYhmE0xpSJYRiG0RhTJsZAIyKP8FY2/U8R+bK3/S+rLEvhqsL58RNy+d7m7bsxXz3XyfzI4Jplq+6KyLCIvEuylYI/KyK/m++PrgibH/vLfL9bQmSsQP7nichrC47Nx/Z3i4h8TERO7mWaxmBjocHGmkFEfh+YV9U39On+zwBuUdUFEZkGzlHVF3jH30I2cewBVX1Fvu9G4LdU9bZIeg8D/gF4CPAKVb1NRH6ObPb9C0VklGxewjnAF8kmxs7nqzN8Enilqt4sIidoviCkiLwJuF9Vd0bu9y952gcix+ZVNaqEuiFftnyTqv5Rr9I0BhvrmRhrFmdNi8g5IvJPInJN/l2GnSLy4tySv0NEvic/7xQR+YCI3Jr/Pb3O/bR4VWFE5CnAo4CP1kgytuquAsfny1+MkK3/9i3NiK0Ii6dIJL9mhYUoIhPAd5wikewbF3vyfNjhnTeW97r25nl3Qb5/h+Tf9cm3/0hEfk1EHi0iN+W9rjtF5IfzU64lW4PK2CCYMjHWC08iW6HgB4BfACZU9SxgF9m3GiCbTf5nqvqDZJ8P2NXgfkurCucLH74R+O2Cc6/KG9vX5A1+2aq77ydbGeCrwH8Ab1DVB/Jrlq0Iq6q3uItE5CrgP4EnAH8ekeHpwF5v+y1AJ8+L//T2/zfwU6q6lez7Lm/MZf5L4GLveV8IvBv4OeB6zRYXfBLZTH00W7r8oSLyiII8MdYZpkyM9cKtmn2H5jtkS3m4HsIdZMvdAPwY8La8Qb4WOCF3NdVCVq4q/HLgOlX9UuT0F6vqDwA/nP/9gpSvunsW8F3gMWRrW71KRDbDyhVhReT73UWqekl+zWeBF6xINfvg29e87adzdHG/v/IfD/hjEbmdbHmZU4FHqeoXgK/nSvA84N9U9etka89dkrsgf0Cz7/847s9lMjYApkyM9cJ3vN+L3vYi2fLfkJX37Zqt7jypqqcGjR+SfX55VkSui91E4qsKbwdeIdkKum8ALhKRnQCq+uX8/7eBvyZTFmWr7v4c2dcUD2v2AaVPkSmuJfToirDnB/u/C7yPox9t8zlEtlbbsksi572YbNznKbni+i/vul1kXzK8hOyzC6jqTcCPAF8G/kpELvLSOi6/r7EBMGVibCQ+CrzCbYjIZHiCql6SK5pnh8dyq/z/kSmS+71rXqyq/1tVTwN+C7haVS+X7GNp4/m1w2QruN6pqgdVdVxVT8uvuTlP8zYy19a5knE8maL5XD7ec1Ke1ghZL+tz+Xnfm+8XspWFPxd59s8C3+ttf4rMVQWZAnGcSDaAfzgPOPAj1j5IpsB+ELg+v+fj8vP/gswVttWT5X+RrUhrbACOrT7FMNYNvwa8PXfhHEv2CeZfqXH964Ex4G/yoY//UNWfLDn/ocD1uSI5hsxt9BcV93g7cBXZqq8CXKWqt4vImcC7ROQYMiPwGlX9cO4ye5eInJCf/xmy79eH3EQ+/qFZCOcrgb/OB9U/4J33buBDInIb2fjHkmJS1f8RkU+QfbnPfUjpHOC3ReQwMA+4nslTyJZ/P4KxIbDQYMPYIOShyx9S1Y91ef0Q2SD+z6rqPQn3ulZVZ7q5l7H2MDeXYWwc/pjsQ261EZEnkn3bYqZKkeTcaYpkY2E9E8MwDKMx1jMxDMMwGmPKxDAMw2iMKRPDMAyjMaZMDMMwjMaYMjEMwzAaY8rEMAzDaMz/B45lQLbgaWtXAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.plot();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This light curve is of the object [K2-133b](https://exoplanetarchive.ipac.caltech.edu/cgi-bin/DisplayOverview/nph-DisplayOverview?objname=K2-133+b&type=CONFIRMED_PLANET), which has a period of 3.0712 days. Let's plot the folded version of it to see what it looks like." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJztvX94XVd55/t9iS3JiqJIPkJKZFtJUZKx7ELsxA1OmUwTUiDQknRIOy1DS0iZ5t4WbpmB9pny0FvukLbDFCilhQduhib8eFpoC5cZl4YxqS8hgco0PyyCY7nYysXHjmILHf04OTmRZMN7/zj73Xr3Omvvs6XzS0d6P8/jxzr757v2Xmu96/2x1iZmhmEYhmGslpc0WwDDMAyjtTFFYhiGYVSFKRLDMAyjKkyRGIZhGFVhisQwDMOoClMkhmEYRlWYIjEMwzCqwhSJYRiGURWmSAzDMIyq2NRsARpBX18fX3nllc0WY0VcuHABmzZtiNcTYmXeGFiZW4cnnnhimplfWum41ivZKrjyyivx+OOPN1uMFTE9PY2+vr5mi9FQrMwbAytz60BEp9IcZ64twzAMoypMkRiGYRhVUTdFQkT3E9EUER2N2U9E9OdEdJKIniKi69S+u4joRPDvLrW9jYjuI6LvE9FxIrqzXvIbhmEY6ainRfIZALcl7H89gKuDf/cA+CQAENFWAO8H8EoANwB4PxH1Bue8D8AUM18DYBeAb9ZFcsMwDCM1dQu2M/MjRHRlwiF3APgclz6IcpiIeojocgA3A3iImWcAgIgeQkkhfQHArwPYGVz/xwCm6yW/YRiGkY5mZm1tA3Ba/T4TbPNuJ6Ke4Pe9RHQzgAkA72Tmc76LE9E9KFk62L59O6anW0vnzM/PN1uEhmNl3hhYmdcfzVQk5NnGCds3AdgO4NvM/G4iejeADwP4Nd/Fmfk+APcBwL59+7gVU+9aUeZqsTJvDKzM64tmZm2dAbBD/d4OYDJhew5AEcBXgu1/B+A6GIZhGE2lmYrkAIC3Btlb+wHMM/NzAA4CeC0R9QZB9tcCOBjEUv4epRgKANwK4FgT5DYaQDabRTabbbYYhmGkoG6uLSL6Akqdfh8RnUEpE2szADDzpwA8COANAE6iZGncHeybIaJ7ATwWXOoDEngH8J8BfJ6I/gzAD+UcwzAMo3nUM2vrzRX2M4B3xOy7H8D9nu2nAPybmghoGIZh1ASb2W4YhmFUhSkSwzAMoypMkRiGYRhVsSGWkTdaj6GhoWaLYBhGSswiMQzDMKrCFIlhGIZRFaZIDMMwjKowRWIYhmFUhSkSwzAMoypMkRiGYRhVYYrEMAzDqApTJIZhGEZVmCIxDMMwqsIUiWEYhlEVpkgMwzCMqjBFYhiGYVSFKRLDMAyjKkyRGIZhGFVhisQwDMOoClMkhmEYRlWYIjEMwzCqwhSJYRiGURWmSAzDMIyqMEViGIZhVIUpEsMwDKMqTJEYhmEYVWGKxDAMw6gKUySGYRhGVZgiMQzDMKrCFIlhGIZRFaZIDMMwjKowRWIYhmFUhSkSwzAMoypMkRiGYRhVYYrEMAzDqApTJIZhGEZVmCIxDMMwqsIUiWEYhlEVpkgMwzCMqqibIiGi+4loioiOxuwnIvpzIjpJRE8R0XVq311EdCL4d5fn3ANx1zUMwzAaSz0tks8AuC1h/+sBXB38uwfAJwGAiLYCeD+AVwK4AcD7iahXTiKiNwEo1EdkwzAMY6XUTZEw8yMAZhIOuQPA57jEYQA9RHQ5gNcBeIiZZ5h5FsBDCBQSEXUBeDeAP6yX3IZhGMbK2NTEe28DcFr9PhNsi9sOAPcC+AiAYqWLE9E9KFk62L59O6anp2sgcuOYn59vtggNx8q8MbAyrz+aqUjIs43jthPRHgBXMfN/IqIrK12cme8DcB8A7Nu3j/v6+qoQtTm0oszVYmXeGFiZ1xfNzNo6A2CH+r0dwGTC9hsBXE9EPwDwLQDXENHDDZHUMAzDiKWZiuQAgLcG2Vv7Acwz83MADgJ4LRH1BkH21wI4yMyfZOZBZr4SwL8G8H1mvrlZwhuGYRgl6ubaIqIvALgZQB8RnUEpE2szADDzpwA8COANAE6iFPO4O9g3Q0T3AngsuNQHmDkpaG8YhmE0kbopEmZ+c4X9DOAdMfvuB3B/wrk/APCT1chnGIZh1Aab2W4YhmFUhSkSwzAMoypMkRiGYRhVYYrEMAzDqApTJIZhGEZVmCIxDMMwqsIUiWEYhlEVpkgMwzCMqjBFYhiGYVRFRUVCRP2ebf+qPuIYhmEYrUYai+RRIvp38oOI3gPgK/UTyTAMw2gl0qy1dTOA+4jolwAMABhH6RO4hmEYhlHZIgmWdv9fKH0P5EqUPo9r30w3DMMwAKSLkTwE4JUorbb7BgAfJaIP11swY2ORzWbx7LPPNlsMwzBWQZoYySeY+a3MPMfMRwH8NID1/QFiwzAMIzUVYyTM/D+c3xcA3Fs3iTYQ2WwWADA0NNRkSQzDqIS113gqKhIieh4ABz/bUPrKYYGZL62nYIZhGEZrkMYiuUT/JqJfgGVtGYZhGAErntkeuLpeXQdZjHVENpsNXQGGYaxv0ri23qR+vgTAPiy7ugyjJgwNDaGzs7PZYhhGBIuLpCPNhMQ3qr8vAPgBgDvqIo1hGIbRcqSJkdzdCEE2IjbKMYzWwdprPLGKhIj+AgkuLGb+7bpIZBiGYbQUSRbJ4w2TwjAMw2hZkhTJXwWTDw1jxZgbwDA2Dknpv/8sfwRuLsMwjA3F0NCQDYpSkKRISP39qnoLYhiGYbQmSYrE5ooYhmEYFUmKkewkoqdQskyGg78R/GZmfkXdpTMMwzDWPEmKZKRhUhiGYRgtS6wiYeZTjRTEMAyj2diSKKtjxYs2GoZhGIbGFIlhGIZRFaZIDMMwjKpIWmvre0hea8uytgyjCrLZLGZnZ9HX19dsUQyjKpKytn4++P8dwf+fD/5/C4Bi3SQyDMNoEmmC7BaQL6di1hYRvYqZ9cz23yOibwP4QL2FMwzDMNY+aWIkFxPRv5YfRPTTAC6un0iGYRhGK5HmC4lvB3A/EV2KUsxkHsCv11UqwzAMo2WoaJEw8xPMfC2AVwDYw8x7mPnJSucR0f1ENEVER2P2ExH9ORGdJKKniOg6te8uIjoR/Lsr2NZJRP9ARMeJ6Gki+mD6YhqGYRj1oqIiIaIBIvpLAH/DzPNEtIuI3p7i2p8BcFvC/tcDuDr4dw+ATwb32wrg/QBeCeAGAO8not7gnA8z804AewG8iohen0IOw1iTDA0NYdu2bc0Ww0ApgC5BdGPlpImRfAbAQQCDwe/vA/iPlU5i5kcAzCQccgeAz3GJwwB6iOhyAK8D8BAzzzDzLICHANzGzEVm/kZw7SUATwLYnkJ+wzDWKK3YgSd9o6QVy1ML0sRI+pj5b4novQDAzBeI6Ec1uPc2AKfV7zPBtrjtIUTUA+CNAD4Wd3EiugclSwfbt2/H9PR0DURuHPPz880WoeFYmTcGusyzs7MAgM7OzmaJU1M54q6z3t9zGkXyAhFlEExOJKL9KAXcq4U82zhhO4L7bwLwBQB/zszPxF2cme8DcB8A7Nu3j1tx0lcrylwtVuaNgZS5WCxGfjeLWsmRdJ1ml7GepFEk7wFwAKVvknwbwEsB/FIN7n0GwA71ezuAyWD7zc72h9Xv+wCcYOY/q4EMhmEYRpVUVCTM/AQR/QyAf4WStfAvzHy+Bvc+AOCdRPRFlALr88z8HBEdBPDHKsD+WgDvBQAi+kMAlwL4DzW4v2EYBgCbpV4tFRUJEU0A+BAzf0pt+yoz/3zCaSCiL6BkWfQR0RmUMrE2A0BwrQcBvAHASZSWXLk72DdDRPcCeCy41AeCbdsBvA/AcQBPEhEAfJyZP52+uIZhNBJbTmRjkMa1dR7ALUT0SgD/W5AxVTFnkZnfXGE/Y3kdL3ff/QDud7adgT9+YhiGsSbYqAozTfpvkZl/GcA4gEeJ6AokrApsGMb6pB6prUmptEbrkMYiIQBg5j8hoidQmlOyta5SGUYFzGViGGuHNIrkD+QPZj5ERK8DcFf9RDIMwzBaiaQPW+1k5uMAntXrYAV8tb5iGYZhNB+zfNORZJG8B8BvAPiIZx8DeHVdJDIMD9agWxN7XxuDpA9b/Ubw/y2NE8cwDMNoNZJcW29KOpGZ/5/ai2MYxlrFrIt4NrrFnOTaemPCPgZgisRoGhu1wRqtRzabxezs7MZca4uZ726kIIZhGEZrkib9F0T0cwB2A+iQbcz8gXoJZRhGbdnorpc4Kj0Xe17pSLPW1qcAdAK4BcCnAfwigH+us1yGEcEadGtgCmtjkmaJlJ9m5rcCmGXm/wLgRkSXfzeaxEb4Glsjy5h0r43wrA1jtaRxbb0Y/F8kokEAOQA/UT+RDMNYS2xEK2OlZa7Hs2ml555GkXw1+LTth1D6Tjqj5OIyPLTSyzfqR63rQb3q1Vqrr2tNnlowNDTU9E8J15s0H7a6N/jzy0T0VQAdzLy+P0BcA9Zjg2gkjXQjrQeXldW31qXSu2uFd5sm2H4RgJ8DcKUcT0Rg5j+tr2itz+TkJID6VoBG3KORrIdOfS2yHupHPTrU9fBcgOYrmzSurb8HsADgewB+XF9xDCMdcQ2nHg0qm81icnISg4ODVV8HqL1sa4lmdGSN6ER1HdD3Sbr3auVaa+80DWkUyXZmfkXdJTFWTL0aTrNHN5q0MtSi8SXdy+1AVku1FmS172YtvVuXOIVtbs5kajXQqYY0iuRrRPRaZv563aVZJ0hnsV6oZ6DXN8prBSpZRHHHTU1Nec9b7f1Xymrq5lp5N41046a5R706cCnn/v37w/ukpVkDhTSK5DCArxDRS1D6fjuh9Mn17rpKtsaJe2Hyu1WUSRrTfCXn1JtWGTHKs6n0DJOyefRzrnS9NJ1sNpvF1NQU+vv7y+RcLa48rpxpO2RBBhUrqZfrLU7YiqRRJB9BaRLi95h5w3+rPW1HVqtRSjWddtpzm90Q5f66U6000pfz3BHhasqyHjqgOHfQWirb4cOHAQBbtmxBsVgMt8s7kzKsdqRfi5iE3Hel11qt4lzJPSrhPsdGkkaRnABw1JRIMr6KVE3nPzs7G2lslSygSsfF3cdHXIVMq0Sr6cSa2fGtxc53bGwMk5OToZvDtXQrdRr1LFMjffO16CTdZ+Grz+J61ErN3VdN21oJWl5ZWSHu3iJvszwhaRTJcwAeJqKvAViUjZb+m0yzR/lpmZycDBsIsPbkjXNjxHUo9RjF+pS1O3JNw9DQ0KobutxHvysgfXlF5lwut6LzVoNPwWi32rlz5/Diiy+WyeC+26RnpY9dTVlcGSvda3BwEHv27Cnb57ri0sb9kupxMwcHqyWNIvn/gn9twT8jJVqZxI2G3FTCuBFerRWTq0D0PYBShyUNv9GmcpqRrnTKU1NTkUYr59Y7iO+OXAX3Pbm/V/Ms9b0kvjE1NYVcLoexsbGK70i/60wmE25bax2WlkeXuZZxR7nG5OQkxsfHw3vI/VZS31fi9nI9DHHHpcF9DpXiZ40gUZEEkxG7mPl3GyTPmiatKS8vdmxsLPw9NjaGXC6HW2+9tSYNs9pGPj4+jrm5OQwPD4fbRHm4uBW1VkrN7SCeffbZWN95pU7EbUSuG8J1E/ju715LH+vKLB2625G7dcMdlcdd14dWlu570QFzzYEDByL75f5iiYyMjITb086HcImLV+nzfW4gkaGjoyNy3vj4ODKZTNmzy+VymJycrJsim5ubC+XM5XKhkk0ijZVUab/b9uLc1XHWuCD9i7u/UrJCPUhUJMz8IyK6riGSrEGSGrvvBenRjm742Ww2bERJ99AV5dy5czh16lSqEZIoOGG1o163MYnM7nPwBUVrUWHjXB4iH4Ay94J0NtpF4T57H6sZ3aZ1tQDLgWX9e2pqKpRfntf09LT3enK+xqfkBbneoUOHYs+TzjqpHHGWsrtNcC0GOTau3uZyORSLRezYsSOULZPJxCpGF3nXcm3XEk1T93UdyeVyyOVyEbnjztFIJ+5zd1UiTmH5rNmpqalwe9L7bzZpXFtjRHQAwN8BeEE2btRvtq+m05ycnMTc3Bx6enpWfb+4xi8NSJvm+lg3WOvroFxEgWg3CICIq8s3il9t4HVwcBBjY2N4+umnccstt4Tl6e/vx/j4eCRl1b2vjLBdy0CP4l2rQK7hjuikXHI//dzjyiUdEYDIqDrO9aXlAJLTf937aOI6o7m5uXC/nKM7yjg5NFJnqs2Gc88ZGRnB1NQUzp49CyD6bNznVEmx+OSJa5uuUtfodinPzvfOdfvyySLKLa5d5PN5LCwsRM6X+lqpTboKRCtBjdy/GS6uNIpkK0pLx79abbNvtq+Aqakp5PN5APCa6nHWRFxj0p22dHxSqeI6XBfdEck5bqfoojtq8S/rjjyXy4Wd8P79+1ec5lgoFEI3h4vbkN1y63OmpqYwMTER/k5SAr6y6gaq7yH31fukI8pms5ibm1vRCFU6mJ07d5Y9Gynv7bffjmw2G1oacc9GnnU+n0d3d/VTvJIUoTthzrfPt72aEbV2Fyd16oKbbSUDAy1jmoGddtX19/djamoKo6Oj4fnyXuJiZcKJEyewtLQEAGX1SQ+U5F3KfV2XZBqakQacZvXfDf/tdl+eO5BudOZ2RHK+HvFpt5Lv5euMF3GRxDUit5H5XGrio+3p6Skbxce54ASp+HL+1NRU5HpxZUy6nsj5wguhwRsrR1yWzcTEBHK53Ioaj34Gck0pCwA88MAD3pRL1+qohA7sasbHx1EsFrFz587IdnGFZjKZ8DmKApHOTP5271MoFCKKZGJiIlQumUymKhekT0nouJko8J6entD6EPQApFAohO9LK0bXPSl1U99DzpGyu+9HW5y6LvgC69pqA0pKeHR0NJRfWz26nsrzlHqi631c/czn86EikXonlngul4u8S+0BkHvoNi//y2DJHdyInCJ7I0iz+u92AH8B4FUoWSLfAvAuZj5TZ9nWJJV84RpdSbu7u8sqXFLQdXp6GgsLCwCWK4cvjiH4XBpamfjSR13Xld4mI+vJyUmMjo4in89jaGgovH4+n0dPT0+koxKqGXnOzc0hk8kkWgnS+GSUFmcx5PP5sLHp8kijFCtRyy3PEVhW4Lox6mw2HzpWkc1mw85QFJRrkebz+cSAsjuS1+9Cl0uTz+fL6p4EdnV99Sn7sbExHDlyBACwd+9er2Wp0e5Buefc3FzYcYurVZehq6sr0hZ0SrKrRDTyfty6kc1mU2WvSYaffndSP6T+6joRl3yi0Z4GKYNG18dCoRD+PTc3F9ZjN+gvSDtzt8s7l2eYxkKrN2lcWw8A+GsAvxT8/tVg22vqJdRaJSkN0R116pGTHuWMj4/Hum+Acr/9apDOsL+/P/xbTHDpzKRB6QwSbV7HXVd3tECpIU1NTYWNxHUxAZXjBIODgxF3gdDf3x9xt4ns+vrSEfT09CCTyYTKtqenJxyhy7VkNOp2OEkWRmdnZ2zsQ1/j2LFj3jLGlVvKce7cOYyOjpb59/U71OTz+UiZhoaGIj72Xbt2Rc7Xxw8PD0cGJcByR3rnnXcCWO6wxGoQWXVmlS63dITDw8PIZDJl9UNbq7puzM3NIZ/PI5PJxLpt3NifvmfcNrmnMDg4GMqYyWTKFJWOj0inLXVarisdvqAHTYVCAYVCAQcPHizLWNPod+bK7Q5odALA1NRUWRxN6o47aNQDTXnWSZMYa0kaRfJSZn5A/f4MEf3Hegm0lojLzJKRTTabDRuYrhzSYeVyuXDELm4gIBrU0yMk7bbo6+sLXR5uxdSdq6AbqtxPyOfzkQ6lWCyGlVcavxsf0QFE15esr9Xd3R0ZaSXhpr+6Cvniiy/2WlTAcqBWnqvbYclzkUYk2UGdnZ1eX7i8W7cRu/7zYrEYcRNoSwNY7ohnZmZQLBZx+PBhDA4OlnWOIr8MJABEyuAGj7Vi1+UtFAro6uoKlaZraUodkOOl05Py6NiR3Lerqyv8LfWhUCiE1pJ2B/lG3mIBjYyMIJfL4dixYxGrSZBydHd3Y2lpKVJvdB2Xuui6nmSAoBWS7lTlOu5z9c1O1/LLuUePHsXi4iLa29sBlN6FW89E4fgGldqy0PfS9VhvE1lFaYsyljJrhaU9DoLPfS2DPV3n14oimSaiXwXwheD3m1EKvm8IDhw4EHGjSEdWKZit3UlSIcRlIGjT2W0wxWIR3d3d4YgKKO8kgPLgvKushoeHMTExgcnJSRSLRVx11VVlnacPuc7+/fsjVowoEel85FpdXV1h55fkDtAp0qKEdePWHZ1Yckn5/T09Pd7GAyxnREkgXD8nYLmjEJnk3oODg2F8orOzM3RDAMvPVWQW+eQ5SJbZnj17Iu9Xx16ENO9B7imdlzx/iUVoZJs8r+HhYe+I14dOXhA3mLay5b1KeeTdZLPZUCFIxyhzgXz1dXJyElu3bkVHR0dESWt8lou+hlhAbpaUey0pT3d3NyYnJ3HkyBHs3bsXuVwutNCljY6MjODYsWN47rnnIvfxBeYLhUJY32SwIm1cK33xQgDACy+8gM2bN3vldF19grxraeNiVelyu1a7z2XYCNIokl8H8HEAH0UpRvJPwbYNg/b5aqThuKMFSfcVpHFrd4MPXQnb2qKLCIjVU6lj6O/vD+MxejKaVDytCKRzE4sFWB69SKeazWYjJr90GkNDQxgZGYlYGeK7HR8fDzs1N8tKPzuRTa7d19cXKj65l9zbNzrXlqDuqMUq0YoPKCn0fD5f9i5Flnw+H7HWCoVC2Wxkaag6VpXL5com3ungrA+5x8DAQKrsIem8XDl0jMONn2mlLJ2Q3Nd1z7jPT56Fz3+fyWTKZJ6bmwuVysLCAhYWFsKBU6FQiCi4QqGApaUlLCwshPVHrGLXksvn88hms+jo6EBnZycGBwcjik7XKd/zkWvMzMxEyiDHi3KUc9vb29HR0RGx3k+ePBm6EXt6esKySJBc6pXcV5dV5BsYGMDS0lL4TlzLQ+6v437d3d2R2I12j7lxOl9/s6YUCTNnAdzeAFnWLFIp3Q5Id+pSQdyRomTLAMu+VxnBSMOXDnhsbMzrstEd5sjISDhSckf+MsKXUbqMKKWxLiwshBVQV16poNrFIR3hnj17UCgUQteNNGag1Ei6u7tDpQIgDMgL2h2lG5eOewg6SCwNSP6XeIhWorqh+Ea/gsgu19Hk83kcOXIkbLRyv/7+ftx0002h4pqYmAhjTq47Rf4WK7S7u7ssA84nX3d3d9mAQZ6hdiFJvMe1BuPcexqxSvSkWPHr6/ct5dLb5PnodyfWndRxKaskXBQKhdA9JM9SlKB+Xt3d3ZEOW2SQ+iqKuVAohEpEj9z1HBnXGtUdubiOtm7diq6ursiM8rm5uTChRZBZ98ViEYODg5HBhLxbSTLp6emJlFG8CO67keOmp6fD56Trijxn7f6W+Ih4EHzuYzcVX56B9Cm1iLemJVaRENEfJJzHzHxvHeRZ07hBQ+2H9o303EwtF1nRUwJ/YnFks9lIByP7JbagMzncdEh3DoSLKBORTVwRnZ2dKBQKOHbsGAYHB8OR+eTkJHbt2hWONqUxyzPQf8tENl0OYNnvrk1+1+0nuM9Wyir+d9f3Lm476XBFHhm1yah4YWEBhUIhUk4JTEvKrDwPSQMVWXVW2kpcB9od6isbAPT29pZ13i6ZTCY8ZteuXd6EBn0PXyadvr90Ul1dXejq6gotv5mZmfA+cf5514Lt7u6OBIO7urpw6aWXRu7b1dUVyaADEE7Qk2QGue7w8HBYVyYnJ/Hcc8+hvb09vIdYPq48omxdxKXb2dkZDmp0ndXLtchxQGnwoa0+IJp1FecV0NYssOwmbW9vx9LSEiYnJ8PBpVaA+Xwex44dQ1dXV6iQ4yar+tzrOmlA3ldchmc9SLJIXvBsuxjA2wFkAGw4RTI5ubwYoLZGtFvI56+XTl6PeH24gUKd4SXBfT0hS3A7NKmketQkSkgsEt2Z9/f3h1bH1q1bASwrnNHR0cQJTq4rQnznUl5dkbX5726T0blWQMLJkyfDv3VjFhdWXKPes2cPcrlc+Lykk5TRozxrUZpablEqWk793HVdAJY720KhEI7Qdbq0flailLX1k7RGl0a709xjpaPWz9kXuNV/i1I+duwYFhcXI9ceGhoKs9G0O0l3uPIsffL6Bleajo4OdHV1hW3JtTpdpSTvRb9vGa3LcykWi6GimJiYCGXt6uoK/3Y7V7FKpExyrDxDqX/iKtXtXqxxsaTiytrX14dz586VuUq1cpfB0NzcXKhQxNWrB3HanSttwE2ikGs3KiU4VpEw80fkbyK6BMC7ANwN4Isofexq3SPuAN2hyMsUk9kdTYp7QFwKOnVUznODoL6Au+sT16MPaRjAcuBRRuBSeR599FEUi8WwImkzXDo46SREPqnk+Xw+HKlNTk5iZmYGHR0d4SjXHWn7LJ+5uTkcPXo09C+7LhnXnZbP59HW1hY7CtPmvZ6cJxPcRGmJlSbzF/Qz1HLOz8+Hz1+sFWC5w3SVc6FQwKOPPorBwUHceOON+PKXvxw+V3dkqXHdcSMjI5HOZnZ2NswCjFvxQM/K1skM8hwOHToUiVtIh+vzk7udqM5wk0wqPeAoFothZw8svy+pK52dnejs7AzrydatWyNKRuSRujM4OIhTp05haWkprJuSyTU1NRW2Hzmns7MTxWIRMzMzkSQIUTo6yC3yiPttcnKybFmSrq6u0Irt6upCZ2cnFhYWkM/nQ2svl8uFz0FbXPo9Sn3WMSTp+KVty3vI5XJhJpibRZhk1YpbUVvFPgtDPAD6fbp9T715SdJOItpKRH8I4CmUlM51zPyfmXntrh5WY2699VbceOON6OnpCSu4IG4B138pv/XyIXItIGpyuiMY1zcORFMb3QCdbzQuabu6U5aORisoaaC+a0gHUSwW0dHR4R3Z6NG+Sz6fx+LiIubn5yPHSAchHRRQUmx79+4tGykPDQ15XSe6nLJNlKlYhNqNIee7I095T+I/d0dwki00NDQUua92Cwly36GhoUjjFQUunYC816GhocjKyz7EItTgStOfAAAgAElEQVSxOa1EpMMQ61PmZPgy3dzYgqDL68ZrdHxKYhxiuUr90M9hcXExfK9a4Ug8T1xXAwMDZW0pDunsOzo6Ioqys7MzktYuo3gAoatO6i5Qqk9XXXVVxOUk70qOkfcv70mu09XVFVH03d3dYdt3n6Gc77NqxZpzk3J0DFUUoB4cify++pLUBvUApN4kxUg+BOBNAO4D8HJmTjdZYPn8+wH8PIApZv5Jz34C8DEAbwBQBPA2Zn4y2HcXgN8PDv1DZv5ssP16AJ8BsAXAgyjNsG/YlxtdpaHnUOhGIWmNlcxKPWIBln3h3d3d6O3tLXM1iNtHOk6dTqwn5/lmbQMlC0tMf58suiw67iPbpOFWmu0rrrX29vawDFIu1zUijUWslsXFxVBJDg8PY2RkBAcPHox1oWjrJi4VulgsYmFhIQzazszMRJZjETm0i0lPdBTFtLCwEBmdy/0FeT5xLkzxa/vmBLhIwF1n+egFM8X6dGXQ9/IhCkhbwzpjSSyN/v5+DA8PR7L99DVnZmawuLiIyy+/PNx2+eWXl807kmtr66pQKODcuXNhTEuIs0Z9CQCCO4hzXUdibfqyJmXeRpzideNM+jlrVxOwXOdlsOfO5XLL4ab766C79gzItYFoCrZYuqIotPUprq7x8fHYhSxrTVKM5D0ofRHx9wG8r9TvAwAIpWB7peHEZ1BKG/5czP7XA7g6+PdKAJ8E8Eoi2grg/QD2oZRu/AQRHWDm2eCYewAcRkmR3AbgaxXkqAqdSy+N2J3cB0Rzxufm5jAyMhKZsb1nzx7vqCKuUzl16lSYsqizaSSvXKcdPvrooxF3gLh7gNIyF0B0NrDcUyqmboxyT+lQ5feuXbvKUhA1MjLW80REicjI6NixY+H1AJSN8mdnZ7G0tFTWgHXwURqnTseUQK5OJdbKrr+/H9lsFvPz81hYWMDzzz+PzZs3o6urK2xkkjHkU5Lyzo4ePZo4+VLLrIOfcQkX0qn6kMmu0iGJ9eFOaNNKQQ9efOmgcoxPAYmFMTMzE1puuVwuUu9k1C+KWZC6cvLkydDqAFAW9Jf4kSBySCxPrqPPB5Y7V0lOEbebnjCo37mWTawNuZarlPRkYD0/Q1xo8vxkm9RzbVHp1RT089fPd3Z2NqzHkiXoWq56RWhxocrz1fN8XHyJNmnnKNWKpBhJoturEsz8CBFdmXDIHQA+F1gUh4moh4guB3AzgIeYeQYAiOghALcR0cMAupl5NNj+OQC/gDorEo1edkNG1dIZiDnsm28Sh7hjRkdHw8YkHUtfX1+4yJsg+6Riy4hOV3g9I9sd4WnfOYBwwp0EEcUlobPNOjo6sLCwEDYKtyHqNF43TnDppZeGjVqsoa1bt4auEekIdAMQi2VoaCiyjIU78Uw61bm5uXCSmQ7S6wC3uGSee+45LC4u4pJLLglTT8XykYbvurfcGez6WerZ6FqJ6gC8O6FtfHw81WKW+j6+WIe2QOU42S6z6HV8xmf9aCtDd4xuLEvqhZuG2t7eHm536e/vx8TEhHcmuUYUoU7JlfMlBiaxDpHFTV8XtAvTTesFEHbMBw8eLHt+QHSlX4mTxKXzyv9uCrhbzomJCczNzaGtrS3SPl3Z3aC5O+B0lUhSNmgjXVpCmgmJ9WIbgNPq95lgW9L2M57tXojoHpSsF2zfvj32A0JJHD9+HNPT0+H3E66//nqcPXsW586dw4ULF3D+/HksLS3h/PnzKBQKOH78OIrFItra2vDMM8+EKY5AaUQiqadASVEAiBwj5/b19WFhYQFtbW1hhyWumIGBAWzbtg1btmwJrSFRONqs37x5cziDeHp6OhxFFotFtLe3o6+vD319fWhra0NbWxuWlpbCDmRhYQFnz55FoVDAhQsXcMkll2BgYCB8hk8//TROnz6N3t5eFItF9Pb2oq2tDe3t7RgYGMCpU6cwOzuLjo4OXLhwAbOzs6FMbW1t4bOT8ra3t4cWmHRKZ8+eRXt7O6666iq88Y1vDDOr5NhCoRC+F/l7dnYWAwMDKBQKOHnyJLq6usIgp6S3btmyBRdddFEog3RQ8n9bWxs6OjrCkazI09bWhk2bNoVyDwwMhOe0t7dj9+7dWFhYwPHjx3H+/HkMDAygra0N58+fD8s4OzuLfD4fbvvGN76BtrY2vPzlL8e3vvUt7N69G+fOnQvfoZTv9OlSc+jt7cXTTz+N7u5uXH311ZG6tbS0hIGBgTApQsrS19eH3t5eFAoFnDp1Ct3d3Whvbw/rgK5zvb29mJ6eDmMlCwsL2Lp1a8TNuXXrVrS1tYWLiko9PnHiBC5cuBAG7OWZbt++PXL+6dOn8cMf/jCc+Cf37+jowLZt28LtW7ZsCbOc9LOXetrW1hZaOyK3PHe515YtW/Diiy/iwoULmJ6exjPPPBM+L5FP3vOWLVvw3e9+F+fOnQvrUz6fx9atW0O5pB5Lfdi2rdT9TE9Ph3IfP34cJ06cCN/X9PQ0XnjhBfT09OCaa64BgLD+6jiSXH96ehpdXV3he5Q6ePr06TKFfcUVV6Cvry8sl8jwzDPP4OzZs+ju7g77nXrTTEVCnm28iu1emPk+lOI72LdvH0uFXwk7d+6MuGhkEtX58+exadMmXHHFFaFFIL7966+/PhyF7dy5EwsLC5iYmMCpU6cicQ+59tTUVBhIzGQyWFpawo4dO3D69OnQnaOtE7EQnnzySSwtLaGnpydUOLrCXHHFFbjxxhvDpaoXFxcjqY/d3d3YsWNHaKq3tbVh9+7dAJbdVB0dHeG1pfHKCE0CrYuLi5H/+/v7ceLECZw/fx779u0riyEAwKZNm9Db2xve/9y5c1haWgpTj7u6unDZZZeFwduhoaHwi3r9/f3h8dLYFhcXw/ciSiCfz2P37t3hxKwLFy6gvb0dW7ZswaZNmzA4OBiWRc5fWlrC0tISnn32WQDLcYmpqSksLS1h586d4UhQGvvS0hJe9rKXASgp9E2bNmHz5s1YXFzE0tJSWM49e/aEbtKOjg7s2LEDExMTaGtrC+tCb28vnnzyydC6k05O3rPcUya0Sr257LLLQpmkTAsLC6FSkjon5ZM4z/T0dJhppJ+tOxdnYmIifKZ6Qu2mTZvCOnfttdcCKI20X/ayl2FwcDAST9J185prrkFXV1c4j0Qyp/r7+yNL6vf29mJhYQFLS0uRgZc8P7m2POvNmzejo6MDO3fuDOva888/j8XFRQwODmJxcRGZTAbXXntt+Izlvtdeey2y2SxefPHFsL7L85fY2OLiYiTdXayPQ4cOheWWutLT04PFxUXMzc2hs7MT58+fD5/t0tIS2tracNlll0Us4NHR0bBeSxabvJfOzs6Ih0LqfXd3d6S+SnxOlFRvby9W0/etlGYqkjMAdqjf2wFMBttvdrY/HGzf7jm+YUjmiUYH8Xzo1DxgefE+d9apuGXE3SWjUEGyRNxl0oHkiVLyDQctj/5b3AduMFXQbi53RrSY46J4tJmvOw43e0Uairij9LpVcs+kAHKSWe/ukwmSQ0NDYccwPz+PQqEQCVTqiXz6Ger0aVcW/d0Icbf40M9Vu9/y+Xw4SS0Ondqt3Yjj4+MR14oOxIrMYpnIcUnxGgBlgXL9XmTui47NSaabuGNuvPHGSAzEnT8jy7AMDAyU1TX57a6f5SY96Pqo09fjXLl6wmE2mw3fNRCdZyHf+ZFnKeWR5+BOAwCWFx/11VXZJgOutra2UGFLWrFeiFTuKYk0etZ6mqVO9CREec+yrdnB9npzAMA7ieiLKAXb55n5OSI6COCPiag3OO61AN7LzDNE9DwR7QfwHQBvRek7KXVDB9oleK6XmQD83wBIOwlIXrp0XNIZDw4O4oknniibSyJ5+jJJyU2XzWQyYSaXr1MbGRnB/v37I41G0KmGcTOjZSZ43JLpErsAlmf06tnsuoOO6zxFSUnnJxP/dEBRritJCHNzczh58mSZX1xkkSVoxFWjl3nxoZWlG/PS5XEDru5xbu6/toJ8K8Lq6wPLMQxR+Hp+gm9iol6eJ5fLhbEb6TB1Vp/cSweI9QTKuMw0sWC0/DrWoN+/nOt+NkFccbqzk2vqZd7lOQPRQYvOYNOKT8/fknPj2qKOO8Wltss+vT9p4qVe80zev9QRcRe6sRd39rlkD0rWnlb++nnt3bu3TG6dRCHPq+kTEquFiL6AkmXRR0RnUMrE2gwAzPwplLKu3gDgJErpv3cH+2aI6F4AjwWX+oAE3gH8JpbTf7+GOgfa3Ylf+m8J0vo+vqMnA8nnSN2go14kT6fdykj36quvDleelc5D5jsA0cUX9bLl7keYtNwyO17LIHMP5FjpdPXo15297qZYijzA8ryJ4eFh7NmzJww26gmRUl535H/11VeH7hhZCkM6N92ZyjPSclx11VUReXQHIwkSHR0duOmmm8pG+FJusc7ESnInnPqSFfSzctHp3ZUC7HHWYCaTCTtW9z0kdRJxqcV6Ap1kIMr8JgARF5yUsVAohOuOCVo56MmGQlya8/DwcJjm7bMMBG0FiEJwrRX3WUiChpRTZ3gJroKstISIO0CQ96g/hSDvWa+RphNPenp6wlhdJpOJLHOjyy6K4uDBg+EKEb7UZElQ8dUZGWw0SoEIdVMkzPzmCvsZwDti9t0P4H7P9scBlM1JqSc6jU+vhwVEsyPiFkjLZrNhCrD7cnWjkGDy8PAwxsfHcfbs2chsco1bgSR10reYoJsRpRu77qx9HxeSUb9ujLlcLvRtx6HXKBKZ7rzzThw6dChcywso/66JnOtOtNTPW7uGtJvFncmtM5dyueVP8Golol0u8oU9SfuUzB+Z1Q8sr2Ir9xAlpzNtfK5PTaVJiGmQ+ihKXytF37t3R/duppauh/K+tfIvFAoRF5M7SpfnoL/RIYoobqUA3YHGWQR61QY9H8K3tLs761vv27t3b2hlawvIp8D0c9RzkXSZJyejX7T0dfayXfYlpY3rDDs3A1HKIe/60KFDkUGbSyVXe71opmtrzePGN9yP/ADLXyATRSIuBJ3GCJTPbdDohfLknt3d3WFwLW75DS2DLNnh4qtUeskWF62MZLKWbmxzc3PYunVr+O0Q7X+X/WKNuB2U67bRz6Onpwe7d+/GqVOnIpaaPkfHIbSv3CWTyWD//v04cOBApCyuewxAZFkZn8UggWpdfpFXztMTQ2XejnQMuv7oj5jJsRJ3SBpB6tG4tojjBi++eQX6Y2tCnFvORS+tr3FXn9X4OtZMprQq7ezsLF588cXIPl9sJQ3aYteDBD0Y0atBaPeQz72oY05uW9eINaJn/QPRr21qN1ycIvGtkSVp19r9NTY2hj179pS1tbjJlO5gqd6YIklAj8w0UgGlkmmTWvy2sk5Xf39/2fnuwnRAubtIVoWNW4zNrWy+hi4NzPVpy4hRd9JuoNHFnfENRDtTXzBSl89Fd2KSTfb000+XddxJoy8gankkzRjX5dL73eeiJ2h2dHSE8uh0YHcJHCmPuMa0vK4sYt3J85RJeiKH/ra8uC9kzoxWONlsNhITkY5TT0p07+2bH+OuiJCEWK6i1Fw3k8Te4lxJLr5OV8riznqX++l6qpWIO4lRxxVkUUX3Oekya1eu3q7btcjhIgMn/eloLefg4CC2bNmC3t7eWNecLq+e3a4Vm06c0e9dx4+aRVWTDjcSMrkqTQaEpOHp3wAiisWNF4irSJSXzlkHliu4fEBK3FluB6tnz1dSDgAi5r5UeneELIvpAfDGhQYHB8OAt67s2h2hR3Y6M6qzsxO7du0q6/zE5y6uFdfU17LpxiX79fF6RriUM5fLYXR0FFNT0c/f9vf3Y9euXRgaGkJnZ2eYmy9IamVPTw/2798fyqjl0u5Q7cuWAYc8z8HBQezevTsiC1D+wSYX7dKrlDGoFZMP/ezkWcjvTCYTDkL0d2ZEPvmXFNTVz0K7gvVaV1pefW+xYHyWoqyjBiASE3An/ElsUGdruej9cr6u88By7HF0dDSyIrYu19jYWMQzUSn+opFjjxw5Enn3Ervz1Qdtde3ZswcjIyMYGRkJY15imTcCs0gScF1bwPLSFUDyWklA+ZcVdaATiLq0ku4v99GZTFKJRGnIyExG6HoykztKFNl9oxg9O7uSy8NFZ0ONjY2FVkp3d3dohcVlfskSKXr9ItcaFJlck97XSUiDljKKktTb5RkDy59BlpnwYvHJvA0gmhXkC4Dq9a/cfRKPkucjHens7Gx4vyNHjoTLZwg6C0k/W22pyfn6vP3795ctxx+Hnv3uugCB8vrppiLL4EKera47cfEDvSqA/qCba83EJSloC1gGHUlp1D7LQ5SMb792J05NTZWlbfvalJwv8TYdU8nn8+jt7Y1YYfJ8JHgudU7QSk3Q1lU+n498BMyVo5GYIknA7Uz13AEdjJNKJi9YvoPhjhalQrkjWP1ddRdJ2RV/P7A8Guvp6Qkn3MkHr0ReHUvQ5wHRwGFShZPj9EhV5AWW3QiTk5NlI3JBB8/FshgaGsKuXbvC43TWmG6c4jqSZAW9323E+st6roLUHYYvKO1jeHg4dFHGLcDodtS+zkq/b9/7PXfuHHK5XGR0LWnNep6LuDY0EkPSI2ldVhdxZ8r7k9GqT+Ho5+zGVnS90McKIqf7XkWubdu2eTP/XPedL24l99V1qlAohJ03UP6lQABlHbgrty9DU87TFqQ7sBH0Fyjd/WNjY6GC0ItXCtpC7+rqimR1yfFawWkZ3diL9gLoZJp6Y4qkAjL61J2Ta8prV4Q7ipfKLKNId7Qp93AVi/jOfbnu4+PjYfaTuJLcCqV9xLpBijtCAnG33nprrCtFKuXhw4cj232mtvZv699y7JEjR7yuOHHR6fRod3Qpykgas5t+mYT49F1Zfdl3buxEN1pXMemRpy9+o0fLPrTcbpaaew0tk6BXGNZJD/ra+mNm+nq+wYMbN5D5Rvq+vrroWt3A8jsElpWD62Jx106T/3XWXZr1yMSqF6vTfYfyDPS2OFeRWzaNLOiZNADRdV/3BxMTExgYGCiL8cg5MtCQVOs4OfScotHR0chilvrajV5nCzBFsiLEfeGOlNzRjDQknVUjfu2kmcVCLpfD2bNnsXv37kil1JVLXECyT0b4upOQzkU63kqmvy6TVEbdEenjtI9aFJP+prU7octXuV23h4zG9+/fH1EUcY1Xpya7aY/i/opTIvK3yC9IWeIsH92p6oFFLpeLzMtxr+feS64j3/HWH/qS0XqcEtKIknXnRPg6al/d05MkJZXZHRS5c4tctLJdSVxABjXaVSn3SOtWdef2aCUh717XZb1d4yYiVFJibp2WZABfcklPT09kmRKd6OCmg8dZh67y023Pp3y0NdgITJEkIBUlroG4Szv7jpGRq3aLydILLjqtz7UwfJli4i/VHahUIC2zm/cuFdP11wMIR6E6Zz+u/G7HXQlfh+GOstPgU4iVMo4ARJSpfm8+t422inzJBUBplKoVi3al6Lohy59rV5UPn7tHX1MHcnV9SHr2leYWAZXntojcruWhrToZrYus+hh5XnqZdI1kVknd9qWsunVUVnZwXU1iQcu1pD6LRaSvo8uWVHb3uzDC4cOHMTo6GiYkyLwceSaVFLCv35D5MElZYhpZgkUPcN3BXyMwRZIC3+xuQc++TSKTyYRLbsvIVc+10HGVnp6eyDchgOhIWEb+OmAncmrrSEYxOo4jI2ltqfjQIyNd0XWnqkfDEhT2dWpuZ6dH+xIU7u9fzrXX99adS9xI0hdo1aNUt4zS6OT56ziADpC65RUqZVXp/drlp2NC+h7ipjp69Cjm5+fR3t4efntFx0ZcZemmZMcpU10XBKlH0lHKTGrpFMXnn8lkylyIPmvch3zlUhSpS1J8Lk3GYU9Pj9c161ttIm3iQZJiyWQyuP3228NriWKXNuW6R9Og341Ys1qJxAX1gaiF6XuWMni1GMkaQFwmOlApoxtpXHHffHBHjnrJDW2huK4WH9rnrJeN8FWguIoXN3p1G1kls146H7m3+NfF1y+duHzM68CBA965OPo+g4OD4XN03TKZTCb0UesRb5yc2hXgTtpLOgcon+fhui98xFksWpmIFSbP33VTZDKZMFliYWEhsuyN26n6yuCubiADCKD8a3wigxuHcqkUp9BWqygi+Z10nl6ttlKyh5RFX1fcqu7grtLo3ZX98OHDoSs1CTe2J+fqOVXuIEEC9KIgBgYGyiZhuhlyvkQJd6AkMZd8Po/Xve51AErJKO5ALY2FXktMkawA6bzd4JaQ1IHncrlIhZRKoyu/XgYhaUmFONwJaVo+d+FD2eaTFSj/bK2vY9WKT//tdkC+LJyVNPpKjUKudfvtt3snxWncZTY0kgIsDVSQ961l19fXy+Toe4sy0xMHXeTbM67l6XM3ybtJyrYCyrOzBF2/tDKpZF2tBB3od7e5naReFiRNYF1wFaobk9QJGXLvJAvTlVeQ55vG0vCl6rryAuVrmaWJx2i5dWxS0DG2ubm5SPLFSp5rNZgiqYA76ve5B4DliuSapZq4Wdo6CwkodcrFYjG2M5SZtK48leIUOnaSpgJXMq19Qdy0o1L3PpVIO8JyM4Yq3XclvmRXzrhn6Fo1lXL6ZWTrrnuVFv3MXXenVmpSx/Ty5bJPd86+Wev6mepUW6Ck2Nz4ie6EBZFPvhniyih/64wxwc3qc6mUTOLiW9DRh3Y9u/K4WZiuC1HXg+PHj2Pnzp1lE3Rdq1wnmbjuOp3dKHN3hoaGMDU1hUcffRRA9JO9Ir+5ttYQupK4y2JrtAkKlBqGGxCvRH9/Pzo6OsKKkzQxy7U+ZDarL+AW12gqjdb0MYI7x8bNLEmKv8TtHxgYiCwj4d4jDt8x0rG4VtlqcS2puKCtHm3q56nnlLjl0qm/MoDQ7ypOEYksEifQ6CwxnxtR9mmmpqYi6dd6HgSwPMnRt8ZX0qRDff1cLof29vbYbLypqamyb4Mkpev6zq9U95LSxt1O2Nfh++aM6EQD16VYKUvMxbcaRKWMy7iliBqFKZIK+CYlSqql+EN16mWSP1sakutW8PnYt2zZUiaLr3PUjebAgQNlQVS5b5rRl8it5ZE4hztC1GUUOXyVPU0GidxPvgbpLj8h116J9VDpvr73kiZrbCWuAt2ZycfD9Nwh8Z3LcieyHIy7JHullRTcTlNno9XLveF2sJXqlh7oxLltfe5QTdKAyiUp1Tfu+KQBlU8xDQ0NeYP9rgs56T7uPJckpC9yU72B6CKmenHKRmVvmSJJQVylFvL5vHdxRd9LFAWkO3eppG6+u+/e2v/rk1NnaiRlOVXCvb57ju/+cW4HrVTTmNk+cz+NNeHGMuLwPee4Y1a6D/B/XyOOc+fOYXKy9G0V6WCHh4cjytvtNLRSX0lGkg+pxzobCYh2cmnrTBw6XRoofVtcKxZ9rzTZWq58bnmkDsj/YknFxQR9g6Q0ZDKZiPWwEnduGpLcjW5Glusu1wk5jcAUSQrc0atO5fTl57uzoeVlun5WubbPn+ySpsOo5ApyA+buqG21vlT3vJVmgcVdT19HdwZaCQ8OLi/N7hsZplVeev5NLTqDSu9LEhKOHz8e2T44OOidiyHXrITUAe0ei1P67vXT4ju+UgKApq+vL5K+LvjmLMUlfci9dOed1hUqVFrQMgk9WIgbdKX1AkiadBK+eBZQ/rzdibmNwhRJBbQZKhXdnUkLJM/qXYl5DSBceVYmqGnrp9JI2dcoKmWq6MpYr8BcnLJJO/LXI8q4NY0km0oyeOTdpF1vyB2Vu3JU6iSlQ9A+ftkmAxDdCUxOTmJ6ehrA8sSyJOLmkbjoCYxJcTFf/MRXpmqJm9DoG1glnS+jcm2NuJZJpWy9laAVgluH4p6br67oTDYtc9r7r0TORqf9CqZIKuCala4FokcIblDQdXGk7ZDc+8tIeaW491tpY0oz63c158UdPz09XVHW1TyPuBGnvsdK3H1pZEwjz+DgIJ555plwRWRZeDPJNanP98kmpKljujPUo/tKyl133O6xlSb0AaX3LHNlNEkWRVLa9mpIG5/wWUD6/Djcfb6YZxLu+Tq1uZKcbuyyEZgiqYK4CURA9AVWO6pzvzZYa9LI14gUwpXiuq4qPaO4gGmlY1aKK8/QUPTDZ1pOcfNoKsWj0taFlcY24jrHWrk85Rl0dHSEZfYdU03bWU38D4iuK1bpnmktitUMHONolqWRFlMkKdD+SXfdo7jjknDTRmtR0YRGjUCESqPzWiqglQTeV9KJNkJJ6u+DJMkl6bWV4mVx1/EFzFfDaq3YWp3vw+fqrOS2rYTr4kuytmpBtddMUk6+jNFGDQBNkawS12Wls2x08K/RrGREWc8G0wya8bzT4D5n3QkMDAzg2muvBYCy787EvZe0sYUk4hIv0rIa96V00vl8PraTq7bzcwdpcfv0/WpRb+rdhuKsVF+Gp5sN2ghMkaQkblKXJk2KrI80jbmRo4u1TlzjX0vPpxpZ9FygSqy0XlSqa7W0jvX16uE+bbU2Uc3ArZLF7yqURg+qTJFUIK0rpVb3Wsn2ZlJth1OpTG6jSxpprub6jSLNqgH1YLWJEmudOEVYz/Ks1mqLy8KsB81+n6ZIUlCtH9ZIT1yjXcloLkkJ1YJa+bmB8m9zNCuoWsvAcKX7zM7O1vUe9ULcc61CI13XpkhWQaNGe/UKdNa7szBWjp5dnfY7Eq3UqWm2bdsW+WLgWqARz9K3MvJ6wRRJCtKO1qqdt9GKbIQyrpa0acl6Jdo0s7NXMn9htWyU91qLdOe10t4bZVX6MEXSJJpd6daKDJVYD9ZTmqCwTLDcv39/w2Jy9aSendpa6bjXat1shkIxRWI0lLXSCaxlmpnGXI94kvvN+lal2rq7VtPTa4EpkhXQ6g2hltT7WTTTTK8lVmc2Lhvp3ZsiMVqC1cyVMFqftWjBtsrgxiYkGlWxXrN91oNpohYAAAzSSURBVDP6O/et/H7WYsffbLLZLGZnZ8MVDPR2YH08K1MkdaAWFePZZ59FsVhs2Uq2nhpJvbFn1NrY+zNF0nSsw01mozyX9VbOtbzwYa1YK3KsBUyRNJBGK421qKTWkixGY1gv73y9JIDUA1MkDWAjV7y1qMyM1sTq0trFFMkGwBqg0aqsxTq7FmVqNqZIaoh12PHYs6nMenhG5v4pR1b/9W1fL5giWaPUY2G7RnZU1qG0LutBoRmNxRRJk7HGaqxHZWv1emPxknpenIhuI6J/IaKTRPR7nv1XENEhInqKiB4mou1q338joqPBv19W228loieJaIyIvkVEV9WzDLWk0d8zaLXvJxiG0ZrUzSIhoosAfALAawCcAfAYER1g5mPqsA8D+Bwzf5aIXg3gvwL4NSL6OQDXAdgDoB3AN4noa8ycB/BJAHcw8zgR/RaA3wfwtnqVIw1pl5ffiJiLy6gVG7kdrXXqaZHcAOAkMz/DzEsAvgjgDueYXQAOBX9/Q+3fBeCbzHyBmV8A8F0AtwX7GEB38PelAFp/zW0gXCG1HphlUlvq+a4MoxWpZ4xkG4DT6vcZAK90jvkugDsBfAzAvwVwCRFlgu3vJ6I/BdAJ4BYAYsn8BwAPEtGLAPIA9vtuTkT3ALgHALZv347p6elalMmLfDp027ZtALCqe8k1JLtjfn6+RtLF36ORyD2Tnk09ylwPavkc5+fnUSgUItua8X409a4nrfKea8l6L3M9FQl5trHz+3cAfJyI3gbgEQDPArjAzF8nop8C8E8AfghgFMCF4Jz/BOANzPwdIvpdAH+KknKJ3oj5PgD3AcC+ffu4np/2LBaLAFBVlpXvGrWWea193tRHK8hYi/etGR4ersl1akUj3kErvOdas57LXE/X1hkAO9Tv7XDcUMw8ycxvYua9AN4XbJsP/v8jZt7DzK9BSSmdIKKXAriWmb8TXOJvAPx0HctgGIZhVKCeiuQxAFcT0U8QURuAXwFwQB9ARH1EJDK8F8D9wfaLAhcXiOgVAF4B4OsAZgFcSkTXBOe8BsB4HcvQcMz/bqwVrC4aaamba4uZLxDROwEcBHARgPuZ+Wki+gCAx5n5AICbAfxXImKUXFvvCE7fDOBRIgJKcZBfZeYLAEBEvwHgy0T0Y5QUy6/XqwxpsUC2YRgbmbpOSGTmBwE86Gz7A/X3lwB8yXPeAkqZW75rfgXAV2orafOxNNnWwQYOhhGlrhMSDcMwjPWPKRLDMAyjKkyRGIZhGFVhizauMcT/Xs8JlIaRBosFGWkxi8QwDMOoClMkhmEYRlWYIjEMwzCqwhSJYRiGURWmSAzDMIyqMEViGIZhVIUpEsMwDKMqTJEYhmEYVWGKxDAMw6gKYnY/Wrj+IKIfAjjVbDlWSB+AjTa93cq8MbAytw5XMPNLKx20IRRJK0JEjzPzvmbL0UiszBsDK/P6w1xbhmEYRlWYIjEMwzCqwhTJ2uW+ZgvQBKzMGwMr8zrDYiSGYRhGVZhFYhiGYVSFKZI1AhFtJaKHiOhE8H9vwrHdRPQsEX28kTLWmjRlJqI9RDRKRE8T0VNE9MvNkLVaiOg2IvoXIjpJRL/n2d9ORH8T7P8OEV3ZeClrS4oyv5uIjgXv9RARXdEMOWtJpTKr436RiJiI1kUmlymStcPvATjEzFcDOBT8juNeAN9siFT1JU2ZiwDeysy7AdwG4M+IqKeBMlYNEV0E4BMAXg9gF4A3E9Eu57C3A5hl5qsAfBTAf2uslLUlZZmPANjHzK8A8CUAf9JYKWtLyjKDiC4B8NsAvtNYCeuHKZK1wx0APhv8/VkAv+A7iIiuBzAA4OsNkqueVCwzM3+fmU8Ef08CmAJQcYLUGuMGACeZ+RlmXgLwRZTKrtHP4ksAbiUiaqCMtaZimZn5G8xcDH4eBrC9wTLWmjTvGSgNBP8EwEIjhasnpkjWDgPM/BwABP/3uwcQ0UsAfATA7zZYtnpRscwaIroBQBuAiQbIVku2ATitfp8JtnmPYeYLAOYBZBoiXX1IU2bN2wF8ra4S1Z+KZSaivQB2MPNXGylYvdnUbAE2EkT0jwAu8+x6X8pL/BaAB5n5dKsMVmtQZrnO5QA+D+AuZv5xLWRrIL6X5aZLpjmmlUhdHiL6VQD7APxMXSWqP4llDgaCHwXwtkYJ1ChMkTQQZv7ZuH1EdI6ILmfm54JOc8pz2I0AbiKi3wLQBaCNiArMnBRPaSo1KDOIqBvAPwD4fWY+XCdR68kZADvU7+0AJmOOOUNEmwBcCmCmMeLVhTRlBhH9LEqDip9h5sUGyVYvKpX5EgA/CeDhYCB4GYADRHQ7Mz/eMCnrgLm21g4HANwV/H0XgP/pHsDMb2HmIWa+EsDvAPjcWlYiKahYZiJqA/AVlMr6dw2UrZY8BuBqIvqJoDy/glLZNfpZ/CKA/5dbe5JXxTIHbp7/G8DtzOwdRLQYiWVm5nlm7mPmK4M2fBilsre0EgFMkawlPgjgNUR0AsBrgt8gon1E9OmmSlY/0pT53wH4NwDeRkRjwb89zRF3dQQxj3cCOAhgHMDfMvPTRPQBIro9OOwvAWSI6CSAdyM5a2/Nk7LMH0LJsv674L26yrWlSFnmdYnNbDcMwzCqwiwSwzAMoypMkRiGYRhVYYrEMAzDqApTJIZhGEZVmCIxDMMwqsIUibGmIaKMSvs9G6x6LL//qcGyJK5W61uVmYgeDlaDFZn7nXMiq8AS0WYi+iwRfY+IxonovcH2DiL6ZyL6brAS8n9R1/jLYPtTRPQlIuqKkf8XiOgPYvYVVv9kvNf7R0pYwdpYX1j6r9EyENH/BaDAzB9u0v1vAfAdZi4S0W8CuJmZf1nt/xhKC0rOMPM7g20PA/gd36SzYBXYf0Bp/bB3MvPjRPTvUZqk9itE1AngGICbAZwCcDEzF4hoM4BvAXgXMx8mom5mzgfX/FMAU8z8Qc/9/im49rRnX4GZvQpoNRDRXQC2M/Mf1eqaxtrFLBKjZZFRNBHdTETfJKK/JaLvE9EHiegtwQj+e0Q0HBz3UiL6MhE9Fvx71Urul7RaLa1uVWbfKrAM4OJgmZQtAJYA5LmEWA2bg38cyCVKhIJzykaHRHQNgEVRIsHs69HgOdyrjusKrK0ng2d3R7D9XiJ6lzruj4jot4nociJ6JLC2jhLRTcEhBwC8eQXPwmhhTJEY64VrAbwLwMsB/BqAa5j5BgCfBvB/BMd8DMBHmfmnANwZ7Fst4Wq1VHlV5geCjvb/DDr7pFVgvwTgBQDPAcgC+DAzzwTnXEREYyitSfYQM4ffsyCiBwCcBbATwF94ZHgVgCfV748B+GTwLM6q7QsA/i0zXwfgFgAfCWT+SwRLuATl/RUAfwXg3wM4yMx7UHoHYwDAzLMA2omolVcwNlJiisRYLzzGzM8FC/9NYNky+B6AK4O/fxbAx4PO+ACA7sC9tCJoebXaDwWbwlWZPYe/hZlfDuCm4N+v0fIqsO/xHH8DgB8BGATwEwDeQ0QvAwBm/lHQYW8HcAMR/aScxMx3B+eMA/B9RfJyAD9Uv18F4AvB35/XxQPwx0T0FIB/RGkZ9AFm/gGAXKAAXwvgCDPnUFpf6u7A7fhyZn5eXWsqkMlY55giMdYLeuXYH6vfP8byKtcvAXAjM+8J/m1zOj4QkVgPD/puQsur1d6uVqu9EcA7iegHAD4M4K1E9EEAYOZng/+fB/DXKCkKvQrsDwDsR2kV2H0ojfD/FzOfDxYy/DZKSiuEmecAPIzSFyP19h8B+BuUrC2XFwF0ONt8AdK3oBTnuT5QWufUeZ9GaQn0uwHcH9zzEZTWQnsWwOeJ6K3qWh3BfY11jikSYyPxdZQW1QNQ+h68ewAz3x0omTe4+yhmtdq4VZmJaBMR9QXnbgbw8wCOVlgFNgvg1VTiYpSUzPEgvtMTXGsLStbV8eC4q4LtBOCNAI57yj4O4Cr1+9souaeAkvIQLkUpWH8+SC7QmWlfQUl5/RRKCxOCSplrU8z831Fyf12nZLkMwA88shjrDPseibGR+G0AnwjcNpsAPALgf1/B+Xq1WgDIMnPSqq7tAA4GSuQilFxF/73CPT4B4AEAR1FyMz3AzE8R0SsAfJZK3wV/CUory341cJN9lkrfbCEA3wXwm57rPoIg3hEsT/8uAH8dBNC/rI77KwB/T0SPoxTvCJUSMy8R0TcAzAXWD1DKKPtdIjoPoABALJLrARwOVsQ11jmW/msYG4QgPfnvmfkfV3n+S1AK2P8SM59Ica8DzHxoNfcyWgtzbRnGxuGPAXSu5kQi2gXgJIBDlZRIwFFTIhsHs0gMwzCMqjCLxDAMw6gKUySGYRhGVZgiMQzDMKrCFIlhGIZRFaZIDMMwjKowRWIYhmFUxf8PNJTmyKYPp+0AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "lc.fold(period=3.0712).plot(alpha=0.2);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can see the hint of an exoplanet transit close to the center, but the motion of the spacecraft has made it difficult to make out above the noise. We can use the `correct` function on the lightcurve to remove this motion. An in-depth look into how the algorithm works can be found [here](http://lightkurve.keplerscience.org/tutorials/motion-correction/replicate-vanderburg-2014-k2sff.html). You can pass the following keywords to the `correct` function:\n", "\n", "* `polyorder` : *int*\n", " Degree of the polynomial which will be used to fit one centroid as a function of the other.\n", "* `niters` : *int*\n", " Number of iterations\n", "* `bins` : *int*\n", " Number of bins to be used to create the piece-wise interpolation of arclength vs flux correction.\n", "* `windows` : *int*\n", " Number of windows to subdivide the data. The SFF algorithm is run independently in each window.\n", "\n", "This [tutorial](gully-tutorial) will teach you more about how to tune these parameters. For this problem, we will use the defaults, but increase the number of windows to 20." ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "100%|██████████| 20/20 [00:00<00:00, 25.46it/s]\n" ] } ], "source": [ "corr_lc = lc.correct(windows=20)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now when we compare the two light curves we can see the clear signal from the exoplanet." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzsvXmUXFd56PvbNY9d3dWjulujBbJlS3geCdjxw4TBBAjgQC52IAmwAoFLSC7cdd9iznSTEOKQFS+/e8FAiOEFwl2OQ0IIwTaenm2BkGzJkttqtbp6rKpT46k659Sp2u+PU1WqbvVQUneru6X9W6tX1xn3t8/Z5/v293377COklCgUCoVCca641lsAhUKhUGxulCFRKBQKxYpQhkShUCgUK0IZEoVCoVCsCGVIFAqFQrEilCFRKBQKxYpQhkShUCgUK0IZEoVCoVCsCGVIFAqFQrEiPOstwPmgp6dH7tixY73FOCts28bjuShuTxNV54sDVefNw4EDB1JSyt7l9tt8NTsHduzYwXPPPbfeYpwVqVSKnp6e9RbjvKLqfHGg6rx5EEKMtbOfCm0pFAqFYkUoQ6JQKBSKFbFmhkQI8VUhxKwQ4vlFtgshxL1CiBEhxCEhxNUt2+4RQrxU/7unZb1PCHG/EOK4EOJFIcSvrZX8CoVCoWiPtcyRPAB8BfjGItvfALyi/ncD8HfADUKIOPAZ4FpAAgeEEA9JKTPA/wBmpZSvFEK4gPgayq9QKIBKpUIikcAwjFU5X7VaJZlMrsq5Ngsbvc6BQIDh4WG8Xu85Hb9mhkRK+ZgQYscSu/wq8A3pfBDlaSFEpxBiC3Ar8CMppQYghPgR8CvAg8D7gUvr568BqbWSX6FQOCQSCaLRKDt27EAIseLzVSqVc1ZYm5WNXGcpJel0mkQiwc6dO8/pHOs5amsIGG9ZTtTXLbheCNFZX/6CEOJW4GXgI1LKmYVOLoT4APABgOHhYVKpzWVzcrnceotw3lF13pjous7Q0BC2ba/K+arV6qqcZzOx0evc0dHB9PT0OevJ9TQkC3Vt5BLrPcAw8ISU8veFEL8P/AXw3oVOLqW8H7gf4Nprr5WbcejdZpR5pag6bzySySQ+n29Vz7lRe+dryUavs9vtPue2uJ6jthLA1pblYWByifVpoAR8v77+H4GrUSgUG44TySInksX1FkNxnlhPQ/IQcHd99NaNQE5KOQX8ELhDCNElhOgC7gB+WM+l/DNODgXgduDIOsitOA8oRaQA+PjHP86Xv/zl5vLrX/96fvu3f7u5/IlPfIIvfelLTE5O8o53vAOAgwcP8oMf/KC5z2c/+1n+4i/+YtmyduzYwS/90i/NWXfllVdyxRVXLHnc/PIeeugh/vRP/3TZ8i4k1nL474PAU8AeIURCCPFbQogPCSE+VN/lB8AJYAT4f4DfBagn2b8APFv/+3wj8Q58EvisEOIQTkjrE2slv0KhWH9uvvlmnnzySQBqtRqpVIoXXnihuf3JJ5/klltuYXBwkO9+97vAmYr9bCgUCoyPOynao0ePtnXM/PLe8pa38KlPfeqcyt+srJkhkVK+W0q5RUrplVIOSyn/t5TyPinlffXtUkr5YSnlJVLKfVLK51qO/aqUcnf972st68eklK+RUu6XUt4upTy1VvIrFIr155ZbbmkakhdeeIErrriCaDRKJpPBNE2OHj3KVVddxcmTJ7niiiuwLItPf/rTfOc73+HKK6/kO9/5DgBHjhzh1ltvZdeuXdx7772Llveud72recyDDz7Iu9/97uY2wzB43/vex759+7jqqqv4yU9+smB5DzzwAB/5yEcAGBsb4/bbb+fqq6/m9ttv59QpR2X95m/+Jh/96Ee5+eab2bVrV9MIblbUm+0KhWLDMjg4iMfj4dSpUzz55JPcdNNN3HDDDTz11FM899xz7N+/f85AAJ/Px+c//3nuuusuDh48yF133QXAiy++yA9/+EOeeeYZPve5z1GpVBYs7x3veAf/9E//BMA///M/c+eddza3/e3f/i0Ahw8f5sEHH+See+6hVqstWF6Dj3zkI9x999387Gc/4zd+4zf46Ec/2tw2NTXF448/zsMPP7zpPRhlSBQKxYam4ZU0DMlNN93UXL755pvbOseb3vQm/H4/PT099PX1MTOz4FsDxONxurq6+Pa3v81ll11GKBRqbnv88cd573udQaKXXnop27dv5/jx40uW+9RTT/Ge97wHgPe+9708/vjjzW1vfetbcblc7N27d1F5NgvKkCgUig1NI09y+PBhrrjiCm688UaeeuqpZn6kHfx+f/O32+1e8p2Yu+66iw9/+MNzwlrgvLi3Ulpf6GyVaTXOvZ4oQ6LYkOzqjbCrN7LeYig2ALfccgsPP/ww8Xgct9tNPB4nm83y1FNPcdNNN52xfzQapVAonHN5b3vb2/hv/+2/8frXv37O+te85jV861vfAuD48eOcOnWKPXv2LFnezTffzLe//W0AvvWtb/HqV7/6nOXayChDolAoNjT79u0jlUpx4403zlkXi8UWfIHutttu48iRI3OS7WdDNBrlk5/85BkvYf7u7/4u1WqVffv2cdddd/HAAw/g9/uXLO/ee+/la1/7GldffTXf/OY3+eu//uuzlmczIDa7S9UO1157rVQfttr4tNbZHB0FwH+Oc/9sFjbDfT569CiXXXbZWR3TeAdoIa9yI887tVZshjovdJ+FEAeklNcud+xF8YVEhUJxflFhyYsLFdpSKBQKxYpQhkShUCgUK0IZEoVCoVCsCGVIFAqFQrEiVLJdsSG50EdrXehcLKPuFA7KI1EoFBsaIURzahIA27bp7e3lzW9+85LHPfLII80JHwHuu+8+vvGNb6yKTMePH+eNb3wju3fv5rLLLuNd73rXeZ3mZH7d2mXHjh1r8rVY5ZEoFIoNTTgc5vnnn6dcLhMMBvnRj37E0NDQssc98sgjRCKR5nxcH/rQh5Y5oj0Mw+BNb3oTX/rSl5qTOv7kJz8hmUzS39+/7PG2bePxnFa91WoVt9t9VjLMr9t6ozwShUKx4XnDG97Av/zLvwBnTu+uaRpvfetb2b9/PzfeeCOHDh3i5MmT3HffffzVX/0VV155JT/96U/nfODq4MGD3Hjjjezfv5+3ve1tZDIZAG699VY++clPcv311/PKV76Sn/70p2fI8g//8A/cdNNNc2YGvu2227jiiisWnGoe4Bvf+AbvfOc7ufPOO7njjjt45JFHuO2223jPe97Dvn37APj7v/97rr/+eq688ko++MEPNr/z/m//9m9cffXVvOpVr+L2229fsG7JZJJf+7Vf47rrruO6667jiSeeACCdTnPHHXdw1VVX8cEPfnDN5vRShkShUGx4fv3Xf51vf/vbGIbBoUOHuOGGG5rbPvOZz3DVVVdx6NAh/viP/5i7776bHTt28KEPfYiPf/zjHDx48IwvH95999382Z/9GYcOHWLfvn187nOfa26zbZtnnnmGL3/5y3PWN3j++ee55pprFpRzoanmDcMAnJmAv/71r/Of//mfADzzzDP80R/9EUeOHOHo0aN85zvf4YknnuDgwYO43W6+9a1vkUwm+Z3f+R2+973v8Ytf/IJ//Md/XLBuH/vYx/j4xz/Os88+y/e+973mVyQ/97nP8epXv5qf//znvOUtb2l+D2W1UaEthUKx4dm/fz8nT57kwQcf5I1vfOOcbY8//jjf+973APjlX/5l0uk0uVxu0XPlcjmy2Syvfe1rAbjnnnt45zvf2dz+9re/HYBrrrmGkydPnpWcjz/+OL/3e78HnDnV/Ote9zri8Xhz3+uvv56d9cEIP/7xjzlw4ADXXXcdAOVymb6+Pp5++mle85rXNPdrPb6V//iP/+DIkdNfHs/n8xQKBR577LHm91Xe9KY30dXVdVb1aRdlSBQKxabgLW95C3/wB3/AI488Qjqdbq5fKFzTOl372dKY3n2x6eYvv/xyHn300QWPXSp0FA6HF12WUnLPPffwJ3/yJ3P2eeihh9qqS61W46mnniIYDJ6xbSXXol1UaEuhUGwK3v/+9/PpT3+6mVNo0Dq9+yOPPEJPTw8dHR2LTu8ei8Xo6upq5j+++c1vNr2TdnjPe97Dk08+2czZgJPHOHz48KJTzS/H7bffzne/+11mZ2cBJ+8zNjbGTTfdxKOPPspofTi1pmnAmVPl33HHHXzlK19pLh88ePCMa/Ov//qvzVzQaqMMiUKh2BQMDw/zsY997Iz1n/3sZ5uf3f3Upz7F17/+dQDuvPNOvv/97zcT0q18/etf5w//8A/Zv38/Bw8e5NOf/nTbcgSDQR5++GH+5m/+hle84hXs3buXBx54gL6+vkWnml+OvXv38sUvfpE77riD/fv387rXvY6pqSl6e3u5//77efvb386rXvWq5qd859ft3nvvbV6DvXv3ct999wFO/uixxx7j6quv5t///d/Ztm1b2/U8G9Q08huUzTC9+Gqj6rwxOZdp5Jd6IXEzTKm+2myGOqtp5BUKxYZCvdF+caFCWwqFQqFYEcqQKBSKZbkYQuAXMyu9v8qQKBSKJQkEAqTTaWVMLlCklKTTaQKBwDmfQ+VIFArFkgwPD5NIJEgmk6tyvnOZW2qzs9HrHAgEGB4ePufjlSFRKBRL4vV6m29WrwabYaTaanOh11mFthQKhUKxIpQhUSgUCsWKUIZEoVAoFCtCGRKFQqFQrAhlSBQKhUKxIpQhUSgUCsWKUIZEoVAoFCtCGRKFQqFQrAhlSBQKhUKxIpQhUSgUCsWKUIZEoVAoFCtCGRKFQqFQrIg1MyRCiK8KIWaFEM8vsl0IIe4VQowIIQ4JIa5u2XaPEOKl+t89Cxz70GLnVSgUCsX5ZS09kgeAX1li+xuAV9T/PgD8HYAQIg58BrgBuB74jBCiq3GQEOLtQHFtRFYoFArF2bJmhkRK+RigLbHLrwLfkA5PA51CiC3A64EfSSk1KWUG+BF1gySEiAC/D3xxreRWKBQKxdmxnt8jGQLGW5YT9XWLrQf4AvCXQGm5kwshPoDj6TA8PEwqlVoFkc8fuVxuvUU476g6XxyoOl94rKchEQusk4utF0JcCeyWUn5cCLFjuZNLKe8H7ge49tpr5Wb8qMxmlHmlqDpfHKg6X1is56itBLC1ZXkYmFxi/U3ANUKIk8DjwCuFEI+cF0kVCoVCsSjraUgeAu6uj966EchJKaeAHwJ3CCG66kn2O4AfSin/Tko5KKXcAbwaOC6lvHW9hFcoFAqFw5qFtoQQDwK3Aj1CiATOSCwvgJTyPuAHwBuBEZycx/vq2zQhxBeAZ+un+ryUcqmkvUKhUCjWkTUzJFLKdy+zXQIfXmTbV4GvLnHsSeCKlcinUCgUitVBvdmuUCgUihWhDIlCoVAoVoQyJAqFQqFYEcqQKBQKhWJFKEOiUCgUihWhDMk6Yo6OYo6OrrcYCoVCsSKUIVEoFArFilCGRKFQKBQrQhkShUKhUKyIZQ2JEKJvgXV71kYchUKhUGw22vFIfiqEeFdjQQjxCeD7ayeSQqFQKDYT7cy1dStwvxDinUA/cBTnE7iKFeLfuXO9RVAoFIoVs6xHUp/a/d9wvgeyA+fzuOqb6QqF4qJCDddfnGU9EiHEj4ApnNl2h4GvCiEek1L+wVoLp7h4OJEsksmUuYA/IqdQXLC0kyP5Wynl3VLKrJTyeeBm4ML+ALFCoVAo2qad0Nb/mbdsSym/sHYiXTycSBY5kVRRQoViMzCulRjXSustxoakndBWAZD1RR/OVw6LUsrYWgqmUCgUis3BsoZEShltXRZCvBU1akuhUFwENJLraoTl0pz1p3allP9HCPGptRBGceEwcuAFAHZfc/k6S6JQrA6ubdvXW4QNSzuhrbe3LLqAazkd6lIoVoVdvRFSwlhvMRSKOTRyIruVQ7Ik7Xgkd7b8toGTwK+uiTQKhUKh2HS0kyN53/kQ5GJkV29kvUXYMJijo9iZDOpFEsVGRT2vi7OoIRFC/A1LhLCklB9dE4kUCoVCsalYyiN57rxJoVAoFBsQlWBvj6UMybeklPZ5k0RxQaFGaykUFw9Lvdn+TONHPcylUCgUFxW7eiMqN9IGSxkS0fL7lrUWRKG42DBHR7FPnVpvMRSKFbNUaEu9K6I4b/h37sQTjS6/o0Kh2HAsZUguFUIcwvFMLqn/pr4spZT711w6hUKhUGx4ljIkl503KRQKhWIDoObWOjcWNSRSyrHzKYhCoVCsN2pKlHOjnQ9bKRQKhUKxKGc9+69CoVgd1ACDzYkKf52J8kgUCoVCsSKWmmvrMEvPtaVGbSkUK+BEskgmU76o56lUvfsLg6VCW2+u//9w/f836/9/A1AfLlYoFBcc7UztoxLyZ7LsqC0hxC1SytY32z8lhHgC+PxaC6dQKBSKjU87OZKwEOLVjQUhxM1AeO1EUigUio2La9t2NSvwPNoZtfVbwFeFEDGcnEkOeP+aSqVQKBSKTcOyHomU8oCU8lXAfuBKKeWVUsqfLXecEOKrQohZIcTzi2wXQoh7hRAjQohDQoirW7bdI4R4qf53T31dSAjxL0KIF4UQLwgh/rT9aioUCoVirVjWkAgh+oUQ/xv4jpQyJ4TYK4T4rTbO/QDwK0tsfwPwivrfB4C/q5cXBz4D3ABcD3xGCNFVP+YvpJSXAlcBtwgh3tCGHArFhmRXb4Tt8eB6i7Gu+Hfu3BAjtkYOvMDIgRfWW4xNSzs5kgeAHwKD9eXjwH9d7iAp5WOAtsQuvwp8Qzo8DXQKIbYArwd+JKXUpJQZ4EfAr0gpS1LKn9TPbQE/A4bbkF+hUGxQTiSLnEgW11uMs2Kpb5RsxvqsBu3kSHqklP+vEOK/A0gpbSFEdRXKHgLGW5YT9XWLrW8ihOgE7gT+erGTCyE+gOPpMDw8TCqVWgWRzx+5XG69RTjvqDpfHLTWOZMpA5ASxnqJA0A+n3fkWKGeWKw+F/p9bseQ6EKIbuovJwohbsRJuK8UscA6ucR66uV7gAeBe6WUJxY7uZTyfuB+gGuvvVb2bMK3vjajzCtF1fnioFHnvCzWl9f3K4TZjpm6HCu7F0vV50K+z+0Ykk8AD+F8k+QJoBd45yqUnQC2tiwPA5P19bfOW/9Iy/L9wEtSyi+vggwKhUKhWCHLGhIp5QEhxGuBPTjewjEpZWUVyn4I+IgQ4ts4ifWclHJKCPFD4I9bEux3AP8dQAjxRSAG/PYqlK9QbGjU9CHnj63x0HqLsKlZ1pAIIV4G/lxKeV/LuoellG9e4jCEEA/ieBY9QogEzkgsL0D9XD8A3giM4Ey58r76Nk0I8QXg2fqpPl9fNwz8D+BF4GdCCICvSCn/V/vVVSgU55NG4nmx5PRGQRnrldFOaKsC3CaEuAH4YH3E1NAyxyClfPcy2yWn5/Gav+2rwFfnrUuwcP5EobigaHgiis3HRjeYa0U7w39LUsq7gKPAT4UQ21liVmCFQnFhshZDW5caSqvYPLTjkQgAKeX/FEIcwHmnJL6mUikUy7BZQiYKxcVAO4bk040fUsofCyFeD9yzdiIpFApQcXvF5mGpD1tdKqV8EZhonQerzsNrK5ZCoVCsP8rzbY+lPJJPAL8D/OUC2yTwy2sikUKxAOqB3pyo+3VxsNSHrX6n/v+28yeOQqFQIa2FUe/VbFyWCm29fakDpZT/tPriKBSKjYryLhbnYveYlwpt3bnENgkoQ6JYNy7WB3azobwIx8hkMmUu4Km2lgxtve98CnIhoR4exUZDtUnFWtLO8F+EEG8CLgcCjXVSys+vlVAKheI0q2EExrUSALuVHZnDciEp5fm2Rztzbd0HhIDbgP8FvAN4Zo3lUijmoB7ozcHFniu4WGlnipSbpZR3Axkp5eeAm5g7/btinbgYvsZ2Puu4VFkXw7Xe6KzFZ3lrp8aonRpb1XNejLQT2irX/5eEEINAGlAOsmJdUTH/88fF6GWcbZ3X4tpspuvejiF5uP5p2z/H+U66xAlxKRbgRLJITSup7xtc5Ky2EhjXSrgixVVXKmutrM7W0G8m5dluZ2aomCRUzHAhB3La+bDVF+o/vyeEeBgISCkv7A8QrxDXtu1MACRX/8G/WDifYaQLIWS1nAJ2bdt+PsVRnAXL3bvNYFzbSba7gTcBOxr7CyGQUn5pbUXb/DRGyqxlAzgfZZxPLgSlvtr4d+7EFVnZdbkQ2sdaKNQLJXLQuDZDxSRw/kO+7YS2/hkwgMNAbW3FUSjaY7HhrGuhbE4ki4yvQrhyrWTbSKyHwVrJdW1X4ba2gdZylgo5NpX7Wcq00e5pO7RjSIallPvXXBLFWbNWD+1GcqXblWE1Hr6lypqvQM6VlXqQK703G+nezmcxg63CnEuzWh2dldCOIflXIcQdUsp/X3NpLhAayuJCYa2Uz2K9vHZY75j/YtdkviKav990zljwuHMt/2w5l7a5UYzO+QzjtlPGRKSXCau06in0Rj1fu6cPOLt7vV4vnrZjSJ4Gvi+EcOF8v13gfHK9Y00l2+Aspkgay5vFmCxlJJZ6p2KxY9aazdJjbFyb5a5hh1j8HK3XebnztaNkTySLTOcMBmLNCSpW3aDNl7Od87eeo9GpOJt2udHzhP6dO/FEo+stxprSjiH5S5yXEA9LKS/6b7W3q8hWy81cidJu99j1fhAb5bcq1eV6+o3j5rv051KXjaqAzobFwkEbqW6PHpsFICJM8vL0/Wzcs0YdzjVUc651bm1bjXLP9lznajjPpozlGNdKuOd1FM4X7RiSl4DnlRFZmoUa0kqUfyZTnvOwLecBLbffYuUsxPwHu51jzkWGhVhPxbcRle+zoxrjWqkZ5pjv6S6nbNeyTuczNr9Ymzwb5l+LhdpzI/TYatTmb1vJs3U2tMrbmFlhsbIb8o6zFdc65EraMSRTwCNCiH8FzMZKNfx3ada7l98u41qp+YDAxpN3sTDGYgplLXqxu3ojmKOjmMUk/p075yjQs1Eau3oj5xzybJTTeq+g/fo2ZE4WzLM67mwxR0cZ10pMRHrnlNEaVpvMmRTlmQZo/r1d6lq17tsY8kq8/bzZfCO4XFlb4yGu2xk/Y9v8UFy7eb+l2vF6dg7OlXYMyWj9z1f/U7RJqzFZrDfU2hiW6uGttmGab0BaywBHYTUe/PM9GqSdnm5DKU/njDkP7ciBF5jOGYwPbV21kVYLMb/n2mD+fWpdNkdHGSqWFh0osNib0q1lNcIW0zmDZMHk2VFt2XvUeq97o/7murVSWHJqEnfIgPieszquVZ7WOq9m3rFxjnGtxOFErllGo7yzae9nE/ZqjTDUTo3hzhlUh7aesd/Z1GG+HOuZP1zSkNRfRoxIKf/wPMmzoWnXlW/c2GdHtebys6MayYLJG/dvWRXlttJeyeFEDk232DNwOgnYMB7zmd9QV8uozVcQY1p50dj5ckpk/kM0PwwxP0zQ+F07Nebsu29Pc/K+E2w/Y9/avNEwDYU+X5HPbxuL5W/aeehbjeX8+7JYHPw7z5yas71RfsMT2Tcca66f34lplLkci+WrAIbqXc1kwaSqleacryFDSMw97nAiR2/Uf8a1SxZMxuedYzXRdKt5XZMFs2lkl6IdL2m57SPHkuTKFYYivcDi4erGdTZHR6nVvbxWGvql0W4nWs53vr2WJQ2JlLIqhLj6vEiyAVnqYV/oBrX2dlof/BPJYvMhWqqM1oY4mTMZyS7f22yco/XYc/EgGj3c1ocpWTDxTSUYKSbn9KIXSoquRoNdLOTRkA84I7yQLJhMHj7GQCyAa9v2Zj2oi7PcKKezYTpn4FrmfA3ZG4nl1mUtkWsq8sb1SqXONNytx7eivTgCgLVl+IxtjfP94NDUgnIDTWU9rpWonRrDLIbO8H4W85Tnr2sw32NoXIfeqJ/BBe5jsmBSKpXYLn1N2Xqj/rYTxI321riuu3ojmMXQggMvWmkoYyK9jGsltBdHiKR0kuEdJAsmoykdWNxAz28vDSU+vz26J8adH/E9TQ+TeQYgHl44sLOQNzudM5gsOmVNVxduKxuBdkJbB4UQDwH/COiNlRfrN9vPRWmOayU03Vq0AbVT3mKKr6HMW13zxr61U2McTuQYv3R3M1k7X0H5phLAXOXUMHq9UT8UnYe9Wi+jYdjmy7OSxOvWeIhfPHmIFyfS3Hz7Dc36DMQCjBw4ghb2MXDpbmDuA+2eGOfqWUdxHi7EiUdKDMUCuGam0A5No/UOwKW7z4jVN87x7KiGb8oJb1S1EkM4dZ2oexmt132xt5OTBXPO9WqUtVjoq1UOWHr47/xygroFQHmJ3rNW36c36m/K1VCUwZkJau7CHMW2kKFtJPhXMhrONTPFJMC+PQwVk9S0EvuGe5nOGfwinWM7c6/N/Ou0nGFZSJ7Fns1Hj83intAWPGfrc9m4dgvlL1qfr4VkaRi3yYkzt01YJfJ5k5KUc45vPEsLdRrAad++gkmysSIyV5Zm57TTCScObXG8kvUIcbVjSOI4U8f/css69c32s2A6Z5ArVwAWdNUX8yYWe5halXZD6TYaVeOYxlDAxWgqoqLzuxoLnKEU59PqHTXiy40ediN01zAAr93T19Ywx9ZvQRTNajPMMZ/5D/K4VkIrmJTqPcnQttM9Q023SGTKlISOVc+hLERDOc838K3eo/biCJpu4a6HAOfnkRqKaDSlo+nWggnZxRjXSuTzeS6XgTOuTaO+d12/jRPJYtPTiId9WC3XxzeVQJuCE/U65soVYkFv2zIsxlKGcP4Lc/OPcy3gfU/njLZ61M32MK8X3xouXkqpw+mEv2vbdmqnxtBaQmeuyOlOUCzoxbdM526omKRWLDFNlIFYgOmcwaPHHNUeD/ua92WhXFmteLqtHJ3WMaRzXYp1g9XYr/XdnumcMWd0VsO47RmIUl3GuDZyU+P1VyTPZ26zndl/L/pvty80zh3a6521KvkG83vvrWGlhW5+64iXR4/NnvFS2fzyfFOJpjLWdAtrnkE5nMhhnTxJLOild8B5QFz1hONiIbgGjYbfUKDaiyOM6BYTkd45CrldD6WhFDTdomRVgbpy1C20sI9AegbCp5OSi42yGU3pTCRyDHUuWdwcNN1iCOgvadQmxpmzupLKAAAgAElEQVQGjk0X0CLOtf3Kj18incgx3BWckxid73UsRMPTY2e8mQcItigQcO5DqVTi8h1zj22EQnuj/uZ17I368RV99Eb9NDq9A7EA7uJpozKulSgYcw3JselC07jEpc+RXUsipyYhvrupuNuZKWAhr7g1bzadMzhWDDAkI1wmis08Sa1YYiAWoBoJcDiRQ0wkmChkSBZ2zOk0tF7X6epp77e1jMZ1abT/E8niHG/xRLLIZD3sNNiy/th0wVm3rxftxRG8yWnivQPNa9mRmkICj+IYiCtFnskJqNVlmLCc8pIFs3k9G0q+td035JrnmJAtVzDrb1A0jt8aDzU7YK3Pc+M6N8pq/G50URr7Nuo0LU4b72TVxGp5jjdEjgRACDEM/A1wC44n8jjwMSllYo1l25A0Qx3FJGYxxNNWGDjT+jeUYbnfaeaxoHdOg5s/bHA+swULe3IEd9TfTKItlMdolGWldIq9g/RG/XPCG4H0DEXgRL2n3Kr8GvI0QleNeoWOPc+uvTtx33ATk8Ukz4xqTCTd7OwJN0MmuXKl2ROXU5NEtlQgvKN5PvfEOG5oe+ROID1DLJdF0y12JKfpByrhgeb24sgJAukZau4rOFCNkiyYXBn1kw16yZUrhLRZvEEvE8N7qfTmIDNKwbCbD1vDU2jtAebKFSLlivMhaVoe3PqzN5rSsYoGw13Bphyto9kWovHg/qTuFcXrxvRwIudMqRHpZXfL/tlyZcmE8vx82zOjGhMRN3urWbSwr+kRtpIrV+Z0CGJBrzOoYirHgadfAJy2Mb83DE6P/5m6Ir5+Z3xBz7KBOTrKkR//nErvANaWYZIFk75jB6kAowEPoBO/dHdLB8UxphG/m1jQ2/xiXnHkBFrYx9ab9zf3HRg6MzfY6MDMb/8TkV4OJDUGlum4ZE+MMRAL0Bv14yr5oKQxWn9uogFHFU7UFTeRuWFLIuEFz9kaaQAnBNoctVb3YBrtqWA4+xZ7Byly2oufn/QHJ6xVOjVOweVMINLaIWy05cbzmyzk8NY7dt7kdF3euR7dWtNOaOtrwD8A76wv/5f6utetlVAbhfludmOUz+zxUWTUD/HTKqF1OCGcjr3Gh2M8eixJvP7QH07kFg3fwOkk3rkQSU7ik+Fm4+6N+ikYNqVypemCO2GfEpfWtzd6rEOR3mYPVwdkcvaMT5BquoWmW0SSk3QYNvTspmDYVEoWVq5MZN5D3uiRbo2HFkzuNhTUQCzA84YNzO3d9Ub9jHHavS8Ydv3BdMJM2RNjpE8mKMe6iQY8hMI+pnMGRd0iGvBwwqzAyMtMJIMMFeNOaKIaZaiYZMAdoCM1BQEPfa90ZJqI9FIUhTl1DnrdxMO+Ob1l31SCIXeMrTtP9+JP/eIo0u/FHA1hJRKOJ+U6PYtQcGaCvokJ3JFLOdGSY5rMmTx6LOm8DzHvWoNj2H0FE2vLMNaWYQqJF3DnTsFAR9ODfXZUa3YEbvDqUIUJ3fFKcuVKU4HtazmvpltUqrlmkvmtdY/ktCI7xUQxyYm6RzVy4IhzHW7ef1pZF3WMqSlSmTL9OEp0Bgjm0vgNN2ZnT9371c7oSFUmJkhnyvRed/mckYOtzM/9NZgfxoXTxsgd9TP2wghGdz/VISfvVTo1TsWw8Uqd2eM+kqE4/TjtKyad8JaVdNoXkbohroehNN3ixZ8fYyaUnNMpbFAwnOv70MFJdvaEYbgXrZqDeZ593rABNwOxQPO9l2qsd47X0VqXsfQM0zmDlNeHOTBMd39nwx7NGdCQrLeN3qgf7cSY06nidGh6qZcYV5N2DEmvlPJrLcsPCCH+61oJtJFoPDCuyNykW7auMMe1Eoez9hmjToaKSQoljRniHD6WRI6fItgVZLol3ASnk3SNY4eKSbInxqj0DtA32MdlMae3OVF1vA5vcppiKE5xBkYA90AU78wUXt0iFoo3H9SCYc9r7Db5lN5UKLppU6icVtyNHmyyYBKcKTBcl7EHGNy3B5+cItYSlikYNrpZYXhmgh1TL5Epm0zSv+z1HNdOT7ndCNl4Z6aQJT/BXJpuwxnBY0xN4ctncLUonj0DUfSETvbEGMmdHWi6Rde8B1DTLSL1OpXSUK5U6dc1+r0RYG7uor+kEcnMUuiaG+dPZEp0pqe4xFsm7eogVami6Ra1U2NsjYf4QcEmXe8IVLUSE4ePAZDRLcqVKs+OavSXDHb2hPGPnMD7syr8X7eyZyBKthRk7IURjk0XKPcPNdtB49rUiqd71HvqIUcm8mi6Ra5ejntmkn6fh/6Sjf7ccU6VdkIo3lS4jTbQCK3Egl4GX/oFr3h+HH9vH6PC6VmX4n0URYxiLkXYf7qtNDzaYv0+N9qoVj/vQolo3bTRdIvLtVHcsxOcsP2ULQ89pAjOTGDheJylNORCcToCHtBBN53755qZQtMtqkNbceMofk23iOze1TQavqkERd0i3j80p7PSkOd079xsdp7gtAdZ7Opj2B90rnmIZueBfJLeqJ9i0MvUiy+TnC1Q7B2kQxQY7gqSK1cI1L3WgSf/Hd2sMrrvNeytZp32gnM9O1JTWEUP09E9zZftGpGEhl/dmZ5Ce/F0J6/hYTWM9tC+PY4sM87+Yb+HnlCAyPgRoi/lKO7dT2T3rmbdG3ojWTDxFU1CWgo7bzAaiuOLnDZkG8WQpIQQ/wV4sL78bpzk+0XBDw5NMRGxmyGEZMGkVFe8tZDBYu9o5soVEmaZvLvCUN1tfmaet9H6fkCjMXhxHqRsOkd3sOr0OoI7COpWU1U3QjkuEaPHyOGZmoYt0DmwneyJMSKZApWSBwYupf+VOymldCaMCrppc4272DQirpkpCMUp9g4yP0WcK1XgqWe4vKSh6XESmRLDXSFy5Qpus9JUPuV6XmPX9MsMd9hE97x6ydzB/BfCXjE+SQEIFLJERJWxw8cYSqUpA4cSObxDQ87DU0ziyWUACIYmiJQrGN390N3vKA1tloY6tLYMk4n0Yiay8LKTnJ0Z1Qhps5TifQzVQwPbZZFUycN0zjEmx+oRnIFYkIETR7EzJWL+CBXXIIfriso6OUn/2Itk/WX0SG+zpy2nJqECYy+MoJdz7Nq7k4LfiaNXtBIj0wUqmTLBXJqKYWOVK0jDBt/Cc5/6phK4i35cM1NUJnIUQnGGu4JIn4eyZZPIlOkx5xrSicPH2DZ5nMCWLdAXY5903hXKtuwTzDmPbik+14A2hsQWdYs9u3cxkQxSGjvFyIEgXLOX/pLmhGDqo/sePZZkqJikUrAoB2w6AX1qlmomixnsoxjtYHhLF4GSRk5EMbr7HWMydgqPX+D2+SlXqiQLJp3jkwRzGdy7tuOamWJg+iTloBedXWi6hRw/haekOcav7hXvGYieMXpQ0y0C6Rm6Ei8D9RDrE0+za2yUsXQ/Ca+H4VduR9MtHj2WJJCeIZszKPcPceVAlKmglxmgOzNNROjkAoNM58qEqjbR9Ay+fAbd30HBqJCznTyHbtqE/R7HI9Zmsf8/DWNmhnKsGzvS29QbRbPaiKCe7kDU9Xt/KUO/t8x8ZkJxJiI9XF3SwHCe+win36tyT4yzdWaKooww0T9EVzBG2KzgKWlwRKPY3c90i+FZS9oxJO8HvgL8FU6O5Mn6uosGTbcYOXDE8QCS0/izKegJkz0xhiWiHOsddMIdAMSYzhmOmxxyeir9tTwx6SUnKhR7B5uubMPLqPQOQKSX0qlxQr94Fk+4A151XbP8riM/o6eco8ssYDBFIr4Ladp4Ro5Rs3TwhfFqKVxhH/Gwr9mwGz0WTbfY9uIhDLtGYWgYvW4IRlM6MyE/0sgRfeFpAMJb+vAYRconR6lKi6mwD1kqsm1mkmr/IIVID9uySS71GgzGt/NSMYvLrhLfuoU9A1Hc8RAjB47w/MQE0YDHSWS/cmcz9m1tccIYDWNTqCvDjD9Khx929oSxAKNSJVUwqGbKdBdMhgS4Czl0q0Ku2+lNxsO+Zu90rKQxU640h8b2Rv1camcJGXlSdBCcmiRo5NBNm7GShqsnjC+foSuVRpdOeEeGt9I3kWDny88RcFcwI4OYbklrNiQW9OLxuXEfeJqe0ZcIWo78e4t5PKUi5b5BCh1dpE8m8OQyFMwKlRdHCKRn8NSVeDCXRjdt6OpjMOZfdFi458ghXBkNAgPERp7HE/YTtqq4KzbR3m1gOB2QgS1OIrsY9BL1e2mYl9GUTkibZXjkF7gqJejto2xVCRazhGKzTs8764RZRg4EqUxMUIr3EcGJ83szaXLlbUznDMK6hTE1RTF0gv6SxiV1Qx4sZunLpqlEAyQyJVxGBV8lg2VXOeZz1EvZyrC7PEXY5yVZchGezlLp7MGIbmU0pTNgVqgWTewXRrhEOIbBq6XwTSXIGSFcJ8foTo/jik0SuQKG9u3BN5Vgsphkut7rDyenqWTKVIBgtBPvkPM0mp09eF5+GalpGPG48+zVNXrBsLFnZjj53PNYPY5Wf8Xsy3SZRexYJ7muPiYivZQzVXbk8kSBLjPPDV6deGeY0qlxfml8Enr7mAnFkaZNNOAh7HfT3RVEr3sdAHsqGr5ShtFU0MnndTneUaRcodvIYRs5MgWTK0WebEkjh+M1RwMex/hraabjZTh8jKFinOlI7xmhUI9hY5tVwtkUUb+XRBtRgtWinVFbp4C3nAdZNiTe5DRDIacHMfbCCNFTLxGbncSbn8Ho6KcQD5IvVxgtOspo9PFfsG3yOEGzSsS06TFyDMgydqyLQjjI8ek8Yb+nOeqjH2cI7au2beHIqKMEfHqertEXcfV3QnSYnnKO7tGjiJLOtt17SJU0jPQshVANYdkII4fHZzMajNFTzhENxpgJxZk9PspYSiebM+ibHmXGHWY21kNPJIDYuo3CsYOEjFmm3VFMu0pvLklAlgnNjCNSaYq9A/TihCBkuYJRdBKHPdEAGI5bTrybuN/NNgxkcpZdvRGOlDQMI4tvNkN52svhUNypa8nxyLSp+hDI/iF0s9q81uGXjiAM5zoGvG6q/U4SNDgzQbakQVojE/NSMGwimVlKpQDJbc5oqljQCyWNmZOQ3LFjzj3MlSt4M2mo6thug0LeyyiDRAFyWaYrVYh30y80gq4yXo8g4PFwSW+UTCBMNfEy+Vya6egt5MoVXLrJwOw0gVwKK9ZDJdJB0OuhApSTGlg1cn4vplUhrI3h8nvxxnuomXn8J44DUL7kKsJ+D9UiyG6n1+rath1/b4RdwGTdQHZYOpeMPEFGt5CRDnqmxrA6eyjt2kY1lcbVFQJ24p4YJ/byC5j5DIVgjL5jBwnm0gx0BHBJg6ItMepG27QdIx0+9Bzugsbk1j2EyhVKBYOgOUlYFgjnHWUWOvY87snjeP1eDJy8lMfI4envb96ngmFjmjZdXjcWUK06H1IdTjo5tpdiQ5StGmEfhPR6HrG3j/5MBmviFL58hqqex5iZIeF3o2dKDCbGcCdncW/Zj8/rxucReLweeiaPE5BOXi8b9DbDnJeUK01vK+x34w77YCqBMTNDOD3DkD+I1t9HwbDprxvuGcCyqwRSM1A3JLumnZc+jw7swLd1G+HMaY8no1fosHQC6RlGy06eJehzU8YJOdqZNLYRwFPIkjarpEXU8fbDPjoCHgwc4xWqt8ns9Bhhf8PYOs9BsuCEqAopnYQnwjCzoKUx7SqD48fo0vy4QrvwFUz63TrSyIHfkT0a8OA33ET9Xna4DSolDZ119kiEEJ9e4jgppfzCGsizobASCbxaCkJxp0c5cgyvVSTgcSM7YqQCnfSPvUh/SSM46iiI9DWvBqAc66Ynl6YvnYB4NwCRzCw7SwkCHhc90QCV5CwJoMfI0fXIj+mfTpGtlCgaVbzVMejvRCZniScThO0yQZeNGfLCzCwhPYfuCWAVDYJeD9EOD0yfxNLS+Ae3MzRyjGoxSTTeT0GECZYKEA1jWFUqEwniAcdrKRg2vtQMndo0vkKGlDeAER8ikM7jKpfRdItXTx2inNQ4sGUQX8DL4OwoHWaJ4pYthP1uuowiwtIZH7GZ/c8Dzmiuzh5SBYtByzEMlYkJ9MTLGNEEnv5+Eq4OyIzQW8gSTk1SsSVes0xUm4F4Pxm/E/IpGDbdo0fpGYgyOzgMeoXLtRP48hks2cVoKE5p7BRbazphv5tt2RQBWXCGJkuLBNBRzNCZnsTlc+HvjeDKptGzaeJuizKOYvUDPUaWwOgRgjOTZAMByvYoMY+LolHFp+dxzewiFozgCfvxmAaULDxex1/x5lK4dJ2w8FHxuPDk/MSxcJWylI4dprRnHx6fF5fbhVVXtF2JlwkaGuUrLj+j7XmOHMI9NoqslAhVDGqlAhXbpNNVBXeF2W1bqWZT5EoVasdHydbDJWGfl0nDRjdtgoAnl8FVLoHwE62HI/WiG8Ou0WNo9E86irNa6iM8NkKwvw+vz6Rfm6C/UMAKhqmanZg+F9HePsLCA5NpClaVGVcHlwIxq4ier0/X0dFJwLCpetwY0TDdo0fYWqniS09T9ggI9zgVTM7iLxgEPC7CPi+pcAcdL71AT8THZKQX065hjJ4kmBOUYt3N6+IeG6Xy/CE8wTDevZei9VnkyhX82RS1ky9TDkWhrw+7MUBDN+kqlwlbBoXxUcqpFFojXHp8lC6zyGB2Co/fQqsF8bqdN0TNSpWBsI9IcpKSroHPgyufhYqBPTNDyJrE9rmbckUKGtFUgkDBUakZlxMajAC9A3uwg15sI0d5xiZY0xkan8YjAti79+A33KClGX/0MUyfhyt+/jgDVhVj95UEXWeGvMrPv0A1ECHdFXI81kuGml6tblbRzSoeWSbnq7QVcloNlipHX2BdGPgtoBu44A1Jg96fP07n8cPI5CzFaJzU1p0EJmfoKjspzQwQKBWISwNGj6ID0U4P0dk8gfQskUgAkjp7JmepBCNY4Q6MWidloGt2HE/hFLKnC0PL4J2aoMs0EUgKNYvZvkG6JxJgFikDqSMvEdVKGIEwXrdFeHKMqj+IqRdwh52wUf/JF8kZFUpGiaBdZdiuEa4Y1EQnXaM/p7dWxt7SQWF8EoBhIFzIIA3HDfeWi7hsC9PyM3XwBXzjo0SAXZcXKBWTeIt5/DPjuPUCfTt24NKmESUdhi/h5f98gr50gmikg1AmS0fIR7E+JDHjjxLU0vj9biLCoGfyBOHJEwCE3IJQKY/LLFItWgSCYfqtKjPbL2U6Z9A9eoRwTsNdqeLT+/GUS/g8LrqMIolMie7cDH6Pi3LfIF6/F5fhZffMODXpIufx0V1M445G8WWm8Kdn0D0BqiEvPlsSzuUI5ZIE4jEik2ME8hnsrjhuy8TrFpj+KBWgMD5JIBij5+hz+JNT2G4PLrcPV7GIYTrXzhY1Jmp+rIkZXCWNgG1ieyzM5CxdHUGqtolIZ3B1pyHow5vP0vPyISf52zJ7bZdZoJadhXAEl1EmlNMoxeK4gkHsoBNWiU+O4Yl307t3J650kqmxUfSuToKjx+lPT+PethUd6DTLdJUzuI8dxoz3U/W66Qr5CE5m6TayRLRRRvUc+bqBK5gVAuUiYWmSsAPUKlX6pUE1l8GVPEU4cQJzeBf6QIVc2SKsF+nX6ul9u4Zt2ujhmNPL1nU6gjohuww29OVHcSHJeDx0mTZDx04SyM8Q9kepFXVynj4Mf5Wwrjujj+wSuuFtXpPOfArdquI1S0QOZglbfhJb9xL1O8OJM/4Ohvv7GUgcZzpv4MllCdQspC9IORjFjnU6SncqQTCXxi4XiGSTUOhhW3mWntlTSF+AHiNP15Gf4c+meDQ8SNDrZlfAi7eogZamMSB8UE/SafZwfNfleDJThCuOB7MlM0XG6wFtlvChArGTR8mNjjMd3wWxAMGXj9IjvBDxOyFbvGyt6XTl8nhMA4/HT1fYz7Ask7d0ikYFv8dNOB4iigfb7yWRKdGTTDPsP84L5Qr2zAzdo0eoBCPQ7+Tvoufp2ySLGhIp5V82fgshosDHgPcB38b52NUFz3TOwEzO4ppMUCsbCCBQrRDKp9C9ASzhp6uo4SvXbW7ATUY3qXb3sSfoxV3ME9JzRLKziIFBsh4X3swsfSWNvNWDoeXwmiXK7k6kYWPYVVw+P9Qk3lKewKlRvFMT4BaUoxFkMEQ1HMWbK+NpuMQePx22ib9gkw9HCXo9hEom+VQasypxVZ2XoNzSYvvIQao1SW1omFo5h551QgGDeoqOYhZpmfgmX8btclEDLF+AslFhMJfGW7Xpf+4xiEbpz0/jKhcpdPUg8jnMnBOucBdydJSSBPIpjEoNTr6MLGaozKQxBnfTZRbpSE9jed1Ey1MENcfAeC2Dii9ApSYJ1ua+tNc/9iLeVIJKTSIqjvsfqhh49BzecJhSYgxX1YvXKNFhFPBn05iz3bh2X4J+KkG5c5ieVIJYLknZMuioJPDYFjV/BK/pIWiZgJtKZzeRUpFQpYybGpbXB1WJL58l4jHIdnRj//iHDHjddHY7nYBQMYfb1DF9QVw+p0foxo3lDxKyTfC5KVlVbEOnemqcXMBNt8eN7RL4PW58eh6RKlDo7afvlTvnjK6R2SxmYgJRsQgBpaCfoNdFsFzE1lJsTc3iDoepWWVG/iVN0OshUilhF110VGp0G1m8L+cxPE6IzIzG6JxO0D+dIN+zBSPaScDrxm9bkE0Tdpm4K1W6cwmCVQtXLovwuAmGLdwVg3I8hh6IELNm6DDykDjBNlvSmZ/BZxtY2SwZVxDRGSNoF+gy85Q8MexqjVBmlkpVQjhM1NSwCzqx8RG63S68lkHNMOkqFqGko9csOmxJ0DbBMIhFnBEQspjHV/MTsE0sy8a0qhTdEXqMHIwfwczN4jd0BjMT2Ed1SpZOMJtFVGr4qjZVyyRYLhCplck/+xwZ3WT75AiiVHS+jZHLslMadLgkptfVPI/H6+YSV4liNM5MfIguu0ZHepqiN0Cto5NasYhVMYgXDaqVEl3SoqrrVAKSfm0Ct9TJBV5Bx/HjxJKzdIV76DItXJaFW9SoJMaoAYSdXE9wdpJATsPt8sDzz+LyuQl095H3ddJh6XTUvf+Ol5/D5/HDwCBdZoHLnvsRnlKRkLQoVUrUEmNECjmGrtoDnDkDwWqz3Oy/ceD3gd8Avg5cLaXMrLlUG4hXBmqUKkXMWgWXWcZllrG8LqxQB11VHVEu4SuXCFQrBLwuBnWLSSAye4zOk0epFgrkgQ5NY9g0SekWBPwEewZwFTT8toUr5MU1axKancXUS1QiMQLVAC63wFXWCZk6rlIeyxfEOztLuOZG+j24CiUolyHgQZgG3tQ0tWCYkG0Qz6eo1mrodfc7VMpSlZAKdOApFgkdO4y7JMj7wlRsx9h4bQuXnsMOhPFIm1ooxGTXENurB/Ai8ZolvHYZYRuIioVdKKCHowRzGSeEYpTxBmOEZiaoaRp2qYzHMqkWnHczOtLTdOgZUrkwbtOmOzuLr6zjCfoh4MF2CWpuH1YwTF9Jo++khlWtYRd1AtUK1ZqjFHxJg5rPR7Bq4c+m6RI+OitF/DUbK28QyOepVAxEcoY4Hiy306P1lwsErTIul8B0+7CrLvwVi3jFRFbKBMtBPLKGBKxoJ1ZRJ1ipEKnauLUKZX8QvCGsQgGvXcFVf1PZVbVp2L+gz8eW9ASucglXJIyrUsFyeSAYJOCJEYxFKBeK9MkyHSWNslmhypnI6UmqhQIlb5CAVcYD1Lw+7LKBu1QkbFcwrE5cGQ2f20cxFqdLWvgiPRTHThFOz1ALhnGHInj1Am6vicuuULMsPOESgUKWzsQIAbOE6fHhN8vEijmidgmXcOEyShAKE9ZzoOfIDwzR6fciTh0nmEvjsS2CEybhek7LcnuoeSPobkFXOUNXUSMTiWNXJXZRx20YuCyDiLCp2GX8+Sr4/VRrzjWsebxIlxu/oeNOT+F2CahYdBU1iICrUiJS70j4bQuvVaZSdkJCffmX8ZczdE1PkI908lLvMB3paVylEp2WgduuIPsHiPT14D9xDJ/HBQPbCMRjVD0CO5MnkE3SUSniNXSEXWHQ7SNt6MhIB12hKl1hP+GoD70YpbOQwmuWKBk6IS1JxF2jquvkI134hoaxJhJ4hcDtdWFksvg7UngLWaqlIgMnjrBFmLjKOpYvgmXXiAW8BEoFKBWoFYvY1So14SJom4QqTl6rr7uPYMg3p634ClmqZhlTLxAoFHCbZTy2BT1b8JRLdJTzzbzkWrNUjuTPgbcD9wP7pJRnNROYEOKrwJuBWSnlFQtsF8BfA28ESsBvSil/Vt92D/B/13f9opTy6/X11wAPAEHgBzhv2K/ZlxtlchY5PYlbS+Kq1KgKF26zhG92Gn+3hGgMG/BWTNweR2F3atNUqhLprSIqFtJTH5RaKEAhj68KNXcMf2IUd7mI7OiAjEbNruDx+TB8Afy2hd8qUwmG8ZhlahJsXxDL46VSMug0y3g8NSYjfRAN0lmcxbYsip4QmDms3Aweo4R0e/HrBedBCgUoW1V8xTwhVxVfPk2v9BNH0lUpY3V1I+y6NvT7saMdlPuGGMxM4HIJcLkJFnO4KxZWJIhLOE2nGnUSxaJq4ynkCVpVp3fuD1CpVPn/2XuXH8uy68zvtx/nfc59xb3xjsyMrMjKYhWraIoU9bCEbjXkltyAoQaMNuCJhw0Y9swjwwYasAFP/Ae04UHD0MDtcU/aDQFt0daDLZLdVBVJKVnJyszKiMyIuBE34j7O++x9PDg3IjOLRUpqiZQNcAExiHPvPWef/VjfWmuvvT4K6M+nDD/6f/DLjMYYRldX7DUFylrQGrcwKEdR2gbKhq0XXfqmF4cMTM3M1OhsiRUa5Wr0fIYwBpTEUR6h46Fthedomus5rVSoyQZtntJ7+jE0YIGgygmLFbXUhI3FLSVCSBopwfPwbINrDdVwyHVvTG82Q7cWW7e0BlokjR9QG4ty1mcZlMIqjXL7SVEAACAASURBVNUOtez6xGlq3LqEXGC1g3FcvOUcmS3Ik5jatLjPn+KtrkAo6uOndIUjOimfPCHyHIqNMZWRBPkS/IDQ1fhFS+U4WLqU6VY7FLFLMJ/hX5+jrqa4KChLZFniVRXlYITxAvwXT9FVSbO8hvCa2rQ0RYmTnrLRCirloMsCx9QIIZBtQ9lKVFPjfu9bOE8D2iJFWkNbVbQUiKZEmQahNInjMayuCdpuHg0vn8FqhdUOCz8kSJf4VNhegliluGmBXa8P2dS0dUHT2g7IpIC6JljNcaoCZzXHK0J0UyNnlxghqeM+4x9+l3g5o417SNOQzM7Yf/Tv6K8uMKal1C62tfhNSXTxgnA+JVMSv7GMpk+QVYWtKlQGoumUtg0iHCXoZ3NqT+NWEk5TqqiHny3xlnOcIicWEhcD1hBkKxbxkLRskGkKosavC3ovjzHT58i6AmvpZ1cIV9Moh3gxw7SG4s4h48vu8EiuPeKywM1TSiE7j8Qr8c+eE16eUcR9qv4Q+fIFtOB5LjJbdbqmrhGOQ+xdUa4ySlvz8vGnJL/+N6wYP0d+kkfy39AxIv73wH/X6X0ABN1m++cnwL+S/40ubfh3f8zn/zHwYP33S8A/BX5p7QX9E+CrdOnG3xZC/Iu1J/RPgX8MfIMOSH4b+Jd/QTv+WpJVFpuXSADZ9YEFUqFhlZLMzjBK4wwHxHWONpZW1oxmZ7iLOZnyoJwBDbQWzxhMkdE2NTLLqABZZCgMwgvwiwobhFSmpc5yHNsimpqsMqyiEWGaIQV4tmF8/imVciiyJW5d4QnNKuxR5SWiqZFCErkKXEUexlDN8bGwWiGLlKEbduEsU5M6LkZp2romFRrt+V0YarFAVCWO5+Kv5oimxs3nyNbiJz2a02NUuqJtDFVRotOMSimMETh1p1BCU6FrMHlGKyRStThF5xk0eg1IyyWN41N7IZRdGEJVIfT79K7OEVWBChM8FLKuQSuEo3GMoTGGpq7JyhJhLUYqmsbgVyW6LenZlgoBLbQtWCQCgdPUKEcjw4imrqjTHFGUNMrByVe4eYpsaqTjENYlTlVQtQ2BaKmAyo9w6wJpGnRVopWi8QLcbIXCYoDGtni2RudZZ3VvbsL1kvj8BEe2tF4IF6fdORRebbqnZU11PQcnoByOkXVFnWZUYYJqVzjpgta21H6If3GGVZpGOyhAZymNkLB+vqpKgmyFMA3UNeHlGfvAKuwhqxK3rBBK48mK1hpMC4oW20Id95HpEp0tcdI5NA1WKty6RDUNDhYlBUIpZFGgTU0gDDJPKYSmlhqnqQnMEtHU1OuUYNnUaFNBXVIoFykFQmmE4xAKg/Uj0iDCrSv82RxjDNYYCsBtGpTWuHFEOD3DL3OauoKmpvJCVJXjpZ0XHAJCCJqra/pX18iqpK46Ra3SFU5rMFWFkhLXdal8j2qyhVPk2Os5BAFRXaDKHNFuUrQGJ0+RRQ5ao5uuz6TM2WhbnF5ArSRBVWIb0GWOvp5hrUUKcKoSU+ZY16MNQtx0gf/4+zRBSFlbYkC1FuE4DKmggjjYIK8ttVA0tqVqLKFpkMagMDhNd/iwqWoQoJcLmqygcX92PIQ/aY9E/nVu3Lbt/y2EuPcTvvI7wO+uPYpvCCEGQogd4O8Cv9e27QxACPF7wG8LIX4f6LVt+8fr678L/EN+ikDSzi7RTx9jsxVNGCMdh7aQCNMgiwKlBA4trRTY3QPcYk79w08IZ1Oo8s5a1Q66zGGd4NEqRasdWsdFhdBTlrrIqB2XJoho8wKqCscahDHI1iKEQDU14ewMP1vhaIm4uiQ0LW4YY+uqc2ltTd1PYApGako/7BZSXXUbx1WO21QIrTor06a0bYtuLY3bWdhunuI7CqkktRb42QKA0vFx1x6J8DxwNLIsaC/OqbWDti1eXSCFQq9WCJXRCIlsO2c8R9Fql9Z1UXmGsk0XB65KqmF3QFCWJTpLEU2DlqDrEjdPSR1NW0t0XSGthKbBSEleG2RekLhVBx62RQGqqbvDlnWN9lyUNRihqYREtBaJRUloTUtbN7hFhqxqLGB8H5UkDNsKXxikVuA4NKZF1wUqWyJ8H21b6tZiW7pMn9pghCR1fLTKka3A1jWtMJRBhGMaJFCcnyNWXZ+iNNb5MZV6o26/REmBl6e0dY3VDk66QBVFp9ARQIGUCq/MkFJgqxy/qTuQQWJtS6Vc3ObNw4uyrtDza9aRJVRrUa2gbS1CCqxUZFGfXLnE5TpzyDSougQhkdYipaH1Aiqlwbb4dYrAB1ehlQQLdWNwTI0ThKAUtE1nnZsGmXf3tZGLFNA6nWdH2+BoyQBLs1jSVDVWSFSZY7wA63pIAZWxuE2NqkrQmrZtibMFfr5CWkujNdJatO3AztsYgd+jqErU5ZTStKAUbdxHtQ2sFri2xRUaW1W4aYqSAul72KKk8QI04NoGK0BKgbIG0TQIDVoIIk+Ra4UsC6znY5WGqkICjqlxbZfQUOB177ruWy0EXrXCVA1FEOMOerDqgJv5HLvKEOkKdx3q0q2ltQbQUGRgLUJqWsCtcoSpqNVfTNT1NyU/q+ywz5M94Plr/x+vr/2k68efc/1zRQjxj+m8F/b397m4uPgrN/D5YJs+EmEtyyAh2hrT/uARqihQTUWLC63FT5fIx3/GajzGCoVbZNgqp7ItwlO0bUvVtrhhgCkNRmucpkJlK0wUd5ZkXaKvr8g8H9EbkBYZom4I8xWNEBhjkMZgPR8RBZDOsbalkQqkoFEaayz+y2Ooa6yUmDBGeBrTWvIwwV8sqPwQLQWqKNDWYGmRtiW4niLXVqyxFrNaocwC1VRY7RLYGlV1mUkVAuME1FEPm2XgOKBUp2T8kFx7VFJjtENgDBgDZYFfZti2wRqDsBaLQAC6rlg6AW0rEULh2xJrLYYW2oaoqciUogwiNAaTp9QtFEGPsK4RtNQIDN3ibpRDoRwiqVAIJC21kJQofKHI/BijHZyswigFSmOkpW0tVTJC7x3gnb5ACrB+0AFcVWKtoWk1eW+EVgKTFYgiRdCiHU0lNW1RYusaPAeEpmkFwfIaXebk8YC5F6OcmqjMKU1DpT10GPAnS8HiB895MS9xXp6QtIq2MbjlFUYIrNZoOmvUtS1KSoyUtLTQ1LRtC1IhgVY72KbByXOkNbiBSxUltFlO1HTtdVdzbJBghUTo9fi1BmMAWgQt1lp6+QK/zDBSrtcVoBTCWqSxaNvQOJqmKPGaikaCdQLqxtACKojQaUVdlRjtUEpFYSXROiJduz7paEKQz5ECXAnSgnAdKu3C5Qzd1FilaaVCNiXKdqG0wvGItUZqiWlbhLEISRcypfNcpKmpHQ8jBGqxoNUO/ipFCUXruYggoE0SyJadJx0EVOH6gG9jsFVJVCyhyKnTBcVgg0YqNF34snVdjJQIa1BNxWJ6SZ51iQtueY2tKzIvxm1rHCmBFtMYrnsbeKbGr0saz0e2Foui0WC17sBrYZF1TZ5mmLJCNjUGQbtY0FZdOF1aQ2tbRGNQssU4Dk2eIdIUL8/Jf/9fc/Gf/NZfWff9VeVvE0g+j9Kn/fe4/rnStu3/Sre/w1e/+tV2PB7/lRt4d2tI5iiM1gQY3NWSSgislJ11IVu052HbFpmlqCxAuApTVQgErqOofQ+TGdq6gdpQBxF1lNBbXFJJgW0711tJiVQC4Yd4rqYkZBn7uI5CLua46/0Lp7bIOu9cfeXgVgW2aVDWQL6iVj1wHLQxjBzwsGRFjhop7HiCl62QeQpaoRvTWfcIhFA4gx5Nf4jpj3GfPyXIU3SVYf0A00tgtQIpaZVGSonTVFhTI6VESIEwAkcr6qrpfjfapRgOCS/P8KsCYQ21bbGCTuk5GiMVlQVZ11jfR2XdPoWSgtZaWB9YVFLghQFaS8xqReMGSCmxUmOlQFmLQdD4EVY7uFWFS4snIHc9ChnQmJbCj3CbEr9Ku81rL0CXBW5Zdm2pa66diEi5RC1o38M0DVZrWtdDS4FMF9jBkNmD91m4Ee9f/BDOXpA7MXVtEIXAKgdrW5RSOLbEcxza0ZDA1QhPo1cW0xikVAx3tri7NWQ4HPLHz1+iHx3z1vWS2BqaIKJyXKJ0jmyh8kPK/oDAUdjLbg9HOgphDEKAWCt8tKb2uiRVi8AgWEz2UNYQVDktAlvVVF6A3zaoqkQqSakjpOsiXZf+7g72xQktdGNeV1gLVeDjWQttiw26fQsX04Ub2xYf2wFCU6HwUa1BVxWiWGH8ENnr4+YStAZHE7UN5XCMpxXx9SXUXftdLKWpaa3BW6+7QnsYqfBNxd7TP0OVBaIxICQIEIjO2wRaKWltgxESqz0wJVIIPEeB40Mc00iBozSR51JHCaXSxOkcqySV0jh1iahLqGuUVMSBi7QlWIUrNTQNGQKDwrWGvO36OihzrONQC4XXFORRDxP3Sa7OEbbBNQ1RXeBai2tryGuaLMU4Hk3SR2zu0CyW6LJE5WnnPdPN89YajHYRUiBNTYvA+CFCK0Sv3+mCcl2jrKn599F9f1X5a4Wv/ppyDBy89v8+8OIvuL7/Odd/qtI6XZxR5xlZ1GPV6w5HuXU3UMJxsK6HdV+5kbKqaJWCKEYHPgQhTdKH/gA/CgmyFdX1Nda+wsFSOlR+SOBIwsUMZx3jLUabXB+8RT4YU4QxpXYwdUPbGBqpKdyAwgsovRBjW8raUqBI4wH1F79MabrwS7GxifEC3Crvwg5hBJ4LsrPijXbIoj4mTYledCeSdZF1mVye27npQtAIiW2hkl1oy8tTSu1yuXe/e8emIetvkCcDjOux2L6DlgJd5p0Vaw3adiE728LKi8msQGYrdLbsNm6BRmpq5UAYInwf67pU2iEdTrBxghRgHJdmveltVRc7lK6LGW1gg7D7TduiWkuIRXoetRdghOruHYRIzwU/gP4A63ZZRMVsjr2eY1uwUUI9miC0RmG7sdYOpXZZFTVycU2uXdowQklBUKaUfshq7x6NH5L1R5RRD5RCXs8QZ6cY29IqRaMdyFbkH37YheJek3xzl2xjCz0ckOxsQRjRKAeruzBbisL2+tjBiDbpo91uHtaOB0GAiftUjkctJLKuuzTfKKLe3Ma6bpescDtfS4QxuI6DdF2aMKJJehR1Q5P0ELoDESMkTRDRjCa0SoNef9+2NEKiXQfHc7FRAkFI2xhkXaOHQ3S/h3Xcbr5u7pIPxtSe390HqI2lbAxMJmR37pMrlzLLu/Cc0uC6WKWQUiBdtxuHoqCtShqpqByPMohoPJ9WSox2EIDRLquov37GOstO0M1JKfDLDLmcsyobVq0kHU4wowkmjPEDD99RHSApjUEir2ddpqTW6CSGpEcWD7rQcFmg51c0nk9jW2xZoVuLspa6N8SLQ9ooZjnaAt/Hx4CUmHXSSSslqrXEbU20uMRdXnVr1+sKQ8qmA3w3WyGtoeoPIe7Ryi6EaaTm2ovJKoORksLxmLU/G1/hb9Mj+RfAfy2E+D/oNtvnbdu+FEL8K+B/EkIM19/7+8B/27btTAixFEL8MvBvgP+CjiflpyZn4QiBwnNcnF5CHLk8ffAOK0BnS6R2cKoUXRWYIMQmfYTrIcsCoXQX8llL67iwswPnF6jVvJvsUlGXFc36dLQsK6Tn4W5tYafXXZG4nX3qXkxvfk0mHJrIxQ98mnVISTZdbF/GXe2oqqnJ+iMIAuLVskvTDULquMdmdk0YBWRVQ9102THCCzBSoaTAbO8i0iXV2VmnHPw3uR2k7yOBqrEY26LSFbJ8dfK21g5tXaPqCqe1uNkKJ19R9Tq2qWq9kCsExvUQrodxXOZRTJTO0QocY6ikREuBYxsqqTpr0PORQUDRWOTWDpdBnxfRmLc++RDKrAMjWWPvPaAJY3qzf0sZ9ZCeR68uMGXLIuxRbWwymh7jZHPKQbc3o7MuIbERmlV/zMrx6TkuNAVLFCaKid0rhGmwdU0jNeVgAxt1ZzHqNiIvclgz4hntkDk+gXLwmgq3yMAPcEcjnPMZSA929lilFdYaiHpcrYs7ztKKeLRJVVzj2ZbaWFzA0ZJUR2TCIcxXuE2N9Txy7aKqEuF4zOMhMgwJ0zmldiGMsXlB1DbodVah8QLkYIhZLpFCQBzDaABVieNprPRZbB3c9oluLGiXxg3JvQBXK9T9IwKTI7KUZneP7PIafz4Daam1Q60cfNelCWNWyZBSCTypsSgqP0ZdXnQb7HEPW1YUykXdfQunzkmBlRPQyz+ldT3wfSw+zdZON8nyHBsEXcZjmoIfUA3HZHmNDAI8YWA5pwoS/CqHokAHAaaukKah6Q/R8yuUgKvxDly8xF33jZWGeryNP7+g1g404NAgpAApEY7GlRJcFxwHZ2uLRrtcLFqGp88Y5CWtbXCVxGqHSjnkGzv05xe4StAMxtg8pyobFGC2djGNRdmGNowweYmqK5RWDN77AufHxzR5SROuuVDCmMa2XTKHq3EcRetEmMZ2wHX3PjroPPJWCTSgnf+fA4kQ4p/TbZyPhRDHdJlYDkDbtv8LXdbVP6CriJ7RHXZkDRj/I/DN9a3+h5uNd+C/5FX677/kp5yxNV2WbLgeYp1ZRBQjd4+Ql+e46RytJe5om+r0tEsTLQuqxiCVix0McXe3GHkO02VJsa6462tJJSW14956N+k6lTZqG7LxDir0cFE0dZcddDbaw8lXiIsprZK0dUsjK1LHR0mJkgJ3ZwuZr/AxNHFEsbHJfLRB7D2FsqRXZoTDHs3lGW4/ookH5Ofnt1YZTcM86bPd97nMUzLHR57X2KYmGI1o10BpXZeLUmCDgOSqO81svIDrcIQuC+K67trjesSbG3wy3u0OC9YV1rTIpqZ2fDLHwzevNoBlFOLFIXWWU1Yei6SPd3FOgMHu3kW/9QANFKuCJltxZRyqpqVQDkpqnGDtmQD+5TmubcijDcR4k3x5xXXesoo3uHd3l2F1iS0XXA02aMIY5+w5A2FQyuXq7kPc6YyqNiyVy9wfUAURVT8lXl7h2gYO7iD37+KfnEF/wMILcbYOMMECyoZsskM73iWYz2jDkEIJpKvxRyPqZclqskO0u0n+8RO8xQxZ5BQvX5KuC/ztffjH9K9PKYcbhL5DkacESR8/CKnyCtnWVMMRzhfew04vuwNtF6fkwsXTEtuUKD9AFjlWCax0sYMR6WCXFLAPvkD8nW/QXF0z623xTk/AbIZ0FIOkj9rdYv541WWazWZUQuJubWK2DqizFQ6wcEM8qcn7Y/J5jg+4nocZjLAbE7IsZbWuO+woSRMN4WpG6CqsMThVBb0+5Dmhu07XXnQHW8vJPmm6RM9+iJUaM9og3NomanLyucJuTLCXUzg97byc0SYr4RH7Cn19jqkrstGELM+xOmOxf4Sbp/QuX+JoF5QDLayckLA/RPX7XWWGiykqXSLyDOu4GM9DTVcgJDYKaFuolMLt9wl6Mbl2WQ02oc1o2230SlE4Hk0Qono9Gu2RJmN6bYWrJDJPkUGA2Bzi1wU0FdXd+wRNRVsWnTHnuJyKgJPvPWMY9iAEGQRdeMtxQTvYMse6AVmygSpznKCrIACwcgNcLySMXGS6pPmbV4ufKz81IGnb9j//Cz5vgf/qx3z2z4B/9jnXvwX8yJmUn5ZUO/s8ef8/ZP87f0Dc1nh0ZRouNzZxXx7jVjlBbwuWEflySd0Ymv4QbAv9zqFKy5pPnR5Os8LTASIeoHo5JutAJ4v6nO8copYL9ucv8foDXnghZdESUJA6IUw2Scsl82hMv8oIrk64HO0z7U/YPX+G10+I7mxTnjyDvMTVnfXEck4W9VE6xyR9lr4mSVcQxZiNHVarAvKcIBrQqzJ6VQrrkvf55i52tUKUDfX2HeI6p0l6rJIhL3NF6Era5fltXxVhQt4fdSVFBhugFYyGiCef4kxPGfQCZkbQzAtUP6Tqj8jTlEx71KYlkZpyY5NmQyDPT8l0AFGC09ZU2kOUDeFkg+/33yJazJDXl+jBBqvkS6zm12xmM8J0Tk0XxusFLkJ7ZNdzVpu7BOen3Fuekvh3ybULQcQw9Ghjn+tsSKMVgauJXIerMIGwy25/4Y9wpKTxeszDPhPZ0Dg+cz8mmNRErsOxCMg3D3lwfQxnL6hVV7yw8XzCJEE2JTZJENu7KOsQbO/T+pp0MOmApPzRsvtNGJO6MXFbUG3v4/gasZizxCWvLIEr+cKXP8A+PSZ9ec7FZI8Lf8CD62PKZecB1vNrtPYIaLBhROBKnNWCarai1C5NmLB48B5Lr0bqY8R4yKpouPISMidk4/oSbWpq7ZBVhqK2+GFMf7KB8DRitSScbNCevsD1PYK9LeiPWDkhUxzkqiAc9CmjhLxu8JvuPAcAwxHiV3+d4uU5g+MfYgBz9xB5NUPPr7oqD36MvxkTOZqwyUk2ujW1ihMsIJ78EKs0y7tHLLwevTanPnqL8ulT3OUSKwX5cET0a7/G8VWO84f/JzQlQmqCtmHclviOhqxTwvlbX+DkwZe498f/qgt9jSc0VYWoK3TgowDrehSNvS2RIlZLfMdBJwlKlljvFUmXjGN6gUZVEabpvJ1y/x7LssaZnuDPL6BuKBtDmdeYeETYi4i0Ik3TzmsLQ3SSMAi6MO9xJXGyAjyPYmMTna26TDY/Q423CYCyP8DGPjJ9k6Ttpyl/m6Gt/8/L3mrKqaeoDu5RVTkrJ6QfOlwdHmFePKOtMnLtEr19xPWLLivMdzQXYY8sHFEHXUn3qtIs3Ih+OUMvl+jG0oQR19EYt8zpVxkfP/gS8XzM4TjiAr87KzEtKMKExNeo8QZLEdN//KcAOGVGkC1JBxMuw4TVywXbaUHjh5x94StEL56hsxWrXrfRZvwYNT1Bz1eE27vY4QgZx1igjGPqzYeYu4esvvchNojwJpuEszOKOMZ/+wHVxoTl8xdcn5wivIrYFuRCo30Hsb0LvQ1SpwsRxP0Bk77fldIHGu2ShwlXToywirY/4rK/RWOvuBh3JUxcs2D43tsc+xssPvouRW1YOT7kV9SOj/W6VObI02QPv0g4OyfxNQ3QnJ1RXGp0kpCPtuiVHcf5JzKhLmqC0QbK0fROP8Us5mSOTxx0qiAtG5z+AG+yQVylWFeTxwM+vfcBcbbAGMsq7FH3NnGUJCwvcVZd2Ec+fJcgcnk66977wfUx/uYYDu9zLAKukzHW1dRvvYd8+C5Hz75L2JR4vsYOR4z3GvyrY1wsdyjI13WRxrGHezol3b3L1QdfZpJ4DLMZl0+P8bIapufkle1ohP0+Tr+mDPpk4YimTbt+OXrIyUXK6PvfJqEmHg+JrMPV7l3KwRj/8hxMS97fYHl9jFs1tEWDWC3R7hVV1MMsAmzUcuokNDv75PGAYblkOBwh6arNyv1t5o9/SNP00FtbCCfAbB3gvjxHXnaGhkm68ipxlTNIxlRXV9Rhj2pjgvrBx8C6pHtWo4uGtKpvFbUNIuxoCOsaVmJ7F7N1gP7eh90+k20xnz6HB+8R7WyyfHlOOtoijycdFcHufbzD99FXf0I93sKrM+b9EbKtCFwNa2Ar6gaVLtk+fdo955d/jd6XP+D8f//n3QHFL/8C6tkT5PkpIs/I+yOWGztk00uKrCNU0+Mx0XBCe3FFWUfkm7tErkNw+YJ5kFDt3+s4dfwYk/Qxkw0elw7bf/YtXk/UTR2fJl/vkW505U1knjKKYnxbko03qf0QFgvo9ZD7d6lml7hJH/f0GLu5zcXhETYcMHzvC389JfiXlJ8DyU+Q6bLkwh/gxldU7gbqi+8hE49kWWLGm6R5Sj3YJFE1+q3OWgrrDEkHHMpz8O/usRodUp+cYJ5/D5mnLCpDrj2uR9t4WpHHA86iEQc2hTWbROhogo0RqVA0wPyt93ihOktznG9Rn03JpU/64Atko02eFg3OxRnBao65uOR8Yx/EJb6jGOxtkxzsoikwnuZ664DjqwzlhFT726jxBtZzqP3+LZtbPRpTbe93lve7H9xym1yfnHK+fYhNfAaLC0gS/LcfEF1lmCXUQUy5fw85iWiPT3HUGelkj8EX3+aydJDihwzbkl6VUnmaC8AkPbRpefnslPmdHmcf/ErH8vj4u5xv/AKr4SbvzT7BDkOy4FXdoGJjiy+IFcvc50n/Sx2hUWWJkj7+/h1WzibpR9/F6Q9JfE0778C+/PXfZDXZZuf5I3ygefcDKuD6+x8yz2qKWRdiqXpDXjox6XCTo/NPqIBH999GXZ7j+wM+WFdcXYx3ugJ5s+9z2R8z277XheGAcZsj+0OyyTY8+y6Z9nCGIzbu7dO8+wHm4+/TXF0gFnPUyXN+8f2HNEf7zKuUcjDmKp4wOz1hGkVM3v2g4+P71re5enrC1VXe8W6EHc/GSTxhK5vRDx3s1g799Bny8AhbLhGhQ+QlNFlNks/pBS4ELvvDgIvjEuIJd+8eoJ49oekPKURBUW1S9CxX1iMBkoNdDsYRG4nHJbC4AptW5Ju7yP27HO9sMs9qko0JTVZjNzap4x7L7XtcvfsLqA++zNbsCfzJv6GNEyaJx9Vkg4v9e3z8ld9gbzXFmZ7y4ipnF3DCmObwCHWwyxVdcdTl8xeQ1fTf+4CXtaJqLMX+PTbu7WMjl+PTJXgB+s4WedEw2NvjaL/P8fc1ajhgVvUoakN08QxbNQRhhMi6UG7jBpyLAH9yj7vDEe30vAsfDkeMf+Vr2MhleX5KG4SIpPNYy/17hMCmKsiurxGDAb37b/H88TFUdUeDsLuNcjyCScQsdMj9Ps3sguZgl2dRx62zeXmMzlY0WpEah0bX6DDBOTyiHo1xXj5hPhxxdpWhnj1BDTtdcOn12O75BMs5F36MO+rK+8vRmHnQJ3n3FR34T1N+DiR/CVneecDjcMRhGPGV93Y4ezKj2t5HLee4gOr1ML1t5NUMUWns7tvkIkL7mtG7R4yWJWdAXjfE/T6Xe/d5sXnI3zHL1gAAIABJREFUxbIgXs4wvRGBo8j7G5wPAz4VCdvVd3FdF1Yl3vUFZ+GIXvGSwf27HP5H/yn/1+99C7smj7phTEgHEyK3A44nIsL54WOKt45IDnaRl1PscIQ9eojd2oE//AZV1ON4crcjwmrTjjfha1+lWfPNzw4fUOY11588Q0Yu8nLKYG+bcOMOzM4pNjZxRxts39vnO/k5c5viJQGToM/9yYgBMH/8Q15EE6rf+kc8O10iZYS9+JSL3fushpuosxdEk7uocUR6uSA8OGAl+gjA+hoBpPOcZyJmd2OjIwjLa7LRJu7WHmerKbP9mCzvFq13fQHFims/AQNqvMFF0GXtjDa3MXcPGdy/y6PTJYvTJRv39nna9gjOTvjCxoRkA/bLmmci5iwckqYVvqs43n8AwHsyI93YxNnb4yu/fJfns4y977ygbxyao4ekvPLE6G9wVVzjTs9xRqfs/f3f4PpPH3el39eEVOn77zH/1rdZvDinXpbs0pEUFeWPj277v/MPufyDP6V9+QKR12/QKgPYjQlm74DR2UvGKuL87lc43TvA+bd/AtkL6tGYNoxo44Snh++jr3ISX/Py4ICsDTvyrqBG9nx68ys+CbYZ9n36QQdQZ8DJ1x4Ca3I2v08QuTRbO9SfdBl/9VqRLYsG1v0xSTzGlcfzuPOyN/I5NnS4WieQTRKP8k8+Zreq8SabMNnk/rudor18evzGO86zGnP3kGOxJqRag3rkKdLSdHS7o032ju5zMApZBg7OZJPHG4e0L19w56Jrp93dp1elXFsP2R8SlTVm2d3fbkyQwy7t/mAU8ulkk0+dHlXU4+7WQVfCHeiHDpOtEX/U75JKH8yOO6+wP6Q52EWGDpyevdH+7OEX8d85gkdT8v4GbpuTxj1+KCKaiwsYBiy/8C6/+ff+A07nBVezC/p0RF1Lf0Dia/TjR+TCUI/GHE4iXL+Pflww2N/m3733VcKXx/zi4We5T3868nMg+QkySTxme3sUAGtWw+ezjO+0Pbb9mPFyTv0LX2P0y+/x5Bvf64BkMCA7+iLZRQqz81t+7mqUsjiOkXvbTPffJstrru5MMBcvSXyNiSecBw5b7ZzwB08o+iOCvQ3U42cMixVlNqO/scUocnk+y5jnNW3ZsPX2IW48gT9/TLF/j7iMmdGFbJrxPsXOAzbSnDGdgrFbO5i9A8xvTYg+/Bbba6a7el3c7cLvg99xWp8dz+HT59STbeyaaCk7GLESnWKe+V+hHziInV2yM6iXBf3Ep9jY4rvGZebtEQ+28YFHpx0Z0bvvvIW+jDmrA67jMcRjdpyUWTpjUdRsronE8q09ZnSkQHFaIaYv2Hx7l+cmYf7RI/qBQ5x4sOqYEkeRizNd4uQOFPADE1BPtlGTba7iCW/92R+Q+JpnWU22LImP7rMCBtsJk9jjz/88x7cph+OIxHOIhOZ6Y4csbhiuOhbByHPY2t9itgba12We1zx//wPccISTdu0HmJ+ddMyDIuHhvOgs5zWL48Eo5M8BFUXUQcyHx3MW/+R/5u7j71AmfZYH73VUy+xRJR7bqgt3fHtecBJP6A0b9gOHrWyGfnoMh6/a83cebvJ4dbguvdJJ/QtfIwsf40xPaeOEtteNY7oGLb2zTwiM1nNheV3jVjXRQFNsbBG+xos+bbuQUHWRcjjZZvtXP+CbT2a4k5LN/T7PTULx+BOSyzPCOwdU67Fytrdo795jcbVArkomd3b5eC47zvaHB0Q7m/DynMXDLzJLK5Ll6g2O9no0pp5sdzz2ec2qaNjKZl0/LUvKwZi8+FEQHty/i4zczhDZ2eX66Zhh6DL6B78NgDcv8IDlH3ZMof3QYfz2Iec/6LygG0rfxfYdLvw+7sEDtjZmsKZicLa3qF3FLK1wZhcMqTsqgMspl7v3qawmoJur1ne5gZVR5JKNNjn1Nc3ZGRcywolqorSjmZZ37jL9ow9hDTx3T57z6HTJ1dYel7LH8SylOXyfdwYV/OAJzdFD5NuHTK9LghtK35+B/BxIfoJs9/01Pew++bJkdnayps50UR98Gf95yEVa8XyWYbd2+KiNebid8KX3H3Ly4UsyoL++12qyy8f9PYJKk8YTRlsu/bQi9jWH44iV6PgF4qbjYchHm9itTfYvp9xXitxxbhfy6bzorND7d/nSr35A9WTGJR0IzHo+9WhMdJVD4OD4mnq8jc1cls9fUKcVA7qSADJymec1+dYe9cqFbNZ5K2uZJB5yHHFGBwT+ZUqxEfPwKIHthzw67QikkrhPP5gSJl3464ZKGDqPIC0N87xmWdQcr3KOgIfbPU7iLr34YSyZLl1yGTAKDHXoUiXeLd/11x9N3xgTdxwxSTxM30etOiXz6HTJ9emS7X6f/mjMWTiC9UKa7HR0tu1gwIXbp0grXhVs72Qx3qEwAbCiHo3ZmmzzQTzoqGBVcDuG1XYCjz9hllacxBNkDHznxW3b1Mq75eSepRVzNWBv1HkN02W3t2K3dhB9H+/wkOnuPeoXF6i4R+Lrrr5ZY7hOxiTrxIdJ4rHd9zm5If2dF/QuXqLOXjD377BFZ4k3gUO/dd5QvK/LR8dzqouUt+k2tt/70hH5Zp/fH24ifM3wte8O7t9lOtlmnM0wa4bL1+fFaK+jlF4FP/q803kBccLD7YRp5FLt7LPd99k9fIi7mjJZLMjWLJbRe0ekUwWnS07nBXdZ0z2/c8QI2FxN3wDDwf27nMQT5h89Ylk0pGXNxbK4nYvAbb9tjSO+tLbIt/s+9A95brrxU8MBYS+4ve902QGjv7PD66kP27/5dzn7ow85nRds7+zSvOuQxhOq7YRN1c1Py5pfvXI7HvbjjuSuV6WI0uH74UNMsGBvMrl9ziTx4OQ5wdmSJ/GkY26sXhD0NINf+Rpb2Yx6bTTezKVqXrBHBz5V4hEMA9I166W8cxc7L9jKZhyMQrhuOIknnMSTnwlH4s+B5C+Q2Vrp3FiGZq/PZF4w2j+ifOeIdF5wOp/eTpDpssSsB3812WVvO+E3Dkd888mMY1cTefqWMezhdoLbRp2Fu7W2cKfX7A8D5oHDSTyhaUPkxYzRkXvrFXxnWfLp6YLtfsDzWcZ0WXISjxmUNQ97DfVkGyesqOkoaCeJx0R5yNDhYt3GWVrxcHuHX3rb5yTuczr3sKs3FcJ232fL6WPr9fXoAL1WCjdc1DdyY+HXJ11K8M2iHv1n/wgzL+DJjP1hyGj/iOLs5PZ3o8hlu88bC2xzp89J7L/xnRm7nMQTDkYhB7/6AfcnMZ9MV7x4dStWw02KYUB/bVHnW3u8v9/vFtbh72A/fUbxaEp8dJ/tvs9HW914btNxxY9WGTdxwhtlOUk83Jfd+Q735p0il1la8c0nM7bXIZ/PSvX0aQdWk134nLrZp/MCOV3RTLYwiwXy5UtCEdG7fImXzknzlKuioVn3y8143Ejia8LEJwN+IBK2BjWjyCU9+iov+j5quurCMS+7vq3ma2rkyS5nq2nH/PmaLIuG6bKzYrcE/OLhiG9UETUwAqp5wcmyhOX0DeCY5zWPTpeYva7jqp19dg9H7AL202XHKd/3+TsP13tbkxh1dcX0UQcOYpYxXb7ydJq2My726cJJdgViZ5cm7ADB7B3Aa3Mv8hwCVzPPawb37xKG5W07Rn3/dp680YdfeZdlNqOJXOSdjkysMjPcl8eMIpezD75M8xlwfHTavcsocsnX80DeuUsdT1AnXUWnG7CtR2OuZiAWcxaznCcXKVsbm0yXBcvnL8gefpHR3gHq5DnzvKY/cRi9c8TH4x2WRc0HWwMGosv+er1s/HRZsie6ufmvT5e0Vzl5bXhykfJ8ljFblpCMeHB4CJ98/KOT7qcoPweSnyDyzl3ya/fWep2lVWeZ9P1OOa1l+rKbvDfK82AUcjAK+eaTV1wA232fk51dnMB5Iyxy4+28LvO85tMq5z4dEFzlCVk4YvLad8zeHdxxdKvQP9gfMBuGnK2msA6t3LjSHSD4XH7m/czeAXIUwiy7tRhv5IZk6fG3v8f0eP7GexyMQk7nBaPIJTg7YbaaMo0nXF7lbK3bfLOoRuvfHI6jH+nfUeR2Xs+dHSozY3FyyXRp+Mp7IUeHm3wyXfF8lv1YC/tm8UwSj/joPuJ0Cesw1DyvcddjcSPyzl1ik9wq5Jv7vg6KNyGTwWvvS/+Ial4Qv/b/yfH8dtxugOWj4zl7q5KH269SQE/ihNWqG3Ob7HRKPZ7ctuEyrQmbrk5U9OIp/rpQ38rxu32gz33zTrLRJu69e1RPn95ygb/eN4wmqF+acPLhS1grQQBWXYze3d9hbzZ9w/N5lE4Yta88X+A2zn7jnZvXAO3GWLr5/MeN1eui79zBXrdMlyXPTcIsXd7OhZug1O2YxN2sV/Pnb9xjNdllsd4f8t45IlyP+02IC16N/f1JDJP3OkCZZRyMQl7cv9udv1nLzbjeAOf2aspnpdpZF9ZYg+MtEdnoIc9nGcdXz0h8h3rSpay05ZLGS3Dv3aPXE/D4z6lHY0bvHN321+E46rzctTy56MxMs/equMck8ZgkHiexd2uUjCKX2cEdZqcLoqLzSkbvHN0Cz/v7/R8x9n6a8nMg+QlyM1F+3AJ5Pss4nRdM15N973Mm33RZ8s0nM6bLktN5dwr8K/v9H/kevLJozsIRJ07CfbidlDfPny5LiDsL+v39Ph8dz98IOxx95V1O5wX5suxIbTI4TTy2XjOat7IZWwK21wvtYBTevuvXH52/8e6na0v2s3sC7stjqouUmyfPRMVivEOy0vThdo/gs+83ilzio/u3i/11K1sd3MFGn7/JHJydoMTidtHaT5+h5sWtYrvJNnNfzm+9is/K/Ul8Gyo4nRdveJGwtmLfOWJw8nwdSuK2f16Xm/8/Op6/Ec6aJB7vD/pv3H/allxe5Rxf5ewPu1DK62G1YmsXHUTE1xfIdNVV3wXyMGE12eVh4t0q8hvDZLvvs7e+dhJ7vOaU/YjcvO/r8/dmz+XG3t17/9XGeXyRwmdA/2YufPSaQQHdfF+tlkDC6bxg8lpbb2SSeOyu++vGM+h9pmreKHLfUHyH44j6MxSxNwaMd7jJ1x91e4/914yyG0/o/iTm64/Ob+91M59vPKLbub5+59dZKV8Xu3r17s/fOeKjR1NGn9EDX390ztcfTXm43Rknie90RtSqgnDEwy8nb3iUn5Wbdft6C24iCDdz7PksewUqnwMMkadJfOcNA/eT6epnCiLwcyD5S8ksrW7j9Z+VmzDRaB0L/nEySTye7N3hBHCP57eL7mah/+LhqLP4pt1k6q1P05/EE4L1vobZO+DkeM6Ti5R+4NxOlpvFdOMx3TxvcP/uOuw1YfdwhD2eM79IO3c6cNjmxy+km0Vv9g6I4+IWMF5Xqv2gy6Ca5zUPj5Ku/Z8DpjdgdCM34HET8vhkumK777NaOUDT9cnse91vTWet3SzWz8r/2975B9d11Qf+85Vi/ZYVK1YiW5JtuU7MhjrJYscLy7IBCh1Kd0Nb2g0t5dd2t9PSLMy2MENnu93dMrS0pVCmMJ1SfkzKlDZtGKZpmxZKliRLISEJTeKExETYsWXZip+RefKT9cvSd/+45zydd3Xfe/fp/dCT9P3MaPTefefee849P77n++Oc6zuPf46+DvpmFiJnfEIZfafzz9+b6rwZZmlopKgQgZUJhGdpaIRZzdJDpPWMT11m6swY4AbwlxyINJ1MZM454Mrh8zw3vI/rzn6Pq86fRubn89uTem3Sp4uXf7Cvg6WeFYHn/UZxvI8lzHPm0jzjrQuM9HcxsdDBg8czDM0sMOr8T+NTl8mNZejvbuNEZqWdhs9kOed8iERCCIBAkHjNfAIKzEv+mU2cyUJsoL1448sYDoSKL3PLnr20x+qxv7uNN960a1V9xOtsMju36v7FyLeVgZeu+m2gt507ju7JX8sL9syl+bx2EQkap2GcP15w/vJ1u1iMlbe/u41+VzcDzi8YChHfVj15oXImWzC5yyQIrMylecanLhft47XEBEkZBnrb82aLyK/Qm5/d+M7lnctxrWWwryPfCCazc1yaW6S3I1INfCX7dMunTzGUyzKxbx+zAN9fmf3NXjfE7kO78tfxgssHA8TJNzzXd9rOnWG59RJTMwsroakB8U6WNHiG7B/o4WzvilPZ29enZhbQuSv5WdWtboZ497dOr9Jo4ubBkf4upqfbWNq+kwnI250zGpk8vI365j2F12hJyGu/c9ROZucKzIulyhbWU5i3sBOGg2mpa470d9EamFi8aSLTvS//zMKBr7u9lW69wralBZa2tTHd1cfy9f8qP3mJD6qR76DQvNTXua1Ai5rMzuXPC9uIryeASfHmo+6ofmJjbTHhHS/3RE8Hjx+byfuKxp35qH3UhZHF2tepqVmgvcAMlURYJ+F1p54bY2hmgQMHbyxIX6x+kvCay/jU5RX/TYlyhu3Xn3t8cmXluH+mYZ292NVPZinSYnb3CbOxhebenLfEitUj3pbiE6XJ7Fw+AvL2W3YDUTBKfKI22NdB27kztE5MQ5ny1QITJBXgB29vu4+rj8U6nZ/tDu9YaZBDuQxDuSjyx3NouI+b9+zivqfOMZ0w4JcjO7uYn1VD4QAymZ3j6v17WXQDyRzJHc+XKWzQxQbWsHN5ITo1socDBwdWnR+/TiWdfrCvg6nu4t4Cf607ju7JmzWK1YXPf9L9v+WEju+gHl/fYd7D6/trPnpyKn/vVqJnMgEFEUVxru1to619G1yeoX02x2zP1SwODufNYCG+bm47eG2Bzb8ncxbcIN46Mc4Tuj0/6QnxZqf7njpXINAmlwrbcd53duw4k9m5fPBI54sTtOXaCrSOOP65Jgng+CDpB8Sjo/1lJy8hvq35c+I+Sf89NGUWy4MnadbuJ1j7B3rK+n58/cY11YXsHEO5DDKzAEQmvAmcmbBngP6hvlRl92m8bzLER0RCJNB8exui+JhUa0yQlCE+64+bB7xzLHS0w2rzB6xoLvGGMNLfxVDbyqxjoLedyzs6izaCg4O9+UEh7iguRWiGC23XxUhSrUP8/cKoJV+2cOAtRxqhksaJCyuDCLDq+SfdtxJbcjyfca3KM5mdY3BohJH+LpZcVF0pTko3Vy8u0rZ0hWVdJretg56BwZLnhPjnP5md4wcnTtHZ1U+/81GFQtxHMHl/ROZcpPUODkUz2FaZzms141OX/Us9AadlDfbmfX6DfR0su92Obzt4LWOPF0auhYOwxz/vLlngpft25NtK2L68udP7NjzhBCepPYZCIw3ej1BuoA1Nz/H8ZGOLQeMmRH/9wdYOxiZy3PCKQ+wf6KHF3bvt3BmGcgvQvzdfNl9HSeY6H+QC0Rjg009m57j/2WhlSn93W35M6nSm7kasbTdBkpKwkRxzPo4kQhUUoo4Rd4iXY7Cvgy7pzjccbwZLGrTi6v9Ab3ve7zAVS9v/kij6iHOFq4SLzdZaJ8ZpuzQPfQcSB9FwgOx36nTY4UvZZpN+393Xzo4dKwJ1sK+DgaXyAiRJGPiBJfytkoEmKb/hNeL5DzUd75QOn6ePBuvr3FYwIEC0IHBZYVlayHbtYGb3PkadyQOS7d9hXm47OLDqGYQRaUlmRIBJCn16Xpv1phTfth5313658yk9vsQqkjSu+PP2bblTFtleJBpvMjtX4CwfjLWrckxm58q2PT9QF/stzH+Shn4ik1tVJ2GgQdykONTKKnyQRhJhoEuYl2J4R/96YoKkDPFO72f13jEG5L9Dsj3b4zvS1MxCfrYwuTSXaGO/KKsHj6TBMew0d3/rNJlL85zI5PLRJgBTz41FMfDunHgIZ0h80Dk03MdI2wyPnpwqCEls2bOX/p5CB2xSY08TQeLvt11a2LmzcGaZH9gYia5z7DitQMuhg0Wulu6+SfWStB6k1HnlCAezsce/g16c5eCtkRPX28R397Xz/K7dzHf1su3KAh3LC+yYv1Tg//GCqdjkpSDEubuNF9370orZ3WvBZHYO+lZMmD5vCTJmJb1Ll8stJqZJMoeG+AlV0fdrB3iBkLbsoXBJOidJMI3u7E509ieZlpPu49t2qai7ED8WheONJ4xg85+zmcWy2nCtMEGSgmKN2uMXZcVJGsi8AJp6LhvNAIcik8N8LhM5KJMWUAWE9t+kfIaRGr5DtAadvJg5Js7+gR6WYyaH0LGddP/QrBTPl/89TQTJSH8X9EeCYoRks8Xy6WivJL+gDFb7Morh81AqkqdUPsuVYSiXYUDagT0l0wGczUZh4Re6+5HFeS5393H10CAj/V08nBC6CxSYfSazc4y3rl3T8u34jqN7OBHTWvYfjPwwg/6Z5maKXaZgkhFnJCYQO6VQsIQmnfhErBihqS5eHt8G/H+voRfzCT54/DytE+PRy6ZKlCPOQG97gfZQiTk3DWG/GZ+6XJD/eERWXCOcmllgKGHXgXphgiQF8dnr1MxC3lRVEJ8fbn8QfPcNIrSz+sF9t4v7pmeg5ACVJoQxrj15fOeIO8zjs7a1hgnGz6s0CqzY9U5kciyfPsV8rovbDo5GK9mPZfLhlgO97exm5fkmzQzTCK+2c2e4pbedjEba4kj/3pLp0+A1w2L4MM9nXphk9/hx2hbmmO/oprOtlX2tyUEcxTQsL0yXT59ioLedF5fa8qZGKC7049dPS1J6H/mUJsz22t42tm9fba7yuzSEfSox6MPVT3ziUKz9F8P3x7WYPL2pG1Y/39BHN9jXAbnS7wXZP9BTto3GTX2e8Hl3vjgRLVR0e+HNup0bGoEJkjKE6q5v6N5ZDYVqbBrTgydp5hM6KC/IHNOBicJTbqac1CnKRaqEjbFeMefFhE258sznVhazeQdw6G9aSoim8hE8vm68qa8c3kxRLO/lBkk/IHx9YjxfZ/7Y1HNjZGcX83lqnRjn7AScd2/eWO7sZKG1G29cax8dhZjDOUk7DWnZE63WPgS0xJ5BnDCQopSgD59b++goLT1RpFglFFvQmHZn2tD/5AWH10bCz1BcKK6lXYcCId6Gij23pLYy0TPA4q72VXmeP1lcwwvvnzafU7lo+5pDl9xrEHqLr2urNSZIyhBXK+MaSDhDiDsFw1kLpB+Q4vcPV2FXQvx+lXamtCaitMfLXefChbmiefUr2b35rpLghWIzTn+P1lg9lctnueeZJhKo1aU79d2z9PV3MfOjtzOZnWN47EkuurfsFXv+pYSdX0iXpo2Fg2E4uy8n3MOBO5621Lm+75y/lOOyJvsGi2kU+wd6WD59Kvq9v7R/LA1xAVQqnSfuR0lr/lw+fYre3DTE4qfy62zKnA8UhDYXy6e3cuT9VSnN2LXABEkVFFtABIUVXu0sP1wBXw/S+i3iK4vXGx/CXGptSEhSOeO+oFqUMx7hFa5DiAclXNO9jat3dLF8TTdz11xHZ3Z8lb07nu+0bSGtPyzMd7FntBbi5/m+MS0LbN+ebBaK53kol2E+lyk56Iak9f95/ATFT9RuTmEKTevr82nGTpdNWpZGrQdZKyZIUhDaJ+P7HhVLFzJ/8iSwMgOJh41WoqGUo5YzkDSdt9zsvJamssG+Dg4cfikTPefhubGyaSs1DdSTf/fawwU+qXBleshZ6aZlqZMbSlwrKSTVE7altANeEqXqdf7kSZanLhcEOlRy/lpJ0g4qWWCYhNf2vfWglLZVC6q9ZimrRtxHW039V4oJkjUSN1n5bROgcMUtUDLapdZUMqOsZ4dZDxqlxpcj7kuIP+dwELi2t4191/WybbCPw8PDPHDmu/jg2GL1UhC0MTEeCaWEvaFKUSzwopLzK9HcQpPYM9PTRQe5+PH5hL3byt0n6XPSd1hxXpcSipXe1+MnkLUg3raT6i0M8gkjyhqBCZKUFFvUFVKssn1niGsm8XRho7hy+jTL2p5v4I2cXTQT+WflNnEEF2b8b29alXatzyet2aQSKslL95HDtI+OMn/yJIvupVxpAwQqNV+VExy11I7D61ViCkrLRusT1Uzcymn8JzK5fADPgNu9upGYIClDsUioet2rkuPrSbUDTrkyxTtdpT6nZnlmZ49FO8Duji2gzM+Anz6WP9Y+OsrMTa2r3k+zFvKvNGiiNuVNYrD2tQ1Ja62gvuVZq9YWOuYvynyxtxvUhILy1yAYoVJMkKSgWjtsKYppKVuVYp3Wb26Z5rWhcSFU60GmVnZugAuHrqd958789/Vyqq51sKyUkf4usotrG3bWu4+kWe+RRPvoKFf19oLWIVMlaKTp2gTJGihXMbWquLVGEJVrQPUeLIzKCVerp32PRLjOZiOxt7+TnTubQ2P0mmEjBtuknZE3CyZIUlDJgrQwnf/uZ1K1dL41C8U6hWlaK6HF5VZbhzvRplmdnWb9QqVO6kru0cg6rfe9ahHu3CxBK43SKpMwQVID4oOm3weqVCRNrTp8NayXnRzSDxDNqD1VWgYfVloKH7p528FrK/LJNaugruegtl4Dd7zey5Wt2kiwtbIeAsUESQOppNNftWdPge18oxLvfH6QPNCc419TsJ5hzPXwJ53I5DiRyaXatbeZqbbtNkt4ej0wQVIB6626NhP1fhbrqabXkgOHK1vfUSuaUVPxebp04cI656Q2JO1AHbKVxgsTJOtMM3b4ZqSSTrmVOnA9aCb/VrP4H0Lik5s0puz1wBYkbkFq2XkbFVVm1A6/rTxs7PppxoG/1uQFR0pOZHJcvDjL4Rt2rjoOm+NZmSCpIV4Y1MIWempqlmlNt7q5GfHrPobWMZhgo9CsddxMmkkzs5l9H2kxQVID4uG9lXQ8PyvZ6I7ItFQaydKMg+x6b6myEWjmjQ/XSvgCMTABEmKCpIE0WpVtBtU5LmT3N2h2a7Pp5qHcdi0bhVpuC7/ZMEHSADZ65FE1bDWNa6tTLpKpGpphYhSymRcaV4oJki1As3VAo7lppgGyGdts3KRlWq8JkpptDqqIAAAMP0lEQVRSzYKlzd4YbSFieTaDwPeD7MQ656OeVNpX9w/0cEFWb3uzkes5jgmSGlILdd7PAvf299Z8Y7tGDlSbZUHhViTfTjb55MaoHSZI1pnNNCsx1sZmFLbWrrcWLfW8uIi8QUSOi8iYiHwg4fe9InK/iDwlIg+IyHDw2++KyNPu747g+I+IyLdF5AkR+bqIHKhnGWrJWt9nsFHuV4qWPXsbuold++jopjcXNjPzJ08WRM5ZXWxu6qaRiEgr8Eng9cAZ4FERuVdVvxMk+wjwZ6p6l4i8Fvgd4G0i8uPAy4BbiF6n9qCI/IOqTgN/DLxJVZ8VkXcDvwG8s17lSEPa7eW3Imbi2lrU01G/lftRs1NPjeQoMKaqJ1R1AfhL4E2xNDcC97vPXwt+vxF4UFWvqOoM8CTwBvebAtvd5z7gbJ3y31D8Dqn1oJk0k81APevKMDYi9fSRDAHjwfczwL+JpXkSeDPwceAngV4RucYd/18i8lGgC3gN4DWZ/wLcJyKzwDTw8qSbi8gvAr8IMDw8zIU67jh68eIsEL35DeDChdIvJip1ja7cRQBmlpdrlLvV90iKIKk32yX6X+rZZLPZBuWmOmr5HLPZLD9Ymi04th71E1KL8l25GLXjq3p7V/22Ueq5lmz2MtdTkEjCsfhbi98HfEJE3gk8RBQ1eEVVvyIitwLfADLAN4Er7pz/DrxRVR8RkfcDHyUSLoU3Uv0U8CmAI0eO6M46vttjWqPZaTVRVv4agzeMAHDhwgVqneeN8HqTetZTrahFfYf8UJOVuRbZmb90CaDoO3U2Qj3Xms1c5noKkjPASPB9mJgZSlXPAj8FICI9wJtVNet++xDwIffbF4DnRWQAuFlVH3GXuBv4xzqWwTCMNWDO9a1FPX0kjwLXi8ioiLQBbwHuDROIyE4R8Xn4deCz7nirM3EhIjcBNwFfAS4CfSJygzvn9cCzdSxDwzH7u9EsWFs00lI3jURVr4jIncCXgVbgs6r6jIj8FvCYqt4LvBr4HRFRItPWr7jTtwH/T0Qg8oP8vKpeARCR/wp8UUSWiQTLf65XGdJijmzDMLYydV2QqKr3AffFjv1m8Pke4J6E8+aIIreSrvkl4Eu1zen6Y2GyGwebOBhGIXVdkGgYhmFsfkyQGIZhGFVhgsQwDMOoCtu0scnw9ve1LGo0jFpiviAjLaaRGIZhGFVhgsQwDMOoChMkhmEYRlWYIDEMwzCqwgSJYRiGURUmSAzDMIyqMEFiGIZhVIUJEsMwDKMqTJAYhmEYVSGq8ZcWbj5EJAOcWu98VMhOoH7vB25OrMxbAyvzxmGvqg6US7QlBMlGREQeU9Uj652PRmJl3hpYmTcfZtoyDMMwqsIEiWEYhlEVJkial0+tdwbWASvz1sDKvMkwH4lhGIZRFaaRGIZhGFVhgqRJEJF+EfknEXne/d9RIu12EZkQkU80Mo+1Jk2ZReQWEfmmiDwjIk+JyB3rkddqEZE3iMhxERkTkQ8k/N4uIne73x8RkX2Nz2VtSVHmXxWR77h6vV9E9q5HPmtJuTIH6X5aRFRENkUklwmS5uEDwP2qej1wv/tejA8CDzYkV/UlTZkvA29X1ZcCbwD+UESubmAeq0ZEWoFPAj8G3Aj8rIjcGEv2C8BFVT0AfAz43cbmsrakLPO/AEdU9SbgHuD3GpvL2pKyzIhIL/Ae4JHG5rB+mCBpHt4E3OU+3wX8RFIiETkMXAd8pUH5qidly6yq31XV593ns8B5oOwCqSbjKDCmqidUdQH4S6Kyh4TP4h7gR0REGpjHWlO2zKr6NVW97L4+DAw3OI+1Jk09QzQR/D1g07xP2wRJ83Cdqp4DcP+vjScQkRbgD4D3Nzhv9aJsmUNE5CjQBnyvAXmrJUPAePD9jDuWmEZVrwBZ4JqG5K4+pClzyC8A/1DXHNWfsmUWkX8NjKjq3zUyY/XmqvXOwFZCRL4KDCb89D9SXuLdwH2qOr5RJqs1KLO/zi7g88A7VHW5FnlrIEmVFQ+XTJNmI5G6PCLy88AR4La65qj+lCyzmwh+DHhnozLUKEyQNBBVfV2x30TkRRHZparn3KB5PiHZK4BXici7gR6gTURyqlrKn7Ku1KDMiMh24O+B31DVh+uU1XpyBhgJvg8DZ4ukOSMiVwF9wFRjslcX0pQZEXkd0aTiNlWdb1De6kW5MvcCPww84CaCg8C9InK7qj7WsFzWATNtNQ/3Au9wn98B/E08gaq+VVX3qOo+4H3AnzWzEElB2TKLSBvwJaKy/nUD81ZLHgWuF5FRV563EJU9JHwWPw38X93Yi7zKltmZef4EuF1VEycRG4ySZVbVrKruVNV9rg8/TFT2DS1EwARJM/Fh4PUi8jzwevcdETkiIp9e15zVjzRl/k/AvwfeKSJPuL9b1ie7a8P5PO4Evgw8C/yVqj4jIr8lIre7ZJ8BrhGRMeBXKR211/SkLPPvE2nWf+3qNS5cNxQpy7wpsZXthmEYRlWYRmIYhmFUhQkSwzAMoypMkBiGYRhVYYLEMAzDqAoTJIZhGEZVmCAxmhoRuSYI+510ux77799ocF5K7labtCuziDzgdoP1eb42dk7BLrAisk1E7hKRYyLyrIj8ujveISLfEpEn3U7I/ye4xmfc8adE5B4R6SmS/58Qkd8s8ltu7U8m8XpflRI7WBubCwv/NTYMIvK/gZyqfmSd7v8a4BFVvSwivwy8WlXvCH7/ONGGklOqeqc79gDwvqRFZ24X2L8n2j/sTlV9TER+jmiR2ltEpAv4DvBq4BTQrao5EdkGfB14r6o+LCLbVXXaXfOjwHlV/XDC/b7hrn0h4becqiYKoLUgIu8AhlX1Q7W6ptG8mEZibFj8LFpEXi0iD4rIX4nId0XkwyLyVjeDPyYiP+TSDYjIF0XkUff3ykruV2q3WlnbrsxJu8Aq0O22SekEFoBpjfBawzb3py5fXoiIO2fV7FBEbgDmvRBxq6+/6Z7DB4N0PU7b+rZ7dm9yxz8oIu8N0n1IRN4jIrtE5CGnbT0tIq9ySe4FfraCZ2FsYEyQGJuFm4H3AoeAtwE3qOpR4NPAf3NpPg58TFVvBd7sflsr+d1qpfyuzJ9zA+3/dIN9qV1g7wFmgHPAaeAjqjrlzmkVkSeI9iT7J1XNv89CRD4HTAIvAf4oIQ+vBL4dfP848MfuWUwGx+eAn1TVlwGvAf7A5fkzuC1cXHnfAvw58HPAl1X1FqI6eAJAVS8C7SKykXcwNlJigsTYLDyqqufcxn/fY0UzOAbsc59fB3zCDcb3AtudeakiZGW32t93h/K7Mickf6uqHgJe5f7eJiu7wP5aQvqjwBKwGxgFfk1E9gOo6pIbsIeBoyLyw/4kVX2XO+dZIOktkruATPD9lcBfuM+fD4sH/LaIPAV8lWgb9OtU9QXg+04A/ijwL6r6faL9pd7lzI6HVPVScK3zLk/GJscEibFZCHeOXQ6+L7Oyy3UL8ApVvcX9DcUGPkTEaw/3Jd1EVnarvT3YrfYVwJ0i8gLwEeDtIvJhAFWdcP8vAV8gEhThLrAvAC8n2gX2CNEM/x9VddFtZPjPREIrj6r+AHiA6I2R4fEl4G4ibSvOLNARO5bkIH0rkZ/nsBNaLwbnfZpoC/R3AZ9193yIaC+0CeDzIvL24Fod7r7GJscEibGV+ArRpnpA9D74eAJVfZcTMm+M/yZFdqsttiuziFwlIjvduduA/wA8XWYX2NPAayWim0jIPOf8O1e7a3USaVfPuXQH3HEB/iPwXELZnwUOBN//mcg8BZHw8PQROesXXXBBGJn2JSLhdSvRxoRIFLl2XlX/lMj89bIgL4PACwl5MTYZ9j4SYyvxHuCTzmxzFfAQ8EsVnB/uVgtwWlVL7eraDnzZCZFWIlPRn5a5xyeBzwFPE5mZPqeqT4nITcBdEr0XvIVoZ9m/c2ayuyR6Z4sATwK/nHDdh3D+Drc9/XuBLzgH+heDdH8O/K2IPEbk78gLJVVdEJGvAT9w2g9EEWXvF5FFIAd4jeQw8LDbEdfY5Fj4r2FsEVx48t+q6lfXeH4LkcP+Z1T1+RT3uldV71/LvYyNhZm2DGPr8NtA11pOFJEbgTHg/nJCxPG0CZGtg2kkhmEYRlWYRmIYhmFUhQkSwzAMoypMkBiGYRhVYYLEMAzDqAoTJIZhGEZVmCAxDMMwquL/A2Zyzqc41X+KAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = lc.fold(period=3.0712).plot(color='C0', alpha=0.2, label='With Motion')\n", "ax = corr_lc.fold(period=3.0712).plot(ax=ax, color='C3', alpha=0.2, label='Motion Corrected');" ] } ], "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 }