{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# GTSPP Data" ] }, { "cell_type": "markdown", "metadata": { "colab_type": "text", "id": "cTOnSg9Sumfc" }, "source": [ "This notebook graphs temperature and salinity against time, as recorded in a dataset in the RIDDC ERDDAP server. I believe this is at one specific location (can be confirmed by checking the map at the ERDDAP server)." ] }, { "cell_type": "code", "execution_count": null, "metadata": { "colab": {}, "colab_type": "code", "id": "sBojlXO9ZZGu" }, "outputs": [], "source": [ "!pip install eofs\n", "!pip install netCDF4\n", "import requests\n", "from netCDF4 import Dataset as NetCDFFile\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "#To use cartopy in a Colab notebook, I need a particular sequence of installs (as shown at this notebook: https://colab.research.google.com/github/adamlamee/CODINGinK12/blob/master/notebooks/quakes.ipynb#scrollTo=3LkZkXvnMAr4)\n", "!apt-get -qq install python-cartopy python3-cartopy;\n", "!pip uninstall -y shapely;\n", "!pip install shapely --no-binary shapely;\n", "import cartopy.crs as ccrs\n", "import cartopy.feature as cfeature\n", "from cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatter\n", "from matplotlib.axes import Axes\n", "from cartopy.mpl.geoaxes import GeoAxes\n", "GeoAxes._pcolormesh_patched = Axes.pcolormesh\n", "import calendar\n", "import datetime" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 606 }, "colab_type": "code", "id": "vt0D0queFqUW", "outputId": "117f68f4-2318-4acb-8ad5-786b273d6952" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "[6. 6. 6. ... 5. 5. 5.]\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEWCAYAAABhffzLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO2dd5gUVdbG3zORnIeMjGQwgAiIKAKCAcy6xjWurmHV1V3DZ1zZVXddV/1Maxb1E8TVNUdQVIISBEGyxCFnyXEYzvdHV9PVNRVuVd1b1T19f8/Tz3SovnVruurUveee8x5iZmg0Go0md8iLuwMajUajiRZt+DUajSbH0IZfo9Focgxt+DUajSbH0IZfo9Focgxt+DUajSbH0IZfo8lwiOgeInol7n5oqg6k4/g1cUNEO0wvawDYC6DCeH0dM4+IvlfxQET9AQxn5pZx90VTdSmIuwMaDTPXSj4nojIA1zDz1/H1KBqIqICZ98fdD03uoV09moyFiPKI6C4iWkxEm4joHSJqYHxWSkRMRFcR0Qoi2kxE1xNRTyKaSURbiOhZU1tXEtH3RPQsEW0lovlENNBl352J6DujnTlEdKbx/jFEtJaI8k3bnkNEM330+WoiWg7gG8s+awL4AkBzItphPJoT0VAiGh7kuI3v/I6I5hnbjiKi1iF/Gk2Wow2/JpO5GcDZAPoBaA5gM4B/W7Y5BkB7ABcCeBLAvQAGATgMwAVE1M+y7WIAjQA8AOD9pFE2Q0SFAD4BMBpAY6MfI4ioIzNPBrATwImmr1wC4C0ffe4HoDOAU8xvMvNOAIMBrGbmWsZjtcP/Rui4iegsAPcAOBdACYDxAEY6tKnJFZhZP/QjYx4AygAMMp7PAzDQ9FkzAOVIuChLATCAFqbPNwG40PT6PQC3Gs+vBLAaxrqW8d4UAJfZ9KEvgLUA8kzvjQQw1Hj+EIBhxvPaSNwIWvvocxuX4+8PYKXlvaFI+P0R4Li/AHC16bM8ALuS/dWP3HzoEb8mk2kN4APDfbEFCaNaAaCJaZt1pue7bV7XMr1exczmaIZlSIzKrTQHsIKZD1i2bWE8fwvAuURUjMRI+idmXuajzyucDtgHosfdGsBTpv78CoBMx6LJQbTh12QyKwAMZuZ6pkc1Zl4VsL0WRESm14cgMQuwshpAKyLKs2y7CgCYeS4SN4LBSHfziPbZLZROdpjdCiQio8z9qc7MP0jejyaL0IZfk8m8AODh5GIkEZUYPuugNAbwRyIqJKLzkfCzf26z3WQk3CF3Gtv2B3AGgLdN27wF4BYAJwB4V2Kf1wFoSER1fXzHjRcA3E1Ehxn9qWscuyaH0YZfk8k8BeBjAKOJaDuASUgsagZlMhILohsBPAzgN8y8yboRM+9DwtAPNrZ9DsDlzDzftNlIJBZpv2HmjbL6bOxjJIAlhnvGzhUlDDN/AOCfAN4mom0AZiNxXJocRidwaXICIroSifyA4+Pui0YTN3rEr9FoNDmGNvwajUaTY2hXj0aj0eQYesSv0Wg0OUZWiLQ1atSIS0tL4+6GRqPRZBXTpk3byMwl1vezwvCXlpZi6tSpcXdDo9FosgoiWmb3vnb1aDQaTY6hDb9Go9HkGNrwazQaTY6hDb9Go9HkGNrwazQaTY6hDb9Go9HkGNrwazQaTY6hDX8OsmTDDvyweKP3hhqNpkqSFQlcGrmc+PhYAEDZI6fF3BONRhMHesSv0Wg0OYY2/BqNRpNjaMOvkc55z/+A81/Qtbw1mkxF+/g10pm2bHPcXdBoPJmwcCMqmNGvQyXxyiqPHvFrso6HPp2La97Qaq1xs21POXbu3R93NwJz6auTccWwKThwIPeKUekRvybreGXC0ri7oAFw5NDRALI/OmxfxQFUy8uPuxuRomzET0StiOhbIppLRHOI6Bbj/aFEtIqIZhiPIar6oNFoNF58vyj3clpUjvj3A7iNmX8iotoAphHRV8Zn/8vMjynct6YKsmzTTuwpPxB3NzRVjBGTl2Ng5yZxdyNSlBl+Zl4DYI3xfDsRzQPQQtX+NFWffv/6Lu4uaKoI05enAhC+mb8+xp7EQySLu0RUCuAoAJONt24ioplENIyI6jt851oimkpEUzds2BBFN11Zu3UP1m3bE3c3NBphznv+B5zxzIS4u5GRbN61L+4uxIpyw09EtQC8B+BWZt4G4HkAbQF0Q2JG8Ljd95j5JWbuwcw9SkriD7fq/Y8xOObvY+LuhkYjzLRlmzFr1da4u6HJQJQafiIqRMLoj2Dm9wGAmdcxcwUzHwDwMoBeKvsgg1fGL4m7CxqNL5hzL0TRD2Ubd8XdhVhRGdVDAF4FMI+ZnzC938y02TkAZqvqgywe+mxe3F3QaDQSORDBjXHt1j34dOZq5fsJgsqonuMAXAZgFhHNMN67B8DFRNQNAAMoA3Cdwj5oDL77ZT2GfV+GN67qqWwfd78/E4c1r6usfY04e/fr6Cc3Nu1M9/Hv3V+B4gK5sfwXvTQRZZt24eQuTVFUkFm5sp6Gn4gaI2HEmwPYjcQIfarhqnGEmScAIJuPPg/QT01IrnztRwBAeYW6kc7IKSsArFDWvkacOauj8+2v3LwLLevXiGx/Mvjgp1Vprzve9yV+eehUqcZ/xebdAKKZXfjF8TZERAOIaBSAzwAMBtAMQBcA9yExiv8rEdWJppsaWZzy5Li4u6BRzKotu3He8xMj29/x//w2sn3JosLGGO/ZJ2+WtGvfflQYUhCfzVwjrV1ZuI34hwD4PTMvt35ARAUATgdwEhKLtxpJ3PbOz3h/+kos/Ye8NPitu8sPPl+6cae0djWZx/Y95TjukW/S3hs+aRku7d1a2j527t2f9bHvFTb6PAx5I/O3p6Rmvht27JXWriwcDT8z3+Hy2X4AHyrpUY7z3k8rpbe5yphyWmFmJNbgNVUFu9HlfR/Olmr4rx8+DeMXZrfMwf6KyqN7mR4Zs3snAz093lE9RHQLEdWhBK8S0U9EdHIUncs1du2LVulQpb8/CuxGbWHYuGMvdu3bjxW/Zm+o313vz1K+jx/LflW+D9XYnToyz6b//Jga8f/zy/kYMXmZxNbDI7LU/Dsj8epkAPWRiNR5RGmvcpTligzO+IX2mc9Tlv6aNRfxDhv5325/HS11Hz0e+hpd/jIKfR/9FhMXb5LadlXCTi9p255ymy0zF7vzSWbuw8L1O9Je3/tBZkWtixj+pC9gCIA3mXkO7KN1NBnKS+PsE9AufXUyzn9BziLgPsXhg/fYjGS3K9SCX7h+u7K2qyIyDdv+igP4zfM/YIIid9LfPplr+352z3/9IWL4pxHRaCQM/yhDaVMHCStAlS/QGrOsglFz1iptP+qqXrLdSHbs3LsfY+atU76fKNgkcQHz1537MHXZZvzpnRneGwdg2Pf29RxyqSCLSALX1Ujo6ixh5l1E1BDAVWq7lTkssoz8yjbuRGmjmsr3u3V3OepWL1S+H1ncPHK60vZ3l1cobd/KXz+Zi6uOO1R6u4vWb8cr45eCCPh5xVbMXbMN39zWD21KaknfV5TIjFVPBhwwM54Y/Qt6t22IPm0bSWvfiYL8zEqyUomI4T/e+HtkLkaAXPTSpLTX/R/7LpKKQ2u37skqw6+aXyOYtUTBdW9Ow+IN6SG1u/ZFe1NTgczBcn5ews5UHGA8/c0iPP3Nokiuufwcsm8iht8c1lkNCVG1aQBOVNKjDGPjjugMjnnQtHjDDnRsWjuyfWcr2/aUo0618DfIqETNrEYfALLJ3nw0Y5X9BzINP6UMf5TIjOPPdDznNsx8hulxEoDDAUTrcM0Rhjw9/uDzTIz9jYtym5jrJP/68pcIe6K55W17v7tMo5lnWKVte6INb86lay6IU2slgM6yO6JJJ5dGH144RWEA8mLK47zoqQoEyS2yhC+GIS6Xci5dcSIibc8g9T/JQ2Kh9yeVndJk1+hjj+KF13EOeQgAMH+tnLDLv31a+eayY+9+1CjMR16eHEN0zRs/SmknE9m8S14c/49L48kt2bs/+9daRBEZ8U9Fwqc/DcBEAP/DzJcq7VUOMtoSDplFdl+5+mAUN8HXfyir9N7hD4zCU2MWStvH1/Ps9W2iGOBu2K5eL0bWOslVr8dzgzz2H994b1RFEPHxv2F6jGDm76PoWC6xdONOXPvmtLT3ZF1EizfIm4I74dTVqlAF6u0fK2kUZiU9H/5a+T4+/jkzi45oKuMmy/yO8XeWURg9+ZhFRDOj62LVx06jZ8aKLVLa3h1BqKDTiN+sChqGONc71m1TP1KWNeKX9f8OyvPfLY51/5mC04Dnw+kOEVEx4Objv8X4e3oUHcll7Bb3xi5w9mtnGu9MtVcU3bB9L+rVKArdfhWYOLiug/yydjs6NQ1f2sJJhTUqcjHPx45HR9lHmr05aRnOPqpFxL2xx3HEz8xJfdeNAFYw8zIAxQC6AtBzOonk2fwKe22EsIIQhdF80GZhFADWS/IrV4VU+uGTnNUZh31fJmUf5nDgOJAh27A9y8Te7HBanI5adsQNkcXdcQCqEVELAKORUOd8XWWnMp0Xxy52jS33i92If9UWOaM3J12SKJC1oLh66x4p7cSJmwT2z5LcenFTXBhe8uD64dO8N8pw7Kp7ZRpC6pzMvAvAuQCeY+bzARymtluZzT++mI/XJY3SALVRHVNiCo0DgPXbs99ga8QpkqB18/2ieOWwZchLZ8MEVcjwE9GxAH6LRP1dAJBbjj4LefjzedLaemty1YgcsRJFCGG2IHOGmKkUKhQ5+2jGqkiixFjCz5QNrkmRX+pWAHcD+ICZ5xBRGwDZV105g/k0A4sxy0CWj78qUBWE2LyoV0OdqOAtb8/AqU+qX8OQ4aaZtWqrhJ6oxTNzl5nHAhhrer0EwB9VdirX2JiBxZhloEf8KVQnuWUCdgJ0MvllnfriOCt+3YUGNcNHomU6joafiD6BSwIpM5+ppEcaqazeGl+Inx7xp4haaTIOjm3TMO4uhGav4kpymYLbiP+xyHqhUUacA0094k+RCyP+Pm0z3/Bb1wme/2133DAiJT02es5a9Dq0QdTdihy3OP6xyQeAKQDWWt6r8sxZnfm+ukxm6+5y5QJu2UIO2H3HxKVMwvo71LWsSxQVqK3C9UmGyFp4HiURnQFgBoAvjdfdiOhjge+1IqJviWguEc0holuM9xsQ0VdEtND4Wz/sQagi7kzIqoCK9YtixRenCuav3RZ3F5TRuVki6zgbqqTNtgzmrDeCnqVqR/vjMiQjX+QKGopE1a0tAMDMMwCIFCPdD+A2Zu4CoDeAG4moC4C7AIxh5vYAxhivNVWMmwa0A6DGzz/uzgFpr6PQIwrLpCXR51OcdkQzqe3ZjVa7tqyLh885XOp+VHLms+kak41rF6e9Vq0M+u40e3mTqBEx/OXMbPV5eE5cmXkNM/9kPN8OYB6AFgDOAvCGsdkbAM4W7260eGmPZOuCXZdm4XVhkrw8bont+6ce3hSAGj9/kzrV0l5v3xsu6WbRevXRInHQrVU9qe39sLhyctVbv++NQjvNkQAsEIjakR3L375J5fKmsrLmMxmRX2wOEV0CIJ+I2huFWX7wsxMiKgVwFIDJAJqYdIDWAmji8J1riWgqEU3dsCGe6ZHXSbZpZ3YuXsocoTklsiVHUlFE9qzZEi5DePe++CM5VLhJrjquNO31rJXh1qxGTqmcaFizuADtm9QK1W6SZZt2eW4TRT7EcY9UfV1+EcN/MxISDXsBvAVgKxJJXUIQUS0A7wG4lZnTHJ2csKy21pWZX2LmHszco6SkRHR3kTJq9lrvjTKQWsWe6RuhaFSrGA1qFoEI2LBNjWzD5HsGHnw+MwsSZrx4/yf5LoCC/Dx0PyQ16p+0RI0cQrXCRCL/OSGVJ5f/6m34F0os8ZjLuOnxVyOiEmbexcz3MnNPZu4J4GnRxomoEAmjP4KZ3zfeXkdEzYzPmwGwL0uUAXh5cpwqKvnBzT9detdnjp+FIXmhquKvZx6Ggvw8NKxZjA2KktPM7p63bUai2cYnirK37zy108HnMkTUnGhcuzj0oruIG2ezgplR2SOnSXeLuZEJJR7dfqmnAfS1ef84AP/r1TAlHOSvApjHzE+YPvoYwBXG8ysAfCTW1ejxUgqUoZkfh5BZqwY1lLbfumGi/ca1i7E+gkImYRd3M0FGPqxC5zfz19m+b85CVamlU5if56pAKsJ+gTWzT2aqCYc87+iWStq1IxOin9zOhKNNo/SDMPMHAE4QaPs4JCScTySiGcZjCIBHAJxERAsBDDJeZxxRxZ/bSTLLIu6RRUlt+SP+3x1XOaBsycZwUgEiBifTWeIgl1BgKhTfWvIN/4yuzQ8+31NeEVpLXyTJ7f2f1FSxOvWwpkrataNso7dLSzVuzl63s0SkVu8EwNGqDXR4P+fwGm0yc+DKRkOeircwR+PaxfhlrdyImVsGtZfaHuBd/OPTmatxbJuGaFir2HW7OHEymm1KUguv+XlyBxldW9Y9+HzTzn0YPdd+1iHK/DXxRVeV1I7ut3132gocG3OWs5sBX09EvaxvElFPAJmRhZABhA0v+37RRtfPvwpxMYmKZsmWDE66ekpqF2Pjjr2hZGqt02InP3KY0Waex431prem44bhP7luEzfPCdS7zfQCIWti1JVSTWnD1Dh6x57KNbajxs3w3wHgHSIaSkRnGI+/AnjH+KxKIzrIDuPX3F9xwNNnGUUUw+Zdcn2Otasl0uAb1y7G/gMcqv3HR6fLAJiN9H2ndT74fEyIhXaRn3qdorWY04+Uk2S1ZZfzje/d648FABwIcX/fYvMbdrCJgQ/Dj2VipQlV6fIf0aKu90YeOOX2tDXNvOpWVydfLYqbVs8UJDJ2CcCVxoMAHMPMk6PoXJyI+t53h1gLeOabRZ4Vh6Io4OFmNMJQUjsReRPGz7/IcuMza6mc0CEV5qtaslckxjwIF/RopaRdM8mbZZgR/2WvTqn0nmzDL8r9H81W0u5NJ7YL3cbL4+0TGgd2TqUrDepim7oUKa6+emZez8wPMPN5xuMvzJyx4ZcysY74WzWobrtdGB+2SNxymBuLKBsVJVk1rmMkcYWI7JnsUjrSPIpasiH4zCjOmr7Vi9QXs0u6UCaHiOO3Fhfp1LQ2mtatVmk7u5lBWKz5AWsV/V592zcK3cZSB/fq2UelFsIzQWLELY7/E8O9U2leQkRtiOhvRPQ7td2LD+vJdXY3++SUhYrT/etVV1MUYtbQkw8+VxVrX2IshqqSZzYvVi5YF9zw3/7uzzK6EwirVowsmpmM8rRlCRfK8EnLpLX/5a32gX1PjVkobR9J/vfCbmmvT+4iJwKnd5t0QbbqEvJb/jN1he37BSZZC1UzFj+4jfh/j0Qc/3wi+pGIPieib4hoCYAXAUxj5mGR9DIGtu5Od3/8aVAH2+1kR60AwNg7+h983qmpmul00g8PqDPMJSFlG/xoIa0QmD1lEo1rF+PIlnXRqr78nIoP/tAH35nOoVOMUMWrj28jfV9WXvu+TEm7o0w3miY2M40gDOqc7nIJGj0ngjmgansGLO46hnMy81oAdwK409DaaQZgN4AFzJxdV1kArO7QPIdQOBXx/uakmyiE4FSN+GsWF6BmUX7gG8uPZeKKltkUi3/PkE649oS2B183qFkkNamnRb3qKC5IjV4bGTOv0kZqE/cAYGCnxtLaOrlLE/z55MSAq6NpADRuwQb06xBexiXKICevyLGoEUrlY+YyZp7IzDNywegDAHsLkAIAdpfLX3w1j8ZnrgyX0Wnlhv5tK703VTCaIgiN61QLnJ281Edi1l9O7xJoH9aZHQA0r1sN397eP1B7IlhdCke3ro+mdaphn6Syf/VqpLsHky6xKKqAydQDurJPKTo1rawk++qEpVLaj3KRNS+P0nST4ib7KlpEhFMmpJXCfLV38qe/WSStrRM6lOB/TNotSYLKBYjMdkpqFQce8YvYqZcv7wHA3+zAjF2OwTFtGuLQRjUDtWeHV47BV3PXYe22Pehw3xdS9metIpVcTBw2oUxK+270lFi2sLfkGr7Wc0TmbyzCq1f0jHR/bmjD78A8j4pJvzG0PWRMOc3Iyq7csbeyH7FvO/uohaBuEpHkspI6wQ2/CKs2JyagXwRUSrU78icu6BqiR5W5eeR0qe35JXkuWCNzZNLc8Lt/94u83E6rezVsLdynvpa/8OyHwgyqHCfUEyKqTkQdVXdGNrNWbg1eN9fDFh5jnISi2bGiLP77EADAD3edCAD47TGHBGrnNZvp8GHN06fNAzqGu2mJGLQwI34RCkIKj822MYayF/msUgRRr0aoriMLyFtwdePo1uGqtMpOVBShSZ1ilD1yGgAgP4P8/Mpq7mYCD302F1e/PjXQd+0uznuHpDJFZxjukbcmBw+RW+1S6ScpTTBicjDJ4TKbhCOrQYsiuqBxnWJs37tfSuyynZ5KUUjDf/mw9MSki3vJT6haq6gmgSiq3ZEA0L5xKqciGT4qG3PeRhDmrFZb99jO9Xlu95TqZwbZfaU1d2OnpHYxtgXUcNm1r7JRvLR3awBAtcK8g6qHYbJe3ZKTkmndQcM57UY3Vi9SMsFKJu9cd2za62Qs/+oAOizWBfYxt/WrtE1S7OpCSRmwXZrbp+2rCquNgnYmoxxWQdOJi3ulZqYy8wXMmBfFy0IqsoadPdhhdZm+e/2xuOPklKMkkyJ7lNXczQSa1gk+/Rw+yX2knfQ/ygwjbG0ScirIz0ONovzA2YR2p5h1/aD7IXJP/iV/H1LJD5vMTP1DAJEz6+JunWqVNU6S71WTVWTEYUW5fg01iXRAeL0erxqx5tBO0aAFv7Sol8pstxs0ycA8uQs7ej+3u31C5qDOjdE5YE1qq4ZQ9cL8tHWKKFxuokRSczcu8vNJWXx3gWSJWwCoViAvfd9ucGF9r39HeTHXgH2uQ82iRKqIKi2dZELkGxPljDKd/iflFQeU6SY9fM4Rob5/1rPfC2+r6hjKTddZ2HUXJ3qUpgYVtaqFKx/qNPouKsjDPkl1LEojjhryg/Kau3Gyb/8BabHRVlo3lP+jtpNUtDpB5RO7rkX+oYkCV48V1Se/bI15p+pkU5dtxrH/GCNlHy3rp+s+mdcp7PIKvNhoSsB7/Hz3iKSVm/273DZZEvzq1ag88zKHxe60iSiTQSNTPYR6IRUunc6aovw87JN0c1Rd2zoMroafiPIBfGauucvM9zFzvKtVgiTTx+0iN8IiQ9DJyiW9gkXw2GE3oKlkcCKYeppdAH4Rkd9V6Te95vj0payNO8JHhbx2ZU+c2Ck9ccj8O1iNrF8GdnafxbmtKzlhldywk5loWCs1qNgW4Obll49mqCnBWFyQH3iwmBX+bwMvdc4KAAeIKLxQdYx8KqmQtXmxUYWUgkydbhFzWCzRteREmBH5D4u9s0Blj/jNXBwwlNaN/jYhtOZjCDva9DotZRQ7sXPp1ShKjW6jiBYbu8CfSLC1SLs1uzlJwtWjXgp9oyKZFFFEhnw7AMwioleJ6OnkQ3XHZLJecjgdgdLS32XJrB4uoRBEEruRcBxBBWHsskhSlsrY6GYKYtO9cgQ+kzRIsZIsMhJkhmSdePVwiIi5oEcidDEKw+83f8a6xtS9tb18QhjDv/JX8ZuqSpkUEUQM//sA7gcwDsA00yPjSSpqXhJy5GYnc1Bo8suWhyltpIgZNjIMcVTeU6l4CKSPPrdKLihjHsVGxUvj7At5iOI0A7quX0KZ87Qj/EcQWcNqj3LQnNllDIBk5C1cdVxp6DbMjJ6TnmXuVGipqCC4j//0Z8RrXMddZtLzzGbmN6LoiAo6NUvEXoctdpEcyZhRsbgrk7iThuJg7pptsRexDktpyPPKyV14iLFobbcw60UlpVqHm3mYEMv9FmNrN9gKw7Dv0zPZazjYhKL8PJRXMA4cYEdFXif8eH/DhJrLQCRzdykRLbE+ouhcWJIREtbyfX5paEQTFBfko2HNIjx49uEAUienqpH07vIKqeGoUagzqsKpHoKZKMpUqkZV2OseQ0XWbiboxRJLspSTOawTIsTSnPR1drfmqCahKIobNR0ibpIL7bIie6w8eNZhAIDDHBIFo0LE1dMDQE/j0RfA0wCGq+yULL6dn1gA+ve3chQu8/MI0+4/6aBAW4Xh4lFV/JlZbmGLTMoc9EttAaPiJoERlJG/7y29TTvO7pYozde1pRqDkBRoeyaA2utL4xanvXbS5akTIjhhrak8p9t5ai2eIptiSYb/xgGV5c+B1P8obvewp+Fn5k2mxypmfhLAaRH0LTTJ8yfsiN+Jx0YvACCnCtcfbHTyZeM2igp785p+/0mhvu/F8QLhs3e9P0v6fqNyHV3XL/H795OcVJekKIRej9VGHeEQhOAUKSPCC2NTNxe3daHj26n9PQ6O+ENG9tzQ375we7IE4/6KeGffIq6e7qZHDyK6HgJrA5nA57MSUSGqizNVBDCaVs0U1QbmrsHuPtOwMer1a3pf9GEupg5N5GvlqFQN9UsyseppBTVrAWBFgMStJKIutKBKslbcXOtOCXaySCbThTX8TslbBcYNWHaBJb+IuHoeNz3+AaA7gAtUdipbSEozT1i40fd37/zvzLTXqiVbr+/nPqOIosTjXkmp8LIo26RGt8YJNx940CplovipZmZloeCMuaHAzV8Et9yMpOCcqqxYWSN+J5K1oe+wXP9RI2L4r2bmAcbjJGa+FoDn8JCIhhHReiKabXpvKBGtIqIZxmNImM57capRZFoVycXjxRv8u5ImWkrUqQ579GJ/BD7HuI/RStS9KXZxtXVtqbYsn8zFUiedKlmZ4Gd0be74WeuGNVGvRiHO6BpO2M4J1Yu7O/dmxuBH5Jf6r+B7Vl4HcKrN+//LzN2Mx+cC7QSmg1Ttm8okp21B/HVWOee4beKUAKn8mYhdKUUngvzPl/mcJSRHeIC7C8MseaCCahLlOZxu4EGDB8YtSK/a5aWEWqMwH+UhfORuSq6yXD2O7Zt+h1FzglWNk4Hjf4CIOhHReQDqEtG5pseVADyDUJl5HIBYrUldhVK6QGqhRsbowFodSyYixUpmrlRXli+JqnDSZJQV4E8vRVQQ7faTU6Gk69Mq+v4AACAASURBVLb5Wxcwp+Y7JQ0B6pPF+iheFAWCD14WWEJYvWQ4CgvyQoXuurlVk4Z5bwSyDQsVhe6K4GYROgI4HUA9AGeYHt0B/D7EPm8iopmGK8hREJ6IriWiqUQ0dcOGYHU8ZS02ObHSqPc6PoCP30ptG615WSSzNt14/YcyZftP4iel3Q8PnnX4wefz1ognEYn+btf0Tf3/dgsUmDeTjJ8H3A1jsWLBvENMi6Iqwl6B4K4868jXK8msII9CRcW4yUar9vG3MAklJqMC48DxP8DMHzHzVQBOZ+arTI8/MnNQPf7nAbQF0A3AGiQWjJ32/xIz92DmHiUlwWrDqk4CkaHWmEus2Fy5HKQMzJnZExaJ34TNbhg3zOfRPT5DRm9/9+eDz91GsqoNv9kNo8rwrzG1u2i9+Gj2R4tuTROPrNbCkNLJpxzmnAsgK47fiUKFooJ+EDnbphPRjUT0nDFKH0ZEw4LsjJnXMXMFMx8A8DISJR0zDlUJWXHh5mKRVrlKgCj+r8/6SFD6el66wqO5prITXtWu3La/zCjdaYfqhW/zTeetKcHqOHthzjLfoXARszA/r5LEgx/cxBCL8hM3eVUjfpVqsn4QuerfBNAUwCkAxgJoCSCQc4qIzEvx5wCY7bRtnHwiqJDYq7SB90YZgNt6Z0cF8fFORHE/DbOOcLVFf182YTWjwmAe8b//0yol+zAvXKqaVQCJ4vFhFnfbN3Y+54O6ekSDIwoiKHwvgojhb8fM9wPYaQi2nQbgGK8vEdFIABMBdCSilUR0NYBHiWgWEc0EMADAn0L0XRi/xUBWCSa73CAp29ZJ7RCoXDwlCG62cF8EGYS/75swqOai36oIY/j9inL5Jc7RXhT7NtdEVhmxUpAfbnG3Z6lzrelUOKe/GYtoMuCxbeQXcAqCULF14+8WIjocQF0AnnnlzHwxMzdj5kJmbsnMrzLzZcx8BDMfycxnMrMa8XETLetXP5hoJYpoxa5m9eQo7DmFINauVoCTuoTXJnEzhi0CHoOofxxITa1VG1YgHulpUY5q5V7c/khFOj1ANDpNbUtSyqK926iLIioKafhFFnf3lvtrX/TUrl6Uj+tO8A62UI2I4X/JiL65D8DHAOYC+KfSXkkkSA3NlaZFyNKGzinisgquH6HwggfcY7iPDJg4tMFHBaFlmxL/z3VbxbNT568NJvGbyYa/i0fIbvvGtX3PTn/dKRZgYD1VVfiwzesUM5arkyQoyCepqrVmDsbx+7QZfgY1QwLURJCNV83dPADbmHkzM49j5jbM3JiZX4yof6EpyA8Q+mU6gbu2cjaMskZRQY2vKH8YYC8YBQQr7g2kSzz06+AedZWU3P1slvgE79v5wUJ4s1l6OghvCIbhWgXUdu1TWyXrP1NXKGu7MF9decSgPn4/tkCG+zYsXjV3DwC4M6K+KKEwP8+3HMFikzaJ2w/aVFJpvjOOdE5Rl4FbWKs1g1iU+SZF0ucv7e66bdIUj5gsHk1irfokShDBvEyByH/k01OCom4ltYvTXu+UVC40DhKLu2oMf3HABC6BHMmDFCoO3RVBpAdfE9HtRNSKiBokH8p7JomC/DzfC5jmk8ptQVJWtqXKaI+kzrsTRQXBZi3bTDMFr/9DEBVMs/271odPNKjdP7RR/BXVCP4yj/0y5rZ+B5/vVjziD8rx7bwXPz+ftRaLN+xUoq6adPWMX+hvxjn2F/HtaysSmPODiOG/EMCNSK+5O1Vlp2RSmEe+Y37NYWmiIX4fTF/pax9R4VWCsU2jYJE2ftxn5ptnkFj+UxSI7VkX1D+9+XjX7WcOPVl6H6zs2leBNVv3YHfA0biXKGHbktTvEEVB9CC8ckUP4W1VhIwmffWTlvyKL2eLRyZ9NTdV07d/R3fXZyaIFYoUYjnU5hH/srQgiWQPf8bGfFGIZlT+6T8/e29kcMQDo3z1Z8WvuwL74r0u8At6tArUboUP91lPU76DqOTBzr2pfjsV/gjDdEsJQq/fuY5CSY0kyTWQ/5tYFuj7IrLXSWloq+snU/CTba86Jv764dOEt11tCly4qKdaqRgZiBRiqUFE9xHRS8br9kR0uvquyaEgn0KlX6u4O2/f62+01ffRb3HGMxMC7cttcRoA6gYovg0A5T6iKsyus7WCkT3PfZeqyCQi9ytSmtGMdfFOhYvlk59XB/retGWbvTeywS0jNcl9p3cBkNnRT9lOc0lh3ioRGc6+hoT+fh/j9SoADynrkWSYgU07M6fSUlCW+4ibN3OsonhqP+4zs5Ed8vR4Fd3x7bqwRt+pKIRjVZ0UpVHA0biIpn/yOKtC9NPExZu8N4oBp0LumYSI4W/LzI/CSORi5l2IvoZFYCYs2ogVilQhswFViTt+UubNI/49PhNjVGGdyalILjPr9JzgEfJq5i0f0U9mRH6RJRsTEWsvjF0SaB8f3Xic6+d1QxRc94toLL9oQqYsssE4ihj+fURUHcZ5RURtAWTdEFpV+Fem48ee+SlK78d95nd06UfZMQ5EF6jNbq0zXapKRcm/v0240EYKCrXttLglvVyHk+8ZGKxjARBNoFSlTeREFFnSYREx/A8A+BJAKyIaAWAMsjC2/5e1mW1MnDDbGFGtebNglJ81io9niF0g+ysO+BqV+vUqDHpinL8vALhGscCamSBFOqJI2nHTfArKQ5/NPfj88BbexYL8SqGHUdl86LN5QjfhYd8vPfi8uaTcGzeywO4LRfV8BeBcAFcCGAmgBzN/p7Zb8hGZFu7cux/3fOBPbz0sNTxi+M2RGoOfEvOPv/5D6kRv11g8Pv1pQUlj0e2SWBcdg7oy3EjmQohedCtD1AYIMqJrFoHBUTHS/GpuSrr6uLb+BMYqBK65Z3yeS1ZeGb/UeyMTUYRSulVasxKXBLxoClk/AAORUNTsq6476piy1Hsh6NUJS0MZpSA/4vBr3IVOgxRnNod+ehW1CMLTgtmiSS7smR4yquLmmnQnif4Ef35HPPw2KGYbo7q0IqDGt5zmTfG5AxGp4jKfNYwB4PJjU3UN/AY9qFiDsEaItfAxu4trjV0knPM5ANcDmIWEfv51RPRv1R2Tzbw13q6esMJPQTTCux/irti4K2RqfSb4Gwv95LMHREWsv0yiiJtX8VOb2/QzkgXERvxBDJ+5F37lpgdJULu1sm1Peo6Nnz7FFV0lckWeCOAUZn6NmV8DMMR4L6sQGlmE/BGmLpNfWz6ImNbG7SnFxkww/HYXQtBQRyd6Haq+mHgS0YVtv4YyCvyuA5jPn7Vb/UXH7djrnXQ4fYX/nAVz5TS/57db2cUkgzp7qs5LY3NAraywiBj+RQDMqWitjPeyim0Cma9hlV6LC+Rr7gQR0zILlWWA3bedXt844iep+4iyxoloElrUiNxo/P6b1piO9axuLXx9NynH7UaQUGuzJLjfyeRhzb1nho+cd6SvNsOMF2esUCdf7YbIv602gHlE9B0RfYeEHn8dIvqYiD5W2juJbN3tPXIOIvo0wuSjV+E/DCKmZV5rKBK4MsJEOtx5asdA31voI3RUhChH15lwM7VFoF+XHONc99cLvxIJqhKZzPIaSzcGX6R3IlN1jGQiYvj/AmAwEmGdDyDh6vkLgMeNR1awUaBwSBAN8eNMaoIiWilzVvtLJgmyuLt4Q8qtJZKYdNVxwUMheyquOzz0jC5C21GIZQS/ETeZmvQqcpM/r7u/UbuZBjWLvDcyocrNeH2/VMnTr+etc9kyGNUKU//HTQJ2I6iEOADc92G0UYRJPG/JzDwWAIiojnl7Zpbv0M5yROK7//erBb7aFBU1C0NxYXCr6VQ2UhZXCt6UwhiZoWce5mv7Wau2CNUPjuK3MyMi7x0mnNFvzWQ/Qn5+qONTl8kvTU2RcGu27kHDWu4L834kma2s2xZPLqxIVM+1RLQWwEwk5JizSpY5SvYIXOjmhSkRrJmTKhBVILXDrX6pmb+d5c+4+iWMj7+NTy1+USXWoGJrQbigR0vhbft1KEHXAOU+/eoZ+Y1yE11UbeRhiMNivjmKRPr5KUOaKYhctXcAOJyZS43Si1klyxwlu3y6ZUbdeoLnNkGyRP1iXpvwm4tQXTBT87chfMsiiLg5nGjfpLbEngRnYKeU4ROZSe0wDQq8woLNEAVTIxW9ySdRNePp3Mw7g1gWIrMW8yXz93OOUNgbeYj8kosByF9BqYJ8KCh5kKRjU3eDYze+WhFQpdON49unBMQWrPO36NpWMDNY9Xqo2SiNmLxM8d7U8OJlRx98LlJC0ux/7nWo+FrLd79swMyVW7HOo0iPaCF3J/zqY3URNOhRhiiL1PL416hfDj4f0ElcjC9ORAz/3QB+IKIXiejp5EN1x7KRzbvCXShOmF0xL48PpqroRi1T9MVrJl0TEURDWKOMhLn3g9nR7Uwi5puXSPLT/JD6U09+7b7eJNIH9/b9ZXhXEyxBmhdhyVq/tTxqFGa+JDMgZvhfBPANgElIlV4UL00TM2EWggYf7q/k38yVauRfzWFx/zfReTQ7d/U2nPVssIItSWQVkLcisqj4ztQV+NN/ZijZf7bhN+s1iIne5hG2GLWOjKi7zm+2bhhGTPIn4RK0sFHUiFjFQmb+s/KeKOK6fm3TpmJ+8DtKVRX/W6MoH78KJB7//fN5+DnkzUdFEpood/53Zmz7VsFhzetgzmoxRVUrIouKYQ3zXg8ffNQRJ5f2FlsHUlE0xwmvkqdhZ0VxIXKL/cKI7GlGRA2SD+U9k4RosXQ7MiVe20vBU5OZBDX6gNjirnmTxgG0gLwizFZtiXZpT1TS2U9I6taQkgheay0iuTuZiIjhvxiGnx8pN49nOCcRDSOi9UQ02/ReAyL6iogWGn/FQxECEmZaeNqRzST2JDg1igrSQtguemmi7XYyBkL//HJ++EayiDDRQCoRWdw1C3zVVlAM/vrhcmU1ZFHNR97J3opwhtnrBhxEmDETENHjP9TmIRLO+TqAUy3v3QVgDDO3R6Kgy12+e+yTMLbwtCMyxfDno1WDlNTrpCU6d04WV/eNroCLH8RG/Kltgrp9dodUf7VDtQ/e100upF02l8+0I1sr+4kkcNUgovuI6CXjdXsiOt3re8w8DoDVQp0F4A3j+RsAzvbZX98EPQmv7FMaSdEGEWoUFWD6cm8xpzD97eQRWiqD+Q9axwHxU9/HYlyU9WTLfS7uBv3tBz0xNtD33Aji9z7GRziqFTfjG3Y8vmbrHrw6wTnSTYbhV539bofInOk1APsA9DFerwLwUMD9NWHmNcbztQAcNVKNdYWpRDR1w4bgKdFBL4h7hnQOvE/Z1CzOTwu5dCLMbcpv2FoQ/Jbl88tzv+0utJ15dNyvg7gErx8XQ1hE6h9/NTeYTo1ZmtlrRJtk+NXuBYPCIqKa6YRbzYpRc9YGbjfJg5/OdfzsxYBF682DiHJF0hZuiJzJbZn5UQDlAMDMuyAhH4cTV5/jrY6ZX2LmHszco6Qk+qSIohAyBrKpUZSPNiX+ZAX84ifzM1Pp0VrsGGatSkU++fmdj28X3Xl4xbApntuYj8MP1kgVEamRPm2D1TuwFilxomZx8EGB201y5BT/wot+CFqD47p+8YofiJz1+4ioOgwjTURtAQSN81pHRM2MdpoB8CdcE5K46luGpUZRAS7o0cp7wxCY5QKylTqCrphJS1JlOEsb1hBuv3Ed9VW0/BA0lLB1g/RjnlrmrSkkovJqx4SFGx0/MxcZunVQh0DtA+4KmuZ9qFh62FsebLR+3QkphdE4zJKj4Sei0cbToQC+BNCKiEYgsSh7Z8D9fQzgCuP5FQA+CthOIB4bHSyeP25qFOULxTiHWZKwFkTPBK7sU+pre9HjN+c6ZMo6ThCCDmTutrgxL311sozuHOTwFinphRfHLnbcbpupRkaYBeGNO5wz5s03R79aQyIErSthPt44yi+6/SdKAICZRwM4F8CVAEYC6MHM33k1TEQjAUwE0JGIVhLR1QAeAXASES0EMMh4HRnDfWbhBSFIqUQvRAt1Wy8dP5nHqjJ2w+BXLlk0NPPrgL5xK5t37kPpXZ9h5BT155UdIiGfdhzSQHyWEwSzlo5bQuEbE8uk7G+fSyz9ys2pNYxCxdFGQdeAVguus8jErad1iehcIjoXQD8AxQCKAJxgvOcKM1/MzM2YuZCZWzLzq8y8iZkHMnN7Zh4UtaZ/FJneXf4ySnqbov7P5Oi10KiU9MQF3YT3kQm1ecNiHr1PLXM+tWQpnq7YnEhwemtyTIY/oKunQLUBFMz+fv4759mAH0Rj6VWHmV7ZJ1ho8KAnxknuiTeuhh/A6QDOsHl4hnNmInEVNg6LqPRx8rQur2Ac0aKuUGEO63f9cHEvtesOYfjNC/ZJbmGwup6SdleVPfHKOnVzcbihwuVhxhwhVltx0RTr/ty4wqfr0C8iVf4yBbczYBkz/46Zr7J5/C6yHmqEa5eaB+2FPuujmnHzHZsjQprU8e8eGtQ5FcH7yvgl2J9FCTDW4036ZlWtE7w7TW1EiirMhW2iqF+7T3AGp9qdGUc8flDcDH/2z/0jIkzyiQiiI3ez9kqhz1Gd2Xa5+WW3m8Lz+rZv5LidEzf0T0UzPPTZPLz+Q5nvNjKBnXv3H7xB2tn9l8eFl892q72gYi1JFred0lH5Pm4+sd3B56IjflJs0vwWo48TN+twWWS9yHKWbhSQzgxBTcHFXTNh8hDcyj2+O3XlweeNa/sfQR1tibWfsMg+3O/EDA8vnbps88EwPOv6yKotu/Hw5/Ok7MOJuSEE4ACg7JHTQn3fjRb1qntvFJLbTk7dXMoFR/x+Jmb9OvjP2bimb/YUJnS0DsycndUsYiDIBM9P+n8QdU6/4mNmd4VbJuQ381Ozioa1inz3y4rTAuXj53cN3bZIYlJQfli08aCP32pPZFVJC6PuKZP+HTO7qpToYr2ftZggETp+azfHSeakp8bM4g2peNzuppR2EUQXX834ycQNYvj9unrMzFzprAtk1iYRDTN1w8nw15NQ0EKl4Xxx3JKDPn7riD+T3TBBOKpVZmd1vzlJrNSmH1ePSMlFK9kUGSdkHYioOhGpd9wp4pObjvfc5p2pqYW0hrX8ZWhaDfOi9d4l8c7u1kK4fdHFXTNhXD1uxSdECoT4wak9GQumQTTq/XBwDdzS1c07szN6zIw5ouiPA9u5bJk9THNxnVmxFor/cvZafDl7Tdp75iCI+07rHDi7OQ5E1DnPADADiexdEFE3IvpYdcdkIjJ63GQKjbvQpzyC1TBvEQgbvfxYsWpDgPjirpkwI3638o6yZWhFxOeC8u60ld4bhWC/Ia5lvd5lDvyc1ltUDy4/mJ763wW5CZcovukGwU8hdOvo/frh0yrVJ9iwPRW+ecyhwbSM4kLEOgwF0AvAFgBg5hkAMlPE3AG30e+bE8uwfvuetLTpZvX8LVpaR/wiqol+LqYaNq4kr8o/RQVqLMOyTXKrMv26U02BegB4eoy/Yt9+Sa6FWI2EzCn/mq17pLXlh7A3zaFn+Mu6joIBPgIGRAY45tkqhxaAjhYRw1/OzNb4vqw6SqfRb9nGnbj/ozn4w/Cf0nRD2jSq5at9qziYzKLrxQV5tgk3XlmbYUb8xREqk85YkVhPMF9oUcofh4EdfPwycXIbumnEi/LlrX0dPwu7PpKJoY1+fieR6+eL2SnJ5yAu0NoKZ7teiFxhc4joEgD5RhGWZ5Aow5g1OJ2EyfjfLbvL0wyPX9eKVbJ2okn90UyQaA/zbMJPFmSYkoKyJA38YC40M//BwZHvPwgHo3o87MlZ3ZoH3sfYBfa1KD6fFV5n3hxZZk6CkpGIlHlm31+R9j8ObO+5zdqtKY2dIPIZnZvXcf187/4KfDB9pRJVYRHrcDOAw5CQYn4LwFYAt0rviUKctEnYFI7ndIGJUCyoTdL30W99t22OnPGjsVKYQfUERFi2SU4uxINnHy6lHRGSrh4vt93Vxwf3jAaJLhHFbAiThouZcet/ZoRuOxNVT/0svvYSSMp8eXxq1hXE8Juv5x2WtZwde/ej431f4k//+RnjXaStgyJSc3cXM9/LzD2Nx33MHI/jMSBO4kxJv1zYc1RkcO0WKeOGecTfqan7CMFMEFfPv35zJADghUvFKlnJYvWW3YEFx6xc1rs1PvhDH+8NJfCtkdPgZU/CiIOpXKA2L8Au2ZgKZ/7459Wh2/ZT0rIqECSGv2X9VKLb0Q9+lfbZ5zNTEUQqZuAiUT1fEVE90+v6RCRfglIhBXn2h5kc8S9YF0xTO4mIkb3j3Z8DtV3D5Af0I1McxE9fapy8QcJHw9DnkW+kln6MQhgMSA0crL7jBevS/fJ+DH+UyVLmUfksietSANAhghrOmUStAOfchT0POfjcatxHmwJEVKx5ibTYiJkPOmCZeTOAzM6nt+B03clynTUVECsLGp1R0zTi32+qzemVjBJEpC35jbUxRJL8KFAFSpSo3AzJn8O6t+cscsOiMsUA0D+AVIAMZIfpZmmxO0+cRAWDJDN2buZ8czTPwH5Z650X5BcRw3+AiA7emoioNbIsqsfJEMgKwRK52wetmWt29aza7FywYbMlLDKIqyf5b7rjvzM9t5UtyfzdfHlVOEXNfhjfO5AqhOJ1n2nuQ7tGtWa8E/skryVErVTpVTT+quNKpexHptvFujZoXmdcsiG15vXl7PAL+VZErMO9ACYQ0ZtENBzAOAB3S+9JRJjD42SNStqWVA7/tI4MrHVORTGPJNx8/LdZXElBMnd37BXXtmnkM7vZiwYSdH+SiLruwoqJbTPWbbxmGH5+i1MOaxpIoiMsyaieoAXc3Yii1vWvHrUJgoj+HdtGbVKW9R5/xbApttupqNAlsrj7JYDuAP4D4G0ARzNzVvn4zbzpkpUalML8PLx9be+09+5+f1ba66ASCmYj4FaJyxoqGmTEP2e1+0VvlmQefHgz3+27ITcxLGVozH220q6xv3wNADjMFIK3wSi8IXOQ3rhONcz926lp763cLDdpzo5/fjkfAPDY6AVS2qtfM3Ujt1u/kX0zsAue+HlFKkQ4SCax3dKgzF7bDRjsEjNXK3C9ilqHCgDrAWwD0IWITpDek4gwK0/KPPe6WGJyP5oRPjKie+v66NoqJRhnXkS0njPWos9B4vjNF8r8tZUTeLaZimqoSrKqI2FhNt90xbopjVrzL0Q4wlSUPjkdV63z7nUuBTkOJ6yKpjJcJP9j4zqU7Qnasbey4TdfE35i+JPY2YcZy50FDGXQ/1/fKW0/iUhUzzVIuHdGAfir8Xeo2m6pw6z/bufj//yPztmMfjgg4a7y+lW9cHGv1Mq/nzjkIDMMs5H8z4+Vqz+Z/bbFARRJRZDh4zbLRdtd7+0a18KQI5oGKkF4ik0Be4egMWn8a9Qvrp9baxyEwWog7z+tS+g2P5yxulLNClnhu0msomoAMHpOyjceZHbXrG5lV+BOxcqrdkEgPz9wsvT9iJyytwDoiUQpxgEAjoKh25ONmEc0drZ5085gdTOtIyUVXk0/RjGIq8fcvt0o1vz/CpMZ7MYRLf1JYtvRsUkqWqLXw2Mqrbds212OOtWCxZl3sglTjDtZSeb+rVnnshQnJ1nalTEwMs90rMKIFQc4LSQyyP9oYOf4gxc/uel4X7U7RBG5evckE7aIqJiZ5wPIWolmc8F1u1Ov+yHBRk8Na6b7EGWPaIB0X3LSJ+tEkHBOtwXbYROW4oR/pTKPw8g+u/EHU2nGoFhvkDeMSFdV3LanvJK+kih2o0A3DZggI03AX9SUDPeYap78OrF28MRXC1B612dSro9tpvUb6wKo9UYThHoKDK5fzEleMhG5elcaCVwfAviKiD4CIH+FNAbsFpiCJi9FEYZnNjCvfV/mum0Qw3xm15SmjHWR6Ymv0hf9ghSfESGMuFwS629hVkvdu78Ce8oPSDWWbr/8F7cEcx3+49wjhbftrTj6RAbrtiVm0knF1AoJI/5lG1OL3mb5BEBO9bU+7Rrh7BA6S0G4oEfLtNcyChLZ4XiVEdGhAMDM5zDzFmYeCuB+AK8COFtJbyIiOfXPqmQE+FMXDOKKMd8sRkxeLrxtpuGmafSvLxP+8qAjfjvc7vkybmRRcf4L6dqLlxxziMOW4flSgsjclS4Lz1PKfg3dPgA8edFRtu+bB0alDYOFagPAuDsGpL1+Z+pK3Dxy+sHXqtyIbmflf40dj0m+wcxjmfljZlYnoh4B7e79AoD87ELVAmG+RNoCGBy371hFpFTht+ylHW4XyyuGnHFQH78dcZfca1BTTg6ENXv61kHeCpVBufM97yRBL67sU+r42fY98s7XQ006PMlfeoJJOO29G4JrQ9nlr3wiQSvJC7f5bh4R3QOgAxH92fohMz+hrlvqSbh50i3/NSEzOa1hdcws9Y6tOqonSh1+J6JaKK1TXaJfPGYhSj+ZwX6QvajoJEt+z5BOgdpzi8raJzHDtqR2caWoJHO9AictMBHiOnXcenwREvH7BQBq2zwCQ0RlRDSLiGYQ0dQwbQXl3g9nVxrxy7Y5b/xQJrdBHwQZ8Wey+0Y2Mkf8quP440JUblwUJ1nyM7oG86O73ZjGLwwus27lvO6p+thJk/HCWJMeU4ifP67JouOwh5l/IaJ/AVjOzCMV7HsAM8sXmhbkrcnLcc5R6QXP69UIN2W2/oZDP5mLK4+Lp0plEB+/VUxs2rJfcXRrb13ybCSMj79RrSJsNEkEZFGN7YzELlIqLMnFZBn0LE1dA+9NW4n9FZyW85KBpQc8cbUOzHwAwG0R9SVyrCP+hiF9pU5uigMHOHINlCCjd2ulsvOen+i7jWwhjHRz+8bpE964ffxhcHLvjb9zgO37mQ4zY+bKLWkLrubs9yC0MWlxjZ67Dje+lR4eHKaEYlyzRRHr8DUR3U5ErYioQfIRcr8MYDQRTSOia+02IKJriWgqEU3dsEHetC2tZR3XcAAAE59JREFUExbDeUJISVynYisvjFuMUXNSIYWPnd811H6S2GXXJgkSx+9n8Tjb+HD6qrTXYVw91oxv1Zm7VrbuClbUxw/1JS0YR817P63Cmc9+jzKT9tNNA9op3WeYdam4xgwip+yFAG5EQrZhmvEI65c/npm7AxgM4EY77R9mfomZezBzj5ISNRrl1vFyWFXE3Q66MJOXpELLxt85AL85uqXtdn5xU6EMUnoxLkngJCpjpq3lBMP81pUnWtH+3+zkCYJiFRdMUivGQuBhuN2m4FGzut71MnINEXXOQ20ebcLslJlXGX/XA/gAQK8w7YlgJyp20UuTLNuEM/xOd2+zznZTiSdhMu3drjhEEB9/kGISsjjtiGaOMdMqCDNK62YJOY36filzlHhUwEz1bOJwk7CebD656fhQ38/YET8RXW73CLpDIqpJRLWTzwGcDGB20PZEGXqGd9nCsOGMIr+hTH9w0lU1fFLlROogUT3VBUfBKrRDorgAnjGyRsNycc/0xCZVPv5+JtfjbofFxEl3D1Sy76pCEJenHzqGLDGZyT7+nqZHXySUOc8Msc8mSBR2+RnAFACfGZr/Sjm3u7d7JWwMucj3ZbpTknIndskqKt02h7cQL/ouillYTRWPfyVHa96K+V+9XGJdgfYmnZ/pK1LJVWZBMpkzyKrIkxeqnUWGDYHOuHDOJMx8s/m1odvzdtAdMvMSAHJWN33gtXA5+k/hSwxE/SMmXT0qIobcOOcoOWsUSd674Vh0a5W9LgfzDd8sZNelWbgbpDlh7/NZa9CnbSMAlYu5ZzP/PO8Ipe2fdqTcgkFVhSC3q50A4glOD4FX1msHCSNOO/eKzAxCK9OXb6kkP6uCnRa5BlnT52qFeahfoxBHt24Q+8KyH6w3eKcb79vX2S+cimIW6Bo+KaWddNNb0+02z0rOP1pe7eaflqdLTgwMUG4xV/Ac8RPRJ0gFwOQD6AzgHZWdylbs4uof+cJdPtkvfds3wnhDJ2Tumm146ms17gsz67enJ8MElRq2MvSMw3BRL3VCYKqw+mWnLdtsu13Y7GDVNV+dsAqHBaVDk1qe9Y9l6f0DwLnPpYvMbdqZ+ZJiQSqDyUBkxP8YgMeNx98BnMDMdyntVZbS2WZqP+z7pTZbBuf1q9IDoMyaIaqw5gu0rBdcjbAqYHetlt71GWZLdrnFlRgmSwP+sz/2xft/CC5gFpYZKzK/XlReHuHnB07GgocGV/rsneuOVbZfER//WAAgooYATgCwB8BKZT3KYsKGg4pgdYm41ZSVwdbd5TZVmeS0nW2y2F6c/swEqe3FZfhljcIL8/Ncw4pVu/eyxXvoFCXX61B1ciluevyfEtHhxvNmSIRc/g7Am0R0q7IeaXxhNcqy6frX0ZUWxuOM988EorLHUWcER83Llx+ttP3bT5FXKDBK5drTj2yGri3V5R4A7iP+Q5k5GV9/FYCvmPlyIwb/ewBPKu2ZJmNoW1LzoB/7L6d3yYqF2OKCPOxVtLAelXR0NmsAJXGrIFVSS20oaqHEO2dJ7WKs3Lzbe0MJPHtJd+X7cPvPmAVBBgL4HACYeTsAdaEqCrHL3tWkc7HNYqvZnfS7kDULzMguhGPGzT867b5BodoWMccySlNWBe2klvWd14OsooCyCau9ZSaKEf/U+wZh7B39le8HcDf8K4joZiI6B0B3AF8CABFVBxB/FeIAzH+w8gKKJp1/nFs5rnqZkZQUNi7dilXsTCZdW9VD2SOnVXq/7JHT0NClqLwIIiNxuxuoX6KaWbVQVMjFizAKqUlu6N/W9v2yR04LnVVrJkyxFVEa1SpG64Y1vTeUgNvRXA3gMABXAriQmZNL5L0BvKa4X5EiK3wtbhrZlHGTQTJOff8BuRM9lSN+lYh4YA5IOLhscKmFQcYNp76NK0nGTdeKzLDTTMDR8DPzema+npnPYubRpve/ZebHouleNBwSoliylUGdmzh+9uwlctLHuznoi6uOKvKKyfZLltp9IVfPMRIiMqwzi/Xb94Ru0w7V6pULHhqMKfeq0RSym33ddKJ8GeYGNbPSyeGIdnpLxu2kO/1IObLDTpmzqhcDB3RUI4+ddQj8m2UluZnp9fCY9NelcsL9XrzsaOluPDNFBXloXLvyzeWAhDu/3azogIyGLVxxbGml94Io4GYK2dvzkBypKFwqigmhk4G/pq+chVenikKDujjPZnIJkWzL9hIkQLx2c3mf1qH3AQANaxXj3O4tvDeUjAzvibk6VhIZbjYr5hDmpIvq5MOy93rIWcN/6uFNlbQbhfvC6by+rLccQzDrr6fYvl8SclG0qtCwVjHuP72L8v14FUNpWie7lTllhMXa/Y9UzHz7tE3JZyQTq+JaFJeBiB5/NSK6kYieI6JhyUcUnVPJH/onXDL9s9B90a6JvRtBdXz5ce0ayW0wW1d3AVztEtZ6UU85wmP1alRerN9hEsyTWWAkWxeS7dyeKi4D8+Ju0mYMPiJ7lT9F4qneBDAfwCkA/gbgtwDmqexUVNiF+4VFRoiaFx0U+I9FqCm5HF/UZv+CHnLlpJ24a3AnZW0Pm5DSfpJprKNIFkteb+MWbMDiDXICBezWJlSPxDs0qa3EdkSJyJXcjpnPJ6KzmPkNInoLwHjVHVPF2Dv6Y0+5uvyztiW1UJSfh3025RBlcXRrdRoeVZEmdYrx93OOwICO0cj02o3UZfGEqZhMkCprTkQ54D+hQ4m05Cq7G5bqmW+FgsXjqBE5c5IZvFsM7Z66ALJW6Lp1w5pSEzvssFsokzX9BxLl3mQpKDpxjcQMXSei8PQsengwJt41EAM7N6lysdgy6dchOy9pOxtvJ48uExURW1EjYvhfIqL6AO4H8DGAuQAeVdqrLMdu4U/myVJUkIe3rkkv8jH86mOktQ8A90WweBkFBfl5WW3wn7qoWyT7qWW4KN20dTIRu9G96hF/VZDS8DT8zPwKM29m5rHM3IaZGzPzC1F0LlupYVO0XOa0HACa1E2PsDm2rdqiHXdIVDrUiHNWt2jCLJMhqipi4KPk+n72Eg6adESiepoQ0atE9IXxugsRXa2+a9kLEVU6AdvaxBuHobgg/eaiehBy4wD52ZCazKHYEDBUPYBQjcqF9eRienbfGhOIDENfBzAKQDLtdAEArcfvQfIEvHFAW3z2x+NxfHvJoZAWopIKlolqX2zUdG5WR7n8gSqqFeZjzG398NRFcmRFqiLJK6wqnLYiUT2NmPkdIrobAJh5PxGpLftURVAd8nVIgxpY/usupfsAgMGKkt2qGl/c0jeyfanQvpE9K61qJMdWKlVlo0JkxL/TKLvIAEBEvQHILS6qCcTNhi7QyN/39tgyHKoqbmX75TP+zvhUXe20bzRqueq4RKRbFBLNqhG5ov+MRDRPWyL6HkAJgN8o7ZVGiPN7tMJ53Vsqj1r500ntlbafrbRqEE3R+cJ8QnlFtt8m1dNH8frE3YM74a5TO2V1lFgSkWLrPxFRPwAdkXBz/cLM5R5f00SEypNwxDXHoEW96q5VlMJQFXylUdCmUS38sm573N3IeGRW3LKDiCKrt6wat2LrPYmoKZDw6wM4GsDDAB4nolCpo0R0KhH9QkSLiOiuMG1p1HFcu0YobaSuIpC2+2K0llgvoirToKa6jOmqhpuz6kUA+wCAiE4A8AiA/0PCv/9S0B0SUT6AfwMYDKALgIuJqGpkC2k0Csjy0PrIcJIT11TGzfDnM/OvxvMLAbzEzO8x8/0AwgR19wKwiJmXMPM+AG8DOCtEexpNbNx5aiKxTaWExlnd5BTwqepUt0mc1NjjdovMJ6ICw80zEMC1gt/zogWAFabXKwHI1RvQZAWqiuFEyfUntMXWXeW4ok+psn1Uhf9TFMhWj63KuP2nRgIYS0QbAeyGochJRO0QQTgnEV0L42ZzyCHyiydr4qcKBEcgL49w95DOSvehfddiVAUNnahwNPzM/DARjQHQDMBoTqVZ5gG4OcQ+VwEwS1W2NN6z7v8lGGsJPXr00F5OH0RRE0ATHbWrZZdwWtQU5BH2H2C0aaQT0ERxtRDMPMnmvQV22/rgRwDtiehQJAz+RQAuCdmmxuDNq3vpDExNTjH5noH4afkW1M0yZdE4iXxoaEg+3ISE/k8+gGHMPCfqflRV+rbPvlKSGn+U6vDONBrWKsZJXbK38HkcxOITYObPAXwex741mmxHeu1jTc6R/aITmixGL8YFQY9uNWHRhl+jyTL6KZYm0FR9tOHXaLKMbKy9oMkstOHXaDSaHEMbfo1Go8kxtOHXaDSaHEMbfk1sNK5dHHcXNJqcRBt+TWxEVcGqKjDujkSZx5cv7xFzTzRVAS3qoomcL27pi4mLN8XdjazikIY1UPbIaXF3Q1NF0IZfEzmdm9VB52Z14u6GRpOzaFePRqPR5Bja8Gs0Gk2OoQ2/RqPR5Bja8Gs0Gk2OoQ2/RqPR5Bja8Gs0Gk2OoQ2/RqPR5Bja8Gs0Gk2OQcwcdx88IaINAJZJbLIRgI0S24sDfQyZgT6GzEAfgz2tmblS5Z6sMPyyIaKpzJzVoif6GDIDfQyZgT4Gf2hXj0aj0eQY2vBrNBpNjpGrhv+luDsgAX0MmYE+hsxAH4MPctLHr9FoNLlMro74NRqNJmfRhl+j0WhyjCph+ImoFRF9S0RziWgOEd1ivN+AiL4iooXG3/rG+0RETxPRIiKaSUTdTW09arQxz9iGMvQYOhHRRCLaS0S3W9o6lYh+MY7vrij6L/MYnNrJpmMwtZdPRNOJ6NNsPAYiqkdE/yWi+cY1cWwWHsOfjDZmE9FIIqqWocfwW8MezSKiH4ioq6ktudc0M2f9A0AzAN2N57UBLADQBcCjAO4y3r8LwD+N50MAfAGAAPQGMNl4vw+A7wHkG4+JAPpn6DE0BtATwMMAbje1kw9gMYA2AIoA/AygS5Ydg2072XQMpvb+DOAtAJ9m8PXgeAwA3gBwjfG8CEC9bDoGAC0ALAVQ3Xj9DoArM/QY+gCobzwfjJRdkn5NV4kRPzOvYeafjOfbAcxD4gc/C4kTF8bfs43nZwH4P04wCUA9ImoGgAFUQ+KfWwygEMC6TDwGZl7PzD8CKLc01QvAImZewsz7ALxttKEcWcfg0k7WHAMAEFFLAKcBeCWCrh9E1jEQUV0AJwB41dhuHzNvyaZjMCgAUJ2ICgDUALBacfdh9MnvMfzAzJuN9ycBaGk8l35NVwnDb4aISgEcBWAygCbMvMb4aC2AJsbzFgBWmL62EkALZp4I4FsAa4zHKGaeF0G30xA8Bidsj01yFz0JeQxO7USKhGN4EsCdAA6o6J8IIY/hUAAbALxmuKteIaKaqvrqRJhjYOZVAB4DsByJa3orM49W1lkHAhzD1Uh4JQAF13SVMvxEVAvAewBuZeZt5s84MWdyjV0lonYAOiNxp20B4EQi6quou059CHUMmYCsY3BrRzUSzqXTAaxn5mnqeumOhN+hAEB3AM8z81EAdiLhmogMCb9DfSRGx4cCaA6gJhFdqqi7Tn3wdQxENAAJw/8/qvpUZQw/ERUi8c8dwczvG2+vM1w4MP6uN95fBaCV6estjffOATCJmXcw8w4k7riRLGYZffRzDE44HVskSDoGp3YiQdIxHAfgTCIqQ2JqfiIRDVfU5UpIOoaVAFYyc3K29V8kbgSRIOkYBgFYyswbmLkcwPtI+NIjwe8xENGRSLgGz2LmTcbb0q/pKmH4iYiQ8EPOY+YnTB99DOAK4/kVAD4yvX85JeiNxPRvDRLTwX5EVGD8YP2Q8Mtl4jE48SOA9kR0KBEVAbjIaEM5so7BpR3lyDoGZr6bmVsycykSv8E3zBzJSFPiMawFsIKIOhpvDQQwV3J3bZF4PSwH0JuIahhtDkSGXtNEdAgSN6bLmHmBaXv513SYleFMeQA4Honp0kwAM4zHEAANAYwBsBDA1wAaGNsTgH8jsVI+C0APTq2ev4jEiTEXwBMZfAxNkRiRbQOwxXhex/hsCBIRBIsB3Jttx+DUTjYdg6XN/og2qkfmudQNwFSjrQ9hRJ1k2TH8FcB8ALMBvAmgOEOP4RUAm03bTjW1JfWa1pINGo1Gk2NUCVePRqPRaMTRhl+j0WhyDG34NRqNJsfQhl+j0WhyDG34NRqNJsfQhl+jsUBEFUQ0w1BU/JmIbiMi12uFiEqJ6JKo+qjRhEEbfo2mMruZuRszHwbgJCSUEh/w+E4pAG34NVmBjuPXaCwQ0Q5mrmV63QaJ7MlGAFojkQSUFCu7iZl/IKJJSOg8LUVCcfFpAI8gkbxVDODfzPxiZAeh0bigDb9GY8Fq+I33tgDoCGA7gAPMvIeI2gMYycw9iKg/EjrwpxvbXwugMTM/RETFSNR5OJ+Zl0Z6MBqNDQVxd0CjyTIKATxLRN0AVADo4LDdyQCOJKLfGK/rAmiPxIxAo4kVbfg1Gg8MV08FEiqKDyBRnKcrEmtke5y+BuBmZh4VSSc1Gh/oxV2NxgUiKgHwAoBnOeEXrQtgDTMfAHAZEsJ+QMIFVNv01VEAbjBUXkFEHeIoYqLR2KFH/BpNZaoT0Qwk3Dr7kVjMTcrqPgfgPSK6HMCXSBQnARIKjBVE9DOA1wE8hUSkz0+GPO8GpEp/ajSxohd3NRqNJsfQrh6NRqPJMbTh12g0mhxDG36NRqPJMbTh12g0mhxDG36NRqPJMbTh12g0mhxDG36NRqPJMf4fhvzep8Q6eEEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light", "tags": [] }, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEWCAYAAABhffzLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dd5gUVdbA4d8h5zwkCUNQECQKiGIgKCqYP8Oa13UX113TqruCuiu6q2IOi3HNCXUNa0ARBERBJClBBBFlBCRLzul8f1T10DPTobqnqsP0eZ9nnumurq66NdN96ta5t+4VVcUYY0zuKJfuAhhjjEktC/zGGJNjLPAbY0yOscBvjDE5xgK/McbkGAv8xhiTYyzwm6whIgUicrz7+GYRecbj+zyvm21EZKuItE53OUx2EevHb1JJRI4G7gU6AvuABcB1qjrDw3sLgN+r6qel2H8+sASoqKp7k91OOojIZ8ArqlomT2ImdSqkuwAmd4hILeBD4ErgTaAScAywK53lykQiUl5V96W7HKZsslSPSaVDAFR1lKruU9UdqjpWVecCiEgbEZkgIr+KyDoReVVE6kTakIgMF5FX3Mf5IqIicqmILHXfe0ukdYHP3d8b3TTJcSKyXkQ6ha3fUES2i0hehP2WE5FbReRnEVkjIi+JSG33tY9F5Kpi688RkbPcx+1FZJy7v+9F5Nyw9V4QkSdE5CMR2Qb0K7adO3FOkiPdco90l6uItA3bxuNuObaKyBQRaSwiD4vIBhFZKCLdwrbZVETeFpG1IrJERK6J8b8zZYgFfpNKi4B9IvKiiJwsInWLvS7A3UBT4FCgOTA8ge0fDbQDBgD/EJFDI6xzrPu7jqrWUNVJwOvARWHrnA+MV9W1Ed7/W/enH9AaqAGMdF8b5b7XORiRDkBLYLSIVAfGAa8BDYHfAI+764RcANwJ1AQmh+9UVW8BvgCucstd5AQT5lzgVqABzpXUVOBr9/lbwINu2coBHwBzgINw/mbXiciJUbZryhAL/CZlVHUzTnBW4D/AWhF5X0Qaua8vVtVxqrrLDboPAsclsIvb3auIOTgBrYvH970InC8i4j6/GHg5yroXAg+q6k+quhUYBvxGRCoA7wJdRaRl2LrvqOou4BSgQFWfV9W9qvoN8DZwTti231PVKaq6X1V3eix7ce+q6iz3/e8CO1X1JTdt9AYQqvH3BPJU9Q5V3a2qP+H8T36T5H5NFrHAb1JKVReo6m9VtRlwGE7t/mEAEWkkIq+LyC8ishl4Baem6tWqsMfbcWrjXso0zV2/r4i0B9oC70dZvSnwc9jzn3Hayhqp6hZgNAeC5/nAq+7jlsARIrIx9INzYmgctq1lXsobx+qwxzsiPA/9TVoCTYuV52agkQ9lMBnOGndN2qjqQhF5AbjCXXQXztVAJ1VdLyJncCCN4ttuoyx/ESfdswp4K0aNewVO0AxpAezlQIAdBdwmIp8DVYCJ7vJlwCRVPSGJsnl9PRHLgCWqerCP2zRZwmr8JmXcxs0bRKSZ+7w5Tq34K3eVmsBWYJOIHAT8NYBirAX24+Tnw70CnIkT/F+K8f5RwF9EpJWI1MA5Wb0R1jX0I5wTwx3u8v3u8g+BQ0TkYhGp6P70jNIOEc3qCOVO1nRgi4jcJCJVRaS8iBwmIj192r7JYBb4TSptAY4Aprk9V74CvgVucF+/HegObMJJmbzjdwFUdTtOA+oUN8XR212+DKcRVHEaUaN5Dif//znO/QA7gavDtr/LLffxOA25oeVbgIE4aaAVOFcW9wCVEyj+I8DZbg+dRxN4Xwluzv8UoKt7HOuAZ4DapdmuyQ52A5cxLhF5DlihqremuyzGBMly/MZQeEfvWRzo9WJMmWWpHpPzROSfOCmn+1R1SbrLY0zQLNVjjDE5xmr8xhiTY7Iix9+gQQPNz89PdzGMMSarzJo1a52qlhhzKisCf35+PjNnzkx3MYwxJquIyM+RlluqxxhjcowFfmOMyTEW+I0xJsdY4DfGmBxjgd8YY3JM4IHfHfXvGxH50H3eSkSmichiEXlDRCoFXQZjjDEHpKLGfy2wIOz5PcBDqtoW2ABcnoIyGGOMcQUa+N1x1wfjDPeKO7Vdf5y5P8GZ/OKMIMtgMs9bs5bz8lcRuxcbY1Ig6Bu4Hgb+hjPBBkB9YGPYpBXLcSZ6LkFEhgBDAFq0aBFwMU0q3fjfOQBc3LtlnDWNMUEIrMYvIqcAa1R1VjLvV9WnVbWHqvbIyytxx7EpA75bsTndRTAmJwVZ4+8DnCYig3DmHq2FM4NQHRGp4Nb6mwG/BFgGk8F27NmX7iIYk5MCq/Gr6jBVbaaq+TjTzU1Q1QtxJp8+213tUuC9oMqQ7R4Y+z1j569KdzGMMWVMOvrx3wRcLyKLcXL+z6ahDFnhhS8LmPrTr+kuhjGmjEnJ6Jyq+hnwmfv4J6BXKvZrMp1NAmRMOtidu8YYk2Ms8BtjTI6xwG/SxrmfzxiTahb4yxhVRTU7cudNaldJ2b7271eenPQjW3ftjb+yMWWcBf4yptWwj/jrW3PTXYyMM2b+KkZ8vJARHy+Iv7IxZZwF/jLorVnL012EjPPClAIAXvlqaXoLYkwGsMBvcsL0gvXpLoIxGcMCvzHG5BgL/CZtBOvVY0w6WOA3Jofs2L2P/fuzo9eXCY4FfmNyxM49+zj0H2P45+jv0l0Uk2YW+MuQbOm/b9Jj8449ADzv9nAyucsCfxliV/Amlp179qe7CCZDWOAvQ/Zbjd/EsM8+H8Zlgb8MscBvYtlnl4TGFeScu1VEZLqIzBGR+SJyu7v8BRFZIiKz3Z+uQZUh12Rb3Fcbjz+lrGJgQoKciGUX0F9Vt4pIRWCyiHzsvvZXVX0rwH3nJPtim1isxm9CAgv86nQx2eo+rej+2CcvQBb3TSxL129PdxFMhgg0xy8i5UVkNrAGGKeq09yX7hSRuSLykIhUDrIMuSTb4n66TlRbdu7h3jEL+faXTekpQJrs2mu9eowj0MCvqvtUtSvQDOglIocBw4D2QE+gHs7k6yWIyBARmSkiM9euXRtkMTNXKQPj9t17+et/57Bp+x6m/fQr+UNHM3f5Rn/KlsU6DR/L45/9yCn/npzuoqRU41qpm//AZLaU9OpR1Y3AROAkVV2pjl3A80SZeF1Vn1bVHqraIy8vLxXFzEilGc/m1a+W8t9Zy/n3hB847+mvADht5BS/ilZqNgFXau3dZzV+4wiyV0+eiNRxH1cFTgAWikgTd5kAZwDfBlWGXBdq7C1XziKsgc8W5eiVsykhyBp/E2CiiMwFZuDk+D8EXhWRecA8oAHwrwDLkNMyvRNHKnP8Z3U/KHU7y1CZkubbs28/n32/xoYYSaMge/XMBbpFWN4/qH2aokKNl09//lOaS5Jeyzds5+dfrUdLhXKZcb/m4Ee/YNHqrdxwwiFcPeDgdBcnJ2XGJ8EE4se1W+OvlAOOvmcis37ekO5ipN3kxet829bni9byj/eSy9IuWu18Lj/+dpVv5XG2u8WuIjwK8gYuk0Y79+zL+O579hXNTj3v/JS1W3YBcMfphyW9ne9WbvalPItWb2HgQ58D8M8zDuPi3i192W5ZZjX+gOzau49Ln5vO96u2pGX/G7fvYe/+zA78qbDHerL4LhT0k/Hlj/5ddYSEgj7AvAxpx8h0FvgTsGfffhZ4rKXMXrqRSYvWcuv/5gVcqsg2bN/Nvn3e6tTPT1lC/tDRKb9Mnrc8+Buoflq7LfB95LKJC9eUWPbtL5vIHzqaY+6dUOK1Jev8/X9s3rmnyPM3Zy73dftllQX+BPzuhRmc/MgXfLcifvAPhVC/5pXdt1+Zt3wTv2zcUaTGtXvv/ojl2bBtNzs9pnpu/8CZkan/A5N8KatXqzbtCHwfNj6Nv7bv3lvk+WUvzOC6178pfL533/7CG+OWrd9B/tDRzChYX/j6jt37fC1P5+FjSyyzPH98FvgT8MUPzmXqoEe/YNP2PXHWduzxId2yeece2tz8EaeOnEyfERPoeeen/PHlWdwzZiGDHv2CQY9+weeL1hb5wK/fvpsN23cntB+/a2PxJPL1VFW+WRq9gXb15p38srHkiWT1lp1JlCzz/Lh2q+fPXCKGvTOXN2Ys9bx+pBPp/2avKHy8JkIa6JwnpxY+TkXq89HxiwPfR7azwO/B3OUb2bW3aE2lyx0laxrhQgNifbO09DnHvvd9VmLZmPmreOKzH1m8xukhcclz04u8vmH7HprXrVbqfWeKez/5njMf/5JP5hftCbJmy06W/rqdI+4aT58RE0rU9t79+peo2+x0UO1AypqIxWu2ctrIyWzeuYe//+9brg2rPQOs27qLgnXbGPDAJPo/8Jlv+31s4mLWbN7JqOnLuOlt7+nIeBdQ73wdOdWyaPUWtu7ay39n+ZOKWbnJuZqI5KFPF/myj7LMAn8cyzds57SRU2h365gSrxWE1ZBXbNzBsrDRD7fvOnBJ/POvpatJr9+WWM0dnFTP+b1aRH196a/befyzxVG/POHWbtnFT2u3MnHhGs59air705A+eeKzHwEYO391keW97hzPsfdNLHx+zeuzmRmWWnh/zgoiGXJsa9668sgSVxHrt+1m1HTvNeDSOv7BScxdvonOw8fy8lc/897sFQx5aWbh673vGk/f+z8D4Ff3c7Bvv3Luk1P57PuS+fWQid+v4dwno/+v7vvke3rdNd5TGa97/Rvyh45m5559US/Tvv1lE1/+uI77x0YOugMf+pwrXp4Z8bVkvPLVzzFf37nH35TSLxt3RLyizFZlOvDf98lC8oeO5r5PFia9jVj9v68e9Q1rt+xiw7bdHDViAsfceyAADXfz5gD/eG9+0vv3qlKFov/K9dt2c3DDGlHXP/a+idw75ntP2+5556f0f2ASl70wg+lL1rOtWJ43Eas3R069PDRuUZFccLjw4FUQ5yT6wZwVnP3kVO76aEHMk9rNgw7lDy/N4szHv+Qn936HLxevo/s/xzHsnXksXLWZ/KGj+fOrXzPr5w3kDx3N/BUHGqNXbdqZUIokEWO/W03+0NHkDx3N3mKBO3/oaK5/czbTC9bz2+dn8Mn8VbzgNs6H/k4F67Zx2fMzmF6wnqNGTGDhqvhtUgtXbWZDlApGKJXT/u9jGDFmQcR1Tvn3ZC74z7SIr4VMWfxriWXR/ufxTI6wrXB/eKl0J5lfNu5gadhNf31GTKDPiJKN1cnasXtf4f84FZ0ciivT/fgfm/hj4e+/ntg+qW1s2hE9rzrvl00cc+8E/nVGp8Jlx9w7gTevOLLIehXCxsqZWbCeTs1qU7lC+aTKE83uvfuLVMY2bt/Nsg3B3K26YuNOtu/eSrcWdT2/55P5q7ji5VlFli3fsIPr35xNXs3KPDXpJx4Z/wMFIwaXeO/GsP9BrJNZOC93K3/ujl3T/4FJ3DjwkCK11ZMe/gKA0fNWMnreSgAGP1pyNM+e+fWoWaUiZz0xhRcv60XrvAPlW7NlJ31GTGCP27uq+LG9nWTa472wnHr437T1zR8x7OT23P3xgYrOqs07C48llpMe/oJyAmP/chzzftlI41pVObJNfaYUu+lr1PRlSZU5mnOenMqj53cjr0Zlygl0b1mXzTv2sG7rbto1rlm43rlPTWX6kvWFf8M5y2KnUEPtcckKBflIn8fSCFVGHr+we+GyU0dOZvJN/di/H1rUT016tkwH/tL6ZumGuLX1nXuKNt4uW7+Dd4rllRes3My9YxbSv31Dzn5yKhce0YI7z3ROFvv2Kw+NW8QNAw9B3OEq12/bzYokLiufdNMhAOu37ykSIPx04TPTWLd1FzNuOZ68mt6mUyge9ME59uJ/K3ACcoXywlFtGgBF+36/PmMZtww+lJpVKjLtp9i1vkRES1HE89e35nJK5yYsW7+D/3yxhLvPOlAJ6HVn0VTKo+N/4MFxi/jib/1Ysm4bN/x3TqnKHEl40E/UfnVSTyEfXn00Fz4Tuxbvh2tGHWjX6NO2Pt+t2MyG7XsYfc3RdGzqtMNMX+JcGZzx2BT+9+c+nra7ZN02Rs9dwR+Pa0OF8tGTG/OWb2Lp+u0M7tykxGv5Q0dTvdKBStqi1VuoWrE8zeuVLkD/6dWvizw/+h4nW/DaH45gzz7luEOCHZG4TKd6vJq9bGORHhP79yvPT1nCmY9/6en9780uGrzu+6RoCmXFpp08/tmPnO32bnh1mpMiKFi3jTY3f8TIiYtpNeyjwvWPuOtTTvn3ZLbsSiyl8nhY4P980Vp+e1R+xPVipUCuHvUNqsqgR76Imh9ft9XpubHNY/kSGRwsf+hoLnluOhf8Zxrbdu3l+1VbePWroimVTsPH8vsXZxYONZ1Os37eUNgddtT0pTw49vuo3QkfHOecXM54bEqJxvhMlI75CqYs/pUN7ndx8KOTmbJ4XZGb8GbHqemH63f/Z9w/dhE3vT2P/fuVL35wer7NXb6R/KGjeXPmMpZv2M6pIyfz59e+LvxeTP2xaIViW1gX1IEPfV6Y0h3z7UoGPjSJjdt3x+xxFuLl5rUL/jONS5+bHjUl6pcyXePv374hExauoWPTWjHXO+OxKXRuVpv3rzoagEfG/8Aj43/wvJ9kLyujNSLuCbvxakcpGqni5cMj+WDOCj5wA/41o76hXrVKUdf98sdfWbR6C/3bN4xZo0p2DoBrRn3D+Ag3CAF8umB1xOVenfTw5/FXSsKjExbz7OQlNK1TNeo6vybRWJ+rIl1x3P1x5HaGaD5dsJrnpizhX6OLvu9vb80tse7OPfs8tQ+EV56Of/Bz1m3dxa2DD+XQJrV4+NNFnNqlKZccmV+4zuxlG+O2gYQ74q7xzBs+kJpVKnp+TyLKdI0/VPOav2Iz05es540ZS6P2cpjrNrBMWrQ2oaCfrC8Xr6N8hHHyiw8xUKEUY+k//Gnpj+OiZ6N/WG9+dx5DXp6VdHrhhS8LYr4eLej7YWGA/cm37d7HD2tsgLw3hvQOZLtPTUpstNlNO/aUCPrRtP/7GLYmeKUdugL+1+gFXPjMNGYUlEwRn/FY4pWf16cvY9WmYGr+ZTrwhzv3qanc9PY8Bj16oKFr6669fOQ23AFMWbyOS1N0CX7BM9OKpGZCDr7l4yLPV6bg7tbSenbyEsDp9hnqqRBid1FmpioVg//qS45PsTZl8TqenPSjpy7Tkdz50QJ63z0+kBsry3SqJ1LIWbjKGbpVRLjs+enMKDiQmws1IKVLpA/IpwuCq/X6KdqH20ZMyEzXDDjYc3dekxy/Gsaveu1rRl9zjC/bCgly6sUqIjJdROaIyHwRud1d3kpEponIYhF5Q0SiJ5FLKVplc6h7p2J40Af43+zod3maxITu4MyUsXKODbiXRLbxawypeJ69tEdK9lOW+TV8dbggr/d2Af1VtQvQFThJRHoD9wAPqWpbYANweYBliOiNmZH7IifThdJEdv2bTlfF/RmS6rnoiOh3MWcSv/uNR3PPmOS7fSbimIPthBvukEY1mHHL8Qm9J4ivUGCBXx2hFq6K7o8C/YG33OUv4ky4Hojf9Gye0Pp7PA5jnGrn9mjG0xcfnu5iJCXVNf4/92vD0W0blFg+sGPjlJYjUa8P6c27fzoq5jrzhg/0dZ9X9Wvr6/YiyZY0/8mHpebzMfYvx1GxfPr/KIG28IhIeRGZDawBxgE/AhtVNdRsvhyIOAu2iAwRkZkiMnPt2rVJ7f/kTk1Y+M+Tknpvpvjrie249+wuDOzYmFd/f0Rg+zn+0IZcfnQrX7c5ceGapMbDL01PprrVKvFKsb/TEa3qJb29SApGDPa9Zt67df24d0L73bXvxhPbJbR+l2aJD2qX/hDnzSO/KTE9uO9C99UkWsH0epNkIgIN/Kq6T1W7As2AXoDncRNU9WlV7aGqPfLykr9cjNRlMhl/OMbfoOjF/ed04c9htbJaPnzxf18suI+8oBuP/KYrz1zak0vD+h374bIXZnDqyMRvAjqyTX3P63ZpXoe7z+rEgPYNAefeDaDIFdKj5wfzpW7XqGb8lcJMvLFv0vuKd9J/ImwIgGhOiXBnaiIeOLdrwu8p3rOnbcMavH3lUVSKcd9HOhQf6yoIobkMEg1Jse6lSVbcoxWRB0SkY2l2oqobgYnAkUAdEQn1JmoGBNqiWs6na81bBnfwZTuJOPvwZkWea1g/pb8mWFsLObhRDfq2O3AirVe9Eqd3dS66wscJOaFDo6S2H+7tK4/kgXO6JPy+/+veLP5KrpqVK3B+rxY8+9ueFIwYXDhWzvGHHih/o1pVEi6DF4mcUBK5QqhaseQ4Tn0ipK9Cpg7rz8mdDgT1xrWq8IF7M2K4Q8JOVOP+cqzn8nRpXoeCEYOT6gJa/Nv36fXHcXjLutx6yqEJbyvbVavkhL36NSrz7/O7ec71B5Eu8/KfXAA87fbE+aOIeLreE5E8EanjPq4KnOBuayJwtrvapcB7iRfbO58q/BkhPF3+51LkZ1+4rBdHtnZq1VF7Pp3cnhm3HF+qMUMOb1mP/zvcWxAPT8md0S1i9i+iclH+weXKCU9dfDhfDRvgeVuJSrTh2us9DdNuSazMxTfbvknNIpWESA5O4GrlPY9j40QSHrTCT36Z0tsrFapWLE//9g0ZevKBhMepXZp6TuH4lbUIFzfwq+ozqtoHuATIB+aKyGsi0i/OW5sAE0VkLjADGKeqHwI3AdeLyGKgPvBsaQ4gnrJ0E0nxQFOaXHO5cpG3GaKq5NWszIu/60VbjyNilkay/6ZY7WQndmxM49qxa/ulObE1qxt9WIZIKkeoyUeSaEqvXnUnFRBqrG2TVyPiCT0dHayiff+CKMtBMYbJSKdJf+vLc7/tSRWP///i6qYj1QNOIy1Ofr49sA6YgxO8X4/2HlWdq6rdVLWzqh6mqne4y39S1V6q2lZVz1HVknO1pUA23lHqZ5lDKbBoFa/w5RemqCvk60N6M+GG46K+PuvWkpfGpT2xj7wgcrrmvrM7x31voo2t0QLTK5eXrtE+FFBCV0pndPV+xZSIRP/WsVavWcX/e0evGRB8L6VkNKxZulTj6V2b+lSSA7zk+B8CFgKDgLtU9XBVvUdVTwWCbwpPk7vP6sT9SeSng+TnuepA4C+60UMa1Sixr14ee8U0LkUuXRB6t65fmKMPz3PfcMIhPHXx4dSvUfLS+LwEu+yGhK6WogXvc3o0Z8x1x5Af8Pjofz+lA0cfXDJ/P//2ExPeVtuGNSgYMZhOzWrTpXkdFtyRfI+2awccnPR7vfDSgB9qcL2ybxsa1Ihf6/WSPXoqhd2iX/v9EUy+KV5iJL60pHqAuUBXVb1CVYsPZNPL9xKlSKwgWrdaRc7v1YIZHoZwePvK6H2vh52c3OQv0RziTkxx15md4qzp6N++IRf1LlpbDzUShj5Lxa8iQnd0hp8QKpTz1qhXmiuD4r0qFoTl/K/q35YTo/TDP7yl98lgEtW+cS0Kfg1mMpuQaF/pilF6vZzaxXvtr2ql8hSMGMyJHRNvqA9N2/nkRQcCpZ/hp0GEk3hxE2/sy2t/OIKbTmrPzFtPiLu+l3aDaJ8jgGMinICTdcWxrTmqbQOaJTHvdeUK5SgYMZgz3Jp+uhp3L1LVIp2xRWQ8gKqmfs6wJLwT4caY0MBikYR6gXQ8KPZwzuBcsvbMr8t5PUrWPM/q3szXy7RaVSpSMGIwF3gMsP3a5RWZHQwo/CD+4ZjWAHRuVqfI6385wanptQibaKJaJW+5yWgNrdHccXrHIr+jiZViSKSBNTQ8dyK9hhLx9pVHMumvfRN6T6KVuQfPTc1VaOPaVSgYMZiTArqxKVq+O/yzfVCdqoWT8XhR2slLRl4Qv0usF5Nv6sewQd56LYX3lDq1S1Mn3el2+734yJaAc4+H36IGfnesnXpAAxGpKyL13J98otx0lam6t6hLwYjBRboo3vlR9GFaQ7nYQ4r1fHg9wjCzdatV4r9/PIp7IuSE82pW5pHfdGPRv04uMtVaKoy8oBsXHtEy6utHtW1AwYjBJWpeJx3WhIIRg6le+UAO1mujWe2qieW8Lzkyn4IRg4uMWx7usQu6xw90CaS/BrldHhvW8v+GGHB6MbWsX73E8jZ5JZeFRDupRTvXRbsSiCXRFGG0G96C7ifRJq+656vZSJrXq8aPdw1i5AXdWHL3IHokeDWYyOc32gm+X7u8hGr5D5934N6If5/fjd6t6xd+3w5vWY+CEYNpUtv/RutYn6IrgFk4Dbpfu49n4XS/HOl7SVLgTI/dBP9xqtNnv/iZtnuEOyu9fFgqVShXGHS86plfuhTGKZ2bJlwDj8brF75d45p0aV4n6utL7h6U0H4Hd27CWcVq5z/eVXQbicS00P/A6+fAi0jDQxRX/KoqXCq6G//hWOfqLlJabMx1RUd9LBgxmDeKzRmdKn70wCtfTjilc1NEhGcCHCAu0gke4KmLE9tn+8bxswpBiBr4VfURVW0F3KiqrcJ+uqhqVgZ+r+NxRPunRhLUHX/nRkgdRRP0wF7hX8hINxeFu/TI6FcZfn2xHwu7JE+kNtuqQXUKRgwucSVXGq/8/gja5FXnppOit+fEuuciWl9uP88HPfOdmmP9CA2k6Qo84T693rmZzO+ednWqVSrSvuF3T6JuLUqe0BONB/kNqvP8ZT2TaswvjVipnv7uw19E5KziPykqn6+C7tMfb0yYWJf82aJhrcqc28OphRfPlatGzt1+8bd+fHh1yTtJkxU+KXYmDHg1/oa+XNm3TcTXCkYMjnkfRKzGRi9qVPYezEobV7fvTn4a0FhC38tkyhfvbuLwHjEnJfm3HhjlLvZ3/9Qn4h3SierXrmGR1GoqxNrbccAE4NQIrynwTiAlymLxPrj//eNR/LR2a+Gk69no3+d3o3OzOgzq1ITerevztjvufizN61UjuU6X0V3WJ5+Xpv4csYtntqhRuUKpKyPxblALF+9u3nh27dkff6UkFO9afMERLdiy09v0h1UqlmdnjHJd0KslH81bBcBdZyXXfnDfOV04r2A9l79Yci7eIO5HSIWopVbV29zfl6WuOOn391O8jcnz5dD+bNhedNLsfXEif73qlahX3d+RIlMtlAbr265hidcSuVQvbf/4207tyG2nlmoIqeY6bNUAABqDSURBVLSLFfK9/iVTeSNis3rxGxlrVK6Q8Jy1ob9DqDdmIg28d53ZiT+9+nXU10NXW41qVU6qYRwAhQGHRq715zc4cBUfVG+xIHi5gauyiFwgIjeLyD9CP6koXBA6NImd0/xdn3xP22lapyodmxYdtih085MfsnGoCSX+VU8Tt4b63G97Bl+gALRu4F+6Lta/2K/BBcPF+9/Ea2iuVaVi3Brut1Fy1c1jnDRCx5rIFckdp3ekW4s6cTtNhP6Mkbr4v/g7f29D8mNgw1Txcgp8Dzgd2AtsC/vJStceH/uOxNIEXK810EiNQkH643GR889BiPflDdW6gghsQWuTV52Pr4s/92kjj91FY33WypeTIl39UsFL2J03PPFGyA+uOpr//ankQG9Th/Vn7F+OPRCcE8gkXXJkPu9G2GZxoW1HOul56ffftmENanhM52RTG56XwN9MVc9T1XtV9YHQT+AlC0isgHPDCYeUattVKpbn9tOc4F+/evRbzN/9Ux8GdYrd0ORnWGzXOPhB1qDolytad9RO7mQeXm8KyxQFIwYz/oa+VK4Qv9yNPfa7jnfu6xqja2wQ6lcPpr2kU7PaEdtimtSuyiGNahZ2Ow4ibVWvWiXaN67JPf8XOX10+2kdY54APr3+uMIG4rlxZkBrk5ea75kfvJzKvhSRTqo6L/DSpECs79rVPoxPEhrrPl4rfSK1m9Ly8/sUL1iF9hVtYKr7z+7C5Ue3omFAY+Qno2aVClEbE6ffMoCN2/cktkGPf/B4J/fQ3zr8LurSiDfKauckZtjyQ/Ecv58qlC/HmOuizz1w6VH5XOrOjBXuq2EDSrRVxBs11a/7ZlLBS43/aGCWiHwvInNFZJ471HJWCuquzUSFejCc0KFRkYlRQvzMhAQ59Hn4hDAanuiJUv6qlcpHvBEuna7uH72ffcOaVRLu99/xIG8BNF5asWmdqnQ6qDZ3J9kbpbh4wwIn0jXUT9EGDPRq2s0DYo7qmozGtaukZDjydPES+E8GDgYG4nTtPIXIXTyzQudmdWIOrFZaXj+7odX+r3szXrisZCOTn7dp+3EJHWrUrFGpaHCINhZR9tR9nFvj/XTbqd56hsWrIFYsX44Prj465uxbXWLcFezVNe6V7kEJzi/glwNzQyT3/ka1qhSO6mq8iXUDV6j7y5YoP1nr8JZ1+V8pZhXyIn5KxPmUR/rydzqodkLzzsbjR4V/wo19KRgxuMTlbPXwE0H2TXEAJF/TjMZLO4AjudNjaLC5kRd0S7pverjKPt99nugYOaERYbNxjoxsFes//pr7exYwkwNj9YSexyQizUVkooh8JyLzReRad/lwEflFRGa7P4kN4OKTrs3rMP3m4KbliydUuwld5n4WNhF3aQd0u6T4sAnuvuLdWZyMutUrUafagdxnNn5596dpGsBk03lvX3kUH11zDKd0bpr0rE7hzu3RnE4H1S75uYnj+csid8kNjU1UL0YHh3CFQ4QntPfkXXFc65hDbOSCWDdwneL+bpXktvcCN6jq1yJSE6edYJz72kOqen+S2y0TCmv87qm3adgImM2TaMyrUrFc4R2MxSdOCWXe/WokLO7QxrWY+tOvRZZl030I6Zr+Ndm/UJWK5enQ1L8xdvJqVuaDJIbU6BfhJj6AYYPac1b3gzz3ciltjj9Rw07OvYnei4sa+EUkZrVTVaPfLue8vhJY6T7eIiILyODhnNs39m/gLi9CwUZ8yoZXLF8u6q3rWuzqwm/DT+vI7R/Mp3vLuoz51rk9PnvCfvquUrLo3JiQiuXLcZjHBm440Ojcu5X/486byGI148fqq69A/xivF+GO4d8NmAb0Aa4SkUtwUkY3qOqGCO8ZAgwBaNEiNXO+plK83i+JOrJ1fcZ+t9rZdrE4VniSCSjQtGtck9f+UHSugmwKatFGyAyaXyf9bFe1Unk+vf7YpGarSrWpw/pTKdmhHzJIrGGZ+8X4SSTo1wDeBq5T1c3AE0AboCvOFUHEE4yqPq2qPVS1R15e6WbWyUTDT+1A//YNOdKn2XW8NPKlIhif0KERxx6Sx40D28VfOUMc3Kgmn1x3bNKjNyYifJTSxy9K7eQ8maxtw5q+tFf4IdYct01qV83qgQFDPHXcFZHDgA5A4V03qvqSh/dVxAn6r6rqO+77Voe9/h/gwwTLnNG8Jg1a59UoMl5NaUdODO+DXTzAn9a1KZ8uWM21A0p3Z7IX1StX4CWfx0BJhXaNawYyqXVx4SmQTLufwTj3BPjdyykTxQ38InIb0Bcn8H+E069/MhAz8IvTuvcssEBVHwxb3sTN/wOcCXybVMkzXKov4mPVlmpUrpC1g6IZk0qNMuiO8iB5qfGfDXQBvlHVy0SkEfCKh/f1AS4G5onIbHfZzcD5ItIVp3JcgDPFY1pkX8dDb7KwR6UxJoW8BP4dqrpfRPa6N3WtgfjzaqjqZCJXfD9KsIzGpERp023h5vxjIHtTOSBTBvjvH9MzV69JnJdk1kwRqQP8B+fmra+B7J1CKgdkU4+asqp2tYplohEwET3zs2uSoTWbd3LuU1P57Ps1pdrOm1ccyVUx5lXORHFr/Kr6J/fhkyIyBqilqlk7SFumsvSMyUYdm9Zi/orNKd3nC1HuGE7Urr37mb5kPef2KN3EoL1a1Stx02Smi3UDV0tgo6pucp/3A84AfhaRhaq6O9p7TXrZSSQ59ndLXConCX99SG++Xroh4rSfJjGxUj1vAtUB3MbY/wJLcRp6Hw++aMaYTPT+VX241oe5KxLVu3V9/tQ3u1IqmSrW6bqqqq5wH18EPKeqD4hIOWB2jPfltGwcpMyYRHRuVqdwIDaTnWLV+MObCPsD4wFUNbe6KiQpnYOUWeOuMSaWWDX+CSLyJs6wCnWBCeDcgAWUifx+JlXO/RyZMJOOK5tUKANjsHhxYsdG6S6CSbNYgf864DygCXC0qoYmHm0M3BJ0wXKNBev0G35qBz6YsyL+illszm0Ds26ie+O/WOPxK/B6hOXfBFqiNEl34E3VWOQmulzod1+7auwJw01uyI1r2yxgYd9ko0wZUdMkxgK/K90NolbhN9no/nM6p7sIJgmeAr+IVBWR7BlgPY2Sjt8W+E0WalgzN0azLGu8DMt8KnA/UAlo5d7MdYeqnhZ04bJZwhcQ1gUzIzx18eFUtfSFKeO83G89HOgFfAagqrNFJNkJ2E0U1uiWGU5MwSxcxqSbl1TPntB4PWHKRGIifBhey7EbY3KFl8A/X0QuAMqLyMEi8m/gy4DLZZJw8mFWWzXGxOcl8F8NdAR2Aa8Bm3Bu7opJRJqLyEQR+U5E5ovIte7yeiIyTkR+cH/bxKM+KZeCOWONKSuqV67AeT2ak1+/WrqLknJxA7+qblfVW1S1p/tzq6ru9LDtvcANqtoB6A38WUQ6AEOB8ap6MM74P0NLcwB+SXd3Tj/UquI02eTCZNHZbGAHGzIhE9SrXol7zu5MjyybQMYPXnr1jAPOUdWN7vO6wOuqemKs97kTqq90H28RkQXAQcDpOJO3A7yI02h8U5LlzzilaSsYdnJ7+rRtkPT7bxncgfz61Tn+UAssmezpS3qkuwi+6pVfj/J2tZlVvPTqaRAK+gCqukFEEpoJQUTygW7ANKCRe1IAWAVEjFIiMgQYAtCiRYtEdpcZkvgeXHFcm1LtskblCqXehjGJ+k8ZO5HlAi85gf0iUhh53Zm5PNdrRaQG8DZwnaoWmaPNHQ8o4rZU9WlV7aGqPfLy8rzuzhiTYrWrVaR2NeuOnE281PhvASaLyCSceuwxuDXxeESkIk7Qf1VV33EXrxaRJqq60h3iuXQzHfvEunMaY3KFl8bdMUB34A2c0ToPV9VP4r1PnJlIngUWqOqDYS+9D1zqPr4UeC/RQvvFgr0xJhd5nSl5H07NvArQQURQ1c/jvKcPcDEwT0RCUzXeDIwA3hSRy4GfgXMTL7b/ykKvHmOM8cJLr57fA9cCzXDm2u0NTMWZjjEqVZ1M9CbOAYkV0xhjjF+8NO5eC/QEflbVfji9czbGfksus/yRMSazeQn8O0M3bIlIZVVdCNgQzXFY5sgYk6m85PiXi0gd4H/AOBHZgJObL1OsodcYkyuiBn4RaaWqS1T1THfRcBGZCNQGxqSkdMYYY3wXK9XzFoCIjA8tUNVJqvq+qu4OvGTGGGMCESvVU05EbgYOEZHri79YrG9+VgrP7lh3TmNMrohV4/8NTv/9CkDNCD/GGGOyUNQav6p+LyL3AUtVdVQKy5TVrJHYGJPpYnbnVNX9wA0pKkuZIpY7MsZkKC/9+D8VkRvdGbXqhX4CL1mKWU3dGJMrvPTjP8/9/eewZQq09r84xhhjghY38Ktqq1QUxBhjTGp4GaTtkkjLVfUl/4uTPpaSN8bkCi+pnp5hj6vgjKz5NZD1gV8tsW+MyUFeUj1Xhz93x+15PbASZTk7lRhjMp2XXj3FbQPi5v1F5DkRWSMi34YtGy4iv4jIbPdnUBL7D4TflX/LHBljMpWXHP8HHKjIlgcOBd70sO0XgJGUTAk9pKr3J1DGtDr78GY0r1st3cUwxhjfeMnxhwfpvTgTsiyP9yZV/VxE8pMsV8a4/5wu6S6CMcb4ystk65NUdRLwLdAQaFTKfV4lInPdVFDdUm7LGGNMgqIGfhH5UEQOcx83wQn8vwNeFpHrktzfE0AboCuwEnggxv6HiMhMEZm5du3aJHfnnXXnNMbkilg1/laqGmqYvQwYp6qnAkfgnAASpqqrVXWfOwbQf4BeMdZ9WlV7qGqPvLy8ZHZnjDEmgliBf0/Y4wHARwCqugXYn8zO3CuHkDNxriLSJrwnj1+9euzWAGNMpovVuLtMRK4GlgPdcadbFJGqQMV4GxaRUUBfoIGILAduA/qKSFecXkIFwBWlKXwms9SRMSZTxQr8lwN3AMcD56nqRnd5b+D5eBtW1fMjLH424RIaY4zxVayJWNYAf4ywfCIwMchCGWOMCU4yd+6WSZaaMcbkCgv8xhiTYyzwG2NMjvEyVk8VnIbejjjDMgOgqkn15c9UvnXntPE5jTEZzkuN/2WgMXAiMAloBmwJslBlgdj4nMaYDOUl8LdV1b8D21T1RWAwzt27xhhjspCXwB+6g3ejO3ZPbZzB2owxxmQhL8MyP+2Oovl34H2gBvCPQEuVBtad0xiTK7xMvfiM+3AS0DrY4mQHO0kYY7JZ3FSPiDQSkWdF5GP3eQcRuTz4oqWWDa5mjMkVXnL8LwCfAE3d54uAZMfjL/PsBGKMyXReAn8DVX0TdyhmVd0L7Au0VGWApYOMMZnKS+DfJiL1cSdcF5HewKZAS5UiVjs3xuQiL716rsfpzdNGRKYAecDZgZbKGGNMYLz06vlaRI4D2gECfK+qe+K8LeskkpqxLI4xJpvFmmy9p4g0hsK8/uHAncADIlIv3oZF5DkRWSMi34Ytqyci40TkB/d3XR+OwReW9jHG5IpYOf6ngN0AInIsMAJ4CSe//7SHbb8AnFRs2VBgvKoeDIx3nxtjjEmhWIG/vKqudx+fBzytqm+74/a0jbdhVf0cWF9s8enAi+7jF4EzEixvxrMrB2NMposZ+EUk1AYwAJgQ9pqXRuFIGqnqSvfxKqBRtBVFZIiIzBSRmWvXrk1yd8YYY4qLFfhHAZNE5D1gB/AFgIi0xYfunKqqEH3welV9WlV7qGqPvLy80u7OV2Kd9I0xWSzWZOt3ish4oAkw1g3U4Jwsrk5yf6tFpImqrhSRJsCaJLfji/BJUyyWG2NyRcyUjap+FWHZolLs733gUpyG4kuB90qxLWOMMUkIbM5dERkFTAXaichyd2C3EcAJIvIDcLz7PCNYo6wxJlck20gbl6qeH+WlAUHt0xhjTHyB1fhzlU22bozJdBb4k+ClHdh6/hhjMpUFfmOMyTE5HfjDG3Stgm6MyRU5HfjDWa8eY0yusMBvjDE5xgJ/EiwtZIzJZhb4fWYpI2NMprPAHxC7KDDGZCoL/MYYk2Ms8Lssb2+MyRU5HfjD0/FecvM2HIMxpizI6cCfLLEMvjEmi1ngN8aYHGOB3xhjcowF/oBYY7ExJlMFNhFLLCJSAGwB9gF7VbVHOsqRKLs5yxhTFqQl8Lv6qeq6NO4/eVabN8ZkMUv1GGNMjklX4FdgrIjMEpEhkVYQkSEiMlNEZq5duzaYQljuxhiTg9IV+I9W1e7AycCfReTY4iuo6tOq2kNVe+Tl5aW+hMYYU0alJfCr6i/u7zXAu0CvdJQjCHYRYYzJdCkP/CJSXURqhh4DA4FvU12O0vA22XrgxTDGmKSko1dPI+BdcSJjBeA1VR2ThnIYY0xOSnngV9WfgC6p3q8xxhiHdec0xpgcY4HfGGNyTE4H/kQ74FiPHWNMWZDTgT9ZsXrs2GQtxphMZ4E/IDZZizEmU1ngN8aYHGOB3xhjcowF/gRY/t4YUxZY4E+C5e+NMdnMAr8xxuSYnA78QfTLt77+xphMl9OBP0g2OqcxJlNZ4DfGmBxjgT8BlsYxxpQFFviTYGkcY0w2s8BvjDE5Ji2BX0ROEpHvRWSxiAxNRxmMMSZXpWPO3fLAY8DJQAfgfBHpkOpyFFepgj9/iia1qwBwYsfGvmzPGGP8lo45d3sBi90pGBGR14HTge9SXZD61SsVPn7iou6e31e1YvmorzWsVYW5wwdSs3I6/rTGGBNfOqLTQcCysOfLgSOKryQiQ4AhAC1atAikIHWrV2Le8IFUr1SBcuXit9hWr1yBoSe354QOjWKuV6tKRb+KaIwxvsvYxl1VfVpVe6hqj7y8vMD2U7NKRU9BP+SPx7WhTV6NwMpjjDFBS0fg/wVoHva8mbvMGGNMCqQj8M8ADhaRViJSCfgN8H4aymGMMTkp5Tl+Vd0rIlcBnwDlgedUdX6qy2GMMbkqLV1PVPUj4KN07NsYY3JdxjbuGmOMCYYFfmOMyTEW+I0xJsdY4DfGmBwjmgWDzIvIWuBnHzfZAFjn4/bSwY4hM9gxZAY7hshaqmqJO2CzIvD7TURmqmqPdJejNOwYMoMdQ2awY0iMpXqMMSbHWOA3xpgck6uB/+l0F8AHdgyZwY4hM9gxJCAnc/zGGJPLcrXGb4wxOcsCvzHG5JgyEfhFpLmITBSR70Rkvohc6y6vJyLjROQH93ddd7mIyKPuZO9zRaR72LbudbexwF3H+ywtqT2G9iIyVUR2iciNxbaVlsns/TqGaNvJpmMI2155EflGRD7MxmMQkToi8paILHS/E0dm4TH8xd3GtyIySkSqZOgxXOjGo3ki8qWIdAnblr/faVXN+h+gCdDdfVwTWIQzkfu9wFB3+VDgHvfxIOBjQIDewDR3+VHAFJzhossDU4G+GXoMDYGewJ3AjWHbKQ/8CLQGKgFzgA5ZdgwRt5NNxxC2veuB14APM/j7EPUYgBeB37uPKwF1sukYcKZ6XQJUdZ+/Cfw2Q4/hKKCu+/hkDsQl37/TZaLGr6orVfVr9/EWYAHOP/x0nA8u7u8z3MenAy+p4yugjog0ARSogvPHrQxUBFZn4jGo6hpVnQHsKbapwsnsVXU3EJrMPnB+HUOM7WTNMQCISDNgMPBMCopeyK9jEJHawLHAs+56u1V1YzYdg6sCUFVEKgDVgBUBFx+3TIkew5equsFd/hXO7IQQwHe6TAT+cCKSD3QDpgGNVHWl+9IqIDRLeqQJ3w9S1anARGCl+/OJqi5IQbGL8HgM0UQ8Np+LGFcpjyHadlLKh2N4GPgbsD+I8nlRymNoBawFnnfTVc+ISPWgyhpNaY5BVX8B7geW4nynN6nq2MAKG0USx3A5TlYCAvhOl6nALyI1gLeB61R1c/hr6lwzxey7KiJtgUNxzrQHAf1F5JiAihutDKU6hkzg1zHE2k7QfPgsnQKsUdVZwZUyNh/+DxWA7sATqtoN2IaTmkgZH/4PdXFqx62ApkB1EbkooOJGK0NCxyAi/XAC/01BlanMBH4RqYjzx31VVd9xF692Uzi4v9e4y6NN+H4m8JWqblXVrThn3JQ0ZrllTOQYoknrZPY+HUO07aSET8fQBzhNRApwLs37i8grARW5BJ+OYTmwXFVDV1tv4ZwIUsKnYzgeWKKqa1V1D/AOTi49JRI9BhHpjJMaPF1Vf3UX+/6dLhOBX0QEJw+5QFUfDHvpfeBS9/GlwHthyy8RR2+cy7+VOJeDx4lIBfcfdhxOXi4TjyGatE1m79cxxNhO4Pw6BlUdpqrNVDUf538wQVVTUtP08RhWActEpJ27aADwnc/FjcjH78NSoLeIVHO3OYAM/U6LSAucE9PFqroobH3/v9OlaRnOlB/gaJzLpbnAbPdnEFAfGA/8AHwK1HPXF+AxnJbyeUAPPdB6/hTOB+M74MEMPobGODWyzcBG93Et97VBOD0IfgRuybZjiLadbDqGYtvsS2p79fj5WeoKzHS39T/cXidZdgy3AwuBb4GXgcoZegzPABvC1p0Zti1fv9M2ZIMxxuSYMpHqMcYY450FfmOMyTEW+I0xJsdY4DfGmBxjgd8YY3KMBX5jihGRfSIy2x1RcY6I3CAiMb8rIpIvIhekqozGlIYFfmNK2qGqXVW1I3ACzkiJt8V5Tz5ggd9kBevHb0wxIrJVVWuEPW+Nc/dkA6Alzk1AocHKrlLVL0XkK5xxnpbgjLj4KDAC5+atysBjqvpUyg7CmBgs8BtTTPHA7y7bCLQDtgD7VXWniBwMjFLVHiLSF2cc+FPc9YcADVX1XyJSGWeeh3NUdUlKD8aYCCqkuwDGZJmKwEgR6QrsAw6Jst5AoLOInO0+rw0cjHNFYExaWeA3Jg431bMPZxTF23Am5+mC00a2M9rbgKtV9ZOUFNKYBFjjrjExiEge8CQwUp28aG1gparuBy7GGdgPnBRQzbC3fgJc6Y7yiogcko5JTIyJxGr8xpRUVURm46R19uI05oaG1X0ceFtELgHG4ExOAs4IjPtEZA7wAvAITk+fr93heddyYOpPY9LKGneNMSbHWKrHGGNyjAV+Y4zJMRb4jTEmx1jgN8aYHGOB3xhjcowFfmOMyTEW+I0xJsf8P9b9miXYleqfAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light", "tags": [] }, "output_type": "display_data" } ], "source": [ "#making the url, opening the dataset:\n", "url = 'https://pricaimcit.services.brown.edu/erddap/tabledap/erdGtsppBest.nc?station_id%2Clongitude%2Clatitude%2Ctime%2Cdepth%2Ctemperature%2Csalinity&longitude=-71.321&latitude=41.579&time%3E=1985-03-31T13%3A15%3A00Z&time%3C=2020-03-30T10%3A08%3A00Z'\n", "r = requests.get(url, allow_redirects=True)\n", "open('test.nc', 'wb').write(r.content)\n", "nc = NetCDFFile('test.nc')\n", "nc\n", "#defining dataset variables:\n", "time = nc.variables['time'][:]\n", "time = nc.variables['time'][:]\n", "mapvar = nc.variables['temperature'][:]\n", "salinity = nc.variables['salinity'][:]\n", "depth = nc.variables['depth'][:]\n", "lat = nc.variables['latitude'][:]\n", "#making a list of times for graphing:\n", "timearray = []\n", "for i in time:\n", " timearray.append(datetime.fromtimestamp(i))\n", "plt.plot(timearray, mapvar)\n", "plt.xlabel('Date')\n", "plt.ylabel('Sea Surface Temperature (Celsius)')\n", "plt.title('Temp over time')\n", "plt.show()\n", "plt.plot(timearray, salinity)\n", "plt.xlabel('Date')\n", "plt.ylabel('Sea Surface Salinity')\n", "plt.title('Salinity over time')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": { "colab": {}, "colab_type": "code", "id": "Bb-lHx0vdjRR" }, "outputs": [], "source": [] } ], "metadata": { "colab": { "authorship_tag": "ABX9TyN7cEddVmr2Qj86AL+GbxTw", "collapsed_sections": [], "include_colab_link": true, "name": "GTSPP Data", "provenance": [] }, "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.7" } }, "nbformat": 4, "nbformat_minor": 2 }