{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import netCDF4 as nc\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "f=nc.Dataset('/data/eolson/results/MEOPAR/SS36runs/CedarRuns/t15r15DiatAlpha0/SalishSea_1h_20150220_20150410_ptrc_T_20150312-20150321.nc')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "mesh=nc.Dataset('/ocean/eolson/MEOPAR/NEMO-forcing/grid/mesh_mask201702_noLPE.nc')" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "<matplotlib.collections.QuadMesh at 0x7f3f51ff0c50>" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD9CAYAAACyYrxEAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnX/QJVV55z8PgzKOiuOg4DiMC0Q2G0tBkWV0TVlsiAFm3YxY/MpWDLBmpxLF0rWsZUyqlE1tsri1m92xTDBvgjpjXAFHLKgKyagoa1IlKIwO4hLCiBN4nQkjP6NFFp3k2T+6L3Pnzr19T3ef7j7d9/upeuvet2//eO6553z79HOe8xxzd4QQQgyXo7o2QAghRLNI6IUQYuBI6IUQYuBI6IUQYuBI6IUQYuBI6IUQYuAECb2ZvdfM7jWz75rZ+/Jta8zsS2b2QP764ny7mdlHzWyPmd1jZmc0+QWEEEIUM1fozezVwH8AzgJOB95qZqcCW4Db3P1U4Lb8f4DzgVPzv83AtQ3YLYQQIpCQHv3PAXe4+9PufhD4P8AFwCZgW77PNuBt+ftNwHbPuANYbWZrI9sthBAikBChvxd4s5kdZ2argI3AeuAEd98PkL8en++/Dnh47PjlfJsQQogOOHreDu5+n5l9BPgS8GNgN3Cw4BCbdpojdjLbTObaYQUrXr+KY4MMFkIIkfEjnnjU3V86b7+5Qg/g7tcB1wGY2e+R9dIfMbO17r4/d80cyHdfJuvxjzgR2DflnEvAEsCxtsY32DkhprTGzn27g/Y79+WnN2yJCGHW76Xfpxv0e7TDl33H34bsFyT0Zna8ux8ws1cAbwfeCJwMXAZck7/enO9+C3ClmV0PbACeGrl4+sS5Lz89SOx37tutyhuZ0JvsJPoduqPqbybawUKyV5rZXwLHAT8F3u/ut5nZccCNwCuAh4CL3P1xMzPgY8B5wNPAFe5+V9H5U+zRT6NuZZYQHYl6fv0hlpjrt43Hl33H3e5+5rz9goS+afoi9NNQ77M8k2W2yGXRFk33uPUbdkOo0Ae5bsRsQl08k8RseH1pZNO+c19sF4fQb9Y/1KOPTBuuiJCbREqNcdzelOxaROQqGxZy3XRMFw1qdM2UGq1EPg2KOgf6XfqLXDeJ0mRvvKobqQ4KQ+0HRXVDYybDRz36lpnlp47pv24r5FNugH4y7+as368/yHXTE5oKWYst9hpI7Tdl65l+234g101PKGpQdR+pq4jzrGNSafhyM4Sj0F8xQj36npCKm6QrodWg7nzKCLvKcBjIdTNwUnsUjy3EchWVQ4Pii4lcNwOnrAumyQG4stcS7SOBX2wk9AMlpGGPC3IdcW5jbkDVfYYqcIv4nUV15LoRz5LKOMA0huh/buLJpy/fXcRBrhtRmnGRmNbbT3Hgt83B4Vghq7EFXuIu5qEevZhLzJ5+7J5504O2bY83SLRFGdSjF9EYiU9Tg651xG2WbXVoU9wl7KINQhce+Y/Ar5Ot/fod4ApgLXA9sAbYBbzD3X9iZscA24HXA48Bl7j73qLzq0ffP+r08puMiS9yM81zQQ1xHEAMm2hx9Ga2Dvgr4FXu/g9mdiNwK7ARuMndrzezjwO73f1aM3sXcJq7/4aZXQpc4O6XFF1DQt8fYufkqXP8vPNOEuLjL0LiLlIjtuvmaOB5ZvZTYBWwH/gF4N/ln28DrgauBTbl7wF2AB8zM/MUBgNELfoyiDg5qFzlOqHHKCWD6ANzhd7df2Bm/51sXdh/AL4I3A086e4H892WgXX5+3XAw/mxB83sKbL1Zh+NbLtomTKpbutepy5tRQpNlokWixcpMlfozezFZL30k4Engc8B50/ZddRjt4LPxs+7GdgMsJJVgeaKFGk6/UHMWbvzqJs4TiIvUiTEdfOLwPfd/YcAZnYT8K+A1WZ2dN6rPxHYl++/DKwHls3saOBFwOOTJ3X3JWAJMh993S8imickl34M0S8657R959k475iqSORFXwgR+oeAN5jZKjLXzTnAXcBXgQvJIm8uA27O978l///r+edfkX++/xSJWtOTiMqcf9akr3nnmneD0mIdos+E+OjvNLMdZCGUB4FvkfXE/wy43sz+S77tuvyQ64BPm9kesp78pU0YLtqhrSySs+Lh24rICd1Xgi76iGbGipl0lfsm9oIrdW8ifXbR6GY1bJSPXtQiBYGYFTkzGdlSJMRVXC59XJmpis0S/P4joRfBTApqV0nMJonxRBEj7LPMco/z9q9KrIldVc/TVFpkpVuuh4RelCK1wcYY0Tttr7rUxYIrXY5hxEbCXh4JvahFUz3V8aeFkQumbM+9aoTMtB5q14umzEPiJ4qQ0Ita1HWbNDVRqegGVPappM1ZrFVuLH0eBBbtIKEXtagz8SiGP3me7zbG4GMbSdWqCnuIf1zCLyT0IhpVBsz64uuNJfZlxb1smUrsxTS08IiIQtdhlmVEs+mcOKE3tWkuoa7LUSw2EnpRinni1PfIkzopjYvy8zRxE9KNQoQi140opGq+mLLMc1eMR+lUpa0Q0hgirwVRRAhy3YjatCXyVW2JSdOJ2WKcQ4IuqqIevZhKqM+6TwtpV/XDdzmwPEIiL6ahHr1onLYEMMaC46G2StTFEJHQiyQJjSOf9lmXA8Kp5AkSYhwJvUiSELGOOfg7KzVDmZDKaecVIgXkoxeHkYLrYpIyaYernHPy3BJr0Rei+ejN7GeBG8Y2nQJ8CNiebz8J2Atc7O5PmJkBW4GNwNPA5e6+q+wXEN0QmjemK6qsETvvHJPbJfhiaJTq0ZvZCuAHwAbg3cDj7n6NmW0BXuzuV5nZRuA9ZEK/Adjq7huKzqsefbqkJPgxb0IScTEEmoq6OQf4nrv/rZltAs7Ot28DbgeuAjYB2/MFwe8ws9Vmttbd95e8lhCHkdJNR4g+cVTJ/S8FPpu/P2Ek3vnr8fn2dcDDY8cs59sOw8w2m9ldZnbXT3mmpBli0Tn35ac3nttGiKEQ3KM3s+cCvwx8cN6uU7Yd4R9y9yVgCTLXTagdol1S9dmHzNqdZ3Ob+eiF6JIyrpvzgV3u/kj+/yMjl4yZrQUO5NuXgfVjx50I7KtvquiStmfBVqGKfRL7YaOB9Ywyrptf4ZDbBuAW4LL8/WXAzWPbf80y3gA8Jf/8MKjrLkmVnft2J38TE+Upyia6aARF3ZjZKjK/+ynu/lS+7TjgRuAVwEPARe7+eB5e+THgPLLwyivc/a6i8yvqpl8MvdEM8Wa2yAx5SUatMCUaZehiD8MShEVnqDn9JfSicST2om8MLQ20hF60Ql/FvkqGyz4KgZjNEHr5oUJfNo5eiF4zbUA5dDC2rzc1MZ1ZwQVD/J3VoxfRGGIDmUafenyimL67ctSjF63Th4YhxDhFIcND6rgoH704jHkrNc0j5YlVMR7TdTMbDqHuuiH85hJ6cRjjQh0jzUAqxGisQ2jwqZB6vZlcA6Hvv7189GImqTfGWZRplDHz2Yv5VI10qTPpqe9++CK0OLgQAVRJ2lbXvbWo1BHrOuWs30hCLwZIlURWVZcrHMJjfQx080sbCb0YLFXFp+yA8pBzqUyjiiukr27AoSAfvZhKSGPuY+NNZbGSvtwMYs4WDhncF+VQCgRRmTINso9iPyKGyDSxfm1XTwizXF5NDGaGlJtuAvOR0ItKVBWZPgs+pCH6VZ+UmhTbJq5bdH2Jezkk9KIWZfzbEvlD9L0sQqlSZkoaF5+o4ZVmthr4E+DVZOu//nvgfuAG4CRgL3Cxuz+RLzyyFdhItvDI5e6+q8J3EB0x2SCHvhxbyt+vq956EaFuvHmD2imW91AJjbrZCvyFu1+YLxK+Cvgt4DZ3v8bMtgBbgKvI1pY9Nf/bAFybv4qeUDXUcNY5prEoPd+6VJ1gNO03LCuss64d+tvJNZMOc103ZnYssJtsGUEf234/cPbY4uC3u/vPmtkf5e8/O7nfrGvIdZMeZQSm7uBhSqKfgq9+nDaFsWyPfBoS8naJ6bo5Bfgh8EkzOx24G3gvcMJIvHOxPz7ffx3Z+rIjlvNtWiB8AEzz3U8OIs7rQU6eI6VEaF1PgGr62kVjL9OuXRSBI1HvDyFCfzRwBvAed7/TzLaSuWlmYVO2HfHYYGabgc0AK1kVYIZok6JGXJT0bNp+dfdpmy7FLEZoZdE56n6fGC4h0T4hQr8MLLv7nfn/O8iE/hEzWzvmujkwtv/6seNPBPZNntTdl4AlyFw3Fe0XHRDDh98Xqoh+zCeUUOHv4uaU0pOYKCYovNLM/hL4dXe/38yuBp6ff/TY2GDsGnf/T2b2b4AryaJuNgAfdfezis4vH32/WaTGnvL4g3rWi0fsFabeA3zGzO4BXgv8HnAN8BYzewB4S/4/wK3Ag8Ae4I+Bd5W0XfSMRRKYVAcnF+k3EOXRhCkRjT707GOncuh65rB85ouNZsaKzkhd8GOLc9diX9YOMRwk9KJThiL2I8bDLkfvy05m6mK2qhg2EnqRFKkKf6zwxdColybKQSK/uGgpQZEUqS4qXiXFwLzw0irpfoVoEvXoRef0SQBDUkB0gXr1i4lcN2IQpCCi00hxlS2J/eIh140YBKm6P7qwQ3lnRFUk9KJXLFL6hTIs2gLlohwSetFbFi3XSuy4frE4SOhFrxniouV16TrVskgPCb0YJIvW259E67OKcST0YrCkOpCbAroRNEtqq4xJ6MXCkGoMfIrI/VONMouhtzmALqEXC42ieKYjkS9HlbrTZhlL6MXCMK0xTltmT4I/v6zapo+/SUo3yyChN7O9wI+AfwQOuvuZZrYGuAE4CdgLXOzuT5iZAVvJVph6Grjc3XfFN12Icsh1U4/U8t03sdpX3UH8VMpmkjI9+n/t7o+O/b8FuG1sKcEtwFXA+cCp+d8G4Nr8VYjkaDvT5BBoa0ZulYRzRVRJLR1yzj5Qx3WzCTg7f78NuJ1M6DcB2z1LonOHma0eLSJex1Ah2kZRO/NpMw1DFzeVNq7dBqFC78AXzcyBP3L3JeCEkXi7+34zOz7fdx3w8Nixy/k2Cb3oNRq4Laau6IdErFR1v1R1yfRd4EeECv2b3H1fLuZfMrO/LtjXpmw7IkWmmW0GNgOsZFWgGUKkgXr7xcQuj1CBn/a+KimKfNWw1yChd/d9+esBM/sCcBbwyMglY2ZrgQP57svA+rHDTwT2TTnnErAEWZri0pYL0TES97SI9XukKPB1mSv0ZvZ84Ch3/1H+/peA3wFuAS4Drslfb84PuQW40syuJxuEfUr+eTFE1KtvlxABrvsbDFHkIaxHfwLwhSxqkqOB/+3uf2Fm3wRuNLN3Ag8BF+X730oWWrmHLLzyiuhWC5EIEve41BHaIYZFxmKu0Lv7g8ARpeDujwFHLAuVR9u8O4p1QgwEZdmcjcS9eTQzVoiKlB0gnPX5oot9SGx72TIaL9dFEvRZSOiFqEhILz1EbJR6YTryt8dDQi9EZEJn207up8HdeKg3fzhJCP0/P+1pdu7UmpdiuMwS7aK4aE3Qqk9qydm6Igmhn4Z+IDEUikS+6HPRDG2mbUiFZIV+GloMQQyNWdP1x9GN4EhiPe0siounV0IPEnsxLKaJ1GQdV2TO4cQa21gkHemd0IPEXgyf8ToukT+ceb1wacOR9FLohUiRpgRZPXpRl94K/aL41kR6zBrMiynIqtfzUcBGOL0V+hFy44guKSPsZZYyVEdmNmVWBQsZ7B4/bqjl3XuhB4m9aJcqER9V66fqdkZRmc9LdxBSfuPjIUMs70EIPQz/jiz6TZ0Vj0bis4h++lDhjtXuR+U873xF+6QYp39U1wbEZue+3c/+CdE0ZRvyuS8//dljQo8tu/9Q6PL7TstXNM6sG+/ktlT0yLKswt1y5ukr/Rs7X9HoNRatkYj2qOqnb+O4oZBS+w3Jtlm0f920zOPHf9l33O3uZ847Lth1Y2YrgLuAH7j7W83sZOB6YA2wC3iHu//EzI4BtgOvBx4DLnH3vWW+TBPELGwhxinrhql7jqEz7aZWNyX0PHdKGb9/2RTIZQbhm9Kl4B69mb0fOBM4Nhf6G4Gb3P16M/s4sNvdrzWzdwGnuftvmNmlwAXufknRudvo0c9Dwi/qUFWkQyJAZonKUG8MZcpk8pgyAjqv/OaVdwqaEdqjDxJ6MzsR2Ab8LvB+4N8CPwRe5u4HzeyNwNXufq6Z7czff93Mjgb+DnipF1woBaGfJIUfUfSTGAuS1BWhPpNC2+tLjH5sod8B/FfghcAHgMuBO9z9lfnn64E/d/dXm9m9wHnuvpx/9j1gg7s/Ouv8KQp9DFKsGKIdYoh91fNI9BeHaD56M3srcMDd7zazs0ebp+zqAZ+Nn3czsBngFev6FeWpyifmEZKnJiSML8R3P+nWkL9fTBKisG8CftnMNgIrgWOB/wWsNrOj3f0gcCKwL99/GVgPLOeumxcBj0+e1N2XgCXIevR1v0ibpBgnK9IjVj0JFe5xwR9ClE5d91boNWb9TkOaPFUqvDLv0X8gH4z9HPD5scHYe9z9D83s3cBrxgZj3+7uFxeddyium6FUClGeItEpWy+qDr6GTPbpi+jHcmvFItW2HdVH/+zOhwv9KRwKr/wW8Kvu/oyZrQQ+DbyOrCd/qbs/WHTeY22Nb7BzDtvWlwpZRKqVQzRP3QiNOvW/70Lfpcj3rc02IvRNMU3oR/ShYpahbxVJlHfBTOuRV3EDhA7Elh2w7WubamPsoW/tU0KfKH2rSOIQ86bCN/XbhrqFNNO2OVJtt4MR+iL6XkFTrTxiPm3PbCy6ZghtDGwuCim124UQehhGBU2p4ohwYoh9VbdQEWWyKg6h/XRBKm12YYR+nKFU2lQqkSgmVqTNpADPm+4/j9DkWqKYee6tFNrpQgr9iKFV6BQqlDiSFOuZXDTlCW1fKcbVL7TQjxhipU6toi0C4w08VhbFovPVdatMi/ap+/QxxLY0jb61Lwl9ztAraN8qZh+J5TqJfd7QxGfTeqJ1XBFDb1OTpNzGJPRTGHIFTbky9p2qYYtlztdGPHzZVL1lv0vfBniH0GYk9AX0pSLWYQiVuGtiZaAse07RLn1uKxL6Egy98fW5IndB2Vzwdc4Vmz71qlPKstnXNiKhr0Eqla8J+lqhY1LFFVNl9muZVZBiCvR4crMh1+Um6Fv7iL5m7CIx5AbS9HT9vlD2+1cpr6rL3tWljbobu41MGzgeahvsAvXoKzDkCjjUG0BXMdChdaUvq0QNdTJWX+u9XDct0fcKXkRfK3/XpFonuuwpx5ob0MRTRJ+R0LdMqo07BkNoEG1RpR5UcYMUhVmWDaOsS6z6kUob6lN9jyb0+UIiXwOOIfPp73D3D5vZyRxaeGQX8A53/4mZHQNsB14PPAZc4u57i64xBKEfJ5UK2yZ9ahxN0OVvXjZVctWbUR3KTNDquiz7REyhN+D57v5jM3sO8FfAe4H3AzeNLSW4292vNbN3AaeNLSV4gbtfUnSNoQn9iEUU/En61nDq0tZvXrVcq8wNmDVAWmdyWAh1y3Lek80Q6mZTSwmuIhP63wT+DHiZux80szcCV7v7uWa2M3//9Xxx8L8DXuoFF5LQD48+NKImI5D6sBLSvFmxTcwnaPIGNYs2o5/arvdRwyvNbAVwN/BK4A+A7wFPuvvBfJdlYF3+fh3wMEB+E3gKOA54tNQ3GACLvIpPm2GcoblchkCbQhKr/o4fWzX6qQlRnib8Q603ZXv0q4EvAB8CPunur8y3rwdudffXmNl3gXPdfTn/7HvAWe7+2MS5NgObAVay6vU/bxtjfJ9eMNTKVERoQ00hzW4dMW0jDjy22JdZQGVW7vwiu6qePxazOgFNdMRS7dGXjroxsw8DTwNXIddNFBZR+EekIOzzCEk73DZNu24mrzEujrFm3bbZg27iiS8F92TMwdiXAj919yfN7HnAF4GPAJcBnx8bjL3H3f/QzN4NvGZsMPbt7n5x0TUWXeinkYLAifo0LWZ1kqpVTXs877whx3aVAyjWtVMQeYgr9KcB24AVwFHAje7+O2Z2CofCK78F/Kq7P5OHY34aeB3wOHCpuz9YdA0JfRgSf1FEk08eMUW+rE2hcwa6vHl0RbTBWHe/h0y0J7c/CJw1Zfv/Ay4KtFOUoOtekUibGAPgsdMu1xl4jXH9mKQg7FXRzNiBo5uBGDHL516WMoOr844pe57Jc7Xdo09N7JUCQQShG4EYUWUiU1HkTdtjE6lPVmsCCb2Igm4Ei0eoiLadJC1meoe6g+SpiL2EXkRFgr84jIdQFtFU7HmbGS6rPIGkIvIgoRcNIcFfXJoKUQzJvjmNpnLhpCTk89AKU6IR2mqEIj1i/6ax0irEpE8iXwYJvWiEIecNEXGpGtJZx+UzVEGfhYReNEaKqQNEGrQZ6x/r2D4jH73oHIn+4jDrSW/SXx/SWw/pSIyfd4gir8FY0Wsk/sMixmzXIQp1XTQYK3pNisvMiXBii7JEvh4SetErmlgTVQPH1WgrL76oj1w3QgSyyELUdI9aETPVkOtGCFGJLoS0zGSsoQ6sNomEXohA2s7t0jR1xbJo0DRkicG615bYhyPXjRAV6LPYhyQta2IsJPS8Vc69qKIfc4Wp9cB24GXAPwFL7r7VzNYANwAnAXuBi939CTMzYCuwkWxt2cvdfVfRNST0ok/0WeTH6cN6vWVZNMGPKfRrgbXuvsvMXgjcDbwNuBx43N2vMbMtwIvd/Soz2wi8h0zoNwBb3X1D0TUk9KIPlO319kUoYy4+kgKLJPYxlxLcD+zP3//IzO4D1gGbgLPz3bYBtwNX5du3e3YHucPMVpvZ2vw8QiwMZVNA1BGougnCyl67bEhq1VWjJj8LIcaSikOj1GCsmZ1Etn7sncAJI/F29/1mdny+2zrg4bHDlvNthwm9mW0GNgOsZFUF04Xojphrq8YQpDophJsUxNBzF+1XNAguMQ8jWOjN7AXA54H3ufvfZ6746btO2XaEf8jdl4AlyFw3oXYI0RVNiEpXvfgy56lqYwyX0Pg5xp8iJPDlCBJ6M3sOmch/xt1vyjc/MnLJ5H78A/n2ZWD92OEnAvtiGSxEH5nV4y7jNknJL97EE8q0c6oXH4e5Qp9H0VwH3Ofuvz/20S3AZcA1+evNY9uvNLPryQZjn5J/XohDTPq36/RSp/mzYy3NN8umFBcMEcWE9OjfBLwD+I6ZfTvf9ltkAn+jmb0TeAi4KP/sVrKImz1k4ZVXRLVYiB5Tddm8KueetW0eddYRaHoNAvXmq6EJU0L0mGlPA7HzxsR6ImhygtWiEhpeeVQbxggh4rFz3+5oKQZC8spMcu7LT3/2r8z5Q48pu6+Yj3r0QvSIaa6fqj3uqhO+qvjty4wbSODDUY9eiIERywcP5UV+Xi++yvlEeyh7pRA9oa6/e5xZYZ1NDxZL4LtBQi9ET4gdmjhtILeOEMe8EYm4yEcvRA+IKZxt96oVbdMc0bJXtoGEXoj5hAyOxu5JN7U4SRvXXgQ0GCtEDUYhjLOEquizpuwJ+Sx2WGJTs2BDQzPlAoqDfPRCjBESn96Ur3zedWeFUk7bP0YahCZ71GXsKRPmKaYjoRciZ1b4YhmBjXXN0H2bIoWsmuNI1OshoRdiBiPBKpN0rK1B06ZDFbv6zpNI4OMgoRcLQ1Pui9DzxhzYbCu18bRrtvV9RTwk9GLQtC3uKazWBM30sptaoEQ0j4ReDI6m1jItc1yXtOXLT7kMxOFI6EVvKBoQbToSZvxaqdKGwKdeBmI6IStMfQJ4K3DA3V+db1sD3ACcBOwFLnb3J/LVqLaSLTzyNHC5u+9qxnQxFJqI1Q7NsCh3xCEW6bsuGiE9+k8BHwO2j23bAtzm7teY2Zb8/6uA84FT878NwLX5qxC1e5x1hahPQtZGtsc+lYeox1yhd/evmdlJE5s3AWfn77cBt5MJ/SZgu2d5Fe4ws9WjBcRjGSzSJ0TQF0lkZoUqzts+b1sVFqncxSGq+uhPGIm3u+83s+Pz7euAh8f2W863SegHSh8HK9tk2iLg48yb1aoUACIGsQdjbcq2qVnTzGwzsBlgJasimyGaQKsChRHL7ZJagjLRX6oK/SMjl4yZrQUO5NuXgfVj+50I7Jt2AndfApYgy15Z0Q4RGYl5feos7wf9Tkks0qSq0N8CXAZck7/ePLb9SjO7nmwQ9in559NFbpdmKJM+YPKYqug3E0WEhFd+lmzg9SVmtgx8mEzgbzSzdwIPARflu99KFlq5hyy88ooGbBYV6GNceF8o66qZTCsQc4KXENPQwiMDRcJejiYWyJjXs29qYW+xOIQuPKKZsQuExH82scuhyWiZKsnNxGKjHv2AmZfKNsSXHMt3vCgzULsIhxxK2YnyaM1Y0RmL2OOcN7Adyw+vAXQxjoReiI5QBI1oC/noRS+oEoqYKrOibxRVI7pGQi86oa8CX+Q6Cc2YOet4IZpCQi9ao895WxZx3EEMBwm9CGJeaOa81AltpN2NTV2by9zYml7sWyw2Enoxk3HxqROC2fai1rGSiLWRjEziLtpAQi8OY9x3XnU91XFi53tpShhjuGbK3pwk8qItJPTiCLdBmwtmpyLuRdcJ3a+p9WuFqIuEfoEpWuGoTH6W1HvtVVIv1+nha+EQkRqaMNVz6rgLJgU9VNxi5cwpI/pVByunfacyA8tNibbcNiIGmhk7MNoaKAxZv7SuHXVvTmVtKLM2a5NI3EVsJPQLTNnFqEPONTqujNCWdX+U8f1P+y51bgJljwtFYZOiSToVejM7D9gKrAD+xN2vKdpfQl+NUN9znfDGkF5w7KeAsteoK6RdDggLUYfOhN7MVgB/A7yFbA3ZbwK/4u7/d9YxEvrZxPBNjyjTq667WEaTqY/nnb8qfU3LIBaXLpOanQXscfcHAfL1YzcBM4VezCbWAGST1ytzbBcToVI5vxBd0YTQrwMeHvt/mWyhcNEiTYlWUeRKURRP0fnUkxaiWZoQepuy7Qj/kJltBjbn/z7zZd9xbwO2xOYlwKNdGzGPFWt5CTzQqJ0r1s765IGAfZ7lJfDAo4f2e6Bo3y7pxe+O7IxNH+z8ZyE7NSH0y8D6sf9PBPZN7uTuS8ASgJndFeJn6hoxa7JgAAAErElEQVTZGRfZGRfZGZe+2BnCUQ2c85vAqWZ2spk9F7gUuKWB6wghhAggeo/e3Q+a2ZXATrLwyk+4+3djX0cIIUQYjeS6cfdbgVtLHLLUhB0NIDvjIjvjIjvj0hc755LEzFghhBDN0YSPXgghREJ0LvRmdp6Z3W9me8xsS9f2jGNme83sO2b2bTO7K9+2xsy+ZGYP5K8v7sCuT5jZATO7d2zbVLss46N5+d5jZmd0bOfVZvaDvEy/bWYbxz77YG7n/WZ2bot2rjezr5rZfWb2XTN7b749qTItsDOpMjWzlWb2DTPbndv5n/PtJ5vZnXl53pAHa2Bmx+T/78k/P6ljOz9lZt8fK8/X5ts7a0u1cffO/sgGa78HnAI8F9gNvKpLmybs2wu8ZGLbfwO25O+3AB/pwK43A2cA986zC9gI/DnZ/IY3AHd2bOfVwAem7Puq/Pc/Bjg5rxcrWrJzLXBG/v6FZCk8XpVamRbYmVSZ5uXygvz9c4A783K6Ebg03/5x4Dfz9+8CPp6/vxS4oaXynGXnp4ALp+zfWVuq+9d1j/7ZdAnu/hNglC4hZTYB2/L324C3tW2Au38NeHxi8yy7NgHbPeMOYLWZzZ/K1Jyds9gEXO/uz7j794E9ZPWjcdx9v7vvyt//CLiPbIZ3UmVaYOcsOinTvFx+nP/7nPzPgV8AduTbJ8tzVM47gHPMbNrEy7bsnEVnbakuXQv9tHQJRRW3bRz4opndnc/kBTjB3fdD1vCA4zuz7nBm2ZViGV+ZP/p+Ysz1lYSdudvgdWS9u2TLdMJOSKxMzWyFmX0bOAB8iexp4kl3PzjFlmftzD9/CjiuCzvdfVSev5uX5/80s2Mm7cxJoS0F0bXQB6VL6JA3ufsZwPnAu83szV0bVIHUyvha4GeA1wL7gf+Rb+/cTjN7AfB54H3u/vdFu07Z1pqtU+xMrkzd/R/d/bVkM+PPAn6uwJZk7DSzVwMfBP4F8C+BNcBVXdtZl66FPihdQle4+7789QDwBbIK+8jocS1/PdCdhYcxy66kytjdH8kb1z8Bf8whV0KndprZc8jE8zPuflO+ObkynWZnqmWa2/YkcDuZT3u1mY3m7ozb8qyd+ecvItzlF9vO83IXmbv7M8AnSag8q9K10CebLsHMnm9mLxy9B34JuJfMvsvy3S4Dbu7GwiOYZdctwK/lEQNvAJ4auSO6YMKneQFZmUJm56V5BMbJwKnAN1qyyYDrgPvc/ffHPkqqTGfZmVqZmtlLzWx1/v55wC+SjSd8Fbgw322yPEflfCHwFc9HPzuw86/Hbu5GNo4wXp7JtKVSdD0aTDaS/TdkPrzf7tqeMbtOIYtY2A18d2Qbme/wNrJUi7cBazqw7bNkj+g/JetlvHOWXWSPm3+Ql+93gDM7tvPTuR33kDWctWP7/3Zu5/3A+S3a+fNkj+D3AN/O/zamVqYFdiZVpsBpwLdye+4FPpRvP4XsRrMH+BxwTL59Zf7/nvzzUzq28yt5ed4L/CmHInM6a0t1/zQzVgghBk7XrhshhBANI6EXQoiBI6EXQoiBI6EXQoiBI6EXQoiBI6EXQoiBI6EXQoiBI6EXQoiB8/8Bz4OMQpRel0gAAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 432x288 with 1 Axes>" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(mesh.variables['tmask'][0,0,:,:])" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "<matplotlib.collections.QuadMesh at 0x7f3f5162a588>" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD8CAYAAAB5Pm/hAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAG6BJREFUeJzt3X+sZHd53/H3Z+beXf/A1Ha4Rlv/qI3rQGjU2ObWOLiJqIHUuCk2bVLhRInVWNpECpVR0hYTpAbUVgptwWmliGipXTYVPwNYtmjS4jqmFkkw3TXrX13AxjGO8ca7GLs2BdbcOU//ON9zznd+7Z29M7Mz9/B5SaOZ+c6Zc545d/z47Dzf8xxFBGZm1l6dRQdgZmbz5URvZtZyTvRmZi3nRG9m1nJO9GZmLedEb2bWchMnekldSV+W9Nn0/AJJ90p6RNInJO2YX5hmZrZVx3NEfyNwMHv+PuDmiLgIeBa4YZaBmZnZbEyU6CWdA/wD4D+n5wKuBD6VFtkLXDuPAM3MbDorEy73u8C/BE5Lz38EeC4iNtLzJ4GzR71R0m5gN8Cpp576mle96lVbj9b42v7HFh2C2bb2o695xaJDOG779+//VkSsbfX9myZ6ST8LHI6I/ZJeXw2PWHRkL4WI2APsAVhfX499+/ZtMVQDeFPn5xcdgtm2due+P1x0CMdN0jemef8kR/RXAG+RdDVwEvBSyiP80yWtpKP6c4CnpgnEzMzmY9Pf6CPiXRFxTkScD7wN+JOI+EXgbuDn0mLXA7fPLUozM9uyaebRvxP4DUmPUv5mf8tsQjIzs1matBgLQER8Hvh8evwYcNnsQzIzs1nymbFmZi3nRG9m1nJO9GZmLedEb2bWck70ZmYt50RvZtZyTvRmZi3nRG9m1nJO9GZmLedEb2bWck70ZmYt50RvZtZyTvRmZi3nRG9m1nJO9GZmLedEb2bWck70ZmYt50RvZtZymyZ6SSdJ+pKk+yU9LOm9afzDkv5C0oF0u3j+4ZqZ2fGa5JqxR4ErI+I7klaBL0j64/Tav4iIT80vPDMzm9amiT4iAvhOerqabjHPoMzMbHYm+o1eUlfSAeAwcGdE3Jte+reSHpB0s6Sdc4vSzMy2bKJEHxG9iLgYOAe4TNKPA+8CXgX8HeBM4J2j3itpt6R9kvYdOXJkRmGbmdmkjmvWTUQ8B3weuCoiDkXpKPBfgMvGvGdPRKxHxPra2trUAZuZ2fGZZNbNmqTT0+OTgTcCX5G0K40JuBZ4aJ6BmpnZ1kwy62YXsFdSl/J/DJ+MiM9K+hNJa4CAA8CvzTFOMzPboklm3TwAXDJi/Mq5RGRmZjPlM2PNzFrOid7MrOWc6M3MWs6J3sys5ZzozcxazonezKzlnOjNzFrOid7MrOWc6M3MWs6J3sys5ZzozcxazonezKzlnOjNzFrOid7MrOWc6M3MWs6J3sys5ZzozcxazonezKzlnOjNzFpu00Qv6SRJX5J0v6SHJb03jV8g6V5Jj0j6hKQd8w/XzMyO1yRH9EeBKyPiJ4CLgaskXQ68D7g5Ii4CngVumF+YZma2VZsm+ih9Jz1dTbcArgQ+lcb3AtfOJUIzM5vKRL/RS+pKOgAcBu4Evg48FxEbaZEngbPHvHe3pH2S9h05cmQWMZuZ2XGYKNFHRC8iLgbOAS4DfmzUYmPeuyci1iNifW1tbeuRmpnZlhzXrJuIeA74PHA5cLqklfTSOcBTsw3NzMxmYZJZN2uSTk+PTwbeCBwE7gZ+Li12PXD7vII0M7OtW9l8EXYBeyV1Kf/H8MmI+Kyk/wN8XNK/Ab4M3DLHOM3MbIs2TfQR8QBwyYjxxyh/rzczsyU2yRG9WTvoGL9URnHi4pjUseIdtIzx29JwCwQzs5ZzojczazknejOzlnOiNzNrORdjrfXU7QIQxciTt9NC2+eYRx2NGO2OWnDofcWLL84nKFtq2+fbbWZmW+JEb2bWck70ZmYt50RvZtZyTvRmZi3nWTfWSp2dJ5UPstYAGpiYMnYWziLaCUw562f0TJzhdXdOPjmNNcsX3/3uVNu25ecjejOzlnOiNzNrOSd6M7OWc6I3M2s5F2OtVbpnnlE+6KWCalFAlEXXKNJYei6NK8Yeo7B5LFspqB6riAp9RdPxmx2z3QneC80+63372YmWt+3HR/RmZi03ycXBz5V0t6SDkh6WdGMaf4+kb0o6kG5Xzz9cMzM7XpP8dLMB/GZE3CfpNGC/pDvTazdHxH+YX3hmZjatSS4Ofgg4lB6/IOkgcPa8AzMzs9k4rmKspPOBS4B7gSuAt0v6ZWAf5VH/UDVH0m5gN8B55503Zbhmw1YuvKB88P2jZfEVoNcr7yOgOgO2OuO1fp4VY4stnA07qgg6WADtqNneJoVXDRZzxy1fbWNUsXWSmOyHzsTFWEkvAT4NvCMingc+CFwIXEx5xP/+Ue+LiD0RsR4R62trazMI2czMjsdEiV7SKmWS/0hEfAYgIp6OiF5EFMCHgMvmF6aZmW3VJLNuBNwCHIyID2Tju7LF3go8NPvwzMxsWpP8Rn8F8EvAg5IOpLHfAq6TdDEQwOPAr84lQjMzm8oks26+AIyq5vzR7MMxM7NZcwsE27b02r8NQPy/o+XAjlVUzaTZ6DX3RflYvf4WCOWMnIExaGapxIgWCfkMlmqGy7FmwYwMfNRsmRFj0nBbhY6G369sLF/PBPGEhKp9Za3lFghmZi3nRG9m1nJO9GZmLedEb2bWci7G2rb0wi/8JKccerF8csoqACqCztGysNh5caMc2yjqdgjqpeLqDzaaFfUGWh/kFwYf1V++kxVeRxVhO8NjMdQWYcyHKka8PtjSQBouEkvD2xi3vcFtSIR2lEOv+fFykf0+JaZtfERvZtZyTvRmZi3nRG9m1nJO9GZmLedEb2bWcp51Y9vK4X/2OgBOejY4emaabZNm03SPFnTTrJdYLY9htFGgHxRpuTT7ZrVbriwCev0zWBQxegZLdUiUzYKJwRkx3WYmTt86Rh1OjdrG4GyaYvR7q89bbSOPuWoBMXL7BbA6ELuaWCIbsnbxEb2ZWcs50ZuZtZwTvZlZyznRm5m1nIuxti0VK7Bxclk27NQdDTrl9c6A6KbXukLpsTZSgbabirNFMVQADRhqbdBX2OxmrQ2qdgfVKjqiWBl470r23tR+YFTxtHx9sDBMHZ/yTg3qb9sQWfm0Xt2otggd6kprtX8i71+fHrsY2z4+ojcza7lJLg5+rqS7JR2U9LCkG9P4mZLulPRIuj9j/uGamdnxmuSIfgP4zYj4MeBy4NclvRq4CbgrIi4C7krPzcxsyWya6CPiUETclx6/ABwEzgauAfamxfYC184rSDMz27rj+o1e0vnAJcC9wMsj4hCU/zMAzhrznt2S9knad+TIkemiNTOz4zbxrBtJLwE+DbwjIp7XhFe8j4g9wB6A9fX12GRxs2OqTtPv7RChgVYAG9BZLR/X81LUofqudqrZKtVXt8i+w9V0lWy2St+MmWqxbqdZfmDWTbHaqeOrZrPEiprtVZvI2i7kE2g0oh1D3e4gm5HTIYuhWm+9jWrWDyPbLER3cCZONjvHUzNaa6I/raRVyiT/kYj4TBp+WtKu9Pou4PB8QjQzs2lMMutGwC3AwYj4QPbSHcD16fH1wO2zD8/MzKY1yU83VwC/BDwo6UAa+y3gd4BPSroBeAL4+fmEaGZm09g00UfEFxh/stwbZhuOmZnNmlsg2LbwjX+d+tA/Uz4vVqgLj3UP9hUoeqnwWhVIaX6fjKhaGqRWCEUMdhMgOlkhtTv8y2ZdoA3qdgeK5rXBYmyxMlzwVVYE7uSF2aqVQ1W0LbLfVvPe8526z0G9XLOS5nM0YyPaNlTPu1nxORVqu1jbuM5uZtZyTvRmZi3nRG9m1nJO9GZmLedirG0PA63fowuRipBFKm4WXVCqJOYFyuoE0mIlFWGrAmgnhq7HjZozbYvVrHA5cAHt/HXV61dTyE1vLXY066hPzO01H6S+IHdeeK1i6jVn+Nb11B7NGblkb6gWyM7wrV/NFx94aygrIKfC8ffe+lpOvu1erD18RG9m1nJO9GZmLedEb2bWck70ZmYt50RvZtZynnVj20s1uaQL0UuPO819pFk31WvE6L7y5Uuq2yfUY/kslNVONl5tt1lXMbDeYiWblZNmvfRWm1k5VesDdWjaIVRtHDaa9SgPqXpv1eahq/o9ZJ0Qmjb0/bN+8liqz9f3eTvN6zF6N1kL+IjezKzlnOjNzFrOid7MrOWc6M3MWs7FWFt6j/7Bpaw+UT6ui51FXvhs7quxotv0pY+q4No0d6/XPdgCITqC7nAv+cH1KqC3o/+9RXYh8LoFwkpW8M36zNeF2bobQ9T98ima5ZW2V3c2KKK5GHpaYd5jfmRBNR8bPLSTmvdky73wCz8JwGkf/fMRK7Ttxkf0ZmYtN8nFwW+VdFjSQ9nYeyR9U9KBdLt6vmGamdlWTXJE/2HgqhHjN0fExen2R7MNy8zMZmXTRB8R9wDfPgGxmJnZHEzzG/3bJT2Qfto5Y9xCknZL2idp35EjR6bYnJmZbcVWE/0HgQuBi4FDwPvHLRgReyJiPSLW19bWtrg5s1J0sls33frGlG7V66JIt1gZcauWX+2UtxVRHOu2CsVqOeOmWBXFqujtKG/FKvRWRS+NlzfGrIu+W3SVfY50AZM0I2bo1lXfDXHMW0jHuOXLNTdrly0l+oh4OiJ6EVEAHwIum21YZmY2K1tK9JJ2ZU/fCjw0blkzM1usTU+YkvQx4PXAyyQ9Cfw28HpJF1OeAvI48KtzjNHMzKawaaKPiOtGDN8yh1jMzGwO3ALBtoW6QFj1is/G+vrRd4bHqtYDRTrHX3WDeGU93Zue7FXP+WKVobHIfuzsrdZN4MvXugy1Eyi6eT/6tFyAeuob00rUsXQ2hmOvWyGQ9aanf/tj5T/QDi4roGqhIFdh28otEMzMWs6J3sys5ZzozcxazonezKzlnOjNzFrOs25sab3sz8oWSo8+2oxVs14EzQyc+l4U3XJGSjWrpVgJOvTPmKlmugAoVL8Xylkz9WybbKZN0aVvOVS2Lujbfnc4ztDwrBuimehCr7wriv6Yy+XzC6SkMdRcLKUaC5qLkVQzcqS+x/V6BifWdJoLj4xqffDsr5QXIPn+GeWLu97/Z8ML2dLzEb2ZWcs50ZuZtZwTvZlZyznRm5m1nIuxtvyyImHegqB+nBdDq+Jqt2lzUNU0VaSx3nA7g/qxVBdZ8wJsdAeXy4q2WeGVgcJmdLI2B73m4xQDXQdUqIm5Lh5H9rgKNJrDs+qzVm+jvy1CjGppMDCWx9y3n8d0Q3jiva/jvN92QXa78RG9mVnLOdGbmbWcE72ZWcs50ZuZtZwTvZlZy3nWjW0PAxcZQRBpNkuRxjodIM2OKaqZLp289UBqE5C1KhhcbwiKlWrmDvV9NQOHbLnq9fy1wQuP5HGqasEQw0dY+TY62Vj9et+FV9Ln6EXzWtUCoYj+5Qe3M+LCI0MtEEbO1hm9PtseNj2il3SrpMOSHsrGzpR0p6RH0v0Z8w3TzMy2apKfbj4MXDUwdhNwV0RcBNyVnpuZ2RLaNNFHxD3AtweGrwH2psd7gWtnHJeZmc3IVouxL4+IQwDp/qxxC0raLWmfpH1HjhzZ4ubMzGyr5l6MjYg9wB6A9fX12GRxs5H6i5Hlaf+DLRCKLnTSN6wqqCqy1geDfeHzQmTVez4vslb95lf6WxpU91G9PtgeIa27ftjrj5MeVCFVQ0UXlGIoUqCdjuqgo1oHWZ/5bhVUtt3OcNW0il0BdDX8Wt2HQX3LW3ts9Yj+aUm7ANL94dmFZGZms7TVRH8HcH16fD1w+2zCMTOzWZtkeuXHgD8HXinpSUk3AL8DvEnSI8Cb0nMzM1tCm/5GHxHXjXnpDTOOxczM5sBnxtrS+adfewKA//bM6c3g4JmxMXxmqrpQ1VmzxeoLgA8WY6MzvN7oDK+3WKG/5/2Y5frOHq0eR3NGbLX9zkZTNK56ycdKf8G13MaIfvQERXVmbKroKiuwRi/yVfQJDRdaY0Tx1mfBto973ZiZtZwTvZlZyznRm5m1nBO9mVnLOdGbmbWcZ93Y0unWc2dKUqTz9yGqvuv094uH/hYERTbrpZpF0vSlz5avWyCk+7wFwmrzWgzMuoGsRUI966bpRVC3HShUb7dqhVBksVStEIimT34ni69eT73dpi1C3bKgiOYzjpgy09cCYUQ/+sHlBscHPfofLwfgb974xfEL2VLxEb2ZWcs50ZuZtZwTvZlZyznRm5m1nIuxti3EQNE078Ge94/vbPQvD/S1I+ij4SLrqD7zRZemBUJfP/roGxvZdqCIugir1G6gQ3Yh8KodQ9bSoWqLUHQH2jUMfqC6aBtNT/nYpAXCQMuDUfsp31a4R30r+IjezKzlnOjNzFrOid7MrOWc6M3MWs6J3sys5TzrxpbK3r/8U/70+7sA6JDNIBmcHZNNNKlnphQMNE9g9Kn82eySvPVB2mjT2qDa5kr0zbaplovO4KybGJqeooIm0NTvIAo1M3vyWTdpPfWMnG7T7qCZLZNNHare2837GFQbHj1NJr9wS7VcDM62GXGBEtvepkr0kh4HXgB6wEZErM8iKDMzm51ZHNH/vYj41gzWY2Zmc+Df6M3MWm7aRB/A5yTtl7R7FgGZmdlsTfvTzRUR8ZSks4A7JX0lIu7JF0j/A9gNcN555025OfthUPWj76T+7uo0xVC6w8XIopcqh6swWLPM+8dX+toYDLQviA4Uq+nNeTFWg8tFc5jUyYrGA9Xg6GVPqmJrAUoxj4pFeexRLVdXT5sPV8VSjCrGDn/ucj39rRL6WitUL41qH+Hi7LY21RF9RDyV7g8DtwGXjVhmT0SsR8T62traNJszM7Mt2HKil3SqpNOqx8DPAA/NKjAzM5uNaX66eTlwm8r5uivARyPiv88kKjMzm5ktJ/qIeAz4iRnGYmZmc+DplWZmLecWCLYUPvvN/QA8U8BqulJHJ/U26HYLNqp2A9mp/fX8m2ymSwweuuQTUrI2B9Xz+uIheQuE1f4rlMRKNOvpNK9Ft3/2i1ReaKR/8536WiB1wN3+VgpV3PUMnGjGBi94Um5LfatrNgCdY0yPCWWxFtmsn3omUfbewRk4I2bkFHedC8BZJ3+Hb73u2bHbtcXzEb2ZWcs50ZuZtZwTvZlZyznRm5m1nIuxthSKVFpcldihjfQ4tULoRN2jXfWp+00BUr3yeKWQ6JRv7SsiVo+b1gJNAbbqC5+3NqgKtHmLg+j0F14BVBVjq1A6oMFiaNYqoWpVoLzwWrU9UDZWNK9psBgb2fbqbgbNNgv6i8GQ9evPCq+VciwVZo9VuM6Lx26HsO34iN7MrOWc6M3MWs6J3sys5ZzozcxazsVYWwpfOroKwPrOgr++8jwAL139PgA7d/yA7+2sKorpLj9Ds1f1rR/Rf17D78nPnh08uzW60fS8z8bqE3KVFTtH9WpXcwHw6rXh9w5XMyO7AHrfxcvri6KngmmANFBwzUNKG1PQnDGbn9VaF16zi46POOt1sFhcXTC9t6MpVq90hi7FbkvKR/RmZi3nRG9m1nJO9GZmLedEb2bWck70ZmYt51k3thSeK04BYKcK/tbqTgAuOumvANh/yjm8cOrJABQb5bGJoJ4mEi9m00WqQ5ds1kvdr37ULJkkpjmtP59MM2pFdT969T8/nk1kH6eOP+tbPzihJ2jaMUS1oJStp1nhqJ7zVW/6qvVE1Sqi2BnEjtSaIm30xcJpZNn5iN7MrOWmSvSSrpL0VUmPSrppVkGZmdnsbDnRS+oCvwe8GXg1cJ2kV88qMDMzm41pjugvAx6NiMci4kXg48A1swnLzMxmZZoqytnAX2bPnwReO7iQpN3A7vT0qKSHptjmifIy4FuLDmICrYnzf144avSRdH/PrOMZpTX78kR54tgvL02cg6S+gvnSxjngldO8eZpEP2qewtB8gojYA+wBkLQvItan2OYJ4ThnazvEuR1iBMc5a9spzmneP81PN08C52bPzwGemiYYMzObvWkS/f8GLpJ0gaQdwNuAO2YTlpmZzcqWf7qJiA1Jbwf+B9AFbo2Ihzd5256tbu8Ec5yztR3i3A4xguOctR+KOBWxhdP0zMxs2/CZsWZmLedEb2bWcick0S9rqwRJ50q6W9JBSQ9LujGNv0fSNyUdSLerlyDWxyU9mOLZl8bOlHSnpEfS/RkLjvGV2T47IOl5Se9Yhv0p6VZJh/PzOMbtP5X+U/q+PiDp0gXH+e8lfSXFcpuk09P4+ZK+l+3X319wnGP/zpLelfbnVyX9/QXG+IksvsclHUjji9yX4/LQ7L6fETHXG2Wh9uvAK4AdwP3Aq+e93Qlj2wVcmh6fBnyNsp3De4B/vuj4BmJ9HHjZwNi/A25Kj28C3rfoOAf+7n8F/I1l2J/ATwOXAg9ttv+Aq4E/pjxX5HLg3gXH+TPASnr8vizO8/PllmB/jvw7p/+m7gd2AhekfNBdRIwDr78f+FdLsC/H5aGZfT9PxBH90rZKiIhDEXFfevwCcJDyjN/t4hpgb3q8F7h2gbEMegPw9Yj4xqIDAYiIe4BvDwyP23/XAH8QpS8Cp0vatag4I+JzEbGRnn6R8pyVhRqzP8e5Bvh4RByNiL8AHqXMC3N1rBhVnh77T4CPzTuOzRwjD83s+3kiEv2oVglLl0wlnQ9cAtybht6e/ll066J/EkkC+Jyk/SrbSgC8PCIOQfllAc5aWHTD3kb/f0TLtj9h/P5b5u/sr1AezVUukPRlSf9L0k8tKqjMqL/zMu7PnwKejohHsrGF78uBPDSz7+eJSPQTtUpYJEkvAT4NvCMingc+CFwIXAwcovwn3qJdERGXUnYL/XVJP73ogMZReQLdW4A/TEPLuD+PZSm/s5LeDWwAH0lDh4DzIuIS4DeAj0p66aLiY/zfeRn353X0H4gsfF+OyENjFx0xdsz9eSIS/VK3SpC0SrlzPxIRnwGIiKcjohcRBfAhTsA/MzcTEU+l+8PAbZQxPV39ky3dH15chH3eDNwXEU/Dcu7PZNz+W7rvrKTrgZ8FfjHSD7Xpp5Bn0uP9lL99/+iiYjzG33mp9qekFeAfAZ+oxha9L0flIWb4/TwRiX5pWyWk3+luAQ5GxAey8fz3rrcCC+24KelUSadVjymLcw9R7sfr02LXA7cvJsIhfUdLy7Y/M+P23x3AL6fZDZcD/7f6J/QiSLoKeCfwloj4bja+pvK6EEh6BXAR8Nhiojzm3/kO4G2Sdkq6gDLOL53o+DJvBL4SEU9WA4vcl+PyELP8fp6gqvLVlJXkrwPvPhHbnDCuv0v5T54HgAPpdjXwX4EH0/gdwK4Fx/kKylkL9wMPV/sQ+BHgLsp+vncBZy7BPj0FeAb4a9nYwvcn5f94DgE/oDwiumHc/qP8p/Hvpe/rg8D6guN8lPI32eo7+vtp2X+cvg/3A/cB/3DBcY79OwPvTvvzq8CbFxVjGv8w8GsDyy5yX47LQzP7froFgplZy/nMWDOzlnOiNzNrOSd6M7OWc6I3M2s5J3ozs5ZzojczazknejOzlvv/uhvSJ5xhN+8AAAAASUVORK5CYII=\n", "text/plain": [ "<Figure size 432x288 with 1 Axes>" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(np.mean(f.variables['diatoms'][:,:,200:400,100],0))" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "<matplotlib.colorbar.Colorbar at 0x7f3f425303c8>" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV0AAAD9CAYAAAAf46TtAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXu0XPdV5/nZv3PqXkmWrZctWZYsybKetmX5oUiyGwIhCRAIMY8A4Rkg3e41i9cMaxpCzwywWMys0KtnIDPMBNydxUroBgdCh4Q36RAIrMQxfiTx+xnHlt+WLUXSvbpV5/z2/PH7/c75nVOn7q0rlXTrXp3vWqWqOnVedVX1rX2++7v3FlWlRYsWLVqcH5iFPoEWLVq0uJDQkm6LFi1anEe0pNuiRYsW5xEt6bZo0aLFeURLui1atGhxHtGSbosWLVqcRwxFuiLy8yLyoIg8JCL/o1+2VkQ+LSJP+Ps1frmIyP8tIk+KyFdE5KZz+QZatGjRYjFhTtIVkeuAfwMcBPYD7xSRncD7gc+o6k7gM/45wDuAnf52O/Chc3DeLVq0aLEoMUykuxe4S1WnVDUD/hH4HuA24CN+nY8A3+0f3wZ8VB3uAlaLyMYRn3eLFi1aLEoMQ7oPAm8WkXUisgL4DuBKYIOqvgjg79f79TcBz0XbH/HLWrRo0eKCRzrXCqr6iIj8JvBp4CTwZSCbZRNp2k3fSiK34+QHLrroopv37Nkz1AkPg8fvfXpk+2qx8Nh18/aFPoU+jOtnbBz/VqPGvffe+5qqXnY2+/i2t1ykR1/PhzveV2b+VlW//WyOF2NO0gVQ1Q8DHwYQkf8DF72+LCIbVfVFLx+84lc/gouEAzYDLzTs8w7gDoADBw7oPffcc8Zvoo63m+8f2b5azB/Hf/wWeisEBGwKamDDBz9/xvv79D1/MsKzGw3O5WdMDl0PqqgxYARU4QtfHmrbcfxbjRoi8rWz3cfR13Pu/tstQ62bbHzi0rM9XoyhSFdE1qvqKyKyBfhe4BbgKuC9wAf8/Sf96p8CfkZE7gQOAceDDDH2EIMYATGgFgDNh/s1vBDwxP93CE0VsUAumJ5gThtQMD0QhWyFu6jZ/u/OnGQvRNg334gojmCtX6gKVhHV/kvFFmcFBWzxhz5ziMhu4GPRou3Ar6jqbw/aZijSBf5URNYBPeCnVfUNEfkA8Mci8j7gWSD89P8VTvd9EpgCfnJ+b2MBoRa1BjHWES9gJjtgBElTMAaZnIA0RdetIn/g0QU+4XOLp/7PW8CAJopdkSNpxq6fuHehT2tRI3vbAUxmwToaFQXJFdO1/rk6sgWw7rlKk2LX4mygKD09+4BKVR8DbgAQkQR4HvjEbNsMKy98Y8Oyo8BbG5Yr8NPD7HcsoRbNQZLwXMAaNMtABEQQEXrrVjD1nsO8vtcwsy5n5898cUFPe1R4+j/eQr7cwjILZOx8378s9CmNLdI9OwEcKRpABE0S99iYcjnuNYDO12eoB1iiir33Qbev2jHaCPfcYRSRbg1vBZ5S1Vnlj2Ej3QsOQVaQxLOvCGJM8XzimaNocimXnU45vcZw5H+9lWxF2BhMBpKDWHePOG0zvK4J5MsgW2nZ8XN3ndf31oQnPnIzIookpzHA1T9y/0Kf0lijTriappBUI1Kx1hFsQ89qe//DxeOWWM8/FCUfvpf4pSISJ53u8DmpOt4D/NFcO2tJdw7YbhcAs3y5I+KZGegK+cuvkDzzNZYDy4E189zvGz95C0kXRAUO78d0M8iV3trlJJ89v5fwz358H5d0psit4Yrveei8HnuxInv0CQCSfZHrRhVEsPe1f8PFADv8z91rqnpgthVEZAJ4F/DLc+2sJd0hYaenMZPLUHGXJMnFF4Mx5MePn9H+1vz+F6r7B9JtW0m+/HDzBucIT//RDaxePgXAhu985Lwe+3xjlI4Dc+M1gI90Iw3Wfrkl3MUABfLRXmO8A7hPVV+ea8UlSbqftqOzzcRfVDtzGjO5DAD1UU2ycmWh1wFgvIYgQn7s2PwO5qPqUeOJ3zuIJkqyPCNJLZesnKZjLCLKqmyatd/5+Dk57lJGkAfk4L4i4YUoZr8jY8lz8gcfG2pfo/y8thge84h0h8EPMYS0AEuUdM8l7MzpynOzYkX19RMnXBTM/KPh7IUXSXdsJ3tyxMb7VJFEMYmL0mOSXTnaI11w0LsfgMP7kV6OGoNgIW9V2nGHAr0RzYf0lbpvB/7tMOu3pHuWsFNTxWOzYgVmxYqKt1dUHfl2UkQMpImLhpPEGd9FIEnInn4GgOzJp0l37yB77MkzPqfn338rmjrvrCaw832tZ/ZcIBQxkHm7l7Wl1CBCcu0u8ofaq4hxhKIjkxdUdQpYN+z6LemOEDEBx0hWroRuDxUhP3pi7h2JkO7e4SqSEnHR8lea9VY5uA+ZcSSvy1L0i19h0wfOH8k++2u3kpyG5a/C2g8vfXLXb7wRyV3RQiBYUQXr7UeF1CAt4Y4zdOEuSFrSPQ/IT54EHPkmK/0FfaQD5yeqRBwy48NA735g4O91uv4yECF7+ZUBa5w5nvyDm0ifW8aGu3OW/9nS8CgPgn3zjYhVb/+r/bWtS6bp/ec3Adri7OAq0hYGLemeRwTyBUryDY+bqo5UK9vMF9krr1aepxsvB2PInu9rhTE0Tr7nMPmEsOd/eY7s2SNnvJ/Fgu473kTStWiIYK26ct1EUED+qfUzL04IeWNvrnOPlnQXCLORaUjEIVI+DvCXrtAfIc+JS1ZCbkm3bS16SxTwZc+FzhxcGOGyeaYL1rL6M0+ivcy5Ny4AmEyRnkV82a4aQQXM51qyXcxwibSWdFt4zEamycUXF1piPUKWZZOuL8TqS1zNflwhZSF7ZDQaY+zOqJ9r5Uci9KpQJXv1tTM+3vnuGtd9x5tcP4RcMT3r9NnlCemn71mg2KjFqOF8ui3pjiVG7aE8WwKJSS5ZubIS+erpGfITr8GLL53VMYY5hyCPxDKJOwktIuUzLRwZNWa+6yD4xjI4O20z/I+Zydx6NhHsioTJP7+7sdt/669d3LBtpNtivogTdEApR5xDwisibZH+KDech+pIj5+9/QDpp+ffb3nmuw665JcnUQRIBGtcHyNRd1/CPxEnISz1BOGFjDbSbXFWqCToLr4YrK1c5s9H++3TkGFoHTk/ebIvAj7TRGD329+EWH+J37XOqpWV2mo4JxUpuDLIKZq6Dl/JtPWNhsQ1G/J67OSf331G59Ri6UAR8uGGoY8cLekuMeRRRVxAiH6HQkicGQOJIX/9jfkd/+TJokpvkG951sO/+UYAJl87jeS5872GYpOo8ABj0MS12Sz8zOLIVXKnIWiSoAKdv23bU7box1jLCyLyPwH/GheVP4BrTL4RuBNYC9wH/JiqdkVkEvgocDNwFPhBVX1m9KfeYhAquu/q1eULnrRmu/RPL7vU2crOwttrp6YKT3I90k1WrQJr+yLmdMtmMAk8f6zoW1w/72HKoyW6X5g4psVigCJ0NZl7xXOAOUlXRDYBPwdco6rTIvLHuL6R3wH8lqreKSK/C7wP+JC/f0NVd4jIe4DfBH7wnL2DFgORrF7tLGB2brINCC6D4z9+C9mksO4/nVmVWaE3h6jb9yGuNwEq7GtJUhKtKmQhurVkX3uOFi1GCVccMd7yQgosF5EesAJ4EfgW4If96x8Bfg1Hurf5xwAfB35HREQvFGPnGCBZuwbyUJaqjnjniVUf/cLcKw2BofTkmHBF5oxo03VlmXt29OjZnN5ZYRwHoLaOiuExtok0VX1eRP4jbg7aNPB3wL3AMVUNo9iPAJv8403Ac37bTESO45pBnLlRs8W8EHTYZNUqR2aeeJO1a1zTHVMSHHE/2BgNFXJnIjmk27eBVbJnBkwwSUy/nDAHsqNHC+JN161bUOJtsTihKuQ6ppGuiKzBRa9XAceAP8E17K0jfHObvj1932oRuR24HWDLluFGIbc4A/gsf9HhbOVFZE999Yx2lV52adkVbVgvsAgkQnr1VRUXRPFadJ7DtrQcl0i3xeKGHddIF3gb8FVVfRVARP4bcCuwWkRSH+1uBkJB/xHgSuCIiKTAKuD1+k79jKE7AA4cONBKDyNEsmqVf2AQMQUxHf/xWzi9RpC3TnLo8mf56sH5uQuyV19zTXSsJb1io3M4RO0pgaJFZXr1VVWNNsgdfoDjGRP/+suKfbWE2+JM4RJpC2PeGuaozwKHfaPeadzEy3uAzwLvxjkY3gt80q//Kf/8C/71v2/13AWCVbJjJTGt+ugXWAXwQRiW8o7/2C2sefgE5rnqFJLshRdn3S4mVRfletLNaUyMpesvc9vVm/RsWF8+Cf1qk9aX0OLsMNaJNFX9ooh8HGcLy4D7cRHqXwJ3ishv+GUf9pt8GPgDEXkSF+G+51yceIt+1P258x4X1IBVf/AFLGAi8pu3tmsMpKl3JWSkmzeVfuCAmq6bbrzcPfDrnYv2lC0ubOTj7NNV1V8FfrW2+GngYMO6p4HxS+teSLB2+GKIIZG9/IojQmtJN6wfmgSTvbtcyeUjTwGQ7ry6sIQNTK4B2YsvlRJGnpOuv6wvCl4omOXLKz8QZ1IE0mJh0VaktTgr9JXudtJ5V5INg+zFlwoZIEa6Z2c5kvy63Y5ULeSPPE7uO5sle3e5frS59RqwId22tY94082b+vYffkDSDevdD4r3HZ8PTbc+A68PqmdVgddi4WDH1b3QYjwx/T2HWP6JakOWeffXPROkqRsntH1bEe3Fky5mm4CrIkiwq3l1Id2y2T+PZH+fnBtUFBGSacnaNWD1rGWUeqc0jaQPMWbWvhMxKZsVKxrtd3Z6ulzHT5OuvF4bdjoXzMRE43J7jqZJL0W4hjct6baYB5Z/4ouk6y9Dswx6zi59rjuMATAxAYmZl/sguXaXe2BxBJb7hJhPrs13AkWQGZyFTQvyLQivrhc3YRb5RZKksVS5CYOi24KMVZ0cMdupxK/b4XPO8yXrFiUUoTeuZcAtxhexxlnYxM4Syd5dhSQAnjAjDmtqhJ5u3+b6NTT4bAvCFYEEsgfd9unWK5039yxG/tQbo6eXXQrgpjEH8oqq8SqFIf6czpVO3EoN4w1VRlYcISKrgf8MXIcLon9KVQeWdLakuwSQrF5duawNBDxXxJvs29O/UITk+r1lz9wBU4gD0p1XF3JATNgVso39uh7Z154rpYURIZBwsnYNGJlT1043rC9dEjB8wUeLJQAZZXHEB4G/UdV3i8gErlXCQLSkuwSQHztWJV7vXkjWrkGSpC8iNDdeA/ire0PjWFQ7gGyTfXvKPrZ+U/vlh0mu3+t62F6/160Ya5uz2LTTzZvIjjw/zNsciNBJLWi7s5FtfCkfHBjp5k2g6twSzO1BbrH4oYwm0hWRS4A3Az8BoKpdYFZxvSXdJYJAOMmqVdFQSUXJ+9a19z+Mufk67wTo35eoVqJdoCDZyufUgljriFjVuRMigp0tqQZOy003XUF6xcYzIro62TZhkJ4aJ7dIUxet57kj3y2bL4hJxxc6RpRI2w68Cvy+iOzH9aX5eVU9NWiDlnSXEOpZ+Fmz+lkD2/qoV0UcgXqoSKU5rb3/4YG7NTdd69a576Ghzpk0cUS36QovtNUaltNfpQZVDTtZtaoipfT9HeaaXhEq3GRA2N9iyUGR+TQxv1RE4nlRd/g2BuA49CbgZ30h2QeB9wP/26CdtaS7lDDRcfdWZ72kL6JYygjWbRetFLZXdcpX7tdNZ48OhibbgE7HHcvYWYsl6siPH4/6BVd17bNyb4iZ13m0WJxwI9iHpr/XVPXAgNeOAEdUNfg3P44j3YFoSfc8Y5T9Tuv9XEUcIapmA1slFsULYRu3gZMaGjqA5Q88Wtne3HQt5qZr50+uHunWK11zHD8OKBRLnAnRFZLK2jXulJN0XqPe06uvqi7wDd/TrVcCzT0iYrS9axczZCT9dFX1JRF5TkR2q+pjuN40gy8FaUl3aSEisSaku3eguYIJCTd3J6o+4o0LFLR5P17nnbeMECDGReKds9fTgqYrBYkP57tMt211BFtrNRk6pLVY+lBGWpH2s8B/9c6Fp3HjzAaiJd0lgHTLZkdk3a5r5xgvF4Flk278jVVEtJAK4sSX1OUIX8hQWL9Cci5XYHb5YlYE/TQqB8Yq6fZtfaRXNF2Ptw0kGSrFwj788nTb1urx6g3b6z8k8+jj22JpYVSTI1T1S8Ag+aEPLekuAWTPHiHZtwd5+agjlaDtZrl7fnLK3S9f5saUW/qcBn0I5FTovJHg6/sqDEJoWF4n0XTPzvJJYW/z90lCumN7RUsuziNEpGFZSLbJgEjFSD/JmvJxXLbc4sKEqrS9F1qcHbrrL2Ly66cc0QZyCk3Gfb8EVJFuv4VM0/KyvB7xFtKDH3seiiUq2rBIYQ9Lrt1V+H9DNVtfFB0/T4zbfyJFo5y+84iboKtWGqan27ZCmrTRaot5wSXS2jLgFmeIN37qFi56OXcyQi+DLCv7G+QKtudWLCJA/wvvL8tLTbffHlbA27mKsuAaKSbX7S6fp1GEKdatGsg3TH144qnK7u2bb8R87n6S6/cOrIJLd+/wK7tjpzu2t2Tb4gwxxjPSWow3nn//ray/v8eyF0+WY8sjrRQoL8+t+nE5RGN2jJMcoGoZi6bz9sGUroZC863LAWGZj17nlDOA/C03w2fvHfh69tiTs25/oWPU04mXsjvDJdLGtIm5iOwGPhYt2g78CvBRv3wb8AzwA6r6hogIrhb5O4Ap4CdU9b7RnnaLgE0f+Dyn33WQ3uplTEx3nd4pxhOrX8nUCNTLBZomzVHtIDQQcP5QrTkOlGQbYICcWf3D5nP3Ay7ijZ/D0v7yt1g4jG1rR+89uwFARBLgeeATOAPwZ1T1AyLyfv/8l3CTgnf62yHgQ/6+xTnCsk/dDYDu3eXzsXlztt5Ht30SQhzh+sGRc41Ej2WAgmyh1GWbPs++kGG21osx2bZoca4wz4q0kWK+8sJbgadU9WsichvwzX75R4B/wJHubcBH/TDKu0RktYhsVNW2i8i5hgGdSJDQbqMe4TJAs42fDyLbWg8GALP/mqpsEB6b6jYAmoDkAAP8vy1anGcs1GDK+R71PcAf+ccbApH6+zC5cBMQl/Ic8csqEJHbReQeEbnn1VfHY/bVYodkFslspKVWS337CNc7Eiqe13qUO1vUa0ATQTuJ03hjm1k92WYMmpYj2/uqwVq0OI9QhZ41Q91GjaEjXV9t8S7gl+datWFZn5DnG0bcAXDgwIF2RPsIEBJNyXW7S/KlwbIVNN3a9IT6enFU27ePOJpVxdx4DWrKpJze/UCxarJ3l3M0GIOm6s7N0lgQ0aLF+YCTF8ZU043wDuA+VX3ZP385yAYishEI42GPAFdG220GXjj7U20xLPIHH8PsvwYzNeOXKBhBALVuggMh6o0iWW0iZ/cCSuRKqEe+0XrSc3qt+DJhRNDMFiQtJK4azrp+D03DKVucP4SZbRfi6J9RVaTNF/Oh+h+ilBYAPgW81z9+L/DJaPmPi8Nh4Hir555/2C/7nhtBavC+XYn1VygSa5oImhoXrab+lvjXgvxgDJok7rVO4m51aSLIFQGqrvCh2F7cPiYS6Ljf/FFPkGgxHOIhmU0DM5cygmVsmNuoMVSkKyIrgLcD/zZa/AHgj0XkfcCzQDAJ/hXOLvYkzjI2a/OHFucQRlxxRCDYLC/1W6NILigWOknpakjLqBVwvRpCpBtNjKg0jPFFGBK6lQXUXRFhfaNomoKxIJNILytGrzdNkRil/3Tc7Gej9tbOByG6DYRrJpdVIt6l/Hdn3OUFVZ0C1tWWHcW5GerrKvDTIzm7FmeF7LEnXU+GIDP4qjJRhTRBxSIaPLsuyWU96Ur4hQ+kHTsRDEhcrhtH0IGE666wYn33oCgtTgQ1KWItZLOTb4tzAztzukK8YUrzUh/pPsIZafNCW5G2xJE/8KhrNBN6MlgFUcicohW6jUnm3QWkTgIIdjMDYqt6rmRRgi5XR5hFAxvK16DqYIiq1aTu1Q0krO3khoVATLyIAbWYiQn32hIkX+deaHsvtDhH0E5S+mnTIC9EUkGIgK3F5OqtYEmZaKvsrCofBMKVvObVna3kt6FLWbpnp9N4uz0X8fohkai2U3rPE2JpIY54lyL5LqbiiBaLEPkDjzpfbOi1YEyVUOteXAtCDjaylQUrmDEuMu7lJdkGLdcOiHADBrRWTHdeXa4/0fFdxGaKXhLpZZfOayJEi7NHEfku4SuPVl5ocU6RPfVVwI9fbyrTrTW4CYUUzp0QrZb7CcJBVshstSrNRmQ7oAlOpTCiTwMObSkT6IhzXRhDuukKsudb5+F5RY1wzcTEkol2x7rhTYulBXv/wxx77y0sO2ZZ9uoMkttSZhAGVp9J7qrdpJuBtVGEO0BOiKULa/16QfA15bQHL3UU5Bu0Z7VV0m5x3mG73UJaCFhKxDvW7oUW4wv9xhtBld7KDhN/8y9DbTN1udBbmYBOkE5ZTM/2kaZYRWyNbFXLSreYbGNyjMm2WK/cf/bskcpxkr2uWY7E+7Q1AlcdWZS7kBatxYIiobaEpQVVIWtJt8WZQP7pfrrveBOaCFPfdxgV0BTyjrDqD75QrHf8x28hnXZEuuaxHLGK6SlqwHaMJz0wvRzJFdPNnXSQ5/1RLFT9uANmj81FuEAhc6gVJDEUHc9ziqq57NnWPgYjKmAYhkgHrBMi3KUS7bbyQoszxsRfuwh36nsPIQhkkObqnqvz1E4et2X0mrmI0vTU28WcXCC5dQmyyJUQW8GKyQ+DpIQYVsHmzUTbtG0qTs1NjdOLe1lRSRcq1ubc1xKGmVxWdo3zsNPTg9cN6wxR3luXEIptBxDrUiDcVtNtMRLYjhBfMYkRTAYms6XvFi8dZIrxXclMNy+rykJkGweyta5lAzuPGZm1UXkTitlq+/YUrgoxLgKXzCXRyN1+021bQS3Z156bY69LCzHhDiLaGPPto7AUSPRMMCrSFZFngBO467NMVWedDNyS7hJCIN1wS7oAisnEN6NRTK4k05kj3W7WR7QDh0jWdVbT8IEtWkf6pJgR0s2bhqouC+N/zI3XuCKNRCDx1W7dmMwT0quvKtwYSwVm+XJgAKk2/a1bnBXOgU/3Lao6lK+xJd0lBJs6srWJYDInJSRdxXQtKCQzFpNZR7q9HHp5P8kGxMtjwq0TwKCOY2HCsAjp9m3Q7Q1Fvvb+h90kiiQppgRrOuHkj25WrBcPqawPuVwMCCQbo4lwi/VaB8fI0fp0W5wxjv/YLYDTbk0GiSoIXPyHLpHWfcebSHqWdKrn9NupbjEBWGtf5kYSnoVo69sXq4QHMflu2zqUzhvPXQM3oUJTg51Y1j+anWiChbWlf9i3jhyXiHgQedqpqVm3s9PTFYJOVq4kP3ly5Od3oUEVstE1KFfg70REgd/zvcIHoiXdJQAThgBbRY0UZBsw8df/Aof3I6czR06JFF3DgMjeVfZiGCayKjuO1V5omircSSFVyJ08QJYNrc0WbSpxn+7k2l2uTeRE6s6hI2WX/FAlZ11SMNm7q/iBiUuPk7Vr/OkZsqNHhzqP+aBRLoibwpvhv/BhH2b5ctRakpUr+/YHkJ84cYZne2FiHvLCpSJyT/T8jhqx/itVfUFE1gOfFpFHVfVzg3bWku4ix7H33kLSVURh5Z13DV7xri/DTde6hmGY/qkQ6hvhhB4Kua1GuPUZaU0VbSEKNd4C1uRySJNivXTrle441pK9MHzL5RAJp3t2+r6/iZch3ElpalBJkE7iknFeu053Xl1IEZIkLkmXJKQbL3c7jm1wqmSvzG+MVJ9kIIJZsaJ55dB7eB6w09MkF19c2S4/fnxe+2jhME9N97XZkmOq+oK/f0VEPgEcBFrSXapY/ZEv8PUfuQUFTvzQYS7+o8HEaydcC0XtxAtxzW5yXJQYiMe7GaiRbOhMVkSz4fvfMM5HbY3Yw6rGIEYhTYHsjNs4Zo8+QXKtj2Rz1/tXEylbTVpfxpwkiFiXXNy7q9rhzGvC6fZtkGWFLAGQblhfXAVonpfadpAy4vc0YMJxJaIt/pZmXmRZEG09sj12rG/d8etbO77QESTSROQiwKjqCf/4W4Ffn22bYZuYrwb+M3Ad7grvp4DHgI8B24BngB9Q1TdERIAP4hqZTwE/oar3nckbajEc1LggFeDkew4PjngF7ETifLpCMc9MbFJ6dIMmGiLX4uY/oIkp5IcK+Q5Ak0bsJA5PRml6ViN7JFfXF9gYbJM+Xe/1a9QRsQXJ8iLyHXZWW7rx8pJw4/eWe41ngMtDxEfhAwoPklWrohOP3kds9TtHUsiFihEl0jYAn3C0Rwr8oar+zWwbDBvpfhD4G1V9tx9QuQL498BnVPUDIvJ+4P24EezvAHb62yHgQ/6+xTlCqDw7+Z7DoDD1fYdZ8af9xFtcfneEfDIh/e+lTKUUrWZ8NJiXG1Y6kJUSQtNHNsgPTWOBBjXAOWtLVGPFHEXns8ZSZT/JOLl2lx+UWUa4syXfQpvJdMN6SJKh2k6mG/ygbGOKv1m6YX1VzgB3DgOSlvOVOlrMDmc7P3vSVdWngf3z2WZO0hWRS4A3Az/hD9IFuiJyG/DNfrWPAP+AI93bgI/6CRJ3icjqMMByPifWYnh8/UdvcaW7Gb5fgnL6XQdZ9qm7K+uFAgnzj/fN+h8fJ5zSHdujy+Iy+RYSbsVY97BOrs1abnwfYMR73Czp1ivPqOghSAyVYxUOBqod0DziKchxkm5ex335lVlfT6++ijBxI3vqq06/rhC/qTb14cKuuDv/EPJzMF59GAwT6W4HXgV+X0T2A/cCPw9sCETqJwL7n3M2AfG354hf1pLuOYLpKZ1TFlFfGCGOfLvf/ibEKp2/cxGt+cf7sN98E/lbbwYLyWfv7d/X/mvAOL9sudAnqaLkj4pzDFQ0YCidEQWZ1IjW1pZG+zxjmcEC4lpMany8kNCLI2vOnGibkO7Y7pJutei473n0gzJf90aLc4NRaLpngmFINwVuAn5WVb8oIh/ESQmD0HjV2beSyO3A7QBbtmwZ4jRaDMLKjzVruNPfcwixoN950C0Jo03EAAAgAElEQVTInctBcoUEsrdHCdkQhIVI9c03AtCzWvRlqPRkCFfGgdD6olhT1YRrKKxpxuvCvhAi3b4NGF5jhXBOgnbcvZu7Ztz7FEUxlYg3uW538R7qEywGoRLxQ6kFP/n00OdZbOsJOd28iVF2UGsxPMa998IR4IiqftE//ziOdF8OsoGIbAReida/Mtp+M9D3qfI+tzsADhw4MHyxfouhsfwT7r/sxA8dJulCetqiViGOREVcZKz4ZeI+kYEUFRp+M6uI/b4BdR01LjOO+zg09czVss8CqnNedmePPUm6e4fTZqEYNeQSheJcDaG3RGGPE0DdmCB/zOyxJ/v2ne68ujzOiCvfsiPPl3a1FucXzbHAecGcpKuqL4nIcyKyW1Ufw00Aftjf3osbxf5e4JN+k08BPyMid+ISaMdbPXdhEWxkp991EDECmbphkzg5gmBO8ESEFQQtubap1DfBWaliP27seIAq0eaRDFHpXtbwya9pnemWzUMRb7J3l7vMMqYs/vC6s3+31fPER9zBIhYRLFAktZrIGFp71mLHuJcB/yzwX71z4WngJ3Hpkz8WkfcBzwKhO/Rf4exiT+IsYz850jNucfYooloK8nXLfevHhlLbwm7ljbnlOrWItU64UUku4FtDDvCZ1Yk4yuQP0zgnf+Rxl1RTRfLS5Vaxjg2SPOLIO0yzYDDhtljc0DFPpKGqXwKaKjLe2rCuAj99lufV4hxg2afu5vRthzB1q5Kqlxcoigpc711HwuCdD+ou1QFXjACOiG3paqi4BUyk/frjAI0JpHTrlX2FGBUYGaqvbv7Q4y4ZmFRta6E02K2k1fOsSx01wk2vvgp6vdZdsMQwtvJCi6WFZZ/8IjPvPOhI1mqh36KOaGP5VgXnxw2JroopIUgJUtFq1U98KFwNwVKmJaE2uRTqRBwItg8ic8oN9ssPk1y/t0zWRcm/gUM0Kztw66Q7trtzO4NkWYvxxzi7F1osMUz+hfPvdt/xJsBHtT7SLSD+HwH1+m0g36JRuvXEa3CjfcBpvVBUfbkn3sgVeuLK3CFGINU6waZbNhfEWxQONFmv8rywilWGaA6D2DvbEu6ShPvtXRjSXRhRo8VYYOKv/8UNswyRrrqJEmJ90YM4WSHcbOpuRYFEXBQREGvAiaCdBE1Tf5+48t/UdRpLt22d8xwrhBskCDFuuRh3a0D+0OOujWVmfVOdwfa1AjUHRUu4SxtWZajbqNFGui3o/G11irC56VpHkkawk4nP8PsXI9KSXL0kUdNIQzKquCmhQYRa3GV+Chjj9NJ4v3k+sGggLA9e3jDGvUkGKKxgAYOSg/VzricEWyxZtJpui7GBve+h4rF8442ONOvNZCxOjqgnpQagsG0ZA1hXsCA1IlR1rRa3b+tzMMTFEtnTzxRE607GrVNMk4h9w00z25qINba7+SnGc5Umj2qc+zhaz5b6qHpFsOPsXmgxnhjll3XQl0z+6f5+N+OtNzhpwXccU8W1hoRSekh8pZmfmSaiqLVuPWP6eteWB4wcDBHRx3PRCsId0I2r6AVsaxHtgCg2e/oZR9hZXkwxDtudaU+IFuOPhbqWaUm3xfzx+S8VRJy95WbSqR6240tt8+ijHBGnCogxKJ54E/E5NSk7gcWoN8kR6S9eKI7ToLvVbWABlWo4XxSxbSv0siKaLpur5/Nqrt5iEWEBE2kt6bY4KySfvde5zg7sA3CThcE5G3yJ8awIlrLQMH02KSCObBvbTUr13q0YuSaiKLruZKjs2xQVd+kVG6tuhiFaObZYJGg13RaLGXrPAwDI/msK4tUkgU7iSY7CKxtf6hdab2KqRNgkOVQO2PB6vR9tWJ5INcIFf06hvaJrlFMk6MJI+ggt2S49tJFuiyWBettEc+M1KEnRnayvSxl4e5qLePuan88Hg+a4eX9xkbgLBCy+K5lVR7Sx/SycE0CeF6N72mbiSwMKWNuSboslCHv/w056aGhuXjRBj6BxdDrbKKCmxHODDuz2hauWC43Ni+1NVKxhokq1MErIR8IhCre2nAJRc2xkr742y8m2GDsozK19nRu0pNvinEPveaAcBRSmPITBlXE0GmNYwo38wH1TLGLrmCqS+6Y9hYuidpygKSfxl9FUo+4gPQSd2Lppxulll/rz9rKJT9Llr78xyxtxWMr2rMrsNyhdK9HfVKMhn5L4kkaRcz5SvvXptlhQnA/7GZTj0wPS3Tv6ZYFBGIZsa0Tb2HNBBE3oJ9442aY1KaKpM1ohi5hom+o3OVm7ZijiXYpIVq/uW1Yn0mTlSsQY8pMnz9dplWhJt8WFiKKT147tiK8w69NkG2SDgmjj5RHJAoNbOdbbVxLtp0LWnmyb2kAOGLle8QczXKS71JCuW0fj1OOG3hcLQrYASJtIa3FhIFjLwBGfvfdBoCzhTffsLLXeoriCShTbFxGHS9O8RrCFRjsgQVf38g4T3dYr3JpcFIl4/7GSrF1T2TY/dqx/nyNGsnLl3CtVCktM82vxD0iDF1oa+l6o2j7CHdsfnnGOdEXkGeAEkAOZqh4QkbXAx4BtwDPAD6jqG+IGwH8Q18h8CvgJVb1v9KfeYjEiWMvAf+YPu+nVkrnBktZaR5553rwDygSc5LORZrlu5bUmou17rV93rGw3TNMcdfY18Uk5VVcll6xdc87JV2tRuJ2aGmq75OKLyyfGVApd8tePlXKBJ+DGaBbO24/LWUF9TmFEEJEEuAd4XlXfOdu684l036KqcYr2/cBnVPUDIvJ+//yXgHcAO/3tEPAhf9+iRT/u+jLgk8nR4mJ45ID5a5WvyyDpICbTxgi4IZptIlSrzSQN/TJDfXtPUI3ky2ijQLNihTuWj1znfeke/t4hmQWN5Bki3IJ0ox+k/PjxeZ71QmKk8sLPA48Al8y14tnIC7cB3+wffwT4Bxzp3gZ81E+QuEtEVocBlmdxrBYXGPIHHyseJ9ftdg/qpBlrvYOiz5hw67LBXOnrgWQ9m7WC6rkVb8KTb01HDK6HUVjOAtmeVeY/vEcjiBjSdev6pQVPsmMrGwyLEckLIrIZ+E7gfwd+Ya71hyVdBf5OnLv89/wk3w2BSP1EYG9gZBMQdwg54pe1pNvijBAIONm3ZzDh1Ym3buuCwRVvAfVET4O9CSgb9sRaaDwpYzYMeD32/w5bgJGsXBnprxHh+ui0Ihf4fQfpwU5NYZYvd4+np/0OE1dEkiSOaOtar/87ZEeXiCd5eNK9VETuiZ7f4Tkw4LeBXwRqf/BmDEu6/0pVX/DE+mkReXSWdZs+VX1vT0RuB24HWMaKkXsVx7FdXouzQ/5A9WOX7NtTXUHV93GInkNEdA3yQhzFhcf1Xg1NRBkTkmrlknzgNoNQkywKAh5mu/h4/nl62aWVHxDVsom7jSSHgmzDLiY67n1FXtmqTm7JXn5l+Pc1zlCYR3HEa6raNCMSEXkn8Iqq3isi3zzMzoYdTPmCv39FRD4BHAReDrKBiGwEwv/GEeDKaPPNwAsN+7wDuAPgElm7QHnEFucC59Lza/ZfUz4ZIDdoAuAiNVF1LRvBVZapUswUCp3Z61l4sVVingtJQ3lcQdjBwzuHJDGsU7/JPdEYtZfHEzUFgabrLyu9xfGPSjQBuSBbqxXiXTKE6zGi4oh/BbxLRL4DWAZcIiL/RVV/dNAGc5KuiFwEGFU94R9/K/DrwKeA9wIf8Pef9Jt8CvgZEbkTl0A73uq5LUaFuLdDQcBhAGVt3A7ie/0Wg9tqRBrLB3HUa6OKiaQWDcUE3dRScmB0nMz+LW8i+UEtK4N3uEaK8X5mI8gikk4TwvijyrGsAraw62XPPj/4vBczRuBeUNVfBn4ZwEe6//NshAvDRbobgE84Jxgp8Ieq+jci8i/AH4vI+4BngRCS/BXOLvYkzjL2k/N/Ky1aDIa5+Tr3INJRK34GVUDcfDRA09C4xvV/0NRHwMGHam3Z7jGOSOsRsJF+Yi/WHdD3IV53Nsmh2E4HLJdSPgmTQVUrEzWGQbp9W7Vyrw5V35XN/W3qU5uXEoaYj3pOMCfpqurTwP6G5UeBtzYsV+CnR3J2LVo0oBit7nXVUJ2mXk7QL36lb5vkut1lv4eAeLIEVNs9QmWCBOKb3iRSEuCgGWtNhDYfjbdpu6SqIYdKvmGR7ry62EcoLqlY62Id+wzIfNGh7lEcxS5V/wHn4poVbUVai0UFc9O10MsRa92Y90TcFApj4K4vD/we5Q8+VjbbKXYm/cQbojzoT7QFNEW4Tc/rkXF92ybUB4DWyDy0wUz2uvcSlzNnjz4BlEM543LqxqPmnBGBLw1E/8/nGS3ptlhUiIdmBsih612Ee8t+NBHkn7/UuG3+0OMu4s01ipZlcEezoOcOKvstTqCWkIKCcPuiylnfXMM6hcdYvJ2rGq3GPSiSa3d5qWU4ZI8+PvdKSxnjKi+0aDHuiOWESuxy6w2OyHzVG1SLLorpwdBcBBHQlNyaS8P1+xlItoP6PsTnMwCVnsNNxx6G4FVJ9u4if+QCJt4halzOBcaOdCXtoFlvoU+jxWLGrTe4e08+oclO3PcBqFxWFxOGmwgrlgnqSah61n9A9VaBOsnORpDx/ga1sZxLK67vP34fJio4sSBZTvbEU7Pvb6lgfj7dkWLsSBec6G8vWc7py1fQvdiw8s67FvqUWiwS6Dfe6B8oYsWFvolArk4PBmSmh0x3IcvcuiLOy1snzKjKq1HbrSMurvDnUJ7YEFF0fV8mItbElBp2OKe4DSYNEkYToTcl+kTAKCSmHEUPxcj7pYqFci80DT1ZUIQSRNtJkEzpnLK8+tO3LvRptVgkkMx1KRPrnFVq3OW43vMA0sudnjvZIXv6GbJnj5A9e8RtGKq2oCweSKJbHFnWbwFWZ+9iVmmUo/2EG4g2kK0xxXloahoj3Eai9dMsUCV/8LHihqrrHJbZ6i3X8lLblqScbtt65v8RiwE65G3EGL9IVy3M9DDTPeTiDqYLk8eVF37xVq74D59f6LNrMeYQq87GKuIKzoKuenBfpW2DHNyH3u3lhlAcEEWViKBxYUTT4Ih6mXElhJEKgRUI28T77tvP3FJC45SNoqmPO7fgZgiIp3YEh0NhF4vOI/Q2bnFuMH6kKwayDHO6i5mZxE4mpFOWi14au6C8xZhCrDq5LpgPgjyQ0By5JJGMUL9sj7249f4v1JwPdfKtN2Gf88Sl73HfFI0Y9Z7B4GUVF7WHkfJNnts6IV+IGNviiPOKWjmi6Vm0Y0hnFDUL9BdqMbYwN13bT2aZdXYwtIgmVcSRbSBidURVfKIadE5R9YUUNcKtl93WZ7PVyXIuh0FfUx6qz+vrRURb2N6K11xCrKiyC0M0W/RDGUkZ8JlgvEg3wJvWJbdIz8KEIZ2xvP6vb+XEFtj6K63McCEi3bwJEkP2Ndc51N73UFkSTENU6CcfNDlXY9LNnn6miArddo7ZpOa1LRB3GKv7ZeuRclIj4YbzqBNqsbyOWoP2SkVZJCsUyHKwtphWodOnQS222208lwsObaTrEScJcsVk6sayGGHFqzmTx4f7dRrnsdZt28kzhG89mO7Z6fRWY6Cbuay+wWm4BuhVG+Mk1+9t3F1y/V7yrzzinlj1fQ3iZJc7hiSmHBdfJ9XAv8b0j4DHk61Ul8XHUKLPu0QEHLenjKcZx1F5pYzXWb5Q9WSbkz3f19yvRYRWXsBFFmKMkxmCzqaKyZXcKKSCyZSp7z1EPilc/EetlexCQtAmg/WrD76ZmP1ytWot/8ojzo/acEmf7t4Bp6Zdn4Xg2wyNbhIpiE7AEa+U28b3OvA5NEoHxTnQGHEV0kFcFuzdBzJoQnFIhPkmNemG9ZDlaLcLqtjpaXTQwM0LES3pglrX2V4SgyZl9lhyxYiiHS2SIkkXvv6jt3DJf/nCAp91iwVBvcQ2ukyPCTZ/8LFKAQCGcmpw1Gc3+9pzpFdf1Vj8oEFOCFGtv7SP3Q3FBGMfiRbP/fr9mi6VbQsEzdZbviTLXLe0LHeyR5ZFBBwN7zSG7MWXKodYav1vR46WdD1UYaaLiCBpCh2gYzA9r7Olgp00qIF0Gk69+zAmV1BY/mdfXNhzb3FeYO97CHPTtWXJb13T9LdgkcofeNQ1iCku/avRatCIK1ptvVlNU5/yXF3DdE+2RTRsIuIlEHBtW0+0Gn48wushUvYHlCRBrEW6nnx7zt1TEHBu+yZAtJgboq28UMIqzMyAtZg0QfMcEsF2jPse+E+vpoLtAALqHQ/T330IkyuTf373Ar6BFucFqkhwKsxRUlt0F7OUl+YNlWDZE0+5FoiRgyZuKNNUduvGwIcrMEXE28hqckNI5hURcKEPU1mvsLmFKNoIag1ijCdfg2RpWUGXWyTPIc8xExPubbaJsuEw7u6F+lx3EbkKuBNYC9wH/JiqdkVkEvgocDNwFPhBVX1mPielee603dNdpJNijEGyhOyijje9K2T+uyCQdwARbOpu+fcdRqySdJWJv2wJeCnC3v8w6Z6dLrosCKrqKkiu3eWINlfIrSPVq6+qrIM6P2vFyxr1PCisY3GSzUensSwgQb6oJdEqfdAjOUQxLiJW8YFDdR01brmdcAk68Q3YTbfjvbiRfOEb+kjacfdJ4rTb4kRaDbcJi6EMOMx1D/hN4LdUdSfwBvA+v/x9wBuqugP4Lb/e/KFe8A8ljZnFdC0ms8WlgfgyveBwEOtuADYRessNU993mJnvOkj3Ow+e0Wm0GGPk1kV63dyV//Zydxk+E926Pej2ikYu2VNfLXouZE8+XfQXCCWvlYYvIYmm3kFjo2Vec+272fIzK9aXJIdtc+3bp5NDKMa1aXTlpiLY1GA7QrYyoXtJyszaCbqr3A0j7jzC6Wa9olmUGGnJdi7okLcRY6hItz7XXdzsnm8Bftiv8hHg14APAbf5xwAfB35HRMRPlJj9OL7uXMRUIobiQ2ot5FJWHPlrMbE4svUfWDVayRrnE+63Jf/uQ56YYfIv2gj4XOF82fUCQaY7r67ODgvygc3L3grxdl7DTbdsrnzO0q1XVvYLvv2jBr8vEEwOoVgnLtFt6oMQtg1Ws+KS1rqrtsy6JuwYBHERryqaOBuciovk845LLKspCXr6sg62A6trH2WNE2wt8TZjEWi69bnu64BjqurbNHEE2OQfbwKeA1DVTESO+/VfG+pIofY9SZx1J03QidT5IIOrIXxRrCJGyiSFUjY7iS/ZgFCN5FwRTv8tlgP5hGA75YfaRSEgGQVRayKs+NPWpjZuyJ54qmzNCBW9Nt28iexI82DF7NkjruACiHWAdNtWMOKa4vj2j07rpUyQBckhRhPxBkTWM7de+Ny6z5YY54YwJL53REjC4a7mclBVbOqJOWi/Aq/83K3YBC7/rfEtGhqlN12GKakeBuNKugPmuje9ax3itXi/twO3AyxjhV9ovE9Xyi5PYQdJ+KUPGpigqTRkhaP7WNdqKMMM66rAxEnXKKVNwi1OZE8+7Yi3sGn5CG+OlozZkeddxFvvaBN6FxjjZIgQVe/eUSNerRJtvXChCcFWFtZXcVbJYs6mOGUtFUzmHBEqbuaOCIh1SWR3E+z4pcMXBWSBLgKG+e/qm+uOi3xXi0jqo93NQCh/OQJcCRwRkRRYBbxe36mq3gHcAXCJrPX+GeuSaKGyxnd7Kqw5xcb+rsF4rkJ5SRd/7itNQdxr6r3uEi3vfesBTE9JPnvvEH+aFmMHa+c9wTZ79kghLTTtzw21NORfeaSMfHfv8HpszekwbPgU5IlAwJ58IXef6VxJcnGqSSglNuIiXuN6BavfRKyQLZvlUJOzvDgk7Mxpt6+JidYdcZYYZhpw01z3HxGRPwHejXMwvBf4pN/kU/75F/zrfz+Mngtei+p2XRI4SZDcwkQHAUwvBybRTmBJ4yLeNEo8hGZR1l/mVQzo/kEU/YpSMq6PIjQRciPYbz2AKCTT2cCZWy0WFpV+CeAId4gptmbFCsQY8pMni2XZ154riVdtKTeI0Fu7AvO5+yv7KMg3iq7dx7JqFRuIPLKjxf0U3AJIFKxgVF1xhndUaObkNWNAbEIyo+STCRt+ewhpYa5G7IMaquOJ2189LBlr2gjkBRFZBnwOmMTx6cdV9Vdn2+ZsLkx+CbhTRH4DuB/4sF/+YeAPRORJXIT7nvnsNBAvacddcnW7SJ4jaUpirdN3dRLpKKIGmwvaMT4BQdUfKVoK5kXPUP+y1SLaLRBVFDlbEOTLU3jbAUzPugKNz7cEPC6ICTbdvs1JAluvLIsdBsBOTQFg9l/jPlun/XiobrcmTySgysSzr5M174qisq0gNH9J1qT31hcVncwoCiqw1t+7qFnEOSIMONJNpLjCC43aL/1vg4uCklWryuOlqTvPyclSusute9+5RbMMrK0l4tx5hb/ZksHoEmkzwLeo6kkR6QD/LCJ/raoDkz/zIt14rruqPg30+bBU9TRwVulrzXM0z0mu3YW8dgydmkaWL4Pp00gvJVF1l3srJjCpkC9LsalA6i/vKkRK9RctPBdBfRKjeElt4yWiqKJGsBMGvumm8EZBBPOP953NW20xC8z+awDQTtI33yxGesVGmOkWDV6GIV4Ac+q0I0wtxFQX4ap1kW+YnKDW6b5ZXim9zV5+heypr7rjB4/wMp+TqM9Lq1/sxSXMULaR9JFv0HvVWic9ZLjXrHGbRn5e+003uY58XgcGMJ+7n+T6vYWfWMUTdmrIl6X+2O5Y5nTufnx81ZvM9NzfJC97PCw50oWRRLr+Kj5cMnX8bdY9j7UEb1cuw8gaZPoimD5dRKDSy1GjJFNdN3I7U9STbxEJgCPY4G5IxFcP1RB4OpSRhg8+FGWldSM8hCB6gdKfFwiKTmGH97tbIJvPf4l029ZCu81eeLGy3TCEC+WEhIp1LJANUeOYrVe619OkuBrKXnjRka3/TAaHRFF80US8YXkDquQLIbEnJD73kIMaJM+dxmtx0kNS3Z/JfGLv8H6y1Cee/TnY1JRJ6NjNY3xuw3uJkxnrSDzOnVy9Ae0Y952wujSCjRF9fX3h2L3ADuD/VdVZ+xGMNenqF79CXl/2jTeSfH3GkWHq69Mz68qDJa80ylEj2AQkqX/wQ3JCUG8NIwcSSiO7cTcfc6DWV8JpsJK5CCR7+wHshJlX5dsofawXRJvIMEL91huQ3GKu3+uuhEY0Qjx4eQtN10hpJYOChPu2q5E9uOKLdNvWkqAD5mtzKg6Zly4H8dJD5qNka6vDKuuBgZdc1TiylcwWshkSCNkVErnnoKnQm0zRYjnlrDnvHJJcSb73EJPHM5KpDL7wZRYbhHm5Fy4VkXui53d4IwAAqpoDN4jIauATInKdqj44aGdjTbpNkH+6n+fefysTJ2D9PWUiBF/hUzgavMZrcnWf3/gzrxQNrsH/kvtfe02kEt0WnaxwEoNYH+Em7rGZyZFMmf6eQ5y6POHSD42vV3LR4/Nfwly/t9TdjZJct9sNXZwDM+88CDKHJbBInuG+kbUeDU1+33TLZjBJGRVv2Qy9XqmZWvqi3UENzRsRAoBQaFE4HZQiRLU1co/KldHyyk18BBNIuJA4Uke6dsKgubocief4UJxhU/d3CX8NTYR8wmBmpNEjOvaYn6b7mqoemHOXqsdE5B+AbweWDukCbPqAI7aXf/ZWVj+dkZ7oucsqoqSZul9tV2LpPqzq9S58scNsXsq+zv2mzDZL7n/1c3XJtUxZ9mpOeipl6nsd+R7bl7Pzf2i7np0LhFaL4f8o2beH/IFHZ91m8i/uRr/hBvQbbkD++UvIgX2FTpzuvLokNeiPbFUHNwRP3Veo0H+tZ7Ysc7rwju2OOJvaQMLgQnw74HGoTTW4UmMoSbfW9yFQpKi3oanvGeF3WP4dY0lOyqIjVfzgI/d9EeXiO11+KFm1CrloBXRS2Hh5X1vJRYHRuBcuA3qecJcDb2OO1geLknQD1v8/jnyP/ptbufiFjPSkzzGr84IZ34UqJM0KrSIODJqG+8XJtWhGVpGQSL13OBUnXfieD52TGZ2TGStegjWPp5x+10GOb085uUW5+hfavr+jgHZCp3L3Y0dmEdVymgQUOmu5kf8/POZbIF63G6ZnwDfEKXrtxsQLRWPwWScwhGOFaDQxPp1S+2rF0ScR8c5xiTuwg1p9O4mmB0NRduw++2X0K1BU01V6AFMGKOEHYdmnBl8V5MePk1hLfuLE7G9gnDEaTXcj8BGv6xrgj1X1L2bbYFGTbsC6/+TId/q7D5FOu6Y48bwoydQlpYNFrMFWVvZiHXAQC8VgwwDjtV7j4w9D0XAnmclJTmdcdqzHZfdB/pabmV7f4dUbDFf9citBDAP75hsLf6x+ww3OM+2b2pDlRRPyQcUQ6e4dZI89WVSRFehlpfYKVW9tg3sl3Xi58wD7puCnvv8wF/3JXa5Hr5+8Gy796+PL4+dFKbE/VvxD35esDdLGcBZ3t794rNAQqIwS8v1NFJe7kFzJ3n4AM5P3eZQDFjXhMi95YSBU9SvAjfPZZkmQbsCrNyZc/DXDyuczR34hGogjVwvgathFvUzgk2MuWWYLK43bIL5kq+VDatGGRD7h4JQIyTcy5aLnc1Y+B9nbDnD8qgl6F8Hlw5jaL1CIghzchzk1g7x2Crq9amObyNLUhFDA0LeeiHMZFI1xtLCJgU+oxZ3DwM1m27Ce7OVXuOhP7oJb9pNHCaQgUYT9No49jxvp7Lza77eMfIGSaOsRbly4UC/AiEqPy0g2jtilOrutIH4qpfXFn8cHJ2oEe1GKebvzqacnushMbygNfVFggcxHS4p0t/yqI7CZdx5EQr9dpfCrA1HVmrfohImvwS4W/LtFxqAk35iM6z0fKqhHLH47k7tLt86JLmsfyciXOQlCjbRTLxqQvjHlZoJNzzitNJBnaKHI7CNpzE3XItM9p0mqkj36RNXSVRBazSPTi0ohinXc8cYKEdIAACAASURBVNL1l7ltn3oB1l9G9sqr7jxCbwZP2OkVG4sR6GGdGMX6NZIuq+JqjFAQrCkJuDYEs2/9MCwzuBOkJN/Ctxt6mYTXcd+Jzt/+S98ul1S/smCRWwAsKdINCG0bs7cdIDmdYTsR6xa6X5U4NXyAs1JxgMGXIHE1W+xnDNGuq48PK1PssEj+TGckUz06xwyaGvQbbqC7ZoKTV6Ss+73ho99xnnp8Nkg3XQFvnOj/AYOhkzYy03PtQHtZ2cAm7Ke4fPf/SVbdMeMo09a+lcE1EZ1TQcIBp2cq54oq6WWXlp7dKGJVtXDs624/69aV2ydJtSG7W9tHuFGJcoymWW51+A59mgaLmJSNo+JufbjAZcm3P72QI91dN2/n0/eMzm8aiCj9785aFyg3f+vNANHsKsBHt6GFnqaCaqmvaaQN41cv4Mk0RMb16NhpvP5LKni9zJYhg3g5o2uhC8tPdVn+guHUdx/i9BrDmt+/8JJv6RUbqwtCaWykqQ7cNlyyB5zuVck1CxFi4ggsDHmsk2uoOgsk/2rZlTTdeHn/ga0tSbK+r0DIYbJwJA1ImfUqty/uaxGsl0GyI0dKLbqGchSQq4rTpEymlStRFExYbw2rywvh2FPfe4gVs5QYL3aMez/dJYHkM2XnsOztBwg+3qL3bhGpUhRFkPjHQTYLubfw3Qnuh2gKQNUD2B8dYaBS9RH0X1XIcy569iQrXjDIjddgl3Xorp5g4q/7L/cWMypdvWp+WCg9semG9QMv0RvhSTR7+pnqpbpIrQjCE288YmfAJX162aUF8WYvvuS8uDNRs5fKQMvocUzGg9BHitJMuB5OvoiiXV8S7yyRtftw7NCD2mu4hbwQCDdc/YXj5IoIJF04fdshjl2VLM3cQ0u65xfpp++h921vctGnb3xeEKspv4RxUiF+vWiQHgg2VLaFijVbFmuI4o3prmOUWNufrfbHA1yGvudY2Uz1WH6qi+6/piyLXQIYulR3HmPEi2bmnsT6EmMBSVJGsybSmRqSp2FZetmlbp+vvuZK0gNBpXHP0RrqhBvtN1SzFZJGYT2LkmNFEi8mZeMTiZbs6eew33RTUZKb7NvjiDbWcqFs/h/mrnXKCFsUCJO2w9v3tjFrhO4KWbqE25Lu+UdTsqD3bW8qHgedK4wDKkqMhaoGJoBKjXQBq25+mydgN9nCgk0qWWoRd+9a+PldBgK2tjhEsm9PeaKqaCfB3r94iHjYRjRnvP89OwsbWbpju/v7ZpknqqhRTUhabdlc3UE9KdVAwoF8y53lqPaQTqeIOBvhm/LrxSsq2f+iSU/o/1AQsKkUVBSkHJ/L7h3w8gl07y4XzYYI1yfMArlqYqruhCy6Kovea5iQYlMhu6TDKzelXPnrS5BwcX+OVl4YEwQiDoMsXRkkSIYjzUz95ZpWvgBqwme4XCYqTuJTxfTchzzp+cY70RBNCdV0mUX8JAJVJ1kUX8S6cR8XESfX7fZRtsvOjzPOFeGevu0QE2900Tem/Rwz64oUepmLaq2CSVw/hfhvmNvGqLbyOLZsxRGnf5y9+pqTQJJkMOGG/7vcIm+ccNFtQNIgP6i6Hwkb9YWIOp4V0sJkp+JOwNsh9T7fke0bbig+P2Ld56ue0A3l82UewmBImPyLu7lyVov/4kdLumOG0MDm9G1ullroXmYyF8manptmoYkvkyRyQ/gvnzXh+yHY1G2XLzOuS1PuIw4bRcLWv9bz1ihPwO7LEjLuDSdrQMUM3YdgKeHo7bdy0ct54YfWxCeo6g6F+IerTrSxWyHSj5/7lVtZ+7Bl+Sszrsy81ks5vexSF/l6LX5gpNukGYdlVgfOcIvRVACS7tlZ3ZcBTRLk0PVOvvLVdsZGxULR50iKycV+24kU082Q10/A1ivh9Aza66GnZ0AVOz0953kuKowr6Q7qjC4iV+GmRqwF7gN+TFW7IjIJfBS4GTgK/KCqPnOOzv+cY9knXfZ2+nsc+eYTgskUYylbRfoOTAXxqvZNsyg0YYUQDTvZomzHV0bDbnS36bpqJ0nd81DyGsihWvyhYCC5dlfxJdROUiRUtGMWZTeo2XDyBw+z9pFpklNd98MUkplx7wQjLmkZV3fViTcs80ivvgq9aBlX/vrnmXnnwYEVWSG5FpJ9hU48IOKdjz5dbPO2A3RXpUyvNWgCl/5uebkfrmySvbt8TxABcpSkKF8vPi/4SBdKog3wMob0cpfUtQq9DPWkuyQJF8aXdBnQGR34BeC3VPVOEfld4H24EezvA95Q1R0i8h5c84cfPEfnf96w/BOOfE/80GFXqdNxOq/TbL3ma0Ak6lVqStIFCCFx0VovwX2gO26h6fjnuSlI1+SK9Kw7TuK/FFlIwuVlxAIVK5qKIL4VoC5P0ESYeddBsuWm+OFIuizq6cYrP3YX+VtuJu26YoYikit+lHyU21CRVqCezMS1ZwyYy6taWNzq3t7gmPBE+/gdB9l1u3v8/C/fyuQbVQKtw37zTdjUMLMmJVsuRal7E0J7y+TaXUgOkgQXQu1HIOQJcvfDFFfIgW8a33F/x/nOmVt0CAnwBcAwM9KU5s7o3wL8sF/+EeDXcKR7m38M8HHgd0REhp2TNgqMssdsvfjg1ZuEZFq4+GsUMkHSVUzXk29wLiiVRFu9gq0gYymfq5ci8gRPwj6qnhG0pxjfiD2QrwBl/z2NrEX+NU8AyUyPRIT02DSauGz7ODghZr7rIDaV4gftTJB89l5y/OidLCI8Xw1WQVN0G5NlGJUzDxQuhCs2Uti3wn68eyHZu4s9v3sKCxx77y2s+qpl4li9U3QV5h/uwwCn33O4aMk4F/KHHi+O57qQOQIthrrGEkNiXJOgTlJ2aPOfyRBBu+btLgHofriiv6+NftygfOxfy44eHe6kFxLjSroA9c7owFPAMT8JGNwE4ND1eRPwHICqZiJyHFgHvMYSwPZ/5woWpv/uKo48toHJ1w3JtLDsdUe+nVOK6ZXE6xIZpeZb/LrG/XxNGf2WC53UoMbNgMuN0ytNz2LU/0xb17zdTRLwLfviy+goQQIgmSCSQWIwvgnMqJNvT/+HW9j+i2VRR7Jvj9OlvQWuiDqNIX3iDfKVy0by2a9M+Ig9v+FxbMGqJyWDlvvqmX9Ei0kSYb9p+dWKG62vu+doQcZzcemp7z/MxPEhGTdCEfnu3QVEMhhEP86OeMkinSF2RuzY7h8k2MmOu5ry0XuQrcgskmXITFYkCun1IMtIoxLpccVYlwHXO6MDe5tW8/dNoULf90pEbgduB9iyZctQJztOWP6tX2Un7jJ02T9ezgNPbyZ9pcOKFw3pNEx83WIyl3CrjHuvDcgsxgWFx5RRcbCgqcH1NbXqbEDe3pNkBk29swIf8fpouUI89UvfqGoqvfqqyuV0um2rszatmOT0xpW88A0TiJZ9Leow+6/xVXWu1HbX770E8ZTeU6crc8Uq5yRCMn0a/DFJk7JJzZCY/u5DLP+zL5J/5RFXkWZtmTir67iD/h6RF/esiXfDev/EOSeytx0oKiPBRaNP/V+3cMXnLMsfgOS63ehEir3vob79Tb6R+blmwis/eytTG2Hbv6/+P1Q8yEVPBv8BOjXlqu8mO+554QX2bzt6DGDvewhz83WYm68ju9f14DY3XweZdXPqOknZmQxgIsHMeD3YT6VwfSzSRtlm3DC28kKMqDP6YWC1iKQ+2t0MhKajR4ArgSMikgKrcFOB6/u6A7gD4MCBA+P/PzQLTn/TS+zE9QM4+he7OfmldXS+bjBdMBmkp6FzyvoiCakkzwpJQgF8MQXRByImT9/TwfkuFdtJkDC6JRdUbFndJlq1RNWtUgFqa9VhFno58vWMzgOPsvXvGt7w4f2kb5xy+zl+qrzsDMb9GE1fvmC/Ks7DusGLPSlG3ejyiYGNyeXgPpKjJ0GVi+9/ATZvcgQ6NR1NbIgy9fVzqJfqFjuWs4rQiiY3AXnOsoeOFJOEk7270GVp2Vv5lv3k6gjL7L+mkAQwblny8kkfTeaseBwXjVNrihPr1fXKtjR1995CVzRO937ewtMb/2nufRA5sA85uA81BmsEXeGr1xJTIWzTzZHcgKbkDy28XDUv+GB/ITCMe2FQZ/TPAu/GORjeC3zSb/Ip//wL/vW/P5967kJj3TsfY130/PW/dJd4rz6yjmUvCytf0KIng0ukeU04oyikKD4M9UGaofAiqlQKc+IwPsEm1kUmIkjcYCVcXheR9oC2iJEuXJkTFmf8X3y9+Yuu1cixf5/RvpqGhOblJaqcnnHHj8k58WTx/NFIP/SvGeMizU1X9Hubm0i2ycHgkW5Y7/Y3RGMd++YbmXjF/QDEUXrhaMCTcZrA1DT5I5FX2btJFODQ9a4SMc+gq+5yPrfVfrxbNrsiijrRhh8aE8a/R1Vt/vWmSRXB42tuuraItM3N15XNn6LKy6IJlE8aG2sxXZfItfc/TLpn59j7xPswrqTLgM7oIvIwcKeI/AZwP/Bhv/6HgT8QkSdxEe57zsF5Lxqs/U6nr631z5/9+D5mTk5y0SMTSObq2yePgyaQTzgyTmYc+ZogE2iQJkIihOhLhyvWUBCM88570lFL5Mf0SR4JROTDqoq22XApHj9v8rpGTblnJ/EIg9atH8vrhsWxsmojGnf8kkXSjZcXl/V9EkLTOTQ9Ls7buv2JNA6gBC/FPFZ6bAtpobYfjE9ETU407ueZj+1nwx8tY8WRKfd/Z9QF//WG6M8eqRZWFH/HslItdh0UbSypORwyW8wCDMJyUUJcDHkNPmLx/l+JpDF3RWY7BhFFDu4ju/uBxvc2rqhcTZ5nDONeaOyMrqpPAwcblp8Glma/wRFgy7vLD+eln1/D17vLePbYGqamJ0gfX0E6BSteUtIZJZ2yiPEj5oPoa3H/a1aw4iPj4K0sxGNHppK7x4opSo0BLztoKR6HS+Ik+hIH2MERYfG6iaLoerOW8LxJY44fxyWwYXmTHhuv24R6v4PgJpgvYjLfdEXfyB45uA9eeqP/PXgM68md+r7DbPlQl/TEKT8Rw3mxh9a2o0Y+fQh/w1zKnhRd938tKn2ThCXLnOSAvwRLQxlx2QbSbesPHaZ6DHemYwdpCgjOA9qKtDkw6hHnsQXttVvdl/ZymiOpF37pVtJpmDymSLCm9bw3uKhmc9Gv5OVNE+N6PHg7mRhFo4KKghqiMTNAc9RrqOq0TbJCQNGruKG0VRr2ER8n3uegS/+613auaDmg3vmrjlka0wQUEWZ47ZWvw4rlZZmuLwUeFBXXkb39AMl0xspnnDTRlEgDH02HBu71c6v/wDXBKpCXkywO7y+cHnYiaLvupb4J2NHfsZJ0ywGlMk7p/2/v3GPluO77/vnNzN5LSqRIkaIkihRFUaJIvSzzIUpUGtmuklhSmqgtjNQpasuuWwGNDTRAU1hpgDb9I4CToi5cNHWg+CU7iR+xXUgB5MqOI8c1bL3fit7Wi3pbIilS97Uz8+sf55yZs7Oze3cvd/cu7z0f4OLunZ2Zc3bu7nd/8zu/x3HHOPt0A4vHGX9cHzFw6Pr9JLNKPGMF2FZJi9K8EN7Cn5dpKb6ZmlKVrsav+xz7QuSEuKiipvTcdKuO4oMbewLrZ4xV9s/V+Ch7sUJ84e0kwt0E2B3X6dyd9lE1oVFzc2V9XbtPsv0cdDKZNx1bBTSJyO97lOYHL2Xqo/tZ89Uy1C7Zud20J0rT8m6jxVcbFUXR2+rrevUZqkkO+YTNTsxBstz+lF9yGnvRDZG0xPuS0uJiiKeaSA+dmMeVsXUvBMaPtTfXFzfPrtpjY3bVWp2R9Q/64qsgZsHFBcqrloJVxHQ6d4H7sPd6K9bN+nTnabGEK5ans5brdN7f1Rfbqmui+rxPLyLbiapP2yVARGIWsfIM5uZgIiZ/326yFXFbJbu3/80VJuElhZP+wvwfk6mUtU/MgIvKaDbh6FTr2BON1vGKa6Lt11yStmwzR/TjB9Bf3mVivlNF0syUG21mRNOzSJrzzoXrkRzipsmKjOcy634wdUCi2aaJEU/z9nb1xxNBdAPHiivSnr9vd1kkB3vbGIPGXpylGgvX1XQA87c6ZcsFEW1b9e7YEhw6W6m+IMbSemxeUVetuY0uXmCNxVkg9a4Kn7ri4C21aituljbrWVqtTCe8xbFJEWEhWU7j9vs58tuXs/rrZar1ui+03r3EF+2AwzPGnzvXND/+fKPYhNA1EnPrb8O9fOvSxUq7qBDXkr4aq6yTCdljTyH/7wHk8kuImplxt0yYKAbiGM3hpHtfNpET1o8dXXKBcWHNpUZk5+bM7zyvj8E+ThiEpSsiZ2JqzZyOMQtuUtXPdTsmiO4SxBW1zj6whyi1oT2ZtX4F8okYTYR4OgMyk1wRYe55o0JdTTJxRVzV71gApYtCtbUlTVUAq/5YR1WEXXOluu63PrnSloZTtY47xQxXM9biqFVIq8wX6eCfO82QNKPxWhPO3cbJd70GrvOwSJGIUrgEXEsh14nCvVbXUqiRlGGBqRZui2jPRab7w0RMXql+1sb+S0xM7dScKYgURTA1V36hZupFKAgkCemLz5cveS49LkqH9s1gLN0U+A+qer+IrAbuE5EfqGrHwOUgukuY+A5j+VbfW06bDn9kPysOZUwcnPMKWxtfbzSXlhZq3cJWYeFZK6pw10opwO6JXDsLr39ufx/fp1znU61avS2TcHOp+m8rF6ClPq63b3W8Lu1zus6hiOhwx1j/q+bGFZGr8dn6Yzixta4ElwlGFJHHYqImRMjjiGwyJlsRUR+IZkh2nEuWK+mJDTjRZKbFMykymxnL1flwM8+1VPF7+2nMSwZtXdtd8GlUXwWzEq6qR0TkcUwphCC6gXbSFcLsmhhpNkw1M1fXF8izxPrvUnPLWikG07bCXaqusZCdBQytAlnna63ewkN9zYRiF2kXtzp6Wd2vns9PrPAtwbpjepmTb/VXLe+6hTGg6DLhRDeOiz5n6aoGGkHj+/cSUzZdrSNauZL8pVfQJ6fKwIMP7DHjJZFZGEvz1roYtMcHL0WEwS+kichWTHht1wpOQXSXMZ3KBc7+xj7imRw0JpqIkSwnnjK3wabvWxlapL6VCia4H4oC7J4Wm+ftny3pp7X+1xpB95M3qkLpqBPgOuFusVZrjolqjvHH8hfT/HHr4o2LWGXfmq6M1eJyidCJxIiuJ7jc+VBfH9i6GriTT75i/LB5jmaZzTIzc88OHerj7EuATl/W7ZwiIvd6f99kyxgUiMgq4DvA76rqO91OFkR3xAw67ndQ+PHDk39T1pA98i/3m+I6M5O2q4WNDc6UaC4jmsvLBTlfZCWqLJiZX8XTbb7iSqNO1VIKqx+OzHeFVFbvOyVXtCRndEhX9oWv6pqIOohr3bG9LOJ5+xYtd1xDSZeG61na8axtTrqyYb7wIoGfPcTcNZdCJEWnEwC1bXqi2RRpZi199Pw44mjlSvLpaaKVK2unNq7v1UHRh6X7C1Xd2/E8ps74d4C/VNXvzneyILqBrqz+q9bwtFf+4xXEcxA1Icpg8nBOPKsk72YmcaOZ2cQNFwVhD7RJ5BpF9t5Oyhx/S9Gry67+F+LqOh8Ut/m5qxjctqAHtIl3mYHn9q0U5vETQYoMO68AuK1B3CasFesUKha8T1VI/fO4qfoJCC5u9l6TwejHCPhaMfG9e5j79X3MXX1p2XMvh+Rv7503/t9ZwkuyK8R8zHdxekTMrd0XgcdV9bO9HBNEN9AXZ/y3VpfEs5/dj8Yw+fYE8TS4jhmrDijJtLLiYLPofAEUKaXZipi8IcytNu3BswZFHzqAxlSZcRfPKlGmJFNG0OMZ62duZl7xHy80Smm/9Vfbay63KygqXht2zHYnvn72l+s84cK3Kq4AdS3YK2LZMavOT0IQaS2X6KznPtsqFWniYuoj+KUkA50ZUD3dXwI+AjwiIi6M5D+p6m2dDgiiGzgmijKFNbx+y/kcyiJmZybIDk+iiclBjlem7Nj0MqdMvsslJ73E7ReexOa7VjGTNbjnxS2kswkrn54ETCW2E182Qj5xJKLxbk4jEhMKF2FijFVN1ES3GgtRhDbEiq+UwporJlwDauOL3W9X32JyohTTyEubLbaVbgm/U28+mZA3ItJVCdmE+bLJYzj5y52v33zMXX0piK1OBy0LoYH5GVD0wk9oD17sShDdwNA47brHOz6XAq8Br3ESAAcuMx2hzuZg275Pf3kv8k7CqhciVrwtZpGPCJWEqJmXRX9c2FPFwlQrhEUGXmIX/DKxJTGhEF77sPUEnpClqfFXTzbMea2lqo0YnYjJVsbksZCeYCz52ZMiZtbDie9/kzXXPl2Ed9V7UfsjSpV4trUtkWRqirlPTbcV6Ql4KP0spA2UILqBsWfL5l/w6ttreDc/geZqofFuQjKT0zhiOmlIribjLo/Mol6zkiXlbultsodkphg8bvHO3da7DKu2FORKaFcSG8FNIvKVDTSOyFbG1m0SMbM25q1LYNvv/YwT3Sn+ZPDXxbkRol0XGFeFzVZL16+C9atgy4Yl1wF6kITaC4FAByZ/9Xm2As//0RUkR2FulZiuGTlEszbjDjUtjcTEtxbRDVAWbIFyAU+1LPsYm7hiSbNyW64mgcFHrStCTUnEfLJhXAYnNjhyZoN1X/gpCTAJrPlab6+t6AJxdArNMrJDh0yLHCC3LXPmQ2ZSE2KXxjblWwrrnr0XF4txgQpBdAOB7jSOYvrPHa3xXYq1Xm19hFxswEJkmnkW3RZ810OkRkBtyJu6ilppVrodXOZYcQyQpsjRnHimCWtPpPGTB4si9QDJqRuKVjndbvGTndsL61TXrTZhXmwq5ulSfYGuwpk9/pQ5F5nxU0bWfeJe60U7QLXoFhwYTnJEr/TSrqe2oIOIrAO+CcYIAX5LVQ/aEIrPAdcCU8DHVPX+4Uw/MCiG2bZ+UEwegnim/KRoBPlkVLSUiTLPss2VuZPs21sAq69RBtioCBNjbOKNJVcktR0n0tgIcTOCLDYdK5z4uk66kcJMSvxmE7ZsbikkXu2xlmzeZKqERbZ+rctSa6Zljzv/AJGixY5Yn7HsuoCiCwhYa9zublv7SFqJG27JhtOiNm/HoufLCdWxLmJeW9AB+BjwQ1X9jIjcCNwIfBq4Bthufy4DPm9/BwIL4uC/3s8Jb2acetfhjsW+wWTSudZGuLrCkW2qGBmrN4sBFdPBfkJMPYPcfAzimZwoU+LpzNSZbWb2PFlZxnB2thRfFwMMXVv7pAdebtuW7DjXhKDVlUbMMiPIkcC0ti+N+41A65IuCh81rYkhaQp5fly0Rx8J42rpdinocB3wfrvbzcCPMKJ7HfBV24zyThFZKyIb7XkCgb455a63kCNTMDvXFljg42fSOab++WUt2WMnfLdrWnyBi53P37ebeLppRVeJIoHZpql5C231Fupa+9Th2vEk289pOb4om2h90t2aY7oeaK5yGZjmlUWCRzU92i4UapoSr127/NJ+K4yte8GnUtDhNCekqvqqiLiufJsAr+UpB+y2FtEVkRuAGwC2bNmygKkHlgvy7qyp79pj37GDH9+PqE2qcP3PEU786zu7HVZL9Pf3o0B00Y7ORbu9CIi+w7ScOM7MGUvUinkvr9UX22Kbax+06YyWlOUQPlZB6b0w/4DpWXSrBR1Eam5r7K4129penS0YcRPA3r17Q0R3oCPVljPzobFZSDOt7mHiiHLCd3oX3OyqPYVfWHJFmjnazIAmMjPXWrvBLVjZx8kZG1sX7OzvOtFracnTNLf+/caOxqtXFwkZ2eHDxfb05VdI1q8nfeutvs63rBhnS7dDQYfXndtARDYC7qv5AHCmd/hmIHzNBkaG687w1r+9AhowG0HzI/tZ87Xu2V/v/Kv9AEUqMpiFt8a7OfF0RtScZCJTZK5pFteqtRh6rXqGdSs0U1O8PM8Lwe3H1xqvO7ml+0W87mQ7vO2PFgS3K4vlXpi342CXgg63Atfbx9cDt3jbPyqGy4HDwZ8bWAzW//lPkQyyCSGbEI7+i8u77p8nkE1A80QhmxTWfO1nrP6rn9E4mpoFtjvuI3v8KXNbX1e8xsePE642wwSYmobpmWJxaz4LN9l0BsnG01t+pNFAJifM7yRB4hiJY1RzNE27ni/g3cnM8zNoerF0aws6AJ8BviUinwBeBFyc0G2YcLFnMCFjHx/ojANjz7iEn735O1cAxsWw7ov1tYNbsLVw4lltqYkw+cwbpC+UyxTJ5k3tItlWzJ2O0QzJ5k3tQqtK+uYvaqeVnLGx7DBRO28XKuf9jo+hg/NyYEBVxhZCL9EL3Qo6XFWzvwKfPMZ5BQLHzIb/3YPQelS7LEcnnEA+NVUIbrJtq/HB1gmujQ7otABWRBV4rgQzSNlzLjl1Q+EPbnEzuOw5KMpcFtuh9BsHd0LPmOSIxVHd8HUYCHhkV+0pHudTU/U7VX22Xs3fjriykABZRvrqa63hYEli6vbaTLbktFOLp1r2q3NVqPYc2RHwyHv8GTAhDTgQsORX7iI52kRsS3OZs7f0frvxaqdjT3C7WrlNayFnrVXBnKC65ArAiK8qycbTi+fTV18rIyNcfQh3jpDosCCCpRsIjAEaCdFsE5mxEQpNGzubpkZ8vcSFFvHsJLiuoI1Ldnj9jd5E0i7UJRtPL8d45dX2iInAwtA+fgZMsHQDAUv04wfMg53b7QZrUUpE0b7HUW1IWSHZttU8sN0p6hbUas9XI6TJGRs7Hx+Ed4GMd+2FQGBZELusM7+Db+Q1OZdKRpq3GFatvZD+/Pm+xvZFNdm8qTx/XSdjOyYwvy850JlQxDywVBhWlbFhkV+5C4D44LQptVjpYix+W8iWAjOds82OBVcgJ9my2TbVlNa03grJaaeGhbR+0YH1SOubILqBZU/hVrjwvJbtrv5sctaZna2iODLi6CziukIz9G/5WX8klAAAEDVJREFUmljeiqVbV8DG7R+Et3/CQlogYIhXrSI64YT27atXm1oDwxjzwvOQrLU4eny+EeH0hZdIXzxgisnktpyj3z0YTD0GrTGdFuo39C3qOII4Mha1y4Tzfcp5DmlWpAEHeiQspAWWO76gdoyRPQZk38UAppXNna29w/JVK4iOzpRxmbnpyFDFr42bbDoD16nCDFCxYaxvuF8r15zLdhJ259W8tLhzWkU3LUPbimuoimYZZBn53Fz/4y8DunaPHiJBdANjQdWCjVevJjtyhHjNmpbb6OzIkQWdP9m5HQ5OlRZktY7tW0crrXy6RwUUoWBQiGILC3UrbD2r5lze3yLe/alrmeG6HU/MW4M3YFEGlvggIl8C/gnwhqpeNN/+QXQD40FNplW8ejXkeVehTXacax64Tg7OB5pmRqyc9Tk902qJLtQCreLO6Z/bK+2YnLvNRDU8/Wzn17D1rPZz+qJfdWO41+gSNbqUjwzUI+ggkyO+AvwvTFuzeQmiGxgp8dq19U9UPgB+bdg6knPONrfV0zNlGJcXwlXXImfeuV2804z9yBPdx956Vr0oOhGsWsmqppSjKukzPy/nX8VPthjEF0KgOwMSXVX9sW3w0BNBdEfM8RZONSiKRR6vUWId3cS2CJtShaPvAr11WOiV+cQWOohlROmaSJzlW9+7LDl3W+15nRgHRkjvonuKiNzr/X2TbcKwIILoBoZGsn49qjkiEer8kpVb4k59upLTTm1PtV3k2+fCDVAVVD+aoM4X7LsHKsd2czsshEGW1VzS9OfT/YWq7h3U0EF0A8eEvwDmfK+uo4Ha4i4q2m5VdElfTTacAmlWCHX29sHBT7wPOlmnRZvz+agUyRm00AYWxthGL9StzInIOuCbwFbgeeC3VPWg7TLxOUwR8yngY6p6/3CmHhgF88bFeoH7xb7NsuB2dREsXrOm5bg6H6820/HqVJvnrbURRMoi4R2+TNInnh7hBAP9U2MIjIheLN2v0L4ydyPwQ1X9jIjcaP/+NHANsN3+XAZ83v4OjDkdxbXLGzM7erTvceZbIBslrv15T5anE1s/48x3fdiW6oHjBGVgoisiXwfej/H9HgD+i6p+sdP+vXSOqFuZu84OAnAz8COM6F4HfNV2j7hTRNa65pX9vYzAMJnPel1oLOw44CxpX9zd6y3cHxftMG6P3NRULcLOoDUSAUAEjaXw/xVhRi4ON4jt8cuAvAuq+tv97L9Qn+5pTkhtN2BX5n4T8JK33wG7LYjuIlMV2uNZWDvhp8H6bouqqyJ79Mlyv4t3et0fqE2Ml6zmVnSRygIGBsdiFTEf9EJa3apC7SsTkRuAGwC2bNky4GksbzpZsktRaIGytxh0bO5YR3zhebZnWev2qjVbyyJ9YAMD5DgT3ded20BENgIuWPIA4OVHshmojfOxcW43Aezduze8gxdIsTDlMU5+01GQvvGmEd4eC3q7QjZAi+C2WT51tWxzDdEHSwHV1trII2ShonsrcD2mDfv1wC3e9k+JyDcwC2iHgz93sNSt9o/VSv8i4Jo4zpcoEbvSjf5nLTJ/ix8TnGbmA9lIugpsiIk9zhlXS7duZQ4jtt8SkU8ALwIuzeo2TLjYM5iQsY8PYc7LiqrILneBdSQbTy9amXfsT+Zlj6XPPleIq8Ze+BfgPGBhUWyZMa6i22Vl7qqafRX45LFOarlTCK21vpaDuyBeuxYiQfxiL0VqbVLe6tv6sXWVtJJTN7Q0fUyffc5s37ndVBlTRf3sMR9bHyG4DpYJyqIthoaMtHGk8g28HKzdY31NyaYzuvt0VdGJpOVvd53FVQWLbOyus6BD0ZkljNYXnR8BQXTHkOzwYeI1a2ot3Hjt2s7PVRfVqmmOFTHXPCefmqrt0uD2z6en+5r7qEl2nNsSSxuff14Z9mXLReZR1Oq7zbLyMfSVohxYIijH3UJaYMh0cil0swh7dUPEq1a1/D2MLg2jINp9IbkI+aR9G//0wWJ7S1+xTMvFM7u9ritEcu42G0JmrKBk61mkz78w7JcRWCzG1acbWHosJH13HMnvfwyA7AN7iGcz5LL3mBbqae6Ff/X+wQrlFZcZQXQD44TL7lrsCl/daP7a3rJE5GTExP+9B9l1AUSguSu16B2Qm+wd1+U3sJwZ74I3gWVGvO7ksRPb2d/YVzxWAURYcctdLfvEF++EzIaExSZKQT3/rD7wyKimGxh3lPY1jxERRDdQUCzE2a6y45A2PPObpdiuuPXuefdXEYgiNBbT9ZfFy7EPjDnB0g0MimT9+vKPyIhQ+vobpji4w4tRdMXGzfbujSBHRfODl6IRRE1l4nv39HRMPtmAWNBqYXH7WqPdFxZ+4MBy5/hLAw6MMelbb9U/UVMHVjW3yQimu8F8gutX8ira8LR0q12YeM/9+j7IlLhpPgiN23sTWoBoz0XGso2NhWvcChClajpY3G3cCsHeDRQoZQupERNE9zgnXrt23sSCZOPp5oGzaHPtqxVOsnlTKaZxBFFctKDJBhRSlRxNkSxHfvJgz8dEl1xgLNskQiMhn4jRSMo03ygHra3WGAiEjLRAfyTr1xuLtkuPrmTLZhtz6jKuop6bOybbtprkgrgiWfaNOoh0Wdl7MWITFfj79q5OybatpXWe5W0xs/lD/2DOs+9is7omWF8uIEJOZGrhBgJ1BJ/ueDJuLdPjdScjEhUuBInj+Q9KzD7pCy/Ns2NZiavtzis6tlCrwtp2Qj4xgb47W5ukEO26wHxB3PdoT+d27gP55V2IKHkSmegFQCPo4QoFlhu2hsdiEER3DCgyxLqknzr/qMRxUai7qLTViTQrH/fStZZSWF0LG3UptJgUW1ElfeJpY0WL9CTk5RxKSzvacxGaRHDFe81zNpssu2oPWQ7xHfcVh81dc2kxf7XuZ/Vfji2GI9VsXqDx/XtbtoVyjIGCYOkuffwaB37qbS8ZYslppxZvkuTUDWZjmrZ0SmiJWgBoNs1vkZbqW92QfRcjmaJzadkF17WywYRk+dW4kh3nolGErkjQiQS96+H6E881W4Q/v+9RuwBm6yNcuQtEzNgipL+ytxBXyUFVQUAik9YrLW3NtTx3rkje30JcYDmiJmpnEQiiO0J6qXHgbsP90oXJhlNM1av5hHPtSS3VsVxxbzBVuKr+3OwDexBVoh+V/lS9+5G2Vf74oh2AloW+Peva1aBNtp4FU1O45ust5Smpj6jI73uUaPeFaBQR5ZExTUUKq7df8vftJqrxDQcCbSy10o4icjXwOYw77Quq+plhjHO8k2w4pWzt7Ra7KAW36P2VK+lbnXt/xReeZ6xZr+EiAI2G8Z8mSeHXBbuApUr09jSiiuzcXrgi6uoPZI8+adwNuUIzhTRt20ffecfUwnXH9FiqMb//MeL3nA8zxqrWRozuvdhEJYjAnQ/1dB4gCG6gPwYUMtav3g1cdEUkBv4U+FVMz7R7RORWVf2HQY91vNOpiaJzE6gVt25hXcn2c0g7LXBFAhIVQeBFGcQOC1hA4attQRWmZ+ykM8iy0sVhvyx6dV/UkT38OMlZZ0KWkx14ecHnCQR6RQEdgKW7EL0bhqW7D3hGVX9uJ/UN4DogiG4PxOtOLsW2h1KNXUO3Gg27Uwppj4W5RYxIN5tlmUNH5HVdEJm3J1k/9LUgFwgcKzqwIuZ9690wRHcT4H+CDmCaVAZ6IcsH1p4nffY50ycszayAzpJsPL10Xzj3hs/UtHVnGB9stQVOILBUGNBCWt96NwzRrYtNarPjReQG4Ab756yI9BaUOVxOATo7T0fBoQHPoVuvxe6GqpnH6wObyUIZ6PWQhXWFWPz3hSHMw3DWsZ7gCAdv/1v99inz7wnAChHxYw9vUtWb7OOe9M5nGKJ7ADjT+3sz0JYGZSd9E4CI3Kuqe4cwl74Yh3mMwxzCPMZvDmEeg0VVrx7QqXrSO59hpKXfA2wXkbNFZAL4MHDrEMYJBAKBxaZvvRu4pauqqYh8CrgdE0LxJVUN9fQCgcCSYyF6N5Q4XVW9Dbitj0Numn+XkTAO8xiHOUCYh884zAHCPMaSfvVONFTVDwQCgZERSo0GAoHACFl00RWRq0XkSRF5RkRuHOG4z4vIIyLyoAsHEZF1IvIDEXna/j55vvMsYNwvicgbfohcp3HF8D/ttXlYRHYPeR5/KCIv22vyoIhc6z33+3YeT4rIBwc0hzNF5A4ReVxEHhORf2+3j/R6dJnHyK6HiKwQkbtF5CE7h/9qt58tInfZa/FNu1iDiEzav5+xz2891jnMM4+viMhz3rV4r90+tPfokkVVF+0H43h+FtgGTAAPAReMaOzngVMq2/4EuNE+vhH44yGMeyWwG3h0vnGBa4HvYWIBLwfuGvI8/hD4vZp9L7D/m0ngbPs/iwcwh43Abvt4NfCUHWuk16PLPEZ2PexrWmUfN4C77Gv8FvBhu/3PgH9nH/8O8Gf28YeBbw7oWnSax1eAD9XsP7T36FL9WWxLt0ihU9U5wKXQLRbXATfbxzcD/3TQA6jqj4G3exz3OuCrargTWCsiG4c4j05cB3xDVWdV9TlMysW+eY7pZQ6vqur99vER4HFMhs9Ir0eXeXRi4NfDviZX47NhfxT4x8C37fbqtXDX6NvAVbLAzI8e59GJob1HlyqLLbp1KXTd3uyDRIHvi8h9YrLjAE5T1VfBfBCBUzsePVg6jbsY1+dT9jbxS557ZejzsLfHuzCW1aJdj8o8YITXQ0RiEXkQkyv4A4wFfUhVXVk3f5xiDvb5w0CloPJg5qGq7lr8kb0W/0NEJqvzqJljoIbFFt2+U+gGyC+p6m7gGuCTInLliMbth1Ffn88D5wDvBV4F/vso5iEiq4DvAL+rqu9023XE8xjp9VDVTFXfi8lq2gec32WcoV2L6jxE5CLg94GdwKXAOuDTw57HUmWxRbfvFLpBoaqv2N9vAP8H8yZ/3d0a2d+DK6PVnU7jjvT6qOrr9gOXA39Oecs8tHmISAMjdH+pqt+1m0d+PermsRjXw457CPgRxke6VkRcPL0/TjEH+/waencX9TuPq60LRlV1FvgyI7oWS5HFFt1FSRkWkRNFZLV7DPwa8Kgd+3q72/XALcOei6XTuLcCH7UrxJcDh91t9zCo+OL+GeaauHl82K6Ynw1sB+4ewHgCfBF4XFU/6z010uvRaR6jvB4iskFE1trHK4FfwfiW7wA+ZHerXgt3jT4E/J2qHrOF2WEeT3hfgoLxK/vXYmTv0SXBYq/kYVY/n8L4r/5gRGNuw6w+PwQ85sbF+MR+CDxtf68bwthfx9yqNjFWwic6jYu5dftTe20eAfYOeR5fs+M8jPkwbfT2/wM7jyeBawY0h3+EuRV9GHjQ/lw76uvRZR4jux7Ae4AH7FiPAv/Ze6/ejVms+2tg0m5fYf9+xj6/bUDXotM8/s5ei0eBv6CMcBjae3Sp/oSMtEAgEBghi+1eCAQCgWVFEN1AIBAYIUF0A4FAYIQE0Q0EAoEREkQ3EAgERkgQ3UAgEBghQXQDgUBghATRDQQCgRHy/wFYurQYCvUwhgAAAABJRU5ErkJggg==\n", "text/plain": [ "<Figure size 432x288 with 2 Axes>" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(np.mean(f.variables['diatoms'][:,0,:,:],0))\n", "plt.colorbar()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "<matplotlib.colorbar.Colorbar at 0x7f3f419f4198>" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWYAAAD9CAYAAACP8N0iAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzsvXmcJFd15/s9EVHdLQmBdhCSeq9Wt1YkNWqp/QyyAQ9gG4EHG8EzA1i27PfQeHt+D+EZA4P93hOzMbyBD7YGsBEztli80LZlY4yNlwfCakkGJLS1Wi2pxaK90dLdVRFx5o97b+SNyIjIyKqsrsys+/186lOZEZE3bmVl/uLEuWcRVSUQCAQC40O03BMIBAKBQJkgzIFAIDBmBGEOBAKBMSMIcyAQCIwZQZgDgUBgzAjCHAgEAmNGJ2EWkV8SkdtF5A4R+WW77QQR+aKI3Gt/H2+3i4j8fyKyR0S+ISIXLuUfEAgEAtPGQGEWkXOAnwMuBs4HfkxEZoFrgC+p6izwJfsc4DXArP25CvjoEsw7EAgEppYuFvM24CZVfU5VU+DvgDcAlwOftMd8Eni9fXw5cL0abgKOE5FTRzzvQCAQGAtE5NUicrf1ElxTs/+DIvLP9uceEXlq0JhdhPl24GUicqKIHA28FjgDeKGqfgfA/j7FHn8a8JD3+v12WyAQCEwVIhIDH8F4Cs4C3iwiZ/nHqOqvqOpLVPUlwH8F/mjQuMmgA1T1ThH5APBF4Bng60DaNte6YfoOErkK4+rgmGOOuWjr1q2DptKZe27ZO7KxAsvPlos2LvcU+hjXz9g4vlej5pZbbnlMVU9ezBj/4oeO0cefyLqd7xuHv6Cqr27YfTGwR1X3AojIDRivwbcajn8z8N5B5xwozACq+nHg4/bE/w/GCv6eiJyqqt+xropH7OH7MRa143Tg2zVjXgdcB7B9+3bdvXt3l6l04lXRT45srMDwfPeXd6KVT9ap//ErCx7vi7s/u8gZjZ6l/Ixlr7gIcshnhEMnJmgsPP+/f7XTa8fxvRo1IvLAYsd4/ImMf/rC2k7Hxqfeu1VEfIG6zuoX1HsIdtSNIyLrgA3A3ww6ZydhFpFTVPUREVkL/ARwqT3B24Br7e/P28N3AVfbK8cO4IBzeYw1EiGRIMkMRD2jPz94cBknNV488Js7Uef8EkBB1P7OzW+NzLa17124EK9EDr5hB5KD5IoKaCxoJOb9zZd7dtOHeVs7v7GPqer2hn2dPASWK4DPqepAU72TMAN/KCInAvPAO1X1SRG5FviMiFwJPAg4E+JGjB96D/Ac8I6O51heNEczjEATIzMJiJCcfBJEEaxeBUmCrl4FSYTMpaR371nuWS8pD797JxqDxpDHZtv6fxMEdzE8+TOXgkKU0buwQfFVNhc6JZ4zIp0cBJVQAXLUKMr8YH3sQicPgeUK4J1dBu3qyvjBmm2PA6+o2a5dTz6OaDqPpvNEugriGLIEEYE0g1yNoRit4vGLT+bZ153CoQsOctrJT7L6VfuWe+oj4bu/spNsDWRrAIW17wtC3MRTb7+0bBtZ20nFE9wK4rRAe8ci8ILru7kqAqNjCIu5jZuBWRHZADyMEd+3VA8SkTOB44FO/+iuFvOKI5+bAyBevdpsEDEujtiYjsfd/QzZqmOJbj6Kb59yFPoHx3H0UeY1uQppFjF/OCGfjyD37nZU0BwkVuI1GUcffZgXvb5pneDI8dBv7ITIWsYC634jCHIbVVHWiMK9U6JBoF/wqSDEy4miZCOoRa+qqYhcDXwBiIFPqOodIvJ+YLeq7rKHvhm4QTsWwA/CPIDsKRNymLzwFCCB+RSZT9E77uGEf1r4uA/+u51E85DqUTzxszuJUkVyOPwC4ZT/emRF8YHf2kk+o6DChmuCIHfhuN8zwvrkOy4t7xA4/hNBdCeBvNEVPByqeiPGhetve0/l+fuGGTMIc0fS7z1CcsrJxnKOI5KN60GE9L77FzRe3eLYoct38ILrv7bImQ7HQ+81oqwCG9813aI8ykiKJ67ciaiiYu6GxBpCQZQnAwWyEQnzUjCVwvzFfHQhQ/6XOX3k0d5iYBybxcH164yLQ6QQbQBNIrI77hnqXPGhpVl+f/jXd6IC+Qxgf6v1rkgOG94VxGRYTvi4uYg98bM7AVDrYH7yZ4wFLXnPqh7EKD+vge6MymJeCqZSmJeS9NHHSs+TdWeAt7ib7nmIZPNGBEjO3DyUQM984WYOvn4HR/3JaK1myYCot/C0PoSyjYwTPvYVI85u4Q+xlvTyzivQjgLzY9zvNAjzIkkf6MWWJxvXGws6zSA31q9kEcmZm8mPXoXOxOSrEzQSstURiJhFo0hY9efGYX3Un3yNQ6+7mDW7Fu7AfuwXTLyxWAP8xR8IQrwUOGsZqF3ke/Idl3L874a7kXFE0eDKWCmke/fVbk82byR65jDEEXrXvQjtRUo0Eg7/+MXmsQAirPl8vRV96PIdSKb2dbBm1z9x0m8fOSG+9+MvhXlBDkds/qWbjth5l4vHr7JiXBuD7H3RwyLgeKOQja8uB2E+EqR7TF2FZHYTyZmbzUaRYuEou7Ps6hjGldEk2AD5yy4AIPr724aZbifu/eRFJN9ZzYv+OufYP5huQX7sF3b2iXCBfe58zoHJYNwTKoMwH0HSe+8rHidnbjYWlgjJ1tlCpKtURXsYqoKc/dBFEEH8pVsWPObjP2fqYGz9T0+SfWPh40wKj/7vRpSVfgtZBU78WBDkyUTIarOpx4MgzMtEWzp3vG1L7eMqw4r24RNmiFIle81LC/eH5BDNe7aDqhEiMSnBiIAq0VyGZMopX30cmc8gG2d7Y3REGbWm1Yn/LQjyJGMW/4IwB4agTXDbRDs/ZjX5mphnT1+DRtIrOIQR4GN/v7vPs/qRtbWKAOMzJ89JNm3oi+NO1p0BYk/sZUr6dwvDcqSrBT7+8zt7YuyJchDj6cHEMQdhnlhGHWO6WJHxRTvZOlvaFx08jO6+h+ct6gyDSffsNdEnqPmtFZNSc5CocTH0SHPgrTa22Lkgol6NChWz2Fr6jto/x1V7a4qsCPHHk00eLObAUpDedS/gCXTes6IX45tuI9m43jyoyXpM1p7em9sDiy6ZWzD3oxcX4YTD8OwbLwFg5jkbteLEOHK/zV1FHmu5AJF9fQh1m16CxRxYcpxAQ0+kfWva3z+IZHYTaM+3DPSe04swqZ3Hg/sLcU7Wnm7E+4GHGo9vY/5HtiMKkirJ0/PoD15g/OK5Fxbh3DUiRVihRkI+Y7IwZ541GTUaly3i1X+6iCIngalAEbJOnfWWhyDMU0Z61719Lo7ieZdMp1xNWrkN5xvW8vbFeUGifOn5IMLqxw9BmiN5DmnFVWLF2P1oYucbCxpFxIeNYGscoYksKgolML1MvCtDRH4F+FnMHcA3McXvTwVuAE4AbgXeqqpzIrIauB64CHgceJOq7hv91ANN9FnQFUFujQg5dysaReRfX3gpUifOyemnke5/uLQvOf00UCV9uFxL3BWF4tGnm8dtsdarWK9FIFCLIsxpvNzTaGSgMIvIacAvAmep6kER+QymGPRrgQ+q6g0i8tvAlcBH7e8nVXWziFwBfAB405L9BYFGkjM3l0S5S8eV7Jt3ASZeOVsNp3x4YZEI6YP7zRyc39m5QipWdOGzduT9sWnjsogYmB5MgsnkuzIS4CgRmQeOBr4D/DC9Sv2fBN6HEebL7WOAzwEfFhHpWiA6sHjibVuMC8AnH+7tH1VomBPoVvzkmg6lVJNNG3rjL7Ds6igYx6a/IVKkOxO9+KeqD4vIf8T09TsI/BVwC/CUqqb2sP2YbrHgdY211f0PACcC5bJsgSXD+YWL9G+ASMzCnk/btbJm30JE0Ilo42sbMh7bSO+7vxi3LpY6EBiEqpDpBFvMInI8xgreADwFfBZ4Tc2h7pvcqWusiFwFXAWwdm23NuKBIXFWsq0XrasidFVCftvC/McLslSt8CabN/ZHe1ToOua4WMyBySafZIsZeCVwv6o+CiAifwTsBI4TkcRazX5nWNc1dr+IJMALgCeqg6rqdcB1ANu3bw9ujhFSWMZWkN1i4NyPXky6JuLx1+3k8Ik5m39xuOJDVUu16RjwhNjhh96xcEEtzqsafM+BBWMW/8Y3KK3LzB4ELhGRozGujFcAu4G/Bd6Iicx4G/B5e/wu+/yrdv/fBP/yMpEr6b29Bb9Vf/5PrAKO/sPuQ0QXnEV0cB7m5s0GJ6wDRNGPoEg2byzvqxHlJpdH0wVgIS6QQMAx8Yt/qvo1EfkcJiQuBW7DWLp/DtwgIr9lt33cvuTjwKdEZA/GUr5iKSYe6KcQQCeeQ4SXNZHf9i1yytETC7JUPSHts6ZrqApycFkERk026XHMqvpe4L2VzXuBi2uOPQSM33L1SqBUqH20H7p07z4jzrmpj5Hu65Zy7S4WRU1qz3puE1vfbQLmwjAurouiPohEEMnYzCvQnVFm/onIq4EPATHwMVW9tuaYn8JEqynwdVV9S/UYn/F1sgQ6k2zaUBblJO4Uszws6d59tnhR5fyzm4rqcYV/W5V0z95aQQZMHeq66nRN1rQ93o29mA7lw1A6Z4m475hg1U8W+QiiMkQkBj4CvAqzvnaziOxS1W95x8wC7wZ+QFWfFJFTBo0bhHlCiS46h/yW282TRS6oDX9yKRcz8kp6tpb39K14O+c+4fPirxtbddkLkZvDYi3WRj+2ajlBp+Y8ycb1vb+lmixT87q6C1vXu49ijBee0puf97vaKDjQjCliNBKL+WJgj6ruBRCRGzBRbH7o088BH1HVJwFU9ZFBgwZhnlDyW27vs5RHJVStuEiPIS4CyeaNPVH2xaTjQmKVIvLD/r0lQWxy4QwZs911Xo0XD29OTYJd2l9NAHKlVIsLlhf6CIsqEBUwroz50aRkF3kblv3AjsoxWwBE5P/H3Gq9T1X/sm3QIMwTjC+OhSVWrY08JMnmjeWIig6LcFU/ct2+6oKkm++w1mJpLhVR9K3XTuS6qPO3EfzO440qwySYnCQiu73n19lwX+iWt5EAs8BlmNDifxCRc1T1qaYTBmGeApKN6yHLis4hXUWv9hZedagEDt937It6yadck1SS7nug9rZ+MTgx7Pz3b1xfslaXSqQD44gMk2DymKpub9jn8jYcfk6Hf8xNqjoP3C8id2OE+uamEwZhngLSvftMSydnLeeZWSyzbZ6qgpNs2lCfgefXXW4Q5OK1/vn37C2EuG+Rz43bwDARHm1jQE9Y28YrRYW4i4hzidhoi+AimH6UoSzmNm4GZkVkA/AwJjy4GnHxJ8Cbgd8TkZMwro3WWNYgzFOCE5Oimpsnsn3H3nd/v4BWFpIK/3XdGBWh7YsK8c7TOmdrNS9UnLv41GsvFFTcLu5vjIA8Ill7erfiS4GJZhSLf7Ye0NXAFzD+40+o6h0i8n5gt6rusvt+RES+BWTA/6mqj7eNG4R5ikheeIrpXu2s3v0t4jLIF1vTvcTRJrhDh47ZxaxGt0aDBesfX41xLrloVIePVonGt05vYDQoMrJC+ap6I3BjZdt7vMcK/Kr96UQQ5mlizWojpLm2uw8aY3MbcMd1SFoZWgQ9a3WYBbN03wONkQ9Dz6HG/x2YbhSYn/BaGYERMsp6ubX1gEWoKeZXMGzkQl/kw6YNfZEbw5CsX1eE3PkhdAuJYigW+5woDxll0V+k37px1p1h5jTAnRFqH08yMtn1mAMThLOWa7qAQIu7oK6IftTyofUiN4Z3E/TicAsWmD7eJ6xtc/Zf5y4OFYKlvHJQRpP5t1QEYZ4CipCvPDOCoybiwnWqRiIvYkPLotTU2cRZj16GX4mFFgysJpp49ZqrVnhrYkZDXeeBr/FZYIJLYDoIFnNgSUn37jMWrA2TI/YsAVXQzPPldhDlOnwRHOAyaIojrgu1K+1viKCovQhUFyWHKOAUaloEVCVYzIEjQBKbiIxq5p8TKZt8Mmzvv+o4RTp0Q0ZgqSJcNYGlai1Xt7edu+Y1RQW6I1TQKDA9mMW/8Y2+CcI8BcRnb4E0N9ZwHgFenYUi7C3rF2mfSPrrMfhUC/VUXRFNIuwf6xcEqlaVsxXq/Ep1VaoW9WIWIQMrnQnv+RcYb+JztyJzKaRZvzXcdVEtqnFz1OGNVy0kNAofdLJ5Y2t1uiDC7Yy6a/c0R52Yxb8J9jGLyJnAp71NG4H3ANfb7euBfcBP2Vqjgika/VrgOeDtqnrraKcdcGTfvItk6yxEUS8qo84irtIxggFoFfhSYsfG9f3HNrkvquO4es7rTNkBP6lkmgUisHyMqlD+UjBwZqp6t6q+RFVfAlyEEds/Bq4BvqSqs8CX7HMwHbRn7c9VwEeXYuKBHuld9/YsTd/6bfsZhIszbhDlqi950MJeV9IHHgq1KgJLjsv86/KzHAzryngFcJ+qPiAil2PK2AF8Evgy8C5MkejrbRriTSJynIicqqrfGdGcA024ULhhrOE66gS2Zlu1pOfAMUNP3sAYMc7NWIed2RXAH9jHL3Ria3+7dil1haNPqw4kIleJyG4R2f3oo48OOY1ALZp3i7rwreG6n7pj214bRf1RETUF8UsNWQcUjw8ElhJVmM+jTj/LQeezisgq4HXAIIdfnbnWpxaqep2qblfV7SeffHLXaQRaSB94yMQOt7ggOtMk0g0ks5uMn9v+pHv3kd53vxHsBks5iHNguTCujKjTz3IwjCvjNcCtqvo9+/x7zkUhIqcCro9Vl8LRgSUkve9+I5RZ1n5gm3hX9w3pHknO3Nx74lvxFZfGOHW/XokkmzdCHC1J895xZ5wz/4a5HLyZnhsDYBfwNvv4bcDnve3/SgyXAAeCf/nI09oU1dHk8/VFueqzHtb14Y/hfiqvDZbz8hBv24LOJOhMTHzetuWezhHFhctN9OKfiByNac/9897ma4HPiMiVwIOAC6K8ERMqtwcTwfGOkc02MBy+deonj1R9vlUhbhoL0IaFQfFFftAiXyRFDowL83P1PupSvUcZnztuoXejjj0ehuzOewCIzzkTVIkuOIv8tl5z52l+32EKUrJV9TngxMq2xzFRGtVjFXjnSGYXWBRFp5LcZQVqfxGjuhoTNeJbEmQvuVAWGmnh5uBnkEejadIaGI7s9ruJz90KOUTnn1XcR/siPY0M0fPviDO+l4zASEj37O0X2rwSMeH/5FqOpLCIau8n6z2uHad6rmqkiNuWVyz3wLKRffOu3hN7sYzOP8sI9RRiojLiTj/LQUjJXgn4Pe3qqLaRymmO7GhydQwK06vsr/rAk80b7TkVsrwoVg+E/ntHiPzrPQvZt5zj87ahIqX9k84oW0stBcFiXgGke/b2uywGxTCrLbhfZwVD2erNayzvlp9aUfbnFke9tPIsJzn1RSN8NwJd8EW4dl1hCsiRTj/LQbCYVwhFGyZnmQ7yDbctCPquEP93R/oq0dWVAY0EoqT4hIawuuUnPndr2eUxwYx7EaNgMa8w0j17yZ+3Gj1qlbVMpWelNv04kc6Nm4E877em2/zMg/ZXLfaQuj0WFIt/TiVEzCLhlDAtCSaBMaTwxa5e3S12GSCKMA2CE0jz9siKzLWk8mo8d6Xm2L7mrr4bo2XsUVnLyxmeNimUFvz8qJmmi+4EoiqkIxJdEXk1pqJmDHxMVa+t7H878B+Ah+2mD6vqx9rGDMI84aQPPGSy7ERMPCoY4a0s1kQXnAW5DW+bS81GVYhA894tnTgBrgrysHQQ5T4aXCzBhWFIZjehSVy+zx3G/9u1AW7VRWUXhZ0bY1pcGqNwZYhIDHwEk+exH7hZRHapanWl9NOqenXXcYMwTwEunTY+e4vZkJsli+K2UxXmvfTsIunECnVdmNyILKOhRbWm84nLDFzJAh1v29IrOGP/b+ld99YfaxNGECG7/e7BY/tZf97/vUl8p0GUR+hjvhjYo6p7AUTkBkyFzUWFsARhnibqQtyc2JayAJdIkBdQ2tN1JSlcGm6MINAF8bYtxWOXrddGFzEuHf+NO4ee0zQwImGuq6a5o+a4fykiLwPuAX5FVVuLjofFv2kiisoLaVXBzTGinOfGVdG0kLcQ/KzBBZT3TPfsHRzWN+SYk0TRbKCFBWdZBvoYslD+Sa5Esf25yhuqSzXNPwXWq+p5wF9j6te3EizmaaWwiHubCit5sf7jIRnG0k337O3vhlJXpH/Thl4CygRa0NXGsqjWdvp21nIQ5dEzRIzyY6q6vWHfwGqatnyF478BHxh0wiDMU0DJl6wKmfkSF0Vqtm0xX+zMs4x9d8EovvRthZHs/q4C3ddB23dz+LjqdNVOKoWLJh+bNlV9QmwZ1GA2u/OekiujrYt4oDuqkI6mCP7NwKyIbMBEXVwBvMU/oNLB6XXAQN9REOZpIM9LlnHVD5ndeY+pz+xbyEuVzdWhlKizdrtaulXxKm75+8qLel+02MzFF0R/HL+7d52luljceUfR2dv9P5OtsxBHvVrXlUSf0EV8OEbhY1bVVESuBr6A+dR9QlXvEJH3A7tVdRfwiyLyOiAFngDePmjcIMwTTnzu1iIWuWmVHkxtimTzxnrhPBKJHTW+Y99XPIw7wglpqSu3v/DotsWuTJpdPNy8sSde1QtFzYVqWKHrs4pFzAXRZ5jyqNX53HWvEWd/W7CeF8Qoa2Wo6o2Ycsf+tvd4j98NvHuYMYMwTzjZN+8yX1ZVkq2zreJcF+nQt7+I3KjsX2yDVx/nRrHnW6iPON27r2c9VyM5fIq51wgl9qI1u6nvPWlyP3Sibi7e+MOIfqkbjPf/qRPl8at7PL7oGKdkdy2UfxzwMeAczIrjzwB3A58G1gP7gJ9S1SdFRDBZMK/FFMp/u6reOvKZBwo0FshNsfpWca6GxFXToNssOL/Q/qjwwuAWvIBXFWT/8aCKd6qFQHa1PPv82U0Xg6a087p9UL5g1L3H9m8JFvLoGOd6zF0t5g8Bf6mqb7RNWY8Gfh34kqpeKyLXANcA78L0Bpy1PzuAj1If1xcYEdkdng8SY2E19nDzBKXJajsi4WjLGWVQZ0F3tGaLuGvrRuli+RYRJL77pMZyL8j73R3BfzxazPrwBAuziDwfeBnWYa2qc8CciFwOXGYP+yTwZYwwXw5cbzuZ3CQix1VWJQMjpoi68G5zG8W5g5j41muyzkYCyYAV7EX6qRdqNZfcGY5B1qy3Pb13YU1IB83Vdz+kd+/pX3ytWYhdikXIQBNCNpqojCWhi8W8EXgU+F0ROR+4Bfgl4IVObG2n7FPs8XWZMKcBQZiXCMlzmLf1L9wtc96zytztb7pnL8nmja0RA07kCpFwglxtyOoYYQr3gl0aTaF/de4EaPfDD0ly5uZaf2/1oujvT9avgzwLDQCWmUn3MSfAhcC/VtWviciHMG6LJrpkwmCzZ64CWLt2bYdpBJpoclsUi4L+4hGUohRqUW3PQltCN8SC064HWcneMX5cdGffcuU9dO95o8uoBdfPMDntxcbP/e1gsxxpRlgrY0noIsz7gf2q+jX7/HMYYf6ec1GIyKnAI97xrZkwAKp6HXAdwPbt20Na0xLgLMM+V4ejyf1QV2+j7vESMkwoXbp3X+/4tuiMKi1xzsU8PF/wQkS4jfThb5O88JTBBwZGzyIrECw1A4VZVb8rIg+JyJmqejemM/a37M/bgGvt78/bl+wCrrZVlnYAB4J/eXkpEhSqljMM/+nsGoXQxKBIiRq6uDgKca4rguTTkgBTF4dMJI2CHELTJptpiMr418D/sBEZe4F3YAogfUZErgQeBFwF8hsxoXJ7MOFy7xjpjAOjYRhB7rKwNyhGukmQfd91S+z0SMS5bf4NYXejtpID44FOweIfqvrPQF0Rj1fUHKvAOxc5r8ASkN69xy5W5QtLyR6QNNF3XIMQOh+rT7Jx/cAY6S7+576kk4XgRNn6n52vPoj0dDHRrozAdJHevae8sDfIvzxou0/dWG6bE91cSdav6xPnvpZTLbHUg6zn9L77mxcvuywSuvPYzjCjjOIIjA+THpURmDKKWhN1HbOb/LNNglZXo6LueF+gO/iZi67eFRF2gr2gOhvDmki5LjjOOTDemLXw8RXm8XWyBJacojh9tfNJW3q2Xza0Wj7Up07A3bbIFDBK1q8bPMcaUa5ur32dn6xR/Xua5lZ5H0L683QzRKH8I06wmAP9ZTXXng5RbKzbpvrKjrrqbtCp/CdRvcuiSXR9K7o03xrfc60ro8sFJLBiCD7mwEThZ6Q1JqF06bI8iIZIjkFiXYpb9saqFeOuPuW6LikDitK/KvrJxn3DMI5hd6P628YVRcgnPSojMJ6M8gvd9EWsTbqoE+uFpmZ3PN73NZdEuckfXrX0m1wZtqVTqaiQ9/rQMWR6GWODOQhzYHhKnUC2zvZ6CC6EIYofDWpW2jdm07gVkU42bYA8Ly+Khqpu082YL/4FYQ4sChdK1ldofoiwtL7juobwdR23S2JJNbmkWlfDHhsqwE0RY2wyB2EOjIQiGcNZm7maxcNhalc4hsnYazu+aWGygdrwQUuwnKePYDEHVgyNjVMXw0KSW+rGGBQC2JLK3VcONTDRKJDnQZgDK5TaLLzFuCqGtZ6r+wa5NerGqjZurSEI9oShQLCYAysZX7RKHamHdXN0iY1eyOubjhtifosR7GkOTUtOP61+7SHLzO88h1zRLANVJElMd/MkWfI61SGOOTD2HInQO+gXqqFcHcuRDLLIhchk04YVa00na0/vvV+JkZq+//8pJ0MSkT3++JGe3sgW/0Tk1Zi+qDHwMVW9tuG4NwKfBV6qqrvbxgzCHFhWihC1Ufii2xiFeTRs+VNWpoujaD7rFoBFIIpqqwemjzy6DDMEkJEs/olIDHwEeBWmScjNIrJLVb9VOe5Y4BeBr/WP0k8Q5sARpaiPYb+kLmmkUaAHWahdBPdIibK/3yWoVP6eIyHUrlt6H8OGMDa93ifX8r6iWFUEcXujgWVnNBbzxcAeVd0LYBuEXI5pJOLzm8C/B36ty6CdhFlE9gFPAxmQqup2ETkB+DSwHtgH/JSqPikigjHrX4splP92Vb21y3kC00+13GfnmhZNdAyFWxSL7ADuc0SjO1ySTEdhrO3mYWwRAAAgAElEQVRwY9/f9O49ZcGvlnN13UC848caBR1NVEZd8+kd/gEicgFwhqr+mYiMTpgtP6Sqj3nPrwG+pKrXisg19vm7gNcAs/ZnB/DR6kQDAUeTQJUSO4aNax4VixHkAdbpUgh0IaxDCnIbjbWom+K9J6p2dWdhPklEfJ/wdbZnadMgxRsjIhHwQeDtw8xsMa6My4HL7ONPAl/GCPPlwPW2k8lNInKca9q6iHMFVhiltG+/jkWVYcLfujIopXux+OF39gI00gSWXEdW38M1C6hj7K3iQXT/1z6mqnUdnGBw8+ljgXOALxtnAi8CdonI69oWALsKswJ/JSIK/I69WrzQia3tlO3a/daZ9qcBQZgDC6LIKqwKdFvCyDAMin1ezNh11IThDerWXUcyu6l/Qc0T5WLMGss/3bN38EWhSZAnyipuYTTX3JuBWRHZADwMXAG8pTiF6gHgJPdcRL4M/NqoojJ+QFW/bcX3iyJyV8uxraa9N8GrgKsA1nD0yGM5x7GUYmBxVK3AgYkrS8Gw3bebXBotHWFKF6C2Poi59hbfvM4wtXcYNUWZBgmyVuYmqtMlyiOIylDVVESuBr6ACZf7hKreISLvB3ar6q6FjNu1Geu37e9HROSPMSuR33MuChE5FXjEHj7ItHdjXgdcB/B8OWGMQ70Dw7KUMdGlkp8L6fTd5ZiFCPwgwW4pO1raXhVsF3LmqGvL5fY3iXiuJmkj13oruzKnqiA7pkaULaO6jqvqjcCNlW3vaTj2si5jDhRmETkGiFT1afv4R4D3A7uAtwHX2t+fty/ZBVxtw0Z2AAeCfzkwKppaTbWy0CJKC+gA3nruYUTa3++Ls9fUtqBGrNt8zIU4N1jGVbI772ndP7FMeK2MFwJ/bB3XCfD7qvqXInIz8BkRuRJ4EHCmzY2YULk9mHC5d4x81oEVTbLuDJAl6D4xROLIwDG6hvG1iXTbOV0c+JALcPG2LcavGAFNZbS9fVMryoCM8X36QGG2gdPn12x/HHhFzXYF3jmS2QUCdThR9v2rHnU9A4tsNFhYgsWggvvD7Ovymqbu4/bxsG6FeNuW8oaqKFeuc9MsyID1MS/3JJoJmX+BiSJ58akwN2efxPZW36T7tnXOLqrcjSqqYqmpFv5Rl9qMWSCsZPepSCGmfSI8gKkX4VokVJcLBEZFXcUx52selLRRW4K0C8tZhqxJoKF3pyCCqBZirZgICuc7lsr83faVKcgewWIOBJaOJku5TqhLJUi7Lh62Mcq6E1AfdQF9rgZzbP0+J8RVQcYTahUh3rZlZYvzIlpVLjVjJ8ySzKDp/HJPIzDBVON4myzpUoRHnZtjmMSTKsMmvzQJMpQs46FfWz55MZ7YxcXCyl5pVvSI4piXirETZjA+svzY1Rw+YQ3ZURFH/XGnSnmBQH0hnig2+5xgu1ZS7nGVmsW2JWUYUW6bT9Pf5L+2RuT9cLl425bC0p62uOUq4xyVsQQxR4tDYvMl0jhCFKJ55am3XbrMswpMDA0il969pyfIIqT33W9+nNXcEgFRdNFu+qljFILelvXn5ux+8rz3O8shTSHLzd+4Z6/J8ovKc1aR2hhmt20krp5xRjv+LAPjZzFrjhyeIzqUIM+bQVKYeVZ5/Od3cuLvfGW5ZxeYJPxiQb4lnSvJmZt7McA1adG12zucp6BrWyrVZgHu6i5xv91FKc9MzYz9D5UO92OeS5Eb1VC5O1aIO2OMGT9hBpibJzo4TzSXo6sjokxZfWCMPfWB8ccXuSNxn1gVVVfLYpiiSINip6uCrNZqzsx3JXnxqUB9JMuK8SW3MM6ujPESZol6yQOqRPM5+UxENK/j7KcPLBPJ5o31VmWReGKfN1nEC2GQe6Ft/LiS2ecLdZOvu02Ii+15eUywTU7z5Q31G2eUsU7JHjsfM5EU1oVkOdF8bnzNqXLgrZfyvV/eudwzDCwTyaYNpTjk2upo7vMTVUSwKlB+Ccxq3HPdIlp1TB/P36xRVPhu635Kr/HHrPNX+/Nwoptl1jL2f7Qsyu51aYamKfGxxxIdfTTRqlVIMlP/N6xEgo95CPwPaGbEWeZzIiJWPZOTHOp2lRvnlvChJOkCsYJTShKpCq41Hv0iPoV/OSu7w5LZTc3FfpyP2P2upn9XRLQQXd/UGeR9c+N7JTtL5y+JbYv161vJ6sRbSZej8/QEEVwZHZFIkMi6M+IIUUUyJcqULFJAkBye/clLyGeEY3//q8s95cARxEVQNLow3HEVsU3v3tPYBSXZvLHeonaUfNND3Pr6ouwVBVKbpWeetMQ6O1H2BTm3URdNr7G4ztPxC14AWYZmZjEwn5vruT0CIfOvK5ormudIJGhsb/0UI86i6LySJ6CRILlxbbzgU0GcVyQtUQ/J7KZeWNy995Xjl5t61e3dN3yHbugPN8sxQhx5z/P+14gvsHklrtoX5cz6itO02K6al19j55o99VTpPNmBAwPnv6IJwjwEqjA3j4gQxTGqMToTEdkrfRQLmUZoBJLBM1dcUrz0eTfctFyzDhxB/LZIffiWo7Wc03vvaz6eBTREHVTH2Ilz3UubFvCq/mdRM4jkJvxNxMYoZ0gGGuWQGbeFuqJOgc6IBlfGcGQZHD4MeY4kMWQZUSxobD7pmtkvRWJ/50KeAGJEWjI45rNBoKce9W7rB8Qdl0S5LsrB0ljkyPczV9wZfj2KklDnlf3+uRrKlZaiNJxlH0c2MSTDWCMCsSJZhubzEMfmJ9declaW1b4PgQpjHJXRWZhFJAZ2Aw+r6o/Z5oM3ACcAtwJvVdU5EVkNXA9cBDwOvElV9w0zKZ1PTePAQ4eRmRmiKEJnTInHPBEkMx/amAiNFFdxK09AY3j6zZeYK2IKx3wuiPQ0kt53P8na061QxeWdLnPNCbJLMW4Q3mTThnqrua7ovWcNV7/W0uT/hfqMxKrQ9/mwvedxbOaQ5IWbQ45aA0D2gEkkccLcCzkN/uQ2xtliHiZc7peAO73nHwA+qKqzwJPAlXb7lcCTqroZ+KA9bjiK0KAc0gyZz5C5FJnLiOZzojQnyhRJFZlX4jkT8xzNK9G8ufsDyGeMSD/7k5fw7BsvaT9nYPLw++JBv5/WW9Rzwpved39pm9veWjK0Ol6W9z6f/o/DNUn1f6rjVY/1/6Zq+F0SoatidHWCrl4Fa1bBmtVGgL3zqlvoC3Rj0sPlROR04EeB/xv4VTF9pn6YXpvuTwLvAz4KXG4fA3wO+LCIiO1s0uFkUTno3ovTlDyHDMgjBJBIjVGRml+SU7KgNTJ3fxqbTc+86RIj2hos6aXkSIUqptZSbGvQWlcStIju2Li+nLZdI84DO3GX3Bc1H/G2yIuWxcXCLRJXanLEVitU0aNXoXEM+x+uvDhYygOZEh/zfwH+L+BY+/xE4ClVTe3z/cBp9vFpwEMAtrX3AXv8Y4NOIn5yQORiSCM0iiDqGfdiDRcpfHL2d6ombj8zAq2RoGLE2Qk1GKF+5k3WgnbXgEjI3Z1gbhJa/N+Sm+zDNbv+qeNbFjhSpHv39YksAKokG9c31mtuep0v0IVFXT3Ot86HySSsc4+AWchzHUowH8vq4qJ6+0HMdyIC/cELyJOI+G9v6T6PI8woY/dlVBX/JlmYReTHgEdU9RYRucxtrjlUO+zzx70KuApgDUfXnbgkxqVdmUvRNqcSBUURxJwoN49RxR2mkbGai44y9ouh1q8X5TmxCEf/YbCkJ5G+cLeON2iNoo4VaBFTnc23sn0awu8GUmd5q5aK3wugGRApROYuUd3Xy7k4XGhpPCKxWkHIiG4sROTVwIeAGPiYql5b2f8LmD6oGfAMcJWqfqttzC4W8w8ArxOR1wJrgOdjLOjjRCSxVvPpwLft8fuBM4D9IpIALwCeqA6qqtcB1wE8X07ofUqd705s/GaeIXmOumIZOdatYe7YJHKSbMXZjZRj/RuYwuDW4i0sZDEWirrbTzHjHrp8B5Irq/80WMYTiWpr7786amOYvfGS2U2mp+Dde/oFuq5c6EKozWBUJI7QXIrbRCGylnQOmbuzpL9biUd83rbu06j7OyLIbzM6El1wVvE4UARFfAR4FUb7bhaRXRXh/X1V/W17/OuA/wy8um3cLl2y3w282w56GfBrqvq/ishngTdiIjPeBnzevmSXff5Vu/9vuvqXi4WLQ6DxPJJlyEwCIkiaGAs6VnM7pxGR5CbZREyrHIn8D1bPWga1bmcXf2q/R3VWNJDHwqHLdwAQzeWs+oubu0w/cISpE9MuouwSTvwMwbZ+gJrEfdXYGgW6yrCCXU2CyUzCFeoSrjLzObahXkUkyN99vXlMd4xvkddhrfS+8EMV4vO2FVZ5dMFZZthJF+jRuDIuBvao6l4AEbkBs85WvDmq+n3v+GO6nHkxcczvAm4Qkd8CbgM+brd/HPiUiOzBWMpXDDOoZhmaq/kwZjk6FxHlCklsPjSrTBEWidQYxfZWTtCSSFcRgMgKeLGh58pQm62lYo1n549OhMM/frHxM6fKzBeCSI8Ldf372nzKxeusICdbZ3tRFmAErO6zMzeg1VlNfHPBQgTb91s7MXWujdyuqUSKuBWeW+9oHCretsXUnPHjrStuD8nsvrRmrnYeU9fNZHSLf8WammU/sKN6kIi8E/hVYBUmcKKVoYRZVb8MfNk+3ou5WlSPOQQsblleczSD+JST0YOHyA8eRNasRg4dRubnkXyN+cCsnoEkIl8VG3F2Qf1oyQIuyM0+YyFLcSyFpa191nPJ/RHB3Gteap8AsbDqz4PLY6lw9YSJItJq5IF/3LozTGH4fQ+Y5x3EGYA0Kwun93kpWdB2EbFUUEgi0n0PGDfI+nXNSSNNDHKB1Ikz1mARgQwUcyGJzzmzd7wNscu//q2SC8P9leUFRG9f7lnfeW6et9UQmQa6/2knichu7/l11hULHdfUVPUjwEdE5C3Av8V4FRoZv8w/n6PWIDMJMp+Cb7WkJkU1AjQVojQpBJpCcG3oHCbiolz31tYqKOJgQezls+cSMbu8u7m+L5FmU/yhHQNcgfdk3RlGfG1xq/S++0nWryuE2IXNFa/r6GN2ZUObQuJKERnQE93ciHKyfl1vLO+iMBRtFnXdPivQAPYbUHLfGfHOic/e0p8ZiVWRzHP52YglYtAk6c1JtRR1F51/VknQJ96NAcMI82Oqur1hn1tTc/jrbXXcgAkrbmWshbnuCxaftw2ZS80HK/E+Kaq9BT6LKNaKoPQhV7TfUsndh10hl5JQF743PPG2ER2HXncxGstQDWNHGee7EkqIFvHKNpPPCWljxt6w41eSTGCAwFqBdmJcGstFeSyWqiVfF55nDRSpliK1SzUqxuqtK0la+JIzyt8FkX7Rru6LQC4+F0lNvY7865Mn0sLIojJuBmZtJvTDGNftW/wDRGRWVZ0v6EeBgX6hsRbmOrJv3MnB1+8gSpWjvv0sGkWlegSCZw2YSPyizIDDhdLVojbdR6S3SOgiPHoDmF+5khzMyWPhuX95CYeOizjh46Ev4VKR7tlbtm5F2msqe0QXnQNp3k1ESndX3vlrDAUXatdnXY+SJpGuCbHzcd+L4qNbCLZvYXu/XQx1cSfZE3y1VnXh544w372BRafHlBH5mG2uxtXAFzDhcp9Q1TtE5P3AblXdBVwtIq8E5jFZ0q1uDJhAYQY46k+MdfrMT+xgzRPzRIez+tuS4vOnSF6J9WzwLfr7xQp7yd9cOVzmlEggPpyTPBdx8A07OHhCzNMbYN2/DSI9cvwsOEsXcc5vuZ347C3EZ28hu+Mekq2zxYJWtaZG6VzOrdHkHnF1Oaox1FasRy7UVd+zT1Ujo8pxzs3XVHjJHeN81e77EkUmVjqK4CYT/ZGcuRlNTJZuvG3LZPYQHJEnUlVvBG6sbHuP9/iXhh1zIoXZcfQfGYE++IYdrPp+SnQ4synXbdZwj+JYX6+j/lsc56suXCNuUcRZ1XbVPHkuI3kuY/UT8zz/gYi5H72YQ8fHHDpOOOXDQaRHgojn6+3dKTUVwi8xZ0zGZHYTZHm5TnMdg0S5bY6jyk6ro7pwWDd/kf4FycrrSqLsuS00MZm2+S23N07BNR/ocrcytozxEtFEC7PD+Xfn/8VLiQ+bbg2iXhhQbhdFmsr8qRb/pDq/k/iujJKI9xJaonntuUtETFGluYznPT3P8x6E/OUXkh4d8+yLEo7/3VDcvwvxOWeS3X63ebxtSy9KwBWJryzSVUk2b6yv3azaHLdcs+jmrF4n0NGFZ5PfekdtjedqH0Jf1JfEeoayQNc9rhPopguHiMmszXPk4nPN4f/0zdpDJ1qUGVm43JIwFcLsePbUhNUHIlZ9PwWbti1lg7j/w+y+iC5arilEyPnZvH0lwVZP1EV7ySyF30+ZeTrn+APz5C+/kMPHz5CtFp736ZAC3ogI8dlbzCKTH9rWNd3aiWSNi6I2OcVfBKw5h4sEyW+9g3jbFlLv9r1YmLS/6xrFLkikq4WRBoXjNYlzG7H3XfAW/yRT48645Hwkz02VxzQn++Zd3eY+7gRhPjIc93vGEk1fub23UKGeb1gBV0PDt4ztbV8prbXPX2f3V28D6/65SvHF7rsqixAdzjjqkZw8EbJXXARA/KXxLUCzbKS5ec9dAkiNWLZFZSRnbi61YCoWD+uE3d/W57vtPS9C5A4eLofsudA7L9HFva4pesMfrwi380LwamlLZqlSRFfU+JGht1hYGBG2KFIco7vLVrJdEp8efENqDJkqYXYkf21iwfOXXUA0n5PPRM1xoQ7B1NxoKnZerHwbcXYi3theyKEU1rh5bqxzyXLiVOGwtVAuPZ/s6IT5Y+KhKtiNczfwxZBs3gjz6dBiXKLi8mirhwGeFVtXvrNpnlUh9ZNQ3DHrzqARe1yy7gyKMrcOifqPHyTKbZ9H36hwPmXfSvYWx+XicxtdGFPDGF9pxkKYt1y0kS/uHl08rhOr6O9vM7/t9vyyC3uXfmtO9y/02fC7yPmn7Q7/O+JZzq0r3MWO3kNx/uyiFoK6mngkB+ZIDkD6iovI1sQrMquwsWu1SK17oO+1PtWLbNOiWZsQa15KYBlo0S4GkXox9uiLk67EGfdtq8Q4l0TZVm+sZgOKqgmJu+T8IgpjGgk+5jEh+vKtxeP85RdiTFctwumgtwjYu92rD6vTGmFurPDVENbU1MZ+5qnDJCLIReegsZCvTpB/uK3Lnzgx1JbPhL4FvTafbSNuca8iYIX7oGo5V6M8qnNdd0Yhzum+B8zrswH3wQMEdqH0vW9VP7IftWKft93V1aVoF+Vy0xx2nMfh41ez6i+nsEZMEObxI/q7W8lffqFZ3FCxIXFClJoyoi5Sw/mVSx9uPzKjwddc34SzYV9VuJ0lb0P34vkMpqzc4rBp012PTTZvLISqJMQ173EfkfSLs0SgeeGOSB94CNK0vH8Y7Dlafcodfcjp3n2lGOJSyGA1FM7f1pf16h74f7vCTEy+Kp5eUQ7CPJ5Ef3dr37b8ZRcA2K4loBIVAq0uvRvK1rMf+1z1Pav0Ij0iIKdI6e5FhDT4vytGWXS+KbdY+LdjmSixTk59Eel3vrt043uujL6kEe//1Zql54tzTYumPn9xbtPporjv2L5xAeKYdK9XbnRQjY0WazfZtAHm097fWnVbxA1CXMmUrX7+dCZGVyU89+Kjhio1MEkIwZUxUTi/tLGmMe2rwNzW+VdZ34DO/bC48pdAcnoWOJgaHGqt8pzeQqKfuOK/vsXSkwzic7cW1cDGPftqqUQ5Pner+f/keW/BryUeuZEGa7mN9MH9ZavXt3b98Vw0xXzafHzd+NU6Hv7f5jr8VITXZTTG27bY2uWV9ZI+7P4IiGMTlfG1b3BU68wmnyDME4izprNXXGSqz0UutjM31nGu9rIrpTocBfbLkhfvcFWwewuProuKs6wl0+JTo0ovUaYq3O5LasOc4nO3Tk+MaUfyl11A9Mycee9yrMVbc+Cg2GdPRJ0Vqz94AfH3DyNpWnvhK6re2TKg1XFqz1EjxF3cOnWRKHUJLoiQbJ3t3bHlNRd3H++OQuMEFSE6NAdnb+l97twcprAm87jSpeffGuDvgdX2+M+p6nttNaUbgBOAW4G3quqciKwGrgcuAh4H3qSq+5Zo/kuOiy9OX2mq/kkiSKo9gfYpLQ5WtheeDXubOePt9jIUi4zFXE1iBfT8JPa4Ukhf5UsXn7u1bz4uFCpvKag+iUQXnUP8xHPmfXIZgdAswm01kz0XRrJpg6kB/Q+3wYVnk329/k6kqHrXMVKjLp55EPnLLyRfFZEnpnbF6j/rReoUsdNOoO26hB/OCTS+HxrZMNIkKupeACixaeeW58U4UyfKMNnCDBwGflhVnxGRGeAfReQvMNX4P6iqN4jIbwNXYuqMXgk8qaqbReQK4APAm5Zo/kcMFxs99+qXIrEieYRkSuREwVIIsu+PxvM5R5XngM6UFxAlMyVMJRHb2zAymVeFG8QMIn7LoGKwulV6+/iS83vdx91LJjgcKr/ldnMhqhOhum11i3vQt823TgddzPrcI9WxrBgfunwHa+zjg2/Yweon50tRQn3sfIm5E1Mlms9Jvth8bEmgs8r/v/p58C5K1TsA2W5SsEWU7JbJWbtYEL2b0rGkS88/xXR2BZixP4ppj+Lqjn4SeB9GmC+3jwE+B3xYRKRr379RMMoaxdUEjmdPTYgyWP1URpRBlkVEmfnylOKToRcvjXc7adeKTJNYd5gVcXf3OSOIignhyzHin4mtzaGI5IVAm1v4ikBXfrsaIfFcWiQUZN+4cxRvz6LIL7sQFSH+24VnPTrXTXLm5t7GttTtqrU8RDJJHX29/2qs8WTrLMfc9xQZ5m9e89gc8bNz7QUzv/LP5iOx8yVFEa2Bc/EFunFBuXeRSs7cbMIIrTXsFrhd8aLogrPKYYzVriZV15oXdTQRFvYkCzOA7QR7C7AZ0xH2PuAp2yEbTBX/0+zjogeWrVV6ADgReGyE8142XAGih35jJ0c9BslzSjwPM89kRPNqftxCIZTE2Udyz2p1n327ii7udWBTZMX0M8y1cGeYPm4Zmkuvfi7UflH6wqAiTIcLkaJI0Kj4/k9fyvP/e69Ik+mkURMiKALffRpdlRTXqiNOg3W7EGoTP1zYnidSM49aGyeqW5iocMn5RIdN9uMw2cN9Lg7/PXcWtPN31yw2K7Z+NZTj9SNQL3y0EGnXUzAGctMWzi+rOq5MfEq2qmbAS0TkOOCPgbp+6DXxCn37CkTkKuAqgLVr13aa7Dhxxm/2yng++L6drH4qYeYZWPNkTjSfEx+27ohce76/wuKwL7Si3WdNV3uyifT6GNq46140h5rn1S+fP46f5avaOw/a9wVyVpQmEfmaVTy96Riy1cILrq+viOf6zTk/+Qk3PwounjYSGy3R7O8UV37Ttv8a9svsFjzTu/f0sgYH0ZRI4tW+WAglcbbz8CvkAWR33MPB1+/gqO8dKubfdAcTHZwv/of5yy8kfV7Slw3aVwfax7VPq8sOrHme33I7sv1cZPu55LZWhlx8rvnfVlwizrouBDv3/NqxuB71Y81EuzJ8VPUpEfkycAlwnIgk1mr2+1y5Hlj7RSQBXoDpll0d6zrgOoDt27eP8Vs0mLXv64n0w9fs5OjvCTMHTecUySFKlWjOWNEu/E6cSHu3f30CDX2XNCfOkmMWbDKFWE07rKYUZHrnqNtXuAFyt4BmFjf19rt53m76iM/diszbSR5O+w9w43fpieisLhs66ERak4jsjvpFt2TrbK/a3LOHegWD0qx+Ya+ruyKS7o1c6+ZVTWRRRQ7OFU/js7dAHBeNHqILz8Z5+IrGqb57wEaDiCqrDs6z6rtKSo0lXHcxahJkr/9lXUag7v4m0UXnIBedY8YtFqBd4fzya5w7Lb+tuXbzWKJMtitDRE4G5q0oHwW8ErOg97fAGzGRGW8DPm9fsss+/6rd/zdH0r+83Jx2bbkg/sPv3kmewFGPxsw8q6x5KiuK62tkrtqSGis6ypQi29AJt1D6AIn2bsFM41gtBip8ztW321rrtVRFyxdsPxzLH/PQfP+iUsMYQ+9zf3ualWN3oXc+15i3pmv1oMW4RjxBd+ftUiwpOv8sZC7t86v683BRHsxlpPf2Ljb+wqJsP7corenf9fgLdMmmDWXxr0meKT1vqi5XTN4c41v10fln9cTYEbv5VP3z9tcttxOdf9bk9f4bY1XqYjGfCnzS+pkj4DOq+mci8i3gBhH5LeA24OP2+I8DnxKRPRhL+YolmPfEcNr/Wxbq+z9wKasOCMd821jU8bySHDQujywRT6hN/Q7x3ATug1TKNBRBY4zFTG6L17lVRe392NtNaIlrbUlmKT13t8htwjBorEE44fd/F+evjGVjiJP16xa2mFcVLnvOZNMG2oonJWduhmcOFvNqTP1284/r/cpP/OxOjttziPiZ+aK8rEZRn9D1Ffev+R/4c/VTtNO79/TmXPiZwX2o4rO3mAOzrLe/CJ+LWi/u0UXntHY7GUeECXdlqOo3gAtqtu8FLq7ZfgiYzlqUI2DDu3q+2vuv3Uk0LyTPQTQHa55UohRWfT83kR6Hc1OK1M8OdIuGLqojA9R09dYosi6B3Fatq2BX9zX3fNSxL6IN4WSOhVjI/mvrjutiXffdAQyIrGgKi6ujzpqsLKa6Tig+8blbe3cO/nydQHYsTZq+cjvH3/Uc0aEMybLCHZV3jZoRaT5XcSHzRLqIyqD3mYLSugM4lxkl67vUcRuKms1jrG+tyCIjcpaSkPk3gFGG3kE5/G7DNe19AA+89VLieSU+1Es8kVSLrEHJFEmsUKsSpc4VIEaEc5d44YU5qZYFmsoNalWoHQu1RP3XNYnzIAa9pilhZBB9lnfzwljJrSMCcym6Kia7w9S9GGRd9536sguJ5nKSZ+aN4O3+Zq3AJWduNpXsmpuckt4AABqtSURBVPzIXVAt5lWk8LuMUf+wpvHiuveXya7ZPOk+5sDy8YJP1UdCHHz9DqK0ItKq5EUJU0ViIMWIshgXB74rw7o3UC1/5+0iXHGco+7L2UYRGVIRdz8uFrpZ4U1FeBaLtxDWeC63zXdJgMnMPJz2pUV3DUPMk6iI4Z5/1XbSy3ew5vO9gkHFAqvfUqthjrWp2fa4vr58ImhivvbV9Qjx/ielcqDe315Y2RmQ50TnbRuLmPiFMNGujMD44Vb1q+Qvv9B+aa1FnMQ9C9rzVWvJX+t9IVVx9amBXqGlISndIjtUy3WPB4ltnUAOEvOhJ9pBlOv21c3D1VJJc3QmNn7Xmbgvs/L7P32pjdLRost7fCgjeXoe2XEe8YGDxkI+ONd/cer6d0dS+JSrZN+4swjRM5+N3nqBzJmLQHbcMeU7LOeLLvzMWt4/qYxo6iLyauBDGOfPx1T12sr+XwV+FmMqPQr8jKq2xmUGYZ4iXOEl/V9e4vlnbfQGoEmNGOZ4Lg8pi7Sq+ahVX9OEjddWrxu5VEXMF2jouROqlnQdoxLkQWO1ukEqEQ9QxPOq/zeowk1fJ/uhi0qZjX7yDdgoiGfnTKNTZyFX5yjSV9LTjwCJt5mFO/9/VwqB9P6u9O49ZN+8y8Sfu0Vll1loC0DFT5okGOe7duOX3oVJFmTLKCxmGxTxEeBVmFDhm0Vkl6r6K7e3AdtV9TkR+d+Af8+AMhVBmKcQ+cd/Ng8uPd88T3NTlCaKTPB/HJnIjlQR8p77ItKy37nOJ1zd5j+OKr9z6Ov0UrU8m1wdS0VdqFnUL7aNL685xv1tfhy6PDdHtHWW5Dvfh62zxbFOUBO7TdMUOaxlQS4NLqU7jfReYwkXYhnRGO/tiM/bBpnpcO0SiPyQuOJC6twXeXnxUmzvxc69FieF0VxbLgb22GAIROQGTFmKQphV9W+9428CfnrQoEGYp5mvmtvounUOAeZ/ZDvxXG4iAlxyinVtqOtQPegWvotQ5/2CVuvuiD2LcyHURWf0MSAKY0gGWo6V/YUla1tTSTpgUc9PCEmikiDjqsO14I7XJDGJSKqm+7irXR2JLQGgvVrWFYZq6zUpKKNKyS5KUFj2Aztajr8S+ItBgwZhXuFoLMbFkcelBUFiMbe6tgh/rZB2OoH295Sj7O6o0rhnIdEZg8hbludrxqu1/OuozrUpqmWQS6WapZeYN7NL3e1k66yxdL3Fufi8beauKYp6qdYu6sPOeeos4xqEoVwZJ4mInwN7nc1cdkNVqR1ZRH4a2A68fNAJgzCvYGb+qibfGsh/8AKziBVRVLgjs7fa1b6FMFTYVkFLDZ+mRcfSWVp93TWWc50wNo3h/p62EMFhe60OE1vtjq+8r8N2qKmrOyLPHe6NG/m/mxcMp5bud2aPqer2hn2uBIXDL09RICKvBP4N8HJVPTzohEGYjzCjjoseFX58td+RO7vsQnvbl5QiO8Te9vaK1LsXVxb+/AWxLgLuviw1LfT8Fl19C5U+1a7Xgwro1y1QFn9EZc6VDtRD0SbmDWMVPvrc8w9b0rvuLTL2fD+z2+b+V75A+24Jl4odn3Nm7bnH9bM6KkYULnczMGsbhzyMyXR+i3+AiFwA/A7walV9pMugQZgDrVSLuc+95qW2Ua2aCKo0N4ku88ZvWdweU9XBBkGrRjLYyI5SOy2HL+6xv8tfQPQnX4nV9ffViXrTQmdTbHMNtb70Aa9pei0MtpKzO+4hPnsL8bYtvSa9Hfs/uljrUZd+nQhaPFhDDWNKG18NfAFjTnxCVe8QkfcDu1V1F/AfgOcBn7Uleh9U1de1jRuEOTAUq/6i3Mr+4Ot3gEA079X0iITk2dR0eLGxsYUo2sI5LjJEbaKFxgKCKfCk2Gp8ZjyX0SjzdkHSLUzWiXZESUzV7YNe3WA3n2rSiD9WdVzMvIu+gnZ7Y7ace3nlDqLvePt+FOPTzX9conLByca8DvK4MKp6zKp6I3BjZdt7vMevHHbMIMyBRdGU7AJw4F9digpEmRFXjUzj2jyG9GghTyBbBaf+p6/wyNU7kRxWHzCp5aue6UUNxM8Z/3aU5Cb0z4lpmtdHRTQljjgrOze1RYpIBOhmQddEqbhn1dZhfee3Iq5JhEauxx5F5MxCKAoP+Sn32RhXfx8zJr5QfiCwEJqK69dxyoeb64Y89xM7iDKID2bEc0KcK5qBJFEvizGXfiu6TqBVQbQs0E7MBlnQbl+mplKcn7psQ9eM5S+lmPF8VcShE1dx9B/e1Pn96ES11ZP1+yennIweniM7cGC055smlIWHZR4BgjAHxp7Dx0XEc5DEkM9HtphTDvN5LxFGtFgcLFLPG3zDhZg5gZbYdnZpKBbUgIoU4Wtqf2Ot4fljVxUNfCPg6IX/+Y04P7Jf3pMogucfi4iQvOjklRdpMQShVkYgsAiO/4SxvA++YQeSKvmqCIls1enMLjh6C3VaE9HhY8LxvKgOV+TJF2joNbktvdhzZxT+ckFXJejXvlEc1vWLVepGomraU9mswM5ttqp3CV7d5+TMzUGcmwjCHAgsnijVXgSI7woWKYfm+YtyLsKjagn7UR1OkKsCDRSFfhyFfzsz7bdmIP9GOQIi2bi+mEtbskbi90Z021whe/e4pjZGlXTP3rLAu0VY+1rXC3ElJI50ZcgEkyNOl9ZSZwDXAy/CrP9ep6ofEpETgE8D64F9wE+p6pNi4kE+BLwWeA54u6reWjd2YHwYZcyqHxM9SuJDptGtX+Bcowjx3Bg+mhjTOV9Fr5WXi8bAi8XO1PqqTdSF5DmuG7SJAnED9keCyOG5vj6B1Z6BrlZzgV/ZDeor+FUKGCVnbm5talu7rfC/W6FfezpkOenDffkPKw/ViS+UnwL/h6reKiLHAreIyBeBtwNfUtVrReQa4BrgXcBrgFn7swP4KO2544FAK9kPXUR8KGX1d56ujbktvl62aFOxMVcQEJu0KLmJCgHbzDaWcpuuXE0NiSzqWZ4uPd0ttNmCPiVyNW2toLbLdp2lmmzeWL/QCD1BLQS84Y1pij4pWftqAputlU+WkZx8EumjjzUMuoIYX13u1FrqO8B37OOnReROTOGOy4HL7GGfBL6MEebLgettA9abROQ4ETnVjhMIDM2qR55GDqf1Pl+futCznS8piV9Rec89rxxezTuQ7eci81JY0UJifNBZVrFgzdyS9etqxbmKy8Br7Xhtfc5NuNeW+vz5PQF9crVzztG5eeJjjyV7+umB85xmJtqV4SMi6zH9/74GvNCJrap+R0ROsYfVVVs6DSvu3lhXAVcBrF27dgFTD6wUnCi3iZRP/vILTReXNC+sZgC+8s+tr6vDtXwquRJaIje6iHLzyXpuh2HEvbTNWufJxvVereuIdP9DfceuaNwd1ZjSWZhF5HnAHwK/rKrfl+aMp07Vlmx1pusAtm/fPr7vUGDZWUjZSdNMVMxiYZYP1Zsu8eonm8GUUknM2qSWCDQnWXdG/3bqhTbZtMHzMy/sK+BbyL7LJN27j+TFp5J+O9yoNjLGqtNJmEVkBiPK/0NV/8hu/p5zUYjIqYArztGp2lIgsFS4Ti7sfEmR/i07ziuFs9XhivmoV1ujiM6w3T0GulOqaP3xycb15cQQyzAWdyHKLorD+rlNsYycdP/Dw811hTHRrgwbZfFx4E5V/c/erl3A24Br7e/Pe9uvtpX8dwAHgn85sCx85Z+NILsKbdvPRXcPsJy11/uu8APb4vZ+Y9OSCFpruRYZUB+0o6Xc54v25ttHNXokUMukR2X8APBW4Jsi4px0v44R5M+IyJXAg4CLkboREyq3BxMu946Rzjgw9oxL6J1sP9cULoq0kytDDs4Vj0uRFGlWep6sO6MsxC2iXOvC2Li+98RvG9VgLSebNrS7UQrfdx7qZXSluso7ZnSJyvhH6v3GAK+oOV6Bdy5yXoHAohloHVeohrXFZ28hu+Oe/gW1JiEG0gfqF9lKYlxHJKVjGhc660LsrLine8MCX1dMgsn4KvOwfRgCganGrztR2+S0RZSHId27rz0Uzhdp/4LRICZdI1YCHnnHn2UgpGQHApZk80bI875Ftb5OKA0sxFp2glp3jJ9RmN53fy+DsDKvRYXorWCCxRwITBLOT5tr58W5rqI8yFJue31trYsRWfArDh3iZxkIFnMgYCmiMDau7+8bCEOJYDUDb+gCQl6eQKkWR6n2chDlhTP5tTICgRVBrcvB727twuJqwuNcYomznIcVYt+KLmKcm2pp+Izx7fjYM8bvXRDmwMhZqupyS0URI+ySR/qyWr0vsItLdr81b3RjLJQ6v3PxuKbCXHL6aSGZZFg0tJYKBMaawoXR4H4okkmgZyl7Al3aD6X6ysVYQ0ZNDAyvc/Nw7aROe3Eo5zksY2wxh8W/wNiRrDuD5LQX92+f3VRuozTKc/op0m6bFep03wPFT4Hm9T7eqigPsYA4iJ6fuSzK5DmkKcnJJ43kPCuGsPgXCAymsFhF6t0DNZboUOM7lwU1hZGiqE+Y6/zE/rz6ChZB0dm79JoRhbMVlrnzcbv6HdJbqIxPOL4QbJ2fhywnP3xoJOefNmTYuidN44i8GtMcJAY+pqrXVva/DPgvwHnAFar6uUFjBmEOjAXJ7Kbegpft2JHevaeoU+FYaP+6UlqzSEmkgdK+TuNV3ReOSjTHsKLc6sIoCX5sLwLlQ2QB51yRKCNJHhGRGPgI8CpMAbebRWSXqn7LO+xBTGORX+s6bhDmwPhgRdlRFA9qEeM+gR2ChZQT7aPOdeFtd0LbNctvMecO2X/dEXRUCSYXA3tUdS+ALd52OVAIs6rus/s6XwqCMAeOKFULuIlBlvEgQV6I6HYRUbDW8iC3ygCBHiTGQWSPAKMR5rrGIItupReE+QgzaaFko6Lrol0n69hzN/ilOBdLFzHsUoxoIa8PQrwMdBfmk0Rkt/f8OtvoAzo2BhmWIMyBJaOvOzT0RSg0CetSi/BCaPQrL4JRC/IoS65ONcP5mB9T1e0N+5akMUgQ5sCicIV/iKL+BqN+12a3sNd1zHES5IX6gBsI1vF4MKKojJuBWRHZADwMXAG8ZbGDdulg8gngx4BHVPUcu+0E4NPAemAf8FOq+qTtdvIhTKH854C3q+qti51kYPno5IKI48Zjq6JaPaZp/OUW48USxHfcGU2XF1VNReRq4AuYcLlPqOodIvJ+YLeq7hKRlwJ/DBwP/LiI/DtVPbtt3C4W8+8BHwau97ZdA3xJVa8VkWvs83cBrwFm7c8O4KOMwBEeWHoaF9NaLN2FiOc4CW7Xxb4uBCGeMJSRZf6p6o2Yzk3+tvd4j2/GuDg606WDyd+LyPrK5suBy+zjTwJfxgjz5cD1tovJTSJynGvYOsykAktLsnljs+CKjJV4Dkuyfl1f/Yoig8+lWFebmC7CVREEeYKZwloZL3Ria7tkn2K314WOnAYEYV5mqhbx0GUoJ4CewMatdYxL/fv8+hhjXDshMHrGuVD+qBf/OoeOiMhVwFUAa9euHfE0VjZNbomRJFSMIZ165TW9boy/nIElZoz/9wsV5u85F4WInAo8Yrd3Dh2xcYDXAWzfvn1836Expy61eFoFuIl0776h3BF9raN82goOFU1P93WfXGA8GfNu4gsV5l3A24Br7e/Pe9uvtmmJO4ADwb88Wuqs4ZUmxFU6Z+y5ehltgqw5pJk5ZiZprbUcYoYnnEm2mEXkDzALfSeJyH7gvRhB/oyIXIkp0OHS2W7EhMrtwYTLvWMJ5ryi6PMNr3ARdnRxX/gJIem+B5q/iKGp6cpkkoVZVd/csOsVNccq8M7FTmqlU0RNgHFNTOFCXZVGV4TrbecK07e4E5L160ri6h4nG9d3cnWUeusFphtlZHWyl4KQ+TeO5Hmprm81tXkarebFCuKC0qWrX8xc+63swJQy3s1sgzCPIenefSSnvoj0O9/t25ds2kCyaUOtFZ2cubm2UHtBTXPP9O49PXeJ8786wcoz0gf3L+pvWWqqlnCnRUDfnyyVgsZj/GUNjBBlKhf/AktMnShDe/xx1yLyydbZ8usm1AJv7NHnF8X38SzkOms4WXt678uqGpqcTjuT7GMOTB/pXfcu9xRGQjWTr7CWq1+4jr7Ecb87CIyYIMyBScP5WsfZz1qXzVhyZTQI8jj/TYEjxWiKGC0VQZgDfVSjG8aBTpXrNq5vtY7H7W8KLCOKWWQfU4IwBwp8a3NcQsd8QR5YXGmMw58CY0iwmANHkrrIhGHSlsdCkD03xaiq3SXrzmjN5AusJKYzJTswxixGWAemNQ9wF5TGGsJ1UFjG1ooZJlKkNoY5kpDRF2hGQcc4NDII84TTxeXQZikvqAHpUoiejaEeJsuxManExXEP6mQdWNmMsesrCPOEslBB7mpNt7o9ch2JIBfnyOsr4lXnUJ17kXLtBDoIcWAYgo95cnlV9JODDzqCJBvXmy4jzrLs2OAUFmgdD/n6xnFPP8088MUzy2p9vsO2fCoJdCTlJrCBQB2qISoj0E4yu6ksIjWC4t/i97VI6kJHAXdiWHITuJodnmAuNM7ZJXG411fHqRPl4mIxxEXIHV91jYRSnYGCMb5wB2E+gsTbthTtbPzsuy5RB9VQNgCyvCSMfT5X50OL+gWq8Tx1flvfFxdJyY1SsrBrhLCgWqNj3wM9C7c6Dr07g1qq2wdcyAKBfhTNsuWeRCNBmI8g2Z33DDymzWIceGtfKYk5KC65btzaGhItldtKFrb2PujJi08tH/dwfyObQpyhVFBoof7rcYm9DkwAK7Hsp4i8GvgQEAMfU9Vrl+I8k06fK8JZ0zXC2iY4zvJc0MKeK3VpQ4fqfL4lAW36MKdp+em3uzWuSfc94BUPykDEPBcBiYYS6SDKgaEYUbjcIL0TkdXA9cBFwOPAm1R1X9uYIxdmEYmBjwCvwvQAvFlEdqnqt0Z9rkmn6XZ7UDRC6dh1ZzSLV6XMZ6euH+u8lo3VkphQfJiTtaf3zkG9oHclfXA/yYknopqTPfHkgscJBLqigI7AYu6od1cC/7O9swuxqori+O9PfkVKZlaIVo4llETYYCIUPlSU+jIFQvOUD4FQCfUQpAhiQQ8GFUSRTGSaRVpWNA+GWRo9pVmNOmLqmEKWOGRq9pA1M6uHva+embn3zh3v+bg06weHs885e87633XPWbPP3nude8bMbpXUCqwBHq123ixazHOBLjP7OQrfBLQAHphrYNT0m6G3t+bWYk0BMfErIEOesxSMrQ/6qvTBSalm0fWcPp3auRxnSCy1F+XXEu9agNWxvAV4XZLiLz6VJYvAPBVI3rEnCD/M6tRCT09q7wDuOXosBPq+3tgF0ceom6Zdmhkxbergub+xbklDsr7j/J9IafCvlnh3sY6Z9Ug6B1wL/F7ppFkE5nJD6YP+M0haCiyNmxckdWagZbhMpoqzcuGXlDVUm5xQvcEbdBSfzZyqPzTcKXcZaKgD1xG4jN8R6895zmz70rZMrrH6OEl7EtttZtYWy7XEu5piYpIsAvMJINFRyTRg0JB8/GBtAJL2mNmcDLQMi0bQ0QgaXEfjaXAd6WJmC1I6VS3xrlTnhKRRwNXAH9VOWmZ0p26+A2ZKapI0BmgF2jOw4ziOUzS1xLt2YEksLwZ2VOtfhgxazLEPZRmwjTB9ZJ2ZHUjbjuM4TtFUineSXgD2mFk78DawUVIXoaXcOtR5M5nHbGZbga3D+JO2oavkQiPoaAQN4DqSNIIGcB0NSbl4Z2arEuW/gWG9dEdDtKgdx3GcnMmij9lxHMepg8IDs6QFkg5J6pK0PEe7xyXtl9RRmgojaZKk7ZKOxPU1GdhdJ6k7OT2wkl0FXou+2SepOWMdqyX9Gn3SIWlR4tiKqOOQpIdS0nCjpJ2SDko6IOnpuD9Xf1TRkZs/JI2TtFvS3qjh+bi/SdKu6IvNcYAJSWPjdlc8Pr1eDUPoWC/pWMIXs+P+zK7REY2ZFbYQOsuPAjOAMcBeYFZOto8DkwfsewlYHsvLgTUZ2J0PNAOdQ9kFFgGfE+ZBzgN2ZaxjNfBsmbqz4nczFmiK39kVKWiYAjTH8gTgcLSVqz+q6MjNH/EzjY/l0cCu+Bk/BFrj/rXAE7H8JLA2lluBzSn5opKO9cDiMvUzu0ZH8lJ0i/liOqOZ/QOU0hmLogXYEMsbgIfTNmBm3zB4DmMluy3Auxb4FpgoaQopUEFHJVqATWZ2wcyOAV2E765eDSfN7IdYPg8cJGRJ5eqPKjoqkbo/4mf6K26OjosB9xHSeGGwL0o+2gLcL11e9kyNOiqR2TU6kik6MJdLZ6x2Q6SJAV9I+l4hCxHgBjM7CeFmBa7PSUslu0X4Z1l8JF2X6MrJXEd8FL+L0EIrzB8DdECO/pB0haQOoBvYTmiJnzWz0qv7knb6pfkCpTTfuhmow8xKvngx+uJVhTem9dNRRqNzmRQdmIedqpgi95hZM7AQeErS/JzsDoe8/fMmcAswGzgJvJyHDknjgY+BZ8zsz2pVc9aRqz/MrNfMZhOyx+YCt1exk5kvBuqQdAewArgNuBuYBDyXtY6RTNGBuab07Swws9/iuhv4lHAjnCo9hsV1dx5aqtjN1T9mdirelH3AW1x6PM9Mh6TRhGD4vpl9Enfn7o9yOorwR7R7Fvia0Gc7USGNd6CdixpUY5pvHToWxO4eM7MLwDvk5IuRStGBuZD0bUlXSZpQKgMPAp30T51cAnyWtZZIJbvtwGNx5HsecK70iJ8FA/oGHyH4pKSjNc4EaAJmArtTsCdCVtRBM3slcShXf1TSkac/JF0naWIsXwk8QOjr3klI44XBvhhWmm8dOn5K/KMUoZ876YvcrtERQ9Gjj4RR3cOE/rSVOdmcQRhV3wscKNkl9NF9BRyJ60kZ2P6A8Fj8L6G18Xglu4THxDeib/YDczLWsTHa2Ue44aYk6q+MOg4BC1PScC/hsXcf0BGXRXn7o4qO3PwB3An8GG11AqsS1+puwgDjR8DYuH9c3O6Kx2ek5ItKOnZEX3QC73Fp5kZm1+hIXjzzz3Ecp8EouivDcRzHGYAHZsdxnAbDA7PjOE6D4YHZcRynwfDA7DiO02B4YHYcx2kwPDA7juM0GB6YHcdxGoz/AOJq5TpS/OwjAAAAAElFTkSuQmCC\n", "text/plain": [ "<Figure size 432x288 with 2 Axes>" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(np.mean(f.variables['flagellates'][:,0,:,:],0))\n", "plt.colorbar()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true }, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [conda env:py37]", "language": "python", "name": "conda-env-py37-py" }, "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.2" } }, "nbformat": 4, "nbformat_minor": 1 }