{ "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/bathymetry_201702.nc')" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARkAAAKvCAYAAABea7viAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2wXPV93/HPJygmCXXARjJgRItMZSikPCiqGoaEjg02\nD3UjEjcdZZwUnMzQdgh2iycU6iaTJkNM8JgET2scFQzMlIQwSoiV1MLGOIlbJkaRsawgJCwZlOoy\nCAQ2cUISEfC3f9xzzWq1D+ecPb+z5+H9mtHo7t7dvUd32De/3+88rCNCAJDKd817AwB0G5EBkBSR\nAZAUkQGQFJEBkBSRAZBUrsjYPs72Jtu7be+yfb7t37G9Pfuzz/b2gcffaHuv7SdtX5Ju8wGkkr2v\n/zx7j2/L7nuz7Yds78n+ftPU18lznIzteyT9n4i4w/YbJH1fRLw08P2PSfrLiPhl22dK+m1J6yS9\nVdLnJb09Il4r9S8FMBe290laGxEvDNx3i6RvRMTNtm+Q9KaI+M+TXmfqSMb2sZIulHSnJEXEK0OB\nsaR/o8WwSNJ6SfdFxKGIeFrSXi0GB0D7rZd0T/b1PZKumPaEZTledJWkg5Lusn2OpC9L+mBEvJx9\n/0ckPRcRe7LbJ0v60sDzF7L7DmP7aklXS9Ixxxzzg2eccUaOTUFX7dr33BH3vfp94x9/9Pe+8p2v\n/+prz78QEStSbFcK6/7FMfGX35zvwP5rf35op6S/G7hrY0RsHHpYSPq87dck/Wb2/RMi4tns+wck\nnTDtZ+WJzDJJayRdGxGP2r5N0g2SfiH7/k/q9VFMbtkGb5SktWvXxrZt24q+BDpk3VW3HnHfC2vG\nT+VPO3vhO19//h2/8RdJNiqRv/zma7p98z+a6zZctOprfxcRa6c87Icj4hnbb5H0kO3dg9+MiLA9\ndb0lz8LvgqSFiHg0u71Ji9GR7WWSflzS7ww8/hlJpwzcXpndB6BFIuKZ7O/nJT2gxWWP52yfJEnZ\n389Pe52pkYmIA5L22z49u+siSU9kX18saXdELAw8ZbOkDbaPtr1K0mpJW3P9qwA0gu1jbL9x6WtJ\n75b0uBbf31dmD7tS0qenvVae6ZIkXSvp3mzP0lOS3p/dv0FDU6WI2Gn7fi2G6FVJ17BnCWidEyQ9\nsLhfR8sk/VZEPGj7zyTdb/tnJf2FFnf6TJQrMhGxXdIR87eIuGrM42+SdFOe1wbQPBHxlKRzRtz/\nohZnM7lxxC9aZ3DRF81HZAAkRWQAJEVkACRFZAAkRWTQSJOO9kW7EBkASREZAEkRGQBJERkASREZ\nAEkRGbQKpxS0D5FBI2y9+7p5bwISITJoDELTTUQGjUJouofIoHEITbcQGTTSUx/40Lw3ARUhMmgs\nQtMNRAZAUkQGjcZopv2IDFrl6ztWznsTUBCRQesQmnYhMgCSIjJovOWP+Yj7GM20B5FBKxCa9iIy\nAJIiMmgNRjPtRGTQKqNCg2YjMmi84RMmCU27EBm0AmdmtxeRQWsMhobRTHsQGbQKoWkfIoPWITTt\nQmQAJEVk0EosBLcHkUFrEZp2IDJoNULTfEQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBS\nRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJE\nBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUrki\nY/s425ts77a9y/b52f3XZvfttH3LwONvtL3X9pO2L0m18QCab1nOx90m6cGI+Ne23yDp+2y/Q9J6\nSedExCHbb5Ek22dK2iDpLElvlfR522+PiNcSbD+Ahps6krF9rKQLJd0pSRHxSkS8JOk/SLo5Ig5l\n9z+fPWW9pPsi4lBEPC1pr6R1KTYeQPPlmS6tknRQ0l22v2L7DtvHSHq7pB+x/ajtP7H9z7LHnyxp\n/8DzF7L7DmP7atvbbG87ePDgjP8MAE2VJzLLJK2RdHtEnCfpZUk3ZPe/WdIPSfp5Sffbdt4fHBEb\nI2JtRKxdsWJF8S0H0Ap5IrMgaSEiHs1ub9JidBYk/V4s2irp25KWS3pG0ikDz1+Z3Qegh6ZGJiIO\nSNpv+/TsroskPSHp9yW9Q5Jsv13SGyS9IGmzpA22j7a9StJqSVsTbDuAFsi7d+laSfdme5aekvR+\nLU6bPmX7cUmvSLoyIkLSTtv3azFEr0q6hj1LQH/likxEbJe0dsS3fmrM42+SdNMM2wWgIzjiF0BS\nRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJEBkBSRAZAUkQGQFJE\nBkBSRAZAUnkvv9k6BT44AQ2yeAVXdAkjGQBJERkASXV2ujSIIXizMbXtNkYyAJIiMgCSIjIAkiIy\nAJIiMgCS6sXepWGXrb5+6mO27Lmlhi0Buq93I5k8gSnyOACT9S4yRRAaYHa9ikyZaBAaYDa9iQyx\nAOajF5GZNTAECiivF5GpAqEByiEyBRAaoDgiUxChAYohMiUQGiA/IlMSoQHyITIzIDTAdEQGQFJE\nZkaMZoDJiEwFCA26yvZRtr9i+w+z22+2/ZDtPdnfb5r2GkSmIoQGHfVBSbsGbt8g6eGIWC3p4ez2\nRESmQoQGXWJ7paR/KemOgbvXS7on+/oeSVdMe51eXrQqpctWX88FrzDRs4eO1Uf2XT7nrfjactvb\nBu7YGBEbhx70G5Kul/TGgftOiIhns68PSDph2k9iJJMAIxq0wAsRsXbgz2GBsf0eSc9HxJfHvUAs\nfqDZ1A81IzKJEBq03AWSftT2Pkn3SXqn7f8l6TnbJ0lS9vfz016IyCREaNBWEXFjRKyMiFMlbZD0\nhYj4KUmbJV2ZPexKSZ+e9lpEBkARN0t6l+09ki7Obk9EZBJjNIO2i4g/joj3ZF+/GBEXRcTqiLg4\nIr4x7flEpgaEBn1GZGpCaNBXHCdTgxcvOFGStO6qWw+7//hHDhx2m+Nr0EWMZBJbCsy47y39kRjt\noJuITEKTAjPusYQGXUNkEikSmOHnEBp0SS8iU/daR5nADD+X0KArehEZqb7QzBKY4dcgNOiC3kRG\nSh+aKgKT4rWAeepVZKR0oUkVBUYzaDuOk6lAnsC8sCa0/DEXft3jHznQ+2vUENp2691IRqp2NJM3\nMIN/l9HXN1pf/91d0svIVGHwILpJhsNSNDSDP6Nvb7i+/Xu7isiUkHf9ZVxQGNFM15d/Zx8QmYJm\nDUwVuv4G7Pq/r2+ITAFVBmbWCHX1jdjVf1efEZmcUoxg8j6WY2bG43fTfESmQmVGJ6zPlJN34R3z\nx3EyOZTZi1SXaaHp2vE1hKV9GMl0XBdGO8PX3UG79DYybfo//PAV9IpqcmimbRthab/eRkZqRmiK\nnmpQVhND08RtQvV6HRmpmtCkXI+ZdRQzqElv6jzbwiimG3JFxvZxtjfZ3m17l+3zbf+S7Wdsb8/+\nXD7w+Btt77X9pO1L0m1+NeY1oqlrFDOoCaGpMjDzWnBHfnlHMrdJejAizpB0jqRd2f2/HhHnZn8+\nI0m2z9Tix1qeJelSSZ+wfVTF2125JkydhlU5ihk0z9BUFZgX1gSBaYmpkbF9rKQLJd0pSRHxSkS8\nNOEp6yXdFxGHIuJpSXslratiY1GdeYSmysCgPfIcJ7NK0kFJd9k+R9KXJX0w+961tv+tpG2SPhQR\n35R0sqQvDTx/Ibuv8bbsuaW2N9+0qVKqUcygJkydiiAu7ZRnurRM0hpJt0fEeZJelnSDpNslvU3S\nuZKelfSxIj/Y9tW2t9nedvDgwWJbnVBd06ZpbxiOCzkcgWmvPJFZkLQQEY9mtzdJWhMRz0XEaxHx\nbUn/U69PiZ6RdMrA81dm9x0mIjZGxNqIWLtixYry/4IaTBtVlF3AzfPG4UC0yb+n085eqHFLUMbU\nyETEAUn7bZ+e3XWRpCdsnzTwsB+T9Hj29WZJG2wfbXuVpNWStla4zXMxz9AsGQxOX6JDYNov77lL\n10q61/YbJD0l6f2SPm77XEkhaZ+kfydJEbHT9v2SnpD0qqRrIuK1qjcchy+S1rGGk8KLF5y4uGtg\nhHGBIS7tkisyEbFd0tqhu396wuNvknTTDNvVSMc/cmDiCGL5Yy59Jvasx8wsXXS86whM+/T+iN+m\nqGJhs21TqEnbO+r3MSowN576mUq3CdUjMgWlWpuRqjnArC2hKRKY085eIDAtxvVkanLa2Qv6+o6V\nuR476/Rp6Q3c1OlT0RHMMOLSLkRmhEnHyly2+vpSazN5A7Nk6fmzxqZpoSk60mINpv2YLhVU9zlO\nfZk+jcIUqRuITAlb9tySdG1m2NJaTdk1m6aEZtbtIDDtRGRKmudZ22WiM+/QzDpNIjDtRWRmsPXu\n6yZ+P+9oZmnvSdn1h7aEJi8C0y1EZkbTQjPN8BuqbGyafALhrHuTJnnnqU/O9HykR2QSK7s2UyY2\neU+4rBPTJBCZCpQdzeSJSJnQ5LmMRB2m/Zxp2zktMIxi2oHIVGRSaGbd01RmzaYpoUmFwLQHkWmZ\nIsGZZ2jaHjFUh8jMUdGjgIfliU3XRzRoPiLTAdNiQ2gwT0SmQybFps7QEC0MIjJzNuuUaZRxsenK\niIZF33YhMhVq2hnPXQ4N2oPIdFzZ0JSNDZHCMCLTACmmTFUgGKgCkalJlZd+KIoLP2GeiEwDzCMC\nqc7cbtq6FOaPyPTELCFj2oRZEBkASRGZGtW1LnPjqZ/5zp9Bw6OZoh+RC5RBZDpmOCxcjwXzxkei\ndETZmBT5jKcmfsQKmo+RTA8MBojd2agbkalQk9ctJoWmydcHRvsRmQaYxxG/KUc0TKkwiMjUaN4j\nhknrNvPetrw4A7t9iEzP1DVtyjOameepFqgPkWmIOqdMLASjTkSmJk2ejgyGpu7RDLqPyPRU2fWZ\nJu9BQzMRmRrkHR3M87oyqaZNjGZAZCApXeAY+YDIdMRH9l0+7004AoGBRGR6q4lRGtaGbcR0RCax\nOvcq5X1TTntcFcevMIrBEiLTM00fHTT1ouooj8j0yLjAVP3GZhSDQUQmoXkcgDcuJHWNYFIH5gv7\nTk/6+qgekemBpk+R0G1EpoOKRKVtayCchd0+RKZBOFmR30EXEZmOY6qEeSMyDVH1/8HzxGXUVIlr\nvKBqRKbDmj6KafLlL1AdIoPG4jOjuoHINMA8FjvbtlcJ7UVkACRFZAAkRWQAJEVkEmJ3MEBkeosj\na1EXIjNnvNnHa/pxPsiHyABIisgASIrIJNbkxV+maqgDkcFhOJ8IVSMyKKyqS2w2eZSH6hAZlJLq\nWr6cU9U9RAal8akEyIPIYCaEBtMQGcyM0GASIoMjlNnDVDY0LP52H5GZMxY60XVEpgb837o8zl9q\nPyKDShz/yIHD/gaWEBnM5PhHDhwRFkKDQUQGR8g7vZsUE0KDJUQGhY0avVSJxfBuyRUZ28fZ3mR7\nt+1dts8f+N6HbIft5QP33Wh7r+0nbV+SYsPRfIxmIOUfydwm6cGIOEPSOZJ2SZLtUyS9W9L/W3qg\n7TMlbZB0lqRLJX3C9lFVbjTag9BgamRsHyvpQkl3SlJEvBIRL2Xf/nVJ10saPHprvaT7IuJQRDwt\naa+kdZVuNVqF0PRbnpHMKkkHJd1l+yu277B9jO31kp6JiK8OPf5kSfsHbi9k9x3G9tW2t9nedvDg\nwbLb3wqTjqBt24WjygaD0PRXnsgsk7RG0u0RcZ6klyX9kqT/IukXy/7giNgYEWsjYu2KFSvKvgxm\nMGqBtekHDn5h3+nz3gQUlCcyC5IWIuLR7PYmLUZnlaSv2t4naaWkx2yfKOkZSacMPH9ldh+AHpoa\nmYg4IGm/7aX/hVwk6bGIeEtEnBoRp2oxRGuyx26WtMH20bZXSVotaWuazUdZ7CZGXZblfNy1ku61\n/QZJT0l6/7gHRsRO2/dLekLSq5KuiYjXZt7SlmrTNXObPlVCO+WKTERsl7R2wvdPHbp9k6SbZtoy\nNM6si7fHP3KAa8/0EEf89lDbpkqcid1uRAZAUkSmIuuuurXwc+ZxjMy4UUzT1mPaNtrCeEQGQFJE\npgKXrb5+3pswkzyjGI7YRVlEZkaTAtOm3ddAKkRmBm0bwZQ9jYBRDGZBZEoaDEyZYz/admIk+sX2\n99jeavurtnfa/m/Z/W+2/ZDtPdnfb5r2WkSmhLaNYMZp2h4lNMohSe+MiHMknSvpUts/JOkGSQ9H\nxGpJD2e3JyIyBXUlMHkxVeqnWPTX2c3vzv6EFq8XdU92/z2Srpj2WkQGtSJajbF86XpO2Z+rhx9g\n+yjb2yU9L+mh7EoMJ0TEs9lDDkg6YdoPynuCJMbo+rk4L15wImGo2KG/fUMTDjZ8ISLGno8oSdmJ\nzefaPk7SA7Z/YOj7YXvqLlRGMgVt2XPLvDehtOHFZnaxI4/scrt/pMVrdj9n+yRJyv5+ftrziUwJ\nS6GZNIrp0iU3uz5aw5Fsr8hGMLL9vZLeJWm3Fq8XdWX2sCslfXraaxGZkto8ohnEaAZjnCTpj2zv\nkPRnWlyT+UNJN0t6l+09ki7Obk9EZGaw9e7r5r0JhY0aReUJDaOZfomIHRFxXkScHRE/EBG/nN3/\nYkRcFBGrI+LiiPjGtNciMjNqY2hGYUSDVIhMBfoSGkYzKIPIVKRNoWnbwjPajcjUrO9vcEZD/UNk\neqrvsUN9iAxqwyimn4gMgKSITAJt2R3c1ikTn4fdLkQGtWCq1F9EpkZtHTnMisD0G5EBkBSRAZAU\nkQGQFJEBkBSRqUlfF30BIgMgKSIDICkikwAfmga8jsgASIrIJNCWc5eAOhCZHqvrA8b4cLh+IzJo\nhY/su3zem4CSiAwaqwEf5YoKEJmeGvcGTrWexJSpv4hMDzFCQJ2IDICkiEzPMIpB3YhMRdZddevY\n73Fy5CLWZfqJyFRgMDBNPhAvzyimaadEMPJqPyLTE1UEpoqRSFWjGT6xoD2IzIzaMIpp4wgG3UFk\nZjBpHWZJV9Zj6lpPmRQ7jvptJyJT0nBgUo5ibjz1M6Wf28RRDAvA/UJkSsgzgqlK1wJTBou/7UZk\nCupjYBh5YBZEpgLjpkpNXo+Z9wiGcPXHsnlvAKo3bRRTJDCjYrBlzy1jH3/Z6utzv/asvrDvdL3z\n1Cdr+3koh5FMg6XYm5IyMHm+j/4hMh0zaRSTOjBFH4d+IDId0oTAlH18XsOjO478bT4i03Af2Xf5\nzNOmugMz6/PQLURmRk3ZszRuFDOvwFT1fLQfkWmAG0/9zEzHxIwz78BU/TpoJyJTQOoD8crGZtQo\npimBAYhMTnUe6VunOg6KG3fszLifPS2QnCjZLkQmgSLrMbNOk2YZxYx7kzOKQZWITA7jRjFNvH5M\n0wJT1RHAnCTZXkRmipTTpFGjmCJTgbJvvLYFZhSmTO1BZCYo8yapa9d12WlSWwIz7xM4UR0iM0ad\nJ/rVpS2BGWc4rIxm2oHIjLD0JnnxghOT/YyqrxVT9iLgbQkM2ovIlDRq0Xde149pSmCqxpSpG4jM\nHKS+4l2XMWVqHyIzRqqp0qTAlHnDdHUUg+4gMhXJM1VKceBdGQQGdSIyJZQ5CC/FCZBltC0wrMu0\nX67I2D7O9ibbu23vsn2+7V+xvcP2dtufs/3WgcffaHuv7SdtX5Ju89NIuVdpnGlTpVlPgpTaFxh0\nQ96RzG2SHoyIMySdI2mXpI9GxNkRca6kP5T0i5Jk+0xJGySdJelSSZ+wfVTlW94gTf1UgsH1GAKD\neZkaGdvHSrpQ0p2SFBGvRMRLEfGtgYcdI2lpDrFe0n0RcSginpa0V9K6aje7XaZNlVKPYpoamLIj\nxr7vYWubPCOZVZIOSrrL9lds32H7GEmyfZPt/ZLep2wkI+lkSfsHnr+Q3dcJVZ8UWSYwRTQ1MEWw\nLtNueSKzTNIaSbdHxHmSXpZ0gyRFxIcj4hRJ90r6uSI/2PbVtrfZ3nbw4MGCm90e87riXRs+PK0N\n24jZ5YnMgqSFiHg0u71Ji9EZdK+k92ZfPyPplIHvrczuO0xEbIyItRGxdsWKFcW2ukEmrcekmCbl\nsfTmnccopgsjJ1RramQi4oCk/baXPnviIklP2F498LD1knZnX2+WtMH20bZXSVotaWuF2zw3RaZK\nqXZZ5z34bp5v9iI/m9FM9+Xdu3StpHtt75B0rqRflXSz7cez+94t6YOSFBE7Jd0v6QlJD0q6JiJe\nq3zLWy7lYm8TRhOpt4HF3/bI9VnYEbFd0tqhu9876rHZ42+SdNMM29U4RU6InNfepOMfOdCIwCzZ\nsueWys7KXv6YG3klQkzHEb8Vm9eRvU0LDLCEyFQoT2BSjGKaHJiU28WUqR2ITA7zHKa3OTBLmr59\nSIvIlFTmVIKqd1m3ac8MoekvIjNF3gXfuj8/qU2BQb8RmZoUvSAVgUFXEJmC6pgmERh0CZGZoI4F\nXwKDriMyBVQ9iiEw6AMiU4Fxi75VXUmfwKDNiMwYw1OloqOYqtZhCAzaLte5S32y7qpbk/+MItMk\niWNM0G6MZAaUDczwqIV1GOB1RCYz6whmKSxVB4ZRDNqOyKi6KVJVgVlCYNAFvY/MqMDM+7olbTjp\nsSk4E7v5eh2ZvCOYWT9Xqeg0icCgS3odmTqwJwl919vINHFX9da7r0u5OcBc9DYy41S1HkNggEVE\nZooqPueaT0BcxFSwn4jMgFSjmGn6dNDdLKEh1u1EZAqadgW8otOkPgUG/URkJhieKhW9xCaBGW3S\naKavv5MuIzKZKj6ClgPD8is7bWLK1D5EZoyioximScWNCw2/m24hMqr+NAICkx97nLqPyFSgyDSJ\nwMyOKVO79D4ydY9ikA8x7o7eR2aUUQfgVXW93lGYMqDLeh2ZoqOYURemmnWPEoEZb9JohhFje/Q6\nMmVUOaIhMEfid9I9vY3MuFFMFecqTfLiBSdK4s00yeDvhrWZ9uttZKpQ5pKaEoHJI8/viClTOxCZ\nmnFJh/yWQsNopt2IzIDUUyUUx6iv/YhMhabtrWIUUz2mTM1HZGpCYNBXRKYGBAZ9RmQyZdZjRj1n\neMpEYNB3RCYhAlMd9jC117J5b0AXvbAm9NQHPjTvzQAagZFMAgQGeB2R0WzHxww/l8AAhyMyaDwO\nyGs3IoNW2LLnFhZ/W4qF3xyGLyKe8gJWQNcQmRGKfr4S6rFlzy1ad9Wt894MFMR0Ca3CsUftQ2RK\nYKQzX4SmXYjMEALSDoSmPXoZmbd9/GPz3gRUgNC0Q+8iw8IhUK/eRQZAvXoVGUYxQP16ExkCA8xH\nLyIzHBiuCwvUpxeRAVCM7VNs/5HtJ2zvtP3B7P43237I9p7s7zdNey0iA2CUVyV9KCLOlPRDkq6x\nfaakGyQ9HBGrJT2c3Z6IyAA4QkQ8GxGPZV//laRdkk6WtF7SPdnD7pF0xbTX4gRJoGbL/mb+64JP\nS8ttbxu4a2NEbBz1WNunSjpP0qOSToiIZ7NvHZB0wrSf1cuRzKQPYeMyDuiJFyJi7cCfcYH5B5J+\nV9J/jIhvDX4vIkLS5E80VA8jM+1THgEssv3dWgzMvRHxe9ndz9k+Kfv+SZKen/Y6vYsMuoWr5aVh\n25LulLQrIgaPAdks6crs6yslfXraa/UiMksn0o0axcxyEXGgwy6Q9NOS3ml7e/bnckk3S3qX7T2S\nLs5uT9Sbhd+td1/H2ddAThHxfyWNW52+qMhr9WIks4SPKwHq16vISIQGqFvvIiNND03R3dis6wDj\n9TIyVeAynUA+RAZAUkQGQFK9j8zXd6ys5HVYlwFGyxUZ28fZ3mR7t+1dts+3/dHs9g7bD9g+buDx\nN9rea/tJ25ek23wATZd3JHObpAcj4gxJ52jxtO+HJP1ARJwt6WuSbpSk7JoTGySdJelSSZ+wfVTV\nGw6gHaZGxvaxki7U4nkMiohXIuKliPhcRLyaPexLkpbmHesl3RcRhyLiaUl7Ja2rftPnjz1M8/fi\nBSfOexMwRZ6RzCpJByXdZfsrtu+wfczQY35G0pbs65Ml7R/43kJ232FsX217m+1tBw8eLLHpaZW5\n5MNpZy9w6gIwJE9klklaI+n2iDhP0ssauOSe7Q9r8VJ99xb5wRGxcelaFitWrCjy1MpVtfi7hNAA\nr8sTmQVJCxHxaHZ7kxajI9tXSXqPpPdlF7CRpGcknTLw/JXZfUDluNRD802NTEQckLTf9unZXRdJ\nesL2pZKul/SjEfE3A0/ZLGmD7aNtr5K0WtLWire7cqNGM2WnTBKjGWBJ3r1L10q61/YOSedK+lVJ\n/13SGyU9lF1r4pOSFBE7Jd0v6QlJD0q6JiJeq3zLG2LS4i+hAXJGJiK2Z+snZ0fEFRHxzYj4xxFx\nSkScm/359wOPvykiTouI0yNiy6TXnpdRF3Kuem1GIjRA74/4nWaWKdMSQoM+63Vk6hrNSHwWN/qr\n15FJaXg0w6ckoK96H5k8o5lJUyY+p2l+Llt9/bw3ATn0NjJLn2AglQ8NgZkfAtMevY2MdHhoisoT\nGC7/kAaBaZfefCTKOFvvvk7rrrpVyx/zEesmX9+x8rBQMHKZPwLTPr0eySxZGtGk2NvEaKY6BKad\niEymitBw6QfgSESmIkuBGRUaRjOzYxTTXkRmQMppk8SRv2URmHYjMkNShwboGyIzwqTQlMUlINBX\nRGaMccfQMG0CiiEyE2y9+7pKRzODCA36gshMMSo0Va3NEBr0AZHJIc+IpuzRwFwCAl1HZHIaXqMZ\nN5qZFJtxx8sQGnQZkSkg1foM0GVEZgZVrM1wMSt0HZHJaWlKM2k0w+UfgCMRmRxYMwHKIzJTpAgM\n1/9FnxCZCfIEhstAAJMRmTGqHMGMugwEazPoCyIzQpWXFmAEg74jMgmNCsy40QynGKCriMyQqkYx\njGCARURmwFJgXrzgxLGPSbkniNEMuojIZMqOYEYt4M4yiuGYHHQNkVH+wNR1PAuhQZf0PjLDgZk0\nVRpW9ShmEKFBV/Q6MkWmSPM4KpfQoAt6G5lRgWnCKGY4ZoQGbdfbyBQx73OLCA3ajMhkioxiZsXp\nBegTIjMsisW3AAASPklEQVTFPEYx8x45AVUiMqp3FJMXoUFXEJkJ6nqjM2VCl/U+MlWOYsp+LMo4\njGbQBb2PzDh1v8HHjWYIzWRNnOricL2ODP+BAun1NjKznmld5sPdBo06cI/RDLqot5FpKhaB0TW9\njczxjxwY+73lj3mmT4ucZTSD/JjutkNvI7Nlzy0TQyNNj82kTypIMW3C6whMe/Q2MlK+0EiTY0No\n6kdg2qXXkZHyh0Ya/xG1qUKDIxGY9ul9ZKTmhGYYo5nDEZh2IjKZJoSG0cx4BKa9iMyAJoaG0Qza\njsgMWQpN3gXhUfJ+PjbyYRTTbkRmghShKbs+w2cyoa2IzBR5p091IDRoIyKTQ56D9kZJMW0iNGgb\nIpNTlaGZ9bozhAZtQmQKaNKIhk8wWMQZ6s1HZAoqGxpUj8C0A5FJYFRoioxm8h6Ux2gGbUBkSmjS\nHqe+hoZRTHsQmZLKTJsGRzNFFn9HHfXLmwxtQWQK2LLnFm3Zc8t3bqcc0XB6wXgEtl2ITAmDoZkk\n9SIwbza0AZHJaTgsS7frXJ9h2tS/f28XEJkcxo1c8oRmltEMV8xDFxCZKaZNjfJOnZaUXfwFo5i2\nIjIT5A3ItOvQVL02w2gGbUJkxig6Qin6+Ly4Wt4iRjHtRWRGKBuMrXdfV/GWLGJtZjx+D81HZFqC\nEQ3aishULNVoZpS+/F+cqVK7EZkEioSmyB4mpk2H6/O/vU1yRcb2cbY32d5te5ft823/hO2dtr9t\ne+3Q42+0vdf2k7YvSbPpzTYpNLNcX4ZpE9om70jmNkkPRsQZks6RtEvS45J+XNIXBx9o+0xJGySd\nJelSSZ+wfVRlW9widU6dgKaaGhnbx0q6UNKdkhQRr0TESxGxKyKeHPGU9ZLui4hDEfG0pL2S1lW5\n0W2SIjSMZpgqtUmekcwqSQcl3WX7K7bvsH3MhMefLGn/wO2F7L7emiU07zz1Sb3z1FEtB9ohT2SW\nSVoj6faIOE/Sy5JumPUH277a9jbb2w4ePDjryzVeyjOy3/XH/ynZawOzyhOZBUkLEfFodnuTFqMz\nzjOSThm4vTK77zARsTEi1kbE2hUrVuTd3k6q4hwmQoOmmhqZiDggab/t07O7LpL0xISnbJa0wfbR\ntldJWi1p68xb2jF8JhP6Iu/epWsl3Wt7h6RzJf2q7R+zvSDpfEn/2/ZnJSkidkq6X4shelDSNRHx\nWvWb3i2TRjNf2Hf62O/1EYu+7ZIrMhGxPZvanB0RV0TENyPigYhYGRFHR8QJEXHJwONviojTIuL0\niNiSbvO7r2hgGM2gaTjid45STJkkQoNmITINMuspBoMIDZqCyHQYoUETEJkaTPoAtuEp0yy7s1kQ\nRRMRmcRSfsIjpxegDYhMjVJ/DhPQREQmoTo+p5rRDJqOyCRSR2BGYV0GTUNkEigamFmPl2E0gyYj\nMhW7bPX1E7/Pugz6hshUaFpgyip6PRmmTGgSIlORVIHJiykTmorIVGDegZH4XG00F5HpoFQnXgJl\nEJkZFR3F8EFl6BsiM4Oqpkks1KLLiExJ4wLz4gUnJvl5fGIB2orIlJByoZe9REcaPraINad2ITIF\n1b0nKc8Ihj1LaDIiU8CsgWHRF31EZHJqwrEwQBsRmRxSBoY9S+g6IlOhWfcslV30ZbEYTUZkpmAU\nA8yGyEzQtHUYjpVBXWx/yvbzth8fuO/Nth+yvSf7+015XovIjFE0MNOmStP2LI2a8hAVzNHdki4d\nuu8GSQ9HxGpJD2e3pyIyI9QxgmGqhCaLiC9K+sbQ3esl3ZN9fY+kK/K81rIKtwtAeyy3vW3g9saI\n2DjlOSdExLPZ1wcknZDnBxGZBmDvUL8s++u/1/GPHJj3ZrwQEWvLPjkiwnauo0uZLo2wZc8t894E\n1mPQRM/ZPkmSsr+fz/MkIjNGytCwHoOW2izpyuzrKyV9Os+TiMwEVYVm0p4l9iqhiWz/tqQ/lXS6\n7QXbPyvpZknvsr1H0sXZ7amIzBRNmDrhSFzuIa2I+MmIOCkivjsiVkbEnRHxYkRcFBGrI+LiiBje\n+zQSkcmhztCUGcVwqQc0GZHJqarQDK7HsFcJfUBkEuMaMug7IlOBqq7rW3bBd9KI6G0f/1jZzWmM\nrXdfN+9NwAyITE5NOVmyTIjWXXVrgi2p11Mf+NC8NwElEZkcmhKYWXQhNGgnIjNFXYGp8tiY4YP9\nltaF2h6a4dEMu7HbgchM0KaPPsn7el0LDZqPyIxRRWBS7VmaddRDaFAnIjNCkcAU2bM0z3OWurYr\nndC0B5HpqbaPZiRC0xZEpqVGTZn6eAQxi7/NR2RmVOTiQ1W/ISatzXA5CTQFkUlo+IPiB3FSI/qC\nyDTEF/adXup5w6OZPk6Z0GxEpgINuF5rLl3bw4R2IDIdwJX00GREJrFJ6zKpsfiLJiAyFZn3lGlw\nNNOndZl5Rhz5EJmaTdqNXXbxt6wuHJDHtWaaj8jUYNz/bavcjZ03UMOLv10IDZqNyFRo0pRpMDSD\no5nh0FQ9msmzLkNokBKR6aFRu7IJDVIhMnPStNGMRGiQBpGpWN4pk1TduUxlosSBeagLkWmIqs9l\nGt6NzWgG87Js3hvQN8sf82GjiK/vWDkyAF/Yd3ryI3lfWBMj93yNCs3gCI2P7kURRGaEsm+ipSvq\nHf/IgVKfxfSRfZcnPZDutLMXjpiijQvNsBcvOPE7obls9fWEBrkxXarQ4BuvirWZWReA8wYr7/rM\nixec+J0/TKuQF5GpWN7QjDOvPU1lFoIJDfIgMgnkmUqk2tOUB6FBnYhMIkuhKTJtWlLlaKboGg+h\nQdWITEJ5QpNX1aGZtEv7hTVRODaEBuMQmZqMC03ec5qkxdBUuUYz7dgZQoMqEJnEUuzqLRObcdOm\nqkMDDOM4mRps2XOLLlt9/djjZwYP0Bs8OG/acTN1HLBX1LTRDNd/6R9GMjVJdfBa3hHNuNMWpu3V\n4spzmBWRqdGk0BRZm6lK1YFZWjBmioVBRKZmeacLVYZm1PNTBGb4NrGBRGTmosy6RNnQFA3M8sc8\nc2CGv0ds+i1XZGwfZ3uT7d22d9k+3/abbT9ke0/295sGHn+j7b22n7R9SbrNb69RoZl2FHDqj7Yt\ns/6SNyDEpr/yjmRuk/RgRJwh6RxJuyTdIOnhiFgt6eHstmyfKWmDpLMkXSrpE7aPqnrD+2KW0w3y\njmLKjF4kdm8jn6mRsX2spAsl3SlJEfFKRLwkab2ke7KH3SPpiuzr9ZLui4hDEfG0pL2S1lW94V2Q\ncn0mT2BmiQuBQV55RjKrJB2UdJftr9i+w/Yxkk6IiGezxxyQdEL29cmS9g88fyG77zC2r7a9zfa2\ngwcPlv8XtNxwaPK86aeFJk+I2DWNuuSJzDJJayTdHhHnSXpZ2dRoSUSEpEL/a4uIjRGxNiLWrlix\noshTO29UAPKsz3xk3+Wlj4cpYpZAve3jH6tsO9AOeY74XZC0EBGPZrc3aTEyz9k+KSKetX2SpOez\n7z8j6ZSB56/M7kMBw5fplI68VGfeqdOoaVKK7cvrbR//mJ76wIdm3ga0w9SRTEQckLTf9tKhpRdJ\nekLSZklXZvddKenT2debJW2wfbTtVZJWS9pa6Vb3RJ4RzTQpAlPFazGi6Y+8e5eulXSv7R2SzpX0\nq5JulvQu23skXZzdVkTslHS/FkP0oKRrIuK1qje8L2YJTcrADL5m2dclNP2Q6wTJiNguae2Ib100\n5vE3Sbpphu3CgKU38bRPOahi3WX4khR5L4hedvpEaLqPs7AbIO+nG4z6OJW88ow2Rl3zZum+MtsH\nSJxW0Bh5r55XZmoy7TnHP3Jg6s/P85g8Pwv9Q2QaJMUbOU9gisizjYQGg4hMA1UVmkmP2Xr3dZVc\ne7jMz0a/EJmGKhOApT09eff4zHIhrVRTJxaCu4fINFieaUmRqCyp6hKYqULDBcm7hcg0XNVTmuHA\n1PGZ1mWuT0NouoPI9Mi4EUzqaVMZhKY7iEwLVPFGnjZFauL6jMTUqQuITEvMEpq8azCpp05lL+tJ\naNqtF0f82h3Znbq33NN8Tw1nPOfdtnumP2TQC2tCyx8zoWkxRjJoPE5VaDcig1YgNO3V2enS4sX6\nUMZlq68v/dw8J1KWDcbS1AntwkgGrcKIpn2IDICkiAwqVfXBeaedvXDExbkYzbQLkUEjDceF0LQX\nkUHtpi3eDgdl3P2Eph2IDI5Qx0mTZRGa9iEyqFzKi2FJ40c6aCYig5GaPJqRCE2bEBmMlfLM7CoO\nqiM07UBkMBEjGsyKyGCqsqFJvTaDdiAyyKXpIxo0F5FBbmVCM2k0w8mO/UBkUEjVoUH3ERkUVuXU\nidFM9xEZ1KLIaObrO1Ym3BLUjcigFEYzyIvIoLSioWFtpp+IDGZS1YhmeDTDlKk7iAxmViQ0jGb6\nh8igEhysh3GIDCoza2iYMnUTkUHtmDL1C5FBozCa6R4ig7lgNNMfRAZzMy40jGa6hcigFQhNexEZ\nzFXe0YxEaNqKyABIishg7hjNdBuRAZAUkUGlqr7oOKOZ9iMyqFwd5zERmvYgMkiCi1phCZFBMqkv\nasVoph2IDJJKdVErtAeRQXIp12gYzTQfkUEt8oaGD4PrHiIDICkig9owmuknIoNacS3g/iEyqF2e\n0DCa6Q4ig7mYNTRoDyKDuZll6sRopj2IDOZqWmgYzbQfkcHclR3RMJppByKDxmM0025EBkBSRAZA\nUkQGrca6TPMRGQBJERkASREZAEkRGQBJERkASREZAEkRGbQCR/22F5EBkBSRAZBUrsjY3mf7z21v\nt70tu+8c23+a3f8Htr9/4PE32t5r+0nbl6TaeADp2L40ew/vtX1D2dcpMpJ5R0ScGxFrs9t3SLoh\nIv6ppAck/Xy2YWdK2iDpLEmXSvqE7aPKbiCA+mXv2f8h6TJJZ0r6yey9Xdgs06W3S/pi9vVDkt6b\nfb1e0n0RcSginpa0V9K6GX4OgPqtk7Q3Ip6KiFck3afF93Zhy3I+LiR93vZrkn4zIjZK2pn90N+X\n9BOSTskee7KkLw08dyG77zC2r5Z0dXbzkO3Hi29+MsslvTDvjRjA9kiL/7sa7fQat2Jm3zr03Gcf\n3PvR5XPejO9ZWvrIbMze10tOlrR/4PaCpH9e5gfljcwPR8Qztt8i6SHbuyX9jKSP2/4FSZslvVLk\nB2f/oI2SZHvbwDRs7tieyZq4PfPehiIi4tJ5b0Odck2XIuKZ7O/ntbj+si4idkfEuyPiByX9tqSv\nZw9/Rq+PaiRpZXYfgPao7H08NTK2j7H9xqWvJb1b0uPZqEa2v0vSf5X0yewpmyVtsH207VWSVkva\nWmbjAMzNn0labXuV7TdocWfO5jIvlGe6dIKkB2wvPf63IuJB2x+0fU32mN+TdJckRcRO2/dLekLS\nq5KuiYjXpvyMjVO+Xze2ZzK2p+Mi4lXbPyfps5KOkvSpiNhZ5rUcEZVuHAAM4ohfAEkRGQBJ1RIZ\n25+y/fzgsTC2f8L2Ttvftr124P5Tbf9tdgrDdtufHP2qlW/PR23vtr3D9gO2jxv4XtLTJIpszxx/\nP7+Sbct225+z/daB7yX7/RTZljp+NyghIpL/kXShpDWSHh+4759o8SCqP5a0duD+UwcfV+P2vFvS\nsuzrX5P0a9nXZ0r6qqSjJa3S4q76o+a4PfP6/Xz/wNcfkPTJOn4/Bbcl+e+GP8X/1DKSiYgvSvrG\n0H27IuLJOn5+zu35XES8mt38khaPC5BqOE2i4PYkN2Z7vjVw8xgtHgUuJf79FNwWNFBT12RWZcPd\nP7H9I3P4+T8jaUv29ajDq484TaLG7ZHm9PuxfZPt/ZLeJ+kXs7vn8vsZsy3S/P/bwZAmRuZZSf8w\nIs6VdJ2k3xq8jERqtj+sxeN77q3rZ04yYnvm9vuJiA9HxCnZtvxcHT+z4LbM9b8djNa4yGTD7hez\nr7+sxTn+2+v42bavkvQeSe+LiKUh+NxOkxi1PfP8/Qy4V6+fdT/v00i+sy0N+d1gSOMiY3vF0vVn\nbL9Ni6clPFXDz71U0vWSfjQi/mbgW3M5TWLc9szx97N64OZ6Sbuzr2v//Yzblnn9bjBFHavLWjyB\n8llJf6/FOfvPSvqx7OtDkp6T9Nnsse/V4mUktkt6TNK/qml79mpxbWF79ueTA4//sBb/r/ikpMvm\nuT1z/P38rqTHJe2Q9AeSTq7j91NkW+r43fCn+B9OKwCQVOOmSwC6hcgASIrIAEiKyABIisgASIrI\nAEiKyABI6v8DkaCEZ4oMBIgAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1,1, figsize=(16,12))\n", "bathy = ax.contourf(grid.variables['Bathymetry'][:], levels = np.arange(0,60, 10))\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))\n", "fig.colorbar(bathy, ax=ax)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": true }, "outputs": [], "source": [ "result = nc.Dataset('/ocean/vdo/MEOPAR/ariane-runs/Apr17/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": 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(720)\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": 9, "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": 12, "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": 14, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf4AAAGDCAYAAADK03I6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XecFdX9//HXZ3sv7MKyLEtfgaVKExEV7BoNxCRGExWN\nLYktiUks3yQmUROT/GKiMbaowV5i7IrdtVIEFOm9t6XDUhaWPb8/ZsDLuuUCe3fu3n0/H4/72HvP\ntM+ZO3s/M2fOzJhzDhEREWkZ4oIOQERERJqOEr+IiEgLosQvIiLSgijxi4iItCBK/CIiIi2IEr+I\niEgLosQvmNlYM7s1oGWbmf3HzDaZ2aSAYqgwsy4NjNPJzJyZJURg+Z+Y2ZH++5vM7MEwpwt73Dqm\nH2FmKw51+gbmHbH1dRAxFJjZbDNLrmec+8zsN00ZV6SEsx03FTM71szmBh2H1E6JPwqZ2RIzKzez\n9JCyS82sLMCwImU4cDLQ3jk3pOZAM7vIzPb6P2pbzewLMzvzUBdmZmVmdmlomXMuwzm36FDneTjM\n7Cxgm3Pucz+WPzrnLm1gMmqOG3Si9bfZkwJa9klmNtXMtpvZCjM7B8A5txZ4H7i8rmmdcz9yzt1y\niMstM7NdZrbN3zanmNkN9e1oRFKktuNwDgz8ba9bSCwfOee6N3Ys0jiU+KNXPHBt0EEcLDOLP8hJ\nOgJLnHPb6xlnvHMuA8gBHgKeNbPcg4zLzCwat/cfAY8FHUTQzKy1mdkhTFcKPAn8H5AN9AOmhIzy\nBHBFowRZu6ucc5lAIXAdcC7w+qHURaTJOOf0irIXsAS4AdgI5PhllwJl/vtOgAMSQqYpAy71318E\nfAL8HdgMLAKG+eXLgXJgTMi0Y4H7gLeBbcAHQMeQ4T38YRuBucA5Naa9F3gd2A6cVEt92gEv+9Mv\nAC7zyy8BdgF7gQrg97VMexHwccjndL/ug4Bc4FVgHbDJf9++xjq5zV8XO/GSwF5/mRXA3f54Dujm\nv08F/gYsBbYAH/tlB6xzvCTzELAaWAncCsT7w7r563ALsB54po7vOcmPKzTm3wGP1/iexwDL/Hn9\nXx3jLvPHrfBfR9eyvFT/+9oEzAJ+Cayo8T39z1+fi4FraizrOeAZfxuZCvTzhz0GVPt1qQB+1VDs\ntcR2Pd52/3ug80H8rzwJ3FLP8ARgByHbc43hY4Fb/fcjgBV4Cbzc/24vrmfeZfj/cyFlHfzlnel/\nTgb+AazyX/8Akmss71chyxsNnAHMw/t/uSlk3kOA8Xj/06uBu4GkkOGh2/FY4F/Aa/73NRHoWk9d\n/guswdtmPwR6+eWXA3uA3f53+0ot037oL3u7P8739tWtxm/aL4Ev/fEeAgqAcX587wC5IeMPBT71\n6zoNGNFYv6/+/If522Sx/7kf3v9Fj8ZcTrS+Ag9Ar1q+FO+f5CTg+ZAfpYNN/FXAxXgtB7fi/fj+\ny/8hOsX/Z8vwxx/rfz7OH34nfrLFS7TL/XklAEf6/zClIdNuAY7Ba0FKqaU+HwL3AClAf7zEckJI\nrB/Xsy4uCoklAa8VZBte4s0Dvg2kAZn+j9eLNdbJMqCXP20itf9Yh/5g/ssfp8hfd8P8dXLAOgde\nAO73108bYBJwhT/sKbwj0Di/zsPrqFsvYHuNst/x9cT/b7yk3Q+oBHrWM25CPevyduAjoBVQDMzA\n/3H2Y50C/BZvh6QL3g7jqSHL2gN8x1+Pv8DbOUgM3WZDllVv7HXENxRvJ3IDXhP9BUBaA/8ri4Bb\ngOl4yfBxoFWNcb4EvlnH9GM5MPFXAX/w63gGXhLPrWPar21LIdv7n/33fwAm+NtIa7xkdkuN5f3W\nX95leP8bT+Jtz73wdqY6++MP9NdRgr9+ZwM/rWM7HuuvxyH++E8AT9ezHn/oL3PfjsoXta2jeqbf\nv+yQutVM/BPwkn0R3o7OVLzfkxTgPeBmf9wiP/Yz8LbLk/3PretY9qt4Owi1vV6tJ+bb/OWm+tvP\nVeH8PsfCK/AA9KrlS/kq8ffGS6qtOfjEPz9kWB9//IKQsg1Af//92NAfBSAD78i4GG/v/aMa8d0f\n8k86Fni0nroU+/PKDCn7EzA2JNaGEn+V/0+83v/x+Fqrgj9uf2BTjXXyhxrj7F9PIWUO7yg9Du+H\ntl8t896/zv0fr0ogNWT4ecD7/vtHgQcIOZKvI95jgDU1yn7H15N5aIvAJODcesatL/EvAk4L+Xw5\nXyX+o4BlNca/EfhPyLImhAyLw0u0x4Zus7Wsr1pjb2C9JAPn4LUibQQerGfc3f6yj8Dbbv8HPFFj\nnE+AC+uYfiwHJv6dHPh/VQ4MrWPar21LfvnTwL/99wuBM0KGnYp3ait0eftaijL9dXZUyPhTgNF1\nLP+nwAs1t+OQej0YMuwMYE5D694fN8efV3bNdVTPNOEk/h+EfP4fcG/I56vxd9rxWn8eqzH/Nwlp\npWyMF97O1hS8pP8GYI05/2h+ReM5T/E552bg7c3ecAiTrw15v9OfX82yjJDPy0OWW4H3g9sO7xz8\nUWa2ed8L+AHQtrZpa9EO2Oic2xZSthRvrz5cE5xzOc65fOfcUOfcOwBmlmZm95vZUjPbineklVOj\nn0F9sdWUj3f0sbCB8Tri/WisDlkn9+Md1YHXdGvAJDObaWY/rGM+m/B+7BuyJuT9Dg783g5GOw5c\nH0tD3ncE2tX4nm/C28nZJ3QbqcZrpm7XwDIPOnbnXCXeUfoXeIm9dz2j78TbOZnnb7d/xEtyoTLx\ndhzDscE5V3WwMddQhPf/A976CV3PSzlwnW1wzu313+/0/9b6f2pmR5jZq2a2xt/e/4i3zdYlrHVv\nZvFmdruZLfTnu8QfVN+8D0XNetX1e9QR+G6NbXE4Xj+KRuOc24O3U9Mb+Jvz9wZaAiX+6HczXhNg\naKLc1xEuLaQsNBEfiuJ9b8wsA685eBXej/0HfuLd98pwzv04ZNr6/mFWAa3MLDTBdcA7L364rgO6\n4x0hZeGdqgAv6dYVW32xrsc7/9+1geUuxzvizw9ZJ1nOuV4Azrk1zrnLnHPt8DqW3RPa4znEArx+\nhwezE1SXcH60VhPyPeN9D/ssBxbX+J4znXOhSTR0G4kD2uN9v+Euv15mlmdmV/mXdb6Hd6plpHNu\naD2TfVlj2QfE4V/l0A3vPHHEmVkxXpP8R37RKrxEtk8HvlpnB+teYA5Q4m/vN3Hgtn6ovg+Mwmtl\nzMZrrSFk3k2dEJfjHfGHbovpzrnbaxvZzMb5V/3U9hpX10L8/7ubgf8AfwvqaowgKPFHOefcArwO\nVdeElK3DS5zn+3vrP6ThZNWQM8xsuJkl4Z0zneCcW47X4nCEmV1gZon+a7CZ9Qwz/uV45zX/ZGYp\nZtYXr1Pf44cZL3hHcjuBzWbWCu+fuCFr8c5f1xZrNfAwcIeZtfPX7dE1fxCcc6uBt/B+LLLMLM7M\nuprZ8QBm9l0za++Pvgnvh7O6luXtxuvUdHw4lW3AOn8Z9V3H/Sxwo5nl+vFdHTJsErDNzK43s1S/\n7r3NbHDIOAPN7Gw/mf4Ub+dngj+szvUaDjO7BO9I83i8Dn7FzrnrnXOzG5j0P8DFZtbFzNLwWsde\nDRk+BK9pfWmtUzcSv/XpeOAlvHX5uj/oKeDX/lUL+Xjn8w91288EtgIVZtYD+HED4x/MfCvxTv+l\n4bUkhArnuz2s77+Gx4GzzOxUfztMMe+eE+1rG9k5d7p/MFLb6/TapvGvuhiL18nwEryd4kO6rLM5\nUuJvHv6A14ks1GV4vWQ34HUC+vQwl/EkXuLciHfEcj6A30R/Ct5lSqvwmg//jHceNlzn4R1FrMLr\nFHfzvub6w/QPvI45+879vxHGNHcC3zHvhkF31TL8F3jn/D7DWxd/pvb/kwvxOsHNwkvuz/FVU+Rg\nYKKZVeBdzXCtq/v66vvxOrEdFufcDvwrGPzm0dqOkn+P19S8GG/H5bGQ6fcCZ+L1k1iMt04fxDsC\n3OclvD4fm/yYz/abS8Hrt/Frf9m/OIQqjMfref9d59xrIc3f9XLOPYzXp2KiX7dKQnaS8U5L3XcI\n8YTrbjPbhpf4/oF37vo0fycSvI61k/FaJqbjdWg71Jtl/QLv6HwbXqfJZw4j7lCP4q27lXjb84Qa\nwx8CSv3v9sU65vE74BF/nHMOJxj/YGEUXovGOrwWgF/SuPnqGrxTc7/xm/gvxtuBPLYRlxG1rAWd\n1hCJSmb2CV6P4s+DjqUuZvY7vM5b5wcdS7jMrA3eZZVHOud2BR2PSLQI7HaaIuJxzh0TdAyxyDlX\nDoR1SkqkJVFTv4iISAuipn4REZEWREf8IiIiLYgSv4iISAsSs5378vPzXadOnRptftu3byc9veYV\ndbEj1usHsV/HWK8fxH4dY71+EPt1DLJ+U6ZMWe+ca93QeDGb+Dt16sTkyZMbbX5lZWWMGDGi0eYX\nbWK9fhD7dYz1+kHs1zHW6wexX8cg62dmYd2oSk39IiIiLYgSv4iISAuixC8iItKCKPGLiIi0IEr8\nIiIiLYgSv4iISAuixC8iItKCKPGLiIi0IEr8IiIiLYgSv4iISAuixC8iItKCxOy9+hvT1GWb+HJd\nFW5uOblpSfRrn42ZBR2WiIjIQVPiD8Ptr89h0pJKmPIZAEcUZPDQmMEUt0oLODIREZGDo6b+MNz2\nrd78ZmgKL/xkGH/5Tl9Wb97FKX//kAc/WsTWXXuYu2Zb0CGKiIiERUf8YSgpyGRlTjxHdsj1XsU5\n/P6VWdz62mxufW02AHd//0jO7Nsu4EhFRETqpyP+Q1BSkMljlwzh/gsG8o0+hQDc/NJMtu3aE3Bk\nIiIi9VPiP0Rmxqm92vKvHwzgpSuPYcP23fz5jTk454IOTUREpE5K/I2gX3EO5w3pwOMTllE2d13Q\n4YiIiNRJib+R/ObMngDMWLkl4EhERETqpsTfSNKSEijMTmHx+u1BhyIiIlInJf5G1KV1OouU+EVE\nJIop8TeijnnpLN2gxC8iItFLib8RFeWksmnHHnbt2Rt0KCIiIrVS4m9EbbNSAFi9ZVfAkYiIiNRO\nib8RFWbvS/w7A45ERESkdkr8jagwJxWANTriFxGRKKXE34jU1C8iItFOib8RpSbFk5GcwIaK3UGH\nIiIiUisl/kaWl5HE+orKoMMQERGpVcQSv5k9bGblZjajRvnVZjbHzGaa2V9Cym80swVmNtfMTg0p\nH2hm0/1hd5mZRSrmxpCfkcyG7Ur8IiISnSJ5xD8WOC20wMxGAqOAfs65XsD/88tLgXOBXv4095hZ\nvD/ZvcBlQIn/OmCe0SYvPUlN/SIiErUilvidcx8CG2sU/xi43TlX6Y9T7pePAp52zlU65xYDC4Ah\nZlYIZDnnJjjvebePAqMjFXNjyM9MVlO/iIhErYQmXt4RwLFmdhuwC/iFc+4zoAiYEDLeCr9sj/++\nZnmtzOxy4HKAgoICysrKGi3wioqKsOZXsX43Gyr28N777xMX3WclDhBu/ZqzWK9jrNcPYr+OsV4/\niP06Nof6NXXiTwBaAUOBwcCzZtalsWbunHsAeABg0KBBbsSIEY01a8rKyghnfitTl/Lywhn0HDCU\nwuzURlt+pIVbv+Ys1usY6/WD2K9jrNcPYr+OzaF+Td2rfwXwvPNMAqqBfGAlUBwyXnu/bKX/vmZ5\n1Cryb+KzcpPu3iciItGnqRP/i8BIADM7AkgC1gMvA+eaWbKZdcbrxDfJObca2GpmQ/3e/BcCLzVx\nzAelfa6f+Dcr8YuISPSJWFO/mT0FjADyzWwFcDPwMPCwf4nfbmCM32lvppk9C8wCqoArnXP7HnH3\nE7wrBFKBcf4rarXzj/hX6IhfRESiUMQSv3PuvDoGnV/H+LcBt9VSPhno3YihRVRaUgK5aYms2rwT\n5xxjP13C6b0Laes/wEdERCRIunNfBLTOTGZDxW6mLtvE71+ZxYUPTww6JBEREUCJPyLy0r279701\ncy0A89ZWMPpfn/B/L0xnd1V1wNGJiEhLpsQfAa0yvLv3TVi0gR5tM8lMSeCL5Zt5YuIy3p29Nujw\nRESkBVPij4D89CRWbt7JzFVbGdG9DW//7Hjuv2AgifHGtBVbgg5PRERasKa+gU+LkJeRTKXfpD+g\nQw5ts1Nom92WnoVZfLF8U8DRiYhIS6Yj/gholZ60//2Ajrn73/cvzmH6ii3srXZBhCUiIqLEHwk9\nC7P2v8/PSN7/vl/7HLbv3svCdRVBhCUiIqLEHwkDO+Zy5ciu/PbM0gPK9x39vzJtVRBhiYiIKPFH\nyi9P7cEPh3c+oKxzfjrf6FPIQx8vpmpvNdNXbGG6OvuJiEgTUuJvYqf0KmDH7r3MXbuNs+7+mLPu\n/hjvrsUiIiKRp8TfxPq1zwHgov98tr/sD6/OCiocERFpYZT4m1jHvDRO7VXAzt17+UafQgAe+XQJ\nW3ftCTgyERFpCXQdfxMzM+6/YND+zxcu2sD3HpjApws2cFrvtgFGJiIiLYGO+AM2oGMumckJfDCv\nPOhQRESkBVDiD1hifBzHdMvng7nr1MlPREQiTok/CgzrlseqLbtYuXln0KGIiEiMU+KPAv2LvZ7+\n05brmn4REYksJf4o0KNtFkkJcXy5YnPQoYiISIxT4o8CSQlxFOWkskJN/SIiEmFK/FGibVYKa7bs\nCjoMERGJcUr8UaJtthK/iIhEnhJ/lCjISqF82y6qq3VJn4iIRI4Sf5QozE5hz17Hhu27gw5FRERi\nmBJ/lGiTmQzAum2V+8ue/Ww5H81fF1RIIiISg3Sv/iiR7yf+Ddu9xH9v2UL+/MYcAE7v3Za7zjuS\nxHjtp4mIyOFRJokSeelJAKyvqKSisoq/vDmHfsU5nD2giHEz1vDkxGUBRygiIrFAiT9K5GX4R/wV\nu5m9eivOwbUnduOOc/rTrziHJyYuDThCERGJBUr8USIrJYGk+DjWV+xmxkrv1r2922UDcGqvAuat\nrWDLjj1BhigiIjFAiT9KmBl5GUmsr6hkfnkFuWmJtMlKAaBPkbcDMHOV7uUvIiKHR4k/iuRnJLNu\nWyXLNuygQ176/vJe/pH/58t1L38RETk8SvxRpDA7hdVbdrJs4w46tErbX94qPYk+Rdm8NWttgNGJ\niEgsUOKPIu1yUlm+cScrN++kQ6vUA4ad0aeQacs3s2LTjoCiExGRWKDEH0Xa5aSwc89e9lY7OoY0\n9QN8o08hAOOmrwkiNBERiRFK/FGkXc5XR/klbTIOGNYhL43eRVm8Nn11U4clIiIxRIk/inTO/+oo\nv6Qg82vDz+hTyBfLN/O7l2cecGtfERGRcEUs8ZvZw2ZWbmYzahl2nZk5M8sPKbvRzBaY2VwzOzWk\nfKCZTfeH3WVmFqmYg1ZamLX/fUby1++mfO7gDrTPTWXsp0s4+95P2FCh5C8iIgcnkkf8Y4HTahaa\nWTFwCrAspKwUOBfo5U9zj5nF+4PvBS4DSvzX1+YZK8yMsl+M4JWrhtc6vFV6Eq9dcyz/uWgwa7dW\n8puXvrZPJSIiUq+IJX7n3IfAxloG/R34FRD64PlRwNPOuUrn3GJgATDEzAqBLOfcBOecAx4FRkcq\n5mjQKT+dPu2z6xyenZrIyB5tuHJEN16fvoY5a7Y2YXQiItLcNek5fjMbBax0zk2rMagIWB7yeYVf\nVuS/r1ne4l14dEeSEuJ49rMVDY8sIiLia7LH8ppZGnATXjN/pJZxOXA5QEFBAWVlZY0274qKikad\nX2MoyYa3pi3huMzyw55XNNavscV6HWO9fhD7dYz1+kHs17E51K/JEj/QFegMTPP757UHpprZEGAl\nUBwybnu/bKX/vmZ5rZxzDwAPAAwaNMiNGDGi0YIvKyujMefXGGa6Bfz1zbn0GzyMXP+xvocqGuvX\n2GK9jrFeP4j9OsZ6/SD269gc6tdkTf3OuenOuTbOuU7OuU54zfYDnHNrgJeBc80s2cw643Xim+Sc\nWw1sNbOhfm/+C4GXmirmaNejrXfJ3+IN2wOOREREmotIXs73FDAe6G5mK8zskrrGdc7NBJ4FZgFv\nAFc65/b6g38CPIjX4W8hMC5SMTc3xf79/Jdv1G18RUQkPBFr6nfOndfA8E41Pt8G3FbLeJOB3o0a\nXIxon+vd6W/Fpp0BRyIiIs2F7tzXjKUlJZCXnqTELyIiYVPib+Y65qWxcF1F0GGIiEgzocTfzJW2\ny2L26q149zcSERGpnxJ/M1damM22XVVq7hcRkbAo8Tdzvdp5D/aZuWoLAM45ZqzcohYAERGplRJ/\nM9e9bSZxBrNWeffsf2rScs7858fc+e78gCMTEZFopMTfzKUkxtO1dQazVnuJ/4mJSwH4xzvzmbho\ng478RUTkAEr8MaBnYRazV29jfUUlM1dtZVjXPAC+98AEzrjrY9ZXVAYcoYiIRAsl/hjQvW0mKzfv\n5N3ZawH45andefiiQdwyujeL1lVwzn3j2bJjT8BRiohINGjKh/RIhJS0yQDgvg8WkZOWSJ+ibBLi\nvX26znnpnP/QRG5/Yzan9S5kQXkFp/VuS1FOapAhi4hIQJT4Y0DPQq9n/+L12/n2gPb7kz7A8JJ8\nvjOwPU9NWs5Tk5YD8P6cch6/9KhAYhURkWAp8ceA4lZp/PzkI3hn9lquPbHka8NvHd2bvu2z6ZiX\nzuQlG/nnewt45rNlfG9whwCiFRGRICnxx4hrTizhmlqSPng9/y88uhMAw7vlM37hBu56dwHfGVhM\nfJw1YZQiIhI0de5rYeLjjIuP6czKzTv5aP66oMMREZEmpsTfAp1U2oZW6Uk8O3l50KGIiEgTU+Jv\ngZIT4hnRvTVTl24OOhQREWliSvwt1BEFmazZuoutu3R9v4hIS6LE30Ltu/Z/QXlFwJGIiEhTUuJv\noUraZAKwYK0Sv4hIS6LE30IV5aaSkhjH/PJtQYciIiJNSIm/hYqPM7rkZzBfTf0iIi2KEn8LVlKQ\nwXw19YuItChK/C1Y+9xU1mzdRXW1CzoUERFpIkr8LVjrjGT2Vjs27dgddCgiItJElPhbsNaZKQCU\nb6sMOBIREWkqSvwtWJusZADWKfGLiLQYSvwtWOuMrxL/5DVVvPTFSpzT+X4RkVimx/K2YK0zvcT/\n5sw1vDWrEr74gqzUREZ2bxNwZCIiEik64m/B0pMT6JKfzluz1hJvXtnbs9YGG5SIiESUjvhbuGO6\n5bNo/XaGFyWQnJXHB3PX4ZzDzIIOTUREIkCJv4W76oRudMpPp7hyCesyWvPWrLUsXFdBN/9e/iIi\nElsOqqnfzOLMLCtSwUjTK8hK4ZLhnUmKN0b45/bPuX8Ct702i2nLNwccnYiINLYGE7+ZPWlmWWaW\nDswAZpnZLyMfmjS1opxUrj2xhB27q/j3R4u57r/T1MtfRCTGhHPEX+qc2wqMBsYBnYELIhqVBOZn\nJx/BnFtO59bRvVlQXsHMVVuDDklERBpROIk/0cwS8RL/y865PYAOA2Pcqb3aAvDxgvUBRyIiIo0p\nnMR/P7AESAc+NLOOQIOHgWb2sJmVm9mMkLK/mtkcM/vSzF4ws5yQYTea2QIzm2tmp4aUDzSz6f6w\nu0zdzZtE68xkurXJYPzCDUGHIiIijajBxO+cu8s5V+ScO8N5lgIjw5j3WOC0GmVvA72dc32BecCN\nAGZWCpwL9PKnucfM4v1p7gUuA0r8V815SoQM7pTLF8s36zy/iEgMCadzX4GZPWRm4/zPpcCYhqZz\nzn0IbKxR9pZzrsr/OAFo778fBTztnKt0zi0GFgBDzKwQyHLOTXBe9nkU75SDNIHeRdls2bmHFZt2\nBh2KiIg0knCa+scCbwLt/M/zgJ82wrJ/iNdZEKAIWB4ybIVfVuS/r1kuTaB3u2wAJi/d2MCYIiLS\nXIRzA59859yzZnYjgHOuysz2Hs5Czez/gCrgicOZTy3zvRy4HKCgoICysrJGm3dFRUWjzi/a1Fa/\nqmpH2zTjZ89M4+YXpnFl/xS6t4qvfQbNQEv8DmNNrNcx1usHsV/H5lC/cBL/djPLw+/Jb2ZDgS2H\nukAzuwg4EzjRfXXyeCVQHDJae79sJV+dDggtr5Vz7gHgAYBBgwa5ESNGHGqYX1NWVkZjzi/a1FW/\nvG6b+N79E9i6u5q7p1Xx/E+GUFLQPO/q11K/w1gS63WM9fpB7NexOdQvnKb+nwMvA13N7BO88+xX\nH8rCzOw04FfAN51zO0IGvQyca2bJZtYZrxPfJOfcamCrmQ31e/NfCLx0KMuWQ3Nkh1w+ueEEPvrV\nSJIT47nx+elBhyQiIoehwSN+59xUMzse6A4YMNe/lr9eZvYUMALIN7MVwM14vfiTgbf9q/ImOOd+\n5JybaWbPArPwTgFc6ZzbdzrhJ3j9DFLx+gSMQ5rUvsf3njekmHvKFrJt1x4yUxIDjkpERA5FnYnf\nzM6uY9ARZoZz7vn6ZuycO6+W4ofqGf824LZayicDvetbljSNo7vk8c/3FjB5ySZG9mgTdDgiInII\n6jviP6ueYQ6oN/FL7OnT3uvlP2v1ViV+EZFmqs7E75y7uCkDkeiXmZJIUU4q89ZuCzoUERE5ROHc\nwOePNW6tm2tmt0Y2LIlWJQUZzFtbEXQYIiJyiMLp1X+6c27/g9mdc5uAMyIXkkSz7gWZLCyvoGpv\nddChiIjIIQgn8cebWfK+D2aWitczX1qgkoJMdu+tZunGHQ2PLCIiUSecG/g8AbxrZv/xP18MPBK5\nkCSadfdv3jN/7Ta6ts4IOBoRETlY4Tyd7894l9n19F+3OOf+EunAJDp1a5NBfJwxY6X3ZOblG3dw\n66uzWK4WABGRZiGcI36cc7pxjgCQmhRP76JsJizaAMDvXp7Ju3PKSUtO4OcnHxFwdCIi0pA6j/jN\n7GP/7zYz2xry2mZmW5suRIk2R3fJ44vlm9mycw/TVnj9Pmev1iYhItIc1Jn4nXPD/b+ZzrmskFem\ncy6r6UKUaHNa77ZUVTsuf3Qy6yt2AzBnjRK/iEhzEM51/I+FUyYtR7/22XTKS2Pi4o10zk/nsmM7\ns3zjTtZu3RV0aCIi0oBwLufrFfrBzBKAgZEJR5oDM+POc49kRPfWPHXZUL51pPfk5A/mrgs4MhER\naUh95/gWnGzpAAAgAElEQVRvNLNtQN/Q8/vAWvRo3BavX3EOYy8eQtvsFHoWZlLcKpXHJixl0/bd\nQYcmIiL1qO8c/5+AbODRGuf385xzNzZdiBLtzIzrTu7O9JVbGHb7e8xapfP9IiLRqt6mfudcNTC4\niWKRZmz0kUW8fs2xJMQZ/ypbEHQ4IiJSh3DO8U81MyV/aVBpuyx+MLQj46avZumG7UGHIyIitQgn\n8R8FjDezhWb2pZlNN7MvIx2YNE8XDetEtYPXp68BYENFJZt36Ly/iEi0COfOfadGPAqJGW2zU+ic\nn86UpRt5fXoa1z79OW0yU3jxymNonalnO4mIBC2ce/Uvdc4tBXYCLuQlUqvBnXJ5Z3Y5P3liKmlJ\nCazZuournpwadFgiIkJ4N/D5ppnNBxYDHwBL0H37pR4/GdGNNpnJ5Gck8erVw7npjJ5MXLyRGSu3\nBB2aiEiLF845/luAocA851xn4ERgQkSjkmatU346n9xwAp/ecCLFrdL4zoD2JCfE8cxny4MOTUSk\nxQsn8e9xzm0A4swszjn3PjAownFJM5cYH0dSgrd5ZaclckafQl78YiV79lYHHJmISMsWTuLfbGYZ\nwIfAE2Z2J6BrteSgnNCjDdt2VTF3zbagQxERadHCSfyjgB3Az4A3gIXAWZEMSmJP/+IcAL5Yvjng\nSEREWrZ6E7+ZjQZ+DJzsnKtyzj3inLvLb/oXCVv73FTyM5L5dOF6AJZu2M4r01axoFwtACIiTanO\n6/jN7B68J/N9CtxiZkOcc7c0WWQSU8yMM/sW8sTEpfzjnXn8870F7K32rgr97Zml/HB454AjFBFp\nGeo74j8OOMF/IM8IYHSTRCQx66JhnchJS+If78wnJSGO+84fwDHd8vjDq7N4cuKyoMMTEWkR6rtz\n327n3F4A59wOM7MmikliVKf8dD745Qhmr95KamICpe2y6N42i+/c+yn/9+J0+rbPpndRdtBhiojE\ntPqO+Hv49+b/0symh3zWvfrlkKUlJTCwYytK22UB0Dk/nfeuG0GcGW/MWBNwdCIisa++I/6eTRaF\ntGjZaYn0Kcpm/CL1GRURibQ6E79/f36RJjGkcyvGfrKEqr3VJMSHc5WpiIgcCv3CSlQ4oiCT3Xur\nWbpxBwA7dlfxr/cXMHv11oAjExGJLUr8EhW6F2QCMM+/s9/1/5vOX9+cyzn3jWfn7r1BhiYiElPq\nTPxm9q7/989NF460VN3aZBBnMHPVVhauq+CVaavoXZTFtsqq/Tf9ERGRw1ffEX+hmQ0DvmlmR5rZ\ngNBXUwUoLUNqUjx9irKZuHgD46avBuDeHwwkPSmed2aXBxydiEjsqK9X/2+B3wDtgTtqDHPACZEK\nSlqmoV3zeOijxSzfuJMjO+RQ3CqN445ozXtz1uJcb3QrCRGRw1fnEb9z7jnn3OnAX5xzI2u8Gkz6\nZvawmZWb2YyQslZm9raZzff/5oYMu9HMFpjZXDM7NaR8oH/vgAVmdpduJBS7vtmvHVXVjjVbd3He\n4A4AnNizgLVbK5m+ckvA0YmIxIYGO/c5524xs2+a2f/zX2eGOe+xwGk1ym4A3nXOlQDv+p8xs1Lg\nXLxnA5wG3GNm8f409wKXASX+q+Y8JUb0apfND4/pzJijOzLqyHYAnNijDfFxxpszdXMfEZHGUF9T\nPwBm9idgCPCEX3StmQ1zzt1U33TOuQ/NrFON4lF49/0HeAQoA673y592zlUCi81sATDEzJYAWc65\nCX4sj+I9M2BcQ3FL8/Tbs0oP+JybnkTf9tl8tmRTQBGJiMQWc87VP4J3e97+zrlq/3M88Llzrm+D\nM/cS/6vOud7+583OuRz/vQGbnHM5ZnY3MME597g/7CG85L4EuN05d5JffixwvXOu1lYHM7scuByg\noKBg4NNPP91QiGGrqKggIyOj0eYXbaK5fg9Nr+TL9Xu5c2TaYc0nmuvYGGK9fhD7dYz1+kHs1zHI\n+o0cOXKKc25QQ+M1eMTvywE2+u8b5SkqzjlnZvXvdRz8PB8AHgAYNGiQGzFiRKPNu6ysjMacX7SJ\n5vrNZiEfvTGHgUOPITMlkXXbKslITiA1Kb7hiUNEcx0bQ6zXD2K/jrFeP4j9OjaH+oWT+P8EfG5m\n7wOG97jeGw5xeWvNrNA5t9rMCoF912mtBIpDxmvvl63039cslxakc346AIvXb6ddTiqn/P0DStpk\n8swVQ9XTX0TkIIXTue8pYCjwPPA/4Gjn3DOHuLyXgTH++zHASyHl55pZspl1xuvEN8k5txrYamZD\n/VMDF4ZMIy1Ej7beXf1mrtrKra/OYtOOPUxaspFxepqfiMhBC6up30/ALx/MjM3sKbyOfPlmtgK4\nGbgdeNbMLgGWAuf4859pZs8Cs4Aq4Ern3L77tP4E7wqBVLzz/urY18J0zEsjOzWRN2as4YN567ji\nuC58MG8dN70wnbbZKQzokNvwTEREBAj/HP9Bc86dV8egE+sY/zbgtlrKJwO9GzE0aWbMjCM75FA2\ndx0A3x1UzCm9Cvj2veM59/4JjPvpsXRtHbudhUREGpMe0iPNwo+O74oZjOzemm5tMhjYsRXvXnc8\nAE9MWBZwdCIizUe9R/z+pXsznXM9migekVoN7ZLHJ9efQEFWyv6yrq0zGNK5FR8vWBdgZCIizUu9\nR/z+efa5ZtahieIRqVO7nFTi4w7sxX9Mt3zmra2gfOuugKISEWlewmnqzwVmmtm7ZvbyvlekAxMJ\nx7El+QB8vECP7hURCUc4nft+E/EoRA5RaWEWuWmJfLxgPWcPaM+aLbtIiDfyM5KDDk1EJCo1mPid\ncx+YWUegxDn3jpmlAQd3yzSRCImLM4Z1y+fj+etZUF7BSXd8QGZKApNuOumg7+wnItISNNjUb2aX\nAc8B9/tFRcCLkQxK5GAc2y2f8m2V/PK5aQBs21XFuBmrA45KRCQ6hXOO/0rgGGArgHNuPtAmkkGJ\nHIzh/nn+z5dtZlT/dhTlpPLzZ6dxx9vzAo5MRCT6hJP4K51zu/d9MLMEoFEfriNyONrnptEu27vM\nb/SRRbx2zXCGd8vn6UnLaOjpkyIiLU04if8DM7sJSDWzk4H/Aq9ENiyRg/PgmME8eelRjOzehpy0\nJEb1b0f5tkpmr94WdGgiIlElnMR/A7AOmA5cAbwO/DqSQYkcrNJ2WQzrlr//89Fd8wD4bMnGuiYR\nEWmRwunVX21mjwAT8Zr45zq1n0qUK8pJpW1WCpOXbmLMsE5BhyMiEjUaTPxm9g3gPmAhYEBnM7vC\nOaen5EnUMjMGdcplin/E75xj227tr4qIhNPU/zdgpHNuhHPueGAk8PfIhiVy+AZ1zGXVll2s3LyT\n656dxtXv7WD6ii1BhyUiEqhwEv8259yCkM+LAPWYkqg3qFMrAB4dv4TnP18JwH+nLA8wIhGR4NWZ\n+M3sbDM7G5hsZq+b2UVmNgavR/9nTRahyCEqLcyiS+t07v9gESmJcfRsFccr01axu6o66NBERAJT\n3xH/Wf4rBVgLHA+MwOvhnxrxyEQOU1yc8etv9KQoJ5WbzujJqZ0S2bRjD397ay4TFm0IOjwRkUDU\n2bnPOXdxUwYiEgkn9CjghBsKAHhv12JO6tmG+z9cxP0fLuKjX42kuFVawBGKiDStcO7V39nM7jCz\n5/VYXmnO4sy45wcDueL4LgC86J/3FxFpScJ5LO+LwEN45/Z1clSataSEOG48vSefL9vMC5+v5KoT\numFmQYclItJkwunVv8s5d5dz7n3n3Af7XhGPTCSCzj6yiEXrtzPNv7yvorKKD+at0739RSTmhZP4\n7zSzm83saDMbsO8V8chEIuj0PoUkJcTxwtQVVO2t5vJHJzPm4Ulc/dTnLN2wPejwREQiJpym/j7A\nBcAJfNXU7/zPIs1SdmoiJ/cs4KnPljN7zTYmLfbu8Pfql6uZvXor7/z8eJ0CEJGYFM4R/3eBLs65\n451zI/2Xkr40e98Z1J7dVdVMWryR28/uw5Lbv8Eto3uzcN12Zq7aGnR4IiIREU7inwHkRDoQkaY2\nsnsb7jt/IA+NGcS5QzoA8I0+hcTHGa9+uTrg6EREIiOcpv4cYI6ZfQZU7it0zn0zYlGJNJHTerc9\n4HOr9CSGdc1j3IzV3HB6j4CiEhGJnHAS/80Rj0IkihxX0prbXp/Num2VtM5MDjocEZFG1WDi16V7\n0tL0K/bObH25YjNFualc9PBnPDhmEL2LsgOOTETk8DWY+M1sG14vfoAkIBHY7pzLimRgIkHpXZRF\nYrwxcfFGZny8hTVbd1E2t1yJX0RiQjhH/Jn73pt3fdMoYGgkgxIJUlpSAkd3zeeBDxcBEGcwddnm\ngKMSEWkc4fTq3895XgROjVA8IlHhwqEdiY8zju6Sx+j+RXy+bJPu6iciMSGcpv6zQz7GAYOAXRGL\nSCQKnFRawJRfn0RaUgL/m7qC5z9fyeL12+nSOiPo0EREDks4vfrPCnlfBSzBa+4XiWk5aUkADOiQ\nC8DExRuV+EWk2QvnHP/FTRGISLQqaZPBEQUZ3Pj8dNKS4hnVvyjokEREDlmdid/MflvPdM45d0sE\n4hGJOnFxxi2jevPzZ6fx82ensb1yL+cNKda9/EWkWaqvc9/2Wl4AlwDXH85CzexnZjbTzGaY2VNm\nlmJmrczsbTOb7//NDRn/RjNbYGZzzUwdC6XJHdUlj7d+dhzDu+Vz0wvT+f6/J/LhvHV8smC9Ov2J\nSLNS5xG/c+5v+96bWSZwLXAx8DTwt7qma4iZFQHXAKXOuZ1m9ixwLlAKvOucu93MbgBuAK43s1J/\neC+gHfCOmR3hnNt7qDGIHIr05ATuO38gv3huGq99uZrxizYAcN6QDhzZIYdTSgv29wsQEYlW9Z7j\nN7NWwM+BHwCPAAOcc5saabmpZrYHSANWATcCI/zhjwBleC0Lo4CnnXOVwGIzWwAMAcY3QhwiByU1\nKZ5/fK8/PQoyKSnI4KP563li4jKemrSMhwoyee2a4STEH9RVsiIiTaq+c/x/Bc4GHgD6OOcqGmOB\nzrmVZvb/gGXATuAt59xbZlbgnNv3SLQ1QIH/vgiYEDKLFX6ZSCAS4+O4+sQSAE4pbcuxJa0Zv3A9\nj4xfysvTVnH2gPYBRygiUjer6/ykmVXjPY2viq9u2QtgeJ37DumWvf65+/8B3wM2A/8FngPuds7l\nhIy3yTmXa2Z3AxOcc4/75Q8B45xzz9Uy78uBywEKCgoGPv3004cSYq0qKirIyIjdS7livX4Q2To6\n57jho53kphg3DEll465qNu5ydM2Oa7JOgPoOm79Yrx/Efh2DrN/IkSOnOOcGNTRefef4I9VeeRKw\n2Dm3DsDMngeGAWvNrNA5t9rMCoFyf/yVQHHI9O39stpifgCvhYJBgwa5ESNGNFrQZWVlNOb8ok2s\n1w8iX8cLWcBf35zLe1vyeX7qSioqq7hgaEd+cWp3slMTI7bcffQdNn+xXj+I/To2h/oFcTJyGTDU\nzNL8e/+fCMwGXgbG+OOMAV7y378MnGtmyWbWGSgBJjVxzCINuuzYLhzVuRWPjl9KdmoiJ5cW8NiE\npQy//T1mrNwSdHgiIkB4d+5rVM65iWb2HDAV7zTC53hH6RnAs2Z2CbAUOMcff6bf83+WP/6V6tEv\n0SgpIY4HLhjEq9NXcUppW1pnJvPlis386LEpXPLIZ7xx7XHkpqvXv4gEK5Dux865m51zPZxzvZ1z\nFzjnKp1zG5xzJzrnSpxzJznnNoaMf5tzrqtzrrtzblwQMYuEIzstkR8c1ZHWmckA9G2fw53nHcna\nrZW8O6e8galFRCJP1x2JRNjADrnkpiVSNleJX0SC1+RN/SItTVyccUppW56ZvJwz+66mICuFP42b\nQ3JCHJcM78yI7m2CDlFEWhAlfpEm8IfRvZi6bBPXPPUFVdXVVDuIM5iwaANPXDqUIZ1bBR2iiLQQ\nSvwiTSA5IZ6/ndOP/01ZQWZKIj8Y2oG0xAS+de8nXPHYZF688hg65qUHHaaItAA6xy/SRPq2z+H3\no3rzi1O7U5idSnZaIg+PGYwDfvrMF+yt1sN+RCTylPhFAtQpP52bzyrl82WbeXzCUrbs3KOn/YlI\nRCnxiwRsdP8ijumWx80vz6Tf79/i1tdmK/mLSMQo8YsEzMy44bSe+2/r+9DHi/nbW/MCjkpEYpU6\n94lEgT7ts5l28yk457ju2WncU7aAuWu3MbhTLpcf1zXo8EQkhijxi0QRM+OW0b3B4PmpK3l71lrS\nkhI4f2jHoEMTkRihpn6RKJOenMAd5/Rn4R/PYET31tzy6izmrd2m8/4i0ih0xC8SpeLjjFtG9ebE\nOz7glL9/SPvc1P3PAOjYKo3ju7dmdP8ivIdcioiER4lfJIoVt0rjlauG8/7ccj6Yu451FZU453h3\nTjkvfrGKVunJHH9E66DDFJFmRIlfJMp1b5tJ97aZ/Oj4rzr57a6qZsRf3+fBjxYp8YvIQdE5fpFm\nKCkhjm/2L2L8wg1s2bEn6HBEpBlR4hdpps7qV0hVteOMuz7itH98yPhVVUGHJCLNgJr6RZqpXu2y\n+fGIrtz/wUI2bI9j3tpq2o9fQlZqImf1bUdcnDr9icjXKfGLNGPXn9aD604+go3bd3PM7e/ym5dm\nArB0ww6uObEk4OhEJBqpqV+kmUuIj6NNVgo/H5jCX77dlxN6tOFf7y9g6YbtAMxYuYUP560LOEoR\niRZK/CIxomdePOcMLubms0pxDr51z6eMm76aM//5MRc+PIlXv1wVdIgiEgWU+EViTMe8dMb+cDAb\nt+/mx09MJS89ifSkeK568nO+d/94rnxyKusrKoMOU0QConP8IjFoWNd8/vLtvrw/t5w/jOpNVmoC\nN780k4XrKnhzxhqyUhL509l9gg5TRAKgxC8So84ZXMw5g4v3f779230B+N3LM3lswlK+P6QDG7ZX\ncmSH3P2PBBaR2KfEL9LCXH1CN/47eTln3f0xAF3y07nje/3pX5wTcGQi0hR0jl+khcnLSObXZ5Yy\npHMrrhzZlU07dvPtez/l5pdmUF2tJwCKxDod8Yu0QOcN6cB5QzoAcMXxXbnt1dk8Mn4ppe2y+N7g\nDgFHJyKRpCN+kRYuKyWR27/dhw6t0nh71tr95RWVVezaszfAyEQkEnTELyKYGSO6t+bZyct5atIy\nPl6wnnHTV9M6M5n3rhtBerJ+KkRihY74RQTwmvyzUhK58fnpvPblajrlpbN2ayV3vjs/6NBEpBFp\nN15EACjKSeXd647nsyUbWbphB98bXMwtr87mgQ8XMWPlFk4pLeDCozvp4T8izZwSv4jsl5mSyAk9\nCvZ/vnV0b4pyUnhq0nJ+98osFq3fzp69jvKtu+jbPodzhxTTJjMZM+0MiDQXSvwiUqf4OOOqE0q4\ncmQ3bn55Jo+OX7p/2Ltzyvn7O/PoXpDJrd/qzeBOrfYPK9+2i8zkRFKT4oMIW0TqocQvIg0yM37/\nzV4cW9KaxHhjRPc2TF22iclLNvLo+KVc89TnvH7NseSmJ/Hwx4v5w6uzOLNvIXd/f0DQoYtIDUr8\nIhIWM+Pk0q9OAwzokMuADrkM7ZLHd+4bz80vz+TKkd249bVZALz65Wp27fmMozrncemxnXU6QCRK\nKPGLyGHp2z6H7w/pwKPjlzBl6SYyUxJ57kdHc+3TX/DunHLemV3O4g3bKS3MYsrSTeSkJfKrU3vo\nNIBIQJT4ReSw/eCoDjw3ZQUrN+/kD6N6UVKQyevXHkt1teMPr85i7KdLDhj/hc9X8tyPhtGtTUYw\nAYu0YIEkfjPLAR4EegMO+CEwF3gG6AQsAc5xzm3yx78RuATYC1zjnHuz6aMWkbqUFGQy/sYT2Lh9\nNx3z0veXx8UZN59VSvvcVLq2zqB/cQ4fzl/HzS/P5Ox7PuG5Hw/jiILMACMXaXmCuoHPncAbzrke\nQD9gNnAD8K5zrgR41/+MmZUC5wK9gNOAe8xMbYQiUSYzJfGApL+PmXHpsV0Y2aMNuelJjOpfxMtX\nDic5MZ6rn/yc3VXVVO2t5pVpq7js0ckc95f3uevd+TinBwaJREKTH/GbWTZwHHARgHNuN7DbzEYB\nI/zRHgHKgOuBUcDTzrlKYLGZLQCGAOObNHARaTQd8tL407f6cOmjk7nvg4VUO8c/3plPUkIcpYVZ\n3PH2PP7+zjyOLWnNmX0KOb1PW1Zs2snGXdVBhy7S7FlT71WbWX/gAWAW3tH+FOBaYKVzLscfx4BN\nzrkcM7sbmOCce9wf9hAwzjn3XC3zvhy4HKCgoGDg008/3WhxV1RUkJERu+cjY71+EPt1bI71u2/a\nLias9h4EVJoXx08HpJAQBx+uqGLRlmomr6liR9VX4yeY44p+KQxu+/Vjll1Vjm27Ha3Tmu+dyJvj\nd3iwYr2OQdZv5MiRU5xzgxoaL4hz/AnAAOBq59xEM7sTv1l/H+ecM7OD3iNxzj2At1PBoEGD3IgR\nIxohXE9ZWRmNOb9oE+v1g9ivY3Os35BhVfzo8anMW7ONOy4YTI+2WQCc4A93zjFuxho+X7aJlMR4\nXpm8iEdn72XMN4bTJjNl/3x+9/LM/R0Irzv5CK4+saSJa9I4muN3eLBivY7NoX5BJP4VwArn3ET/\n83N4iX+tmRU651abWSFQ7g9fCRSHTN/eLxORZi4tKYFHfziEvdWO+FqeAWBmnNGnkDP6FAJQtGcl\nv/l0F397cx5//k5fAD6Yt46xny5hWNc8AO58dz6DOrXiaP+ziByoydvEnHNrgOVm1t0vOhGv2f9l\nYIxfNgZ4yX//MnCumSWbWWegBJjUhCGLSITVlvRr0zY9jouGdeLZKcuZsXILO3ZX8X8vTKdL63Qe\nvmgw954/kI55afzo8Sl8sXxzhKMWaZ6COhl2NfCEmX0J9Af+CNwOnGxm84GT/M8452YCz+LtHLwB\nXOmc2xtI1CISuKtOKCEvPZkxD0/i4v98xopNO/nTt/qQkhhPdmoi/7loCJkpCZz/4ERe+3I1yzfu\nCDpkkagSSOJ3zn3hnBvknOvrnBvtnNvknNvgnDvROVfinDvJObcxZPzbnHNdnXPdnXPjgohZRKJD\ndmoif/lOHzbu2M20FZv53VmlHNXlq2b9DnlpPHvF0eSkJXLlk1M5/q/v88aMNV+bT9Xeavbs1VUC\n0vLozn0i0uyc0KOAt392HG2yUshKSfza8HY5qYy79limLd/CX9+cw4+fmMI/zzuSb/QpxMxwznHZ\no5OZuWor15/Wg7MHFOlZAtJiKPGLSLPUrU39d/zLTElkeEk+AzsezQ8enMBVT37OHa3ncdGwTkxY\ntIH3564D4Lr/TuNf7y9gcKdWXHRMJ3oWZjVF+CKBab4XvIqIhCE1KZ7HLz2KG0/vQcWuKn770kzG\nzVjDL045gkV/PINbR/emuFUar3y5im/e/TG/f2Umu/aoG5HELh3xi0jMS0tK4Irju/LdQcXMXr2V\nkoKM/fcBOH9oR84f2pGN23fzp9dn859PlvDZko3cOroP7bJTaJOV0sDcRZoXJX4RaTFapSdxTLf8\nOof99bv9OLVXW3727BeM/tcnAIzu347rT+9BYXZqU4YqEjFK/CIiIU4qLeD9X4zg4/nrmbpsE09/\ntpxPFm7ghZ8Mo31uGkvWb2fH7r2UtlNfAGmelPhFRGrIz0hm9JFFjD6yiO8f1YHv3jeeS8ZO5oKj\nO/L7V2ayZ6+jpE0Gvzy1O13bZNApLz3smxCJBE2JX0SkHj3aZvGv7w/gwocn8esXZ5DnP1r4jRmr\nufyxKQCc1LMNv/5GKcmJcUxbvoWF6yoY0rkVgzrm6jJBiTpK/CIiDTjuiNb89sxS9uyt5pxBxeSm\nJ3HD6T14b045788p55nJy3lndvnXpjultIDfj+pF+dZKZq/eynFHtKZdjvoKSLCU+EVEwvDD4Z0P\n+JyUEMdpvdtyWu+2fGtAEas272RDxW7aZCXTr30Ot742i7dmreWtWWsPmK5H20z+feEgilulNWX4\nIvsp8YuIHKahXb7+JMB/XziIKUs38fmyzbTOTGbP3mo+mr+e9+eU8617PqFv+xza5aRw9Qkl7K12\nzFi5hYkr9pC0cD0lbTJpnZkcQE2kJVDiFxGJADNjUKdWDOrUan/ZdwcV8+WKzfzx9dms3rKLsrnl\nPD5h2QHTPTRjIknxcfx7zCCOP6I16ysq+WTBejrlpdOvOKepqyExSIlfRKQJ9W2fw9OXHw3A/LXb\neG9OOfFxRu+ibOZM/4Jupf249bVZXPXEVP75/SP51XNfUr6tEoDzhnSgKCeFQZ1a1drKIBIOJX4R\nkYCUFGRSUvDVMwd2LYtneEk+D100mFF3f8JF//mMhDjjoTGDePGLVTw16avWgc756SQnxLFi007y\nMpLo3S6bIzvkcMHRHUlOiA+iOtJMKPGLiESZopxU/vujo3ln1lqGdcujV7tsTuxZwB3n9GN3VTUP\nfrSYmau2ANC9bSbbK6uYvHQjr01fTdncdfzp7D6Mm7GaODMuPLoTSQl6LIt8RYlfRCQKdc5P57Lj\nuhxQlhgfR2J8HNeeVFLrNI98uoTbx83h2L+8v7/s7vcXULXXUdwqjV+d1p1ju+WzbOMOMlIS9j+v\nQFoWJX4RkRgxZlgnTujRhoc+Xky7nBRK2mTyyper9l81cNkjk2mTmfz/27vz6CrrO4/j729WskFI\niJCQAAGiAkpYwioquKCCRc5oKT3UsUcsQ6VVj8fTkXG0tUePjh0dxy5TLdpji4iKWtFaFS2ojFBl\nEWTfV9mMgBCWEPKdP+5DuFAWJ433Js/9vM655z7P73ly+X7I8n2We5+Hz/ceAuCG3sXcO7wrLTJT\n41y5xJIav4hIiJTkZfKzEd3q5oecfw4Au6uqeeDPy9n0ZRVjL+7I9r0HmTR7PdPmb2FEeRGPjSon\nJVmnBBKBGr+ISAJomZXGo6PKTxgb2q0N0z/9nD/O3ciqHfu4eVAp0+Zv4UB1Ddf3KmZURQlZ6WoT\nYeLpQpoAAA22SURBVKPvqIhIgurTIY8+HfLoVtScn05fyk+mLSYzLZn2+Vnc//oynnx/HaWtsujQ\nKpM7rzxPFxUKCTV+EZEEN7pvO/qW5vG7D9czpl87Lmjbgo/WfMGjM1YxZ10lc9ZVsrHyAD3b5bL5\ny4OMqihhUFmreJct9aTGLyIidCzI5qF/urBufmDnVgzs3Ap356G/rOCpD9bx0dpK0lKSeH3x5wy7\nsJCc9BRGlBcxsLM2ApoSNX4RETktM+PfhnVh4jXnA3DoSC0/f2MZM1fsZO/BI7w0fwuv/HAgRbkZ\nrNy+j17tc8lMU2tpzPTdERGRszIzADLSkuuODOw9cIShj7/Pt5+cQ22tU1PrNG+WQrv8THq1a8mt\ngzvTpoWuFdDYqPGLiEi9tMhMZcoP+vObmWspyEmn7Jxs3l+1i217D/KHORv5aG0lk/65gsLcZiz7\n/Cs6FmTj7rh73YaExJ4av4iI1FunguwTPiZ4fe9iAOasreSHz83nW7+aTWpyEl9WVdetc/6iDxna\nrQ1JBuUluaSnJNEhP4ui3IyY15+I1PhFRKTBDeiUz+s/GsRPpy8lPSWJi8sK+GzrXmYv38KufYd5\n4r3Vf/c15+Sk8+sxvegTdStjaXhq/CIi8o0oycvkme/3OWFs1qxKBgy6mKrDR9lYWcW2vYf4cPUu\ntuw+yIKNu5nw3AKm/KA/nc/JjlPV4afGLyIiMZWekkx6SjJ5WWn0BIZdWAjA8m1fMWbS37h96kIm\nj+3Hiu37mLuuki6FORTlZtC9ODe+hYeEGr+IiDQKXQqb88j13bn1uQX0efBdamr9hOVXdWvNmH7t\nObd1DvnZaaTq3gL1osYvIiKNxhVdW/PS+AH84u2VlLbK4o4ryti8+yBvL93O8x9v4u2lOwBITTZ6\nlrQkIy2ZrXsOUtIyg4KcdB4YeSFpKdogOBM1fhERaVTKS3KZfEu/uvn87HR6lORy++VlvLNsBxu+\nqGLvwSPMXLGTbXsP0bEgi5krdwEworytLid8Fmr8IiLSJDRLTWZEeVHd/L3Xdq2bPlBdQ4/7Z/DW\n0m1UVh3mkbdWcuuQTozp1z4epTZqavwiItLkZaalMLJnEZPnbmLy3E0A3PPqEhZv3ssFbZtT1jqH\n/h3z41xl46DGLyIiofDz6y6gfX4W1TW13NC7mF/+dTUvzNvMC/Miy++4oowkM9rlZfKt8iKSkxLz\n6oFxa/xmlgzMA7a6+7Vmlge8AHQANgCj3H13sO5EYCxwFLjN3d+OS9EiItJoNUtNZsKQznXzj9xQ\nzj3Du3KguoZ7/7SUx989ftGgB99cTv+O+dx8UQd6tmsZj3LjJp5vfbwdWB41fzfwnruXAe8F85hZ\nV2A00A24GvhNsNEgIiJyRi0yUilskcET3+3B1d3a8MgN3Xn02+WUF7fgrSXb+P7vP6Fy/+F4lxlT\ncWn8ZlYMDAcmRQ1fBzwbTD8LjIwan+ruh919PbAG6BurWkVEpOnLTEvhtzf2ZlRFCdf3LmbSTX14\n87aL2XfoCPdNX0rtSdcMCDNzj31YM5sGPATkAHcFh/r3uHtusNyA3e6ea2a/Aua6++Rg2dPAX9x9\n2iledxwwDqB169a9p06d2mA179+/n+zs8F5CMuz5IPwZw54Pwp8x7Pmg8WWctqqaN9Yd4cc905m3\nvYYkM27smkazlPqd/49nviFDhsx394qzrRfzc/xmdi2w093nm9ngU63j7m5m/+8tEnd/CngKoKKi\nwgcPPuXL18usWbNoyNdrbMKeD8KfMez5IPwZw54PGl/GgYNqef+BGfxy4fHD/eurUqjokEdykjGg\nUz7Xdi/kSI1TU1tLfnb6GV+vseU7lXi8ue8iYISZDQOaAc3NbDKww8wK3X2bmRUCO4P1twIlUV9f\nHIyJiIj8Q9JSkviXSzry2IxVXNWtDaP7tuOZ2etZvGUPB6qP8urCrfzry4txh+Qk48K2LUhLTuLK\nrq255eJSIgeom5aYN353nwhMBAj2+O9y9++Z2S+Am4CHg+fXgi+ZDkwxs8eAIqAM+DjWdYuISDj9\n6LIyxl/aiZTg2v+XnlsAgLvzyYbdzF69CzPjQHUNS7Z+RWXVYR58czlFuRkM717IpsoDLNqyh+HB\nzYYau8b0Of6HgRfNbCywERgF4O5LzexFYBlQA0xw96PxK1NERMIm5RQ3/DEz+pbm0bc074TxmqO1\nDH38AyZMWcCzc/JYsHE3NbXOE++t5vI2NQyOUc31FdfG7+6zgFnBdCVw+WnWexB4MGaFiYiInEZK\nchKTx/ZjwpQFzNvwJUO7tqFVThofrankt4sPcyBrCROv6UJGWuP85Hlj2uMXERFpEopyM5g6rj/7\nD9XUveHv0JGjjHvyXf44dyPrv6jiqRsr2Lz7AC8v2MJ5rXOoaJ9Hu/zMOFeuxi8iIlIv6SnJpGcf\n36tvlprMzRekM7xfF+5+5TNG/vp/+WL/YSqrqoHIrYSnjutP7/Z5p3vJmNBNi0VERBrQd/q045mb\n+rD34BEAXho/gNcmXERBdjr3vbaUeFw/J5r2+EVERBrYkPPPYc7EywDqPvJ311XnceeLi3jgz8v5\n9+Fd4vZRQDV+ERGRb8DJjX1kj7Ys3LSHp2ev58PVu/hOn3aMHVQa87rU+EVERGIgKcm4f0Q32rbM\n4KO1laSnxOdsuxq/iIhIjCQlGeMv7cT4SzvFr4a4/csiIiISc2r8IiIiCUSNX0REJIGo8YuIiCQQ\nNX4REZEEosYvIiKSQNT4RUREEogav4iISAJR4xcREUkgavwiIiIJRI1fREQkgajxi4iIJBA1fhER\nkQRi7h7vGr4RZrYL2NiAL9kK+KIBX6+xCXs+CH/GsOeD8GcMez4If8Z45mvv7gVnWym0jb+hmdk8\nd6+Idx3flLDng/BnDHs+CH/GsOeD8GdsCvl0qF9ERCSBqPGLiIgkEDX+r++peBfwDQt7Pgh/xrDn\ng/BnDHs+CH/GRp9P5/hFREQSiPb4RUREEoga/1mY2dVmttLM1pjZ3fGup77M7Bkz22lmS6LG8sxs\nhpmtDp5bRi2bGGReaWZXxafqr8/MSsxsppktM7OlZnZ7MB6KjGbWzMw+NrNFQb77g/FQ5ItmZslm\nttDM3gjmQ5PRzDaY2Wdm9qmZzQvGQpMPwMxyzWyama0ws+VmNiAsGc3svOB7d+zxlZnd0eTyubse\np3kAycBaoCOQBiwCusa7rnpmuQToBSyJGnsEuDuYvhv4j2C6a5A1HSgN/g+S453hLPkKgV7BdA6w\nKsgRioyAAdnBdCrwN6B/WPKdlPVOYArwRgh/TjcArU4aC02+oO5ngVuC6TQgN2wZg9qTge1A+6aW\nT3v8Z9YXWOPu69y9GpgKXBfnmurF3T8Avjxp+Doiv6QEzyOjxqe6+2F3Xw+sIfJ/0Wi5+zZ3XxBM\n7wOWA20JSUaP2B/MpgYPJyT5jjGzYmA4MClqOFQZTyE0+cysBZGdjKcB3L3a3fcQooxRLgfWuvtG\nmlg+Nf4zawtsjprfEoyFRWt33xZMbwdaB9NNOreZdQB6EtkrDk3G4BD4p8BOYIa7hypf4HHgJ0Bt\n1FiYMjrwrpnNN7NxwViY8pUCu4DfB6drJplZFuHKeMxo4PlguknlU+MXILJHSeSPUpNmZtnAy8Ad\n7v5V9LKmntHdj7p7D6AY6GtmF5y0vEnnM7NrgZ3uPv906zT1jMCg4Ht4DTDBzC6JXhiCfClETin+\nj7v3BKqIHPquE4KMmFkaMAJ46eRlTSGfGv+ZbQVKouaLg7Gw2GFmhQDB885gvEnmNrNUIk3/OXd/\nJRgOVUaA4NDpTOBqwpXvImCEmW0gclrtMjObTIgyuvvW4Hkn8CqRw76hyUdkj3ZLcDQKYBqRDYEw\nZYTIhtsCd98RzDepfGr8Z/YJUGZmpcEW3mhgepxrakjTgZuC6ZuA16LGR5tZupmVAmXAx3Go72sz\nMyNyXnG5uz8WtSgUGc2swMxyg+kM4EpgBSHJB+DuE9292N07EPld+6u7f4+QZDSzLDPLOTYNDAWW\nEJJ8AO6+HdhsZucFQ5cDywhRxsB3OX6YH5pavni/u7CxP4BhRN4hvha4J971/AM5nge2AUeIbJWP\nBfKB94DVwLtAXtT69wSZVwLXxLv+r5FvEJHDa4uBT4PHsLBkBLoDC4N8S4D7gvFQ5DtF3sEcf1d/\nKDIS+XTQouCx9Njfk7Dki6q5BzAv+Fn9E9AyTBmBLKASaBE11qTy6cp9IiIiCUSH+kVERBKIGr+I\niEgCUeMXERFJIGr8IiIiCUSNX0REJIGkxLsAEWm8zOzYx5QA2gBHiVySFeCAuw+MS2EiUm/6OJ+I\nfC1m9jNgv7v/Z7xrEZH606F+EakXM9sfPA82s/fN7DUzW2dmD5vZGDP7OLj3fKdgvQIze9nMPgke\nF8U3gUhiUuMXkYZQDowHugA3Aue6e18it9f9cbDOfwP/5e59gOs58da7IhIjOscvIg3hEw9uS2pm\na4F3gvHPgCHB9BVA18htFQBobmbZ7r4/ppWKJDg1fhFpCIejpmuj5ms5/ncmCejv7odiWZiInEiH\n+kUkVt7h+GF/zKxHHGsRSVhq/CISK7cBFWa22MyWEXlPgIjEmD7OJyIikkC0xy8iIpJA1PhFREQS\niBq/iIhIAlHjFxERSSBq/CIiIglEjV9ERCSBqPGLiIgkEDV+ERGRBPJ/57LbHzM0BQ8AAAAASUVO\nRK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig,ax=plt.subplots(1,1,figsize=(8,6))\n", "time = range(720)\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": 20, "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": 21, "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": 23, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAf4AAAGDCAYAAADK03I6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcHHWd//HXp3vuI5lMjiF3QhLCfSXhFhJRQUTCssrG\nBUSXwxPx57ECu4oXir/1/rmoLCAoQkRUYBFERAJyhJAQICQQct/3NZlJ5v78/qia0JnM0ZPMTHdV\nv5+PRz+muqq66/Pp7ulPf7/1rSpzd0RERCQ3JDIdgIiIiPQdFX4REZEcosIvIiKSQ1T4RUREcogK\nv4iISA5R4RcREckhKvyCmd1tZt/O0LbNzH5lZjvMbE6GYqgxs8O7WGeMmbmZ5fXC9p83s5PC6ZvM\n7I40H5f2uh08fqqZrT3Yx3fx3L32enUjhioze9PMCjtZ5xdm9tW+jKu3pPM57itm9i4zW5zpOKR9\nKvxZyMxWmtlmMytNmXe1mc3KYFi95SzgvcAIdz+l7UIz+5iZNYdfatVm9qqZXXiwGzOzWWZ2deo8\ndy9z9+UH+5yHwsw+COx29/lhLN9x96u7eBht1810oQ0/s+/J0LbfY2avmFmtma01s0sB3H0T8DRw\nbUePdfdPuvu3DnK7s8yszsx2h5/NeWZ2Q2c/NHpTb32O02kYhJ+98Smx/MPdJ/Z0LNIzVPizVxK4\nPtNBdJeZJbv5kNHASnev7WSdF929DKgA7gQeMLMB3YzLzCwbP++fBH6T6SCymZlVdbLsaOA+4D+A\n/sAJwLyUVX4LfKIXw/usu5cDQ4EvAjOAx8zMenGbIofG3XXLshuwErgB2A5UhPOuBmaF02MAB/JS\nHjMLuDqc/hjwPPAjYCewHDgjnL8G2AxcmfLYu4FfAE8Cu4FngNEpy48Ml20HFgOXtnnsz4HHgFrg\nPe3kMwx4JHz8UuCacP5VQB3QDNQA32jnsR8Dnku5XxrmPhkYADwKbAF2hNMj2rwmt4SvxV6CItAc\nbrMG+Fm4ngPjw+li4AfAKmAX8Fw4b7/XnKDI3AlsANYB3waS4bLx4Wu4C9gK/K6D97kgjCs15q8D\n97Z5n68EVofP9R8drLs6XLcmvJ3ezvaKw/drB7AI+DKwts379Ifw9VwBfK7Nth4Efhd+Rl4BTgiX\n/QZoCXOpAf69q9jT+B/IB/4p/Nzs7mS9+4BvdbI8D9hDyue5zfK7gW+H01OBtQQFfHP43n68k+ee\nRfg/lzJvVLi9C8P7hcCPgfXh7cdAYZvt/XvK9i4GLgDeJvh/uSnluU8BXiT4n94A/AwoSFme+jm+\nG/hv4M/h+/USMK6TXH4PbCT4zD4LHBPOvxZoBBrC9/Z/23nss+G2a8N1/qU1tzbfaV8GXg/XuxOo\nAh4P4/sbMCBl/dOAF8JcXwOm9tT3a/j8Z4SfyZHh/RMI/i+O7MntZOst4wHo1s6bEvyTvAf4Y8qX\nUncLfxPwcYKeg28TfPn+d/hF9L7wn60sXP/u8P7Z4fKfEBZbgkK7JnyuPOCk8B/m6JTH7gLOJOhB\nKmonn2eB24Ai4ESCwvLulFif6+S1+FhKLHkEvSC7CQrvQOCfgRKgPPzyeqjNa7IaOCZ8bD7tf1mn\nfmH+d7jO8PC1OyN8TfZ7zYE/Ab8MX58hwBzgE+Gy+wlaoIkw57M6yO0YoLbNvK9zYOH/H4KifQJQ\nDxzVybp5nbyWtwL/ACqBkcAbhF/OYazzgK8R/CA5nOAH43kp22oEPhS+jl8i+HGQn/qZTdlWp7F3\nEuNxwA8JCuGLBK31ik7WXw58C1hAUAzvBSrbrPM6cFEHj7+b/Qt/E/DNMMcLCIr4gA4ee8BnKeXz\n/r1w+pvA7PAzMpigmH2rzfa+Fm7vGoL/jfsIPs/HEPyYGhuuP4mgIOaFr++bwOc7+BzfDWwj+LGQ\nR/Cjd2Ynr+O/hdts/aHyanuvUSeP37ftlNzaFv7ZBMV+ePj+vkLwfVIE/B24OVx3eBj7BQSfy/eG\n9wd3sO1HCX4gtHd7tJOYbwm3Wxx+fj6bzvdzHG4ZD0C3dt6Udwr/sQRFdTDdL/xLUpYdF65flTJv\nG3BiOH136pcCUEbQMh5J8Ov9H23i+2XKP+ndwK87yWVk+FzlKfO+C9ydEmtXhb8p/CfeGn55HNCr\nEK57IrCjzWvyzTbr7HudUuY5QSs9QfBFe0I7z73vNQ+/vOqB4pTlHwGeDqd/DdxOSku+g3jPBDa2\nmfd1DizmqT0Cc4AZnazbWeFfDpyfcv9a3in8pwKr26x/I/CrlG3NTlmWICi070r9zLbzerUbezux\nvRuYS/Aj8zvAEWn+rzSE2z6C4HP7B+C3bdZ5HvhoB4+/m/0L/172/7/aDJzWwWMP+CyF82cC/xNO\nLwMuSFl2HsGurdTttfYUlYev2akp688DLu5g+58H/tT2c5yS1x0pyy4A3krzNa0In6t/29eok8ek\nU/gvS7n/B+DnKfevI/zRDnwF+E2b53+ClF7KnrgR/NiaR1D0/wJYTz5/Nt8yNuJWuubub5jZowTd\n/m928+GbUqb3hs/Xdl5Zyv01KdutMbPtBF2/o4FTzWxnyrp57L9feg0dGwZsd/fdKfNWEXTVp2u2\nu5/VdqaZlRDszjifoNsfoNzMku7enEZsbQ0iaH0s62K90QRfGhtSduUmUrb17wSt0DlmtgP4gbvf\n1c7z7CD4su/KxpTpPez/vnXHMPZ/PValTI8GhrV5n5MEPQStUj8jLeERAcO62Ga6sQ8h+PH1IkHX\n7uounrfVXoIfJ28DmNl3CLqNU5UT/HBMxzZ3b0oz5o4MJ2jZQ/D6pL7Oq9j/NduW8lndG/5t9//U\nzI4g6A2ZTNDLlcf+4xnaSuu1D8fl3AJ8mKCR0RIuGkTQ8OgpbfPq6PtoNPDhcOBrq3yCgZo9xt0b\nzexu4KfAFzz8NZALsnGwk+zvZoIuwOEp81oHwpWkzDvsELczsnXCzMoIuoPXE3zZP+PuFSm3Mnf/\nVMpjO/uHWQ9UmllqgRtFsF/8UH0RmEjQQupHsKsCIHVgVdvYOot1K8H+/3FdbHcNQYt/UMpr0s/d\njwFw943ufo27DyPoqr4tdcRziqUE4w6Ht7Osu9L50tpAyvtM8D60WgOsaPM+l7v7BSnrpH5GEsAI\ngvc33e13yN1nEnyGf0Mw9mO9mf2PmR3wg6+N19tse784wqMcxhP8mOh1ZjaSoEu+9QfTeoJC1moU\n77xm3fVz4C1gQvh5v4n9P+sH61+B6QS9jP0JemtIee6+LohrCFr8qZ/FUne/tb2Vzezx8Kif9m6P\nd7SR8P/uZuBXwA8ydTRGJqjwZzl3X0owoOpzKfO2EBTOy80saWb/RtfFqisXmNlZZlZA0Fqd7e5r\nCPafHWFmV5hZfnibYmZHpRn/GoLWz3fNrMjMjif4Yr/3EOOFoCW3F9hpZpUE/8Rd2USw/7q9WFuA\nu4Afmtmw8LU9ve0XgrtvAP5K8GXRz8wSZjbOzM4BMLMPm9mIcPUdBF+cLbTh7g0ErdNz0km2C1vC\nbXR2HPcDwI1mNiCM77qUZXOA3Wb2FTMrDnM/1sympKwzycwuCYvp5wl+/MwOl3X4uqbL3evc/T53\nfx/BmICVwK/MrLMemF8BHzezw8MeoBsIPrOtTiHoWl/V7qN7iJmVhO//wwSv5WPhovuB/zSzwWY2\niGB//sF+9suBaqDGzI4EPtXF+t153nqC3X8lBLtaUqXz3h7y+5/iXuCDZnZe+DkssuCcEyPaW9nd\n3x82Rtq7vb+9x4RHXdxNMMjwKoIfxQd1WGcUqfBHwzcJBpGluoZglOw2gkFAL7R9UDfdR1A4txO0\nWC4HCLvo30dwmNJ6gu7D7xEMAkrXRwhaEesJBsXd7O5tu2MPxo8JBua07vv/SxqP+QnwIQtOGPTT\ndpZ/iWCf38sEr8X3aP//5KMEg+AWERT3BwkO6QKYArxkZjUEo9Kv946Pr/4lcEUacXfK3fcQHsFg\nZjvN7LR2VvsGQVfzCoIfLr9JeXwzcCHBOIkVBK/pHQQtwFYPE4z52BHGfIm7N4bLvktQ4Haa2Zd6\nIJ817n6Lu08gODKgo/XuIhhT8VKYWz0pP5KBywiOWOktPzOz3QSF78cE+67PD39EQjCwdi5Bz8QC\nggFtB3uyrC8RtM53Ewya/N0hxJ3q1wSv3TqCz/PsNsvvBI4O39uHOniOrwP3hOtceijBhI2F6QQ9\nGlsIegC+TM/Wq88R7F76atjF/3GCH5Dv6sFtZC3Lod0aIlnJzJ4nGFE8P9OxdMTMvk4weOvyTMeS\nLjMbQnBY5UnuXpfpeESyhQb3iWSYu5+Z6RjiyN03A2ntkhLJJerqFxERySHq6hcREckhavGLiIjk\nEBV+ERGRHBLbwX2DBg3yMWPG9Njz1dbWUlra9oi6+Ih7fhD/HOOeH8Q/x7jnB/HPMZP5zZs3b6u7\nD+5qvdgW/jFjxjB37twee75Zs2YxderUHnu+bBP3/CD+OcY9P4h/jnHPD+KfYybzM7O0TlSlrn4R\nEZEcosIvIiKSQ1T4RUREcogKv4iISA5R4RcREckhKvwiIiI5RIVfREQkh6jwi4iI5BAVfhERkRyi\nwi8iIpJDeq3wm9ldZrbZzN5ImVdpZk+a2ZLw74CUZTea2VIzW2xm56XMn2RmC8JlPzUz662YRURE\n4q43W/x3A+e3mXcD8JS7TwCeCu9jZkcDM4BjwsfcZmbJ8DE/B64BJoS3ts/Z69ydNbtbaG7xvt60\niIhIj+q1wu/uzwLb28yeDtwTTt8DXJwyf6a717v7CmApcIqZDQX6uftsd3fg1ymP6TM3/nEBX31+\nLy8s29rXmxYREelRfX11vip33xBObwSqwunhwOyU9daG8xrD6bbz22Vm1wLXAlRVVTFr1qxDDrjF\nnRfeqgNg9rzXaF4Xzwsa1tTU9Mjrlc3inmPc84P45xj3/CD+OUYhv4xVMXd3M+vRvnN3vx24HWDy\n5MneU5dGHHH0bt73o2eZeNTRTD1hWI88Z7aJ+6UyIf45xj0/iH+Occ8P4p9jFPLr61H9m8Lue8K/\nm8P564CRKeuNCOetC6fbzu9T+cngZWpuaenrTYuIiPSovi78jwBXhtNXAg+nzJ9hZoVmNpZgEN+c\ncLdAtZmdFo7m/2jKY/pMXiI4kKCxWYP7REQk2nqtq9/M7gemAoPMbC1wM3Ar8ICZXQWsAi4FcPeF\nZvYAsAhoAj7j7s3hU32a4AiBYuDx8NanWlv8TSr8IiIScb1W+N39Ix0sOreD9W8Bbmln/lzg2B4M\nrdvykkGLv0ld/SIiEnE6c18aWrv61eIXEZGoU+FPQ15rV79a/CIiEnEq/GnQ4D4REYkLFf40aHCf\niIjEhQp/GsIGv47jFxGRyFPhT4OZkTRo1EV6REQk4lT405RMQFOzWvwiIhJtKvxpSho0qcUvIiIR\np8KfpqRpcJ+IiESfCn+akgnTcfwiIhJ5KvxpSpqO4xcRkehT4U9T0qBZ+/hFRCTiVPjTFLT41dUv\nIiLRpsKfpuBwPrX4RUQk2lT405Q0De4TEZHoU+FPUzKh4/hFRCT6VPjTpOP4RUQkDlT406TBfSIi\nEgcq/GlSV7+IiMSBCn+agsF9KvwiIhJtKvxpSpiuziciItGnwp+mPB3HLyIiMaDCn6akQaOO4xcR\nkYhT4U9TQufqFxGRGFDhT1PSTF39IiISeSr8aVKLX0RE4kCFP00Jg2ZX4RcRkWhT4U9T0qBFLX4R\nEYk4Ff40menMfSIiEn0q/GlKqMUvIiIxoMKfpqRa/CIiEgMq/GlKmGlwn4iIRJ4Kf5rU1S8iInGg\nwp+mhLr6RUQkBlT405Sw4K9a/SIiEmUq/GlqLfzazy8iIlGmwp+mfYVfLX4REYkwFf40qfCLiEgc\nqPCnKWlB5VdXv4iIRJkKf5paX6hmXZpXREQiTIU/TabBfSIiEgMq/GlK6nA+ERGJARX+NLUO7tNJ\nfEREJMpU+NOkUf0iIhIHKvxpUuEXEZE4UOFPU0KH84mISAyo8KdJg/tERCQOVPjTZBrcJyIiMaDC\nnybt4xcRkThQ4U/Tvq5+7eMXEZEIU+FPk47jFxGROMhI4Tez/2NmC83sDTO738yKzKzSzJ40syXh\n3wEp699oZkvNbLGZnZeJmBMa3CciIjHQ54XfzIYDnwMmu/uxQBKYAdwAPOXuE4CnwvuY2dHh8mOA\n84HbzCzZ13HvO5xPhV9ERCIsU139eUCxmeUBJcB6YDpwT7j8HuDicHo6MNPd6919BbAUOKWP49Xg\nPhERiYW8vt6gu68zs+8Dq4G9wF/d/a9mVuXuG8LVNgJV4fRwYHbKU6wN5x3AzK4FrgWoqqpi1qxZ\nPRZ3/d69gPHKq6/RsLbPOxx6XU1NTY++Xtko7jnGPT+If45xzw/in2MU8uvzwh/uu58OjAV2Ar83\ns8tT13F3N7NuN63d/XbgdoDJkyf71KlTDz3g0NKHngLqOOa445g6cUiPPW+2mDVrFj35emWjuOcY\n9/wg/jnGPT+If45RyC8TXf3vAVa4+xZ3bwT+CJwBbDKzoQDh383h+uuAkSmPHxHO61Ma3CciInGQ\nicK/GjjNzErMzIBzgTeBR4Arw3WuBB4Opx8BZphZoZmNBSYAc/o45n3H8Wsfv4iIRFkm9vG/ZGYP\nAq8ATcB8gu75MuABM7sKWAVcGq6/0MweABaF63/G3Zv7Om6N6hcRkTjo88IP4O43Aze3mV1P0Ppv\nb/1bgFt6O67OtHaN6Op8IiISZTpzX5p0OJ+IiMSBCn+aVPhFRCQOVPjTpMIvIiJxoMKfpoSuzici\nIjGgwp8mXZ1PRETiQIU/Ta2H8+kEPiIiEmUq/GlKqsUvIiIxoMKfJg3uExGROFDhT1Ne+Eo1NLdk\nNhAREZFDoMKfptbC39SsFr+IiESXCn+aEmaYQaNa/CIiEmEq/N2Qn0zQqBa/iIhEmAp/N+QnTC1+\nERGJNBX+bsjPS9Ckwi8iIhGmwt8NeYkEDerqFxGRCFPh74aCpKnFLyIikabC3w15yYT28YuISKSp\n8HdDftJo1Jn7REQkwlT4uyE/maCxSS1+ERGJLhX+bshPJnSRHhERiTQV/m7IT+o4fhERiTYV/m7Q\n4D4REYk6Ff5uKNApe0VEJOJU+LshT8fxi4hIxKnwd0N+UmfuExGRaFPh7wYN7hMRkahT4e+G/KQu\n0iMiItGmwt8NeQkN7hMRkWhT4e+Ggjx19YuISLSp8HdDvo7jFxGRiFPh74a8RIImdfWLiEiEqfB3\nQ36e0aAWv4iIRJgKfzfkJ9TVLyIi0abC3w3lRXm0OOyua8x0KCIiIgelW4XfzBJm1q+3gsl2h/Uv\nAmBTdV2GIxERETk4XRZ+M7vPzPqZWSnwBrDIzL7c+6Fln6p+QeHfuKs+w5GIiIgcnHRa/Ee7ezVw\nMfA4MBa4olejylL7Cr9a/CIiElHpFP58M8snKPyPuHsjkJPHtB3WT139IiISbekU/l8CK4FS4Fkz\nGw1U92ZQ2aq4IMnA0gKWb6nNdCgiIiIHpcvC7+4/dffh7n6BB1YB0/ogtqx06uGVvLBsK+452ekh\nIiIRl87gviozu9PMHg/vHw1c2euRZakzxg1iw6461mzfm+lQREREui2drv67gSeAYeH9t4HP91ZA\n2e6ooeUALN2yO8ORiIiIdF86hX+Quz8AtAC4exPQ3KtRZbFxg8sAWLZZ+/lFRCR60in8tWY2kHAk\nv5mdBuzq1aiyWEVJAYPKCli6uSbToYiIiHRbXhrrfAF4BBhnZs8Dg4EP9WpUWW5YRbGO5RcRkUjq\nsvC7+ytmdg4wETBgcXgsf84qK8yjtr4p02GIiIh0W4eF38wu6WDREWaGu/+xl2LKeqWFeWyv3ZPp\nMERERLqtsxb/BztZ5kDOFv7ywjxq1OIXEZEI6rDwu/vH+zKQKClVV7+IiERUOifw+Y6ZVaTcH2Bm\n3+7dsLJbqVr8IiISUekczvd+d9/ZesfddwAX9F5I2a+8KI/GZqe+KWdPZyAiIhGVTuFPmllh6x0z\nKwYKO1m/S2ZWYWYPmtlbZvammZ1uZpVm9qSZLQn/DkhZ/0YzW2pmi83svEPZdk8oLUgCUFuvwi8i\nItGSTuH/LfCUmV1lZlcBTwL3HOJ2fwL8xd2PBE4A3gRuAJ5y9wnAU+H91msDzACOAc4HbjOz5CFu\n/5CUFeUDUFOn7n4REYmWdI7j/56ZvQ6cG876lrs/cbAbNLP+wNnAx8LnbwAazGw6MDVc7R5gFvAV\nYDow093rgRVmthQ4BXjxYGM4VGWFwe8O7ecXEZGosb6+vKyZnQjcDiwiaO3PA64H1rl7RbiOATvc\nvcLMfgbMdvd7w2V3Ao+7+4PtPPe1wLUAVVVVk2bOnNljcdfU1FBWFpyn/42tzXx/bh03nlLExMqM\ndj70mNT84iruOcY9P4h/jnHPD+KfYybzmzZt2jx3n9zVep2dwOc5dz/LzHYTnqe/dRHg7t7vIGPL\nA04GrnP3l8zsJ4Td+q3c3c2s279I3P12gh8VTJ482adOnXqQIR5o1qxZtD7fwLW7+P7c5xg78Rim\nHnNYj20jk1Lzi6u45xj3/CD+OcY9P4h/jlHIr8N9/O5+Vvi33N37pdzKD6HoA6wF1rr7S+H9Bwl+\nCGwys6EA4d/N4fJ1wMiUx48I52VMZVkBANtrGzIZhoiISLelcxz/b9KZly533wisMbOJ4axzCbr9\nHwGuDOddCTwcTj8CzDCzQjMbC0wA5hzs9nvCwNKg8G9T4RcRkYhJ5+p8x6TeMbM8YNIhbvc64Ldm\nVgAsBz5O8CPkgfDIgVXApQDuvtDMHiD4cdAEfMbdM3ocXVF+kpKCpFr8IiISOZ3t478RuAkoNrPq\n1tlAA+F+9IPl7q8C7Q1AOLedebj7LcAth7LNnlZZWqDCLyIikdPZPv7vAv2BX7fZvz/Q3W/suxCz\n08DSAnX1i4hI5HS6j9/dW4ApfRRLpAQt/vpMhyEiItIt6Zy57xUzU/Fvo7K0kB21jZkOQ0REpFvS\nGdx3KnCZma0CannnOP7jezWyLDewrIBtavGLiEjEpFP4M35RnGxUWVpAXWMLexqaKClI52UUERHJ\nvC67+t19lbuvAvYSnMGv9ZbTKluP5a/RAD8REYmOdE7gc5GZLQFWAM8AK4HHezmurNd6Eh8d0ici\nIlGSzuC+bwGnAW+7+1iCY+1n92pUETBAhV9ERCIoncLf6O7bgISZJdz9ado/+U5OGVRaCMCWGg3w\nExGR6EhnVNpOMysDniU4ze5mgtH9OW1oRRF5CWPVtpx/KUREJELSafFPB/YA/wf4C7AM+GBvBhUF\n+ckEoweWsHRzTaZDERERSVunLX4zuxgYDyxw9yeAe/okqogYM7CUJxZu4rEFG7jguKGZDkdERKRL\nHbb4zew2glb+QOBbZvbVPosqIqZOHAzA5+6fT019U4ajERER6VpnLf6zgRPcvdnMSoB/EIzwl9AV\np4+hpCCPL/7+NbburqesUCfyERGR7NbZPv6G1uveu/seglP1Shv9i/MBqK7TeftFRCT7ddZEPdLM\nXg+nDRgX3te5+lP0LwkK/669KvwiIpL9Oiv8R/VZFBHWr0iFX0REoqPDwh+en1+6sK+rf68G94mI\nSPZL5zh+6US/4uC3k1r8IiISBSr8h6g4P0l+0jS4T0REIqGz4/ifCv9+r+/CiR4zo39xvlr8IiIS\nCZ0N7htqZmcAF5nZTNoczufur/RqZBFSXpRPtQq/iIhEQGeF/2vAV4ERwA/bLHPg3b0VVNSUFCTZ\n29Cc6TBERES61Nmo/geBB83sq+6uM/Z1orQgjz0q/CIiEgFdnmPW3b9lZhcRnMIXYJa7P9q7YUVL\nSWGSHbUNmQ5DRESkS12O6jez7wLXA4vC2/Vm9p3eDixKSgvyqFWLX0REIiCdq8p8ADjR3VsAzOwe\nYD5wU28GFiXF2scvIiIRke5x/BUp0/17I5AoKy1IUtugM/eJiEj2S6fF/11gvpk9TXBI39nADb0a\nVcSUFOaxp14tfhERyX7pDO6738xmAVPCWV9x9429GlXElBYkaWhuobG5hfykToYoIiLZK50WP+6+\nAXikl2OJrJKC4GXc09BM/2IVfhERyV6qUj2gpCAJwB7t5xcRkSynwt8DSgqDFn+t9vOLiEiW67Tw\nm1nSzN7qq2CiqqwwaPHv1hX6REQky3Va+N29GVhsZqP6KJ5IGlZRDMC6nXszHImIiEjn0hncNwBY\naGZzgNrWme5+Ua9FFTEjB5QAsGrbngxHIiIi0rl0Cv9Xez2KiCstzGNQWSFrtqvwi4hIdkvnOP5n\nzGw0MMHd/2ZmJUCy90OLltEDS3hj/S7cHTPLdDgiIiLtSuciPdcADwK/DGcNBx7qzaCi6J9PHsEb\n66qZ9faWTIciIiLSoXQO5/sMcCZQDeDuS4AhvRlUFF1y8nASBvNX7ch0KCIiIh1Kp/DXu/u+i82b\nWR7gvRdSNBXlJxk3uIxFG6ozHYqIiEiH0in8z5jZTUCxmb0X+D3wv70bVjQdPawfb27YnekwRERE\nOpRO4b8B2AIsAD4BPAb8Z28GFVXjBpexbude9jboDH4iIpKd0hnV32Jm9wAvEXTxL3Z3dfW3Y8yg\nUgBWba/lyMP6ZTgaERGRA6Uzqv8DwDLgp8DPgKVm9v7eDiyKxg4MCv/KrbVdrCkiIpIZ6ZzA5wfA\nNHdfCmBm44A/A4/3ZmBRNHpQcAa/lTqDn4iIZKl09vHvbi36oeWARrC1o19RPiUFSbbsrs90KCIi\nIu3qsMVvZpeEk3PN7DHgAYJ9/B8GXu6D2CJpcHmhCr+IiGStzrr6P5gyvQk4J5zeAhT3WkQRN7hM\nhV9ERLJXh4Xf3T/el4HExZB+hby9qSbTYYiIiLSry8F9ZjYWuA4Yk7q+LsvbvsFlhTy3ZGumwxAR\nEWlXOqP6HwLuJDhbX0tPbdjMksBcYJ27X2hmlcDvCH5grAQudfcd4bo3AlcBzcDn3P2Jnoqjpw3p\nV0R1XRN7G5opLtBFDEVEJLukM6q/zt1/6u5Pu/szrbce2Pb1wJsp928AnnL3CcBT4X3M7GhgBnAM\ncD5wW/j5pjNyAAAeZUlEQVSjISuNCY/lX6Fj+UVEJAulU/h/YmY3m9npZnZy6+1QNmpmI4APAHek\nzJ4O3BNO3wNcnDJ/prvXu/sKYClwyqFsvzeNGxIU/qVbtJ9fRESyj3V19l0z+y5wBcHZ+1q7+t3d\n333QGzV7EPguUA58Kezq3+nuFeFyA3a4e4WZ/QyY7e73hsvuBB539wfbed5rgWsBqqqqJs2cOfNg\nQzxATU0NZWVlXa7X0Ox84sk9XDQun3+aUNBj2+9t6eYXZXHPMe75QfxzjHt+EP8cM5nftGnT5rn7\n5K7WS2cf/4eBw1MvzXsozOxCYLO7zzOzqe2t4+5uZt2+HoC73w7cDjB58mSfOrXdpz8os2bNIt3n\nGznvaZpL+zN16iF1jPSp7uQXVXHPMe75QfxzjHt+EP8co5BfOoX/DaAC2NxD2zwTuMjMLgCKgH5m\ndi+wycyGuvsGMxuasr11wMiUx48I52WtcYNLWbZF+/hFRCT7pLOPvwJ4y8yeMLNHWm8Hu0F3v9Hd\nR7j7GIJBe39398uBR4Arw9WuBB4Opx8BZphZYXho4QRgzsFuvy+MH1LG8i01rN+5l6cXb6auUZfp\nFRGR7JBOi//mXo8icCvwgJldBawCLgVw94Vm9gCwCGgCPuPuWV1Jxw8po76phTNu/TsA3774WC4/\nbXSGoxIREUmj8PfQoXsdPfcsYFY4vQ04t4P1bgFu6a04etr5xwzlrudWUpif4PW1u1imEf4iIpIl\n0jlz326Ci/MAFAD5QK279+vNwKKsf0k+f/7cWZgZH/jpP1ity/SKiEiWSKfFX946HR5mNx04rTeD\nioO8ZDB8YlRliU7mIyIiWSOdwX37eOAh4Lxeiid2RlWWsHr7Hro6X4KIiEhfSKer/5KUuwlgMlDX\naxHFzODyQuqbWqipb6K8KD/T4YiISI5LZ1T/B1OmmwguoDO9V6KJocrS4Ox922sbVPhFRCTj0tnH\n//G+CCSuBpYFhX9bbQOjwwv4iIiIZEqHhd/MvtbJ49zdv9UL8cROZWkhANtreuSMxyIiIoeksxZ/\ne0PRS4GrgIGACn8aBrZ29e9R4RcRkczrsPC7+w9ap82sHLge+DgwE/hBR4+T/aXu4xcREcm0Tvfx\nm1kl8AXgMuAe4GR339EXgcVFSUGSwryECr+IiGSFzvbx/xdwCcFlbo9zd5139iCYGQNLC9imffwi\nIpIFOjuBzxeBYcB/AuvNrDq87Taz6r4JLx4qywrYXluf6TBEREQ63cffrbP6SccqSwvV1S8iIllB\nxb0PDCwtYJsKv4iIZAEV/j5QWVqgFr+IiGQFFf4+UFlawJ6GZuoamzMdioiI5DgV/j4wsrIEgHte\nWJnZQEREJOep8PeBDxw3lA8cP5TvPv4Wdz23grrGZl5fu1M9ACIi0ufSuTqfHKJkwvjhpSfw9sbd\nfPPRRXz7z4tocThqaD/+9OkzKMpPZjpEERHJESr8faQwL8mDnzqDvy7cyML11TS3OL+ZvYpXVu3g\njPGDMh2eiIjkCBX+PtS/OJ8PTx7Jh4GNu+r4zexVLNtaq8IvIiJ9Rvv4M6SqXyGlBUmWbdaZkEVE\npO+o8GeImTFuSBnLtqjwi4hI31Hhz6Bxg8vU4hcRkT6lwp9B4waXsn5XHbX1TZkORUREcoQKfwaN\nG1wGwIJ1u3D3DEcjIiK5QIU/g44d3h+AGbfP5qY/vaHiLyIivU6FP4NGVpbwi8tPZlRlCffPWc0d\n/1iR6ZBERCTmVPgz7Pxjh/LMl6dyxriB3P3CSlpaglZ/fVOzegBERKTHqfBnATPjslNHs27nXq67\nfz5zV27nxG88yW2zlmU6NBERiRkV/izx/mMP4z1HVfHnBRv40C9eZG9jM//v70toblGrX0REeo4K\nf5ZIJIw7rpzMe46qoiCZ4J9PHkFdYwtvbqjOdGgiIhIjOld/lvnZv55EY3MLO/c08odX1jJ/9Y59\no/9FREQOlVr8WaYoP0l5UT4jBhQzuLyQ+at3ZjokERGJERX+LGVmnDSygldW78h0KCIiEiMq/Fns\n5NEDWLltD1t212c6FBERiQkV/ix2zhGDSSaMz/9uPnWNzZkOR0REYkCFP4sdNbQf//Wh43l+6TYu\n/eWLfOXB16nRBX1EROQQqPBnuUtOHsHHzhjD62t38bu5a7jgJ/9g+RZdyldERA6OCn8EfP2iY1jw\n9fdx39WnUlvfxOV3vMTWGu33FxGR7lPhj4jyonzOGD+IOz82hfW76vjf19ZnOiQREYkgFf6IOXFk\nBaMHlvDckq2ZDkVERCJIhT+C3nd0FbPe3sITCzdmOhQREYkYFf4I+uy0CYwdVMonfjOPxxZsyHQ4\nIiISISr8EdS/JJ9HrzuLE0dWcNOfFvDo6+tpCa/it6m6jlod8iciIh3QRXoiqig/yQ8vPYFP3juP\nz943n3OOWMuwimIemLuGSaMHMPOa00gkLNNhiohIllGLP8IOH1zG49efzZfPm8gzb2/hwXlrOHVs\nJXNWbGf+Gp3jX0REDqQWf8QlE8anzhnH8IpiJo0eQFF+kim3/I35q3cyaXRlpsMTEZEso8IfA4mE\ncfFJw/fdH9a/iFfX6HK+IiJyIHX1x9DJowfw8srtuHumQxERkSzT54XfzEaa2dNmtsjMFprZ9eH8\nSjN70syWhH8HpDzmRjNbamaLzey8vo45as4aP4hN1fUs3axz+ouIyP4y0eJvAr7o7kcDpwGfMbOj\ngRuAp9x9AvBUeJ9w2QzgGOB84DYzS2Yg7sg4Z2JwOd+f/n0pjy3YwPbahkyHJCIiWaLPC7+7b3D3\nV8Lp3cCbwHBgOnBPuNo9wMXh9HRgprvXu/sKYClwSt9GHS1D+xcz/cRh/O9r6/n0b1/hktueZ3N1\nXabDEhGRLJDRwX1mNgY4CXgJqHL31tPQbQSqwunhwOyUh60N50knvvNPxzF14hD2NjTxjf9dxLk/\neIZhFcV8+5+OZcoYjfYXEclVlqkBYGZWBjwD3OLufzSzne5ekbJ8h7sPMLOfAbPd/d5w/p3A4+7+\nYDvPeS1wLUBVVdWkmTNn9li8NTU1lJWV9djz9aVV1c08uaqJVzY1MbI8wcXjCzhq4P57S6KcX7ri\nnmPc84P45xj3/CD+OWYyv2nTps1z98ldrujufX4D8oEngC+kzFsMDA2nhwKLw+kbgRtT1nsCOL2r\nbUyaNMl70tNPP92jz5cJ33v8TR/9lUd99Fce9V89t3y/ZXHIrytxzzHu+bnHP8e45+ce/xwzmR8w\n19OowZkY1W/AncCb7v7DlEWPAFeG01cCD6fMn2FmhWY2FpgAzOmreOPkC+89gqe+eA7vmjCIH/1t\nCTU6p7+ISM7JxKj+M4ErgHeb2avh7QLgVuC9ZrYEeE94H3dfCDwALAL+AnzG3ZszEHfk5SUTjBtc\nxvXnTmDX3kYe15X9RERyTp8P7nP354COrh5zbgePuQW4pdeCyjGTRg9gxIBiHn19Ax+ePDLT4YiI\nSB/SmftykJlx4fHDeG7pVh3jLyKSY1T4c9QHjhtKc4vzzNubMx2KiIj0IRX+HHXk0HLyEsaSTTqt\nr4hILlHhz1H5yQSjB5awbIsKv4hILlHhz2HjBpfpQj4iIjlGhT+HTR4zgGVbalmwdlemQxERkT6i\nwp/DZpwyigEl+Xz2/ld4Y2szO/dohL+ISNyp8OewfkX53PWxKWzZXc/359Zx3o+fZcOuvZkOS0RE\nepEKf447adQAnvriOVx3UiHVe5u4+eGFmQ5JRER6UUYvyyvZYWj/YiZV5WGVo/jpU0tYta2W0QNL\nMx2WiIj0ArX4ZZ9LJ48gL2FM/f4sTv3O33hw3tpMhyQiIj1MhV/2GTGghPuvPY1r3nU4wyuK+dLv\nX2POiu2ZDktERHqQCr/sZ8qYSm664Ch+e/VpDCor4CdPvU1wmWcREYkDFX5pV3FBks9MG8/zS7fx\nZ12+V0QkNlT4pUMfPX0Mxw3vz80PL+SNdTrJj4hIHKjwS4eSCePHM07EDC657QXe3FCd6ZBEROQQ\nqfBLp8YNLuPPn3sX5UV5XHHnSzz86rpMhyQiIodAhV+6VNWviF9eMYmmFuf6ma/ynw8tyHRIIiJy\nkFT4JS2Tx1Qy60tTOWpoP+6dvZp5q3SYn4hIFKnwS9oqSgp48JOnU1GSz89nLc90OCIichB0yl7p\nltLCPD52xhh+/LclvLFuF4+/sYF5q3ZwwXFDmTFlFAV5+i0pIpLNVPil2648fQy3P7ucC//fcwCM\nGVjC1x5eyEPz13HWhMFcddZY+hfnZzhKERFpjwq/dNuA0gJ+c9Up3PfSGs4+YhAXnTCM3728hv96\nYjHz1yzhsQUbuO2ykzmiqjzToYqISBsq/HJQJo2uZNLoyn33Z5wyihmnjOKZt7fwpd+/xmV3vMRV\nZ43FgLxkAndn2pFDGDe4LHNBi4iICr/0rHOOGMx9V5/K5Xe+xK2Pv7Xfsv96YjH3Xn0qU8ZUdvBo\nERHpbSr80uMmVJXz/FfeTUNzC00tzpvrqxlYVsiVd83h+08s5nefOB2Amvom8hJGUX4ywxGLiOQO\nFX7pFXnJBHnJYIT/qYcPBOAjp4zk+399mxv+8DrVdY08tmAjBXkJPnraaD4zbTwDSgsyGbKISE7Q\nsVfSZy4/bTQDSwuY+fIa/vLGRi4/bRQXHj+Uu55fwdn/92nu+MdyXQJYRKSXqcUvfaaipIA/fvoM\nGptbGFlZQmFe0MX/yXPGcevjb/HtP79J/+J8Pjx5JO6OO+xpbKa0IImZZTh6EZF4UItf+tTogaWM\nH1K+r+gDHFFVzh0fncz4IWX88ZV17NrTyFX3zOXwmx7j2Juf4Ed/W5LBiEVE4kWFX7JCImG8/9jD\neHH5Nk745l/5+1ubOfKwcoryE9zzwko2VddlOkQRkVhQV79kjU9NHUd9Uwt/W7SJ733oeKaMqWT2\n8m1ccedLXHLbC/zrqaNYtKGa/sX5/NuZYxk/ROcEEBHpLhV+yRolBXncdMFR3HTBUfvmnXb4QP7w\nqTO45tdz+a8nFjOsfxE79zbyu5fXcM27DufUwyuZNnFIBqMWEYkWFX7JesePqODZf5/Glt31DK8o\nZnttA9fdP59fPLOMXzyzjItPHMaFxw/jPUdXUV3XyP88u5w5K7Zz22UnM7CsMNPhi4hkFRV+iYTC\nvCQjBpQAMLCskN9efSo19U1857E3eWj+eh56dT3DK4qpbWhi555GEgY/ePJtvvNPx2U4chGR7KLB\nfRJJZkZ5UT7fveR45n/tvRw9tB/VdY1MHl3Jo9edxYcmjeC+l1bzzNtbANhd18iWPS0ZjlpEJPPU\n4pfIK8pP8sdPn0HCjIK84LfsZaeO5oG5a/m3u19mSHkhG3YFRwUsaHyDb04/NpPhiohklAq/xELb\n8/2fMLKCV7/2Xn7296Ws3FZL/+J8Sn0Pv35xFSu37aGsMMmAkgKue/cEfv3iSn71/EoSFpxq+PLT\nRnHt2ePoX5wPQHOLk0wYTc0tJBOmkwmJSKSp8EtsVZQU8J8XHr3v/pN/f5pnqwcze/k2AFZuq+W3\nL60G4ORRFZQW5lGUn+S/n17GH+at40f/ciLVdY187v75nDFuIAvW7WJo/2K+Mf0YTh41ICM5iYgc\nKhV+yRn5CeNbF7/Tzf/qmp38+fX1FOYl+cJ7jyCRsH3z/+3ul/nI/8wGYGBpAXNWbGdvYzP1TS3M\nuH02wyuKmVhVzlsbqynISzBlTCUFeQkmjR6wr6dge20Dr67ZydkTBnP0sH4s21LDpNEDWLdjLzv2\nNHLSyIp92xQR6Ssq/JKzThxZwYkjK9qd/9QXzmH28m20OEwZO4DKkgLMjG219eHugz3MWbmd40f0\nZ3ddE4++voG9jc386vmVBzxf6rxBZQVsrWkAIJkwkgljUGkBx4+owAzW7tjLks27OWFEBQPLCpgw\npJxPTR1HUX6S5hZnd10jZYV5+658KCLSXSr8Iu0YUFrA+48besD8IeVFHQ4OrGtsZuH66n1XGEwk\njCMPK+fR1zawftdexg4q5eFX1zPxsHLGDixl+dZaHGfR+mqWbalhyeYaAC4+cRhvbdzNlt31PLZg\nIw+/uo6BZYWs2raHrTX1HD6olFv/+fh9PQut1u5uYfHG3eQnjbGDSjUWQUTapcIv0kOK8pNMGn3g\nvv9Lp4zcNz39xOEdPn7p5t3s2tu033M88tp6fj93De7BOIQTRlbwi2eWcekvX2z/SZ5/FoCTRlVw\nzbsOp7wojyOqyqnqV3SQWYlI3Kjwi2SJ8UPKD5h30QnDuOiEYfvN+5cpI3l5xXa8zboLFy7kmGOO\nYVN1Hbc/u5xP//YVAAqSCc6ZOJgzxw1kw646Vm/fw5hBpYwbvP+1DoZXFHPkYeW0dhQkEkZhXgL3\nA4+aEJHoUuEXiZhBZYXt7oYo2baYqeH8j5wyioXrq2lsbuGh+et46q3NPLloEwmDUZUl/GXhRrzt\nL4cOFCQTHDWsH8mUPQfHj6hgcPk7p0PuV5zPmeMGUpCXoLa+udPnc5zX1+5iy+56ivOTjKwsYVRl\nCf2L8ykpTLJh54FXYkwmjMMHlaYXsIh0SoVfJIZSdzucdvhAmluc9Tv3UlqYR2VpAVtr6tnb8E6B\nbnHnldU72Lmncd+87bUNbKttoKXFWbdz7775exqauefFlWn/cOiOgmSChub2z7A4sLQAa2mk6KW/\nk59MMGXMANxh7qodmMEpYyopyk9y9hGDmDZxCGbGiq21DC4vpKxQX3UirfTfIJIDkgljZGXJvvuD\n2rl40eiB6beoG5tbaEmp/Es21TB7+TYam53RA0s6eWRgSHkhEw8rZ9H6amrqm6hrbOGtjdXs3NPI\n5DEDyG9z1MLWmnpeW7OLjRs3UnVYJZur6/nrok0AHDe8P03Nzl8WbqShqYW7X1hJaUGScUPKeH3t\nLgCK8hOUFeZx7PD+FKfsthhVWcJZEwZx1vhBGgwpOUOFX0S6rW1hPnZ4f44d3r/bz3Pq4QP3TX/g\n+AN3X+zndJg1awdTp57Y4SqNzS38af46nn17C8u31HLVWWMpyk+wevteVm+rZX1Kz0VzS/Bj4ZfP\nLufwwaUM7V/ElDGV++U2fkhZWj9kAAaUFOwbRLlu51721DcxfkiZflBI1lHhF5HYyE8muHTySC6d\nPLLrlQl2Zzy5aCOPLdjIym21PL90yUFv2wyG9ivCYd+1IaafOIwpYypZsmk31XVNzF21nZYWKMhL\ncMqYSoYPKOaoof14eeV26hubWbuunqd3vdHFdozjR/Tfr9dmYzigc1hFMYcPKqUoP8mEIWUMH1DM\n5NEDdN4H2Y8Kv4jkrMrSAv5lyij+Zcoo3H2/8QXNLc5LK7ZT19D5YEUABxatr95X8EdWFrOpuo77\n56zh4VfXU1qQpDA/yUkjK6goKWBj9V6eWLRx35iKvIRRVpRHY2MT+VvXd7qtusZm7n7hwHEQg8sL\n2VPfxOzl22hoeme5WfD8B6s4P8kpYyvJSyTY09jMyyu209SS3pUuh1cUM/Gwcox3tr9lax0z18zb\nb70BpQWcNLKCnugcOWpoP4ZXFOPAWxuqWbdzL+OHlHH44LIDzn1R19jM3oZmCvISlObQOJDcyVRE\npBNmRmHe/octTps4JO3HX9DmSAt357JTR2MGRw/t126X/5rte9i8u45xg8uoKClg1qxZTJ06tdPt\nNDa3sHB9Nc1tiu8xw/pTlJ+krrGZbbUNbN1dz6rte3hrQ3XaObRn1bY9LA1PLgXBLpkh5QeOEWmr\nxeGV1TtYuXXPfvNra1uooXbffcdZtXgP989ZfUhxpmP8kDLKi/JoaXEWbaimsTkYp2IGxwzrR14i\nwbItNfQvzt/vqJXuqN61l58sej6tdS88fhhXnTX2oLZzKCJT+M3sfOAnQBK4w91vzXBIIiIdMrMu\nxz2MrCzZb9BlOvKTiXZPNd2qKD/J8IpihlcUc8LIigPOA5FpwY+bs/ebV1PfxI7ahkN+7vqmFuas\n2E5DU9BLM6C0gGOH92feyh2s27mXV9fs3Dco9V0TBnNY/yImDCljzfbgVNkAk0YHR4u0HORhK415\nlvZRJIV5mdkFE4nCb2ZJ4L+B9wJrgZfN7BF3X5TZyERE5FCVFeb12CGX44eUHTCv7cmqelPww+bU\nPtvewYjKiI9TgKXuvtzdG4CZwPQMxyQiIhI5USn8w4E1KffXhvNERESkG8x74/RbPczMPgSc7+5X\nh/evAE5198+2We9a4FqAqqqqSTNnzuyxGGpqaigr67vuor4W9/wg/jnGPT+If45xzw/in2Mm85s2\nbdo8d5/c1XqR2McPrANSD8wdEc7bj7vfDtwOMHnyZO9qdGx3pDPaNsrinh/EP8e45wfxzzHu+UH8\nc4xCflHp6n8ZmGBmY82sAJgBPJLhmERERCInEi1+d28ys88CTxAczneXuy/McFgiIiKRE4nCD+Du\njwGPZToOERGRKItKV7+IiIj0ABV+ERGRHKLCLyIikkNU+EVERHKICr+IiEgOUeEXERHJIZE4Ze/B\nMLMtwKoefMpBwNYefL5sE/f8IP45xj0/iH+Occ8P4p9jJvMb7e6Du1optoW/p5nZ3HTOgRxVcc8P\n4p9j3POD+OcY9/wg/jlGIT919YuIiOQQFX4REZEcosKfvtszHUAvi3t+EP8c454fxD/HuOcH8c8x\n6/PTPn4REZEcoha/iIhIDlHh74KZnW9mi81sqZndkOl4DpaZ3WVmm83sjZR5lWb2pJktCf8OSFl2\nY5jzYjM7LzNRp8/MRprZ02a2yMwWmtn14fxY5GhmRWY2x8xeC/P7Rjg/FvmlMrOkmc03s0fD+7HJ\n0cxWmtkCM3vVzOaG82KTH4CZVZjZg2b2lpm9aWanxyVHM5sYvnett2oz+3zk8nN33Tq4AUlgGXA4\nUAC8Bhyd6bgOMpezgZOBN1Lm/V/ghnD6BuB74fTRYa6FwNjwNUhmOocu8hsKnBxOlwNvh3nEIkfA\ngLJwOh94CTgtLvm1yfULwH3AozH8nK4EBrWZF5v8wrjvAa4OpwuAirjlGMaeBDYCo6OWn1r8nTsF\nWOruy929AZgJTM9wTAfF3Z8FtreZPZ3gn5Tw78Up82e6e727rwCWErwWWcvdN7j7K+H0buBNYDgx\nydEDNeHd/PDmxCS/VmY2AvgAcEfK7Fjl2I7Y5Gdm/QkaGXcCuHuDu+8kRjmmOBdY5u6riFh+Kvyd\nGw6sSbm/NpwXF1XuviGc3ghUhdORztvMxgAnEbSKY5Nj2AX+KrAZeNLdY5Vf6MfAvwMtKfPilKMD\nfzOzeWZ2bTgvTvmNBbYAvwp319xhZqXEK8dWM4D7w+lI5afCL0DQoiT4Uoo0MysD/gB83t2rU5dF\nPUd3b3b3E4ERwClmdmyb5ZHOz8wuBDa7+7yO1ol6jsBZ4Xv4fuAzZnZ26sIY5JdHsEvx5+5+ElBL\n0PW9TwxyxMwKgIuA37ddFoX8VPg7tw4YmXJ/RDgvLjaZ2VCA8O/mcH4k8zazfIKi/1t3/2M4O1Y5\nAoRdp08D5xOv/M4ELjKzlQS71d5tZvcSoxzdfV34dzPwJ4Ju39jkR9CiXRv2RgE8SPBDIE45QvDD\n7RV33xTej1R+KvydexmYYGZjw194M4BHMhxTT3oEuDKcvhJ4OGX+DDMrNLOxwARgTgbiS5uZGcF+\nxTfd/Ycpi2KRo5kNNrOKcLoYeC/wFjHJD8Ddb3T3Ee4+huB/7e/ufjkxydHMSs2svHUaeB/wBjHJ\nD8DdNwJrzGxiOOtcYBExyjH0Ed7p5oeo5Zfp0YXZfgMuIBghvgz4j0zHcwh53A9sABoJfpVfBQwE\nngKWAH8DKlPW/48w58XA+zMdfxr5nUXQvfY68Gp4uyAuOQLHA/PD/N4AvhbOj0V+7eQ7lXdG9cci\nR4Kjg14Lbwtbv0/ikl9KzCcCc8PP6kPAgDjlCJQC24D+KfMilZ/O3CciIpJD1NUvIiKSQ1T4RURE\ncogKv4iISA5R4RcREckhKvwiIiI5JC/TAYhI9jKz1sOUAA4DmglOyQqwx93PyEhgInLQdDifiKTF\nzL4O1Lj79zMdi4gcPHX1i8hBMbOa8O9UM3vGzB42s+VmdquZXWZmc8Jrz48L1xtsZn8ws5fD25mZ\nzUAkN6nwi0hPOAH4JHAUcAVwhLufQnB53evCdX4C/MjdpwD/zP6X3hWRPqJ9/CLSE1728LKkZrYM\n+Gs4fwEwLZx+D3B0cFkFAPqZWZm71/RppCI5ToVfRHpCfcp0S8r9Ft75nkkAp7l7XV8GJiL7U1e/\niPSVv/JOtz9mdmIGYxHJWSr8ItJXPgdMNrPXzWwRwZgAEeljOpxPREQkh6jFLyIikkNU+EVERHKI\nCr+IiEgOUeEXERHJISr8IiIiOUSFX0REJIeo8IuIiOQQFX4REZEc8v8B/kA7Yo0QzBkAAAAASUVO\nRK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig,ax=plt.subplots(1,1,figsize=(8,6))\n", "time = range(720)\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 }