{ "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": "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": 5, "metadata": { "collapsed": true }, "outputs": [], "source": [ "result = nc.Dataset('/ocean/vdo/MEOPAR/ariane-runs/Jan17/ariane_trajectories_qualitative.nc')\n", "latt = result.variables['traj_lat']\n", "lont = result.variables['traj_lon']\n", "bathy, lons, lats = tidetools.get_bathy_data(grid)\n", "mask = lont[:].mask" ] }, { "cell_type": "code", "execution_count": 6, "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": 10, "metadata": { "collapsed": false }, "outputs": [], "source": [ "index_deep_particles=[]\n", "d = result.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": 13, "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": 15, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf4AAAGDCAYAAADK03I6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd4HNXZ9/Hvrd4ly7Lk3rsNGGzAgDE2mOYABgIE3tAT\nHNJJeR4CCQl5CCUFEgghQCiGUAyhhRCaQ2zAgDG2Affee5MtyUWypPP+MWOzFiprWavZ8vtc11ya\nnXqf2dXeO2fOnDHnHCIiIpIYkoIOQERERFqPEr+IiEgCUeIXERFJIEr8IiIiCUSJX0REJIEo8YuI\niCQQJX7BzCaY2W8C2reZ2eNmVmpm0wOKocLMejaxTHczc2aWEoH9f2BmR/vjN5vZI2GuF/ayDaw/\nyszWNnf9JrYdseN1CDGUmNkCM0tvZJkHzeyW1owrUsL5HLcWMzvZzBYFHYfUT4k/CpnZSjPbbGbZ\nIdO+aWZTAgwrUkYApwOdnXPH1Z1pZlebWY3/pVZmZp+Z2TnN3ZmZTTGzb4ZOc87lOOeWN3ebh8PM\nzgXKnXOf+rHc4Zz7ZhOrUXfZoBOt/5kdE9C+x5jZLDPbZWZrzewSAOfcJmAyML6hdZ1z1zvnbmvm\nfqeY2V4zK/c/mzPN7GeN/dCIpEh9jsM5MfA/e71DYnnfOdevpWORlqHEH72SgR8GHcShMrPkQ1yl\nG7DSOberkWU+cs7lAAXAo8DzZtbmEOMyM4vGz/v1wN+DDiJoZtbOzKwZ6w0EngF+DuQDRwEzQxZ5\nGvhWiwRZv+8553KBDsBPgEuB15tTFpFW45zTEGUDsBL4GbAdKPCnfROY4o93BxyQErLOFOCb/vjV\nwAfAH4EdwHLgRH/6GmAzcFXIuhOAB4FJQDnwLtAtZH5/f952YBFwSZ11/wq8DuwCxtRTno7Aq/76\nS4Hr/OnfAPYCNUAF8Ot61r0amBryOtsv+zCgDfAasAUo9cc71zkmt/vHYg9eEqjx91kB3O8v54De\n/ngmcDewCtgJTPWnHXTM8ZLMo8AGYB3wGyDZn9fbP4Y7ga3Acw28z2l+XKEx3wo8Ved9vgpY7W/r\n5w0su9pftsIfTqhnf5n++1UKzAf+B1hb53160T+eK4Af1NnXC8Bz/mdkFnCUP+/vQK1flgrgf5uK\nvZ7YbsT73P8a6HEI/yvPALc1Mj8F2E3I57nO/AnAb/zxUcBavAS+2X9vr2lk21Pw/+dCpnX193eO\n/zod+BOw3h/+BKTX2d//huzvfGAssBjv/+XmkG0fB3yE9z+9AbgfSAuZH/o5ngD8Bfi3/359DPRq\npCz/ADbifWbfAwb508cD+4Aq/739Vz3rvufve5e/zNf2l63Od9r/ALP95R4FSoA3/Pj+A7QJWX44\n8KFf1s+BUS31/epv/0T/M9nFf30U3v9F/5bcT7QOgQegoZ43xfsnGQO8FPKldKiJvxq4Bq/m4Dd4\nX75/8b+IzvD/2XL85Sf4r0f68+/FT7Z4iXaNv60U4Gj/H2ZgyLo7gZPwapAy6inPe8ADQAYwBC+x\nnBoS69RGjsXVIbGk4NWClOMl3rbAV4EsINf/8nqlzjFZDQzy102l/i/r0C/Mv/jLdPKP3Yn+MTno\nmAMvAw/5x6cYmA58y5/3LN4ZaJJf5hENlG0QsKvOtFv5cuL/G17SPgqoBAY0smxKI8fyLuB9oBDo\nAszF/3L2Y50J/BLvB0lPvB+MZ4bsax9wkX8cf4r34yA19DMbsq9GY28gvuF4PyK34VXRXwFkNfG/\nshy4DZiDlwyfAgrrLDMbOK+B9SdwcOKvBv7PL+NYvCTepoF1v/RZCvm8/9Yf/z9gmv8ZaYeXzG6r\ns79f+vu7Du9/4xm8z/MgvB9TPfzlh/rHKMU/vguAGxr4HE/wj+Nx/vJPAxMbOY7X+vvc/0Pls/qO\nUSPrH9h3SNnqJv5peMm+E94PnVl43ycZwH+BX/nLdvJjH4v3uTzdf92ugX2/hvcDob7htUZivt3f\nb6b/+fleON/P8TAEHoCGet6ULxL/YLyk2o5DT/xLQuYd4S9fEjJtGzDEH58Q+qUA5OCdGXfB+/X+\nfp34Hgr5J50APNlIWbr428oNmXYnMCEk1qYSf7X/T7zV//L4Uq2Cv+wQoLTOMfm/OsscOE4h0xze\nWXoS3hftUfVs+8Ax97+8KoHMkPmXAZP98SeBhwk5k28g3pOAjXWm3cqXk3lojcB04NJGlm0s8S8H\nzgp5PZ4vEv/xwOo6y98EPB6yr2kh85LwEu3JoZ/Zeo5XvbE3cVzSgUvwapG2A480smyVv+++eJ/b\nF4Gn6yzzAXBlA+tP4ODEv4eD/682A8MbWPdLnyV/+kTgb/74MmBsyLwz8S5the5vf01Rrn/Mjg9Z\nfiZwfgP7vwF4ue7nOKRcj4TMGwssbOrY+8sW+NvKr3uMGlknnMT/9ZDXLwJ/DXn9ffwf7Xi1P3+v\ns/23CKmlbIkB78fWTLyk/yZgLbn9aB6i8Zqn+Jxzc/F+zf6sGatvChnf42+v7rSckNdrQvZbgfeF\n2xHvGvzxZrZj/wB8HWhf37r16Ahsd86Vh0xbhferPlzTnHMFzrki59xw59x/AMwsy8weMrNVZlaG\nd6ZVUKedQWOx1VWEd/axrInluuF9aWwIOSYP4Z3VgVd1a8B0M5tnZtc2sJ1SvC/7pmwMGd/Nwe/b\noejIwcdjVch4N6Bjnff5ZrwfOfuFfkZq8aqpOzaxz0OO3TlXiXeW/hleYh/cyOJ78H6cLPY/t3fg\nJblQuXg/HMOxzTlXfagx19EJ7/8HvOMTepxXcfAx2+acq/HH9/h/6/0/NbO+ZvaamW30P+934H1m\nGxLWsTezZDO7y8yW+dtd6c9qbNvNUbdcDX0fdQMurvNZHIHXjqLFOOf24f2oGQzc7fxfA4lAiT/6\n/QqvCjA0Ue5vCJcVMi00ETdHl/0jZpaDVx28Hu/L/l0/8e4fcpxz3w5Zt7F/mPVAoZmFJriueNfF\nD9dPgH54Z0h5eJcqwEu6DcXWWKxb8a7/92piv2vwzviLQo5JnnNuEIBzbqNz7jrnXEe8hmUPhLZ4\nDrEUr93hofwIakg4X1obCHmf8d6H/dYAK+q8z7nOudAkGvoZSQI6472/4e6/UWbW1sy+59/W+V+8\nSy2jnXPDG1ltdp19HxSHf5dDb7zrxBFnZl3wquTf9yetx0tk+3Xli2N2qP4KLAT6+J/3mzn4s95c\n/w8Yh1fLmI9XW0PItls7Ia7BO+MP/SxmO+fuqm9hM3vDv+unvuGNhnbi/9/9CngcuDuouzGCoMQf\n5ZxzS/EaVP0gZNoWvMR5uf9r/VqaTlZNGWtmI8wsDe+a6TTn3Bq8Goe+ZnaFmaX6w7FmNiDM+Nfg\nXde808wyzOxIvEZ9Tx1mvOCdye0BdphZId4/cVM24V2/ri/WWuAx4B4z6+gf2xPqfiE45zYAb+N9\nWeSZWZKZ9TKzUwDM7GIz6+wvXor3xVlbz/6q8Bo1nRJOYZuwxd9HY/dxPw/cZGZt/Pi+HzJvOlBu\nZjeaWaZf9sFmdmzIMkPN7EI/md6A9+Nnmj+vweMaDjP7Bt6Z5il4Dfy6OOdudM4taGLVx4FrzKyn\nmWXh1Y69FjL/OLyq9VX1rt1C/NqnU4B/4h3L1/1ZzwK/8O9aKMK7nt/cz34uUAZUmFl/4NtNLH8o\n263Eu/yXhVeTECqc9/aw3v86ngLONbMz/c9hhnl9TnSub2Hn3Nn+yUh9w9n1rePfdTEBr5HhN/B+\nFDfrts5YpMQfG/4PrxFZqOvwWsluw2sE9OFh7uMZvMS5He+M5XIAv4r+DLzblNbjVR/+Fu86bLgu\nwzuLWI/XKO5X+6vrD9Of8Brm7L/2/2YY69wLXGReh0H31TP/p3jX/D7BOxa/pf7/kyvxGsHNx0vu\nL/BFVeSxwMdmVoF3N8MPXcP3Vz+E14jtsDjnduPfweBXj9Z3lvxrvKrmFXg/XP4esn4NcA5eO4kV\neMf0EbwzwP3+idfmo9SP+UK/uhS8dhu/8Pf902YU4SO8lvcXO+f+HVL93Sjn3GN4bSo+9stWSciP\nZLzLUg82I55w3W9m5XiJ7094167P8n9EgtewdgZezcQcvAZtze0s66d4Z+fleI0mnzuMuEM9iXfs\n1uF9nqfVmf8oMNB/b19pYBu3Ak/4y1xyOMH4Jwvj8Go0tuDVAPwPLZuvfoB3ae4Wv4r/GrwfkCe3\n4D6iliXQZQ2RqGRmH+C1KP406FgaYma34jXeujzoWMJlZsV4t1Ue7ZzbG3Q8ItEisO40RcTjnDsp\n6BjikXNuMxDWJSmRRKKqfhERkQSiqn4REZEEojN+ERGRBKLELyIikkAi1rjPzB7Duz1os3NusD/t\nObwOV8DrFnKHc26ImXXH63d6//ObpznnrvfXGYp3v2Um3r2xPwynh6WioiLXvXv3lioOu3btIju7\n7h118SFeyxav5YL4LVu8lgtUtlgUa+WaOXPmVudcuyYXjFRfwHi9qB0DzG1g/t3AL/3x7o0sNx3v\nwRSG9ySns8PZ/9ChQ11Lmjx5cotuL5rEa9nitVzOxW/Z4rVczqlssSjWygXMcEH21e+ce48v+qs+\niN9r0iV4vVo1yMw6AHnOuWl+oZ7Ee2yliIiINENQ1/hPBjY555aETOthZp+Z2bshvSd1wnsQyH5r\nObSHu4iIiEiIiN7O51+7f8351/hDpv8VWOqcu9t/nY73bPht/jX9V/C6oe0L3OWcG+MvdzJwo3Pu\nnAb2Nx7vUaOUlJQMnThxYouVpaKigpyc5j4ULbrFa9nitVwQv2WL13KByhaLYq1co0ePnumcG9bU\ncq3ec5//gI8L8fqDBw48hrPSH59pZsvwkv46vCeA7deZRp7q5px7GO856AwbNsyNGjWqxeKeMmUK\nLbm9aBKvZYvXckH8li1eywUqWyyK13IFUdU/BljonDtQhe8/uSrZH+8J9AGWO+8paGVmNtxvF3Al\n3oNCREREpBkilvjN7Fm8J271M7O1/mM3wXvKW91GfSOB2Wb2Gd5Tzq53zu1vGPgdvKeELQWW4bXs\nFxERkWaIWFW/c+6yBqZfXc+0F/EeZ1nf8jOAwfXNExERkUOjnvtEREQSiBK/iIhIAlHiFxERSSBK\n/CIiIglEiT8Ms1aXMntLNR8u20pNbeQ6PBIREYm0Vu/AJxbd9fpCpq+shJkf8/Xju3L7BUcEHZKI\niEiz6Iw/DLdfMJhbhmdw4TGdeO6TNeyuqg46JBERkWZR4g9Dn5JcehUkc+6RHamudXy2ZkfQIYmI\niDSLEv8hOKZbGwBmrSoNOBIREZHmUeI/BPmZqXQtzGLBxvKgQxEREWkWJf5D1L99Lgs3lB14Xbqr\nivPun8rMVdsbWUtERCQ6KPEfov7tc1mxdReV1TUAPPnRKmav3cnD7y0PODIREZGmKfEfos5tsqh1\nsGlnJQBTl24BYGtFVZBhiYiIhEWJ/xB1KMgAYMPOPdTUOuat96r9Z64q5XvPzKK6pjbI8ERERBql\nxH+IOuRnArBh515WbK1gd1UNv/jKAAqz03ht9gb+9v6KgCMUERFpmBL/IeqQ753xr9+5hznrdgIw\nok8Rs245nbMHt+eP/1nM0s0VQYYoIiLSICX+Q5SdnkKbrFRWb9vNnLVlZKQm0btdDgC/HjeIjJQk\n/vSfxQFHKSIiUj/11d8MAzrkMX9DGRkpyQzokEdKsvf7qTg3g9MGlPD+kq045zCzgCMVERE5mM74\nm2FghzwWbCjj87U7OKJT/kHzhnVvw9aKSpZt2RVQdCIiIg1T4m+GQZ3y2FfjqKyuZXCdxH9a/xKS\nk4z/e20+O3brFj8REYkuSvzNMKjjF8l+SJeCg+a1z8/gupN78v6SLXztoWm6vU9ERKKKEn8z9CzK\nPjDepzjnS/N/dnZ//vS1ISzaVM7PX56r5C8iIlFDjfuaISU5ifEje9K1MKvBBnznHdWReevLePi9\n5XQryuI7o3q3cpQiIiJfpjP+Zrp57AAuH96twflmxs1jB3BCz7Y8O301zrlWjE5ERKR+SvwRdtbg\n9qzZvofN5ZVBhyIiIqLEH2k9/PYAq7btDjgSERERJf6I69Y2C4CV23Rfv4iIBE+JP8I6FWSSkmSs\n1hm/iIhEASX+CEtJTqI4N52NZXuDDkVERESJvzUU5qSxrUKN+0REJHhK/K2gMDud7bvUfa+IiARP\nib8VtM1OY5sSv4iIRAEl/lZQmJ2mM34REYkKSvytoDA7jd1VNeypqgk6FBERSXBK/K2gXU46AJvL\n1bJfRESCpcTfCgZ0yANgzrqdAUciIiKJLmKJ38weM7PNZjY3ZNqtZrbOzD7zh7Eh824ys6VmtsjM\nzgyZPtTM5vjz7rOGHocXxfp3yCU9JYlPV+8IOhQREUlwkTzjnwCcVc/0PzrnhvjD6wBmNhC4FBjk\nr/OAmSX7y/8VuA7o4w/1bTOqpSYncXTXAj5cti3oUEREJMFFLPE7594Dtoe5+DhgonOu0jm3AlgK\nHGdmHYA859w05z3X9kng/MhEHFmj+hWzYEMZSzeXc9tr81m5VX33i4hI60sJYJ/fN7MrgRnAT5xz\npUAnYFrIMmv9afv88brT62Vm44HxACUlJUyZMqXFgq6oqDis7eWW1wIw5p73ANiyYS0X9klridAO\n2+GWLVrFa7kgfssWr+UClS0WxWu5Wjvx/xW4DXD+37uBa1tq4865h4GHAYYNG+ZGjRrVUptmypQp\nHM72nHP8ec5/D/TZX5tdxKhRx7RQdIfncMsWreK1XBC/ZYvXcoHKFovitVytmvidc5v2j5vZ34DX\n/JfrgC4hi3b2p63zx+tOjzlmxq3nDeLdxVuYv34nSzZVBB2SiIgkoFa9nc+/Zr/fBcD+Fv+vApea\nWbqZ9cBrxDfdObcBKDOz4X5r/iuBf7ZmzC3prMHtufPCIziuRyErtu3Ca7YgIiLSeiJ2xm9mzwKj\ngCIzWwv8ChhlZkPwqvpXAt8CcM7NM7PngflANfBd59z+bu6+g3eHQCbwhj/EtLY56VRV17Krqoac\n9CCaWYiISKKKWNZxzl1Wz+RHG1n+duD2eqbPAAa3YGiBK8z2GvWV7qpS4hcRkValnvsC0NZP/Hpi\nn4iItDYl/gC08RP/9l2VAUciIiKJRok/AAfO+Ct0xi8iIq1LiT8AB67x71biFxGR1qXEH4Cc9BTS\nkpN0jV9ERFqdEn8AzIw22alsV1W/iIi0MiX+gBRmp6uqX0REWp0Sf0DaZqepql9ERFqdEn9ACrPT\n2K7ELyIirUyJPyCF2Wm6xi8iIq1OiT8ghdlplFdWU1VdG3QoIiKSQJT4A6J7+UVEJAhK/AFR730i\nIhIEJf6AFB7or1+JX0REWo8Sf0AOJH5V9YuISCtS4g/IgcRfoSf0iYhI61HiD0hBVhpmquoXEZHW\npcQfkOQkoygnnTWle4IORUREEogSf4BO6tWWKYs2U12je/lFRKR1KPEH6IxB7SndvY8Zq0qDDkVE\nRBKEEn+ARvZtR1pKEm/P2xR0KCIikiCU+AOUk57CiN5FvD1/I865oMMREZEEoMQfsDMHlbC2dA/z\nN5QFHYqIiCQAJf6AnTagBDNU3S8iIq1CiT9gRTnpHNu9kJc+Xasn9YmISMQp8UeBb4/qxZrte5j4\nyeqgQxERkTiXEnQAAqP6tmN4z0J+89oC3l+ylQHtc7lhTF+Skizo0EREJM7ojD8KmBk/O3sA+2pr\nmTR/E/f9dylPfLQy6LBERCQOKfFHiSFdCpj0o5G885NTGN2vHb99cyGrtu0KOiwREYkzSvxRpHdx\nLr3a5XDnhUeSkpTEXW8sDDokERGJM0r8Uah9fgZXntCNt+ZtZM323UGHIyIicUSJP0pdcUI3ksyY\n8OHKoEMREZE4osQfpTrkZ3L6wBJe/Xw9tbXqzldERFqGEn8UO31gCVvKK5m7fmfQoYiISJxQ4o9i\np/Rthxm8s2Bz0KGIiEicUOKPYm1z0jm6SwGTFynxi4hIy4hY4jezx8xss5nNDZn2ezNbaGazzexl\nMyvwp3c3sz1m9pk/PBiyzlAzm2NmS83sPjNLqO7sThtQwuy1O7nvnSWs2Kr7+kVE5PBE8ox/AnBW\nnWmTgMHOuSOBxcBNIfOWOeeG+MP1IdP/ClwH9PGHutuMa6f2LwbgnkmLueShj6hRQz8RETkMEUv8\nzrn3gO11pr3tnKv2X04DOje2DTPrAOQ556Y55xzwJHB+JOKNVgM65HHXhUfQoyibLeWVfLhsa9Ah\niYhIDDMvn0Zo42bdgdecc4Prmfcv4Dnn3FP+cvOAJcBO4BfOuffNbBhwl3NujL/OycCNzrlzGtjf\neGA8QElJydCJEye2WFkqKirIyclpse0dqn21jh/8dzfDO6Rw1aD0Ft120GWLlHgtF8Rv2eK1XKCy\nxaJYK9fo0aNnOueGNbVcIE/nM7OfA9XA0/6kDUBX59w2MxsKvGJmgw51u865h4GHAYYNG+ZGjRrV\nQhHDlClTaMntNUfPue9DVjqjRh3XotuNhrJFQryWC+K3bPFaLlDZYlG8lqvVE7+ZXQ2cA5zmV9/j\nnKsEKv3xmWa2DOgLrOPgywGd/WkJqSQvg4079wYdhoiIxLBWvZ3PzM4C/hc4zzm3O2R6OzNL9sd7\n4jXiW+6c2wCUmdlwvzX/lcA/WzPmaFKSl87m8sqgwxARkRgWydv5ngU+AvqZ2Voz+wZwP5ALTKpz\n295IYLaZfQa8AFzvnNvfMPA7wCPAUmAZ8EakYo52xbkZbNtVyb6a2qBDERGRGBWxqn7n3GX1TH60\ngWVfBF5sYN4M4EuNAxNRSV4GzsHWiko65GcGHY6IiMQg9dwXQ4pzvdb8m8pU3S8iIs2jxB9DSvIy\nANhcpgZ+IiLSPEr8MaQkzz/jVwM/ERFpJiX+GNI2J50k0xm/iIg0nxJ/DElOMopy0tmkxC8iIs2k\nxB9jSvIydC+/iIg0mxJ/jCnJS1erfhERaTYl/hhTnJeha/wiItJsSvwxpiQ3g227qqiqVu99IiJy\n6JT4Y0yxf0vf6u27m1hSRETky5T4Y0x7vxOfMfe8qyp/ERE5ZEr8Meak3kVcd3IPAJ7+eHXA0YiI\nSKxR4o8xaSlJ/PwrAzm+RyHPfbKGzeU66xcRkfAp8ceoG8/uz/ZdVdz7nyVBhyIiIjFEiT9GHdO1\nDeOGdOSlWevYU1UTdDgiIhIjlPhj2HlDOrJnXw0n/24yFzzwAff/dwnOuaDDEhGRKHZIid/Mksws\nL1LByKE5vkdbSvLSKduzj12V1fzh7cU89sHKoMMSEZEoltLUAmb2DHA9UAN8AuSZ2b3Oud9HOjhp\nXFpKEu/972icg/SUJL7xxAzufnsR44Z0pCgnPejwREQkCoVzxj/QOVcGnA+8AfQArohoVBK29JRk\nMlKTMTNuHjuA6hrH956ZFXRYIiISpcJJ/KlmloqX+F91zu0DdCE5CvUuzuGHY/owbfl21qhnPxER\nqUc4if8hYCWQDbxnZt2AskgGJc13zpEdAPjX7PUBRyIiItGoycTvnLvPOdfJOTfWeVYBo1shNmmG\nbm2zObFXW578cJUe5CMiIl/SZOI3sxIze9TM3vBfDwSuinhk0mzjR/ZkY9lefvqPzznl95OZu25n\n0CGJiEiUCKeqfwLwFtDRf70YuCFSAcnhO6VvO/q3z+XVz9ezattufvfWoqBDEhGRKBFO4i9yzj0P\n1AI456rxbu2TKGVm3HhWf47olM8pfdvx2epSdewjIiJAeIl/l5m1xW/Jb2bDAdUdR7nR/Yv51/dH\ncNbg9pTtrWa1WvmLiAhhdOAD/Bh4FehlZh8A7YCLIhqVtJiBHbyOFhdtLKdb2+yAoxERkaA1mfid\nc7PM7BSgH2DAIv9efokBXQqzAFhbuifgSEREJBo0mPjN7MIGZvU1M5xzL0UoJmlBbbJSyUpLZk2p\nqvpFRKTxM/5zG5nnACX+GGBmdGmTpTN+EREBGkn8zrlrWjMQiZxObTJZp8QvIiKE14HPHWZWEPK6\njZn9JrJhSUtql5POlorKoMMQEZEoEM7tfGc753bsf+GcKwXGRi4kaWlFuWls31VFba3u5RcRSXTh\nJP5kMzvwcHczywT0sPcYUpSTTk2to3R3VdChiIhIwMK5j/9p4B0ze9x/fQ3wRORCkpZWlOP9Ttta\nUUXbHP1mExFJZOE8ne+3wO3AAH+4zTn3u6bWM7PHzGyzmc0NmVZoZpPMbIn/t03IvJvMbKmZLTKz\nM0OmDzWzOf68+8zMDrWQia5d7v7Er+v8IiKJLpyqfpxzbzjnfuoPb4W57QnAWXWm/Qx4xznXB3jH\nf73/iX+XAoP8dR4ws2R/nb8C1wF9/KHuNqUJ7fMyAFi+dVfAkYiISNAaTPxmNtX/W25mZSFDuZmV\nNbVh59x7wPY6k8fxxWWCJ4DzQ6ZPdM5VOudWAEuB48ysA5DnnJvmvKfMPBmyjoSpW9ssehfn8MzH\nq9m7T89XEhFJZA0mfufcCP9vrnMuL2TIdc7lNXN/Jc65Df74RqDEH+8ErAlZbq0/rZM/Xne6HAIz\n44YxfVi4sYybX57DvpraoEMSEZGANNm4z8z+7py7oqlph8o558ysRe8vM7PxwHiAkpISpkyZ0mLb\nrqioaNHttbYcYEzXFF6atY7V6zby3SHp7G8uEetla0i8lgvit2zxWi5Q2WJRvJYrnFb9g0JfmFkK\nMLSZ+9tkZh2ccxv8avzN/vR1QJeQ5Tr709b543Wn18s59zDwMMCwYcPcqFGjmhnml02ZMoWW3F4Q\nRo50/O6tRTz47jLa9T2GwZ3ycc7x7rvvxnzZ6hMP71lD4rVs8VouUNliUbyWq7Fr/DeZWTlwZOj1\nfWAT8M9m7u9V4Cp//KqQ7bwKXGpm6WbWA68R33T/skCZmQ33W/NfeRj7TnhJSca3RvYkJcm4840F\nnP+XDzjnz1OpqFLHPiIiiaKxvvrvNLPfAo8456491A2b2bPAKKDIzNYCvwLuAp43s28Aq4BL/H3N\nM7PngflANfBd59z+VmjfwbtDIBN4wx+kmdpkp3HNSd352/srDkx7sjaZM06tJS0lrJs8REQkhjVa\n1e+cqzU8msCZAAAgAElEQVSzY5uzYefcZQ3MOq2B5W/H6y+g7vQZwODmxCD1u3nsAC49risleRn8\n4a1FTPhwJeP/PoPHrz72wHV/ERGJT+Gc4s1qbvKX6GRm9GqXQ056Cr88ZyAX9E5lyqItfLZmR9Mr\ni4hITAsn8R8PfGRmy8xstt+L3uxIByatIynJOKN7KhmpSbw4a23TK4iISEwLp1X/mU0vIrEsM8U4\nc1B7/vX5Bm45ZyDpKclNryQiIjEpnL76VznnVgF7ABcySBy58JjO7Nyzj3cWbGbZlgrG3PMuM1bW\n7XhRRERiXTgd+JwH3A10xLvvvhuwgDr390tsG9G7iJK8dJ7+eBVby6tYurmCSQs2Max7YdChiYhI\nCwrnGv9twHBgsXOuB16r/GkRjUpaXXKScf7Rnfhg6TYWbSqnKCeNT1ersZ+ISLwJ5xr/PufcNjNL\nMrMk59xkM/tTxCOTVnfVCd1ZvW03o/sXM399Gc/PWENtrSMpSbf4iYjEi3AS/w4zywHeA542s82A\nnu8ahzoWZPLXy73emCdOX83uqhrWlu6ha9usgCMTEZGWEk5V/zhgN/Aj4E1gGXBuJIOS4PVrnwvA\nwo1NPoFZRERiSKNn/GZ2PtAbmOOcewt4olWiksD1KfES/+JN5ZwxqH3A0YiISEtp7CE9D+Cd5bcF\nbjOzW1otKglcTnoKXQozWbixPOhQRESkBTV2xj8SOMo5V2NmWcD7eC38JUH0K8ljkRK/iEhcaewa\nf9X+J+Q553YDatqdYAZ0yGX51l1sq6gMOhQREWkhjSX+/n7f/LPNbE7Ia/XVnyDOO6ojNbWOJz5c\nGXQoIiLSQhqr6h/QalFIVOpTkstXjuzAQ+8t54oTutMuNz3okERE5DA1eMa/v4/+hobWDFKC86Mx\nfaisruUlPblPRCQuhHMfvySw3sW5HNejkMc/WMnefTVBhyMiIodJiV+a9I0RPdhYtpfP1qjvfhGR\nWNfYffzv+H9/23rhSDQa2CEPgOVb1FOziEisa6xxXwczOxE4z8wmUud2PufcrIhGJlGjU0Em6SlJ\nLN9SEXQoIiJymBpL/L8EbgE6A/fUmeeAUyMVlESXpCSjR1E2K7Z6Z/ylu6qYuaqUU/sX68l9IiIx\npsHE75x7AXjBzG5xzqnHvgTXs102CzZ4vfj99B+f887CzYw9oj33X3aMkr+ISAxpsnGfc+42MzvP\nzP7gD+e0RmASXXoW5bB6+27eX7KFdxZuJjstmdfnbKT/L99kq3r2ExGJGU0mfjO7E/ghMN8ffmhm\nd0Q6MIkuvYqzqal13PTSHIpz0/nkF2PoXZxDVXWtevYTEYkh4dzO9xXgdOfcY865x4CzAJ31J5jB\nHfMBWFu6h4uHdSYrLYVJPxrJyX2KeH3OhoCjExGRcIV7H39ByHh+JAKR6NanJJd7Lx3Ct0b25Nuj\negNgZpzQqy3Ltuxi+66qgCMUEZFwNNaqf787gU/NbDLeLX0jgZ9FNCqJSuOGdGLckE4HTTu+R1sA\nHnl/OT85ox/JaugnIhLVwmnc9ywwHHgJeBE4wTn3XKQDk9hwTNcCLjy6Ew9MWcaAW97k8Q9WBB2S\niIg0IqyqfufcBufcq/6wMdJBSewwM+6+5CgevPwYehfncMfrC9hTpT79RUSilfrql8NmZpw1uAPX\nj+rFvhrH6u27gw5JREQaoMQvLaZbYRYAq7apT38RkWjVaOI3s2QzW9hawUhs69Z2f+LXGb+ISLRq\nNPE752qARWbWtZXikRhWkJVGm6xUlm7Ww3xERKJVOLfztQHmmdl04EAdrnPuvIhFJTHryM4FfL52\nR9BhiIhIA8JJ/LdEPAqJG0d3LeDed5ZQUVlNTno4Hy8REWlN4dzH/y6wEkj1xz8BZkU4LolRR3dt\ng3Mwe43O+kVEolE4D+m5DngBeMif1Al4pbk7NLN+ZvZZyFBmZjeY2a1mti5k+tiQdW4ys6VmtsjM\nzmzuviXyhnT2end+8qNVzFu/M+BoRESkrnBu5/sucBJQBuCcWwIUN3eHzrlFzrkhzrkhwFBgN/Cy\nP/uP++c5514HMLOBwKXAILwHBD1gZsnN3b9EVn5WKkO6FPDmvI185b6p3P/fJUGHJCIiIcJJ/JXO\nuQNPYDGzFMC10P5PA5Y551Y1ssw4YKJzrtI5twJYChzXQvuXCJhwzbH87cphjOzbjgemLKOqujbo\nkERExGfONZ7Dzex3wA7gSuD7wHeA+c65nx/2zs0eA2Y55+43s1uBa4CdwAzgJ865UjO7H5jmnHvK\nX+dR4A3n3Av1bG88MB6gpKRk6MSJEw83xAMqKirIyclpse1Fk0iVbdamau77tJIbj81gQNvWr6TR\nexZ74rVcoLLFolgr1+jRo2c654Y1uaBzrtEBr1bgOuAfeNf6r8P/wXA4A5AGbAVK/NclQLK/v9uB\nx/zp9wOXh6z3KHBRU9sfOnSoa0mTJ09u0e1Fk0iVbcfuKtftxtfcXyYvicj2m6L3LPbEa7mcU9li\nUayVC5jhwsi/Td5v5ZyrNbMngI/xqvgX+Ts4XGfjne1v8vezaf8MM/sb8Jr/ch3QJWS9zv40iXL5\nmal0yM9g8cbyoEMRERFfOK36vwIsA+7DO/teamZnt8C+LwOeDdlPh5B5FwBz/fFXgUvNLN3MegB9\ngOktsH9pBX1Lclm8ST35iYhEi3B6WLkbGO2cWwpgZr2AfwNvNHenZpYNnA58K2Ty78xsCF6twsr9\n85xz88zseWA+UA1813ldCUsM6FOcw7Tl23DOYWZBhyMikvDCSfzl+5O+bzlwWHW3zrldQNs6065o\nZPnb8a77S4zp1jaLyupatpRXUpyXEXQ4IiIJr8HEb2YX+qMzzOx14Hm8s/GL8XrvE2lSF/9Rvau3\n76Y4L0Nn/iIiAWvsjP/ckPFNwCn++BYgM2IRSVzpWvjFo3qnLt3KI++v4Gdn9+fy4d0CjkxEJDE1\nmPidc9e0ZiASnzq1ySQtOYmf/OPzA9N+8cpcehRlc1LvogAjExFJTOG06u9hZveY2Utm9ur+oTWC\nk9iXnpLMOUd5N2yMPaI9//7BCACuefwT9eUvIhKAcBr3vYLXac6/APW9Kofsl+cM5PQBJZwxqD3J\nScbUG0dz7p+n8ts3F/Hktep9WUSkNYWT+Pc65+6LeCQStwqy0jj7iC+6aejcJotrT+rB3ZMWs7Z0\nN53bZAUYnYhIYgnnIT33mtmvzOwEMztm/xDxyCSujRvSCYDnZ6wNOBIRkcQSzhn/EcAVwKl8UdXv\n/NcizdK1bRZjj2jPQ+8uA+CUvu0Y2q1NwFGJiMS/cM74LwZ6OudOcc6N9gclfTlst5wzkCQz7ntn\nCdc/NTPocEREEkI4iX8uUBDpQCTxdMjP5F/fH0Gngky2lFdStndf0CGJiMS9cBJ/AbDQzN7S7XzS\n0noX53DPJUcB8O6iLQFHIyIS/8K5xv+riEchCe3Y7oW0z8vglU/Xce5RHYMOR0QkrjWZ+J1z77ZG\nIJK4kpKMcUM68ujUFWzYuYcO+eoRWkQkUsLpua/czMr8Ya+Z1ZhZWWsEJ4njoqGdAbj8kY8p17V+\nEZGIaTLxO+dynXN5zrk8vIfzfBV4IOKRSULpU5LLY1cfy7Itu3hq2uqgwxERiVvhNO47wHleAc6M\nUDySwEb2bceRnfN5a97GoEMREYlbTV7jN7MLQ14mAcOAvRGLSBLaGQNL+MPbi9lUtpeSvIygwxER\niTvhnPGfGzKcCZQD4yIZlCSuMwa1B2DS/E0BRyIiEp/CadV/TWsEIgLQpziHHkXZvD1/E5cP7xZ0\nOCIicafBxG9mv2xkPeecuy0C8UiCMzPOGFjCYx+soGzvPpyD/33hc356Rj/6lOQGHZ6ISMxrrKp/\nVz0DwDeAGyMclySwMwaVsK/GMWXRFv71+XremreJrz08jZ27dZufiMjhavCM3zl39/5xM8sFfghc\nA0wE7m5oPZHDdXSXNpTkpTPhgxWU760GoGzPPs69fyr//sEIcjNSA45QRCR2Ndq4z8wKzew3wGy8\nHwnHOOdudM5tbpXoJCElJRnjR/Zi1uodLNlcwePXHMuEa45j9fbdnPPnqergR0TkMDSY+M3s98An\neK34j3DO3eqcK221yCShff34rgB0zM9gdL9iRvQp4u6Lj2LVtt08+O6ygKMTEYldjbXq/wlQCfwC\n+LmZ7Z9ueI378iIcmySwjNRkPvzZqaSnfPHb9KtDO/Phsm08MGUZSzZVcGr/Yi4e1oXkJGtkSyIi\nEqqxa/yH1KufSEvrWPDlh/XcfsFgctKTeeKjVbw9fxOfr93BnRceGUB0IiKxScldYkpGajK/HjeY\nZXeMZfzInjw7fQ3TV2wPOiwRkZihxC8xKTnJ+NGYvrTPy+CWV+ayuVy9SIuIhEOJX2JWZloyvzl/\nMCu37eK6J2dSU+uCDklEJOop8UtMGzOwhNsvOILP1+zgw2Vbgw5HRCTqKfFLzDv3qA7kZqTw+pwN\nQYciIhL1lPgl5qWnJNO/fS7Lt+xqemERkQSnxC9xoWthNqu27Q46DBGRqKfEL3Ghe9ssNpbtZe++\nmqBDERGJakr8Ehd6FecAMGu1epUWEWmMEr/EhVP7F9MmK5WnP14ddCgiIlEtkMRvZivNbI6ZfWZm\nM/xphWY2ycyW+H/bhCx/k5ktNbNFZnZmEDFLdMtITeaMge15b/EWqmtq2VS2l5mbqqnVvf0iIgcJ\n8ox/tHNuiHNumP/6Z8A7zrk+wDv+a8xsIHApMAg4C3jAzJKDCFii2yn92lG+t5rfv7WIayd8wp8/\nreSBKUuDDktEJKpEU1X/OOAJf/wJ4PyQ6ROdc5XOuRXAUuC4AOKTKHfGwBLOGtSeh95bzrz1ZQA8\n/fFqqmtqA45MRCR6mHOtXxVqZiuAnUAN8JBz7mEz2+GcK/DnG1DqnCsws/uBac65p/x5jwJvOOde\nqGe744HxACUlJUMnTpzYYjFXVFSQk5PTYtuLJvFUtp2Vjgc+28vAtsmUpFXx4Hzjsv5pnNk9NejQ\nWlQ8vWeh4rVcoLLFolgr1+jRo2eG1KI3qMHH8kbYCOfcOjMrBiaZ2cLQmc45Z2aH/IvEOfcw8DDA\nsGHD3KhRo1okWIApU6bQktuLJvFWtnF+K5DJkyezqCabV5dv59TjjmTMgGK835SxL97es/3itVyg\nssWieC1XIFX9zrl1/t/NwMt4VfebzKwDgP93s7/4OqBLyOqd/WkijTIzfnPBEXQpzOK6J2dw3v0f\nsGrbLmprHUHUdImIRINWT/xmlm1mufvHgTOAucCrwFX+YlcB//THXwUuNbN0M+sB9AGmt27UEqs6\nFWTy6vdGcMcFR7CmdDen/H4K/W55gwffXR50aCIigQjijL8EmGpmn+Ml8H87594E7gJON7MlwBj/\nNc65ecDzwHzgTeC7zjl1zyZhS0tJ4v8d35UnrjmOjNQkzIxHpy5nzXZ18SsiiafVr/E755YDR9Uz\nfRtwWgPr3A7cHuHQJM4d1aWAT34+hjnrdnL1459w2t3v8sK3T+DIzgVBhyYi0mqi6XY+kYjLzUjl\nxF5FvHXDSLLTk/nNvxfoer+IJBQlfklIPYqy+fEZ/Zi+YjufrFT//iKSOJT4JWGdd1RHzGDa8m1B\nhyIi0mqU+CVh5Wem0q8kl09Wbg86FBGRVqPELwnt2O6FzFpVqm59RSRhKPFLQhvWvQ27qmpYuLEc\ngKWbK/homar+RSR+KfFLQhvesy1mMGn+Jnbu2cc5f36fy/42jfcWbwk6NBGRiFDil4RWkpfBiN5F\nPPbBCsbdP5W9+7wq/5tfnsPLn67VrX4iEneU+CXh/fq8QXQtzGL9jr3cccERTBw/nMzUZH703Odc\n9+RMqmtq2VdTS0VlddChiogctqCezicSNXq2y+HV742gfO8+CrLSAHjrhpE8MnU5d7y+kBtfnMOn\na0pZs303N4zpy3Un9yQtRb+ZRSQ2KfGLAMlJdiDpAyQlGeNH9mLltt088/FqzKBH22x+/9Yi5m8o\n4y//75gAoxURaT4lfpFG3DZuMMd1L6RvSS4DO+Zx3ztLuGfSYi4aupnR/YqDDk9E5JCpvlKkEclJ\nxvlHd2JgxzwArj+lFz2LsrntX/Opqta9/yISe5T4RQ5BWkoSt5wzkOVbd/HEhyuDDkdE5JAp8Ysc\notH9ixndrx33vbOEvftqWLq5gtdmr1fvfyISE5T4RZrhyhO6U15ZzYfLtnL5Ix/zvWc+ZYJqAEQk\nBijxizTD8T0LSUtJ4toJM9hYtheAxz9YqQ5/RCTqKfGLNENWWgo3jOkDQPe2Wdx2/mDW7djD8q27\nAo5MRKRxup1PpJm+fUovBnXMp1NBJhmpSZjBOfdN5eavDOCK4d2CDk9EpF464xdpJjPjlL7t6F2c\nQ+c2Wfx87ACv1f8rc7nyseks21IRdIgiIl+ixC/SQr55ck9m/GIMlx3XhZkrt3PWn95j7L3v67Y/\nEYkqSvwiLSg1OYk7LzySyf8ziqtP7E6tc/zq1XncM2kxldU1QYcnIqJr/CKRUJybwc+/MpCKymq+\n98ws7ntnCZ+uLuWxq48lNVm/t0UkOPoGEomgnPQUJlxzHL/96hG8v2QrJ9z5X9bt2BN0WCKSwJT4\nRVrB147tyl0XHsHWikqueXw6m8v38va8jdz7nyVBhyYiCUZV/SKt5NLjulKSn8G3n5rJaX94l/LK\nagBO6deOIV0KAo5ORBKFzvhFWtHofsU8dtWxjOpfzFUndCPJ4LlPVgcdlogkEJ3xi7SyE3sXcWLv\nIgBKd+/jjbkbuW3cYFLU6E9EWoG+aUQCdOag9uzYvY+Zq0qDDkVEEoQSv0iARvhn/l97eBo3TPyU\nPVW6119EIktV/SIBys9K5fYLBvP5mh38Y+Zapi3fzrDubaiqruXiYV3okJ/B4E75QYcpInFEiV8k\nYF8/vhtfP74bx3Yv5LdvLuS12RsAeHv+JgB+ekZfvndqnyBDFJE4osQvEiUuHtaFswZ71/wLs9NY\nuLGcp6at4g9vL2ZQx3xG9y8OOkQRiQNK/CJRJDcjldyMVACGdmvDkZ3zmbFqO/dMWsyofu0ws4Aj\nFJFYp8QvEsVSk5P4wal9+J8XZvOv2Rv4cOlWPly2jRHF1YwKOjgRiUmt3qrfzLqY2WQzm29m88zs\nh/70W81snZl95g9jQ9a5ycyWmtkiMzuztWMWCdIFR3eiR1E2P3j2U56bsYZa53h2YRVz1+0MOjQR\niUFB3M5XDfzEOTcQGA5818wG+vP+6Jwb4g+vA/jzLgUGAWcBD5hZcgBxiwQiJTmJP192NF8b1oXH\nrjqW1394MlmpcMlDH/Hwe8vYu0+3AIpI+Fo98TvnNjjnZvnj5cACoFMjq4wDJjrnKp1zK4ClwHGR\nj1QkegzulM9vLzqS0f2LyctI5SfDMhjSpYA7Xl/IFY9+zKayvUGHKCIxItAOfMysO3A08LE/6ftm\nNtvMHjOzNv60TsCakNXW0vgPBZG41zM/mWeuG869lw5h/voyTv3DFI769dv8cOKnfLB0K9U1tUGH\nKCJRypxzwezYLAd4F7jdOfeSmZUAWwEH3AZ0cM5da2b3A9Occ0/56z0KvOGce6GebY4HxgOUlJQM\nnThxYovFW1FRQU5OTottL5rEa9nitVxwcNlW7qzh9RX72L0P5m7zqv1zUuGSfmmM7JwaZJiHLFHe\ns3gTr2WLtXKNHj16pnNuWJMLOudafQBSgbeAHzcwvzsw1x+/CbgpZN5bwAlN7WPo0KGuJU2ePLlF\ntxdN4rVs8Vou5xou27aKSvfGnA3u4gc/dN1ufM19sGRL6wZ2mBLxPYsH8Vq2WCsXMMOFkYODaNVv\nwKPAAufcPSHTO4QsdgEw1x9/FbjUzNLNrAfQB5jeWvGKxJLC7DTOGtyeJ689ji6FmdzxxgKqqmv5\n8fOfcdNLc6ipDaaGT0SiRxD38Z8EXAHMMbPP/Gk3A5eZ2RC8qv6VwLcAnHPzzOx5YD7eHQHfdc6p\nGbNIIzJSk7nxrP5875lPGfqbSZTvrQZgeM9Cxg1RExmRRNbqid85NxWor/ux1xtZ53bg9ogFJRKH\nzjmyIxt37uWteRu54OjO/GXyUm57bQF5Gan0apdDl8JM9QQokoDUc59IHPvmyT355sk9AUhPSeIn\n//icayZ8AsB5R3XkB6f1ISM1iTlrd9KnJJfexbHTkElEmkeJXyRBfHVoZ3q2y2bJpgrmbyjjiY9W\n8urn6w/MTzL4yRn9+M6oXqoJEIljSvwiCeTorm04uqvXRcb1p/TimY9XUV5ZzRkD2/PzV+bw+7cW\n0b1tNl85skMTWxKRWBVoBz4iEpz2+Rn8+Ix+/OrcQZzQqy2TfnQKfYpz+MvkpTjneOT95Yy7fypv\nzt24/1ZaEYkDSvwiAkByknHdyJ7M31DGT/8xm9tfX8Dna3dy/VMzGXvfVP752Tr9ABCJA6rqF5ED\nLjqmMy/NWsuLs9bStTCLl75zIne+vpAXZ63lhxM/497/LOGiYZ1JS/bOGXZX1XBs90KO71FIUpLa\nBYjEAiV+ETkgKcm455IhPDZ1BVed2J2inHT+cPGR/HrcIJ75eBUTP1nD795c9KX1ju9RyF8vH0ph\ndloAUYvIoVDiF5GDdCzI5BfnDDzw2szISU9h/MhefHNET9bt2HNgXkZqMo+8v5yH3lvOdU/O4JqT\nunPfO0sozE7jocuHkZ8VW88KEEkESvwiErakJKNLYdZB024aO4Ds9BTumbSYmatK6dY2i2nLt3Pq\n3VMozsuga2EmPz69H/3a5wYUtYiEUuIXkcN29UndmbtuJ6kpSfzhoqOYuaqUpz9exb6aWqav2M5l\nf5vGpB+NpCArjY1le+mYn6G+AkQCosQvIoctLyOVh6/84mmgI/oUMaJPEQDz15cx7i9TOfNP79E2\nO51Fm8rpVJDJoI55jOhTxJUndGfnnn1UVdcCsLlsL5lpyeRm6DKBSCQo8YtIRA3smMfL3zmJO15f\nwN59NXxrZE+mLd/G2/M38fb8Tbzy6To+X7uTmlpHWjJUvfkOKUnGeUM68oeLjqKqppYnP1rJ2tI9\nnDWoPSf2Lgq6SCIxTYlfRCJucKd8nrlu+EHTKqtrGHvv+8xavYOzB7dnSJcCZi1YyrGD+rBsyy6e\nnb6aGStLKchKZfbanWSkJvHUtFXcd9nRnHNkx4BKIhL7lPhFJBDpKck8+Y3jmb++jDEDijEzprg1\njDq5J845uhZm8dwnq1lbuod7Lx3C6QNLuOqx6fz4uc/JTE3mtAElVNfUkmSmPgREDoESv4gEplNB\nJp0KMr803cz49qhefHtUr4OmP3zFMC772zS+8cQMehfnsGrbLkb0LuLxa45rrZBFYp667BWRmNEm\nO40Xv30iw3sWsquymn01jsmLtrBsSwW1tepOWCQcOuMXkZiSnZ7CxPEnAPDJyu1c/OBHnHb3u/Qp\nzuGvlx9Dr3Y5ulVQpBFK/CISs4Z1a8ODlw9l/Y49/PE/ixlzz3v0b5/LWYPb07s4h6y0ZPbuq2V4\nz7bqTljEp8QvIjHLzDhrcHsAzhhUwgsz1/LGnI386T9LvrRsarKRlZZCm6xUzhzUnm+e3JN2uemt\nHbJI4JT4RSQudG6TxQ1j+nLDmL6s2b6bPftq2LBzL0kG01dsZ1+NY/2OPZTuruKRqSt4+dN1XHhM\nZ75yRAcGd8rT5QFJGEr8IhJ39j9PoG+J93yAk/u0O2j+rNWl/HHSYv72/nIefHcZndtk8q2RPbni\nhO6tHapIq1PiF5GEc0zXNvz9G8dTuquKSQs28Y8Za7jln/P4x8y1ZKQkA3B01wK6tc2mQ34GQ7u3\nIcXvKyA5yUj3lxGJRUr8IpKw2mSnccmwLowb0pGbX5rLvPU7yU5LoaKymofeW17vOslJxtmD2/OV\nIzpwXI9C2uak45yjbG81eRkpumQgUU+JX0QSXnpKMndfctRB07bvqqKqupYFG8tYvLH8wPRNZZU8\n/fEqXpu9gey0ZC4e1oXV23fz34Wbyc9Mpaq6lqHd2tCxIIPtu6rIy0jltAElrN1ew7GV1WSnf/G1\nW1PrqHWOir3VAKSlJJGdnkJVdS1pKepmRSJDiV9EpB77b/9rn5/B6H7FB8370el9mL++jLveXMgz\n01djwLUn9WBrRSX5mal8vGIby7ZUkJaSROmuKl76dB0Ad8+cRNe2WfRul0Ov4myen7GWLeWVB7ab\nZNA2J52tFZUc07UNpw0opkN+BpmpyZzcp91BPxpEmkufIhGRQ5SbkcrxPdvy8ndOorK6BucgI7X+\n6/57qmp4Z+Emliycz56cTizeVM6cdTt5c95GhnZrw3lHdaRDfgYpSca2XVVsKa8kIzWZN+du5Hdv\nLjqwnZ7tsvnqMZ1ZW7qHNdt38+GyraSlJLF3X+2BZfIzUzmhZ1vSU5PITE3mmyf3pHdxTsSPh8QW\nJX4RkcPQVEO/zLRkzjmyI1O2L2bUqAEAOOco3b2v0U6Fbj1vENsqKinfW83CjeXc+uo8fv/WInLT\nU2iXm85FQzvjHHQoyGR/q4LV23fz2ZodOOfYWLaXf8/ewE/P7EfXwizW7djDzj376FuSy+y1O+iQ\nn0mfkhzyMlJxOPoW5+phRwlCiV9EpJWZWVg9CbbNSadtTjrdi7I5c1AJldW1pCUnhZWg15bu5gfP\nfsqvXp0XVkw9i7I5pV87vnpMZwZ3yg9rHYlNSvwiIjHAzBq8nFCfzm2yeOH6E3n503UkJcHxPdpS\nWV3L+h17OLJzPlsrqli6uYIt5ZWYwcTpq3nm49U8NW0V147owfdP7UN2WjIvzlrHPz9bx8AOeXTI\nzziw/Q4FmeRnpgKwcHsNbdfu5P+3d6cxdpV1HMe/v7kzd9ZOp9NlCi3QKVSgGixQyippBRXQWAm+\nKFHDCw0xEYX4wkAwRn2lxri8QAxBjIlKE0GQEAybFiIuhUKrpWydWmhLy0yHbjOdff6+OKftnYFO\nF+8IUfIAAAnrSURBVIfennN/n+TkPue5Z26e32z/e5+znTWrifqiT3U82bnwm5nlVFWVuOHCuWP6\n2mc0AslxCgfaADcuOZ1dvYN8++H13PPsJh56cRszp9Ty8tt7OXVqHX/v6GZkojsgrv4bxeoqrl80\nh5bGGjr3DjClrpqtu/oYjeCUqXW0NhZpn9FETUFcMn86bc11Y15ieGSUrp4BZjfX+bTID5ALv5mZ\nAcl1De76wgWsefNd7l7VQde+AW69agG3XrWAgeFR+odGABiNYHN3LwPDyYGF69au49Qzz+F3/3qL\nh9YmZzDU1xQYGB7htGkN1NUUWLN5F72Dw5S+d6itruLSM6czraFI174BnuvYSaTPz26uo6WhhmkN\nRWanMw3FQhXzZjRSW13FJz/cxuad+xmJ938zUl9ToKWhhghoaahBwIym2oO7SYZHRunuHUy2LRZo\nrquZ7G/nScuF38zMxrjwjFbuval1TF99sTBmGn9606EbHA1uKbB00RyWL5oz4esODI+wbVcfu/uG\neO6NnXT1DPCPjm42dfVSJVhx0emc3trA/sFhNnfvp39ohO17+tj2Zh8Ae/qG2NM3BMD3H91wzLna\nmmu5/vy57Ood5IkNO9i1P3ktKTkjYk5LPWfNamJOSz31NQX6Oodp2bKbYqGK+TMb2dkzwN87uon0\nzUZNoYr2GY1USUjJJaKPZXdMubjwm5nZCVFbXWD+zOT0wgtOn3bMXz86GvQPj7B+217+vH47F7e3\nMnNK3ftu+27vILv2DyJgaCQYGhnl6Vc7+eUzHTTVVnP1ubO4cF4rBYnOff3s7Bngvzt7+UdHN50l\n11b4xbrnjnp8U+trOG/uVKbUVY/rL3L1ubO49MzpNBTLX3bLPwIzM7OjUFWV3Fp5SXsrS9pbj/wF\n49x02Ty69g3QXF894WmY+/qHGA148PFnmXfOR9jbN8y23X0UqsQVZ804eEbG3v4h3t6dzEbs6h3i\nmde7eP2dfezY03/wtQJ4q3s/969+i+oqUayuYlpDkVnNtXzmvFP58hXtx5zj/+XCb2ZmFWPmlNoj\nbjMl3d8/v6XA0nPaDrvdqdRzzuzmg+vjD6Q84N3eQV5+ew//3NTNwNAob+/pY1//MLVluiyzC7+Z\nmdkHqLWxyMcWzHzP7aHLJTN3gZB0jaTXJG2UdHu5x2NmZpZFmSj8kgrAXcC1wELgRkkLyzsqMzOz\n7MlE4QeWABsjYlNEDAIrgeVlHpOZmVnmZKXwzwG2lKxvTfvMzMzsGCgOc9Wjk4mkzwPXRMRX0vUv\nARdHxC3jtrsZuBmgra3twpUrV07aGHp6emhqyuftLfOaLa+5IL/Z8poLnC2LspZr2bJlayJi8ZG2\ny8pR/duA00rW56Z9Y0TEPcA9AIsXL46lS5dO2gBWrVrFZL7eySSv2fKaC/KbLa+5wNmyKK+5sjLV\n/zywQFK7pCKwAnikzGMyMzPLnEx84o+IYUm3AI8DBeC+iDi6m0ybmZnZQZko/AAR8RjwWLnHYWZm\nlmVZmeo3MzOzSeDCb2ZmVkFc+M3MzCqIC7+ZmVkFycQFfI6HpC7gzUl8yRnAzkl8vZNJXrPlNRfk\nN1tec4GzZVHWcp0REUe8BWBuC/9kk/TC0VwRKYvymi2vuSC/2fKaC5wti/Kay1P9ZmZmFcSF38zM\nrIK48B+9e8o9gA9QXrPlNRfkN1tec4GzZVEuc3kfv5mZWQXxJ34zM7MK4sJ/BJKukfSapI2Sbi/3\neI6VpPskdUpaX9LXKulJSW+kj9NKnrsjzfqapE+VZ9RHJuk0SX+VtEHSy5JuTfvzkK1O0mpJ69Js\n30v7M58NQFJB0kuSHk3X85Jrs6T/SFor6YW0Ly/ZWiQ9IOlVSa9IujQP2SSdnf68Dix7Jd2Wh2wT\niggvh1lI7gTYAcwHisA6YGG5x3WMGa4ELgDWl/T9CLg9bd8O/DBtL0wz1gLtafZCuTMcJtcpwAVp\newrwejr+PGQT0JS2a4B/AZfkIVs63m8CvwcezcvvYzrezcCMcX15yfYb4Ctpuwi05CVbScYCsAM4\nI2/Zxi/+xD+xJcDGiNgUEYPASmB5mcd0TCLiWeDdcd3LSf6QSR8/V9K/MiIGIuK/wEaS78FJJyK2\nR8SLaXsf8Aowh3xki4joSVdr0iXIQTZJc4FPA/eWdGc+1wQyn03SVJIPEL8CiIjBiNhNDrKNcxXQ\nERFvkr9sY7jwT2wOsKVkfWval3VtEbE9be8A2tJ2JvNKmgecT/LJOBfZ0unwtUAn8GRE5CXbz4Bv\nAaMlfXnIBcmbs6ckrZF0c9qXh2ztQBfw63QXzb2SGslHtlIrgPvTdt6yjeHCX+Eimb/K7KkdkpqA\nB4HbImJv6XNZzhYRIxGxCJgLLJH0kXHPZy6bpM8AnRGx5nDbZDFXiSvSn9m1wNckXVn6ZIazVZPs\nLrw7Is4Hekmmvw/KcDYAJBWBzwJ/GP9c1rO9Hxf+iW0DTitZn5v2Zd07kk4BSB870/5M5ZVUQ1L0\nfxcRf0y7c5HtgHRK9a/ANWQ/2+XAZyVtJtlt9nFJvyX7uQCIiG3pYyfwEMkUcB6ybQW2prNOAA+Q\nvBHIQ7YDrgVejIh30vU8ZXsPF/6JPQ8skNSeviNcATxS5jFNhkeAm9L2TcCfSvpXSKqV1A4sAFaX\nYXxHJEkk+xxfiYiflDyVh2wzJbWk7XrgE8CrZDxbRNwREXMjYh7J39JfIuKLZDwXgKRGSVMOtIFP\nAuvJQbaI2AFskXR22nUVsIEcZCtxI4em+SFf2d6r3EcXnuwLcB3JEeMdwJ3lHs9xjP9+YDswRPLO\n/cvAdOBp4A3gKaC1ZPs706yvAdeWe/wT5LqCZPrt38DadLkuJ9nOA15Ks60HvpP2Zz5byXiXcuio\n/sznIjnzZ126vHzgf0UesqVjXQS8kP5OPgxMy1G2RqAbmFrSl4tsh1t85T4zM7MK4ql+MzOzCuLC\nb2ZmVkFc+M3MzCqIC7+ZmVkFceE3MzOrINXlHoCZnbwkHTitCWA2MEJy+VaA/RFxWVkGZmbHzafz\nmdlRkfRdoCciflzusZjZ8fNUv5kdF0k96eNSSc9I+pOkTZJ+IOkLklan96c/M91upqQHJT2fLpeX\nN4FZZXLhN7PJ8FHgq8C5wJeAD0XEEpLb73493ebnwE8j4iLgBsbemtfMThDv4zezyfB8pLcxldQB\nPJH2/wdYlravBhYmt1kAoFlSU0T0nNCRmlU4F34zmwwDJe3RkvVRDv2fqQIuiYj+EzkwMxvLU/1m\ndqI8waFpfyQtKuNYzCqWC7+ZnSjfABZL+rekDSTHBJjZCebT+czMzCqIP/GbmZlVEBd+MzOzCuLC\nb2ZmVkFc+M3MzCqIC7+ZmVkFceE3MzOrIC78ZmZmFcSF38zMrIL8D2+kE8hpjA2PAAAAAElFTkSu\nQmCC\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": 16, "metadata": { "collapsed": true }, "outputs": [], "source": [ "index_shallow_particles=[]\n", "for n in range(2694):\n", " if result.variables['init_z'][n] < 6:\n", " index_shallow_particles.append(n)" ] }, { "cell_type": "code", "execution_count": 19, "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": 21, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf4AAAGDCAYAAADK03I6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XmcHHWd//HXZ7rnvmcSJpNMyM2RcJoQCJfhkkMEdBVx\nF0VWQF1cYHVXQddb8NiVH7iKiqKgoBFRkWU5BCREkBBIuHKQg9z3ncxMkjk/vz+qJnQmc/RMZqam\nu9/Px6MfU11VXfX5dPf0p7/f+la1uTsiIiKSGbKiDkBEREQGjgq/iIhIBlHhFxERySAq/CIiIhlE\nhV9ERCSDqPCLiIhkEBV+wczuNbNvRbRvM7NfmtkOM5sTUQx1Zja2m3VGm5mbWbwf9v+CmZ0YTn/R\nzH6e5OOSXreTx083s7W9fXw32+6356sHMVSZ2SIzy+1inZ+Y2ZcHMq7+ksz7eKCY2RlmtjjqOKRj\nKvyDkJmtNLPNZlaYMO8aM5sZYVj95XTgPKDG3ae2X2hmHzezlvBDbbeZvWZmF/d2Z2Y208yuSZzn\n7kXuvry32zwUZvY+oNbdXw1juc3dr+nmYbRfN+pCG75nz41o3+ea2TwzqzeztWZ2OYC7bwKeBa7r\n7LHu/il3/2Yv9zvTzPaZWW343pxrZjd39UWjP/XX+ziZhkH43hufEMvf3P3Ivo5F+oYK/+AVA26M\nOoieMrNYDx8yCljp7vVdrPOiuxcBZcA9wINmVt7DuMzMBuP7/VPAr6MOYjAzs6oulk0EfgN8CSgF\njgfmJqzyAPDJfgzvM+5eDFQDnwOuAB4zM+vHfYocGnfXbZDdgJXAzcB2oCycdw0wM5weDTgQT3jM\nTOCacPrjwAvA/wN2AsuBU8P5a4DNwFUJj70X+AnwFFALPAeMSlh+VLhsO7AYuLzdY38MPAbUA+d2\nkM9w4JHw8cuAa8P5nwD2AS1AHfD1Dh77ceD5hPuFYe5TgHLgUWALsCOcrmn3nNwaPhd7CYpAS7jP\nOuCH4XoOjA+n84HvA6uAXcDz4bwDnnOCInMPsAFYB3wLiIXLxofP4S5gK/C7Tl7nnDCuxJi/Btzf\n7nW+ClgdbutLnay7Oly3LrxN62B/+eHrtQNYCPwHsLbd6/SH8PlcAdzQbl8PAb8L3yPzgOPDZb8G\nWsNc6oDPdxd7Ev8D2cD7w/dNbRfr/Qb4ZhfL48AeEt7P7ZbfC3wrnJ4OrCUo4JvD1/bqLrY9k/B/\nLmHe4eH+Lg7v5wJ3AOvD2x1Abrv9fT5hf5cBFwFLCP5fvpiw7anAiwT/0xuAHwI5CcsT38f3Aj8C\n/i98vV4CxnWRy++BjQTv2VnApHD+dUAT0Bi+tv/bwWNnhfuuD9f5cFtu7T7T/gN4I1zvHqAKeDyM\n72mgPGH9U4C/h7m+Dkzvq8/XcPunhu/JkeH94wn+L47qy/0M1lvkAejWwYsS/JOcC/wx4UOpp4W/\nGbiaoOfgWwQfvj8KP4jeE/6zFYXr3xvePzNcfidhsSUotGvCbcWBE8N/mIkJj90FnEbQg5TXQT6z\ngLuAPOAEgsJydkKsz3fxXHw8IZY4QS9ILUHhrQT+ASgAisMPr4fbPSergUnhY7Pp+MM68QPzR+E6\nI8Ln7tTwOTngOQf+BPw0fH4OA+YAnwyX/ZagBZoV5nx6J7lNAurbzfsaBxf+nxEU7eOBBuDoLtaN\nd/Fcfgf4G1ABjATmE344h7HOBb5C8IVkLMEXxvMT9tUEfDB8Hv+d4MtBduJ7NmFfXcbeRYzHArcT\nFMIXCVrrZV2svxz4JvAmQTG8H6hot84bwCWdPP5eDiz8zcA3whwvIiji5Z089qD3UsL7/bvh9DeA\n2eF7ZChBMftmu/19JdzftQT/G78heD9PIvgyNSZcfzJBQYyHz+8i4KZO3sf3AtsIvizECb70zuji\nefzncJ9tX1Re6+g56uLx+/edkFv7wj+boNiPCF/feQSfJ3nAX4GvhuuOCGO/iOB9eV54f2gn+36U\n4AtCR7dHu4j51nC/+eH75zPJfD6nwy3yAHTr4EV5p/AfQ1BUh9Lzwr80Ydmx4fpVCfO2ASeE0/cm\nfigARQQt45EE397/1i6+nyb8k94L/KqLXEaG2ypOmPdt4N6EWLsr/M3hP/HW8MPjoF6FcN0TgB3t\nnpNvtFtn//OUMM8JWulZBB+0x3ew7f3Pefjh1QDkJyz/CPBsOP0r4G4SWvKdxHsasLHdvK9xcDFP\n7BGYA1zRxbpdFf7lwAUJ96/jncJ/MrC63fq3AL9M2NfshGVZBIX2jMT3bAfPV4exdxDb2cArBF8y\nbwOOSPJ/pTHc9xEE79s/AA+0W+cF4GOdPP5eDiz8eznw/2ozcEonjz3ovRTOnwH8LJx+G7goYdn5\nBIe2EvfX1lNUHD5nJyesPxe4rJP93wT8qf37OCGvnycsuwh4K8nntCzcVmn756iLxyRT+P8p4f4f\ngB8n3P9Xwi/twBeAX7fb/pMk9FL2xY3gy9ZcgqL/BGB9uf3BfItsxK10z93nm9mjBN3+i3r48E0J\n03vD7bWfV5Rwf03CfuvMbDtB1+8o4GQz25mwbpwDj0uvoXPDge3uXpswbxVBV32yZrv76e1nmlkB\nweGMCwi6/QGKzSzm7i1JxNbeEILWx9vdrDeK4ENjQ8Kh3KyEfX2eoBU6x8x2AN939190sJ0dBB/2\n3dmYML2HA1+3nhjOgc/HqoTpUcDwdq9zjKCHoE3ie6Q1PCNgeDf7TDb2wwi+fL1I0LW7upvtttlL\n8OVkCYCZ3UbQbZyomOCLYzK2uXtzkjF3ZgRByx6C5yfxeV7Fgc/ZtoT36t7wb4f/p2Z2BEFvyBSC\nXq44B45naC+p5z4cl3Mr8CGCRkZruGgIQcOjr7TPq7PPo1HAh8KBr22yCQZq9hl3bzKze4EfAJ/1\n8NtAJhiMg53kQF8l6AIckTCvbSBcQcK8YYe4n5FtE2ZWRNAdvJ7gw/45dy9LuBW5+6cTHtvVP8x6\noMLMEgvc4QTHxQ/V54AjCVpIJQSHKgASB1a1j62rWLcSHP8f181+1xC0+IckPCcl7j4JwN03uvu1\n7j6coKv6rsQRzwmWEYw7HNHBsp5K5kNrAwmvM8Hr0GYNsKLd61zs7hclrJP4HskCaghe32T33yl3\nn0HwHv41wdiP9Wb2MzM76AtfO2+02/cBcYRnOYwn+DLR78xsJEGXfNsXpvUEhazN4bzznPXUj4G3\ngAnh+/2LHPhe761/BC4l6GUsJeitIWHbA10Q1xC0+BPfi4Xu/p2OVjazx8Ozfjq6Pd7ZTsL/u68C\nvwS+H9XZGFFQ4R/k3H0ZwYCqGxLmbSEonFeaWczM/pnui1V3LjKz080sh6C1Otvd1xAcPzvCzD5q\nZtnh7SQzOzrJ+NcQtH6+bWZ5ZnYcwQf7/YcYLwQtub3ATjOrIPgn7s4mguPXHcXaCvwCuN3MhofP\n7bT2HwjuvgH4C8GHRYmZZZnZODN7N4CZfcjMasLVdxB8cLbSjrs3ErRO351Mst3YEu6jq/O4HwRu\nMbPyML5/TVg2B6g1sy+YWX6Y+zFmdlLCOpPN7ANhMb2J4MvP7HBZp89rstx9n7v/xt3fQzAmYCXw\nSzPrqgfml8DVZjY27AG6meA922YqQdf6qg4f3UfMrCB8/f9M8Fw+Fi76LfCfZjbUzIYQHM/v7Xu/\nGNgN1JnZUcCnu1m/J9ttIDj8V0BwqCVRMq/tIb/+Ce4H3mdm54fvwzwLrjlR09HK7n5h2Bjp6HZh\nR48Jz7q4l2CQ4ScIvhT36rTOVKTCnxq+QTCILNG1BKNktxEMAvp7+wf10G8ICud2ghbLlQBhF/17\nCE5TWk/QffhdgkFAyfoIQStiPcGguK+6e/vu2N64g2BgTtux/yeSeMydwActuGDQDzpY/u8Ex/xe\nJnguvkvH/ycfIxgEt5CguD9EcEoXwEnAS2ZWRzAq/Ubv/PzqnwIfTSLuLrn7HsIzGMxsp5md0sFq\nXyfoal5B8MXl1wmPbwEuJhgnsYLgOf05QQuwzZ8JxnzsCGP+gLs3hcu+TVDgdprZv/dBPmvc/VZ3\nn0BwZkBn6/2CYEzFS2FuDSR8SQb+ieCMlf7yQzOrJSh8dxAcu74g/BIJwcDaVwh6Jt4kGNDW24tl\n/TtB67yWYNDk7w4h7kS/Inju1hG8n2e3W34PMDF8bR/uZBtfA+4L17n8UIIJGwuXEvRobCHoAfgP\n+rZe3UBweOnLYRf/1QRfIM/ow30MWpZBhzVEBiUze4FgRPGrUcfSGTP7GsHgrSujjiVZZnYYwWmV\nJ7r7vqjjERksNLhPJGLuflrUMaQjd98MJHVISiSTqKtfREQkg6irX0REJIOoxS8iIpJBVPhFREQy\nSNoO7hsyZIiPHj26z7ZXX19PYWH7M+rSQ7rmlq55Qfrmlq55gXJLRamW19y5c7e6+9Du1kvbwj96\n9GheeeWVPtvezJkzmT59ep9tbzBJ19zSNS9I39zSNS9Qbqko1fIys6QuVKWufhERkQyiwi8iIpJB\nVPhFREQyiAq/iIhIBlHhFxERySAq/CIiIhlEhV9ERCSDqPCLiIhkEBV+ERGRDNJvhd/MfmFmm81s\nfsK8CjN7ysyWhn/LE5bdYmbLzGyxmZ2fMH+ymb0ZLvuBmVl/xSwiIpLu+rPFfy9wQbt5NwPPuPsE\n4JnwPmY2EbgCmBQ+5i4zi4WP+TFwLTAhvLXfpoiIiCSp3wq/u88CtrebfSlwXzh9H3BZwvwZ7t7g\n7iuAZcBUM6sGStx9trs78KuExwyYeat3sHRHC0EIIiIiqWugf6Snyt03hNMbgapwegQwO2G9teG8\npnC6/fwOmdl1wHUAVVVVzJw5s0+CvnX2XpbubGXFrqd5z+jsPtnmYFJXV9dnz9Vgkq55Qfrmlq55\ngXJLRemaV2S/zufubmZ92oR297uBuwGmTJniffWrSiOOruXae57nyTXw9Y+eSXYsvcZEptovUCUr\nXfOC9M0tXfMC5ZaK0jWvga5gm8Lue8K/m8P564CRCevVhPPWhdPt5w+oCVXFXDo+m231jTy/bOtA\n715ERKTPDHThfwS4Kpy+CvhzwvwrzCzXzMYQDOKbEx4W2G1mp4Sj+T+W8JgBdXRFMNZw4frdUexe\nRESkT/RbV7+Z/RaYDgwxs7XAV4HvAA+a2SeAVcDlAO6+wMweBBYCzcD17t4SbupfCM4QyAceD28D\nLi9uVBbmsHbHnih2LyIi0if6rfC7+0c6WXROJ+vfCtzawfxXgGP6MLReq6koYM32vVGHISIi0mvp\nNUqtn40sz2eNWvwiIpLCVPh7YEhRLtvrGqMOQ0REpNdU+HsgLzvGvuaW7lcUEREZpFT4eyAvO4um\nFqelVVfwExGR1KTC3wN52cEpffua1OoXEZHUpMLfA3nx4OlS4RcRkVSlwt8DbS3+hubWiCMRERHp\nHRX+HlBXv4iIpDoV/h7Iy27r6leLX0REUpMKfw/kxsMWv07pExGRFKXC3wO52RrcJyIiqU2Fvwf2\nD+5TV7+IiKQoFf4eyIu3jepXi19ERFKTCn8PaHCfiIikOhX+HtDpfCIikupU+HtAhV9ERFKdCn8P\n5LZdsldX7hMRkRSlwt8DbYVfo/pFRCRVqfD3QDyWRSzLaGxRV7+IiKQmFf4eyoll0aiufhERSVEq\n/D2UE1fhFxGR1KXC30M58SwaW1T4RUQkNanw91BOLIsGtfhFRCRFqfD3UG62Cr+IiKQuFf4e0uA+\nERFJZSr8PZSrwX0iIpLCVPh7SKP6RUQklanw95BG9YuISCpT4e+h3HhMLX4REUlZKvw9FJzOp0v2\niohIalLh7yEd4xcRkVSmwt9DKvwiIpLKVPh7SIP7REQklanw91BuXFfuExGR1KXC30M5KvwiIpLC\nVPh7KDe8ZK+7Rx2KiIhIj6nw91BOPHjKmlpU+EVEJPWo8PdQbjwGoAF+IiKSklT4e6itxa9T+kRE\nJBWp8PeQCr+IiKQyFf4eyokFT5ku2ysiIqlIhb+H1OIXEZFUpsLfQ7nxtha/Cr+IiKQeFf4e2t/i\n16h+ERFJQSr8PaSufhERSWUq/D2krn4REUllKvw9lBMLL+ATFn5duldERFKJCn8P5WYf2NV/5T0v\nMeVbT0UZkoiISNLiUQeQatrO429sCc7jf2HZtijDERER6RG1+HsocXBfS6u6+UVEJLWo8PdQYuFf\ns33P/vk61i8iIqkgksJvZv9mZgvMbL6Z/dbM8syswsyeMrOl4d/yhPVvMbNlZrbYzM6PIuY2OQmj\n+rfUNeyfv7dJl/AVEZHBb8ALv5mNAG4Aprj7MUAMuAK4GXjG3ScAz4T3MbOJ4fJJwAXAXWYWG+i4\n2ySezlff0Lx/fl3CtIiIyGAVVVd/HMg3szhQAKwHLgXuC5ffB1wWTl8KzHD3BndfASwDpg5wvPvt\nH9zX3Mqexnda+fUNavGLiMjgZ1EcmzazG4Fbgb3AX9z9n8xsp7uXhcsN2OHuZWb2Q2C2u98fLrsH\neNzdH+pgu9cB1wFUVVVNnjFjRp/FXFdXR1FREQDXPFnP+aOzGVZo3DO/EYCvTctjdGlkHRGHJDG3\ndJKueUH65paueYFyS0WpltdZZ501192ndLfegJ/OFx67vxQYA+wEfm9mVyau4+5uZj3+RuLudwN3\nA0yZMsWnT59+6AGHZs6cSdv28p59kqrhNYysyIf5CwE48pgTmDauss/2N5ASc0sn6ZoXpG9u6ZoX\nKLdUlK55RdHVfy6wwt23uHsT8EfgVGCTmVUDhH83h+uvA0YmPL4mnBeZ3HgWjS0t1B/Q1a9j/CIi\nMvhFUfhXA6eYWUHYpX8OsAh4BLgqXOcq4M/h9CPAFWaWa2ZjgAnAnAGO+QB52TEamjS4T0REUs+A\nd/W7+0tm9hAwD2gGXiXoni8CHjSzTwCrgMvD9ReY2YPAwnD969090pF0ufEs9rUb3Jc4LSIiMlhF\ncsled/8q8NV2sxsIWv8drX8rwWDAQSEnnkVDUwv1Dc3kxrNoaG6lqUW/1iciIoOfrtzXC7nZseA8\n/sZmygtygHd+tEdERGQwU+HvhaCV30J9QwtlBdkANKrFLyIiKUC/ztcLufEs6hqaaWn1dwq/Wvwi\nIpIC1OLvhbzsGPuaWqlvaKEoN5vsmKnFLyIiKUEt/l5o6+pvaXUKc2PkxLLU4hcRkZSgwt8LufHg\nPP6G5lYKcuJkx7M0ql9ERFKCuvp7ITc7OIVvT2MzRWrxi4hIClHh74XceBZ7G5vZ09hCQU6cnLgK\nv4iIpAYV/l7Ijcf2X6d//zF+dfWLiEgKUOHvhdz4O09bYa5a/CIikjpU+HshLzu2f7qwratfLX4R\nEUkBKvy9kNjiL8iJka3BfSIikiJU+HshN7tdV39Mp/OJiEhqUOHvhWElefundYxfRERSiQp/L4yq\nLNw/PaaykOxYcF6/iIjIYKfC3wsjK/L3T5cWZJOrK/eJiEiKUOHvhdx47ID7GtUvIiKpQtfq76XH\nbjiD4rzg6cuOmY7xi4hISlDh76WJw0v2T+fEs2hq8QijERERSY66+vtAYW6c2n1NuKv4i4jI4KbC\n3weGFuXS1OLs3tccdSgiIiJdUuHvA5VFOQBsq2uIOBIREZGuqfD3gcrCXAC21TdGHImIiEjXVPj7\nwJCisPCrxS8iIoOcCn8fGBJ29W+tU4tfREQGNxX+PlBRmENuPItlm+uiDkVERKRLKvx9IB7L4vTx\nQ3jmrU06pU9ERAY1Ff4+csaEIazZvpeNu/dFHYqIiEinVPj7yMThpQC8taE24khEREQ6p8LfR46q\nLgZg4YbdEUciIiLSORX+PlKSl011aR7Lt9RHHYqIiEinVPj70IiyfNbt3BN1GCIiIp1S4e9DI8rz\nWbdzb9RhiIiIdEqFvw+NKMtnw859tLTqlD4RERmcVPj70PCyfJpbnS21unSviIgMTj0q/GaWZWYl\n/RVMqisryAZg976miCMRERHpWLeF38x+Y2YlZlYIzAcWmtl/9H9oqacoNw5ArQq/iIgMUsm0+Ce6\n+27gMuBxYAzw0X6NKkUV5wUt/tp9zRFHIiIi0rFkCn+2mWUTFP5H3L0J0Oi1DhTnBS3+ugYVfhER\nGZySKfw/BVYChcAsMxsF6PJ0HWjr6q9Ti19ERAapeHcruPsPgB8kzFplZmf1X0ipq0gtfhERGeSS\nGdxXZWb3mNnj4f2JwFX9HlkKKswJCv9utfhFRGSQSqar/17gSWB4eH8JcFN/BZTKYllGUW5cXf0i\nIjJoJVP4h7j7g0ArgLs3Ay39GlUKK8qNU9eg0/lERGRwSqbw15tZJeFIfjM7BdjVr1GlsJL8ODv3\nqPCLiMjg1O3gPuCzwCPAODN7ARgKfLBfo0phw0rz2bBrX9RhiIiIdCiZUf3zzOzdwJGAAYvDc/ml\nAyPK8pm/Th0iIiIyOHVa+M3sA50sOsLMcPc/9lNMKa2mPJ/t9Y3saWymICeZDhUREZGB01Vlel8X\nyxxQ4e9ATXk+AOt27GVCVXHE0YiIiByo08Lv7lcPZCDpYsyQQgCWba5T4RcRkUEnmQv43GZmZQn3\ny83sW/0bVuo6oqqYLINFG3RVYxERGXySOZ3vQnff2XbH3XcAFx3KTs2szMweMrO3zGyRmU0zswoz\ne8rMloZ/yxPWv8XMlpnZYjM7/1D23d/ysmOMGVLIwg21UYciIiJykGQKf8zMctvumFk+kNvF+sm4\nE3jC3Y8CjgcWATcDz7j7BOCZ8H7bJYKvACYBFwB3mVnsEPffr04YWc7cVdtpbdWPGIqIyOCSTOF/\nAHjGzD5hZp8AngLu6+0OzawUOBO4B8DdG8MehUsTtnsfwc8AE86f4e4N7r4CWAZM7e3+B8Jp4yvZ\nsaeJheruFxGRQcbcu2+VmtmFwDnh3afc/cle79DsBOBuYCFBa38ucCOwzt3LwnUM2OHuZWb2Q2C2\nu98fLrsHeNzdH+pg29cB1wFUVVVNnjFjRm/DPEhdXR1FRUVJrbu70fnszD2cPjzOx4851M6R/teT\n3FJJuuYF6ZtbuuYFyi0VpVpeZ5111lx3n9LdekmdaO7ujwOPH3JU7+zzXcC/uvtLZnYnYbd+wv7c\nzHrcT+7udxN8qWDKlCk+ffr0Pgg3MHPmTHqyvTl73mTGnDV88UOTOWKQj+7vaW6pIl3zgvTNLV3z\nAuWWitI1r067+s3s+fBvrZntTrjVmtmh9GGvBda6+0vh/YcIvghsMrPqcJ/VwOZw+TpgZMLja8J5\ng9pnzzuSorw4X3tkQdShiIiI7Ndp4Xf308O/xe5eknArdveS3u7Q3TcCa8zsyHDWOQTd/o8AV4Xz\nrgL+HE4/AlxhZrlmNgaYAMzp7f4HSkVhDv908uG8uHwbDc36MUMRERkckjmP/9fJzOuhfwUeMLM3\ngBOA24DvAOeZ2VLg3PA+7r4AeJDgy8ETwPXunhKVdPxhRbjDmu17ow5FREQESO4Y/6TEO2YWByYf\nyk7d/TWgowEI53QwD3e/Fbj1UPYZhVGVwVX8Vm2rZ/xhqTNARERE0ldXx/hvMbNa4LjE4/vAJt7p\nhpcujA4L/8pteyKOREREJNDVMf5vA6XAr9od369091sGLsTUVV6QTX52jPU71dUvIiKDQ5fH+N29\nFThpgGJJO2ZGdWkeG3fvizoUERERILkr980zMxX/XqoqyWPjLhV+EREZHJIp/CcDL5rZ22b2hpm9\nGY7GlyRUl6rwi4jI4JHMqP5B/Wt4g11VaR6bdu+jtdXJyrKowxERkQzXbYvf3Ve5+ypgL+AJN0lC\ndWkeza3O1vqGqEMRERFJ6gI+l4QX1VkBPAespO+u25/2qkryANi0S4VfRESil8wx/m8CpwBL3H0M\nwUV2ZvdrVGmkujQo/Bt26ZQ+ERGJXjKFv8ndtwFZZpbl7s/S8VX3pAPD2lr8OqVPREQGgWQG9+00\nsyJgFsH19TcD9f0bVvqoLMolnmVs0Mh+EREZBJJp8V8K7AH+jeBHct4G3tefQaWTWJZxWHGuTukT\nEZFBocsWv5ldBowH3nT3J4H7BiSqNDNMV+8TEZFBoqsf6bmLoJVfCXzTzL48YFGlmerSfLX4RURk\nUOiqq/9M4OzwB3mmA5cNSERpqKokaPG76/IHIiISra4Kf6O7twC4+x5Al53rperSPPY0trB7X3PU\noYiISIbr6hj/UQnX5DdgXHjfAHf34/o9ujRxWEkuAFtq91Ganx1xNCIiksm6KvxHD1gUaa6iMAeA\nbXWNjD8s4mBERCSjdVr4w+vzSx+oLAxa/NvrGyOOREREMl0y5/HLIaosClv8KvwiIhIxFf4BUF7w\nTle/iIhIlLo6j/+Z8O93By6c9JQTz6I4L852/TSviIhErKvBfdVmdipwiZnNoN3pfO4+r18jSzND\ninLZqq5+ERGJWFeF/yvAl4Ea4PZ2yxw4u7+CSkcjyvJZu31P1GGIiEiG62pU/0PAQ2b2ZXf/5gDG\nlJZGVRbw6BsbAJi3egeThpeQG49FHJWIiGSabn+W192/aWaXEFzCF2Cmuz/av2Gln9GVheza28SS\nTbV84K6/M3VMBQ9+clrUYYmISIbpdlS/mX0buBFYGN5uNLPb+juwdDOqsgCAWUu2ADBnxXZaWnXt\nfhERGVjdtviB9wInuHsrgJndB7wKfLE/A0s3w8vyAXhl5Y798zbs2ktNeUFUIYmISAZK9jz+soTp\n0v4IJN0NK80D4JVV7xT+lVs12E9ERAZWMi3+bwOvmtmzBKf0nQnc3K9RpaGKghxyYllsrXvnXP6V\n2+o5fcKQCKMSEZFM022L391/C5wC/BH4AzDN3X/X34Glm6ws2/8rfaeMraA4L86C9bsijkpERDJN\nMi1+3H0D8Eg/x5L2ygqyWbtjL5NHlVOUm80Ly7ZFHZKIiGQYXat/AH3uPUdSlBvn7KOqOGl0Oau3\n72HXnqaowxIRkQySVItf+sZZRx7Gm197D2bG+p17AdhUu4/SguyIIxMRkUzRZYvfzGJm9tZABZMJ\nzIKfPKgqCUb5b9q9L8pwREQkw3RZ+N29BVhsZocPUDwZoyoc6Ldpt36xT0REBk4yXf3lwAIzmwPU\nt81090upRkP4AAAgAElEQVT6LaoMcFixWvwiIjLwkin8X+73KDJQfk6Mkrw4G3btjToUERHJIMn8\nSM9zZjYKmODuT5tZAaCflesD4w4rYsmmuqjDEBGRDJLMj/RcCzwE/DScNQJ4uD+DyhRHV5ewaMNu\n3PVjPSIiMjCSOY//euA0YDeAuy8FDuvPoDLFxOoSavc1s2yzWv0iIjIwkin8De7e2HbHzOKAmqh9\n4PxJw8iJZ/Hr2auiDkVERDJEMoX/OTP7IpBvZucBvwf+t3/DygxDi3M5ffwQXlq+PepQREQkQyRT\n+G8GtgBvAp8EHgP+sz+DyiTHjihl6eZa9jQ2Rx2KiIhkgGRG9bea2X3ASwRd/Itdo9H6zHE1pbQ6\nLFy/mymjK6IOR0RE0lwyo/rfC7wN/AD4IbDMzC7s78AyxbEjSgF4Y61+oldERPpfMhfw+T5wlrsv\nAzCzccD/AY/3Z2CZ4rCSPIaV5PHmOhV+ERHpf8kc469tK/qh5UBtP8WTkY6tKeWNtTujDkNERDJA\npy1+M/tAOPmKmT0GPEhwjP9DwMsDEFvGOG5EKU8v2kTtviaK8/QTvSIi0n+66up/X8L0JuDd4fQW\nIL/fIspAx9aU4g7z1+1m2rjKqMMREZE01mnhd/er+3PHZhYDXgHWufvFZlYB/A4YDawELnf3HeG6\ntwCfAFqAG9z9yf6MbaCdOLIcM3hpxTbeNaqM3Lh+CkFERPpHMqP6x5jZ7Wb2RzN7pO3WB/u+EViU\ncP9m4Bl3nwA8E97HzCYCVwCTgAuAu8IvDWmjtCCbY4aXcsfTSznyP5/gn34+m+VbdBlfERHpe8kM\n7nuYoAX+PwQj/NtuvWZmNcB7gZ8nzL4UuC+cvg+4LGH+DHdvcPcVwDJg6qHsfzC69syxnDK2go9N\nG8X8dbv59P3zaG5pjTosERFJM8mczrfP3X/Qx/u9A/g8UJwwr8rdN4TTG4GqcHoEMDthvbXhvLRy\nyfHDueT44QBMHVPBZ37zKs8t2cI5R1d180gREZHkJVP47zSzrwJ/ARraZrr7vN7s0MwuBja7+1wz\nm97ROu7uZtbjqwOa2XXAdQBVVVXMnDmzNyF2qK6urk+315W8Vqc4G+5+8lVim/L6fX8DmdtASte8\nIH1zS9e8QLmlonTNK5nCfyzwUeBsoK3v2cP7vXEacImZXQTkASVmdj+wycyq3X2DmVUDm8P11wEj\nEx5fE847iLvfDdwNMGXKFJ8+fXovQzzYzJkz6cvtdee8La8xc8kWzjzz3WRlWb/ua6BzGyjpmhek\nb27pmhcot1SUrnklc4z/Q8BYd3+3u58V3npb9HH3W9y9xt1HEwza+6u7Xwk8AlwVrnYV8Odw+hHg\nCjPLNbMxwARgTm/3nypOHT+E7fWNLNmsayWJiEjfSabwzwfK+jsQ4DvAeWa2FDg3vI+7LyC4eNBC\n4AngendvGYB4InV0dTD8YfmW+ogjERGRdJJMV38Z8JaZvcyBx/gvOdSdu/tMYGY4vQ04p5P1bgVu\nPdT9pZKasgIA1u3YG3EkIiKSTpIp/F/t9yjkICX5cYpy46zbqcIvIiJ9p9vC7+7PDUQgciAzY0RZ\nPmvV4hcRkT7UbeE3s1qCUfwAOUA2UO/uJf0ZmMCYIYW8uW4X7o5Z/47sFxGRzNDt4D53L3b3krDQ\n5wP/ANzV75EJp08Ywrqde1m+VQP8RESkbyQzqn8/DzwMnN9P8UiCdx8xFIDnFm+JOBIREUkXyXT1\nfyDhbhYwBdjXbxHJfiMrChg7tJDnlmzhqlNHc9PvXmPq6HI+Om101KGJiEiKSmZU//sSppsJfrDn\n0n6JRg5y6rhK/jRvHZ978DX+9/X1/O/r6znrqMOoKS+IOjQREUlByRzjvzrhdq273+rum7t7nPSN\nidWl1De28PBr63nPxCriWcYvX1gZdVgiIpKiOm3xm9lXunicu/s3+yEeaaftCn4A//OPJ3L9A/N4\nauEmvnzxxAijEhGRVNVVi7++gxvAJ4Av9HNcEjq6uoTJo8r5xcenkBuPcdr4IazevodnF6vTRURE\neq7Twu/u32+7EfziXT5wNTADGDtA8WW8vOwYf/j0qZx9VBUAFx1bzYiyfD7167ks3aQf8BERkZ7p\n8hi/mVWY2beANwgOC7zL3b+gY/zRqSrJ40/Xn0pjSyuPz98YdTgiIpJiOi38ZvZfwMtALXCsu3/N\n3XcMWGTSqcOK8zh6WAkvrdgWdSgiIpJiumrxfw4YDvwnsN7Mdoe3WjPbPTDhSWemjqlg7qodNDa3\nRh2KiIikkK6O8We5e37iJXvDW7Gu0x+9k8dUsK+plTfX7Yo6FBERSSE9umSvDB5Tx1SQZfCXBRtx\ndzbu2kdzi1r/IiLStWSu3CeDUGVRLhceW81PZy3nsfkbWLN9L2OHFHLx8cO5+tTRlBfmRB2iiIgM\nQir8Kewbl0yisjCHzbsbOGPCUJ6Yv5EfPLOUvY3NfOm9usCPiIgcTIU/hVUW5fKNS4/Zf/+29x/L\nv/72VWa8vIbPvedI8rJjEUYnIiKDkY7xp5kPTxlJ7b5mnl60KepQRERkEFLhTzPTxlVSVZLLn+at\nizoUEREZhFT400wsy7jsxBE8t2QLr63ZyQ2/fZXX1+yMOiwRERkkVPjT0AdOrKG51bnsRy/wyOvr\nef9dL/DE/A1RhyUiIoOACn8aOnJYMaeNrwRg7NBCWh0+++Dr7NrTFHFkIiISNRX+NPXfHzqeL150\nFH/6l9N4/MYz2NPYwveefIvWVo86NBERiZAKf5qqLs3nujPHUZqfzdHVJVx5yuE88NJqfvvy6qhD\nExGRCKnwZ4hvXHIMI8ryeXqhTvMTEclkKvwZIivLOPuow3hpxXa21TVEHY6IiEREhT+DfGzaKBqa\nW7nn+RVRhyIiIhFR4c8gE6qKmTKqnFlLtwDg7vzmpdU8tKSRNdv3RBydiIgMBF2rP8OcNn4Itz+1\nhNseW8TOPY08+MpaAB793rNcf9Y4/uP8oyKOUERE+pNa/BnmylNGcdSwYu6etZwHX1nLlFHlfPv0\nfKaOqeDHM99m7Y49PLNoE5f+8Hnmrd4RdbgiItLH1OLPMBWFOTx+4xnsa2pl0cbdTBpewovP/43v\nf+h4zvjes1z3q7ks2VRLc6tz+1+WcP81J0cdsoiI9CG1+DOQmZGfE+Ndh5eTGw9+undkRQEXHjOM\nhRt2c0RVMf8yfRzPL9vK/HW7Io5WRET6klr8st/tl5/A+47fzPQjh9LU4vx69iq+/fgiHrjmlKhD\nExGRPqIWv+yXnxPjomOrKciJU5qfzWfOGs8Ly7axYmt91KGJiEgfUeGXTl10bDUAf5q3NuJIRESk\nr6jwS6fajvv/4K/L+OFfl0YdjoiI9AEVfunS//vwCbz7iKHcNfNt/ayviEgaUOGXLuVlx/jseUew\np7GFJxdsjDocERE5RCr80q3jakoZWZHP/725IepQRETkEKnwS7fMjIuOqeaFZVvZuacx6nBEROQQ\nqPBLUt53/HCaW527Zy3n78u28qsXV7KnsTnqsEREpId0AR9JyjEjSpk6poK7Zr7NXTPfBuAP89bx\no388kZrygoijExGRZKnwS9J+cuVkXl+7kwXrdrF8Sz1/fHUdZ3//OR761DSOqymLOjwREUmCuvol\naRWFOZx15GF85uwJ3P7hE3jshjMoycvme08sjjo0ERFJkgq/9NrE4SVcPqWGF5dv0zn+IiIpQoVf\nDsn5k4bR0ur8evbKqEMREZEkqPDLITl+ZBnvPbaa//7LEkbf/H9c/cs5tLZ61GGJiEgnVPjlkH33\ng8dx7IhSAJ5dvIVV2/dEHJGIiHRGhV8OWVFunN9/aho/uXIyAIs27I44IhER6cyAF34zG2lmz5rZ\nQjNbYGY3hvMrzOwpM1sa/i1PeMwtZrbMzBab2fkDHbN0Ly87xvQjh5JlsGD9rqjDERGRTkTR4m8G\nPufuE4FTgOvNbCJwM/CMu08AngnvEy67ApgEXADcZWaxCOKWbuRlx5g8qpwn5m/EXcf5RUQGowEv\n/O6+wd3nhdO1wCJgBHApcF+42n3AZeH0pcAMd29w9xXAMmDqwEYtybrsxBG8vaWepZvrog5FREQ6\nEOkxfjMbDZwIvARUuXvbz79tBKrC6RHAmoSHrQ3nySB05oShAMyYs4aXlm+juaU14ohERCSRRdUl\na2ZFwHPAre7+RzPb6e5lCct3uHu5mf0QmO3u94fz7wEed/eHOtjmdcB1AFVVVZNnzJjRZ/HW1dVR\nVFTUZ9sbTPo6t1v+tocN9cH76rghMT47Ja/Ptt0Tes1ST7rmBcotFaVaXmedddZcd5/S3XqRXKvf\nzLKBPwAPuPsfw9mbzKza3TeYWTWwOZy/DhiZ8PCacN5B3P1u4G6AKVOm+PTp0/ss5pkzZ9KX2xtM\n+jq3P0/ex+pte/jDvLX8ds4aJpxwMiPK8vts+8nSa5Z60jUvUG6pKF3zimJUvwH3AIvc/faERY8A\nV4XTVwF/Tph/hZnlmtkYYAIwZ6DilZ47rDiPKaMruO7McQA8MX9jxBGJiEibKI7xnwZ8FDjbzF4L\nbxcB3wHOM7OlwLnhfdx9AfAgsBB4Arje3VsiiFt6aMyQQo4aVsyTKvwiIoPGgHf1u/vzgHWy+JxO\nHnMrcGu/BSX95oJjhnHnM0vZvHsfr67ZSUleNtPGVUYdlohIxorkGL9kjguOGcYdTy9l6m3P7J93\n/VnjuHzKSEZVFkYYmYhIZtIle6VfHVlVzAWThgFw7RljeO+x1fzo2bf5hx//nX1NOmIjIjLQ1OKX\nfmVm/PAfT2RTbQMjyvJxd06bM4Qv/ulNPnX/XO656iRiWZ0d+RERkb6mFr/0u3gsa//pfGbGR6aO\n5N/fcwQzF2/hjqeX8MbanRFHKCKSOVT4ZcCZGdefNZ6Ljh3G//x1GZf88AV+9eLKqMMSEckIKvwS\nCTPjzitO5I4Pn8Cwkjxue2wRSzfVRh2WiEjaU+GXyGTHsrjsxBH86fpTKc7L5mO/mMPaHXuiDktE\nJK2p8Evkqkvzue/qqdTta+b07z7Lo2+sjzokEZG0pcIvg8LE4SU8+Klp5GfH+J9nlhHVj0eJiKQ7\nFX4ZNI6uLuFL7z2axZtqWbRBx/tFRPqDCr8MKhccM4wsg8fe3BB1KCIiaUmFXwaVIUW5TBtXyf+9\nuYHWVmfWki1c9qMXuOCOWfxs1nIdAhAROUQq/DLoXD5lJCu21nPF3bO5/oF5bKltoCAnxq2PLWLi\nV57k248tor6hmcbm1qhDFRFJObpkrww6lxw/nM27G/jNnNWMGVrInVecyKiKAu57cSVPL9rET2ct\n56ezlpMdMz56ymhuOm8CJXnZUYctIpISVPhl0DEzrj1zLNeeOfaA+VefNoaPnzqapxdtZtnmOp59\nazO/eGEFM5ds5rv/cBwnja6IKGIRkdShrn5JKWbGeROr+PT0cfzuk6fwpYuOZvW2PVxx92xWbK2P\nOjwRkUFPhV9SVlvPwPNfOJuYGffPXhV1SCIig54Kv6S8YaV5vGtUGa+s2hF1KCIig54Kv6SF42rK\nWLRhN00tGukvItIVFX5JC+86vJzG5lZ+9/IaADbu2sdLG5r1RUBEpB2N6pe0cN7EKiaPKuc/H57P\n5toGnpy/kcWbGljaOI+fXPku4jF9xxURARV+SROxLOOBa07mxhmv8oNnlgJQlA1PL9rExK88iVnw\nM8Anj6lg/GFF1Dc287elW/nUu8fxkamHRxy9iMjAUeGXtJGXHePOK07k/tmrOLq6hL2r32RP5ZEs\nWL8LgC27G3j4tXU889Zmsgwqi3K55Y9vMnZIISePrYw4ehGRgaHCL2klLzvGNWcEF/6Zuda45Pjh\nXHL88P3Lb/vAsTQ0tdLiTkFOjNO/+1f+56/LmDqmAjOLKmwRkQGjA5+SUfKyY5QWZFNRmENedoxP\nTx/P88u28uzizQDs3tekHwISkbSmwi8Z7WPTRjGiLJ87n17Ktx9fxHFf+wv/fO/LOhtARNKWCr9k\ntOxYFp+/4EheX7uLnz63nDFDCnl28RbufWFl1KGJiPQLHeOXjHfpCSMozImzeFMtn3r3OK771Svc\n+tgiXlqxnXiWcfHx1azatofRlYW897jqqMMVETkkKvwiwLkTqzh3YhUAX790EnUNzfz1rU0U5sZ5\nYsFGALIMttVPIi8e4+SxFYyqLIwyZBGRXlHhF2mnpryA331yGs0trTQ0t/LBn7zI1NHlLNywm6/8\neQEAQ4tzeeELZ5MT19EyEUktKvwinYjHsojHsnjshtMxM5pbWnlx+TbeXLeL7z2xmGcWbeLCY9X1\nLyKpRc0VkW60nd8fj2VxxoShXHfGWEZVFvBff1lMS6sza8kW/vnel3li/saIIxUR6Z5a/CI9FI9l\n8dnzjuDGGa/xyV/P5elFmwB48e1tjD/sNFZu3UNreC2Agpw408ZVEsvSxYFEZHBQ4RfphfMnDaM4\nL87TizYxsiKfL793Ip9+YB7n3j7roHUvn1LDTeceAUBhTpzSguyBDldEZD8VfpFeyMuO8chnTmfO\nim18cPJIYlnGg5+cxotvb+Vdh5dTkh8U90ff2MBPnnubB19ZCwRnBtxwzgRuPGcCa3fsZURZPlnq\nDRCRAaTCL9JLY4YUMmbIO6f0TR5VzuRR5QesM2l4CSeNLmdrXQMAs5Zu5Y6nl/KzWcupb2zhQ5Nr\n+K8PHT+gcYtIZlPhF+lHZsY5R1ftv/+hySM5Y/wQFqzfzfPLtvL7uWspK8jm6OoSpoyqICscbptl\nxrCSPPUGiEifU+EXGUBZWcYVUw8HoLG5lU/c9zI/+9uKDtcdVpJHaX42Jx5exk3nHsGw0ryBDFVE\n0pQKv0hEcuJZ/Oqfp7JrbxOLN9ayavue/ct2721i3uodNDa38sd563h28WYevv40qkvzI4xYRNKB\nCr9IhMyMsoIcTh5bycljKztcZ9GG3bz/rhf4t9+9xgdOrOHVNTs5s0Q/HSwivaPCLzLIHV1dwtcv\nmcSX/7yA2cu3A7CuJs6F50YcmIikJBV+kRTw4ZMO54wJQ3l55XaefWszD7+2ntseW8RN504AICeW\nRXOrkxvP2n+lQRGRjqjwi6SI4WX5XHrCCE4dN4T5Kzdy96zl3D1r+QHrjB1aSFVxHtnxLL72vomM\nHVoEBAMJzSA7pqt0i2Q6FX6RFDO0OJcvTM1ja9E4du9twoElm2opzo2zdHMdzS3OvOXb+OyDr3PS\n6HIamlt5+NV15MRjXHxcNdmxznsEzIwPTq7hiKrigUtIRAaUCr9ICsrOMj4SnhbYkR88s5SfPPc2\nSzbVAnDK2Er2NDbz4CtrutzuvqYWXlqxnStPPpzZy7czb/UOttQ24OFvD0wdU8Htl59AWUE29Y0t\n+x8XzzLysmN9kJmI9DcVfpE0dMM5E7jhnAk9ftwvX1jB1/93Ia+v2UlhTozqsnxOGl3O2KFF7Gtq\n4fevrOWCO2dRmBNn+db6/Y+LZxnH1ZQSDw8lVBbmkJ8TY+2OvZ3ua2J1CR8+aSTZsSzmrdrBw6+t\no6CpgSVZbwPBRYwuPm448Zjx6OvriWUZ739XDUW5+tgSORT6DxKR/T5+6mjOOaqK5tZWqkvzyc85\nsBV/8XHDuXvW2zS1OJecMJzCnOAjZMW2epZvqQPAHf7+9jb2NbVwfE3Z/qsRJmpthftnr+Lev6/c\nP6+mPJ91O5p5evVb++d994m3iGdlsbcp6F349uNvUVWSx4kjyygryDlou+UF2ZxxxFCGl+axbHPd\nAcsKc+McV1PKjj1NPDB7FTv2NAFwzIgShpXkEY9lceLhZRoHIWlPhV9E9jMzDq8s6HT5tHGVTBvX\n8fUGEjW3tNLqwUWKOrNqWz2vr90FBD0Ep46r5Km/zuS0088A4K2Nu3ngpdUU5MT46CmjWbG1nsfn\nb2D33iaeXrQJb3cpAwfqGpr5/lNLOt1ndsxodWh1pygnTkNLK43NrfuXx7KsyzEQbXLjMU4dV8mp\n44dQmNP1IY6Jw0uoLs1nb7OuvSCDgwq/iPS5eBKt5lGVhYyqLDxgXk7MKAy78iePqmDyqIr9y44c\nVswFxwzrcptrtu/hb0u3snH3PqaOrjigiC/fWs/KrfXEY8b7TxzB+MOKaW5p5fW1u2huaWVTbQML\n1u1KKr+1O/fy4tvbeHz+xqTWb3P0gr8xsjyfM44YypRR5cTC32LIMhhdWUh9Qwu79wU9ESV52foJ\nZ+kXKvwikjZGVhTwjyd3POixoysjxmNZB/yi4iXHD096X00trazf2fkYBoB9Ta3MWbGNxhZn7oKl\n7Ipl8/LK7fxl4aaD1i0ryGZPQwuNLUEPRH52jPcdX01BTpz3nziC40eWJR2bSFdU+EVEeiE7lnVQ\nj0VHjhwWnBo5rnkV06efQkurs3D9blYn/DbDrr1NzF21g6LcGMNK83GCdf6ycBM79zRx799XktOu\nF6W8MJsTRpaRFV6waURZPkcMK2bRht3U7ms+oEehTVFunCmjK4j38a8+1jU6O+ob+3Sbg0FXeeVm\nZ1GQ030JbWhuYU9Dy0Hz83NikZ0JkzKF38wuAO4EYsDP3f07EYckItJjsSzj2JpSjq0pPWB+Zz0V\nm3fv4zdzVtOQMBbBPfgNhxXhmRUtrc6TCzbS6sH2c2JZPDR3bf8l0ZG/PjWw+xsoneQVyzImDS85\n6MtVIndYuqn2gFNf2+Rnx/j8BUdy9Wlj+izUZKVE4TezGPAj4DxgLfCymT3i7gujjUxEpH8dVpLH\nTece0e16O+obqWtopjgvTm48xta6hoPWeXtLHSsTTsPsK0uXLWPC+PF9vt2odZXXiq31B5zS2pn3\nTBrGcTWlJH49cODlldspzotmDEdKFH5gKrDM3ZcDmNkM4FJAhV9EBCgvzKG88J1THEdWHHx2xsiK\nAjiy7/c9s2kV0yNoufa3/swripZ+G/P258QMQmb2QeACd78mvP9R4GR3/0y79a4DrgOoqqqaPGPG\njD6Loa6ujqKioj7b3mCSrrmla16Qvrmla16g3FJRquV11llnzXX3Kd2tlyot/qS4+93A3QBTpkzx\n6dOn99m2Z86cSV9ubzBJ19zSNS9I39zSNS9QbqkoXfNKlUtUrQNGJtyvCeeJiIhID6RK4X8ZmGBm\nY8wsB7gCeCTimERERFJOSnT1u3uzmX0GeJLgdL5fuPuCiMMSERFJOSlR+AHc/THgsajjEBERSWWp\n0tUvIiIifUCFX0REJIOo8IuIiGQQFX4REZEMosIvIiKSQVT4RUREMogKv4iISAZJiR/p6Q0z2wKs\n6sNNDgG29uH2BpN0zS1d84L0zS1d8wLllopSLa9R7j60u5XStvD3NTN7JZlfPUpF6ZpbuuYF6Ztb\nuuYFyi0VpWte6uoXERHJICr8IiIiGUSFP3l3Rx1AP0rX3NI1L0jf3NI1L1BuqSgt89IxfhERkQyi\nFr+IiEgGUeHvhpldYGaLzWyZmd0cdTw9ZWa/MLPNZjY/YV6FmT1lZkvDv+UJy24Jc11sZudHE3X3\nzGykmT1rZgvNbIGZ3RjOT4fc8sxsjpm9Hub29XB+yucGYGYxM3vVzB4N76dLXivN7E0ze83MXgnn\npUtuZWb2kJm9ZWaLzGxaOuRmZkeGr1fbbbeZ3ZQOuXXJ3XXr5AbEgLeBsUAO8DowMeq4epjDmcC7\ngPkJ874H3BxO3wx8N5yeGOaYC4wJc49FnUMneVUD7wqni4ElYfzpkJsBReF0NvAScEo65BbG+1ng\nN8Cj6fJ+DONdCQxpNy9dcrsPuCaczgHK0iW3hBxjwEZgVLrl1v6mFn/XpgLL3H25uzcCM4BLI46p\nR9x9FrC93exLCf6RCf9eljB/hrs3uPsKYBnBczDouPsGd58XTtcCi4ARpEdu7u514d3s8OakQW5m\nVgO8F/h5wuyUz6sLKZ+bmZUSNCDuAXD3RnffSRrk1s45wNvuvor0y+0AKvxdGwGsSbi/NpyX6qrc\nfUM4vRGoCqdTMl8zGw2cSNAyTovcwu7w14DNwFPuni653QF8HmhNmJcOeUHw5expM5trZteF89Ih\ntzHAFuCX4SGan5tZIemRW6IrgN+G0+mW2wFU+DOcB/1XKXtqh5kVAX8AbnL33YnLUjk3d29x9xOA\nGmCqmR3TbnnK5WZmFwOb3X1uZ+ukYl4JTg9fswuB683szMSFKZxbnOBw4Y/d/USgnqD7e78Uzg0A\nM8sBLgF+335ZqufWERX+rq0DRibcrwnnpbpNZlYNEP7dHM5PqXzNLJug6D/g7n8MZ6dFbm3CLtVn\ngQtI/dxOAy4xs5UEh83ONrP7Sf28AHD3deHfzcCfCLqA0yG3tcDasNcJ4CGCLwLpkFubC4F57r4p\nvJ9OuR1Ehb9rLwMTzGxM+I3wCuCRiGPqC48AV4XTVwF/Tph/hZnlmtkYYAIwJ4L4umVmRnDMcZG7\n356wKB1yG2pmZeF0PnAe8BYpnpu73+LuNe4+muB/6a/ufiUpnheAmRWaWXHbNPAeYD5pkJu7bwTW\nmNmR4axzgIWkQW4JPsI73fyQXrkdLOrRhYP9BlxEMGL8beBLUcfTi/h/C2wAmgi+uX8CqASeAZYC\nTwMVCet/Kcx1MXBh1PF3kdfpBN1vbwCvhbeL0iS344BXw9zmA18J56d8bgnxTuedUf0pnxfBmT+v\nh7cFbZ8V6ZBbGOsJwCvhe/JhoDyNcisEtgGlCfPSIrfObrpyn4iISAZRV7+IiEgGUeEXERHJICr8\nIiIiGUSFX0REJIOo8IuIiGSQeNQBiMjg9f/bu39dCMIoDOPPqxadRC+R0FAgQkPiDpSi13ANSoVE\nXIhCu52WjkTHVWyEgqPYEXSy2cwW3/NLJvMnU5xp5s18M9+cJN/TmgAWgA9Gv28FeK2q7akUJmls\nTueT9C9JzoBhVV1MuxZJ43OoX9JYkgy79W6S2yQ3SV6SnCc5THLX9adf7M6bT3Kd5L5bdqZ7BVKb\nDCxj9GwAAAC3SURBVH5Jk7AKHAPLwBGwVFWbjNrvnnTnXAGXVbUBHPC3Na+knviOX9Ik3FfXxjTJ\nMzDojj8Ce932PrAyarMAwFyS2aoa9lqp1DiDX9IkvP/a/vy1/8nPfWYG2Kqqtz4Lk/SXQ/2S+jLg\nZ9ifJGtTrEVqlsEvqS+nwHqShyRPjL4JkNQzp/NJktQQn/glSWqIwS9JUkMMfkmSGmLwS5LUEINf\nkqSGGPySJDXE4JckqSEGvyRJDfkCnkB6JSwWcPAAAAAASUVORK5CYII=\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 }