Skip to content

Instantly share code, notes, and snippets.

@cmitR
Forked from ecomodeller/period-kdtree.ipynb
Created July 6, 2022 06:31
Show Gist options
  • Select an option

  • Save cmitR/1a8f02688759a2771e4ca9c7282bc27e to your computer and use it in GitHub Desktop.

Select an option

Save cmitR/1a8f02688759a2771e4ca9c7282bc27e to your computer and use it in GitHub Desktop.
Period KDtree.ipynb
Display the source blob
Display the rendered blob
Raw
{
"nbformat": 4,
"nbformat_minor": 0,
"metadata": {
"colab": {
"name": "Period KDtree.ipynb",
"provenance": [],
"collapsed_sections": [],
"authorship_tag": "ABX9TyPfRqMNxFGNBh8AU+vjuAYM",
"include_colab_link": true
},
"kernelspec": {
"name": "python3",
"display_name": "Python 3"
},
"language_info": {
"name": "python"
}
},
"cells": [
{
"cell_type": "markdown",
"metadata": {
"id": "view-in-github",
"colab_type": "text"
},
"source": [
"<a href=\"https://colab.research.google.com/gist/ecomodeller/628c519a420df79f811ab99f0dd3365e/period-kdtree.ipynb\" target=\"_parent\"><img src=\"https://colab.research.google.com/assets/colab-badge.svg\" alt=\"Open In Colab\"/></a>"
]
},
{
"cell_type": "code",
"execution_count": 43,
"metadata": {
"id": "ce_-7psZ529s"
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"from scipy.spatial import cKDTree"
]
},
{
"cell_type": "code",
"source": [
"n = 1000\n",
"lat = np.random.uniform(0,180,n) # Has to be positive\n",
"lon = np.random.uniform(0,360,n)"
],
"metadata": {
"id": "UAnSIkNa59lw"
},
"execution_count": 65,
"outputs": []
},
{
"cell_type": "code",
"source": [
"X = np.vstack([lon,lat]).T"
],
"metadata": {
"id": "vkQb23nC6VhW"
},
"execution_count": 66,
"outputs": []
},
{
"cell_type": "code",
"source": [
"tree = cKDTree(X,boxsize=[360,10000])"
],
"metadata": {
"id": "i6ctzjM_AmG2"
},
"execution_count": 67,
"outputs": []
},
{
"cell_type": "code",
"source": [
"d,idx=tree.query([355,20], k=100)"
],
"metadata": {
"id": "FI5-IAjBA11j"
},
"execution_count": 74,
"outputs": []
},
{
"cell_type": "code",
"source": [
"plt.scatter(lon,lat, marker=\".\", label=\"All\")\n",
"plt.plot(355,20,'k*',markersize=10, label=\"POI\")\n",
"plt.scatter(lon[idx],lat[idx], label=\"k nearest\")\n",
"plt.legend(loc='upper left')\n",
"plt.title(\"Finding nearest points fast with periodic boundaries\");"
],
"metadata": {
"colab": {
"base_uri": "https://localhost:8080/",
"height": 281
},
"id": "xivnap_iBKea",
"outputId": "b782757a-0c3a-4f9e-952c-d2f1624392fd"
},
"execution_count": 75,
"outputs": [
{
"output_type": "display_data",
"data": {
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
],
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOy9e5wcVZn//z7dM5NkYkiGBCFhMgmRi5oomMmSsMgCBv2KvyiaiIC3RYXIS7ywsruoYMSg631hfcl3MSgvllXCJQkXg3wXg4CATCQTEZNFbiGTTBIumUwuJJPMdPf5/VFVPadr6nKq6lR1z9Cf1yuvTHdXnTp1Lp/znOd5zvMIKSV11FFHHXWMLOSqXYE66qijjjrMo07uddRRRx0jEHVyr6OOOuoYgaiTex111FHHCESd3Ouoo446RiDq5F5HHXXUMQIxoshdCPG6EGJGzHsfFkJcZP/9CSHEA2ZrVwdk27ZCiO8IIXYKIV7O4nlhEELcL4T4x4DfbxZCfCfLOinPbrPnTz7GvWcIIbqVzxuFEGdELGO6EEIKIRqiPt80VC6IeX/k908Dw5LchRCbhRB99mB0/k2RUr5JSrkpaflSyl9LKd9noq7DEXb7npVG2VHaVghxoRDisTjPEUK0AZcDb5dSHhWnDLucCuJKAinl2VLK/7LLjf1uaUBKucWeP0UDZc2UUj5soFrDErXy/sOS3G180B6Mzr/t1a5QLaEWJKAqow3okVK+Wu2K1DrqY8UMaq4dpZTD7h+wGTjL43sJHGv/fTNwPXAfsA9YC7xFufa9wN+APcDPgEeAi+zfLgQec5V7CfA8sNsuV9i/5YGfADuBl4Av2tc3BNT9n4Gn7WffDoxWfl8APGU/54/AO5Xfvga8aL/P/wIfUX67EHgcuBboAb4DjAJ+DGwBXgFuAMbY108CVtvP2QU8irXY/zdQAvqA14F/9XiHM4Bu4Bv2e28GPqH8Ph64BXgN6AKuAnJR2hZ4G3AQKNr12G1f/wH73fcB24B/9qjfWXb9S/a9N9vf3wm8bLf7H4CZyj1DygXGusp5HZjietYxdr2d97sReFX5/b+By+y/HwYuCni3mwkYs67nTrfbbjGwHdihtoXdl8546QHuAA533fs5e2z8Qfmuwb5mCnCvPTZeAC5Wyh5j17XXbrN/Abq95ifW/PgGg+O2E5ga431GAdfZv223/x7lNaYMc8FbgN/bbbgT+DUwwfWuV2DN50NAg+v9g/phNPAr+/vdwJPAkcZ4sloEnajS+uTeA5xsN/ivgdsUYtsHfBRoBP4JKBBM7quBCVgS4WvA++3fLrEHeCvQAqwhnNz/hDV5DgeeAS6xf3sX8Cow154U/2hf7wzic+37csB5wH5gslLnAvAl+33HYBH9vfZzxgG/Ab5nX/89LLJvtP+dxuCC5dm+yjucYT/r37Em3el2XU6wf78FuMd+5nTgOeBzMdq24lr7ux3AafbfLcDsgDp2u777rF0nhyieCivXqxyPZ20B2u2/nwU2AW9TfnuX/ffD+IyxsDHr8czpdtstx1qE3mG3nUMqXwE6sMblKODnwHLXvbfY945hKLn/Afi/WAR0kl32e+zfvo8lDBwOTAU24E/u/wL8FTgBa9E+EZgY432W2u/zZuAILMHnmoC2NMUFx2KR/yj7uX8ArnO961N2O4zxeP+gfvg81pxsxprv7cBhxnjSNPFm8c9uvNexVrvdwN0+HfoL5Z4PAH+z//400KH8JrAk0SByf7fy+Q7ga/bfvwc+r/x2FuHk/knl8w+BG+y//9MZsMrvzwKn+5T1FHCOUuctrnfaT6WEcgrwkjJZ7nHay6OOOuQ+1tUm37QHaT+Wrtv57fPAwzHatuJa+7stdnmBk4AQUsZaTCQwPqjcsHLsa/4b+CpwlN1fP8Ra9N1S/cN+YyxszHo8c7pd/7e6xtIv7b+fAeYrv00GBrDIzbl3hkd5DVhEVQTGKb9/j8Ed0CbsBdj+vBh/cn8We4yGtGHY+7wIfED57f8AmwPa0ggXeNTzw8CfXe/6Wb/5E9IPn8W1Ozf5bzjr3D8spZxg//uwzzWql8QB4E3231OArc4P0mr1rQRDqyyNcoLKmgZcLoTY7fzDmmhTAIQQnxZCPKX8NgtL8vB69hFYEkGncv3/s78H+BHWdvsBIcQmIcTXNOqtoldKuV/53GXXcxKWBNTl+u3ogLL82sMLi7AmZ5cQ4hEhxCk6lRVC5IUQ3xdCvCiE2Is1AWGw/WKVa+MRrEXgH7Aku4exdjOnA49KKUsRyorSFlDZ504fgDWW7lL6/hkswj7S514VU4BdUsp9rrKPVn53P9cPU7GIWRd+7zOFoWNqCvqIxQVCiCOFELcJIbbZ4+ZXVM45d53dCOqH/wb+B7hNCLFdCPFDIURjhHcKxHAm9yTYgTXoABBCCPVzjLJalc9xywFrkHxXWbQmSCmbpZTLhRDTsPS5X8Ta1k7A2g4L5X6p/L0TS188UylrvJTyTQBSyn1SysullDOADwFfFULM9yjHDy1CiLHK5zYsXehOLMlkmuu3bZptoGJIPaSUT0opz8Hant+NJenr4OPAOVg7q/FYkiLY7RdQrk5bPIKl1jrD/vsx4FQscn/E5x6dcnWgjjenD8AaS2e7xtJoKaXaD3512A4cLoQY5yrbubdi/ti/+WErlt5aF37vs52hY8r5bT+WIAOAECKKd1QYF/wbVju9Q0p5GPBJKuccBPelbz9IKQeklN+WUr4d+Hsse9unI9Q9EG9Ucr8PmCmEWGhbuL+MtaWOgzuArwghjhZCTMAyrsTFjcAlQoi5wsJYIcT/Z0+ysViD6DUAIcRnsCR3T9jS4o3AtUKIN9v3HC2E+D/23wuEEMfag3kPljThSJivADrnBb4thGgSQpyGNTDvlJYr3R3Ad4UQ4+xF6atYEk9UvAK0CiGa7Do32X7y46WUA8Bepc5hGIdl8OrBIoJ/c34IKfcVYKIQYrxfwVLK57EW0k8Cj0gp99r3LcKf3CveLQG+KYRoFkLMBD6DZaAHy57yXbv9EUIcIYQ4R6dAKeVWLHXB94QQo4UQ78Qyvjp9eAfwdSFEixCiFcvO44dfANcIIY6zx/Q7hRATY7zPcuAq+z0mAUuU+vwFaz6fJIQYDVyt8542wrhgHJYKeI8Q4mgsG0IU+PaDEOJMIcQ7hHW2YC+WUBRllxeINyS5Syl3Yhknv4812Y/D8jSJgxuBB7Cs5X8Gfoulj47sLyylXAdcjGWx78VSm1xo//a/WF45T2ARwzs06nyFXUaHvaVcg2XYAuud12AN3CeA/yulfMj+7XtYE2m3EOKffcp+2a7jdiwD1SVSyr/Zv30JS5rahCXF3grcFN4CQ/B7YCPwshBip/3dp4DN9vtcAnxCs6xbsLby27AM4B2u3z3Ltd9pObDJbg8/VcAjWK6XW5XPAlgf4d3i4BGsPn4Q+LGU0jkg9h9YxvQHhBD7sN53boRyL8Da3WwH7gK+JaVcY//2bay2fAlr7P93QDn/jrUYPIBFYL/EMuBGfZ/vAOuw5tlfsdr1OwBSyuewbEhrsLyutM8PaHDBt4HZWALQfcAq3bJtBPXDUcAKrHZ5Buvdg9oyEhzviDoMQQhxNpaBdFroxcMUwjp99yspZWvYtXWkAyHEdCxybZRSFqpbmzpqEW9Iyd0khBBjhBAfEEI02Nu2b2FJOnXUUUcdVUOd3JNDYG3derHUMs9g6QPrqKOOOqqGulqmjjrqqGMEoi6511FHHXWMQNREoJtJkybJ6dOnV7saddRRRx3DCp2dnTullEd4/VYT5D59+nTWrVtX7WrUUUcddQwrCCF8TwfX1TJ11FFHHSMQdXKvo4466hiBqJN7HXXUUccIRE3o3L0wMDBAd3c3Bw8erHZVhjVGjx5Na2srjY3Ggs3VUUcdwwA1S+7d3d2MGzeO6dOnY8W2qiMqpJT09PTQ3d3NMcccU+3q1FFHHRmiZtUyBw8eZOLEiXViTwAhBBMnTqzvfuqo4w2ImiV3oE7sBlBvw5GPzq5ern/oBTq7eqtdlZrCG71dalYtEwd79uzhwgsv5Oabb2b8+MHw2/sPFdh/qMDYUQ2MHRX8ylGuraOOaqOzq5dP/KKD/kKJpoYcv75oHu3TWhKV17Gph3kzJkYuJ8m9pmG6XYYjalpyj4p7772Xu+++m9/85jfl7/YfKvDSzv28svcgL+3cz/5DldFR9x8q8Oreg2VSV6/N5/OcdNJJzJo1i3PPPZcDBw4Alj3gnHPO4bjjjuMtb3kLX/nKV+jv7wfg4YcfZsGCBdm9dBWQpkT0Rpe2oqJjUw/9hRIlCQOFEh2bemKX5RDiTx54lk/8oiNSHyS5Nw1EbZeROO5GFLnfdNNNFf+DRd5SSiSWgVEldzeZ7z7QX3Ht6DFjeOqpp9iwYQNNTU3ccMMNSClZuHAhH/7wh3n++ed57rnneP3117nyyivL5fYXSmx67XV6Xj+U1atnhrQmcWdXL1fe9VcuWPZEzRBEraOzq5ftu/toyAnyAhobcsybEZTkKBhJFgqTi4wJzJsxkaaGXEW7+BF4GmO6FhaLUL2DEOImrBRqr0opZ9nf3c5gRp8JwG4p5Ul2AoFnsDKeg5VV/BLTlXZw1lln8eCDD5Y/NzVZGcsef/zxCl3z3Hefzo3L70YIUaFqUYkf+38hBEiJEKIiUeJpp53G008/ze9//3tGjx7NZz7zGQDy+TzXXnstxxxzDN/+9rfZ09fPwUKR1w8VeN1eSCa+aVRaTZA5vCZx3O2us41vaW5i6eqNHBoolZNRJi3bFGpJ1aBCVTs05HOcd/JUFs1uTVRHhxAHCqXIC0Wce9Ns2/ZpLfz6onnl8oEKNc2SBTPpPdDPvBkTjY5pqB2VkI5S+WastG+3OF9IKc9z/hZC/AQrBZWDF6WUJ5mqYBCuvPJKnnjiibK6xFGNOP8DNDc3861vfpMjDxs9RI8+dlRDBZm3NDfR0txU1rk7KBQK3H///bz//e9n48aNtLe3V9TjsMMOo62tjRdeeIH9hyqz6+3pG9Ai9+Gi6486if0msDoBckJQchZZrAD5SaVQE8RRK5PUCyohFYsljp4wJnHd3IQYpbyo92bRtu3TWsplXv/QC+X26i+UWHLPBkpSlok+7qLmBdOLRVyEsoiU8g+2RD4EdnLljwHvMVstPZx55pmsXr2aBQsWlAleRXNzM/fddx9nnHGG5/1jRzVwzKSxQ0jV+b+vr4+TTrLWqdNOO43Pfe5z3HDDDYF1GjsqX/F5/Jjww0OOekjai8wxk8bWLMFHmcRBE1idAEhJLicQSPL5HB9tb00khZoijqwnaZQFyb3ItjQ3cf1DLySWglVCTPPerNtWbS8hBMWSJUwMFEr0Huj3HdNxhIQgASjLnWBSBjkNeMXO/u7gGCHEn7GSvl4lpXw04TMCceaZZ3L77bdz7rnnVvhzjx49mttvv92X2B0EScpjbJ27ire//e2sWLGi4ru9e/eyZcsWjj32WHbt2sXohjxvGtXA+DGN2lK7qh5y7xxqDbqTOGgCuyeAuk1OOuhNEUcSNUVURF2Q1EXWUWu5761VlRJk27bg3V7qs73GdBwhwWlzr/HsVqUlFWLCkJRBLsDKDO9gB9AmpewRQrQDdwshZkop97pvFEIsBhYDtLW1JarE7t27aWhoIJfLMWrUKA4dOkRDQwO7d+9OVK4X5s+fz9e+9jVuueUWPv3pT1MsFrn88su58MILaW5uBqCpIceMI96kXaZbPVTLxB4FLc1N5IQA5JAJnEQFEIYgqRbQfmaadXQjzoLkEJKqclCNmUl3L1nqxJOWr1NXlcBPOGpc6PVR+yRsMVDL6y+UWL52C6vWd6em7ovNIkKIBmAhUFZASykPAYfsvzuFEC8CxwNDgrVLKZcBywDmzJmTKNffL3/5Sw4cOMCJJ57ID37wA6644gr+8pe/cNNNN/HJT34SiK7T3n+ogLT/V68XQnDXXXfxhS98gWuuuYZSqcQHPvAB/u3f/i12/f3UQ9WAqQnd2dXL0tUbKZYk+ZxgyYKZQ8pLogIIgp9U25DPgZQUSnKIUc1Lalu5vhsBLIwoXcXdyjfkBANFq72SGjOT7l6y1onr1inMfqNbV51nR91dhLW5U57jOOCohdJSSSVhkbOAv0kpu50vhBBHALuklEUhxAzgOGBTwjqGYvz48fzoRz/isssuI5fL8Z73vIfrrruORx+1NEI9rx9i++6DQLBO21kA8jnBjj0H6fhbNy/t3D/k+qlTp1b40qs444wzQlVBXqg2qYPZCe0MdAkUS5IN2/eE3mMSflItWJPKbVRT37Wzq5cLlj1Bf9GSOe7s7Gb5xXptkagN7V0OEU8V+0nBTQ05+gdKZWeBKNBZHLJU++jab0ySZdTdRdhi4JS3an03d67bSrE0dEdrEjqukMuBM4BJQohu4FtSyl8C51OpkgH4B2CpEGIAKAGXSCl3ma3yUNx9990Vn/P5PJdffjmXX345+w8V2L77IGVfDB+dtmrUBBF6vQ6GiweMA5OTxJFE+4tWS67o7E5VvwjeZKNOuLwtuRdLcohRTX3Xjk09DBQHN5NR2iJuG3Zs6qFQtBfDYvS2d0ui7dNaWLJgZnkBW7p6IyccNc6IURCGkq2fjtkU+Uex35gkyyi7C53FwClv4ezW1BdGHW+ZC3y+v9Dju5XAyuTVMgdLvaJqfbx12qpRUyCxZSiwdeBx1DrDxQPGgclJ0j6thXPnTOXWtVtiE1YU+El2Xv7OfkY1B/NmTKQxL8qSe5S2iNuG82ZMpCE/uAiZIKjeA/2UpIy1WIcRVYX+eGDoLgjM6vyD1FZZ2kbCoLsYpKWSVFHbbGMAY0c1kBOiLJFPmTDak2TzOQFY7njlbbG0JvfBgSI79hyMRNTDzQMGzE+ShbNbWbm+O7FPfBg6u3q5bs1zZV2mm8i8pFrwN6q1T2th+eJTYuncE7WhlJX/J0TSMwlBBOR2LXQvIoBRnf+SBTMD1VZZkOVwQ22zjQHoGCv3HyqwY4+lkwfBuFEN7Ds4YP0oJXv6BiIT9XDwgPEiU5OTJArRxdVVq/dJrHgauhJz0LvGbYe4C1THph4KtpqoWJJGdjlptn+YayGQaBfoVsPcv2FHJLVVkPHV/b2pA2+1sHNQUXuMkwLCVClulUxDXlQQ8/gxjRzoL0Yi6qgeMFnr57M6falLkkl01c59OQGnHjuJy846vioTLEmbpqU3TrP9w1wLk+wC3e1x9qzJPLl5l1b7+PWD1/dAxUnppefM4uNzo7lmh/V7lIXGJN4Q5B5GnEFhCPI5y/A2efxoiiUZiXxrWT9fK0ekHcQ9cem+r1rEDsnatNp646SLi9cikmQX6NUeYb7pDllu293n2Q9+wc2c70pSsuSeDZEMzxDc71EWGtN9PuLJXYc4g6TsX/zqdi676JPc8/CfmD/vXby2o5u5CxawYcMGHn74YX784x+zevXqxHXMWj+f9QnBMHht83VOXKZBinElqgrPnJxg++4+Ort6IxF8tRamLBaXqO2qei85n/3uc5/+bLCFMnVs+415J7YRQCmGSixoLvkRfxbC1RuC3HWI00vK3n+owP13r+RdfzeP3969gnnveofRelUsJvbOgYz0815eJCZikyStU/u06Ccu1UmfdKubRKJy2vSGR17k9397lVvXbmFlwAnEWtPTprm4xD3Kr3uPSpbFYonzT25jyoQxWoLA0nNmWd4+JUlTY7xdi9/C6Ef8WQhXI57ckxg2+w7sZ/2THfzijnv58mcu4LvXLPW9Nmq2J3U3MXn8aO06mYQzmU35LJsiq7gnLk1sdU1IVA/97VWKJUsS7B9Ir64mkNUCE+dQVJS+cI8ZPw8nrwXs43PbtMIRBMFvYfQj/ix2SiOK3L0GapKj/fetvpdTz5jP9BnHMqGlhfWdnRw3bcqQ66LqzN27iT19AxVucEnVMlGNsyZ8lpN4u+gO/DBJxwQxJ5WoOjb1lLf4ADmfUAJZ2TyCyDvLBSbqoahfXzQvUl8kJUuTuxZdl9K01XAjhtyDBmpcD5TfrFrBRz51MQI4+0OLWH33Cv7pK18ecl1Unbl7N+HnjRPHg8ZroQmDCZ/ligWiUOK6Nc+FGjeD+sw98HUmr4mtblKScOqgel+kVdcwRAlklbZRPaxdvepy6ZnHRuqLatosHNTKjgxGELl7DY63HjUutnvhrl27eOThh9i4YQMIKBVL5HKCy778pSHXRlX9eO0mRjfmKz7H9aDxWmjC4GXMjOqzrJJaScJjz+9k7Uu7hoQ1VaWaqOQSNnlNbXVNe3mkWdcghLWv1wKTdiRIXbWKM8ZqgbCjoJa80EYMubsHx0lTJyRyL1yxYgWf+tSn+PnPf17+7vTTT2fr1q1Dro2j+vHKChWUAvCVvQfL2aTCynUvNDrZLFX9+8LZrUNOZurEzPj1RfO4bs1zPPb8znJwLjWsKQxNdWYiq5PXe+ggLSLTrUPaxBW2O/AyqldD6gyKgT7cUEteaCOG3N0DdWrLGF7ZezC2e+Hy5cu54oorKlQjixYt4nvf+57n9aYPH6kkLYHXDxU44BGh0oFaz7g2BveWcuHs1vJvOkTUPq2Fy846nic37/IMawqV6p2gDDhhdUtKPLeu3eIbFTIpasULRmd3oParn5dS2u6RtaLGMAGdNs9qfIwYcofKgbr/UCHR8f+HHnpoiGrkc5//Al/+8qDOPW54Xx04JP3K3oPlRNtSI6Kls0t582HRPXBMbCmdwe0X1tQt1ehKrya3u51dvSy5ZwMFx6PF4Pa51sgqyu7A6yDZBTd2lD/rhj2OgrB+rZWFMgqC2jzL8TGiyF1F3OP/zolUR+9dzeBfY0c1cORhozmgEHdYRMsk9TS1pXQGt1dY07h6ZpPb3Y5NPWVXRbAOsZjaPqetc80yO9Kq9d30F6wY+P2FEqvWdweSlul+rbWF0gSy1MmPWHKH6Mf/HRc2AWX/82oH/9JZpEwFKTNt5POSYOLqmU3Wbd6MiYxqtBJZ5HL+Hi1xy05L55p1dqSV67srfvOLVWni4JdXvyYlQpNZxUzNiSx18iOa3HXhSL4OHAm4WJJlYs3nRNnzJGuSD1ukTKbpq2XvBFN1S9NTJc2ys/bEWDS7lRXrtjJQlDTmBYtsG4zOYSOnvklcGJMQoamF0HRS6yy8pByMOHKP4xvuNl46krtD6E7avTQCe5mKBjlcsj1lhTBpK81FLK2ys/bEaJ9mxbVX21HnsFFLc1OsUAMmT3GaWgjd5zdMJLXOSoAaUWwQ5hvuRaTOd07UR0fnrhI6dl4m07r34ZityQTSNpLFjWNSS0GzvJCE7OI+301EOoeNohJrlMNsujC1EDrlpJHUuh7yNwL8DvC4pW+EoKW5kTGNeV+J/NW9BytivDtZmkzq3qttsK0GstAbmyQX53evSag7OU2+cxyyM/l83cNGUYg1DXWTKfVHmPdXXNREyF8hxE3AAuBVKeUs+7urgYuB1+zLviGl/K3929eBzwFF4MtSyv8xWuMAuA2L+ZyoSHqtSt+79vcjfCTy/YcKvPjSS3zqvEXcteaPZeNq1HjuUesbt9ynnnqK7du384EPfMBIvdJEFnrjqFJbUJ1MxOOu9qlFd+wgndAQftD1nY9CrKbVTeqie+mZx2pf61fPIO+vuMhiTOiwyc3Az4BbXN9fK6X8sfqFEOLtwPnATGAKsEYIcbyUsmigrqFwGxbdGZYYTHttY6hE7qhK9hyw0uwdPraJCc1NqUnUE5obEVDxjEKhQEOD/vOeeuop1q1bNyzIPQu9sUly8ZuEUfyz47xzGh4a/QMlSsDjL+zkyc27AhekoOfr7B6i7DBMGhmjLLpxUgsOB68qB6EMIqX8gxBiumZ55wC3SSkPAS8JIV4ATgaeiF1DXTx9Bzy4lLF7uhk7vhXmL4ETFlZIxpPHj+bgQJFdBwYqvlMlclUdg5Q05nO8sm0LH1m4kB9e+zPe/ffzyiT88MMPc/XVVzNp0iQ2bNhAe3s7v/rVrxBC0NnZyVe/+lVef/11Jk2axM0338zkyZO58cYbWbZsGQcPHeKo1ul89z9uoLl5LFd8+RLeNLaZP//5z5x66qlceumlXHrppbz22ms0Nzdz44038ta3vpU777yTb3/72+TzecaPH8+aNWtYsmQJfX19PPbYY3z961/nvPPOG9I8tXIYJCtvAVPk4jcJgyanc/q1WJKMarRIY8mCmdy/YQdnz5o8pF7uvjG9ZXfe77o1z/H4CztDpcVq+JebIs4oEnE1d1RZzIMk4ugXhRCfBtYBl0spe4GjgQ7lmm77uyEQQiwGFgO0tUXLWTgET98Bv/kyDPRZn/dshd98mbHAMScsHGJEndDs76HiqEocj5ltXS9y4ac+wZIf/Ywpx83iJSUEQF9/kfV//jPr1v+FY49p49RTT+Xxxx9n7ty5fOlLX+Kee+7hiCOO4Pbbb+fKK6/kpptuYuHChVx88cW8uvcgV111Fatu+xWf+MxiCiVJd3c3f/zjH8nn88yfP58bbriB4447jrVr1/KFL3yB3//+9yxdupT/+Z//4eijj2b37t00NTWxdOlS1q1bx89+9jPP5ukvlGLplNOC8wzHZU7nmWnU0Ss0a2dXb0XSEr9J6Pf9kNOvA9bhn5X2gaAnN++qSOPmRaRp6aCd0BBh0qL6/EMDwYeXag1RJOJqx4FJ22smLrn/J3ANll7jGuAnwGejFCClXAYsA5gzZ47f+Qg9PLh0kNgdDPRZkvw7P+ZJ4H5qFke1c3DXKPb09nDBuYu48ZblTJhyzBBD7ct7+ph54mz6R7fQN1DipJNOYvPmzUyYMIENGzbw3ve+F4BiscjkyZMB2LBhA1dddRW7envZs3cff3/6exBC0JATnHvuueTzeV5//XX++Mc/cu655w6G3+3vB+DUU0/lwgsv5GMf+xgLFy7Uap5DhWJknXKaiPrMNLxfwpIlq8/xm4Re33dsGhrP3Qmi5kWWXkSepg5aR1qcN2MiDXkrwqcE7ly31Tf5Ra0hikSc1S6yWohF7lLKV5y/hRA3Ak4S0W3AVOXSVvu7dLGnO9r3IRg7qoFJ40YxYfx42l3Yr3IAACAASURBVNraeOrJDs788IwKw+f+QwUk0NjUVI75ks/n2d93iJ37DvK2t7+dtR0dQ8q+8MILufvuuznxxBP5+Y2/5KGHH+aYSWNpyOcYO9aKvV4qlRg/YQKrH/ojvYoKaf+hAjfccANr167lvvvuo729nc7OztD3GdWQj6xT1kUcidotGa4MkQxNe7/4lQnhcevDUNZvFwbjuZ9w1DhWdHZ7kqUXkaehgz40YOV1XXrOrFAjY/u0Fj7a3srytVuQQDFGXtFqIqpKbri8V1Tk4twkhJisfPwIsMH++17gfCHEKCHEMcBxwJ+SVVED41ujfa+JpqYm7rrrLm5f/mueXHMvRx42uqySydsSGVjbl3xOMFAs8drrhzjsyDZ2vPwqv3/kUQAGBgbYuHEjAPv27WPy5MkMDAyw4o7bGN2YH7KLyI9q5qij27jt9juQUlKSkr9tfNry4nnxRebOncvSpUs54ogj2Lp1K+PGjWPfvn3+72ET3Fffd8IQonPIJS+IfQrwJw88yyd+0UFnV6/WffNmTKQhJ8ptt6KzO/DeqHX0I+6wMpO0hQOHmC9/3wnc/vlT+PjctjJZCvsahyzV69190z6thUvPPNaIV4bjo10oSZbcs0GrnxbNbmVUY7K2GMlw1He6Y74a0HGFXA6cAUwSQnQD3wLOEEKchDU3NwOfB5BSbhRC3AH8L1AALs3EU2b+kkqdO0DjGOv7hBg7diyrV6/mve99LxNbxvOhD30IoCLwlLA/F4olkNDQ1MRPfv5fXPmNr3Pg9X0UCgUuu+wyZs6cyTXXXMPcuXM54ogjmDt3ricp7z9U4Hs/XcZ3vnE5N/70xxQKBd7/oYWcffopXPIv/8Lzzz+PlJL58+dz4okn0tbWxve//31OOukkX4NqkITiFb/dC0lyXLrrcu6cqdxqS4aFYomlv9nIzKPHex7tbp/W4mmQ9Ns16Kg1/KTjtJJ9LJrdyqr13Z51MiE9BrVFPifKNoCS1JPCR7rKIgmGS0AzocZUqRbmzJkj161bV/HdM888w9ve9jbtMg6uv43Gh64ht28bwvGWeefHTFe1DL90dnFPnKqnZ9VyEILDmxt93TG9olm6r3O3pUMETtaluG5jYOmoHcKKMshVdYE6ApvyguWLTxli8NXVj7vfsZaIKa06hZFNmrHr34i4/qEX+MkDz1KSkBfw1fedEKrqSgtCiE4p5Ryv30bEccj9hwq8dOTZyPPeH+sYv258F/d1XsG64gTw8loodMpR73MIMhfy/ioR5DzypUZxG9PJcakuJGqWHUcav+ruv6LKFwPFoZJlHP14LepS06pT2A7q43PbOOGocTW32JlElot52M6wVgSLEUPucY/x68Z3Ua8DwZQJo5n4plGRPHGi1P/NGin1Ku6z4ZfQw4FKBEhJLmcd4vKKpa0OUHVA53OC7bv76OzqDSQst3QuoGzU+/jcNnoP9A+5pzE/NLa632RKEjGwFiafKeiqobxcPUcCslaTBKmsakllU9Pk7hBuGPI5QdzYL7oLQyWRSrbvtvT7E980SvtZfogbhsAvmqV6v1vt5iYCr7yVfgP01xfNY6UdY+PWtVu4s7M7MDuPs5CohmfHqHfCUeMqPEuEELznrW/mktPf4qlTj6IfD4oF48QIKZSqp6IwvbgE+dyneTiqVlCNw0h+Qk21Q02oqFlyHz16ND09PUycODGQ4PcfKrBjz0GcUAKTx4dLvCp0idX6fjB8gQS27z7o6e0SBrd6J248dvU+L527lJKenh5Gjx5MuadjKPMboO3TWiyjYHEwPV3QARf3sXcHJdtbJEit40VEbr2m1wTzcv37+Nw2Tx1/NSaf1+lVSJ6n1N0WWR2OqgVU+zCSX13y+RzblB1u1qhZcm9tbaW7u5vXXnst8Lp9BwfY21coS659rzXw6ujGSM8qFkocKhQZ1ZBny15/79BDhwrsPjBQJgfneeMiPK+/UGLn64fKeuYJzY0VZD7UaS8a3PePHj2a1tZKl9Aw3W/QZHGb34PM8epCsq9vgF889hKlkqSp0Yr57agHvIxRcYnIy/XP0TerxC6gHHs8KzWF1+nVleu7y+nsgqTpqNJ+FoejagWmzwUkKcepi7NDvO1PyeO/x0XNkntjYyPHHHNM6HWdXb180eWx8bYUG7HseWCTVNTnWZb2zTielA05we2fP6WmJKigyeKXnSeoLOf+9848SttLJy4R+bn+tTQ3VSxEc6a30NLcxNW/2UihmI2aYuX67srcrTkrzEXYIhYUmdKPiKIejhrudghT7qQm1FaOm3ChpOeskBZqltx1EXXVTjqIk3oezJsxseylAoMqiiw7PkqYU6/v3dl5omLj9j2hhBZXGmuf1sLSc2ZVuP7Nm2ElkMgJKElLav/zlt0USzIVFY1X+3Z29bKis7v8vLygfHp1pcv/XTeNXRAR+bVfkCprpOnio8Kk2qoWdknDntxBf9U2uTIn0Y2WycdDRZGFaiBpG8R5/4pclDlBQz5HsRg88NXnRFmU/RZgZ7IJe3F1q2hMTEC/9nUOMDnPO//kNj4+1wqYp5IwDCVtL6LQISLdfhqpuviocBv4W5qbYpdlUlUUFyOC3HVRK4NYJZ8oB4lMQG2D/gzbQH1usSQ57+SpHD1hjPZuK+qC5CY2dbI5be64dZ47Z6qxwFh+Uvad67aWF5PGhhwLFXWW89yOTT1s292nfZ7AlGRYC1JmLaB9mnX2wjF4L129sSKCZ5zyqrlIjghy15Xqgny1s4bT8dc/9EKmC05Lc1NZ31+SJJJOosBNIFEyyJtalNXJ5qda8zt0pQsnbs5A0crD60jZjg1AAB9tbx3yzPKuJp+jwfZ8CtvVmJIMa0HKhMp5DMk9iOKg90B/eVeX1nzMyr4x7Mnd72i6V+M5g3jl+m5WdHaz/E9bWFklS7aDrKWm3gP9ZYfOnP05CyQhkCRt5DeRwnTPjm7ecVeMND6E3cK2C6/XwqaiYldTLHH+yW1MUXY1fjsXk5JhXFWbKZJyq+0QIjNDt4q052OW9o1hT+5uqU5NjOBnaOrY1EOhGC4JZrHC6pCeSYlm3oyJjGqszhY8LhnFXRiiTqSK07vEy3TvjC2JRdReahWgwsbiJhS3iqhW1IkqTJNUxTsWJU5u45GWISnLvhz25O6eGGpiBL/G01mds1xhg0jPvWVHykSnK00N3mpkb4r6nKgTSTWolaS1s4m6AM6bYSW6cA6xOPc69Q86/RvFrRGq675omqTcKlOECDW4QzptkKauPMud+rAnd8cI4oSDPeGocb6hVdV73JPJVDhb03DXA+JJlCqieBfpnh516lpLftJRJ5Lb6BpH5w5QPqHmEXE16PSv33P8xmuUqJimSdA0SbnfEfzHU5yoprUC9xhzDO5p1HvYk3tnV2+5g52M7jqSqTqZvCZKrXgQuI8zI2WosS0q/Pyy/cjDTVC6pyyzRpxdiiphOxMvChzjqV8Go6BxFUTAbvIPEj7cfbdkwUyWrt44JCRDEqShvnC/o5+KUo1q6pxVCBN2aumQlvP8tDUDw57cvQZ5VNLzKkMnnG1UxBlgUSSaOPCL9R1EHu4t9MZte8pH+7Pa5ei0ZdwJnUQlFyYU+JFi1GcGPcft7nr/hh2eIRmS9lGa6gs/VNpFLI8kKYOFnVo8pKW+R1pJyIc9ubt1nC3NTZE70m+imBy8SQaYjkQTt04VsU4UYg4iD4egHK+jv27bU/a+yWKXo9OWSdo7iUrOaZtV67t94+6ouwPHsBr1mUGSs9vddebkw3jixZ7I2ZiiICvJ2D0uvaKauut13ZrnMhc+wuDwVppJyIc9uQMVOk6do+1upG0hh9r0eOjY1FMZ60SIioUtqE0c6d7xOsoJOPXYSVx21vGpv5dOWyZp77gqOVUX7Hhs+QWN6uzq5YJlT5Rj9Fz9oVlD/OPD4Cd8uN1dx41p9AzJYApZOx/ozlW1XlkKHzpon5Z+EnKdHKo3AQuAV6WUs+zvfgR8EOgHXgQ+I6XcLYSYDjwDPGvf3iGlvMRYbT3g1nFK4p3cS3uLWSs6fHedRjVaIXlzti7WvUOIohrIgtg7u3p5autuhBDkPJKM+NUtSnvHWex1M1w5C8BTW3fT74ROLkoeevbVIf7xceHl7to+rSW1bExZCy66c1WtV5bChy6C8uqagI7kfjPwM+AW5bvfAV+XUhaEED8Avg5cYf/2opTyJKO1DIDXAZFFs1trxnjiIIvdQdZ1yvqdHGnXIcV8TrBkwUzP5yb1Soi62Ffogn0yXKkLgDtHwat7Dw7xj4/bnn79kpYAU4uCC1RH+IiCtOdPKLlLKf9gS+Tqdw8oHzuAjxqtVQSo+l+hfFdLnejAS9da7XombSvVg0b9nAY6NvWUE4WAFVEz6IRtVl4JoKcLrpAkbWNgyfZ8Ou/v2nj2lY1GXQuzGlu1KLhA7dXLyy6RZj+Z0Ll/Frhd+XyMEOLPwF7gKinlowaeEQrHFa/a4QT84KR4e3XfIR557jUKReuYtcmgVdVAlvrWeTMm0pgXZcldhwSjqAySGAV1iCRsARjOSaxrXaCqNqrhsZOI3IUQVwIF4Nf2VzuANilljxCiHbhbCDFTSrnX497FwGKAtrb4PreONbzWjJUqOrt6ueBGq2NV9Bclt66tfnybqMGyVBLMUt/aPs2KJe/s0nQWRV2VQRahkHWM1Lr6fV2D4nBdLGoFptqwGg4VscldCHEhlqF1vrSzMEspDwGH7L87hRAvAscD69z3SymXAcsA5syZE5StzRfuvJg54S3NZTnIvZ7leE54odruWV5tGJbuzX1AJkt9a1RJzI9QdRNiRDWqRjk8Fwe6i1At+nYPN5hsw2rYJWKRuxDi/cC/AqdLKQ8o3x8B7JJSFoUQM4DjgE1GauoBZ0I6bk5e1vAsB7nfsfwVnd0V1zXk4D1vPZKHn3tNK35GmlDbEAiVLNwk2Hugv6b0ml5wE6rOieSo5yWyGme6EmAtut4ON5hsw2ro/3VcIZcDZwCThBDdwLewvGNGAb+zrf6Oy+M/AEuFEANACbhESrkrpbr7WsOroTbwUw8BFIqDGXje2TqeJR+cOaSeSesUtyynDfsHSpQYuvtxl+slgdSKXlMXOieSo44b96nQ69Y8l4p3hq4EWKseLMMJacTPcY+HNLUKOt4yF3h8/Uufa1cCK5NWShdeq2E11AZh6qEKI9oHB133TJFi0tOvfsGy4kQwzBJJF7SgE8nPvryPnO13rjNuyoukTfCPPb+zHOvIZBvptn8t9dNwgXs8pd2Gae/2hv0JVTdBVkNtEKYeqsXTr14DWbfcWpDUTS1oXn3iBKMrlmSgL71XmdeteY7Hnt+Zqi1Ft/1roZ+GC0wlRIkicKStVRj25O5GNdQGYYclsn5+mJSpG5tl2+6+IYmsa8UDI+nECOoTdbGWMtiX3l3mZWcdz5ObdyXaKdZKG6uoxTqZRNB4iuKdZCr4mwmMOHKvxna02lvgqM8PI0Z3yrPzT24rJ3SuFQ+MNCdG1qELVKS9VY9D0m8Ezxu/Po/y7uq80on0mDZvjDhyrxaSSOcmpKIozw8jL3WQFkuSKRPG0D4t+2TeEJwDVU3SUg29dtD9js0i6klk3R1JliT9RvC88evzKO8+b0b0SI9p7upHHLkPNynDHXDKL5GCyW1xGHn5kX/WHhhBfelO0mIiPrmKpJMurO5R21637CDEJWndOtWq2ka3bl59HmXMt09LP9JjFIw4ch9uUoZa35L0TqSQxoIVRF5+5J+1+imoL3VUS9UkG7/6hfWlThtHHePqCeS4EVPDDNC1KlDp2peCTg1HGfNpR3qMghFH7nGMi9UkgXkzJpZDxIIVDMs9WdNYsMLe24/8s/TACOrLoN9qgWzmzZjoGZ9dpy/D2jjKGPdyDY6TG1bHAJ3EWystRLEv+Y2VKGO+2vY3FSOO3KM0bi2QQPu0lsFECiVJU+PQyWpaHVIL762DoL4M+i2rA0Wh8IjPbqIvo9gbvFyDLz3z2Dhv44uo7+SX2jEKdBeHKPalfkOCU5YCUBBGHLlDvGD+1VThnHDUOD72d1N9g2GZlgZ03rvaOxoHYeojr9+yOlAE/u3UsanHMz67ib5U7Q1rX9rFhu17WORjuIu7mETpf7936uzqHRLkLSi1Y5T31xVOwtrbnZKwpblJux61jhFJ7rpQB34+J9i+u4/Orl5j3gk6cA9Ux+XQ67mmJK6wCT9cJPsgLJzdysZte3i628rvmkYS4iCVR1AbJ5Xs3NLm8rVbfNP5xVlM4vS/+53ciVXu7Oxm+cVWPfxSO+oiqlDm1d7OvNq2u68iJaHfmYZaEXai4A1N7s7AdxI9L/+Td/jdNMku7PBEGs8Nm/DV3tEETaSwSVbho5/Pkc8LCkWZShJiN8m6VQ1p6V6dhcMJdxF2GjbqYmKi/zs2VSZWccqZNyM4taMOkqq23GOkMS8olqTvQb3hKuy8ockdhiZ69hrMaZJd0EBN87lBEz5rl0cVYS6EYZOswke/WGLW0ePL0rtp1zS1nYQQ5Ry+ajCyNEjAWZxXre/mznVbK4jJD1EkTxP9P2+Gd2KVIBVOUjWQLlat7y4vjMViifNPbmPKhDHl93SPsZqx4UTEG57cIdzzYvvuPhpyQmsSRUXQQK0WyVbT4p/E/RGGtpnp9HUqVMPmzMmHcfMTmzOPa79QI19wVMkzisE2qAy/xCpeKpykaiBddHb1cue6reUQ1/l8rqJu6kG9/gGLyM+eNTkzG45J1MmdYIOQun077+SpvoarpM/3206nTbJBJ0DDnpeGHtI55TdQKJHP67s/qvV2t1la6evcB6niuhkmgU4/RZU8TR0Qq0XHho5NPWVjrgA+2l45n8vGeDsE9uMv7Cz37f0bdqQeFM4k6uRuQx2IDmlt391XscU/2j6GX6166UBHJ+38DmidjvV7TlRpS3sxsH3+y//b0FnsvJ6R1IDphyzcDE0gqvdQVmSb9HBVHLgFhEUuBwZnjF235jkef2FnRd9GDQpXbSPsG4bc/Ro6yHjSkBNDoiJWq563rt1S3ib7EXAY4bp/XzS7NfR0rB+8CMD5PskpRkey8tORBxF11oavatomokAlLB3JM4v3MnW4Kip0BIT2ad7RPaPspGvBCPuGIHe/hvb6vsIgV5Kcd/JUjraNLWl3jl89b127hW/c9VcAHn1+J4AnwYdJXO7frcQiladjV67v1tLhusMB+6Wl89oFBfk2JyGWrL18qqk2i3q9H2F5Qfe9kkim1dz1qDYc9bP7Gr8QHLWmavLDG4Lc/Rra63uvbVuanaKTEvD+DTsq7rl/ww5Pcg8jRq93mzVlfPl0bENesKKzm0IxOA6HVzhgP0m+rPbJicDDImo7xCXMNCTOOGEaTG3Ho0p/JuLWBL1Xkrq5kebuIMyV1nF9do9zrwQ2cfuvFnZ1I5bc1Y7ya2iv77P0FPHamnrV8+xZk8sSu/PZC2F19/pdNThu393H8j9tCZQ2VBIvKOGAgSF1V6+Vis+z+7CIe8E4d87UWP7opvsurl3B1HY8qvSnc70p+4Nu3YJ2EmnMMx1XWscNEqzDbSvXWwnsTapRsuQRP2iRuxDiJmAB8KqUcpb93eHA7cB0YDPwMSllr7AyZv8H8AHgAHChlHK9+ar7w6uDo0Q5NDEBvI5eu+G1NXUOVQnlOkdKD9O569Td63fnO6fOQdKG33Ftv7YsnwDO50BKT3fSCm+OouTWtd6HyXRg0ngaZ2sddE9UiT6q9JeltKjzLJ2dhGnS03GlVc30EljRac0302qUtAz5utCV3G8Gfgbconz3NeBBKeX3hRBfsz9fAZwNHGf/mwv8p/1/ZvDqYL8DJWl0gN/Ra/dz/CbIqvXd9BdKFQT38blt2p4scaEjbfQe6Pc9ru1uS3d54G1wjXriMip0jeluqP0TFp7C6x71xKNz4KhQ0g+WFVX6y1Ja1HlWNfTOQYuO6ubojDOwPOEkQ3eewx1a5C6l/IMQYrrr63OAM+y//wt4GIvczwFukVJKoEMIMUEIMVlKuYMM4GXsy6qjVOOh19FrHTVJNbIduesURl6jGvUngRfhe12jhoEw2W9RjOl+/RMWnsLrHre7qaoKiNKvUYUPHWOhKeiMFccFUwiRSVCuoEWnfZp1OGvJPRvK8W1yULY/LdI4EDackETnfqRC2C8DR9p/Hw1sVa7rtr+rIHchxGJgMUBbmxmJ1M/Yl0VHuZ/txDQBQj0TvKTYWs16oyOxedVP1zBpeoJFMab71atjU3B4Cr93gcETjw6xC4LHQ1LUggueAzeZLl290XjGLL/n+j2j90A/JSltTzE49dhJFQe6ouzsah1GDKpSSimEkOFXVtyzDFgGMGfOnEj3+sHtxjglw0NHQ10orQUrSOfuwD2Iaj3rTdDk8arfsy/vK0/wUY3e/vduI68p+KlWouink+iyK56fz3HG8UcwadwoE6/mCZ1FK0vSUsm0Fk51zpsxeAK6IZ/zPambZJ7VyqKQhNxfcdQtQojJwKv299uAqcp1rfZ3qSMtg5JOZ8V1ofQbRH731oL/bBDc9Vu5vps7ntw6GL97IDzypVOOickRpFqJ4hYYV5et3tvS3FQ+1u8XojcpwuZA/aAXviegVcSdZ2HeOlmSfhJyvxf4R+D79v/3KN9/UQhxG5YhdU9W+vakBiU/dYLOZIj77KiDqCYniwJ3/QSUD0kB5HKV8bu9FgMnal8+Fy0kgh/8VCtRojYm2VE492ZhTwkbh1kKB858uvCU6WzcsTd2EDKT6NgUfALaQdx55te+Ko9EDfURF7qukMuxjKeThBDdwLewSP0OIcTngC7gY/blv8Vyg3wByxXyM4brnAr8SDzKZFCNWc++vE/rOHXQIPKLk1Jt/1k/OPV1jpK3NDexcfseGnKCQkmWB3WQjUFA2fhYKEULiRAEdTvuDkiWFbJamIMWoqzq4PYpzwkSBSEzBd33jzvP/MpXeSRqqI+40PWWucDnp/ke10rg0iSVioskW04/Eo8yGbwGtF+MlzAde9C7mNZLm4DXgSxHBdGQz3G+T0RNL++S2xU1TkkajMGusR0PQtJtdS0szKbr4Ncmbp/yWlEj6r6/GtQsiueRX/nzZkwcEuoj7bYYUSdUk2w5/Ug8ymTQGdC6OvZa16274a7v/Rt2aEfUdL//Re8+hmWPbkJKa3E0IV3qbsf9YEpXrbMwewWzM7kgmBIOgtpE9SkvYQk6UXcKUd47arKPMAcH571K0nKK8HIEiFJ++7QWlp4zazDUR0OObRrnJpJgRJF7ki2nF4m7QxiEreA6A1qXtGtZt+41kdz1PXvWZM8gVTohiW9+YjNSQj4nWLJgppHBn7Q9s1psg3ZAWXlH6RJlUJu4DclRIz5GWUx1oqFGWRzV9wJzh+o+PreNE44ax6r13dz+5BaWr93CinVbWb74lFT6dESRe9Itp7riuhN1IGXo6UKdAZ22zi9tBIVqDUuSESVNngSklIkSFkdxMQ1DVott0A4oix1cFG+PsDZJskOIspgGXRvHkKkevirJwYNOJvq8fVoLK9d3UyhZn/uLViTWOrlnCPeAAb0VPGxARyEZ3cmRpYuV2i79A0OTQqthW93115mwuoe4whaJW9duGVK3JGST1WKruwPSQZxxoePtobanKsxE0U2H1TlKAo+gMRPHkJl01xEGEfLZFEYUuZv04XUfPvELehUHJg2iaSc9cBNES3MTOWFFl3EMRKbcON3eNn7vELZIdHb1suSeDYO+9baLZVJidvdbGjphr0XEaweks2uJMxf8vIr82twp08S8izuWgxbeuIZMHb183PG0cHYrd3YOBuhb6MoGZQojitxN6kXdA8Ypv5ZUJDA0P6aXtBoXfvrfYkmSzwkuevcxkZJCB03CKGQUtkh0bOopxw4BSzIKi1MfFSZ0wkHhcP0+6z430Vzw8CrSlY6TzDt3OVESePiRsduQ2dSoJ5yFxYRPspi1T2th+cXp7wJHFLmb1ot6TbKkMK1CUd9ZCEGxZO6ot5/+19GHjxvTGFlV4TcJo54nCHruvBlWcLP+AStJyHve+mbWPPOKkTC8cerrdS3Ek3bTNsirXkUFRcp12twdjjrJs0zVOQyOITPKLiuob0wsZiZ3734YUeSepREyDimkcfTbrR9cunqjscmho/81NUijTmznuZ1dvVz/0AsV/eC16/rD8695eu7E7Y8o9fW6Ni5BxDXIA0PayQt+8fodeIWjbp9mBQhz8g3Uol0jyjgN65ta9mRTMaLIHbJZEeOSQlrudOo7R5FQdMp1E8Si2a1I+3/Tky/qxI5y0Mur7CT9EaW+ftfGIQjd56rCB+jvEoLi9QcZWx13Tfcp1KhCUBbzNww6XkBB6sVaUd+OOHLPAmlLXUlgenKoErJKEIs0jEBpT+yoqhwv9U2S/ohSXy8VX5JgZGHGPrWvFs5u1W4nR6UVlOwi6Gi9jodNUqRNoDp949UHtRCxVcWIJfc0B0AUUjDpa20SSQ526CxoWQx0E+Rsoj/inihNS0p195VAf5cQ1CZ+v0Ul/SRIc1y5dztRob7voQEr8med3A0jbWKJsjWOEs43K8Rpn6hEmpYKSoUJck7aH+7DbqcffwSPPPdaOX58FtH/3HD31cLZrSyMkAQlqE28fotK+pC+ITtq+Z1dvVxwY4cd510goHxoUdcl03EjdZwO7ly3NbNkQV4YkeSeFbGElZlmPZL4WMepV1QizcroVO3F0u2K+rv/faX8m8mollEQ1FeOp46pXYr6TF3Sd2ctO3fOVG0S1BlXcYQXx1AMVKTI9DqoF7TwfbS9leVrt8SOYWQSI5LcTRNLXCkjLYJL6mMdt15Rdcy1ooJKE05bqjlSVRiNamlDZzy6+8rEbjZOGV5jpmJBLEpuXRuem1YtL2xcxRFe3H2Xt/09RcSDeotmt7JqfXfqQo0ORiS5PltOjAAAIABJREFUmyAW9Rh03MBNaRGcn9+0n6TmvvbSM4/NhHirLVVnAaePV63v5s51WymWJEKARCBtac/kBK+mp5apnah7QYx6LiNsXMURXmZNGU9eQFFCU15w9YdmlXMSRHEvriWhZkSSO4Tn+QxqfHewoSgrd5R6xIV78LY0N5X1hY0NOZZfPDT0qlco45FKvHF3WmGnEv1+c9pS1WtDOieaq+mpZWon6hCgk/qwWDQr5aoEqxPzxnHllEBDziJ21U7iRHLUzQJQK3NrxJK7H6JEJixZRzHJ5QSC+HFl0ojHrUoHqr6wv1BppW+fZuaASVz4pS5MS7KJK9kG5XLV3b25J3UabZ1EpWbC+GxKKnXaalEEQ2/U8kHPv1+d7wLvSKQr7TmWVu7bNPCGI3cdycc9gZIE40rLc0clkpXruyt+UyWMoAMmzu9RvQocKSbsIJMfYYa1RxLyjyvZuu9btb67PKFzhsM6RIUpd9okEqVaB914L7p1An9DbxZjIWzBzMJBIw3EJnchxAnA7cpXM4AlwATgYuA1+/tvSCl/G7uGhqEj+ZiUUIIGho56SKcOi2a3smLdVgaKksa8qDhgFPb8KAuP4y7m7BLCEg342QbCIjomWQzjSrbu+ySDuVyltAKlSWkmKmgU+MUjz3px0ekT0yE5shoLYfO9HAk1Qv/rCEFpH8aKTe5SymeBkwCEEHlgG3AXVkLsa6WUPzZSQ8PQJW5TejO/AaaTPUZ3YLdPa2H54lM83ylogEeVSDo29ZRj24PlMhZ22tHr2WlKSXEXZvd9z768r7wDklip/8aNaczcSKa2RxaJlb0IR6dP0jD0ZjkW/Oa7s/Mt2erZC0+ZrqXD9xKCnHdyxnzah/xMqWXmAy9KKbuESCv0vDnEJe44K63fAAsbuFEHtt87BQ3wqFLuvBkTaWzIlQdtY17EOu0YFtGxIScYKFrSchwpOa4xXb2vY1MPOWEFz8oJGDem0ag6QhfzZmSXWNmPoHXGiTpeDw2UtLMLBQk/23f30ZATifIoJBXS3Pa3Xzz2UqjPu5cQtHJ9d9k2FjUkRFyYIvfzgeXK5y8KIT4NrAMul1L2GnpO1ZBki+gMMDWCYdiEMekjH4f4/cpZfvG8IdtN97FttTyvZ3t9p7qeYicDIaag4GfEddwVw9IlQvaR/4IOB3nFI09jS+8nUOiME2dR7i9atokVnd1aweW8ynaf+j3v5KnGA9XpQh0Huj7vXkLQzn2Hymq+qCEh4iIxuQshmoAPAV+3v/pP4Bqs3ew1wE+Az3rctxhYDNDWlu3xbBW6kyTpFtFrcQiaMCb1/kGIKtm4r3dPRJ1cs254uZ5KoFg0085gbYEPDijSVEgfmm5/dfFyG+fDBAd3PHLA1/U1Sl2CdnL5fI5tu/vo7OodsqNx2kdF+7QWzp0zlVud05kR+s49ptS5ViyWOHrCmKoZMNVxoOvz7haCZk0Zz9X3biir+fL56CEh4sCE5H42sF5K+QqA8z+AEOJGYLXXTVLKZcAygDlz5ui6kBqFn7HKC0klOb/DRHFULWHvlOUBCvd7gbWq9xdKXLfmOS476/jQeiRxPdUJrQCWqkCFo/LRVdMkgTPOHMktJ6ggcR3BQa3LlXf91df1VbcufmGS1QNZt/1pS9n179mX94Uew184u5WVBk5nZr1rCoPa9rohtdV7rn/ohXLKRwF8tL21os3TgglyvwBFJSOEmCyl3GF//AiwwcAzUkEUY1VSSS6LAVuNkKNuac+R3EsSHnt+J09u3hVaDx3XUz9Vi25ohXxOVEywc+dMBdI3asHgOHMkGDeJRx0bbkkoimQUtpA4i43ThwMFS39+x5NbK3LS+i1AJnY7We1a4yBswfcap+7+1QmXbQKJyF0IMRZ4L/B55esfCiFOwhpzm12/1RR0jFXuzqrlAWvKHzeK9O9+L4Dr1jzHY8/vrPALd+oXRwXlt2j57Ya8yirrrW3J09kSZ+G/7Ezu/oESJSzJ3X1aWHdsdHZZ5quGvKDo4fqqW5eghcR9jYCKnLQ54W3oNrlrNLVryhJBUWCrsVglIncp5X5gouu7TyWqUYZwG6sa8oLtip7RtCSc9oA1sTswERzqsrOOr0jH19LcFFpmUNv4kbDf+3qV5ZdHM6i9TJGVW2+rsyvxgjua4nlz2yKHlNUhGq8Fe+X67nJO2qXnzNJegGsZplWYQcJCNRarN9wJVTecSe/EuVj+p8EIdavWd1dYuNOS7EzBhIRgwnDcsamnQrWStMwgEo8SQ8Q9wbzay0TAOJ1nO4hCihWGxpJkSoqGRnd9o0ZiXLm+27ikapKMTS5GTr329Q1EPuyUJkYUucft/LKesVh5/PzOdVsrLNxRjXvVQFIJIYn0HzRh0sqY5Pyd1E3VXf+sQg5EWfiqtTODaJEY8znBis5uCkVzUvyta7dox1XXgU6768xnt7FcYNl3liyYGV6/p++AB5fCnm4Y3wrzl8A7Pxb7ndwYMeSeJGCUI6m5j5/7WbhNPr/WkET6X+mz0zGxo9BV2yRJb1bhtUM2IQeC3A/dqIWdmR/Uum3f3cfyP20x9ozOrl6W3LMh1KAbBWELpe58dhvLnXAVXsHHKvD0HfCbL8NAn/V5z1brMxgj+BFD7nEGrbsDVVUCUBF0P8xolZVxLgvEkf47u3pZ0TkYFtV9ujRNneO8GdHSm+n4eScNGKcLhxS93A/9CD7LnVlU47qj2tJxiYxyxkTHoBulzLCFUnc+hxnLfaXzB5cOEruDgT7r+zq5VyLOltXdgb0H+iuOmNdiWrkwVEs15Ki1YNDVMKvnt0/TT2+m4+eddfuV1YKl+HkDdOBlD0ljN6rTjlHKnjdjIqMac4EG3Tj1DVoodedzoLE8SDrf0+1Znu/3MTC8yV1ZFdvHt/LbMy/nfnGa9sQM68AoUpLugE6TOKqlGurs6mXb7j4a8rly4oWFGfnyOtBNb6bj512NHVfawkFY9EWvcZlkN+rVjupzopStu+iqasH+gaGH6KLkFoiy0PuOGT/p/P4rLCl+z9ah94w3N2+GL7l7rIoznvg6l37wpzBNb1tjWlILIoYsiLcaqiG3e978tx3JpHGjjD8jrI90+7JWdlhupL1r8BsbQapJv7YyEdp3yYKZkfohbNF1qwVLwOMvWIfoliyYyYbte4YYeSHYEJ94ofeTwvt2wcyPwF9urST/xjGW2sYQhi+5B62KESzQOh1oQuI2ZZ0PQlLiivN89b0KJcnv//YqJSmNZayJGvpYRy+cRsyYahzcifJsv7Gh9l//QGmIR0pQUK8koX17D/Qb7QdVLeigJC3j65J7NpQ9n0A/t0Bi+EnnAM8/AB/8ad1bxhNBq2LfLvua5BZoUxJ3FOt8WJwbP6jGuajBeuK+p/peIgX3wTR2I6ZUL9X0kLp17Ra+aZNWk0bwML9Fzd1/7qiH7vhHcfvDa/ybVIG5XTERgmKxckyCZQ/SzS2QGPOXwKqLvX/b021xkkEyd2P4knvQqqgioQU6bDDHsc57Hbjp2NRT1hcmTcoQJ9+j7qT1CscQNWpeFOjsRrKwZZjQSSfJvOUOq/zNu/9K0WYs3eBhXmQapf90Y6x7vUfaaid3+cCQd8rnBOfOmVrhSeV3iC1xHR17oB8M6tb9MHzJff6SSp17EBJYoIPIJY51Hgb1fA35HB9tb2XR7FZampsqpO0gjw834hqqdN8z7H1V0tCNmqeLMFJIW3oOKt9xwXR81IMWsySZt9y/LZzdWiZ2B0nCqur0n26M9SAdfphNKum4cZevMybVe9QDSfmc4L/mdHFq1/XR1SZue6AbhnXrfhi+5O408oNLwyX4BKtkELnEIdIKHWehxPK1ll/zotmt2CkqgKF+4n5IaqjSec8o75uGt0lQmWkbkUPLt4POlf+PWU7Q7+7fBNCUtxJjgBVEzFSkQb+2VusQFGM9TIfvt2NJc4HWHZPq7vkDPMq7/vILEPZhpCgqXi97oIPxU43r1v0wfMkdBnVWV4/3v8bAKuk3OOIYMJ17ysmXsSasBMuXV9G5qxKF7mGLJIaqsElQi54madcpqHzHNz3Mt16nnuphGIC/bN1dPqnqvtdJ9LByfTcCIgcPc2DCIBt0nW7moiy8vHTedd6MieQEFCX8a8MdNAvXKVNdFa+vpkDAP2UXAX14kztYW6AKmVeByFsWacOrpDpQohKpavS8c93Wsu5y0WxLPePlhxu0zU3bUOVV92ockopbpygEFlVXHGVhCatn+7QWliyYWTaSPvC/r/Dws6+yfPEpvvcmaf84KkWdvvfS4fcPWETf0tykFe88Sx9/N4TN7lPETu/CdFS8Gfiw62D4k/uDS/HWOAr4yA2pELt7oERNnOyQr1earSDJ3HHrcm9zs3Tt0104TOhQoxirk2734+iKo7Z9WNv1HuinpByzHygO7gbC7o3a3nGkZd2+V6/b0rOfZY9uoliSXH3vBhBiSECxoHbMyg3Zuc5p/+1yEq1eBK9D0F72wIz07CqGP7n7rqQyFb2WyS2kzmTRcTU0Jamb0n2aKMdEGe6FMSjtXxxdMZi1McybUZlYOSewEoaHIE5bZaFi6+zq5RePvYSzXlk2AuuDTrzzrNyQva67lvP5Yf6X5IsxCLrCHpiOD7sOhj+5+26BpqbyuKz1zu5t7tX3bmCgKLUNrlFgyrXPxAJoooyyHtsuJyjtn3otkEm4Xzfap1mJlX/+yIs8+MwrlCQsXb0x1CU2rhSeRWYwNeCXgyHBtQLuNyFIRVUp9TzxK07fehf5/X2WalcWoxtCU/Zh18HwJ/eMt0BRJ0WULDth7lqdXb0gbPuCECZepwJRFq4wF8GkC6CJMpy+8kr756fzVkPL6pKQSbRPa+HEqRNY88wr2otL3LZKyzaj1mtUY46DSnLyHHDqsZO0Eqe3NDdZyS9iJktXoa1S2vM7ePG7g3wii4N8UmWyjorhT+4mt0CawfOj6J1Vv1nn1Kl7EOr4QHdssuJkF4qWZ02xaF6ijLJwBUlVJqRCk5Jl2+HNNOaDD96ArfO2XRqjkJBpRCXrWjR0g7/zgE6bdnb1snT1RoolqZX8Iqp9xbetMgjFmxUSk7sQYjOwDygCBSnlHCHE4cDtwHSsJNkfk1L2Jn2WL0xsgVIInq/6zRZK1qlTYEj6tiCidAfmUiMvpiFR6ixcXlEg3XUxIRUmLUPn4I060d2kWg1ih3hknbYUHhdBzgM6Lr66yS90VTihGZ18Q/FuTT1zkmmYktzPlFKqpuWvAQ9KKb8vhPia/fkKQ8+qhKkGT2HFnjdjIvmcKG/ziyXJsj+8OCRbUZgvtZo387yTp3L0hDFVk9Dci835J0dP0pxGncJ0/14Hb7ykvWpLwOq7RPXCqhX4uZSGtb36e9Tdi+4J69CMTmNaBmNTuXH3F6A0YP2dQuYk00hLLXMOcIb9938BD5MGuSeRtt2Lgt8p1wShC9qntbD0nFnlASWBzT0HgMoARkGSmnvQeh351oVptzJTSZqT1CuJ7t9L2nMHysoSIyFVY9T0dEGS9qLZrUj7fxMGYreBNyijkyccYncQJvxVWdI3Qe4SeEAIIYGfSymXAUdKKXfYv78MHOm+SQixGFgM0NYWLfphGXGlba9Fwe8gVMKDB05kx6vu/iuq44AQVOgR/VzBOjbpZ84JQlpuZS3NTVz/0Aux65a0Xkl0/7V24jaLk5ppQ/cdgtrePSZ0QyuEqaYcA29gRqe+iNpjP+EvgxypYTBB7u+WUm4TQrwZ+J0Q4m/qj1JKaRM/ru+XAcsA5syZEy/uUdxUVZ6xHzyqYMjrpvdAv2fokSA9omkpzm/SRZWa3a6ZbvtB1i6PYQQdNOFrzRBZa4uNLoLsFn7vENT2aS1yWv2tG21Wvd4LNWCYTUzuUspt9v+vCiHuAk4GXhFCTJZS7hBCTAZeTfocT8Q95htE/mMOt1Zvg9sod8RHsOJIB01e0wNcnXT5fI5tu/u4de2WWOTsEOb1D71gzBc9CqG5F6QkBK2SvzusblonJoPqUkuLjQ6S2C38Fl6Ti5xXiOrAdg2Kwe5GkPAXZJi9dlYmqppE5C6EGAvkpJT77L/fBywF7gX+Efi+/f89SSvqibg+7kGrc9NYuOIlc3XEktAdpY8A3n1cuIudaSnO7ZZ225+2kEuYXMOkL3qUcwN+YYeTwG0o9jomb+rkrYnwDkmeYRJp2C10x4ROfPzI/fXOj1nZ3LyMqmMOt/hBh5gD7XhbB/9PUVWTVHI/ErhLWAdqGoBbpZT/TwjxJHCHEOJzQBeQztIU18c9LEOKYTi6PreLXdihizhSXFiZHZusSIaW/t/yIZZS75CIG6YkzSiEFke9pEN2FeXax+Tdi556zaGBEis1EmS465FkcdB5D2P+3ppIS5UUNibci7E7CQck2P2e/QNvofHsH+iT8PwlcM+lUAx24UxTVZOI3KWUm4ATPb7vAeYnKVsbcXzcg1bnGAbUzq7ewPCrXiSoMwmjSnE6ZbonY1JjbVJJMyrBeJFJlEQXfmRXobZS0rSphDVvxkQaclYcdQms6OyO5L3kFevm7FmTQ9u/s6u3vOMqlILj3egSmimbTrVUSRVtWZTcunYLK12Zx9Q+/XDjH7lo3b/AI9vDhUC30DjGfqdVi63vdATIII5xIwWBEkbCCdW48FudIxpQO7t6uWDZE+XECbc9uZVrPPKfukkwDaORTpm1pNfVJZgwHfv1D71QPjvQPxCc6MKvnd3lOve6fbXPnTOVW9duiXVKWI1fU5Lw6PM7efT5neQEvu+vnnJ2J3hOoq9Ougtx6hbFJ9+kush5TzUvgjs4nNOnPU/8ivnP/9KKFQN66hBHaEzi9aJD7JBaKOA3LrkbClvQsanH3sZbKJb08p+msZ2N4qlQC8Y6HeLV0bGrBusSg5EUdXN+OnCX69VGTpKMOP3mkI0a6wYIfH/1pCYMTfDs94wwEk26C4kq+Zv2/nLec+X6blZ0dlMoegeHa5/WAqt+DsWYnitJ3K393KtVpBgH641L7hA/bIFyOOGisVN4seHDrCqcWv65JMPzn6YhQdeSVK4DncVIZwHoPdBPTlgkmRPWZ92cn1FhwjvnsrOO58nNu+gfKFEiOECZ28vJybkbxwvFfU2SXcjK9d1DTloH3ZvGTtV5z0WzW4ODw8V1mU5yr2+eCQVjDo+mx4+INza5x4FrmzZq/zZ+NOom3jb5ML6/7Z1IO2aFjkSXhgQdt0wdN0DT0CFKnQXA65qw0ANJ6520rIV2ztyZU8YH6tzTXLDj7kI6u3pZ0dkdKd9vmj78FQtmYTDrUxl+nisiB1dPCN61p+Fu7aDgk2fVEIQMSeybBebMmSPXrVtX7WpY8DsyXP7ex71p/FQ6Fz5aM3FJ4hpi/dwAs4Bf/XW9RLwM1g6Z1MpR/jjqjDTHVJzyr3/oBX7ywLOUpKV4+PjcNr77kXek8qwocIKCFUuSUY1K27r15l5oHOOdktPrXr9rVVw7S+8w1PipifKqCiE6pZRzvH6rS+4q/IwnWzrgL7cGD4493VXVZSfRaeq4AaaNoPrrqhnc+nLTEm+SxdO5L4p6IotYM3HGrFsKX2goPEBSOCGbh4xdt31N5Kw47Sr89OhJ3K3DFhRIzVMG6uReCT/jSefNQweDGxknv3WjY1OPr8cIBBOTjhtgXHhJ1FkdOTdJJnGJ1is3q656olZjzdSqbSdQ9aPa166e4F2AH9HGdbeG4AUFUuWNOrmr8OvcMGKvQvJbN/w8RiCcmHTcAOPAi9j8wh34TcwsT1sGPS8u0bp92+/fsKN8tqCluYmOTT2At2dOEFll3S5u1IrHlQrtRSeuHj0MXipdR+Xip95JkTdGBrmbCq0ZNWgQRM+tGBG6emg/j5GOTT1s292n5f8e5gYYVic33IR4/4YdkSI4Zh0CNyiRQ1yDoJMqTkpZ4aoXtNA58COrkRAaOC2Ux/HTd8C1PpyQRmrOMH/4KiTNHv7kHveQgWogdZLgjjkc8k3hR4ZBz6iSEH6T2Ot7N/m0NDdVuALq+nrHrZMX3HU6e9Zknty8SzuCY5ZqibBEDnFUEWqqOCf1raMPDlroVHhJyFm1S7V3B7HhxQn3XGqfGO21T5zmKu9pGJPsmTr+8BknzR7+5B7nkIG78x21S98uyDUORob005OJvBFiD5s8fpPYvdV3gjWp5ON2BTz/5DamGMjgFIVYvAjxhKPGaRNGmu5zbqxc310mdvBO5OBFtLqp4oSkIo5P2EIXBK+FPElMfS/Uwu4g9uLixQnF/sETo14nR/t2JQvilcSXPiUMf3KP06ie8dxtONlWglQ0smiVsWpx7O1VnDgwDgG0NDeVE3+U5KB+3U0+bo8GE5MzKuF6qXt065Gm4c7t17+ic3C85AWeiRyiJjYPi+MTZaFTobaLiZj6Xoi6OzAt5SdaXOISapIgXmnp8RNg+JN7nEYN6/y+XSFxIUTisJ1J4sCoIYRzeCf9CNLXrlrfrZ2+TLfctJCG4c5NHAtnt1IolgDLb/v8k9uGxAbyIpuwPgxrqyTv5txrIqa++z2dRUN3EU9Dyk+keopjO3MQlFnJHUhMzfsQpMf3UgGnbKuDkUDucYwjSTrfK15EjBU/SRyYeTMqQwjr3tvZ1csFN1qTEGDFuq0sX3xKLIIfVjpYF9zEIRi6ywm6R43oGNaHabeVCdWVSujqLkA3YmgaNoBE76XrY+4JaR1AUonXrcZVBT9HuPvgT61/boMpeKuAM0i7N/zJPY4VOm7nj59qLJF2Egk47r0dm3oYsIkdrANLteI77UaaxjyvQzgLZ7dqhUFQIzr+8YWdXHzaDMaNaayqS2LYWAiLde9I3TkhKNlePQOFEr0H+rWiPaZhG0m0Q3RzQlMz9O/Xv99NvEFqXBgU7v5pw1DeuXaW/70DfXDXJZV1Nog3VvgBz63VrsqtUv9+nzjv9jFhv2PFGseIq+194Jbcm/LCV3KvZl2zMObFiaXT2dXL0t9s5C/de8rf5QXcccnf1+QCCeFtqYYSyAE5xegbpd2rPbYDoRsKwA1nTl89gdAgYAi4evfQr3XuTeB5Vw8/AN5bq8YxsPDGykYNO2ygoQbyMrypoUmr5X3QPq2F5RfPC9W5V9tTIgtXP0ddEuVd26e1MOvo8RXkXpLU7O4HwtsyzOirCx31U9UWgLgGVuc+HTWun41P596UsjGNDHLXOcTk5zJ5/xWV1+qoeRrGDJblCtvpdypTN9lC2tCZhNU+9p6lC2TQu3qR0cLZrdyxbms5hn9jPjwiYjWh05ZOhEpTHlVeqKrAEESwjWOs+RyUlS1MjRtk46tijJnY5C6EmArcgpVHVQLLpJT/IYS4GrgYeM2+9BtSyt8mragvgoJ9Pf/AoArGz/ulb5dVhpvgvVZRL6neFbbT71SmbrKFWkCW5OogLNtSWggKe+CXJOS2xack8jjKEkFt6eU1lBaqKjD4EawjmEHwbtwv7Z7qLaOVsi9Agk/BZTKJ5F4ALpdSrhdCjAM6hRC/s3+7Vkr54+TV04CfRL7uJsq6rrB0V7pbIo0DU0GnMvM+iXx1EXVbGzdGe9bujjrZltKC37sGkdFw8xbyq2+WhJtEYEikznF29QN99qFE26HAK1GG6q7ozGsYFPbiqk3CCD6lGDOxyV1KuQPYYf+9TwjxDHC0qYppw3c7E8FQrLsl0jgw5UUWcQ+rqEiS1ixOjPYsCazaaiA/d9Osdy9ZI8t3jCsw+MZo94PbaeLQvsGDiXLQU2zIiVSHgOPmSw1CYDx5ASd+PBVvGSM6dyHEdOBdwFrgVOCLQohPA+uwpPtej3sWA4sB2tra3D/rI5HPulKGF9y6fD/1jut+N1mYIMqoccCvW/Nc1WO066LWiNSRFOMaF5M+N6vdUtbv6DUPVB97dz2GxPvxCGVdgSB/dC+4DZlx86UG4ek7LHdH38iy0lIfp4DE5C6EeBOwErhMSrlXCPGfwDVYovM1wE+Az7rvk1IuA5aB5QoZuwKe+jSNxLQO/LZEXrr8XOPQwGIZhfvVJUBHYncMuDlBWXI3FaPdNLJWAwWhWoa/NJ8bNWxCVvAaq2p9Ojb1UFJctXNh6fzC/NG9oO7GTceHefoOK2CZQux7DkouvLuPmz88hvGjRbLyQ5CI3IUQjVjE/msp5SoAKeUryu83AqsT1TAMXt4tx70vJHOSTf5BR4C9BkppwNLVNY1NLWynn/SmS4BqsKoccOqxk7jsrOPLv1WbPFW437UW6lUtFVFaz40TNiErqGMVGFIf9eBYTogh8X6GzJU4JKnuuk3Hh7n/iiERZu99doC7ny3wm+cG+OQ7m5KVH4Ik3jIC+CXwjJTy35XvJ9v6eICPAPETBMZF2zzrXxzrtgO/gdLXC1e8NPj56TvsQxLJyV4nqUbYJHRL+JeddXyFIbBWUCvSoxvVUhGl9VwvIk/rWVFz3ZbJe6BECUtyV+sTJNB46uKjqmjdu27Tcd491EI3/Xmg/P8n39mU6s4/ieR+KvAp4K9CiKfs774BXCCEOAlLL7IZ+HyiGobBzxXygz9NlHhWaxWPG0veB34SVRRdbC2pOMB8ZqO0Ua32S+u5XkSexrN0Fmuva9Toll66f9V7yfnsq4s/7n2w7pd6FfbataeQUOOsW/bz4EuDapmmvPX/41uLiG/vBfbCVecxf/4y1qxZE/s5XkjiLfMYln7DjfR82r2QhhEE9FZxw8/2moheEwKCVSy1ouIImvC1ZkRVUSvtZwJZLVY6i7XXNZeeeay215eqVvLUxa/SNEyqoUKCUuPFhVMmcOVpo3ii+wAHbIed/mLl/wDNzc1cddVVyZ7pgeF/QtXXCLJ16OGkKNBZxX2eLfd0sz6G54PXRHSHdF25vptV67trTp3hhTBf8VraYVQbQQthUi8a92KVhkpMZ7GOs6AHqZXcuni5p9tT2qyAKqAZ3nl7lXnmMQ3HkXnPAAARA0lEQVSsvqCZBcsHCV5Fc3Mz9913H2eccUa85wVg+JN7kJ4tbkf9/+2da4wc1ZWAv9M9Dw8PmfEQOQY/4kmcF5bDerJhWLFISwIERGQCXssQtNFuAlGykRKRSICiIK+l7G5eu4iIBJnHkk0AG4gDJIq1BGIl+cEYMIudAdbYO2BjY2DxKyTYHs/03R91q7umum69uqqrunU/aeTuqnLV6Vt1zz117rnn+EfzK9Y55/D71w2hka+pIa66Y4ypaafE3cqR+Sw9Y3askDN/R/R3CIFSujOCiOrMrVjIeYYN5nXuOFWbglxypcqVbiDOYO09Jqo4uEtct9LW3Yd4txriTHnLLKS/gloeb/0B5/ybxT1sWDnA3z54jGNTjTeOWbNmsWHDhlwUO3SDcg/L3eBdZRbXjxaWzsAbgWMIjXxH9fGdE6s4UXOiACanaty3ZQ/QHOrlJ6jz+x9kcErCJXVnFJG0KS/rPEuF164wQTcjp3vf1nxqZnz56PAQPRXhxLSi6gn5K12u9BDiDNYjiwbZ8frbxkLkQccHPUP+a41NHGDn1Cr+uedOThJDDWRVi/Xm3VJoouH/Hj4GPX39VKaP0V+F49PQU4HDhwMySWZE5yt392ZtvDZ4v6uc4756mUbzrfc0L0TwhEaqI3t5TQ3x3alVbJK/prdHOOEJ84LmUC8vYUrF/yCnKdRcVGRK0IRYq2x8dm89NroVhdfOMEHXlQbOgP/NR8ZRHuUGUK+iLQ3nQulypbeIfzL0eNTCJOJHif2gcj5Mwfd6bqdHas0H+UMO8yiNZzjnXdtqvHP0OB+ZW+Hbn5jFDY8fY9sbx7j71n/lmmuuSX+9EDpfuUMjoX7QjXLzRHgJe/UyjdqmFWY6NFKA13cfYsnEAe7XHXDjs3t58JlXmZpWgaFeXpIolaTujCIjU7IeWLbuPsSDz7xaHzSr1fQKr51hgv4lddNaubnXBZiadgasqakaa3/xPGedOZsrl8/PRREXNWk8NnFgRiFyRaMGcCs0Bqwl7FHvZ/jJm6JDGk0LIJdclF4QQyDG7P4TfPfCfr462kdFhAsWn8wtY5P8/rU96a8VQXcodzBHt5gWMnmVuNfHLpVgRe4W9PDjGeX9HWZk0WC9wo8b6mXyNeYZPVJkZErWA4tXOQiwciR9VsakYYKtuLauXD6fh3Sq4GoFKpVK04phb8z3tr1H2Lb3SL0UYpyKSJ2A637y3sOgGsBpaPS/98HpJ0e7Ypetctyt3iSDKMf9unA0nd/dEIjx8Inr8A7x1Yrwtb/q52vgFPTIYUFkd1ViCgprMln0bjhUaFIfTe+Ak9zHv+o1YQWVKCu2EycJ41z3M3c2fM1ZWO5JzucWSjHlK4/bLt6l8tWKE6HhL6IdR/awLJ1uTqDf72xMCgrw9Ys/EKjcS139KAR3AVIcn3vutFBZLZPreElRkSmsElN3Kfcgtj8AD3+pkRkOnInQy3/oNGJUo3tTg8YpChKA2wn3HT7K+qf2UFNOebbrLwrutGWkFUWStRJKopCvWvckk7qwRl9PhfuvTadIbtu8i+/914667dVTETZ8IXlx8SjilkLMex4lb0OjNPnwjWXwDGXz0rL9Adg403oPJOGgYsvsiZi/R82Me4txpMjpPCP1brVCT0WYrqkm90hWnSmPTtmqIsnavxv3fGMTB+oVk6A1t9Do8BBVjzuhpvIpLj6yKF4pxCB3l7s9i2cozyRm3nNf2UKBkEye9TwmVYNYtsoc9OElwyRi3a/cn1jblLyH6UmnoZ9YG16lCVqOe/V2wunpGqs/tpAzThtoeh3PojPl1SnTpEUog8tgdHiI3qrULfdW5htGFg2ydsXSGe6EoHNl8bvjRod45wsGT+rL7N5nMU+Sd9qJzJ5106TqnOHMckbVmb0gfS3WFHS/cg8bCU1pfJOcIwJ/Jwzy+2ZlheUVFRM3LULWg1WrjCwa5P7rzg31uSfh6nMWhhZeyXOVqf86YxMzc7Fnee9bnYBvR9qJzH6vaVL15d82jsmqaEcrtVhT0P3KPSpT3Iw0vobjpAK/vL5Rk9UdySHSBx8npjgrKyzPqBh/EWV/WgRv54rT8dpl2Xut4K27D3Hb5l25WdVRq0xbmYx1CVOc/mco7W8NWjiX5FztSDuR6bO+8zEifeFZFO2ol/vTkXcDc5x9SbLVJqC7lHvQhGec6uNuGl9T5Iyanplt7sirsPELQG3mNsPoHvWa7X/g01oleSxOMRVRDutcUR2vCMs+z1WtLqbfPTZxoL7oaqqmuPmRcT7w7lNTXd/0bHjv/eBJfaz95fMt/Vb3nGnaLc+0E95z3HzZWWwa388lS+e1dr6kZTaTBlb49Yqadqx0fw3XjOke5R6UNmDjtc7o+JGrtdVtsMxdP5fb0KFlsVwCVsDp0X3r7AtTJQ3zHpvWKsl68jKOMglK0xo2yLTySp3W4s/D1+smrXKtcNPv9k/GTtfST8aGKU73voS9VSUlTbu1amTEzQvvDmBPv3Iw9WAJxC/VOXt+umRjeWWujaB7lLupxNbRg058+qdudb5HpfFdtkqHLKVDHdnbsoWYhwWeljjKJIiwfWlfqVuxvvPw9dZUsxUe9LtHFg3y+fMWc/vvJoDWVmWmcfO14rJIe660RkbcvPC3PP5SJikogHhv966eSKOo88hhE4PuUe5hDeU2vhs/GvRKFWeVagze7p/L5B9bt5qytsCT4Lecsh5o0p4zrfWdZTHo0eEhKiL1fOK1mFb4qQO99cq+FVpblZnUzdfKPWu3oRF1j73K323L1AOYt88PDELPQMP/veSi5jm2MMMvSP/U87ob/Pk5lddz6R7lHvVq5TZ+UKx6kE8sFcJb59xI3+ZyFqGIg8lyyrpTJzmnq5wHT+ozWpGmV/ms/fszQiJrir7eePd4dHiI/t58ngtTNtGs5lz8WTPzDH+NelPwKv+KNGoEJ76ev88fPehY5256bxNx4+KjVr7nWF7PpXuUe9SrVdgoaXLpSNVJEzpjJA8ZQD76Dwxf8Pfc+96GMsoyG2I7KDLJWBB+5RxkfYcp8Dx+T1RIZBB5WcDtXHB082VnGSdqs5Ijqp38yj9UsYdNfKb1gwfpmUovTP55Zo4Yk06B4BJ/OdA9yt1tqE03NC9KiholjZkga81LkI3LlYHLnDrh7sNWhljvpHg7T7VaYd/ho2zdfagw2f3K+dA7k00pG8IU+OjwED3Vxu/JymL2W8ZROWy8/yeLkEyXPAdj/7k3je9vKfw1LlFzObEGyaiJz7R+cH9isIFBmPxTQ+f4U4w3IdnmrAmhkteJReSTIrJDRHaJyI15XWcGy1Y5IY1X3NGIIQXHlxaGyaoP2m48dsGMr6aFSZ3AFcvn8/EPzQWlWP/UHj5z5xhbdx8qRBZ3sKmGpEuOPMbNn5RxHiVXSd+3ZQ9XrXuS+7bs4d4te7jqDnN7uRbu9x/bkUm7xmmfrM59ydJ5xmvlKYefkUWDkbVXQy1zSNbn/Sxb5SjoNYed9TH+BZBuLHvU+d3KbmtOc/7d/kD0tROQi+UuIlXgNuBCYC/wtIg8qpR6IY/rBeLNCXP0YHi4Upxi2AmPzXNBUV74w/xqShW+ECmOpRZ2zNiEkyJY0VoIoh9/W3lzlIdZrllb2nlOeAad2+SOKlOEFxBtmSfp82mu48aym86fR/1WH3m5ZT4G7FJKTQCIyHpgBdAe5Z7UnxanGHbCY0v3sMfAq3hQikpFEJqTnEF7FyLFmRw0HZPXIDujrVBUBdwcZWHXyUOePCOr/OdOG/7abo6ffAb9f97XvMO/piVFltem85lSBtdTjmfo809AXsr9TMD7i/cC53gPEJHrgOsAFi5MtxTbSBp/WpKMjzGPLdPDHge/4gkLHSzbxKuJvAbZoLYaf+1IZA6bThz0O5HfLvgi5724dmY91aA1La0q0rA3gLDztyH2vbAJVaXUOmAdOPncMz15u9J4dhlJFE8nuZ3yCuVMq6Q7bdDvRIbOvYZvvvA617OeeRzgxCln0H/xmuwjVNK+AbRBR+VSrENEzgXWKKUu1t9vAlBK/UvQ8ZkX6wiKMU1R5cQSThnS+losJkr9fGako4oo1vE0sEREFgP7gNXA1Tldq5ms/GmWUKwFaikzpX4+26CjciuzJyKXArcAVeBupdS3TMfmWmbPYrFYupRCyuwppX4F/Cqv81ssFovFTG6LmCwWi8VSHFa5WywWSxdilbvFYrF0IVa5WywWSxeSW7RMIiFE/g/YnfK/nw68laE4edIpslo5s6dTZLVyZk+esi5SSr0raEcplHsriMgzplCgstEpslo5s6dTZLVyZk9Rslq3jMVisXQhVrlbLBZLF9INyn1d0QIkoFNktXJmT6fIauXMnkJk7Xifu8VisVia6QbL3WKxWCw+rHK3WCyWLqSjlXshRbhjIiKviMgfROQ5EXlGb5sjIr8WkZ3630LykYrI3SLypoiMe7YFyiYOt+o23i4iywuWc42I7NPt+pzOPuruu0nLuUNELm6jnAtEZLOIvCAiz4vIV/T2UrVpiJxlbNNZIvKUiGzTsv6T3r5YRLZomTaISJ/e3q+/79L731OwnPeIyMueNj1bb2/fvVdKdeQfTirh/wWGgT5gG/DhouXyyPcKcLpv23eAG/XnG4FvFyTb+cByYDxKNuBSYBMgwCiwpWA51wBfDzj2w/oZ6AcW62ej2iY55wHL9edTgZe0PKVq0xA5y9imApyiP/cCW3RbPQCs1ttvB76oP38JuF1/Xg1sKFjOe4CVAce37d53suVeL8KtlJoE3CLcZWYF8GP9+cfA5UUIoZT6HXDQt9kk2wrgP5XDGHCaiMwrUE4TK4D1SqnjSqmXgV04z0juKKX2K6We1Z/fBl7EqSNcqjYNkdNEkW2qlFJ/0l979Z8CLgAe0tv9beq29UPAx0VECpTTRNvufScr96Ai3GEPartRwGMislWcYuAAc5VS+/Xn14G5xYgWiEm2Mrbzl/Ur7d0e11Yp5NTugL/AseBK26Y+OaGEbSoiVRF5DngT+DXOm8NhpdRUgDx1WfX+I0BbCvv65VRKuW36Ld2m/y4i/X45Nbm1aScr97JznlJqOXAJ8I8icr53p3Le0UoZh1pm2YAfAe8Fzgb2A98vVpwGInIK8DPgq0qpP3r3lalNA+QsZZsqpaaVUmcD83HeGD5YsEiB+OUUkaXATTjy/iUwB7ih3XJ1snLfByzwfJ+vt5UCpdQ+/e+bwM9xHs433Fcw/e+bxUnYhEm2UrWzUuoN3ZlqwB003ASFyikivTgK816l1Ea9uXRtGiRnWdvURSl1GNgMnIvjxnAryHnlqcuq988GDhQk5ye1C0wppY4D/0EBbdrJyr1ehFvPmK8GHi1YJgBE5GQROdX9DFwEjOPI91l92GeBR4qRMBCTbI8Cf6dn+UeBIx5XQ9vx+Sc/jdOu4Mi5WkdNLAaWAE+1SSYB7gJeVEr9m2dXqdrUJGdJ2/RdInKa/jwAXIgzR7AZWKkP87ep29Yrgd/ot6Ui5Pwfz6AuOPMC3jZtz73Pa6a2HX84M88v4fjivlG0PB65hnGiDLYBz7uy4fgAnwB2Ao8DcwqS736c1+8TOD6/z5lkw5nVv0238R+AjxYs50+0HNtxOso8z/Hf0HLuAC5po5zn4bhctgPP6b9Ly9amIXKWsU2XAf+tZRoHbtbbh3EGmF3Ag0C/3j5Lf9+l9w8XLOdvdJuOAz+lEVHTtntv0w9YLBZLF9LJbhmLxWKxGLDK3WKxWLoQq9wtFoulC7HK3WKxWLoQq9wtFoulC7HK3WKxWLoQq9wtFoulC/l/bD1B3hzUJkkAAAAASUVORK5CYII=\n"
},
"metadata": {
"needs_background": "light"
}
}
]
},
{
"cell_type": "code",
"source": [
""
],
"metadata": {
"id": "YBf6X09TpAjG"
},
"execution_count": null,
"outputs": []
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment