{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# How to cut out Target Pixel Files from Kepler Superstamps or TESS FFIs?\n", "\n", "You can use `lightkurve` to cut Target Pixel Files (TPFs) out of a series of standard astronomical images, such as [K2 Superstamp Mosaics](https://archive.stsci.edu/prepds/k2superstamp/) or TESS Full-Frame-Images (FFIs).\n", "This brief tutorial will demonstrate how!" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's assume you have downloaded a set [simulated TESS FFI images](http://archive.stsci.edu/tess/ete-6.html) to a local directory called `data`. `lightkurve` will assume that the files are given in *time order*. So we'll sort the filenames first:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from glob import glob\n", "import numpy as np\n", "fnames = np.sort(glob('data/*.fits'))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now we use the `KeplerTargetPixelFile` class and its function `from_fits_images()` to create the new TPF. This will cut out around the position keyword. You can pass a pixel position in units of the original image or RA and Dec coordinates." ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "scrolled": false }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "100%|██████████| 3/3 [00:00<00:00, 19.98it/s]\n" ] } ], "source": [ "from lightkurve import KeplerTargetPixelFile\n", "from astropy.coordinates import SkyCoord\n", "tpf = KeplerTargetPixelFile.from_fits_images(images=fnames, \n", " position=SkyCoord(257.13700, 24.48958, unit='deg'), \n", " size=(9,9),\n", " target_id='MyCutOut')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We then have a fully functioning `KeplerTargetPixelFile` object! You can read more about such objects in the [tutorial on their use](https://docs.lightkurve.org/tutorials/1.02-target-pixel-files.html)." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGRCAYAAAA5GvouAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xtc1GX6+P/XcBIEBIFQUxFSSUPUBI0y8USah8pIVzFb67eFu+uWa1rbx07Wlu2nWreDSW627noAzUrTSiBFRfHAkiiSJhQCmil4muF8mJnfH3x5f0QOzrxnRga4nj3m8Yh5X3PfN4rMNff7vq9bYzQajQghhBBCWIFDaw9ACCGEEO2HJBZCCCGEsBpJLIQQQghhNZJYCCGEEMJqJLEQQgghhNVIYiGEEEIIq5HEQgghhBBWI4mFEEIIIaxGEgshhBBCWI0kFkIIIYSwGqfWHoAwzQMPPKD6tQsWLCAqKsqKo7FfP/74I0eOHKFr165MmjRJdTtfffUVq1evBiAhIQEPD48G12fNmkVZWVmD55ycnHBzc8PDw4OAgAD69u3LyJEjCQgIUD0OU+Tl5bFgwQLl61GjRvH888+3+JqkpCRWrFihfP3kk0/y0EMP2WR8JSUl7Ny5k6NHj1JYWIhOpwOgS5cuBAUFMWzYMEaNGoWXl5dV+92xYwdXrlxh2LBhDBgwwKTXFBUVkZyczNGjR7lw4QKlpaW4u7vTrVs3hg4dysSJE/H397fqOK/1/fffc+rUKXr16kVkZKTN+hHCliSxaCO8vb2bfL6yspLKysoWYzp16mSzcdmbU6dOkZCQQFBQkEWJhak6deqEm5ub8nVZWRklJSX8+uuvHD58mPj4eEJDQ5k/fz49e/a0+XgADh06RGlpaaNk6FrffffdTRnLtm3bWL9+PRUVFcpznTp1wtHRkYsXL3Lx4kX++9//8p///IfHHnuMBx980Gp979ixg9OnT+Pu7n7DxMJoNLJx40Y2b95MTU0NAA4ODnTu3JmSkhK0Wi05OTls2bKF6dOnExMTg0ajsdpY6x05coRt27Zx1113SWIh2ixJLNqIdevWNfl8fHw8CQkJLcYI25k4cSJPPfVUg+dKS0vJzc1lz549pKamcvz4cZ555hmWLl1KaGioTcfj7+9PUVERe/fuZcqUKU3GnDlzhlOnTtGlSxeqq6uVxNTa4uLi+PbbbwG47bbbeOSRR7jzzjvx9PQE6pLi48ePs2vXLg4cOMDBgwetmliY47333iMlJQWAQYMGMWvWLEJCQnByckKv1/PDDz+wadMmsrKySEhI4MKFCyxcuLBVxiqEvZM1FkJYmYeHB3feeScLFy7kb3/7G15eXlRXV/Pmm29y6dIlm/Y9fvx4AHbu3NlsTP21MWPG4OjoaJNx7NixQ0kq7rvvPpYvX05kZKSSVAC4uroyfPhwXnjhBd5991169Ohhk7HcyLZt25SkYtKkSSxbtowhQ4bg5FT3ucvR0ZHBgwfz5ptvMnXqVABSUlLYvn17q4xXCHsnMxYdiF6v58SJE6Snp/PDDz9w6dIldDodbm5uBAUFMXbsWMaOHdvkm8219/ETEhK4fPkyX3zxBVlZWVy9epXevXvzwQcfKPEGg4HExES+++47zp49i6OjIwEBAdx///2MGzeON954g8OHD/Pggw82+sRfr7CwkO3bt5OVlaW8IXfr1o3w8HCmTZtG165dldjS0lJiYmKUr0+fPt1oXYot1xE05/bbb2fx4sW8/PLLlJWV8fnnnzNv3jyb9RceHs63337LTz/9RH5+PoGBgQ2u6/V6du/eDUBUVBS7du1q1MbFixf53e9+h8Fg4LXXXmPYsGHN9rdq1Sq+/vprgoKClL//iooKZfYsODiY+fPn3zCBCQ4Opl+/fg2ee+aZZzh9+nSLf2+ffPKJcuvgpZdeAhqujwFYvXp1g6/h/9bNlJeXEx8fD0C/fv2YN29ei7c4nnzySXJycsjJySE+Pp6oqKgGt8KaGs/1Dh48yLJly3B3d2fjxo1A43Uyhw8fbvTzu2TJEu6+++5mxyaEvZAZiw6koKCAJUuWsHXrVnJzc9HpdDg7O1NSUkJWVhbvv/8+r7/+OrW1tS22k5mZycKFC0lJSaGsrKzRm0ZNTQ1//etfiYuL46effqK6uhqNRsOpU6f4xz/+wcqVK2841i1btvD000+TmJjIuXPnMBqNGAwGCgsL+fLLL5k/fz4nTpxQ4h0cHPD29lZ+yTs6OuLt7d3g4erqquJPzXJDhw5l4MCBAOzevRuDwdDgel5eHg888AAPPPAAX331lUV9OTk5MWbMGKDpWYvvv/+eK1eu0LdvX4KCgppsw8/Pj+HDhwN1izybU11drSQpEydOVJ5PTU2lpKQEgJkzZ5o8K+LgYJ1fR66urnh7eyv9urm5NfpZqO9r7969yiJcU8bq6OjIb37zG6Aumd27d69Vxlz/81q/HsrZ2bnRmF1cXKzSlxC2JjMWHYiTkxMjR45kzJgx3H777Xh7e6PRaCgrKyM1NZX169dz5MgRNm/e3ODT//U+/PBDBgwYwJNPPqm8Of3yyy/K9Q0bNpCRkQHU7Z546KGH8PDwoLS0lC+//JLNmze3uLBw165d/Otf/8LFxYVHHnmEiRMn4uvrqyQW//nPf8jIyGDZsmV89NFHeHl50blzZ9atW6d8Wg0ICGgwg9Lahg8fzsmTJykrK+P06dP07dvXZn3dd999fPXVV+zZs4fHH39cmdKH/1u0ed9997XYxuTJkzl8+DDp6elcvXq1yYXB+/fvp6ysDFdXV8aOHas8f+zYMQBcXFxanO2wlYkTJzJx4kRlxuPRRx9tdsYjKysLqBtreHi4Se2HhYXh4uJCdXU1x48f5/7777d4zH369GHdunXKjMewYcOanfEQwt7JjEUHEhAQwAsvvEBERARdu3ZVpnzd3d2ZNGkSzz33HADffPMNRqOx2Xb8/PxYunRpg0+89TsedDod27ZtA2DatGk8+uijShLh4eHBb3/7Wx588EFKS0ubbLuiokKZtl60aBGzZ8/G19cXqPtEGxgYyMsvv0xoaCharZZvvvnGkj+Sm+baP6vz58/btK8+ffrQr18/tFot//3vf5XntVotGRkZODs7M3r06BbbuPPOO+nevTu1tbXNrteon80YNWoUnTt3Vp4vLCxUxnFtUmOP6scaEBBg8lidnJyULcQFBQU2G5sQbZUkFkIxZMgQnJ2d0Wq1nD17ttm4hx56CGdn5yavpaenU1NTg6OjI9OnT28yZsaMGc1Oex84cIDS0lJuvfVW7rnnniZjHBwclOn+I0eOtPAd2Y9rZ2jqbxPUu+2229i+fTvbt2+32hqQ+hmJa5OC3bt3U1tbS0RERIszRgAajUb5JJ6cnNwo0Txz5oxyK+r6T+z139+1CzXtVf1Yu3TpYtbr6uOv/7sUQsitkA6nsrKS5ORkDh8+TGFhIaWlpU2uqbh06RK9e/duso077rij2fZ//vlnAHr37t1swSNvb28CAgLIz89vdK3+zaq4uJjHHnus2X7qaw0UFxc3G9ORRUZG8umnnyprKrp27aokGaYWS4uKimL9+vX8+uuvZGVlMWTIEOVacnIyUDcTExwc3OB1Lc12tRcd4XsUQi1JLDqQ4uJiXnzxRX799VflORcXF7y8vJTbIlqtFqPR2GJtg5YqJGq1WgB8fHxaHIuPj0+TicXly5eBusTh6tWrLbYBUFVVdcMYe3DtrZ+b8Unew8ODiIgIUlNT2b17N4MGDaKgoAA/Pz+GDh1qUhteXl7cc889pKamkpSUpCQWNTU1yvbMptYXdOnShStXrrSJT/Oenp5cuXJFqQZqqrY0KyPEzSaJRQeycuVKfv31V3x8fHjiiSe48847GyUJv/nNb6ioqGjxE1lLq/frX6e2KmH9jonIyEhlzUd7cPr0aeX/b1a9hqioKFJTU9m5c6eyrmPcuHFm7b6YPHkyqampHDx4EK1Wi5eXFwcPHkSn0+Hq6qrckrpWQEAABQUFFBQUUFtba9frLHr37k1hYSGFhYUmj7W2trbB2gwhREOyxqKDKC0tVdYjPPPMM4wZM6ZRUlFZWdmg9LIa9bsH6mcemtPc9frXNzWb0ZbVL6L08PBoVFvCVoYMGcItt9zCmTNnlN0g9QW0TBUSEkJAQAC1tbXKLEVzizav7RfqtqNasgamfutndXV1szHXn9dirmvHeu1C15ZkZGQoYxo8eHCDa/VJW/2tuqaUl5erGaoQbYYkFh3E5cuXldmA2267rcmYo0ePWtxP/TbKwsJC5bbI9a5evcqZM2eavFa/fqOwsFDVivv6mRJ7ugd+9OhRTp48CZg/Y2AJBwcHxo0bB9R9yg4JCeHWW281u536M1eSkpI4d+4cx48fB5q+DQI0qLC5adMm9Hq9Sf1cX9+jfoHpxYsXm31NTk5Os9dM+VkYPXo07u7uAHz22Wc3HKter+ezzz4D6nZTXX+eR/2YW1r7c+rUqRuOWYi2TBKLDqL+lyc0nJavV11drVQgtMTw4cNxdnZGr9fz5ZdfNhnzxRdfNPsL/N5771XGGhcX1+InP6DRttX6T9CWfpK1lpycHN59912g7k3nkUceuan9T5gwgWnTpilbf9UYO3Ysrq6u/PLLL3zwwQcYjcYmF23Wc3NzY86cOUDd9//RRx/d8A37p59+anDaKqDM7KSnpze5wPjgwYPNJqhg2s9C586dmTVrljKGVatWtZiIrF69mtzcXABmz57daMamflvx2bNnmxxbUVERe/bsueGYm9uOLURbIIlFB+Hr60ufPn0A+Pjjjzl58qTyC/Snn37ipZde4vz58xbfD/fy8lJKEW/ZsoWEhATlF3tpaSnr169n69atzW53dHd3V0p8//DDD/zP//wPWVlZyhuT0Wjk3LlzfP311zzzzDONfknXf4/FxcWtthW1rKyMzMxM3nvvPf7yl7+g1Wrp1KkTL730UpOLWq1ZefN6/v7+/O53v+N3v/ud6gPQ3N3dGTVqFFD3dwLNz1bUmzx5slKN87vvvuPZZ59l3759DRZ0VlZWcuTIEd555x0WLVrUYFExoMwGXLx4kX/84x9cuXJFed2OHTv4+9//3uLiyfqfhcOHD7e4EHjatGlKXY8dO3bw4osvcuzYMeVnTq/Xc/z4cV566SW+/vproG6mo6kD04YOHUqXLl0wGo2888475OfnYzQa0ev1fP/99yxZsqTFf2P1Y87NzSUvL6/ZOCHsmf2uqhJWN2/ePF599VV+/fVXnn/+eVxcXHB0dKSiogJnZ2eee+453n///RuW9L6ROXPmkJ+fz5EjR4iPj2fjxo24u7tTVlaGwWBg6tSpFBcXc/jw4SbrYYwfP57q6mr++c9/curUKV588UWcnJxwc3OjoqKiwfiunzru378//fr146effuLVV1/F3d1dmQGZNWvWDStOmispKYnU1FTl6/Ly8kZrAgYPHsz8+fNV3YawF5MnT1bWaTS3aPN6f/rTn+jVqxfx8fHk5eXx9ttvA3UzGhqNpsFag86dOyvJS73+/fszdepUvv76a1JTU0lNTcXd3Z2KigoMBgPjxo3D3d292cPA7rvvPhITEzl9+jRz587Fy8tL+Xn78MMPG8w2LFq0iO7du/PFF19w/Phxjh8/joODQ4OfW6grjvXII480O/vj4uLCH//4R95++21Onz7N008/jaurK3q9npqaGgIDA5k1axbvv/9+k68PCwvDz8+PixcvsmDBAjw9PZUy9fPnz2+VSqZCmEsSiw4kNDSUd955h40bN/LDDz9QUVFBly5duOuuu4iOjiYoKKjZX3jmcHZ25pVXXmlwCJnBYOD2229n8uTJjBkzhhdeeAFoeIvmWpMmTSIsLIxvv/2WzMxMLly4QFlZGW5ubgQGBhISEkJERESTNTWWLl1KfHw8mZmZXLx4UZkxscWiuaqqKmXLa33y4+vrS0BAAH379uXee+9tth5IW9KvXz+6devGhQsXml202ZRp06Yxbtw4du7cydGjRyksLFRmLfz8/AgKCiI8PJzIyMgmZ7FiY2MJDAwkMTGRwsJCDAYD/fv3Z9KkSYwfP55PPvmk2b779u3L66+/zpdffklubi5arVZJEK5fz6HRaJgzZw4TJkwgOTmZzMxMzp8/T1lZGR4eHnTr1o2hQ4cyYcIEunfv3uL3PHLkSN544w02b95MTk4ONTU1dOvWjcjISKKjo1ucSXN1deWtt94iISGB48ePN9i221a2VguhMdrTKjfRIdTW1jJnzhzKysp4+eWXGTFiRGsPSdzAuXPn+P3vf4/RaOTvf/97s+srhBBC1liIm+7bb7+lrKyMTp06ERIS0trDESaoPz+mb9++klQIIVokiYWwiTfffJPDhw83WKx36dIlNmzYwL/+9S+g7nZHc7dChP3Izs5mx44dADz88MOtPBohhL2TNRbCJg4dOsShQ4eAuvvGDg4ODdY4DBs2rMWzQETrKi8v5+mnn6a6ulrZURESEtJogaUQQlxPEgthE/Pnz+fo0aOcPn0arVZLVVUV3t7e9O3blzFjxhAZGXnTCkUJ8xkMBoqKitBoNPj4+DBixAjmzp0rf2dCiBuSxZtCCCGEsBr5+CGEEEIIq5HEQgghhBBWY7drLMrLy9m0aRN5eXnk5eWh0+mIiYlh9uzZquKgrnTyunXryM/PR6fT4eLiQs+ePZkyZQpjx45tdixJSUmsWLECV1dXNm/ebPXvVQghhGgv7HbGoqSkhKSkJGpqaoiIiLA4DurOcPDz8+O3v/0tr776KgsXLqRbt24sX76cTZs2NfmaS5cusWbNmibPeBBCCCFEQ3Y7Y+Hv709CQgIajQatVktycrJFcVBX0vr6g5hGjBjBhQsXSExMZObMmY1e89FHHxESEoKHhwcHDhyw7JsSQggh2jm7nbHQaDSNDpiyJK4lnp6eODo6Nnp+9+7dZGdn84c//MGi9oUQQoiOwm4TC1syGAzo9Xq0Wi3ffPMNmZmZPPLIIw1irl69yieffMLcuXPx8/NrpZEKIYQQbYvd3gqxpbi4OBITE4G6EyljY2OZNGlSo5hevXoxefJkk9s1GAxcuHABJycnk2ZRnJ2dmzw2XAghbhaj0UhFRQU+Pj6tWgCturqa2tpaVa91cnLCxcXFyiMSanXIxGLGjBlMmDABrVZLeno6q1atorKykujoaADS0tJIT0/n/fffN+s2y+XLl4mNjbXVsIUQwmbWrFnTarOz1dXVPPn/PcQVrbq3pK5du7J69WpJLuxEh0ws/P398ff3ByA8PByAtWvXMn78eFxcXPj444+ZOnUqPj4+lJaWAiiZdGlpKU5OTri6ujZq183NDYDKPS6gNyEhMQBGy9aHCBNo7OCOn9HQ2iOoI4V2xfUcjbhGVSu/v1pDbW0tV7ROfPK/uXR2M+/fSnmFA0/9pT+1tbWSWNiJDplYXC84OJgdO3Zw/vx5vL29uXr1Klu3bmXr1q2NYmNiYrjrrrt46aWXGl1TZjeqHaBWEga7YeHiXquwlwRS8grRDEsXwVtDJ9daOrmal1jojXbwwUE0IIkFkJWVhYODA927d8fNzY1ly5Y1ivn888/Jzs5m6dKldOnSpRVGKYQQQtg/u04sMjIyqKqqoqKiAoDCwkLS0tIACAsLU25HmBq3YsUK3NzcCA4OxtvbG51OR1paGvv27SM6OhovLy+ARrUuAHbu3ImDg0OT14QQQljOgBGDmdNq5sYL27PrxCIuLo6ioiLl67S0NCVhWL16tZIwmBo3YMAAdu7cSUpKCmVlZbi6uhIUFMSzzz7bYklvIYQQtmf4f/+Z9xphb+w6sfj000+tGhcVFUVUVJSqsSxcuJCFCxeaFOtydyXQ8H6lvtAJfYFd/3ELIUSr0huN6M1cYGxuvLA9eaezgeqDrrJ4UwghzGRUcSvEKLdC7I4kFkIIIeyCHiN6MxMFc+OF7ck+HSGEEEJYjcxYCCGEsAuyK6R9kMRCCCGEXZDFm+2DJBZCCCHsggHzt4/KdlP7I4mFDch2UyGEMJ8s3mwf5J3OBmS7qRBCmE9vrHuY+xphXySxEEIIYReMmH9rQ/IK+yPbTYUQQghhNTJjIYQQwi7o0aDHvNvI5sYL25PEQgghhF0wGOse5r5G2BdJLIQQQtgFmbFoHySxsAHZbiqEEOaTxKJ9kHc6G5DtpkIIYT6DUYPBaN7vTnPjhe1JYiGEEMIuyIxF+yDbTYUQQghhNTJjIYQQwi4YcECv4jXCvkhiIYQQwi7IGov2QRILG5BdIUIIYT5ZY9E+yDudDciuECGEMJ/e6KDiEDK5FWJv7DaxKC8vZ9OmTeTl5ZGXl4dOpyMmJobZs2erigPIy8tj3bp15Ofno9PpcHFxoWfPnkyZMoWxY8eaHSeEEMJ6DDiYfQiZrLGwP3abWJSUlJCUlERgYCAREREkJydbFAdQVlaGn58fkZGR+Pr6UllZyd69e1m+fDlFRUXMnDnTrDghhBBCNGS3iYW/vz8JCQloNBq0Wm2zCYOpcQChoaGEhoY2eG7EiBFcuHCBxMREJWEwNU4IIYT1yBqL9sFu55A0Gg0azY1/YEyNa4mnpyeOjo5WixNCCGE+vVHz/9ZZmPOQxMLe2O2MhS0ZDAaMRiOlpaXs37+fzMxM5s2bpzquEUcTVx8ZAIP8oxBCCAAjGrPXWBhlxsLudMjEIi4ujsTERACcnJyIjY1l0qRJquOu5zq+yqRx1OY6UZvrbMbIhRCi/dLjgB7ztoXIrRD70yETixkzZjBhwgS0Wi3p6emsWrWKyspKoqOjVcVdr3JXJ9Cb8MNubmpuKo193OHSONjJP3h7+POwk+lao97cuoa2GISZ+wlFh6HHAb2ZPx+SWNifDplY+Pv74+/vD0B4eDgAa9euZfz48Xh5eZkd14heI3UshBDCTHXbTc1LLAySWNgdO/go1/qCg4PR6/WcP3/eKnFCCCFER9UhZyyul5WVhYODA927d7dKnBBCCPPV7Qox/zXCvth1YpGRkUFVVRUVFRUAFBYWkpaWBkBYWBiurq5mxa1YsQI3NzeCg4Px9vZGp9ORlpbGvn37iI6OVm5vmBonhBDCemTxZvtg14lFXFwcRUVFytdpaWlKwrB69WolYTA1bsCAAezcuZOUlBTKyspwdXUlKCiIZ599tkGpblPjhBBCWI/B6IDBzMWbcrqp/bHrxOLTTz+1alxUVBRRUVFWi2uOnG4qhBDmM6iYsZDFm/ZH3ulsQE43FUII8+mNGhxkjUWbJ7tChBBCCGE1MmMhhBDCLkgdi/ZBEgshhBB2oe5WiJmnm0ohV7sjiYUQQgi7YFBxCJmtTkYQ6kliIYQQwi7ojQ4qFm/aZixCPUksbEC2mwohhPn0OJi9o8DcY/WSkpJYsWIFrq6ubN68ucG1iooK1q9fz/79+ykpKaFXr15Mnz6dyMjIRu2YGmuLNu2dvNPZgGw3FUII8xmNGgxmzkCYU0/r0qVLrFmzBh8fH8rLyxtdX7ZsGbm5ucydO5eePXuyd+9e3nnnHQwGA2PGjFEVa4s27Z0kFkIIITqEjz76iJCQEDw8PDhw4ECDaxkZGRw9epTFixczevRoAAYPHkxRURFr1qxh1KhRODo6mhVrizbbAqljIYQQwi7UnRVi/sMUu3fvJjs7mz/84Q9NXj948CBubm7ce++9DZ6Piori8uXL5OTkmB1rizbbAkkshBBC2IW6s0LMf9zI1atX+eSTT5g7dy5+fn5NxhQUFNCrV69GswKBgYHKdXNjbdFmWyCJhRBCCLugR6PqcSNxcXH06tWLyZMnNxtTUlKCp6dno+frnyspKTE71hZttgWyxsIGZFeIEEKYr24GwtzXtHw9LS2N9PR03n//fTQaWVR/M8g7nQ3IrhAhhDCfHo3ZBbpb2m5aUVHBxx9/zNSpU/Hx8aG0tBSA2tpaAEpLS3FycsLV1RVPT88mZwXqn7t2NsHUWFu02RZIYiGEEKJd0ul0XL16la1bt7J169ZG12NiYrjrrrt46aWXCAwMJDU1Fb1e32CdQ/3ahj59+ijPmRprizbbAlWJxenTp3FwcGhT36gQQgj7Zu1bIV27dmXZsmWNnv/888/Jzs5m6dKldOnSBYCIiAiSkpI4cOAAo0aNUmJ37dqFj48PwcHBynOmxtqizbZAVWKxYMECQkNDefPNN609HiGEEB2UwahBb+YhZIYWKmS5uLgQGhra6PmdO3fi4ODQ4Fp4eDhDhw5l5cqVlJeX06NHD1JTUzly5AiLFi1qMItgaqwt2mwLVCUWHh4edO3a1dpjEUII0YEZ0Ji9yqLuEDLrHBiyZMkS1q1bx4YNG5SS2s8991yTJbVNjbVFm/ZOYzSaUxC1zmuvvcbFixf58MMPbTGmNqu8vJyZM2dSmdzKizc19rGLWONgJwtY7eHPw2gfZzAa9eaerGCLQcipUXbFyYjrxCo2bdpE586dW2UI9b87Q19xxNHVzGPTK40cf13fquMXDan6jRsTE8PZs2fZsmWLtcfTLrjcXYnLqIYPxz61rT0sIYSwawajRtVD2BdVt0LOnj3L2LFj+fe//83u3bsZPnw4t9xyCy4uLk3Gjxs3zqJBtjWy3VQIIcxXV57bzBkLK90GEdajKrF477330Gg0GI1G8vPzyc/Pb7LwiNFoRKPRqEosysvL2bRpE3l5eeTl5aHT6YiJiWH27Nmq4gDy8vJYt24d+fn56HQ6XFxc6NmzJ1OmTGHs2LFK3LFjx9izZw8nT57k4sWLuLu7079/f2bNmkW/fv3M/l6EEEKIjkJVYjFr1iybVzArKSkhKSmJwMBAIiIiSE5OtigOoKysDD8/PyIjI/H19aWyspK9e/eyfPlyioqKmDlzJgA7duygpKSEBx98kN69e6PT6diyZQuLFy/mtddeY8iQITb5noUQoiMzGDVozN4VYqPBCNVUJRZNzQZYm7+/PwkJCWg0GrRabbMJg6lxAKGhoY22Ho0YMYILFy6QmJioJBa///3v8fb2bhA3bNgwYmNj2bx5syQWQghhAwYcVOwKkczC3tht5U1TZ0SsMXPi6enJ1atXla+vTyoA3NzcCAgI4OJIAXamAAAgAElEQVTFixb3J4QQojGDUQMyY9HmWZxY5OXlkZubi06nIyAggLvuuguAmpoaampq7HL7j8FgwGg0Ulpayv79+8nMzGTevHktvqasrIyff/6ZwYMH37gDRxN/0g2AQRZ5CiEESGLRXqhOLM6cOcP7779Pbm6u8ty4ceOUxOK7775j1apVvPLKK4SFhVk+UiuKi4sjMTERACcnJ2JjY5k0aVKLr/n444+prKzkN7/5zQ3bdx1fZdI4anOdqM11NinWHPZSP0LjZB8TYho7OLxH42QfVfMMV67eOMjWY6gy7d+HzUk9DbtjMDqoSCzk79HeqPrNX1RUxAsvvEBJSQkREREMHDiQNWvWNIiJjIxk9erVHDhwwO4SixkzZjBhwgS0Wi3p6emsWrWKyspKoqOjm4xfv349e/bsYd68eSbtCqnc1Qn0JvzjsI+aSUIIIYTVqEosEhISKC0t5c9//rOylfT6xMLDw4PevXvz448/Wj5KK/P398ff3x+oq88OsHbtWsaPH4+Xl1eD2ISEBDZt2sRjjz3G1KlTTetAr5E6FkIIYSY9GoyqSnoLe6Kq8uaRI0e47bbbblifwt/fn8uXL6sa2M0UHByMXq/n/PnzDZ5PSEggPj6e2bNnm3QLRAghhHpSebN9UDVjUVJSwh133HHDOI1GQ3V1tZoubqqsrCwcHBzo3r278tzGjRuJj49n5syZxMTEtOLohBCiYzCgYo2FbDe1O6oSiy5dunDhwoUbxp05cwZfX181XQCQkZFBVVUVFRUVABQWFpKWlgZAWFgYrq6uZsWtWLECNzc3goOD8fb2RqfTkZaWxr59+4iOjlZug2zZsoUNGzYwbNgwwsPDG93OGTBggOrvSQghRNPqbmtIYtHWqUosBg0axL59+zhx4kSzMxfp6en88ssv3H///aoHFxcXR1FRkfJ1WlqakjCsXr1aSRhMjRswYAA7d+4kJSWFsrIyXF1dCQoK4tlnn21Q0js9PR2ou+Vz5MiRRuPavn17i+N2ubuS6/9x6Aud0BfYxy4JIYSwR3qjBqPZu0LkVoi9UfVON2PGDA4cOMBf//pXnnjiCWWLKUBFRQWHDh1i9erVdOrUiYcfflj14D799FOrxkVFRREVFXXDuLfeesuk9pojh5AJIYT5jEaHui2nZr3GRoMRqqlKLPr06cOiRYt47733+Oijj/joo4/QaDSkpKSQkpICgLOzM4sWLaJHjx5WHbAQQggh7JfqufmRI0fSr18/tm3bxtGjRykqKsJgMODr68vQoUOZNm0at956qzXHKoQQoh1TcwiZubdOhO1ZdNO/W7duPPXUU9YaixBCiA7MgMbsQ8jMrXshbE9WEwohhLALMmPRPliUWNTU1HDo0CFOnDjBpUuXAPD19WXgwIHcfffdODtb/xwMIYQQ7ZPB6IBGFm+2eaoTi2PHjvHee+9x+fJljNf9zX777bd07dqVBQsWcOedd1o8yLZGtpsKIYT5ZMaifVD1Tnfq1Clee+01amtrCQ4OJjIykm7dumE0GikuLiY1NZVTp07x17/+lbfeeovbb7/d2uO2a7LdVAghzCdrLNoHVYnF+vXr0ev1/OEPf2jyuPEHHniAxMREVq5cyYYNG3j99dctHqgQQggh7J+qQ8hycnLo169fk0lFvfvvv5/+/ftz6tQp1YMTQgjRccghZO2DqsRCo9GYVPiqR48eaDTyly6EEOLGJLFoH1TdCgkODiY/P/+Gcfn5+fTv319NF0IIIToYo4qzQsw9DVXYnqoZizlz5nDu3DnWr1+PwWBodN1oNLJhwwbOnTvHnDlzLB6kEEKI9k9mLNoHk2Ys6s//uNa4cePYvHkze/bs4Z577sHf3x+AoqIiDhw4QHFxMRMmTOCXX37pcLtCZLupEEKYz4AGc49NB426T8jCZkx6p3vvvfeaXCthNBopKipi69atyvVra1okJSWRnJzMuHHjrDTctkG2mwohhOioTEosZs2aJYswhRBC2JTBqDF/zYTcCjFZbW0tv/zyC1qtlrKyMtzd3fHy8qJnz544OVlvRt2klmbPnm21DoUQQoimSGJhfVqtll27dvHf//6XnJwcamtrG8U4OzsTHBxMeHg448ePx8vLy6I+5aa/EEIIu2AwoiKxsMlQ2rxz586xYcMGDh48qCQTXbp0oWfPnnh6euLm5kZ5eTmlpaWcPXuW7OxssrOzWb9+PXfffTePPvoot956q6q+JbEQQghhF2TGwjpWrVpFYmIiBoOBwYMHM3r0aAYNGkT37t2bfc358+fJyspi79697N+/nwMHDnD//fczb948s/tXnVhotVq+/fZbsrOzuXz5MjU1NU3GaTQaPvnkE7XdCCGE6CCMklhYRXJyMpMnTyY6OhpfX1+TXtO9e3e6d+/OhAkTuHTpEl988QXJyck3L7HIz8/nxRdfpLS0tNHJpkK2mwohhBpqt5uKhlavXk3Xrl1Vv97X15fY2FhmzJih6vWq3uk++eQTSkpKGDt2LA8//DDdu3fH1dVV1QCaU15ezqZNm8jLyyMvLw+dTkdMTEyjhaSmxgHk5eWxbt068vPz0el0uLi40LNnT6ZMmcLYsWNVtdkU2W4qhBCitViSVFijHVV1RU6dOkVgYCALFy4kMDDQ6kkFQElJCUlJSdTU1BAREWFxHEBZWRl+fn789re/5dVXX2XhwoV069aN5cuXs2nTJlVtCiGEsA6pvNk+qJqxcHNzU71a1FT+/v4kJCSg0WjQarUkJydbFAcQGhpKaGhog+dGjBjBhQsXSExMZObMmWa3KYQQwjpkjUX7oGrGYvDgweTl5Vl7LA1oNBqTinKZGtcST09PHB0drdqmEEII8xgxf7bCKGssTHbmzBmefPJJm/ej+hCy8vJy1qxZ0+QhZPbOYDCg1+vRarV88803ZGZm8sgjj1ivA0cjOJnwcJCFr0IIUa/+dFNzH8I0tbW1FBcX27wfVbdCevTowdtvv80bb7zBoUOHCA0NbXZLi0ajYdasWRYN0tri4uJITEwEwMnJidjYWCZNmmS19l3HV5kUV5vrRG2us9X6FUKItsygIlHQSGKhSEhIaPH6lStXbso4VCUWtbW1fPbZZ/zyyy8YjUZ+/fXXZmPtMbGYMWMGEyZMQKvVkp6ezqpVq6isrCQ6Otoq7VemuIHehB92A2CLWy4a+zjrT+Pm1tpDAEAXeVtrD4ELj5iWbNqa84nW/7Pos/xoaw8BAGN1dWsPAaNe39pDkN2a7UhCQgJdu3Zt9tyPpsp524KqxGL9+vWkpKTg7e3N6NGjbbLd1Jb8/f2VY97Dw8MBWLt2rVVqpAN1SYVsNxVCCLMYjXUP815kk6G0SbfccgtPPPEE9957b5PX8/LyWLhwoc3HoSqx2LNnD15eXnzwwQd4e3tbe0w3XXBwMDt27OD8+fPWSSyEEEKYzYD5izE1MuWiCAoKIi8vr9nEQqPR3JSilqoSi9LSUoYNG9YukgqArKwsHBwcWqyjLoQQwrZULcaUNRaKhx9+mIqKimav9+jRgzfffNPm41CVWAQEBHD16lVrj6WRjIwMqqqqlD+owsJC0tLSAAgLC1Nuv5gat2LFCtzc3AgODsbb2xudTkdaWhr79u0jOjq6wWyFqW0KIYSwDlm8aZmQkJAWr7u6ujaq5WQLqhKLhx9+mOXLl3Py5EkGDhxo7TEp4uLiKCoqUr5OS0tT3txXr16tvLmbGjdgwAB27txJSkoKZWVluLq6EhQUxLPPPtugpLc5bQohhLAOWWPRPqhKLG6//XamTJnCa6+9xkMPPcTQoUNbPEGtfqGkuT799FOrxkVFRREVFWXVNpvicnc5jQ4hK3BGXyhbS4UQojlyK8S6tFot//73v1mwYMFN7VdVYvHkk08qi0A2btzIxo0bW4z/6quvVA2urao+2Fl2hQghhGhV5eXlpKSk8PTTT+PgcPPKEKhKLEJCQqTktRBCCKuSGYv2QVVi8dZbb1l7HEIIITo4WyzezMvLY926deTn56PT6XBxcaFnz55MmTKl0dq6iooK1q9fz/79+ykpKaFXr15Mnz6dyMjIRu2aGmuLNu2dqsRCCCGEsDZbLN4sKyvDz8+PyMhIfH19qaysZO/evSxfvpyioiLlVGuAZcuWkZuby9y5c+nZsyd79+7lnXfewWAwMGbMmAbtmhprizbtnSQWQggh7IItboWEhoY22mI5YsQILly4QGJiopJYZGRkcPToURYvXszo0aOBupO8i4qKWLNmDaNGjVJOwTY11hZttgWqEosbHXRyLXs8K0QIIYQduol1LDw9PRvUYzp48CBubm6NqlZGRUXx7rvvkpOTo5RXMDXWFm22BaoTi5ZKg9Yv7DQajR0ysZDtpkIIYV8MBgNGo5HS0lL2799PZmYm8+bNU64XFBTQq1evRrMCgYGByvX6N3ZTY23RZlugKrFobk+s0WikuLiYzMxMfvzxR6ZMmUK/fv0sGmBbJNtNhRDCfEbMr3dlanxcXByJiYkAODk5ERsby6RJk5TrJSUlTR7r4OnpqVw3N9YWbZrrZpwNcj1VicX48eNbvB4TE8PmzZv57LPPmDhxoqqBCSGE6FjUrLEwNX7GjBlMmDABrVZLeno6q1atorKykujoaDVDbRM8PT2JiYm5qTUsAGzW24wZM/D19WXt2rW26kIIIUR7YlT5MIG/vz/9+/cnPDycP/7xj0ycOJG1a9ei1WqBujfhpmYF6p+rnzkwJ9YWbZrDw8ODmJgYs19nKZumMYGBgZw4ccKWXQghhGgn6rabasx8qOsrODgYvV7P+fPngbr3q7Nnz6LX6xvEFRQUANCnTx/lOVNjbdFmW2DTxOLXX3/FYDDYsgshhBDtRH0dC3MfamRlZeHg4KCsa4iIiKCiooIDBw40iNu1axc+Pj4EBwcrz5kaa4s21Th79iw//PADVVVVqtswh03qWJSWlrJp0yZOnz59U45otTeyK0QIIcxnizUWK1aswM3NjeDgYLy9vdHpdKSlpbFv3z6io6Px8vICIDw8nKFDh7Jy5UrKy8vp0aMHqampHDlyhEWLFjXYrWFqrC3aVGPr1q189913vPPOOw0SlCtXrrBr1y6MRiPDhw9XdqBYSvUhZM2prKykpKQEo9GIi4sLc+fOVT24tkp2hQghhH0YMGAAO3fuJCUlhbKyMlxdXQkKCuLZZ59tVNJ7yZIlrFu3jg0bNigltZ977rkmS2qbGmuLNs118uRJunXr1iCpqKmp4bnnnqO4uBij0cj69et57LHHmD59ukV9gcrEoqioqNlrjo6O+Pn5MWjQIB555BECAgJUD04IIURHolFxqFjL8VFRUURFRZnUkpubG7GxscTGxlot1hZtmuvSpUuNamCkpqZSVFRE3759iYyMZMeOHaxbt46BAwcSEhJiUX+qEott27ZZ1KkQQghxPTVrJlqhTEObU11dTefOnRs8l5aWhkaj4S9/+Qvdu3dn5MiRzJs3j+3bt7dOYiGEEEJYnS0rZHVgPj4+FBcXK19XVVWRlZXFwIEDlcWr/v7+3HHHHZw8edLi/m5u1QwhhBCiGeZvNVVxaFkHFBoayqlTp5StqykpKVRXVxMWFtYgzsfHB51OZ3F/Js1Y5OTkWNSJJdtkhBBCdBAyY2ET0dHRpKam8j//8z8MGjSI77//HgcHB0aNGtUgTqfTNbplooZJicXixYuVg8XU+Oqrr1S/ti2S7aZCCGE+W5b07sh69+7NkiVL+PDDDzl06BAajYZHH320wdkkBoOB3Nxc/Pz8LO7PpMQiJCTE7MQiJyeH6upqixKStkq2mwohhLAnYWFh/Otf/+LcuXO4u7vTtWvXBtePHj1KaWlpo2Pb1TApsXjrrbdMbjAjI4P4+Hiqq6sB9bdBysvL2bRpE3l5eeTl5aHT6YiJiWH27Nmq4gDy8vJYt24d+fn56HQ6XFxc6NmzJ1OmTGm0n7miooL169ezf/9+ZT/x9OnTLd5PLIQQohlyK8SmHBwc6NWrV5PXjEYj48aNY+TIkRb3Y7VdIUeOHCE+Pp7c3FyMRiP9+/cnJiaG8PBwVe2VlJSQlJREYGAgERERJCcnWxQHUFZWhp+fH5GRkfj6+lJZWcnevXtZvnw5RUVFzJw5U4ldtmwZubm5zJ07l549e7J3717eeecdDAYDY8aMUfU9CSGEaImGG9WlaPo1wlJhYWGNFnOqZXFikZmZSXx8PDk5ORiNRvr27cvs2bMZPny4Re36+/uTkJCARqNBq9U2mzCYGgd1K2OvLzE+YsQILly4QGJiopJYZGRkcPToURYvXszo0aMBGDx4MEVFRaxZs4ZRo0ZZVF5VCCFEE2TGol1QnVgcO3aMDRs2cOrUKYxGI7fddhuzZ89mxIgRVhmYqWszrLGGw9PTk6tXrypfHzx4EDc3t0b3mqKionj33XfJyclpVMVMCCGEhSSxsCqtVsu///1vFixYcFP7NTuxyMrKIj4+npMnT2I0GgkKCiImJoaIiAhbjM8mDAYDRqOR0tJS9u/fT2ZmJvPmzVOuFxQU0KtXr0azEvUHtBQUFLScWDia+JNuAAwyjSeEEEBdOW9zd3nIrpBmlZeXk5KSwtNPP42Dw80rW2VyYnH8+HHi4+M5ceIERqORwMBAYmJiuPvuu205PpuIi4sjMTERACcnJ2JjY5k0aZJyvaSkpME2nHqenp7K9Za4jqswaRy1uc7U/uRi6rBNpnG2k4Kqfl1vHHMT7PpgRWsPgU4a+9hqPND5sdYeAg5eXVp7CADoL19p7SGAwQ4+bmvsYAyiXTHpHejFF18kOzsbgICAAGJiYrjnnntsOjBbmjFjBhMmTECr1ZKens6qVauorKwkOjraKu1XpriB3oQs2mCV7oQQol2Qs0LaB5MSi+PHj6PRaHBxccHHx4fk5OQWF0leS6PR8Oqrr1o0SGvz9/fH398fQNm1snbtWsaPH4+Xlxeenp5NzkrUP1c/c9EsvUbqWAghhLlkjUW7YPKcudFopKqqiszMTLM6aAsFsoKDg9mxYwfnz5/Hy8uLwMBAUlNT0ev1DdZZ1NdZ79OnT2sNVQgh2i9ZY9EumJRYvPnmm7YeR6vKysrCwcFBWVcRERFBUlISBw4caFBLfdeuXfj4+MjZJ0IIYQMao/lLPmSJiP0xKbG4vvbDzZKRkUFVVRUVFXWLIQsLC0lLSwPqinm4urqaFbdixQrc3NwIDg7G29sbnU5HWloa+/btIzo6Gi8vL6Du9sjQoUNZuXIl5eXl9OjRg9TUVI4cOcKiRYukhoUQQtiC3AppF+xk+0DT4uLiKCoqUr5OS0tTEobVq1crCYOpcQMGDGDnzp2kpKRQVlaGq6srQUFBPPvss41Kei9ZsoR169axYcMGpaT3c889JyW9hRDCVuRWiNUZW2F1q10nFp9++qlV46KiooiKijIp1s3NjdjYWGJjY02Kv5acbiqEEKK1eXp6EhMTc1NrWICdJxZtlZxuKoQQKsitEKvy8PAgJibmpvcriYUQQgj7IIlFuyCJhRBCCPsgiYXNlJSU8OWXX3L06FFKS0txd3dn3LhxPPjgg1bvSxILIYQQdkLF4k05Nv2GLl26xPPPP8/FixcxGo24ublRVFTE6dOnlZiTJ09SWVnJoEGDcHa2bD3gzV3RIYQQQjSjvo6FuQ/RsrVr11JcXExUVBTx8fFs2rSp0W6R6upqli5dyt69ey3uT9WMxZkzZ+jdu7fFnbdXsitECCGEvThy5Ai33norf/rTn5qthj1kyBC8vLxIT083efdkc1QlFvPnz8fLy4uQkBBCQ0MZNGiQlLm+huwKEUIIFWSNhU2UlZUxaNCgGx6x0aNHD/Lz8y3uT1ViERYWxsmTJzlw4AAHDx4EoEuXLoSEhDBo0CBCQ0Ml0RBCCCHsgI+PD5cuXbphnK+vb+slFq+++ioGg4G8vDyys7PJyspSEo0DBw6g0Wjw8PAgJCSEwYMHM3XqVIsHKoQQon2Ts0JsY/DgwaSkpFBQUNDih/6Kigpqa2st7k/1rhAHBwf69etHv379mDZtGkajkby8PI4fP052djaZmZkcPnyYw4cPS2IhhBDixqSkt0089NBD7N69m//93/9l6dKl+Pv7N4qpqqrip59+wtfX1+L+rLYrpLi4mIKCAgoKCsjPz6empgaj0SgHdgkhhDCNUeVDtKhPnz48+eST/PLLLyxYsIC1a9c2uF5VVcXKlSvR6XQMHTrU4v5Uz1gUFRWRnZ3N8ePHOX78OMXFxRiNRpycnOjfvz9jxowhNDSUAQMGWDxIIYQQHYAs3rSZKVOm4OPjw8qVK/n8888B2LdvHz/88APFxcXo9Xo8PT35zW9+Y3FfqhKLp556SjlN1NHRsVEi0alTJ4sH1pbJdlMhhBD25u677+bOO+8kMTGRQ4cOkZ+fz/nz53FxcSEsLIzHH38cPz8/i/tRlVhcuHABjUZDQEAA06dPJywsDA8PD4sH017IdlMhhDCfLN60PVdXV6ZNm8a0adMAqK2txcnJukW4VbU2efJksrOzKSwsZPny5QAEBgYqNS0GDRokiYYQQgjzyK2Qm87aSQWoTCx+//vfA6DT6ZRdINnZ2Wzfvp1t27ah0WiURCM0NJS77rrLqoMWQgjRDkli0S5YlKp06dKFkSNHMnLkSKDu9LTs7Gy+//57du/eTX5+Ptu3b+err76yymCFEEK0X3IrxDqsdeyG2nasMgdSU1PDqVOnlNmLH3/8kZqaGoAblhAVQgghAKljYSV/+tOfGDVqFDNmzFBVBTsvL4/PP/+cAwcOsHXrVrNfryqxaCqRqK2tVU5L8/X1VUp7h4aGqulCCCGEECrMmjWLLVu2sG/fPgIDAxkzZgyDBg0iKCioyTUVNTU1/Pzzzxw/fpy9e/dy5swZOnXqxKxZs1T1ryqxiImJUQpgAdxyyy3Kos3Q0FC6d++uajDthWw3FUIIFWSNhVXExMQwadIkPvvsM1JSUlizZg0ajQZHR0f8/f3x8PDAzc2NiooKSkpKKCoqwmAwYDQa6dy5Mw888AAzZszAy8tLVf+qEgsvLy8liRg0aJBNEony8nI2bdpEXl4eeXl56HQ6YmJimD17tqo4gGPHjrFnzx5OnjzJxYsXcXd3p3///syaNYt+/fo1iM3JyWH9+vX8+OOPGI1G+vfvz5w5c7jjjjtuOHbZbiqEEOaTNRbW4+3tTWxsLHPnzmX//v3897//5eTJk5w7d65RbNeuXbnjjjsYPnw49957Ly4uLhb1rSqx+PTTTy3q1BQlJSUkJSURGBhIREQEycnJFsUB7Nixg5KSEh588EF69+6NTqdjy5YtLF68mNdee40hQ4YAdUnFCy+8QHBwMAsXLgTgiy++4KWXXmLZsmVSTVQIIWxBZiysrlOnTowfP57x48cDoNVquXr1KuXl5XTu3Blvb2/VMxPNsf4GVivx9/cnISEBjUaDVqttNmEwNQ7qtsl6e3s3eG7YsGHExsayefNmJbHYsGED7u7uLF26FFdXVwCGDBnCU089xb/+9S/efvttK32XQggh6smMhe15eXlZPZG4nkWJRUFBAd988w0nTpzg8uXLQN257yEhIUyePFnVatR6pu4mMWfXyfVJBYCbmxsBAQFcvHhRee7kyZOEh4crSQVA586dCQkJ4eDBg1y+fBkfHx+T+xVCCGECmbFoF1QnFtu2bWPNmjXKgo96paWlFBYWkpyczBNPPMGDDz5olYHaSllZGT///DODBw9WnqupqcHZufFCy/rn8vPzW04sHE38STcABlmLIYQQgCQW7YSqxCIzM5PVq1fTqVMn7r//fsaNG4e/vz8ajYYLFy6we/duEhMT+fTTT+nTp49yi8Eeffzxx1RWVjY40S0gIIBTp05hMBhwcKg7WV6v15OTkwPUretoieu4CpP6rs11pvYnyxbJCCGEEPZEVWKxdetWHB0def311xk4cGCDa0FBQQQFBXHPPffwwgsvsGXLFrtNLNavX8+ePXuYN29eg10hU6dO5YMPPuDjjz9m5syZGAwGEhISlBNd65ON5lTt7gx6E2YiDKBxsP6MhcbFPra1GjzcWnsIAHTS2Mefhz1wSfNs7SFgrDzb2kMQdkqDijUWNhmJsISqxCI3N5dBgwY1SiquNWDAAEJDQ5VP+fYmISGBTZs28dhjjzF16tQG1+677z60Wi2fffYZO3bsAOq+n4cffpgvvvjixusr9BrTEgshhBCinVGVWFRVVdGlS5cbxnXp0oWqqio1XdhUQkIC8fHxzJ49u8EtkGtNnz6dhx56iHPnzuHm5oa/vz8rVqzA1dW1Uc0LIYQQViBrLNoFVYmFn58fP/74I3q9HkdHxyZj9Ho9P/74I35+fhYN0No2btxIfHw8M2fOJCYmpsVYZ2dnZWdLUVER+/fvZ8KECXTq1OlmDFUIIToU2W7aPqhKLO666y62bt3Khx9+SGxsLJ07d25wvby8nH/+859cvHiRadOmqR5cRkYGVVVVVFTULYYsLCwkLS0NgLCwMGU7qKlxW7ZsYcOGDQwbNozw8HB+/PHHBv3VF74qKCjgwIED9OvXD2dnZ06fPs3nn39Ojx49mDNnjurvRwghRAtsMGNhTsXliooK1q9fz/79+ykpKaFXr15Mnz6dyMjIRu2aGmuLNs114cIFunXrZlJseno6I0aMsKg/VYnFjBkzOHjwILt37+bQoUOEh4c32BWSkZFBeXk53bt3Z8aMGaoHFxcXpyyYBEhLS1MShtWrVysJg6lx6enpABw5coQjR4406m/79u0AODk5cezYMbZv305FRQW33HILkyZNYvr06Q1qWwghhLAiGyQWplZcBli2bBm5ubnMnTuXnj17snfvXt555x0MBgNjxoxp0K6psbZo01wLFiwgNjaWcePGNRtTXV3N6tWrSUpK4quvvlLdF6hMLDw9Pfnb37DgG+kAACAASURBVP7GRx99REZGBqmpqY1iwsPDmT9/Ph4eHqoHZ2rpcFPj3nrrLZPievbsyd/+9jeTYpviHFFGo0PICp0xnJGtpUIIcTOZWnE5IyODo0ePsnjxYkaPHg3A4MGDKSoqYs2aNYwaNUq59W9qrC3aVKOmpob333+fjIwM/vjHPzZ6X/7pp5/4+9//zi+//GKVs79UF8jy9fXllVde4fz5840qb95xxx0d+oTTmkPusitECCHMpWKNxY1mLEytuHzw4EHc3Ny49957G8RGRUXx7rvvkpOTo+yENDXWFm2q8d577/Huu++yf/9+Tp48ycKFC5WikJs3byYhIYHa2lqioqKIjY1V1ce1LD4rpHv37h06iRBCCGElN2lXSFMVlwsKCujVq1ejWYHAwEDlev0bu6mxtmhTjd69e7N8+XLWrVvHli1bePnll5k6dSo///wzJ06cwNPTkz/96U/cfffdqtq/XsuVnix09epV/v3vf9uyCyGEEO1E/a4Qcx/maqricklJCZ6ejQvI1T93bcVlU2Nt0aZajo6OPP7447zxxhu4ubnx9ddfc/LkSYYMGcKKFSusllSAjU43LS4u5ssvv+S7776jpqaGxx9/3BbdCCGEaE9uwoxFcxWXO4Ly8nKSk5MpLy9Xnjt79iyFhYV07drVav2YnFgYDAZSU1PJzMzk6tWreHt7ExYWxr333quUuC4uLiYhIYHdu3djMBgAiIiIsNpghRBCtGM2Tixaqrjs6enZ5KxA/XPXziaYGmuLNtXKzs7mH//4B8XFxdx22238+c9/JjU1lS+++IJXXnmFBx54gLlz5zZ5AKe5TEos9Ho9S5cuJSsrq8FJpnv27GH//v0sWbKE7777jn/+859UV1cDdbUuYmJiCAoKsniQQgghhCVuVHE5MDCQ1NTURoUfCwoKAJRiiebE2qJNNf7zn/+wZcsWjEYj0dHRzJkzBycnJwIDAwkLC2P58uVs376dY8eOsWjRImVdh1omrbH4+uuvOXbsGE5OTtx///3MmzePxx57jIEDB3L48GFWrFjBihUrqKqqYujQobz33nssWbKkwyYVzhFlOI9s+HDoXd3awxJCCLtmqzUWplRcjoiIoKKiggMHDjR4fteuXfj4+BAcHGx2rC3aVOOLL76ga9euvPHGGzz++OM4Of3fnEJISAgffvghkZGRFBQUsGjRItX91DNpxmLfvn04ODjw1ltvNfjmZsyYwcqVK0lMTESj0fD4448THR1t8aDaOtluKoQQKtjgVoipFZfDw8MZOnQoK1eupLy8nB49epCamsqRI0dYtGhRg1kEU2Nt0aYa9957b5P1K+p17tyZRYsWMXz4cOLi4lT3U8+kxOLs2bMMGDCgyYwpOjqaxMREevbsKUmFEEII9WyQWJhacRlgyZIlrFu3jg0bNigltZ977rkmS2qbGmuLNs31/PPPmxQXGRnJHXfcYVFfYGJiUVFR0Wyd8frnO+ptDyGEENahUZFY3OhWiKkVl6GucFZsbKxJRaJMjbVFm7ZkjYNDTUosjEajsvPjehpN3ZS/i4uUqxZCCGEBOTbdJrKzs82KHzRokEX92aSOhRBCCGEuW8xYiLpbLPWTAKa4aYeQpaSkkJKS0uQ1jUbT4nVLB9nWyCFkQggh7MXYsWObTCyMRiMXL17k559/pry8nBEjRlh0cGg9kxOLa+tXiJbJrhAhhFBBboXYxMKFC1u8XlJSwocffkhhYSHvvvuuxf2ZlFhs27bN4o6EEEKIFkli0So8PT159tlniY2N5T//+Q/z58+3qD2bHkImhBBCmEqj8iEs5+rqSnBwsLI91xKyeFMIIYT9kBmIVlNRUUFpaanF7UhiIYQQwi7IrpDWk56ezg8//EDv3r0tbksSCyGEEKIde//995u9VlFRwblz5ygoKMBoNPLwww9b3J8kFjYg202FEEIFWbxpE7t27bphzC233EJMTAzjxo2zuD9JLGxAtpsKIYQKkljYxJtvvtnsNWdnZ7p27drssR1q2G1iUV5ezqZNm8jLyyMvLw+dTkdMTAyzZ89WFQdw7Ngx9uzZw8mTJ7l48SLu7u7079+fWbNm0a9fvwaxP//8MwkJCeTm5lJaWsott9zC6NGjefjhh3F1dbXp9y6EEB2RrLGwjdDQ0Jvan91uNy0pKSEpKYmamhoiIiIsjgPYsWMHRUVFPPjgg7z66qvExsZy9epVFi9ezLFjx5S4wsJCnn/+eYqKinjyySd55ZVXiIyMZOPGjVYpHiKEEKIJ/397dx4XZbk3fvwzbLKDGG6JIu6aoIhmnmMmmT25oJUI+GTlr1xOm3qyXm6VZi7lcTvp0Uf0aCmMHH2k8jGXk4oamkikqKEiJmoqCKgMO8L8/uAwOQIyq9zC9/16zau455rr/t6XA/Od674WrYkPoSiK7bFo2rQparUalUrFnTt32Lt3r1nlACZNmoSnp6fescDAQCZMmMDWrVsJCAgA4ODBg5SUlDBjxgxatGgBQEBAADk5OezZs4e8vDyLLHsqhBDiD9JjUT8YlFjUtAeIoUwZDGLohinGbKxyf1IBFdvUtm7dmqysLN0xO7uKZnF2dtYr6+rqio2Nje55IYQQFiRjLCzizTffNPm1KpWKyMhIs85v0Cfk8uXLjfoAr6TValGpVBYZZWot+fn5pKWl4e/vrzsWHBzMt99+y+rVq3n99ddxd3fn9OnT7N69myFDhtQ+xsLWwHd6OaCVQZ5CCCEsJzMzs07Pb1BiER4eblJi8ShYs2YNRUVFjB49WnesWbNmLF68mAULFjB+/Hjd8eHDh+v9XJNGAwsMOvfdC/aUpTUyPuhaaEvvWrxOU9hcu1nXIQDw7Ng36joEbrdXxlTjx9Wn6zoEyouL6zqECmVldR0BaMvrOgJQ0AaTcivEMup6fy+DEovqZljUB5s3byYuLo6JEyfqzQrJyMhg3rx5eHp6Mn36dDw8PDh//jwxMTEUFRXx3nvvPbDe4gPOhk03VcDfFCGEUAy5FVIvNNjBAmq1mpiYGMaOHcuwYcP0nvvqq68oLCzk73//u+62xxNPPIG7uzsrVqxg4MCBD56+U6aSdSyEEMJYklhYxKxZswgMDOTll1+u8lxBQQF2dnY4OFivF9Xs6aYajYZffvmFgwcPkpKSYomYrE6tVhMdHc2YMWP0boFUunjxIj4+PlXGUnTo0AGomI4qhBDCslRa0x5C36lTp7h69Wq1z0VERLBmzRqrnt/kHotbt26xdu1ajh49ivY/9+iCg4Pp0qULADt37mTz5s3Mnj2bbt26WSZaC9iyZQvR0dGEhYURERFRbZkmTZqQnp5OYWEhTk5OuuNnz57VPS+EEMLCpMfC6rRare4z21pMSizu3LnDhx9+SEZGBn5+fnTp0oWdO3fqlenbty+RkZHEx8ebnFgkJiZSXFxMYWEhUNFTEB8fD0CvXr10PQqGlouNjSUqKorAwECCgoJ0iUKlzp07AxASEsL8+fP56KOPGDFiBO7u7pw7d45t27bh4+NDr169TLoeIYQQor4zKbGIiYkhIyOD//7v/yYsLAygSmLRpEkTfHx8OHPmjMnBrV69Wm/aTHx8vC5hWLdunS5hMLRcQkICAElJSSQlJVU5344dOwB48skn+eyzz9i2bRuRkZHk5+fj7e3N888/T2hoKPb29iZfkxBCiOqptFqjZ6moFDSrRVQwKbE4duwYrVq10iUVNfH29ubcuXMmBQawfv16i5ZbuHChwef29/fXW9vCGLK7qRBCmEBuhdQLJiUWt27d4sknn6y1nIODg+72REMiu5sKIYTxZB2L+sGkxMLFxYXs7Oxay127do3GjRubcgohhBANjfRYWMz+/fur3Y5DpVLV+Fylb7/91qxzmzTdtHPnzqSmppKenl5jmV9//ZVLly4pakaIEEII5ZLpppZTOfvDlIe5TOqxGDlyJAkJCXz22We88847VRaLOnPmDMuWLcPW1pYRI0aYHaQQQogGQhIFsz0SS3rfr1u3brzxxhusX7+ejz/+GCcnJ1QqFUePHiUhIYG8vDwAJkyYQLt27SwasBBCCCGUy+QFskJCQujYsSPbtm0jOTkZrVZLQUEB9vb29OjRg9GjRzfY2yAyK0QIIYxnym0NuRWiPGbtFdK5c2dmz56NVqslNzeX8vJy3N3dsbW1tVR8jySZFSKEECYwJUmQxEJxLLIJmUqlwsPDwxJVCSGEaKCkx6J+MGlWSGxsrEHlNBoNixYtMuUUQgghGprKlTeNfQhFManHYsOGDSQlJTFlypQaN+Q6efIky5cvJycnx6wAhRBCNAzSY1E/mNRj0b17d06ePMm7776r25Oj0t27d3WzRXJycnjxxRctEqgQQgghlM+kHov58+ezfft2Nm/ezBdffMHAgQOZOHEimZmZLFmyhEuXLuHt7c2UKVOqrHEhhBBCVEsGb9YLJg/efOmll+jZsyeLFy/mwIEDnDx5ktzcXEpLSxkwYACTJk3CxcXFkrE+MmS6qRBCGE9Vzv1/Omt/jSQWimPWrJC2bdvy8ccfM2XKFLKzs1GpVAwcOJCpU6daKr5Hkkw3FUIIE0iPRb1gVmJx9OhRVq1aRUFBAe3atePy5cvExcUBMHHiRJydnS0RoxBCiIZAa3SHhSQWCmRSYlFUVMTatWvZt28ftra2vP7667z44oukp6ezZMkSDhw4wJkzZ5g6dWqDXX1TCCGEkbQmbG8q000Vx6RZIZMnT+aHH37g8ccfZ8mSJbz00kuoVCp8fX1ZunQpw4cP5+bNm8yaNYuvv/7a0jELIYSoh2R30/rBpMTi+vXrDB06lOXLl9O2bVu95+zt7Rk/fjxz5szBw8OD//3f/7VIoEIIIYRQPpNuhXz88ccEBQU9sEzPnj358ssvWbVqlUmBCSGEaGBk8Ga9YFJiUVtSUcnd3Z0ZM2aYcopHmkw3FUII45l0W0MSC8WxyCZk1lBQUEBMTAwXL17k4sWL5ObmEhERwZgxY0wqBxXLjMfFxZGSkkJWVhYuLi506NCB8PBw2rdvryu3bNky9u/fX2NsixcvpnPnzjU+L9NNhRDCBDJ4s14wKLGo/JDt27cvzs7OD/zQrU5wcLDRgWk0Gvbs2YOvry99+/Zl7969ZpUD2LVrFxqNhpCQEHx8fMjNzSU2NpZp06Yxd+5cAgICAAgPD+eFF16o8vp58+Zhb29Phw4djL4eIYQQDyY9FvWDQYnF8uXLUalUdOrUCWdnZ93PtdFqtahUKpMSi6ZNm6JWq1GpVNy5c6fGhMHQcgCTJk3C09NT71hgYCATJkxg69atusSiRYsWtGjRQq/cqVOnyM3NJSwsDFtbW6OvRwghRC0ksagXDEoswsPDUalUuLu76/1sTYbWb0wc9ycVAE5OTrRu3ZqsrKwHvvbf//43KpWKQYMGGXw+IYQQhlOBJAr1gEGJxf3jFaobv/Coys/PJy0tDX9//weWiY+PJyAggObNm9deqa2BvxnlgFbGYgghhKg/FDt482FZs2YNRUVFjB49usYyhw4doqSkhOeee86gOhsNLDCo3N0L9pSlNTKorBBC1HvlJgzeRPvAFZmMGeBfWFjI5s2b+fHHH9FoNLRq1YpRo0bx9NNPm1zWGnUqnVGJRWJiIj/99BM3b97E3t4eX19fBg0aZNi3eAXavHkzcXFxTJw4UW9WyP327t2Lm5sbTz31lEH1Fu13MmxWSDn/+UWysJJSy9dpCk1eXUcAQKPjqXUdAk1/KqvrEAAoV8B7Q1umjLZAW17XEYj7WeHPoTED/BcsWEBqaiqvvfYajz/+OAcPHmTx4sWUl5fzzDPPmFTWGnUqncGJxd/+9jcOHz4MVAzKBDh+/DixsbF8+OGHPPnkk9aJ0ErUajUxMTGMHTuWYcOG1Vjut99+48KFC4SEhGBvb29Y5WUquCu3OIQQwhjWWJ7b0AH+iYmJnDhxgmnTpjFgwAAA/P39yczMZMOGDfTv3183cN/Qstao81Fg0JLee/fu5dChQ9jY2PDss88yYcIExo4dS6dOnSgtLWXZsmXk5+dbO1aLUavVREdHM2bMmAfeAoGKQZsAgwcPfhihCSFEw6XVmvZ4AJVKZdAg/6NHj+Lk5MSf//xnveODBg0iJyeH8+fPG13WGnU+CgxKLPbv349KpWLOnDm89957DBs2jNDQUL744guCg4MpLCzk6NGj1o7VIrZs2UJ0dDRhYWFEREQ8sGxpaSlxcXF07NiRNm3aPKQIhRCiYarLTcjS09Np1apVlV4BX19f3fPGlrVGnY8Cg26FXLp0iU6dOunWebjX6NGj2b9/P5cuXbJ0bCQmJlJcXExhYSEAly9fJj4+HoBevXrh6OhoVLnY2FiioqIIDAwkKCiIs2fP6p3v/tU0f/rpJzQaDa+99prFr00IIYRyaDSaascLurm56Z43tqw16nwUGJRYFBYWVlkwqlJlQxQUGDYTwhirV68mMzNT93N8fLwuYVi3bp0uYTC0XEJCAgBJSUkkJSVVOd+OHTv0ft67dy+Ojo7079/fglclhBCiWrKGRb1gUGKh1Wqxsan+rknlca0V1mtfv369RcstXLjQqPPPmzfPqPKVHJ4qoMomZOn2lF02cPCnEEI0QCpT9gqxUDbi5uZWba9A5bHKngNjylqjzkeBQWMshHFKjjpTclj/IUmFEELUotzEhwX4+vpy9epVyu6bDl05tuHecXaGlrVGnY8Cg6eb7t+/v8bNx1Qq1QOf//bbb02LTgghRINRlz0Wffv2Zc+ePRw5ckTv9ve+ffvw8vKiY8eORpe1Rp2PAoMTC2vc6hBCCCF0rPQxY8gA/6CgIHr06ME//vEPCgoKaNGiBYcOHSIpKYn3339fb7aGoWWtUeejQKWVjMFiCgoKCAsLo+jfznW6QJZKIW9Alb0yVoxX2dV9HEpZbVIrK2/+QVberGCnxXFwETExMTg7O9dJCJV/Ows0vTF+p4m7OLsdf2D8b7zxht4A/3utW7eOZs2aARUTFTZt2qS3pHZoaGiNy28bUtYadSqdJBYWJImFPkks/qCUD1NJLO4hiUWFBpBYiIer7v/iCiGEEFhnSW/x8EliYQUy3VQIIUxQh4M3heVIYmEFJUfr9laIEEI8ilRyd6pekMRCCCGEMkiPRb0giYUQQgjlMDZPkM5hxZGVN4UQQghhMdJjIYQQQhFUWv5zO0Q8yiSxEEIIoQxarQmJhSQiSiOJhRXIdFMhhDBBOTLGoh6QxMIKZLqpEEIYTyU9FvWCJBZCCCGUQRKLekESCyGEEMogiUW9INNNhRBCCGEx0mMhhBBCGWTwZr0giYUVyKwQIYQwngzerB8ksbACmRUihBAmkMSiXlBsYlFQUEBMTAwXL17k4sWL5ObmEhERwZgxY0wqB3Dy5Eni4uJISUkhKysLFxcXOnToQHh4OO3bt69S/syZM2zdupWzZ89SWlpKkyZNCA4OJjw83GrXLYQQDZYkFvWCYgdvajQa9uzZQ2lpKX379jW7HMCuXbvIzMwkJCSETz75hAkTJnD79m2mTZvGyZMn9crGxcUxc+ZMXFxc+Otf/8onn3zCyy+/jFaWmxVCCOuoTCyMfQhFUWyPRdOmTVGr1ahUKu7cucPevXvNKgcwadIkPD099Y4FBgYyYcIEtm7dSkBAAADZ2dmsWrWK559/nrfeektX1t/f3wJXJoQQQtRfik0sVCrDxigYWg6oklQAODk50bp1a7KysnTH9u7dS1FREaNGjTK4biGEEGaSWSH1gmITi4clPz+ftLQ0vd6I06dP4+bmxtWrV/nss89IT0/Hzc2Np556inHjxuHs7PzgSm0N/M0oB8rlt0IIIUBmhdQXDT6xWLNmDUVFRYwePVp3LDs7m+LiYhYtWkRoaCjjx48nNTWVqKgo0tPT+fzzzx/YU+IYXGjQue+m2nP3goPZ1yCEEPWCJBb1QoNOLDZv3kxcXBwTJ07UmxWi1WopKSnh1VdfJTQ0FIDu3btjZ2dHZGQkJ0+epEePHjXWW7SvEZQZ0BNRDmjLzb2MKrRlFq/SNFa4NlNoS0rqOgS05Qr546eQfxNFkEF/FZTUDlotGPu7YqOg+AWg4Fkh1qZWq4mJiWHs2LEMGzZM7zk3NzegYmDnvXr16gVAWlragysvU1WsY1HbQ26DCCHEH2RWSL3QIBMLtVpNdHQ0Y8aM0bsFUsnX17fa11VONTVmwKgQQggDSWJRLzS4xGLLli1ER0cTFhZGREREtWX69esHwM8//6x3PDExEYBOnTpZN0ghhBDiEaXoMRaJiYkUFxdTWFgxGPLy5cvEx8cDFbclHB0djSoXGxtLVFQUgYGBBAUFcfbsWb3zde7cGai4BdKnTx+2bNmCVqulU6dOpKamsmXLFnr37k23bt2sf/FCCNHQmNIDIT0WiqPoxGL16tVkZmbqfo6Pj9clDOvWrdMlDIaWS0hIACApKYmkpKQq59uxY4fu/z/88EPUajW7d+9GrVbj5eXFiBEjauzlEEIIYaZyEwZvyqwQxVFpZY1qiykoKCAsLIxyDVTZ3fSyHWXpDymPUynjDpfKRsaiVJJZIQokf/oq2GlxfL6YmJiY2tfosZLKv513f/WDclvjXmxThl3Xi3Uav9Cn6B6LR1XJUUfZ3VQIIYylxYRbIVaJRJhBEgshhBDKILdC6gVJLIQQQiiDDN6sF5RxM14IIYQQ9YL0WAghhFAG6bGoFySxsAKHp4qo01khQgjxKJLEol6QTzorkFkhQghhgvLyis0ZjXuRNSIRZpDEQgghhDJIj0W9IImFEEIIZZDEol6QWSFCCCGEsBjpsRBCCKEMWhMWyFJJj4XSSGIhhBBCEbTacuMX0pT9bxRHEgsrkOmmQghhAlOW9JYeC8WRTzorkOmmQghhAhm8WS9IYiGEEEIZTFnHQiW3QpRGZoVYUGlpacX/SNcc2GixbVcsbQGg0mLXvgRspC2w0WLXoVTaAqQtqlPZY2HsQyiKJBYWpEsspFXBBuzal0pbQEVbdJC2AP7TFnelLaCiLTqWSVuIekduhQghhFAEbbm24naIMaRXVHEksRBCCKEMWq0J000lsVAaSSyswOHJYtDKdFMhhDBKudaEwZuSWCiNYj/pCgoKiImJ4eLFi1y8eJHc3FwiIiIYM2aMSeUATp48SVxcHCkpKWRlZeHi4kKHDh0IDw+nffv2unKnTp1i5syZ1ca1ePFiOnfu/MDYS441gmLTb5zatrlrVhJi27qUssv2Jr/eUnVYgo1PCeVXHOq0DkvEYAkW+Xc1971l5ustxRJxKKUOcykhBovRllf5Ulb7aySxUBrFvhs1Gg179uzB19eXvn37snfvXrPKAezatQuNRkNISAg+Pj7k5uYSGxvLtGnTmDt3LgEBAXrlX331Vbp37653rE2bNuZfXC1sW5v7x98SH0DKSCxsW5ea/aFubh2WiMESLPLvau57y8zXW4ol4rBtU6aIOsylhBgsRSs9FvWCYt+NTZs2Ra1Wo1KpuHPnTo0Jg6HlACZNmoSnp6fescDAQCZMmMDWrVurJBYtW7astXdCCCGEEH9QbGKhUhnWHWZoOaBKUgHg5ORE69atycrKMrgeIYQQViC3QuoFxSYWD0t+fj5paWn4+/tXeW7NmjV88cUXNGrUiM6dOxMWFka3bt3qIEohhKj/tFoTboXIAmOK0+ATizVr1lBUVMTo0aN1x5ydnQkJCeGJJ57A3d2d69evs337dmbOnMknn3xCYGBgtXVpKzNnBwN/M8qpITvXgp05vyzmvt4Cddhq9f9rThx1XoeZr7dkW9T1v6uS3hd13hZm1qGUtvjP+bVK+OZva8LuprZWiUSYoUEnFps3byYuLo6JEyfqzQpp164d7dq10/3crVs3+vbty7vvvsuGDRtqTCzu3r0LgGP/UrNjcxxcZN7rnyswPwYL1NFooAXqGJRf53VYIgbH4ELz67DEv6u57y0zXw/g+Gyx+XVYIo7nLRCHmXU4Diqp8xgACgsLcXFxMbseU9jZ2dG4cWNuPXvLpNc3btwYO7sG/XGmKA32X0KtVhMTE8PYsWMZNmxYreVdXV3p3bs3u3btori4mEaNGlUp06xZM9auXYudnZ1BYz/s7e2xt6/7mRdCiIZLq9VSWFiIl5dXncXg4ODAunXrdF/OjGVnZ4eDQ93P3BIVGmRioVariY6OZsyYMXq3QGpT2VVYU9JgY2NDixYtLBKjEEI8LHXVU3EvBwcHSQ7qiQa3/c2WLVuIjo4mLCyMiIgIg1+Xl5fH8ePH8fPzkze/EEIIUQNF91gkJiZSXFxMYWHFvenLly8THx8PQK9evXB0dDSqXGxsLFFRUQQGBhIUFMTZs2f1zle5ZsXixYvx9vamQ4cOuLu7c+3aNWJjY7l9+zZTpkyx/oULIYQQjyiVVhFDgav3xhtvkJmZWe1z69ato1mzZkaVmzFjBqdPn67xfDt27ABg69at/Pjjj2RkZFBYWIibmxtdu3Zl1KhRdOzY0ZxLEkIIIeo1RScWdSktLQ21Wk1qaip5eXl4e3szYMAAXnzxRV0PiDH7lNxvz549rFy5EkdHR7Zu3WrtyzGbtdrjzJkzbN26lbNnz1JaWkqTJk0IDg4mPDz8YV2a0azRFobUqUSGxG3oHj2VCgsL2bx5Mz/++CMajYZWrVoxatQonn766Yd9eUaxdFsY225CKIXtnDlz5tR1EEpz+fJlpk2bhq2tLWPHjmXw4ME4OzsTExNDenq67g9cTk4OK1euxN3dna5du5KWlkb37t2r7C9yv+zsbObPn4+bmxulpaVGDSCtC9Zqj7i4OObPn0/79u2JiIhg0KBB+Pj4UFhYWGsb1hVrtIWhdSqNoXFv3LiRW7du8dxzz+nWhzlx4gRbtmyhS5cuNG/eXK/eTz/9lMTERF555RVCQkLQaDRs2rSJli1b4uvrOULhzAAAE75JREFUWwdXWjtrtIWx7SaEUih6jEVdOXjwICUlJcyYMUM3yyMgIICcnBz27NlDXl4erq6uRu1Tcq9Vq1bRrVs3XF1dOXLkiDUvxSKs0R7Z2dmsWrWK559/nrfeekt3vLoVUJXEGm1haJ1KY2jcxuzRk5iYyIkTJ5g2bRoDBgwAKt4TmZmZbNiwgf79+2Nrq7wVkazRFsbubSSEUjS4WSGGqFxoxdnZWe+4q6srNjY2uudVKpVRe5UAHDhwgNOnT/OXv/zFMsE+BNZoj71791JUVMSoUaMsG6yVWaMtDK1TaQyN25g9eo4ePYqTkxN//vOf9Y4PGjSInJwczp8/b8lLsBhrtIXsbSQeVZJYVCM4OBgXFxdWr17NjRs3KCgoICEhgd27dzNkyBCT73nfvn2byMhIXnvtNR577DELR2091miP06dP4+bmxtWrV3nvvfcYMWIEr7zyCqtWraKgwPzVJa3FGm1hrfebtZkTd+UePa1bt9Y7np6eTqtWrar0SlTeAklPT7f4dViCNdrC3LJC1BVlfhWqY82aNWPx4sUsWLCA8ePH644PHz5c72djrV69mlatWjFkyBBLhPnQWKM9srOzKS4uZtGiRYSGhjJ+/HhSU1OJiooiPT2dzz//3OjeoIfBGm1hrfebtZkTd3V79ABoNJpqxw64ubnpnlcia7SFuWWFqCuSWFQjIyODefPm4enpyfTp0/Hw8OD8+fPExMRQVFTEe++9Z3Sd8fHxJCQksGLFCkV+YD6INdpDq9VSUlLCq6++SmhoKADdu3fHzs6OyMhITp48SY8ePSx9KWazRltYo86HwdS4a9qj51H2MNqiPrabqJ8ksajGV199RWFhIX//+991XZiVO52uWLGCgQMHGjVrobCwkDVr1jBs2DC8vLzIy8sD/ti0LC8vDzs7O8V2eVu6PeCPb6D3b+jWq1cvIiMjSUtLU2RiYY22sEadD4Mpcde2R4+bm1u1vRKVxyrfN0pjjbYwtawQdU3GWFTj4sWL+Pj4VPmg79ChA1AxtcwYubm53L59m2+++YaIiAjd49ChQxQVFREREcHf/vY3i8VvaZZuD6DGaYO17cdS16zRFtao82EwNm5D9ujx9fXl6tWrlJWV6R2vHFvRpk0bS4VvUdZoC1PKCqEE0mNRjSZNmpCenk5hYSFOTk6645VLgDdp0sSo+ho3bsyCBQuqHN+2bRunT59mzpw5uLu7mxe0FVm6PQD69evHnj17+Pnnn/W2qE9MTASgU6dOZkZtHdZoC2vU+TAYE7ehe/T07duXPXv2cOTIEfr37687vm/fPry8vBS78q012sLYskIohSyQVQ1XV1d++OEHkpOTcXJy4vbt2xw6dIioqChatGjB//t//083aj0xMZHffvuN9PR0jh07hoeHByqViitXrtC0aVPs7OywtbWlWbNmVR4nTpzg6tWrvPPOO9VOLVMKS7cHQIsWLUhLS+OHH34AKm4LHT58mOjoaAIDA3XjLpTGGm1hTJ1KYmjcsbGxfP311wQGBvL888+TlZWl97h3hlTLli1JSUlhz549uLm5UVBQoFti/6233sLPz68Or7hm1mgLY8oKoSSypHcNkpOT2bZtG5cuXSI/Px9vb2969+5NaGioXu+CofuUVGfZsmUcOXLkkVjS2xrtUVxcjFqt5uDBg9y6dQsvLy+eeeYZIiIisLe3t/o1mcoabWFonUpjSNyG7tFTqbCwkE2bNukt6R0aGqrYFUgrWbotjG03IZRCEgshhBBCWIwM3hRCCCGExUhiIYQQQgiLkcRCCCGEEBYjiYUQQgghLEYSCyGEEEJYjCQWQgghhLAYSSyEEEIIYTGypHcDM3z4cL2fVSoVzs7OtGnThuDgYAYPHqy3T0d0dDRqtZrJkyczaNCghxqrOecuKipi9+7dJCQkcOXKFfLy8mjUqBGtWrWiR48eDB48mKZNm5oVX+UCRrUthFYf3LtY04IFC6rdFO3UqVPMnDmTwMBA5s6d+7BDrNXw4cNp2rQp69evr+tQhKjXJLFooIKDgwEoLy/nxo0bpKSk8Ouvv5KcnMwHH3xQx9GZ5+zZsyxcuJCcnBwaNWpEp06d8PT0pKCggNTUVGJiYti+fTsff/yxIndQVbqoqCgWLVpU12EIIRRKEosGaurUqXo///LLL8ydO5dDhw4xYMAA+vTpA8DQoUPp378/Xl5edRGm0X777TdmzZpFSUkJL7/8MuHh4Xo7TpaXl/PTTz+xceNGsrKy6jDSR5ODgwNnzpzh5MmTBAQE1HU4QggFkjEWAoCePXsycOBAAH766SfdcQ8PD3x8fHBxcamr0Aym1WpZunQpJSUljBkzhtdff73KNtY2Njb069ePZcuW6ba0FoYbMmQIUHGbSgghqiM9FkKncufIe7/JVzfOITExkblz59KiRQtWrFiht020Vqtl9uzZJCcnM27cOF566SW9c5w5c4ZvvvmGlJQU8vPz8fLyok+fPoSHh+Ph4WFW/ElJSVy6dInHHnuM0aNHP7Csi4tLlWSpqKiIb775hsOHD3Pjxg3s7Oxo27YtQ4YMMXgDrIyMDN58802eeOIJFi5cWOX5msaNVG5YtmPHDnbu3Mn333/PjRs38PT0ZMiQIbz00kuoVCouXLhAVFQUZ8+epaysDH9/fyZMmFBlvMiyZcvYv38/CxYsQKVSoVarSU1NBaBbt26MGzeO1q1bG3RN9/rTn/7EiRMn+PXXXzlx4oRBt5JqGytz77VXqhyvERwczLhx4/j66685fvw4RUVFtG3blnHjxtGlSxcAdu3axffff8+1a9dwd3dn8ODBhIWFYWNT/fem0tJS/vWvfxEXF0d2drZu87vRo0fj4OBQbfldu3Zx4MABfv/9d8rLy2ndujX/9V//xXPPPac3Jgn+GMuxZs0atm3bxsGDB8nIyKBXr17Mnj271vYS4lEnPRZCp7CwEKDWnUWDgoIYOnQo169fZ+3atXrPxcbGkpycjL+/PyNHjtR77rvvvmPGjBkkJCTQokULnnzySRwcHPi///s/3n//fXJycsyKPzExEaj48DN2m/GCggJmzJhBVFQUd+7coXfv3nTp0oXz58+zePFiIiMjzYrNUJGRkfzzn//Ew8ODgIAANBoNGzduJDo6ml9//ZXp06eTkZGBv78/np6eHDt2jNmzZ1NcXFxtfQkJCcyaNQuNRkPPnj3x8vIiMTGR6dOnc+vWLZNiDA8PBx5Or0V+fj4ffPABP//8M507d6ZNmzakpKTw0UcfkZ6eztq1a1m3bh2urq4EBASQn59PdHQ0mzdvrrY+rVbLokWL2L59Oz4+PgQFBZGXl0dMTAyffvopZWVleuWLior46KOPiIyMJDMzky5dutC9e3euX7/Ol19+yT/+8Y9qz1NeXs78+fPZvn277r3+qNxOFMJc0mMhgIo/uMePHwfA19e31vLjxo0jOTmZH374gd69e9OvXz9+++03Nm3ahIuLC1OmTNH7xnj27FnWr1+Pt7c3s2fPpm3btrrzxsTEEBUVxdq1a5k+fbrJ15CWlgZAu3btjH7tpk2buHDhAj169GDmzJm6XpgrV64wc+ZMvvvuO3r27ElQUJDJ8RkiPj6epUuX0qZNG935J0+eTGxsLPv372fs2LGMGDECqPgmPWfOHJKTkzl8+HC1vQHfffcdf/3rXxkwYAAAZWVlfPHFFxw5coSdO3fyyiuvGB1jv379aNu2LSkpKSQlJREYGGjGFT/YsWPH6N+/P1OmTNH1JlT2gHz++ecUFBTotdfly5eZPHky3333HaGhoXq9aQA3b95Eq9WyatUqmjdvDsCdO3eYNWsWJ0+eZOfOnYSEhOjK//Of/+TMmTMMHDiQv/zlL7r67ty5w7x589i9ezd9+vShd+/eeufJysrC3t6eNWvW0KRJE6u1jxBKJD0WDVxZWRnXrl1jxYoVnD17Fnt7e4OmdjZq1Ihp06ZhZ2fHypUruX79OosXL+bu3bu8/fbbeHt765Xftm0b5eXlvP3227qkAiqmu4aFheHn58fRo0e5c+eOydei0WgAjL6lUlRUxN69e7GxsdH78ADw8fHR3Va5t6veWl555RXdh2Tl+YOCgiguLsbb21uXVEBFz1Llh+CpU6eqre/pp5/WJRUAtra2hIaGAhW3pUyhUqmIiIgArN9r4eLiwltvvaV3i2LkyJGoVCquXLlSpb1at25N7969KS4u5sKFC9XWGR4erksqoOL9Mm7cOAC+//573fHbt2/z73//m2bNmvHuu+/qvS88PDx4++23Adi9e3e153nttdckqRANkiQWDdTw4cMZPnw4I0eOZOLEiezbtw8nJyc++OADWrRoYVAdfn5+jB07Fo1Gw5QpU7hy5QoDBw6kf//+euXKy8tJTk7Gycmp2pkEKpWKrl27Ul5erut1MIVWqzXpdRcuXKCkpIQOHTrQsmXLKs9XDmpNSUkx+RyGqm7MQuWHYM+ePWt8rqbbGtW95vHHH3/gawzRt29f/Pz8OHfuHD///LPJ9dSmffv2uLq66h1zdnbGzc0NeHB71XRr7f73J0CvXr1wdXXl999/1yW3p0+f5u7duwQGBlZ7e7Bt27Y4OTnpxq7cS6VS6WZWCdHQyK2QBqpyHQsbGxvdAln9+vWr8ke8NiNHjuTw4cNcuHCBJk2aMGnSpCplNBqNbvzG/eMu7pebm2vU+e/l7u6u98FgqMoPoJoWuXJ1dcXFxYX8/HwKCgqsOkOmum+4lTNbHvRcaWlptfU99thjVY5VfvOu6TWGqOy1mD9/PtHR0fTq1cvkuh6kpm/8jo6O5ObmGt0mrq6uODs7V1tn06ZNycvLIycnBw8PDzIzM4GKwaG7du2qMcaSkpIqxzw8PGodqyREfSWJRQN1/zoWprp8+TLp6elARVKQmZlZZYxGeXk5UPGB9tRTTz2wvvtvoRjDz8+PlJQU0tLSdL0Mlnb/DABjVbaFteq3dn336tu3L+3ateP8+fMkJibSqFEjk+qprU0exJLXd39vVOVATj8/P4PGHd2rutklQjQUklgIk5WWlrJkyRJKS0t55plniIuLY8mSJSxdulTv25q7uzv29vbY2dlZLKGpTlBQEDt37iQ+Pp5x48YZPDOkcrR+RkZGtc/n5+eTn5+Po6NjlcGA97Ozq/iVKioqqvb5+rYo15gxY5g3bx7R0dG6cQr3e1CblJWVcfv2bavGeK+8vDwKCgqq7bW4efMmAI0bNwb+6O3p3r07b7755kOLUYhHnYyxECbbuHEjly5d4plnnuH9999nwIABXLp0iY0bN+qVs7W1pXv37mg0Gt1+E9bQq1cvWrduTVZWFv/6178eWLagoEDX09K+fXscHBxITU3l2rVrVcrGxcUB0LVr11q/Ibu7u2NnZ0dGRkaVqYulpaVWvf660KdPH9q3b09qaioJCQnVlqlM3H7//fcqzyUnJ3P37l2rxni/w4cPVzmWlJREXl4eLVu2xNPTEwB/f39sbGw4fvx4lX9LIUTNJLEQJvnll1/YsWMH3t7eunEVkyZNwtvbmx07dvDLL7/olQ8NDcXGxoZly5ZVOxshOzubnTt3mhWTSqXi/fffx8HBgejoaL766qsq35K1Wi3Hjh1j6tSpukF3jo6OPPfcc5SXl7N69Wq91/z+++/ExMQAMGzYsFpjsLe3p1OnTmg0Gr3ruXv3LuvWrauxV+RRNmbMGEB/RsW9nnjiCaAiQbv3+m/cuMH//M//WD/A+2zZskUvjjt37rBhwwbgj5VFoWJ8x7PPPsu1a9dYunRptWN3UlJSdOunCCEqyK0QYbTc3FyWL1+OSqVi6tSpusGMrq6uTJ06ldmzZ7N8+XK+/PJL3N3dgYoPl/HjxxMZGcn06dPx9fWlZcuWlJSUcPPmTa5cuYKTkxNDhw41KzY/Pz/mzZvHwoUL2bZtGzt27KBz5854enqSn5/PhQsXuH37Ng4ODnrjOV599VXOnTvHiRMnGD9+PN26daO4uJjk5GRKSkoYPnx4lbUKahIeHs4nn3xCZGQkhw8fpnHjxly4cIHi4mKCg4PZv3+/WdeoNL1796Zjx46cP3++2uebN2+uu+7JkyfTrVs3ioqKOHfuHEFBQZSWluoGSlqbt7c3vr6+vP322wQEBGBra0tycjL5+fn4+/tXSR4nTJhARkYGhw4d4vjx4/j5+eHl5cWtW7e4fv062dnZhISEWH19EyEeJZJYCKOtXLmSnJwcXn755SrbZ3fv3p2RI0eyfft2Vq5cycyZM3XPDRs2jM6dO/Ptt99y+vRpEhIScHJyokmTJrzwwgv86U9/skh8Xbt2Ze3atbpt0y9dukReXh6Ojo60atWKF154gcGDB+vNmHB2dmbhwoXExsZy+PBhEhISsLOzo3379gwZMkRvLYja9OjRg9mzZ6NWq0lLS8PR0ZGAgABef/119u3bZ5FrVJqIiIgHbpX+zjvv4OXlRVxcHElJSXh7exMaGsqoUaOYMGHCQ4tTpVIxY8YM1Go1Bw8eJCcnBy8vL4YOHcro0aOrjMtxdHTk008/Zd++fRw4cIBLly5x7tw5PD09ad68OSEhIQYv9y5EQ6HSWntivhBCCCEaDBljIYQQQgiLkcRCCCGEEBYjiYUQQgghLEYSCyGEEEJYjCQWQgghhLAYSSyEEEIIYTGSWAghhBDCYiSxEEIIIYTFSGIhhBBCCIuRxEIIIYQQFiOJhRBCCCEsRhILIYQQQliMJBZCCCGEsJj/D7fIEo2fAK18AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "%matplotlib inline\n", "tpf.plot();" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.6" } }, "nbformat": 4, "nbformat_minor": 2 }