{ "metadata": { "name": "", "signature": "sha256:6fc408129f5b4b64743709010c48b9bc4b997eeaf88da894bb35e6808e5ea09d" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "In this Notebook the findpeaks.py script is illustrated with a time series of daily discharge. \n", "\n", "This function finds peaks in a time series, e.g. a hydrograph. A peak is a data \n", "sample that is larger than its two neighboring samples AND the largest sample \n", "within a moving time window.\n", "\n", "Three methods are available: \n", "peaks over threshold, annual maxima, selection of n peaks" ] }, { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline\n", "from hydrotools.statistics import findpeaks as fp\n", "import pandas as pd\n", "import matplotlib.pyplot as plt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 34 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's load some sample data\n", "try to load the data using a relative path \n", "OR edit the directories to subdirectory data to find the hydrograph data" ] }, { "cell_type": "code", "collapsed": false, "input": [ "try:\n", " fname = r\"data\\hydrograph.xls\"\n", " sample = pd.read_excel(fname,index_col=0)\n", "except: \n", " OETdir = \"d:\\\\scripts\\\\OET\\\\\"\n", " DATAdir = \"python\\\\applications\\\\hydrotools\\\\notebooks\\\\statistics\\\\data\\\\\"\n", " fname = \"hydrograph.xls\"\n", " sample = pd.read_excel(OETdir + DATAdir + fname,index_col=0)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 35 }, { "cell_type": "markdown", "metadata": {}, "source": [ "First, we will find all peaks given a 30 days period between peaks and a minumum value (threshold). The distance betweeen peaks is important in order to be able to treat them as independent events. This yields only the largest peaks" ] }, { "cell_type": "code", "collapsed": false, "input": [ "peak_dist = 30\n", "threshold = 200 #m3/s\n", "POT2 = fp.peaks_over_threshold(sample, threshold=threshold, peak_dist=peak_dist)\n", "\n", "fig = plt.figure()\n", "ax = fig.add_subplot(111)\n", "sample.plot(ax=ax,style='-b', legend=False)\n", "POT2.plot(ax=ax,style='og',markersize=6, legend=False)\n", "ax.legend(['Discharge','peaks above threshold'],fontsize=10)\n", "fig.show()" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEPCAYAAACwWiQoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnWmYFNXVgN8GWRSJ4IaI6CAqhkVZBERFGgUHMyxGRBQ1\njBr9FFnUqIAxLokaNUYFoiYxChojRsHgMooQpFVcgwIqyKqDIvsiyDBsM/f7cbumq3uquqt6qruq\np8/7PP3U0rdunVN169xTp+4CgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIWaUNsMD02w6M\nBg4FZgPLgVlAE9Mx44EVwFLgvGwKKwiCIKSmDrAOaAk8BNwW3T8WeCC63hZYCNQDCoCV0eMEQRCE\ngHAe8H50fSnQLLp+VHQbtFc+1nTMTOD0rEgnCIKQp7j1mC8BpkbXmwEbousbiBn2o4E1pmPWAC3S\nFVAQBEFIjRtjXh8YALxs8Z+K/uxI9p8gCIJQQw5wkfZ84DNgU3R7Azq8sh5oDmyM7v8BHVM3OCa6\nr4qjjz5arV27Nh15BUEQ8plVwAlWf7jxzC8lFmIBeA0YHl0fDsww7b8E7cm3Ak4EPjVntHbtWpRS\nKKXo1atX1bqT31133ZXT6UXf3JI/k/oGTXa5t8HXF2htZ6CdeuaNgD7ANaZ9DwAvAVcDpcDF0f1L\novuXAPuBESQJszRs2NChCJpwOJzT6UVfb/MPWno3+gZNdrm3yQmavkFBGQwfPlzlE6Jv7Saf9M0n\nXZUKhr4kcYx9b/9dXFzstwhZRfSt3eSTvvmkKwRf35BP541WMoJQezj00EPZtm2b32IItYSmTZuy\ndevWuH2hUAhs7LbvnnkkEvFbhKwi+tZetm3b5uqDl/zkl+zn1jHw3ZgLgiAINUfCLILgEaFQCCnX\ngldYladAh1kEQfCGunXr0qlTJ9q3b0/Hjh155JFHqozBZ599xpgxY1znGQ6H+eyzz7wWVcgAvhvz\nfIqpgugrZI6DDjqIBQsW8NVXXzF79mzeeust7rnnHgC6dOnChAkTXOcZ9QTTYv/+/WkfK7jHd2Mu\nCIL3HHHEEfz973/nL3/5C6Ar1QEDBgDw7rvv0qlTJzp16kTnzp0pKysD4MEHH+SUU06hY8eO3H77\n7VV5vfzyy3Tv3p02bdowb948AEpLSzn77LPp0qULXbp04aOPPqo6T8+ePRk0aBDt27dHKcWIESP4\n+c9/znnnnUdRURHTp08H9NtCOBzmtNNOo1+/fqxfvz5r10fwDiUItQ2/y/XBBx9cbV+TJk3Uxo0b\n1dy5c1X//v2VUkoNGDBAffjhh0oppcrKytT+/fvVm2++qc444wxVXl6ulFJq27ZtSimlwuGwuuWW\nW5RSSr355puqT58+Simldu3apXbv3q2UUmr58uXqtNNOU0opNXfuXNWoUSNVWlqqlFLq5ZdfVr/4\nxS+UUkqtX79eNW3aVE2fPl3t3btX9ejRQ23evFkppdSLL76orrrqKu8vSg5jVZ5I0mnIzUBbgiCk\nSQ2iFXj9TfXMM8/kpptu4rLLLuPCCy+kRYsWzJkzh6uuuqqqy3qTJrGJwy688EIAOnfuTGlpKQB7\n9+5l5MiRLFq0iLp167JixYqq9N26deO4444D4IMPPuDii/VIH82aNaN3794ALFu2jMWLF9OnTx8A\nKioqOProo71VNM/wPcySbzFV0Tc/USr9X7p888031K1blyOOOCJu/9ixY3n66acpLy/nzDPPZNmy\nZVEZrU/WoEEDQH9gNeLgjz76KM2bN+eLL75g/vz57Nmzpyp9o0aNEnS3zrddu3YsWLCABQsW8MUX\nXzBz5sz0FBWAABhzQRC8Z9OmTVx33XWMGjWq2n+rVq2iXbt23HbbbXTt2pVly5bRt29fJk+eTHl5\nOUDKDis7duzgqKOOAuC5556joqLCMt2ZZ57J9OnTUUqxYcOGqsq9TZs2bNq0iY8//hiAffv2sWTJ\nknTVFQhAmCVXRwpLF9FXyBTl5eV06tSJffv2ccABB/CrX/2Km2++GdCtUoyWKRMmTGDu3LnUqVOH\n9u3bc/7551OvXj0WLlzIaaedRv369SkqKuLee++tdg4jjxEjRjB48GCee+45+vXrx8EHH1wtDcDg\nwYOZM2cObdu2pWXLlnTu3JlDDjmEevXqMW3aNEaPHs327dvZv38/N910E23bts3kJarVSKchQfAI\n6TRkTVlZGY0aNWLLli10796dDz/8kCOPPNJvsQJPznUayreYqugr5Bv9+/enU6dOnH322dx5551i\nyDOE72EWQRBqN3PnzvVbhLxAwiyC4BESZhG8JOfCLIIgCELN8d2Y51tMVfQVBCET+G7MBUEQhJrj\nuzHPt3bIoq+QCxQXF1cNiOUFU6ZMsezAlCnuv//+qvXS0lI6dOjg+Tnuvvtu/vznP7s6xtwe34wX\n19t3Yy4IQaNkdgmFVxYSLg5TeGUhJbNLsp6XlzKkQ02Gvs1Gfqn44x//6PoYu16sdqSjk90x5k5d\n6eK7Mc+3mKroG2xKZpcw5vExzCqYxbut3mVWwSzGPD4mLWOabl7pHldaWsrJJ5/M5ZdfTtu2bRky\nZEhV93y74WafeuopunXrRseOHbnooouq0kPM8Pzud7/jyiuvpLKyknHjxtGuXTtOPfVUbr311moy\nfPrpp5xxxhl07tyZM888k+XLl1f99/3339O7d29OOukkfv/731ftf+SRR+jQoQMdOnSoGnN93Lhx\nPPHEE1VpzF7wn/70J7p168app57K3XffXU2GcePGVfWGveKKKwiFQlRUVHDttdfSvn17CgsL2b17\nN6DfHG+66Sa6du3KxIkTba/TxIkTq/QeNmxY1bmWLFlC7969ad26NZMmTUqqkxmlFCNHjuTkk0+m\nb9++bNy4MWstoZoA04CvgSVAd+BQYDawHJgVTWMwHlgBLAXOs8ivakjHuXPnZnYcyYAh+gab84rP\nU9xNtV/hlYUpjyVhyNJ080r3uG+//VaFQqGq4W2vuuoq9fDDD6t9+/bZDje7ZcuWquPvuOMONWnS\nJKWUUsXFxWratGnqlltuUddff71SSqnNmzerNm3aVKXfvn17NRl27Nih9u/fr5RSavbs2Wrw4MFK\nKaUmT56smjdvrrZu3arKy8tV+/bt1fz589X8+fNVhw4d1K5du9TOnTtVu3bt1IIFC9SCBQtUr169\nqvJt27atWrNmjXr77bfVtddeq5RSqqKiQvXv31+999571eQwDwf87bffqgMOOEAtWrRIKaXUxRdf\nrJ5//nmllB7i94YbblBKqaTX6eijj1Z79+6N0/uuu+5SZ5xxhtq7d6/avHmzOuyww9T+/fstdVq4\ncGGcXNOnT1d9+/ZVlZWVau3atapJkyZq+vTpcToklidjn52RdtppaALwJnBR9JhGwG+jxvwhYCww\nLvprCwyNLlsA/wVOAiqtMs63mKroG2z2qD2W+3dX7s5aXjWRoWXLlvTo0QOAyy+/nIkTJ9KvXz/b\n4Wa//PJL7rjjDrZv387OnTvp168foD3HP/zhD3Tv3p2//e1vgB4Wt2HDhlx99dX079+f/v37Vzv/\njz/+yK9+9StWrlxJKBSKm23ovPPOo2nTpoAeVnfevHmEQiEuvPBCDjzwwKr977//PqNGjWLjxo2s\nW7eOjRs30rRpU1q0aMGjjz7KrFmz6NSpE6CHCli5ciU9e/ZMel1atWrFKaecAuhZl4yhfAGGDh0K\nwNKlS22v0ymnnMKwYcO44IILuOCCCwD95tK/f3/q1avHYYcdxpFHHsn69euZN29eNZ3ee+89Tj31\n1KpzvvfeewwbNoxQKETz5s0555xzksrvBCfG/BCgJzA8ur0f2A4MBHpF9z0LRNDGfBAwFdgHlAIr\ngW7AxzWWVhAyTINQA8v9Des0zFpeNZHBHHdVSlV1PGnXrh0ffvhhtfTFxcW89tprdOjQgWeffbYq\nLBYKhejatSufffYZ27Zto2nTptStW5dPP/2UOXPmMG3aNP7yl78wZ86cuPx+97vfce655/Kf//yH\n1atX21bmhmzGutX+IUOGMG3aNNavX88ll1xSlWb8+PFce+21Ka+FGWMYX9BD+RphFogN2ZvsOpWU\nlPDee+/x+uuvc9999/Hll18CUL9+/bh89+/fX62zj1kng0x0MHMSM28FbAImA58DT6E982bAhmia\nDdFtgKOBNabj16A9dEtyLaZaU0TfYDN62GhaL2gdt6/1560Zdan7lhjp5lUTGb777ruqYWVfeOEF\nevbsmXS42Z07d3LUUUexb98+nn/++Tij069fP8aNG0dRURE7d+6krKyMH3/8kfPPP59HHnmERYsW\nVTv/jh07qrzZyZMnx/03e/Zstm3bRnl5Oa+++ipnnXUWPXv2ZMaMGZSXl1NWVsaMGTOqvOyhQ4cy\ndepUpk2bxpAhQwAoLCzkmWeeqZrq7ocffmDTpk3V5KhXr17SOUgTjS3YD8urlOK7774jHA7zwAMP\nVL3FWBnjUCiUVCeDs88+m3//+99UVlaybt06T4Y8cOKZHwB0BkYC/wMeQ3vgZpLGcqz+Ky4upqCg\ngNLSUhYuXEjHjh2ranHDAMi2bGd7u6hvEV8s+oLbH/gPZxQeROMGDenVsxeN6sUmXLA7PpGivkUA\nTJo6id2Vu2lYpyGjRo6q2m9HuseBNkiPP/44V111Fe3ateP6669POtysEUo54ogj6N69Ozt37qzK\nKxQKMXjwYH766ScGDhzICy+8wKBBg9i9ezdKKR599NFq57/tttsYPnw49957L0VFRVWVQygUolu3\nbgwePJg1a9ZwxRVX0LlzZ0Dbgm7dugFwzTXXVIUj2rZty86dOznmmGNo1kz7in379uXrr7+uCiU1\nbtyY559/vtoEHNdeey2nnHIKXbp04d5777X0jBPX69evb3mdTjrpJK644gq2b9+OUooxY8ZwyCGH\n2LZA6dSpk61ORvpf/vKXvPPOO7Rt25Zjjz2WM844w/J+RiIRpkyZAkBBQYFlmio9kv6rOQr4CO2h\nA5yF/sB5PNAbWA80B+YCJxMz9A9ElzOBu4BPTHkqr18xBMFLQiHYvBkOO8zNMf6OzVJaWsqAAQOq\nQgBCbpOJsVnWA9+jP2IC9AEWA68Ti6MPB2ZE118DLgHqoyuAE4FPnSogCEL6ZLs9txAcnLYzHwX8\nC1gEnALch/a8+6KbJp5DzBNfArwUXb4FjCBJCCbXYqo1RfQVMkVBQQFffPGF32IIPuG0aeIioKvF\n/j426e+P/gQhZ5FIoJBLyHjmgmBBKASbNsHhh7s5RsYzF7xDxjMXBEHIQ3w35vkWUxV9cwe3TnbT\npk2rmqvJT341/Rm9ZZ3iuzEXhNrC1q1bmTt3LkqpvPjlk65+6Lt161ZX5U9i5oJgQSgEGzaATCQv\nBIlQSGLmgiAItRrfjXkux1TTQfSt3eSTvvmkKwRfX9+NuSAEFYkECrmExMwFwYJQCNavh2bNUqcV\nhGwhMXNBSAPxN4RcwndjHvQ4lNeIvrWbfNI3n3SF4OvruzEXhKASkgEIhRxCYuaCYEEoBOvWwVFH\n+S2JIMSQmLkgCEItx3djHvQ4lNeIvrlDOi+PuayvW/JJVwi+vr4bc0EQBKHmSMxcECwIhWDtWmje\n3G9JBCGGxMwFIQ3E3xByCd+NedDjUF4j+tZu8knffNIVgq+v78ZcEIKKeOZCLiExc0GwIBSCNWug\nRQu/JRGEGBIzFwRBqOX4bsyDHofyGtE3d5B25snJJ10h+Po6NealwBfAAuDT6L5DgdnAcmAW0MSU\nfjywAlgKnOeFoIK/lMwuofDKQsLFYQqvLKRkdonfIgmCYMJpzPxboAtgnmH0IWBzdDkWaAqMA9oC\nLwBdgRbAf4GTgErTsRIzzyFKZpcw5vExrOq0qmpf6wWtmXDDBIr6FvkoWeYIheC776BlS78lEYQY\nXsXMEzMYCDwbXX8WuCC6PgiYCuxDe/QrgW4uziMEjIkvTIwz5ACrOq1i0tRJPkkkCEIiTo25QnvY\n84FrovuaARui6xui2wBHA2tMx65Be+iWBD0O5TW5qO8etcdy/+7K3SmPzUV9a0I+6ZtPukLw9T3A\nYbozgXXAEeg4+dKE/1X0Z0e1/4qLiykoKKC0tJSFCxfSsWNHwuEwELtosh2M7bL1ZQDQSi/4Vi8a\n1mno2fk+mv8RkcUR9qg9lK0v48I+FzL+lvG+6q9Uds+Xa9sGQZGnNuobiUSYMmUKAAUFBSQjnXbm\ndwE70R56GFgPNAfmAiej4+YAD0SXM6PHfGLKQ2LmOYRlzPzz1kwY6U3MPIgx+VAIVq+GY4/15fSC\nYElNY+YHAY2j643QrVO+BF4Dhkf3DwdmRNdfAy4B6qN9uROJtYARcpCivkVMuGECnecXwuReFK4u\n9MyQg8TkBcELnBjzZsD7wEK0d/0GuiniA0BfdNPEc4h54kuAl6LLt4ARJAnBJL7C1HZyVd+ivkXc\n838zYXWEmc/MdGzInehbk5h80MjV+5sO+aQrBF9fJzHzb4GOFvu3An1sjrk/+hOElDQINbDcb8Tk\nBUFIjYzNIjjmjTdgwADvB6DKdEw+HUIhKC2F447z5fSCYEmymLnT1iyCkDEMg/2HJyfxyee7KTyn\nIaNGjqq1HZIEIRPI2CxZRvS1pqhvEY/e6j4mn0lkbJbk5JOuEHx9fTfmgmAQ8ivoJwi1AImZC47J\nVMzc4NZb4eGHgzEpRCgE33wDrVqlTisI2ULGMxdygsWL/ZYgniBUKoLgFN+NedDjUF4j+tpTx/fS\nWHPy6f7mk64QfH1rweMj1BaCFjMXz1zIJSRmLjgm0zHzAQP0OYJQNEIhWLECTjjBb0kEIYbEzAVP\nyLSRDZpnLgi5hO/GPOhxKK8Rfe0JmjGXdubJySddIfj6+m7MBcGgNhhzQfALiZkLjnn9dRg4MHNG\n7pe/hBkzgmFEQyFYtgxOOslvSQQhhsTMhZxAPHNBSB/fjXnQ41BeI/raUxuMeT7d33zSFYKvr+/G\nXBAMgmbMBSGXkJi54JhMx8yHDIFp04IR3giF4Ouv4eST/ZZEEGJIzFzICXr29FuCeIJQqQiCU3w3\n5kGPQ3lNLuub6Rjy8ce7zz9o5PL9dUs+6QrB19d3Yy4IQUU8cyGXkJi54JjXXoNBgzJn5DI99osb\nQiH46ito185vSQQhhsTMhZxAWrMIQvr4bsyDHofymlzWN9Mx86AZc7f6lswuoWu/roSLwxReWUjJ\n7JLMCBYQcrksp0PQ9T3AYbq6wHxgDTAAOBT4N3AcUApcDPwYTTseuAqoAEYDs7wTVxCyhxtjXjK7\nhDGPj2HVUasgOtXcqsdXAQRicmqh9uPUF7oZ6AI0BgYCDwGbo8uxQFNgHNAWeAHoCrQA/gucBFQm\n5Ccx8xzk1VfhggsyF9N+6y34xS+CEzP/4gvo0MFZ+sIrC5lVUN1vKVxdyMxnZnosnZCv1DRmfgzw\nC+AfpkwGAs9G158FLoiuDwKmAvvQHvtKoFsaMgt5SC6HWfaoPZb7d1fu9kgaQUiOE2P+KHAr8d51\nM2BDdH1DdBvgaHQoxmAN2kO3JehxKK8RfXMHN8a8QaiBXvk2fn/DOg29Eyhg5PK9TYeg65sqZt4f\n2AgsAMI2aVT0Z4flf8XFxRQUFFBaWsrChQvp2LEj4bA+hXHRZDtY20YRyNX8MylPuF2YVfNWsarJ\nqiqD3npba0aNHBUYfby/PgRKntqobyQSYcqUKQAUFBSQjFQvtvcDVwD7gYbAz4BX0DHxMLAeaA7M\nBU5Gx80BHoguZwJ3AZ8k5Csx8xwk0zHzt9+Gfv2CEzNfsAA6dnR+TMnsEvoXT4J6uyk8pyGjLh0l\nHz8FT0kWM3cTpewF3IJuzfIQsAV4EG3AmxD/AbQbsQ+gJ1DdOxdjnoNky5hXVvofPw+F4PPPoVMn\nd8fVravll+ItZAIvOw0ZRfQBoC+wHDiHmCe+BHgpunwLGEHyEEy1V5jaTi7rm63xvXPZECoV8VuE\nrJHLZTkdgq6v03bmAO9GfwBbgT426e6P/gQhp8nlSkXIP2RsFsExM2boeTozHWbZv1+HK/wkFIL5\n86FLF3fHSZhFyCQyNouQExgGMJcNod+xfiF/8d2YBz0O5TWib2qCYszTkUNi5rWXoOvruzEXcods\nGdmgGHNByCUkZi445j//gQsvzJyxnTkTzj8fdu+GBg0ycw6nhELwv//Baae5O65ePR3zl+ItZAKJ\nmQs5QdAMYNDkEYRk+G7Mgx6H8hrRNzVBMaLpyFFZGfFcjqAiZTlY+G7MBSGRoBhzQcglJGYuOCbT\nMfM334SiIigrg4MOysw5nBIKwccfQ/fu7o6rXx/27ZMKScgMEjMXcoqgGMKgyCEITvDdmAc9DuU1\ntUHfNWtSpzHI5Zh5Okg789pL0PX13ZgLuYNhZFu29FeObJHLlYqQf0jMXHDMK6/A4MF6PRO3r6QE\n+veHHTugcWPv83dDKAQffABnnOHuuIYNYc8eqQiEzCAxcyGnEEMoCO7x3ZgHPQ7lNaJvaoJizKWd\neXKkLAcL3425IBgEbdTEoMghCE7w3Zgbk5jmC6Jvaj5JnDE2h6hbN+y3CFlDynKwcDPTkFDLKJld\nwsQXJrJH7aFBqAGjh40OxATEQ4boj6B+I565kEv47pkHPQ7lNUHRt2R2CWMeH8Osglm82+pdZhXM\nYszjYyiZXeLpeTIdMy+ZXULhlYWEi8MUXlnoqfwSM09OUMpytgi6vuKZ5ykTX5jIqk6r4vat6rSK\nSVMn2XrnmfZUjfwrK52lNyoksx6rHtfrQXjDEIRs4rtnHvQ4lNcERd89ao/l/t2Vuz09Tzr6Oq00\nklVIXpBO5SUx89pL0PUVzzyLBClG3SBkPftDwzoNsyxJdZwa0UxXSBIzF3IJ3z3zoMehvKIqRq0y\nG6N2yuhho2m9oHXcvtaft2bUpaM8PU8693e3Q1scxApJYua1l6Drm8qYNwQ+ARYCS4A/RvcfCswG\nlgOzgCamY8YDK4ClwHmpBPho/kcZ+4AVJDIdEnBLUd8iJtwwgZ+9WAiTe1G4upAJIyf4Gmt26wln\nq0IShFwgVZhlN9Ab2BVNOw84CxiINuYPAWOBcdFfW2BodNkC+C9wEmD5SatkdglPz3s6Lz5gVYUE\nWsXv9zpG7YaivkWc2KiIz5bCzGdSp08n7JDJOKNRRvoXT+K0M3ZzWOOGjBo5yrOyk46+BxwQZo91\n9KfWEfQYstcEXV8nMfNd0WV9oC6wDW3Me0X3PwtE0MZ8EDAV2AeUAiuBbsDHVhmn06IiVwliSKA2\nUNS3CNYW8ew90Lat39IIgju8/I7mJGZeBx1m2QDMBRYDzaLbRJfNoutHA+bRrtegPXRL9qg98G31\n/X56q5miKiRg0jcIIYFQhsfNdBNnDNoHx3TkqaiIeC5HUAl6DNlrvNbX674eTjzzSqAjcAjwNjrs\nYkZFf3ZY/ldcXMz3C7/XG+uBo6gKQezasItIJFL1WmNcxFzeblSvERNumED/S++Benvp2rE+d918\nF0V9i3yXDyJEIqnTg7v0brfTzR8ifPoptG3rrTxKpXd8pq5P0LYNgiJPrulbFZkwHLxWOjJxzyP3\n0KheI8LhMJFIhClTpgBQUFBAMtz6Zb8DyoFfo5+89UBztMd+MjrUAvBAdDkTuAv9EdWMUkpZdvpo\n/Xlr3z/EZZJQCBo0cN5iI9N07Qrz5zvzQqdN013tITNe9KuvwgUXuM8/FILFi70Ns4RC8N//wrnn\nujvu4IP1HKZBe8sQgke4OMy7rd6ttr/Xt72ITIlYHpNsPPNUnvnhwH7gR+BAoC9wD/AaMBx4MLqc\nEU3/GvAC8Ag6vHIi8Kld5obBHnrDJMr27qbwHG8/YAUVedBzg3TuU6bDVkLtwevvaKli5s2Bd9Ax\n80+A14E5aM+7L7pp4jnEPPElwEvR5VvACJKHYGhUrxEt686E1RFmPjOz1htyiDjurl4bSHxFTUZN\nKrlMVJDp5Ll/f8RzOYKKm3tbG/BaX6+b1qbyzL8EOlvs3wr0sTnm/uhPEPIO8cwFpxiO6+U3TuLH\nsppHJnzvzh/7aJQvhHM2zBK0duaZJt125vlCLt/bdMiEvkV9i2jbpIgPlzjr65EM37vzQ/55M7lq\nzPMNiZkLuYTvxjzf4m66f1X+kMsx83SQmHntJej6+m7M85GgGB4QTzIZQbpPgpAK3415vsXdYh1j\ngkGmDVa27m8mKiWJmScn357doOvruzGH/PMOxeOzpjaEWfKtLAs1x6sy47sxD3ocKhXu56CMZEMs\nxwRpbJagIe3Mk5PL9zYdgq6v700Tc5naMAelG4MVFO83W+SbvoJ3+DGrmO/GPBwO5+yraXpD+IYz\nLlemyLd25ulQr17YbxGyRr7dW6f6+uXk+R5mgdyNM2ZrUuRMEqRrXxti5oLg16xivhvzoMehkpHe\nQDmRjMiSDdIxmLl8f2fPdn/Mvn0Rz+UIKrl8b9PBqb5+OXm+G/NcJt/moAyy95uJN4zHHw+GHEJw\ncNLgwa9ZxQIRM89VzHNQHnrUbrqemnqgnDp1wnk1aqKb+xu0MEs6hlli5rWXsn1ljmLho4eNZtXj\nq6rN0zBqZGadPN+NOeS2N2PMQdmjE7xRw4Fygk6QPfNMsG+f+2NyuSwLyXHa4MFYv2vSJD77IvVo\niNLOPGeJ+C1AVsm3+ysx85rjvu9GdtiwboPlfqtYeFHfIu4dkd15GsQzTyDd9qH54LW60dG4jhvW\nbaDZs80cXcdt22ooYAAIUlnORYLcd6Me9Sz328XCs20TfDfmQYq71aQgOb1xdeqEqTyghMIrs9uh\nwAuc6hh3HQv0PifX8dprMy9bpqlfP+y3CFkjE89uen03ssPdv7nbcs7imsbCvXIAfDfmEBxvJhsF\nSdUrgePHMKsgeJ6HVwT5gXSLHz358pkg990wN3g4s/duDq6fnVi4U3w35jruFvZZCk1NCpJTz7Ci\n0d0wJDcNnVMd125ZW+WR8y3QSq8G4YF0RX33b2o6Zh7OvGwBIBKJeO6d+9WszwmRSKSqwcNbT0Lj\nxt7kW2tk0ndXAAAgAElEQVQ+gAaJmhQkx6/5B1g3kcg5Q2dDyewSVn23yvK/IDyQrjjcfU++oLxl\n5iqjh42m9efB77vhZViv1hjzIMXMs9EJKHRAM8v9uWDonBTgiS9MpLxTOcyJ7oh65QfOPDBwD2RK\n6rl/U2vQIJwhYYJHpubEHDd4Avy9ECb3onB1IRNGTgjEW6tZ36B8ozHje5gFguPNmGNidQ/cTZ+z\nnc+W7fTm1t02mv0vr4oLtWSjQ4Edbq69Ex33qD1VBpx3gBCg4PimxwfigXTFPvdvakEpy7lMrx46\nlAE1n+TYT5yWhVrjmUcikUA9ADomNpPGWzLTPjS0txGs0J5H6/f99zy89jCqQlUFwDmx5TFHHePt\nibLBZvdvanv3RjIsVHDIVDvzINkDM0HX14kxbwnMBRYDXwGjo/sPBWYDy4FZQBPTMeOBFcBS4Dxv\nRA02To1iKATs1RXGqAHZ61DgBU509Gu8moy89u4tYsINuuI97A1nFW9QDVEukQvXMFfDLPuAm4CF\nwMHAZ2gjfmV0+RAwFhgX/bUFhkaXLYD/AicBliOSBClmbsbtzXKSfscO2Ls3nJY8mcLrMIs5VEW9\n3RzRtCETHnIWqgoiRuuFXqfDdAev/BIzr70EXV8nxnx99AewE/gabaQHAr2i+59F91MfBwwCpqIr\ngVJgJdAN+NjuBGVlruXOOJmoeU87zfs8g4hhAAF6nApFfX0WyANcvXkJNSIXrmEQPXO3MfMCoBPw\nCdAMMAYr2BDdBjgaWGM6Zg3a+FsSiURYvtylFDnKihUQtLFZMj1t3ObNEfcHpcHKlVk5TUr27In4\nLULWCHoM2WuCrq+b1iwHA9OBMcBPCf+p6M+Oav8VFxdTUFBAaWkpOoLTEaOzhXHRjNeabG9DhP37\ncSwPRNiyJXX6WGeS+O0g6BuJpE5vlt9NeqfyuM3fLP9FF4FS7s7nRp5Nm5zpox9Md/Ln6raB1/l/\n8omRfzD1hQjvvw8DByZPHwo5y//HH+31jUQiTJkyBYCCggK8oB7wNnCjad9S4KjoevPoNsRi5wYz\nge4J+Skz2udTgQGUatzYXfrevZ2lM/8efTR9Gb2ia9eYPCtXJk/7j384v1dGuoEDnctivjZuyET5\nMecJSv3yl86OO/bYYJXlXOTbb4NnE8yAUlu2pE43a5YzHc4917muJHGanYRZQsDTwBLgMdP+14Dh\n0fXhwAzT/kuA+ugWxycCnzo4T04TxBiaE8xyl5Y6T5tvSMw8e+TCNXRSHoLYzvxM4HKgN7Ag+usH\nPAD0RTdNPCe6DdrovxRdvgWMIEltUv0VJjdxbugiGZQieGzaFEnruKBVHE4fuN27IxmVI0hk6tkN\nqjE36+ukfNZx+EUymzHzedgb/T42+++P/pJSXg6ff+5AAh9wa0xycSq4sWNh40bn6TNtYAcPhunT\n9fqSJdCuXWbP54agVS61maAac7dkWw9fe4D+/e/wm9+Eq7YrKvyTpaY4f9jDaRyTGR56CL77znn6\ndOQ99NCw47Snnhpb/+EH9+fKJE51P+igcEblCBKxD9D5QdD19dWYJ3qzjz1mnS4X8NswB5USFzN+\nma9h0K6nxMwFM16GWbzC97FZzDHkDdZT7PlC5sIsEZeSBIf0DGwki+fKHE7lcRsznz8fPGpxlnUk\nZp6cIH4AzRpBe4DdkMuyB4Uge+aZYsYMWL3abymCRVCNOcTKpRjzBLQS4artXH6A04mZ5xrp3Z9w\nhvPPDq+/rj/KpqJRo7CrfO+7Lz15gkDQY8heY9bXS2PuFb4a8/ffj992+jDv2JH5j6WZGGjLi2P8\nJNPyBvF6jBgRW7/11tTpg+xV5gq5cA2dhFWz3TTRV2P+yitgjqk6fZgPOQQefDATErnn9df1Mh9i\n5ukRSeuooBj2J590lz5bY9EEgXyMmWcizOIVORsz//bbzMkBzmX5/nvn6YNaSJ2STc88KMbczJtv\npk6zdWvm5ajt5MJzIjFzS8JpHRWUh924EbnYztwt2YyZ59q1AR3+y+VvIm4Jasw8UyFYs75O3sTz\n0JjHyPRwrG5w2644F3uAJhIkAxokWZxSXOy3BMIHH8ABWZjZOIhv4gEw5pGqtVw25s5vbiRNifwn\n0+3MzfnnYm9gHWKJ+CxF9shUzLwmz/aaNanTpIs5Zi6eeQqC5I059bTdh1li3HILvPii++Myhd/X\n33z+li39kyNdst3jT6iOm3tQXp7+Nw7pAWpJmCbRqaD9NiZm3Bpz52GWcNxWkAaNTHX9Mx0zz3X0\nwxv2WYrsEaSY+d69OsTixoAOGwaHHeY8vbQzd0GQXq0zEWapyXmyQWaMeXrnD1JZcIp45v7x3HNw\n1lnu7kE6PW8lzJKUSNWamwc4F2Pmmkha58ldIoC+Tl26JE9pvha5+EG5bl0Iasx8zx7429+8zTNI\nMXM9zaNxD5zh1oimOzZLtp7xABjzYCKtWaoze3bN8nczdr145t7y8cdw3XV+S5F53NyDmnjEL7+c\nOo3TDka1yDMPu+pVZRCUjgXuW7OE4/YFyTNPJcsbb6STazit8+di5ai9wrDPUmSPIMXMjecwk8bc\nrO+kSc6PS1WWa5Exh+3b9TIoTRPTKRDpfhDJJWOezfNn0phPnQoHHuh9vkH2zGs72TDmECujbobr\nzqMwS6RqPOd//MNXQarIrDGPxO2zO27XLli/3rkcwSWS1lGZDLN8/DHs3u19vrrcRLzP2AMy8SYb\npJi5QSaNuVt98zDMAvXq+S1BPG4urlF4vPbMr7kGmjd3LocX5JJnvnMnnHlmZuVxi3jmmWHVqtgH\nzlS4+QDqlJkzYdOm9I/PVsgwAMUvzIoVmcu9Vy949113x7h5KI20zpsqheP22RnQdeucy+AVmTHm\n4bTOn8oz/+EH+PDD9CTK1PcWczvzoHzTySThcJiuXasPZe01J5wATzyRPI1xvd1cd6dpzz8f7rlH\n6xvkoa4DYMzTw+lr+HvvQTisPTmnpGPMve4R5ocx8KrQbdlS8/Onqhwz4YHVlKB65jfeqA1SJpg/\nH95+OzN5m0kVcsxGzDwUcuehO22TXovCLJG0jvrnP2HpUufp3QyZm07t7vXYLLns2X35pXkrklYe\nuWjM//1vcKPv3r2ZkiSeCRP0NxivMWLIXnueVvk5lT9T7cy1TBHWrnV+TPyxmceJMX8G2ACYH9FD\ngdnAcmAW0MT033hgBbAUOM8bMa2pSRwrGenEzNPtEWZ3o9N5bawpXhW6RJlXrnR//lRvXl54wf/8\np/tjvIx/ZnJQqGySDWOVqjxkMswC7ub+BO3QPPusXg+SZz4Z6JewbxzamJ8EzIluA7QFhkaX/YAn\nUp8jHLeVzmuM16QTZvnuu9RprWLmdjfaD888WaEzmo+6J8yllzpLma0wi9F56Ve/cn9s6vBe2HFe\nufz2BbF219kw5k4NYqbCLPoDbNhxePf222Nt0YPUmuV9YFvCvoFAtN7hWeCC6PogYCqwDygFVgLd\n3AgUhIGnMlW7p+OZZ5NkD0y7ds7zSXyg5s93L0uqh8Yw5m4NyaZNNftg52WTyVw15i+8EC97Nlpr\nOD1HOo6YE556yp0c5hZ6QW/N0gwdeiG6bBZdPxowvzyuAVokzyrCUUfFttL3AL0j3QLh7FUwErfP\nqWf+7beZf/iTFboffnCeT7ycEcfHufHMjbRujeuRR7pLn0jq80Uc52UuO999lztTzi1apJfZjJmn\nKg/GtXTzxuY+VBdxrKvZmGfLM/diTg4V/SX7vxrFxcVAAVDKzp0LgY5AmMrKWCExXuMSt40HRinr\n/63Sz5kDHTo4S19REYkelzq9uaNIRUWYunXt04dC8fIbr21W6fWDHdvWHxWdye90OxYS0NuVlc7T\nv/46DBgQZu5c2LcvQv36sfQLFljnn0p+pcKceSZ8/73Wd8gQ+/Q6HGe89jrLP7H8pEpvJb/V/Xri\niUi0tYW7/Fu3jqU/7jjo0SPMhx96d3/D4XDUCDqTx+n2998b+WlWr44QiXgjL8BHH8XLC5Ho9wX7\n45cv1/+/+qpO70SeOnWcyWOuoJ1ez82bY/+//36EQw6xz3/TJvv8IpEIU6ZMAaDA6F1ZQwqI/wC6\nFDD86ebRbdCx83GmdDOB7hb5KaWU0nWWUkceGVt/4w2VEiPtO+84Txs9pSMOO8x5+ldfjeW/c2fy\ntE2bxssDSg0ZYp32/PPj5X7/fXc6OCFRlqlTnaft0SO2/4kn4tN+8EH19E7uwZgxSj36qFKXXKLU\nv/6VPO2aNTq/sjJ399epPHbyb9lSPe2pp7rXVSmlvv8+Pn2bNs50cMPLL6f3DCRj3LhYXqDUzTd7\nk6/Bd9/FywtKXX118mNeeCFez8rK1Oc591xn18Sc79y5zq7lZZfF0m3cmDztRRe5Kb/2jnO6YZbX\ngOHR9eHADNP+S4D6QCvgROBTJxmuXw9Dhzr7kGiQ7NWrshJKSpznZSbdMEuqXmpWr1NOwiynnZZa\nDqXctaW3wk1szzwWdGLoId3XRuN1tEEDPWSrk7TZHl3R6nzp6puNdumZaJKYiNcx4XTCLInhFScy\npXP904ndBylmPhX4EGgDfA9cCTwA9EU3TTwnug2wBHgpunwLGEGSmuTccwEi7Nunt994A0aMcC58\nsos0fz707+88LzOZNeaRuH1OjPlnn6WWY9o0aNw4dbpkuCl0ZiOR+PDFX7+I4zyV0no3aABLlqRO\nC0E05hFH+fz0E3z9tRcSJScThsTQ1wgH+P0BdPp07Qg6TZ+YxulQAW5i5m6+pXmFk5i5XcOyPjb7\n74/+UmIovG2bLiAHHQRlZU6O1GTqIqXbQiUdz9xOB7fe3vffu0tvxe9/D5dd5uzcP/4YW08s4DX5\noBMKwd//rtcbNIDevY1KPx4rY75xY80/cKbCK8987Fh48sn4fZlo4nf11fHblZXWzsqCBXD44fqa\nN2qkf3Yk6psNY5XMOFu1gKuocD7mU3m5fqu9/nqYMSN52nQ8c7d2YetWXRbcTGkHPvcANbe7DoXg\nrbfcHZ/swtbEoKT7Rdxs4OzThuP22emQ+MBl44v4ihXp9UpMbszDrvM59VS9vO8+eOSR5MeYH5Rm\nzbTHm0msHsz4exV2lE8mRm20IrF82RmWzp1h0CB9DS++OHmexv01PuD57ZlbVU6VlamNqFHe9u7V\n4/zoj6fJCKfUdd8+3RnMbEOcdnjSH3GhWzd3TYENAmDMY+vt20P9+s49lGTpamLcDnDRxsdckE4+\nOXlaq8JlVzjcdorxqnlTOg+mV555WZk+9sor7fNO3J/4oCTrbu30rS9ZuXLjmSfLJxvxciuSGRaj\nfKZ6y0vU1++YudX1//FH55751Klw0UXO0s6alfz/Bx6Ali3deeYGp5+ul2vWuBsv3SAAxjxStd6g\nARx8sI7/btyY+vhknT+y5Zm7OY9ubhiJ2+fUmKeq4LwyDs7jh/Zcckls/bLLIo6O2bgRJk/W19PJ\ndxM7Y57sOvznP45ESYrVR+b//c+8FalaKy+3z8evsWX++lf7/4xr+dNPyRsiJMbM/Q6zWD2Dbj78\nzpvnNGUk7k0xfgwijW6SmF6YxdAx3bLhqzE3YyjUsKF+zbv11tTHPPig/X81MW5uPPOaesR2D4Jb\nY+6VZ27XiuTQQ6vvM5q9Jsq2alVs/ZRTnJ3X/PA58aiMcyY+KMmug9MyYXetDzzQ3dDEyd4ErGTJ\nRrf4m2+2/88wJqWlcNxxzvP0O8yS7gxeRpp0J4GxKtvGfU3HM68Fxjwct2W8Jjst2F59QDRYt87d\nzXV/nnDcltOYufdyWGMXy23QQC/r1NHjTkCseWKye9W+fdhRz0a7gZJSlQM3XmFN314OO8zJ4Fjh\nqrVkTUX9MubJcGqUE2Pm6Xjmc+a4S+82Zu4G5/MdhFOmsDLmTmPmRrqcNObmwptuHM7r7v9dusCO\nHd7mmQw7PRPfDrIVZjG6aifSsKFejh8PF1wQ/1+ywnreefCzn8W2jWaoiSQa81Gj9NJtzDyZF+S0\nwrM7Z7167kY6TOaZu3n7yxbphtjcGvOffoI+dm3h8CZmbuRhd9y778LcufZ5pouVMe/VK/kxRlkw\nrn8OG/MIUP2GOP1YZfeApuupbttmvW6HU28qllckbr9dgXJ7Q430v/61u+NAG+hzztHrdm3zDc+8\nokJ/qDZjZ6AB5s2LULdurKXPffdZpzMKv3HfjArA7mE0vN7EWaqSdTaqSYXXtq1utpm6CWikam3Z\nMvtUZmPuZkz8TJLsPpqpaTvzdJ5Nt8bcMIx2Fc0777iXwUkfAqvRG3fsSN5SxggrGi3JjGfNLQEI\ns2isBpZygt0rfLoPrlkOqw8ciTh9AN32zrSLmRsfWBIxPgY//bS784Au+NOmJU9jeOb791ef2d6u\nOaM5zn7IIXr51VfWaRM91RYtYuez4qqr9HLAgPj9Xhhz8z19+WXtSS5eDD16WBvzYcOs80nWQsJ8\nf3v0qH5eP3BrzI2mdOl+AHVTCbgNsxjlxu65Mz/nXbs6lyMVVp45xEKTViSOsJjOBBgQCM88DOgm\niWbsXvch/kL97nfWadJ9jTXn7aQtsPseYWHXckCsYrFrsvSvf8XW7YyLHXv3QtOmydOYPfNE7IxA\nnTrmwYo006dbp038on/MMfbnA/tKrXt3WLjQXh637NqlW1iBbnJmFWaJ96TCHHRQ6nzN5fP66/XS\ni5ZENcGpMTdo1SoMuO8fkk7vXbfNPI28O3ZMfUx8a6RkhFOmMAxzokzJejSbbZ/be2AmMJ651auF\nlcf39df6gZ86VW936mSdn1WYwonhNdfYTpqypdO9N3GGG6vhZRPlHzNGL514M1Onak8yFErePM7A\nyWBshmfuxpgnesnnnBPrEJSIoZeRV8uWemln4BKvz/DhsXW7bvLme2C0F7bSx3xPzdf7mGOsPfNE\nR8RJpWGW//jj9TKT08iddVbqNE4rE0N249q4/cZk1xIp8X8zycq91bNulD27JpaZaudvVNJu8jd7\n5ubwsl1FZIfvnvnhh0e0ICZJDONy443VjzEmOjj9dLjtNnc1mZOPV2Y57LxIM0bBMyoXu/atsXwj\nhELxPRutei3axcydxhuNuPbkyanTjh6dOk0yz9zuDeann8zDycLIkXDssdZpE8fJMDzzDz+0Tp94\nHR5/PLZuV4GZr2nXrrr7up2Hb3DZZbH1Jk20TkZ4wcB4GLVjEWH8eIiOWmqL2TNv1UovnRrzhx/W\n16uiAiZOdHbMiBHQoUPyNE4bExiyx4Y6dkdixe3mGCusnpVU9zU9Yx5JmcLOM3dyDMSHhZJFJ6zw\n3ZjffHP1qeLq1tWG+ptvtHLmGtzcRK5NG3svzKp2dzKLuFEw2rTRcjkd8MnwXPslTrBnQ/PmsXWr\n2eztjHmq0QQTsco7EaPg/eY39mnMMXOIfxDtHpzEt6369eH1163TGpWEkX+qcSkSDV+jRrEOMXYV\nqtmDbtAATjzRuvyYy475QQuF9PeCxPJqGDel4G9/0/HRs8+2r7gg/v4ak7M4nZzi1lt1pfLDD7E3\ntkSeey5ezksvhU8+0XrXtJOPcU3SbV9uXF8nb42pWqWAteG8/PLk+dqFYZ2+af/2t9b7jWtj5UTZ\nPYvmMmZ2Htzie5ila9cwhx8e265bV38su+MObXwbN47v6PDNN3oZCmlP4/nnrVu+WN2Ua65JLY/x\nwLdurZepvnob5zFqVLteqTF5woRC8Q/6Qw9VT++VMU+WvlmzWLwW4IYb7A2QYZiNe2V+GOy+1H/1\nVXzM/Oyz9dLKmCR6a4bnbfWquXatdbjDCM3YtYQyVy7162tj/qc/WaetU8f6Q1RRUfVzGw9jeTn0\n7BkGdGW9YYO94bRqzeKGrVuTD408fHj1QccOPFCXB7cx7kQMfdu0Cad1vHGvnRhzo/zOmqW7ylth\nZcxTvWUkhsYMEu9X/NhA4aq1P/xBDwyYeB7jvhrftvr2jenwwQfW5zQ/6857o1bHd888kYULY8O5\nGsbX8J5KS3UzOgMjlGA1YXBi3iedpJcvvQT/+Ie9TMZNNgrIqFHJ5yU1zpOq8Jg9i82b9cBGBosX\nV09vLqBmj9ntjO7JhkVo2FC/ARkcc4zuMGVVARiG8I47YvsMY7Frl7VXecIJ8duNG+tjrGRK9MwN\nrD5m2vUONppX2jUpNR60sWP1A7RnD7z5ZvV0Sun/zW9PBm+/Xb28GR/Qli2LfSxt2FBXLkuXYkm6\nH+iNcjRrVszzTuVp33lnbP3EE92NTGqFYcxvusn5MaNH65ZB114bK1924Tnzs2t+c/rvf63TJ8bs\nzWMk2bVUMXQYO1afz9Al8Tk2P3v3R8eC/eMfY2NJJc5vm9h7ef58bVMeftj+rTTZ24CTN2sD3435\nokWRuH3t28diiEataNxEs9ING8aayFldpMQLNHu2Xg4dmtxDN4zoe+/F9vXubf/AGOfZtSs267vV\nDYjJE+HEE7X8lZXawJm7v1vJ//DDscpo6NDqr/lmmjaNb2b4yivOR+irV0+P1mb1dmG0pjF7NOZY\nu11sO5JQE55wgh5uNRHj+pr1vuEG6zwTR6Qzh7oGDLA35kY642Ppww/rpdXgSXbesvF9wyynWUVz\nee7a1d4bM4z5v/+tly++qMtzqo+QxnW67rrYvlThGXO380GDUocODVJViubp41KNzT5pkh6m46mn\nYs1OnYxwaR6S9pNPrNMkhtzuv19/nwHtDFi19zeMrtHI4cEH9bNo1xIKYMUKPZ75uOhcahdeqO+b\nVb4Gho6DB2snMpnT98orsfV27fS9uusu+/SJ+B5mSfaKefDB+uE0alfzhTriCL20GuvaCisvywrD\nmJ97rjZcxrZduEIpLef118da1lh9uDAe/vPP1z/Qui9apEMXiQY6sTL63/9ihTJxIH4zixbpgnP/\n/bpiPP10+49xVh7BddfBY49Z52vIbPDb38ZeJ1MNVWtwzTXw5z9XP7dhpMxe1p//rN8WPk2Yq8pc\nDhJfsR98UFfuVq/wifFXo0wktnhI5imtXauvgV0rCXMb/Cuv1BWG1YdNQ4eBA/Vy6FBdzlN9sLbq\nBp/YcSoR86Qlw4bpvgjm7wqGQevWLf64SZOs80s0WF26JB/Ay47HHkt+rQsK9Bj7Bjt3WlcaxmiD\nBpWVunPaK6/obylWo5kaOhgfs+vV08c88YS9PInfgIqLdSMJc3PhxDcuo3IuKNDNZq2iAkrpZ7V/\nf50n6HIwfrz+sO9m9jU/UEop1atX6nk8p01T6oIL9HqfPtXn39u7V28vWBB/3FdfxdL+4Q963403\nxvbNmWN9vpNO0v8vWaK3GzbU22vXWqd/7TWliopi22PHKnXVVdXTrV6t87n00ur/3XijUsOGKVVR\nEdv38MM6/e9/H592yxalmjfXc4KasZuX8JNPlDr6aKV27Kj+37HHKlVaGr9v1y6lCgqUmjWrev4/\n+1n1PJTS80z27RvbHjTIfm7TPXuU+vnPlXrppfj9X36pz7FhQ/z+p55S6rTT9HEGo0bptB06WJ9j\n8GB9HxKZMKH6/JC3367U5Zfr+VvLyvS+8nKlGjSwzlsppYYOVerOO2Pb5jkizVRWKtW/v1J33109\nj8mTdXrzPV+0SKnDD9dl1w7jPP/4h1LXXafXf/GL6nNemmVatSr+v0suUWr06Nj2+PGx63L66Uq9\n+aY+BpR65ZXqMvzzn/H5z5kTkymV3Im/Rx6pnnblSv3f3Xcr1bOnXl+0SG9366bUvn3x6YcOjeXX\nsaM+3mDnTv28fPBB/DHGnKEPPlg97YcfWsv900/VZR0zRqlrroltP/tsdR0N5s9XqkWL6vnccUf8\nc/7YY9quKaXL8bXXmuWxn7nNL1RFhVb044+rXyAzkYi+oeYLu3p1fJphw5Q65hid1uDLL5U6+eTq\nRmz7dj1hc58+8Ybjhx/0RR45Uqmzz47tNyZmvf12a/nuvFOpc86JbW/bplTbtjq9+UEtLdX5fPRR\n9Tx27FDqzDOVKi7WxlQpbXhGjrQ+5xtvaAO9eHFsHyh1333W6X/9az1prPkh2L1bH5NozJVS6q23\ndP4rVujtLVt02m7drPMvK9PX+pZbYhMtP/mkdVql9DU44gil3n03tm/hQmvjXFmp1IABekLf/fv1\nPmNC4a1brfNft05XSH/+sz6+rEypefOsDe6OHUp16aL3H3qo3rd9u67E7fjhB6WOO06p++/X19TI\nN9GhUEpP2nzkkdopMfPEE0oNH149/dNP62u5fn38/sWLrQ3E7t16QumHHortM54tUOqvf61+ji1b\ntNNyzz36+tg5As88o/cb5cDgxRf1/nff1RX5nj1a92bN9DVJrFjWro2XvaxMqWXLdEV97LHVjdvK\nlUq1aqXXKyv187Vhg9arTx+l7r03Pr1dZWrw8stKnXiiUps2xfYZFUAiL72k1AknxGyDkW9iZWCw\nbZtSLVvG7u/zz8fLk3gtrrhCqRtuiN+XaMzNbNqkr+t//mPIE0Bjbsyq/pe/zLXWIsrSpUo1aRKr\nof/v/6qnqazUD0erVkqddZZSf/+79ogOOsg6z717lbrtNm3Ur7pKe5HGxT/sMD1DvDnt7bfrG9ar\nl541vrRUG5aXXtLHmL1SpfRs3GecoVQ4rI1IZaX2dA4/XKm5c6313blTF7Djj4955Yn5mnnuOS3r\nH/8Y/5BYUV6uVL9+SvXuHfP6vvlGH2P2Ysw89ZQ2Qk8/ndxYGWzapD3yI4+Mr3Dt9J01SxfSUaO0\nl5jsYdyxQ1+LHj30Q3X//Tqt2VtPZOVKXfl07Kg9eyP/Tp2qp/3pp9hM7eZZ1ZPx/ff6/hpvckaZ\nsdL388+1R/Z//xdv2MxOgJl77tHpp0yJ6fjOO7Hjysvj05eW6redwkI9e/ydd+p0yTz8tWv19ezQ\nQaft08c63V//qivG+fNj+158UVeAibr+8INS7dvre7pvny4TlZVKLV8ek/355+Pzv/pqpQYO1JWS\nwfKvPS8AAAgVSURBVMqV+jmw4rvv9LUxV1JG3omVjpmxY5Xq3FlXiskqMKX09WvTRr+dG+kqKuzL\n8vz5+tl+8814z7y4uHraLVt0ZTFxYmzfb39rb8yV0s4V6GeGIBrzjRu1gG+9ZX2BzMydq9RNN2lP\nK/EVy8y+fdrADhumPb9kr31K6Ydg0iSlLr5YP0CffqofykSvSClt1F98UVcSzZsrVa+elr9NG3tZ\nnnxSP+zHHqsrmkMPtS8QBrNnxwzL558nl//rr3W+dl6Amf37dSVx5JFah5EjUx/zwQfaiBr5W12X\nRD76SKl//zu2nUzfDRviQ19mTz2RigpdwRj69uiRWpaKCu2VXX65PuaOO5KXny1bdDkApWbOTJ1/\nZaV+Gxw1KmZ07fTdulWX4aZNY/q+8YZ93vPm6criyCOV+s1vdLls184+fXm5Lu/HHx/L3+qty8z+\n/doonnVW9QrCzIsvajkuuUR766NH63Chla7btmnHwXjbuesupT77zL6s7dmj1IUXaudnxQpdaX3y\nSfLKdPly/fZy+eX6LcxJ5VtZqQ2m2Wu2Cn8ZPP20DrVBzJNOVpY/+EDbhUMOSf28rFql7cLVV+sK\n77e/jYWC7YhEjHwDaMznz9evh7lKebl1DC2Rykod3/7Xv6qHh7ygsjK5h5pIWZl+6H/xC+35OWHz\n5vRkc8quXfEhqWTs36/Ue+9pzzgX2bBBqddfj4WMUrFsmfYqjzoqufEx2LdPG0Sv79n27TqWW1Sk\njcqECfZpKytjIZqjj9bLk0+2T19RocNEdrFmK3bs0GG9ww7TaW+80Zkeq1frGPesWamfmxUrdEXq\ntGxu2KAdyccfT532xx9j334gtTFXygiN2htzj6Y0cI166SXFiy866zIvCPlOZaVuRePVJCTZQind\nZHfv3tjokHasW6f7XCxbpvthpEoPuqXHSy/BFVfoTnC5xoYN8MYbukPdiSemTh/SBcCyFGSqaPQD\nHgPqAv8AEid4U0VFii5doHfvSLWR9WozkYjoW5vJJ33zSVcIhr7JjHkm2pnXBf6CNuhtgUuBnycm\nOvJIPYbFwmSt9Gshom/tJp/0zSddIfj6ZsKYdwNWAqXAPuBFYFBiomee0Y3wfzSmoMkTRN/aTT7p\nm0+6QvD1zYQxbwGYhyJaE91nSWlpqavME7uH51p60dfb/IOW3o2+QZNd7m1ygqZvIpkw5q6azrh9\ndQnaDXabXvT1Nv+gpXejb9Bkl3ubnKDpm0gmPoCeDtyNjpkDjAcqif8IuhJonYFzC4Ig1GZWASek\nTOURB0RPWADUBxZi8QFUEARBCD7nA8vQHvj4FGkFQRAEQRAEQcgEzwAbgC9N+04FPgK+AF4DjBGW\nLwMWmH4VgDGUfpdoHiuACRmXOn280jcCLDX9Z5pML1C40bchMDW6fwkwznRMbby/yfSNUPvub31g\ncnT/QqCX6ZhcuL9e6RohN+6ta3oCnYi/QP+L7ge4Evh94kFAe3RYxuBTdJt1gDeJfVANGjXR1zyt\nwFygs0W6oOFG32K0cQM4EPgWMGYZrY33txh7fWvj/b0BeDq6fgRgnkQtF+6vV7oG4t5momni+0Di\nhFMnRvcD/BcYbHHcMGIPQnN0jWjMMfMccIG3YnpGTfRNmHTKt7Fy3OBG33VAI3Sv4EbAXmAHtff+\n2ulrUNvu78/RhgxgE/Aj0JXcub811dU8pbbv9zZb08YtJtYLdAjQ0iLNxcSMeQt0ZyODH0jS8SiA\nuNXX4Fn0a9od1ZMHGjt930Ybs3XoHsF/Qj8EtfX+2ulrUNvu7yJgILryaoUOrRxDbt9fN7qan2vf\n7222jPlVwAj0q8nBaI/FTHdgFzrOWBtIR9/L0KGXntHfFZkX0zPs9L0cHW5ojn4Abokuc5109K2N\n9/cZtNGeDzwKfIj+DhS4Mbdd4FZXyO17m5IC4uNQZk4CEufZfpT4j0XNAfPUrZcCaUwZmzUKqJm+\niQwHbKbTDQQFJNf34+j6E2gDZ/A0cBFwFLXr/ibTd4jFMbl+fxPLs8EHwMnk1vNbQM10TcS3e5st\nz/wI0/nuAJ5MkGEI8fHjdejX1e7oWNQVwIzMi+kZbvWtS+wLeD1gAPYFLIgk6ms8uEuBc6LrjdC9\ng5cC66ld9zeZvl9T++6vUZ4PROsJ0Bc9sN5Scvv5datrrt/bpEwF1qJfT75Hv7aMRnciWgbcn5A+\njH5lScRo2rQSmJghWb3AC30PQr++LQK+Qnvuvn9QscGNvg2A59H3cTHwG9N/tfH+2unbiNp5fwvQ\nBm0JMIv4GHIu3F8vdM2leysIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIQowK9LgZX6GH\nL72Z1G2Bj0P3XBQEQRACwk+m9SOA2eh5apMRBl7PkDyCIAhCGvyUsN0K2BxdLwDeAz6L/npE93+M\nHvFwATAG3b37T+hhXRcB12ZUYkEQBKEaicYc9DjWR6DH3GgQ3XcienIC0DPJmD3za4HfRtcbRNMV\neC2oIKTLAX4LIAg+Ux/4C3q6sAq0QYfqMfXzgA7oUR8BfgacgB63XBB8R4y5kI8cjzbcm9Cx83Xo\nkf3qAruTHDcSHW8XhMCRrSFwBSEoHIEestYYc/pn6CF5AX6FNuigQzONTce9jZ60wHCATkKPdikI\ngiBkif3YN008Af1BcyHwALE5Ow8A5kT3j4mmvw89O/uX0f9+lh3xBUEQBEEQBEEQBEEQBEEQBEEQ\nBEEQBEEQBEEQBEEQBEEQBEEQBEEQhLzj/wFlVSW9MegoZAAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 47 }, { "cell_type": "markdown", "metadata": {}, "source": [ "To find a good value for the threshold we use the find_n_peaks function, with n equal to the number of years of the hydrograph. This function finds the the n largest peaks in the hydrograph." ] }, { "cell_type": "code", "collapsed": false, "input": [ "n_peaks=round(len(sample)/365,0)\n", "peak_dist = 30\n", "POT3 = fp.find_n_peaks(sample, n_peaks=n_peaks, threshold=None, peak_dist=peak_dist)\n", "threshold = POT3.min()\n", "\n", "print 'the found threshold = ' + str(round(threshold,2)) + ' m3/s'\n", "\n", "fig = plt.figure()\n", "ax = fig.add_subplot(111)\n", "sample.plot(ax=ax,style='-b', legend=False)\n", "POT3.plot(ax=ax,style='og',markersize=6, legend=False)\n", "ax.legend(['Discharge','n peaks equal to no. of years'],fontsize=10)\n", "fig.show()" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "the found threshold = 138.98 m3/s\n" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEPCAYAAACwWiQoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXl4FEX6+D8DBBBQwQNERcMliOEIl4ICI4JRw+HqeuCx\nBK9VJKDueqD+FPerrrqsLrLqeoMXLB7LCkEMqwyLKLrIoYJAiARB5BBQJBAISf3+qOmZnkn3TPek\nZ7pnUp/nmaePqa6ut6v6rbfergMUCoVCoVAoFAqFQqFQKBQKhUKhUCgUCoVCoVAoFAqFQqFQKFJK\nJ2CF7vcLMB44BlgArAeKgea6ayYCJcBa4PxUJlahUCgU8akH/Ai0AZ4A7gqevxt4LLjfBVgJZAHZ\nwIbgdQqFQqHwCOcDi4P7a4FWwf0TgscgrfK7ddfMB85KSeoUCoWijmLXYr4SmBHcbwVsD+5vJ6zY\nTwS26K7ZApyUaAIVCoVCER87yrwhMBx42+A/EfyZEes/hUKhUNSSBjbCXgh8CewMHm9Hule2Aa2B\nHcHzPyB96honB8+FOPHEE8XWrVsTSa9CoVDUZUqBDkZ/2LHMRxF2sQC8D4wO7o8GZuvOX4m05NsC\nHYEv9BFt3boVIQRCCAYNGhTat/J78MEH0zq8kje90p9Meb2WdpW33pcXaG+moK1a5k2BIcCNunOP\nAbOA64Ey4PLg+TXB82uAw8BYYrhZGjdubDEJEr/fn9bhlbzOxu+18Hbk9VraVd7GxmvyegWhMXr0\naFGXUPJmNnVJ3rokqxDekJcYhrHr/b8LCgrcTkJKUfJmNnVJ3rokK3hfXp9L9w1WMgqFMxxzzDHs\n2bPH7WQoFI7RokULdu/eHXHO5/OBid52XZkHAoG09RUlgpI3Ofh8PpSBoMgkjMp0LGXuuptFoVAo\nFLXHdctcoXACZZkrMg1lmSsULlC/fn1yc3PJycmhR48ePPnkk6EX8csvv2TChAm24/T7/Xz55ZdO\nJ1WRobiuzAOBgNtJSClK3sykSZMmrFixgm+++YYFCxbwwQcf8NBDDwHQq1cvpkyZYjvOoBWWEIcP\nH074WkV64royVygyjeOPP54XXniBv//974Cs0IYPHw7AokWLyM3NJTc3l549e1JeXg7A448/Trdu\n3ejRowf33ntvKK63336bM888k06dOvHJJ58AUFZWxsCBA+nVqxe9evXis88+C91nwIABjBw5kpyc\nHIQQjB07ltNPP53zzz+f/Px83n33XUC2Fvx+P7179+aCCy5g27ZtKXs+iszCnR73iozF7TLVrFmz\nGueaN28uduzYIRYuXCiGDRsmhBBi+PDh4tNPPxVCCFFeXi4OHz4s5s2bJ/r37y8OHDgghBBiz549\nQggh/H6/+OMf/yiEEGLevHliyJAhQggh9u/fLyoqKoQQQqxfv1707t1bCCHEwoULRdOmTUVZWZkQ\nQoi3335bXHTRRUIIIbZt2yZatGgh3n33XXHo0CHRr18/8dNPPwkhhJg5c6a47rrrnH8oilphVKaJ\nMWjIzkRbCkVaUgtvBU5/Uz377LO5/fbbufrqq7nkkks46aST+Oijj7juuutCw8WbNw8v2nXJJZcA\n0LNnT8rKygA4dOgQ48aNY9WqVdSvX5+SkpJQ+L59+3LqqacCsGTJEi6/XM6y0apVK84991wA1q1b\nx+rVqxkyZAgAVVVVnHjiic4Kqkg5rrtZ6opPVUPJm3qESPyXKN999x3169fn+OOPjzh/99138/LL\nL3PgwAHOPvts1q1bF0yj8c0aNWoEyA+smh/8qaeeonXr1nz11VcsW7aMgwcPhsI3bdo0SnbjeM84\n4wxWrFjBihUr+Oqrr5g/f35igio8g+vKXKHINHbu3MnNN99MYWFhjf9KS0s544wzuOuuu+jTpw/r\n1q1j6NChvPrqqxw4cAAg7kjWvXv3csIJJwDw2muvUVVVZRju7LPP5t1330UIwfbt20MVa6dOndi5\ncydLly4FoLKykjVr1iQqrsIjuO5mqUujIUHJm6kcOHCA3NxcKisradCgAb/73e+44447ANkrReuZ\nMmXKFBYuXEi9evXIycnhwgsvJCsri5UrV9K7d28aNmxIfn4+Dz/8cI17aHGMHTuWSy+9lNdee40L\nLriAZs2a1QgDcOmll/LRRx/RpUsX2rRpQ8+ePTn66KPJysrinXfeYfz48fzyyy8cPnyY22+/nS5d\nuiTzESmSjBo0pMgI1KAhY8rLy2natCm7du3izDPP5NNPP6Vly5ZuJ0thgbQbNOQFn2oqUfIqUsmw\nYcPIzc1l4MCBPPDAA0qRZzCuu1kUCkXyWLhwodtJUKQI5WZRZATKzaLINNLOzaJQKBSK2uO6Mq9r\nPlUlr0KhSAauK3OFQqFQ1B7XlXld6YesoeRVOEFZWRldu3Z1OxmmZGdn11jyDODRRx91ITXxGTVq\nFN27d09odkuvoHqzKDKWogVFPP3W0xwUB2nka8T4q8aTPzQ/adcpwphN3/vnP/85YlZIL7Bt2zaW\nLVsWMcdNMjl8+DANGjivel23zOuaT1XJmxqKFhQx4ZkJFGcXs6jtIoqzi5nwzASKFhQl5bqysjJO\nP/10brrpJnJycsjLy6OioqJGuIKCAm6++Wb69OlDp06dKCqS8VZVVXHnnXfSt29funfvzgsvvADA\nvn37GDJkCL169aJbt268//77NeL87rvv6NmzJ19++SWrV6/mzDPPJDc3l+7du7Nhw4Ya4YuLi+nf\nvz+9evXi8ssvD03DO3/+fE4//XR69erF+PHjQ9P2Tpo0ib/+9a+h63Nycvj+++8B+M1vfkPv3r3J\nycnhxRdfjPmM7rnnntBI2WuvvRaAJ598kq5du9K1a1dTq7hZs2bcf//99OjRg379+rFjx47QMx88\neDDdu3dnyJAhbN68Oeb9KyoqGDNmDN26daNnz56hsnn++efzww8/kJubG5pmGODXX3+lXbt2oTlx\n9u7dS7t27aiqqqK0tJQLL7yQ3r17M3DgwNAcO3PmzOGss86iZ8+eDB06NJTWSZMmce2113LOOecw\nevRoVq9eTd++fWPmU7JoDrwDfAusAc4EjgEWAOuB4mAYjYlACbAWON8gvtCUjgsXLkz2TJKeQsmb\nHIiaLvT8gvMFk6jxyxuTFzOeRK/buHGjaNCggVi1apUQQojLL79cvPHGGzXCFRQUiAsvvFAIIURJ\nSYk4+eSTRUVFhXj++efFww8/LIQQoqKiQvTu3Vts3LhRHD58WOzdu1cIIcTOnTtFhw4dQvfLyckR\na9euFbm5ueKrr74SQghRWFgo3nzzTSGEEJWVlaFpdTV27twpBg4cKPbv3y+EEOKxxx4Tf/rTn8SB\nAwdEmzZtxIYNG0LpHz58uBBCiEmTJonJkyeH4sjJyRGbNm0SQgixe/duIYScljcnJyd0nJ2dLXbt\n2lVDfv1UwcuWLRNdu3YV+/fvF/v27RNnnHGGWLFiRY1rfD6fmDt3rhBCiLvuuiv0nIYNGyZee+01\nIYQQr7zyirj44otrXKtn8uTJ4vrrrxdCCLF27VpxyimniIMHD4qysjKRk5NjeM2YMWPE7NmzhRBC\nPP/886EpiQcPHixKSkqEEEIsXbpUDB48WAgRnr5YCCFefPFF8Yc//EEIIcSDDz4oevfuHZq6OF4+\nCZG8KXCnAPOA3wavaQrcF1TmTwB3A/cEf12AK4Lbk4D/AKcB1UYR1zWfqpI3NRwUBw3PV1TXtJad\nuA6gbdu2dOvWDZCrC2lT1kajTUvboUMH2rVrx9q1aykuLubrr7/mnXfeAaQVuGHDBk4++WQmTpzI\n4sWLqVevHlu3bg1Zezt27ODiiy/mX//6F507dwagX79+PPLII2zZsoVLLrmEDh06RNx76dKlrFmz\nhv79+wNyOt3+/fuzbt062rZtS/v27QG45pprQq2DWEyZMoXZs2cDsHnzZkpKSujbt2/c6wA++eQT\nLrnkEo444ghATve7ePFievToERFOm68G5HNdsGBBSBbt3tdccw133XVXzPstWbKE8ePHA3KysVNP\nPZX169dHzG0TzQ033MATTzzByJEjmTZtGi+99BL79u3j008/5bLLLguFO3ToUOgZXH755Wzbto1D\nhw7Rrl07QLqdRowYEZoFM14+JYIVN8vRwADgleDxYeAXYAQwPXhuOnBxcH8kMAOoBMqADYC13FUo\nHKKRr5Hh+cb1GiflOghPVwuRU9bGQ/Mv//3vfw9NS1taWsqQIUN44403+Omnn1i+fDkrVqygZcuW\nIfdN8+bNOfXUU1m8eHEorlGjRjFnzhyOOOIILrroIsMRoEOHDg3dZ/Xq1YbuEaEbrNKgQQOqq8O2\nmHb/QCDARx99xNKlS1m5ciW5ubmGrqVYcuvvI4Qw9LVnZWWF9uvVqxfxXIXNgWJ2w/fv35+ysjIC\ngQBVVVV06dKF6upqWrRoEXqG2nMEKCwsZPz48Xz11Vc8//zzoZkwQS4tqGEln+xiRZm3BXYCrwLL\ngReRlnkrYHswzPbgMcCJwBbd9VuQFrohyoec2bgl7/irxtN+RfuIc+2Xt6dwVM1paZ24zipCCN5+\n+22EEJSWlvLdd9/RuXNn8vLyePbZZ0OKav369ezfv5+9e/fSsmVL6tevz8KFC9m0aVMoroYNG/Le\ne+/x2muvMWPGDAA2btxI27ZtKSwsZOTIkXz99dcR9z/rrLNYsmQJpaWlgJyIq6SkhM6dO1NWVsZ3\n330HwIwZM0KKNTs7m+XLlwOwfPlyNm7cCMjWQ4sWLWjcuDFr164NTakbi6ysrJCMAwYMYPbs2Rw4\ncIDy8nJmz57NgAEDLD/L/v37M3PmTADefPNNBg4cGDP8gAEDePPNNwH5fL///ns6deoU9z6/+93v\nuPrqq7nuuusAOOqoo2jbtm2oFSWE4KuvvgLkM9EW+pg2bVoojuhKJF4+JYIVN0sDoCcwDvgf8Dek\nO0VPTF+O0X8FBQVkZ2dTVlbGypUr6dGjR6hJrikAdayO7Rzr0XqfTJ0xlYrqChrXa0zhuMK4vVIS\nvQ5q9uAwsjJ9Ph+nnHIKffv2Ze/evTz//PM0bNiQG264gbKyMnr27IkQgpYtWzJ79myuvvpqhg8f\nTrdu3ejduzenn356RFxNmjRh7ty5DB06lGbNmrFmzRpef/11srKyaN26Nffdd1/E/Y877jimTZvG\nqFGjQotaPPLII3Ts2JEXXniB/Px8mjRpwoABA0IKX5tuNycnJ7QeKcAFF1zAP/7xD7p06UKnTp3o\n169f3Gd000030a1bN3r16sXrr79OQUFByC1z44030r1795jPVT+d8NSpUxkzZgx/+ctfaNmyJa++\n+ioAzz//PAC///3vI+IZO3Yst9xyC926daNBgwZMnz49ZPXHWjz7qquu4v7772fUqFGhc2+++Sa3\n3HILDz/8MJWVlYwaNYpu3boxadIkLrvsMlq0aMHgwYNDla8+3QCzZs2KmU8agUAgVClkZ2ebphGs\nzc1yAvAZ0kIHOAf5gbMdcC6wDWgNLAQ6E1b0jwW384EHgc91cQq7zR2FIhbpMjfLmDFjGD58eGg5\nOK+yaNEiJk+ezJw5c9xOiuu88847zJkzh+nTp8cP7CB252axYplvAzYjP2KuB4YAq4O/0cDjwe3s\nYPj3gbeAJ5HulY7AF/bEUCgUbhPLWq0rFBYW8uGHHzJv3jy3kxIXq7nVHXgJaAiUAmOA+sAs4BTk\nh87LgZ+D4e8FrkN+LJ0AfBgVX8gyDwQCdaqHh5I3OaSLZa5QWCUZljnAKqCPwfkhJuEfDf4UCoVC\nkQLUfOaKjEBZ5opMQ81nrlAoFHUQ15W56ned2aRK3hYtWoS6f6mf+mXCr0WLFrbeAdeVuULhBLt3\n70YI4dgPBDt22L9u4cKFjqbDy7+6JKsb8hpNIRwL5TNXKAzw+WD7dlCL2Su8hM+nfOYKhUKR0biu\nzJUPObNR8mYudUlW8L68ritzhcKrKE+gIp1QPnOFwgCfD7Ztg1at4odVKFKF8pkrFAmg7A1FOuG6\nMve6H8pplLyZTV2Sty7JCt6X13VlrlB4FZ+aNFCRRiifuUJhgM8HP/4IJ5zgdkoUijDKZ65QKBQZ\njuvK3Ot+KKdR8qYPiTQe01leu9QlWcH78rquzBUKhUJRe5TPXKEwwOeDrVuhdWu3U6JQhFE+c4Ui\nAZS9oUgnXFfmXvdDOY2SN7OpS/LWJVnB+/K6rswVCq+iLHNFOqF85gqFAT4fbNkCJ53kdkoUijDK\nZ65QKBQZjuvK3Ot+KKdR8qYPqp95bOqSrOB9ea0q8zLgK2AF8EXw3DHAAmA9UAw014WfCJQAa4Hz\nnUiowl2KFhSRNyYPf4GfvDF5FC0ocjtJCoVCh1Wf+UagF6BfYfQJ4Kfg9m6gBXAP0AV4C+gDnAT8\nBzgNqNZdq3zmaUTRgiImPDOB0tzS0Ln2K9oz5dYp5A/NdzFlycPng++/hzZt3E6JQhHGKZ95dAQj\ngOnB/enAxcH9kcAMoBJp0W8A+tq4j8JjPP3W0xGKHKA0t5SpM6a6lCKFQhGNVWUukBb2MuDG4LlW\nwPbg/vbgMcCJwBbdtVuQFrohXvdDOU06yntQHDQ8X1FdEffadJS3NtQleeuSrOB9eRtYDHc28CNw\nPNJPvjbqfxH8mVHjv4KCArKzsykrK2PlypX06NEDv98PhB+aOvbGcfm2cgBoKzdslJvG9Ro7dr/P\nln1GYHWAg+Ig5dvKuWTIJUz840RX5RcitfdLt2MNr6QnE+UNBAJMmzYNgOzsbGKRSD/zB4F9SAvd\nD2wDWgMLgc5IvznAY8Ht/OA1n+viUD7zNMLQZ768PVPGOeMz96JP3ueDTZvglFNcub1CYUhtfeZN\ngCOD+02RvVO+Bt4HRgfPjwZmB/ffB64EGiJtuY6Ee8Ao0pD8oflMuXUKPZflwauDyNuU55giB+WT\nVyicwIoybwUsBlYireu5yK6IjwFDkV0TBxO2xNcAs4LbD4CxxHDBRDdhMp10lTd/aD4P/X4+bAow\n/5X5lhW5FXlr45P3Gumav4lQl2QF78trxWe+EehhcH43MMTkmkeDP4UiLo18jQzPaz55hUIRHzU3\ni8Iyc+fC8OHOT0CVbJ98Ivh8UFYGp57qyu0VCkNi+cyt9mZRKJKGprD/77mpfL68grzBjSkcV5ix\nA5IUimSg5mZJMUpeY/KH5vPUnfZ98slEzc0Sm7okK3hfXteVuUKh4XPL6adQZADKZ66wTLJ85hp3\n3gmTJ3tjUQifD777Dtq2jR9WoUgVaj5zRVqwerXbKYjEC5WKQmEV15W51/1QTqPkNaee66Wx9tSl\n/K1LsoL35c2A10eRKXjNZ64sc0U6oXzmCssk22c+fLi8hxeKhs8HJSXQoYPbKVEowiifucIRkq1k\nvWaZKxTphOvK3Ot+KKdR8prjNWWu+pnHpi7JCt6X13VlrlBoZIIyVyjcQvnMFZaZMwdGjEiekvvN\nb2D2bG8oUZ8P1q2D005zOyUKRRjlM1ekBcoyVygSx3Vl7nU/lNMoec3JBGVel/K3LskK3pfXdWWu\nUGh4TZkrFOmE8pkrLJNsn/lll8E773jDveHzwbffQufObqdEoQijfOaKtGDAALdTEIkXKhWFwiqu\nK3Ov+6GcJp3lTbYPuV07+/F7jXTOX7vUJVnB+/K6rswVCq+iLHNFOqF85grLvP8+jByZPCWX7Llf\n7ODzwTffwBlnuJ0ShSKM8pkr0gLVm0WhSBzXlbnX/VBOk87yJttn7jVlblfeogVF9LmgD/4CP3lj\n8ihaUJSchHmEdC7LieB1eRtYDFcfWAZsAYYDxwD/BE4FyoDLgZ+DYScC1wFVwHig2LnkKhSpw44y\nL1pQxIRnJlB6QikEl5orfaYUwBOLUysyH6u20B1AL+BIYATwBPBTcHs30AK4B+gCvAX0AU4C/gOc\nBlRHxad85mnIv/8NF1+cPJ/2Bx/ARRd5x2f+1VfQtau18Hlj8ijOrmm35G3KY/4r8x1OnaKuUluf\n+cnARcBLukhGANOD+9OBi4P7I4EZQCXSYt8A9E0gzYo6SDq7WQ6Kg4bnK6orHEqNQhEbK8r8KeBO\nIq3rVsD24P724DHAiUhXjMYWpIVuitf9UE6j5E0f7CjzRr5Gcmdj5PnG9Ro7lyCPkc55mwhelzee\nz3wYsANYAfhNwojgzwzD/woKCsjOzqasrIyVK1fSo0cP/H55C+2hqWNvHWtFIF3jT2Z6/Gf4Kf2k\nlNLmpSGF3n5PewrHFXpGHuefD55KTybKGwgEmDZtGgDZ2dnEIl7D9lHgWuAw0Bg4CngP6RP3A9uA\n1sBCoDPSbw7wWHA7H3gQ+DwqXuUzT0OS7TP/8EO44ALv+MxXrIAePaxfU7SgiGEFUyGrgrzBjSkc\nVag+fiocJZbP3I6XchDwR2RvlieAXcDjSAXenMgPoH0JfwDtQE3rXCnzNCRVyry62n3/uc8Hy5dD\nbq696+rXl+lXxVuRDJwcNKQV0ceAocB6YDBhS3wNMCu4/QAYS2wXTI0mTKaTzvKman7vdFaEQgTc\nTkLKSOeynAhel9dqP3OARcEfwG5giEm4R4M/hSKtSedKRVH3UHOzKCwze7ZcpzPZbpbDh6W7wk18\nPli2DHr1snedcrMokomam0WRFmgKMJ0Vodu+fkXdxXVl7nU/lNMoeePjFWWeSDqUzzxz8bq8ritz\nRfqQKiXrFWWuUKQTymeusMy//gWXXJI8ZTt/Plx4IVRUQKNGybmHVXw++N//oHdve9dlZUmfvyre\nimSgfOaKtMBrCtBr6VEoYuG6Mve6H8pplLzx8YoSTSQd1dUBx9PhVVRZ9hauK3OFIhqvKHOFIp1Q\nPnOFZZLtM583D/LzobwcmjRJzj2s4vPB0qVw5pn2rmvYECorVYWkSA7KZ65IK7yiCL2SDoXCCq4r\nc6/7oZwmE+TdsiV+GI109pkngupnnrl4XV7XlbkifdCUbJs27qYjVaRzpaKoeyifucIy770Hl14q\n95ORfUVFMGwY7N0LRx7pfPx28PlgyRLo39/edY0bw8GDqiJQJAflM1ekFUoRKhT2cV2Ze90P5TRK\n3vh4RZmrfuaxUWXZW7iuzBUKDa/NmuiVdCgUVnBdmWuLmNYVlLzx+Tx6xdg0on59v9tJSBmqLHsL\nOysNKTKMogVFPP3W0xwUB2nka8T4q8Z7YgHiyy6TH0HdRlnminTCdcvc634op/GKvEULipjwzASK\ns4tZ1HYRxdnFTHhmAkULihy9T7J95kULisgbk4e/wE/emDxH06985rHxSllOFV6XV1nmdZSn33qa\n0tzSiHOluaVMnTHV1DpPtqWqxV9dbS28ViHp5Sh9Ru57oYWhUKQS1y1zr/uhnMYr8h4UBw3PV1RX\nOHqfROS1WmnEqpCcIJHKS/nMMxevy6ss8xTiJR91I5/x6g+N6zVOcUpqYlWJJrtCUj5zRTrhumXu\ndT+UU4R81CK5PmqrjL9qPO1XtI841355ewpHFTp6n0Tyt8KiLvZihaR85pmL1+WNp8wbA58DK4E1\nwJ+D548BFgDrgWKgue6aiUAJsBY4P14CPlv2WdI+YHmJZLsE7JI/NJ8pt07hqJl58Oog8jblMWXc\nFFd9zXYt4VRVSApFOhDPzVIBnAvsD4b9BDgHGIFU5k8AdwP3BH9dgCuC25OA/wCnAYaftIoWFPHy\nJy/XiQ9YIZdA28jzTvuo7ZA/NJ+OTfP5ci3MfyV++ETcDsn0M2plZFjBVHr3r+DYIxtTOK7QsbKT\niLwNGvg5aOz9yTi87kN2Gq/La8Vnvj+4bQjUB/Yglfmg4PnpQACpzEcCM4BKoAzYAPQFlhpFnEiP\ninTFiy6BTCB/aD5szWf6Q9Cli9upUSjs4eR3NCs+83pIN8t2YCGwGmgVPCa4bRXcPxHQz3a9BWmh\nG3JQHISNNc+7aa0mi5BLQCevF1wCviTPm2nHz+i1D46JpKeqKuB4OryK133ITuO0vE6P9bBimVcD\nPYCjgQ+Rbhc9Ivgzw/C/goICNq/cLA+2AScQckHs376fQCAQatZoDzGdj5tmNWXKrVMYNuohyDpE\nnx4NefCOB8kfmu96+iBAIBA/PNgLb/c40fghwBdfQJcuzqZHiMSuT9bz8dqxhlfSk27yhjwTmoHX\nVnomHnryIZpmNcXv9xMIBJg2bRoA2dnZxMKuXfb/gAPADcg3bxvQGmmxd0a6WgAeC27nAw8iP6Lq\nEUIIw0Ef7Ze3d/1DXDLx+aBRI+s9NpJNnz6wbJk1K/Sdd+RQe0iOFf3vf8PFF9uP3+eD1auddbP4\nfPCf/8B559m7rlkzuYap11oZCu/hL/CzqO2iGucHbRxEYFrA8JpY85nHs8yPAw4DPwNHAEOBh4D3\ngdHA48Ht7GD494G3gCeR7pWOwBdmkWsK+4pbp1J+qIK8wc5+wPIq6kVPDxLJp2S7rRSZg9Pf0eL5\nzFsDHyN95p8Dc4CPkJb3UGTXxMGELfE1wKzg9gNgLLFdMDTNakqb+vNhU4D5r8zPeEUOAcvD1TOB\n6CZqLGpTySWjgkwkzsOHA46nw6vYydtMwGl5ne5aG88y/xroaXB+NzDE5JpHgz+Fos6hLHOFVTTD\n9ZrbpvJzee09E64P5w9/NKor+NPWzeK1fubJJtF+5nWFdM7bREiGvPlD8+nSPJ9P11gb6xEL14fz\nQ92zZtJVmdc1lM9ckU64rszrmt9Njq+qO6SzzzwRlM88c/G6vK4r87qIVxQPKEsyFl7KJ4UiHq4r\n87rmdwsPjPEGyVZYqcrfZFRKymcem7r27npdXteVOdQ961BZfMZkgpulrpVlRe1xqsy4rsy97oeK\nh/01KAOpSJZlvDQ3i9dQ/cxjk855mwhel9f1ronpTCasQWlHYXnF+k0VdU1ehXO4saqY68rc7/en\nbdM0sSl8/UlPV7Koa/3MEyEry+92ElJGXctbq/K6ZeS57maB9PUzpmpR5GTipWefCT5zhcKtVcVc\nV+Ze90PFIrGJcgJJSUsqSERhpnP+Llhg/5rKyoDj6fAq6Zy3iWBVXreMPNeVeTpT19ag9LL1m4wW\nxjPPeCMwK6EpAAAgAElEQVQdCu9gpcODW6uKecJnnq7o16A85oQK+nSPP1FOvXr+OjVrop389Zqb\nJRHFrHzmmUt5ZbklX/j4q8ZT+kxpjXUaCscl18hzXZlDelsz2hqU/XJhbi0nyvE6XrbMk0Flpf1r\n0rksK2JjtcODtv/g1Kl8+VX82RBVP/O0JeB2AlJKXctf5TOvPfbHbqSG7T9uNzxv5AvPH5rPw2NT\nu06DssyjSLR/aF2wWu3IqD3H7T9up9X0Vpae4549tUygB/BSWU5HvDx2I4ssw/NmvvBU6wTXlbmX\n/G61KUhWM65ePT/VDYrIG5PaAQVOYFXGiOeYLc9ZeY433ZT8tCWbhg39bichZSTj3U1s7EZqmPSH\nSYZrFtfWF+6UAeC6MgfvWDOpKEgiqwjaTaA423uWh1N4+YW0ixsj+eoyXh67oe/wcPa5FTRrmBpf\nuFVcV+bS7+Z3ORWS2hQkq5ZhVdNJcFl6KjqrMm7dtTVkkbMRaCt3vfBC2qKh/Zaa9Jn7k582DxAI\nBBy3zt3q1meFQCAQ6vDwwXNw5JHOxJsxH0C9RG0KkuVmfgPjLhJpp+hMKFpQROn3pYb/eeGFtMVx\n9kfyeaWVma6Mv2o87Zd7f+yGk269jFHmXvKZp2IQkK9BK8Pz6aDorBTgp996mgO5B+Cj4ImgVX7E\n/CM890LGJct+S61RI3+SEuM9krUm5j2XToEX8uDVQeRtymPKuCmeaLXq5fXKNxo9rrtZwDvWjN4n\nVv+ICoYMtL5attXMrb9nPIffLo1wtaRiQIEZdp69FRkPioMhBc7HgA8Q0K5FO0+8kLaotN9S80pZ\nTmcG9ZOuDKj9IsduYrUsZIxlHggEPPUCSJ/YfI7clZz+ob5DTaFEWh7tF7tveThtYYRcVdnA4PD2\n5BNOdvZGqeAn+y21Q4cCSU6Ud0hWP3Mv6QM9XpfXijJvAywEVgPfAOOD548BFgDrgWKgue6aiUAJ\nsBY435mkehurStHnAw7JCqNweOoGFDiBFRndmq8mKc3eQ/lMuVVWvMfOtVbxelURpRPp8AzT1c1S\nCdwOrASaAV8ilfiY4PYJ4G7gnuCvC3BFcHsS8B/gNMBwRhIv+cz12M0sK+H37oVDh/wJpSdZOO1m\n0buqyKrg+BaNmfKENVeVF9F6Lww6C9610ORXPvPMxevyWlHm24I/gH3At0glPQIYFDw/HTlO/R5g\nJDADWQmUARuAvsBSsxuUl9tOd9JJRs3bu7fzcXoRTQEC9OsO+UNdTpAD2Gp5KWpFOjxDL1rmdn3m\n2UAu8DnQCtAmK9gePAY4Ediiu2YLUvkbEggEWL/eZirSlJIS8NrcLEIADa3NhZFIAf7pp0Ct0meV\nDRtScpu4HDwYcDsJKcPrPmSn8bq8dnqzNAPeBSYAv0b9J4I/M2r8V1BQQHZ2NmVlZUgPTg+0wRba\nQ9OaNak+hgCHD2M5PRBg16744cODSSKP3ZT35/1FcOJNFIutoV4o3zz0DeNWjWPiHyeapj8QiB+/\nXfnsxq9//r/9LQhR++dhlp6dO63JI19Me+lP12MNp+P//HMtfm/KCwEWL4YRI2KH9/msxf/zz+by\nBgIBpk2bBkB2djZOkAV8CNymO7cWOCG43zp4DGHfucZ84Myo+IQeafMJzwBCHHmkvfDnnmstnP73\n1FOJp9Epjjr9fMEkavzyxuTVCPvSS9bzSgs3YoT1tOifjR2SUX70cYIQv/mNtetOOcVbZTkd2bjR\nezpBDwixa1f8cMXF1mQ47zzrshLDaLbiZvEBLwNrgL/pzr8PjA7ujwZm685fCTRE2nodgS+iI9VP\nc8mJedDQG9NcJooXfWix0J7/Pt/nhv8bDYxJNxmdRPnMU0c6PEMr5cGL/czPBq4BzgVWBH8XAI8B\nQ5FdEwcHj0Eq/VnB7QfAWAxqkwnPTKA4u5hFLIKbiqHjBM/MW5wI1hVdIImpsIY2q2FxdjHVx/5i\nGMapEak7dwYSus5rFYfVF66iIpDUdHiJmu4HZ/CqMtfLa6V81rP4RTKVyvyTYLgeyI+fuUjXyW5g\nCLLb4fnAz7prHgU6AJ2R7pkaRM95wWXJX73aDnaVSTotBRcxq2F7wkPvg5j1C0+2gr300vD+mjXJ\nvZddvFa5ZDJeVeZ2SbUcro8ADQ39Bg5Upe9kU9Zfdn8C1zhLxOyQ2UiF/jE0mXN0zIExiaT3mGP8\nlsN27x7e/+EH+/dKJlZlb9LEbzlOr66oY5XwB+i6gdfl9cTcLBrbNnt/sikz0slyqzE7ZLb8dVl+\nFvNfme/ovYps6Cf9M/Ta83TaZ+7lFXUU8XHSzeIUrlnmoSHfG4MnZrUnp5V3ZtVLnpslYDMlzmM0\n5J5Z7bl4QOznn5iCDSRyUdoqc6s+81gLeKQLymcem1R/AHXNMp9y6xT+8NhU1q3eBvVPgJ8KaXdR\nfIvEqyu/eE35xCJ6yD2VjeGnQs7Kdfc5etkydxovr6jjNl5V5hAul0qZ68gfmk/J6nxu/zh8Lt4D\n8nLTNBGfuZvoh9xbJTEF609y/Klhzhz5UbZLl9jhmjb1W4rPyyvqWMXrPmSn8fv9SVHmTuHqB9DF\niyOP4z2gVDZNkzHRlhPXuEmy0+vF5zF2bHj/zjvjh7f6Ars1u2Q64GXLXMOKW9WLXROTxnvvgd6n\nGu9l9mLTdM4cuU0nn3lqCSR0lVcU+3PP2QtvdS6a/KHh6XW9tqKOVbzmM09276BAIOBpy9xTvVni\nPaBUNk2tKpPNm62H9/m8o6QSIZWWuRef07x58cPs3m09Pr2rK51X1HGaRJRgql2wXvSZu9/P3IZP\n1YtNUy0j0qmfeaKk0meebs8G5Hz1Xvkmkgq85DNPhQtWL6+Vlnid+QBqRLwXWN8Lo9UpFfQ43foa\nnU6nRUPLCNsjQBsW8cLcp5nztXd65XhJgXopLVYpKHA7BXUXuy7Y2vaKU24WQwJo1ozllWy25nNR\nHrySxKapXWVu3c0SgIbl0HEC6/2laFO5e6VXTiwS72futx1/VVUi93IX6WIJUFes80AgkBTr3G45\nK1pQxDerv4kYTa5h5IJN1CUTCAQYMMAPeNMy94CbJYyXrDGrlrZ9Nwtw3NNwmfcGjCTr+Vv9MKW/\nf5s2yUlLMkn1iD9FWDHvytlleY4hJ1wyXhwB6gHL3E/z5vDzz+mtzK27WfyQNcnwH7cHjMR7/gnl\nT8Nyz44NcBr58vpdTkXq8ILPvIZi/hjwwbH7j2XKI8a9gxLtFef3+0MtRi+6WTxlS3ipaZ0MN0uI\nSm8OGEmKMj/OuhWU7m4WZZmnnhqTxg0GzoWcLjmmxkJtesVpZVS5WQwJhPbsvMBeGcBiX5kH4Kfx\n8La3euUkjarthqfjLX6RTlMKa9SvD14dR3DwIDz/vLNxJqufuZ13OxHFnGivuETnZkmVx8EDytyb\nJLU3y6F8KJEDRk6Y750BI/FkXrAggUgPZxmejmcFKcvcWZYuhZtvdjsVzmM2aVwsxawN2DpqZuID\ntt5+O34YqwOMMqhror35DjS8MuTXbm+WkE/1kOyVM+wiePHFJCXOJvFkmDs3gUh/nUT7FZE+8/bL\n21M4LvbiF+lrmftdTkXq8ILPXFPAdz0xlTUlFfTs2pjlG+J3V84fms9pzfJZttb6gC2/3x9c6B2m\nToX77rN2XXV17Io+g5Q5/BJcucyOMk9m06VeveT0ZjHKNC999E1KWg7lM+XW8AyNeYPNxwakSpnP\nmAHXXQcHDjgbr5ct80wmf2g+P2zM5/e/h0f+ARdaNDoSUaJaGd1u7D2MeU2y8UDxC5CdLfdeesnV\nhISw81LaV+aBiHNm1+3fD9u2WU+HdwkExwbMh00B5r8y31JzNplulqVLoSIJHYdkuQk4H7EDJKMl\n6wWfeTSJvLtWsStvqt0sHlDmkGXsVnUNOw9XKzxOW+Y33gitW1tPhxO43UqwY5nv2wdnn53c9NhF\nWebJobSUkHsjHtLV5Szz58POnYlfnyqXoQeKn5+SkuTFPmgQLFpk7xo7L6UW1npXJX/EOTMF+uOP\n1tPgFMlR5v6E7h/PMv/hB/j008RSlKzvLfp+5l75ppNM/H4/ffrUnMraaTp0gGefjR1Ge952nrvV\nsBdeCA89FDmfuR3qkJslMaw2w//7X/D7pSVnlUSUudMjwtxQBk4Vul27an//eJVjMiyw2uJVy/y2\n26RCSgbLlsGHHyYnbj3xXI7a+5LMd8zns2ehW+2TnkFulkBCV73+Oqxdaz38xo3xw2gkUrsn6jN3\nIg1e4+uv9UeBhOJIR2X+z3+CHXkPHUpWSiKZMkV+g3EazYfstOVpFJ/V9NspF3beMZmmAFu3Wr8m\n8trkY0WZvwJsB/Sv6DHAAmA9UAw01/03ESgB1gLnO5NMY2rjx4pFIj7zREeEmWV0Is3G2uJUoYtO\n84YN9u8fr+XlhBX8+uv2r3HS/7lli3NxuUkqlFW88pBMNwvYW/sTpEEzfbrc95Jl/ipwQdS5e5DK\n/DTk9Db3BM93Aa4Ibi8Ano1/D3/EUSLNGKdJxM3y/ffxwxr5zM0y2g3LPFah07qP2sfPqFHWQqbK\nzbJ8udz+7nf2r43v3vNbjiudW18Q7meeCmVuVSEmy80iP8D6Lbt3771X9kUHb/VmWQzsiTo3AgjW\nO0wHLg7ujwRmAJVAGbAB6GsnQUnq7WSLZNXuiVjmqSTWC3PGGdbjiX6hli2zn5Z4L42mzO0qkp07\na/fBzskuk+mqzN96KzLtqeitYfUeiRhiVtAG9llNh76Hntd7s7RCul4IblsF908E9I3HLcBJsaMK\ncMIJ4aPELUDnSLRAWGsKBiLOWbXMN25M/ssfq9D98IP1eCLTGbB8nR3LXAtrV7m2bGkvfDTx7xew\nHJe+7Hz/vb0l59xk1Sq5TaXPPF550J6lnRabfVddwLKsemWeTsP5RfAX6/8aFBQUIKc5K2PfvpVA\nD8BPdXW4kGjNuOhj7YURwvh/o/AffQRdu1oLX1UVCF4XP7x+oEhVlZ/69c3D+3yR6deabUbh5Ysd\nPpYfFa2l3+px2CUgj6urrYefMweGD/ezcCFUVgZo2DAcfsUK4/jjpV8IP2efDZs3S3kvu8w8vHTH\nac1ea/FHl5944Y3Sb5Rfzz4bCPa2sBd/+/bh8KeeCv36+fn0U+fy1+/3B5WgtfRYPd68WYtPsmlT\ngEDAmfQCfPZZZHohEPy+YH79+vXy/3//W4a3kp569aylR19BW32eP/0U/n/x4gBHH20e/86d5vEF\nAgGmTZsGQLY2urKWZBP5AXQtoNnTrYPHIH3n9+jCzQfONIhPCCGErLOEaNkyvD93roiLFvbjj62H\nDd7SEsceaz38v/8djn/fvthhW7SITA8IcdllxmEvvDAy3YsX25PBCtFpmTHDeth+/cLnn302MuyS\nJTXDW8mDCROEeOopIa68Uog334wddssWGV95ub38tZoes/Tv2lUzbPfu9mUVQojNmyPDd+pkTQY7\nvP12Yu9ALO65JxwXCHHHHc7Eq/H995HpBSGuvz72NW+9FSlndXX8+5x3nrVnoo934UJrz/Lqq8Ph\nduyIHfa3v7VTfs0N50TdLO8Do4P7o4HZuvNXAg2Rizh1BL6wEuG2bXDFFdY+JGrEanpVV0OR8YI2\ncUnUzRJvlJpRc8qKm6V37/jpEMJeX3oj7Pj2Nm0K70e7HhJtNmrN0UaN5JStVsKmenZFo/slKm8q\n+qUno0tiNE77hBNxs0S7V6ykKZHnn4jv3ks+8xnAp0AnYDMwBngMGIrsmjg4eAywBpgV3H4AjCVG\nTXLeeQABKivl8dy5MHas9cTHekjLlsGwYdbj0pNcZR6IOGdFmX/5Zfx0vPMOHHlk/HCxsFPo9Eoi\n+uWLfH4By3EKIeVu1AjWrIkfFryozAOW4vn1V/j2WydSFJtkKBJNXs0d4PYH0HfflYag1fDRYaxO\nFWDHZ27nW5pTWPGZm3UsG2Jy/tHgLy6awHv2yALSpAmUl1u5UpKsh5RoD5VELHMzGexae5s32wtv\nxJ/+BFdfbe3eP/8c3o8u4LX5oOPzwQsvyP1GjeDcc7VKPxIjZb5jR+0/cMbDKcv87rvhuecizyWj\ni9/110cem03HumIFHHecfOZNm8qfGdHypkJZxVLORj3gqqqsz/l04IBs1d5yC8yeHTtsIpa5Xb2w\ne7csC8cea+1eoXvaC+4s+n7XPh988IG962M92NoolES/iOsVnHlYf8Q5MxmiX7hUfBEvKUlsVGJs\nZe63HU/37nL7yCPw5JOxr9G/KK1aSYs3mRi9mJF55bcUTzJmbTQiunyZKZaePWHkSPkML788dpxa\n/mof8Ny2zI0qp+rq+EpUK2+HDsl5fuTH01j448paWSkHg+l1iNUBT/IjLvTta68rsIYHlHl4PycH\nGja0bqHEClcb5dbARh8ffUHq3Dl2WKPCZVY47A6Kcap7UyIvplOWeXm5vHbMGPO4o89Hvyixhltb\nbfXFKld2LPNY8aTCX25ELMWilc94rbxoed32mRs9/59/tm6Zz5gBv/2ttbDFxbH/f+wxaNPGnmWu\ncdZZcrtli7350jU8oMwDof1GjaBZM+n/3bEj/vWxBn+kyjK3cx/Z3TAQcc6qMo9XwTmlHKz7D825\n8srw/tVXByxds2MHvPqqfJ5WvpuYKfNYz+Ff/7KUlJgYfWT+3//0R4HQXqzFL9yaW+Yf/zD/T3uW\nv/4auyNCtM/cbTeL0Tto58PvJ59YDRmIaClGzkEkkV0SE3OzaDImWjZcVeZ6NIEaN5bNvDvvjH/N\n44+b/1cb5WbHMq+tRWz2IthV5k5Z5ma9SI45puY5rdtrdNpKwyvE0a2btfvqXz4rFpV2z+gXJdZz\nsFomzJ71EUfYm5o4VkvAKC2pGBZ/xx3m/2nKpKwMTj3Vepxuu1kSXcFLC5PoIjBGZVvL10Qs8wxQ\n5v6II62ZbLVgO/UBUePHH+1lrv37+COOrPrMnU+HMWa+3EbBRdDr1ZPzTkC4e2KsvMrJ8Vsa2Wg2\nUVK8cmDHKqxt6+XYY61MjuUP7cXqKuqWMo+F3aUSw4Ps7N/ro4/shbfrM7eD9fUO/HFDGClzqz5z\nLVxaKnN94U3UD+f08P9evWDvXmfjjIWZnNGtg1S5WbSh2tE0biy3EyfCxRdH/hersJ5/Phx1VPhY\n64YaTbQyLwyu92zXZx7LCrJa4ZndMyvL3kyHsSxzO62/VJGoi82uMv/1Vxhi1hcOZ3zmWhxm1y1a\nBAsXmseZKEbKfNCg2NdoZUF7/mmszANAzQyx+rHK7AVN1FLds8d43wyr1lQ4rkDEebMCZTdDtfA3\n3GDvOpAKevBguW/WN1+zzKuq5IdqPWYKGuCTTwLUrx/u6fPII8bhtMKv5ZtWAZi9jJrVG71KVazB\nRrWp8Lp0kd0243cBDYT21q0zD6VX5nbmxE8msfJRT237mSfybtpV5ppiNKtoPv7YfhqsjCEwmr1x\n797YPWU0t6LWk0x71+ziATeLxGhiKSuYNeETfXH16TD6wBGN1RfQ7uhMM5+59oElGu1j8Msv27sP\nyIL/zjuxw2iW+eHD0nesx6w7o97PfvTRcvvNN8Zhoy3Vk04K38+I666T2+HDI887ocz1efr229KS\nXL0a+vUzVuZXXWUcT6weEvr87dev5n3dwK4y17rSJfoB1E4lYNfNopUbs/dO/5736WM9HfEwsswh\n7Jo0InqGxUQWwABPWOZ+QHZJ1GPW3IfIB/X//p9xmESbsfq4rfQFtj8izG87HRCuWMy6LL35Znjf\nTLmYcegQtGgRO4zeMo/GTAnUq6efrEjy7rvGYaO/6J98svn9wLxSO/NMWLnSPD122b9f9rAC2eXM\nyM0SaUn5adIkfrz68nnLLXLrRE+i2mBVmWu0besH7I8PSWT0rt1unlrcPXrEvyayN1Is/HFDaIo5\nOk2xRjTrdZ/dPNDjGcvcqGlhZPF9+6184WfMkMe5ucbxGbkprChefY1tpStbIsN7o1e4MZpeNjr9\nEybIrRVrZsYMaUn6fLG7x2lYmYxNs8ztKPNoK3nw4PCAoGg0ubS42rSRWzMFF/18Ro8O75sNk9fn\ngdZf2EgefZ7qn/fJJxtb5tGGiJVKQ5/+du3kNpnLyJ1zTvwwVisTLe3as7H7jcmsJ1L0/3pilXuj\nd10re2ZdLJPVz1+rpO3Er7fM9e5ls4rIDNct8+OOC8iE6FKiKZfbbqt5jbbQwVlnwV132avJrHy8\n0qfDzIrUoxU8rXIx698ajjeAzxc5stFo1KKZz9yqv1Hza7/6avyw48fHDxPLMjdrwfz6q346WRg3\nDk45xThs9DwZmmX+6afG4aOfwzPPhPfNKjD9M+3TRw5fN7PwNa6+OrzfvLmUSXMvaGgvozQsAkyc\nCMFZS03RW+Zt28qtVWU+ebJ8XlVV8PTT1q4ZOxa6do0dxmpnAi3t4amO7RFdcdu5xgijdyVeviam\nzANxQ5hZ5laugUi3UCzvhBGuK/M77qi5VFz9+lJRf/edFE5fg+u7yHXqZG6FGdXuVlYR1wpGp04y\nXVYnfNIs1wuiF9gzoXXr8L7RavZmyjzebILRGMUdjVbw/vAH8zB6nzlEvohmL050a6thQ5gzxzis\nVklo8ceblyJa8TVtGh4QY1ah6i3oRo2gY0fj8qMvO/oXzeeT3wuiy6um3ISA55+X/tGBA80rLojM\nX21xFquLU9x5p6xUfvgh3GKL5rXXItM5ahR8/rmUu7aDfLRnkmj/cu35Wmk1xuuVAsaK85prYsdr\n5oa12tK+7z7j89qzMTKizN5FfRnTGw92cd3N0qePn+OOCx/Xry8/lt1/v1S+Rx4ZOdDhu+/k1ueT\nlsYbbxj3fDHKlBtvjJ8e7YVv315u43311u6j1ahmo1LD6fHj80W+6E88UTO8U8o8VvhWrcL+WoBb\nbzVXQJpi1vJK/zKYfan/5ptIn/nAgXJrpEyirTXN8jZqam7dauzu0FwzZj2h9JVLw4ZSmf/lL8Zh\n69Uz/hCVn1/z3trLeOAADBjgB2RlvX27ueI06s1ih927Y0+NPHp0zUnHjjhClge7Pu5oNHk7dfIn\ndL2W11aUuVZ+i4vlUHkjjJR5vFZGtGtMIzq/IucG8of2/u//5MSA0ffR8lX7tjV0aFiGJUuM76l/\n162PRq2J65Z5NCtXhqdz1ZSvZj2VlcludBqaK8FoweDouE87TW5nzYKXXjJPk5bJWgEpLIy9Lql2\nn3iFR29Z/PSTnNhIY/XqmuH1BVRvMdtd0T3WtAiNG8sWkMbJJ8sBU0YVgKYI778/fE5TFvv3G1uV\nHTpEHh95pLzGKE3RlrmG0cdMs9HBWvdKsy6l2ot2993yBTp4EObNqxlOCPm/vvWk8eGHNcub9gFt\n3brwx9LGjWXlsnYthiT6gV4rR8XFYcs7nqX9wAPh/Y4d7c1MaoSmzG+/3fo148fLnkE33RQuX2bu\nOf27q285/ec/xuGjffb6OZLMeqpoMtx9t7yfJkv0e6x/9x4NzgX75z+H55KKXt82evTysmVSp0ye\nbN4qjdUasNKy1nBdma9aFYg4l5MT9iFqtaKWiXqhGzcOd5EzekjRD2jBArm94orYFrqmRP/73/C5\nc881f2G0++zfH1713SgDwukJ0LGjTH91tVRw+uHvRumfPDlcGV1xRc1mvp4WLSK7Gb73nvUZ+rKy\n5GxtRq0LrTeN3qLR+9rNfNuBqJqwQwc53Wo02vPVy33rrcZxRs9Ip3d1DR9ursy1cNrH0smT5dZo\n8iQza1n7vqFPp15EfXnu08fcGtOU+T//KbczZ8ryHO8jpPacbr45fC6ee0Y/7HzkyPiuQ414laJ+\n+bh4c7NPnSqn6XjxxXC3UyszXOqnpP38c+Mw0S63Rx+V32dAGgNG/f01pat1cnj8cfkumvWEAigp\nkfOZ3xNcS+2SS2S+GcWrocl46aXSiIxl9L33Xnj/jDNkXj34oHn4aFx3s8RqYjZrJl9OrXbVP6jj\nj5dbo7mujTCysozQlPl550nFpR2buSuEkOm85ZZwzxqjDxfay3/hhfIHUvZVq6TrIlpBR1dG//tf\nuFBGT8SvZ9UqWXAefVRWjGedZf4xzsgiuPlm+NvfjOPV0qxx333h5mS8qWo1brwR/vrXmvfWlJTe\nyvrrX2Vr4Yuotar05SC6if3447JyN2rCR/tftTIR3eMhlqW0dat8Bma9JPR98MeMkRWG0YdNTYYR\nI+T2iitkOY/3wdpoGHz0wKlo9IuWXHWVHIug/66gKbS+fSOvmzrVOL5ohdWrV+wJvMz4299iP+vs\nbDnHvsa+fcaVhjbboEZ1tRyc9t578luK0Wymmgzax+ysLHnNs8+apyf6G1BBgewkoe8uHN3i0irn\n7GzZbdbIKyCEfFeHDZNxgiwHEyfKD/t2Vl9zAyGEEIMGxV/H8513hLj4Yrk/ZEjN9fcOHZLHK1ZE\nXvfNN+Gw//d/8txtt4XPffSR8f1OO03+v2aNPG7cWB5v3Woc/v33hcjPDx/ffbcQ111XM9ymTTKe\nUaNq/nfbbUJcdZUQVVXhc5Mny/B/+lNk2F27hGjdWq4JqsdsXcLPPxfixBOF2Lu35n+nnCJEWVnk\nuf37hcjOFqK4uGb8Rx1VMw4h5DqTQ4eGj0eONF/b9OBBIU4/XYhZsyLPf/21vMf27ZHnX3xRiN69\n5XUahYUybNeuxve49FKZD9FMmVJzfch77xXimmvk+q3l5fLcgQNCNGpkHLcQQlxxhRAPPBA+1q8R\nqae6Wohhw4SYNKlmHK++KsPr83zVKiGOO06WXTO0+7z0khA33yz3L7qo5pqX+jSVlkb+d+WVQowf\nHz6eODH8XM46S4h58+Q1IMR779VMw+uvR8b/0UfhNMVLd/TvySdrht2wQf43aZIQAwbI/VWr5HHf\nvkJUVkaGv+KKcHw9esjrNfbtk+/LkiWR12hrhj7+eM2wn35qnO5ff62Z1gkThLjxxvDx9Ok1ZdRY\ntmTpDF8AAAt6SURBVEyIk06qGc/990e+53/7m9RrQshyfNNN+vSYr9zmFqKqSgq6dGnNB6QnEJAZ\nqn+wmzZFhrnqKiFOPlmG1fj6ayE6d66pxH75RS7YPGRIpOL44Qf5kMeNE2LgwPB5bWHWe+81Tt8D\nDwgxeHD4eM8eIbp0keH1L2pZmYzns89qxrF3rxBnny1EQYFUpkJIxTNunPE9586VCnr16vA5EOKR\nR4zD33CDXDRW/xJUVMhropW5EEJ88IGMv6REHu/aJcP27Wscf3m5fNZ//GN4oeXnnjMOK4R8Bscf\nL8SiReFzK1caK+fqaiGGD5cL+h4+LM9pCwrv3m0c/48/ygrpr3+V15eXC/HJJ8YKd+9eIXr1kueP\nOUae++UXWYmb8cMPQpx6qhCPPiqfqRZvtEEhhFy0uWVLaZToefZZIUaPrhn+5Zfls9y2LfL86tXG\nCqKiQi4o/cQT4XPauwVC/OMfNe+xa5c0Wh56SD4fM0PglVfkea0caMycKc8vWiQr8oMHpeytWsln\nEl2xbN0amfbyciHWrZMV9Smn1FRuGzYI0bat3K+ulu/X9u1SriFDhHj44cjwZpWpxttvC9GxoxA7\nd4bPaRVANLNmCdGhQ1g3aPFGVwYae/YI0aZNOH/feCMyPdHP4tprhbj11shz0cpcz86d8rn+619a\nejyozLVV1f/+94XGUgRZu1aI5s3DNfTvf18zTHW1fDnathXinHOEeOEFaRE1aWIc56FDQtx1l1Tq\n110nrUjt4R97rFwhXh/23ntlhg0aJFeNLyuTimXWLHmN3ioVQq7G3b+/EH6/VCLV1dLSOe44IRYu\nNJZ33z5ZwNq1C1vl0fHqee01mdY//znyJTHiwAEhLrhAiHPPDVt9330nr9FbMXpefFEqoZdfjq2s\nNHbulBZ5y5aRFa6ZvMXFspAWFkorMdbLuHevfBb9+smX6tFHZVi9tR7Nhg2y8unRQ1r2Wvy5uTXD\n/vpreKV2/arqsdi8Weav1pLTyoyRvMuXS4vs97+PVGx6I0DPQw/J8NOmhWX8+OPwdQcORIYvK5Ot\nnbw8uXr8Aw/IcLEs/K1b5fPs2lWGHTLEONw//iErxmXLwudmzpQVYLSsP/wgRE6OzNPKSlkmqquF\nWL8+nPY33oiM//rrhRgxQlZKGhs2yPfAiO+/l89GX0lpcUdXOnruvluInj1lpRirAhNCPr9OnWTr\nXAtXVWVelpctk+/2vHmRlnlBQc2wu3bJyuLpp8Pn7rvPXJkLIY0rkO8MXlTmO3bIBH7wgfED0rNw\noRC33y4tregmlp7KSqlgr7pKWn6xmn1CyJdg6lQhLr9cvkBffCFfymirSAip1GfOlJVE69ZCZGXJ\n9HfqZJ6W556TL/spp8iK5phjzAuExoIFYcWyfHns9H/7rYzXzArQc/iwrCRatpQyjBsX/5olS6QS\n1eI3ei7RfPaZEP/8Z/g4lrzbt0e6vvSWejRVVbKC0eTt1y9+WqqqpFV2zTXymvvvj11+du2S5QCE\nmD8/fvzV1bI1WFgYVrpm8u7eLctwixZheefONY/7k09kZdGypRB/+IMsl2ecYR7+wAFZ3tu1C8dv\n1OrSc/iwVIrnnFOzgtAzc6ZMx5VXSmt9/HjpLjSSdc8eaThorZ0HHxTiyy/Ny9rBg0Jccok0fkpK\nZKX1+eexK9P162Xr5ZprZCvMSuVbXS0Vpt5qNnJ/abz8snS1QdiSjlWWlyyReuHoo+O/L6WlUi9c\nf72s8O67L+wKNiMQ0OL1oDJftkw2D9OVAweMfWjRVFdL//abb9Z0DzlBdXVsCzWa8nL50l90kbT8\nrPDTT4mlzSr790e6pGJx+LAQ//2vtIzTke3bhZgzJ+wyise6ddKqPOGE2MpHo7JSKkSn8+yXX6Qv\nNz9fKpUpU8zDVleHXTQnnii3nTubh6+qkm4iM1+zEXv3SrfescfKsLfdZk2OTZukj7u4OP57U1Ii\nK1KrZXP7dmlIPvNM/LA//xz+9gPxlbkQmmvUXJk7tKSBbcSsWYKZM60NmVco6jrV1bIXjVOLkKQK\nIWSX3UOHwrNDmvHjj3LMxbp1chxGvPAge3rMmgXXXisHwaUb27fD3LlyQF3HjvHD+2QBMCwFySoa\nFwB/A+oDLwHRC7yJ/HxBr15w7rmBGjPrZTKBgJI3k6lL8tYlWcEb8sZS5snoZ14f+DtSoXcBRgGn\nRwdq2VLOYbEyVi/9DETJm9nUJXnrkqzgfXmTocz7AhuAMqASmAmMjA70yiuyE/7P2hI0dQQlb2ZT\nl+StS7KC9+VNhjI/CdBPRbQleM6QsrIyW5FHDw9Pt/BKXmfj91p4O/J6Le0qb2PjNXmjSYYyt9V1\nxm7TxWsZbDe8ktfZ+L0W3o68Xku7ytvYeE3eaJLxAfQsYBLSZw4wEagm8iPoBqB9Eu6tUCgUmUwp\n0CFuKIdoELxhNtAQWInBB1CFQqFQeJ8LgXVIC3xinLAKhUKhUCgUCoUiGbwCbAe+1p3rDnwGfAW8\nD2gzLF8NrND9qgBtKv1ewThKgClJT3XiOCVvAFir+0+3mJ6nsCNvY2BG8Pwa4B7dNZmYv7HkDZB5\n+dsQeDV4fiUwSHdNOuSvU7IGSI+8tc0AIJfIB/S/4HmAMcCfoi8CcpBuGY0vkH3WAeYR/qDqNWoj\nr35ZgYVAT4NwXsOOvAVI5QZwBLAR0FYZzcT8LcBc3kzM31uBl4P7xwP6RdTSIX+dktUTeZuMromL\ngegFpzoGzwP8B7jU4LqrCL8IrZE1orbGzGvAxc4m0zFqI2/UolOuzZVjBzvy/gg0RY4KbgocAvaS\nuflrJq9GpuXv6UhFBrAT+BnoQ/rkb21l1S+p7XrepmrZuNWER4FeBrQxCHM5YWV+EnKwkcYPxBh4\n5EHsyqsxHdlMu79mcE9jJu+HSGX2I3JE8F+QL0Gm5q+ZvBqZlr+rgBHIyqst0rVyMumdv3Zk1b/X\nrudtqpT5dcBYZNOkGdJi0XMmsB/pZ8wEEpH3aqTrZUDwd23yk+kYZvJeg3Q3tEa+AH8MbtOdROTN\nxPx9Bam0lwFPAZ8ivwN5bs5tG9iVFdI7b+OSTaQfSs9pQPQ6208R+bGoNaBfunUUkMCSsSkjm9rJ\nG81owGQ5XU+QTWx5lwb3n0UqOI2Xgd8CJ5BZ+RtL3ssMrkn3/I0uzxpLgM6k1/ubTe1kjca1vE2V\nZX687n73A89FpeEyIv3HPyKbq2cifVHXArOTn0zHsCtvfcJfwLOA4ZgXMC8SLa/24q4FBgf3myJH\nB68FtpFZ+RtL3m/JvPzVyvMRSDkBhiIn1ltLer+/dmVN97yNyQxgK7J5shnZbBmPHES0Dng0Krwf\n2WSJRuvatAF4OklpdQIn5G2CbL6tAr5BWu6uf1AxwY68jYA3kPm4GviD7r9MzF8zeZuSmfmbjVRo\na4BiIn3I6ZC/TsiaTnmrUCgUCoVCoVAoFAqFQqFQKBQKhUKhUCgUCoVCoVAoFAqFQqFQKMJUIefN\n+AY5fekdxO8LfCpy5KJCoVAoPMKvuv3jgQXIdWpj4QfmJCk9CoVCoUiAX6OO2wI/Bfezgf8CXwZ/\n/YLnlyJnPFwBTEAO7/4LclrXVcBNSU2xQqFQKGoQrcxBzmN9PHLOjUbBcx2RixOAXElGb5nfBNwX\n3G8UDJftdEIVikRp4HYCFAqXaQj8HblcWBVSoUNNn/r5QFfkrI8ARwEdkPOWKxSuo5S5oi7SDqm4\ndyJ95z8iZ/arD1TEuG4c0t+uUHiOVE2Bq1B4heORU9Zqc04fhZySF+B3SIUO0jVzpO66D5GLFmgG\n0GnI2S4VCoVCkSIOY941sQPyg+ZK4DHCa3Y2AD4Knp8QDP8IcnX2r4P/HZWa5CsUCoVCoVAoFAqF\nQqFQKBQKhUKhUCgUCoVCoVAoFAqFQqFQKBQKhUKhUCgUCoWizvH/AYcaTtDIXbAVAAAAAElFTkSu\nQmCC\n", "text": [ "" ] } ], "prompt_number": 48 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Alternatively, peaks can be selected based on the maximum peak of a hydrological year." ] }, { "cell_type": "code", "collapsed": false, "input": [ "first_month = 'OCT'\n", "peak_dist = 15\n", "AM = fp.annual_max(sample,first_month=first_month,threshold=None,peak_dist=peak_dist)\n", "\n", "fig = plt.figure()\n", "ax = fig.add_subplot(111)\n", "sample.plot(ax=ax,style='-b', legend=False)\n", "AM.plot(ax=ax,style='og',markersize=6, legend=False)\n", "ax.legend(['Discharge','annual maxima peaks'],fontsize=10)\n", "fig.show()" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEPCAYAAACwWiQoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnWmYFNXVgN8GYYi4gQsgooMIKKJhWCWItMqmI2AkiqKG\nEWNiWNVowCVuiUb9NArEJa6giSgIcQFBEGgiuLGjINvoIAiMCCgCM2xzvx+3a7q6p7q7qqe7q7r7\nvM/TTy1969Y5VbdOnXvqLiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgiAIgpBWWgLLTL+fgBFA\nfWA2sA6YBRxnOuZOYD2wBuiZTmEFQRCE+NQAtgJNgMeAPwf3jwIeCa63ApYDtYB8YEPwOEEQBMEj\n9AQ+Cq6vARoE1xsGt0F75aNMx8wEzkuLdIIgCDmKU4/5amBicL0BUBpcLyVk2E8GNpuO2Qw0TlRA\nQRAEIT5OjHltoA8w2eI/FfxFI9Z/giAIQjU5wkHaS4AlwPbgdik6vLINaAR8H9z/HTqmbnBKcF8l\nJ598stqyZUsi8gqCIOQyxcAZVn848cyvIRRiAXgXGBRcHwS8bdp/NdqTbwo0Bz43Z7RlyxaUUiil\n6NatW+W6nd99992X0elF38ySP5X6ek12ubfe1xdoFs1A2/XM6wLdgZtM+x4BJgE3AiXAVcH9q4P7\nVwOHgCHECLPUqVPHpggav9+f0elF3+Tm77X0TvT1muxyb2PjNX29gjIYNGiQyiVE3+wml/TNJV2V\n8oa+xHCMXW//XVRU5LYIaUX0zW5ySd9c0hW8r6/PpfMGXzKCkH3Ur1+fXbt2uS2GkOHUq1ePnTt3\nhu3z+XwQxW677pkHAgG3RUgrom92EwgE2LVrl6MPX/KTn9XPqUPgujEXBEEQqo+EWQQhyfh8PqR8\nC9XFqhx5OswiCEJyqVmzJgUFBbRu3Zo2bdrwj3/8o9IoLFmyhJEjRzrO0+/3s2TJkmSLKiQR1415\nLsZUcwnRN/0ceeSRLFu2jC+//JLZs2czY8YMHnjgAQDatWvHmDFjHOcZ9AgT4tChQwkfK9jHdWMu\nCELqOPHEE3n++ef55z//CeiXTZ8+fQCYP38+BQUFFBQU0LZtW/bu3QvAo48+yrnnnkubNm246667\nKvOaPHkynTp1omXLlixYsACAkpISLrjgAtq1a0e7du345JNPKs/TtWtX+vXrR+vWrVFKMWTIEM46\n6yx69uxJYWEhU6ZMAXRtwe/30759e3r37s22bdvSdn2E6qMEIVtxu3wfddRRVfYdd9xx6vvvv1fz\n5s1Tl112mVJKqT59+qiPP/5YKaXU3r171aFDh9T777+vfvWrX6mysjKllFK7du1SSinl9/vV7bff\nrpRS6v3331fdu3dXSim1b98+VV5erpRSat26dap9+/ZKKaXmzZun6tatq0pKSpRSSk2ePFldeuml\nSimltm3bpurVq6emTJmiDhw4oDp37qx++OEHpZRSb7zxhho8eHDyL0oGYlWOiNFpyMlAW4IgVJNq\nRCtI9jfVLl26cOutt3LttddyxRVX0LhxY+bMmcPgwYMru64fd1xoArErrrgCgLZt21JSUgLAgQMH\nGDZsGCtWrKBmzZqsX7++Mn3Hjh057bTTAFi4cCFXXaVH/GjQoAEXXnghAGvXrmXVqlV0794dgMOH\nD3PyyScnV9EcwfUwixdijOlE9M1u4umrVOK/RPn666+pWbMmJ554Ytj+UaNG8dJLL1FWVkaXLl1Y\nu3ZtUEbrk+Xl5QH6A6sRB3/yySdp1KgRK1euZPHixezfv78yfd26dSN0t8737LPPZtmyZSxbtoyV\nK1cyc+bMxBTNcVw35oIgpI7t27dz8803M3z48Cr/FRcXc/bZZ/PnP/+ZDh06sHbtWnr06MErr7xC\nWVkZQNyOK7t376Zhw4YAvPrqqxw+fNgyXZcuXZgyZQpKKUpLSytfei1btmT79u18+umnABw8eJDV\nq1cnqm5O43qYJVNHCksU0Te78YK+ZWVlFBQUcPDgQY444gh++9vfcttttwG6VYrRMmXMmDHMmzeP\nGjVq0Lp1ay655BJq1arF8uXLad++PbVr16awsJC//e1vVc5h5DFkyBD69+/Pq6++Su/evTnqqKOq\npAHo378/c+bMoVWrVjRp0oS2bdty7LHHUqtWLd566y1GjBjBTz/9xKFDh7j11ltp1apVKi9RViKd\nhgQhyUinIWv27t1L3bp12bFjB506deLjjz/mpJNOclssz5JxnYYkpprdiL6CwWWXXUZBQQEXXHAB\n9957rxjyJON6mEUQhNxg3rx5bouQ1UiYRRCSjIRZhGSQcWEWQRAEofq4bsxzLcYo+mY3uaav4B1c\nN+aCIAhC9ZGYuSAkmVyPmRcVFdGnTx/69++flvN16dKFhQsXpuVciZCfn8/SpUupX7++o+Ocxsyl\nNYsgRDB99nTGvj6W/Wo/eb48RgwcQWGPQlfyTIUsqcbcMSkdeNmQQ/WGD3aC62GWXIsxir7eZvrs\n6Yx8eiSz8mcxv+l8ZuXPYuTTI5k+e7qt4630TTTP6sjy61//mvbt29O6dWteeOGFyv1HHXUU99xz\nD23atKFz5858//33gPamR44cSZcuXWjWrFnl8LTmIXMBhg0bxoQJEwB48MEH6dixI+eccw5/+MMf\nws5vVTPx+/3cdtttdOjQgbPOOotFixbx61//mhYtWvCXv/wlpuwbN26kRYsW7Nixg4qKCrp27cqH\nH35YqZMha7du3bj88stp1qwZo0eP5rXXXqNjx46ce+65fP311wC89957nHfeebRt25YePXpUXgMz\n48ePp1+/flx44YW0aNGCBx98sPK/f//733Tq1ImCggJuvvlmKioqAN0btkOHDrRu3Zr777+/Sp5l\nZWVccsklvPTSS+zbt4/CwkLatGnDOeecw6RJkyzvoxPsGvPjgLeAr4DVQCegPjAbWAfMCqYxuBNY\nD6wBelZbSkFIE2NfH0txQXHYvuKCYsZNHJf2PKsjy8svv8zixYtZtGgRY8eOrRxjZd++fXTu3Jnl\ny5dzwQUXhBn6bdu2sXDhQqZNm8bo0aMt8zV7mcOHD+fzzz/niy++oKysjGnTpsWUyefzkZeXx6JF\ni/jjH/9Iv379eO655/jyyy8ZP358pYxWsp922mmMGjWKP/7xjzzxxBO0bt26cqRFs0wrV67kX//6\nF1999RWvvfYaxcXFfP755/zud79j3Dh93bp27cqnn37K0qVLGTBgAI899pilvIsWLWLq1KmsXLmS\nyZMns2TJEr766ismTZrExx9/zLJly6hRowb/+c9/AHjooYdYtGgRK1asYP78+Xz55ZeVef3888/0\n7duXa6+9lhtvvJEZM2bQuHFjli9fzhdffEHv3r1jXjs72DXmY4D3gbOAc9FGejTamLcA5gS3AVoB\nA4LL3sAzsc7jhbEs0ono6232q/2W+8srym0db6VvonlWR5YxY8ZUet+bNm2qHJrWGG8F9KxDxlC2\nPp+Pyy+/HICzzjqL0tLSuOeYO3cu5513Hueeey5z5861NUBW3759AWjdujWtW7emQYMG1K5dm9NP\nP51NmzZVkX3z5s2sW7cOgBtvvJGffvqJf/3rXzz++OOW+Xfo0KEyzzPOOINevXpVns/QddOmTfTs\n2ZNzzz2Xxx9/nFWrVlnm1bNnT+rVq0edOnW44oorWLBgAXPmzGHJkiW0b9+egoIC5s6dyzfffAPA\nm2++Sbt27Wjbti2rVq2qvB5KKfr168fgwYO57rrrADj33HOZPXs2o0ePZsGCBRxzzDFxr1087Bjz\nY4GuwMvB7UPAT0BfYEJw3wTg8uB6P2AicBAoATYAHastqSCkgTxfnuX+OjXqpD3PRI8LBALMmTOH\nTz/9lOXLl1NQUEB5uX4B1KpVqzJdjRo1wqZ0q127duW6ESY54ogjKsMIoEMFPp+P8vJyhg4dypQp\nU1i5ciU33XRT5Tli6hQcRrdGjRqV62ZZImVv06ZN5bC6+/btY/Pmzfh8Pn7++eeY+Ueew6zr8OHD\nGTFiRKUXbyV3ZJxbKVW5b9CgQZVD9q5Zs4Z7772Xb775hieeeIK5c+eyYsUKCgsLK/P1+Xycf/75\nzJgxozK/5s2bs2zZMs455xzuuece/vrXv8a9dvGwY8ybAtuBV4ClwAtAXaABYLy+S4PbACcDm03H\nbwYaR8s802Kq1UX09TYjBo6g2bJmYfuaLW3G8GuqDiFrhZW+ieaZ6HG7d++u9CjXrFlTObxsIpx2\n2mmsXr2aAwcO8OOPPzJ37lyASkN1/PHHs2fPHiZPnpzwOQyUUjFlHzVqFNdffz0PPPAAN910U8Ln\n2b17d+UEGOPHj48qy+zZs9m1axdlZWW88847nH/++Vx88cW89dZbbN++HYCdO3fy7bff8vPPP1O3\nbl2OOeYYSktLwww36O8L9erVY+jQoQBs3bqVOnXqcO2113L77bezdOnShPUxsNOa5QigLTAMWAQ8\nRSikYhBzOiOr/4qKisjPz6ekpKTyDWxUUY0HQrZlO93bhT0KWbliJXc98l9+1etIjs6rQ7eu3ahb\nKzTRQrz8IjFan4ybOI7yinLq1KjD8GHD47ZKSfS43r1789xzz9GqVStatmxJ586dK/8ze5yRrU6s\n1ps0acJVV11F69atadq0KW3btgX0DEQ33XQTrVu3pmHDhnTq1ClMhngtOKxavPh8vqiyz58/nyVL\nljB27Fh8Ph9TpkxhwoQJDBo0KKoO0c53//33c+WVV1KvXj0uuugiNm7caJm+Y8eO9O/fn82bN3P9\n9ddX6v63v/2Nnj17UlFRQa1atXjmmWfo2LEjBQUFnHnmmTRp0oTzzz+/Sp5jxoxh8ODBjBo1iosv\nvpg77riDGjVqULt2bZ599llLuQOBQOULJz8/P/Y1jfmvpiHwCdpDBzgf/YHzdOBCYBvQCJgHnEnI\n0D8SXM4E7gM+M+Up7cwFT+PzwQ8/wPHHJ3JsbrczzwbGjx/PkiVLKj+aukEqxmbZBmxCf+gE6A6s\nAt4DBgX3DQLeDq6/C1wN1Ea/AJoDn9tVQBAEwW3S3VY+GdhtzTIc+A+wAt2a5SG0590D3TTxIkKe\n+GpgUnA5AxhCjBBMpsVUq4vom93kmr7ZyqBBgxg7dqzbYjjCbg/QFUAHi/3do6R/OPgThIxFIiVC\nJiFjswiCBT4fbN8OJ5yQyLESMxeqj4xnLgiCkIO4bsxzLcYo+mYOiTjXgUCAevXqVX5Ak5/8Ev3V\nq1fPUdmTURMFIcns3LnTbRHSQiAQyLjhGqqD1/WVmLkgWODzQWkpyATygpfw+SRmLgiCkNW4bswz\nOaaaCKJvdpNL+uaSruB9fV035oLgVSQSKGQSEjMXBAt8Pti2DRo0iJ9WENKFxMwFIQHE3xAyCdeN\nudfjUMlG9M1ucknfXNIVvK+v68ZcELyKL7MGzRNyHImZC4IFPh9s3QoNG7otiSCEkJi5IAhCluO6\nMfd6HCrZiL6ZQ6Jjs+QKuaQreF9f1425IAiCUH0kZi4IFvh8sGULNGrktiSCEEJi5oKQAOJvCJmE\n68bc63GoZCP6Zje5pG8u6Qre19d1Yy4IXkU8cyGTkJi5IFjg88HmzdC4sduSCEIIiZkLgiBkOa4b\nc6/HoZKN6Js5SDvz2OSSruB9fe0a8xJgJbAM+Dy4rz4wG1gHzAKOM6W/E1gPrAF6JkNQwV2mz55O\nrxt64S/y0+uGXkyfPd1tkQRBMGE3Zv4N0A4wz1T7GPBDcDkKqAeMBloBrwMdgMbAh0ALoMJ0rMTM\nM4jps6cz8umRFBcUV+5rtqwZY4aOobBHoYuSpQ6fD779Fpo0cVsSQQiRrJh5ZAZ9gQnB9QnA5cH1\nfsBE4CDao98AdHRwHsFjjH19bJghByguKGbcxHEuSSQIQiR2jblCe9iLgZuC+xoApcH10uA2wMnA\nZtOxm9EeuiVej0Mlm0zUd7/ab7m/vKI87rGZqG91yCV9c0lX8L6+R9hM1wXYCpyIjpOvifhfBX/R\nqPJfUVER+fn5lJSUsHz5ctq0aYPf7wdCF022vbG9d9teAGiqF3yjF3Vq1Ena+T5Z/AmBVQH2q/3s\n3baXK7pfwZ233+mq/kql93yZtm3gFXmyUd9AIMD48eMByM/PJxaJtDO/D9iD9tD9wDagETAPOBMd\nNwd4JLicGTzmM1MeEjPPICxj5kubMWZYcmLmXozJ+3ywcSOceqorpxcES6obMz8SODq4XhfdOuUL\n4F1gUHD/IODt4Pq7wNVAbbQv15xQCxghAynsUciYoWNou7gXvNKNXht7Jc2Qg8TkBSEZ2DHmDYCP\ngOVo73oauiniI0APdNPEiwh54quBScHlDGAIMUIwkVWYbCdT9S3sUcgDf5gJGwPMfHmmbUNuR9/q\nxOS9Rqbe30TIJV3B+/raiZl/A7Sx2L8T6B7lmIeDP0GIS54vz3K/EZMXBCE+MjaLYJtp06BPn+QP\nQJXqmHwi+HxQUgKnnebK6QXBklgxc7utWQQhZRgG+6/PjuOzpeX0uqgOw4cNz9oOSYKQCmRsljQj\n+lpT2KOQJ+9wHpNPJTI2S2xySVfwvr6uG3NBMPC5FfQThCxAYuaCbVIVMze44w54/HFvTArh88HX\nX0PTpvHTCkK6kPHMhYxg1Sq3JQjHCy8VQbCL68bc63GoZCP6RqeG66Wx+uTS/c0lXcH7+mbB4yNk\nC16LmYtnLmQSEjMXbJPqmHmfPvocXigaPh+sXw9nnOG2JIIQQmLmQlJItZH1mmcuCJmE68bc63Go\nZCP6Rsdrxlzamccml3QF7+vrujEXBINsMOaC4BYSMxds89570Ldv6ozcr38Nb7/tDSPq88HatdCi\nhduSCEIIiZkLGYF45oKQOK4bc6/HoZKN6BudbDDmuXR/c0lX8L6+rhtzQTDwmjEXhExCYuaCbVId\nM7/ySnjrLW+EN3w++OorOPNMtyURhBASMxcygq5d3ZYgHC+8VATBLq4bc6/HoZJNJuub6hjy6ac7\nz99rZPL9dUou6Qre19d1Yy4IXkU8cyGTkJi5YJt334V+/VJn5FI99osTfD748ks4+2y3JRGEEBIz\nFzICac0iCInjujH3ehwq2WSyvqmOmXvNmDvVd/rs6XTo3QF/kZ9eN/Ri+uzpqRHMI2RyWU4Er+t7\nhM10NYHFwGagD1AfeBM4DSgBrgJ+DKa9ExgMHAZGALOSJ64gpA8nxnz67OmMfHokxQ2LITjVXPHT\nxQCemJxayH7s+kK3Ae2Ao4G+wGPAD8HlKKAeMBpoBbwOdAAaAx8CLYCKiPwkZp6BvPMOXH556mLa\nM2bApZd6J2a+ciWcc4699L1u6MWs/Kp+S6+NvZj58swkSyfkKtWNmZ8CXAq8aMqkLzAhuD4BuDy4\n3g+YCBxEe+wbgI4JyCzkIJkcZtmv9lvuL68oT5I0ghAbO8b8SeAOwr3rBkBpcL00uA1wMjoUY7AZ\n7aFHxetxqGQj+mYOTox5ni9Pr3wTvr9OjTrJE8hjZPK9TQSv6xsvZn4Z8D2wDPBHSaOCv2hY/ldU\nVER+fj4lJSUsX76cNm3a4PfrUxgXTba9tW0UgUzNP5Xy+M/2U7ygmOLjiisNerNdzRg+bLhn9En+\n9cFT8mSjvoFAgPHjxwOQn59PLOJVbB8GrgcOAXWAY4Cp6Ji4H9gGNALmAWei4+YAjwSXM4H7gM8i\n8pWYeQaS6pj5Bx9A797eiZkvWwZt2tg/Zvrs6VxWNA5qldProjoMv2a4fPwUkkqsmLmTKGU34HZ0\na5bHgB3Ao2gDfhzhH0A7EvoAegZVvXMx5hlIuox5RYX78XOfD5YuhYICZ8fVrKnll+ItpIJkdhoy\niugjQA9gHXARIU98NTApuJwBDCF2CKZKFSbbyWR90zW+dyYbQqUCbouQNjK5LCeC1/W1284cYH7w\nB7AT6B4l3cPBnyBkNJn8UhFyDxmbRbDN22/reTpTHWY5dEiHK9zE54PFi6FdO2fHSZhFSCUyNouQ\nERgGMJMNoduxfiF3cd2Yez0OlWxE3/h4xZgnIofEzLMXr+vrujEXMod0GVmvGHNByCQkZi7Y5r//\nhSuuSJ2xnTkTLrkEysshLy8157CLzweLFkH79s6Oq1VLx/yleAupQGLmQkbgNQPoNXkEIRauG3Ov\nx6GSjegbH68Y0UTkqKgIJF0OryJl2Vu4bswFIRKvGHNByCQkZi7YJtUx8/ffh8JC2LsXjjwyNeew\ni88Hn34KnTo5O652bTh4UF5IQmqQmLmQUXjFEHpFDkGwg+vG3OtxqGSTDfpu3hw/jUEmx8wTQdqZ\nZy9e19d1Yy5kDoaRbdLEXTnSRSa/VITcQ2Lmgm2mToX+/fV6Km7f9Olw2WWwezccfXTy83eCzwcL\nF8KvfuXsuDp1YP9+eREIqUFi5kJGIYZQEJzjujH3ehwq2Yi+8fGKMZd25rGRsuwtXDfmgmDgtVET\nvSKHINjBdWNuTGKaK4i+8fkscsbYDKJmTb/bIqQNKcvewslMQ0KWMX32dMa+Ppb9aj95vjxGDBzh\niQmIr7xSfwR1G/HMhUzCdc/c63GoZOMVfafPns7Ip0cyK38W85vOZ1b+LEY+PZLps6cn9TypjplP\nnz2dXjf0wl/kp9cNvZIqv8TMY+OVspwuvK6veOY5ytjXx1JcUBy2r7igmHETx0X1zlPtqRr5V1TY\nS2+8kMx6FD+t171QwxCEdOK6Z+71OFSy8Yq++9V+y/3lFeVJPU8i+tp9acR6ISWDRF5eEjPPXryu\nr3jmacRLMeo8n/XsD3Vq1EmzJFWxa0RT/UKSmLmQSbjumXs9DpUsKmPUKrUxaruMGDiCZsuahe1r\ntrQZw68ZntTzJHJ/y23aYi++kCRmnr14Xd94xrwO8BmwHFgN/D24vz4wG1gHzAKOMx1zJ7AeWAP0\njCfAJ4s/SdkHLC+R6pCAUwp7FDJm6BiOeaMXvNKNXht7MWbYGFdjzU494XS9kAQhE4gXZikHLgT2\nBdMuAM4H+qKN+WPAKGB08NcKGBBcNgY+BFoAlp+0ps+ezksLXsqJD1iVIYGm4fuTHaN2QmGPQprX\nLWTJGpj5cvz0iYQdUhlnNMrIZUXjaP+rco4/ug7Dhw1PWtlJRN8jjvCz3zr6k3V4PYacbLyur52Y\n+b7gsjZQE9iFNubdgvsnAAG0Me8HTAQOAiXABqAj8KlVxom0qMhUvBgSyAYKexTClkImPACtWrkt\njSA4I5nf0ezEzGugwyylwDxgFdAguE1w2SC4fjJgHu16M9pDt2S/2g/fVN3vpreaKipDAiZ9vRAS\n8KV43EwncUavfXBMRJ7DhwNJl8OreD2GnGySrW+y+3rY8cwrgDbAscAH6LCLGRX8RcPyv6KiIjYt\n36Q3tgENqQxB7CvdRyAQqKzWGBcxk7fr1qrLmKFjuOyaB6DWATq0qc19t91HYY9C1+WDAIFA/PTg\nLL3T7UTzhwCffw6tWiVXHqUSOz5V18dr2wZekSfT9K2MTBgOXlMdmXjgHw9Qt1Zd/H4/gUCA8ePH\nA5Cfn08snPplfwHKgN+hn7xtQCO0x34mOtQC8EhwORO4D/0R1YxSSll2+mi2tJnrH+JSic8HeXn2\nW2ykmg4dYPFie17oW2/prvaQGi/6nXfg8sud5+/zwapVyQ2z+Hzw4Ydw8cXOjjvqKD2HqddqGYL3\n8Bf5md90fpX93b7pRmB8wPKYWOOZx/PMTwAOAT8CvwB6AA8A7wKDgEeDy7eD6d8FXgf+gQ6vNAc+\nj5a5YbAHDB3H3gPl9LoouR+wvIo86JlBIvcp1WErIXtI9ne0eDHzRsBcdMz8M+A9YA7a8+6Bbpp4\nESFPfDUwKbicAQwhdgiGurXq0qTmTNgYYObLM7PekEPAdnf1bCCyihqL6rzkUvGCTCTPQ4cCSZfD\nqzi5t9lAsvVNdtPaeJ75F0Bbi/07ge5Rjnk4+BOEnEM8c8EuhuN63S3j+HFv9SMTrnfnD300yhX8\nGRtm8Vo781STaDvzXCGT720ipELfwh6FtDqukI9X2+vrEQvXu/ND7nkzmWrMcw2JmQuZhOvGPNfi\nbrp/Ve6QyTHzRJCYefbidX1dN+a5iFcMD4gnGQsv3SdBiIfrxjzX4m6hjjHeINUGK133NxUvJYmZ\nxybXnl2v6+u6MYfc8w7F47MmG8IsuVaWheqTrDLjujH3ehwqHs7noAykQyzbeGlsFq8h7cxjk8n3\nNhG8rq/rTRMzmWyYg9KJwfKK95suck1fIXm4MauY68bc7/dnbNU0sSF8/SmXK1XkWjvzRKhVy++2\nCGkj1+6tXX3dcvJcD7NA5sYZ0zUpcirx0rXPhpi5ILg1q5jrxtzrcahYJDZQTiAlsqSDRAxmJt/f\n2bOdH3PwYCDpcniVTL63iWBXX7ecPNeNeSaTa3NQetn7TUUN4+mnvSGH4B3sNHhwa1YxT8TMMxXz\nHJT1G5bT4ZfxB8qpUcOfU6MmOrm/XguzJGKYJWaevew9uNdWLHzEwBEUP11cZZ6G4cNS6+S5bswh\ns70ZYw7KzgUwrZoD5XgdL3vmqeDgQefHZHJZFmJjt8GDsX7fuHEsWRl/NERpZ56xBNwWIK3k2v2V\nmHn1cd53Iz2Ubi213G8VCy/sUcjfhqR3ngbxzCNItH1oLnitTnQ0rmPp1lIaTGhg6zru2lVNAT2A\nl8pyJuLlvhu1qGW5P1osPN02wXVj7qW4W3UKkt0bl8kxc7s6hl3HfL3PznX8/e9TL1uqqV3b77YI\naSMVz25ifTfSw/1/ut9yzuLqxsKzJswC3vFm0tE+1Cu6phK32tkKmY+X+24U9ihkzNAx8Hwvuqzv\nRq+NvWJOPp/uZ911Y+6lmGp1CpJdz7CiIgDA2WfDsmV2JfMGdnUMu47fhFa98ECmGomZVw+3mvXZ\nIRAIBBs8zGTGs8mLhWeVZ+4VqlOQ7Bo648atXg0LFtiVLLPw8gOZCF9/DXv22EubCzWvVDJi4Aia\nLfV+341khvWyxph7KWaejk5ANWv6K9e9Eue1i115RwwcwamfBa9jU73w4gNpl2bNYORIe2nz8vwp\nlcVLpGqMyGNFAAAgAElEQVROzNH9dSiDV+KHMtKJWV8vPruufwAF73gz5k5ANX9RTvcL7M+W7cWb\nawcn196ujoU9Cvn2MmDFOOo3KqfDudWbddwL/PijvXReKcuZTLfOuu8GVH+SYzexWxaSVWZcN+aB\nQACfz++2GJUYnYCOPi41BUnHzP3JzzhBUvYSOqCv45n5AWa+7E/RSbzHgQMBvHR/U0kgEEiJd+7V\nF6LX9bUTZmkCzANWAV8CI4L76wOzgXXALOA40zF3AuuBNUDP5IjqbZzGzJ0c4xW8LK9XZPOqIcok\nMuEaeqW8mbFjzA8CtwJnA+cBQ4GzgNFoY94CmBPcBmgFDAguewPPxDqPl2LmZpzeLDvpd++GAwf8\nCcmTKlIRZjFzwgl+5wdlMBIzz168rq8dY74NWB5c3wN8BTQG+gITgvsnAJcH1/sBE9EvgRJgA9Ax\n1gn27nUicnpIxZu3ffvk5ymkh0RqXkJiZMI1zFTP3Ew+UAB8BjQAjMEKSoPbACcDm03HbEYbf0sC\ngQDr1jmUIkNZvx68NjZLqqeN++GHgPODEmDDhrScJi779wfcFiFtpKqPiFeNudf1dfIB9ChgCjAS\n+DniPxX8RaPKf0VFReTn51NSUoJ2/NtgfDgyLppRrUn3NgQ4dAjb8kCAHTvipw99GAvf9oK+gUD8\n9Gb5naS3K4/T/M3y/+Y3oJSz8zmRZ/t2e/roB9OZ/Jm6bZDs/D/7zMjfm/pCgI8+gr59Y6c3GnbE\ny//HH6PrGwgEGD9+PAD5+fkkg1rAB8Atpn1rgIbB9UbBbdCx89GmdDOBThH5KTPa51OeAZQ6+mhn\n6S+80F468+/JJxOXMVl06BCSZ8OG2GlffNH+vTLS9e1rXxbztXFCKsqPOU9Q6te/tnfcqad6qyxn\nIt984z2bYAaU2rEjfrpZs+zpcPHF9nUlhtNsJ8ziA14CVgNPmfa/CwwKrg8C3jbtvxqoje4y0hz4\n3MZ5MhovxtDsYJa7pMR+2lxDYubpIxOuoZ3ykO525naMeRfgOuBCYFnw1xt4BOiBbpp4UXAbtNGf\nFFzOAIYQ421StQqTmdg3dIEUSuE9tm8PJHSc114cdh+48vJASuXwEql6dr1qzM362imfNWx+kUxn\nzHwB0Y1+9yj7Hw7+YlJWBkuX2pDABZwak0wc1nbUKPj+e/vpU21g+/eHKVP0+urVejAyr+C1l0s2\n41Vj7pScGjXx+efhT3/yV24fPuyeLNXF/sPuT+CY1PDYY/Dtt/bTJyJv/fp+22l/+cvQ+nffOT9X\nKrGr+5FH+m3n6dUZdewS+gCdG3hdX1e780d6s089BX/6kzuyVBe3DbNXme7APpmvodeuZ7Jj5l6e\nUUeITzLDLMnC9VETzTHkUusp9lwhdWGWgENJvENiBjaQxnOlDrvy2I2ZZ8MEHhIzj03ODbRlxs4F\nSnSOzlTjNeOTiXjZM082Xp5Rx228aswhVC7FmEeglfBXbse7QF6umiYSM880EjOw/hTnnx7ee09/\nlG3VKna6unX9tvLLhgk8vB5DTjZ+vz8lxjxZuBpm+eij8O14FyidVdNUDLSVjGPcJNXyevF6DBkS\nWr/jjvjp7T7A6ZgIJVPxsmduYCesmu6mia4a86lTwRxTjfcwe7Fq+t57epkLMfPECCR0lFcM+7PP\nOktvdywa8+TAXptRxy5ei5mnunVQIBDwtGeeUTHzdFZN7RqTTZvsp/f5vGOkEiGdnrkXr9P778dP\ns3On/fyMiVAgs2fUSTaJGMF0h2C9GDP3QGsWv+2UXqyaGjciE9uZOyWdMfNMuzagx6vP5G8iTvFS\nzDwdIVizvnZq4jn1ATSSeA+weY7OBqeW0+as1M0t6bRdcSb2AI3ESwbUS7LYpajIbQlyF6ch2Oq2\nipMwiyUBDG/GzgUyqqaX9oKXU1g1dWrM7YdZAmSq95Z4O3O/4/wzsTewDrEEyNT765RAiubETKSc\nOQnBJhqSCQQCdO3qB7zpmXsgzBLCS96YXU/beZglxO23wxtvOD8uVbh9/c3nb9LEPTkSJd09/oQQ\nTkKwyQjJSA9QS/wcF5wK2m1jYsapMbcfZvGHbXlp0Mh41z/VMfNMRz+8fpelqEqqWnl4KWbeo1sh\nv79wDG0+j986KNFWcWZ9JcwSBy9VrVMRZqnOedJBaox5Yuf3Ulmwixc9cy93tEsmr74Ko24pZOrU\nQq64InbroOq0ijPKqIRZLAlUrjl5gL3SgcW5MQ8kdJ7MJQDo69SuXeyU5muRiR+Ua9YEr/UjSGUr\nj1S1M0/kmdDTPBr3IDaJtopLdGyWdD3jnvLMvYS0ZqnK7NnVy9/J2PWp8sxTObaPFz1zL3a0SyV2\n7oFxvweOGMfusnJ6XeS8VdzkyXD33bHTmDsYxfK+s6hporPxDgy80uXXaWuWyJiq1XFuDSYWT4dp\n0xLJ1Z/Q+VPxckx1yEF7hf5q55NMUtnRzksxc+M5tPtCLexRSIujClm8xn6HLb/fX1kDGDcuvjE3\nqKiILVcWGXP46Se9dGLMU1l1qVEjNa1ZrG5a5HFuxjjdDvmk2pjHCjkk49p60TMfMXAExU8Xh+nd\nbGkzhg/LrjFgnBpz8zFOMMqok+G60/VceaD4BcjP12svvpjcnBP9ip9IgbBvzANh+yKPy4ZxrsMJ\nJHRUKsIsqQ456HITSEpeycIYA6bj8uSPAfP3x/+eklYy1TF+qTTmTr8R2I04ZJVnXqtW8vOsjofr\n5OIahSdZnrmbMc5M8sz37IFevWDhQvv5p3psHy965qDL+9F5hXTrlrwxYKbPns4/3/4nW3psqdyX\nqhpkcTGcdhocYcNa2fkA6pSZM+N/vI9Fur6neaD4+Vm/Pvm5VsfDdfJQGmntN1Xyh+2LNKBujnOd\nGmPuT+j88Tzz776Djz92Jkmqx/YxtzP3yjedVDH29bFhhhxSV4M84wx45pnYaYzr7eS62017ySXw\nwAPh45k7QVqzxCHew14dDzcRY56sHmFuxjiTVeh27Kj++eO9HBPxwAyPceh949i4JbFWDLHwqmd+\nyy3wwgvJzTPdNcht22L/n46Yuc8HmzfbT2+3TXpWtjN3wmuvwZo10f+vjoebyNs90Zh5JEaMs/5b\nmTvO9RdfmLcCCeWRCmMO+vpe3nEmbAww8+WZSb2ub74JTvQ9cCBpp47JmDGwb19y88zz5cE3Vfcn\nowZp9SzZld9JuXDynGuZAmzZEi9ltGNTjx1j/jJQCpgf0frAbGAdMAs4zvTfncB6YA3QMzliWrN9\ne/T/qlOlTiRmnmiPMKsbXdijkLYNtcH54JXkGpxYJKvQReq5YYPz88ereSXDC37tNefHJDP+6cTL\n8xojBo7g5A0nh+1L5XDU8cpDKsMs4GzuT9AOzYQJej1dnrmdMMsrwDjgVdO+0Whj/hgwKrg9GmgF\nDAguGwMfAi2AGOr4w7a2b4cTT7QnfKwLax4ul1rl/KpjHe6yWaVOJMzy7bfx01rFzKPdaDdirrEK\nndF81Dl+rrnGXspUh1kMjM5Lv/0tXH+9s2MPH45XPvy288rkuHphj0Ke53kuKxpHk2bltDojdcNR\ng32DmKowi25f7rfdyuquu0L9MrzUmuUjID9iX1+gW3B9ArpuORroB0wEDgIlwAagI/CpXYECAbjy\nSrupY2OeyeW59+Ccc+wdl6q3u13P3Gm+ySLWA3P22fbziXygFi92Lku8h8Yw5k5rE9u3V5171gmH\nDyev9VWmGvPXX4drrwWl9PN1zXXw6KOpPafdGlEijpgdjG8OduUwlxGvt2ZpgA69EFw2CK6fDJgr\nj5vRHnoMAjRsGNpK3ANMHokWCHtVwUDYPrue+TffpP7hj1XovvvOfj7hcgZsH+fEMzfSOm2PftJJ\nztJHEv98Adt5mcvOt986m3LOTVas0Euj3XWyY8JW+cUrD8a1dFJjcx6qC9jW1WzMveSZx0MFf7H+\nr0JRURHa4S9hz57lQBvAT0VFqJAY3YUjt40HRinr/63Sz5kD55xjL/3hw4HgcfHTmzuKHD7sp2bN\n6Ol9vnD5jWqbVXr9YIe29UdFa3n+/vjfmfrhVOo2rEueLw//2X46t+8cVb9Q54dweSoq7Kd/7z3o\n08fPvHlw8GCA2rVD6Zcts84/3vVUyk+XLrBpk9b3yiujp9ffS4xqr738I8tPvPRW8lvdr2eeCQRb\nWzjLv1mzUPrTToPOnf18/HF8+Z1sayNoTx6725s2GflpNm4MEAgkR16ATz4JlxcCwe8L0Y9ft07/\n/847Or0deWrUsCeP+QVt93r+8EPo/48+CnDssdHz3749en6BQIDx48cDkG/0rqwm+YR/AF0DGP50\no+A2hGLnBjOBThb5KaWU0u8spU46KbQ+bZqKi5F27lz7aYOntMXxx9tP/847ofz37Imdtl69cHlA\nqSuvtE57ySXhcn/0kbVM02ZNU836NVPcT+WvWb9matqs+BcyUpaJE+2n7dw5tP+ZZ8LTLlxYNb2d\nezBypFJPPqnU1Vcr9Z//xE67ebPOb+9eZ/fXrjzR5N+xo2raX/7Sua5KKbVpU3j6li3t6eCEyZMT\newZiMXp0KC9Q6rbbkpOvwbffhssLSt14Y+xjXn89XM+Kivjnufhie9fEnO+8efau5bXXhtJ9/33s\ntL/5jZPyG91xTjTM8i4wKLg+CHjbtP9qoDbQFGgOfG4nw23bYMAAex8SDWJVvSoqYHqCvYsTDbMY\ng/BEw6o6ZSfM0r599DyT2f3fSWxv48bQemToIdFqo1EdzcuD/dbNmKukTfe451bnS1TfdLRLt9Ok\nr7qTVyQ7JpxImCUyvGJHpkSufyKxey/FzCcCHwMtgU3ADcAjQA9008SLgtsAq4FJweUMYAgx3iQX\nXwwQ4OBBvT1tGgwZYl/4WBdp8WK47DL7eZlJrTEPhO2zY8yXLImeZzI7bzgpdGYjEfnwhV+/gO08\nldJ65+XB6tXx04IXjXnAVj4//wxffZUMiWIT754aw17Myp/F/KbzmZU/i5FPj4xp0A19jXBAOoxV\nrHNMmaIdQbvpI9PEe25D2I+ZO/mWlizsmK1r0B82awNN0E0VdwLd0c0OewI/mtI/DJwBnAl8EPPk\nwbPv2qULyJFHOhM+VRcp0RYqiXjm0XSwK0Myu/8/+KD9j1k/mu545DHV+aDj88Hzz8Pjj8M998Cc\nOdbprIz5998nfl67JMszHzUKuncP35eKziU33hi+HWnkEqnZReqbDmMVyzhbjX/lRKayMti6FS6/\nvHpymKlOjX3nzsR6UbvaA9Tc7trngxkznB0f68JWx6Ak+kXcbOCip/WH7YumQ2TtINqDnszxRtav\nT6xXYmxj7neczy9/qZcPPQT/+EfsY8wPSoMG2uNNJVYPZvi98tvKpzxNc0NElq9I+ROp2Rn31/iA\n57ZnblWTrqiIb0SN8nbggB7nR388jYU/rq4HD+rOYGYbYrfDk/6ICx07OmsKbOABYx5ab90aate2\n76HESlcdY25ndDYDc0E688zYaa0KV7TCYfeFYnT/b/m/5HT/T+TBTJZnvnevPvaGG6LnHbk/8kGJ\n1d167157csQqV04881j5uDWOS6T8idTsIvV1O2Zudf1//NF+f4CJE+E3v7GXdtas2P8/8gg0aeLM\nMzc47zy93LzZ2XjpBh4w5oHK9bw8OOooeOste1XmWJ0/0uWZOzmPbm4YCNtn15jH6+06pF9yxhux\nHz+MztVXh9avvTZg65jvv4dXXtHX0853k2jGPJaR/O9/bYkSkz17qu5btMi8FahcKyuLnk8qhmq1\nw3PPhW9b1eyO+G8zBnSPXrOLjJk7CWkk+rHVqTF3MhbNggV2UwbCaorhYxBpdJPExMIsho6Jlg1X\njbkZQ6E6deCqq+COO+IfE6vXWXU8HyeeeXUb/Ed7EJwY82TIYRCtFUn9+lX3Gc1eI2UrNoVgzz3X\n3nnND58dj8o4Z+SDEus62C0T0a71L36hY6t2iVUTsJIlHQMy3XZb+LZRs+P5XvzizW7wfC8OfTWG\nwdfadwjseuaJfGy1cw4nPaut0sQbkTEaVmXbuK+JeOZZYMz9YVtGNdluwa7uB8RItm51dnOdn8cf\ntmU3Zp58OayJFsvNC9bGa9TQ405AqHlirHvVurXfVs/GaAMlxSsHTrzC6oY2jj/ezuBY/so1Ky8+\nlizpGl0vEj3sxUyaHA7AlplwwN7kLaFOdvbOU51mtE5j5k6YP99uSn/cFFbG3G7M3EiXkcbcXHgT\njcMlu/t/u3awe3dy84xFND0jawfxHvRkxWCNrtqR1AmGUO+8s+pX/1iFtWdPOOaY0LbRDDWSSGM+\nPFjLdxozj+UF2X3hRTtnrVrORjqM5Zk7qf2li0RDbHaNud2PrcmImRt5RDtu/nyYNy96noliZcy7\ndYt9jFEWjOufwcY8AFS9IXY/VkV7QBP1VHftsl6Phl1vKpRXIGx/tALl9IYa6X/3O2fHgTbQF12k\n16O1zTc888OH9YdqM9EMNMCCBQFq1gy19HnoIet0RuE37pvxAoj2MBpeb+QsVbE6G1XnhdeqlR5c\natOmeCkDlWtr10ZPZTbmTsbETyWx7qOZRNuZV6cZrVNjbhjGaC+auXPjntKCQNwUVqM37t4du6WM\nEVY0WpLlWV+muHggzKKxGljKDtGq8Ik+uGY5rD5wRGL3AYxV5bYiWszc+MASifEx+KWXnJ0HdMF/\n663YaQzP/NAhHTs2E605oznOfuyxevnll9ZpIz3Vxo1D57Ni8GC97NMnfH8yjLn5nk6erJs7rloF\nnTtbG/OBA63zidVCwnx/O3euel43cGrMjaZ0dj3z6jSjdRpmMcpNtOfO/Jx36BD39Lax8swhFJq0\nInKExUQmwABPeOZ+QDdJNBOtug/hF+ovf7FOk2g11py3nbbAznuE+R3LAaEXS7QmS//5T2g9mnGJ\nxoEDUK9e7DRmzzySaEagRg3zYEWaKVOs00Z+0T/llOjng+gvtU6dYPny6PI4Zd8+3cIKdJMzqzBL\nuCflt9X5zVw+//hHvbQT5qhu1/tY2DXmBk2b+gH7/UOMj60Xf62b0fYssd+M1mkzT6PctGkT/5jw\n1kix8MdNYRjmSJli9Wg22z6n98CMZzxzq6qFlcf31Vf6gZ84UW8XFFjnZxWmsGN4zW9sO03ZEune\nGznDjdXwspHyjxypl3aqtBMnak/S54vdPM7AzmBshmfuxJhHeskXXRTqEBSJoZeRV5MmehnNwEVe\nn0GDQuvRusmb74HRXthKH/M9NV/vU06x9swjHRE7Lw2z/KefrpfxOmxVpzXI+efHl8luzNyQ3bg2\nTr4xFfYoZOo43Yz27Wesm9E6jZlbPetG2Ys21lOq2vkbL2kn+Zs9c3N4OdqLKBque+YnnBDQgpgk\nMYzLLbdUPcaY6OC88+DPf3b2JrPz8cosRzQv0oxR8IyXS7T2raF8A/h84T0brXotRouZ2/0WYMS1\nX3klftoRI+KnieWZR6vB/PyzeThZGDYMTj3VOm3kOBmGZ/7xx9bpI6/D00+H1qO9wMzXtEMHOOGE\n6B6+wbXXhtaPO07rZIQXDIyHUTsWAe68E4KjlkbF7Jk3baqX8Yx5dVqDDBkSf3IWu40JDNlDQx07\nI/LF7eQYK6yelXj3NTFjHoibIppnbucYCA8LxYpOWOG6Mb/ttqpzedasqQ31119r5cweg7mJXMuW\n0b0wq7f7BzFHigmdG3Te27fbH/DJ8Fx7945/DoBGjULrVuMwRDPm8UYTjMTOGA9GwfvTn6KnMcfM\nIfxBjPbgRNa2ateG996zTmu8JIz8jz8+uixQ1fDVrRvqEBPthWr2oPPyoHlz6/JjLjvmB83n098L\nIsurYdyUgn/9S8dHL7gg+osLwu+vMTlLvCacTrrev/pquJzXXAOffab1ru5YKsY1SbTnp3F97dQa\n47VKAWvDed11sfONFoa1W9O++27r/ca1sXKioj2L5jJmdh6c4nqYpUMHPyecENquWVN/LLvnHm18\njz46vKPD11/rpc+nPY1//9u65YvVTbnppvjyGA98s+B3mnhfvY3zGG/UaL1SQ/L48fnCH/THHqua\nPlnGPFb6Bg1C8VqAoUOjGyDDMBv3yvwwRPtS/+WX4THzCy7QSytjEumtGZ63VVVzyxbrcIcRmonW\nEsr8cqldWxvz//s/67Q1alh/iCosrHpu42EsK4OuXf2AflmXlkY3nFatWeLhpDXIoEFVZ1X6xS90\neXA6BlIkhr4tW/oTOt6413aMuVF+Z83SXeWtsDLm8WoZkaExg8j7FT42kL9y7a9/1QMDRp7HuK/G\nt60ePUI6LFxofU7zs26/N2pVXPfMI1m+XLesOProkPE1vKeSEt2MzsAIJVhNGByZd4sWejlpErz4\nYnSZjJtsFJDhw61HZYs8T7zCY/YsfvgB2rYNba9aVTW9uYCaPWanM7rHGhahTh1dAzI45RTdYcrq\nBWAYwnvuCe0zjMW+fdZe5RlnhG8ffbQ+xkqmSM/cwOpjZrTewUbzymhNSo0HbdQo/QDt3w/vv181\nnVL6f3PtyeCDD6qWN+MD2tq1oY+lderol8uaNViSyAf6EQNH0GxpeGsQJjVjyFWxW4Pce29ovXlz\n+81+o2EY81tvtX/MiBG6ZdDvfx8qX9HCc+Zn11xz+vBD6/SRMXvzGEnRWqoYOowapc9n6BL5HJuf\nvYcf1su//z00llTk/LaRvZcXL9Y25fHHo9dKY9UGnIye6LoxX7EiELavdetQDNF4Kxo30ax0nTqh\nJnJWFynyAs2erZcDBsT20A0j+r//hfZdeGF0D8s4z759oVnfrW5ASJ4AzZtr+SsqtIErLo6VXhcE\n42U0YEDVar6ZevXCmxlOnWp/hL5atfRobVa1C6M1jdmjMcfao8W2AxFvwjPOgGXLqqYzrq9Z76FD\nrfOMHJHOHOrq0ye6MTfSGR9LH39cL60GT4rmLRvfN8xymlU0l+cOHaJ7Y4Yxf/NNvXzjDV2eY32E\nLOxRyBM36673vKK73rNhDJ3bxW4NYu523q9f/NChQbyXonn6uHhjs48bp4fpeOGFULNTOyNcvv12\naP2zz6zTRIbcHn5Yf58B7QxYtfc3jK7RyOHRR/WzGK0lFMD69Xo889HBudSuuELfN6t8DQwd+/fX\nTmQsp2/q1ND62Wfre3XffdHTR+J6mCVWFfOoo/TDabxdzRfqxBP1Uk9wER8rL8sKw5hffLE2XMZ2\ntHCFUlrOP/4x1LLG6sOF8fBfcon+gdZ9xQoduog00JEvo0WLQoUyciB+MytW6ILz8MP6xXjeedE/\nxll5BDffDE89ZZ2vIbPB3XeHqpPxhqo1uOkmeOKJquc2jLnZy3riCV1b+DxiripzOYisYj/6qH65\nW1XhI+OvRpmIbPEQy1PaskVfg2itJMxt8G+4Qb8wrD5sGjr07auXAwboch7vg3WeT3e9Z2Ogsut9\nZMepSI4+OrQ+cKDui2D+rmAYtI4dw48bF+W7aqTBateu6gBednjqqdjXOj9fj7FvsGeP9UvDGG3Q\noKJCd06bOlV/S7EazdTQwfiYXauWPuaZZ6LLE/kNqKhIN5IwNxeOrHEZL+f8fN1s1ioqoJR+Vi+7\nTOcJuhzceaf+sO9k9jU3UEop1a1b/Hk833pLqcsv1+vdu1edf+/AAb29bFn4cV9+GUr717/qfbfc\nEto3Z471+Vq00P+vXq2369TR21u2WKd/912lCgtD26NGKTV4cNV0GzfqfK65pup/t9yi1MCBSh0+\nHNr3+OM6/YMPhqfdsUOpRo30nKBmos1L+NlnSp18slK7d1f979RTlSopCd+3b59S+flKzZpVNf9j\njqmah1J6nskePULb/fpFn9t0/36lzjpLqUmTwvd/8YU+R2lp+P4XXlCqfXt9nMHw4TrtOedYn6N/\nf30fIhkzpur8kHfdpdR11+n5W/fu1fvKypTKy7POWymlBgxQ6t57Q9vmOSLNVFQoddllSt1/f9U8\nXnlFpzff8xUrlDrhBF12o2Gc58UXlbr5Zr1+6aVV57w0y1RcHP7f1VcrNWJEaPvOO0PX5bzzlHr/\nfX0MKDV1alUZXnstPP85c0IyxZM78vePf1RNu2GD/u/++5Xq2lWvr1ihtzt2VOrgwfD0AwaE8mvT\nRh9vsGePfl4WLgw/xpgz9NFHq6b9+GNruX/+uaqsI0cqddNNoe0JE6rqaLB4sVKNG1fN5557wp/z\np57Sdk0pXY5//3uzPNFnbnMLdfiwVvTTT6teIDOBgL6h5gu7cWN4moEDlTrlFJ3W4IsvlDrzzKpG\n7Kef9ITN3buHG47vvtMXedgwpS64ILTfmJj1rrus5bv3XqUuuii0vWuXUq1a6fTmB7WkROfzySdV\n89i9W6kuXZQqKtLGVClteIYNsz7ntGnaQK9aFdoHSj30kHX63/1OTxprfgjKy/UxkcZcKaVmzND5\nr1+vt3fs0Gk7drTOf+9efa1vvz000fKzz1qnVUpfgxNPVGr+/NC+5cutjXNFhVJ9+ugJfQ8d0vuM\nCYV37rTOf+tW/UJ64gl9/N69Si1YYG1wd+9Wql07vb9+fb3vp5/0Szwa332n1GmnKfXww/qaGvlG\nOhRK6UmbTzpJOyVmnnlGqUGDqqZ/6SV9LbdtC9+/apW1gSgv1xNKP/ZYaJ/xbIFSzz1X9Rw7dmin\n5YEH9PWJ5gi8/LLeb5QDgzfe0Pvnz9cv8v37te4NGuhrEvli2bIlXPa9e5Vau1a/qE89tapx27BB\nqaZN9XpFhX6+Sku1Xt27K/W3v4Wnj/YyNZg8WanmzZXavj20z3gBRDJpklJnnBGyDUa+kS8Dg127\nlGrSJHR///3vcHkir8X11ys1dGj4vkhjbmb7dn1d//tfQx4PGnNjVvV//nOetRZB1qxR6rjjQm/o\nP/yhapqKCv1wNG2q1PnnK/X889ojOvJI6zwPHFDqz3/WRn3wYO1FGhf/+OP1DPHmtHfdpW9Yt256\n1gN8ls8AAAkvSURBVPiSEm1YJk3Sx5i9UqX0bNy/+pVSfr82IhUV2tM54QSl5s2z1nfPHl3ATj89\n5JVH5mvm1Ve1rH//e/hDYkVZmVK9eyt14YUhr+/rr/UxZi/GzAsvaCP00kuxjZXB9u3aIz/ppPAX\nbjR9Z83ShXT4cO0lxnoYd+/W16JzZ/1QPfywTmv21iPZsEG/fNq00Z69kX9BQdW0P/8cmqndPKt6\nLDZt0vfXqMkZZcZK36VLtUf2hz+EGzazE2DmgQd0+vHjQzrOnRs6rqwsPH1Jia7t9OqlZ4+/916d\nLpaHv2WLvp7nnKPTdu9une655/SLcfHi0L433tAvwEhdv/tOqdat9T09eFCXiYoKpdatC8n+73+H\n53/jjUr17atfSgYbNujnwIpvv9XXxvySMvKOfOmYGTVKqbZt9Usx1gtMKX39WrbUtXMj3eHD0cvy\n4sX62X7//XDPvKioatodO/TLYuzY0L67745uzJXSzhXoZwYvGvPvv9cCzphhfYHMzJun1K23ak8r\nsopl5uBBbWAHDtSeX6xqn1L6IRg3TqmrrtIP0Oef64cy0itSShv1N97QL4lGjZSqVUvL37JldFme\nfVY/7Keeql809etHLxAGs2eHDMvSpbHl/+ornW80L8DMoUP6JXHSSVqHYcPiH7NwoTaiRv5W1yWS\nTz5R6s03Q9ux9C0tDQ99mT31SA4f1i8YQ9/OnePLcviw9squu04fc889scvPjh26HIBSM2fGz7+i\nQtcGhw8PGd1o+u7cqctwvXohfadNi573ggX6ZXHSSUr96U+6XJ59dvT0ZWW6vJ9+eih/q1qXmUOH\ntFE8//yqLwgzb7yh5bj6au2tjxihw4VWuu7apR0Ho7Zz331KLVkSvazt36/UFVdo52f9ev3S+uyz\n2C/Tdet07eW663QtzM7Lt6JCG0yz12wV/jJ46SUdaoOQJx2rLC9cqO3CscfGf16Ki7VduPFG/cK7\n++5QKDgagYCRrweN+eLFunqYqZSVWcfQIqmo0PHt//ynangoGVRUxPZQI9m7Vz/0l16qPT87/PBD\nYrLZZd++8JBULA4dUup//9OecSZSWqrUe++FQkbxWLtWe5UNG8Y2PgYHD2qDmOx79tNPOpZbWKiN\nypgx0dNWVIRCNCefrJdnnhk9/eHDOkwULdZsxe7dOqx3/PE67S232NNj40Yd4541K/5zs369fpHa\nLZulpdqRfPrp+Gl//DH07QfiG3OljNBodGOepCkNHKMmTVK88Ya9LvOCkOtUVOhWNMmahCRdKKWb\n7B44EBodMhpbt+o+F2vX6n4Y8dKDbukxaRJcf73uBJdplJbCtGm6Q13z5vHT+3QBsCwFqSoavYGn\ngJrAi0DkBG+qsFDRrh1ceGGgysh62UwgIPpmM7mkby7pCt7QN5YxT0U785rAP9EGvRVwDXBWZKKT\nTtJjWCyP1Uo/CxF9s5tc0jeXdAXv65sKY94R2ACUAAeBN4B+kYleflk3wv/RmIImRxB9s5tc0jeX\ndAXv65sKY94YMA9FtDm4z5KSkhJHmUd2D8+09KJvcvP3Wnon+npNdrm3sfGavpGkwpg7ajrjtOri\ntRvsNL3om9z8vZbeib5ek13ubWy8pm8kqfgAeh5wPzpmDnAnUEH4R9ANQMTwb4IgCEIcioEz4qZK\nEkcET5gP1AaWY/EBVBAEQfA+lwBr0R74nXHSCoIgCIIgCIKQCl4GSoEvTPt+CXwCrATeBYwRlq8F\nlpl+hwFjKP12wTzWA2NSLnXiJEvfALDG9J9pMj1P4UTfOsDE4P7VwGjTMdl4f2PpGyD77m9t4JXg\n/uVAN9MxmXB/k6VrgMy4t47pChQQfoEWBfcD3AA8GHkQ0BodljH4HN1mHeB9Qh9UvUZ19DVPKzAP\naGuRzms40bcIbdwAfgF8AxizjGbj/S0iur7ZeH+HAi8F108EzJOoZcL9TZaunri3qWia+BEQOeFU\n8+B+gA+B/hbHDST0IDRCvxGNOWZeBS5PrphJozr6Rkw65dpYOU5wou9WoC66V3Bd4ACwm+y9v9H0\nNci2+3sW2pABbAd+BDqQOfe3urq2Nx3n+r1N17Rxqwj1Ar0SaGKR5ipCxrwxurORwXfE6HjkQZzq\nazABXU27p2pyTxNN3w/Qxmwrukfw/6Efgmy9v9H0Nci2+7sC6It+eTVFh1ZOIbPvrxNdzc+16/c2\nXcZ8MDAEXTU5Cu2xmOkE7EPHGbOBRPS9Fh166Rr8XZ96MZNGNH2vQ4cbGqEfgNuDy0wnEX2z8f6+\njDbai4EngY/R34E8N+a2A5zqCpl9b+OST3gcykwLIHKe7ScJ/1jUCDBP3XoNkMCUsWkjn+rpG8kg\nIMp0up4gn9j6fhpcfwZt4AxeAn4DNCS77m8sfa+0OCbT729keTZYCJxJZj2/+VRP10hcu7fp8sxP\nNJ3vHuDZCBmuJDx+vBVdXe2EjkVdD7ydejGThlN9axL6Al4L6EP0AuZFIvU1Htw1wEXB9bro3sFr\ngG1k1/2Npe9XZN/9NcrzL9B6AvRAD6y3hsx+fp3qmun3NiYTgS3o6skmdLVlBLoT0Vrg4Yj0fnSV\nJRKjadMGYGyKZE0GydD3SHT1bQXwJdpzd/2DShSc6JsH/Bt9H1cBfzL9l433N5q+dcnO+5uPNmir\ngVmEx5Az4f4mQ9dMureCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIAiCIIQ4jB4340v08KW3\nEb8t8GnonouCIAiCR/jZtH4iMBs9T20s/MB7KZJHEARBSICfI7abAj8E1/OB/wFLgr/Owf2fokc8\nXAaMRHfv/j/0sK4rgN+nVGJBEAShCpHGHPQ41ieix9zIC+5rjp6cAPRMMmbP/PfA3cH1vGC6/GQL\nKgiJcoTbAgiCy9QG/omeLuww2qBD1Zh6T+Ac9KiPAMcAZ6DHLRcE1xFjLuQip6MN93Z07HwremS/\nmkB5jOOGoePtguA50jUEriB4hRPRQ9YaY04fgx6SF+C3aIMOOjRztOm4D9CTFhgOUAv0aJeCIAhC\nmjhE9KaJZ6A/aC4HHiE0Z+cRwJzg/pHB9A+hZ2f/IvjfMekRXxAEQRAEQRAEQRAEQRAEQRAEQRAE\nQRAEQRAEQRAEQRAEQRAEQRAEQcg5/h+xqrmaBSJpLwAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 49 }, { "cell_type": "markdown", "metadata": {}, "source": [ "The following figure compares the peaks selected with all three method" ] }, { "cell_type": "code", "collapsed": false, "input": [ "fig = plt.figure()\n", "ax = fig.add_subplot(111)\n", "sample.plot(ax=ax,style='-b', legend=False)\n", "AM.plot(ax=ax,style='*y',markersize=10, legend=False)\n", "POT2.plot(ax=ax,style='or',markersize=10, legend=False)\n", "POT3.plot(ax=ax,style='og',markersize=5, legend=False)\n", "ax.legend(['Discharge','AM','peaks over threshold method','peaks equal to no. of years'],fontsize=10)\n", "fig.show()\n" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXMAAAEPCAYAAACwWiQoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXl4FEX6+D9DDjDACkFFZIGgAgsJkHBnERyQQFRAjgT0\niysBhB+yguKBoCuwyrq4uuLBuoK64AVIgnKsgiDLcCgoICiCXMoYQUEgyJEAuer3R0/P9Ezm6JnM\n0ZPU53nmmemequqq7uq333r7rbdAIpFIJBKJRCKRSCQSiUQikUgkEolEIpFIJBKJRCKRSCQSiSSs\ntAJ2aT5ngUlAIrAOOAisBepp8kwDDgH7gb7hrKxEIpFIfFMD+AVoAvwDmGLb/xgw2/a7DbAbiAOS\ngMO2fBKJRCIxCH2Bzbbf+4GGtt/X2rZB0cof0+RZA3QLS+0kEomkmuKvxnwnsNj2uyFwwvb7BA7B\nfh1wVJPnKNA40ApKJBKJxDf+CPN4YACQ6+Y/Yft4wtt/EolEIqkksX6kvRXYCZy0bZ9AMa8cBxoB\nv9r2H0Oxqav83rbPznXXXSd+/vnnQOorkUgk1ZnvgRvd/eGPZn4XDhMLwEpgpO33SGC5Zv+dKJp8\nc6AF8KW2oJ9//hkhBEIIbr75ZvtvPZ8ZM2ZEdXrZ3uiqfyjba7S6y2tr/PYCN3gS0Ho189pAH2Cs\nZt9sYCkwBrACw2z799n27wNKgQl4MbPUqlVLZxUUzGZzVKeX7Q1u+UZL7097jVZ3eW29Y7T2GgWh\nMnLkSFGdkO2t2lSn9lantgphjPbiRTGOuP93Tk5OpKsQVmR7qzbVqb3Vqa1g/PaaInRc20NGIqk8\niYmJnDlzJtLVkEiCSv369SkoKHDaZzKZwIPcjrgwt1gsUWsrCgTZ3uBjMpmQyoGkquGuX3sT5hE3\ns0gkEomk8kRcM5dIKovUzCVVEamZSyRhJiYmhrS0NFJSUkhNTeWFF16w34Q7d+7kgQce8LtMs9nM\nzp07g11VSRUm4sLcYrFEugphRba36pGQkMCuXbv49ttvWbduHatXr+avf/0rAB07duSll17yu0yb\nBhYQpaWlAeeVRC8RF+YSSVXi6quvZv78+cydOxdQHmYDBgwAYOPGjaSlpZGWlkaHDh0oLCwE4Nln\nn6Vdu3akpqby+OOP28vKzc2la9eutGrVii1btgBgtVrp2bMnHTt2pGPHjmzdutV+nB49enDHHXeQ\nkpKCEIIJEybQunVr+vbty+23386yZcsAZbRgNpvp1KkTmZmZHD9+PGznR1L1iIzHvaRKEun+VKdO\nnQr76tWrJ3799VexYcMG0b9/fyGEEAMGDBCff/65EEKIwsJCUVpaKj7++GPxxz/+UVy8eFEIIcSZ\nM2eEEEKYzWbxyCOPCCGE+Pjjj0WfPn2EEEIUFRWJS5cuCSGEOHjwoOjUqZMQQogNGzaI2rVrC6vV\nKoQQIjc3V9x2221CCCGOHz8u6tevL5YtWyaKi4tFenq6OHXqlBBCiCVLlojRo0cH/6RIKo27fo2X\nSUP+BNqSSKKOSlgrCPY71e7duzN58mRGjBjBkCFDaNy4MevXr2f06NH2qeL16jkW7BoyZAgAHTp0\nwGq1AlBcXMz999/P119/TUxMDIcOHbKn79KlC82aNQPgs88+Y9gwJcJGw4YN6dWrFwAHDhxg7969\n9OnTB4CysjKuu+664DZUEhEibmapDjZVLbK94UWIwD+B8sMPPxATE8PVV1/ttP+xxx7jzTff5OLF\ni3Tv3p0DBw7Y6uj+YDVr1gSUF6yqHXzOnDk0atSIb775hh07dnD58mV7+tq1a7u03X25ycnJ7Nq1\ni127dvHNN9+wZs2awBoqMRQRF+YSSVXi5MmTjB8/nokTJ1b47/vvvyc5OZkpU6bQuXNnDhw4QEZG\nBgsWLODixYsAPmeynjt3jmuvvRaAt99+m7KyMrfpunfvzrJlyxBCcOLECftDtVWrVpw8eZJt27YB\nUFJSwr59+wJtrsRARNzMUp1mQ4Jsb1Xk4sWLpKWlUVJSQmxsLPfccw8PPfQQoHilqJ4pL730Ehs2\nbKBGjRqkpKRw6623EhcXx+7du+nUqRPx8fHcfvvtzJo1q8Ix1DImTJjA0KFDefvtt8nMzKROnToV\n0gAMHTqU9evX06ZNG5o0aUKHDh248soriYuLIy8vj0mTJnH27FlKS0uZPHkybdq0CeUpkoQBOWlI\nEvXISUPuKSwspHbt2pw+fZquXbvy+eefc80110S6WhKdRN2koUjbVMONbK8kXPTv35+0tDR69uzJ\n9OnTpSCv4kTczCKRSELDhg0bIl0FSRiRZhZJ1CPNLJKqSNSZWSQSiURSeSIuzKubTVW2VyKRhIKI\nC3OJJJS89tqcsOaTSCJFxIV5dfBD1iLbGz6OHj3KCy/8lWPHjoUl3/Lly6lRo4Z9ZqfVaqVGjRo8\n+eST9jSnTp0iLi7O7aQiiaQyRFyYSyShYsmSV5k8+SxLlrwalnyLFy+mf//+LF682L6vefPmfPzx\nx/bt3NxcUlJSKhXiViJxR8S9WeSamFWbcK4BOmfOU+zZ8zF16ypBq4qLjzF8+GHef/9G4uMbA3D+\n/CXatr2NyZOn2/MHmk/LhQsXSElJYdOmTfTr14/vvvsOq9XKgAEDaN++PZMnT6Zjx4706tWLvn37\n8vPPP/PKK6+E8rRIohx/vVn0+pnXA94AklFCMI4CDgHvA80AKzAM+M2WfhowGigDJgFr9TdBIgmM\nUaMmMXHiKjIyNqKZ5c7w4YeBw5w/D++/34lRoyYFJZ+WFStWkJmZSdOmTbn66qv56quvSExMBODO\nO+9kyZIlNGzYkJiYGK677jp+/vnnILZcItFvZnkJ+BhoDbQD9gNTgXVAS2C9bRugDTDc9p0JvOrt\nONVJSwXZ3lBSr149XnllHe+/34kLF5z/UwXy3LnrnMLMViaflsWLF5OdnQ1AdnY2ixcvtptS+vXr\nx7p161iyZAnDhw+vfEMlEjfoEeZXAj2A/9i2S4GzwEDgLdu+t4BBtt93AIuBEhSN/TDQJTjVlUi8\nowrmvLyWTvuXLWvpVSAHmg+goKCADRs2MGbMGJo3b85zzz1Hbm6ufYgcFxdHx44deeGFF8jOzpYT\nnCQhQY8wbw6cBBYAXwGvA7WBhsAJW5oTtm2A64CjmvxHgcaeCq9ufsiyvaHnwoULXHXVOazWWN58\nswU//hjLVVedsy/TFux8eXl53HPPPVitVo4cOUJ+fj5JSUnk5+fb0zz88MM8++yzXh8KEkll0CPM\nY4EOKOaSDkAhDpOKitfljNz9l5OTw8yZM1m4cCEvvvii001vsVjkttzWve3KokVzOX68jOPHH+SN\nN/bwyy8Pcvx4GYsW/ctjnsrkW7JkCYMHD3baN3ToUGbPnm03tbRp04Y//elPgHNYXInEGxaLhZyc\nHLu89IaeHnUtsBVFQwe4CeUF5/VAL+A40AjYAPwBh6CfbfteA8wAvtCUKWOzSIKG61v/uXOfwWy+\nnZSU9vZ93377NZs2rWbCBFc9hErnk0hCgb/eLHrVg03AvcBBYCaQYNt/GngWRYDXs323ARah2Mkb\nA58CN+KsnUthLgkaMtCWpCoSqkBbE4H3gK9RvFn+hqJ5Z6AI+N44NPF9wFLb92pgAl5MMN6Gy1UR\n2V6JRBIK9PqZfw10drO/j4f0z9g+EolEIgkDEZ8BKpFUFmlmkVRFZDxziUQiqYZEfNk4i4xVUqWJ\nRHuLiorYtmkTW/PyuJifzxVNm5KelUW3nj1JSEgIej6JpDojVDZs2CCqE7K9wUfbn57MzhazkpLE\n+rg4UQhCgCgEsT4uTsxKShJPZme7LSPQfMFm5MiRIi8vLyzH8her1SoWLVpk316wYIG4//77g36c\nm2++WezYsUN3+g0bNoj+/fu7/a9Zs2bi9OnTwaqaECL458FTfm2/1u7zJFQjbmapTloqyPaGkqKi\nImpu384TViu9S0rs/rMJQO+SEp6wWonfvp2ioqKg5AsFRppMVFpa6rR95MgRFi1aZN/WW9fy8nK/\njhvMcxCK8xnoeQg1ERfmEkmw2LZpE+k+FpRIP3aMbZs2BSWfitVq5Q9/+AN33303bdq0ITs7m4sX\nLwKwc+dOzGYznTp1IjMzk+PHjwPw+uuv06VLF1JTU8nKyrKnB4dwePLJJxk1ahTl5eVMnTqV5ORk\n2rdvz6OPPlqhDgUFBQwaNIj27duTnp7Onj17KC8vp3nz5pw9e9aerkWLFpw8eZKTJ0+SlZVFly5d\n6NKlC59//jkAM2fO5E9/+hM33XQTI0eOdDrG1KlT2bx5M2lpabz44osA/Pzzz9x66620bNmSxx57\nzJ62Tp06PPLII6SmprJ161beffddunbtSlpaGuPHj6e8vJyysjJycnJo27Yt7dq146WXXrLnz83N\npWvXrrRq1YotW7YAcOnSJUaNGkW7du3o0KGDW7fX06dP07dvX1JSUhg7dqzHF+N16tRhypQppKSk\nkJGRwbZt27j55pu54YYbWLVqFQBlZWU8+uijdOnShfbt2zN//ny/z8PixYtp164dbdu2ZepUx8Sz\nBQsW0KpVK7p27Wo/99GKfdggzQ5Vm3CaWWaNGWM3kXj6FIKYNWaMU/5A86kcOXJEmEwm8fnnnwsh\nhBg9erR4/vnnRUlJiUhPTxenTp0SQgixZMkSMXr0aCGEcBr6/+UvfxGvvPKKEEKInJwckZeXJx55\n5BFx3333CSGEOHXqlGjVqpU9/dmzZyvU4f777xdPPfWUEEKI//3vfyI1NVUIIcQDDzwgFixYIIQQ\nYtu2bSIjI0MIIcRdd90ltmzZIoQQ4scffxStW7cWQggxY8YM0alTJ3Hp0qUKx7BYLE7mjAULFojr\nr79enDt3Tly6dEk0a9ZMHD16VAghhMlkErm5uUIIIfbt2ycGDBggSktLhRBCTJgwQbz99tti586d\n9vpo22U2m8UjjzwihBDi448/Fn369BFCCPH888+LMbZrsH//ftG0aVNx6dIlJzPLxIkTxdNPPy2E\nEOKjjz4SJpPJrZnFZDKJNWvWCCGEGDx4sMjIyBClpaXi66+/tp+7efPmiVmzZgkhhLh06ZLo1KmT\nOHLkiO7zcOzYMdG0aVNx6tQpUVpaKnr37i2WL18ufv75Z/v+4uJi0b17dzFx4sQKdcRPM0vEX4BK\nJMHiYn4+vl5TJtjSBSOfliZNmpCeng7A3Xffzcsvv0xmZiZ79+6lTx9lOkZZWRnXXXcdAHv27OEv\nf/kLZ8+e5cKFC2RmZgIghODpp5+ma9euzJs3D1AiOtaqVYsxY8bQv39/+vfvX+H4n332GR988AEA\nvXr14vTp01y4cIHhw4fz1FNPkZOT4xSC99NPP+W7776z5z9//jyFhYWYTCYGDhxIzZo1KxxDuHGT\nu+WWW6hbty6gxJ/58ccfady4MTExMQwdOhSA9evXs3PnTjp16gTAxYsXadiwIQMGDOCHH35g0qRJ\n3H777fTt29de9pAhQwDo0KEDVqvV3sZJk5SY8q1ataJZs2YcPHjQqU6bN2/mww8/BOC2226jfv36\nbq9XfHw8/fr1A6Bt27bUqlWLmJgYUlJS7Mdbu3Yte/bsIS8vD4Bz585x+PBhYmOdxaa782C1Wjl1\n6hRms5kGDRoAMGLECDbZRnfa/cOHD6/QjkCIuDCXNuSqTTjbe0XTphSBV8FcZEsXjHxatHZTIYTd\nRzg5OdntMDonJ4eVK1fStm1b3nrrLbvJwGQy0blzZ3bu3MmZM2eoX78+MTExfPnll6xfv568vDzm\nzp3L+vXrK5TpKmwBunXrxuHDhzl16hQrVqxg+vTp9rRffPEF8fHxFfL447mjFfoxMTF2O3utWrWc\nzsnIkSN55pmK8wi/+eYb1qxZw2uvvcbSpUt58803ncrVlumuje7s1e7OgytxcXH23zVq1LCfhxo1\najgdb+7cuWRkZDjldWfecXceXOvmqV566qsHaTOXVBnSs7LYprlJ3bE1Lo70rKyg5NOSn5/Ptm3b\nAFi0aBE9evSgVatWnDx50r6/pKSEffv2AUq43WuvvZaSkhLeffddpxs/MzOTqVOncvvtt3PhwgUK\nCwv57bffuPXWW3nhhRf4+uuvKxy/R48evPfee4AibK6++mrq1KmDyWRi8ODBTJ48mTZt2tg11b59\n+/Lyyy/b87sr05Xf/e53nD9/3r6tVwjdcsst5OXlcfLkSUCx7+fn53P69GlKS0sZMmQITz/9NLt2\n7fJajraNBw8eJD8/n1atWjml6dmzp/3l5OrVqzlz5oyuOrqjX79+vPrqq3bhfvDgQYqKinSdB5PJ\nRJcuXdi4cSOnT5+mrKyMJUuWYDab6dq1Kxs3bqSgoICSkhJyc3MDrqOWiAvz6ha7Q7Y3dHTr2ZOt\njT2Gzgdga+PGdOvZMyj5tLRq1Yp//etftGnThrNnz3LfffcRFxdHXl4ejz32GKmpqaSlpbF161YA\nuynlpptuonXr1k5lmUwmsrKyGDt2LAMHDuT8+fP2tUR79OjBnDlzKhx/5syZ7Ny5k/bt2/P444/z\n1ltv2f8bPnw47733ntMqRy+//DI7duygffv2JCcn20066vHd0a5dO2JiYkhNTeXFF1/0GspXu791\n69bMmjWLvn370r59e/r27cvx48c5duwYvXr1Ii0tjT/96U/8/e9/91rWhAkTKC8vp127dtx55528\n9dZbxMXFOdVjxowZbNq0iZSUFD788EOaNWvms36u2+rve++9lzZt2tChQwfatm3LfffdR1lZme7z\ncO211zJ79mx69epFamoqnTp1YsCAAVx77bXMnDmT9PR0brrpJpKTk4PiERPx6fxyEk3VJhzt1U57\nnj5sGPHbt5N+7BjpNjfDIhTNemvjxhR37sxTS5dWKCPQfIB94eY9e/aErI2S6keoQuAGGxEsO5FE\n4trpwz0D1Gq1MnDgQL755pugtktSvZHCXFLtkIG2JFWRqAu0JW3IVZvq1l6JJFJEXJhLJBKJpPJI\nM4sk6pFmFklVJOrMLBKJRCKpPBEX5tXNpirbK/FGTk4Oy5Yti3Q13LJw4UImTpxYYf/GjRvt/vNG\nYvPmzSQnJ9OhQwcuXboU6eqEnIgLc4kkFOSuyKVvTl/65vQld4X+GXaB5gsWRgmn6g5PdduwYYMh\nI/+99957PP7443z11VfUqlUr5MfzN9RvVaFCNDCJJFBc+9PS5UtF4p2JghkIZiAS70wUuStyfZYT\naL4jR46IVq1aiREjRojWrVuLrKwsUVRUJIQQYseOHeLmm28WHTt2FP369RO//PKLEEKI+fPni86d\nO4v27duLoUOH2tPn5OSIZcuWCSGUaIo5OTmirKxMPPbYY6JNmzaiXbt29oiCWi5cuCBGjRolunTp\nItLS0sSKFSuEEEIUFRWJ4cOHi9atW4vBgweLrl27ip07dwohhKhdu7Y9f25ursjJyRFCCLFy5UrR\ntWtXkZaWJvr06SNOnDghhBBi4cKFFRZROHLkiLj22mtF48aNRWpqqtiyZYs4cuSI6NWrl2jXrp24\n5ZZbRH5+foX6zpgxQ4waNUqYzWZx/fXXi5dfftn+3z//+U+RkpIiUlJSxIsvvujz/H/66aciLS1N\ntG3bVowePVpcvnxZvP766yIxMVE0b95cjBgxwin99OnTncp9/PHHxUsvvSSEEOIf//iH6Ny5s2jX\nrp2YMWOGPc2gQYNEx44dRXJyspg/f759f+3atcXDDz8s2rdvL7Zs2eLzOvmDOzmJl6iJkaJSjZSE\nl8LCQrF+9Woxa8wYcc8fWopZY8aI9atXi8LCwkhXTQhRsdNnjMxQBPJM22cGImNkhofclc9nhBC4\n06ZNE++++64QQogzZ86Ili1bisLCQvHPf/7THjb2m2++EbGxsXZhXqdOHXv+vLw8uzA/c+aMff/r\nr78uHn74YSGE5xVxZs6cKf75z3/at/v37y/efvttIYQQ//nPf8SgQYMq5JkxY4bo3r27KC4uFqdO\nnRINGjQQpaWlYseOHaJt27aiqKhIXLhwQSQnJ4tdu3ZVyK9y8eJF0aRJE3Ho0CEhhBD33HOPXVBr\nH4xarFar6NChgxBCiLKyMnHDDTeIgoIC8cknn4hx48bZ9/fv319s2rRJCCFEQUGBEEJ5OKakpNi3\ntaF+9Vwnf3AnJzHySkPVzaYabe2dPmwYc5KTYeBAJr/5Jm/tP8jkN9+EgQOZk5zM9GHDvOaPtvYG\nimsI3C1btnDgwAF7CNy0tDT+9re/ccy2CMaePXvo0aMH7dq147333rMH4BK2ELjnzp3j1VdfBZxD\n4H744YdcccUVFY6/du1aZs+eTVpaGr169eLy5cvk5+ezefNm7r77bgD7IhC++Omnn+jbty/t2rXj\n+eeft9fNG0LjdbFt2zb+7//+z+lcuGIymbj99tuJi4ujQYMGXHPNNRw/fpwtW7YwZMgQrrjiCmrX\nrs2QIUPYvHmzx+MeOHCA5s2bc+ONNwJKdMZNmkVEhBsvp2bNmtGgQQN2797N2rVr6dChA/Xr12ft\n2rWsXbuWtLQ0OnbsyIEDBzh8+DAAL730EqmpqaSnp/PTTz9x6NAhAKdQv1deeaXP6xRK9IbAtQLn\ngDKgBOgCJALvA81s/w8DfrOlnwaMtqWfBKwNVoUl4UO7nFpuPAxqpOwfexqyi0vobbUyy5bOSAse\njx08lp1LdlLQqgCAxAOJjLtrXMjygTFC4H7wwQe0aNGiwn53As21ztqVjiZOnMgjjzxC//792bhx\nIzNnztR1DvQcU4s2/K42bKw2r3ouPaE3zKwr9957LwsWLODEiROMHj3avn/atGmMG+d8zS0WC+vX\nr2fbtm3UqlWLXr162V+oakP9xsbG6rpOoUKvZi4AM5CGIsgBpgLrgJbAets2QBtguO07E3jV23Gq\nU9ApiK72TnlwHJ3zFUE+3gzrximf8WbIs92HnfOtTHnQs8CLRHuz78hm3l3zyLBmkGHNYN5d88ga\n6Dl8bWXzQeRD4Pbr188ppK0aTlYbEvbbb791ih/TsGFD9u/fT3l5OR9++KG9DufOnbMvorFw4UKf\nba9bt65TSNg//vGPLFmyBFBeQvb0Em1Si8lkokePHixfvpyLFy9SWFjI8uXL6dGjh8c8LVu2xGq1\n8v333wPwzjvv6OpzgwcPZs2aNezYscO+SEW/fv34z3/+Q2FhIQDHjh3j5MmTnDt3jvr161OrVi32\n799vv56u6LlOocSfxSlcH48DgZttv98CLCgC/Q5gMYoGbwUOozwA3J8BiWFpcNHETeUwqAEUpGPv\nAQXpMH8PZP0CPcrh84uV88DQBrg6+NlmWnbvoSvAlTeyBmbpFsTByKeGwB09ejTJyclOIXAnTZrE\n2bNnKS0ttccVV0PgXn311XTt2pULFy7YyzKZTAwdOpTz588zcOBAFi1axB133MGlS5cQQrgNgfvk\nk0/y4IMP0q5dO8rLy7n++utZuXIl9913H6NGjaJNmza0bt2ajh072vPMnj2b/v37c/XVV9OpUye7\nEJs5cybZ2dnUr1+f3r178+OPP9rr5U5LHjBgAFlZWaxYsYK5c+fyyiuvMGrUKJ577jmuueYaFixY\n4PacuSsrLS2NnJwcunRRdMaxY8fSvn17+3+uMc9r1arFggULyM7OprS0lC5dujB+/HivxwBlcYre\nvXtTv359e5qMjAy+++47u7msbt26vPvuu2RmZvLaa6/Rpk0bWrVqZf/ftfzz58/7vE5G4AdgF7AD\nGGvbp436btJsvwKM0Pz3BjDUpTy7QV+uiWlcnsjIEAJERiMqvhhs5Fgf84kMzy8JfbX3yexsMSsp\nSayPi7Ovw1kIYn1cnJiVlCSezM72WU8i/EL9yJEjIiUlJaJ10IvZbLa/AK3OlJWVidTUVHH48OFI\nV8Uj7vo1QVgDtDvwC3A1imllv6tw9nYQd//l5OSQlJSE1Wpl9+7dpKam2odHqv1Qbkd2W11O7Y8n\nYdtHcP52AKj7EXRXFo2hCPglLs4pbrne8rt06ULN7dvpbrViiYXZNpv8H0+CuaSEJ2w2+TVr1lCr\nVi2P5RkBI/uHS5zZt28fAwYMYMiQIdxwww2Rro5XLBaL3dSVlJTkNW0gPXAGcAFFQzcDx4FGwAbg\nDzhs57Nt32tseb7QlGF7yEiMzP/WrIGBA+ldUkJePMxX1p9l3GnIKlZ+r4+Lw7RyJb1tCxL7w/3j\n7mbgm+9xNlaxwxfYRq+JW2GeRTnGJzVg1ZgRzJ3/rsdyZGwWSVUkFLFZEoC6tt+1gb7AHmAlMNK2\nfySw3PZ7JXAnEA80B1oAX+ptgMQ4aJdTyyqGtb8oH1WQg+/l1Lyh2uRf19rkTTabvO3B0aNcSSeR\nSLyjR5g3BDYDu1G06/+iuBrOBjKAg0BvHJr4PmCp7Xs1MAEvJpjq4oesEk3tTUhI4HLnzkxJrMe6\n2FiKbPuLgHWxsUxpUI/izp29vqT01t6yEyfw9XozwZZOIpF4R4/N/AiQ6mZ/AdDHQ55nbB9JlPPU\n0qV2b5PJL+Sxa10+d4xRllObWQlvE3DY5Meehp1bnc0s404rv4ts6SQSiXdkPHOJbv77XxgwQHE5\nCQbBssknJiZy5swZj/9LJNFI/fr1KSgocNrnzWbuj5+5RBJUuvXsyZzGjelttZJVrPitu7K1cWMe\n8mGTd+3wwcBkAqsVmjULetESSUiQsVnCjGyvg2DY5ENJICOQ6nR9q1NbwfjtlZq5JKJobfLTXs1j\n66rg2eQlkuqEtJlLdBNsm7krjz4Kzz8fuvL9wWSCH36A5s0jXROJxIFcA1QSFezdG+kaOGOEh4pE\nopeIC3Oj26GCjWyvZ2pEvDdWnup0fatTW8H47a0Ct4+kqmC08CZSM5dEE9JmLtFNqG3mAwYoxzBC\n1zCZ4NAhsC1gI5EYAmkzlwSFUAtZo2nmEkk0EXFhbnQ7VLCR7fWM0YS59DP3TnVqKxi/vREX5hKJ\nSlUQ5hJJpJA2c4luVq2CgQNDJ+QGD4bly40hRE0mOHAAWraMdE0kEgfSZi6JCqRmLpEETsSFudHt\nUMFGttdsawP0AAAgAElEQVQzVUGYV6frW53aCsZvb8SFuUSiYjRhLpFEE9JmLtFNqG3m2dmQl2cM\n84bJBN99B3/4Q6RrIpE4kDZzSVTQo0eka+CMER4qEoleIi7MjW6HCjbR3N5Q25Cvv97/8o1GNF9f\nf6lObQXjtzfiwlwiMSpSM5dEE9JmLtHNypVwxx2hE3Khjv3iDyYTfPstJCfrS68usLE1L4+Dn22m\nZfcepGdl0U0usCEJInINUElUEK3eLNOHDaPm9u2kHzvG5JISEoCi/QfZ9vbbzGncmMudO/PU0qWR\nrqakihNxM4vR7VDBJprbG2qbudGEuZ72FhUVUXP7dp6wWjltKqFnIvRtBB/FQ++SEp6wWonfvp2i\noiLfhUUZ0dyXA8Ho7dUrzGOAXcAq23YisA44CKwF6mnSTgMOAfuBvsGppkQSfvQI8ykPjqNzvpXc\neBhvhp0DYN045XdevJKmc76VKQ+OC2VVJRLdwvwBYB+gdu+pKMK8JbDetg3QBhhu+84EXvV1DLPZ\n7FeFox3ZXs8YTTPXQ4OLJm4qh9cbQEE60BwwKb/nN1DS9ChX0lU1ZF82FnqE+e+B24A3cBjeBwJv\n2X6/BQyy/b4DWAyUAFbgMNAlSHWVVHGMJsz1aOZlJ07g6/Vmgi2dRBJK9AjzOcCjQLlmX0NA7Z0n\nbNsA1wFHNemOAo29FW50O1Swke2NHvQI8yuaNqUIGHsaErcCRwCh/B53WklTZEtX1YjmaxsIRm+v\nL2+W/sCvKPZys4c0Aof5xdP/FcjJySEpKQmr1cru3btJTU21D2PUkya3jbWtdoFoLT8U9UnPyuK1\nhQvpUFzGPAvMrqPkHHAOskqV3/+OieEKTSxdo7Sv8uenarXHiO21WCwsXLgQgKSkJLzha2D7DPAn\noBSoBfwO+ADojNLTjwONgA3AH3DYzmfbvtcAM4AvXMqVfuZRyIoVMGhQ6PzAP/kEMjON42e+axek\npnpPV1RUxJzkZJ6wWj2mmZWUxEN790p/c0mlqUxslseBJiivde4E/oci3FcCI21pRgLLbb9X2tLF\n2/K0AL4MvOqS6ogRhDnoq0dCQgKXO3dmSmI91sXGojogFgHrYmOZ0qAexZ07S0EuCTn++pmr3Xs2\nkIHimtgbhya+D1hq+14NTMC7CabCEKaqE83tDVd8b6MIc708tXQpM386RsyqVdzAbXQhgzljxhCz\nahUz849V2QlD0dyXA8Ho7fVnBuhG2wegAOjjId0zto9EEtX481BJSEigd2YmJ0y1OC7MfPlG6Ool\nkbhDxmaR6Gb5cmWdzlDbzEtLISYmNMfQi8kEO3ZAx47+5YuJgfLy6BtdSKIDGc9cEhWoAjCaBaHR\nfOUl1YeIC3Oj26GCjWyvb4wizAOphxCWoNfDqMi+bCwiLswl0UO4hKxRhLlEEk1Im7lENx9+CEOG\nhE7YrlkDt94Kly5BzZqhOYZeTCbYvh06dfIvX1ycYvOX3VsSCqTNXBIVGE0AGq0+Eok3Ii7MjW6H\nCjayvb4xihANpB7l5Zag18OoyL5sLCIuzCUSV4wizCWSaELazCW6CbXN/OOP4fbbobAQIj373WSC\nbduga1f/8sXHQ0mJfCBJQoO0mUuiCqMIQqPUQyLRQ8SFudHtUMGmKrT36FHfaVSi2WYeCNLPvOpi\n9PZGXJhLogdVyDZpEtl6hItofqhIqh/SZi7RzQcfwNChyu9QXL6PPoL+/eHcOahbN/jl+4PJBJ99\nBn/8o3/5atWCy5flg0ASGqTNXBJVSEEokfhPxIW50e1QwUa21zdGEebSz9w7si8bC3/imUuqEEVF\nRWzbtImteXkc/GwzLbv3ID0ri249e0ZsVRyjRU00Sj0kEj1Im3k1ZPqwYdTcvp30Y8foVlJCAsoy\nZ9vi4tjauDGXO3d2uzpOqG3m//0vDBigxGjp1893+lA+kEwm2LwZbrrJv3xXXKHElpHdWxIKvNnM\npWZezSgqKqLm9u08YbWSGw+DGin7x56G7OISelutzLKli5SGnp2tvAT1hvaBNFl9IO0/yLa332aO\nlweSP0iBLIkmpM08zES6vVMeHEfnfEWQjzfDunHKZ7wZ8uKVNJ3zrUx5cFxQjhcKm7n2gXTaVMKg\nRtC3EXwUD71LSnjCaiV++3aKioq8F1TJerhD2syrLkZvb8SFeXWhqKiI/61ZwzvPPcfI1q342733\n8r81ayotcPylwUUTN5XD6w2gIB1lwGZSfs9voKTpUa6kcyXUmqpafnm593ThfiBJJNFAxM0sZrM5\n0lUIOVqTwCshMgnopezECXwZTxJs6YJBINfX10NDfSAN0j6QsD2Q9kDWL8oD6XM3D6Rg1sMdMTHm\nSh0zmqgO964Wo7c3opq5qq3+7d57I6qthpJwmQT0ckXTphSh2MgTtwJC+SRuhXGnbXW2pYsUvoRo\nuB5I0mYuiSYiJsynDxvGnORkdvfvz+Q33+St/QeZ/OabMHAgc5KTmT5sWKSqFlQqmAT6RtYkkJ6V\nxba4OLKLYZ4FMuYrn3kWyCpW0myNiyM9KysoxwvEznjpkvf/jfxAkjbzqovR2+tLmNcCvgB2A/uA\nv9v2JwLrgIPAWqCeJs804BCwH+jrqWBVWz1nKou4thpKKtioQbeNOhR069mTrY0bA4rwXvuL8lEF\nOcDWxo3p1rNnWOqjRa8mHO4HkkQSCoJtmfBlM78E9EJRdGKBLcBNwEAUYf4P4DFgqu3TBhhu+24M\nfAq0BCq80lK11VfMDiG3cyuYLMoNqWqrc+e/G1DDjEIFk0DzimmCaaP2RUJCApc7d2bKud/IOHeB\n7qWldj/zz2JjWXdlHWp17uzWLTEQs0Mo7IzdevZkTuPG9LZaySpWbOSubG3cmIcq+UAKpL2xsWYu\nX67UYaMGo9uQg00w2xsK11o9Zhb1MREPxABnUIT5W7b9bwGDbL/vABYDJYAVOAx0cVdooB4V0YYR\nTQJPLV3KzJ+OEbNqFTc3GEMXMpgzZgwxq1YxM/9Y2F7GBor9gZRYj3WxsfYOWgSsi41lSoN6FHt4\nIEkkkSZU79H0eLPUAL4CbgD+DewFGgKqKnnCtg1wHbBNk/coioZeAfttZiXi2mooSc/KYtvbb5Nd\nXILJArN3QGJNRZBH0iSQkJBA78xMaJ7J9tPw5RuhOY7FYtGt0fijCT+1dKl9BmjyrXnckZ7P1W2a\nkp6VxcwghSQIRDMvK7MA5kofOxrw59pWBYLV3ikPjmOg5j1asCwTeoR5OZAKXAl8gmJ20WLTNT3i\n9r+7gStPQPxaKG4HXAuJvyhCzoJi31G1VfXFg3oio2m7W8+eTExMpMaJE2QVw1UFFc/FosREsjXb\n4a4vWLBYfKd3CCl96f3dDqT83pmZWKlF6jjIyQlufYQILH+ozo/RtlWMUp9oae/5Iyco11omrEq5\nqmvtVb+AsFkmLBYLCxcuBCApKYlg8iTwCMrLzWtt+xrZtsFhO1dZA7hbRVGsj4sTAkRuPCKjkfLJ\njUcIRSESn8bFifWrV4uqwJPZ2eLRxHpibWysKLS1rxDE2thY8WiDeuLJ7OyI1a1TJ+WU6yE31355\nQsLy5YGVD0Ls3RvcuoAQn37qf77atUN3fiRVgycyMoRAkXnMQDDT9pmh7FNvgicyMirkxYvi7Mtm\nfhUOT5UrgAxgF7ASGGnbPxJYbvu9ErgTxb7eHGgBfOmuYKN6VIQCrY06mTF0M0WXjbq6EoiZxRT9\nr3kkISZU79F8mVkaobzgrGH7vAOsRxHoS4ExKIME1Sl8n23/PqAUmICHJ4nqUXHNb+eYUF6u26Mi\nWlFt1FZqcTTGzLYQ2aiNhiVENvNg5g1mmaWlFqTNvGoSrPa6vkebv0fZX9n3aL6E+R6gg5v9BUAf\nD3mesX28or7Aem3uXLo/c5C4s/ncMSa4L7AkknAjNXOJL0LlWhvR2CwJCQk8NGUKry+A3WdD51Fh\nLMxRO03cKH7m4SJQP/PqQjRf20AIVnsrM9fDGxEPtAXVT5uJVmFe3ZA2c0mo0LrW9h6ZR/mvlbdM\nRDwErqvbT9XHEukKhBV/rq/RbOaBoNjMqwfV7d4NdnvV92gxN77BdtbyxBtv0DszM2ATc8SFeXXE\nKIIHpCbpDSNdJ4nEFxE3s0Sz3S2wNSjN4ayiT0ItsMJ1fUPxUJI2c+9E870bCEZvb8SFOUSndliZ\nQDlS43NPVTCzRGNflgQXf5W8YPWZiJtZotHuVrlAOZZwV9croRY+0Xh9VQL3M68eRPO1DQQ97VXX\naWDgwLCv0xBxYR6NVKU1KP0RWEbRfsNFdWuvpHJEelWxiAtzs9kcdUPTyiyKbDSbuT9UNz/zQIiL\nM0e6CmGjul1bX+2NtJIXcWEO0WdnDPeiyKHESOe+KtjMJVUXXysDVU7JqzwRF+bRaHerXKAcS7iq\nGXQCEZjReH1V1q3zP09JiSXo9TAq0Xxt/WX6sGFMvP56r7bwSCt5ERfm0Uh1XYPSyNpvKEYY//qX\nMeohiSyqLfxPJ054tYXHNWoU0VXFIu6aGI12t8oEyqlRw0x5hRVRqy7+XF+jmVkCEczSZl45Apu7\nEVrUlYFO+lgZaP6J5jYlL7jREPUScWEOxtFm9HakUAXKMTp6BaYRb8hAKCnxP49R+nI0EopFjoOB\nagsfpLWF41gZKOsXxRa+sX5Dtgag5AWrz0RcmCt2N3OEa+F/R3Jdg7L1Vfl0v8N3oBwh1hAH/O3e\nPC7m53NF06ZRKeg84XoevwS6GOCGDBeKzdwc4VqEh2DGM9e69eXGw6BGyv6xpyG7uITeViuzbOnC\nfZ+otvCCy57TJAA1Tp+OqJIXcWEOkddmAu1IjgUnMmnTBZ7wEcJ3+rBh3Co2MZECer5ZYr/Q0SLo\nfGnm7s5jwWV47IK+G/LMmdDUO3CK+N8aZYSh98Eb6b4crYRqkeNgoDo89D8HR7Y66ubOFv7EG2/Y\nlbyH5uTx1dqqv06D07p27dpFdt3EP48dIT6pgVgaj0jsa1uXb4byW12XdE0NxJ/HjnCbH4S49Vbv\nxygsLBSzkpKEQDmOuu7pUs26p08nJYnCwsIQtNAznTvrP/fvvON9jc4/jx0h1pgqdx4DXQP066/9\ny+OLlmSLMSSJ9XFxTuu2ro+LE7OSkjyu23rVVXIN0ECYfvfdotDHupiFIKbffXfY67Z+9eqA1iz+\n+GN9faFHD/19hkqsARoWIq3NhMM/dNumTaQfO+Z1QkH6sWNs27Spkq2JHEc/2UwPETk/22BRVFRE\nD7bzBpGZyVcdibRbnze69ewZ0JrF4ZZrERfmRvBVDUZH8mWC2JqXR7eSEp6t41nQpZeUsDUvz4+a\nhxdvbSwqKiLh119JAE5oe9UR53TRMJlq26ZNZBHYg1f6mQdGqBY5Dgaqw8OddWuzLjYW9RFeBKyL\njWVKg3oUV8IWXmVegBoBbUfa6cMm5glfwvxifr6uB8bF/HydtTYW2zZtov7ly7wbD99fB+wGUpX/\nahvghvSHrXl5TKbEq/dCekkJc/Ly6J2Z6ZQ30qPMaCVUixwHi6eWLmXNmjXEAMm35jG8Vz51rw+O\nLbzKCHMj+JmHoyOpD4zHLsD4AB8YkcbbA2trXh63C8H0BlB4G/A9sF7574ZDjvO4pUYNw0+mqsyD\nt2ZNcyiqZEiCee+GapHjYJJpe3BbyWTah3DllRGrilsibmaByGszgdrEtPjSzI06a9Sfc++tjRfz\n8zEDZ9UdNwJ9gFugoaaXfRIf7/U8GoHKDPkj3ZejFfvcjcR6rI0JvikjEujtC1VGM7dYLJhM5ojW\nIRyTgFTNo4bBNI9gzZxUBVv703BsK1y0jTzqfgT3nFaU9M+AH6+91vA3ZHpWFpveDGykVlxsQfqZ\nB4Y6d+ODJZtIHpPH1RjLrS/Y7VUJpzBvArwNXIOio8wHXgYSgfeBZoAVGAb8ZsszDRgNlAGTgLXB\nqW7ocJ0EdF1MPrfl6O9IvoSi+sB4znqRixRwMw4/861xcWxt3Njwmoe3NqqmqrziEt61wIt7FC29\n0Un4sRSuAzrGxnLTv/8d9rr5S7eePZlEYzLx/8ErNfPKkZCQQPrNytwNK/Clj7kbkcLIcYq8cS32\nV1nUAQ4ArYF/AFNs+x8DZtt+t0F5/RUHJAGHqWjOcfKdTEszlm8uCHHllf6l79PHd7qzZ4WAQhHH\napHEGHF3ywwxa8wYsX716rD7l6t06aL/3P/nP579wLV+9J4+vvzoK+Nnvnu3f3l80ZJsMZok8amL\nn/mncXHiaS9+5tddZ6y+HI0cPhxYPwgXIMSZM77TrVmjrw29egXHz1yPZn7c9gG4AHwHNAYGAjfb\n9r+FEtt1KnAHsBgoQdHYDwNdgG2eDlBYqKMWYSYUT95OnQASKEHRPB6cAA88EPzjRAJ15DELxW0v\nvST6Rh5aDrKUgxQxYqUyUku+Op/0gcpI7SE5AzSkRMM5NKJm7u8L0CQgDfgCaAioDsMnbNugjKiP\navIcRRH+brFYLBw86GctopRDh8Bo8cyFAOJz6ZvTl745fcldkes9rReeWrqUh/buxbRyJcmMoQsZ\n/PGa2zCtXMlDe/eGNFTB4cOhKFUN1/AG8Tet5Yk33qB3ZqbXB9Lly5ZQVMSQhGqOiFGFudHb688L\n0DrAMuAB4LzLf17Vf3f/5eTkkJSUhNVqxeGUbAYcJ0192RDubbBQWoru+oCF06d9p3e8GHPejmR7\nCy7kQrsxrOM8JMHOJTvZ9+0+bu5+s9f6Wyzuy0tISKBGrVpYuRsrZv54o4UateDLL7/0WR895bvb\nBgtZWSBEcM+Ptj4nTzq2veVXbkz/6h+t2yrBLv+LL9TyjdlesLB5Mwwc6D296tjhq/zffvPcXovF\nwsKFCwFISkoiGMQBnwAPavbtR7GnAzSybYNiapmqSbcG6OpSXgUblJHsYyBE3br+pe/VS1867WfO\nnMDrGCx+1yqjYiyMkRlu077xhv5rpaYbOFB/XSpjMw92/9GWCUIMHqwvX9OmxurL0ciRI8aTCVpA\niNOnfadbu1ZfG265JXyxWUzAm8A+4EXN/pXASNvvkcByzf47gXigOdAC+NK10NwVjqE98Z6H9tGC\nEW1o3lDPf2HZTt15oq2NwURv241qIogmouEc6ukP4fYz1yPMuwN3A72AXbZPJor3SgZwEOiNw5tl\nH7DU9r0amICbp8n4JeNZl7SOdawD83iIN25MEj3oF3SWENZCH7krcu3nv6xrAewwOSbGHEhk3JDg\nrR5+8qQloHxGe3DoveEuXbL4XfZrr83xO48RqGh+CA5GFeba9urpnzV0vpEMpzDfYkuXivLyMw3F\ndFKAMsevJdAXh485wDMocwD/gGKeqUBBqwJ7zAvSC6DB/EDqHzL8FSbRtBTc6x++7jj/LYB6At5K\npONXGcy7ax5ZA93PQg21gB061PF7377QHstfQtX2o0eP8sILf+XYsWOhOUAUYlRh7i/VLmoizR0/\ny8oiV43Kov9mNweQJ8TcCFzuyLMPrPUoyCGw+iYmmnWnbd/e8dtosk1v2xMSzLrLzF2Ryy139yCm\nyVmmzbw/sIpFEMcL6OqB0dsbsen8iQcSFe0QYGsinB7Hiy/Cww/ry//aa3MYP35y6CroJ4YRzDoY\nO3gsO5fsrHD+Q8FHH+lPqz2HRjufwbKZz5nzFHv2fMwvZ86xofwQl82lABzZtZJTg9vQouk1nD9/\nibZtb2Py5OmVrLUkVATTzBIsIqaZz7trHi03ZsArHcEyD4qz0BvmOhxD09CZWSx+1iT4ZN+Rzby7\n5sH8DOVjO/++CEzAWgLJFLXC3JfNfNSoSZSUlHG8+Dsup5Xa49pfTivnZNl39OmzkeLiMkaNmlTp\nOocaaTP3jhFfgIaErIFZjB+4FgqetwsSPSfIqENTowkfX2QNzIJf1iofHYI8HBhZMw8W9erV45VX\n1nHi+O8q/FdaBu+/34m5c9dRr169CNTOGBhVmIOjXxpRmEc0aqLSCLN929MJioahaSA282gjMAFr\nDnH54WHVKuWlbJs23tPVrm32WVa9evV45pGX+H+Lx1DcQRnS1d0LV5Y2iipBblQbcqhMsGazOSTC\nPFhE9AXo5s3O255OUCSGpv4KlkAEkZGFlztCXV8jno8JExy/H33Ud3q9N3CfHn3oG1+X1p/B796/\ngj81qUGntoJCIwYqCjOVEYLh8g7SY1Y1omtiyPjgA9DaVD3dzEYemq5apXxHk808vFgCymUUwe5v\nxN5Tpyy60i1aNJcra8Xz665HOLf/DE2ueYjjx8tYtOhf/lcyQhjNZh5qE6zFYjG0Zh7xxSm0eDtB\n4R6a6hUmP/2kP73JZBwhFQjh1MyNeJ4+/th3moICfWUlJPyOqVPX8d57ij/m1KnP8e23d7Np0+pK\n1LBq4I8QrKwJNlCTjBFt5pH3M/fDpqodmiYuM8bQVL0QUe1nrpNw2syj7dwAnDsHett7//2Pk5LS\n3mlfSkp7JkyY6iGH8TCCzbwyJlh/TTLa9uoZiVdDYe7A1w2sHZoO6Bjaoam/fsV+zwCNz2X+Kt9h\nZ8OJkQSokeqil5ycSNeg+hGoCbayJhlpZnGLBVWb8b30mmNoWqNGaIem/gpz/WYWC8SfBPN4DqYX\ncBAl7KzJZPI6+9IIBO5nbva7/GicDayYWCwY1WMp2F4elhCtielvP1u3cR2J9ZtwfPteRGdlnzsT\nbGVNMhaLhR49zIDUzH3i6yKGc2iqV9P238wCNHhdiUdjGxIWtCpg/geRj00TKm1YGyFT7+IXTZqE\npi6hJNwz/vyhqsaAUYPG7b1pL6I+1FgOVyyKd2uCDaZXnJwB6hYz6gjISENrf4W5fjOLOYDahAdf\n5z+g6xN/0hEhM2kd45eMJ29ldEfI9IRy85ojXAv3LFnyKpMnn2XJkleDVqYRbOauQePKB8HFCz3d\nmmAr6xWnba8RzSwGEOYOjDS0DoWZxc7psUo8lBCFnQ2UkAjzBpqbzccoJJxmllCEnTWSZj5nzlOM\nHt2NBx4w88ADZr7/PpfWreHw4aX2faNHd2POnKciXdUQYGLq1OeYOnUdtWs7C27VKy7+qxr2+88f\nrzi1j0ozi1ss9l/+3MBGmcDivzC3QHG2Eg9lfgaNP/YcdjZa41w7cVmnrx7O5zCUIYVDZXKIiQGj\nzCNQTQoZGRsZPHgjw4cri6QOH36YwYM3BmWiXaj8zP25t8cOHkviAYdiVOdbR9A4TyZY1SuujZ8T\ntgKNzRIui4MBhLkxCbk3S7ESG+W2VPdhZyNh4/TV5nXrAij0XH+nm03vKCSUmnkoTA5gLM1cNSm8\n/34nLlxw/u/8+aoTA0YNGtd6ixI07s/JvoPGqV5xpl8Cn7CVq8MBTe8EoyqkmfsX70DFnxMQSg3X\nX28WV5uqp3yhEjje8NWG//43gEJLZzpFaNS7+EUwNfNwmRwUzdwcjCoHBVWg5+W1dNq/bFnLoAhy\n1YYc6RFk1sAsHshWgsb17OrbI0z1irviiueAmh5NMq5obeavvKK/fr76cpUItKVy9qzy7Y8w15tW\n1XAHDBhG48aNdeWpUSM03izuLpqaT3Wbqlu3FgDFxccYPhzef38pDzzwGUDIg4mFajiYNTAL7siy\n/dZ3/GAK81GjJjFx4ioyMjZSp45jv2J6OGzXVCsb28dImrnKhQsXuOqqcxw5Esu8ec25774jXHXV\nOQoLC4OilQdyf4UC9d7Scw3uv/9xpzygmGRcPeXcofZRveG6tXlCjQG6n4WkJOXXG28Ev/RANFx/\nbkr/hbnFaZ+aLxw2zshgCShXMM0s4TI5KP3GUqkygs2iRXM5fryMnTsf5Isv9vDLLw8GbaKdxWIJ\nyQiyMsIvkHtXL/6+Iwi3mcUQmnlcXPDKCoaG68/JVTtPZTVzVeBMnJjB8OE7nDTIcNk4I+0a6o9m\nfuEC9OsHn32mr2z1/D74YFdycg7a9wfL5ADG1MxVk0JBQXv+9a/KT7TT3l9Hj/7GNdcUhnwE+f33\n0KwZxOqQVoqpK7isWQMdOwaeP1zrAxtAmJs5dCh4pQVjSO3PTamm1e+qZHbapxVg4RA43giNMDcH\ndHxfmvmxY/D55/7VRDU5/PBDLPPnB9/koPUzN0pQNdWksGmTY59ek4I7wmWy0nLjjfDSSzDJS5Gq\nouSPIqY37a23wp//DHPnmiku1l++SjUyswSGp5s9GEPqQIR5sGaEqQLn++9jmTq1BT/+GGsXOKEm\nWJ3u9OnKH9/XwzEQDUw1Obz/fvBNDmBMzRzgwQcVgRQMIuUlc/y49//9sZm75vEn/cmT+tPr9Umv\nQt4sloByvfMO7N/v/r/KvsUP5OkeqM3cFVXgvPpqaAROONizR7tlCaiMUAhz1eTQurV/Xgx6ef99\n8Ke9gWh5gfDSS1BUFLzy1Ptr7tzfO+0P1gjS3b2kt/7+9At/7nOlThZ+/ll/Hue8oUePMP8PcALQ\n3qKJwDrgILAW0F69acAhYD/QNzjVdI+3p6Sq4R46pGi4Vqt+DTcQm3mgM8JcL7QqcBITFYEzbVpw\nBY43gtXpXNt5+LD/x/dlZglEC3aN7fPOO/7H9gmm/fPo0eCVFW4uXLjAlVcWcuhQLM8+G/oRpK/+\nEEozC/i39icoCs1bbym/jaSZLwAyXfZNRRHmLYH1tm2ANsBw23cm8KrvY5idtgIZxrhD1XCnTVM0\n3G+/1a/hBmJmyc/3ndadzdz1QqsCx9VtKhxxrr11OtV91H/M3HWXvpShNrOofPWV8n3PPf7n9e1l\nY9ZdlpEXLvbFokVziYuLZdq0B0lJCf0IUq9ADJWZpbQUwKzby+rxxx2+6EaaNLQZOOOybyBge+7w\nFjDI9vsOYDFQAliBw0AXfyoUrBnCqoZ7+rSi4Y4cqV/DDdXTXY9mHki5KpWdvOHthklO1l+O6w21\nY4f/dfF106jC3N/RxMmTFdee9YdgukxGqzBftAhmzXLcXzVqBN9k5YreEVEgipgeXn/dv3poPfTC\n5SYc4YUAAB6aSURBVM0SqM28IYrpBdt3Q9vv6wDt4PEo4GMmgYVrr3VsBa4BOlOZcLmBdgh9Q0GL\n0z5PF9r1Rj9yxPvNH4zp/946nT/FOtfTojufP5q5mtZf4XrNNf6ld8X38Sy6y9L2nfx8/UvORZqv\nv4bz5x/n1ClFx1OvRbBGkO4e0L76g3ou/Rmx+W+qs+hWHrTCPJr8zG1RN7z+X4GcnBwgCbBy4cJu\nIBUwU17ucM5Xp8+6bqs3jBDu/3eXfv16aNtWX/qyMostn+/02okiZWVmYmI8pzeZnOuvDtvcpVdu\nbMe28lLRc33ef/91++SNjh0zvLbPMfnBuT7l5frTr1oFAwaY2bABSkosxMc70u/a5b58X+dTCDPd\nu8NPPyntzc72nF4xx6nDXn3lu/YfX+nd1d/d9Xr1VYvN28K/8m+4wZG+WTNITzfz+ee+6+/PtiIE\n9dVn//5djB8/2Wf5P/2klqfw448WLJbg1Bdg61bn+oLF9n7Bc/6DB5X/V6xQ0uupT40a+uqjfUDr\nPZ+nTjn+37zZwpVXei7/5EnP5VksFhYuXAhAkjq7spIk4fwCdD+g6tONbNug2M61j+Y1QFc35Qkh\nhFCeWUJcc43j93//K3yipv3f//SntR1SFw0a6E+/YoWj/AsXvKetX9+5PiBEdrb7tLfe6lzvzZud\n6/TCC38Vo0Z1FZMm3SwmTbpZjB9/o9iwATF+/I32faNGdRUvvPBXj/VxrcvixZ7r7po2Pd2x/9VX\nndN+9lnF9HquwQMPCDFnjhB33inEe+95T3v0qFJeYaF/11dvfTzV//Tpimnbt/e/rUII8dNPzulb\ntdLXBn/IzdVXp59++km0aHGlOHr0qM8yp051lAVCPPRQkCprIz/fub4gxJgx3vMsWuTczvJy38e5\n5RZ910lb7oYN+q7viBGOdL/+6j1tVpY//dez4hyomWUlMNL2eySwXLP/TiAeaA60AL7UU+Dx4zB8\nuL4XiSrehl7l5fDRR/rL0hKomUV5SeIZd8MpPWaWTp0q/h+K6f/+2PZ+/NHx29X0EOiwUR2O1qwJ\nly/rSxvuGPjujhdoe8Phl67Xpa8y0/KDbRMOxMzial7RU6dAzn8gtnsj2cwXA58DrYCfgFHAbCAD\nxTWxt20bYB+w1Pa9GpiAlyfJLbcAWCgpUbb/+1+YMEF/5b2dpB07oH9//WVpCa0wtzjt0yPMd+6s\n+H8oJm/40+m0QsL15nM+fxbdZQqhtLtmTdi3z3daMKIwt+gq5/x5+O67YNTIO56uaWUiSartVc0B\n4RBW3o6xbJmiCOpN75rG133rQL/N3J93acFCj83ck2NZHw/7n7F9fKI2+MwZpYMkJIA/bqqhOkmB\neqgEopl7aoOeOgR7+v9TT8GIEfqO/dtvjt+uHbwyL3RMJphvW4ioZk3o1Ut96DvjTpj/+mvlX3D6\nIlia+WOPwb//7bwvFJNLxoxx3i4vV+67ykzLd21vOISVN+HszgOurEx/zKeLF5VYP/fdB8uXe08b\niGbur1woKFD6QoMG+o5lP6Z/yYOL1u/aZILVfsb+8XZiKyNQAn0jrhVwntOanfZ5aoPr6MDTja6N\nN1LZ6f+HDgU2K9G7MDf7XU57mxPS3/4GL7zgPY/2RmnYUNF4Q4m7G9P5Wpl1lXPpUjBq4xvX/qXW\nvzIjO/X6qi/wAtHM/XWj9XYMdyPp8nLfQlTtb8XFSpwf5eWpN8w+21pSokwG08oQvROeDtr0sS5d\n/HMFVjGAMHf8TkmB+Hj9Goq3dJUR5nqis6loO9If/uA9rbvO5alz6H2gqJOj1qwJzvT/QG7MYGnm\nhYVK3lGjPJftut/1RvE23Vrv881bv/JHM/dWTqTiuGjr7ynsxQsveB/ZubbX3z7jy43WX5u5u/P/\n22/6NfPFiyHL95oWAKxd6/3/2bOhSRP/NHOVbt2U76NH/YuXrmIAYW6x/65ZE+rUgbw8ZcjsC2+T\nP8KlmftzHMXd0OK0T68w9yQY1MlRN90UnHgj+u2HnrnzTsfvESMsuvL8+issWKCcTz3vTTwJc29C\n8sMPdVXFK65aLMD27doti/3XxYueywlFqFY9vPaa87Y27MX06S04fDiWpk3PcfCg5yefq83cXzNL\nIC9b/RXm/sSi2bJFb0qL00jROQaRguKSGJiZRW1joH0josJci9qgWrVg2DB49FHfeZ591vN/ldF8\n/NHMK+vw7+lG0CvMgz3935MXSWJixX2q26tr3b7/3vG7XTt9x9XefHo0KvWYrjeKt+uht094OtdX\nXAG//KKvDPA+EnBXl3AEZHroIedtbdiLEyf2MHWqMrLLzNQ/svOlmQdj2T5/hbmec6mm8RWR0RPu\n+rZ6XQPRzKuAMDc7banDZL0duzIvEN3xyy/+XVz/j2N22tJrMw9+PdzjyZZbs6byXaOGEncCHO6J\n3q5VSopZ18xGT4GSfPUDf7TCypo2GjTQExzLbP/lTov3VpdwRdfT4hr24vTp5/j003WUlnoe2bna\nzH1dg2C40fprM/eHjRv1pjT7TOFOmOu1mavpolKYaztvoHa4YE3/V+nYEc6dC26Z3vDUTtfRga8b\nPVg22K+/dr+/lrJwE9OmwaBBzv9566x9+8LvNHJBdUN1xVWYT5yofPtrM/emBel94Hk6Zlycf5EO\nvWnm/oz+Qok27IV67srL23P+vP6RnS9h5e/L1mDYzNUyPOXbuBE2bPBe70BwJ8xvvtl7HrUvqOc/\nioW5Bah4QfS+rPJ0gwaqqZ454/63J/RqU46yLE77PXUofy+omv7ee/3LB4qA7t1b+e3JN1/VzMvK\nlBfVWjwJaIAtWyzExDg8ff72N/fp1M6vXjf1AeDpZlSFgusqVd4mG1XmgdemjeK2+dNPvlJa7L8O\nHPCcSivM/YmJr1LZoGru8HYdtQTiZ17ZNQb8FeaqYPT0oPnf/7wezgMWnyncRW88d867p4xqVlQ9\nydR7zV8MYGZRcBdYSg+ehvCB3rjaerh7weGK3hvQ25DbHZ5s5uoLFlfUl8FvvunfcUDp+Hl53tOo\nmnlpqWI71uLJnVFrZ7/ySuX722/dp3XVVNWF3j1p2qNHK98DBjjvD4Yw117T3FxFe9y7F9LT3Qvz\n//s/9+V485DQXt/09IrH9UYwgqq5w19hrrrS6TV1qS9brVb/3Wj9NbOo/cbTfae9zzt39nl43bjT\nzMFhmnSHa4TFQBbAAENo5mZAcUnU4mm4D84n6skn3acJdBirLVuPL7D/M8LMftcDHA8WTy5L773n\n+O1JuHiiuBjq1/eeRquZu+JJCNSooQ1WpLBsmfu0rm/0f/97z8cDzw+1rl1h927P9fGXoiLsk2qa\nNHFvZnHWpMwkJPguV9s/77tP+db7oqwyU++9oVeYqzRvbgb0zw9RX7b+9JPiRnvsmH43Wn/dPNV+\nk5rqO4+zN5I3zD5TqILZtU7eZjRrZZ+/10CLYTRzd0MLdxrfd98pN/zixcp2Wpr78tyZKfQIXu0T\nW48rWyDTe995x/k/dwqWa/0feED51jOkXbxY0SRNJu/ucSp6grGpmrk/wtxVS+7d2zEhyBW1XWpZ\nTZoo354EnOv5GTnS8dvTNHntNVD9hd21R3tNtef79793r5m7KiJ6Hhra+l9/vfLtaYQTDG+Qm27y\nXSe9DxO17uq50fuOSX3Z+sADysvWyZPdu9H6azN3d6+rfc9TrKdQ+fmrD2l/ytdq5tpBiqcHkSci\nrplfdZVFqYimJqpwefDBinnUhQ66dYMpU/x7kul5eaWthyctUova8dSHiyf/Vke5Fkwm55mN7mYt\nerKZ630XoNq1Fyzwndbbqucq3jRzTyOY8+e14WTh/vuhaVP3aV3jZKia+eefu0/veh7+pVHuPD3A\ntOe0c2e46irPGr7KiBGO3/XqKW06eNA5jXozKoqFhWnTwBa11CNazbx5c+XbkzAPhjfIhAnQtq33\nOul1JlDr7gh1rA/1Zav2wa3XjdZfYe7rugYmzC0+U3jSzPXkAWezkDfrhDsiLswfeqjiUnExMYqg\n/uEHpXFajUHrIteqlWctzN3T/ZNPfNdJ7RitWin10hvwSdVcM10X2PNAo0aO3+5Ws/ckzH1FE3TF\nXdmuqB3v4Yc9p9HazMH5IerpxnEdbcXHw6pV7tOqDwm1fF9xKVwFX+3ajgkxnh6oWg26Zk1o0cJ9\n/9H2He2NZjIp7wtc+6sq3ISAefMU+2jPnp4fXOB8fdXFWTy9/wlk6v3bbzvX86674IsvlHZXNpaK\nek4CDbClnl89o0ZfXingXnDefbf3cj2ZYfWOtJ94wv1+9dy4U6I83YvaPqZVHvwl4maWzp3NXHWV\nYzsmRnlZ9pe/KMK3bl3niQ4//KB8m0yKpvHuu+49X9xdlLFjfddHveFvuEH59vXWWz2OepN5mpXq\nqI8Zk8n5Rv/HPyqmD5Yw95a+YUOHvRbgz3/2LIBUwaxeK+3N4OlN/bffOtvMe/ZUvt0JE1czi6p5\nuxtq/vyze3OHaprx9D5N+3CJj1eE+XPPuU9bo4b7F1G3317x2OrNePEi9OhhBpSH9YkTngWnO28W\nb3jyBvn73917g4wcWTHo2BVXKP3B3xhIrqjtbdXKHFB+9VrrEeZq/127Vpkq7w53wtzXKMPVNKbi\ner2cYwOZ7b+efloJDOh6HPW6qu+2MjIcbfjsM/fH1N7r+mejViTimrkru3crnhV16zqEr6o9Wa2K\nG52Kakpwt2Cwa9ktbffA0qXwxhue66ReZLWDTJzofV1S9Ti+Oo9Wszh1Cjp0cGzv3VsxvbaDajVm\nf1d09xYWoVYtZQSk8vvfKxOm3D0AVEH4l7849qnCoqjIvVZ5443O23XrKnnc1clVM1dx9zLT0+xg\n1b3Sk0upeqM99phyA12+DB9/XDGdEMr/2tGTyiefVOxv6gu0AwccL0tr1VIeLvv345ZAXtBrp96/\n8IIy9f76689x7px3b5Dp0x2/W7TwLzKpO1RhPnmy/jyTJimeQePGOfqXJ/Oc9t7Vjpw+/dR9eleb\nvTZGkidPFbUNjz2mHE9ti+t9rL33nrHFgv373x2xpFzXt3WdvbxjhyJTnn/e86jU22hAz8haJeLC\n/OuvLU77UlIcNkT1qaheRG2ja9VyuMi5O0muJ2jdOuV7+HDvGroqRDdtcuzr1cuzhqUep6jIseq7\nuwvgqI+FFi2U+peXKwJOO/3dXf2ff97xMBo+vOIwX0v9+s5uhh98oD9CX1ycEq3N3ehC9abRajRa\nW7sn27bF5Ul4442wa1fFdOr51bb7z392X6ZrRDqtqWvAAM/CXE2nvix9/nnl213wJE/asvp+Q1tP\nbRO1/blzZ8/amCrM339f+V6yROnP3l5Caqfer1qlTL0/ebKMBQu8e4Nop53fcYdv06GKr4eidvk4\nX7HZX3lFCdPx+usOt1M9ES61IWm/+MJ9GleT2zPPKO9nQFEG3Pn7q0JXdXJ49lnlXvTkCQVw6JAS\nz3yqzcQ/ZIhy3dyVq6K2cehQRYn0pvR98IHjd3Kycq1mzPCc3pWIm1m8DTHr1FFuTvXpqj1RV1/9\n/9s78+CoqiyMf1OIKIyoLFFAMTCsKjUoDoKCgOKogGjpKBBhjCCLI5ujJcVSCC44qCgIAiUGNyxj\ntHBKKVAcJDggixHIsMgSnKBDMGJYAiFAkj7zx9e33uu3dNPQId2d86tKpd/t827f8+69555733v3\n8L/XXtdeeHlZXhhjftttNFzm2G+5QoTlfOwx68karxsXpvPfdRf/AOqem8ulC6eBdg5G331nNUrn\nRvx2cnPZcKZN48DYqZP/zTgvj2DECGDmTO98TZkNEyda08lIW9Uahg4FZsxw/7Yx5nYva8YMzhY2\nOGJV2duBc4o9fToHd68pvHP91bQJ5xMP4TylggJeA7+nJOzP4D/yCAcMrxubRoe+ffm/Xz+283A3\nrGvXrouuXd2v3h8/Hn5TtYsusj6npfFdBPt9BWPQOnYMPW/2bO/8nAarQwf3Bl6nw8yZ4a91air3\n2DccO+Y9aJjdBg2BAF9OW7yY91K8djM1Opib2TVr8py5YZ72dN4DSk/nQxL2x4WdMy4zOKem8rFZ\nr1UBEfbVPn2YJ8B2MH48b+xHE32tKhARkW7dIsfx/OQTkXvv5eeePd3x906d4vGmTaHnbd1qyT73\nHNPGjrXSVqzw/r1Wrfj99u08vuACHhcUeMt/9plI797W8bhxIoMHu+X27mU+Awa4vxs7ViQtTaSi\nwkp75RXKP/tsqGxRkUijRowJascvLuH69SKNG4sUF7u/a9pUJD8/NO34cZHUVJHly935163rzkOE\ncSZvv906vuce/9imJ0+KtG0rkpUVmr5lC3+jsDA0fcECkRtu4HmGUaMo266d92/cfz/rwcmsWe74\nkBMmiAwcyPitJSVMKy0VqVXLO28RkX79RCZPto7tMSLtBAIiffqITJnizuPttylvr/PcXJEGDdh2\n/TC/89ZbIiNG8HOvXu6Yl/Yy7dkT+l3//iKjR1vH48db16VTJ5GlS3kOILJ4sbsM778fmv+KFVaZ\nIpXb+ffqq27ZvDx+N2WKSNeu/Jyby+OOHUXKykLl+/Wz8mvfnucbjh1jf1mzJvQcEzN0+nS37Lff\nepf76FF3WceMERk61Dp+9123joacHJEmTdz5TJoU2s9nzqRdE2E7HjbMXh7/yG1VhVRUUNF169wX\nyE52NivUfmH37g2VSUsTueIKyhq2bBFp08ZtxI4cYcDmnj1DDce+fbzII0eK3HKLlW4Cs06Y4F2+\nyZNFbr3VOj50SOTqqylv76j5+cxn7Vp3HsXFIjffLJKeTmMqQsMzcqT3by5ZQgO9bZuVBoi88IK3\n/KOPMmisvROcOMFznMZcRGTZMua/ezePi4oo27Gjd/4lJbzWTz1lBVqeN89bVoTXoGFDkVWrrLTN\nm72NcyAgcvfdDOhbXs40E1D44EHv/Pfv54A0YwbPLykRWb3a2+AWF4t06MD0evWYduQIB3E/9u0T\nueoqkWnTeE1Nvk6HQoRBm1NS6JTYmTtX5OGH3fIZGbyWv/wSmr5tm7eBOHGCAaVfeslKM30LEJk/\n3/0bRUV0WqZO5fXxcwQWLmS6aQeGzEymr1rFgfzkSep+2WW8Js6BpaAgtOwlJSI7d3KgbtrUbdzy\n8kSaNePnQID9q7CQevXsKfL886HyfoOp4eOPRVq2FDlwwEozA4CTrCyRFi0s22DydQ4GhkOHRK68\n0qrfRYtCy+O8FoMGiTz+eGia05jbOXCA1/XTT0154tCYm6jqc+as9NYiyI4dIpdcYo3Qw4e7ZQIB\ndo5mzUS6dBF58016RLVre+d56pTI00/TqA8eTC/SXPz69Rkh3i47YQIrrFs3Ro3Pz6dhycriOXav\nVITRuG+6SaR7dxqRQICeToMGIitXeut77BgbWPPmllfuzNfOe++xrC++GNpJvCgtFbnzTpEePSyv\n78cfeY7di7GzYAGNUEZGeGNlOHCAHnlKSuiA66fv8uVspKNG0UsM1xmLi3ktOndmp5o2jbJ2b91J\nXh4Hn/bt6dmb/K+7zi179KgVqd0eVT0cP//M+jUzOdNmvPTduJEe2fDhoYbN7gTYmTqV8u+8Y+n4\n9dfWeaWlofL5+Zzt3HEHo8dPnky5cB5+QQGvZ7t2lO3Z01tu/nwOjDk5VlpmJgdAp6779olcey3r\ntKyMbSIQENm1yyr7okWh+Q8ZItK3LwclQ14e+4EXP/3Ea2MfpEzezkHHzrhxItdfz0Ex3AAmwuvX\nujVn50auosK/LefksG8vXRrqmaenu2WLijhYvP66lTZxor8xF6FzBbDPIB6N+a+/soDLlnlfIDsr\nV4o88QQ9LecUy05ZGQ1sWho9v3DTPhF2gtmzRR58kB1owwZ2SqdXJEKjnpnJQaJRI5GaNVn+1q39\nyzJvHjt706YcaOrV828Qhq++sgzLxo3hy//DD8zXzwuwU17OQSIlhTqMHBn5nDVraERN/l7Xxcna\ntSIffWQdh9O3sDB06cvuqTupqOAAY/Tt3DlyWSoq6JUNHMhzJk0K336KitgOAJEvvoicfyDA2eCo\nUZbR9dP34EG24UsvtfRdssQ/79WrOVikpIg8+STb5TXX+MuXlrK9N29u5e8167JTXk6j2KWLe4Cw\nk5nJcvTvT2999GguF3rpeugQHQcz23nmGZHvv/dvaydPitx3H52f3bs5aK1fH34w3bWLs5eBAzkL\nO53BNxCgwbR7zV7LX4aMDC61AZYnHa4tr1lDu3DxxZH7y549tAtDhnDAmzjRWgr2Izvb5BuHxjwn\nh9PDRKW01HsNzUkgwPXtDz5wLw/FgkAgvIfqpKSEnb5XL3p+p8Nvv51Z2U6X48dDl6TCUV4u8s03\n9IwTkcJCkc8/t5aMIrFzJ73Kyy8Pb3wMZWU0iLGusyNHuJbbuzeNyqxZ/rKBgLVE07gx/7dp4y9f\nUcFlIr+1Zi+Ki7msV78+ZceOPT099u7lGvfy5ZH7ze7dHEhPt20WFtKRfOONyLKHD1v3foDIxlzE\nLI36G/MYhTSIGsnKEmRmnt4r84pS3QkE+BRNrIKQnCtE+MjuqVPW7pB+7N/Pdy527uR7GJHkAT7p\nkZUFDBrEl+ASjcJCYMkSvlDXsmVk+d+xAXi2gspqGncCmAmgBoC3ADgDvEnv3oIOHYAePbJdO+sl\nM9nZqm8yU530rU66AvGhbzhjXhnPmdcAMAc06FcDGACgrVMoJYV7WGwO95R+EqL6JjfVSd/qpCsQ\n//pWhjHvCCAPQD6AMgCZAO5xCi1cyIfwD5sQNNUE1Te5qU76ViddgfjXtzKMeRMA9q2I/hdM8yQ/\nPz+qzJ2vhyeavOob2/zjTT4afeOt7Fq34Yk3fZ1UhjGP6tGZaKcu8VbB0cqrvrHNP97ko9E33squ\ndRueeNPXSWXcAO0EYAq4Zg4A4wEEEHoTNA/AHyrhtxVFUZKZPQBaRJSKEecFfzAVwPkANsPjBqii\nKIoS/9wFYCfogY+PIKsoiqIoiqIoSmWwEEAhgC22tD8CWAvgPwA+A2B2WH4IwCbbXwUAs5V+h2Ae\nuwHMqvRSnzmx0jcbwA7bd7ZgenFFNPpeAODDYPp2APbIvclYv+H0zUby1e/5AN4Opm8G0M12TiLU\nb6x0zUZi1G3UdAVwHUIv0HfBdAB4BMCzzpMAXAsuyxg2gM+sA8BSWDdU442z0Xe37XglgOs95OKN\naPRNB40bAFwI4L8ATJTRZKzfdPjrm4z1+ziAjODnhgDsQdQSoX5jpWtc1G1lPJr4bwDOgFMtg+kA\n8C8A93uclwarIzQCR0QTY+Y9APfGtpgx42z0dQSdqrK9cqIhGn33A6gDvhVcB8ApAMVI3vr109eQ\nbPXbFjRkAHAAwGEAf0Li1O/Z6nqD7bwqr9tzFTZuG6y3QB8AcKWHzIOwjHkT8GUjwz6EefEoDolW\nX8O74DRtkls8rvHT90vQmO0H3wh+GewEyVq/fvoakq1+cwH0BQevZuDSyhVI7PqNRld7v67yuj1X\nxnwwgL+BU5Pfgx6LnRsBHAfXGZOBM9H3IXDppWvwb1DlFzNm+Ok7EFxuaAR2gKeC/xOdM9E3Get3\nIWi0cwC8BuBb8D5Q3O25HQXR6gokdt1GJBWh61B2WgFwxtl+DaE3ixoBsIduHQDgDELGnjNScXb6\nOnkYgE843bggFeH1XRf8PBc0cIYMAH8BcDmSq37D6fuAxzmJXr/O9mxYA6ANEqv/puLsdHVSZXV7\nrjzzhrbfmwRgnqMMDyB0/Xg/OF29EVyLGgTgn5VfzJgRrb41YN0Brwngbvg3sHjEqa/puDsA3Br8\nXAd8O3gHgF+QXPUbTt8fkHz1a9rzhaCeAHA7uLHeDiR2/41W10Sv27B8CKAAnJ78DE5bRoMvEe0E\nMM0h3x2csjgxjzblAXi9ksoaC2Khb21w+pYLYCvouVf5DRUfotG3FoBFYD1uA/Ck7btkrF8/fesg\nOes3FTRo2wEsR+gaciLUbyx0TaS6VRRFURRFURRFURRFURRFURRFURRFURRFURRFURRFURTFogLc\nN2MruH3p3xH5WeCrwDcXFUVRlDjhqO1zQwBfgXFqw9EdwOeVVB5FURTlDDjqOG4G4Lfg51QA3wD4\nPvjXOZi+DtzxcBOAMeDr3S+D27rmAhhWqSVWFEVRXDiNOcB9rBuCe27UCqa1BIMTAIwkY/fMhwGY\nGPxcKyiXGuuCKsqZcl5VF0BRqpjzAcwBw4VVgAYdcK+p/xlAO3DXRwCoC6AFuG+5olQ5asyV6khz\n0HAfANfO94M7+9UAcCLMeSPB9XZFiTvO1Ra4ihIvNAS3rDV7TtcFt+QFgL+CBh3g0sxFtvO+BIMW\nGAeoFbjbpaIoinKOKIf/o4ktwBuamwH8A1bMzvMArAimjwnKvwBGZ98S/K7uuSm+oiiKoiiKoiiK\noiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKoiiKUu34P1AOzvG+vTgKAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 55 } ], "metadata": {} } ] }