{ "metadata": { "name": "poking_eyeball_rps" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "# import some modules using standard naming conventions\n", "from pylab import *\n", "import pandas as pd\n", "\n", "# find the .csv files and read them into pandas data sets\n", "import glob\n", "flist = glob.glob('/usgs/data1/csherwood/poking_eyeball/*.csv')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "def Phi(mm):\n", " return -1.44296504*log(mm)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "def my_read(csv):\n", " phi = float(csv[-6:-4])/10. # extract phi size from CSV file name\n", " df = pd.read_csv(csv) # create Pandas DataFrame from each CSV file\n", " mg = df.geom_mean.mean() # mean of geometric means\n", " m = df.arith_mean \n", " ma = m.mean() # mean of arithmetic means\n", " sa = m.std() # std of arithmetic means\n", " return phi,ma,sa,mg" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "# list of stats for each file\n", "stats = [my_read(f) for f in flist]\n", "# convert to array\n", "stats = array(stats)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "# create a data frame from the stats array\n", "df_stats = pd.DataFrame(stats,columns=['Phi','Ma','Sa','Mg'])\n", "df_stats" ], "language": "python", "metadata": {}, "outputs": [ { "html": [ "
\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
PhiMaSaMg
0 0.0 0.557322 0.067476 0.468722
1 0.5 0.543682 0.134539 0.430470
2 1.0 0.487743 0.082410 0.381277
3 1.5 0.373238 0.048243 0.267245
4 2.0 0.326144 0.031638 0.209683
5 2.5 0.301667 0.034709 0.199168
6 3.0 0.235943 0.032090 0.134345
7 3.5 0.238555 0.015709 0.135005
8 4.0 0.204032 0.015977 0.106308
\n", "
" ], "output_type": "pyout", "prompt_number": 6, "text": [ " Phi Ma Sa Mg\n", "0 0.0 0.557322 0.067476 0.468722\n", "1 0.5 0.543682 0.134539 0.430470\n", "2 1.0 0.487743 0.082410 0.381277\n", "3 1.5 0.373238 0.048243 0.267245\n", "4 2.0 0.326144 0.031638 0.209683\n", "5 2.5 0.301667 0.034709 0.199168\n", "6 3.0 0.235943 0.032090 0.134345\n", "7 3.5 0.238555 0.015709 0.135005\n", "8 4.0 0.204032 0.015977 0.106308" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "# list of full DataFrames\n", "df_list = [pd.read_csv(f) for f in flist]" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "# 1st column of stats array is the phi size\n", "phi = stats[:,0]" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "# make a dictionary of phi_size: full data_frame\n", "d = dict(zip(phi, df_list))" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "plot(phi,phi,'-k')\n", "plot(phi-0.5,phi,'--k') # one sieve larger\n", "for p,df in d.iteritems():\n", " plot(p*ones_like(df.geom_mean),Phi(df.geom_mean),'o',color='0.8')\n", "plot(phi,Phi(df_stats.Mg),'or')\n", "xlim([-.5,4.5])\n", "xlabel('Sieve size (phi)')\n", "ylabel('DGS geometric mean (phi)')\n", "title('MOF Poking Eyeball Tests, May 3, 2013')\n", "grid()\n", " " ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEXCAYAAABF40RQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXlYE/fWx78B3LAg4FKsqCxW3EpZRKCKIr2vqBSwLlgR\nFJcWN7witAWXamuvYt1FoQq2KoKVClfADbQmVgFBqigi2FYLqLUq4IIbhGTeP7yZJEJIAslMQn6f\n5+F5nP3kODkzOSuHoigKBAKBQNAp9NgWgEAgEAjMQ4w/gUAg6CDE+BMIBIIOQow/gUAg6CDE+BMI\nBIIOQow/gUAg6CDE+BOaJTg4GCtXrmxy2/z58/Htt98yLFHLKS8vh56eHoRCYYuO19PTw61btwA0\nrxcCQRsgxp9hLC0t0aFDB1RXV0utd3BwgJ6eHiorK+l1JSUlGD9+PMzNzdG7d2/4+/vjzp079HYe\njwc9PT0YGRnRf35+fk1eNzg4GB06dICRkRH69euHhQsXory8XK68HA4HHA6nyW1xcXFYsWKFAp9a\nOURGWvJzGRkZ4eeff1b5tVqKLL2sXbuWlrdTp04wMDCgl9977z2lr8Pj8dC7d29ViIzVq1dDT08P\n27dvl1q/bds26Onp4euvv1bJdST54osv0KdPH3Tp0gVubm74z3/+o/Cx+/btw9ChQ2FiYoKRI0c2\nkvv58+cIDAxE9+7d8d577+Ho0aP0Nj6fj8mTJ8PKygp6eno4e/as1LFbtmyBjY0NjI2N4eDggPDw\ncAgEgtZ9WC2DGH+G4XA4sLa2xsGDB+l1xcXFePnypZQxuXLlClxdXdGzZ0+cPXsWqampEAgEcHBw\nwP379+n9evXqhdraWvovPT1d5nW//PJL1NbW4uzZsygvL0d0dLRCMrNVB/jkyROpzzZlyhRW5JBF\nU3pZtmwZLe/333+PDz74gF4uLi5mQUoxHA4H/fv3x/79+6XW79u3D7a2tjIf8q1hzpw5uH79Op48\neYIff/wR33//PU6cOKHQsS9fvsS2bdtQVVWF7du3Y8eOHTh58iS9ffHixaisrEROTg4+++wzBAQE\n4K+//qK3jxw5EgcOHIC5uXmjz+bn54fCwkI8ffoU//3vf8Hj8bB7927VfGgtgRh/FggMDJT6Au7b\ntw8zZsyQMiYREREYNWoU9uzZA1tbWwwbNgypqano3r07vvrqqxZdV3T+Xr16ITAwkH5Tunv3Llas\nWAFLS0vMmjULly9fbvL42tpajB49GkuWLAEg7frg8XiwsLDA7t27YW1tjeHDh0t9yZ89e4a1a9fC\nwsICnp6eiI6Ohru7u9Kf4eLFizA3N5fSVVpaGuzt7enlEydOwNfXF7a2ttiyZQuePXsmdY6ff/4Z\ntra2cHFxkTImBQUFcHNzg6mpKdzc3LBjxw40NDQoLaMIiqKk5Lx37x7WrFmDfv36YerUqcjPz6e3\nnT9/HuPGjYOZmRksLCywefNmvHjxAuPGjcPff/8NIyMjGBsb459//sH169cxceJE9OjRA+bm5ggP\nD1dYJmdnZ7x48QLXr18H8PrXZV1dHYYOHUrL+ujRI3z00Ufo0aMH3n33XXz11Vd48OABrbuhQ4dK\nnXPz5s2YMGFCk9eztbXFW2+9RZ/bwMAAnTp1UkjWefPmYfjw4TAwMIC9vT1mzpyJvXv3Anj9YPjp\np5+wbt069O/fH6GhoXBzc6O3t2vXDosXL8bw4cOhr6/f6NzW1tYwNTWllw0MDGBoaKiQXG0FYvxZ\nwNXVFU+fPkVZWRkEAgEOHTqEwMBAertQKEROTg68vb0bHevt7d3oJ6yyVFZWYv/+/Zg4cSIAICAg\nAE+ePMHFixfh7u4ODw8PvHr1it6fw+GguroaH374Idzd3bF161Z6veQb1YMHD1BUVIT8/HzMmjUL\nn332Gb1t9erVOH/+PM6fP4+lS5di+/btct80m3qzdnZ2RteuXZGVlUWvS0xMxMyZMwEAGRkZiIyM\nxLJly8Dj8XDhwgWsW7dO6hwHDhzAyZMnERUVhYCAANy4cQPAawMgetPctGkTNm7cCB6Pp4hKFcLb\n2xsGBgYoLCzEjBkzMG7cODx//hwAEBYWhpCQEFRVVaGkpASjR4+GoaEhTp48iXfeeQe1tbV4+vQp\nzM3NsWrVKowePRp3797FrVu34O/vr5QcQUFB9MvHvn37EBQUJLWdoijMmTMHlZWVOHnyJAoKCmiX\ni6+vL/766y+UlZXR+0vqvymio6Px1ltvYdCgQVi5ciU8PDyUkldEXl4e3n33XQCv7+GXL19KudLe\ne+89KbnkkZycDCMjI1hbW2Ps2LHNfoa2CDH+LCH6Ap46dQqDBg1Cr1696G1VVVV49eoV3NzcGh3n\n6uqK27dv08t///03TE1N6b/Dhw83eT2KorBx40aYmppi/PjxGDJkCCIjI1FVVYWLFy8iOjoa3bt3\nx+zZs/Hee+9JvbXfvXsXHh4emDp1Kr755ptG5xUhFArxzTffoHv37ggODsbjx49pw3rixAlERETA\n0tISH330Ef71r3/JdSd169ZN6rOJzjVjxgwcOHAAAFBTU4Ps7GwEBAQAAA4dOoQvv/ySdplFRUXh\nyJEjUucNCwuDlZUVJkyYgLFjx9K/gBwdHTFs2DDo6+vjgw8+QGBgoEw3mrL88ccfePHiBaKiomBi\nYgJvb2+MGjUKx48fp3V38+ZNPH36FF26dIGDgwOAph+AQqEQlZWVqKmpgaGhIVxcXBSSQXSuwMBA\nHDx4EA0NDY1ePADAzMwMH3/8MTp27AgbGxtERETQeujQoQP8/f1p/ZeUlKCiogIfffSRzOtGRkbi\nyZMnSE9Px4oVK3D69GmF5JVk9+7dKC4uxhdffAEAqK6uhomJCYyNjel9rKysGsXSmiMgIAC1tbU4\nffo0Dh06RP9q0BWI8WcBDoeDoKAgJCUlNeny6datGzp27Ijc3NxGx+bl5aFPnz708jvvvINHjx7R\nf5MnT5Z5zc8//xyPHj3CtWvXsGHDBrzzzju4cOECrK2t0blzZ3rfoUOH4vz58wBeG4xjx47h1atX\nCAkJafZz9ezZE926dQPw+i26W7duuHv3Lp4+fYrS0lLaoAGvDa08qqurpT6bra0tAGD69OnIzMzE\nixcvkJKSgpEjR+Ltt98GAJw+fRrz58+nHxijR49GeXk57bYAIOUicnBwQF5eHoDXD7l58+bBzs4O\nxsbG2LJlC65evSpXTkU4ffo0/vrrL6mH2S+//IJz584BAPbv348rV67A2toaU6ZMwZUrV2Sea8uW\nLXjx4gWGDBmCsWPHKvVLkMPhoHfv3ujXrx+ioqLQv39/WFhYSO0jFAqxbNkyuLu7w8TEBJMmTcL1\n69fpe3TmzJlITk4G8Pqtf+rUqWjXrl2z1zUwMICPjw8++eQTqXiXIvz3v//FypUrceLECRgZGQEA\nunbtisePH+Pp06f0frdu3aLvP2Xw9PTEggULkJiYqPSx2gwx/izRp08fWFtb48SJE7T7RYSenh6G\nDx+OY8eONTru2LFjGDlyZIuu2dRbpKurK27dukW7HwDQ7h/gtbH49NNP4eXlhfHjx+PFixdSxysS\nJDQ2NsaAAQOkYgmXLl1q0WcAAAsLC7i6uiItLQ0HDhyQclt4enoiPj5e6qHx/Plz9OjRg97nTTk+\n+OADAMC3334LPp+P48eP48mTJwgLC2txWuibeHp6wsbGRkqup0+f0u6UwYMHY//+/bh37x7ee+89\nzJ07FwCgr6/f6P+tT58+2LlzJ/755x/4+/tj2rRpCsspOteMGTOwefNmzJgxg94m+r/8+eefcezY\nMfz444+oqqpCamqqVPzC1dUV7du3x6+//oqDBw82chs1x/Pnz9GzZ0+F98/KykJISAiOHTuGQYMG\n0ev79OmDTp06ST2ci4uLMWDAAIXP3Rq52gLE+LPInj17cObMmSYDYBs3bsSvv/6KOXPmoKysDAUF\nBZg4cSIePnzYopQ8WS6Wbt26wdnZGcuWLcODBw+wd+9elJSUwMvLS+q4HTt2wNbWFj4+PnQ84M2A\nZnOMHz8emzZtQkVFBY4fP45ffvmlRT5/ETNmzMD69etx7do1qYdnUFAQvvvuO5w/fx4CgQAPHz5E\nRkaG1LHbt2/HX3/9hczMTGRnZ9Mui7///htmZmbo2rUreDxeo6wYRWVrClHgc+PGjfjnn3/A5/Nx\n8eJFlJWVgc/nIykpCU+ePAEAGBoa0m+477//PqqqqnDv3j36XAcOHMDDhw9BURQ6d+6Mt956i95m\naWnZrNwipk6dilOnTtEZVJL/l3///TdMTEzQrVs3/P7771i/fn2j44OCgrBo0SK0b9+efni+CUVR\n2LVrFx4/foznz58jLS0NP//8M+bMmaOQvFwuFwEBAUhLS2sUZO7UqROmTZuGFStW4I8//kBMTAwu\nXLiA4OBgep+6ujr6XpX8NwAkJCTg4cOHqK+vx8mTJxEfH08/cHUFYvxZxNraWsr9IWkM33//feTl\n5eHevXsYNWoUPv74YxgYGODy5cswNzdv8pjmaC5fPykpCYaGhnB2dgaPx8Mvv/xCP5Akj9u9ezcs\nLCwwYcIE1NXVNTpnc7KsWrUKbm5u+OCDD7Bx40bMnj1byl/bFCYmJlJ5/qJAMwBMnDgRlZWVtG9a\nxLhx4/DNN99gx44d6N69O9zc3FBQUCAlY0BAAMaOHYtvv/0WiYmJ6N+/P4DXQemioiJYWFhgw4YN\nWLRokczP15w+Ze1z5MgR8Pl8fPjhh3Q8or6+HsBrg25lZYU+ffogPz+f/qzGxsb44osvMHLkSJiZ\nmeHevXvIysrCkCFD8Pbbb+PAgQPYvXs39PT0UF9fj5qaGri6usqVp2PHjvD09KR1J7lt9uzZ6NWr\nF/r374+goCDMnj270WcNCgpCSUlJo3jBmxw5cgQ2NjawsbFBRkYG/TkByJV3zZo1qK2txbhx4+h7\nQDIJYvv27bCwsICbmxt2796NgwcPwtLSkt5ua2sLQ0ND/P333/Dy8kLnzp3pOprc3Fy89957MDc3\nxw8//ID//Oc/LQ5EayscMsyFwAZTpkyBm5sbli5d2qLjKYpC//79sWvXLnh6eqpYOu0kJycHsbGx\nSEpKUvu1Xrx4AXNzc1y+fBk2NjYtOgeT8hIao/Y3f1Fhko+PT5Pbo6KiYG1tDScnJ6XStAjaxY0b\nN3D16lXU1dUhKSkJp06dwpgxY1p8vuTkZJiYmBDDL8Hw4cMZM6Rbt27F//3f/7XY8APMyktojIG6\nL7Bt2zYMGjQItbW1jbYVFBTg3LlzKCwsRFZWFiIiIqRKtAlth9raWkybNg337t2Dp6cn9u/fjyFD\nhrToXB4eHvjnn3+wZ88eFUtJUARLS0uYmZnhp59+YlsUQitQq9vnzp07CA4OxvLly7F582ZkZmZK\nbY+JiYFAIKArRm1sbHDz5k11iUMgEAiE/6HWN/+wsDBs2LBBKhdXkoKCAqk0se7du+PmzZtSPyXV\n0W+EQCAQdIHm3u3V5vM/evQoevToAQcHB5kCNJUq2JSxF+3H1t+qVatYl0FT/oguiC40WRdlZWVo\naGhg5doCgQBffvklTE1Nce3aNdZ1IQ+1Gf/c3FxkZGTAysoK06ZNw5kzZ6QKSgDAxcWFbjAFAA8f\nPoS1tbW6RGoxirQ+1hWILsQQXYjRFF3Y2to22chN3Tx79gwTJ05EXl4evLy8MHjwYMZlUBa1Gf+1\na9fi9u3b+Ouvv/DTTz/RQT5JXFxckJqaiurqaiQnJ2PgwIHqEodAIBDUQkVFBYYPH45u3brh1KlT\n6NChA9siKYTas31EiNw5u3btAgCEhIRg2LBhGDFiBIYOHQozMzO6WZSmIVk1qOsQXYghuhDDtC74\nfD7Onj2Lf/3rX4xe903y8vIwadIkfP7551iyZAk4HI7W3BcaX+TF4XAU8l8RCATdoLq6GlOmTIGR\nkRH++9//Qk+PnUYFiYmJCA8Px759+zBu3DhWZGgOebaTtHdQAFX2dNd2iC7EEF2IYUoXZWVlcHV1\nhZOTE9LS0lgx/EKhEJGRkVi9ejW4XG4jw68t9wVjbh8CgUBoDdnZ2QgMDMR3333Hmmvl2bNnCAwM\nxKNHj5Cfn9+iFtKaAnH7EAgEjefly5f0aM0RI0awIkNFRQV8fX3h7OyM2NhYtG/fnhU5FEWe7STG\nn0AgaAVCoZA1/35TgV1Nh/j8VYC2+PCYgOhCDNGFGCZ0wWZg18/PD3v27EFYWJhcw68t9wXx+RMI\nBEITiMZZ/vzzz+ByuVpRuKUMxO1DIBA0iuzsbBQXFyM8PJw1GSQDu6mpqVoZ2CVuHwKBoBVQFIWY\nmBjMnDkTLi4urMnxZsWuNhp+RSDGXwG0xYfHBEQXYoguxLRWF3w+H/Pnz8fu3buRm5vLWkZPXl4e\n3NzcEBwcjPj4+BZl9GjLfUF8/gQCgVVqamowZcoUdOrUCTk5OXJnO6sLTa/YVTXE508gEFjl3r17\nSEhIwLJly1jpyCkZ2M3IyGgzgV2S508gEAgyaAuBXVmQgK8K0BYfHhMQXYghuhCjjbpQV2BXW3RB\njD+BQGAMPp8PgUDAthgqCexqO8TtQyAQGEHUijkwMBCzZ89mTQ5dCezKs50k24dAIKidsrIy+Pj4\nYMKECZg5cyYrMgiFQixfvhwpKSltsmJXWYjbRwG0xYfHBEQXYoguxDSni+zsbIwcORLLli3Dhg0b\nWJ2xm5ubi/z8fLUafm25L4jxJxAIaiMtLQ0zZ85EamoqZs2axYoMulKxqyzE508gENRGVVUVamtr\nYWVlxcr1tbEVs6pgLdXz1atXcHFxgb29PVxdXbFly5ZG+/B4PHTp0gUODg5wcHDAt99+qy5xCAQC\nC3Tr1o01w69sK2adg1Ijz58/pyiKol69ekUNHjyY+uOPP6S2c7lcysfHp9lzqFlEheByuWyLoDEQ\nXYghuhCjSboQCARUZGQkZW1tTV27do3x62uKLuTZTrVm+xgaGgJ4HWxpaGhAhw4dmnr4qFMEAoHA\nEBcuXICTkxPatWvHmgxtacauulGr8RcKhXBwcEBJSQm2bt2K3r17S23ncDjIzc2Fvb09PD09sXDh\nQtjY2DQ6T3BwMCwtLQEAJiYmsLe3h4eHBwBxZF3dyyKYup6mLovWaYo8bC57eHholDxsLVMUheLi\nYqxbtw7r169Hnz59WJGnoqICH374IWxtbXHq1Cm0b9+eNf2IYPL6PB4Pe/fuBQDaXjYHIwHf8vJy\njB8/HklJSXBwcKDX19bWQl9fH+3atcO+fftw5MgRHD16VFpAEvAlEDQWPp+PRYsWITc3FxkZGSSw\nq0FoRG8fS0tLjB8/Hvn5+VLrjYyMYGhoiHbt2mHOnDm4ePEi6urqmBBJKd58musyRBdidF0X1dXV\n8PLywt27dxEdHU0Cu/9DW+4LtRn/qqoqPH78GMDrmyQ7Oxt+fn5S+9y/f59+MmVmZsLOzq7JuACB\nQNA8li9fDicnJ6Snp6Nz586MX18oFCIqKgqrV68Gl8tt060a1IHa3D7FxcWYOXMmBAIBzM3NMX36\ndMyYMQO7du0CAISEhGDnzp2Ii4uDgYEB7OzsEBERATs7O2kBiduHQNBI+Hw+a8FdTWzFzOVykZGR\nAT09PQiFQvj6+mL06NGsyUP6+RMIhDZFRUUFfH194ezsjNjYWI3oyMnlcpGSkoK5c+fS6xISEuDv\n78/aA0AjfP7ajrb48JiA6EIM0YUYpnShqa2YMzIyaMP/22+/AQDmzp2LzMxMNsVqFtLVk0AgNEtN\nTQ2+/fZbrF27Fh07dmRNDk1uxaynp4dLly7hwoULePDgAfLz8+Hq6sp68Lk5iNuHQCDIRLIVc3R0\nNGszdkWtmDV1xm5AQABevHiBrl27wsDAAA0NDaiurkbnzp2RlJTEikyknz+BQGgR2dnZCAwMxPr1\n62V25FR3kFNbKnZra2vRo0cPLFiwgF4XGxuLhw8fsihV8xCfvwIQ364YogsxbVUXFEUhJiZGbitm\nUZAzMDAQtra2CAwMpAelqAJtasUsEAhowy/y+S9YsAANDQ1sitUs5M2fQCA04vbt28jNzW22cEsy\nyCli7ty5SEpKavXbv7ZV7LZv3x5F586h6NAh3K+pQbGZGeynTtWYgHRTEOOvAJJ9bXQdogsxbVUX\nHA4H3333ndz99PTEjgMnJyep41uDKLC7d+9ejB8/vlXnYorHt2+jdNMmxNy5Q6+LuHMHT95+m0Wp\nmocYfwKB0CKEQmGT61uaoCEUCrFs2TL8/PPPWjdjt/2tW9j4v44GIjbeuYMxz5+zJJF8iPFXAJ5E\nF0tdh+hCTFvRRUNDAwwMlDcFvr6+WLduHbp06YIHDx6gR48eePz4MRYtWqT0ubQlsCsLY4ksKB4A\nD9F6Pc0Nq2quZAQCQa2IArtv9txShrfeegsLFiyAj48PFixYACMjI6XPoU2BXVnIer9/rsGxCpLn\nTyDoIHw+H6GhoTh//jwyMzNb1JEzLCwMgYGBjdYnJSVh8+bNCp1D2wK7slj66ad4cvAg9ki4eWZ3\n7gyTadOwOT6eFZlInj+BQJCipqYGU6ZMQadOnZCbmwtjY+MWnUdPhktDUQOujYFdWWyOj8dSAOOO\nHEFnisJzDgcDJ0xgzfArAjH+CtBWfLuqgOhCjDbqoqqqCm5ubiqp2JUM+P722290xo+8X+rqCOxq\nQkfNzfHxQHy81twXco1/bW0tSkpKcOPGDejp6aF///4YNGhQi3x7BAKBXbp27YqEhASMGjWq1ed6\nM+Cbn58vN+CrjsCurI6aAFhtqazpyPT5X7x4ETt27MClS5dgZWUFGxsbUBSFW7du4datW3ByckJo\naCiGDh2qXgGJz59A0Ei4XC727duH0NBQep2oMrgpo6uuVsyqiD20RVrs8//xxx+xaNEiODs7N7m9\noKAAP/zwg9qNP4FA0EwyMjKkDD8AhIaGNlnhq87AbmtjD7qKzFTP2NhYmYYfAIYNG4bY2Fi1CKVp\ntNUeLi2B6EKMpuuipqYG5eXlaju/pNEV9bMBGhtd0YzdhIQEtczYVXWxWWvR9PtChMw3/wMHDiAw\nMBCbNm1qtI3D4WDp0qVqFYxAILQcUSvm+fPnq+27Ks/oMlWx27dvX6xbtw5RUVH0urVr18Ld3V0t\n12sryDT+z/+Xr1pbWyv1pKYoSud+TmlD5J4piC7EaKouFGnFrAp8fX2RkJCAuXPn0pk+8fHxmDp1\nKqMVuxUVFfDy8kJcXBz09fUhEAgwduxYlJaWqu2azaGp98WbyDT+ISEhAIDVq1e36MSvXr3CqFGj\nUFdXh44dO2Lq1KkICwtrtF9UVBQOHToEU1NTJCUlYcCAAS26HoGg64gqdtetW4fU1FS1v/mK/PpJ\nSUl0cHHq1KmwtrbG8OHD4ezsjJSUFLV3tnz69CkcHR3h6Ogotf7ChQtqve6baEK6qTLITfWsqanB\n0aNHkZeXh1evXgF47fb54Ycfmj2uY8eO4HK5MDQ0RF1dHZycnODj44N+/frR+xQUFODcuXMoLCxE\nVlYWIiIicPTo0VZ+JNXDdt6uJt1UbOtCk9A0XfB4POzevVtuK2ZVMnr0aIwePZrWhWjGLpMVu/fu\n3VNqvTqQTDcV1TxoerqpXOO/aNEidO7cGZ6enmjXrh0AxaPohoaGAF7n9jY0NKBDhw5S2/Pz8zF5\n8mSYmZlh2rRpWLFiRZPnCQ4OhqWlJQDAxMQE9vb29JdOFFxR53JRURGj15Nc3rJlC3g8Hr766isA\nrwNrW7duBQD6S8ekPEVFRYxejywrvuzh4YFNmzahoqKCNv5MXR94HdgNDQ1FZGQk/SufievX19cj\nNjYWCxYsoAPPFy5cgKmpKWOfXzTb4LfffsONGzfg5OSEuXPnIjo6GhwOh5H/fx6Ph7179wIAbS+b\nQ25vn8GDB6OkpETuiZpCKBTCwcEBJSUl2Lp1a6Pij6CgIAQFBWHMmDEAAFdXVyQlJcHGxkYsoI7n\n+ZMcZoKmIwrspqSkIDMzk/FWzDNnzsSIESOQn59P+/xdXFyQk5NDG0N1Ex4ejoCAgEbrk5OTm0ya\nYYJW9/b55JNPsGfPHkyfPh0dO3ZU6uJ6enq4cuUKysvLMX78eAwfPhwODg70doqiGgmna8FkeZAc\nZoImIxnYLSgoYKUjJ5/Px4ULF6Tm5+7cuRN8Pp8xGTQt3VQRZOb5v/XWWzAyMsJ3332HTz/9FMbG\nxjAyMoKRkZHSjaAsLS0xfvx45OfnS613cXHB9evX6eWHDx/C2tpayY+gfiR/2jKNpt1UbOpC02BL\nFxRFITY2Vuq7wwYVFRVwGzIEKCyE5ePH2Dp9On49doxxOXr27AlXV1fExcVh9+7diIuLg5ubG8zN\nzRmTQZT5BIhrHuLj4+Hj48OYDMoi883/2bNnrTpxVVUVDAwMYGJigurqamRnZyM8PFxqHxcXFyxd\nuhQzZsxAVlYWBg4c2KprtkUk0+lEiNLpCLqHqBVzTk4Oxo0bx5oceXl5+MTbGx/p6WFndTV4d+/C\n4+pVLL95EwAw0tubMVmEQiH0nj9Hp5IStOfzUd+uHfSGDGH0BUky8+nOnTsoKyvD1KlTNTbYCyjY\nz/+3337DyZMnweFwMHbs2EYpVU1RXFyMmTNnQiAQwNzcHNOnT8eMGTOwa9cuAOJU0sjISBw6dAhm\nZmY4cOBAoweArvv8gdeZBJmZmbQufHx8NPqmIqiH6upqTJkyBYaGhkhOTsZvv/3GShaYqBXzxxYW\n2HX5cqPtK728sObkSbXLIWLHunW4Eh2N+KdP6XWfGhvj/chILJIo/NI15NpOSg5bt26l3NzcqA0b\nNlAbNmyghg8fTm3dulXeYSpDAREJhDZPaWkp1a9fPyoiIoJqaGigzpw5Q82bN48qLCyk/+bNm0ed\nOXNGbTIIBALqyy+/pKytralr165Rq0aNoiig0d+qUaPUJkNTLB8zpkk5Vnh5MSqHpiHPdsod4xgf\nH4+TJ08iIiICEREROH78OOI1eECBOiB+bjFEF2KY1MXFixcRFRWFDRs2QF9fn04tlGTu3LnIzMxU\ny/WfPXtfT7mpAAAgAElEQVSGiRMnIi8vD/n5+Rg8eDAaJFK3eRL7CpRMDGktBnV1Ta7X/19dEtNo\ny3dErvHv168f/vjjD3r55s2bUoVaBAJB/QQFBWH27Nn0MpNZYLJm7I5ZvBjLJdKyAWCZjQ3+741O\nn+rmwfOmJ+g+ePGCUTm0Dbmpnk+fPsWwYcPo3N2SkhJ4eHjAx8cHHA4HGRkZaheSbUQFFQSiC0nY\n1AVTWWDNtWIWBXVXxsRA/9Ur/NKxI8aGhjIa7AWAp++8gwgLC2y8c4deF25hgdp33mFUDhHa8h2R\na/xFlaWSiAIJJNecQFA9QqFQ5pu9CCaywBSZsTvS25txY/8mHbp1w8DwcCxOSUG7+nrw27fH+/7+\neFxczKpcmo5M4y8y7s09xVT9lqGp8DSshwubEF2Iey3dvXsXvXr1UmmWTVlZGT755BMcP34c7zTz\n5iqrqZoq5FCmFbM6daEo9+7dg/2CBbB/o5FdanY2o3KI0JbviEzj7+7uDg8PDwQEBMDW1pYe9NzQ\n0IAbN24gOTkZPB4POTk5jAlLILCNOht4SbZibs7wixA1VVMlyrRi1pRmZqampnRvHxE7d+6Eqakp\nYzJoIzLz/AUCATIyMhAfH4+rV69CX18fFEVBIBDAzs4On332Gfz8/OT+PG21gCTPn6BBqKPXEkVR\n2LFjB9auXYuUlBTWhpAoO2NXU/pOhYWFYdCgQY16+5SWlup0/6sW9/bR19fHxx9/jI8//hjA68Av\nh8OBkZGR6qUkELQEdWTZhIeH49SpU4y2Yn6TlszY1ZS+U76+vkhJScH8+fPpdaQKXj4Kv7aLevvo\nItqSt8sEuq4LySwbybm1rfl16ufnh5ycHNYMf0tn7KpDFy1h9OjR8Pf3R1JSEpKTk5GUlMRqawVt\n+Y7IzfYhEAhi1JFlM2rUKFWIpjStnbGrSX2n1BH/aOso1NuHTYjPn6BptIVeS5KB3dTU1Ba3Ym4L\numirtLqfP4FAkKalb5kUReH333+Hra2tGqRSHMnAripm7IrSwslLmnYh1+d/+vRpeHp6wsTEpMX9\n/LUdbfHhMQHRhRhldMHn8zF//nxMnz4dAoFAJdfncrkICwtDeHg4wsLCwOVy5R4jmrEbHByM+Pj4\nVhl+UapnYGAgbG1tERgYiJSUFIXkaMtoy3dE7pt/ZGQktm3bBjc3N7WndRIIbRHJVsxnzpyha2Za\ng2SOvQh5OfaKVOwqg6zmcklJScT1owXItebt27eHk5OTTht+bajWYwqiCzGK6KKsrAyurq5wdHRE\nenq6yn41K9PVUygUIjIyEqtXrwaXy1WJ4Qdep3+LcHJyov/95MkTlZxfW9GW74jcN393d3dMmDAB\nU6ZMgYmJCYDXgYSJEyeqXTgCQZupq6vDRx99hOXLl2PWrFkqPbeiOfbKVOwqy71795RaT9As5L7O\n379/H+bm5jh//jyOHj2Ko0ePqq1nuKaiLT48JiC6ECNPFx06dMDly5dVbvgBxbp6ymrFrCpEbRUA\ncZ4/W20VWhL/UBfa8h2R++a/d+9eBsQgENom6iqMlJdj35KKXWXp0aMHBg0ahLi4ONy/fx8FBQVw\nc3NDaWmpyq/VHC2JfxAUyPOvr68Hl8tFVlYWHj16RN9EP/zwQ7Mnvn37NmbMmIEHDx6ge/fu+Oyz\nzxAQECC1D4/Hg5+fH6ytrQEAkyZNwooVK6QFJClkBEKTyMqxFwV29+3bp9Yh700ZXdEDiEmjqyk9\nhjSNVuf5r1ixAi9evMDx48fx73//GwcPHlToP7Zdu3bYsmUL7O3tUVVVhWHDhsHHx6fRm9CoUaN0\nYiAMoW1TVlaGS5cuNXrBUTeSOfaiwG5LK3aVRZ1tpZVBU3oMaRtyff6//PILYmJi0KlTJ/z73//G\niRMn8Msvv8g9sbm5Oezt7QEA3bp1w+DBg1FYWNhoP214q9cWHx4TEF2IEekiOzsbI0eORH19PWPX\nlsyxDwgIQGBgIHbs2IEzZ87QM3aZYPTo0di8eTN8fHywefNmVtwsTE01UxRt+Y7IffPX19cHh8OB\ng4MDTp06hXfffRcvlJyN+eeff6KkpATDhg2TWs/hcJCbmwt7e3t4enpi4cKFsHljJigABAcHw9LS\nEgBgYmICe3t7Op1KpGh1LhcVFTF6PU1eLioq0ih52FymKAqhoaFITk7GkSNH4O7uztj1RameokCr\nk5MTVqxYgejoaFy7do1xfYhg4/+jT58+dPxDpI9Lly5h6tSprMjDlr3g8Xh0jFZkL5tDrs8/Pj4e\nkyZNwq1btxAREYG///4ba9asUbh5U21tLTw8PPDVV1/Bz8+v0TZ9fX20a9cO+/btw5EjR3D06FFp\nAYnPn6BhcLlcpKen4+bNm7h//z7CwsIwbdo0RmUIDw9v0sWUnJyMTZs2MSqLJkB6DDVGnu1Ua2M3\nPp8Pb29vjB8/HkuWLGl2X4qiYG5ujsrKSnTo0EEsIDH+BA1CVmaJv78/CXISNAp5tlOuz7+6uhrr\n16+Hr68vAOD69evYs2eP3AtTFIU5c+ZgyJAhMg3//fv3aeEyMzNhZ2cnZfg1hTd/2uoyuq4Lycpa\nkYtBVmWtuhAKhaipqcHKlSul1sfHx8PHx4cxOSTR9ftCEm3RhVyf/1dffYXBgwejvLwcAPDuu+/C\n398fc+bMafa4nJwcHDhwAHZ2dnBwcAAArF27FpWVlQCAkJAQHD58GHFxcTAwMICdnZ1O/lwlaBds\nZ5ZIVux6eXlh5cqVaN++Perr6zFmzBidd3UQFEeu28fFxQX5+flwcHDA5cuXQVEU7O3tceXKFWYE\nJG4fAsuI7j8Oh8Oqu0WyFfPUqVORlpbGuvuJoLm02u3j6OiI27dv08tpaWmsDZgmEJiGz+dj3rx5\ntKtTVFkrCRPuljdbMR8/flzhxm4EQlPIdfssWbIECxcuRHl5Ofr16wcrKyu6n4euwOPx6NQqXUeX\ndCHZinnDhg0AXue1X7lyBStXrsTTp09hbGysdndLU62Y9fT0UHTuHIoOHUJ7Ph/17drBfupU1gqb\ndOm+kIe26EKu8be1tUVGRgYePnwIgUAAc3NzJuQiEFiltLQUvr6+mDBhAqKjo+ke/FwuFzdu3MCa\nNWvw22+/wcnJCQkJCeByuSp/AAiFQixfvpwekCJZuHX/zz/BSUlBzJ079LqIO3fwQKK1MoHQHHJ9\n/hRF4dKlS8jLy0N9fT1dTr506VJmBCQ+fwLD/Prrr5g8eTLWr1/fqCMnUz5/eTN2Fzg7I7aJivkF\nzs6ILShQmRwE7aXVvX0WLFiAP/74Ay4uLq2e9UkgtBYul4uMjAzo6elBKBTC19dX5W/clpaWSEtL\nw4gRIxptYyLbR5EZuz06d27y2B6GhiqTg9C2kWv8z507h+LiYp1ukqQtPjwmYFMXTLXu7dOnD/r0\n6dPkNqFQSPva79fU4G0zM9hPnaqyX6eKtmJukFEPI+jYUSVyKAv5jojRFl3Izfb5+OOPsX//ftTV\n1TEhD4EgE2VGF6oLmx49cGHVKsRcuIAFv/+OmAsXcGHVKlh3797qcycmJsLPzw8JCQkICwtr9oVr\nzOLFWP5GH6xlNjb4v9DQVstB0A0UGuM4adIkfPrpp3T1LYfDkZrf2dbRhqc4U7CpC3W4XCorK2Fh\nYaHwjOp/eDzE/+/e9/jfuvinT7Hy7FkgKqpFMogCu4cOHVK4FfNIb28AwMqYGOi/egVBx44YGxpK\nr2ca8h0Roy26kGv858+fj6NHj8LNzY34/AmsourWvdnZ2QgMDMTJkyfh6Oio0DEGMn4B67961SIZ\nJAO7BQUFSo1aHOntzZqxJ2g/cl93evXqBUdHR502/NrSq4MJ2NSFr68vYmJipNZt375d6QIriqKw\nfft2zJw5E2lpaQobfkDa186TWN8SX7u6Z+wyCfmOiNEWXch98x84cCDc3d3h6+uLLl26AACjqZ4E\ngiTPnj1DXFwc9PX1IRAI8Pz5c6WO5/P5WLRoEXJzc5GbmwsrKyuljh+zeDGW37yJ/9y8Sa9bZmOD\nsUr62pmYsUsgNIfcPP/Vq1e/3vGNm3PVqlVqE0oSkudPEKGKHPuwsDD88ccfSE5OhrGxcYvk2LFu\nHXi7dqGTQICX+vrwCAnBIiX8/U1V7BIIqqbVef4i408gsI0qAr5fffUVjI2N6YpdZeFyuSiprERU\naiq9TtEK35YEduXJou6aB0LbRbEUBx1HW3x4TMCmLlQR8DU1NW2x4Qda3s//2bNnmDhxInJyclBQ\nUKASw//mDF9RGwg2IN8RMdqiC2L8CVoDWx01JWnJr4/KykqMGDECXbt2xenTp1US2NWEmgeCdiPX\n7UPQnrxdJmBTFyKXRlJSEu3PnDp1apOuDj6fj/j4eISEhLTqTf9NJH99OEk0UZP160MU2I2IiJBb\nuKUMsupsnjx5opLzKwv5jojRFl3IffOfMWMGHj16RC/X1NRg9uzZahWKQJDF6NGjsXnzZmzatAmb\nN29u0vBXV1fDy8sLx48fx6sW5t/LQplfH5IVu0uXLlVpRs+9e/eUWk8gvIlc43/16lWYmprSy2Zm\nZrSvU1fQFh8eE2i6LsrKyuDq6gonJyekp6ejs4wGaC1l9OjR8Pf3R1JSEtavX4+kpKRGvz6EQiGi\noqKwatUqcLlctWT0mJqaNpqrsXPnTqnvKpNo+n3BJNqiC7lun759++KPP/7Au+++CwD4/fffYWFh\noXbBCARlEVXsNtWKWdWI3DxvuntEFbs1NTVKV+wqQ48ePTBo0CCpmgc3NzeUlpaq5XqEtofcPP+s\nrCwsXLgQ//rXv0BRFE6fPo24uDiMGTOm2RPfvn0bM2bMwIMHD9C9e3d89tlnCAgIaLRfVFQUDh06\nBFNTUyQlJWHAgAHSApI8f4ICUBSF6dOnY8GCBU22YlYVsjqL+vv7w8bGBr6+vnByckJcXJxaq+Kb\nkiM+Pl5mDISge8iznXKNPwC8ePECx44dAwCMHz9eoZ/S//zzD/755x/Y29ujqqoKw4YNw5UrV2Bk\nZETvU1BQgKVLlyIjIwNZWVlISkrC0aNHlfoABAKTyCo027lzJ06ePKnywG5zcLlcZGZm0t8RHx8f\nYvgJNC0u8hLNJ62pqQEAeHp6AgDq6upQV1cHMzOzZi9sbm5Oj3zs1q0bBg8ejMLCQqmbMz8/H5Mn\nT4aZmRmmTZuGFStWKP7JGERb+nMzga7rQjLVUzTGEXg99jE+Ph7eDDZaGz16tMYYe12/LyTRFl3I\nNP7Tpk3DsWPH4Ojo2OgthsPh4NatWwpf5M8//0RJSQmGDRsmtb6goABBQUH0cvfu3XHz5k3YvNGn\nPDg4GJaWlgAAExMT2Nvb08oVBVfUuVxUVMTo9TR5uaioSKPkYXr52rVrUkZflPzQp08feHt7sy4f\nW8siNEUeNpfZshc8Hg979+4FANpeNkezbh+KonD79m2ZU40Uoba2Fh4eHvjqq6/g5+cntS0wMBBB\nQUHw8vICALi6uiI5ORnW1tZiAYnbh/AGWVlZ+PHHH3Hw4EHGG6IFBASgU6dOWLBgAb1u586dqKur\nQ1JSEqOykPYOhOZodW8fb29vFBcXt+jifD4fkyZNQlBQUCPDDwAuLi64fv06bfwfPnwoZfgJmoMm\nGBqKohATE4N169YhJSWFlU6YPXv2xIABAxpl2ZSVlTEqB1MjLQltl2bz/DkcDtzc3JCenq70iSmK\nwpw5czBkyBAsWbKkyX1cXFyQmpqK6upqJCcnY+DAgUpfhwne/Gmra0j2kbG1tWWljwyfz8e8efMQ\nt2EDXA0MsDsgAJMtLbFj3TrGZABeFzk6Ojpi/vz5cHJywvz58+Ho6Mj4r1NNa++g698RSbRFFwoN\ncE9ISEDXrl3pAC6Hw8HVq1ebPS4nJwcHDhyAnZ0dHBwcAABr165FZWUlACAkJATDhg3DiBEjMHTo\nUJiZmeHAgQOt/TwENZCRkQFHR0fExsbiwYMHyM/Ph6urKzIzMxl5y3z69CkmTJiAlw8eYMSTJ4iv\nraW3fRodjR2AUi2VW0piYiJOnjwJoVCIxYsX0+tFKZZMoo6RlgTdQm6qZ0VFRaO3Gg6Hg759+6pV\nMMlrEZ8/u0yfPh0dO3aU8nPHxsbi1atXjPi5BQIB9u7dixPffIPD/3t5kGRy3744XF6ututLtmLO\nzMzEgwcPWE+xVMVsA0LbptV5/kFBQUhMTJS7Tl0Q488+48ePx5o1axqtX7lyJY4fP86YHEG9eyPx\nzp3G6y0skHj7tlquKVmxm5aWpjGjFkmRF0Ee8myn3N4+165dk1p+8eIFrl+/3nrJtAht8eGpi549\ne9L/luzrJLmeCV7K6M4pa31rkdeKmc37QrLHUHJycpM9hphE178jkmiLLmT6/NeuXYt169bh5cuX\nUlW53bt3R6iS80oJ2o2scYeimc6qhM/no7a2tskiQo+QEHwaHY14iXbGnxobwyMkROVyqKsVsyrR\npCIvgvYh1+0TGRmJ6OhopuRpBHH7sA9TLobq6mpMmTIFLi4uWCcji6e183MVITExEUuXLsXevXsZ\nrdglEFRJq/P8o6OjkZ+fDy6Xi8jISFRWVuKff/5pVK1LaLsoM0SlpZSVlcHHxwcTJkzAt99+K3O/\nwa6uuPngAV1vMNjVVWUySAZ2eTxeq0ctEggaDSWH//znP9Qnn3xCDRgwgKIoiqqurqacnJzkHaYy\nFBBR7XC5XLZF0BjUoYusrCyqe/fu1A8//NDsfmfOnKHmzZtHFRYW0n/z5s2jzpw502oZamtrKT8/\nP8rd3Z16+PChQseQ+0IM0YUYTdGFPNspN+CbmZmJpKQkdOzYEcDrYS719fVqfiQRdIULFy5g5syZ\nSE1NlduDX12FTRUVFRg+fLhKZ+wSCJqOXLePhYWFlLEvLS1F//791SqUpiFqoqTLSLZ3SE9PV1l7\nh2HDhuHSpUsKZQ6po7CpNYFdcl+IIboQoy26kGv8Q0JC4OPjgwcPHmDWrFk4d+4c4uPjmZCNoCGo\ns4+Mnp6ewimjVVVVuHTpEi5cuAADAwM0NDTA1dUV1dXVLbp2YmIiwsPDsXfvXrWMWiQQNBmFh7mc\nOHECQqEQPj4+tAuICTQh24enJf251YVkNalkO2Omq0nHjRsHQ0NDLFu2jF63du1avHz5Uqliszcr\ndlsa2NX1+0ISogsxmqKLVhd5AUCnTp1ga2sLa2trXL9+HZcuXVKZgATNR1XullOnTtG9nVoCh8OR\nMvwAGi3L49mzZ5g4cSJycnJQUFBAMnoIOotct8/333+PtWvXonfv3lIzSZns6Mg2mvAUZxOhUEj/\nW/TWDzQeXi4LSqIVc3p6eovnQ8iaiavorNyKigr4+vpi6NChSElJafWMXV2/LyQhuhCjLbqQa/x3\n7NiB69ev46233mJCHoIG4uvri4SEhCaLvOTB5/OxaNEi5ObmIjc3F1ZWVi2WQ1ZFsYmJidxjtaFi\nl0BgErlunwEDBuD+/ftMyKKxaEuvDnUh2Udm/fr1CveRqa6uhpeXF+7evYucnJxWGX7g9TjPmJgY\nqXXbt2/HzJkzmz0uMTERfn5+SEhIwNKlS1Vm+HX9vpCE6EKMtuhCoQpfZ2dnDB48mH7D4nA4yMjI\nULtwBM1B1EdGmWDWzz//DCcnJ0RHR0NfBc3Xmqo0Dg4OlvkQkgzscrlc4t8nECSQm+3j4OCACRMm\nwM3NjfaRcjgcjBo1ihkBNSDbh6B9aGorZgKBKVrdz9/e3h6XLl2SmfGhbojxJyiLZGA3Li6u1YFd\nAkEbaXWq50cffYSQkBCcOXMGly5dov90AS6Xi7CwMPj7+yMsLEynMpxkoen+zLy8PLi5uWHmzJlI\nSEhQq+HXdF0wCdGFGG3RhVyf//nz58HhcBpNcmrrhlCyqlVU2KSqqta2RnV1NUJCQrBlyxb07t2b\nNTlIxS6BoDgKVfi2lNmzZ+PYsWPo0aMHiouLG23n8Xjw8/ODtbU1AGDSpElYsWKFtIAsuX3IjFTF\nKC0tha+vLyZMmKCywK6yqKpil0BoS7S6n39tbS0OHDhAZ/f4+fkhMDBQobz/WbNmITQ0FDNmzJC5\nz6hRozQyc0gdTcS0GcnGbkKhEL6+vqivr0dQUBDWr18vtyOnupAM7BYUFJDALoGgIHJ9/hs3bkRx\ncTG++eYbfP311yguLsaGDRsUOrm7uztMTU2b3UdTg7mSVa2Sc2s1VV51InKBBQYGwtbWFoGBgYiP\nj8fs2bMVasWsLthuxawtvl0mILoQoy26kPvmn56ejsLCQhgYvN7V0dERQ4cOxddff93qi3M4HOTm\n5sLe3h6enp5YuHAhbGxsGu0XHBwMS0tLAK+rOe3t7elcc5GiVb0sqmp1cHDAjRs34OTkhPj4eNja\n2krluqvr+pq0HBsbi8jISADAjRs3AADh4eHo0qULBAIBK/ro0KEDJk2ahAkTJmDKlCl0YFcT9KWL\nyyI0RR42l4uKili5Po/Hw969ewGAtpfNIdfnP3v2bHh5ecHf3x8URSE1NRUnTpzADz/8IPfkAFBe\nXg4fH58mff61tbXQ19dHu3btsG/fPhw5cgRHjx6VFpDFVE8ul4vMzExaBh8fH50M9oaHhyMgIKDR\n+uTkZGzatIlxeUSB3R9//JHM2CUQZNBqn39UVBS+/PJLREREAHg9fENVA92NjIzof8+ZMwfLly9H\nXV0dOnTooJLztxZRVauuI+kCk4TphzKp2CUQVIdc4//uu+8iLS0NfD4fFEWpNG/6/v376NGjBzgc\nDjIzM2FnZ6cxhl8SSbeGLuLp6YnvVq5Ej0ePcL+mBm+bmeG+iQnmf/EFYzJoYmBX1+8LSYguxGiL\nLuQa/+rqaiQkJCAnJwcZGRm4fv068vLyMGfOHLknnzZtGs6ePYuqqir07t0bX3/9Nfh8PoDXE8IO\nHz6MuLg4GBgYwM7OjhUXAkE2FEVhx44d2BkdjXEUhS337oEHwANA2DvvQP/FC0bkEFXsOjs7q6QV\nM4FAUMDnv3DhQgwePBjff/89rl69Cj6fDwcHB1y7do0ZAUl7B1bg8/kIDQ1FTk4OxpiaYtO5c432\nWenlhTUnT6pVDlEr5s8//xxLlizR2VRbAkFZWt3eobCwEAsWLKCLdwwMDFgp5CEwh6gV8507d5CT\nkwMjvaZvE/1Xr9Qqh2QrZtKDn0BQLXKNv6OjI27fvk0vp6Wlwd3dXa1CaRpvprO1ZRoaGjBq1Cg4\nOTkhPT0dxsbGaJCIw/Ak9hWoaZazUChEVFQUVq1aBS6Xq7GtGnTpvpAH0YUYbdGFXJ//kiVLsHDh\nQpSXl6Nfv36wsrJCbGwsE7Kxjqiq9e7du0hPT4evr2+bz/4xMDBAenq6VL2FuYcHPr1wAfFPn9Lr\nPjU2xvtqaOutiYFdAqEtonBvnwcPHkAoFMLc3FzdMknBls9fsrGbiISEBPj7+7f5B8CbhIWFYUjf\nvriSkoJ29fXgt2+P9/39UVJZqdI+R5KB3djYWBLYJRBaQavz/FNTUxv5Wq2srODg4NB66TSYjIwM\nKcMPAHPnzkVSUpLOGX89PT3Yu7vD/g133/XkZJVdgwR2CQRmkWv8Dxw4gJMnT8LV1RUAkJ+fDxcX\nF9y9exfbtm3DuHHj1C4kG0g2dhO1dAbYaezWVFM1VTyAqqurceXKFXh6eja735t9jkS6UNUvMm1t\nxawt+dxMQHQhRlt0ITfg+/LlS1y+fBlcLhdcLheXL19Ghw4dkJ2djW3btjEhIytoSlWrZFO1gIAA\nBAYGIiUlpdXzFMrKyuDq6oozZ87I3VfU50iS+Ph4+Pj4tEoGbQnsEghtEbk+fzs7u9fpfv9rxVBb\nW4sPPvgAxcXFsLOzw9WrV9UroAb5/OPj4zF16lRG3T4if3vRoUNoz+ejvl072E+d2ip/e3Z2NgID\nA5VqxazqPkeiwO6jR4+QmppKArsEgopptc9/7ty58Pb2xscffwzgdZfPOXPm4Pnz5zAzM1OdpBrG\n6NGjceXKFaxcuRLt27dHfX09xowZw7i/v7q8HKWHDyPmzh16XcSdO6gaOlTpc4kqdteuXYu0tDSM\nGDFC4WNV2eeIVOwSCOyjULbPxYsXkZWVBQ6HAy8vLzg5OTHm+9aEN3/JMY5MZ/tMtrTE4YqKxuv7\n9sXh8nKlzlVRUQF/f3/89NNPsLKyUupYybTXXr16tTju0JYCu9ri22UCogsxmqKLVr/5A4CzszOc\nnZ1VJpQ2oCnZPiYGTf8XyVrfHH379sWFCxeUNriqmmesrYFdAqEtIjfgq6tIZvuIslsA5rN9Hjc0\nKLVeHi2RX/JBKNLF3LlzkZmZqdDxbTWwqwlvd5oC0YUYbdGF8q+POkJVVVWT66urqxmVw9DBAREC\nATZK+PzDLSxgyGCdRWvmGUsGdknFLoGgOSj15n/r1i0IBAJ1yaJR8Pl8uo2FaIbvzp076ZbUTNHV\n0hIDw8Ox2M0N4U5OWOzmhkHh4ejWjM+eoihs374dSUlJKpGhpfOMKysrMWLECHTr1g2nTp1qc4Zf\nW3q4MAHRhRht0YVM4x8SEkK3bW5oaICrqys+/PBD9O/fH8ePH2dMQLbo2bMn3jY0xBo/PxyKisIa\nPz+Yd+7MeHsLX19fFN64gZkxMQjYtQszY2JwsaxMZo49n8/HvHnzsHv3bnzwwQcqk0HZPP+8vDy4\nurpi5syZiI+PJxk9BIKGITPbZ9CgQbh+/ToAYN++fdi/fz9++eUX/Pnnn1iyZEmjWbtqE5ClbJ9A\nPz+YX7ok5W6JsLDAfScnJB45wqgsiubYV1dXY8qUKTA0NERycjKMjY0ZlwEggV0CQROQazspGbz/\n/vv0v8eMGUP9+OOPTW5TN82IqFbmDx1KUUCjv/nOzqzII4/S0lKqX79+VEREBNXQ0MCKDAKBgIqM\njKSsra2pa9eusSIDgUB4jTzbKdPtY29vjx9++AG//vorLl26hMmTJwN47QJ6peYhHppAj86d6X/z\nJNJfMO0AABW+SURBVNcbGjIuiyKIhptv2LBBrcN2ZPkznz17hokTJyI3Nxf5+fk6MVxdW3y7TEB0\nIUZbdCEz2+frr7/Gd999h8TERMTHx+Ott94CAJw+fRre3t6MCcgWDTIGyatrgElzKNLYbdCgQRg0\naJDaZWhqtkFlZSV8fX0xdOhQUrFLIGgJCvfzV5bZs2fj2LFj6NGjB4qLi5vcJyoqCocOHYKpqSmS\nkpIwYMCAxgKy5PP/9dgx/Pezz7Dl77/pdUveeQcTd+/GSAYffpowV6A5GTp27NhmKnYJhLZEi2f4\nXrt2DampqfTykiVLMGvWLMyePRuXLl2Se+FZs2bhZDPDvQsKCnDu3DkUFhYiIiICERERcs/JJAJD\nQ5RaWuLjXr0Q0LMnPu7VC2WWlhAw7PaRVWmsaIGVOmXYtWsXmbFLIGgpMo1/ZGSkVLZIdnY2Pvro\nI3h4eOCbb76Re2J3d3eYmprK3J6fn4/JkyfDzMwM06ZNQ2lpqZKiq5f4+Hj0GjAAK9LT4bF6NVak\np6PXgAGNUh7VjawCq9raWlZkkMzzv3PnTpuq2FUWbfHtMgHRhRht0YVMn/+ff/6J//u//6OXjYyM\nMGnSJADArl27Wn3hgoICBAUF0cvdu3fHzZs3pWbHiggODoalpSUAwMTEBPb29nQJtUjRql5+/Pgx\nwsPD8dtvv+HGjRtwcnLCggULsGDBAqnGTeq6vmi5srJSaoCKyPiKWmyr+/o8Ho+u9wCAGzduAHjd\n5sHS0hIPHz5kVB9kWTOXRWiKPGwuFxUVsXJ9Ho+HvXv3AgBtL5tDps/fysoKJSUlMHzDzfH8+XMM\nGjQIFU10mnyT8vJy+Pj4NOnzDwwMRFBQELy8vAAArq6uSE5OhrW1tbSALPn8J0yYgJUrVzZav2bN\nGhxhMM9/69atOH/+PKKiouh1a9euhbu7O5YsWcKIDAEBAejUqRMWLFhAr9u5cyfq6upUVkVMIBBU\nS4u7etra2uLgwYOYM2eO1PpDhw5h4MCBrRbMxcUF169fp43/w4cPGxl+NunSpUuT601MTBiVo6Ki\nAl5eXoiLi4O+vj4EAgHGjh3LqJusZ8+eGDBggJQMbm5uKCsrY0wGAoGgWmQa/82bN8PX1xdZWVkY\nO3YsKIrC8ePHUVRUhIyMjFZf2MXFBUuXLsWMGTOQlZWlkgeKKgkODkZMTAxCQ0Npt8v27dsRHBzM\nqBx6enpwdHSEo6Oj1HomDe+dO3cQEBAAR0dHKReUpsVpmEbS3aXrEF2I0RZdyDT+gwYNwo0bN3D2\n7FkcPnwYHA4HCxcuxKhRoxQqIpo2bRrOnj2Lqqoq9O7dG19//TXdFC0kJATDhg3DiBEjMHToUJiZ\nmeHAgQOq+1QqQJRGmZSUhDt37qCsrAzBwcGMT/Jic5awqHCssLAQsbGxUm4f0UhLAoGgnSiU519b\nWwsOh0MXejEJWz5/TYGtWcJvztgtLi5W6QxfAoGgXuTZTpnGn6IobNu2Dd999x0ePXoEiqLQrVs3\nfP7551i8eHGbH+PIBhRFITs7G2PGjJHSr6qHp8tDsmI3NjaWVOwSCFpIixu77dmzh3J1daXS0tKo\n2tpa6unTp9Thw4epDz74gNqzZ0/Luw0pSTMiMgaXy1X7Nerr66nPPvuMGjx4MFVTU6P268kiNzeX\n6tmzJ7V582ZKKBQ22s6ELrQFogsxRBdiNEUX8mynTJ9/XFwcvv32WzobBwAmTZqELl26IDIyErNn\nz1bhM0q3kWzFnJubq9JWzMpAWjETCLqDTLdP//79UVpa2ii4KxAIMHDgQPz+++/MCNjG3T5l/xvM\nMmHCBERHR6u1I6csRIHdlJQUZGRk6ERHTgKhrdPiPP/OnTs3aYj09fXRWaLdMaF1LFmyBMuXL282\nhVSRrp4tRTKwm5+f3+ZGLRIIhKaR2dvn6tWrMDIyavJPVpfOtsqbJeyq5NixY3INf0pKCgIDAxEQ\nEIDAwECkpKSAy+W2+toVFRUYPny4UjN21akLbYPoQgzRhRht0YVM4y8QCFBbW9vkX0NDA5Mytmnk\nuXnU1dUzLy8Pbm5uCA4OJjN2CQQdRKbbhyCGzWo9WV09W5Nq25rArjZULjIF0YUYogsx2qILmW/+\nBNVSVlaG8PBwpYPXqqzwFQqFiIqKwqpVq3S6FTOBQCDGXyFa68PLysrCyJEjMWTIEKXf2H19fRvN\nEIiPj4ePj49S55GcsVtQUNDijB5t8WcyAdGFGKILMdqiC+L2aYbm5tYqAkVRiImJwbp165Camgp3\nd3elZZDsMSRK3VK2tUNFRQV8fX3h7OxMZuwSCAQAapzhqyrYyvNv7excPp+PRYsWITc3FxkZGbCy\nslKnuDLJy8sjM3YJBB2kxTN8dZ3WZtkIBAKYmZkhJyeHNcOfmJhIZuwSCIQmIcZfBrLm1ipqQDt2\n7Ih169ax0qpBnYFdbfFnMgHRhRiiCzHaogvi85cBm330W4OoYrempgYFBQWkYpdAIDQJ8fnLQJk+\n+hRFob6+Hh06dGBaTCkkA7ukFTOBoNu0uJ+/psBmYzdF+ujz+XyEhoaiffv22L59OytyAiSwSyAQ\npGlxP39NQRNElNWfu6qqiho9ejTl7e1NPXnyhFmhJNi/fz/VvXt36tixY2q/lqb0KtcEiC7EEF2I\n0RRdyLOdJODbQsrKyuDq6gonJyekp6e3ucAugUBo26jV7fPrr78iJCQEDQ0NWLx4MUJDQ6W283g8\n+Pn5wdraGsDrYTErVqyQFlAD+/kXFxfjww8/xPr16zFr1ixWZJAM7KalpZHALoFAkIJVn7+DgwO2\nbduGvn37wsvLC+fPn5cyUjweD5s3b0ZGRoZsATXQ+PP5fBQXF8PR0ZGV65PALoFAkAdrRV5PnjwB\nAIwcORJ9+/bFmDFjkJ+f32g/TTPsTfFm3m67du1YM/xst2LWlhxmJiC6EEN0IUZbdKG2PP+LFy9i\nwIAB9PKgQYNw4cIFeHt70+s4HA5yc3Nhb28PT09PLFy4EDY2No3OFRwcDEtLSwCAiYkJ7O3t6bap\nIkWrc7moqIjR68laTkxMRGhoKCIjIxEWFsaKPEVFRax9frKsucsiNEUeNpfZshc8Hg979+4FANpe\nNofa3D6nT5/Gnj17cPDgQQDA999/j7t372LNmjX0PrW1tdDX10e7du2wb98+HDlyBEePHpUWkGW3\nT1lZGbp3746uXbuyJoNoxu6hQ4eQmZlJZuwSCAS5sOb2cXZ2RllZGb1cUlICV1dXqX2MjIxgaGiI\ndu3aYc6cObh48SLq6urUJZLSiFoxN+WuYgpRK+acnJxWtWImEAgESdRm/Lt06QLgdcZPeXk5Tp06\nBRcXF6l97t+/Tz+ZMjMzYWdnx3qVLPA6DrF9+3YEBwcjLS0NhoaGrMghOWP39OnTGpHR8+bPfF2G\n6EIM0YUYbdGFWnv7bN26FSEhIeDz+Vi8eDG6deuGXbt2AQBCQkJw+PBhxMXFwcDAAHZ2dti0aZM6\nxVEIyVbMubm5sLKyYuU/k1TsEggEdULaO7zB999/j6NHjyI5OZmVwi2gdTN2CQQCASC9fZRGIBAA\nAPT19Rm7pggS2CUQCKqCDHNREn19/UaGnwm3j7YEdrXFn8kERBdiiC7EaIsuiPHXADQxsEsgENo2\nOuv24fP5WLVqFRYvXgxzc3OVn19RSGCXQCCoA3m2UycneVVXV2PKlCkwNDRkLY0TIIFdAoHAHjrn\n9iktLVW6FbOqfXja3IpZW/yZTEB0IYboQoy26EKn3vyzsrIQFBSkMa2YyYxdAoHAFjrl89+8eTOc\nnZ3h7u6ukvMpC2nFTCAQmILk+WsIJLBLIBCYhOT5q4DW+vASExPh5+eHhIQEhIWFabXh1xZ/JhMQ\nXYghuhCjLbposz7/hoYGGBiw+/EkK3a5XK7GFm4RCATdo026fbKyshAREYHCwkLWuoSSGbsEAoFN\ndMrtI9mKOS4ujjXDTyp2CQSCptNmjD+fz8e8efMQHx+P3NxcjBgxQmXnVsaHx/aMXXWjLf5MJiC6\nEEN0IUZbdNEmfP4CgQDjxo1Dp06dkJubCyMjI1bkIBW7BAJBW2gzPv+zZ89ixIgRpBUzgUAggOT5\nqx0S2CUQCJqITgV81YUsH54uBna1xZ/JBEQXYoguxGiLLrTO+PP5fJSUlDB6zaKiokbr2npgVxZN\n6UJXIboQQ3QhRlt0oVbj/+uvv2LgwIF49913ERMT0+Q+UVFRsLa2hpOTE8rKypo9X3V1Nby8vLB+\n/Xp1iCuTx48fSy23pYpdZXlTF7oM0YUYogsx2qILtRr/f//739i1axdOnz6NnTt3oqqqSmp7QUEB\nzp07h8LCQkRERCAiIkLmucrKyuhWzD/++KM6xZaJNrdiJhAIBEnUZvyfPHkCABg5ciT69u2LMWPG\nID8/X2qf/Px8TJ48GWZmZpg2bRpKS0ubPFd2djZGjhyJZcuWYcOGDYxn9JSXl2vNjN3/b+/+Y2r6\n/wCOP4tRyRDC/P4x8yPVne69FK6F3Ya0Yf2akvzI/Ghshj9s4Q+b2Yr5vTAj+dnWQsuv5HeW0BST\nO+4wFbFZofTjfv/w6XwUN77IuX3O67GdrXPO6757ddZe99xz3+d1WprValU7BYchx+Jfciz+1WqO\nha2FXLhwwRYeHq6s796927Zu3bpGMXPmzLGdO3dOWTcajTaLxdIoBpBFFllkkeUXluaoepOXzWb7\nZipS0+vnTfcLIYT4fS122Uev1zf6AreoqIgxY8Y0ijEajTx8+FBZf/PmDYMGDWqplIQQQvyjxYp/\np06dgC8zfqxWKxcuXMBoNDaKMRqNpKWl8fbtW1JTUxk+fHhLpSOEEOIrLXrZZ+vWrcTFxVFTU0N8\nfDzdunVj7969AMTFxWEwGBg3bhx+fn54eHiQkpLSkukIIYT4h8O3d1Db1atXiYuLo7a2lvj4eJYv\nX652SqqJjY3l7NmzeHp68uDBA7XTUdWLFy+Ijo7m9evXdO/enUWLFhEZGal2Wn9dVVUVJpOJ6upq\nXFxcCAsLY+XKlWqnpaq6ujr8/Pzo06cPp0+fVjsdu6T4/4BOp2Pbtm30798fs9nM9evXNdHG4Xuu\nXbuGu7s70dHRmi/+paWllJaW4uvrS3l5OQaDgYKCAtU6yqrp48ePuLm5UV1dzejRo0lPT2fIkCFq\np6WaxMRE8vPzqaioICMjQ+107Gp17R3+pp+5V0FLxo8fT5cuXdROwyH07NkTX19fALp168bIkSO5\nc+eOylmpw83NDfjS5LC2tla1hyg5gpcvX5KZmcmCBQscfqaiFP9m5OXlMWzYMGV9xIgR5ObmqpiR\ncEQWi4WioiIMBoPaqaiivr4eHx8fevTowbJly+jbt6/aKalm5cqVbNmyBWdnxy+tjp+hEA6soqKC\nsLAwkpKS6NChg9rpqMLZ2ZmCggIsFgu7du3i3r17aqekijNnzuDp6YlOp3P4s36Q4t+sn7lXQWhX\nTU0Ns2bNIioqipCQELXTUd2AAQOYOnWqZi+N3rx5k4yMDAYOHEhERATZ2dlER0ernZZdUvyb8TP3\nKghtstlszJ8/Hy8vL1asWKF2OqopLy9Xuli+ffuW8+fPa/aNcNOmTbx48YJnz55x7NgxAgMDOXTo\nkNpp2SXF/wca7lWYPHkyS5Ys0exMH4CIiAj8/f0pLi6mb9++qnVXdQQ3btwgJSWF7OxsdDodOp2O\nrKwstdP660pKSggMDMTHx4fIyEhWrVpFr1691E7LITh6q3eZ6imEEBokZ/5CCKFBUvyFEEKDpPgL\nIYQGSfEXQggNkuIvWo3k5GRMJhPe3t7odDry8vIAWLhwod1HgLaU06dPs3nz5t8e5/3790yePLnZ\nGKvVyqhRo767LyEhgezsbABCQkIoKSn57ZyENshsH9EqvHr1iqCgIHJzc3Fzc+Pdu3dUV1e3+mmF\nSUlJdOzYkQULFtiNsVqtBAcH/7CZ3smTJyksLGTDhg1/Ok3xHyRn/qJVKC4uxtPTU2ki5uHhoRT+\niRMnkp+fD3zpxxQdHY3RaGTt2rVUV1eTlZVFaGioMlZOTg7BwcF245tKTU1l7Nixylx2gIMHDyrt\nvX19fZW5/m5ubly7do2qqioSExMxmUxMmzaNnJyc7/5d+/fvJzw8XBkzNDSUwMBAfH19OXLkiBJX\nX1/P0qVLGTFiBIsXL6ampgaAmJgY0tLSAJgxYwapqam/doCF5kjxF62CyWSivr6e/v37Ex8fj8Vi\nUfY5OTkpN9SsXr2a7du3c/v2bWw2G+np6UyZMoXbt2/z6dMnAI4fP05ERITd+KY2btzIpUuXKCgo\nYM+ePcrvbHD//n3u3bvHxo0b0ev1jB07lmPHjtG2bVuuXLnCgQMHWLNmzTfjfvr0icrKStzd3ZVt\nFy9eZP/+/WRlZbF+/XrKy8sBePToETNnzqSwsBCr1cqtW7e+yaN9+/a4urpSWlr6awdZaIoUf9Eq\nODk5kZ2dzalTp3B1dSUgIIDMzMxGMfn5+RQWFjJx4kR0Oh1nzpzh6tWrtGnThqCgIDIyMqitrSUz\nM5OQkBC78U35+fkRERHBqVOn7DZve/LkCatXr+bEiRO0bduWtLQ0kpOT0el0BAUFUVZWxrNnzxq9\nxmKx0K9fv0bbxo8fz8CBA+nZsydms5lz584B0Lt3byZNmoSzszMmk0kp/k0NHjyYx48f//RxFdrV\noo9xFOJP0+v16PV6hg8fztGjR5k6daqyr76+Hi8vLy5fvvzN68LDw9mxYwceHh74+fnRoUMH6urq\n7MZ/LSUlhZs3b5KSksKWLVuUTwkNKisrCQsLY9++ffTo0UPJZefOnUyYMKHZsZt+5dZ0veFTTefO\nnZVt7dq1o7Ky0u54raGdsFCf/JeIVqG4uJgnT54AUFtbS25uLv7+/o1i9Ho9ZWVlyjMXPnz4oLzG\nZDJx9+5dkpOTlUs+BoPBbnwDm82G1WrF39+fxMRESkpKqKqqahQTGxvLvHnzCAgIULZFRkayd+9e\nKioqAL7b5njIkCE8f/680bbr169jtVopKyvj/PnzmM3m/6s98NOnTxk6dOhPxwvtkuIvWoXKykpi\nYmIYOXIkAQEBuLi4MHfu3G/iDh8+zO7du/H29sbf31+5BOLs7Mz06dPJyspi+vTpP4xvUFdXR1RU\nFN7e3kyaNIn169fj4uKinJE/f/6ctLQ0Dhw4oHzpe/fuXWbPno3BYMBsNuPl5UVCQsI3ubq6uuLu\n7q68QTg5OTFlyhRiY2Mxm80kJCTQtWtXZd/Xvl5v+Pnz5898/PhR+fQhRHNkqqcQKkpKSsLd3Z2F\nCxdy8OBB8vPz2b59+y+NdeLECYqKimSqp/gpcuYvhIrmz5/P8ePHgcazln7FkSNHWLx48Z9KTfzH\nyZm/EEJokJz5CyGEBknxF0IIDZLiL4QQGiTFXwghNEiKvxBCaJAUfyGE0KD/AUmTbdQ4J8sJAAAA\nAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 14 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }