{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Quickstart" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If you have a working version of Python installed on your system, it is easy to install Lightkurve and its dependencies using the ``pip`` package manager. In a terminal window or Jupyter notebook cell, type:\n", "```\n", "! python -m pip install lightkurve --upgrade\n", "```\n", "\n", "See our [installation instructions](about/install.rst) page for details and troubleshooting information.\n", "\n", "With Lightkurve installed, it is easy to extract brightness time series data (astronomers call this a *light curve*)\n", "from the tiny images of stars collected by NASA's Kepler and TESS planet-hunting telescopes.\n", "\n", "For example, let's download and display the pixels of a famous star named [KIC 8462852](https://en.wikipedia.org/wiki/KIC_8462852), also known as *Tabby's Star* or *Boyajian's Star*, which is known to show unusual light fluctuations.\n", "\n", "First, we start Python and use the `search_targetpixelfile` function to obtain the Kepler pixel data for the star from the [data archive](https://archive.stsci.edu/kepler/):" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:10:22.272397Z", "iopub.status.busy": "2023-11-03T14:10:22.272256Z", "iopub.status.idle": "2023-11-03T14:10:24.822998Z", "shell.execute_reply": "2023-11-03T14:10:24.822597Z" } }, "outputs": [], "source": [ "from lightkurve import search_targetpixelfile\n", "pixelfile = search_targetpixelfile(\"KIC 8462852\", quarter=16).download();" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Next, let's display the first image in this data set:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:10:24.825049Z", "iopub.status.busy": "2023-11-03T14:10:24.824920Z", "iopub.status.idle": "2023-11-03T14:10:25.094527Z", "shell.execute_reply": "2023-11-03T14:10:25.094188Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiEAAAGZCAYAAABfZuECAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8pXeV/AAAACXBIWXMAAA9hAAAPYQGoP6dpAACDeklEQVR4nO3deVxU9f4/8NcMMKyyo4KCICoQuKAm6rgkSrmSWaaWlphW11yStGtmpqlp6nX5drVV83JRMxfMLUxFUMHELcEFAUFEEEF2GNaZ8/uD35zLOAPMDHPOMOP7eR88bpzzPufzOeAw7/msAoZhGBBCCCGE8Eyo7woQQggh5PlESQghhBBC9IKSEEIIIYToBSUhhBBCCNELSkIIIYQQoheUhBBCCCFELygJIYQQQoheUBJCCCGEEL2gJIQQQgghekFJCCHE6GzZsgUTJkzAli1b9F0VQkgzTNUNnDBhgtaFLFy4EKNGjdL6+raqoqICR48eBQCEhobCxsZGq/s8efIEs2fPBqD6Z7VlyxbExMQoHBMKhbCwsICVlRU6duyIrl27om/fvggMDIRQyH1uWVNTgz/++AOXLl1CdnY2JBIJzM3N4erqir59+2LChAlwcHBQ+3719fVYtGgRHjx4AAAIDg7GokWLWrzu/v37OHXqFJKSklBYWAiBQAAHBwd4e3sjKCgIw4cPV7omMzMTiYmJuH37NrKyslBWVgYzMzO4uLigV69eGD9+PDp16sTZ83/22We4detWs/d3cnLC7t27lY5LpVLcunULV69eRUpKCnJyciCRSGBhYQF3d3e8+OKLGDt2rNb/FjWVmZmJuLg4JCUloaCgABUVFRCJRHByckL37t0xaNAg9O/fH6amav+pIW2cNq+51lx/5swZbNu2Te36ffzxxxg5cqTCsezsbKSkpCA9PR33799HZmYmamtrAQDHjh1T+96N5eXlYd68eaipqQFgvO9zXFP7L4O9vb3K49XV1aiurm42xtzcXOOKGYLKykrs27cPADBy5EjO//ALhULY2tqy39fU1ODp06d4+vQpbt26haNHj8LZ2RmzZ8+GWCzmrB75+fn44osvkJubyx6ztrZGVVUV7t+/j/v37+PkyZNYvnw5AgIC1Lrnb7/9xiYg6mAYBrt378aRI0cgk8kAAJaWlpBKpcjNzUVubi4ePXqk9ActNjYW//rXvxSOyev+8OFDPHz4ENHR0Xj//fcxZswYTp/fwsICFhYWKs/Z2dmpPL5jxw78+eef7PdCoRCWlpaorKxESkoKUlJScOzYMXz++efw9fVtsuzWkkgk2LFjB86fPw/59lMCgQBWVlaora1FdnY2srOzERMTA1dXV3zyySfw8fHhrD6Ee9q+5lp7vbm5eZPvLXISiYRNKrp37650fseOHS0m/ppgGAbffvstm4AQ7amdhPz3v/9VeXzv3r3sG3FTMUQ3nJ2dsXPnToVjdXV1ePDgAa5evYqTJ0/i6dOnWL9+PSZPnox33nmHk3ps3rwZubm5MDU1xezZsxEcHAxLS0vU1dXh+vXr+O6771BYWIhvvvkGP//8c4tJ6IMHD3DgwAF07NgR1dXVKCkpabEOP/74I44fPw5LS0tMnToVI0aMYFseSktL2VaOZ9XX18PMzAxisRjDhw/HCy+8ACsrK9TV1eHOnTv46aefkJWVhe+++w6urq7o06cPZ8//2muv4a233mrxWZ+tv729PUaMGIHBgwejW7duMDU1RVVVFS5cuIDdu3ejpKQEX331Fb777rsmk5nWqKiowD//+U88fPgQAoEAw4YNw5gxY+Dj4wMzMzMAQGFhIa5du4Zjx47hwYMHSElJoSTEwGn7mmvt9UOHDsXQoUObrdu8efOQlZUFHx8feHh4KJ03MTGBu7s7vL294e3tjcLCQhw5ckSDp1cUHR2NpKQk+Pn54e7du1rfh2iQhJC2yczMDN27d0f37t0xbtw4fPPNN0hKSsKBAwfg4eGBl156Safl5efn4/bt2wCAyZMnY9y4cQp1CQoKgoWFBZYvX46SkhLcvn0bffv2bfJ+UqkU27ZtQ319PebOnYt///vfLdbh2rVrOH78OExNTfHVV18pfeK3s7PD4MGDMXjwYKVrfX198dNPP8HJyUnhuJmZGXr37o3169dj7ty5KC4uxsGDB5WSEF0/v6bGjh2LuXPnKiU2lpaWePnll+Hh4YElS5agvLwc0dHRmDJlis7Kltu0aRMePnwIExMTLF68GEOGDFGKcXJywssvv4yQkBD88ccfEAgEOq8H4U9rXnO6uL459+7dY5OXl19+WWXMqlWrYGJiwn5/5swZjcuRKygowO7du9GuXTvMnj0bn3zyidb3IhwnIfX19bh27RquXLmC+/fvo7CwEOXl5bC2tkbXrl0xcuRIDBs2TOUfqOTkZCxbtgxAQ5/d/fv3ERUVhVu3bqGkpAR+fn5Yt24dG//gwQPs378ft27dQmVlJRwdHfHiiy9iypQpyM7OVriXKhKJBCdOnMDly5eRk5OD6upq2Nvbw8/PD6GhoUovmmf79eVjOuQCAgIU6scHW1tbfPbZZ5g3bx4KCwsRGRmJIUOG6LQ/vqioiP3vbt26qYxp3BxaVVXV7P2OHDmC9PR0BAcHIzAwUK06yFvexo4dq3GXQ+fOnZs9b2Njg0GDBuHkyZNIS0tTOq/r59dUS60Jvr6+cHd3R3Z2tsr6t9bVq1dx7do1AMDUqVNVJiCNCQQCjB07lm1+l8vLy8PFixeRnJyMvLw8dmyAi4sLAgMD8eqrr6J9+/bN3js2NhYnT57EgwcPIBAI0LlzZ4SEhOCVV15R61mysrJw7NgxdmyCUChEx44d8eKLL+LVV19V2Yokb/mVv75v3ryJI0eOIDU1FVVVVejQoQOGDRuG119/HSKRqMmyy8rKcOLECVy9ehW5ubmoqamBg4MD3NzcMGjQIAwfPhzW1tY6qbMutOY1p4vrm3P69GkADYl4Uy0mjROQ1tq+fTskEgkWLlzI2c/7ecJpEnL37l2sWbOG/d7KygpmZmYoLS3FjRs3cOPGDVy6dAmffvpps4Mp4+PjsWnTJtTX18PKykrpH9SlS5ewYcMG1NfXA2j4x1hcXIzjx48jISEBM2bMaLaeGRkZWL16NZ4+fQqgoZ/d3NwcT58+xYULF3Dx4kXMmDEDkydPZq9p164dbG1tUVZWBqAhAWj8DO3atVPzp6RbNjY2CA0NxS+//IInT57g9u3b6N27t0LMe++9h/z8fK0SpY4dO7L/nZ6ejgEDBijFyN/8hEIhvL29m7xXTk4O9u7dCzs7O7z33ntqlZ+bm4t79+4BAEaMGKFJ1dUmf/N49o0T0O3zc0Vef6lUqvN7nzhxAkDDGJiJEyeqfd2zr+9t27axSbypqSk7rkU+luTs2bP44osv4O/vr3QvhmHwf//3f+ynWYFAAGtra6SnpyM1NRXJycktJt6HDh1CREQE+zs2NzdHfX09Hjx4gAcPHuDs2bNYsWJFs7+/w4cPs4OHraysUF9fj0ePHmHv3r24desWvvrqK5VvftevX8fGjRtRUVEBoOEN0srKCkVFRcjPz8fff/8NBwcHDBo0SGd1bs1rvrWvOS5fs9XV1Th//jyAhm4bS0tLnd7/WTExMbh27Rp69eqFUaNG4cmTJ5yW9zzgNAkxNzfH6NGjIRaL0aNHD1hZWQEAysvLce7cOezZswfx8fE4fvw4QkNDm7zPtm3b0KdPH8yaNQvu7u4AwA4KzMvLw7/+9S/U19fD29sb8+bNQ7du3cAwDG7evIlvv/1WaRxFY0VFRfjyyy9RUlKCQYMG4c0334SnpydMTU1RUlKC48eP4+DBg4iIiEDnzp3ZPwzLli1TmNWyefNmdOjQQSc/t9Z68cUX8csvvwAAbt26pZSEtIa9vT0GDRqES5cu4cCBA7Czs2PHRMhbvr777jsAwMSJExXetBuTv5HU1tZi/vz5CgNum3Pnzh0ADW9cXl5eSExMxNGjR5Geno66ujo4Ozujb9++mDhxota/j+TkZABAly5dlM7p6vmBhk/yZ8+eRVFRETuzJjAwEOPGjVPqLlJXaWkp2zTt6emp1T2aIpVK2a6oPn36NDmoVh1du3bFkCFDEBgYiI4dO0IoFEIqleL+/fvYs2cPrl+/jg0bNuDHH39U6no6duwYm4CMHz8eU6dOhZ2dHSorK3H06FHs27eP/Vujyp9//ondu3fD0tISb7zxBkJCQuDg4ACpVIrMzEz88ssvSEpKwpo1a7Bjxw6Vb2yZmZm4c+cO3njjDbYFQiKRICoqCr/++iuSkpIQExODkJAQhevu37+PtWvXora2Fh4eHggLC0OfPn1gamrKlh8XF6dUf13UWVutfc1x+Zq9ePEi29rYVFeMrhQXF+Pnn3+GSCTCRx99xGlZzxNOk5AePXqgR48eSsfbtWuH0NBQODk5Yf369S0mIe7u7li+fLnCpwo3NzcADbMqampqYG9vj9WrV7MtEAKBAH369MGqVauwYMGCJu8dGRmJkpISDB8+HIsXL1Y4Z29vj+nTp8PGxgY7d+7Evn37lD6dtEWdO3eGqakp6uvrkZeXp/P7L1iwADU1Nbh+/Tq+//57fP/99+zsEJlMBk9PT0ybNq3ZZvHjx4/jzp076Nu3r0bjVnJycgA0tPhERETg8OHDAMD+0ZaPsj979iz++c9/ol+/fho924ULF3D//n0AaLL+unh+AHj8+DFMTU1hYWGByspKpKenIz09HSdOnMDHH3+s1b+1PXv2oL6+HiYmJkrTFFsrPz+f/YPftWvXVt1rzpw5SsdMTEzQo0cPrFixAh9//DEePHiAhIQEhU/PtbW1+PXXXwE0fKr+4IMP2HPW1taYNm0aamtrcfDgQZXlSiQS7Nq1CwCwdOlShfE6JiYm6NatG7766issXrwY6enp+PPPP/Hqq68q3aeyshLTpk1TGFhsZWWFt99+G1lZWbh06RLOnz+vlIT8+OOPqK2thZubGzZs2KDQ5SIv/9luPl3VWVutfc1x+ZqVzxTr0qUL5wOfv//+e5SXl+Odd95h339I6+l1sbL+/fsDaPhjXFxc3GTcpEmTVDZrMgyDhIQEAMCYMWNUdoF07ty5yX7r2tpaxMXFAQBef/31JssPDg4G0PDpp7l6thUCgYCdLlxeXq50fufOnTh27JjWY1ZsbGzw2WefYeLEiex4nsrKSraZuKqqCmVlZU12Bzx58gQREREwNzfHP/7xD43Kljdhl5WV4fDhw+jVqxd27NiB/fv348CBA/jqq6/g4uKCqqoqfPPNN8jPz1f73jk5Odi+fTsA4IUXXmjyTby1z9+zZ08sXLgQu3fvxuHDh7Fv3z7s27cPCxcuhL29PSQSCTZs2ICUlBS16w40JFB//PEHgIbXTEvjXzTV+N8Sl92NJiYm7BuR/FO03I0bN9h6TJ06VeX1b7zxRpPjMRISElBZWcmuq9NU+cOGDQPQ0HWiipmZGV577TWV5wYOHAgASlPOc3Nz2ed55513VI754KrOrXnNt/Y1x9VrNjs7m52Z8myyp2sXL15EQkICPD09MWnSJE7Let5wPjtGIpEgOjoaiYmJePToESorK9mxG409ffq0ycWd/Pz8VB7Py8tDZWUlADS7HkPPnj1x7tw5pePp6ens3PIVK1a0+CxAw8hoTRbhMkb379/HmjVrUFRUhLFjx2LMmDHo0KEDSkpKkJiYiL179yIiIgK3b9/GihUrlMYD/Pvf/0Z1dTVmzZrVbHeFKvI1KWQyGRwdHfHFF1+w3QJCoRCBgYFYunQpFi9ejKqqKhw5cgTvv/9+i/ctLi7GqlWr2EHNS5YsaXKcUmufX9W0XGtra4waNQr+/v5YtGgRKisrsXv3bqxfv16tn8vt27exdetWAECvXr3w9ttvq3WdPt2+fRt//vkn7t27h8LCQna9ocbk47Tk5ONtnJ2dm/w0am1tDW9vb5VTJ+XHHj161OxYMfnfhYKCApXnPTw8muzycHR0BKD8AUBetlAo1OjTvq7qrK3Wvua4es3KB6SamZlxNj4MaEiefvjhBwiFQsyfP1+ng1wJx0lITk4Oli9frvCHxNzcHNbW1uwnSPmaEKr+AMk1NQK5tLSU/W/5C1+VpvrXG890UGdtCgAGsTgNwzBscqbrT6wSiQQrV65ESUkJpk6dqvBm17FjR4SGhqJLly744osvcO3aNZw5c0ahr/bUqVP4+++/4e3t3WwXXFMa/+EfO3asynEJPXr0QK9evXDz5k3cuHGjxXuWlJTg888/x+PHj2Fvb481a9bA2dlZZWxrn78lrq6uGDduHH777TfcuXMHZWVlLY6XSUlJwapVq1BbWws/Pz+lrktdafxvSVULmyZ2796NQ4cOsd8LhULY2NiwA0rliyA++3qTv+ZbGjPT1O9P/pqvra1l37Sb09TrvbkxF/Kf/bMtYfJWVFtbW43G0+iqztpq7WuOi9dsfX09+8Fy4MCBao8p08ZPP/2EkpIShIaGqhxeQFqH0yRk27ZtePr0Kdq3b49Zs2ahV69eCn/IpFKpWiPs1fmDqs06BI1nPxw6dKjZKXWG5NGjR6irqwPQ8KamS7GxsWzC1lRzdO/evdG1a1fcv38fCQkJ7JtwZWUlfvnlFwiFQsyZM0flH9TGn5rk4w/Mzc3Z1oTGbz7yQcqquLu74+bNmy027coTkOzsbNjb22Pt2rXN3rc1z68u+RRGhmHw5MmTZv/ApqSk4Msvv0RVVRV8fX2xcuVKzmYItG/fHpaWlqiqqkJGRobW97lx4wabgIwdOxZjx45F586dFV7nkZGR2L9/P/vvQVfkr/mhQ4fi008/1em9W6LtWin6rDPQ+tecrl+zAJCYmMi+DrkckJqcnIzY2Fg4OjrijTfeUJpy3zjhq6urQ1VVFQQCQasGbT9vOEtCCgoK2GbEJUuWqJwbrm7rQ1Mat5AUFRU1ud9HYWGhyuONu1Xy8/N13oeuL1euXGH/u2fPnjq9d3Z2NoCGn31zMxDc3Nxw//59hSlsFRUVbAvN0qVLmy0nNjYWsbGxABqSWflASE1nfDT3h7+kpATLli1jE5A1a9aoXG2xsdY8v67dvXtXIQFZtWpVs3VqLRMTE/j7++Pq1av4+++/UV1drdUf2wsXLgAA+vbt2+SYoKbGXslf8029puWaOi9f/luTsUK6Iv97U1ZWptHPTp91Blr/mtPla1ZO3hXToUMHnc7+e5b89VtUVNTiCtQ7duzAjh07YG1tzQ6eJi3jbGBq4y6Ypuba//33360qo2PHjuzgLvm0SlWaOte9e3e2+TcxMVHj8hu/WHT9iU1bFRUV7IJsrq6ueOGFF3R6f/kzy/+QNkWeYOr6U7mfnx/7x1ueEKjy8OFDAGhyyl9xcbFSAqJqSu6z+Hh++ZoKAoGgyQW7+E5A5OQrxFZWVmq07HXjVkf5mIWmZtgwDIOkpCSV5+QLwT19+hSPHz9WGSORSJCenq7ynPz1cP/+fYXuWD7IP4jJZDJ2wTd16LPOQOtfc7p6zco9ffqUHXw7atQoWo3XwHGWhDT+g5iZmal0XiKRYP/+/a0qQyAQsMv8/vHHH+wo7MZyc3Nx8eJFlddbWFiwmyUdOnSoxU8az/aDN35G+Sd8fSovL8e6devYBHDGjBk6HxsgTygZhkF0dLTKmKysLHYWQOMWsA4dOuDYsWPNfsnfdIODg9ljjd+szM3N2VURT548qTIRkC9YBUDlYmKqumDUSUBa+/zy65qTl5fHLgjm6+urcjxU4wTEz8+PtwQEaJjRJl/Z9tdff0V8fHyL10RHR+PUqVPs9/IPDqr+LgANr+Wmppb36dOHnfnV1KfNQ4cONTl2QiwWw9raGvX19di5c2ezvw+ZTKbyb4q23Nzc2MXXIiIiIJFI1LpOn3UGWv+a08VrtrGzZ89CJpNBKBRyvmvtqFGjmv179fPPP7OxCxcuxLFjx6gVREOcJSHu7u5wcXEB0NCc3viTSUpKCpYtW6aTF4t8Ol5JSQm++OILdo0H+WJlK1asaHYDtXfeeQeOjo4oKyvDkiVLEBMTo/DHobS0FPHx8Vi7di02btyocK2NjQ3b33nmzBlOVqhsSX19PdLT07Fv3z7MnTuX/QQ5ZcqUJpcwfu+99zBhwgR89tlnGpcnFovZQcARERHYv38/u2qsfPXCFStWQCqVwtTUFOPHj9fyyZr29ttvw9raGkVFRVizZg376Uomk+Hvv//G+vXrwTAM7O3tlcYclZaWsgmIg4MDvv766xa7YBpr7fMfPHgQW7ZswdWrVxX+/UskEsTExODTTz9FRUUFTE1NMXPmTKXyG48B8fPzw8qVKzVOQD777DNMmDBB7VVqn7V48WK4u7tDKpXim2++waZNm3D79m2FWW9FRUU4e/YsPv74Y2zfvl0hKZBPM7127Rp+/fVX9k2poqICv/32G3744YcmB1Sbm5uz++HExMTgp59+Yn/+EokEv/76Kw4cONDk9FcbGxt2jZLz589j1apVuHfvHttSI5PJkJ2djaioKHz00UcKXZu68P7770MkEiE3Nxeffvoprl27xv7cpFIpUlNTsX37doVWYl3UuTWveaB1rzldXC/HMAzbFdO3b98mByA/q66uDqWlpexX40So8fHS0lKVKyUT7nA2JkQoFOLDDz/E119/jYcPH2LRokVsMlBTU8Nu8rV8+fJWlePm5obw8HBs3LgR6enp+Pjjj2FpaQmZTIaamho4OTnhvffew7Zt29gdPhtzdHTEmjVrsHbtWuTk5GDLli0QCoWwtrZGXV2dwj9WVTuqjhkzBpGRkTh+/DhOnToFOzs7CIVC+Pj46HwQ2dOnTxWm6NXW1qKqqkrhk5GLiwvmzJnD2aJqFhYW+Pzzz/HVV1+htLQUkZGRiIyMhKWlJaqrq9m6mJubY9GiRZws6uPk5IQvvvgCq1evxs2bNzF37lz29yV/s7O3t8cXX3yh1JJw8uRJttm3qqqK3VOoKZs3b2aTaaD1z19XV4eYmBjExMQAaOiuMTU1VVhnxNraGgsWLFDZlfbf//6XHRyXnZ2tsFjXs5ydnbFly5Zmn08btra22LRpE7799lvEx8cjLi4OcXFx7PLpz87icHd3V1h+PTg4GDExMbh9+zb27NmDvXv3wtraGhKJBDKZDP3794e3t3eTLaWhoaHIyMjAuXPncPToURw/fhxWVlbs9cOGDYOpqSn7M37WyJEjUVtbix9//BHXrl3DtWvXYGZmBgsLC1RVVSkkU7pu6u/atSs+//xzbNiwAVlZWVi5ciW7bH3jsp9dD0SfdQZa95rTxfVySUlJ7DgNTQakxsXFYdu2bSrPTZ8+XeH7n3/+uc2sfv084HR2zIABA7B+/Xrs378fd+/eZTdp6t27N15//XWdDQQVi8Vwc3NT2MDOyckJQUFBePPNN9mm8aY+Hbm7u+Pbb7/F2bNnkZCQgMzMTJSXl8PU1BSurq7o2rUrAgMDIRaLla6dPHkyLC0tce7cOeTk5KCwsBAMw7S4+ZY2ZDIZO9ZAPgLbyckJHTp0gLe3N/r164c+ffo0uw+PLvTo0QM7duzAH3/8gStXriAnJwcSiQQWFhbsQLHx48drvAaIJvz9/bFjxw5ERUXhypUrePr0KQQCATw9PZvdzKtxwiafBtocVZ+KWvP8YrEYDMMgJSUFjx8/Rnl5OSQSCaytreHu7o7AwEC88sorTa5F07g+LbUkNjXbSz5oszUrTFpZWeGf//wnMjMzce7cOSQnJ6OgoAAVFRUQiURwd3dH9+7dMWTIEPTt21ehW1C+k+rBgwdx/vx5PHnyBAzDoHv37hg5ciReeeWVZpu0hUIhwsPD0adPH5w8eRJZWVmQSqXw9vZGSEgIRo8eza6Z0pQxY8agb9++OHHiBP7++288efIElZWVsLKyQseOHeHr64ugoCD06tVL659RU/r27YsffvgBR48exdWrV5GXl4eamho4OjqyW0OoGmypzzoD2r/mdHU98L8VUu3t7VvstiGGQcC0lRGVHIqIiMCBAwfQq1cvrF27Vt/VIURvnj59irCwMAiFQuzYsaPJGWWEEMIHvS7bzofS0lK2D1HTfUQIMTbyMUMjR46kBIQQonecL9vOh6NHj6KmpgZisRgdOnSAiYkJ6urqcPPmTezcuRMlJSWws7PjfCQ1IW1dUlISzMzMmtx3hRBC+GQU3TE//fQTjh49CgDsoFKJRMLOVrG2tsbnn3+u84W7CCGEEKI9o2gJCQ4OhlAoxK1bt1BUVISysjKIRCJ06NABffv2RWhoaIt7TRBCCCGEX0bREkIIIYQQw2P0A1MJIYQQ0jZREkIIIYQQvaAkhBBCCCF6QUkIIYQQQvSCkhBCCCGE6AUlIYQQQgjRC0pCCCGEEKIXRrFYWVskk8lQVFQES0tLTrbWJoQQvjAMg6qqKjg6OnK+S3dzamtrUV9fr9W1pqamTe4sTfSHkhCOFBUVISwsTN/VIIQQnfnll1/g7Oysl7Jra2sxe9arKC7V7m3LwcEBP//8MyUibQwlIRyxtLQE0PCitbKy4qXMV+3f4aUcFs8tPAKhCa/lQcDvJz5GWsdreXw/H98tgrwvBi2T8lrc7yURvJUlkUgQFhbG/l3Th/r6ehSXmuKnb9JgZSnT6FpJlRBz/tkd9fX1lIS0MZSEcET+B9fKyoq3JAT1PHf78N3NZGLkz2fsvz8hz+Vp9j6lg/L4fT7e/q400ha6ls0t6mFuodkvV8rQ8Me2ipIQQgghBkMGBjJo1sqlaTzhDyUhhBBCDIbs//9Ps2tIW0VtVIQQQgjRC2oJIYQQYjCkDAOphoOONY0n/KEkhBBCiMFgtBgTwtCYkDaLkhBCCCEGQwoGUg2TCk3jCX8oCSGEEGIwaHaMcaEkhBBCiMGgMSHGhWbHEEIIIUQvqCWEEEKIwZBB83U/aJ2QtouSEEIIIQaDi4GpqampiImJQVJSEvLz89GuXTv4+PhgxowZ6NSpk0JsXV0d9uzZg3PnzqGiogKenp6YPn06AgMDle6rbiwX9zQU1B3DsfDwcMydO1fh68SJE/quFiGEGCQpo91Xcw4dOoSEhAT07t0bc+bMwejRo3H79m18/PHHyMrKUojdunUrjhw5guHDh2POnDkQCoVYtWoVbt++rXRfdWO5uKehoJYQjm3evFkvG00RQogxYqB590pL7SYTJ07E4sWLYWZmxh4bOnQo5s2bh4MHD+KTTz4B0NBicv78eYSFhWHSpEkAgODgYMybNw+7d+/Gxo0b2evVjeXinoaEWkIIIYQYDCkEWn01x8/PTyEBAQA3Nzd4eHggOzubPRYfHw+hUIjRo0ezx0QiEUJCQpCSkoKCggKNY7m4pyGhJIQQQgh5BsMwKCkpga2tLXssIyMDnTp1Umrd7tGjBwAgMzNT41gu7mlIKAkhhBBiMGSMdl+aio2NRWFhIYYOHcoeKyoqgoODg1Ks/FhhYaHGsVzc05DQmBBCCCEGQ53uFVXXaCI7Oxvff/89fH19ERwczB6vra1V6rYBGrpE5Oc1jeXinoaEkhBCCCEGg+skpLi4GF999RWsrKywdOlSmJiYsOdEIhHq6uqUrpG/+cuTAU1iubinIaEkhBBCiMGQMQLIGM2SEHXjKysrsXLlSlRWVmL9+vVwcnJSOO/o6Kiyy6O4uBgAFOLVjeXinoaExoQQQggxGFzMjgEaWhNWr16NnJwcrFixAh4eHkoxXl5eyMnJgUQiUTh+79499rymsVzc05BQEkIIIeS5JpVKsWHDBqSkpGDp0qXw9fVVGScWiyGTyRAdHc0eq6urw5kzZ+Dj4wMXFxeNY7m4pyGh7hhCCCEGQwYhpFpc05xdu3bh8uXLGDBgAMrLy3Hu3DmF8yNGjAAA+Pj4QCwWIyIiAqWlpXB1dUVMTAzy8/OxYMEChWvUjeXinoaEkhBCCCEGg4sxIRkZGQCAxMREJCYmKp2XJyFAw1YckZGRCnu3rFixAgEBAUrXqRvLxT0NhYBhGC1mUJOWSCQSTJkyBfv37+dt2fYQ4WReymEJNPtD0OriGo1S56dAfnsrmXrlUe+c4vn5BEJ+/70w2iwO0RoyTT+ft85p2QHeytLH37Om6hC+8S7MLTVbuL2mSojNS/z0Wn+iGrWEEEIIMRhSRtjihnSqriFtEyUhHAsPD4dQqPgCGDduHMaNG6enGukQz41oTH09r+VByG/Li8BUeREiTssT8VyeGb9/bmRV1byWx9Rquq0a0YYMQo03sGtpTAjRH0pCOEa76BJCCCGqURJCCCHEYPCxbDvhDyUhhBBCDIaUEWg8xkOq4Wwawh9KQgghhBgMBgKNx4Qw1BLSZlESQgghxGBIIYQUmg2Kp+6YtouSEEIIIQZDCiGkGs7MoySk7aJ5S4QQQgjRC6NqCamqqsLhw4eRmpqK1NRUVFRUYOHChRg1apRCXFZWFvbt24f09HQUFxfD3NwcHh4emDRpEgYMGKDy3unp6di3bx/u3LmD2tpadOzYEa+88gpCQ0P5eDRCCCGQrxOiWUuIjFpC2iyjSkLKysrw66+/wsXFBV5eXkhOTlYZV1BQgKqqKowcORKOjo6oqalBQkICVq9ejY8++gijR49WiL9+/TpWr14Nb29vTJkyBZaWlnj8+DEKCwv5eCxCCCH/X8PsGM2vIW2TUSUhjo6OiIiIgIODA9LS0hAeHq4yrn///ujfv7/CsXHjxmHRokU4cuSIQhIikUiwZcsWvPjii1i6dKnS6qeEEEL4QwNTjYtRvaOamZnBwcFBq2tNTEzg7OyMyspKheNxcXEoKSnBjBkzIBQKUV1dDZmMlmcmhBB9kDFCrb5I22RULSGaqq6uRk1NDSQSCS5fvoxr165h6NChCjF///03rKysUFhYiLVr1yInJwcWFhYYMWIEZs+eDZFIpKfaE0LI80emRUsIjQlpu57rJGTnzp2Ijo4GAAiFQgwaNAgffvihQkxubi6kUinWrFmDkJAQvPPOO0hOTsbx48dRWVmJJUuWNFuGRCJRqy5mZmYwM+N3QzFCCCFEn57rJCQ0NBRisRhFRUW4cOECZDIZ6urqFGLkrSVjxozBBx98AAAYPHgw6uvrER0djbfffhtubm5NlhEWFqZWXaZNm4a33npL+4chhJDngJQRQEgDU43Gc52EuLu7w93dHQAQHByML774AqtXr8a//vUvCAQN/2jl3S3Dhg1TuHb48OGIjo5GSkpKs0nIL7/8otYuutQKQgghLaMpusaFRus0IhaLkZaWhpycHPaYo6MjAMDe3l4h1s7ODgBQUVHR7D2trKzU+qIkhBBCWibfwE6zL0pC2ipKQhqpra0FoDiOo1u3bgCgtCZIUVERAMDW1pan2hFCCJFBoNUXaZueyySkpKRE6Vh9fT1iYmIgEonYLhoAGDJkCADg9OnTCvF//vknTExM0LNnT07rSgghhBgroxsTIp+1Im+5SExMZP97/PjxsLa2xvbt2yGRSBAQEABHR0eUlJQgNjYWjx49wnvvvQdLS0v2ft7e3ggJCcHp06chlUoREBCA5ORkxMfHY/LkyXByctLLcxJCyPNIygi1GJjKTV1I6xldEhIVFYX8/Hz2+0uXLuHSpUsAgJdeegnW1tYYOnQoTp8+jZMnT6K8vByWlpbo1q0bZs6ciaCgIKV7zp07Fy4uLjhz5gz++usvuLi4YPbs2Xj11Vd5ey5CCCENK6Zq2oQv5aQmRBeMLgnZuXNnizHDhg1Tmu3SHFNTU0ybNg3Tpk1rTdUIIYS0EsMIINOwZYNpIV7dzU+3bNmCmJiYJu+ze/dutnU8OTkZy5YtUxm3ceNG+Pr6st/X1dVhz549OHfuHCoqKuDp6Ynp06cjMDBQ6VpNYg2B0SUhhBBCjJcUQo2HmbbUEqLu5qdjxoxBnz59FI4xDIMdO3agffv2KrvnJ0yYgO7duyscc3V1Vfh+69atiI+PR2hoKNzc3HD27FmsWrUKa9euhb+/v9axhoCSEEIIIQajYS8YTa9p/ry6m5/6+voqtGAAwO3bt1FTU4OXXnpJ5TX+/v4Qi8VNlp2amorz588jLCwMkyZNAtCwbtW8efOwe/dubNy4UatYQ/Fczo4hhBBC5Fqz+WlcXBwEAgGGDx/eZIxEIoFUqro9Jj4+HkKhUGH3dpFIhJCQEKSkpKCgoECrWENBLSGEEEIMhhQCnXfHaKu+vh4XL16Er68vOnTooDJm27ZtqKqqglAohL+/P8LCwhS6ZzIyMtCpUyellbV79OgBAMjMzISLi4vGsYaCkhBCCCEGg4vuGG1dv34d5eXlKrtiTE1NMXjwYPTv3x+2trZ4+PAhoqKisHTpUmzYsAHe3t4AGha+VNUKIz/WeKFMTWINBSUhhBBCDEZbagmJi4uDqakpu6hlY35+fvDz82O/DwoKglgsxvz58xEREYFVq1YBaFipW9W2HfJ9y+QreWsaaygoCeFYeHg4hELFoTfjxo3DuHHj9FQjojZGxmtxQkd+F74TWFm2HKTTAnleOjvnMa/FCUz53f8pxORN/gozZWDxMn/FNaettIRUVVXh8uXLCAwMVHv7Djc3NwwcOBAJCQmQSqUwMTGBSCRS2r0d+F9CIU8w5P+tbqyhoCSEY5s3b1ZrF11CCCEtkzECjTekk7W0UIgW/vrrr2ZnxTTF2dkZ9fX1qKmpgZWVFRwdHVV2oxQXFwOAwrRfTWINBc2OIYQQQjQUGxsLS0tLDBgwQKPr8vLyIBKJYGFhAQDw8vJCTk6OwsapAHDv3j32vJwmsYaCkhBCCCEGoy3soltaWoqbN29i4MCBbDKhKuZZmZmZSExMRGBgINtNLxaLIZPJEB0dzcbV1dXhzJkz8PHxUZjtokmsoeC1O+bnn3+GjY0Npk6dymexhBBCjISUEQIadsdIGQYtDU9VZ/NTuQsXLkAqlTbbFbNhwwaIRCL4+vrC3t4eDx8+xKlTp2Bubo53332XjfPx8YFYLEZERARKS0vh6uqKmJgY5OfnY8GCBQr31CTWUPCahJw4cULlBnGEEEKIOmSMAAKNx4S0HKPO5qdysbGxsLe3R+/evZu8X1BQEOLi4vD7779DIpHAzs4OgwYNwrRp0+Dm5qYQGx4ejsjISIX9YFasWIGAgACl+2oSawh4TUKcnJwgk/E744AQQojxkEIIaNi9IkXLWYg6m5/Kbdq0qcWY0NBQhIaGqnU/kUiEWbNmYdasWTqNNQS8jgkZOHAgbt26pTSohhBCCFGHjBFo9UXaJl6TkLfeegsuLi5YtWoV7t+/z2fRhBBCCGljeO2OWbt2LczMzHD37l2Eh4fDwcEBLi4uTS6wsnbtWj6rRwghpI2TQajxmqkyNbpjiH7wmoQkJyez/80wDIqKilBUVKQyVsD36oqEEELaPBkj0Hh2DFd7x5DW432KLiGEEKItSkKMC69JSPv27fksjhBCiJGRabFOCBfLthPdoL1jCCGEGAwpBGA0HhNC2iq9JCHylePS0tJQVlaGoKAghIWFAQDu3r2LtLQ0jBgxAu3atdNH9XSKdtElhBBCVOM9CTly5Aj+85//QCptWEJXIBCgrKxMIWbnzp0wMzPDmDFj+K6eztEuuoQQojs0JsS48LpOyJUrV7Br1y44Ozvjs88+w3//+18wz/TV+fn5wdbWFpcvX+azaoQQQgyADELIGA2/aK/WNovXlpAjR47AwsICq1evRseOHZuM69q1K3JycnisGSGEEEPQML6D1gkxFrymh/fv34ePj0+zCQgA2Nraori4mKdaEUIIMRRSRqDVF2mbeG0Jqaurg6WlZYtxJSUlMDEx4aFGhBBCDAnz/7tYNLuGo8qQVuO1JaRDhw7IzMxsNqaurg4PHjxQ2uqYEEIIIcaF1yQkKCgI+fn5OHLkSJMxhw8fRllZGQYPHsxfxQghhBgE2kXXuPDaHfP6668jNjYWv/zyC+7du4dBgwYBaOh+uXTpEi5duoS4uDh06NCB1tEghBCiRAaBxhvYabq4GeEPr0mIjY0N1qxZg/Xr1yM+Ph4JCQkAgOvXr+P69etgGAbu7u74/PPPaW0NQgghSmSMAAINWzYYaglps3hfrKxTp07Ytm0bEhMTcePGDeTn50Mmk8HZ2Rl9+vTB4MGDaVAqIYQQlWSMEAIamGo09LJsu1AoxMCBAzFw4EB9FE8IIcRAUUuIcaFl5AghhBCiF3ppCcnIyMDJkydx+/ZtFBUVAQAcHR3h7++P0aNHo1u3bvqoFiGEkDaOBqYaF96TkH379mH//v2QyRQ3V87JyUFOTg5Onz6NKVOm4K233uK7apygXXQNl4DnsUl13fldG8e0opbX8gSPnvBanomjA6/lyYpLeC2PqeNxoEMbGlRB3THGhdckJCYmBvv27YOFhQXGjRuHYcOGoUOHDgCA/Px8nD9/HidOnMD+/fvRsWNHBAcH81k9TtAuuoQQojuUhBgXXpOQo0ePwsTEBF9//bVSl4unpyc8PT0xePBgLFmyBEePHjWKJIQQQojuMIxA86SihfiqqiocPnwYqampSE1NRUVFBRYuXIhRo0YpxCUnJ2PZsmUq77Fx40b4+voqHKurq8OePXtw7tw5VFRUwNPTE9OnT0dgYKBWcZrGGgJek5Ds7Gz06tWr2TEf3bp1Q69evXDr1i0ea0YIIcQQyBhBi0mFEkbQ7CyMsrIy/Prrr3BxcYGXlxeSk5Obvd2ECRPQvXt3hWOurq5KcVu3bkV8fDxCQ0Ph5uaGs2fPYtWqVVi7di38/f01jtM01hDwmoRYWVnBxsamxThra2vqwiCEEMILR0dHREREwMHBAWlpaQgPD2823t/fH2KxuNmY1NRUnD9/HmFhYZg0aRIAIDg4GPPmzcPu3buxceNGjeI0jTUUvE7R7du3L27duoWampomY2pqanD79m307duXx5oRQggxBDIItPpqjpmZGRwcNBvILJFIIJVKmzwfHx8PoVCI0aNHs8dEIhFCQkKQkpKCgoICjeI0jTUUvCYhM2fOhKmpKb7++mvk5uYqnX/8+DHWrVsHU1NTzJw5k8+qEUIIMQBtYQO7bdu2YcqUKZg0aRKWLVuGtLQ0pZiMjAx06tRJqVW/R48eAMDuKK9unKaxhoLT7pht27YpHfPy8kJiYiLmzp0LLy8vtG/fHkDD7JjMzEwwDIMXX3wR//3vf7FgwQIuq0cIIcTAaDsmRBdMTU0xePBg9O/fH7a2tnj48CGioqKwdOlSbNiwAd7e3mxsUVGRytYV+bHCwkKN4jSNNRScJiFnz55t8pxMJsP9+/dx//59pXOJiYkQCAQaJyHqjnDOysrCvn37kJ6ejuLiYpibm8PDwwOTJk3CgAEDmi1j//79iIyMhIeHB7Zv365R/QghhLSOjIEWSYhuyvbz84Ofnx/7fVBQEMRiMebPn4+IiAisWrWKPVdbWwszMzOle4hEIva8JnGaxmqrrq4Ojx49QmlpKSorK2FtbQ07Ozt07txZZdmtxWkSsnbtWi5vr0TdEc4FBQWoqqrCyJEj4ejoiJqaGiQkJGD16tX46KOPFPrbGnv69CkOHDgACwsLLh+DEEJIE/TZEqKKm5sbBg4ciISEBEilUnYDVpFIhLq6OqV4eaIgTxzUjdM0VhOlpaU4e/Ysrly5gtTUVNTX1yvFmJqaokePHnjxxRcxcuRI2NnZaVWW0n11cpcm9OzZk8vbK1F3hHP//v3Rv39/hWPjxo3DokWLcOTIkSaTkF27dsHHxwcymQxlZWU6rz8hhBDD4+zsjPr6etTU1LDjNRwdHVV2jxQXFwMAnJycNIrTNFYdubm52LNnDy5dusQmHra2tujUqRPatWsHS0tLSCQSVFRU4NGjR7h9+zZu376NyMhIDBo0CG+//Tbc3Fq30rNe9o7hijYjnOVMTEzg7OyscoARANy6dQvx8fHYtm0bfvjhh9ZUkxBCiJaYNtYSAgB5eXkQiUQKreReXl5ISkqCRCJRGEh679499rwmcZrGtuT777/HqVOnIJPJ0KtXLwwfPhwBAQHo2LFjs8+ZlJSEuLg4XLx4EQkJCRg9ejQ++OADtct9lt520ZVKpSguLkZ+fn6TX1yrrq5GaWkpHj9+jCNHjuDatWvo3bu3yrr+8MMPePnll+Hp6cl5vQghhKjGxRRddZWWliody8zMRGJiIgIDAxX2CROLxZDJZIiOjmaP1dXV4cyZM/Dx8YGLi4tGcZrGtuT06dMYO3Ysdu3ahdWrV2PUqFHNJiAA0LFjR7z88stYu3Ytdu3ahTFjxuD06dNql6kK7y0hN27cwIEDB5CSktLsHGsA+P333zmty86dO9lfplAoxKBBg/Dhhx8qxUVHR6OgoABr1qzRuAyJRKJWnJmZGSeDfgghxJhwNSbk+PHjqKysZLs7EhMT2f8eP348rK2tsWHDBohEIvj6+sLe3h4PHz7EqVOnYG5ujnfffVfhfj4+PhCLxYiIiEBpaSlcXV0RExOD/Px8hUkX6sZpGtuSn3/+WeueA6Ch6+f999/H5MmTtb4HwHMSEh8fjw0bNoBhGNja2sLFxQWWlpZ8VkFBaGgoxGIxioqKcOHCBchkMqVBP2VlZdizZw+mTJmi1UCcsLAwteKmTZtmNDsHE0IIV7jqjomKilJogb906RIuXboEAHjppZdgbW2NoKAgxMXF4ffff4dEIoGdnR0GDRqEadOmqRwbER4ejsjISIV9XlasWIGAgACt4jSNbU5rEhBd3ofXJGTfvn0AgPnz52PkyJFKW9zzzd3dHe7u7gAalr794osvsHr1avzrX/+CQNDwjzYyMhI2NjYYP368VmX88ssvai1BT60ghBDSMgaab2AnUKM7ZufOnS3GhIaGIjQ0VO1yRSIRZs2ahVmzZukkTtNYQ8BrEvL48WMEBAQgJCSEz2LVJhaLsX37duTk5KBz587Izc3FqVOnMHv2bBQVFbFxdXV1kEqlePLkCaysrNCuXbsm72llZUX74BBCCCEq8JqE2Nvbw9bWls8iNSKfay0fx1FYWAiZTIYff/wRP/74o1L87NmzERoaijlz5vBaT0IIeV4xjOYtIVzPjjEm165dw3fffYeff/6Zl/J4TUKGDBmCc+fOoba2VutFVXShpKQE9vb2Csfq6+sRExMDkUjEdtF4eHhg2bJlStdHRkaiqqoKc+bMUbl9MyGEEG7ItEhCBJSEqK26uprXjfB4TUKmTZuG5ORkrF69GnPnzuXkDVydEc7bt2+HRCJBQEAAHB0dUVJSgtjYWDx69AjvvfceO1hWPujoWUePHgUAlecIIYRwh2EavjS7iJOqGJTIyEi14h49esRxTRTxmoRYWFhgzZo1WLJkCf7xj3+gffv2cHJyanKAqjbLvqszwnno0KE4ffo0Tp48ifLyclhaWqJbt26YOXMmgoKCtHs4QgghnJNBAEbDdT/UGZhq7H777TdYW1u3OEZRF/vPaILXJKSwsBDLly9Hbm4uGIZBXl4e8vLyVMbKZ6doSp0RzsOGDcOwYcO0uj8ArFu3TutrCSGEaI/GhGinY8eO8Pf3x8KFC5uNky+lwRdek5Cff/4ZOTk56N27NyZMmIAOHTrodZ0QQggh5Hng6+uLu3fvqhXLaNzfpT1ek5CbN2/Czc0NK1euZHcaJIQQQtRFA1O1Exoaijt37rQYFxAQoNVQCG3xmoTIZDJ4e3tTAkIIIUQrNDBVO926dUO3bt1ajLOzs0PPnj15qFEDXpMQHx8fPH78mM8iCSGEGBEaE2JceF03fcaMGcjKysKpU6f4LJYQQoiRkCchmn4R1QoKCvD555/rrXxeW0IePnyIkSNHYseOHYiNjUWfPn2anaIbHBzMZ/UIIYS0cTQmRLdqampw69YtvZXPaxKydetWCAQCMAyD27dvNzlIhmEYCAQCo0hCwsPDlZKscePGYdy4cXqqEVGXUEe7TKpL4szvKsJmtzJ5LU/m1ZnX8kzyi3ktj6mv57U8QowBr0nI1KlTtV7/w1Bt3ryZNrAjhBAdoYGpxoXXJOStt97iszhCCCFGhgamGhdekxBCCCGkVWhMiFGhJIQQQojBYKB57wr1xrRdvCYhmk4D4nPVNkIIIW2fNt0xNEW37eI1CUlOTm4xRj575nkbwEoIIYToA597xTyL9w3sVJHJZHj69Clu3LiBY8eOYezYsRg7diyfVSOEEGIIqD9GpxwcHPCPf/xDb+XzmoS0b9++yXMdO3ZEQEAAevXqhS+//BI+Pj7NxhNCCHn+NEzR1bQ7hqPKGAFra2uMGTNGb+Xzumy7Onr37o1u3brh4MGD+q4KIYSQNka+ToimX6RtapOzY5ydnXHt2jV9V4MQQkgbw8XA1KqqKhw+fBipqalITU1FRUUFFi5ciFGjRinEpaamIiYmBklJScjPz0e7du3g4+ODGTNmoFOnTgqxycnJWLZsmcryNm7cCF9fX/b7uro67NmzB+fOnUNFRQU8PT0xffp0BAYGKl2rSawmHj58iJKSEvTo0QMWFhatupcm2lwSUlNTg7S0NJiZmem7KoQQQp4DZWVl+PXXX+Hi4gIvL68mJ1EcOnQId+/ehVgshqenJ0pKSnD8+HF8/PHH2LRpE7p06aJ0zYQJE9C9e3eFY66urgrfb926FfHx8QgNDYWbmxvOnj2LVatWYe3atfD399c6VhOHDx/GuXPnsHHjRvTo0YM9XlxcjNOnT4NhGAwYMABeXl5al6EKr0lIfn5+k+eqq6uRk5ODI0eO4OnTpxg2bBiPNSOEEGIYBFqsgNp8vKOjIyIiIuDg4IC0tDSEh4erjJs4cSIWL16s8CF56NChmDdvHg4ePIhPPvlE6Rp/f3+IxeImy05NTcX58+cRFhaGSZMmAWjYvHXevHnYvXs3Nm7cqFWsplJSUuDq6qqQgNTV1WHJkiUoKCgAwzDYu3cv3n33XbZsXeA1CZk9e3aLU28ZhkGnTp0QFhbGU60IIYQYCm3GeLQUb2ZmBgc1Nqz08/NTOubm5gYPDw9kZ2c3eZ1EIoG5uTlMTEyUzsXHx0MoFGL06NHsMZFIhJCQEERERKCgoAAuLi4ax2qquLhYqSXl/PnzyM/PR/fu3TF8+HCcPHkS//nPf+Dr64sXXnhBq3KexWsS4u/v32QSYmpqCgcHB/Ts2RPDhg2DSMTvjqJcoV10CSFEh9rYFF2GYVBSUgIPDw+V57dt24aqqioIhUL4+/sjLCxMoXsmIyMDnTp1UtroVN4ikZmZySYWmsRqqq6uDpaWlgrHEhISIBQK8emnn6JDhw4YNGgQ3n//fRw9etQwk5B169bxWVybQLvoEkKI7rS1FVNjY2NRWFiIt99+W+G4qakpBg8ejP79+8PW1hYPHz5EVFQUli5dig0bNsDb2xsAUFRUpLIVRn6ssLCQPaZJrKacnJzw5MkT9vvq6mrcvHkTvr6+6NChAwDAxcUF/v7+uHv3rtblPKvNTdElhBBCmsRo+cWB7OxsfP/99/D19UVwcLDCOT8/P3z22WcICQlBUFAQJk+ejE2bNgEAIiIi2Lja2lqVEzHkvQG1tbVaxWqqZ8+eSEtLQ2ZmJgDg3LlzqK2tRb9+/RTiHBwcUFZWpnU5z2pzs2MIIYSQtq64uBhfffUVrKyssHTpUpXjPZ7l5uaGgQMHIiEhAVKpFCYmJhCJRKirq1OKlScUjYcmaBKrqddeew1xcXH4/PPP4e/vj+vXr0MoFGLo0KEKcWVlZTpt3ec0Cdm3b1+rrp82bZqOakIIIcQYtIXumMrKSqxcuRKVlZVYv349nJyc1L7W2dkZ9fX1qKmpgZWVFRwdHVV2oxQXFwOAwr01idWUu7s7li1bhm+//RaXL1+GQCDA22+/jY4dO7IxMpkMaWlpWo87UYXzJES+IZ26Gg9cpSSEEEKIAj0PTK2trcXq1auRk5ODNWvWNDkgtSl5eXkQiUTsgmBeXl5ISkqCRCJRaGG4d+8ee15Ok1ht9OvXD7t27UJubi6sra2Vxp/cuHEDFRUVSq0jrcFpEvLOO+9oFF9YWIjTp0+jtraWdtElhBCiggAtrfuh+prWk0ql2LBhA1JSUrB8+XKFVU+fVVpaCjs7O4VjmZmZSExMRL9+/dhZk2KxGFFRUYiOjmbX36irq8OZM2fg4+Oj0OqgSay2hEIhOnfurPKcQCDAqFGjMHjw4FaXI8dpEvLGG2+oFVdcXIwDBw7gzz//RF1dHaysrBAaGspl1QghhBgijlpCjh8/jsrKSra7IzExkf3v8ePHw9raGrt27cLly5cxYMAAlJeX49y5cwr3GDFiBPvfGzZsgEgkgq+vL+zt7fHw4UOcOnUK5ubmePfdd9k4Hx8fiMViREREoLS0FK6uroiJiUF+fj4WLFigcH9NYrnQt29f9O3bV6f31OvA1NLSUhw8eBB//PEH6urqYGFhgddeew0TJ06EjY2NPqtGCCGkLeIoCYmKilJY1fvSpUu4dOkSAOCll16CtbU1MjIyADQkKImJiUr3aJyEBAUFIS4uDr///jskEgns7OwwaNAgTJs2DW5ubgrXhYeHIzIyUmE/mBUrViAgIECpDE1iDYFekpCysjIcOnQIf/zxB6qrq2FhYYHXX38dr732Gtq1a6ePKhFCCHmO7dy5s8UYTda6Cg0NVbtFXyQSYdasWZg1a5ZOY9VRUFCArVu3Yu3atTq5n6Z4TUIqKipw+PBhHD9+nE0+Jk2ahEmTJsHW1pbPqhBCCDFEjBZ7x3C4WJmhq6mpwa1bt/RWPi9JSGVlJaKionD8+HF2Df2JEyfi9ddfVxq4QwghhDSFi71jiP5wmoRIJBIcOXIER48ehUQigUgkwquvvorXX38d9vb2XBZNCCHEGLWxvWNI63CahMyaNQtVVVUwMzPDhAkT8MYbb6i1UyEhhBCiEnXHGBXOW0IEAgHq6upw8uRJnDx5UqPro6KiOKoZfyaPnYJn56hLs0wgzaIV8zUl5HnQctribryW53yT349rNX35fb66di0va61LNveb3lqdC4yM39/fadkB3sqSSCSYMmUKb+U1R8A0fGl6DWmbOH8nlK+WKpVKuS6qTapNMAfqKQsnhBBCnsVpEnL06FEub08IIeR5Q2NCjAr1CRBCCDEcNCZE5zTZ303XKAkhhBBiOKglRKccHBzwj3/8Q2/lUxJCCCHEcFASolPW1tYYM2aM3sqnJIQQQojhoCTEqFASQgghhBCUlpbiwIEDuHnzJsrKymBjY4OQkBBMnDiRszKNLgmpqqrC4cOHkZqaitTUVFRUVGDhwoUYNWqUQlxWVhb27duH9PR0FBcXw9zcHB4eHpg0aRIGDBigEJuamoqYmBgkJSUhPz8f7dq1g4+PD2bMmIFOnTrx+XiEEPKc02JgKmhgaksKCwuxePFiFBUVAQBsbW3x6NEjZGVlsTG3b99GVVUVevXqBZFIpJNyhTq5SxtSVlaGX3/9FdnZ2fDy8moyrqCgAFVVVRg5ciTef/99TJ06FQCwevVqREdHK8QeOnQICQkJ6N27N+bMmYPRo0fj9u3b+PjjjxV+QYQQQrglX6xM0y/SvIiICBQWFmLMmDHYv38//vvf/yrNmpFKpVi9ejXi4uJ0Vq7RtYQ4OjoiIiICDg4OSEtLQ3h4uMq4/v37o3///grHxo0bh0WLFuHIkSMYPXo0e3zixIlYvHgxzMzM2GNDhw7FvHnzcPDgQXzyySfcPAwhhBBFNCaEE9evX4ebmxs++OADCASqW4569eoFe3t7XLlyBSEhITopl9eWkISEBJSVlXFahpmZmdb705iYmMDZ2RmVlZUKx/38/BQSEABwc3ODh4cHsrP5XRqaEEII0bXKykp4eXk1mYDIeXl54cGDBzorl9eWkPXr10MgEMDDwwMBAQHsl52dHZ/VUFBdXY2amhpIJBJcvnwZ165dw9ChQ1u8jmEYlJSUwMPDg4daEkIIAWjvGK44OjqisLCwxbh27dqhtLRUZ+XymoS8+uqruHXrFjIyMpCVlcVuaNe5c2f07NkTPXv25D0p2blzJzsGRCgUYtCgQfjwww9bvC42NhaFhYV4++23mw80UfNfvwyAjAZPEUII4V+fPn1w5swZZGVloUuXLk3GVVRUoL6+Xmfl8pqEvPfeewAadmS8desWkpOTkZycjMzMTGRnZ+OPP/4AAHTq1Am9evVSKxlordDQUIjFYhQVFeHChQuQyWSoq6tr9prs7Gx8//338PX1RXBwcLOxFqNq1apHfaoJ6tPMWg4khJDnGS3bzolXX30VZ8+exTfffIOVK1eiffv2SjE1NTVIT0+Hk5OTzsrVy8BUKysrDBgwgJ0KK09Krl69irNnz+LRo0fIycnhJQlxd3eHu7s7ACA4OBhffPEFVq9ejX/9618q+8aKi4vx1VdfwcrKCkuXLoWJSfPbk1efEQFSNV4AMq2qTwghzxcamMoJd3d3vP/++/j++++xYMECpVVUq6ursX37dpSVlWHQoEE6K1evs2Nqa2tx9+5dJCcn49atW0hLS2NbIZydnfVSJ7FYjO3btyMnJwedO3dWOFdZWYmVK1eisrIS69evVy8blAqAesrCCSFEJygJ4cyYMWPg4OCA7du349ChQwCA8+fP4/bt2ygoKIBUKkW7du3w5ptv6qxMXpMQVUlHfX09GIaBi4sLhg4dyo4NUdUUxFcdgYbWmWePr169Gjk5OVizZg0NSCWEED2ggancGjhwIPr06YPo6Gj89ddfePDgAfLy8iASidCvXz/MnDlTp40EvCYh06ZNYwe0ODs76zXpKCkpgb29vcKx+vp6xMTEQCQSsV00QMMCLRs2bEBKSgqWL18OX19fXutKCCHk/+OgJUTdlbYBoK6uDnv27MG5c+dQUVEBT09PTJ8+HYGBgVrHcnHP1rCwsMDEiRPZ5drr6+thaspNusBrEiLvaunSpQtGjBiBnj17olu3bi3OS9bU8ePHUVlZyU43SkxMZP97/PjxsLa2xvbt2yGRSBAQEABHR0eUlJQgNjYWjx49wnvvvQdLS0v2frt27cLly5cxYMAAlJeX49y5cwrljRgxQqf1J4QQwh/5StsuLi7w8vJCcnJyk7Fbt25FfHw8QkND4ebmhrNnz2LVqlVYu3Yt/P39tYrl4p66xFUCAvCchMyePRtJSUm4c+cOdu/eDYFAAEtLS7zwwgtsi4i3t3erk5KoqCjk5+ez31+6dAmXLl0CALz00kuwtrbG0KFDcfr0aZw8eRLl5eWwtLREt27dMHPmTAQFBSncLyMjA0BDMpOYmKhUHiUhhBDCEw5aQtRdaTs1NRXnz59HWFgYJk2aBKBhQsO8efOwe/dubNy4UeNYLu6pjpam4vJ1H16TkNDQUISGhoJhGGRmZiIpKQm3bt3CnTt3cPXqVQgEAlhZWbFJibY79+3cubPFmGHDhmHYsGFq3W/dunVa1YMQQohucTEmRN2VtuPj4yEUChW29RCJRAgJCUFERAQKCgrg4uKiUSwX91THggULMGTIELzxxhvN7rPWlPv37+PAgQO4dOkSfv/9d42vl9PL7BiBQICuXbuia9eumDhxIpuUnD17FtHR0bh69SquXr3K6fbBhBBCDJAe1wnJyMhAp06dYGVlpXC8R48eAIDMzEw2CVA3lot7qmPq1KmIiorCxYsX0aVLFwwbNgwBAQHw9vZW2qYEaJickZGRgeTkZMTFxSE7Oxvm5uaYNm2aWuU1Ra9TdPPz8xUWLSsoKGB37eOyD4oQQoiB0uMU3aKiIpUtJvJjjZc9VzeWi3uqY9q0aRgzZgx+++03xMTEICIiAgKBAEKhEC4uLrC2toaVlRUkEgkqKirw9OlTyGQyMAwDKysrTJgwAZMnT271Cue8vtOrSjqAhn1YTE1N4efnh4CAAPTs2ZNmoBBCCFGizym6tbW1KlsJRCIRe17TWC7uqS57e3u8//77ePfdd3Hx4kVcuXIFd+7cQV5enlKsg4MDXnjhBbz44osYMmQIW2Zr8T4wVSAQgGEYmJmZwc/Pjx2Q6uvrq7OHIoQQQnRNJBKp3NZD/ubf+D1M3Vgu7qkpc3NzjBw5EiNHjgQAlJaWoqSkBBKJBFZWVrC3t+dsTzdekxB/f3+FpENVRkcIIYQ0SY/dMU3tNFtcXAwACqtoqxvLxT1by87OjreNZHlNQmiWCSGEkNbQZ3eMl5cXkpKS2BYCuXv37rHnNY3l4p6GRO+jPysqKgAANjY2eq4JN0SDawAojsyWZplAmqX3H32rCXgePPx0cgCv5aXN+I7X8jCD3+LeeziE1/LypjryWp70ma0XOCeT8lve80qPLSFisRhRUVGIjo5m1+moq6vDmTNn4OPjozAzRd1YLu5pSPTyTnj16lUcPXoUd+/eVejLeuGFFzBhwgT0799fH9XiRG2COW1gRwghusJREqLOSts+Pj4Qi8WIiIhAaWkpXF1dERMTg/z8fCxYsEDhfurGcnFPQ8J7EvLTTz/h+PHj7FRcKysrCAQCVFZW4saNG/j7778xYcIEzJ49m++qEUIIaeME0KI7Ro0YdVbaBoDw8HBERkYq7N2yYsUKBAQot9SqG8vFPQ0Fr0nIhQsXcOzYMdjZ2WHKlCkYMWIE+4uVSCQ4d+4c9u/fj2PHjsHHxwdDhw7ls3qEEEKeU+qstA00tNrPmjULs2bN0lksF/c0FEI+Cztx4gTMzMywfv16tnlLzsrKCuPGjcO6detgamqKkydP8lk1QgghhPCM1yTkwYMH6NWrFzp16tRkTKdOndCrVy9kZmbyWDNCCCEGgdHyi7RJvCYhdXV1sLCwaDHOwsJC5YIshBBCnm/yKbqafpHmPX78WO3Yv/76S2fl8pqEuLq64tatW6iurm4yprq6Grdu3YKrqyuPNSOEEGIQqCWEEwsXLsSZM2eajampqcG3336r0zW/eE1ChgwZgtLSUqxduxa5ublK5x8/fox169ahrKyMBqUSQghRRkkIJxiGYROM8vJypfOpqalYuHAhTp8+rdNGAl5nx7z22mu4fPkybt68iblz58Lb2xvt27cHABQUFCA9PR0ymQzdunXDxIkT+awaIYQQ8tzaunUrNm3ahEuXLrEJR58+fcAwDH777Tf8+uuvkEqlePnll3W6hAavSYi5uTm+/vprRERE4PTp00hLS0NaWhp7XiQSISQkBO+88w7Mzc35rBohhBBDoM0YD2oJaVGnTp2wadMm7NmzB4cOHcKXX36JMWPGICMjAykpKbCzs8P8+fMxYMAAnZbL+2JllpaW+OCDD/Duu+/i/v37KCoqAtCwMY+3t7daA1cJIYQ8p/S4bLuxMzExwTvvvIN+/fph9erV+OOPPwAAffr0QXh4OOzt7XVept42MLGwsIC/v7/KcyUlJThy5AhmzpzJb6UIIYS0aQItkhCaHaO+yspKnDx5EpJGey9lZWUhMzMTgYGBOi+P14GpLSkoKMAPP/yA2bNnIyoqSt/VIYQQ0tbQwFTOJCcnY/78+bhw4QK8vLzwf//3f3jjjTdQWlqKlStX4qefftL58hmct4TIZDKcP38eN27cQElJCezt7dGvXz8MGTIEQmFDDlRQUIB9+/bh3LlzkMlkAICBAwdyXTVeGPMuuoQQwjvqjuHE7t27ceTIETAMg0mTJmH69OkwNTWFp6cn+vXrh82bN+P48eNISkrCJ598Ak9PT52Uy+k7oVQqxcqVK5GUlMRuWAcAsbGxuHjxIpYtW4azZ8/ihx9+QE1NDRiGwcCBAzFt2jR4eXlxWTXe0C66hBBC2rrDhw/DyckJ4eHh6Nmzp8I5f39/fPvtt/j+++8RGxuLTz75BIcOHdJJuZwmISdOnMDNmzdhZmaGkSNHokuXLpBIJLh27RouX76M7du3488//wTDMAgMDMS7776Lrl27clklQgghBozGhHBjyJAhmDt3LmxsbFSet7KyQnh4OAYMGIAdO3borFxOk5Dz589DKBRi3bp16NGjB3t88uTJ2LFjB6KjoyEQCBAWFobXXnuNy6oQQggxBtQdw4lPP/1UrbghQ4bAz89PZ+VyOjD10aNH8PX1VUhA5CZNmgSgYW4yJSCEEELUQgNT9c7JyUln9+K0JaSqqgodOnRQeU5+3FjGfhBCCOEedccYF06TEIZh2BkwzxIIGgZrikQiLqtACCHEmFB3DCdu3bqlUXxAQIBOyqV5ooQQQshzbtmyZWzjgDp+//13nZTLeRISExODmJgYlecEAkGz53X1kIQQQowDdcdwY8SIESqTEIZh8PTpU9y/fx8SiQRBQUGwtrbWWbmcJyGN1wchhBBCWoW6YzixaNGiZs9XVFTg22+/RVZWFjZt2qSzcjlNQo4ePcrl7QkhhDxvKAnRCxsbGyxatAjvv/8+/vOf/+Cjjz7SyX3b1N4xhBBCSHMEWn6R1rOwsECPHj2QmJios3vSwFRCCCGGRcctG1u2bGlybCLQsK+Kk5MTkpOTsWzZMpUxGzduhK+vr8Kxuro67NmzB+fOnUNFRQU8PT0xffp0pd1o1Y1rC6qqqlBRUaGz+1ESQggh5Lk2ZswY9OnTR+EYwzDYsWMH2rdvr7Q414QJE9C9e3eFY66urkr33bp1K+Lj4xEaGgo3NzecPXsWq1atwtq1a+Hv769xnL4lJibi9u3bcHd319k9KQnhmDHvoiu00d0IaXWYl8l4La/Xprm8lpe0WHf7MajjekQvXsvrkHOV1/KY+npey4PQhN/ynlNczI7x9fVVasW4ffs2ampq8NJLLynF+/v7QywWN3vP1NRUnD9/HmFhYewK4cHBwZg3bx52796NjRs3ahTHtW3btjV5rqqqCrm5ucjKygLDMDpd5dzw3wnbONpFlxBCdIingalxcXEQCAQYPny4yvMSiQTm5uYwMVGdfMbHx0MoFGL06NHsMZFIhJCQEERERKCgoAAuLi5qx3Ht7NmzLca4uLhg2rRpCA4O1lm5lIQQQggxHDwkIfX19bh48SJ8fX1Vbj2ybds2VFVVQSgUwt/fH2FhYUrdMxkZGejUqROsrKwUjsv3UsvMzISLi4vacVxbu3Ztk+fMzMzg4ODQ5DYsrUFJCCGEEIPBx2Jl169fR3l5uVJXjKmpKQYPHoz+/fvD1tYWDx8+RFRUFJYuXYoNGzbA29ubjS0qKoKDg4PSveXHCgsLNYrjWs+ePXkp51mUhBBCCCGNxMXFwdTUFEOGDFE47ufnp7CNfVBQEMRiMebPn4+IiAisWrWKPVdbWwszMzOle8v3S6utrdUozlhREkIIIcRwcNwdU1VVhcuXLyMwMBC2trYtxru5uWHgwIFISEiAVCplx4iIRCLU1dUpxcuTCnmSoW6crs2ePVvrawUCAX766Sed1IPTJKS5edfq0OXgF0IIIYaP6+6Yv/76q8lZMU1xdnZGfX09ampq2LEdjo6OKrtSiouLAYCd9qtunK7l5+dzcl9NcZqEbN26VaNd+eQYhoFAIKAkhBBCiCKOW0JiY2NhaWmJAQMGqH1NXl4eRCIRLCws2GNeXl5ISkqCRCJRGHR679499rwmcbrWVrZV4TQJmTp1qlZJiLaqqqpw+PBhpKamIjU1FRUVFVi4cCFGjRqlEJeVlYV9+/YhPT0dxcXFMDc3h4eHByZNmqTyH54hrWZHCCHGjMuWkNLSUty8eRPDhg1TSCgan7ezs1M4lpmZicTERPTr1w9C4f92QhGLxYiKikJ0dDS7/kddXR3OnDkDHx8fdsaLunHGitMk5K233uLy9krKysrw66+/wsXFBV5eXkhOTlYZV1BQgKqqKowcORKOjo6oqalBQkICVq9ejY8++khhvjZgOKvZEUKI0eOwJeTChQuQSqVNdsVs2LABIpEIvr6+sLe3x8OHD3Hq1CmYm5vj3XffVYj18fGBWCxGREQESktL4erqipiYGOTn52PBggUaxxkroxqY6ujoiIiICDg4OCAtLQ3h4eEq4/r374/+/fsrHBs3bhwWLVqEI0eOKCQhbWU1O0IIIdyKjY2Fvb09evfurfJ8UFAQ4uLi8Pvvv0MikcDOzg6DBg3CtGnT4ObmphQfHh6OyMhIhVb0FStWICAgQKs4Xfr888/Rt29fvP7660rn8vPzYWlpiXbt2nFWvpzekpDy8nKkp6ejrKwM7du3V5j2pC35giraMDExgbOzM9LS0hSOt5XV7AghhIDTlpBNmzY1ez40NBShoaFqFysSiTBr1izMmjVLJ3G6lJycjPbt26s8N2fOHIwcOZKXlhjek5DS0lL8+OOPSEhIgEzWsBdIcHAwm4ScOnUKu3fvxvLlyznv6qiurkZNTQ0kEgkuX76Ma9euYejQoQoxbWU1O0IIIfwsVva8YxgGDMPPD43XJKS8vBxLlixBXl4evLy88MILL+DEiRMKMYMHD8Z3332HhIQEzpOQnTt3Ijo6GgAgFAoxaNAgfPjhhwoxrV7NzkTNX6QMgIz2mCGEkGbxtHcM4QevSchvv/2GvLw8TJ06lR20+mwS0q5dO3h6euLWrVuc1yc0NBRisRhFRUW4cOECZDKZ0qIxrV3NzmKUeqvd1aeaoD5NuRxCCCH/I2AYQMNP6QKePtUTzfGahPz1119wc3NrcdZMx44deUlC3N3d4e7uDqChS+iLL77A6tWr8a9//YudWtza1eyqz4gAqRotHPzuUk8IIYaJWkKMirDlEN0pLCxUa+EVgUAAiUTCQ40UicVipKWlIScnhz3m6OjIrlzXmNqr2UkFQL0aX9QVQwgh5DnDa0uIlZWVyjf0Z+Xl5SktCMMHeetG4wRIX6vZEUIIUUYDU3UnJiZG5fYqAoGgyXNyv//+u07qwGtLSPfu3ZGWloa8vLwmYzIzM5GRkaGTKbtNKSkpUTpWX1+PmJgYiEQitosGaGgdkclk7ABW4PlazY4QQtoURssvokQ+C0abL13htSVk/PjxuHbtGr7++mssWbJE4c0eAHJzc7F582YADYuHaeP48eOorKxkZ60kJiay/z1+/HhYW1tj+/btkEgkCAgIgKOjI0pKShAbG4tHjx7hvffeg6WlJXu/5301O0IIaUuoJUQ3nou9Y57Vr18/TJo0CYcPH8a8efPg6uoKgUCA69evY/78+cjOzoZMJsObb76p9fTcqKgohd0BL126hEuXLgEAXnrpJVhbW2Po0KE4ffo0Tp48ifLyclhaWqJbt26YOXMmgoKClO6pj9XsCCGENIGSCqPB+2JlM2fORLdu3fDbb7/hwYMHABoGeRYXF6Nz586YMmUKhg8frvX9d+7c2WLMsGHDMGzYMLXvqY/V7AghhCjTplWDWkLaLr0s2z5kyBAMGTIEpaWlePLkCRiGgbOzc8szTQghhBBiNHhNQioqKmBjY8N+b2dn1+QsmLS0NHTv3p2vqhFCCDEE2rRqUEtIm8Xr7Jh58+bh5s2bzcYwDIP9+/fj008/5alWhBBCDIWA0e6LtE28JiHFxcVYsWIFdu7cifr6eqXz+fn5+Oyzz7Bnzx5YW1vzWTVCCCGGQL5su6ZfpE3iNQn5+uuv4ezsjKNHjyI8PBxZWVnsuXPnzmHBggW4c+cOAgMD8e233/JZNUIIIQaAWkKMC69jQvz9/fHtt99ix44dOH/+PMLDw/H222/j/v37uHDhAszMzDBnzhxMmDCBz2oRQggxFDQmxKjwPjvGysoKixcvRv/+/bFjxw785z//AQB4e3sjPDxcaQEzQycaXANAcV8YaZYJpFl6mZikU9KyCl7Ls43L4Le89o68lhdgNpfX8rpE5/JaXr1Uymt5ENB+TIS0dXp5J5RKpXjw4AFqamrY5V/r6+shkxnfVrK1CeYNG9QRQghpNYEMz36ua/kaaglps3gdEwIAOTk5WLx4MaKiouDo6Ijly5dj2LBhyMrKQnh4uM42xSGEEGKEaO8Yo8JrS8gff/yBXbt2oaamBmKxGB999BFsbGwwYMAAvPjii/juu++wa9cuXLt2DR9//DEcHfltDieEENLGMRo3hFAS0obx2hLy3XffQSgU4uOPP8Y///lPhYXLhg8fjn//+9/w9/fH33//jfnz5/NZNUIIIYaApugaFV6TEF9fX2zbtg3BwcEqzzs7O2Pt2rV49913UVVVxWfVCCGEGACaomtceO2OWb9+PYTC5vMegUCA119/HX379uWpVoQQQp5nycnJWLZsmcpzGzduhK+vL/t9XV0d9uzZo7Cr+vTp0xEYGKh0rbqxmtzT2PCahLSUgDTm5eXFYU0IIYQYJA7XCZkwYYLSnmWurq4K32/duhXx8fEIDQ2Fm5sbzp49i1WrVmHt2rXw9/fXKlaTexobw1+sghBCyHNDq64VNa/x9/eHWCxu8nxqairOnz+PsLAwTJo0CQAQHByMefPmYffu3di4caPGsZrc0xhxmoTMnj0bAoEAq1evRseOHTF79my1rxUIBPjpp584rB0hhBCDw2gx51aDgakSiQTm5uYwMTFROhcfHw+hUIjRo0ezx0QiEUJCQhAREYGCggK4uLhoFKvJPY0Rp0lIfn4+gIbFyRp/TwghhGiDy5aQbdu2oaqqCkKhEP7+/ggLC1PonsnIyECnTp1gZWWlcF2PHj0AAJmZmWzCoG6sJvc0RpwmIUePHm32e0IIIUQjHCQhpqamGDx4MPr37w9bW1s8fPgQUVFRWLp0KTZs2ABvb28AQFFRERwcHJSulx8rLCxkj6kbq8k9jRGNCSGEEPJc8/Pzg5+fH/t9UFAQxGIx5s+fj4iICKxatQoAUFtbCzMzM6XrRSIRe15O3VhN7mmMeElCrl69ir/++gsFBQUwMzODp6cnRo0ahY4dO/JRPCGEECMhAHhZAdXNzQ0DBw5EQkICpFIpTExMIBKJUFdXpxQrTxTkiYP8v9WJ1eSexojzJGTTpk24cOECALCb1V25cgVRUVH49NNPERQUxHUV9MqYd9ElhBDeybTZDIbRamlOZ2dn1NfXo6amBlZWVnB0dFTZPVJcXAwAcHJyYo+pG6vJPY0Rp++Ef/75J86fPw8TExOMGDECXbt2RVVVFa5cuYKUlBRs2bIFO3fuhLW1NZfV0CvaRZcQQnSIx9VP8/LyIBKJYGFhAaBh/aqkpCRIJBKFgaT37t1jz8upG6vJPY0Rp8u2x8TEQCAQYOXKlViwYAHGjx+PyZMnY8OGDQgODkZVVRUuXbrEZRUIIYQYES6WbS8tLVU6lpmZicTERAQGBrILbYrFYshkMkRHR7NxdXV1OHPmDHx8fBRmsagbq8k9jRGnLSEPHjyAj48PevfurXTuzTffRExMDB48eMBlFQghhBgTbdYJaSF+w4YNEIlE8PX1hb29PR4+fIhTp07B3Nwc7777Lhvn4+MDsViMiIgIlJaWwtXVFTExMcjPz8eCBQsU7qlurCb3NEacJiFVVVVKS97KyQelSiQSLqtACCGENCsoKAhxcXH4/fffIZFIYGdnh0GDBmHatGlwc3NTiA0PD0dkZKTCPi8rVqxAQECA0n3VjdXknsaG0ySEYZgm94uRH2doi2VCCCFq4mJH3NDQUISGhqoVKxKJMGvWLMyaNUtnsZrc09jQFA1CCCGGgz63GhXOk5CYmBjExMSoPCcQCJo9//vvv3NZNUIIIQZGwMGYEKI/nCch1N1CCCFEZ2T6rgDRJV73jiGEEEJag1pCjAun64QQQgghhDSFBqYSQggxHNSoYVQoCSGEEGI4qDvGqFASQgghxGBwsU4I0R9KQjjmPctZacG2cePGYdy4cTovK0Q4Wef3bJZMymtx0qdPeS1PUFbGa3ld9lbyWp6sqJjX8sDwPK2B55l5QktzXst7blFLiFGhJIRjmzdvVtgZkRBCiPYENEXXqNDsGEIIIYToBbWEEEIIMRzUHWNUKAkhhBBiWDTNKQSc1ILoACUhhBBCDIaAAe+Djgl3KAkhhBBiOBhGiySEkpa2ipIQQgghhkMG6o4xIkaVhFRVVeHw4cNITU1FamoqKioqsHDhQowaNUohLjU1FTExMUhKSkJ+fj7atWsHHx8fzJgxA506dVK6b25uLiIjI3Hnzh2Ul5fDxcUFw4cPx2uvvQYLCwu+Ho8QQggxKkaVhJSVleHXX3+Fi4sLvLy8kJycrDLu0KFDuHv3LsRiMTw9PVFSUoLjx4/j448/xqZNm9ClSxc2tqCgAOHh4bC2tsa4cePQrl07pKSkYO/evbh//z6WL1/O1+MRQshzT0DdMUbFqJIQR0dHREREwMHBAWlpaQgPD1cZN3HiRCxevBhmZmbssaFDh2LevHk4ePAgPvnkE/b4uXPnUFlZiW+++YZNTkaPHg2GYRATE4OKigrY2Nhw+2CEEEIaUBJiVIwqCTEzM4ODg0OLcX5+fkrH3Nzc4OHhgezsbIXjEokEAGBvb69w3MHBAUKhEKamRvUjJISQto2SEKNCK6b+fwzDoKSkBLa2tgrHe/bsCQD49ttvkZGRgYKCAly4cAF//PEHxo8fT2NCCCGETzItv0ibRB/j/7/Y2FgUFhbi7bffVjjer18/TJ8+Hb/99hsuX77MHn/zzTcxY8aMFu8rb0lpiZmZmUL3ECGEEGU0JsS4UBICIDs7G99//z18fX0RHBysdL59+/YICAjA4MGD0a5dO1y9ehUHDhyAg4MDxo8f3+y9w8LC1KrDtGnT8NZbb2lVf0IIIdpTd8ZkcnIyli1bpvIeGzduhK+vr8Kxuro67NmzB+fOnUNFRQU8PT0xffp0BAYGahVnjJ77JKS4uBhfffUVrKyssHTpUpiYmCicP3/+PP7973/jhx9+gLOzMwBg8ODBkMlk2L17N4YNG6bUhdPYL7/8otYuutQKQgghauCgJUSTGZMAMGHCBHTv3l3hmKurq9J9t27divj4eISGhsLNzQ1nz57FqlWrsHbtWvj7+2scZ4ye6ySksrISK1euRGVlJdavXw8nJyelmJMnT8Lb25tNQOSCgoJw9uxZZGRkoE+fPk2WYWVlpVYSQgghRA0cJCGazJgEAH9/f4jF4mbvmZqaivPnzyMsLAyTJk0CAAQHB2PevHnYvXs3Nm7cqFGcsXpuB6bW1tZi9erVyMnJwYoVK+Dh4aEyrqSkBDKZ8qim+vp6AIBUKuW0noQQQhqRJyGafjXDz89PqTW6qRmTchKJpNm///Hx8RAKhRg9ejR7TCQSISQkBCkpKSgoKNAozlg9ly0hUqkUGzZsQEpKCpYvX67Uj9eYm5sbbty4gZycHIW+wfPnz0MoFMLT05OHGhNCCAHA27Lt8hmTqj6gbtu2DVVVVRAKhfD390dYWJhS90xGRgY6deqk1BLeo0cPAEBmZiZcXFzUjjNWRpeEHD9+HJWVlSgsLAQAJCYmsv89fvx4WFtbY9euXbh8+TIGDBiA8vJynDt3TuEeI0aMYP970qRJuHbtGpYuXcqumHrlyhVcu3YNL7/8ssouHEIIIdzga3aMqhmTpqamGDx4MPr37w9bW1s8fPgQUVFRWLp0KTZs2ABvb282tqioSOW6VfJj8vcldeOMldElIVFRUcjPz2e/v3TpEi5dugQAeOmll2BtbY2MjAwADQlKYmKi0j0aJyEBAQHYuHEj9u7di5MnT6K8vBwdOnTAjBkz8Prrr3P8NIQQQvjW1IxJPz8/hcUug4KCIBaLMX/+fERERGDVqlXsudraWpUTDkQiEXtekzhjZXRJyM6dO1uMWbdunUb37NGjB1auXKlljQghhOgMxy0hLc2YfJabmxsGDhyIhIQESKVSNl4kEqGurk4pXp5UyJMMdeOMldElIYQQQowYwwAyDZMQoXrx6syYVMXZ2Rn19fWoqalhx3Y4Ojqq7EopLi4GAPbe6sYZq+d2dgwhhBADxMHsGED9GZOq5OXlQSQSKWzj4eXlhZycHKVVs+/du8ee1yTOWFESQgghxHBwkIQ0njG5dOnSJmdMlpaWKh3LzMxEYmIiAgMDIRT+7y1VLBZDJpMhOjqaPVZXV4czZ87Ax8eHnfGibpyxou4YQgghzzV1Z0xu2LABIpEIvr6+sLe3x8OHD3Hq1CmYm5vj3XffVbjGx8cHYrEYERERKC0thaurK2JiYpCfn48FCxZoHGesKAkhhBBiOLQZmNpCvLozJoOCghAXF4fff/8dEokEdnZ2GDRoEKZNmwY3Nzel68LDwxEZGamwJ8yKFSsQEBCgVZwxoiSEEEKI4ZBpMTC1hdkx6s6YDA0NRWhoqNqlikQizJo1C7NmzdJJnDGiJIRj4eHhCv2EADBu3DiMGzdOTzUyYBpPy2tlcTU1vJZX/yiH1/J4x/PvD8Lmp1bqmsDKktfynluMDGA0XAKVUd56g7QNlIRwbPPmzbSBHSGE6AoDLbpjOKkJ0QFKQgghhBgODrpjiP7QFF1CCCGE6AW1hBBCCDEcHMyOIfpDSQghhBDDQUmIUaEkhBBCiOGgJMSoUBJCCCHEcMhkgMYzbmmKbltFSQghhBDDQS0hRoVmxxBCCCFEL6glhBBCiOGglhCjQkkIIYQQw8FosViZgJKQtoqSEEIIIQaDYWSaL4BKe8e0WZSEEEIIMRzaLNtOLSFtFiUhHKNddAkhRIdoTIhRoSSEY7SLLiGEEKIaJSGEEEIMhzaLlQloTEhbReuEPO+EDEy71wFCI22upOczbMb+fAIGwk5lNGZBE/LuGE2/SJtEScjzTgiY9pAa77+E5+H5utcb9/MZ9e+PgUnncuNNsjjAyBgwMpmGX/TzbauoO4YQQojhYBgtpuhSEtJWURJCCCHEcMgYLcaEUBLSVhlrIychhBBC2jhKQtqgEydOaHWdSZd6HdeEm/L4vk5b9Hxt4zptmXjU8XqdsEOFVtdpS9u/E9pe12YwMu2+SJtESUgbpH0SItVxTbgpj+/rtMX783nw/SatZTKh7XW8//60TEK0vE7YoVKr67T1vCYhDQNTNf8ibRONCSGEEGI4GBnACDS8hpKQtoqSEEIIIQaDYbQYmEpToNssSkIIIYQYDmoJMSqUhHCE+f//6CUSicbXymQyra4DGMBUwxebCaP4/1yXx/d19Hy6LY/v6wzt+Uw0/IgulCn+vwZkMv7+vsjjmbbwZm4KaLxQCL3TtVkCpk38qzI+T58+RVhYmL6rQQghOvPLL7/A2dlZL2XX1tZi9uzZKC4u1up6BwcH/PzzzxCJRDquGWkNSkI4IpPJUFRUBEtLSwgEGjYdEkJIG8IwDKqqquDo6AihUH+TKmtra1Ffr93sLFNTU0pA2iBKQgghhBCiF7ROCCGEEEL0gpIQQgghhOgFJSGEEEII0QuauGQkUlNTERMTg6SkJOTn56Ndu3bw8fHBjBkz0KlTJ43j5HJzcxEZGYk7d+6gvLwcLi4uGD58OF577TVYWFi0uefLysrCvn37kJ6ejuLiYpibm8PDwwOTJk3CgAEDlO5bV1eHPXv24Ny5c6ioqICnpyemT5+OwMBA3p5Njotn1PT3bUjP9qz9+/cjMjISHh4e2L59O9ePpIDL50tPT8e+fftw584d1NbWomPHjnjllVcQGhrK1+MRwhkamGok1q1bh7t370IsFsPT0xMlJSU4fvw4qqursWnTJnTp0kWjOAAoKCjA/PnzYW1tjdGjR6Ndu3ZISUnB2bNnERQUhOXLl7e557t69SqOHTsGX19fODo6oqamBgkJCbh9+zY++ugjjB49WuG+GzduRHx8PEJDQ+Hm5oazZ88iLS0Na9euhb+/P2/Px9UzavL7NrRna+zp06f48MMPIRAI0L59e96TEK6e7/r161i9ejW8vb0xZMgQWFpa4vHjx2AYhpYAIMaBIUbhzp07TG1trcKxnJwc5rXXXmM2bdqkcRzDMMz+/fuZ8ePHMw8ePFA4vnnzZmb8+PFMeXm5jp+iaZrU+1n19fXM/PnzmQ8++EDh+L1795jx48czhw4dYo/V1NQwc+bMYRYvXqy7yquJi2dszT11iYtna+ybb75hli1bxixdupSZO3euTuqsCS6er7Kykpk+fTqzdu1aRiqV6rzOhLQFNCbESPj5+cHMzEzhmJubGzw8PJCdna1xHPC/VRLt7e0Vjjs4OEAoFMLUlL/ePE3q/SwTExM4OzujslJxl9P4+HgIhUKFT58ikQghISFISUlBQUGB7h5ADVw8Y2vuqUtcPJvcrVu3EB8fjzlz5uisvpri4vni4uJQUlKCGTNmQCgUorq6GjIZbUlPjAslIUaMYRiUlJTA1tZWq7iePXsCAL799ltkZGSgoKAAFy5cwB9//IHx48fzOiZEleaer7q6GqWlpXj8+DGOHDmCa9euoXfv3goxGRkZ6NSpE6ysrBSO9+jRAwCQmZnJXeXV1Npn1PSefNLFs0mlUvzwww94+eWX4enpyUOt1dfa5/v7779hZWWFwsJCfPjhh5g8eTKmTJmCHTt2oLa2lq/HIIRTNDDViMXGxqKwsBBvv/22VnH9+vXD9OnT8dtvv+Hy5cvs8TfffBMzZszgpM6aaO75du7ciejoaACAUCjEoEGD8OGHHyrEFBUVwcHBQela+bHCwkIOaq2Z1j6jpvfkky6eLTo6GgUFBVizZg3n9dVUa58vNzcXUqkUa9asQUhICN555x0kJyfj+PHjqKysxJIlS3h5DkK4REmIkcrOzsb3338PX19fBAcHax3Xvn17BAQEYPDgwWjXrh2uXr2KAwcOwMHBAePHj+fyEZrVUr1DQ0MhFotRVFSECxcuQCaToa6uTiGmtrZWqQkdALu0s74/beriGTW9J1908WxlZWXYs2cPpkyZAjs7O76qrhZdPF91dTVqamowZswYfPDBBwCAwYMHo76+HtHR0Xj77bfh5ubGy/MQwhXqjjFCxcXF+Oqrr2BlZYWlS5fCxMREq7jz58/j3//+N+bPn49XXnkFgwcPxoIFCxAcHIzdu3ejrKyMj8fRuN4A4O7ujj59+iA4OBhffvklqqqqsHr1aoVdQEUikco3bXnyoc99JnT1jJrekw+6erbIyEjY2NjoNRlWRZf/PgFg2LBhCtcOHz4cAJCSksLhUxDCD0pCjExlZSVWrlyJyspKrFq1Ck5OTlrHnTx5Et7e3kq7ZgYFBaGmpgYZGRmcPENz1H2+Z4nFYqSlpSEnJ4c95ujoqHJHTvkxde+ta7p8xtbeU9d09Wy5ubk4deoUJkyYgKKiIjx58gRPnjxBXV0dpFIpnjx5gvLyci4fRSVd//sElAeGy1t9KioqdFNpQvSIumOMSG1tLVavXo2cnBysWbMGHh4erYorKSmBjY2N0nH5LpZSqVR3lVeDuvVu6lrgfzN+AMDLywtJSUmQSCQKg1Pv3bvHnuebrp+xtffUJV0+W2FhIWQyGX788Uf8+OOPSvGzZ89GaGgorzNmdP2769atG/7++28UFhaic+fO7PGioiIA0PvAYkJ0gZIQIyGVSrFhwwakpKRg+fLl8PX1bVUc0DDF8MaNG8jJyVFY9fH8+fMQCoW8zkZQt94lJSVKnxzr6+sRExMDkUgEd3d39rhYLEZUVBSio6MxadIkAA0rqJ45cwY+Pj5wcXHh7HlU4eIZNfl9c0nXz+bh4YFly5YpXR8ZGYmqqirMmTMHrq6uOn+OpnDxuxsyZAgOHjyI06dPK8yc+fPPP2FiYsLOXiPEkFESYiR27dqFy5cvY8CAASgvL8e5c+cUzo8YMUKjOACYNGkSrl27hqVLl2LcuHFo164drly5gmvXruHll1/mtUlf3Xpv374dEokEAQEBcHR0RElJCWJjY/Ho0SO89957sLS0ZK/x8fGBWCxGREQESktL4erqipiYGOTn52PBggW8PZscF8+oye+bS7p+Njs7OwwaNEipnKNHjwKAynNc4uJ35+3tjZCQEJw+fRpSqRQBAQFITk5GfHw8Jk+erLcuNUJ0iZZtNxKfffYZbt261eT5Y8eOaRQnl5qair179yIjIwPl5eXo0KEDgoOD8frrr/M6sFHdep8/fx6nT5/GgwcPUF5eDktLS3Tr1g3jx49HUFCQ0nW1tbWIjIxEbGyswt4xffv25exZmsLFM2r6++YKV78/VeWUlZXxvmw7V89XX1+PAwcO4MyZMygqKoKLiwvGjRuHV199lbNnIYRPlIQQQgghRC9odgwhhBBC9IKSEEIIIYToBSUhhBBCCNELSkIIIYQQoheUhBBCCCFELygJIYQQQoheUBJCCCGEEL2gJIQQQgghekHLthMAwIQJExS+FwgEsLKyQpcuXRAcHIyXX34ZAoFAIb59+/bYuXMnr/Xcu3cv9u3bh4ULF2LUqFEaXVtdXY3o6GgkJiYiOzsbFRUVMDc3R+fOndGnTx+8/PLLaN++favqJ1858+eff0aHDh1ada+2rvEqoV9//bXKvUxSUlKwZMkSBAQEYN26dXxXsUX6+ndMCGlASQhREBwcDACQyWTIy8vD3bt3cefOHSQlJWHJkiV6rp327t69i3Xr1qG4uBjm5ubw8fGBvb09JBIJ0tLSsH//fhw+fBgrVqxAnz599F1dg7Nnzx6sX79e39UghBgYSkKIgkWLFil8f+PGDaxatQrnz5/H8OHDMWDAAADAjh07YGpqGP98MjIysHz5ctTW1uL111/H1KlTYWFhwZ6XyWT466+/sHv3bjx9+lSPNTVMIpEIt2/fxs2bNxV2eyWEkJbQmBDSrMDAQHYH0L/++os97u7uzutW6dpiGAabN29GbW0t3nrrLcycOVMhAQEAoVCIwYMHY8uWLejevbueamq4xo4dC6ChNYQQQjRhGB9liV517doVABRaCZ7tS5dKpfjss89w9+5dfPjhhxg3bpzCPW7fvo1ly5bB3t4e3377LWxtbdlz165dw/Hjx5GamgqJRAInJycMHDgQb775pkKcNq5du4asrCw4OzvjzTffbDbW2toa1tbWCseqq6tx5MgRXLhwAXl5eTA1NYWXlxfGjh2LYcOGqVWHJ0+eYPbs2U2Oi2hqnMt7772H/Px8HDt2DCdOnMDJkyeRl5cHe3t7jB07FpMmTYJAIEB6ejr27t2Lu3fvor6+Hr1798b777+vNL5ly5YtiImJwddffw2BQIB9+/YhLS0NAODv74+wsDB4eHio9UyNDRw4EElJSbh79y6uX7+u1g7ELY3tafzscsnJyVi2bBmCg4MRFhaGiIgIXLlyBVVVVejatSvCwsLg5+cHAPjjjz9w8uRJ5ObmwtbWFiEhIZg6dSqEQtWfu+rq6vDbb78hNjYWhYWFcHR0xEsvvYQ333wTIpFIKV4qleLUqVOIiYnBw4cPIZVK0alTJ4wcORLjx49X2mFa/jxHjx7F8ePH8eeffyI3NxedOnXC//3f/7X48yLEWFFLCGlRVVUVAMDMzKzJGBMTE4SHh8PS0hK7du1CdnY2e66yshKbN28GwzD4+OOPFRKL3bt3Y+XKlfj777/RqVMnBAUFwcTEBL///jsWL16M4uLiVtX96tWrAACxWKz0xtASiUSCzz77DHv27EFpaSlefPFF+Pn5ITU1FRs3bsSPP/7Yqrqp66effsKuXbvQvn179O7dG+Xl5di9ezf27t2LO3fuYOnSpSgqKkKfPn3g4OCAy5cvY/ny5aipqVF5v8TERHz++eeoqalBv3794OjoiKtXr2Lp0qVa/7zfeustAA3JBdcqKyuxZMkS3Lx5Ez179oSnpyfu3r2LL774AllZWfjxxx/x888/w9nZGb1790ZlZSX27duHyMhIlfdjGAbr1q3D4cOH4e7ujv79+6OiogL79+/HV199BalUqhBfU1ODFStW4LvvvkNubi58fHzQp08fFBcX4+eff8a6desgk8lUlrV9+3bs2rUL9vb2CAoKQseOHXX+8yHEkFBLCGkWwzC4cuUKAMDT07PZ2I4dO+KDDz7A1q1bsWnTJmzatAlmZmb47rvvkJ+fj9DQUAQGBrLxFy9exKFDh9ClSxcsW7YMbm5ubJl79+7Fr7/+ip9++gmffvqp1vXPyMgAAHh7e2t87X//+1+kp6ejV69e+Pzzz2FlZQUAyM7OxrJly3Ds2DH06dOHHSfDlYsXL+Lf//432/2VnZ2NhQsXIioqCjExMXjvvfcwZswYAA2f6FeuXImkpCRcuHBBZSvD0aNHsXTpUgwaNAhAw6f6DRs2ICEhASdOnMD06dM1rmNQUBC6deuGe/fu4erVq+jfv38rnrh5ly9fxksvvYSFCxey45LkLSvffPMNKisrFX5eDx8+xMKFC3H06FFMnjwZlpaWCvcrKCgAwzDYvn07mxSUlpbi888/x82bN3H8+HG8+uqrbPyuXbuQlJSEoUOH4qOPPmJbzyQSCTZu3IjLly/j1KlT7O+ksUuXLmHr1q3o0qULJz8bQgwNtYQQlaRSKXJzc7Ft2zakpKTAzMxMrSmxI0eOhFgsRkZGBiIjIxEbG4u4uDh06dIFM2fOVIj97bffAABLlixhExCgYXrwW2+9ha5duyI+Ph6lpaVaP0d5eTkAwM7OTqPrqqur8eeff0IoFOLDDz9kExCgYTyMvGuncXcBV95++22F8TfyT+s1NTVwdnZWeLMzMzNDaGgogIbuC1WGDRvGJiBAQyvW5MmTATR0m2lr2rRpALgfG2JlZYUPPvhAYWD0q6++CoFAgOzsbKWfl4eHB1588UXU1NQgPT1d5T2nTp2q0CphZ2eHsLAwAMCJEyfY4yUlJfjzzz/h7OyMhQsXKnTfWVlZYcGCBTA1NcXJkydVlvP6669TAkJII9QSQhQ8u14IAFhaWmLRokVqD0SdN28e7t27hyNHjsDc3BxmZmb45JNPFLpzSkpKkJmZCTc3N5V/lAUCAfz8/JCRkYH79++rNc5Al9LT01FbW4tu3brB3d1d6fyIESPw448/4s6dO5DJZE2ONdCFxq1HcvI1SFSdk7+ZNtW1ouoaeRLYmu6vAQMGoHv37khLS0NiYiJnLUTdunWDjY2NwjFra2vY2NigvLy82Z9XUVGRynsOHTpU6Vi/fv1gY2ODx48fo6ioCI6OjkhOTkZ9fT369esHc3NzpWscHBzg5uaGrKws1NTUKMUEBQWp/ZyEPA8oCSEK5OuECIVCdrGywYMHK/3Rb46NjQ0+/PBDrFmzBlVVVZg5cya8vLwUYvLz8wEAubm5KhOfxsrKyjR8iv9p164dAGjcmiJ/s2pqwTEbGxtYW1ujsrISFRUVrR5A2xwnJyelY/IuBVXn5LN/6urq1L6fvKWnqWvU9dZbb2HVqlXYu3cvZ0mIqvoDDT+T8vLyZn9eqp7PxsZGoaWrsfbt26OiooJNQuT/bk+dOoVTp041W0/5YniNubi4NHsNIc8bSkKIgmfXCdHWhQsX2P9W1QQuH7jn4OCg8pNrY635w921a1fcvXsX9+/fZ6catzVNDWKUa66VpfEqturistWmf//+8PHxwb1793Dp0iU4ODhodZ/mfiYtPTOXzyevV9euXVscI6VqHR1VM20IeZ5REkJ0Li4uDnFxcfDw8ICpqSkuXryIF198kW1lAQBnZ2cAgK2trc4SH1X69++PEydOID4+HmFhYWrPkHF0dATwvxabZ1VWVqKyshIikajFViL5m1F1dbXK88a2QNpbb72FL7/8Evv27cPcuXNVxjT3M5FKpSgpKeGyigoqKiogkUhUtoYUFBQA+N+/B/m/2xdeeAEffPABb3UkxFjRwFSiUwUFBfjuu+9gZmaGxYsX45NPPoFIJMIPP/yAvLw8Ns7Z2RmdO3dGdnY2cnJyOKtPv3794OHhgadPn7IDYZsikUiQlZUFoGHcgUgkwv3795Gbm6sUGxsbC6DhzailT962trYwMTHBkydPlKZ71tfXs/uvGIu+ffvCz88PmZmZiI+PVxkjf1NX9buXj7vg08WLF5WOXb9+HeXl5ejYsSNb3169ekEoFCIxMZH3OhJijCgJITojk8mwZcsWVFZWYsaMGfDy8oKHhwfeffddSCQSbN68WeFNeMqUKZDJZFi3bh07lbaxsrKyFvvdWyIQCNhEaO/evfjPf/6j9OmbYRhcvnwZixYtYhfvsrCwQEhICGQyGb777juFa3JycrB//34AqgfyPsvMzAy+vr4oLy9XmGkhlUqxc+dOPHnypFXP2BbJ1w1papaIv78/gIZkrvHz5+Xl8bb+SmP79u1TqEdpaSl++eUXAFBYeM/JyQkhISHIz8/Hxo0bVQ7kzc3NbTL5IoQoou4YojNRUVFITk5G7969MXHiRPb4hAkTcPXqVdy4cQMHDx7ElClTAAAvvfQSHj58iAMHDmDRokXw8vJiZ3Y8fvwYDx48gKWlJV555ZVW1atr165YvXo11q1bh4MHD+LYsWPw9fWFvb09KisrkZ6ejpKSEohEIoXxJ++88w7u3buHv//+G3PmzIG/vz9qamqQlJSE2tpaTJgwQe3Bl1OnTsWXX36Jn376CRcuXICDgwPS09NRU1OD4OBgxMTEtOoZ25o+ffrA39+/ySm/rq6u7HMvXLiQ/dneu3cP/fr1Q01NTZNdYbrm4uICT09PfPTRR+jduzdMTEyQlJSEyspK9OrVSynRnDNnDp48eYKEhARcv34dXl5ecHFxQU1NDR4+fIjHjx8jKCgIYrGYl/oTYsgoCSE6IV8XxMbGBosWLVIYPCgQCPDxxx9j3rx52LdvHwIDA9GjRw8ADW/0ffv2xfHjx3H37l1kZWXB0tISTk5OGDt2rM7+kL/wwgv48ccfER0djcTERDx48AAVFRWwsLBA586dMWbMGLz88stsnz/QMGNk3bp1iIqKwoULF5CYmAhTU1N069YNY8eOxfDhw9Uuv0+fPli+fDn27duH+/fvw8LCAr1798bMmTNx9uxZnTxjW/PWW2/h888/b/L8vHnz4OjoiNjYWFy/fh0uLi5444038MYbb+D999/nrZ4CgQCfffYZ9u3bh7i4OHYmzLhx4/Dmm28qjSMyNzfHypUrERcXh7NnzyIzMxNpaWmwtbVF+/btMWLECLWX9CfkeSdgGIbRdyUIIYQQ8vyhMSGEEEII0QtKQgghhBCiF5SEEEIIIUQvKAkhhBBCiF5QEkIIIYQQvaAkhBBCCCF6QUkIIYQQQvSCkhBCCCGE6AUlIYQQQgjRC0pCCCGEEKIXlIQQQgghRC8oCSGEEEKIXlASQgghhBC9+H/ixAwPSvj44wAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "pixelfile.plot(frame=1);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "It looks like the star is an isolated object, so we can extract a light curve by simply summing up all the pixel values in each image:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:10:25.115768Z", "iopub.status.busy": "2023-11-03T14:10:25.115587Z", "iopub.status.idle": "2023-11-03T14:10:25.137086Z", "shell.execute_reply": "2023-11-03T14:10:25.136779Z" } }, "outputs": [], "source": [ "lc = pixelfile.to_lightcurve(aperture_mask='all');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The above method returned a `LightCurve` object which gives us access to the number of photons received by the spacecraft over time (known as the *flux*). The time is an [AstroPy Time object](https://docs.astropy.org/en/stable/time/) in units of *days*:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2023-11-03T14:10:25.138959Z", "iopub.status.busy": "2023-11-03T14:10:25.138841Z", "iopub.status.idle": "2023-11-03T14:10:25.141387Z", "shell.execute_reply": "2023-11-03T14:10:25.141142Z" } }, "outputs": [ { "data": { "text/plain": [ "