{ "cells": [ { "cell_type": "code", "execution_count": 21, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import matplotlib.dates as mpl_dates" ] }, { "cell_type": "code", "execution_count": 31, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAGdCAYAAAAR5XdZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB1a0lEQVR4nO3de3xcdZ0//teZSya3mdyvTdqmLfdCgVLkUmhxoQWxwuJXWfnuAi76xbWtsuj+dru6C+u6W9wVZVcQ70VUKl5QURAsAi2KlbaktBQovSRN0txvk8lMMtfz++PM58xMMveZk7mc1/PxyAOSTDInnSTzzvv2kWRZlkFERESkU4ZcXwARERFRLjEYIiIiIl1jMERERES6xmCIiIiIdI3BEBEREekagyEiIiLSNQZDREREpGsMhoiIiEjXTLm+gGQEAgH09/fDarVCkqRcXw4RERElQZZlOBwOtLa2wmDI3/xLQQRD/f39aG9vz/VlEBERURp6e3vR1taW68uIqSCCIavVCkD5x7TZbDm+GiIiIkrG1NQU2tvb1efxfFUQwZAojdlsNgZDREREBSbfW1zyt4BHREREtAAYDBEREZGuMRgiIiIiXWMwRERERLrGYIiIiIh0jcEQERER6RqDISIiItI1BkNERESkawyGiIiISNcYDBEREZGuMRgiIiIiXWMwRERERLqWUjC0fft2rFmzBlarFY2Njbj55ptx9OjRuB/zhz/8AVdeeSXq6upQVlaGs88+G1/96lczumgiItInf0DGd145icN99lxfChWRlIKh3bt3Y/Pmzdi7dy927doFn8+HDRs2wOl0xvyYiooKbNmyBXv27MHbb7+Nz3/+8/j85z+Pb33rWxlfPBER6cueYyP44jNv4/O/ejPXl0JFRJJlWU73g0dGRtDY2Ijdu3fj6quvTvrjbrnlFlRUVOAHP/hBUrefmppCVVUV7HY7bDZbupdLREQF7usvH8d/PXcU1lITDt+/MdeXQwkUyvN3Rj1DdruSpqytrU36Yzo7O/Hqq69i3bp1MW/jdrsxNTUV8UJERHR8eBoA4Jj1we7y5vhqqFikHQzJsox7770Xa9euxcqVKxPevq2tDRaLBZdccgk2b96Mj33sYzFvu337dlRVVakv7e3t6V4mEREVkRPBYAgAeidcObwSKiZpB0NbtmzBoUOHsHPnzqRu/8orr2D//v34xje+gYceeijux23btg12u1196e3tTfcyiYioSMiyjBMjoR7VPgZDlCWmdD5o69atePrpp7Fnzx60tbUl9TEdHR0AgPPPPx9DQ0O4//778ZGPfCTqbS0WCywWSzqXRkRERWpwahbTbp/6eu/4TA6vhopJSsGQLMvYunUrfvGLX+Dll19WA5xUybIMt9ud1scSEZE+HRuajnidZTLKlpSCoc2bN+OJJ57Ar371K1itVgwODgIAqqqqUFZWBkApcZ0+fRqPP/44AOCRRx7B4sWLcfbZZwNQ9g59+ctfxtatW7P5dRARUZE7PhwZDPVNMDNE2ZFSMPToo48CANavXx/x9h07duDOO+8EAAwMDKCnp0d9XyAQwLZt29DV1QWTyYTly5fjgQcewN13353ZlRMRka4cH1GCoVVtVXijz47ecWaGKDtSLpMl8thjj0W8vnXrVmaBiIgoYyIztP6sRrzRZ0ffxAxkWYYkSTm+Mip0PJuMiIgKghirv/rMBkgSMOP1Y3Tak+OromLAYIiIiPLeuNODMacS+JzTYkWzrRRAdsbrZVnG/7xwDL86eDrjz0WFicEQERHlPVEiW1RdhvISE9prygEAvVloou4adeKrL7yL+58+kvHnosLEYIiIiPKeCIaWN1YCANpqlQnmbDRRDzuUVS9Ts76kemOp+DAYIiKivCeCoTOCwZDIDGWjTDYRLL/5AzLcvkDGn48KD4MhIiLKe2KsfoXIDNUomaFs7BoSvUgAMOPxZ/z5qPAwGCIiorwnJslEMNReG+wZykKZbDwsGHJ5GQzpEYMhIiLKa063D6cnlQzQiobIYOj05Az8gcz6fCKCobCzz0g/GAwREVFeOxEskdVVlKCmogQA0GwrhckgweuXMTQ1m9HnjwiGWCbTJQZDRESU147PKZEBgNEgobU6O31DEy4GQ3rHYIiIiPJatGAIANqzNF4/Nh0eDLFMpkcMhoiIKK/FDIbUxYuZBUMskxGDISIiymtzx+qF0ERZ+mUyWZYx7uJovd4xGCIiorzl8QVwakzJ/MwNhkK7htLPDDk9fnjCFi06WSbTJQZDRESUt7rHnPAHZFRaTOrhrEKbuoU6/czQhDPy1HuWyfSJwRAREeWt8DPJJEmKeJ9ooB6wz8DrT+8YjbE5wRDLZPrEYIiIiPKW2jzdUDnvfQ2VFlhMBgRkoH8yvezQuNMd8TrLZPrEYIiIiPJWrEkyAJAkKeMzysad3ojXmRnSJwZDRESUt+IFQ0DmZ5TNzQyxZ0ifGAwREVFe8gdk9SiOM2IFQxnuGhKZoYoSIwAuXdQrBkNERJSXTk/MwO0LoMRkUDNAc4kyWbq7hkRmSEymMTOkTwyGiIgoLx0fcQAAltVXwGiQot5GBEnp7hoS26cXBYMqJ4MhXWIwREREeSl8rD6WUJks3cyQEgyJDNMMy2S6xGCIiIjyUryxekHsGhpxuDHrTT2rMzcYYplMnxgMERFRXhLB0BlNsYOhqjIzKi0mAOmVykLBEHuG9IzBEBER5R1ZlhOO1QORu4ZSLZV5/QFMzSplsVBmiGUyPWIwREREeWfE4cbUrA8GCeior4h7W7WJOsVdQ+JcMoMEtFQpwdCsNwB/QE7jiqmQMRgiIqK8I7JCi2vLYTEZ49423SbqcZcSDFWXl6ilNgCYSaP3iAobgyEiIso7x0cSl8iE0K6h1DJD49NKMFRbUYJSswHiHFiWyvSHwRAREeWdUL+QNeFtQ7uG0ssM1ZaXQJIklJuVDBTPJ9MfBkNERJR3kmmeFsR4fapHcohJstqKEgBAebBU5nQzGNIbBkNERJR3UgmGxFj8pMsLx6w3wa1DxkSZrDIYDAXPJ5vxskymNwyGiIgor9hnvBh2KGeGLW+IP0kGAJUWE2rKzQBSO6NsIqxMBgBlZnFYKzNDesNgiIiI8orICjXbSmEtNSf1MemcUTY2p0xWwTKZbjEYIiKivCICmiV10U+qjyad8fqJuT1DLJPpFoMhIiLKK2IrdE2wfJWMttrUx+vnNlCzTKZfDIaIiCivTM0oTdC2MlOCW4aIJupslMlcLJPpDoMhIiLKKyIYqipLrl8IANqDixeT3TUky/K8MllZCTNDesVgiIiI8spUcDzelmTzNBBqoO4dd0GWE58tNjXrgy94BpnaM6SWybLTMxQIyPD6A1n5XKQtBkNERJRXpmaUYMSWQmZoUbWSGXJ6/JhwJd41JLJC5SVGlAaDILF0MRuZoUBAxqaH/4CND+3BLM86y3sMhoiIKK+omaEUeoZKzUY0Wi0AkmuintsvBISmybIRDI1Ou3GkfwonR5x4o3cy489H2mIwREREeUVtoE6hTAakdkbZ3EkyIDwYyrxMFj7iv697POPPR9piMERERHlFjNanUiYDQk3UyZxRNrd5GgDKS7JXJgufavtzF4OhfMdgiIiI8kqmmaFMy2TZOLU+PDv1+qkJ+NhIndcYDBERUd6QZTmtniEAaFMzQ4nLZHPPJQNCo/XObJTJwgIyp8ePtwccGX9O0g6DISIiyhuz3gC8fmXkPeXMUE0KmaE5J9YDQEWwTJbNzJDRIAEA/tw1lvHnJO0wGCIiorxhD5bIjAZJLVsla1lDJQCgZ9wFjy9+WWrc6QYQmRnK5jSZ6Bm6+ox6AGyizncMhoiIKG+IEllVmRmSJKX0sU02CyotJvgDMk6NOePedjy4iyi8ZyhbZTJ/QMbpSSUz9JcXtwEA9nVPJLUMknKDwRAREeWNUPN0av1CACBJEpY3VAAAjg9Px72tyAzVaVAmG3bMwuuXYTJI2HBuEywmA8adHpwYiX9NlDsMhoiIKG+EmqdT6xcSljcqpbJEwdCEU7mfmigN1L6AnLDMFk/vuJIVaqkuRanZiIsWVwMAXuuaSPtzkrYYDBERUd5Qj+JIsXlaWBEMhuJlYdw+P6bdyv3UVVjUt4f3KGWyeFH0C4mG7kuX1gJg31A+SykY2r59O9asWQOr1YrGxkbcfPPNOHr0aNyPeeqpp3DdddehoaEBNpsNl19+OZ5//vmMLpqIiIpTumP1wvJgE/XxOMGQ2D5tNEiwhpXjzEYDSozK02ImTdRikkyM+l/aUQcAeI3LF/NWSsHQ7t27sXnzZuzduxe7du2Cz+fDhg0b4HTGblTbs2cPrrvuOjz77LM4cOAArrnmGmzatAmdnZ0ZXzwRERWXdBcuCmpmaNiJQCB6w7IIhmrKS2AwRDZpl2VhokyM9rcFM0MXLa6G0SDh9OSM2lhN+SWl0Pu5556LeH3Hjh1obGzEgQMHcPXVV0f9mIceeiji9f/8z//Er371K/z617/GRRddlNrVEhFRUUv3KA5hcW05TAYJM14/BqZm1dPsw4lgqC5skkwoLzHCPuPNsEymBDzttcp9V1hMWNlqwxt9duzrGseiixal/blJGxn1DNntdgBAbW1t0h8TCATgcDhS+hgiItKHTKbJAKXUtbQ+/kSZmhmqmB9wZWPXkDgbTWSGAODSDuU57zX2DeWltIMhWZZx7733Yu3atVi5cmXSH/fggw/C6XTiwx/+cMzbuN1uTE1NRbwQEVHxy3SaDABWNIhSWfxgKLx5WijPcLze5w9gwD4LINRADQBrgk3U7BvKT2kHQ1u2bMGhQ4ewc+fOpD9m586duP/++/Hkk0+isbEx5u22b9+Oqqoq9aW9vT3dyyQiogIiNlBXZRAMLW8MZoZiNFFPJJEZSnfx4uDULPwBGWajhEZrKNgSwdDx4Wk1GKP8kVYwtHXrVjz99NN46aWX0NbWltTHPPnkk7jrrrvwk5/8BNdee23c227btg12u1196e3tTecyiYiowGQ6Wg+EN1FHD4ZCJ9ZHywxlViYTO4YWVZdFNGfXVJTgzCblujhin39SCoZkWcaWLVvw1FNP4cUXX0RHR0dSH7dz507ceeedeOKJJ3DjjTcmvL3FYoHNZot4ISKi4pfpaD0ArGiwAoi9ayh+A7Vyvy53epkhdcdQbfm897FUlr9SCoY2b96MH/7wh3jiiSdgtVoxODiIwcFBzMyERgW3bduG22+/XX19586duP322/Hggw/isssuUz9GNF8TEREJmY7WA8Cy4JEco9MeTLrml6RCDdTRp8kAwOVNMzM0Z8dQONFEzcxQ/kkpGHr00Udht9uxfv16tLS0qC9PPvmkepuBgQH09PSor3/zm9+Ez+fD5s2bIz7m05/+dPa+CiIiKniyLGc8Wg8oo+ytVaUAomeHEo3WA+k3UPdFmSQTRGboSP+UugGb8kNKechkTtx97LHHIl5/+eWXU7kLIiLSKZfHD39wUWImmSFAOaOs3z6L48PTWL0kcpXLhCu0dHGusmCZzOlONxiKnRlqrS5DW00Z+iZm8PqpCVx9ZkNa90HZx7PJiIgoL4h+IbNRQqk5s6cncSzHiZHIExICARkTLuV+wk+sFypEZsibZs/QeOzMEMBzyvIVgyEiIsoL4ZNkkiQluHV8K2KcXm+f8arZp+iZofSnyTy+AAangjuGaudnhoCw5Ytsos4rDIaIiCgvZGPhoqAe2DonGBoPlsisFhNKTPOfAsszKJMN2mcRkAGLyYCGyvlj+wCwJhgMHeydhNuX/pZryi4GQ0RElBfsrsyO4ggnMkO9Ey7Mhk2Giebp2iglMgCosKRfJgsdw1EWM7O1rL4C9ZUlcPsCONzHqep8wWCIiIjyQjYzQ/WVJbCVmiDLQNdoqG8o/MT6aMrM6ZfJ4k2SCZIkhfYNsW8obzAYIiKivKDuGMpCMCRJUmgTddh4fbyxeiB86WLqwZDYPh1tkiycCIb2sW8obzAYIiKivKDuGMpwrF6I1kQdb+EiAJRbxNLF1Mtk8bZPhxNN1Pu7J9RmbsotBkNERJQXQpmhzHuGgOhN1IkzQ+kvXYy3YyjcOS02WEwGONw+nJ6YiXtbWhgMhoiIKC+oPUNZzgyF7xpSG6hjBUPm9KfJRAN1e5yeIQAwGiTUB6fNxqMcF0ILj8EQERHlBXXPUBZ6hoBQMHRyZFotRyVbJpvx+hFIoYTl9vkxNOUGkDgzpNy/8jVOOBkM5QMGQ0RElBdCmaHslMnaaspRYjLA7Qugf1IpRyVbJgOA2RT2AIlyV5nZGDPrFE5Ms00wM5QXGAwREVFeyOZoPaCUo5bVKyfYi76hRJmhUlMoGEqlVCb6hdprY+8YCieCoXFmhvICgyEiIsoL9pns9gwB85uoE2WGDAYprSbqUPN0/H4hQWSPmBnKDwyGiIgoL4ieoaosZYYA5fR6QNk1NOPxYya4jTpeKUsEQ6mM14dvn05GKDPkTfo+SDsMhoiIKOcCARmO2eyO1gORu4bE5JbZKKHSEvs+xGGtaZXJks4MsYE6nzAYIiKinHN6fBDDW9ktkyk9QydGpjE+HRqrj9fXUxHcQp1Kmax3PMXMUDAztVCj9ZMuDw6cmliQ+ypEDIaIiCjnxPbpEpMBpWZjglsnb1l9JSQJmHB5cXzEASD2uWSCmhnyJF8mCzVQJ5kZEtNkC5QZ+sxP3sAHH30VnT0MiKJhMERERDk3pUHzNKAENouqlWzNa8GzwOpinFgvpJoZmvX6MTqd/I4hIJQZWqgG6qNDSiB4dNCxIPdXaBgMERFRzmX7KI5wom9IBEO1FZa4txeZoWRPrhdnklVaTEk3f4emybyQZW3PJ5NlGcMOJVgbnJrV9L4KFYMhIiLKuWwf0hpuRUPksRy15fHvQ50mS7JM1ht2JlkyO4YAoDp4Df6ArH7tWpma9cHjCwAABu0MhqJhMERERDkXygxlPxgS4/VCosxQebBMlnRmSG2eTq5fCAAsJiMqgkGX1n1DI45QAMTMUHQMhoiIKOdCCxe1K5MJYqw9lvKUy2Sh7dOpWKiJsuHgmWkAM0OxMBgiIqKcy/ZRHOFEmUxInBlKrUyW6vbp0HUszETZyHRYMMTMUFQMhoiIKOe02D4t1FSURGycTnSQaqplslS3T6vXtUDnk4VnhiZdXsx6k9+fpBcMhoiIKOdCJ9ZnPxgCIrNDiYOh1M4mS3X79Nzr0Hq8PjwzBLBUFg2DISIiyjktR+sBYHljhfr/iYKhVJYuOt0+NbPTlmrP0AKdTzY8pzTGUtl8DIaIiCjntM4MLQ/LDFUnGK2vSKFMJrJCVWXmlK99oc4nEzuGBGaG5mMwREREOSd6hrRooAZCE2VVZWaYjfGf+lIpk/Wl2S8ELNwW6pFgMNRoVRrHmRmaj8EQERHlXCgzpE2ZbPWSGiyrr8D7zm9JeNtUymSpHtAaTpTJtA6GRGbogrYqAMwMRaPNdx0REVEKtFy6CADWUjN+/5l1SW2ITuVssnSbp4GFmSab9frVHU4XtFXjhbeHGQxFwcwQERHlVCAgw+HW7jgOIdmjMlI5m6wv7CiOVIWfT6YVcYBsidGAM5usAIABlsnmYTBEREQ55XD7IM4q1WqaLBWpLF0UO4baa9PIDAUbqCddHvgD2hzWKkpkDVYLWqpKAQBDzAzNw2CIiIhySpTISs0GWEzGHF9NqEzm9cvqAaexpLt9GgiVyQJy6N8g28TCxQarBc3BYGjYMQufP/7XpTcMhoiIKKe0HqtPlSiTAfH7hhyzXrUfZ1EaZTKz0QBrsGFcq/PJxMLFRqsF9ZUWGA0SAjIwOq1t03ahYTBEREQ5pfVYfapKTAaYDEp/kcsbu1QmGpFtpSZUWtIr72l9PtlIsD+owaoEQhyvj47BEBER5ZTWY/XpSObkehFQiPJTOrSeKAtlhpRrFNc6aJ/R5P4KFYMhIiLKKa3H6tOhHtbqjh0MDdhFMJR6iUyoKRdN1Nr3DAFAs00EQ8wMhWMwREREOTU1q/1YfaqSmSgTU1nNNkva9yO2UGvVMzQ8Z/u0yAxxvD4SgyEiIsoprQ9pTUe5JRgMeeNkhqYyzwzVlmvcMySCIVtkZojj9ZEYDBERUU7l2zQZAJSbE5fJQpmhDHqGKrTrGQoEZHXpYsPczBCDoQgMhoiIKKfs+dgzZElcJhMBRUsGDdS1KR7WOjbtxjuDU0nddtzlgS8gQ5KA+so5mSGWySIwGCIiopxSR+vzKTMkTq6PUyYTAUVTJpmhFKfJ/vb7+/G+/3kFPWOuhLcVJbLa8hKYjcrTfUuwpDdgn4Usa7P1uhAxGCIiopwSZbKqPMoMlQXLZM4YZTK3z4+xYACTyWh9KueTBQIy3uq3IyADb/bbE94+/CgOQfQOuX0BNSNHDIaIiCjH8rGBuiJYJpuJUSYTI+slJoM6Hp+O2uD5ZMlkhoYdbnj9Sjana9SZ+PZhCxeFUrNRDcDYNxTCYIiIiHLKkYej9YlOrlcXLtpKIUlS2vdTHSyTTc16E54X1jcRKo0lEwzNXbgoiLIet1CHMBgiIqKcysuli6JMFiMYCi1cTL9EBgDVwa9ZlpGwbCUOhQWSzQzNL5MBob1IHK8PYTBEREQ54w/IcLhFZqhwymTZGKsHAJPRoPZKJZooC88MdSeTGZqzcFFoDmuiJgWDISIiypnp2VCwYS2gMlk2xuqFWnXXUPKZoTGnJ2Emae7CRYHj9fMxGCIiopwRk2RlZiNKTPnzlJTooNZsjNULogE7URN1eDAEJM4ODTuCDdSVkcFQCxcvzpM/33lERKQ79jycJAPCDmqNUSYbCJ76ns3MUKIyWW+wTGYJBo3dY/GDoVBmaE4DdRUzQ3MxGCIiopxRm6fzqEQGJJMZUgKNpiwEQ8ksXvQHZPRPKgHYpR21AICTI7GDIafbpzZ/z+0ZYmZoPgZDRESUM/m4cBGIHwwFArKaVclqZihOMDTsmIXXL8NkkHDZsjoA8TNDYuFieYkRFZbIrJso7dlnvJiJEezpTUrB0Pbt27FmzRpYrVY0Njbi5ptvxtGjR+N+zMDAAG677TacddZZMBgMuOeeezK5XiIiKiLqURx5FwyJMtn8YGHU6YYvIMMgze/HSYfYNTQep0wm+oVaqkuxvKESQPzx+liTZIAytSeCPe4aUqQUDO3evRubN2/G3r17sWvXLvh8PmzYsAFOZ+wHxO12o6GhAZ/73OewatWqjC+YiEgPfP4AHLPFf1xC6MT6fOsZin1Q65A9tL/HZMy8wCK2UE/GOZJDjNW3VZdjWUMFACUYinW+mNo8HSUYkiRJnSgbZKkMAJDSd99zzz0X8fqOHTvQ2NiIAwcO4Oqrr476MUuXLsX//M//AAC+973vpXmZRET68sVn3sZjr3bj+vOa8dmNZ2FFY2WuL0kT+bhwEQiN1s94/ZBlOWLLtGieznTHkJBMz1DfuHKfbTVlWFxbDklSNnePOz2oi5KdEgsX526fFpqrSnFy1Mkm6qCMQlq7XTkorra2NisXI7jdbkxNTUW8EBHphSzLePqNfgDAc0cGseGru/GPPzukPgkXk6k8PIoDACqCZTJZBma9kcdkZHOsHkhumkyUydpqylFqNqI1uDgxVqlMHMURLTMEhAI5NlEr0g6GZFnGvffei7Vr12LlypXZvCZs374dVVVV6kt7e3tWPz8RUT47OerEuNODEpMB157ThIAMPLm/F+v++2X857Nvx220LTT5eEgroOw9EpxzSmXZXLgIADUVSWSGJoNlsholCOqoD5XKool1FIfQzPH6CGkHQ1u2bMGhQ4ewc+fObF4PAGDbtm2w2+3qS29vb9bvg4goX+3vHgcAXNheje/ccQl+/ndX4NKOWnh8AXxrz0lc/V8v4ZGXjic82LMQhHqG8iszZDBIKDUrT5FzJ65E03E2xuoBoDZYJnPM+uCN8ZiKzFB7bTmAJIKhYM9QtAZqIBQMFWO2MR1pBUNbt27F008/jZdeegltbW3ZviZYLBbYbLaIFyIivdjXPQEAWLO0BgCwekkNnvx/l2HHR9fgnBYbHG4f/vv5o3j2zcFcXmZW2PO0ZwgIlcrmTpQNZjkzZCszwxBsSYpWKgvfMSQyQ0uDwVCs8fpYCxeF0Mn17vQvvIikFAzJsowtW7bgqaeewosvvoiOjg6trouISLdEZuiSpaF+TEmScM1ZjXhm61psPK8JANA77or68YVEHa3Ps8wQEGqinlsmG8xyz5DRIKnj9RNRzicL3zEk7rOjXskQxVq8KIKhWKP/IpAbZGYIQIrTZJs3b8YTTzyBX/3qV7BarRgcVP4qqaqqQlmZEq1u27YNp0+fxuOPP65+3MGDBwEA09PTGBkZwcGDB1FSUoJzzz03S18GEVFxGHbMonvMBUkCLl5cM+/9BoMU3DMzpD7hFTK1TJZnPUNAaLw+vEwmy3JYZqgsa/dVU27GuNMTtW9IlMhaq8tgDKaQOuqV6cJTYy4EAjIMhtC0m9cfwFjw88w9pFUQDdQjDjd8/kBWVgQUspS++x599FEAwPr16yPevmPHDtx5550AlCWLPT09Ee+/6KKL1P8/cOAAnnjiCSxZsgTd3d2pXzERURE7ECyRndVkjbmVWTTFiomhQiYaqPNtAzUQffGiw+1TX8/WaD0gxuudmIxSJlN3DNWEgq+2GiUwmvH6MeSYjQjMxqaVz2E0SGo/0lx1lRaYDBJ8ARkj0+6sBnaFKKVgKNZyp3CPPfZYWh9HRETh/UKxV5aowVCBZ4Z8/oB6flY+lsmiLV4cCmaFqsrMahktG9SJsmjB0HhkvxAAmI0GLK4tR9eoE12jzohgRjRP11eWRGSMwhkNEhqtFvTbZzFon9V9MKTvvBgRUZ7Zf0r0C80vkQmiD2S0wIMhx2woyLDm2QZqIPr5ZGKsPptZISA0URZtbUL4jqFwS+uU17tHI3vHQkdxxL/GZrVviOP1DIaIiPKE0+3DkX5lyaweMkOiX6iixJiXPSvRymSiebo5S5NkQmjX0PwG6rk7hgTRN9Q1Oh3x9uE455KFU4Mh7hpiMERElC8O9k7CH5CxqLoMrdWxyxYiGHK4fQV96ni+HtIqqJkhdyiDNahVZih4Plm00fpYmSExUdY1JzOUaOGi0GxTvscYDDEYIiLKG/uCI/Vr4pTIAKDSYlIXAo4WcBN1vi5cFNTMkHcBMkMxzieLtmNIWKouXozMDI1Mx1+4KDRXKe9nmYzBEBFR3tgfbJ6+JE6JDFB2Dom/+ocLuFSWr0dxCNFG69XMUJaDoVjnk0XbMSSILdQ94y74A6FBpaQzQ8GmaQZDDIaIFsSfT47h7h/sx+lJLjij6Hz+AF7vSTxJJogm6kLuG1K3T+dpZkhdurgAZbJY55NF2zEktFaVocRkgNcv4/RE6HeLCJAbEjVQ29gzJDAYIloAj73ajeePDOG3hwdyfSmUp94amILL44et1IQzGisT3r4Ydg2FFi7mZzBUIXqGFrBMNumKbKCOtmNIMBgkdaKsK+xYjtBRHPEzQy1h02R6X4HDYIhoAYhfTqIsQDTXvrASWazdMOGKYaIsdBRHvpbJgj1DwcyQ2+dXMzdajdZPu31w+0LBV7QdQ+GW1gXPKAse2CrLcsKjOAQRLLl9gXlBmN4wGCJaAKLJdWrWl+CWpFeh88jiN08LDZWh4xQKlcgM5eP2aSBUJhOj9aIXx2IyoLo8u9dsLTWpZbDwwCTWJJnQ0RB5er19xguPPwAgcc+QxWRUe5X0XipjMES0AEaD6/EdDIYoClmWk9o8Ha44MkN5XiazBBuog2WygbDmaUlKnL1LhcEgoSYYYIX3DcXaMSR01EUGQ+L7oarMjFJz4g3Zat+QzpuoGQwRaWzW68d0MM3umNV3KpqiOzXmwui0GyVGA85fVJXUxxRHz1D+nlgPAGVmpUwmGqjVfqEsl8iEmihbqBNlhkLj9UowFGqejp8VErh4UcFgiEhj4X+5MzNE0Yj9Qhe0VSX11zwQerIr5CM5Cm20ftCuBCbZbp4W5p5PFm/HkLAsGAz1Tbjg8QXCjuJILRgaYGaIiLQUvhTP4WZmiOZLdr9QuPAyWaFOAuX70kVRJhPTZIN25WdZq2Bo7vlk8XYMCQ1WCypKjAjIyr4hcUhr0sFQ8PMOMRgiIi2JfiGAmSGKbt+p5DZPh6uvVJ44Pf6AOpVVaPL9OI4ydZosGAxNBTNDWpXJ5pxPFm/HkCBJkloq6x51Jr1wUWCZTMFgiEhjEZkhBkM0x9i0GydHlH6P1UuSD4YsJqM6hSWOXyg0+b50sTxYsvT4A/D5A2qTcYtWZbLyyPPJ4u0YCqcGQ2NOtYcs0Yn1AhuoFQyGiDQ2FhEMeQu2pEHa2H9KKZGd2VSJ6mCZJFmFfCSHxxdQp7TytmfIEurfcnn9asAQq2SVqblHciTaMSSIvqGTYZmhRAsXhRZmhgAwGCLSXHiZzOuX4fYFcng1lG9C+4WS7xcSCvlIjvDJSmueZoZKjAa1PDU961ODzpaq+MFJuuYe1ppokkwIX7woeoYSLVwUmoLBkH3GG3EGm94wGCLS2NzR5ymO1xedr/3+GP72sX2Y9ab+ZCL2C12aTjBUwLuGxFi91WKK2Q+Ta5IkqaWy3nEXfAEZBinUr5VtczNDvSmWybpGnUkfxSFYLSb12BE9Z4cYDMXBcgZlw9zRZ/YNFZcJpwcP/f4YXnxnGAeCJa9kzXj8ePO0HUDym6fDFfKuoXxfuCiIUtnJ4B6fBqsFJqM2T52igXpiTgN1osyQKJMN2GfVIDPRIa2CJElqdmjArt+DpBkMxfCLzj5c9O+78KcTY7m+FCpwo9MMhorZrreG4A8ofziJhtdkdfZOwBeQ0VJVikXVqZdeCjszpDzhW/P0XDJBnE92cmQaANCsUYkMCI3Wjzs9Se0YEmoqSiKONCkxGVI6700dr2dmiOb6+YHTmHR58eqJ0VxfChU40TNkNiqlAG6hLi7PHB5Q/793PLW/rMP3C6VzvEOh9gx5/QHs61J6pfI9M1QWLJOJDc/NSZaf0lFTofxbzHj96B5zwheIv2MoXEcwOwQoO4ZS+X667T2Lcf+mc3FBW3XK11ws8jskz5FAQMYbvZMAoPuTfCkzHl9AHR9eXFuOEyNOTDMzVDQmXR788XjoD6beFDNDYpIslf1C4QotM2R3ebFzXw++/2q3uvG4PUEJKNfE4kWx/kCr5mkAqLSYYDZK8PpltXwab8dQuI76ChwMPm8lu2NIeP8FrSlfa7FhMBTFyVEnHMGzaCZnGAxR+sacypOU0SChPRgMsUxWPHa9NQRfINRb2DueWjB0fMgBADiv1ZbW/atHcuR5z1D3qBM7/tiFnx7oU0+Ar68swd9cthR/u3Zpbi8uAbF48VTwsdVqrB5Q+ndqyksw7HDjUJ8SDLXXJhd8zc0MUWoYDEUhskKA8pcfUbrGgiWyuooSdbEcp8mKx7PBEtm15zThhbeH1IbXZLh9fgwEezSW1FUkuHV0Ihgac3rg8wcSNvaeHJnGztd68P+uXp5y9iAd404Ptj11CL97awhiHuXsZivuWtuBD1zYCospuXPYcklMk4m+MK0WLgoiGDoczAy1VSeXOVsaEQxpe43FiMFQFAfDgiE7M0OUATHlU19pURtFmRkqDvYZL/4QLJF9Yt0yvPD2EIYdbsx6/Ukdtto7PgNZBipKjKirSG9Uu6a8BEaDBH9AxpjTkzBr8ejLJ/DTA32osJhwz7VnpnWfqfjVwdN4/sgQAOCasxpw19pluHJFXVr9UbkSvngR0DYzBIT6ho6IYChB87TQERZQL0SgW2zYQB3FG32T6v8zGKJMiLH6eqtFXSzHYKg47HprCF6/jDObKrF6SQ0qLUqwm2x2qGdc6UFZXFeRdnBgNEhqIJVM39C7w8pE1LHgf7V2Ovhv8dErl2LHRy/F2jPqCyoQAkIn1wtaZ4bEriFnsJzYlmSZbGl9KIPEMlnqGAzNMev14+2BKfV1NlBTJsQkWX1lSVhmiN9TxUCUyN53fgskSVL/gk+2ifrUmHK7JbWZNRAn20Qty7I6Hi6agbVWKE3S8VSURBZQtDqxXqiZcyRLoh1DgrXUjPrgdGGyCxcphMHQHG8NTMHrl1FiUv5ppma9aq2YKFWjYWUyG8tkRWNq1otXjo0AAG48vwVA6Ekr2cyQGgzVLUwwNDrtUb/3ukanEViA32tiiZ/W2RQtlYVlhqrLzUmVQDNRWzE3GEp+eu3/vmcxzm62YvXi1LeZ6x2DoTlE8/R7OpRvJlnmX/KUvlAwVBIqk7n5/VToXgiWyM5orMQZTVYAoamfviQnyk6NiTJZhsFQZXJbqEVWCABmvYEFOXpBZIZa0lgomS/Cy2TNGvcLAZGZIbNRSqkZ+u+vOxPP3XM1qsrze3dTPmIwNIdonl6ztFY9r4V9Q5SuUTZQFyVRIrshmBUCQqWgpMtk46JMlt4kmZBsZkgcJ6G+rnGpzOcPqBuNCzkzVB5WJtO6RAZEZoaS3TFEmWMwNIfIDK1qr0Z1MEJn3xCla9QheobYQF0spma92POuMkV2Y1gwpPYMJbGF2h+Q0Re83UKVycIzQwBwclTbJuqRaTcCMmAySGovSyFa8MxQWDCUSomMMsNgKMyE04PuYB1/VVuVetYLFy9SusTSxcjMEL+fCtmLbw/D4w9geUMFzmyqVN/eXit6hhJnhganZuHxB2A2SmjNsISUfDCkZIJE75rWmaH+SSUr1GQrLejsRngwpPVYPQDUhJW4kt0xRJljMBRGjNR31Fegujx08B0XL1I6/AEZ485gZsgamiabYmaooImzyG4MTpEJIhiacHkx7Y7/GIt+obaa8owDhaR7hoJlsvee3QgAODGibWZo0F74JTIgsky2EF9LeM8QM0MLh8FQmDd6lSVXq9qqACiTAwB7hig9404PAjIgScpp1KJM5vEF4Pb5c3x1lA7HrBe731WmyN53QUvE+yotJvWv+kTHcohJssUZjtUDQGMwWxEvM+TxBdATvKZrz20CoH1mSJ0kK+DmaWBOZmiBe4aS3TFEmWMwFOZgr3Jo4oXt1QDCgiH2DFEaRPN0bXkJTEaDupQPYN9QoXrxnWF4fAEsq6/AWcEpsnBivD7ZYCjTfiEgVCabdvvg8kT/vuqdcMEfkFFRYsRly+oAAP32Gcx6tQvKB4okMxQ+Wr8QX0t5iVFd7ZLsjiHKHIOhIFmW8UbwYLxVwWCoqizYQM3MEKUhfJIMULYFiwlFBkOFae6ixbnU8foEu4bU7dNZyAxVlBhRFtx9Ixr25xJZoI6GCtRVKC0Asgx0jWqXHSqGHUNA5NLFhWigliQJV62oR7OtFOe0pHeAL6WOwVBQ7/gMxp0emI0Szg2eIB3qGWIwRKkTwVBdZSjtHZoo4/dUoXG6fXj5aLBEdn5L1NskO14vMkNL0zygNZwkSaEm6unou4PEJFlHfSUkSVJPONeyVFYsmaF6qwWlZgOabaXqc4LWvnPHJfjDP14TkU0mbfFfOuhgsHn63BabepJyqGeIDdSUuvCxesFaasLgFDDNzFDB+f07w3D7Auior8A5LfNLZADQVivKZLEzQ7IsoyeLZTJAKZX1jLti9g2JoGdZMAha1lCBg72T88bts2lgUgRDhd33Umkx4Tdb16LUbFywc9UkSYLJWLgTeIWIwVBQ+H4hoZqZIcrAqDOyTAaAE2UF7NlDokTWHPNJUUz/xBuvn3B54QhOm7VnoUwGhE2UxQqGgjuFljUowdDyhsrg27XJDPn8AQw7iiMzBAArGqMHv1Q8WCYLEpunLwwLhqo4TUYZUDNDVpbJisGrJ5RFixvPa455m/aw88lkOfrZX2KsvtlWmrVzrhLtGhKZIREELVPLZNpkhoYdxbFwkfSDwRAArz+AN09HNk8DQDUbqCkDcxuoAfBIjgI1NetVs3krGitj3k5khqbdvpgZZXWsPkslMiAsGIqya8ju8mIsuO+qQy2TBTNDI86YQVsmRPN0k60UhgJeuEj6wWAIwNFBB9y+AGylJnSENTRWhY3Wa/ELg4qbCIYaIoIhHslRiPonlSf3mnJzxBK+uUrNRjQGA5NYTdTqWH2WSmRA/MzQiWCJrNlWiopgQ+6SunJIEuBw+xIua0yHaJ5urS78EhnpA4MhhEpkq9qrI/6KET1DHn8AMxru46DiFC0zZOORHAVJBEPJHJ2R6IyyU8Gx+mw1TwPxe4bU5umG0B96pWajep1aTJQVS/M06QeDIYQ1T7dVR7y9vMQIc7Cjn03UlIpAQMbYdLSeIZbJCtHpSZHpSPzknuiMstAkWeZj9UK8zJDoCwoPhgBgWX2oVJZtxTJWT/rBYAjRm6cBZbxRLF5kEzWlwj7jhS+glFbrKqKUydz8fgKU3prrvrIbn//l4VxfSlwiM7QomWAowa6hU+PZHasHInuG5pb0Q2P1kb1OIjjSoom6WBYukn7oPhhyzHpxPPjL4IL2qnnvF7uGmBmiVIgSma3UpK7WB5gZmutgzySODU/jqddP53VfXqhMlvjJXWyhjlYmc3l8avZmSW32MkNisafXL8/7w23uWL2wTMPxepEZamaZjAqE7oOhw6ftkGXlL75G6/xfdGLjKBcvUipG1RJZ5FixyAxxz5BCZE9cHr/6b5aPUusZip0ZEoelVpWZ1QGNbLCYjOofbuGlMn9ARnewLCfG6oXlwckyLY7kEJkhNlBTodB9MBSrRCYs1OLF48MO/NW3/oSX3hnW9H5oYURrngbCM0PMNAKRB5r2JDjcNJf6U+kZCgZDp6PsGuoezX6JTIjWRH16YgYeXwAlJsO8axeZoZ5xFzy+QNauw+sPYDh4Dc0sk1GB0H0wFNo8Pb9EBoTG67XcNRQIyPiHnx3C3pPj+MHeU5rdDy2caGP1AMtkc/VEBEPanZOVCZ8/gMEpJRhKpmeopboUBglw+wLzGpqzeUDrXNF2DYmx+o66Chjn7PtpsllQUWKEPyBnNRAddrghy4DZKKG+ggsXqTDoPhgKZYZqor6/egEaqH/+eh86e5TrGLRHP2iRCksoM1QS8XarhRuow/WGne7eMxb/pPdcGXa44Q/IMBulecFtNGajQR0pn1sqy+YBrXNFmyiLNlYvSJKEDg2aqAcmuXCRCo+ug6FB+yyGptwwSMDKRbaot9G6gXpq1osvPfdO6JqmGAwVg2iHtAKhzNCsNwCvP3uliULVVwBlMtEv1FyV/JN7rF1D4mvM5vZpIVqZLNZYvaCO12exb0hduMjmaSogug6GRFbozCZrzK2yWjdQP7TrGEanPWr6fdzpwSwXPBY8NTM0p4G6sjT0fab3UpnT7VOPiQDyt0x2WjRPp/DkHmvXkBbbp4W4maH66EeIaDFeL5qn2S9EhSSlYGj79u1Ys2YNrFYrGhsbcfPNN+Po0aMJP2737t1YvXo1SktLsWzZMnzjG99I+4KzSQRDFy2ujnkbLTNDRwcd+P6fugEA2285H6Vm5eEYnsr+enxaWLEaqM1GA8qCh3PqvVTWNxE9a5JvRPN0Mv1CgrprKCwz5PUH1MAqmwsXhWg9Q7HG6oXwM8qyRV24yEkyKiApBUO7d+/G5s2bsXfvXuzatQs+nw8bNmyA0xn7B6mrqwvve9/7cNVVV6GzsxP//M//jE996lP4+c9/nvHFZyrW5ulwVRpNk8myjPuefhP+gIzrz2vG1Wc2oNmm/PIQf1lR4VJH6+f0DAFsohbEJJkIMoam3HmZFU1lrF5Qy2RhmaHTEzPwB2RYTAb1/LJsmpsZmnb7MBT8w2pZQ4zMkDi9PptlskmWyajwxD5xMIrnnnsu4vUdO3agsbERBw4cwNVXXx31Y77xjW9g8eLFeOihhwAA55xzDvbv348vf/nL+OAHP5jeVWfJN/56Nd7om8TZzdaYt6ku16aB+pnDA9h7chwWkwGfu/EcAEpauXvMxb6hAifLcszMEKAEQ8MON4OhYKBw/qIqTM164Zj1oXfchTOaYv885kI6wZAok4UHQ2Lz9OLack0ai+cGQ13BbE99ZYn6R91c4hT7cacHky6P+vsuEwNTYuEiM0NUODLqGbLb7QCA2tramLf505/+hA0bNkS8bePGjdi/fz+83ugBhtvtxtTUVMSLFqrKzbj6zAY02mL/0FarPUPZC4acbh/+45m3AQCfXL9C/cUpJlAGOFFW0KbdPriDe1uiB0OcKANCZbH22jJ11DwfS2WnU9g+LYgt1AOTs/AFG+V7xrJ/QGs40UA97vLA6w+ESmQx+oUAoMJiUjPSJ7JUKhtIo8eKKNfSDoZkWca9996LtWvXYuXKlTFvNzg4iKampoi3NTU1wefzYXR0NOrHbN++HVVVVepLe3t7upeZMfEX1bTbl7Xpn0deOo4B+yzaa8tw97pl6tvFX1Icry9sokRWUWJEWYlx3vtZJlOIfpr22nI1QMjHYCiVc8mEJmspzEYJvoCsZnpPaXBAa7ia8hIYDRJkWcn0nIgzVh8um03UXn9A7VliZogKSdrB0JYtW3Do0CHs3Lkz4W0lKTIlLLayzn27sG3bNtjtdvWlt7c33cvMmC0svZyN7FDXqBPfeaULAPCv7z8PpebQk6X4C43BUGGLNUkm2JgZAhCatGqvKVezo/kWDDlmverRKS0pBEMGg6QGTyLo0+KA1rn3KXrURhxuNbgRpbBY1GAoC31DQ1OzkGWgxGhAXUXmJTeihZJWMLR161Y8/fTTeOmll9DW1hb3ts3NzRgcHIx42/DwMEwmE+rq6qJ+jMVigc1mi3jJFaNBgi34l3ymTdSyLOPffn0EHn8A685swLXnNEa8X/wlNcCeoYI26ojdLwQwMwQoPwu90cpkY/kVDImSdVWZGZWWlFos5/UNia9Ni+3TQnjfUGjhYuwyGRC2aygLmSHxh1xTlYULF6mgpBQMybKMLVu24KmnnsKLL76Ijo6OhB9z+eWXY9euXRFv+93vfodLLrkEZnP2DirUUraaqH//9jBePjoCs1HCfZvOnZcZawkGQ0PMDBW0WNunBTUYcus3GJpweeH0KJNjbTXledszdDqN5mlBHNjaFzyj7NS46BnSpkwGhPqGhqZm1QNYky+TZZ4Z6hdj9ewXogKTUjC0efNm/PCHP8QTTzwBq9WKwcFBDA4OYmYmNAq+bds23H777errn/jEJ3Dq1Cnce++9ePvtt/G9730P3/3ud/HZz342e1+FxsSuoUwWL856/fjCb94CANy1dlnUv9ZEmWzYEWq6pMIzMh19+7TABurQWH2j1YJSsxFLapUn5J5x17zDTXMp1C+Uev+LaKLuG3dh2OHGrDcAg5Ra71GqRGbo8Gk7Zrx+mAxSwkyUOM3+1JgL/kBm//aiebqF/UJUYFIKhh599FHY7XasX78eLS0t6suTTz6p3mZgYAA9PT3q6x0dHXj22Wfx8ssv48ILL8S///u/43//939zPlafimzsGvr2npPoGXeh2VaKre9dEfU2dZUWmAwSAnLk4jQqLPHG6oFQZmhKx2UyUTpSJymrS2E0SHD7Qiee54N0xuoFkRnqnXCpzdOt1WUoMWm3+F8EQ3/uGgeglOTMxvj3J67J4w/M25idqgFmhqhApVQET+Yvtscee2ze29atW4fXX389lbvKK5kGQ30TLjzy8nEAwD/feA4qYvQeGA0SmmylOD05gwH7LH+hFKixBA3UocyQjoMhMUkWXE5oNhrQWl2K3vEZ9Iy70BRn3cVCEtun0wmG2sPOJzsVHKvX4oDWcKJMdnw4/ubpcEaDhI66ChwdcuDkiDOjMp5YGMvMEBUaXZ9NlqxQmSy9YOg/nnkbs94A3tNRi00XtMS9bTP7hgqeGK1vSNQzpOMyWWjHUKiEo5bK8qiJOpOeIfG1DTlmcTzYnKzFAa3hGqyRQUii5unQ7ZR/+xMZNlEPqpkhBkNUWBgMJaG6LP0G6j8cG8Vv3xyE0SDh3246L+Y6ASF0JAeDoUIlymR1scpkFk6ThY/VCyJ4OJVHTdSZ9AzVVZSgzGyELAN7T4wB0OaA1nANc7KRyxKM1au3CwZDXRmO14sG6nSCR6JcYjCUhNBhrak1UHt8Adz39JsAgL+5bAnObk68IkBdvMjx+oKVeLSeDdSigbqtNvSkKRp9e/MkGPIHZDXTkc6TuyRJ6hllh04r2/q12jEkzAuGks0M1Wd+YKvHF1D/EODCRSo0DIaSIBYvTqaYGfr+q904MeJEXUUJ/v66M5P6GJFeZmaoMM14/OrIeMLRep1mhvwBWS0/hU865dsW6hGHG76ADKNBQqM1vSd3ke0S7ZaLazXuGZoXDKWWGRJHeKSDCxepkDEYSkJ1Gg3Uw1OzeOiFdwEA/3j92TEPSpxLNI6yZ6gwib+MLSZDzCV9Ihhyefy6XKEwNDULr1+GySBFDAmIwOhUnvQMiYCt2aZMuqVDNFELWvcMVZQYURbcam8rNSUdlIgM0tCUG9Np7r8aDDugNVE7AFG+YTCUBLF0cSqFzND2374Dp8ePVe3V+D+r42/pDqdmhqZmEtyS8tFI2Fh9rCcEUSYDAKfbvyDXlU9EGay1uiwiyBBZlNFpN1ye3GfN0jmTbK7wBvH6SkvKW6xTJUmSmh1a1lCZdFBSVWZWM5ldaZbK+rljiAoYg6EkqD1DSQZD+7rH8YvO05Ak4AsfOC+ltfShaTJ3Xi2fo+So/UIxxuoBoMRkgCW4a2ZKh31DvRPigNbIIKOqzKz+rInR+1zqT+O0+rnawjJDWvcLCaFgKLWSnNo3lGapjJNkVMgYDCUhVCbzIJBgQ6s/IONff3UEAHDrJe1Y1V6d0n01WkshSYDHH8C4M/2N15QbY874Y/WCnncNqWP1NfODg1CpLPOjITKVycJFoa0mfHXAwgRDTTYlGFqeZPO0oI7XD6cXDKkLFzlJRgWIwVASRAN1QAamE6Tvf9F5Gm8PTMFWasI/bDwr5fsqMRlQV6H8MmMTdeFJNEkm2HS8a6gvyo4hIZ9Orz+dwcJFIfxr1LpfSLhrbQc2rWrFBy9OvjwPhIKnEyyTkQ5pW8AuEqVmI0rNBsx6A7C7vLCVxm6G3t+trMH/v5ctiblnJpGWqlKMTrsxaJ/FykVVaX0Oyo1ER3EIep4om3sUR7gleTRen42eoaoyM2ylJkzN+hasTLZ6SS1WL6lN+eNWNCrB0PE0M0OigZqb86kQMTOUpGQXL3YH0/tnNqWWog7XrDZRMzNUaMT26bpky2Ru/WWG5h7FES6fTq/vDx4tsSjKdabiwsU1kCRgVVt1Fq5KOyIYOjk6ndaUozi6hJkhKkTMDCWputyMwanZhOP13aPKL/FMzvdp5nh9wRphZigut8+PIYfyfR0tM7Q4T7ZQO90+9Wc90yf3b/71aoxOu6N+vflkUXUZysxGzHj96Bl3Jb2wEYhcuMhgiAoRM0NJUg9rnYnd1Dzj8aup4o5MgiEuXixYLJPFd3piBrIMlJmNUXfgiIChb3wm4bCClsSBo7ZSU8QqhHSUlRjzPhACAINBUpuoUy2VDQV/75WYDKjlwkUqQAyGkpTMyfWnxpUSma3UpI4Ip6NFPZIj9+PFlBrRQN1gTa5MprfR+vCx+mg7cFqry2AySPD4A2oGKRey0TxdiNS+oRQPbB0IG6vnwkUqRAyGkpTMyfWiRNZRX5HRLwT1fDJmhgqK2+fHVDDTw8xQdL1xxuoBwGgIneeVy03U2WieLkQrGtJrohaZNJbIqFAxGEqS2EIdNxgKNk9n0i8ERJ5cz8WLhUPshTIbpYTHr+h1z1BvnLF6IR/G67OxY6gQpTtRFsoM6evfi4oHg6EkVZUlPrleLIpbWp9hMBT868rl8cOR5jlBpI3H/tiFa778Mn5zqH/e+0YdwUmyithHcQhWiz73DMUbqxfECHoux+tP6zQYOiM4BXtieDqlP8QGuGOIChyDoSQl0zPUNRoMhjLcJ1JeYlLvj6Wy/PKLztPoGnViyxOd+MxP3ogIZtTm6QT9QoCey2Sxx+qFfDiwNRtHcRSiJXUVMBkkOD3+lAY4+nkUBxU4BkNJSuZ8MvHLO9MyGRAqlTEYyi9DU271/3/+eh9u/N8/4MCpCQDJj9UD4WUyZobmyoddQ2Jnjt56hsxGg5qZS6VUNsgyGRU4BkNJUpcuxsgMzYT9JdWRYZkMYBN1PvIHZDXgefi2i7Cougw94y58+Jt/wv+8cAzDwfFicZxKPHrMDDlmvWpmNX4wpPz85KpMFgjIakOw3spkQHp9Q+Lfq5mZISpQDIaSlGiaTPwVay01oSaDsXqhhbuG8s7YtBv+gAyDBFx/XjN+e89VuPnCVvgDMr76wrv4398fB8AyWSyiRFZTbkalJfa+V3Ga/ZjTg+kc9MyNTrvh9cswGiQ0WtM7UqeQpTpe7/b51c3regweqTgwGEpSoqWLol8o07F6Qc0M8UiOvCFKZA1WC0xGA2ylZjz0VxfhoVsvhNVigid4hEFDCmWyabcvp8sFF1IyJTJA+bcRi/t6ctA31BfsF2q2lcJk1N+vyDMarQCA40PJBUNDduXnwmIyZOUPQaJc0N9Pepqqgj/ks94AZr3+ee8/laWxeiHUM8TFi/lCbNltskWWAm6+aBGe/fRVWLO0BkaDhIsWVyf8XCIzBADTHn1khxLtGAqXy74hvTZPC6lmhsJ3DHHhIhUqnk2WJKvFBKNBgj8gwz7jRanZGPF+sWOoI0snU/NIjvwjNiI3Wuc/SbbXluMnd1+OabcvqeMbSs1GlBgN8PgDcMz6YMvwyIdCkMyOIWFxbTkO9k6iJ7jVfSHpdceQII7kGHd6MO70JDxegzuGqBgwM5QkSZLijtdn44DWcOIXC8tk+UMcnNtki14GkyQppXOsQn1D+pgoCz+KI5HcZob0eRSHUF5iUqfokmmiFlO0LTrNpFFxYDCUguqy2E3U3VlauCiIMtmkyxu1LEcLT/QMzS2TpUtvTdQplcnqRDC08GVivS5cDJfKRNnek2MAgIvaq7W8JCJNMRhKgegbmruFetYbGqvPdOGiYCszoSxYiuN4fX4QZbLmrAVD+tk1JMsy+tTMUPI9Q7kYrw+dS6bfTIcIho4NO+Lebtbrx4EeZc/WFSvqNb8uIq0wGEpBaKIs8slLpImtpaaE9fVkSZLE8fo8IzJDjTHKZKnSU2ZodNqDGa8fkpRcY7IIhvomXPAv8LSd3nuGAOCMJDND+7sn4PEF0GwrxbIsZcWJcoHBUArUMtmcniG1RFaXnbF6ITRez4myfDAcY5osXSIYmtJBMCTG6pttpbCYjAlurdyuxGiA1x9agLgQXB4fJoI/33oOhkRm6ESCYOiPJ0YBAFesqOMkGRU0BkMpECfXz9011D2a3X4hITRe705wS9Ka2+fHWPBU+uwFQ/opk6XSLwQABoOEtmCj9UI2UYvmaavFpIsJv1hEMNRvn4UzzuLLV48rwdCVy1kio8LGYCgFVTEaqLuDZbJs9QsJoSM5mBnKtRGHEpCWGLO3WE5PZbJUxuoFdaJsARcvskSmqC4vQX2l8sffiRj7huwzXhw+bQcAXMl+ISpwDIZSoB7W6prbMxQqk2UTe4byR3i/ULbKAVaLfkbr1dPqkxirF3IxXq/3hYvhljcEm6hjbKLee3IMAVnZS8QzyajQMRhKQczMkFomy25mSJRjhrhrKOey3S8EhJfJdJAZmkitTAbkOhjSd2YIAM5oir+JmiUyKiYMhlIQLTM06/WjXx2rz3ZmSPmFzMxQ7oWO4sjewZ26KpMleS5ZuFwEQ6d1vnAx3IqG+BNlfzyh7Be6ckXdgl0TkVYYDKWgqkypoYdnhtTT6i3ZG6sXROp5ZNoNb/AQUMqNQVEmi3IUR7r00kDt8wfUxuSUymR1C79rKLRjiMHQiuCBrdEmyoamZnF8eBqSBFy2jMEQFT4GQymojrJ0MXySLNujpXUVJTAbJcgyMOzgRFkuaVMm00dmaMA+C39ARonRgKYUgkkxTTnh8sLtW5gt7P12lskEMVF2atw179//1eBI/crWKnXKlqiQMRhKgdgzNDXrUxfBdaun1We3XwhQxotFJoJbqHNL3T5dxTJZqsTm6UU1ZTAYkv+DoarMjBKT8itqZAH+GPAHZAyoZTI2BDfZLKi0mOAPyOrZi8IfjyslsitYIqMiwWAoBbay0Ej1VLBUJsbqOzTavtpSxWAoH6jnkrFMlrLhNI8xkSQJDZWW4OfQPhganJqFxx+A2SjxBHYo//7Lo2yilmWZzdNUdBgMpcBsNKAyOA4tjuQQZbJsnVY/V7M6Xs9dQ7kkGqgbs1gmswUzQ9NuH2R5YY+cWEiZNJ+Ljxme0j4YOhX8WW6vKYcxhQxWMYt2LEf3mAv99lmUGA1Ys7Q2V5dGlFUMhlI0d7z+lJoZyn6ZDAj9Nc3x+txxeXxqKSu702TK91JABpyehemJyQWxQT2dfitRJhbZJS2JLK8WJe9CpZ5eHzZe/8dgVuiixdUoK0l8tApRIWAwlKLwJmplrF7J2GifGWIwlCuiRFZeYlQzg9lQajbAFMxAFHOpTPRbpZNVa1zIzNC4tlneQhRtvF40T3PrNBUTBkMpEsGQfcaL3nEXZFkZq6/L8li9IHoX2DOUOyIr12wrzerEoCRJumiiHp5Kr2cIABqtomdI++//U6PMDM2lHtg6Mg1/QEYgIONP3C9ERSh7f+bqhCiTTbq86BL9QvXlmp3YzMxQ7oX6hbJXIhOspWZMuLzFnRkSzedp/PuFymQLkRkSZwwyMyS015ajxGSAxxdA34QLjlkfJlxeVFpMuKCtOteXR5Q1DIZSJBYvTrq88PiURYha/vIUwdCwYxaBgJzSaDJlx5AGO4YEkRmaKtLMkCzLGf37NSxQmUyWZfWMwcXMDKmMBgnL6ivwzqADx4en1XLZezpqYTaysEDFg9/NKQovk3VpdEBruEarBZIEeP0yxpyexB9AWRfKbGgXDBVrmcw+44U7+EdDgzWdzNDCjNaPTLvh8vhhkIC2Go7Vh1sRNlEmjuC4gv1CVGQYDKVILF6cnPGETqvXaMcQoIzzi10r7BvKDbVMlsaTeSLFvmtIBJLV5WaUmlOfPBJlsjGnGz4Nj6TpCU6StVSVwWLihFQ4EQy9PTCFfV3jANgvRMWHwVCK1MyQy6tuZV2qcVqdu4ZyS5RoxOOQTVZLcWeGhjJongaUI2mMBuVIGi0zo2KsfqlGKzIKmQiGnj8yhBmvH/WVJTiryZrjqyLKLgZDKRIN1EOOWXWsXsvMEMBdQ7kmRsO1LZMVa2Yos2WVBoOE+kqlT0/LvqGeMY7VxyKCoRmvsgvr8uX1mg2MEOUKg6EUiQbqo4MOyDJQqeFYvdDCibKckWVZLU9m8ygOIVQmK87MkOj1acqgxLgQixfVhYu1zAzN1VFfgfC5jSuXs0RGxYfBUIpEmczrV45PWKrhWL3QzF1DOTM141MbgLUZrS/uMpkaSGaQVVuIJupTzAzFZDEZI/5duGyRihGDoRSJYEhYiF+e4qT0QZbJFpwokVWVpdcAnEjxN1Cnfy6ZsBBbqMWOIS5cjG55cBN1e20Z2pk9oyKUcjC0Z88ebNq0Ca2trZAkCb/85S8TfswjjzyCc845B2VlZTjrrLPw+OOPp3OteaG6LLIk1rEAwVB9cJpsdFr7xXMUKdMG4ESKfc/QkCPztQQNGpfJ7C4vJl1KMMpgKLrzWm0AgKvPaMjxlRBpI+Wli06nE6tWrcJHP/pRfPCDH0x4+0cffRTbtm3Dt7/9baxZswavvfYaPv7xj6OmpgabNm1K66JzqdRsQInRAE9wzHchfnmK/Syj09wztNDEaLgWJTIgFAxNF2kwNJyFhZVal8nEmWSNVgvKS7iHNpqPXdUBa6kJH1rdnutLIdJEyj/5N9xwA2644Yakb/+DH/wAd999N2699VYAwLJly7B371586UtfKshgSJIkVJWbMRL8xdyh8SQZAHXP0LjTA68/wM2vC0jL7dNAWJnMXXxlskBADjVQ53EwxNPqE7OWmvGxq5bl+jKINKP5s6rb7UZpaeQvwrKyMrz22mvweqM/AbjdbkxNTUW85BOxeBFYmJ6hmnJl1wqgBESUOZ8/gBmPP+HtstHzEo+tiBuoR51u+AMyJAnqeHw6RCA1rFHPHMfqiUjzYGjjxo34zne+gwMHDkCWZezfvx/f+9734PV6MTo6GvVjtm/fjqqqKvWlvT2/UrOiibrSYsrol3yyDAZJHd8fWYADK/Xg3p+8gTX/8QJ6g42zsWjfMxQarZdlWZP7yBXR8FxfaYEpg2ymKFGOONwIBLL/b8SxeiLSPBj6l3/5F9xwww247LLLYDabcdNNN+HOO+8EABiN0adztm3bBrvdrr709vZqfZkpEbuGltRpP1YviCZqBkOZ8/oDeP7IIKbdPvz+7aG4tw31DGnbQO0PyOpSu2KRrUCyvlI5n88XkDHhyn5mVB2rX4CSNxHlJ82DobKyMnzve9+Dy+VCd3c3enp6sHTpUlitVtTXR99XYbFYYLPZIl7yidhCreUBrXOJJuoRTpRl7OigQ90dtO/URNzbZqMBOJ7yEqNaAi22UlnogNvMSoxmowG15cEt1Br8MXBqbGGO1SGi/LVgnbhmsxltbW0wGo348Y9/jPe///0wGAqzEVican1u68IFaWowxMxQxt7om1T/f3/3eMzyVGQDsDY9Q5IkodJSnEdyZHoUR7gGjZqoXR6f+jmX1DIzRKRXKU+TTU9P4/jx4+rrXV1dOHjwIGpra7F48WJs27YNp0+fVncJvfvuu3jttdfwnve8BxMTE/jKV76CN998E9///vez91UssI9fvQxnNFXivWc3Lth9skyWPYd67er/D0250TcxE3WR3JjTA5/aAKxNMAQopTL7jLfodg2pzedZOMak0VaKdwYdWW+i7gn2jFWXm1E1Z6EqEelHysHQ/v37cc0116iv33vvvQCAO+64A4899hgGBgbQ09Ojvt/v9+PBBx/E0aNHYTabcc011+DVV1/F0qVLM7/6HKm0mPD+C1oX9D5Du4YYDGVKZIaMBgn+gIx93eNRgyHxZF5fadF0nYHSRD1ThGWy7E3iaTVe3z3K5mkiSiMYWr9+fdypl8ceeyzi9XPOOQednZ0pXxhFElNrzAxlxuXx4d0hBwDgxvNb8PQb/djXPYFbLm6bd9thh7Zj9YK1aMtkwRJjVRYyQxqViXkmGREBPJusYDAzlB1H+qcQkJUA58YLWgAofUPRqE/mGpxWH65YD2tVg8lslMnUzFB2y2TiTDI2TxPpG4OhAqHVX8Z680bvJABgVVs1LllSAwA4NjyNiSjLLMWJ61qN1Qvq+WQzxZMZ8voD6vExWSmTqYsXtckMLWZmiEjXGAwVCNHAOzXrw2yR7aNZSG/0Kc3Tq9qrUVdpwfIG5UnwQJQR+4Uqky0O9qscPm1PcMvCIXp7zEYJNeWZLybVqmeIY/VEBDAYKhhVZWaYjco+mjEeyZE2kRm6oK0KALBmaS0AYN+p+aUyUSbTavu0cPWZykngrxwbhV+DDcu5oI7VW0thMGS+mLQx7OT6bG3qdvv86J+cAQAsZjBEpGsMhgqEJEnqga0slaVnwulRR6kvWFQNALgkGAzt756fGdL6kFbhwvZqdbw+fAdSIRvO8plu4kiOWW8ADnd2eqv6JmYQkJXFlw0ark4govzHYKiA1LNvKCOHgmWojvoKdafMmqVK39Chvsl55cfQ0kBtnyhNRgOuOkPZxr776Iim97VQQtunsxNIlpqNam9VtvqGetTT6isW7FgdIspPDIYKiPjrlRNl6ZlbIgOUfp0GqwVev4xDfaGencgGYG0zQwCwLlgq2/1usQRD2c+qZXuirFuM1XPHEJHuMRgqINxCnZlDwRLUqrZq9W2SJKnZoX1hI/YjYQ3AtVloAE5E9A290TcZdbKt0AxqkFVT+4aylBkSzdNL6hkMEekdg6ECwl1D6ZNlGQd7xSRZVcT7Llki+oZCwVC2G4ATaakqw1lNVsgy8MrxUc3vT2vDGuxoEv1H2coMqQsXeSYZke4xGCogPKw1fQP2WYxOu2E0SDivNTIYEhNl+09NIBCc5hI9L1r3C4Vbd1awVFYEfUMimGzOwvZpIdu7hrhwkYgEBkMFhGWy9IkS2VlNVpSajRHvO6fFivISIxyzPrw7rBzVkc1DRpMV3jcU0HDEftrtwwO/fQfHg1+rFrJ5LpmQzV1D/oCM3nFRJmNmiEjvGAwVEJbJ0herRAYo01wXLxZ9Q8qIvRZP5olcsrQG5SVGjE678fbglGb38/MDffjG7hP4ux++rsleoxmPH1PBo0Wyub27IYsN1P2TM/D6ZZQYDZrvkSKi/MdgqICwTJa+aM3T4S4JNlGLvqFQmWzhnigtJiOuWF4HQNupMtE4fGx4Gr851J/1zy8CyTKzUT2ENhtCixcz//4X+6baa8tgXICeMCLKbwyGCog4ud7p8cPlKa5DPbUUCMg4HBybvyBGMLRmzvJFkX1Y6KyBWirTsG9owD6j/v9DLxyDzx/I6ucPz6plc3+P6N8ayULPUDdPqyeiMAyGCkilxYRSs/KQjToKf/x6oZwcdcLh9qHUbMCZTZVRb3NhezWMBgmnJ2dwenJmwbZPz7XuzEYAyllpjlltDm7tt4fKTF2jTvyi83RWP/+QI7sLFwXRM+Rw+zDjyex8vtDCRTZPExGDoYIiSVKoiXo6O+PFeiBKZCtbq2AyRv+Wr7CYcF6rDYBSKhMn1i9kzxCgnJHVUV8BX0DGqyfGNLkPcR7XplWtAID/ffEYvFnMDg1rFEhWWkwoCza/Z9o3JDJDS5kZIiIwGCo4ob4hZoaSFdo8XR33dqJU9sqxUU0agJOl5TZqt8+v9pz94/Vnob7Sgt7xGfx0f1/W7kOr5nNJktRSWaZ9Q6Jvige0EhHAYKjgqIe1cqIsaW/0xZ4kCyc2Uf/uyCAAoNRsgK00ew3AyQrvG8rWCe3CkF35vrGYDFhUXYZPrl8OAHj4xWNw+zIrPQmDWT6XLJw6Xp9B35Asy2owxMwQEQEMhgoOD2tNjccXwFsDyph6rEkyYXVwE7XICjXbSnNygOd7ltWixGTA6ckZnBiZzurn7g82T7dWl0GSJNz2nsVotpWi3z6LH7/Wm5X7CB1wq0UwJCbK0i+TjTjcmPH6YZCARdVl2bo0IipgDIYKDA9rTc3RQQc8vgCqyswJm2UbrBZ0hC3gy0WJDADKS0x4T4cSmL0cZ6rshbeGcPMjf8Tek8n3FolJspbgZuhSsxGb37sCAPDIS8cx6808OyR6hrSYxGvIwuJFsXl6UU0ZSkz8FUhEDIYKDncNpeaNYPP0BW1VSWV5LllSo/7/Qk+ShUvUN/TMoQF84ocHcLB3MqV+n/5JJVBpDcuI3HpJOxZVl2HY4cYP957K4KqVEtSQWibLfvO52jOUQZmse5TN00QUicFQgeGRHKkRzdOJSmSCaKIGgCbrwk6ShRPB0J+7xueNkf/q4Gls3fk6fMHt0eJYiWSISbLWsDPDSkwGfOovlOzQoy+fgNOd/g4rh9uHmWB2qVGDo0yyUSZTm6dr2TxNRAoGQwWGR3Kk5pDaPF2d1O3FJmogt5mhFY2VaK0qhccXwN6uUBns5wf68PdPHkRADmWxTo07k/68IhhqmdMrc8vFbVhSV44xpwff/1N32tc9FFxJYCs1oazEmODWqctGA3XogFZmhohIwWCowDSEZYayPWlUbJxuH44FDyNd1RZ/kkzoqK9AXYWy6bspiyeup0qSJKw7S1nAKLZRP7mvB5/92RsIyMBHLm3HN/5mNQDl6JBke30G7PPLZABgNhrw6b84AwDwrT0n0174OKThJBkQVibLKDOkBI8cqycigcFQgam3Kk/Ubl8A0xmUM/TgzdN2BGSlkTfZZmhJkvC3azuwrKECl3XUJv4ADYlS2Z53R/DDvafwjz8/DFkGbr98Cf7j5vNRV1Ginv3VN5FcqSxamUy46cJFWN5QgUmXFzv+2J3WNYtJsmaNAsmmYJlswuWFx5feokiO1RPRXAyGCkx5iQmVwSdA9g3FdyjJ/UJzbb5mBV78zPqcTZMJV6yog8kg4eSoE5//5ZsAgLvWduDfPnAeDAYJkiShPdj30pNE39C026euDZhbJgMAo0HC1vcq2aGfHehLK/M4FMzYaNEvBADV5WaUBLeIp7Nra8A+A/uMF0aDxKM4iEjFYKgAiQNbGQzFd1CdJKvO6XWky1ZqxsVh022fWLccn7/xnIipONEELLId8QwEs0K20lBAPdeG85pgMRnQM+7C2wOOlK95WMNJMkDJ3Knj9VOpl8o6eyYBAGc3W1Fqzn5PExEVJgZDBSjURM0jOeJJdZIsH334knaYDBI+/Rdn4B+vP2veegCR3UgmM3R6MrRwMZbyEpNannvuzYGUrzd0ppt2WbVMdg0dDH5PXJhkQz0R6QODoQIU2jXEw1pj6ZtwoW9iBkaDhIsWV+f6ctL2f1a34a0vXI+/v+7MqHuSRJksmfH6WM3Tc91wfjMA4LngsSSpEGUyLQ+4bcwgGOrsmQAAXLS4JsEtiUhPGAwVoHqeT5bQn0+OA1CWLVbEKAkVinhbklMpk6lj9Qmam997dhNMBgnvDk2nfBzIsMbTZEBoomwkxTKZ1x/A4dNKHxkzQ0QUjsFQAVKP5ODJ9TGJIyouW1aX4yvRVniZLFHDc7Tt09FUlZlxxYp6AMBzbyafHQoEZHXkXdNgSF28mNofA0cHHZj1BmArNWFZPSfJiCiEwVABUg9rZWYoJrGosNiDodbqMhgkZdVCoob6AfWQ1sSByg0rlVLZ8ymUyiZcHnj9SkDWoOH27nTLZJ2ih6y9GgbDwh/AS0T5i8FQAeJhrfH1TbjQO670C4WfNVaMzEaDmuk5laBvKLRjKPFJ7ded2wRJUtYTJLvDaDBYtqqvLIHZqN2vlnQXLx4MTpKxX4iI5mIwVIB4WGt8e4uoXygZaqksTt+QLMvoT7KBGlD60sQ5bc8fGUrqOkS/kFY7hgS1TJbikRydvcHmafYLEdEcDIYKUH3Y+WQ8kmM+vfQLCYuTWLw45vTA4wtAkpLv51FLZUn2DWm9fVpoDPv+9weS+/63u7w4OaIcw5HsOXVEpB8MhgqQWLro9cuwz6R3hlQx01swlMwW6oFg83RDpSXudFq4jecpwdC+U+NJlaSGNF64KNRVWmCQgIAMjDmTyw6JBZxL68pRGzx7johIYDBUgCwmI6rKzABYKpurdzy0X6jY+4WEJbXKZFS8YKjfHv20+nhaq8uwqq0KsgzseitxqUzrozgEo0FCXWVqp9eLfiGO1BNRNAyGClSxHcnh9QeSLnnE8+cuffULAcmVyUTz9KIkJsnCXb+yBUByI/ZDC7B9WghNlCXXRK32C7F5moiiYDBUoBqKaLx+0uXB9Q/twV88+DJmvf6MPpfeSmRAKBgacbjh8vii3kZsn25JYpIs3MbzmgAAfzoxBrsrfkl2IbZPC2owlERmSJZlHsNBRHExGCpQ6hbqAs8MybKMz/zkDZwYcaJ7zIU/HBvN6PPpMRiqKjerZdPe8Zmotzmd5PbpuZY1VOKsJit8ARkvvB2/VDa0ANunBXEfyewa6h5zYdLlRYnJgHNabFpfGhEVIAZDBapYDmv9zitd+P07w+rrv01h4/FceuwXEhKVygbUMllqmSEAuH5l4rPKfP6Auvcq38pkB4MlspWttqSbx4lIX/iboUAVw66hA6cm8KXn3gEA3HLRIgDAC28PwesPpPX59NgvJITOKHNGfb84iiOVBmpBBEN73h2B0x29DDc67YEsB5ubF2Baq8GW/K4hLlskokQYDBWoQj+sdcLpwdYnXocvIGPTqlb894dWoa6iBPYZr3rIaqr0WCIT4p1e7/MH1AxKMkdxzHV2sxVL68rh9gXw8tGRqLd5Izi63mi1LMhRF6kcydHJfiEiSoDBUIFSy2QFmBkKBGR85qdvoN8+i476CvznX66E0SBhQ7BZ97dvDqT1efUcDIUf2DrXkMONgAyYjRLqK1JvbpYkCRtjlMpOjEzjkz86gLt/cAAA0LFAB6A2JpkZnfX68Vb/FAAGQ0QUG4OhAtVQwJmh7/zhJF58ZxglJgMevu0iWEuV5l+x5O/5I0Mpj9mLfiGTDvuFgLAyWZRgqF9tni5LO2tzffCxefHtIcx6/Ri0z2LbU4ew4at78OzhQUgScMvFi/CVD1+Y3heQIjEVNzg1qzaHR3Ok3w5fQEZ9pQVtNamXCIlIH/TVWFFERGZoLHgkgTEPTuGWZRk/3teLsWk3rlhRjwsWVcE058DOA6fG8aXnjgIA7tt0Ls5rrVLfd8XyelhLTRiddqOzZwKXBM/GSoae+4WAUDDUNz6DQECOCHr605wkC7eqrRotVaUYsM9iyxOdeOXYCNw+pbfr2nMa8dmNZ+Hs5oWb1GquKsVly2qx9+Q4Htr1Lv77Q6ui3q4zbNmiJOX+Z4SI8pP+njWKRG1FCaTgkQQTLo/aQ5RLe46NYttTh5VXfvcurBYTLlteh7Ur6nHlinrUVZRg6xOd8Af7hG67dHHEx5eYDLj2nCb8ovM0fvvmYErBkJ5LZIAS6JgMEjz+AIYcsxH7hETzdDIHtMZiMEjYeF4zHnu1Wx2xX7O0Bv94/dkpPU7Z9I/Xn42//Pqr+Pnrffj41ctwZpN13m1Ev9BFi6sX9uKIqKCwTFagzEYDasrzZwt1ICDjS79VJsPObKpEVZkZDrcPu94awn1PH8G1X9mNyx/4vdontP2W86P+pa6Ocb85mNIhtHoPhkxGAxYFy0BzT68fCB7FkU7zdLhb17SjzGzE2c1WfO/OS/CTuy/PWSAEKNNh15/XjIAM/Fcw2ziXOknGfiEiioOZoQLWUGnBuNOj7nfJpV8f6sdbA1OwWkx48v9dDluZGUf67Xjl2Cj+eHwU+7snMOsNoMRkwCO3XYzKGKWsq89oQJnZiNOTM3jz9BTOb6uKertw4f1Cq3XYLyQsri3HqTEXTo278J6woFAdq09x+/Rc57TYcPj+DfNKn7n0D9efhV1vD+GFt4ewv3s8Ijgbdij9RJIEXMBgiIjiYDBUwBqsFhwdcuQ8M+TxBfDg794FANy9bhlqgntmLmirxgVt1dh8zQrMePx4vWcCjVYLzohSzhDKSoxYf1YDfvvmIJ47MpBUMKT3fiFhcYzx+v4MFi7OlU+BEAAsb6jEhy9pw87XevGl597BT+6+XM04iqzQmY3WmME3ERHAMllBy5fDWne+1oOecRcarBb87dqOqLcpKzHiyhX1cQMhIbxUlgy9l8iEWFuoB9QT67XfDJ0Ln/6LM2ExGbCvewIvhm0zZ78QESWLwVABCx3JkbtgyOn24WsvHgMAfPovzkB5SeZ/gb/37EaUGA04MeLEsSFHwtszGFKEtlCHgqEZjx8TwQNWMy2T5avmqlJ89EolCP/Sc++oaxkOhk2SERHFk3IwtGfPHmzatAmtra2QJAm//OUvE37Mj370I6xatQrl5eVoaWnBRz/6UYyNjaVzvRQmHw5r/c4rXRid9qCjvgK3rmnPyue0lppx5QolsEmUHWK/UMjiuvllsv5gVqjSYoKttHhLRX+3bjmqysx4d2gav+g8DX9AxqHgVuwLmRkiogRSDoacTidWrVqFhx9+OKnb/+EPf8Dtt9+Ou+66C0eOHMFPf/pT7Nu3Dx/72MdSvliKpJ5PlqPM0Oi0G9/acwIA8JkNZ8KcxX6SG1a2AIh/OCjAfqFw4kiOMacH08EzxAbU5unSot6zU1VuxifXLwcAfHXXuzjSb4fT40dFiRFnNCYuzRKRvqX87HXDDTfgi1/8Im655Zakbr93714sXboUn/rUp9DR0YG1a9fi7rvvxv79+1O+WIoUOpIjNyfXP/zicTg9fpy/qArvCwYv2XLtuU0wSMCR/ql5o+LhWCILsZWaUVOubPMW/2aieTqTHUOF4o4rlqKlqhSnJ2fUfVcXtFXnxUJSIspvmvcMXXHFFejr68Ozzz4LWZYxNDSEn/3sZ7jxxhtjfozb7cbU1FTEC82Xy8Nae8dd+NGfTwFQlt9l+3DO2ooSvKdDCXCej5Ed2vXWEJ4PltEYDCkW1ylng4km6v4s7RgqBKVmI/7+2jMBKEE0wOZpIkrOggRDP/rRj3DrrbeipKQEzc3NqK6uxte+9rWYH7N9+3ZUVVWpL+3t2elFKTYiMzTh8sDrDyzofT/4u6Pw+mVcdUY91p5Rr8l93HB+9MNBZzx+fO4Xh/Hxx/fD4fZhVXs13rMsd8v/8snc8frwc8n04JaLF2FFY6X6OpuniSgZmgdDb731Fj71qU/hX//1X3HgwAE899xz6Orqwic+8YmYH7Nt2zbY7Xb1pbe3V+vLLEg15SUwGiTIMjDuzF6pzOcP4D+eeQv/38/ewPf+0IVXj49GfP63+qfwqzf6AShZIa1sOFcJhg6cmsDQlNL78uZpO2782iv40Z97AAAfv6oDP7n7MlhMRs2uo5AsrlWCnlPjTgDAgD3zozgKiclowD9sPEt9nc3TRJQMzTtOt2/fjiuvvBL/8A//AAC44IILUFFRgauuugpf/OIX0dIyv9fEYrHAYsn9WVv5zmiQUFtRghGHGyMON5ps2SmF/PbNQXz7la55b2+wWnB2sxXDU27IMvD+C1qwclHipYjpaq4qxcWLq/F6zySee3MQs14/vhzMSDXZLHjwQxdqlpUqVKFdQ0pGSJzo3prBIa2FZsO5TfjUX5yB8hIjGq36+bqJKH2aB0MulwsmU+TdGI3KX/GpnD1F0TVUWpRgKErf0N6TY/jWnpO4/rxmfDiFsfef7FcycVeuqENFiQlHhxw4NeZSgy4AMBkkfHbDWfE+TVZcv7IZr/dM4t9/8xZ8wf0xG89rwgO3XKBuuqaQxbVKz1DvuAuyLKvTZHrJDAGAJEm497ozc30ZRFRAUg6Gpqencfz4cfX1rq4uHDx4ELW1tVi8eDG2bduG06dP4/HHHwcAbNq0CR//+Mfx6KOPYuPGjRgYGMA999yDSy+9FK2trdn7SnSq3moBBiJ3DY07Pdj+7Nv46YE+AEpQtHFlM6rKzAk/X++4C68cGwUAPHDLBeq4ttPtw7tDDrwz6MCxoWmsXlKDpfUVGnxFka4/rwX/+ew78AVklJmNuP8D5+LDl7QX9Zh4JsSuob4JF8adHsx4/QCULBsREUWXcjC0f/9+XHPNNerr9957LwDgjjvuwGOPPYaBgQH09PSo77/zzjvhcDjw8MMP4zOf+Qyqq6vx3ve+F1/60peycPnUELZ4UZZl/OxAH/7z2bfVrcO2UhOmZn34yb5efPzqZQk/nwigrlxRpwZCAFBhMeGixTW4aPHCLjZcXFeO2y9fgu4xF+7fdC6WNVQm/iAda7aVwmyU4PXLOHBqAoBybEupmT1VRESxpBwMrV+/Pm5567HHHpv3tq1bt2Lr1q2p3hUlQUyUdfZM4q++tVddQnh2sxX/8Zfn49iQA//01GF8/0/d+Nu1HXF3rvgDMn4WLJHdumax9hefpC/ctDLXl1AwjAYJ7TXlODnqVL8X9DJJRkSULp5NVuDEYa0vvD2EP3eNo8xsxLYbzsavt67F6iU1uPmiRagpN6NvYga73hqK+7leOTaCfvssqsrM2HBu00JcPmlAZPTEQsoWlsiIiOJiMFTgwhtj33t2I37391fj7nXL1aMxSs1GfORSJcuz44/zJ8TCicbpv7xoEcsqBUxMlL01oCwe1FPzNBFROvR9mFMReO/Zjfjk+uVY1V6NDec2RW0s/pvLl+Cbe07iz13jONJvx3mt88fhx6bdauYoWweuUm4sCTZRi2q2HrZPExFlgpmhAldqNuL/u/5sbDyvOeaEVUtVGW5YqSww/P6r3VFv84vO0/D6ZVzQVoVzWmxaXS4tgPDGd4CZISKiRBgM6cRHr+wAAPzyYD/G5uwkkmUZT+5TSmQfvoRZoUK3eE4wxAZqIqL4GAzpxMWLq7GqrQoeXwA7X+uJeN/rPZM4NjyNUrMBH7iQu58K3dxgiGUyIqL4GAzphCRJanboB3tPRRzs+pNgVuh957fAVpp4MSPltwqLSZ0yNBokHklBRJQAgyEded/5LWiwWjA05cazhwcAANNuH359SDl09a/yaLcQZUb0DTXbSuPuliIiIgZDulJiMuCv37MEALDjj90AgGcO9cPl8WNZfQXWLF3Y7dKknSXBYIglMiKixBgM6cxt71mMEqMBB3sn0dkzoTZOf4jnfRWVxXXKuXGcJCMiSozBkM40WC3YtEppkv63X7+F13smYTRI+ODqRTm+Msqm/3NxG25Y2Yw7r1ia60shIsp7DIZ06KNXLgUAHOydBKAsbmSTbXFZXFeOR/969YIfrEtEVIgYDOnQykVVuHRprfr6rdwtREREOsZgSKf+du1SAECTzYL1ZzXk9mKIiIhyiGeT6dTG85rx1VtX4cwmK0xGxsRERKRfDIZ0SpIk/OVFbbm+DCIiopxjSoCIiIh0jcEQERER6RqDISIiItI1BkNERESkawyGiIiISNcYDBEREZGuMRgiIiIiXWMwRERERLrGYIiIiIh0jcEQERER6RqDISIiItI1BkNERESkawyGiIiISNcK4tR6WZYBAFNTUzm+EiIiIkqWeN4Wz+P5qiCCIYfDAQBob2/P8ZUQERFRqhwOB6qqqnJ9GTFJcr6HawACgQD6+/thtVohSVKuLycjU1NTaG9vR29vL2w2W64vR7f4OBQ3Pr75g49F8UrmsZVlGQ6HA62trTAY8rczpyAyQwaDAW1tbbm+jKyy2Wz8xZAH+DgUNz6++YOPRfFK9Njmc0ZIyN8wjYiIiGgBMBgiIiIiXWMwtMAsFgvuu+8+WCyWXF+KrvFxKG58fPMHH4viVUyPbUE0UBMRERFphZkhIiIi0jUGQ0RERKRrDIaIiIhI1xgMERERka4VZTC0fft2rFmzBlarFY2Njbj55ptx9OjRiNvIsoz7778fra2tKCsrw/r163HkyJGI23zrW9/C+vXrYbPZIEkSJicnI97/8ssvQ5KkqC/79u2LeX2zs7O48847cf7558NkMuHmm2+ed5unnnoK1113HRoaGmCz2XD55Zfj+eefT+rr//rXv46Ojg6UlpZi9erVeOWVV+Z97o0bN6K+vh6SJOHgwYNJfd506Pmx2LNnDzZt2oTW1lZIkoRf/vKX825z5513zrveyy67LOHnzhcL9fgCwLvvvoubbroJ9fX1sNlsuPLKK/HSSy8lvMbDhw9j3bp1KCsrw6JFi/CFL3wh4pykgYEB3HbbbTjrrLNgMBhwzz33JP3158vPmp4fB/6cKfL98c3357SiDIZ2796NzZs3Y+/evdi1axd8Ph82bNgAp9Op3ua//uu/8JWvfAUPP/ww9u3bh+bmZlx33XXqOWgA4HK5cP311+Of//mfo97PFVdcgYGBgYiXj33sY1i6dCkuueSSmNfn9/tRVlaGT33qU7j22muj3mbPnj247rrr8Oyzz+LAgQO45pprsGnTJnR2dsb92p988kncc889+NznPofOzk5cddVVuOGGG9DT06Pexul04sorr8QDDzwQ93Nlg54fC6fTiVWrVuHhhx+Oe7vrr78+4rqfffbZuLfPJwv1+ALAjTfeCJ/PhxdffBEHDhzAhRdeiPe///0YHByM+TFTU1O47rrr0Nrain379uFrX/savvzlL+MrX/mKehu3242GhgZ87nOfw6pVq5L+2vPpZ03PjwN/zhT5/vjm/XOarAPDw8MyAHn37t2yLMtyIBCQm5ub5QceeEC9zezsrFxVVSV/4xvfmPfxL730kgxAnpiYiHs/Ho9HbmxslL/whS8kfW133HGHfNNNNyV123PPPVf+t3/7t7i3ufTSS+VPfOITEW87++yz5X/6p3+ad9uuri4ZgNzZ2Zns5WZMT49FOADyL37xi4zusxBo9fiOjIzIAOQ9e/aob5uampIByC+88ELM6/n6178uV1VVybOzs+rbtm/fLre2tsqBQGDe7detWyd/+tOfTuprzeefNT09DuH4c1YYj6+QT89pRZkZmstutwMAamtrAQBdXV0YHBzEhg0b1NtYLBasW7cOr776atr38/TTT2N0dBR33nlnRtcbTSAQgMPhUL+GaDweDw4cOBDxdQHAhg0bMvq6skkvj0UqXn75ZTQ2NuLMM8/Exz/+cQwPD2fl8+aCVo9vXV0dzjnnHDz++ONwOp3w+Xz45je/iaamJqxevTrmx/3pT3/CunXrIpbCbdy4Ef39/eju7k7xqwvJ9581vTwOqeDPWWIL+fjm23NaQRzUmglZlnHvvfdi7dq1WLlyJQCo6b6mpqaI2zY1NeHUqVNp39d3v/tdbNy4Ee3t7elfcAwPPvggnE4nPvzhD8e8zejoKPx+f9SvK16Kc6Ho6bFI1g033IAPfehDWLJkCbq6uvAv//IveO9734sDBw4U3FZXLR9fSZKwa9cu3HTTTbBarTAYDGhqasJzzz2H6urqmB83ODiIpUuXzrtv8b6Ojo6kryFcPv+s6elxSBZ/zpKzkI9vvj2nFX1maMuWLTh06BB27tw5732SJEW8LsvyvLclq6+vD88//zzuuuuuiLefd955qKysRGVlJW644Ya0PvfOnTtx//3348knn0RjYyMA4JVXXlE/b2VlJX70ox+pt8/m15VNenwsErn11ltx4403YuXKldi0aRN++9vf4t1338UzzzyT1vXlkpaPryzL+OQnP4nGxka88soreO2113DTTTfh/e9/PwYGBgDEfnyj3Xe0t8dSaD9renwcEuHPWXIW6vHNx+e0os4Mbd26FU8//TT27NmDtrY29e3Nzc0AlIi1paVFffvw8PC8CDRZO3bsQF1dHT7wgQ9EvP3ZZ5+F1+sFAJSVlaX8eZ988kncdddd+OlPfxrR4HvJJZdEdMw3NTXBYrHAaDTOi5gz+bqyRW+PRbpaWlqwZMkSHDt2LO3PkQtaP74vvvgifvOb32BiYgI2mw2AMmGya9cufP/738c//dM/RX18m5ubo/48AMk/ToX0s6a3xyFd/DmLbiEe33x9TivKzJAsy9iyZQueeuopvPjii/NSdB0dHWhubsauXbvUt3k8HuzevRtXXHFFWve3Y8cO3H777TCbzRHvW7JkCVasWIEVK1Zg0aJFKX3enTt34s4778QTTzyBG2+8MeJ9ZWVl6uddsWIFrFYrSkpKsHr16oivCwB27dqV1teVDXp9LNI1NjaG3t7eiF9o+WyhHl+XywUAMBgif2UZDAYEAgEA0R/fyy+/HHv27IHH41E/5ne/+x1aW1vnpfVjKYSfNb0+Duniz1l0Wj++ef2clnLLdQH4u7/7O7mqqkp++eWX5YGBAfXF5XKpt3nggQfkqqoq+amnnpIPHz4sf+QjH5FbWlrkqakp9TYDAwNyZ2en/O1vf1vtsO/s7JTHxsYi7u+FF16QAchvvfVW0td45MgRubOzU960aZO8fv16ubOzM6ID/oknnpBNJpP8yCOPRHwNk5OTcT/vj3/8Y9lsNsvf/e535bfeeku+55575IqKCrm7u1u9zdjYmNzZ2Sk/88wzMgD5xz/+sdzZ2SkPDAwkff3J0vNj4XA41M8FQP7KV74id3Z2yqdOnVLf/5nPfEZ+9dVX5a6uLvmll16SL7/8cnnRokURX3s+W6jHd2RkRK6rq5NvueUW+eDBg/LRo0flz372s7LZbJYPHjwY8/omJyflpqYm+SMf+Yh8+PBh+amnnpJtNpv85S9/OeJ24nFavXq1fNttt8mdnZ3ykSNH4n7t+fSzpufHgT9ninx/fPP9Oa0ogyEAUV927Nih3iYQCMj33Xef3NzcLFssFvnqq6+WDx8+HPF57rvvvoSfR5Zl+SMf+Yh8xRVXpHSNS5Ysifq5hXXr1kV9/x133JHwcz/yyCPykiVL5JKSEvniiy9Wxy+FHTt2RP3c9913X0pfQzL0/FiI8dVYH+dyueQNGzbIDQ0NstlslhcvXizfcccdck9PT0rXn0sL+fju27dP3rBhg1xbWytbrVb5sssuk5999tmE13jo0CH5qquuki0Wi9zc3Czff//988Z9o933kiVLEn7ufPlZ0/PjwJ8zRb4/vvn+nCbJctiKSCIiIiKdKcqeISIiIqJkMRgiIiIiXWMwRERERLrGYIiIiIh0jcEQERER6RqDISIiItI1BkNERESkawyGiIiISNcYDBEREZGuMRgiIiIiXWMwRERERLrGYIiIiIh07f8HIIsQ7FQFIecAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "data = pd.read_csv('/ocean/cstang/general/7120-01-JAN-2017_slev.csv')\n", "date = pd.to_datetime(data.iloc[:,0])\n", "ssh = data.iloc[:,1]\n", "\n", "fig,ax = plt.subplots()\n", "ax.plot(date.iloc[330:365+30],ssh.iloc[330:365+30])\n", "\n", "locator = mpl_dates.AutoDateLocator(interval_multiples=True, maxticks=10) #auto chooses ticks based on data\n", "formatter = mpl_dates.AutoDateFormatter(locator) #auto formats based on how many ticks\n", "ax.xaxis.set_major_locator(locator)\n", "ax.xaxis.set_major_formatter(formatter)\n", "ax.tick_params(axis='x')" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3.10.6 ('analysis-camryn')", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.6" }, "orig_nbformat": 4, "vscode": { "interpreter": { "hash": "45037fb5df905fe0ec215a1f079d5355b535b26f32fb585b4920c91cfb7d4d84" } } }, "nbformat": 4, "nbformat_minor": 2 }