{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import matplotlib.colorbar as colorbar\n", "import matplotlib.colors as mplcolours\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import xarray as xr\n", "\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "mesh = xr.open_dataset('/home/sallen/MEOPAR/grid/mesh_mask201702.nc')\n", "land = 1-mesh.tmask[0, 0]" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "imin, imax = 100, 400\n", "jmin, jmax = 100, 350" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "watercolour = 'lightskyblue'\n", "landcolour = 'papayawhip'\n", "waterland_cmap = mplcolours.LinearSegmentedColormap.from_list('mycmap', [(0, watercolour), (0.85, watercolour), \n", " (0.850001, landcolour), (1, landcolour)])" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "array(6, dtype='timedelta64[ns]')" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "results = xr.open_dataset('/data/sallen/results/MIDOSS/AKNS_dispersion19/Lagrangian_AKNS_crude_15jan18-22jan18_AKNS.nc')\n", "scale = 86400 * 1e9\n", "(results.Beaching_Time[:].max()-results.Beaching_Time[:].min())/scale" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAADxNJREFUeJzt3XuMXOV9xvHn2d0ZE8DEXu+aEGyySVQS0VQN7tZKBaTCrVJCE9JKVZWql6SNZCl/RCRqG6BIrfJHpUCq9CJVrdzQhio0KU2gpUikkBSSRqqNxsRcjB1CqRGuTb2+kNgFeW+//jFn7LU9l2Myl/3Nfj/SyDNn3h3/3tdnHr/7zjlnHBECAOQxMugCAADnh+AGgGQIbgBIhuAGgGQIbgBIhuAGgGQIbgBIhuAGgGQIbgBIZqwXLzoxMRFTU1O9eGkAGEo7d+48HBGTZdr2JLinpqZUq9V68dIAMJRsv1i2LUslAJAMwQ0AyRDcAJAMwQ0AyRDcAJBMqeC2vcb2V23vtb3H9s/0ujAAQHNlDwf8c0lfj4hfsV2VdGEPawIAtNExuG1fIum9kj4qSRExK2m2t2UBAFops1TyNkkzkv7O9ndtf8H2RWc3sr3Vds12bWZmpuuFAgDqygT3mKRNkv4qIq6W9H+Sbj27UURsi4jpiJienCx11iYA4HUoE9z7Je2PiB3F46+qHuQAgAHoGNwR8bKkl2y/o9j0c5Ke7WlVAICWyh5V8glJ9xRHlLwg6bd7VxIAoJ1SwR0RuyRN97gWAEAJnDkJAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQzFiZRrb3STouaUHSfERM97IoAEBrpYK7cH1EHO5ZJQCAUlgqAYBkygZ3SHrY9k7bW3tZEACgvbJLJddExAHb6yU9YntvRHx7aYMi0LdK0hVXXNHlMgEADaVm3BFxoPjzkKT7JW1u0mZbRExHxPTk5GR3qwQAnNIxuG1fZHt1476k90l6pteFAQCaK7NUcqmk+2032v9DRHy9p1UBAFrqGNwR8YKkn+xDLQCAEjgcEACSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSIbgBIBmCGwCSKR3ctkdtf9f2g70sCADQ3vnMuG+WtKdXhQAAyikV3LY3SPpFSV/obTkAgE7Kzrj/TNKnJS22amB7q+2a7drMzExXigMAnKtjcNv+gKRDEbGzXbuI2BYR0xExPTk52bUCAQBnKjPjvkbSTbb3SfqKpC22v9TTqgAALXUM7oi4LSI2RMSUpA9L+veI+I2eVwYAaIrjuAEgmbHzaRwRj0l6rCeVAABKYcYNAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQTMfgtn2B7cdtP2l7t+3P9KMwAEBzYyXanJS0JSJO2K5I+o7thyJie49rAwA00TG4IyIknSgeVopb9LIoAEBrpda4bY/a3iXpkKRHImJHb8vCcvD883dq+44bOra7ofac/uC5/X2oCIBUMrgjYiEi3i1pg6TNtt91dhvbW23XbNdmZma6XScGIvTqq/tU/6WrtX2vndRinyoCcJ5HlUTEK5Iek3TONCwitkXEdERMT05Odqk8DFKlOq6IOS0snGjZZm4x9Mr8gtZVRvtYGbCylTmqZNL2muL+GyT9vKS9vS4Mg1etjEuSZmePtGxzbG5ekrSuUuZzbgDdUObddpmku22Pqh7090bEg70tC8tBpQjuubljkqaatjnSCO4qwQ30S5mjSp6SdHUfasEyU62ukyTNzh1t2eYIM26g7zhzEi2dmnHPEtzAckJwo6VqtVjjbjfjnq0H9wRLJUDfENxoaXT0Qo2MXKC5Nh9ONmbca8cIbqBfCG60Va2Md1jjXtDasVGNjbiPVQErG8GNtirVcc11WCrhiBKgvwhutFWtjGu2w4eTfDAJ9BfBjbY6zrgJbqDvCG60Va2saz/jZqkE6DuCG21VKuNaXHxNCwuvnfPcYoSOzTPjBvqN4EZbp47lbjLrfmV+QQvByTdAvxHcaKtSWStJmps791juxsk341wZEOgrghttVdqcPckFpoDBILjRVrXN9Uq4TgkwGAQ32qpUWl8hsLFUQnAD/UVwo62xsdWyK+1n3CyVAH1FcKMt26pU1rZc4754dESrRtiNgH7iHYeOqi3OnjwyyzHcwCAQ3Oio0uJ6JUfmOGsSGASCGx1VK+PNj+PmOiXAQBDc6KhSbTHjnl0guIEBILjRUbUyroWFE1pcPHlqW0SwVAIMCMGNjiqnvu392KltxxcWNRfBjBsYAIIbHTU7e5KTb4DBIbjRUaVy7vVKOPkGGByCGx01Lu26dMZ9lOuUAANDcKOjxox76Uk4h08FN5d0BfqN4EZHlcoaSSNnLpXMslQCDArBjY7sEVUqa878cHJuXm8YsS4aZcYN9BvBjVKq1XXnfDg5zvo2MBAEN0qpVMbPORyQZRJgMAhulFKtjJ8z4+aIEmAwOga37Y22H7W9x/Zu2zf3ozAsL5WzLu1KcAODU+adNy/pdyPiCdurJe20/UhEPNvj2rCM1K8Q+IoiFmSP1i8wxVIJMBAd33kRcVDSweL+cdt7JF0uieBejh66VXr56dOP3/QT0vs/+7pe6jP/ulvPHvihJOnkySlNjP2yrrv2mOZHx/Xa4qImOs24z66lRD1L/86Gq958if7ogz/+uvowcF389wAazmuN2/aUpKsl7Wjy3FbbNdu1mZmZ7lSHZWNkpB7Ss7NH+XZ3YMAcEeUa2hdL+pakP46I+9q1nZ6ejlqt1oXysFzMzx/X3NwPtGrVpVrUmA7OzumS0RG9kfAGusL2zoiYLtO21LvOdkXS1yTd0ym0MZzGxlZrbGy1pPqvaRsvqA62IGAFK3NUiSXdJWlPRHy+9yUBANops8Z9jaTflLTF9q7idmOP6wIAtFDmqJLvSHIfagEAlMCZkwCQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMkQ3ACQDMENAMl0DG7bf2v7kO1n+lEQAKC9MjPuL0q6ocd1AABK6hjcEfFtSUf7UAsAoISxQReA1r7xN/+oVYfd9LkLKxdq4+ormj5XffNFWvPBt5+x7T/ufU6HXzpx6vHExot13a9eKUl69IvbdOjFF049t/4tb9P1H936o5bfU3c8fof2Ht176vE7x9+pWzbfck67Vv3O2GegoWsfTtreartmuzYzM9OtlwUAnMUR0bmRPSXpwYh4V5kXnZ6ejlqt9qNVBgAriO2dETFdpi2HAwJAMmUOB/yypP+U9A7b+21/rPdlAQBa6fjhZET8Wj8KAQCUw1IJACRDcANAMgQ3ACRDcANAMgQ3ACRT6gSc835Re0bSi11/4TNNSDrc479jOVvp/ZcYA4kxkIZnDN4SEZNlGvYkuPvBdq3sWUbDaKX3X2IMJMZAWpljwFIJACRDcANAMpmDe9ugCxiwld5/iTGQGANpBY5B2jVuAFipMs+4AWBFWpbB3ewLim2/2/Z227uKL2zYXGy37b+w/bztp2xvGlzl3WN7o+1Hbe+xvdv2zcX2cduP2P5+8efaYvtQjUOb/n/O9t6ij/fbXrPkZ24r+v89278wuOq7o9UYLHn+92yH7Yni8VDtA1L7MbD9ieLferftO5dsH6r9oKmIWHY3Se+VtEnSM0u2PSzp/cX9GyU9tuT+Q5Is6T2Sdgy6/i6NwWWSNhX3V0t6TtJVku6UdGux/VZJdwzjOLTp//skjRXb71jS/6skPSlplaS3SvovSaOD7kcvxqB4vFHSv6l+vsTEMO4DHfaD6yV9Q9Kq4rn1w7ofNLstyxl3NP+C4pB0SXH/jZIOFPc/JOnvo267pDW2L+tPpb0TEQcj4oni/nFJeyRdrnp/7y6a3S3pl4r7QzUOrfofEQ9HxHzRbLukDcX9D0n6SkScjIj/lvS8pM39rrub2uwDkvSnkj6t+vuiYaj2AantGHxc0mcj4mTx3KHiR4ZuP2hmWQZ3C5+U9DnbL0n6E0m3Fdsvl/TSknb7dXrnHgrFV8ddLWmHpEsj4qBU36klrS+aDe04nNX/pX5H9RmmNMT9l84cA9s3SfqfiHjyrGYrZgwkXSnpOts7bH/L9k8XzYZ6DBoyBffHJX0qIjZK+pSku4rtzb4GfWgOlbF9saSvSfpkRPywXdMm29KPQ6v+275d0rykexqbmvx4+v5LZ46B6n2+XdIfNmvaZNvQjUGxH4xJWqv6ktDvS7rXtjXEY7BUpuD+iKT7ivv/pNO//uxXfb2vYYNOL6OkZrui+s56T0Q0+v6/jV9/iz8bvyIO3Ti06L9sf0TSByT9ehQLmxrC/ktNx+Dtqq/dPml7n+r9fML2m7RyxkCq9/W+YlnocUmLql+zZCjH4GyZgvuApJ8t7m+R9P3i/gOSfqv4RP09kn7QWErIrJg93CVpT0R8fslTD6j+n5iKP/9lyfahGYdW/bd9g6RbJN0UEa8u+ZEHJH3Y9irbb5X0Y5Ie72fN3dZsDCLi6YhYHxFTETGlelBtioiXNWT7gNT2ffDPqueAbF8pqar6haaGbj9oatCfjja7SfqypIOS5lTfMT8m6VpJO1X/xHiHpJ8q2lrSX6r+6fHTkqYHXX+XxuBa1X/Fe0rSruJ2o6R1kr6p+n9c35Q0Pozj0Kb/z6u+htnY9tdLfub2ov/fU3EEUuZbqzE4q80+nT6qZKj2gQ77QVXSlyQ9I+kJSVuGdT9oduPMSQBIJtNSCQBABDcApENwA0AyBDcAJENwA0AyBDcAJENwA0AyBDcAJPP/erRdaI8bazsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "days = np.ma.masked_array((results.Beaching_Time[imin:imax, jmin:jmax]-results.Beaching_Time[:].min())/scale,\n", " results.Beaching_Time[imin:imax, jmin:jmax] == results.Beaching_Time[:].min())\n", "#plt.plot((results.Beaching_Time[imin:imax, jmin:jmax]-results.Beaching_Time[:].min())/scale);\n", "plt.plot(days);" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xl4VEW6+PHv23tiQtgi+6agIqCoYVEWERABGZk77o46gwKOo3NdEK8zjqPeER33nzCKOhcZXBkVFWQTZRVZwyrIjuwEwhay9V6/P7rTJhAgQJKT7ryf58mT09V1Om+R5k11nTpVYoxBKaVU4rJZHYBSSqmKpYleKaUSnCZ6pZRKcJrolVIqwWmiV0qpBKeJXimlEpwmeqWUSnCa6JVSKsFpoldKqQTnsDoAgLp165rmzZtbHYZSSsWVZcuWHTDGpJ+qXpVI9M2bNyczM9PqMJRSKq6IyPay1NOhG6WUSnCa6JVSKsFpoldKqQSniV4ppRKcJnqllEpwmuiVUirBVYnplUqpimfyDxOe/28oOIxcOgDbeR2tDklVEu3RK1VNhOf+C/auh5wszA/jMEf3Wx2SqiSa6JWqLvIOANE9om12yD9saTiq8miiV6qakDbXgt0FDjckpUF6C6tDUpVEx+iVSmAm6Mf8OB1TcATbxT2Req2gMAcaXIQ4XFaHpyqJJnqlElj4+/dg12oIBQlvy8T2m+eQus2tDktVMh26USqR7dsEoWD0gUDOXkvDUdbQRK9UImvUFuxOQEAEajW2OiJlAR26USqB2brcjUlvgSnIwdbqKsR9jtUhKQtoolcqgYnNjlzUw+owlMV06EYppRKcJnqllEpwmuiVUirBaaJXSqkEp4leKaUSnCZ6pZRKcJrolVIqwWmiV0qpBKeJXimlEpwmeqWUSnCa6JVSKsFpoldKqQR3ykQvIh4RWSIiq0RkrYg8Gy0XERkhIhtFZJ2I/Hex8pEisllEVovI5RXdCKWUUidWltUrfUBPY0yeiDiB+SIyDWgNNAEuMsaEReTcaP1+QKvoVydgdPS7UkopC5wy0RtjDJAXfeiMfhngfuAOY0w4Wm9/tM5A4P3oeYtEpKaINDDG6NY2SillgTKN0YuIXURWAvuBb40xi4HzgVtFJFNEpolIq2j1RsDOYqfvipYd+5pDo+dmZmdnn10rlFJKnVCZEr0xJmSMaQ80BjqKSFvADXiNMRnAv4D3otWltJco5TXfNcZkGGMy0tPTzyx6pZRSp3Ras26MMUeAOUBfIj31CdGnvgQuiR7vIjJ2X6QxsOesolRKKXXGyjLrJl1EakaPk4DewHrgK6BntNrVwMbo8STg7ujsm85Ajo7PK6WUdcoy66YBME5E7ET+MHxqjJksIvOBj0TkESIXawdH608F+gObgQJgUPmHrZRSqqzKMutmNXBZKeVHgOtLKTfAA+USnVJKqbOmd8YqFQ8K9kHuDojMZlbqtJRl6EYpZaXd82D394BASkNofTeI9tFU2em7Ramqbvd8CAcg7Ie83VB4wOqIVJzRRK9UVec855djY8CRZF0sKi5poleqqrvwdkiuB640OP+/wJVqdUQqzugYvVJVXfK5cMn9Vkeh4pj26JVSKsFpoldKqQSniV4ppRKcJnqlqrj9O7bw1iN38PrQX7Hmh++sDkfFIU30SlVxHz//KPu2beRI9l6+HPk0uYd1Hr06PTrrRqkq6Mj+vUx+5wW8+bnHJfbPXv0zoUCQ3nc+QIt2GRZFqOKJRNYgs1ZGRobJzMy0Ogylqox//vctHNj1M8YYbHYHNrsdm82OwRDwecEYnG4Pj7zzNeek1bI6XGUREVkW3fzppHToRqkq6PC+3RR1wuwOBzc+8hx3PvVG5MlouYiNowf3WRWiiiOa6JWqgtpfMwCnJwmXJ4naDZpwYYduNLv4Mq7o/Wtc0fK09PqkNznf6lBVHNChG6WqIGMMm5cvwFeYzwUduuNye2LlW1YtpjA3hws7dMfl0XVvqrOyDt3oxVilqiARodUVXUotb9m+swURqXimQzdKKZXgNNErpVSC00SvlFIJThO9UkolOE30SimV4DTRK6VUgtNEr5RSCU4TvVJKJThN9EpVUeFwmFAwaHUYKgFooleqCtq4bD7P39Gd527twuzx71odjopzmuiVqoK+eONpAj4v4XCI+V/8m5wDWVaHpOKYJnqlqiATDpd4HD7msVKnQxO9UlXQr+7/C3anC5vdQYe+N1Hr3IZWh6TimK5eqVQly9v7Mz99/gbhgI8LfjWU2q0uO65O2y7XckFGN0LBAEnnpFoQpUok2qNXqpKtHPsMeXu2UJC9i9UfjCDoKyi1nsvt0SSvyoUmeqUqkTGGQH5O8QKChXnWBaSqhVMmehHxiMgSEVklImtF5Nljnh8lInnFHrtF5D8isllEFotI8/IPW6n4JCI06twfm9ON3eWhZos2uNPSrQ5LJbiyjNH7gJ7GmDwRcQLzRWSaMWaRiGQANY+pfy9w2BjTUkRuA14Ebi3fsJWKX60GDKHepd0JB/zUbNEGEbE6JJXgTtmjNxFFPXZn9MuIiB14GXj8mFMGAuOix58DvUTfyUrFiAhpTS+i1vmXIDa71eGoaqBMY/QiYheRlcB+4FtjzGLgQWCSMWbvMdUbATsBjDFBIAeoU8prDhWRTBHJzM7OPps2KKWUOokyJXpjTMgY0x5oDHQUke7AzcCoUqqX1ns3pbzmu8aYDGNMRnq6jlEqpVRFOa1ZN8aYI8Ac4BqgJbBZRLYBySKyOVptF9AEQEQcQBpwqJziVSouBYNBZk75imlf/gdvYenTKZWqKKe8GCsi6UDAGHNERJKA3sCLxpj6xerkGWNaRh9OAn4HLARuAmYZY47r0StVnTzz8BAyF3wPGL74YAzvfD4dm01nN6vKUZZZNw2AcdGLrzbgU2PM5JPUHwN8EO3hHwJuO/swlYpf4XCYH2bNoKi/s33rZg7sz+Lc+rqsgaocp0z0xpjVwPH3aJesk1Ls2Etk/F4pBdhsNuo3bsq+3TsJh8O43G7SatW2OixVjehaN0pVoIN7dhAKBnjtvU8Z/dKz+Hw+hj76Z9xuj9WhqWpEE71SFWTmR6NZMOlDRIS2Xa7lf0f+n9UhqWpKrwYpVQHCoRDffzGWoN9HwOdl9bzp5OcctjosVU1poleqAojNhsud9MtjERwut4URqepME71SFUBEuP0vr1Kj9rmck1abGx8dgTsp+bRewxijO0upcqFj9EpVkBZtMxg2ZuoZnZv5w1yeemgwfp+Pe/40nN8O/VM5R6eqE+3RK1UF/X34gxTm5xMKBhn35mvsz9pjdUgqjmmiV6oSbP9pBa8NGcDLg65j7cKZp6zv93ljx+FwmFAwWJHhqQSniV6pCmaM4ePnHyXnQBZ5Rw7yxetP4SvMP+k5Dscvo6rhcAi3J+kktZU6OU30SlUCf7GFzIwx+L2FJ60fLNaDd7pc5OcerbDYVOLTRK9UBRMRut10Dw6XG6fbQ5su15Jaq+5Jz/n9g8NwuT14kpLIuOpqGjc/r5KiVYlIqsLCkhkZGSYzM9PqMJSqUAd2bycUDOBOPocZ/34Dv7eQ3nc+QP0WF5Raf8/O7RTk53HeBa11pUtVKhFZZozJOFU9nV6pVCWp26gZACMf+A2HsnZhwmF2rl/FY+9Nx1nK2jcNmzSr7BBVgtJuglKVyBgTS/IAwWCAvBzdl0dVLE30SlUiEeGijlfj9CThdHuo06AJaXXqWR2WSnA6dKNUJbv5sRf4acF3BHw+2nbtg81uP65OwO9n6hfjKcjLpe9/3UqtOie/eKvUyWiiV6qS2e0O2nXre9I6zzxyH5kL5hEKBfniw/f4aMYCXLoomjpDOnSjVBW09Ie5+LyFBAMB8nKPsmfHdqtDUnFME71SVVCr1m1xOJ2ICHaHg3MbNLI6JBXHdOhGqQpijGHRokX4/X66det2WnPh//H2+/z7zdfIO5rDHUMeJPmccyowUpXoNNErVUEeeOAB3n//fUSEXr168eWXXyIiZTo3Na0mf/rL/1ZwhKq60ESvVDnbtm0bfr+fd999l1AoBMD06dPZv38/9erpVEpV+XSMXqly9Mwzz9C6dWsuvfRSnE5nrNzhcJCammphZKo60x69UuUkGAwyYsSI2MqTDoeDiy++mFAoxFtvvUVy8ultJVgVBPx+vIUFpKbVtDoUdRY00StVTmw223EXXJcuXRqXCR5g7YpMhg+5A7/PR8du1/D3UWOwl3Jzl6r6dOhGqXJy7Gbe4XC4xLry8eb//f0vFOTnEQwGWLH4B1Ysmm91SOoMaaJX6ix9/vnnpKamHjcG73A4yjzLpiqy2Yr33g2iSyXHLf3NKXUWQqEQd911F3l5eRQWFmKMweVy4XQ6efXVV+P6AuywZ18kNa0mNpuNLj37clmnLlaHpM6QjtErdRbC4TCBQKBEWVZWFikpKSVm3cSjC9pcwqSFawkE/LrOTpzTHr1SZ8Fut5e4QCkiJCcnx32SLyIimuQTgCZ6pc7AypUr6dKlC1dddVWJC64iQmHhyTf+Vqqy6dCNUqcpHA7Tq1cvDh2K7AzlcrliF15vueUWatbUOeeqatFEr1QZbdmyhSeffJJgMEhOTk6sPBQKMWvWLJKTk2nfvr2FESpVulMmehHxAPMAd7T+58aYp0XkIyADCABLgPuMMQGJzCd7A+gPFAC/N8Ysr6gGKFUZwuEw3bp1Y9++fUCkF5+UlARAjx496NJFZ6SoqqssPXof0NMYkyciTmC+iEwDPgLujNb5GBgMjAb6Aa2iX52iZZ3KO3ClKlNubi7Z2dmxG6L8fj8ffvghHo+Hvn1PvluUUlY7ZaI3xhggL/rQGf0yxpipRXVEZAnQOPpwIPB+9LxFIlJTRBoYY/aWb+hKVZ7U1FTsdnvswqvL5eI3v/lNXN8QpaqPMo3Ri4gdWAa0BN40xiwu9pwTuAt4KFrUCNhZ7PRd0bISiV5EhgJDAZo2bXqG4StVsfLy8hg3bhw+n6/E7JpAIEB+fj4pKSkWRnd6wuEws6Z+xcHs/fTsP5D0eg2sDklVkjIlemNMCGgvIjWBL0WkrTFmTfTpt4B5xpjvo49L6+KYUl7zXeBdgIyMjOOeV8pqxhi6d+/O+vXrY3e8FiX7Vq1acU6c7fr05j+eZsrnnxAMBvj4X//kw+k/kFojzeqwVCU4rVk3xpgjIjIH6AusEZGngXTgvmLVdgFNij1uDOw5yziVqjTBQICFc74jvyCfNWvWxO58FRGeeOIJnE4nDz30UNwN28yZ/jXewgIg8olk008/cnnnrhZHpSpDWWbdpAOBaJJPAnoDL4rIYOA6oJcxJlzslEnAgyIynshF2Bwdn1fxwhjD8CF3sP7HlRhjcDsdhEIhRITGjRszYsSIuEvwRS5qdxlL5s8m4PcTDoVo0uJ8q0NSlaQsPfoGwLjoOL0N+NQYM1lEgsB2YGH0jf+FMeZ/galEplZuJjK9clCFRK5UBcg5fIg1y5cSCPgBaNOwBud36InD6eS5556L2yQP8ORLo3hv5Mvs27uL2wbdr2P01YhEJsdYKyMjw2RmZlodhlIE/H76XdEyNhbv9niYvnxLXCd4lbhEZJkxJuNU9XStG6WKKcjPKzFzwO/zEYrjzUOUAk30SpXgSU7G4fhl5cnklFTsDl0pRMU3TfRKFeN2exjx5lgaNG5C42YteOmdD3XYRsU9TfRKHePnTevZn7WXrD272L51k9XhKHXWNNErVUxBfh7vvDKCUDBIMBDgtWefIBQKWR2WUmdFE71SJ2HCYX7VqTUDOrXmh1nfWB2OUmdEE71SxSSfk8KQR/+M3eHA4XBiiPTy847m8Oyjf9DevYpLmuiVOsatg/7A1MyNTJi3okR5MBAgrIlexSGdN6ZUKdxuDzaxIcXW6LPZ7TrVUsUl7dErdQI+n5fiMyvDoRAP/+4mHv7dTWzduN66wJQ6TZrolTqBlNQa9Oh7A56kZDxJSdjtDlZnLmLlkgU88vubYrtNKVXV6edQVarNKxay9JsJ1GvWku433YvD6Tz1SQnoyZdG8dOq5fh8Xh6757ZYee7RHHxeL0nJyRZGp1TZaKJXx9m3fTPjXxxOwOdl84qFePNz6T94uNVhWUJEaNP+CgCuuLIra1ZEFt+7JKOTJnkVNzTRq+Nk/byBolVNg34f239aaXFEVcMLb3/AwtnfAtDu8g78+5+vYozhN3feQ1qt2hZHp9SJaaJXx3GlpePzRi5EGgNem8fqkKoEh8NBt2v7YYxh0A092bV9KxiYPX0S4ybP1TVxVJWliV7F7Nm5nYVzviV7317WZuWR4jAU+IOkhrcw4YP/o3Gz8+jUvafVYVqusKCAHT9vjs2p371jG3m5R0mtkYbP5+Pjjz8mGAxyxx13xN2+sioxaaJXAOzbs4vB/3UtgYAfERvhcIjDgQBut4ec7dt4+5UR2O127vnTcG4ZdN+pXzCBJSUn06hJM/bs2gFAvYaNSUmtAcD111/PwoULAXj77bfJzMzUnr6ynCZ6BcDyxT8QDocI+CNb6DU7ryWXdrwKv7eQOd9MxltYSACY/tWn1T7RiwijPp7I+DGjMcZw2z33IyIEAgFmzZoVu76xdu1aPvvsM5KTk+nbty8OvdlKWUTfeQqABg0a4vd5AbCJ0Kx5Mx59+h+sWrqQ6V99FqtXo2ZNq0KsUmrWqsMfHvtriTKn00mzZs3YsWNHLNkPGjQIm81G586dmTFjhvbulSX0hqlqzltYwJoVS8ndt4O2TetRK9lFo1rJNPBEEtLB7P243O5Y/aNHjlgValyYM2cOffr0oUePHvj9fgoKCsjLy+P7779n3759Voenqint0VdjuUdzGPKbPhw9cphkh3BJ0zrU8tTF7nBSr9l5ADRo3BSbRPoDDqeLZue3sjLkKm/kyJHMmzcPAI/HQ2FhIQAul4ua+mlIWUQTfYIJhULsz9pD7Tp1cThdJz1eOn8ORw4dxFtYQAFwVJrSvElj6jVrSf/BjwHQ+pLL+OP/PMOn/36H5q0u5NFnXrS2gVVYIBDgjTfeiC1l7HA46NSpE6FQiJEjR+Lx6DRVZQ1N9AkkPy+XP942gKzdu3A6naTVqs2B/Vk4nC5q1a5D9r692B1O6qSfy/49u4ms2BUZS7Y7HNRo3JIHXx1d4jX9fh/fTPyU/Vl7OJi9j+1bNtL2sg4WtK7qczgcOByOEmvWz549m6SkJAujUkoTfVwL+P04nM7IjA+/n2+/nkDWrp34fF583kIK8vMIh8P4vF4Ko8fgxVuQH1uQq37DxhQWFNDs/Fb895PPHfczFs+dxdaN62MXat/8x7OM/s/kymxmXAiFQgQCAYLBYKzMGIPP59NEryyniT4O+byFDLv3NtauyKR+oybUrF2H9T+uJDWtFkXLp4sIRWvsnujYZrPR6uJ2/H3UmBP+LJfHg4n2+kVE13cpxfvvv8+QIUMIh8PYbLZYj95ut+OspovBqapFZ93EoWlf/IdNa3/EGEPW7p1sWLsaYwy5OUdITauFzWajSYuWXNG5GzabjcbNz6PDVd2x2Ww0atqcjl2vwWaz0bBJcx544tmT/qwOXa6mx3W/wma3U69hYx59+h+V1Mr4EAwGGTJkCH6/n2AwSDgc5pxzziEpKYl33nlH74xVVYIUzfe1UkZGhsnMzLQ6jCon4Pfzv4/dz9L5c7ng4nak1qzFsgXzqF03nex9WQT8PhDBFr2TVWw2evYbyF9f/mdsvrYx5pTHZXG69auLQCBAUlJSiV784cOHSUlJ0X8vVeFEZJkxJuNU9XTopgqbOuETlnw/B5+3kLUrMwEhFAqyb+9uUlJqUIAhJbUGNWrWYs+O7aTUSOOe/x5eIsGU5bgsNGmVzmazlRiugciNU/rvpaoSTfRnwe/38crfhrNq6SK69OzDA088i91uP+3XCYVCjBrxFAvnfkf7Dp1JSk5h4dzvSE1LIxgMxOoUJY9wKESHrj148M/PklIjDbvdTs7hQ6TUSNPb7CtZKBSi+KdiEaFDhw6Ew2Heeustrr76agujUypCs8JZ+OT/3mTO9K/x+3xMmfAJLS9qQ/8bbz/t15k64ROmffkffN5CZk75CgOEgkEOHcjG5XIj7sj8a5fbg9/nRWw27hjyIDVr14m9RvFjVXlcLhfDhw/n9ddfR0QIhUKsWbMGgAEDBnDgwAHcxe4sVsoKmujPwu4d2/H7fAD4vF6mf/kpkz//mC49++D3+Vj6w1yuuqYPvx3y4HEf5Y0xfPSvf7Jg9gzsdnts+mLx6XkBv48e1w3g1nv+QIPGTXE4HOzctpX6jZrEVktU1nv++ecZPHgwPp+Ptm3bxsp9Ph833HADxhhefvllLr30UgujVNWZJvqzUKNGsVvajeGnVcsJBgNsXLMqtprh1g0/cW69BvQZeFOJc2dM/JwP3/5/eAsLcbncOJ0unC4XgWAAu82GzWYnFApx492DaXlRm9h5xY9V1XHeeZElIwYNGsT48eNjvfvvvvuOcDhMjx49yMrK0t69soQm+jIwxvDVJ/9m9bIlXN3nenIOH2Ll0oUc3J9Vol7ReHrxXrm3MDIcs3DeTDp378l1A28GYOvGdXij66D4/T7aXt4Bm83OlVf3os/Am9m8bg0tLriIc+s3rKRWls2qpYuY9J/3aXp+K+649wGcLpfVIVUp//rXvxg0aBA+n49evXrFygsKCjh06BANGjSwMDpVXZ0y0YuIB5gHuKP1PzfGPC0iLYDxQG1gOXCXMcYvIm7gfeAK4CBwqzFmWwXFXym+/PjfvPvqc3gLC5k3Ywo2m42A3x/rhdvtdkKhMCKRWRjhUAiDYLfbCAVDrFjyA36fj4WzZuD2JNHjugFc0+8GJn4yDgRCwRAb1qwi4Pezce0qatVJp+9/3WJ1s4+z8+ctPD70t/i8hbg9Hg7t38cjOq++BBGhS5cuAKSmppKbmwtEOgt16uh1FGWNsvTofUBPY0yeiDiB+SIyDXgUeN0YM15E3gbuBUZHvx82xrQUkduAF4FbKyj+cvf9d9PYsGY1V/bozeGD2az/cRVrViyN9b5DwSBFE+mCwQA33HoXF7S5hEuu6ITf72Pd6hW0u6IjoWCQn1YtZ/nC75k1bRIAXm8hMyd/yeZ1a2nf8UoeemoEc76ZQjgcYun8OZE6hYWsXLqwSib6zevXxmYV+bxeMhfM472RL9GgSTOuG3gzNlti33935MgR3n77bUSEP/zhD6SlpZ2wbigUiiV5iPwBOHDgAA0bVq1PaKp6OK0bpkQkGZgP3A9MAeobY4IiciXwjDHmOhH5Jnq8UEQcQBaQbk7yg6rKDVNTPv+YUc8/hbewELvDEb1I6osc2+wEAn7sdkekRx/w43K7ee29T2lz2YnvV1gyfw5P/ele/D4vDqcTDAQC/sit8dE1apwuF5jIHw6X28PTr43mqmv6VGLLy2b/3t38/lfX4C0swOlyEQqGCIWCuNwefnPnPdw37EmrQ6wwxhguueQSNm7cCMBFF13EypUrTzpfvkaNGrFk73Q6ycvLw6VDXaoclesNUyJiB5YBLYE3gS3AEWNM0WD0LqBR9LgRsBMg+kcgB6gDHDitFlSwH5cvYeuGdXTo2oPDB7LZvH4ts6d/Heu5h0MhQsXG2nte/2saNGpCh649yMs9yrpVy8m4qvtJkzxAx649+Mfb77M6czFbNv7EvBlTgcgdlUUCfj9XXdOHC9tcQpvLMsi4qnsFtPjsndugEW9/NpU50ydz9Mghpnz+CcGCAD5vIXO/mUKDxk2p36gJHbv2sDrUcnf06FHWr18fu/6ydu1axowZg9vt5pZbbjnuIuuxPXqbzaY9emWZ0+3R1wS+BP4GjDXGtIyWNwGmGmPaicha4DpjzK7oc1uAjsaYg8e81lBgKEDTpk2v2L59e3m0p0xmTvmKl/46DDCEw+For0wIh0PYbXZ8Pm+sF+/3+3B7knjxnQ9p3/HKs/q582ZM5fn/+RNebyEOpwsh0rv3JCUx/O+v0uv6X5dH8yrFnp3buWdgz8isIbcHYww2myBi4677H+a3Qx60OsRyFQ6HSU5OxhedTmuz2fB4PIgIl19+OXPnzj2ud689elXRKmQJBGPMERGZA3QGaoqII9qrbwzsiVbbBTQBdkWHbtKAQ6W81rvAuxAZujmdOE7H9q2b2LVtK+0u70iNmrUAmPL5J/i8kZ67iMTubHQ4nXTu0ZukpGSuvm4A2fuyWLNiKT2uG3DWSR6ge5/+BIJ+Fs+bTceuPbDZbCyaN4uMq7rHVZIHaNikGa++9ylff/ohoVCQ+TOnU5ifD8C0CZ8kXKLPzc0tMZsqHA5TUFAAwMKFC5k8eTLJyclcc801sSURju3RHzx4UGfdKEuUZdZNOhCIJvkkoDeRC6yzgZuIzLz5HTAxesqk6OOF0ednnWx8viItmP0tzz56H3a7A5fLzXuTZlG7bjoXtr2EtSuX4vN6sdntsVk0DoeDX91yZ4mhkxtuvbNcY+rV/9f06v9LUu/Zf2C5vn5latP+Ctq0v4LtWzfFhqScLhctW7c9xZnxJzU19bhNRYp68Ha7ndtvvx0RoXfv3nz55ZfY7XZSUlLIy8sDIn8YatWqZUnsSpWlR98AGBcdp7cBnxpjJovIT8B4EXkOWAEULWo+BvhARDYT6cnfVgFxl8n4997C543ccRoKhVgwewYDbv4t9/xpOMYY1q1aTp+BN7N7x8+sXZFJn1/fXGXHx6uyZue14pnX3+E/Y9+mSfPzuP/xv1kd0mkzxvDTTz/hcrlo1apVifL169dTUFBQokcvIvTv3x9jDNOmTYsN6UyZMoWlS5fi8XjIj37CKXLkyBHq169fOQ1SqpiEXqb4paeG8e2kCQT8fjxJyTw3agwZXXSRKXW8e++9l/Hjx2OMYdiwYfz9738H4IEHHmDs2LEYE7me4/f7AUhLS+PQoUOICHXq1OHw4cNAZO0bm80WGxL0Rjsa55xzDtnZ2brblCpXZR2jj+tE7/N5ycvJoXb6uaVOc8vPy+XFJx9ly/qfuP7mO7hj8APlEa46Q8YY9u/fT2pqKslVZKeqQ4cOceTIEVq3bh1L4na7nb0ZsfmbAAAO5ElEQVR791JYWEjLli1jM6RsNhvt2rXDZrPxzjvv0KFDZO/cZcuWcd999xEIBFizZk1sm0aHw0Hbtm0JhUK8+eabdOvWzZpGqoSV8Il+/ZpVPDroFgJ+H63bXcar7/1Hb8evwsLhMDfddBNTp07FbrczceJEevfubWlMTz75JK+88krsYnxRQne73bFkXbw8OTmZo0ePnnAp6nA4TFpaWmxc3uPxcOjQIe3FqwpT1kQft7cyvvnC0xTk5RLw+9m0bg0L5nxb4vni/1GV9RYsWMC3336Lz+ejoKCAP/7xj4TD4ViSLX5c0cLhMNnZ2bzyyiv4/X4CgQDhcJh69erRsGFDgsEggUDguPKJEyeedL8Bm83GxIkTadSoEfXq1WP8+PGa5FWVELeJ3ulylRiuKdqEefeObdzaqyO92jbhbw8NKTFLQlnH5XKVSOSHDh3C5XKRnp7OnXfeicvlom7duixatKhC4xg5ciRJSUk0a9asRGfA7Xazd+9eNmzYUCKZO51Odu/eze7du8v0CaRnz57s2rWLrKwsBg6M3xlVKrHEbaJ/+G/PUye9HiJC5+696Hx15D/hqOf/xv6sPRhjWDp/DvNnTrc4UgXQoUMHBg0ahN1ujw1vhEIhDh48yMcff0woFOLQoUPcfffdFRbDwYMHefzxx/H7/RQWFmKz2XC5XCQlJfHBBx8gIqSkpDBq1ChcLhcej4exY8ee0a5hSlUlcZvoa9aowbUd23LtZa24sv3Fsd693+/FRHtqxpjIBtrKciLCqFGj8Pl8TJs27YRbHh44cIC6dety0UUX8dRTT1G3bl0uvPBCfvzxx7OOIRgMlvhUYYwhLS2NunXrkp6eHisfOnQohYWFFBQUcNttls0OVqrcxG2i//qdF9i7dT2+gjyWTP+MrauXAHD/Y0+RkloDl9tNi1YX0e3a/hZHqoqz2+106tSJ6667DrfbTVJSEl27dsXtdsfmnh88eJANGzYwYsQIDh48yMaNG7n55pvP+mfXq1ePhx56CJfLFRtKys7OZufOnbGdoIoUTZFUKhHEbaLPP3Io1nMXoDA3B4BWF7djwvcr+eibBbw1/mvc0f1WVdVhs9mYMGEC27ZtIzs7m3nz5rFt2zamT59+wrVg9u7dS5s2bejRowc///xzqXVmzJhBu3bt6Nq1a2yVyWP179+fli1b0rp16xKJPC8vr9IuBitV2eJ2euWWlYv45B+PITYbNWqnM/TlD3AnVY252erMhEIhevfuzdKlSwkGg5x33nns2LGDYDBIKBQiGAxis9lo27Ytq1atKnHuwYMHadq0KQUFBYgI559/Pps2bSpRJy8vj/r168fuWE1KSord2DR8+HCeffbZSmurUuWhQhY1q0rOb9+Zh0Z/xdED+6jXvBUOp86hj3d2u53vvvuOVatWUadOHZo0acKqVavYvXs3t99+O3l5eYTDYXbs2BE7Z93qFbw36mV8QSDaaTHGsHPnTvr370+NGjW45557eOONNxCREssY+Hw+VqxYgcfj4YILLqjQtu3YsYNhw4ZRWFjICy+8QLt27Sr05ylVXNz26FX1EQgE6NSpE5s2bSIcDvP444/z9NNPk3s0h1uuyaCwIB+bzcbarDwK/KHIRfhAAL/fH9v1KhwOx5YWttlsGGO46667GD16dKW0oWXLlvz8888YY6hZsyZ79uzB49FhRXV2Er5Hr6oPp9PJggULmDlzJrVr1+bKKyNLRmdn7YnVCYfDXNqkNo+88CaHDx9myJAh+P3+EnPlw+EwPp+PyZMnk5SURPfuFbuA3cGDB3nhhRfwer1s3bo1dg2gsLCQRx55BLvdzvDhw2nWrFmFxqGU9uhV3Ar4/fz+hms4sC8LEaFH31/xxPOvEwwGadeuXWyIx+l0EggEEBH69evHZ599VinxXXrppaxbt45wOIzD4cDhcGCMwRgTu+6Qnp7Ojh07dEMSdUa0R68SntPl4p3PpjFn+tekpNage5/rgchiYkuWLOHTTz8lJSWFvn37MmHCBJKTk8tlmubJ5OfnM2bMGPx+Pz/++OMvm9o4HLz44oskJSUxePDgWHlubi67d++mRYsWFRpXovv2229ZvHgxvXv3pnPnzlaHU+Voj16pctSlSxeWL18e3VrRRjgcJhwO07Rp09jyCj169GDx4sWEQiEaNGjA5s2bY0t4qNM3YcIE7r77brxeL263m5kzZ8aG9xKd9uiVqmSBQICFCxfGeusul4sbbrgBEeGVV15h4sSJ+Hw+JkyYwNixYyksLOS+++7TJH+WJkyYENvW0ev18s0331SbRF9WmuiVKidOp5NWrVqxdevW2CyfqVOnIiLMmjULr9eLiHDxxRezcOHC2IwgdXauvvpqJk6cSEFBAUlJSXTq1MnqkKocTfRKlaO5c+fyzDPP4Pf7GTt2bKy8+LaCq1evZs+ePTRu3NiKEBPO0KFDCYfDzJ49mxtvvJF+/frFntuzZw+rV6/msssuo169ehZGaS0do1eqgjRt2pTdu3cTDocREUQktjnJ3r17da36CrZy5Uq6deuG3W6PrGa7dGmF3xhX2RJ+4xGlqrpZs2bRt29fevXqxaRJk+jXrx89e/Zk1qxZmuQrwejRo8nLyyMnJ4fc3FzGjRtndUiW0aEbpSpIy5YtmTJlSuzxgAEDLIym+mnWrBkejwev14vH46Fp06ZWh2QZTfRKqYQ0bNgwNmzYwJw5c+jfvz+DBw+2OiTLaKJXSiUkt9tdrYdritMxeqWUSnCa6JVSKsFpoldKqQSniV4ppRKcJnqllEpwmuiVUirBaaJXSqkEp4leKaUSnCZ6pZRKcJrolVIqwZ0y0YtIExGZLSLrRGStiDwULW8vIotEZKWIZIpIx2i5iMhIEdksIqtF5PKKboRSSqkTK8taN0FgmDFmuYikAstE5FvgJeBZY8w0EekffdwD6Ae0in51AkZHvyullLLAKXv0xpi9xpjl0eNcYB3QCDBAjWi1NGBP9Hgg8L6JWATUFJEG5R65UkqpMjmt1StFpDlwGbAYeBj4RkReIfIH46potUbAzmKn7YqW7T3LWJVSSp2BMl+MFZEUYALwsDHmKHA/8IgxpgnwCDCmqGoppx+3X6GIDI2O7WdmZ2effuRKKaXKpEyJXkScRJL8R8aYL6LFvwOKjj8DOkaPdwFNip3emF+GdWKMMe8aYzKMMRnp6elnErtSSqkyKMusGyHSW19njHmt2FN7gKujxz2BTdHjScDd0dk3nYEcY4wO2yillEXKMkbfBbgL+FFEVkbL/gIMAd4QEQfgBYZGn5sK9Ac2AwXAoHKNWCmlEkAwGGTLli00bNiQ1NTUCv1Zp0z0xpj5lD7uDnBFKfUN8MBZxqWUUgkrNzeXTp06sWPHDux2O7NmzeKKK45Lp+VG74xVSqlK9tlnn7F9+3by8/M5evQof/3rXyv052miV0qpSpaSkkLk8ifY7XZq1KhxijPOjiZ6pZSqZDfeeCMDBw7E7XbTunVrXnvttVOfdBYkMqRurYyMDJOZmWl1GEopFVdEZJkxJuNU9bRHr5RSCU4TvVJKJThN9EopleA00SulVILTRK+UUglOE71SSiU4TfRKKZXgqsQ8ehHJBraf4el1gQPlGE48qG5t1vYmturWXii/NjczxpxynfcqkejPhohkluWGgURS3dqs7U1s1a29UPlt1qEbpZRKcJrolVIqwSVCon/X6gAsUN3arO1NbNWtvVDJbY77MXqllFInlwg9eqWUUidR5RO9iLwnIvtFZM0x5X8SkQ0islZEXipW/mcR2Rx97rrKj/jslNZeEWkvIotEZKWIZIpIx2i5iMjIaHtXi8jl1kV+ZkSkiYjMFpF10d/lQ9Hy2iLyrYhsin6vFS2P6zafpL0vi8j6aJu+FJGaxc6J9/d0qW0u9vxjImJEpG70cUL+jqPPWZO3jDFV+gvoDlwOrClWdg3wHeCOPj43+v1iYBXgBloAWwC71W0oh/bOAPpFj/sDc4odTyOyp29nYLHV8Z9BexsAl0ePU4GN0d/jS8AT0fIngBcToc0naW8fwBEtf7FYexPhPV1qm6OPmwDfELmPpm6C/44ty1tVvkdvjJkHHDqm+H7gH8YYX7TO/mj5QGC8McZnjPkZ2Ax0rLRgy8EJ2muAor3G0oA90eOBwPsmYhFQU0QaVE6k5cMYs9cYszx6nAusAxoRadu4aLVxwK+jx3Hd5hO11xgzwxgTjFZbBDSOHifCe/pEv2OA14HHibzHiyTk7xgL81aVT/QncAHQTUQWi8hcEekQLW8E7CxWbxe/vKHi2cPAyyKyE3gF+HO0PKHaKyLNgcuAxUA9Y8xeiPzHAc6NVkuYNh/T3uLuIdKjhQRqL5Rss4jcAOw2xqw6plrCtPmY37FlectRni9WiRxALSIf6zoAn4rIeUQ+6h0rEaYV3Q88YoyZICK3AGOA3iRQe0UkBZgAPGyMOVq0cXJpVUspi7s2H9veYuVPAkHgo6KiUk6Pu/ZCyTYTaeOTRIasjqtaSlnctbmU97RleStee/S7gC+iH+2WAGEia0fsIjLmV6QxvwxzxLPfAV9Ejz/jl491CdFeEXES+Q/xkTGmqJ37ij6uR78XfcyN+zafoL2IyO+AAcBvTXTwlgRoL5Ta5vOJjEevEpFtRNq1XETqkwBtPsHv2LK8Fa+J/iugJ4CIXAC4iCwQNAm4TUTcItICaAUssSzK8rMHuDp63BPYFD2eBNwdnaXQGcgpGu6IFxLpuo8B1hljXiv21CQif+CIfp9YrDxu23yi9opIX+B/gBuMMQXFTon793RpbTbG/GiMOdcY09wY05xIsrvcGJNFgv6OsTJvWXFV+nS+gE+AvUCAyJvh3ug/0IfAGmA50LNY/SeJXLXeQHSmSjx9naC9XYFlRK7MLwauiNYV4M1oe38EMqyO/wza25XIx9TVwMroV3+gDjCTyB+1mUDtRGjzSdq7mcg4bVHZ28XOiff3dKltPqbONn6ZdZOov2PL8pbeGauUUgkuXodulFJKlZEmeqWUSnCa6JVSKsFpoldKqQSniV4ppRKcJnqllEpwmuiVUirBaaJXSqkE9/8BGyLNoAmg3CcAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "xs = range(jmin, jmax)\n", "ys = range(imin, imax)\n", "xx, yy = np.meshgrid(xs, ys)\n", "plt.scatter(xx, yy, s=10, c=days, cmap='copper');" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAADxCAYAAADr5V2tAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAACfdJREFUeJzt3V+o7llZB/DvM8eZJstB6NgfVNBQvBkw46DFQFhmTDXYTRcKCdVF3SRKRjRCSEFQN9FAXrTxT0lmhDUgkZY3wyDU4D4yVM4Y1FB4nKlxX6Rjg9mc83RxdrBx3O9+19v7smYxnw8czv7tvdaP39XD4lnPs1Z1dwBYxy2zPwCAMQI3wGIEboDFCNwAixG4ARYjcAMsRuAGmKiqXlNVD5/595WqetfGOeq4AZ4bqupSki8meUN3/9t546y4AZ473pTkXzYF7SR5wcgb7/6BV/XJl5/+f30VwPPN1c8/8bkkXzvzq6PuPvomQ9+a5KMXvW8ocJ98+ekcf+gXRqYAPO/VD/7G17r7ysYxVbcleUuSey9631DgTpLIiQMcwo8n+Wx3/8dFA8cCd0fgBjiMt2WLNEmy04p7eAYAG1TVC5O8OckvbjNeqgRgsu5+Osl3bDteOSDAYgZX3G3FDTCZHDfAYsarSm6I3AAzja+4LbkBppIqAViMckCAxVhxAyxGyzvAYtRxAyxGqgRgMTYnARbjrBKAxVhxAyxGyzvAYmxOAizGWSUAi7HiBliMBhyAxWh5B1iMVAnAYtRxAyzGihtgsqp6cZL3J7kzN6Psz3f335433oobYL77knyyu3+6qm5L8sJNg3eo4wZgX6rqjiQ/lORnk6S7v57k65vmaHkHOLzLVXV85vmou49Of/7eJF9K8qGqem2Sq0ne2d3/dd7L1HEDHN5Jd185528vSPL9Sd7R3Q9V1X1Jfi3Jr5/3Mse6Asx1Lcm17n7o9PljuRnIz2VzEmCi7v73qvpCVb2mu/8pyZuSPLJpjnJAgPnekeQjpxUljyX5uU2DtbwDTNbdDyc5Lwf+LFIlAIuRKgFYjBU3wGKsuAEWM7g5qQEHYLbhFXcL3ABT7RC4D/EZAGzLihtgMUOBu2NvEmC24c5JK26AuYZTJY7jBphrOFVyQ+QGmGp8c/IQXwHA1oZvwLkhxw0w1Q6bkwf6EgC2Ml4OKHIDTCXHDbCY4RX3dVUlAFNpeQdYjM1JgMWMd07KcgNMtUNVyYG+BICtjDfg2JwE2Kuq+tckTyW5nuSZ7r6yafwO5YACN8AB/HB3n2wzcPyQKXEbYCrncQMc3uWqOj7zfNTdR2eeO8nfVFUn+YNv+NuzuHMS4PBOLshb39Xdj1fVdyb5VFV9vrsfPG+wBhyAybr78dP/n6yq+5O8Psl+AncnuS5wA+xNVX1bklu6+6nTn38syW9umiNVAjDXdyW5v6qSmzH5T7r7k5smDG5OtlQJwB5192NJXjsyRzkgwGI04AAsRo4bYDEuUgBYjFQJwGJcpACwGJ2TAIuxOQmwmMHNydbyDjCZVAnAYmxOAixmh5Z3kRtgph3quAGYSY4bYDHjVSVa3gGmsjkJsBhnlQAsRuckwGKUAwIsZnjFLXADzCVVArCYsVRJq+MGmG2HVMkhPgPg+a2qLiU5TvLF7r5n09jBwN1W3ACH8c4kjya546KBzioBmKyqXpbkJ5P8VpJfvmi8W94B5vu9JL+a5EXbDN6h5V3gBhh0uaqOzzwfdfdRklTVPUme7O6rVfXGbV6mHBDg8E66+8o5f7sryVuq6ieS3J7kjqr64+7+mfNeNpwqcVYJwP50971J7k2S0xX3r2wK2olyQIDluEgB4Dmiux9I8sBF44YvUnBWCcBcw1UlUiUAc0mVACxGOSDAYsbLAUVugKmGV9zXBW6AqdzyDrAYqRKAxTjWFWAxwxcpaMABmGuHVMmBvgSArYx3TmqdBJhq/HRAWW6AqaRKABbjrBKAxTirBGAxctwAixnLcXerKgGYTKoEYDE2JwEWM1wOeONAHwLAdqy4ARaj5R1goqq6PcmDSb4lN2Pyx7r7vZvmjKdKxG2AffrvJD/S3V+tqluTfLqqPtHdf3fehB3O4xa5Afalb+afv3r6eOvpv42BVjkgwGRVdSnJ1SSvSvK+7n5o0/jhixRsTgIMu1xVx2eej7r76P8euvt6ku+rqhcnub+q7uzufzzvZXLcAId30t1XLhrU3f9ZVQ8kuTvJfgJ3Oq4uA9ijqnpJkv85DdrfmuRHk/zOpjnquAHm+p4kf3Sa574lyZ91919umiBVAjBRd/99kteNzLHiBljMDnXcAMw0vDlpxQ0w12COu3NdkhtgKqkSgMXYnARYjLNKABYzeFmwzkmA2YZX3AI3wFxSJQCL2aHlXeQGmMmKG2Ax4xcpqOQGmGqHlvcDfQkAWxnOcWt5B5hL5yTAYpxVArCY4VSJFTfAXDtcFnygLwFgK+Mt7yI3wFTDFymo4waYa4dDpg7xGQBsSzkgwGJ0TgJMVFUvT/LhJN+d5EaSo+6+b9Oc8XJAOW6AfXomybu7+7NV9aIkV6vqU939yHkThlMlWt4B9qe7n0jyxOnPT1XVo0lemmR/gVuqBGDY5ao6PvN81N1H3zioql6R5HVJHtr0Mp2TAId30t1XNg2oqm9P8udJ3tXdX9k0drwccHQCABtV1a25GbQ/0t1/cdH4waqStuIG2KOqqiQfSPJod//uNnPkuAHmuivJ25P8Q1U9fPq793T3X503wWXBABN196eT1MicHVreBW6AmXaoKjnQlwCwlR1a3kVugJlsTgIsZvyWd5EbYKrBFbc6boDZpEoAFrNDy7vIDTDTWI7bRQoA07m6DGAxLlIAWMwOV5cBMJNUCcBi3PIOsBgNOACL2eE87gN9CQBbcR43wGJsTgIsZrhzUqoEYK7xFbdKboCpnA4IsJjxixTkSgCmUscNMFlVfTDJPUme7O47Lxq/Q44bgD37wyS/n+TD2wx2yzvAZN39YFW9YtvxOicBDu9yVR2feT7q7qNdX6YBB+DwTrr7yr5epuUdYDHquAEWs0PLu8gNsE9V9dEkb8zNXPi1JO/t7g+cN96KG2Cy7n7byPjxBhyV3ABTWXEDLMZZJQCLcawrwGLc8g6wmB1a3kVugJlsTgIsRss7wGIGUyUtcANMZnMSYDGOdQVYjM1JgMUMlwNqwAGYa3jFreUdYC6bkwCLkSoBWMwODTiH+AwAtjV+kYJcCcBUygEBFjN+kYLIDTDVDlUlAjfATFIlAIsZLwcUuQH2qqruTnJfkktJ3t/dv71p/Hg54I4fBsCzVdWlJO9L8uYk15J8pqo+3t2PnDdnPHAr5AbYp9cn+efufixJqupPk/xUkr0F7r9OcnnnzwN4frq9qo7PPB9199Hpzy9N8oUzf7uW5A2bXjaW4+6+e2Q8ABeqb/K7jamNWw70IQBs51qSl595flmSxzdNELgB5vpMkldX1Sur6rYkb03y8U0ThjcnAdif7n6mqn4pN/cQLyX5YHd/btOcUpcNsBapEoDFCNwAixG4ARYjcAMsRuAGWIzADbAYgRtgMf8Lkyb8poksUwkAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(1, 1)\n", "cnorm = mplcolours.Normalize(vmin=0, vmax=7)\n", "colorbar.ColorbarBase(ax, cmap='copper', norm=cnorm);" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "[2214051840 values with dtype=float64]\n", "Coordinates:\n", " * time (time) datetime64[ns] 2018-01-15T00:30:00 ... 2018-01-21T11:30:00\n", " * grid_y (grid_y) int16 0 1 2 3 4 5 6 7 ... 888 889 890 891 892 893 894 895\n", " * grid_x (grid_x) int16 0 1 2 3 4 5 6 7 ... 388 389 390 391 392 393 394 395\n", " * grid_z (grid_z) int16 0 1 2 3 4 5 6 7 8 9 ... 31 32 33 34 35 36 37 38 39\n", "Attributes:\n", " standard_name: OilWaterColumnOilVol_3D\n", " long_name: OilWaterColumnOilVol_3D\n", " units: m3/gridcell" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "results.OilWaterColumnOilVol_3D" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "oilwaterhours = np.ma.masked_array(results.OilWaterColumnOilVol_3D[:, 39, imin:imax, jmin:jmax].sum(axis=0)/(440*500)*1000, \n", " mask = results.OilWaterColumnOilVol_3D[:, 39, imin:imax, jmin:jmax].sum(axis=0)/(440*500)*1000\n", " < 0.004)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "beachvol = np.ma.masked_array(results.Beaching_Volume[imin:imax, jmin:jmax],\n", " mask = results.Beaching_Volume[imin:imax, jmin:jmax] == 0)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0.18833146875494242\n" ] } ], "source": [ "print (beachvol.max()/440/500*1000)" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA30AAAFJCAYAAAAv/kQrAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzsnXecXGX1/99nZlvapm0CqSQkpEDoLQSkVxWkiWDnJ+JXBRRrBEVA0GAXUBARAUVQkWroCERaQigBQhIISUgjyWbTt8/M+f1xZzeTybQ7c+/Mndnzzuu+snPL85x757mfec5TziOqimEYhmEYhmEYhlGZhEptgGEYhmEYhmEYhuEf5vQZhmEYhmEYhmFUMOb0GYZhGIZhGIZhVDDm9BmGYRiGYRiGYVQw5vQZhmEYhmEYhmFUMOb0GYZhGIZhGIZhVDDm9BmGYRiGYRiGYVQw5vQZhmEYhmEYhmFUMOb0GYZhuERE5ovI0fG/l4nI8X7mkZxPIXmKyO0ico1HZhoVhl/lOUj0hHs0DMNIxpw+w+ghiMhFIjJXRNpF5PYs5w4SkftFpFlEPhCRT+d6PNu1bvIqxI40ef1CRN4RkaEZzvmiiLwlIi0iskZEbhKRAYnnqOpeqvpsprziaR0hIi+KyGYR2SAiL4jIwdmuc5NHUBGRZ0XkAh/SzViO4/m2ici2+LYox3RrReTP8XK0VUReF5FTks7JVFaz2TVGRB4RkY3xcnWjiFS5fwLZnZZCyl1QiN9ja/w73CgiM0VkVKntMgzDKFfM6TOMnsNq4BrgthzO/T3QAewCfAa4SUT2yvF4tmvd5FWIHTsgIhOBo4B3gIY053wbuA74LtAfmArsBjwpIjUZ7iFVWvXAf4AbgEHACOAqoN1NOsZO5FKOL1LVvvFtYo7pVgErcMpIf+BHwD9FZEzCOZnKWza7/gCsA4YB+8Xz+VqOtuVMhZW7U1W1L84zW4tzT4ZhGEYemNNnGD0EVb1PVR8AmjKdJyJ9gLOAH6nqNlV9HngI+Fy249mudZNXIXakubXTgbuBicB7KWypx6kcX6yqj6lqp6ouA87Bcfw+m3BuLsPDJgCo6t2qGlXVVlV9QlXfTErnB/Hex40i8hcRqXORByIyXET+LSKNIrJURC5JOr6/iLwW7736B1CXIa3J8Z6yTeIMLz0tydbviMib8R6kf3TZmoON00Xk/bgN74jIGUnHc04713LsFlVtVtUrVXWZqsZU9T/AUuDAuI0Zy1sOdo0F/qmqbaq6BngMSNdAkfZ5ichfgdHAw/FesO8lXZ5Lucv4fSTZskxEvhv/bprjvaG7iMij8eufEpGB8XPTlp+EtFyXIVVtA+4F9sz1HnLNS0Qmxd+bc7PZYRiGUc6Y02cYRjITgKiqvpuwbx7bK6iZjme71k1ehdiRimNxKuTzVbUzxfFpOA7RfYk7VXUb8ChwQpp00/EuEBWRO0TklK6KcQo+A5wEjMO5px/mmoGIhICHce57BHAc8E0ROSl+vAZ4APgrTq/Pv3Acl1RpVcfTegIYClwM3BXvIe3iHOBkHAdmH+CLOZr6PvARnB60q4C/iciwpHPyTTsVPxOR9eIMazw6nwREZBec72N+fJfb8pbM74BzRaS3iIwATsFx/FKR9nmp6ueA5cR7wVT150nX5lLucvk+EjkLp/xPAE7FeR8uw+kxDwGX5Fh+II/vWUR6A58CXnZ5DxnzEpED4vZerKr3ZLPDMAyjnDGnzzCMZPoCm5P2bQb65XA827Vu8irEjlTsAxyCM8wuFQ3AelWNpDj2IWmGhKZDVbcARwAK/AloFJGH4s5EIjeq6gpV3QBcC5znIpuDgSGqerWqdqjqknheXb0WU4Fq4Lfxnst7gVfSpDUV55nOiKf1X5xhgon2XK+qq+O2PowzTDErqvqv+HUxVf0HTk/rIUmn5ZV2Cr4P7I7jBN+C0yM2zk0CcQfmLuAOVV0Y3+22vCXzHI6DuAVYCczFcch3IsfnlZJcyl0e6d+gqmtVdRXwP2C2qr6uqu3A/cD+5FZ+wN33/ICIbMJ5ZicAv3B5D5ny+ghOT+0X4r26hmEYFY05fYZhJLMNqE/aVw9szeF4tmvd5FWIHTsgIn0BAXZT1VlpbFkPNEjq4BrD4sddoaoLVPWLqjoSmAIMB36bdNqKhL8/iJ+TK7sBw+PD6TbFK8iX4cw5I57WKlXVpDxSMRxYoaqxpHNHJHxek/B3C04lPysi8nkReSPBxins7ETnlXYyqjpbVbeqaruq3gG8AHw01+vjvad/xZm7d1HCIbdlOznNx3F6kfvg3PtAnPmjqc7P5XmlJVu5yyP9tQl/t6b43Jfcyg+4+55PV9UBQC3Od/GciOzq4h4y5fV/wIuq+kyG/A3DMCoGc/oMw0jmXaBKRPZI2Lcv24e5ZTqe7Vo3eRViRzL1OBXCTEMnX8IJdnFm4s74XK5TgKczXJuVeI/R7TiV00QSIxKOxgkIkisrgKWqOiBh66eqXU7Oh8AIEZGkPFKxGhgVd1ASz13lwp6dEJHdcHqcLgIGxyvxb+M44cVAc80r/pz+jOM0n5U0DNht2U5kEM73fGPcGW0C/kIKZzTH56XJ16Ujudz5+H34Un4A4nMT7wOiwBEe3cP/AaNF5DeF2mcYhlEOmNNnGD0EEamKBzIIA2ERqevq1RJn7bbbwQlogdMjcbWI9BGRw4FP4PR+ZDye7Vo3eRViR4rbPwjHAXpTRE5IcoKIp7cZZ27QDSJysohUixO58V84w/FSpZvpeU8SkW+LyMj451E4Q91eTjr16yIyUkQG4fTS/cNFNnOALSLyfRHpJSJhEZki28PzvwREcOZcVYnImaQfxjcbaAa+F7/3o3HmbxU616kPjpPSCCAi57Oz45szWcrxABE5qWufiHwGOBKnl22HspeGm4DJOPPlWhMPZCtvmexS1fU4QWG+Gj9vAPAFnDmByeTyvNbiDGFN9XyylTtPv48E/Co/iMMncHpHF+DNPWzFme93pIjMKNRGwzCMoGNOn2H0HH6IMxRrOk4kyla293yNwhkG18XXgF44IebvBr6qqvNzPJ7tWjd5FWIH0D0/61M4PTjPAROThjt2Ew+KcRnwS5x5RLNxetOOi89fcsNW4FBgtog041S63wa+nXTe33GCSSyJbzkvnK6qUZyK9X44TsV64Fac4BaoagdOz+UXgY04z+G+NGl1AKfh9Gqux5n7+PmEOW35oKr6DvArHAd0LbA3O37/bslUjqtxnl8jzj1cjDNEsGutvuSy10289+grOM9yjWxf5+8zCadlKm+Z7ALnezg5bttiHGf80mQ7cnxePwN+GB/a+J2kYxnLnQ/fR5fdfpSfh0VkG867eC3O/Lv5Xt2Dqm7CmSt4ioj8pAA7DcMwAo+kqfsYhtFDECfC4zxgnzRRLcsyr3JARJYBF6jqU6W2xWtE5DXganWWMSg5VvYMwzCMnkyqgAWGYfQg4i30kystL6N0iLNg+WTg9VLb0oWVPcMwDKMnY8M7DcMwDM8Qketwhqt+X1XTRQo1DMMwDKOI2PBOwzAMwzAMwzCMCsZ6+gzDMAzDMAzDMCoYc/oMwzAMwzAMwzAqGHP6DMMwDMMwDMMwKhhz+gzDMAzDMAzDMCoYc/oMwzAMwzAMwzAqGHP6DMMwDMMwDMMwKhhz+gzDMAzDMAzDMCoYc/oMwzAMwzAMwzAqGHP6DMMwDMMwDMMwKhhz+gzDMAzDMAzDMCoYc/oMwzAMwzAMwzAqGHP6DMMwDMMwDMMwKhhz+gzDMAzDMAzDMCoYc/oMwzAMwzAMwzAqGHP6DMMwDMMwDMMwKhhz+gzDMAzDMAzDMCoYc/oMwzAMw/ANETlZRBaJyGIRmZ7i+LdE5B0ReVNEnhaR3RKOfUFE3otvX0jYf6CIvBVP83oRkWLdj2EYhleIyEQReSNh2yIi3/QlL1X1I13DMAzDMHo4IhIG3gVOAFYCrwDnqeo7CeccA8xW1RYR+SpwtKp+SkQGAXOBgwAFXgUOVNWNIjIH+AbwMvAIcL2qPlrMezMMw/CSuF6uAg5V1Q+8Tt96+gzDMAzD8ItDgMWqukRVO4B7gE8knqCqz6hqS/zjy8DI+N8nAU+q6gZV3Qg8CZwsIsOAelV9SZ2W6zuB04txM4ZhGD5yHPC+Hw4fQJUfiRqGkZqTTt5Hm9Zvc33dq68ufQ2nlfthVX3Yc8MMw+jRHH9ib21qirq+7o3XOuYDbQm7blHVWxI+jwBWJHxeCRyaIckvAV09dqmuHRHfVqbY3+NoGDxQx4zukbdeUta016Y9tmtte97Xpksn12uCSLbn4RevvjF/vaoO8Sv9k6eO1/WbW7KfCLy68MNsOtnFucDdHpiXEnP6DKOINK3fxuy5P3F9XZV8VlX1Qh9MMgzDoKkpyrMvuXceBtQubVPVgzKckmquXcp5JSLyWZyhnEdluTbnNCudMaNHMPeZf5bajB7HjKV7pD02fex7eV+bLp1crwki2Z6HX8jAvXzpLeti/eYW5v4lt2qZHHZVNp1ERGqA04AfeGBeSlw5fdaiVFwKaUnKdn26dMq1NalcWpIUJRZz35puZMa0qfhYS3dulIs2gUAs7IcpK4FRCZ9HAqt3yl3keOBy4ChVbU+49uika5+N7x+ZtH+nNA2jmHjt3MxYukd3mtPHvleWjl++z6TrXkvlMOaMt3FRTgFeU9W1XiaaiCunz1qUiku6FzzXl6CntCaVUlTctyQpqpG88zNSY9pUfFJphdfalJxmT9GnkmiTgsR8CYD5CrCHiIzFCVBwLvDpxBNEZH/gj8DJqrou4dDjwE9FZGD884nAD1R1g4hsFZGpwGzg88ANfhhvGNlw+566ceAqwfHLRLb7yXY88E6hO87Dx6GdYMM7yw43BTxXgTBR2X7cdwFRULWePqPy8EObvLouqLit8PivT947faoaEZGLcBy4MHCbqs4XkauBuar6EPALoC/wr/jKC8tV9bS4c/cTHMcR4GpV3RD/+6vA7UAvnDmAFrnTKApe1BlMA72jqHW4RBTPevpEpDdOhOOveJJgGszpKxMqrDWjaLgVR78rWYoSs54+o8IwfSoOflZuBN96+lDVR3CWVUjcd0XC38dnuPY24LYU++cCUzw00zBc4afumabmRnKdregOsUejO+PRiwd7k1p6bMmGgFLsgluylpKA0/VcZizdo3srDGd4p9vNMIKEVy3dbvMytuO5NilIzP1mGJVO4jtW7Llmqd5tq6MFCNXctoBgPX1lQDFak1KJilW2HLxtSbI5fUblYNpUWjxv5TYnzjA8J52jWIiWpdNHw8iEOX0BodCgLV7knTi3z/AJVTRmTp9RXhQStMWr/BPzM8fPBxQkOA3ShhEIvNQZN3Ws5FgLRhAJVi9eLtjwzgqikGEA08e+l7IVKuhiE3T7UqIR95thlAgvKj3phiBmen8TW8dTnVcO73452JiIDe80DO9JpwPZnMByb9jyar3CQKM5bgHBnL4KItULVumiUn7YnD6j55Gp0pOJctansqzwxNT9ZhgVTFejU/L77Pb9TTccM5A6YORGV/TOMprTZ05fQPBSWJJb1E1UgoRCrNP9Zhglwu9KjxEQ4sM73W6G0VPwqufeTTCYchstUGrS/V75hvX0GYXipbCYqHhPIaKimndPX38RuUVETvX8hgwjR0pR6TFyp+AKTyyPzTAM1+QacbdcG8nKcqRDXpSX12dOXxlQyMtRyaJSnijEIu432KyqF6rqw6W+A8PoolBtMu0JDs46fep6M4yeQqFDM/NtjMmUT1Abzrxb5srwEoveWQYUK4JnUMUjHT2nJckwgolpU2pMmwyjPCiXUQflvFZfuTzjvCizdi/r6QsgpZrom2nZiCC+rGXZkqR59/QZRskJWku3aZOHKDa80+hRFLJGnl9zlIOoZ16RuDxYxVBmgVysp69IlEtLR+LLmGhr0F/Scnm+gC3BYASKdO980EjX0m3a5B0WmMXoqSSOKOj6O5u25LNeqFdrjAZd99JRrnanpcw003r6ikA+hdzPCERBbR33gqC3JAmKxCKuN8PwA2vpLh5l0ftnPX1GDyaX9zPVesZBICh2lJqi6msZLtlgPX1FJlWruluRyYZfLUlepes3fttY0ByjruGdhhFAkvUp1buUWPbzeQ96uj75SWHaZIutG0YXpZxDl6+WlcvoB7/p6fefCevpKzHJlSyveuH8WqqhknsJ3ZB/i73N6TPKg1SVnmK9+/nmY9pUYFTUXFutA9iCbRhuyaRphTZuFUqmPIPa29jzKD+dNKevCOQiLPbS9hQU0YjrzTD8wCo9RiISc78ZRrkTJP3IdY51umPWy1VkymuZPnP6eiImKiVEgVjU/WYYPhKkSk8ibvXJtKkALHqnYXSTqCVB1RUbeRUArKfPSEfyy5k4DKfYopLrPEI/A8r0TCyQixFMUulTqr/9xs0cZ9MmbxF1vxlGJVDsOXuZ8isk6JPpoZEJc/qKTLGFJRsmKsVGrafP6PHkqk3FWgfQiGM9fUYPptj6YXpVAVhPn1EuZKpQeR1R1Ihji7MbPYBoVHn33a1s3NiR9pxsLd35HMuWrpEBG95pGAUTbW1m45yn2DjnKaKt2zxLN3HUVeCXfukp5DqfLzg+nzl9pSBovX1GcZFY1PVmGMXAC73o6Ijx8VNf5KSTX2D/A57m2Wcbfc3P8AYBRMX1ZhiVRCGapKqsnXknW9+Zy6b5c1k786+oj7085vgFAOvpM/LF7QvsR2tP8qLL6dI3sckTteGdRrAp1BF78qm1zHt7Ky0tUVpbY1x11QLP7DJt8hnr6TOMvIls2Uhk60bQGGFidGzbgna2Z7wm36WzujC9KzHW02cEDS9ExYTFO/Ls6esvIreIyKmltt8w0qGqzJjxLtGO7Q0V9f2rPc0jV+fPcIkN7zSMtOSiM03Pz+z+W4FtdYOR6tqC805Vh7Po6kGhvLw+c/oCRq4vbK+hJ3PVoeM8zTvXta9MVAoh756+zap6oao+XOo7MHomubz3a9a0s2xZ8w77fvXLvQvOO9dKj2lTgbiZoxKsuoxhlJzoto3b/5Ywr+1+HiKFD4HOpmupluGyofNFosx00pw+o5tcK0wmKvkjqjanz6hYBg6sprY2DICI0G+Xfowf37fgdHOt9JR6QflyR2LiejOMSqKQhqN+kw9GqqqRqmpaawexdtDkoo1GsFEPRi6Y01emtK57jNZ1jwGFrenihsQ8rEJVADanzwgwhehIXV2Y+/49lWnTBrHboaM5+4azilrpAdOmvMmnly9ALdiG4QX5RhUGqN/nMIaeeB6DP/JxHj3oKqLh7UM7k691U29zo2mJaz+bFvpNjkFcLJCLkQk3FaSrDh23wzDPdMKSDTeLHJuoFID19BkBp5BKD8A++/Tn/vsO44zfnMGAUQPSXuunNuVqq5FETNxvOSAiJ4vIIhFZLCLTUxw/UkReE5GIiJydsP8YEXkjYWsTkdPjx24XkaUJx/bz7DkYRp7U7jKS3mMmE6nq1b0vcQ5y4v9+Y3U0nynDxrGqUhtgeEOyqCT/7TXJjp9VsAzDSEfi2lJ+60Xy/BbTptIiImHg98AJwErgFRF5SFXfSThtOfBF4DuJ16rqM8B+8XQGAYuBJxJO+a6q3uuf9YaxI/noSb71smwOm+nbjiQ/rx8UJdcAeXQ5YE5fifDqZf3x7PcBmLF0+z4TFf/wRFSs584IOJne+3xaj/2q9Bjb8UabPDElmUOAxaq6BEBE7gE+AXQ7faq6LH4skwVnA4+qaosvVhpGBdJT6m8l+60oL5/PnL5S4SYakxEMvBEVRWIW59wINqZPPZD8hyE1iMjchM+3qOotCZ9HACsSPq8EDs0jn3OBXyftu1ZErgCeBqarauZF0QyjjPBieGZPaLgvaeNggObr5YLN6evhVLoYeIVnoqJYIBfDyBHTp+y4mfOYmTzm8zlz+tar6kEJ2y07J7wTrmpKIjIM2Bt4PGH3D4BJwMHAIOD7btI0jFSUk+bkYqvN6/OZMpvTZ06fkRETFa/Je50+wygK5VTpyUaXNpk+5YiK+y07K4FRCZ9HAqtdWnYOcL+qdnabqvqhOrQDf8EZRmoYeZNN+7q0pBA9KWRofK5ppbKzkjWwpL9ZHkXvFJEBInKviCwUkQUicpgf5prTVwJyFZZ8MVEJNqIx15thFIMgVnrySb8nOXqeVXgUJOZ+y4FXgD1EZKyI1OAM03zIpXXnAXcn7oj3/iHO6tenA2+7TNMw8mYHjclSqS/l2sY9pY5WAY2VvwMeU9VJwL7AAj8yMacvwLitvJTqhe4JouJdxcp6+ozKoEub1i9p4t6L7uVfF93L+sXr056b+H+x6AnaBB7qkw9LNqhqBLgIZ2jmAuCfqjpfRK4WkdMARORgEVkJfBL4o4jM77peRMbg9BQ+l5T0XSLyFvAW0ABcU/D9Gz2KgqOdb2lk+oofMn3Jdzmu8QGmj3k3pc4kjoYq9F1NjJyeT1qmgx6ieNLTJyL1wJHAnwFUtUNVN/lhsjl9RcCLZRRyqbwUEg49lY0mKjviXcUq5n4zDB8otEzHYsrMi+9h5WsrWfnqSv5zyT/47uh3U+aTb6UnUxRRt1SqNnmGT4uzq+ojqjpBVcep6rXxfVeo6kPxv19R1ZGq2kdVB6vqXgnXLlPVEao7DnlQ1WNVdW9VnaKqn1XVbYXdvNETKWiB9PnPQEcroOy7dTa3L6zLqFdeOXxp7TE8ec6uyF0nG0RkbsJ2YUIquwONwF9E5HURuVVE+vhhrjl9RaIgYQHqhp7IlYeO7U4rUz6FYKKSmYIFRdWcPiNQdJXpXIZ2JtPaGmXz5kh3Q+amLRGuWbBbxnwKtTWTPblQqbqW6/eYFsW3xdkNoyJJGIWjQJVG8k4qky51vdPJ55gGpqc4jl+OvXzOD2SmgFdVwAHATaq6P9AMTPfDYnP6AoTbF9HN+bmcm6o13kQlNYVUriQWdb0ZRhDp06eKo49uoFevMFV1VYw6cBS1fWrz0qZcKz1e/Jj3BH3KC38CuRhG2ZPynZp0BDEJEUOoGjCUz0527/TlOoUnnQamG07qVb5GUVgJrFTV2fHP9+I4gZ5jTl+AyPSitq17grZ1T+zwkub6w17ofJpCRMVIxnr6jMpi+PA6OjqiaGeE+qH9APdOR77Bp/LVJiMNsTw2wyhTCg1698jMp/jd/XO5eeYbXP/XR9m2ZYuv9qTTtkJ0sKfMe/YND+b0qeoaYIWITIzvOg54xw9zzekrAm5epHQ9SEcN/x1HDf+d77akW37Bq9akHi8qijl9RmAoVJvuuOMD/vrXFUSjEI3C2/+Zz5v/fNMXe9JpY6oW61yHq5o2JZJHL5/19Bk9jC5dWfjqK7z5/CzaOzrZ0txGe0sLN177i+4Kvl89aX44fonnmw66xLu5zxfjBKd6E9gP+KkP1prTV26ctdvN3X+bqJQj1tNnlDeJ7/zvrl+84w+awry/vlqSZRMKzavHa5OCxsT1Zhg9iS6NeGfO7J0Pbl4DbVtdpZfvyIRcr+tpdbRyXbJHVd+Iz/XbR1VPV9WNfuRjTl+ZMDf2FA3VtTvs8+ul9zv9chaVREHJ6x4UW7LBKGsSy/2wYXWI0N2aWSMh+g7tV5Qhlrnk4XYYejkv5l6wNoH19Bk9Cjc6lfxejRw3nnBV1Y4n9RkAtX1cp52PPelszCfdVJpRLhroie7li0dLNhSTquynGIWSz9jqZBqqxzE3tpCvzpqZV9r5np9MukAK2Ry5VOHavQrK4Cdei4igiPXcGQHBbaUnmZtv2p+v/N/rLFvcQjRazYDRAzn+Jyd2p13onBkvyCUoVaIWlUujlC82mjQZRjeZ3rGDjz8RBZa8/RbRSCcra8fCkRdAKNx9TjEa5nON85Ap+EtiY1fRlzwod4Ljz+WEOX0BIJcf7w9W/LagF9FEJUCY02eUCdm0adSo3jwy83AgtTZ41RvnBen0KVVUvB6pT4r13BlGjkgoxKEnnsyhJ54MlE4z8u1RTFVHg9QN9EYmysvrs+GdZcRPDpvETw6bVPR886kIpasAltvwKc+FzwK5GBVILFaaH75EbfJCn8ppPogvlTJbp8/oIXjx/kQ6O9nc1ESsDH+j02lmYj3NyAHvArkUBevp85lKaS0pZHx6phalSnk+htETmTdvM+d9eg4bN3awz1mrOPrSoxEpriPg1SiGQtIyjFKxpr02+0mGa7I5PRvXreXOGdfS2d5Gv4GD4DO3Q12/4hjnIek0r1y0sOQOaoDm6+WCq54+Exd35DN5NhM/emkhP3ppYaFmFZ1KEBXvsOidfmDa5J5c9CkdsZgyf/4Wvvq112lq6iAWg7cfns/aBWu9NtNXykWD0uGt/bZkg2FAZu1TVZqWLebpe/5G67ZtRDo72di0Ad58pIgWGoHBevqMXHDj7B0w4mcAnPP8bX6ZYxQLBUo0FM4wciWbw/fpz7zC7NkbaGnZHllWRIh0WKTZskWx4ZplSi7BioztuFmWKpl5D93DmgVvsX7d5u07JQRVNV6ZZ7ikZL19AYvMmQuu5/SVe8toEMi1cI4adQmjRl3iszVG0QlAT5+I9BGRO0TkTyLyGc8zKAE9PkCQCzJVerLp04KFW3dy+MLVYUYdNIoR+wz31E4jO56W+VxbrQPYgm0YuVDI+xKLRlj11mtEOzvYtb6W2uoQIDBiCuzzUe+MTMAcecNL8grkYhWr7Hj5jM55/jZfe/mCLipBt88d/g3vFJHbRGSdiLydtP9kEVkkIotFZHp895nAvar6ZeA0b++xtJg+ZabQ59O/vnqHwC31w+r5f/f9P0677lQk1PN6iipJn2xxdsNIj4TCVNf1AqA6HGLcyF3hu0/Dp38HVf5NMQi6xgTdPj9R1Zy2oGDROwPM12c9Sphqrp66B1dP9bciG+SIdSWfqOslXcM73W65cTtwcuIOEQkDvwdOAfYEzhORPYGRwIr4aTYmz8iZ2toQ++/fn3B1mH679OPUGR+nb0Mf3wK4BFmboMIaGWxOn2GkpbW1lTVVg2iJhWjuNYRlx15ZtGGdQdbBIERlL5UOl9na7Ob0BZmbjvxY0fMMuqiUMn/P0JigVsgXAAAgAElEQVT7DRpEZG7CduFOyarOAjYk7T4EWKyqS1S1A7gH+ASwEsfxA9MBwwWfPGc2r7yyiWgkSqQ9wqAxg4qSb6krFZkotW2e6FPXnD5bssGoUAp9T/7469/x9H//xwOzF/HQ86+xKdTXI8tyI8gNTEGYXlGK/Mutp88CufiAFwXvgBE/gxAoUWIa8cCq3Cj1S5uJUtvmjePpqucukfWqelAe141ge48eOM7eocD1wI0i8jHg4XwMMsqPQt+hdWvbWLBga/fnztZOWppaqB9WX6hpOVFqDUhHUO1yh/XclTsW0MXBj/exo62N+a+9QSTi1MckFKZ97TJqGkalvcYWOq9synFac95On4mL/y/ykTWTAdhr1syM59l6d2VE8aN3pqrFqao2A+cX05BiUuqe4VLjlx784NvvEEaIqoJAXf86eg/qTaQtQlXdzj8nVukpM6znriLoyfWzXLUmWZuyPasXZj5EfW2I9nbnZ1yAutFTiLZu5bLJH3LdsgmFmJ2Wrt+yxP+NgKAEqhcvFwoe1tVTC6AbYUkUEzcifO/W+7l36/0p0/GSxHvpiT8SbvCkvPs3py8VK4HEpsiRwOqC7E9DPCJo2I+08yEIw01KgZt7dvu+r1rZRo2GqSZEbbiKKafsxc0f/SN/OO5GHv/J477+APZ0Rz4bhZb1XOemBHWuimG4dfjcsG3TJob0rWXUoL4M69+H3Y/5JItnnMk735vG5ZdcSrStueA8smEaGDz8r8J5iydzeXpaxSpfYXHzsk6uHcSZfc/gx7Pfd5VHoZio+ItqvlP66C8it4jIqS6zfAXYQ0TGikgNcC7wkBf3IiIhEfm0iMwUkXXAQuBDEZkvIr8QkZ4lDGVEqkakXN77i7+1O9VVIaoIIRF49Z7X6GzpJBqDBY8tZPWb29sT/Kr0GD5igVyMCqFL41Jtqc7LxsEnnER1TQ31vWoY3LeW9XP+Q2TTWohF+eD9pTQ9d9cO53utVYk9k6WeQ5yJoNrlF+XWOGZz+gok1wLu9kXY2LHjPL6+u36CbWse3Ok8r5y0xHSS/w4aZT/EIb9mn82qulPwlkRE5G7gaJygLyuBH6vqn0XkIuBxIAzcpqrz8zEgBc8ATwE/AN5WddxTERkEHAPMEJH7VfVvHuVnuMAPbfroqbvwla+9RhSIoUgb3QOINaY0vdfEiH1HAN5rU+IQp8RjhofY8M6Koqc14qbSCK/YdfRu9B0wgE2NjQBs27xjzLSRjS/xrbHH+KJJiUNQ/RjB4ledquzrahWIZ05fTxIXP4Wli1mRJwA4Mv65oXr8Dnn7QeKY8aCKSqpWumKKSsHfuQLer7XuJK16Xpr9jwCP+JDl8aramSK/DcC/gX+LSLUP+bqiJ2lTF37eb21diG0RZ5WPag3Toc7fNeEQ00/o8KWgJVd6vMTLBq5E+4rdcOaFNqn13BlljN/vWl3vPoDj9PWurmJzpMM5ILDXvvukvc4r3fLq/lJphR+NaZXu8CnlN6fPevryoBgFuatX78ARF7Bf7wEcV5d5krBXFYygi0pFtPT75PQVm0SHT0QGAnsAdQnHZ6VyCo3yRUS4/LJJ/OAyp7O4kyh1hAmLcPanRrDPfv15ZOnO1/WUSk/ideWnTbYEQyXS0xq9/LzXUz73Rf5yzZUADB/Qm7bOKO2RGHtMmsirk7/O60uLs2ZfofhZJkoZcKYU5by8XD6P1+cqvx+5/ClG4eo37IwdPpdTYAo/7SzVj5hneWoeW4ARkQuAWThDSK+K/39lKW1KplzeGy/w+9049dRh9OrlxOtRgbZQlBaJ8PcHV/OjucN9zdsr/CwPZV3Btjl9FUlX3aGc6hBuKcZ9DR01mt79+gFQFQ4xbmg9grJ48TJal7+T8dqg6UKllYPSPF8lprltQcHzRZkrXVyKfU971g3gH1sf4R9bcxs0FSRhqaTv30uHT2PieiP/QC7F4BvAwcAHqnoMsD9dY2ACRE/QJr/va+PGDqZ/bz5DBtVQP6yecE04PizQ6QVcv3h92muDpE2ZKMey4dWzVRXXm1F+lGMZz0QxgptEIhGuv/4W5q1ro6Z+IO1VvVnauMVpl4100PzubN/yLjcqrXylpQyjHHvu9FUyxYya9MkxN3Nuv5N5p20TE8PT+OFLmVuRehplLSqxPLZ4IBdVDeJC6m2q2gYgIrWquhCYWGKbslLWZSiJYqzL1d4W46SjXuTRR9ayamU77RvaGDN1N6pqq7qDuQzZY0jKa1W17OY+5EvZlquu+cbutSkrInKyiCwSkcUiMj3F8SNF5DURiYjI2UnHoiLyRnx7KGH/WBGZLSLvicg/4tGJjRwp23KaAr/rZKrK1772XW699W/MeWMBdz77NsuHHEJbvMhJVQ29xx3gqw2J+Hm/5VguStWg2DWnL5ctKPju9FVSq3qxC1ZMQYmhlTIJLE6llAejm5UiMgB4AHhSRB7Ep7UAvaZSymIxGqNu+91yPlzb1v1Zo8qk4yfxkYuO4MBPH8B5t55L70G9d7ruD79fwphRjzNpj6dY8coKT2wpxpql5ULQe1Dja3f+HjgF2BM4T0T2TDptOfBF4O8pkmhV1f3i22kJ+68DfqOqewAbgS+5sClQa4qWinIs76Vg1qyXeOmlud2fo+2ttNUNYvg5P2TgtLMYfcFv6LvHwTtdt3L5cr7z5a9x4Tmf4du3vuiZPT1tnmYmSv0cym2GTtF6+kxccufiyX9kQI0QEjjn+ds45/nbcr7Wy+dc6pcpKHj+HGLifgswqnqGqm5S1SuBHwF/Bk4vrVW5Y9qUnUhnjLtuXklvrQIFiW+jDx3Nfmfvx5EXHcmgMYN2um716lZ+cd17dHYqW7dGeOxHj3lij1V6HDx/Bv7M6TsEWKyqS1S1A7gH+MQO2aouU9U3ybHvUEQEOBa4N77rDjJojq0pahTCjTfeutO+kVNPYdBhZzDyMz+hfspRKa/78dU38OGqVTRv28bah39L25rM6y67Jd/3P1UjYTE6aCpRs8ttTp9F7wwYMw74AxBmU0fuhcTvoAT5pJ/q5S7HyrX3lSq65uhVDPEK2GeA3VX1ahEZDewHzCmtZYZX/O+hDbQ0R+lNNVWEiRJj7EkTqKuvy3jd9e+PJpLgGEQ7o36bmhOl0CevI9r5UYHKU5saRGRuwudbVPWWhM8jgMQu3pXAoS7Sr4unHwFmqOoDwGBgk6p2LWi7Mp5POmxN0TRYA0pmliz5gIULd3w+9aMmMHDM5KzXRlo2d/8toTDRli2e2ZW8fqnba/OlkEjxpYzs6QfO8M5SW+GOojp9Ji65MaJXlJnblnLhrAcynnfNYc4ImULm+2V6AQtZG8sLUfEirULxfq5U3j13/UXkFuDhAM7r+wNOC/2xwNXAVpx1+nYe7xJQijEnrlyJRWM8eOUyBsRqWBlqRnB+7A7+2tSs1/Yf0Z9JH53EwpkLATjq26lbxFORrXJRSAUiX00pxFn0Wsc8L7P5R+Ncr6oHZTieKlE3VaXRqrpaRHYH/isibwGpas+Z0iyLNUVLhdXN0vPzn1+/076pF/0qp2t3Pe1SVt39Y5AQdaP2pPdue+ecby7OVSm/t3yXuPFDB0tZdoM0Xy8Xit7TV+ovKMi89alvMLRuDM+u7c0RNeO5dtoUhDATQ4ewIPoCsN3Bu3balO7rrjlsT9eOXzpBSVWJ8rt1Jlt5KPYix6nwtNzmV7HarKoXemOA5xyqqgeIyOsAqrqxXIMqmD7tzMyvL6PvlggxlN1ifelEqdtvEL0H7zx/LxkR4bjpxzH1y1Opqqmitl9t1msyVSaSdagYLce56FOptanLjkLxKRrnSmBUwueRuJjzq6qr4/8vEZFncaID/xsYICJV8d6+jGnmsl5oT19T1LRvZ+6/fyYvvLDjgJWq3v0YMjl7e+b0se8xg9PoM+Fgots2UTdiAhLKPI00m/alOhaEIZmlrqMl5lvUMhywyJy5YNE7A8J/T7kCgLaoIwr3t6Se+/LTafvy02n7AhCSakLxxsmufdnIJeBDUIdmlvoHyZMx7/kv2RBkOuNBERRARIZQMUvQ92yiEWXRzM2MqAoxpk7ZIh1U18Q45JqjM163+o3VLL/uX7x868tEOiL0GdzHE4cvFUFYc6/U2gQePQd/5hu/AuwRj7ZZA5wLPJTlGgBEZKCI1Mb/bgAOB95Rp3n9GaAr0ucXgAddpHlIPGLokSJyZC7XGT2P667buZfviO/cnPGaWHsLTc/9nScenkmsvYWagcPoNWqyZw5fqs9e0lVHzCc4WMXooAti3eEWM29Bweb0BYimrfXU13Sw94Awz3eMByBEiI00MZnDOev5P3U7d5e/+DbgvGT1w85Om2YmMglLKXv2sl1bage04BbRWMW1tVwP3A8MFZFrcSpiPyytSUahqCq/OXMJrRGhiioOqRUOrhE2fGcq1bv0TXvdphWbeOCbDxBpixCuXUfzumaOu+y4rPnlWunxs1W5UG2CMh+R4ENPn6pGROQi4HEgDNymqvNF5Gpgrqo+JCIH42jIQOBUEblKVfcCJgN/FJEYTiP1DFXtGtbyfeAeEbkGeB0ngFRGROQCnHVFRwJvAFOBl3CGphtGN3+7426am1t22Dfm2E8x4qD0WqaqLPnd+bStfpc7gbphLzDue//EmfaeO4n1nHTvs5d641YzbDqEg83pM/KmPVJNTIU+VZ10xHozRXanNRrh/dByBse2r32l7BgIYcbSPYjFR6XcfORp/N+snBpQge3CEjRRKVlXfTHQ4EfjdEM8iMss4FXgOJz5O6er6oKSGmYUzEv3bWLh/5pRgWEaoioSYo2EGHXkyIzXrV+8Hgk7ZTzaHmXVG6tc551rpccrBysffao0bfJzEWFVfQR4JGnfFQl/v4LjiCVf9yKQciKUqi7BiQzqhm/gzDV+WVWPEZFJwFUu06hY0g2nrrSyno3lS5dxy01JbQgiTDn7oozXxdpbaF3xDsScelrrykXE2poJ90rfSNZFqhgKyXWwVN+Dl8H28iG5fFREw3zOBGsNvlwoidOXqkD0NFFJ5L5jZlAbFja09WLmqr6M6au0RjsZXFPDiPB47mt+Etg+j+/yF9/eIYiLSIjBVbvnlFfiy5D8YqR7UUotKqnSLbWoFIJP82ZKgqqqiDygqgfihEEve6zi4/DU/Y10EqNWhZXtCgKjbphGuE/meBe7TtkVAAkJ4Zow448Zn1N++VR6CtECL77PIJaLgmyqoAapNLSpapuIICK1qrpQRCaW2qig0vVuVWpDRzpeeOY5eleH2Ziwb6+zL6b/yMxaE6rtTXX/IXRuWgdAdf1gQnV98rbDrwZwLx2+dOmXso5WzHJaZj5fcHr6cp2/Uaks2jSAftVO9Ok1rcLgmhom1ce4c8MSGqrGMjg2hIWxl9Je3xRZ0j3kMxeSxdxrvBYVLyt7JSe/4Z1Bjt75sogcHG+trzh6asVHhihtKFGUMMLGybVMmjYs63V9h/TlvL+cx3tPv0f98HomnDjBVb7pghQErdKTWC6C1NKd/71JRTVIpWGliAwAHgCeFJGNuAgqY/QMqqqrGNS3FkXZ1h5h4OAh7Pe5y7JeJyLs/q2/sXbmjaDKLh+/2PXQzuTgVOX6e1MqHSzX5yUiy3Ain0eBSJaIyHkTGKevJ3PP0gEc2hAhqkJjRwcHD6qmOSLcuWEJY2OjmRP7L+tZirB9InBitE43zh7s3KJeri9JOYqK5r9OX5Cjdx4DfEVEPgCacYZ4qqruU1qzjEKoawjxVvVG6rWaZumk88MQB7d2Ut0re2T7AaMHcPD5+a/Y4Welpxh6V47aBPgypy9IqOoZ8T+vFJFngP5A6qhpRk5U4nSMhqFDEBEa+vWioR9QV0Wss4NQdfag1DWDhjPqcz/138g4ueqM37176fL0cgh+EBv61dsgLceo6novE0ymrCJKdEVPDOIXXwhLO7fQKxzjtQ11nLBLmF5hpT0Kh1SNZwutgDOXr6FqLA1VY3dYrgGcyJ1dwz1zJd85dtnSzCfiUy7pBoXCbckjOl7wh1ydAozDCYZwKvDx+P89hkrUpv4DqmmXKI2hNlokSkdzJwsfyn9N0FzxQ5u60i30/XUzNaHYPZOFO3wVGVl4B8ThsyJyhao+hxPMZb9S21WOVJre7YCChBKqx23bWP/s34pqQi7vc6bnnyr4lds0CjnXS7rsz3QfftQ9s6FATHPbgkJgnb7EoTNdm19ORal5ZdX32NjhdLou2RZmVUuI95s7aOzoYGHsJSaFDgOgMbKYxsjina6PuVheKJeJwYlkE5VsL6PbHwa3ouJFWcilXHlZ7lTF9RZkVPWDxA0nIMP3Sm2XXyQPP0x8l4KkTW0tUf5z+1oevm0trc3R7Bckcd7nRhAKJbZkKi//eTaxiH+rcST38HkRTTP57y7cOur5zEUvtDzkorGelruuBdrdbOXFH4DDgPPin7cCvy+dOeVFpjlcQdI+VWXN8g9oWvNhXtcfdNhUkkt24xN/QqORwo3LQKppLF6ml0wpIn/6kX7Jyp46ZS2XDWgQkbkJW/LILQWeEJFXUxzzjEAP76zYVqQUzG4SeoehsaODaYOrWdNWxaCaEFtaHYdPCHPZi/NSrseX78LsuZ6b7nvItD+fl7AUAX6yzRX0PH+lEpdsQET2Az4NnAMsBe4rrUX+Ug7a9O3TFvD+282own/uXMfNz0xxNb+kX301R5w4kGcf2wDxuX0d29pYs2ANw/ce7rm9ubTmuiFVw1Y5aVM2/Mg/6A1MHnCoqh4gIq8DqOrG+NqBRo4EXftUlftuupElb7+JAod/7DQO/5i7gSf9+tdz4GFTeeWFF7v3xVq3sm3RS/Tb8yMeW+yQj/5Vev0sF0ptg4tALuuzzNM7XFVXi8hQnPnGC1V1VsEGJlF5tc8ic9V+N3HVfjcVnM7buoThvWBcnxpebOqkORqhLQpNoUYWxeZ0n3fZi/O47MV5BefnBjeteF7+IOSap59RQv2gUoZQicgEEblCRBYANwIrAFHVY1T1hhKbV1RK/cOTzPq1bbz7RjPRiBM9fMnbLaxb2e46nV/+eh9CvZVWIrQTJRaL8eA3H2Tev+YFIlR1LtrkR9CqoH3fnqAuhpuXz9DzZDpFJIzT/IaIDAH867qucILWwwewYd1a3n3jNaKRCLFIhP89eH9e6XzuK18iFE5YUF1jLLvrStbOf9kjS70h12GgyaMa0o1YyTUfN/WzfMtI0IZzpsJFT1+2dFbH/1+Hs2ap2+VocqIsnL6gfLl+s3ib89szpKaGObH/MieymHGx0RwaOmqH866dNoVbjjydW448vftzqUklGoWKihvyuTbdkK9CbcmM+6GdAW59X4izNt+pqnpE3NFzP46wjAmqNl319Xd3nGCu8IevLnOdzuDBNfzt7wchYWeYZwihs7mTWb+exawZz3pmr99k056gaVO263xrkFL3W5lxPU6FaqiIXAs8DxQv6kYFEUTdA/jLTbfsUMnujEZ57EH3Aa8HNzRw9a9/vsO+2MbVPHfFJ5l318/TXJU/hTRKeVX/8eOaQshUb8y0v1jzTRWnxSiXLRMi0kdE+nX9DZwIuIvQmCOBd/qCKixdHDh4AwcO3sDfP/KLgtJZsuIXNEcjPNW6mMaODiaFDmNwbAgjetUwdXA1Z/Y5gaePuISnj7jEI8v9J+iikkkcMglGQaLSNbzT7RZfskFEghQg5SxgDfCMiPxJRLoWZzdKRCwWo3lLJ6/M2uzsUGcTYNGc5rzSPOCAAewypI4qCSHxr1dVWfTwApoWNnpjOMWv9AQxj0SKXeFRKm++cSLijG2ehTPf+GfAh8DpqvqvkhpmeEY0GuX1Oa+yYsM2OqMx2iNRPli/lQf/kd9XPGb8OOr79+/+rAqRzk7e/udvaFw41yuzgdLoS775B7VeXoqhxx719O0CPC8i84A5wExV9SWqcOCdvqDT0HcLWzq8mRLw9KqLWbLiFywNLacxtJZh4X6sau1gcK0z1LMLIUwo/g+2L9buNoKn31SaqJQ4WtlmVb0wSGv0qer9qvopYBLwLHApsIuI3CQiJ5bUuB6EqtK2fAXtj15K+2MXcNVhs2jRCNvopFaFOoR6rWJLde4BnxKprQ3zyGPTmDixb7cTiUKNhNi2eqtn9xHEd76cCPpcq1KjTs3rAVVdqKq/V9UbVXVBqe0qV4JS3jZt2MjPr7iKL531Kb7wibOprQqzpbWThR9u4t01m2ntjBLpzC8Ai4gw/dqrGNjQQGtHhKWNW53xE6qsecPz6Va+kaytha6Lnc9376e+l7Ju5sWICFVdoqr7xre9VPVav+wNvNMXFGFJxXUH/p5/L96dtmgVbVHvYuKMi40mRIjmaIT3Q8u5fs0HzI0sYWloJUtDK3c49+kjLkkZ3KXY9ARR8YpKmdPXhao2q+pdqvpxnMidbwDTS2xWj2Dryhae/OxLLPvzHYSHrOa9J/dhXVNfagixWTpYHWphEx2sDDWzqZ/7OX1dDGmoYf2KTqoQ6gjThyrChNj1oBEe3o1/eK1PhebvJZ5WeCp/Tt/LIpL/ApJGYIjFYtxx85/4+me/yLy5r9HW2obGYuxSX0ff2h3rY/0H9E+TSnZ2230soQFDeW/tFlo6tjuPDRMPzDvNZPzWHz/WYy5mHS3byKuSkWMvXxDmv3cReKcPguH4Lbvg//HSadNpu+EjtP7mKJ486Upnf7Py6OpeXL1iZeYEcmT63jcDMDY2miE1Tg/i4NgQBsYGszbyLmsj7wJwwaz7uGDW9gCJIalmcvhw/n3El9On7XPFIwii4nVeXrcgVfoQKlXdoKp/VNVjS21LMShGGY11RGl8pZHN727u3rd5XSePf/tdXv2/mSyYo1RVRVGFaCTM+Fpo0DqGaB1VhGgKtbM51EnvPvk3TD07swndBhGUTmJ0EmO3k8ZTW1/rxS0WpdKT6m8v0ssVrys8nrdua2VrU5xjgJdE5H0ReVNE3hKRN0ttVLlSqrrZ3Bdf5vwzzuGJh/6z07GVG1voW1dNQ99auoIVh6vy1761axt59tkXdthXO2AIw/Y/Ks0VwSQI9eh8ybQkWCkpx3X6Ar1kQyJ+OBW5sPwrXyRcFQEcB2zT4pF0ttdw5+Jh1IUhHgSMZSt+5VmeUwdX82JTJ40dHayPLSVUFWJd5H0uf/Ht7uAtjR//DgCzY8+h8dgZC2MvAdsDu1z+4vZ5oH4+u1TrYZV6EU8v1vvz5R5UKnLJhp6Mn9oU64zx/FdfYNuKZogqLWOGMvz03bnrgpUcs9s6Bu3eigg88sBhXLDvCsYf8Q6jX9qT4+jHllAdL4/rRceitSDwwxkT87Zj0dvb6K3VRONLN4QRdjt8jCf3WCxt6qLYvXzJeZViLdKcKT8nzi2nlNoAIzc6Ojr41VXX8M68txgxehTfuuJyRIQrv/19Nm3YkPa6mCprNrci1XVQVUttCI45Of/ZBu+8s2innpo9z7wo7/SKTdHqNjmSjwYGmR0CppUBVvuMs/j8C1l8/o7rIa75xnnb/24cwqMfjGH2W3uzZu1Q6qvhteZNvBJ9j1ei3v0o14WUsChbaGULrd37P1p3AgAhQnxm4DjP8nNDMVpbgiQGXi+U2kWlDe80/Cm3sZjywC/fZ+PirYQiMUKqhN5bx2+/sJTaTuiMhOnVp5WJu62jcVUDV3/5q/zp9s/RdsPhhP96An0ePZtjbzmNX90yhc9dMIpCfpsOOLw/iNCPGvpTSzUh5lzxLO0bWrNfXCQyLcCe6bxCFmevLCoqsnBKVPWDxA1nOPr3Sm1XOePXe/Hvv93N26/PIxaLsWLZB1z6/y7kWxd8Jb3DFwrTZ+JhDPrIuYyffi8TfjyTXT/xLT73lQs4tgCnb/DggTvte+ueX9LZ4t18Zi/oafUz8K+O5oZyi3JsTl8SXVE47zj81wDUD93A7Hem8L+Vo2lqF2qrIjS313HSiPVM7TcAgKbIEvoNO6M7jUJejCvnfZWtkRD19Open29d5H0eaXuSW448nXWR9+lUWNUCb27q4NDQUexSNQHV+HIPVePzzjsX/F6bxauX1st0vA7LX+kVK6NwXp+1mU/u/Sp/uGEV4fjXH1PY2hkiCkRR3vuwPx8sHcbwUY2ccdpbjL3rQDr+8AU27Lo3HaOGQm0VTQsaufTCBdx243Iu/sKbnPHNda5t0fY29gu/xJ6TWuggRgcxIih0xljzzHJvb7xIVJI2eZlmD5jTh4jsJyI/F5FlwDU4S88YAaFp/Xp+dtkV/Ofe+3Y6FoumCn4vDDz8HCb/7H/sfsmfGXHuFfQa5QS1W/f4Lfz55tv4+mfPZ/WKVXnVzfr168cuuwzZYV+kZSsLHvyj67T8pqfVzxLTKlUgMJvT5yN+tTLcOu03RDqrWL5mVwC+tafzMjeta+DR5z7Cy40DWdbsjIStDkXp16uFySOXc974Dzih14TudM4Zc7MnCwH/av5XeCU+bHNi6BAmhrav0dhQNZY/rH+KuZElzI491x3YJSTVCGEunPVA3vnmSj5ROEsRmCWfxd0zXeNJ+esa3ul+yYZAISInxJdp2C/++cJs1xjZWbe6nS8dP4/vnPkOjWs72NJZxUOratjaHmZ1S5j7l/WlnRjriLKhPcwzb4zh+reP5YkfXMPGcTuXzzWvrCIaDz4Qi8Lbf5/Hp3/prndu83230/b68/zyjP8SFiVClGpCVGuYyOoWT+7bK/KNEFzsIZf5aJNftnSTxxp9AarLZEREJojIFSKyALgRWAGIqh4TX1/UKAAv5p22NDdz5823csnnv8Tbb8zL6ZrqwaOYcPUTjPz0lVT1HbDDsaZn7yK6bQPa0UpLczPfvehbfG2Qu2UWZs58grPP/iLr1zftdGzd/Nmu0ioW+dbPihUB020+mXrzktdVLnZPpM3pKyH5zKu5eerv2NpZxcBaWN3UQHukmrlN9arZM6oAACAASURBVAB82FrHu2uG8+v3labQYvZiLAcOEh5fOZxDGjYxrKOGvr1aOX10I5M3HcfL651Hee20KWz98H7X9ncNJd31d3cDcPXIo5m5WtlEK3v16k9T+1gaI4sRwoDT+9dQNZa1kXfZpWrCDmklf4bSzYlMtgFyF6V8bS70xe8ac+7XM6uQ4ZpfA84Hfigig4D9SmxPYMml3EejytWXL+SJO5qo7gxRhSAihBFe3VTLh5tD9KGKNlXaQhE6VVlPJ80SYZ9zp6RNt2HKLs4iferMPYgQ44XfPkfj549hyNDcgrBEPlxBTZUSDsO4qmrWdYSJEqO3VlGrhTdIBEGb3NpRiM2FrkfoV8WmK8hUhbIQ+B9wqqouBhCRS0trUmWR+E4kl9HEz8nvTSwW439PP8Nf/3grrS3uGpGGnfldagenjiAc7tVvxyHtnW38+NLv8qs/3UQoHM4p/euuu4GOjtRL3UQj+S2BUwrS6VUqpymTtqX7XvMdMp8LuWptyZZsKLM5fWXn9OUrLIlcvs/NtEaECY5/x7DezcxZuyvXfvgk39v1RDa2OxWZXtWd7NdnEJs76mmNxljVGqImBOcvepnLh53A2H7b2HPXVURVGFTbi4uX/K/g+9t85cfof+VMPj7lTfpWT+ZbS+ewqBliGhcYgWOqT+TZzqcBZzjn2si7TAodljK9a6dN2SGgi1/4KSp+T0TORVS8yq+AilV/EbkFeDgga/U1quom4DsiMgPo8WHQk8tRrvoUiyl/v2MF996zmn6RWiIoNYSIqlKNEJIwus9Qdr/8MJ7/8qN0bu6gXZxhTg0H7spup+/cyNPFrgcOZ+zpk1h433wi8aGZvapCLF3anNXpe+P3H7D0vvdp2Hcg+45bBwpH7fsBj8+Z2F2Zqi4TJyHVO578fSS2QOeiORVZ4SmT7zMPzgLOBZ4RkceAe3CaQwwPyVYOU5XlO2+6hWefeIrOztROlPTuT/2+J7D5pXt32F87YiL1+6QPEt1wzOdY/987iTZv7N63rnEjaz9cw7CRmZeaWb16Df/854Ns3Lgp7TlSAaUnXSNStt+yTBTb+QraPMOgU3ZOHxT2JR844jpOHuys2VITijFrnbCseQhbO2FCaCp/Wb+ExshiDgodw4qWXRCgKS5GbzRvY2x1PafUnsSL65Xd+ykbt/VlQF0LHdEqjq5ygq3sPuq7LFnxC1d29Rq4pfvvNz55KX37Qm0oxjeGHA/A3RuXdh9vjUa47EVn+MO106YghFkUm8PR1ccBECNGlE5uOvJj+T0kDyl3UfE8L827p2+zqgZpCOXMrj9UdbqIXFxKY4JCLuUlsdy/8HwTX/z8q3Q2K2FCVBOhTqpo1ShVhOi1ax/2veoIBh3gDD0/4f6zee/2N9n07gaGHTWaMWdMRLLUPo74/lFsWLeZ5bOXEw6FiUWU8097nVGjenP97XszYUrfna6Z95tFvH17I4d89FUuv+NwfnXT3dT3a2N883Lmzp1Ee0ToRYiqsf1cPqHik06Duvanc/AKGf5ZDLy3qXKDRqnq/cD9ItIHOB24FNhFRG4C7lfVJ0pqYA9k8cJFPPrAQ7z28pzUDl+4mhHnXcWgw5yI5et2Hcu6/9yAxiL0m3w4o7/0GySUfqRBqLY3e/zwQRZdcQLa2bVGqfKlK//Clz+2P2eckbp+tHz5Ss444wt0dHRktL+uf0NO95kLXo12SKcJmdJOnt7idopSMXQyKKNBdqKMhrh3UZZOXzayFY4XmzrYv38t16xYTWNkMUOax7MXYxlX3Z/3O2FgqIFjh1ZTF4rw2JoYc2L/3X5x52E0hRoZHBvCvR/058ThNUwZsoa+Ne18akyYK5atBeAXB/6e77769ay2vvyJ7wPQp248gxqaaFznTBh+Y+k4lmzrTVXIKVHrIu8zMXQI42r60xaNdS/IfvmLb3PrkWeyp4zJmpcXL04pRCWX/AuleKIiqAdD4kqNqj6YtOvmkhhSZiSXsW9e8ibNzVEECCFskQ4iYRj/kd0ZeexYRpwwFgltr4hX961hz4sOcp3vqb88lWUvLOO9/yxk7awVoMqKD1r5zDGvMm3aYC67aQ92Ge70/K18ejVL/rmCmNZR26uDT3xkEVd84xzOPGcetFXRoTFCEmZrb2HssSMLeh5dlKrSk06bMqWV7vqy1yat6OGdAKhqM3AXcFd8WPongemAOX0+k1iGN6xv4qeXXUF7W9tO50ltbwbsfzLDz7uCUFVN9/6hx5/P0OPPd5VndX0De1z+IKvuvorWdct5f8kSti57nCtmP86NN/6Z3/72Gvbee8/u87du3cYnP3l+VocPYOyxn3Rli99k0h+/RgwUq2Es1by+oDTEBSlISy5UnNOXqmD33fUTAEwMT2NR9EUmhqfx3OZWBtLAwFADTazjLd6nqX0Jp9SexJi+wtZOeGubsIVmhlSNpzGymAmhqSyMvcRBHANAXVgYVNvOys2DCEuMmlCUn49rYOlWp+X8tmm/YWivFj7+9OU72fTqWd+O/8A6Y8v7D9jMW+9NIIQiomxsr6VvVZS3NlUDMLRqHMTg/Y7N3VE9AWcx9hC8o8tYF3mfZ6ft271un2oMEe8cjJ4kKr5SYa3pInIrcKaINAOrgTeBNy1Awo6kKmNbtjhBVhRoJYII1I3ux9TrjvM0bxFh7BFj2bKwiXXPLO+ehaDAvBe2cuoBs/n+L8bTb0AVHTcvIhyK0qnw1gtTOO2sWUyZsJolHw7jfy/tSbgvbIuF2esvxxPuW+2pnYWQKZgEZH/Hy0WbfK3wVLjTl4iqbgD+GN8MH0kuw6tXrOiOON6NhNj1zO8x5NjPe5p37ZDR7H7Jn5lz0/fZ+uZb3fvXrWvk05/+Cvvvvzc33DCD/v3rueKKGbS07OyIJts56tCTGXXISZ7aWShudCFREwvRErdxGgrJJ5A9fVhPX0lJVygGV49lHyawOraVhurxNLGO9Z2LmRieBjhLLoAzP+7R9sc5JnYizVGnMnbMwP70DvdjS+d4ntvsjO/+ILSC3WKjWNHWztMfDmRU7ygT+m9BEYb22czAuhb61rXSuNUZRjp9b6cTZMZb/7eDXZ2RKjqiVSzfPJB5a4ezsqWWtW3CoFrloMGbuWV5B8PCzvCpI2p35/n2JcSI0RBygrrsUjWBRZE5DA2NY13kfQAue3Eetx55JhfMuo9rDttzZ2Et8Plmc/xyOc8NxWrNKaaoVGBr+keAXVS1U0RGAPsC+5TYprIgJrHuieAKxFBq+uUWYCUfJn1yCovvfIuWtg5CQC8cpy0ahZ98+11iVco3R0YYUgP9+7Tw5PxdeXjBWQwYo/SbcQYNl/am97IttL2zkeZXG6kd2tszx6/QdzAXffJam4rR4mzaZFQa9f36Ek0a0inhqqxD1Qth1NRTeO+xO3eqpb/++lucfPKnCIWE5ubmnS8MhRl/4v9n78zj5Kiq/v2cniUzSWayTkJICNmAsK9ZWAwYtoAgQUBBFFAEUVSU11dF/amgICCCwIv6RjZBFAVk82WT3RC2EJZACJAFSEIgk21mMslMprvP74/qnnQ6vVR1V3UtfZ986jPT3VW3bk/3/eaec88953TGfOoE3v/P/Qwauwfjjzht2/PKxOtxnstAc0O7grTqVmkUJRkyqy/8cWYZFEs7u32sifOefaDXyHsnMYd1sdXAtvXtXko+yXo28oc1j3Db2iXc17GI/fsN7D3v5cQTfBBbxsKOzXQlhPbN9bR2NfDix9uzpG0QPfFaRg5aQ1dii109aeSVHD7yeg4feT1PfzCWjZv78O5aKy78lbUNfNIlPL5pET1J4fnWgYyoaWJwvfUR7dDX+mKtji9ldXwpU2PTOab/OL7YfCRH9B3Xe49fHbgbH/cs5FcH7sZPn18AWPv+Mv9G5WCnbl0pm7ndODcUaCSLs78ADAJQ1RWq+pCqXu5znwJHLm3abkQDSZQkiorSb2g/jvzxEZ71oc+ABj57+8kc0DScodpAjQqbJGElGAISceX+j+rZEBdqa5T6/h28llQG3HgGtSMGEOtTw4e/eJllv3mVD38zj9dO/Bc9ne5lsfNan4w2FUaTzg+DoRjZ2qfxHnbboYWBjVscRlJbR9Meh3rWhxH7HsbE488hV/6eDRs20N7eQSK7DqAIU799NVO+eSXD9zyY8Yefyuu3X8adJ+/I89d+1/XQPjM/Cx9q83ALEZkvIm9kHf8RkWtEZEix6yO10pePD5b9bqsvemY2y7RBNChpGV8H9BnPw92PMjk2nfVY6YNb44v4TsvR9K1JsEfXOAbWxZgN7CHjeDO2hCFt44FGtmvczJA+m6mvSfDm6uHsP2I5dbEENQLz1+f+2Ic1bkRE+dyYTupr4vR53zIq31i/mZGNVjz7zFFd3Le8gSP6jmNg/VjWbxb2GthNRzzO3R9v8Uz1ZvjMIFd4p1OPkptL+E7biJoHSSOypy+LWcAzInIT8CJWaGebz30KLJnjb+rFnyN28b0sa01w0NcPYo/j85dfcIu6HZvY9YmZDHn0feZe/B/icWu1sYcEgvBRd4xrFjdaBdgFDrx0GrG6GKrK3PP/TWxJG4IgQGLdZh774UN85n9OcKVva5as4ZxfzaN//1p+/JNdimYYtZux2Q7Vrk3W7CTwDiZDiElr38jxE2DQKMbX1bK+s5uuCdPof+Q3qc9TfsEt9v/aJYycdBSzf/sNute3Fj1/8Pi9GDf98wC0rVjEoz883ip4Cix54k4kVsPUb//Wtf51t6/lB4938atpMerr64tfgJmf+Y0Pe/oeBhLAX1OPT039bAduBY4vdHFkjT6nX8ZTZt/IrGkzWZzK4plO3nJi3xl0JZTWLmv/3Q59YWHHZmYOmMB/2togBkt72lm8Wvnu2FoG9elieWd/epIxVm+wakJMbO5kU7w/u9U001SbpCMeY3jjBj7p7M/2Te3MXz2Mf39cy5HbxelOwMA65avju1iwvo73OoT7ljcA8NjGRXy+fjwD65WOeA2Pf5zcan9fTOqYFDuUl5PP8OM5rzNr2kxaaie4VrDdqbFYqgh5mQI9CIRg5c4pfwFuw9KTbwJ7iUiDqo73t1vBJf2dHjIGjrnlbF/6MOzoMcw4fDSvXPc8b9/zFvF4klpibKCbRuqoIcbEc/dlhyOsSIK3H3qbjlda6U8tMSx/eRJl2VsfudKfnq4e7v763Wze0E0sJsx/o53HnzrE9vV+6JMf4ede3dNySEVOmxCRI4HPAzeo6msicq6qzvK7X9VK7/f3rJvYvOw1+vYfQt9hEwpf5CLb7X0In7t5Hu89dgfzbv45yc3dOc+r6z+IIy67FxFh84Y2Hvne0b0GX5olT9zJ/mf/grq+5Wcxbl04lyd/9nkQ4c2hA/j732+iqWnbzMr5qJT+ZV8bJQ10jD/ZOw9W1YMzHs8XkedU9WAR+VKxiyNr9BVj1rSZtMYXAXBY7ZE8N+17nPvsfVuFQgLM3mydM6SzhWYaGVJXx8iGekY0xpmqA3ihHaY2D+CF9jZWbGyktasP2zVuoiamrO9uoDtRQ0vjRnYdUENPMkZdLElNdz3D+26gdWM/ehI1bEzEmDGih7fb+zCyb5KVm2Ks6+5Da7fQ0gCtXTEWb25jPKNZ2KZ0JuK8kHwSoYZdYpOJpVaOPq5ZDgXCbX51oJWlSiTGpWCriLzdunlOqVQcuBMDsiJ7Z6KZIW+5qv488wkR8W5jmsE1YrUxJl14MBNP34s5v5/DpmUdxNd2o7XCnl/Zj50/M7H33E/e+oQ2NlOLsJE4MeCDmg2M2MedDJ6drZ0kehKoWsXqFyzo4PxvvMY11+5Ffb2z1XE3wqTylXmwc56T+2Titz55pU0iMgO4Fitz2Y3Z4d8iMg34HdZe4FNV9e7U8/sAfwCasbzbl6rq31Ov3QocCqSjCs5S1ddy3P6bwFeAn6Yydu7j7rszlERdHxg3xZdbx2rr2OXYsxh32Eksf/hG1j1xI4s/Xs/muGXUNQxs4dhrn6CuoS8AHz73IPHuTds2JAIuJct7444riXdZ0WUrVvZw5i//zl2/PouaAoXkC9WCLYRbBl/6sZ8a6OeeQgU/9vT1F5EpqvoigIhMBtLegXixi6vW6MvFXYd8bavHE2MHQhJOmv0nIBUK2m3t/5u9Cs599j72Pfr7fNCZZGrzAFZsgjH9kqzc1Mhug9axYN0g+tUmUBX2aPmYWxaOpakOJjZv4q21QxnX3EZCY4zp30mNKOs211IbU+piMd5u78P6zfBSfBFDki0MpJGpQ+pYtlF4svMZJJX1c11sDUmSrI4vhTgcOef/OJItYauFErmUOolwmsq81DYy2yp1ULsxUXJbVCJo9L0mIheo6rXpJ1Q1t/vUEEiahjdx9MWFs9FNmD6Be/85n65EglqE1tgm4jFl2nenudKH5hHN9G/pT/tH7SRTe2sefugT9tn3A845d6ytNtya9FRKn8rF9Xt7oE0iUgPcABwJLAdeFpEHVHVBxmkfAmcB38+6fCNwhqq+JyLbA6+IyKOqmq6a/d9pA7EAranzvy8ilwOTynxLhohQ17eJsSd9j9HHnsOe8R42d21Ekwn6Dh1JLMPYqu8/IOf1Ox/7Feoa+7nSlz4DhiA1tWgiTjK+mcWP/43v9yzlmmt+lfN8L0vE5MOtEHo3MoYGBR/SuHwNuFlE+mMF3LQDX0vVIP11sYur1ujLDHm89KA9egurT45N37ouXw5a44toqZ3ArGkze5/bs248Q+otI7tvbYJFbQPpV5tgQP1mWpraSCZjHNTSyfrN9SRVGNKnm4aaOEOa2hnY2MmC1u24afVSTh80np4kvLKhjfH1A9hLxrHnEOGRNe0s21jPPZ2P0FI7gdXxpQyrHU8y9W+X2GQgVcIhg3Qyl+zf7VCNohI0QQkJw4EjROSHwDzgdeA1Vb3L324Z3GSHA3Zg8pmTmHvb3F7v5vDdhjNw9EBX2o/VxvjCzV/gwR88yEevWiGj3d1JPl6Z23/gxVgNsjaVe62PTAYWqeoSABG5EzgB6P0PSVXfT722lZdSVd/N+P0jEVkFtADrsc//ZbTxIxH5dgnvwRBhahr7UwPUNQ3K+foOB36GgWN2Y/3SLfkgpKaWCUd90bU+7H/2JbQvf491S6x7JHu6eeKJZ0kkEgVX+8JCoeRa5Wir35pY6T19qvoysKeIDAAkwwEG8I9i10cuo0Qp7BybysfJDVZxZDbRUjuBltoJLEw+z5pY61Yhn+nX0oxNjuaIxgn0JIWB9T1sTsaYt7aBNd21jGteT3eihs7uBjq6Gpm/vh9rumt5p72BoY0b2X7wGgBWbWgmJspFo0YyqWUN+w7qIqYx9hyYYId+Mf62bikDaWTFps1MjB1oreoBn8TfpTW+iJ2SY1iYfJ6Fyed7+zW8duetEtakyQ5fjRKlZuwrhtuiUmL2zgEiMktECm7S9QNV/byq7gqMBX4GvIs10TNEjKnfOJDP33gKk844gKN+fhQn//FkV9Os92nqw/QfTKe+Xz31/eppaqrl9C/v4Fr7fhEKbSo9s/BQEZmbcZyb1fJIYFnG4+Wp5xyRCmOqBxZnPH1pKnvdNflCylX1/qyn/uj03obqRmIxjr7yQUZNPZZYXR/q+jYx6bzLGbjjrq7do3FQC0f/5v+o69ucvimNQ7bPafCVMu6LZfXMPM9N7NzXb8OtHFTtHW4hIn1E5IvA+cB3RORnIvIzu9dX7UpfJu8mX6CldgJDki2MqGkimUh9QrX07vtLG0uH1FsG326PXwXAnEMv5G8djwDwmc4ZtDTAuP4JxjZtYPjAdQzqt4Hhw1oRUfbbez6JeC2bOvqyvm0AIsrmnjpGD17NgI39qJEk765t4fnVDezat4H7VlmF2JUEn/Auw2t35pP4u+wmB/O2voBQw7Da8SzmQyRpCcNJs//EjdM+B8COZ57L863NvNjWAcCJz/1vryF42UF7W+9xzuuuZn2yG/rp1iDP1/fsewRFVMrI3tmmqtmTqUCRCumclzoMEWXY7tsxbPftPGt/8NjBnPnPM1m7ZC1Dxg1h3Lhl25xTSiF1O1nuqlmboOTQ89WqekCB13M16mgaJCIjgNuBM3XLnoWLgI+xDMFZwA+BS4q0cyPwORHpBD4C3sDKNny9k/4Yqo/aPo0c+uObPb1HTV0fjrj0HubdcjE1dX044JxfYWObVlH8SoBiVwMhWDpoF5/q9N2PtY/5FcDxNhpj9KVojS+ilUXAgayLrWYPGcc6Xb3NefdufGSbFbRP1x3FQc9czaUH7cFJNTOQvkq/2h6eWjqB4Y2bWNk+kHiyhmH92lm7qR8DGjbR1tXI2q5GVmxsoF9tkt0GrSGhMdp76hhcD/d1LIIYTMLKxrlLbDLrWMPpzTMY3S/B2yute6+KL0ZJMDk2HYAhXzyfPWRMqmcbtupn8qTvEbvnGnf/cGw7gEv16rjpYQqsgGgks3caDK7SOKCRkfuWn749e7JjtKkI3sxflgOZy7WjsAwuW4hIM1aI5k9V9YX086qa+l+QbhG5hW33A+biU8BwVe0RkZHA3ljJYwyGQDB4/J4c8avMbaru5iQwGugi/mTvHKWqM0q92Bh9bFu3bzLT6SROa9Ja5ZsYO3Cr0MlZ02b27glcmejYqq17Oh+hpXsCJ3ZNYK+BnWyK1/LJpka279tJd7yODZvrSSRjbE7W8EFnI6rw4EdKc10znfFapm6/jJqVozhRJvDs+jZiwDmDj+T+9iV8tWVHIMFLq4VdZSpvp/7/+3TdUQB8Eu8Earn7illoXGh7cyx9X96fgfWDaaixnKP3Hvx1Tnzuf7eq6VdOTHX6OrvXR2awl0EEE7kYDIHGjj4ZbfKsZMPLwE4iMhZYgVVXytZmKBGpB+4FbsveIywiI1R1pVjxxTOBbfczbMsLwCBglaquSPXnIdvvxGCoMGZ+Fmy08qlc5ojInqo6v5SLjdGXxU/mvMk9h5zDmlgrn645ijd1CSQtw++k2X/qTd6SDvfMrIM3MXYgAAvjzzNr7SKGt+/MhduPZGF7Hc119fxn1UBuXf8IE2MHsibWCrT2hpQu3VDH0aM/ZNWGZiYMWM+qrhamDRzAmm6oiymnDRrL/R9v4G19ASXRG9oJ8FTPY+wmB3PyiEZO++z9rH9tHP+afTAb4rWMb25j+sgVLFo/mKa6Hvg4xmUH7b1NUpdShMXz8gYpKlXaoVKCFxWjT0QOBF5QH6qTGqqHUic9ldCnSpad8fpeijfapKpxEfkW8ChWyYabVfUtEbkEmKuqD4jIJCzjbhBwvIhcrKq7Y9XXmwYMEZGzUk2mSzPcISItWOGjrwHn2ejOLOAZEbkJeBErtLOtyDUGQ+iI2vysUvcKAYcAZ4nIUqzwTgFUVW1FLBijLwcLk89DElbGWtijZhyp6gjMOfRCWuOL+MmcN3uNvuzafmcNnMGaDVZ2zVXxxfx5WQsjautorI3z2dEr2blpOmOb2lndNYL1m+uZOGgtr7T2Z1zTBj5sG0RnvI5R/Tt4r8P6z3fx5ja6EgN4qOvfW/VxaO1YhiRaaK35BICxfZroiCuvPXcAB57/BKeN/zsbP9iOfz0zjdruBoY1bqQrXtu7OpgLL+tAOR2skR3cKlEK7zwTuEFE3gUeAR5R1Y997pMhgtgNzyy3fTfPr0jdTzfxMPRcVR8ia0VNVX+W8fvLWGGf2df9BfhLnjanl9CVvwC3Yc19vgnsJSINqjq+hLYMhooQRv1LnxcqDSwBH1zex5RzsTH6cpAO97z0oD1YmLTCJ5/qeaw3jDKT7BIJ/7fBCgkdVjueI/qOo6WPcm3r4zz5Hnyx+UiSCiP61vJBZyMAdy4ZTnMd3LJ6NYOSQ9i3XyOruupTRdkhprFegy+9wgewl4yjoU+Mo5r609ptJZdZ3ilM3P0dNs4bSHxjAysW78igPl28v6GJhApdiRqOaziSl+JL8r73Sg/SzBU2Nw09pyt3FVvlIzorfap6HoCITMQSoltTaYSfwjICn1PVhI9dNESMtD5lju9EPEGsJrZVFlEnGlZsD8vlS3fyXJsK6VVFHWAR0aYCLFfVn2c+kS/rp8EQNHLpXy7cmsNVpQY6oJLF2UWkWVXbgY6iJxfAkdG3XZ/qqrecudfvzWkzeVOXbGX4pcM5M0ln+7QYB8A3hx7BDasf5Y52q8be44uhNW5l31QSvW29k3yJ4/oezj5D1vK7twcwpL6eEbX9mNLvSNq3bMGjoQZGNCaJJ5W2HmFQvXLM6GXssvNiujr60bF6ILPf2Y2X1/TnkGHr6V8bZ0NPHd3JWK/B98KhF/JU/CkALnru1a3eQ6ZXyY1NvbnayLy+Ekv3QQoPKDF7Z2BR1YXAQuAaEWkEPg2cAlwNFMrq5xrVpk3VTHri84Md3+Xws5by7mPv0jCggROvP5GhOw11VZvSj0tNke7kPD+LvKeJikOqAK+JyAWqem36iVTW4dCyXZ/uQP3/ZvCWzM/4sE/+xgcfLOPNsV+kX4u1UO7U4LOrgYXOz4WTfgRJA51SwZW+vwLHYWXtVLbOiqykDY4iOF7pC+OHUgrpkM204Zc25g6ac1/v83MOvZA3dQktsQlZxp7FHe1Wps/kSd/jsh2ssM532gYwceB67lg8g+Z6+GhTkoNbktSKMrbpAPrWtVIjSaYOqWNjAupiwsK2JG/q+5zQPI7OOIztH2f95hoU+LAzyXGjuljf1Zf2tQN4+t1dOXDHJUwd/y7771jHB6uGA7BkQ1+WdwpCDbuLvbpX2at+XocYOP2PK6xhA8kIT6xUdRNWGFfFkyOYiU/1cPnSnfjwxQ9ZNnsRqsqm9Zt48oonmffEvq60nW/SE21t8iyRS5AYDhwhIj/EKi3zOvBadpKYsGI0sHp49+FbufPmi6klQbLmRo77n6fpO3R7z+7n5DvlpC5pKQZlkKhUIhdVPS71c2w57ZQUTYU1xQAAIABJREFU3hn2D6kUsjN8Ar21/aykLBYTYweSTH0JBtKXWdNm0hpfREvtBE7vGs/A+jitm/py9PadDO+3gZamNj5uG8TKjf14e30zT69KcNiwGvrXJbi29Sn+a/jhnDexlRdXjWFTwjLul26oZfGGzYzvX88x23cz81OzaW8dzLPv7Mp/vT+bXyYO5fwfWPVkls+aydruBgCG9IHBGwezkk7W1qxlWM0EPom/y68P3neb1b406QlQJT5nr0UlEN/VaO3pCySB+rwNnpFMJLdydSbjybznlktVTHj8ST9eUVT189Ab0rk7sCcwGYiE0ZcmUN8rgycsfuyvJLo3kQBqG+r4eP5zjPv0KX53Ky+VruFcKfzQzFS5mR3JsOFU9Vk715a9py/Xf3Bh+9BykV2LL99rlx60ByTpTe6yMPk8J/adwfqeJE/1PEZLzNpv1xpfxO9aF/HpuqOY0L8fr7dt5oOYVWLhkD4tdCVgpyZlfP8a3uuAXQfEuW+f3RjU7z1isSRTgDHDPubjtUPYa5/5fPjeWPo2biKZjPHyvH24fdEwDhjSw292/BTPfFLPvy74KtOG1dJcF2dCcxtDG7p4fOVA1tas5ZP4uxCH4bU7I2KFGl520N78eM7rOd9vviKa5Xizy9076KTuTOZ5fn83o7SnL+g4mXwbwsfoKaMZutcoNs1dRk1dDYd+/1AuXzoCKH+lrRx9yrVXxQ5+61M1aVMqpHNe6og0xgCMJoN32pu2Ze+S2NyFapIBO+wMuJOXwYv5md/65gWqUOnk5SJyBfAFYAGQzpmgQGWMvmonl3F478ZHAJgcm85L8SdRTSISY2LsQFYmOnhq3fNbrRI+272Iz/SfwCNr2tm9cQCTh3Tz4cZ6rntrBE/2PMH0usM5ecdOEh9vz4aeei69/xhmr+tgRG0/OhNx2rUb6KBhfRN7D1YWd3cwfUh/zpz8Ao1NnbSuHM5zSyawZEMPg3UwY2OHsTRm1ca96LlXueygvSv298o38MsRg0psaHaTaplYGQxeEquJ8dmrP8um9Zuo71dPbf2W/85KLUHjxaSkEsWOXcNEIRgMoWD/sy+hpq6B9e8vYKdjz2LIhC3zuHLK3FT7/Mwpycqv9M0Edil1L7InRl/2Bxk16z4f2QZgZikH1WRvgfeWWmv1b2HyeSZjZZ3+vw2LaE0uYmEnTPz3m/z+oL1REkyMHciKRAffXPRSb1s/nvM6zYd8h85EHKC3DMPbnfBs93iogSfXQN0LU/n+SQ8w78UxbEzUcMroOFcuX8ualHNgSKKltz2g1/jLt+LnBrlWCb36fuQSFDsTO689UlEx+kSkg203FKdRVW2ucJdsUa36FEVEhL6D+rrSVi7vdCW1ye51RpucY2qKWhjtiw61fRo54Jxfutqm3/Mzu9cF53urrq70iUgNMBdYkd7Dl4MlQB1WjT7HVGSlL1gfUuXILP2QSbqY+1avZWxHST83vHZn1rGGrz37z21W415MPtNbnH0Kh7KeTcQ0Bglolj68qM/QvPpQzrrpeMb3jzG+qYvG2jgzm3dk1ronmMhkgIL7+cql0Gfu9fchsB4kFZIRyd6pqk1+98FgyIWd8e+XPgVXm6Jr9GFqihqqiCDrHwRYA0vAZS/SBcDbwDYOcxG5PnW7jVhZiJ8gw/BT1e/YuUHFwjurMblCtrFnZ59gZtbQG6d9rvd1JbHV9ZmrcTdO+xxH9LWytf61/d+938IXk88gPTFWrpvCAUMTNNTEWZP6iqytWWv9kthi+Dld4auGwptuo3hXALnSiMhsVT0kY8VvK4K60pdNNWpT1PE607CfeOdEjW72TlNTNDdG+6JJlPUPgrOQZO2DdsfsE5FRwGeAS4ELc5wyN/XzFeCBUu9j9vRVgExj7VcH7gbAT59fsM1r2Y+/9uw/cz6fTeZKYHr1L/3cSX0P578mvcnw0R8Rq0nwypoZHN7zaQDmJ5f1tvHrg61U505X/fKJS6FkBn57kYoJhtcZS6MysVLVQ1I/I7HiF5T/SAzukc8xlU+f/P787WiTnfNKJSralI8g1BQNIkb7oklU52d2zqsUDvb0DRWRuRmPZ6nqrIzHvwN+AOScT6nqnwFEpB/QlXZSpUJC+9jtRMXjzKLqecjFT+a8WdBYc5Mfz3mdH895na89+89eY/HHc15n539fxdsrRzFw0hKSPbXsM3g92zXEtjL4Mkkbf2Dt8Ss1yUvQP+di/UtPFgP0PgaIyCwROd7vjmQjIg0isoeI7C4iDX73pxwC9HkbXMJO8eEgYTc0K8jvIQyo6iZVfUhVv62qVWvwpTHfp+hixzi6fOlOgSnNZWfulTlH8+27m8reaecAVqvqARlHr8EnIscBq1T1FRt3fQJozHjcCDxut8u+bC6qZnH56fMLelf5KsWiW27kqv03cN1jRzC8fxvHj17FD0YN41N9dmRYzQSG1Uzo3eMHcNO0U3jpsC2ry5mGYD7ylXXIfi4oouIEt7+vquL4ANpU9VxVfdDVzpSBiNSKyJXAcuA24C/AMhG5UkTq/O2dwbAFu/oUFJyUeXBVn1ScHwaDIbDki2Lw0mAq1m56Lph9OKWU0l1uo2rvKMLBwGdF5H3gTmC6iPwlz7kNqrphy/11A2A7m5kJ76wi+t73W9Z97Sz2HL+Io055hxMeHsXsV/djl+1WUF+/ho2btufiF3dhdL8Y2zUkeKN1PK2JpUXbzbXM7tUAtJt9zI006ZlhBK6FEigkozNR+g1WKMJYVe0AEJFm4KrUcYGPfTMYCk54vLpXGqfalNmGXc+8mzWwVEG9q29vMBh8JFOfKjU/c0qmpgXJCZcPBdSFVC6qehFwEYCIHAZ8X1W/lOf0ThHZT1Xnpc7fH9hk917G6Ksy3rzxVt7EWs0DGJqo46/jNjL6hNdY95+x8OIuTB+xmsbaHl5aMwZiY3ii56mcbVVCRIphZ79O9mtOJ0huio9GK1nCccDOmWnQVbVdRL6BtW/GGH0GX8jWBT/0Kd89i+mJU21yc29LhLTJYKhagqB/+e5r16FlF69zMBTDhzp93wXuEpGPUo9HYBVrt0U0cscbHHP2s3dx9rN3sbpmDfPfH0diVQ33zj6E28+/mw2b6zl5/gKmDUv0GnxpIzGTzEFmZ3neDeFxEg+eL3zBqTi4HQJRYnhnENFcda9SG4yruh6WwV+y9ahSYeWl6FN2G07IXOlzgwhp01aISIeItOc4OkSk3e/+GQxukkv/ojw/y/wZdlT16QI1+lDVl4GJwDeAbwK72twLCPho9EXlAwo7Zz97FwtuvJWWn43h+dZG6oa0c/eH/RlWM4FR/Tp79/wBzD/8+9tcb3eiU85+vmxxcNKWm5M9t76zSRXHR0BZICJnZD8pIl/CWukLJUabokMldCJ74hKGsKTcODf4wmL0qWqTqjbnOJrCUlqmEhjtixaVmJ8BJc/P3MaP76+DRC5usguwG7AvcFqueVg+fA3vDErKVcOWUg1rrnmGo0ds5s6vP8Z3Lz+Hs1uaaK6Lc81H9urYerGxttxMdeWKWaHnLyqhzbBMlGxwPvBPEfkqVu0YBSZhZZM60c+OlUuY9hUYCpP+DL0q01COPhltMgQNo33RJQiJT7Lx4ntWjgY6xX17rjAi8nPgMCyj7yGsuqOzsZLpFcX3PX3G8AsWE++Oc9ce64nttR2LNvTw9b0XsaptIEeu35HTd36fX37quxz2n99VrD/p70dQviPlCqWVySkaEytVXQFMEZHpwO6AAA+r6hP+9sxg8J5yPduByzKsRDYbp4jMVtVDRKQD651K5k+z2rctZm5mKEa5DgIvDb5KoKl/FeZkYG/gVVX9iogMB260e3Eg9vTl+pACViOtarjouVc5/+1ulv1+OA/feDsTD3uJl1YN4/LT7+eXr44C4IqD9+eKg/cv2lY+D3t14zy0M8DhnQCo6pOqer2qXhc1g6+QNhl9Cj9Gm7agRHdPn6oekvrZlBnWacI7nWG0LzqY+ZkLqJXIxc7hIptUNQnEU9nSVwHj7F7s+0pfmkKhKuZLWVnOfvYu/g4c9K3Pc/A5D/Oj3/2d5PG/5Zev/p3XP9meF9+335YboVNRIywTJYOFmeBEEz/DOoNK1LVJRBqACVg27mJV7fK5S4Gm0PfbzM3Cj9HA8vFgv14x5orIQOBPWNtqNgAv2b04MEZfIYy4VJ4Vx36fgy59mO6DpiPxzSQ/fIAFq0Ywst8GhtbYdioYchD1iVU1YbSpuonUZx/hOn0iUgtcBnwV+BArrHOUiNwC/ERVe/zsn8EQVnJpoFNDMMyJsCpt86nqN1O//lFEHgGaVfUNu9eHwugz+IMM6Et86ET6PfNnksO247Dd+wOw+i3bdSA9oxKi4oUHS4lUcXZERIBRqrrM774YDEEgvBOe8IRrlsBvgCZgrKp2AKRCo65KHaamqMEQMvxeZbRC4itr9aXmXKcD41T1EhEZLSKTVdXWal8g9vTZwe8Pt9oY+dBV/O+PP0/dtbdz9y+OJ9E0hOE/7cOIC7v55MkYP3zOdlkQ17H7XUhPoIJi8AGWNz1C+2ZSdfru87sffmK0yeCUXNrkBK++c1HSpiyOA85JG3wAqtqOVevqWN96FXKM9hmyKWd+5sV9vEZtHi7ye+BA4LTU4w7gBrsXm5U+Q17a77yea7keZsLvEtZz394u9T1bV7l+ZIbQBWWgl07wE7OUwAsiMilVNNRgqCrS+hR2bUoncokoqjlc8qqaEJGKb8oxGKJI2DWwFJKV39M3RVX3E5FXAVR1nYjU2704VEaf2T/jP9ev+5sv942KmFhen8hNrD4NnCci7wOdbEmDvpevvaogRpuqm0rqk5f3irDRt0BEzlDVrWpZiciXgIU+9SkSmLp+hkrPzwIzH9TK7+kDekSkhtQCooi0ALZ3Y4cmvDONSRdscEIQvysRDKE6Bitl8HTgeKxQquN97ZEPGG2qPiL1eXsYei4iM0TkHRFZJCI/yvH6NBGZJyJxETk567UzReS91HFmxvP7i8j8VJvXpfa65ON84HwReVpEfisiV4nIM8B3sEI8DWUSqbFgsE01f+7pPX12Dhe5DrgXGCYil2IVZr/M7sWhWunLJDib3w1BJaix4hEM7/wYOAkYw9aacokvvfEZo02GYgRTm7xxMKW80jcARwLLgZdF5AFVXZBx2ofAWcD3s64dDPwcOABrjvVK6tp1wB+Ac4EXgIeAGcDDufqgqiuAKSIyHdgdKxrh4ajVFfWb7O+n0T9DPoKpgc6p9EKfqt4hIq8Ah2Pp2ExVfdvu9aE1+gwGQ2C4H2jDqhnT7XNfDAZDsJgMLFLVJQAicidwAtBr9Knq+6nXssOUjgb+raprU6//G5ghIk9jpSp/PvX8bcBM8hh9Gfd5Eniy/LdkMBgM6seePlR1oYh8COwGrHZyrTH6DIYMvPcihSJc0ymjVHWG350wGAxlkixJm4aKyNyMx7NUdVbG45FAZkmX5cAUm23nunZk6lie43mDwWCoCFZ4Z2XuJSKfxQrtXAv8FCt64hNgjIj8UFX/bKed0O3py0XQlnsN/nL50p3yhpUU+654HY6iaoV3Oj0CzhwR2dPvTgQRs88v2pSiF0EOeStxT99qVT0g45iV1WwuAbM7Vcp3bTltGgwGlwiynlWCCu7p+yVwFPB14B/A4ao6FdiLrLD4QkTC6ANj+BksMgUon/Hn93clKolcUkkU3gAOAealEjW8kfG8IYXf3zlDMEjrURC1Sb1L5LIc2CHj8SjgI5vdynft8tTvttoUkckicqGITE/tMUw/P1BE+tvsi8EhxullyCatfUHUwICTVNV3U6WxlqbD5VV1FRC320ikwjtN2vRokf48Mz/T7PTQxQQi3+u5av9VqtZWMjolG47zuwMGgx9k6ke2PuXSq+zzirWdPr/S6fA9ClV6GdhJRMYCK4BTgS/avPZR4DIRGZR6fBRwkaquFZEOEZkKvAicAVxfoJ0vA/VYK4Rni8hC4NrU4xtTfTJ4hJmbRY+oaqAjKluyIZbSwRiQTP2enkzaXsCLlNEHRlyiQnqgZw/47M+2nM+6kDh5RZQKIKvqB373IUwYbYoWuSYzTvTJicOqUt5vL7RJVeMi8i0sA64GuFlV3xKRS4C5qvqAiEzCSkM+CDheRC5W1d1Txt0vsQxHgEvSSV2wSi3cCjRiJXAplMTlh8B5WBk+a7D2/50H3EOEIp6CjNG/6JBvfpb9e67Hudqx83pQVwCTlYsqH4CVLC8t0vMyXrPdicgZfWDEJerY8fzYFYjKC0ko9ugZPCLQXkuDKxT7/yeokxcvk0yp6kNYZRUyn/tZxu8vs3W4ZuZ5NwM353h+LrCHzftvBK4GrhaRPsBOQAuwPbCfvXdhKBczN4s+bs7Pgk4lE7mo6hg32omk0WcIN4UEwY5YBF1QgrzSJyLjgJ8AA1T15GLnGwzVRqGQ8bATZG0qBxERTWVTUNVu4M2M1+7MPsdgMDgnCvMzp4RNMiJr9JkCydEiKkKhQHYhKrcQkZux9tmtUtU9Mp6fgbV/pQa4UVUvz9s/a3Pw2SJyt437TcZK4vIa8IyqJlLPDwTiqrqhnPcTZcyKXziJtsEXXaMPeEpE7gHuV9UP00+KSD1wiIicCTyFFS5q8BCjfeEml9ZFQf9KJWQ2X3XEslfzFzJsRP6z8i5DHlgTlq3q5aUy1d0AHINVyPM0EdlNRPYUkX9lHcMcvpsvA7sA+wK3icj/E5FmtiRHMBTBZLcLD9XwOWlSHB8hYQaQAP4mIh+JyAIRWQq8B5wGXKOqt/rZwWrDaF/4MJ/XtiRRW0dQiOxKXzZm5S/4uCEoYRClEvf0FSuAjKo+KyJjsq6bDCxKp/dNhTKdoKq/pvzsmyY5gksY73ewcUtXgq5PUV3pU9Uu4PfA70WkDhgKbFLV9f72zGC0LxwEXbv8oJJ7+jIRkb2BT6Ue/kdVX7d7bVVOzMyXN1hUm8dPEccHxQsg52MksCzj8fLUczkRkSEi8kdgXxG5qOD7UN2oqler6lHAscC3gC9hpVBfVuhaQ26qbSyEger5PJxHIITRSFTVHlVdaQy+YFE94yx8mM8mDzYLs7u5709ELgDuAIaljr+IyLftXl81K33ZGO9SMKg2MdHKZ+/MdbO8CqSqa7BW64o3bJIjeIbJcuc/1aZNRHtPnyEEmIisYFF1GlgCPsxuzgamqGongIhcATxP4TqlvVSt0Zcm+0vtldCYSdzWlCsm+f6WlSqwXg7JyorEcmCHjMejgI9catskR/CQXN9jLzTEaNO2lKMhhf6WQdanKNUQNYQfMzfzD6/mZ1HEB5+2YO1PTpMgt3M/J1Vv9FUCIyqVI6gTqkzU/vjMZICIzAIeVNUHHVz3MrCTiIwFVgCnAl8spQM5mAF8FSs5wlhgPVah5BjwGFZyhNdcupfBA4w2bYuXGhJ0fTJGn6GaMPpnKAcvs7EX4BbgRRG5N/V4JjnqmObDGH1ZlBr2WUg8jKhYuDXhCfPfUyk5kUubqp5b6AQR+RtwGFbSl+XAz1X1JhH5FvAoVrKVm1X1rVI6kI1JjlB5SpmkGG2yR9ANskoQVaNPRDrIHdYugKpqc4W7ZHCImZt5i5mflUalV/pU9WoReRqrXJYAX1HVV+1eb4y+PNgRmOxBYvYJ5scNQYnE31W9iwFX1dPyPP8Q8JA3d+29Rw+w0st7GCyMNrmL0aY04UzMYgdVbfK7DwZ3MPrnPkYDS6fS0Z0icruqfhmYl+O5ohijrwilDAYjMFtw03tuQjEMhi2Uo01g9Mlok8EQXoz+lY/RwDJRX/b07Z75IFWLeX+7Fxujz0OqWWD8CpUKcqKENMnS9vQZDK5Szc4po085UMJUbN0RIjJbVQ/JCPOUzJ8mvLP6MPpnCBOpElo/BhpFpD39NLAZsFvCqzrr9PmBGWSV4/KlOwVWyNMZ8kqohTVARGaJyPE+vwVDxKg2bfLq/dptN4LaFHhU9ZDUzyZVbc7+6Xf/DP5RbfrnJXb+lkGenznFytFg7yiEiDSIyEsi8rqIvCUiF29zL9Vfp8LUf5PSrrR+DVHVgjWVMzErfQZX8VpA7YQQBFvES67TVzSRi1+IiACnA+NU9RIRGQ1sp6ov+dw1g6GXSuhC2PUp6hU1RaQBmIA1X1ucSkZlMFQNQdGfQEc9OEDzlz12QjcwXVU3pJLizRaRh1X1hW3u58DAy4VZ6TO4wo/GvlexAVzsPkH3ImkJR8D5PXAgkE4k0wHc4F93DE6Jwn++hYj6+3OLpIrjIwyISK2IXIlVt/Q24C/AMhG5MjXJMhgiTxB0MD1XDEJf3EDV3lG4DVVV3ZB6WJc6PJn6mZU+Q9n4MXiLedTTrwVNWMoo2RBkpqjqfiLyKoCqrksVaDcYfMVokwM0uiUbgN8ATcBYVe0AEJFm4KrUcYGPfTP4TJSTkFRaZ4rtlcx8PnAa6Bh1kshlqIjMzXg8S1V79+KlErK8ghWJcIOqvuheP7dgjD5DyYR/wPqDD8U8vaYnJVgKICItRPJtGsKE0SdnaIRLNgDHATtrxgxNVdtF5BvAQozRZ4gYRv+8J72nzyarVfWAvG2pJoB9RGQgcK+I7KGqb2afJyLjgeWq2i0ihwF7AbfZrY9swjsNhgoTwUQu1wH3AsNE5FJgNnCZv10yVDNhmPAEcWUhqolcSEVQ5XgyQSgi6A0G+4RB/yCYGugUtfnPdnuW8fY0MCPPKfcACRGZANwEjAX+ard9Y/RVkLAMxGIEJR7bbqaoIKFa8r6ZNlU9V1Uf9Ps9ZJJK4vIs8APg11gF2meq6l2+dszgmCCM6XIJkzYFkQgbfQtE5IzsJ0XkS1grfQZDJAiK9gSlH55icz9fsQhQEWlJrfAhIo3AEeTXpaSqxoETgd+p6veAEXa7bMI7K0zYY8fLGci53nf236MahCJKbmVVVRG5T1X3x0yeDCEmW5+y96ZEXpsiXKcPOB/4p4h8FWvfjAKTgEasyZOhygn73Ay8mZ9lvhZ5DSwBl4qzjwD+nNomEwP+oar/ynNuj4icBpwJpCO/bCejMkafwTZuDvh8bVWDuEQwkcsLIjJJVV/2uyOG6sUtzchsJ3MiVA3OqRCt3DlCVVcAU0RkOrA7VlHjh1X1CX97ZjC4g1fzs1waGFX9c4rDPX3521F9A9jX5ulfAc4DLlXVpSIyFisbsS1MeKcPhHHAeNHnQkU63fS4BakYqGJlOHF6BJxPYxl+i0XkDRGZLyJv+N0pg3PCqk1e6VM5hK1QcTqRS0TDOwFQ1SdV9XpVvc4YfIZswqh/4F2/y52fhXELjlNU1dbhIuOA76rq31L3X6qql9u92Kz0GUrGiec7e2CHfaCXTvgmSjY4xu8OGAyZ5AvVtHO+H9oUlELFEdQmg8FQApXUwaDoXyn4sF3nVOBaEbkHuEVV33ZysVnpM5REJYy4qBqGJa70BTl755l5DoOh4uTSDbe1JLLaFNHi7AZDNZGOILAbSRAEp3woNVUrv9Knql/CCgVdDNwiIs+LyLki0mTnerPSZ7BFticm16bnfCt/pQ5mtz0/QYhHV0r2prep6rkud8ctOjN+b8CqieXI+2QwlIodbUqfl349+zmneBG25Lu3O8LF2UXkwkKvq+rVleqLwVBpCkU+eGlshU4DS8DdyE2799T21EpfI/BdrGRU/y0i16nq9YWuNUafwRPcEJJCRloovUIp3Nj4GyRU9beZj0XkKuABn7pjqEKcTBaMNuWmDIdUUURkBnAtUAPcmL0HRUT6ALcB+wNrgC+o6vsicjrw3xmn7gXsp6qvicjTWFnvNqVeO0pVV+XpQtoLvgtW1s60Ph2PVXLGYAgtTo0lo4HuoCjJClt9qWivrwLjgduByaq6SkT6YjnbjdFnCDdRE5CI2Xy56Iu12dhgiDSR0yYPskal0pDfABwJLAdeFpEHVHVBxmlnA+tUdYKInApcgWX43QHckWpnT+B+VX0t47rTVXVusT6o6sWpNh7DMho7Uo9/AZiaogZDiURNA53iw3zuFOAaVd3KWaWqG1MlaQpijD4fiMogyfbwhOV9hTGEIMiIyHy2aF8N0AJc4l+PDKUSljFcjLB6nyOqTZOBRaq6BEBE7gROADKNvhOAX6R+vxv4HxER3XozzGnA38rsy2hgc8bjzcCYMts0RIQwaIQdwqp/EFkNdA1VPaPAa0UzEhujL2QEuYBopeLEw4xqJOv0HZfxexz4RFXjfnXG4A9B1ibYttCwIZuSMwsPFZHM1bZZqjor4/FIYFnG4+XAlKw2es9R1biItAFDgNUZ53wByzjM5BYRSQD3AL/S4hkTbgdeEpF7sRxVJwJ/LnKNwVCUsOgfGA10E5fLMRRFRKZihXDuCtRjOdo7VbXZzvUme2eFKXew+TlYndzbq9pZUSCCdfomA2tV9QOswqH/EJH9fO6TocL4PZFwUjvKaFMOtOTsnatV9YCMY1ZWy7ksyeyZUsFzRGQKsFFV38x4/XRV3RP4VOr4ctG3qHoplkatA9YDX1HVXxe7zhB9wjw3c3p/o4HuoWrvcJH/wYp6eA8rkcvXKLKPLxOz0ldByqlrFxScZsAMqnfdzxCCCGbI+3+qepeIHAIcDVwF/IFtvfmGABMVfTLaVBoeJnJZDuyQ8XgU8FGec5aLSC0wAFib8fqpZIV2quqK1M8OEfkrlvPptmKdUdV5wDyH78FQZYRx+0pUQiPD8j6syK3K7+pT1UUiUqOqCaxohzl2rzUrfT5RTEDSnpgwfPEN9lEiWacvkfr5GeAPqno/VtiBISQ4rdNktCmaqIrjwwYvAzuJyFgRqccy4LKz+z7AltqeJwNPpkM1RSSGlbzgzvTJIlIrIkNTv9dhhZi/SRHE4ksi8rPU49EiMtnOmzBUD7m0LSxzMrPiV1l8WOnbmNLR10TkShH5HtDP7sVmpa9ClOMlKidWPOjoaoqjAAAgAElEQVRx5oUI82bkQpQoAEGu07dCRP4XOAK4IpV+3TiUQoJf2uTG9W5QSh+iq03ur/Sl9uh9C3gUa//Jzar6lohcAsxV1QeAm4DbRWQR1grfqRlNTAOWpxPBpOgDPJoy+GqAx4E/2ejO77H8aNOxkk11YO0HnFTOezSEm1xOr0LGUFCjBAyVRyufv/PLWPOrbwHfw4qQOMnuxcboCwlOJybZxYrTuCVSXi+/hz1EqxDJnNtXQs3ngRnAVaq6XkRGsHVtLUOEKWUMeqlPQQkNKkeb/HkP4lmSKVV9CHgo67mfZfzehbWal+vap4GpWc91YtX0c8oUVd1PRF5NtbMu5TU3VCnV7vQqlag6vpxS6ehOVf1ARFpSv1/s9Hpj9FWI7MFdyn/o+SYRTtryK07dibg5Nfictp/Gj4mVEsni7BtFZDFwtIgcDfxHVR/zu18Ge+QaO6XqU64xWKo+VXIC4YY+ZbZRrjZB5fXJgzCkINKTqhuYDh1tIRS5sgxeUa7+FdOsYm3lez2o+pc+v1g75c41neaP8ANrPlcZ0RQRAX6OtcInQExE4sD1qmq7RJYx+kKEFyJQqZUyt9ovFnLhhnB5TdQmViJyAXAO8M/UU38RkVmqajujlCE4lDouyjX48l1bicmP1/rkVgip12gyclEI2VwH3AsMF5FLsfYP/tTfLhmCRND0L1/bbuJm+2HXQKdUcD73XeBgYJKqLgUQkXHAH0Tke6p6jZ1GjNEXErwe9E4nWGEYjPnwt+8SxfDOs7HCpjoBROQK4HkcpBE2hBsv9SlI2uS1dvitqxHMLLwVqnqHiLwCHJ56aqaqvu1nnwzhpxLzs6CETkZdA52hldzTdwZwpKr21i5V1SUi8iXgMcCW0WeSLRgAS7SCICpeT9h8FxSbmZ48zPzkBcKWDJ6kfo/27DFilDMuKuGFjoI2lRrmVSmsUKWS6vSFhlSSqf2wSkIMAU5JZ/I0VC9B1r9K3cMOUddAx1R2DleXafD1dkG1Faiz24hZ6Qs4bocM+EUxT5Vb7ynX3smgJHaALSUbIsYtwIsicm/q8UysbHyGkFDqnrNcBGWsOcEvfXKz7bIJh4OpXO4H2oBXgG6f+2IICKUaVZXQjEpgZyUxzBEUXlHJPX3A5hJf2wpj9PlAuV/wIBgxTicufm5K9vtvlU2JiVwGiMgs4EFVfdDdHpWHql4tIk8Dh2Ct8H1FVV/1t1cGJ7jpHPFbn0pxlPnhSQ+aLkH0wzuBUao6w+9OGMJPlJxexYjie3KLCjrK9haR9hzPC9BgtxFj9FUQNzPi+T2xysaP/gTZsCtEiRoR5Dp9qOo8YJ7f/TCUjtMxFKR9JsWotD6FU5vCF65ZAnNEZE9Vne93RwzBIUpOr3QfMgmS06tQhs8wohWy+lS1xo12jNEXAtycXHk5uCs5kMMy2cwmvW8mSohIA/BNrJU+BWYDf0jV3jJEGLcNP6/1yWhTfpTohneKyHyst1gLfEVElmCFdwqgqrqXn/0z+I+bTq8gGH6ZGKe8d4RNMo3RF2KcDOQgT0SiLAi5iODE6jaggy3ZOk8DbidPsWVDtMg1+QmTNvm5lyVoRDi88zi/O2CIHmFzeoFxylc7xuirAqphAOYrihxEIpjIZRdV3Tvj8VMi8rpvvTEYXMKtcK/stgKpTRq9KIQ0qvoBgIj8GbhAVdenHg8Cfgt81cfuGUKMG4ZfNczRoogVHREuL74x+gyRIujiGdEQqldFZKqqvgAgIlOA53zuk8EQKIKuTQAaQY9UFnulDT4AVV0nIvv62SFDdRN0XXDT8eV2u76jJSfm8w1Tpy8kBKLGnAeU8p6CLpJVyBSsBAnvi8j7WIXZDxWR+SLyhr9dM1QCo02GkBBLre4BICKDMc5vQ5lEVSvM/Kw4qmrrCApG7EJGmDLmuUFmaFRU3ncEnekmBbohUmPUDlF7r1YUQjTDOzP4LZaD6u7U41OAS33sjyEiVJv+QTTnZ04JkD1nC2P0hZxyszKlr83VTimD2E2PV+b9IyMoIQwHsEE/VV2Q+YSIHKaqT/vUH0MAcCNxQL5U6kHSpugQ/ZINqnqbiMwFpmNl7vxctnYZDG7gZsZMN/QvVzvFKPQeIjk/c4im/oUJY/RVKZkDOV8dQLeoVkHIhRK+FL82+IeI3A5ciVUk9ErgAOBAX3tlCCWVrClltCkDDZ/XukRWAi9hadVQEZmmqs/63CdDBHHLKe8Xft8/DITNiW+MvirHzUmP8aTbI2wiYYMpwBXAHKAJuAM42NceGQJBuZoQZR0IGlGsIZqNiHwNuAAYBbwGTMXagzzdz34ZDJlkO+W9NL6MxpZHkPbr2cEkcgkZ2QO0XDHIdX0UPOlB9lCpOj8CTg+wCWjE8p4vVa2CPICGrXBbm9xqIx9Gm7YlgtqUzQXAJOADVf00sC/Q6m+XDFElSE4v45T3gBDqpDH6qhAvJh2ltOn1hC6IwqRYiVycHsAAEZklIsdXvNPFeRnL6JsEHAKclpEowWDwnSAZWkHVJrASuTg9QkaXqnYBiEgfVV0I7OJznwwRoFJOL+OUDw7Wdh17/wohIjuIyFMi8raIvCUiF3jVZxPeGSK8EBU3ki0EkSC/nxLDO9tU9VyXu+IWZ6vq3NTvHwMniMiX/eyQwV/cGn9uhDYFTQuC1p9MIhh6ns1yERkI3Af8W0TWAR/53CeDoZdqccp7fY9K4ZJmxoH/UtV5ItIEvCIi//YiyZRZ6QsJpRp8Tq+rVFY8rwd7UD3psCWZi5MjiIjIDwBUda6InJL18q4+dMngE6WOt0L65EaCqVJ1JophpcUoJbQzSGFLdlDVE1V1var+Avh/wE3ATH97ZQg7XszPckUERGV+Vql7VAI36vSp6kpVnZf6vQN4GxjpRX+N0RcCylnh83pgFUrnW0igvJ74BFVQrGQJzo+AcmrG7xdlvWZq91UB2eP8R2PfK1mfKjlmi+mPl/oUVG0CK5GL0yOsqOozqvqAqm72uy+GcBJ0/TPzM+9x4BwbKiJzM46ckVsiMgZrr/GLXvTXhHeGADfCm/J5zb3whNsVDC+yUgVeTELoHS+A5Pk912NDBPFqvOUK/7GrK25okxcEXpuIlDblREQEOB0Yp6qXiMhoYDtVfcnnrhlCjteJUtyKcvBzfgbh0EG7KI6yd65W1QMKnSAi/YF7gO+qanuZ3cuJMfqqBDvFPd2YVJWCW0Wcw0KE0lpqnt9zPTYYcpJPm9w2+Eqh2rQp6kYf8HssCZ4OXAJ0YE2yJvnZKUM4cStZSz59C2J4ejkEtV/l4NZ8TkTqsLToDlX9p0vNboMJ76xSgjL4MsWuVC98UN6LHSIW3rm3iLSLSAewV+r39OM9/e6cIZwEZTy7oU1hIl2nL+LhnVNU9XygC0BV1wH1/nbJUO3k0rxSnF752ioHN7INOw17DQ029/MVWw1MRSDcBLytqld72WVj9BkcU4lY8EgKRIqoJHJR1RpVbVbVJlWtTf2eflznd/8M1Yef+1SiQFS0qQA9IlJDqusi0kKkgi8MBveoNqe8jxwMfBmYLiKvpY5jvbiRMfqqCLc810HNBhUWIrTSZzC4QqY3OfunE0rVJqNPQAm6ZFebRGSGiLwjIotE5Ec5Xu8jIn9Pvf5iKpkBIjJGRDZlTIT+mHHN/iIyP3XNdSlveTGuA+4FhovIpcBzwGX23oXBEHyMU76yuJHlWFVnq6qo6l6quk/qeMiL/hqjr4opdeCWuqHY6f3snB9G8Yl6WnSDoVxKNcSCVtIhTPqkSElHMVIrazcAxwC7AaeJyG5Zp50NrFPVCcA1wBUZry3OmAidl/H8H4BzgZ1SR9GMwap6B/AD4NdYIZ4nqOpdRd+EweAhbpVnAOOUryTpRC7lhndWEmP0VRleF8Us1m61e9UVK5bI6WEwVANeZYwDo0128GilbzKwSFWXpMoj3AmckHXOCcCfU7/fDRxeaOVOREYAzar6vFozqtsoUG9PRKaKyNMi8k+gETgPOB94RkRMeRlDoChHi4LklK8GwhYGb4y+KsHNulRmsBsMBreIqjYFqS8eU6z+1EhgWcbj5WxbeLj3HFWNA23AkNRrY0XkVRF5RkQ+lXH+8iJtZvI/WGGcfwOeBM5W1e2AaVirfgaD7/hdH884vpyTVLV1BAVTssHgqtCU25aXnv6gECQBKIdUls5cb0YAVdXmCnfJYMiLG9oE0TbmSlSmYvWncq3YZd8q3zkrgdGqukZE9gfuE5HdbbaZSa2qPgYgIpeo6osAqrrQ3lZAg8Ef3AxZr7QGRlkrIR3e6XcvnGGMvirAjVjxYoM3V5tBrWHjNyHTiLyoapPffTBEDy/25XmhT1HTJqtkgydNLwd2yHg8CvgozznLRaQWGACsTYVudgOo6isishjYOXX+qCJtZpIZJb8p67WoSLIhAniVZM+LVcSoaaBj1FFx9kBgwjsNtigmGFUxwF1APcyQZzBUI3YmM070qZq1zKOSDS8DO4nIWBGpB04FHsg65wHgzNTvJwNPqqqKSEsqEQwiMg4rYcsSVV0JdKT26glwBnB/gT6YmqKGQOKGMVYsLDOf06uatc4twpaMzxh9Bk8oR0zsimA4625pSf+CiIjMTv3syJxEpR/73T+DIR9u6VO+dsKpTd44pFJ79L4FPAq8DfxDVd8SkUtE5LOp024ChojIIuBCIF3WYRrwhoi8jpXg5TxVXZt67RvAjcAiYDHwcIE+mJqihqrFGHfekURtHUHBhHcaXMGPSU4YhayMEKoBIjILeFBVH3S1UyWiqoekfh0GTMB6e4tUtcu/XhkM21JpfQqjNoF3cY6pmlMPZT33s4zfu4BTclx3D3BPnjbnAnu421ODIfq4qU/5QjyrIT9DGPf0mZU+QyAJq6fcDiWWbGhT1XODYvABiEitiFwJfAjcCtwOLBORK0XEeM8NgSTqE5FSMeVkDIZoYpzyXmGvRl+Q9v0Zo68KMLHbwcKuSARRMLL4DTAYGKeq+6vqvsB4YCBwla89M4QGr7Upyg4kt/FoT5/BYIgQRlO3YPb0GQKLV5OrbAEwBmZ+IuZNPw44R1U70k+oajvWXptjfeuVIXQYbQoGEdImgyEUGE0KN2Hb02eMPkNZ+DmpCqu3KUIrfao5OqeqCcwigKGCuJ3Bs1zCqE3p/SlOD4PBUB5eapNxfHlHGHXSGH1VhBuFOXMJSOZhl2LnRlmYIuRNXyAiZ2Q/KSJfAhb60B9DSHHDSCqUltypnkRZfwoRIW0yGKqGfPppnPLeEzbHvTH6qgQ3B6BbbRWajEVVMKzsner4CCjnA+eLyNMi8lsRuUpEngG+gxXiaTAUJYjaBNVp+Jk9fQZD5XHL6WWc8oZiGKOvCoiqARVWolKnT1VXqOoU4BLgfawsnpeo6mRVXeFr5wwGF6imyU3E9hsbDFWJccpXFhPeaQgUXgxOtydC1SYsUZtYqeqTqnq9ql6nqk/43R9DeHB7jHthpOVq006/w6hfUdMmgyHoBF0nqsnxVQphi9YyRp8h0ARdEJ2iNjM9BTXzk8FgsMgVThV2THinwVA5wuCUz0dUHV9OCKNOGqPPYBsv6/2V2m7oREUjtafPYCiZMKzylUuY9MmEdxoMBidE0fHlFJPIxRBJKjGhKuUeQZzoFSMqe/oMhiDgpTMq8x6VvM4fSlEmo00GQyl44fQKmlM+8piSDYYgkPa+XL50p5JTludqsxI47afxNBkM4SJzzEZZmyBcK30GgyGcBNUpH33CV2u51u8OGNznR2Pf6zX40riVErhQm36KQnbfgooVQhUcATAYKk32OM2nTWkds0OmEZmvXb/0ISzaBCZc02DwmkxdSutCWJxDTjQ5TZj0zylWCS6/e+EMY/RFFDsGXymC47U4OW0/l4EbdIzRZzAUJj2enU4yvNSnUtoOky5BsBIOGAxRxK4mRMXpBVE3/MKlmia80+Aa5U64Srk+fGJi9s0YDGGjVG0LiwcfTCIXg8EPcmlEptPLaVtebnkptd1wzdGcYfb0GQJFoUHqhTD4MckJ38TKlGwwGCD/2I2KNoUNFXV8GAyG4OOHUx6ibfCByd5psInXHpmoYFcwQvN3FEhK0vFhMFQSo03FsZuEJkx/R7PSZzBUjko6vbxsN2j3rBSq1p4+O0dQMEafxxSbPAXBCxKEPgT9b+QmZqXPEBSCPO7CMlnw++/kFia802CwME4ve0TR8eWUsG3RMUafh+TaXBtUKlHrqhSiJhbW8Hdu9hkMbpJZ0iXXa0Egc1+L0abKYPYbG6qFoDvkgzR/jKLWuYXZ02coSpAHULkC43eGqDBgVvoMQSXIY6hc489oU3HMSp+hGihkUAVlrGb2LQiGXy6C8rfyE7OnzwDkF5VKD5JSJkrlCExQRCDI4RlmT5/BTyo5Lordq5KTmaAkhwnqfm7FrPQZqpugjclM3Ih4MI4vd0nX6QvTnj5Tp89lchXezCSz9oobA7BYLZZCJQ0K3b+UIpxeCUo5dbqCVtLBCu40RpzBH4rpkxf3KsWrXkybCl1b7Dq3cVpLK8gYZTJEnUo65IvNP0rRslLmZnb7U2qb5VwTpPlZNWBW+iqM13HauQZgNafaDeJ7SJbwz2Bwk2Ljwotxk8shUw5BHNtOCGL/VZwfBkNY8CPKwc49s7XAK23wW3OCGOFQHvZCO+2Ed4rIzSKySkTe9LLHxujziGJebTe++IXS/boxuJx6nypxHyf4LXC5CXYiFxGZKSJ/EpH7ReSoit3YUDHyjQs39clrbSp0j1wETZ+CqE2mhqghyvgZal6sNEMlohaCqDlRQG0eNrgVmOFBF7fChHe6SJD2ylQzQRY3Bc/26InIzcBxwCpV3SPj+RnAtUANcKOqXp63f6r3AfeJyCDgKuAxTzprqDiFNCPKehI0PQhafzIxMQWGqGPH6eUFla7Jl0kQNSeIfXKM4lqSFlV9VkTGuNJYAcxKn0dke3CK7b1zSiQGjAcE/++iXoZ33kqWp0hEaoAbgGOA3YDTRGQ3EdlTRP6VdQzLuPSnqesMESN7pS3XpCPo46icfb5e3ccOwf67lpLGxd6ER0RmiMg7IrJIRH6U4/U+IvL31Osvpic/InKkiLwiIvNTP6dnXPN0qs3XUsew7HYNBju4Pc6j7EQrl6CW4CkFh4lchorI3IzjXD/6bFb6PKTUpXs7bZZLUAadm3+boLynYiiJUi4bKiJzMx7PUtVZW7Wb21M0GVikqksARORO4ARV/TXWquBWiIgAlwMPq+q8UjpqCAelJFEp1l45SQbKvb/b7VWbNqWLs7tNhuPpSGA58LKIPKCqCzJOOxtYp6oTRORU4ArgC8Bq4HhV/UhE9gAeBUZmXHe6qmbqosFQkOwkIl44vdzQwUKYqLLg4GClb7WqHuBlX+xgjD4fKEVQyk3GEsSBG8Q+eU0Z2TtLFYyRwLKMx8uBKQXO/zZwBDBARCao6h9LuKfh/7d3PyF6XWUcx39PEm03JS3VjU0pBePCiCAt7VYQpbsgVppGtAs3otG1WVRLySIroaVuoqnQQltKRZhiIVDqSiompW7iokSjOOnGGgxFaGhmHhf3neTtO+877/1zzrnn3PP9wNCZ6fveufdmzjP3ee5zzynU0NjU5zm7GuNAriI9o7e08CRpPuk7Kump2eevSXrOzMzd3517zUVJt5vZbe5+PcaOoi4hE75QRa+c1Fb46qO0p5pJ+gIKMeBjDbL5C6wcBl9dbVOjWja/3so45e7PSno23u5gLEPjU4rYFEPXmEds6mxdF0KbwtPN17j7DTO7JuluNXf6dnxL0rsLCd9vzGxL0m8lnfKcVkFGccYoenV9Xyo57lOOSgs5JH0BDZ0dKkXrZuylItpsv/ZgkngJhk1J9859fUjS+yl3AHnIZaH0seITsWm9nkswrOtCaFN42vM1ZnZETcvn/IzC33H3K2Z2h5qk77uSXmi3y0A3Me96hV6/eSgKX+2FyvnM7GVJX1VTRNuU9HN3Pxtm67eQ9A0wP0jHSvjaDqZUSdmQhU+7nKNyg4j3fabvoJmdkfS6u7/e4X3nJR02s/slXZF0TNLxPjuAcoRcziBVbGpryMLEErFplZ0lGyJoU3jaec2mmR2QdFDSVUkys0OSfifpe+7+t5v7635l9t8PzewlNW2kJH3oZepFr7bbr73w1YW7tB1u9s7Hg2xoDZK+HsYcFKVeUIQwhWNvLqx63em75u57zva0qlJkZifUTICwX9Lz7n6xzw6gDGPFp1LGZ4wW91KOfZ22s3F21KbwtCHpCUlvS3pU0lvu7mZ2p6TfSzrp7n/cefEsMbzT3T8ws0+pmZTqzRg7DwyJqTklWqtiX9ciYao1UktQWHcnSV9Xy37Z9xoAoQZ8bosLt/mZXZ7VmXpgmBdrsfVVlSJ3f0PSG1F+KLISepHfttsLMX5TxoApTroQQozINHtGb1fhycyelnTB3TcknZX0opldUnOH79js7SckfV7Sk2b25Ox735D0P0nnZgnffjUJ368i7D4KVFIXVlch5o0Yurj7usSxnuu59svW5IKkL6JUz/CFNMbyEvUECKlJ+Xq1dwJ7Cl19zXH8Ep/iaWYWjnMBs6zw5O4/m/v8I0nfXvK+U5JOrdjsAyH3EeWjC2s8tR4/d/rQSqrb47lXX3Lcp5hc8e70AaGkSohymrxgUY77FFth1y/ATV26sHIqeo3xvF3ojpBaNYuzlxU1SfpaGKN6FGqQ5dzCVGcgcW170olcMHGpx3fIcdslPuUax6Zk28q6gAFiyDHhG0Ou+5WTwnI+kr4Y+g6UWAMsx4Gb4z6l0vNO39qJXIB1hlSTiU/TFnH2TiCKkgvyO9vKtehVYwysAUnfGiFvgw+d8CW2VEElh2MdT+8lG4BdUrT75DBeiU1pkPJhyih6tZPjPuWKiVwmJNQt/twvqKQ0F1W5HOuYmh5wnunDMKGeCSY23ZLL8Y6JO30oBQX5OHI43mI47Z2TEbuCntvAyvnZv2lxJnLBILHHKbGpTrR3okbElltyi/25YyKXCUhRPcp1YMW8uMr1mJNzyZnIBT3Fjk+5jtPYiV+ux50a5SiUoKYurFRqO95QCsv5SPpSY2DVrbnPx0QuiI9Yg27KW2gY9amtCwt588KyPpK+SEqqosdW63Gv4jzThx642AmvxmNehfZO5CzF2nK5xgM6HfJVWsTcN/YOjOn05cO0M86JcS5KPA9xNbN3dv1AXULEpi7JXe7jlNiUgDXr9HX9AEqyKrbWGg9qPe4QXK5tb/eRiyqTvpDJXukPAe+ci9KPoyTu250/UI/5sRgr8SsFsSmd5k5f9w8gtNR3tXJPfLg5kanZ7J1tPnIx2fbO05cP3/yFXjdgQl9YLbsVn+Pgol0AGMdY8YnYlOfxAthd8BoyVoe+PwcUvPLHM30ZGWvAlDJQQ/eJlx5g0+i9ZAOzd07MGHGi1ti0s03sjWf6EBsF+fVqW5anZIXlfNNL+mJfKKR4mLgUUz62WFy9l2xg9s7CpUhiphSfhiZ+OR9brpi9EylQkN8bRa8yuMqLmZNJ+mKsXzV0kEx1kE31uNJwntGrUOj4NPWET+p3gZb7MeWs6UEo6wIGGKKWeFHLcY5hu7CQOYmkL8cK+lQRPIbr2d6JAqWYdZLYhFBI+hALXVjd0OlQBp7pSyxVNbikARTqTuWiks5Btpx1+mrR9w9213FW0riMFZtibbM2JH0ILccurBKk+vuBATKbmbONYpM+qtvrTWH2qunpPZELCkF1dr02sSn2gsT4JBZnR2h0YaVVy9+PXPBMX2RtB/fUfvG7JG8EwLwNmMgFmWsz9pa1JJUcr7pU3hfPz7pjJ5alt21j7wGmgi6s3eh0wNiKSfpqTfikWxXvPhdITEqTm50lkDtjyYZMDYlNpY+rndg0JD4Rm/LAnT701ScGll7wkvoXvWpqU506JnKJoE8FfUpCtEItBpkuVbgpn9sx9HymjyUbMlRzMUoKc/FCbMoFs3eimyFjdQpjN0TRa2c7Q/cDY3AmcgmNhG/38c9/r+ux0zI1Np7pq0ltsWn++8SmsrikLZI+BFZDDBx6jItJI4WvchSW8+Wb9K37pa/hF31ZO8C8VOdgCm0YeSHpKxmxabfFc0JsKhN3+tAWBXmKXrVr5mgoK2ZmmfSt+uWfcgDpgvNQOJZsKNZef5gZlw3OQ7lI+rBO7ckeMK+0q7nskr4x72iVYlV1u297wDqc/5Bo75wSxsZuy+JTrAXlOf/huFxbRmzCanQ5LDfmXTq6HUbk3OkLquZf5FzuKNT8bxAPF1ZTUPPYyCU+IZyYz/SZ2SOSnpG0X9Kv3f30wv+/TdILkh6Q9B9Jj7n7P2b/76Sk70vakvQTdz/XZptIhzF/C+eiLoXlfLIuWaqZ/VvSP+PtDlCc+9z9s21fbLbPrUetxfXxO+7+YOc3VoLYBOzSKTYd2HfI7/j0jzv/kP9e/+mescnM9kt6T9LXJW1KOi/pcXf/69xrfijpy+7+AzM7Jumb7v6YmX1R0suSHpL0OUlvSvrC7G17brMWD37lS37hD69G/zl0YX3SkMIXnVhp2F1Hol43mZkf2NducdMb257FNVynq88uf0AALONynpsJjtgEDOOStixKbHpI0iV3/7skmdkrko5Kmk/Qjkp6avb5a5KeMzObff8Vd78u6bKZXZptTy22iUhqTzRy6XSo/d8hB6H6tlJ1LmTd3glM0Dnpxmd6vM9YnB1ALNt+5dyHH53sE5tuN7MLc1+fcfczc1/fI+lfc19vSnp4YRs3X+PuN8zsmqS7Z9//08J775l9vm6bVXjnLxc/sLuOJO1yOJnyhxUm5bnh32Gt+2L/gBDP9M26IX6puc4FM9uI0blA0gck5O6PjL0PALAoYmxa1v+0eKW06jWrvr+vxTarQJcDMJ5Az/S16YYIgqQPAADEsinp3rmvD0l6f8VrNigmAosAAAENSURBVM3sgKSDkq6uee+6bQJATOckte2O2Ksjok03RBAkfQAAIJbzkg6b2f2Srkg6Jun4wms2JD0h6W1Jj0p6y93dzDYkvWRmv1AzkcthSX9Wcwdw3TYBIJqA3RFtuiGCIOkDAABRzJ7RO6GmKr5f0vPuftHMnpZ0wd03JJ2V9OJsoparapI4zV73qpo2pxuSfuTuW5K0bJupjw0AAmjTDRFEpyUbAAAAAADDzVra35P0NTWdC+clHY9RyOJOHwAAAAAktqobIsbP4k4fAAAAAEzYsmmPAQAAAAATQdIHAAAAABNG0gcAAAAAE0bSBwAAAAATRtIHAAAAABNG0gcAAAAAE0bSBwAAAAAT9n+q4X1ZOEVzjQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, axs = plt.subplots(1, 3, figsize=(15, 5))\n", "for ax in axs:\n", " ax.set_aspect(500/440.)\n", " ax.tick_params(axis='both', which='both', bottom=False, left=False, labelbottom=False, labelleft=False)\n", " ax.pcolormesh(land[imin:imax, jmin:jmax], cmap=waterland_cmap)\n", "axs[2].scatter(xx-jmin, yy-imin, s=10, c=days, cmap='copper', vmin=0, vmax=7);\n", "ax_cbar = fig.add_axes([0.9, 0.12, 0.01, 0.76])\n", "cb = colorbar.ColorbarBase(ax_cbar, cmap='copper', norm=cnorm);\n", "cb.set_label('Days to Beaching')\n", "colours = axs[0].pcolormesh(xx-jmin, yy-imin, \n", " oilwaterhours, cmap='inferno', norm=mplcolours.LogNorm(0.004, 4));\n", "cb = fig.colorbar(colours, ax=axs[0])\n", "cb.set_label('Oil Exposure ($\\ell$ hr Oil / $m^2$ Area)')\n", "colours = axs[1].scatter(xx-jmin, yy-imin, s=10, c=beachvol/(440*500)*1000, cmap='plasma', vmin=0, vmax=0.2);\n", "cb = fig.colorbar(colours, ax=axs[1])\n", "cb.set_label('Beached Oil ($\\ell$ Oil / $m^2$ Area)')\n", "fig.suptitle('1,000,000 $\\ell$ Oil Spilled on Jan 15, 2018 at Salmon Bank');" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "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.7.1" } }, "nbformat": 4, "nbformat_minor": 4 }