{ "cells": [ { "cell_type": "code", "execution_count": 34, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import netCDF4 as nc\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from salishsea_tools import viz_tools, tidetools, geo_tools\n", "from matplotlib import animation, rc\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 38, "metadata": { "collapsed": true }, "outputs": [], "source": [ "particles = nc.Dataset('/ocean/vdo/MEOPAR/ariane-runs/NanimoOilSpill/ariane_trajectories_qualitative.nc')" ] }, { "cell_type": "code", "execution_count": 39, "metadata": { "collapsed": true }, "outputs": [], "source": [ "rc('animation', html='html5')" ] }, { "cell_type": "code", "execution_count": 40, "metadata": { "collapsed": false }, "outputs": [], "source": [ "Bathy = nc.Dataset('/data/vdo/MEOPAR/NEMO-forcing/grid/bathymetry_201702.nc')" ] }, { "cell_type": "code", "execution_count": 118, "metadata": { "collapsed": true }, "outputs": [], "source": [ "test3 = nc.Dataset('/ocean/vdo/MEOPAR/ariane-runs/test3/ariane_trajectories_qualitative.nc')" ] }, { "cell_type": "code", "execution_count": 122, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "" ] }, "execution_count": 122, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmwAAAG7CAYAAAB+RskvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlYlOX+BvD7YRGRxUTQJBc4Ku4IiYqiIcvgvpSpmaSl\nkQKVnexYmrmXLaf0VOLayXKXtNJOJsOiuCsKuO9Z7oqgLLLP9/eHOr9MUVTgHeD+XNdcDTPv+z7f\nmXOU2+d9FiUiICIiIiLTZaZ1AURERER0fwxsRERERCaOgY2IiIjIxDGwEREREZk4BjYiIiIiE8fA\nRkRERGTiGNiIyGQppcYrpRZqXQcRkdYU12EjIq0opTL/8mM1ALkACm/9PFJElpZ9VQ9HKRUAYDaA\n+gB2AnhZRP7QtioiqmjYw0ZEmhER29sPAH8C6P2X18pDWHMEsAbABwAcACQAWKlpUURUITGwEZHJ\nUkpNVkotufXcRSklSqlXlFJnlFJpSqlRSqm2Sql9SqlrSqmv/3b+cKXU4VvHblBKNSjhEp8DcFBE\nIkUkB8BkAK2VUk1LuB0iquQY2IiovGkPoDGAQQBmAXgfQCCAFgAGKqV8AUAp1RfAeNwMVU4ANgNY\nXtRFbwW+oh7vFXFaCwDJt38QkSwAJ269TkRUYiy0LoCI6CFNu9WbFaWUygKwXEQuA4BSajMATwCb\nAIwCMENEDt967yMA45VSDe41xkxEnniEWmwBXPnba+kA7B7hWkRERWIPGxGVN5f+8jz7Hj/b3nre\nAMB/bveSAUgFoAA8VYK1ZAKw/9tr1QFklGAbREQMbERUYZ3BzZmmT/zlYS0i2+51sFIq8z6P8UW0\ncRBA679cwwZAw1uvExGVGAY2Iqqo5gIYp5RqAQBKqepKqQFFHfzXGav3eHxUxGk/AmiplOqvlKoK\nYBKAZBE5UuKfhogqNQY2IqqQRORHAJ8AWKGUSgdwAED3Em7jCoD+AD4EkAagHYAXSrINIiKAC+cS\nERERmTz2sBERERGZuGIv66GUMsfNVbzPiUgvpVRr3BwjYgvgNIAhIpJ+j/NO4+aMqUIABSLiVQJ1\nExEREVUaD9PDNhrA4b/8vBDAeyLSCjcH3v7rPuf6iYgHwxoRERHRwytWYFNK1QXQEzdD2m1uAOJv\nPdfj5sBbIiIiIiphxb0lOgvAWNy5evdBAH0B/ARgAIB6RZwrAKKVUoUA5onI/HsdpJR6DcBrAGBj\nY9OmadPH24ovOzsb6enpSE9PR2ZmJgwGAwDA1tYW9vb2sLe3x5kzZyAiaNas2SO1ceTIzZn7j1vr\n/Rw7dgyZmZlwd3eHmZkZkpKSUKtWLdStW7dU2ktLS8OpU6fg4OCA2rVro1q1aiV27cLCQiQlJaFO\nnTpwdnYusetWdhcvXsS5c+fg6OiIGjVqwNbWFmZmd/5bLD8/H4cOHYKFhQWeeuop2NnZwdzc3Pj+\n0aNHUVhYiObNmxtfMxgMyMjIwLlz55CXl4dmzZrBysqqzD7Xo/rzzz9x9epVeHh4QCmldTlEVMHt\n378feXl5KSLiVKoNich9HwB6AYi49bwLgF9uPW8KIArAHtxce+hqEec/deu/tXBzz71nHtRmmzZt\npCTl5ORIbGysjBs3Try8vEQpJbgZJOX9999/pGumpKSIUkomTZpUorX+3dSpU0UpJZcvX5b//e9/\nAkCioqJKtc133nnH+P3UrFlTBg0aJAsXLpQ//vjjsa67atUqASDbtm0roUpJRKSwsFBGjBghVlZW\nAkCqVKkifn5+8tFHH8nu3bulsLBQRERiYmLE3t5eAIi5ubl06NBBJk6cKHq9XiwsLOTdd9+95/V/\n//13eeKJJ6RNmzaSk5NTlh/tkTRu3FiCgoK0LoOIKoFDhw7d/n2ZIA/INo/7KE5gmwHgLG5OLLgI\n4AaAJX87xg3ArmJcazKAdx50XEkHtr9LSUmRlStXypgxY+TMmTOPdI3ly5cLANm+fXsJV3cnHx8f\n8fLyEhGRN998U6ytrSU7O7tU2xQRuXDhgixevFiGDh0qderUMQa4Jk2ayOuvvy5r166V9PT0h7rm\n8OHD5YknnpD8/PxSqrpyy8rKkg0bNsg777wj7u7ud4TugQMHyoIFC+T48eOyceNGef/996Vdu3Zi\nZmZmPC42NrbIa//0008CQMLDw8vwEz2a559/XiwtLWX37t1al0JEFdyKFStMJ7DdcfCdPWy1bv3X\nDMD3AIbf43gbAHZ/eb4NQLcHtVPaga0kjB07VgDI2bNnS62NtLQ0MTc3lwkTJoiISPPmzSUgIKDU\n2iuKwWCQ/fv3yxdffCHdu3eXatWqCQCxsLCQTp06yZQpU+T48eMPvIazs7MMGDCgjKqmixcvypIl\nS2TYsGF3hG43NzcJDw+XU6dOSWpqqvzwww8yc+ZMY09cUd5++20BIKtWrSqjT/BoUlJSpF69euLq\n6ippaWlal0NEFVh+fr4sW7bM5APbaADHbj0+xv8vwusM4Ndbz/9x6zZoMm6OeXu/OO2Uh8B24MAB\nsba2Fj8/PykoKCiVNn744QcBIJs3bxYRkWeffVasrKwkMTGxVNorrtu3mN977z1p06aNKKXE0dHx\nnr2VGRkZMm/ePPHw8BAAsmjRIg0qJoPBIAcOHJCZM2dKjx49xNraWpo2bSoZGRnFvkZeXp54e3uL\nnZ3dAwO61rZv3y4WFhbSr18/MRgMWpdDRBWcyQW2snqUh8AmIvLtt98KAJk8eXKpXD8kJETs7e0l\nLy9PREQuX74szs7O0rhx44e+HVmaDhw4IDY2NtKpUyfj7c6DBw/K66+/bhwz5e7uLvPmzXtgLw6V\njdjYWDEzM5Pg4OCHCjR//PGHODg4iJubm8ycOVMOHDhgsoHo888/FwAyc+ZMrUshogqOga0cGDZs\nmCilJDo6ukSvazAYpF69etK/f/87Xo+PjxczMzN54YUXTOoX5ZIlS4y323x9fY2D34cMGSJbt241\nqVrppqlTpwoAWbBgwUOdFx0dLW5ubsZbrM7OzjJ06FBZvHixXLlypZSqfXgGg0H69u0rFhYWpT7W\nlIgqNwa2ciAzM1OaNWsmtWvXlgsXLpTYdQ8ePCgAZP78+Xe999FHHwkAmTt3bom1VxLCwsIEgLi6\nusonn3wily9f1rokuo+CggIJDAyUqlWrSnJy8kOf/8cff8jChQtl0KBBUrNmTeP/9qYkNTVVXFxc\npH79+pKVlaV1OURUQZVFYDPJzd+9vLwkISFB6zKK7eDBg2jXrh3at28PvV5/x/pWj+qLL77AmDFj\n8Mcff6B+/fp3vGcwGNCjRw9s3LgR27dvh6en52O3V1LS0tJQvXr1u9YBI9N0+fJleHh4wM7ODgkJ\nCbCzs3vwSfdgMBjQsmVLODo6Ij4+/sEnlKFPP/0U7777Lv7880/Uq1fUcpFERI9OKbVHSnk3J/5W\nLQEtWrRAREQE4uLi4OXlhffffx8bN25Ebm7uI19zw4YNsLW1xd69e3H9+vU73jMzM8PixYvh6OiI\ngQMHIj39ri1cNVOjRg2GtXKkVq1aWL58OU6cOIGQkBBkZmY+0nVSUlJw+PBhdO3atYQrfHzx8fFo\n1KgRwxoRlWv8zVpChg0bhtmzZ8PGxgaffPIJ/Pz84ODggB49emDWrFkPHarat28PEcGzzz6LmjVr\nwsfHB5MnT8bWrVtRUFAAJycnrFixAr///jtCQkJgij2lVD74+vpi2rRpWLlyJRwcHNClSxd8+OGH\n2L17NwoLC4t1jaioKAAwucCWk5ODuLg4dOvWTetSiIgeC2+JloL09HRs3LgRUVFR0Ov1OHbsGPz9\n/REVFfVQt0vz8vKwY8cO6PV6REVFISEhAQaDAXZ2dvD394dOp8PevXvx3//+F7Nnz0ZYWFgpfiqq\nyEQEcXFx+O2336DX65GUlATgZo9pnz59sGDBAlhaWhZ5/ksvvYTffvsNly5dMqke1ujoaOh0Oqxb\ntw69evXSuhwiqqDK4pYoA1sZWLhwIUJCQjB58mRMmjTpka+TlpaGmJgY6PV66PV6/P7773e8f/jw\n4VLd15Qqj8uXLyMmJgaffPIJ9u/fj2vXrhU5vs1gMKBOnToIDAzE0qVLy7jS+5s6dSomTZoEX19f\n9OjRA0FBQcZ9eYmISgoDWwUhIhg2bBiWLFmCtWvXomfPniWyKfXJkycRFRWFFStWID4+HoWFhfxF\nRCWqffv2MDc3x7Zt24o8Zu/evWjTpg2+//57vPTSS2VY3YP9+eefmDVrFvR6PQ4cOAAAcHJyQmBg\nIIYNG2Zyt3CJqHxiYKtAMjMz0a5dOxw+fBh169aFTqdDUFAQAgIC4OTkpHV5RHdJSUlBrVq1MHny\nZEycOLHI42bMmIHx48fj4sWLqF27dhlW+HDOnz+P6OhoYw91SkoK4uLi0LlzZ61LI6JyjoGtgklN\nTcXq1asRFRWF6OhoXLt2DQDg6emJoKAg6HQ6+Pj4oGrVqhpXSgSsWLECgwcPxo4dO9C+ffsij/P1\n9UVGRgb27t1bhtU9nvT0dHh5eSErKwtJSUn8RxMRPRYu61HBODg4ICQkBJGRkUhJScHOnTsxbdo0\n2Nra4vPPP0dgYCAcHBzQrVs3fP7559i3bx9nf5Jmbt8GHTduHD755BMkJibCYDDccUxBQQG2bduG\ntm3balHiI7O3t8eqVatw9epVBAcH3/W5iIhMDXvYTERGRgY2bdpkvF1z+PBhAMCTTz6JwMBABAUF\nITAwEHXq1NG4UqosTp48iYiICERFRd0x/isgIACjRo2Cr68vACAoKAibN2/Gzp074e7urmXJD23+\n/PkYOXIkpk+fjvfff1/rcoionOIt0Urs7NmzxvB2e7wNALRq1co4/q1z586oVq2axpVSZXDhwgXj\n+K/ly5ejS5cu0Ov1AIBLly7Bw8MD9vb2j7VbghZEBMHBwVixYgViYmLQpUsXrUu6LxHB/v37Ua9e\nPdSoUUPrcojoFgY2AnBz2YTk5GTjum5btmxBbm4uqlSpgk6dOhnHv3l4eHCWKJWqq1evwsnJCZMm\nTbpjiZpNmzbB398fgwYNwtKlS0tkFnRZyczMhJeXF65fv45t27bB1dVV65Lucu3aNXz33XeYO3cu\njhw5AjMzM3h5eRn/7Ht7e6NKlSqa1CYi2LhxI/773//CzMwMgYGB0Ol0ePLJJ+84rqCgAP/73/8Q\nFRUFT09P6HQ6NGjQQJOaiUpaWQQ2zTd6v9ejPG3+roWsrCz57bffZMyYMeLu7i4ABIA4OjrKN998\no3V5VIEtX75cAMiOHTvuem/69OkCQObOnatBZY9n3759Ym1tLQCkZcuW8s9//lPWr19fahvGHz9+\nXMaMGSMtWrSQIUOGyKJFi+TcuXN3HZeQkCAjRoww1ta+fXuZN2+eTJo0STp06CDm5uYCQGxsbKR3\n795y6tSpUqn3Xq5duyZffvmlNGvWTABIjRo1pGbNmsa/j1q1aiVvv/22/PLLLzJt2jSpW7euAJAq\nVaoYj3Fzc5Pw8HD56aef5Pr162VWO1FJQxls/q55OLvXg4Ht4Vy4cEEWL14sPj4+YmFhcc9fpkQl\n4eWXXxYHBwcpKCi4673CwkLp2rWrWFlZSWJiogbVPZ6jR4/Kxx9/LAEBAWJlZWUMF35+fjJjxgxJ\nSEiQwsLCR75+fn6+/PjjjxIUFCQAxMLCQvz8/KRWrVrGANO8eXMZPXq0fPXVV9KuXTsBINWqVZNX\nX31V9u7de9c109LSZM2aNRIaGir29vbi4eEh2dnZdxxjMBgeueZ72bNnj7z66qtSrVo1ASDt2rWT\nRYsWyY0bN6SwsFD27NkjM2bMEH9//zvCmU6nkzVr1kheXp4cOHBAZs6cKd27dzdex9zcXHx8fGTy\n5MmydetWyc/PL9G6iUoTAxs9lNTUVHFxcZEGDRrI1atXtS6HKqAGDRpIr169inx/3759YmFhIW5u\nbiUeFMpSVlaWbNiwQd555x1p3bq1MXTUrFlTBg0aJAsXLpQ//vijWNc6f/68TJ061djDVLduXZk6\ndaqcP39eRG4G3aSkJPn0009Fp9NJ1apVBYA0bdpUvvzyS0lLSytWO7/88osAkJEjR4rBYJDY2FgZ\nMGCAWFlZiYeHh/zrX/8SvV5/R6C7HSK7du0qLVu2lDfffFPWrVsn6enpd1z7xo0bsmjRImnfvr0A\nEGtraxkxYoQkJCQ88HuMiYmR48ePF3lMTk6OxMXFyfjx46Vt27ailBIAYm9vL3379pVvv/22XP9/\niSoHBjZ6aDt37hRLS0vp06cP/5KjEjdw4MC7enENBoPExcUZ3wMgzz777D174cqrixcvypIlS2TY\nsGHi7Ox8xy29119/XX7++ec7bun9NTDd/k6CgoLkxx9/fGDP0Y0bN+Tw4cOP9Od37NixAkBcXFyM\ntylfffVV8fX1FUtLSwEgVatWFZ1OJ2+99ZYxRD711FMSGBhovPVqYWEhnTt3lmnTpsmYMWPEwcFB\nAEiTJk1k1qxZxQ6RjyIlJUVWrVolISEhxs8xa9asUmuPqCQwsNEjmTVrlgCQzz//XOtSqIK53Ytb\nv359OXXq1F1jmN5++205duyY1mWWKoPBYLyl16NHD+MtPQsLC/Hx8ZG3335bmjZtqsl3kpeXJz17\n9hRvb2/jbcrbMjIy5JdffpHRo0dL8+bNRSklOp3ujhCZnZ0t0dHR8u6778rTTz9t/FwDBgyQ2NjY\nMv9HoMFgkD59+oilpaXs3LmzTNsmehhlEdg4S7QCEhH0798f69atw+bNm+Ht7a11SVSB7N69G+3a\ntTP+3LZtW4SFhWHgwIGVcpmZ3NxcbN++HXq9HlFRUdizZw/atWuH0NBQDBw4ENbW1lqXeE95eXkP\nnFl6ezkhR0fHsijpntLS0uDp6Qng5r61Dg4OmtVCVJRKu6yHlZWVDBs2DDqdDgEBAfwD+giuXbsG\nT09PFBYWIjExETVr1tS6JKpA9u3bhy1btqBdu3bw8irdmezlTXGCED2cXbt2oVOnTujevTt++umn\ncrVsDFUOlTaw1ahRQwwGA9LT06GUgpeXF3Q6HXQ6HTp27Kj5X4b79+/HqlWr4OrqCp1Oh3r16mla\nT1ESEhLQsWNHdO3aFT///DPXaCOicus///kP3nrrLfz73//GmDFjtC6H6A6VNrB5eXnJjh07sGvX\nLuNisTt37kRhYSGqVasGX19f44KRzZs3v+e/tkQEiYmJAFAiC8rm5uZi9erVmDNnDrZs2XLHe02a\nNDHW06VLF5Na6f2rr77Cm2++iU8//RT/+te/tC6HiOiRiAief/55rF27FvHx8ejQoYPWJREZVerA\n9vcxbNevX8fGjRuN40SOHz8OAHB2dr5jr01bW1ssW7YMERERSEpKAgDUrFkTAQEBxlBVv379Ytdy\n+vRpzJs3D9988w2uXLmChg0bIjQ0FMOGDcPFixeNgXLTpk3Izs6GhYUFvL29jW15eXnBwsKi5L6c\nhyQiGDhwIH788Uds2bKF49mIqNy6du0a2rRpg7y8PCQlJXGoB5kMBrb7+OOPP4z7bEZHRyM1NRUA\nYG1tjezsbLi7uyM0NBS2trbG4y5cuAAAcHNzM95i9fPzg729/R3XNhgM2LBhAyIiIvC///0PSin0\n7t0boaGh0Ol09+yty83NxbZt24xt7dmzByKC6tWrIyAgwNhew4YNS+hbKr5z586hbt26mDJlCiZO\nnFjm7RMRlZS9e/eiQ4cOCAwMxLp16zjUg0wCA1sx3R5Yr9frcfbsWbz44ovo2LHjHbdKRQQHDx40\n9tDd7hEzNzeHt7e38Xbmzp07MXfuXPz++++oXbs2QkJC8Nprrz30OLWUlBTExsYiKioKUVFROHPm\nDADgH//4hzG8+fv7l8kGzmvXrkXfvn0RGxsLPz+/Um+PiKg0RUREIDw8HB9//DHeffddrcshYmAr\nTUX1iAGAr68vwsLC0K9fvxKZ4CAiOHbsmLGt2NhYZGZmwszMDG3btoVOp0NQUBC8vb1haWn52O39\nXXh4OL777jukpqZqPmGDiOhxiQgGDx6MH374AXFxcejcubPWJVElx8BWhq5evYotW7agcePGaN68\neam2lZ+fj507dxp7+3bt2gWDwQBbW1t06dLFOP6tSZMmJTJ9vVGjRmjevDnWrl1bAtUTEWkvPT0d\nXl5eyMrKQmJiImrVqqV1SVSJMbBVEteuXUNcXJyxB+7EiRMAgLp16xp73wICAuDk5PTQ1z579izq\n1auHGTNm4L333ivp0omINJOcnAxvb2907twZ69evh7m5udYlUSXFwFZJ/f7773dMqLh27RoAwNPT\nE/369cP48eOLPfM0NzcX7u7uyMzM5L9CiajCWbhwIUJCQjB16lR88MEHWpdDlVSlDWxOTk5ye7Hc\noKAgtGjRotKubF1YWIg9e/YYJy9s3rwZ48aNw0cffVTsayQnJ6N9+/bw9fXF+vXrOauKiCoMEcHQ\noUOxbNky6PV6+Pv7a10SVUKVNrA5OjrK1atXjT/XqVPHuNbagAEDYGVlpWF12goJCcHChQuxfv16\ndOvWrdjnzZ8/HyNHjsS0adMwYcKEUqyQiKhsZWZmom3btkhLS0NiYiLq1KmjdUlUyVTawPb000/L\njRs3cPr0abz77rv46quvkJaWBuDm9iRvvvmmxhVqJzs7G97e3jh37hySkpJQt27dYp0nIggODsaK\nFSsQGxsLX1/fUq6UiKjsHDx4EG3btkX79u2h1+s1XbCcKp+yCGwmeW/MzMwMkZGRKCgowNSpU5GW\nloZmzZrhyy+/REhIiNblacra2hqRkZHIzc3FCy+8gIKCgmKdp5TC3Llz0ahRIwwePBiXL18u5UqJ\niMpOixYtMGfOHGzcuBFTpkzRuhyiEmeSgQ0AWrVqhd27d2P27NmIi4vDwYMH8cYbb8Da2lrr0jTn\n5uaG+fPnY+vWrQ91e9POzg6RkZFIS0vDa6+9VooVEhGVvWHDhuGVV17Bhx9+iKioKK3LISpRJnlL\ntLLPEi2uUaNGYd68eVi8eDEGDRpU7EV3g4KCcObMGRw+fLiUKyQiKls3btxA+/btcfHiRZw/f75U\nFiMn+rtKe0uUimfmzJnw9PTESy+9BAcHB/Tp0wdff/01jh49iqKCeE5ODrZu3YrAwMAyrpaIqPRV\nq1YN/v7+yMzMRGFhodblEJUYjsosx6ytrbF582Zs2LDBuGvCunXrAAD169c37lkaEBAAR0dHAMCW\nLVtw48YNdO3aVcvSiYhKzYYNG9ClSxdUrVpV61KISgxviVYwp06dMoa32NhYXLt2DUopeHp6QqfT\nISoqCvv27cO1a9dga2urdblERCXq9OnTcHV1xcyZM/HWW29pXQ5VEmVxS5Q9bBXMP/7xD4wcORIj\nR45EYWEhEhISjAHu888/R0FBAUaOHMmwRkQV0oYNGwDcHKt7P6mpqahevTq3s6Jyg2PYKjBzc3O0\nb98eEyZMQHx8PFJTU7Fr1y7MnTtX69KIiErFoUOHAAB+fn4YMmQIFi1ahHPnzgEACgoK8NNPPyEo\nKAg1a9aEo6Mj+vfvj7lz5+LUqVPGa+Tk5OD7779Hhw4dYG9vj549e2LWrFk4ePBgkeODiUobb4kS\nEVGFkZGRgTVr1hj3Y7695mSzZs2QkZGBs2fPom7duhg2bBguXryIqKgonDlzBsDNOxRt27ZFdHQ0\nrl69iiZNmqBTp06Ij4/H8ePHAQDOzs7G8cGBgYGoXbu2Zp+VTEel3emAgY2IiB6XwWDA/v37jeHN\n0tISr776Knr16mXcCUFEcOzYMeMxO3bsQOfOnREWFgY/Pz/jPtanT582HhMTE4PU1FQAQOvWrY0B\nrnPnzlwrtJJiYCMiIjIxhYWFSExMNI4P3rp1K/Lz82FlZYUXX3wRCxYs4Ni4SoaTDoiIiEyMubk5\nvLy84OXlhXHjxiErKwvx8fFYs2YNFi5cCBcXF0ycOFHrMqmCYWAjIiJ6DDY2NujevTu6deuGnJwc\nTJ48GT4+PggICNC6NKpAOEuUiIioBCilMGfOHDRt2hSDBw/GZ599huTkZBgMBq1LowqAgY2IiKiE\n2NraYvXq1ahTpw7Gjh0LDw8P1KlT564lRogeFicdEBERlYJz584hOjr6riVGmjdvbpxZ6uvry4XM\nKwDOEiUiIqoAbi8xEhUVBb1ej82bNyMnJweWlpbo2LEj+vTpg9GjR3N2aTnFwEZERFQB5eTkYMuW\nLcalQZKSkjBp0iRMnjxZ69LoETCwERERVXAigpdffhmLFy9GVFQUAgMDtS6JHlJZBDZOOiAiItKQ\nUgoRERFo1qwZhgwZgh07dqCwsFDrssjEMLARERFpzMbGBpGRkcjOzkaHDh1Qq1YtDBw4EAsWLMDp\n06e1Lo9MAAMbERGRCWjevDlOnTqFpUuXok+fPti2bRtee+01uLq6ws3NDfHx8VqXSBriGDYiIiIT\nJCI4fPgw9Ho9vvrqK2RkZCApKQl16tQxHnPhwgU4OjrC0tJSw0qJY9iIiIgqKaUUmjdvjtGjR2Pt\n2rXIzMzE4MGDkZWVhSVLlsDHxwfOzs6oWbMm+vXrh9mzZ+PYsWMwxY4YenzcS5SIiMjENW/eHHPm\nzMGwYcPg5OSE7OxsuLm5Yfr06Thz5gyioqLw888/AwDq169vXJg3ICAAjo6OGldPJYGBjYiIqBwY\nOnQojh49iuPHj2PkyJHw9/eHUsr4/smTJ427KqxevRrffPMNlFJ4+umnjQHOx8cHVlZWGn4KelQc\nw0ZERFTBFBQUYM+ePcadFbZv346CggJYW1vjmWeeMQa4Vq1a3RH66NFw4VwiIiJ6bBkZGdi0aZMx\nwB05cgQA8OSTTyIwMNAY4P46oYEeLDU1FYsWLcKYMWMY2IiIiKhknTlz5o6N6VNSUgAALVu2NIa3\nZ555BjY2NhpXappOnjyJadOmYeXKlcjJyQEABjYiIiIqPQaDAcnJycbwtnnzZuTm5qJKlSrw8fEx\nBjhPT09uTn+Lj48PkpKSMHToUISGhqJ169YMbERERFR2srOzsXnzZmOAS05OBgA4ODggICAAOp0O\n3bp1Q73fwTGRAAAgAElEQVR69TSuVBupqalwcnLChAkTMGXKFABlM4aNs0SJiIjIyNraGkFBQQgK\nCgIAXLp06Y7bp5GRkbC0tIRer4evr6/G1Za96OhoGAwGdOvWrUzbZQ8bERERFYuI4NChQ3juueeQ\nkZGBxMRE1K5dW+uyytSIESOwZs0aXLlyBRYWN/u9uNMBERERmQylFFq0aIHIyEikpaVhyJAhKCws\n1LqsMpGfn4/IyEisXbsWOp3OGNbKCgMbERERPRR3d3d8/fXXiImJQbNmzfD666/j559/Rnp6uvEY\nEcHWrVsRHByMZs2a4dVXX8WqVatw9epVDSt/eGlpaZg4cSLq16+PgQMHwsbGBmPGjCnzOnhLlIiI\niB6aiGDRokWIjIzEpk2bcOPGDZibm8Pb2xsdOnRAVFQU9u3bB3t7e3To0AE7duzA9evXjbsvBAUF\nQafToWPHjia9+8LQoUOxZMkSdO/eHaGhoejevftds2W5cC4RERGZvNzcXGzfvt24MO+ePXvQunVr\nhIWFYfDgwbC1tUVBQQESEhKg1+sRFRWFHTt2oKCgANWqVbtj94WWLVuazO4LhYWFqF27Nnr06IHv\nv/++yONMKrAppcwBJAA4JyK9lFKtAcwFYAvgNIAhIpJenHMf1BYDGxERUfmVk5MDKyur+wavjIwM\nbNy40Tj79K+7L7z44ov47LPPYGam7cit3bt3o127dli6dClefPHFIo8ztWU9RgM4DMD+1s8LAbwj\nIpuUUsMB/AvAB8U8l4iIiCqoqlWrPvAYOzs79O7dG7179wZwc/cFvV6PdevW4YsvvoCjoyPGjRtX\n2qXe14YNG6CUgk6n07QOoJiTDpRSdQH0xM2QdpsbgPhbz/UA+j/EuURERERG9erVw/Dhw7FmzRq8\n8MILmDBhAjZt2qRpTXq9Hh4eHnByctK0DqD4s0RnARgLwPCX1w4C6Hvr+QAARS15fK9ziYiIiO6i\nlML8+fPRqFEjDB48GJcuXdKsljp16uDQoUNITEzUrIbbHhjYlFK9AFwWkT1/e2s4gDCl1B4AdgDy\nHuLce7XzmlIqQSmVcOXKleJVT0RERBWOnZ0dfvjhB6SlpSE4OFiztd6+/vprODo6YsCAAbh+/bom\nNdxWnB42HwB9lFKnAawA4K+UWiIiR0QkSETaAFgO4GRxz71XIyIyX0S8RMTLFLoeiYiISDutWrXC\n119/jejoaHz44Yea1ODo6Ihvv/0WJ0+eRHh4uCY13PbAwCYi40Skroi4AHgBQKyIBCulagGAUsoM\nwATcnDFarHNL8gMQERFRxTR8+HC89NJLmDx5MmJiYsq07SNHjuCtt97CwIEDAQDXrl0r0/b/7nHm\nyw5WSh0DcATAeQDfAoBSylkp9WtJFEdERESVl1IKc+bMQdOmTfHiiy/iwoULpdpefn4+fvjhB/j7\n+6NZs2aIiIhA9+7dER8fj3Xr1pVq2w/ChXOJiIjIpB06dAht27ZF27ZtER0dXeL7eJ49exbz58/H\nwoULceHCBTRo0AAjR47E8OHDi7W5PTd/JyIiokqvefPmmDNnDjZt2oS3334bJTE50WAwQK/X49ln\nn4WLiwumT58OT09PrFu3DidPnsS4ceOKFdbKCnvYiIiIqFwYOXIk5s+fDwDw8PAw7kfaqVOnYi3W\nCwBXr17FokWLMHfuXJw4cQKOjo4YMWIERo4cCVdX10eqy6S2pipLDGxERET0dyKC3bt3G7ez2rZt\nG/Lz81G1alUMHDgQ33zzzT1vl4oIdu3ahTlz5mDlypXIycmBj48PQkND8fzzzz/25vMMbERERERF\nyMzMxKZNm/Dzzz9jwYIFGD9+/B1LgGRlZWH58uWYM2cO9u7dC1tbWwQHByM0NBTu7u4lVoep7SVK\nREREZDJsbW3Rs2dP9OzZEwaDAR999BE6d+4MFxcXzJkzB9999x2uX7+Oli1bIiIiAsHBwbCzs9O6\n7EfCwEZERETl3ldffYXdu3ejX79+yM3NhaWlJQYMGIDQ0FD4+PhAKaV1iY+FgY2IiIjKPWtra0RG\nRiIsLAyBgYEYPnw4atWqpXVZJYaBjYiIiCoENzc3REdHa11GqeA6bEREREQmjoGNiIiIyMQxsBER\nERGZOAY2IiIiIhPHwEZERERk4hjYiIiIiEwcAxsRERGRiWNgIyIiIjJxDGxEREREJo6BjYiIiMjE\nMbARERERmTgGNiIiIiITx8BGREREZOIY2IiIiIhMHAMbERERkYljYCMiIiIycQxsRERERCaOgY2I\niIjIxDGwEREREZk4BjYiIiIiE8fARkRERGTiGNiIiIiITBwDGxEREZGJY2AjIiIiMnEMbEREREQm\njoGNiIiIyMQxsBERERGZOAY2IiIiIhPHwEZERERk4hjYiIiIiEwcAxsRERGRiWNgIyIiIjJxDGxE\nREREJo6BjYiIiMjEMbARERERmTgGNiIiIiITx8BGREREZOIY2IiIiIhMHAMbERERkYljYCMiIiIy\ncQxsRERERCau3AW2nJwcbN++HZmZmVqXQkRERFQmyk1gO3XqFMaOHYu6deuiY8eOqFGjBnx9fTF9\n+nTs2rULhYWFWpdIREREVCqUiGhdw128vLwkISEBhYWF+O233xAREYH169fDzMwMffv2xfPPP4/k\n5GRERUUhMTERAFCjRg34+/tDp9MhKCgIrq6uGn8KIiIiqgyUUntExKtU2zDFwObs7CxvvPEG5s+f\nj9OnT+PJJ59ESEgIXnvtNdStW/eOY69cuYKYmBhERUVBr9fj7NmzAICGDRsiKCgIOp0O/v7+qF69\nuhYfhYiIiCq4ShvYlFICAF26dEFYWBj69esHS0vLB54nIjh69KgxvG3cuBGZmZmwt7fHli1b0KpV\nq1KvnYiIiCqXShvYPD095ccff4SLi8tjXScvLw/bt2/HCy+8gCeeeAK7d++Gra1tyRRJREREhLIJ\nbCY56cDc3PyxwxoAVKlSBb6+vli2bBmOHTuG0NBQmGJAJSIiIrofkwxsJc3Pzw+TJ0/GkiVL8N//\n/lfrcoiIiIgeSqUIbAAwfvx4BAQE4I033sCBAwe0LoeIiIio2CpNYDM3N8fSpUthb2+PAQMGcOFd\nIiIiKjcqTWADgNq1a2PZsmU4evQowsLCOJ6NiIiIyoVKFdgAwN/fH5MmTcLixYuxaNEircshIiIi\neqBKF9gAYMKECfD390d4eDjHsxEREZHJq5SBjePZiIiIqDyplIENAJ588kmOZyMiIqJyodIGNoDj\n2YiIiKh8qNSBDeB4NiIiIjJ9lT6wcTwbERERmbpKH9iAm+PZli5d+tDj2QwGA6KiotCvXz8EBQUh\nOzu7lCt9OMeOHcOiRYtQUFCgdSlERET0GBjYbgkICMDEiROxePFiTJo0CSdPnizy2NTUVHzxxRdo\n0qQJunbtiq1bt0Kv12P06NFlWPG9FRQU4KeffoJOp0OTJk3wyiuvYOLEiVqXRURERI+Bge0vPvjg\nA/Ts2RPTpk1Do0aN8I9//AOjRo3C6tWrkZaWhoSEBAwfPhxPPfUUxowZY5xpevbsWbz33ntYsGAB\nli5dqkntFy9exPTp0+Hq6opnn30WR44cwbRp0zB06FDMmDEDv/32myZ1ERER0eNTprichZeXlyQk\nJGjStojg2LFj0Ov10Ov1iIuLQ0ZGBpRSEBHY2NggODgYoaGhaN26tfG8goICdOnSBUlJSQgJCUFQ\nUBCeeeYZ2NjYlGq9Z86cwdixY7F69Wrk5+dDp9MhNDQUvXv3hoWFBbKzs+Ht7Y1z584hKSkJdevW\nLdV6iIiIKhul1B4R8SrVNoob2JRS5gASAJwTkV5KqdYA5gKwBXAawBARSf/bOVUBxAOwAmAB4AcR\nmfSgtrQMbH+Xn5+PnTt3IjY2FjVr1kRwcDCqV69+z2PPnTuHkJAQxMbGIjc3F1WqVEHHjh0RFBQE\nnU4HT09PmJubl2h9wcHBWLNmDUaNGoVRo0bBzc3trmOOHj0KLy8veHh4IC4uDhYWFiVaAxERUWVm\naoHtbQBeAOxvBbbdAN4RkU1KqeEAXEXkg7+dowDYiEimUsoSwBYAo0Vkx/3aMqXA9iiys7OxZcsW\n6PV6REVFITk5GQDg4OCAwMBA6HQ66HQ6NGjQ4LHaMRgMqF27Nrp164bFixff99ilS5ciODgY48aN\nw0cfffRY7RIREdH/K4vAVqwxbEqpugB6Alj4l5fdcLP3DAD0APr//Ty56fY6GZa3HqZ3D7aEWVtb\nQ6fT4dNPP0VSUhIuXryIJUuWoHfv3tiyZQtCQkLg4uICNzc3vPPOO8jJyXmkdvbu3YuUlBR069bt\ngccOGTIEISEhHM9GRERUDhV30sEsAGMBGP7y2kEAfW89HwCg3r1OVEqZK6WSAFwGoBeRnUUc95pS\nKkEplXDlypVillU+1K5dG0OGDMGiRYtw9uxZHDhwADNnzkSjRo3w+eef45133nmk627YsAEAoNPp\ninX8f/7zH7Rq1QrBwcE4e/bsI7VJREREZe+BgU0p1QvAZRHZ87e3hgMIU0rtAWAHIO9e54tIoYh4\nAKgLoJ1SqmURx80XES8R8XJycnqoD1GeKKXQokULvPXWW/j1118xZswYzJ49G5GRkcW+Rl5eHlat\nWoWFCxeiTZs2qFWrVrHOs7a2RmRkJHJycjB48GCuz0ZERFROFKeHzQdAH6XUaQArAPgrpZaIyBER\nCRKRNgCWAyh64TIAInINQByAB9+/q0RmzJgBb29vjBgxAmvXrkVWVlaRx549exYTJ05EgwYNMGjQ\nICilHno8WpMmTTBv3jxs2bKF67MRERGVEw8MbCIyTkTqiogLgBcAxIpIsFKqFgAopcwATMDNGaN3\nUEo5KaWeuPXcGoAOwJESrL/cs7S0xIoVK2BjY4O+ffuiRo0a8PPzw0cffYTdu3ejoKAA0dHReO65\n5+Di4oLp06ejTZs2+PXXX3HixAkEBQU9dJtDhgzBq6++yvFsRERE5cRDrcOmlOqCmzNDeymlRgMI\nv/XWGgDjRESUUs4AFopID6WUO4DvAJjjZjhcJSJTH9ROeZ8l+iiKmllqZWWF3Nxc1KxZE6+++ipG\njhwJV1fXEmmvffv2uHDhAhITE7k+GxER0SMyqWU9ylJlDGx/d+nSJcTExGDbtm1o3749BgwYgKpV\nq5ZoG0eOHIGXlxc8PT25PhsREdEjYmCjUnd7fbaGDRsaF/j18/PDE088AeDmzg/btm1DREQEdu/e\njY4dO0Kn0yEwMBC1a9fWuHoiIiLtMbBRmVi8eDFWrVqFuLg4ZGVlwczMDO3atYOPjw+io6ORnJyM\n6tWro1OnTtixYweuXr0KAHB3dzeGvM6dO8Pa2lrjT0JERFT2GNioTOXl5WHHjh3GfVR3796NVq1a\nITw8HC+++CJsbGxgMBiQmJhoHGu3detW5OXlwcrKCp06dTIGuNatW8PMrLjL/BEREZVfDGykqdv7\nod7cYezesrKyEB8fbwx5Bw4cAAA4OTlh8ODB+OKLL0p8/1QiIiJTUhaBjaPMqUhWVlYPPMbGxgbd\nu3dH9+7dAQAXLlxAdHQ01q1bhy+//BJOTk6YMGFCaZdKRERUobGHjUqFiOCll17C8uXLERsbC19f\nX61LIiIiKhUms/k70cNSSmHu3Llo3LgxBg8ejMuXL2tdEhERUbnFwEalxtbWFqtWrUJaWhqCg4Nh\nMBi0LomIiKhcYmCjUuXu7o6vvvoKer3+ofc9JSIiops46YBK3YgRI7Bx40ZMmjQJtWvXxrPPPgtH\nR0etyyIiIio3OOmAykRmZiZ8fHywb98+KKXg6elpXLPNx8enWDNSiYiITBEnHVCFYWtriz179mD7\n9u2YMmUKbGxs8O9//xsBAQGoUaMGVq5cqXWJREREJouBjcqMhYUFvL298cEHHyA+Ph6pqamYM2cO\nsrOzcenSJa3LIyIiMlkMbJXQpUuXjPuBasnOzg63b8l369ZN42qIiIhMFwNbJWEwGBAVFYW+ffvC\n2dkZTk5O8PLywvjx4xEXF4fc3FxN6tqwYQNcXFzQuHFjTdonIiIqDxjYKri0tDR88cUXaNKkCbp2\n7Yrt27dj7NixmDJlCqytrfHZZ5/B398fDg4O6NGjB2bOnImDBw+iLCaj5ObmIiYmBt26dbvvfqVE\nRESVHZf1qKD27t2L2bNnY/ny5cjOzkbHjh0xZcoU9O/f3zgj84MPPkB6ejo2btwIvV6PqKgorF+/\nHgDg7OyMwMBA6HQ6BAYG4sknnyyx2i5duoSFCxdi3rx5yMzMRJ8+fUrs2kRERBURl/WoQHJycrBq\n1SpERERg586dqFatGoKDgxEaGgoPD49iXePPP/+EXq+HXq9HdHS0caybu7s7dDoddDodOnfujGrV\nqhnPOXjwIObMmYOtW7fC29sbOp0O/v7+eOKJJ4zHiAi2bt2K2bNnY/Xq1cjPz0dgYCBef/119O3b\nt2S/CCIiojJUFst6MLBVECtWrMDrr7+Oq1evokmTJggNDcWwYcPuCE0Py2AwICkpCVFRUdDr9diy\nZQvy8vJQpUoVdOrUCZ07d0ZcXBzi4+NhZWUFb29v7NmzB5mZmTAzM0O7du2g0+ng6OiIhQsXYv/+\n/ahevTpefvllhIaGokmTJiX4DRAREWmDgY2KrWHDhqhWrRr+85//wM/Pr1TGhN24cQObN282Brj9\n+/fD1dUVoaGhGD58OGrWrIn8/Hzs3LnTeMyuXbtgMBjg4eGB8PBwDB48GDY2NiVeGxERkVYY2KhY\nTpw4gcaNG+Prr79GeHh4mbV7/fp12NnZwcys6LkraWlpuHjxIpo2bcqJBUREVCGVRWDjpIMK4Lff\nfgNQ9muZVa9e/YHH1KhRAzVq1CiDaoiIiCouLutRAWzYsAENGzZEw4YNtS6FiIiISgEDWwWwa9cu\nDuAnIiKqwBjYKoDg4GD8+uuvWL16tdalEBERUSlgYKsAZsyYgXbt2mHEiBE4deqU1uUQERFRCWNg\nqwCqVKmClStXQikFb29vvPzyy1i2bBkuX76sdWlERERUAhjYKggXFxf8+uuv8PPzw7p16zBkyBDU\nrl0bHh4eGDt2LPR6PbKzs7Uuk4iIiB4B12GrgAoLC5GYmGjcYmrr1q3Iy8tD1apV0alTJwQFBUGn\n08Hd3f2+a6hpraCgAGvXrsW3336LKlWqGLfGetTZsCkpKfjuu+/g4uKC/v37l3C1RERUWXHhXCoR\nWVlZiI+PN27wfvDgQQCAk5MTAgMDjQHuqaee0rjSmy5cuIAFCxZg/vz5OHfuHOrVqwelFP78808A\ngKurqzG8BQQE3HedNxHBzp07ERERgVWrViE3NxcWFhbYvHkzvL29y+ojERFRBcbARqXi/PnziI6O\nNm4fdXusW7NmzYzhzdfXF7a2tve9joiU2O4FIoL4+HhERERgzZo1KCgoQNeuXREWFoaePXvCzMwM\nx48fN/YaxsbGIiMjA2ZmZvDy8jLW7e3tjSpVquDGjRtYtmwZIiIikJiYCDs7OwwdOhRDhgzBkCFD\njL2QDg4OJVI/ERFVXgxsVOoMBgP2799vDELx8fHIycmBpaUlOnbsaOzJatOmDczNzWEwGBAdHY2I\niAj8+uuvcHd3h06nQ1BQEDp27AgrKysAQF5eHtasWYO5c+ciJSUF/v7+0Ol06NKlC+zs7Iztp6en\nY/HixYiIiMChQ4dQo0YNvPLKKxg1ahQaN25cZN35+fnYtWuXsddw586dMBgMsLGxQYcOHZCQkIBr\n166hZcuWCAsLQ3BwsLHd3bt3w8fHB35+fvjwww/h6ekJc3Pze343mzZtgqurK1xcXEr2iyciogqD\ngY3KXE5ODrZs2WIMcImJiQBubjHVpUsX7N+/HydOnICTkxOee+45HDx4EDt27EBBQQGqVauGZ555\nBo0bN8bKlStx+fJlNGzYEI0aNUJ8fDyys7NhYWEBb29v6HQ6XLhwAYsXL0ZWVhbatGmD8PBwDBo0\nCNWqVXvouq9du4a4uDjo9Xps3rwZLVq0QHh4ODp16nTPXsC5c+ciLCwMIgIHBwcEBAQYe+lsbW3x\n7bffYs6cOTh16hRcXFyQmJiIJ5544rG/XyIiqngY2EhzV65cQUxMDPR6PeLi4vDUU08hLCwMzz33\nnLE3LSMjAxs3bjSGvGPHjqFXr14ICwuDTqeDmZkZcnNzsW3bNmOP2N69e2FlZYUXXngBYWFhaNu2\nbZl/tkuXLiE6OtpY9/nz5wEA5ubmKCwsROfOndG7d2+MHz8evXv3xurVq7mBPRER3YWBjcql/Px8\nWFpa3veY1NRUmJubF2sD+bIgIjh8+DD0ej0uXLiA4OBgtGzZEgDw+eef45133sHHH3+MMWPGwMLC\nosjrXL16FdWqVYO1tXVZlU5ERBpjYCMyASKCvn37Yt26dbC3tzeOx9PpdGjUqBFEBHq9HhEREfjl\nl19gaWmJzp07G8f2mfryKURE9HgY2IhMRG5uLtauXWu8fXr69GkAQIMGDWBhYYGTJ0/CyckJw4cP\nR15e3l3Lp9wOeKa0fAoREZUMBjYiEyQiOHnypHE8XmZmJl555RX079/fOK4P+P/lU26HvEuXLgG4\nuXzK7d634iyfQkREpo2BjaiCEBHj8ilRUVF3LJ+yYMECDBs2TOsSiYjoETGwEVVQOTk52Lp1K6ZM\nmYJdu3Zh586daN26tdZlERHRIyiLwMaR0EQaqFq1KgICAvDDDz/AwcEBAwcOREZGhtZlERGRiWJg\nI9JQrVq1sHz5cpw4cQJdunTBp59+isTERBgMhjuOy8/PR2RkJIYMGYLDhw9rVC0REWmFt0SJTMB3\n332Hf//73zhw4ACAmzNLAwMDodPpcPr0acyfPx8XL14EADRt2hS7d+/mZAUiIhPBMWxElcxfZ5ZG\nR0fj4sWLUEqhR48eCAsLg6WlJbp27Yrg4GB899133HmBiMgElEVgK3rJdiIqc87Ozhg6dCiGDh0K\nEcGhQ4dgZ2eH+vXrG4+ZNGkSJk+eDEdHR7zwwgto06bNPTevP3z4MNavXw8PDw/4+PjcseTIbSKC\nXbt24erVq+jcuTPs7OxK9fMREdGjYQ8bUTlTWFiIF198EatWrQIA1KhRw7j7gr+/P5KTkxEREYG4\nuDjjOdbW1vD19TWu/+bq6oply5YhIiICSUlJAAALCwt06NDBuMCvl5fXfbfhIiKim3hLlIiKdOXK\nFcTExBgX5j1z5ozxvQYNGmDUqFEYNGgQDh48aFz/7ciRIwAAMzMzGAwGuLu7IywsDI0aNTLeit27\ndy9EBNWrV8dbb72FyZMna/QJiYjKBwY2IioWEcHRo0cRFxeHevXqoXv37ve8TXrmzBlER0fj0KFD\n6NevHzp27HjXOLiUlBTExMRg6dKlWLduHZYtW4bBgweX1UchIip3GNiISDP5+fno0qUL9u3bhz17\n9sDNzU3rkoiITBIXziUizVhaWmLFihWoUqUKBg4ciOzsbK1LIiKqtBjYiKhI9erVw/fff4/k5GT8\n85//1LocIqJKi4GNiO6rZ8+eGDt2LObNm4cVK1ZoXQ4RUaXEOftE9EDTp0/Hli1bEBISgtzcXAQF\nBaFOnTpal0VEVGmwh42IHuj2eDZHR0e8/PLLcHZ2RqtWrfD2229j/fr1yMrK0rpEIqIKjbNEiajY\nDAYDkpOTjWu/bd68Gbm5uahSpQo6duyIoKAg6HQ6PP300zAz478Hiahy4LIeRGTSsrOzsXnzZmOA\nS05OBgDUrFkTAQEB6NatG4KDg2FpaalxpUREpYeBjYjKlUuXLhl3TNDr9Th//jzCw8Px9ddfa10a\nEVGpYWAjonJLRDBmzBjMnDkTkZGReP7557UuiYioVHDhXCIqt5RS+OSTT+Dt7Y0RI0bg5MmTWpdE\nRFRuMbARUakoLCzEhg0bYG5ujvT0dMydO1frkoiIyi2uw0ZEJerKlSv45ptvMG/ePJw+fRp16tTB\n5MmT8cYbb2hdGhFRucXARkSPTUSwY8cOzJ49G5GRkcjLy0OXLl3w2WefoW/fvpwlSkT0mBjYiOiR\nZWVlYdmyZYiIiEBSUhLs7e0xcuRIjBo1Cs2bN9e6PCKiCoOBjYge2tGjRzFnzhwsWrQI169fR6tW\nrTBnzhwEBwfD1tZW6/KIiCocBjYiKpaCggKsXbsWERERiImJgaWlJfr374/w8HD4+PhAKaV1iURE\nFRYDGxHd14ULF7Bw4ULMmzcP586dQ7169fDhhx9ixIgRqF27ttblERFVCgxsRFSkMWPG4Msvv0RB\nQQGCgoIwe/Zs9OrVC+bm5lqXRkRUqTCwEdE9nT9/Hl988QWee+45fPzxx2jcuLHWJRERVVpcOJeI\n7mnDhg0AgIkTJzKsERFpjIGNiO5pw4YNePLJJ+Hu7q51KURElV6xA5tSylwplaiU+uXWz62VUtuV\nUvuVUuuUUvb3OKeeUipOKXVIKXVQKTW6JIsnotJRWFgIvV6Prl27cvYnEZEJeJgettEADv/l54UA\n3hORVgB+BPCve5xTAGCMiDQH4A0gXCnF1TSJTNzvv/+O1NRUNGrUSOtSiIgIxQxsSqm6AHriZki7\nzQ1A/K3negD9/36eiFwQkb23nmfgZuB76nEKJqLS5+rqCl9fX8yYMQOHDx9+8AlERFSqitvDNgvA\nWACGv7x2EEDfW88HAKh3vwsopVwAeALYWcT7rymlEpRSCVeuXClmWURUGszNzbFs2TLY2NhgwIAB\nuHHjhtYlERFVag8MbEqpXgAui8iev701HECYUmoPADsAefe5hi2A1QDeEpH0ex0jIvNFxEtEvJyc\nnIr9AYiodDg7O2Pp0v9r777jc7zeB45/ThIkYkbVnlU1S6wgNklQtCgVjVGUyret0tbot4P+Sumi\nVGJWtbQRq2pmETMpiaCUKqVij4iVhCTP+f2R5PkaQSTPSnK9X6/nlWfc59zXfWRczn3GUv7880/e\nfMtIXWsAACAASURBVPNNa4cjhBD5WlZ62NyBHkqpk0AA0EEptURrfURr7am1bgz8AhzPrLBSqgBp\nydpSrfUqE8UthLAADw8PPvzwQxYtWsSPP/5otvPEx8cTHx9vtvpzoxMnTjB+/Hhq1apFr1698Pf3\n5/jxTH/NCiHyAaW1zvrBSrUD3tNad1NKPa21vqiUsgN+AMK11t/fd7wCFgNxWut3snqeJk2a6Kio\nqCzHJYQwn9TUVDp16sTu3buJioqidu3aJqlXa01kZCR+fn4EBgaSkpJCkyZN8PT0xMPDg+bNm1Ow\nYEGTnCu3MBgMBAUF4efnx/r161FK0b59e/7++29OnToFpI0v9PDw4IUXXqB79+4yi1cIG6CUitZa\nNzHnOXKyDpu3UuoocAQ4CywCUEqVV0ptSD/GHRhAWq/cvvRH1xxFLISwKHt7e5YuXYqzszMdO3bk\n7bffZt26ddy4cSNb9d26dYsFCxbQuHFjWrZsyZo1axg+fDgffvgh9vb2fP7557Rt2xYXFxe6detG\nRESEia/I9ly5coUvv/ySZ599lq5du7Jnzx7++9//cvLkSUJDQzl58iR//fUX3333HfXr1+eXX37h\nxRdfZP/+/dYOXQhhKVprm3s0btxYCyFsS0REhPby8tJOTk4a0A4ODrp169b6//7v//SRI0ceW/7v\nv//Wo0aN0sWLF9eArl+/vp4zZ46+cePGPcfFx8fr1atX62HDhmlAv/XWW+a6JKsyGAz6999/14MG\nDdKFChXSgG7Tpo0OCAjQt2/ffmTZUaNGaUdHR52QkGChaIUQjwJEaTPnRk90S9RS5JaoELYrKSmJ\nXbt2ERwcTEhICHv37sXZ2ZmoqChq1aqVaRmtNWXLluXSpUv069cPX19f3N3dH3k7b9u2bbRt25ZV\nq1bRs2dPc12OxSUkJLBs2TJmz55NdHQ0RYoUYcCAAfj6+lKvXr0s1VG7dm0qV65s3D5MCGFdtn5L\nVAiRDzk6OtKhQwemTp1KdHQ0//zzD05OTo9c/kMpxeDBg9Fa06NHD1q1avXYsVebNm3C3t6eDh06\nmOMyLO7vv//m3XffpWLFigwZMoTExERmz57NmTNn8PPzy3KydurUKY4cOYKXl5eZIxZC2BJJ2IQQ\nOVKtWjWWLFnCwYMHefvttx963GeffUbLli15/fXXOXr06GPrDQoKomXLlhQvXtyU4VpUSkoKa9as\nwcvLi5o1azJz5kw8PT0JDw/n4MGD+Pr6UqzYA7v6PVJwcDAAnp6e5ghZCGGjJGETQuSYl5cXH3zw\nAQsXLmTJkiWZHlOgQAECAgIoWLAgffv2JSkpKdPjMnqi9u7dm2t7kVJTU5k2bRrVq1fnpZde4tCh\nQ3z66afExsYSEBBA27Ztsz2786mnngJg9erVpgxZCGHrzD1ILjsPmXQgRO6TnJysW7durZ2dnfXh\nw4cfetz69es1oN944w3jeykpKXrNmjXa09PTOKGhb9++Oi4uzhKhm9zBgwc1oMuVK6dXrVqlk5OT\nTVa3wWDQPj4+2s7OTm/ZssVk9Qohsg+ZdCCEyE3OnDlDw4YNKVeuHL///jtOTk6ZHjdu3Di++OIL\nZs6cyY0bN5g7dy6nTp2iQoUKjBgxgmHDhlGuXDkLR29aw4YN4/vvv2fjxo0m7ym8efMmTZo04dq1\na+zbt48yZcqYtH4hxJOxxKQDSdiEECa1adMmunTpwrBhw5g/f36mxyQnJ9OuXTt27doFQMeOHfnP\nf/5D9+7dcXBwsGS4ZpOQkICbmxvnz59n3759VKhQwaT1//HHHzRr1oxWrVoZJ2gIIaxDZokKIXKd\nzp07M2HCBBYsWMDSpUszPaZAgQKsX7+eqKgoDh8+TGhoKD179swzyRpA4cKFWb58OYmJiXh7e5OS\nkmLS+uvXr893331HaGgon3/+uUnrFkLYHulhE0KYXEpKCu3btycmJobo6Giee+45a4dkNUuWLGHA\ngAH897//5bPPPjNp3VprBg4cyM8//8z8+fPp2bMnJUuWzPTY+Ph4du7ciaurK+XLlzdpHELkd3JL\nVAiRa50+fdqYHERGRj50PFt+YO7xbO7u7hw4cAA7OzuaNm2Kh4eHcT/WP//8Ez8/P5YuXWpcJ69u\n3bp4eHjg6elJmzZtcHZ2NmlMQuQ3krAJIXK1jPFsw4cPZ+7cudYOx2oSEhKoX78+Li4u7Nmzx+T1\nJycns3v3buPuE7t37yY1NRVHR0eSkpJwcnKif//+9OnThwMHDhASEsK2bdu4ffs2BQsWpEOHDvz4\n44+ULl3a5LEJkR9IwiaEyPU++OADPv/8c5YuXUr//v2tHY5VGAwGypcvT7t27QgICDD7+a5du8aW\nLVvYsmUL1apVY9CgQQ/cKk1MTGTHjh0EBwcza9Ys2rVrx4YNG7Czk6HNQjwpSdiEELmejGeDffv2\n4erqyqJFixg8eLC1w3nAnDlzGDlyJFOmTGHChAnWDkeIXEdmiQohcj0HBwcCAgJwcnKib9++JCYm\nWjski8vYpN1Wt5MaMWIE/fr148MPP+Tzzz9n//792OJ/5oXIzyRhE0KYXYUKFfjpp584cOAA77zz\njrXDsahz586xfPlynn/+eZudnamUYu7cubRs2ZIPPvjAuPixj48Pixcv5uzZs9YOUYh8TxI2IYRF\ndO7cmfHjxzNv3jx+/vlna4djVlprtm3bxiuvvELlypWJjo5m+PDh1g7rkYoVK8b27ds5ffo0P/zw\nAx07diQkJITBgwdToUIFpk+fbu0QhcjXZAybEMJiMsaz7du3j6ioqDw3nu369essWbIEPz8/Dh06\nRIkSJRgyZAhvvPEGzz77rLXDe2IGg4E//viDjz76iA0bNhAeHk6rVq2sHZYQNkcmHQgh8pzTp0/T\nsGFDKlSokGfWZ0tNTWXcuHHMnTuXmzdv0qhRI3x9ffH29qZw4cLWDi/Hrl+/TqNGjUhKSmLfvn08\n9dRT1g5JCJsikw6EEHlOxYoVjePZxo4da7U4vLygcOG0rzl1584d5s2bR2pqKjt27CAqKoqhQ4fm\niWQN0m6XLl++nEuXLtGiRQvGjh1LaGgoSUlJ1g5NiHxDEjYhhMV16dKFZs2aERkZaZXze3lBcDAk\nJqZ9zWnS5uTkxLx580hMTGTt2rUopUwTqA1xdXVlxYoVVKxYkRkzZuDh4UHJkiXx9PTkyy+/5PTp\n09YOUYg8TRI2IYTF3bx5k5iYGNq1a2exc97do7Z9+//eL8s5PgxpC+fP56j+fv368cYbbzBt2jTW\nr1+fw2htU/fu3dmyZQtxcXGsX7+eN954g7NnzzJ27FhatGjB5cuXrR2iEHmWJGxCCIsLDw8nOTnZ\n5Ptq3i8jSStV6t4etbuHzX3E/9FS72Bt009zfL7p06fToEEDBg4cSGxsbI7rs1VFihSha9euTJ8+\nnYMHD/L7779z8eJFBg4ciMFgsHZ4QuRJkrAJISwuKCgIJycns844vPu2Z1zcvZ8lJkKSnRMahS/+\n2GOg+2l/UOrebO4JOTo6snz5cu7cuUO/fv1ITk7O4VXkDs2aNWPGjBls3LiRqVOnyqK7QpiBJGxC\nCIvbvn07dnZ2vPfee6xZs4br16+b4RwP/8zJCWoX/Ief6c8t0iYG3KIwv9i/CidO5Oi8zz77LPPn\nz2fXrl28++67Zrk2W/TGG2/Qt29f/vvf/1KpUiVee+01fv75Zy5evGjt0ITIEyRhE0JY3Mcff0zb\ntm1ZtGgRL730Ei4uLrRq1YpJkyYRERFBSkpKjs/RuvWD79nZgYtLWo/biaRyXKMYjiSRiCOOJHE1\ntRheg8rm+Nz9+vVj5MiRzJo1CxcXF9zd3Zk4cSK7du0yybXZIqUUixcvZsGCBbi7u/Pbb7/x6quv\nUqZMGVxdXRk7diwhISH5cmsyIUxB1mETQljN7du3iYiIICQkhODgYKKjo9FaU7x4cdq3b4+Hhwee\nnp4888wz2Zp56eUFoaFw97AqO7v/vV5BL85TjgVqOMP0PMpyjpdZhacnpG//mW0Gg4Ht27cTEhJC\nSEgIUVFRGAwGihUrRp8+ffD396dAgQI5O4kNS01NJSYmhuDgYEJCQti5cyfJyck4OjrSqlUr47/t\n888/j52d9B2I3E0WzhVC5CtXrlxh8+bNxj/y//77LwBVq1blxx9/pHVm3WaPUbhw2pi1DErB3b/2\nlEp73J3UOTlBQkJ2ryJzcXFxbN68mfXr1/PDDz8wfvx4Pv/8c9OexIbdvHmTbdu2GRPYQ4cOAVC6\ndGk6deqEh4cHHh4eVKxY0VgmMTGRwMBAfvjhB4oVK2Y8pmbNmnly6RSRe0nCJoTIt7TWHDt2jIUL\nFzJt2jQWLFjA0KFDn7iejMkHGTw9077e3/N2N1P0sD3KiBEjmDdvHhs2bKBLly7mO5ENO3v2LKGh\nocYE7sKFCwDUrl0bDw8PHBwc+OGHH4iLi6NmzZokJydzIn18YeXKlY3JW8eOHR+584LWmmvXrlGi\nRAmLXJfInyyRsKG1trlH48aNtRBCaK31119/rQF96tSpbNfh6am1k5PWLi5pXzNep/W1pT2U+t9n\n5paQkKCff/55XapUKT1v3jx94sQJ85/UhhkMBr1//3791VdfaS8vL+3o6Kjt7e31yy+/rDdv3qwN\nBoPWWutjx45pf39/3bNnT128eHENaKWUbty4sR4/frwOCwvTSUlJWmutr169qr/99ltdq1YtDeha\ntWrpt956S//222/6+vXrxnOfOXNGT5w4UVeqVEl/9NFHVrl+kfsBUdrMuZH0sAkhbJqnpyenT5/m\nzz//zFE99/e0ZUw++N95zNurdr+jR4/SpUsX/vnnHyBtdmlGr1H79u0pXry45YKxMUlJSSQlJT2y\nVywlJYU9e/YYe+giIyNJSUnBycmJpk2bEhUVRUJCAm5ubnTp0oXIyEi2bt1KYmIiDg4OtGjRglKl\nSrFu3TrjRJDJkyfzwQcfWOoyRR4it0SFEPlaQkICLi4ujBw5kunTp+eorvvHsjk5pc0k3b497asl\nk7UMWmsOHz5sTDrCw8O5desW9vb2uLm5GRM4Nzc3HBwcLB9gLnLjxg3Cw8MJCQlhx44duLq64uvr\nS+PGjY3HJCUlsWvXLmN7nzlzhldffZWiRYsyceJEoqOjadSokRWvQuRWkrAJIfK1TZs20aVLFzZu\n3Ejnzp1zVFdmY9mskaQ9yp07d4yzZu+fWdq+fXs8PT3x8PCgRo0aMujehHx8fAgODub8+fMyY1Vk\niyRsQoh8bfTo0cyZM4e4uDiccrADQYaMfUSt1aP2pDJmlmYkcBmD7qtUqWJcFqNjx464uLhYOdLc\ny2AwULZsWTw9PVmyZIm1wxG5lCRsQoh8rXbt2lSuXJmg3JBdWcCxY8eMyVtYWBjXr19HKUXjxo2N\nt09btmxJoUKFrB1qrpGYmEixYsXo3Lkzv/32m/RcimyxRMImfb9CCJt06tQpjhw5kuNboXlJjRo1\nGDlyJKtWreLKlSvs2rWLiRMnUqhQIb744gs6dOiAi4vLPRuz2+J/ym2Jk5MT06ZNY926dcycOdPa\n4QjxUNLDJoSwSfPnz2f48OEcOnSIOnXqWOScly9fxsnJCWdnZ4ucz5SuX79uHHQfHBzM0aNHAShX\nrpyx961Tp06ULZvzrbfyGq01L774Ips2bWLHjh00a9bM2iGJXEZuiQoh8q3evXuze/duTp06ZZHb\nVFprateuTbFixdi+fXuuv6146tQp4+3T0NBQrly5AsDzzz9vHP/WunVrk4wNzAvi4uJwdXXFzs6O\nmJgYWWhXPBFJ2IQQ+VJKSgpPPfUUffr0Yf78+RY559GjR3nuuecAGDVqFDNmzLDIeS3BYDCwb98+\n45ZfO3bs4M6dOxQqVOiefT0bNGiQr2dJRkZG0rp1a7p3787KlStlPJvIMhnDJoTIl37//XeuXbuG\nl5eXxc6ZMbGhV69efPvtt6xevdpi5zY3Ozs7GjVqxPjx4wkLC+Pq1ats3LiR//znP1y4cIHx48fT\nqFEjypYti7e3N99//z2xsbHWDtvktNbs2bPnodfWvHlzpk2bxurVq5k1a5aFoxPi0aSHTQhhcz7+\n+GMmT57M5cuXKVmypEXO+cILL/D333/zxx9/0KpVK6Kjo2nSpImx96lFixYULFjQIrFY2rlz5wgN\nDSU4OJjQ0FDOnz8PQK1atYzX37ZtW4oWLWrlSLMnPj6exYsX4+fnZxzbl3FtHh4etGvXznhtd49n\n27lzJ02bNrVm6CKXkFuiQoh8qU2bNkRHR/POO+8Yl6owZ7IUHR1NmzZtGDJkCLNmzeLixYv4+/sb\ntzxKTU3F2dmZtm3b4uPjg7e3t9lisTatNQcPHjSOf7t/O6eMxXubNGmCvb29tcN9pH379uHn58fS\npUtJSEigRYsWDBs2jPj4eIKDg9m2bZvx2lq2bMn8+fOpWbOmjGcTT0w2fxdC5Es///yzdnd31/b2\n9hrQhQsX1l27dtXTp0/Xly5dMsk5EhIS9A8//KCbNWumAe3s7KyjoqIeOC4+Pl7/+uuv2tfXV1ev\nXl0D+tdffzVJDLlBYmKiDgsL0+PHj9eNGjXSSikN6BIlSuhevXppf39/ffz4cWuHaWQwGPQvv/yi\nW7RooQHt5OSkhw4dqqOjox84NuPavL29NaBXrlxp/CwiIkI7ODjonj17GjefF+JhsMDm71ZPzjJ7\nSMImhNBa62vXruk1a9bo//znP7pmzZoa0A0aNNCJiYnZrvP48eP6vffe0y4uLhrQtWrV0jNnztTx\n8fGPLZuUlKQbNWqkS5QooU+cOJHtGHKzS5cu6YCAAD1kyBBdqVIlDWhAV69eXY8YMUKvXLlSX7t2\nzWrxrV+/XgO6Zs2aesaMGfrq1auPLfPxxx9rOzs7feXKFa11WtIXHh5u/J4LDQ01d9gil5OETQgh\n7rJmzRoN6JEjRz5RuZSUFL1u3TrdpUsXrZTS9vb2unfv3josLOyJe0+OHTumixUrpps1a6aPHj36\nyPJnz57VBw4c0KmpqU90jtzCYDDoI0eO6FmzZunu3bvrIkWKaEBXrFhRX7hwwSoxjRw5Ujs7O+vb\nt29nuYybm5tu3ry5vnbtmv7uu+90nTp1NKBLliyp3333XZ2QkGDGiEVeIAmbEELc57333tOAXrJk\nyWOTrUuXLulp06bpatWqaUCXK1dOf/LJJ/r06dM5imHFihXGW4NVq1bVr7/+ug4MDNRXrlzRBoNB\nh4WF6d69extv6ZYuXVp7e3vrRYsW5fjctuzOnTt6w4YN2tHRUXt4eFg8UTUYDLpatWq6R48eWS5z\n5coVYy9hRsLZpEkTvWjRIknURJZJwiaEEPe5c+eObt68uTFZGj58uF6+fPk9t7MiIyP1gAEDdKFC\nhTSg27Vrp5cvX67v3LljsjiOHz+u/fz89EsvvaSLFSumAa2U0mXLltWALlWqlB47dqxetGiR9vHx\n0WXKlDEmBrVr19Zvv/22Xrdunb5x44bJYnqYO3fu6FmzZulz586Z/Vxaaz137lwN6M8++8wi58tw\n9OhRDejZs2c/cRlAv/baa3r37t1mjFDkVZZI2GSWqBAi14mPj+fnn38mODiYLVu2GDdBb9KkCamp\nqezdu5eiRYsycOBARo4cSd26dc0aT0pKCrt37yYkJIQ//viDHj160LdvXxwdHY3HaK35448/jFtH\nbdu2jaSkJAoUKEDLli2NS0w0btzYZLMv//rrL1q1asXly5cB6Nq1K+vXrzdJ3Y+itebVV18lICDg\nnoV5755ZmpKSwm+//cacOXO4evUqHTp0wNPTE3d393vaLavi4uKYMGEC8+bN4/jx41SvXj3LZU+c\nOEHx4sVxcXF54vMKAbKsh7XDEELkAhnJUsYq/klJSQwdOpQBAwbY9LphSUlJ7Ny50xh3TEwMACVL\nlqRjx47GBK5atWrGMkePHsXf359ff/31njXS6tata1yVX+u0xWFbtGiBwWAwlu3Tpw9Lly6lQIEC\nFrm+mzdvMm3aNDZu3MjevXvRWlOiRAk6dOhAjRo1WLp0KWfOnKFSpUpUrVqVyMhIkpOTcXR0pE2b\nNsZrq1+//iN3HNizZw9+fn4EBASQlJTEK6+8QkBAgEWuUYgMkrAJIUQ+cenSJcLCwow9cKdPnwag\nRo0adOrUiWPHjhEaGkqBAgXw9PTk+PHjHDlyBEjb4L1Tp07Ur1+fZcuWER0dbax32bJl9O3b1yrX\nlOHy5cuEhYUZk9PY2Fg8PT3x9fXlhRdewMHBgZs3b7J161bj9R8+fBiAMmXK0KlTJzw9PenUqRPl\ny5cnMTGRgIAA/Pz8iIqKwtnZmQEDBuDr60v9+vWteq0if5KETQgh8iGtNX/99ZcxwQkPD6dkyZKM\nGDGCoUOHUrZsWQBiY2ONC9yGhIRw5coV6tatS+/evfn000+ZPHkyH3zwgcnju3PnDqtXr8bPz49z\n586xfft2ypQpk+Vru3XrFkWKFHnkcadPnyY0NNR4bZcuXQKgTp06nDt3jqtXr1KnTh1GjhzJwIED\nKVasWI6vS4jskoRNCCEEqamp2NnZPfLWoMFg4MyZM1SsWJElS5YwcOBAoqKiaNy4scniMBgMTJ06\nlZkzZ3LhwgWqV6/O2bNnadWqFZs2bTLbzgcGg4EDBw4QEhJCWFgYJUuW5I033qBNmzayQbuwCZKw\nCSGEeGI+Pj4EBwdz/vx57OzsTFbvrVu3qFatmrG3q3Xr1sTGxnLy5Ek++eQTJk6caLJzCZGbWCJh\nM91PshBCCJtw8OBBnn766XsmHZiCs7Mzp06dIiQkhLFjx3Lz5k1OnjwJwKRJk0hKSjLp+YQQ/yMJ\nmxBC5DFjxozh0KFDTJo0yeR1Ozo60qlTJ6ZNm8bevXu5cOEC8+bNY82aNdlajkMIkTVyS1QIIfKg\noUOHsmjRIoKCgvDw8LB2OELkaXJLVAghRLbMmjWLOnXq8Oqrr3L27Flrh5NrnDhxgtu3b1s7DCEe\nIAmbEELkQYULF2b58uXcunWL/v37k5KSYu2QbNadO3cICAigTZs2VK9enT59+mCLd59E/iYJmxBC\n5FG1a9fG39+frVu38umnn1o7HJsUEhJCpUqV8Pb25uzZs/Tv35+1a9cyffp0a4cmxD0kYRNCiDxs\n4MCBDB48mM8++4xjx45ZOxybc/78eS5evIiPjw9Hjx5lyZIl9OrVi3HjxhEREWHt8IQwkoRNCCHy\nuFq1aqG1xtnZ2dqh2BwfHx969epFQEAAv//+O0opFi5cSKVKlWjdujXu7u5MmjSJXbt2yW1lYVUy\nS1QIIfK49u3bExcXx/79+60dik2Kj4+nUaNGpKSkcPjwYZydnTl58iTz588nODiY6OhotNYUK1aM\n9u3bM3z4cLp27WrtsIUNkVmiQgghcuTmzZvs3LmTzp07WzsUm1WiRAnefvttYmNjjTNqq1atyuTJ\nk9mzZw+XLl0iMDCQV155hejoaF588UUiIyOtHLXIbyRhE0KIPGzLli0kJyfj5eVl7VBs2o4dO6hQ\noQI1atR44LNSpUrRp08f5s2bx4EDB6hYsSL9+vUjLi7OCpGK/EoSNiGEyMNiY2MBZNuoR0hJSSE0\nNBQvL6/HbiZfsmRJli1bxtmzZxk4cCD//vuvhaIU+Z0kbEIIkYcNHjyYunXrMnjwYFlA9yF2797N\ntWvXstwL2axZM77++mvWr19P1apVee6553jzzTdZs2YN169fN3O0Ir+ShE0IIfKwwoULExgYKAvo\nPsLOnTsBcHV1zXKZt956iz///JPp06dTo0YNFi1axEsvvYSLiwuenp5yu1SYnCRsQgiRx9WpU8e4\ngK45NoTPqTFjxtCrVy/8/f05fvy4xc/v4eFBwYIFGT16NAaDIcvlateuzTvvvMP69euJi4tjy5Yt\njB07lvDwcF577TXZLUGYlCRsQgiRDwwcOJDXXnuNyZMnExISYu1w7nH9+nVWr16Nr68vNWrUoHr1\n6owYMYIVK1ZYZOxdw4YNjbc4v/nmm2zVUahQIdq1a8eUKVP44osv+O2335gxY4aJIxX5mazDJoQQ\n+URCQgLNmjXj4sWL7Nu3j/Lly1s7JABu3LhBkyZNOHnyJBMmTGDfvn1s2bKF69ev4+XlxYYNG7Cz\nM2//gtaaPn368Ouvv9KnTx88PDzw8PCgUqVK2aqrZ8+erF+/nh07duDm5maGiIUtkXXYhBBCmEzG\neLa4uDi++uora4djVLRoUQIDA1FKERERwapVq7h8+TLTp08nKCiIadOmmT2GjB0OBg0aRHh4OEOH\nDqVy5crUrl2bt99+m02bNmX5FqdSikWLFlGhQgVeeeUVrl69auboRX6Q5YRNKWWvlIpRSq1Lf91A\nKRWhlPpDKbVWKVXsIeW+V0pdVEodNFXQQgghssfJyYnU1FSqVKli7VDu0aBBA7799luCg4OZOnUq\nBQoUYNSoUfTr148PP/yQ7du3mz2G4sWLs3DhQs6ePcuBAwf4+uuvqVq1KgsWLKBLly58/fXXWa6r\nZMmSBAYGcvbsWRnPJkxDa52lBzAG+BlYl/56D9A2/fkQ4P8eUq4N0Ag4mNVzNW7cWAshhDA9f39/\nDejDhw9bO5QHGAwG3a9fP21nZ6e3bt2qtdb62rVr+tlnn9Xly5fXFy9etEpciYmJulevXtrBwUFH\nREQ8Udnp06drQM+aNctM0QlbAETpLOY42X1kqYdNKVUReAFYcNfbNYFt6c9DgN4PSQi3ATK/WQgh\nbEBQUBCVK1fmueees3YoD1BKMXfuXJ555hm8vb356aefSEhIIDAwkCtXrjBgwIAnmsVpKo6OjsYN\n4V955ZUnWrJj1KhRlCtXjtDQUDNGKPKDrN4SnQGMBe7+STkEvJj+vA/w5CMzhRBCWExycjJhYWFZ\nWtHfWooVK8by5cuBtJmt5cqVY+DAgTg4OBAUFMQnn3xilbhKlChBYGAg586dY/DgwVm+xXn6edNV\naQAAHmdJREFU9GnOnTtHmzZtzByhyOsem7AppboBF7XW0fd9NATwVUpFA0WBOzkJRCk1XCkVpZSK\nunTpUk6qEkIIkYmUlBRSU1P566+/bHoB3QYNGhAbG0t0dDRTp06ldOnSJCcnA/DZZ58xZcoUq8TV\npEkTvvrqK9auXcv06dOzVCYoKAhA9nIVOfbYZT2UUp8DA4AUwBEoBqzSWvvcdUxNYInWutlD6qhK\n2ti3elkJSpb1EEII8/jxxx8ZNGgQH330EZ9++qm1w8myhIQE1q9fT9++fQH4888/qV27tsXj0FrT\nu3dv1q5dy/bt22nevPkjj3/55ZeJjIwkNjbWZns1Rc7ZxLIeWusJWuuKWuuqQD9gs9baRyn1dHqQ\ndsCHwBxzBiqEECLnBg4cyODBg/nss88ICwuzdjhZVrhwYfr06cP3338PpO3ekJqaavE4lFJ8//33\nVKxY8bHj2VJSUmz+FrTIPXKyDpu3UuoocAQ4CywCUEqVV0ptyDhIKfULEAE8p5Q6rZQampOAhRBC\n5Mx3332Ho6MjS5cutXYoT+y1116jQYMGAPzyyy9WiSGr49l2795NfHy83A4VJvFECZvWOlxr3S39\n+bda65rpj/Hp01rRWp/VWne9q4y31rqc1rpAek/dQtNeghBCiCdx+vRpEhMTc+0K/DExMWzfvh0f\nH5/HH2wmTZs2fex4tqCgIOzs7PDw8LBwdCIvkp0OhBAin9m0aROQewfCK6Vo1aqVtcPgrbfeomfP\nnowbN47IyMgHPt+0aRNubm6ULFnSCtGJvEYSNiGEyGeCgoKoWbMmVatWtXYouVrGeLbM1me7cuUK\ne/bsoUOHDlaMUOQlkrAJIUQ+kpycTHh4OAULFmTVqlXEx8dbO6RcrUSJEixbtuyB8WwnTpxAa82U\nKVNo2rQpH3zwAVu2bOH27dtWjljkVpKwCSFEPuLg4EDPnj05efIkvXv3plSpUrRo0YKPP/6YnTt3\nyp6X2XD3eLZvvvkGSFuz7ffff+eTTz6hUKFCfPHFF3To0AEXFxd69eolG8KLJ/bYddisQdZhE0II\n80pOTiYyMpKQkBBCQkLYvXs3BoOBiRMnWm03gdzsceuzXb9+nfDwcIKDg5k3bx5du3Zl9erVstxH\nHmGJddgkYRNCCMHVq1d58803+eWXXwgJCaFjx47WDinXiY+Px9XVFa01e/fuxcXFJdPjvv32W955\n5x2mT5/OO++8Y+EohTlIwiaEEMJibt26RdOmTYmLi2Pfvn2ULVvW2iHlOnv27MHd3Z2mTZsyYsQI\nOnXqRPny5e85RmtNr169WLduHTt27Mi1y6uI/7GJnQ6EEELkD87OzgQGBnL9+nX69+9vlZ0Ecrum\nTZsyb948/v77bwYNGkSFChWoV68eo0ePZsOGDSQkJDywW4KMZxNZIQmbEEIIo3r16vHdd9+xZcsW\nPvvsM2uHkysNHjyY8+fPExMTwxdffEG5cuXw9/fnhRdeoEGDBly9epWSJUuybNkyzp49y2uvvSaT\nPcRjScImhBDiHq+99hoDBgxg0qRJbN682drh5Ep2dnY0bNiQ999/n5CQEK5evcry5cv5999/GTJk\nCFprmjVrxpdffsmaNWuYMWOGtUMWNk7GsAkhhHjAzZs3adq0KfHx8ezbt48yZcpYO6Q8YcaMGYwe\nPZoZM2YwatQo43i29evXs2PHDpo1a2btEEU2yKQDIYQQVnPw4EGaNWuGu7s7mzZtwt7e3toh5Xpa\na3r27MmGDRuMCdrVq1dp1KgRWmtiYmJkK6tcSCYdCCGEsJp69eoxa9YsQkNDmTJlirXDyROUUixa\ntIjy5cvTt2/fB8azDRo0iLNnz1o7TGGDJGETQgjxUEOGDMHHx4eJEyeyZcsWa4eTJ5QsWZLAwMB7\nJhxkjGdbu3btPTNLN27cyK1bt6wdsrABcktUCCHEI928eZMmTZpw7do1Gc9mQhnj2Tw8POjWrRue\nnp7cvn2boKAgQkJC2L59O7dv36ZgwYJ06tSJJUuWyO1SGyW3RIUQQlhdkSJFWL58OfHx8fj4+Mj6\nbCYyatQoPv74Y/755x9GjRpF7dq16datG3/99RdDhw7l4sWLBAcH8/bbbxMcHGycXSryJ0nYhBBC\nPFb9+vX57rvvZDybCSmlmDRpEseOHeP48ePMmTMHNzc3Vq1ahbe3Nz/++CMeHh58+eWXfPHFF/z6\n66/MnDnT2mELK5GETQghRJYMGTKEV199lYkTJxIeHm7tcKzmxo0b+Pv74+rqyrhx40xSZ/Xq1Rkx\nYgQrVqxg8eLFANStW9f4+TvvvEOPHj14//332bNnj0nOKXIXSdiEEEJkiVKKOXPm8Oyzz+Lt7c2F\nCxesHZJF3bx5kzfffJPy5cvj6+vLvn37uHPnjsnPExQUhLOzM+7u7sb3MmaXlitXjldeeYX4+HiT\nn1fYNknYhBBCZFmRIkUIDAwkPj6eAQMGYDAYrB2SxURERDB79mwqV67Mp59+CoCXl5fJz7Np0yY6\ndOhAwYIF73nfxcWFZcuWERsbK+PZ8iFJ2IQQQjyR559/nlmzZhESEsLnn39u7XAsplOnTvTq1Yuj\nR4+yYcMGChUqRJs2bUx6juPHj/PPP//g6emZ6efNmzdn6tSprF69mvfff58jR45I4pZPSMImhBDi\niQ0dOpT+/fvz8ccfs3XrVmuHYxFKKRYuXEilSpWIjIykfv36FC5c2KTncHBwwMHBgbCwsIcmYmPG\njOGVV17h66+/pnbt2lSpUoWhQ4cSEBDA5cuXTRqPsB2SsAkhhHhiGePZatSogbe3NxcvXrR2SBZR\nokQJli1bRrly5fj4449NXn+VKlUeOyNUKUVAQIBxZmnTpk2NM0uffvppuV2aR8nCuUIIIbJt//79\nuLm50bZtWzZu3IidnfQD5JTWmpdeeomNGzcSGRlJo0aNHlsmNTWVqKgofvjhB+bMmcPMmTN56623\nLBCtAFk4VwghhI1r0KABM2fOJDg4mKlTp1o7nDxBKcU333xDcnIyoaGhWSpjb2+Pm5sbfn5+dO/e\nnXfffRfp+MhbJGETQgiRI6+//jre3t589NFH+WY8m7lFR0cD3LO0R1Yopfjhhx8oV64cffv2leU/\n8hBJ2IQQQuSIUoq5c+dSo0YN+vfvn2/Gs5lTUFAQxYsXx83N7YnLuri4EBAQQGxsLMOGDZPxbHmE\nJGxCCCFyrGjRogQGBnLlypV8tz6bqWmtCQoKolOnTjg4OGSrjhYtWjB16lRWrlzJ7NmzTRyhsAZJ\n2IQQQphEgwYN+PbbbwkODmbatGnWDifXOn78OGfOnHni26H3GzNmDN27d2fMmDEyni0PkFmiQggh\nTEZrTf/+/QkMDGTLli0mX1g2P0hISKBOnToopdi7dy8lS5bMdl1xcXG4urpib2/P3r17KVGihAkj\nFRlklqgQQohcJWM8W/Xq1fH29ubSpUvWDinXKVy4MMuWLeP06dM5XlNNxrPlHZKwCSGEMKlixYoZ\nx7MNHDhQxrNlg5ubG1OnTn3kArpZdfd4tu+++85EEQpLk4RNCCGEybm6ujJjxgw2bdok49myKWMM\n2vvvv8+ePXtyXFe3bt147733ZDxbLiVj2IQQQpiF1pp+/fqxcuVKwsPDadWqlbVDynUyxqDZ2dkR\nExOTozFoV65cwdXVFQcHBxnPZmIyhk0IIUSupZRi/vz5VK1alX79+snG5NmQMQbNFOPZSpUqxbJl\ny4iNjWXo0KEyni2XkYRNCCGE2WSMZ7t06RKDBg2S8WzZkDEGbfXq1Tkeg9aiRQs+//xzVq1aJeuz\n5TKSsAkhhDCrRo0a8c0337Bhwwa++uora4eTK2WMZzPFHqEZ49lkv9HcRcawCSGEMDutNX379mX1\n6tVs3bo1x4vC5kd3r6k2ffp02rdvT7FixbJVl4xnMy0ZwyaEECJPUEqxYMECqlSpQo8ePRg+fDgr\nVqwgLi7O2qHlGi4uLixbtoyrV6/y0ksv4eLiQqtWrZg0aRK7du0iJSUly3VljGc7deqUrM+WS0jC\nJoQQwiKKFy/O2rVrad26NcuWLaNPnz6ULl0aNzc3PvzwQ/bv32/tEG1e8+bNuXDhAuHh4YwbN447\nd+4wadIk3N3dKVeuHNHR0VmuK2M828qVK/Hz8zNj1MIU5JaoEEIIi0tOTmb37t1MnjyZjRs3AlCx\nYkViY2OtHFnuExcXx+bNmxkzZswT3+I0GAx0796d0NBQIiIiaNSokZmjzZvklqgQQog8KTk5mQkT\nJrBx40aKFCmCr68vISEh1g4rV3JxceHll182LtnxJLc47ezsWLx4MU8//TR9+vTh2rVrZo5WZJck\nbEIIISwuODiY7du38+mnn3LmzBlmz55NrVq1rB1WrtaiRQumTJnCypUreemll1iwYAGnTp16bLmn\nnnqKgIAA/v33X15//XUZz2ajJGETQghhcUFBQRQpUoRx48Zle6ajeNC7775r3Mrq9ddfp0qVKjz3\n3HO89dZb/PbbbyQnJ2dazt3dnSlTprB8+XL8/f0tHLXIChnDJoQQwqK01jzzzDPUr1+fNWvWWDuc\nPElrzaFDhwgJCSEkJIStW7eSkJBAnz59WLZsGUqpB8rIeLbskzFsQggh8pxjx45x4sQJvLy8rB1K\nnqWUol69eowePZoNGzYQFxfHxIkTWb58+UNnhN49nq1v375cv37dwlGLR5GETQghhEVcvnyZL7/8\nEi8vL5RSdO7c2doh5RuFChXio48+omvXrowZM4a9e/dmelzGeLaTJ0/KeDYbIwmbEEKIbDl69Cjv\nvvsuffr0Yfbs2Rw9evSBP/Baa3bv3s2gQYOoWLEiY8eOpXLlyqxdu5bq1atbKfL8KaszQt3d3Zk8\neTKBgYEyns2GyBg2IYQQWZaamsq6deuYPXs2ISEhODg4UK5cOeP6aVWqVMHDwwMPDw9u3ryJn58f\n0dHRFClShAEDBuDr60u9evWsfBX5286dO2nbti29evV65Hi2bt26ERYWRmRkJK6urlaINPewxBg2\nSdiEEEI81oULF1i4cCFz587l1KlTVKxYkREjRjBs2DDKlCnD8ePHjQPcw8LCjOOf6tati6+vLz4+\nPjIb1IZMmzaN8ePHM3v2bHx9fTM95vLlyzRs2BBHR0f27t0r/36PIAmbEEIIq9Fas2vXLvz8/Fi+\nfDnJycl07NgRX19fevTogYODQ6blUlJS2LNnD5C2lVJmPTjCuu7uQXvUjNAdO3bQrl07evfuTUBA\ngPxbPoQkbEIIISwuOTmZRYsW4efnx/79+ylWrBiDBw9m5MiRsrhtHpLVHrSM3jh/f3/eeOMNC0eZ\nO8iyHkIIISzuyy+/ZMSIEWitmTt3LmfPnuXbb7+VZC2PyeqM0Pfff58uXbrwzjvvEBMTY+EoRQbp\nYRNCCHEPNzc3lFJERETILbB8YOrUqUyYMAE/Pz9GjhyZ6TGXLl3C1dWVwoUL89dff8n3xX2kh00I\nIYRFXblyhT179tClSxf5o5xPjB07ls6dOzN69OiH9qCVLl2aF198kdOnT5OSkmLhCAVIwiaEEOIu\noaGhaK1lF4J8xM7Ojh9//JFSpUo9coeD0NBQ2rVrR4ECBSwcoQBJ2IQQQtxl06ZN2Nvbc/v27Ydu\nFC7yntKlSxMQEMA///xD165dWbhwIadOnTJ+fvLkSY4ePSqJvBVlPidbCCFEvlS+fHm01rRr146i\nRYvSvn1740K4NWvWlNukeVjr1q2ZM2cOn3zyCcOGDQOgZs2aeHh4cOLECQBJ2KxIJh0IIYS4x9Wr\nV9m8eTMhISEEBwcb/1hXqlQJT09PPDw86NixI0899ZSVIxXmoLXmzz//NP77b926lYSEBFq1asW2\nbdskac+ErMMmhBDC6u7fxeDatWsopWjUqBEeHh4MHz6catWqWTtMYSa3b99m7969tGjRwtqh2CxJ\n2IQQQtiUlJQUoqKijAlcREQEFSpUICYmhpIlS1o7PJO7dOkSCxcu5MCBA/j7+1O8eHFrhyRskCRs\nQgghbNrvv/9Oq1at6NatG6tWrcozt8tiYmL45ptvCAwM5M6dOyilePnllx+6WbrI32QdNiGEEDbN\nzc2NadOm8euvvzJz5kxrh2MSV69exd3dnTVr1jB8+HAOHTrE1KlTWb58Of7+/tYOT+RTMktUCCFE\njowePZqtW7fy/vvv07JlS5o2bWrtkHIkNDSUxMREduzYgbu7OwC1atVi69atjB49mlu3btG5c2fq\n1av3QG+bwWAgJCSEU6dOMWzYMOmNEyYjt0SFEELkWFxcHK6urtjZ2RETE0OJEiWsHVK2DR06lJUr\nV3L58mUcHP7Xr3H58mW6du3Knj17AChbtiydOnXC09OTZs2asW7dOvz9/Tl+/DiAbJaej8gYNiGE\nELlGREQEbdq0oUePHqxYsSJX9i5pralUqRLNmzdnxYoVmR4TGxtLaGgowcHBhIaGcvnyZeNnrVq1\nwtfXl8WLFxMeHk5kZCQNGza0VPjCSiRhE0IIkat8/fXXvPfee8yaNYs333zT2uE8sUOHDlGvXj3m\nzZvH66+//tjjDQYD+/fvJyIiAnd3dxo0aACkzS5t2LAhhQsXZtmyZTRs2BA7Oxk2nlfJpAMhhBBP\n5Pvvv6dKlSr4+PiwePFizp49a9Hzjxkzhm7duvHuu+8SHR1t0XObQsGCBbG3tycsLIysdGjY2dnh\n6uqKr6+vMVmD/231FBsbS+PGjSlTpgz9+vXj+++/JzY21pyXYJNSUlJ47733+OKLL4iJicFgMFg7\npFxHetiEECIPad26NYcPH8bOzo5Lly4BULduXTw8PPD09KRNmzY4OzubNYYrV67g6upKgQIF2Lt3\nb65bu2zq1KlMmDDBJGPQLly4YFyzLjg4mPPnzwNp/yY///wzzz//vClCtnnh4eG0b9/e+Lp06dJ0\n6tTJuO1ZxYoVrRhdzlmihw2ttc09GjdurIUQQjyZ+Ph4bW9vrz/44AOdmpqqY2Ji9LRp03SnTp10\noUKFNKALFiyo27dvr6dMmaL37NmjU1NTzRLLzp07tb29vX755Ze1wWAwyznMJTU1VXfu3FkXKlRI\n792712T1GgwG/ccff+hvvvlGly1bVj/77LP6+vXrJqvflo0fP147ODjoI0eO6MWLF2sfHx9dpkwZ\nDWillF64cKG1Q8wRIEqbOTeyenKW2UMSNiGEeHKrVq3SgN66desDnyUkJOigoCD93nvv6QYNGmhA\nA7pUqVL6/fffN0tSNW3aNA3o2bNnm7xuc7t48aIuX768fvbZZ/W1a9dMXv/WrVu1nZ2d9vb2znUJ\nbXa4urrqNm3a3POewWDQBw4c0B06dNCOjo56//79Voou52wqYQPsgRhgXfrrBkAE8AewFij2kHKd\ngb+AY8D4rJxLEjYhhHhyw4cP10WLFtV37tx57LHnz5/XS5Ys0T179tSAnjVrlsnjSU1N1V27dtUF\nCxY0aU+VpWzbtk3b2dnpfv36mSWpmjx5sgb0nDlzTF63LTl//rwG9OTJkx/6edmyZfVzzz2nb9y4\nYeHoTMPWErYxwM93JWx7gLbpz4cA/5dJGXvgOFAdKAjsB+o87lySsAkhxJMxGAy6SpUq+sUXX3yi\ncqmpqfqFF17QBQsW1NHR0SaP69KlS7pixYr6mWeeMUtPlblNmTLFbElVamqq9vT01IUKFdL79u0z\nef224scff9SAjoqKeugxmzdv1nZ2dtrHxydX9jhaImHL0ixRpVRF4AVgwV1v1wS2pT8PAXpnUrQZ\ncExr/Y/W+g4QALyYlXMKIYTIuqNHj/Lvv//SuXPnJypnZ2fH4sWLefrpp3n55Zf56aefjAPjTeGp\np54iICCAkydP8vrrr2f8Zz7XGDduHF5eXowaNYr9+/ebtG47Ozt++uknXFxc6Nu3Lzdu3DBp/bYi\nKCiI0qVL4+rq+tBj2rdvzyeffMKSJUtYtGiRBaPLPbI0S1QptQL4HCgKvKe17qaU2gV8obX+VSk1\nBpiktS56X7mXgc5a62HprwcAblrrBxbnUUoNB4anv6wHHMzBdeVFTwGXH3tU/iPtkjlpl8xJuzxI\n2iRz0i6Zk3bJ3HP350Cm9ti9RJVS3YCLWutopVS7uz4aAsxUSn0E/AbcyUkgWut5wLz0c0Zpc0+P\nzWWkTTIn7ZI5aZfMSbs8SNokc9IumZN2yZxSyuxrkWVl83d3oIdSqivgCBRTSi3RWvsAngBKqZqk\n3TK93xmg0l2vK6a/J4QQQgghsuixY9i01hO01hW11lWBfsBmrbWPUuppAKWUHfAhMCeT4nuAZ5VS\n1ZRSBdPL/2ay6IUQQggh8oGcbE3lrZQ6ChwBzgKLAJRS5ZVSGwC01inAm0AQcBgI1FofykLd83IQ\nV14lbZI5aZfMSbtkTtrlQdImmZN2yZy0S+bM3i42uTWVEEIIIYT4H9n8XQghhBDCxknCJoQQQghh\n4yyWsCml+iilDimlDEqpJne976GUilZK/ZH+tUMmZX9TSmW6LptSqplSal/6Y79Sqqc5r8PUzNgu\njy1vy8zYLqWUUluUUjeVUt+Z8xpMzVxtkv75BKXUMaXUX0opL3Ndgzlkp12UUpvSf18cUkrNUUrZ\nZ1JvQaXUovTy++9b1sjmmbFdCiilFqeXP6yUmmCpazIFM7bLq3f9LdqXXn9DS11XTpirTdKPe14p\nFZF+3B9KKUdLXJMpmPF7papSKvGu75XMJm0+yNxbKWQ8gNrAc0A40OSu912B8unP6wFn7ivXi7Qt\nsQ4+pN7CgEP683LAxYzXueFhxnZ5ZHlbf5ixXZyBVsAbwHfWvk4baZM6pG0bVwioRtp2cvbWvl5z\ntgvpex8DClgJ9Muk3v8Ai9KfPw1EA3bWvl4baJf+QED688LASaCqta/X2u1y3znqA8etfa3WbhPS\nlg47ADRIf11KfrdogKoP+338qEdW1mEzCa31YQCl1P3vx9z18hDgpJQqpLW+rZQqQtoepsOBwIfU\nm3DXS0cgV82iMGO7PLS8KeM3FzO2yy1gh1KqhlkCNyNztQlp28UFpH9vnFBKHSNtW7kIE1+CWWSn\nXbTW19PfdyBtn+PMfm/UATan13VRKRUPNAF2m/YKzMOM7aIBZ6WUA+BE2qLp1zM5ziaZsV3u5k3a\nVoy5ghnbxBM4oLXen17fFVPHbk4W+l7JMlsbw9Yb2HtXUvF/wNdAwsOLgFLKTSl1CPgDeEOnLSeS\nl2SrXR5RPq/IabvkRdlpkwpA7F2vT6e/l5c88DOglAoirUf+BrAikzL7SVs03EEpVQ1ozL0LgecF\n2WmXFcAt4BxwCvhKax1ngVgtKTvtcrdXgF/MF55VZKdNagJaKRWklNqrlBprmVAtKrvfK9XSb4du\nVUq1zsqJTNrDppQKBcpm8tF/tdZrHlO2LjCN/+2e0BB4Rms9WilV9VFltda/A3WVUrWBxUqpjVrr\npGxcgllYq10yK29LrNkutkraJHOmbJcMWmuv9PE0S4EOQMh9Rb8n7ZZIFPAvsAtIzdYFmImV2qUZ\nae1QHigJbFdKhWqt/8neVZieldolo7wbkKC1tqn9sK3UJg6kDUFpStp/GsOUUtFa67DsXYXpWald\nzgGVtdZXlFKNgV+VUnXv6p3LlEkTNq11p+yUU0pVBFYDA7XWx9PfbgE0UUqdJC3Op5VS4Vrrdo84\n/2Gl1E3S7imbfV+vrLJWuzykvM2w9veLLbJSm9j8FnImbpe7601SSq0h7bZwyH2fpQCj76prF3A0\nO3GYizXahbQxbJu01snARaXUTtJuFdtMwmaldsnQDxvsXbNSm5wGtmmtL6fXtQFoBNhMwmal3y23\ngdvpz6OVUsdJ6418ZN5i9VuiSqkSwHpgvNZ6Z8b7Wmt/rXV5nbYlVivg6EOSkmrpYylQSlUBapE2\nCDZXM0G7ZFo+t8tpu+RFJmiT34B+SqlC6bf+niWXjNN6lIe1i1KqiFKqXPpzB9L2QT6SSfnCSinn\n9OceQIrW+k+LBG9GOW0X0m6Ddkg/zhlo/pDjchUTtEvGVo19yUXj1x7FBG0SBNRP/1lyANoC+f5n\nSClVWqXPHlVKVSftd+7j/8OjLTfboidp2fZt4AIQlP7+h6SNh9h31+Pp+8pW5a4ZFUAP4NP05wNI\nG/S3D9gLvGSpa7LxdnlseVt+mKtd0l+fBOKAm+nnqGPt67WBNvkvabND/wK6WPtazdkuQBnS9jk+\nABwEZvG/meZ3/wxVTW+Pw0AoUMXa12oj7VIEWE7a790/gfetfa220C7pr9sBkda+RhtrE5/075WD\nwBfWvlZbaBfSxr3dnbd0z0o8sjWVEEIIIYSNs/otUSGEEEII8WiSsAkhhBBC2DhJ2IQQQgghbJwk\nbEIIIYQQNk4SNiGEEEIIGycJmxBCCCGEjZOETQghhBDCxv0/OWOUrnwkrxEAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize = (10,10))\n", "def animate(i):\n", " ax.clear()\n", " for n in range(45):\n", " ax.plot(test3.variables['traj_lon'][i, n], \n", " test3.variables['traj_lat'][i, n], 'b.')\n", " viz_tools.plot_coastline(ax, Bathy, coords = 'map')\n", " viz_tools.set_aspect(ax)\n", " ax.plot(213, 418, 'r*')\n", " ax.plot(-123.89095676900132, 49.16340592936349, 'r*')\n", " ax.set_xlim(-124.3, -123.5)\n", " ax.set_ylim(49, 49.5)\n", " ax.set_title('Time = ' + str(i))\n", "interval = 0.25#in seconds\n", "ani40 = animation.FuncAnimation(fig,animate,\n", " frames=44,interval=interval*1e+3, repeat=False)\n", "ani40" ] }, { "cell_type": "code", "execution_count": 123, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "array([ 59. , 59.25, 59.5 , 59.75, 60. , 59. , 59.25, 59.5 ,\n", " 59.75, 60. , 59. , 59.25, 59.5 , 59.75, 60. , 59. ,\n", " 59.25, 59.5 , 59.75, 60. , 59. , 59.25, 59.5 , 59.75,\n", " 60. , 59. , 59.25, 59.5 , 59.75, 60. , 59. , 59.25,\n", " 59.5 , 59.75, 60. , 59. , 59.25, 59.5 , 59.75, 60. ,\n", " 59. , 59.25, 59.5 , 59.75, 60. ])" ] }, "execution_count": 123, "metadata": {}, "output_type": "execute_result" } ], "source": [ "test3.variables['final_t'][:]" ] }, { "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 }