{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 1. Preprocess seismic waveforms\n", "\n", "An essential step before detecting and locating earthquakes is the adequate preprocessing of the seismic data. Preprocessing addresses specific needs for a given data set and task and, thus, cannot be straightforwardly generalized to new data sets. Preprocessing workflows must remain as simple as possible. The preprocessing steps taken in this notebook are given as an example and should not be blindly reproduced on other data sets.\n", "\n", "### Contents\n", "\n", "* [Read and destination paths](#read-and-destination-paths)\n", "\n", "* [Show an example raw trace](#show-an-example-raw-trace)\n", "* [Process](#process)\n", "* [Compare raw and processed trace](#compare-raw-and-processed-traces)" ] }, { "cell_type": "code", "execution_count": 43, "metadata": {}, "outputs": [], "source": [ "import obspy\n", "import os\n", "import tqdm\n", "\n", "import matplotlib.pyplot as plt\n", "\n", "from glob import glob" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read and destination paths\n", "\n", "First, we define the path from which we read the raw data and save the preprocessed data. Also extract the path to all input waveforms." ] }, { "cell_type": "code", "execution_count": 44, "metadata": {}, "outputs": [], "source": [ "DIRPATH_RAW = \"../data/raw/\"\n", "DIRPATH_PROCESSED = \"../data/processed/\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Again create directory if needed and copy the metadata therein." ] }, { "cell_type": "code", "execution_count": 45, "metadata": {}, "outputs": [], "source": [ "# Create directory\n", "os.makedirs(DIRPATH_PROCESSED, exist_ok=True)\n", "\n", "# Copy meta to destination\n", "!cp {DIRPATH_RAW}*xml {DIRPATH_PROCESSED}" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Show an example raw trace\n", "\n", "We collect all available waveform files to manipulate them afterward." ] }, { "cell_type": "code", "execution_count": 46, "metadata": {}, "outputs": [], "source": [ "filepaths_raw = sorted(glob(os.path.join(DIRPATH_RAW, \"*.mseed\")))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Using `obspy`, we can visualize the waveforms. In the following cell, the selected waveform has not been preprocessed yet: the amplitudes are in counts (signed integers) and the waveform exhibits low-frequency components, especially with a day-long period. Note that, here, we do not show the earthquake itself to see the ambient noise better." ] }, { "cell_type": "code", "execution_count": 47, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAj0AAADtCAYAAABUMMQHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAABU3klEQVR4nO3deVwU9f8H8Beici0siHgg4JWRIGLemqaVWp6peSKmmFemZeaR5q2oWWb1MzXLpBQ8MjUrzzQlr8BUTg9EEFPBA9h1BRZY3r8/fOx89z6AhYV9Px+PeTx25jPzmc98ZnbnvZ/PHHZERGCMMcYYq+ZqVHYBGGOMMcYqAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwi5DL5QgLC4OPjw/EYjF69uyJhIQEtXnWrFkDLy8v1KlTB3PnzoXyjSjFxcV466230KhRI9jZ2SEzM1NtuVmzZqFZs2ZwdXVF+/btER0dbbAssbGxCA4OhrOzM3r06IHbt28LaYGBgRCJRMJQo0YNrFu3Tmc+169fx4ABA1C3bl14eXkhNDQUOTk5avNs27YNLVq0gIuLC1q2bIkbN25g1apVQv4ODg6oVauWMD516lQAwOHDh/Hcc8/BxcUFb775plq+Dx8+RP/+/eHs7Ax/f3+cOHHC4Pbqq1djdaEpPz8foaGhcHV1hZ+fH3bu3KmWHhERAR8fH7i5uSEsLAyFhYUGy1Ua58+fxyuvvAIPDw80bNgQ06dPV1uPobqJjo7Gyy+/DBcXF7zxxhtq+ZqyLzXpq9e///5b7RhydnZGjRo18PDhQ535REREoE2bNnB1dUWzZs2wefNmIc1QXqrHao0aNeDk5CSM//333ygpKcHMmTPh7u6O+vXrY/369WrrNXSMaTJ2zBk6xjQZOuaMHWOMlTtizAJkMhktX76c7ty5Q8XFxbRu3Tpq1qyZkP7HH3+Qn58fpaam0r1796hly5a0detWIiIqKiqiL7/8ks6fP08A6P79+2p5L126lFJSUkihUNCePXvIw8ODpFKpznIUFBSQj48Pbd26lfLz82nu3LnUvXt3nfM+evSIateuTSkpKTrT//nnH/rxxx8pNzeXZDIZDR8+nCZMmCCkHzx4kAIDAykxMZFKSkro5s2b9PjxY7U8Vq9eTePGjVOblpWVRe7u7nTo0CGSyWQ0ZswYGjt2rJA+fPhwmjhxIj19+pT2799PHh4elJ2drbOMhurVnLogIpozZw717duXJBIJnT17lsRiMV2/fp2IiOLj48nDw4NiY2MpNzeXevbsSYsWLdKbV2kdPnyY9u/fTzKZjB49ekTdu3en5cuXC+mG6iY2NpYiIyNpxYoV9Prrr6vla2xfajJUr5q++uor6tatm968Nm/eTOfPn6eioiJKTEykevXq0enTp83Ky9/fn/766y+1ad988w29+OKLlJWVRdeuXaMGDRrQiRMniMj4MabJUL2aUxfGjjlDxxhjlsBBD6sQcrmc7Ozs6NGjR0RENGrUKFqzZo2QvnXrVnrllVe0ltMV9Gjy9vamixcv6kw7cuQIvfDCC8K4TCYjJycnSk9P15p348aN1LlzZ5O2h4jo2LFj1KpVK2G8Q4cOdOzYMYPL6Ap6Nm/eTG+88YYwfvPmTXJ0dKSCggJ68uQJ1a5dm+7duyekd+/enX788Ued+RuqV3PqgoioQYMGdOHCBWF87NixtGzZMiIi+vjjj2nq1KlC2okTJ6hp06YGt708bNmyhQYMGEBEZHLdbNu2TSvo0aS5LzWZerwSEXXs2JE2b95sdFuUQkJC6PPPPzcrL11BT+fOnWnXrl3C+KJFiygsLIyIDB9jmozVqzl1YeyYM3SMMWYJ3L3FKsT58+dRv359eHp6AgCSk5MRFBQkpAcHByMpKcnsfNPT05GdnY3nnntOZ7rmelxcXNC8eXMkJydrzRsZGYkxY8aYvO5z584hMDAQAKBQKHD58mUkJCTAx8cHTZs2xfLlyw02++srY/PmzVGzZk3cunULKSkpEIvFaNiwoZCuWldnzpyBu7u73rxU5zVWF2vWrMGAAQMAADk5OcjMzDQ5r+DgYKSlpSE/P9/o9paFap0bq5vS5quLqcfrzZs3ceXKFQwfPtyk9SoUCsTExOhct7l5mbPvVY8xAJg2bRqmTZsGwHi9GquL1q1bIyoqSue8qsecsWOMMUuoWdkFYNWfRCLBlClTEB4eLkyTyWRwc3MTxt3c3CCTyczKt6ioCOPGjcOcOXMgFot1zqO5Hn3rSk9PR0xMDPbt22fSuq9cuYKvv/5auJ4oKysLxcXFOHHiBBITEyGVStG3b1/4+voiLCzMYF4ymQxeXl46y1hQUKCz/Lm5uQCAbt26CZ91ba/qthqri48//lgtH3t7ezg7O5uUl/KzTCaDk5OTwe0trcOHD+Pw4cOIi4szuD2q9WEKzX2pi6nHa2RkJN544w3UqVPHpHUvXLgQjRo1wuuvv17mvIzte33HGABs3LhRbz7KeZX1aqwu4uPjjeYlk8mMHmOMWQK39DCLKigowODBg9G/f39MmDBBmC4SiSCVSoVxqVQKkUhkcr5EhPHjx6NevXpYunSpMF31Ys+MjAyt9ehbV1RUFHr16oV69erpzUspLS0NAwcOxNatW4V/6MoT/bx58+Du7g4/Pz+89957OHTokNFtMVRGU8uvLy/Vec3JSyQSQaFQIC8vz6S8lJ+N7UPVC3X79u1rcF5VsbGxCAsLw4EDB1C/fn2zt0cfXftSVxlNPV6joqLUWgsNbe/mzZuxb98+7N27F3Z2dkbzMqY8972hec357ho7tg0dY4xZAgc9zGKKi4sxatQoeHt74/PPP1dLCwgIULubKy4uzmD3gqYZM2bg3r172LFjB2rU+N9hnJSUJPyL9PPz01rP06dPkZqaioCAALX8dJ1gNPMCgMzMTPTu3RuLFi3C4MGDhXk9PDzg7e2ttrwpXVuAdl3cunULxcXFaNasGVq0aAGJRKJ2B5uhujJUr6bWhXJ7GjRoYHJecXFxaNq0qdFWnu7duwt1evjwYYPzKl29ehWDBg3CDz/8gM6dOwvTza0bTfr2pa4ymnK8xsbG4v79+xg4cKDR7d29ezfCw8Nx9OhR1K1bV6tsuvIyxpx9r3qMaTJWr+Z8dw0dc8aOMcYsopKvKWLV2Pjx46lPnz5UWFiolfb7779T48aN6datW3T//n0KDAxUuwOkoKCA8vPzCQClp6dTfn6+kLZ48WJq06YNSSQSo2VQ3j2ybds2KigooI8//ljrjqXLly+Ti4sLyWQyg3nl5uZS69at9V5ouWDBAurfvz9JpVK6e/cutWrVirZt26Y2j6G7t44cOUJPnz6l0NBQtTtrhg0bRpMnT6a8vDz69ddfDd69ZaheTakLVbNnzxa25/z58yQWi+natWtE9OzurTp16tC///5Lubm59Oqrr1rk7q2MjAzy9fWliIgInemG6kahUFB+fj5t2bKFevfuTfn5+cKxaGxfajJ2vBIRffDBB1r7VpejR4+Sl5cXxcXF6Z3HWF66LmTesGEDtWvXjh48eEDXr1+nhg0b0p9//klExo8xTYbq1ZS6UDJ2zBk6xhizBA56mEWkp6cTAHJ0dCQXFxdhiI6OFuZZtWoVeXp6kru7O82ZM4dKSkqEtMaNGxMAtUEJADk4OKjlu2PHDr1liYmJoaCgIHJ0dKTu3btr3a00e/ZsCgkJMbpNERERBEBtvS4uLkK6XC6niRMnkpubG3l7e9PixYvVtolId9BD9Ow24GbNmpGTkxMNHDhQLah58OAB9e3bl5ycnKhFixZ0/PhxIS06OlqtDESG69VQXYSHh6vd4ZOXl0chISHk4uJCPj4+FBkZqbaebdu2kbe3N4lEIho3bpzOO4HKaunSpWRnZ6dW3wEBAUK6obr566+/tI4hZd0b25e6GKrX4uJiatCggdG794iIevbsSTVr1lRb75QpU8zKS1fQo1Ao6IMPPiCxWExeXl60bt06tXRDx9iUKVPUymCoXo3VRUBAgNr30dAxZ+wYY6y82RGZ2AbPGGOMMVaF8TU9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsAgc9jDHGGLMJHPQwxhhjzCZw0MMYY4wxm8BBD2OMMcZsQs3KLoA1KCkpwb179+Dq6go7O7vKLg5jjDFmE4gIT548gbe3N2rUsHw7DAc9AO7duwdfX9/KLgZjjDFmk+7cuQMfHx+Lr4eDHgCurq4AnlW6m5tbJZeGMcYYsw1SqRS+vr7CedjSLBL0iEQitfGnT59i7969eOuttxAREYGJEyfC0dFRSE9OToafnx8AIDY2FhMnTkRKSgo6dOiAn376CY0bNwYA5OfnY9KkSfj111/h4eGBTz/9FKNHjxbyiYiIwMKFCyGVSvHWW2/h22+/Re3atY2WV9ml5ebmxkEPY4wxVsEq6tISiwQ9MplM+BwfH4/OnTujT58+wrRevXrhyJEjWsvJ5XIMHToUy5YtQ0hICJYsWYKxY8ciOjoaALBkyRJkZ2fj7t27SExMRL9+/dCuXTs8//zzSEhIwKxZs3Ds2DG0aNECgwcPxsqVK7F8+fJy2aa0tDRcvXoVRUVF5ZKfratZsyaaNm2Kli1b8nVUjDHGKoTFu7ciIyPx5ptvmtR0derUKYhEIkyYMAEAsHjxYnh5eeH27dto3Lgxtm/fjgMHDsDNzQ1du3bFoEGDsGvXLixevBhRUVEYOXIk2rdvDwBYtGgRJk6cWOagRyqVYtasWbhy5UqZ8mG6NWnSBF9//TW8vb0ruyiMMcaqOYsGPUSEnTt3YuPGjWrTz549C09PT9SvXx/vv/8+pk6dCuBZN1dQUJAwn4uLC5o3b47k5GS4ubkhMzNTLT04OBgxMTHCsq+//rpaWlpaGvLz8+Hk5KS2frlcDrlcLoxLpVK927Bs2TLcunULa9euRZcuXbTyYqVTWFiIuLg4hIeHY+bMmdi9eze3+DDGGLMoi94fFh0djby8PLVgpEePHkhISMDDhw+xbds2LF++HPv37wfwrFtM85oaNzc3yGQyyGQy2Nvbw9nZWStN17LKz6pdbUqrV6+GWCwWBn13bkmlUpw5cwZTp07Fq6++WqqAJy8vD4mJicjJyTF72eqsdu3a6NChA+bNm4dbt24hJSWlsovEGGOsmrNo0BMZGYkRI0agVq1awrSmTZuiSZMmqFGjBjp16oT3339fCHpEIpFWq4tUKoVIJIJIJIJCoUBeXp5Wmq5llZ81L6oGgPnz50MikQjDnTt3dJb/v//+g0KhQOvWrUtZA0BqaioKCgqQmppa6jyqs+DgYABAenp65RaEMcZYtWexoKewsBB79+7FmDFjDBdA5WFEAQEBSEhIEMafPn2K1NRUBAQEwMPDAw0aNFBLj4uLQ2BgoM5l4+Li0LRpU52tMw4ODsKdWobu2FIoFACgFrSZq6SkpNTL2oKaNZ/1sBYXF1dySRgrH0SETZs2CTdgMMash8WCnkOHDkEsFqNr165q048cOYKHDx8CAC5duoSvv/4aAwYMAAD07NkTMpkMERERkMvlWLlyJdq3by/csh4aGooVK1bgyZMnuHDhAg4ePIiRI0cCAEJCQrBnzx5cunQJEokE4eHhCA0NLfftGjZsGD788EO1aZ9//jl69OiBxo0b48KFC2ppK1euxJYtW3Tm1bNnTzg6OsLV1RXu7u7o1q0bIiIitOb77rvvEBQUBBcXF/j5+WHcuHFCy0h+fj5CQ0Ph6uoKPz8/7Ny5U23Za9euoXfv3nB1dUX9+vXx1VdfmbSdS5cuRa1atSASieDq6op27drhzJkzaunKa7GULly4gCZNmgjjTZo0gbOzs9BSJxKJsHXrVpPWz1hVderUKUybNg09evSo7KIwxjRYLOiJjIxESEiI1sWpx48fR2BgIEQiEUaPHo158+ZhxIgRAJ61wOzbtw9ffPEF3N3dcfbsWWzfvl1Ydvny5RCLxWjYsCGGDx+OjRs3wt/fHwAQFBSEdevWYeDAgfDx8YGvry8++eSTct+u//u//8P27dtx6dIlAEBGRgZWr16Nb7/9tlQX4kZERODJkye4ffs2Zs6ciQULFuDjjz8W0leuXInFixfj008/xePHj5GcnIyXXnoJJ0+eBKB+G/+uXbvw7rvv4saNGwCeBUR9+/bF22+/jUePHuHmzZtqjw4w5p133oFMJoNUKsW7774r7CdznDx5UrgmSyaT4Z133jE7D8aqkrS0tMouAmNMD4vdvfXzzz/rnL5u3TqsW7dO73IdOnRAfHy8zjQnJydERkbqXXb8+PEYP368WeU0V8OGDbFy5UpMnjwZ//zzD2bMmIH3338fL7zwgs75TQ2ExGIxhg0bhtq1a2PYsGH46KOPUKtWLaxatQpRUVHo16+fMO/kyZOFz4Zu49+2bRu6du2KsWPHAngWVLZs2dLsbbazs0NISAgmTZoEiUQCsVhsdh6M2YqKeH8QY6x0+NtZClOmTIGjoyOGDh2KlJQUzJ8/v9zy7tevH0pKShAbG4vz58+jsLBQ6P7TlJOTo/M2/qSkJADPnm5dp04ddOrUCV5eXhg0aJDei7YNUSgUiIyMhL+/Pwc8jBnBQQ9j1ovfvVUKdnZ22LJlCwIDA3Hy5Em1V1307t0b9vb2wnheXp5ZrU81a9ZE3bp1kZOTAyJC3bp1hYt9NRm7jf/u3buIiYnB8ePH0bp1ayxcuBChoaE4ffq0SWXZtm0bdu3aBblcDiJCVFSUznQlhUIBT09PtXk06+PIkSPo3LmzSetnrCpSPd4ZY9aF/5KUUkBAAABodRcdP34cubm5wjB48GCz8i0uLsajR4/g4eEBT09PPHr0SO+dTcZu43dycsLQoUPRoUMHODg4YOnSpfj77791PrtIl7CwMOTm5iIvLw9//fUXJkyYgMuXL2ulK4fjx49r5aFZHxzwsOqOW3oYs1787bQyhw8fRo0aNdCxY0d06dIFtWrVwh9//KFzXmO38bdq1UrnckRkVpns7OzQpUsXtGjRAidOnDBrWcZsDQc9jFkv/nZaCalUin379mHq1KmYNWsW6tatC3d3d3zyySeYNm0ajhw5Arlcjry8PHz//ff44YcfABi+jT80NBQHDx7E5cuXUVRUhBUrVuDll1826T1ommJjY5GcnFyqC6EZsyUc9DBmvfjbaWGqd29FRkYKrTBK48ePh0gkgq+vL9atW4eVK1dizZo1QvrChQuxZMkSzJkzBx4eHvD398fff/+N1157DYDh2/hbtmyJb775BkOGDIGXlxcSEhLUHgEQGBgo3A2XkZEBkUiEjIwMIX3r1q0QiURwcXHB8OHDsXjxYvTv39+s7X/11VfVntOzYsUKs5ZnrKrha3oYs152ZG5fRzUklUohFoshkUjUns6ckJCAsLAw7N69G82bNy9V3snJycI1N8o3wLP/KSwsRNeuXbF8+XK12/IZq6oOHDiAIUOGADC/K5kxW6Pv/Gsp3NJjYfzmcMZsC3dvMWa9+NtpAv63xhgzFQc9jFkv/nYa4OLiAgDIzs4udR7c0mOYsm6Vdc1YVcfX9DBmvTjoMaBJkybw8vLCb7/9VurWHg56DDt48CBq1qyJF198sbKLwli5UG3p4VZixqwLP5HZgBo1amDKlClYuXIl7t27hy5dusDJycmsPO7fv4/8/HwAwPXr1y1RzCqpsLAQcXFxOHPmDMaNG1chF7AxVhFUgx6FQqH3ieqMsYrH30YjBg8eDLFYjP3792PXrl0oKioya/mUlBQ8efIEANC2bVtLFLFKqlmzJpo2bYoFCxYId7owVh2odm+VlJRUYkkYY5o46DHBK6+8gldeeaVUy77++us4duwYAJj8zivGWNWl2tLTv39/bNiwQXh2FmOscvE1PRbGFzUyZltUg54///zT7Ad6MsYsh4MeC+OghzHbonnLempqaiWVhDGmiYMeC+NndjBmW/iPDmPWy2Jn5J49e8LR0VF451Lfvn2FtDVr1sDLywt16tTB3Llz1W7rjI2NRXBwMJydndGjRw/cvn1bSMvPz0doaChcXV3h5+eHnTt3qq0zIiICPj4+cHNzQ1hYGAoLCy21eSbjH0DGbAv/0WHMeln02xkREQGZTAaZTIbDhw8DAA4dOoRNmzbhn3/+QVJSEn7//Xds27YNACCXyzF06FB88MEHyM7ORufOnTF27FghvyVLliA7Oxt3797Frl278O677+LGjRsAnr0na9asWThw4ADu3LmD9PR0rFy50pKbZxIOehizLRz0MGa9KvzbuX37dkybNg3NmjVDw4YNMXv2bOzYsQMAcOrUKYhEIkyYMAGOjo5YvHgxLl68KLT2bN++HUuWLIGbmxu6du2KQYMGYdeuXQCAqKgojBw5Eu3bt4dYLMaiRYuEfDXJ5XJIpVK1wVL4B5Ax28Lfecasl0W/nTNmzICXlxd69+6N+Ph4AM/eOh4UFCTMExwcjKSkJJ1pLi4uaN68OZKTk5GTk4PMzEyTlw0ODkZaWprwYEBVq1evhlgsFgZfX9/y3XAV3NLDmG3hp7AzZr0sFvSsXbsWaWlpyMjIQO/evdGvXz+hq0v16btubm6QyWQAoJWmmi6TyWBvbw9nZ2eTllV+Vqarmj9/PiQSiTDcuXOn/DZcA//rY4wxxqyDxc7IHTt2hEgkgpOTE+bOnQuRSISYmBiIRCK17iSpVAqRSAQAWmmq6SKRCAqFAnl5eSYtq/ysTFfl4OAANzc3tcFSuKWHMcYYsw4V1gyhbPEICAhAQkKCMD0uLg6BgYE6054+fYrU1FQEBATAw8MDDRo0MHnZuLg4NG3a1Ox3ZZU3DnoYY4wx62CRoCc3NxfHjx+HXC5HYWEh1q9fj+zsbLRv3x6hoaHYtGkT0tLSkJmZiS+++AKhoaEAnt3mLpPJEBERAblcjpUrV6J9+/Zo3LgxACA0NBQrVqzAkydPcOHCBRw8eBAjR44EAISEhGDPnj24dOkSJBIJwsPDhXwrE3dvMcYYY9bBImfkoqIizJ8/H56enmjQoAF+++03HDp0CG5ubujfvz+mTJmCDh06oGXLlujXrx/CwsIAPOt22rdvH7744gu4u7vj7Nmz2L59u5Dv8uXLIRaL0bBhQwwfPhwbN24U3mkTFBSEdevWYeDAgfDx8YGvry8++eQTS2yeWbilhzHGGLMOdqT6ZEAbJZVKIRaLIZFIyv36nmnTpmHTpk0AAK5qxqq/S5cuoV27dmrT+LvPmG6WPP/qwn0vFsbdW4wxxph14DOyhXH3FmOMMWYdOOixMG7pYYwxxqwDn5EtjFt6GGOMMevAQY+FqQY9fDEjY4wxVnk46LEw1e6tkpKSSiwJY4wxZts46LEw1ZYehUJRiSVhjDHGbBsHPRbGQQ9jjDFmHTjosTDu3mKMMcasAwc9FsYtPYwxxph14KDHwlRbejjoYYwxxioPBz0WptrSw91bjDHGWOXhoMfCuKWHMcYYsw4c9FiYnZ2d8JmDHsYYY6zycNBTgbh7i1V1mZmZSE5OruxiMMZYqdSs7ALYEm7pYVVdw4YNAQDp6elo3LhxJZeGMcbMwy09FYhbelh1ERMTU9lFYDasoKAABQUFlV0MVgVZJOiRy+UICwuDj48PxGIxevbsiYSEBABAREQEatasCZFIJAwZGRnCsrGxsQgODoazszN69OiB27dvC2n5+fkIDQ2Fq6sr/Pz8sHPnTrX1RkREwMfHB25ubggLC0NhYaElNq/UuKWHVRfW9t1itkOhUMDLywtisRjFxcWVXRxWxVgk6CkuLkazZs1w4cIFZGdnY9CgQRg8eLCQ3qtXL8hkMmHw8/MD8CxYGjp0KD744ANkZ2ejc+fOGDt2rLDckiVLkJ2djbt372LXrl149913cePGDQBAQkICZs2ahQMHDuDOnTtIT0/HypUrLbF5pcZBD6suOOhhlUV53igsLERWVlZlF4dVMRYJelxcXLBo0SL4+PjA3t4e06dPR1paGh4/fmxwuVOnTkEkEmHChAlwdHTE4sWLcfHiRaG1Z/v27ViyZAnc3NzQtWtXDBo0CLt27QIAREVFYeTIkWjfvj3EYjEWLVqEHTt26FyPXC6HVCpVGyoCd2+x6oKDHlZZatWqJXzm45CZq0Ku6Tl//jzq168PT09PAMDZs2fh6emJgIAAbN68WZgvOTkZQUFBwriLiwuaN2+O5ORk5OTkIDMzUy09ODgYSUlJOpcNDg5GWloa8vPztcqzevVqiMViYfD19S33bdaFW3pYdSGXyyu7CAwAEVV2ESqc6mNAOOhh5rJ40CORSDBlyhSEh4cDAHr06IGEhAQ8fPgQ27Ztw/Lly7F//34Az5ot3dzc1JZ3c3MTmjPt7e3h7OyslaZrWeVnZbqq+fPnQyKRCMOdO3fKd6P14KCHVRd8sql8v/32Gzw9PfHHH39UdlEqDR+HzFwWvWW9oKAAgwcPRv/+/TFhwgQAQNOmTYX0Tp064f3338f+/fsxZMgQiEQira4mqVQqXPCsUCiQl5cnBD7KNABayyo/K9NVOTg4wMHBoXw31gTcvcWqCz7ZVL5BgwYBAAYMGGCTLT4AH4fMfBZr6SkuLsaoUaPg7e2Nzz//XH8BVF7TEBAQINzlBQBPnz5FamoqAgIC4OHhgQYNGqilx8XFITAwUOeycXFxaNq0KZycnMpzs8qEW3qYMQsWLMC0adMquxhGcfcWswYc9DBzWSzomTRpEvLz8xEREaHWB3vkyBE8fPgQAHDp0iV8/fXXGDBgAACgZ8+ekMlkiIiIgFwux8qVK9G+fXvhIWihoaFYsWIFnjx5ggsXLuDgwYMYOXIkACAkJAR79uzBpUuXIJFIEB4ejtDQUEttXqlw0MOMWb16NTZt2iTclWit+GTDrAEfh8xcFgl6bt++jYiICERHR8PDw0Ponvr7779x/PhxBAYGQiQSYfTo0Zg3bx5GjBgB4Fm30759+/DFF1/A3d0dZ8+exfbt24V8ly9fDrFYjIYNG2L48OHYuHEj/P39AQBBQUFYt24dBg4cCB8fH/j6+uKTTz6xxOaVGndvMVNJJJLKLoJBfLKpfKp/Jm0VH4fMXBa5pqdx48Z6+5i7d++OdevW6V22Q4cOiI+P15nm5OSEyMhIvcuOHz8e48ePN6usFYlbepiprP1ps3yyqXwODg5Wf5xYGh+HzFz8GooKxC09zFTWfjLja3oqn6OjY2UXodJx0MPMxUFPBeKWHmYqXc+XsiZ8sjGPJe6uqow7UK2NNR+HmZmZePr0aWUXg2ngoKcCWUvQk5WVhQ0bNlj9dSO2TLWlJz4+HpmZmeWW9+3bt5GTk1OmPKz5ZGONLPHd55Ye6z0OHzx4gIYNG6J+/fqVXRSmwaLP6WHqrKV7680338Q///yDEydOCA+GZNZFGfTcuHEDwcHBAMqnteD+/fto0qRJmfOz1pONtVIoFKhZs3x/bjnosb7jsLi4GP/++6/wJ4VbeqwPt/RUIGtp6fnnn38AAAcOHKjcgjC9lN1b//77b5nykcvlak8l13eTQGnyLU8ymaxCHrCXl5eH9PR0i69HU1m++5MmTcIrr7yilQd3b1lf0DNjxgx07twZK1asEKZZy59d9gwHPRXIWoIeZv2ULT2qTxQvTVDg7+8PV1dXnU8oL8vxWJ4nm5SUFLi6umL48OHllqc+rVu3RtOmTREXF2fxdanSrOuSkhKcOHHCpG7G77//HqdOncLff/+tNp1beqwv6FG+S1L1z0peXl5lFYfpwEFPBeKIn5lK2dLj6uoqTCtN68rt27cBABcuXACgHvSUpem9PE82mzZtAgD88ssv5ZanPqmpqQCAPXv2GJ03Li4OkyZNwt27d8u8Xs2gZ+vWrejVqxe6dOlich6aJ09zW3quXLmCXr16ITY21qzlrJm1BT266Hr/I6s8fE1PBTp58iQGDhxY2cWAvb09tzpZOV0tPU+ePCn1v3tlS0/t2rWFabpe8Guq8jzZqAZ2FcWUE1G7du2gUChw/fp1REdHl2l9xcXFauN79+4FAFy/ft3kPDSDXnOPhf79++PevXvo2LFjtXlXV1UIevi6HuvCLT0V6Msvvyzzj2d5qIyTDDOPMuixt7cXpj158sTk5TMyMrBw4UJhXBn0qD7Ftyz/QMvzmh7V4zEjI0MtTaFQmLXdpjJl25V/DJTXwJWF5p+M0nwHyxr0ZGVlmb1Oa2dtQY+up2RbW0tPdnY2/vjjD5v948tBTwU7dOiQzumHDx/Gxx9/XCEHYnUJembMmIHWrVtXyz5zXc/pMefkP3jwYISHhwvjuh5PUJYfY0u19DRu3FhoBQGAXr16wc3NzWAXU3R0tNB9Zypz6rI8ttVQ0PPTTz+Z1PWt+TR6c7u3StuqZ8107ZuHDx+afTyUF2dnZ61p1hb09OrVCwMGDMD69etLtTwR4dGjR+VcqorDQU8Fy87O1jm9X79++PTTT/H9999bvAyqXSZV2YYNG5CQkKD2fjZjpk6ditmzZ1uwVNpKSkpw7tw5s5q5dT2R2ZwT9eXLl9XGlS09qqwx6AGAuXPn4qOPPsLx48dx6tQpAEBUVJTOZWUyGXr06IEuXbqY1fpU0SciQ0HPuHHjsHXrVqN5/P7772rjqi09pnRXadazTCbDiBEj0KVLlyrbCrRmzRqtac8//zy6dOmCZcuW4fjx45g7d67wkmuliIgIfPPNN+VenooMevLy8jBjxgycOHFCmJabm4tvvvlGa3tVKX8btmzZYvY6U1JSUKNGDXh5eVXZu3856Klgjx8/NpheEXeVWEtLDxHh7bffxtKlS8uUj6kP7rt37x6+/fZbrFu3Ds2bN6+wvvbvvvsOL730Evr27WvyMrqCHnN+PDWfCWMo6Nm5cyfOnDmjNy9dJ9Ty7N7SDMLT0tLwxRdfoE+fPsI0fXc5qZbD0A+9psoOejS3+ciRIyblo7ovVFt6ioqKjM6v2dIzZcoU/Pzzz7hw4QKmTZtm0vqrgtzcXADA0qVL0adPH3z22WcYPXq0kF5SUoKwsDBMnz693B9fUBFBzw8//ICgoCDMmzcPGzZsQK9evYS06dOnY/r06Xj99deN5mPO90Xp7bffFj5PnTrV7OWtAQc9FUxfS4+SrqCovJ+JYi0tPUlJSdi+fTuWLVuGgwcPlrqbylhTq1QqxerVq9GoUSNh2q1btzBixIhSrc9cypYozVuODSlr91bdunXVxvUFPdevX0dISAi6d++uM5/CwkIEBQVp1VV5tvSY8tA+fUGP6jVPxr5bqio76NH842Hq09FVjwvVlh5dQfLRo0fRsGFDoYVIM+hRbT1LTk42af2GpKSkYMiQIYiJiSlzXuXtxIkTGDx4ML766iu16eYEPfv27dNqQdXk5OSkNU3Xn6uioiJ07969VMHmO++8g8TERGzYsEErTXn5hLFyAv8LDs2hGihV1dZBDnoMyM3NRc+ePbFx40aTlzHWN2+spUcz+j5z5gzc3NywcuVKvcsQEVJTU7Fq1SpMnDgRISEhOH36tNY8yqBC9Qf33LlzSEtL05lvSkoKLl26ZLC8muvQ3H4iwtmzZ5GZmQkiUvvnqfqv6M033zT5OS2aJxB9/1i+++47hISEQCwWY8GCBVrphw4dgp2dHezs7GBvb4/Tp09rBZgnT54s88XnDRs2NHuZn3/+WWuaJYIe1WBCVyBz8eJFJCUl4eeff1bbdxV9Aakpz7NRfreuXr1q9K4oU+pSNaAqq1GjRqnVn2bQk5SUZFKZVPejakuPrj9GgwYNQlZWlnDHqGbQoxo0ZWRk4KOPPhK6E3W5fPkyOnfujL/++gvAszvSTp06JQRiYWFhOHDgADp16qQ3j6SkJMyfP79cX4Fjaovtr7/+ipkzZ6pNM7W1Iz4+Hm+99Rbatm1rcD5TW3piY2Nx5swZbNq0qVyu41T+7jZo0KDMeRkiFostmn+FIEYSiYQAkEQiofj4eBo1ahQdPnyYxo8fTwAIACkUCpLL5VRSUkJERBkZGXT27Fl69OgR7d+/n3Jzc+nPP/8kANSvXz/Kz8+n+/fvC8urDgsWLKAzZ85QSUkJ7dq1i2bNmqU1T/v27cne3l5t2r///ksPHz6kvLw8On78ON2/f5/27dtHo0aN0rmeQ4cOUUlJCX366afCtNGjR1OrVq205j148CAdPHiQ3nvvPWrSpAn17dtXSDt8+DD99ttvlJ+fL9RZcXEx/fHHHxQZGUl///03ffDBBzrLoG9o3bo1rVmzRmu6h4cH/fLLLySTyejJkydUUlJCFy9epOzsbJLJZJSWlkYeHh5qddalSxciIiopKaH4+Hjq1KkTffHFF2aVR3WYOXOm1rTly5fTd999Rz/99BONHDmSkpKSiIjo5s2bdPPmTTp79ixFRUXR/fv3iYgoOzubdu/eTQUFBTR9+nQhH4lEQrm5uXqPRdV1/vfff2rHYFhYGBUXF9Pt27dp165d9PDhQ/r5559pw4YNFBERQRkZGUI+PXr0UMurV69eRER09epVYdpXX32lNn737l2t8iQkJAjpY8aMET7b29vTrl27SCKRGP1+FRYW0qeffkqxsbHCfvrrr78oMzOTiIj++OMPo/ukT58+OvPOyckR5tm7dy/l5+cL46rHq2b9NmjQwGi53dzchPmLi4t1zvPtt99Sr169SCqVCtP+/fdfndsQExMjzLN161atdBcXF53rUJ3n4sWLwvSPPvpImK6674mIkpOT1ZYjIhoxYoTauKenp85y6hMYGKg2z/LlywkADRkyhIiIWrRoYTQPV1dXAkCjRo3SO48xeXl5auVNS0tTSzd2LD19+lT4/M0336gte+fOHVq1ahU9fvxYmCaXy2nLli3CMosXLyaFQkFFRUV04MAB4TgmIurWrZvW+j799FOtbVDdP//9959Z269rmx48eEBERK+++qrRfWDKvtbnlVde0bl8cXExzZw5k5YvX05Ez77zSjk5OcJ5U9WDBw/ozp07tGLFCuG3sSJw0EP/C3p44KGihjfffJMmT55MdnZ2BEA4GVhyKCwspCNHjgjj4eHhdP36dWH88uXLWt+NmzdvGszz9ddfN/r92rZtmzA/EdHx48cJADk7O1NJSYlJZe/QoYPOvFWDni1btlBubq4wfuPGDa35lWkikchouRs2bCjMv2nTJp3zKNPfeecdYVpMTIzObdi3b58wzw8//KBzHoVCoXcdqnVIpB70qG7rgwcPdC4zceJEYby4uJj8/Px0lkGf5557Tm2eRo0aqY2rBtr6KNOdnZ31zmOMZtCjDKY116FvqFevnvC5R48east26NCBAFC7du0oIyODhg4dSrVr1zaYn7u7O+3YsYPCw8N1pnfp0oU++OADysvLowsXLlCnTp1o0aJFQvoXX3xBjx8/puLiYiE4SExMpBMnTtDmzZvpww8/pLi4OPrhhx/o559/1rmOUaNGUXZ2tlpQsnTpUtq/fz/t3r2bpk6dKnyPNJf95JNPSCwWU3BwMF26dIn+7//+jyZOnEijRo0iX19f6tSpk95tf+2110r9e6Q5cNBTgTjo4cEWhxo1aqgFW0ePHqWioiLKy8ujc+fOUZcuXXS2QmoOWVlZdPHiRYqKiqLevXvThQsX1L5fyhYBAFRSUkLLli0Txq9du2ZSWb29vWnnzp10+/Zttbznzp0rzLNmzRq17/KpU6e0vuuqeRqjepIfOHCgznlU85NKpUbrSyQS0YULF/QGPQ8fPjS4DtVyqwY9CQkJwvQrV67oXObDDz8UxpOTk9WCOtVB2UIWGxtLbdu2pT///JOI1P/lX7lyhRo3biyMFxUVqbU4P336lAoLCyk+Pp4KCwupsLCQFAqF2npGjx5Nu3btolu3btG+fftIJpNRQUEBKRQKSktLo6SkJMrLyyOiZy0JS5cupbVr19LRo0fV8hkxYgRdvHiRVqxYQQsXLizVd+GNN94gR0fHSv9O2vLAQU8pPHjwgPr160dOTk70/PPPC19WYzjo4YGHqjUsWbKEbt26pTX9jTfeUBufNWsWbdiwgZYtW0aTJ09WSzty5AiFh4fT+++/T5cvX6bNmzfTf//9R/n5+XTv3j0KCAhQmz8xMVE4eT9+/JgKCwvV0mfPnm1y+b/55hud0+Pj44noWfdAZGQkZWRkaM2zc+dOWrNmDTk4OAjTtm3bRg8fPqTHjx/TTz/9VOn7hwcezB0qKuixI6omzyMHMGLECIjFYnz11Vc4duwYJkyYgNTUVHh4eBhcTiqVVo8LtBhjjLEqSCKRVMgDNKtN0COTyeDp6Yn09HThbpmXX34ZEydOVHu2APDsTgfVux2kUil8fX0rtLzVmY+PD5o3b651BxljjDGmS0UFPdXmhaMpKSkQi8VqtwcHBwcjKSlJa97Vq1dj2bJlOvPp06cPFAqF2lMuAeDBgwdQKBRYtWoVXn75ZWRnZ8PNzQ2NGzdGfHw8evXqBalUip9//hl16tTBnDlz9JZ1+PDhqFOnDurXr48LFy7g2LFjaumtWrVCkyZN4OvrK7yBGnj2LA57e3s8fPgQBw4cQI8ePXD27FkUFxfj999/R48ePVBUVITffvsNDg4O6NGjBxYuXIjatWvj9u3b+Oijj/D777+rBXyNGzfG7du3MXjwYHh4eODVV19FcXExunbtimvXruH06dNo27Yt5HI5Nm7ciH///ReNGzfG2rVr8dNPPyE3Nxd//fUX0tLS0KJFC53vngGAxMREHD58GDNmzBBulSUizJ8/H59++ikAoFmzZujbty88PT3RoUMHREVFYefOnQCePQhrypQpiIyMxI8//liqB2sZ06tXL/z5558605ydnfH+++8LT4Bt06YN/P39sXv37nIvh9Kbb76JX3/9VWfaN998A39/f1y7dg3Tp09XS6tVq5beh9UxxqofOzu7Mr9E1tPTE3v37sXJkyexYsUKAM9+51599VWIxWKsWrVK7dwxevRoNGnSBCdPnkTbtm1x4sQJPHr0CAcOHEB0dDSysrLQsmVL3LhxA+PHj0erVq1Qo0YNXL16FYmJiTh+/DjS0tKwYMECvPbaa2Uqu1kqpBOtAkRHR1Pz5s3Vpi1YsICmTZumNW9BQQFJJBJhuHPnjlafYkFBAZ08eZIA/XeOGHL9+nVavXq1Vr/l33//rTWv5rUBqndwpKSkUNeuXem3334zuwz6KNczdOjQcsuzLDIzM2np0qV0584dk5dRfRxATEyM2i2RRUVFlJ+fT++++67RfuTnn39e5zGi6t69eySXy4Vxzdsv09LS6MGDB3Tq1CnhDoyEhAQqKioiIu2LZ7dv305r166lPXv2kFwup7/++os+++wzGjZsmNa8qrcB37p1i/bu3at1h49CoaBffvmFsrKyiOjZRZ+a2/niiy9Wep89DzzY4iCTyejs2bP03Xff0f79+7XuhlJ+nyUSCeXl5ZFMJqPExEQiIjpw4AB17NiRzp07p/e7DUD4XVKO9+nTh7Kysujjjz+m6Oho4ZEL6enpwjybNm2i9evX06VLl0z6zVW9UxIAjR071qTljFF9ZExFqDZBz6VLl8jLy0tt2vTp02nu3LlGlzVU6QkJCfTkyZNSlSkpKUnr4Dx58qTOeZXpX3/9danWZQ7lurp162bxdVmK6p0gv/76q975SkpKaOPGjTp/KLp3714hZdX8cdJn/vz5wnyrVq0iIqL27dubtKyhdQKg9PR0rVt9N2zYYPEf/AEDBqiNx8XFlettrjxYfhgwYIDeZ/ooB83ntyiHzz//nCIiItQuulYdtmzZQllZWXTq1Clq3bo1NWrUiL788ktKTEykkydP0qlTp+jPP/+k4uJiOn36NF28eFErj/Pnz5d62zIzM4WAID4+nh48eECTJk0iV1dX4YLwrl27qj3j6OjRozR58mS6evUqyeVyiomJofDwcLXHJagOuowdO9ZgujnfbdU8lOPKZ5dpksvlwjxr1qwp07oHDRpk9vK6cNBTSk+ePKHatWsLD4cjIurevTv9+OOPRpe1VKXruh139+7dOuddtGgRtWrVinJycsq1DLooy/L8889bfF2W9PLLL5O7u7tJdbZnzx4CQHXq1KElS5bQmjVr6NGjR5YvJBGJxWKTfuBWrlwpzKd82JhqgGAOzeMuPT2diIief/55tfzefvtttXGFQiHcslxSUkJLlizR+SOrfDhjTk6OUG4PDw9h/Y8ePRL+wRoK3Ex5MCEAmjt3LoWGhlKHDh3o2rVrJJfLtW6BvnbtGg0cOFAY/+WXX+jJkyeUlpZGUqmUPvvsM618z5w5QxkZGXTz5k1KSUmh9evX0//93//R3r17DZZHtR6tZdB3y7Wfnx+dOnVKGH/ttdeE74PmEBYWRnPmzKEpU6boPaEqpxkbV6VQKNRaKYYPH673gY/G9OvXT2td8fHxdPnyZcrKyjJaT6qtnroeEUBEQsuu8pb5efPmGf0e6nvmlC6qjxswl6F1qB6fxpb/6KOPyrRuzWcclRYHPWUwbNgwmjx5MuXl5dGvv/5KHh4elJ2dbXQ5S1V6SkqK1sG5cePGcl1HaSjL4u7uXtlFKROFQiH8KFkz5UPhvL29Dc735ZdfCvtGGfQMHTq0VD+OmsedMuhp166dWn6qz27RZfPmzTp/ZK9evSrMo1AoaPv27XTr1i2deQwePFjvOvT9O9YcVJ9mrG875XK5WhO8sntRaffu3ULa4cOHdT4IUEnzdvgZM2bQoEGDaPbs2TRz5kySy+VUUFAgpNerV48OHTpU5sBFWSZlF6quLtqOHTtSdHS01nRl96auJ4qr1lVwcLBW3Wnun7y8PJo0aRIdPnxYmPbVV18RAFq7di0RETk5Oakt16VLFwJA/v7+eut137591LVrV0pJSdE7jzGq3wnNY0pf949mPSk/m9qlrtx2Y99DU4Me1afkm2vnzp308ssvC8fbjBkztNZfp04do2UsTffUyy+/LCzfpk0bs5fXhYOeMnjw4AH17duXnJycqEWLFnT8+HGTlrNUpaelpakd/A0aNKiw1gVDlCc6awjAbIFMJqO1a9dSamqqwfm+++47raBH9TUU5tAX9PTs2VMtP+Uj4PXlr/oE544dO+oMeoxJTk4mOzs7mj59ukll1TWcOXNG57Kqdab8dx4XFyc870bV4cOHhXlPnz5tsMzZ2dlq61+8eLHO+RYsWEBNmjQRXkWg2tKkOURHR9P48eN1Pl9I3z5QbREAnr0ihYiE6xAB0MiRIykuLk5YZu3atWrLKP/RK8cbNWqkt94NKSkpof/++08IyGJjYyk4OJiOHj1KRM+uzfvoo490Pgm7PI0ePdpgmXVtV1BQEA0aNIjWr19PJSUl1LlzZ2rVqpXJrU2PHz8mDw8Po9dBmlqnERERpQ56VMlkMrVrDE3Zl8r0N954w+z1ZWZmCssfOnSoVGXWxEFPJbBUpas+WGzs2LEG/1lWpJKSEr3/ylnl2blzp1bQo/pOM3No/vAq97fqNTZEpHa9ky4KhYJGjx5NixcvVguQzAl6iJ61HOh6/w6R/ndAqQ5//PGHzmX37dunFfToc+bMGZODHs0WA13vJdNFVxfHsGHD6Pz582rzmRr0aD5hWFmHjx490gpqlHbs2KEWaCnfg6Sc5uDgoLcMVYHqHwFlMK9Kc5s6d+6s9sRqomf7ydzuNdV3L+pjap2qBr7lyZR9GRwcTADo+++/L9U6SkpKTOpBMVVFBz38lnUL0nxLc40a1lHddnZ2aNq0aWUXg2lwcXHRmlbah2aKRCK1ceUbvDXf7u3l5WUwnxo1aiAqKgrLli1Te6u3uZycnPQ+zkDXdmvKzc3VOd3Yg0dVmfMMENXv7tixY+Ht7W3Scrq2MSgoCJ07dza4XKtWrXD16lWt6Zpv7Vbm7+TkJExTvuVc6dVXXxWW7d69O2rVqqWWrrztWDUPAIiKijJYRmuh+nZ4T09PrfSMjAyhDgDg/PnzaNWqldo8dnZ2Wr/PxtSuXVvvMWyupk2bIj4+Hv/991+55KfUvXt3AEC7du30znP69GmcOnUKYWFhpVqHnZ2dWd87a2MdZ+FqqmbN/z0GiarHMyCZBWkGKkDpg56YmBiMHz9eGHd3dwegHfQMHDgQ7dq1w4QJE4zmWZagxxDNoKdLly4AgEmTJgnT9D081JwfX81tN1VeXl6pllMqKSkxOk9ERAReeOEFremaQY+S6olfM+hp2LAhMjIytE6oQ4YMAQAMGzYMAHD27FkhbfTo0Rg9erTRcloDY8ehr68v9u3bh2bNmmHGjBkVVCrzBQUFoVGjRuWa5y+//IIVK1bgwIEDeucRi8Xo0aOH1fwJr2jV5uGE1kj1n4QpP3zMtukKejp27FiqvFq2bIlt27Zh0qRJePToEfz8/HSuw8HBARcvXjQpT0sFPZplqlu3LvLz8+Hg4IB33nkHSUlJwj9YTaotlsb+uau29BQWFppcPolEYvK8uigUCoPpd+/e1duSpC/oUT1haQY9gO4g8aeffsJvv/2G/v37AwBefPFFIe3p06cGy2hNTDkOxWIxbt68WW4tM1WFl5cXFi5cWNnFsGoc9FiQ6o+wsR8+xnR183Tr1g27du1CixYtSpVn165d1cZL29oBPGvetwTN7c7NzRVaMjp16oROnTrpXdbNzQ2JiYmwt7c3GvSobruyu88U+rrWTFVcXGww3VDXmb6gR5WuoEcXkUiktzVH2bpWFai2chlSGQGPvb09/9ZbOQ56LEi1e4u/CMwYXS09ADBy5MhyW0dZgp6K6t76+OOPzVo+MDDQpPlUr22pyKBH13ffwcFB7ZH++mhed6OLqUGPLjdv3hRezlxVWOo4LA8XL17EihUr4OTkhMjIyMouDtPBNjv1KojqP09j//YYUz35W+p40RdYmaIiurfu3buHfv36WWQ9qqRSqdF5lF1IQUFBZVqXrn1pavBZni09ujRv3hzvvvuuVQcSmqy5rG3atMEvv/xi8EJiVrk46LEg7t5i5lA9+VvqGgtrb+lRfWGwJZkS/CUmJmLmzJlqL/01hWa3iq7vvqnBpylBT0FBgWkFqyZM7d6qTKbckcgqBwc9FsRBDzNH7dq10a1bNwQFBVnskQIvvfRSqZe11DU97733HgBUyJuWd+zYgfHjx2PMmDFG523ZsiXWr1+P+vXrm7UO1W5twPJBT1laeqoia27pUeKgx3rxNT0WpHqHBXdvMWPs7OwQHR0NIrLY7aSNGzfGtWvXUKdOHbOXtdTJpm3btsjKyipVmcw1ZswYkwKesrC3t0dRUZEwruu7/8ILLyAxMdFoXqpBz+nTp3XOw0GP9XnllVcA/O9REcx6cNBjQarN3NzSw0xhZ2dn8btO/P39S7WcJU829erVs1jeFc2Ulp4NGzaguLgYU6ZMMZiX6oXMyscOaLK1oKcqdG95e3vj3r17ZepOZpbBQU8F4aCHVXVV4R+2NdAMenS19NSvXx/79+83mpdqS4++B5zaWtBTVY7Diro+jZmHr+mpIBz0sKrOUtf0VDemtPSYSjXo0XfBsq0FPVWhpYdZLw56Kghf08OquqryD7uyaT4ksSzffdXuLX2vw7C1u7f4OGRlwUFPBeGWHlbV8cnGNOXZ0qOal62+K0kTH4esLPiangrCQQ+r6lRPNvwCXf1MuabHHJ9++inS09PRpk0btektW7bE1atX0aFDhzLlX9Vw9xYrCw56Kgh3b7GqTjXoMeeFnbZGs3urrH945s6dq3P60aNHsXnzZkybNq1M+Vc13NLDyqJc20uvX7+OAQMGoG7duvDy8kJoaChycnKE9PHjx8PBwQEikQgikUjrnTkRERHw8fGBm5sbwsLC1H5YU1NT8dJLL8HZ2Rlt27ZFXFyckFZSUoKZM2fC3d0d9evXx/r168tzs8oFt/Swqk71QmZT3htlq8qze8sQX19fhIeHo1GjRhbJ31pxiyMri3INeiQSCUaMGIHU1FSkp6ejsLAQs2fPVptn2bJlkMlkkMlkSEpKEqYnJCRg1qxZOHDgAO7cuYP09HSsXLlSSB89ejT69OmD7OxsTJgwAUOGDBFaTzZv3ozo6GjcuHED0dHRWLt2LU6ePFmem1ZmHPSwqo5bekxT3t1bTJ1q9xbXLTNXuQY9HTt2xNtvvw2xWAwXFxdMmjQJMTExJi0bFRWFkSNHon379hCLxVi0aBF27NgB4FkL0vXr1zF//nw4Ojpi+vTpUCgUOHfuHABg+/btmDdvHurVqwd/f39MmjRJWNZa8JeTVXWqF9JyS49+5d29xdSJxWLhc61atSqxJKwqsujtAOfOndPqwvrss8/g6emJrl27Ijo6WpienJys9jbj4OBgpKWlIT8/H8nJyfD391drXm/durXQUqRrWdVWJE1yuRxSqVRtsDT+4WPVCQc9+nFLj2W5uLjg999/x8GDB01+hxljSha7kPnKlSv4+uuv1QKbDz74AOvXr4eLiwt+/vlnDBw4EImJifD19YVMJoObm5swr/KzsitMNU2ZLpPJhHk0l1Wm6bJ69WosW7asXLbTVBz0sOrEUi9ErQ4q6poeW9a/f//KLgKrosxq6enTpw8cHR11DqrX36SlpWHgwIHYunWrWkvPiy++CA8PD9SuXRtjxoxBly5dcPz4cQDP3jqs2uKi/Ky86FmzNUYqlQpRvq5lDf0DmD9/PiQSiTDcuXPHnGooFf7hY9XBv//+i99++w0tW7as7KJYLW7pYcx6mdXSc+zYMaPzZGZmonfv3li0aBEGDx5scF7VawQCAgKQkJAgjMfFxaFp06ZwcnJCQEAArl+/jqKiIqEPNz4+HnPmzFFbNiAgQFhWs1tNlYODQ4Xf9sg/fKw6aNu2Ldq2bVvZxbBqfE0PY9ar3O/eev311/H2229j8uTJWum//PILnj59iuLiYuzevRtnzpzBq6++CgAICQnBnj17cOnSJUgkEoSHhyM0NBTAs7dC+/v7Y82aNZDL5di4cSPs7e3RtWtXAEBoaCg+++wzPHz4EDdu3MD333+PMWPGlOemlRn/8DFmG7ilhzHrVa5Bz4EDBxAfH4+1a9cK3VKq3Uzr16+Ht7c36tatiy+++AL79+9HkyZNAABBQUFYt24dBg4cCB8fH/j6+uKTTz4Rlo2KisKRI0fg7u6O7777Dvv27RN+XN59911069YNLVq0QLdu3TB79my89tpr5blpZcZBD2O2ga/pYcx62RE/3QlSqRRisRgSiUTrgumysrOzA/Ds4mqJRFKueTPGrM8bb7yBo0ePCuOtWrVS67pnjP2PJc+/uvAb7CoI/9tjzDZwSw9j1ouDngrCP3yM2QYOehizXhz0VBD+4WPMNvCFzIxZLw56Kgj/8DFmG/iWdcasFwc9FYSvF2fMNnBLD2PWi4MexhgrR3xND2PWi4MexhgrR5rdW9zSw5j14KCHMcbKEbf0MGa9OOhhjLFyxNf0MGa9OOhhjLFyxHdvMWa9OOhhjLFyxC09jFkvDnoYY6wcaQY9Q4cOraSSMMY01TQ+C2OMMVOpBj3ffvstRo0aVYmlYYyp4qCHMcbKkeo1PZMnT67EkjDGNHH3FmOMlSPN7i3GmPXgoIcxxsoRBz2MWS8OehhjrBxp3rLOGLMe5R702NnZwcXFBSKRCCKRCKtWrRLS8vPzERoaCldXV/j5+WHnzp1qy0ZERMDHxwdubm4ICwtDYWGhkJaamoqXXnoJzs7OaNu2LeLi4oS0kpISzJw5E+7u7qhfvz7Wr19f3pvFGGMm4ZYexqyXRVp6UlNTIZPJIJPJsGDBAmH6kiVLkJ2djbt372LXrl149913cePGDQBAQkICZs2ahQMHDuDOnTtIT0/HypUrhWVHjx6NPn36IDs7GxMmTMCQIUOE519s3rwZ0dHRuHHjBqKjo7F27VqcPHnSEptmtho1uDGNMVvCQQ9j1qtCz8jbt2/HkiVL4Obmhq5du2LQoEHYtWsXACAqKgojR45E+/btIRaLsWjRIuzYsQMAcP36dVy/fh3z58+Ho6Mjpk+fDoVCgXPnzgn5zps3D/Xq1YO/vz8mTZokLKuLXC6HVCpVGyyFfwAZsy3cvcWY9bJI0NO2bVs0atQI48ePx+PHjwEAOTk5yMzMRFBQkDBfcHAwkpKSAADJyclaaWlpacjPz0dycjL8/f1Ru3ZtIb1169YGl1Wm6bJ69WqIxWJh8PX1LZ8N14F/ABmzLfxHhzHrVe5BT3R0NG7fvo0rV64gLy8PEyZMAADIZDLY29vD2dlZmNfNzQ0ymUxId3NzU0tTTtdMM2VZZZou8+fPh0QiEYY7d+6Ucav146CHMdui+hvHGLMuZgU9ffr0gaOjo85Bef1N9+7dUatWLXh5eeHrr7/GoUOHUFhYCJFIBIVCgby8PCE/qVQKkUgEABCJRGrdTMrPyguiNbugjC2rTNPFwcEBbm5uaoOl8L8+xmzLqFGj0LVrV7XrGRlj1sGsM/KxY8fMylx5ES8RwcPDAw0aNEBCQgI6deoEAIiLi0NgYCAAICAgAAkJCcKycXFxaNq0KZycnBAQEIDr16+jqKgItWrVAgDEx8djzpw5assGBARo5VvZuKWHMdvi6OiIs2fPVnYxGGM6lGv3VlJSEuLi4qBQKJCTk4OZM2eid+/ecHBwAACEhoZixYoVePLkCS5cuICDBw9i5MiRAICQkBDs2bMHly5dgkQiQXh4OEJDQwEA/v7+8Pf3x5o1ayCXy7Fx40bY29uja9euQr6fffYZHj58iBs3buD777/HmDFjynPTSm3IkCEAgBYtWlRySRhjjDEbR+XoxIkT1KJFC3J2dqb69etTaGgoZWVlCel5eXkUEhJCLi4u5OPjQ5GRkWrLb9u2jby9vUkkEtG4ceOooKBASEtJSaGuXbuSo6MjtWnThi5fviykKRQK+uCDD0gsFpOXlxetW7fOrHJLJBICQBKJpHQbboBMJqMtW7bQvXv3yj1vxhhjrCqz5PlXFzsiosoOvCqbVCqFWCyGRCKx6PU9jDHGGPufij7/8pPzGGOMMWYTOOhhjDHGmE3goIcxxhhjNoEfIoNnt9QDsOjrKBhjjDGmTnnerajLiznoAYRXZVjydRSMMcYY0+3x48cQi8UWXw8HPQDq1KkDAMjIyKiQSmcVSyqVwtfXF3fu3OG786oh3r/VG+/f6k0ikcDPz084D1saBz3435OjxWIxf6mqMUu/coRVLt6/1Rvv3+pNeR62+HoqZC2MMcYYY5WMgx7GGGOM2QQOevDsretLliwR3hHGqhfev9Ub79/qjfdv9VbR+5dfQ8EYY4wxm8AtPYwxxhizCRz0MMYYY8wmcNDDGGOMMZvAQQ9jjDHGbAIHPYwxxhizCaUKeuRyOcLCwuDj4wOxWIyePXsiISFBSF+zZg28vLxQp04dzJ07V3iRWHFxMd566y00atQIdnZ2yMzMVMt31qxZaNasGVxdXdG+fXtER0cbLEdsbCyCg4Ph7OyMHj164Pbt20JaYGAgRCKRMNSoUQPr1q3Tmc/169cxYMAA1K1bF15eXggNDUVOTo7aPNu2bUOLFi3g4uKCli1b4saNG1i1apWQv4ODA2rVqiWMT506FQBw+PBhPPfcc3BxccGbb76plu/Dhw/Rv39/ODs7w9/fHydOnDC4vfrq1VhdVBeG6uvChQvo3LkzRCIRfHx8sGfPHr35GKqr/Px8hIaGwtXVFX5+fti5c6dFt4n9z5IlSxAQEIAaNWpg165dwvSIiAi0adMGrq6uaNasGTZv3mwwH96/1knf/pXL5XjnnXfg5eWFunXrYuzYsXj69KnefHj/WidDv8+Gzl2aLL5/qRRkMhktX76c7ty5Q8XFxbRu3Tpq1qwZERH98ccf5OfnR6mpqXTv3j1q2bIlbd26lYiIioqK6Msvv6Tz588TALp//75avkuXLqWUlBRSKBS0Z88e8vDwIKlUqrMMBQUF5OPjQ1u3bqX8/HyaO3cude/eXee8jx49otq1a1NKSorO9H/++Yd+/PFHys3NJZlMRsOHD6cJEyYI6QcPHqTAwEBKTEykkpISunnzJj1+/Fgtj9WrV9O4cePUpmVlZZG7uzsdOnSIZDIZjRkzhsaOHSukDx8+nCZOnEhPnz6l/fv3k4eHB2VnZ+sso6F6NacuqjJ99XXv3j3y9vam33//nYqKiujRo0d08+ZNnXkYq6s5c+ZQ3759SSKR0NmzZ0ksFtP169crahNt2vbt2+nYsWPUqVMn2rlzpzB98+bNdP78eSoqKqLExESqV68enT59WmcevH+tl779+9lnn1GHDh3o8ePHJJVKqXfv3rRw4UKdefD+tV76fp8Nnbs0VcT+LVXQo0kul5OdnR09evSIRo0aRWvWrBHStm7dSq+88or2inUEPZq8vb3p4sWLOtOOHDlCL7zwgjAuk8nIycmJ0tPTtebduHEjde7c2dTNoWPHjlGrVq2E8Q4dOtCxY8cMLqMr6Nm8eTO98cYbwvjNmzfJ0dGRCgoK6MmTJ1S7dm26d++ekN69e3f68ccfdeZvqF7NqYuqylB9zZkzhxYsWGBSPsbqqkGDBnThwgUhfezYsbRs2bJy2gpmih49eqidFDWFhITQ559/rjON96/109y/06dPp8WLFwvj33zzDQ0dOlTnsrx/rZOh32dTYwKiitm/5XJNz/nz51G/fn14enoiOTkZQUFBQlpwcDCSkpLMzjM9PR3Z2dl47rnndKZrrsfFxQXNmzdHcnKy1ryRkZEYM2aMyes+d+4cAgMDAQAKhQKXL19GQkICfHx80LRpUyxfvtxg85y+MjZv3hw1a9bErVu3kJKSArFYjIYNGwrpqnV15swZuLu7681LdV5z6qKqMlRfsbGxsLOzQ2BgIBo2bIixY8eqdSO2bt0aUVFRAAzXVU5ODjIzM8vl+GWWoVAoEBMTI3w/Ad6/Vd24ceMQHR2Nhw8fIjc3F7/88gt69+4tpPP+tX6Gfp+NxQQVvX/L/JZ1iUSCKVOmIDw8HAAgk8nU3oTr5uYGmUxmVp5FRUUYN24c5syZA7FYrHMezfXoW1d6ejpiYmKwb98+k9Z95coVfP3118L1RFlZWSguLsaJEyeQmJgIqVSKvn37wtfXF2FhYQbzkslk8PLy0lnGgoICneXPzc0FAHTr1k34rGt7VbfV1LqoyvRtY25uLu7evYvIyEgcPXoUjRo1wsSJEzFz5kz8+OOPAID4+Hij+chkMshkMtjb28PZ2VkrjVmHhQsXolGjRnj99deFabx/q7bmzZvD3d0d9evXh52dHV577TVMnDhRSOf9a/0M/T4biwkqev+WqaWnoKAAgwcPRv/+/TFhwgQAgEgkglQqFeaRSqUQiUQm50lEGD9+POrVq4elS5cK01UvTM7IyNBaj751RUVFoVevXqhXr57evJTS0tIwcOBAbN26Vfgn6eTkBACYN28e3N3d4efnh/feew+HDh0yui2Gymhq+fXlpTqvuXlVRYa20cnJCWFhYXj++efh4uKCRYsW6d0/xvaJQqFAXl6eVhqrfJs3b8a+ffuwd+9e2NnZ6ZyH92/VM23aNLi4uEAikSAnJwd169bFhx9+qHNe3r/WyZxznTnnOc18ymP/ljroKS4uxqhRo+Dt7Y3PP/9cmB4QEKB2J1dcXJxaU7QxM2bMwL1797Bjxw7UqPG/4iUlJQnRnp+fn9Z6nj59itTUVAQEBKjlFxUVpdW1pZkXAGRmZqJ3795YtGgRBg8eLMzr4eEBb29vteVN6doCtOvi1q1bKC4uRrNmzdCiRQtIJBK1O9gM1ZWhejW1LqoyQ/XVqlUrtXkN7R9DdeXh4YEGDRqU6fhllrF7926Eh4fj6NGjqFu3rt75eP9WPfHx8QgLC4Orqyvc3NwwYcIEnDx5Uue8vH+tk6HfZ3NiggrZv2ZdAaRi/Pjx1KdPHyosLFSb/vvvv1Pjxo3p1q1bdP/+fQoMDFS7UrugoIDy8/MJAKWnp1N+fr6QtnjxYmrTpg1JJBKj61de5b1t2zYqKCigjz/+WOuOpcuXL5OLiwvJZDKDeeXm5lLr1q31XhC1YMEC6t+/P0mlUrp79y61atWKtm3bpjaPobu3jhw5Qk+fPqXQ0FC1u7eGDRtGkydPpry8PPr1118N3r1lqF5NqYvqQF99HTt2jJo0aUKpqamUl5dHY8aM0doXSsbqavbs2cK+Pn/+PInFYrp27VoFbaFtKywspPz8fOrevTv99NNPlJ+fTwqFgo4ePUpeXl4UFxdnNA/ev9ZL3/59++236e2336anT58Kd7mOHDlSZx68f62Xvt9nYzGBqorYv6UKetLT0wkAOTo6kouLizBER0cTEdGqVavI09OT3N3dac6cOVRSUiIs27hxYwKgNgiFAcjBwUEtzx07dugtR0xMDAUFBZGjoyN1795d626l2bNnU0hIiNHtiYiIIABq63VxcRHS5XI5TZw4kdzc3Mjb25sWL16stk1EuoMeome3mjdr1oycnJxo4MCBakHNgwcPqG/fvuTk5EQtWrSg48ePC2nR0dFqZSAyXK/G6qI6MFRfX331FTVo0IDq1KlDo0ePVqvngIAAtePIUF3l5eVRSEgIubi4kI+PD0VGRlbMxjEaN26c1m/DX3/9RT179qSaNWuqfTenTJkiLMf7t2rQt3+zsrLorbfeIk9PT6pTpw4NGjSI7t69KyzH+7dqMPT7bOjcVdH7147IxL4axhhjjLEqjF9DwRhjjDGbwEEPY4wxxmwCBz2MMcYYswkc9DDGGGPMJnDQwxhjjDGbwEEPY4wxxmwCBz2MMcYYswkc9DDGGGPMJnDQwxhjjDGbwEEPY4wxxmwCBz2MMcYYswn/D/tGie9tCb4JAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "stream = obspy.read(filepaths_raw[0])\n", "stream.plot(size=(600, 250), show=True)\n", "plt.gcf().set_facecolor('w')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Preprocess\n", "\n", "The next cell loads every waveform and applies the preprocessing workflow, which includes detrending individual segments, filling gaps over the entire day, synchronizing the traces on the requested start and end times, resampling at the target sampling rate and removing the instrument gain. This cell runs in about 3 minutes on a desktop machine." ] }, { "cell_type": "code", "execution_count": 48, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "Processing data: 0%| | 0/24 [00:00" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAvIAAADtCAYAAADdhJErAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAABf20lEQVR4nO3dd3hTZf8/8HeA0pVOKKOUFsqmTGVIoQwfQQFBFEFlgwiIgyEgILKUoYCK+vAgigxtAQeoDyDjAWUWQfYSKLRQSlldabrbfH5/8Ov5Jm2SJl1J4P26rnNdyTl37vM5JyfJJ/e5z31UIiIgIiIiIiKHUsHWARARERERkfWYyBMREREROSAm8kREREREDoiJPBERERGRA2IiT0RERETkgJjIExERERE5ICbyREREREQOiIk8EREREZEDYiJPREREROSAmMgTERERETkgJvJEDiorKwsjR45EQEAAvLy80LVrV5w9e9agzOLFi+Hn5wdfX19MmzYNIgIAyM3NRf/+/VGrVi2oVCrcvn3b4HWTJ09GcHAwPDw80KZNG+zfv99sLMeOHUPLli3h5uaGLl264Pr168qykJAQqNVqZapQoQKWLVtmtJ5Lly7h2WefRdWqVeHn54chQ4YgKSnJoMyaNWvQoEEDuLu7o0mTJrh8+TIWLlyo1O/s7AwnJyfl+bhx4wAAv//+O+rXrw93d3c899xzBvXeu3cPvXv3hpubGxo1aoQ9e/aY3V5T+7WofVFQRkYGhgwZAg8PDwQGBmLDhg0Gy9euXYuAgAB4enpi5MiRyM7ONhtXcURGRqJbt27w8fFBzZo18eabbxqsx9y+2b9/Pzp37gx3d3c888wzBvVa8l4WZGq/HjhwwOAYcnNzQ4UKFXDv3j2j9axduxatWrWCh4cHgoODsXLlSoPlGo0Gr776Knx9feHt7Y1BgwYBMDxWK1SoAFdXV+X5gQMH8NFHH6Fx48bw8PBA06ZNsXnzZqXOa9euoW3btvDx8YGvry/69euH+Ph4k9vatWtXuLi4KPX37NlTWZaeno6xY8eiWrVqqF69OpYuXWp2v6lUKri7uyt1LVy4UFkWExODp59+Gl5eXmjUqBF2795tsp779+9j4MCB8PX1RWBgIMLDww2Wx8fH48UXX4SXlxeqVKmCadOmmY2LiMqBEJFD0mq1Mn/+fImNjZXc3FxZtmyZBAcHK8u3bdsmgYGBcvXqVbl165Y0adJEVq9eLSIiOTk58tlnn0lkZKQAkPj4eIO6586dK1euXJG8vDz54YcfxMfHRzQajdE4MjMzJSAgQFavXi0ZGRkybdo0CQsLM1r2/v37UrlyZbly5YrR5X/99ZesW7dOkpOTRavVyoABA2TUqFHK8t9++01CQkLk3LlzotPpJCoqShISEgzqWLRokQwfPtxg3p07d8Tb21u2b98uWq1WBg8eLEOHDlWWDxgwQEaPHi1paWmyZcsW8fHxkcTERKMxmtuv1uwLEZGpU6dKz549JSUlRQ4dOiReXl5y6dIlERE5c+aM+Pj4yLFjxyQ5OVm6du0q77//vsm6iuv333+XLVu2iFarlfv370tYWJjMnz9fWW5u3xw7dkzCw8Plgw8+kKefftqg3qLey4LM7deCli9fLp06dTJZ18qVKyUyMlJycnLk3LlzUq1aNdm3b5+y/IUXXpC3335bkpOTJTs7W06cOFGojkaNGskff/xhMG/JkiVy5swZycvLk3379omnp6dERUWJiIhGo5Fr166JTqeTrKwseffdd2XgwIEmY+zSpYts2LDB6LKZM2dK9+7dJSUlRW7cuCH169eXHTt2mKzL2Gc4X2hoqLz33nuSk5Mj+/btE19fX7l3757RsoMGDZJhw4ZJZmamst/Onz8vIiI6nU7atGkjH3/8sWi1WsnIyJDTp0+bjImIygcTeaKHRFZWlqhUKrl//76IiLz88suyePFiZfnq1aulW7duhV5nLgnI5+/vL3///bfRZTt27JDGjRsrz7Varbi6ukpMTEyhsitWrJAnnnjCou0REdm1a5c0a9ZMed62bVvZtWuX2dcYS+RXrlwpzzzzjPI8KipKXFxcJDMzU1JTU6Vy5cpy69YtZXlYWJisW7fOaP3m9qs1+0JEpEaNGnLkyBHl+dChQ2XevHkiIjJ9+nQZN26csmzPnj1St25ds9teGlatWiXPPvusiIjF+2bNmjWFEvmCCr6XBVl6vIqItGvXTlauXFnktuQbNGiQLF26VEREzp07J0FBQZKbm2v2NcYS+YJCQ0Plp59+KjQ/KytLZsyYIW3btjX5WnOJ/GOPPSbbt29Xni9cuFAGDRpksi5Tn2GNRiMqlUrS0tKUeT169JBVq1YZrcfX11cuXLigPB8zZozMnDlTRB780QoNDTUZAxHZBrvWED0kIiMjUb16dVSpUgUAcOHCBTRv3lxZ3rJlS5w/f97qemNiYpCYmIj69esbXV5wPe7u7qhXrx4uXLhQqGx4eDgGDx5s8boPHz6MkJAQAEBeXh5OnjyJs2fPIiAgAHXr1sX8+fMNurWYUjDGevXqoVKlSrh27RquXLkCLy8v1KxZU1muv68OHjwIb29vk3Xply1qXyxevBjPPvssACApKQm3b9+2uK6WLVsiOjoaGRkZRW5vSejv86L2TXHrNcbS4zUqKgqnTp3CgAEDLFpvXl4ejh49qqz777//RsOGDTFkyBBUqVIF7dq1w4EDB6zcGiA1NRXnz59H06ZNDeZ7e3vD1dUVS5cuxTvvvKPMHz9+PMaPH29Q9q233oKfnx+6d++OM2fOGCzTP65FxGBftGjRAhEREQblH3vsMdSqVQsjRoxAQkKCxXUVZKrssWPHULduXTz99NOoWrUqunXrVqzjgIhKFxN5oodASkoKxo4diwULFijztFotPD09leeenp7QarVW1ZuTk4Phw4dj6tSp8PLyMlqm4HpMrSsmJgZHjx7FwIEDLVr3qVOn8Pnnn+P9998HANy5cwe5ubnYs2cPzp07h3379mHTpk1Yu3ZtkXWZi7Go+Dt16oTk5GSTdemXLaqu6dOnY+vWrUrZihUrws3NzaK68h9b+x5a4/fff8fvv/+OSZMmWbQ9lir4Xhpj6fEaHh6OZ555Br6+vhate9asWahVqxaefvppAEBcXBx2796Np556Crdv38b06dPRr18/JCYmWrVNY8eORd++fdGkSROD+cnJyUhKSsKiRYtQp04dZf6KFSuwYsUK5fnHH3+M6Oho3LhxA927d0evXr2U7e3Rowc+++wzJCcnIyYmBmvXrkV6erry2jNnzij9+oEH1ypcv34dp06dQnp6OkaNGgUA8PDwQPv27bFw4UJkZ2fjzz//xL59+wzq0tejRw989NFHyMjIwNmzZ7F582albFxcHDZt2oSJEyfi1q1beP755/Hcc88hJyfHqv1GRKWLiTyRg8vMzES/fv3Qu3dv5QccANRqNTQajfJco9FArVZbXK+IYMSIEahWrRrmzp2rzNe/IPDGjRuF1mNqXREREXjqqadQrVo1k3Xli46ORp8+fbB69WqlJdXV1RUA8O6778Lb2xuBgYF44403sH379iK3xVyMlsZvqi79stbUpVarkZeXZ5BUmasr/3FR76H+haH6F1AW5dixYxg5ciR++eUXVK9e3ertMcXYe2ksRkuP14iICIOzOua2d+XKldi8eTN++uknqFQqAA+Oo7p16+LVV1+Fk5MTXnjhBdSrVw+RkZEWb9P06dNx48YNfPXVV0aXe3p6Yvjw4Xj++edNnjFq164d1Go1XF1dMW3aNKjVahw9ehTAgz8fderUQZMmTdCzZ08MGDAAtWrVMhlPWFgYnJyc4Ofnh88//xzbt29XLlgODw/H33//DX9/f3z44Yd48cUXTdb1+eefIz09HUFBQRg1ahReeeUVpayrqyvCwsLQs2dPVK5cGW+//TY0Gg0uX75s8X4jotLHRJ7IgeXm5uLll1+Gv79/oZEtmjZtajCKzenTp812bSjorbfewq1bt/D999+jQoX/+6o4f/680pIdGBhYaD1paWm4evVqoS4HBRMwY3UBwO3bt9G9e3e8//776Nevn1LWx8cH/v7+Bq+3pFsNUHhfXLt2Dbm5uQgODkaDBg2QkpJiMHKPuX1lbr9aui/yt6dGjRoW13X69GnUrVtX+UNjSlhYmLJPf//9d7Nl8128eBF9+/bFt99+iyeeeEKZb+2+KcjUe2ksRkuO12PHjiE+Ph59+vQpcns3bdqEBQsWYOfOnahataoyv1mzZoXitPQ4AoAlS5bgv//9L/773/+afS90Oh3i4+NNtn4XpP8Zc3d3x6pVqxAfH4+LFy9CpVKhTZs2VtWTv03BwcHYuXMn7t+/j//973+4fv26ybr8/Pzw448/4u7duzh27BiSkpKUssb2GxHZAdt0zSei0jBixAjp0aOHZGdnF1q2detWCQoKkmvXrkl8fLyEhIQYjAKSmZkpGRkZAkBiYmIkIyNDWTZ79mxp1aqVpKSkFBlD/kgta9askczMTJk+fXqhkVpOnjwp7u7uotVqzdaVnJwsLVq0UC74LGjmzJnSu3dv0Wg0EhcXJ82aNZM1a9YYlDE3as2OHTskLS1NhgwZYjBqzYsvvihjxoyR9PR0+fXXX82OWmNuv1qyL/RNmTJF2Z7IyEjx8vKSf/75R0QejFrj6+srx48fl+TkZHnyySfLZNSaGzduSO3atWXt2rVGl5vbN3l5eZKRkSGrVq2S7t27S0ZGhnIsFvVeFlTU8SoiMmHChELvrTE7d+4UPz8/o6OqZGdnS506dWTt2rWSm5srv/76q1SpUqXQ6EfGLnZdvXq1BAYGSmxsbKF69+3bJ8ePH5fc3FxJTEyUoUOHSrt27YzGl5SUJLt27ZLMzEzJysqSTz75RKpXr6583mJjYyU+Pl5yc3Nl165d4u/vL3FxcUbrOnfunJw6dUpZ7yuvvCI9e/ZUll+4cEG0Wq2kp6fLZ599Jq1atRKdTme0rqioKElMTJTs7GzZsGGDBAYGKhfK3r9/X3x9fWXXrl2Sm5srX375pdSrV8/odw8RlR8m8kQOKiYmRgCIi4uLuLu7K9P+/fuVMgsXLpQqVaqIt7e3TJ061eAHPCgoSAAYTPkAiLOzs0G933//vclYjh49Ks2bNxcXFxcJCwsrNErLlClTzI66kW/t2rUCwGC97u7uyvKsrCwZPXq0eHp6ir+/v8yePbtQUmIskRd5MOpGcHCwuLq6Sp8+fQwS9bt370rPnj3F1dVVGjRoILt371aW7d+/3yAGEfP71dy+WLBggcHoOenp6TJo0CBxd3eXgIAACQ8PN1jPmjVrxN/fX9RqtQwfPlwyMzOL3IfWmjt3rqhUKoP93bRpU2W5uX3zxx9/FDqG8vd9Ue+lMeb2a25urtSoUaPIUYtERLp27SqVKlUyWO/YsWOV5adPn5Y2bdqIu7u7tGrVymBoynzGEvk6deqIk5OTQb0LFiwQkQfHV+PGjcXd3V38/PzkxRdfNHjvx44dq8Rw9+5defzxx8Xd3V18fHykW7ducvz4caXs3r17JSAgQFxdXaV169YGn2kRkaZNmyqfxz179kiDBg3Ezc1NqlevLkOGDJE7d+4oZZcsWSK+vr6iVqvl2WeflRs3bijLCh7b4eHhUq1aNXFzc5POnTvL2bNnDda7d+9eady4sajVaunUqVOh5URU/lQiVpxTJCIiIiIiu8A+8kREREREDoiJPBERERGRA2IiT0RERETkgJjIExERERE5ICbyREREREQOiIk8EREREZEDYiJPREREROSAmMgTERERETkgJvJERERERA6IiTwRERERkQNiIk9ERERE5ICYyBMREREROSAm8kREREREDoiJPBERERGRA2IiT0RERETkgJjIExERERE5ICbyREREREQOiIk8EREREZEDYiJPREREROSAmMgTERERETkgJvL/3/79+9GnTx/4+/tDpVLhl19+KfN1xsXFYciQIahSpQrc3NzQqlUrHD9+vMzXS0RERESOj4n8/5eWloaWLVviyy+/LJf1JSUloWPHjnBycsLvv/+OCxcuYNmyZfD29i6X9RMRERGRY1OJiNg6CHujUqmwZcsW9OvXT5mXnZ2NWbNmITw8HMnJyWjWrBk++ugjdO3atVjrmD59Og4dOoQDBw6UTtBERERE9Ehhi7yFRo4ciUOHDmHjxo04c+YMBgwYgGeeeQZXrlwpVn2//fYb2rRpgwEDBqBatWpo3bo1vv7661KOmoiIiIgeVmyRN6Jgi/zVq1fRoEED3Lx5E/7+/kq5p556Cu3atcPChQutXoeLiwsAYPLkyRgwYACOHj2KiRMn4quvvsKwYcNKZTuIiIiI6OFVydYBOIITJ05ARNCwYUOD+VlZWahSpQoAICYmBnXr1jVbzxtvvKH0wdfpdGjTpo3yJ6B169Y4f/48/vOf/zCRJyIiIqIiMZG3gE6nQ8WKFXH8+HFUrFjRYJlarQYA1KpVCxcvXjRbj4+Pj/K4Zs2aaNq0qcHyJk2a4Oeffy6lqImIiIjoYcZE3gKtW7dGXl4e7t69i7CwMKNlnJyc0LhxY4vr7NixIy5dumQw7/LlywgKCipRrERERET0aGAi//9ptVpERUUpz6Ojo3Hq1Cn4+vqiYcOGGDx4MIYNG4Zly5ahdevWuH//Pvbu3YvmzZujV69eVq9v0qRJCA0NxcKFCzFw4EAcPXoUq1atwqpVq0pzs4iIiIjoIcWLXf+/P//8E926dSs0f/jw4Vi7di1ycnLw4YcfYv369YiLi0OVKlXQoUMHzJs3D82bNy/WOrdu3YoZM2bgypUrqFu3LiZPnozXXnutpJtCRERERI8AJvJERERERA6I48gTERERETmgR7qPvE6nw61bt+Dh4QGVSmXrcIiIiIjoESciSE1Nhb+/PypUMN/m/kgn8rdu3ULt2rVtHQYRERERkYHY2FgEBASYLfNIJ/IeHh4AHuwoT09PG0dDRERERI86jUaD2rVrK3mqOY90Ip/fncbT05OJPBERERHZDUu6fT/Siby18vLycPr0acTFxUGn09k6nIeCSqWCr68v2rRpAxcXF1uHQ0REROQw7DqRj4yMRMeOHbFw4UJMnz690PIRI0Zgw4YNcHJyAgAEBQXh/PnzZRLLhQsXMGnSJCQkJJRJ/Y86V1dXzJw5Ez179rR1KEREREQOwW4TeZ1Oh0mTJqFt27Zmy82bN89okl+aMjIy8OabbyIoKAiffPIJGjdujIoVK5bpOh8VOp0ON2/exKpVqzBnzhzUq1cPDRs2tHVYRERERHbPbseRX7VqFdq3b48mTZqUWp1ZWVnQaDQGkyUOHz4MjUaDDz74ACEhIaWaxN+5cwfnz59HTk5OqdXpSCpUqIDAwEDMmzcPHh4e2LVrl61DIiIiInIIdpnIJyYm4rPPPsPcuXOLLLtkyRJUqVIFoaGh2L9/v9myixYtgpeXlzJZOvRkdHQ0fH19ixwCqDhiY2ORkZGBW7dulXrdjqRixYoICQlBTEyMrUMhIiIicgh2mcjPnDkTEydOhI+Pj9lyEyZMQFRUFOLj4/HGG2+gT58+iI2NNVl+xowZSElJUSZzZfXl5eWhUqWy7YXEi2cBJycn5Obm2joMIiIiIodgd4n8yZMncfToUbz22mtFlm3dujV8fHxQuXJlDB48GB06dMDu3btNlnd2dlaGmizpkJMvvvgiJk2aZDBv6dKl6NKlC4KCgnDkyBGDZePGjTN5hmHs2LFo3LgxPDw84O3tjU6dOmHt2rWFyn399ddo3rw53N3dERgYiOHDhyst2BkZGRgyZAg8PDwQGBiIDRs2GLz2n3/+Qffu3eHh4YHq1atj+fLlFm3n3Llz4eTkBLVaDQ8PDzz++OM4ePCgwfJx48YZvObIkSOoU6cOACA8PBxqtdpgcnNzg0qlKvIMChERERGZZneJ/L59+3D58mXUqlULNWrUwKZNm7BgwQKLEvuibmNbmr744gt89913OHHiBADgxo0bWLRoEb766iuLxv0s6OOPP0ZqaiquX7+OiRMnYubMmQYX8X744YeYPXs2PvroIyQkJODChQvo2LEj9u7dCwCYM2cOEhMTERcXh40bN+L111/H5cuXATxI8nv27Ilhw4bh/v37iIqKQo8ePSyO7dVXX4VWq4VGo8Hrr7+OgQMHWvzawYMHQ6vVGkyjRo1CWFgYOnbsaHE9RERERGTI7hL5MWPGICoqCqdOncKpU6fQt29fTJgwAUuWLClU9ueff0ZaWhpyc3OxadMmHDx4EE8++WS5xFmzZk18+OGHGDNmDPLy8vDWW2/h7bffRuPGjUtUr5eXF1588UWsXLkSn3zyCe7du4fk5GQsXLgQ//nPf9CrVy+4uLhArVZjzJgxGDVqFADgu+++w5w5c+Dp6YnQ0FD07dsXGzduBACsWbMGoaGhGDp0KJydneHh4VGsi4hVKhUGDRqE+Ph4pKSkFGv7Nm/ejE2bNiEiIoIj/xARERGVgN0l8m5ubqhRo4Yyubq6Qq1Ww9vbG+Hh4QgJCVHKfvrpp/D390fVqlXxySefYMuWLUqXjvIwduxYuLi44IUXXsCVK1cwY8aMYtclIgbPe/XqBZ1Oh2PHjiEyMhLZ2dl49tlnjb42KSkJt2/fRvPmzZV5LVu2VMbUP3bsGHx9fdG+fXv4+fmhb9++Fl8foC8vLw/h4eFo1KgRvLy8rH59TEwMRo8ejbVr15bJhcNEREREjxK7HUc+n35f8cGDB2Pw4MHKc/2+2ragUqmwatUqhISEYO/evahcubKyrHv37gYtzhkZGVYl+pUqVULVqlWRlJQEEUHVqlVNXnCr1WpRsWJFuLm5KfM8PT2h1WoBAHFxcTh69Ch2796NFi1aYNasWRgyZAj27dtnUSxr1qzBxo0bkZWVBRFBRESE0eX58vLyUKVKFYMyOTk5ePnllzFq1Cj07t3bovUSERERkWl21yLvaJo2bQoAhbqq7N69G8nJyco0cuRIs/UUbJHPzc3F/fv34ePjgypVquD+/fsmR3RRq9XIy8tDenq6Mk+j0UCtVgN4cNfUF154AW3btoWzszPmzp2LAwcOKIl+UUaOHInk5GSkp6fjjz/+wKhRo3Dy5MlCy/MnYxccz5w5EyKCRYsWWbROIiIiIjKPibyd+v3331GhQgW0a9cOHTp0gJOTE7Zt22a0rI+PD2rUqIGzZ88q806fPq10Q2rWrJnR1xX881AUlUqFDh06oEGDBtizZ4/Fr9u+fTtWr16NjRs3wsnJyap1EhEREZFxTOTtjEajwebNmzFu3DhMnjwZVatWhbe3N9577z2MHz8eO3bsQFZWFtLT0/HNN9/g22+/BQAMGTIEH3zwAVJTU3HkyBH89ttveOmll5Rlv/32G06ePImcnBx88MEH6Ny5Mzw8PKyO79ixY7hw4YLFF8vGxcVh+PDhWLVqFerWrWv1+oiIiIjIOCbyNlDwol0AmD59OtRqNWrXro1ly5bhww8/xOLFi5Xls2bNwpw5czB16lT4+PigUaNGOHDgAP71r38BAObPnw8vLy/UrFkTAwYMwIoVK9CoUSMAD7r9/Pvf/8bzzz8PPz8/nD17Ft99951Sd0hICMLDwwE8GEZTrVbjxo0byvLVq1dDrVbD3d0dAwYMwOzZsy3u57569Wrcv38fI0aMKDSefP46iYiIiMh6KrG2f8VDRKPRwMvLCykpKWZvDvXVV1/h119/xfbt20s9hr///hvAg37uJR260tG98847yMvLw2effWbrUIiIiIhswtL8FGCLPBERERGRQ2Iibyce4RMjCu4DIiIiIssxkbeAu7s7NBqNyeEfqXQkJSXB3d3d1mEQEREROQQm8hZo164dMjMzjY6PXloe9dboq1ev4ty5c2jfvr2tQyEiIiJyCHZ/Z1d70KBBAzz55JOYN28eDhw4gEaNGpm8y6q1rl27BgCoXLmywU2WHhV5eXm4efMmdu3ahXr16imj8BARERGReRy1xsKrgnNzc/H999/jf//7H27duoW8vLxSieHEiRMAHtx91dKx2R8mKpUKvr6+6NKlC4YPHw5vb29bh0RERERkM9bkp0zkLdxRZUWlUgF4MNb7hQsXbBIDEREREdkHDj/pgEqrhZ+IiIiIHg1M5O0EE3kiIiIisobdJvJjxoxBzZo14enpiebNm2Pr1q1Gy2VkZGDIkCHw8PBAYGAgNmzYUM6Rlg4m8kRE9qM0r4UiIiordpvIT548GTExMdBoNPj2228xZMgQJCUlFSo3Z84cJCYmIi4uDhs3bsTrr7+Oy5cv2yDikuEPBhGRfdi7dy9q1aqF3r172zoUIiKz7DaRb9y4MZydnQE8uCA0MzMT8fHxhcp99913mDNnDjw9PREaGoq+ffti48aNRuvMysqCRqMxmOwFE3kiIvvw5ZdfAgB27txp40iIiMyz20QeAMaPHw9XV1e0bdsWzzzzDJo2bWqwPCkpCbdv30bz5s2VeS1btsT58+eN1rdo0SJ4eXkpU+3atcs0fmswkScisg+ldZ8QIqKyZteJ/IoVK6DVarF792506dKl0HKtVouKFSvCzc1Nmefp6QmtVmu0vhkzZiAlJUWZYmNjyyx2azGRJyKyD05OTrYOgYjIInadyANAxYoV8dRTT2HPnj2FTnOq1Wrk5eUhPT1dmafRaKBWq43W5ezsDE9PT4PJXjCRJyKyD2yRJyJHYfeJfD6dToerV68azPPx8UGNGjVw9uxZZd7p06cREhJS3uGVGBN5IiL7wESeiByFXSbyWq0W4eHh0Gq1yM3Nxc8//4w//vgDYWFhhcoOGTIEH3zwAVJTU3HkyBH89ttveOmll2wQdckwkScisg/6XWuGDx+OW7du2TAaIiLT7DKRV6lUWL16NQICAlClShUsXLgQERERaN68OcLDww1a3OfPnw8vLy/UrFkTAwYMwIoVK9CoUSMbRl88TOSJiOyDfov8+vXrMWzYMBtGQ0RkmkpExNZB2IpGo4GXlxdSUlJs1l9epVIBeNAClJ2dbZMYiIjo/0ycOBHLly9XnlevXh23b9+2YURE9CixJj+1yxb5RxFb5ImI7EPBUWtycnJsFAkRkXlM5O2ETqezdQhERITCF7vaSyJvL3EQkf1gIk9ERKSnYCKfm5tro0j+z3vvvQdXV1eDUdqIiJjIExER6bHHrjULFy5EXl4epk+fbutQiMiOMJEnIiLSY48t8vk4KAIR6WMiT0REpMeebwhlD2cHiMh+MJEnIiLSU7BrjT1hizwR6WMiT0REpMeeW+SZyBORPibyREREeuw5kWfXGiLSx0SeiIhID7vWEJGjYCJPRESkx55b5JnIE5E+JvJERER67DmRZ9caItLHRJ6IiEhPxYoVbR2CSWyRJyJ9dpnIz5kzB02bNkWFChWwceNGk+VGjBgBZ2dnqNVqqNVqhISElGOURET0MFKpVLYOwSQm8kSkzy4T+QYNGmD58uVo165dkWXnzZsHrVYLrVaL8+fPl0N0REREtsGuNUSkzy4T+SFDhqB79+5wcXEp1XqzsrKg0WgMJiIiIkfBFnki0meXibw1lixZgipVqiA0NBT79+83W3bRokXw8vJSptq1a5dTlKbZ8ylcIiKyL0zkiUifQyfyEyZMQFRUFOLj4/HGG2+gT58+iI2NNVl+xowZSElJUSZzZcuLPV9URURlY/LkyXjvvfdsHQYR/X95eXkYM2YMvvnmG1uHQmQVh07kW7duDR8fH1SuXBmDBw9Ghw4dsHv3bpPlnZ2d4enpaTDZGhN5okfLnTt38Omnn2LhwoVIS0uzdThEBGD79u34+uuv8dprr9k6lDL18ccf46effrJ1GFSK7Hew3GKoUMHx/pcwkSd6tOh3p0tLS4O7u7sNoyEiAEhPT7d1CGXu+PHjePfddwEAImLjaKi02GXmm5OTg8zMTOh0OoPHBf38889IS0tDbm4uNm3ahIMHD+LJJ5+0QcTFx0Se6NHi5OSkPM7IyLBhJESUz9XV1dYhlJnw8HA89thjOHHihK1DoTJgl4n8a6+9BldXVxw4cADDhg2Dq6sr9u/fj/DwcIOx4j/99FP4+/ujatWq+OSTT7BlyxbUqVPHdoEXAxN5okeLfos8E/mHQ3Z2Nn766Sfcv3/f1qGYZaxBjB5wc3NTHufl5dkwktI3ZMgQnDx5EhMmTFDm2XOLvIjg0qVLyM3NtXUoDsEuE/m1a9dCRAymrl27YvDgwQZjxR88eBApKSlITk7GX3/9hX/96182jLp4mMgTWS4zMxMbN260+4TJUkzkHw4LFy7EgAED0KlTJ1uHYtLmzZvh7e2NrVu32joUu6TfIv+wfi71t8ue70ewZs0aNG7cGIMGDbL4NTqdDtHR0WUYlf2yy0T+UcJEnshyc+bMwSuvvIIuXbqUy/qio6MRHByM5cuXl0n9D2vC8KjZvHkzAODSpUs2jsS0/v37IzU1FX369LF1KHZJ/741j0J/+czMTFuHYNLSpUsBAD/++KPFrxk5ciSCg4Oxbt26sgrLbjGRtzEm8kSGRASrVq3CyZMnCy375ZdfAAAXLlwol1hmzZqF6OhoTJw4sUzqZyL/cFCr1SWu48iRI7hz547J5bznSNnSHyzjUfhcZmVl2ToEk/S7OVlq/fr1AB409jxqmMjbGBN5IkO//PILxo4di8cee6zQstJImKxR1p9PW7SKXbx4ESNHjsTVq1fLfd0Pq5KOPBQZGYkOHTqgRo0aJstUrly5ROsgy7FF3raKk8jnexSH9GUib2NM5B8dD9sFVGXln3/+MbmsNBP569evY82aNWb7ipb1HwdbtPx1794da9euxVNPPVXu6zYlf4QyR1XSRH7fvn1FltEf7ag4SpIcPWrYIv9/4uLiyj3pZyJvHSbyNqafyNvzVeRUMmlpaQgODsbLL79s61DsnrmkqDQT69atW2PUqFH4+OOPTX72zK2vND6vtkgY4uLiAAAxMTHlvm5TunXrhho1ath1S6i597ukibwlry9pi7ypY3nLli0IDAzErFmzsGTJEv4OwfoW+StXruDzzz+36+4qBVkSa1RUFAICAtC4ceNyiOj/lOTzVNrfqY7weWAib2P6ifzDMDTYyZMn8e233zrEwV+edu7ciRs3bmDTpk22DsXulVcin5SUBACIiIhA48aNMXjwYIvXd/bsWVSrVg2ff/55iWKwdcvfxo0bbbr+fPv370diYqJFLdO2Yu6MmiMn8i+88AJiY2OxYMECTJs2Dbt27SrReh4G1n4umzRpggkTJmDBggVlFFHps6SVfefOnQAenL0sLf/88w/WrVsHEUFKSgqmT5+OU6dOGZSxl7NH//zzD2rVqlXi7/myxkTexvQT+Yeh68Vjjz2GV199FVu2bLF1KHbF09NTefww/GErS/pJTcF9VRZdXS5cuIDLly8jIiKi0DL99en/OZ0yZQru379vMC5zcdgikdffv6+88oryY20P7Pm0uLmuPwUT8YyMDKtadfVfb6oRpKRdayz97OgP8fyosvS9u3v3LsaNG6f8du/YsaMswypVlrTIe3t7l/p6mzRpghEjRmD9+vWYN28ePvroI7Ru3dqgjL0k8tOmTUN8fHyJv+fLGhN5G9NP5AcOHIhVq1bZMJrSs3//fluHYFe8vLyUxxqNxoaR2D/9pKZgq5GHh0eZrjs7O9vguX7yox+LtT80N27cwIgRIwqNxGOLRF7/WASAw4cPGzw/duyYVWeOpk2bhqFDhxpNQC9duoS7d+9aXFdxE/nTp0/j008/LdPGEEsT+by8PNSuXRtqtdrirhbmjvl8ZdUiX1D+mapHWVGfy1u3biEqKgoTJ07EV199pcxPSUkpVFZE8PLLL2P48OGlEltcXFyp3CjJkhZ5/US+tLsN7dq1y+RwrfaSyOt/5j766CO0adMGycnJtgvIBCbyNqafyP/6668YO3asDaMpPQkJCbYOwa7oj1Fsj18E9kT/xiz37t0zWKafjJTFXf8Kvjf6PyharVZ57OPjY1W9o0aNwrp16wqNxGOLRF7/7BBQeJvbtWuHl19+GZGRkRbVt2TJEnz//fc4duyYwfzY2Fg0btwY1atXtzi24ibyrVq1wuTJk/H5559b3K1Pp9PhpZdesrg7hKWJfGpqKhISEiAiuHLlikV16x/zpvaBNYm8iBR6X5nIW66oFvlatWqhQYMGhT4jGo0GIoK//vpL2Y+JiYnYtGkT1q9fX+j7zFqRkZEICAhAjx49SlQPAKxYsaLIMvrfFaX9u5WYmIgqVaoYXab/vWvLbrr6jR7Tp0/H8ePHsXjxYpvFYwoTeRszNmrNw9C//GG582ZZeFQS+cjISEyaNMnq5Ez/MxEWFmawTD8ZKYtuGAXfG/1Y9BN5X19fi+u8du0a9uzZY3SZPSTyxloRAeDvv/+2qt6VK1fi2rVryvOLFy8qjy3tTqa/j4tj8uTJCAsLs+g79ODBg/jhhx8wa9Ysg/mffPKJ0RuAmUvk9RMP/ePS0u9B/ePM1HGt37XG1P68e/cu+vfvj/bt28PHxwf//e9/lWXWJPLx8fGYOnUqoqKiLHpNedi7dy+WL19eLr+PlnatKXjBeEpKCnbs2IEnnngCzZo1A2D4viUmJpYorm+//RYA8Mcff5SoHgD4+eefiyyjP7Z+af9uJSQkWJTIWzpiTlm04hc8ewmU7vUCpYWJvI0ZS+RN/bDmy8nJQceOHfHGG2+UVVglxkTetKLe34dFaGgoPvvsM0ybNq3YdcTGxho812+VLGnSZ4y51khTLfJFdecw14JjD11rTP1AW/vDvWbNGtSrV0957ufnpzy29AzdgQMHzC4/dOgQJkyYYPZP3KFDhyxan36Cld9tIC0tDe+88w6OHDlSqLy5RF4/4dGPrVu3blYf/8a2LS8vz6AbgqlY5s+fj82bNytnR1577TVlmf5ZA/3XV6pUyaCOy5cvY8SIEVi6dCkef/xxq2Iv6H//+x86dOiAc+fOGczfvHkzWrRoYdWN3f71r39h4sSJ2Lp1a4liskRxP5cZGRnYvn07gAfdbwoqaSJvKvEtD2XRIq/fIKL/51T/DJWlf6pKesG5McYSeWu6CpYXJvI2ZiyRL+pAiYyMxOHDh7FixQqldeLy5csGrX537tzBmjVrbHbTB3atMe1RaZHPZ6o1uqRKmsgbu1OmuffGVCJf1B8zcz8w5Z3IZ2ZmYvfu3Qbzbt68ibfeeqvQ/NmzZ5doXfqJcnx8vEWv2bp1q9nW+06dOuHzzz/Hu+++q8z7/vvvC5W7efNmkevatm2b8jj/+8pca29OTg5WrVpV5JmKgon4kiVLioxFn7ERMsaPH2/QR9lUIp+ammrwXP9OsabOZhVsyfz777+VkWtKej1P9+7dceTIkUL3LOjfvz/Onj2Lfv36WV2npV2+SqIkw6DqJ6cF/5iaSuQTEhLwxRdfFPm7WZxE/ujRoyaXWdM9sbS7XCUkJBjsK/3vUf0/l5autywS+YJnLwGYvfuyrTCRtzH9lpx8RSXy+heg5B/kjRo1wlNPPaWcSh04cCBGjRqFt99+u/SCtYKjtsjrdDr88ssviI2NhYiUyYWpjpTI63S6Ep/KtiSpKg5LE3kRMfpjYGxEBktb5PXPDBT1Q/P7778bPNdPVFevXo3Lly+bfX1RNBoN/vOf/xT5A6PT6QxayfMdP34cX375pdF+t5Ym4PoKXjAMALdv37b49cZaMgvS/9MxdOjQQsuLOuZ0Op1B3/hJkyYVuc569eph7NixaNu2LaKjo7F8+XKjDSXGWtSt+QzpXzyZr+AgCMb2MQCL7wyr1WqRnJyM7Ozscrmw0NSxaek1BPry74NQlkryB1s/+evcubPBe68/tOeFCxdw8OBBAA/OnLz99tuoWrUqbty4oSTYOp0OFy5cwJEjR7B+/XqDRP7YsWMYN24cdu7ciT///FPp5hMZGQmVSoX27dsjOzsb3bp1MxmrNX32S/t3Kzk5Gc7OzspzU39i9M8qmaOfyG/evLlkwf1/TOQfISKCX375BZcvX8bt27cRERGB3r17o0OHDoVaTnQ6HSZPnoz169ebrK+oRF7/4I+Pjzf4V923b19kZGQoo8Z8/fXXxdkkxc6dO/Hqq69a3fqZ3zJ04cIFkz869mjbtm14/vnnERgYiLfffhteXl6FLuIrKUfpWqPT6dCuXTt06NChRENmWtuX/cSJEwbP+/bta/SPYUJCApKTk5Genm4yvsTERFSvXh2+vr6FTskbS+Sjo6MhIrh58yaOHz9usMzUnzpTifzmzZvRtGnTQiMzFGxVa9SoEX788UcMGzYMGRkZuHDhAqpVq4YuXboo69y4caPJbiezZ8/G+PHjC11PkJaWho0bNyp1/Pzzz1Z/jq3p+pDPWCJuLpEv+N7p97M3JTo62uzyK1euYMuWLdBqtdDpdOjduzfGjx9vsvwPP/xQ5Dr1de7cGRMnToSrqytUKpVBn/rS6FKRnZ2Ne/fuYenSpfjpp58KLf/yyy/RpUsXhIeHG8yvWbOmyTr1z0B17NgRPj4+cHZ2tupPVmZmJrZu3Wpxi7V+FwlLHD58GGPGjFH66RtLHm/evInr169jxIgRhcYfLy3mtq/gH/OCCiaR+t8z+WdbcnJyEBISgrCwMDz11FMGwzUHBQXByckJiYmJaNSoEUJCQtChQwcMHz4cr776qlKuXbt2+Oqrr/DMM8+gW7duqFu3Lr766iuEhoYCeNAS36BBA7PbYs17X9YNUAkJCcjKyir05/jPP/80+Zrly5fjm2++AWB4Zmn16tWlEpOxHhN22Ugpduru3bvSq1cvcXV1lYYNG8r//vc/o+XS09Nl8ODBolarpXbt2hIREWHxOlJSUgSApKSkmCzz+++/S6tWreTEiRMiIrJixQpxdXWVVatWKWUiIyMFgNFp7ty5cvz4cXnhhRdkypQpcvDgQWXZ0aNHjb5mxYoVheK4ceOGzJgxQ6Kjo2XHjh1K2datW8vgwYMNXt+3b1+pXLmy8txS+/btky+++EJ0Op2kpqbKp59+qtQxcuRIi+rQj2Pz5s0CQNq0aSPfffed5OXlWRxLadLpdKLT6Swqu2DBgkLvR9WqVYu13qSkJGW9U6ZMUerz9/cvVDYvL08uXrxYKM67d+8WmpeUlCQLFy6Ua9eumVy3TqeTu3fvSkREhFy5csVgWWpqqvTv31++/fbbQnVnZ2fLjh07JCMjQ5KTk5WY849/rVYrUVFRIiISHx8vWq3WZAz6+1DfTz/9JP3795fk5GSJioqShg0byrhx46R79+7y7LPPGv1M+Pn5yZUrV2T06NFGl3fp0kV0Op1s3rxZ3nnnHcnMzJS7d+8WKpeQkCBnzpyRtLQ0ad26tdG6Jk6cKH5+fgJA6tevr8wPCAiQoUOHFir//fffF7n9JZlq166tPPbw8JDnnntO0tPT5dq1a/Lll18alE1PT5e7d+9KQkKCTJo0SZn/1VdfWbSuY8eOmfwuys7OlsTERMnKyhKRB9+9hw8fLlTHwYMH5ebNm7JmzRpl3kcffWR0H/373/8WV1dXg9evWbNGUlJS5MiRIxYdU+a255lnnpGrV68avP9//vmnDBs2rFDZ5ORkSU1NLbX3TX86fvy4iIj89ttvsnPnTsnIyJBt27Yp+/Ktt94qdt3Lly+X27dvS0JCgjz99NOFlv/444/yxx9/yMSJE4tV/4oVK+TmzZty+fJladu2rTL/+eefl+vXr4uIyHfffSdNmjSRgwcPSlZWlnz33XeyZs0acXd3V8pnZWVJZmam5ObmilqtVuZ36NBBXn31VenRo0ehdbu4uIhOp5OsrCyT8enLyckp8rtep9PJmTNn5OTJk8r+FxGZNm2aUufo0aOV+V9++aXMmzdPzpw5I9u3by/xsfDaa6+VyTFWnOnXX381u6/085VFixaZLHf27FlZuXKlxMbGGsxPTEw0+J3atm2byVhmzJihPH7zzTcNlh0+fLjQOvPzNwBy+vRp6dq1q/K8du3aotPp5NKlS5KTk2N2G83R/w4zdcxZw9I8RMSy/DSfSsQ+h0gZOHAgvLy8sHz5cuzatQujRo3C1atXCw37Nm3aNJw7dw4bN27EuXPn0KtXLxw9ehQNGzYsch0ajQZeXl5ISUkxegpFq9XCx8dHafF2dXUt1z6tISEhcHJywuXLl0vUZ2/MmDG4e/cukpOTERwcjCpVqsDd3R0VKlTAX3/9hZycHBw6dMig5bRSpUqF+s9NmTIFaWlpaNy4MTQaDVxcXBAQEAARQW5uLry8vPDcc8+ZjMPb2xvz5s3DrVu3ULlyZXh6eqJy5cqoXLkyTp8+DScnJ1SrVg0uLi5ITU2Fj48PcnJyEBMTg6ysLKjVavj7+yMoKEgZEcPT0xMxMTG4ffs2OnfujKtXryIlJQVNmzZFxYoVERcXh2XLlgF4MARgQECAMkTY9evX4ebmhqCgINy6dQu1atXCr7/+avSGKD4+Pujfvz/UajVOnTqFI0eOwNPTEx07dkS1atWQnZ2NQ4cO4fLly/D29kZGRobSp7VXr17KBVD6WrdujejoaKSkpBicfm3Tpg0SEhIQGxurvAft27dHcnIydDpdodPRrVq1QlRUFNRqNQIDA3Ht2rVCrQYVK1a0eoxtPz+/Eg+XVhZ12TtPT8+H/l4BtWvXLnQhsiWcnZ0REhJS6KxLcdSqVQvJycmlOnpR3bp1i2ztJ/sSFBSEhIQEg7NNfn5+yMjIKJML4sl22rZtCw8PD+zdu9eq1zVo0ACZmZkIDQ3FnTt3oFKpoNPp4OvrCycnJ+WMXL9+/aDRaODq6oqLFy+aPDvYokULZGRkoFGjRnBzc0NSUhJyc3OhVqvh7e2N7OxsJCcno3bt2nBzc8PNmzexZ88epKSkoHfv3pg1axbq168Pb29vVKxY0ej1WkXlp/rsMpHXarWoUqUKYmJilFOFnTt3xujRozFs2DCDsjVr1sQvv/yC9u3bAwCGDRuG+vXrG71QKysry+CCIY1Gg9q1ayMgIABOTk6oWLEiKlasaPImBUREREREpa1atWpQqVRQqVRKtydLEvlKZpfayJUrV+Dl5WXQ369ly5aFWkqTkpJw+/ZtNG/e3KCcqau0Fy1ahHnz5hWaX1YX4xERERERFaW4Q1vaZSKv1WoL/QPx9PQsdLGFVqtFxYoVDS5y8PT0NHk6bcaMGZg8ebLyPL9FfteuXfDw8EBeXh7y8vKQmZmJ27dv49KlS1i4cGHpbdhDwtfXF0lJSahcuTIqVaqEzMxMVKtWzeIRLtRqtfIeqVQqoyM6mOsKkt81KCMjA6mpqTYbYpOIiIjIGn369EGjRo3g5+cHHx8fNG7cGF5eXsjKyoKzszNEBKmpqYUGMDCpGP31y9yJEyfEz8/PYN6bb74p06ZNM5iXmJgoACQtLU2Zt3TpUhk4cKBF67HkYgKVSiUAZOrUqaLT6eQ///mPdOnSRZ5//nn54osvJDMz0+AimZdffrnQhRFvvPGGuLi4SNu2bQ0uIA0KCjJ5IcX9+/fl6NGjcvbsWYmKihKtVivLly+Xbdu2SXBwsNkLWN5//33l8XvvvWdy2/QvvEhPT1dec+nSJbl9+7bBBSerV68ucn/OmjVLKb9v3z557rnnBIB4enrKrFmzJDMz06L3pTTpdDr54Ycf5MyZMxaV9/T0NLpPk5OTC9Vb8MKVvLw80Wq1ygVXBw4ckNTUVBERg4smw8LCJCEhweC10dHRMmLECLl48aIy79SpU1K7dm359ttvJTc3V+7duyepqamSkpIiixYtkh9//NHkBcR79uyRxx9/XIAHFyqnpqZKTk6OZGdny3//+18lDpEHF63mX4wnIpKZmSl5eXmyZ88eg2Pivffek59++kkyMjLk+vXrcuvWLUlKSpLIyEiZO3euaDQaycvLkxs3bsju3buV137//feSmpoqx44dUz6zAOTtt9+WNWvWiIuLi8GFTsamIUOGyOeff25y+dy5c+XixYuyefNmOX78uKSkpMjZs2cLlbt69arExsaW6OLCgtPFixclNjZWoqOjlWPjxo0bpVZ/wem5556TqKgoiY6OloEDBxosO3funFy4cEF27twpQ4YMKfG60tPTJTc3V0RE/ve//8mSJUskMTFR2c7XX3+90GsiIiIkPT3d4CLdM2fOKJ+N3NxcycjIEJ1OJ40bNzb6/Zf/ecr/jOl0Orl8+bJSZt26dcp8c/F7eHgYXLCXkpIix48fl6pVqxpdb//+/cvkPfvmm28kPDxchgwZIjt27JDx48dLWFiYnD59WmJiYkpc/7fffmv0Qt2PP/5YNmzYIHPmzFHm9ezZ06q6BwwYIEuXLjX4XQEefH9s2bJFpk6dKqGhoQJANmzYIOvXr5eOHTvKY489ppT18fGR+vXrS8uWLWXnzp0Wr7tu3boSERGh/JYAKHRB74cffqh8Z/30008SGxurHGNarVYyMjJERCQrK0vi4+Pl5s2bcu3aNXnnnXckMjJSdDpdofdAo9FIQkKCHD9+XJnXuXPnMjk2bDlptVqD36L09HTls5eXl2dQNicnR5KTkyUhIUHOnj0rBw8elH/++Ud27Nghy5Ytk61btyq/I1qtVuLj4yU9PV00Go3yG2cqjuDgYDl//rzJ5R07dlRiy5eSkqIM6lFwEIL58+fL1atXZc+ePaLRaCQ3N9fgwmZL6OdE+pOpgVdKmzUXu9plIp+amiqVK1eW+Ph4ZV5YWJjy5a2vRo0aBqMbDB06VObNm2fReizZUVevXpXZs2fLvXv3TJZJS0uTgQMHyurVq0Wn04lGo5FOnTpJ06ZNjV4xvWbNGtm3b59MmDDB6IFSFP0k5MsvvzRIGjp06GDw49a/f3+L9oWIyPXr1+XAgQMG8/bt2ydTpkwx+LNkyrJly5T17tu3T7RarSxevFhJcBxBly5dCr0fY8aMKXG9+j8C9uDo0aOF/kwUlJeXJy+88ILFf4z13bx5U9negp/bPXv2yMSJEwsdUzk5OZKQkGDy85Cdna3M69mzp+Tk5BT5JZc/clKdOnWUBFTkQVJq7LNXo0YN+fnnn8XJyanQMv0/ovnzmjZtanLd8+fPL/KHtG/fvrJp0yaJj49Xkor9+/fL66+/Lp06dZKzZ8/KqVOnlPJ//fVXoR8kU98dubm58vfffyvLGjRoUGQ8Bf/gFPzxLOiLL74w+X5169atyGM+ISHBIOEzVzY3N1cps3DhQqPbnz/l19m7d28ReTDa2A8//KC8Rj9B01+v/kg/xZn8/f2Njg5TlILljx07powqlZGRYbDs3//+t0GDw/37903W8+GHH4qIyMmTJ5V5M2bMkK+//tqi7fH09DSIU6fTye3bt41uQ8FjRT8RfP3115X5+p/jd999V65duya3bt2SrKwsmTlzpnTv3l1iY2MlNjZW+czu3bvXYN9cunRJevfubXQ0k+Iy9X7dunVL7t69KyKGx2Bxp2XLlklOTo4sW7ZM+SPerl07GT9+vEGjoP7vqbFJf0QgU5OxkYCKe1yWBv3GIf0RgAIDA0XkQQPmkSNHDNbr7Oxssr7k5GRJSUkxaDgCIB988EGpxGtsnzGRt8KLL74oY8aMkfT0dPn111/Fx8fH4Ic435QpU6R3796i0WgkMjJSvLy85J9//rFoHdbsKGvpdLoih1ucPHlysT5c+h9wffp/GvKXd+vWrXgbUAzffvutst59+/aV23pL0z///CNqtVqmTp2qDE9Y8M9NcQwYMKBUvxDtnf4P3uLFi616rbnPQ/68Dh06WBVLQfqJjbF1JSQkyLVr15T5tWrVMhpHwfkF3bt3z+wPaXZ2tkXbsHPnTjl58qTRZVOnThUAMmrUKKPLY2JiZNu2bQbDtZmaCg5nWdQf+B9++MHkPtRvJTOn4DrNyS8zduxYZZ6xBpHbt2/LggULJC4uzmg9BZPj/CT/o48+MrpfCg5X+dtvv8lTTz0l169fN0hYu3XrJosWLSpRwmSsvP6yiIgIkw0D+sM6ApBx48aJiBgcyzNmzJDs7Gzp3bu3RWcgSiK/jj59+hid/8orr1hUz6VLl5TXHDt2rEQxmWLtMWhsMnaW/cSJE8rjU6dOFarv/v37Bt8Dr776qvTq1Ut0Op0kJiaaPOuUmpoqa9eulZCQEOnYsaO0aNHCYCjszz//vMgzVkXRP55Kw+nTp5X69BN5d3d3g/2hH+OgQYOKrLfgUN5JSUmlEq+xfbZr165SqbsoD0Uif/fuXenZs6e4urpKgwYNZPfu3SIi8v333xu0gqWnp8ugQYPE3d1dAgICJDw83OJ1lGUib4n8H2BrP1yHDh0qsuy6devE39+/zL70jNmyZYsSl6Mm8iL/94coLS3N4j+FRdHvpvSoyN/e8ePHF+t1xvZX/jxzLeGWMNX1xVQsXl5eRud7eHgUuS4fHx+DdegnnpYm8uZkZ2fLH3/8oXQhMMfX17fQNn/22WfKD2Z4eLjBsqLi27dvn8l9qN9dypyNGzda/P03ZcoUcXd3l3Pnzhlsf8EYLLlnhX43yHzr1q0zelx88skn8vHHHxd5nLRp00a++eYbq7/TLT0OgQeJ/C+//CLA/3WPy1ewq9Kzzz4rImJwpuudd95Ryt++fdtsote1a9ciY7dkux577DGj80NDQy2qR6PRKK/Zu3dviWIyxdL3y9j7mz+99957BmfzvvnmG4MuGsYaFSwxYsQIqV27tkX3h8lfnp8L5Xd7Kk6uERgYWKq/W/pnaguOKZ+vYJcWS/Iz/T96X3zxRanEKmI8kd+8eXOp1W+ONfmp3d7Z1c/PD9u3b0d6ejouX76Mp556CgAwePBgg9FrXF1dER4eDq1Wi9jYWAwaNMhWIVutUqXiXWscGhqKb775xuB2zwUNGzYMcXFxaNOmTXHDs1rBMf4dVf774ubmhkaNGpVKndWrVy+VehyRtduuVquVx5GRkUbLlPQugwXvrloUUxfQ59/B2JyCd9v08PCwat1FcXJyQteuXeHi4lJk2WeeeabQvLZt2+Lu3btYv369wV0of/vtNzg5OZmtr0aNGiaX1a5du8h4AOveiyVLliAxMREhISHKvKJiNMXYPS/M3Rn1zTffRM+ePZW7cxqj1WpRtWpV5bmHhwcWLVpkVVyW3GX2ueeew6lTp7Bjxw6D+W+99ZbB8/w7q+oPHqF/Z2lzw9p98cUX2Lhxo0UxF6XgwAWPPfYYAKB///4WvV7/M1PckT1KyyuvvGJyWWpqqsHx/Oqrr8LV1RX37t3D/fv3jd4p1BJr1qxBTEwMgoODiyy7b98+zJkzBwMHDgQAbNq0qVjrBB4MbFGa9PeNqTuc63+PdevWrcihFwHAy8tLeVySu5BbwpLv/PJml6PWPCr0P9R16tRBTEyMxa/Vv1WzvXhYEvmyEBAQYOsQyt2hQ4fw448/YtKkSVa9Tn9UoyeeeMJoGVM/ApbSH+nKEtbeTEuf/ueiYcOGpZ7IW2PFihWIiIgwmKfRaODn5wfA8HbygYGBRdZn7k/aoEGDsGDBApPvYT79xNcSlStXNrv8pZdeQoUKRbdR9evXD82bN1eSSsD8HxNXV1ejN3bTl5qaquxLALhx44bBPjXF2dlZucfJgAEDiiwPPBhquaBx48ahQYMGSEhIwLp167BkyRIAho1G+p8dY3/+6tWrhzfeeANvvvmmRXGYs2XLFrz77rv45ptvDObv2bMHhw8fRo8ePSyu6+WXX8aff/5p9M9oaTB2E0Rj3NzcULt2bWg0mkLfQ1qtFr6+vsoY4PmsPcaNqVChgkX1dO7cGZ07d1aeW9tooa8krzVG/3greNPCfPo3R7L0Rpj6iXxZJ9r2eLM/JvI2pJ/IW/Jlb+/0Exaxv/uM2dQrr7yC1atXo0OHDrYOpdyEhoYiNDTU6tfpt8ibUpp38yxr+p+L3bt3Y9u2bTaLRf8HL5/+D6q1f8YLtpaNHDnSYF3Xr18vMqnWTxb++9//WrX+giZOnIhPP/3UorIuLi44ffq0QeJQr169Eq2/YIu8pd+D7u7uBjcrLGjOnDlG74FSUIUKFdC9e3cADxJfY/QTEWN3lFy/fn2xPrfG9OvXD/369Ss039vbG7169bKqrg0bNiA3N7fYZ7KL4uzsbFEiDwBXr15FTk4OgoODcefOHWV+/s0sy4olf64LsuRMnSml3SKvLyEhAT4+PkhKSjJZxtJEXn8bS9rIo8/Ynzt7bJG32641jwL9RL5nz542jKR06CcBjpRolQdnZ2ccOHAAH3/8sa1DsXuWJPKO5PXXXwcAtGnTBoGBgTZtkTdGPx79BgVLklCVSoX4+HhER0cjPT0dq1evNlhu6vbj+vQTn+IkKvqKav0vqGBsbm5uJm/LbomCLfKWJoZFnSFq2LBhsWMqqKhEx9LkyRbKKokHHnxHW8rJyQlubm6YOHGiwfycnJwyTeSXLl2KFi1a4N///rfFrzH1+fvzzz+LfG1ZJvKNGzcucl8V51gszUTe2Pbrd+uzF2yRtyH9L6X27dvj8OHDCAoKsmFEJePu7q48Nvcvm8ic8k7kW7VqhVOnTln1Gk9PT4tPsfbq1QtnzpxRWnvtKZGfNm0a+vTpozzXT+QtbXky1x3FEvpJrLVdnvJFR0fj2LFjFve5Nqdu3brFfq1Op4OPjw8mTJiA9PR0g6TeHP3vTmNKM4F15ES+LFmTyOebOnUq2rVrhyNHjuA///kPFi9ejOjoaPzyyy94/PHHSz3GmjVr4vTp06VSV5cuXYosUxaJ/NmzZxEZGYmBAwfi1KlT+Oijj0x+TmydyPv4+CjXZJw7dw4nTpww+L60F0zkbajghS+O3u1C/58/E3kqLnOJvJOTE3JyckplPdu2bcM333yDlStXYuvWrUYTOHd3d6Nnl7y8vKzqK9m8eXPlsb0k8t7e3vjoo48M5umfoi7pBcXW2LFjB+Lj41G/fv1ivb5OnTqoU6dO6QZVAp999plV5e0pkX9Uz6YWpwtKxYoV8eSTT+LJJ5/EzJkzATw4e3Lu3LkS/SG0Fz179rT6Yu2iNGvWDM2aNQMAzJs3D0FBQSZ7JBQnkS/NPuytWrXCpUuXADxoibfH1niAibxNFfcKdkfQoEEDW4dADspcIh8cHKx8sZZUr169lH66o0aNMhmLscTGkpEUTLGXRN7YSDz6f8YtbU0uDU8//XS5rcse2VMi7+gNSsVVnBZ5U+w14bNWWFgY9u7dW+JrR0xxdnZWuh4aU5w/laXZIv/ll19CpVLZ5eAi+thH3oYexkT+1KlTWLt2rVWjERDpM5fIb9myBW3btsXWrVttGsvUqVMBFO/aFv0/Aba4KDz/e8dUN74//vgDK1eutLq/+aPMklFyzCkqkS/uEJvGmEqOVCoVrl69aldnNspTaSbyD5Nu3bqV+NqV4rL0GhN9pdk1s2rVqtiwYYMy/Lm9YiJvQ2V54Y6ttGzZEsOHDy/yAjciU8x9ETdp0gRHjx5F7969bRrLsGHDcObMGfzyyy9W19mgQQN069YNvXv3LnIoxbLw999/o1+/fvjtt9+MLu/atSvGjh1bzlE5tpImD+XRIj9r1iwAKDSqT36L6PLlyy0ap/xhVZLRXcj2tm7dis6dO2PFihW2DqXcPXyZpAN5GFvkiUqqVq1atg5BYSrBUqlUBv3erVGhQgXs3bu3JGGVSKtWrbBlyxabrd8RFXXmxMPDo0R9c8sjkZ8/fz7GjBlT6EZdX375JSZPnlxm3SccxcPcIv/SSy9h06ZNCA4OLtGoTPasd+/e5dbAY2/YIm9DTOSJCpswYQL69OmDr7/+2tahlKgvPD08irrAOn8Ywvwx3K1VHom8SqUyerfdChUqoH79+o/8WdSHOZGPiIhAfHy80bsZk+Nji7wNMZEnKszV1dVkt4/y9sUXXyAsLAzvvPOOrUMhG8rOzja7fNKkSXjiiScM7hJrjfK82JWMe5gT+QoVKqBGjRoO031IpVLxppJW4LeDDfHLmci+1a9fH7du3XrkWysfdUUl8hUrVkSnTp2KXX/79u3NLudvRdlzlCS3JBxlG93c3B7ZYVCLg11rbIgt8kT2j0k8FZXIl1T//v3x9ddf48SJE0aXM5Evew9zi3w+V1dXW4dgkaLOUJEhfjvYEBN5IiL7V9aJvEqlwujRo00uZyJf9h6FRN5RWuSbNm2q3FGVimZXLfK5ubno378/atWqBZVKhdu3b5stX6dOHbi5uUGtVkOtVmPcuHHlFGnpYCJPRGT/yjqRLwoT+bLHRN5+rF+/Hv3798f+/fttHYpDsLtvh86dO2Pq1KkW311u7969DnvjEn45ExHZP1sn8qV5QygyzlGS3JJwlK41tWvXxk8//WTrMByGXWWSlSpVwoQJE8qs/qysLGRlZSnPSzLub2lgizwRkf2zdSLPRp+yxxZ5clR21bWmOPr164fq1avj+eefx/Xr182WXbRoEby8vJTJ2Ji65YmJPBGR/WMi//B7FBL5xx9/3NYhUBlw6EQ+IiICMTExuHLlCgIDA9GvXz+zY4/OmDEDKSkpyhQbG1uO0RbGRJ6IyP4xkX/4PQqJfN26dXH8+HHExMTYOhQqReWayPfo0QMuLi5Gpw8//NDq+kJDQ+Hi4gJPT0988sknuHLlCqKjo02Wd3Z2hqenp8FkS/xyJiKyPwX7pNtTIq/T6WwYycPrUel28thjjyEoKMjWYVApKtdMcteuXWVWt0qlcrjxntkiT0Rkf5ycnJCTk6M8t6dEPi8vz4aRPLwehRZ5ejjZXdearKwsZGZmFnpc0I0bNxAZGYmcnBykpaVh6tSpCAoKQp06dcox2pJhIk9EZH/suUVe/w8GlZ5WrVrZOgSiYrG7RL5Ro0bKEEl16tQxGC5p3LhxyljxqampGDNmDLy9vVGnTh1ERUXh119/RYUKdrdJJjGRJyKyP5UrVzZ4bk+JfG5urg0jeXiFhobixx9/NHl3XSJ7ZXedtM1dhLFy5UrlcUhICM6ePVsOEZUd9pEnIrI/9tYirx8PE/my8+KLL9o6BCKrOU7z9UOILfJERPbH3hJ5/TPNTOSJSB8TeRtiIk9EZH/srWuN/kAOTOSJSB8TeRtiIk9EZH/srUVeHxN5ItLHRN6G2EeeiMj+FEzkZ86caaNICmMiT0T6mEnaEFvkiYjsj37XmitXrqB+/fo2jMYQE3ki0scWeRtiIk9EZH/0W+TtKYkHOI48ERliIm9DTOSJiOxPwa419oQt8kSkj4m8DbGPPBGR/Sk4ao09YSJPRPqYyNsQW+SJiOwPW+SJyFEwkbchJvJERPZn7NixAIAOHTrYOJLCmMgTkT727bAh/UReRGwYCRER5evbty8uXryIunXr2jqUQpjIE5E+JvI2pN9HPi8vz4aREBGRvsaNG9s6BKOYyBORPnatsSH9Fnl+ORMRUVH4W0FE+uwqkddqtejUqROqVKkCHx8f/Otf/8I///xjsvzVq1fRsWNHuLm54bHHHsPp06fLMdqS00/k2SJPRERFYSJPRPrsKpF3dnbG119/jXv37iEhIQEvvPAChg8fbrL8K6+8gh49eiAxMRGjRo3C888/71BfckzkiYjIGrwhFBHps6tE3snJCU2aNEGFChUgIqhQoQKio6ONlr106RIuXbqEGTNmwMXFBW+++Sby8vJw+PDhco66+NhHnoiIrOFIjVVEVPbs8mLXFi1a4OLFi9DpdPj444+Nlrlw4QIaNWpkcOOOFi1a4Pz58+jcubPR12RlZSErK0t5rtFoSjdwK7GPPBERWSMgIMDWIRCRHbGrFvl8Z86cgUajwcqVK9GkSROjZbRaLTw9PQ3meXp6QqvVmqx30aJF8PLyUqbatWuXatzWUqlUymO2yBMRkSm7du3C8OHDMW/ePFuHQkR2pFwT+R49esDFxcXo9OGHHxqUdXV1xejRozFq1CgkJSUVqkutVhdqUddoNFCr1SbXP2PGDKSkpChTbGxs6WxYKWAiT0REpnTv3h1r166Ft7e3rUMhIjtSrl1rdu3aZVV5EYFWq0V8fDx8fHwMljVt2hSXLl1CTk6OcjvtM2fOYOrUqSbrc3Z2hrOzs/WBlwMm8kRERERkDbvqWnP69Gns378f2dnZSEtLw8yZM+Ht7Y0GDRoUKtuoUSM0atQIixcvRlZWFlasWIGKFSsiNDTUBpGXHPvIExEREZE17CqRz8nJwYQJE1ClShUEBgbi1KlT2L59u9LiPm7cOIwbN04pHxERgR07dsDb2xtff/01Nm/ebDASjCNhizwRERERWcOust42bdrg5MmTJpevXLnS4Hn9+vVx6NChsg6rXDCRJyIiIiJr2FWL/KOMiTwRERERWYOJvJ1gH3kiIiIisgYTeTvBFnkiIiIisgYTeTvBRJ6IiIiIrMFE3k4wkSciIiIiazCRtxPsI09ERERE1mAibyfYIk9ERERE1mAibyeYyBMRERGRNZjI2wkm8kRERERkDSbydoJ95ImIiIjIGkzk7QRb5ImIiIjIGkzk7cQTTzxh6xCIiIiIyIFUsnUAj7pLly7h8OHDGDp0qK1DISIiIiIHwkTexho2bIiGDRvaOgwiIiIicjB227Vm8eLFUKlUOHLkiMkyderUgZubG9RqNdRqNcaNG1eOERIRERER2Y5dtsjHxcUhIiICNWrUKLLs3r172b+ciIiIiB45dtki/84772DevHlwdnYu1XqzsrKg0WgMJiIiIiIiR2R3LfJ//vkn7t+/j+effx6TJk0qsny/fv0gIggNDcVnn32GoKAgk2UXLVqEefPmFZrPhJ6IiIiI7EF+XioiRRcWO5KTkyOtWrWSs2fPiohIUFCQREZGmix/6NAhycjIkJSUFHn77belVatWotPpTJbPzMyUlJQUZTp16pQA4MSJEydOnDhx4sTJrqbY2Ngic+dybZHv0aMH9u/fb3TZrFmz4OHhgU6dOqFZs2YW1RcaGgoAcHFxwSeffAIvLy9ER0cjODjYaHlnZ2eD7jr5rfc3btyAl5eXNZtCVG40Gg1q166N2NhYeHp62jocIqN4nJIj4HFKjkBEkJqaCn9//yLLqkQsabcvH/369cP+/ftRuXJlAMC9e/fg7e2NpUuXYuTIkWZfq9Pp4OXlhdOnT5tM5AvSaDTw8vJCSkoKP9Bkt3ickiPgcUqOgMcpPWzsqo/82rVrkZmZqTxv27YtvvrqK3Tt2rVQ2Rs3biAuLg5t2rRBdnY2Zs+ejaCgINSpU6f8AiYiIiIishG7SuS9vb0NnlesWBG+vr5wc3MDAGWc+JUrVyI1NRVjxozBtWvX4ObmhtDQUPz666+oUMEuB+IhIiIiIipVdpXIFxQTE2PwfOXKlcrjkJAQnD17tkT1Ozs7Y86cOaU+zCVRaeJxSo6Axyk5Ah6n9LCxqz7yRERERERkGfZDISIiIiJyQEzkiYiIiIgcEBN5IiIiIiIHxESeiIiIiMgBMZEnIiIiInJAViXyWVlZGDlyJAICAuDl5YWuXbsaDAG5ePFi+Pn5wdfXF9OmTUP+gDi5ubno378/atWqBZVKhdu3bxvUO3nyZAQHB8PDwwNt2rTB/v37zcZx7NgxtGzZEm5ubujSpQuuX7+uLAsJCYFarVamChUqYNmyZUbruXTpEp599llUrVoVfn5+GDJkCJKSkgzKrFmzBg0aNIC7uzuaNGmCy5cvY+HChUr9zs7OcHJyUp7nj3X/+++/o379+nB3d8dzzz1nUO+9e/fQu3dvuLm5oVGjRtizZ4/Z7TW1X4vaF2RfzL3vR44cwRNPPAG1Wo2AgAD88MMPJusx955nZGRgyJAh8PDwQGBgIDZs2FCm20QPnzlz5qBp06aoUKECNm7cqMxfu3YtWrVqBQ8PDwQHBxsMB2wMj1MqS6aO06ysLLz66qvw8/ND1apVMXToUKSlpZmsh8cpOTyxglarlfnz50tsbKzk5ubKsmXLJDg4WEREtm3bJoGBgXL16lW5deuWNGnSRFavXi0iIjk5OfLZZ59JZGSkAJD4+HiDeufOnStXrlyRvLw8+eGHH8THx0c0Go3RGDIzMyUgIEBWr14tGRkZMm3aNAkLCzNa9v79+1K5cmW5cuWK0eV//fWXrFu3TpKTk0Wr1cqAAQNk1KhRyvLffvtNQkJC5Ny5c6LT6SQqKkoSEhIM6li0aJEMHz7cYN6dO3fE29tbtm/fLlqtVgYPHixDhw5Vlg8YMEBGjx4taWlpsmXLFvHx8ZHExESjMZrbr9bsC7I9U+/7rVu3xN/fX7Zu3So5OTly//59iYqKMlpHUe/51KlTpWfPnpKSkiKHDh0SLy8vuXTpUnltIj0EvvvuO9m1a5e0b99eNmzYoMxfuXKlREZGSk5Ojpw7d06qVasm+/btM1oHj1Mqa6aO0yVLlkjbtm0lISFBNBqNdO/eXWbNmmW0Dh6n9DCwKpEvKCsrS1Qqldy/f19efvllWbx4sbJs9erV0q1bt8IrNJLIF+Tv7y9///230WU7duyQxo0bK8+1Wq24urpKTExMobIrVqyQJ554wtLNkV27dkmzZs2U523btpVdu3aZfY2xRH7lypXyzDPPKM+joqLExcVFMjMzJTU1VSpXriy3bt1SloeFhcm6deuM1m9uv1qzL8i2zL3vU6dOlZkzZ1pUT1HveY0aNeTIkSPK8qFDh8q8efNKaSvoUdKlSxeDBKmgQYMGydKlS40u43FK5aXgcfrmm2/K7Nmzlef//ve/5YUXXjD6Wh6n9DAoUR/5yMhIVK9eHVWqVMGFCxfQvHlzZVnLli1x/vx5q+uMiYlBYmIi6tevb3R5wfW4u7ujXr16uHDhQqGy4eHhGDx4sMXrPnz4MEJCQgAAeXl5OHnyJM6ePYuAgADUrVsX8+fPN+jWYkrBGOvVq4dKlSrh2rVruHLlCry8vFCzZk1luf6+OnjwILy9vU3WpV/Wmn1BtmXufT927BhUKhVCQkJQs2ZNDB061KArVosWLRAREQHA/HuelJSE27dvl8rnkMicvLw8HD16VPm+BHickn0YPnw49u/fj3v37iE5ORk///wzunfvrizncUoPm0rFfWFKSgrGjh2LBQsWAAC0Wi08PT2V5Z6entBqtVbVmZOTg+HDh2Pq1Knw8vIyWqbgekytKyYmBkePHsXmzZstWvepU6fw+eefK/3z79y5g9zcXOzZswfnzp2DRqNBz549Ubt2bYwcOdJsXVqtFn5+fkZjzMzMNBp/cnIyAKBTp07KY2Pbq7+tlu4Lsj1T71VycjLi4uIQHh6OnTt3olatWhg9ejQmTpyIdevWAQDOnDlTZD1arRZarRYVK1aEm5tboWVEpWnWrFmoVasWnn76aWUej1OyB/Xq1YO3tzeqV68OlUqFf/3rXxg9erSynMcpPWyK1SKfmZmJfv36oXfv3hg1ahQAQK1WQ6PRKGU0Gg3UarXFdYoIRowYgWrVqmHu3LnKfP2LV2/cuFFoPabWFRERgaeeegrVqlUzWVe+6Oho9OnTB6tXr1ZamFxdXQEA7777Lry9vREYGIg33ngD27dvL3JbzMVoafym6tIva21dZDvm3itXV1eMHDkSDRs2hLu7O95//32Tx1lRx1ZeXh7S09MLLSMqLStXrsTmzZvx008/QaVSGS3D45RsZfz48XB3d0dKSgqSkpJQtWpVTJo0yWhZHqf0MLA6kc/NzcXLL78Mf39/LF26VJnftGlTgxFsTp8+bXDatShvvfUWbt26he+//x4VKvxfWOfPn1f+HQcGBhZaT1paGq5evYqmTZsa1BcREVGoW03BugDg9u3b6N69O95//33069dPKevj4wN/f3+D11vSrQYovC+uXbuG3NxcBAcHo0GDBkhJSTEYucfcvjK3Xy3dF2R75t73Zs2aGZQ1d5yZe899fHxQo0aNEn0OiczZtGkTFixYgJ07d6Jq1aomy/E4JVs5c+YMRo4cCQ8PD3h6emLUqFHYu3ev0bI8TumhYG2n+hEjRkiPHj0kOzvbYP7WrVslKChIrl27JvHx8RISEqKMriLy4OrwjIwMASAxMTGSkZGhLJs9e7a0atVKUlJSilx//lXma9askczMTJk+fXqhkVpOnjwp7u7uotVqzdaVnJwsLVq0MHnxysyZM6V3796i0WgkLi5OmjVrJmvWrDEoY27Umh07dkhaWpoMGTLEYNSaF198UcaMGSPp6eny66+/mh21xtx+tWRfkP0w9b7v2rVL6tSpI1evXpX09HQZPHhwoWMqX1Hv+ZQpU5RjNjIyUry8vOSff/4ppy2kh0F2drZkZGRIWFiYrF+/XjIyMiQvL0927twpfn5+cvr06SLr4HFKZc3UcTps2DAZNmyYpKWlKaPGvfTSS0br4HFKDwOrEvmYmBgBIC4uLuLu7q5M+/fvFxGRhQsXSpUqVcTb21umTp0qOp1OeW1QUJAAMJiUIABxdnY2qPP77783GcfRo0elefPm4uLiImFhYYVGaZkyZYoMGjSoyO1Zu3atADBYr7u7u7I8KytLRo8eLZ6enuLv7y+zZ8822CYR44m8yINhI4ODg8XV1VX69OljkKjfvXtXevbsKa6urtKgQQPZvXu3smz//v0GMYiY369F7QuyH+be9+XLl0uNGjXE19dXXnnlFYPjpWnTpgafB3PveXp6ugwaNEjc3d0lICBAwsPDy2fj6KExfPjwQt/Vf/zxh3Tt2lUqVapk8F05duxY5XU8Tqk8mTpO79y5I/3795cqVaqIr6+v9O3bV+Li4pTX8Tilh41KxML+IkREREREZDdKNPwkERERERHZBhN5IiIiIiIHxESeiIiIiMgBMZEnIiIiInJATOSJiIiIiBwQE3kiIiIiIgfERJ6IiIiIyAExkSciIiIickBM5ImIiIiIHBATeSIiIiIiB8REnoiIiIjIAf0/1ZANBPbtAvEAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Get a file to read\n", "filepath_raw = filepaths_raw[2]\n", "filename_raw = os.path.basename(filepath_raw)\n", "filepath_processed = os.path.join(DIRPATH_PROCESSED, filename_raw)\n", "\n", "# Loop over cases and show\n", "for filepath in (filepath_raw, filepath_processed):\n", " stream = obspy.read(filepath)\n", " stream.plot(figsize=(12, 6), show=True)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Data are never perfect... Check out station DE07, which came with a huge anomaly shortly before 4am." ] }, { "cell_type": "code", "execution_count": 59, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxUAAADtCAYAAADJNFjjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA2AUlEQVR4nO3deXwV9b3/8XcWyXY4ScCwEwIC0URAqSBVEXBBKFCDilfZAqiYa90uF1D6K4qIgF61Fi0Pri0FK1D1WhAXNotWLGJJUXbFEAigiICQhAPJIcvn94c3cznJyQInCQl5PR+PecjM93u+85mZb8b5nJnvnCAzMwEAAADAOQo+3wEAAAAAqN9IKgAAAAAEhKQCAAAAQEBIKgAAAAAEhKQCAAAAQEBIKgAAAAAEhKQCAAAAQEBIKgAAAAAEhKQCAAAAQEBIKgAAAAAEhKQCaOC8Xq/Gjh2rNm3aKDo6Wn379tW2bdt86syePVtxcXFq0qSJJk+eLDOTJBUWFur2229X69atFRQUpEOHDvl8bsKECerQoYMaN26sq666SuvWraswlvT0dHXr1k2RkZHq06eP9u3b55QlJyfL5XI5U3BwsF544QW/7ezatUuDBw/WxRdfrLi4OI0cOVLHjx/3qbNgwQJ16tRJUVFRuuyyy/TNN99o5syZTvthYWG66KKLnPm0tDRJ0sqVK9WxY0dFRUXp1ltv9Wl3/Pjxatmypdxut7p06aL333+/3G1duHChQkNDfbZp//79VdoXpeXl5WnkyJFq3Lix4uPj9Ze//KXMutq0aSO3262xY8fq9OnT5bZ1rjZs2KB+/fopNjZWLVu21IMPPuizniNHjmjQoEGKjIxUYmKi1q5d65StW7dO119/vaKiojRgwACfdqtyLEsrr79++umnPvs7MjJSwcHBOnLkiN92Fi5cqCuuuEKNGzdWhw4dNG/ePJ/y3Nxc3XPPPWrSpIliYmI0fPhwSb59NTg4WBEREc78p59+quLiYj366KOKiYlR8+bN9dvf/tZpc8+ePerRo4diY2PVpEkTpaSk6Pvvvy93W/v27avw8HCn/YEDB1Y5/tIq6nOV9TEAkAFo0Dwej02fPt0OHDhghYWF9sILL1iHDh2c8g8++MDi4+MtMzPTDh48aJdddpnNnz/fzMwKCgrspZdesg0bNpgk+/77733anjZtmmVkZFhRUZG99dZbFhsba7m5uX7jyM/PtzZt2tj8+fMtLy/PJk+ebL179/Zb9+jRo9aoUSPLyMjwW/7Pf/7TXnvtNcvOzjaPx2PDhg2zcePGOeXvvvuuJScn2/bt2624uNh2795tP/74o08bs2bNstTUVJ9lP/zwg8XExNiKFSvM4/HYiBEjbNSoUU75V199Zfn5+WZmtnHjRouOjrZjx475jXHBggV2yy23BLwvzMwmTZpkAwcOtJycHFu/fr1FR0fbrl27zMxs69atFhsba+np6ZadnW19+/a1qVOnltvWuVq5cqUtW7bMPB6PHT161Hr37m3Tp093yocNG2b33nuvnTx50pYtW2axsbHOvklPT7fFixfb008/XWafVHYsS6uov5b2u9/9zq677rpy25o3b55t2LDBCgoKbPv27dasWTP75JNPnPLbbrvNHn74YcvOzrbTp0/bF198UaaNxMRE+/jjj32W/f73v7crr7zSfvjhB/v666+tRYsWtnbtWjMzy83NtT179lhxcbF5vV577LHH7M477yw3xj59+thf/vKXc4r/TJX1uYr6GACYmZFUAPDh9XotKCjIjh49amZmd911l82ePdspnz9/vvXr16/M5/wlFaW1atXK/vWvf/ktW7VqlV166aXOvMfjsYiICMvKyipTd+7cudarV68qbY+Z2Zo1a+zyyy935nv06GFr1qyp8DP+kop58+bZgAEDnPndu3dbeHi4k0icKT093cLCwmzHjh1+268oqTibfWFm1qJFC/v888+d+VGjRtlTTz1lZmaPP/64paWlOWVr16619u3b+22nOr366qs2ePBgMzM7ceKENWrUyA4ePOiU9+7d21577TWfz1S0T0qUPpalVbW/mpn17NnT5s2bV+m2lBg+fLg9//zzZma2fft2a9eunRUWFlb4GX9JRa9eveyNN95w5qdOnWpjx44t81mv12tTpkyxHj16lNt+RUlFRfGXVlmfq6iPAYCZGY8/AfCxYcMGNW/eXE2bNpUk7dy5U126dHHKu3Xrph07dpx1u1lZWTp27Jg6duzot7z0eqKionTJJZdo586dZeouXrxYI0aMqPK6P/vsMyUnJ0uSioqK9OWXX2rbtm1q06aN2rdvr+nTpzuPyFSkdIyXXHKJQkNDtWfPHmfZAw88oIiICPXo0UMDBgxQUlKSJOkf//iHYmJifNpbv369mjZtqqSkJJ9HUyrbF7Nnz9bgwYMlScePH9ehQ4fKPUb+jt/evXuVl5dX6fYG4sx9npGRoejoaLVs2dJvjOfarj9V7a+7d+/W5s2bNWzYsCqtt6ioSBs3bnTW/a9//UudO3fWyJEj1bRpU/Xs2VOffvppldqqSowxMTGKiIjQ888/r//8z/90lj/wwAN64IEHfOo+9NBDiouL080336ytW7dWKX5J6tq1q5YsWeI3pjP7XGV9DAAkKfR8BwCg7sjJydH999+vZ555xlnm8XjkdrudebfbLY/Hc1btFhQUKDU1VZMmTVJ0dLTfOqXXU966srKytHHjRi1durRK6968ebPmzJnjjOf44YcfVFhYqLVr12r79u3Kzc3VwIED1bZtW40dO7bCtjwej+Li4iqMce7cuXr55Zf18ccf+4xNue6665Sdne3M9+nTR9u2bVN8fLzS09M1dOhQNW/eXEOHDq10Xzz++OM+MYWEhCgyMtJvXX/Hr2R5REREhdt7rlauXKmVK1dqy5YtfmMoiePM/VEVpY+lP1Xtr4sXL9aAAQPUpEmTKq37N7/5jVq3bq1bbrlFkvTdd9/pww8/1B//+Ef9+c9/1nvvvaeUlBRlZGRU2mZVYszOzlZubq7+8Ic/KCEhwVk+d+5cn3rPPfeckpKSFBISopdfflm/+MUv9PXXX8vlclUYvySfBKSiPldZHwMAiYHaAP5Xfn6+UlJSNGjQII0bN85Z7nK5lJub68zn5uaWuWCpiJlpzJgxatasmaZNm+YsP3Mw6/79+8usp7x1LVmyRDfddJOaNWtWblsl9u7dqyFDhmj+/PnON7QlF9KPPfaYYmJiFB8fr1/96ldasWJFpdtS1RhDQkJ00003ae3atVq9erXfttq3b6+EhAQFBwfr6quv1sMPP6xly5ad1XpK6hYVFenUqVN+6/o7fiXLK3LmoOYzB/9WJj09XWPHjtU777yj5s2bn/X2lMffsfQXY1X765IlS3zudlW0vfPmzdPSpUv19ttvKygoSNJP/ah9+/a65557dNFFF+m2227TJZdcog0bNlS6LVWN0e12KzU1VUOHDi33TlrPnj3lcrkUERGhyZMny+VyaePGjZXGX1lMZ8ZVWR8DAImkAoB+eovTXXfdpVatWun555/3KUtKSvL5xn3Lli0VPn5S2kMPPaSDBw9q0aJFCg7+v1POjh07nG9B4+Pjy6zn5MmTyszMdB4fKlH6YtBfW5J06NAh3XzzzZo6dapSUlKcurGxsWrVqpXP56vy6JNUdl/s2bNHhYWF6tChg9/6xcXFyszMrFLbZ+6bqu4L6aftadGiRbnHyN/xa9++faV3KXr37u3s05UrV1ZpG7766iv98pe/1J/+9Cf16tXLWd6pUyfl5OT4vB3sbPpRecfSX4xV6a/p6en6/vvvNWTIkEq3980339Qzzzyj1atX6+KLL3aWX3755WXiPNd+VNG+KC4u1vfff+9zQV+RM/uRVH78lcV0Zp+rrI8BgCTe/gTAbMyYMda/f387ffp0mbL333/f2rVrZ3v27LHvv//ekpOTfd6mk5+fb3l5eSbJsrKyLC8vzyl74okn7IorrrCcnJxKYyh5+8yCBQssPz/fHn/88TJvPPryyy8tKirKPB5PhW1lZ2db165dyx1I+utf/9oGDRpkubm59t1339nll19uCxYs8KlT0dufVq1aZSdPnrSRI0c6b386ceKELVq0yE6cOGEFBQX29ttvW3h4uG3dutVvDCtXrrTDhw+bmdmmTZusdevW9uabb1Z5X5xp4sSJzvZs2LDBoqOj7euvvzazn97+1KRJE9u0aZNlZ2fbDTfcUCNvf9q/f7+1bdvWFi5c6Lf8jjvusPHjx9upU6ds+fLlPm9/Kioqsry8PHv11Vft5ptvtry8PKcvVnYsS6usv5qZPfLII2WOrT+rV6+2uLg427JlS5my06dPW0JCgi1cuNAKCwtt+fLl1rRp0zJvEfM3UPuVV16xn/3sZ3b48GHbtWuXtWzZ0v72t7+Zmdknn3ximzZtssLCQjt27JiNGjXKevbs6Te+48eP25o1ayw/P9+8Xq+9+OKL1rx5c+fvraL4S6usz1XUxwDAjLc/AQ1eVlaWSbLw8HCLiopypnXr1jl1Zs6caU2bNrWYmBibNGmSFRcXO2Xt2rUzST5TCUkWFhbm0+6iRYvKjWXjxo3WpUsXCw8Pt969e5d529HEiRNt+PDhlW7TwoULTZLPeqOiopxyr9dr9957r7ndbmvVqpU98cQTPttk5j+pMPvplaUdOnSwiIgIGzJkiHNh7PF4rF+/fhYdHW1ut9u6d+9uS5cudT63bt06nxgmTJhgcXFxFhUVZZ07d7Y5c+ZUeV8888wzPm+hOnXqlA0fPtyioqKsTZs2tnjxYp+2FixYYK1atTKXy2Wpqal+31YVqGnTpllQUJDP/k5KSnLKDx8+bAMHDrSIiAjr1KmTffjhh07Zxx9/XKYPlez7yo6lPxX118LCQmvRokWlb/8yM+vbt6+Fhob6rPf+++93yrds2WJXXXWVRUVF2RVXXOH3da3+koqioiJ75JFHLDo62uLi4uyFF15wyj744AO79NJLLSoqyuLi4uyOO+7wOfb333+/E8Phw4ftZz/7mUVFRVlsbKz169fPNm3aVOX4k5KSfP4eK+pzlfUxAAgyq+L9WgAAAADwgzEVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVZ2HdunUaMmSIWrVqpaCgIL3zzjtn3cbq1avVq1cvNW7cWHFxcbr99tu1d+/e6g8WAAAAqCUkFWfh5MmT6tatm1555ZVz+vyePXt066236oYbbtDmzZu1evVqHT16VLfddls1RwoAAADUniAzs/MdRH0UFBSkZcuWKSUlxVl2+vRp/eY3v9HixYuVnZ2tyy+/XM8++6z69u0rSXr77bd19913y+v1Kjj4p3zuvffe06233iqv16uLLrroPGwJAAAAEBjuVFSjsWPHav369XrjjTe0detWDRs2TAMGDFBGRoYk6aqrrlJISIgWLFigoqIi5eTk6PXXX1f//v1JKAAAAFBvcafiHJW+U5GZmalOnTrp22+/VatWrZx6N910k3r27KmZM2dK+mlcxrBhw/Tjjz+qqKhIP//5z7VixQrFxMSch60AAAAAAsedimryxRdfyMzUuXNnuVwuZ/rkk0+UmZkpSTp06JDuvfdepaamKj09XZ988okaNWqkO+64Q+R2AAAAqK9Cz3cAF4ri4mKFhIRo06ZNCgkJ8SlzuVySpN///vdyu9167rnnnLJFixapbdu2+uc//6levXrVaswAAABAdSCpqCZXXnmlioqKdPjwYfXu3dtvnVOnTpVJOErmi4uLazxGAAAAoCbw+NNZ8Hg82rx5szZv3ixJ2rt3rzZv3qz9+/erc+fOGjFihEaPHq2lS5dq7969Sk9P17PPPqsVK1ZIkgYNGqT09HRNnz5dGRkZ+uKLLzR27Fi1a9dOV1555XncMgAAAODcMVD7LPz9739Xv379yixPTU3VwoULVVBQoBkzZujPf/6zvvvuOzVt2lQ///nP9dRTT6lLly6SpDfeeEPPPfecvvnmG0VGRurnP/+5nn32WV166aW1vTkAAABAtSCpAAAAABAQHn8CAAAAEBAGaleiuLhYBw8eVOPGjRUUFHS+wwEAAEADZ2Y6ceKEWrVqpeDgunGPgKSiEgcPHlTbtm3PdxgAAACAjwMHDqhNmzbnOwxJJBWVaty4saSfDprb7T7P0QAAAKChy83NVdu2bZ3r1LqApKISJY88ud1ukgoAAADUGXXp0XySihrw3Xffadu2bfJ6vec7lHovKChILVq0UPfu3RUaSncFAACoi+r1VdqRI0c0ZswYffzxx2rbtq3mzp2rG2+8sUy9vLw83XfffVq+fLliY2P17LPP6u677672ePLz8zVlyhR9+umn1d52Q+d2u/Xss8+qR48e5zsUAAAAlFKvk4pf/epXatWqlY4ePao1a9Zo2LBhyszMVGxsrE+9J598UseOHdN3332n7du36xe/+IV+9rOfqXPnztUaz4svvuj8YnafPn0UGRlZp25L1UdFRUXKzMzU7373Oz366KN677331KRJk/MdFgAAAM5Qb3/8zuPxqGnTpsrKylLLli0lSddff73uvfdejR492qduy5Yt9c477+jqq6+WJI0ePVodO3bUE088UaZdr9fr89hSyUCYBx98UGFhYeXGU1xcrOXLl+uee+7RqFGjqmMTcYacnBzdcccdSkxMVKdOnaqlzXra9fG/So5fcHCwgoKCZGYqLi6WmVX4er3U1FR17dq12uP585//rK1bt5aJLyQkREVFRQoODnbiKxEUFORMpZXUKyoqkpk521jy35LtPlNd+xKjof6NVXQ8a1td6xNAXXKuf5d14e/K6/XqlVdeUU5OTp0Z81tv71RkZGQoOjraSSgkqVu3btqxY4dPvePHj+vQoUPq0qWLT72NGzf6bXfWrFl66qmnyixftWqVQkJCyo2nsLBQHo9H8fHxysnJOdvNQRU0a9ZMn3/+ub755hu/5edycqgLJwacm5JjV3KhHhQU5CQTxcXFfo/tnj17VFhYqN/97nfVHs+ECRMUFBSkiy++2FlmZioqKlJoaKiKi4t9EgEzc6bytk1SmeShqKhIISEhPp+rqxfvJcdF8t2muhovqteZx7+62ivd7vk6h1f3ttVVtbmdJN5np7Cw8HyHUEa9TSo8Hk+ZzMztdis7O7tMvZCQEEVGRvrU83g8ftudMmWKJkyY4MyX3KnYtGlThZngwYMH9ctf/lKdO3fW5Zdffg5bhMq0adNGgwcP1uTJk893KKinunXrVmP/4zIzTZ48mf4JAKhxubm5io6OPt9h+KgbP8F3Dlwul3Jzc32W5ebmyuVylalXVFSkU6dOVVivRFhYmPP62EBfI3vHHXfoP/7jP3yWPf/88+rTp4/atWunzz//3KcsLS1N06ZN89tW3759FR4ersaNGysmJkbXXXedFi5c6FMnISFBkZGRcrlczjR//nxJP31z++ijjyomJkbNmzfXb3/7W+dzixcv9vlMyXqqYuHChQoNDXU+e9lll+mdd97xKR8wYIDPZw4dOuTzzUDJtp0Zw5NPPlml9QMAAOD8q7dJRadOnZSTk6NDhw45y7Zs2aLk5GSferGxsWrRooW2bdtWYb2a8PLLL+v111/XF198IUnav3+/Zs2apf/+7/8+p9ttCxcu1IkTJ7Rv3z49+uij+vWvf63HH3/cp85HH30kj8fjTPfcc48kad68eVq3bp2++eYbrVu3Ts8995w++ugjSdKIESN8PvPAAw8oJSWlynHddNNNzmdfeukljRgx4qwfAVu4cKFPDP4eQQMAAEDdVG+TCpfLpV/+8pd68sknlZeXp3fffVfbt2/XkCFDytQdOXKknn76aZ04cUKff/653n33Xf3bv/1bjcfYsmVLzZgxQ+PHj1dRUZEeeughPfzww7r00ksDajc6Olp33HGH5s2bpxdffFFHjhyp9DOvv/66HnvsMTVr1kyJiYm67777tGjRojL1iouL9cYbb2jkyJHnFNstt9yiqKgoZWVlndPngZpWk48/AQDQUNXbpEKS5s6dqwMHDqhp06aaOHGi3nrrLcXGxmrx4sU+dyKmT5/uDOoeNmyY5s6dq8TExFqJ8f7771d4eLhuu+02ZWRkaMqUKdXW9i9+8QsVFxcrPT290ro7d+4sM1i99KB26ac7HUVFRbrpppvOOh4z06pVq1RcXKyOHTue9ecBAABQP9XbgdqSFBcXpxUrVpRZPmLECI0YMcKZj4iI0OLFi2szNEdQUJBeffVVJScn66OPPlKjRo2csptvvtnnjVJ5eXlnlXSEhobq4osv1vHjx8ttc9WqVerVq1eZge3lDVZfvHix7rrrrgrfdFXa2rVrFRMTo4KCAuXn5+ull15SVFRUmfIS/r7Rveeee5SWlubMv/rqq7rzzjurHANQFTX9lo/6+hYRAAACVa+TivoiKSlJknTZZZf5LP/www/Vq1cvZ/7Mi+qqKCws1NGjR31+7K90myVKD2z3N1g9Pz9fS5cu1YcffnhWcdx4441atWqVJOmrr77SjTfeqISEBOdRtDPLpZ8Gap/5KmBJmj9/vu66666zWi8AAADqhnr9+FNDt3LlSgUHB6tnz56V1k1KSqp0sPp7772nZs2aVam98lx22WW69tprtXr16nNuAwAAAPULSUU9lJubq6VLlyotLU0TJkzw+bGt8owcOVL/9V//pSNHjuibb77RH//4R59HxKSfHn0qvexs7d69W+vXry9zVwaoKxioDQBA9SOpqKNKDzaXpDFjxsjlcqlt27Z64YUXNGPGDM2ePdunzg033ODzew9PP/20JOnf//3fdd1116lTp0667rrrNHHiRN14443O544fP66VK1f6TSoGDhyomTNnOvMul0uffvqpM/+3v/1NLpdLUVFR6tOnj+6+++6zfpSrZNtKpnvvvfesPg8AAIDzJ8j4eq1CJb9YmJOTU6Vf1J47d25Ajw+hfMOHD9cVV1zBLxbjnF155ZW69tpr9corr1R727Gxsfp//+//aeLEidXeNgAAZ6rq9Wlt4k4FAAAAgICQVKDe4KYaAABA3URSUU1Kfpfh2LFj5zmSC5OZ6dixYz6/fwGcCwZqAwBQ/Ugqqkl0dLQ6d+6s999/X8XFxec7nAvOhg0b9OOPPzJeBQHhx+8AAKgZ/PhdNbrvvvv02GOPafTo0erdu7eioqK4yAhQUVGRdu/erb/97W/q2bOnunfvfr5DAgAAQCkkFdWoX79+evnll/U///M/+utf/yqv13u+Q6r3goOD1bx5c40ZM0ajR49WSEjI+Q4JAAAApZBUVLOrr75aV1999fkOA0A5GFMBAED1Y0wFAAAAgICQVABoMBioDQBAzSCpAAAAABCQeptUpKenq1u3boqMjFSfPn20b9++cusmJCQoMjJSLpdLLpdLaWlptRgpAAAAcGGrl0mF1+vVbbfdpkceeUTHjh1Tr169NGrUqAo/89FHH8nj8cjj8WjevHm1FCmAuoaB2gAAVL96+fanv//973K5XBo3bpwk6YknnlBcXJz27dundu3anefoAAAAgIalXt6p2Llzp7p06eLMR0VF6ZJLLtHOnTvL/UxKSoqaN2+uoUOHVviolNfrVW5urs8E4MLAQG0AAGpGvUwqPB6P3G63zzK32y2Px+O3/pIlS5SVlaWMjAzFx8crJSWl3EcVZs2apejoaGdq27ZttccPAAAAXEjqZFLRv39/hYeH+51mzJghl8tV5g5Cbm6uXC6X3/auueYahYeHy+1268UXX1RGRob27t3rt+6UKVOUk5PjTAcOHKj27QMAAAAuJHVyTMWaNWsqLF+9erVeffVVZ/7kyZPKzMxUUlJSpW0HBQVV+IhCWFiYwsLCqh4sgHqFgdoAAFS/OnmnojJ9+/aVx+PRwoUL5fV6NWPGDF111VV+B2nv379fGzZsUEFBgU6ePKlJkyapXbt2SkhIqP3AAZxXjKkAAKBm1MukIiwsTEuXLtWLL76omJgYrV+/Xq+//rpTnpaW5vwWxYkTJzR+/HjFxMQoISFBu3fv1vLlyxUcXC83HQAAAKhz6uTjT1XRo0cPbd261W/Zmb9DkZycrG3bttVWWAAAAECDw9f1ABoUxj4AAFD9SCoAoBqQrAAAGjKSCgANBgO1AQCoGSQVAAAAAAJCUgEAAAAgICQVABoUfvwOAIDqR1IBAAAAICAkFQAaDAZqAwBQM0gqAAAAAASEpAIAAABAQEgqADQoDNQGAKD6kVQAaDAYUwEAQM0gqQAAAAAQEJIKAAAAAAEhqQDQoDD2AQCA6kdSAQDVgGQFANCQ1cukorCwULfffrtat26toKAgHTp0qML6mZmZuvbaaxUZGanu3btry5YttRQpgLqEgdoAANSMeplUSNL111+vv/71r1Wqe/fdd6t///46duyYxo0bp6FDh6qwsLCGIwQAAAAahnqZVISGhuqRRx5Rr169Kq27a9cu7dq1S1OmTFF4eLgefPBBFRUV6bPPPvNb3+v1Kjc312cCAAAAUL56mVScjZ07dyoxMVGNGjVylnXt2lU7duzwW3/WrFmKjo52prZt29ZWqABqAT9+BwBA9bvgkwqPxyO32+2zzO12y+Px+K0/ZcoU5eTkONOBAwdqI0wAAACg3qqTSUX//v0VHh7ud5oxY8ZZteVyuco8wpSbmyuXy+W3flhYmNxut88E4MLAQG0AAGpG6PkOwJ81a9ZUW1tJSUnatWuXCgoKdNFFF0mStm7dqkmTJlXbOgAAAICGrE7eqagKr9er/Pz8Mv8uLTExUYmJiZo9e7a8Xq/mzp2rkJAQXXPNNbUZLgAAAHDBqrdJRWJioiIiIiRJCQkJzr8lKS0tTWlpac78kiVLtGrVKsXExOgPf/iDli5dqtDQOnmTBkANY6A2AADVr95eWWdlZZVbNm/ePJ/5jh07av369TUcEYC6jjEVAADUjHp7pwIAAABA3UBSAQAAACAgJBUAGhTGPgAAUP1IKgCgGpCsAAAaMpIKAA0GA7UBAKgZJBUAAAAAAkJSAQAAACAgJBUAGhR+/A4AgOpHUgEAAAAgICQVABoMBmoDAFAzSCoAAAAABISkAgAAAEBASCoANCgM1AYAoPqRVABoMBhTAQBAzSCpAAAAABCQeplUFBYW6vbbb1fr1q0VFBSkQ4cOVVg/ISFBkZGRcrlccrlcSktLq6VIAQAAgAtfvUwqJOn666/XX//61yrX/+ijj+TxeOTxeDRv3rwajAxAXcbYBwAAql/o+Q7gXISGhuqRRx6pkba9Xq+8Xq8zn5ubWyPrAXBhIVkBADRk9fZOxdlKSUlR8+bNNXToUO3bt6/cerNmzVJ0dLQztW3bthajBFCTGKgNAEDNaBBJxZIlS5SVlaWMjAzFx8crJSWl3G8Vp0yZopycHGc6cOBALUcLAAAA1C91Mqno37+/wsPD/U4zZsw46/auueYahYeHy+1268UXX1RGRob27t3rt25YWJjcbrfPBAAAAKB8dXJMxZo1a2qs7aCgIB5RABowxj4AAFD96uSdiqrwer3Kz88v8+/S9u/frw0bNqigoEAnT57UpEmT1K5dOyUkJNRitAAudCQrAICGrN4mFYmJiYqIiJD00+9QlPxbktLS0pzfojhx4oTGjx+vmJgYJSQkaPfu3Vq+fLmCg+vtpgM4RwzUBgCgZtTJx5+qIisrq9yyM3+HIjk5Wdu2bauFiAAAAICGia/rAQAAAASEpAJAg1JTYx8YUwEAaMhIKgAAAAAEhKQCQIPBQG0AAGoGSQUAAACAgJBUAGhQGPsAAED1I6kAAAAAEBCSCgANBmMqAACoGSQVAAAAAAJCUgEAAAAgICQVABoUBmoDAFD9SCoAIEAkKgCAho6kAkCDwUBtAABqBkkFAAAAgICQVAAAAAAISL1MKnbt2qXBgwfr4osvVlxcnEaOHKnjx4+XWz8zM1PXXnutIiMj1b17d23ZsqUWowVQl9TE+AfGVAAAGrp6mVTk5OTozjvvVGZmprKysnT69GlNnDix3Pp33323+vfvr2PHjmncuHEaOnSoCgsLazFiAAAA4MJVL5OKnj17avTo0YqOjlZUVJTuu+8+bdy40W/dXbt2adeuXZoyZYrCw8P14IMPqqioSJ999lktRw3gfGOgNgAANaNeJhWlffbZZ0pOTvZbtnPnTiUmJqpRo0bOsq5du2rHjh1+63u9XuXm5vpMAAAAAMoXer4DCNTmzZs1Z84crVu3zm+5x+OR2+32WeZ2u+XxePzWnzVrlp566qlqjxMAAAC4UNXJOxX9+/dXeHi432nGjBlOvb1792rIkCGaP39+uXcqXC5XmbsNubm5crlcfutPmTJFOTk5znTgwIHq2zAA5x0DtQEAqH518k7FmjVrKq1z6NAh3XzzzZo6dapSUlLKrZeUlKRdu3apoKBAF110kSRp69atmjRpkt/6YWFhCgsLO6e4AdRtjKkAAKBm1Mk7FZXJycnRLbfcotGjR2v8+PEV1k1MTFRiYqJmz54tr9eruXPnKiQkRNdcc00tRQsAAABc2OplUvHOO+9o69ateu655+RyuZypRFpamtLS0pz5JUuWaNWqVYqJidEf/vAHLV26VKGhdfImDQAAAFDv1Msr69TUVKWmppZbPm/ePJ/5jh07av369TUdFoB6gPEPAABUv3p5pwIA6hISFQBAQ0dSAaDBYKA2AAA1g6QCAAAAQEBIKgAAAAAEhKQCQIPCj98BAFD9SCoAAAAABISkAkCDwUBtAABqBkkFAAAAgICQVAAAAAAICEkFgAaFgdoAAFQ/kgoADQZjKgAAqBkkFQAAAAACQlIBAAAAICAkFQAaFMY/AABQ/UgqACBAJCoAgIaOpAJAg8FAbQAAaka9TCp27dqlwYMH6+KLL1ZcXJxGjhyp48ePl1s/ISFBkZGRcrlccrlcSktLq8VoAQAAgAtbvUwqcnJydOeddyozM1NZWVk6ffq0Jk6cWOFnPvroI3k8Hnk8Hs2bN6+WIgUAAAAufKHnO4Bz0bNnT/Xs2dOZv++++zRhwoRqadvr9crr9Trzubm51dIugLqB8Q8AAFS/enmnorTPPvtMycnJFdZJSUlR8+bNNXToUO3bt6/cerNmzVJ0dLQztW3btrrDBXCBIVEBADR09T6p2Lx5s+bMmaOpU6eWW2fJkiXKyspSRkaG4uPjlZKSUu5FwJQpU5STk+NMBw4cqKnQAdQyBmoDAFAz6mRS0b9/f4WHh/udZsyY4dTbu3evhgwZovnz51d4p+Kaa65ReHi43G63XnzxRWVkZGjv3r1+64aFhcntdvtMAAAAAMpXJ8dUrFmzptI6hw4d0s0336ypU6cqJSWlym0HBQXxbSIAAABQjerknYrK5OTk6JZbbtHo0aM1fvz4Cuvu379fGzZsUEFBgU6ePKlJkyapXbt2SkhIqJ1gAdQpNTH+gTEVAICGrl4mFe+88462bt2q5557zvntCZfL5ZSnpaU5v0Vx4sQJjR8/XjExMUpISNDu3bu1fPlyBQfXy00HEADGVAAAUDPq5ONPlUlNTVVqamq55Wf+DkVycrK2bdtWG2EBAAAADRJf1wMAAAAISL28UwEA5+rw4cP6+OOPq7XNvLy8am0PAID6hqQCQIPRpEkTrVixQjfccEONtB8TE1Mj7QIAUNcFGa8tqVBubq6io6OVk5PDb1YA9ZzX69W3335bI22HhoaqXbt2NdI2AABnqovXp9ypANBghIWF6ZJLLjnfYQAAcMFhoDYAAACAgJBUAAAAAAgISQUAAACAgDCmohIl49hzc3PPcyQAAADA/12X1qX3LZFUVOLHH3+UJLVt2/Y8RwIAAAD8nx9//FHR0dHnOwxJJBWVatKkiSRp//79deagAaXl5uaqbdu2OnDgQJ15tRxQGv0U9QH9FPVBTk6O4uPjnevUuoCkohLBwT8NO4mOjubkgjrP7XbTT1Hn0U9RH9BPUR+UXKfWBXUnEgAAAAD1EkkFAAAAgICQVFQiLCxMTz75pMLCws53KEC56KeoD+inqA/op6gP6mI/DbK69C4qAAAAAPUOdyoAAAAABISkAgAAAEBASCoAAAAABISkAgAAAEBASCoAAAAABOSskgqv16uxY8eqTZs2io6OVt++fbVt2zanfPbs2YqLi1OTJk00efJklbxYqrCwULfffrtat26toKAgHTp0yKfdCRMmqEOHDmrcuLGuuuoqrVu3rsI40tPT1a1bN0VGRqpPnz7at2+fU5acnCyXy+VMwcHBeuGFF/y2s2vXLg0ePFgXX3yx4uLiNHLkSB0/ftynzoIFC9SpUydFRUXpsssu0zfffKOZM2c67YeFhemiiy5y5tPS0iRJK1euVMeOHRUVFaVbb73Vp93x48erZcuWcrvd6tKli95///1yt3XhwoUKDQ312ab9+/dXaV+UlpeXp5EjR6px48aKj4/XX/7ylzLratOmjdxut8aOHavTp0+X2xbO3pEjRzRo0CBFRkYqMTFRa9eulSQtX75cl156qaKjo9WyZUtNmDBBRUVF5bYTFBSkqKgopz/MnDnTKavsGJ/p888/V9euXRUTE6NmzZopNTVVHo/HKc/MzNS1116ryMhIde/eXVu2bKmGvYC6rjb66dy5c3XFFVcoNDRUs2fPrjSmis5N9NOGqa71U86n8Kc2+unZXkPX6PnUzoLH47Hp06fbgQMHrLCw0F544QXr0KGDmZl98MEHFh8fb5mZmXbw4EG77LLLbP78+WZmVlBQYC+99JJt2LDBJNn333/v0+60adMsIyPDioqK7K233rLY2FjLzc31G0N+fr61adPG5s+fb3l5eTZ58mTr3bu337pHjx61Ro0aWUZGht/yf/7zn/baa69Zdna2eTweGzZsmI0bN84pf/fddy05Odm2b99uxcXFtnv3bvvxxx992pg1a5alpqb6LPvhhx8sJibGVqxYYR6Px0aMGGGjRo1yyr/66ivLz883M7ONGzdadHS0HTt2zG+MCxYssFtuuSXgfWFmNmnSJBs4cKDl5OTY+vXrLTo62nbt2mVmZlu3brXY2FhLT0+37Oxs69u3r02dOrXctnD2hg0bZvfee6+dPHnSli1bZrGxsXbs2DH79ttv7fDhw2ZmduzYMbvxxhtt7ty55bbj72+oREXHuLQjR47Yt99+a2ZmJ0+etJEjR9rkyZOd8h49eti0adMsLy/PXn75ZWvfvr0VFBSc6+ajnqiNfrps2TJ77733bOjQoTZr1qwK46ns3EQ/bZjqWj/lfAp/aqOfns01dE2fT88qqSjN6/VaUFCQHT161O666y6bPXu2UzZ//nzr169f2RVWsGNKtGrVyv71r3/5LVu1apVdeumlzrzH47GIiAjLysoqU3fu3LnWq1evqm6OrVmzxi6//HJnvkePHrZmzZoKP+MvqZg3b54NGDDAmd+9e7eFh4c7icSZ0tPTLSwszHbs2OG3/YqSirPZF2ZmLVq0sM8//9yZHzVqlD311FNmZvb4449bWlqaU7Z27Vpr376933Zw9k6cOGGNGjWygwcPOst69+5tr732mk+9kpPLpEmTym2ror+hio5xRTwej40aNcqGDRtmZmZff/21ud1u83q9Tp34+Hj75JNPKm0L9Vdt9dMSqamplV6sVXRuop82THWxn56J8ynMar+flqjoGrqmz6cBjanYsGGDmjdvrqZNm2rnzp3q0qWLU9atWzft2LHjrNvMysrSsWPH1LFjR7/lpdcTFRWlSy65RDt37ixTd/HixRoxYkSV1/3ZZ58pOTlZklRUVKQvv/xS27ZtU5s2bdS+fXtNnz7deaSrIqVjvOSSSxQaGqo9e/Y4yx544AFFRESoR48eGjBggJKSkiRJ//jHPxQTE+PT3vr169W0aVMlJSVp3rx55a6n9L6YPXu2Bg8eLEk6fvy4Dh06VO4x8nf89u7dq7y8vEq3F5XLyMhwbnOWOHP//+Mf/1B0dLSaNGmibdu2ady4cU69rl27asmSJT7tde/eXa1bt9aYMWP0448/Sqr8GO/fv18xMTE+j8+VLHO5XFq2bJkefPBBST/1h8TERDVq1MgnjnP5m0b9URv9tDKl+2lF5yb6acNUF/vpmcs4n0I6P/209DV0bZ9PzzmpyMnJ0f33369nnnlGkuTxeOR2u51yt9vt8zxhVRQUFCg1NVWTJk1SdHS03zql11PeurKysrRx40bdeeedVVr35s2bNWfOHE2dOlWS9MMPP6iwsFBr167V9u3b9cknn+jNN9/UwoULK22rKjHOnTtXHo9HH374ofr06eMsv+6665Sdne3M9+nTR9u2bdORI0e0YMECTZ8+XcuWLavSeh5//HFnvIbH41FISIgiIyP91vV3/EqWI3CVHavrrrtOOTk52rNnj9LS0hQbG+vU27p1q4YPH+7Mr1u3Tvv27dPmzZt16tQp50RU2TGOj49Xdna24uPjnfKSZT/88IMmT57snPyq+neGC0tt9NPKlO6nFZ2b6KcNU13sp2cu43wKqfb7qb9r6No+n55TUpGfn6+UlBQNGjTI2TCXy6Xc3FynTm5urlwuV5XbNDONGTNGzZo107Rp05zlZw683r9/f5n1lLeuJUuW6KabblKzZs3KbavE3r17NWTIEM2fP9+5UxERESFJeuyxxxQTE6P4+Hj96le/0ooVKyrdlqrGGBISoptuuklr167V6tWr/bbVvn17JSQkKDg4WFdffbUefvhhJ6mo6npK6hYVFenUqVN+6/o7fiXLEbiqHqv27durS5cuevTRR8ttq3fv3rrooosUFxenOXPmaMWKFTp9+nSlx7gizZo108CBAzV69OizihcXltrop4HGdOa5iX7aMNXFfnomzqeQareflncNXVlM1X0+PeukorCwUHfddZdatWql559/3lmelJTk8yaoLVu2OBfoVfHQQw/p4MGDWrRokYKD/y+sHTt2OBlUfHx8mfWcPHlSmZmZzuNDJZYsWVLm0afSbUnSoUOHdPPNN2vq1KlKSUlx6sbGxqpVq1Y+n6/Ko09S2X2xZ88eFRYWqkOHDn7rFxcXKzMzs0ptn7lvqrovpJ+2p0WLFuUeI3/Hr3379k5yhcB06tRJOTk5Pm8+K+9v5Fz6g5lVeowrc+Z6k5KStGvXLhUUFDjlW7duPau/adQ/tdFPz1ZF5yb6acNUF/tpReulnzZMtdlPy7uGLq3Gz6dVHn3xv8aMGWP9+/e306dP+yx///33rV27drZnzx77/vvvLTk52Xn7k9lPbyrKy8szSZaVlWV5eXlO2RNPPGFXXHGF5eTkVLr+kjceLViwwPLz8+3xxx8v88ajL7/80qKioszj8VTYVnZ2tnXt2rXcgay//vWvbdCgQZabm2vfffedXX755bZgwQKfOhW9/WnVqlXOWyBK3v504sQJW7RokZ04ccIKCgrs7bfftvDwcNu6davfGFauXOm8IWDTpk3WunVre/PNN6u8L840ceJEZ3s2bNhg0dHR9vXXX5vZT28EaNKkiW3atMmys7Pthhtu4O1P1eyOO+6w8ePH26lTp2z58uXOWyDeeust27dvn5mZZWRk2JVXXunz1pAzbd++3TZv3myFhYV27Ngxu/vuu23gwIFOeUXHuLQPPvjAvv76aysuLraDBw9a//797c4773TKe/ToYdOnT7f8/Hz7/e9/z9tKGoja6KcFBQWWl5dnI0eOtKefftry8vKssLDQb1uVnZvopw1TXeunnE/hT23007O5hq7p8+lZJRVZWVkmycLDwy0qKsqZ1q1bZ2ZmM2fOtKZNm1pMTIxNmjTJiouLnc+2a9fOJPlMThCShYWF+bS5aNGicuPYuHGjdenSxcLDw613795l3nY0ceJEGz58eKXbs3DhQpPks96oqCin3Ov12r333mtut9tatWplTzzxhM82mflPKsx+OsF06NDBIiIibMiQIc4rYz0ej/Xr18+io6PN7XZb9+7dbenSpc7n1q1b5xPDhAkTLC4uzqKioqxz5842Z86cKu+LZ555xuctVKdOnbLhw4dbVFSUtWnTxhYvXuzT1oIFC6xVq1bmcrksNTXV79uqcO4OHz5sAwcOtIiICOvUqZN9+OGHZvbTcWrdurVFRkZa27ZtbcKECT5Jd1JSkvP3sHbtWuvUqZNFRkZa8+bNbeTIkfbDDz84dSs6xvv27bOoqCjnRPanP/3J2rdvb5GRkdayZUsbN26czyuTMzIy7JprrrHw8HC74oor7Msvv6zJ3YM6ojb66ZNPPlnm/wclX9iU7qdmFZ+b6KcNU13rp5xP4U9t9NOKrqFr+3wa9L8BAQAAAMA5CeiVsgAAAABAUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAAJCUgEAAAAgICQVAAAAAALy/wGorkmG2bQ4RQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxUAAADtCAYAAADJNFjjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA9YElEQVR4nO3deVgUV7o/8C+LAt1Nd6OCgCyKCwqixonLYzRq4hqXaEYzRjEu4xiTaJLxqpH5jUtwIzPqJPEO12viSEbRxMnVqIlrcCZmjEayqKAREUWJihtL2wIty/v7w9t1aWg2m7bB/n6ep5+Hqjp16lTVobrePqdOuYiIgIiIiIiI6CG5OroARERERETUuDGoICIiIiIimzCoICIiIiIimzCoICIiIiIimzCoICIiIiIimzCoICIiIiIimzCoICIiIiIimzCoICIiIiIimzCoICIiIiIimzCoICIiIiIimzCoIHJyJpMJ06ZNQ1BQEHQ6HQYMGICUlBSLNHFxcfD19UWzZs2wYMECiAgAoKSkBL/+9a/RqlUruLi4IDs722K9uXPnIiwsDN7e3njyySdx5MiRasuSnJyMrl27QqVSoX///rh8+bKyLDIyEhqNRvm4urpizZo1VvNJS0vDyJEj0aJFC/j6+iI6Ohq5ubkWaTZt2oT27dtDrVajU6dOOH/+PFauXKnk7+HhgSZNmijTs2bNAgDs27cP7dq1g1qtxvPPP2+R78yZMxEQEACtVouoqCh88cUXVe5rQkIC3N3dLfbpypUrtToWFRUWFiI6Ohre3t4ICQnBtm3bKm0rKCgIWq0W06ZNw/3796vM62EdO3YMAwcOhI+PDwICAjB79myL7dy6dQsjRoyASqVCeHg4kpKSlGVHjhzB008/DbVajWHDhlnkW5tzWVFV9fWbb76xON4qlQqurq64deuW1XwSEhLQrVs3eHt7IywsDOvXr7dYbjAY8Nvf/hbNmjWDXq/HxIkTAVjWVVdXV3h5eSnT33zzDcrKyvDWW29Br9ejZcuW+Mtf/qLkefHiRfTo0QM+Pj5o1qwZxowZg+vXr1e5rwMGDICnp6eS//Dhw2td/oqqq3M11TEiIggROTWj0SixsbGSlZUlJSUlsmbNGgkLC1OWf/nllxISEiIZGRly7do16dSpk2zcuFFERIqLi+W9996TY8eOCQC5fv26Rd5Lly6V9PR0KS0tle3bt4uPj48YDAar5SgqKpKgoCDZuHGjFBYWyoIFC6Rfv35W096+fVuaNm0q6enpVpd/99138vHHH0teXp4YjUYZP368TJ8+XVm+e/duiYyMlNTUVCkrK5MLFy7InTt3LPJYtWqVTJkyxWLejRs3RK/Xy969e8VoNMqkSZNk8uTJyvKff/5ZioqKRETkxIkTotPpJCcnx2oZN23aJEOHDrX5WIiIzJ8/X4YPHy75+fly9OhR0el0kpaWJiIip0+fFh8fH0lOTpa8vDwZMGCALFq0qMq8Hta+fftk586dYjQa5fbt29KvXz+JjY1Vlo8fP15mzJgh9+7dk507d4qPj49ybJKTkyUxMVGWLVtW6ZjUdC4rqq6+VvT+++9L3759q8xr/fr1cuzYMSkuLpbU1FTx8/OTr7/+Wln+wgsvyBtvvCF5eXly//59+fHHHyvlER4eLv/85z8t5v31r3+VJ554Qm7cuCHnzp0Tf39/SUpKEhERg8EgFy9elLKyMjGZTPL222/Liy++WGUZ+/fvL9u2bXuo8pdXU52rro4REYmIMKggIgsmk0lcXFzk9u3bIiIyYcIEiYuLU5Zv3LhRBg4cWGk9a0FFRYGBgfL9999bXbZ//37p2LGjMm00GsXLy0syMzMrpY2Pj5fevXvXan9ERA4ePCidO3dWpnv06CEHDx6sdh1rQcX69etl2LBhyvSFCxfE09NTCSTKS05OFg8PDzlz5ozV/KsLKupyLERE/P395fjx48r05MmT5Z133hERkYULF8qsWbOUZUlJSdKmTRur+dSnDRs2yMiRI0VE5O7du9K0aVO5du2asrxfv37y8ccfW6xT3TExq3guK6ptfRUR6dmzp6xfv77GfTGbOHGirF69WkREUlNTJTQ0VEpKSqpdx1pQ0bt3b/nkk0+U6UWLFsm0adMqrWsymSQmJkZ69OhRZf7VBRXVlb+imupcdXWMiEhEhN2fiMjCsWPH0LJlSzRv3hwAcPbsWURFRSnLu3btijNnztQ538zMTOTk5KBdu3ZWl1fcjlqtRtu2bXH27NlKaRMTEzFp0qRab/vbb79FZGQkAKC0tBQ//fQTUlJSEBQUhDZt2iA2NlbpIlOdimVs27Yt3N3dcfHiRWXea6+9Bi8vL/To0QPDhg1DREQEAODf//439Hq9RX5Hjx5F8+bNERERYdE1paZjERcXh5EjRwIAcnNzkZ2dXeU5snb+Ll26hMLCwhr31xblj3l6ejp0Oh0CAgKslvFh87WmtvX1woULOHnyJMaPH1+r7ZaWluLEiRPKtr///nt06NAB0dHRaN68OXr27IlvvvmmVnnVpox6vR5eXl5YvXo1/uM//kOZ/9prr+G1116zSDtnzhz4+vpi8ODBOH36dK3KDwBdunTB1q1brZapfJ2rqY4REQGAu6MLQEQNR35+Pl555RWsWLFCmWc0GqHVapVprVYLo9FYp3yLi4sxZcoUzJ8/HzqdzmqaitupaluZmZk4ceIEduzYUattnzx5Eh988IHyPMeNGzdQUlKCpKQkpKamwmAwYPjw4QgODsa0adOqzctoNMLX17faMsbHx2PdunX45z//afFsSt++fZGXl6dM9+/fHykpKQgJCUFycjLGjh2Lli1bYuzYsTUei4ULF1qUyc3NDSqVympaa+fPPN/Ly6va/X1Y+/btw759+3Dq1CmrZTCXo/zxqI2K59Ka2tbXxMREDBs2DM2aNavVtv/4xz+iVatWGDp0KADg6tWrOHToED766CP8/e9/x549ezBmzBikp6fXmGdtypiXlweDwYAPP/wQrVu3VubHx8dbpPvTn/6EiIgIuLm5Yd26dXjuuedw7tw5aDSaassPwCIAqa7O1VTHiIgAPqhNRP+rqKgIY8aMwYgRIzB9+nRlvkajgcFgUKYNBkOlG5bqiAimTp0KPz8/LF26VJlf/mHWK1euVNpOVdvaunUrBg0aBD8/vyrzMrt06RJGjRqFjRs3Kr/Qmm+k3377bej1eoSEhOD111/H3r17a9yX2pbRzc0NgwYNQlJSEg4cOGA1rzZt2qB169ZwdXVFr1698MYbb2Dnzp112o45bWlpKQoKCqymtXb+zPOrU/6h5vIP/9YkOTkZ06ZNw+eff46WLVvWeX+qYu1cWitjbevr1q1bLVq7qtvf9evXY8eOHfjss8/g4uIC4EE9atOmDX7729+iSZMmeOGFF9C2bVscO3asxn2pbRm1Wi2mTJmCsWPHVtmS1rNnT2g0Gnh5eWHBggXQaDQ4ceJEjeWvqUzly1VTHSMiAhhUEBEejOI0YcIEBAYGYvXq1RbLIiIiLH5xP3XqVLXdTyqaM2cOrl27hi1btsDV9f8uOWfOnFF+BQ0JCam0nXv37iEjI0PpPmRW8WbQWl4AkJ2djcGDB2PRokUYM2aMktbHxweBgYEW69em6xNQ+VhcvHgRJSUlCAsLs5q+rKwMGRkZtcq7/LGp7bEAHuyPv79/lefI2vlr06ZNja0U/fr1U47pvn37arUPP//8M0aPHo2//e1v6N27tzK/ffv2yM/PtxgdrC71qKpzaa2MtamvycnJuH79OkaNGlXj/n766adYsWIFDhw4gBYtWijzO3fuXKmcD1uPqjsWZWVluH79usUNfXXK1yOg6vLXVKbyda6mOkZEBICjPxGRyNSpU2XIkCFy//79Ssu++OILCQ0NlYsXL8r169clMjLSYjSdoqIiKSwsFACSmZkphYWFyrLFixdLt27dJD8/v8YymEef2bRpkxQVFcnChQsrjXj0008/iVqtFqPRWG1eeXl50qVLlyofJP3DH/4gI0aMEIPBIFevXpXOnTvLpk2bLNJUN/rT/v375d69exIdHa2M/nT37l3ZsmWL3L17V4qLi+Wzzz4TT09POX36tNUy7Nu3T27evCkiIj/88IO0atVKPv3001ofi/LmzZun7M+xY8dEp9PJuXPnROTB6E/NmjWTH374QfLy8uSZZ56xy+hPV65ckeDgYElISLC6fNy4cTJz5kwpKCiQXbt2WYz+VFpaKoWFhbJhwwYZPHiwFBYWKnWxpnNZUU31VUTkzTffrHRurTlw4ID4+vrKqVOnKi27f/++tG7dWhISEqSkpER27dolzZs3rzSKmLUHtf/zP/9TfvWrX8nNmzclLS1NAgIC5KuvvhIRka+//lp++OEHKSkpkZycHJk8ebL07NnTavlyc3Pl4MGDUlRUJCaTSdauXSstW7ZU/t+qK39FNdW56uoYEZEIR38icnqZmZkCQDw9PUWtViufI0eOKGlWrlwpzZs3F71eL/Pnz5eysjJlWWhoqACw+JgBEA8PD4t8t2zZUmVZTpw4IVFRUeLp6Sn9+vWrNNrRvHnzZOLEiTXuU0JCggCw2K5arVaWm0wmmTFjhmi1WgkMDJTFixdb7JOI9aBC5MGQpWFhYeLl5SWjRo1SboyNRqMMHDhQdDqdaLVa6d69u+zYsUNZ78iRIxZlmDt3rvj6+oparZYOHTrIBx98UOtjsWLFCotRqAoKCmTixImiVqslKChIEhMTLfLatGmTBAYGikajkSlTplgdrcpWS5cuFRcXF4vjHRERoSy/efOmDB8+XLy8vKR9+/Zy6NAhZdk///nPSnXIfOxrOpfWVFdfS0pKxN/fv8bRv0REBgwYIO7u7hbbfeWVV5Tlp06dkieffFLUarV069bN6nCt1oKK0tJSefPNN0Wn04mvr6+sWbNGWfbll19Kx44dRa1Wi6+vr4wbN87i3L/yyitKGW7evCm/+tWvRK1Wi4+PjwwcOFB++OGHWpc/IiLC4v+xujpXUx0jInIRqWV7LRERERERkRV8poKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzCoIKIiIiIiGzi7ugCNHRlZWW4du0avL294eLi4ujiEBEREZGTExHcvXsXgYGBcHVtGG0EDCpqcO3aNQQHBzu6GEREREREFrKyshAUFOToYgBgUFEjb29vAA9OmlardXBpiIiIiMjZGQwGBAcHK/epDQGDihqYuzxptVoGFURERETUYDSkrvkMKuzg6tWrSElJgclkcnRRGj0XFxf4+/uje/fucHdndSUiIiJqiHiXVo+KiooQExODb775xtFFeexotVq8++676NGjh6OLQkREREQVMKioR2vXrkVycjJiY2PRv39/qFSqBtUs1RiVlpYiIyMD77//Pt566y3s2bMHzZo1c3SxiIiIiKichjEG1WOgpKQEBw8eRHR0NJ577jmo1WoGFPXAzc0NHTp0wIoVK1BcXIx//etfji4SNWIffPAB7ty5Y5e8jUYjjEajXfImIiJq6BhU1JM7d+7AaDSiS5cuji7KY0mv1yM0NBSZmZmOLgo1Urdu3cKbb76JOXPm2CV/rVYLHx8fu+RNRETU0DGoqCelpaUAwIeJ7ahJkyYoKSlxdDGokSorKwMAFBYW2iV/EWH9JCIip8Wgwo7GjRuH3//+9xbzVq9ejf79+yM0NBTHjx+3WDZr1iwsXbrUal4DBgyAp6cnvL29odfr0bdvXyQkJFikad26NVQqFTQajfLZuHEjgAc3VG+99Rb0ej1atmyJv/zlL8p6iYmJFuuYt1MbCQkJcHd3V9bt1KkTPv/8c4vlw4YNs1gnOzvbomuYed/Kl2HJkiW12j4REREROR5/VrejdevWISoqCpMnT0b37t1x5coVrFq1CkePHq10o10bCQkJmDBhAvLz83Ho0CG88cYbOHfuHOLi4pQ0hw8fRu/evSutu379ehw5cgTnz59Hbm4uBgwYgK5du+KZZ57BpEmTMGnSJCXt3LlzcevWrVqXa9CgQdi/fz8A4MCBA3jhhRdw7do16HS6Ou8bERERETU+bKmwo4CAACxfvhwzZ85EaWkp5syZgzfeeAMdO3a0KV+dTodx48Zh/fr1WLt2ba0CgM2bN+Ptt9+Gn58fwsPD8bvf/Q5btmyplK6srAyffPIJoqOjH6psQ4cOhVqt5rMPRERERE6EQYWdvfLKK/D09MQLL7yA9PR0xMTE1Fvezz33HMrKypCcnFxj2rNnzyIqKkqZ7tq1K86cOVMp3eHDh1FaWopBgwbVuTwigv3796OsrAzt2rWr8/pERERE1Dix+5Odubi4YMOGDYiMjMThw4fRtGlTZdngwYPh5uamTBcWFtYp6HB3d0eLFi2Qm5tbZZ779+9H7969YTQaodVqlflardbq8JeJiYmYMGGCRR41SUpKgl6vR3FxMYqKivDee+9BrVZXWm4mIpXy+O1vf4tZs2Yp0xs2bMCLL75Y6zIQERERkeMwqHgEIiIiAACdOnWymH/o0CGL5x/K31TXRklJCW7fvm0xjGXFPM00Gg0MBoMybTAYoNFoLNIUFRVhx44dOHToUJ3K8eyzzyrPVPz888949tln0bp1a4waNarScuDBg9oBAQEWeWzcuJHPVBARERE1Uuz+1Ijt27cPrq6u6NmzZ41pIyIikJKSokyfOnUKkZGRFmn27NkDPz+/WuVXlU6dOuGpp57CgQMHHjoPInuw1kJGRERE9YNBRSNkMBiwY8cOzJo1C3PnzkWLFi1qXCc6Ohp//vOfcevWLZw/fx4fffSRxYhPwIOuTxXn1dWFCxdw9OjRSq0yRI5mfk8FgwsiIqL65/Cg4tatWxgxYgRUKhXCw8ORlJRkNV1hYSGio6Ph7e2NkJAQbNu2TVl2/PhxdOnSBXq9Hn5+fpgyZYrF8wIZGRl46qmnoFKp0L17d5w6dcru+2WrxMTESi0JU6dOhUajQXBwMNasWYPly5dbDCcLAM8884zF+x6WLVsGAHj11VfRt29ftG/fHn379sW8efPw7LPPKuvl5uZi3759VoOK4cOHY+XKlcq0RqPBN998o0x/9dVX0Gg0UKvV6N+/P1566aU6d+Uy75v5M2PGjDqtT1QTBhNERET24yIO/qZ98cUXodPp8P777+PgwYOYPn06MjIyLJ4TAIAFCxYgNTUVn3zyCVJTU/Hcc8/hxIkT6NChA27fvg2TyYRWrVqhoKAAr7zyCgIDA/Huu+8CAHr27IkRI0bg7bffxkcffYS1a9fi/PnztXr7tcFggE6nQ35+vsWDzhVdu3YNo0ePRnx8vE3dh6hqEydORLdu3bBgwQJHF4UaoaysLISEhOD555+3eEFjfTG/0JHBCxER2Vtt708fJYe2VBiNRuzatQuxsbFQqVQYM2YMOnfujD179lRKu3nzZixZsgRarRZ9+vTB6NGj8cknnwAAWrRogVatWgF48IXu4uKCS5cuAQDS0tKQlpaGmJgYeHp6Yvbs2SgtLcW3335rtUwmkwkGg8HiQ0SNH2/2iYiI7MehQUV6ejp0Op3FSEDW3p+Qm5uL7Ozsat+zcOXKFej1emg0GuzcuROzZ88G8OD9DOHh4RZDuXbp0sXqOxoAYNWqVdDpdMonODi4XvaVbMebQrIF6w8REZH9OLylomKTjbX3JxiNRri5uUGlUlWZLiQkBHl5ebhx4wYWLFigBCq13YZZTEwM8vPzlU9WVlat9sX8XoacnJxapae6ERHk5ORYvP+CqC7MD2oTERFR/XNoUFHx3QmA9fcnaDQalJaWoqCgoNp0AODn54fhw4fj5ZdfrtM2zDw8PKDVai0+taHT6dChQwd88cUXvHmxg2PHjuHOnTt8XoUeGlsqiIiI7MehL79r37498vPzkZ2dDX9/fwAP3p9QceQfHx8f+Pv7IyUlBb169VLSVRwdyaysrAwZGRkAHryfIS0tDcXFxWjSpAkA4PTp05g/f36978/vfvc7vP3223j55ZfRr18/qNVq5eFNejilpaW4cOECvvrqK/Ts2RPdu3d3dJGokWKwT0REZD8ODSo0Gg1Gjx6NJUuW4L333sOhQ4eQmpqqvIm5vOjoaCxbtgzbtm3DmTNnsHv3bnz33XcAgL1796Jt27bo0KEDsrOzsWjRIgwcOBAAEB4ejvDwcMTFxWHBggXYuHEj3Nzc0KdPn3rfn4EDB2LdunX4xz/+gf/5n/+ByWSq9204G1dXV7Rs2RJTp07Fyy+/DDc3N0cXiRoptlQQERHZj0ODCgCIj4/HlClT0Lx5cwQFBWH79u3w8fFBYmIiVq5cqTxQHRsbixkzZiAgIAA+Pj6Ij49HeHg4AODGjRuYPXs2bty4AZ1Oh+HDh+PPf/6zso2tW7diypQpWLlyJTp27IgdO3bUajjZh9GrVy+lNYWIGg62VBAREdmPw99T0dA1xHGAiajuzp07h06dOvE9FURE1Og1xPtTh79Rm4joUWBLBRERkf0wqCAip8AWBCIiIvthUEFEToEtFURERPbDoIKInAJbKoiIiOyHQQUROQW2VBAREdkPgwoicgpsqSAiIrIfBhVE5BTYUkFERGQ/DCqIyCmwpYKIiMh+GFQQkVNgSwUREZH9MKggIqdgbqlgiwUREVH9Y1BBRE6BwQQREZH9MKggIqfA7k9ERET2w6CCiJwCWyqIiIjsh0EFETkFtlQQERHZD4MKInIKbKkgIiKyHwYVROQU2FJBRERkPwwqiMgpsKWCiIjIfhweVNy6dQsjRoyASqVCeHg4kpKSrKYrLCxEdHQ0vL29ERISgm3btinLzp49i8GDB0On06Fjx46V1m3dujVUKhU0Gg00Gg1mzZplt/0hooaJLRVERET24+7oArz++usIDAzE7du3cfDgQYwfPx4ZGRnw8fGxSLdkyRLk5OTg6tWrSE1NxXPPPYdf/epX6NChA5o0aYKJEyfiN7/5DVavXm11O4cPH0bv3r0fxS4RUQPElgoiIiL7cWhLhdFoxK5duxAbGwuVSoUxY8agc+fO2LNnT6W0mzdvxpIlS6DVatGnTx+MHj0an3zyCQCgffv2mDZtGtq1a2dzmUwmEwwGg8WHiBo/c0uFi4uLg0tCRET0+HFoUJGeng6dToeAgABlXteuXXHmzBmLdLm5ucjOzkZUVFS16aozZswYtGzZEmPHjsXly5erTLdq1SrodDrlExwcXIc9IqKGqrS0FADg5ubm4JIQERE9fhzeUqHVai3mabVaGI3GSunc3NygUqmqTVeVrVu3IjMzE+np6QgJCcGYMWOq7AoRExOD/Px85ZOVlVXHvSKihsgcVLi6OvxRMiIioseOQ5+p0Gg0lboXGQwGaDSaSulKS0tRUFCgBBbW0lWlT58+AABPT0+sXbsWOp0Oly5dQlhYWKW0Hh4e8PDweJjdIaIGrKSkBABbKoiIiOzBoT/ZtW/fHvn5+cjOzlbmnTp1CpGRkRbpfHx84O/vj5SUlGrT1YaLiwv7VBM5IXZ/IiIish+HBhUajQajR4/GkiVLUFhYiN27dyM1NRWjRo2qlDY6OhrLli3D3bt3cfz4cezevRu/+c1vADwY1aWoqAj379+3+BsArly5gmPHjqG4uBj37t3D/PnzERoaitatWz/KXSUiB3tULRUcZYqIiJyRwzsXx8fHIysrC82bN8e8efOwfft2+Pj4IDEx0aIlIjY2Vnmoe/z48YiPj0d4eDgA4PLly/Dy8sLQoUNx/vx5eHl5YciQIQCAu3fvYubMmdDr9WjdujUuXLiAXbt2sV81kZN5VC0VfB8GERE5I4e/p8LX1xd79+6tNH/SpEmYNGmSMu3l5YXExESrebRu3brKXwcjIyMtuk0RkXN6VC0VJSUl7GJFREROhz/XE5FTeFQtFcXFxXbNn4iIqCFiUEFETsHcUmHvgRrM2yEiInImDCqIyCmYWyrsjS0VRETkjBhUEJFTMLcg2Ht0JrZUEBGRM2JQQUROwdxSYe+ggi0VRETkjBhUEJFTeFTdn9hSQUREzohBBRE5BfPN/sWLFzF27FgYjcZ6yXfnzp0YPnx4pe0QERE5E4e/p4KI6FEwt1QkJSUBAM6dO4cnn3zS5nzfeOMN/PLLL5W2Q0RE5EzYUkFETqFiC0J+fn695KtSqSymGVQQEZEzYlBBRE6h4s1+Xl6eXfJlUEFERM6IQQUROYWKLRUMKoiIiOoPgwoicgoVb/Zzc3PrJd+KwUpZWVm95EtERNSYMKggIqdgr5aK+/fvW0yzpYKIiJwRgwoicgr2eqbi3r171W6HiIjIGTCoICKnYI/uTyLCoIKIiAgMKojoMTd//nz8+c9/tkv3p7t371aaV1JSggEDBuDIkSM2509ERNRYMKggosfa6tWrsWDBgkoBQH0EFdZaO/77v/8bX3/9NRYuXGhz/kRERI2Fw4OKW7duYcSIEVCpVAgPD1fedltRYWEhoqOj4e3tjZCQEGzbts1ieUJCAoKCgqDVajFt2jSLhyczMjLw1FNPQaVSoXv37jh16pRd94mIGobygcS1a9csluXm5qK4uBgGg+Gh8s7NzbUamFy5cgUA4O/v/1D5EhERNUYODypef/11BAYG4vbt23j33Xcxfvx4q7/+LVmyBDk5Obh69So++eQTvPrqqzh//jwAICUlBXPnzsXnn3+OrKwsZGZmYvny5cq6L730EoYMGYKcnBxMnz4dY8eOrdQVgogeH2VlZSgrK0NycrIyr+J1JS8vD9HR0dDpdPj555/rlP/IkSPRrFkzHD16tNIyo9EIADh27JhSFiIiosedi4iIozZuNBrRvHlzZGZmIiAgAADw9NNPY8aMGXj55Zct0gYEBODzzz9Hr169AAAvv/wy2rVrh8WLFyMmJgZ5eXn4r//6LwDA4cOHMWPGDFy8eBFpaWno2bMnbt26haZNmwIAQkNDsXnzZjz99NOVymQymWAymZRpg8GA4OBgrF27Fl5eXtXuj4uLi/K3Aw+rUyl/zB+1x/kclz+utvxdn3nV9HdBQQGWLl2KX375RVn2hz/8AStXrgQAtG7dGpmZmQCApk2bwtXVFUVFRUra999/H2vXrsXly5cBPAgc+vfvjw0bNiA9PR0AMHToUIwcORJz5syx2Ec3NzflAe2goCClDOW3DwDdu3fHpEmTUFpaivT0dLi4uMDf3x9arVa5vogI3Nzcqr2eVFf3XFxcKi2veE7KL7eWl5ubmxKYlefq6lrt+bVFTWWq7+0RmdXmWs56Rw1NYWEh5s6di/z8fGi1WkcXB4CDg4qffvoJQ4cOxc2bN5V5c+bMgUqlwrvvvqvMy83NRbNmzXDv3j2oVCoAwJo1a3DixAl8+umneP755zF06FC89tprAIA7d+6gRYsWKCgowP79+7Fq1SqcOHFCyW/UqFF47rnn8Oqrr1Yq09KlS/HOO+9Umm/ty9QaEanzxcfaOrXdVnUe94tgQ7ipfxyPcVU3d7W96WtoAgICcP36dWi1WqWrU6tWrXD16tV620bLli1x48YNAIBarYZGo1GmrfHw8LD48QIAmjRpogQEpaWlynXBXMdqc40QEZSVlVmsZ55f22DPvH3zNc/8sRZkPMz1rjo15SUi9b5Nql8N7dxUda2qWI+qK3djut6R8zBf7xtSUOHuyI0bjcZKB0Kr1Vbqp2w0GuHm5qYEFOZ05m4GFfMx/200GqvchnndimJiYjB37lxl2txSkZub22BOGlFDVd0v6vb+G3jwzojDhw+jW7duWL58ObZs2YLf/e53WLZsmcWzE35+fkpQsWfPHvzjH//A2LFjMXr0aPzyyy9o2rQp9Ho9PDw8UFhYiOzsbGi1Wuj1eixfvhwdOnTA4cOHsXHjRmg0GotWzHv37mHUqFH46quvcPv2bXz33XfIzc2Fi4sLnn76aXh6egIAiouLATwIJqwdx4Z2c0ZERA2HwWCATqdzdDEsODSo0Gg0lR6SNBgM0Gg0ldKVlpaioKBACSzKp6uYj/lvjUZT622YeXh4wMPDw7YdI3JS9uoaU1sqlQq/+c1vAAB///vf8eqrr6JHjx7YsGEDsrOzlXS+vr7K3wMHDsTIkSOV6ZCQkEp5hoWFKdNLly4F8H8PZBcWFipdK8vnkZKSgnv37qFt27ZWy2otmDBjQEFERI2NQx/Ubt++PfLz8y2+7E+dOoXIyEiLdD4+PvD390dKSorVdBEREZWWtWnTBl5eXoiIiEBaWpryqyAAnD59utI2iOjx4uLigj59+qBJkybw8fGxWNaiRQvlb7Va/VD5Dxo0CMCDl925u1v+PqPT6eDv719lQEFERPS4cWhQodFoMHr0aCxZsgSFhYXYvXs3UlNTMWrUqEppo6OjsWzZMty9exfHjx/H7t27lV8kJ06ciO3bt+PHH39Efn4+VqxYgejoaABAeHg4wsPDERcXB5PJhPj4eLi5uaFPnz6PdF+JyHHMXY7M/Pz8bM4zMDBQ+dvV1fJSWnF7REREjzuHDykbHx+PrKwsNG/eHPPmzcP27dvh4+ODxMREi9aE2NhY6HQ6BAQEYPz48YiPj0d4eDgAICoqCmvWrMGoUaMQFBSE4OBg/L//9/+Udbdu3Yr9+/dDr9fjww8/xI4dOyr9skhEj6+KI7e1atUK+/fvx/Xr1x86z8DAQBw4cACXLl2qFFTUNFIcERHR48ahoz81BuYHYRrS0/VEVDeDBg2yeLHm2rVr8fvf/77e8n/iiSdw8uRJZXrTpk2YOnVqveVPRERUXkO8P3V4SwURkb1V7I7k5uZWr/mz+xMRETk7BhVE9Nir2B2pYhBgK3Z/IiIiZ8eggogeexWHkLZ3S4W3t3e95k9ERNTQMaggosdex44dLabtHVR06tSpXvMnIiJq6BhUENFjr2XLlhbT9u7+VP7lekRERM6AQQURPfYqvr26vlsqKubHIauJiMjZMKggIqdT3y0VLi4u9ZofERFRY8Oggoicjr2fqSAiInI2/CYkosdecHCwxbS9uz8RERE5GwYVRPTYe/rpp3H69Gn06dMHgP0e1F63bh0uXbpUr3kTERE1BgwqiMgpREVFoaysDID9uj/5+PigdevW9Zo3ERFRY8CggoichogAsF9QwQe2iYjIWTGoICKnUVpaCsB+3Z8YVBARkbNiUEFETsPe3Z+IiIicFb8JichpmIMKtlQQERHVLwYVROQ07NVSwSFliYjI2TGoICKnwQe1iYiI7MOhQUVycjK6du0KlUqF/v374/Lly1WmzcjIwFNPPQWVSoXu3bvj1KlTyrLPPvsMvXr1goeHB2bNmmWxXmZmJlxcXKDRaJRPYmKi3faJiBouez+oTURE5Kwc9k1oMpnwwgsv4M0330ROTg569+6NyZMnV5n+pZdewpAhQ5CTk4Pp06dj7NixKCkpAQA0a9YMCxYswIwZM6yu6+HhAaPRqHwmTZpkl30ioobN3g9qs6WCiIiclcOCin/961/QaDSYPn06PD09sXjxYnz//fdWWyvS0tKQlpaGmJgYeHp6Yvbs2SgtLcW3334LAHjmmWfw61//Gr6+vjaXy2QywWAwWHyI6PHAoIKIiMg+HBZUnD17FlFRUcq0Wq1G27ZtcfbsWatpw8PD0bRpU2Vely5dcObMmVpt6/79+wgMDERoaCjeeustFBQUVJl21apV0Ol0yic4OLgOe0VEDZm9R38iIiJyVg77JjQajdBqtRbztFotjEajTWkratGiBb7//ntkZWXhyJEjOHnyJGJiYqpMHxMTg/z8fOWTlZVVyz0iooaOD2oTERHZh92CiiFDhsDT09PqZ/ny5dBoNJW6FhkMBmg0mkp51SWttXW7d+8ONzc3hIaGIi4uDjt37qwyvYeHB7RarcWHiB4PfFCbiIjIPtztlfHBgwerXX7gwAFs2LBBmb537x4yMjIQERFRKW1ERATS0tJQXFyMJk2aAABOnz6N+fPn17lc/PIncl58poKIiMg+HHaHPWDAABiNRiQkJMBkMmH58uV48sknERoaWilteHg4wsPDERcXB5PJhPj4eLi5uaFPnz4AHvz6WFRUhJKSEou/AeDEiRNIT0+HiODatWuIiYnByJEjH+m+ElHDwKCCiIjIPhwWVHh4eGDHjh1Yu3Yt9Ho9jh49is2bNyvLZ82aZfHOia1bt2L//v3Q6/X48MMPsWPHDri7P2ho2bx5M7y8vLBixQp89NFH8PLywvLlywEAFy5cwODBg6HRaNCrVy9ERETgT3/606PdWSJqEPigNhERkX24iPnJRbLKYDBAp9MhPz+fz1cQNXKBgYG4fv06zpw5Y7Wr5cN69dVXsX79euzYsQNjx46tt3yJiIisaYj3p/x5jYichvlB7fru/lTf+RERETU2DCqIyGnYu/sTn6kgIiJnxaCCiJyGvR/UJiIiclb8JiQip8GWCiIiIvtgUEFEToNDyhIREdkHgwoichrs/kRERGQf/CYkIqfB7k9ERET2waCCiJwGWyqIiIjsg9+EROQ02FJBRERkHwwqiMhp8EFtIiIi+2BQQUROg92fiIiI7IPfhETkNNj9iYiIyD4YVBCR06nvlgpzfiJSr/kSERE1FgwqiMjp2Kv7U2lpab3mS0RE1FgwqCAip2Ov7k/m7lVERETOhkEFETkdewUV7P5ERETOikEFETmd+n6gmt2fiIjI2Tk0qEhOTkbXrl2hUqnQv39/XL58ucq0GRkZeOqpp6BSqdC9e3ecOnVKWbZ+/XqEhYVBq9UiJCQEq1atslh33759aNeuHdRqNZ5//nnk5ubabZ+IyPmw+xMRETk7hwUVJpMJL7zwAt58803k5OSgd+/emDx5cpXpX3rpJQwZMgQ5OTmYPn06xo4di5KSEgDAsGHD8OOPP8JgMOC7775DYmIi9u3bBwC4efMmJk6ciHXr1uHmzZvw9vbGm2+++Uj2kYicA4MKIiJydg4LKv71r39Bo9Fg+vTp8PT0xOLFi/H9999bba1IS0tDWloaYmJi4OnpidmzZ6O0tBTffvstAKB169bQ6/VKehcXF1y6dAkAsHPnTvTu3RvDhw+HWq3GO++8g3/84x8wmUxWy2UymWAwGCw+RETVYVBBRETOzmFBxdmzZxEVFaVMq9VqtG3bFmfPnrWaNjw8HE2bNlXmdenSBWfOnFGmt27dCm9vbwQGBqKoqAjjxo2zup22bdvC3d0dFy9etFquVatWQafTKZ/g4GCb95WIHm8MKoiIyNk5LKgwGo3QarUW87RaLYxG40OlnThxIu7evYuUlBRER0dDpVLVeTsAEBMTg/z8fOWTlZX1UPtHRM6DQQURETk7uwUVQ4YMgaenp9XP8uXLodFoKnUtMhgM0Gg0lfKqS9rOnTtDrVZjxYoVdV4XADw8PKDVai0+RETVYVBBRETOzm5BxcGDB1FUVGT188c//hERERFISUlR0t+7dw8ZGRmIiIiolFdERATS0tJQXFyszDt9+jQiIyOtbrusrAwZGRnKuuW3c/HiRZSUlCAsLKy+dpWInByDCiIicnYO6/40YMAAGI1GJCQkwGQyYfny5XjyyScRGhpaKW14eDjCw8MRFxcHk8mE+Ph4uLm5oU+fPgCAzZs34+bNmxARnDx5En/9618xcOBAAMDYsWNx/PhxHDhwAAUFBViyZAnGjx8PDw+PR7q/RPT4cnNzA8CggoiInJfDggoPDw/s2LEDa9euhV6vx9GjR7F582Zl+axZszBr1ixleuvWrdi/fz/0ej0+/PBD7NixA+7u7gCAEydOICoqCt7e3hg3bhxmzpyprOvn54fExES89tpraNGiBfLz8/H+++8/2p0loscaX35HRETOzkVExNGFaMgMBgN0Oh3y8/P5fAVRI2d+k3Z9X/b+/ve/Y8qUKYiPj8err75ar3kTERFV1BDvTx36Rm0ioseBuaWCv9EQEZGzYlBBRGQjdn8iIiJnx6CCiMhGHP2JiIicHYMKIiIbMaggIiJnx6CCiMhGDCqIiMjZMaggIrIRgwoiInJ2DCqIiGzEoIKIiJwdgwoiIhsFBgYCAEJCQhxcEiIiIsdgUEFETqN37952ybdnz5748ccfMWHCBLvkT0RE1NC5O7oARESPyjfffIOSkhK75P3EE0/YJV8iIqLGgEEFETkNd3d3uLvzskdERFTf2P2JiIiIiIhswqCCiIiIiIhswqCCiIiIiIhsws7FNRARAIDBYHBwSYiIiIiI/u++1Hyf2hAwqKjBnTt3AADBwcEOLgkRERER0f+5c+cOdDqdo4sBgEFFjZo1awYAuHLlSoM5aUQVGQwGBAcHIysrC1qt1tHFIbKK9ZQaA9ZTagzy8/MREhKi3Kc2BAwqauDq+uCxE51Ox4sLNXharZb1lBo81lNqDFhPqTEw36c2BA2nJERERERE1CgxqCAiIiIiIpswqKiBh4cHlixZAg8PD0cXhahKrKfUGLCeUmPAekqNQUOspy7SkMaiIiIiIiKiRoctFUREREREZBMGFUREREREZBMGFUREREREZBMGFUREREREZBMGFUREREREZJM6BRUmkwnTpk1DUFAQdDodBgwYgJSUFGV5XFwcfH190axZMyxYsADmgaVKSkrw61//Gq1atYKLiwuys7Mt8p07dy7CwsLg7e2NJ598EkeOHKm2HMnJyejatStUKhX69++Py5cvK8siIyOh0WiUj6urK9asWWM1n7S0NIwcORItWrSAr68voqOjkZuba5Fm06ZNaN++PdRqNTp16oTz589j5cqVSv4eHh5o0qSJMj1r1iwAwL59+9CuXTuo1Wo8//zzFvnOnDkTAQEB0Gq1iIqKwhdffFHlviYkJMDd3d1in65cuVKrY1FRYWEhoqOj4e3tjZCQEGzbtq3StoKCgqDVajFt2jTcv3+/yryo7m7duoURI0ZApVIhPDwcSUlJAIBdu3ahY8eO0Ol0CAgIwNy5c1FaWlplPi4uLlCr1Up9WLlypbKspnNc3vHjx9GlSxfo9Xr4+flhypQpMBqNyvKMjAw89dRTUKlU6N69O06dOlUPR4EaukdRT+Pj49GtWze4u7sjLi6uxjJVd21iPXVODa2e8npK1jyKelrXe2i7Xk+lDoxGo8TGxkpWVpaUlJTImjVrJCwsTEREvvzySwkJCZGMjAy5du2adOrUSTZu3CgiIsXFxfLee+/JsWPHBIBcv37dIt+lS5dKenq6lJaWyvbt28XHx0cMBoPVMhQVFUlQUJBs3LhRCgsLZcGCBdKvXz+raW/fvi1NmzaV9PR0q8u/++47+fjjjyUvL0+MRqOMHz9epk+frizfvXu3REZGSmpqqpSVlcmFCxfkzp07FnmsWrVKpkyZYjHvxo0botfrZe/evWI0GmXSpEkyefJkZfnPP/8sRUVFIiJy4sQJ0el0kpOTY7WMmzZtkqFDh9p8LERE5s+fL8OHD5f8/Hw5evSo6HQ6SUtLExGR06dPi4+PjyQnJ0teXp4MGDBAFi1aVGVeVHfjx4+XGTNmyL1792Tnzp3i4+MjOTk58ssvv8jNmzdFRCQnJ0eeffZZiY+PrzIfa/9DZtWd44pu3bolv/zyi4iI3Lt3T6Kjo2XBggXK8h49esjSpUulsLBQ1q1bJ23atJHi4uKH3X1qJB5FPd25c6fs2bNHxo4dK6tWraq2PDVdm1hPnVNDq6e8npI1j6Ke1uUe2t7X0zoFFRWZTCZxcXGR27dvy4QJEyQuLk5ZtnHjRhk4cGDlDVZzYMwCAwPl+++/t7ps//790rFjR2XaaDSKl5eXZGZmVkobHx8vvXv3ru3uyMGDB6Vz587KdI8ePeTgwYPVrmMtqFi/fr0MGzZMmb5w4YJ4enoqgUR5ycnJ4uHhIWfOnLGaf3VBRV2OhYiIv7+/HD9+XJmePHmyvPPOOyIisnDhQpk1a5ayLCkpSdq0aWM1H6q7u3fvStOmTeXatWvKvH79+snHH39skc58cZk/f36VeVX3P1TdOa6O0WiUyZMny/jx40VE5Ny5c6LVasVkMilpQkJC5Ouvv64xL2q8HlU9NZsyZUqNN2vVXZtYT51TQ6yn5fF6SiKPvp6aVXcPbe/rqU3PVBw7dgwtW7ZE8+bNcfbsWURFRSnLunbtijNnztQ5z8zMTOTk5KBdu3ZWl1fcjlqtRtu2bXH27NlKaRMTEzFp0qRab/vbb79FZGQkAKC0tBQ//fQTUlJSEBQUhDZt2iA2Nlbp0lWdimVs27Yt3N3dcfHiRWXea6+9Bi8vL/To0QPDhg1DREQEAODf//439Hq9RX5Hjx5F8+bNERERgfXr11e5nYrHIi4uDiNHjgQA5ObmIjs7u8pzZO38Xbp0CYWFhTXuL9UsPT1daeY0K3/8//3vf0On06FZs2ZISUnB9OnTlXRdunTB1q1bLfLr3r07WrVqhalTp+LOnTsAaj7HV65cgV6vt+g+Z56n0Wiwc+dOzJ49G8CD+hAeHo6mTZtalONh/qep8XgU9bQmFetpddcm1lPn1BDrafl5vJ4S4Jh6WvEe+lFfTx86qMjPz8crr7yCFStWAACMRiO0Wq2yXKvVWvQnrI3i4mJMmTIF8+fPh06ns5qm4naq2lZmZiZOnDiBF198sVbbPnnyJD744AMsWrQIAHDjxg2UlJQgKSkJqamp+Prrr/Hpp58iISGhxrxqU8b4+HgYjUYcOnQI/fv3V+b37dsXeXl5ynT//v2RkpKCW7duYdOmTYiNjcXOnTtrtZ2FCxcqz2sYjUa4ublBpVJZTWvt/Jnnk+1qOld9+/ZFfn4+Ll68iFmzZsHHx0dJd/r0aUycOFGZPnLkCC5fvoyTJ0+ioKBAuRDVdI5DQkKQl5eHkJAQZbl53o0bN7BgwQLl4lfb/zN6vDyKelqTivW0umsT66lzaoj1tPw8Xk8JePT11No99KO+nj5UUFFUVIQxY8ZgxIgRyo5pNBoYDAYljcFggEajqXWeIoKpU6fCz88PS5cuVeaXf/D6ypUrlbZT1ba2bt2KQYMGwc/Pr8q8zC5duoRRo0Zh48aNSkuFl5cXAODtt9+GXq9HSEgIXn/9dezdu7fGfaltGd3c3DBo0CAkJSXhwIEDVvNq06YNWrduDVdXV/Tq1QtvvPGGElTUdjvmtKWlpSgoKLCa1tr5M88n29X2XLVp0wZRUVF46623qsyrX79+aNKkCXx9ffHBBx9g7969uH//fo3nuDp+fn4YPnw4Xn755TqVlx4vj6Ke2lqm8tcm1lPn1BDraXm8nhLwaOtpVffQNZWpvq+ndQ4qSkpKMGHCBAQGBmL16tXK/IiICIuRoE6dOqXcoNfGnDlzcO3aNWzZsgWurv9XrDNnzigRVEhISKXt3Lt3DxkZGUr3IbOtW7dW6vpUMS8AyM7OxuDBg7Fo0SKMGTNGSevj44PAwECL9WvT9QmofCwuXryIkpIShIWFWU1fVlaGjIyMWuVd/tjU9lgAD/bH39+/ynNk7fy1adNGCa7INu3bt0d+fr7FyGdV/Y88TH0QkRrPcU3KbzciIgJpaWkoLi5Wlp8+fbpO/9PU+DyKelpX1V2bWE+dU0Osp9Vtl/XUOT3KelrVPXRFdr+e1vrpi/81depUGTJkiNy/f99i/hdffCGhoaFy8eJFuX79ukRGRiqjP4k8GKmosLBQAEhmZqYUFhYqyxYvXizdunWT/Pz8GrdvHvFo06ZNUlRUJAsXLqw04tFPP/0karVajEZjtXnl5eVJly5dqnyQ9Q9/+IOMGDFCDAaDXL16VTp37iybNm2ySFPd6E/79+9XRoEwj/509+5d2bJli9y9e1eKi4vls88+E09PTzl9+rTVMuzbt08ZIeCHH36QVq1ayaefflrrY1HevHnzlP05duyY6HQ6OXfunIg8GBGgWbNm8sMPP0heXp4888wzHP2pno0bN05mzpwpBQUFsmvXLmUUiO3bt8vly5dFRCQ9PV2eeOIJi1FDyktNTZWTJ09KSUmJ5OTkyEsvvSTDhw9Xlld3jiv68ssv5dy5c1JWVibXrl2TIUOGyIsvvqgs79Gjh8TGxkpRUZH89a9/5WglTuJR1NPi4mIpLCyU6OhoWbZsmRQWFkpJSYnVvGq6NrGeOqeGVk95PSVrHkU9rcs9tL2vp3UKKjIzMwWAeHp6ilqtVj5HjhwREZGVK1dK8+bNRa/Xy/z586WsrExZNzQ0VABYfJRCAOLh4WGR55YtW6osx4kTJyQqKko8PT2lX79+lUY7mjdvnkycOLHG/UlISBAAFttVq9XKcpPJJDNmzBCtViuBgYGyePFii30SsR5UiDy4wISFhYmXl5eMGjVKGTLWaDTKwIEDRafTiVarle7du8uOHTuU9Y4cOWJRhrlz54qvr6+o1Wrp0KGDfPDBB7U+FitWrLAYhaqgoEAmTpwoarVagoKCJDEx0SKvTZs2SWBgoGg0GpkyZYrV0aro4d28eVOGDx8uXl5e0r59ezl06JCIPDhPrVq1EpVKJcHBwTJ37lyLoDsiIkL5f0hKSpL27duLSqWSli1bSnR0tNy4cUNJW905vnz5sqjVauVC9re//U3atGkjKpVKAgICZPr06RZDJqenp0ufPn3E09NTunXrJj/99JM9Dw81EI+ini5ZsqTS94H5B5uK9VSk+msT66lzamj1lNdTsuZR1NPq7qEf9fXU5X8LRERERERE9FBsGlKWiIiIiIiIQQUREREREdmEQQUREREREdmEQQUREREREdmEQQUREREREdmEQQUREREREdmEQQUREREREdmEQQUREREREdmEQQUREREREdmEQQUREREREdmEQQUREREREdnk/wNPaP4b4Xr3fgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Get a file to read\n", "filepath_raw = filepaths_raw[12]\n", "filename_raw = os.path.basename(filepath_raw)\n", "filepath_processed = os.path.join(DIRPATH_PROCESSED, filename_raw)\n", "\n", "# Loop over cases and show\n", "for filepath in (filepath_raw, filepath_processed):\n", " stream = obspy.read(filepath)\n", " stream = stream.slice(\n", " starttime=obspy.UTCDateTime(\"2012-07-26T03:50:00\"), endtime=obspy.UTCDateTime(\"2012-07-26T03:52:00\")\n", " )\n", " stream.plot(figsize=(12, 6), method=\"full\", show=True)" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.10.5 ('py310')", "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.10.4" }, "orig_nbformat": 4, "vscode": { "interpreter": { "hash": "312801dd0aa5461639464a6df6a6c417c90a144a1d808d7cfd3dc2851df10201" } } }, "nbformat": 4, "nbformat_minor": 2 }