{ "metadata": { "name": "", "signature": "sha256:f7ab8af70315cd2eed5a951ff27cdc4d3439222a3fe3e4a74f18309e36598d32" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "heading", "level": 1, "metadata": {}, "source": [ "obtain PDO data from JISAO PDO page and plot" ] }, { "cell_type": "code", "collapsed": false, "input": [ "#plot since\n", "tsince=\"2000-1\"" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "# python packages\n", "%matplotlib inline\n", "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "#obtain PDO data from\n", "ds=np.DataSource()\n", "f=ds.open(\"http://jisao.washington.edu/pdo/PDO.latest\")\n", "data=np.empty((0,))\n", "for line in f:\n", " if line[0:2]==\"19\" or line[0:2]==\"20\":\n", " indata=np.array((line.split())[1:],dtype=np.float32)\n", " data=np.append(data,indata)\n", "f.close()" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "# make time series using pandas\n", "date_range=pd.date_range(start=\"1900-1-1\",periods=len(data),freq='M')\n", "pdo=pd.Series(data,index=date_range)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "#plot form tsince to the latest\n", "plt.fill_between(pdo[tsince:].index,pdo[tsince:],0.,where=(pdo[tsince:]>0),color=\"red\")\n", "plt.fill_between(pdo[tsince:].index,pdo[tsince:],0.,where=(pdo[tsince:]<0),color=\"blue\")\n", "plt.xlabel(\"Year\")\n", "plt.ylabel(\"PDO\")\n", "plt.savefig(\"PDOsince\"+tsince+\".png\")" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEPCAYAAABBUX+lAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXuYHHWV9799ne6Z6bnkRjK5kAzXQCBBQO4wKkhQQBS5\nrBEVL+vj5RXddb2wYsKu+vrq47q66HpbL+urqz64+AKuCouZBVcUXUIggYRAIoGQGZLMTDKTmenp\nme73j9NlVXdXV1d1/aqrqvv7eZ55pqenpvpMd9X3d37nd875AYQQQgghhBBCCCGEEEIIIYQQQggh\nhBBCCCEk5KQA/B7AYwCeBPC//TWHEEKIl7QXv8cB/A7AhT7aQgghTU3U59efLH5PAogBGPHRFkII\naWr8FvwoJKQzDGAzJLRDCCGkiemGhHQGfLaDEEKalrjfBhQ5DODnAM4CMKg9uXbt2sLWrVv9sokQ\nQsLKVgDryp/0M6SzAEBP8XEawGUAthgP2Lp1KwqFgvKvjRs3enLeRnzRdtpO28Px5aftANaaia6f\nHv4SAN+DDDpRAN8H8ICP9hBCSFPjp+A/AeBlPr4+IYS0FH5n6fjCwMCA3ybUDW33B9ruD7RdLRG/\nDahBoRiPIoQQYpNIJAKY6HtLeviEENKKUPAJIaRFoOATQkiLQMEnhJAWgYJPCCEtAgWfEEJaBAo+\nIYS0CBR8QghpESj4hBDSIlDwCSGkRaDgE0JIi0DBJ4SQFoGCTwghLQIFnxBCWgQKPiGEtAgUfEII\naREo+IQQ0iJQ8AkhpEWg4BNCSItAwSeEkBaBgk8IIS0CBZ8QQloECj4hhLQIFHxCCGkRKPiEENIi\nUPAJIaRF8FPwlwPYDGA7gG0APuCjLeHgzjuBBx8EDh/22xJCSAiJ+Pjai4tfjwHoBPA/AK4B8JTh\nmEKhUPDBtAAyPg709ACZDFAoAAcPAomE31YRQgJIJBIBTPTdTw9/CCL2ADABEfo+/8wJOMPDQDot\n3n0+D/zpT35bRAgJGUGJ4a8EcAaA3/tsR3AZHgbicXkcjQLPPOOvPYSQ0BEEwe8EcCeAWyCePjFj\neFhCOQCQzQK7dvlrDyEkdMR9fv0EgJ8C+L8AfmZ2wKZNm/78eGBgAAMDA42wK3gMDQEzM/I4mwW2\nbfPXHkJIYBgcHMTg4GDN4/xctI0A+B6AQwA+VOUYLtpq3HYb8KlP6T+fdx7w29/6Zw8hJLAEcdH2\nAgBvBvAKAFuKX+t9tCfY7N1b+vPu3f7YQQgJLX6GdH6DYKwhhIPnny/9+eBBIJdjaiYhxDYU3LAw\nNFT6cyoFPPecP7Z4xeHDwH33+W0FIU0LBT8sHDxY+nM83nyZOr/9LfCxj/ltBSFNCwU/LIyNlf48\nPd18ufijo8BLL/ltBSFNCwU/DExM6Dn4Gs2Ymjk2BoyM+G0FIU0LBT8MDA9LzL6cZhP8kRFgakoW\nowkhyqHgh4GhISAWq3z+yScrPf8wc+CAfKeXT4gnUPDDgLGtgpHp6cp0zTCjxe8PHfLXDkKaFAp+\nGBge1tsqGEkkgIcfbrw9XqFlIpVnJBFClEDBDwP794s3X874OGCjf0ZoGBmRTqD08AnxBAp+GLAq\nsNq8uXF2eM3YmKxVUPAJ8QQKfhiwitPv3i2ZLc3A+DgwN0fBJ8QjKPhhoLytgpF0Gvif/2mcLV4y\nMSG7ee3f77clhDQlFPwwYOXxTk83R5vkfF5fp6DgE+IJFPwwUN5WwcjMTHM0HDtyRO/8OTzsry2E\nNCkU/KAzMwPMzlof88QTjbHFS8bGdMFnWiYhnkDBDzoTE7V73mstCcLM6KheTTw66q8thDQpFPyg\nMzEhrZCtSKfD3yp5bAyIFHdkO3LEX1sIaVIo+EFnYkKKkWqxY4f3tnjJ6KjePmJyUhZxCSFKoeAH\nnfHx2oJ/9CiwfXtj7PGK0VF9rSIel92vCCFKoeAHnYmJ2sfk8+HMxX/hBeDGG+Xx2JjeFjmZ5MIt\nIR5AwQ86ZpufmPHkk97bopqnngLuukuvrtUEn+0VCPEECn7Q0apPa/HCC+GLew8NSdrp3r2VufcU\nfEKUQ8EPOhMT4gHXIpEQ4QwTL74o33fs0Dc/ASSWT8EnRDkU/KAzPm5vy794PHyZOloX0KeeKt3l\nKpdjDJ8QD6DgB53x8dqVtoD0oXnqKe/tUcmf/iTfH3usVPBnZvTdrwghyqDgBx27VaczM8CWLd7a\nopp9++T71q2VxVZsoEaIcmqUcBLfsWqcVk7YBF/z4nfvrqw1oOATohwKftBxUoD0zDOSqWOnMjcI\naLOXmZnKdQrN+yeEOMMitBsSZWhhnAh+NKrHxYOOMfsona4cpNgimRDnFArAq19d9dd+C/63AQwD\naIL+vh4xPm7/2Hhc4uFhYP9+IJWSx7OzUl1rZGwsfHUFhPjNo49KTU4V/Bb87wBY77MNwcZOawXj\nsY8+6p0tKtm/X2+HfPRopYefSDA1kxCn/OhHlr/2W/AfAsDm51ZMTto/Np8Pz3aH+/eXevDZbOnv\nk0nncfzf/Q548MHax73qVXoNACHNQqEA/PCHlu3U/RZ8Ugsngg8Ajz/ujR2q2b+/VOTNag2cCv5P\nfwp84xu1j/vjHyn4pPnYvl1CoR0dVQ8JfJbOpk2b/vx4YGAAAwMDvtniC9rG3nY5fFi+urv15woF\nfXORoPD885KdU41czrngDw8DO3daH5PNSs4/w0WkiRgcHMTgJz8pemHRiiVUgt+SlIc6apFOi5d/\n0UXy8113AXfcAdx/f7DSNffssf791JQMCk44cKB2PyEt+4eCT5qIgYEBDOTzEiZNpXB7Fd0IkAKQ\nCmZm7LVGLv8bYwHWL38J/PrXwFe+otY2t9gR82eecXbOkRERcqtmc5rgG5u1EdIM2Ehl9lvw/w3A\nbwGcCOB5ADf7a07AsLOBeTnT07J4qfHww/L9Yx+TitagYCfP3mmcXUvltEhLw9CQfGclL2k2bHSY\n9Vvw/wJAH4A2AMshaZpEw84G5mb893/L93weePppeTw9Ddx2m3ubrrgC+P3v3Z/HTvtjrX2yXY4c\nkQHSamDTBJ+VvKSZKBQq+1GZ4LfgEyvs7GdrxtCQhDaefVYfMMo932xW8t+dUChI2qeZB33jjfa3\nWcxm7S1GO42zT0xIbr/V+oAm+KzkVcf0NLub+s3oqF7XYgEFP8hMTNQn+KmU5KNv2VL698ZGbF/7\nGvChDzk779CQeBHlDd0KBeBnPwM+8AH759GqbK3IZmXx1g6FgqSwTk8Du3ZVP04LE3GDFXX8wz8A\nixcDCxZIq2vSeIaHgba2modR8IOMkypbI0eOAL/6FfCHP5Sewzjle/HF0li/HZ4odsAoF/yRERHc\nxx6zV/j0/PP21iZSKfthHWO17vbt1Y/TsnicdCGth4mJ8PQ1csvzz+shhQce0J/fscM69ZaoY3jY\nlnNIwQ8ydjcwN+O++4Df/Kb07419eYaHJeTj5Pya4Bs3KwEkZp5KiYd9yy21z7l7t71tG+Nx+7H2\nsTG9H49Vdo82gNiId7rizjuBm1skB0ELk+VywLZt+vPXXgv867/6Y1OrMTRkq/cUBT/I2N3A3IwX\nX6xspGas2j1wQH7WblY7PPKIfC+P1+7Zo4v8U0/V9sqffdbe+kE+b1/wR0f19Qqrv9HWBWZmvPU+\nDxyQvkb1DthhwrjWos2uCgUZeL/9bX9sajWGh21dzxT8IDM+bs8TNiOVqhSb6Wl9ANFi2E8+af+c\nWny2fDH12Wf1WHs8XnuQ2rbNnhBOTzvz8LVq4snJ6gOKFsppa/M2jq+td7RC+qfxfdQWzPftk8/j\n0Ue5XtII9u2zVaRJwQ8yExP2NjA3Y3y8ctU+Htdj+tqqvnEKbkU+r8eky0M6Tz5pb99dDatFVSO5\nXO2KXI3RUX0QSafN4+dHj+oDqNfdOLVZTth2IasH43rI6KgIz86d4nTE48C//7t/trUKNmtWKPhB\nxu4G5mYUCpVebiKhb6hy5IiIn91Uyj179IXW8n12d+xwZptVYVQ5zz5r77jRUV3MIxHzgWJ4WM8O\nikS8FfyhIXkNu+9vmDGuDWmD7c6deurvv/yLb6a1DDZnwhT8IGN3A3O7xGK64Guevl0PdNs2fcZQ\nvgtXrf41RqannS2Y2vXwx8b02dDsrPnaxNCQ/j/k894K/oEDMug+9JB3rxEE8vnStaFoVGL3jz+u\n11ps2cJWFl5js66kNQV/Zgb4/OeDt6PSHXeUiqFqwY9GRay1nHWgdqbO1VcDGzZIda02YzCmes7N\nOYvRPveceIF2sRvDHxnRF62mp6sLvva/5nLeCr722YVlB7J6OXKktBo8mxXBN/7fqZR1qixxj817\nMPDdMj1h40bgs58FbrgBWLHCb2t0Pvc5oLdXBBZwtp+tXcbG9Jz1fF4EcP9+oK+v8tg9e6TLJlBa\nGWsMFb3wgiyA2u3bv2ePrYrAP6O1M+7qsj7OmDmUz5s3Zxsa0mcB2ay3Xqf22Y2NyVdPj3ev5Scj\nIxLq0wbbbFZmg8bU2KC15m42bLZVAFrRw3/0UeBLXwLa253Hnr1mZkZ2rNFQnSs+NydCNDqq56wn\nk9W9rzvuEPEsb4MwM6PHy3fvdtbgbfduZ+mQ6bS9sE55qqjZOsH+/Xo2UaHgXT+duTl9AEynm7v6\n9NChyn5Pjz7qfWEb0bHZVgFoRcG/4Qa56efmJGc8aDzwgC6ITjYwt8PsrAj+2Jh+k2azekGVkWwW\n+OY3zcU5mdQHoz17nC0sP/20/XYJQPUF2HLKp7Rm6ZDlmTtepUwaB9Tp6fAv3D7+ePVMm/KMLUDC\nOU7CdsQdNtsqAK0o+Jp4ZLPB9LzyeeC//kse19taoRq5nO7ha9PsbBb4+c8rj73zzuqx/URC9+B2\n7XLWhM1pLHd62l5b53LhMQvXlId5vGr4dfCgLvgzM3rBWli5917g7/7O/HcjI5VrYUFbG2t2bLZV\nAFpR8I0YF5aCUhE5Owv85Cfy2Gk3y1rMzYkYGXPWAemAafTkn3lGWiRUG3CiUV3wnQq40578MzP2\nZmLl6x1m/X7KM5K8WrQ9eLD0BnQyowkiO3dKrYVZTcjISOXz7e3O92Im9WOzrQLQ6oKvic8jjwBB\n2Su3UJBtCXfv9uamOXCgNGcdEG9U66H/7LPAueeaT9U1IhFdUJ3uSlVPGMXOWkt5+GtysjTU9M53\nVq5FeBVnPngwOA6ECnbtquyTo3HwYOX7msvVXyFOnGOzrQLQ6oI/OSme4be/7Xz/VC+ZmgJOO80b\nQRoZKc1ZB8STv/tuEanrr9e7X1ajUNBtc1JEBdTXv8ZODL98cGxr0z34u+6S7qHlr+1VA7WDB+sv\nmAsie/dKGM9s4xuz9Fd2yGwsNtsqAK0u+Om0hAt+/ONgeWSTk/LlhU0jI7LAabwp83npZ/+LX8j0\nvdbrzs2J4E9MNCZc8dJL1jblcpUik0joxSi33mo+W8rnvbH/wAHnm88HlUJB3v9cTvZGLsdJ8z3i\nDQ62Ag2f4D/9tDrvKZ8X797rVrlBYmzMvCpv/37gXe+yt26Qy8l59uxpTDZGPG4tLMbWyBqRiD5Q\nVOtL71UDtaGh5glpjIzo6xHa/shGWEHrLy+9ZJ50UYVwCf7vfiehjt/8Rs35JiaA7343WN6914yP\nm9+k8bj9Qq+ZGRHZ3bsbU1STTFqHdcbGKmsB5uZkYLPKxInFvBnsm2m/3Oee01P+hocrrxF2wvSX\nD3zAUYPF8Aj+7t3A+vXilatM+8rlWkvwJybMb9KjR51lBb30kizY2tmb1i1zc9bZPaOjlWlp2azY\n+PTTtnOUldFM+7sawwXptOyiZsSLanBij4cfBu65x9GaSXgE/6//Wryxjg61561nz9gwMzlpnYFj\nlwMHZP2jEQt0k5PWXTONvfA1cjnxtFWGAO3STGGO557TB/XJycqwTiuFQ4PGQw85vv/Co3ZHj3rj\niXtVJLJpkzfndcvMjJqmbIcONa5SOZ+37ts/Omr+Oe7dK/njqusZaqG66Z2f7Nqli8rsbOln7tWi\nN/GM8Ah+2Lj9dr8tMCeRUOOBjo42dpNuq9DB2Ji5F79vnz/V1M3k9e7cWf135Z0ySeCpJfivAfAg\ngEPFr/8C8FqvjSIekkjUv4uWkbGx4MSqx8bMp7YHDkhIp5Hkco1Z12gUVoP6yEhldhQJNFbD87sA\nvBvARwBo3Z/OBPBZAMsAfN1b04gnxGJ6a2Q37NsnN3sQCoxmZ83/nwMHGh/OGRmRReJmEX2rdNhD\nh5y1uia+Y+Xh/xWAywH8GsDh4tevAVwB4EPem0Y8IRJRk7UyO+usLbIfHD6s5n/9wQ+s14+M6zUH\nDwb/fbHL0aPWBWQqFv9JQ6kV0jFLsj0EoIXyGJuMQkFd3DUI3r0VyaSa//Wmm4DNm6v/3rheU944\nLcw895y+B7AZZp0ySaCxujKPAFhn8vxaAKoata8HsAPALgAfVXROYkU+r0aQEonGh0uckkyqa3Hw\n+c/bO+7gweYRwb17rUM2hw6pWQ8iDcPqzv9rAP8PwCYAVwG4GsDtAO4u/s4tMQB3QET/FAB/AWC1\ngvMSK1TdoGHwYgsFdbH0zZuBF1+sfVwzNQ6bmrIOZU1MUPDd0uD0bau79jcAzoEI89sAvLV4/DkA\nHlLw2i8H8AyAPwHIAfgRgNcpOC+xwrg9oRvCcKOrbi/9la8A3/oW8CEuYRFFNDh9u1aAcwjAbQAW\nFX9WmYe3FICxJ/ELkMGkhHvvBV7xCuDP9bVFoXnqKWet2Lu6gEsusT4mm5UdBt3o4ctfDhxzjDze\ntcteK/dkErjssgYVRRQKakIx1cIWxbj+H/9Y2vrezvtfi8OHgQcfLH3u7LOBxdX+oFahXtHWLVvM\nuzx3dwMXX1z8IZsFPvMZedzfD3zxizXPa2R0VF0LqEZw8snACUDtdRo7xZDFe3brVokSEWDVKmDN\nGnk8NSWNSJ1EAs88E+gDzP/IwhmzEvwIgI0A3g/x8gFgDsA/Afg7uF+4tfX311yzCW94A3BKPI6B\n007DwKJFwHHH4S1vFNG3mxU2Pi6z++R115k3t7rgAjz8MPCGN9Sf2JHNAh/9aHHQ3rgR73mPVKLX\nWjecmpKNo0648UbZANpLDh8WQxctqn1sPXR2AvPnY/0ZMpnQOh6Mj4vDnbr+euc99AHg6qvxk59I\nrygt9TubBT74QeCzl54DXHih83O2twNLluCadaVNITXGx+Xeib3xjZKeuH27KGFx5Nq0qTgj37hR\n//naNcA555R6DVdeiR/8APibvwlH2vrsLHDBBcB9X14tHozxf1m/HkDxf7345XJgLZJJoL8fb3iF\nlG6EIRroJXNz4jM8/jiAjRuxeTNw7bX2dSebBd73PuALV+nv/+DYGAa14sRotK6iyL8CcD+AVYbn\n+gHcV/ydW84F8EvDzx9H5cJtIZMpFH7yk0IFa9cWCuJe2PuKxQqF6enK8xh54IFCobvb2XnLvz7x\nCf18F19s728ymUJhxw5r28JGT0/p/xiPFwqTk+7O+fWvFwrpdOl5P/xh97b29Zl/LpFIoTA7W/3v\nAOufjXz5y4VCW5u7a6uRXwMD1u+Z1f9ajRUr/P+/gvJ18sn6+3LPPc5155Zban0+5g611Vj7FgBv\nArDH8NxuABuKv3PLHyGzxpUAkgBugCwI26K7W4EFhJjQ6h4oaV6sLu04ALOmKwdQO/Zvh1lIuOhX\nAJ4E8GMAtrtxnXKKAgsIMSGfb57aKUKMWAm3VRqGqhSNXxS/HLN6tcS8mmUnuWYj7D212tqc15UV\nQ/mEeEYy6S7z1+q2PB3VC6wasK+dNf39FPwgs3Chvod42Jg/X7z8gmkUtDpB7YhNmoNYTPagcSP4\nViGdGIBMlS/f/bf+fuc3JGkcS5eW/qyqwLcRZDLANdf4bQVpZtrbnf9NMum+tMTqFkxDmqTdAema\n6bvIG1m1insvBJljjy39OZ1u/E6D9dLVBVx/vd9WkGam/P6wg4rO5laC/z1IO+RtkL74X3D3UmpJ\npyXlmwST/v7SGgmtGC0M9PRIsd+llzZmj/ZWodU6KUciEto08+ZPOKG+86VdBtOtBH81gDcD+BqA\nawFcbHGsLyxb5rcFpBpLlpQ2Wlyxwj9bnDJ/vsxG7r8/PGGoMJDJ+G1BY+noAK6+ujLjK50GVq50\nfr5CwV6dmxVWl/NslceB4cQT/baAVKOvrzRT5/jj/bPFKfPn+21BczJvnt8W2EfFQB+LyVpQeeg5\nkajPWZ2bk5mnG6z+LS1LR/s6zfA4EJt2rlmj1gOLx6WcPharb1GlXmZmwp/GWM6SJfqieiIRrsHZ\nq64TYcDLENaCBd6dWzUq6jCmpoDzzzdvp1GP4Ody7vtROcnSiRsed7l7WTUcf7xaYT7/fOlvcfiw\nNxe+WQwznQY+8QmJeTcTfX16ymwqFZ6QTiIB9Pb6bYV/eFlwFqaBVIUjmUzKrOa440qfz2aB5cud\nn292VpoFulkLCXWEsnxh0C3xOHDqqRJ7K/aIUkYsJoOIcSBJpYB77hHBb7bFwfnz9cKlaDRcgu9V\n244wfMZetvNftCgc7wGgZg+bVcUuZOvKtpGam6svbNjeLoNIPYOFRugF36sL9MYb1S4yJZPSLbGj\nQ38uGgVe+Up1rxEkIhHJdgHkMwqL4Mfjut21cFJZ294e/CyVZNJbQe7tVd8tNJXyZjE4m3Xv5Wvt\nj1/2stIEhnnz6nuftazE8gHECaEW/MWLvdtW9fLL1VXxRiLSvfe660ovopNPDo/HUw/aFH5mRj6r\nMBCJ2PfwnVTWZjLBz/hpa5Pr1M41WU8bid5etSGjdBr46lfV73Oj4cbWREKEHtDbwGj09dV3Tu26\nvOqq+vL4gZALfiTinZBkMrLJgKpzffzjIvDGi/Pcc9WcP6hoF3Zvb/C9WyN2PXwnhEHwo1EpOLNT\nwV5PG4meHrXXwZVXAjff7F1WlZtEinQaOOkkebx2bWmmTj0pmYC+tvT2t9e/6VrAL8HaXHSRd+d+\n05vcFzoAsvHK+efLuTSvt7NT9pZoZrQLe8kSNeeLxbwXzXzemxh+GIoE83lgYMA7W1UuhqdS+r3v\nVQaYm8Epn9cFf9GiUq++3hTlhQvrt0cj9IJ/+eXeXaDHHacm5njGGfqU7tRT5XskIiN/M6NlHq1a\nZX2cXbq7vRf8XM47Dz/ovZ+00Nt553lz/p4ede9BMin3FQCcdZaacxqJRt2FW1Op0sy7q66Sc6ZS\nzjx8Y+xfRZZT6AX/kkvsxfHrudBUeCSRiHzYGueeKx/81FTz9/TXPPuTT1Zzvp6exgi+Vx6+iswP\nL8nlZEHxyiu96XvU26tuzW1qSneY1q4tTYZQgRtvuqNDQl7GkNBrXiPXgJMsm0ikdB1BRXuS0Jf7\nHHusvMHT07WPdeqt9/S4v0kLBenJorF2rUwVly4Nx/6mbtCmsao8fJUeohVGr0oVmYx3CQaq6OqS\nAfUtb3Ffwm9GT4+6rLoFC/TsnNWr1a8RrVgB7NhR398mk8A73lH63MUXi0Y5qbJNp/XYfzKppnAt\n9B4+YK/6LJ12PkXr6XHfna6trXTxd80aOae2gt/MaB6+qpRMFQNwLbyqsA6D4GutD3p61CUsGGlv\nVzdga+EcQARfdabOqlX1fV4dHcBtt1U6De3tcu8fPWpf8ONx3QZVBYFNIfhXXFH7Rq0nzq/CIzn7\n7NKp3XHHyc9eeFBBQ/PwVQq+16KpOjSgkU5LwU05qZQ3M4p68LqjqYpuj4CInzFZo6tLbS5+NCr3\naT3OXi4HvPvd5r97/evl3rcbLopE9ChAPK6mF1FTCL4dD7+ehbhUyl3MOJkEXve60uficeD002Ug\naHbmz5eLVJXgd3d7WwkKiHh4QSRiHhdPJoOzT0D5pjVeoEKY29srZ8gqm/O1tYkXbjZA12LduurO\n52tfKwu5diMNuVxpkRU9/CLHH1+7SKLeN8vNFP9jHwNuuqny+QcflAKXZicSAfbuVSeiiYQ3TeaM\nayleCT5Q3ZMPymJuvfnhTlDx/k5NlYZ0ALUhUs2bdjoQmzl4Rs44A3j0Ufvnm56W7pjRqITC6OEX\niURqX6z1Fme48Uhe9SrzaXJHR3NX2BpRHSLxIsaupcoC3rbwNbM9l/OuUtQJyWT9FaBOUOGldnRU\nhkXOOEPdtRGNip1Oz9fWBlx2mfUxTu6HZFLi/u3tcp3QwzdQK9Wp3jQrrxppkfrwom/KeefpM0Qv\nBd/sZp+ZMY9rNzqDq62tMbuSua2KjUYlNFLO4sVqZ3+9vc6dldlZtTMNrdNmNCrXCQXfQC0Pv95q\nTy+KcEj9eBFyWbtWF10ve7abCUh3d2XWRirV+IyeaLQx/Y7cvr/pNPC+91U+39mpbtY8Nyfi6jTR\n49xz1aaHLlwogj85qW7Bu6kEv1ocP5Go37Pg7kfBwote9aefLgIbiXjbs91ssDrmmNJ+6fG4LOz9\n5V96Z4cZ+XxjPHy3rzFvHnDOOZXPq6y218InTpyLVEp2t1JJX5+Ifiym7v9rGsHv66u+KJZM1i8U\n5R5JUFLoWhUvQi7Ll8uimJe98AFzAVm+XPLINdrbgW9+UwoKG9lwbmamcSGdekMv6TTw/vebe/Id\nHeoWv2dmZGbv5FpIJOrfjaqapqxYIf9rX5+6UGZTCX61FMpYrH7BN94ExrxY4g9ehFwiEeCEE2Qq\n72UIz+zc/f3y2toC4dSUxIF7ehp7rWltFbzGTk/8aqGZfB5461vNf6eydUUyKYOSE82Yna1/drhs\nWaXtsZjeAvn449Vdl4EX/PJ+EtXo66ueNxuJ1C/48+frr1+teIY0Dq+80LPOks/WSw+//BpMJCSc\ns2qVfo3194vH19PT2H2OtbYKXmOnRbJZNW48Lo0Sq33+nZ3q7k0tfOJkAJybqz9LaPnyyo3OUyl9\nTeW009SFlkMh+HZWy5csqd5Px00Oa0+PfjPG4/Z69hDvmDfPGyHUWlV76eGXC34qJcVOK1fqYjUw\noB+rQoDtLmQ2wrsH6v+/4nHgC1+o/vvOTndtUDIZfeahhU+Ms8la7+PsbP2C39FROeuJx3XBf9nL\n9FbLbvFuoZqiAAAUHUlEQVRL8K8DsB3AHICaiUx2BL+zs7oQzM7W7+GXe1r08P3Fq1DH6afLd69j\n+MbZajQq0/kVK8TD6+zU48CqBh67vWtU9Fq3Q60GeKlUZcFTMgls2GBdTdvW5q5Pj7ESWnvvu7v1\nz8vOeoqbHbLKw0GFgj6bedObgH/+5/rPbcQvwX8CwOsBPFjrwELB/shZzUtxk8Pa26uP7rkcF239\nxqtQh7b/qJeCX+6UzM6Kh9/WJoPB1JSegdLb21jnolHZaL291rH2ZFLy7I3iWSgAn/609Xndrq/N\nzuoFcNp7YRyga6XJut26sbyOKJfTPfxIRF3KqV+CvwPA03YOzOftF0BUyyN2co5yjB5JNmud79/e\nbr8THqmPnh5vqpQzGeD2272tNs1kSgV/elp/vWXL5PrRWkk3olGcES/TUY3U6kA7NQVce21pzvmS\nJfbWbtw4Y9PTejhNC+V0ddnPlHLbD6l89pLNejPrCnwMf27OvlhXq7Z108rA6GllMtYFXIlE829q\n4jde9sT/5Ce9bWSWyZReh+m0LlLHHy+xWu33Kjq1Otm1SdU2lLWo9X8dd1zlPVTeN6cabgQ/nQau\nu04ea4NL+T7EVjNLt2HGE04oPX97u9oN3zW8zAO4H4CZz30rgHvsniSb3YQvf1k+zIGBAQxow7AJ\nxgIWI25yWI2e1oIF1h0FW2EXK79pRE98ryivBjV6cJdfXiouyaR4l27+144OuSe2bbOeLajaXMMO\nsZgMqtWSHy65RK8u1Wyz20rcTS+defP0PjjabMdYN9HVJY5ftffRbah3xQo5x8SEbo8TBgcHMTg4\nWPM4LwW/RhshexQKm/CpT9mL2y5fLhdIuQfhJi5rbMm7ZIl1q9+2Nlbmek13d/A3EqlGueNhdB7M\neqh3dABjY/W/XiQifeOfeUYXEjMSicZl6QDVd6jLZGRnqExGPO7xcflubBFc67z1smiRpMQuWaK/\nF11dMptsa5O1lUceqf73btseLF9eGj5yGmIrd4Zvv/120+OCENKxnHRGo/YX6fr6zKfkbsrxNY8E\nkEKIJUuqT/urzTCIOrq7Jb4ZRsqLg2pt/ei2unJ6WgS0Vhw6Fmus4FeLx+fzenqsFp417l1bCzft\nB7TB99Zb9U3RNa9eE3wr3Hr4y5frjkw0Wvv16sUvwX89gOcBnAvg5wB+Ue1AJzHVatW2bqer2hrA\nqlUy8laL19n1REj9aFWQYSST0deDYjGJ21rhtlFcPC7ZR7XCQpFIYwX/ssvM79OpKb269MQT5Xs6\nbd/bdTNAaoPv+9+vb07U1SUiXCiIAFvNLN22Zl66VC++6uwE1q93d75q+CX4dwFYDiANifNfUe1A\np4Jv9qG4zUDQsiuWLZNzmV2sqVRr7FMbBLzahtBrjPvatrfXFny3IrxsmdwTtYoF8/nGCv6ll5p7\n48aFytNOk+9O1sTqHSATCfOEj0xGwrmTkzLzsJpZur0mk0n9PZmaKt3CUSVBCOlY4iQ2Vq3a1m3b\n1+5u8cyWLBHBN/OY2tq4YNsoVHZGbCQdHXpKYiRSe1s+t+tBJ55YmRlkhpvCxHq44ILKVgJA6Uxc\nGwzPP9/+eestVmtrM9cIbYe1zs7aKZIqnBAtU6q/37t6kKYS/HTavErPrfeiFYv09Yngm6WVZbOl\nXQ+Jd3i5DaGXxGJ6OGp6uvaaj5s87EhEYt92wjXZbGMFv7fXvF7FKLrae6PF0+1Qr+DHYtWdwnRa\nvP9IxDpkpMIJ0RJCrrzS/bmqEXjBdxob6+8vXaRSIfiap6Wt4JsVjkQijctlbnUaKU6q0Rb37HRj\ndNMorqND779S6zyxWOMryF/96sqZhzGsos1+7C7YAuIV11Nvk89bC74227D6vFS0L9Zep9Y2iW4I\nvOA7HTnvvVe8cM2T0vandIN2wxxzjJzPzKZVq1pnn1q/aWS8WTXajLXWlpyAvVbC1YjFdNGs9Vp+\nhMguu6xSJI1ZSwsXAh/8YO11DiOdnfUVK83MVHfWMhl9PcFqLVDFrFPbzvDCC92fqxpNJ/jHHgv8\n4Q8yYmuVhm4Ff+FC+UC1QcTsfNpFQbynUUVCXqDFerUsFCt6e+uvtjSGjGpt/+nHvs0XXVS63tbW\nVlqXEIkAX/yis4wsqwaKVlitYSxcqA+cVmuBKjz8/n5ZHPYyKSHwgl/PG7l0KfDQQ3IhT0y4F/ze\n3tJpcbngdHYCV13l7jWIfRrV98ULtJvZTqjCTu/4aszN6WsAVtt/aq/TaBYtKl2UbmtzHxLt6Kjv\n/bIKBX3kI5JVBFSvso/H3adlAhK7//nP3Z/HisALfr1TpZUrgV//Wqb/brN0entLvaTyBlvpNHD9\n9e5eg9jHzTZ5fqNtNGLHw3fTKM7Y5dVq+0/Av+pwY1abijWwejcyt5oxXnONbtfSpeYDSjyuZoPx\neNz7cGXgbxs3081164DhYfex9de+Vm+fC5TGRDs6gNtu86bRETFHEzCrdgFBJZORRcJaKZmAiHa9\njeKMsyCr7T+BxvXCL+eMM8QpKxQkrOKX4Nt93YULZSai9fnRiMXUePiNoKkFH1CzkDp/fqkXtGyZ\n3ED5vHzY73iH+9cg9rn0Une7G/mJFj6xI/g9PaU98Ts6JH/dTkM1Y9qj1fafgH/ZZaeeKkJ59KjE\n81UIvhVtbebFU3YW0AERfLOZZTQaHsEPfEgniEU2l14qaWWXXw584xvh+bCbhWOO0Uvww0Zvr8xO\n7Cw89/bqA1t7O/C1rwHvfKe91zFmvFht/xmL+bcmctJJeogkkXB/H3V0VB8M02nzjKdIpPaitsbC\nheYOZCSiJqTTCALv4QexjP7ss4FfVO3+QxrBG98IfOYzflvhnN5eiQXbmXkaWzEsWCDb/E1MiDCW\nhxWMJBKlgt/ZKffR4cOVx7a1+ZfmetJJ+kCkwoby5nRGqqW3trfbb3q4cKH5TCkSCY/TF3gPP4iC\nT/zn6qv9tqA+uruBk0+2d2w0KrOBZFLaJ0cistBZa72ora0yseB1rzMfZBrdKdPIvHm6ELspMtPo\n7Kze4KzaQBCP1+5aqrFwYfXNWyj4iqDgEzPOOiuc18all0pHRrt0dIiIbdggP69eXbsZWixWKfgb\nNpinOKsoTHSDFk6x2ljILsZeReVMTZkvgOdy9gVf66lVTqHAkI4ywnhTE++JxaRaM2zVzWee6ez4\nTEbEUFuzWLhQPHyrzo1a3ycjr3iFufdbKPhbubxmjezIZVd0rbCqSp6dNffyZ2bsDzZaP53yTWkK\nBXr4yqDgk2p8+tPAm9/stxXesnQp8N73lj5XK+Zs3BxdI5GQ9OJy5ub8FXytpbjVTnJ2iUSqt1Pv\n7TUf8ObPd1asZTYbyucp+Mqg4JNqnHIKcPrpflvhLT/6EXDzzaXP1arSrdbZ8aabKu+nmRl/BV9r\n8KYqNbRagdlJJ5mHe+xm6GiYpXDOzoYnpEPBJyTA9PVV5n6feab1xkDz5pmHuq64Anj960u90bk5\nf1OfVQt+NeE99VTzGL72+nYxa4E+N0cPXxkUfEJKOeWUUk9WK/yplfESjwPf/77UjmjCqG3f6Rf9\n/WKXKsE3E95EQoS93PuPx2UgcMLq1ZWD7exs49tL1wsFn5CQsXp1aXpgRwewaZMuOrUqRzdsAL7z\nHRFCFV0e3ZBIAF/6kv1c+FqY6UUqJYve5d5/e7vzxeJVqyoFPx63bl0RJAJvJgWfkFL6+krDE/E4\ncMst+iDQ31/7HDfcIIOEk37zXvHe99bf978cswEsEpFFYTMtsfNelR9fHhpSZXsjCLzgh+nNJKQR\nRCKSz9/RIZ7lhg1yn5x1lrNWAbfeKs3Lmgmz9YhcTjz88t9NTzv38FeurNyP12o9JWgEXvDDlmdN\nSCO4+27gW98Sj/Nd75Lnrr5avM/ylEwrwhKKsItZs8VcTtY1ylut11N01tlZGa+n4BNCPCUSAW68\nEdi1S09NfeUr5bsTwW82zAR/3jwR9/LfLVlSn0NZvgF7WBZsAQo+IU3DunVSSKSiiCmsaO2ntayl\naFSvpC3f2ctp/F6jfIE5LDn4AAWfkKYhFhOPP6yto1XQ3a2HqbZuFW9c212sPHxj3HHLCWvWlM4M\nwiT4ge+lQwixj5+N0IKA1iL57W+XrqSPP673HTLuQZBM1u/hH3+8iLzWojpMmYQUfEJI06Btc/jh\nD8vPxrh9V5eksM7Outs0vb+/tEV1mASfIR1CSNNwwgnAe95jHtbq6tLTvKPR+he3V60qbcQWJsH3\ny8P/PIArAcwAeBbAzQBM9uMhhBD7XHSRfJmRyeidMd1smr58eemeBH5XKzvBLw//PgCnAlgL4GkA\nH/fJDkJIi9DVpS/outk0PR4v7TBKwa/N/QC07Qh+D2CZxbGEEOIaozAnEu5CMcaWFOUFXUEmCDH8\ntwP4D7+NIIQ0N5mM3gfH7R4A558vs4VolDF8jfsBLDZ5/lYA9xQf/y0kjv9DD+0ghBB0denbHLrd\nNP3MM0Xos9nw9MIHvBX8y2r8/m0AXgPgVVYHbdq06c+PBwYGMDAw4NIsQkgrksno2TXl7RGcsm6d\nzBbi8WAUXg0ODmJwcLDmcX61JlsP4AsALgFw0OK4QsFsmxpCCHHI1JR4+bOz0k76H/+x/nPNzYnQ\nJxLA178evL2VI1IKXKHvfsXw/wlAJyTsswXAV32ygxDSIqRSEtJJJNz3G4rFpACLIR17BGDbBUJI\nKxGJSIVtJKKmo+jZZwM7dwYjpGOXIGTpEEJIQ2hvl3CMij10zztPP2dYoOATQlqG9nYJw6jw8Neu\n1c8ZFij4hJCWQSu+UuHhaxvPMKRDCCEBJJOROL7Z3rf1nOvKK4HFZtVGAYXtkQkhLUN3t/sqWyP3\n3FP7mCBBD58Q0jL09ACLFvlthX9Q8AkhLUNvr7Q3blUY0iGEtAwLFuiboLQifrVWsAtbKxBClLFv\nn1TbNruXX621AgWfEEKajKD10iGEENJgKPiEENIiUPAJIaRFoOATQkiLQMEnhJAWgYJPCCEtAgWf\nEEJaBAo+IYS0CBR8QghpESj4hBDSIlDwCSGkRaDgE0JIi0DBJ4SQFoGCTwghLQIFnxBCWgQKPiGE\ntAgUfEIIaREo+IQQ0iL4Jfh/D2ArgMcAPACgyXeYJIQQ//FL8D8HYC2AdQB+BmBjI198cHCwkS+n\nFNruD7TdH2i7WvwS/HHD404ABxv54kH8IOxC2/2BtvsDbVdL3MfX/jSAmwBMAjjXRzsIIaQl8NLD\nvx/AEyZfVxV//7cAVgD4LoAvemgHIYQQABG/DYCI/n8AWGPyu8cgsX5CCCH22QpZIw0EJxge/y8A\n3/fLEEIIId5yJyS88xiAnwJY5K85hBBCCCGEkECxHMBmANsBbAPwgeLz8yCLx08DuA9Aj+FvPg5g\nF4AdAF5teP7TAPaiNHXUS1Ta/kvIrGk7gH8BkPDScKi1fbD43Jbi1wIP7QbU2Z4x2LwFwAF4n4Sg\n8n2/ARLv3Qbgs55aLTi1fV7x+HEA/1R2rqDfq1a2N/pebSoWQ1+g6ASwE8BqSIHXR4rPfxT6BX0K\n5M1OAFgJ4BnoC9gvL56vUReRSts7Dee9E8CbvTK6iErbNwN4mcf2GlFhu1mW2x8BXOiJxTqq3vf5\nAJ4rfgckY+6Vnlru3PZ2ABcAeDcqRTPo96qV7Y2+V5uanwG4FOLNHFN8bnHxZ0C8nY8ajv8lKmsB\nGnURlaPC9gSAuwGs985MU9zYvhnAmQ2wsRoq3vcTIR5no6nX9rMB/Kfh+ZsAfMVTSyupZbvG21Ap\nmhpBvVc13obqtjf0Xm3G5mkrAZwB4PeQD2G4+Pww9A+lD8ALhr95AcDSBtlnxUq4t/1XxeOnIDd2\no1iJ+mzvM/z8PUhY5BNeGmrCSqi5Zm4E8CPPrDRnJep/33cBOAnAsZAizGvQ2L5WK1Hbdo1C48yy\nxUq4t73h92qzCX4nJOvnFlSO+gVYXzR+X1CqbL8cwBIAbQDeqtJAC9zYrrEBUotxUfHrJpUGWqDy\nmrkBwL+pM60mbt/3MQDvAfBjAA8C2ANgTrGN1VBxzfiFKtsbfq82k+AnIB/C9yFTLUBGz8XFx0sA\nvFR8vA+lnsyy4nN+odr2bPF8Z3thbBmqbH+x+H0CwA8h8VmvUfm+r4V4yVu8MrYMVbbfCwnvnA9Z\ndNzpncl/xontQUO17Y28V5tG8COQle4nAfyj4fm7oY+cb4X+Ad0NmX4nAayCFII90hBLK1Flewfk\nYgNEeK6E9+KjyvYY9KycBKT9xhNeGg7118xfQAaqRqDSdq0Gphfi7X/LM6sFp7Yb/85vVNnux73a\nVFwIIA/JRNDS49ZD0qL+E+ZpardCshV2QKZWGp8D8DyA2eL3T4bE9kWQm3grgMcBfB7e3ySqbO+A\nZLdo6YFfDJHtGs9CFm0bgUrbfwhJDdwO4HqvDUd9tv8JwCFI+OR5ACcXnw/DvWpmux/3KiGEEEII\nIYQQQgghhBBCCCGEEEIIIYQQQgghhHhFBMBDKG1mdR2AX/hjDiGEEC85FVJN2Qbpm/I0pEK1HuKq\njCLELazuIsSc/wNgElIFPAHpKLkG0vphE6ScfiWAfy0eAwDvB/AwgAEAfw9gBFJZeVLDrCaEEOKY\ndkgbgscBfAbSzROQsvmdxd+nIbMAQPrT/KH4eAD6IEFIYOB0kxBzJiFtgycgPWauAvDh4u/aIN0n\nhwDcAemUOQcRfY1HILtJERIYKPiEVCdf/IoAeANkwxAjmwDsh/TujwGYNvzuaAPsI8QRzdIemRAv\n+RX0DasB2ekIALogXj4AvAUi+oQEFgo+IdYUIAuwCUg8fxuA24u/+yqk//ljkIXZibK/I4QQQggh\nhBBCCCGEEEIIIYQQQgghhBBCCCGEEEIIIYQQQkiQ+f+CgJfllCnJrAAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }