Skip to content

Instantly share code, notes, and snippets.

@phobson
Created December 15, 2015 02:28
Show Gist options
  • Save phobson/794918e83668502f44b0 to your computer and use it in GitHub Desktop.
Save phobson/794918e83668502f44b0 to your computer and use it in GitHub Desktop.
Fetching data from hydra
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Portland Rain\n",
"### imports, etc"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import warnings\n",
"warnings.simplefilter('ignore')\n",
"from functools import partial\n",
"\n",
"import numpy as np\n",
"import pandas\n",
"from matplotlib import pyplot\n",
"import seaborn\n",
"seaborn.set(style='ticks')\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### helpers"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def water_year(date):\n",
" fmt = '{:d}-{:d}'\n",
" if date.month >= 10:\n",
" return fmt.format(date.year, date.year + 1)\n",
" else:\n",
" return fmt.format(date.year - 1, date.year)\n",
" \n",
"def cumulative_water(df, aggcol='wateryear', valcol='precip_inches'):\n",
" cumsum = (\n",
" df.groupby(by=aggcol)\n",
" .cumsum()\n",
" .fillna(method='ffill')\n",
" )\n",
" return cumsum"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### fetch, prep data"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"url = \"http://or.water.usgs.gov/precip/arleta.rain\"\n",
"\n",
"df = (\n",
" pandas.read_table(url, skiprows=12, header=None, sep='\\s+', na_values=['-'], parse_dates=[0], usecols=[0, 1])\n",
" .rename(columns={0: 'Date', 1: 'precip_inches'})\n",
" .assign(precip_inches=lambda df: df['precip_inches'] * 0.01)\n",
" .assign(wateryear=lambda df: df['Date'].apply(water_year))\n",
" .set_index(['Date'])\n",
" .sort_index()\n",
" .assign(cum_precip_inches=lambda df: df.pipe(cumulative_water))\n",
" .loc['2003-10':]\n",
")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### plop the plot"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x8b535f8>"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfQAAAEJCAYAAABi2tVNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXeYHFeZr99JmijNSBrlHI9tSZYsJzknjDHY2LAYLmEJ\nXtJilmXvwiUtYBYWdi9p7wILGJaMCbYJi73GOEcZW7ZkyQpHWaORJueZzlV1/6iqVmvo6anuSqel\nep/HjzUzmq6j7qrznS/9vgrDMIiIiIiIiIgobyrDXkBERERERESEeyKDHhERERERcQoQGfSIiIiI\niIhTgMigR0REREREnAJEBj0iIiIiIuIUIDLoERERERERpwDVfl9ACPECMGR9eQj4IvAjQAdellLe\n5vcaIiIiIiIiTnV89dCFELUAUsqrrf/+Bvga8Ekp5RVApRDiJj/XEBERERERcTrgt4e+HmgUQjwA\nVAGfAjZKKZ+0fn4/cC3we5/XERERERERcUrjdw49BnxZSnkd8LfAz4GKnJ+PAM0+ryEiIiIiIuKU\nx28PfS+wH0BKuU8I0QdszPn5VGCw2BcVQlQDC4F2KWXGi4VGRERERESUM34b9FuBdcBtQoj5wDTg\nT0KIK6SUjwPXA48UegEhxO3AZ/P97OGHH/Z2tREREREl8PKBXjTdYP2qWWEvRVlu/Eczs/o3r13L\nzVesCHk1ZU3FRD/w26D/F/BDIcSTmFXt7wT6gO8LIWqA3cDdhV5ASnk7cHvu94QQSzEr5iMiIiJC\nJaPpfOI/nwbgD1+Nanwno7pqQnsU4RJfDbqUMg28Lc+PrvTzuhERERFB8cz242EvoayoqorkT/wi\nemcjIiIiXPCrh/aGvYSyoroy8tD9wndhmYiIiIhTgR37e3n25Q7SGZ1kWqOrP0Y8meFo10jYSysr\nIg/dPyKDHhERETEJfUNxPvWdpzGMk79fX1tFc1MtgyPJcBZWhkQ5dP+IDHpERETEJHT3xzEM2LBq\nFu+5eS0VFRXMbK6joa4GgI9/6yl2HuwLeZXlQeSh+0dk0CMiIiImYTSeAmDD6lksnjst5NWUN1EO\n3T+io1JEREREAXTdYPOODgCaGqaEvJryJ/LQ/SPy0CMiIiIm4LEXjnLH73YwEksDsG7FzJBXVP5U\nRR66b0QGPSIiIiIPB9oH+eqdLwKwaE4T121ayvxZTXn/bkVkoxxTHXnovhEZ9IiIiIhxdPSOccfv\ndgCwUczmM+/eFHmWHlEZvY++ER2VIiIiIsbxH7/eyq5D/TTWVfPRvz4vMuYeEkUz/CMy6BERERHj\n6B6IA/DlD11OU31NyKuJiHBGZNAjIiIichiNp+kfirNiYTOL5kwNezkREY6JDHpEREREDi/t6yGj\nGWyIRqH6wni1vQjviAx6RERERA5dfTEAxJIZIa/k1CGezIS9hNOCyKBHREREWIzGUvzyQQnA3JkN\nRf++EbmfeRkeS4W9hNOCqG0tIiIiwuKBZ49kvcl5Mxsd/14FwZdu67rBsZ5RNN1A03Q03UDXDfqG\nEwwOJ9AMg01r5jF7RvEHE68ZHouG1wRBZNAjIiIiLAZHTcPz6VsvpK5Wze0xlkjzhR88x44DvZP+\n3YPHhvjw/9oYwKoKMzQaeehBoOYdGxERERECsYTpnc+f5dw7D5o9hweyxnzjGbOZP7ORyqoKqior\nqaqsoKa6kqkNU7jjdztIprSQV2umIe5+ZF/YyzgtiAx6REREhMWeI/0ANNap23v+p+eOAPChN27g\n2guX5P07gyPJrNJd2Hz7nu3RaNmAiIriIiIiIjBD2W2dIwA0N9WGvJqJeXFPFwDnnzU35JU44+Dx\nIQBWL24JeSWnPpFBj4iIiAC27e0B4OyVrcrqjRuGQUYzK+lbpk5+6FCh5n4snqa5aQobVs8Oeymn\nPJOG3IUQTcBVwCpAB/YDD0kpEz6vLSIiIiIwOq3+81dfsizklUxMZ1+MdEbngjLxzsE06I11NSH0\nAZx+TGjQhRANwGeB1wPbgSNAGrgY+LoQ4jfA56WUo0EsNCIiIsIvBkeS/PDenQDMmV56m5dh+Dd8\nZGg0yXu/9BAAc1sLr1GlASij8TStLfVhL+O0oJCH/jPgDuATUko99wdCiErgBuvv3Ozf8iIiIiL8\n50Vp5qUXzGpi2YLmon8/CAP65LZj2T+/9boznP1SyDH3VFojndFpjAbcBEIhg/5XUsq8t4Nl4P9b\nCPEHf5YVERERERzt3Wag8bZb1is7KnVwxOyR/8L7L6ZB4Sr8XEZiZv95NLEuGCY06LYxF0K0AG8F\nZsCJNIiU8p8nMvgRERER5cRYPA3AtIYpIa9kYoYs+VQnxXCqcM+j+wGYW4TqXkTpOKlyvwuzKK4K\n06Db/0VEREScEozFTUEZlUPDL1jtas2N5WPQdx82+/pvunxFyCs5PXAiLDNXSnmt7yuJiIiICInn\ndnUA0NSgrkGvm2Ju19ManUcRjICS6Hc9vJdDx4fRdB1NM9B0g7F4mv1HB1k+v7msogrljBODvlUI\ncbaUcrvvq4mIiIgImN89vp94UqO1uS5rNFUkmdZobalXrke+uz/GT/5nd96f1U2p4h2vOeuk70UT\n6fzDyd27FtOodwEJzHC7IaVc7uvKIiIiIgLgaJdZEPfBN24IeSWFSaU1x5K0FQH2re1rHwTgza8U\nvOaSZVRVVVJZAVVVlVRXVZ4oMlTrHHJK4sSgv873VURERESERDpjDjBZNGeq69fy0/dMpjRmTKsr\n6neCcIaP95gHoiVzpyktmXs6MGFRnBDiBuuPV0zwX0REWWMYBj/5n13ssQp3Ik5P0hlTZqOmWl0l\nbMMwiCcz1NZUhb2UkzjQPpgNt8+aHonHhE2hO/h86/9X5fnvSn+XFRHhP/LIAHc9vI+PfuPJsJcS\nESInDLpaxjIX22jOnVm6ip3XxBJp7n3qEAArF7WwalE0fCVsCoXcvwQgpXzXRH9BCFEXabpHlCuD\no8mwlxChAOXgoW/eYVbhX33eIkd/P4gU+mfv2MyeIwOAOco1yLx9RH4KGfSfCyH+CPxSSjmS+wMh\nxFTg7cArmCTHLoSYDWyx/q4G/AhzyMvLUsrbSl96RIQ7Yol02EuIKIF4MsOh40OIJTM8UXXLGvQq\ndQ36SCzFwtlNykws03SD/e3mWNR/fMtGls0vXi43wnsK3cG3AFOA54UQzwoh7hZC/FII8SzwImbN\n4i2FXlwIUQ18B4hZ3/oa8Ekp5RVApRDiJtf/goiyYN/RAb551zZ+9/iBsJeSJZbIhL2EiCLRdIO3\n3/5HPvbNp3hia7snr9nRN0ZVZYVy7WA2um4wGksxVSEVu9FYioymc9G6eVx5rrOoQYT/FJJ+1YFv\nAd8SQqznxPjUA1LKlxy+/leAbwOfwDwAbJRS2gnL+4Frgd+XuPaIMuIXf5I8v8tUurr2gsWhK3IZ\nhsHuQ1ExXLnx5LZjJFJmVfqIJYXqBk036B+KU61w/jyWSKMbxQnK+H00GbOiW07b6HKJutD9w5GK\ngmXAnRpxAIQQ7wS6pZQPCiE+aX07NyIwAkRxmtOEvsETpRZ6yMIS9z11kO/8dkf264Y6dcVEIk6m\nuz924gsPrFbfUBzdgGXzp7l/MZ8Ytg4uKnnoMUsqt6He+bNTETWiT8rmHR1s3dvN2151ZlEHOBs/\nd7J3AboQ4lpgPfATYFbOz6cCg5O9iBDidsy57BFliq4bHDw+FPYyyGg6//enW9i+rwcwK4Y7+2Js\nFGrkJSMm56S6Bw/Ohc9sPw7A8hJGpubFMPDaP7aLN0uRT/VLlW2/JSaj8jCbcmPL7i6++KPnADhr\n2Uyu3Liw6NfwrQpESnmFlPIqKeVVwDbgr4H7hRCXW3/lemDSfiEp5e1Syorc/4Blfq07wnsOjTPm\nYTnox3pG2byjg2Ra46J18/ji314KQLXC1c0RJ+N13YM9NrWUzTMXPwu8hyyD3tykjvHMjkWNDLon\nHO0a4UuWMYfSD2KOPHQhxM3A1UAGuF9K+WBJV4OPAN8TQtQAu4G7S3ydiDIhndGzvaphY2+Mt1yz\nmrdcdwY9A/GQVxRRLENj3rYa2uHsBbOaPH1dLzAMg9u//ywv7ukGKE6FzecWsnjSPFgtnaduqqJc\n6Owb48Nfe4yU1W3hhkkNuhDiK8BFwC8xPfrPCyHOk1J+yelFpJRX53x5ZbGLjFAXwzDYKnvoH07Q\nP5zgSMcwmn7idPnSvh5GrVnT1VWVZDT3N22xvLSvh/buUQ4eMyMFkTxl+ZJbi+EFz+/qBKBJwbGp\niZTGi3u6qa+tYvmCFs4poWXNr2iYXZgY1Z+4Z8/hflIZnUvXz2deayN3Pbyv5M/NyadxI7BGSpkB\nEEJ8F9iKJTwTcfqSzmh87vvP8tK+3oJ/b8a0Ot5w9Sp2HOhl846OQKctbd7Rkc1L2cyZoY7aVoRz\nnn25A9k2kP3a7V1kGAYZzXyVKgV70BMp0wveeMYcPv728yf528GSsDx0lafTlQtdA2ah5zXnL6Z/\n2N2B1cmn0Q20APauXZPz54jTmJ/evydrzN91w1lMa6xl4Zwm5s1sPOnvTW2YQmVlBS8fDP62efSF\nowC8+uKlrF3RSmNdDRtWz5rktyJUZMcB8/5pqKv2JJd+uGMYgHMUvR8SSdMLrlfQaNoeet0Uddv9\nVGV/+2A2/dfVH+P5nWY77+zp9YEY9H7gJSHEf2Pm0K8HuoUQPwCQUt7qagURZYtdIXz7ezZx7hlz\nQl7NyXT0jvG93+/I9r6/7fozlWr7iSieuGXE33jNan503y7Xr9fZZ3pGqqqc2R56XW3xRtPvBrFU\n2jToNaUMizmNG9GPdo3wD19//C++XzelitnTG7JSuqW+SU4M+m+s/2y2lHSliFOK4z2jdPXHaG2p\nV86YJ5IZvvDDP9PWaSoWX7dpSWFjfhpvMOWE7ZXbeVu3mZseK9S5arGaQ0UGR0wvrqlevYOoLZc7\npYgOkdNR6n1gOMG//3Jrtt3SrifatHYuqxdPp6qyknmtjSybP4262mrXB7EJDboQYq6UshN4NN/P\npZRtLq8dUca0WzOQz1o6I+SVnMzetgE+9s0ns7nRn9x+HdOn5p8hfTpuMOXA0a4RRmMn6+x3D8R4\n2ooI1ZegTpaPbqvLYfZ072oqvDwbHu02D6SL5qhXgZ+yZsirPNBGBZ7f3cWLspvKygpsZeGpDTW8\n5bozfIkMFfLQvw/cAOTGB2zVBANY7vlqIsqGviEz17PxDOeVt0EoRb18oJeMZrBm+UzeeM3qCY15\nhJocaB/kw3lCkjarFrV4lrftGTQ9dC/mePtxb9sRpkVzphb9u34fVlNpjZrqymjC2gQMjCT40b27\nePlgHwBf+sAlnLVspu/XLaTlfoP1/2VCiBopZdrqH6+VUo76vrIIpWm3vIcFs9XyHl6wenZvvXEN\nqxdPD3k1EcXwyJY2vv6LrYBpxDatnXvSz+trq3nt5SvYKrut77jzh7sH4tRUV9LcqGYb49a9PVRX\nVbJwdvEG3W96B+ORd16Ax19s55EtZkHuojlTWV6kN+5b25oQ4hbgM8A6YDHwmBDig1LKaKjKaUy/\n5aGXEq70s2vNVrBauVDNvGjExPzgDzsB2LBqFp945/k0TBBa98In1HWDjp5R5sxoUHLK2pNbj9Hd\nH2PVohZXhtOPZ03XDfqHk0wvQYr2dMEuuPzahy9n1SLnjoXbgIeTO+XTmLPMkVIeAM4FPufushHl\nTkkDI3zcN0djKf7znpc43jum7CYdMTGdfWMMjaZYuaiFz7//4gmNuVe0dY0wlsgoG8XZfcScBPiq\ni5aGu5BxpNIaP7zXPHipGDlQgXse2cd9T5vqmF7WZzjBSZX7FClll/2FlLJbCBHtlqc5w2MpGuqq\nlQm7Pb39OPc/cxiAM5epVagXMTlPv2QWvBUj+uPG+9y21xzQs25Fa+kv4iP2aFiVNBOe29XJ5//r\nz9mv3//6dSGuRk2OdA5nWypvuHRZyaqUfirFPSWE+AXwc+vrNwKbS7tcxKnAjv29HO4YZmGJ+XPD\n4z6xwx3DfPMuc7rv595zEeeI4jbBkKe5RnCineeGSyafu+RFIdbRLrMGRCxR00O3U0dup5l5+ax9\n5zfbAfPQdcs1q1k8tzQdd6+ff5Xo6B0D4G3Xn8GbXiFKeAV397YTg34b8HfA+4A08ATwn66uGlG2\n7D86yCe//TQAc8cpwk2GH2GdZFrj77/2GADzWhvZsHqW4w0/KtBVB9uABaWzn7SUzupr1VNhA/P9\nqK6qpFYRJbb27hF6BuKsXNTC1/7+8pIOVaf645bOaHzrbtOxmNXivnOiFCa9m6WUSSHE3ZjT0R4A\nFkkpU76vLEI5EskM//DvZkvRuhWtfOhNG0JeEXT2jqHrBs1NU/j6h6+IcudlStYjbQxGRMXuo55S\nitJZAbyK9ozE0kxrrFGmLWz7flN297L185VZk2oc7hjOigGdvdJtqqS0G2nSBKgQ4k3AH4D/B8wA\nNgsh3lbS1SLKlrF4mjd96r7s1x97+3lK9HjbodpXXriERgUnZkU4o28wQWVlRVFTz9wYT9tD98wD\n9tDGZTSd/uGEq2iF10a3d9AU4VmlaBGhChzvMcPt7715Ha0leuhBVLl/DLgYGJFSdgPnAJ9wd9mI\ncuPp7cfRDXMD/M7Hr3EXGvUwhdbVbz5EKspjRjhnYCTBjGl1zqaeeWCrkpYWeTHSpUFx+PgwyZTG\nGUvcF3d6ETEwDIM/7+ykuqqCJSXmzU8Hjlv58/mziktFeomTu1mTUo7YX0gpO4Dgh1pHhIp9Qv/0\nrReyYFZpxXBeew3HekazQiTTp50aPbHb9nbzyJajaCHMjQ+TZFoLdHJXIpVhSk2VkuHjFy3hHFVE\nmzr7YrR1jnD+WXMDS4mUI+1WoeX8Vvefm59V7juFEB8EaoQQG4APANtKu1xEuWLnhloCKlpywr1P\nHgRg4ewmLlo7r+TXUaXq9njPKJ/+rtlAMnNaHesValnym0RKY1ZLsQa99M9taCRJS5Oaxske+Vvq\nwRm8LUAbjZv1DcW0FJ6O7DrcT2NdNbNdSAm7PV86rXL/JyAO/AB4BPhHd5eNKDfs2dEtHqhDeWU+\ne6yowSfecT51ilYrO6F7IMZ/3v1SNgoCJ0LCpwO6bpBMaY4/Qy+M1Vgiw9yZahqowZEkdVOqOLeI\nOQl+Ys9lr1NwLrsqHGgfpHcwzkXr5jlLG/mEk0/om1LKdxHlzU9bkmmN3Yf7mTez0VXIzevg5vBY\nisrKipKGV6hCRtP5my88GPYyQsU+yGQywaQZDMMgkcr4rkZXCsm0xmg8XVT7pd/Erbns9SXMZc/H\nqaj78PxuU3ttk4tIoRc4OUqsFUKokcyJCIVRq6Vo1aIWZTYZwzDoHojR0jRFmTWVwgu7syKM/PqL\nr+HWG9eEuJpwGLbur1lFymSWahgSKQ3DULMHvW/IPNzMbA6/g8Qmbs2hr3XroZfxczoZtlCRV6Nu\nSz3zOPmEdKBNCCExw+4ASCmvLvGaEWXGiDWbuqnBG4/G8OCIfqxnlL6hBJesn+/BisJj92FTs/tv\nXrtGSQMTBImkaTDmtzqrDnZ7gIsnbY/Tj/fb3b2998gAAHPcaoB7aDtlm7kmL9Jtpyqbd3QAsHSe\n2xnn/ivF/R9XV4goe4bHzIK4ogax5MPDTeb5XaZnu2FVeRWOHWgfpH84QSqjc7RrhHse3Q/A1ect\nDnll4ZGwesKDqoOIJcwDakOdd9fz6tYesIpPlxY5bnMi3B6ef/HAHv5gFZ+6KfY6lTnSMUw6ozNr\nen3osy2cKMU9LoS4GbgayAD3SylP76TfacbOA32Adx66FxzrGQVQdlpWLu3dI7wou3lxT3d2Xnsu\nKxe1/EVtghdRjHIh6zEH1Lbmr4fujrRVR1DrsYJdqdz3jDk17B2vOSsaSTwBOw+Z++ONly737DX9\nnIf+FeAi4JeYOffPCyHOk1J+qbRLRpQbOyyDvna5OpOpbIU4T2Yy+2g70xmdj3/rKYZGT6glX3/x\nUua3NtJQV8OqRS0nFfWdwmnGCcka2CI95lI/NpUNui1JW1MTvuCNpukMjaZYs3wmb7h6VdjLURZb\ntnjJPPeiO0G0rd0IrJFSZgCEEN8FtgKRQT8NMAyD/e0DLJrTxMpF6pzQxyyD7kbuNYhiuj2H+xka\nTTFjWh3/8OZzWL14upLV1WFiG9ig2qLsIi8vQ+5eYVf6uw3denFvx61USGN0vxZEpfvJyQq6gRag\n1/q6JufPEac4nX0x4kmt5FGJuVR4lGnUdINte3uoqa70fLiG1/RaVcv/65WCDavV6CtWjYHhBBBc\n0VVMYQ/dDrlPqfaoRczF79qGSsX3SSX8uZ9K++ScrKAfeEkI8d+YOfTrgW4hxA8ApJS3lnTliLLg\nF3/aA6CUhvMzLx0H3HnnQbHDmlKlqiqZChyyRIuK1RMotc7gRauOQUX9/5RHHrpbEqkMX/zxc4AP\ntTMKlYcMjCSYUl3lai+xVTTdzq4H98WVTgz6b6z/bLa4vOYpx31PH+L5XZ2sXzWL1125MuzleMr+\n9iEArr3Auypst/VeP7l/FwBvv/5MD1bjH7FEmke2mDKeK6KCogk51j3KjGm1jrso3EaT+4bMiIBY\nol5BZdrOoYds0Hfs72X/0UEALjnbm9ZQ1epDjvWM8v5/fZil86bxjY9cVfLr9A3Fqa6qVKKtz0mV\n+4+DWEi5ks5o3PG7Hei6wZ4jA6eUQX9uVydHu0aYO7Oh5HGAuXjxQCeSGTr7YlRWwFXnLXL/gj6y\nv30QTTc4Z/UsZhcrmuLTmlRkaDTpyUALp/QNxWlpqvXknh6P28NqOu1RDt3dMugeMFNFf/+mDaxb\nqU4xrJc8sfUYcELWulT2tg0ys7lOCYGr8Espy5x4UkPXzaf4VGs1stvVXnvZipBXcgI7X3Xx2fOp\nDlEz2QkdvTEALtuwoIjfCn9TCJJkWiOR0pgWYEqifzjh+XQ+rzbztGYbdI9qQ0rckjr7zFGg8wI8\naAWNrcrnZsqf/T55rd/u57S1iAIkLJ3jU5GdB/uorKzw3BN2c+6x328vK6L9OobZvfLzXUzNOtUZ\nttr5mhudG1g3xZWxRJpESmNms5oiKSeK4sI7rGY0nSe2ttNQV82KBd4I3ATFwWND2Zx2ITRd54Fn\njwDuZHYPHDNTkisXevQ+uTwYOtoVhRBLgTXAH4HFUspDrq56CmHLVp5qDI0mkW0DrFk+kyavis88\ncGKStqqYByIkfvrCfUNxfvuYqQK31IP+1FOV49ahpzkgD91ukVOhxWg8um6w+3C/J90bbuzCi3u6\n6R9Ocv3FS8tqiuHTLx3nX3/yfNG/56Y63W4zXK+IYqUTYZk3YY5PbcAUmNkshPiIlPJnfi+uHDje\nOxb2Enyh32olWjZfLWNkG/TagFTFSuHnf9zDLx+UgGnMy6EaPyzstr5SNtVSIj2aZv6Siumabft6\nGIunmT29nspKb46bpaQB91nFcGuXz/RkDUFx4Ji57ovWzWOVA82M6qpKfvCHna4idLr1/lZVehxy\nL/H3nDxFHwMuBp6QUnYLIc4BHgJOa4OeSGb41UN7ufuRfdnvnUopdFv9yLV+u8fsOWIOMylWVSwo\n9h8dzBrzGy5dxjtec1bIK1Kb4THzPitKtMiFrdN0ewNWr1ZhcMQ8RF+xcWGo6+iwnJRyk3qNWX3z\nb36lYJlDLXy7LbdU7AOiV+fDINrWNCnliBACACllhxDC0eBiIUQl8D1AYE5tez+QBH5kff2ylPK2\nEtYdOlv39mSN+dkrWxmNpenoO3W8dVuq1A+Dbrg4E7d3myHaOTOcTeYKmmd3mlOXbrhkGe973dkl\nv86pdDgsRNAHx4xVdOZ1EZMX2ENqwkzRPPx8G49vbaepvoY5M1xOfJsAN8//hK9pGNn0TZARMdtD\nr/TYQy8VJ6vYKYT4IFAjhNgghLgD2Obw9W8EDCnlpcCngS8CXwM+KaW8AqgUQtxUysLDJp40pUff\nc9Na/uVvLznlipN7rLYVL1t73L5FPQPxbCHLBWfNcb8gH7AlaV+5aUlJv69A58tfcLhjmCOd7lp7\nJsL20McPp/ELXWEPPZG060O8iD4V/++LJzP8x6/Nrf26TUs8P/T4+Y4/+3InW/f2UF1VQXNTcR0M\nbg7PdsTHqxRJlhIX5eQTuw1YgDkL/QfAMPABJy8upfw98F7ryyXAALBRSvmk9b37gVcUs2BVsHO5\nJ4sJnDpulV2h7dcpvRS27O4EYPaMBk/10L1sN8xqzJ8i+tc/vX83f/eVR/nglx/15fWzKlslGPRS\nPD2/Q+5u7qSk1cHhZX1IMev5wg/+jK4b3HjZct55wxrP1hAE9n71xleIIifVubsPsgdEj07iQQxn\neQ/w71LKT5RyASmlLoT4EXAzcAtwbc6PR4Dy6ouwSKat4izr5lHRsyoVTTf4059NT9iXGcgl7nrP\n7jQN+j++ZaM36/D8UG1kN5apAXmcftI9EOPXD+319RqDo0mqKiuKCrm7+dh8C7l7cC/FPezgKIX2\n7hEAXnuZd2NAgyCj6fz3EwcAOO/MYOclaLp5P1VWqWEAnBj0BcCzQgiJWQj3GyllrJiLSCnfKYSY\nDTwP5FqIqcBgod8VQtwOfLaY6wWB7aHntpecKnlP28usm1JFk4e5TTfiG8d6RnlxTzctTbWcuXSG\nZ2vykt8+doC9bYOcuXRG2Q+0iCczfOa7m7Nfz53pT6RmcCRJc1Ot9yHLCVA65O6DxkJx19dYNn8a\nc2eqWZ8yEZt3dDAwYh4MF5Si+eBi37bsOZUee3SlLmnSY6qU8qNSymXAvwCbgG1CiJ86eXEhxNuE\nEB+3vkwAGrBFCHGF9b3rgSfz/vKJ698upazI/Q9Y5uT6fjJ+5KN620PpjFqFSsUpnPmLXXm7dN40\nJSQW83Hf0wcB+Mhbzw15Je7Zsb83G23wi8dfbKerP8asUus0Stj1sh66ggZ9NGYepL3o4Mg+Ig7f\nI8MwSKS00A4TbugdNOt93ve6dUWn4txuJdkDokceutuJlI7iTkKICsyxqVMwq9Mnl+Ix+Q1wjhDi\nccx8+Yd8V3GCAAAgAElEQVQwc/KfE0I8bb3m3cUuWgWyxTwKT9E60jGc3cCKYdSDWeOFKOX0aXsv\nF66d6+1iPKKjd4zugTjnnTmH2Z7UHYQb7hlLmPfArTeuodWFktZEHDo+xFd+/gIANxYZ4nWzCWdz\n6ApWuR86PkRDXTWtIajYpTM6um4ore8wEYeOm2ptpY54dlN1nw25K+JkOBGW+QZm/nsr8HPgQ1LK\nhJMXt0Lzb8rzoyuLWKOSxPIMtVcp4v7inm4++73NXH3eIv7hzc5zzppu8OBzbQBFV4v6yYkKYDU3\nnDt+twNwX32vxrZwIgI1Y5ppzL1OJ333t+b7deuNawLtu1a5D713MM78WU2BpR9ysVvmyjFVtLdt\nAChNwtXtO31CWEaN+8nJp7cXszK9x+/FlBPx8UPtFTmh2ew6bA5WeWTLUccG/XjvKP/nG09me9DX\nLFNHKSpbAVyj3oZjGAYvH+ilosIcGlPO7D86yMNb2tjXZpa21NdV+3Jv20Mtbri09OxZKWcMO50U\n9njS8ei6GfL2SpK22E8s4UOF/UR4fTi0zmih5P51zZ+2Nc+Hswgh3iulvAOYAfytLSxjI6X859Iu\neWoQT2aorKDIFongKKX456ltx7PG/KNvO5czl6lTfDZorcuPNIDb/WVgJEkipXHRunlKRTVK4ecP\n7GHL7i4AqqtKLDJyQCyRZvn85pKmirnJM97zqKmvf7ZiI0Ftg+q1h+w0nDy+JsgXfPJ5RmNpFs0p\n/T511YdueGzQfWxbq5jgz6BWdDlQDMPgaNcIuw/301hfc1KBlkpV7noJggfdA2bzwjc/chVLfFSr\nKqXv2y58mdvqXbW12wIUmzZLdGXxnKmevF6Y9A7Gqa+t5t//9xVMa6z1bjBPDppuEE9qNNQHG22R\nR/o5eGyIxXOnssInWdNSNQ2yEb+QitL6h8ws6oyp5XUgNQyDsUSahbNLNOguok+xRJr7njLnlCkf\ncpdSftf642Ep5Y9zfyaEKEu5Vi94+UAfn/z20wDMyJmprMbHeYJSFIx6LKM5y4/ec9xFbgdHzTrM\nFgU94COdZv/ukhKLcvIRxOFwaDRJ72CcnQf7spK6HX1jzGqpZ76Pc7Bt4+VWfKfY90haudZrL1js\n6rr5cPv8hz10yD4we6kMGQTP7exE1w2aGoIXcvroN57MFhB7Hu4vcQMoFHL/MDANeL8QIlfHshp4\nK/Ctkq5Y5tgSmBetm/eX1bkKuehZjeEirOhYPE1VZYWShTGDo0mmVFcqubYjHZaHPs8DDz2gk2Es\nkebd//Jgthgql3yDLby8s+99ymzvKzlfXOJ7FLcKWUuthvaT8UJVriny9FyuBv3JbceBYFMoPQNx\nntx2jDbrIH/7ezZli0fd4qdS3H7gXMzHJ/cySeCd7i6rPoZh8OuH9nLg2BDpjE46o9HdH88OYHnd\nFStPyjErVhNXUsg9nsxQX1utZJ/30GiS5qm1Sq6trWuE6qoKX71arxkcNfP+S+ZOZeMZc1g+f1o2\nDD2v9WRvw8u3vLNvjJ//0Zxwdd6ZwerxjyXsyIB6h8K0NVe7xuOaHCc+xo79vdz5J3NCYLkZdHvM\n842XlqZuV8qt/Z3fbOe5XaZq5V9dtZJzz1BnrkShkPu9wL1CiF9LKXfn/kwIUV6fegkc7x3jZ388\nebRe3ZQqFsxqZO7MRlYuUlOxNpZI89Dzbdxr5XaK8dDjyYyvY0lLzVkbhsHgSNLXvH6pGIZBW+cw\n82c1KVc5XQg7xLtuZSu33hicbveew+b424vWzePyc4IdExqzeuu9nAPgFbaHPqUm+HvoP369FYAF\ns5rKTiVuJJaisa7ala5AsXUP9oTAf73tUs7wSbXSz3noZwkhfgk0YR5oqoAGYFaJ1ywL7CKR6y9a\nyq03rqGmunLSmybsgPuew/189BsnC+8Vkw+PJzLM8EFExC3xZIZ0Rlcyf94zECee1DzNnweBt5O9\nnLN9fy9gPleuKXIjzqcdoQrptOmhB901k0xrdPbFmDGtju98/JpAr+0Fo7GUK3nqUqJPybRGfW0V\na5Z739YbhFLc/wU+DOzGzJ3/EPiVq6uWAfYUqMVzp1JXO/kJ0KuK6WIZiaU42jXC4Y5h/mwNL5lS\nU8XNV6wA4HyHQieGYWRD7n5TbKmB7wVxLk5ibV1mHm3x3PKpcNc0nf/ZbEZwghTqMQyDrbIbgFWL\nSq8yLzUFYHvoKk7Csz30Utr48uH0LeqyUohBpj+8dHxG42mmBlwQl85onn1OXuNk9x6QUj4qhLgE\naJZS3i6EeMHvhYXNwKjpoU+fqp7HCvD8rk7+857t2WKWXL70gUtIpTV+9/gBx09POqOj6YavBr3U\njXhoxAxxNXsss+tFbtieULXI45Y1P6M9P71/N4+90A7ArOkO2wA9KPg82jVC71CCyzYs8HToj1Ni\niQyVlRVKypumM3ZRXHAh92Ra466H9wH+Dd/JxWunJ53RSKQ0murd3UvF3tmptH7SUC6VcLJ7x4UQ\nqzE99CuFEI9QpiNPi8H20FuK6MsMssj9qZeO0zsYp7lpCuesnk1DXTWVlRXMmFbHyoUt7LZylU75\nC+U7hRiJW7r5Co4kte+TmV5VufoY6UmkMnztzhd5+YAZ9v70rRc6iuB4taL97ab63NoV3oQqi33c\nxhJpGvwu+ixxD0h57KE7YfP24zz2onmwW7tcLaEdJ9jzNNy1rBV/L6TSmu/7pOdKcTn8E/AF4K+B\njwPvA/6rtMuVB4c7hrMnV8cGPcCI+wPPHuGRLUcB+ME/vdKT06LKBt0e56piMZPKednxyMMDbN7R\nQUUFbBSzuWBNsINuRqxpYm5TJ6UcenoH47R1jjDbN40FdxtA0uMcupPlDFkG8bY3rFdKFdIpA8Pm\nYdqrljGnpDI6zU0+Hbx8bFsDQEr5OPC49eX5QojpUsoBd5dVG3nE/OfNndnAvKKqPoNx0Z/cZp6q\nb75ihWehH9ugNwSRQy/yfYr5PP3NDWMKV06Pp2fQVAK87Q0buG7Tkkn+tvfY1cFBf45bZTdfvdPM\nEp7lQyGTF9ga840e54MLeXr2YbS4PU4d+kcsdTuXBr1Ybzid1kLpRnBCIWGZR5nAQgkhkFJe7duq\nQuawNY7vfa8723Efd1AOeiyR5qV9vUyfWsvfvHath69reegKepq2EpvXOXQv2Hd0kMrKCqYqmA4Y\nzy+sXuOS54+75EC7+Vx5pUTodCP+919uZWg0xbzWRt538zpPru01dz9iRgSD9DbtIkEVn3kn2J1I\n0128Z8UGVgzDIJXRfU+NlDrStdAneXtJr1jm3L/5MPc+fYiKimAKRYrlWI8p0enUy3F6W6gccj9w\nzMy9nrHEr57P0h6eodEkHb1jLJrTpOyQnlxsOeBScthexJ5GxlLeCPAUsQnrusHgaJLFc6fyzY9c\npaQwUSqtZRX7FswKzlu2D/EqRr6cMDBse+huu1+c3922AJBfz7vbu7OQsMzjAEKIy11eo2zIaDp3\nPmCKyXziHeezcHZxlctBFMXZD+HlGxZ4+rqjQeapi3yf4skMUxumKFdZ+uuH9wKwbJ4PNaIe30uH\nO4bpG0pw5tIZxb+PFRWe3NzpjB74wad3KI6uG8yb2aikMQd4+Pk2ADatnetd25qDf2tWaEfBQ7wT\nBux2VhedSMXeEilb0U9RESknn+Tncv5cA5wNPAk84cuKQuK5nZ08+NwRBkeSXLp+PhetK26udVCb\nRd+Q2abW4PGpuqvf7Ef1q2gISn+P4klNubCgpun86dkjALzhmlWeva5ft9HzllRlmP3yqYzmubTp\nZHT1xQK9XikMWJ0Sm9bO8/y1C0WfYsmQ0mweHFbH4mnuf+YwEOzAplRW0c/n+9ivKncp5VW5Xwsh\nlgFfL+1yapLRdL72ixcZi6epqa7kLdedEfaSJuSpl8xhBM0e52ztqv7ZM9RLM8STmdDyvvnYfaif\nz31/M4mUxprlM/MOM1GNrn7TsP3FQKEAMXOP3nk2TlIltijROauDEbYsZR8esarNly8I9j4ajaWo\nqqwILGri5WHV3gfBfW1NMcGnE+2F/njobt+jolclpTwEqGvxSmD34X7G4mmuvWAxP/vcq0oWCQmi\nxt3Wwt60ztlp3olOsa4b2dzQoiLTDKVQzPvU0TvGWDxdlB6A39zx+x2MJTJsWDWLD71xQ9jLcUTP\ngBnZme1USGYcXqST0mmNKR5shMXseXZlfRhCNk6x1QaD1lnoH0ows6Ve2VREIQasCvd/fMtGV+sv\ntgUy2w2kWMTQZtJVCSF+yIk9uAI4E3jZz0UFzc6DfQCce8YcpduPntvVyUgszYqFzZNqcBdzjw+N\nJdF1g4vWzStqOluxlPLcySPmAcaNVKiXjMRS7D86yIqFzXz+/ReHvRzH9A8naKyrLqno0as7IpXR\naQ5YMtMeQqNisaeNnwa90EEskdKY2azufleI5yyZa68VGifD1sTwWz7Yz+Esj427zl3AQyVeT0n2\nHzWrqMWS6e5eyOequJ/fbxbsXXXuIk9ft6PXzJ/PUTDc3mHlQFVRsrLVqVYsUOOA4ZSgdPoLkc7o\n3vbvOnjc7OpxlbsQNE2nuqoicH3wZFpTUgbXCban7MUExmK27THfNTHcHZ+d5NB/LISYCuRau7lA\nm6srK8LOg328sKeLKks2tVT8jlrFEmkOHh+ivraaV1+81NPXth8O1aRVdd3gme1mrmy+j+08pTzQ\nfofcSm2lm4hkSgu1PckwDM+GWhQTYk2mzHtbVcOVzmiMxNKcvTLYA6tmpdmCnrbnBQPDCdq7R1m7\nYibVLsamAkXbz7GsMqSakQ0nIfcvA+8F+qxvVWCejcOrrvGI+546yHd+uwMw1ZL8DDe7ZdveHgDE\n4umen+SD9mKc5PU37zjOv/1kC5puMKWmypeiuFJyb37P1PbrDkykMsxsCW/QUEYzMIzg233sKWaq\nGvRDx4eB4A/T9kFHtVZQJxzuMN+z5kav6mqcH56zE/vq1TwIOVnVzcACKeWo34sJGtlmSry++uKl\n3HT5Ctev51fAXdN0/ucZc9zlqy9Z5vnrn/Bi1LhJY4k0X/7ZC2i6wRlLpnPdpqWTjq8NCjsFMCtE\n41gshmGQTGvuvDGX6aTuAfN90zTvnhInr5RUPORuFyvObPbjwDrxz1QWkpqMHmvC5Hlnzg782oHl\n0H0czrIdqAVOOYPeO2hWSr77pnWuPQc/K0U3v9zBS/vMCVlL5jkrAimmetP20IOcjV2I+54+RDqj\ns3xBM1/+kFq6RvZ0tVaF2ugmI5nWMIzSvVQvbm37fQtafVF1D73X0pU4K+DhKKPWoBx3k8pKw206\nKasQ58EhqNhbe8xndT23z5oTg/5TYL8QYgeQsb9ZzlruhmHQ0TvGywd7aZlaq6zqj0231UP8ples\ndi+bmYdkwAZ9stPnywfM7M4/v/eiAFZTHENWX3NzgGIWbonbOv0hRmDaOs0wadDCNnZ4Nqjok5N0\nUi62lHPQ+g+2MmRTGcq+2oc0r/arYj6yh54zxaTKtm0NU0Tm74EjPq8lMO58QPLLB81BFQtmeWcg\n/SpyHxw1K6v9Gnd5wkP39yZ1EsUwDIMDxwaZ2VynpNEsR4P+0n4zuuPVUJRieWJrO/95z3bA6+dt\n8gfO/iuNPm/ApXpWh48PU1lZwRIfDjqFluR3LYif2JoZUzwpsHT+d//1J88zEktTN6WqZD0H53g/\nnMVmSEr5k5JeXVF+89h+KirgdVes5LWXq1/b57cRsSMAXo9uLIWxRIah0RTnnTkn7KXkZcg6XPld\nxOTl4bDPyjmuC7iS2ma7daC45vxFnnyuxWzCyVSGVkXFUwzDoK1rhPmtjb61rE106ElZ89dVHQNa\niKxaW8Brf2F3FwDf/OjVvhUT+jacJYenhBD3APcDKfub5Wrkj/eMkkprnL2ylXfduCbs5TjClq8s\nRe7ViWHotHTcVZAwtbXqVc1RD44mmdpQ475dZiJ8MDy2ZrebQh435ws7X/uO15wVuGFNpDTl2jFt\negbjjMWDb1mDnNoCRYsFC5H2eECK03s7oxmsXtyipF6HjROD3ggMA5eM+35ZGvStVvuX18pjfu5T\nQ6NJ6qZUUVdERWox60mkTIGJKp/b9py8esKWVgyg+raUz2xoNFlW4fZkWssOkSm1orlYeczxZOVX\n64M3rMNjKV8HDrlh1yFTBTHogjg44eWWs0H3IuRejE+s6zpVlcFEBXyrcpdSvqu0l1aTHqt9RtWQ\nbj6GRvw1IsmUpsyDne2JV7AqWdMNRmKpwOUm3fCje3dmIzzTXc+NLp4/bj7M9v29LJjVGHjxqS0b\nPFVRHXc71VXsmGbHVFRMaBgCmxrmA6mMvwNS8qHrBrqBf5E5C7cRLCfCMofIE5WQUqqffM6DLd05\n3YUqXCEMw/A0rGgYBoOjKVb4OIlJJQlIrytYvaR/KIFhBDuu0S2Pv3gMgH961wVMdzE3uhQMw+C7\nlnDTNecv9ux1nUQMNE3nM3dsBuC6TUs9u7ZXGIbBI1uOAuEUKybL2KDbXRvFRCwL4sAb1nQzKlBV\npV4tRi5O3pErc/5cA7wOsy+9LLENutfjR92GJScilsiQ0XSfPfRMoHnGQtXJdshdFZGbXHoGTY9q\nXqt/MrReks7ojMRSzJpez4U+zNqeDPveXb24hVuuWR3otQ8cGyJmbfznnaVeNO6+pw9xrGeUiopw\nZijY743fAin5cFvwOTCSpL62ypOoolPfK2MJIvntobvFSch9fLval4UQW4Av+LMkf2nvHqG6qqJs\n2jWyMoclzvx10tqTTGnUtgRwUnfw8NzzyH4ApipQcT+eWCLA0YkeVLkfOj4EeFMvUsombOfO/Qop\nF1qT3Zb1V1etVCadlMuRTnPC2ntuWheKnnqXFe4PUljGi8hlIpnhcMewp4dqJ0I3mmZ56AHJg/uW\nQxdC5Ep1VQBrAEcxIiFENfADYCkwBfgXYBfwI0AHXpZS3lbUil1woH2QYz1jLJ471TfddsPwtkDu\nmR3mcJKib2CHa7BlQVXxiLtCqHFwKgZit634OeTEy3vHlgueN9Pl5lfimuxomNfRHyfvkV2LEVQB\nY7ERur1HTNnpyzYs8GM5QOGPbY9VX1BOBZ4A3/6NqWfgZpBWLk4/tcA89ACU4j6X82cD6AXe4fD1\n3wb0SinfLoRoAV4CtgGflFI+KYT4thDiJinl74tadYkcOGZ6LBet8z786FeVu50v2uRTyDSV0V3J\ngpbCRObTMAzG4mnWLJ8ZSASlmI/soeeOcO/TpoH0ajPxk2Ray+ZowyoAPdBujiUOoygtK5akqFZ5\nu6UQF0YkKpnWGBlLUVGh3oTFybBnAtz2hvWBXveUyaFLKa8SQsyWUnYLIRqA+VLK/Q5f/9eY89MB\nqjClYzdKKZ+0vnc/cC0QiEG3Q6Z+zrL2WiwubYV6/DK4QQ6vmOxRONGOol6eyp6KdfMVKzj/LH8U\n+7yko3cMw4AL18xl7Yrg+5wNw8iqw83wqbq+UGTFHjikYnFlOqORSmusWNgcytChfW0DJFIar/Fh\n0JPfxOIZ6murPO00cRKgCz6HXpolmXR1Qoi/A/5ofTkL+IMQ4r1OXlxKGZNSjlnz1O8CPsXJ+/oI\nEJiaSTzAHmevSKe97Ln8S5IKtYmp3EpjK8S99rIVgeXR3PDyAVOdLawWuwFrGAvApes9Dis7ePvj\nSXW7JexUhB9zGXKZKGp4xKrLWb3YP8fGL0biKW+jdw5Dq2mrVc73tjWXv+9kde8DLoNsgdy5wN85\nvYAQYhHwCPBjKeUvMXPnNlOBwUl+/3YhhJH7H3DI6fVzyY4MVFRYPx9eqyKNJ5m2qsqDNKIT9cZa\n/1YVi5hOyO+WR4jSnm62dsVMj16xOI+ho9dUH/yrq1aGEvY+4aGr96z71WmTj3xRDPuwNct3PXJv\n2ba3m56BOHPd1oSUQMznKWte4eRurwGSOV+ncPh0CyHmAA8At0kpH7W+vVUIcbmU8gngekxjPyFS\nytuB28e97lJKMOp25auvM4ANA/fnrBPYIgrFai47XUFQg1lg8irXsDSandA1EGNqQ01g0QO3Iyaf\neskspvRiiEQpd3Nnn2nQw9h8Idj7uljsw2FY+WvVxiU7ZZul8nnlxoWevq6TJ61vyBzZ2lgf1NS+\n0n7Pyep+BzwihPi19fXrcZ7z/gTQAnxaCPEZzPfu74FvCCFqgN3A3cUtuXR6rSEVfrQd+VUUZ3vo\nfoV67M3F72lUTlBZ7GJwJOm+WtwBXtxGo7EUPYNxKiu8nW5WDHa7pR/vWTESwnW16t1L2er/kCrM\nT4g3hf/MF4M98nXNcq+iTs72bV03+On9uwBobvT3M/NdKU5K+TEhxBuAK4A08B9Syt85eXEp5YeB\nD+f50ZXFLNILntvVyQt7upnWOKUsqpTBvJF2HuyjvrbKt6EWIz4r5xVDWPrSk52Gdd0gkcqUTarm\nJ/fvJpXW+F/XCt/aMwuh6Qa/e/wAEJ4Iz2hI40GdeFZZg+579X9FXu8zkbLFm4J9ztzeifaQH897\n5yf50PYc6edol9mV4GeboRc42qGklHcToCftB4+/2A7ADZcu93Xik5dV7gctYRA/JTuT2TGK4bet\npYJei8P7IJHKYBg+p2o8pM0SLbnlmlWhXP+Q1R66bP40ZoeggmYYBm0d5nvgV4X9eIrZUuz7vDak\n6EFS4XREIboGYtRUVzLNQy/Zycf23M5OAP76+jMDy6GXakfUS1b6gKYbvLCnGzDbjvzAD+nXb921\nDYDXX1X6xjyZxxCkVzzZO3Qi5K7WbZktpiwTgx5LpGmsq/b0YFRMTm/vUVM05bWX+TvuYaI1HTo+\nzMHjQ8yb6d+ccTdkrFbU6oAmd40nm44osxx6MpWhvrY60C6TdEbnnkfNLu0zQ5iKVyxO2tbK61PP\nw5GOYcbiaa48d6Hvm7JbnWKb4bEU+9uHWDZ/Gq84f1HRv+80CmG3Y6hQiDYSWCiyOLLtjmUSco8l\nMtR7GGouxvs0DIMHn2sDYNWi6Z6t4eT1FF5Q75BZK3P5OWqGR7MGvTqc9sdESqOiIthpZV6QSuu+\naFQU2rK3WOqQ81obWeth7t4vnLw7z/u+Cp+xi+GWzp3m30U8fjb3tplezvpVs3wVn7BD7rUBejIT\nCYLYucWpiqlX9Q+bFa5NAbasuDkYxhKZ0A4fx3pG2X/U7ERdGFIP/JDVljV3ppptWZmsLrjfozjJ\na62SKY26Kf7V5fhFKq15n46b5D0Yi5t70isvXFIW75eTO6pLCHGZEKK8RH9zGB4rrx5igJ8/sAeA\nJXP93RRtkYlAPPRJnge/tL/d0tlnz60OoGLc5aaRTGuMxFKhiSeNjJmFS6+6aKnvodGJzjzP7OgA\n1NUpt1XHwvKQk+mMMrMbiiGV0QPvgLFb/ILocMllIqdnszXbYyKc3FHnAY8DcSGEbv2nFb3CEAlW\n19mbmHuPpVl82Tnuei4nW83AiOl9LvYzeuEQezqYl0UvXmD3oM5sDn5udbE89oJZ/BlW618saRr0\n2SHM+Lax60KWLwhMhLIo7MldYY3iTFgeermRTGm+1NcUioYlAlbSLHSeNwyDr/zshYK/76RtbVbR\nq1KMIEQmvPRF0hmNodEUZ69s9b1YbXgsxbTGKYGGk/ORTGv82aombVHMs7JD7jOaw2/tm4wnt5kG\n/d03rfX0dZ0eU1WQVx6NpamvrVL2AGbPZ/B70IcZcf/LTy6R0kIZmGPjdLphLqPxNBlNp6HW231q\n0kJdhaSxM5qRVdOcCCfjU2cDbwWaMP/9VcAyKeXbvVhkEATZd+lFUZwtzeh3yFDXDXoG44F5U4Vy\nUDFLNKKmupKWqUGNvHSGvbawDz2ToesGe9sGWDRnKsvme+mdOjc8L+w2u0lafGy1tJnIMIzGUzSF\nZLCcPP6aHXIPyUNPpjKhGCg32SQ7bTrLj72qwKadtR0KiF3ZkadCOLmjfgNswByF2gi8lpP12JVn\nzBIk8DPM5GXBxGFrspffhU2JVIZkSqNVAU/GPnmqKNwQK5O2tcMdw8STGss9NebFsX2/Kc951nL/\nWnwmC0v2DyeVzZ9DTpV7CAY9o+lkNKPsQu4Ja9iO18/gZNt20BruhdqfvTLorVLKdwB/wDTuVwJr\nnC1PDXYd7gdg4exwqm6LZdCSY/U7B2jLygY9PCPfgVjlSWv72wepqAhWiKOUQM9Dz5vtYutWhtNe\nk85o9AzGWbN8pq9iSIVIpjUymq5cYWUumYBC7vlQWeO+EPGslG+w6+7qN2cSTA8oaliIpEcGfcD6\nvwTWSymHMAe2lA39wwnmzmwIxMPyoiRuwMrZzp3hvrKyUL7K70lu4ym0fak6Cz2RyjA4kqSmqjIQ\nGVU3V7AHolxy9nxvFlMko7E0hkFgKZN8ZHOeCh4MbQLz0Csq/uLwfLzHlDAtF00Fm7GEf1HWQnt2\nPJmhprpSCQlhJx66k0/1ESHEXcBHgD8JITYCiSLXFyqJZIZWvwuaPNzre6y++ekuZCudZACCNuiF\nsKfKqbCWXOwxoBetC8dIFsPwaIrqqgp/woMOTqrxVPgFcckymCSWyZhvZhgh9ye2HgPU6ySZjK2W\n0qfXw4YmU/hMpQNulSuwHFsyuBCT3lFSyk8BH7dmob8Z01N/veMFhow5WEMLZSZzqezY3wtAa4u/\nue2sSpwC8pjpEDTlbQoVMtoFiovmhjO1zCmptIZsG2Ba4xTPBTCcvlxQFcGFNmE7LKlyn3VGt5Xi\ngjfo9gjpay9YHPi13dDebUYWvC32nJxkSlMm2uNVyB3gQiHEvwB7gD4pZeHudoXo6jf7uYMqaCql\nJSOXvW0DHO8dY8XCZt9bS8IKc+d7j8Lw0J0YquHs7Gp1PRrDMPjMHZsBmONBmqZU7MNPmMWDZRFy\n93kkciHsZz7oELJbtu0ziy39mN5XaMtOpjUlWtbghANWCCda7v8KvBrTK68G3iWE+Krr1QXE/nZT\nhl5v+eEAACAASURBVNL3ObYevc4fNx8G4LWX+TNEJhe7HzYwI+ognKRCtCAXW1RmRtB54SIOht0D\ncXYe7KOiAv73Wzb6uKjC2GJIlQFJZObNM2bUHPCTi6YbVFTgu5JeRcVfZkqyh/gQ359iXR5bzdKX\nyX2TfATJdDge+kT6AZPh5FO9DvhrICGlHAauBa4vcn2hsc/Slb5w7dxgLuiyKu5olzn2cf2qVg8W\nUxj74Q4j9DeetKIbsR3hmROQ9GOxtjCd0fjyz7YA8O7XrmVuwBKVuYxZ/fqrFrX4e6EC75FdcBZW\nj7eTg1gmo4emEhdurUppB5hjViHfNecVP6TKGfk/M003iCfSgbWsQeF3yK70L4ST2Jidibf/1bWU\nUR/6Dqsvdr4PoZpcvMpbdg/EmTW93rXKlZNxrmEVxeV7fLKz0BU4XORiV9eqKCpzpGOYb9y1DXlk\ngMrKCi5Z71/hXj6PYTyHLU+qOcQq97AOqcU8/hk9PIN+4plXKxJWCDt3PNOHwuZCH9toLIVuqDNb\nIuGRQf818CtghhDiw5je+p3ulhYM8WSG/e1DVFVWKKFV7oR0RmNqQzA9vHYur6YqIJ3iAo/PifC/\nWhtNMtA5AM4ZGk3yd199FMMwc7Hf/tjVvkmdOjVWtkEPStgm3yEjzPy0U0wPPaAWyHERg+yBJ4Qe\n+FIJqy7ix/ftAvw5SExKnvOzJx66lPLfhBDXAUeAxcBnpZT3Fr/C4NluFVKsCXCOrds+9FRGD8xj\nVqltLZ1WM+SeFbRQpDDG5qV9PRgGtDbX8Yl3XhBqqN3m0PFh5rU2+t6pUFDPIOi6kBII8hkfTzxp\nyr6WwyhQmxMDUvw5VE+UJWmz0p83Xe5/PZNNoY9lMh13cKblfh/wI+CTUsq086WFzy8elAC8+pJl\nIa/EGYZhkE5rgXmp6bDyaXmLmcLz0AuFk5MpjeqqSuU8Pju3/77Xn83qxdNDXg0MWhXuydTkXoSf\nlIOHbmqphxPx6RuK+6/J4THJtH+zOAoZ0EQyw9SGGiUOy3BiBkAhnNz1/wa8CtgrhPiWEOJ8twsL\ngq7+GAfah6isgPWrymNgnKYb6Ia3XmqhGp3AleIKPDyqeuiybYDG+uA338keXbt4MpAZ7Uxe6zVk\nDc84/6yAik8h75ukUtQpH5puMBpLh6bUNjKWUroFMx9hhdzjIeqX5HvcdAcFl05C7k8ATwgh6oE3\nAL8RQgwB3we+LaVMFrnWQOi2WmhuuWZ1oAVNbvrQPS1YcaIUp1B4MpXtiVcntG1HMOwBDcHgLBTa\n3j1KdVUF81qDMOiTr2lkLAUEM/q2ULg4HfKs8cm47+mDpDI6y+b5X2cwvm3NMEyHIQwNeTf4rf43\n0Y6dTGUCL4grVGek69546AghrgS+CXwRuB/4e2Au8N9Ofj8Mdh7sA2BaUzAfiBcpqVTAXqpK3oz9\nb1ehhc5m1GrDuiBIr9MBmm6w7+ggM6bV+d7L7JS9bebIh6khVwT7NZXLK+Rh833aeMbswK9tG4Sg\ndAImpEifJ+nr4KaJ34t4UlNqiI3mwKA7yaEfAQ4CPwQ+KKWMW99/DHje3RL9w7D+8UF4DF6RDrjq\nPKwWlvE56+37e/jvJw8CKDHK1cb2zFUbZGEffuLJyYUmgqJ7wJw/4Le64WSMxs1IQVgGvdCWaxgG\nLx/so7Kygk1B6WLkYIdsgxgylI9SzxFJvyfE5fnQNN0gldbCu4/yrMmJh+5ktVdLKQ+M/6aUUgPC\nk6WahLh1EwRV0OCk73sygla5OnGACNcrfmqbqSR81rIZoU7qGs9ozDQOQQpLOME+qZ+1zL+Z48Vi\n30tiSXAFeuO3t2Ra466H9wEwsyXYwi8nVeMv7euhfzjBuWfMDqX40y6qUiWq45REyh6dGlxRXFun\n2YLpZMKZpxT4aJzk0J0MZ/kLY14O2E34qobe8jEwbJYjeBkCL1TBbeeIww5zb9/fQ31tNV94/8Wh\nrmM8tuyrX/3dpaKFME97sr0k2zERwOFwok243/q8qqsqmKdIZXIuL1gTw8SSoA5iJ49PDdtDL5Uw\niuI6+8warAvWqJNus5/7QqiTsPQYe5Rj0ML6bmaz/M/ThwBvBmw4eWSHrUKmoFXQxr9HfUMJ5rU2\nhiYqM9FntnWvqWMwe3rwBr3QfZTdmAPKhZbLKF7bk3vVpqVK9lnbmgaXn7MglOsrk0MvkmRao6LC\nv3srn9MTT5r1MyqlbB1E3B0NZ3m/F4sJmqA9dC+ekUFrstdV5y50/2IO2HWoH4BZARmsfJvsWDxN\nIqWFIq9YaNPfebAvOyhnpd/a5Dk4uY/sjbmqUp3zeCaEjonxHSVBjW8tlbBFiuxUTblVuSesEaZ+\nHNImesm0PbM+rANqnhO9V1XuHyxhOaFjh0xUqlKcjO37e2msr2H6tGDyf6m0RmVlRaijFI/3jmbX\nohK2jOnNV6xg9nQfpjy5IJsLVWhjTisg6JL0WVHMLfEw0oC5Ifcy9dATyUzgqdNMSKJbhT4ar4ri\njgohHgH+DMTtb0op/9nB74aGXdCkwuCRQrx8oJf7Nx9my+4uwIcCrAILSqQyoeca7QldG1YH38ZT\nCFu4xc+BJ6WS9bQUyoWqEHK3B+moJtNrY7fUBXXgGG8c7FSNSpEdJ8QSGV87TfKlt1TUM9B0D6Rf\ngWdz/qzODjIJ8ZTGojnBqGiVykPPtfH/frU1+/XyBc289VVnePLaTsJTiZRGa0twm1++FY3FTa8l\nDDW2QtjelEo5NBv7wQ62uKnwUTWd0amsgKogN8BxSzrSGax6XrHEU6aOelgHMTuyE7Y9L0Z8S9MN\nhsaSLPDtM83/WYR9QM2rFOdgxqkTpbjPCSEagRXAy0C9lHKs2AUGiWGYc2yb5kwN4+KO/+q9T5u9\n1+eeMZtbrlkd6BAZXTdIpsIXTrDFW1QbT2qHt1Q6odtkQ6cBGYbJrpLRdHYf7g+sXz/fYdUwDO58\nYA8ASwNQYZuIQo9/IpkJNXoQdpV7KZH+sXgaw4CpDf7tD/k+s7BmAhRUivOibU0IcTXwEvB7YA5w\nWAjxSudLDJ6BkSS6EawgSLE3azqjcfj4MFNqqvjsuzcFaszhRM467AKimBUmbQwxj58Pu8grvLD2\nxA+vaiH3Xz24Fwi3XsXu2KiuqqQ14B50p6Qyuu+T6HKp4OQKbhWLKSfDFgpqqvenaHbCojiFZLFt\nnITcnaz2S8ClwKCUsgO4Aviyu6X5y7a9Zr9nkNXJxfLte7aj6QZnr2wNpcUm4bM+ciFyQ252Dj0M\n8ZZC7/qJiuCgT+iToym2MR/tNkPd7715XaDXzT3y9A+bPeiv2rREyZY1MAutwhRx0gKO7HjBaMyK\n4PnooefjwT+3AX7JzU5OPmfcq2lrlVLKTvsLKeWuYhYmhLhQCPGo9ecVQognhRCPCyG+VczrFMOz\nL5vLPXNp8EpaTgLuhmHwzHZTHe2Dt6wPZT1Z9aUAvap8+2yYBr0QumJecC6qrW1wJElFBYHJmeb7\nV/cOmvW6QXWIlEIqHayHPp4TVe6hLaFogknJ/eUuOWIVVS+ZO83H6+ahwGdj1/UUwolBbxdC3AAY\nQogWIcSngDYnaxNCfBT4HmBXFn0Nc676FUClEOImJ69TLFul6aEHWTldjFdw5wOSsUSG+a2NvqmQ\nTbacMD30XEYTahr08EPuE6OH0E88UfrOMAw6+8ZobqoNPJqRi60lP2eGWi2GuaQyeqiqjBkFK7cn\n4+UDvYB/Bj3fE5TOaGi6wYZVs5QKuY8l0pOmSJ2s9n3AW4FFwAFgA/Aeh2vYD7wu5+tzpZRPWn++\nH3iFw9cpilRaY+XCZiU3YzA1nQH+4c3hSeGH4aHb5NqG4dFw1OomI6yQuxMGRkwBoqDqDgodVtu7\nR+kbSgReAzIe21iF6QEXwjAMMhmNKUEaiIqTpV9TafM9CrtuxilHOodztPmDU2u0hx7VKzaUKZbI\n0DjJmpyseL2U8s253xBCvB74zWS/KKX8rRBiSc63cneGEcDzclRNN6yCuHAMhJMi967+GLNnNHBG\nCCkBm2QyfA/dMAxelN1UVqgnAKRaWDsXO18clMJfIeww4KwAN1yb3FqME1XJ6n1eYB44dCPcIqtk\n2vysVD30jGf/0UEAXnnhEs4/c45v1xm/Zw9Zip1hzgHJK0ebyEwayZxwxUKIN2GGyv9ZCPGZcb/z\nSRwY9DzklulNBQYn+wUhxO3AZ51eQJWBIxOxeUcH/cMJ1q1oDXENx7n3KVM3PkxVLbsyuXV6g3KF\nOhnN7KsOa12FDobD1oYThlzueOx+3UCNRJ6PJKPbA2vUeu413eD5XZ387P7dADSHqGsQxpCTfDht\n7B20IlGb1s717XPNF316ZodZ3zQ3hPRNod1mLJFm9ozCB+dCu/k04GJMw3tVzvczwKecLnAcLwoh\nLpdSPgFcDzwy2S9IKW8Hbs/9nhBiKXAo39/PWJWAYVWTFhJNSKY1vvyzLQCsXxWMQc+3nh/ft5tj\nPabk6sIwxHesJe1tGwBg1cJwuxHyvUeablAZQhW5k1IM+yAUpnGwsQ/QgYaS86CF/NzbjL+Xvn3P\nSzzw7BEAls9v5s2vFGEsCzhRYBZeyN3Z4VjXDb565wts32fmz6dP9bfQcfzT3zdoRsDOP0udSWvp\njEY6o9NQW6KHLqX8HvA9IcQ1UsqHPVrXR6zXrAF2A3d79LpZQhMEcHCv3vfUIdIZnaXzpvGma8N7\nsMcSaebMaOArH7o80Pnj40/DsYQZAgxyhvZJFPjMNE1XSis9lyHLoAfpoU90TE2FqKiVu6ZMCCNl\nc5no+bdHpn7s7edx0br5gaZwxl/piW3HgPALYSejdzDOE1uPUV1VwfIFzYEr/9mTOkONgI174Oy9\nsmESRU0n8dakEOL3QBPmPVIFLJFSLnWyLinlEUxPHynlPuBKJ79XKipXcnb1mwJ7t964JtR1xJMZ\nZjbXBWrM82GLysxQsNVI0w2qFUsD2NihSCVC7mnboIcrIWxH5lR67nXdoHcwjlg8nUvXhzMyNRfb\n2blwzbyQV1KYpCV6dc35i/ngLRsCv352UmcIRXETFaBmDfokHrqTu//7wO8wjf+3gH3Ab4tYY6Bk\nDXq1epuxXd2+dH7AvY05aJbka5gFHzajivagQ3ghdyf0DsVpmVob2vz4XFJ2yL0m7JC7em2G+9vN\nEiFV6kPGEmnqplQp+bzlklWxDCrXPy5NcmLMbfh7pI09eGgyD93JUxiXUv4QeAwYwGxZu8Ld8vwj\nrDGOhTR4AfqG4hzrGaOyItyBH0HPic+H/fzYvfAqHC7Go2m6shXTsURamTY/e/MNpXI6Zx9WMTLX\nZY1wvnRDOBP7KipOzumPxtLKG3M44aEHcU/lc4gTSY3qqspQOxLGp7gef7Ed8MZDTwghZgAS2CSl\nNIBwZ24WoMcSmPC7kKJYXrIKPK48d1Gw0pTj7oxQZjJbjP9n273wKvbFaroRqrdXSMhlaDSlhPZ9\nPJnhQPsQEGxeNt/zE7Pua5XqHtKa3Roa/oHVMAy6+mNlYdDDnDMxPJZCtg1QX6vWnmQb9DOXFW51\ndnKnfQ34FfB64HkhxFuBLS7X5xs9lgTk3JnhKEZNtBE/uuUoAKsD0pef6NBwrNusbg87fw4n2mhU\nLNLRdIPKULy9wgbp0RfM+2hq0PnzcTd2Kq3xni8+yJAlDNTSFO4BurPXrE+ZPV0dpTi78l6FNMCu\nQ/2AmtGw8bRbe1RQnRO5d/aP7t0JwLxWdXxWwzAYGEmyYmEzG0Vh9dMJ3zFbEEZKeRfwSinlCHAu\n8Dbgrz1cr6fY7ViBS0AWeGb3tg2wzcqfX3P+4oAWlJ/ndpk694FrFOdgiybsOWxuMmF7MPnOYJqm\nK1cUp2k6X//FVgBuvnxFYNfNdzYcHElmjfkH/upszlgaUqeCRWdfjDkzGpQSTckopDb4iOVQvPri\nZSGvZHLxLTttGkTRZ26qtK1zmAefM1XNw1TxhJNTJdv2mrbDSdS50E76jBBiFHgQ+JMQ4lHLqG91\ntVKfOdIxDMCy+eHMRM6n8GMbrvfevI66kE/IdnFFGINrxlNjbb5hVWtPFMWIJzMMjqZYtVAdbw/g\n4HEzvF1dVcn61bNCXYt9H73mkmVcH5KRsJ+0/uEEg6NJLlyjTt8wgB56oV5F9j2yHZ2L14VX4e40\n02jfW3NmBuMl27azrcucGnjlxoUsnD01kGuPJ997dNz67JwI3Ux4dJRSLgBejRlevxl4zpqU9hkh\nxKbSlus/AyMJ6murA88VFbpXuwbM4pjQ+q1z6B8yRRNmNodfY6BrBlOqK5XyqgC6+2PousGyBeEc\nCvNhGAaf+vYzALzh6lUhrwYOWwfnBgX0rjv7zHB70P3KueQrirU9dBWKK/uG4syYVhu6Q+EEu0Ur\nkDqRnI9mYNhsB71AsYNh3EpNnutA/rbgpyulPIA5kOVHQogW4CbgHzGV4sJPwuZhNJ4OfHbuZHT3\nmwY9jPxebrzga3e+wNa9PdROqVKiyj2taUpuMAMj5qFnhgJ1Bja/fHBvtqDxdVcGF26fiCNZgx7+\ns2a/L6oVfJ3IoYcfcu8bSrBkXnhptmKw9SmCPizaz/10hZ57ONGZ5KTWqJCWezVwKfAq4DqgHngI\n+AwOJFvDYjSWDq0gDsibkO3oHaOhrprmpuBCy/l8gi27TdWqt73qzGAr7S3GXzOTMZRqM7KxhVvC\nLRw8cSO1d49w5wN7qKqs4BsfuSoUIzr+traLTy/fEJ5gip1nDLNzoxCaHq56HZiH53RGJ53RJ53U\npQp728z+/eDuc/M+sjukZilUWAknlOuc3N+F/sYAsBm4C3idlPKwB2vzlUQyQzyZCUXjupCBHB5L\n0dxUG4oRtUmkMozEUmxYPYubrwjfwwNTn1ilecNgGomf/nEPEI5Bz3eLHO0yc2gbz5jNojnB5/by\nhZO7+2NUVlaEkroZ/x7ZIVr1DHq4Ve72+5QMcVRysfz8j3s42jVCa3NdIFoLufdSz2CcigpoDTEd\nme9ZSySd63UU2k2/C8wCbgXeJYS4RAih1u47jj57rGQIoxwnwjAMRmIppjWEK9Npe50qyaxmNEMN\ng57jfvYNJbIpksUhdgLk8kOrlSZMb3g83QNxWlvqlajgtj8v1bQnsiH3kN8jW8BJRb2HXOLJDHc9\nvBeAD7xhfeAKe2PxNA211aF/XuPJhtzdGHQp5UeklOcAf4U52eyDwF4hxG+FEO/3ZqneYisztYZo\n0MeHJkfjaTKaEXrf9/b9prCNChO6bNIZLdSQe77twg7fXnz2PBbMCq/IKhc7FCiWhN+ZAGb73MBI\nQpmDsy0qE2RKywnZkHvI7Y+JMvHQt+3tQdMNXn/lykAnndk1PbFEmnoFakLg5Na+UbumwKWHDoCU\n8jhwJ/BN4A5gFWYeXTns4SfzQxAFmOiR7beiBtND8oz/f3tnHlxXeR3wn2Rbcmwt3vd94bMBY7Cd\nBBMgYMCFlBZapkOH0gUS2gZKyNBppqQlkGkD2YaGFpppk2GgyWQoZNomoQnQhrgYAjRJ2ZcP8AZC\nXiVrfZb0lts/7v3eu3qWrPf8pPsdPZ3fjMfS85V83nn3fuc75zuLO2fcE5U8nb56thc54HgdpTM5\nGR56DNd20uemME46kyWTzXHm2rliml109PQTBP6Sh4qPrsrxYJIk47mxTKimIO+hy2ngNNjtCYKA\nHz69i7se/F8ANqxJZrQ0DF6TUn0Z73kGQx25dfcOMGVybUkRlhMlxV1BOCXtXGAV8DzwU+Bqa+3r\nJynvmOKGfTR4DG8Xz0N+/Od7AQ8lNUU3Rk8q1M3KhXJKsTLZnLikuMQHQ4xA/nzY80ITv61f29UG\nyOmmJanjYFxPrr+877LMfiEh9+G2NXtau/jWD14D4My1c0fshjbaBITrdqo/w1IhHrojCAL2Heim\ncVpdSTlYJ1olbiI04J8FfmWtzY2SjGNGrzPoPspXhmryP5DhsWf3ALBty9KEBRpMdyrs6iWhpC8I\nArLZHLnAzxztE5Ff/Dwvws4wuCMAr/XeRfe2mxr4kQTDokMRBOG9VJjzLUdHUNgc+r7HpYfc3dp0\nxpo53HnD2cmenUdGsm8gSy4XiOirAAXHcO/+LgbSWZaW2IlxWOmttdtHR7TkaI36OUupR7334bCp\n3rYtS2n0nBTXk0ozeVKNXy8m9pym82NuZRl0n4Mh4PiQW8Ggy7inIey1XVtbw+olycwlKKY4TOqQ\nshg7DkbJenUex9wGAQJD7oNxx1yb1833lpDmat+9Dz0qev5/+eZBgJJzCmStphXizomlJH69te8o\nANf/xmneZHARwO7UAA0lhm2SIBP1a54iLOSe5OjGUpBYktV6pIf5s6Z59zwBOnvC6o2LP7xMzL3t\naD3Sy+RJfkr7oFAC1Z+faijnHorT3+93Ex0EhedsmhBn0HE0qk46bVVpuU/+n8hRxHlXPnqDF9cP\ndqcGONJxjA2r5/ipiy/6vjslZ4Z2QMFDl2AU4nkPUkLujryHLsSgu/GtUsofXY/yWQJaGRfTeyzN\nknmN3sugpHvoqf7QO/Yhn1sne/MeuoznzLGrJWyyU+ra7X81HSX2HeiivUvGgIae1ABf+044YXb1\nEhlJaP0DGe9eXnzT4yYq+Qy5D+XQDQjz0J3B8jXApkC46cknegnZiLkqksVzZSToOXK5gGN9aRHH\nAJ2RlyfRoL+66wj3Pfoy4PeotPVweFwrZaMKYfWGG3s7o0Sn0P9TOUq8814Y3t5SQgP7sSQIYOdL\nH/BiNPIuyRKM4QiCgHQ2J8ZIgfyQu++MYMfb0bGNz/sovvEZSAvI3I7J40a4Nk2Xcczm6BvIkAv8\n5/MEAew9EPbdn5f0SOkSeGN3WDExb+aH2LjW0wTBIOCNPaEcpYa2x4q4j+Fywj62cVHJJZmyVtMK\nGIgMhK+khvii1xVlbd56zSbvmcAQ1sMGgSzjKTUpzo1QrPeYyASF3Ad3hialLn4gIyNz29HZG+pH\nWlOZ3mMJTgwbjmhNcl0ilwvpfBgv7XPdPW//5NleIohu3XaZ9lKeMyhEw+aX0VtexlM5CjjPYcoU\nv28pCCAVPcwi6nSDsDkJ+NdNnqAQcpdiGBx7Pgi9mUXeusQNPgc42t1H47QpYur10xI89Bi7PwgT\nYaUsxEG0FXNNrnw1lIrTnRqgaXpd4q1Uj2OI/95GESifE84C5LTpdQSxNbKc4T7+D3hGCWe0vJWI\nxHTuGtz43J3HM37zYVLPXmc8iuF2n1KiBm2dx/jRzt3sbu1kzZJm5s6UYSDaOvuYLyBU6rwqCR56\nPBejo7uf5oY6733ci5fcN/eGZ5+rFvn3its7+5ghrM89hBUKe1s7WTy3wVuOiLuX3Ho02fOmJ75u\nu9bB5Wzmq8ag5z10AR7fgbZwdy7FKOS9YSkeOjKS4uLcdv+z7I8+t7NPX+hZmpAgCIQkMxZwk5+k\nNCkJ2wfLiBY4Onv6+ZcfvwnAqsV+k2JTfWl6+zIsnifj84rzvSfeIheEcxN8lxxK89Ah1jp4Invo\nvo1WNpfjtV1HmNVUL2jRi3Qj5GYNCMhmwptVwgasrauP/W291NbW8OUbz2XtMj8NU4rJZANyAdQJ\n2ogd7gibpUg4sw6i4yRp9dUvRQmxS+c3ehl366ihkIMhbcY3hF3QAM45Y5FXOYIAMnlvWEYvg4CT\ni2LKWSkqpKs3TGrwFVZ2t8E773WQCxAV4nIJg77PPeOPSlpEyD2UyLUM3v7R5axfOUvMebW0EjqA\n7z/1DuD3zDruzKUzOREldHE6omY3v3fpOu+ep8PHwKqR6E4N0DhtCms8dRwE8otSNhtQW1sj5vOC\ncKohlBc1kPUkVIAbMen7xnUP83lCZlcHyOkpHUdSyF1i440gkGnQWw/3UjdlEls3yDiWGBA4sc+1\n6/TfO6CAJFkc3b1p7y2xHZlszvv5eTEu5F6OXLKehApo6zpG47Q67yMUX3k3DLeV2gggCVy3MSkL\nXxDIqkOXMo0KBnuf/YImvwWEtdU9x9KsXzFTxHFSQCDyDN2Via1Y6DchLn4v+eq7PxxBEIQeuoiN\nRkAmmxN1fu4GWEF5Hrr/p7JCOnv6efhJy/sHe1i+wON5VfT0HOkIIwVSOsQBtEbdxhZ7K8WKiK0w\nksrW3FmVBCMVx1VLeO82Fn1uLQfD+2j+LBnh22wuIJcLRNxDju/+5C327u9i9ZJmMd4neJpAOQwB\nYe/0bC7wriO3ImWygYjz86EqgcqRS9YKViYH21N8+u8ez39//llLPEoT8m5LJ7Obp3rfncc50BYm\nMi2a49mgxyjcrHIWY0khdwjLjQBmN8uolnBNXBbMlpFgVZiDLmcZe/KFfQBsWee3Y2UxvislYHC5\noWvkIuEoIAjgWF9alIcOkIoiq+VEnf1/yhXw7zveBcJQ6b23XuDfA41YuahZVHLFgajJxXwhCzHE\nPXQ5RlSSQW851M19j74E4G1aVzEvvH4AkLEIAzz36n4Aptb7/9ziz/tX/+w81q+c5VGa45G06YGC\nQffuodfUkM0FHOnsY+1SWccSbR3hhn5OGRt6WZ9ymdh97UAj3/78Jczw2GmoGEkN/ntSA/z8lf3U\nTa712o1pEAH5oRo+w1zFey7f+Rdx/uN/duW/Pn2V/3kAAL2p8Ahg9WJZC9+Fm5f6FmEQa5ZKOW4r\n3OASNj1xunvDe6lxupyjAEnrNhQ2rOVE6OSsYCdBdyrNsqXTRRlzgOUL/Z3lO5yxeublViBscOE7\nauD++6d+9T5P/fJ9QJYRXTJPwOcW+3rypFoe+OtLRLQPDYKwXr+mBlZ67n4Wv4//4BPrvQ9kKkZS\n1AnC505CYmUcGw3TahKUZyBhkFb8WMJVAM2ZUfrzn/hqaoypAf4R2Aj0AZ+y1u4+2d8nJbSVi00c\nuPxjqzxKMjSfufos3yLkccb8snNWcOpKv9ONHM0NdV6TKofiwS9sp1lAtUQNYQOX13e3MaOxtsxa\nSgAACOJJREFUXtRZo4SzYQjHpUplat0k75v5OF9+6Bf5ryU1vJHmofceS7Nw9vSyPjsfT+aVQL21\n9hzgNuCeSn7ZOULqYbujxjZrljT7H4JQxJlr53rtWDUUs5qmcuNVG8VkKF93+WkiFr1szDBIMOYO\nt1+VMgTFISXvYX806tL3mfBQSOmk9/6h7vzXyxY0csGmJWwy8zxKNJh1y2U4hwA/eHoX7V19TJ9W\n3pGEj0/6XOBxAGvtC8aYLZX8Mik7PNeprknIIux62wP87nbjUZICcXMpxbNyLJgtoxTLJQtJRcI4\n4Ph9JOXIxo3d3bZFznm+259K2fS8HYXZAe7/i20eJTme2hqYJSH5NPrM3Ab6ws3lVW75cI+agM7Y\n9xljzEnLMU/IAJTOyKA3C8kAXjhnOjMa6jn/rMWctkpGWDvOh0+Vde65foWM3blrSiKhJtbxzvsd\n+a+3rJfjUYGc7G1XhrlM0LGNMwpSdHTVBWsB+Js/2epZkgJu09PUUM8kYZHVqy8+hd88b3VZP+Pj\nk+4C4nd9rbU2N9zFwzAJYMbUATraD9HRPmqynTQfqukhnepgKjNoaWnxLQ4Ad99wOoAYeTraD5FO\nRR9W/2zvcuWCIC9Pa+sHXmVxHOtpI51qZ9HCJu/6cSxsSvPegW42rJnDVHpoaen1Kk97Z1/+c+vp\nOExLS9qrPAAz6/s4dDTFzPo+MZ/bQKqddKqPXF9WhEwr58J9t2xiUu2ACHkAUt1tpFNHqW9oECHT\n4YPd+Xt7+ezckDJddNFFK4AWa22m+N9qgiDZZA5jzG8Dl1trrzfGnA3cbq399RNcfydwR1LyKYqi\nKIpwVlpr9xa/6MOguyz3M6KXrrPWvl3m76gnzJBfA2RHuHwPsLJcOceQzwLf8C1EEdJ0BPL0pDoa\nGdXRyKiOSkOaniTpaA8wZSgPnSAIxuWfU045JRjN6xKU+07fMkjXkUQ9qY5UR6qjiasnSTo6kW5k\n1AxNLHb4FmCcsMO3AOOAHb4FGAfs8C3AOGCHbwHGATt8C1AKatATxlq7w7cM4wHV08iojkZGdTQy\nqqORGS86UoOuKIqiKFXAeDboXxzl6yYyqqORUR2NjOpoZFRHpaF6Gp5hdZN4lruiKIqiKKPPePbQ\nFUVRFEWJUIOuKIqiKFWAGnRFURRFqQLUoCuKoihKFaAGXVEURVGqABlz9crEGDMZeABYAdQBXwLe\nAB4EcsBr1tqbomtvAP4YSANfstb+Z+z3rAOeB+ZZa2UPoi6TSnUUjbS9B9gM1AN3Wmt/nPDbGFNG\nQUdNwMNAA+FsgWuttYcSfhtjSjk6iq6fCzwDbLDWDhhjpgLfBeYRTlr8Q2ttW5LvYawZBR01Eeqo\nCZgC/Lm19vkk30MSVKqn2OtVu25Xynj10K8FjlhrzwcuBe4jND6ft9Z+HKg1xlxhjJkP3Axsja67\n2xgzBcAY0wh8nXAhrkYq1dHvA5OttecBVxIOwqk2KtXRHwGvRD//CPA5D+9hrClJRwDGmO3AE0B8\n2P2nKejoO8DtSQqfEJXq6Fbgv621FwDXAfcnKHuSVKqnibBuV8R4NeiPUFgYJgEZYJO1dmf02k+A\nS4CPAM9YazPW2i7gHQpT3v4ZuA1IJSZ1slSio43ArwGtxpjHCHX1oySFT4hK76NXCb0qor+r0Vso\nRUcXR19ngYuA9tjPnws8PsS11USlOroH+Kfo6ynAsTGV1h+V6gmqf92uiHEZcrfWpiC/W3sU+CvC\nXZujm3CBbQQ6Y6/3AM3GmDuAx6y1r0bjXKuOCnXUBMwBVltrLzfGnE8YFvv42EueHJXeR8ARYLsx\n5nVgJnBeAmInSok6ao6u/Wl0bfyZaqKgO6fPqqJSHUWbRIwxCwijGJ9JRPCEqVRPE2HdrpTx6qFj\njFkKPAU8ZK19mPAMxtEIdBCe2TUN8fq1wCeNMT8DFgBPJiJ0wlSoozbgMQBr7dPAKUnInDQV6ugO\n4CvW2tMIIxr/lojQCVOijuLE2092RdcMd21VUKGOMMZsAP4L+Etr7TNjKatPKtTThFi3K2FcGvTo\nTPMJ4HPW2oeil1+MPEmAy4CdwC+Ac40xdcaYZmAdYeLFWmvtNmvthcABwrBqVVGpjgiTUT4R/a6N\nwL4k5U+CUdBROwXv8zAFw1U1lKGjOHHv6Vmi+yj6u/jacU+lOjLGnEoYjr7GWlu1RqpSPU2EdbtS\nxmXInfAMZQZwuzHmC4S7uFuAf4iSld4Evm+tDYwxf09onGoIky+KzzkDBi9A1UJFOjLGfAv4pjHm\nuej3/Wnyb2HMqVRHXwC+bYy5ifBZ+pSXdzG2lKSjop+Je1XfBB4yxuwE+oFrxl7kxKlUR3cRVpLc\nG4WSO6y1vzX2YidOpXoqfr0a1+2K0OEsiqIoilIFjMuQu6IoiqIog1GDriiKoihVgBp0RVEURakC\n1KAriqIoShWgBl1RFEVRqgA16IqiKIpSBYzXOnRFUUYZY8xy4G3gdcIa36nAK8DNJ5oiZ4x5ylq7\nLRkpFUUZDjXoiqLE+cBau8l9Y4y5i7DZx/nD/wgXjLVQiqKMjBp0RVFOxB3AgajX+M3A6YSzzS1w\nFfAVAGPMc9barcaYS4EvEq4te4AbrLVHvUiuKBMMPUNXFGVYrLVp4F3gCqDfWnsOsBaYBlxmrb0l\num6rMWYOcDew3Vq7mXB4xlf9SK4oEw/10BVFGYkAeBHYY4y5kXA4zRqgIfbvAB8FlgE/i3qS1xJO\n7VMUJQHUoCuKMizR0AwDrAb+FvgG8AAwh+OHY0wCdlprr4x+to4qnECnKFLRkLuiKHHiYz1rCM/D\nnwNWAf8ajb08RJgkNym6NGuMqQVeALYaY9ZGr98BfC0pwRVloqMeuqIocRYaY/6P0LDXEobarwGW\nAN8zxvwO4RjU54CV0c/8EHgZ2AxcDzwSGfgW4NpkxVeUiYuOT1UURVGUKkBD7oqiKIpSBahBVxRF\nUZQqQA26oiiKolQBatAVRVEUpQpQg64oiqIoVYAadEVRFEWpAtSgK4qiKEoVoAZdURRFUaqA/wde\nlEtIzyKueQAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x8a5fa90>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig, ax = pyplot.subplots(figsize=(8, 3.75))\n",
"df['cum_precip_inches'].plot(ax=ax)\n",
"seaborn.despine(fig=fig)\n",
"ax.set_ylabel('Water year cumulative precip (in)')"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.4.3"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment