{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "import netCDF4 as nc\n", "import numpy as np\n", "import matplotlib.patches as patches\n", "from salishsea_tools import viz_tools, geo_tools, tidetools\n", "from bathy_helpers import *\n", "import matplotlib.path as mpltPath\n", "from fast_ll2ij_SalishSea201702 import *" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "grid = nc.Dataset('/data/vdo/MEOPAR/NEMO-forcing/grid/bathy_meter_SalishSea2.nc')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "# October 2016" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAMEAAAKvCAYAAADTI6ufAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2MXNd93vHvs7vcpURJlALTbkTSMINSUqU2spmtUAG1\nAEMVqRRpGMdSvLXgFzkB60JhDTSIIcF1ECAg0NQNAqaA7bBSBKOgLdNsnDAtTIkK0OSfiCpt0a74\nVrOSXS5h2eva9ItE7nK5v/6xl/aImnvm7szs7ix/zwcYcOaee2fODvnw3L3nnnMUEZhlNrTcFTBb\nbg6BpecQWHoOgaXnEFh6DoGl1ygEkm6UtF/SSUknJN0t6QuSjlaPb0o62rL/Y5JOSzoladviVd+s\ndyMN99sNHIyIBySNAtdGxHsvF0r6I+CH1fPbgQngDuBm4FlJt0TEpf5W3aw/OrYEktYC9wBPAETE\nTEScaykX8BvA56tN24GnImI6Il4GTgN39bviZv3SpCXYBEwBT0q6E/gK8NGIeLUqfyfwnYj4RvV6\nPfBcy/GT1bbXkbQD2AGwZs2aX7rtttu6+wnsqveNF75ZXzgkfjT7ve9FxLpu379JCEaALcDOiDgs\naTfwKPCJqvxf8rNWoLGI2APsARgfH48jR44s9C0sifvXfri2TNdew8FXPvWtXt6/yS/Gk8BkRByu\nXu9nPhRIGgF+HfhCy/5ngY0trzdU28wGUscQRMQrwBlJt1ab7gWOV8//GXAyIiZbDjkATEgak7QJ\n2Aw838c6m/VV06tDO4G91ZWhl4CHq+0TXHEqFBHHJO1jPiizwCO+MmSDrFEIIuIoMN5m+4dq9t8F\n7OqpZmZLxD3Glp5DYOk5BJZe01+MzZbP8PCivr1bAkvPIbD0HAJLzyGw9BwCS88hsPQcAkvP/QQ2\n+BZ5qlC3BJaeQ2DpOQSWnkNg6TkElp5DYOk5BJae+wls4MX58/WFc3M9v79bAkvPIbD0HAJLzyGw\n9BwCS88hsPQcAkvPIbD03FlmA+/pC3try7atfqjn93dLYOk5BJaeQ2DpOQSWnkNg6TkElp5DYOk5\nBLailfoQmnIILD2HwNJzCCw9h8DScwgsPYfA0nMILD2HwNJzCCw9h8DScwgsPYfA0nMILD2HwNJz\nCCw9h8DScwgsPYfA0nMILD2HwNJzCCw9h8DScwgsPYfA0nMILD2HwNJzCCw9h8DScwgsPYfA0nMI\nLD2HwNJzCCw9h8DScwgsPYfA0msUAkk3Stov6aSkE5LurrbvrLYdk/QfWvZ/TNJpSackbVusypv1\nw0jD/XYDByPiAUmjwLWS3gVsB+6MiGlJbwaQdDswAdwB3Aw8K+mWiLi0CPU361nHlkDSWuAe4AmA\niJiJiHPAvwb+fURMV9u/Wx2yHXgqIqYj4mXgNHDXYlTerB+anA5tAqaAJyW9IOlxSWuAW4B3Sjos\n6W8k/eNq//XAmZbjJ6ttryNph6Qjko5MTU31+GOYda9JCEaALcCnI+IdwKvAo9X2nwP+CfC7wD5J\navrBEbEnIsYjYnzdunULr7lZnzQJwSQwGRGHq9f7mQ/FJPDnMe95YA54E3AW2Nhy/IZqm9lA6hiC\niHgFOCPp1mrTvcBx4C+AdwFIugUYBb4HHAAmJI1J2gRsBp5fhLqb9UXTq0M7gb3VlaGXgIeZPy36\nM0kvAjPAByMigGOS9jEflFngEV8ZskGm+X+3y2t8fDyOHDmy3NWwFUrSVyJivNvj3WNs6TkElp5D\nYOk5BJaeQ2DpOQSWnkNg6TkElp5DYOk5BJaeQ2DpOQSWnkNg6TkElp5DYOk5BJaeQ2DpNR1eueIs\nYOILGyDLMdLRLYGl5xBYelft6VCrdk3sfUMPdv1+h+a+2Et17ArLferqlsDScwgsPYfA0nMILD2H\nwNJLcXWolytBdvVzS2DpOQSWnkNg6TkElp5DYOk5BJaeQ2DpOQSWnkNg6TkElp5DYOk5BJaeQ2Dp\nOQSWnkNg6aUYT9Bv3Y5P8CwVgylFCNr94/NAG7vMp0OWnkNg6TkElp5DYOk5BJaeQ2DpOQSWXop+\ngkFx3/B7224fGh0tHjc3M1MsP3TpC13XyRyChVN946mh+nn2Y27plyGyZnw6ZOk5BJaeQ2DpOQSW\nnkNg6TkElp4vkV4F6vofYOX1ISx0nIeGh3v+zLQhKI3yKv1FlPoCSn8hGi2UjZT/GoZirrZs7uJs\n8VjrzKdDlp5DYOk5BJaeQ2DpOQSWnkNg6aW9RJrFIPUhDOpcTw7BApXGBXTdFxAdxhqU+h8uXap/\n26tkDINGVtWXrRqBHrtKfDpk6TkElp5DYOk1CoGkGyXtl3RS0glJd0v6fUlnJR2tHv+8Zf/HJJ2W\ndErStsWrvlnvmv5ivBs4GBEPSBoFrgW2AX8cEf+xdUdJtwMTwB3AzcCzkm6JiPrf4MyWUceWQNJa\n4B7gCYCImImIc4VDtgNPRcR0RLwMnAbu6kdlzRZDk5ZgEzAFPCnpTuArwEersp2SPgAcAX4nIn4A\nrAeeazl+stpmA6a2D6Fw6zZ0fxv6oGryO8EIsAX4dES8A3gVeBT4NPALwNuBbwN/tJAPlrRD0hFJ\nR6amphZW62Wk4eHaB5cu1T8iah8xN1d8MBf1j5VEQ+0frbu0+V6Hrlld/7jh+p6r1SQEk8BkRByu\nXu8HtkTEdyLiUkTMAf+Zn53ynAU2thy/odr2OhGxJyLGI2J83bp13f8EZj3qGIKIeAU4I+nWatO9\nwHFJP9+y27uBF6vnB4AJSWOSNgGbgef7WGezvmp6dWgnsLe6MvQS8DDwJ5LeDgTwTeBfAUTEMUn7\ngOPMd2g/4itDNsgahSAijgLjV2x+f2H/XcCuHupltmTcY2zpOQSWnm+ltgVbiX0BJQ5BG113BpXu\ne1dhvqLVY+UKzVysLxsurJdQGKcQhYU/lvMyhtosWFL8ftZc2/Nn+nTI0nMILD2HwNJzCCw9h8DS\ncwgsPYfA0nM/wQKV+hC2jr6vtixm6yfHGerUT3Bt/V9TqS+gNN6gOGym1IfQYQxDaf2G2jWgWyrT\nrp+Aa6+pf89Vvf8Tdktg6TkElp5DYOk5BJaeQ2DpOQSWnkNg6TkElp47y1aCDot91xqq/z9uaLR+\nANDcD39UWxbnL5Q/s65DjEKnYMuYobaLmRTeMwo/Y1NuCSw9h8DScwgsPYfA0nMILD2HwNJzCCw9\n9xMMgCgsyA2g0rXwVfXX+4v9C2NtBq9c/rzp6dqyocLgICj/LJ1+zrp9NFc4brgwiKchtwSWnkNg\n6TkElp5DYOk5BJaeQ2DpOQSWnvsJ+qm0EnxhkY7ePrMwGVapPhcKC3+UJtgq3NsPoOFicWft+hJK\n/R0eT2DWO4fA0nMILD2HwNJzCCw9h8DScwgsPfcTDAB1uPbetUv1/QRx/nx9WXGh7+5X+hYNjm23\n6EahLyD60P/ilsDScwgsPYfA0nMILD2HwNJzCCw9h8DScz/BUind919aALuX9y0o9U1Et2MUOhmu\nn+vosrb1ulg/9kHRYSH0BtwSWHoOgaXnEFh6DoGl5xBYeg6BpecQWHruJ+inbscFlOb5ge7/q+qy\nb0Kjhev5He7fj8LaBo3mXhpe2A8bXp/ArHcOgaXnEFh6DoGl5xBYeg6BpecQWHruJ+inLu+1jw7H\nicKk/6V5gEpz95fKSu/ZaR3jUp/HTGFNhMrcT159w7ah1fVjBmK41wUR3BKYOQRmDoGl5xBYeo1C\nIOlGSfslnZR0QtLdLWW/Iykkvall22OSTks6JWnbYlTcrF+aXh3aDRyMiAckjQLXAkjaCGwF/u/l\nHSXdDkwAdwA3A89KuiUiup/O2GwRdWwJJK0F7gGeAIiImYg4VxX/MfAxoPW62HbgqYiYjoiXgdPA\nXX2ttVkfNTkd2gRMAU9KekHS45LWSNoOnI2Ir12x/3rgTMvryWrb60jaIemIpCNTU1Pd1t+sZ01O\nh0aALcDOiDgsaTfw+8y3Dlu7/eCI2APsARgfH+9uBqkVRIUFqRdvkY76M9CYLizEsViLdDSYZEzt\nOsZKi3kv0aCaSWAyIg5Xr/czH4pNwNckfRPYAHxV0t8DzgIbW47fUG0zG0gdQxARrwBnJN1abboX\n+GpEvDki3hYRb2M+KFuqfQ8AE5LGJG0CNgPPL071zXrX9OrQTmBvdWXoJeDhuh0j4pikfcBxYBZ4\nxFeGbJA1CkFEHAXGC+Vvu+L1LmBXTzUzWyLuMbb0HAJLzyGw9Dyopp+6vd6/wAmnmorCAJhiX0Bh\nYEzHfoLSd9Dk+2m3AHnh5xh6tTDZV0NuCSw9h8DScwgsPYfA0nMILD2HwNJzCCw99xP0UfF++VWL\n9FWXFuJod839ssKCGSpMaNWpn6B0bKNFyxfaZ1KaRKzpW/T8DmYrnENg6TkElp5DYOk5BJaeQ2Dp\nOQSWnvsJBkHpWj+UF8EuXe8vLG5Rmgdp7rXX6o/rsChGzBYW4miwiElceOP4gFJdY2Rtx/fsxC2B\npecQWHoOgaXnEFh6DoGl5xBYeg6Bped+gn4qXUMvLXI90vsc+20V+h+K4wJK6xr0sD5Bk3mHtGrV\nGzdet6Z2/7nVbfZfILcElp5DYOk5BJaeQ2DpOQSWnkNg6TkElp77Cfqp1BdQmnOnNF6gB6X1CUp9\nAeU37TAmoNc1mdvVq/Rz9IFbAkvPIbD0HAJLzyGw9BwCS88hsPQcAkvP/QQLtG31Q/WFvV4jt2Xh\nEPSRCgtxlCaQ6kmXi3SUBsf0NHCm0JkWDd623Wdrun7h8aHz9WVN+b8uS88hsPQcAkvPIbD0HAJL\nzyGw9BwCS8/9BP002vtEUAtVGjgTFwsLZhSO0/AwT1/Y21V97ht6sPi+bbV0LbTtTyks2K3Zzgt/\ndOKWwNJzCCw9h8DScwgsPYfA0nMILD1fIm1j66qJ2rJO6/jayuMQLFRp4Exp8q3SBFudxhrMFa6F\nXyxNsFV/3FzpuMXSZDHv6Tcu5k1pnMZM7z+HT4csPYfA0nMILD2HwNJzCCw9h8DScwgsPfcTLFBx\nbqFCWU9KcwSVxhMUjtPwMM/MfK6narVzaO6LtWX3Db+38xu06Ywsztl0/kKTahW5JbD0HAJLzyGw\n9BqFQNKNkvZLOinphKS7Jf2BpK9LOirpGUk3t+z/mKTTkk5J2rZ41TfrXdOWYDdwMCJuA+4ETgCf\njIhfjIi3A/8N+D0ASbcDE8AdwP3ApyT51ksbWB1DIGktcA/wBEBEzETEuYj4Uctua4DLt1BuB56K\niOmIeBk4DdzV32qb9U+TlmATMAU8KekFSY9LWgMgaZekM8BDVC0BsB4403L8ZLXNbCA1CcEIsAX4\ndES8A3gVeBQgIj4eERuBvcBvL+SDJe2QdETSkampqQVWe3FpeLj2gVT/GBqqf4yM1D8uXiw+4vyF\n+sfF2frH7MXiYyBdutTmMVf/mO1hLYVKkxBMApMRcbh6vZ/5ULTaC7ynen4W2NhStqHa9joRsSci\nxiNifN26dQurtVkfdQxBRLwCnJF0a7XpXuC4pM0tu20HTlbPDwATksYkbQI2A8/3sc5mfdW0n38n\nsFfSKPAS8DDweBWMOeBbwEcAIuKYpH3AcWAWeCSiyUI9ZsujUQgi4igwfsXm97Tbt9p/F7Crh3qZ\nLRn3GFt6DoGl5xBYeh5P0E5pbqHCXPnF+YMK6w3Ha+eL1Zm70GYunhaLMS5gUTSZd6jN3E0xU79W\nsVaP9VQlcEtg5hCYOQSWnkNg6TkElp5DYOk5BJae+wnaGS6sm9tlX8DcuR8WP/Lpn3y2Y7VWuro5\niVRau6GT0toNDbklsPQcAkvPIbD0HAJLzyGw9BwCS88hsPTS9hNsu+b99YWFfgK7+qQNQUlx1frC\nII65H5wrvm+GDrFetZsULGYKfx8dBiQ14f/yLD2HwNJzCCw9h8DScwgsPYfA0nMILL28/QSlgTM3\nXF9/XB8GcVi9dgNvtq35QP0BHlRj1juHwNJzCCw9h8DScwgsPYfA0nMILL20/QTFMQOlyaBmCotg\nFybfsu7FxdnasuJkaA25JbD0HAJLzyGw9BwCS88hsPQcAksvxSXSdnMMaXTVMtTEBlGKECyYxwwM\nlNJi5cVJ1Bry6ZCl5xBYeg6BpecQWHoOgaXnEFh6DoGll6OfINpc9+/y3v9o916XXbrU1Xta96IP\n37lbAkvPIbD0HAJLzyGw9BwCS88hsPQcAksvRT/B0xf2vmHbtus+WLu/LtbPLVSaryhKcxnZomi3\n7vFCuSWw9BwCS88hsPQcAkvPIbD0HAJLzyGw9BwCSy9FZ1lbhUE1pYEauv66+rIL0z1VyZaHWwJL\nzyGw9BwCS69RCCTdKGm/pJOSTki6W9Inq9dfl/QlSTe27P+YpNOSTknatnjVN+td05ZgN3AwIm4D\n7gROAIeAfxgRvwj8b+AxAEm3AxPAHcD9wKck+fZKG1gdQyBpLXAP8ARARMxExLmIeCYiLq+t+Ryw\noXq+HXgqIqYj4mXgNHBX/6tu1h9NWoJNwBTwpKQXJD0uac0V+3wY+HL1fD1wpqVsstr2OpJ2SDoi\n6cjU1FQXVTfrjyYhGAG2AJ+OiHcArwKPXi6U9HFgFnjjyJWCiNgTEeMRMb5u3bqFHNofly7VP7qk\nm9YWHzaYmoRgEpiMiMPV6/3MhwJJHwJ+BXgo4qe9T2eBjS3Hb6i2mQ2kjiGIiFeAM5JurTbdCxyX\ndD/wMeBXI+K1lkMOABOSxiRtAjYDz/e53mZ90/S2iZ3AXkmjwEvAw8D/BMaAQ5IAnouIj0TEMUn7\ngOPMnyY9EhGepNMGVqMQRMRRYPyKzX+/sP8uYFcP9TJbMu4xtvQcAkvPIbD08o4nKIjCuACNFL4y\nT7619DQE3a238lNuCSw9h8DScwgsPYfA0nMILD2HwNJzCCy9tP0E7RbuuGzrqonassJS3rYMNDzc\n81+KWwJLzyGw9BwCS88hsPQcAkvPIbD00l4iXQ5bR99XW/bMzOeWsCYrS+mSNer9/3GHoJ3CF1ta\nzLtk7vyFbmtjBRoehh7X8/bpkKXnEFh6DoGl5xBYeg6BpecQWHppL5Fuu+b9y10FGxBpQ1AytHqs\ntkxrb6g/cLi+YR36/rleqpRaaV3pfnSW+XTI0nMILD2HwNJzCCw9h8DScwgsPYfA0nM/QTujq+rL\nxkbry6Zn6suix0n0bdG4JbD0HAJLzyGw9BwCS88hsPQcAkvPIbD0HAJLL21n2dPn/0tt2bY1H6g/\n8OL3uvq8ucIC4dZBp8nQPPmWWW8cAkvPIbD0HAJLzyGw9BwCSy/tJdJBU7osGzOFcQrAMxef6nd1\nls19Qw++cWMf5hYqcQjaKP6jK0wE1e0kUVrlv4aSocL3o9Vj0OP6Jz4dsvQcAkvPIbD0HAJLzyGw\n9BwCS88hsPR8gbqN0vX+Ullpoe+ha1bXf2CnBcLn8kzcpZE2E5+NFPoJShOlNeSWwNJzCCw9h8DS\ncwgsPYfA0nMILD1fIr0KbLvug7VlT//ks0tYk5XJIWjj0NwXuzpu66qJ+sLCIh0q9SEADKn+ba+y\n+YzU5mct9b8wWlg0pSGfDll6DoGl5xBYeo1CIOlGSfslnZR0QtLdkh6UdEzSnKTxK/Z/TNJpSack\nbVucqpv1R9NfjHcDByPiAUmjwLXAOeDXgT9t3VHS7cAEcAdwM/CspFsiojAK3Wz5dAyBpLXAPcCH\nACJiBphhPgRIb/htfjvwVERMAy9LOg3cBfxd32pt1kdNToc2AVPAk5JekPS4pDWF/dcDZ1peT1bb\nzAZSk9OhEWALsDMiDkvaDTwKfKKXD5a0A9gB8Na3vrWXtxoYpbEGpfUJNDfX9WdqKMG1jdK8TB0m\nJmuiyTc4CUxGxOHq9X7mQ1HnLLCx5fWGatvrRMSeiBiPiPF169Y1ra9Z33UMQUS8ApyRdGu16V7g\neOGQA8CEpDFJm4DNwPM919RskTS9OrQT2FtdGXoJeFjSu4H/BKwD/rukoxGxLSKOSdrHfFBmgUd8\nZcgGWaMQRMRRYPyKzV+qHu323wXs6q1qZksjwW9VZmUOgaXnW6mvch5r0JlDsFSivi+g07V+rbm2\nvrAwJ8/cD851rNagefrC3jdsu/+m36rdvx8zMvl0yNJzCCw9h8DScwgsPYfA0nMILD2HwNJzCCw9\nd5b1U2nB7naLT1wu6zT5VqGzLL5f6BCbnS2/7wqhQocgq8d6fn+3BJaeQ2DpOQSWnkNg6TkElp5D\nYOk5BJae+wn6qLSYhErXs9eUJvSjvNh3YQGPq0WUBiQVFj9pyi2BpecQWHoOgaXnEFh6DoGl5xBY\ner5EukBbR9+33FWwPnMI+kijhTEDpQm0xuqPA2C2flJvrSoc24eFrgeBCuM0KCyM0pRPhyw9h8DS\ncwgsPYfA0nMILD2HwNLzJdI27ht6sLasNHWKrUwOwQJ1PWagtCD1XPme+Cj1P1xXGItw8WLxfVeK\nL099prbsl9d9pOf39+mQpecQWHoOgaXnEFh6DoGl5xBYer5Emti21Q/VlrVbT/hq5RAsUMwWrr1f\nKKxPULjvP64vjDUAojDeQNMz9WWF+/Dnfvzj4meuFNGHNRh8OmTpOQSWnkNg6TkElp5DYOk5BJae\nL5FaW6X5lZ6Z+dwS1mTxOQQLFKV7/2fq+xDitfO1ZfpBh7+Gm66vLyuNUyitXVAYFxEXV876xwd/\n8DjSEz29h0+HLD2HwNJzCCw9h8DScwgsPYfA0vMlUluw0jiETgZxnIJD0MahuS/2/T3vv+Hh+sIO\n8wPpx6/VF5bWICiUabTQp1HoJ9BVuG6yT4csPYfA0nMILD2HwNJzCCw9h8DScwgsPYfA0nNn2RKZ\nuzBdW9bpfyKp0EEVhUE+hYmpYqZ+0q6iwmAcAEp1HVBuCSw9h8DScwgsvUYhkHSjpP2STko6Ielu\nST8n6ZCkb1R/3tSy/2OSTks6JWnb4lXfrHdNW4LdwMGIuA24EzgBPAr8dURsBv66eo2k24EJ4A7g\nfuBTkjr8NmW2fDqGQNJa4B7gCYCImImIc8B24LPVbp8Ffq16vh14KiKmI+Jl4DRwV78rbtYvTVqC\nTcAU8KSkFyQ9LmkN8JaI+Ha1zyvAW6rn64EzLcdPVtteR9IOSUckHZmamur+JzDrUZN+ghFgC7Az\nIg5L2k116nNZRISk8orUV4iIPcAegPHx8QUduxKVFveYu1A+tvg/VWFSr9ICFt1OsNXpuJU46KZJ\nSzAJTEbE4er1fuZD8R1JPw9Q/fndqvwssLHl+A3VNrOB1DEEEfEKcEbSrdWme4HjwAHgg9W2DwJ/\nWT0/AExIGpO0CdgMPN/XWpv1UdPbJnYCeyWNAi8BDzMfoH2SfhP4FvAbABFxTNI+5oMyCzwSEZf6\nXnOzPmkUgog4Coy3Kbq3Zv9dwK4e6mW2ZNxjbOk5BJaeb6W2JVVa/GO5OARLpDSh131DDxaPjZn6\na+/FRUNirr7oUnfXKtRhPEHxEkhhcfHlNJi1MltCDoGl5xBYeg6BpecQWHoOgaWX4hJpccqSlaB+\ntpalt8hLHN9/02+9YVucL6wBfc01PX+mWwJLzyGw9K7a06EozcxmRaUe7E49xiotH1U4LdWqVR3r\ntVjcElh6DoGl5xBYeg6BpecQWHpX7dUh614vYx8ozEs0dM3q2jKNFP4pXrem/rixMThXrlInbgks\nPYfA0nMILD2HwNJzCCw9h8DScwgsPfcT2IKU+hAA7ht+b21ZTBdGB62q/6dYvHO1y/mTWrklsPQc\nAkvPIbD0HAJLzyGw9BwCS88hsPQcAkvPnWXWV4cufaG2rLRKjWbqFztnTWE6ltVjTapV5JbA0nMI\nLD2HwNJzCCw9h8DScwgsPYfA0nM/gS2Z0gLiUZi0i4v1fQhyP4FZ7xwCS88hsPQcAkvPIbD0HAJL\nzyGw9NxPYEum27EGnL9QW9SPpV/dElh6DoGl5xBYeg6BpecQWHoOgaXnEFh67iewwVcch1CYr6gh\ntwSWnkNg6TkElp5DYOk5BJaeQ2DpOQSWnvsJbCCU5iTS8Gj9gZfmev5stwSWnkNg6TkElp5DYOk1\nCoGkb0r6X5KOSjpSbbtT0t9V2/9K0g0t+z8m6bSkU5K2LVblzfphIS3BuyLi7RExXr1+HHg0Iv4R\n8CXgdwEk3Q5MAHcA9wOfkjTcxzqb9VUvp0O3AH9bPT8EvKd6vh14KiKmI+Jl4DRwVw+fY7aomvYT\nBPCspEvAn0bEHuAY8//g/wJ4ENhY7bseeK7l2Mlq2+tI2gHsqF5OS3px4dVfNG8CvrfclWiRuz71\n0w5dLru1l7dvGoJ/GhFnJb0ZOCTpJPBh4E8kfQI4AMws5IOrIO0BkHSk5TRr2bk+ZYNYn16Ob3Q6\nFBFnqz+/y/z5/10RcTIitkbELwGfB/5PtftZftYqAGyotpkNpI4hkLRG0vWXnwNbgRerVgFJQ8C/\nAz5THXIAmJA0JmkTsBl4fjEqb9YPTU6H3gJ8SdLl/T8XEQclfVTSI9U+fw48CRARxyTtA44Ds8Aj\nEVF/Y8i8PV3VfvG4PmVXVX0UEf2qiNmK5B5jS88hsPSWJASS/kzSd1v7AiQ9KOmYpDlJ4y3b3ybp\nfHWLxlFJn2n/rn2vzyclnZT0dUlfknRjS9mi3gaykPos4/fzB1Vdjkp6RtLNLWWL9v0spC5dfzcR\nsegP4B5gC/Biy7Z/wHwnx/8Axlu2v611vyWsz1ZgpHr+h8AfVs9vB74GjAGbmL8UPLyM9Vmu7+eG\nluf/BvjMUnw/C6xLV9/NkrQEEfG3wPev2HYiIk4txec3rM8zEXF5WfXnmO/fgCW4DWSB9Vl0NfX5\nUcvLNczfRQCL/P0ssC5dGdTfCTZVzdnfSHrnMnz+h4EvV8/XA2daytreBrKE9YFl+n4k7ZJ0BngI\n+L1q87J8PzV1gS6+m0EMwbeBt0bE24F/C3yu9TbtxSbp48z3b+xdqs8saVOfZft+IuLjEbGxqstv\nL8VnLrAuXX03AxeCqln9f9XzrzB/jnnLUny2pA8BvwI8FNVJJst4G0i7+izn99NiLz+7a3i5b5P5\naV26/W45jipCAAAAs0lEQVQGLgSS1l0efyDpF5i/7eKlJfjc+4GPAb8aEa+1FC3LbSB19VnG72dz\ny8vtwMnq+ZJ/P3V16fq7WcyrDC2/wX+e+abqIvPnjL8JvLt6Pg18B3i62vc9zN+mfRT4KvAvlqg+\np5k/tz1aPT7Tsv/Hmf9f5RTwy8tZn2X8fv4r8CLwdeCvgPVL8f0spC7dfje+bcLSG7jTIbOl5hBY\neg6BpecQWHoOgaXnEFh6DoGl9/8BNUMwjKXSGukAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1, figsize=(16,12))\n", "ax.pcolormesh(grid.variables['Bathymetry'][:])\n", "ax.set_ylim(590,670)\n", "ax.set_xlim(115,135)\n", "viz_tools.set_aspect(ax)\n", "ax.add_patch(patches.Rectangle((118,598), 15, 59, fill=False, linewidth=3))" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "Oct = nc.Dataset('/ocean/vdo/MEOPAR/ariane-runs/Oct16/ariane_trajectories_qualitative.nc')\n", "latt = Oct.variables['traj_lat']\n", "lont = Oct.variables['traj_lon']\n", "bathy, lons, lats = tidetools.get_bathy_data(grid)\n", "mask = lont[:].mask" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": true }, "outputs": [], "source": [ "def still_inside(time, number):\n", " number_of_particles = np.zeros(time)\n", " for n in range(time):\n", " for m in range(number):\n", " if (mask[n,m]) == False:\n", " y,x = fast_ll2ij_SalishSea201702(lont[n,m],latt[n,m],lons, lats)\n", " if (598" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig,ax=plt.subplots(1,1,figsize=(8,6))\n", "time = range(744)\n", "ax.plot(time, op)\n", "ax.grid('on')\n", "ax.set_title('Number of Particles in Domain', fontsize=16)\n", "ax.set_ylabel('Number of Particles', fontsize=14)\n", "ax.set_xlabel('Time (h)', fontsize=14)\n", "ax.tick_params(labelsize=12)" ] }, { "cell_type": "code", "execution_count": 17, "metadata": { "collapsed": true }, "outputs": [], "source": [ "index_deep_particles=[]\n", "d = Oct.variables['init_z']\n", "for n in range(2694):\n", " if d[n] > 6:\n", " index_deep_particles.append(n)" ] }, { "cell_type": "code", "execution_count": 18, "metadata": { "collapsed": true }, "outputs": [], "source": [ "mask = lont[:].mask\n", "def deep_particles(time, index):\n", " number_of_deep_particles = np.zeros(time)\n", " for n in range(time):\n", " for m in index:\n", " if (mask[n,m]) == False: \n", " y,x = fast_ll2ij_SalishSea201702(lont[n,m],latt[n,m],lons, lats)\n", " if (598" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf4AAAGDCAYAAADK03I6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd8XnX5//HXldHsZjRNOhKadNLBbCllScsWweJCUFmC\nfFXc4sCFylC+/nB9FaSCgqIWRAREpkDYpbTM0r13OtI2o23aJNfvj3NabmJW09w5yZ338/G4Hzn3\n56zrc+4793XO53zOOebuiIiISN+QFHUAIiIi0n2U+EVERPoQJX4REZE+RIlfRESkD1HiFxER6UOU\n+EVERPoQJX7BzO40s+sjWreZ2R/NbJuZzY4ohlozG97ONGVm5maWEof1v2hmR4XD3zGz2zs4X4en\nbWX+qWa2trPzt7PsuG2vA4ih2MwWmFlaG9P8zsy+351xxUtHvsfdxcxOMrNFUcchLVPi74HMbKWZ\nbTKzrJiyK8ysIsKw4uVE4HSgxN0nNx9pZpeaWWP4o1ZtZm+Y2TmdXZmZVZjZFbFl7p7t7ss7u8yD\nYWbnAjXu/noYy43ufkU7s9F82qgTbfidPS2idZ9mZq+ZWZ2ZrTWz8wHcvRJ4BriytXnd/bPufl0n\n11thZrvNrCb8bs41s2+3taMRT/H6HnfkwCD87o2MieV5dx/T1bFI11Di77mSgS9HHcSBMrPkA5xl\nGLDS3evamOZld88G8oA7gHvNLP8A4zIz64nf988Cf446iKiZ2UAzs07MNw74K/BdIBc4ApgbM8lf\ngP/pkiBb9gV3zwEGA18HLgAe6UxdRLqNu+vVw17ASuDbQBWQF5ZdAVSEw2WAAykx81QAV4TDlwIv\nAr8AtgPLgePD8jXAJuCSmHnvBH4HPAnUAM8Cw2LGHxqOqwIWAec3m/dW4BGgDjithfoMAR4K518K\nfCYsvxzYDTQCtcCPWpj3UuCFmPdZYd0nAfnAw8BmYFs4XNJsm9wQbotdBEmgMVxnLfCbcDoHRobD\nGcDNwCpgB/BCWPaebU6QZO4ANgDrgOuB5HDcyHAb7gC2APe08jn3C+OKjfmHwN3NPudLgNXhsr7b\nyrSrw2lrw9dxLawvI/y8tgHzgW8Aa5t9Tv8It+cK4EvN1nUfcE/4HXkNOCIc92egKaxLLfDN9mJv\nIbZvEXzvfwSUH8D/yl+B69oYnwLsJOb73Gz8ncD14fBUYC1BAt8UfraXtbHsCsL/uZiyQ8L1nRO+\nTwN+CawPX78E0pqt75sx6zsPOBtYTPD/8p2YZU8GXib4n94A/AboFzM+9nt8J/Bb4N/h5/UKMKKN\nuvwd2EjwnX0OGB+WXwnsBfaEn+2/Wpj3uXDddeE0H99Xt2a/ad8A3gqnuwMoBh4N4/sPkB8z/RTg\npbCubwJTu+r3NVz+8eF3sjR8fwTB/8WhXbmenvqKPAC9WvhQgn+S04D7Y36UDjTxNwCXEbQcXE/w\n4/vb8IfojPCfLTuc/s7w/fvC8b8iTLYEiXZNuKwU4KjwH2ZczLw7gBMIWpDSW6jPc8AtQDpwJEFi\nOSUm1hfa2BaXxsSSQtAKUkOQeAcAHwEygZzwx+uBZttkNTA+nDeVln+sY38wfxtOMzTcdseH2+Q9\n2xz4J3BbuH2KgNnA/4Tj/kZwBJoU1vnEVuo2HqhrVvZD/jvx/54gaR8B1ANj25g2pY1t+VPgeaAA\nKAXmEf44h7HOBX5AsEMynGCH8cyYde0FPhpux6sJdg5SY7+zMetqM/ZW4ptCsBO5laCJ/iIgs53/\nleXAdcDbBMnwbqCg2TRvAR9sZf47eW/ibwB+HNbxbIIknt/KvP/1XYr5vt8UDv8YmBV+RwYSJLPr\nmq3vB+H6PkPwv/FXgu/zeIKdqfJw+onhNkoJt+8C4CutfI/vDLfj5HD6vwAz29iOnw7XuW9H5Y2W\ntlEb8+9fd0zdmif+WQTJfijBjs5rBL8n6cDTwLXhtEPD2M8m+F6eHr4f2Mq6HybYQWjp9XAbMd8Q\nrjcj/P58oSO/z4nwijwAvVr4UN5N/BMIkupADjzxL4kZd1g4fXFM2VbgyHD4ztgfBSCb4Mi4lGDv\n/flm8d0W8096J/CnNupSGi4rJ6bsJ8CdMbG2l/gbwn/iLeGPx3+1KoTTHglsa7ZNftxsmv3bKabM\nCY7Skwh+aI9oYdn7t3n441UPZMSMvxB4Jhz+EzCDmCP5VuI9AdjYrOyH/Hcyj20RmA1c0Ma0bSX+\n5cBZMe+v5N3Efyywutn01wB/jFnXrJhxSQSJ9qTY72wL26vF2NvZLmnA+QStSFXA7W1Muydc92iC\n7+0/gL80m+ZF4OJW5r+T9yb+Xbz3/2oTMKWVef/ruxSWzwR+Hw4vA86OGXcmwamt2PXtaynKCbfZ\nsTHTzwXOa2X9XwH+2fx7HFOv22PGnQ0sbG/bh9PmhcvKbb6N2pinI4n/kzHv/wHcGvP+i4Q77QSt\nP39utvzHiWml7IoXwc7WXIKk/xhgXbn8nvzqiec8JeTu8wj2Zr/didkrY4Z3hctrXpYd835NzHpr\nCX5whxCcgz/WzLbvewGfBAa1NG8LhgBV7l4TU7aKYK++o2a5e567F7r7FHf/D4CZZZrZbWa2ysyq\nCY608pr1M2grtuYKCY4+lrUz3TCCH40NMdvkNoKjOgiabg2YbWbvmNmnW1nONoIf+/ZsjBneyXs/\ntwMxhPduj1Uxw8OAIc0+5+8Q7OTsE/sdaSJoph7SzjoPOHZ3ryc4Sn+DILFPaGPyXQQ7J4vD7+2N\nBEkuVg7BjmNHbHX3hgONuZmhBP8/EGyf2O28ivdus63u3hgO7wr/tvh/amajzexhM9sYft9vJPjO\ntqZD297Mks3sp2a2LFzuynBUW8vujOb1au33aBjwsWbfxRMJ+lF0GXffS7BTMwG42cO9gb5Aib/n\nu5agCTA2Ue7rCJcZUxabiDujdN+AmWUTNAevJ/ixfzZMvPte2e7+uZh52/qHWQ8UmFlsgjuE4Lz4\nwfo6MIbgCKk/wakKCJJua7G1FesWgvP/I9pZ7xqCI/7CmG3S393HA7j7Rnf/jLsPIehYdktsj+cY\nSwn6HR7ITlBrOvKjtYGYz5ngc9hnDbCi2eec4+6xSTT2O5IElBB8vh1df5vMbICZfSG8rPNpglMt\n09x9ShuzvdVs3e+JI7zKYSTBeeK4M7NSgib558Oi9QSJbJ9DeHebHahbgYXAqPD7/h3e+13vrE8A\n0wlaGXMJWmuIWXZ3J8Q1BEf8sd/FLHf/aUsTm9mj4VU/Lb0ebW0l4f/dtcAfgZujuhojCkr8PZy7\nLyXoUPWlmLLNBInzU+He+qdpP1m152wzO9HM+hGcM53l7msIWhxGm9lFZpYavo4xs7EdjH8NwXnN\nn5hZupkdTtCp7+6DjBeCI7ldwHYzKyD4J25PJcH565ZibQL+APzczIaE2/a45j8I7r4BeILgx6K/\nmSWZ2QgzOxnAzD5mZiXh5NsIfjibWljfHoJOTSd3pLLt2Byuo63ruO8FrjGz/DC+L8aMmw3UmNm3\nzCwjrPsEMzsmZpqJZvbhMJl+hWDnZ1Y4rtXt2hFmdjnBkebJBB38St39W+6+oJ1Z/whcZmbDzSyT\noHXs4Zjxkwma1le1OHcXCVufTgYeJNiWj4Sj/gZ8L7xqoZDgfH5nv/s5QDVQa2aHAp9rZ/oDWW49\nwem/TIKWhFgd+WwP6vNv5m7gXDM7M/wepltwz4mSliZ29/eHByMtvd7f0jzhVRd3EnQyvJxgp7hT\nl3X2Rkr8vcOPCTqRxfoMQS/ZrQSdgF46yHX8lSBxVhEcsXwKIGyiP4PgMqX1BM2HNxGch+2oCwmO\nItYTdIq7dl9z/UH6JUHHnH3n/h/rwDy/Aj5qwQ2Dft3C+KsJzvm9SrAtbqLl/5OLCTrBzSdI7vfx\nblPkMcArZlZLcDXDl73166tvI+jEdlDcfSfhFQxh82hLR8k/ImhqXkGw4/LnmPkbgXMI+kmsINim\ntxMcAe7zIEGfj21hzB8Om0sh6LfxvXDdV3eiCi8T9Lz/mLv/O6b5u03u/geCPhWvhHWrJ2YnmeC0\n1O86EU9H/cbMaggS3y8Jzl2fFe5EQtCxdg5By8TbBB3aOnuzrKsJjs5rCDpN3nMQccf6E8G2W0fw\nfZ7VbPwdwLjws32glWX8ELgrnOb8gwkmPFiYTtCisZmgBeAbdG2++hLBqbnvh038lxHsQJ7Uhevo\nsawPndYQ6ZHM7EWCHsWvRx1La8zshwSdtz4VdSwdZWZFBJdVHuXuu6OOR6SniOx2miIScPcToo4h\nEbn7JqBDp6RE+hI19YuIiPQhauoXERHpQ3TELyIi0oco8YuIiPQhCdu5r7Cw0MvKyrpseXV1dWRl\nNb+iLjEkat0StV6QuHVL1HqB6tYb9bZ6zZ07d4u7D2xvuoRN/GVlZcyZM6fLlldRUcHUqVO7bHk9\nSaLWLVHrBYlbt0StF6huvVFvq5eZdehGVWrqFxER6UOU+EVERPoQJX4REZE+RIlfRESkD1HiFxER\n6UOU+EVERPoQJX4REZE+RIlfRESkD1HiFxER6UOU+EVERPoQJX4REZE+RIm/A15bvY23NjewcGN1\n1KGIiIgcFCX+DvjpIwv5+dx6zvrl8zy1oDLqcERERDpNib8DbvjQBL57bDopScbj72yMOhwREZFO\nS9jH8nalUcU5rMtP5tSxRcxaXhV1OCIiIp2mI/4DcHhJHqurdlJb3xB1KCIiIp2ixH8ARhVlA7Ck\nsibiSERERDpHif8AjC7OAWBJZW3EkYiIiHSOEv8BKC3IpF9yEsu2KPGLiEjvpMR/AJKTjJL8DFZv\n3Rl1KCIiIp2ixH+ADhmQyeoqJX4REemdlPgP0LCCTFZv3Ym7Rx2KiIjIAVPiP0Al+ZnU1DdQo0v6\nRESkF1LiP0ADc9IA2FxTD8D89dX88/W1NDWpBUBERHo+3bnvAMUm/iG5GVx252wqq+up2d3AxceV\nRRuciIhIO+J2xG9mfzCzTWY2r1n5F81soZm9Y2b/G1N+jZktNbNFZnZmTPlEM3s7HPdrM7N4xdwR\nsYn/qYWVVFYHR/6PvL0hyrBEREQ6JJ5N/XcCZ8UWmNk0YDpwhLuPB/5fWD4OuAAYH85zi5klh7Pd\nCnwGGBW+3rPM7jYw+93E/8Q7leRnpvK5qSOYvaKKbXV7ogxNRESkXXFL/O7+HND8iTafA37q7vXh\nNJvC8unATHevd/cVwFJgspkNBvq7+ywPutH/CTgvXjF3RF5mKqnJxotLt/DwW+v50FElvH/CIJoc\nnl64qf0FiIiIRKi7O/eNBk4ys1fM7FkzOyYsHwqsiZlubVg2NBxuXh4ZM+OQgkyeCpP8VdNGcNjQ\nXAb1T+eJ+Xpkr4iI9Gzd3bkvBSgApgDHAPea2fCuWriZXQlcCVBcXExFRUVXLZra2tr9yxuQvJtl\nwMi8JN6e8zIA43IbeGZhJf95+hlSkiLthnDAYuuWSBK1XpC4dUvUeoHq1hslar26O/GvBe4Pm+1n\nm1kTUAisA0pjpisJy9aFw83LW+TuM4AZAJMmTfKpU6d2WeAVFRXsW95zNfOZvXEFV5wyganHBGFX\nZq3m6X+8zZgjj6W0ILPL1tsdYuuWSBK1XpC4dUvUeoHq1hslar26u6n/AWAagJmNBvoBW4CHgAvM\nLM3Mygk68c129w1AtZlNCXvzXww82M0x/5fPTxvB988Zx0cmvrtPUpIfJPs123aypbaeC2a8zLx1\nO6IKUUREpEXxvJzvb8DLwBgzW2tmlwN/AIaHl/jNBC7xwDvAvcB84DHgKndvDBf1eeB2gg5/y4BH\n4xVzRxVmp3H5ieUkxzTpl+RnALB22y4efGM9s5ZX8aWZr7NrT2NrixEREel2cWvqd/cLWxn1qVam\nvwG4oYXyOcCELgwtLgbnZmAGa6t2MmtFFQOy+rF8cx03PbaQH35wfNThiYiIALplb5fpl5LEkNwM\nFmys4Y3V2/noxBLOn1TCX19Zzd7GpqjDExERAZT4u9To4myenF/JnsYmjh1ewHEjBrCnsYlVW+ui\nDk1ERARQ4u9So4tz9g8fUZK3//2ijbVRhSQiIvIeSvxdaPzQ3P3DA7LTGDEwmySDRZU1EUYlIiLy\nLiX+LnTGuGIABuemA5CemkxZYRaLNyrxi4hIz6DH8nah9NRknvjq+8hOe3ezji7KYbGO+EVEpIfQ\nEX8XG12cw5C8jJj32azcWkd9g67nFxGR6Cnxx1lJfiZNDpuq66MORURERIk/3or6pwFQWb074khE\nRESU+ONuUNjRb6MSv4iI9ABK/HE2qH+Q+CvV1C8iIj2AEn+c5WakkpaSpKZ+ERHpEZT448zMGJqf\nweqtO6MORURERIm/O5QPyGKl7tcvIiI9gBJ/NygvDBJ/U5NHHYqIiPRxSvzdoKwwi917m9SzX0RE\nIqfE3w3KC7MAWLlFzf0iIhItJf5uUBYm/hXhef6FG6tpVLO/iIhEQIm/Gwzun05aShIrt9Tx55dX\nctYvn+eul1ZGHZaIiPRBSvzdICnJKBuQxYotddwzZw0ALyzdEnFUIiLSF+mxvN2krDCT11ZvZ0tt\ncAe/2SuqqNm9l5z01IgjExGRvkRH/N2krDCLzTX1uMMPzhlHbX0DH731ZerqG6IOTURE+hAl/m5S\nPiBr//Clx5dxw4cmsKiyhtufXxFhVCIi0tco8XeTE0YWctjQXH73qYkkJRmfPHYYk8sLePydjVGH\nJiIifYgSfzcpLcjkX188kbMmDNpfdvLogczfUM3WWj25T0REuocSf4SOKs0D4J311RFHIiIifYUS\nf4TGDMoBYNHGmogjERGRvkKJP0IDstMYmJPGgo064hcRke4Rt8RvZn8ws01mNq+FcV83Mzezwpiy\na8xsqZktMrMzY8onmtnb4bhfm5nFK+YoHFGSx2urtkUdhoiI9BHxPOK/EzireaGZlQJnAKtjysYB\nFwDjw3luMbPkcPStwGeAUeHrv5bZm00ZXsDKrTvZsGNX1KGIiEgfELfE7+7PAVUtjPoF8E0g9ik1\n04GZ7l7v7iuApcBkMxsM9Hf3We7uwJ+A8+IVcxQmlRUA8Oaa7RFHIiIifUG33rLXzKYD69z9zWYt\n9kOBWTHv14Zle8Ph5uWtLf9K4EqA4uJiKioquiZwoLa2tkuXt8+eRseAB154i18+8gZnlaVyZFH3\n3kk5XnWLWqLWCxK3bolaL1DdeqNErVe3ZRgzywS+Q9DMHxfuPgOYATBp0iSfOnVqly27oqKCrlxe\nrJFvPstjK2sBWFhVz5vXnkxuRvfdwz+edYtSotYLErduiVovUN16o0StV3f26h8BlANvmtlKoAR4\nzcwGAeuA0phpS8KydeFw8/KEMqks/z3v73xxZTSBiIhIwuu2xO/ub7t7kbuXuXsZQbP90e6+EXgI\nuMDM0sysnKAT32x33wBUm9mUsDf/xcCD3RVzdzl5dBEAHzhsMKePK+b255ezdFNtxFGJiEgiiufl\nfH8DXgbGmNlaM7u8tWnd/R3gXmA+8Bhwlbs3hqM/D9xO0OFvGfBovGKOymlji7j23HH85COH8b0P\njMWBXz21JOqwREQkAcXtHL+7X9jO+LJm728AbmhhujnAhC4NrodJSU7ishPKAeifnsrUMQOZs7Kl\nCyJEREQOju7c1wNNGpbPhh27dW2/iIh0OSX+HujQwf0BWLapLuJIREQk0Sjx90DDC7MAWL5FHfxE\nRKRrKfH3QANz0sjql8zyzTriFxGRrqXE3wOZGeUDs1ixRYlfRES6lhJ/D1VemK3ELyIiXU6Jv4cq\nL8xi7bad1Dc0tj+xiIhIBynx91DDC7NoclhTtTPqUEREJIEo8fdQ5ft69quDn4iIdCEl/h6qLEz8\nOs8vIiJdSYm/h8rNSKUwu9/+xF9X30BTk0cclYiI9HZK/D1YeWEWy7fUsXtvIyf/7Bmm//ZFPbVP\nREQOihJ/D1ZeGFzLP3fVNrbU7mFxZQ3n3/ayOvyJiEinKfH3YOWF2Wyuqeefr68jNdm4//PHs7ex\niRv+vSDq0EREpJdS4u/B9vXsv2/uWj5w2GDGD8nlnMMH8+KyLTTqfL+IiHSCEn8PduignP3Dl55Q\nDsCU4QOo2d3Agg3VUYUlIiK9WErUAUjrygqz+Pb7D2X55lqOLM0D4PCS4O/89dVMGJobZXgiItIL\nKfH3cJ89ecR73h9SkEl6ahILN9ZEFJGIiPRmaurvZZKTjDHFOSyqVFO/iIgcOCX+XmjMoBwWbtAR\nv4iIHDgl/l5ozKD+bK3bw+aa+qhDERGRXkaJvxcaG/b2X7Sxhurde3ls3gYaGpsijkpERHoDJf5e\naEyY+BdurOZ/H1vIZ+9+jW/e91bEUYmISG+gxN8LDchOozA7jVdXVnH/a+sAuP/1dTz+zsaIIxMR\nkZ5Oib+XGjs4h8ffqWTnnkYeuOoExhTncPMTi3DXHf1ERKR1Svy91KRhBQAMyOrHkaV5XHz8MBZX\n1ur6fhERaZMSfy/1oaOGAnD+MaUATByWD8CyzXpsr4iItC5uid/M/mBmm8xsXkzZz8xsoZm9ZWb/\nNLO8mHHXmNlSM1tkZmfGlE80s7fDcb82M4tXzL3JIQMyef6b07j6jDEADM7NAGDjjt1RhiUiIj1c\nPI/47wTOalb2JDDB3Q8HFgPXAJjZOOACYHw4zy1mlhzOcyvwGWBU+Gq+zD6rtCCT5KRgP6h/egqZ\n/ZJZv12JX0REWhe3xO/uzwFVzcqecPeG8O0soCQcng7MdPd6d18BLAUmm9lgoL+7z/Kg19qfgPPi\nFXNvZmYMzk1nY/WuqEMREZEeLMqH9HwauCccHkqwI7DP2rBsbzjcvLxFZnYlcCVAcXExFRUVXRZs\nbW1tly4vHtKbdrFw9a4DjrM31K0zErVekLh1S9R6gerWGyVqvSJJ/Gb2XaAB+EtXLtfdZwAzACZN\nmuRTp07tsmVXVFTQlcuLhwc2vs5rq7cfcJy9oW6dkaj1gsStW6LWC1S33ihR69Xtid/MLgXOAU71\ndy86XweUxkxWEpat493TAbHl0oL8rH5sq9sTdRgiItKDdevlfGZ2FvBN4IPuvjNm1EPABWaWZmbl\nBJ34Zrv7BqDazKaEvfkvBh7szph7k4LMftTUN7CnQfftFxGRlsXzcr6/AS8DY8xsrZldDvwGyAGe\nNLM3zOx3AO7+DnAvMB94DLjK3RvDRX0euJ2gw98y4NF4xdzbFWT3A2DbTh31i4hIy+LW1O/uF7ZQ\nfEcb098A3NBC+RxgQheGlrAKMoPEX1W3h+L+6RFHIyIiPZHu3JdA8rPCI36d5xcRkVYo8SeQgjDx\nb1XiFxGRVijxJ5B9iV/n+EVEpDVK/AkkLyMVCM7xi4iItESJP4GkJCeRm5Gqc/wiItIqJf4EU5DV\nT+f4RUSkVQeU+M0sycz6xysYOXgFWf10jl9ERFrVbuI3s7+aWX8zywLmAfPN7BvxD006Iz+zH1V1\ne6MOQ0REeqiOHPGPc/dqgsfhPgqUAxfFNSrptAFZ/dhSWw+Au7M1HBYREYGOJf5UM0slSPwPufte\nwNuZRyIyKDedLbX17Glo4on5lUy8/j88u3hz1GGJiEgP0ZHEfxuwEsgCnjOzYUB1PIOSzhual4E7\nVFbv5pmFmwD47dNLI45KRER6inYTv7v/2t2HuvvZHlgFTOuG2KQThuRlALB++y7mrd8BwOyVVSzd\nVBNlWCIi0kN0pHNfsZndYWaPhu/HAZfEPTLplMF5wcN5VlftZPHGWo46JA+Aj982i/qGxrZmFRGR\nPqAjTf13Ao8DQ8L3i4GvxCsgOThD8zIwg2cXb2ZPYxOfOnYYPzx3HFvr9vDYvI1RhyciIhHrSOIv\ndPd7gSYAd28AdOjYQ6WnJlM2IIuH39oAwOjiHC4+rozcjFReXrY14uhERCRqHUn8dWY2gLAnv5lN\nAXbENSo5KIcOygGgX0oSI4uySUoyJpcX8MqKqogjExGRqHUk8X8NeAgYYWYvAn8CvhjXqOSgTC4v\nAGDSsHwy+iUDcNjQXFZurWPXHjXWiIj0ZSntTeDur5nZycAYwIBF4bX80kNdclwZyUnGiSML95eN\nKsrGHZZtrmXC0NwIoxMRkSi1mvjN7MOtjBptZrj7/XGKSQ5SUpJx8XFl7ykbVZwNwNJNSvwiIn1Z\nW0f857YxzgEl/l5k2IAsUpKMJbqeX0SkT2s18bv7Zd0ZiMRXanISZYVZLKmsjToUERGJUEdu4HOj\nmeXFvM83s+vjG5bEw6iibJZuUuIXEenLOtKr//3uvn3fG3ffBpwdv5AkXkYVZbNya53u4Cci0od1\nJPEnm1navjdmlgGktTG99FAji3NoclixpS7qUEREJCIdSfx/AZ4ys8vN7HLgSeCu+IYl8TCqKOjZ\nr/P8IiJ9V0eu47/JzN4CTg2LrnP3x+MblsRDeWEWyUnGG2u2M2dlFTnpqXz9jNFRhyUiIt2o3cQP\n4O6PAo/GORaJs/TUZE4ePZA7XlixvywvM5WREcYkIiLdq9WmfjN7IfxbY2bVMa8aM6tub8Fm9gcz\n22Rm82LKCszsSTNbEv7Njxl3jZktNbNFZnZmTPlEM3s7HPdrM7POV1e+ctoojhs+gOvOm8Cphxbx\ny/8sYXeDRx2WiIh0k1YTv7ufGP7Ncff+Ma8cd+/fgWXfCZzVrOzbwFPuPgp4KnyPmY0DLgDGh/Pc\nYmbJ4Ty3Ap8BRoWv5suUA3B4SR5/u3IKF00ZxuenjaS2voFZGxqiDktERLpJR67j/3NHyppz9+eA\n5o+Dm867HQPvAs6LKZ/p7vXuvgJYCkw2s8FAf3ef5e5O8ICg85AucfQheYwpzuG5tUr8IiJ9RUd6\n9Y+PfWNmKcDETq6v2N03hMMbgeJweCiwJma6tWHZ0HC4ebl0ATNj6piBrK5uItivEhGRRNfWQ3qu\nAb4DZMSc0zdgDzDjYFfs7m5mXZptzOxK4EqA4uJiKioqumzZtbW1Xbq8nqJ2014aHP71ZAX9+yVW\n94lE/cwgceuWqPUC1a03StR6tXWv/p+Y2U3A7e7+6S5aX6WZDXb3DWEz/qawfB1QGjNdSVi2Lhxu\nXt5azDOzdnShAAAgAElEQVQId0omTZrkU6dO7aKwoaKigq5cXk+xe94G/rLwNYaPPzrhntqXqJ8Z\nJG7dErVeoLr1Rolarzab+t29CTimC9f3EHBJOHwJ8GBM+QVmlmZm5QSd+GaHpwWqzWxK2Jv/4ph5\npAsMys0AoLJ6d8SRiIhId+jIOf7XzOyAk7+Z/Q14GRhjZmvDu/79FDjdzJYAp4Xvcfd3gHuB+cBj\nwFXuvu+G8p8Hbifo8LcM3U+gSw3OTQdgww4lfhGRvqAjN/A5Fvikma0C6gjO87u7H97WTO5+YSuj\nTm2p0N1vAG5ooXwOMKEDcUonFGankWSwUYlfRKRP6EjiP7P9SaS3Sk4y8tJMR/wiIn1ER+7VvwrA\nzIqA9LhHJN0uP83YWL0r6jBERKQbdOQGPh8Mz8mvAJ4FVqLz7AklP93U1C8i0kd0pHPfdcAUYLG7\nlxOco58V16ikWxWkB039+27i8+zizdz02EJq63VHPxGRRNORxL/X3bcCSWaW5O7PAJPiHJd0o/z0\nJHbuaaQmTPS3PLOUWyuW8dFbX2LHrr0RRyciIl2pI4l/u5llA88BfzGzXxH07pcEUZAe3LFvX3P/\nvmv6F26s4Tv3vx1ZXCIi0vU6kvinAzuBrxJcY78MODeeQUn3yg8T/4Ydu2lqctZv382Fkw9h6piB\n/PvtDTrqFxFJIG0mfjM7D/gccLq7N7j7Xe7+67DpXxJEflqQ+Ct37GZLXT17GpsYOziHS44rA2Bx\nZU2E0YmISFdqNfGb2S0ER/kDgOvM7PvdFpV0q31H/Ot37GLV1p0AlORnMGZQDhA0+YuISGJo64j/\nfcAp7n4NMBU4r1sikm6XkmQU5aSxbtsuXl62FTM4sjSfwbnp5GWm8vba7funffCNdZz5i+d4cemW\nCCMWEZHOaivx79l3v3x330lwq15JUKUFmazZtpPnl2zmsKG5FGT1w8w4pqyAV1ZUAbCtbg/fe2Ae\niyprePit9RFHLCIindHWnfsONbO3wmEDRoTvO3SvfuldDinI5KkFldTtaeSzJw/fX37c8AE8Ob+S\nZZtruenRhdTVN1CUk8bb63ZEGK2IiHRWW4l/bLdFIZErzc+gendwHf9JowbuLz99XDE/fng+F93+\nCut37ObqM0ZTW9/IHS8sp76hkbSU5KhCFhGRTmi1qd/dV7X16s4gJf6mDB9Av5QkTjm0iMllBfvL\nSwsyGT+kP+t37OaYsnyumjaSw0ty2dvoLFKnPxGRXqcj1/FLH3D8yELeuvYMZlw0kaSk93bnuOKk\ncgC+fOpozIzDhuYCqLlfRKQX6shjeaWPSE9tudn+Q0eVcNzwQgblBg9nLMnPIDcjlXlK/CIivU5b\n1/E/Ff69qfvCkZ5qX9IH9h/1v7VWiV9EpLdpq6l/sJkdD3zQzI4ys6NjX90VoPRMh5XksriyhvqG\nxqhDERGRA9BWU/8PgO8DJcDPm41z4JR4BSU932FD3+3gd3hJXtThiIhIB7Wa+N39PuA+M/u+u1/X\njTFJL7Cvg99ba3co8YuI9CLtdu5z9+vM7IMEt/AFqHD3h+MblvR0JfkZ5GW+t4PfpurdrNhSx7HD\nB0QYmYiItKXdy/nM7CfAl4H54evLZnZjvAOTns3MGFOcw5JNtQCs3baTyTc+xcdnzOKd9er0JyLS\nU3XkOv4PEDyW9w/u/gfgLOCc+IYlvcGIomyWbqrF3fntM8v2l3/3n/NoaGyKMDIREWlNR2/gE3sS\nNzcegUjvM3JgNjt27WVL7R5mr9jKKYcW8fPzj+CNNdt5cn5l1OGJiEgLOpL4fwK8bmZ3mtldwFzg\nhviGJb3B+CH9geBRvcs21zFxWD7TjxzK0LwM7pu7NuLoRESkJR3p3Pc3M6sAjgmLvuXuG+MalfQK\nk8oKGJKbzvX/XkCSwQcOG0xyknHq2CLum7uWPQ1N9EvRXaFFRHqSDv0qu/sGd38ofCnpCwDJScat\nn5rIxyaW8KPpEygrzAKCp/vt3NPIa6u3RRyhiIg0p3v1y0E5ojSPI0rfex3/lOEFpCQZzy/ZzJTh\nA3h6YSXPLd7CF08ZyYDstIgiFRERiCjxm9lXgSsI7gD4NnAZkAncA5QBK4Hz3X1bOP01wOVAI/Al\nd3+8+6OWjspJT+WYsgLunbOWwbkZfO+BeQAsrqzhL1cci5m1swQREYmXNpv6zSzZzBZ25QrNbCjw\nJWCSu08AkoELgG8DT7n7KOCp8D1mNi4cP57gUsJbzKzlx8hJj/GDc8dRv7eR7z0wj0H90/nKaaN4\nadlW3llfHXVoIiJ9WpuJ390bgUVmdkgXrzcFyDCzFIIj/fXAdOCucPxdwHnh8HRgprvXu/sKYCkw\nuYvjkS42dnB/nrl6Kn++fDJPfO19fGJy8BV6fsmWiCMTEenbzN3bnsDsOeAoYDZQt6/c3T/Y6ZWa\nfZngksBdwBPu/kkz2+7ueeF4A7a5e56Z/QaY5e53h+PuAB4NnyXQfLlXAlcCFBcXT5w5c2ZnQ/wv\ntbW1ZGdnd9nyepLuqtt3X9hJQXoSX5+U3v7EXUCfWe+TqPUC1a036m31mjZt2lx3n9TedB05x//9\nLohnPzPLJziKLwe2A383s0/FTuPubmZt75G0wN1nADMAJk2a5FOnTj34gEMVFRV05fJ6ku6q2/Gb\n3+TphZs4+eSTu+U8vz6z3idR6wWqW2+UqPVq93I+d3+WoLNdajj8KvDaQazzNGCFu292973A/cDx\nQKWZDQYI/24Kp18HlMbMXxKWSS9zWEkuW+v2sH7H7qhDERHpszrykJ7PAPcBt4VFQ4EHDmKdq4Ep\nZpYZNumfCiwAHgIuCae5BHgwHH4IuMDM0sysHBhFcNpBeplxg4M7/S3aqA5+IiJR6UhT/1UEnele\nAXD3JWZW1NkVuvsrZnYfQatBA/A6QfN8NnCvmV0OrALOD6d/x8zuJXgyYANwVdjpUHqZEQODc2XL\nNtVxyqERByMi0kd1JPHXu/uefedkw574B3z+PZa7Xwtc23w9BEf/LU1/A3o+QK+Xn9WPAVn9WBo+\nyldERLpfR27Z+6yZfYfg8rvTgb8D/4pvWJKoRhRls2yzEr+ISFQ6kvi/DWwmuMPe/wCPAN+LZ1CS\nuEYMzGbp5lrau4xURETioyNP52sKH8f7CkET/yLXr7Z00siibLbv3EtV3R7dt19EJAId6dX/AWAZ\n8GvgN8BSM3t/vAOTxDRiYPAEP53nFxGJRkea+m8Gprn7VHc/GZgG/CK+YUmiGj8kF4DX12yPOBIR\nkb6pI4m/xt2XxrxfDtTEKR5JcANz0hhZlM2s5VujDkVEpE9q9Ry/mX04HJxjZo8A9xKc4/8Ywd37\nRDplyvAC/vnaOhoam0hJ7si+p4iIdJW2OvedGzNcCZwcDm8GMuIWkSS8KcMHcPes1cxbX02TO4cU\nZFKojn4iIt2i1cTv7pd1ZyDSdxw3fAD9kpP4yszXWbl1J2OKc3j8q++LOiwRkT6hI736y83s52Z2\nv5k9tO/VHcFJYhqQncZPP3IYNbsbAFhUWcM/X1/LQt3DX0Qk7jpyy94HgDsI7tbXFN9wpK/48NEl\nfPjoEqp37+W4G5/iq/e8CcAPzhlHkzsnjCwks18yGanJFPVPjzhaEZHE0ZHEv9vdfx33SKRP6p+e\nym8+eTQ3/nsB67fv4scPz3/P+JL8DF741ikRRScikng6kvh/ZWbXAk8QPEgHAHd/LW5RSZ8ybUwR\n08YED3xcUlnD2m27uLViGbNXVrF22y5eXVnFMWUFEUcpIpIYOpL4DwMuAk7h3aZ+D9+LdKlRxTmM\nKs5h2qFFbKrezYdueYlrH3yHR758UtShiYgkhI4k/o8Bw919T7yDEYlV1D+dTxx7CD97fBFVdXso\nyOoXdUgiIr1eR+6eMg/Ii3cgIi2ZMnwAAM8v2RxxJCIiiaEjR/x5wEIze5X3nuP/YNyiEgkdVZpH\nSX4GM2evYfqRQ6MOR0Sk1+tI4r827lGItCIpybhwctDcv3xzLcMHZkcdkohIr9ZuU7+7P9vSqzuC\nEwH42MQSAB6dtzHiSEREer+O3Lmvxsyqw9duM2s0M91iTbpNUf90RhZlM2dlFQDPLt7Mmb94jmWb\nayOOTESk9+nIEX+Ou/d39/4ED+f5CHBL3CMTiXFMWT5zVm7jnldXc8kfZrOosoa7XloZdVgiIr3O\nAT0T1QMPAGfGKR6RFn10Yik19Q186x9vM3Zwf045tIg/vbyKMd97lLmrtkUdnohIr9Fu5z4z+3DM\n2yRgErA7bhGJtGDisHyuO28C89dX86VTR5Kf2Y/zb3uZt9bu4LN3z+WmjxzGtDFFmFnUoYqI9Ggd\n6dV/bsxwA7ASmB6XaETacNGUYe95P/PKKTy3eDM/fGg+n75zDpPLCzhtbBGXHl9Ov5QDaswSEekz\n2k387n5ZdwQicqAy+6Vw1oTBnDq2mLteWskfX1zJjY8s5MZHFnLiyEI+XKKHSYqINNdq4jezH7Qx\nn7v7dXGIR+SApSYnccVJw7nipOH89pmlzFlZxUvLtrKz2vjwWVFHJyLSs7R1xF/XQlkWcDkwAFDi\nlx7nqmkjAfi/p5Zw85OLmbOyikl6sp+IyH6tngh195v3vYAZBJfyXQbMBIYfzErNLM/M7jOzhWa2\nwMyOM7MCM3vSzJaEf/Njpr/GzJaa2SIz0xUF0q4rThpOXprxv48vwt2jDkdEpMdoswdUmIyvB94i\naB042t2/5e6bDnK9vwIec/dDgSOABcC3gafcfRTwVPgeMxsHXACMB84CbjGz5INcvyS4jH7JnDsi\nldkrqvjds8tZvXUnH7/tZeauqoo6NBGRSLWa+M3sZ8CrQA1wmLv/0N0P+oJpM8sF3gfcAeDue9x9\nO8GVAneFk90FnBcOTwdmunu9u68AlgKTDzYOSXynlKbwvtED+b+nl3D5Xa/yyooqLvz9K7yyfGvU\noYmIRMZaawY1syaCp/E1ALETGUHnvv6dWqHZkQSnDuYTHO3PBb4MrHP3vHAaA7a5e56Z/QaY5e53\nh+PuAB519/taWPaVwJUAxcXFE2fOnNmZEFtUW1tLdnZiPiAmUetWW1vL8l3p/Hxu8FDJ80am8tza\nBpIMvntsOvnpvfeSv0T+zBKxXqC69Ua9rV7Tpk2b6+6T2puu1c597h6vX8UU4Gjgi+7+ipn9irBZ\nP2bdbmYHfGLW3WcQ7FQwadIknzp1aheEG6ioqKArl9eTJGrdKioq+NI5U5k8cSvpqckcWZrH66u3\n8cnbX+FHsxv40QfH99pH/SbyZ5aI9QLVrTdK1HpFccizFljr7q+E7+8j2BGoNLPBAOHfff0I1gGl\nMfOXhGUiHTJl+ACOLM0D4KhD8nnoCydQNiCLL898gwde11dJRPqWbk/87r4RWGNmY8KiUwma/R8C\nLgnLLgEeDIcfAi4wszQzKwdGAbO7MWRJMCOLcrjvs8dxZGkeNz22kIZG3ehHRPqOqE5yfhH4i5m9\nBRwJ3Aj8FDjdzJYAp4Xvcfd3gHsJdg4eA65y98ZIopaEkZKcxBemjWTDjt385NGFANTVN+jSPxFJ\neB25V3+Xc/c3CB7209yprUx/A3BDXIOSPufUsUVcenwZd7ywgk019Tw2bwMXH1fG988ZF3VoIiJx\n03u7NYscJDPjux8YyyePPYTH5m1gb6NzxwsrWLChOurQRETiRolf+rTU5CRu+NBhvP3DM5l1zakU\nZqfx43/NjzosEZG4UeIXAdJTkxmUm87lJ5bz8vKtLNtcG3VIIiJxocQvEmP6kUMAeHJ+ZcSRiIjE\nhxK/SIwheRmMG9yfpxcc7OMoRER6JiV+kWZOG1vEnFVVbKvbE3UoIiJdTolfpJlTxxbT5FCxWEf9\nIpJ4lPhFmjlsaC5FOWn8fc5a5qys4rfPLGWv7u4nIgkikhv4iPRkSUnGp6YM4+dPLualZS8DsKW2\nnmvPHR9xZCIiB0+JX6QFXzp1FMePGMB/FmzijTXbuPOllQD8z/tGMCg3PdrgREQOghK/SCsmlRUw\nqayAXXsaueGR+fzxxZX86eVVfPioofzsY0dEHZ6ISKfoHL9IOzL6JXP9eYfxzNVTGTEwi7/PXcua\nqp1RhyUi0ilK/CIdVF6YxZ2XTSY5yZjx3PKowxER6RQ19YscgCF5GXxi8iH8bfZqLj2hjA3bd3P3\nrFWkJBs3nHcYuZmpUYcoItImJX6RA/Tl00bx4BvrOO+3L1Jb30BuRirbd+5l7OD+XDVtZNThiYi0\nSU39IgeoMDuNP19+LEPzMphSPoBnr57GlOEF/H3OGtw96vBERNqkxC/SCUeU5vHYV97H366cQm5m\nKh+dWMrKrTuZvaIq6tBERNqkxC/SBc4+bBCF2Wn85NGFNDXpqF9Eei4lfpEukNkvhe+cfShvrNnO\ng2+uizocEZFWKfGLdJEPHTWU8sIs7nppFZtr6qmrb4g6JBGR/6LEL9JFzIzzjhzKG2u2c8wN/+Ho\n657k73PWsHRTjZr/RaTH0OV8Il3oqmkjGFmUTVVdPX95ZTXfuO8tAE4cWcjvLppIdpr+5UQkWvoV\nEulCKclJfODwwQB8ZGIJs5ZvZeHGGm5+YjGn3fwsnzj2EKYMH8DYwTnkpOtmPyLS/ZT4ReIks18K\npxxazCmHFnNMWQE/e3wRP39yMQD5malMG1PEqOIcPnz0UIr764l/ItI9lPhFusExZQXcc+UU1m7b\nxVtrd/DP19fy3JIt3P/6Ol5bvY0ZF03EzKIOU0T6ACV+kW5iZpQWZFJakMkHDh+Mu3PzE4v5zTNL\n+citL5Gf2Y+UZON7HxhHaUFm1OGKSIJS4heJiJnxtdNHU5DVjz/PWsW89dXsbWzi2cWbKcpJZ2RR\nNkeV5nHx8WXkZqg/gIh0jcgSv5klA3OAde5+jpkVAPcAZcBK4Hx33xZOew1wOdAIfMndH48kaJEu\nlpRkfPrEcj59YjnuzvItddz10kpWbKljyaYanl64idkrq/jz5cdGHaqIJIgoj/i/DCwA+ofvvw08\n5e4/NbNvh++/ZWbjgAuA8cAQ4D9mNtrdG6MIWiRezIwRA7P58fQJ+8tmPLeMGx9ZyBV3vUpyklFX\n38gdl06KMEoR6e0iSfxmVgJ8ALgB+FpYPB2YGg7fBVQA3wrLZ7p7PbDCzJYCk4GXuzFkkUhccnwZ\nD7+1gf8s2LS/7H3/+wz19XsomFPBdedNICsthbtnreKCY0qZVFYQYbQi0htEdcT/S+CbQE5MWbG7\nbwiHNwLF4fBQYFbMdGvDMpGEl5aSzN8/exwbtu9mT2MTr6/eRsWizVRt2cza3U1c/fc3KS/M4qVl\nW3nozfV88thDuGjKMIYPzI46dBHpoay7nx9uZucAZ7v7581sKnB1eI5/u7vnxUy3zd3zzew3wCx3\nvzssvwN41N3va2HZVwJXAhQXF0+cOXNml8VdW1tLdnZi/pgmat0StV4Q1O2t6nRmvFUPwGmHpFC7\n15mzsZG0FPjhcRkMzOx9d+RO9M9Mdetdelu9pk2bNtfd2z0XGMUR/wnAB83sbCAd6G9mdwOVZjbY\n3TeY2WBgX9vmOqA0Zv6SsOy/uPsMYAbApEmTfOrUqV0WdEVFBV25vJ4kUeuWqPWCoG5fPeMk3qyZ\nzfade/nFp48jJz2V1Vt3cvovnuXVnQP437OPiDrMA5bon5nq1rskar26/ZDA3a9x9xJ3LyPotPe0\nu38KeAi4JJzsEuDBcPgh4AIzSzOzcmAUMLubwxbpcTL6JXPP/xzHY185af/tfw8ZkMn0I4fwyNsb\n2dvYFHGEItIT9aS2wJ8Cp5vZEuC08D3u/g5wLzAfeAy4Sj36Rd7V/I5/p44tpra+gVdXVkUUkYj0\nZJHewMfdKwh67+PuW4FTW5nuBoIrAESkHSeMLCQ12ahYtJnjRxRGHY6I9DA96YhfRLpAdloKx5YP\n4JmFm3ht9Ta+du8bzFq+le0790Qdmoj0ALplr0gCmjpmINf/ewEfvuUlAO5/bR2pycZ3zx7LpSeU\nRxydiERJR/wiCeiUQ4sAyEhN5oGrTuD3F09iyvAB3PDIAtZU7Yw4OhGJkhK/SAIaPjCbv3/2OGZ9\n51SOLM3j9HHF/L+PHUGSGb94cjEAd89axQk/fZoVW+oijlZEupOa+kUS1DHNbt9b3D+dS08o47Zn\nl/P6mu37E/43/v4mXztjNMcNH/BfVwiISOLREb9IH/K5k0cwqH86K7fWccWJ5fzgnHHMXb2NT/z+\nFf744sqowxORbqAjfpE+JC+zH49/5X1sqN7FoYOCB2N+8MghfPGvr/OL/yzm48eUkpWmnwWRRKYj\nfpE+JjczdX/SByjMTuPqM8dQs7uBv81eDUB9Q6Pu/CeSoLRrLyIcfUgeJ40q5OYnFjNr+VZeWraV\nw4bmcvcVx5KarOMDkUSi/2gRwcz42UePYGRRNvPWVTOprIBXVlQx6ruP8oW/vkZTU/c+xVNE4kdH\n/CICwKDcdP71xRMBcHcef2cjj79TyT9fX8eehiZ++8mjSU1Owt3V+1+kF1PiF5H/YmacNWEwZ4wb\nREFWP+54YQX/9/RSUpKM3z+/nLMnDObUsUWcMX5Q1KGKyAFS4heRViUlGd8/Zxwbq3fz66eWADB+\nSH/umbOGe+as4f8uPIpzjxgScZQiciCU+EWkXTd95HAOH5rLyKJsTh1bzJqqnVzyx9ncUrGMcw4f\nrKZ/kV5EnftEpF3ZaSn8z8kjOHVsMQClBZlcedJwFmyoZvaKqoijE5EDocQvIp0y/cihFGT14+Yn\nFuOuXv8ivYUSv4h0Ska/ZL56+mhmr6zi5WVbeXJ+Jbv2NGonQKSH0zl+Eem0jxw9lJ88soBP3P7K\n/rLSggw+fUI5H5lYQv/01AijE5GW6IhfRDots18K1547jvzMVC6aMoyxg/uzpmoXP/rXfD79x1dp\n1I1/RHocHfGLyEH5+DGHcP6k0v09+xubnL/NXs33HpjHUwsqOWP8INZv38XuvY0MH5gdcbQiosQv\nIgct9nK+5CTjgmNKubViGT97fBGPvL2BB95YD0BuRiqTyws45dAizp9USnKSLgMU6W5K/CLS5VKS\nk/juB8by+b+8xpJNtVx6fBnZaSms3baTZxZt5sn5lWypqeeLp46KOlSRPkeJX0Ti4uzDBvOfr50M\nwMiid5v43Z3P3f0aNz+5mFkrtrK5pp4ffXACx40YEFWoIn2KEr+IxE1swt/HzPjCKSP5z4JK3llf\nTe3uBi78/SymHzmEMamNnKyHAInElRK/iHS7CUNzeeuHZ5CRmsziylp+88xSnpxfyYN7GtmRsZBr\nzh4bdYgiCUuJX0Qikdkv+PkZMyiH/7vwKLbW1vOFO57htueWc/Logextcp5ZuIkzxw/i6GF5pKUk\nRxyxSGJQ4heRHmFAdhqXTkhj417jU3e8wr5bANz50kpy0lK4+PhhXH3GGJ0GEDlIuoGPiPQYacnG\n7z41kfLCLCaXF/DyNafw/z52BO8bPZDfPrOMH/1rftQhivR63X7Eb2alwJ+AYsCBGe7+KzMrAO4B\nyoCVwPnuvi2c5xrgcqAR+JK7P97dcYtI9xgzKIenvj51//uPTizhI0cPZdC/07njhRVMGV7AWRMG\nRxegSC8XxRF/A/B1dx8HTAGuMrNxwLeBp9x9FPBU+J5w3AXAeOAs4BYz08k+kT7EzLjm/Ycyqiib\nb/3jbT72u5f47TNL9UAgkU7o9sTv7hvc/bVwuAZYAAwFpgN3hZPdBZwXDk8HZrp7vbuvAJYCk7s3\nahGJWkpyEjd99HBGFmXzzvpqfvb4Ii7946u8urKKeet20NDYFHWIIr2CRbnHbGZlwHPABGC1u+eF\n5QZsc/c8M/sNMMvd7w7H3QE86u73tbC8K4ErAYqLiyfOnDmzy2Ktra0lOzsx7zOeqHVL1HpB4tat\no/Vqcufp1Q3cu3gPexqDsuG5SXxsdD92NzrD+idRkN6zujAl6mcGiVu33lavadOmzXX3Se1NF1mv\nfjPLBv4BfMXdq2N76rq7m9kB75G4+wxgBsCkSZN86tSpXRQtVFRU0JXL60kStW6JWi9I3LodSL1O\nAa7asYuFG2pYtrmW6/+9gJte3Q1AemoSpx5azNjBOVw4+RAGZKfFL+gOStTPDBK3bolar0gSv5ml\nEiT9v7j7/WFxpZkNdvcNZjYY2BSWrwNKY2YvCctEpI8bnJvB4NwMph1axBGleWypqWdAdhr/fH0t\nzyzczL/f3sBdL6/i5+cfwUmjBkYdrkiPEEWvfgPuABa4+89jRj0EXAL8NPz7YEz5X83s58AQYBQw\nu/siFpHe4Jiygv3Dk8uD4QUbqrnsj69y0R2zmVxWwNptOxman0FR/3QmlxWQl5nK+CH9GVmUE1XY\nIt0uiiP+E4CLgLfN7I2w7DsECf9eM7scWAWcD+Du75jZvcB8gisCrnL3xu4PW0R6m7GD+1Pxjalc\n/+/5/7+9ew+Oq7zPOP797Wq1q7ssWZaE5YssYyxjEUwM2OBgE64mDIRJLwSaYRoIZCZJk2GaFiZt\nh6YNvdOm005SSkihbfAMl3IxyYDB3OIwYBtfseS7wDdJlmxJlmTJq9XbP/ZYXslYvkTW6px9PjM7\ne867Z4/fx17rp/Oes+flve2t1FYWsuaTw6xuPMyrGw8Mbje5OIeinAj50SzKCqJ8fdF0Pj+tZIQ9\ni/jXmBd+59yvgVPdeuu6U7znR8CPzlunRCSwYpEwf/3lusH1jp44rd195EezaOns4+2tLWxr6WJ/\n+1F64v2sbOigq6+fp76uLw9JMOmWvSKSUYpyIxTlRgAoL4xRV1U05PW/+VU9P3tvNy2dvUwqjKWj\niyLn1fj6vouISJrddcVUHPDYim0A7G8/yhtbmonrPgESEDriFxFJMa00j/sWVfMf7+5iSkkur2zY\nT0PTET43pZjnvrmQSFjHS+Jv+gSLiAzzxzddxIIZJfzDa1tpaDrCopkT2bCnnT9/cTMvrd/HsX4d\n/QxLPwMAAAzTSURBVIt/6YhfRGSYSDjEL+5bwIr6ZvoTjlvqKvjOM+tYtnoPy1bv4dHCeh64poZw\nyOiNJ1jdeBiAb3yhmitnlKa59yIjU+EXEfkMoZBx08UVg+v/dtdlPHJbH+s+becnb+/gh8tPTBEc\ni4TIiYT59Y6D/NXtc7nt0gt47eNm/mvVbuZOLuJPb57NgCYUknFChV9E5AxNzI9yw5xyrq+dxO7W\nbkJmZGeFqCyK0dzZx31Pr+b7z23k+89tHHzPR5+28/T7n1CWY/x4aisLZ5Syv6OXifnZRLM00aiM\nPRV+EZGzZGbMKBs6eUtFUYxnH7iKf3p9K8cSAyyeVcbiWWWs2NLM7rZunnp3G3f95wcU50Zo74mT\nFTIKYllcX1vO3QumUVOWx/s729i8r4P7F9eQH9WPZzk/9MkSERklOdlh/uzWOUPaltZVAlCT2EOD\nq2J7yxEun17CgY5etjcf4aX1+3l27d4h7/nXlTu4emYpN11cwR3zJlMQi4xZBgk+FX4RkTEQDRvf\nXXLhSe0dPXFe3XSAzt44tZWFtHT28tGnyVsKr9rRxlO/aeShpbVcN3sS7Ufj7G7tprwwStWE3DSk\nkCBQ4RcRSaOi3Ah3XTl1SNvvzp/Co3fU8dbWFh55eQvfeHoNE3IjdPb2kxhwhEPGLXWV1JTlcUlV\nETPLCijNzyZPpwfkDOhTIiIyDpkZX5xdzqKZZby6aT+rdrRRURhj3tRi3t12kOUbD/DKhv1D3pOX\nHaauqoiZk/KZWZbP3QumDd5wqDeeIBbRxYSiwi8iMq5lZ4W4Y14Vd8yrGmy7rracR267mL7+Ad6s\nb6GrL86Huw8TjYRY23iY5RsP0N4T55ebmqgqyWHPoR5WNx5mWmku86YUc+cVU6mtLKQo5+yvHXDO\nsb+jl5LcbHKy9YuEH6nwi4j4kJkRi4T50iXJiwd///Khpwsee30rP1/VyKZ9HUyekMP1tZP4YPch\nVmxp5sX1+zGD3EiYaCRM9cQ8YpEQS2ZNIjeaXO/oibO9pYvmzt7BffYnHL/Z1cqeQ0eJRULUVhYy\nMT/KtHCcuV19lOZlY2Y45zA71SSskm4q/CIiAfTgjRfx4I0XndTe2RvnvW2tbG85wrpP23l/VxuR\nsLHrYDerdtSftH2ymJ9Yv/iCIu69upqdB7vZ3drNx/s6WNFxjCc2vcEFRTGmlOSyYW87N8yp4KaL\ny/nChWXnNLIg548Kv4hIBimMRbxRguRIQV9/gmhWmMSAo627j774AI1t3eRmZzFzUv5pi/bAgOPJ\nl1fSP2E6axoP0dTZy8IZpbzV0MIrG/ZTGEvup7Ioh+LcE/uaVppLTnayBO093MOaxsM455hemsfc\nyUUsraugOCeb+qZOZpUX6L4Go0h/kyIiGez43QPDIWNSQQyAKSVn/lXBUMiYWRxmyeIaWFwz2N4b\nT/BmfQvPf7SXvv4E6/e009efACCecHQcjQ/ZT93kIvKiYd7dfpAX1u3jh8u3kBMJczSeoCgnwsNL\nZwNw8EgfV84opaIwRmVxjK7efqKREIe6jzG5OEenGM6ACr+IiIy649cfHL8GIZVzjsM9cRIDyfkL\nskLGhLzswdd3HeziV5ubaGztZv70Cfz0nV089MKmk/YTMhhwJ56nlOSwcEYpWeEQhbEId185lZKU\n/R7oOMqHuw8ztSSXBTNKGHAQTyRnWjzUfYzmzl7KC2OUDDu9ETQq/CIiMqbMbEhBHm5GWT7funbm\n4PpXLqtiV2u3NyoRZdWONtp7jrG7rZtoVpiBAUdZQZSVDS2sbDgIOFq7jvHTd3ae8s8oiGXRG08Q\nT5x68qS6iWH6ypoIm9FxNE5bdx/dfQlGmm4pPxpm8axJRLNOzHp/QXEO2Snr6abCLyIi41pWOMSs\n8oLB9ZvnVnzmdvdcNX1w+ZO2blY2tAwe0QPkRyMsrClla9MR3tl2kLzsMJMKowDkRMJMLc3j07Zu\njsYTtPfEeeK9nTzw32vPur+P/rJhyHpRToTqiXknjSLceskF3Luo+qz3/9tS4RcRkcCZVprHH179\n2UW1emLeKX95gLLBpdrQAUpr6ghb8lREXjSLsvzoiEfvDU2d1B/oHFyPJxwf7DpEy5Hek7aNpmkU\nQIVfRETkMxRkG1fVTDyr98yuKGR2ReGQtt+bP2U0u/VbGz8nHUREROS8U+EXERHJICr8IiIiGUSF\nX0REJIOo8IuIiGQQ3xR+M7vZzLaa2Q4zeyjd/REREfEjXxR+MwsD/w4sBeYAXzWzOentlYiIiP/4\novADVwA7nHO7nHPHgGXA7Wnuk4iIiO+YcyPddXh8MLPfAW52zt3nrX8NuNI59+1h290P3A9QXl7+\n+WXLlo1aH7q6usjPzx+1/Y0nQc0W1FwQ3GxBzQXK5kd+y3Xttdeudc7NP912gbpzn3PuceBxgPnz\n57slS5aM2r7ffvttRnN/40lQswU1FwQ3W1BzgbL5UVBz+WWofx+Qes/DKq9NREREzoJfCv9q4EIz\nqzazbOBO4OU090lERMR3fDHU75zrN7NvA68BYeBJ59zHae6WiIiI7/ji4r5zYWYHgU9GcZcTgdZR\n3N94EtRsQc0Fwc0W1FygbH7kt1zTnHNlp9sosIV/tJnZmjO5WtKPgpotqLkguNmCmguUzY+Cmssv\n5/hFRERkFKjwi4iIZBAV/jP3eLo7cB4FNVtQc0FwswU1FyibHwUyl87xi4iIZBAd8YuIiGQQFf7T\n8Pt0wGb2pJm1mNnmlLYSM1thZtu95wkprz3sZd1qZjelp9enZ2ZTzOwtM9tiZh+b2Xe99iBki5nZ\nh2a2wcv2l16777NBcrZNM1tnZsu99aDkajSzTWa23szWeG1ByVZsZs+ZWYOZ1ZvZwiBkM7OLvH+v\n449OM/teELKNyDmnxykeJG8WtBOYAWQDG4A56e7XWWa4BrgM2JzS9vfAQ97yQ8DfectzvIxRoNrL\nHk53hlPkqgQu85YLgG1e/4OQzYB8bzkCfAAsCEI2r78PAr8Algfl8+j1txGYOKwtKNmeAu7zlrOB\n4qBkS8kYBpqAaUHLNvyhI/6R+X46YOfcu8ChYc23k/yPjPf85ZT2Zc65PufcbmAHyb+Dccc5d8A5\n95G3fASoByYTjGzOOdflrUa8hyMA2cysCvgS8ERKs+9zjcD32cysiOQBxM8AnHPHnHPtBCDbMNcB\nO51znxC8bEOo8I9sMrAnZX2v1+Z35c65A95yE1DuLfsyr5lNB+aRPDIORDZvOHw90AKscM4FJdu/\nAH8CDKS0BSEXJH85e8PM1npThEMwslUDB4Gfe6donjCzPIKRLdWdwDPectCyDaHCn+FccvzKt1/t\nMLN84Hnge865ztTX/JzNOZdwzl1KcibKK8xs7rDXfZfNzG4FWpxza0+1jR9zpVjk/ZstBb5lZtek\nvujjbFkkTxf+xDk3D+gmOfw9yMfZAPAmf7sNeHb4a37P9llU+EcW1OmAm82sEsB7bvHafZXXzCIk\ni/7/Oude8JoDke04b0j1LeBm/J/tauA2M2skedrsi2b2P/g/FwDOuX3ecwvwfySHgIOQbS+w1xt1\nAniO5C8CQch23FLgI+dcs7cepGwnUeEfWVCnA34ZuMdbvgd4KaX9TjOLmlk1cCHwYRr6d1pmZiTP\nOdY75x5LeSkI2crMrNhbzgFuABrweTbn3MPOuSrn3HSS/5dWOuf+AJ/nAjCzPDMrOL4M3AhsJgDZ\nnHNNwB4zu8hrug7YQgCypfgqJ4b5IVjZTpbuqwvH+wO4heQV4zuBH6S7P+fQ/2eAA0Cc5G/u9wKl\nwJvAduANoCRl+x94WbcCS9Pd/xFyLSI5/LYRWO89bglItkuAdV62zcBfeO2+z5bS3yWcuKrf97lI\nfvNng/f4+PjPiiBk8/p6KbDG+0y+CEwIULY8oA0oSmkLRLZTPXTnPhERkQyioX4REZEMosIvIiKS\nQVT4RUREMogKv4iISAZR4RcREckgWenugIiMX2Z2/GtNABVAguTtWwF6nHNXpaVjInLO9HU+ETkj\nZvYI0OWc+8d090VEzp2G+kXknJhZl/e8xMzeMbOXzGyXmf2tmd1tZh9689PXeNuVmdnzZrbae1yd\n3gQimUmFX0RGw+eAbwK1wNeAWc65K0hOv/sdb5sfA//snLsc+ApDp+YVkTGic/wiMhpWO28aUzPb\nCbzutW8CrvWWrwfmJKdZAKDQzPKdc11j2lORDKfCLyKjoS9leSBlfYATP2dCwALnXO9YdkxEhtJQ\nv4iMldc5MeyPmV2axr6IZCwVfhEZK38EzDezjWa2heQ1ASIyxvR1PhERkQyiI34REZEMosIvIiKS\nQVT4RUREMogKv4iISAZR4RcREckgKvwiIiIZRIVfREQkg6jwi4iIZJD/B1RNWLcomLpMAAAAAElF\nTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig,ax=plt.subplots(1,1,figsize=(8,6))\n", "time = range(744)\n", "ax.plot(time, sa)\n", "ax.grid('on')\n", "ax.set_title('Number of Particles (init depth > 6) in Domain at time = x')\n", "ax.set_ylabel('Number of Particles')\n", "ax.set_xlabel('Time')" ] }, { "cell_type": "code", "execution_count": 22, "metadata": { "collapsed": false }, "outputs": [], "source": [ "index_shallow_particles=[]\n", "for n in range(2694):\n", " if Oct.variables['init_z'][n] < 6:\n", " index_shallow_particles.append(n)" ] }, { "cell_type": "code", "execution_count": 23, "metadata": { "collapsed": true }, "outputs": [], "source": [ "mask=lont[:].mask\n", "def shallow_particles(time, index):\n", " number_of_shallow_particles = np.zeros(time)\n", " for n in range(time):\n", " for m in index:\n", " if (mask[n,m]) == False: \n", " y,x = fast_ll2ij_SalishSea201702(lont[n,m],latt[n,m],lons, lats)\n", " if (598" ] }, "execution_count": 25, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf4AAAGDCAYAAADK03I6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XecXXWd//HXZ2bu9JZkkknvoYUSIISqJoBSFGFx5Ydr\nARZldbHt2sBdxYZld2XVtaysSlBKRERhXUFpAQQhhZ5GQnpIb1Mz9fP745wZbiZT7iR35tx75v18\nPO5jTj+f77137ud8v+d7zjF3R0RERIaGnKgDEBERkcGjxC8iIjKEKPGLiIgMIUr8IiIiQ4gSv4iI\nyBCixC8iIjKEKPELZjbfzL4R0b7NzG4zs71mtiiiGOrMbGofy0w2MzezvAHY/9NmdnI4/EUz+1mK\n66W8bA/rzzWzzYe7fh/bHrD3qx8xVJvZCjMr6GWZ/zazLw1mXAMlle/xYDGzt5jZqqjjkO4p8Wcg\nM1tvZjvMrCRp2ofNbGGEYQ2Uc4C3A+PdfU7XmWZ2tZm1hT9qNWb2opm963B3ZmYLzezDydPcvdTd\n1x7uNo+EmV0C1Lr7C2Es33T3D/exGl2XjTrRht/Z8yPa9/lm9ryZ1ZvZZjO7AsDdtwOPA9f1tK67\nf9Tdv36Y+11oZgfMrDb8bi41sxt6O9AYSAP1PU6lYhB+96YnxfKUux+d7lgkPZT4M1cu8Kmog+gv\nM8vt5yqTgPXuXt/LMn9191KgEvg5cI+ZDetnXGZmmfh9/yjwq6iDyGRmVt3LvOOAu4B/ASqAk4Cl\nSYvcCfzDAIb3cXcvA8YAnwGuBP5oZjaA+xQ5Mu6uV4a9gPXADcAeoDKc9mFgYTg8GXAgL2mdhcCH\nw+GrgaeB/wT2AWuBs8Lpm4AdwFVJ684H/ht4GKgFngAmJc0/Jpy3B1gFXNFl3Z8AfwTqgfO7Kc9Y\n4IFw/TXAR8Lp1wIHgDagDvhqN+teDfwlabwkLPtsYBjwB2AnsDccHt/lPbk5fC8aCZJAW7jPOuCH\n4XIOTA+Hi4DvAhuA/cBfwmkHvecESebnwFZgC/ANIDecNz18D/cDu4Bf9/A554dxJcf8FeCOLp/z\nVcDGcFv/0sOyG8Nl68LXmd3sryj8vPYCy4HPAZu7fE6/Dd/PdcAnu+zrXuDX4XfkeeCkcN6vgPaw\nLHXA5/uKPYX/gQTwN+H3praX5e4Cvt7L/DyggaTvc5f584FvhMNzgc0ECXxH+Nle08u2FxL+zyVN\nmxju713heAHwPeCN8PU9oKDL/j6ftL/LgIuB1wj+X76YtO05wF8J/qe3Aj8E8pPmJ3+P5wM/Av4v\n/LyeA6b1UpbfANsIvrNPAjPD6dcBLUBz+Nn+bzfrPhnuuz5c5v91lK3Lb9rngJfD5X4OVAMPhvE9\nAgxLWv4M4JmwrC8Bc9P1+xpu/6zwOzkhHD+J4P/imHTuJ1NfkQegVzcfSvBPcj5wX9KPUn8Tfytw\nDUHLwTcIfnx/FP4QvSP8ZysNl58fjr81nP99wmRLkGg3hdvKA04O/2GOS1p3P3A2QQtSYTfleRL4\nMVAIzCJILOcmxfqXXt6Lq5NiySNoBaklSLwjgPcAxUBZ+OP1+y7vyUZgZrhugu5/rJN/MH8ULjMu\nfO/OCt+Tg95z4HfAT8P3ZxSwCPiHcN7dBDXQnLDM5/RQtplAfZdpX+HQxP8/BEn7JKAJOLaXZfN6\neS+/DTwFDAcmAK8S/jiHsS4FvkxwQDKV4IDxgqR9tQB/G76PnyU4OEgkf2eT9tVr7L3EeAJwC0Ei\n/CtBbb2yl+XXAl8HXiFIhncAw7ss8zLw7h7Wn8/Bib8V+FpYxosJkviwHtY95LuU9H3/Tjj8NeDZ\n8DsykiCZfb3L/r4c7u8jBP8bdxF8n2cSHExNCZc/lSAh5oXv7wrg0z18j+cDuwkOFvIIDnoX9PI+\n/n24z44DlRe7e496Wb9z30ll65r4nyVI9uPCz/d5gt+TQuAx4KZw2XFh7BcTfC/fHo6P7GHffyA4\nQOju9YdeYr453G9R+P35eCq/z3F4RR6AXt18KG8m/uMJkupI+p/4VyfNOyFcvjpp2m5gVjg8P/lH\nASglqBlPIDh6f6pLfD9N+iedD/yyl7JMCLdVljTtW8D8pFj7Svyt4T/xrvDH45BWhXDZWcDeLu/J\n17os0/k+JU1zglp6DsEP7UndbLvzPQ9/vJqAoqT57wMeD4d/CdxKUk2+h3jPBrZ1mfYVDk3myS0C\ni4Are1m2t8S/Frgwafw63kz8pwMbuyx/I3Bb0r6eTZqXQ5Bo35L8ne3m/eo29m5iOxdYQnCQ+U3g\nqBT/V5rDfR9F8L39LXBnl2WeBj7Uw/rzOTjxN3Lw/9UO4Iwe1j3kuxROXwD8Tzj8OnBx0rwLCE5t\nJe+vo6WoLHzPTk9afilwWQ/7/zTwu67f46Ry/Sxp3sXAyhTf08pwWxVd36Ne1kkl8b8/afy3wE+S\nxj9BeNAOfAH4VZft/4mkVsp0vAgOtpYSJP2HAEvn9jP5FVmPW+mbu79qZn8gaPZf0c/VtycNN4bb\n6zqtNGl8U9J+68xsD0HT7yTgdDPbl7RsHgefl95Ez8YCe9y9NmnaBoKm+lQ96+7ndJ1oZsUEpzMu\nJGj2Bygzs1x3b0shtq6qCGofr/ex3CSCH42tSadyc5L29XmCWugiM9sLfNfdf9HNdvYS/Nj3ZVvS\ncAMHf279MZaD348NScOTgLFdPudcghaCDsnfkfbwioCxfewz1dhHERx8/ZWgaXdjH9vt0EhwcPIa\ngJl9k6DZOFkZwYFjKna7e2uKMfdkHEHNHoL3J/l93sDB79nupO9qY/i32/9TMzuKoDVkNkErVx4H\n92foKqX3PuyXczPwXoJKRns4q4qg4pEuXcvV0+/RJOC9YcfXDgmCjppp4+4tZjYf+AHwzx4eDQwF\nmdjZSQ52E0ET4LikaR0d4YqTpo0+wv1M6Bgws1KC5uA3CH7sn3D3yqRXqbt/LGnd3v5h3gCGm1ly\ngptIcF78SH0GOJqghlROcKoCILljVdfYeot1F8H5/2l97HcTQY2/Kuk9KXf3mQDuvs3dP+LuYwma\nqn+c3OM5yRqCfofjupnXX6n8aG0l6XMm+Bw6bALWdfmcy9z94qRlkr8jOcB4gs831f33yN0XEHyH\nf0XQ9+MNM/sfMzvkgK+Ll7vs+6A4wqscphMcTAw4M5tA0CTfccD0BkEi6zCRN9+z/voJsBKYEX7f\nv8jB3/XD9XfApQStjBUErTUkbXuwE+Imghp/8nexxN2/3d3CZvZgeNVPd68He9pJ+H93E3Ab8N2o\nrsaIghJ/hnP3NQQdqj6ZNG0nQeL8gJnlmtnf03ey6svFZnaOmeUT1FafdfdNBOfPjjKzD5pZInyd\nZmbHphj/JoLaz7fMrNDMTiT4Yb/jCOOFoCbXCOwzs+EE/8R92U5w/rq7WNuBXwC3mNnY8L09s+sP\ngrtvBf5M8GNRbmY5ZjbNzN4GYGbvNbPx4eJ7CX442+nC3ZsJaqdvS6WwfdgZ7qO367jvAW40s2Fh\nfJ9ImrcIqDWzL5hZUVj2483stKRlTjWzy8Nk+mmCg59nw3k9vq+pcvcD7n6Xu7+DoE/AeuA2M+ut\nBeY24Bozmxq2AN1A8J3tMIegaX1Dt2uniZkVh5///QTv5R/DWXcD/2pmI82siuB8/uF+98uAGqDO\nzI4BPtbH8v3ZbhPB6b9iglMtyVL5bI/4809yB3CJmV0Qfg8LLbjnxPjuFnb3i8LKSHevi7pbJ7zq\nYj5BJ8NrCQ6KD+uyzmykxJ8dvkbQiSzZRwh6ye4m6AT0TNeV+ukugsS5h6DG8gGAsIn+HQSXKb1B\n0Hz4HYJOQKl6H0Et4g2CTnE3uXvX5tjD8T2Cjjkd5/4fSmGd7wN/a8ENg37QzfzPEpzzW0zwXnyH\n7v9PPkTQCW45QXK/l+CSLoDTgOfMrI6gV/qnvOfrq38KfDCFuHvl7g2EVzCY2T4zO6Obxb5K0NS8\njuDA5VdJ67cB7yLoJ7GO4D39GUENsMP9BH0+9oYxX+7uLeG8bxEkuH1m9tk0lGeTu9/s7jMIrgzo\nablfEPSpeC4sWxNJB8nA+wmuWBkoPzSzWoLE9z2Cc9cXhgeREHSsXULQMvEKQYe2w71Z1mcJaue1\nBJ0mf30EcSf7JcF7t4Xg+/xsl/k/B44LP9vf97CNrwC3h8tccSTBhJWFSwlaNHYStAB8jvTmq08S\nnF76UtjEfw3BAeRb0riPjGVD6LSGSEYys6cJehS/EHUsPTGzrxB03vpA1LGkysxGEVxWebK7H4g6\nHpFMoc59IhFz97OjjiGO3H0HkNIpKZGhRE39IiIiQ4ia+kVERIYQ1fhFRESGECV+ERGRISS2nfuq\nqqp88uTJadtefX09JSVdr6iLh7iWLa7lgviWLa7lApUtG2VbuZYuXbrL3Uf2tVxsE//kyZNZsmRJ\n2ra3cOFC5s6dm7btZZK4li2u5YL4li2u5QKVLRtlW7nMLKUbVampX0REZAhR4hcRERlClPhFRESG\nECV+ERGRIUSJX0REZAhR4hcRERlClPhFRESGECV+ERGRIUSJX0REZAgZsMRvZr8wsx1m9mrStOFm\n9rCZrQ7/Dkuad6OZrTGzVWZ2QdL0U83slXDeD8zMBipmERGRuBvIGv984MIu024AHnX3GcCj4Thm\ndhxwJTAzXOfHZpYbrvMT4CPAjPDVdZsiIiKSogFL/O7+JLCny+RLgdvD4duBy5KmL3D3JndfB6wB\n5pjZGKDc3Z91dwd+mbTOoNnX0ExNs7NlX+Ng71pERCStBvscf7W7bw2HtwHV4fA4YFPScpvDaePC\n4a7TB9Un7n6BTz7WwKU/fHqwdy0iIpJWkT2dz93dzDyd2zSz64DrAKqrq1m4cOERb9PdmVnUxlPA\nrrom/vzo4+TnxqubQV1dXVreq0wT13JBfMsW13KBypaN4lquwU78281sjLtvDZvxd4TTtwATkpYb\nH07bEg53nd4td78VuBVg9uzZnq7HKc4DzB7hJy81Mf7YUzlubHlatpspsu3Rk6mKa7kgvmWLa7lA\nZctGcS3XYDf1PwBcFQ5fBdyfNP1KMyswsykEnfgWhacFaszsjLA3/4eS1hlUY0uDt2r1jtoodi8i\nIpIWA1bjN7O7gblAlZltBm4Cvg3cY2bXAhuAKwDcfZmZ3QMsB1qB6929LdzUPxJcIVAEPBi+Bt2w\ngqB5f3ddcxS7FxERSYsBS/zu/r4eZp3Xw/I3Azd3M30JcHwaQzssxYng7/7GlmgDEREROQK6c1+K\ncswoK8xT4hcRkaymxN8PFUUJJX4REclqSvz9oMQvIiLZTom/HyqLE+xrUOc+ERHJXkr8/aAav4iI\nZDsl/n6oKMpX4hcRkaymxN8PHTX+4HlBIiIi2UeJvx8qihK0tDmNLW19LywiIpKBlPj7oTK8i8++\nBjX3i4hIdlLi74eKoiDx6zy/iIhkKyX+fqhU4hcRkSynxN8P5UVq6hcRkeymxN8PHU39Narxi4hI\nllLi74fOzn2NunufiIhkJyX+fijOD55i3NCsy/lERCQ7KfH3Q26Okcg1DrS0Rx2KiIjIYVHi76fC\nvFyaWlXjFxGR7KTE308FiVzV+EVEJGsp8fdTQV4OTbplr4iIZCkl/n4qTOTQ1Koav4iIZCcl/n4q\nTORyQDV+ERHJUkr8/VSYyOWAOveJiEiWUuLvp+Acv5r6RUQkOynx95Nq/CIiks2U+PupMJGjy/lE\nRCRrKfH3U4Fu4CMiIllMib+fVOMXEZFspsTfTwV5upxPRESylxJ/PxXoBj4iIpLFlPj7qTAvl+bW\ndtrbPepQRERE+k2Jv58KE7kAqvWLiEhWUuLvp8JE8JapZ7+IiGQjJf5+KsgLavzq2S8iItlIib+f\nOmr86tkvIiLZSIm/n3SOX0REspkSfz8V5KnGLyIi2UuJv586avxK/CIiko2U+Pup4xz/+t31fPRX\nS9m4uyHiiERERFKnxN9PHb36n1y9i4eWbePq2xZFHJGIiEjqlPj7qaPGv23/AQBqm1qjDEdERKRf\nlPj7qaPGv2VvIwDuunWviIhkDyX+furo3LetJqjx72toUfIXEZGsocTfTwWJg9+y1nan5oCa+0VE\nJDso8fdTYdjUn2xvfXMEkYiIiPSfEn8/JXKNHAuGywryANitxC8iIllCib+fzKyzg9/MceUAbA/P\n94uIiGQ6Jf7DkMgNqvznTK8ix2DlttqIIxIREUmNEv9hKMoPavyTRpQwpaqEFVtrIo5IREQkNUr8\nh6EkPzi3P7qikBPHV7J4/R4am3XvfhERyXxK/IfhLTOqABhVVsD75kxkX0MLDy3bGnFUIiIifcuL\nOoBs9C/vPI6/OWU8k0aUMK6yiESu6Ty/iIhkhUhq/Gb2T2a2zMxeNbO7zazQzIab2cNmtjr8Oyxp\n+RvNbI2ZrTKzC6KIOVl+Xg6zJlQCkJebw6QRJazdWR9xVCIiIn0b9MRvZuOATwKz3f14IBe4ErgB\neNTdZwCPhuOY2XHh/JnAhcCPzezQu+hEaNrIEtburIs6DBERkT5FdY4/DygyszygGHgDuBS4PZx/\nO3BZOHwpsMDdm9x9HbAGmDPI8fZqwrBi3tina/lFRCTzWRQPmDGzTwE3A43An939/Wa2z90rw/kG\n7HX3SjP7IfCsu98Rzvs58KC739vNdq8DrgOorq4+dcGCBWmLua6ujtLS0m7nPfB6M/etbuF/3lFM\nouO2flmkt7Jls7iWC+JbtriWC1S2bJRt5Zo3b95Sd5/d13KD3rkvPHd/KTAF2Af8xsw+kLyMu7uZ\n9fuIxN1vBW4FmD17ts+dO/fIAw4tXLiQnra3qWA9961exkmnncmossK07XOw9Fa2bBbXckF8yxbX\ncoHKlo3iWq4omvrPB9a5+053bwHuA84CtpvZGIDw745w+S3AhKT1x4fTMkZFcT4ANY0tEUciIiLS\nuygS/0bgDDMrDpv0zwNWAA8AV4XLXAXcHw4/AFxpZgVmNgWYASwa5Jh7VVGUAGC/Er+IiGS4QW/q\nd/fnzOxe4HmgFXiBoHm+FLjHzK4FNgBXhMsvM7N7gOXh8te7e0bdJq8j8e9rUOIXEZHMFskNfNz9\nJuCmLpObCGr/3S1/M0FnwIxUqRq/iIhkCd2yNw3U1C8iItlCiT8NytXULyIiWUKJPw1yc4yywjzV\n+EVEJOMp8adJRVFCiV9ERDKeEn+aVBYr8YuISOZT4k8T1fhFRCQbKPGnSUVRgn0NzVGHISIi0isl\n/jSpKMpnf2Nr1GGIiIj0Sok/TSqKEtQ0thDF0w5FRERSpcSfJuVFeTS3tXOgpT3qUERERHqkxJ8m\npQXB3Y/rm9XcLyIimUuJP01K8sPE36TELyIimUuJP01Kwhp/nRK/iIhkMCX+NCkpyAWgoTmjnhgs\nIiJyECX+NFGNX0REsoESf5p0du5T4hcRkQymxJ8mxflhU3+TmvpFRCRzKfGnSama+kVEJAso8adJ\niZr6RUQkCyjxp0kiN4f8vBzq1atfREQymBJ/GpUW5FF7QI/mFRGRzKXEn0YVRQn2Nyrxi4hI5lLi\nT6OOxN/Q3MranXV6Up+IiGQcJf406ng070d+uYRzv/sEf319d9QhiYiIHESJP40qixPsa2xhxdZa\nAFbvqIs4IhERkYMp8adRR1N/jhkAW/Y1RhyRiIjIwfKiDiBOOhK/heOb9zZEGo+IiEhXSvxpVFGU\nwB06uvRt3qsav4iIZBY19adRVWlB53BujrGjpinCaERERA6lxJ9G00aWdg5PH1nKnoZmXdInIiIZ\nRYk/jaaNKukcPm5sOc2t7TToFr4iIpJBlPjTqDj/zS4TZ04dAcCe+uaowhERETmEOvel2SP//DZK\nC/J4dct+IEj8E4YXRxyViIhIQIk/zaaPCs7zv7E/6NG/p0E1fhERyRxq6h8gw4vzAdhVq579IiKS\nOfqV+M0sx8zKByqYOBk3rIiywjwWr98TdSgiIiKd+kz8ZnaXmZWbWQnwKrDczD438KFlt0RuDnOP\nHsXjq3ZGHYqIiEinVGr8x7l7DXAZ8CAwBfjggEYVEyeNr2BnbZN69ouISMZIJfEnzCxBkPgfcPcW\n3rwrrfRi6sjguv61O/WUPhERyQypJP6fAuuBEuBJM5sE1AxkUHExtSro4b92Z33EkYiIiAT6vJzP\n3X8A/CBp0gYzmzdwIcXH+GFFAGzW43lFRCRDpNK5r9rMfm5mD4bjxwFXDXhkMZCXm0NpQR61B1qi\nDkVERARIral/PvAnYGw4/hrw6YEKKG7KCvOoO9AadRgiIiJAaom/yt3vAdoB3L0V0JNnUlRWmEet\nEr+IiGSIVBJ/vZmNIOzJb2ZnAPsHNKoYKS3Io7ZJTf0iIpIZUrlX/z8DDwDTzOxpYCTwtwMaVYyU\nFSbYp/v1i4hIhkilV//zZvY24GjAgFXhtfySgtLCPDbtaYg6DBEREaCXxG9ml/cw6ygzw93vG6CY\nYqW8MI8aneMXEZEM0VuN/5Je5jmgxJ+CssIEdTrHLyIiGaLHxO/u1wxmIHFVWpDHgZZ2mlvbyc/T\nU5BFRCRaqdzA55tmVpk0PszMvjGwYcVHSUFwbNXYrCsgRUQkeqlUQS9y930dI+6+F7j4SHZqZpVm\ndq+ZrTSzFWZ2ppkNN7OHzWx1+HdY0vI3mtkaM1tlZhccyb4HW0l+LgD1zTrPLyIi0Usl8eeaWUHH\niJkVAQW9LJ+K7wMPufsxwEnACuAG4FF3nwE8Go533CL4SmAmcCHwYzPLPcL9D5risMbfoMQvIiIZ\nIJXEfyfwqJlda2bXAg8Dtx/uDs2sAngr8HMAd28OWxQuTdru7QSPASacvsDdm9x9HbAGmHO4+x9s\nHTX+BjX1i4hIBjB373shs4uA88LRh939T4e9Q7NZwK3AcoLa/lLgU8AWd68MlzFgr7tXmtkPgWfd\n/Y5w3s+BB9393m62fR1wHUB1dfWpCxYsONwwD1FXV0dpaWm/11uxu43vLD7AF04r5NgRmdlQcbhl\ny3RxLRfEt2xxLReobNko28o1b968pe4+u6/lUrlzH+7+IPDgEUf15j5PAT7h7s+Z2fcJm/WT9udm\n1vcRSRfufivBQQWzZ8/2uXPnpiHcwMKFCzmc7Q3btA8WP81Rxx3P3GOr0xZPOh1u2TJdXMsF8S1b\nXMsFKls2imu5emzqN7O/hH9rzawm6VVrZjVHsM/NwGZ3fy4cv5fgQGC7mY0J9zkG2BHO3wJMSFp/\nfDgtK5QUdHTuU1O/iIhEr8fE7+7nhH/L3L086VXm7uWHu0N33wZsMrOjw0nnETT7PwBcFU67Crg/\nHH4AuNLMCsxsCjADWHS4+x9sRfkdl/Opc5+IiESvz6Z+M/uVu3+wr2n99AngTjPLB9YC1xAchNwT\ndiDcAFwB4O7LzOwegoODVuB6d8+a6nPn5XxNWROyiIjEWCrn+Gcmj5hZHnDqkezU3V8EuuuAcF43\n03D3m4Gbj2SfUSnO1+V8IiKSOXo7x3+jmdUCJyaf3we282YzvPQhPy+HRK7pcj4REckIvZ3j/xZQ\nAfyyy/n9Ee5+4+CFmP2KErlK/CIikhF6vYGPu7cDpw1SLLFVUpBHfZOa+kVEJHqp3LnveTNT8j8C\nxfm5NLSoxi8iItFLpXPf6cD7zWwDUA8YwT12ThzQyGKkpCCPBtX4RUQkA6SS+LPqaXiZqCiRqxv4\niIhIRuizqd/dN7j7BqAR8KSXpKikIE+X84mISEboM/Gb2bvNbDWwDngCWE/67ts/JBTnq1e/iIhk\nhlQ6930dOAN4zd2nENxk59kBjSpmSvLzaNCd+0REJAOkkvhb3H03kGNmOe7+ON3fdU96UJSfS72a\n+kVEJAOk0rlvn5mVAk8S3F9/B0HvfklRSUEujc1tuDtmFnU4IiIyhKVS478UaAD+CXgIeB24ZCCD\nipvi/Dxa253mtvaoQxERkSGu1xq/mV0GTAdecfc/AbcPSlQx0/GEvoamNgryciOORkREhrLeHtLz\nY4Ja/gjg62b2pUGLKmY6ntCn8/wiIhK13mr8bwVOcvc2MysGniLo4S/9VBjW+A+0qKlfRESi1ds5\n/mZ3bwNw9waCW/XKYShKdCR+XdInIiLR6q3Gf4yZvRwOGzAtHNe9+vupMBEcXzUq8YuISMR6S/zH\nDloUMacav4iIZIoeE394f35Jg8Iw8Tfqtr0iIhKxVK7jlyPUkfgPtKpzn4iIREuJfxAUdfTqV41f\nREQi1tt1/I+Gf78zeOHEU2GeOveJiEhm6K1z3xgzOwt4t5ktoMvlfO7+/IBGFiOdNX4lfhERiVhv\nif/LwJeA8cAtXeY5cO5ABRU3heFtelXjFxGRqPXWq/9e4F4z+5K76459RyAnx8jPy9Gd+0REJHJ9\nPpbX3b9uZu8muIUvwEJ3/8PAhhU/RYlcNfWLiEjk+uzVb2bfAj4FLA9fnzKzbw50YHFTmMjRdfwi\nIhK5Pmv8wDuBWe7eDmBmtwMvAF8cyMDipiiRy4HWNhYs2siKrTV85d0zMdPjD0REZHClkvgBKoE9\n4XDFAMUSa6WFeextaOGG+14BYO4xo5h39KiIoxIRkaEmlRv4fAt4wczmh7X9pcDNAxtW/IwqK+TJ\n13Z2ji9dvzfCaEREZKhKpXPf3Wa2EDgtnPQFd982oFHFUHV5wUHja3fVRRSJiIgMZSk19bv7VuCB\nAY4l1kaVFXYOnzi+gmVv1PDipn3MmlAZYVQiIjLU6F79g6S6/M3Ef870KjbsbuCyHz3NnvrmCKMS\nEZGhRol/kIypfDPx//05UzqHX968L4pwRERkiOo18ZtZrpmtHKxg4uzsaVUATB9VSlVpAS/d9A4A\nXt68P8qwRERkiOn1HL+7t5nZKjOb6O4bByuoOMrPy+GlL78DxwGoKEowrrKItTvVyU9ERAZPKp37\nhgHLzGwRUN8x0d3fPWBRxVRFceKg8TEVhWyrORBRNCIiMhSlkvi/NOBRDFGjKwp5dYua+kVEZPD0\n2bnP3Z8A1gOJcHgx8PwAxzUkjC4PavzuHnUoIiIyRKTykJ6PAPcCPw0njQN+P5BBDRWjKwo50NLO\n/saWqEOrLf1+AAAgAElEQVQREZEhIpXL+a4HzgZqANx9NaCbzKfB6IrgEj+d5xcRkcGSSuJvcvfO\nu8yYWR6gtuk0GBMm/q37lfhFRGRwpJL4nzCzLwJFZvZ24DfA/w5sWENDx938tinxi4jIIEkl8d8A\n7AReAf4B+CPwrwMZ1FAxqqwQMyV+EREZPKk8na89fBzvcwRN/Ktc3dDTIj8vh6rSAiV+EREZNH0m\nfjN7J/DfwOuAAVPM7B/c/cGBDm4oGFdZxIY99X0vKCIikgap3MDnu8A8d18DYGbTgP8DlPjT4Ojq\nMh5esR13x8yiDkdERGIulXP8tR1JP7QWqB2geIaco0aXsae+mV11ejyviIgMvB4Tv5ldbmaXA0vM\n7I9mdrWZXUXQo3/xoEUYc8eMLgNg8fo9vP2WJ3jwla0RRyQiInHWW43/kvBVCGwH3gbMJejhXzTg\nkQ0RR4eJ/7t/XsXqHXV87M7nWbhqR8RRiYhIXPV4jt/drxnIHZtZLrAE2OLu7zKz4cCvgckEzwa4\nwt33hsveCFwLtAGfdPc/DWRsg6mqtICq0nxe31lPfl4OE4YV8YXfvswTn5tHYSI36vBERCRmUrlX\n/xQzu8XM7jOzBzpeadj3p4AVSeM3AI+6+wzg0XAcMzsOuBKYCVwI/Dg8aIiNGaOCWv9ls8byjctO\nYHtNE/cu3RxxVCIiEkepdO77PUEN/L8Ievh3vA6bmY0H3gn8LGnypcDt4fDtwGVJ0xe4e5O7rwPW\nAHOOZP+Z5hPnTgfgkpPGcsbU4UytKuHxlWruFxGR9Evlcr4D7v6DNO/3e8DngbKkadXu3tGzbRtQ\nHQ6PA55NWm5zOC02zppexYqvXUhRftCQMWtCJU+t2aVL/EREJO1SSfzfN7ObgD8DTR0T3f35w9mh\nmb0L2OHuS81sbnfLuLubWb/vDmhm1wHXAVRXV7Nw4cLDCbFbdXV1ad1eb4oaW9hZ28zv//Q4wwpT\naZQ5MoNZtsEU13JBfMsW13KBypaN4lquVBL/CcAHgXOB9nCah+OH42zg3WZ2McEVA+Vmdgew3czG\nuPtWMxsDdLR1bwEmJK0/Ppx2CHe/FbgVYPbs2T537tzDDPFQCxcuJJ3b603+ml3cufI5Rk47kXNm\nVA34/gazbIMpruWC+JYtruUClS0bxbVcqVQn3wtMdfe3ufu88HW4SR93v9Hdx7v7ZIJOe4+5+weA\nB4CrwsWuAu4Phx8ArjSzAjObAswAFh3u/rPB1JGlAKzdVRdxJCIiEjep1PhfBSp5swY+UL4N3GNm\n1wIbgCsA3H2Zmd0DLAdagevdvW2AY4lUdXkBJfm5rN2pe/iLiEh6pZL4K4GVZraYg8/xv/tId+7u\nC4GF4fBu4LwelrsZuPlI95ctzIzJVSWs26XELyIi6ZVK4r9pwKOQQ4ytLGLDbiV+ERFJrz4Tv7s/\nMRiByMHGVRbx7Ou7ow5DRERips/Eb2a1BL34AfKBBFDv7uUDGdhQN6aikNqmVmoOtFBemIg6HBER\niYlUavydN9mx4G4ylwJnDGRQAmMqg+cgbd13gPLRSvwiIpIe/bo7jAd+D1wwQPFIaNLwYgDW6ZI+\nERFJo1Sa+i9PGs0BZgMHBiwiAYLH9eYYLN9ay4XHj4k6HBERiYlUevVfkjTcSvDAnksHJBrpVJjI\nZUpVCSu21kQdioiIxEgq5/ivGYxA5FDTRpayYXcD9U2tXH3bItbtauAvX5hHYSJWTyUWEZFB1GPi\nN7Mv97Keu/vXByAeSTKqvIDF6/fw30+8zuL1ewHYsLuBo0eX9bGmiIhI93rr3FffzQvgWuALAxyX\nACNLC9nb0MJfk67n37SnIcKIREQk2/WY+N39ux0vgifeFQHXAAuAqYMU35A2qrwAgCUb9nLxCaMB\n2KjELyIiR6DXy/nMbLiZfQN4meC0wCnu/gV3H+gH9ggwqqygc/jcY6opyc/lkRXb2Vnb1MtaIiIi\nPesx8ZvZvwOLgVrgBHf/irvvHbTIhOryws7heUeP5ANnTOKZ13fzrv96ij31zRFGJiIi2aq3Gv9n\ngLHAvwJvmFlN+Ko1M11jNghmji3n4/Omc/28aYwoLeCGi47hv953Mttrmnhq9c6owxMRkSzUY69+\nd+/XXf0k/cyMz15w9EHjFx0/moK8HF7evJ9LZ42LMDoREclGSu5ZJi83h5ljy3lp076oQxERkSyk\nxJ+FTpsynJc276OhuTXqUEREJMso8Wehs6dV0dLmPPmazvOLiEj/KPFnoTOmjmBqVQnfe2R11KGI\niEiWUeLPQvl5Obzn1PGs3FbLrromnl27m5oDLVGHJSIiWUCJP0udOL4CgOt+uYQrb32Wm+5fFnFE\nIiKSDZT4s9SJ4yrJz83h+Y1B7/6HXt3GXt3UR0RE+qDEn6UqihP85AOn8N5Tx/P9K2fR2NLGZ37z\nUtRhiYhIhuvxBj6S+c47tprzjq0GYPnWGm59ci3b9h9gdEVhH2uKiMhQpRp/TFxy4ljcYcmGPVGH\nIiIiGUyJPyamjyolx2D19rqoQxERkQymxB8ThYlcJg4vZvWO2qhDERGRDKbEHyPTR5Wydmd91GGI\niEgGU+KPkfHDitmytxF3jzoUERHJUEr8MTKusojaplZqGvXwHhER6Z4Sf4yMH1YEwOZ9DRFHIiIi\nmUqJP0bGhYl/y97GiCMREZFMpcQfI1WlBQDs0a17RUSkB0r8MVJRlABgf6Oe1CciIt1T4o+R4vxc\nErnGPiV+ERHpgRJ/jJgZFUUJ1fhFRKRHSvwxU16UYH+DEr+IiHRPiT9mKlXjFxGRXijxx4ya+kVE\npDdK/DFTWZzP3gZdziciIt1T4o+ZUeUF7Khpor1d9+sXEZFDKfHHzMThxTS3tbO99kDUoYiISAZS\n4o+ZCcOKAdi4W/frFxGRQynxx8zE4WHi36PELyIih1Lij5lxw4rIz83h9Z31UYciIiIZSIk/ZhK5\nOUwfVcqKrTVRhyIiIhlIiT+Gjh1TzvKtNfzwsdX816OraW1rjzokERHJEHlRByDpd8qkSn77/Gb+\n48+vAbCnoZmbLpkZcVQiIpIJVOOPofOOqe4cPnliJfOfWc+za3dHGJGIiGQKJf4YGl1RyLcvP4F/\nfeex/PLv5zCqrIArb32WPy3bFnVoIiISMSX+mLpyzkQ+/JaplBUmuOsjZ1CQl8MPH1uDu+7oJyIy\nlA164jezCWb2uJktN7NlZvapcPpwM3vYzFaHf4clrXOjma0xs1VmdsFgx5ztpo0s5aZLZvLKlv38\nZsnmqMMREZEIRVHjbwU+4+7HAWcA15vZccANwKPuPgN4NBwnnHclMBO4EPixmeVGEHdWu/yUcRxd\nXcaNv3uFNTtquXfpZppb1dtfRGSoGfTE7+5b3f35cLgWWAGMAy4Fbg8Xux24LBy+FFjg7k3uvg5Y\nA8wZ3KizX2Eilx9/4BTa2p3zb3mSz/7mJb7+h+VRhyUiIoPMojzna2aTgSeB44GN7l4ZTjdgr7tX\nmtkPgWfd/Y5w3s+BB9393m62dx1wHUB1dfWpCxYsSFusdXV1lJaWpm17Ufnd6mZe2tlGbbOz+4Dz\nvmPyObuqKRZl6youn1l34lq2uJYLVLZslG3lmjdv3lJ3n93XcpFdx29mpcBvgU+7e02Q6wPu7mbW\n7yMSd78VuBVg9uzZPnfu3DRFCwsXLiSd24tKRxG27m/k0h8+zd0rmyicWcBN75obZVgDIi6fWXfi\nWra4lgtUtmwU13JF0qvfzBIESf9Od78vnLzdzMaE88cAO8LpW4AJSauPD6fJERhTUcRjn53L7EnD\nmL+smVc27486JBERGQRR9Oo34OfACne/JWnWA8BV4fBVwP1J0680swIzmwLMABYNVrxxVlqQxy+u\nOY3CPLjtmXVRhyMiIoMgihr/2cAHgXPN7MXwdTHwbeDtZrYaOD8cx92XAfcAy4GHgOvdvS2CuGOp\nvDDBaaPzeOjVberlLyIyBAz6OX53/wtgPcw+r4d1bgZuHrCghrgTqnJ5cnMTy97Yz8kTh/W9goiI\nZC3duU+YURl8DZZu2AtAS1s72/YfiDIkEREZIEr8QmVhDlNHlvDEazsB+Mw9L3HGtx7lOw+tjDgy\nERFJNyV+AeDtx1bz19d389CrW3ngpTcA+MnC1/n14o0RRyYiIumkxC9A8FCfHDM+esfzANz54dM5\nY+pwvnz/Mm5/Zj0HWoL+lPsamvWgHxGRLBbZDXwks0ypKuHu605n5bZa3jJ9JBNHFDNtZCkfv+t5\nbnpgGf/xp1WcOnkYT6/ZxXFjK/jPK05i6sjsuaOViIgEVOOXTqdOGs77T5/ExBHFAIyuKOSOD5/O\n96+cxfnHVbNk/V4qivLZsLue6+96gR216gAoIpJtVOOXXhUmcrl01jgunTWO/Q0ttLS3s3TDXv7x\nzuf54n2v8rOr+rwttIiIZBDV+CVlFcUJqkoLuGDmaK6fO41HVmxn056GqMMSEZF+UOKXw/Le2cHj\nE/74ytaIIxERkf5Q4pfDMmF4MSeOr1DiFxHJMkr8ctguPmEML23ez4ub9rFqWy3t7brMT0Qk0ynx\ny2F776njGVdZxGU/epoLvvckP3x8TdQhiYhIH5T45bCNKC3glitOIj83+Brd9vS6zhv9iIhIZlLi\nlyNy+tQRvPrVC7jrI6ezt6GF07/5KE++tpPm1nbd4U9EJAMp8csRy8/L4cypI7jhomMYXpLPh36x\niKP+9UE+esdSADbubuBXf13P5r269E9EJGq6gY+khZnx0bdN44rZE/jt0s28uGkf//fKVs6/5Qk2\n7K6npc0Z/fjr3PuxMxk/rDjqcEVEhiwlfkmr4SX5fOStU6k90EJxfi71za287aiRnDVtBP/06xe5\n/q4X+N3HziInx6IOVURkSFLilwFRVpjg39970kHTvnrpTP7p1y9x9+KNvP/0SRFFJiIytOkcvwya\ny2aNY86U4XzvkdU0tar3v4hIFJT4ZdCYGZ84dzo7a5t48JVtUYcjIjIkKfHLoDp7WhXjhxVxz5JN\nUYciIjIkKfHLoMrJMa46czLPvL6b+1/cEnU4IiJDjhK/DLqrz57MaZOH8akFLzL5hv/j1K8/zLb9\nB6IOS0RkSFCvfhl0idwcfnbVaXzy7hd44rWd7K5v5u9+9izvPXUCl5w0hvqmNh5evo3GljYWr9tL\nU2sbx4wu5+9On8jwknyWvVHD2MpCThxfGXVRRESyjhK/RKKiKMHtfz+HAy1tfPvBlcx/Zj3feWgl\n//6nlSQ/5G/GqFLKCvP49ZJN/DqpX0BejvHdK07ibUeNZFddE9NHlUVQChGR7KPEL5EqTOTypXcd\nxyUnjaGqtIC7nttIYSKXq86aTEVRgtzwRj8vbtrH9poDrN1ZT16O8dMn1/KpBS92buf8Y6vJMWgN\njxpyzJic18IpB1ooL0xEUjYRkUykxC+Ry80xTp00HIAbLz6222VmTTi4Wf9DZ03iu39+jR01B2hu\na+ePr2xjdHkhI8sKANjf2MIje5q5ffnDXH7yeEaU5rNk/V6GlST49PlHceyY8oEtlIhIhlLil6xU\nkJfLF5MOEvbUN1NZlOi8FXB7u/ODex9la84ofvfCFlra25k0vJiV22p4avUzPPX5eYwoLYgqfBGR\nyCjxSywML8k/aDwnx5g1Ko9Pzz2R7/ztiZ3T1+yo5fxbnuRjdzzPNy8/Xn0DRGTI0eV8MqRMH1XG\nd95zAq/tqOWKnz7LwlU7og5JRGRQqcYvQ87/O20isycP56O/WsrVty1mzuThjB9eRF6O8b45Ezl5\n4rCoQxQRGTCq8cuQNG1kKfd+7CyOHVPOS5v38djKHdyzZDN/8+NnuOXh16g90BJ1iCIiA0I1fhmy\nKooS/O4fz6K13SnMy2HJhr1898+r+MGjq3no1a3c8w9nUlmc3/eGRESyiBK/DGmFidzO4TOmjuA3\nHz2LJ1/byYdvX8Ksrz1MSX4w/3MXHM3VZ0+JKkwRkbRR4hfp4q1HjeRX185hweJN7G9sYev+A3zl\nf5ezblc9oyuKAMgxePWNGhqbW/n0+Udx/LiKiKMWEUmNEr9IN06fOoLTp44AYNkb+/nqA8v51bMb\nDrqdcGlBHvl5Obzrv/7CtJElvGXGSD7y1qmMqyyKKGoRkb4p8Yv0YebYCu756Jk0t7bT7o47rNpe\ny6ThxTS0tPHjx9ewZkcd859Zz/xn1nP5yeP41ntOYPX2OgoTOUwbWUpNYyvlRXmYWdTFEZEhTolf\nJEX5eW9eBNNxC+FhwM1/cwIAK7fVMP/p9SxYvIknV+9iV10TZjCusojNexuZPWkYN158LKdO6vly\nwT31zYfcjEhEJJ2U+EXS5JjR5Xz7PSfytqNG8pulm6kuL2BHTRMOvPPEMdz3/Bbe85NnGF1eSEfF\nP8eMkydWcvqU4Ty3bg9/eHkrHzhjIvOOHoUZbNjd0Pmwohc27uPqsyYzuaqEHQ3t7Kxt6nw2gYhI\nqpT4RdLsohPGcNEJYw6Z/qnzZvDLv25g7c66zmn7G1t44rWd/OHlreTnBi0Kdzy7kTue3djttu9/\ncQszx1bw9JpG/MlHAMjPzeHTb5/BxcePYdH6PRw3plydDUWkR0r8IoOkOD+Pj75t2iHT29qdFVtr\nmDiimPLCBOt21VPTGNxAqKQgjw2762ltdyaPKOFbD65gZ20TM6tyKS6r4JSJw3h+w17+7aFV/NtD\nqzq3+cWLj2HOlBG8smU/T722kzf2N1KQl8vZ00Zw1vQq2tqdJev3UpjIYfbk4eSFDzcyg6Oqy2hu\nayfHjNIC/USIxI3+q0UilptjB9XQp1SVHDR/+qjSzuH518wBYOHChcydeyYQPInw0ZU72FF7gGNG\nl/GDR9fwzT+u7FynoijBKRMreWPfAX7w2Bp+8NiaXuMpSuTS0taOGVQU5ZNjcOqkYZQXJg5Z1gze\nMmMkF58wGnfYVd/E/oYWpo4s7SybiGQWJX6RLJeTY7z9uOrO8duuPo1F6/fQ0NzKqLLCgw4qdtY2\n8cqWfRhB34J9DS2s3fXmqYfaA60sWb8XJ7husd1hX0MzL27cd9CljB0OtLaxYPEmplSVsGVfI82t\n7QfNv+j40RxVXcaLm/Z1TisrzGPe0aN496yxJHJ113CRwabELxIzOTnGGeE9CLoaWVbAuce8eZBQ\nWZzP5C4tDJfOGpfyvtranfue38w9SzZxysRhTKkqBoIDhu01B7jzuY08+Oo2jhldRkF4l8QXNu7l\nDy9v5bt/XsX4YcWUFebx8XOn6+FIIoNEiV9EDltujvHe2RN47+wJ3c6fObYCM3jfnImd05pb27l7\n0UZ+s3QTi9bvAeDRlTuoLE4wurCNZb6GC2aOZvqoUtxd9z4QSTMlfhEZMH93+sRDpuXn5XDVWZP5\n0JmTWL61hvGVxfx6yUbW7qznsWWb+fc/BQ9KOnt6FYvX7+H4sRXk5MDR1eVcOWcC+bk5vL6zjjuf\n28jkESWMG3bwnRJPGFfBqPAyx0Reju6kKNKFEr+IRMLMmDk26H9w3VuDqx0WjtjDyKNO5ku/f5XH\nVu6gtCCPlrZ22lqd2/+6nl88va5z/USu8fiqHXg3fQ+SnT19BGdNq+LkiZXkmJHINU4aX0me+hfI\nEKXELyIZZebYCn77sbNoaG6jJOlywvW76nlh014ADOPs6VWUFOTS0vZm5m9sbmPx+j20tgedDBet\n28MDL77B02t2H7SPqSNLOD486Jg5tpzRFYXk5eRw2pRhFIV9EUoLdItliSclfhHJOGZ2UNIHmFxV\nckhHxK4qihJcctLYzvG/OXk837r8RDbsrmfLvkYguLJh/jPreWXLfhqb23jgpTe63dakEcVUlxem\nHHNFUYILZ45mV11wt8Zl4dMbT5k0jFwzXl/XzCp7vXP5Y8eUM7aykNycHCYNLyZHlz7KIFHiF5HY\nmzSihEkj3jxo6Lhywd3ZvLeRlrZ2ttc0seyN/UDQcrBo/R5a2tq73V5X7Q6Pr9zBw8u3d04rLcij\nMJHDIyt2vLngqpXdrA2jygo4Z3oVFcWH3iuhJ8OK85k1oZI2d5au30t9c2vK6yZyczht8vDO1o2S\nglxmTahUC8cQkTWJ38wuBL4P5AI/c/dvRxySiGQ5M2PC8OASxKkjSzlzWveXQaZi4+4GVmyr4fQp\nw8nPyyE/N4ccMw60tgHw1JNP8Za3vgWAppZ2nlsXHFjsqmviL6t38ejKHbT31WGhg0Nt05uJ3ox+\n3WWxobmNW59ce9C0vBw77FaH9vZ2ch55kJL8XGZPHk4iN9jO8JJ8Thp/6AFFfVMrSzbsJZFrzJk8\nfMD6W5w4voKjqssOmlZzoGXI35UyK0puZrnAj4C3A5uBxWb2gLsvjzYyEZHAxBHFTBxRfMj04vzg\nZ7YgzzqHi/PhwuNHdy5zzdlT+r2/TXsa2FZzAAhOS4wqS/20RM2BFlZtq+0cX7uzjnW7GvodQ4eN\nGzcyceJE1u2qY+3OegCc4GCop+dOVJUW0NTSxn3Pbzns/fbFLLjKo+MOko3NbazcVotZ0LejrxtI\n1exv5PvLnx6w+N514liuPaf/n/2RyorED8wB1rj7WgAzWwBcCijxi8iQNGF4cWdrRX+VFyY4bfLw\nzvHk4cOxcOE25s495pDpNQda2N/Q0u06YyuLaGkLnjI5EFra2rl36WZe2bK/c1ppQR7nHjOKptZ2\nXtte28va4TbyBrZloCAvmitLsiXxjwM2JY1vBk6PKBYREUlBeWGi22c8dMjNyT3sg5dUfP7CQw9G\n+iN4Jkb8Uo15queUImRmfwtc6O4fDsc/CJzu7h/vstx1wHUA1dXVpy5YsCBtMdTV1VFaWtr3glko\nrmWLa7kgvmWLa7lAZctG2VauefPmLXX32X0tly01/i1A8j1Bx4fTDuLutwK3AsyePdvnzp2btgCC\nI7/0bS+TxLVscS0XxLdscS0XqGzZKK7lypZbVy0GZpjZFDPLB64EHog4JhERkayTFTV+d281s48D\nfyK4nO8X7r4s4rBERESyTlYkfgB3/yPwx6jjEBERyWbZ0tQvIiIiaaDELyIiMoQo8YuIiAwhSvwi\nIiJDiBK/iIjIEKLELyIiMoQo8YuIiAwhSvwiIiJDSFY8pOdwmNlOYEMaN1kF7Erj9jJJXMsW13JB\nfMsW13KBypaNsq1ck9x9ZF8LxTbxp5uZLUnlqUfZKK5li2u5IL5li2u5QGXLRnEtl5r6RUREhhAl\nfhERkSFEiT91t0YdwACKa9niWi6Ib9niWi5Q2bJRLMulc/wiIiJDiGr8IiIiQ4gSfx/M7EIzW2Vm\na8zshqjj6S8z+4WZ7TCzV5OmDTezh81sdfh3WNK8G8OyrjKzC6KJum9mNsHMHjez5Wa2zMw+FU6P\nQ9kKzWyRmb0Ulu2r4fSsLxuAmeWa2Qtm9odwPC7lWm9mr5jZi2a2JJwWl7JVmtm9ZrbSzFaY2Zlx\nKJuZHR1+Xh2vGjP7dBzK1it316uHF5ALvA5MBfKBl4Djoo6rn2V4K3AK8GrStH8DbgiHbwC+Ew4f\nF5axAJgSlj036jL0UK4xwCnhcBnwWhh/HMpmQGk4nACeA86IQ9nCeP8ZuAv4Q1y+j2G864GqLtPi\nUrbbgQ+Hw/lAZVzKllTGXGAbMCluZev6Uo2/d3OANe6+1t2bgQXApRHH1C/u/iSwp8vkSwn+kQn/\nXpY0fYG7N7n7OmANwXuQcdx9q7s/Hw7XAiuAccSjbO7udeFoInw5MSibmY0H3gn8LGly1perF1lf\nNjOrIKhA/BzA3ZvdfR8xKFsX5wGvu/sG4le2gyjx924csClpfHM4LdtVu/vWcHgbUB0OZ2V5zWwy\ncDJBzTgWZQubw18EdgAPu3tcyvY94PNAe9K0OJQLgoOzR8xsqZldF06LQ9mmADuB28JTND8zsxLi\nUbZkVwJ3h8NxK9tBlPiHOA/ar7L20g4zKwV+C3za3WuS52Vz2dy9zd1nAeOBOWZ2fJf5WVc2M3sX\nsMPdl/a0TDaWK8k54Wd2EXC9mb01eWYWly2P4HThT9z9ZKCeoPm7UxaXDQAzywfeDfym67xsL1t3\nlPh7twWYkDQ+PpyW7bab2RiA8O+OcHpWldfMEgRJ/053vy+cHIuydQibVB8HLiT7y3Y28G4zW09w\n2uxcM7uD7C8XAO6+Jfy7A/gdQRNwHMq2GdgctjoB3EtwIBCHsnW4CHje3beH43Eq2yGU+Hu3GJhh\nZlPCI8IrgQcijikdHgCuCoevAu5Pmn6lmRWY2RRgBrAogvj6ZGZGcM5xhbvfkjQrDmUbaWaV4XAR\n8HZgJVleNne/0d3Hu/tkgv+lx9z9A2R5uQDMrMTMyjqGgXcArxKDsrn7NmCTmR0dTjoPWE4Mypbk\nfbzZzA/xKtuhou5dmOkv4GKCHuOvA/8SdTyHEf/dwFagheDI/VpgBPAosBp4BBietPy/hGVdBVwU\ndfy9lOscgua3l4EXw9fFMSnbicALYdleBb4cTs/6siXFO5c3e/VnfbkIrvx5KXwt6/itiEPZwlhn\nAUvC7+TvgWExKlsJsBuoSJoWi7L19NKd+0RERIYQNfWLiIgMIUr8IiIiQ4gSv4iIyBCixC8iIjKE\nKPGLiIgMIXlRByAimcvMOi5rAhgNtBHcvhWgwd3PiiQwETlsupxPRFJiZl8B6tz9P6KORUQOn5r6\nReSwmFld+HeumT1hZveb2Voz+7aZvd/MFoXPp58WLjfSzH5rZovD19nRlkBkaFLiF5F0OAn4KHAs\n8EHgKHefQ/D43U+Ey3wf+E93Pw14Dwc/mldEBonO8YtIOiz28DGmZvY68Odw+ivAvHD4fOC44DEL\nAJSbWam71w1qpCJDnBK/iKRDU9Jwe9J4+/9v3w6NEIqBKIq+pSIEgqoo8JeBoBZcEIg/NBDEO6eC\nVbmTnSTnOXNJcl9rvXcOBvyy6gd2OXKu/TMz1z/OArWEH9jlkeQ2M8+ZeeX7JgDYzHc+ACjixg8A\nRR7Wv8YAAAAiSURBVIQfAIoIPwAUEX4AKCL8AFBE+AGgiPADQBHhB4AiH9wXGekTp8IhAAAAAElF\nTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig,ax=plt.subplots(1,1,figsize=(8,6))\n", "time = range(744)\n", "ax.plot(time, sb)\n", "ax.grid('on')\n", "ax.set_title('Number of Particles (init depth < 6) in Domain at time = x')\n", "ax.set_ylabel('Number of Particles')\n", "ax.set_xlabel('Time')" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.0" } }, "nbformat": 4, "nbformat_minor": 2 }