Skip to content

Instantly share code, notes, and snippets.

@georgebv
Last active June 26, 2022 17:48
Show Gist options
  • Save georgebv/0329912d7a46628226f7558d124ae321 to your computer and use it in GitHub Desktop.
Save georgebv/0329912d7a46628226f7558d124ae321 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,
"id": "68dfa34b-ac80-4882-8d30-ea55e74c79f4",
"metadata": {},
"outputs": [],
"source": [
"import pandas as pd\n",
"import pyextremes"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "c0406a47-1138-418d-b707-d8637c09627c",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<AxesSubplot:>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD4CAYAAAAXUaZHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAABek0lEQVR4nO29d5gkV3nv/z0VOk8OG2Zmo3aV0a4SQglkSWABtgTYWJgfxgZ+2AY/Bl8n8DXYxhbYvra5P4zxvbJJxgZMFmARhAhCAcmrsFppc97ZnRw7d1fX+f1RdapPVVd1V3dX90zPnM/z7LMz1am6puqt7/me97wvoZRCIBAIBGsLaaV3QCAQCATBI4K7QCAQrEFEcBcIBII1iAjuAoFAsAYRwV0gEAjWIMpK7wAADA4O0m3btq30bggEAkFH8fTTT89SSofcHlsVwX3btm3Yt2/fSu+GQCAQdBSEkDNejwlbRiAQCNYgIrgLBALBGkQEd4FAIFiDiOAuEAgEaxAR3AUCgWANIoK7QCAQrEFEcBcIBII1iAjuAkEHMLGUxcOHplZ6NwQdhAjuAkEH8IUnz+K3/v3pld4NQQchgrtA0AHkNB3FEoWui+Y6An+I4C4QdABayQjqJdE5TeATEdwFgg5A03Xj/5II7gJ/iOAuEHQARTOosyAvENRCBHeBoAPQSkZQLwnPXeATEdwFgg6ABXVNBHeBT0RwFwg6gCIL7sJzF/hEBHeBoANgtozw3AV+EcFdIOgAmB0jPHeBX0RwFwg6gLJyF8Fd4I+awZ0QEiGEPEUI2U8IeZEQ8hfm9u2EkCcJIccJIf9JCAmZ28Pm78fNx7e1+DsIBGseTXjugjrxo9zzAH6OUnoVgD0Afp4QcgOAvwHwUUrpRQAWALzdfP7bASyY2z9qPk8gEDSBJvLcBXVSM7hTg5T5q2r+owB+DsBXzO2fBXCP+fPd5u8wH7+dEEKC2mGBYD3Cgrrw3AV+8eW5E0JkQshzAKYBPATgBIBFSqlmPmUcwIj58wiAcwBgPr4EYMDlPd9JCNlHCNk3MzPT1JcQCNY65RWqIrgL/OEruFNKS5TSPQBGAVwP4JJmP5hSej+l9FpK6bVDQ0PNvp1AsKYpCc9dUCd1ZctQShcB/AjAywD0EkIU86FRAOfNn88DGAMA8/EeAHNB7KxAsF4pijx3QZ34yZYZIoT0mj9HAdwJ4BCMIP9L5tPeCuAB8+dvmr/DfPyHlIo6pQJBM4g8d0G9KLWfgk0APksIkWHcDL5EKf02IeQggC8SQv4KwLMAPmk+/5MAPkcIOQ5gHsC9LdhvgWBdIfLcBfVSM7hTSp8HsNdl+0kY/rtzew7ALweydwKBAIDIcxfUj1ihKhB0AFYnJuG5C3wigrtA0AFYnZiELSPwiQjuAkEHICZUBfUigrtA0AEwW6YoPHeBT0RwFwg6gKLVZk947gJ/iOAuEHQAos2eoF5EcBcIVjmUUuG5C+pGBHeBYJXDq3XhuQv8IoK7QLDK4dW68NwFfhHBXSBY5bDJVEB47gL/iOAuEKxy+JIDJWHLCHwigrtAsMopclZMUSh3gU9EcBcIVjnCcxc0ggjuAsEqh7dlhOcu8IsI7oKWMLmUw/X3/QAnZlK1nyyoCj+hKjx3gV9EcBe0hHMLGUwn8zg5k17pXel4eFtGKHeBX0RwF7QEpjZ51SlojKLNlhHHU+APEdwFLaFcxVAEo2bhA7ooPyDwiwjugpbAAlJeE8G9WXgrRrTZE/hFBHdBSyi2QLl/4sfHcWwqGdj7dQoiW0bQCCK4C1qCZcsEpNwLmo6//e4R/NeBiUDer5PQRPkBQQOI4C5oCcyWCaqKIRsBrEdbQhOLmAQNIIK7oCWwoF4IyJaxRgLrMLjxE6rr8eYmaAwR3AUtgVkJhaBsGZZaqa2/4MZulIQIW0bgHxHcBS2BFbgKakKVqdf1mOfN1HpEkUVwF/hGBHdBSygFrNyZYl+PnYjYDS0akoXnLvBNzeBOCBkjhPyIEHKQEPIiIeQ95vY/J4ScJ4Q8Z/57Nfea9xNCjhNCjhBCXtXKLyBYnWgBK3fmtWvrcFFUWblLwnMX+Ebx8RwNwO9TSp8hhHQBeJoQ8pD52EcppX/HP5kQchmAewFcDmAzgB8QQnZTSktB7rhgdVOeUA04W2Yd2hJMuYdVYcsI/FNTuVNKJyilz5g/JwEcAjBS5SV3A/gipTRPKT0F4DiA64PYWUHnoAVcW6Zsy6xD5W4G9LAiieAu8E1dnjshZBuAvQCeNDf9DiHkeULIpwghfea2EQDnuJeNw+VmQAh5JyFkHyFk38zMTP17LljVsAnVwDx3fR3nuZvfWXjugnrwHdwJIQkAXwXwXkrpMoB/BrATwB4AEwD+vp4PppTeTym9llJ67dDQUD0vFXQAwSv39Zstw45hRJHX5c1N0Bi+gjshRIUR2P+DUvo1AKCUTlFKS5RSHcC/oGy9nAcwxr181NwmWEcEPaHK3i8oD7+TYJUgw6qwZQT+8ZMtQwB8EsAhSuk/cNs3cU97HYAXzJ+/CeBeQkiYELIdwC4ATwW3y4JOgAX1oKpCFkrrOFtGL+e5i5K/Ar/4yZa5CcBbABwghDxnbvsTAG8ihOwBQAGcBvCbAEApfZEQ8iUAB2Fk2rxbZMqsP4Ku587ebz3aEpYto0rr0pYSNEbN4E4pfRQAcXnowSqvuQ/AfU3sl6DDaVXhsPVYW6akUxAChBRJ9FAV+EasUBW0hKDrua/nqpDFEoUqSZAlycpCEghqIYK7oCUEXTisFc0/OgWtpEORCRSJCM9d4BsR3AUtwcpzD8xzX88rVClkiUCWyLqcUBY0hgjugpYQeJ77us6W0aHKElRZKHeBf0RwF7QE5o0Hb8usv+CmlSgUiQjPXVAXIrgLWkK5nnvA2TLrULkXSxSqLAnPXVAXIrgLWoJlywSm3Nev517SdctzL+kUlK6/YyCoHxHcBS1BC7iH6nrOlinqFIpMoMrGchOh3gV+EMFd0BLYYqNCSQ9Eaa7nPHetpFt57sD6HL0I6kcEd0FLYEGY0mCUJgto63H5fclMhVQkQ7mL4C7wgwjugpbA2ydBTKqyrJtiaf15zsaEquG5AxAlCAS+EMFd0BJ4dRlEOiSv2Neb56zpOhQzzx1Yn/V1BPUjgrugJfCLjYKYVGVt9oD1l+teLLEVqsblut5uboLGEMFd0BJ45R5EhguvVtebci3phi0jPHdBPYjgLmgJWokiJBunVxC2DK/W13rGTK5ob3+glXQokiQ8d0FdiOAuaAmariMakgEEo9x5m2ct15d54fwSrviz7+HcfMbaxiZUFeG5C+pABHdBSyiWKGJmcA/Ec+ezb9awLXFuPgNNp5hYylnbyqmQwnMX+EcEd0FL0Eq8cg8gFdJmy6xd5cp6zvLWTNHMlmG2zFq3pQTBIIK7oCUUdYp4yOjiGEgqZMB586sVFtT5xuJaiULlFjEJ5S7wgwjugpZgV+4B2zLrTLkbtowkPHdBXYjgLggcXafQKQL23NdHtoybci+WdDMVUnjuAv+I4C4IHKYsg7RliiXd8pzXsnJ1U+6aWRVSeO6CehDBXRA4LPgEmwpJEVNl2/uvRdw9dyPPXRElfwV1IIK7IHBY8I0F7Lmzm8VazpbJFT2UOzehupZHLoLgEMFdEDgs+FjKXWteaRZ1vXyzWMPKNa+5Z8soslT23NfwyEUQHCK4CwKHKXfmuecDKhwWUdePcs878tz5kr+itozADzWDOyFkjBDyI0LIQULIi4SQ95jb+wkhDxFCjpn/95nbCSHkY4SQ44SQ5wkhV7f6SwhWF8yGsZR2QCV/YwEuilqtOJW7rlNQCmOFqvDcBXXgR7lrAH6fUnoZgBsAvJsQchmA9wF4mFK6C8DD5u8AcBeAXea/dwL458D3WrCqYcqSKe0gPPeCFmze/GrF6bkzi0uVJa4q5Nr9/oLgqBncKaUTlNJnzJ+TAA4BGAFwN4DPmk/7LIB7zJ/vBvBv1OBnAHoJIZuC3nHB6oXZJvGwmeceUFXIqGrYPGs5uDmVO7O4FK62zFrOFhIER12eOyFkG4C9AJ4EsIFSOmE+NAlgg/nzCIBz3MvGzW3O93onIWQfIWTfzMxMvfvdUrSSjlv+9of41v4LK70rHQmzTSKKDEICSoXUg61Vs1rJO5Q7GwXJEoEsbBlBHfgO7oSQBICvAngvpXSZf4waTS3rOuMopfdTSq+llF47NDRUz0tbTjKn4dx8Fgcnlms/WVABU9ZGazjJVvSrESilKJYo4qF1kOdeodzdbJnO+f5aSceTJ+dWejfWJb6COyFEhRHY/4NS+jVz8xSzW8z/p83t5wGMcS8fNbd1DKm8BgBYyhZXeE86E6asFZkgJEtNK3cWzKw89w6zZYxJUX8B2Uu5K7ZOTJ3z/X9waBq/cv/PcHo2vdK7su7wky1DAHwSwCFK6T9wD30TwFvNn98K4AFu+6+ZWTM3AFji7JuOIJkTwb0ZLLUpGU2dm/Xc2c0hqnamLfNrn3oKf/ntQ76eW6Hc9c723GdTeQDAQqawwnuy/lB8POcmAG8BcIAQ8py57U8A/DWALxFC3g7gDIA3mo89CODVAI4DyAD4jSB3uB0w5b4sgntD8GozpDSv3IuOFa+dlud+eDIJyVTdtahQ7uZ3VSSpIz13JpSyhVKNZwqCpmZwp5Q+CsDrzLzd5fkUwLub3K8VJS1smaYoWj4xMT33YJS7tYipg4IbpRSLmYJvoeBU7rzF1YmeeypvfO+0CO5tR6xQdSEpgntTlNP3JIRkqWlbxlmILIjUynaRzGvQdIrlnM/gbip29j9T6bYJ1Q4auTDlniloK7wn6w8R3F1ICc+9KYK3ZYzXhxUZEumsCcWFtOE1+1HulFJLsZeVu/G/LHVm+YFycBfKvd2I4O4CG0ouZ4vQO+hCWi1o3KpKVZaangDlbR5FljpqQnHeCu5azYyZQkkHe4ozW0aVCQgxAnwneu4iuLcfEdxdYMpdp0BKDCfrhl9VGUy2TNmaCAVws2gnLEukUNJtlR7dYI+rMqnIc2eZMrJEOky5G0IpkxfXUbsRwd2FVL6sMpYywpqpl2LJrtyDmlBVzOJZnWXLlM+fWtYMU+s9URUFTYeuU1sqJACoEulMz70olHu7EcHdBWbLAMJ3b4RWee6qYtQ070TlDqDmpCpLg+yOqgAMtW+NguTOVO4srVgo9/YjgrsLKe5EFLnu9cNbCUGsULVsGXNRlB/lenw6hX98+JjvlaGtgnnuALCUrR7gWNGwHjO454olqyokK/eryFKHee6mLSM897YjgrsLyZyGiGocGqHc66fskZt57k2nQvITqv6U64MHJvD3Dx3Fcm5lFWM9yp2V+2XBPa/pVtclZst0knKnlIoJ1RVEBHcX0nkNm3ujAERwbwS+cJhhyzQXjAqcLaNK/jx8poJXeuQ1ny4gZFoqtfbFTblbx9KcUO0kzz1X1K0bUafluc+l8tj7oe/jmbMLK70rDSOCuwupvIYREdwbpmjLlgluEZMqSYZy9xHcmAr2u3ioVSxkihjrN86l2hOqxj73csqdHwUBgCx3Tipkkpu76rQVqmfmM1jIFHGogyvDiuDuQiqnYbgrAlkiIrg3gManLioksGwZVTFuFn7y3MvKvX7F+JOjM1gMqNDVQrqArQNxY19qWERuyr3E1XMHDAXfKbZMkvu+nVZbhv39Fzs4W04EdxeSeQ1dEQU9UVUE9wbQdB2EGAFJDWJC1UoHlKDIkvV7NRpV7qm8hl//9FP40r5ztZ/sg4VMARt7Iggrkm/l3m1T7uW0UgAdtYiJBfeoKiPdYbYMC+pB3eRXAhHcHVBKkRbBvSmKJQrV9IhDstR0g2z2+pAs+fac2SKgej33hXQBlNpVZ6PoOsVCpoj+WAjdUdXHhKqb515OKwUMq6tTesiyxYAbusMdqNyNv9WCUO5rh2yxBJ0CibCCbhHcG0Ir6VYwUgOYUNW4dEDDc/ej3I1gUu/fjz0/iOyOZE5DSafojanojig1LSJ2Q7KyZbgJSTahqnSS527ezIa7I1al1U5hMcuUe+de/yK4O2BqIx42lPtKZ1t0IppOLY+YrVBtJt+8wHn4qixZud/VsJR7nQqcBfdsACsqWRpkfzxknEv1KnetxK0ZYKmQnee5b+iOBHI820nZcxe2zJqBlfsVtkzjFEu65RGHFcnc1nhAsuW5S/Up93pvzkyp5QJQ7vNmYOiLm7ZMzVRIU7nHOOVeqrRlOka5m9fSxu4wiiXaUaWaLc+9g69/EdwdMOWeCCvoiSoiuDeAVqLlWihmUGrGJ+YnFRWfE7R5FtzrnFANVLmbq1P7YyF0R9Saowg35V7U7ROqneS5W7ZMVwRAZ2XMlG0ZodzXDMwbTDBbJle7VKvATlEvK3f2fzOqrZGG241OqC5mjYu5WnAv6RSnfDR8ZqUH+uMhdPsQCnlNR0iWrF6x+WJ5haqVCtlRnruGWEhGImI0fMsUO8d3X+JSITv1+hfB3QEbSiZMW6akU1utGUFttBK1bISQZcsEoNzZIiZfqZCN5bn7mVD93ouTuOMffmI1f/aCDe2NCVXDlqkWKHLFEsKKZLUTNJS7s/xA53juqZyRdcZ636bznaPcWZaM1sHXvwjuDuy2jDE8FtZMfWi6ztkypnJvMrjLEoEkEWMRj69FTI3lubMSz7kqyn1yKYeSTmsO2eczBagysTKvNJ1WHRHkijrCqmyVKzA8d+NYEtKJnnsRXREVsZCh3DvKlskUEDdvSp2aMSOCu4NUXgT3ZimWqBXUWaBqbkKVWt69KvvznJudUK0WiFidFLboyIuFdAF9sRAIIeiOqOb+eKvAvGYod0kySiXnNGOFKrNkgE7z3DUkwpxy75CFTCWdYjmnWSuLRXBfI6Q4W6ZbBPeGsOW5B+C5F7jsG7+2TCtTIVkzl1qTrvPpAvrjIQBAd1Qx98f7XMoXdasaaViRkC/qthsl0HmeO2/LdIpyZ4Jg+5AR3Bc6dFJVBHcHqbyGkCwhrMiWche57vWh6dRadBOE565xAc5o1lH9vSilyBVLIMT4e9ZTRZFlSVQLRGzSvVawWsgYyh2Ar3Mpr5Usvz2iyshrRlVIdqMEDM+9c4J7Ed2cLdMpyp2dA9uZcu/Q618EdwepnGbN7gtbpjGMPHd7KmSznjv/frU8d02n0KmRggjUV0pg2YdyZ8G9mi8PGJNylnJntkwV5Z4r6ta6AKbc+RslYNgynTKh6rRlOqWmO1PqWwdiADo3HVIEdwepvIZ42DgZRXBvjBKv3ANKhSwvv6+t3FnQHeoKA6hvUpVdyNliyTOzhVl3tWyZhXQBveaCJD8Wn1O5sxWqisQr986p5+60ZTql1R6bVN82aNoy6c68/kVwd2CoDeNCTIQVyBLp2AmVlaLYglRI9j6qbKQCVkspZH77cLexeMZvOmSxpCNdKCGqyqC0/D5O0oXayt0oGsZ57uZosNq+uCp37lgC5silA5S7VtKRLZbQFVERD7M8985Q7mytw2AijK6wYv3eadQM7oSQTxFCpgkhL3Db/pwQcp4Q8pz579XcY+8nhBwnhBwhhLyqVTveKtJ5DV3myWhkOYhVqvWiuSxiaspz51Mrzf+rBTgWdIdN5e7378eet6mn+opKlq9dLVtmOVeETmF57l2R2p57ruii3HX7hGqnlPzlExPCigRCgEyH5Llb6xOiKnrjaseKOz/K/TMAft5l+0cppXvMfw8CACHkMgD3ArjcfM0nCCFyUDvbDlL5sucOQNSXaQB7+YEAsmU0bkLV/L+a786Cbr22DLuIN7Lg7qE00z5sGX51KmCMYKKqXD1bRtOt4F723HVHKmR9i5h0neLDDx7CmbnaK2qDhM1zdEUUEEIQDykd5Lkbf6PuqIreaGjteu6U0kcAzPt8v7sBfJFSmqeUngJwHMD1Texf20nljUkghgju9cMXDgspbEK1cbXpNkFbrTIk62jElLvfbCf2d/Yd3KsEKxYg+szgDhjpkNVtmZJlyzDlXuRulED9nvvkcg73P3ISPzg07fs1QcCCO7OjoiG57j6q5+YzK9ImcSlTQHfEsGR7Y2rH1nRvxnP/HULI86Zt02duGwHAt7AZN7dVQAh5JyFkHyFk38zMTMXjBy8s4/h0sonda4xkTrM8QgC+qvkJ7Gg657nLhhJtpmEHb/OwQOdHubOCVX4DxJLprdayZZjlkNOqBHdTufeZE6oAapb9dVXu3I0S8J/nz2BBtt311FnRMDZ/FQ/JdSv3N//rk/j77x0JfN9qsZgtWjfl3tgaVu4e/DOAnQD2AJgA8Pf1vgGl9H5K6bWU0muHhoYqHv/AAy/gQ98+1ODuNQ7rwsQQyr1+NC67RVUCqAqplW8WZVumtnLvi6uQJeJ7QrVsyxgNrd2UO6XUavZcrSywVe43xin3SPXg7qbcNZcVqvV47izItjvHnLdlACDagC0znczh0GT7Bd5ipmg1Ke+Lqesrz51SOkUpLVFKdQD/grL1ch7AGPfUUXNb3Sxni22/Y7IZfmHLNIfdRgmgtoxtgpbZMlWyZUzlHlHluibELVum21u55zXdCq7VPPcFh+cOsFGge5CllCKvGbVlAHu2jOpYxFQrW4hnpZR7Km8P7vE6bRmtpCNX1Ns+VwAY6bA9sbJyX8oWO2IS20lDwZ0Qson79XUAWCbNNwHcSwgJE0K2A9gF4KlGPiNTKFlFvNoFy4JwBndR9rc+bLaMEkSeu0v2TZX3Y8o9osi+epcymHLf0G149W7Bmw+S1bJlFjJFhBTJyvEGUPVGw9IuWfmBiCqbPVT1ikVMAHwHG1bltN2ZKpYtY/Pc/e8DGx1NLefr9uqbZTFbtOy03qhq9tTtPIHnJxXyCwCeAHAxIWScEPJ2AH9LCDlACHkewG0Afg8AKKUvAvgSgIMAvgvg3ZTShs6qbLFUd12QZknm7SckAFH214M/+PJ+fO5nZ1wfK5b0ikVMQRUOs2yZKhOquWI5ULJSu35YyhbRFVGsm7tbHjtftraWcu83i4Yxqt1o2GgjrHDKXdNtN0qg3JHJr+/OglK7z99la0KVee5KXUGa39/Ts5lgd64GNlsmbvzfiZOqSq0nUErf5LL5k1Wefx+A+5rZKQBtv1sD9oqQDH6VKstVFgA/ODSFvKbjLTdsrXjMXsUx2MJhLM+92s2CKfewKlsjLz8sZYvojamIVlkun7Ipd+/gPsetTmXwNd35oM/vs1U4TJWR13TzRmn33IE6lLv53dudhpjKa1BlYs0hxOpV7nxwn0vjss3dge+jG0ZFyGLZloka/y9kCtiOeFv2IShW5QpVXafIFQ3PrZ3lTdM1grvAQNcplrJFZD1uwEaxK+PUkiUCiQRYOKyOPPeIIvnqgMRYzBTQGw1ZnZDcPHc2MUlI9eB+bDpp1SZhdEcV6LRsObjtM6/cAcNOUWyLmNjIZXUr92SuiERYsW5isXB9wZ3fXz9dr4LCuPnCUu7sBr3Ugcp9VQZ3PsWsnb47UzlOWwYQwZ1nOWdcAF7ZJM7c7JDirzWeF/wEreIjz50F3bAq123L9ERVKx3R7fuxoNMXC3naMnOpPM7MZXD1lj7b9u4qq1Sdyp3tA1PAjHIqqL/juXKpkJptpBur05axKfc2BneWGcPsmN5YWbl3GqsyuPN3+Hoq+jWLNcPvyHMHRNlfHjbx6KbEmF3ATwKqstRkVUhqqVdVqq3crclJRapvQjVbRE9MNRpmeChzFnQG4iHPPPjnzi0CAPaM9dq2W2V/XfbHS7mn85ql1oHyzc2vLZNaKVsmZ08pjoVk5Iq67/1Oc9k2p9uYMcMy9JgdwyZW6ylB8K8/PYl/+tHx4HeuTlZlcOcvGjbJ2Q7YhcAvYmJL2GeS1ftlrieq1TxndgE/CRiSpaazZUJcsw6gunLNFUuQJQJFltATVZEr6pYyrsaSOZFGCEFUdbcRWNbJQCLkmS3z7NlFyBLBlaM9tu1loVApWHIeyj1dKFnzDACn3H0GSTbf0H5bxr7S22rY4bN4GGuIcvnmbpyea9+EKju3e8yg3hVRQUh9ZX+/8vQ4Hjww0ZL9q4fVGdyLK6vceVtmKBFGSJYwvpht236sdviyuE6sZtZysLZMuVaNjzx3TUfEVL5s+Xut84hSatkygJG6V82WGUiEPT33Z88t4JKNXVaTCkZVW4bLzQfKyh1ARbMOoJ4JVTbKane2TLHClgH8l/1lyv2KzT2YSebbdnNadCw+kyWCnqj/hUyUUpyZy7TdBnNjVQZ3XjG103NnJ1CcuygliWBzbwTjCyK4M6rZMswucdoyTadCciV/gep57nx1Rb+tEtMFYzUom0CLqLLrClR20Q55BPeSTrH/3BL2bumteIy12nPbF2ueQLErdwB2W8bKFqrPcy+WqK/RS1Akc5p1YwVQd8MOdi1eMWKMftrlu/MVIRl9sZDvVMiZZB7ZYskaeawkqzS4lwN6O4eTqZyGeEi2LfcGgNG+mAjuHEzduAU/NtGpOmqQN2rLUEqNVEjzb6JY2SLVFjGV66JXU8s8LOAy5R7zUu4FDSFFQiKsuDb0ODGTQiqvYe9YX8Vrq3VjKi9iqlTuqkueu/9FTOXPSrcx4Dirq9bbao9NJO/e0AUAbfPdF7mKkIyeqOrblmEWklDuJs7FnzbPvY0rw5wnJGO0L4rzC+1dSLGaYUPUjEtws5S7HMyEKgtiFeUHqqZCVir3Wrnu7OLtMSfSvDz3tFk1NBqSodPKsgrPnl0AAFfl3lWlYYdTudtsGRfl7tdzT+U0K7WzXQGHUmPRn3NCFfDfJDudNwr4bRs00knbpdyXskWrIiSjL+a/pjsrl5Atlla8ZMGqCO4HJ5Ztv9uyZdp4B0zm7RUhGaN9UcymCjV7Zq4X2Ile0mlFcCvbMuWLI9yE51503Cz8rFDNa+XOTT1RFlBrKHc2HOdsGffyAyXEw7J183BOqj57dhE9URXbBysXvCiyhHhIdrVlnMqdt2WctWUAf8qdUopkTrNKGLereFimYAQ23nNnrSvdcvzdSOU1xEMKYiEFG7rDONWmVaoLmYKtTDNgpEP6TYU8w03+rnRD8FUR3HVqL4S0UhOqfBcmntE+Qz0Ia8aAH6LmCvbgVrZlnJ57g8HdYfMoPlao2pR7xJ/n7rRloiHZMxUyHlKsrBbnc549u4i9W3orVqAyuqOq62iUvU+EpUKq9u5LjHo895zZYJvVymmXLeO20juqGj+7LXz78r5z+PgPj9m2sRESAGwbiLfVluH9dsC44ftdxHRmngvuK2zNrIrgDtgDOhu6SaTNE6o5d1tmpM8oATsurBkAsGUOZIr2v0/ZluE998ZTIdnEqbNwWK08d8tzr5JbzsO+E1PusZDsuUI1HlZcV7Emc0UcnU5W5LfzeJX9zVklE8wJVaWs3HmLqx7Pnd1ENpkljNsVbNjn8rYMU+5uVtc391/Al/aN27axERJgBPd2VYc01jo4lHs0hGRe83VD5fdTBHcTfuKUnQADiXD7PXcPWwYQyp3B+4/OAMguAN4nDilSw52YNIfnXi6cVcWWcfQiDSlSzZruTLmzxSsRD889lS/Zgju/mvr58SVQCuzdUjmZyvDqxlQuHMZqy3AjH0cnJsCf587mGTaYJYzblQ657KjlDsCq1+Nmy8ynCxUTlinOIt02GMdsqtCWWLCUKdgarADl1ap+fPfTs2mM9BrxIqiMmeVcEQcvLNd+ooPVE9w5hc6GbgPxUFuzZYyFF5XFwYa7IlBlIoK7yWKmYNkSzgBYngC1K/dGOzEVNLsto/qoMpkr6rYJyVpNMgDjwg3JkvW9oqq3LZPgPHf+5sYmU/eM9np+TjXlHlIky87hlbtsKz9Qe+TCYNcO6yzVrvS8pKMiJFDOlnGzZRbSBSznNNtohLdltluTqq0fOS+42jLGDZ916vJiMWN8j0s3GUXOglLu//jwMfzS/3m87rLjqya4O0upRlVWi7s9wV3XKeZdKvkBhlra3BvF+SYWMp1fzGIutTZWuS5mi9ZQ3xkAmaLmrYSQQhrOlnEq93K2TPV67vyEpKGWa3nuBfTEVCu4Rj0nVJnnXrni8uRsGpt7ItbqRje8yiHki+WFV4BTuVf679VGLgymdJlyb78tUz4O5YwdF+Vuqnb+b8SPoreZk9OnWmzNOCtCMliwn09XP4dYGuTlZgXLoITpf59eMPpb1Pl+qya48/m4mUIJsZDRRaddnvvZ+QyyxRIuNvNqnYz2RZvy3H/rc0/jz791sOHXrxZYRUimBp3KnSlq3koINTOhWrJ77mXlWr2ee0S1K3c/E6q8YmN57k61xOwCZjPkuWyZ5axWERicdEc8bBmtZHVhAso1ZgDn/EU9nrvxOe3OlmHfjy3aAoybUkSVKm6Y2ULJyjhadAR3Zsts7TeC+9kWB/eFTAGUAv2OmzPrplUrY4b57Sy4B3EzzWsly5Kpt3jhqgnuNuVeKCEakpEIK22rLXPY7NV48UaP4N7b3EKm8YUMTkynGn59NT7y4CH8xbdebMl7O0nmNFBanqRzBnevPPdGbRnLw7dsGX/13Png6GcEuJgplx4AgEhIBqXlFEXASC3MFFgqpPH9+GBlLLmv3iKBZcvojuCcd9yQZIlUZAix7YA/z50p6MFECLJEVlS5A8bKb+c+zHMBk/nulFKbLcNiwWyqucqMJZ1WFRlTyzkA5ZshgwX3+XSt4G6IvyBtmUMTSWvU28HB3T6hGgvJ6IqobVPuhyeXQQisFXFORvqimEnmG8p1L+kUi9liU7ZONZ44OYdHj8225L2dMPWyude4AJzHo6jbgzEAqErji5hYEGeFwwghkCVSsxMTHyh7oiqSfpR7rNJG4D111j/VK1tmOVu0+cxudEdUs6a7/bzOaSWbzw6UfXdbtkwdnnu5SbWKWEhuWypkMqdBIkbfVJ6oSwbSAhcwmXLPFXXo1F7AbyARqhlca/Gn3ziA3/zc056Ps+DObCyG3+B+ei6Njd0RDCSM5wcxx/GcOY8D1F9TftUEd36xUsb03BMRxVSK1U/k5VwRf/yV520nSr0cnkhi+0DcGm47YRkztQL0Y8dn8X9/csK2jQ33lrLFlkwQL2QKmGmTn88uwFrKnfeJm6kK6VTugKFkvYKb0Wjartz7Yioml3MVx/7CYhYPPHce//bEaZxfzNqWnEddPHU+f9stWyaZ02xWhBs9Hitmc0Xd5rMDZd/dXbn78dzL+5sIV6rmVuFs1MGIh5SKmxofMJlyLx/n8t+wP958cD85k8ahCe+sk6ll4xpyBveIKiMWkmvGl7NzGWwdiCGqypBIMMr92XOLYIdxTSj3nGnLdEUUaDq1DY3deOzYLP5z3zl8/+Bkw59/ZCrpackA5YVM52tYM1/edw4fe9i+IIM/KWu9vhEWM0UsZoptKQzFLkDmuTs9VM0lGBtVIRtLhXR67uxnr/crlih0Cptyf/3Vo8gUSvjkT0/ZvscvfvxRvOeLz+GDD7yIxUwRV3FZLlGXErXsHI2HFMsftyn3nA/l7rFiNu+i3NkNyv7d7Z77Bx94AX/6jQOun8XK7soSqbvNXTMs5zTbjZLh1iR7wWbLGMfEOs68co+HMNdkcE/lNcwk8xWWGGNyKQdCymW+efpitW8up83gTghBPKQEIuSeO7eIl5jF0zo2uPP2S6aoIRZSrNWitdLYjk4ZXjZrklAvmYKG03NpXLLRu0+j31z3+UwR6ULJllM8x3mF5xeDTefSSrql0Oaa9CT9wC7ATaYt40xtK7qmQpKmyw84389LuVr9U7lAuWesF3ddsRH3P3LCylj68IOHsJAp4t/f/lLs+9M7cOy+u/DWG7dZr3GzXVJc0LGUuxn8dbOJultQ4/EqZFZNuctVPPenTs3jmTOLrp/FFDRgqPd2pRUnc+79huMurfbsyt3eEpAP7v3xUNPZZum8Bk2nnhOjU8s5DMTDtpspYyARss0PuL33bCqPrQNxa9+bVe7z6QLOzGXw8t1DAOC77DBjVQR3iZCKRUxR03MHaq9SPTZtTIY+d26poc8/OpUCpcAlm7yV+4buCBSJ1MyYYcp2Nlk+EZpR7tPJHD7z2ClPa4r/g8+2wZph328oEYYsEW/l7ij5q+nUUzE5+evvHLY62Wguyl2potyt/qmOQPn7r7wY2WIJ//SjE3jixBy+tG8c77hlO27eNYjBROUF7a7cjZ/jYRmqbHj/7POSeWOiudvHhCpQacu4KfeIpdy989ynk3lP8ZPkuiEFEWz8suzowsSIqkqlck8XQIhx82HK1K2X8UAibNqbjRfjYjFm2qPxztRyDht7KlU7UFu5s8lU1jc3Hpabzk7ab4rVGy8ahCqTzlTuzpn8bMH03MP+Gi0cM5X7kcnlhlbhHZk0fLhLqtgyLNe9pnI3TwDeA59Pl3+ut+nHRx48jD//1kHPz+VXzbWjW9QiV4Ml5rKK0638ACvi5XdS9Vv7L+D7LxoWm9uKV1XyHgm4KXcAuGg4gTdeO4Z//9kZ/NFX92OsP4r33r7bcx/cqhiyizVu+skRpZzax5S4nwlV/vkMYxLYYctYnntlnntJN5rHz6cLnhd9Ml/O3omFFN9Fu5rFa2LZUO6V2TJ9sRD64yEXz91uyxRLtKl1L7WC++RyHhsdfjujlufP0iC3mcrdGCk1d7yfPbsAiQBXjvSYZYc7MLhLpLL8QMz03IHqiwG0ko6Tsyns3pCAToED4/Wr90MTScRCMsb6YlWfN9JbO9edTbrwgZalcBmlg/0H91OzaTzw3HkA8PQb+WXbbQnuGSNgKLLkmv3gVjgsZK0qrR3ctZKOyeUcJs3MBStbRuFrmkueee5WL1K18tR+zx27AAKcm8/ivnuu9Jw8B+C6SMmpKPniYtaqzBoTqpbnnnMG95JtVS1QVu6yS567plNrpJZ0rO5kpLgm1Ymw3MYJVc11BOPm+y+ki+iLqeiNlbsdedkyQO2MFS+0km6dG9PmueVkejmH4SrBvdqEKisYtsVS7s2PlJ49t4iLN3YjHlbQE/Xf6J2xSoK73ZbJFs08dx8t0k7PZVAsUfzyNWMAGvPdD08u4+KNXZAcTTqcGAuZvINzXitZ6mjWptwL6Imq2DpQX678x394HOya9Tqx2q7cM+VVvG6t6Jhy531iFuj9ZMxMLOVQ0ilmknloJd1VuSsy8WyzZ1VXVCsD96aeKP7q7ivwh6+6GLeaPqYXTk8dqJzoCyvl78+CdS3lzm4MzoVMeU23LWICyjcotxWq7Bgx3Oqu8LZMLKy0rbZM0iPfPxZSKtrszacL6I+HbMqU2V8J1+De2DnOp4G6Kfe8VsJculBVuacLJc9U6LPzGfTFVOvv32xw13WK/ecWrSJ0PdHaC/GcrIrgLkvl4F7SKQqajpiqWAeqWsGg46bf/tId/djSH6s7uFNKcWQyWXUylTHaF8N0Mu+ZleIVaOfTBQzEQxipo4TBmbk0vvHcedxx6QYA3sqdnxxqi+eeLVr9Jd0aWljZLQ7P3Xistl96zhwZ6dT4zlZqJd+ZSPJW7iyzyqmCGW+8bgzvvu2imvvBVL2t5SMLOqFK5c5UldtEIo8iG12c6lHuiovnXizZg7vbhc973+2aUNV1iqTHxHIsJFc0eFkwbZneWKjCc49zqZADccMLb3QhU4q7sbmJoGkrDdLdc6+1SnUxU7CeAzR/vE/OprGc07CXC+6LNWrbOFkVwV0iZc+dqYtYqOy5VztILFNm51ACe8Z66w7u08k8FjLFqn47g2XMXFh0H9bxQ0Y+0M6l8+iPhzDSG/O9EOqffnQcskTw/ldfAqC2ch/uCrcl151fyelW85xlcbh57n5sGX5kM7mUs3x6Z6s5rzz3fBXlXg8xq/54pXKPmUHHKC5m7N+yT1sGYCUInKmQ1Tz3ymyZkq7bFKhbSQM+ayUWMva1WtmGIEgXjIllL+VOqd3qYsq9l2tl59bLmC0MatSW4ZMyppOV1y/b5sxxZzBB45WRluQsMMA9M6genh9fBADsMTt68Tc/v9QM7oSQTxFCpgkhL3Db+gkhDxFCjpn/95nbCSHkY4SQ44SQ5wkhV/vZCWNC1TgQ7GKK+LRljk2nMNoXRTysYM9YLyaWctZKMz+wRQ31BPdz8+6+Ox+AK5R7ImTVhZ9Yqr5/5xez+Noz5/Gr12/BjsE4VJl4e+7ZAmSJYPtgvI22jHGiu3mobnnuLDD7mVAd547t5HKunC1js2UkT1umlnL3SyRUWV4gnTf6p7KRSESVrPM16dOWASqLh+nmaNWZ4eO+QrXsuVdT7gVNR17TrXRiJpQyLe4mxq+KdcIyUZg4opRanY96Y4btoOtG6YFYSLbZpM167uyGQUhZpfNMLrkvYHJ+vpdyd2YIxZtU7mfnMyCkPEHbE/XfMITh5wr4DICfd2x7H4CHKaW7ADxs/g4AdwHYZf57J4B/9rUTpHxxsIsppspQzRKs1Q7Ssakkdg0nAJTvcs+eXfTzsQCAI2ZNGV+2TH/1jkysQ/pgIlzhuffHw1ad51qTqk+enIOmU7zp+i0ghFSdzGElSoe7I23LlmEFtqKqy4SqywpV9r2fPrOAWowvZC2/e2o5V34/hZ+gJVUmVINR7iFZgkQcnnvBXu+fb8XHlHOt2jKAWfaXU9rspufM8HFT7pJEIBHDvuQVqDO4s2uGz5YBWl8ZstrcwxZWAGw+be1jsUTRb9oyOjVuDqwhCk9ElREPyQ2v5WDHY6Q36uq5W3VlagR3r5tL0rGALRFSUND0htd3TCzmMJgIW6NeVh+pnr6sNYM7pfQRAPOOzXcD+Kz582cB3MNt/zdq8DMAvYSQTTV3QiJIF0pWYSagnIrWFXFvSwaYmTIzaasezGWbuqHKpC5r5vBkEptqlGllbOw26rqf88iYYYscdm9IWBYJKyU8EA9xJQyqZ9ycmk1DIrD6cPbFvFfnLWWMmiiDiVDLgzurCNlnTagqlROqug6JwKa6rtnah90bEvjMY6dr5imPL2RxxUg3FInYbBk+wCmS5G3LBKTcCSEVcwp8dyDAXvN9OVdELCTbVLYX3VG7516+ITmUu1q5QhUwvj/z3FlWijO4s2smYU3weZfcDZKycq+8ybEccJYTvmCW0O0zbRnAGImylbVO+hMhzDU8oWrs146hBKaTuYrzcGo5h5AiuZb8BvwE90rlzn9uvVxYymIzV8CMHR8+FrJMOi8avQI2UEonzJ8nAWwwfx4BcI573ri5rSoyISjpFLmibl1MbEKrK6x42jJn5jMolHTsMoN7RJVx2aZuPHeutkJkHJtOWq+vuZ81ct2Zut69oQuzSWPBxWK2CJ0aJ8fGnggkUlu5n5xNY6w/Zt21BxLeDXoXTJtkqCtcsTK2UZ47t1jR0xIoV4RkZW2NPPfKNnvOAEcIwa/fuB0HJ5bx1CmnTrAzvpDBWH8Mw11h05ax13MHDMvHy+IJSrkDlTcv1rSZEVHtE6p+LBmgsmGHlb5ZUX6gcoUq+5157heZo9bK4G4Psokmg41f3FrsMQbiIcRDshXcmRjqj6tWUF3MFG0VIe2vDzftue8YjCNX1G21rADDAtzQHfbsfdsTVSER77kvZ4aQn/nCakws5az6TezzAfvf+dOPna76Hk1PqFLjFlj3sjFCyDsJIfsIIfsy6fIwjQ3z2dC8K+Id3NniJWbLAMBVY704ML7ke/gyk8xjg0stCS9G+6LennumgK6wgk09EWSLRlokS90aSISgyhI2dEdqLmQ6NZO2VDtgKPdqtkxfTMVQwswmSDZ28vN8ed85/N33j1aMBNgNppebUHWzZVSXlNLX7R1Bb0zFZx4/7fm5BU3HxHIOo30xbOiJmLaMMRJwplZ6lx8IRrkDQDQkIeeYULU3fS7bMn6KhjG6o3ZbhmVfeSt3+/FUZGJ57lsHjDkZZ/bNsiPIWrZMi9Mhy/n+lTc6Qgi2DMRxdp4pd+N8MrJlmHIvVoyQGAPxUNO2zI4h47py+u5TyzlPSwYwzr9ejxF00cyht0+osptpYyOlyaWcVeIDKAd3PiNvYql6HGn0Cphidov5/7S5/TyAMe55o+a2Ciil91NKr6WUXtvTbQTndF7jsmVMxRHxnpg4NmX45RdxwX3vll6kCyUcnqzdc5BS0zJJ+A/uY30xz4VMC+kCeuOqVXhoNpm3TkY2rBvprb6QiVKKU7P24N5fpWjSUqaAnmjI+syZlP/JZC8mzQnfZ87aR0DOJtKuee667mpNREMy7r1uC7734qTn8ZtYyoJSYKwvio3dEUwt51HUdRdbwjtbJlDl7ujGlM5riPHBPcRny9Sj3BVbTfesxz6zG5TzeLLvP5PMY7gr7JoD7Wx1l2gy2PilnBLqfqPb2h+zVnMyFW7kuRvXx2Km4NnLuJnKkCyGsOvKmTEztZz3XMDE6IupriNoNyuK3ZwaUe7LOaN67CbelonZlXuxpHuutGU0Gty/CeCt5s9vBfAAt/3XzKyZGwAscfaN906YQ6FUXrNO9LIt4+25H5tOYaQ3apt8een2AQDAz05WH/4DMDuaUwxw+am1GOuPYTZVcLU/5jNF9MdCGEywQJu3ncCAURe+Wq771HIe2WIJO4bKN6z+uJEG5TaJyJS79ZkB+O4XzODunABlqWq9XJ57sWRvgFAs0QqlyXjLy7aCEILP/eyM6+PM7hrti2FDdwRTSzkUNVoR3NUqnZ1yxQCVu8NzT5n9UxlhrrOQn0YdjO6ovaY7y55y5lizwBxxfBdZkrCQKSCv6RjqCrt2mko5Ak7M8txbPaFqv6k42ToQw7mFLHSugBfLlgGM4OU2oQoYnvt8urH6Mum8hrAiWVYHf51QSmsqd8DbFnJrTtKMDTZhZhNVs2Wmk3nUOgx+UiG/AOAJABcTQsYJIW8H8NcA7iSEHANwh/k7ADwI4CSA4wD+BcC7/HwZNvnG2zJsQjVRpdXe0akkdm9I2LZt7o1i60AMT5yo3bxi3qGq/WBNirqo70Uzrcum3M2TgQXfkd4oJpdyntkeJ2cMq2mHQ7kD5WwcRq5YQrZYQl88hGFLuTdvy0yawz1ncF9yKPeYS3EtraTbVpPyjPRG8arLN+CLT51zXa3K7K7Rvig2dEeQzGtYyhY9bQk38loJikR8TWzWwjkyyRRKNs89qsoomA08lrO1K0IyrPoy5nnN39R4fuGqzfj8O15aMbJUJGLdEIa6wqbN4zGh6kiFbIctE5Ilz5vrloEYCppRYmI+XYAiEXSZy+uBsufuFtwH42EUSpV+ua/9yhsTnsPmDZS3ZZJ5DZlCyXMBE6MvrnoEdzfl3nhwv2Bef5t5W4azrQBgwsdiSD/ZMm+ilG6ilKqU0lFK6ScppXOU0tsppbsopXdQSufN51JK6bsppTsppVdSSvf5+TIyU+45zSVbxt1zZ5kybpOhN+4cwJOn5mv67nOcH+4XdgG6ZczMp43Vdm7KnS2CGO2LQdMppjwU9slZY8jqtGWAyhxbFmx7oir64yEQ0rxyzxZKWMgUEVIkHBhfsq3GZR4p89yZjcD70ppObTnuTl51+UYsZYs47dIPc3whC1ki2NQTsXKixxcyFYFalb2zZdwKcDUKnw0D2Pt6sscB44ZSly3jqOk+vpBBSJaseRPr/UMybrxosOL1ikysi3uohi3DL2IC2pMK2RWpbNTBYP1Qz8xlrBx3QghUc+XuYqZY1ZYByqKsHtgNoyusIKJKNltmaqn6Aqby54ddm2Q75zeA8gKsRmwZZou6KXd2zlyosVYGWC0rVE3lni5oFf5jV1hBqqBVlItlmTK83864YccAkjkNL16oXkSM+eFsabMfxvrZQqbKO+dCulzhTiKGcp9PF9AVUazMl5Eqyh8w0iAjqmQbIvZ7rI6zhrWxEBRZwkC8+XRIVrDrtouHUCjpeOF8+RjyFSGBcsDIOIK7Wz1sBvt7sclwnvGFDDb1RKCYE8/GtqxVeIxRrT680YUpmNOanzB29vUEyudoplCqb0LVURlyfD6Lkb5ozdpGDEUiljhgnnuFcjdtCHbexUPNee7/+tOT+NjDx2paIs6UQCcsHfLsfBoL6aJ1bgPGeTWXziNX1G0jJEa/KcIaadqRMtMrCSEY6grb/GrWgamWLdMfNzx3Zyxyzm8AfOppI7ZMFhKBNRoHjEyqiCpZN/FAlHs7sHnuhRIkUvZMuyIqqEvPyaPW4qNK5f6yHYbv/sSJuaqfa/nhdSj3oUQYYUWqmBRkRcP64ypkiWAgYZQDmDNz3BnWQqbFDE7MpPDpx07ZVMSp2TS2DcTtq/MS7sqdzZwzm2QwEW46uLOT5jUv2QzAbs3wFSEBj+Be0m056U52DiVASLkGP8+5haxle7ELbWIpWzESMPK8vT33oJR7hPPc3fp6MuU+ny6gpNM6lLvTlslY39sPRiqkEWCGuiIeyt3eMEMyuzE1qtw//dhp/MNDR/GJH5+o+rxkrljVntrUY/RFODOXMcr9xsvP7YurluhJeKRSAo2tUuVHXcNdEZstM+nRO9VJfzyMkk4rnISqtkwDJQguLOUw3BWpGLH2RstlkSeWcq6jG55VEdzZd2C2TCxUHtaxP7JzeHN4MglCgF3DlcF9uDuCnUNxPF4juDMFUM+EKiHETIe03zlZoO0z34sF2rlU3ubps+D+Zw+8iNv//if4i28dxCd+VL5gTs2msXPIPhqxlHvaGdzZBKdxgQx1hZsuHsa83CtHerB1IOYI7gXbIg+3srhFlzx3nogqY0t/DMem3ZU7K7vMOtDr1GURT7XaMgEq9xhXO6dchrZ844iYNzcWKGoVDWNUKPeFbIXfXg02pxFSJHRHFHNRlL3XsFvZ3UZruidzRnP3/ngI/+t7R/CFp856Pnc5W31iWZEljPZFcWY+g4W0vdhWbzRkzT8kXFIhm6kMmcprVimG4a6w3ZbxHdyNv5uzI5PbhGpYkaBIpCFbZmIpa0uDZPA38YmlrC2bxo1VEdwlQkDMhrLZomars+3VsOPIZBLbqjS0vnHnIP779HzV5b9zqQLiIblupTfWH6vw3J3e+lBXGDOpglV6gBENybjpogHsGErgA6+9DDfuHMBDB6dAqZF1cnY+Y/PbgfINw5nrXlbu5mcGoNwnuWXY12zpw9NnFqygYZQeKF+MLF01a7NldM9sGcau4QSOO2yZvFbC1HLeCnKxkGIFCddsGY88d6NdXXCeO7tx8f1TGSyLhQUK/3nu5ZrumYKGuXShLuXORjLDXcaim56oipLZ5o+RzGkV6rfRmu6sON+HX3cFXnHxEP7n1w/g4UNTrs9N5jR0havf5LYMxHGWee68LRNTMWUeS7cJVb+VIQuajtOz9jmdtE25O22ZHLojStX6/kD52nbeXNyUOyGk4bK/E0s5bO6pPB/swT2HTb3Vz5lVEdwB46Jhs9ZR7uJkB8wZ3N0yZXhetnMAmUIJz1dp3jGfztdlyTCMXHe7cuf9bwAYTISsbBnnyOA/3nEDvvHum/D2m7fjnr0jOL+YxYsXlnF2PoOSTiuCuypL6IooFcNRlj3Txyn3mVS+qVZkFxaz6IupiIZkXLOtD7OpAs7OZ/Ds2QU8fmIOO4fK+xZ1Ue5aiVa1ZQBg53ACJ2dTtowhVkyKD3JMSVVky1TJcw/UczeDO6XU1oXJepwpdzNQ+LVlmGBZyhYtG6Ku4G4eX5aV5bZ60a2meizUWLBh9Zcu39yDT7z5auwa7sKHHzzk2jbRz9zDlv4oTs+msZApOpS7aqX3uQX3aEhGLCRXtWVKOsVv//vTeNX/fqRiMpzd7Ia7I0jmNOvxyaWcNVKsBru5OCdVk7kiIqpUIUIaKftLKcXEovv+9MTKNe8vLOZs5QncWDXBPWHe5VgXJkY5uNtrcZyeS+PiKsW+btjB8t29rZk5h6r2y1h/FEvZom1VIKuTwU5WFmgX0oWqN5DbLxmGRIDvvTiJUzNmpsxQvOJ5Ay4LOBazBYQUyQqyg4kwCpreVCsy40Q3As01W/sAAN9+fgLv/NzT2NAdxgdee5n13HLN8/LnFUvui5h4dg13oViiVvcaoJwGOdZftic2WsHdacsYPVndbmL5YmV1xUaJhGRQaqx6dWsgwRc4A9xXZbph1XTPatYIsB5bhq3WZdk15UwKu3J3KuhEWLFuUpmChs8/edZXA5Ujk8uIh2SM9EYRCyn4zZfvwImZNB49XpluvOzRHJtna38cybxRBItX7rzl5+Un11rI9HffP4KHD08jr+k2i5LPwGE3RWanTSXzNS0ZANb8gJty92oI7nUzTeaKrufvUraIbLHkarmwifO8VsJsKl/zhrRqgrtxIIxOJ1FbcDebZHMH6dhUCjqtXqa3Px7CJRu78HiVfPe5VAGDdfjtDCsdkgtO85ZyN1W0GWg1vfoiqYFEGNdt6zeCuzmU3DFYGdz73IJ72qjQyOYnrPz6Jnz3C0tlRbBruAtdYQX/63tHkCuU8Mm3XmfLuWZ/J14hGdkytW0ZADjO+e7jLgqWXXDOkQArb+CW657XSoFNqMbYyKRQwlOnDJHAsj2A8pxD2XP3Z8sAZk33XNH63iwLyw/Mc2c5292uyr0yayVmXmMA8NVnzuNPvn4A9z9SfYIUAI5MJbGb61T2mpdswmAijE8/dsr2PK1k1IaqdRy2cMfQ6bkz3LJlALMEgUdwf+C58/jnH5+wRvTMvmEt9th7siyU6aTR9ev8QtZXcPdW7u4ZQoYtUznHcWQyievu+wG+9kzl4n02gt3sYrkwW2bKLE/sZt3wrJrgnoioli0Tc/Hc+YVMR8yyA7trFPx62c4B7Du94Nk5ad4xoeOXMSu4l62ZRZYDznnujFp59K+6fCOOTqXww8PTRuOCWOXz+126rzs9S6sEQRO+++RS1lIEskSwd2sfJAJ87Ff3VhzvmOqVLVP9tNrpGtwzUCRiu8hYrnvIYbOw8r9u1kyuqAeaCgkYmVpfeOocbtw5YBtZWMGdee4+bRkA1sKj8YUswkpljns1mOc+lDCOlZct4/Tc45xy/+nRGQDAx3543Fo450a5U1n5bx9WZLz5pVvwoyMzliABygKs1nHgb5B9cbvnzvBS7gOJMOZcxMuxqST+6CvP4/rt/fjI668EYKQiA1zbPmbLdBnHbTqZx/2PnMRsKo9XXFy97SJgnA8RVarIWvMarbjZMrpO8SdfP4BcUbfiGA+rF+Om3HujKtKFklWbx23SlWf1BHdzCGN47uU/rJvnfmRyGSFFwraB6kPZG3cOIq/peObMYsVjrK5MQ567qbL4dMh5s2gYC0SD3MVay/p55eVGUc0nTs5V+O3l93CzZYq2C6LZEgRsAROvGj7wmkvx2bddj9suHq54vlsrumrlBxiJsIKR3qhVGwgATsyksLk3aisQttFDubPf3SZVcwEqd/Y+Dx2cwvnFLN50/RbH42xCtRHlrprKPYORvqjnoh832DGylLsj+2Y5V0S6ULIJDACIm6mQxZKOJ07M4Y5LNyCsSPifX3/Bc55mxuxU5ryxv/mGLVBlgs9yheCqlfvl2cLdIPk8917O1nIrHAZ42zLMivnHN+21bEU2gmUt9lgGDjtuPz02g48+dBR3XbERr7myZmVya3+d6028GoLHXeY4Pv/UWTx9ZgESMea3nEy4LGBisGud1c1yew7Pqgnu8ZBRZiBbsGfLxEMKCIFtyfGRqRR2DSdqersv3dEPWSJ4zMUbTOY1FEo6Bhvw3HuiKhJhxTapupAu2FSITbnXGB2M9sVw+WZj/qBqcM/Y62osZgqWDcR/ZqPBfdKlYcGuDV24ZZe7qgkrEghx2jK1lTtgLGZi6ZDJXBE/OTqDW3fbV2Nu8PDcrZ6sLn5xPkjlbgb3zzx+Gv3xkHUTdj4+tZxDWJHquql0Rw3Pvd40SICbUGWeu6Oo1GkPe4/ZBPvPLSKZ1/CGq0fw/rsuxRMn5/CVp8cxl8rj+fFFnOCUPFOXFzss0OGuCF77ks348r5z1nwY+/xannsspFjnKp/nzo9Y3fLcgbIt47wZHZ1KYkN3GBu6I9b1ZgX3HAvuxmf1x0JQJIIvPHUO3VEFf3XPFb5vrv0u5be9GoI7s2WmlnP4m+8cxo07B3Ddtn7XjmwTS1koEqm4MQPlEdphc4K7I1IhgXL1x0yhZA33AWPxRSKk2CZUj0wu42IfNdi7IyquGu3BYy6+eyN1ZRjlXPeycmcFvBhDNuVe+zNedflGANWDe0HTbSp5MVO0T0hFVSgSadhzZwuYag33GIQQs6a7PVtGrqHcASO4H59OoaRTfOfAJHJFHW+4etT2HK/gzmwJN889SOXOUj3PzGXwS9eMVtRbL8856L5z3BlMuZ+bz2CsjkwZwCgcBpRv5omQAomUl8GfskpY2LPJ2ITqI0dnIBFjZHvvdWO4dmsf/vArz+Oav/oBfvHjj+G1H3vUWkPBMmXcrrdfv3Eb0oUSvvGs4R2Xy/3WHsFsNdW7zXM3rx9VJhXHmsGuA2e+vpE9V+7r0BVRLM/duUZBkog1yv3I619SV1VYt8Y5XumfibBss2U+9O2DyJd03Pe6KzHSF3VdZTqxmMOG7khFDX+AD+7L6I4orhlFPKsnuJsnXtYxoQoYwWbfaSPfejFTwNRyvkJJeHHzRYPYf26xot71XAOrU3nG+u3pkKxOBqPHDLSAv+D+mpdsQkiWrG7nTvqsBRzGfhvHwm7LsJPWq5lILaoNCb1wFtcq6rprPXcnu4YTyGs6zi9k8dVnxrFjMI49ju/OvH+nzcNa+LmtYQhUuYfK73PvdWMVj0e4AOQ3x738fNVqzl6vcmfHg9kLkkTQzeVAs+C+1WFbsgbV3z84havGetETUyFJBB/9lT34rZfvxAdfexnue90VyBZL+MrT4wAMlTiYCLsGwKvGejHSG8WTZgOWevrIbhmI2TK9gHJwrxa03OrLlHSKY1Mp2w1oKFFuGO9sOQgAt186jHfcvB13XmYfjdViwKXlZdUJVbPDXCqv4cEDE/iNG7dh+2Acm3uimErmK+pfXaiyOIkF96NTKdcJVyerKrgbtox9QhUA3nHzDhw4v4SHDk6VlYTP4H7jRYPQKfAzx2rV+QZWp/KM9RkLmdjwcD5dsPmHLND6XSS1cyiBZz94p2uhKH4/2X5nCiUUSrpNuQPAbZcM4TsvTNgmuvzCJnNq1djgcTbscOvE5MYuM6PhR0em8eSpebz+6pGKofFgIgzZpcKjpdwdE6qU0pZ47jfs6LeVYGZIEin3uKxbuStWGmI9Oe5A2XPnayLxC1xOz6Yx0hutOA7Mcz48mbRZbWP9Mbzvrkvwtpu3480v3Yprtvbh80+eBaUUR6eSVbPSXjLaY60lWfbpuQPA227ajg/94uW2vzkLXl6ZMkB5XmmWS0c8N59BXtOxm9vPwUSYm1CtXKNw3+uuxJ9yab1+cWatFUs6ssWSRyqkgpJOkdd0HBhfAqXADTuNFO1NvZGKPrhA9Zx7ZlsVNL2mJQOsouAeDyvQdApNp7a7OQC8/uoRbBuI4R8eOmr5TX6D+94tvYiqcoXvzmbc6xmS8Yz2RZEplCx1sJgpVmS5DHaF6nr/aorFqdytxhmO3Orfu3M3QrKEDz94yPfnMiaWctYCJr84m2T7mVAFgIuGjL/fP/7wOADgdQ5LBjCC2OWbu63sJAYL9s5uTMUSBaXBNOoAjJtcRJXwtpu2ez6Hnat+c9wZ/PPrDe6KZDRN57OI+Jrup2bT2DZYORqIcUHz1l3uIgIA3vzSLTg5a+SxH51KVr3WXjLai7NmKYF6lPsVIz241zFBHVaMRUrVaqa4KXe37LnBrpCL517f6Mr182MhpPKalYHnrJvPw7fa2z++CAC4arQXQDmNkZ9UpZQaq1M9VHkPd87UWp0KrKLgblsc4gguiizhvXfsxuHJJP7lpyfRHVF8q8uwIuP67f14zKHcG6krw8MWST3w7AXktRJSec2qPcHYPdzlWrWyEZzKfcGReskY7orgXbddhIcOTuFxl4nkakxwC5j8Eg0pyDQwodoTUzFs1sJ52Y4Bq+aOk6+/6yb87u0X2baFzJtH0aHcc+YFF5QtM5AI4/k/exVeac6HuMEyZurJlAHsAbBeW+buPSP47ZfvtG1jyt2tkxeDiYeusIKrPOw/AHj1lZvQF1Px1985jFxRrzq/9ZLRHgDAgfNLlufuNRnqh96o6pkpA5QnEfmS0ayIIN9uczARrvDcAwnuCft1WC1DiK/pvv/cIrb0x6ybE5vXYnnt7D3zVVQ5n5GzyUf8W5XBPeYyLPuFqzZj13AC4wtZXLKxu67UsZsvGsTx6ZRVJxkwDmSsgboyjMs2d+P67f34zOOnrdSoPseN4iNvuBKfePPVDb2/E6dyZyqNn8RlvP3m7RjpjeJD3z7ou5cswGpa+LdkACPXPcutUDVsGX9/G2bNvOGaStXOkCVS8bdmNw+nLZNnXZgCUu5AZY69E0u512vLmB59RJUwWOe8z22XDOP/vXWHbRtbvTifLmA5p2HbgFtwN/b1ZTsHqpZljqgyfvnaMbx4wUi5q6bcrxgxgvvz44tYzhYRVeWq712L/kSo6ihouDuCHUNx2+rYo9MpjPXbO7INJsJYMldzplxsmUYpK24jliTz3hlCCa7V3v5zi7YbKpvX4vug1przUmTJKn7WUco9bgvulRenLBH83p27Afi3ZBg3mT42b83MpfJ1Nelw4203bcf5xSy+tO8cAHvOLsBqMAcTaLrCClSZWCthFzLuyh0wLs73v/oSHJ5M4i2ffBJ/+93D+Ob+C56LuRgT3AImv1RMqJYqe556ccXmHnSFFfz8Fd7K2A1283DmubOUzKCUux8ili3TmHIf7YvVJVQ83y+qYimrWYp2h0sJCxaAbtlde8HOr5qWCSHlm7AbPVEV2wfjeH58qa6a9l586O4r8EevuqTqc27dNYSfnZyzzuejk8mK0QXz5udSBavFXjM3HYbVj8G0U8q13L2V++nZDC4s5XCVOcphz4+HZJtyZ+tmvEaxQNnO8yPCVk1w54c1Xp7vz1++Eb/18p1447WVmQvVuGRjF/rjIVtKZKN1ZXjuvGwDRvui+ORPjWXYboE2KAgh6IuFLK/RWTTMyWuu3ITffPkOzKUKuP+Rk/jdLzyL//Gf+z0Xq2QLJSw6FjD5IRqSK5p11CocxnjPHbvwnffeUvdw2SvPPW/+HtQN1Q+RhpU7C+71HW8vmHI/adYnclPuV4704AOvvQxvuHqk5vttG4zjtouHsGs44TqS5mGTqsl87boytbh6Sx8u2+xdMwoAbtk1iFxRx9OnF1As6Tg5m6royMZGQ7OpvNViLwhY4GWB2NnxiocF9ydOGnGHzwYjhGBzb9Sm3A9OJCERVLVyme/uR7kH840DwK0JghNJInjfXdXv6l6vu3HnAB49NgtKKQghmE8XfNWTqIYsEfz6jdvwV/9lTF42kjNfD2whEwAsmf/3eAR3Qgjef9eleP9dl6Kg6bj/kRP4u+8fxd5He/GOW3ZUPL+RTBnAbEXHdSsq6f6yZQDDfqsVONxgNw9nnvtKKPeGJ1QjwQb37qiCQknHoYkkZInYyiQwZIng7Td7Tw47+f/etNe2QM2LK0d68MBzF3BiOh1YEK3GDTsGoMoEjxybxVBXGMUSrVTuXJ0lr56sjRAPK+iPl+vOJ11a7DGYaHn8+JyZHNBje3xTb9S2kOnghWXsGEpUTWhg6aJ+rtNVo9wTNWyZZrll1yCmk3mrNnWjdWWcvPG6McTN/e2LN6daasEvvV7IFBELyZ6LPXhCioR333YRXnnZBnzkO4fxlJmXzGP1bfS5gIkRC8nWhCqb4PST594M7ObhzHNnw/R2Knd2IboNy6vRa3bsci40ahSm6PaPL2KsLxqIBdEdUa06LNVgXvKRqWTdI5hGiIcVXL2lDz89NuNZZ4otIpxNFqwWe0Ex2hflgnvtCdWTs2ns3tBVEbQ390Rs2TKHJpZx2abqo5aeqOo7o21VBvd6UvH8wvJ6f3psBpRSzKUKTXvugHEBvPG6MURV2VbVrhWwHNtvPHse33j2fF0WCiEEf/fGqzDWF8XvfP6ZiuJLFxpYwATYUyFZaqJf5d4oqkeeuzWh2lbPvdE8dxVf+a2XWd52s7Dg/sL5Jc9Vzq3i8s3dYPfzdih3ALh19xBevLCMx08Yqtg5x8A3qXc2Nm8WI7gzW6bKhCo3Kt0z1lPx+KaeKGZTBeS1EpYyRrerWpbUG68bw+/evsvXfq6a4M6nPzUyVK/F5t4oLhpO4CdHZ5Ay68o0mgbp5H13XYJv/+7NNTMrmmUgHsKp2TTe+5/PYbQvio/du7eu13dHVHzizddgOpnHF//7nO2xR4/NoCeqVp3McSMakpHXdJR0WlbuPrNlGkX1yHPPrYByb3RCFQD2bukLTMiw4J7XdGxrc3CPhRSr3WWznrtfbjHz9L/2zDi2DsQq/ubRkIx4SMasGdy7AgzuI71RnF/IglKjnyrfiJwnxsU0lt/Ow0bJk0s5HJwwMpMuraHcb7t4GL9RZd0Fz+oJ7qHannuz3LprCE+dmrdmqJudUGWEFbmi72kreMloL4a6wrjvdVfga++6qeZd3g2WwvnVZ8atydV0XsP3XpwySiDUeYNiFlq2WLI6K/mdUG0U1SPPnSn3oJp1+KHRVMig4Re4uPUDaDVXmpkgzWbL+OXyzT3oi6lV8/CNnsKFQD13wMhwMpqBGGmnXjc0VS4Hfbd1BXxa5SEzuNeyZeph1QR3SSKWd90KWwYAbtltlAD+7guTAGrXWV9t/NI1o/jv/3kH3vzSra6FhfzyhqtHcHImbS0b//7BSWSLJbxub+0sCidWq71CyQq2rbZlrDx3h3Jnk83Vlq8HDVOM7VKsXvA3l3YrdwBWml+7bnKyRKwUZ6++DqwEQStsGcDImEnmilXnWxJhBVFVti2wYmw2lfvEUhYHJ5Yx1BV2rQbZKKsmuAPllW2tmFAFgBu2DyAkS/jGc0YVu6BsmU7jris3IaxI+NozRnGorz97AaN9UVyzpa/u94pyTbLZEmu3HOsg6U+EIBHg1Ky9SfmTJ+cx1BUOLAPFD9EmbJkg4ZV7uz13wBhVAu3z3AFjJA7UCO7Mlglwv9iK4vGFrGfRMEYirODKkR5XwVNeyJTDwQu1J1PrZVUF97i5UCeImX43oiEZ123vs4pqtTp1cbXSHVFx52Ub8M39F3B+MYtHj83gnj0jVhu1euCbZP/o8DS6wgqu29Yf9C7b6I6ouGqsF4+Y3YQAIw3z8RNzuHHnQCCLgvxy2yXDeMsNW1tmJfqFpWKGFKlm+7VWcMVID957xy688rL6FqQ1w6tfYqzl8OqiNNgVwuRyztZiLwhGLOWeNWu5e49W/sedu/GeO9wnQKMhGX0xFWfm0jg2nazpt9fLqgrubAjTSm7lquENBOS5dyJvuHoUC5ki/vDL+6FT4J69mxt6HzbKShc0/PDwNG7dPdSymzPPrbuG8Pz4olV3/Nh0CrOpPG7a6V0QqxVcs7UPf1lHs4dWIUsEXWEF2wZiDd2kg/j8996xu+4Vzs2QCCt4/12Xeloug4lwIPVu3D63N6aatkx15X7P3hHLPnJjU08Ujx6bRbFEG5pDq0ZTVyEh5DQh5AAh5DlCyD5zWz8h5CFCyDHzf99j/URYaZnfzmApkbGQ3PLPWs3csmsQg4kQHj8xhytGunHRcH0lHRjMc3769AKmk3ncdkllO75WcOvuIegUVo0RVlrixosG2vL5q5H+RAg7AsqbXwvwrS4TVYqRNcJoXxTnF2vbMrXY3Bux0pBXoy1zG6V0D6X0WvP39wF4mFK6C8DD5u++SISVlk+GXbqpC4OJ8Lq1ZBiKLOHuPcYE6j176p9IZTDl/u0DEyAEvhoNB8FVoz3ojiiWNfPY8TlsHYjVXWFxLfHRX9nT0ArutYo9uAc70TvaG/Nly9SC+e4RVQp8rqQVkfRuAK8wf/4sgB8D+GM/L/ztV+ysaD4bNIQQ/D83bGn553QCb33ZNpyZy+D1LrXU/cKC+/5zi9gz1mu7oFqJIku4edcgHjk6C62k48mTc3jtVf6aHK9Vrm5gQnwtM9RVFnDVygg3wmhfFD86YjTlbk65G8H9ko3dTWXAudFscKcAvk8IoQD+L6X0fgAbKKUT5uOTAFz7WBFC3gngnQCwZYuxSm9vm07O996xuy2fs9rZMhDDv7712tpPrAK/eOT2NlkyjFt3DeHBA5P42rPnkcxruLHNfrtgdcMLjaCzeEb7olahumaUO0uHDHoyFWg+uN9MKT1PCBkG8BAh5DD/IKWUmoG/AvNGcD8AXHvttf6LjgtWFXzaarv8dsatZunajz50FABw487167cLKuGDe5B57gAwwtl/zdw4mC0T9GQq0KTnTik9b/4/DeDrAK4HMEUI2QQA5v/Tze6kYPXCJqU3dIdxeQtO0GqwkhITSzlcsrGr4ZaJgrVJnMu+C7JwGGCv5llv0Tieq8Z68Bs3bcOr6+xp4IeGgzshJE4I6WI/A3glgBcAfBPAW82nvRXAA83upGD1ElFkqDLBz12yYUXSAVlqa7V0M8H6ZdD03YMO7iNccG/GlgkrMv7sFy5viTBpRrlvAPAoIWQ/gKcA/Bel9LsA/hrAnYSQYwDuMH8XrFEkieAzv3E9/uCVKzOPcfulhhX0ch/dhQTrD2bNBG3LdEdUa1VwO1fl1kPDe0UpPQngKpftcwBub2anBJ3FSqrmmy4axHfecwsuqbP1omB9MJgIB9Ziz8loXxRL2ea7T7WKVbVCVSBohEs31dcwXbB+2D4Yx6YWrZpl5bHXnHIXCASC1c5779iFd9TRWrAe2II5EdwFAoGgzTTap9cPv3LdGIa7w75aXa4EIrgLBAJBA1y8sQsXr+K5HuG5CwQCwRpEBHeBQCBYg4jgLhAIBGsQEdwFAoFgDSKCu0AgEKxBRHAXCASCNYgI7gKBQLAGEcFdIBAI1iCE0pXvk0EImQFwpsbTegAsNfD2jb5upT5zEMBsGz9THJ/WvVYcn9X3mWvt+GyllLqXRKWUdsQ/APe383Ur+Jn7Omhf1/zxWYnvKY5PSz9zXRwfSmlH2TLfavPrVuozG0Ucn9Z9Zru/pzg+rX1tuz9vRc6DVWHLCMoQQvZRSpvrWr2GEcenOuL4VGc9HZ9OUu7rhftXegdWOeL4VEccn+qsm+MjlLtAIBCsQYRyFwgEgjWICO4CgUCwBhHBvcUQQj5FCJkmhLzAbbuKEPIEIeQAIeRbhJBuc/udhJCnze1PE0J+jnvNNeb244SQj5E10jS0nuPDPb6FEJIihPwBt+20+fznCCH72vkdWkmd58+bze/P/umEkD3mY+v+/CGEhAghnza37yeEvIJ7zY8JIUe4Yzfc/m8TMM3kUYp/vnJVbwVwNYAXuG3/DeDl5s9vA/CX5s97AWw2f74CwHnuNU8BuAEAAfAdAHet9Hdr9/HhHv8KgC8D+ANu22kAgyv9fVbD8TG3XwnghDh/bNfXuwF82vx5GMDTACTz9x8DuHalv0+Q/4RybzGU0kcAzDs27wbwiPnzQwDeYD73WUrpBXP7iwCihJAwIWQTgG5K6c+ocSb+G4B7Wr7zbaCe4wMAhJB7AJyCcXzWPPUeH443AfgiAIjzxzo+lwH4ofm6aQCLANZsWqQI7ivDiwDuNn/+ZQBjLs95A4BnKKV5ACMAxrnHxs1taxXX40MISQD4YwB/4fIaCuD7pp31zrbs5crh5/z5FQBfMH8W54/BfgC/SAhRCCHbAVwD+7H7tGnJfGAt2FYiuK8MbwPwLkLI0wC6ABT4BwkhlwP4GwC/uQL7thrwOj5/DuCjlNKUy2tuppReDeAuAO8mhNzalj1dGWqdPy8FkKGUvuD24nWA1/H5FIwb2z4A/xvA4wBK5mNvppReCeAW899b2rnDrUBZ6R1Yj1BKDwN4JQAQQnYDeA17jBAyCuDrAH6NUnrC3HwewCj3FqPmtjVJlePzUgC/RAj5WwC9AHRCSI5S+nFK6XnztdOEkK8DuB7lofmaotr5Y3IvyqodEOfPa8ztGoDfY88jhDwO4Kj5GDt/koSQz8M4f/6tvXseLEK5rwBsJp4QIgH4UwD/x/y9F8B/AXgfpfQx9nxK6QSAZULIDeZw8dcAPNDu/W4XXseHUnoLpXQbpXQbDOX1YUrpxwkhcUJIl/maOIwLe82qVq/jw217I0y/HRDnD8rXV8w8P0AIuROARik9aNo0g+Z2FcBrsQbOH6HcWwwh5AsAXgFgkBAyDuDPACQIIe82n/I1AJ82f/4dABcB+CAh5IPmtleakz/vAvAZAFEY2Q7facsXaDF1Hh8vNgD4ummTKgA+Tyn9bmv2uL00cHxuBXCOUnrS8Vbi/DEyZL5HCNFhjFyY9RI2t6sAZAA/APAv7fkGrUOUHxAIBII1iLBlBAKBYA0igrtAIBCsQURwFwgEgjWICO4CgUCwBhHBXSAQCNYgIrgLBALBGkQEd4FAIFiD/P9vlAYq1mJCoQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"extremes = pd.Series(\n",
" data=[\n",
" 129.7, 51.4, 148.8, 76.3, 68.3, 91.8, 116.5, 97.4, 81.3, 65.8, 50.9, 100.2,\n",
" 123.2, 199.4, 137.7, 104.7, 72.6, 63.2, 116.5, 71.9, 121.6, 87.8, 89.2, 76.2,\n",
" 103.7, 112.0, 60.7, 63.2, 119.7, 75.5, 68.3, 67.6, 62.3, 51.1, 62.5, 99.4,\n",
" 81.4, 86.3, 143.8, 44.5, 58.4, 270.4, 103.5, 94.0, 68.9, 80.0, 112.5, 73.4,\n",
" 262.8, 58.1, 43.3, 59.0, 181.8, 117.1, 109.4, 115.9, 116.2, 175.8, 73.4, 83.7,\n",
" 165.6, 56.8, 123.4, 142.2, 195.2, 130.1, 64.5, 97.1, 219.4, 95.7, 319.4, 72.5,\n",
" 115.5, 53.0, 81.0, 152.0, 68.5, 83.5, 77.0, 77.5, 67.5, 91.0, 86.5, 77.0,\n",
" 57.0, 125.5, 197.0, 105.5, 75.5, 128.0, 70.0, 150.0, 121.5, 72.0, 57.0, 81.0,\n",
" 86.0, 40.5, 118.5, 76.5, 82.5, 122.0, 179.5, 123.0, 71.5, 38.5, 108.5, 138.0,\n",
" 57.0, 76.5, 83.0, 47.0, 69.0, 105.5, 124.5, 64.0, 100.5],\n",
" index=[pd.to_datetime(f\"{i:d}-01-01\") for i in range(1897, 2014)],\n",
" name=\"precipitation\",\n",
")\n",
"extremes.plot()"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "0c6749bc-1a39-4854-b1b1-75ef4d8fb4bd",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
" Univariate Extreme Value Analysis \n",
"========================================================================================\n",
" Source Data \n",
"----------------------------------------------------------------------------------------\n",
"Data label: precipitation Size: 117\n",
"Start: January 1897 End: January 2013\n",
"========================================================================================\n",
" Extreme Values \n",
"----------------------------------------------------------------------------------------\n",
"Count: 117 Extraction method: BM\n",
"Type: high Block size: 365 days 05:49:12\n",
"========================================================================================\n",
" Model \n",
"----------------------------------------------------------------------------------------\n",
"Model: MLE Distribution: genextreme\n",
"Log-likelihood: -588.265 AIC: 1182.743\n",
"----------------------------------------------------------------------------------------\n",
"Free parameters: c=-0.197 Fixed parameters: All parameters are free\n",
" loc=77.842 \n",
" scale=28.110 \n",
"========================================================================================"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"model = pyextremes.EVA.from_extremes(\n",
" extremes=extremes,\n",
" method=\"BM\",\n",
" extremes_type=\"high\",\n",
" block_size=\"365.2425D\",\n",
")\n",
"model.fit_model(model=\"MLE\", distribution=\"genextreme\")\n",
"model"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"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.5"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment