{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example of using the Radiation class, retrieving the data via ftp server" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import BCO\n", "from BCO.Instruments import Radiation\n", "import matplotlib.pyplot as plt\n", "from matplotlib.dates import HourLocator, DateFormatter, DayLocator" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Set up ftp-usage:" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Successfully loaded username and password\n" ] } ], "source": [ "BCO.settings.set_ftp(True)\n", "BCO.settings.path_to_ftp_file(\"/home/tobias/Documents/ftp_access.txt\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Initiate a new instance with data from the 1st January 2018 to the 3rd January 2018:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Downloading Radiation__Deebles_Point__DownwellingRadiation__1s__20180101.nc.bz2\n", "Downloading Radiation__Deebles_Point__DownwellingRadiation__1s__20180102.nc.bz2\n", "Downloading Radiation__Deebles_Point__DownwellingRadiation__1s__20180103.nc.bz2\n" ] } ], "source": [ "rad = Radiation(\"20180101\",\"20180103\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Load some variables:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Scattering needs to be one of: direct ,diffuse ,global\n", "Longwaveradiation at the surface is only measured as diffuse radiation. Setting scattering to diffuse!\n" ] } ], "source": [ "lw = rad.getRadiation(\"LW\")\n", "sw = rad.getRadiation(\"SW\",\"direct\")\n", "time = rad.getTime()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Plotting the just loaded data:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtwAAAEVCAYAAADTpjKIAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzs3Xd8VFXawPHfSSP00FuAUEQE6SggxYZKca3ouq4dF111VxddxY4dd197R7FXdFVUBEVE6UjvHQKEGtJJL+f9485MZpJpd2buzCR5vp8PzMydO/eeJLc8c85zzlFaa4QQQgghhBDWiIl0AYQQQgghhKjNJOAWQgghhBDCQhJwCyGEEEIIYSEJuIUQQgghhLCQBNxCCCGEEEJYSAJuIYQQQgghLCQBtxBCCCGEEBaSgFsIIYQQQggLScAthBBCCCGEheIiXQArtGzZUqekpES6GEIIIYQQohZbvXr1ca11K1/r1cqAOyUlhVWrVkW6GEIIIYQQohZTSu3zZ72IpJQopf6llNqslNqklPpMKZWolOqilFqhlNqplPpCKZVgW7ee7fUu2/spkSizEEIIIYQQgQh7wK2U6gD8ExistT4ViAWuAp4FXtBanwRkARNtH5kIZGmtuwMv2NYTQgghhBCiRohUp8k4oL5SKg5oABwGzgG+sr3/AXCJ7fnFttfY3j9XKaXCWFYhhBBCCCECFvYcbq31QaXU/wH7gULgZ2A1kK21LrOtlgZ0sD3vABywfbZMKZUDtACOO29XKTUJmATQqVMnq38MIYQQQgjhp9LSUtLS0igqKop0UQKSmJhIcnIy8fHxAX0+7AG3UqoZRq11FyAb+BIY62ZVbf+Il/cqF2g9HZgOMHjw4GrvCyGEEEKIyEhLS6Nx48akpKRQ0xIVtNZkZGSQlpZGly5dAtpGJFJKRgN7tdbpWutS4GvgDCDJlmICkAwcsj1PAzoC2N5vCmSGt8hCCCGEECJQRUVFtGjRosYF2wBKKVq0aBFU7XwkAu79wFClVANbLva5wBZgATDBts71wCzb8+9sr7G9/6vWWmqwhRBCCFG7aG38q6VqYrBtF2zZwx5wa61XYHR+XANstJVhOnAfMFkptQsjR3uG7SMzgBa25ZOBKeEusxBCCCGE5aafCY8lRboUwgIRmfhGa/0o8GiVxXuA092sWwRcEY5yCSGEEEJEzOH1kS5BrdaoUSNOnDgRkX1HalhAIYQQQggh6oRaObW7EEIIIYSITo99v5kth3JDus1e7Zvw6J96+7XubbfdxpgxY7jooou49NJLadasGe+++y4zZsxg7969PPnkkyEtG0gNtxBCCCGEqENGjRrFokWLADh48CBbtmwBYPHixYwcOdKSfUoNtxBCCCGECBt/a6KtMnLkSF588UW2bNlCr169yMrK4vDhwyxbtoyXX37Zkn1KwC2EEEIIIeqMDh06kJWVxdy5cxk1ahSZmZnMnDmTRo0a0bhxY0v2KQG3EEIIIYSoU4YNG8aLL77Ir7/+SkZGBhMmTGDChAm+PxggyeEWQgghhBB1ysiRIykrK6N79+4MHDiQzMxMy/K3QWq4hRBCCCFEHeA8BvfEiROZOHEiAPHx8eTn51u6b6nhFkIIIYQQwkIScAshhBBCCGEhCbiFEEIIIYSwkATcQgghhBBCWEgCbiGEEEIIISwkAbcQQgghhBAWkoBbCCGEEEIIC0nALYQQQgghhIXCPvGNUupk4AunRV2BR4APbctTgFTgSq11llJKAS8B44AC4Aat9ZpwllkIIYQQIiAZu6FhS0hsGumSRI85U+DIxtBus20fGDvN6yr5+flceeWVpKWlUV5ezrXXXsuKFSv4+uuvmTVrFldddRU5OTlUVFTQq1cv9uzZE7LihT3g1lpvB/oDKKVigYPAN8AUYL7WeppSaort9X3AWOAk278hwBu2RyGEEEKI6PbKQGhxEvxjVaRLUufNnTuX9u3bM3v2bABycnJ48803AVi0aBGnnnoqK1eupKysjCFDQhtqRnpq93OB3VrrfUqpi4GzbMs/AH7DCLgvBj7UWmtguVIqSSnVTmt9OBIFFkIIIYQwJWNnpEsQXXzURFulT58+3HPPPdx3331ceOGFjBw5ku7du7N161b++OMPJk+ezMKFCykvL2fkyJEh3Xekc7ivAj6zPW9jD6Jtj61tyzsAB5w+k2Zb5kIpNUkptUoptSo9Pd3CIgshhBBCiJqmR48erF69mj59+nD//ffz+OOPM3LkSObMmUN8fDyjR49m8eLFLF68mFGjRoV03xELuJVSCcBFwJe+VnWzTFdboPV0rfVgrfXgVq1ahaKIQgghhBCiljh06BANGjTgmmuu4Z577mHNmjWMGjWKF198kWHDhtGqVSsyMjLYtm0bvXv3Dum+I5lSMhZYo7U+ant91J4qopRqBxyzLU8DOjp9Lhk4FMZyCiGEEEKEz9SmMDUn0qWodTZu3Mi///1vYmJiiI+P54033qB3794cPXrUUaPdt29fWrdujTFmR+hEMuD+C5XpJADfAdcD02yPs5yW36GU+hyjs2SO5G8LIYQQQggzLrjgAi644IJqy4uLix3Pp0+fbsm+IxJwK6UaAOcBtzgtngbMVEpNBPYDV9iW/4gxJOAujGEBbwxjUYUQQgghhAhKRAJurXUB0KLKsgyMUUuqrquB28NUNCGEEEIIIUIq0qOUCCGEEELUTgWZkS5BVDHqUGumYMsuAbcQQgghhBVKCyJdgqiRmJhIRkZGjQy6tdZkZGSQmJgY8DYiPfGNEEKI2m7XL/Dx5XD7SoiJhRbdIl0iIcKjBgaXVklOTiYtLY2aOldKYmIiycnJAX9eAm4hglWSD8tehxH/gtg4KC8DpYzAQghhBNsAr51mPD6abZwjQog6Iz4+ni5dukS6GBEjKSVCBGvuFFjwJGz6ynj9RAt444zIlkmIaNJpmOvr98ZC5h7I3BuZ8ggRNlLDLQwScAsRrCMbjcdD6yqXpW+LTFmEiEb5x11f718GLw+Al/vD4fWRKZMQQoSRBNxCBOvQWuNxxRuw9fvK5U+0gpf6wcavIlMuIaJFxk7P7701KnzlEEKICJGAW4hQ+uKayuflJZCVCv+bGLHiCBEV2vWLdAmEECKiJOAWIhyKciNdAiEiJ7ae5/dOOj985RAi3GSUEmEjo5QIEQ7TOsKFL8LGL6FRG9j8NTycYYxqIkRt520s4oSG4SuHEEJEiNRwCxEuP9wF+5YYwTYYHceEqAtKCz2/t/kb97WAJflQUWFdmYQIC6nhFgYJuIUIVts+gX3ugwthalMJKkTtpjVk7va+TkaV98vL4On28N0d8NODRvAthAitigp45zzY9mOkS1InSMAtRLDKioP7/OF1vtcRoiZa+Q68ObLy9flPuV+vzKkGPO8ofHSJ8XzdJ7DsVVj6inVlFKKuKi2AtD8qO/Z/PMGoBNo1P7LlqqUk4BYiWKVFwX0+dVFoyiFEtJl9NxzdWPn6jDvcr/fmiMrnz/Wofk789gxsnxP68glR2xRmmwiY7ekutllfd80zHj++DAoyQ12yOk8CbiGCUVoEOfuD28a8R+DZFJjlIRgRoi5Y9pr39z+7KjzlECKUwj1KyczrjIA5P8P3uvayKVX9vbfPDm25hATcQgTF3gEyWIVZsPYjeGe0DCMl6qafHvC9jgyvKYR3x3cYj2VOLa+rP4Bvb3NTa12lhttZYlMrSlenRSTgVkolKaW+UkptU0ptVUoNU0o1V0rNU0rttD02s62rlFIvK6V2KaU2KKUGRqLMoo5Y9Dx8f1f1oPfb22DHz9XX1wF0eKzXxPN7aSth/3Lz2xQiEkqLgu/D4Gz9597fP3EsdPsSIizCXYHiJnj+/p9Gf4ivbnJdbr/PleRV/8zh9UY+t3zJDZlI1XC/BMzVWvcE+gFbgSnAfK31ScB822uAscBJtn+TgDfCX1xRZ8x/DFa/Byveqly29BXjYvXpFbDmI9f1cw+HvgzvjYGlr8roJSL6FGbDlu8qXz/VBp5sbTz/9nbjBh2Mb27x/v6rg4LbvhB1WU5alQVOXwbKy9x/5rmTLStOXRP2gFsp1QQYBcwA0FqXaK2zgYuBD2yrfQDYuqlzMfChNiwHkpRS7cJcbFHXbJlV+fznhyqff1clz7phC2v2//ODRvAvRDT5ZALMvBbWfQZlJZXL59wH6z42nu9fARXl1pWhKMe6bQtRm2XsdH3t3JK78L/uP+Nt0iphSiRquLsC6cB7Sqm1Sql3lFINgTZa68MAtkdbtQkdgANOn0+zLXOhlJqklFqllFqVnp5u7U8gar/9S+H/TnafRuIsJoCZIv1NQ1nyImz9wfz2hbBK2krj8dtbYdbtlctXvFn5/N3zYe79kL7DmjJs+Q4Ob5C+DkJ45cf54ZwOlr7NuqIIIDIBdxwwEHhDaz0AyKcyfcQdNwlJ1Y8krfV0rfVgrfXgVq1ahaakom5ZOcP19YkjRhqJNxtmWlcegC/+Cp9eFfzQg0KEwoBrK59v9HLs//EWvHaaNWVIXQxvjYQXA5xwSojazN2II54839O/9VIXw8sDrUmhrEMiEXCnAWla6xW2119hBOBH7akitsdjTut3dPp8MnAoTGUVdcnsyeY/E9AY2iYuiAA75hi5slKjJyKtwkOeZzj85QvjcYOtY2WOU8PngT+MsYetTGURIhDRet2u2kco96Dndd8fb8wWu+ELa8tUy4U94NZaHwEOKKXsmfjnAluA74DrbcuuB+xJtN8B19lGKxkK5NhTT4QImdJC3+vY2W/qx3d6X88TMzUQzn5+yOiUlrknsM8LEYyKclj/WeT2Hxvv+b0Z5xljD783NnzlESKa+Qr0Fzzp+tqeLuZNeYnvdYRHkRql5B/AJ0qpDUB/4GlgGnCeUmoncJ7tNcCPwB5gF/A2cFv4iytqPTPDjT3d3ng8siHAnQUYcC971XjctzTA/QoRBOeUq9FOHXq9DXMZSu76S2yfA5u/qXx9YEX1dYSIpEArWAKVd8T2xEfAfXC1+W0veEpmoAxCAD2+gqe1XgcMdvPWuW7W1cDtbtYVInRe6uv/umVFsDeI6diDvf6WST63iIA5/658PuIuaNfPmKzp8hnGjdjTKAehEhNbfZnMPimiXbhTSrStBTZrHyR18rJegOX6Txe45mvoXi1cEz7ITJNCBOKDC10vWA1ahm/fs+8O376EqOq+fcZjt7NhwrtGDd4pF1m/XyW3KyG8OuzU6prU0fN6gN8T8px0QfVl39/ld5FEJbmCCREK9ZMiXQIhrGMfc7teE/fHuvajs2LQU0WHuWleiJpm5dv+r1vi5/jajdtWX5azH47JMIJmScAtRKAjGxQ7TYfbvJuJD0rgIGqY3542Hos9TPPcupcfG7Ed92f8M8BCROloD0JEC+dxtX2JTfBvPU9fpguz/N+XACTgFiLwntc/ODWruasF8CTcnWiECNbiF7y/H1fP9zbsx/35T8Alb3pfVwhhnsuwfT7uM+76RLhTdfhAuxNH3C8XHvkVcCulmvvxT9rURc0U9qGOQhBwR+vYrqL2ydhd+fy2GjAKiHPLkxB1VojuER0Gul/+5Q2h2X4d4u8oJYds/7xFCrGAly6xQkSZshKjI1YoZs8yU2sdihruQ2s9XwiFCCXnL6QtTwpiQ2Fq2XkmGabmhGdfQtR2DWXm7lDxN+DeapuG3SOl1NoQlEeI8CjMgmdTQrjBMKeJvH023LEaWnYP735F3VN8wnhUMf43Q7sTyVSqHT9B9n44/W+RK4Oom7L3R2a/vlpBJbUx7PwNuIeFaB0hokOoB+83NWRZiC50BRmABNzCYlu+NR5vmB3khpyPe4tToioqIMZ2Tk51Gh2lXX/oeJq1+xbCmfaQA239jq3fxZKXYXignaDrHp9RglLqPOAVpVR/2+tJ7tbTWstsHKLmKPVzSCR/RaK2YMssmHE+bPwq/PsWdYd9htOKssiWw0y/BU/n44zR5kZyECJYTZ3Gw34ijOkZwfTzufQt/9ab97CcTyb4Uy13G/Bv4Bql1DkYU7ELUbO9OSJy+w5VcL78NWMq6/9NDM32hPCmZY/gPh8tTdhfXBPpEoi6Kuwd9D3IOQh7Fwb++aTOlc+P7wi+PHWEPwF3utY6W2t9D3A+IO1xQgTFgsAjdTHkHw/9doWwMzP0pVtBHveBBOzuhjTb+XNw5RCipns1yDCuzCmhIZKVVzWMPwG3I3FPaz0F+NC64ghRU5kIBrqPDv3u3x8P/zUz+Y4QYRZsDXcgTeRLX3a/XNKwRG3n7XwpzQ/8swAnjpovj/AdcGutZwEopVraXr9idaGEsNTqD0K/TTPBxPA7Q79/IawQ0vHeI5BS8suj7pdLGpYIm0jNmSBzNUQbM0MrvGtZKYQIp++t6FVtZhxumeBV1BD2zsWjp4Z2u1ZP3CQTQ4m6LlTngD+VSSU+aswFYC7gjpIeL0IEYcdPkS6BEDWHvV9A/eaRLYcVyqKkA5uomV4eALPuMIad/OnBSJcmxEwG60tesqYYtYyZgFuqDETNFw0jFETLaA1C+PLTA8ZjWQhGfY22437dJ5EugajJMvfA2o+M5/ahM6OKhSHbyeNdXx/dbN2+ahGp4RZ1S7QMyyRETbDtB+Mx91AINhZlt5DD6yNdAlEXRCq9KVT77Ti0+rLmXVxf268TwiszAff9odqpUipVKbVRKbVOKbXKtqy5UmqeUmqn7bGZbblSSr2slNqllNqglBoYqnKIOibUs0s6i5bau/TtkS6BqI3OeTj4bQR9joQ4cFn9Xmi3J0RUCeZ8cTpXG7cJuiTC4HfArbXeFOJ9n6217q+1Hmx7PQWYr7U+CZhvew0wFjjJ9m8S8EaIyyHqirrQ7PXbtEiXQNRGsXGRLoE1ds2PdAlErVcTs3FrYpmjn6nhEpRSg5VS3yil1thqmzcqpTaEqCwXA/bx2j4ALnFa/qE2LAeSlFLtQrRPUZdoN5NgmHHe46Eph5U2fx3pEgjhgXMNt9U3dD+3v2OutcUQIlLCncribpIp4cLs+GSfAO8BlwN/Ai60PZqlgZ+VUquVUpNsy9porQ8D2B5b25Z3AA44fTbNtsyFUmqSUmqVUmpVenp6AEUStV+0fGu3OP0kdbG12xciENGSduXs0NpIl0DUdhEbojLM+5Wx7X0yG3Cna62/01rv1Vrvs/8LYL/DtdYDMdJFbldKjfKyrrurdLUjSWs9XWs9WGs9uFWrVgEUSdR6Zmq4mySb3LiZYMLiC+H7432vI0TYhTHgrijzb72WJ1tbDlE7lYZg1B6rhTvQl9ZVn8wG3I8qpd5RSv1FKXWZ/Z/ZnWqtD9kejwHfAKcDR+2pIrbHY7bV04COTh9PBkLRZV7UNWYuQN3ONvf5aKy9EyIYWYHUpXgRzlOkvNS/9UpOQM5B+OZWKCu2tkyi9njKTEfCGljDLRNHWcJswH0j0B8Yg5FKYk8r8ZtSqqFSqrH9OXA+sAn4Drjettr1wCzb8++A62yjlQwFcuypJ0KYYuoiYuEFRy5moiaIttQkK86bknyYex+s/0zyuYUIxmk3R7oEUc9s1/N+Wus+Qe6zDfCNMmoE44BPtdZzlVIrgZlKqYnAfuAK2/o/AuOAXUABRtAvhBDCSvbhwM57IrLlsFLJicqx+WXmSeGPvQvNrV/Tx+H2V2MZy8IXswH3cqVUL631lkB3qLXeA/RzszwDONfNcg3cHuj+hKgU7AXI2+ejLKXk6BZo0yvSpRA12e4FxmN+qDqhO50jgQQDZtK2YmL9W6/kBBzZaDz/+mboe4X39YUwPZSkn8d67iH45THTxQl6v6Gy/UcYdU9491nDmE0pGQGsU0ptt2BYQCGstTzIIdxDVmMQhgthdojzb0XdY5+uuu+VkS2HnZnzL8bPuqSS/MDKIuowk9fvqsftfA8tRnOnwIbPAyuSP/s192Hjob2JeQYPrg5if3WD2RruMZaUQgirrf0YdpuomXB7rZJOk6IOSmgYmu1E4zkiAbcwa8lLwX1+0f/BuSGYuTUYcfWhrND3eq1Psb4sdYhfAbdS6i5gCbBWa+3neEtCRJFZIchKClUNdzhy66RjpgiVhEYh2pAE3EKEj5d7QEJD/wJusyrK/U/nqoP8TSlJBl4CjimlflNKPa2UGq+Uam5h2YQIjf+eZP4zURgbmJKfLjN/icBoDdM6Vb6Ob+Df5065yPv74azhzvNzIKuSE9aWQ4hIDc/n7bNWBcWLnrdmu7WEXwG31voerfUZQFvgASATuAnYpJQKuAOlEGGRf6z6srs2wWXveP6M2ZSSaPP9P+HxZv6PRyyE3ab/QVFO5Wt/A+6Rk60pTwCyP/tbYB/MSg1pOYSIWGtjTlrgnw20zMc2B77POsBsp8n6QBOgqe3fIWBFqAslvJu58gDrD2RHuhg1W1JHaNPbywpuLjg1KN52mHt/pEsgQmjupiP89Z3l1u6kQQvX1zFmbxOeONdwW3sybT+czca0HN8rVlPTm7Zql/IKzdJdxyNdjCBFqIZ75rWBf9bBw/nQtKP75VLB45VfV1Kl1HSl1BLgC2AYsBS4wjaVuoyLHUblFZp7/7eBi19bYvm+cotKOZhtQZ5XTWBprUQYI/eVb4dvX8Jyt368miW7MqiosPAY+tWicbfDmFKi0Hy5+oD5D0r+aVSZsXgPV7+zgs//2B/pohjyjpr/TDT2pyn0t8LOTdn/+j9o0t796hXSxc8bf6suOgH1gCPAQYzp1qWKNQIO57gGwNrHyZxXVMo3a9PQWpNd4Dqxw1u/7+bdxXs9frb/Yz8zfNqvHt/PKShlzf4st2UoKi2n32M/syGtDh4mvmbcCvcFuJbWOizdfZxVqZmRLkZElFVoMvNLKCwpD/3GAx7ey+qA2tx58+GyAIbG9DsQCa2Zqw6QMmU2T82WDE1nr/y6C4ApX2+McElsnusR6RKERly9wD9brxEez8Vaeq8JFb9GKdFaj1HG1JC9gTOAu4FTlVKZwDKt9aMWllE4yS6oPKCvnbGCRTuPM+P6wZx7Shu36z/63Wa+XnOQf32xHoC5d42kZ9smADwzZxsA1wztTEKc63evJbuOY69Em7vpCGNObQvAop3pXDvjD7f7Sp02Hq010+Zu463f9wBww3srWfPweWScKEYDLRsFcaIHoiiQZmVP/Lzhj38OVnrJDw+3E0ehaXKkSxFyV79tZLOlThvv1/pl5RU89eNWbjgjhc4tQjTUXYQUlJQx8Il5gP8/f+T5CMhbngzHt4enKJ6EsRY+K7+EmasOoBR8vNyowX170V4eHB/5Cau01nS5/0cA3r/xNM46uXVEypFXZNSYtmiYEJH9m3E4p5B2TetXf8PfMeHdMu45ZeUV7Eo/4bh3By0uEYpzQ7MtZ8V5od9mLeJ3cp42bMKYan0OxjCB3YA7LSqbcOOr1ZUdIRbtNHLbJn6wiowTxdXW1Vrz9ZqDLsv2pBvDYDk3Sd/71fpqn/3rO5Wp+bd+vJqUKbP5dMV+j8E2wPETxSzfk+kItgEy80tYuCOdQU/+wuAnf2HOxsOs3udfrWRmfgmPfb+ZotIgavHKA2ziclcDEI1Ng/54wX2uenZBiaP1I1psPpTD5JnrKC33PsJKmdP73lppnG07ksd7S1K58/N1AOw8msddn6+t1vITSkWl5eQXmzsGc4tKWbY7w+179WxfjLcdqYE3Nl/BrM9g1/pguKwoj2N5RZbvB+CG91fyzJxtPP3jNpcfPa+oei3hsdwilu62Ppf5gW828uqvO9mfWVBZzvdWWrrPotJyjy01ifExtkcj1SczvyQkraaHcwpDft0b9syv/LHXzb0tiIA7t6iUlCmz6f7gHMa8uIiUKbNJmTKbjWk5XPHmUno8NMf09cUcd+ecl/Mw7Q/Y+r1lpfFk2e4Myq1MswsRf8fh/idGzfZwoBQj2F4GvAtESVtP7ZddUML7S1PdvjfoyV8cz5fdfw7tmtYn381F7LZP1vDRxNNJz6sM0BdX6ZTiqZn+gW8q/9RvXjOQ83u1pbC0nBmL9/L8vB0MfvIXtzXY171bGaT//ZM1ALz+14GM69MOMPLSS8oqqJ9gXFRzCko5UVLGO4v28N6SVNo1TWTSqG5uy+ST15u4lxO0WWfX170uCWz/UWz6wj28/ttuOiQ14PQu1o/wuWZ/FolxsfRq77mW5vp3/+D4iRJuHtHV63qHcyqDosd/2MKNw1NQPgI2e+uQ/Qb10vyd/LDhMBf2bU+F1rRolMCgzsbvYVVqJj3bNaG8QvPNmjSuHZbCugPZnNqhCfXivOf53vn5WhLjYnl2Ql/GvbyIY7nFbHrsAq+fcTbx/ZWsTM1i02MX0Kie6yW6bdNE9mUUcNV0iztORiX/b6gqwH4Sn336Hi/njGL6beMY0KlZQNvwl3PH930ZlQFuVn4pjRPjXdZ9eNYmftp8lO1PjvF5/AWqtLyCT1cYNe0DO1f+7G2aWNsqef4LC9mfWVCtpaa0vIKiUuOLdfqJYrTW3Pj+StYfyGbNw+fRPMBa7+Mnihn2zK/cemY3poztGXT5nW1Iy+b0Ls0557nfHJVbJ6v9/FTlV5gyZXa1z37Y6Dijqiwb8PjPQPW/959eXex43vvRnwB474bTODuo0vvLx7n1xTUwNZQty54dzS3iijeXsT+zgAfHncLfRnUNy34D5W8NdwrwFXC61rqr1vparfXrWuv1WmsZ7NciWmuW7j7OrHUHmbF4L/0fn1dtnQ5J1Zuwhj3zKylTZvOPT9e43e61M/5g8szKWu3jJ1xr+Ca8ucxruVY9NJoxp7YjJkbRsF4ct5/d3Wlb1WvaAa4Z2okvJg11vL7tkzXkFpWiteberzZwyiNzST2ez7LdGZz/4u8Mn/Yr8bHG4Tlr3aFq28srKuWc535j/tYAOrEE4GBpqCb/gGgZ7mSDbRSHnMLA8u5mrjzAAaeasKqKy8o54hQYX/b6Usa9vMjrNu3H4mc+OknlVqkFLPSjFeTXbcbwkPafd8tho0n1l61HmfTRai5/wzjucwpKmfDmMsa9tIgHv9nI1O+38N6SvVz+xlJeW7Db4/aP5RWhtWbWukN8scrosLcnPZ8TJmugVqZmAfCRmxzv1hW5AAAgAElEQVTkNk0Sqy0LqgUolHx84dHONWPuahejoKXl2uLPWZl4G79tT7d0P+5qJUfb0gKz3LS4/LTZuM45pxSGWmZ+5X7t5+2Y3m05llccstrDigrNzJUHXM4J59r01xbsou/Unygtr+CELZ0kpUUDSsoqyMgvcXxJ2XYk8HQI+5ebT1aYzPEvzPK5yoniMlKmzHYE22bkFXv/26ZOG8/Xt53h8f3HPvjOxN58/T2D/HuHIbXknUV7OPe53x3Hz1M/brV8n8HyN4c7egZYjVJaa581bP4qLivnlfm7eHXBLrfvXzawA89d0Q/AsU+tNeUVmrEvLWLnMWMyhwW2m8b7N57ms1mwsKScj5fv453Flekg/5nQlysGJZOZX8L6tGzO6ek+Tzw2RrHigXMZ8vR8x+u3rxvE0K4tiFGKLYdz6Z+cREyMYv2j59PvsZ8B6Dv1Z5ftnPV/v7m8nr7QKMvJbRpX2+ee9Hz2pOcz8YNVTL92EOf3buv+BwvRTfyXrUe5slV73GTomReBwCJlymzeu+E0UHBWj1YopRzxUV5RKb/vSCe7oIR/fbGOeZPPpFur6l8wVqZmcvfM9Xx3x3BKyiu4938b6N2+CU0S45l0ZlfOrpLn+cK8nbz5+262PTHG0STsLK+olOyCUpo1TOD79Ye4YlBlnvlHy/fRJ7kpVw52P/yUPbfT+XWDBO+Xs3eXGKknGuPGaL8pfr7SdTSL7EIj8NifWeC4mD8527iYL9/jPtVj4Y50rnv3D96+brDb98srNLEx7q8P87ceZcXeTB4Y5zqN8rNzt/H3s1xbdrLyqwdjB7ML3f693Nl1LI/Scs0p7UKUC2pCXnEZ3va689gJTgrVCIRBemn+Tsb1acfJbV2vPUWl5W6PZbPswe2I7i0dLYz9kpvyy9ajHistwAjG3X3pCsTOo3lM/X4zr/xlIM0bJrgE3PYKmW6tG6I3Q26hcZ4Ga/X+LO793wZ+2XqU6VXOlcKScv77k5HDf9pTvzi+XPRJTiI1o4BnbX2OwOi/cfd5PfjHueYnNbP/fk13OF7/hc9VTjhdl/p3TGLdgWyfiVA3nJHiseW6qoGdmlVrCbDXlt8S+0O19e3vpQZ6yHgqvK972Htj4dbF3tcJQnZBieOaXJP4Oyyg+6pSk+vUZo/M2kzKlNnc+fnaarlhK/Zk8PHyffy0+QhzNhozoBWVlnP8RDHFZeWOHNLjJ4p55set9Hn0Z5dg+6yTWzH1T72YMrYne54ex/NX9rcFTJVng1KKuNgY5k0+k9Rp45l1+3DHewM7N2PHk2O5a3T1i1OHpPp0SKrPFW8t5akft3I0t/Ji3yAhFqUULRrV8xhs27VpkkhvWwrAgI5JnNOzDQ0S4kiMj2Vgp2bE2IKNpvXjeWCcuWa879YbNdzlFZofNhwiu6CEu7+srKGf9NFqrn57OYt2prvcNIwP32FqX5VcrzQaxfSF7r8AGe9Ht4FqBzd/uIob31vJnE1HABwBYFZBKde/+wd3fr6OCm2MXuOOvenulo9WczzP+D1vPpTLsj0ZPPfzdr5ek0bKlNlM/c6Y/OCthcZ2rn57udsm1Cn/28jI/yzguhkruP/rjfy0+ahLJeme9Hx+3nzE8bqsvIL9ttqp3emuMwTa817nbDzMvgz3tUutGhvturFKOc7DqkrLK/jUW+260x+6rLyC4jLjpv3sXCMY+G175SRL//xsreO5t7zgiR+sYvrCPT7z1sE1lcbOTK3n6OcXMvYl760MVimrQW2hz8W/wSOzNrksm7nqAP0e+7naSFGBsLeypLSsnFTInkLlbSjWnABruPOLy/jbh6vYfKiyqf+TFftZsiuDX7Z6rj1vZUsRPJJbedz9vPkIl76+JKChKe3n6c9bjlJRoV3ulRn5lfce57LY7ytfrnadyOW5eTtM7x8gw9aKVma2/H4Msbpqn1EL3je5Kd/ePpzUaeP54R/Dq633j3O6kzptPKnTxjP1ot7VUsfs/EmNsgfgRURRx9Ij1mYaO6fEAnRrZXSCL4nyi4y/2fynKKU2eHlfYUyEU2fZ/+Cz1h1i/YFsnrmsL4dzCjm9S3P+XCXf8vphnfmgSnNx5xYNHE1dbZskcnLbxvxnQt+AazP6dUyq9k34rtE9uHpIJ8a+uIinL+vDBb3bMvW7zby/NNXtRb5LS3MjOTxxyalc9vpSLh/kfUSMcX3a8fSP21yW/XTXKJo3TOC6d/8gMT6Gtfsr8xvLKjQ9H55D48T4aiea3dLdGSzdnUH/jkl86/Rlg+0/eimJuRYJb5e+0nIdTZe7ar6uN5WUok8BI1XhtJTmjmbiQ1X+9s0aGD/J8j0ZNK0fzx2frmHiiMrcuBV7M3l+nutoEvnF5Y5asfeXpvLAuFMclSBr9rt2ciosKefDZanMtgW99vd3HM2jQXyso+/Bm7bA/4FxPZk0qhv//Xk7b/2+h55tG9Opuevsh8fyiunashF//2QNTevHs/7R8x3vXfzaEo7lFjmOnezCEoptF+a3rxvM3z5c5Vj3pAfneP9FOrlmxgo2H8plw6Pns/mQ0cT9yYrKYN3+RRHgYJb7EQycJ2d59ddd1b4UX/r6Ei4fmMzR3CL+Nqqr2/SUtKwCBnX2nW88a11lB+ovVu7ngt5tSWoQvqO2IsiWHa112KaluTx2EXdX6QD34bJUissq2HQwt9rfcsG2Y6zZn8Xd55/s1/btKVFdWla2TPSwteQ9MmszV5/eibjY6vVh2QGmf731+27mbTlKXIzijWsGAZVfAp/5cSuFJeVuc7Xto/k417pPnrmeE8Vl7M8sIMXEPSIzv4Rdxyq/KP+y9ShndG/peD3i2QVuPzf21LZMm7PN7Xt2mw/loFBe+32AcQy9s6iyFbesvMLt79mt1qdAhudKF4B1tpSXa4ZU9gGKc9Oydd2wFJfXS+47h0XPPOdfOdxIfeJceP85Y9BmH1KmzGZjwzKqtxubEdkqpnSn4/Gcnq05p2drHvp2E9kFJbQOUQuQFfwNuP2pkoySRMLIuGF4F64e0pnRz/9OakYBf3nbc6emqsE2VOaV/d8V/bh8YIeQpadU1bpxIqsfPs/x+tIBHdw2Z618cLSjRtBfAzs1Y+WDo2nZyPtNPLlZAzY9dgFzNh7mkxX7Obdna0fT7Zw7RwJGC8D2I3nEKMWfXl1MUWkFRaWVJ1mMgrNPbs38ba7Ttq9z6ohUUaF9NOH4f9HQKPZXeB4aq6SsPKoDbmfL9mRw2lOVnWz3VKktzioo4ePl+3jo28oaPucOswC/bHX9vZdXaLq2auhI0/CWW3/KI3PdLn9p/k6UMprW1zsFok//uI2bhnfhly3GNrcdyWPbkTySGsTTNzmJhTvSWbjjuCM4yCks5c3fd9MvOYknZ29xBMN2RaUVjgBiRPeWtGpcz+MXuariYivPy+V7jIDs5g9WeVrdYc/xfAanuHZMPZxT6NL56bUFu7hpeBeXddbuz3Z8+Zy3pfrvNCE2hi2Hcrm4f4dq71VNc7OPzgJw3/82Gl+cruwPGOfKB8tS+VPvFrSsuiE/vbN4L95GoM8uLPO6bV+1eQWl5YR7MEd7q8PHy/ex6aBxHO1OP8F5tGHr4VzGvrSILi0bsve4cdz//axuHMou5KNl+3j4wl4egzl7DW5Pp5SVZg0TaNEwgYz8Eg7nFNHR9qXSOUc/0BF17F9iWztd02NUZQvXo99t5qlLT3X5zLe3D3f0EdqQlsMZ3VpyJLeIZg3jOVFcxj4vAfecjYfp0baxS6rT0KfnU+LUijPpo9W8f+Np1T778cQhDO3anNgYRU5hqccvhSktjN/P3uP5jH/ZOI/WP3o+TevHcyi7kGfmbOPJi09l86Ec2ifVJ6VlQ05+eK5LLeixvGLau+kH5ZaJ0Td8dUKvem9t2iCeC/u0gyrDsH8y8XT+PMOP5IEf73Y7fn7qtPFMX7gbqkynUVxWTmOv4wlEb5vt9IW7XSrszu7Zmqb1jU7GOYWlNT/g1loHMIOAZ0qpWGAVcFBrfaFSqgvwOdAcWANcq7UuUUrVAz4EBgEZwJ+11qmhLEsoJcTFsPDes1my67hjWL1ze7bm0oEdGHtqO2JjFMdyiziUU0TPto1DkgsYrH4dk9jx5FjSsgro0rIh87YcJblZA9PBtp2/n2tUL44rBnfkCg85uonxsfTrmITWmkcu7MXz83Zw04guKOCu0Sc5AomqgSG49gD3lLt25n8X8Mo5CfT1UL5n5mzjfteBAqjwEr5nF5bif7fK6LqYLajSQWzmqjRmrkrzsLZ75RXaZSx3+2g0ZmkNlw9Kdgm4Af7703YSqozO0KZxIjOuH8xJD87hzd93O2rEAY81Yg9f2IsnftjCgcxC4mMVifExLJtyjq0VxfWLQO/2Tdh8KJfGiXFc0r8DHy3fx9LdGaRMmc3EEZWBcdUvfe7sPFrZiWjXsRO0aJjAoWzX9JCyCs1eD+kwUDkU4Cc3D+HOz9dx21ndePyHLby1cA9vLdzDygdH06xBPLlFlWN033BGChMGJbvN8XYeGWPN/iwe+34LK1dk8rrPn8a9r9cc5GYvp39puSZlymzWP3I+21IzGWJy+/szCjjF92pA6AYQXLEnky9WHeB7p9aKt37fzS2juvKqbVIWe7ANRmvRI7M2s3R3BuP6tGNI1xZut2v/guecI96oXhzPXdmPG95byUfL93H/2J4opVxaoKp2cPeXvZNmakYB2QUlJDVI4Giu6/GXadv2ZQM6ML5vO/p3THL5mQtLyl1SHQ97SH3RWjvOf+dW1hI3KVPu+hYN6drc8UXFHmzvfWacY1zwD246nbmbDvPZHwfILy7jbKe+PwezCmlaP563ft/N9+sPOf5uHZLqG31PqqQcHMkt8j/gNsHli4jfwWv19YZ0cX/8VHOo+tC+djd1qd7Zs6Xyt9OpyWEBw8DetwuMPPmrTuvIClvlx4GsAk5y0+crWgQzInsw7gS2gqMPzbPAC1rrz5VSbwITgTdsj1la6+5Kqats6/05EgU2Y3j3lh4no2jdJDHqvoElxMXQ1XZD9tj5MEKUUtw0ogs3jeji9v2rTutIfKzion4dqtWcnq48d6rYl1HAv7/aVm24Jruql74Lerdhw0bPF84fNx5hkp9n01/fWc4n/q0aUvGUURrkKT+gU5JLuk+3Vg3ZnZ7vNe8U4Lkr+tGhWX3HcHYPjT+FvslJPD9vO22bJHLtsBQuf2MpYDRjXzEomS9Xp3H2ya1YsD2dt5wusnbtkxIdI9n4q0OSce4t2H6M1o0TbX0fFFVHWptz50jKKzQXvrKY9k3r89hFvfloeWW9www/x/62e3vRXt5eVPmZGAXPXNbH8bpT8wbszyxwae72ZHj3lqx6aDRgDIlo59xqYff+0lQ+X7mfgW6GuLPnnwPcb5vJ747s/5qYncFQVFpOz4fn0svnfdhYod/jP3NVbBpDqnyh7dKqoVGt4sHjP2zhMz+bkQIdFtBZ84YJXDOjcj6CB8b1ZMmuDH7fkc65z//udiSK0c8vdDz/Y28mAzo1IyEuhuyCEh6ZtZknLjmVPekn+HTFfhrXi6N5gwT+cU53R232Gd2MNoDpC/eQGB/L5PN6uKRhVE3/qkprzd8/XsPglGbcPLIyDcye+//7jnQuenUJz13ZzzEajt1z83bQtH48z/+5v8vyKwcnM3NVWrVO/AezCymv0Lz4yw4Gdm7GWT1acSyv2KUzYsaJYh76dlO1DpfPX9nPZaSsHm0asePoCTok1Xd7TiulaJhgpJu1bZLoyHm2D4lnN+7lRXx7+3AOVenrcDC7kAterOy7MOP6wUz8YBUHswrdnhuBaJAQS4EVM7/6w9O5959uxBWEdvz2x77fTIPiY/w7pFv1T0WF5viJEgZ1bsZXtw5zVLy1s13Xb3p/VVRPBBb2gFsplQyMB54CJttmsDwHuNq2ygfAVIyA+2LbczCGJXxVKaV0NM3UISIqLjaGP5/WyWVZ6rTx/L4jnTM/vdrDp8xrn9SAa4d2MtpffEiZMttrr/Cs/BII84SbAAmU0q1tM5dJU7Y9MYZxLy9iT3o+398xgtX7Mpn6vRHEDenSnKFdW/DS/J0ATBiUzP1jezrGfD+pdSPm3DmS2RsPO9IVrhnayTFr3it/GcC9X22gvEJzwaltaVQvjjvPPYn84jJHMPD5pGGAcSG9dEAHDmYXMrxbC/p0aErH5g24eWQXvl17yJHS0rlFA/42siuz1h3kxT8PAGDdI+dx8werHB2WvBnatQVKUb1zrZO/nN6RU9o1oay8gquHdOLSAR2IiVGc36sNP1dJ65h5yzDeXbyXNfuzeORPvXjr9z3kl5SR0qIhGw/meExVqdDwpVMrwl+HdOKZOdv4YUNlZ85BnZux2ulnumJQMveOcc3u+/Cm013GuXenqLSCpbaJdH6ZPIqZq9L4anUamw7mkjJlNvGxitJy45LaK6Z6Y6a9xeiXyaNo17Q+h3OKGP3871736Y6vDOzYMM7y6I9Jo7o6WkqevbwPfz6tE38d0pnej/7kCLaVgkGdmrk99p6bt4Pn5u2gZ9vGZBWUcDS32CWv/28juxATo1zyvhPiYuib3JQNaTm8PH8nE0d04Z+fG51ve7ZtzLwtR7nhjBTOe2EhjevFcWqHpgzp2pwXfzHO0TG92zJ38xHmbj7CN2sPVkunAmP0nVs/MtIPBnZKculjcWaPVtXWf/rSPhSVVriUHeCdRXsdU6974jw/hF3nFg24bGAybZsmOmaLnXH9adSLj6FerOdW31l3DOezPw7QrVVDbhjexeULrPMoWZe8tsSlpfXqIZ34dMV+RxrZgnvOcuSr7/cyrKlZU//UmytPc9diG0y44u9nPZw7AQbb9329kWfjjRGcpixzvZ9tSMumjcrHVw7lrR+tZu7mI8y8ZZjbFButNb9tT2dQSjOaVBl3HozBJto2TXSkCmqtecPWinlRv/Yu6XJdWzZk5EktmXxejwB+2vAxFXDbguO/Al211o8rpToBbbXW3q/4rl4E7gVHzn4LIFtrbe8NlAbYExI7AAcAtNZlSqkc2/rVjiKl1CRgEkCnTp2qvi3qGHc3Dmep08ZTmLYBPMzAPr5vO6MNxsnATs38Crij1ebzt8M5EwBj5rqkBgkkxMXw691nOdbpk9yUG4Z3obisnLiYGGJjFP9yuog5z2g6b/KZgHHxyy0qI7+4jFvP7MaTl/RxNFv/qV97lzL8y8MFMSZG8YJTzVrzhgn80zbk19VDOnH1ENdz+pqhlZ2Skhok8NXfzyC7oISG9eIoLC3nQGYBrRrVQylFakY+7yzaw+UDk0lqkMCzl/Xl3v9t8DiEn31q7bjYGJ6+tLIW+q1rB3Egs5CC0jLG2GrLerdvwpvXDnKsc2Hfyp83M7+EFXsyHM3r9ho6O+cg7ZIBHXimShrMLaO60j6pPhe+sphfJo+ie+vqTaWjerRi11Nj6W7r7Nm2SSIXD2jPiO4tWbs/m+edRnL46S5jGw+MO4UrByc7amPtwfYvk0fhLZ/EufbWHV8Bda92Tdh9yzjS84o5/OtO8NwK7laDhPCm4E0a2ZXBnZvRuUVDRwDXsF4cv91zFv+auY7urRrxnwl9UUpRWl7B9iN53PvVBsf47nbuZgX9fNJQhnpIN3nt6oG8u2Qv7y1JdQyhemaPVpzXqw0PfbuJ814w/g55xWUs25PBMqehKuc6jepTNdi+b0xPx2g6Gfkl3DvmZG47qzvlFZqHvt1IbIziwXHVp5WPi43h5b8M4OW/DEBrTUl5BSv2ZHLHp2tw7sPZL7kpWQWljiB2VI9WLNxhpKtdPaQTnZs34PozUhydCM/o1pK7Rp9EUv14Rw2/N91bN+bhC43ydUiqz+ldmvPH3kw+vXkIbZok8s1tZ3Dp60YrWXpeMX06NOW/V/Sle6tGLNudwd7j+XwxaahjMIAOSfX5ZetRhnZtweVvLHXkz185OJlpl/V1jKxlbHB7tfJU5T7YNiOIL5zKXLNUqY4lXnmujffVQuRPSRdv3gM04Mq3jPkN6sXFUFxWQa92TbjlzK4ufUoAXv7LAM7t2Zqi0nJufH+lY56I7U+OYcuhXMffFuCSKn1WlFJ8NNFsklr4ma3hfh2owKiRfhzIA/4HVO/54IZS6kLgmNZ6tVLqLPtiN6tqP95zXaj1dGA6wODBg6UGvK4r8V1zUT/O80WqX3JStYDbm2uHdjZ6JfjhisEdIzM/68L/wDkPAvhMazIzm51Syvj5nYRz9Iuq+4yPjaF3+8pBk1o1rsdpTh0Wrzyto9ebo6ebiVKKTi18BwZ2zRsmMNY2myrA5sfHVI6Ze2ZXFu04zpe3DqOhrXnc3ehFp3Zo6rOJNC42htRp4zmSU0TrxvUcgcLIk1rxj3O6u+2A3b11Y1KnjWfr4Vz2ZRQwoFOSxxGRUqeNdzus46zbh3Oyc1+UIxvhTS8FVcZQlG2bJuJ2UEYfDZczrhts9OgJk5gYVa2jKxj5ud/c5jrUW3xsDKd2aMqXtw5jya7jnNOzNTFKsTcjnxYNE9Da6BTpz3wNHZs34NE/9aZ140RHgPz2dYNJiIuhcWIc6w/kUD8hhn7JSfRJbsq2I3l0at6AhNgYRv5nATeP6MI7i/dyy5lduXVUN7YezqVj8wZ0bN6Av5/VjaW7jrMvs4CrbOdAbIzimcs89WZxpZSiXlwso3q0Yu0j5zs+H6i7RgdeIznzlmEurwd0asaGqedz5n8W0KZJIp/8bYij5nTuXSMpKatwmcHzljO78siszY5Utgxbq5e9D8sp7ZrQs21j9mXk8/Wxcb4LVFYCcW6ue8E0yPv72SJz093nk0gSnvuLnNuzNex2/yX63J6t2bzdd2vipsSbmTpwqWNQBvvIUFsO51YLtsF1GFVnJz/kmio67bI+NG1QvUa8JjAbcA/RWg9USq0F0FpnKaXM3FmHAxcppcYBiRg53C8CSUqpOFstdzJgb7tKAzoCaUqpOIyhB93POy6EMz9mBTPP88WvvlPtW+q08ZWJUG7ceEZKZAJuFflOusJw/9hTuH+s93XMjlTUtmn1gNnXNk5p16RyEhwvN/fQ50VGX51IBTHEENw4vg3rxbn0g6naWdXM3/TvZ3XjxuEp1IuLcXzu4v4dqo1I4zxEof3v9NCFlTXVzkPv2V97nq/Qf8EE2lZpkhjv+CLgrF5cbLVKhGuHdmbhjnTHiEsPjjuFsgrNs3O30a1VQ7YezmWrrbViWXwvhsVu4fnSCUyO/8r9zgszoXGAfaCKT8CWbwP7LEBWqt+rxuN75tvzerWB3UZ63V8ucr2f3XZ2dxhU30jy9WHqRb2ZelFvisvKiVWKzIISJn+xnvoJsczbcpQJg5KZMrYnmw/lcr1TatzF/dvz5CWncu5zv5NVUMJNI7qQnFSfa6sMp1jTmA24S20jjGgApVQr8P8KpbW+H7jf9tmzgHu01n9VSn0JTMAYqeR6YJbtI9/ZXi+zvf+r5G8Lv1SEegpkHzeXmnBY6uieFKAmCmXaccSPoBd6W7v9Uyd4f9/nL9P/31AgnSbL4hqQUHbC94phFA0jWdVWSinevm4wBzILXVqu7LO7nigu4/UFu3h3yV6GxRr9Woo6jYLDniJNj9My+i7MEW/TnIRWW5Xh//kRovua/ctO68aJfHxz9dSPM3u0cvul/o8HR4dk/9HC7ES6LwPfAK2VUk8Bi4GnQ1CO+zA6UO7CyNGeYVs+A2hhWz4ZmBKCfYm6YNV7od9m43a+14lmLaO7Q4mIsFw/Zs3wykfA3PdKpxcR/3pRTVls6IeHE9HNW5pYo3px3DumJ9ueqGyKeuCi/m7XtU7oz5O4IFtxROBM1XBrrT9RSq0GzsW4ul6itQ5oQnut9W/Ab7bne4DT3axTBFwRyPZFHbfkxdBv02THlKgz7LZIl6BGiLLBMmq1pU0v5IycH4wXPmvTrP3DnKjfngbF6b5XFHVT/eYEdAxGWetnbNU5CjsOhQMeJurzdDHs4L7DufDOdAShtd6mtX5Na/1qoMG2EDWOryjMTJQWsYhOIsmaJGx/La1h7v1h2FHlTxRYDGJuZlizshr7Ny27qIPiG0L/q71fu0tDN8SggwXBehwVrmdHk/aeVvWsWWeYmgM9LwxVseoEUwG3UuoDpVSS0+tmSql3Q18sISKtFgan3/8T3vXRU0+YokJ4nESsImzHT7A80PklgxNNrQnlklIi3CkrgVL7iB5eDtiv/+bhjeiq4Y6r2mnS3UkYZbXytYXZGu6+WmvH+DNa6yxgQGiLJEQ08hUZRFHk4M3+pb7XEXXLZyGavDeaoucAlMb5P+SjqEPSbKNnHNvq/RhPqz5Nvf/Cd+6Yy+EOslx5R3yvU4eYDbhjlFKOeVCVUs2J3PTwQoSX14BCagRqi1DWWgsnYQzIAxmlRGq4hVultmniR91DYDncwawU+vtKtRxuKx1yP7Z2XWU2WH4OWGYbxg+MDo2hGKVEiNA5IR2fRM2jq9xca3iFsXfupy8L2eYD+dWVxUoNt3DDPqdDw1ZQ4XsM62gXRwVNlHO+ucmzpW3l7Ls+L1IlnifXqYtM1XBrrT8ELgOO2v5dZlsmRPT4v+6RLkF0y9gd6RLUGrU6KA4La36BA2J2mf5Meaz32VdFHVVw3HhMbEpgx2sYZpo0IbbqlO5uL2Je9pvQ0P+drf3I/3XrALOdJusB/TFmiGwOTFBKPWJFwYSwnpeLStWLUEgjqwhHafuXRXb/opayuJ+DxR25yiSHW7gz1zb9R2JTC79hh++eEG8mpSTYn3fvouA+X8uYzeGeBVwMlAH5Tv+EqN1qU6/tDPO1f3WJ1FrXTeVxksMtvIirR60chzuoMbV9zcAcxnzxGsBsDney1nqMJSURIup5u7jUoCjNTJOgCJuq9+WwBP7lkclJja4QxCCdJoVPYf82bsU43FWC4ME3wtz7vPPFyXcAACAASURBVH+ozxWw8Uvv6wifzNZwL1VK9fG9mhB1TTSGEB50PSfSJRDRYvV74duXr2AlwjWBFSo2ovsXUSo2AU693PbCqhzu8B371QJudz9T1XPx8ncC3+EP/wr8s7WM2RruEcANSqm9QDHGX0prrfuGvGRCRJSVOdwRlvYHNGxpzBYm6rYf7wndtkycI4EM2ydE2BXnQXkJtOxhvA73fcCimSb9F4IZlpt2NLG/2s1swC3T1Ik6wsKAINLB+9wpxr+pOZEth3BR9YirzeOB24dA1JE+F5xEUVFEtDi6xXhsG0TDvl9BczgnvvFjpslQatjK2u3XIGaHBdwH5AJtgM5O/4SoeczWHtS2O/Ku+ZEugRBRxM3t8OBqOL4z/EUR0eH4duOxTW/jsRbcA2KVmRpuH/dIf+6hZUUm9le7mR0W8GZgIfAT8JjtcWroiyWEsNzGryJdghqvFtx/w8jXLyuyaSZu01zePgc++FP4CyOiQ75tDO5GbWwLwjwOtwXnRPVhAS2+iBVkWrv9GsRsp8k7gdOAfVrrs4EBgEzrJ0RNbP5f/2mkSyAiqcJMTZc/auA54I+8w5EugYiUQ2uMx3jbCDa14Bt2tZQSr0Lw8/4mk5HbmQ24i7TWRWBMgqO13gacHPpiCRFp7i40Xi4+NfFC3HlEpEsgnFRrnbX6kCo5YfEOQs3iGvAaeAoLi239PvhtBNPx0YJOk9VquGvivauGMhtwpymlkoBvgXlKqVnAITMbUEolKqX+UEqtV0ptVko9ZlveRSm1Qim1Uyn1hVIqwba8nu31Ltv7KSbLLIR7Zi40Snlf39SFMUoucPsWR7oEUanO3H/KisO7P+dfrO18iaZfta5ampC3AIgap+XJ0HGI04JoOmKdtDrF71WrTXyDgr/9GljnRn8vluG+1kQps50mL9VaZ2utpwIPAzOAS0zusxg4R2vdD2Oa+DFKqaHAs8ALWuuTgCxgom39iUCW1ro78IJtPSGCF2UzgImap0aPJBIFqRLVgtwIqlaSNe9HoBQiqpQWQrMuQW4kHPcZ//fhdljADoPg384zEIe4zEW5od1eDWW2httBa/271vo7rXWJyc9prbW9LTPe9k8D5wD2XlwfUBnIX2x7je39c5WqM3VQorZp1z/SJRDC8NbISJfA1QXPOHVOiwInjkW6BCLSirKhfrMIFsDPwNdE5VGFmTkmfIZafoZiMksl4GfArZRabHvMU0rlVn00u1OlVKxSah1wDJgH7Aaytdb2bP40oIPteQfgAIDt/RyghZttTlJKrVJKrUpPl36cwoNmXWDI332vF+x3uta9g/t8uBzZFOkSCAfXm2aNq1UI9uZ88hi4Z0fIimOWUgHXP4naqKLCmPgmsUnlskDuC2FpSTVTw101pSQMfro//PuMQn5dYbTWI2yPjbXWTao+mt2p1rpca90fSAZOB9wlINmPIHdHeLWjS2s9XWs9WGs9uFUrGWhdeHDnOhg7LYAPKkx1mux9aQD7iIClr0S6BFGnRqeJ+KvMVMOkZfz6Td88H677ztqCtHU3WXIdOA6EZyV5gIZ6ja3fl6dA3oqZJlU588sHeN+3pFtawq+ZJpVSk729r7V+PpCda62zlVK/AUOBJKVUnK0WO5nKzphpQEeMDptxQFNABnYU4VVbs5g2fA6XveX36nIZdlVjD4vCrEiXwH/Jg41HKydqqrF/SGEZ+zmSmBTkhvy4agYb4Jr4fBzl7secdysEU7vblRZBfKL/69dC/rahNbb9Gwz8HSPNowNwK9DLzA6VUq1sI52glKoPjAa2AguACbbVrgdm2Z5/Z3uN7f1ftZavXyIUQngYmTok5fAV1YX1qnZkQxh3ZuMySkn4d2+ac3ml01fdk3fEeGzc1mlhtH4xsyrgDqGa9CXfIv6mlDymtX4MaAkM1FrfrbW+GxiEURttRjtggVJqA7ASmKe1/gG4D5islNqFkaM9w7b+DKCFbflkYIrJ/QkRgCAvrNF6XRY1hqV9wxcF1Cjpg5nyatv/JofmtJDXzU/rCFn7LN2/iDJ7FxqPgQyX5ywc36RN7COe8vDeni6dbjyWFoRzr1HJr5QSJ50A5+S/EiDFzAa01hswZqisunwPRj531eVFwBWmSilEtJKm66hXJ/5E0VLbFDW/az8mujq+A5p1DktpRBRY8JTx2LCl9fsK40UnNuw13LZ9nTgKLbqFcb/Rx2y37I+AP5RSU5VSjwIrgA9DXywhAtMYq75F+5j4xsx2ajBJ5nJVY/+a3c6OdAnMs/Lg8+fc/vJG6/Yvok+nM4zHps6N+GG+APp9zJur4fZo8EQYdkfl61Dc8w6vNx5XvRv8tmo4sxPfPAXciDExTTZwo9b6aSsKJkQgRraNwJBHUEeqRYWVwvplRobAq0JVP4Wrvi7Jk2+cdYlS0Hm467Jo/fuHqtPkhc/DBU+FqFA2Q241HjsNDe12a6BArrp7gWXAWqCxUmpUaIskRGBWPTSa5yfUhPGvJTgXvll6lBRkhH6bJr50Rl3YonwM+2l3eJ3lRRFRIv84NGheZWEgR24wR3voa7jjVJnvlfzmxznTwDZtyuy7Q7jfmslUwK2UuhlYCPwEPGZ7nBr6YglhXstG9UiMMXExadoxdDuPrWdi5agLN/yma3DZ/VUnvg6V5Ed2/7YaOVO/a0tbkfwMuA+ttbAMIqoc3w65h12XhbuGe9lr/q1nqoa7wncOd8sexmMoZkZOaGg8Jkn/B7M13HcCpwH7tNZnY3R+lGkdRfTQFf6vW632whcvN+QRd5nclqgNLB1JxEpbLZ5Ixh03vytTo5RYyl0A4qZsP/zL8pKIKFBeajwerzLzqel7RpA2zAz5JuMo8x1wdz0T7lgNA64JfodKQaue0LZP8Nuq4cwG3EW2UUNQStXTWm8DTg59sYQIUKQqYO3f4oUIUF1oPYhebnK4Rd2Vf9x4HD3VdXliU3jwqLltBVUrbk2nSb8O9ZbdQ9eqVL85FOWEZls1mNlhAdNsk9Z8C8xTSmVROSOkEJG3+r1Il0CIkKh5AaD5AkfNj1jzftnCSifcTXpjE87ZEv0N1k3E9BGZ+CaxCeRKqOh3wK2MttN/aq2zgalKqQUY06zPtapwQpiWkxaa7VS9AYfqhuxv56xoJZWwImA1+LgXdcvGr4zHhq0jWw6/mR2lJMziG8jEN5gIuLXWWin1Lcbskmitf7esVEIEykwOt1lSC1YnmMnLDuUREa0jjkUPqztNVl0k53udtexV47Fhi8iWo7zE9zomxVGOUiG62Ph7jiQ0hBIJuM3mcC9XSp1mSUmECAXLAu6QhlYh3JaorZRVAaY9su/7Z2u2718hIrhvN5SquR1ghXXqh6KTZBDHesFxP3dhZljACKSUaA15klJiNof7bOAWpdQ+IB8jCtFa674hL5kQZuUdhdRFkS6FENGtrMh4TOoU2u36CljdvR81QW4NT/US1qjXJNIl8JO5TpMlpkO/IG351ngsL4XY+PDuO4qY/a2PtaQUQoTCgRWRLkGtF2X1kiIQ9rGkd/8a2XJYzszR6mbdtZ94WFVH0RcFYYnWvaBZF4iJ8IysDVpCkR/rmajhjqUc86FfkEbcBb8+CcV54R9aMYqY7TRZobU+YGF5hAhcSKertuqGKjVp0c7MXyeUcVfYvszYJ9PoMRYOrg7XXqOcmz9k5m73q5YVh3ekChF+RTlQPynSpcDvq4KJVMp4yikKd9VJ4/bGY3FunQ64/Y5QtNYaYzhAIaKT5bVOtTRQLsqFZa9Lr70oY9nhvO0H4zGxqUU78KTyBwroSDP5CzG1upkc7gVPmiqHqIFyD8I6Dy0cAFNNjCnt13XVw7Hn97CA/gfcjmEBOw+H+w/6/bmg1GtsPOb7mZNeS0mnSVGLWBgQ1+Ym5HkPw0/3w46ffK4qMXktEpcQ6RJYyrKOYUtfsWa7IjrYZ5nseWFkywFYUcNtzDQJxMRBvUYBlaqSyfuiBaOu1CSBdJq8VSmVinSaFMILDxcipYi6TOhSW5JgYWZky1HHhf3LTO7hMO/QSSA/rKW/oFr8hVqYU5hlPHY9K4w79XBs+xtIm0kpUeWoUPVD8HcbDVsZj2X+JKTXXmYD7jFUXpmiLGoQIoRqc412VfZp6YtyI1uOKGHmT1+jh5Irzgvv/tz8rrSFga6pLSuZ2l3YFGQYj9GQa+z3zO7+B9yxVNhaf8J4wNv7PJTW7YDbr5QSpVSeUioX2ARstP3bBGy2PfpNKdVRKbVAKbVVKbVZKXWnbXlzpdQ8pdRO22Mz23KllHpZKbVLKbVBKTXQzP5EHRLSHmxVr3TK5PZr0PdR+wxgi5+PbDmEC0tuh3lHKp9bOUlUFDCZ8S113MJgD7hDMgY3BHcv8DelpMp6N3qeANyRUhLOb5hxtoC7rDB8+4xCfgXcWuvGWusmbv411lqbHaiyDLhba30KMBS4XSnVC5gCzNdanwTMt70GYyjCk2z/JgFvmNyfqDMkhzsg7fobjykjfK6qa9IXif9v78zD7KjKhP97b3cn6ex7IAtJ2A1IWJp9FxkWgYAMAqMQEIYBRSU6jqifgs4ziuKMz4gLgqDgh4gKKktUMA6DyiKBL0DYQ8AkJISQkLWT7nT3+f44dXOr69a9t+p2bbfu+3ue+1TVuWetOm/VW6fe8x6lnBfvK+2b3mjzrkNGwo5Ch8s7xr7a0xVf3kq6dDqmdUOTXGUy4kmTo6f5RuuVFtqKkyaTfMUsKtxpmrFlgMSdTBpjVhljnnb2NwEvAlOAOcBtTrTbgDOd/TnA7cbyODBaRHZOuNqKkl9aHMuySN0qKmFJ/GWmL2KFu+EJoYB0bY6vGkq6bFtvt5lwCxgUz72jwr28jxZaceQ+yUGkbY5Xlwe/mFyZGSTVJ6yIzAAOAJ4AJhljVoFVyoGJTrQpgNv39wonzJvXZSKyUEQWrlmzJs5qK81AnDejrLn6eOAzdrt1fbr1yAgNbZcdhqhHuGuS7HkNNcId1oZ77auh66M0CEXlMHG3mX7U66XEvzMbKSrcCT+DxkxPtryMkprCLSLDgbuBq4wx1WZr+fWcst5ijLnJGNNhjOmYMGFCVNVUGgn1wz0wljyUdg0UN3F0t27XyGyaI9wJvHSGf3EKEf/Wk0LmrTQMT99ut4MG6jIP6O6En541gPQBv6R4Fe4Kfb/PUbjtvwk+z9rHJFdWhklF4RaRNqyyfYcx5h4neHXRVMTZvu2ErwDcBklTgZVJ1VVRlBJZG5xXQvLHa0v7kU+aDP8ANzG+JMdpH67kmKI5RqFl4Hm9uXDgeQShTOGuYFIirbRKHwX60unzM49JvswMkbjC7SwRfwvwojHG7RbhXmCusz8X+K0r/ELHW8lhwIai6Ymi9CNWrwtRjZZl0A+3kg2S7hZJj3AP+AEfr4IQp4tCpYEYvydM2DuavMRHad/1uGjydhPQpKTPUcStWUkE/T3MCMxO74W2oQMvs4EJ64c7Co4ELgCeE5FFTtgXgOuAX4jIJcAy4Bznv/nAqcASoBO4ONnqKg1DI0wCq6V0DBoB3Qn7R3bTux1a2tIrX0mOxG24+xUeewnhvJQIIiHr1NtTmnCs5IfNb0enGBZ8+kchhvur99lXyaQE+wIwiJ6IRrjDyRh/fzSCMhuXxO8Wxpi/UPnV6gSf+Ab4eKyVUvJBpCPcPl100qxoss6yXUb3lgabnZ9fYvcMnQE/3JkZU5Y6znbXxmwsjqJEy/LHo8srCrOUWhhDYC8lYlW+VEa433p24OU1OOoHTMkPcY/YtbXHmz+QurlJDSUsw68KDU/i5zbqL0I1R8wy7KXESRGK4hLgSv6YfEA0+fgq3BFLegg57nNMXKYX3oa+nggKD9GWg//ZSdK8TxFVuJX8kCmTkga9qWTqHCqxEskDt4npStH0S4mPIaNh6sEDy6OvF7Zv9bfh9mMg5h1+clxjhBuIxitVmK9koxxvztubd7VJVbiV/BDnJ/JIZ3RXU8ZT/siuSljzkKJJifFsAxF6pclwsUO7Eeyq5s1WaUj6+koL3wyEey6D/9jJf4Tbb4R3IKO+G98sD6vQl03Ui5uFqXfRr3kU57dBUYVbyQ+1FIgxM4Pn1axuwmqY5Zgm/hyYNLF3waSvpU+DYrdTD0o9J3v9sujroaRLUXl94bfV49Vi8a/s1m/SZNQ88GmfwMp+uCMlzEv7oBF2290ZbR0aCFW4lfxQyxxi+pHJ1KMWWVZadYQ7NZJ/mYm6vLgneYarbyGs15GwrFxUO47SWPRss9sTvhxRhoHW7RsYrT5zi4KYlCTB0HGl/UHD7DboYj45RH0aKfkhVTdnQcnIiF4l1Ia7eciAl5LMUM+n9gl7RV8PJV2Ky7oPi2i16qBfTqL+nFXRpCTqEe4qLw9XPNb/PO5QuLdEW4cGQke4lfxQU1kcyMjCABa++dDt8MGbB1B2yPIGgircmSH2VzP3wzKKZaxr4mpRAoP5obyUGKM23Ap0rrVb98hs1ET9Jas40NTPv3cGTEomzYLhboXbucc08WRjVbiV/BDliN1+H4our1lzYNS00nEcxrlRKUw1bbijKUbJAKmOcDsdKYwoxG/UHi66ugXMH5Er3Al80SzK8U7vdRWblElJCJkZ7DyjmtikRBVuJT9EqUAUZ1THQZa1VrXhTo3Ye8XLv/MUGLHCnbGJxoUwyo5I+Oqrwp0/tqyx22HjYywkYi8lxbTuDpyYl5IQ95DWIXZbtJNvQlThVvJDpm1STYX9OPJXFB/uPK//8UEXlfaz/BKYGCE17s518VRDSY/Nb9tJiImYWFUjRF8sPvfcynRFxTril+Iwz9ziwnFqUqIoOaCW8A9EqYhq9C6uUcCoFCZVvDJDaJvisAyf5DpIwqjabcMdf3nhVpoUQisj25vXvVlu2bYe2kfH+7XGr+97ywszEr30f4qJ3Bn4Fx11u8Ksvtw62G5/f3W0dWggVOFW8kOjKIvV6pn2Z/mX7k+3/CYm+e7bIPKSBPXIXaa/qCl10bUZBo+IMMM6Zaye/hjApCTyEe4xM4LH9XNf2GSowq3kh5oPwCS9lHjKapSXgUV3pl0DJSn69ckoHsRxvyyGzD+sH+6w1e9ThTt3dG+O1pzE/Uy69E8w7TAifdF9/RHXgVvh9lftTJpuaVsHQcsg2OWI9OqQMqpwK/mhIUacJJ5R7H0/GH2eSn547Ps+gXHPK/BS6veZe/00htAad0P4/VdC0bWp5E0jCrwDLYHv/QHj3XZ6af9o14qTFcp5fdzRAcuPid5uWPZounVIEVW4lRxR4zGe6ChzwiMJs86MKKPq56hRBurzQKTvZf/ztfKwVC+mLTvOpd1D5SyhfJpY1Gd9/ti2MVoPVe5BoGIHCyJ39Qj/zrOr/3/+Xaxrnxk+XyUyVOFW8kOcI9xRaj+pz4BXskis6m+3n2cAV4kpKd+RT+Kql7psuFXhzh1dm2BwTAo3UPE1sEz+6pGLGmmkjonBcdGkIzeJK9wicquIvC0ii11hY0XkIRF51dmOccJFRL4jIktE5FkROTDp+ioNREOYlFCyrxu7a7r18KNJb4RNiftaR6H41sojYeW6EPYVJmz9dIQ7f3RtDD5psm1YgEhB+2C9pifuNLXVudg9HwWlST38pDHC/RPgZE/Y1cACY8wewALnGOAUYA/ndxnwg4TqqDQiaSvcVy+HL60NHj9K5bYQ8ZK9Sv7pJy8JP4jr6fqxKgt1GLekfb9RoqW70yrcQVeZvPA3tePEbcPdL4lHnbt6ef/jvp4YZChkfnufZrdN6sM+cYXbGPMI4D3bc4DbnP3bgDNd4bcby+PAaBHZOZmaKg3HC/fWiFD9Kf/MyOMAWMuY+sofMhJaBrh0bktbfemmHjywcndQw4Y7e9PdcoPxPJzjV4HdJiX5Ux5D+eGux4Y7h+esqXn3DbsdG9DOOdBqlD6egOKy4famGTKy/3EWvsjsd67dNukqrVmx4Z5kjFkF4GwnOuFTAPdr2gonrAwRuUxEForIwjVr1sRaWSWjvPXsgJJva7GfErfK4ChqUx8Xza8vXVY+FSqNQ78Hf7JeSjKHSHh7clW488WGFXY7enp0eSZpw12r//Ztry/fKBk2wW43r063HimRFYW7En69w/fJYIy5yRjTYYzpmDBhQszVUpqPhFaanLh3NOXUi9pwNxFRX+tsKdSFuKujspIvtm2w2/bRARME6GC+L2V+/SYBG+7uTiR2oahB0Vxn3evp1iMlsqJwry6aijjbt53wFcA0V7ypwMqE66YoA2eXw+DgS+GsG2MqIKobqSoRTUO/S50tZTkKRrE5VHw1KWlytq2328Ejq8crEkQp7u0JliaKl7daCvemVekufAMwaKjdPvz1dOuREllRuO8F5jr7c4HfusIvdLyVHAZsKJqeKEpoBnJTO/CCcPELHlvuQgt84D9h9C7uCtVfn5TQQb34KPuoHPuzMeGLmbDZ0zmyIHhkkXIPDrW8CKnCnS+KI9xR+uF+6selfQlhwx3FpEkv7XXOTapaZsh6jnQsgpt0obY03ALeCTwG7CUiK0TkEuA64EQReRU40TkGmA8sBZYANwMfS7q+igKEv1m1tdvt7PPL/wt7k7p2Q7j4A0E16uYhcRvuCux9Gux5Snrlg7+yMqyWaaLKSq7Y+i60tkPbkIAJAtzHtwSdTxaBSUm1+ux7Nsw+n0E9W+rI1+HAC+tPW0QEho63HlOakAG6VAiPMcZHAwHgBJ+4Bvh4vDVSlCDUOToX1KeroiROxAvf1DuCfd4dAy/bh3CO/upZ+EZHuHPFtvXhBlYC9fcq09B6t8Pdl8Ax/xa8zGpUcw37j7faKDKAPnvGDfY3UIaOU7eAiqIotdAxvRyR4teMJNxLhrJXFQn/vrBlLbz+SMhESmbpfDd6s4vBnlWF3Z3s7Rfghd/Cr/8F3nx64GVJ7bUYDBlYr0EVbkVpAtL2ANIQqEqdGmWnPmab5+mHx5t/GaX2iKPsxz4KHSLv0Ll3bYDbToee7jgqpCRN96aQXyQD9JhJ+5TH977orl4MD30pRLkVCLL4WRbcxw4dq364FSXXzD4fjpw38HyOvGrgeRRvvN6JlVlAbbibB/fksJQexKl7TSgykParaUk+6NpUPiJdjUp9pttlJ93n6RtvLYYVT9r9qvfaGPxwQ6Dl30MRdFVON+1jYKuOcCtKfpnaAYUIuvvO+w08jzEz4KhPwz/9Inzaoz/jHx6VwpSFERAleSo9/CcfECKTjPWdUNXJWN2V5OnaDINCKNyV+szXJpf2jWd1x63rXGFVFO6g3bEQdmXiiPv5YXVMsRsyquQRpslQhVtpPIyBb+4GT94SccZVbkZRKqIi8P5rYNxu/cNH+i6i2p8TvhxdPXyp3k7v8uNKdHjtmjPx7hP6gV6FTDSoAvWsNLkDlYmGp6cb1r4K2zuDp2kZVDuOdzn1KR2l+7x39LsfQRbVMXb1yMM+Bl9aWzs+0Gq6AsULTEsrjNg5XJrBI+157t0ebV0aAFW4lcZj7RLofAce+HTEGaf84Jw1J93yIdtKkRIfXjddM46227j7Q2a6m09Fgr5c6kto4/PyfLt99cHgaVoCvIy6R7hFYPhEa8MM1U2Rgoy0F18Oho61iq8fnomUU959sna+YfnMS+HiF1fybMJRblW4lcbj5jIPkhFTRQvI/cM1MxqQkiTeCVdRLv4RFXEq/yI+0yaDynre7wlNwIbl4dN4Fe4Rk8vjeF9kpVAa2faam7jZuMJ1UKF/FVdrfPH+yvnMWwz/UvKkY+JS+Q66KHjcIY7CvXV9LFXJMhmctaUoNehqvjfjxNAR7iYlguteq++4/k9CRR1BCPMAhLrPgU6abHw6g5lk9KNlsN3u80FY9QxM3r88jtdsRAql/hKw3xQw9Pn1zVXP2u32rZUTj5xsfw4miCeTgXL5X2Hd0sr/7xjhbj6FW0e4FSUQzg0v9wppDRvuhGrRjHg/niTa00ZNTbK0eNjvvH6H2whgY1tECvWLdu6/ejUBm1aHT9PSCvNegLNutM8Fv37Qz6SkxaNwB+s3Uumu+/r/2m015dZbHbeJyf4fCZyudsaul4ed9oVZZ1SO28Qj3KpwK41Pd2fpbb8iDaAoD+jBHbGXkr5eWPFUNHkq6fHGX+CRb9WON3qX+OsSNx6NOZyTkoHIjyrcDU89I9wAo6ZA62Bsb/PpB70uH+2FFvsL4qXExWBqTC4cNiFQPrZEVz8/7PLA6Wry9O3B4xbN1ZrwS7WalCiNTec6uPtSeG1B2jWJgCw8uJ0b8hM3wh++ABf8BnY7Pt0qKfXzkw/Uly6LX3IGjwwVvUA4U49Cn0exCTxpUk1KGp7OdwaW3j1y7ca9KFKhtX88rweTCoyTDXSaIZUjnPvTwNU0AVajjJ2iwq0j3IrSYHxzZrmyfa3Pm3PUDv/jIAufpot6VnHxBs+5zUIV80rjn9swSnodjd1p3xr/v7ffYViFu22bZ5QzqL1r4184pTjCPW6P+tIXWv0VaPcItxSsWUkxnvcFrwI9psK46L5n2+2UgwJXs5/CnVa/bR9jt2rDrSg5JYnJIvWwUwQL6UTJzGPtdth4u93wZnp1aXIkyVHmshX2kik70lJaB/c7fJTZdqdtWKCalNvKBqydKtyNT7czwbZemSu02JUqX32of3ivy+/1jhFup78EHOGeImv8/ygOioSps3vxt7TmbbQNgdb2plzeXRVupTkYNS18mnnPl/bjGiG//M+ugww8uIujJauesduNqnDnmnG7222lFUyjYvT0/sexKKkVFI+RHndtYRfq8NLa7gnIgNwqA6Oap48gtLTZSYx3/GP/cK9JScFlUhJw4ZeLW3/v/8crFcKr0M8tYNEfeBSMmWG3+3wwWPz2MdYctMlQhVtpDirZIQ/fqXIa9whA2jatly6Ai+YnUJCjPDz1E7sNMQNeaUD2dmy8g6xyOhCGT/QN8DrGKAAAFDVJREFUDr264/u+VPk/T14tRZMS79etqQf7pCX4S4D35bsJV8zLFV2boHuT3a/3RbBSup5tpf1C0UtJ0aSkxz8NwOCSH/w/9B4SvLwa1L+aag2mH2m3uwdcI2PoWLXhVpSmI+iEk+KNKuikk/0+BDOPgaPmBa9LtdnmUztgxpHB86pZlr8CVDbqsGVN6XMroKN58eFd2j0RTrgGPvNKuUI85UDY4x/gtG8Hz6vawzzMF6Kh4yr/99Zz1SrQ72iHwu2V2YPmlic1IFKnwh3QFlfJKL/73MDzWPm0f3jR7AMck5IW16RJH4X7gzfb7fQjdgRtYXB5PHe+IZC4J00GfRFoHwNbdYQ7k4jIySLysogsEZGr066PkiJBR5OO+wIccEHteNNKoweBXv73Pq204EE12sfA3PvKP2f7UVxGO8TklzKCLDPspuCZiFMcUWn1aZuOcqdCIt9UCi0wYlJ5eMtg+PAvYdI+0ZRTQeEuX90R6HUpIjOP8fxZ5YEuAmf+AID7ew8tuVN7z+mlOBf8BnZ/f3navT+AeJWF5Y/3P77w3lI5+5xVCr/vUyUTLKXxWONamnx8nZMmK9G9ubQvBfsrTtD0e5bt9j747FLY75wdQe10l8er0yuXics0svjwDOqxp3202nBnEbGvZN8DTgFmAeeLyKx0a6WkxrLHa8cBOO5zMOe70Zff0gonfiXaPA+5zG6Lyk21Eb5K+E0KnfTe8rAi+3ps7QaPsNtun9X53A8kJR+017DfjPrBHCa/4gTO1iH2pdVN1Zdo2WEis9EM4zq5wQa7Rgt3rHLnZcZR1Px6U1TGdj0W9nK5W3ztT/BD74uB0jC8+/fS/lk/jDbvLpfCXWiFJX+0I9vLHoftPvfalkEwbBzMOguO+AQA+xbeKP2/fSs8/xv4xYX2eL9zQ1UnNoW7OGg0Ye9g8YeMhm3qhzuLHAIsMcYsBRCRnwNzgBdSrZWXxffAqkXOJxXT/9NKv7Aq26pxKR1vWmlnOA8eAeteh5lH20mBpteG9/Xa/e2d8M6rdmJUcSWsHfn1Bdx3yt2+FYaMdHzhuuvq10b61939f+i4rv+XPQbvvh7uuhz5KRi/V9Uoy4fvx7TNFRbOGTUNxu/ZP8z9AI+CWWf0d2V41XN29OMb02Ho+NrpOz7qH37m9+CBf4UVfyv/b/gkmH0+PHOnPR4z3bqwevhr8OZCV7yd4NeX75h9P3FbN99uc1Zlu/vXdlvrOvvu49n36w8+++68erfb35gZnpXeKuRfI+yGtpV2vPWXvyyvp6fu3297y3q1uOtO/8+orrCb296yOz/7KTcXz93Pbi+Ld+Wq9ZzbZr0aCIaZfxgDQ1y3aG+9+3qsXE7psO3v6ynJfl8PLPe57kUO+1h52OfegAVfhYW3lrzUhMH71QRg+lHw97/A2JnB8znwQnj463CGzwvzHifabfvY8k/SwydZhQX4p9Y/lcLdJh+TDyjtH3IZ/O0muz9kJGJqeI0YNgGueAzG7go9PpPsrh1lV7v0XXXQr/+l8Z8nXjUTgMDpAv5XT/2NsV8P9jzZ6ZOu+0HZM7LPJ6zG87Svt+SDe+rB9hlXDyOnwsYVpeOOS+w9c8OyUljrEPt1ZeEtcOtJ/vkUBz4KBTjms/DoDVzeeh+Xt94H1/rEP/07oaoZmx/u2efZr1FBvuiC/QK8aRXceLQd1Y+KKQdVX+UyZRpB4Z4CLHcdrwAO9UYSkcuAywB22SWFVdNe+xM8+wvn04p4tnjCvMfV4laIs+Ude4Pp7YJ3Xqbmp5y2oaX0UijVY8e+VN/v6bJ+M93h/erqqTfFTdC4lc6Ba3+DqxuMmAzn3wk3HQufeNquqHffJ+HIq/q3+8SvVj8vQPdJ18PdJzHu+Cvh7vt466BPs2O65LzF5Ql2ng3H/x944gf1r1BWjUGOG7OrFvuPyF3yR1jzEivu/xqjCl2MOOEaG37Z/9oVAxd8BQ69HCa+B86+GVa/YFdEW/qwfagsWWBtzPt67EvZmlfg7Fvg3Tfgse/Cmpet8nL0v9rJLQu+CsseBYThwIGylZaWAqxYWePa1tgH//5SFu7dd8XbsNy+gO7AGy942FHDt7N+63ZYva48nqfsQ0d08W7ndli7pWo8gP1HdrFuy3bY2M3O4kzO2thTFm/mkB76NncydUw7K97dyrDeAnRK5XqvW2r73+rnSy7HCi0lW9F+58XFsZ+DYz9bHt4+Bk78d6tQBvU24GboWDjzRnjpfvu5+JRvwoS9YNHPYJ8z+0Wdcvg58Np/M+EIn+WlOy6xirMnzQ4u+I0dRFj1DPzxGjj84/ala+ZxUCjQe8YPaLn3ilL8GcfYl8sD55bCrt1gla+iwg1M3OdYKOrpwyZaGVjt2IyP38uabU1yPq62DYFPPQv3/DMsf8JJJPDK70oLe3jNZfyuYxr/ldnOBf2PKv8FzTNk/Y2xL1bP/MwqrJWeiVWflzW2kw+AsbvBSV/zNjA4p/0XPP4DqzDvdSrsfz78/TFY9H+tud7sc6F1EJx6vZ2Mv2WNLf+NR6wr1omz7L3efT6GjIKTr4PfO1a0R3zSTmYfOtY+12bNKavGS22z2DR0F3ymBQMw/fCz4cWv19/OagRVtsHK46Pfgbeetc+bqDjgI5lWuMVk3IeoiJwDnGSMudQ5vgA4xBjziUppOjo6zMKFCyv9nR/6eh1leJtVnsR52BZnQ6ftWUNRFEVRFCXHiMhTxpiOWvEaYYR7BeB2ojwVqDB002QU7XbbvH5hFUVRFEVRlKyQ+UmTwJPAHiIyU0QGAecB96ZcJ0VRFEVRFEUJROZHuI0xPSJyJfAHoAW41RjzfI1kiqIoiqIoipIJMq9wAxhj5gNJLLOnKIqiKIqiKJHSCCYliqIoiqIoitKwqMKtKIqiKIqiKDGiCreiKIqiKIqixIgq3IqiKIqiKIoSI5lf+KYeRGQN8PcUit4FWFYzVmMwCthQM1ZjkJe25KUdoLKSRfLSDshXW1RWskde2gH5aktasjLdGDOhVqRcKtxpISJrgpz0RkBEbjLGXJZ2PaIgL23JSztAZSWL5KUdkLu2qKxkjLy0A3LXlkzLipqURMv6tCsQIfelXYEIyUtb8tIOUFnJInlpB+SrLSor2SMv7YB8tSXTsqIj3BEiIguNMR1p10NRso7KiqIEQ2VFUYKRdVnREe5ouSntCihKg6CyoijBUFlRlGBkWlZ0hFtRFEVRFEVRYkRHuBVFURRFURQlRlThroKInCwiL4vIEhG52gmbKSJPiMirInKXiAyqkPbzTrqXReSkanmm1A4Rkf8QkVdE5EUR+WSFtHOdtr4qInNd4QeJyHNOnt8REUmgHbeKyNsistgVdr2IvCQiz4rIr0VkdIW0vuc96PVMqC37i8jjIrJIRBaKyCEV0mbmmrjKVlnJ0HVRWdkRLzPXxCk3F3JSpS0qK6isRNSW3MjKDowx+vP5AS3Aa8CuwCDgGWAW8AvgPCfOjcAVPmlnOfEHAzOdfFoq5ZlSOy4GbgcKTryJPmnHAkud7Rhnf4zz39+AwwEBfgecksA1OQY4EFjsCvsHoNXZ/wbwjaDnwPmv5vVMsC0PFs8jcCrwcNavicpKNq+Lykomr0ku5ERlRWVFZaW+n45wV+YQYIkxZqkxphv4OTAHeB/wKyfObcCZPmnnAD83xnQZY14Hljj5VcozjXZcAXzVGNMHYIx52yftScBDxph1xph3gYeAk0VkZ2CkMeYxY3v47fifh0gxxjwCrPOEPWiM6XEOHwem+iT1PQfOm3qQ6xk5fm0BDDDS2R8FrPRJmqlr4qCykrHrorICZOyakB85oUq5KisqK1GQJ1nZgSrclZkCLHcdr3DC1rsEsRiGiJwhIl+tkbZSeJxUKnM34FznE9PvRGQPABHpEJEf1Ug7xdn3hqfNR7Fv4IjIZBGZ74RXasc4KlzPlLgKuF5ElgPfAj4PDXFNVFayeV2qobJSHh43eZGTavVRWVFZiYI8ycoOVOGujJ+dUotPmAEwxtxrjPlylbSmSnicVCpzMLDNWJ+VNwO3AhhjFhpjLq2RNo12VEVEvgj0AHcAGGNWGmNOLf7tkySL7bgCmGeMmQbMA26BhrgmKivZakdVVFZSa0te5IQq5aqsZK8dKisZuS6qcFdmBTDNdTwVWAaMFpFWV5jf5xm/tCurhMdJtbrc7YT9GtgvZNqpPuGp4EzwOA34sPPZy0uldrxDsOuZFHOBe5z9X2I/gXnJ4jVRWcnmdSlDZSXVa5IXOalVH5UVlZWBkidZKWESNBhvpB/Qip04MJOSgf0+2A7rNtr/mE/afehvtL8U+3bmm2dK7bgO+KgT5zjgSZ+0Y4HXsZMoxjj7Y53/ngQOozSR4tSErssM+k8IORl4AZgQ9hw4/9W8ngm25UXgOGf/BOCpBrkmKivZvC4qKxm6JnmRE5UVlRWVlTrblWRhjfbDzuh9BTuz9YtO2K7YWbtLnIs/2Ak/AztZpJj2i066l3HN6vXLM6V2jAYeAJ4DHgNmO+EdwI9caT/qtHUJcLErvANY7OT5XZxFlGJux53AKmA79m31Eqdey4FFzu9GJ+5kYH6t817peqbUlqOAp5wbwRPAQVm/Jior2bwuKivZuyZ5khOVFZUVlZXwP11pUlEURVEURVFiRG24FUVRFEVRFCVGVOFWFEVRFEVRlBhRhVtRFEVRFEVRYkQVbkVRFEVRFEWJkaZVuEXkZBF5WUSWiMjVTtiVzrERkfGe+G0i8lSltCHT3yoib4vIYk+c/UXkcRFZ5KzU5ecvU1ESJSZZucMJX+zIQ1uF9CorSsMQk6zcIiLPiMizIvIrERnu+m9nEXnQ2f+9iKwXkfs9ZZwgIk87svIXEdk9znOgKEGIQ1ZccW8Qkc2esPRlJWm3KFn4YX0yvoZ1MVP0xzgLOADrw/INYLwnzfHADZXSOnFqpnf2jwEOxOUr0wl/EMeFDdZ9zcNpnyv9NfcvRlk5FevXVbCurK7wpnf2VVb01xC/GGVlpCv+fwFXu44vBj7j7J8AnA7c7ynjFeA9zv7HgJ+kfa7019y/uGTFidcB/BTY7Emfuqw06wj3IcASY8xSY0w38HNgjjHm/xlj3qiQ5mSs03fftAAB02OMeQRY5xPHACOd/VGkuzqVokB8sjLfOGD9qk71Sa+yojQSccnKRgAREaCd/stRu2VlAbDJpwyVFSVrxCIrItICXA/8W5X0qclKsyrcU7BO7YuscMKqcTzwcJ1p3emrcRVwvYgsB74FfD5AvooSJ7HKimNKcgHwe5/01VBZUbJGbLIiIj8G3gL2xo7yFZWLvYwxL9Qo41JgvoiswMradbUaoigxE5esXAnca4xZ5U6YFVlpVoVbfMIqrgAkIpOBdcaYzrBpfdJX4wpgnjFmGjAPuKVGfEWJm7hl5fvAI8aYP/ukr4bKipI1YpMVY8zF2FUOXwTOdYIPxa4cWIt52OW4pwI/xpqlKEqaRC4rTpxzcF5IPWRCVppV4V4BTHMdT6X6p4NTgD/UmdabvhpzgXuc/V9iP50oSprEJisicg0wAfh0hfTVUFlRskaszxVjTC9wF3C2K737y1AZIjIBu7x6Udm4CziiWhpFSYA4ZOUAYHdgiYi8AQwVkSWu9KnLSrMq3E8Ce4jITBEZBJwH3Fsl/g7bnzrSetNXYyVwrLP/PuDVAGkUJU5ikRURuRQ4CTjfGNNXIX01VFaUrBG5rIhld9hhw3068JKT5gRgQY06vQuMEpE9neMTsaPkipImkcuKMeYBY8xOxpgZxpgZQKcxpuhlJBuykvZs1bR+WM8Gr2Bnu37RCfsk9u2pB/tA/xF2RuyiWmlDpr8TWAVsd+Jf4oQfBTyFnXX7BHBQ2udJf/qLSVZ6nLBFzu/LKiv6a/Rf1LKCHRT7K/AcsBi4AzupawLwJ0/6PwNrgK1OeSc54Wc56Z/B2sDumvZ50p/+4niueOJsdraZkRVxClEqICJHAR8xxlyeRnpFaRRUVhQlGBHIykeAqcYYnQCp5Jo8yYoq3IqiKIqiKIoSI81qw60oiqIoiqIoiaAKt6IoiqIoiqLEiCrciqIoiqIoihIjqnAriqIoiqIoSoyowq0oiqIoiqIoMaIKt6IoSgYRkXEissj5vSUib7qOH42hvItEZI2I/Mh1/F1PnIdFpENEnnDqscxJU6zXDBEZLiI/FJHXROR5EXlERA4VkXYnTreIjI+6/oqiKFmmNe0KKIqiKOUYY9YC+wOIyLXYhRy+FXOxdxljrgxQt0Odel0EdLjTiMjPgdeBPYwxfSKyK/AeY8xWYH9n2WVFUZSmQhVuRVGUBkNENhtjhovIccBXgNVY5fwe7EppnwLagTONMa+JyATgRmAXJ4urjDF/jaFeuwGHAh82xvQBGGOWAkujLktRFKWRUJMSRVGUxmY2VsF+L3ABsKcx5hDsssifcOL8N/BtY8zBwNnOf3GwD3YZ5t6Y8lcURWlIdIRbURSlsXnSGLMKQEReAx50wp8Djnf23w/MEpFimpEiMsIYs6lKvpWWIdbliRVFUUKiCreiKEpj0+Xa73Md91G6xxeAwx076qCsBcZ4wsYC71RJ8zwwW0QKRZMSRVEURU1KFEVRmoEHAffExv0DpHkSOFJEdnLSdACDgeWVEhhjXgMWAl8RZzhdRPYQkTkDqLuiKErDoyPciqIo+eeTwPdE5Fnsff8R4PJqCYwxq0XkU8B8ESkAm4HzA4xcXwr8J7BERDqxI+WfHWgDFEVRGhkxRs3xFEVRmh0/F38xlfOGU0410xRFUZRcoSYliqIoCsBW4JTiwjdRU1z4BmjD2pcriqI0DTrCrSiKoiiKoigxoiPciqIoiqIoihIjqnAriqIoiqIoSoyowq0oiqIoiqIoMaIKt6IoiqIoiqLEyP8Hyvlt4VfJqtwAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig,ax = plt.subplots(figsize=(12,4))\n", "\n", "hours = HourLocator(byhour=[0,6,12,18])\n", "days = DayLocator()\n", "\n", "h_fmt = DateFormatter(\"%H:00\")\n", "d_fmt = DateFormatter(\"\\n%x\")\n", "\n", "ax.plot(time,lw,label=\"lw\")\n", "ax.plot(time,sw,label=\"sw\")\n", "\n", "ax.xaxis.set_minor_locator(hours)\n", "ax.xaxis.set_minor_formatter(h_fmt)\n", "ax.xaxis.set_major_locator(days)\n", "ax.xaxis.set_major_formatter(d_fmt)\n", "\n", "ax.set_xlabel(\"Time [UTC]\")\n", "ax.set_ylabel(\"Irradiance [Wm$^{-2}}$]\")\n", "\n", "ax.legend()\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Deleting temporary files:" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Successfully deleted all temporary files\n" ] } ], "source": [ "rad.close()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.4" } }, "nbformat": 4, "nbformat_minor": 2 }