{ "metadata": { "name": "", "signature": "sha256:fe87cbab75175db4de86cade6dc5fb48ab4c19b65e7db7c62b8bb80c9f9a9fe8" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline\n", "from matplotlib import pylab\n", "import matplotlib.pyplot as plt\n", "import netCDF4 as NC\n", "import numpy as np\n", "import datetime\n", "import requests\n", "\n", "from salishsea_tools import nc_tools\n", "from salishsea_tools import tidetools\n", "import os\n", "import csv" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Grab observations at Point Atkinson from 2013" ] }, { "cell_type": "code", "collapsed": false, "input": [ "yr=2013\n", "\n", "statID = '7795'# station ID number\n", "start = '31-Dec-' + str(yr-1) #start date \n", "end = '02-Jan-' + str(yr+1) #end date\n", "\n", "tidetools.get_dfo_wlev(statID,start,end)\n", "# now there is a csv file in this directory with the sear surface height information. Let's take a look at the plots.\n", "filename = 'wlev_' +str(statID) + '_' + start +'_' +end +'.csv'" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "raw", "metadata": {}, "source": [ "Now go to matlab to run the tidal predictions. \n", "\n", "Tidal predcitions: 'Point Atkinson_t_tide_compare8_31-Dec-2013_02-Jan-2015.csv'\n", "\n", "Harmonics: 'Point Atkinson_harmonics_31-Dec-2012_02-Jan-2014.csv'" ] }, { "cell_type": "code", "collapsed": false, "input": [ "yr=2014\n", "\n", "start = '31-Dec-' + str(yr-1) #start date \n", "end = '02-Jan-' + str(yr+1) #end date\n", "\n", "filename = 'Point Atkinson_t_tide_compare8_' + start+'_'+end+'.csv'\n", "\n", "\n", "import datetime\n", "import pandas as pd\n", "import pytz\n", "\n", "def dateParserMeasured2(s):\n", " #convert the string to a datetime object\n", " unaware = datetime.datetime.strptime(s, \"%d-%b-%Y %H:%M:%S \")\n", " #add in the local time zone (Canada/Pacific)\n", " aware = unaware.replace(tzinfo=pytz.timezone('Canada/Pacific'))\n", " #convert to UTC\n", " return aware.astimezone(pytz.timezone('UTC'))\n", " \n", "ttide = pd.read_csv(filename,skiprows=3,parse_dates=[0],date_parser=dateParserMeasured2,delimiter='\\t')\n", "ttide = ttide.rename(columns={'Time_Local ': 'time', ' pred_8 ': 'pred_8', ' pred_all ': 'pred_all'})" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot and adjust axis. " ] }, { "cell_type": "code", "collapsed": false, "input": [ "\n", "\n", "fig,ax=plt.subplots(1,1,figsize=(10,5))\n", "ax.plot(ttide.time,ttide.pred_all,'-k')\n", "ax.set_xlim(datetime.datetime(2014,9,9),datetime.datetime(2014,9,16))\n", "ax.set_ylim([-3,3])\n", "ax.grid()\n", "ax.set_ylabel('Tidal elevation (m)')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 4, "text": [ "" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAnsAAAE4CAYAAAAjEDXIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXl4FGW2/7+dhZB9IftOQkhYAmFfFCeA4LC7cL3qHUdw\nG0eQq4PjNuN1GxkdBxX14rgNjl51HEUcQFRAE0BWJQlbSMjW2fed7EnX74/8KnSSXqvet6o6ns/z\n8Dx2d1X18U1X1bfO+b7n1QmCIIAgCIIgCIIYkTipHQBBEARBEATBDxJ7BEEQBEEQIxgSewRBEARB\nECMYEnsEQRAEQRAjGBJ7BEEQBEEQIxgSewRBEARBECMY1cReZ2cn5syZg5SUFEycOBGPP/64WqEQ\nBEEQBEGMWHRq9tlrb2+Hh4cHent7cfXVV+Ovf/0rrr76arXCIQiCIAiCGHGoWsb18PAAAHR3d6Ov\nrw8BAQFqhkMQBEEQBDHiUFXsGQwGpKSkICQkBAsXLsTEiRPVDIcgCIIgCGLEoarYc3JyQlZWFsrK\nynD48GGkp6erGQ5BEARBEMSIw0XtAADA19cXK1aswE8//YTU1NSB9yMiIlBRUaFeYARBEARBEDYS\nHx+P/Px8tcMYhmqZvbq6OjQ1NQEAOjo6cODAAUybNm3QNhUVFRAEgf5J/PfUU0+pHoOj/qOxo/Gj\n8XPcfzR+NHZq/SsoKFBDUllFtcxeZWUl7rjjDhgMBhgMBtx+++1YvHixWuGMSPR6vdohOCw0dvKg\n8ZMHjZ88aPykQ2M3MlFN7CUnJyMjI0OtrycIgiAIgvhZQCtojGDWrVundggOC42dPGj85EHjJw8a\nP+nQ2I1MVG2qbA2dTgcNh0cQBEEQBDGAVnULZfZGMNTKRjo0dvKg8ZMHjZ88aPykQ2M3MiGxRxAE\nQRAEMYKhMi5BEARBEAQDtKpbKLNHEARBEAQxgiGxN4Ih74V0aOzkQeMnDxo/edD4SYfGbmRCYo8g\nCIIgCGIEQ549giAIgiAIBmhVt1BmjyAIgiAIYgRDYm8EQ94L6dDYyYPGTx40fvKg8ZMOjd3IhMQe\nQRAEQRDECIY8ewRBEARBEAzQqm6hzB5BEARBEMQIhsTeCIa8F9KhsZMHjZ88aPzkQeMnHRq7kQmJ\nPYIgCIIgiBEMefYIgiAIgiAYoFXdQpk9giAIgiCIEQyJvREMeS+kQ2MnDxo/edD4yYPGTzo0diMT\nEnsEQRAEQRAjGPLsEQRBEARBMECruoUyewRBEARBECMYEnsjGPJeSIfGTh40fvKg8ZMHjZ90aOxG\nJiT2CIIgCIIgRjDk2SMIgiAIgmCAVnULZfYIgiAIgiBGMCT2RjDkvZAOjZ08aPzkQeMnDxo/6dDY\njUxI7BEEQRAEQYxgyLNHEARBEATBAK3qFsrsEQRBEARBjGBI7I1gyHshHRo7edD4yYPGTx40ftKh\nsRuZkNgjCIIgCIIYwZBnjyAIgiAIggFa1S2U2SMIgiAIghjBkNgbwZD3Qjo0dvKg8ZMHjZ88aPyk\nQ2M3MiGxRxAEQRAEMYIhzx5BEARBEAQDtKpbKLNHEARBEAQxglFN7JWWlmLhwoWYNGkSJk+ejNde\ne02tUEYs5L2QDo2dPGj85EHjJw8aP+nQ2I1MXNT6YldXV7zyyitISUnB5cuXMWPGDCxZsgQTJkxQ\nKySCIAiCIIgRh2Y8e9dffz0eeOABLF68eOA9rda+CYIgCIIghqJV3aIJz55er0dmZibmzJmjdigE\nQRAEQRAjCtXKuCKXL1/G2rVrsW3bNnh5eQ37fN26dYiNjQUA+Pn5ISUlBampqQCueAvotenXr776\nKo2XxNfGvhUtxONor2n8aPxo/BzztfieVuLR+mvxv/V6PbSMqmXcnp4erFy5EsuWLcODDz447HOt\npkMdhfT09IEfJmEfNHbyoPGTB42fPGj8pENjJw+t6hbVxJ4gCLjjjjswZswYvPLKKya30eqgEQRB\nEARBDEWrukU1sffDDz/gmmuuwZQpU6DT6QAAf/7zn/HLX/7ySnAaHTSCIAiCIIihaFW3OKn1xVdf\nfTUMBgOysrKQmZmJzMzMQUKPkI+xp4CwDxo7edD4yYPGTx40ftKhsRuZqCb2CIIgCIIgCP5ops+e\nKbSaDiUIgiAIghiKVnULZfYIgiAIgiBGMCT2RjDkvZAOjZ08aPzkQeMnDxo/6dDYjUxI7BEEQRAE\nQYxgyLNHEARBEATBAK3qFsrsEQRBEARBjGBI7I1gyHshHRo7edD4yYPGTx40ftKhsRuZkNgjCIIg\nCIIYwZBnjyAIgiAIggFa1S2U2SMIgiAIghjBkNgbwZD3Qjo0dvKg8ZMHjZ88aPykQ2M3MiGx9zNE\niylmRyE/Px+ZmZlqh0EQBEEQNkOevZ8ZH374IV577TX8+OOPaofikGzYsAH/93//hx9++AHJyclq\nh0MQBEFoCK3qFsrs/Yw4deoUfve73yEvLw+lpaVqh+OQnDx5EuvXr8eKFStQVlamdjgOzdGjR7Fr\n1y61wyAIghjxOITYO3fuHFasWIG+vj61Q3EojL0XlZWVuOmmm/Duu+9iyZIlSEtLUy8wB8CUb6W9\nvR3Z2dn485//jE2bNmHZsmVobm5WPjgHwJrvRxAEbNq0CU8//bQi8Tga5JuSB42fdGjsRiYOIfY+\n//xzfP3119ixY4faoTgkXV1duPHGG3HvvfdizZo1WLRoEb7//nu1w3I4MjIyMGnSJLi7u2Pz5s1Y\nuHAhbr/9drXDckgOHjyIzs5OVFZWorCwUO1wCIIgRjQO4dmbNWsWbrnlFvz1r39Fbm4ufHx81A7N\noXjppZdw5MgRfPnll3ByckJubi6WLFmC4uJi6HQ6tcNzGLZu3Qq9Xo/XX38dQH+mLzAwEM3NzXB1\ndVU5Osdi0aJFWL9+PQ4dOoQJEyZg8+bNaodEEAQhG/LsSaS6uhp5eXkDZbMtW7aoHZLDkZWVhbVr\n18LJqf/PPX78ePT19VFGxU5OnDiBOXPmDLz28PBAREQE8vLyVIzK8Th58iQKCgpwyy234MYbbyTf\nHkEQBGc0L/a+/fZbLF68GK6urnj++efx7rvvkkixEdF7cenSJSQkJAy8r9PpqJRrBVO+lRMnTmDu\n3LmD3ktOTsa5c+cUispxsOT7efHFF/Hwww/D1dUVixcvxvnz51FdXa1ccA4A+abkQeMnHRq7kYnm\nxd6+ffuwfPlyAEBYWBh+97vf4ZFHHlE5KsdBEATk5eVh/Pjxg94nsWcf5eXl6OjoQHx8/KD3k5OT\ncf78eZWicjwuXryIo0eP4q677gIAuLm54Ze//CX+/e9/qxwZQRDEyEXznj1/f3+cP38e4eHhAICO\njg7Exsbi+PHjiIuLUzlC7VNbW4vExEQ0NDQMer+4uBizZ89GVVUV+fZs4IsvvsB7772Hr776atD7\nn332GT766CN8+eWXKkXmWNx9992IjY3FH//4x4H3Pv30U7z//vv4+uuvVYzM8SguLkZ5eTk8PDzg\n4eGByMhIeHh4qB0WQfys0apnz0XtAKwRExMzIPQAwN3dHbNnz0ZmZiaJPRsYWsIViYmJgaenJ7Kz\nszFp0iQVInMsTJVwASrj2ssPP/yAhx56aNB7y5Ytwz333IPm5mb4+vqqFJnjccstt6C9vR0AUFdX\nh0WLFuHDDz9UOSqCILSI5su4YgnXmJSUFJw5c0aFaByL9PR0kyVcESrlmmeob+XkyZMmxd64ceNQ\nWVmJtrY2hSJzDEz5fnp7e6HX64eVwn18fLBgwQLs27dPoei0jzXfVFdXF86ePYtjx47hzJkzOHjw\nII4ePapMcA4A+c6kQ2M3MtG82Fu2bNmw96ZOnaopsXf8+HG88847aodhkry8PJOZPUCbYs9gMKCl\npUXtMAbR29uL06dPY9asWcM+c3FxQWJiIi5cuKBCZI5FSUkJQkJCMHr06GGf3XDDDTQr1w4yMzOR\nmJgIT09PAEBiYiIaGxtpogtBECbRvNgzlU3RktgzGAzYsGEDnnnmGc3V6VNTUy2KvYULF+LQoUOa\nWpnk2WefxerVq9UOA6mpqQP/fe7cOURHR8PPz8/ktlTKHY7x+IlY+i0uX74c3333HeeoHAdT42fM\nUFuBk5MT5syZgxMnTnCOzDGwNn6EeWjsRiaaF3suLsNthfHx8aivr0djY6MKEQ1m586dcHJygqen\nJ06fPq12OMO4dOmS2TJuWFgYxowZg5ycHIWjMk1FRQVef/11ZGRkoK6uTu1wBjBXwhUhsWcblsRe\nWFgYenp6hk0kIkxjykM6b948EnsEQZhE82LPFE5OTkhOTsbZs2dVjaO3txdPPvkktmzZgjVr1miu\nfURaWhry8/PN3mCB/vJPfn6+glGZ5+mnn8Zdd92FxYsXq+7fMvatDG2mPJTJkydrsv1KS0sLsrKy\nVPluU74fS2JPp9MhISFBM79FtbHmmzIl9ubOnYvjx49zjMpxIN+ZdGjsRiYOKfYAbZRy//GPfyAs\nLAxLlizRpNirr6+Hp6enxeXlxo0bp4kb7IULF/Dll1/iiSeewKpVq7Bnzx61QxrgzJkzmD59utnP\ntZrZ++///m/ccMMNmrEXWBJ7gHZ+i1qnsrISra2tw8Zyzpw5OH36NHp7e1WKzHG5/fbb8f7776sd\nBkFwg8SeRDo7O/HMM89gy5Yt0Ol0mDt3Lqqrq1FUVKRaTEMJDAw0W8IV0coN9tFHH8Xjjz8OPz8/\nrFixAgcOHEBXV5dq8Rj7VkpLSxETE2N224iICHR3d6OmpkaByGzj8OHDOHjwIIB+M7/S2OvZA4CE\nhARaeu7/Y8k3dfLkScyZM2dYf0w/Pz9ERUVp8sFDaezxnRUVFWHv3r34/e9/j/Lycn5BOQjk2RuZ\nOLTYU6tEBQBvvfUWUlJSMG/ePACAs7MzVq1apansnrWbK6ANsZeWlobs7Gzcf//9AICQkBBMmDAB\nhw4dUjUuAGhvb8fly5cRGBhodhudTqepUm53dzd++9vf4tVXX8XNN9+MnTt3qh0Senp6UFpaarE3\nphZ+i46AuZ6PAPn2pPDOO+9g/fr1uP/++7FhwwbNZMIJgiUOK/aSk5Nx8eJF1UoWX3zxBTZs2DDo\nPa2Vcr/77jurYi8+Ph4FBQUKRWSa7du344knnoCbm9vAe2qXckXfSnl5OSIiIuDkZPlU0VIp95VX\nXkFsbCxuvPFG3HTTTdi5c6fiN7Chvh+9Xo/w8HCMGjXK7D5azuy1tLQoKqIs+aYsiT3y7fVjq++s\nu7sbf//73/Gb3/wGTzzxBC5duoTPP/+cb3Aahzx7IxOHFXteXl6IjIxEbm6uKt+fnZ2NKVOmDHrv\n2muvRUZGhmZmFJaVlVkt48bExKC8vBzd3d0KRTWcixcvYubMmYPeW716Nfbs2aP6U3ZZWRkiIyOt\nbqcVsafX6/HSSy/h9ddfh06nw6xZs9De3o7s7GxV43KULLM5Nm3ahBUrVqCzs1PVOMSej7Nnzzb5\nOWX27OPLL7/ExIkTkZiYCDc3N7z33nv47//+b81cwwmCFQ4r9gD1fHu1tbXo7e1FaGjooPfd3d2x\naNGiYeunqkVjY6PVG+yoUaMQGRkJvV6vTFBD6OvrQ0FBwbA4J02aBJ1Op5qAEn0rpaWlDiX2Xnzx\nRWzYsGGgXKrT6XDjjTcqXsod6vuxRewFBweju7tbEy2VjPnmm2+Qnp6OCRMmKNb42Zxv6vz584iK\nijLb83HChAmorq5GfX09x+i0j62+szfffBP33XffwOt58+bhpptuwrPPPsspMu1Dnr2RicOLPTV8\nexcvXsSECROGGaQB7ZRy+/r6UFhYiHHjxlndVs2Mil6vR3Bw8MBKACI6nW4gu6cmZWVliIqKsrrd\npEmTcOHCBRgMBgWiMs/Zs2exePHiQe+JpVw1sUXs6XQ6zWX3Wltb8Zvf/AZvv/02HnjgAbz77ruq\nxmOphAv0e4dnzZpF2T0byMnJQU5ODq6//vpB7//617+mUiYx4nBosafWGrnZ2dmYOHGiyc9WrlyJ\n/fv3q37TLy0thZeXFzw8PKxuq+YNNjc3F4mJiSY/U9O3J17sbS3j+vv7w9/fX7UMqUhOTs6w8Zw/\nfz6qq6sV/RsPvVnaIvYA7fn2HnvsMSxevBhLly7F9ddfj7NnzyricTUnNqyJPUD7pdz29nbuvkJb\nxNpbb72FO++8c5iPNDk5Gbm5uaqX7NWChO7IxKHFnlplXEtiLzAwEB4eHqisrFQ4qsHk5eXZJFIA\ndSdpWBJ711xzDc6cOYO2tjaFo7qCrZk9oF80q9l6p66uDgaDAcHBwYPed3Z2xg033KBqdi8vL0/z\nWeahHDlyBF9++SW2bt0KAHBzc8Ptt9+Ov//976rFZIvY0/okjSeffBILFy5ESUmJajH09vbigw8+\nwD333DPss9GjRyMhIUEzs+sJggWqir0777wTISEhSE5OlrR/ZGQkuru7UVVVxTgyy1gSe4A2Zrjm\n5eWZNXEPRauZvVGjRmHs2LEoLCxUOCr7PXsAEB0dreoNTMzqmbIXKF3KNfb9dHd3o7y8HGPHjrW6\nn5bE3jvvvIM//OEP8Pf3H3jvrrvuwo4dO7h3ATDlm2ptbUVpaSkmTZpkcd/Zs2fjp59+Un1ykynO\nnTuHDz74AOvWrcOTTz7J7Xus+c4uXbqEgIAAxMbGmvx8+vTpqvSn1AJyPXsdHR2qV7aI4agq9tav\nX49vvvlG8v46nU6VUq41saeFG5alNXGHolWxB6hf1rO1jAuoL/Zyc3ORlJRk8rPU1FTk5OSoMsuw\nsLAQUVFRcHV1tbqt2n9vY7Kzs4etnDJp0iTExsaqspxfQUEB4uLi4OzsbHE7sSek1ia6CIKADRs2\n4Nlnn8Vf/vIXfPvtt6r1Ss3KysK0adPMfj5t2jRkZGQoGJE0BEFARUWF2mEMkJeXh6SkJLz88stq\nh0IMQVWxt2DBgkFPzVJITk5WNN3e1NSElpYWi6U9rWT2bPWcxMXFobi4GH19fZyjGo41sTdu3DhV\nbv7p6eno7OxEa2srgoKCbNonOjoaxcXFnCMzjym/noiLiwuSkpKQk5OjSCzGvh9b/XqANh6UAMBg\nMCAnJ8fkQ90999zDfaKGKd9UQUEB4uPjre6r0+lUy4hb4sMPP0RHRwfuvfde+Pj44Mknn8QjjzzC\n5bus+c4yMzORkpJi9nNHyex9/PHHiIuLG1gthwVSPXvnzp1Damoq1q5di23btqGnp4dZTIR8HNqz\nB/TfHJS8qFmaiSuiBbFXWFiI8PBwm7YdPXo0goODUVpayjmqwbS0tKC5udli5iwhIUG1m39ZWRnC\nw8OtNlQW0XJmD+hvy6GU2DPGHrEXEhKCzs5ONDU1cY7KMiUlJfD39ze5rvTNN9+MQ4cOKZ45s1Xs\nAf0PcFpaurGxsRGPPvootm/fPpCZvPfee6HX67F//37F47GW2UtJScG5c+c0vc5wb28vnnnmGTz5\n5JO47bbbcPjwYdViOXnyJK699lps3boVW7duRXx8PL744gvV4pFCb28v9u7dq7mMOCscXuwp/QRr\nrYQL9Is9tbMTlZWVWLNmjc3bqxGzWGq2JKbUKuulpqbaVcIF1Bd7ljJ7ABTN7Bn7fuwRe1ppv2Lp\nPPf09ERycjLXEqQp35Q9Yk9rmb1t27Zh5cqVmDVr1sB7rq6ueOGFF/D73/+eeVXBku9MEASrmT1v\nb29ERESo8nBkKx9++CEiIiLwhz/8AZ988gnWrl3LZGKOvZ69zs5OrFixAn//+99xyy23AAAefPBB\nvPLKK7JjUZLPP/8cd955J2JiYrB48WK88cYbIyo76WLrhp2dndDpdIOWtFKCdevWDZho/fz8kJKS\nMvBjTE9PR11d3cBFTUw/G3/O+vW3336LGTNmWNx+8uTJKCgoUCQeU6/nzJmDjo4OnDlzBjqdzqb9\nx40bh3379sHFxUWxeHft2jWojG9q+9ra2oEbv9LjeeDAgUFtGaxtr9frodfrIQgCdDqdovF2d3dD\nr9ejrKxsQKQM3b63txdHjhyx+f+H1eu8vDysWbPG5u1FgX/58mVF4jP1+sKFC/D29kZ6errJz1NS\nUrBz506bzy8Wr3/66adBYs/S9nFxcQPNoNUYv6GvT506hauuumpYPP7+/ujq6kJmZqZif+9x48bB\n2dkZOTk5yM3NNbt9ZGQkPv74Y2zZskXx8bL2uqenB0888QQef/xxAMDixYuxefNmLFu2DDk5OQgN\nDVUsnq6uLiQlJcHT03Pg77tq1Srcf//92L59+8Ca51oaP1Ovn3/+eWzcuBEPP/wwDhw4gM2bN6Os\nrAwvvPCCxf3F/1a77ZZVBDP09fUJO3fuFNauXSuEh4cLoaGhQkhIiBAeHi7cdNNNwhdffCEYDAZz\nu9tMUVGRMHnyZJOfWQhvgLa2NsHNzU3o7e2VHYstLFu2TNi9e7fFbQwGg+Dj4yPU19crEtNQCgsL\nhejoaCEtLc3mfV544QVh8+bN/IIywR//+Efhf/7nfyxu09fXJ4wePVpoa2tTKKp+0tLShC1btgiP\nPPKIXfsFBAQI1dXVnKIyT3Z2tjBu3DiL21y8eNHqNqww/u1FR0cL+fn5Nu/7+OOPC88++yyHqGxn\n3bp1wttvv23287fffltYt24dt+83de7GxsYKly5dsmn/ffv2CUuXLmUclXQiIiKEoqIik5/de++9\nwquvvsr0+yxd+3bv3i1cd911Vo/x4osvCg8++CDDqNjx1ltvCUuWLBn2/m233Sb87W9/k3Vse+4b\ngiAImzZtEp5//vlh77/66qvCzTffLCsWpcjMzBQiIyOFnp6egff27NkjzJo1y+5j2aJb1MBs/Sw1\nNRWnT5/Gww8/jMLCQlRWVqKqqgqFhYV4+OGH8eOPP+IXv/iFLKF56623Yv78+bh06RKioqKwY8cO\nu4/h4eGBgIAAxWYk2VLG1el0qpZyq6qqEBYWZtc+apTOrE3OAAAnJyeMHTtWlbG0t4wL9K81rEYp\n15pfD+gv1ZeWliraLLazsxPV1dWIiYmxeR8tzMi1dp4rvXpPT08PKioqbB7HuLg4zZRx6+vr0dra\najb2q666CkePHlUsnqysLIslXBGtzsjt6urCn/70J5NLuq1evRq7d+9WNJ59+/Zh+fLlw95fv349\nDh48qKq1xVb+93//F/fddx9cXK4UO5ctW4a6ujqcOnVKxcgYYk4FdnZ2WlWKtmwjBwvhDeKqq64S\n0tPTucYiCILQ2toquLu725RFXLt2rfDxxx9zj8kUX3zxhbBmzRq79snKyjKbYeXFlClThJ9++snq\ndqtXrxZ27typQETDv3fXrl127bNmzRpVYv3zn/9sU2Y2MTFROHfunAIR9ZOXlyfExcXZtc/hw4eF\nefPmcYrIOgaDQfDy8rKYmW9raxPc3d2F7u5uRWLKy8sTYmNjbd6+o6ND0YqHJdLS0oT58+eb/Twv\nL08IDw9nUimyhRtuuEH45z//aXW72tpawcfHR+jr61MgKtt5++23hWXLlpn8rKmpSfD29hYuX76s\nSCyXLl0SwsLCzP7tHnroIburI0pTX18v+Pn5mazIvPTSS8Ltt99u1/Fs1S1KYzazZ+zNa2xsxJkz\nZ5CRkTHwb+g2aqLUzLOcnByMHz/eap8roD9TptaM3MrKSoSGhtq1T1xcHAoKChRrhmkwGJCfn281\nsweol+mxp6GyiFqTNGzJ7AHKz8iVMoZqzsAG+jO6Xl5eCAgIMLuNh4cHYmJiFBtLeyZnAP0z7AMD\nA1FeXs4xKts4e/YspkyZYvbz+Ph49Pb2Knbe2JrZCwwMhJ+fn2YypCJpaWm4+eabTX7m6+uL2bNn\n48CBA4rEImb1zHWn+I//+A+mbWF4sGPHDqxcuXLYykNA/8IPe/bsQU1NjQqRscWs2BN58sknMWXK\nFGzatAmbN28e+KcllCpZXLx40WoJV0TN9itiGdfYQGoNb29v+Pj4KLbMW2lpKfz9/eHl5WV1WzV6\n7aWnp9u1VJqIWmLP2kxckaSkJFy8eJF7POJvT0opPCQkBO3t7WhubuYQmXWys7OtrlIB8F2ucei5\na6/YA/pn5Gqh/Yo1safT6ZiXcs1d+5qamlBbW2vT0n2ANku5mZmZFtvGyF1T3J77hrkSrsjUqVNx\n8eJFdHV1SY6HJ319fdi+fTs2btxo8vOAgACsXbuWe19NJbAq9j799FMUFBTg0KFDSEtLG/inJZQS\ne7b49UTU9OxJyewByvr2bPHriaiR6enu7kZzc7PNDZVF1BB7giBoNrMnReyp3X7F1vNcSd+eVLGn\nhazUuXPnLIo9oN+3d+zYMe6xnDlzBsnJyTZVZwDtNVdua2tDcXGxxd/nqlWrsHfvXu5N8tva2nDs\n2DFce+21Zrfx8PBAXFwcLly4wDUWqXz77bcICAiwuLTohg0b8Oabb2q656ItWBV7kyZN0nyTQS2K\nPTXLuGJmT5wibitaFntKZ/bi4+PtaqgsoobYq62tBXBlmSxLKJXZE397UsQe0D+OZWVljKOyjQsX\nLth0nvNcqnHouStF7GmhsXJfXx8uXLiAyZMnW9yOdWbP3LXP1hKuyPTp0zWV2Tt79iwmTJhgcenB\nuLg4BAUFSZ5YYOt94/vvv8esWbNMNh43RmtjaMw333yDm2++2eIiCSkpKYiNjVV84gtrrN7Jnnji\nCUybNg1Lly7FqlWrsGrVKqxevVqJ2GxGSbE3YcIEm7aNiIhAY2Mj2traOEc1nKqqKkmZPSXLPvaI\nvcjISDQ0NCg6llJKuIA6Yk/M6lm6YIkkJiYiNzdXMW+m1HGMiIhQTezZk9k7c+YMBEHgHpOjZvYK\nCwsRFBQEX19fi9tNmzYNeXl5aG1t5RqPtRLoULRWxrU1/tWrV8sq5drCvn37sGLFCqvbzZgxQ1Nj\naExGRgZmzpxpdbv169fj008/VSAiflgVe7/+9a/x2GOP4bHHHtOsZy80NBTNzc1cxUB3dzdKS0tt\n9no4OTkhNjZWlYutWMa1x3sB9IsqpQzdtnrMgP6xFCeQKMX+/fslZaRCQ0PR2NioaHsTe8bS19cX\nfn5+3IXQjsokAAAgAElEQVSU+NuTMkED6Bd7akwuEATBZrEXFhYGQRC4+FyNz11BEFBYWIi4uDi7\njqGFzN7Zs2eRnJxsdTs3NzdMmzYNJ06cYPK95q599mb2wsPDNbF8n4g9Yk9qJsqW+4YgCFb9eiJa\nzez19fXhzJkzNo3n8uXLsX//fodeUcOq2PPy8sKmTZuwaNEipKamIjU1VXZ/PdaIvdh4XthKSkoQ\nHh5uMX0+FDUmaRgMBtTU1CAkJMTufaOiohRbH9dWj5mI0pM0amtrJYkUJycnREREKLrOsL1jqVQp\nF5Bexo2MjFQls1dZWYlRo0bZVBLX6XRcS7kiVVVV8PDwsFouG4oWMnvWJmcYw9u319XVhUuXLlkt\nKRuj0+k0IZpFbBV7s2fPRm1tLbe/f3Z2NpycnGy67mh1neFLly4hODgYfn5+VrcNDQ1FQkKCov0g\nWWNV7C1YsACPP/44jh8/Pqz1ipbgfWHT6/UDy7bZihq+vfr6enh7e8PNzc1uz55SN9iuri7U1NTY\nVd5TepLGqFGjJJUfAeVLufZk9gBl1shNTU0dyIiYamlgDbUye/b4cgF+M3KNz93CwkK7S7hAf1aq\nsbER7e3tDCOzD3vFHqubqalrX3Z2NuLi4uDu7m7XsbQgmoH+xtrZ2dk2jaeTkxNWrFiBr776yu7v\nseW+kZ6ejkWLFtlkHfH29kZkZKRiD5i2kpGRMbD0qS2sWLECe/fu5RgRX6yKvYyMDJw4cQJPPPGE\nZsu4AH/fnhSxp8aMXCmrZ4iIYo+3B6m8vBxhYWE2z4gDlJ+kITUjBSgv9uzN7E2YMEGRC29FRYWk\nSS6Aepk9W9uuiCgxI1eKXw/ov+HHxMSoumanPWJv3rx5OHnyJLdZpGfPnsXUqVPt3k8rq5FcvHgR\n0dHRNrWrAvpF2w8//MAllmPHjuGqq66yeXstlnIzMjIwffp0m7dfuXLlyBZ76enpg1quaLH1CsD/\nhCwqKsLYsWPt2keNMq5x2xV7PXs+Pj5wdnbm7k8pLS1FdHS0Xfso3YojOzvbIcSewWBASUmJXQ8i\nSmT2xD6FUsdQzOwpMfnBGHsze7zKuMbnrlSxB6jr27t8+TKqqqps9jkHBgYiLCwM586dk/3dpq59\n9gp5Ea30K7R3csm8efNw/Phxu7/HlvvG8ePHMW/ePJuPqcVJGvaKvWnTpqG5uVnVhu9yMCv23n//\nfYs19u7ubklr2fJCi5k9Ncq4UmfiiiiRUSkpKbG7RKp0Zq+2tlZyGVfJ9XFra2vh4+OD0aNH27yP\nEmIPkD4TF+h/8HByclK8sbKtbVdEkpKSUFxczLVUKkfsqSlUzp8/j6SkpEHrjVrjqquukiRQbMFe\nIS+ilcyevWJv3Lhx6OjoYH49r6ysRFNTk13WEa1l9gwGg93jKac0rgXMir3Lly9j1qxZuPXWW7F1\n61Z8/PHH+Oijj7B161bceuutmDNnDjo6OpSM1SJKZPbsFXuxsbEoKytTdAaPcRnXXs8eoMwkjdLS\nUrtFgJLtV7q6utDe3i7JawYom9mTkj2LiIhAW1sb1/6Zqampkmfiiqjh27t06ZJdNzFXV1ckJibi\n/PnzTOMwPnflZvbUEir2lHBFUlJScPbsWdnfbc6zZ2vrLGMcNbOn0+kkZfes3TfErJ499oxp06Yh\nKyuLe6NnWykqKoKPj4/dTfNHpNjbuHEjMjIysGHDBvT09OCHH37A0aNH0dvbO/DZ/fffr2SsFhk7\ndiz0ej23so9er7e7jDtq1CiEhYWhuLiYS0ymkLp6hogSmT0pYk/J9isVFRUIDQ2V5DUDlBV75eXl\niIiIsGsfnU6nSHZPThkXUL7XXmdnJxobG+0+f6ZOncpEoJjDUTN7UsTe5MmTmQtnAGhvb0dFRYWk\ncYyNjUVxcbGqQsVgMCArK8susQdIL+Va4vjx45g/f75d+/j7+yMoKEiVNc5Ncfr0absmZ4hce+21\nOHHiBPd+kDyweDfT6XS4+uqr8dhjj2H79u3Yvn07Hn30UVx11VU2zcJREi8vL3h7e6Oqqor5sTs7\nO1FfXy9p4oPSvj3jzJ69nj1AObFnr2cPUM63V1VVBU9PT8n7R0VFoaSkRBG/mVRBJTZX5oVczx6g\nbN9HQPqEkvHjxzP/XYrnbmtrKy5fviz5AU7NmaRyxJ7cc2fotS83Nxfjxo2zq6Qs4uHhgYCAAFRU\nVMiKSQ5iJsqWlkDGSBF71u4bx44ds8uvJ6Il3569fj0Rb29vzJ07FwcPHuQQFV+kpS40Cq+SRXFx\nMaKiouyaPSoSFRWlaHZCbmZPiTKuFM8e0J8xU6J/XVVVFQICAiTv7+XlBQ8PD9TV1TGMyjRSMnvA\nlWwFT1hk9pQUe1Lj5fkQUlBQgLi4OMkP1+IEDaUnugD9JXF7ZokDQFBQENzc3Jj/3aX69UTULuVK\nyeoBwKxZs3D27FlmTd67urqQlZVlcS1Zc2jJtydV7AGOOyuXxJ4NSCnhioSHhyv6RGg8QUOKZ0+r\nZVxAuXYclZWVNnX9t4RSpVypAiUmJoar2EtNTWWS2VPyQUlqvDwmD4nnrtQeeyJ+fn5wdnZGfX09\no8hsQyyJS6mGsCjlDr32Xbx4UZbYU3uShr1+PRFPT08kJSXZJbIs3TcyMzORmJhoc/sXY7Qi9gRB\nkCX2VqxYgX379im25CQrSOzZgJSZuCJKZyfk9NkD+N9gW1tb0d3dLSlzptTNX+6MZkA5sSc1s8db\n7HV3d6O+vl7WODpKZk+0avDInklZJm0oarRfKSsrk9xjkYdvz9Eze1LFHsDWtye1hAtcEXtqZJmN\nKS0thaurq+T7ZHx8PPz8/DQhXO3B6pnY2dmJjz76CM8//zyeeeYZPPPMM3j22WeViM1ueIk9KT32\nRJS8YXV0dKCzs3Ng+Rcpnj2xjMvrhBSzelLKUkqKvZaWFlnH+Lln9r744guEhoZKsj6IOEpmz8fH\nB56enkz9wuK5K+faI6LkMogiJSUlkny5QL/Yk9trb+i1T67YUzuzl52dbdcyb8bMnz/frmXoLN03\njh07ZvfkDJGgoCCMGjWKi6/eHuxdOcMUjjgr16rYW7NmDXbv3g1XV1d4eXnBy8tLlnmdJ7zMyI6S\n2RMzUnImz/Dubya1hAso53+srKyU5dkD+oUK7xusIAgoKyuTlNmLjo5GWVkZt1KE1LWFjXGUzB7A\nz7fHIrOnRgsbOWIvOTmZaWavq6sLxcXFSEhIkHwMNSe69Pb2oqKiQvJ4ipk9uQ/wgiBImolrjNLN\n8U0hp4Qr4oi+PatTk8rLy/Htt98qEYtseDWzldJjT0TJC+3QyRlSPHvAlYyKLQtE24scsSf6Hw0G\ng+S2KLZQVVWFJUuWyDpGeHg413YcANDS0gKdTgcfHx+793V3d4evry+qqqoQHh7OPLagoCDZYi8o\nKAgtLS3o7Oy0q2m0VFiIvQULFjCJRTx3WWT2HE3sTZw4ERcvXkRfX5/kzLDxtS8vLw+xsbEYNWqU\npGMB6q5EUlpaipCQEMnxx8bGDqy0ExMTY3V7c/eNkpIS9PX1Sb4XAuzPEylkZGTg7rvvlnWMq666\nCvn5+aisrJRlm1ISq3fM+fPnc79psSIiIgKVlZUWV/6QgpwJGiEhIWhsbER3dzfTmEzBwmsG8C37\nSG27AgCjR4+Gr68vamtrGUc1GBbjqMTEHNGvJzWTy3PdVLmTM4D+3ophYWGKTXDSWmbPYDCguLhY\n1s0VUEfsyTnPvb29ERISwiyTJreEC/Sfz/X19aosJCCnsgT0t1CbP3++bN+eWMKVUzlSeiUkU0hd\nNs8YV1dXLF26FF9//TWjqPhjVewdOXIEM2bMwPjx45GcnIzk5GS7eycpxahRoxAUFITKykpmx2xv\nb0dLSwtCQkIk7e/s7Izg4GBFfApDJ2dI8ewBfL1SUtuuiPAujxoMBlRXV8tuOBwWFsb0d2gKuYKK\nZ/uVEydOyBZ7gHKNlXt6elBXVydZ5LMWe+np6aiqqhrwA8rB0TJ7QH8pV45vz/jax0LsOTs7c304\nsgSL7O68efNs9u2Zu2/ImZwhonYZt7u7GxUVFbIfoIB+354jlXKtir2vv/4aeXl5OHDgAPbu3Ys9\ne/Zg9+7dSsQmCdbGeL1ej+joaFllw/DwcEUutnJ77InwFFRyyrgAf9N+Q0MDvL29ZZV8AGUze1Lh\nOUmjpqZG1t9ZRKnGypWVlQgODpbUdBfgcxNjcZMHHFPssZyRy0LsAepN0mDxO7B3koYpDh8+jF/8\n4heyjqG22CsuLkZERARcXV1lH2vZsmX47rvv0NXVxSAy/lhVMLGxsWhqasLu3buxZ88eNDc3M1HF\nvGDt25NTwhVR6mI7NLMn1bPHcyKE1sWeWMKVOnYivr6+6O3txeXLl9kEZgK5mT2eYq+rq8uhMnty\nx3LcuHHIy8tjNos9NTUVRUVFsidnAFeuP0q1vBAEQXYGX67YMz5/pa6JOxS1JmmwuAfNnDkTOTk5\nNi3zZeraV19fD71eL3tiA+vzxF7y8/Nl9a00JigoCBMnTsThw4eZHI83VsXetm3b8Ktf/Qq1tbWo\nrq7Gr371K7z22mtKxCaJ6OhopjcwOZMzRJQSeywzezxusIIgyBZ7vGfksjLc6nQ67qVcLWf2WHj2\nAOUye3Lj9ff3h6urK1M/aWFhIZPMno+PD3Q6nex2QrbS0NAAV1dXSROHROSWcUV6enqQn5+PxMRE\n2cdSa5IGi3uQm5sbpk+fLtm3d/jwYcyfP19y5lvE398fbm5uqKmpkXUcqRQUFGDcuHHMjrdy5UqH\nacFiVey9++67OHnyJJ599lk899xzOHHiBN555x0lYpMEjzKuI2X2jMWeHM8ejzJufX09Ro8eLan7\nuohSmT2pY2cM71KuVjN7PT09qK6uZvLgodS5w0KcsixRpaenMyvjAsqWcuWWcIH+tZv1er3kZb7E\n87egoACRkZFwd3eXFQ+gXmaP1e/gmmuuwZEjR6xuZ+ral56eLrvaIaJmKZdlZg+44ttTu1G0Ldhk\nRDP2q/FsecECHmLPkTJ7LLJSvBory83qAcqJPRbwnknKKrPH+u9cVVUFPz8/Jr4YpRorsxB7CQkJ\nTG9iP2exN2rUKMTFxSE3N1fWcVj59QB1MnudnZ2oq6uTdZ6LLFiwQHLJkbXYU2tGLuvM3tSpU9Hd\n3Y0zZ84wOyYvrCq39evXY86cOXj66afx1FNPYe7cubjzzjuViE0SrMWeo5RxDQYDampqBs0alnpy\n8ir7sBJ7PGfjioKZxYUtPDxc02VcX19fuLi4oKGhgWFU/X9nVk/PP9fMnujZc0SxJ6ftijFyVtIQ\nz1+5a+IaI2b2lMziiN5HOSvRiMybNw+nT5+2OqFg6LWvvr4eRUVFsv16IqwfiuwhPz+fqdjT6XS4\n//778fLLLzM7Ji+sir3f/e532LFjB/z9/TFmzBi8//77eOihh5SITRJanaDBe2ZmfX09fH19Zc8i\nBfp/wDxEFYsnft5mc5aZPZ5l3K6uLjQ1NSE4OFjWcXiUcsvKypjMxAX6x7Cqqgp9fX1MjmcOrYm9\nnp4eVFVVMRtHR8vsAWxW0jh//jwzsefv7w8XFxfU19czOZ4tsEg2iPj4+CApKQk//vijXfuJfj0W\nmXpAvTJuX18fs0lPxtx3333Yu3ev4ksS2otZsSdmdRoaGjB27Fj86le/wn/9138hJiaGeSaAJX5+\nfujr62Oy3Fdrayva29sRFBQk6zhi6xWeT4RVVVXDegHK8Z3xmAjBIrPn4eEBT09P1NXVMYpqMOIk\nFxaePZ4TNCoqKhAWFibbVsFL7LFi1KhR8Pf3527o1prY+9e//oXw8HBmN1hHFHtyMnvi+XvmzBlM\nnTpVdiwiSrdfYZndBfpLudZ8e0OvfSxLuIB6Zdzy8nIEBgYy8W8a4+fnh3Xr1mHbtm1Mj8sas3eK\nW2+9FQAwffp0zJgxAzNnzsTMmTMxY8YM2YsI80Sn0zEr5Yp+PTkdwwH+680CQHV1teTGz6bg4ZVi\nIfYAvjNyh7avkQPPzB6r2a68xJ7cByRjeLdf6evrY7JsHEuxV1lZyfQmr1QLG4Cd2JsyZYqsGbkd\nHR3Q6/VM2q6IKN1YmUVlyRhbxN5QDh06xFTsiWVcpSc1sJ6cYcyDDz6IHTt2cL3Hy8Ws2BOnE+v1\nehQVFQ37p2VYiz0W8H6yNiX25JygWi3jAnxN+6z67AF8xZ5cv54IL7Ent/mqMbzbr1RXVyMgIEC2\nBWLMmDEwGAxMynw+Pj7MxZ6SmT0WD3Vij1cplaTU1FScP38eiYmJTKwtIrxbPw2FZRkX6Bd7x44d\ns2iLML72NTQ0oLCwkGmCJyAgAM7OztyqM+ZgPTnDmOjoaCxbtgxvvfUWl+OzwGoNaPHixTa9pyVY\nib3CwkJm9X01xJ4ctFrGBfiJvY6ODrS3t8Pf35/J8XiWcbWe2WMRmwjvrBSreHU6HbPsHuvynVJi\nr6enBzU1NbKzpEB/5wc5/faysrKQkpIiOw5jeK4bbgrWv4OgoCCEhYXZvN49a7+eiBq+PZ6ZPQD4\n/e9/r+lSrlmx19HRgfr6etTW1qKhoWHgn16vV3zpHXthNUnD0cWeHN8Z68yeWCpjkY3iJfbE3nA6\nnY6JZ4/nKhpazuyVlpYy/fvwzuyxFKesbmInT55kaiQPCQlBfX09enp6mB3TFOXl5QgJCWEmDqZM\nmSKprUV6ejqysrKY+vUA5VoBibAWe0B/vz1LLViMr33p6elMs/QiCQkJivv2WM/EHcrUqVMxadIk\nbseXi1mx99Zbb2HmzJnIzc0d8OnNmDEDq1evxsaNG5WM0W5YraLBcuYO7xm5WvfsVVVVMSmVAfza\nr7CciQvwXUVDq5m93t5eVFdXY8yYMcyOyXslEi2KPdaePRcXFwQHB3MdR4Bd2xWRqVOn2pyFGsqZ\nM2ccOrN3+fJltLW1Mb2uA/b59lhPzhBRI7PHs4wrsmPHDq7Hl4NZsffggw+iqKgIL7300iCv3tmz\nZx1C7FFmT55nj7U4ZeXXA/j5ZoyXm2N1gePl22OV2QsODkZ7ezva2toYRIUBobdkyRImxwOA0NBQ\nVFVVMTveULQo9urq6phndJQo5bI8zwHpmb1rrrkGZ8+edejMnl6vR0xMjOwJgkMRmyubmyAhXvty\nc3NRVVWFmTNnMv1+QHmxJwgC9zIuACbXZF5Y9ext2rQJ58+fx7/+9S988MEHA/+0DAuxJwgCs7Up\ngSvtV3jBOrPn7++Prq4uZiKAlV8P4HfBZZ3ZA/iJPZY+M5brSbP26wGOldlj0TCWV0bHEcVecnIy\nsrOz0dvba9d+hYWFCAgIYOa/FQkPD0dNTY3d8UiBRwkX6M/m+/j4WO239+677+KOO+5g7tcDlG+/\nUlNTA3d3d/j6+ir2nVrDqth7+umn8cADD2Djxo1IS0vDI488gt27dysRm2QiIiJQVVUl64Ssrq6G\np6cnvL29mcXkSJ49nU7HVKCyFHuiYZ/HMl9i2xUWnj2Aj1Bh1SpEhGUpVxROrMYP+Pll9oqKihAU\nFMQ8o+OIYs/b2xthYWF2j+lHH33EvIQLAK6urggKCuJeDgfYt10x5t5778X27dtNfpaeno7u7m58\n8MEHuOuuu7h8v+jZU6r9ihJZPa1jVex9/vnnOHjwIMLCwrBjxw6cOXMGTU1NTL78m2++QVJSEhIS\nEvDiiy8yOSbQf0IGBwfLyqiwLOECfC+0BoMBtbW1sldTGArLmFm1YwAALy8vuLm5obGxkcnxRIzL\nuKzgkdmrqamBv78/s5YSLMVeaWkp88xecHAwGhoauGVTWIq94OBgdHZ2yrpGFhUVMRPyxigl9lid\n5yJSSrn5+fnMS7givJdsFGHddsWYu+66C//+97/Ntj/Zs2cPJkyYgPHjx3P5ftHTq9QCDbwnZzgC\nVsWeu7s7nJ2d4eLigubmZgQHBzP5off19WHjxo345ptvkJ2djU8++QQXL16UfVwRuTcw1suqhIaG\noq6ujstsuMbGRnh6esLNzW3Q+3J9Zyx9e6yN2zxKucZlXC179lj59UR4ZPZYmrqdnZ0xZswYLqto\nGAwGpuPJov1KUVERF5+UI2b2AGmTNJqamrhk9gDlJmnwKuMC/WLrhhtuwHvvvTfss9TUVLzzzju4\n++67uXw3cOU8UaqUW1BQQJk9axvMnDkTjY2NuOeeezBz5kxMmzYN8+fPl/3Fp06dwrhx4xAbGwtX\nV1fccsst+Pe//y37uCJyfXusM3suLi4ICgpCdXU1s2OKsPbribC8ObAs4wJ8nq5Zrp4hwqOMy9oX\nx1rssc7sAPx8e3V1dfD29ma6hBILscfjJu+oYk9KZo9Hjz0RpSZp8BR7ALBx40Zs3759WIPl4uJi\n/Pjjj7jpppu4fTfQf54UFBRw/Q4RyuzZIPbefPNN+Pv747777sP+/fvxj3/8g8n04vLy8kE3Bda9\ntLQm9gB+F1tzYk+ub0rLYo/HjFzjMi4rz5kjZPZYzV4H+Hj2AH5ir6ysjPkMOrlir6CgAO3t7Qwj\n6oe32GtubkZfXx/zSRH2Zvbq6urQ2NjIrQSqVGaPp2cP6F8KNTw8HHv37h30/lNPPYVbb72V+Rqy\nQ4mNjWXe49McSrRd0Tou1jZYtWoVbr31VqxZs4bpD89W8/G6desGTlo/Pz+kpKQMlIjEG4qp19HR\n0fj2228H9QmytP3Q14WFhUhOTpa8v6nXbm5u+PbbbzFnzhwmxxNfi2Jv6OdZWVmyjt/Y2DhwDDnx\nzZs3D42NjcjJycGlS5eYjGdkZCR++OEHJCYmMjmewWBAZWUlcnNzERMTI/t44uvLly8PiBRWf++h\ngkru8WJjY1FSUsLkePn5+YiMjGR2PPF1aGgoDh06BC8vL2bnY3p6Oo4dOzYg9ljFO27cOBw5ckTy\n/rm5uVi7di3T8QP6H2BLSkogCMKgxuGsjv/FF18gICBg4NrO6vjXXHMNmpqasHv3bvj4+Fjdvq+v\nD/Hx8Th06BDT/z/xdWRkJI4fP858/IxfNzY2oru7G2fOnMHChQuZH198vWjRIrzxxhtYs2YN0tPT\n0dfXh3379uHAgQNc//8AoLu7G8eOHYMIz+/Lz89HVVUV0hnez8XX4n8ruWayJAQrpKWlCffdd58Q\nHR0t3HjjjcJnn30mdHR0WNvNKsePHxeuu+66gddbtmwRXnjhhUHb2BCeWXbv3i0sX75c8v6RkZFC\nUVGR5P1N8dvf/lZ47bXXmB5TEATh1VdfFTZu3Mj8uEeOHBHmzp0r+zj5+flCbGwsg4iu8N577wl3\n3HEHs+PV1dUJ/v7+zI4nYjAYBHd3d6GlpYXZMW+//XZhx44dzI7X3d0tuLq6Cj09PbKO09fXJ7i6\nujK5PgzliSeeEJ599lnmx/3b3/4m3H333UyPmZ6eLlx11VWS9u3q6hLc3NyErq4upjGJeHt7Cw0N\nDVyO/fXXXwvXXnstl2PPnz9fSE9Pt2nbrVu3Cg888ACXOARBEI4ePSrMnj2b2/EFQRAyMzOF5ORk\nrt8hCILQ2dkphISECH/5y1+EO++8U0hMTBQWLFjA/XsFof/3smTJEu7f09DQIHh5eQkGg4H7dwmC\nPN3CE6tl3NTUVLz55psoKCjAfffdh3/9619MZn3OnDkTeXl50Ov16O7uxqefforVq1fLPq6IHB9S\nV1cXampqmM8qVLqMKxdWEzRYl3AB9k2feczEBa60sGFZgiwvL2f622Qxex3o/x36+flh9OjRjCK7\nAq/2K+Xl5cxnvsop4xYWFiIqKorZTOuh8Czlsv5dGmOPb4+nXw/g19TdGB6zmk3h5uaG5557DllZ\nWZgxYwb++c9/4vvvv+f+vUB/RUGJbJg42ZJ1KyNHw6rYA/rXyd25cyf+9re/4ccff8Qdd9wh+4td\nXFzwxhtv4LrrrsPEiRPxn//5n5gwYYLs44qIjWIFCX18iouLERUVBRcXq1Vuu+C1ZFp1dbVJAW6c\nZpZCeHg4qqqqYDAYZB2Hx4WLtRduaENluWNnDGu/GQ+fGQvfnvHEEZbjB/Dz7LH2PwL9sba2tqK1\ntdXufXNycpCUlMR8/ER4ij0ev0sRe3x7WVlZXJseh4WFoba2lus6w6y7F1jinnvuwUcffYT7778f\nKSkp+OGHHxT5XnENe7n3F2vw9j46ClbF3s0334ykpCR8//332LhxIwoKCvD6668z+fJly5YhNzcX\n+fn5ePzxx5kcU8TX1xc6nU5SvysekzMAfjcsXpk9Nzc3+Pj4oLa2VtZxeFy4WK9IwmMmrghLYSoI\nApcMSnR0tGzTOa+ZuAC/zF5FRQVzgeLk5IT4+HhJ2b3c3FwkJiYyjccY3pk9nmLPlsxec3Mz1/50\ngDLrDPOohmgNd3d3+Pn5cW2YDvSLPZ6/B0fBqti78847UVhYiLfeegsLFy6Ek5NNyUDV0el0iIuL\nQ2Fhod37slwmzZjQ0FAuF4iamhqTYk80ksqBxc2Bx4UrICAAHR0dzGYtDi3jshg7EZZl3ObmZjg5\nOTFb2UWEdWaP5fgBfDN7PBoYSy3l5uTkDJp0xBpHFXuTJ0+2adm0/fv3Y8GCBVi6dCmXOER4N1ZW\nqoxrCl6/PVMoUcrlLf4dBavKbcGCBdiyZQvuueceAEBeXt6wqdpaJT4+XlIfH56ZPR5PMbwyewCb\nDBoPscfaC8fzRhUWFsYss8crTtZijzViZk+KLcMSPDJ7gHSxl5ubi6SkJObxiDiqZ8/b2xvx8fE4\ndeqUxe327duHFStWcInBGN7tV5Qs46qJEmKPMnv9WBV769evx6hRowamSIeHh+MPf/gD98BYIFXs\nsV49QyQwMBAtLS3o7u5mdkxBELj12QPY+Ax5PaWyLOVWVFQMyvCw9EyxLOPyElRRUVGyxZ7xUmms\nPWeenp5wdXVFc3Mzs2N2dnaitbUVgYGBzI4pIkXsCYIwkNnj5dljPVnIGJ6ePQBYs2YNvvzyS7Of\nGwllKxMAACAASURBVAwGfP3111i+fDm38RPh3VhZzTIu77EzRimxR549G8ReQUEBHn300YHZYZ6e\nntyDYoXWMntOTk4IDg5mmt1raWmBq6srPDw8mB3TGFZlXB5PqSwnvPDO7Gk9A6n1zB7A3rcnlu55\nWFOkiL26ujoIgoCgoCDm8YjwKod3dnaipaWFa+zXX389du3aZTa7e/r0aQQEBChyY+eZ2evr6+OW\ncdYavMWeIAjQ6/UDvVN/zli9yrm5uaGjo2PgdUFBwbA1WLWKFLEnCAI3sQewv9haKuFqwbPX0tKC\nnp4e5l31AbYZs6GZPdaePa1n9rTu2QPYnzs8Bb4UsSeWcHU6HTffFI8VXYAr5w9PT/f06dPR1dVl\ndg114xIub98Zz/YrVVVVGDNmjGr32ZHk2WtoaICLiwv8/Py4fYejYPXMfPrpp/HLX/4SZWVluO22\n27Bo0SK8+OKLSsQmGylir6GhATqdjos4AdhP0uDp1wPkiz2xHMGjxxEr/5EgCKisrORi1AfY3mB5\nCZSAgAB0d3dLahcC9JfQeHq2AD5ij9ffPDIyEg0NDWhra7N5H7GEy5PQ0FDU1NQMWw9VLrxLuEC/\nT/f66683W8r96quvsHz5cq4xiPCcoPFzmIkrwlvs0eSMK1gVe0uXLsXOnTuxY8cO3HbbbTh9+vTA\n8i1aJyoqCjU1Nejq6rJ5H55ZPYD9JA1LYo+F90KuL47nhYuViKqvr4eHh8egtSBZ+lZ8fHzQ19cn\nWUgZwyuzp9PpZLVfqaurg7e398AY8vD9sC7j8iyVOTk5YezYsXY9bBpPzuDlmxo1ahR8fX1RV1fH\n9Lg8s6TGmBN71dXVuHTpEq6++moA/H1nPDN7as7EBZT17PHutUd+vSuYFXunT59GRkYGMjIyUFJS\ngrCwMISFhaGkpAQZGRlKxigZFxcXREVFoaioyOZ9eE3OEFGyjMsCFpk9XrPKWIm9oSVc1rCcOczz\npiqnlMvbrwc4VhkXsL+Uy7vHngiPUq5SYm/BggUoKCgYJrS++eYbXHvttXB1deUeA9D/4FFXV8d0\nsp3Izymzx7vXHs3EvYLZJSI2b95ssfSWlpbGJSDWiKVcW9sZKJHZYymWeXv2AgMD0dbWho6OjkGZ\nL1vheeFiVcY1daNi7VsRxd748eNlHYenqJIj9oxn4gJ8fD+hoaE4f/48s+NVVFRwXVbLXrFnXMbl\n6ZsSf4vTpk1jdkyeDbWNcXV1xYoVK7B7927cf//9A+9/9dVXg1qu8Padubi4ICQkBBUVFczFhNpt\nV5T07AFXSrk8Hrj1er3sa+5IwWxmLz09HWlpaWb/OQr2+vZ4NVQWcbTMnpiVkpoJ4FmSEPvXye29\nxjuzB7DptcezVQggr/0KZfaGk5CQYLPY6+7uRklJCeLj47nFI8Ky76OIUpk9YHgpt6enBwcOHMCy\nZcsU+X4RXu1X1C7jKg1P3x559q5g1bPX1taG5557ziGbKgP2i73s7Gyma/QORUmxx8p7ISeDxvMp\n1cvLC66urpKWxDPG1I2KtW+FRelMnFDAa8YjyzKuI3j2tFTGLSwsRFRU1MAMTJ6+KUcu4wLAdddd\nhxMnTqCiogK7d+/GunXrkJCQwG1ta3Pwar+idhlXSc8ewFfskWfvCiO6qTJgn9gTBAHnzp1DcnIy\nt3hY37B4Z/YA+WKP54WLRa89pTJ7ckU+7xuqnAkaSpTxWD4oCYLA/e9uj9hTYiauCI9ee7xnYhvj\n6emJhQsXIi4uDq+88gpmzpyJXbt2KfLdxvCapKF2GVdpeIk96rE3mBHdVBmwT+yVlJTA09MTY8aM\n4RaP2PqA1ewj3p49QPqMXEEQuIs9VhkzJTx7cuPkXSqVk9krLi4e9Hfm4fsZM2YMWltb7Zpdb46m\npia4urrCy8uLQWSmiYqKQm1trU3tV4Yuk8bbs8cys2cwGLi2LjLFjh07UFlZibS0NDz00EPcz19T\n8Gi/0tXVhYaGBu4P8JZQy7PHmrq6OowePRo+Pj7Mj+2IjOimygAQFxcHvV5vk7jindUD+lsf+Pj4\nMGt9oFRmT8rNoa6uDu7u7lwfEFgsmaZEZo/FbFzemT3RgyTlQaSgoADjxo3jENUVxBVoqqurZR9L\nib+5i4sLkpOTbZqQpWRmj/WSaTU1NfD19VX0vhAQEMCtF6qt8MjslZWVITw8HM7OzkyPq2V4iT3y\n6w1mRDdVBgAPDw/4+/vbJAiUEHsAuzLK5cuXIQiC2eyE2p49JcoRrMq4vD17LEzxvEtlo0ePhr+/\nv91iqq2tDfX19YPGkJfvh5UNQimP2dy5c3HixAmr2w1tu8LTN8V6goaSJVxbUcJ3xiOzV1JSonoJ\nV2nPHq9ee+TXG8yIbqosYmsp19HEXnV1NYKDg7msTmGMHLHH28cltyTV09OD+vp6BAcHM4xqOKzK\nuLwFipRSrjiDXYlsBKtzRymxN2fOHJw8edLiNoIgDCvj8oT1KhpK/C61CI/MntqTM9SAV6896rE3\nGJum9QUGBmLlypVYuXIlt7YPPNGa2GOVnbBWwmXlvZAq9pRoISC31151dTWCgoLg4jK45SRr34qP\njw8MBoOsVTSUyKBIab9iqoTLy/fDarlBJcq4gG2ZvYKCAri7uyMoKGjgPZ6+qVGjRsHPz4+ZlUTJ\nmbi2ooTvLDQ0FPX19UwbK2tB7Cnt2QP4lHJJ7A2G36rVGsIWsdfd3Y38/HyubVdEWGb2lDDyilkp\ne/vZOUJmj+f6qMbodDrZf3etZvby8/MV6Q8HsFtuUCmBEhcXh66uLosZoO+//x6LFi3inqE3hmUp\nV4tiTwmcnZ0RGhrKpLG7iBbKuGrAQ+yRZ28wJPb+P7m5uYiJiZG0SoS9sBJ7NTU1FsUeK++FOMnC\n3kyAEieb3Aka5jI8PHwrcoRpX18fqqurERYWxjiqwUhpv2Iqs8fL98Pq3OG5Lq4xOp3OailXFHvG\n8PZNsZyk8XP17AHsS7layOwp7dkDKLOnBGbFXkNDg8V/joQtYk+pEi7A7oZVVVU1qJEoT6ScjHl5\nedyXqgkLC5PlP1IyKyHn715dXY2AgICBFki8kJrZ4z0TV4TlBA2lWoVYKuUaDAaTYo83LNuv/Fw9\newD7SRpaEHtqEBsba9ca9tYQBAHFxcUk9owwuzbu9OnTLZYVWP5heKNFscfihlVRUWFxfUuW3gtx\nDGfNmmXT9oIgIC8vDwkJCcxiMIWrqyv8/f1RW1srSfiay+zx8K3IucEqlT1hVcbl5ftxtAkaQP8k\njeeee87kZxcuXICvr++w0h1v39RIL+Mq5TtjndnTQhlXLc/ezp07mR2vpqYGnp6eXPtoOhpmxR6v\n5UvUYMyYMTAYDGhoaEBAQIDJbc6dO4e77rpLkXhY3rCMF//mib3LzlVUVMDb21uRhpZiKVeK2Csv\nL1fs4iandKZU9sResdfd3Y2KigrFutSzOHd6e3tRW1urWOPa2bNnIyMjAz09PXB1dR30mRpZPaD/\nt3ju3Dkmx9JiGVcpIiMj7bouWqKlpQW9vb2q9w9Ug7i4OKYJJPLrDccmz15jYyNOnTqFw4cPD/xz\nJHQ6nVWxomRmT5xRaO+Eh6FYm1HI0nsRFxeHwsJCm7e/dOkS9xKuiJxee0p69uRkU5TKngQFBaGt\nrc3mWcN6vR6RkZHDyss8PXvV1dWy2oZUV1cjMDBwmPDihZi5O3/+/LDPzIk93r4pVpm9lpYW9PX1\nwdfXl0FU7HBEz55YwlVyoo4p1PDsxcTEoLS0lFk7IPLrDceq2HvnnXdwzTXXYOnSpXjqqadw3XXX\n4emnn1YgNLbEx8ebXaeyqakJ9fX1iIuLUyQWb29v6HQ6WW04AOXaRwD2Z/aUKOGKyCmPKmXUB+TF\nqVRmz8nJCePHj0dubq5N2xcUFCg2Exfobxvi7++PmpoaycdQ8m8uYsq319vbi8OHD6tSNmM1QUN8\nCFFboKhFVFQUM8+eFkq4auHm5obg4GBmwrmoqIgaKg/Bqtjbtm0bTp06hdjYWKSlpSEzM1NzT3G2\nMHfuXLNPLOfPn8fEiRPh5KTc5GS55aienh6rayjy8OzZipKZPTkzcs0Z9Xl59qT+zZU0GyclJSEn\nJ8embc1NzuApYCIiImTdFJScnCFiakZuZmYmIiMjTZ7DvAUgqwkaWvTrAcr5zlhO0NDK5Aw1Hj4A\nYOzYsXZVjyyh5KQxR8Gquhk9evRAO5LOzk4kJSXZ/NSvJVavXo09e/aYXJJFyRKuiNxJGlVVVQgK\nClJsDcXIyEjU1NSgs7PTpu2VzOxJLeO2tbWhq6tLMY+MnNKZkk+qEyZMwMWLF23aVo2LqtxG2moI\nFFOZve+++04Vvx4AhISEMFlFo6ys7Gfr1wP6x7GxsRFdXV2yj1VcXPyzzewBbH17JPaGY1XsRUVF\nobGxEddffz2WLFmC1atXO2QtPCEhAT4+Pjh9+vSwz86dO4cpU6YoGo/czJ4tJVyW3gsXFxdER0fb\nPHFH6cyeFBEljqGpEhQP34qcVTTEJcmUICkpyWaxZ66My9P3w0LsKZ3ZmzRpEsrLy9HY2Djw3vff\nf4/Fixeb3J63b4rVKhpazewp5TtzdnZGWFgYk8bKWik9quHZAyizxxurYm/Xrl3w9/fH008/jeee\new533303vvzySyViY46Y3RuKGpk9ucs+qeE7srWU29vbi6KiIsW8XFLLuEqPodRVNC5fvozLly8r\n1lNxwoQJssu4PJEr9kpKShSbPSzi4uKCGTNmYOvWrdDr9ejq6sLx48dxzTXXKBqHMSxKuT9nn5kI\nq0kaWhF7asEqs9fe3o76+npNPoSoiV1NladMmYKrr74aly9fVjJGZqxevRq7d+8e9F5HR4dqZVw5\nYs+W7ARr74WtM3JLSkoQEhKiyGokgPSbv6Ux5OVbkXKDFdsIKGWCT0hIQFFREXp6eixu19fXB71e\nb3JiE2/Pnhyxp9ZMvb/85S8oKSnBnDlzkJCQgKSkJPj5+ZncVgnfFIs2Nlqd9aik74yVb0+v12tC\n7Dm6Z0+sgihlcXIUrDZVFgQBJSUlA76mxsZGxMTEOFRTZZF58+ahvLx80NPo66+/jsWLFw9ahFwJ\nwsLCcOHCBcn7KzkTV8TWzJ6SJVwACAwMRFtbG9rb2+Hh4WHzfmqMoZRJGkVFRYrNFAf6fboREREo\nLCxEYmKi2e3KysoQGBiomKgXkSv2iouLFc/sAf399j744AMYDAZkZWUp1vrFHCwye2qNpZZgkdnr\n6OhAY2Oj4tcjLcEqs0clXNOYzezp9XoUFRVhyZIl2Lt3L+rr61FfX4+vvvoKS5YsUTJGZjg7O2P5\n8uUDpdyGhga89NJL2LJli+KxyJ2gobRnD7Bd7Ck5OQPobxcSExNjdyNwS34jnn3i7L3BKunXE7HF\nt2fposrT9xMZGSlZ7HV3d6OmpkbVEo+TkxOmT59usZqghG9Kbq89cUkqLYo9JX1nLNqv6PV6REdH\nK9oRwhxqefZCQ0PR3NyMtrY2WcchsWcaq7+s48ePY/ny5QOvly1bhmPHjnENiierVq0aEHsvvPAC\nbrzxRovZC16wKOMqfcOytYyrdGYPkLZ2r1qZPSllXCUze4Btvr2CggJVLqpyMnulpaUIDw+Hi4vZ\nosbPBrm99mhJqn4iIyNlZ/ZoxYf+hyAp1/GhkNgzjVWxFx4ejj/96U8Dmb7nn3/eoY2P1113HY4d\nO4bs7Gy89957eOqpp1SJQ+5TtS1ChYdnr6ioyGT7GmPy8vIcQuxZahvB07Nn7w1Wy5k9c5NwePp+\nfH190dfXh5aWFrv31WomaihK+KbklnG16tcDlPWdscjsaWlyhlqePYCNb0/pRu+OglWx98knn6Cm\npgY33HADbrzxRtTU1OCTTz5RIjYueHt7Y968eVi1ahXuu+8+1TwSY8aMQUdHB9rb2yXtr0ZWysvL\nCz4+PlbFyqVLlxQt4wL9Fwl7/R4FBQWKZ8ykiHy1Mntyyrg80el0krN7WhYoSiO3uuAowpk3LCZo\naEnsqQkL3x5l9kxjVeyNGTMGr732GjIzM5GZmYlt27YhICBAidi4sXr1ajQ1NeGRRx5RLQbxhiUl\n/d/e3o6Ojg6rfwce3gtrpdyuri5UVFQofkO1N7PX2tqKlpYWhIWFmfycl2/F3syeIAiq3AjEVTQs\nrd9sqYzL2/cjR+w5gkBRwjclZ+UZQNvCWUnfWUhICJqammxuOG8KrczEBdTz7AHyM3vi/ccRznGl\nMSv2Nm7cCKDf4zb03+rVqxULkAd33XUXjhw5ovqyb1JncVlqBswba5M0CgsLER0drfhMw9jYWLue\nCAsLCxEXF6e4IdrezF5NTQ08PDzg7e3NMarhBAQEwN3d3WysgiCoWi6RKvaUXHZO64SFhaGurg7d\n3d2S9qfMXj9OTk6yhTNl9vqRm9nT6/WIiopSfaa7FjHrUv7HP/6BN954A5s3bx72maMvej169GhM\nnDhR7TAkp/9tLeHy8F5YE3tqTM4A+p8I7cnsWRMqvHwrxqto2CLg1PDriYjZPVMe3ZycHAQGBsLH\nx8fkvrx9PyO9jKuEb8rFxQVhYWEoKyuTZBPQ6/VYunQph8jko7TvTHxwl/rwoyWx58iePfLrmcds\nWkMsz6Smpg7794tf/ELWl3722WeYNGkSnJ2dkZGRIetYjozUWVxqrJ4hEhcXZ1HsKd12RSQoKAgd\nHR02L0Wmtt/M1r+7Gn49EUu+vf3796t6o5eT2aNs1BViYmJQXFwsaV/Kkl5BziSN5uZmdHd3IzAw\nkHFUjofovbZkH7EE+fXMY1bs1dbW4uWXX8bWrVuH/Xv55ZdlfWlycjJ27dql6lJBWkBuGdcaPLwX\n8fHxFp+81Mrs6XQ6u3x71p4AefpWEhISkJ+fb9O2WsjsmcKa2NOiZ6+3txeVlZVmZ2BrCaV8U9HR\n0ZLEniAImvY/Ku07k9N+RekVcqyhpmfP19cXbm5uqK2tlbQ/iT3zmBV7fX19aG1tHViX0/iflIXc\njUlKSlJFEGgNqWVcNRZyF7Els6fW35al2OPJ+PHjkZuba9O2WszsdXV14ciRI1i0aJEKUfUjpbFy\nWVkZQkJCMGrUKE5ROR4xMTEoKSmxe7+Ghga4uLiYXe7t54aczJ6WSrhaQI5vj8Seecx69kJDQ1Xr\nQfdzQU5mb/r06Va34+G9CA0NRXt7O1paWob5tQRBQE5OjiplXMC+SRrWLgo8fSvjx49HZmamTdsW\nFhbi1ltv5RaLJcxl9o4dO4aJEydanA2uRc+eI5VwlfJNxcTE4MSJE3bvp/WxVNp3FhkZiQMHDkja\nV0szcQF1PXvAFd/enDlz7N6XPHvm4TYVccmSJUhOTh72T1y9gpA3QUMtz55OpzPbfuXs2bPw8PBA\nVFSUCpHZPkmjq6sLlZWVA+sjK8348eNx6dIlm7ZVM7MXFRWFpqYmNDc3D3pfbb8e0P/QUVtbi56e\nHpv3cZTJGUoi1bNHYzkYOevjUmZvMFIze729vSguLlbteql1zGb2Dh48KOvAUp9yhrJu3bqBi4qf\nnx9SUlIGnjxEb4Gjvj5//jxaW1vR3t4ODw8Pm/cXy7jWtn/11Ve5jJdYym1qahr0+csvv4yZM2cO\neE+UHs+2tjacOnUKIua2DwsLQ1RUFI4ePWr2eMa+FdbximLP2vYHDx5EeXn5gHhWejwPHz6M8PBw\n5OTkYM6cOQOf79+/H9u2bbO4P8/xS01NhYuLC3x8fLBr1y7cfPPNNu2flpYGY9Q+/9UcP/F1VVXV\ngNizZ//i4mI4OzsjPT1dE+M19LVS4ye+joqKQkFBgaTxKCoqwqJFizQzfuJ7an3/2LFj8eOPP9q9\n/2effQYfHx+MHj1a0XjF/5a7zBt3BBVJTU0VfvrpJ7OfqxyeIsTHxwu5ubk2b28wGAR3d3ehtbXV\n6rZpaWkyIjPPH//4R+HBBx8c9v7kyZOFo0ePcvlOW/jxxx+FlJQUq9t99dVXwtKlSy1uw2vsBEEQ\n+vr6BA8PD6GlpcXidvn5+UJsbCy3OGzhoYceEu65556B1zU1NYKvr6/Q3d1tcT+e4ycyc+ZM4fjx\n4zZvv379euGdd97hGBE7lBg/QRCEtrY2wc3NTejr67Nrv02bNglbt27lFJV8lBo/EYPBIHh5eQmN\njY127ztp0iQhKyuLQ1TSUHrshrJ//35h0aJFkvZbuHAhh4jsQ6u6hVsZ1xK7du1CVFQU/l979x4U\n1Xm/AfxZkXBVEY1WQSXCAnJbCIKiUamXhJjo4C3GW6xGezGa6IymdtRWjZq2qbGdabWJ0cQS05oW\nL6NW7RRtvCSiIuANAQUVUFS8Al64+P7+8Lc0uLtwdpc9tz6fGWfcs3t2333mLPvd97zve44ePYrX\nXnsNr776qhLNUAV7Z3Hdu3cP7u7uki4+bv4F0tKmT5+OtLQ0VFdXN2wrKChARUUF+vbt65LXlELq\naVwp4zpclR3wdBFWo9GIwsLCJh+n5Excs1/96lfYu3dvQ0//v//9byQnJze7aKkr8zOzd9yelk49\nypEfAHh7e6Nt27a4fv26XftxzF5jBoMBISEhzX6mnyX+/wo5ajou5c7uWY6utdfUtbrJhWP2mjJq\n1CiUlJTg4cOHKC8vx549e5RohirYO9ZDiWviPuuFF15A//79sXnz5oZt27Ztw6hRo2S/IsX3+fv7\no66uruH0si1qmLElZdyekuP1zNq1a4dPP/0UM2fORGVlpSrG65k5UuypuUBRiiMzcrVUOMtFyg+4\nZ928eRMeHh6KX81JTXr06IHy8nI8fPjQrv2aunwjKVTs0X/ZO0nDnmVXvj+moKXNnj0bf/zjHxsW\nv9y6dStGjx7tsteTwmAwSOrdk9Kz58rsAGnFnhp69gAgJSUFycnJWLhwoeRiz9X5AfYVe/X19Sgr\nK1NsUo695MjPzJFJGmrv2ZMzP7PQ0FC7iz21zcQFlF1nDwDc3d1hNBptLuhuixp+xKsZiz2F2Xsa\nVw09ewAwZMgQ1NTU4ODBgygtLcWFCxecvrJKS5Cy/IoapudLWWtPDT17Zh9//DG2b98ODw8PxbMz\ns2etvatXr6Jjx47w8PBwcau0x95iz3zFhw4dOriwVdrjSM8eZ+JaFxMTg1OnTtm1z5kzZ9CrVy8X\ntUj7WOwpzN7FOO1ZdsXV487MvXvbt2/HiBEjVHHx6eYWVq6vr5dURLl63IqWevYAoH379khLS8OC\nBQskrfSvtjF7WjuFK+e4KXuLPfNl0tRyxQdrlBh3ZjQaJS+pZKbGYk/pMXvA06ts2VPs3b17F+Xl\n5QgLC3Nhq7SNxZ7C7O3ZU/LqGc966623kJGRgU8++UTxU7hmzZ3GLSsrQ4cOHeDt7S1fo6wwF3vC\nxjUg6+vrUVBQoNgC1dYMHjwYP/vZz5RuRgN7rjHM67jaZu8l07RWOMvF3LNn6zNtjRqLPTWwt2cv\nJycHJpMJbm5uLmyVtrHYU5grJ2i4euxF27ZtMWnSJBQXF2PYsGEufS2pmjuNK3XGlquz8/f3h7u7\nO27cuGH1/tOnTyMgIECzp8rkHLMn5ctVaxMK1DxmTwuFsxLjzjp27AgAuHXrluR91FjsKT1mD3ha\n7J0+fVry40+ePCnpqlL/y1jsKaxjx46oqqrCgwcPJD2+tLRUsatnWLNgwQL8/ve/h5eXl9JNAdB8\nz56aZmw1dSr30KFDeOmll2Rukba0adMGrVu3bnb2NcDeqKbYOxuXWVpnMBjsHreXl5fHU49WdO3a\nFXV1dZKXBGKx1zwWewozGAySxx4JIVBYWCi5WJFj7EX37t0xY8YMl7+OVOaePVu9PVInZ8iRXVhY\nmM1i7/Dhw5ou9uQa9yN1NrsWeqO+T85xU+3bt0d9fb3FZfFs0UKWSo07s6fYu337Nu7du6e6LNUw\nZs9gMNh1KpfFXvNY7KmA1EkaFRUVMBgMDacLyJKfnx/c3Nxw+/Ztq/eraXq+rZ49IYTmiz25RERE\n4MyZM80+rqioSHVfqmphMBjsOpXLnj3b7JmkkZubi5iYGEXXJlUzqcVedXU1Ll26hIiICBlapV08\nylRA6iSN/Px8hIWFSZ4Fp4axF0po6lSu1J49ObKzVeyZxxyqbSyPPeQ69kwmU7NfCPfu3cP169dV\nU+RLIfdnV+okjfr6euTl5al+iQul/vbZs9Zebm4uTCaTi1tkP7V8b0gt9nJzcxEZGamK1SDUjMWe\nCkjt2cvPz0doaKgMLdK2qKgoHD9+3GK7EEJVl9SxtdaeuVdPzUtbqEVMTAxyc3ObfMzJkyc5U68Z\nUnv2Lly4gE6dOvGKDzbYcxpXrcWeWkhdfoWncKVhsacC9vbsSaWGsRdKSE1NxbZt2yy2FxUVwcfH\nB+3bt2/2OeTILiQkBEVFRaivr2+0XQ+ncOU69kwmU7PFXlZWFuLj42VpT0uR+7MrdZKGVgoUpcfs\nSZkhrtYs1fK9ERkZifPnz6O2trbJx7HYk4bFngpILfYKCgo4c0uClJQUfPfdd7hz506j7X/9618x\nduxYhVplycvLC507d7boUTl8+DAGDBigUKu0pUePHqiursbNmzdtPkaLxZ7cpPbsqbVAUQs/Pz94\neXmhvLy8ycfV1tbi/PnziIqKkqll2uPj44Nu3bo1OwaSxZ40LPZUwFWncdUy9kJuPj4+GDx4MHbv\n3t2wTQiBtLQ0TJkyRdJzyJXds+P2KioqUFZWhujoaFle31Xkyk/KrL2TJ09qrtiT+7Ort2JPyb99\nUk7l5ufno1u3bvDx8ZGpVdKp6XujufX2Hj16hIKCAs3/vZQDiz0VkNKzV1dXh+LiYk0NMlfSqFGj\nsHXr1obbx48fhxACiYmJCrbKUnh4OI4cOdJw+8iRI0hKSuL4Mjs0NUnj/v37KC0tRXh4uMyt0hap\nEzS0UuwpSUqxxxylae6H3JkzZ2A0GuHp6Sljq7SJxZ4KPP/8880urHzp0iV06dLFrsWL1TL2CO4e\nBwAAFOJJREFUQgkjRoxARkZGQ6ZpaWmYPHmy5EkPcmU3b948bNy4Ef/4xz8A6OcUrpzHXlPj9rKz\nsxETE4PWrVvL1p6WIPdnt0uXLrh9+zYePXpk8zG3bt3CvXv3NDFLXMm/fVov9tT0vSGl156ncKVh\nsacCUhZWtndyxv86f39/JCYmYt++faitrcWWLVswefJkpZtloWfPnti1axdmzZqF/fv388oZDmhq\nRi7H60nj5uaGgICAJoeTcF04aaQsv6LmYk9NWOy1HH5qVaJHjx4oKiqyeb8jy66oaeyFEkaNGoVt\n27Zh3759MBqN6Nmzp+R95cwuLi4OX3/9NcaPH4/Tp08jISFBttd2FTnzi4qKsjlrT6vFnhKf3V69\nejX5xaqlAkXpMXvNTSrIyclRbZZq+t4ICgrCnTt3LCbbmbHYk47Fnkq8+OKLyMrKsnk/Z+LaLzU1\nFbt378bGjRslT8xQSnJyMtavX49x48bB29tb6eZoinnWnrU1C7Va7Clh4MCBOHjwoM37tVTsKSkk\nJAQXL17EkydPrN5fXl6O2tpaBAYGytwy7WnVqhWio6OtfjfW1NTg7NmzPCYlYrGnEomJicjMzLR5\nvyOncdU09kIJXbt2RWhoKHbv3o1x48bZta8S2aWmpuKLL76Q/XVdQe78rI3bq6ysRElJiSYvo6TE\n8Tdo0CB88803Nu/XUrGn5N8+X19f+Pn52RyWY85RrYumq+17Y9KkSfjDH/5gsX3Tpk3o168ffH19\nFWiV9rDYU4k+ffogMzPT5mKcvHqGY958802kpqaiQ4cOSjeFXMjajNycnBxERUVpbnKGUuLj43Hx\n4kWrp8y4Lpx9IiIikJ2dbfW+3NxcxMbGytwi7Xr77beRnZ2NEydONGx7/PgxPvjgAyxfvlzBlmkL\niz2V6N69O4QQVpdgqaysxL179+zu9lfT2AulvPvuu0hLS7N7P2bnHLnzszZJQ8uncJU4/tzd3dG3\nb18cOnTI4r7z58+jR48eqlwXzhqlP7+jRo3C119/bfU+tfeQKp3dszw9PbFw4UIsW7asYdv69esR\nExODpKQkBVumLSz2VMJgMCAxMRHHjh2zuK+goABGo5Gz4BxgMBjw3HPPKd0McjFrPXtaLvaUMmjQ\nIKvj9tQ8oUCNxo0bh127dqG6utriPrUXe2o0Y8YMZGdnIysrCw8ePMCqVavYq2cnVg8qYqvYc/QU\nrtrGXmgJs3OO3Pl1794dDx48aHTZNC0Xe0odf7bG7WmtQFH689upUyckJSVh586djbY/evQIFy9e\nVPU4UqWzs8bT0xM///nPsWzZMqxduxZJSUmchWsnFnsqYmuSBmfiEjXNfNk086ncsrIyXL58GZGR\nkQq3TFsSExORl5eH+/fvN9qutWJPDSZOnIivvvqq0bb169cjOTkZHh4eCrVKu2bOnImsrCysWLGi\n0SldkobFnookJCQgKysL9fX1jbY7uqCy2sZeaAmzc44S+ZlMJvzkJz9Bp06dEBkZialTp8Ld3V32\ndrQEpY4/Dw8P9O7du9El/IQQmiv21PD5TU1NxTfffINbt24BAG7fvo0PPvgAH330kcIta5oasrPG\n09MTK1aswKRJkzhRyAEs9lTE398fXbp0QV5eXqPtvHoGUfMWLlyIdevWITs7G3fu3MHatWuVbpIm\nPXsqNz8/H0+ePEFAQICCrdKeNm3aICUlBenp6QCA5cuXY+zYsSxUnDBt2jT86U9/UroZmmQQttb6\nUAGDwWBzKRK9mjx5MgYPHozp06cDePqruk2bNigrK0O7du0Ubh0R6d3+/fuxaNEifPfdd7hz5w76\n9euHd955B7Nnz1a6aZqzY8cOrFmzBp988gleeuklnDt3Ds8//7zSzSIXUmvdwp49lXl23N7Ro0fR\nvn17FnpEJIu+ffvi1KlTuHv3LsaOHYtXXnmFhZ6DUlJScPr0aUybNg0LFy5koUeKYbGnMt+fkVtZ\nWYkpU6bg448/dui51Dr2QguYnXOYn3OUzM/b2xuxsbEYOnQofHx8sHr1asXa4ii1HH8eHh4YPXo0\nbty4oZmCWS3ZUctisacysbGxyM/Px4MHDzB79mwkJyfbfakvIiJnDBkyBEIIfPXVV3Bzc1O6OZq2\nfPly7Nq1izNwSVEcs6dCCQkJSExMREZGBrKysjSzaj0R6cPjx48BgAUKkZ3UWrfwopEqlJiYiM8+\n+wxHjx5loUdEsmORR6QvPI2rQtOmTcOXX36JuLg4p56HYy8cx+ycw/ycw/ycw/wcx+z0iT17KtS7\nd2/07t1b6WYQERGRDnDMHhEREVELUGvdwtO4RERERDrGYk/HOPbCcczOOczPOczPOczPccxOn1js\nEREREemYImP2FixYgF27duG5555DcHAwPv/8c6uXA1PruW8iIiKiZ6m1blGkZ+/ll1/G2bNnkZub\ni9DQUHz44YdKNIOIiIhI9xQp9oYNG4ZWrZ6+dJ8+fVBaWqpEM3SPYy8cx+ycw/ycw/ycw/wcx+z0\nSfExexs3bsTw4cOVbgYRERGRLrlszN6wYcNQXl5usX3VqlUYMWIEAGDlypU4efIk0tPTrTfOYMDU\nqVMRFBQEAPDz80NsbCySk5MB/PcXCG/zNm/zNm/zNm/ztty3zf+/dOkSAGDTpk2qHLOn2KLKX3zx\nBdavX4+MjAx4enpafYxaBzoSERERPUutdUsrJV507969+Oijj7Bjxw6bhR457/u/PMg+zM45zM85\nzM85zM9xzE6fFCn25syZg6qqKgwbNgxxcXGYNWuWEs0gIiIi0j1eG5eIiIioBai1blGkZ4+IiIiI\n5MFiT8c49sJxzM45zM85zM85zM9xzE6fWOwRERER6RjH7BERERG1ALXWLezZIyIiItIxFns6xrEX\njmN2zmF+zmF+zmF+jmN2+sRij4iIiEjHOGaPiIiIqAWotW5hzx4RERGRjrHY0zGOvXAcs3MO83MO\n83MO83Mcs9MnFntEREREOsYxe0REREQtQK11C3v2iIiIiHSMxZ6OceyF45idc5ifc5ifc5if45id\nPrHYIyIiItIxjtkjIiIiagFqrVvYs0dERESkYyz2dIxjLxzH7JzD/JzD/JzD/BzH7PSJxR4RERGR\njnHMHhEREVELUGvdwp49IiIiIh1jsadjHHvhOGbnHObnHObnHObnOGanTyz2iIiIiHSMY/aIiIiI\nWoBa6xb27BERERHpGIs9HePYC8cxO+cwP+cwP+cwP8cxO31isUdERESkYxyzR0RERNQC1Fq3sGeP\niIiISMdY7OkYx144jtk5h/k5h/k5h/k5jtnpE4s9IiIiIh3jmD0iIiKiFqDWuoU9e0REREQ6xmJP\nxzj2wnHMzjnMzznMzznMz3HMTp9Y7BERERHpGMfsEREREbUAtdYt7NkjIiIi0jFFir0lS5bAZDIh\nNjYWQ4YMQUlJiRLN0D2OvXAcs3MO83MO83MO83Mcs9MnRYq9999/H7m5ucjJyUFqaiqWLVumRDN0\nLycnR+kmaBazcw7zcw7zcw7zcxyz0ydFir02bdo0/L+qqgodO3ZUohm6d/fuXaWboFnMzjnMzznM\nzznMz3HMTp9aK/XCixYtQlpaGry9vXH06FGlmkFERESkay7r2Rs2bBiio6Mt/u3cuRMAsHLlSly5\ncgU/+tGPMG/ePFc143/apUuXlG6CZjE75zA/5zA/5zA/xzE7fVJ86ZUrV65g+PDhOHPmjMV9ISEh\nuHjxogKtIiIiIrJPcHAwLly4oHQzLChyGrewsBBGoxEAsGPHDsTFxVl9nBoDIyIiItISRXr2xo4d\ni/z8fLi5uSE4OBjr1q1Dp06d5G4GERERke4pfhqXiIiIiFxH0gSNlStXIioqCiaTCXFxcTh27FiL\nvPiHH34Io9GI8PBw/Otf/2rYvmXLFphMJkRFRWHhwoVW9928eTNMJhNiYmLQv39/nDp1quG+vXv3\nIjw8HEajEb/5zW8atv/9739HZGQk3NzccPLkSYvnvHLlCnx9fbF69eoWeX+Aa7K7ffs2fvjDH6JN\nmzaYM2dOo/uysrIQHR0No9GI9957z+r+WskOkD+/RYsWoXv37o2WB3oW87Oe38OHD/Haa6+hV69e\niIqKwi9+8Qur+zM/28dfSkoKYmNjERkZibfffhu1tbUW+2slP7mzMxs5ciSio6Ot3qeV7AD580tO\nTkZ4eDji4uIQFxeHiooKi/2Zn+38ampq8OMf/xhhYWHo1asXtm7darG/ovmJZnz77bciKSlJ1NTU\nCCGEuHXrlrh69WpzuzXr7NmzwmQyiZqaGlFcXCyCg4PFkydPREVFhejevbuoqKgQQggxdepUkZGR\nYbVdd+/eFUIIsWfPHtGnTx8hhBB1dXUiODhYFBcXi5qaGmEymcS5c+eEEELk5eWJ/Px8kZycLLKy\nsiyec8yYMeKNN94Qv/vd75x+f+Y2uiK76upqcfjwYfHnP/9ZzJ49u9F9CQkJIjMzUwghxKuvvir2\n7NljtV1qz87cTrnzy8zMFNeuXRO+vr5Ntov5Web34MED8Z///EcIIURNTY0YMGAAjz8rmjr+Kisr\nG/4/ZswYkZaWZrVdas9PieyEECI9PV1MnDhRREdH22yX2rMzt1Pu/Gy9v2fbxfys5/fLX/5SLFmy\npOG2uYZ5tl1K5ddsz155eTk6duwId3d3AIC/vz+6dOkC4GkvUnJyMnr37o2UlBSUl5cDePoLYe7c\nuYiLi0N0dDSOHz9u8bw7duzAhAkT4O7ujqCgIISEhCAzMxNFRUUwGo3o0KEDAGDIkCFIT0+32D8p\nKQnt2rUDAPTp0welpaUAgGPHjiEkJARBQUFwd3fHm2++iR07dgAAwsPDERoaavV9bt++HT179kRE\nRERzkUjmquy8vb3Rv39/eHh4NNp+7do1VFZWIjExEQDw1ltvYfv27Rb7ayE7QP78ACAxMRE/+MEP\nmmwX87Oen5eXFwYNGgQAcHd3x4svvoiysjKL/Zmf7ePP19cXAFBbW4uamhqrC85rIT8lsquqqsKa\nNWuwePFimxei10J2gDL5AbCZmxnzs53f559/3uhshrmG+T4l82u22Hv55ZdRUlKCsLAwvPPOOzh4\n8CCAp3+M5syZg/T0dJw4cQLTpk3DokWLAAAGgwEPHz5EdnY21q5di+nTp1s879WrVxEYGNhwOzAw\nEFevXoXRaER+fj4uX76Muro6bN++vdlr527YsAHDhw8HAJSVlaFbt26NntfaF873VVVV4be//S2W\nLl3aXBx2cVV2ZgaDodHtsrKyRpkGBAQ0+97Vmh0gf36OYH7W3b17Fzt37sSQIUOabCPzs/TKK6+g\nc+fO8PLyQkpKSpNtVGt+SmS3ZMkSzJ8/H97e3pLaqNbsAOWOvalTpyIuLg4rVqxoto3M77/MVx1Z\nvHgx4uPj8cYbb+DGjRtNtlHu/JpdesXHxwdZWVk4dOgQDhw4gPHjx+PXv/414uPjcfbsWQwdOhQA\nUF9fj65duzbsN2HCBADAgAEDcP/+fdy/fx9t27ZttkF+fn5Yt24dxo8fj1atWqFfv35NrrV34MAB\nbNy4EUeOHAHg2Bf40qVLMW/ePHh7ezf7y8YecmdnLzVnBzA/ZymVX11dHSZMmID33nsPQUFBNh/H\n/Kzbt28fHj9+jPHjx2PTpk2YOnWq1cepOT+5s8vJyUFRURHWrFkjaVFgNWcHKHPsbd68GV27dkVV\nVRXGjBmDtLQ0TJkyxepjmV9jdXV1KC0tRf/+/bF69WqsWbMG8+fPx1/+8herj1ciP0nr7LVq1QqD\nBg3CoEGDEB0djU2bNiE+Ph6RkZH49ttvJTXs2TcTEBDQqMeutLQUAQEBAIDXX38dr7/+OgDg008/\nRevW1pt56tQpzJw5E3v37kX79u2tPm9JSUmj3i5rjh07hvT0dLz//vu4e/cuWrVqBS8vL8yaNUvS\ne2uKK7KzJSAgoKFbGGic6bO0kB0gb372YH62mQcpv/vuuzYfw/ya5uHhgTFjxiAzM9NqsaeF/OTM\n7ujRozhx4gReeOEF1NXV4caNGxg8eDD2799v8VgtZAfIf+yZix5fX19MnDgRx44ds1rsMT9LHTp0\ngLe3N0aPHg3g6fJyGzZssPpYpfJr9jRuQUEBCgsLG25nZ2cjKCgIYWFhuHnzZsN1bWtra3Hu3LmG\nx23ZsgUAcPjwYfj5+VnMbhw5ciT+9re/oaamBsXFxSgsLGwYa2bu/rxz5w7WrVuHGTNmWLTrypUr\nGD16NL788kuEhIQ0bO/duzcKCwtx6dIl1NTUYMuWLRg5cqTF/t+vhA8ePIji4mIUFxdj7ty5WLRo\nUYsccK7Kztp7AIAuXbqgbdu2yMzMhBACaWlpSE1NtdhPC9kB8ucnFfOzfB9mixcvxv3797FmzRqb\n7WJ+lu8DAKqrq3Ht2jUAT3sKdu3aZXXBeS3kJ3d2P/3pT1FWVobi4mIcPnwYoaGhVgs9LWQHyJ9f\nfX19w+zb2tpa7Ny50+qMZuZn+T6Ap0XhiBEjcODAAQBARkYGIiMjLfZTNL8mp28IIbKyskS/fv1E\nRESEiImJEWPGjBG3bt0SQgiRk5MjBg4cKEwmk4iMjBSfffaZEOLprJ65c+eKuLg4ER0dLY4fP271\nuVeuXCmCg4NFWFiY2Lt3b8P2CRMmiIiICBERESG2bNlidd8ZM2YIf39/ERsbK2JjY0VCQkLDff/8\n5z9FaGioCA4OFqtWrWrYvnXrVhEYGCg8PT1F586dRUpKisXzLl26VKxevbq5WCRxZXY9evQQ/v7+\nwtfXVwQGBoq8vDwhhBAnTpwQUVFRIjg4WMyZM8fqvlrITghl8luwYIEIDAwUbm5uIjAwUCxbtsxi\nX+ZnPb+SkhJhMBhEREREQzYbNmyw2Jf5Wc/v+vXrIiEhQcTExIjo6Ggxf/588eTJE4t9tZCfXNl1\n69at4bNrVlxcbHM2rhayE0L+Y6+6ulrEx8eLmJgYERkZKebOnavZY08IZY6/y5cvi4EDB4qYmBgx\ndOhQUVJSYrGvkvk1W+w5QsoUbrKO2TmH+TmH+TmH+TmO2TmH+TlH7/lJWlSZiIiIiLSJl0sjIiIi\n0jH27BERERHpGIs9IiIiIh1jsUdERESkYyz2iIiIiHSMxR4RERGRjrHYIyIiItKx/wNYYKPhO3xI\nGQAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }