Skip to content

Instantly share code, notes, and snippets.

@snsinfu
Created December 18, 2020 09:02
Show Gist options
  • Select an option

  • Save snsinfu/c2d540d336ffd3b3800d1b9d14af7f00 to your computer and use it in GitHub Desktop.

Select an option

Save snsinfu/c2d540d336ffd3b3800d1b9d14af7f00 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"import pandas as pd\n",
"\n",
"from umap import UMAP"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"random = np.random.RandomState(0)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Load data"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"\n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
" \n",
"time123456789...2122232425262728293000.142917886938691...372109099310.259705986938691...372109099320.341235986938691...371109099430.407438986938692...371109099440.848406986939692...3711090984..................................................................99953240.32000021005559513...531124102015099963240.39000021005559514...531124102014099973240.41000021005559513...531134102014099983240.99000021005559513...531133103014099993242.17000021005559513...5311331040140
\n",
"10000 rows × 31 columns
\n",
"
"
],
"text/plain": [
" time 1 2 3 4 5 6 7 8 9 ... 21 22 23 24 25 26 \\\n",
"0 0.142917 8 8 6 9 3 8 6 9 1 ... 3 7 2 1 0 9 \n",
"1 0.259705 9 8 6 9 3 8 6 9 1 ... 3 7 2 1 0 9 \n",
"2 0.341235 9 8 6 9 3 8 6 9 1 ... 3 7 1 1 0 9 \n",
"3 0.407438 9 8 6 9 3 8 6 9 2 ... 3 7 1 1 0 9 \n",
"4 0.848406 9 8 6 9 3 9 6 9 2 ... 3 7 1 1 0 9 \n",
"... ... .. .. .. .. .. .. .. .. .. ... .. .. .. .. .. .. \n",
"9995 3240.320000 2 10 0 5 5 5 9 5 13 ... 5 3 1 12 4 10 \n",
"9996 3240.390000 2 10 0 5 5 5 9 5 14 ... 5 3 1 12 4 10 \n",
"9997 3240.410000 2 10 0 5 5 5 9 5 13 ... 5 3 1 13 4 10 \n",
"9998 3240.990000 2 10 0 5 5 5 9 5 13 ... 5 3 1 13 3 10 \n",
"9999 3242.170000 2 10 0 5 5 5 9 5 13 ... 5 3 1 13 3 10 \n",
"\n",
" 27 28 29 30 \n",
"0 0 9 9 3 \n",
"1 0 9 9 3 \n",
"2 0 9 9 4 \n",
"3 0 9 9 4 \n",
"4 0 9 8 4 \n",
"... .. .. .. .. \n",
"9995 2 0 15 0 \n",
"9996 2 0 14 0 \n",
"9997 2 0 14 0 \n",
"9998 3 0 14 0 \n",
"9999 4 0 14 0 \n",
"\n",
"[10000 rows x 31 columns]"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"data = pd.read_csv(\"_output.dat\", header=None, sep=\" \")\n",
"data.rename(columns={0: \"time\"}, inplace=True)\n",
"data"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### UMAP projection"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"umap = UMAP(random_state=random).fit(data.sample(1000).iloc[:, 1:].values)\n",
"umap_points = umap.transform(data.iloc[:, 1:].values)"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUQAAAEvCAYAAAA92bhfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAAA1C0lEQVR4nO3dd5wcdf348dd7ZnevpueSXHqBVEqA0DGABCQIX4qigoXiD0RFBRVF9CsqfhXFgiIWFMSvXxDpoYlUCaFfAiQhIZCQdql3uZTrW+b9+2PmWnLpNzu7d+/n47G5ndnZnfdebt/7+cyniapijDEGnKgDMMaYXGEJ0RhjApYQjTEmYAnRGGMClhCNMSZgCdEYYwKxqANob+DAgTp69OiowzDGdDNz586tVtWy3R2XUwlx9OjRVFRURB2GMaabEZGVe3KcVZmNMSZgCdEYYwKWEI0xJmAJ0RhjApYQjTEmYAnRGGMClhCNMSZgCdEYYwKWEI0xJmAJ0Zgc0ZzO8NW3a3l9czrqUHqsnBq6Z0xP82JVkqfWp5lb08i/N4MIvFKT4Y2T+0YdWo9kCdGYiMyuTnHqnDoyqnhxwAX1YFSRQ0NGKXYl6hB7HKsyGxORY/vHiAHqBjvEvz1cneIrC+v3+/UfXNdM36eq6PPURj765mbGzq7i8Fdq+PGy/X/t7spKiMZk2VNrk/zo7TpebwRtVyQpUPj8qAIO6Rvj0N7uzl9gD/34vToagyT77OYUAOvSaX61Ks33xpXs9+t3R5YQjQnZtqRy1et1PLY2Q52nIKACxP3H3Yxw5yFFnDeykMIuqiZvTnrMr1U8B+IJ/9ok+EsONwv8bnU9V46wpLg9qzIbE7KfzG/kvsq2ZAhQXgAn93Oo/Egfkuf158IxRV2WDAGSnvLhATFUab2JBIlR4OWtyS47V3diJURjQrRsW4Zb30v6iTAoflw4HP5yfP/Qzrm0PsP0l7ayIan+ZUkHUP8+ojgiPFaT5PGqJmYOLMARa7xpYSVEY0KS9pRTn6htbSxxPBgv8MOpfUI97wPrkqxPKgo4QZHHLyEqMdf/icCnFm/llsqGUGPJN5YQjQnJvE0ZalKA598mi7Cy2eHNTZlQz3twqf+xLnXAy4CbgURCiQftNC0FQhGYVGKVxPbst2FMSA7sJYgKsQwg8PqFfdjU5DGgMNxyyAub08SAOg/A4bC+Hu+kwPPbVPC8ICkqnNa/INRY8o2VEI0JweZm5YLnGv2GXYUPPtELIPRkCPDjCcWsmtGPa8YVcuogl3dSGVoCEWBQ3KGfAxWH9ws9lnxjJURjQlCxMcXrVR6CgAf9CrJX9og7wpBC4WeTSmjMKENfbUQFHAd+MKKErw3vlbVY8o2VEE1WeJ7yyztS/O8jPWPigpc3+B2hUfjMmOji+NS7Na2fcgG+NLQ0umDyQJckRBG5Q0Q2isjCdvt+ICJrROSt4HZGV5zL5Kd0Gh6fA399WDn54hR1DUoqrXz7xgxPzvaiDq/LzVre1nDyrWnZT0KeKo9WN/JSXRIJSocjEkLcsS42u9JVJcQ7gdM72f9rVZ0a3J7oonOZPJRICGdNDzYcOPPLaTZtVhqblVQq0tBCMaq3/9G65hCH4aXZr4g9s6WZi5ZuRlEcx799sqwo63Hkmy65hqiqs0VkdFe8lum+rr44ziOzU6j4Exp88jsZnr8thtPNSi01jWleXOMxugS+c3g0VdQT+xTgBkMERfx+iJ8ZbEP1difsr64rRWR+UKW2Ji3D6OHtNhw4+Yo0jc3dq8o8/aEGBNjYGF0MnioqiiuKg3JFWZyRBdaGujthJsQ/AOOAqcA64JedHSQil4tIhYhUVFVVhRiOyQXfvSzGzde4JIKJDRRl5lcyLF2t0QbWRRpTHhuCRDh9SHQl3yLXwZW28csPbemG1yVCEFpCVNUNqppRVQ/4M3DUTo67TVWnqeq0srKysMIxOeKAEcLUiQ5P/z6OoqgDGUe55IYUv7s3v1ugX12fZuTf6/wNhYMHxCOJY3Pa46PvrKP9EOUry62rzZ4ILSGKSHm7zXOBhTs71vRMJx3ZMoYMcKA5lb+lxAXVGc56osHv/xyMArni4GhGgbzXmGJuk+dP6gBcMrCYLw6xhLgnuuSigoj8AzgJGCgilcD1wEkiMhX/T2QF8IWuOJfpPq75dIwxQ1Lc8QQg8OgrSrwozRVnuSTi+dXQ8tu5jf6YZQE8uGqKQ7/CaN7D+43NraXDb5WX8plB1vdwT4lq7nwrT5s2TSsqKqIOw2TZ/GVJrvhlMGmqwMM/ijG4f36NGRj6p1oy0Fo6XH5pKcURJfVh89bgAT8a2ovLhvSOJIZcIyJzVXXa7o7Lr7860y0dMi5BUSLYcOCc76dZtjbDxTclmbMw3JlhusKaWg/1wFFwMnD7jKJIkuG/tjTyg9WbaWmzv9hKhnvNEqLJCc/enKC4AP8CiwMX/izN4rXKN25P82RFmjueSpPxcqc2094/FqUQFcjAABHOHBdN95YHaxr4U1W9P98hStyxj/fest+YyRnP/irhrzdC21T7ADfPynD7Uxm25uhicUeVOzgexFT4+KToPlITE4KIIqJML42mhTvfWUI0OeWV3yaI0bG6+YXTHfqWwBnXJzn9e80RRbZzt1QkEaC8BL59XDTD4+6vqedX1fW4DjgC/xw/OJI48p11XTc55/Gfxrl/doZRBcIBk4W0B+mMhypsrocllR4ThufOd/nHJsR4bW2KdfX+9cQD+u//EqJ7I63K1ytraBkBOS6WXy30ucQSoskZr7/h8cc/KWfNhEOHOfz0B8rnLlcOmOjw7x8XUL3VY+5S5cChufOBb04r35+douXq5tx16awnxCPeWQ2ABH1t/jNl+K4ON7tgCdHkjBdfVKqq4Y6/AygUwN//CUVFHn+/02VgH4ePHBF1lB1d9lgTjWlar3l+fFJil8d3tZvWbmKr+tN7eZ7y/Pjy3T/J7FTu1DtMj3fJRUD7IbcONCehfLvPeC71nV2yMegWFIT0twXZGzN85fIN3FpT35qMVx0yknGF1piyPywhmpwxf66DkwG2m/xm0KC2+w2Nyse/lOGuWdHPkLN+m8f6Wr//IQoxhZlZ6nLzYl0Dj9X7DUwOsOKgkVk5b3dnVWaTM8ZN8BdBctRfPhMXPFGenwNzXktz/10OBQnhQ0cJ48fAv55PM6Cfw1FTd/297nlKMg2Fif2/9tiQ9LhuVopN9R4Va4ASQP2k5ClkY2rHWZu3cfW6La3D8z6YYsmwq1hCNDlj2AihyIF0E5QPh41blFoAgWQazrrAoyAB55wBtXXKr+4AEY+n/nfXCfE3s9LcN8ejfy84+1jh86fGdzopredp62PJtJLYrsX2kr81s2CdP8xQRZGgoNoy7LCmPsPgkq6reDV5HvdvrePs3qX0ch0+98E65jT51XIRWDDRGlC6kiVEkzMcR7jwUnhlttJQD/c85LJmrXLFVRkyHqDQnIK7HgaCtouBezDt8PhhDhnPY0Mt/PEp5Q/PJom7Skog5sLQvtCUguMnwP1z4bQp8PZqWFcHCIwrg0tPcHlrtceCtcGLCgiCmxK8eHAZT+GafyU5cbRyzYld07jyRmMzP9xYw1821lDZeunUz74V4wZTYqNRupRN7mByjucpngexdqWzxiaPr13jsXo9pGP4X+UKRx4MP/n27r/XN9UqX7g1yQfVgBOMhnEBAa99ITC47zmwXf9wPNHWkqACuARzOirJBEG9GeIC7319/6frT6py/opK3k2lyXh+txpV8DzhtNJi/jjS5g/dUza5g8lbjiMdkiFAUaHDsUcLOOAEpUWANxZAKr37L/UBvYT7ry2gtYdgy8srBEN/O+qkzUYUSl2YfXWCR68o4LBhEBfB8RwIGoNE/RUGN9TufaOPqnZoQb901RreTaURAdcBL3jfH0weZckwJFZlNnlj+vEOjuNR2hv69FJ+czuMGAHxvRiZcdPFMX79WJpN9UKfYrjmbJfaRmXjNuW0Q2IoMKi3XxJbVaOkM0ptk/LXlzOcNCHGBUf5H5mBveDui4tIZZRZ89N8+4kkGpQqBei7FyP4NmcyfHXdOt5MNuNpUPxU8adXbCmxen7pZWZp8Z6/sNlrVmU2Zj9t2Kbc/mqK2+emUNdPXIWuMO+qIhLu7pP1ssZGzlm3BhHIqKB+iw2q/mp5o0X4ybAhHF5cGP6b6ab2tMpsJURj9tPg3sJ1pyX41ow4R/yikQYPmoHFGz3ujW9kYzpFVTrNFlX+OGgYByYKcEWoyWR4ub6W79ZU47p+8lMPjhaXSoEJ8QTXDSljaNw+ptliv2ljukjMEd7+VjGLN2RIezC3ZCtPbqujtRKmcOEGf9yxKmhwIbOlWqwKZY7D70aNsNbjiFhCNKaLTRrsN900NfmLTLUmvNZ/OuoD9Mbh/D79+Fx/W748SpYQjQnJkYUl/Kd8LFdVr+WtVBMitJYWFTg6XsCfho6INEbTkSVEY0LUy3X5/aBhHL9mmb86ggMvDR1LXBwc2fPWcZMddqHCmJAViMMw1+9JriiLU82WDHOUJURjsqBRNRhwp1xWXcnnN66KOiTTCUuIxmTBU0PHIqLEXHBEmZ9q5L66zVGHZbZjCdGYLCl32pYWEKAynb3JZM2esYRoTJZMihe2tjKLwMeLekUbkNmBJURjsuRnA4fz+dJ+iPgJ8fwau46YaywhGpNFl/ceRPuu1/+p3xZZLGZHlhCNybKf9h+BP3BPuW7b2t0eb7LHEqIxWTbIjXOQG8xcI8px6xeTymSiDcoAlhCNybqhsTjf6jeE9n2zT6x6j6QlxchZQjQmAuPjhfx1wChKg20ROLn6vUhjMpYQjYnMhHgRTw2ZxDGJtlmwT9mwiC1eOsKoejZLiMZE7Ff9R5HAH9bXDJy96T02W1KMhCVEY3LAz/uMQMSfDUcVvrt1NffWV/OPhmqa1cNrt9RHLi370d10yfRfInIHcCawUVUPCvb1B/4JjAZWAJ9QVRu8aXq06xuW0U9ijHAKOTNeRkEwM/aRhaV8orkv9yW3gMCidCML0w24jvKX+g1kAFdgWryEpZkm7uk3nrhYeaarddVv9E7g9O32XQs8q6oHAs8G28b0aGmUinQtdyTX8fXGjo0oX+kzlCNjxYAyxHEY5QqCvx60iOIBryXrSakS237RaNMluiQhqupsoGa73WcDfwvu/w04pyvOZUw++5/iA/hiwTAcYIXXxMNNG1mRbuTF5GYu2LaQBWwl4WaolhTryCCirUP9wP9Zq15ny0abLhDmjNmDVXUdgKquE5FBIZ7LmLyw1UvzSKq6NaHdllwLSX8tZsWfBccVyLS7TCi0Pday+6d1lXyvly0/0NUivwghIpeLSIWIVFRVVUUdjjGhWuzVs9CrB9rWVwFoqQE7DkGJ0C8ZDhAXR2CSFPPt0mHc1fcA+jguk2O2YH0YwiwhbhCR8qB0WA5s7OwgVb0NuA38hepDjMeYyB0T60Mh0IQiEixKD62L0rdfiOr/eo9nsFuAqn9si4f7T4wg8p4hzBLiI8BFwf2LgFkhnsuYnHZNwxL+2lTJDfXLSKI4gIuSQBjjFCIKRQg/LxzHUEnw+YJBDHZbljG1BpRs6apuN/8ATgIGikglcD1wI3CviHweWAWc3xXnMiYfDXEKqNM0b+jW1lLgV50RzCgp2+HYvxZYCTAqXZIQVfWCnTx0Sle8vjH57gsFw7m2fjEStI6IwIhCuw6Ya2xdZmOy4JcNH7CGJE5Q+z3XGcwEtyTaoMwOLCEakwVb1F9QquVy4LRE7wijMTsTebcbY3qCG0smtZYOAR5Jd9rpwkTMEqIxWVAgDg+UHMFNiQkA9CMecUSmM1ZlNiaLxsVL+bN7MP3EEmIusoRoTJYNcBJRh2B2wqrMxhgTsIRojDEBS4jGGBOwhGiMMQFLiMYYE7CEaIwxAUuIxhgTsIRojDEBS4jGGBOwhGiMMQFLiMYYE7CEaIwxAUuIxhgTsIRojDEBS4jGGBOwhGiMMQFLiMYYE7CEaIwxAUuIxhgTsIRojDEBS4jGGBOwhGiMMQFLiMYYE7CEaIwxAVuo3mTHphWQSUHNSqheDuOmQ/mEqKMypgNLiCZ86xfDrOv8+6qgwILH4cQvQbIRpp4daXjGtLCEaMJXMmDHfRkPnvsdODEYdjCIQP+R4LjZj8+YgCVEE75NK3f+WCYF93zNv3/iFXDIR7MTkzGdsEYVE76hU9rui/h/dfHCtn0a3A44LsuBmbxQtQaWvOnfb2rwL7uEJPQSooisAGqBDJBW1Wlhn9PkmEQxFBRDqtFPiADHXgIv/x80bms77qlfw8xvAwoFpZGEanLQXb+Euc+1bU85Bq79YyinylaV+WRVrc7SuUwuOu8XsOQ5aNwC40+BIROh7wh44DtBQ4vCynnwpwvAjcNlf/cTqenZ1q+Eec913FfSO7TTWZXZZEfvcqhdB6vegAGjwcvAgJEw7CCIFYOnkPYgnYbSgRAvijpikwsa6iCx3d/CRdeGdrpsJEQFnhKRuSJy+fYPisjlIlIhIhVVVVVZCMdE5qD/guMug/uvhJf+CPVbYO1iaK7z/0rA/+nRVrU2+2VBcgvn1Vfwlfr5UYeyb375JWhu7Ljv7ZdCO102EuLxqno4MBP4sohMb/+gqt6mqtNUdVpZWVkWwjGRGTQexh4Pkz8Ko4/xS4hfuAfYLvl5mUjC647+k6kBYA1JftawiPfTdfwxuYI/JFdEG9iemhg0OQjgCsQdSDfu8in7I/SEqKprg58bgYeAo8I+p8lxh54LI47w7ycKofcQ/37G82+bKv0O22a/faVoLI6A68Dr2sB3kot5MVO9/VdQbqpe418/dPBvrgMXfw+OPC20U4aaEEWkRER6tdwHTgMWhnlOk4dOvxoKe223M7yuFT3NSfQBVWKOR8xRUihDpSDqsHbN8+B754HjQMwFN7hNPw9K+4Z22rBLiIOBOSLyNvA68LiqPhnyOU2+GX4QfP4O0HZ/js+E062iJ7qoaAx9RXAE/4bySHpN1GG1eeNxuP40+Pt34X/OhZ98DL7zYdB2l04cIJGAOQ+HGkqo3W5U9QPg0DDPYfLc0zfC5pWwptJvYXbxSwWTT4o6sm6jl8T4U+ER/LTpHRZSj+NAnXpRh+VrrIOHfu7fX/Ri0Kjm+Q1rpQV+z4PGlH8NMd0MD98CA8ph0tGhhGPdbky0Bo2HZAMkABQyCqkMVC6KOrJuxRHhu0UHEQuGBTm5ckkilmi77+FfQ/aAkj7+Psfxrx06jt/zoLkebvsGbFobSjiWEE20Dj4bSgZCIt5x/wt/iyaeFqkmePkeqNscbRxdqFEziASFcFEaNB11SPCP73fcbsnT9Vv9BHjU6Ts+nk7DwjmhhGMJ0URM/T/88kkwfGzHh9YvjSYkgA3L4PnbYXlF9s/97hx44ma/tX3Jy132sgU4lEkcEb8G+s/kCqq0qctef58cfjr0K9/5+OT5T0Pv/m3bGc8fyfShj4cSjiVEEy0nBuf8AqacCVIP4w5re2z53OjiGj4FrrgDDpqR/XNvWApLX4eX/wEP/gjSyS55WUeEy9wD0OA76Hmq+Wkq4ksTB50E19wDkz7kbzvtOgQ5Ag1p2FbT8Tm9B/hV6BBYQjS5Id0M6SZo+ACOPAUGlsLog6KNacCIaEbMnHgxfPVuOOVyuOR3Ha+z7an5T8O2HacPmOj2CfogKohSSoyNXsSlRICzr4JhYyHhQixoDh91ECQ7qdafd3VoYVhCNLmhPigFxApgzVtwwR+gfMoun9LtFfeBIQfs/fM2LINHfw73fb/Th38oUxD8avNK6vldZsl+hdkl+pTBplV+A0pBHIriMOXEHY9zBAaNCC0MmyDW5Iap58Hg8TDvn7BhMTTVto1gMXtn0FiYdg5Mmt7pw2MSvZAkIIoAjaSyGV3n/vaNHfeNmgjfvA3KRsCzd8PQA2DqSVAQ3sQfVkI0uaP8IH9KMID5D0YbSz4TgY98GUYevNNDrpfJrcP3quiaa5T7ZfU7O+6b+xiMm+pfMzz3K3D0zFCTIVhCNLlm6X/8n5VvwjuPRRpKd1bt+NfmWqrOjVF3wfnoVTB8Ipx7LXzjXv9yQf2WrIeR11VmL+231rvbdWFr2gaF4c0hCUBjvbKtGgaP6njR/ekHPA46UigfKbz+oscLz8I5n4IDJ9p3zx7JBKUVLw2v3Qnb1sGxl0UaUnd0hNMfCQarlCIUScSp4LDT/VuLS34FBSVZDyOvE+ItU6B6md9zY9Ah0H8CSAzefgCGTvM/RyOmwqCJu24sbNjizzhV2snicNurrVHihXD/zfDiQ1A0yKOmGiQOUgBNCjFVYnGlIQMEney/0fn1bbO9geP8/9QWNRH2RezGHBEGU8gGmhhPn6jD2dGgMZGcNq8TYuMWwPVHe1UuCEZ7BYlvzVy4781gdvoCGHQgXHq3n/S2/+L5wwV+H+CjL1LeeAIKi+H4T0GsEA49GR78LZT0gjmzINnkl0gHj4HmtNJY48egHmSawEtAMgZJBRx/Crczzs3qryV/qcK446HmA/++4/iL2qvCmnkweLLNpN2FNtMMwLHYPKQtRENcwWpvTZs2TSsq9nxkwJPfhJd+C56AF8MfkxTwABzw3OB+y2PilyiduJ/EUs3+MepCqlBBwHP845RgAhYHf0NAxd/nOZBxFW33uh6QTABxWhNzUQFc8yPh4MPzYga6aD3xPSjuB4dfAA9d1ba/tD80b4XDPwcTZ0YWXnfzu9Ri3qSGIVLAl5zJDHO67xo2IjJ3Txa4y+sLW+89Cnggin89pLOJlj12nJA57U+6m2mmNdkpQCo4Hv856gbPbf+d0W6qe1FBtjunZNq9hsLwkVgy3FPDDoMhB0PvoXDRvf5iVHHXT4bxYhhqEyd1pYlOHxRYr8383Hs76nByQl5XmSeeDZt+AcOnwck/hlUVsOAhKJ8KVe/BiGkw8Qx49Lsw6hgo6g8bl0FpGWxYAstf9V+nJV2JEyTBIOmJ51eNNy73r/fO+BwcMLVjcmtuhHfnKSqwaik8eT989mtQ8TpUvAQHTMzWb6MbOPRjHbcP+xQs/49/P9UAK16CQ87Peljd0bOZdczyKgH/zz1npgOLWF5XmVWh+l0om7Rv51s5F955Cl69B7bWQEoVTfgt1FIIE4+D//cLWL0YysdCoshKeln3ws2w5jX/viPwqbsjDae7eM2r5uHMSjbQ6Dc4qnBTbBr9c30m7X20p1XmvE6IXWX5XLj5HDjz23DSZRDvnn8T+clLw72fa9s+6xYo2YPuAGa3qjINXOfNbb0sFBeHW9zjcLvhioc94hpiVxl5CHzyRjju05YMc44TgzN+2bb9WHgD+3uaf+pyfwEq8edHzJDhu5k3og4rUpYQ8bvRHPdpKOkXdSSmU73LYeRxwUYGtoUzW3JPUyQugiJ4uOLhilLUw1NCz373Jn8cd2XbNPJPfivqaLqFvsQRFBcNeqUpKkouXUbLNkuIJn8MP6bt/rM37HyWZbNbdV6SZ1lNTNTvYit+MthAIxu0566JbQnR5I/jroR4MEh90xKYf1+08eSxUifBtGCEiiseDhkclIEU0q+btjTvCUuIJr/812/a7r/3SHRxdAObnc30c9L0kmZ6SZJTnUHc4B5Fgbi7f3I3ZQnR5Bc3DjNuaGkahQc+G3VEeStJGkgRjFTlBIlmQoVcYgnR5J/+Y2BEu9mgH7o0uljyUAaP51lGnTSCECxN6nCjN5f3tPssu7ovLCGa/HTUZX6rs+MAGX+GHLNH/o95PMP7rXOWAByswxlOL/qwDwtadSOWEE3+mtnueuILP4RkQ3Sx5LgqbeCf3ns0aIqV+At6tY5eVniNNSyjhi25sJxAhCwhmvxV2BcGtpsB519fiiyUXFdJHa+zght5nuZ2Uzq1lBALJMM4KWQkvaIKMSdYQjT57UNfh4LebQucP2rXE7e3XuuZyyoKHY80Hmcxie/yYa5nBpdyFAlxQZRzZXL0SwlEzBKiyX9n3ALFwZKl6sGSWdHGk0Nmee/xK32FD6QaUeEIhnIMoygmQYIYYxnARUzDAe7k1ajDjVzP/jow3cepN8Ljl4Gm0fcfhvefQQ44GSaeF3VkWfWIvsdmmjib8dzAbGKipNShQON8VY6mjNIdntMbf1mGDDYnopUQTfdxxm0tKz0g1KFLn4g6oqybw0reYT0/4QUSaNDHUPi2Tu80GQL8lVf95TKAjWzLZrg5xxKi6T5EYPhH2q34kMFbPz/CgLIrjccxDMcJJmtQ/FUxfuXMoNSJ7/R5zbStybzWEqIx3Ycc+gnoO8FPio6DzvsN6SWPRh1WVtytb/I6q1rbl9IIP+e0XT7nXTbQSArwS9YbqOUB3uqQJHsSS4imexEHOeFa9NRb/fXFHEGWzyJdcUvUkYXKQ1lMNY4ES+8Cl3E4zm5mvx4cVKMFKKOEt1jNItbRFCTJnib0hCgip4vIEhFZKiLXhn0+YwCcgmKcmX9pXXRRqt8mveTBqMMKzULWI0EynMlYbpSPMEEG7vZ5cdwgCSgbqaORFBm81lJjTxNqQhQRF7gVmAlMBi4QkclhntOYFiJC7PS/tK0Ku/JJkk9/NdKYwjKG/kyijKvleKbLgXv8vA3UQmvjS9ttA1uppSmcYHNY2CXEo4ClqvqBqiaBe4CzQz6nMR3ETv8LGcfxP+zaRPPTl5Nc+iiaqo86tC7TiwI+y+EM2klL8s6Moj+nMIFR9KeAGA4eDh6P8Ba/4RlWsimkiHNT2AlxGLC63XZlsK+ViFwuIhUiUlFVVRVyOKanih33A7xgBW4BdMVjNM3+Jtq8NdrAIhbD5XgO4CKO5RvMCEqIflIUPB5lHs+wMOowsybshNjZFd0O876r6m2qOk1Vp5WVlYUcjumpnNKhJE77E3rkdeCW+PvUo+nla2l86Vo00xxxhNGL43IKk4JuOwp41NLAwg5lmu4t7IRYCYxotz0csCXTTGQS/UaT+PCvSUz/JTJihr8zuZXG2VfR8PwXaVr1NJm6ytbjVZVkTc8pIR3LAa33Wxae+jhHRhVO1oW6UL2IxID3gFOANcAbwIWq+k5nx0e1UL3pubx0E80vXt06wsUD1BEKD7uGzKYFJLcuIVO/CkdiJMZ9goKBhyNu915zZDGVzKKi9fcxijIu5Piow9ovObFQvaqmgSuBfwOLgXt3lgyNiYITKyR+3E0ArUkRVZrevInk6n9D7UpEFTRN09K72Tr3+zRUPk2mcUN0QVevg99+EzaEU5WdxHDOx1/hUICqHjR6JfR+iKr6hKqOV9Vxqvo/YZ/PmL0VKyjF6TepdTwvtCVGAFHa1irONJJc/ShbF96c/UBbbKmGt+dATXhJeRxDWrvgfJJjdnd4t2Gz3RgDFE79KpnN79I0/7dAMJu0SGurYPsLSwqQaUC9FLKLMcKhOeBguD38qbp6U8xWGigkgvcYERu6Z0zA6Tue2MiZSN/x/lot4s8onQZSiRhpRzokxvS2ZRFFmh0nMBGAV3k/4kiyxxKiMQERh4IxZ1F88JcpPPhreOJ35vZiLiqCxlw8ETwR3H6TifWZEHXIoSoKFpxaTOVujuw+rMpszHbEiUNBKQ0FMWIZr0OxIV52FAX9DiHRbwqym4kT8t1Q+gOQJMMS1jCh45iKbskSojGdiBeV44hLJiZ4qkEDg0PpyLNx4ns3PC5fFZFAgvf+jMzvEQnRqszG7ETGcfx+ia6Dug5O0cBwk+Ebz8Jz94f3+vvA/w141GkTKTJRhxM6S4jG7MTgg7+PxhzUEdQtpmzKt8I94eyH4d93h3uOveTiVyMdlFt0Fk/pvKhDCpVVmY3ZCTdeQu+hZ7J13eMQiyEScvnhK78AL7dKYS2t6jEUD3iH5SzS5ZzGkUyWkVGGFgpLiMbsQq9BJ1LYexJONobrJXJrSOBSXctwyqjEn4XKCdJjjFhrC3R3YwnRmN2IFw6KOoRINJOikWYSuGRI+9OmASdxGGNkSNThhcKuIRpjOjVFRnGRnMqRTOgwj9/TVNCo3XO6NEuIxphdOlomUYRfnRfgFA6nSHKret9VLCEaY3brYMYCcAxTOETGRBxNeCwhGmN2q55GAJpJRhxJuCwhGhOGjZUw59Goo+gym6gFoIot0QYSMkuIxoThmXvhtuuhdnPUkXSJwxgH+C3P3Zl1uzEmDOddASecCb36RR1Jl0gGw/Y2071XKbQSojFhKCyGkeOjjqLLpEkHP5WGbtrlBiwhGmP2QEJjiHqgHlV0j8sAnbGEaIzZred4A/D7IUp4C3VGzhKiMWaXqnQzwxncuujUM7wedUihsYRojNmlxbqcNWxs3U4Qp8brnkuTWkI0xuzScXIovSmBYKHWLWxjOWuiDisUlhCNMbsUE5da6gBaq81DGBBpTGGxfojGdEK9FOnkVuKFA3d7bDK9jXfX/AaAmDuA8UM/T8zpPpMfeOrhqYdDMGGsCLVSH3FU4bCEaEwnaiofpHHLPMQpYtCBV+PGCgE6nSjWbZf80plNLFp9E6PKzqdPcX4vU5pRj2W6mg1a1Tr9lwP00mImyOgIIwuPJURjOlHc73Aat8xDvUbWL/mJXzJyHIZN/B5urATPS9GUXE8yXcem2nmgCq3Lkiorq+5lWP+P0r/0sLxcrtRTj5W6hmf1FYAO8yEqmpfvaU9YQjSmE0W9DgSKgMa2ZOB5LK/8C4nCIRQkBlO1dXbrmiMCiFOEp42t+9bUPM7qmico73cKg3sfm+V3sH+e9uawno20raoCLi7nygxKpSS6wEJmjSrG7ES81yja90FWwK1fT7rmbaq3vOiXCluOjZUxZcTXGTv4ktZ9HqAiVG17I2sxdwVVZTNbadYkMVWc4NabIsqc/t12cliwhGhMpxpqF5Ose3+H/R4gniJeOhi14S/kPnHEF9nauJwlVf+LIigCwSp94uTfgkx96AX4JV8HKCDGMc5hkcaUDVZlNqYT4iRwS4aTrF+Loyk8R1Dxb02J4GMTXDd03d5sbljOss33AkpGpLV7CkAqUxfRu9g3IoIGib5FnDgjKCepSRZkFlKohaRIM8wtp8wpiyzWrmYJ0ZhOFJWMo2jsl1i74i80N3zQmhxcIKXamjQyIjRQR/3me/wDRBCBuAzAdYroXTiGIX2OjOpt7LNpHMwa1gF+Ym+kkZe9CuIqvKtLaLm2OD/zNjM5vdskRasyG7MLbtFwPOhwa0mGadelubAAXLf1+N4FBzKk13SadRN9SkYxvN+HiDmFEUW/7wbFBvIZ5zyKKECCESrrdCPN6k8Q21KVBqUiUxFdoF3MSojG7ILQvjuNzwE8x6Wp0G9c6F94EAMLx9O36AAkuG5YUlBOSbw82+F2qSKnkL5eH9bThIN/XVFa06BPgFpqyWgGV9ydvVTesBKiMbtQ2msynvglwxZpXHoNmI5DguLYIEb3/TD9ise3JkOAPoXjiLnF2Q+4iw2XoYCfANfoaj5gKQB96N16GaGJJu5K38Wy1DIAUpqiWZsiiHb/hVZCFJEfAJcBVcGu61T1ibDOZ0wYiotHMWTohVRXPYG6RWSa1tNQ4FLb+DoT+19AaXwIbh62Iu+pIQz0E18wdA/8/ohnx/6LuZm5LNbF/k5VXmYOrzS/QJwYLg7HutMZFhsZUeT7JuwS4q9VdWpws2Ro8lKv3lMoKZ/JhkQdW0qLycTjjO41kz4FI7t1MgSo0Lda77dUk4cwBEccjowdydFyNL3pTYw44FehM6RJkWR25hmWp5dmPeb9YdcQjdmFZKaWRTX30ZiuwhEh47oI0Ksgv68P7qkj3Km8lVlAkhQeKWqpZazTtlD9hNgEJuCP2X40OYstbMIlSJ6qVKU3kNBCxFGGuiMieQ97Q1TDmQ88qDJfDGwDKoBvqOoOizGIyOXA5QAjR448YuXKlaHEY8y+qFj/B1LUQrt+ecXuIA4ZdMkun9cT1Ws9K5MreJPXWjutt86QA5RQygmxGfSLZX/qMBGZq6rTdnvc/iREEXkGGNLJQ98FXgWq8X8fNwDlqnrprl5v2rRpWlHRfZrwTX6rT23k7U13An5V0O+s7HBs+TWRxpXLVJVV3gqWpt9jo1YCdOikDlBIMX2kLyfETyXmxEhqM/VaRz8nvES5pwlxv6rMqjpjD4P5M/DY/pzLmGwrjPULupl4fk88EcaUzow6rJwmIoxyxzDKHUODV8+yzFJGOCP4ILWEbWyhig0008QGXctLqWfYpltwcKljG2PlQKozKyl3x3Fo4lg8MqzJrODt1Kscm5jBQLezslfXCrOVuVxV1wWb5wILwzqXMWFwJY4jRaS0Hge/lLOqYTblpQdFHVpeKHZKONg5FIDD3bbZftJeig26lmZt4v3MIrZoDQqs9pYj2sQH6QUsTy9AAQ9hoDuUhGSnc3uYjSo/F5Gp+FXmFcAXQjyXMaEYWHQQ6xpfxwsuLXnaEHFE+S/mxBnGKADGxiawJVPDOq2kSIt5LfU8cfXaDZVUTiz4aNbmXwwtIarqZ8N6bWOyRcRDg8kaFBhT+uGoQ+p2+rr9EU+oSM3pcImi5Xe+MrmU0QUHZiUWG6lizC4MLT7Mn+XGccBx6Bsfgaq3+yeavZKkmVrdRgkFQUpsGx30Zvo5PC+TlTisH6Ixu1DVvByCyRwA5tb+gz7xkUztfU60gXUzZc4Qzin8NE1eA28lX+Kg+NEUOcXMbnyUNM0dhkWGyRKiMbuwPrm4daU5TwQhg2AlxLAUOsUcU3hq6/bJxedm9fxWZTZmF47o/cnWiWFbrmmNLjwh6rBMSCwhGrMLTV5t60iLlp/1WrWzw3Pa+vRKnmu4jybPWsp3xhKiMbvwQUOFv1hUsC0I5YlJUYa0zxwcXIl12yVEu4JdQzRmF9ZmFiHi4PnLRtEncQDVmfUMdIdk7UJ/VxkUG8GgWO5PsBAlS4jG7IaIv4qehz9J6prmNQxyRnBc0RlRh2a6WH59xRmTZQpk8PvEpYOGFVTZ7G3kveRbUYZmQmAJ0UTP8+DRO2Dlu1FH0sGmZCUZ18FzxF+G1BHimqE//ejrDKRJ66MO0XQxqzKbaP3xh/D0PyHmwMcbYNTEqCNqtbDpeb9EGDRCOIBLEUcUn06J0yfK0ExILCGa6Jw7GTJBJ2cROO+KaOPZzuHFH+XFhrtb+x8CTCw41pJhN2YJ0UTjv7YrCd63cIflPqNWq5vxoHVsLeLQpI3RBmVCZQnRZFf9NrjgqLZtAR58J+eSIcDA2DB6u0MY6AyjX2wwq9KLGJWYHHVYJkSWEE34arfA7TfC4nmwblXHx2blVkNKewkpZHrJx1u3y+NjI4zGZIMlRBOe9avhP4+CG4PnHt7x8QfmZz0kY3bFEqIJT8UL8I9b4NcPQL+B8PwjMOEQOPOz0L8s6uiM2UFoy5DuC1t1r5vxPNi0HsqGRh2J6eH2dNU965htwuM4lgxNXrGEaIwxAUuIxhgTsIRojDEBS4jGGBOwhNid/ese+MKZsGFN1JEYkxcsIeYyz4MfXg9lfeGwQ+Dxx2DVyl0/RxUuPgU+PBZuug7eXwRfPAsa6rISsjH5zBJiLnppDgzpC0MHwC9+CnX1sGgRnHcOHDgO7r1n58+tXA6rlnfc11QHyxaFGbEx3YIlxFzUvz80NEJDgz/pgbdd5/m77oK33+r8uZec5k+plWm3dnAsDv++L7RwjekuLCHmGs+DU7Zb93f7/6Wn/g2fvnDH527bAqn0jvszaZj/ClSt66oojemWLCHmGseBgoKO++Iu/Pf1kMpAMgNNKXhnMfz1rx2Pq6/r+JxTz4LCmD/F1pZqmPOv0MM3Jp9ZQsxFP/0FuI5/iwe35Uvapm1uceONHbfLh0Ms5lezn1kGl34LhoyAUz8Ov34Qzvps1t6CMfnIZrvJRaNHQ/kg2Fjlb4vA/bN2PG7mzI7bNdVQVg5Fpf724GFwx7OhhmpMd2IJMRe9swA2bQInmEVaxG9p3t7QcnjmEXjhCZj3CtRt9ff/4NbsxWpMN2IJMRdNmOT/bD+t/vP/ga99DWo3wfJKOOxAePxP8EjGf9xtd+zUo7MWqjHdiSXEXHRp0IIstCXFe+6EN9+Ez86AMoVEM2TaXVTMqH9F+Md/9idjNcbsNWtUyUUzPgJHHwsfO9+vNhcXwdXf9h/70e/hExf7fQuPPL7j84aPg+mnZT1cY7qL/UqIInK+iLwjIp6ITNvuse+IyFIRWSIiH9m/MHuYl1+AcQfArX+Fj30SMilYF4xHHjMerr4Bvv9buPL78Lmv+PszHix/3+/HaIzZJ/tbQlwInAfMbr9TRCYDnwKmAKcDvxcRdz/P1XPMngc/+41//xe3wpNzYPqHOx5z4kwYNxGaGjruf/6J7MRoTDe0XwlRVRer6pJOHjobuEdVm1V1ObAUOKqT40xneveBRMK/n0jAxCk7P/acz0D/QW3b5cPDjc2Ybiysa4jDgNXttiuDfTsQkctFpEJEKqqqqkIKpxsbORZSKX+8c0ahemPUERmTt3abEEXkGRFZ2Mnt7F09rZN9nS7vp6q3qeo0VZ1WVmZLU+6TL1/X9tv95qXwtc/Cqg8iDcmYfLTbbjeqOmMfXrcSGNFueziwdh9ex+zO334HtwZD+FT9bjovPQvDx8A1P4o2NmPyTFhV5keAT4lIgYiMAQ4EXg/pXD3bX272J31IZfyWZg/wHPi/2+CH34BUMuoIjckb+9vt5lwRqQSOBR4XkX8DqOo7wL3AIuBJ4MuqmtnfYE0nGhr86rLiJ8N0xp/uK+PBP/8X7vx9xAEakz9EtdNLe5GYNm2aVlRURB1G/nhvEVwww68qtyRF8LdT7b5/ho2EJ16FeDyKKI2JnIjMVdVpuzvORqrks3jcv2ZYNgSGDIVjPwTnXwS33gVjD/SPyai/Dsu7C/3tHPoCNCbX2FjmfDbmQHhgNvTuC/0GtO2f/QysXBaMgw4S4Mdm+IlTYnDCifDN/4aB1qpvTHtWQsx3o8Z1TIYARx4HMz7qj4P+zg3+vozCmrVQuQru+Tv8vwuyH6sxOc4SYndUVAzX/RT++yb47OXwyrsQ364y8NbcaGIzJodZQuyuygbDJz4Hrgv9B8K7G+GC7ZYQ+OQZ0cRmTI6yhNiTHHI4jBzpjyNygIpXo47ImJxijSo9yQUX+zeAl/7jr/1sjGllCbGnOv6kqCMwJudYldkYYwKWEI0xJmAJ0RhjApYQjTEmYAnRGGMClhCNMSZgCdEYYwKWEI0xJmAJ0RhjApYQjTEmkFNLCIhIFbAywhAGAtURnr8r5Pt7yPf4wd5DLtg+/lGqutsZkXMqIUZNRCr2ZN2FXJbv7yHf4wd7D7lgX+O3KrMxxgQsIRpjTMASYke3RR1AF8j395Dv8YO9h1ywT/HbNURjjAlYCdEYYwKWELcjIj8QkTUi8lZwy4uVmETkdBFZIiJLReTaqOPZFyKyQkQWBL/3iqjj2RMicoeIbBSRhe329ReRp0Xk/eBnvyhj3J2dvIe8+RyIyAgReV5EFovIOyLytWD/Xv8/WELs3K9VdWpweyLqYHZHRFzgVmAmMBm4QEQmRxvVPjs5+L3nS5ePO4HTt9t3LfCsqh4IPBts57I72fE9QP58DtLAN1R1EnAM8OXg73+v/x8sIXYPRwFLVfUDVU0C9wBnRxxTj6Cqs4Ga7XafDfwtuP834JxsxrS3dvIe8oaqrlPVecH9WmAxMIx9+H+whNi5K0VkflCVyOnqTmAYsLrddmWwL98o8JSIzBWRy6MOZj8MVtV14H9YgUERx7Ov8u1zgIiMBg4DXmMf/h96ZEIUkWdEZGEnt7OBPwDjgKnAOuCXUca6h6STffnYfeB4VT0cv+r/ZRGZHnVAPVjefQ5EpBR4ALhKVbfty2v0yGVIVXXGnhwnIn8GHgs5nK5QCYxotz0cWBtRLPtMVdcGPzeKyEP4lwJmRxvVPtkgIuWquk5EyoGNUQe0t1R1Q8v9fPgciEgcPxnepaoPBrv3+v+hR5YQdyX4xbU4F1i4s2NzyBvAgSIyRkQSwKeARyKOaa+ISImI9Gq5D5xGfvzuO/MIcFFw/yJgVoSx7JN8+hyIiAC3A4tV9VftHtrr/wfrmL0dEfk7fjVBgRXAF1quQ+SyoFvEzYAL3KGq/xNtRHtHRMYCDwWbMeDufHgPIvIP4CT82VU2ANcDDwP3AiOBVcD5qpqzjRY7eQ8nkSefAxE5AXgRWAB4we7r8K8j7tX/gyVEY4wJWJXZGGMClhCNMSZgCdEYYwKWEI0xJmAJ0RhjApYQjTEmYAnRGGMClhCNMSbw/wF7+/aFhROxogAAAABJRU5ErkJggg==\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(figsize=(5, 5))\n",
"ax.scatter(umap_points[:, 0], umap_points[:, 1], c=data[\"time\"].values, lw=0, s=3, cmap=\"rainbow\")\n",
"pass"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### Trajectories of some species"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAl4AAADCCAYAAACYGsyXAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABdXklEQVR4nO2dd5wbxfXAvyOdrhf77HO3wTbGdAyYDqaDIRBDAqGHAAkkoYVA6Am9hVDCjxaa6SXEFAOhxfRiwKaDbVxwt6/Y1/tJ8/vjaW8lnaTT6U5X3/fz0Ue7s7uzb2dnZ9/OvHnPWGtRFEVRFEVRUo+npwVQFEVRFEUZKKjipSiKoiiK0k2o4qUoiqIoitJNqOKlKIqiKIrSTajipSiKoiiK0k2o4qUoiqIoitJN9AnFa/r06Rbovt+TT1qMaZtujGXrrbvuPC+9ZJk82XLHHZL3559LemmprBtjKSnpeL7vvCPHfvhh95ab/vSnP/3pT3/6s8QhZYqXMeZhY0yJMea7kLRCY8xbxpjFwf/BieRVVlaWKjE7jt/ftfltsQVsuaUs77yz/Hu97vZAoON5OjImc6yiKIqiKCkjKcXLGPNtArs9AkyPSLsEmGOtnQTMCa73Lbpa8QJoaura86nipSiKoii9krRYG4wxv4i1CRjRXsbW2veNMZtGJM8A9g0uPwq8C1zcXl7dzvLl8n/66XD11fDdd/Duu5K2ZEnn8n73Xdh3X1leuxaam+MrXokqT9bCzJmw666qeCmKoihKLyWm4gU8CzxJ9LHKzCTPN9xauw7AWrvOGDMs1o7GmDOAMwDGjRuX5OmSpKZG/h9+GBYuhI8/7pp8AwHYbz/5NwYqKiAzEw4+GO64w92voAAuugjuvDPxHq/SUlEUzz5b8nPOpyiKoihKryHeUOM3wD+stadG/oCKVAtmrb3fWjvVWju1qKgo1acLJyMDrrxSlququi5fY+TfUaaMgc03h7w8OO88dz+vF26+GUaO7LjyZIz2eCmKoihKLyWe4vUnIJbWcVSS5ys2xowECP6XJJlP95Ge3vV5trQktp/HozZeiqIoitKPiKl4WWs/sNaujLFtXpLnmw2cElw+BXgpyXy6nvp6UVSqq6GhITw9ksrK+D1hgYA7XFlfL/vX1EBjo6Rt2CA2WX6//McjnvIUCIiNWOh+NTWu/M42RVEURVF6Bal0J/E08Akw2Riz2hhzOnATcJAxZjFwUHC95/nxR8jOhp/9DPLzZZhv1CjZtmBB2/3HjYPCQnjzzej5XX21DB8uXQpjx8ovLw+ysmT7mDHwr3/BokVQVxdbroqK+Mb8p54K224ry3fdJf8zZ8JJJ8nyf/8b+1hFURRFUbqdlCle1trjrbUjrbU+a+0Ya+1D1toN1toDrLWTgv8bU3X+DuEoP8XF8r/zznDGGdIb9cknkvb997JurfRgHXEE1NZGz8/p7aqrk+WSEpgwIXyf6mpR8iZPji3X1KnhPr0i+ewzUd5Aes9uuEGUR4AZM0SZVBRFURSl19AnPNf3S7rD/qq9YUxFURRFUbqVdhUvY8x5xph8IzxkjPnCGHNwdwjXr0m1UmSMKl6KoiiK0stIpMfrNGttFXAwUAScSm+xzepqvvxS/ocOddMcFxCeiKLy++EPf4C33nLTfvlL2f+222T9+OPFwN2Ytobul14Kd98df9ZiczMceqjYg0Vj4UJXxueeExmdGZPxhigVRVEURekRElG8gpoHhwEzrbVfh6T1D8aPh/vug9mz4b334IEH3G077gj/+x9MmhR+zK23wj77wNy5btrzz7vLJ50EN94Ib78tfsFefVUUrW+/hQMOcPcbOTK2XPfcI/+//33713D77eJA9c475Vy77KI9XoqiKIrSy4jnud5hvjHmTWA8cKkxJg/oXw6iCgrgzDOjb/P5whUlh0mTxPlpLKZMEQN8h223dWcgjh3rpg8aFDuPeIb3IHJPmCA9dc65hg4VuX78Mf6xiqIoiqJ0O4koXqcDU4Bl1to6Y8wQZLhRSZbIYcvOEC/Oo/Z4KYqiKEqvIhENwAJbAecG13NIPlZj/8Ja8ZX19dewbFnix5kuGKn97jtRumI5STUGVq/u/HkURek7lJXBnDk9LYXS3wgE4IUX9J3SRSTS43UPMrS4P3ANUA3MAnZOoVx9g333hXffhQsuCB8WTE+H3XePfdyxx4qPr5aW9ocTY/GXv8DEiRJce+nStttzcmBl1MADiqL0Vx59FC68UHu7la5l6VL4xS/g8svhuut6Wpo+TyKK167W2h2NMV8CWGvLjTEpCGDYBznwQOlZuuGG8HQnNFAsDjpIfp0hLU3Oe8ABsW3Q1IGqogwsVOFSUoEzW17rV5eQyFBjszHGiww5Yowpor8Z1yuKovQHusKMQVEi6Q6H3wOIRBSvO4EXgGHGmOuBD4Eb4h+iKIqidDuqeCmpQBWvLqVdxcta+yRwEXAjsA440lr7XKoF6zPk54uvLsfn1g47pOY8F18s/9ZK4/rqq3LuWOTlwfvvw5Ah7Q99RuPss+U8xojNWno6PPxwcrIrfY/tt4fLLutpKfo/CxbIM/b1153LZ/58yeeCC7pGLqV/8M9/hrfjF12UXD6O4nXDDZLXFlt0nYwDkJiKlzEmP/hfCJQATwNPAcXBNAUkoDbAzJlSOefN69r8nSDY//mP/Dtj7I2N4sA1FjvuKLMeGxpiz3yMR0mJGO+DBPo+91woLe14Pkrf5Jtv4IMPelqK/k9FRfh/spSXw7RpMmlHURxKSuT/P/+RzgFnvaMEAvIx5nzEO+8lJSniGdc/BRwOzCdo3xXEBNcnpFCuvkd6uvtlkQpC/XUZI45d28Pn65w8jr+x9HQNQTQQ0WGrvoXPJ1EyFMXBeYazsqC6Ovl8AgF5B6TrvLquIKbiZa09PPg/vvvEUWKSTK+VoigDC1WWlWh09sPZ7+9ax98DnHbdSRhjjgLettZWBtcHAftaa19MrWh9jFQ1eE6+xcWwahX87W8dO762VmI45uS4aZmZEtsx2tdxdTUccgh88glsuWX4tltukS7mP/0JRo2SWJT/+Ic29n2ZDz6QYfKsLKkToV+0H3wAy5fLUPM//gFTp4pvuJtuSr4RXrQIbr5ZbEVGjOiSS+jTfPWV/J9/PjzxBGy1VeLH3n8/fPyx3IuZM6Gw0H0WTz01/LnMz4cTT4TPPoOzzpK0ujqxCWtshJ//HI48siuuSEkl9fXS/jY3w/TpYtpSWgrDhslzBaIkXXABXHUV/PCDpGUGfZ4/8YT77Hq9Uv823VRsgocNk2c7Gp99JucO5bTT5L+6WurSpElQWSkh8bbbTpyLAxxzDBx6aNs8//xnqXuR7c4AwNh2/HIYY76y1k6JSPvSWpsiK/K2TJ061c7ratupruSVV2CvveLHXUwWvx+23lqUoDPPhF/9Ch57LPFG8pVX2o7r//nP8kCOGtV2/wUL3Mb/++/lRTt9ujjQe+89cdB4wgky3r/77tIApCXiDk7plVx5pTS+77wDixfD8OGS7ry0X3pJGvbf/lbigi5bJg1wZpLBK558UgLIv/OOOCAe6Fx1FTzyiDyLZ54Jp5yS+LE77SSTefbYQz6uQGxAX3wRNm4M3/fMM2HGDJg1y7UTXbwY9twTDj9c0mbO7IILUlLKihWiKJ1xhig5zz4Ld98tNriOr63qalG0v/oKLrkERo6Ehx6CDRtg9mw3r0svdd8N990necSaiHXJJeK1/oknZGLXqlWusnTxxRIxwWHQIFG0cnPlw3/wYLjrrrZ5GiNyLlwoMvY/YvZIJPLGjPZpq2/aUA4/PHV5e73yQHz3nazvsUfHvkyjyXbFFYkdu9VWrhI2caL8Pv008XMrfYMdd4S5c8PTMjLEWNth0iT5COgq1BGjy29+Az/91PHjMjKk18FRvByHydHah7POalvmfr+8FPfeW2ZAK32DceNgv/1EwU5Lk97Nc8+Nvq/PJ/XBGBg61O2lArj1VlG89tpL6mCsPECO33prWf7Zz8K3Pfpo9Pqz777yAbBwYex8B6iT70TGC+YZY24zxkw0xkwwxtyOGNwriqIkhypenaezdjd+v06aUZQeIJGeq3OAvwLPBtffBBLsMlG6jHnz5Eunq/jwQwk1tG6ddA3/+GNiszIDARlydIz9a2uhoECWq6vh88/dfXfbreu/aEpLxV7N64W1a2W2Z7LxLpubxUZmyhT3Gvo7334rw1pDhkjZvfKK2Pds2CD+6JyhRr9fyuftt92hZCcQfHNz8kONjr+qrnDIuGZN+LT2iRNhk03C9yktleHznByxYRk7tvPnTZZPP5XnZexY6UH86ivp6Ro/Xsr7ww9luaVFhnH23DP68+jU2x9+CFec4j27DQ3yjIP4/KqsFJchzhDRF1/Ivd5lFxki6k98/73YyAKMHp18e9FdrFwpNldffin2lVtuKfVl3Tp4/XU3Bu9XX7nPkd8v9w9clyJvvx3fBZBjs+Xk0dQkxxgjz5ITA7i2Vmw9t9suej6RH1EVFSKnM9ry+ediHjNmjNiSFhSEP/8ffwxHHNFxO69PP5VnenwfnP9nrU3oB+Qmum9X/3baaSc7oHn/fWt3393a/faz9uabO5/fmWdaO368tZdcYi1Yu8ce1ublWTtsmLVTpkjanXdGP/aRR1xZwNqHHnK33XOPpO23n7XDh1v72GOdlzUSsPbYY60991xZhuTz+uADOf6WW7pOvt4OWHveebI8a5asP/ywtaeeau1ee8m9228/a488UsrFudd/+5u1e+8t+7/3XufOD9a++Wbnr+WUU6zdckuRb5ttROZILr7YPefhh3f+nMlSWSky7LSTtbvtJmlg7a67Wjt7trX33+/Kud9+1mZlWfvjj9Hzeu89d9958yTtzjslj1gcfbS1m20Wfg5nee5caQNGjrT2wQe79rp7A3l5Uj923tna7bbraWnis3y53JMbb3Tvz+67y7bTTnPTvvxS7tlRR1nb0mLtYYe5z65zb3ff3doDDrD2+++jn+uyy2S/e++11u+39tBD5dihQ60dN849F1g7aZK1r70WPZ8HHnD3O/lka/fcU5Y//1zqana2W+d8PnffHXe09owzrB08uOPtQVWV5HHwwR07rnuJrU/F22hF4doD+AFYGVzfHrinveO68jfgFa9UcN111l54odv4P/OMtSecYO0OO1j75z8nlsfpp8tD5/B//2ftWWfJ8imnWDtzZldL7SqKJ5zQecXr3Xfl+Jtu6jr5ejsgjZ211j75pKyvWJH48dOmSbklS2amKHhvvJF8Hg4nneQq9y+9ZO0RR7Tdx6njIC+nnqK83NqCAms//VQUAGutzc2VF4gDyIvJWmsnT7Z2wYLoeb39tnzYgLVffNExOcaMcZ8Z50Xu8LvfWfuvf3Usv77A2LFSx7/91tqtt+5paeKzdKnck+uuszYnR5Z32UW2nXpq59u8RDjqKFHCQT4AwNoffoh/TF6eK9ePP8pyS4usO4qZtfJx71zD/PmSdsghsZW6WFRUSB4HHNCx47qXmDpNIgYCtwOHABuCPWRfA9PiHqEoqcTaATf9uF9hjMZ+6yo6Y6PVlZMllK6hP7Zroc6/FSAx43qstasikvSJVXoOVbz6Nh6PGtd3FZ0xrnfcDyi9ByciSX96PlTxakMiT+0qY8wegDXGpBtjLgQWpFguJdVkZopTTBAjxcxMUWa+/DJxw+mMDPjd71yj/Msvd4+dP1+mOY8aJb/NNhPnfD6fmxb5++c/45/v+uvlf+5ccR7pcMQRHbt2B8dZ4CWXiMGzY3zcXfzvf9HL4ZprOpbPrbe2zSOaSwHHsP3xx+X/xBPlv6NhZvbd152eHnpOpx7Eur+jRkkvS16eGPQ7+3bUHct118lx//mPW98yM+GNN+Q+Oq5X3nvPreMgDh2TjVXXUfbe2y1fgKeeEoP2zEwxNh41SgzeQ33gZWW5ge8bG2HbbWW/2bPFuNkxnneM9KHjk1caGtzlyZPDDenr6sTf16hR8ePAtkdNjRhnO/f86afb7jN/vhhb33tv8udx+Pvf4eyz3fWFC6UezJgh5b12rdTxzEwxtH/00c6fMxVMnOi6bLj5ZjEcB5n4MGqU+Fnzet06kirS0sSQH1wnx+29E/bZR9yahO7r1NfQSWGhYYuc/TZudJ2thuL3Q1FR9DbGmSDx/vvi0+zAA+W/rxBvHNKK1j0UeBIoBkqBJ4Ah7R3XlT+18UoBLS3Wrllj7cqV1q5da20gYG19vaQ1NSWWR0OD7O8Ycq5ZY21jo2z73e8kfc0a+RUWil3WBRe4aaG/q692jb5j8ZvfWHvFFdauXi3HrFghEw822yy5Mpg2TQzM16wRo+ePPkoun2R55BExeg4th9tvt/bXv+5YPmedZe2117p5vPee2AhF8sYbYsg6dqysg5RlR6itde/5W2+Fyw7WHnNM9Pvr/DZssLauTgx6wdqPP5aJHh3h1FOt/cc/3HprrfyvWyf2f3PmSNpTT4kNSFWVnHfYsNgG610NSJ13uPJKa88/X5aLi92yCKWqytrqand5zRqxlbzrLrl3jp3MTTfJ5JKKio7LVVkp9mbWikF1ba27bc4c2zrZZOjQjuftsH69tUOGiPx/+EN0G8rZs+Vcf/lL8udx2GSTcLun99+X9aFDxZ6utNTd9oc/SFvTGwFrlywRA3unHa6okPpy9NGyfc4ca5ubUyuH8x4oL5fljRvbP6axUd4HDjU17nIg4L5TKiutXbVKrsnhppusPfvs6HmGvkPWrLF2xgw3raLC2ueft3affSRtxIgkLjalxNRp2nUnYa0tA05sbz+lj+H1tvVcn5kZ3Zt9LDIy3P2dL5LQbU66c77Q3pBICgrEpUE8PB7x2jx6tJsWy9NyIvh88vU4alTPDV1mZ4eXx+DByeUzaJCbT1VV7P3y891r9fnki7IjZGe7vSxjxrS9l7HubyTOeUeOTM7ea/DgcG/XxsjXeWTvXVGR9LBB6nsKIokMZO+c3+nJiMSR01nOy3NDfUWWUWZmci5QQsvA4wnvMXPqRVdE4HDal9BrCqUrbfziudEYOlR+oeu9Faf9zcpy05x77Nwnrzf1kUIi3wOJjIBEtp+hIeqMcZ+F/Py2z2FWVux76IyQODjl4aSFtpd9KIJKu0ONQaepLxtjSo0xJcaYl4wxE7pDOEVRejFdoawa07/sWVJFfysjnVyhDGASURGfAu4GjgquHwc8DeyaKqGUPkikrUnkV1JpqTiI3HHH6McbI3HAcnLgwgsl4G8ojY1i6+LYEYRSVib2PkcfLXk0Nrqx62prJZi34/A1lMWL3eW6OolbuOOOEv/uuOMk/eabpSfOGLFnKyoSm5KWFlm+8EI3j1tuEVl+//v2nfr98AP8+9/Rv8BnzZKwHPPnw+abu9cSjXnzxAHhFluEp5eUSAw1EOPWggKJtdbYKPHeLr6488bV0RSvRJUxxyjc44n/Eo52/z7/XMKcxOL++8Xe64cfwm2YysvhssvEpm+nnaIf+/zz4pjyiiuiX8sbb4iTyYMOEruSaDjGxMXFEg/v2GPFni/W/vGwVmLxOWG+rr8e3nxTnBOnkrIyt/6sXCm2R5dfHt4zUVUFt90mAZlDe7ZKS13HrNbCPfe4sSOtlfzmB4OfvPuue56ttpL64NjoARx8sDh6jsfy5fIfCMjxTz4p6zU10fd/803XRg6kjjQ1Sd3IyoI77hAbvdA68t13cr3R2p/Osm6d2LdGa6Mc+vsM1CefDO/pg+jXHO2ZfO89+V+9Wp7Nb74Ru75ovWhz58Kuu8Lf/tajjoITCZL9qbV214i0udbaFD/5Lr0+SPZA5803xdgx1Gv4+vXyAps+XdaffFIejKOPFiPSSO68E847T7qL335bGr5Qli2T41atkiEuh6YmN+j3229LV3wg4PYQfP89bLONa0gfitcrwZ8HDRLjzm+/FaXkyy/hk09kH2Pkxf+f/8Af/iAGz7vuKgbeV18dPmPH65XtF1wAJ58cv8zuuAOee04M40NfomvWiKL13XeyvNlm4QpiJFdcIQrcm2/KMCyIcnXPPSLbvHkiO8gw3HXXiZxNTXLdZ5zRfrSCaLz6Khx2WPixL78simvoUHAsamrgtdck0Ppuu0ndiMZ330mZRt6/E08MrwcOb74p989hr73EAzyIQnvfffLyvPrq6OfbZx8x2P3pJ7c8Qzn+eCnTPfeU4NbRcAIZX3ONGHIfdJBMaPj8c/FC3hHmz4epU0X5OussSdtpJzG03nbbjuXVHtXV8MADEsT82WddD+iXXCL/tbXhH1jz5sHOO0t5T5nipr/+uhi7L1kiz1So4fQrr8gHmINzX9etgzlzpG7us48MJc2bJy/H9oJ3O3Vw40YZeho0CM45R47deWfYf3933x9+kHrq8N//unEGFy2SDx1jpE259VZ3vzFj5HlMRc/jyy/Dr34lSuzvfx/9efzkE/joI3leu3vIPNX89JO0YdEoKgqPL7lsmbTpzoSqykqps2VlUtdGjICHH5Z2MVKx8vvl42H//eUjPnQ4NDXEbljjGYAFlbKbgEuATYFNgIuQEEKFQGF7x3fFT43rBwD//KfrsO/999tuX7o0thH2nDniFdnacEd+1lr73XfWbrVV4nJ8/LHrWdxaySsQEGeTDz7oOsBsbHSdXTp4PNYef3xiHvtvvz32ZILrr7d29GiRo72JA1dcYe0118TePnOm67CwqMjab75pX7buZM0acdYYi652ennVVeKFPxYHHSRl9dNP0bcfd5x4yz7llNh5LF8unr+tFeP+3//e2l/8IlmJxcGxMxkB4nuoTwXOeUMN8a0Vz+SOF/VQXntNnGJG46KL5Jittgp/Tr/5RrzLb7+9m9/DD8uEmvbw+cQ7umMEPn68tBeJ8Je/uNe3aJGkQVsn0uPHh8vblcye3bNRFfoLDz4o7XROjjtJJZTmZrmHL77YXRIlb1wPHBv8PzMi/TTAAmrvpXSeZHpdlL6POlNNjP421NTX7rn6DVS6kERmNfbBCJRKn2bpUndm1dixsuzYjMRixQrpanYCv377rfwvWdLx869a5R7vEAhIWqgdgpPmYK28IJ1g0rGoqhIbhHisXy/XvnZtW1nAtY0qKYk+5BaNUB9OvQWPR4Y9v/1WZjBNnuwq4U1N8YdZk2XBguhlCq6foW+/Fbu44mKxvRk3ToZ4VqwQe6YlS6SeOsPma9e6s3JD761jUxdpv9IR6uu738dcNL79Nnyo0Xm21q1zhxo3bJCho1g4Q3XRhuzWrGk7y3bjRvdeTZzonr+uzg3iHAjIekOD5OvYkyVCqByLFrmzpIuLw+uIc2+dtM03b9//nd8vdS30HIWFUp9Ch7FLStTJaFexZEmf+EiJqXgZY3YGVllr1wfXfw38ElgBXGWt7UDtVpR22HxzmWY/bZprW1FeLvYxM2eKw81YNi2jR0vjdsIJYhtSWirLDrt2YB6I45LAOX5aMDrWlCligzFnjtigeL1imxR6nr32Enuciy6SX6yX7WWXic1OqBPYULbYQq5pv/2kwQ89B8gLZskSUcwKCsRQNBYTJsCQIfLFXljYcfcRqSYvDyZNkmtcvFhsmhwnko88Ig49Q208OsuWW4oNSGSZOng8YuR96qki19y5Uqd+/nOR45NPZMLFPfeIklheLtcwZYoc65SvY4S9xx7wzDNiS5gse+8NL70kEzZycqLbSKaSrbcWu5rf/Cb6lP2TTnIVk5NOkokEf/lL9Ly23Vbq94UXhitoRUVST3Ny3DKcOFE+Yk44QRSh884TGx0Qx7hXXin2m7vtJk5GH31UJsVUVibulmXrraVcBw+W59Lh88/D64jz4XfCCaJYX3ONa3MXi9dfF5tAx/a1ulqUcBBl0fnAuOUW2GWXxORVYjNxopTvDjtEV4o9Htkn1KFrTxFrDBL4gqANFxKbcS2ieF0L/Cfe+GV7P2A6sAhYAlzS3v5q4zVAefxxa088UZYnTrR28eKelSdRsrPDHQhGcuaZYrOTLAsWiK3CrFnJ59EbCbXvsVYch/7xjz0jy5NPWrvttmJf9sADEhD+o4/EUbDDoEHJ2RX1JxYsCHfWe+ih1r76atef58or5ecQadt48cXW3nijOMhN1qFyopx3nthotscLL4jDT4fXX3ftyRzHv9ZKnepu581Kd5CUjZfXur1axwL3W2tnAbOMMV8lq+gZY7yIe4qDgNXA58aY2dbaH5LNU1EURVEUpS8Qz4Gq1xjjKGYHAG+HbOuMi9hdgCXW2mXW2ibgGWBGJ/JTFEVRFEXpE8RToJ4G3jPGlAH1wAcAxpjNgMpOnHM0sCpkfTXqjFWJhs8nfn+2316cLkaGYOmt+HxiV+b1Rt++erX4F+pM/qH//QWfT/wZObZxZWVwzDE9J8u334rtls8njlXffTfcR5nPJ3Z9aWlyT/vb/UgEn0+M0rffXtaXLYPzz0/Nea64Al54QdbXrRNfW6HbL71U7AJTfR98PnGi3J5/scrK8Oc8VK7tt3dtvBYvHph1ZyATbxwS2A3xWJ8TkrY5sGO849rJ8xjgwZD1k4H/i7LfGcC84K86ZHlaxH63h2w7I04et0dsmxay7fXQY4GikG3zgKKIY19PtTwR23q7PH/uZfL0tvKJJc/yXiZPbyufWPIs72Xy9LbyiSXP471Mnt5WPlHlCab3Gnl6W/nEkyd4TK+Rp6fKJ/LXruf6rsYYszsyK/KQ4PqlANbaG7tVkLZyzbPWTu1JGfoqWnbJoeWWHFpuyaHllhxabsmjZReddoNkp4DPgUnGmPHGmHQk9uPsHpBDURRFURSlW+mMkXxSWGtbjDFnA28AXuBha+333S2HoiiKoihKd9PtiheAtfa/wH/b3bF7ieHNUkkALbvk0HJLDi235NBySw4tt+TRsotCt9t4KYqiKIqiDFR6wsZLURRFURRlQKKKl6IoiqIoSjehipeiKIqiKEo30SPG9R1l+vTp9vXXX+9pMfo0z79TzV3PlfP2PeP4++MbeP2TWl74+2gKcl3v6rc+uYFFK5pI95k2xzc0WuoaA5SW+5m8STplFX4O2iWH02cM6sar6D88N6eKe2dV8I9zh7HjFpldkuf+f1zJrw/L5zeHD4q6/a//KmVjpZ9LfzOEMcOie8qubwhwwZ0leAwsXN7E8CFpDM7zULzRz1H75HLC9IJOy3nzYxt4Y24t2ZmGvGwPe26fzdnHDCYQsBx8ziq22DQ9bP/8HA83/HFYWNr/PqvlmbeqyEx362ptg2VwnofMdMNXixtpaLQcc0Aef/jl4NZ9vlvayB3PbGSLTdLJy/bw+YIGRg1NY+maZtaVtfCL/fI4cp9crnmojOraAPvsmM3vfzGYrqCxKcChf1oNwFbjw6/xF/vlkeEzzHylknVlLfz9nGFsPSGDNaXNnHzlOvKyPVx0ciF//VcZF51cyPTdc7tEpmgsXtXEmTeub5WxrNLP/lNzOOPIQQCcd1sxpx1ewPabS719Z14t1z68IeyaSsv9HLpHTsy6GIvijS2c+49inr1hdNTtL71XzcIVTVz86yEx86hvCHD5vaU0tVjWlLZw1tGDOXCXnA7JoXSM0ooWjr1sLQAjh3hZt8HPZmN83HzOMAbnxYjgkQBPv1nFnM9rW5/z8uoAO0zO4MITY99/kLZws7E+qmoDDM714vUG24dcD7edPzxpeTpI2xdpkJQpXsaYh4HDgRJr7TbBtELgWWBTxPv0r6y15e3lVVZWlioxBwwNTYHWZWvhopMLw5QugAviVOYlq5q49J5S8rI93PWXETz9ZhXVtf6Uydvf8fvh2APzukzpAnj7nnFxt197ZlG7eWRlerjnohFt0h9+uQJ/IMoBSVBdJxn99fShrC1tYWVxc+s2C9z1l7bnj+TAXXLavEy/XtzAwy9XYozh8lOHsOd22W2O22ZiBg9ePrJ1/cwY+d9/6Uj+83YV6zd0XR0PhMxjinWNe26fzXm3FuMP7uwPwJhhaYwf5Wst/7KK1D53fr9l8rj0VhmffauKjVUR5wx5pew3NYf9pobfi8dfq6SpueMTt6x1I+lEY8Y+ee0G9s3KdF+uNz+2geYWnUCWavzB6pGf4+H4Qwq47amN3HLusDbvmI5y/MH5HH9wfuv6a5/U8M3ixnaPy8ww/PP84WRluoN63yxp4KGXOhPtsOtI5VDjI8D0iLRLgDnW2knAnOC6oiiKkiAt/v6rSFhLnH6CjuP10mUfDIrSVcRUvIwxY40xzxhjPjDGXGaM8YVse7G9jK217wMbI5JnAI8Glx8FjuywxIqiKN2AgaAm0Lto6ccdzV2sd+E1hkCg991DZWATb6jxYWAWMBc4HXjPGHOEtXYDsEmS5xturV0HYK1dZ4wZ1t4BStfwyCvSxXri39ZSUe1n6pYdG+LypRk2VPopGiRdx2leeOatat77sj5sv5/tmcMJh7h2QOffXkxJefibYtTQNG45d+De+tc+ruH+Fys4cXp++zv3EhqbLE+9XcWsd6rJyfJw9P55HLVvXsz9X3i3mkdeqSQ320O6z/B/Fw4nN8vDVQ+U8vE3Umd8aYbGZsvilU188FUd980qx5eW/GvXl2b4dkkjGT7Dkft0jQ1U5Cv70rtLWLyqiYx095t16/HpXHbq0HbzWlXcktA5A9by59tLGD4kjZYWS262B4/HcNUDYnLx1BtV/HbGIP5w03qq6gLc9qdhDC+M3pRfcV8pC5Y3kpnu4ch9cjnmAKlzK9c3c95txWQHh2JO+VkBmemGf71QQVOzZcQQd4jIl2Z4bk41H34t962sooX0du6TL80w8+VK3p5X15qWk2W49+IReD2Gp96o5PVPannsqlFhx1lr4w41dhSPB+54ppwWv9jRJcO9s8qZv7CBYYO9bWwNAV75sIYHXqwgN7ttP8bJh+a32uN9u6SBP91ewvV/KGK3bbIAGf664t5SXvj7GM68aT019QHS0ww3/rGI0THsMHsTNz+2gbnf1VNZI92KHgO+YFVM83alCi34vIb3v6zjmyWNjB2Wxk1nD+NPtxWzprSF3GwPk8b4+OyHBpqbLR5P22MXrmjkxL+t5Z6Lhnd6GLQzxFO8iqy19wWXzzHGnAS8b4z5OW3boy7HGHMGEu2bcePi264o7eNLM9xybhFDg4rTiCEdM+8bNyKNp68bRUbQ8P6offPYY7ussH0++LKeFeubw9J+XNnEnRcMJzNDjqutt/zlzpJkL6NfsKqkhcP2zOHkQztvqN5dTN0yk+fmVDNmWBo7bZHJquLmuPuvLG7moF1zOGrfXM67tZi6+gC5WR5+XNnEtb8fysTR6Qwv9GItfPhVHWtKW9hhi0xO66AxdihbbprOU9eOwlrCFIdkMVE0gEUrmrj8tKGt+a8qbuGBFysSyq+q1s+Wm6Zzwx/j29rd+MdhVIbYT+ZmefClGX47o4CTr1znyrKyiaJBXiprAgwvjJ7XohVNXHzyEBatbApT/DZW+YOTFop44Z1q1pY2k5XhYYfNMzj+kHzyc9zy+/m0XHbdxv1Q83oMwwvjl+8v98tj7ynh7cNvrl5HIABej7QLq0uiKKK2a3u8nJfvi+9VJ614PTenGoBla6LX+VXFzey7Uza/OjA8/5feqwm7xpJyP9ZC8QY3rbTcT029pcVvWbammYeuGMGNj26gvDrA6D7wbfrjyibO+dVgrnt4A2OGpXHXX4aTk+Vhpy0yycnqekum/aZms+X4dKpqAlxxXykAC1c0cet5w7jy/lK+/DHAmUcNYpets8I+jgC22DSdR68chT9gyU2BbB0h3tvXZ4zJtNY2AFhrnzDGrEdiLCY7RaTYGDMy2Ns1Eoj5BrbW3k8w3MDUqVO1r7iTZGUYxgzzMaQguReSMSbsqzrNaxhdFP5FNiivkRXr2x47cmha65d1lRrkY4BRQ9Kizh7trTg9UUMHeSnI9bYayMfCID2bo4t8eD3udWakexg1NK1V8U/zyoscC3lZHgqTrJ8gdbSjHxTtETnSmJlhGDU0jZFD5TyNTR1rmrIzPe1+aedme6L2nmRneiLWDb526lBmumFUURrFG1sojeh5Lsz3MrrIR16OFyeCSU6Wp81zHe1Zbw9fWttjQvXYaEotBL/ou7DLyxOse6m0izNG2rjI6y3I9VDXkPh5031St9rrTexNpKeJzABFg72tCvvQQamZt+f1SL3KynDrcnaG1HGvx5DmlXYn2nsuFe1DssRT+x4Edg1NsNb+DzgG+C7J880GTgkunwK8lGQ+iqL0A7wegz9gCViL8fS+F45+8XUv7c1q7ChOlWpObJS3SzFo/VGiE1P9s9beHiP9S+Cg9jI2xjwN7AsMNcasBq4EbgL+bYw5HViJKHFKEjQ2BXhujtjbVNaI9l+Y7+Xn0/J4e14tKyOG/Ooau6cJ+PDrOoa/4n5tRJvKXV0X4JFXKlrXp26ZxTYTM7pDvJRSUx/ghXeq2WSkj2k7ZLOmpJm3Pqtts9+3SxvZbeuucyPRE7w7v478nLbfbTtukckXCxtYsLyp1VdYZW2Ax/5bydBBXiqqw3tdPJ6gm4UuHmLqCoyBdaXyxq5rCDDrneqoPX0/rW0Oq88AW2ySwebj0nn5w5rW3qS1ZV339n/3izrqGiwZ6QGeeK2SCaOj90hV1bryvv9lHYX5nqg+wF6fW8vYYb6Y+XQFLX5YXdLM+FHpfLOkARCfW6FT/tdvaKGuvuumIXqDWW+o9LNweSNbbJpYO/PRN3UsXtnUJv3HlU1sPi6d97+sY8/ts/hxZRP//l81Zxw1qG0mhrAu0x+D+c39rp7y4HPw01pppx99tbK1V66x2fLUG5VsPi4dX5ohEIDpe+RQkOPluTlVbD0hgymb91z7UVnj56X3awgELGWV7vPc3cptebW8RxqCvc6VtYE+M+M3Zf1u1trjY2w6IFXnHEisLWvh4ZfFYH6XrTLZbGw69z1fwc+n5XHf8xVsOzGDscPd23vyoQUMykvtuPaUzTP5ZYQdxW9nDCIrw32l5mZ5OOOoQTQ0SuP648omNlT6+4XitXxtMzNfqWTyJulM2yGbzxc08NL7NcyYFv6i23FyBnts39bHVG9m4ph0Tj28gK0mZDC80Bt1yPjDr+v5ZnEji1c38cv98tglqFyef/xg1gWVjl/ul9c6NAHSIxEIiPLV2zq8crM81ASVgJXrm5n5ciW//XlBq50kwKiiNE47oiCswV9b1sL3y5qYMS2Xt+fVst9Ocq9HDU1j/506V8/vvGA4D75UwR1PbyQrw3DOMYNZvi62vd2vDsxj2OA0dtoyiw2Vfj76pp4RQ2So1ynuvadktSqHe0TxfdZVTBjlY96CBsaPSqe8Ssp1VUkLm49zHa+uLW1h9LCuey3tNzWHtWUtLFndzJx5dQkrXk+9XkVutoctN01n+u45NDZb3plXx5zPa9l8XDpXPVDGnRcM5425NQBsuUl6mzyMMWFD1evKWthuswwmh+w7fpQMxaf7DGceNYh0n+HXhxWwdHUTjU2Wx/5bBcDoYWlMHpfOQ7Mr2W9qdo8qXotXNfHoq5WcfGg+P5+Wy7gRPs751eDW4ffuoCDXw++OHERjU4BfH1ZAfo6H848fTMlGPxPGtL0XvY3eMeCpdIpdts7i0D1yeP4dMQLNyfJw0qH5jB/VvRVwxJC0dj1VezyG4w5yZ/O98mENi1a07xCvr+CN0G333TG7w967eyN52R5OPsydDBDtmtZv8FOysYW9tg+/5kN2iz3D0OOR6f5i29N18nYFQwd5w+xtJo9Lb+O5PzPdw0kRkyTmLajn2bfkWRw33Nel93+biRkcuW8edzy9kcnj0tl3p8TMbUcNlWezeGNbhXnC6HQmjE59W7HDFqEG+rDJiOi9axO7UJZJY9O56ndFPDenqs3s6nik+wzHH5wfpuBMGlsV1mMbsLbVVm3s8LbXEjnU6PXKbNv27tneU7LZe0o2lTV+nv1fdZvtKZgs2GF2nJwZVq/jzXBOBV6PCXOsCvHbmd6GxmpUFKXH8Hqkx8v2wh4vRYmHvz09zvRKN3BKL6DdHi9jzHnATKAaMbjfAbjEWvtmimXrNpqaLaXl7gB1Xo4nbDp1NCpr/NSE2HuMGJoWNnsrUSqq/dQGhzOKBstMt+KNLeRmeWJOx21usVG/3hqaLGtKmvtciIyScj9rStzhkqxMD4X5PedjJRGstRRv9OMPGWJy6lB5lVyPY78xUAgELKtKWiganPi983jEpqWqxh91Jl9PsyF4LzvSWwIS33BDZJidLiTUdqsjBAISvzAVPpbaw1oZSlxT0oxF7vuakmZyMl1ZyqtTV2brN7RQWx9oNeCPbF9r6wNU1wUoyPXQGCPc0foNLa1t1bqylrB3QCQGqKkLtO5f35C87Vpzi8Vpasoq3fYyYGVbhs+QneVhcJ6XFr+lvNpPUXBmYWl5C/WNYkTpNZCZ4Ul6drs/YNlQqTPTO0siQ42nWWv/aYw5BCgCTkUUsX6jeK0ta2n1CdLcbCnM93LvJfHjxZ123TqamiwFeV7Kq/ycd1whB+/acS8bZ/+jmBa/paHRcsReufx8Wi7HX7GWPbfL4trfR/f38+J71dw7q4KsDIMvzTByaBq+NMOIIV4uvruUNA8UtKM49haGF3rF19DdUv42YKmuCzD71rE9LFl8Vqxv4Ywb1jEswnHlTltksni1XM/a0hY2G9P7nSB2FcOHpLGhso7sDvjIKcj1kuaFud83cMphBe0f0I0U5nuprg201s0pkxKzDxpS4KWy2s9zc6rZb8eut5kaNthLQa6H8UkYwm86Kp1n3xL7pV27eYLH2OFpzHy5ks8XNLDZmHTGjUjj//5d3qa+/CIFw1ajhqbx0OxKbnpsAwuXN5Gf4+GhK0aG7XPdw2V8+n0DUyZlsHB5U5uPv1FD03jmzarW+vDgSxVkBIM3h9qxOowuSuPlD2v4eonsb5CP60TJTDcMG+ylpNzPsjXNrT1sX/3YyMlXrWNUURprg5M/hhd68Qfg3zeM5tUPa/jns+W8fOsYLHDs5WvdayhKo2RjC6//c2yrq42O8O78Om5+bCNH79+9Q4v9jURqgXN3DgNmWmu/NrGcsPRRNh3p44mrxYPy0tVN3PjohnaPyUz3cOcFRYwu8nHL4xtoSbKXqaXF8s8LhvPO/Doqqv00Bx+u5jizM1r8cNxBeZxx1OCw9KeuHZ2UDD3Jzltl8eItY1rX6xsD/PLiNT0oUWK0tFg2GenjgctGxtxn/z+uDDOk7e8cvX8eT75excG7JP4BUpjvZebfRrW/Yw+wyUgfz/99TPs7RjB+VHpSxyXKVuMzeCHJ/CODDncnM6blMWNaz7yw99w+m6t+a3jp/Wo2VPqpjdL75LS49U2WC08sZFyEDdq0HbKZtkPiivReU7LZa0ryindGuodnrh/NY/+txO+3eL1w8K45vPlpLeNH+XjoipFcdk8Jc79r4J6LR3D6teJc13l3+AMWayVwdVVtgOxMwxNXj+KAs1YmLVOL33Lwrjn88ejB7e+sxCSRT9P5xpg3EcXrDWNMHqBhRxVFURSlB9EA4H2TRHq8TgemAMustXXGmCHIcKOiKIqiKD2Ev4/4rVLCSUTxssBWwOHANUi4oL7t/TEOHo/E5Drr7+GxbxYsbyLNK9OTQQLFOsb0xgO3PbWRVz+qYcHyJrbc1B1eOuVnBeyytRuzbG1ZC1fdX9oagmVjlR+vR2Z3/ft/1cxfKI4FP/u+gftfKG8znAjw1OuVHNFDXfapxmNkksBZf19PXo6HG/9YFDO8SDS+XNTAgy9VYAxc/OshUad5d5avf2zg709sJC8BY/BhHbDp6Os4z4O3N8x3V5QIvF749HtpXxsaLb+6bA05WR4uOKGQO57ZSPFGsZdatKKJI/fpSUnD8XrgkVeqGDk0rdXO0HFeLDamjXiMOKY+6+/rWyd1nHdbCelpct2bjvRRkCvHWAtn3VLcOot4+u45fPZDA6XlfqZNyWrjMsVh/YYWHppdyc5b9tvXf7eRyFvhHmRocX9E8aoGZgE7p1CuHmOTET7uu2REGw+4Z99STIsfzjpGFKE0rxso9nczBnHYHrmUVfi56oEyDtwlh8mbpPP8O9X8tLY5TPEqLW+hvDrA1WcMBSSm2dBBaRy+Zy5bT5CHKi/bw7dLGvnvxzVRZaxtsBydZMDX3k5GuocHLx9BQ5Pl7FuKO3z8yvXNDCnwUlrup3hjS0oUr5XFzQzO83DFaUPj7vfsDaP6zCSHriA328PDfx3JuOEDR9lU+g47TM7k7ouG4/fDguWN3DurgrIKPz+ubKKp2fL3c4aRn+2hsjbAZr3ICeeR++Sxw2RRdkYNTePkwwpajfr/dFwhf/jlYPKyPdx3yYjW2ZihPsTysj3kZXtaA4Y//NeR1AVt3D76up5FK5v46Ot6Tjwkn+9/auut36Gswo8/YPndkYNScZkDikRayF2ttTsaY74EsNaWG2N6T63sYjweE+ZJOZSRQ9PYanzbmU0FuRI42PHOvclIH1uNz+C9L+qi5jOqqG0+WZmesLR408V9afTKqfddRWedORbme2X6dAqZMDq93YCrRSkKFNub2XTkwJnFqfQt0ryGLYOe6yON68cMS2vd1tumKOVkeaK+d0AcvaYHA6VPTFBZDH1Gl61pprLWT06WYezwNJa14wJn1NA0BuUNnI/JVJHI27vZGOOlNXC8KUKN6xVFURRFUTpMIp/kdwIvAMOMMdcDRwNXpFSqXkpkOJhEmDOvtvWLBGBNaeKRRH/4qYm359Wy/1SZnt/YFOCNubUDaibLyx/UMH333LAyjMfydeKcsanF8vIHNawqbmH3bbMYMSSNQMDy+ie1Yc4Rc7M8HBTD/1pZRQsffFXfJv27pY0xndsqitK3eOHd6m6NM9jb+PS7htYA11/92MAL71YzeZP0sF62yho/H3wVfQRH6Tjt1jZr7ZPGmPlIcGsDHGmtXZByyXoZl586pNU4MRZDCrwcc0Bea1fuHttm0dRsWVkc3n37sz3a93M0ZlgaR+6Tyw2PbGhVvJava+aOZ8o58ZB8fAOgnfj1Yfk8/HIlW0/ISLgbfW1ZC1tums52m+Xy3dJGXvu4hsYmy3EH51NRE+AfT25kxj5uTK/Z79ew707ZrZMdQpn7XQMPvFjBIbuH36/8XA97bJvVZn9FUfoGE0b7OHKfXD79voE1pS0csHPfClrfVWw7MYO9ts8iK8Ow3aRMDtkth0+/r+eLhQ1hDry/WNjArHequfDEwh6Utv8Q8/VtjMm31lYZYwqBEuDpkG2F1tqN3SFgb+GAndtXltJ9hj/80p2FuP3mmWyfZBT5glwvZx09mNkfhBvYbz4undNnDEoqz77Gbw4fxIdft+1xikea1zBpXDp7bpfNATvn8K/ny8O2D87zcN6xbuPx6ofRJzA4HLBzdtj+iqL0fYoGpXHusYW89nENtzyxkV23GZgfUpuM9HHecW77du6xhXz4VR1vzK1ts++0KdlM373vBKLuzcTrN3kKcSExn/Ag686EiQkplEtRFEVRFKXfEVPxstYeHvwf333iKJEEAlATDKKd6pl6vZW6hkBrGYDY2mVleLDWUtsQXibRHAo2Nltq6gOtU6hDafFL/tmZbQPjNjQNIGM6RVGUIPWNAeobA2RliHmNPwABOzDfP6mgXUshY8xRwNvW2srg+iBgX2vti6kVTTFGXEccf7kbu9Dx5zJQGDbYy2X3lIal1TdZnrluFO/Mr+PeWRXkZLr2WcbAySHBlocM8vLoK5X8Z04VAGMj4q/lZhl+cfEanDYlNC+AYw7smbh2iqKknsHBQNiJOEMeKAzK8/LFokaOuGA1T183iqJBafzv81p86hi5yzC2HS3WGPOVtXZKRNqX1todUilYKFOnTrXz5s3rrtMpvZwT/7aWv59TxNuf19Hst5x2xKBO5ecEmv3nn4ex7WYDS7FVFEWJxklXruXms4oYPczHtQ+Xsce2WQnZOiutxNRUE1Hzo+0zAObUKYqiKIqidC2JKF7zjDG3GWMmGmMmGGNuRwzuFUVRFEVRlA6QiOJ1DtAEPAv8G6gHzkqlUIoSj9r6ANc8WMZbn9W2BnrtDE4Q7o4E41YURenP1NYHuPqhMi67p4RvFjd2SVurCO3aeLXuaEyutTa+06MUoTZeSig/rmyirEJcLW+xaQaF+Z2LHbaquJk1pS1M3TKTNDUgVRRFCWtnjTFMmZRBVqZOQugAMV8miRjX7wE8CORaa8cZY7YHzrTW/rFrZYyNKl6KoiiKovQhOmVcfztwCLABwFr7NTCta+RSFEVRFEUZOCTUb2itXRWR5E+BLIqiKIqiKP2aRNxCrAoON1pjTDpwLjDggmQriqIoiqJ0lkR6vH6PzGIcDawBpqCzGhVFURRFUTpMuz1e1toy4MRukEVRFEVRFKVf026PV9Bp6svGmFJjTIkx5iVjzITuEE5RFEVRFKU/kchQ41OI49SRwCjgOeDpVAqlKIqiKIrSH0lE8TLW2settS3B3xNAYl5XFUVRFEVRlFYSmdX4jjHmEuAZROE6FnjVGFMIYK3dmEL5FEVRFEVR+g2JKF7HBv/PjEg/DVHE1N5LURRFURQlARKZ1Ti+OwRRFEVRFEXp78S08TLG7GyMGRGy/uvgjMY7nWFGRVEURVEUJXHiGdf/C2gCMMZMA24CHgMqgfs7c1JjzHRjzCJjzJKg/ZiiKIqiKEq/J95QozfEcP5Y4H5r7SxgljHmq2RPaIzxAncDBwGrgc+NMbOttT8km6eiKIqiKEpfIF6Pl9cY4yhmBwBvh2xLxCg/FrsAS6y1y6y1TchsyRmdyE9RFEVRFKVvYK2N+gMuBz4CXgK+RPx5AWwGfBTruPZ+wNHAgyHrJwN3RdnvDGBe8FcdsjwtYr/bQ7adESeP2yO2TQvZ9nrosUBRyLZ5QFHEsa+nWp6Ibb1dnj/3Mnl6W/nEkmd5L5Ont5VPLHmW9zJ5elv5xJLn8V4mT28rn6jyBNN7jTy9rXziyRM8ptfI01PlE/lzlKmoGGN2QzzWv2mtrQ2mbQ7kWmu/iHlgHIwxxwCHWGt/G1w/GdjFWntOMvl1FcaYedbaqT0pQ19Fyy45tNySQ8stObTckkPLLXm07KITd8jQWjs3StqPnTznamBsyPoYYG0n81QURVEURen1JBIyqKv5HJhkjBlvjEkHjgNm94AciqIoiqIo3UpnjOSTwlrbYow5G3gD8AIPW2u/7245otApFxkDHC275NBySw4tt+TQcksOLbfk0bKLQlwbL0VRFEVRFKXr6ImhRkVRFEVRlAGJKl6KoiiKoijdhCpeiqIoiqIo3YQqXoqiKIqiKN1Et89qTIbp06fb119/vafFUBRFURRFSQQTa0PKeryMMQ8bY0qMMd+FpBUaY94yxiwO/g9OJK+ysrJUiakoiqIoitJtxFS8jDFbGGNeM8a8aoyZaIx5xBhTYYz5zBizZQJ5PwJMj0i7BJhjrZ0EzAmuK4qiKIqiDAji9XjdD9wDPAG8jQSOHAxcC9zVXsbW2veBjRHJM4BHg8uPAkd2TFxFURRFUZS+SzzFK89a+7K19mmg2Vr7jBVeRhSwZBhurV0HEPwfFmtHY8wZxph5xph5paWlSZ5OURRFURSl9xBP8fKGLN8WsS09BbKEYa2931o71Vo7taioKNWnUxRFURRFSTnxFK+7jTG5ANbae5xEY8xmwP+SPF+xMWZkMJ+RQEmS+SiKoiiKovQ5Yipe1tp/WWtroqQvsdb+KcnzzQZOCS6fAryUZD6KoiiKoih9jlS6k3ga+ASYbIxZbYw5HbgJOMgYsxg4KLjet6h8ChYa2HBzT0uiKIqiKEofI2UOVK21x8fYdECqztkt+NfLf8vanpVDURRFUZQ+h4YMSgotNkVRFEVROk67GoQx5jxjTL4RHjLGfGGMObg7hFMURVEURelPJNJ1c5q1tgo4GCgCTqUv2mZ1BU1LoPp5Wa6cCevPgo3/17MyKYqiKIrSZ0hE8XICPR4GzLTWfk2c4I/9mvpPwTbA2Neh6EbwDoXyO3taKkVRFEVR+giJKF7zjTFvIorXG8aYPCCQWrF6MemTIecgGHwWFJzU09IoiqIoitKHSGRW4+nAFGCZtbbOGDMEGW5UFEVRFEVROkAiPV4W2Ao4N7ieA2SmTKLeTP0nYP3haYHKnpFFURRFUZQ+RyKK1z3A7oDjl6sauDtlEvVmKu4G3zh33VMAfg3grSiKoihKYiQy1LirtXZHY8yXANbacmNMyoNk91ryQ/zCegcBvp6SRFEURVGUPkYiPV7NxhgvMuSIMaaIgWxcb5t6WgJFURRFUfooiShedwIvAMOMMdcDHwI3pFSq3owqXoqiKIqiJEm7ipe19kngIuBGYB1wpLX2uVQL1q00/gCLh0HTT7H3qZ0j/57ckEQP0AzVL6ZQOEVR+jQN38DSSbB0ovtbf3bn821eDYuyoXFR5/NSFKXbiGnjZYzJt9ZWGWMKgRLg6ZBthdbajd0hYLfgmwgmK2goPz76Ps0rIHt/yNzBTTNpkH8yNC3tFjEVRemDtKyFtJEwcqas138KFf/qfL7+MrD10LIOMiZ3Pj9FUbqFeMb1TwGHA/MJ2ncFMcH1CSmUq3vxZEDa8Pb3823SNi1tOOBvm64oiuLgyYb0ibLcsrqLM9f2R1H6EjEVL2vt4cH/GF1AiuBlIM81UBSlp1HFS1H6Eu3aeBljjjLGFISsDzLGHJlSqXqCQC3UvQ+BBii7HkqvhLJrwF8JTcuh5qXoxxlPW6eqitIbaPxe7Bd7Cy3rofQqebYGkl1Sy0qwjeFp9e9D8fnJT9ZpKXGHKwd6+1P3EVTc39NSKErCJDKr8Uprbat7dmttBXBlyiTqKbL2hPqPoHk5lF0hClXFg9D4LdS9A3XvQcFpUQ70ol+cSq/kp23gp617WgqX+rlQ/QzUvgm1r/a0NN2Hvwo8g9z1jO3kv/wOaF6VXJ4Nn0HFQ84JOiFcP6D8dlh/Zk9LoSgJk4jiFW2fRByv9i1yDhJjeZBA2EOvBN+mtDZqeUdB9l5RDvSC1aFGpTfi7WkB2pK+BWTt3tNSdD+OfReAd3DX5JlzMOQcru0PpqcFUJQOkYjiNc8Yc5sxZqIxZoIx5nbE4H4A4Gm/UTMeBvwXp9I7MRpVod9jtMddUfoaiShe5wBNwLPAv4F64KxUCtVjNK8WtxEOxguBGvAXxznI68Zr9FdCS1lKRexxbLNcp9I7aF4p98ShpRSsFbtE2yBpTUvk1xKlHtsANC2TbV3VcxKolyG0piWyDODfIG4P5KRyzkjbJP9Gkd1fDo3fyfPYel3FEGh0l1tKQo7b4F5joCbi2pYGr31911xbPPzVUu7NK+U6QK7ZH6dNaA7u3xFa1olvMJoBL7SsgUBdcjL3VfzVcl+bV7jDtbalC/Ld4C63lPae3kT/xvC6nSzWir/KZOwC/ZVS5pHtf8s68Fd0IJ+NYks9gEnEgWqttfYSYF9r7VRr7WXW2tpukK178Y2F5sVQfBZkOL66PGLvVXoppG8Z/bj6ua6R64o9YUlRt4jbY5RcAkt1omuvoHkVLN0Eyu9z05YMg8oHYVnIPVo1HVYdJDZfkdS9DcsmwpIRUPlY18hVeiksHQfLJkFZ0Bz0px2g9CLI2FqG8ivubut4ePEQeQaXTYaftoWlm0ojDSLfhmvc5VDbtZX7wsr9YcUeUHyem177GizbDFYdIrKkmhVTpdyXbgINc0URWDIOqv8D6ZuH71vwW/lff6Y4VI3nvDmSpZtC2eVSjhnbQtnVUHpZl11Gn2D97+SeLt0UGj6VtMonOpdnSzEsHioTrUCepaqn4x/TXSweAquP6Hw+9R/BsglQ9UzHj117PCyfCmtPctNsAJaMkjYmURYPgeIucCDch2nXVssYswfwIJALjDPGbA+caa39Y6qF61aydoNJEV+mxitf7MPugMLzoh4mgbKDBKpSJV3voWkBBMp7WgoF3BlxkTPjAtWQuStsOtdN81eKMhQrD+i6+xqap7PsyYKxb7mOPuveQzrSIwg0uD1E3iHhvXmB2ujLtknybvhMDPdD03OPhFHPwOL8zlxRYjT9GH5ux/3hxCVt9x35gPxAFM2OzG509i26QSJppA2THsIBhR9G/RvW/goytoHM3dvOHO0orc9TSM9ZoLpzeXYl0XqsO0qsNiPRYwf9Dhq/CU2UP9NBe9KWNR0/fz8ikaHG24FDgA0A1tqvgWmpFKr3kICPLrWjURSlR+h/c5ySxngI9/OtKL2XRBQvrLWRc54HiDVnIobz2vgpitIDGG17XAzqyFrpKySieK0KDjdaY0y6MeZCYEGK5eodGK8Yv8Zr4NJGyn/JhdAS1E+X7ya/VYekzrnhxn+653F+tf9r/7jiP0HJpR07V/08+Gn74Dlek7RV02H5zmLkGsrK/WDFtI4ZWzqs+hksLoIlY2DlQV1jLNtfqX0bVh0qy03fi9HswnRZL72sbZ01XhkKd+rKsi3hpylQcr67T8mfZVvN68nLVflYiK2YD8qD9bR5RbhMxsC637jyLAy6BFg9ndaeC/9GWLmPbAc3L5AYhc6xrXkbqP63bK94EIrPkXRjZJik6tnkr6ujrD5KbD4TUY4CNWKDF/osL8oUO7flu8GKfSWo9k9TYO1vQg4MDu+YNKieBdVBJ8+2RZ6f5bvBir27ZoiqN1H6V6h52S1bT55cY/FZHZtEsWJPWDxSJkRUPACrDpT0lXu59az4D7K88sD4w3NrTxE7xsVD3boMYsO34dYOXV4bnOH2pgXh7xbnHtd9FEWeX8Oqw8LTqv4N604NLofYrjV8JTL/tKP4nKuaJXWv+vkQGZqgbg6YTBnOL/6TpNd/5OYB0LQYlm0lk2ciqXwMfsyT5drXpayc63HsPf0VsGKftscGauQdtNDI/mXXRy2qvkIin0y/B/4JjAbWAG/QX2c1RjLiPmi+xHV4GI0hl0Hdu1KpM6bAiPtp7SVbuY9UWJPV9bI1fAY50yE3aNS48VZoXAA5B8Y/rvyfgA+G3Zj4uZqXABaG3yGN+toToPYN2RYoB2+eu2/du6KM+jeE278lQu1/3WV/mdhs6Fd9dJoWSJzQnP3cGXQEG+hx70LaqPD9Pbmw6Tdga4Mv5r0hbQyMfg5MDvjXgSdfFPrG79161VEa5kPeDBh8PqSNcOMSmnTwhYR3HXY7DD7XXV+xOww6CwpOArxBX1c2ZJaZpdVfkydXPmhsfXje3kL35Vj/mfi5KrpBtuefAg1fQv6xyV1Xoox+MTzuq6cg5q6tjPsA/CGzNAMNsGo/8I6Aomth9c9FufaXiUPnodfDoNNEoQTIPxnqP5FteTOkDOrfh3HvwdoTRSlJJBZtX6FujvhZzP0ZTFgqcTDX/FK2tRRLvUuE+o+Dx6yG+s/Fifaw26Uegdh3eYJt28oDZOaoNz16XjUvSttfdnX4zMiK+8UGbcgFHb7MVpwZgJt8EkywMpHEZEL2/mJbmL1n+DHV/25r89b4FWTuBDl/gYqZbnpzUElq/FLa88avxWaw4RvI+0VQhuCzNuQKmdSxMahMNi2WCR6OfWPzSmmbWlZDekQ454b54TMz/RtgzCsyOa3xW8g7Uup4/ftty8Bf6dqW5R0jk4L6MO2+1ay1ZcCJ3SBL7yNtRPsPsScDcg+D8nukQmbtHLIxxXEc0yfJpACQRjpRTIzGI+65tnLPlbUrVMfxuG0yO54/iIJq60VpSKbHbKCRsZ00eo1fhac79ymSzG3l3+lJDK0/BLelPddJoQxk7gxZU2XVNyb6bmnD2yoDOQdFkb0DsxFNhvuRYzyQsYv7/GZMlq/5VJM3o+PHpE8AQl5Szksuc4qUh8kATFDhbIaMrcLbJU8WpI2OyNQrx3pyOi5Pb8dkyuQRkx7ycu+kE1WnvuQeFmN7GnFtyDyDpLwzthHFI/zgzskGcu+dZ6P1Q8srEyui7p/vujkKlSNzJ+kg8MRro2PI68mX911GRDSMtE0S7GkM5mvSgx0S2XJNNR2MYpGxRZ9XvBKJ1TjBGPOyMabUGFNijHnJGDOhveMUkOJVg09FSYhkPggUBUi993pDQm15tDrcW3yBtaKe/nuaRMZxngLuBo4Krh8HPA3smiqh+iTRXEkYk5qHrvIpaFoEOaFfZxaqHpdzegqkAfDkQ+6h4rPFpEPL2uCutVB+F5hcKDjFHbJwqP0fNC2UnqdAbduvi1Dbq8pH3K55B3+55J9zMDQvlbTcI8VmwNbJ1xtA1j5uL0z9XHfoyDbJfuV3RXyxeyD/hI4PYcajaTk0/yRDdpG0FEP1c/I1N+gMcdrZ8BnknwTebnBPEI2WEhlGqHsfvEGfcTWvQHpXxWS0UPuKDI00LgjadeRAwW+kfuX+TNw8hFLzBqRvJmFx6t6Hgk2TO7XJ6KzwUofK75Khkowp4duqZ0kPXN6xkDYU6oJDTdl7tJ9v3QcyBNMqazYUnCrPTt2H0PBF52Vvg2PrVhZ85gytw8nR9q1+Vp79yB7Qqqfd4Zv0yeJU1ju4c8OugRqofBzwQ+4vwDcq/v4VD0sItpz9O3aehi+kPXB6e+reD3es20qwD6H8bsgMMQ0J1LptiLXSBvoiegcr7hU7wcj6EkYcZcVaCYQObh0uvyt0B/lrXgU1/4VBpyb+kdH4PdS8FufcARnmzNhK3gn5J0HlTNcNRvm98szWvAwN8yArwiFBw3yRycF5dkDausYfpM2peyv8Hdf4lbQ7dXNk3V8aPHaRrNe8FuF2AtcOzLFZa7WZs2J/6R3sOkYuv0t826VPkP0jncc65Q3u+8qkS/325ENm0Bdn9YuQvXfbNquHSUTxMtbax0PWnzDGDGzvZ5HkHCgPVc5BERtS1OO17kQYdLYM+TnkHSXnalwIFffIctomsNlycXwHUiHzjhYXGI0Lxf4g7yjwRtihlF4iSlftHPF+PugPkH+ku73g1+Jw1pMnD0qoF/EhfwXvUMm7/A550ANBj+QlQV9og86Sh7JxAYy4W9LKroOsPSD/ePAMFmWoeSVhDV7Ny+DbRBqSrqLsMnkxbRHlPtW8BqWXS4PTslbKoupp8E2E3EO6ToaOUPsGlFwABb8Tm4i0UdA4Qxrd3CNg0O8TyMQLQy6XuImR5M2QoMOBeqh8GKqeFFuMnP1h3SniOyn/mPBjVk+XBn/kI2Inkj2zbb7tMezWzsdwNJlQeInU7Ywp0uA6ZB8kL+yNd0gjnH8crDoYsDA5AX/QpZeLUpcZfOYqH4Tcn4sCV/Y3sbHqKkwGFF7sDnsV3SCTfNI3l2vL3LHtMTk/k+H5qsfkBTv0WkkfdJbYzzQulJdV6ZXyPJqMzile9fOg+I9igkAaDG6n3q0/Xf6jPWfxWH2UPHtbBF/WpZdImUcOdw29Ekr+4rZtIM9t1eNiy5g7Q56d5iVSJiDDg74x8p8xOLy+tCFOj5ej5HiHw5BLRcF3ZPAOAV/QmXH1LJnMkr2nDEkmQsX98sFaGDEhashf5R5m7gCrfyYKWO0rkLUXrP+t2FhioOwaUUrK75J2POfgcB9zG/4BjV9I25t7hLzHfBPFd2X1v2WSStUTolhlBT9QfBMAjzh0zT1CPkoztnCvGWDjTZA51X1eIOhvbQepv7bJVT5zfxaMVrFQ5PLki41m1eOIuY5f7tnQ6+QZ9w6B5hA/YMXnyNB7wxdyz32bwcTFsm3NUTD83vbrZzeTiOL1jjHmEuAZ5G4dC7xqjCkEsNZuTKF8fYPMnWDETlE2eEiZjdeI/wtfzzlAfgCVD4Q/XA4Fp8Hw2931qseJyaCzRPEyPhhxV/i2vKOCil4cmhaJQXD+cfLV5ODbTPIrvzs83fig8ML4+a7qgHfvRAm089L1DpOHOVAX9BXUC8j7Vfg9cZTXRDEGiq6Lvi17n3AbvcHnw4YEJmK0GvJ6IHP7jskDUPjnjh8TifHEnjSSNVV+a45z0zoSgMNkwNC/uh9X1SGev026vNTKb49+bEcxHhh2k7teeH7sfR2y95Bf+b+CitcVkj74THef2rdkBlz6lq5BdWfI2kdeuKkkUAGE9LCbTDHwbmMfuD+Mjwgh3LRc2rjCC6DwT1B8vnwMOs9OZLsWFxNiWxUFT57YTaVvBptXuOlZu0NtB22YIs/ryN+aZKDomtCTh09C8g6H4bfJcsNn8p93lHu9Tk8vyGznIX8NTmqJxB/sXQy2e/lBU29vgUzsqv0vDD5Hnom8EK/66RPk4zDn0Ag5Y5C1e/hH14i7pber6nEZ3fBvCL9XTcvCRz08OVB0i3QwBKrajj7ZWL3EPUciipfzWXRmRPppyJtd7b1ikqBdgKIofZyB5vpkgLVrJsm23Hh6oY0X8ZXIfkfvczuayKxGDcyXLMYTnJ0XMiPG+OTrxJm5FLZ/unzRdWXoISewcMztUQJeOzZctpnkXyg2fDnqeRrdgKsJfZX4o/tFCzTEDhfiyQ/aBwRk9pdJl0bHKePWIM7R5AsJPGzrgWCUgkC1HBfq18dkJBayxGS0M6MoiL9aZG6zf4BucRQZqAy/HmdGYKQvI+e+BSqj3+PeSKCuY0Gl/VVt66dtkuv1+2Rbqvz1dTWtdpTBZ8+kuTZQgTq5FpPhKgzx6qpjG2rrJS9PXvC44PPlyZUelbCwVPXyHLaRyx89CLRzbGs70UMKrm2W3jd/OuARFzqBhuBzH68uGbesbUhbE629iXreBMMgtT6Hkb09LW179Y2R3l5/ZXzfZK359eBHRaznKlDXtk60lkF1ePlav1sOnvy2Ns09QEzFyxizM7DKWrs+uP5r4JfACuAqHWJMAN8ECUAMcsNB7IR8E2XKsZMmG4A0GWJYf3b0xqk1303bOe94Geozae75oe2Uc9+m4iivDWlBP1Cd6KJNC7oRSBshL6jSoI2CM5ToHSGG69WO+wIDQy6Jn2egHmpmQ97Pw9N/2gb8xbSZpBuoF6PwygfEQDxzR9jkfRkmWntCUCkLNlTR4hgCZB8odiE1s8VwOW0srA3aOLXe08agz7HsdvyOBWTocuLS+NdZ954EfvbkgnckTAyJAdjwTceUhmTwjZf4gQDD75Z6snIvWV/3G/Hx5RirOkN3tW9K8HRfL/9O840ROxvH5qg9al4ROyNPrtguOniLxAmtJwvwiKF9xT0pEblDpI+XehgN71B5+fsmiv+lpePkZbTZKnk+Fg8VRctkB20Hv2zfJittNJRdJbZXQ6+HIReKndm638gQ2bB/iA8shyVjYPMNbfMpOR/K/y+iTUTaS4DFg4LbvIkbSjv+BR0n1+mTknd1kzYClm0t9ztQC+PegZXTZJhuw7WxfT16R0DdO+Hty8q92l5nPEbcH3+7bzzUfyjLK/YIl8U3BiofgsEhsYY9hWL3t3QcYCT+YjTSRsqQYet6iAsT31ixjQt9JkKPg7a+BDuC4z8tcyrU/a/tNn9xSJmmiWF+2khoWSF1fOk4t2036bD+91IOw+8LH37vKay1UX/AF0BhcHkasBZRvK4F/hPruER+wHRgEbAEuKS9/XfaaSfbZ1l9tLULcNcX5Vu7ZKK19fPC9ws0WrvAZ+2GW61df37nz9u4yNqlm1u7ZLy1jUuTy2MB1i7dqvOydBUb7rB2/blt038cbm3zurbp6/9s7aqfW7t8b2vr51u7bAdJL7/f2rW/TV6OFQeE39ON/yfrdZ/FP655vbU/Dms//+pXrF15mLVNa6xdPDJ8W8mV8uspFo+ytmm1u77ycGurZvecPMmyAGsXbyLXE4/Kp61dfWy3iNQjOO1Dc6m1Pw6x1l9t7cIca5dsFl7HI6mZY+2Kfd310Hq58S5rF3isXXemrC+bYm39F9b6G6xdmB49v7WnW1v+QNv05hKRo/z+ZK4uNaw8xNrq10Sutb9t+4wqvYMFWLswW+rO6mPlWS/7R3dKEFOnifd57rVur9axwP3W2lnALGPMV8kqesYYL+Ke4iBgNfC5MWa2tfaH+EcqiqIovYeeH7JRlLgYb680A4g3TctrTOu4yQFAqDOnzsRx2QVYYq1dZq1tQmZLJuHuWVEURVEUJRZBdxS9jHgK1NPAe8aYMqAe+ADAGLMZ0Bkr2tFAaLyZ1fRnZ6yRzkUDVUEj0kid1wAt4jclVtiKDuEJTgW2Uc7VkWxi2Iv0CB6xtYoMwBooJ+o1Go84Zk3fUrY3L4VVR4gdQGd8RhlvW7mc88XFI/5qVh0Rfzd/sUwJNx5oWRe+f9OiGFO/uwuP+PNyQvM0fAamd/nISRjjFee48e5Hy5rUu0zoUTzi94k0WSZoeN28RDbHKht/WbhzY+MRX3wN88UHn0kTP3irjgg6UfYEjbpboufZ+LXrJyqU1meql7hyATHiXh+0i6p8MEq4JqXX4MkTJ7KBjYBP7Lzq3oVRj3etI+4OYmycaaXGmN2AkcCb1orTG2PM5kCutTYpV83GmGOAQ6y1vw2unwzsYq09J2K/M4AzgquTEZuwVDIUKEvxOforWnbJoeWWHFpuyaHllhxabskzkMuuzFo7PdqGuEOG1tq5UdJ+jLZvB1gNjA1ZH4MY7kee536gnekcXYcxZp61dmp3na8/oWWXHFpuyaHllhxabsmh5ZY8WnbR6Yn+28+BScaY8caYdCT24+wekENRFEVRFKVb6YyRfFJYa1uCsR7fQCzfHrbWft/dciiKoiiKonQ33a54AVhr/wv8t90du5duG9bsh2jZJYeWW3JouSWHlltyaLklj5ZdFOIa1yuKoiiKoihdRy+ao6soiqIoitK/UcULMMZMN8YsMsYsMca0EzBw4GGMWW6M+dYY85UxZl4wrdAY85YxZnHwf3DI/pcGy3KRMeaQnpO8ezHGPGyMKTHGfBeS1uFyMsbsFCzvJcaYO43pBVFdU0iMcrvKGLMmWOe+MsYcFrJNyw0wxow1xrxjjFlgjPneGHNeMF3rXBzilJvWuXYwxmQaYz4zxnwdLLurg+la5zpCvHhCA+GHGPgvBSYA6cDXwFY9LVdv+gHLgaERaX8nGGcTuAS4Obi8VbAMM4DxwbL19vQ1dFM5TQN2BL7rTDkBnwG7I151XwMO7elr64Fyuwq4MMq+Wm5uWYwEdgwu5wE/BstH61xy5aZ1rv2yM4gfTwAf8Cmwm9a5jv20x0tDGCXLDODR4PKjwJEh6c9YaxuttT8hgdB36X7xuh9r7fvAxojkDpWTMWYkkG+t/cRK6/RYyDH9khjlFgsttyDW2nU26MjaWlsNLEAig2idi0OccouFllsQK9QEV33Bn0XrXIdQxSt6CCONARGOBd40xsw3ElEAYLi1dh1IQwYMC6ZreYbT0XIaHVyOTB+InG2M+SY4FOkMXWi5RcEYsymwA9IDoXUuQSLKDbTOtYsxxmuM+QooAd6y1mqd6yCqeEk3ZyQ61TOcPa21OwKHAmcZY6bF2VfLMzFilZOWn3AvMBGYAqwDbg2ma7lFYIzJBWYBf7LWVsXbNUragC27KOWmdS4BrLV+a+0UJOrMLsaYbeLsrmUXBVW8EgxhNJCx1q4N/pcALyBDh8XB7mKC/yXB3bU8w+loOa0OLkemDyistcXBBj4APIA7XK3lFoIxxocoD09aa58PJmuda4do5aZ1rmNYayuAd4HpaJ3rEKp4aQijuBhjcowxec4ycDDwHVJGpwR3OwV4Kbg8GzjOGJNhjBkPTEKMKAcqHSqnYDd9tTFmt+Asn1+HHDNgcBrxIEchdQ603FoJXudDwAJr7W0hm7TOxSFWuWmdax9jTJExZlBwOQs4EFiI1rmO0dPW/b3hBxyGzGxZClze0/L0ph8y2/Pr4O97p3yAIcAcYHHwvzDkmMuDZbmIATRTBXgaGaJoRr7oTk+mnICpSKO/FLiLoKPj/vqLUW6PA98C3yCN90gttzblthcyPPMN8FXwd5jWuaTLTetc+2W3HfBlsIy+A/4WTNc614Gfeq5XFEVRFEXpJnSoUVEURVEUpZtQxUtRFEVRFKWbUMVLURRFURSlm1DFS1EURVEUpZtQxUtRFEVRFKWbUMVLURRFURSlm1DFS1EURVEUpZtQxUtRFEVRFKWb+H88FrBbab+bqwAAAABJRU5ErkJggg==\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"data_slice = data.iloc[:]\n",
"\n",
"fig, axes = plt.subplots(3, 1, figsize=(10, 3), gridspec_kw=dict(hspace=0.2))\n",
"\n",
"ax = axes[0]\n",
"x = data_slice[\"time\"].values\n",
"y = data_slice[1].values\n",
"ax.step(x, y, lw=1, color=\"red\", zorder=2)\n",
"ax.set_ylim(0, 11)\n",
"ax.set_ylabel(\"Species 1\")\n",
"ax.set_xticklabels([])\n",
"ax.spines[\"top\"].set_visible(False)\n",
"ax.spines[\"right\"].set_visible(False)\n",
"ax.spines[\"bottom\"].set_linestyle((0, (3, 6)))\n",
"\n",
"ax = axes[1]\n",
"x = data_slice[\"time\"].values\n",
"y = data_slice[2].values\n",
"ax.step(x, y, lw=1, color=\"royalblue\")\n",
"ax.set_ylim(0, 11)\n",
"ax.set_ylabel(\"Species 2\")\n",
"ax.set_xticklabels([])\n",
"ax.spines[\"top\"].set_visible(False)\n",
"ax.spines[\"right\"].set_visible(False)\n",
"ax.spines[\"bottom\"].set_linestyle((0, (3, 6)))\n",
"\n",
"ax = axes[2]\n",
"x = data_slice[\"time\"].values\n",
"y = data_slice[3].values\n",
"ax.step(x, y, lw=1, color=\"gold\")\n",
"ax.set_ylim(0, 11)\n",
"ax.set_ylabel(\"Species 3\")\n",
"ax.spines[\"top\"].set_visible(False)\n",
"ax.spines[\"right\"].set_visible(False)\n",
"\n",
"pass"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAUcAAAE9CAYAAACY8KDMAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABXqUlEQVR4nO29d3yb5bn//761JUvyzN6bDCALSAmrh9KmZZSWLg7zS0pK6YHCKQfCLIVSWlrK4VdKIaWUFWgPZYS9QsMqCSQhi4QMyPa2bMmytnT//rCfp5YtJ7b02JLt+/16+ZXokXTpeiQ9H93zcwkpJQqFQqFIx5TvBBQKhaIQUeKoUCgUGVDiqFAoFBlQ4qhQKBQZUOKoUCgUGVDiqFAoFBmw5DuB7rBo0SL52muv5TsNhUIx8BBd3dEvWo719fX5TkGhUAwy+oU4KhQKRV+jxFGhUCgyoMRRoVAoMqDEUaFQKDKgxFGhUCgyoMRRoVAoMqDEUaFQKDKgxFGhUCgyoMRRoVAoMqDEUaFQKDKgxDGPNDQ05DsFhULRBUoc84jf7893CgqFoguUOOaRWCyW7xQUCkUXKHHMI9FoNN8pKBSKLlDimCcSiQSqLK5CUbgoccwTjY2NlJaW5jsNhULRBUoc80Q0GsVut+c7DYVC0QVKHPNEMBjE7XbnOw2FQtEFShzzRCgUoqioKN9pKBSKLlDimCeklAjRZW0fhUKRZ5Q4KhQKRQaUOCoUCkUGlDgqFApFBpQ4KhQKRQaUOOYJtTtGoShslDgqFApFBnpNHIUQDwshaoUQWzocv0IIsV0I8akQ4q7eev1CRy3jUSgKm95sOT4CLGp/QAjxZeCbwFFSypnA73rx9QuWaDSq7MoUigLH0luBpZTvCiHGdzj8Y+DXUspo22Nqe+v1e4Pbb7+dr33tawwdOhSbzdbt5wkhEEJgs9mw2Wxs27aNjz76CCklxx9/fC9mrFAosqXXxLELpgInCiHuACLANVLKj/s4h6z55JNPeOmll7j77rs54YQTevTcZDJJPB4nFoths9m46qqrWLhwIQ8//DBTp07tpYwVCkW29PWEjAUoBRYA/wP8n+hi8E0IsUQIsVYIsbaurq4vc+ySa665hgsvvJDPPvusx881m804HA68Xi9ut5s777yTF198UQmjQlGg9HXL8QDwrGxdx/KRECIFVACd1E9KuQxYBjB//vyCWPdSV1fHt771LUaOHJlTnObmZi644ALl56hQFDB93XJ8HvgPACHEVMAG1PdxDlmTTCYNiTN79mw1W61QFDi91nIUQjwFnAJUCCEOAD8HHgYeblveEwMukv1oNbRR4ghqKY9CUej05mz1uV3cdX5vvWZvk0qlDBO1fvSboFAMStQOmR6QSqXynYJCoegjlDj2ANWtVigGD0oce0BLSwsOhyPfaSgUij5AiWMP8Pl8hi2/UWOOCkVho8Sxm0gpDe1WKxSKwkaJYzeJRqNYrVbD4qkxR4WisFHi2E2am5sNLaWqutUKRWGjxLGbxONxLJa+3m2pUCjyhRLHPKG61QpFYaPEsZukUimcTqcyqVUoBglKHLtJXV0dkyZNorm52ZB4UkoCgQB+v9+QeAqFwliUOHaTlpYWYrEYu3btyjlWJBLhxRdf5Be/+AU+ny+nWIlEggMHDhAKhXLOS6FQ/Bs1w9BNfve737FixQq++93v8re//Y1oNEokEiESiegzz12NI3acmV61ahW33nor1dXVjB8/nv/8z/+kvLw8q7yCwSAvvvgiFouFOXPm9NhrUstZK+Hgcrkwm81Z5aJQDCSUOHaTP/3pT6xYsYLdu3ezY8cOhg0bhsPhoLi4GJPp3w3wSCRCU1NTJ5OK9sJZWlpKdXU1AOFwGLfbnXVeX3zxBR988AH/+te/ePXVV7My4k2lUnoJh4MHDzJ27Nis81EoBgpKHLvJ+++/zw9/+EM2bNhAc3MzXq+XcDjc6XF2u52Kiooul/0EAgG2bNmC1Wpl0qRJXHDBBTQ0NACtAtq+FdqxxZnp/uHDh3POOedw7rnnMm3atKzOzWQyYbfbsdvtBAKBrGIoFAMN0R8WI8+fP1+uXbs2rzlUVVWRTCZxu90kk8msu8HQ2s3+v//7PxYsWMCYMWPSWp7Zsm7dOubNm5dznIMHDzJq1Kic4ygU/YQu19SpCZluMmzYMGpqagiFQpSVleUUSwiBy+XCZrMZIoxaTIVCYRxKHLuJyWRCSomU0hAhMrrF3h96AApFf0KJYw8YDAKkWqAKRStKHHuIUeJhMpmUBZpCUcAocewBRraqiouLDZ0ZttlsRKNRw+IpFIMdJY55wul0ZlwKlC1lZWU0NjYaFk+hGOwocewBNpvNMOMJo8f2TCaTqo6oUBiIEscekEwmDVt6YzabSSQShsRSKBTGo8SxB1itVuLxuCGxvF6voWOOg2EmXaHoS5Q49gAjx/XsdjuRSMSQWAqFwniUOPYAI8cJjR5zVOsTFQpjUeKYR4zsCqtutUJhLL0mjkKIh4UQtUKILRnuu0YIIYUQFb31+r2BkQKkWnoKRWHTmy3HR4BFHQ8KIcYApwH7evG1ew23200wGMw5jpFCG4/HM1qcZUsymVS7dxSDnl7zc5RSviuEGJ/hrnuAa4EVvfXavcHbb7/NypUrWbhwIYlEglmzZmUda9OmTXzwwQesWLGC+fPnp91ntVopLS3tURnY3/72t7z44ossW7YsJ7ux6upqVq9ezfr167n88suzjqMhhMDpdFJSUpJzLIWir+lTs1shxFnAQSnlxv7UrUylUmzdupVf/epXAGzdujUnd54//OEPPPTQQwD861//4tvf/ra+fjIWi+Hz+TKugez4eslkkhtuuIHHH38cgAsvvJA///nPaY+x2+2UlJTgcDgwm80ZZ8i1Vmc0GuXXv/410WiUa665JmdRe+edd3C73ZSVlTFy5EjsdntO8RSKvqTPxFEI4QJuBL7azccvAZYAebft37NnD88//7x++8knn2Tx4sWMHz8+q3gXXnghtbW1vPDCCxw4cCBtYbnNZmPo0KHdilNZWcnOnTv12xs2bGDdunXMnTtXF1K/38/evXuJRCKkUqk0D0khBA6Hg5KSEoQQVFVVceSRR/Lwww/zwAMPsHTp0qzOT2PXrl2MGjWKsrIybDZbTrEUir6mL1uOk4AJgNZqHA2sF0IcK6Ws7vhgKeUyYBm0OoH3YZ6dmDhxIn/60594+OGHmTBhAuXl5VkLI8C4ceM48cQTWbx4MYsWdRqW7TbDhw/nyiuvJB6Pc8sttzBixAiEEMybN6/brdpQKERTUxMAQ4YMYdOmTcyYMYNrr70267w0xo4dy7Bhw5gwYULOsRSKvqbPxFFKuRnQm0RCiD3AfCllfV/lkAtjx47llFNO4aijjup2y64rxowZwymnnEJlZWVOLSqTycQPfvADjjnmGCZPnpxVDJfLhcvl0m//v//3/xgyZIgh2yRnzZrFiBEjco6jUOSD3lzK8xTwITBNCHFACLG4t16rL9AKZ0kpcy5davRi8vbiZgSjRo1Ss9WKQU9vzlafe5j7x/fWa/cXCnFSSgiB2+2mpaUFr9eb73QUiryhdsjkCZfLRSgUyncanUgkEng8noLMTaHoS5Q49hCjWnslJSUFaU7r8/koLy9X2xEVgx4ljj3EKNEw0jjXSBKJBFarVYmjYtCjxDFPFOJ4o0Kh+DdKHHuAEMKwWdxoNIrb7S7I1qNCoVDi2COGDBlCfb0xyzKllJSWlhbkuKNCoVDi2CPMZrNhLUchBGazuWCLYhnV7Vdjl4r+ihLHPCKlLEjxMMr+rFCXKykU3UGJYx4pKiqipaUl32n0Gh6Ph+bm5nynoVBkhRLHPKGNOWqmD7lg9HZEo1qzqpa2oj+jxLEHVFVV5TuFjKhaNAqF8Shx7AGvvPKK4eUN/H6/IbFqa2sNXRak1mEqBjt96gTe12zbto2VK1cyb948xo4dm/UFb7fb+fWvf83vfvc7brnlFj7//HMmTJiQdbxPPvmEU045hUAgwLXXXstXvvKVrOJA6zmuWLGCBx98kPvvv5+vf/3rWccKhUJs3ryZESNG0NjYmPGHwGw26zZnh3Inqqurw2azsX79etxutzKxUPQ7RH/oRs2fP1+uXbu2x8877bTT2LlzJ9XV1fzsZz/j5JNPpry8XL+//blrQheNRvH5fGn3rVu3jl/84hdpse+++26KiopIJpOEw2HC4bBe2qCoqEgXDy2uNpYXDoe5+uqr9Ther5empqY0odVKJXRnvO7222/ngQceAGDhwoWsWrWqy8ce7rM+cOAAF198Mdu3b+f000/nL3/5S6fHJBIJwuEwoVAobVlTxx+KRCLBW2+9xbXXXsvUqVN56623cDqdhz0fhaKP6bKFM2BbjqlUCrfbzbRp0xg3bhzf+MY3cLvdJBIJXSTMZjPl5eVpxaxsNhtlZWVpZq8TJkxIE8cRI0Zw4oknIoTAZDLhdDpxuVxYrVZSqRShUIiWlpZOayKFEJ0WkQcCAVasWIHJZNIFxmazUVpaqrfM2teraS9wNTU1aZUQ58+ff9hF5Ydq7e7btw+fz0dNTQ1HH310p/ullLS0tNDS0oIQolMsLTchBF988QU/+MEPuOSSS7jnnnuUMCr6HQNWHE0mE1dffTWhUIht27bhdDqZMWMGVqtVf0w8HqexsTGtmFU8Hqe6Or1qQ1lZGT/96U+59957ufjii7nwwguZMmVKjwtQSSmprKzksssu44MPPmDu3LkMHz6cs88+O6tzlFIyY8YMvvKVr7BmzRquv/56hgwZklUsaB0+uOGGG6irq9OdytuLshCC4uJiRo4cedhYiUSCF154gbq6OlV9UNEvGdDd6srKSvbu3YvP52Pu3Lk5Wfbv3LmT9evXc9JJJzF8+HCklFmXEqisrGTdunWMHj2aYDDIiSeemHVe0WiUN998E5vNxle/2q3aZV2SSqUwmUz8/e9/58QTT+yWCHZFZWUlLpdLCaOi0OmyKzXgZ6sTiQSNjY1UVFTkFKeoqIhJkyZRXl5OQ0NDTjVWTCYTdrudQCCQc3czkUgwbNiwHtW5PlRe0NpSNkLUlDAq+jMDXhwbGhooLi5O607ngs1mIx6P5xSjqKgIi8WC3+/PeclMMpkklUoZdn4A4XA457o0hepXqVB0lwEvjslk0pBWlZG4XC7C4bAhsRKJBCUlJWkTM7lg1H7vgb41UjHwGfDiGAqFDK/Olytms9mwxeTJZBKHw5Fza1ajsbGRsrKynOM4HA7DfgAUinww4MUxGo0a1uXMtJwmW4wSx0QikbaeMlcCgYAhC7bVDhtFf2fAi2OhGsoaJY65zJpnwihxVCj6OwNeHDMtVs4Go+23ysvL8fl8hsUzSmxjsRh2u92QWApFf2bAiyMY08Vzu91p4pirGDmdTiKRSK5pGU4kElHiqFAwCMSxqanJ8DG0iooKGhoaco5pFO23F+aKz+czZEIG1Lijon8z4MWxpaXF8NlqI9Y6ahMpRqDVmjYKJWoKxSAQR6Mu9Hg8buh6Sb/fT3FxsSGxYrEYNpst5ziBQACPx2NARq30h62pCkVXDGhx1LqbRixIbmxsNHQ7XCQSwel0GlLNMB6P43A40gw0siEYDBoqjgpFf6bXxFEI8bAQolYIsaXdsd8KIT4TQmwSQjwnhCjprdcH9L3LxcXFOTtuJxIJgsEgu3btMiS3VCqFEKKTd2Q2xGIxhg8fXlDjoO1R2wgV/ZHebDk+AizqcOxNYJaU8ihgB3B9L74+n3zyCWvWrOHTTz/F7/dTWVmZVZxgMMh7773H//7v/7Jly5bDP+EwJBIJVq5cyTPPPMOyZcvYvXt3TvGqq6upq6tj+/btJBIJUqlUVoIbCoX45z//yfbt21m5ciUHDx6ksrKSmpqaHsc7ePAgkUiEd99917Ba3wpFX9Jrm46llO8KIcZ3OPZGu5urge/01usD7N69m2XLlrFs2TIuuugifvnLX3ZZJEurBtjRJWfHjh2sX7+e9957j1deeYVXXnmFJUuWsGTJkm4LRiQSwe/36yLxxBNP8Ic//AGAoUOHsmfPHjweDyUlJVitVqSUen2ZqqqqtMkfKSWxWIxQKEQ4HObNN9/UY23fvp3GxkZSqZTuIt4+R83lPBQK0djYSF1dXZr576OPPsrq1au59tprue222xg3bhxFRUXEYjHWr1/fKV5paWnaZFd7s9tkMsk555zDhg0beP7551m4cGHOzkgKRV+ST0eGS4C/9+YL7NixQ/9/PB5n2LBhXc7qhkIhdu/enbb2ULvYH3zwwbTyA8uWLWPKlCkZx+faTwBpz7dYLAghiMVifPjhhzz22GP6Y2pra3nsscew2+34/X7i8ThSSmw2G16vlxEjRuB2uykpKUFKqYuf0+mkurpaF0aAadOmcfPNNxOPx/UuttPpxOl0YrfbcTgcFBcXk0wmqaurw+12E4lEiMfjRCIRVq9erce65ZZbOO+885BSYjab07wwtXNsbGxk//79BAIB/ZjZbGbo0KHU1tayYcMGAM4++2zq6uoyvu8KRaGSF3EUQtwIJIDlh3jMEmAJwNixY7N6nfYtroqKCqqrq7FYLBlbfE6nk+nTp2ec3b7lllvSxPG73/0uFRUVjBw5stOawPZ1n7X/m81miouLsdvtGetUX3jhhdjt9rRtgO139oTDYb3OzIwZM/RF2qFQiJdeeolPP/007XlFRUWMHTuWhoYGampq2L59uz7uZzKZcLlcDBs2TN8q2JW92OjRo/VZ8EgkQkNDQ9r5lZSUMGrUqLRZ90QiQV1dHR6Ph29961s899xzjBs3joMHD6qWo6Jf0atO4G3d6peklLPaHbsIuAw4VUoZ6k6cbJ3AX375ZTZu3MgzzzzDunXrevx8jWQyyYMPPshNN91EY2Mjb7zxBqeddlpWsWKxGE8//TQ//OEP+fGPf8zChQs54ogjKC0tzdjq7Ci2HVmyZAkvvfQSAF//+te57LLL9PsqKiooLi6mtLRUPyaEwGw2U1RUhNPpTBPkX/3qV+zcuZPZs2fjdrtZtGiR/pp2u52ysrIeLY3au3cvK1eu5JRTTmHixIndfp5C0YcURoEtIcQi4Drg5O4KYy7MmTOHuXPncsQRR+QUx2w2c9xxx3Hqqady1VVXMW7cuKxj2Ww2vvzlL/P3v/8dq9XK0UcfnVM5guXLl/P8888TCoU466yzcor1ta99jRtuuIE9e/bg9/sZNWpU1rEArFYrRx11lBJGRb+kN5fyPAV8CEwTQhwQQiwG7gM8wJtCiA1CiAd66/U1/H5/WsspW6SUTJ8+Xa/HnCtCiJyKYWl4vV7DfBMtFgvxeJxt27Yxffp0Q2IqFP2V3pytPjfD4c6FkHuZYDCI2+02JJbJZKK2ttbQsbNC2qpXUVFBfX098XjckB8AhaI/M6B3yJhMJn1JS65oNaoTiYQh/omZltpkizYTnqubTvvxTYVisDOgxbGkpMRQo9uSkpKMs83ZYGSLMZlM4vf7c97e2NGWTaEYzHQpjkKIMUKIvwkh3hNC3CCEsLa77/k+yS5HHA6HoZ6JdrudWCxmWGs0Ho8b4qbT1NSE2+3O2eXH6/WmrVk0AtUSVfRXDtVyfBhYBVwBjADeEUKUt92X/XRtH9Iblv9G7NOGVtEworUHxpaCMFIYTSaTEkdFv+VQ4jhESvmAlHKDlPIK4H7gXSHEJKBffOODwSAmk4mioqKcY2kOP06nk1Ao91VILS0tOBwOQ1qOFovFsP3LoVAIh8NhSCwjx3wVir7mULPVViGEQ0oZAZBSPiGEqAZeB3JXmz6iqamJKVOm5BxHW+ZiVKW/ZDJpmNmtkbW5W1paDPOZVOKo6M8cquX4EHBc+wNSyreA7wK5W9P0M7RlLmaz2ZAL3uv1FuTkjt1uN6wGthJHRX+mS3GUUt4jpXwnw/FPpJTZ7Z3rxwghqK+vZ8qUKTn7Jgoh8Hg8BAIBg7IzDiMdz4UQShwV/ZYBu5RHs/w3sti93W7H7Xbn3LIqKirCZDLl7E6uYeSkh5HdfTUho+jPDFhxbG5u1i3FjBDIUCiE1+slkUjkHM/tdhsmjBrZGtx2xMjxS9WtVvRnBqw49sa2PG1Rea4iZDKZdN9Go9ZhlpWV4fP5coohpaSystKwnNqLY7bvWTwe5+DBg+zbt8+wsVCFojsctokghPgp8FegmdZJmjnA0g6u3gVJXV0dq1atYsKECWlmrdng8/l48803eeutt7j44osZNWoUUsqsrLzWrFnDk08+SUNDAxMmTGDWrFnY7XZKS0t7vDXxwIED/OUvf+Hjjz9m8+bN3HbbbXg8nqxmnHfs2MH1119PMpnktttuY+LEifr75na7u4zZXvi0bZbV1dWUlZXxzDPPMHz4cDweT5e+nOFwGJ/Pl/E9tFqtjBw5ksbGRrZs2cKcOXN6fF4KRTYc1s9RCLFRSnm0EOJrwE+Am4G/Sinn9kWC0HM/x1QqxWOPPcaBAwe4+eabgVZvwdGjR2e1L/qLL77g8ssv5/XXXwdgxowZPPLII0Br6YFEIsGYMWNwuVydWkjafmWTyURZWRn79u1LW1r07W9/m4svvlgvpeD1eikvL8fj8WA2mzsJj5RSj+f3+zn11FPTXm/58n/7B3s8Hn19ZiaKi4v1MgexWIwrrrhCL4cAsHLlSn0RfUtLC36/v1Os0aNH4/F4GDZsGG63Wy+RsHnz5jQhu/feexk/frxewkH7i8ViehE0KSWpVIpkMonD4cDlcuku6ldddRXvv/8+L7zwAmeeeeahPzCFovvk5OeoPfkbtIriRlFIVjIZ2L9/P06nUxdGgLPOOou//vWvaYultVaOdkGaTCb94tdERQjB008/rQsjwNatWxFC6N6JgUCAL774Ql8cHolECAQCJJNJEokEgUCAcDhMfX091dXVabk+++yzfOtb32L69Okkk0kCgQD19fXs2LGD6urqtDygVcRMJhNmszljJcQ333yT73znO3rXWHtee2HVbldVVeFwOPQfjPbCCHDPPfdw2223UVFR0WmSxul04nA4qKuro7a2lvXr1+tDDmazmbfeeivt8Vu3bmXChAlYrVZMJpP+p7n/aO+9dlyrd6MNPbz//vv657hq1SpOPvnkTB+9QmEY3RHHdUKIN4AJwPVCCA9Q0KPsqVSKBQsWcNddd3HttdcCcOeddzJs2LBOj9VEUdtlkkn3Tz31VB544AF9XeJFF13E/v3704ROEwWA8vJyJk6cyNChQ/W6LQ6Hg2QyycqVK3n++efT4peXlxMOh4nFYiQSCVwuFx6Ph3nz5unC1b7kgtbCAnjyySfTYl177bVMnTqVWCx22C2FWp1rbZfOcccdx5o1a/T7FyxYwPvvv091dbXeQk4kEkSjUcLhMOFwGLfbzZAhQ5g+fbreBW9ubu60xfKHP/yh/j5p59P+h6g9QgisVitut5tEIsHevXv1+375y18ye/bsQ56XQmEE3elWm4DZwBdSyqa2/dWjpJSb+iA/IPsyCQ0NDaxfv15vjXz5y1/O6vVTqRQmk4kJEyZw9tlnc88993TrOT6fT59E0N7n/fv3s3HjRn70ox8Brd3q++67D5vNpv9pravusHv3biZOnMgZZ5zBSSedxDXXXNPt8c9UKkVDQ4Oe42effaZ307///e/z+9//HiEEFosFq9Wq/6vll0wmaWxspLm5mf379+t1brRzfeyxxygqKuI//uM/OPPMMzvV2+kJH3/8MT6fjxNPPDGt4qFCkSM5daslMAM4A7iN1q2Dxmy+7WXKy8uZOXMmQE7lAzShuuOOOxgzZky3n5PJFHfkyJEcd9xx+P1+5s+fz9SpUxk+fHjWs+sTJkzgvvvuo7y8HJfL1aM4JpMpzY182LBh3HTTTYwZM4bTTz/9sO+Z2WymoqKCiooKJkyY0On+kSNHEgwGMxYi6ynHHHNMTs9XKHpKd5on9wNfAjRn72bgj72WUQHjdrs71bXOFrvdzpAhQ9K6l9liNpuZPHlyznHq6+uZN2+eoW7n4XDYMCMLhaIv6U7L8Tgp5VwhxCcAUspGIcSg9NCPRCKGXejaDh4jiMVihpQ10HbHGOEq3j6mUYvKFYq+pDstx7gQwkybTZkQYggFPiHTnlQqZdiCcKP8F6HVgzHXRdsa0WgUu91ekFv1Cnxhg0LRJd0Rx/8PeA4YKoS4A3gf+FWvZmUgPp+P8vLywz+wG2gTM0ZgtVr12eJc0RzKjUI5gSsU3ehWSymXCyHWAafSOrNztpRyW69nZhBGdTl7g6amJn3CKBc0hx8jRM3v9xvm56hQ9Ge6FEchhFdKGRBClAG1wFPt7iuTUhrTJ+xnGNmqCofDhixLaW/Emysejyfj4vJssFqtaj+0ot9yqJbjk7Qu31lHelkE0XZ7Yi/mZQhGjjdqGNVNNKq4FrTOVldXVzNxYu4fiZE2Y2VlZWSzPlWhKAS6FEcp5Rlt/3ZewNZPMHICBYwdPwsGg7jdbkNiaS4/Ro2HOp1OwuFwzsuWzGazYbVtFIq+5rBXkxDiW0KI4na3S4QQZ/dqVgZhxAWu0XGPc65EIhHDctO61bnklkgk8Pl8DBkyhClTprBz505DctNQEzOK/kZ3mho/l1LqG2WllE3Az3stI4Mx6qL0+XyUlpYaEstohgwZQl1dXU7i2NDQQHFxMVarldGjR3PgwAHD8nO5XITDYcPiKRR9QXfEMdNjBt2q3mg0auhODyklNpuNaDSacyyPx0Nzc3NOPwSpVEqf0LFYLIYtMwLjan0rFH1Jd8RxrRDi90KISUKIiUKIe2idpCl4DuVl2BOSyaQuPEZ2D8vLyw1ZCG6UC3j7f43A6XQSiUTSzCgUiv5Cd8TxCiAG/B34PyBMq+ntIRFCPCyEqBVCbGl3rEwI8aYQYmfbv73aT62vr0dKSW1tbU7W/1JKPvzwQ0Mv8FQqxf79+w0Z27NarWzfvj2nGLW1tQwdOhQhBIlEwpAlOG63m2AwmHMchSIfHFYcpZQtUsqlwClSyvlSyhuklN2pDvUIsKjDsaXASinlFGBl2+1eY+vWrfzwhz/kgQce4ODBg1nH8fl8/O1vf+O6667j3XffzVkk165dyxNPPMFrr73W4/Ks4XCY6upqKisr2bJlC7/4xS8YPXo0b7/9Nn/7298O+/xAIEBlZWWnv23btvG3v/2NSy65hAULFvDyyy+zcePGHp+blJKmpiaqqqp48cUXufHGG1WrUdEv6Y6f4/G01o5xSynHCiGOBn4kpbz8sMGFGA+8JKWc1XZ7O60iWyWEGAGsklJOO1ycbPwcNZ9Djeuuu47bb7/9sGsLo9EoDQ0NpFIpAoEALS0t3HfffTz22GNpj/v444+7fH5LS0taS9XhcGC1WmlpaWH//v1cfnn6W3fddddht9ux2WyUl5fjdDrxer2MGjWq0/IcraSAyWTik08+4Ywzzki7/+abb6a0tJSGhgaCwWCniRC3262XrG3fld6yZQvr169nz549+mPPPfdcxo0bR2lpqX4OsViMSCSi51tSUkIymaSmpobm5magdYwxEolw44036rEeeughFi9efKi3XqHIBzn5Od4DfA14AaCtTMJJWSYyTEpZ1RanSggxtKsHCiGWAEuALgszdUUikeCzzz5LO+ZyudiwYQMmk4lgMNjlBIHdbqe8vByTyYTX66W4uJhvfOMbncTx1ltv1fJMO+7xeCgpKcHr9eqmsNFolHg8TlFRUUZ37nfeeYdoNEokEqGlpQW73Y7JZCIQCDBkyBDi8ThNTU0Eg0HdDDcWi2UcZ7z99tv52c9+RklJCcOGDWPMmDH4/X7279/PgQMHCAQCemu1vRP3hx9+SFVVVVqsp556irPOOouGhgaamprw+/1YLBZsNhuRSIR4PE5zczNlZWWMHz8et9uN2Wxm48aN7Nu3Ly1WXV1dxvdboShUujXrLKXc30EEen1lr5RyGbAMWluOPXmuxWLpJKhDhw5l1KhRABQVFfVo//Dcuem1xCZOnMiVV17JkUce2anLqJURiMViRKPRtPuTySSxWIwnnnhCL7nwve99j9///vf4fD4ikQixWIympiYaGxsJBoPU1dVhtVrxer14vV7i8TixWAyLxcLLL7/Mc889l/b6999/v25Sm0gkqK2tpbi4mFNPPZVRo0bp9mHthdHr9XLvvffy1FNP6eUMlixZwpIlS5BSYrFYcLlcuFyutKJfWuvT7/fr9XO0Y08//TS/+c1v9Pf+ggsu6Pb7rVAUAt0Rx/1tXWvZ5uN4JZCt8USNEGJEu251bZZxDovVauXCCy+krKyMY445hpEjR2btBj5lyhTuvPNONm7cyPe//33mzp3LiBEjst7+9/zzz3PKKafw9a9/nRtvvBEhBFOmTOnRUqFYLMbs2bMpKytjx44dvPfee/zXf/0X5513nl4xsDtoY4ThcJijjjqKX/ziFwD853/+J5MnT+5WjEzv68SJE5k6dSp79uzhhBNOUJ6Oin5Hd76xlwH3AqOAg8DrdGO2ugteAC4Cft3274os4xyWqVOncs455+itxXnz5uUU79RTT+Xtt99m5syZPe7md2ThwoVcfPHFXH755Rx11FFZxbDZbIwYMYLLLruMuro6HA4Hf/jDH3ocRwhBaWkppaWlBAIBjj76aHbv3p3z7p3S0lIWLVrE+vXrKS0tNcxZXKHoK7pjWVYPnNfTwEKIp4BTgAohxAFad9X8Gvg/IcRiYB/w3Z7GzRcmkwm73W7IQnCLxcK4ceNyXoNZVFSk7x/XurW5EAwGKS4uznkrYns08w8jHIMUir7ksOIohJhIa8txAa1uPB8CV0spvzjU86SU53Zx16ldHC9opJRUVFRQX1/f7SJbh4uXKx6Ph5qaGubOnWvY9jytRK1RKCdwRX+lO4vAn6R18fcIYCTwNO28HQuZruoiZxtLCGGocOSaV/ua1kahlVw1CqPFVqHoK7ojjkJK+biUMtH29wTp/o4FjVEFnkpKSvQZZiMwSjSsViuVlZVpJVazJRaL4XQ6Dd1XbTableGtol/SHXH8pxBiqRBivBBinBDiWuDltq2AuRUj7gP27NnDuHHjco6jeRwaRVFRkSFb62w2Gzt37jSkqx+LxXA4HIa3HJU4Kvoj3WlSfb/t3x91OH4JBe4ILoSgoaHBEKuxWCxGUVGRYa2qESNGdFp0nQ1Dhw5l+/bthhTYKioqoqWlOztDu4/VajXEeUih6Gu6M1vdb53AXS5XToYT7QmHw8yZM4dt27axYMGCnONVVFSwZcuWwz/wMGgLxDPtvMk3WoXFVKrfVPJVKHS67FYLIY4RQgxvd/tCIcQKIcT/1x+609DaFTZiiQu0imNpaalhrSCbzWZYdzOVShlm7mBkXsXFxezZs6dgTYIVikNxqDHHB2m1KqNtL/WvgccAP23b+godbQ+yEbQ3gzVCiIyqISOlJBAIGFJOVQiBy+UiFAoZMgNus9mIxWJqOY+iX3IocTS3K7/6fWCZlPIZKeXNQPf2lRUARheoLy0tNaQLa4RBLaALWVFRUU5xgsEgRUVFOBwOQqGQoTZjShwV/ZFDiqMQQhuTPBV4u919/WKjrNF2/1rLyohxTKNqOu/cuZOpU6fmHMfv9+P1eg0tzQqt5hdGlaBVKPqSQ4njU8A7QogVtLp/vwcghJhMa9e64HE4HIYXdiq0ZSmJRAKLxZKzoBlZKrY90WgUm81meFyForfpUhyllHcAP6PV0fsE+e+rz0Rr6YSCx+l00tTUxJo1a3KOVV9fz/79+zv5RGaLEV1qQLc3e++993KaFQ4EArqbj5ElaDX7NkX/RUo5KN3cD9k9llKuznBsR++lYywHDx7k4Ycf5vrrr6elpaXH43Kaoez+/fu5+uqrdWGcP38+5eXl2O32HucUDodpbm7mxz/+MdFolOOPP54jjzwy7TFut7tL2zHNqRxav7R33HGHvpi8vddiNhw4cEB3HBo3bhyXXHIJ0DoEUFFR0W3BTCaTHDx4kBtuuIHly5fz2Wef8cc//jGnlumuXbtoaWlh8uTJPf4cU6kUzc3N+hpOi8WC2+0mkUjoDueZzk1KSTKZJBqN6n/t39/KykosFgu1tbWMGzcOp9OZcTeW5nzUHdOSUCjEhg0b8Hq9lJWVUVFRcdiWdzQaxefzpeXWsbicyWRiyJAh3TIAiUajncyJm5qaKCsrw2w267WGBjqHLZNQCGRTJmHVqlV8+ctfTjv25JNP4nA4SKVSncbBtA+7rq5OH6f0eDxMnDiRXbt2cf7556c9funSpRx77LEkk0m8Xq8eV0pJKpUiEomQTCaxWq1px6PRKN/73vf0OJpDeVFRkZ5DMBikvr6eQCCAEEKvGxONRnE6ncycORO73c7777/PVVdd1enczzjjDIYMGaKvM5RSMmHCBIYNG4bT6cTpdOozyeFwmGg0yttvv817772XVq9648aNVFRUEI1Gqa6uBlqFUvN/hNaLsLi4WBc+7fv06quvcsstt+ix5s+fz5/+9Cf9dmlpaY9s0V577TUWL17MOeecw9KlS/XPrL3xrkZ7YUgmkzQ0NOB2u3G5XAB89NFHfPHFF3i9Xo499lhisRhSSlwulz7uqmGxWPSSEHa7PW3FQmVlJbW1tVxyySUce+yxPPjgg5SUlOBwODrl09jYeNix6oMHD/Lcc89xxBFHcM455+BwOPjss88Ih8MZz9NisVBRUYHD4aCsrAwhBM3NzfqPpVZOWEpJOBymqKhI3/3UPo42wVhSUgK0rjLouB31gw8+4PjjjyeZTLJ//34mTOi3y587klOZhH5HLBbjuuuuSzt2xhlnMHnyZP3LXVVVpXdDtW6DEII5c+YwbNgwACKRCH6/P62uisaqVavweDyYTCb8fj+RSASTyaR3SbVSB8lkErPZrLcAXn/99bQ4oVCIk046iTlz5pBMJmlsbCQej2O32/F6vQghGD58OOPHj8dutxMIBNi8eTOff/45r7zySqe8TjvtNObMmaM7iI8aNYpgMKiXQdC6uG63G5vNhtPpxOFwUF9fnyaMAB9++CHz588H0Ftr8Xic8ePH67c1s9z260mbm5vThBFa69Fo70EqldKHAwD9/emqVSOl1FvpzzzzDPfffz9Dh/67wobf70/b2dOxxTR//vw0wbPb7Wzfvl2vpTN//nxdQJqbm9O2TyaTSUKhUMb1stXV1Tz77LPU1NTw4osv8pvf/AaHw9Ht1QzJZJJ4PE4wGMTn8/Hb3/5W/0xTqRSLFi1i1KhRXa7VTSQS1NTUIKVkz549CCHweDwUFRURjUZ59dVXmTlzJk6nkyOOOILi4mKklFRVVaW1MsePHw+gewfE43EqKyv1+0OhEMXFxezatQuXy2XIsrH+wIBsOTY2NrJv3z5mz56tH9u6dav+oQohGDp0aMaLMRKJ0NjYmPbl2bFjR1or9Mgjj+Sjjz7qtrdjKpWioaGBeDzO6tWrOXjwIFdeeSVFRUVce+21TJgwgeOPPx6v10tpaWla1yyRSBAIBDotPm9sbCSRSHD22Weze/du/fhvf/tbpk1rrVkWDodpamrC4/Hgdrv1c4rFYnrXXKO2tpZ//OMfbNq0CYBXXnklTYA6cqjvjd1u5/PPP2fZsmW8+uqrXHPNNfzXf/1Xl7PWiUSChoaGjHu6tdbhQw89xKhRoygtLWXBggVp7uNer7fHXfZt27YxfPhwIpEII0aM6NFzNSorK/H5fDz++ONYrVbGjBnD/PnzO9Xpbv//9sJtNpux2WwUFRVRWlqKx+Nh9erV7N69m+OPPx6bzZZWEK2n7Nq1i2AwSE1NDR6PR/cQHTZsWI/8NaWU+Hw+du7cic1m61Q2pJ/T5Rs7IMVR6/K88cYbvPLKK/zxj3885IXenXhvvvkmkydPZv369XzpS1/SHcZ7ilYfeu3atUgpGTlyJCNGjOhUZbA7hMNhHnvsMZqamrjrrrvYtGlT1nl9/vnnrFq1iurqatatW8ezzz6bVZxUKkUoFCIQCCClJJFI6K7l2aJ146WU1NfXZ+2enonKysqsy2dUVlZSVFREZWUl06dP13sJ+cypPVVVVQwZMsSwEhVSStavX5+zq36BMbi61VrX9sgjj+Tjjz/OefuaEIJZs2bp3btcBqPbi7QWJxthhNbZ+DPPPJPNmzfz7W9/O+ucACZNmsSkSZNYv359TuJjMplwu92GLgsqKSmhqKiIpqYmQ4UxkUjkLGbFxcV6l77Q3M614mhGIYTQu+CDgeyuyn6Cy+UiFosZZnywe/fugvxyaFZj5eXlhsQrxN6ENrlkJHV1dYb4YA4mjPqO9QcGtDhqtVWM8idsaGgw7MsRCoX02dNc0c7TyAXqg8FJR0qZdatdMfAZ0N8MIYRhLQ6bzUY0GqW5udmQLmN9fb1hFfmcTidjxoxh+/btOcfyer1IKQ0x4h3oGGlsoig8BrQ4amMuRrSCKioq8Pl8BIPBHtWF7opgMIjH48k5DrQuvRg2bFinhbvZ4HQ6MZlMBSmORnf3c13IXF5ebthOp95gMCzU7k0GtDgmEglcLldBFngKhUI5T+5A6wRIfX09I0eONKQVU19fr8bhuonJZNJLzyoGHgNWHJ1OJz6fD4/HU5DiCK0TKdlsQWxPaWkpDQ0N+gL0XLHZbMq9u4cU4gSWIncGrDiWlJTg8/koLy8nEAgYEtPoiyASiXR7IXlX2O12fYG4EeLo9XpJpVIF2a0uZApRIAsxp/7EgBVHbSuV3W7vtBskG+LxOCUlJYaJhhCCSCSSc8uxfTyj4rhcLsMLbRUaRgpHWVlZQdbwUeTGgBVHaL3QvV6vIS3Huro6vF6vod6ERtXU7o0WwkBvdRi5LMvhcBhWyE1ROAxocQQYNWqU7iiTC1otlEJ0ta6rq6O8vBwhhBor7CbxeFyZ8CoOyYAXR7vdbsgsrtbSMKL7mkwmDV18rE3sVFRUUF9fn3M8Nfvac9R7NvDIizgKIa4WQnwqhNgihHhKCJHbrMQhsNlshsxWG7mbwsjdMfBv5xqz2WzIgveORqmFQiHmBIU5BFGon2F/os/FUQgxCrgSmC+lnAWYgR/01ut5PJ6cu5qJRIK9e/cSCoWorq7mzTff1I1OsxFev99PPB6npqYmp7y03IwyFwiFQj0yoO1rjFxwbZSgGTWm3R4jcjP6B3gwki9XHgvgFELEARdQeZjH9xgpJdu3b2fVqlWsWLGChQsXpg3A98SJOhAI8J3vfEe/PWXKFNauXYvT6aSuru6w4muxWCgvL9ddW/7whz9w1113AfDTn/5UN+bV7Ps7OklnQrPG//Of/8zPf/5zHnzwQRYvXsxll13WrXPKxPbt23n88ce55557+PGPf5yVzZi2xTKVSpFIJHj66af5yle+wmmnnZZVTvX19ZSWlrJ8+XJ+9atfcdtttzFr1iymT5+eU8to9+7d+nuXC263m/fff5/nnnuOO++8s8txTM0JXkqJ3+/XBTWVSul/DQ0NLF++nBEjRnDTTTcd8nW1LZ7Nzc0Z76+pqaGiogK/3095eXnOpXv/+c9/4vF4mDVrFjabbVDsSc+Ln6MQ4qfAHbRWNXxDSnneoR6fTZmE119/nUWLFum3b7vtNr7+9a8D/3avTiQSXdZGaX9s7969nHPOOWn3z5o1C7PZjMlkwuv1pn35tOdqDuNaDRbtC/XQQw+lxWppacFut1NbW4vP5yMUCtHS0qKLTEtLC8FgUG9RpFIptm3bhsVi4Z577kmL9cYbbzB9+nT9tbTn2O12ysrKDvmlvuCCC3jiiSf020ceeSS//OUvDylCbreb4uJiSktLsdvt2O12PB4Pa9euZeHChfrjzjzzTBYsWEAsFsPpdOL1ejGbzfqfEIJAIIDZbNZXBZjNZuLxOI8++iivvfaaHuvZZ59lwYIFWbWw6uvree2117juuus4//zz+clPfsKCBQv09yoej3f60+j4PiSTSc477zzeffddoPUH97e//S1DhgxBCIHJZKK4uBir1aqPBQshKC4u1l3ktcX7JpOJcDjMnDlzgNaSCYdCCKFbw2X6fDRPyPr6evbu3Zvxh65jmQltaEb7HLXPAOD666/n3XffxW6389RTT+lu+dni9/t1J/o8Uzh+jkKIUuCbwASgCXhaCHG+lPKJDo9bAiwB9KJPPXydtNu33HILM2fO5Mwzz9Rt7zVHbE3EYrEYoVAorWZHRUVFxi/q448/rn/xg8HgIdc/HjhwgE8//ZRYLNbJ0RtgxYoVugu4VqdEK40QCASw2+0UFRWlCds3v/lNgsEgGzZs4J///Gfaa7W3abNYLJSWlpJKpfD5fF22cjdu3JgmjACXX365XiahPXa7HYfDgcPhIBQK4ff7aWxs1NeTCiF49NFH055zww036DP92nO0FlMymSSVSjFs2DC9ZRWLxUgmk/p97dmwYQPHHXdcl+/3oaioqGDSpEkA+mL3yspK/ftisViwWq1YrVZcLleXxbegVRzbm9KGw2H27dvH/v379ZUDTU1NRKNRKioqsFgs2Gw2bDYbLpdLr/+iCdOuXbv0WH6/n+nTp2d1jtD6o1VZWUkgEOCtt97iggsu6PSYTG7lsViM5uZmvWpkMplECMFTTz3F3r17Wbp0KRs3buS0007LquVeU1NDOBzWt89qn0Uhko9u9VeA3VLKOgAhxLPA8UDalSmlXAYsg9aWY09f5LTTTuO8885j+fLlHHvssXz3u99l7NixfPrppxQVFeF0OrFarQwfPlx/jvaldTgcenF77YJ3Op16UalHH32UGTNm0NzcTCQSOeRsuBCC2bNnc+yxxwKtX74RI0bohbG+973vcfLJJ+uPbf+FtVgsjBo16pDLhx566CHOOOMMxo8fz6WXXsq3vvWttPvj8ThNTU167ZeumDdvHieffDInn3wyNpuNP/7xj5x99tmdnIO0H5FIJEJ9fb3uft1xzeAZZ5zBmjVr2LhxI2vWrNHPPxs0R/FNmzbxj3/8g9ra2pycss855xzWr1/P7Nmzc+qam81mli9fztSpU3G5XJ3qFnVEa5lGo1GCwSCNjY36D7EQgqqqKi699FKWLFnCEUcckXVe0DoW6vV6GTlyJJMmTcppCVoqleLGG2/kq1/9Kh6Ph7KysqxjJZNJxo8fT319fUGbdkAeutVCiOOAh4FjaO1WPwKslVL+oavnZNOtbmho4PPPP2f48OG43W7C4XDWJQQ0/v73vzN06FC9RotGT8u0trS08PrrrxMMBrnwwgtzyqmyspKamhqmT5+Oz+fL2V6/qqqK2tpaKisrWbRoUdbiUVlZSVVVFcFgkFAolFOsdevWUVZWxmuvvcb555/Pjh07BppVP4Be1EqrJlhIVFZWIqXM+RqqqqrSu/jt/59HCqdbLaVcI4T4B7AeSACf0NZCNJLy8nKi0aguFu2rqWXL5MmTGTFiRM4CZLVamTJlimG7KiZNmmSoS/aQIUOIx+M5taq0bmVRURGTJ0/OKVZFRQVNTU2cddZZXY6xDSQKTRgBRowYwfr163MWR41YLFaQGyrak5cpJynlz6WUR0gpZ0kpL5BSdh6IK0BSqZQhs3Rms9kwpyBt146RLuBG4PF4DNuHPmTIEOrq6gzfj67oPkbu3deKsBX6UqOBPx+PcWu+wuGwIbNrmg+gUVgsFiKRSEEVeCoqKtJn2HO9sFwulz7eu3fvXsaNG2dEigVHeXm5ITucegsjhuBKSkr0sdZcHal6m0Ehjk1NTYYUIjdKHLWJF6N+jS0WCzU1NYZ2x3LNTXP3iUajhvppBoPBnNfsFSp2u73gegBG43A4iEaj/aJ+T2FnlyPtL3AjhMjIcRJt+VCuaDs0EomEoWM4RuWmLQkxAovFkvNYqKIw6A+f4YAWx96oOVJoH2pRURGhUMiwePF43LDtiCaTydAfFG3scSBTiPu0e4P+cJ4DWhyNxmhhLDShhdYhiKKiIkMETUpp6NhSIb5fgwkjJ2X6AwNaHI32NzTaK9GoX08jf4Xj8TiBQMCQsrFCCCwWiyHLjFQZVEVfM6DF0WazGTrAXchdASN/jQvR7qqiooJt27apyoh5pBB/zHuTAS2O7YtPGUEhd6sL8QtnZE5CCOLxeMEvHB7IuN1uQ2soFTpKHLuJUQvANYz8cvSHL5pReL1e/H5/vtPoNRwOh76ms9Aw0ruyEH/MOzKgxVFb+mEEzc3NeDweQ2INFowWbSklbrd7QFdGLC0tLdhKhoPpRxgGuDhqS2+M+JVqaWkp+O1ORmJEyQUjxy617YgD/QLVJhFTqVTBVTQ0qgyHUZN0vc2AFkcwbuF2KpVi586d1NbWGpCVsUSjUcMq6cXjcXbt2oXZbM6pC9XU1MTy5ctZt26dIXm53W6am5t1R+1c0KzojGyB5voDrHWlV69ezdKlS9m/fz8tLS26u1FVVRWVlZUcPHiw20NFUkqqq6sN6T1FIhE8Ho8hrdry8vKCtyuD/JVJ6DN8Ph9HHXVUTjHC4TB//vOfue2224BWs9Wjjz46q1ixWIx//etffPnLXwbgww8/1J2ou/v8jqa1K1eu5MQTT8xpD7nP5+ODDz7grLPOwmQyce+99/KDHxy6tI+UktraWmpraztdsA0NDWzdupU33niDkpISfvCDH6R5Z/aUuro6XnjhBdauXUs8HufOO+885ONjsViX+5QbGxv561//Sjgc5oILLtD3pJvNZoYNG5Zxsb9WxgDQPT/D4TDr1q1jzJgx2O125s2bl3XLNhAIcMcdd3DHHXcAcMopp7B48eJO8TRh13Jp/3yHw4HNZkMIwZYtW3jyySd57rnnePXVV/nSl76UVV4aK1asIBwOY7FYOOqoo0gkEsTjcWKxWCfHdI14PE59fX2nH45YLMajjz7K5ZdfnvV1pFFfX8/+/ft1B3UjGfDi2NLSkpMhg5SSTZs26cII8PHHH3drSYnFYqGsrEz3NEwmkzz99NP87Gc/0x+zYsUKfD4fQ4cO7TJO+73YVquV0tJS/ZxaWlq49dZbaWxs5KGHHmLs2LFMmDAhLf+O2Gy2TiUTampqeOedd4BWIbjiiisoLi5mxowZJBIJPvvsMyKRiJ5DcXGxXq522LBhnZxytm/fzhtvvAHA1VdfzZYtW7jtttv0shFWqxWbzab//1CiUltby+zZs/Xb3/jGN1i7di0+n0/PSTtXbW1lWVkZo0ePxmQyUVJSou+JDwQCTJkyRd9VpDmwt7S0EI/HO/UMXC4Xo0ePxuPx6J95MBikpaWFrVu38vrrr/P5558DcNNNN3H22WenvffthxYO9X9oLe2hcemll/L222/rP0DJZFIXpGAwiN/v18W8vLwcq9XK/v37SSaTuN1uqqur2b17NwD33Xcfu3btoqSkBEgv46GVoggGg4TDYex2O263u9M1o/1QHnHEEZSWlnL00Ufr15b2GXbE5XIxZ86cThOZTz75JI8//jjTpk0jkUgQDAYJBAJIKbHb7RQXF6fFtVgsukO7tg5Xez98Pl+vzQUMaHHcvHkzM2fOzCmGEIKZM2dy991366J20kknUVxcTCKRoLa2lkgkotv5t2/RJRIJ1q1bR1FRkW6WcMIJJzBz5kw+/fRTAN2hu7S0FPj3hdJeLDoKXPtxVKfTidlspqWlhXPPPZc1a9YctovX3NzM3r17Oz3u9NNP5+6779Zvf/7557S0tGC1Wpk2bRoOhwMppV5Ea968ebr9f8d6JB1biXfddRdlZWVIKUkmk3qLIxQKEY/H03LpeO6aKWpVVRVLlizhm9/8JiNHjmTmzJkZjUC0iz0SiZBKpfjiiy/0bqtml6UxefJkXC5X2gWm5aKVSqiqqiIQCOittdLSUoYNG4bL5SIQCOjieOaZZ2IymfTSAu3HvA8cOJBWbVIThPafZ0dDDafTyYQJEzjxxBP1WjOaoYfX69UdbhoaGvS65Vp9oY8//lgXx+3bt+P1ejnhhBMoKSlJ+5y0YafRo0fjcDiwWCwZx4r/+Mc/8pOf/IQvfelLjBw5kk8//ZSSkpJOIlpcXKyfRzwep7q6Ou3+qqoqEokEv/zlLykuLta/K5MmTcJkMhGNRgkEAmmO6dqPgraUS6s/ZLFY9No7a9euZcqUKYYYzGjkpcBWT8nGCRxam9xG7PTQKhnecccd1NfX8/3vf59jjz0Wi8XC0KFDcTgcerEt7Ut8KK644greeustqqqqeOedd3LqWsRiMVatWoXf7+f000/PadJISsnixYuxWCxs3ryZN95445C/ytovd6bB9fr6eg4ePEhlZSWXXHKJIRMpO3fu1EtY5GK6unXrVmKxGBMnTsTr9Xb5uPY1htqLv9/vp6WlhaqqKioqKrj99ts56qijOPvss/VeQkNDA4lEQn/uqFGj0noHbre702t//PHHbNmyhWefffawY7WJRIKGhoaMjkdaPZpEIsGBAwf4yle+wr59+zCbzYwePbrLmKlUinA4TDAY7BRXSskjjzzCVVddxZQpU7qMob03HdEKvO3evRubzUZJSQlut/uQ59hHdPnFHNDiaDQ7d+7E6XQe8gvWHb744gu++OILvXtilNFDriQSCT766COOP/54Vq9e3aOx0K7QquAZxf79+7FYLDnZ6x88eJBEIpGzL2RVVRVut5umpibWrFnDtGnT9PpEFRUVWU0Erl69Go/Hk3OPpzc4ePBg1j9KkUgkbRLGyO9EjhROmYT+jFE+gk6nk1AoRGlpacEII7SOkWpjh0YIY2/QcfwxG4QQhhnmejwePB4PZrPZkAteG1MeaDgcjkISxG4x4JfyFCIul4uWlpZ+sUug0EilUgXleD6YGOhrTDuixLEHGLWg3Mg9qor80RtOQYVo+gH9oyCW0Shx7AFG1X4xm82DotXYG+dYSMJRVlZWsFv9jKaxsXFAdvcPhRLHHmBk1cDBII7xeHxAVwo0ulBaIZNMJgtqfLwvUOLYA7Q1bEYwGMQxEokYWmHOqLo7hc5gOMf+gBLHHmA2mw1rKcRisQHdqgIMrzFdV1d3yJ1E+cLIve2KwkGJYw8wsuU4GMZwjO6KGWl2a2TrTNnZDUyUOPYAI1uORpvnDgZisZghLbSSkhKamppyT6iNwdALGIyoq7MHGNVyHIyD20bQ1NRkSGvb5XIVrNs2FNaM/GBGiWMPMKrl6PP5BnyXujdIJpOGLQBXAtQzBuMkkRLHHmBUyzEej6uWo0JR4ORFHIUQJUKIfwghPhNCbBNC5ObE2UcMpnVthUZjY6PuR2gEhdwSKsQlS4OxpZ2vluO9wGtSyiOAo4FtecqjR9TV1VFXV5dznFQqRU1NDXv37qWuro54PF4QNTU0m62tW7caEi8YDNLY2GjIexYKhairq6Oqqkq3ESsEHA4H27dvNyRWIBBg3bp11NTUpHlO5kIhCm1/oc/7dkIIL3AScDGAlDIGGLtB1WB8Ph///d//zaOPPgrAww8/zAknnHBIl57i4mJcLheJRCLNKv6LL77g1ltvZeXKldx6662cfvrpXHnllXz1q1/NOr+9e/fy7rvv8tWvfpWSkpKsZ07r6+t55JFHKC8v58c//nGP7Kk049/23HHHHTz66KPceuut/PSnPzWk5RcMBtP2NLc3btX2JR/OLiwQCLB+/XpDXGLC4TB33303+/btY+bMmcydO7fHqxCam5tpbm7m8ccfZ+nSpcyYMYP//u//ZtGiRV1+BlJKAoFAl3VwUqkUBw4cwGKxMG/evB6fV0c2b97M0KFDOXjwoGGORoVOPga+JgJ1wF+FEEcD64CfSikNr7epmY4ezv9Pc7eORCK6e3R7PvvsM9577720Y6WlpTidTvx+P4FAQDeS0KrHbdiwQW+JlZeX692SvXv3snLlSj3Oyy+/zLJly6isrOyUkxbvcLz33nusWrWK+vp6Zs6cqRv8trfpN5vNlJaWpolG+9ewWCyMHz+e3bt388knn3D++ednfK1oNEp9fX2nvMxmM8OHD9eFYfv27fqPya233srs2bM55phjDnsuXbFnzx4OHjyI3+/nqKOO6tQa0m6nUim2bt2aZjTbnlAohJSSl156iUmTJjF9+nT9fWr/l0qlDnssGAzyyCOP8PLLL/Pyyy/zox/9iNmzZ+PxeAiHwySTSd38VctDE/H2+Xs8HoqKili6dCnQasa7dOlSTj75ZA4ePEgoFNLLCGgxtOe53e5OhihCCBKJBNXV1QwfPpyqqqpOJRky0dV3LZlM8thjj+H1ehk7duzhP6wBQj7E0QLMBa6QUq4RQtwLLAVubv8gIcQSYAmQ9QeidVnbf7E6vIb+r81mw+Fw6N587Rk+fHhaq+ett95i69atHH300Xi9Xrxer+5qLKXEZDIxbtw4pk2bht1uT2tJzJ07lxUrVvDWW28BcP755+fUgtm3bx8LFizA5/NxxRVXdPkFj8fj+P3+tEJI7R8biUSorq7G6XQSi8UoLi7m4MGDneLY7XZGjhx5WNGeMmUKK1as4Jvf/Cb3339/J9frrhyOysvLM7Z8tfdo586dh32/DtXibWho4J133uHkk0+mpaWFF154gfnz56eVIdD+Mh3TfmiEEJSXl3PRRRfxyCOP8Pnnn2M2mznyyCOJx+N6+Qqv16v3MA71nkkpmT59ul547cwzz8Tv9+P3+ykqKtJ/eDKVzOgYWyt3cOyxxxKLxXJuIScSCW6++Wbee+89jjzyyJxi9Sf63AlcCDEcWC2lHN92+0RgqZTy9K6ek4sTeDAYJJVKHdIOvzvs2LGDRCLBxx9/zEUXXZSTw3VLSwtut5v/+Z//4eKLL2bGjBlZ59XS0oLT6WTv3r1phbWypbm5mZ07dzJx4kRDusGrVq3ilFNO6dZjpZT4fL5OXWb498Xf1NTE9OnTc85L2/dthEXYm2++SSgUYvTo0SQSibSaOpnidzzWvnXf2NjIggULMtbG6Q5aL8hkMhm+pTEYDFJUVDTQJmcKxwlcSlkthNgvhJgmpdwOnAoYMwOQAbvdbshuiKlTp1JZWcnMmTNz3jNssVhYt24dyWQyJ2GEf7uTG7VDI5VK6Xb/RtCTOiFaa+xQ5FIeoT2aIYYRF/qMGTNyqmmjIaUkHA5nLYzQej5Gmn20p0BqvvQZ+ZqtvgJYLoTYBMwGftVbL2SxWAwpat4en8932Iv4UNjtdubOnVuQv8C51L4erBj1OWqVBRWFQV5WIkspNwDz++K1ClGANIzKzeihkUJ+zxSKvkLtkBkA+P1+w+r1qjVxCkUrShwHAKFQyLAxQoVC0YoSxx6QSCRobGxkyJAh+U6l11BdaoWilUEhjkZd8MFgEIfDMeiqsCkUg5FBIY5GjaPlusyiP6DGHHuOtitKMbAYFOJoFEYKR0tLi1q2MUAoLi7G7/fnOw2FwQwKcTSqW22kD6ORM8xGosYce06hO4srsmNQiKNRLT6jhUMJkUJRuAwKcTQKi8VCMplU43IKxSBgUIijkS00m81m+HZEhUJReAwKcTQCrbXodDoNcWmuqanJOUahk0wmVflZRb9lUHxzjVhm8frrr7Nq1Spqamp45513su5aB4NBdu/ezeLFi7nyyiv59NNPc8pLSsn777+Pz+dj3759VFZWUllZSXV1dY/POxaL8e6772b0ccyGxsZGampq9JwqKytpamrK6vOQUlJVVaXbmeVCKBRi7969VFZWql6AoksGdAm8tWvXsmbNGjZv3swDDzyQdZxt27Zxzjnn6C3G4447jiFDhmC1WgmFQoRCIRKJBD6fj2g0ipQSq9VKcXExVqs1zYXZ7XZTX1/PJ598wieffML111+f0znW19dzxRVXUFtby7Jly7j44otpbm6mrq6ODRs26Bd/+xw0c95wOExzczPRaJTm5mY++ugj7rvvPgB+/vOfc+mll2Kz2SgqKsLpdPZoeKKuro433niD1atXY7VaOfXUU4HWtaK1tbVdCmQikaCuri6jqevIkSM5cOAAEydO7PH7JKXE7/cTCoWIx+OsWLECj8fD3LlzDyuQ2nk7HI5OHpdCCJqamnSfT8XAYcCKYyqV4t577+WJJ55gxIgROYnjmDFj0kqyxmIxGhoaMJlMuFwuvfzAzJkzdS+9RCJBU1MTsVgszSk8Go3yl7/8RY81b948/vWvf6WJKPy7G28ymSgrK8Nut3cyf43FYrzwwgt67ZYlS5YghMDj8VBRUUFpaSlmsxkpJclkklQqpYv4vn37MJvNOJ1OrFYrI0eO5Gtf+5oujqNGjUIIQUtLC3V1dYTD4UO6qXe8L5VK8f7773PgwAE+//xzZsyYwciRI3G5XLhcLlpaWjKuDbTb7cyZMydjd3z9+vUMHTqUhoYGysvLCYfDNDU1kUgkaGho0CfLAoGAXg6hfWkBra6P1Wpl9uzZVFdXM2TIkIxmtB3PD1qFvaamRr9f+xE0m83MnTu3U76K/s2AFcf9+/dzxhlnYDKZmDJlClVVVT2OoV0E27ZtY8GCBbzzzjsAPPHEE0yYMOGwu2XGjh2r1xoJBoNIKbHb7dx+++08+eSTANx9993MnTu3S7PaZDJJU1MT0Wi0U26aIGsce+yxTJ06Nc2EQkqJzWajoqICh8OB2WzW/ywWC2azWS8HkEgkOOecc3jmmWdwu91pLeBDuYJnKnmQSqW49NJLqampwWQysXPnTmpra4nH47pvYaaYyWSS6urqTsellKxZs4ZPP/2UyZMnM3HiRJxOJ8XFxVgsFsrLy/VtndOnTz+sEUcwGDTEvDUejxMOh4lEIjnHUhQWfV4mIRtyKZOQSqUMmRRYvXo1yWSSuro6zj777Jzj7dmzR299jh49Oiczi40bN2K32wkGg8yf39kmMx6P6y0s6NzKay9uyWSSlStXMmfOHI4++uisc4LWKn8ej0dvIRYVFal96YpCo3DKJPQ1Rs2WTp48mYqKCsPqCdtsNubNm0d9fb1eLTBbhgwZQiKR6LL1YrVaeyS+X/3qVw0pW6rV7TGiFo1C0dcMitlqI9AEzKj90NpYVq7CqLFr1y6mTJliSCyFQqHEccDQ0tKiDG8VCgNR4pgnVN0XhaKwUeKoUCgUGVDiOEBQZqsKhbEocRwgqG61QmEsShwHAA0NDZSVlRkSKxaLqbWICgVKHPOGkS29qqoqRowYYUispqYmSktLDYmlUPRnlDjmCSNnq6PRqGGFvxKJhGGlIBSK/owSxzxhsVgMtcvqD9tAFYr+hBLHPOH1emlubs53GgqFogvyJo5CCLMQ4hMhxEv5yiGf2O32gnRyUS1QhaKVfLYcfwpsy+Pr5xWjJmT27t2rGzwYJWxGmWsoFP2ZvIijEGI0cDrwUD5efyBRU1PDunXrWL58ORs3buzSlLa7CCF44YUXDMxQoeif5Gta8n+BawFPnl6/IIhGo4et15KphakZ2EajUW6//XZeeql1ZGLx4sVceumluvs3oJdzOFxLVavR8sorr7BhwwaOPvpopk+fnuWZKRT9nz4XRyHEGUCtlHKdEOKUQzxuCbAEWh21ByITJkzI6nlSSt2/cf78+bo47t27F5/Px9ChQ/XHbdu2jUQioZdb0IRy6NChmM1mPWY0GuWtt97imWeeYdOmTcyaNYubbroJm82W+4kqFP2QfLQcFwJnCSG+ATgArxDiCSnl+e0fJKVcBiyDVifwvk+zcBFCYLVasVqtHHnkkYwbN46FCxdy1113dWohjho1Cmhdv9jY2Eg8HieVSrFp06a0ujgARx11FKlUiokTJ3LBBRcoYVQMavJaJqGt5XiNlPKMQz0ulzIJA53q6mqGDh1qiON5PB6nrq4Ov9/PmDFjDKmxolAUOIO3TMJAZ/jw4YbF0qoQGlEiQaHo7+RVHKWUq4BV+cxBoVAoMqF2yCgUCkUGlDgqFApFBpQ4KhQKRQaUOCoUCkUGlDgqFApFBpQ4KhQKRQaUOCoUCkUGlDgqFApFBpQ4KhQKRQaUOCoUCkUG8mo80V2EEHXA3hxCVAD1BqVjJCqvnlOouam8ekah5FUvpVyU6Y5+IY65IoRYK6Wcn+88OqLy6jmFmpvKq2cUal7tUd1qhUKhyIASR4VCocjAYBHHZflOoAtUXj2nUHNTefWMQs1LZ1CMOSoUCkVPGSwtR4VCoegRA1ochRCLhBDbhRC7hBBL852PhhBijBDin0KIbUKIT4UQP813Tu0RQpiFEJ8IIV7Kdy4aQogSIcQ/hBCftb1vX8p3TgBCiKvbPsMtQoinhBCOPObysBCiVgixpd2xMiHEm0KInW3/lhZIXr9t+yw3CSGeE0KU9HVeh2PAiqMQwgz8Efg6MAM4VwgxI79Z6SSAn0kppwMLgJ8UUG4APwW25TuJDtwLvCalPAI4mgLITwgxCrgSmC+lnAWYgR/kMaVHgI5r9pYCK6WUU4CVbbf7mkfonNebwCwp5VHADuD6vk7qcAxYcQSOBXZJKb+QUsaAvwHfzHNOAEgpq6SU69v+30zrhT4qv1m1IoQYDZwOPJTvXDSEEF7gJOAvAFLKmJSyKa9J/RsL4BRCWAAXUJmvRKSU7wK+Doe/CTza9v9HgbP7MifInJeU8g0pZaLt5mpgdF/ndTgGsjiOAva3u32AAhGg9gghxgNzgDV5TkXjf4FrgVSe82jPRKAO+Gtbd/8hIURRvpOSUh4EfgfsA6oAv5Tyjfxm1YlhUsoqaP1RBobmOZ9MXAK8mu8kOjKQxTFTPdqCmpoXQriBZ4CrpJSBAsjnDKBWSrku37l0wALMBf4kpZwDtJCf7mEabeN33wQmACOBIiHE+fnNqn8hhLiR1mGm5fnOpSMDWRwPAGPa3R5NHrs8HRFCWGkVxuVSymfznU8bC4GzhBB7aB2G+A8hxBP5TQlo/SwPSCm11vU/aBXLfPMVYLeUsk5KGQeeBY7Pc04dqRFCjABo+7c2z/noCCEuAs4AzpMFuKZwIIvjx8AUIcQEIYSN1oHyF/KcEwBCCEHr+Nk2KeXv852PhpTyeinlaCnleFrfr7ellHlvCUkpq4H9QohpbYdOBbbmMSWNfcACIYSr7TM9lQKYKOrAC8BFbf+/CFiRx1x0hBCLgOuAs6SUoXznk4kBK45tg73/BbxO6xf2/6SUn+Y3K52FwAW0tsw2tP19I99JFThXAMuFEJuA2cCv8psOtLVk/wGsBzbTej3lbeeHEOIp4ENgmhDigBBiMfBr4DQhxE7gtLbbhZDXfYAHeLPt+/9AX+d1ONQOGYVCocjAgG05KhQKRS4ocVQoFIoMKHFUKBSKDChxVCgUigwocVQoFIoMKHFU5AUhxI1tbjab2pZyHGdw/FdydXoRQpwkhFgvhEgIIb5jUGqKfoIl3wkoBh9tdmNnAHOllFEhRAVgM/I1pJRGrBvdB1wMXGNALEU/Q7UcFflgBK0lMaMAUsp6KWUlgBBijxDiN0KIj9r+JrcdHyKEeEYI8XHb38K2424hxF+FEJvbWqHntItT0fb/89tibRBCPNjmV2kWQjzS5sO4WQhxdcckpZR7pJSbKCwTDkUfocRRkQ/eAMYIIXYIIe4XQpzc4f6AlPJYWndR/G/bsXuBe6SUxwDn8G9LtZtpdcM5ss0b8O32gYQQ04HvAwullLOBJHAerbtsRkkpZ0kpjwT+avA5Kvo5qlut6HOklEEhxDzgRODLwN+FEEullI+0PeSpdv/e0/b/rwAzWrcwA+AVQnjajusGs1LKxg4vdyowD/i47blOWs0XXgQmCiH+ALxMq2ArFDpKHBV5QUqZBFYBq4QQm2k1RXhEu7v9Q9v+NQFfklKG28dpM3w41B5YATwqpezkNC2EOBr4GvAT4Hu0+goqFIDqVivygBBimhBiSrtDs4G97W5/v92/H7b9/w1ajUS0GLO7ON6xRspK4DtCiKFt95cJIca1jUeapJTP0No1LwQLNEUBoVqOinzgBv7QttQmAewClrS73y6EWEPrj/e5bceuBP7Y5spjAd4FLgN+2XZ8C63jib+g1VcRACnlViHETcAbQggTEKe1pRim1VlcayBkalkeAzwHlAJnCiF+IaWcacD5K/oBypVHUVC0Ge3Ol1LW5zsXxeBGdasVCoUiA6rlqFAoFBlQLUeFQqHIgBJHhUKhyIASR4VCociAEkeFQqHIgBJHhUKhyIASR4VCocjA/w/GLxdlhcbVUAAAAABJRU5ErkJggg==\n",
"text/plain": [
""
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax = plt.subplots(figsize=(5, 5))\n",
"\n",
"x = data[1].values + random.normal(0, 0.05, size=len(data))\n",
"y = data[2].values + random.normal(0, 0.05, size=len(data))\n",
"ax.plot(x, y, lw=0.1, color=\"black\")\n",
"ax.set_xlabel(\"Species 1\")\n",
"ax.set_ylabel(\"Species 2\")\n",
"\n",
"ax.spines[\"top\"].set_visible(False)\n",
"ax.spines[\"right\"].set_visible(False)\n",
"\n",
"pass"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.9"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment