{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Demo of running AMRClaw Fortran code and plotting results" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This notebook can be found in `$CLAW/apps/notebooks/amrclaw_advection_2d_square`\n", "\n", "The is adapted from $CLAW/amrclaw/examples/advection_2d_square to illustrate the use of IPython notebooks.\n", "\n", "The setrun.py and setplot.py files are taken from that example but some parameters are modified in the notebook below.\n", "\n", "The 2-dimensional advection equation $q_t + uq_x + vq_y = 0$ is solved in the unit square with periodic boundary conditions. The constant advection velocities $u$ and $v$ are specified in setrun.py but can be changed below." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Have plots appear inline in notebook:" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "source": [ "%pylab inline" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Check that the CLAW environment variable is set. (It must be set in the Unix shell before starting the notebook server)." ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Using Clawpack from /Users/rjl/git/clawpack\n" ] } ], "source": [ "import os\n", "try:\n", " CLAW = os.environ['CLAW'] \n", " print \"Using Clawpack from \", CLAW\n", "except:\n", " print \"*** Environment variable CLAW must be set to run code\"" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Module with functions used to execute system commands and capture output:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [], "source": [ "from clawpack.clawutil import nbtools" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Compile the code:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Executing shell command: make new\n", "Done... Check this file to see output:\n" ] }, { "data": { "text/html": [ "compile_output.txt
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/compile_output.txt" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "nbtools.make_exe(new=True) # new=True ==> force recompilation of all code" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Make documentation files:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "See the README.html file for links to input files...\n" ] }, { "data": { "text/html": [ "README.html
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/README.html" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "nbtools.make_htmls()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Run the code and plot results using the setrun.py and setplot.py files in this directory:\n", "\n", "First create data files needed for the Fortran code, using parameters specified in setrun.py:" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [], "source": [ "nbtools.make_data(verbose=False)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now run the code and produce plots. Specifying a label insures the resulting plot directory will persist after later runs are done below." ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Executing shell command: make output OUTDIR=_output_1\n", "Done... Check this file to see output:\n" ] }, { "data": { "text/html": [ "run_output_1.txt
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/run_output_1.txt" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Executing shell command: make plots OUTDIR=_output_1 PLOTDIR=_plots_1\n", "Done... Check this file to see output:\n" ] }, { "data": { "text/html": [ "plot_output_1.txt
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/plot_output_1.txt" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "View plots created at this link:\n" ] }, { "data": { "text/html": [ "_plots_1/_PlotIndex.html
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/_plots_1/_PlotIndex.html" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "outdir,plotdir = nbtools.make_output_and_plots(label='1')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Display the animation inline:\n", "\n", "Clicking on the _PlotIndex link above, you can view an animation of the results. \n", "\n", "After creating all png files in the _plots directory, these can also be combined in an animation that is displayed inline:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "\n", "\n", "\n", "
\n", " \n", "
\n", " \n", "
\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
\n", "\n", " \n", " \n", " Once \n", " Loop \n", " Reflect \n", "
\n", "\n", "
\n", "\n", "\n", "\n" ], "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "import clawpack.visclaw.JSAnimation.JSAnimation_frametools as J\n", "anim = J.make_anim(plotdir, figno=0, figsize=(6,4))\n", "anim" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Illustrate how to adjust some parameters and rerun the code: \n", "\n", "First read in the current rundata. (See the README.html file for a link to `setrun.py`)." ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import setrun\n", "rundata = setrun.setrun()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we adjust the mesh to be $40 \\times 40$ and set two gauges at $(0.6, 0.4)$ and $(0.6,0.8)$." ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false }, "outputs": [], "source": [ "rundata.clawdata.num_cells = [40,40]\n", "rundata.gaugedata.gauges = [[1, 0.6, 0.4, 0., 10.], [2, 0.6, 0.8, 0., 10.]]" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The advection velocities can also be changed (in setrun.py, $(u,v) = (0.5,1)$." ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false }, "outputs": [], "source": [ "rundata.probdata.u = -1.0\n", "rundata.probdata.v = 1.0" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Change the limiter if you want to experiment with other methods. A list of 1 limiter is required since this is a scalar problem with only 1 wave in the Riemann solution.\n", "\n", "Here we set the limiter to ['none'] to see the oscillations that arise if Lax-Wendroff is used with no limiter." ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false }, "outputs": [], "source": [ "rundata.clawdata.limiter = ['none']" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Write the data out (for the Fortran code to read in) and run the code:" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Executing shell command: make output OUTDIR=_output_2\n", "Done... Check this file to see output:\n" ] }, { "data": { "text/html": [ "run_output_2.txt
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/run_output_2.txt" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "Executing shell command: make plots OUTDIR=_output_2 PLOTDIR=_plots_2\n", "Done... Check this file to see output:\n" ] }, { "data": { "text/html": [ "plot_output_2.txt
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/plot_output_2.txt" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "View plots created at this link:\n" ] }, { "data": { "text/html": [ "_plots_2/_PlotIndex.html
" ], "text/plain": [ "/Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/_plots_2/_PlotIndex.html" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "rundata.write()\n", "outdir, plotdir = nbtools.make_output_and_plots(label='2')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Changing plot parameters and plotting inline\n", "\n", "In addition to producing a _plots directory using the parameters in setplot.py, plot parameters can be changed and plots produced inline as the examples below illustrate..." ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [], "source": [ "import setplot\n", "plotdata = setplot.setplot()\n", "plotdata.outdir = outdir" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plot one frame of the solution:" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " Reading Frame 2 at t = 0.4 from outdir = /Users/rjl/git/clawpack/apps/notebooks/amrclaw/advection_2d_square/_output_2\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUUAAAEKCAYAAACFekfkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8FHX+x/HXJiEkEFoISO9ditxRREVBRIoCiniiUhRF\nLCj8bCh4HpyKwqGeihU9KwIiIKKCghoFpYkK0kuAQGiGFkIxZJPfH5vdzCZbJpvs7iR5Px+PeWQ3\n852Zz87sfPY7M9+Zry372NJsREQEgIhwByAiYiVKiiIiBkqKIiIGSooiIgZKiiIiBkqKIiIGxT4p\n7kk+RETVq8nKygpo+plzv6XXDY8VcVTBUZxiFSmuLJMUV6zayCW9xlC5wXVUbTyQy/qM5ZffthXp\nMjwl0Ftv7MHX854r0uWY0a3fQ7zz4WKv48Mda0TVq0naczAkyyqIPcmH6N7/YcrX6UfLziP49odf\n/U6TkXGelp1HULf1LQWa18effkf9trcSV7cf1w+dyPETp1zj/vorgxGjp1Gp/gBqtryJF1+b5zbt\n73/s5O/d76V8nX50uPJe1m/c5Tb+xdfmUbPlTVSqP4A77n+ejIzzrnHHjqdx/dCJxNXtR4N2Q5g1\n7zu3ab/94VdadB5B+Tr9uHLAIyTvP+J3HYh5lkiKaWmnuXbwE4wZdT3Hdy8gZdNs/vXoUMqWjQ7K\n8rKzw99e3WazmSoXzlitsJ7yuvnOyfy9XVOO7ZrHM0/czqDbniL16Emf0/znlU+oXq0yeVe5r3lt\n2rKHux98iZlvPc7hrZ9QLrYs9z78smvaiVM+ZNeeAyT/8THfL/wPU1/5hK+/XQs4kvCAW//FsJuu\n4sTuBQwffDUDbv0X589nAvD1t2uZ8vIcvvtsKns3zCRp70H+9dwHrnnf98grxJQtw5Ftc5n55mPc\n89DLbN66F4DUoye5Yfi/eWbC7RxPmk+Hi5px04inC71eJZclkuL2Xfux2WzcNLAbNpuNmJhoenb/\nO21aNQQcO+fT02bSoN0QLmh+I8PvnUpa2mmP82rQbojbL/7E5z5g6N2O2tXl1zwIQOWG11Ox3gBW\nrd3Mex9/Tde+/+cq//PqTXTscR+VG1xHp6tGs3LNZte4bv0e4snJ73FZn7FUrDeAXjc8xtFjaR7j\nOHEynWsHP0H1ZjcS32gg/W7+JykHUgGY8PT/WL7yD0aPm06Fev154LFX801vJtaIqlfz+v8W0bTD\ncCrWG8CTk99j1+4DdLn6ASo3uI7Bdzzj2hEBvvh6FRddPooqDa/n0t5j+GPzbo+xO5fd7vJRVKjX\nn7mf/eCxXKht37mf3/7YyaTHhlG2bDQD+3Wl7YUNmbdouddpdu89yMy53/H42Jsx5nh/85r56bf0\n79OFyy5uTfnysTw1/jbmf/ETp0+fBeCDOUv558NDqFSxPC2a1eOuYX15b9Y3ACSuWI/dnsWYuwdS\npkwU9991HdnZ2Xy//HcA3p+9lDuH9qFl8/pUrhTHk48McU17+vRZ5n/xE0+Nv41y5WK49OLWDOh7\nCR9+sgyA+YtW0LplA27o35Xo6DJMHDeU9ZuS2L5zf5Gv79LKEkmxeZO6REZGcNt9U1mybK3bYQrA\nuzO/5v3Z35C4aBpJv35AevpZRo+b7nFeNpvNrRZmfL38qxcBOLnnM9KSF3Jxx1Zu0x47nsY1g59g\n7N0DOZY0nwfvuYFrBj/hFs+s+d/z3quPcGT7XDLOZzJt+lyPcWRlZXHHkN4kb5hJ8oaZxMZEu2J+\n5okRdO3Shlen3s+p5M95+bn78k3vL1anb75fx2+Jb7Dqm5eZ8vInjBz7IrPeHk/yhpn8sXk3s+Z9\nD8BvG3ZyxwPPM+O/D3IsaT6jbruW/rc86XbY5vTjly8AsGH5W5xK/pwbr7siX5kVqzZSpeH1Xoef\nV2/yGG9hbNq6h0b1a1K+fKzrf+1aN2ZTTi3Kk/vHvcqzT95BTIz7UYe/eW3aupd2FzZyjWvUoCZl\no8uwfVcKx0+c4uChY7RrnTu+7YWN3KZte2FDt+W1a507fvM293m3vbARh48c5/iJU2zflUJUVCRN\nGtXOnfbCRmzauscVt3G55crF0KRhLTZu8fwDJwVniaRYoUI5Vnz1IjZsjBz7ItWb3ciAW5/kyJ/H\nAcev9kP3DaJBvRqULx/Ls0/ewez5iaYurhgPAf0dDn75zWqaN6nDrTf2ICIigsE3dKdF07p8vngl\n4Eiwt9/SiyaNahMTE80/rruC3//Y5XFe8VUqcv21lxETE01cXCzjH7yZH37a4DU2X3H78uj9/yAu\nLpZWLerTplUD+vToSIN6NahYsTx9rurIbxt2AvDW+18yavg1dPxbc2w2G8MG96Rs2TKs+mWLqeXk\nddnFrTm+e4HX4ZLOFwY0X1/ST5+lUsXybv+rWKEcp06d8Vh+wRcryM7OZkDfS0zNq0JcLOnpjprg\n6TPnPC8r/YyrjHG8Y9xZH3GW51T6mZzx5/JNC3Aq/Szpp8+63ufGZZz3OSpWyB9X+ulzHteBFFxU\nuANwatGsHu+++ggA23bsY8jdzzF2/Ot8PGM8Bw8do36dC1xl69WpTmamncNHjhdpDAcOHaVenepu\n/6tftzoHDh11va9RPd71OjYmmvScw6m8zpw5x/9NeIOvv/vFVdNMP32O7OxsV+3V7HlFXy6oXsUQ\nT1m39zEx0Rz58wQAe/cd5oM5y3hlxkLX+POZmRw8dKzQMYRKXPlY0vIkwBMn0/MlEXAchj46cQaL\nP5lsel4n005TIS7WNf5knlM0J0+dpkJcOeJyyqSdOkNC1Uoep/UUZ4W4cjnjY9zGO5dTIS7Wa1zO\nz1ghLpa0U6fzjXcuWwrPEjXFvJo3rcvwwT3ZuGUPALVqVmXPvsOu8cn7jxAVFemWAJzKl4vh9Jnc\nX81DR3J3en9JqHbNBPbuc7+St3ffEWrXTCjwZ3j+1U/Zvms/a5a9wsm9C/nhi+fJzs521QD9xVIU\nCdNG7jzq1anOhAdvdqvNpe9bxE0DuwU07+Ur/6BCvf5eh59WbSx0/Hld2KIBSXsPumpqAOs3JnFh\niwb5yu5ISmHvviN0veZBara8iRuGT+Lg4aPUbHkTyfuP+J3XhS3qs35Tkmvcrt0HyMjIpFnj2lSp\nXIGaNeLdjhLWb0yidcvcaTdscj+c/WPzbi5sUd/1OX7faJx2FxdUr0KVyhVo1rg2mZl2dial5I7f\ntMtt2vUbc+M6ffosu/Yc9LgOJDCWSIrbduzjhVc/dV2I2Lf/CLPmfU+XnPNoNw/szouvz2NP8iHS\n088y/qn/MXhgNyIi8od/UZvGzJ6fSGamnV9+28a8RStcCaZa1UpERNjYtfuAxzj6XNWJ7bv2M2ve\nd2Rm2pkzP5GtO/Zxba/OrjJmD2vTT58lNiaaShXLc+x4GpOmfuQ2/oJqldm1x3McZmL1xu10AblJ\neOSwvrzx7hesWbeV7OxsTp8+y5ffrHZLCm7xVa/ic9ldu7ThVPLnXodLL25doLjNaNakDhe1bsyk\nqR9y7lwG8xctZ+OWPdzQv2u+sm1aNWT/xo9Z/+MbrP/xDd5+6UEuqF6F9T++QZ1aCX7ndeugHixa\nspIVqzZy+vRZ/jn5PW7od5nrHOSwm3ry9PMzOXEynS3b9vL2R4u57earAeh2WTsiIyN4+c0F/PVX\nBi+/uYCIiAiuvLx9zrRX8c5HS9iybS/HT5ziqWkzuf0Wx7Tly8cy8NpLefLZ9zlz5hwrVm1k0ZJV\nDP3HVQBcf+2lbNyyh/mLlnPuXAaTpn7ERa0b06xJnSJf36WVJZJihbhYVq/bSuee9xNXtx9deo2h\n7YWNeP6pUQCMGNKbof+4isuveZBGfxtGuXJleWXKaNf0xlrVU+OHs2v3Aao0up6JUz7k1kFXusaV\nKxfDhAdv4dI+/0d8o4Gs/mWL24WZqvEV+WLWUzz/6qckNLmBaa/O5YtZTxFfpaLHZeW9qGM09u6B\nnD2XQULTQVzSeyx9enR0Kztm1PV8+vly4hsNZOzjr+Wb3l+seWPxGB+55f9+UTNm/Pf/GD1uOvGN\nBtK04218MHupx9gBJo4byvD7plKl4fV8uvBHr+VCbfY7E/jl9+3ENx7IhKffZd77T1I13rF9nLVX\ngMjISKpXq+IaqlSuQGREBNWrVXH9mPqaV6sW9XnjhTHcetezXNDiH5w9l8Fr0x5wxTHpsWE0blCL\n+m1vpfuARxj3wE1cfWUHAKKjy/DZRxP5YM5SqjQayAdzlvLZRxOJiooEoFePjjx6/z/oPuARGrQb\nQuMGtZj02HDXvF+b9gBnz2VQvfmNDBn1HG88P4aWzR01xYSqlZj3/pNMePpd4hsP5JfftzP7nQlB\nXeeljc3fQ2ZHjJ7Gl0tXUz2hMn/8NMNjmQcee5XFy9ZQLjaG9159hPZtmwQlWBERp2DlJr81xdtv\n7cWSuc96Hf/V0tXsTEphxy/v89aLY7nnoZf8LlREpLCClZv8JsWuXdpQpXKc1/GfL17J8MGO8yGd\nO7TkRNrpIr8qLCKSV7ByU6HPKaYcPErd2tVc7+vUSmD/gT8LO1sRkUIJNDcVyYWWvFdki6I5iYhI\nYQWSmwrdeLt2zarsS8nNvvsPpHps12eL71nYRYlIgLKPeW9pEE4FzQsF+Rxmc1NehU6K/ft0YfqM\nhQy+oTur1m6mcsXyHhtVAzyNtZoOROC4TTCaDOJIpx7JAHRiDVW7HmViFZg4ABgM5L9pwuEM4Lzb\nMIaw3iM0caJjKE4KHXNGzt90IBlYk/P+Mzi8+AJmMxiAVVxMOo7zTwmkUo9kKuO42yeCLDKIzpld\nNFlEuL4bzvEAWTkHVt/yIz24vBBBh9YTPBPuEHyKwG6qXBaRBZpvQXKTkd9d+OY7n+GHnzeQejSN\nuq1vYdJjwzif6XjyyqjbrqVvz858tXQNTf4+nPLlYnh3+sMFClxESjfjD5AveUsFKzf5TYqz3vZf\nu5s+9X5TC5vKo6bKhYqxphjPMZqxHYA9NOCeFa/D5ThqH5ne50EWuVsrsId/S6CM6/4ccAJw3gG3\nAmYzmGk4doT9uN/xEUUm1XHc0hlBFueIyZlNjFtN0dMOew476yz2XfbN6jXFwHacosxNRiE92Euj\nov9CYXKE6hygFgDxHOOq7GV0q7IDUvGd7LL8jA+hbt3CHUHBFTpm57rPxHEqw3GnKIdPXcAqLs6X\nDHEVj3Jt74LrRYaFv8vFTZTPWkeuDP9FioRlnpJjRd0u8F/GSkplUgyLbuEOoEQJtKYYLEqKBs6a\nbDpxrsMpUz9i1tqmpVcWru11jhjXhRWxNiVFC6vBIcBxdTKOdDDbRYwlHqtRSkUY/kbh2mYVSSPB\neSwtlqakaFGVOUGjnLP0ddjv2KGczXCU9KzPmRRzKodVIo9Tz57sOl+V6eGr7hxnbJIjoaekKCJi\nUKqTYhzpoVycR54a5UaTQWVOuA63KnPC/fDZV01RtUhrcNYUc04FEweVT55wNbvJe6U50CY5WcVw\ng4d/r/OtVCdFs5feQ8X4xXcOzvc2W7YSXjFn3Ka+xuV97emvcTopWtbLC+Kdvv/FnqsVQQlyxus9\np8WTsVLiawiVkNYUrbYxnSs6kyhiOOfagTKIxp4dSSR2SzXOFh+Mjbhz/jovnnhKjBlEu/4fQZbr\nu2nmjhYrsNq+VBhWW8chTYpWvcKXRQTniHFLiplEEZll7kZ1i23T0sv4A5bpnvg8MSZF52vnVeri\neO6wuCrVSdGqjJV043uLbSspiCzctqnnIhEeX0toKSmKiBgoKVqQp6vPOS+kuIrwffXZUST/1WfV\nGENPSdFCfDXJsdqGEpMi8r4tWFKU0LNak5xSnRRFJPys9oOkYwUpsay2s4lnhWmnuGTZWlp0HkHT\nDsOZ8tLsfOOPnzjF9UMn0q7rKDpfdT+btuwxEY9ICaXzg8VDoEnRbrczetx0lsydzOaV7zBr3vds\n2bbXrczkF2bxt7ZNWL/8TT54/VHGjH/NRDySj3amYqwAje39NdmR0Ag0Ka5Zt40mDWvRoF4NypSJ\nYvDA7ixcvNKtzJbtyXTvehEAzZvWZU/yYf5MPeEnHhGRMAo0KaYcTM3X2X3KQfdnaLZr3Yj5i1YA\nsGbdVvbuO8z+A76fs6kLLSISVt7OF6axjjR+9TqdmY7tHxszmDGPv0b7K+6mTauGtG/bhMhI33VB\nJUURCStvTXLiaUc87VzvD/C22/i8nd3vS/mTOrXcO7uvUKEc/zN0bdrwoqE0ql/TZzw6fBaRsAr0\n8LlD++bsSEphT/IhMjLOM2dBIv17d3ErczLtNBkZ5wGY8f5XXHFpG+LiYn3Go5qiB2rKIRI6ge5v\nUVGRTJ8yml6DHsduz+KOIb1p2bw+b773BQCjbruWzVv3ctvo/2Cz2WjdogHvvPKQ//kGFI2ISBEp\nTCWkT89O9OnZye1/o2671vW6S6dWbFvzboHmqaQoImFltSMzJUURCSslRRERAz0QQkTEQDVFERED\nJUUREQMlRRERAyVFEREDJUUREQMlRRERAzXJERExsFpN0e9Tcvz1gZB69CS9Bz3ORZePovUlI3nv\n46+DEqiIlEyF6aMlOPH4YKYPhOkzFtK+bRN+//FNEhdN46F/vklmpj2oQYtIyVGskqKZPhBq1ogn\n7dRpANJOnaFqfEWioiKDF7GIlChWS4o+zyl66gNh9bqtbmVGDuvLlQMeoVarmziVfpZP3nkiOJGK\nSIlktXOKPpOimT4QJr8wi4taNyZx0fPs2n2AngPHsf7HN6lQoZyH0hMNr7vlDCJStBJzhuKhWCVF\nM30g/LxmExMevAWAxg1r0bB+Dbbt3EeH9s09zHFioQMWEX+64V7hmBSeMEwqTJOcJcvWMnbC69jt\ndu4c2odxYwa7jU89epIho57j0JFjZGZm8fDoQdx2Sy+f8/R5TtFMHwgtmtZj2Q+/AXD4yHG27dhP\nowa+O4YREXEK9JxisC4E+6wpmukDYfyDN3P76Gm06zqKrKwspk4aSXyVigVdLyJSSgV6+Gy8EAy4\nLgS3bF7fVaZmjXg2bEoCzF8I9tt4218fCAlVK7Fo1lPmP4mIiEGgSTFYF4J1R4uIhJW3pLiPnexn\np9fpiv5CsDMeEZEw8nYOsT6NuJSrXUNeZi8E3zjgcsD9QrDveEREwijQCy3BuhCsw2cRCatAm+QE\n60KwkqKIhFVhGm8H40KwkqKIhFWxuqNFRCTYlBRFRAyUFEVEDJQURUQMlBRFRAzUcZWIiIFqiiIi\nBkqKIiIGSooiIgZKiiIiBkqKIiIGuvosImKgmqKIiIGSooiIgZKiiIiBkqKIiEFhkuKSZWsZO+F1\n7HY7dw7tw7gxg93GT3tlLjM//RaAzEw7W7Ynk7pzHpUrxXmdp5KiiIRVoEnRbrczetx0li2YQu2a\nCXTscR/9e3dx6/f54ftv5OH7bwTgi69X8d835vtMiI54RETCKIpMU0Nea9Zto0nDWjSoV4MyZaIY\nPLA7Cxev9Lqcjz/9jpsHdvcbj5KiiIRVoL35pRxMpW7taq73dWolkHIw1eMyzpw5x9ff/cIN/bv6\njUeHzyISVt4OnzdxhE386XEcgM1mM72MRUtWcdnFrf0eOoOSooiEmbek2IYE2pDbuf2nbHYbX7tm\nVfal5CbNfSl/UqdWAp7MXpBo6tDZEY+ISBgFevjcoX1zdiSlsCf5EBkZ55mzIJH+vbvkK3cy7TQ/\n/ryBAX0vMRWPaooiElaBXn2Oiopk+pTR9Br0OHZ7FncM6U3L5vV5870vgNz+nz/78id6XdmB2Niy\n5uYbUDQiIkWkMO0U+/TsRJ+endz+50yGTsNvvprhN19tep5KiiISVnpKjoiIgW7zExExUFIUETFQ\nUhQRMbBaUvTbTnHJsrW06DyCph2GM+Wl2R7LJK5YT/sr7qb1JSPp1u+hIg9SREquQNspBovPmqKZ\np1CcOJnOfY+8wtefPkud2tVIPXoy6EGLSMlRrGqKZp5C8fGn33FDv8uok3NjdkLVSsGLVkRKnECf\nkhMsPpOimadQ7NiVwrETp+je/2E6XHkvH85ZGpxIRaREKlaHz2aeQnE+M5Nf1+/k28+mcubsX3Tp\n9QAXd2hJ08Z1iixIESm5rHb47DMpmnkKRd3a1UiIr0RsbFliY8tyeZc2rN+Y5CUpTjS87pYziEjR\nSswZiodilRSNT6GoVaMqcxYkMmvGeLcyA/pcwuhx07Hb7fz113lWr9vKg/cO8jLHiUUUtoh41w33\nCsek8IRhUrFKimaeQtGiWT16X9mRtpeNIiLCxshhfWnVor6v2YqIuBSrpAjmnkJh7BxGRKQgil1S\nFBEJJqs9JUdP3haRsCpMk5xg3HGnmqKIhFUw+30O5I471RRFJKwCrSkG6447JUURCatg9vscyB13\nOnwWkbAK9PA5WHfcKSmKSFh5u/r8M3+xkr+8Tlf0d9w56PBZRMLK2+HyZZThEeJcQ15m+n0e0OcS\nVqzeiN1u58yZc6xet5VWzX3fXKKaooiEVTD7fQ7kjjslRSlZCnDsY7U7KUqrYPf7XNA77pQURSSs\nrPbjpKQoUsxEkEVWCbocUKqTotU+vDMeT22hrBarFFAEXtu35RZx3/6+ylhNQeKy5ifIZbV1HNKk\nGE1GKBfnl3NjRJPhGpzvI8gyf36q5PxoF28R5G6LCNy2qSfOcRFkkZmzKzh/Ho0Js7g7F+4A/NAD\nIfKI55hlxlXmhPs/dHKh+MmzzQq0vU1OF+5xxvcFGWdVxaqPlqLmbSOF60tmrCk6/1eZE8SRTmRZ\nu6NwFBb46RC/jLXEKKAcxJ1OBxzb+wjV3Yo7E2JlThBBFjE59alzxLjO18VzzOvOWJwS4wGvc7UG\nq9XGQ5oUW7A1lIvzy5gUK3OCRiQBUI9kSAAP7UU9zcTzawkNw+GyMxkCEA/1/kymGdsBOEAt0qgI\nQA0O0YgkEkjNmTSLc8QAuUmxJB0+KykWTEiTYl++CuXiTIsmgzjSHckQ6MQaaAbE5xTwleyUCMPH\nuO6dCdF5l1cT6LRtDXtoADhqTOk5v3IJpFKH/a7aYgRZZBANQAbRPpNi3qu+VtuhPfku3AH4YbV1\nqF3aF/XSWuzdY3s93CEUudFMD3cIRcpq5xRt2ceWZodkQfE9OUfZUCzKNGNtILKsPbeW0Qy4GLgy\n5/3lkFORyM94cVPnH0PPeeHyDJAKOUfLsApHFWmF4+2O7KauQ+Q40kkglTgc5xxttmzs2ZE5s4ty\nqw1arRYTiBj+IvuY/0dmhYMtvifZ/h924yibTUg+R0gPn8tO8P7Ei7Byno9ynkOMx1FLbJTz3tda\niiK3IZgSYug513kMju3WxPC+Hq4ftqapO3ITaDSOQ+3o3HlEZjkurEVm2R3bsyRty2fCHYAfZrPQ\n+aBG4RLaRiePhnRpBZP3gkkUuWvH3znFkrQDFTfOdR+dM1TMeV8P6ERuIsxb4cu73Yp/hdA7qydF\ni+0/aoknIuFVqpNiRf9FLMt47tBZkxTryVtzFOsr1UmxJMhEa02kKCkpiogYWCwpWiwcC3OemFez\nG5GiFWFy8GDJsrW06DyCph2GM+Wl2fnGJ65YT6X6A2h/xd20v+Junp420284qimapTUlEhwB7lt2\nu53R46azbMEUatdMoGOP++jfuwst8/TBcsWlbfn846eCHU6A0kO6tIIz3kdrvJiimmHxkEVu05os\nHOd/je/z0na1hgC3w5p122jSsBYN6tUAYPDA7ixcvDJfUswu4O0p+lqISHgFePiccjCVurWrud7X\nqZVAysFUtzI2G/y8ZhPtuo6i7z/Gs3nrXr/hhLamaK1nSeby9NMQge5UKc6yyF9zNLM9S3IjbqsK\ncP+y2fzfH/i3tk3Z98fHlCsXw+Kla7hu6L/YvvY9n9OENimW818krPLe1aJkWLzkefK22y2YYl1e\n9rPEc5Do487g2jWrsi/lT9f7fSl/UqdWgluZChVyk06fnp2495FXOHY8jfgq3htN6/KBlGzRYLFe\nMAqvpH0mL0mxWznH4DQpzX18h/bN2ZGUwp7kQ9SqUZU5CxKZNWO8W5nDR45TvVplbDYba9ZtJTs7\n22dChFAnRd1hIKFirDWWxO9dSfpMAR6RRUVFMn3KaHoNehy7PYs7hvSmZfP6vPneF4Cj/+dPP/+R\n1//3BVFREZSLjWH22+P9zDXEjw4r6FUgESk8my00j9wKhC2+J9lNTZbdUQIfHSYiko/Fzt37Dcdf\ni3Gntb9uI6paL+YvWl6kAYpICVeIO1qCFY5XzhbjS+ZOZvPKd5g173u2bMvfzsdutzNu0tv07tFR\nh8giUjDFKSkaW4yXKRPlajGe1ytvLWRQ/65Uq1opaIGKSAlVnJKimRbjKQdSWbj4Z+4Z0Q9wnNQV\nETHNYknR54UWMy3Gx45/jeeevAObzUY22Tp8FpGCsdjlXp/hmGkxvm79DgbfORmA1GMnWbxsLWXK\nRNK/zyX55jdxYu7rbt0cg4gUrcREx1BsWOzqs892ipmZdpp3up1vP5tKrRpV6XTVaGbNGJ/vKRRO\nt9/3H/r1vpiB/brmX5DaKYqEheXbKXY2WXa1BdopmmkxLiJSKMWpplikC1JNUSQsLF9TvNRk2Z8s\nUFMUEQk6i9UUlRRFJLyUFEVEDCyWhSwWjoiUOqopiogYKCmKiBhYLClaLBwRKXUKce9zMB5tqKQo\nIuEVYFIM1qMNlRRFJLwCTIrBerShkqKIhFeUySGPYD3aUBdaRCS8vPX7fAQS//Q8DoL3aEMlRREJ\nL2/9PtdwDE6TNruPL+pHGzopKYpIeAV4Eq9D++bsSEphT/IhatWoypwFicya4d6vc9JvH7peOx9t\n6CshgpKiiIRbgEkxWI821KPDREo4yz86bKjJsh/q0WEiUhpYrA2MkqKIhJfFspDFwhGRUkc1RRER\nAyVFEREDJUUREQMlRRERAyVFEREDJUUREQOLZSGLhSMipY5qiiIiBkqKIiIGSooiIgZKiiIiBkqK\nIiIGFstCFgtHREodi9UULRaOiJQ6AXZxCrBk2VpadB5B0w7DmfLS7HzjF371M+26jqL9FXfz9+73\n8t2Pv/keBXO2AAAKtklEQVQNRzVFEQmvAKtmdrud0eOms2zBFGrXTKBjj/vo37sLLZvXd5W56or2\nDOjr6JPlj827uX7oRHauez8Y4YiIFJEAa4pr1m2jScNaNKhXgzJlohg8sDsLF690K1O+fKzrdXr6\nWRLiK/oNRzVFEQmvAKtmKQdTqVu7mut9nVoJrF63NV+5z778icf//Q4HDx/jm3nP+Z2vkqKIhJeX\npJi4ExJ3eZ/MZrOZmv1111zKdddcyvKVfzD0nilsW/Ouz/KmkuKSZWsZO+F17HY7dw7tw7gxg93G\nz5z7LVNf/oTs7GwqxMXy+vNjaHthI1MBi0gp5yULdWvhGJwmfeM+vnbNquxL+dP1fl/Kn9SpleB1\nMV27tCEz087RY2lU9XEY7bfi6jyZuWTuZDavfIdZ875ny7a9bmUa1a/Jj1++wIYVb/HPh4dw19gX\n/c1WRMQhwHOKHdo3Z0dSCnuSD5GRcZ45CxLp37uLW5lduw+QndO38q/rdwD4TIhgoqZoPJkJuE5m\nGq/wdOnUyvW6c4cW7D+Q6m+2IiIOAZ5TjIqKZPqU0fQa9Dh2exZ3DOlNy+b1efO9LwAYddu1zFu0\nnA9mL6NMmUjiyscy++0J/ufrr4DZk5lO73y4hL49O5n5TCIihWoD06dnJ/rkyTejbrvW9frRB27i\n0QduKtA8/SZFsyczAb5f/jv/m7mEn5b81+P4iRNzX3fr5hhEpGglJjqGYsNiDQP9JkWzJzM3bEpi\n5NgXWTJ3MlUqV/A4L2NSFJHgyFvhmDQpXJGYZLGk6DccMyczk/cfYeCwSXz0xjiaNKodtGBFpAQq\nxG1+weC3pmjmZOa/p37I8ROnuOfhlwEoUyaKNcumBzdyESkZLNZa2pZ9bGl2SBYU35PskCxJRIxs\nNsg+tjTcYXhki+9J9ocmyw4NzeewWI4WkVLHYucUlRRFJLyUFEVEDJQURUQMlBRFRAyUFEVEDCyW\nhSwWjoiUOqopiogYKCmKiBgoKYqIGCgpiogYWCwpWiwcESl1CvGUnCXL1tKi8wiadhjOlJdm5xs/\nc+63tOs6iraX3cWlvcewYVOS33BUUxSR8AowCzn7j1q2YAq1aybQscd99O/dxa2rFGf/UZUqlmfJ\nsrXcNfZFVi19xed8VVMUkfAKsKZo7D+qTJkoV/9RRl06taJSxfKA+f6jlBRFJLwCTIqe+o9KOeg9\n6ZntP0qHzyISXl6qZom/OgZvirL/KCMlRREJLy9JsVsHx+A06R338UXZf5SJcEREQiTAw+dg9R+l\nmqKIhFeAVbNg9R+lPlpESjjL99Gy3mTZduqjRURKA4udxFNSFJHwUlIUETFQUhQRMVBSFBExUFIU\nETGwWBayWDgiUuqopigiYqCkKCJioKQoImKgpCgiYqCkKCJioKQoImJgsSwU2nCyQrq0wFnsl0sK\nobh850ozi+1vFsvRFpOF5TaYCBlAdLiDKEIW28f8huOvX1WABx57laYdhtOu6yh+27DT+8wyLD5k\n5gzO16CaRnGTlTM4t6NzOGcYwv09K+xAAOWtLIj9Pm/dnkyXqx8gpmZfnp8+11Q4PmuKZvpV/Wrp\nanYmpbDjl/dZ/csW7nnoJe/9qh4zFVN4GFd8lGNIXA3drqDY/ConJkK3buGOomCKNGZnMgT3ZIjh\n/+B5RyvAj1/iz9DtkoCjlLwCrCmayU9V4yvyypT7+Oyrn4smHDP9qn6+eCXDB18NQOcOLTmRdprD\nR457n+lWH+8DHedvGWbH5ankJn7ro6yRRWqTiYnhjqDgijzmM3nebze83pBn3FYv5fKOy/OdSVzk\nfVxIvrOB7jNWFcR+n6slVKZD++aUiTJ/ptBnSU/9qq5etzVPmaP5yuw/8CcXVK+Sf4Zb8/zN+//C\njPO2LDPjnCs9CkdiLAdUzBl3Bt81xSzDX4udGylVssg9VDwDnMgZwJHwnMlyA64jAaLInxiNP3BF\n8b1UYvQvwP3GTH4KhM+kaLZf1ew8na94ne4rc0EFlacN4EyI5XLeJwD1cOxIqTgSpLcNZ9yJlBRD\nz7n+M3EcKjv7Qk/G8QOXlPM+ldxD6GggDogxzMfbOeS82zQL2IE1vsslRYCXewvS73NB+AzHTL+q\necvsP5BK7Zr5+14FsD1fmFDDY5KX06NWNWlSuCMouGIZ80r/ZcScQHOb2X6fC8pnUjT2q1qrRlXm\nLEhk1ozxbmX69+nC9BkLGXxDd1at3UzliuU9HjpbtTcxEQmfwuQFM/nJtZwCdCXqMyma6Ve1b8/O\nfLV0DU3+Ppzy5WJ4d/rDBfhYIiKBMZOfDh0+Rsceo0k7dZqIiAheenMBm1e+Q1xcrNf5hqzfZxGR\n4qDILw0UaWPvEPAX78y539Ku6yjaXnYXl/Yew4ZNSR7mElpm1jHA2l+3EVWtF/MXLQ9hdPmZiTdx\nxXraX3E3rS8ZSbd+D4U4wvz8xZx69CS9Bz3ORZePovUlI3nv46/DEGWuEaOncUHzG2lz6UivZay0\n31lZkSZFZ2PKJXMns3nlO8ya9z1btu11K2Ns7P3Wi2O556GXijKEAjETb6P6NfnxyxfYsOIt/vnw\nEO4a+2KYonUwE7Oz3LhJb9O7R0cKcDqlyJmJ98TJdO575BUWffwUG3+ewafvPRmmaB3MxDx9xkLa\nt23C7z++SeKiaTz0zzfJzLSHKWK4/dZeLJn7rNfxVtrvrK5Ik2JQGnsHkZl4u3RqRaWK5QHo3KEF\n+w+keppVyJiJGeCVtxYyqH9XqlWtFIYoc5mJ9+NPv+OGfpdRJ6fNWUIxiLlmjXjSTp0GIO3UGarG\nVyQqKjIc4QLQtUsbqlSO8zreSvud1RVpUvTUmDLlYGqeMp4be4eDmXiN3vlwCX17dgpFaF6ZWscH\nUlm4+GfuGdEPCLzJQ1EwE++OXSkcO3GK7v0fpsOV9/LhnPC2VDAT88hhfdm0dS+1Wt1Eu8tH8dLk\ne0MdZoFYab+zuiJ9Sk6RN/YOsoIs9/vlv/O/mUv4acl/gxiRf2ZiHjv+NZ578g5sNhvZZIf18NlM\nvOczM/l1/U6+/WwqZ87+RZdeD3Bxh5Y0bVwnBBHmZybmyS/M4qLWjUlc9Dy7dh+g58BxrP/xTSpU\nKOd32nCxyn5ndUWaFIu6sXewmW38uWFTEiPHvsiSuZOpUrlCKEPMx0zM69bvYPCdkwFIPXaSxcvW\nUqZMJP37hP4pBmbirVu7GgnxlYiNLUtsbFku79KG9RuTwpYUzcT885pNTHjwFgAaN6xFw/o12LZz\nHx3aNw9prGZZab+zuiI9fDY2pszIOM+cBYn0793FrUz/Pl34IOfwyFdj71AwE2/y/iMMHDaJj94Y\nR5NGtcMSp5GZmJN++5DdvzuGQf268vq0B8KSEM3GO6DPJaxYvRG73c6ZM+dYvW4rrQxPOgk1MzG3\naFqPZT/8BsDhI8fZtmM/jRrUDEe4plhpv7O6Iq0pFrfG3mbi/ffUDzl+4hT3PPwyAGXKRLFm2XRL\nx2wlZuJt0aweva/sSNvLRhERYWPksL60ahG+pGgm5vEP3szto6fRrusosrKymDppJPFVKvqZc/Dc\nfOcz/PDzBlKPplG39S1MemwY5zMzXfFaab+zOjXeFhEx0HNdREQMlBRFRAyUFEVEDJQURUQMlBRF\nRAyUFEVEDJQURUQMlBRFRAz+H24erxeO1RI+AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQcAAAEKCAYAAAAfNZB5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8FHX+x/HXbEnvCUgIEHovohSBoCAioBSlnKA0QUQF\nlR/qgdjg9BA41FPhTuU8FUTkkCYWFFSUKF0EKSLSCSCEJKRny8zvj0k22WRSCMlmIZ/n45EH2eyU\nD7uz7/nOd74zq2hJGzSEEKIQU1UXIITwThIOQghDEg5CCEMSDkIIQxIOQghDEg5CCENXfTgcP3kO\nU+TtqKparvmXrviGPkOmV3BVleNqqlVc/bwmHOK37qNrn8cJq38XkY0GE9dvCjt3H6rQdRgFyX3D\nevHVyjkVup6y6DHgCd5d8mWxz1d1rabI2zl6/KxH1nU5jp88R8+BTxJYZwAtOo/jm+9/LnUem81O\ni87jqNv63sta1keffEts2/sIqjuAu0fNJDklzfVcTo6NcZPnExo7iOgW9/Dav1a6zfvLr39wY89H\nCKwzgA63PsKefUfcnn/tXyuJbnEPobGDGP/oK9hsdtdzScmp3D1qJkF1B1C/3UiWrfzWbd5vvv+Z\n5p3HEVhnALcOeoqTp8+X+hqUh1eEQ2pqBv2HP8vjE+8m+dhqEvZ/zAt/HYWvr0+lrE/Tqn7cl6Io\nZZquKmv1htepsBEPzObGdk1IOrKSvz97P0PHvkjixUslzvOPN/9HzRphFH7JS1rW/oPHeWjq6yx9\n52n+/O1/BPj78siTb7jmnTl3CUeOn+Hkrx/x3dp/MO/N//HVNzsAPYwG3fcCo++5jZRjqxkz/HYG\n3fcCdrsDgK++2cHcN5bz7Zp5nNi7lKMnzvLCnMWuZU966k38fK2cP7SCpW9P5+En3uDAbycASLx4\niSFj/sbfn7mf5KOr6HB9U+4Z99IVv65GvCIcfj9yGkVRuGdwDxRFwc/Ph949b6RNywaAvpG+NH8p\n9duN5LpmwxjzyDxSUzMMl1W/3Ui3PcDMOYsZ9ZC+t735zqkAhDW4m5B6g9i64wDvf/QV3e/4P9f0\nP23bT8dekwirfxedbpvMlu0HXM/1GPAEz89+n7h+UwipN4g+Q6ZzMSnVsI6US+n0H/4sNZsOI6Lh\nYAaMeI6EM4kAPPPSf9m85VcmT1tAcL2BPDZ9YZH5y1KrKfJ2/v3fdTTpMIaQeoN4fvb7HDl2hi63\nP0ZY/bsYPv7vrg0S4LOvtnL9zRMJb3A33fo+zq8HjhnWnrfudjdPJLjeQFas+d5wOk/7/Y/T7P71\nD2ZNH42vrw+DB3SnbasGrFy3udh5jp04y9IV3/L0lBEUzLrSlrX0k28Y2K8LcTe1JjDQnxdnjGXV\nZz+SkZEFwOLlG3juyZGEhgTSvGk9Hhx9B+8v+xqATfF7cDpVHn9oMFarhUcfvAtN0/hu8y8AfPDx\nBh4Y1Y8WzWIJCw3i+adGuubNyMhi1Wc/8uKMsQQE+NHtptYMuqMrS/63EYBV6+Jp3aI+QwZ2x8fH\nysxpo9iz/yi//3G6wl9vrwiHZo3rYjabGDtpHus37nBrvgG8t/QrPvj4azatm8/RnxeTnp7F5GkL\nDJelKIrbXrng75u/eA2AS8fXkHpyLTd1bOk2b1JyKncOf5YpDw0m6egqpj48hDuHP+tWz7JV3/H+\nwqc4//sKbHYH8xesMKxDVVXGj+zLyb1LObl3Kf5+Pq6a//7sOLp3acPCeY+SdvJT3pgzqcj8pdWa\n5+vvdrF701ts/foN5r7xPyZMeY1l/5nByb1L+fXAMZat/A6A3Xv/YPxjr7Don1NJOrqKiWP7M/De\n592as3l++PxVAPZufoe0k58y7K5bikwTv3Uf4Q3uLvbnp237Deu9Evt/O07D2GgCA/1df2vXuhH7\nc/eqRh6dtpCXnx+Pn597K7S0Ze3/7QTtWjV0PdewfjS+PlZ+P5JAckoaZ88l0a51/vNtWzV0m7dt\nqwZu62vXOv/5A4fcl922VUP+PJ9Mckoavx9JwGIx07hhTP68rRqy/7fjrroLrjcgwI/GDWqz76Bx\n0F8JrwiH4OAA4r94DQWFCVNeo2bTYQy673nOX0gG9BR/YtJQ6terRWCgPy8/P56PV20qUydkwaZx\nac3kz7/eRrPGdbhvWC9MJhPDh/SkeZO6fPrlFkAPmvvv7UPjhjH4+fnwl7tu4ZdfjxguKyI8hLv7\nx+Hn50NQkD8zpo7g+x/3FltbSXWX5K+P/oWgIH9aNo+lTcv69OvVkfr1ahESEki/2zqye+8fALzz\nwedMHHMnHW9ohqIojB7eG19fK1t3HizTegqLu6k1ycdWF/vTtXOrci23JOkZWYSGBLr9LSQ4gLS0\nTMPpV38Wj6ZpDLqja5mWFRzkT3q63jLIyMw2Xld6pmuags/rz2WVUGcgaemZuc9nF5kXIC09i/SM\nLNfj/LoKLjubkOCidaVnZBu+BlfCUuFLLKfmTevx3sKnADh0+BQjH5rDlBn/5qNFMzh7LonYOte5\npq1XpyYOh5M/zydXaA1nzl2kXp2abn+LrVuTM+cuuh7Xqhnh+t3fz4f03GZmYZmZ2fzfM2/x1bc7\nXS2P9IxsNE1ztWbK2u9Qkutqhheox9ftsZ+fD+cvpABw4tSfLF6+kTcXrXU9b3c4OHsu6Ypr8JSg\nQH9SCwVByqX0Ih8m0Jvnf525iC//N7vMy7qUmkFwkL/r+UuFDl0vpWUQHBRAUO40qWmZREWGGs5r\nVGdwUEDu835uz+etJzjIv9i68v6PwUH+pKZlFHk+b90VyStaDoU1a1KXMcN7s+/gcQBqR0dy/NSf\nrudPnj6PxWJ2+yDkCQzwIyMzP0XPnc/f+Ev7MMZER3HilHvP74lT54mJjrrs/8MrCz/h9yOn2b7x\nTS6dWMv3n72CpmmuFkFptVREcCjkL6NenZo8M3WE2949/dQ67hnco1zL3rzlV4LrDSz258et+664\n/sJaNa/P0RNnXXtugD37jtKqef0i0x4+msCJU+fpfudUolvcw5Axszj750WiW9zDydPnS11Wq+ax\n7Nl/1PXckWNnsNkcNG0UQ3hYMNG1ItxajXv2HaV1i/x59+53b+b/euAYrZrHuv4fv+wrOO8RrqsZ\nTnhYME0bxeBwOPnjaEL+8/uPuM27Z19+XRkZWRw5ftbwNbhSXhEOhw6f4tWFn7g67E6dPs+yld/R\nJfc4e8Tgnrz275UcP3mO9PQsZrz4X4YP7oHJVLT869s04uNVm3A4nOzcfYiV6+JdH7QakaGYTApH\njp0xrKPfbZ34/chplq38FofDyfJVm/jt8Cn69+nsmqaszf30jCz8/XwIDQkkKTmVWfM+dHv+uhph\nHDluXEdZai2O22EU+WE0YfQdvPXeZ2zf9RuappGRkcXnX29z+3C41VczvMR1d+/ShrSTnxb70+2m\n1pdVd1k0bVyH61s3Yta8JWRn21i1bjP7Dh5nyMDuRaZt07IBp/d9xJ4f3mLPD2/xn9encl3NcPb8\n8BZ1akeVuqz7hvZi3fotxG/dR0ZGFs/Nfp8hA+JcfRSj7+nNS68sJeVSOgcPneA/H37J2BG3A9Aj\nrh1ms4k33l5NTo6NN95ejclk4tab2+fOexvvfrieg4dOkJySxovzl3L/vfq8gYH+DO7fjedf/oDM\nzGzit+5j3fqtjPrLbQDc3b8b+w4eZ9W6zWRn25g170Oub92Ipo3rVPjr7RXhEBzkz7Zdv9G596ME\n1R1Alz6P07ZVQ155cSIA40b2ZdRfbuPmO6fS8IbRBAT48ubcya75C+5lX5wxhiPHzhDe8G5mzl3C\nfUNvdT0XEODHM1PvpVu//yOi4WC27Tzo1oEZGRHCZ8te5JWFnxDVeAjzF67gs2UvEhEeYriuwp2f\nBU15aDBZ2Taimgyla98p9OvV0W3axyfezSefbiai4WCmPP2vIvOXVmvhWgzrI3/6G69vyqJ//h+T\npy0gouFgmnQcy+KPNxjWDjBz2ijGTJpHeIO7+WTtD8VO52kfv/sMO3/5nYhGg3nmpfdY+cHzREbo\n709eawbAbDZTs0a46yc8LBizyUTNGuGunUpJy2rZPJa3Xn2c+x58meua/4WsbBv/mv+Yq45Z00fT\nqH5tYtveR89BTzHtsXu4/dYOAPj4WFnz4UwWL99AeMPBLF6+gTUfzsRiMQPQp1dH/vroX+g56Cnq\ntxtJo/q1mTV9jGvZ/5r/GFnZNmo2G8bIiXN465XHadFMbzlERYay8oPneeal94hoNJidv/zOx+8+\nUymvtVLazV7GTZ7P5xu2UTMqjF9/XGQ4zWPTF/Llxu0E+Pvx/sKnaN+2caUUK4TwnFJbDvff14f1\nK14u9vkvNmzjj6MJHN75Ae+8NoWHn3i9QgsUQlSNUsOhe5c2hIcFFfv8p19uYcxw/Xipc4cWpKRm\nVPhZBCGE511xn0PC2YvUjanhelyndhSnz1y40sUKIapYhXRIFu7Br4jTcEKIqnXFg6BioiM5lZDf\nUjh9JtFwXIAS0ftKVyWEKCctqfgzU8W54nAY2K8LCxatZfiQnmzdcYCwkEDDwUm68t1zofIoQEkn\na2aiKDMxm6F2XZVOXRx88pEvt/axExGlclM3B1vjLfS83UFIqMZrL/uxcVtaCcszkFzca3X55syZ\nydzpj1fY8ipSMsb/zzlzZjJ9+syyLSS8aF/WpRSFNvVC2Z+QQqemodgdkJaqYLddQbGGCm4rM4mM\nep79CanU8g2r6BVVqJatHRzYV76rm0sNhxEP/J3vf9pL4sVU6ra+l1nTR2N36Ff6TRzbnzt6d+aL\nDdtpfOMYAgP8eG/BkyUszdsON0of0KRp4FT1jfCzVT7E773EAyOCiG3gZPdOC2dOm6gbq3JTnIOH\nRweSnKQQHuF9lzpfq7ZsttDuBgc/b7cQFqFy7owJp6P46fOOeAseCSuK++NCcxj8rmC25C3L27Zp\nd5YruOtBqeGw7D+lD7BYMO/RMq2saQtnmaarbJoGmgpOJ2RlKVxKUVBVcDr0v2nglhuaCtlZEFNH\npWUblakzsli0wBdfX0hJVggL17BYoHkrJ4cOmLgpzjv+n9XBTz/oLbedWy2ER2gcP6IUuW8DwAOT\ncpg5Lwt/fzh90sRXn1lZsdTKLzstOFVQNAgO1YiK0ggK1jCZwW4DhwMcDgVN08NAUeBiosqgYXZM\nJu/Zposz9F47e0u/H44hj154tfWA8b0PqpqmwdE/TMya7s/GL6zY7XpIQA/XNE6nQmi4fljU9gYn\nZ0+bCI/UyM5S8PPXk+S6WioXzpuAqtlg4uJ6VMl6r8SV1nzsDxPDRtr4fqOVSykKNrse5gXVjVWZ\ntyDL7fEDk3J4YFIOqgo5OWC1gqWMn4b4TR2J66Evz1u36YKeL6kxXwKvGD5d1RQFGjVReWtJRoFg\ngILhAGC16iEQGKSRna1gt+stjrxLPHz9IKfir5wts+oYDhkZSu77ATabgtlMkZb+DR2LD2uTCfz9\nyx4MAHE94spX7FVGwqGAn7db8PUDpZhXRdP0rc5hVzBbNMxmza0J63SCyeyBQoWLj4+Gw65gsYDZ\nrOl9B4X6D86f9+5+AW8l4ZDr3FmFSWMDcNgxPGZVFLDl6L+nXlIICNDw8dH3PHktjewsfS8kPCc0\nTCMlWW89BARi2Bm5b7ckdnlU63BQVTh+1MQ/5/jStVUIigK+vhoWg23JpOhNV4AL5xWCQ8DPX++4\nyj15Q2aGQkCgnKnwpMgaGhfO653C/v4aVp+irbeMDIWjfxTd1NPT4aFRAXRrE8KQPkGs/p+VLOMr\n2Kslj3ZINqkR6snVGdI0/cfhUMjKhBo1NbrfaufBx3L43xIfFDPYi95WEcUEdpvepLh4wURwqJbb\nlM3vAMvOzu+cFJ4RHqFxKVkhIkrDbNH7hTRNQS3QzaCqsG6llcen5bj+ZrPB8DuDiG2o9zX9ftDE\n4nd8+b8HA2jcTCUgQG+RXEw0kZEOqlNBMYHFomEywT2jbbw4P4umNat+my7JtJnl7wTzaDhs8ZKe\nXUXRN6KAwPyOqI5NQ2jfUT8llpOtD6IpeItKRQFH7gaXkQG+vvrvJpOGqub1RVxex5a4cv4BGinJ\nJqJ9VTQVzAZtYUWBX39xb04c3Gfmwp8mPv0uHZMJ2lzvZMgIOxfOKxw7YiIrUx+vEhGlEhQEZouG\n6lRwOvXDSF8/fSfw0/5Uw8NQbxEQqPHXyaVPZ8Sjm3JUDe/dq6anKVitJQ2GIb+jSytt4IzwNoXv\nRZydDT6+GoVvJlajpkaNmsWd3Sj6hteoee1uBB4Nh9fn+npydcUyWyAoSKNhE5UbOzsIDITRD+aw\nNV5/ORwGhxWQf8rS6gMOu4ZiUtC0/EE3iqnoRigql92m4OOjHyYqCqgGZys0DerVd39jbuzkJOmi\niQP7TLRsrbqmmzPTj0+W+nD+nIkWbZzc0NFBk+Yq/gEa2VmQmamQk63Q5noHve9w8M85fiiK5rWt\nh7ieJQwXLYVHwyE5qer7P/X+Br2l8PFiM8ePmHjnowx69bWz9n8+ZGboYxf0VkGB8fSa3r8Aeqel\nw6m4XjxF0acxmyUcPM1m01sADgegFP/612/o/oTFAn8ZaeOTpT48/3I2mgbTHvVn904LS1anEx2j\nsX+vmZ93mDnwq5msTH0cS0Cghp+fht2up0HqJcV1itsbZWeVvzaPhsPMud7XFfzt1xYevDeQLQdS\nOZNgIjtL33D0TsmCN2vFNbjG6qMHTOHjW0UpOjpPVC6HHXx8c88YacaHeooCNa4r+kTvO+zMnOYP\nZHPgVzOfr/Fhy/5LhOT2Mcb1cBDXo+Q97wtzvG+brihVvyuvYrfe7iAgUOPUCYXMjNxx9KXMoyj6\nWPzCG6LeD+G9e5FrkablHu6V9KYpeuuisBrXqSQl6u/X3t1mbopzuIJBSDiQmQmJ503Uqq3mDmrS\nSr3OTtMApbjBUtduB5W30vKP/Awp6KciC8vKzB+XEh2jcjZBgr2gah8OP35voU17ByaTgn+AhqW0\nMxbkh4Ooeq5OYEV/S4zeFk3TL64qLOmiQkSk/ma3u8HJvj2WAtfViGofDtviLXS/1UF2lt7hZLGU\nEA55eygNTBIOXiHvlLKrFVfM+2I0sC05SSE0XH9TwyM0IqNUjh2p9h8Jl2r/SiReMFErWsXhUDCb\n9fPeRtkgWeCd3MKhhIAw6ii22xR8C/RFRETqoyKFrtqHg9ms36vBW89Ti9KV5b0r7mK6glQ1/3S1\nkHAgKFgjLTXv6koFTZVWQnVhMms4HPnvdmam3u8kdNU+HAKDNDIyFBSTfi8A2TSuYiW8eUb36LBY\n8q+oBcjKhICAii/ralXtw8FiAdVZqIkpTYdqofBhhdOp38RH6Kp9OIBcQFWtacX8Ljw7fNobOZ25\n92TQytaxJa4dhcerKKZyDH+vwO8d8TbVvuUgA5quIeV4HwvuEArem0NIOADSYhA62Q7cSTgIIQxJ\nOCAdkkIYqfbhIE1JIYxV+3AQQhiTcBBCGJJwEFeFs2fP8OwTz9K1dVfqBtWlTb02TBo7ifS0s1Vd\n2jWr2g+CEleHvXt/Zv2n6xk9YTQ3dr6RP8/9ydyZczl/7lbGPbKTQOSiiIom4SCuCl26dGf7oe2Y\nCnzRRLsb2tGxaSd+P7CGznH3VmF11yY5rBBXhZCQULdgAGjUpBFWawDp6eeqqKprm4SDuGrt37sf\nuz2TiMjGVV3KNUnCQVyVVFXl6cefJiKyCY2b9a/qcq5J0ucgrkp/e/pv7Nq2iyH3foXJZEaut654\npbYc1m/cQfPO42jSYQxzX/+4yPOJFy/Rd+jTXH/zRFp3ncD7H31VKYUKkefdf73LgvkLWPjBQqJj\nOlR1OdesEsPB6XQyedoC1q+YzYEt77Js5XccPHTCbZoFi9bSvm1jfvnhbTatm88Tz72NwyE3/xcl\nW7NmBUOG3M6QIbezZs2KMs/36cpPmf7YdGb9YxZ3DburEisUJR5WbN91iMYNalO/Xi0Ahg/uydov\nt9CiWaxrmuhaEezdfxSA1LRMIiNCsFjkFr6ieGvWrGDq1IdISUkCYPfuXSiKwqBBQ0ucL35TPA+N\nfIgHH3uQSVMneaLUaq3ElkPC2UTqxtRwPa5TO4qEs4lu00wYfQf7fztB7Zb30O7mibw++5HKqVRc\nM5YsWeQKBoCUlCQWL36nxHkOHTrIyLtG0qR5E+76y13s2LqDHVt3cOb0dlKSjlV2ydVSiS0HpQyX\nLM5+dRnXt27EpnWvcOTYGXoPnsaeH94mOLjoiLU5M+e4fo/rEUdcj7hylCyqo59/3k5aahr79+6n\nb9e+rr9rGrS5fiTDRr5VhdV5l/j4TcTHb7ri5ZQYDjHRkZxKuOB6fCrhAnVqR7lN89P2/TwzVR+d\n1qhBbRrE1uLQH6fo0L5ZkeVNnzn9igsWV79Royawe/cuV+shLCyC0aMfLHGeESPGMOKRgUX+/vIL\nfvq3bMvZCpe4uB7ExfVwPZ43b1a5llNiOHRo34zDRxM4fvIctWtFsnz1JpYtmuE2TfMm9dj4/W66\n3dSaP88nc+jwaRrWjy5XMaJ6uOuuYSiK4jqUGD36wVL7G4TnlRgOFouZBXMn02fo0zidKuNH9qVF\ns1jefv8zACaO7c+MqSO4f/J82nWfiKqqzJs1gYjwEI8UL65egwYNlUDwcqUOgurXuxP9endy+9vE\nsfkj0qIiQ1m37MWKr0wIUaVk+LQQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDC\nkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwlC1\n/yLd5q2cOBwQEAC39bNz+qSJc2cVkhIVsrMUNMBihpg6Ku07OQCIqavSoYuDoCBIS1UICdVvi37L\nbXYiouQW6Z7UorUTkwl8feHCn06ysyDpokJqioItB0wWCPCHOvXUIvPG1FXp3M3hetx3gJ2AQHn/\n8iha0gaPvBpKRG+StKTSJ6xuksMrdHHhJFfo8ipKMhXw/wz3wv9bBb9/lSEiQkFL2nDZ81X7w4p/\nzvHlHy/6kXBa4YbGITSMCCXaL5RIUxgRSv5P3eBQurUJBmDdKiuDewdyz51B3NYpmIP79Zfx7t5B\nbPtJvifUk175ux+vzvbj7Td8Gdo3kNjQUGr5hrq9d1HmMD5daS0y72errUwcmf/NbDe1DOHUiWr/\nkXCp9ocVdruCzaZ/rVp2loLNpuBw6o8LUp1gy9G/HtDhAJtNQVUhK0tBy22xZhf4XXiGzQYmE9jt\n+nvicICquX+No6bp71lhTidkZ+dPm5VZ9H2vziQmhRCGJByEEIYkHIQQhiQchBCGJByEEIYkHIQQ\nhiQchBCGJByEEIYkHIQQhiQchBCGJByEEIYkHIQQhkoNh/Ubd9C88ziadBjD3Nc/NpxmU/we2t/y\nEK27TqDHgCcqvEghhOeVeFWm0+lk8rQFbFw9l5joKDr2msTAvl1o0SzWNU3KpXQmPfUmX33yMnVi\napB48VKlFy2EqHwlthy27zpE4wa1qV+vFlarheGDe7L2yy1u03z0ybcMGRBHnZgaAERFhlZetUII\njykxHBLOJlI390MPUKd2FAlnE92mOXwkgaSUNHoOfJIOtz7CkuWXf8cZIYT3KfGwQlGUkp4GwO5w\n8POeP/hmzTwys3Lo0ucxburQgiaN6lRYkUIIzysxHGKiIzmVcMH1+FTCBerUjnKbpm5MDaIiQvH3\n98Xf35ebu7Rhz76jhuEwZ+Yc1+9xPeKI6xF3pfULIQqJj99EfPymK15OieHQoX0zDh9N4PjJc9Su\nFcny1ZtYtmiG2zSD+nVl8rQFOJ1OcnLsbNv1G1MfGWq4vOkzp19xwUKIksXF9SAurofr8bx5s8q1\nnBLDwWIxs2DuZPoMfRqnU2X8yL60aBbL2+9/BsDEsf1p3rQefW/tSNu4iZhMChNG30HL5rElLVYI\ncRUo9Qaz/Xp3ol/vTm5/mzi2v9vjJx8dxpOPDqvYyoQQVUpGSAohDEk4CCEMSTgIIQxJOAghDEk4\nCCEMSTgIIQxV+3CQ70YUwli1/yJdTYMil5BIYFQLFbJjCE+ugIV4p2rfclBV/Vua85ThWjNxDSn4\nfiuKvj0IXbUPh+wsBT8/DU3N3UquYG8ihyjeSzP40BfeMfj6Qk6O52rydtU+HM6dUYisoeFwgNms\noZhKzwdNc59G0xTX35GWh1dS1aJvjC1Hwccn/50MDddISqz2HwmXavtKaBr8+ouZ7762EtfDQWaG\ngq8fWCwFP/lKkXlA3wuZFL0Zqiia2/NyWOJZhocChdJd08BmKzpv0kWFsPD8ia+/0cGWzcbdcNWx\nVejRDsm2sSGeXJ2bvL27qoLTAZdSFGrXURk/KYcGjVS+/txCSKhGZoZCiquPqcAWUWAjtNuV3FZG\ngRYDRZupovKZTGC357/uioJh6y01pegff9lpplsPh+vx2Ik5DLo1mPAIlUZNVeI3Wfh2vZVDB8xk\nZCgEBGgEBWv4B2gMGWHjmZeyaRtb9LaIBXcYVW3qM9nlntej4fD5D+meXF0RiqKBAhYLhIRqBATo\nf09Lg1dn+xMSqnHhTw2rj4LTUXSP5HTq/2ZlgtVHQdP0jTJvOqdDX7bwHKsVsjIVvcWH3qIrTFHg\n4H6z29+OHDbx3ddWHpqSzYYvLMTUU2nZRmXxqgze/Zcvyz4w0e0WBzPnZdH2BgdBQZCZCelpClmZ\nCkHB+vq+2JxWqPXoXU3H0DCVqRPLN69HN+UP/+vjydUVS9NAdUJ2tkLyRYVtP1poe4OD+E1WMtIV\nLBZw2N3nMSn5TdPECyYCA1VsNgWLRcGRu/PJylLw9fOevUZ14B+gce6Mgp+/hqqC2QKKQafihs+t\n2O16mGRkwFOPBNCpm4PBvYO5sbOTg7+aadLcyQ2dnHTq6qDbLRoZGQpbNlv4fqPF1SpUTHrYtO/g\npE9/u9ds08Xpebu99ImK4dFw2BrvBbvV3C6FvHx3OsHHV+O7r60EBGr4+kJGRqH0V8Bkhsx0/e+H\nDpiJiFLD5R/aAAAU90lEQVRJS1Xw8QVb7saYkqQQESnh4Ek1aqps/8lCeIQDew5YrOBw6B/ggv0E\nKSkKE+4NpHtPB2/+wxebTT87ERKqsWWzhaBgjYwM2LnVDJj1MDAbtEQUQAMFjdvv9JJtugTRMeU/\nN+vR/9nOLVX3Qhb5yGqgaRqqqu/5g0M0Ei/opypU1X16kwl8fPSNbeXHVjZtsHD/wzn4+KokJymk\npSmkpelNzqgaEg6eFNtQ5dgfJurVV0m8YEJ15gZEgR2mpukBvm6llS/WWlHQdwq+fmC36f+mXVL4\nZacFqw+YzaBvIEr+mam8t1XJX6am5bDDYJv2pgOL2PrOcs/r0U+r3VH6NB6jgZ7/+l4mLVXBbMk/\nH567g8Bk0n/8AzTuvNvOhBFB9L/bxrkzJtq2d3IpWeHCnyZ++t5Km/ZO6XPwsHr1VU4eN9G8lZML\n501omqY3/wu1HPJ+Lxgathx9R5CRnnvmyaT3G6nO3L6DEj7lSRf1eR3etE0bOH++/FHl0U258HG8\n18hNAoc9v7dbyf2xWDUcdoXrovUzGVsPXiI4WOOW9iE8+Ww2OTmw8BU/0lLhhblZVf0/qXZqXqcR\nEAhHfjcx6oEctvxg4fQpE8kXFbfOYiNOJ/p7n/v5UVVQy9jwS0kyoWlevE3nSrxawsFrFd7D5HZK\nmEzg7w+1Gjlp297B6uU+RERqfLHGSr9BdmIbqEx8LIf6DVWO/G7irmFevqVcoyY+ns3jEwL5YGU6\n6z8NxseqYTYrZdura27/lFn2VTKSMjNTwqFC5R1KBAZp+PlB6iWFvgPs3DfOxuerrTz3chb979aD\nwGrF9buoGo/8Xw67tlm4b2AQT7+YxbxZ/oTYNJKSlPzjwwpmNus/vn4Vv+yKdENHJ7/tK9+8ng0H\nb+qpKUQh/1DCatXw9dc7IdNSFWbOy6b3Hfpu6IaO5e/gEZXDZIL/Ls9g6Xs+/PtVP04dNxEUrBIS\nCpkZYLcpxY5wLDyiVSv4Swnbq5+/htkMNa5TOX3CXPyEVeyHb63lntej4eDnjSmr5XdDq5reIWmz\nKYSEqrRsqzJ9Zhadujpdk65abqVlGyctWsnle95EUWDkOBsjx9nYt9fEP2b5sXObhdQUE1YfMCla\ngQaEkjtwKbfXMq91oUBZ92B1Y/Wl1avv3eGg11e+eT0aDi+9kunJ1ZWJqio4nfoAp4x0hf17zBw6\nYGbewky2bLYw+zl/1nyjj+z8+7N+LF7kS2iYRvyvqfj6VnHxwtDsZ/2p30jlvvGZTBoTyP0P6/1C\nwSGaa7h1TrZCVqa+I8jbP+QdTppM7s0Mo1ZH0xb6ziGuh4Ofvi//3rmy6fWVb15FS9rgkRPzSkRv\nkrQkT6yq3JxO2POzmf43BzNspI2HpmTTr1swx1MuYbfDdT7hvPx6JvNf9OPfSzK4rW8FnMdKDr/y\nZRQQjnfefCSZCvh/lvHGKrGhYaxYn8ann/jw1TorX29NIzyi+o4/iVAi0JI2XPZ8Hm05DOkT5MnV\nGXJdWZl7VsLphJwchZRkhTOnTNSqrfLXF7K4qbuDqRMDmTpDv3DFaoVPv0tj2qMBDBtpo1cfLz/B\nXY3NeSOTyWMDeem1TLIyoVOzEHr1tVOnnorVh9yL6xTOn1NITjLppyMV8PXV8PHVOxrzrvYsrq+i\n9x12Hno8h6F9q36bLsnoCeU/reLRlkO3Wwyum60KeYeaiushSu5VfRlpCieOmQgK1rhntI0nn82u\n3IFN0nIou8u4JdsHi3xY9KYviedNNGzqRHWC05l3sZyGyQwWM/kDpsgb9Zi7gNw+CMNuCA3ietr5\n6/M5fPu1d5/wa9xUpX3DMO9vOdzSy8tO+eW+6SaTPs7eP1AjMkqlRWuVxk1Vufz6KjZmgo1R423s\n3W1m1zYzCadNZKQpqKo+vNrXT8PfH6w+mr6T0MCp5o+QLXxvjsL50LyNE0WBHVsszJvl76n/1mX7\n6wvlH5jn0XBYv847r2BTVXI7JfVOqqSLJsxmjb/Nz2LUeC9p7YjL8sVaK4+OC0BRoMZ1GsEhGj65\nQaCqeke03V40DAzvB2HQtjaZ7XDXtX1oKR2SxfjqMyvPPeHP9kOplbsiOawos3eXzefrz79m59ad\npCSn8Ol3n9Ltlm6G0w6+PYh+A+1MmHyVDGWsRFdFh+Tvv3lPOz3v+FK/t4NCdrY+4OncGRN7fjbz\nyVIfBg71ssOgam75kuWYTCZ69e3FymUrUUq4J1+X7g5mP+fH2QSFm3s5aNREJTRcxWzWB0VlZuod\nk/rFV/pyTGbN1RlZmpBQjVrR1/YZEI+Gw6i7vKxnN7cZaTaDn59+C7CatTRatnHy+eY0mjaXgU7e\n5OstXwNwcP9BVi5bWeK0Tz2XzZARNj5814d5s/w4dcLMpWT9egsfX/0uYP4BGr5++o2FNU3vsFRz\nB8DmdUwWFxSD/mJjxt+ymTPT75rtc5DDiqomhxVll3u24sC+A3Rv2511m9bR9eauV77ca1x5Dyu8\np50vhPAqEg5CCEMSDkIIQ6WGw/qNO2jeeRxNOoxh7usfFzvdjp8PYanRh1XrNldogUKIqlFiODid\nTiZPW8D6FbM5sOVdlq38joOHil7/6XQ6mTbrP/Tt1bFafjOQENeiEsNh+65DNG5Qm/r1amG1Whg+\nuCdrv9xSZLo331nL0IHdqRFZ9Nt/hDCyZs0Khgy5nSFDbmfNmhVVXY4wUOI4h4SzidSNqeF6XKd2\nFNt2/eY+zZlE1n75E9+u/Qc7fj4k3xUpSrVmzQqmTn2IlBT91Pbu3btQFIVBg4aWON/unbs5efwk\nCacSAIjfFM+F8xeIbRDL9TdeX+l1VzclhkNJI9DyTJnxL+Y8Px5FUdDQ5LBClGrJkkWuYABISUli\n8eJ3Sg2H/yz8Dx9/oPd7KYrC3JlzARgxdgQL/rug8gqupkoMh5joSE4lXHA9PpVwgTq1o9ym2bXn\nMMMfmA1AYtIlvty4A6vVzMB+RQenzJk5x/V7XI844nrEXVHxonpZ+N5CFr63sKrL8Hrxm+KJ3xR/\nxcspcYSkw+GkWaf7+WbNPGrXiqTTbZNZtmgGLZrFGk5//6R/MKDvTQwe0L3oimSEpLFqOEKy8GFF\nWFgEr732dqkth8u5n4PIVykXXlksZhbMnUyfoU/jdKqMH9mXFs1iefv9zwCYOLZ/+aoV1dpddw1D\nURQWL34HgNGjHyw9GITHefbaiiTpkKhsV0PLQXhWRIQi11YIISqOhIMQwpCEgxDCkISDEMKQhIMQ\nwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQ\nhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISD\nEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCUJnCYf3GHTTvPI4mHcYw\n9/WPizy/dMU3tOs+kbZxD9Kt7+Ps3X+0wgsVQniWpbQJnE4nk6ctYOPqucRER9Gx1yQG9u1Ci2ax\nrmkaxkbzw+evEhoSyPqNO3hwymts3fBmpRYuhKhcpbYctu86ROMGtalfrxZWq4Xhg3uy9sstbtN0\n6dSS0JBAADp3aM7pM4mVU60QwmNKDYeEs4nUjanhelyndhQJZ4v/8L+7ZD139O5UMdUJIapMqYcV\niqKUeWHfbf6F/y5dz4/r/2n4/Jw5M12/x8X1IC6uR5mXLYQom/j4TcTHb7ri5ZQaDjHRkZxKuOB6\nfCrhAnVqRxWZbu/+o0yY8hrrV8wmPCzYcFnTp88sf6VCiDIpvOOdN29WuZZT6mFFh/bNOHw0geMn\nz2Gz2Vm+ehMD+3Zxm+bk6fMMHj2LD9+aRuOGMeUqRAjhXUptOVgsZhbMnUyfoU/jdKqMH9mXFs1i\nefv9zwCYOLY/f5u3hOSUNB5+8g0ArFYL2zcuqNzKhRCVStGSNmgeWVFEb5KSPLKqai2c5KouwVAy\n4VVdQrUVEaGgJW247PlkhKQQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISD\nEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDC\nkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwpCEgxDCkISDEMKQhIMQwlCp37JdocK980tePS5ZvlS2\nXGT78ShpOQghDEk4CCEMSTgIIQxJOAghDHm2Q1KIKmKzQdJFBYcdAoMgJFTDbK7qqrybR8OhVUyo\nJ1dXJooJLBYNXz8IDdOIqatSv6FKi9ZOut9qJ7q2VtUlinKw2SB+k4UfvrHy7VcWDv9mJjRMw2KF\njHTISFcIDdMIDdMwW0BT9XnsNgWnM385ipL/czV68rmscs+raEkbStz612/cwZRn/o3T6eSBUf2Y\n9vjwItM8Nn0hX27cToC/H+8vfIr2bRsXXVFEb/ad9q5TUZqmbxQOh0J2NqQkKZw+aeLEMTO//mIm\n/jsLz8/JYtR4W8VuHJV4KjMc73qN8yRTAf/nMp7K/PF7C4+NDyAiSqNXXzu39HLQvqMDX19QVTCb\nweGAi4kKaakKDgeYTODjA1YfvUWhKLnbR4GfqzEggkM0YkMi0JI2XPa8JbYcnE4nk6ctYOPqucRE\nR9Gx1yQG9u1Ci2axrmm+2LCNP44mcHjnB2zbeZCHn3idrRveNFze3bcFX3aBlU1R9I3F11cjOESj\nZi2N1u0cPDc7i20/bua12f2wWmHEGFtVl1qq+PhNDIhrV9VlXJb4+E3ExfWokGU5nTB3lh9L3/Xl\nlbcy6djFweJFvrw4w4/fD5pJuqigKHoQhIRqhEdoBAbpYaCqYLfrQaHmthzcWgwFgiEzYxMBgRVT\nc2V79Knscs9bYjhs33WIxg1qU79eLQCGD+7J2i+3uIXDp19uYczw2wHo3KEFKakZ/Hk+metqFt1T\nZGaUu85KkxcOdhs4nAp2O5w5beX1uX7UqbeVfy+5mXsHBNF/sI1g78s2N9U5HDQNHhsfwIljJr7d\nlcpnq3yYfH8I9RuppKfphxGBgRqKST90SE9XsFg1V0tC1fRQsNv1oNA0MOWFg5KbDbkBkXrpe+DK\na/aEc2fK39wpMRwSziZSN6aG63Gd2lFs2/VboWkuFpnm9JkLhuGQluqd7TK92ahgtyuoTqhZS8Xp\ngLNnTMx/0Y8mzZ1s+9HCbX0dnP9TISxcw8enqqsWBS37wIe9u818tSWNvz3tz8YvLCgmOHncROol\nBYtFIzNTQdMUzGaNgECNnGyF7CwFNFwf/LIcOuTkKF67LRd2NqH8JyRLDAeljAdZmubebVHcfAvf\nzyxjWZUv7zjS6dT3JBnpkHDaxKjxNnZsNXP9jU4WvmJj5PhsYuqqREbp/8e8Y1PhXfrfbePmW+0E\nBMDIcTamzsji8G8WkpMUbDaoFa3y5ac+BAbq7+PAoTZOnzRxMdGE3a63Hvz89I5pS+6nonCfA+jh\nsepjO4OHe8+2XJLmrZy891b55i0xHGKiIzmVcMH1+FTCBerUjipxmtNnEomJdp8mz6i7vbxdDrzy\nkvvj99+eWzWFlNO8ebOquoTLVhU1z/vblc2/evnLFVOIFysxHDq0b8bhowkcP3mO2rUiWb56E8sW\nzXCbZmC/LixYtJbhQ3qydccBwkICDQ8pytNbKoSoOiWGg8ViZsHcyfQZ+jROp8r4kX1p0SyWt9//\nDICJY/tzR+/OfLFhO41vHENggB/vLXjSI4ULISpXqeMchBDVU4VfW7F+4w6adx5Hkw5jmPv6x4bT\nPDZ9IU06jKFd94ns3vtHRZdwWUqrd+mKb2jXfSJt4x6kW9/H2bv/aBVU6a4srzHAjp8PYanRh1Xr\nNnuwuqLKUu+m+D20v+UhWnedQI8BT3i4wqJKqznx4iX6Dn2a62+eSOuuE3j/o6+qoMp84ybP57pm\nw2jTbUKx01zu565CwyFv0NT6FbM5sOVdlq38joOHTrhNU3DQ1DuvTeHhJ16vyBIuS1nqbRgbzQ+f\nv8re+Hd47smRPDjltSqqVleWmvOmmzbrP/Tt1RGtCtuGZak35VI6k556k3Ufvci+nxbxyfvPV1G1\nurLUvGDRWtq3bcwvP7zNpnXzeeK5t3E4nMUssfLdf18f1q8ovpO0PJ+7Cg2HgoOmrFaLa9BUQcUN\nmqoKZam3S6eWhIYEAtC5Q3NOn0msilJdylIzwJvvrGXowO7UiKza61nKUu9Hn3zLkAFx1MkdLxN1\nFdQcXSuC1DR9VF9qWiaRESFYLFV3JVf3Lm0IDwsq9vnyfO4qNByMBk0lnE0sNI3xoKmqUJZ6C3p3\nyXru6N3JE6UVq0yv8ZlE1n75Ew+PGwBU7TUBZan38JEEklLS6DnwSTrc+ghLllftma2y1Dxh9B3s\n/+0EtVveQ7ubJ/L67Ec8XeZlKc/nrkKvyqzoQVOV7XLW+93mX/jv0vX8uP6flVhR6cpS85QZ/2LO\n8+NRFAUNrUoPK8pSr93h4Oc9f/DNmnlkZuXQpc9j3NShBU0a1fFAhUWVpebZry7j+taN2LTuFY4c\nO0PvwdPY88PbBAcHeKDC8rncz12FhkNFD5qqbGWpF2Dv/qNMmPIa61fMJjysagdylaXmXXsOM/yB\n2QAkJl3iy407sFrNDOzX1aO1QtnqrRtTg6iIUPz9ffH39+XmLm3Ys+9olYVDWWr+aft+npl6LwCN\nGtSmQWwtDv1xig7tm3m01rIqz+euQg8rCg6astnsLF+9iYF9u7hNM7BfFxbnNhtLGjTlCWWp9+Tp\n8wwePYsP35pG44YxVVJnQWWp+ejuJRz7Rf8ZOqA7/57/WJUEQ1nrHdSvK/Hb9uF0OsnMzGbbrt9o\nWeDiPk8rS83Nm9Rj4/e7AfjzfDKHDp+mYf3oqii3TMrzuavQlsPVNmiqLPX+bd4SklPSePjJNwCw\nWi1s37jAq2v2JmWpt3nTevS9tSNt4yZiMilMGH0HLZtXXTiUpeYZU0dw/+T5tOs+EVVVmTdrAhHh\nIVVW84gH/s73P+0l8WIqdVvfy6zpo7E7HK56y/O5k0FQQghDcoNZIYQhCQchhCEJByGEIQkHIYQh\nCQchhCEJByGEIQkHIYQhCQchhKH/B3QQRv9NUZYAAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAQcAAAEKCAYAAAAfNZB5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlYVGX7B/DvKEiIiGwqO66A4vaGC4ampokmmLjvW0ru\nvi4vLj/NpUxJLZXKpZLKJcVdS0rMtdwyxTQlBVc0FXFB1JCB3x/UxMhhzjAcZp4Zvp/r4rpgzvOc\nc58zZ+55thlUuel7ckFE9IIypg6AiMTE5EBEkpgciEgSkwMRSWJyICJJTA5EJEmY5DBi4hK8u3Bt\nodvLOL+OlCu3jBhRQbHrvkeLjv81aQyGOHTkN/g3HWLqMMjMlEhy+GbzPjRtOwYVvMJQxa87mrUb\ng0+/2KmzzqeLxuH/JvUtiXB0GjQqGjPeizX6cQ2lT7wvJtIWwfVw4dgXJRvY33wb9MOPB08Z5VhF\nkX7/Ebr0n4UKXmHwbdAP6zf/qFe9196cjDLOryMnJ0fvfe098Cv8mw6BnWcY2nSejGs37mhtj5q1\nCi41u8KlZldMmf2Z1rYr1/5E6/BJsPMMQ0DTIdh74Fet7es2/Qif+n1RwSsMXfrPwv0HGZptf/2V\nhSGjF8LBpzPcAnriw082a9U9/dslvNx6JOw8wxDUZiQSzybrPHfFk8OimDiMn/4posb1wO0LG3E7\nKQ7LF4/DT8fOISvruWSd/BeelJGba5q1bSqVymTH1mXU5GV4ycYad5LisHbFFIyYuBS/X7iqs87a\nuL3IzlZDpVLpva+0ew/RdeAcvDd9MO6nbEFQw9roOeRdTd0VsbuwffcRnDm0AmcOrcDO+KNYEbtL\ns733W/PwcoNaSE/ejPf+bzC6DZqLtHsPAQDnzl/B2xOWYO3Kqbh9YSPK29pg5KSlmrqzFnyN5Cs3\nce23ddi3/QNEL9uI7/eeAABkZT1H577vYEDPtnhweSsG9nodnfu+g+fPsws9f0WTw8NHmXhnwVf4\ndOFYRIS1gJ2dLQCgYb2aWLNiCsqVswaQ9+43YuISdOwxDRW8wrDv0OkC74gfLN0I9zo94Vm3N75Y\nE6/zuK3CJmLqnM/RtO0YOPh0xpv93tHKqN0HzYFbQE9U8n0Tr3aaoHkiV8Z+i3Wb9iF62UbYe4ej\nc9+ZAIDrN+4gYsAsVK7dHS41u2JMVIzW8SbPXAmn6hGo3qg/4hNOaJ3/0DGLNHHPeC9Wk/gupaTi\n1U4TUMn3TbjW6oZeQ98r9HyKGm9+Ld+YAABo0DIS9t7hiNt2APsPJ8IrsI+mjG+Dfli4LA71Q4bD\n3jscQ8cswu0799Gh+zQ4+HRGuy5RePDwsab80RO/o3n7cXCs1gUNW0biwE+JknH3f3s+rt24g7A+\nM2HvHY6Fy+IKPUdjysx8ii27fsLcaYNQvvxLeKVZIDp3bI6vNyYUWufho0zM+WANomcN00p2cvva\nsvMwAgN80TW8BcqVs8asqP5IPJeCPy7dAAB8uX4PJo3qBnc3F7i7uWDS6G6IXfcDAOCPSzdw6rdL\nmD1lAGxsyiEirAXq162GzTsPAQDWbtqL8A7BCGkWCDs7W8ydNghbdv2EzMynAICvNuzBjEn94FDR\nDv61vTF8QEfErs/b9/7DiVCrczDu7QhYW1thzPA3kZubq7OVp2hyOHL8d/z113N07thctuz6zfsw\nY1I/PL6+EyHNAqGCCv8k6PiEE1j0ySYkbInGHydWI+GFppWUrzckYHXMJNw6vwFWZcti7JSPNdve\neL0pLv0Si7sXN+E/9Wuhb+T7AIDhg95A325tEDW2BzKu7cD2tXOgVqvRqfcMVPOuiquJa5B6bj16\nR7TW7OvYyQvwr+WFe8mb8b+xPTB03CLNtkGjPkC5clZIPvkVTh34FD/sO4nPvtoNAJgxLxahrzXG\ngyvbkHpuPcYOf7PQcylKvC86+O1iAMCZQyuRcW0Hur/5aoEyKpUKW3Ydxt5t0Ug69gV2/XAUHXpM\nw/x3huLOH3HIyc3B0hVbAQCpN9PQqfcMzJzcD/cvb8XCOZHoOnCO5t1M6zlYPgXenpWxa/1cZFzb\ngUljuhcoc+3GHThW61Lozzeb9xV6XQz1R3IqrKzKomZ1D81jDepWx7kLVwqtM23u5xg5JAxVKjsW\naV/nLlxBg8Dqmm3ly7+EmtXcNdt/T7qKBoE1NNvr162Oc38n/3MXrqC6j5vmTRUAGgTWyLf9KhrU\n/Xff1X3dYFPOGn8kp+L+gwzc+jNd69ja+76K+nWraZ1Lg8B/t0tRNDmkpT+Ei7MDypT5d7f/vOOU\n9+iEw0fPah5/843mCG5SBwBgY1NOaz8btx3AkL7tUcffB+XLv4TZUwboPK5KpcKAXm015edOG4iN\n2w5qMv6gPu1hZ2cLa2srvBPVH4lnU5CR8URTP38r+PjJJNy6nY4P5gyHra0NbGzKoXnTuprtPl6V\nMbR/h7xj9myHW3+m487d+7h95z52JxzHh++NgK2tDVxdKmH8iC74Zut+AEC5cta4cu02Um+moVw5\na619vqgo8RpqzPDOcHWpBHc3F7RoVg/BjQPQILAGbGzKocsbr+DUb3n90TVxCejYrglC2zYGALRt\n9R8ENayN7/YcN+i43p6Vcf/y1kJ/enVtLb+TInqc+RQV7ctrPWZfoTwyHj+VLP/LqSQcOXEeYyQS\nuNy+Hmc+Q0V7O63tFe3zb38Kh4p2WtseZ0pvy9u3LR7/XTfzybMC2/P2/URT5sV9F3bcvO12hV4D\nALAqdIsBnB0rIu3eQ+Tk5GgSxM/fLwEAeAX20TSxVVDB09210P3cup2Oxv/x0/zt7VlZ9theHv+W\n8fasjOfPs5F27yGcHO0x/d3V2LTjEO6mPdDElZb+EPYvPMkAcD31Lny8KmsluPyqVnbS/F6+/EsA\n8m6ItHsP8fy5Gm4BPTXbc3JyNLFHzxqGGfNi0aTdaDg62GPiqK4Y3De0wP5zcnIwbe4XesdrqCqu\n/74j2tqW0/r7JZtympvt6vU7iNt+EDvjj2q2Z6vVaNOyoWKxlLQKdrZ4lC+5AnndBvsKtgXK5uTk\nYOTkZfho3gite+CfN5rC9vVPwrCvYItHGZmFHiuvfqbWtgp2+bcVHmcFO1s8fPTCvjMyYV+hPCr8\nXeZRxhO4ODsUctzC45aiaHIIblIHNjbW2PbtT4gIa2HwftyqOGmN8L442ivlxfLW1lZwcXbAmo0J\n2LH7CPZui4aPVxU8ePgYTtUjNO++Lw42eXm44tqNO1Cr1ShbtqzeMXt5uMLGxhr3kjdLJpYqlR2x\n8qO8adCfjp5F24govPpKA1T3ddMqtzZub5HiVUphrRFvz8ro36OtJnY5cuFdu3EHdZu/Vej2lR+O\nR++ubfQ6lr5q1/BAdrYal1JSNd2BxHPJCAzwLVD2UcYTnDx9ET3/HhNSq/Pe0DwD+2DT6hloWK+G\n5L7q+vsAAOr6++LLb37Q7C8z8ymSr9xCXX/fv7f74PRvyQhqlPfml3g2RRNHXX9fpFy9hcePn2pe\n7IlnU9C/R1tN3cRzKZp9J1++iaysbNSu4QE7O1u4VXXC6d+S0bbVfyT27YPFL8xenDmXgjHDOhd6\n3RTtVlRyqIB3JvfHyMnLsHnHIWRkPEFOTg5O/3YJmU+eacrlouCdmItczQ3a481XEbv+B5xPuoon\nT55hdvQancfNzc3Fmo17NeVnvv8lunduCZVKhceZz2BjYw2nSvbIzHyKaXO1p/SqVHZEytV/p/2a\nBvnDrYoTpsz+HE+ePMOzZ1n4+dg52XN3q+qM11u/jAnTl2vOO/nyTRz8+QwAIG7bAdxIvau5TiqV\nCmXKFHwlFTVeKVUqOyL58k3ZmPXRr/tr2Pn9Ufzw4y9Qq9V49iwL+w8nIvVmmvSxXR2RfLnw+Lw9\nKyPj2o5Cf5RODABgZ2eLiE6vYOb7X+LJk2c4fPQsdsYf1bzo8qvkUAG3zn+DxIPLkXhwOb7bkDfT\n8Ou+T9DkZX/ZfXXp9ArOnr+CLTsP4dmzLMyOXoOGgTVQu6YnAGBAr3ZY/Mlm3LyVhtSbaVj8yWYM\n6v06AKB2TU80DKyB2dFf49mzLGzZeQhnz19B1/C8N9q+3V7DzvgjOHz0LDIzn2LGvFh0DQvRjFEM\n6NkO7y5aiwcPH+N80lV8tma3Zt+tQhqgbNkyWLpiK/76KwtLV2xFmTJl0KZlo0Kvm+JTmZPH9sDi\nd99G9NKNqBrQE1X9e+LtCUsQPestBDfOG2PIP/j4j/yPhbZtjPFvd0GbN/+H2o0H47WWDXW+Y6pU\nKvTv+RoGjVoIt4CeyHqejaXzRwIABvRsCx+vKvCo2xuBrwxHcOMArX0N7ReK35OuwrFaF0QMmIUy\nZcpg57q5uHT5Jrzr94VXvT7YuO2A5jgvxpH/768++R+ynmejTvBbcKoege6D5+LP2/cBAL+c/gPN\nXh+bN8vQ7x0sfX8kfL2rFjiXosYrZVZUfwwcFQ3Hal2waftBqFTy7+j5t+edZ97vnh6u2L5mNuZ9\nuB6Va3eHd/2+WPTxpkKnn6f+txfeXbQWjtW6YPHHm3Qf1Ig+WTgWT59lobJfd/SLnI/li8YhwC/v\n3f7ajTuw9w7XJO/Kro6aHxdnB6hUKlSp7AhrayvZfbk4O2DzlzMx/d3VcKoRgV9O/4FvPp+uiSNy\nUCeEhTZDvZDhqN9iOMJCm2H4oDc027/5fDp+Of0HnGpEYPq7q7H5y5lwdqoIAKjj74Pli8eh7/D3\nUcW/B54+y8InC8dq6s6eMgA1fN3hU78vWneejKixPfF6myAAeWNe29bMwlcb9sCxegS+2rAH29bM\ngpVV4a1jldyXvQwZvRDf7jmGyi6V8NtPqyTLjJ3yMXYnHEd525cQ+/FkNKpfU9cuFdc6fBL692iL\nIf0K9uGJyDCyLYfBfdsjPu79Qrd/t+cYLqWk4uIvX2Llh+MxYuISRQPUl4gLb4jMmWxyaBFcD46V\nKhS6fcfuIxjYK69f0zQoAA8eZeL2nfvKRainkhqoIyqtij1bkXrrHrw8/p2W9HR3wY2bdwssHilJ\n+3YsNNqxiEoLRQYkX2zS812cyPwVu+Xg4eaM63+P8gLAjZtp8HBzKVBO5dSuuIciIgPlpu8pcp1i\nJ4fwDsGIWbUdvbq2xtETv6NSRbtCuxTpuenFPZxRzZ81H1NmTSnRYzipnOQLFUFuuvY1Vjk5FXhc\n6rGiPm5IWSWIcA/9c184qZyEiEeOofeYbHLo/dZ7OPDzGaTdewSvwD6YPWUAnmfnfcwzclAndGzX\nFN/tOY6aLw+EXfmXsDpmkkGBEJFYZJPD+s+myxVBTPQYRYIp7dLTc+HkpNL7d6l6IjP0nMzl/CyN\nop+tsDQhrUKMerz8LwB9fpd6rLAmvNTjRSmrxD4MPSfREoOx7wtTEeY7JEVUWm4CKprScl+w5SAQ\nJboVIg9IslthXthyEIgS3QqRWUq3orRgciAiSUwOAklPzy3S77oeE5Gh52Qu52dpZD+yrdiBnNqZ\nxYIRY1N6EZSlE+keMqdFUIaskGTLgYgkcbZCIJytkL8WZDxsOQiEsxXyZcl4mByISBK7FRZG5OXT\nZF7YciAiSWw5CIQDkvLXgoyHLQeBcEBSviwZD5MDEUlichAIl0/LlyXj4fJpE+Py6aIR6R7i8mki\nKpU4WyEQzlbIXwsyHrYcBMLZCvmyZDxsORAVgyWPGbHlIBDOVsiXFY0S52SMeobgbIWJWfI7T0kQ\n6R4yp+eOsxVmTol3iNz0dK0fqceLUlaJfRT3nNhyME3LgclBIByQlC8rGiXOyRj1DMHkQESSmBwE\nwgFJ+bKiseRuBQckTcycBrVEINI9ZE7PHQckiUgxTA4C4WyF/LUQjSV3K5gcBMLZCvmyouFsBRGV\nOkwOAuFshXxZ0Vhyt4KzFSZmTiPeIhDpHjKn546zFUSkGNmPbMcnnMD46Z9CrVbjrf4dEDWul9b2\ntHsP0S9yPv68k47s7BxMGt0Ng/q0L7GASTeRv+yFzIvOloNarcboqBjEx83D70c+x/rN+3A+6apW\nmZhV29Gofk2cPrgC+3cuxMQZK5CdrS7RoImo5OlMDsdPJqFmNXf4eleFtbUVekW0xvbdR7TKuFV1\nwqOMTADAo4wncHaqCCursiUXsQXjgKR8WdFY8oCkzm5F6q00eHm4av72dHfBsZMXtMoMG9ARbTpP\nhnudnsh4/BQbP/+/YgVUmikxny3y/8rkOgfzWuegMzmoVPI7n7d4PRoG1sD+nYuQfPkm2kVEIfHg\nCtjbly9Qdv6s+ZrfQ1qFIKRViAEhE5Ex6EwOHm7OuJ56V/P39dS78HR30Srz8/FzmD6hDwCgRjV3\nVPOpiqRL1xHUyK/A/qbMmqJEzBYrPZ3fPi13LUSjxDkZo54hdI45BDXyw8WUVFy59ieysp5jw9b9\nCA8N1irjX8sbCQdOAQBu37mPpIs3UN3XzeCASjMun5YvK5pS262wsiqLmAWj0b7bVKjVORjaLxQB\nfj5YEbsLABA5qBOmTeiNwaMXokGLSOTk5CB69jA4OVYsVlBEZHpcIWli+VfZsVshfU5avwt0Dzmp\nnMymW2HICkkmBxMzpyW4IhDpHjKn547Lp4lIMUwOAuGXvchfC9FY8iIoJgeBcLZCvqxoLHm2gsmB\niCTxH+laGJGXT5N5YcuBiCQxOQiEA5Ly10I0HJAko+CApHxZ0XBAkohKHSYHgfDLXuTLisaSuxVc\nPm1i5rQEVwQi3UPm9Nxx+TQRKYbJQSCcrZC/FqKx5G4Fk4NAOFshX1Y0nK0golKHyUEgnK2QLysa\nS+5WcLbCxMxpxFsEIt1D5vTccbaCiBTDT2VaGJG/Q5LMi1GTgzk1w0hMvIeMx6jJQaT+oiiU/vZp\nkSnyjcsC3UPm9O3ThuCApInxnbBoRLqHzOm544AkESmGyUEgXD4tfy1EY8nrHJgcBMLl0/JlRcPl\n00RU6jA5CITLp+XLisaSuxWcrTAxcxrxFoFI95A5PXecrSAixTA5CISzFfLXQjSW3K1gchAIZyvk\ny4qGsxVEVOrwU5kWhv8rk5TClgMRSZJNDvEJJ+DfdAhqBQ3EgiXfSJbZfzgRjV59G4HNh6FV2ETF\ngywtOCApfy1EY8kDkjq7FWq1GqOjYpCwdQE83FzQ+LVRCA8NRoCfj6bMg4ePMWryMny/6X14ergi\n7d7DYgVUmnFAUr6saErtgOTxk0moWc0dvt5VYW1thV4RrbF99xGtMus2/YiuYSHw9HAFALg4OxQr\nICISg87kkHorDV5/v+gBwNPdBam30rTKXExORfqDDLQOn4SgNiPx9Yair8SiPFw+LV9WNKW2W6FS\nyTdLnmdn49fES9i7LRpPnv6F4PZj0SwoALVqeBYrsNJIieajyLMV7FaYV7dCZ3LwcHPG9dS7mr+v\np96Fp7uLVhkvD1e4ODnA1tYGtrY2aBlcD4lnUySTw/xZ8zW/h7QKQUirkGIFT0QlR2e3IqiRHy6m\npOLKtT+RlfUcG7buR3hosFaZzh2a4/Cxs1Cr1Xjy5BmOnbyAOvkGLPObMmuK5oeJoSDOVshfC9GU\n2m6FlVVZxCwYjfbdpkKtzsHQfqEI8PPBithdAIDIQZ3gX9sboW0ao35IJMqUUWHYgI6o4y+dHEg3\nzlbIlxWNJXcr+JFtE1P6Y7+W/n8rRLqHLP0j21w+bWFEHpAk88Ll00QkiclBIByQlL8WorHkAUkm\nB4FwQFK+rGgseUCSyYGIJDE5CITLp+XLisaSuxWcyjQxc5oOE4FI95A5PXf89mkiUgzXOQgkPb34\n/2pd5EVQhp6T1HZRiNSSKYyhLRy2HATC2Qr5smQ8TA5EJInJQSCcrZAvS8bD2QoTM6cRbxGIdA85\nqZyEiqcwTionzlYQkXI4W2FhRJ6tIPPClgMRSWJyEAgHJOXLkvFwQNLEOCBZNCLdQ5Y+IGnUMQe+\nEKi4RLuHRItHSUZNDuaQZY0t/83F5dPS56T1u0D3kDm1HAzBMQeBcPm0fFkyHiYHIpLE5CAQzlbI\nlyXj4WyFiVnygFZJEOkeMqcxBy6fJiLFcPm0QDhbIX8tyHjYchAIZyvky5LxMDkQkSR2KywM/1cm\nKYUtByKSxJaDQDggKX8tyHjYchAIByTly5LxMDkQkSQmB4Fw+bR8WTIeLp82MS6fLhqR7iEunyai\nUkk2OcQnnIB/0yGoFTQQC5Z8U2i5E78mwcq1PbbsPKRogKWJEt2K3PR0rR+px4tSVol9FPec2K0w\nDZ3JQa1WY3RUDOLj5uH3I59j/eZ9OJ90VbJc1OzPEPpaY+TyeTQYZyvky5Lx6EwOx08moWY1d/h6\nV4W1tRV6RbTG9t1HCpRbtnI7uoW3gKuzQ4kFSkTGpTM5pN5Kg5eHq+ZvT3cXpN5K0y5zMw3bd/+M\nEUPCAAAqJnmDcbZCviwZj84Vkio9Xunjp32C+TOHQqVSIRe57FYUgxLdCpE/W8FuhXnRmRw83Jxx\nPfWu5u/rqXfh6e6iVeZk4kX0emseACAt/SF2J5yAtXVZhHdoXmB/82fN1/we0ioEIa1CihU8ERV0\neP9hHN5/uNj70bnOITtbDb8mg7F3WzTcqzqjSdvRWL9qGgL8fCTLDx71AcJCmyEirEXBA3GdgySl\n1zmI/NkKJYh0D1n6OgedLQcrq7KIWTAa7btNhVqdg6H9QhHg54MVsbsAAJGDOhU5SKLiEO0eEi0e\nJcmuc+jQrgmSjq/GpZNfYup/ewPISwpSiWH1x5MlWw3/MHQQytCy5lJPyWOITJTrbc7PrzHvC6Mu\nnyYi0+DyaSJSjFGTg6U00URudlr68mlzrCdCbIZgt8LCWPpsBRmG3QoiUgy7FQLUU/IYIhPlepvz\n88tuBREpit0KIlIMv5rewnBAkpTClgMRSeKApAD1lDyGyES53ub8/HJAkogUxQFJIlIMuxUC1FPy\nGFw+LV49EWIzBLsVFoazFSSF3QoiUgy7FQLUU/IYIhPlepvz88tuBREpit0KIlIMuxUC1FPyGJyt\nEK+eCLEZgt0KC8PZCpLCbgURKYafyrQwIv87PDIvbDkQkSQOSApQT8ljcEBSvHoixGYIDkhaGA5I\nkhQOSBKRYtitEKCekscQmSjX25yfX3YriEhR7FYQkWLYrRCgnpLH4GyFePVEiM0Q7FZYGM5WkBR2\nK4hIMexWCFBPyWOITJTrbc7PL7sVRKQodiuISDF6fSozPuEExk//FGq1Gm/174Cocb20tq+N24vo\npRuRm5sL+wq2+HTRONSvW71EAibdOCBJSpFNDmq1GqOjYpCwdQE83FzQ+LVRCA8NRoCfj6ZMdR83\nHPx2MRwq2iE+4QSGj/8QR/csK9HAiahkyXYrjp9MQs1q7vD1rgprayv0imiN7buPaJUJblIHDhXt\nAABNg/xx42aa5L4sZXBH5AErkYlyvc35+RVqQHLT9oP4/sdfsGrJBADAmo0JOHbyApYtGC1ZfuGy\nOPyRfAMrP/qv9oE4IElkMoYMSMp2K1Qqld4723foNL5YG4+f4j8qciBEJBbZboWHmzOup97V/H09\n9S483V0KlDtzLgXDxn+IHevmwLGSveS+LKWJJnKzk8unxasnQmyGkO1WZGer4ddkMPZui4Z7VWc0\naTsa61dN0xqQvHbjDtp0now1y6PQrHEd6QOxW2EUnK0gKSXSrbCyKouYBaPRvttUqNU5GNovFAF+\nPlgRuwsAEDmoE+ZEf437DzIwYtJSAIC1tRWOJ8QUORgiEodei6A6tGuCpOOrcenkl5j6394A8pJC\n5KBOAIDPlk7EveQtOHVgOU4dWF5oYrCUJprIzU6RiXK9zfn5FapboRR2K4hMh8uniUgx/FSmAPWU\nPAZnK8SrJ0JshmC3wsJwtoKksFtBRIrh/8q0MPxfmaQUthyISBIHJAWop+QxOCApXj0RYjMEByQt\nDAckSQoHJIlIMexWCFBPyWOITJTrbc7PL7sVRKQodiuISDHsVghQT8ljcLZCvHoixGYIdissDGcr\nSAq7FUSkGC6ftjBcPk1KYcuBiCRxQFKAekoegwOS4tUTITZDcEDSwnBAkqRwQJKIFMNuhQD1lDyG\nyES53ub8/LJbQUSKYreCiBTDboUA9ZQ8BmcrxKsnQmyGYLfCwnC2gqSwW0FEimG3QoB6Sh5DZKJc\nb3N+fi22W5Gemy5fsBRwUrGZbQjeP4ZxUjkZ1K0w6gev+KKg4uD9Y1xGTQ7M/HlK8ia39AFJ3kNF\nZ+j9xgFJIpLE5GAiIg08GYsSg3BkPByQNAH2nQ3He6joDB2QZMuBiCTJJof4hBPwbzoEtYIGYsGS\nbyTLjJ3yMWoFDUSDFpE4deaS4kFaIi6fZrdCdDqTg1qtxuioGMTHzcPvRz7H+s37cD7pqlaZ7/Yc\nw6WUVFz85Uus/HA8RkxcUqIBG9Ph/YdLbN9OTiqdj+nzu67HRGToOYl2fiV5X4hEZ3I4fjIJNau5\nw9e7KqytrdArojW27z6iVWbH7iMY2Ot1AEDToAA8eJSJ23ful1zERlRabgIqmtJyX+hMDqm30uDl\n4ar529PdBam30l4oc69AmRs37yocpuXhbAW7FaLTuQhKpdKvOZebq/3k6VuvNCupboXIX01vKd2K\n0kJncvBwc8b11H9bAddT78LT3UVnmRs30+Dhpl3mH+Y4hRc9O9rUIVA+otxDpeG+0Jkcghr54WJK\nKq5c+xPuVZ2xYet+rF81TatMeIdgxKzajl5dW+Poid9RqaIdqlR2LLAvQ+ZZich0dCYHK6uyiFkw\nGu27TYVanYOh/UIR4OeDFbG7AACRgzqhY7um+G7PcdR8eSDsyr+E1TGTjBI4EZUso62QJCLzovgK\nSXNbNCUX79q4vWjQIhL1Q4bjldBxOHMuxQRRatPnGgPAiV+TYOXaHlt2HjJidAXpE+/+w4lo9Orb\nCGw+DK3CJho5woLkYk679xCh3aaiYctIBDYfhth135sgyn8NGb0QVfy6o94rwwotU9TXnaLJwdwW\nTekTb3UfNxz8djHOHF6JGZP6Yfj4D00UbR59Yv6nXNTszxD6WmPkmrBtqE+8Dx4+xqjJy7Bz3Vyc\n/XkVNsVu/9IEAAAC70lEQVTONFG0efSJOWbVdjSqXxOnD67A/p0LMXHGCmRnq00UMTC4b3vEx71f\n6HZDXneKJgdzWzSlT7zBTerAoaIdAKBpkD9u3EyT2pXR6BMzACxbuR3dwlvA1dnBBFH+S5941236\nEV3DQuD593oZFzOI2a2qEx5lZAIAHmU8gbNTRVhZlTVFuACAFsH14FipQqHbDXndKZoczG3RlD7x\n5vf51/Ho2K6JMUIrlF7X+GYatu/+GSOGhAEATLnsRJ94LyanIv1BBlqHT0JQm5H4eoNpZ7b0iXnY\ngI44d+Eq3Ov0RIOWkVgyb6SxwywSQ153in4TlLktmirKcfcdOo0v1sbjp/iPSjAiefrEPH7aJ5g/\ncyhUKhVykWvSboU+8T7PzsaviZewd1s0njz9C8Htx6JZUABq1fA0QoQF6RPzvMXr0TCwBvbvXITk\nyzfRLiIKiQdXwN6+vBEiNExRX3eKJgelF02VNH3iBYAz51IwbPyHiI+bB8dK9sYMsQB9Yj6ZeBG9\n3poHAEhLf4jdCSdgbV0W4R2aGzVWQL94vTxc4eLkAFtbG9ja2qBlcD0knk0xWXLQJ+afj5/D9Al9\nAAA1qrmjmk9VJF26jqBGfkaNVV+GvO4U7VbkXzSVlfUcG7buR3hosFaZ8A7B+OrvZqOuRVPGoE+8\n127cQcSA2VizPAo1q3uYJM789Ik55dTXuHw676dbWAt8unCsSRKDvvF27tAch4+dhVqtxpMnz3Ds\n5AXU8fMxSbyAfjH71/JGwoFTAIDbd+4j6eINVPd1M0W4ejHkdadoy8HcFk3pE++c6K9x/0EGRkxa\nCgCwtrbC8YQYoWMWiT7x+tf2RmibxqgfEokyZVQYNqAj6vibLjnoE/O0Cb0xePRCNGgRiZycHETP\nHgYnx4omi7n3W+/hwM9nkHbvEbwC+2D2lAF4np2tideQ1x0XQRGRJH5NHBFJYnIgIklMDkQkicmB\niCQxORCRJCYHIpLE5EBEkpgciEjS/wM6mX4iEn445gAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plotdata.plotframe(2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Adjust some plot parameters:" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "\n", "Current plot figures, axes, and items:\n", "---------------------------------------\n", " figname = pcolor, figno = 0\n", " axesname = AXES1, axescmd = subplot(1,1,1)\n", " itemname = ITEM1, plot_type = 2d_pcolor\n", " \n", " figname = contour, figno = 1\n", " axesname = AXES1, axescmd = subplot(1,1,1)\n", " itemname = ITEM1, plot_type = 2d_contour\n", " \n", " figname = cells, figno = 2\n", " axesname = AXES1, axescmd = subplot(1,1,1)\n", " itemname = ITEM1, plot_type = 2d_patch\n", " \n", " figname = q, figno = 300\n", " axesname = AXES1, axescmd = subplot(1,1,1)\n", " itemname = ITEM1, plot_type = 1d_plot\n", " \n" ] } ], "source": [ "plotdata.showitems() # shows the currently defined figures, axes, items" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Change the color limits for the pcolor plot to show the overshoots and undershoots with Lax-Wendroff. Also suppress showing the other plots..." ] }, { "cell_type": "code", "execution_count": 19, "metadata": { "collapsed": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVMAAAEKCAYAAACvyrZLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8E2X+wPFPSgu9KL04Wu77kEPWIhZBuQoFpUjBFeRS\nFFEB4SdqEVYX1pUVFkWl7Ip4gIjAIiKiUKVKReRGBLmPcpazlFLaUkvT/v5IkyZp0kyOdtLwfb9e\nebXJPDPznUzmm+eZeeaJpihjYxFCCCGc4qV2AEII4QkkmQohhAtIMhVCCBeQZCqEEC4gyVQIIVxA\nkqkQQrhApU+mp89ewiusD4WFhQ7Nv2zVj/QdPNXFUZWPyhSrEHcat0mmW7YfoEvfSQQ3eoSwpvF0\n7TeZ3XuPunQdlhLv8Ed78f3qt1y6HiW6D5jCx0s3WJ2udqxeYX1IPX2xQtZlj9NnL9Ej7iUC6g2g\ndecx/Pjzbzbnyc+/TevOY6jf9nG7lvXFlz/RsP1wAusPYNDIGVzPvGmY9uef+YyZMJcaDQcS0fox\n5v1ntcm8v/9xgnt6PE9AvQFE9XyefQdOmkyf95/VRLR+jBoNB/LUxLfJz79tmJZxPYtBI2cQWH8A\njTqMYPnqn0zm/fHn32jVeQwB9QbQc+DLnD1/xeZ7IMqfWyTTrKwcHh76NyaNG8T1U2tIO7iCv78y\nkmrVqpbL+oqK1L9PQaPRKCqnZqzu8D6ZG/b0LO7p0JyMk6t5829PMuSJN0i/dqPMef49/3/UqhmM\n+Vte1rIOHj7Nsy++x7IPX+Xykf/h71eN51963zDvjNlLOXn6Amf/+IJNa//NnPn/4/sfdwG65D1w\n+N8Z9VhvMk+tYfTQPgwc/ndu3y4A4PsfdzH7/ZX89PUczuxfRuqZi/z9rc8Myx7/8nx8q/lw5egq\nli2cynNT3ufQkTMApF+7weDR/+DN6U9yPfUrou5uwWNj/un0+yqc5xbJ9NjJ82g0Gh6L745Go8HX\ntyoxPe6hXZvGgO6g/ufcZTTqMILaLR9l9PNzyMrKsbisRh1GmNQwZrz1GSOf1dXmHnjoRQCCGw8i\nqMFAtu86xOIvvqdb//8zlN+64yCdeo0nuNEj3Nt7Att2HjJM6z5gCq/PWkzXfpMJajCQvoOnci0j\ny2IcmTeyeXjo36jV4lFCm8QzYNhrpF1IB2D6Pz/hl21/MCEhkeoN4nhh6oJS8yuJ1SusD//9ZB3N\no0YT1GAgr89azMlTF4ju8wLBjR5h6FNvGg5ggG+/387dD4wjpPEg7o+dxB+HTlmMXb/uDg+Mo3qD\nOFZ9/bPFchXt2Inz7P3jBDOnjqJatarED+hG+7sas3rdL1bnOXXmIstW/cSrk4dh/N1ga1nLvvyR\nuH7RdL2vLQEBfrwx7Qm++vZXcnJuAfDZyo289tIIagQF0KpFA54Z1Z/Fy38AIGXLPrTaQiY9G4+P\njzcTn3mEoqIiNv3yOwBLVmzk6ZH9aN2yIcE1Ann95RGGeXNybvHVt7/yxrQn8Pf35f772jKwfxeW\n/i8ZgK/WbaFt60YMjutG1ao+zEgYyb6DqRw7cd7l77ewj1sk05bN6lOlihdPjJ9DUvIuk+YUwKfL\nvmfJih9IWTeX1N8+Izv7FhMSEi0uS6PRmNT6jP//Zf08AG6c/pqss2u5r1Mbk3kzrmfx0NC/MfnZ\neDJSv+LF5wbz0NC/mcSz/KtNLF7wMleOrSL/dgFzE1dZjKOwsJCnRsRydv8yzu5fhp9vVUPMb/5t\nDN2i27FgzkRunv2G998aX2p+W7Hq/bBpD3tTPmD7D+8z+/3/MXbyPJZ/NI2z+5fxx6FTLF+9CYC9\n+0/w1Atvs+jdF8lI/YpxTzxM3OOvmzQv9TZ/9w4A+3/5kJtnv+HRRx4sVWbL9gOENB5k9bF1x0GL\n8Trj4JHTNGkYQUCAn+G1Dm2bcrC41mbJxIQF/Ov1p/D1NW3l2FrWwSNn6HBXE8O0Jo0iqFbVh2Mn\n07ieeZOLlzLo0LZkevu7mpjM2/6uxibr69C2ZPqho6bLbn9XEy5fuc71zJscO5mGt3cVmjWpWzLv\nXU04eOS0IW7j9fr7+9KscSQHDlv+YhQVxy2SafXq/mxZPw8NGsZOnketFo8ycPjrXLl6HdDVEqaM\nH0KjBnUICPDjX68/xYqvUhRddDJuqtpqtn73ww5aNqvH8Ed74eXlxdDBPWjVvD7fbNgG6BLzk4/3\npVmTuvj6VuWvjzzI73+ctLis0JAgBj3cFV/fqgQG+jHtxWH8/Ot+q7GVFXdZXpn4VwID/WjTqiHt\n2jSiX69ONGpQh6CgAPr17sTe/ScA+HDJd4wb/RCd/tISjUbDqKExVKvmw/bdhxWtx1zX+9py/dQa\nq48une9yaLllyc65RY2gAJPXgqr7c/NmrsXya77dQlFREQP7d1G0rOqBfmRn62qeObl5lteVnWso\nYzxdN+1WGXEGcDM7t3h6Xql5AW5m3yI755bheUlcxsvOI6h66biyc/Isvgei4nirHYBeqxYN+HTB\nywAcPX6OEc++xeRp/+WLRdO4eCmDhvVqG8o2qFeLggItl69cd2kMFy5do0G9WiavNaxfiwuXrhme\n16kVavjfz7cq2cXNPnO5uXn83/QP+P6n3YaabXZOHkVFRYbastLzpmWpXSvEKJ5qJs99faty5Wom\nAGfOXeazlcnMX7TWMP12QQEXL2U4HUNFCQzwI8sscWbeyC6VfEDXXH5lxiI2/G+W4mXdyMqheqCf\nYfoNs1NJN27mUD3Qn8DiMlk3cwkPq2FxXktxVg/0L57uazJdv57qgX5W49JvY/VAP7Ju5pSarl+3\nUI9b1EzNtWxen9FDYzhw+DQAkRFhnD532TD97PkreHtXMUkcegH+vuTklnxLX7pSkixsJa+6EeGc\nOWd6ZfTMuSvUjQi3exveXvAlx06eZ2fyfG6cWcvP375NUVGRocZpKxZXJFoNJctoUK8W018cZlJ7\nzD63jsfiuzu07F+2/UH1BnFWH79uP+B0/ObuatWI1DMXDTVDgH0HUrmrVaNSZY+npnHm3BW6PfQi\nEa0fY/DomVy8fI2I1o9x9vwVm8u6q1VD9h1MNUw7eeoC+fkFtGhal5Dg6kTUCTVplew7kErb1iXz\n7j9o2uz+49Ap7mrV0LAdvx8wnvcktWuFEBJcnRZN61JQoOVEalrJ9IMnTebdd6AkrpycW5w8fdHi\neyAqllsk06PHz/HOgi8NF2jOnb/C8tWbiC4+Tzgsvgfz/rua02cvkZ19i2lvfMLQ+O54eZUO/+52\nTVnxVQoFBVp27z3K6nVbDImpZlgNvLw0nDx1wWIc/Xrfy7GT51m++icKCrSs/CqFI8fP8XDfzoYy\nSpvf2Tm38POtSo2gADKuZzFzzucm02vXDObkactxKInVGpPTGpQk77Gj+vPBp9+yc88RioqKyMm5\nxXc/7DBJJibx1Qopc93dottx8+w3Vh/339fWrriVaNGsHne3bcrMOUvJy8vnq3W/cODwaQbHdStV\ntl2bxpw/8AX7Nn/Avs0f8NF7L1K7Vgj7Nn9Avchwm8saPqQX65K2sWX7AXJybvHarMUMHtDVcI51\n1GMx/PPtZWTeyObw0TN89PkGnhjWB4DuXTtQpYoX7y9cw59/5vP+wjV4eXnR84GOxfP25uPPkzh8\n9AzXM2/yxtxlPPm4bt6AAD/iH76f1/+1hNzcPLZsP8C6pO2M/GtvAAY9fD8HDp/mq3W/kJeXz8w5\nn3N326a0aFbP5e+3sI9bJNPqgX7s2HOEzjETCaw/gOi+k2h/VxPefmMcAGNGxDLyr7154KEXafKX\nUfj7V2P+7AmG+Y1rcW9MG83JUxcIaTKIGbOXMnxIT8M0f39fpr/4OPf3+z9Cm8SzY/dhkwtWYaFB\nfLv8Dd5e8CXhzQYzd8Eqvl3+BqEhQRbXZX6xy9jkZ+O5lZdPePMhdImdTL9enUzKTho3iC+/+YXQ\nJvFMfvU/pea3Fat5LBbjo6T8PXe3YNG7/8eEhERCm8TTvNMTfLZio8XYAWYkjGT0+DmENB7El2s3\nWy1X0VZ8PJ3dvx8jtGk80//5KauXvE5YqG7/6GvLAFWqVKFWzRDDIyS4OlW8vKhVM8TwJVzWstq0\nasgH70xi+DP/onarv3IrL5//zH3BEMfMqaNo2iiShu2H02PgyyS88Bh9ekYBULWqD19/PoPPVm4k\npEk8n63cyNefz8DbuwoAfXt14pWJf6XHwJdp1GEETRtFMnPqaMOy/zP3BW7l5VOr5aOMGPcWH7w9\nidYtdTXT8LAarF7yOtP/+SmhTePZ/fsxVnw8vVzfc6GMxtbg0GMmzOW7jTuoFR7MH78usljmhakL\n2JC8E38/XxYveJmO7ZuVS7BCCKFnKzctW/Ujc97/H0VFRVQP9OO/b+sqaQBJybuYPP2/aLVanh7Z\nj4RJQ52Ox2bN9MnhfUla9S+r09dv3MGJ1DSO717Ch/Mm89yU95wOSgghbLGVm5o0jGDzd++wf8uH\nvPbSCJ6ZrOtuqNVqmZCQSNKqWRza9jHLV2/i8FHr3euUsplMu0W3IyQ40Or0bzZsY/RQ3fmezlGt\nyczKcflVdiGEMGcrN0Xf28bQBa1zVCvOF1+T2bnnKM0aR9KoQR18fLwZGt+DtcXdH53h9DnTtIvX\nqF+3puF5vchwzl+46uxihRDCZT5emkT/mHsBSLuYXipnpV1Md3odLulnan6F2xXdeoQQwhU2/fI7\nnyxL4tekd4Hyy09OJ9O6EWGcSyupiZ6/kG6xX6YmNMbZVQkhHFSUYb3nhprszQv2bsf+g6mMnTyP\npFWzCAmuDpTOWefSrlIv0v6+5OacTqZx/aJJXLSWoYN7sH3XIYKDAix2pgf4hCedXZ1LeVFo+OtL\nHuHoqvpNSKVhzTPMCIcZbYD7AH0HBV8gC9DfC5BttEBvVL2nbMZ6mNFfvfU7wqmYC4sfALno9smJ\n4ue/walzjdmK7lbSVJqQhy8AgWQTTjr+6O408qKQguIdV4A3hXgZPhuWfM1eHqGjg0FXvDF8qnYI\nZfJCq6hcIVXsWu7Z81eIHzWTzz9IMBnrIKpjS46npnH67CUi64Sxck0KyxdNs2vZltg89Ic9/SY/\nb91P+rUs6rd9nJlTR3G7QDcS0bgnHqZ/TGfWb9xJs3tGE+Dvy6eJLzkdlBDizlHWF5cx81K2ctM/\n5izleuZNniseOtHHx5udyYl4e1chcfYE+g55Fa1WNyCRvh+vM2wm0+Uf2e4QnDhnoqKVfcvDispV\ntKrk408ukeju+EknnKCrWeB/HdLR7UXf4sL+QD4ll+4sfQ4sTSvETW6R8ADm7ytAAbraafF1hFvn\n/NhKF9ajq/bupz256O5vDyaTSC4QhG74RC8KyUc3qpT+r56lA/0YOW77WbbM3Wumjv1Khq3c9NH7\nU/jo/SkWp/WLuZd+xRekXKVCG6W/c3dFrs4m/U70poBAsskkGAB/cmnDIbpXu65rxheC4RjzpfS7\nVmj2VyXdm6u7fkc4HbPxe58HxfmRS9QhlSbspz0ABzC9vfUSdQilZNwG42Rq3My3dKDfIs/tPsuV\nmTcFtguhq8O4M6krlaF7JRuI545Mpjboa6Ou5Md9Ll+mPbKx3reyMvKiUNHD3VVozdSdPgTGO6cq\n+XhRaDjw8vA1XJAwfGl6Gf2VryD3Vbxb86lquOBkKaHm4k/V4rqOeTPf/AKU/v9CN9rx7nQsOasy\nJEolKjSZGjer3IFxMz+ILILRjf0ZRBa+5GFy+sz43Jxn7HvP44XhE+1PLoFkG/apsarkE0qGYZqt\nZFrZufvP7XnKe12hyVR/gcddmJ8zrVX8sQslg0CyS5KpFyU1VOPuOEJ9xi0GbwwXCvXdn/SfuUvU\nMdRQQ8lw6gJUZXNE7QBs8IT3GNxopH0hxJ1JkqkDLDW53IG+ZhpY3APfn1xdM9/Su2OpZupl9leo\nw6iZ70se/uQaap+hZBjOkQaTaXJap6xmvqcc6O7MU97jO7pmanzAGF8xrCxXD0XZzPejtf1tPo/x\nX1H+lHaNcnd3dDK1STraV3r5VDXUSD2Fp22Tp3xxVWgyNT+5rzbjLi/5VDXEV4B3Sdcoe/azZ3wm\nKi+jUzCFeBn2oX7fGjfljZ/rXwNl9+a7A3c7lpzh7u+1UhWaTAvctCKsT6bmg11YZT7JMz4LnqF4\nX+i/EM0TqPH/+ufGA53IOdOK5ynvcYVmN3fq9AzWd6L+jJqH7OM7lvHnzbBPjV43n25tmihfkkyF\nEMIFJJk6oLK8aYYrvVI5qdTMr+SbN90t9d7wotBmM19fa60sn2d35ynv4x1dM1XSVcYiz9j3nslL\n/8dy4iwp5njXKE85+N2FdI0SQggXcPTLacyEuXy3cQe1woP549dFFsu8MHUBG5J34u/ny+IFL9Ox\nve4nMxp1GEFQdX+qVPEyDBrtLEmm4o4hNUr35Oh+eXJ4XyY+8wijnpttcfr6jTs4kZrG8d1L2LH7\nMM9NeY/tG+cDuh/VS1k3l9CQIIfjNidnBcUdw/iKvnAfjo5n2i26HSHB1oci/GbDNkYP7QNA56jW\nZGblcPnKdcN0sx9VdsF2CCGEisprcOi0i9eoX7em4Xm9yHDSLup+10ajgd6DXiGq5/MsWrLeJdsh\nzXwhhKrK8/RLkZXq55b184iMCOdqeiYx8Qm0alGfbtHtnFqXJFMhLJDuTxXH2nucxR6y+M3h5daN\nCONc2lXD8/MX0qkbEQ5AZPHfmuHBDHqoKzv3HHE6mUozXwgLZOSwiuNtuPnX9BFKBxrxpOFhr7h+\n0Xy2ciMA23cdIjgogNq1QsjNzePmzVwAcnJu8cOm3bRr09gF2yGEECpy9Etr2NNv8vPW/aRfy6J+\n28eZOXUUtwt0fVbHPfEw/WM6s37jTprdM5oAf18+TXwJgEtXrhM/agYABQWFDH+0J316Rjm9HZJM\nhRCqcjSZLv9ous0yiXMmlnqtSaMIft+80KF1lkWSqRBCVZ5yOkWSqRBCVZJMhRDCBSSZCiGEC8hA\nJ0II4QJSMxVCCBeQZCqEEC4gyVQIIVxAkqkQQriAJFMhhHABSaZCCOEC0jVKCCFcwFNqpjaH4EtK\n3kWrzmNoHjWa2e+tKDU9/doNYoe8yt0PjKNtl7Es/uL7cglUCOGZymuk/YpWZjLVarVMSEgkadUs\nDm37mOWrN3H46BmTMomL1tKxfTN+37yQlHVzmfLaQgoKtOUatBDCc9wRyXTnnqM0axxJowZ18PHx\nZmh8D9Zu2GZSJqJOKFk3cwDIuplLWGgQ3t5Vyi9iIYRHuSOSadrFdKs/SKU3dlR/Dh45Q2Sbx+jw\nwDjem/V8+UQqhPBIziRTW6chr2feZNDIGXToNo7OvSdy8PBpxfPaq8wLUBqNxuYCZr2znLvbNiVl\n3ducPHWBmPgE9m1eSPXq/qXKHudzw/+htCeM9g6ELIQoyzX2k8F+tcNQzNFap/40ZPKa2dSNCKdT\nr/HExUbTumVDQ5lZ7yznL+2bsWbpDI4eP8f4V+aTvGaOonntVWYyNf9BqnNpV6kXGW5SZuvOg0x/\n8XEAmjaOpHHDOhw9cY6oji1LLa85IxwOVAihTJhZReUky1SMxjZHu0YZn4YEDKchjRPi4WNnmTp5\nKAAtm9fn9NnLXLl6nZOnLtqc115lNvOjOrbkeGoap89eIj//NivXpBAXG21SplXzBiT/vBeAy1eu\nc/T4eZo0inA4ICHEncXRZr6S05Ad2jbhq3VbANi55whnzl3m/IV0RfPaq8yaqbd3FRJnT6DvkFfR\nagt5akQsrVs2ZOHibwHdj1ZNe3EYT06YS4du4ygsLGTOzLGEhgQ5FZQQ4s5hrZl/nuOkcdzqfEpO\nQ06dNJRJr/6Hjg8+S7s2jenYvhlVqngpmtdeNjvt94u5l34x95q8Nu6Jhw3/h4fVYN3yN1wemBDi\nzmAtmTagKQ1oani+iw0m05Wchqxe3Z9Pin+VFKDx3SNp2iiSW7f+tDmv/dshhBAqcrSZr+Q05I2s\nHPLzbwOwaMl6Hry/HYGBformtZfcTiqEUJWjV/OVnIY8dOQMT0z4NxqNhratGvHx/CllzusMSaZC\nCFU50yHf1mnI6HvbcHTnp4rndYYkUyGEqmTUKCGEcIHKcKuoEpJMhRCqkmQqhBAuIMlUCCFcQJKp\nEEK4gCRTIYRwAUmmQgjhAtI1SgghXEBqpkII4QKSTIUQwgUkmQohhAtIMhVCCBeQZCqEEC4gV/OF\nEMIFpGYqhBAu4EwyTUrexeTp/0Wr1fL0yH4kTBpqMn3u/FUs+/JHAAoKtBw+dpb0E6sJrhFIow4j\nCKruT5UqXvj4eLMzOdGp7ZBkKoRQlaPJVKvVMiEhkeQ1s6kbEU6nXuOJi402GTH/pYmP8tLERwH4\n9vvtvPvBVwTXCAR0P8iXsm6uy34AVH4DSgihKkd/A2rnnqM0axxJowZ18PHxZmh8D9Zu2GZ1PV98\n+RPD4nuYvFZU5MrtEEIIFTmaTNMuplO/bk3D83qR4aRdTLe4jtzcPL7/aTeD47oZXtNooPegV4jq\n+TyLlqx3ejukmS+EUJW1Zv5hLnGYy1bn02g0itexLmk7Xe9ra2jiA/y64V0i6oRxNT2TmPgEWrWo\nT7fodsoDNyPJVAihKmtdo9oRTjtKfsv+a/abTK8bEca5tKuG5+fSrlIvMhxLVqxJKdXEj6gTBkDN\n8GAGPdSVnXuOOJVMpZkvhFCVo838qI4tOZ6axumzl8jPv83KNSnExUaXKncjK4fNW/czsH8Xw2u5\nuXncvJkLQE7OLX7YtJt2bRo7tR1SMxVCqMrRq/ne3lVInD2BvkNeRast5KkRsbRu2ZCFi78FSn7y\n+evvfqVvzyj8/KoZ5r189TqDRs4AoKCgkOGP9qRPzyintkOSqRBCVc70M+0Xcy/9Yu41eU2fRPVG\nD+vD6GF9TF5r3DCC3zcvdHi9lkgyFUKoSu6AEkIIF5BkKoQQLiDJVAghXEBGjRJCCBeQmqkQQriA\nJFMhhHABSaZCCOECnpJMbd5OmpS8i1adx9A8ajSz31thsUzKln10fPBZ2nYZS/cBU1wepBDCczl6\nO6m7KbNmqmTw1cwb2Yx/eT7ff/kv6tWtSfq1G+UetBDCc1SGRKlEmTVTJYOvfvHlTwwe0JV6xeMK\nhofVKL9ohRAex5sCRQ93V2YyVTL46vGTaWRk3qRH3EtE9XyepSs3lk+kQgiPdEc085UMvnq7oIDf\n9p3gx6/nkHvrT6L7vsB9Ua1p3rSey4IUQniuypAolSgzmSoZfLV+3ZqEh9bAz68afn7VeCC6HfsO\npFpMpsf53PB/KO0Jo72z8QshzFxjPxlmAym7szsimRoPvhpZJ4yVa1JYvmiaSZmB/bowISERrVbL\nn3/eZseeI7z4/BCLy2vOCNdFLoSwKMysonKSZSpGY5unJNMyz5kaD77aJvppHhvU3TD4qn4A1lYt\nGhDbsxPtu46jc8xExo7qT5tWDctarBBCGDhzztSZrptK5rWHzU77SgZfNf5taiGEsIejNVNnum4q\nmdf+7RBCCBU52jXKma6bSua1lyRTIYSqHG3mO9N1U8m89pJ784UQqrLWzN9BLju4ZXU+Z7puKpnX\nXpJMhRCqspZMo/ElGl/D80QyTKY703WzXmS4zXnt3w4hhFCRo818466b+fm3WbkmhbjYaJMyA/t1\nYcuOA2i1WnJz89ix5whtWjZUNK+9pGYqhFCVo1fzjbtuarWFPDUi1tB1E3S9joy7bnp5aUy6blqa\n1xmSTIUQqnJmEBNnum5amtcZkkyFEKrylDugJJkKIVQlyVQIIVxAkqkQlYynHLSexlP2iyRTISoZ\nLwop9KBejZJMHeBub5o+HvO43C1O4Rjj/Wj+v63nlv53J+4alyM8ZVsqNJm62++46Hei+WAKXhTa\nF6uX2V+hOn2CrEo+oNvHxv+b7299Tc+8xucpB7o7c7e84CjVD/9Ast1mmj+5pi8Yf9XIMVU5mFUP\njPep+b53p8+ePdOMn9szzV3dEb8B5WrWdq7aH05vCgyv+ZOLL3mGWkyphOr++/TOpa8aeENV8vEl\nD9Dt00CyycWfQLIJIsuwvwPJxotC8qkKYPhrPM0StT+znpRQK0OiVKJCk2kkFypydYp5U4A/udTi\nCgDhpOs+hPrjygsMLZGC4keh0TTP+CxUbl7oPs3F+8yn2m3C/0w3fOYyCSYXfwCCyaQWVwyJxjiZ\nFpgdEp5yoLszT3mPKzSZtnezH/kyvgDlSx7h6MYzbEIqmtpFEGRUWJ9M84z+L1mQKUufDdVPqHgQ\n/XtZiOn56qqU7LNa0ORcKunoRgLyJ5e84hGIgsgilAzDKQAvCg1JtABvCvHymAO8MvCU91oO8bKE\nqh2AcFZvktUOweX68IPaIbiUnDN1gLt+CPRXevXNPk3tImhWPDEU3VeO/jpGAZCPaTMfC/+LiuUN\n+FLyBdgEGuad4frVEHqTTBsOGWqfvuQRSLbhnKpxzVR/6Lr7weuux5IjNJoiZQUVFlNLhSbT2t0u\nV+Tq7GN0vo0gdAelfqxYLyg+7nSJ1PicqTdyzlRNRhedCATqFD+vCoRDSPp1AEKyr5ecnqla/LDW\nW8P49IEn+EXtAGxQmoVul2sUTqvYO6Aetl1EVea1TOMDNdfodfPkKf1M1WH8fntTUjsF3ZdhE0r2\nlbV9didw92TqIfvCQzZDCFFpeSl8WJCUvItWncfQPGo0s99bYXUVu347infNvqz+puSbpVGHEbTv\n+gwdH3yWe3tPcHozKrZmeneFrs0+hZh2f8qlpGmfazRN3wSU2qh78abkij7oaqj+lNRUq2LaC8Az\nbrrxDA4eQ1qtlgkJiSSvmU3diHA69RpPXGx0qRHztVotCTM/IrZXJ5PXNRoNKevmEhoShCtUbDL1\nr9C12cf4ANN3fyrut2/StPe082meRN/XFHTJU59QKf5fv9/057zlXLd7cPB42rnnKM0aR9Koge5E\n+dD4HqzdsK1UMp3/4VqGxHVj129HSy2jyIUXtSQtlEVqL5VfptoBlIMM20UqFQeb+WkX06lft6bh\neb3IcNI0YaUnAAARxklEQVQuppuWuZDO2g1beW7MAACMf+FZo4Heg14hqufzLFqy3unNqNiaqbt/\nCIyb+flG/5t3DhfuqZDSrQlfdAlVfxpAX85azVTJDRjuwN2PJXtYOaZSbkNKGRUajXFmtGLytP/w\n1utPodFoKKLIpCb664Z3iagTxtX0TGLiE2jVoj7dotvZGXyJik2mlyp0bfYz7x5j/Nx8h0tSdS/m\nCTKfMi9cCDdiJQt194buRs9n5plOrxsRxrm0q4bn59KuUi8y3KTMnn3HGfr0LADSM26wIXkXPj5V\niOvXhYg6YQDUDA9m0ENd2bnnSCVKpu7e5DLvGuVtZZpwP/okanyhUM6LVg4OHltRHVtyPDWN02cv\nEVknjJVrUli+aJpJmdS9Sw3/Pzn+3wyIvY+4fl3Izc1Dqy2kenV/cnJu8cOm3fz9lZHObIWMtC+E\nUJmDydTbuwqJsyfQd8iraLWFPDUiltYtG7Jw8bdA6Z98NnbpynXiR80AoKCgkOGP9qRPzyjHAimm\nKcrYWCE3aWlCYyj6V0WsyQGWmvDm3Z+kZurezDvnm9dMlQ4+44E1Wc2rUJSxUe0wLNKExlBUV2HZ\nNPfdDqjommlV20VUYemgknvuKxf9BScvs+d6tpKpBybRSsNDji9p5pdFxiqt/LzxvC5uVSnpteAJ\nJJk6wF1rpmXxkB19R7DW48ITqwyV8ViyxkOOsYr9mFWWN01qo56jsnzm7mQe8mVXsZshSUpUNPnM\nuT8P+cKzuRn2jsry1Tp3H+9LCOFWnBg1yp2UGaJ+VJakVbM4tO1jlq/exOGjZyyW04/K4sqBA4QQ\ndwAPSaZlNvNdMSqLicrS5KoEO04oVFk+c3cyDzneytwMZ0dlEUJ1HnKgmvC0bboTaqbOjspSSiV4\nQ4QH8sTPnSdt051wNd/ZUVnMzTAaMrB7c91DCOFaKcd1j0rDQ74YykymzozKYsmM/i6IWAhRJvOK\nyswN6sWiyJ2QTJ0ZlUUIIRTxkGRasaNGza+INQkhjGkmuu9oS5rQGIruV1j2V/fdDvCYU79CiErL\nQ2qmHrIZQohKy4muUbbu0Fy7fisduo2j44PPck+P5/lp817F89pLaqZCCHU5mIX0d2gmr5lN3Yhw\nOvUaT1xstMlNRb0f7MjA/roL4n8cOsWgkTM4sWeJonntJTVTIYS6HKyZGt+h6ePjbbhD01hAgJ/h\n/+zsW4SHBime15HNEEII9TiYTJXcoQnw9Xe/0rrzGPr9dRrvvzXernntIc18IYS6rFTpUq5BSob1\n2ZTcoQnwyEP388hD9/PLtj8Y+dxsjuz4xIEgbZNkKoRQl5Vk2r2m7qE384TpdCV3aBrrFt2OggIt\nGddvUi8y3K55lZBmvhBCXQ42843v0MzPv83KNSnExUablDl56gJFxQOG/LZPd49tWGiQonntJTVT\nIYS6HKzSKblDc/W6X/hsRTI+PlUIDPBjxUfTy5zXGXIHlBAezu3vgBqssOxq990OkJqpEEJtHnKy\nUZKpEEJdkkyFEMIFJJkKIYQLSDIVQggXkGQqhBAuIMlUCCFcwEOykIdshhCi0pKaqRBCuIAkUyGE\ncAFJpkII4QKSTIUQwgUkmQohhAtIMhVCCBfwkCzkIZshhKi0PKRm6iGbIYSotBwcaR8gKXkXrTqP\noXnUaGa/t6LU9CPHzhLd5wV8I/rzduIqk2mNOoygfddn6Pjgs9zbe4LTmyE1UyGEuhys0mm1WiYk\nJJK8ZjZ1I8Lp1Gs8cbHRJiPmh4UGMX/2eL5ev7XU/BqNhpR1cwkNCXI0chNSMxVCqMvBmunOPUdp\n1jiSRg3q4OPjzdD4HqzdsM2kTM3wYKI6tsTH23K9sciFvzMiyVQIoS4Hk2naxXTq1y35+dJ6keGk\nXUxXvFqNBnoPeoWons+zaMl6JzZAR5r5Qgh1WclCKacg5bT12TQajVOr/XXDu0TUCeNqeiYx8Qm0\nalGfbtHtHF6eJFMhhLqstI+7N9U99GammE6vGxHGubSrhufn0q5SLzJc8Woj6oQBulMBgx7qys49\nR5xKptLMF0Koy8FmflTHlhxPTeP02Uvk599m5ZoU4mKjLa6iyOzkaG5uHjdv5gKQk3OLHzbtpl2b\nxk5thtRMhRDqcrBK5+1dhcTZE+g75FW02kKeGhFL65YNWbj4WwDGPfEwly5n0KnXBLJu5uDl5cV7\nC9dwaNvHXEnPJH7UDAAKCgoZ/mhP+vSMcmozNEUZG114PauMFYXGUDS/ItYkhDCmmei+vzevCY2h\naLbCsgnuux0gNVMhhNo85GSjJFMhhLo8JJkq2gxbt2wtW/UjHbqNo33XZ7g/dhL7D6a6PFAhhIfy\nVvhwczZDVHLLVpOGEWz+7h1qBAWQlLyLZybPY/tGOUEqhFDgTqmZKrllK/reNtQICgCgc1Qrzl9Q\nfheCEOIO58RAJ+7EZoj23rL18dIk+sfc65rohBCez0OSqc1mvj23bG365Xc+WZbEr0nvWpw+w+j2\n1+7NdQ8hhGulHNc9Ko1KkCiVsJlMld6ytf9gKmMnzyNp1SxCgqtbXNaM/k5EKoRQxLyiMnODerEo\n4iHJ1OZmKLll6+z5K8SPmsnnHyTQrEndcgtWCOGB7pRmvpJbtv4xZynXM2/y3EvvA+Dj483O5MTy\njVwI4RkqQbcnJeR2UiE8nNvfTrpUYdmR7rsd4DHfCUKISqsSNOGVkGQqhFCXJFMhhHABSaZCCOEC\nHpJMPWQzhBCVlhNdo2wNwgTwwtQFNI8aTYdu49i7/4Rd89q7GUIIoR4HR43SD8KUtGoWh7Z9zPLV\nmzh89IxJmfUbd3AiNY3ju5fw4bzJPDflPcXz2kuSqRBCXQ7WTJUMwvTNhm2MHtoHgM5RrcnMyuHS\n5QxF8zqyGUIIoR4Hk6mSQZjSLl6zWObCJcuvO0MuQAkh1GWlSpeyX/ewRukgTOa/TFpeJJkKIdRl\nJZl2v1v30Ju5zHS6kkGYzMucv5BOvcia3L5doGgAJ3tIM18IoS4Hm/lKBmGK6xfNZyt1t6Bu33WI\n4KAAatcKUTSvvaRmKoRQl4NVOiWDMPWP6cz6jTtpds9oAvx9+TTxpTLndYYMdCKEh3P7gU5+Vlj2\nQffdDpCaqRBCbR5yslGSqRBCXZJMhRDCBSSZCiGEC0gyFUIIF5BkKoQQLiDJVAghXMBDspCHbIYQ\notKSmqkQQriAJFMhhHABSaZCCOECkkyFEMIFJJkKIYQLeEgW8pDNEEJUWh5SM/WQzRBCVFpO/NSz\nNRnXs4gZlECLTk/QJz6BzBvZFsuNmTCX2i0fpd39Y01en/HWZ9S7axgdH3yWjg8+S1LyLkWbIYQQ\n6imHZPrWuyuJ6fEXju1aTK8HO/LWuysslntyeF+SVv2r1OsajYYXnx/M3p8/YO/PHxDbu5OizRBC\nCPWUQzL9JqnkJ55HD+3D1+u3WizXLbodIcGBFqfZ+0N8kkyFEOoqh2R6+cp1atcKAaB2rRAuX7lu\nd1jzF62lQ7dxPDXxbaunCYzJBSghhLocrNLFDErg0pWMUq+/+bcnTZ5rNBrFPwut99yYh3n9lREA\nvDZrMVP+tpCP508pc56KTaaFFbo2x0l93XNUls/cncxKFkrZontYs3HNbKvTatcK4dLlDOrUDuXi\npWvUqhlsV0i1aoYY/n96ZD8GDHvN5jxSM7XE/ACU5Fp5SPKsfKwcX90f0D30ZlrPnaXExUazZMUP\nJEwaypIVG3mkfxe7Qrp46RoRdcIAWPPtr7Rr09jmPJImhKhsCtQOwMXK4Zzp1MmPsTHlN1p0eoKf\nNu9l6uShAFy4mM5Dj003lBv29Jt0iZ3EsZNp1G/7OJ8uSwIgYeZHtO/6DB26jePnrfuZ9+ZzNtdp\n86eek5J3MXn6f9FqtTw9sh8Jk4aWKvPC1AVsSN6Jv58vixe8TMf2zUqvKDSGordtxuMevKz8Fe6t\n0ML/5n/17qB9qpnivj+RrAmNoShHYdkA990OsNHM12q1TEhIJHnNbOpGhNOp13jiYqNp3bKhocz6\njTs4kZrG8d1L2LH7MM9NeY/tG+dbXqDtC2LqMk6e3pByGro30/1f6p2ydJCqfICmHIfuzdWNwV4u\ni1m/HwqMnuebPTfmxL5KOQXdbbf6hFIe8sVW5mbs3HOUZo0jadSgDj4+3gyN78HaDdtMynyzoaQ/\nV+eo1mRm5ZTdDeFCGc8dnWZrHUqnXTF9mnLM6Ekl2OEpx9WOwH4ujznf7Pklo//PljHNjs9eyh+O\nzeeyz6yjx4y7KodmvhrKrJmmXUynft2ahuf1IsPZseeIWZlrpcqcv3DV0MfLxAWzv+avOzPN2rqU\nTiuujXIFqAr4Fr+eT8m7pN+phVi+0FFIpdjpHsV4P+hroflAbvEDdEkzr/j/s5Tsay8cTqiqf2Y9\nKaF6yDFTZjJV2jfL/E4Bq/PtVxaUKrzRJVGAICAU3UGZjS6xVjUql0/pZiV4zIei0jBPpPr9BZCO\nLlGmFz/PMirvjW6fGn/67Tm3ehn3/ixXNh7Sp6jMzagbEca5tKuG5+fSrlIvMrzMMucvpFM3wrSM\nnmabxZfd2sxf1I7APjM3qB2B/SplzOfVjsBz2Nmf3m2VmUyjOrbkeGoap89eIrJOGCvXpLB80TST\nMnH9oklctJahg3uwfdchgoMCLDbx3fkqnBBCHZ6UF8pMpt7eVUicPYG+Q15Fqy3kqRGxtG7ZkIWL\nvwVg3BMP0z+mM+s37qTZPaMJ8Pfl08SXKiRwIYRwJzb7mQohhLDN5ZdMkpJ30arzGJpHjWb2e5bH\nEHxh6gKaR42mQ7dx7N1/wtUh2MVWvMtW/UiHbuNo3/UZ7o+dxP6DqSpEaUrJewyw67ejeNfsy1fr\n1D3xqyTelC376Pjgs7TtMpbuA8oeUKIi2Io5/doNYoe8yt0PjKNtl7Es/uJ7FaIsYW2QY2PudNx5\nIpcmU30n/6RVszi07WOWr97E4aNnTMoYd/L/cN5knpvynitDsIuSeJs0jGDzd++wf8uHvPbSCJ6Z\nPE+laHWUxKwvlzDzI2J7dcLOYRldSkm8mTeyGf/yfNZ98QYHti7iy8WvqxStjpKYExetpWP7Zvy+\neSEp6+Yy5bWFFBRoVYrY+iDHeu503HkqlybTcunkX46UxBt9bxtqBAUA0DmqFecvpFtaVIVREjPA\n/A/XMiSuGzXDaqgQZQkl8X7x5U8MHtCVesX9lcMrQcwRdULJuqm7DzLrZi5hoUF4e1dRI1yg7EGO\nwb2OO0/l0mRqqZN/2sV0szKWO/mrQUm8xj5emkT/mHsrIjSrFL3HF9JZu2Erz40ZAKjb9URJvMdP\nppGReZMecS8R1fN5lq5U9wqvkpjHjurPwSNniGzzGB0eGMd7s56v6DDt4k7HnadyaXdZl3fyL2f2\nrHfTL7/zybIkfk16txwjsk1JzJOn/Ye3Xn8KjUZDEUWqNvOVxHu7oIDf9p3gx6/nkHvrT6L7vsB9\nUa1p3rReBURYmpKYZ72znLvbNiVl3ducPHWBmPgE9m1eSPXq/hUQoWPc5bjzVC5Npq7u5F/elMQL\nsP9gKmMnzyNp1SxCgqtXZIilKIl5z77jDH16FgDpGTfYkLwLH58qxPWzb0xHV1ASb/26NQkPrYGf\nXzX8/KrxQHQ79h1IVS2ZKol5686DTH/xcQCaNo6kccM6HD1xjqiOLSs0VqXc6bjzVC5t5ht38s/P\nv83KNSnExUablInrF81nxc24sjr5VwQl8Z49f4X4UTP5/IMEmjWpq0qcxpTEnLp3Kad+1z2GDOjG\nf+e+oEoiVRrvwH5d2LLjAFqtltzcPHbsOUIbo5HJKpqSmFs1b0Dyz3sB3e8NHT1+niaNItQIVxF3\nOu48lUtrppWtk7+SeP8xZynXM2/y3EvvA+Dj483O5ES3jtmdKIm3VYsGxPbsRPuu4/Dy0jB2VH/a\ntFIvmSqJedqLw3hywlw6dBtHYWEhc2aOJTQkSLWYhz39Jj9v3U/6tSzqt32cmVNHcbugwBCvOx13\nnko67QshhAvIOEdCCOECkkyFEMIFJJkKIYQLSDIVQggXkGQqhBAuIMlUCCFcQJKpEEK4gCRTIYRw\ngf8H9J1u2ROgApYAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plotitem = plotdata.getitem('ITEM1','AXES1','pcolor')\n", "plotitem.pcolor_cmin = -0.2\n", "plotitem.pcolor_cmax = 1.2\n", "plotdata.getfigure('contour').show = False\n", "plotdata.getfigure('cells').show = False\n", "plotdata.plotframe(2)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Explore the gauge output:\n", "\n", "The overshoots and undershoots with Lax-Wendroff are even more visible in the gauge output..." ] }, { "cell_type": "code", "execution_count": 21, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Gauge 1 at [0.6, 0.4] has a maximum value of q = 1.16098 at t = 1.84500\n", "Gauge 2 at [0.6, 0.8] has a maximum value of q = 1.15636 at t = 1.54062\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 21, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYYAAAEZCAYAAACTsIJzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmYHFW5/z9vT8++JpmQCcmEBNnXBAibgsPVi4jKqrIq\nekVBL4riVdxJrjsoP5VcFq+gIAJeN0AgAoIDCBhAEnaQLSE7CcnsMz3d0+/vj+qe6enMUt1T1V3d\n/X6eZ55Ud586/VbOTH3rXc45oqoYhmEYRpJQvg0wDMMwgoUJg2EYhjEKEwbDMAxjFCYMhmEYxihM\nGAzDMIxRmDAYhmEYozBhMAzDMEZhwmAEHhE5XURWiEiPiGwWkX+IyKfzbZeXiEi7iHwiqP0ZpYUJ\ngxFoROSLwE+AHwKzVHUWcD7wdhGpyKtx3uL1TFObuWpkjQmDEVhEpBFYCnxaVf+oqr0AqrpKVc9W\n1cFEu/eJyEoR6RSRN0TkkpQ+2kRkbVq/q0XkXYnjahG5XkS2icjzIvLl1PYisrOI/EFE3hSR10Tk\nsxPZKyI3JNquFpGvi4gkPlsiIr9OaTtfROIiUiYi3wWOApaJSLeI/CzRJi4inxWRV0Vki4hcOpX+\nDMMt4XwbYBgTcARQCdw2Sbse4GxVfU5E9gfuFZFVqjreecrIE/UlwDxgAVAHLE9+JiIh4M/An4DT\ngFbgryLykqreM0a/VwD1ib6agXuAjcB1jP8Er6r6dRE5Evi1ql6X9vlJwMGJfv8KvARcO4X+DGNS\nzGMwgkwzsFVV48k3ROQREdkuIn0ichSAqj6gqs8ljp8BbgHe6fI7PgR8T1U7VXU98FNAEp8tBppV\n9TuqGlPV14FfAKendyIiZTji8VVV7VXVNcCPgY8km7iwZaw2P1TVDlVdixNSO2OK/RnGpJjHYASZ\nt4BmEQklxUFVjwRIhHuSYZXDgB8A+wIVOF7G/7n8jp2B1FDTupTjXYCdRWR7yntlwINj9NMMlANr\nUt57A5jj0g4Y2wtIte2NhL1T6c8wJsU8BiPIPApEcMIpE3ETcCswV1WbgKsZ+d3uBWqSDRNP9jNT\nzt2IEyJKknq8FnhdVael/DSo6vvHsGErEAXmp7w3jxGhGWUH0JJ2/ng38Xlpx+un2J9hTIoJgxFY\nVLUDJ/l8pYicKiL1IhISkYVAbUrTOmC7qg6KyKHAmYzcGP8FVInI8SJSDnwDx6NI8n/AV0WkSUTm\nABeknPsY0J1ISFcnErv7icghY9g6lOjruyJSJyK7AF8Abkw0WQkcLSKtiaT6V9O62Ay8bYz/hv9K\n2NYKfA747RT7M4xJMWEwAo2qXgZcBHwZ2JT4uTrx+tFEs88A/y0iXcA3Gbl5oqqdic9/gfP03sPo\n8Mx/J95/HSdZ/DtgMHHuEPB+YCHwGrAF+DnQMI65n8V5kn8NeAj4DfDLRF9/Tdj1NPA4TlI79an+\np8AHE9VRP0l5/zbgnzhCcAdOInsq/RnGpIht1GMYIyQmzn1YVY8JgC1xYDdVfS3fthilhXkMRkkj\nIi0i8vZEiGpPHO/kT/m2yzDyiVUlGaVOBU5oagHQAdwMXJlXi0Ywd97ICxZKMgzDMEZhoSTDMAxj\nFAURShIRc2sMwzCyQFUzngFfMB6DqhbtzyWXXMLzbz7PnlfsyTl/Oofrnrwub7bceqtywgnKlVcq\n55/v3fXl+//Yz7FLHp/9x7P51cpf8eq2V9nl/+2SF3vicaWsTBkcVP72N+Xoo23sSvn6sqVghKHY\nGYgNUF1ezbSqaWwf2D75CT7R1wc1NTB9OmzbljczCpLuSDf1lfV5HcNoFEIhKC+HadNge/5+lYwC\nxoQhIPTH+qkKVzGtehodAx15s6O314QhW7oHu6mvqKexqpGewR6G4kM5tyE5fmBjaGSPCUMAaGtr\nczyGcDVNVU1s78+vx1Bb6+1Npa2tzZuOAkjqtXVHummobCAkIRoqG/Ii8MnxA288hmIeOyj+68sW\nE4YA0NbWRn804TEUYSipmP/4RgnDoBNKAvI2jqkeQ20tDA5CJJJ9f8U8dlD815ctBVGVVAoM5xiq\n8ysMqaGkt97KmxkFSc9gD3UVdQDOOObB80v1GESccdy+HVrS1171kcQmc0aOmUqyOR0ThoCQzDE0\nVTXlNcfQ1wezZ0NDg3McjTqJTGNy+qP9VIergWB4DDASTsqlMIC3NyljcrwWYwslBYRkjmFaVX6e\nNJMknzhFrKolUwZiA1SFq4BgeAzgjKEloI1MMWEICMM5hoCEkgBmzLCbiltUlf5YP9XlwfIYkqEk\nw8gEE4aAECSPwcodMycajxKSEOGQE53N1zj29prHYEwdE4aAkMwxVIWrUJSB2EBe7EgNRZgwuCc1\njATkzfNLFXawcKCRHSYMASFZlSQiefUa+vqg2omGmDBkQGriGYLjMVgoaUduueUWDjvsMOrq6pg1\naxaHH344V111Vb7N2oG//e1vHHPMMTQ1NbFgwYKcfrcJQ0BI5hggf0+b4NS8VyUefK1k1T2p+QVw\nxnDbQO5VdSyPwcR9hB//+Md8/vOf5+KLL2bz5s1s3ryZq6++mocffpjBwcF8mzeKuro6zj33XC67\n7LKcf7evwiAi14nIZhF5ZoI2PxORl0XkKRFZ5Kc9QaY/NvLEmc/Zz4ODUFHhHJvH4J4dQkljeAyx\neIyuSJevdoxXrppKdzdcey3EYr6aEjg6Ozu55JJLuOqqqzjllFOoTbhWCxcu5MYbb6Qi8Yt/5513\nsmjRIhobG5k3bx5Lly4d7qO9vZ3W1tZR/c6fP5/77rsPgP7+fs455xymT5/OPvvsw6WXXjqq/YYN\nGzj11FPZaaed2HXXXbniiivGtXfx4sWcddZZOfcWwH+P4ZfAceN9KCLH4+xpuzvwKSB4/lyOGFXq\nWJW/9ZIGB0fmLYwlDAF7qAoMO4SSxvD6vnbf12j5UQv3v36/b3YMDOxYlZQ+hv/5n/DVr8KXv+yb\nGYHk0UcfJRKJcOKJJ07Yrq6ujhtvvJHOzk7uvPNOrrrqKm677bZx24vI8DyCpUuX8sYbb/D6669z\n7733cuONNw5/Fo/H+cAHPsCiRYvYsGED9913Hz/5yU+45557vLtIj/BVGFT1IWCiR98TgOsTbVcA\nTSIyy0+bgsqoUseUm0rvYC8/fuTHPPfmczmxYyKP4aKLYOed4bHHcmJKQZEeSppePX2Ux9AV6eLq\nJ67mivdewXl3nOdbcUF//0goEHb0GAYG4NZb4aGH4Prr4bXXfDFjQkS8+cmUrVu30tzcTCg0cts7\n8sgjmTZtGjU1NTz00EMAvPOd72TfffcFYP/99+f000/ngQcecPUdv/vd7/ja175GY2Mjc+bM4cIL\nLxye7Pf444+zdetWvvGNbxAOh1mwYAHnnnsut9xyS+YX4zP5nvk8B1ib8nodMBfYnB9z8keqxzCj\negZb+7YCcEn7Jfxj3T/4wcM/4Pvv+j4LmhZQFioDnNBEdChKbUUte87Yk1l1U9fUaHRsYejvh5//\nHP7nf+Dkk+FPf4K5c2HrVujqcp5S99+/dGdJp4eSZtbMZEvfFlQVEeHBNQ+yeM5iPnHQJ7jrlbv4\nwd9/wJK2Jd7bMTBaGNLF/e9/h/32gz33hAsvhG99C2680XMzJiRfk6JnzJjB1q1bicfjw+LwyCOP\nANDa2jp8A1+xYgVf+cpXeO655xgcHCQSifDhD3/Y1Xds2LBhVOho7ty5w8dr1qxhw4YNTJs2bfi9\noaEhjj766Clfm9fkWxgA0rV/zF+bJUuWDB+3tbUV3eJXqTeWuQ1zWdu5lqH4EDc9cxP3n3M/XZEu\nLmm/hEgsQlzjKEo4FKY8VE7PYA/Pb3mey99zOR9b+LEp2ZHqMaROcPvHP5wb/znnQDgMp53m3IRm\nzIDGRujshA0bHPH44AenZEJBkh5KqquoQ1XpGeyhvrKeB1Y/wDHzjwHgp8f9lEXXLCIkIaJDUWbV\nzeLTh3x6WPCnQrowzJwJW7aMvF6xAt7xDuf4ootg993hkUfgyCOn/NWB54gjjqCyspJbb72VU045\nZdx2Z555Jp/73Oe4++67qaio4Atf+AJbtzoParW1tfT19Q23HRoaYkvKf/Ds2bNZu3Yte+21FwBr\n144897a2trJgwQL+9a9/eX1pw7S3t9Pe3j7lfvItDOuB1EzO3MR7O5AqDMXI4NAgFWXOHbm1oZUn\nNjzBqk2raKpqYq9m55ds+VnLxz3/pa0vcdgvDuPkvU6msaoxezvGCSU9+CAkH2zOOsv5SefJJ+Hd\n74ZjjnEEo5RIDyWJCC11LWzu3Ux9ZT0rN63ki0d8EXCE/68f+SvXP3U99RX1XP/U9cTiMT5/+Oen\nbEe6MMyY4SSbIxGorISnn4YPfMD5rK7OEfL3vQ8WLIBvf9s5Llaampq45JJL+MxnPoOqcuyxx1Jb\nW8vTTz9Nb2/vcLuenh6mTZtGRUUFjz32GDfddBPvec97ANhjjz0YGBjgrrvu4t///d/53ve+RyRl\n+doPf/jDfP/732fx4sX09vaybNmy4RzDoYceSn19PZdeeimf/exnqaio4IUXXmBgYIBDDjlkB3tV\nlUgkQjQaHT4WkeEk+VikPzSnJs4zId/lqrcDHwUQkcOBDlUtuTASQHQoSnnIicO0NraytmstD699\nmHfMe4er8/ds3pO2+W386cU/TcmO9ORzslw1VRjG46CDnJvOdddNyYSCJD2UBDCrbhabezajqjy1\n+SkWtiwc/uzAlgO5/D2Xs/SYpVx7wrX88OEferKxT7owhEKw006wOfFX9fTTcMABI59/4AOwcSN8\n73vwkY+MtCtWvvSlL3H55Zdz6aWX0tLSQktLC+effz6XXnopRxxxBABXXnkl3/rWt2hoaODb3/42\np5122vD5jY2NXHnllZx77rnMnTuXurq6UaGjb33rW8ydO5cFCxZw7LHH8qEPfWj4Rl5WVsYdd9zB\nqlWr2HXXXZk5cyaf+tSn6Ooau1LtgQceoKamhve9732sXbuW6upqjjtu3Foeb/F5v9GbgQ3AIE4u\n4T+A84DzUtosA14BngIOGqcfLXYO/d9D9dG1j6qq6urtq3Xu5XP1Q//3If3Vyl+57uMX//yFnv3H\ns6dkR3W1ak+Pczw0pFpRodrZqVpXp7p9++Tn33ab6rHHTsmEguTnT/xcP3HbJ0a9d+LNJ+rvn/u9\nrutcpzMvnanxeHzc8/f5n310xboVU7bjmGNU//rX0e8dfLDqihWq/f2qVVWqkcjY537606pLlkzZ\nBC2Fv1e3XHnlldrW1ub794z3f554P+N7t99VSWeo6s6qWqGqrap6napeo6rXpLS5QFV3U9UDVfVJ\nP+0JMqkew871O7Ouax23vXQbx77tWNd97D9r/ylXL6Umn0Mh2GMPJzm5227Q1DT5+Ycf7sSx4/Ep\nmVFwpM5DSZIMJT21+SkObDlwwqWR37vbe/nLK3+Zsh3pHgM4y6hv2gTPP++M43iRiFNOgXvvnbIJ\nJc2mTZt4+OGHicfjvPTSS1x++eWcfPLJ+TYrY/IdSjISRONRysscYSgvK+fI1iMZHBpkdv1s133s\n3bw3L259MeuQRDzuTHoKp2Se9t/fqUR65zvd9bHTTk4I6qWXsjKhYBkrlNRS18KG7g08tekpDpx1\n4ITnHzbnMJ7cOPXnorGEoaXFCRelh5HSOeIIWLXKmT1tZMfg4CDnn38+DQ0NvOtd7+Kkk07iM5/5\nTL7NyhgThoCQ6jEA3PuRe1nz+TUZ9VFfWc/M2pm83vF6djYkvIXUB9t993WeNDMpAttzT3j11axM\nKFj6o6OTz+AI9XNbnmPlppWj8gtjccCsA3jmzXEXCHDNWMIwdy688QY888zEwlBbC3vt5QiIkR3z\n5s3jmWeeoaenh3Xr1nHZZZcRDue7xidzTBgCQqrHAFBTXsO8xnkZ9zO/aT5rO9dO3nAMUhPPSc44\nw1k+4YQTMrBhPqxenZUJBUtkKEJlWeWo9xbNXsTKjSt5cM2DkxYR7DZ9Nzb1bKJnsGdKdowlDAce\nCE89BY8+CmMUv4xizz1Lz9szdqTwpKxISfcYsmVW7Sw292ZXWpJaqppk112dn0yYPx/WZObsFDyR\nWITp1dNHvbfb9N1Y07mG5ppm5jfNn/D8slAZu0/fnZe2vsTBOx+ctR1jCcNBBzkznWMxJ1w0EXvt\nZcJgmMcQGNI9hmyZVTuLN3vfzM6G6PiJyUwoRY9hcGhwB48hJCFuPvVm7j77bld97NK0C290vjEl\nO8YShtZWZ3b64sWj11EaC/MYDDCPITB45jEkauezYSyPIRt22aX0hCEyFBmeoJjK6fud7rqPeQ3z\nfBEGEXjzTWeG+mTsthu88sqUTDCKAPMYAoJXHsNOtTt5GkrKhlmznBtRKREZilAZrpy84QTMa5ya\nMKg6wlA5hhkzZ7ob2zlznKVNjNLGhCEgeJljyDaUNFbyORuam0tvg5+xQkmZMq9xHm90ZS8M0SiU\nlY0uN86UmTOdda9SVnkwShAThoDgWY6hztvkczbU1Tk3qYH8bFudFyKxsUNJmdDa2Dolj2GsMFKm\nhEIj8x6KlULZ2vOyyy5j//33p6GhgV133ZUf/ehHOftuE4YAoKrE4jFPPIbp1dPZ1p/dtmteJZ9F\nSs9r8CKUNKt2Flt6t0zecBwGBkb2654KO+8M68dcyrLwKaStPQF+/etf09HRwV/+8heWLVvGb3/7\n25x8rwlDAIjFY5RJ2YRLJrilvqKe7kh3Vud65TGAs6pnYqXiksCLUFJzTTNb+qYmDFP1GKB48wyF\ntrXnl770JRYuXEgoFGKPPfbgxBNP5OGHH/bs/2MiTBgCgFdhJICGyga6B/MvDM3NpSUMXoSSGiob\niMQiWe/uNl7iOVNmzy5OYSjkrT1VlQcffJD99tsvgyvOHitXDQBeJZ7BmTE9EBsgFo8RDmU2vF4l\nn6EEhcGDUJKI0FzTzNa+rcxtmDv5CWl45THMmOFvGFCWTt0zBtBLMtsKbrytPV944QUikQh33303\nRx11FO9MWRgsdWvPyQQFnK09r776ahobG2lsbOTCCy8c3ksmdWtPYNTWnsceO/Fimck+Pv7xj2d0\nzdliwhAAvPQYRIT6inp6BntoqnKxHGoKXnoMTU3Q0eFNX4WAF6EkgJm1M7MWBq/Gb8YMfye5ZXpD\n94pC3dpz2bJl3HjjjTz00EOU52jvXAslBQAvPQZwFtPrioy9+ceEdniUfAZoaHB2DisVvAglQSLP\nkGUCOrlL25RtKNLCgdStPSfizDPP5KSTTmLdunV0dHRw/vnnE0+sI+92a88kY23tuX379uGfrq4u\n7rjjjnFtue6667j00ku577772HnnnTO+5mwxYQgAXnoMkMgzZJGA9tJjqK93lmEoFbwIJQHDoaSs\nbPBIGPwOJeWL1K09//CHP9Dd3U08HmfVqlWTbu2ZzBOkbu0ZjUb5zne+M+bWnh0dHaxfv37crT37\n+/sZGhri2Wef5YknnhjT3t/85jd8/etf55577mH+/Pn+/ceMgQlDAIgORTPOB0xEfUV2HoPXwlBK\nHoNXoaTmahMGPymkrT2/+c1vsm3bNhYvXkx9fT319fU529vBcgwBIBr3PpSUTWWSl8nnhoYS8xg8\nCiU1VTXRGenMzgYTBleceeaZnHnmmeN+fuqpp3LqqaeO+/k555zDOeecM/z6i1/84vBxTU0NN9xw\nw/Drq666apRwzJ49m5tuusmVna+99pqrdn5gHkMAiA55H0oyjyG3eBVKaqhsoHMg/8JQShVlXlIs\nW3uaxxAAvJr1nCTbSW6WfM4er0JJjVWNvLo9u+3vIhFvxq++3il9jUa98yBLheTWnq+//jpNTU2c\nccYZBbm1pwlDAPAj+RwEj6FUQkmq6lkoqbGyMauxA2f8vPAYREZCgTNmTL2/UiK5tWehY6GkAOB1\nuWpNeQ290d7JG6ZhoaTsiMVjhCREWahsyn01VDbkPccAzt4NndmZYRQBJgwBwGuPoTpcndWyCpZ8\nzo7BoUFP8gvghJLynWMARxhKZfyMHTFhCABeewzV5dX0R/szPs88huwYb/e2bJhKKMk8BsMrLMcQ\nAPzwGPpjmQuDl8nn2lpImSBa1ERiEU8Sz5DwGAIQSmpomJoweLFSsJE/TBgCQDF6DFVVzo0qHnc2\nfylmvAwlTbVctb7eEzOm5DEk1xwyCpci/5MtDILiMXgpDCKOOPRnbkbB4WUoqb6int5oL3GNZ3yu\nl+NnOYbSxoQhAPjiMWQpDF7WrdfUlEY4yctQUlmojJrymqzmoQQplGQUNr4Kg4gcJyIvisjLInLx\nGJ83i8hfRGSViDwrIh/z056g4ovHkOdQEpSOMHgZSgKoq6jLqtzYks+GV/gmDCJSBiwDjgP2Ac4Q\nkb3Tml0ArFTVhUAb8GMRKbm8h9ceQ1W4Ku/JZ3CEwUJJmVNbXkvvYH6Foa4OejM3wSgS/PQYDgVe\nUdXVqhoFbgHSt0DaCDQkjhuAt1Q15qNNgcTrRfSCkHyG0vEYvAwlAdRW1ObdY6irg54eb/oyCg8/\nn87nAGtTXq8DDktr87/A/SKyAagH3G2TVGR4vYheEJLPANXVpSEMXoeSzGMw8o2fwuCmZu1rwCpV\nbRORtwH3isiBqrpD5i255ylAW1sbbW1tXtmZd/zwGPI98xlKyGPwOpQ0BY/By3ko5jEUHu3t7bS3\nt0+5Hz+FYT3QmvK6FcdrSOVI4LsAqvqqiLwO7AnssKVRqjAUG754DBZKyhmeh5Ky9Bi8WkQPLJRU\nqKQ/NC9dujSrfvzMMTwB7C4i80WkAjgNuD2tzYvAuwFEZBaOKORvd4o84UuOISDJ51IQBs9DSQHJ\nMVgoqXTxzWNQ1ZiIXADcDZQB16rqCyJyXuLza4DvAb8UkadwROrLqrrNL5uCSjF7DFaVlDlByTGY\nx1C6+FoaqqrLgeVp712TcrwV+ICfNhQC0XiU2vJaz/pL5hhUNaM1ayyUlB2+hJLy7DFYjqG0sZnP\nAcBrjyEkIcrLyokMRTI6z5LP2eHV7m1JaiuC4TFYKKl0MWEIAF7nGCC7cJKVq2aHL6GkAHgMvb1g\n6+GVJiYMAcBrjwGc2c+ZegxeJ58rK52bVbETiUW8Tz5nWZXk1fiFw85PKYyfsSMmDAHAD4+hoqyC\nSCzzUJKXwlBV5WwqX+x4HUoKwlpJYAnoUsaEIQB4vYgeQGW4ksGhwYzO8UMYSuGJsxhDSQB//KN3\n+zsYhUXJLVgXRLxeRA8SHkOek8+VlaXhMXgdSqopr6EvmnlyxmthOPpo7/oyCgvzGAJANB4lHPJW\noyvLguExlIIweB1KqgpXZVw4EI/D0JC3wm6ULiYMAcCP5HOmOQZViMW8vbFYKCk7slnrKrlOkm21\nbHiBCUMA8CP5XBmuzCiUFI06ouDljaVkQklD3oaSslkd12tvzyhtTBgCQCwe88VjyCSU5MeNxUJJ\n2ZHNfhpe5xeM0saEIQD4kXyuLKvMKJTkdeIZSiiUFPM4lJSFx2DCYHiJCUMACEK5qnkM2eN1KCmb\n5LMJg+ElJgwBIAjlql7PeobSyTEMDg16nnw2j8HIJyYMAcAXjyHDclW/PIZSCSV5mmMIZ1eVZMJg\neIUJQwDwzWPIMMdgoaTs8DqUVFFWQXQoylB8yL0NJgyGh5gwBICgeAxeJ59LJZTktccgIlSFqzLy\nGqxc1fASE4YAEIQcg4WSssdrjwEyzzOYx2B4iQlDAPCrKimTUJIfyedSCSV5nXyGzPfTMGEwvMSE\nIQD45THkO/kcDjvr98Ri3vYbNLwOJYF5DEZ+MWEIAH7lGPIdShIpjXCSL6GkDCuTTBgMLzFhCAC+\nzHwOwAQ3cG5Wg5kt8lpw+OExZDrJzYTB8BIThgDgh8eQTbmqH0s2V1SUgDAEIPlsVUmGl5gwBAC/\n1krKxGPwI/kMxS8MQ/EhVJUyKfO030yTzwMDTtjOMLzAhCEA+OYx5DnHAMUvDElvQTzeCCHTeQwW\nSjK8xIQhz8Q1Tlzjnj9xBiXHUPTC4EN+ATLfTyMSMY/B8A4ThjyTDCN5/cRpHkNu8CO/AJkvmz4w\nYB6D4R2uhEFEjhKRjyeOZ4rIAn/NKh38CCNBMPZjgBIQBp88hmxCSeYxGF4xqTCIyBLgy8BXE29V\nADf6aFNJ4UfiGaC8rJxoPOreDks+Z4Ufs54h83ko5jEYXuLGYzgZOBHoBVDV9UC9m85F5DgReVFE\nXhaRi8dp0yYiK0XkWRFpd2l30eCXx5BcodMtFkrKDt9CSRkuaWIeg+ElYRdtIqoaT8bARaTWTcci\nUgYsA94NrAceF5HbVfWFlDZNwP8A71HVdSLSnOkFFDq+eQyhcks+5wDfks/mMRh5xI3H8DsRuQZo\nEpFPAfcBv3Bx3qHAK6q6WlWjwC04nkcqZwJ/UNV1AKq61b3pxYFfHkOmoSQThuwIksdgwmB4xaQe\ng6peJiLHAt3AHsA3VfVeF33PAdamvF4HHJbWZnegXET+hhOe+qmq/tqV5UWCXx5DNqEkSz5njp8e\nQ89gj3s7LJRkeIibUBKqeg9wT4Z9q4s25cBBwLuAGuBREfmHqr6c4XcVLL55DBmGkiz5nB1B8Rgs\nlGR4yaTCICLdY7zdCTwOfFFVXxvn1PVAa8rrVhyvIZW1wFZV7Qf6ReRB4EBgB2FYsmTJ8HFbWxtt\nbW2TmV4QRIeihEOu9DkjLJSUG4JSlWQegwHQ3t5Oe3v7lPtxc0f6Kc4N/ObE69OBtwErgeuAtnHO\newLYXUTmAxuA04Az0trcBixLJKorcUJNl4/VWaowFBPReHBCSSYMmePnzOdM5jGYx2DAjg/NS5cu\nzaofN8JwgqoekPL65yKySlUvFpGvjneSqsZE5ALgbqAMuFZVXxCR8xKfX6OqL4rIX4CngTjwv6r6\nfFZXUqBEh4IRSjJhyA6/QklV4SrzGIy84UYY+kTkNOB3idcfBJKPMhPmEVR1ObA87b1r0l7/CPiR\nK2uLkFg8FogJbpZ8zg5fy1Utx2DkCTflqmcBHwHeTPx8FDhbRKqBC3y0rSQIygQ3Sz5nR2TIFtEz\nig835aoVKNT1AAAdpUlEQVSvAu8f5+O/e2tO6WET3AobX5PPNo/ByBNuqpJ+mfaWAqjqf/hiUYlh\nE9wKm0jMx3LVDGc+m8dgeIWbHMOdjOQSqnHWTtrgm0Ulhp8eQyweQ1VdLentpzD0uJ+nVXD4Fkoy\nj8HII25CSb9PfS0iNwEP+2ZRieGXxyAihENhovGoq1CHJZ+zIxKLUF1e7Xm/2XgMJgyGV2SzUc8e\nwEyvDSlV/PIYwPEa3CagLfmcHX56DG7nMcRiIAJh7+dJGiWKmxxDDyOhJAU2A2MuoW1kjl8eAziV\nSYNDg9Qy+YK4lmPIDr+Sz1XhKtehJPMWDK9xE0qqy4UhpYqvHkMGCWg/hSHiPiJScAQh+WylqobX\nuHI+RWQazkqow79+qvqgX0aVEn4tiQGZhZLMY8iOICSfzWMwvMZNKOmTwOdwFsFbCRwOPAr8m7+m\nlQZ+LYkBI6EkN1jyOTt8XV3VPAYjT7hJPl+Is+nOalU9BliEs7qq4QG+egwBCSUVtTAEYEkMK1U1\nvMaNMAwklsVGRKpU9UVgT3/NKh389BisKsl//Eo+h0Nh4hpnKD40aVub3GZ4jZscw9pEjuFW4F4R\n2Q6s9tWqEsJPjyHTUJIJQ+b4FUoSkeFwUk2oZmIbzGMwPMZNVdLJicMlItIONAB/8dOoUsJXj8FC\nSb7jVygJRuYy1JRPLAzmMRhe4yb5PD3l5dOJf91s22m4IBqPTvqHny1uV1hVdUJJfiSfKyuLXBh8\n8hjA/VwG8xgMr3GTY3gS2Iqz3ebLieM1IvKkiBzsp3GlgN85BjehpGjUmTXrYkmljDGPIXvcViZZ\nuarhNW6E4V7gvao6Q1VnAMcBdwD/CVzlp3GlQBCqkvxKPEMJCMNQxJfkM7ivTLJyVcNr3AjDEap6\nd/KFqt6TeO9RwKfbSeng9zwGN6Ekv/ILUPzCMDg06Fsoya3HYKEkw2vcVCVtFJGLgVsAAT4MbBaR\nMpx9mo0p4PfMZzehJBOG7PE7+ezGY7Dks+E1bjyGM3FmPd8K/AmYB5wBlOGIhDEF3C6LnQ1uQ0l+\nzXqGEhAGH5PP5jEY+cJNueoWxt/b+RVvzSk9/JogBRZKygXmMRjFSDb7MRgeMjg0GIiqJBOG7PDb\nY3CzJ4N5DIbXmDDkGT89hvKQ+1CSCUPmqKqv41cVrnJdrmoeg+ElkwqDiLxjjPfe7o85pUd0yL8c\nQxBCSWVlEI/D0ORL/hQcg0ODhENhQuLP81Um5armMRhe4uY3+oox3lvmtSGliq8eQ5n7qiS/ks8i\njuhE3a3MUVAMxAaoDnu/33OSTJLP5jEYXjJu8llEjgCOBGaKyEU4paoA9VgIyjOKPZQEI+GkYrt5\nDcQGqC73URgySD6bx2B4yURVSRU4IlCW+DdJF/BBP40qJYq9KgmKN8/QH+unKuyf2lWWWbmqkR/G\nFQZVfQB4QER+paqrc2dSaRGEUJLfN5ZiFYaB2IC/whC2clUjP7gJCfWJyI9E5C4R+Vvi5343nYvI\ncSLyooi8nJg9PV67xSISE5FTXFteJAQhlGTCkB390X5/cwwuPQYTBsNr3AjDb4AXgV2BJTib9Dwx\n2UmJJTOW4Sy6tw9whojsPU67H+Ls8eDD+p7BZnBo0NeNetyEkkwYsiMoHkN/P1T7p09GCeJGGGao\n6i+AQVV9QFU/Dvybi/MOBV5R1dWqGsVZa+nEMdp9Fvg9sMWt0cWEhZIKF7+Tz1XhKlcT3Pr6TBgM\nb3EjDMk/6U0i8n4ROQiY5uK8OcDalNfrEu8NIyJzcMQiuXx3yW0A5OtaSRZK8pWgJJ/7+6HGn72e\njBLFzeqq3xWRJuCLOHMaGoAvuDjPzU3+J8BXVFVFRCjRUFK+q5L8LncsVmHwex6D2z27LZRkeI2b\nRfT+nDjsANoy6Hs9zqqsSVpxvIZUDgZucTSBZuC9IhJV1dvTO1uyZMnwcVtbG21tmZgSXHwPJcXd\nhZL8TF4WqzD0R332GFxOcLNQkpGkvb2d9vb2KffjxmPIlieA3UVkPrABOA1nue5hVHXX5LGI/BL4\n81iiAKOFoZjwvSrJks++4Xvy2eUEN/MYjCTpD81Lly7Nqh/fhEFVYyJyAXA3ziS5a1X1BRE5L/H5\nNX59d6GQXITNzx3c3Cafa2t9McGxo4iFIQhLYliOwfAaPz0GVHU5sDztvTEFIVHtVFLE4jFfF2Fz\nu1FPJALTp/tiAuAIQ2Ty+1vBkZPkswuPwUJJhtdMKgwi8kWcRHIyMTzqWFUv98m2osfPOQxgoSS/\n8X2tJBceQyzmrF7r1yKIRmnixmM4GFgM3I4jCO8HHgf+5aNdJYGf+QXILJRkwpA5/dF+X4XBzfgl\n8wtScvV8hp+4EYZW4CBV7QYQkUuAu1T1LF8tKwH8nMMA7kNJfi+pUFlZnMIwEBtgWrWbKT3Z4SaU\nZPkFww/cBLd3AlLvLtHEe8YU8dtjsFCSv/THfF4ryUUoyfILhh+48RhuAB4TkT/ihJJOAq731aoS\nwUJJhU0QylWtVNXwAzcT3L4rIn8BjsJJPH9MVVf6blkJ4LvHkEFVkglD5vTH/M0xuPEYTBgMP3BV\nrqqq/wT+6bMtJYeFkgqb3sFeasv9mwDixmPo67Mcg+E9tkVnHvFzchsEK5RUjPMYeqO91Fb4Jwxu\nxq+319/JiUZpYsKQR4ISSrJF9LKjL9pHTbl/j+tuQkk9PSYMhveYMOSRXCSfgxBKKtZy1SCEknp7\noa7ONxOMEsWEIY/4vTVkecjdRj1+18IXq8fgdyipMlzJ4NAgquOvYG+hJMMPTBjyiN/ljm5DSX5X\nthStMAz2+hpKCkmIslDZhGPY02Meg+E9Jgx5xO+1dtyGkvyubClWYeiL9vkaSoLJw0nmMRh+YMKQ\nR/xenTOTUJJ5DJmhqvRG/fUYYPLKJBMGww9MGPLIQGyAqjL/hCEcCjOkQxPGqGMx56fCvxx4UZar\nRoYilEmZr+XGMHllkoWSDD8wYcgjfoeSRMSZ5DZBjDqZePZzdc5i9Bj6on2+Jp6TWCjJyAcmDHnE\n7+QzJPZ9niAUkYslFYqxXNXvUtUkbjwGEwbDa0wY8ojfm8nD5Mti5EIYitFj8LtUNYkbj8FCSYbX\nmDDkkVx4DBVlFROGknKx1k4xCoPfs56TWPLZyAcmDHnE783kIRihpGIUhqCEkrq7ob7edzOMEsOE\nIY/kJMcQkFBSsVUl9Qz2BCKU1NkJjY2+m2GUGCYMecTveQxgoSS/6Ip00Vjp/x15Mo+hqwsaGnw3\nwygxTBjyiFUlFS6dkU4aKv2/I0/kMaiaMBj+YMKQR/yexwDBCSUVmzAEwWPo63P+b8v9nWNnlCAm\nDHkkKKEkE4bM6RzIjccwUVVSZ6d5C4Y/mDDkkSCEknKxpEIxCkNXpIvGqhx4DBOEkrq6LPFs+IMJ\nQx7JSbnqJKGkXJQ7FqUwDHblLscwTijJPAbDL0wY8kiuJrhN5DF0d/t/cynGctXOgc7c5RjMYzBy\njAlDHsnJkhiTbNaTC4+hrMypoBka8vd7cklXxDwGo3jxXRhE5DgReVFEXhaRi8f4/CwReUpEnhaR\nh0XkgLH6icf9tjT3dA92U1/p7105CKEkkeIrWc1ZueoEHsO2bTB9uu8mGCWIr8IgImXAMuA4YB/g\nDBHZO63Za8DRqnoA8G3g52P1FYv5aWl+6I50U1/h713ZTSgpF0sqFFueIVfJ54nG7623YMYM300w\nShC/PYZDgVdUdbWqRoFbgBNTG6jqo6ramXi5Apg7VkfFFIYAGBwaJBaPlUQoCYpPGDoGOnKTY5gg\nlGTCYPiF38IwB1ib8npd4r3x+ARw11gfFJvH0B3ppqGyAfFzhxzchZJyEacuJmGIa5xt/duYUeP/\nXXmiUJIJg+EXYZ/7H39PyTRE5BjgP4C3j/X5d7+7hKrEw3VbWxttbW0emJc/cpFfgMlDSV1d5jFk\nSudAJ7XltVSU+bgfaoLJPIbmZt9NMAqI9vZ22tvbp9yP38KwHmhNed2K4zWMIpFw/l/gOFXdPlZH\n//VfS4rqjyBXVS2Tbe2Zy1BSsZSsbunbQnNNbn4ZJ1oSY+tW8xiM0aQ/NC9dujSrfvwOJT0B7C4i\n80WkAjgNuD21gYjMA/4InK2qr4zXUTGGkvxOPEMix5DnqiQoLo9ha99WZtbOzMl3TTTz2UJJhl/4\n6jGoakxELgDuBsqAa1X1BRE5L/H5NcC3gGnAVYl4e1RVD03vq9iSz92D3XlfaycadXYAy8UkqWIq\nV93SmzuPYaLx27IFZuZGn4wSw+9QEqq6HFie9t41KcfnAudO1k+xeQxdka6c5BgmCiVt2wbTpkEo\nB9Mci85jqMmRxzBOKKmvz/kxj8Hwg4KZ+VxswpCrUFJFWcW4oaStW3OXvCwmYchpjmGcUNL69TBn\njjN50DC8pmCEodhCSTlLPk+wuqoJQ3Zs6tnErNpZOfmu8TyGpDAYhh8UjDAUm8cQhKqkXMaoi6kq\naU3nGnZp2iUn3zWex7BuHcwdcyqoYUwdE4Y8sbl3c06eOi2U5D1rOtawS2OOhGECjyE6715i8SL7\nwzACQcEIQ7GFkjZ0b2B2/Wzfv2eiUNKWLblLXhZTVdKazjXMb5qfk++qClcxEBvY4f1Vb7zK76uP\n5eZnbs6JHUZpUTDCUGwew8aejcyu818YJpo5u2YNzJvnuwnA2B7D44/DzQV2X+uOdDMQG8hZ8rk6\nXD2mMDzScyOt1Xvym2d+kxM7jNKiYISh2DyGjd0bc+IxjBWK2Ni9kee3PM9rr8Hb3ua7CcDYwvDV\nr8KZZzqeS6Hw6vZXmd803/c1rpJUhavoj/bv8P6m2Iucvd85rNq0Kid2GKVFwQhDMXkMcY2zuXdz\nTjyGqnDVqOSlqjL3/83lyGuP5NVX8ycMqvDKK7DHHuDB0i45Y+XGlSxsWZiz76su39Fj2LIFhurW\ncNw+7yAaj7Kxe2PO7DFKAxOGPPBW31vUVdRRGa70/bsqyypH3Vg2926mpryGuMbZuL2D1tYJTvaQ\ndGHYtg06OuDss+GJJ3Jjgxes3LSSg1oOytn3VYWr6I/1ozqyHuVjj0HZjNUsmDafRS2LzGswPKdg\nhKGYQknPb3mePWbskZPvqgpXjQolre5Yzd7Ne7Nf0xHMPOQBystzYsYO5aobNjh1+IsXF5YwrFi/\ngoNm504YwqEwIQmNqj665/4IQ5VvsXP9zixsWWjCYHhOwQhDMXkMK9av4LA5h+XkuyrDoz2G1R2r\n2aVpF+o7D2PWQY/nxAbY0WPYuBFmz4aDD4Z//rMwtm5d37Wel996mbfPG3NleN9Ieg3ghOD+/MBa\nZlXPoSxUxqKWRazctDKn9hjFT0EJw0MPwfLlk7cNOo+ue5RD5+ywTqAvpE+QWtOxhvmN83lz1WKG\nWnInDOnlqklhmDkTmpqcfEPQuXbltZy010k52YchldTKpOXLIVS/hdbpOwGwaPYintz4ZE7tMYof\n3xfR84qhIfjOd+Cee+CZZ2C//fJtUXY8ufFJ/v7G37n2hGtz8n3poaTXO1YTXb8vW1YeRs/CjxAd\nilJeVk4kFqFjoIOQhKgpr6G6vJqQePfcUFHhbAqUJCkMAIcc4pSu7pGIrq1e7Uy+q6+HlhZnh7l8\nrQmkqnRGOrnn1Xv42Yqf8fgncyemSarCVbyxsZ8n/gWf/CR86kfbWBGfDsCeM/ZkU88mOgc6h/eg\nXv7ych5Z+wiz62ezV/Ne7DtzX3aq3SlnlVRG4VMwwnDOM7PpPRiqDoUDb3BWBI3reBvEjb9xnIqO\n3UQy72u8z1R07I8S31F5z9W0Lpk+0j6lbfJ4rPeyaaszKuH0gZEb62mbmb313fzhVztxwTNv49F1\nj1IVruLk355MdChKXOP0x/rpj/ZTGa50RCJcTVmobDgBquioY+e7dPg4+Tr1874BZagafnmZ87qn\nXwnVwXWXQv8Byh0vwWd/qAwMwMCAEgo55zohJh2xf3icdOzXKd/ppu24n0nKv4N1lG3dn5qHb+Xg\nHyzIaGzSj7M5Z+CTVRx7/AB7z4Rly6Bn1228/Krz+1MWKuOAWQewctNK2ua3ccuzt3DxXy/mYwd+\njFWbVnHLs7fw3JbnUFWaqpqoLq8eXiYlFo8RHUr8G48yFHcSeSEJISJcsPgCvn701zFKj4IRhh/M\nf5JvfBPuvht6uiEac2bthlKeglL/mAQZfi/1QUmEkSen0f8gIju0hdHfsUNfibOHuxz+V0Y6TvmO\nxspGwl8vG7Ov9OPJPnfTdk1nFf92Q4TXhpz3/+367XzzndM4YgGcOHgiX/nrV1jdsZqr3ncVJ+51\n4vB5qspAbIC+aB990T7iGh+5LpzrTj1Ofiak/l+NfH7DDfD008KPv+y8/uSn4L3vFU49BV5/Hd53\nvHDFFXDRRcI998DOs0f67OuD7m4hHndyEfEhZ5xCIadNKDFuoZTxk1CqnSM2htLeH9VGRvobvh4Z\nudZMxmOi40zPOexX1dzwtX4Wtjivf/qPbUyvHnmwOGreUdz/+v0csvMhXHT3Rdx2+m0snrN4+HNV\n5a3+t+iKdNEX7Rv2EstD5YRDYcrLnH/LxPm9VJS4xqktr8UoTQpGGIY6Z9P/JizazZYazoTqcmeC\nW3LPhe0D25lWNQ2A/zryv1jTsYYLD7twlCiAc1OsLq+muryaGUx93YxZdUAv7JS41/S+CQt2guYa\naN4XDjsA/uMMuPNO2HfX0ec2VQHT03ssHarLRy+L8Vb/W8yoHhmT9+/xfi5YfgGqyjELjhklCuCM\nZXNNc85maxuFT8EIwyuvOKtJmihkRvpaO9sHttNU1QRATXkN156Ym1xHdTX0p0zg3b7d2SQoyR/+\nAJs25W6JjkKiOlw9avbztv5t7Dljz+HXR7QeQW15Ldetuo5HP/FoPkw0ioyCEYaNG53qFSMzKsOj\nq5K2929nWvW0Cc7wh8mEoaLCRGE80sV9W//oUFI4FOa+j95Hf6x/1PuGkS0FIwybNuVmb+JiIznz\nWVUZ0iH6on052QcinXRh6OgwoXdLdXn18DwG2FEYkm2qy6tzbZpRpBTMPIaNG00YsqEsVEZZqIxY\nPEbnQCcNlQ2elqG6JVUY4nHo7DRhcMtkHoNheE3BCIN5DNmTXHp7+0B+wkgwWhi6u6GmBsIF46/m\nl7FyDDNqcrSRhlGSFIwwdHaaMGRL8olze/9IRVKuSRWG7dvNW8iEdI/hrf63zGMwfKVghAFMGLIl\nmYDOt8fQ1+ccd3SMTjwbE1MdrqYv6vznDcWH6I5001hpfwyGf5gwlABB9BhMGNxTV1FHb7QXgI6B\nDhoqGygLleXZKqOYMWEoAVJzDMk5DLnGhCF76irq6BnsASzxbOQGE4YSIBlK6hjoyLvHoGo5hkxJ\nFwZLPBt+Y8JQAlSHnTr4fE1uA6cCKRx2lt62HENmpAqDJZ6NXFAwwnD88XDEEfm2ojBpqGygO9I9\nap2kfJD0GiyUlBmpOQYLJRm5wFdhEJHjRORFEXlZRC4ep83PEp8/JSKLxuvrzjuhqso/W4uZhsoG\nOiOdea1KAhOGbNkhx1BlwmD4i2/CICJlwDLgOGAf4AwR2TutzfHAbqq6O/Ap4Cq/7CllGiob6Ip0\n5bUqCUaEwZbDyAxLPhu5xk+P4VDgFVVdrapR4BbgxLQ2JwDXA6jqCqBJRGb5aFNJMiwM5jEUJCYM\nRq7xUxjmAGtTXq9LvDdZm7k+2lSSpHoM+SpXBWerzq4uE4ZMSU8+W1WS4Td+rlYz0Z6YqaTvsDDm\neUuWLBk+bmtro62tLSujSpGGygbWda1ja9/WvG7WMn26IwomDJlhHoPhlvb2dtrb26fcj5/CsB5o\nTXndiuMRTNRmbuK9HUgVBiMzGiob2NSzicGhwbwupTB9OmzbZvMYMsWEwXBL+kPz0qVLs+rHz1DS\nE8DuIjJfRCqA04Db09rcDnwUQEQOBzpUdbOPNpUkDZUNvLztZVrqWkb2u84D06fDli2OODTbLpOu\nqQpXMTg0SCweM2EwcoJvHoOqxkTkAuBuoAy4VlVfEJHzEp9fo6p3icjxIvIK0At83C97SpmGygb+\n9da/2Lt578kb+8j06fDSS04Yqbw8r6YUFCIy7DWYMBi5wNcV8VV1ObA87b1r0l5f4KcNxkjyuaWu\nJa92TJ8O998Ps2fn1YyCpLmmmc09m+kc6MxrybFRGhTMzGcje5Jbee5Uu1Ne7Zg+HZ57zoQhG1rq\nWnj2zWdpqmqylVUN3zFhKAH2at4LYNS+wfkgWZVkwpA5LXUtPL7hceY1zsu3KUYJYJsrlgDhUJit\nX9qa9yfN6YnQ+E75dVwKkpbaFh5b/5gJg5ETzGMoEWbUzMjr5DaAhQvh9NPhXe/KqxkFSdJjaG1o\nnbyxYUwR8xiMnFFdDTffnG8rCpOWuhZ6BnvMYzBygnkMhlEA7DZ9NwD2mLFHni0xSgFRdbtyRf4Q\nES0EOw3DTwZiA1SFbe15wz0igqpmPKvVhMEwDKNIyVYYLJRkGIZhjMKEwTAMwxiFCUMA8GKZ3CBT\nzNdXzNcGdn2liglDACj2X85ivr5ivjaw6ytVTBgMwzCMUZgwGIZhGKMomHLVfNtgGIZRiBTtPAbD\nMAwjd1goyTAMwxiFCYNhGIYxikAJg4gcJyIvisjLInLxOG1+lvj8KRFZlGsbp8Jk1ycibSLSKSIr\nEz/fyIedmSIi14nIZhF5ZoI2hTxuE15foY5bEhFpFZG/ichzIvKsiHxunHYFOYZurq9Qx1BEqkRk\nhYisEpHnReT747TLbOxUNRA/QBnwCjAfKAdWAXuntTkeuCtxfBjwj3zb7fH1tQG359vWLK7tKGAR\n8Mw4nxfsuLm8voIctxT7W4CFieM64KUi+9tzc30FO4ZATeLfMPAP4B1THbsgeQyHAq+o6mpVjQK3\nACemtTkBuB5AVVcATSIyK7dmZo2b6wPIuIIg36jqQ8D2CZoU8ri5uT4owHFLoqqbVHVV4rgHeAHY\nOa1ZwY6hy+uDAh1DVe1LHFbgPIBuS2uS8dgFSRjmAGtTXq9LvDdZm7k+2+UVbq5PgSMT7t5dIrJP\nzqzzl0IeNzcUzbiJyHwc72hF2kdFMYYTXF/BjqGIhERkFbAZ+JuqPp/WJOOxC9IObm7rZtNVvVDq\nbd3Y+STQqqp9IvJe4FagWHZmKdRxc0NRjJuI1AG/By5MPFnv0CTtdUGN4STXV7BjqKpxYKGINAJ3\ni0ibqranNcto7ILkMawHUje0bcVRtonazE28VwhMen2q2p10C1V1OVAuItNzZ6JvFPK4TUoxjJuI\nlAN/AG5U1VvHaFLQYzjZ9RXDGKpqJ3AncEjaRxmPXZCE4QlgdxGZLyIVwGnA7Wltbgc+CiAihwMd\nqro5t2ZmzaTXJyKzREQSx4fiTEBMjxcWIoU8bpNS6OOWsP1a4HlV/ck4zQp2DN1cX6GOoYg0i0hT\n4rga+HdgZVqzjMcuMKEkVY2JyAXA3TgJlGtV9QUROS/x+TWqepeIHC8irwC9wMfzaHJGuLk+4IPA\np0UkBvQBp+fN4AwQkZuBdwLNIrIWuASn8qrgxw0mvz4KdNxSeDtwNvC0iCRvKl8D5kFRjOGk10fh\njuFs4HoRCeE86P9aVe+b6n3TlsQwDMMwRhGkUJJhGIYRAEwYDMMwjFGYMBiGYRijMGEwDMMwRmHC\nYBiGYYzChMEwDMMYhQmDYaQhIo0i8unE8WwR+V2+bTKMXGLzGAwjjcRCa39W1f3zbIph5IXAzHw2\njADxA+BtiVmyL+Os3b+/iHwMOAmoAXYHfgxUAWcCEeB4Vd0uIm8DlgEzcWbRflJVX8r9ZRhGdlgo\nyTB25GLgVVVdBHwp7bN9gZOBxcB3gS5VPQh4lMR6NMDPgc+q6iGJ86/MidWG4RHmMRjGjsg4x+Cs\nd98L9IpIB/DnxPvPAAeISC1wJPC7xJps4GygYhgFgwmDYWRGJOU4nvI6jvP3FAK2J7wNwyhILJRk\nGDvSDdRneI6As64/8LqIfBCcJZ9F5ACP7TMMXzFhMIw0VPUt4GEReQa4lJHdrpTRO1+lHydfnwV8\nIrHd4rM4e+4aRsFg5aqGYRjGKMxjMAzDMEZhwmAYhmGMwoTBMAzDGIUJg2EYhjEKEwbDMAxjFCYM\nhmEYxihMGAzDMIxRmDAYhmEYo/j/bLh+v/XCom4AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "for gaugeno in [1,2]:\n", " gauge = plotdata.getgauge(gaugeno)\n", " q = gauge.q[0,:]\n", " t = gauge.t\n", " qmax = q.max()\n", " tmax = t[q.argmax()]\n", " print \"Gauge %s at %s has a maximum value of q = %7.5f at t = %7.5f\" \\\n", " % (gaugeno, gauge.location, qmax, tmax)\n", " plot(t,q,label=\"Gauge %s\" % gaugeno)\n", "legend()\n", "xlim(0,3)\n", "ylim(-0.1,1.1)\n", "xlabel('time')\n", "ylabel('q at gauge')\n", "title(\"Gauge output\")" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": false }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.10" } }, "nbformat": 4, "nbformat_minor": 0 }