{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# ``atmosphere.py`` tutorial\n", "\n", "This tutorial needs your help to make it better!\n", "\n", "This tutorial requires pvlib > 0.6.0.\n", "\n", "Authors:\n", "* Will Holmgren (@wholmgren), University of Arizona. 2015, March 2016, August 2018." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n", "import matplotlib.pyplot as plt\n", "\n", "# built in python modules\n", "import datetime\n", "import logging\n", "import os\n", "import inspect\n", "\n", "# python add-ons\n", "import numpy as np\n", "import pandas as pd" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import pvlib\n", "from pvlib.location import Location" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "tus = Location(32.2, -111, 'US/Arizona', 700, 'Tucson')" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Location: \n", " name: Tucson\n", " latitude: 32.2\n", " longitude: -111\n", " altitude: 700\n", " tz: US/Arizona\n" ] } ], "source": [ "print(tus)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ " apparent_zenith zenith apparent_elevation \\\n", "2014-06-24 00:00:00-07:00 124.039086 124.039086 -34.039086 \n", "2014-06-24 00:01:00-07:00 124.065488 124.065488 -34.065488 \n", "2014-06-24 00:02:00-07:00 124.090882 124.090882 -34.090882 \n", "2014-06-24 00:03:00-07:00 124.115267 124.115267 -34.115267 \n", "2014-06-24 00:04:00-07:00 124.138641 124.138641 -34.138641 \n", "\n", " elevation azimuth equation_of_time \n", "2014-06-24 00:00:00-07:00 -34.039086 352.695507 -2.366772 \n", "2014-06-24 00:01:00-07:00 -34.065488 352.970528 -2.366922 \n", "2014-06-24 00:02:00-07:00 -34.090882 353.245771 -2.367072 \n", "2014-06-24 00:03:00-07:00 -34.115267 353.521228 -2.367222 \n", "2014-06-24 00:04:00-07:00 -34.138641 353.796891 -2.367372 \n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAD7CAYAAACfQGjDAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XlclWX++P/XxWEVlF12BRFlFwUVtRQ0l3Lvq2NONVqT1uTU9JtpsfnMZM04leWnmfHTNk4plY65lOWUlUsuuQuKC4sCArIo+74ezrl+f4CECorK4QBez8eDB3Cfe3mfw+F+n/u6rvt9CSkliqIoitIaE2MHoCiKonRdKkkoiqIobVJJQlEURWmTShKKoihKm1SSUBRFUdqkkoSiKIrSJpUkFEVRlDapJKEoiqK0SSUJRVEUpU2mxg7gTjk5OUlvb29jh6EoitKtxMXFFUopnW+2XrdPEt7e3sTGxho7DEVRlG5FCJHZnvVUc5OiKIrSJpUkFEVRlDapJKEoiqK0qdv3SbRGq9WSnZ1NbW2tsUNRugBLS0s8PT0xMzMzdiiK0u30yCSRnZ1N79698fb2Rghh7HAUI5JSUlRURHZ2Nj4+PsYOR1G6nR7Z3FRbW4ujo6NKEApCCBwdHdVVpaLcph55JQHcVoJoqNdRWVqHqbkJNnaWBohKMQb1YUFRbl+PTRK3Qq+XVJXWUVNRD4DUa4wckaIoStdg0OYmIYSlEOKYEOKUECJBCPFa0/IYIUS6ECK+6SusabkQQqwSQqQKIU4LIYYZMj4pJTWV9RTnVFJTUY+ljRmm5ipBKIqiXGHoPok6YLyUcggQBkwRQkQ2PfaClDKs6Su+adn9gF/T12LgA0MF1lCvo+RyNRVFtZiYmmDvak0fRytMTARSGuqoXUtMTAy5ubkGP05ubi5z5swBID4+nu3btzc/9uqrr7Jy5UqDx6AoPUl9TQOHvkgl8aDh/38N2twkpZRAZdOvZk1fNzoFzwQ+bdruiBDCTgjhJqW8dLsxvPbfBBJzy1sJDrT1OjSmJphofm6zbtDqQcobXlEEuvdh2fSg2w2pwzU0NGBqeut/ypiYGIKDg3F3dzdAVD9zd3dny5YtQGOSiI2N5YEHHjDoMRWlJ5J6ybmjlzm8NY3q8nrC7vMy+DENPrpJCKERQsQD+cBOKeXRpof+1tSk9HchhEXTMg8gq8Xm2U3Lrt3nYiFErBAitqCg4DYDAzMLzVUJoqPNmjWL8PBwgoKCWL16NQA2Njb84Q9/YNiwYUyYMIEr8UdFRfHcc88xevRogoODOXbsGADHjh1j9OjRDB06lNGjR3Pu3Dmg8QQ/d+5cpk+fzqRJkwB4++23GT58OKGhoSxbtgyAjIwMAgICWLRoEUFBQUyaNImamhq2bNlCbGwsDz/8MGFhYdTU1FwXf2xsLGFhYYSFhRESEtLcAZyWlsaUKVMIDw/n3nvvJTk5GYCFCxfy7LPPMnr0aAYMGNCcGDIyMggODqa+vp5XXnmFjRs3EhYWxsaNGwFITEwkKiqKAQMGsGrVKoP8LRSlu7uUVsaWt+LY/UkSvR0tmbM0gjFz/Ax/YCllp3wBdsAeIBhwAwRgAXwCvNK0zrfAPS222Q2E32i/4eHh8lqJiYnXLWuvkstVsii38ra3b6moqEhKKWV1dbUMCgqShYWFEpDr1q2TUkr52muvySVLlkgppRw3bpx84oknpJRS7tu3TwYFBUkppSwrK5NarVZKKeXOnTvlgw8+KKWUcu3atdLDw6P5GD/88INctGiR1Ov1UqfTyalTp8p9+/bJ9PR0qdFo5MmTJ6WUUs6dO1d+9tlnzcc8fvx4u57L888/L59//nkppZTjx4+X58+fl1JKeeTIERkdHS2llHLBggVyzpw5UqfTyYSEBOnr6yullDI9Pb35+axdu7b5OUsp5bJly+SoUaNkbW2tLCgokA4ODrK+vv4WXuX2uZP3hKIYU1lBtfx+9Rn57pO75ZoXf5JJh3OlXqe/4/0CsbId5+5OG90kpSwVQuwFpkgprzRC1wkh1gLPN/2eDbS8fvIEDN/oZiCrVq1i69atAGRlZZGSkoKJiQnz5s0D4JFHHuHBBx9sXn/+/PkAjB07lvLyckpLS6moqGDBggWkpKQghECr1TavP3HiRBwcHADYsWMHO3bsYOjQoQBUVlaSkpJCv3798PHxISwsDIDw8HAyMjJu6Xls2rSJEydOsGPHDiorKzl06BBz585tfryurq7551mzZmFiYkJgYCB5eXnt2v/UqVOxsLDAwsKCvn37kpeXh6en5y3FqCg9TV1NA3HfZXDqxyxMhGD4VG+GTuqPmUXnDq4xaJIQQjgD2qYEYQXcB6y40s8gGtsvZgFnmzbZBvxWCPE5MBIok3fQH2FMe/fuZdeuXRw+fJhevXoRFRXV6g1dLcfwXzueXwjBn//8Z6Kjo9m6dSsZGRlERUU1P25tbd38s5SSl19+mSeffPKqfWRkZGBhYdH8u0ajabVpqS0JCQksW7aM/fv3o9Fo0Ov12NnZER8f3+r6LY8l2zkC4Nr4Ghoa2h2fovQ0ep2exAO5HP1vOrWVWvwjXRk50xcbe4ubb2wAhu6TcAP2CCFOA8dp7JP4BlgvhDgDnAGcgOVN628HLgCpwL+Bpw0cn8GUlZVhb29Pr169SE5O5siRIwDo9frmtvr//Oc/3HPPPc3bXGmjP3DgALa2ttja2lJWVoaHR2O3TExMTJvHmzx5MmvWrKGysnGcQE5ODvn5+TeMsXfv3lRUVNzwOTz00EN8+umnODs3zk3Sp08ffHx82Lx5M9CYCE6dOnXD49zKMRXlbiWl5EJ8AZ8vP86+DedxcLNm7ssRTFgYaLQEAYYf3XQaGNrK8vFtrC+BJYaMqbNMmTKFDz/8kNDQUAYPHkxkZOPIX2traxISEggPD8fW1rY5MQDY29szevRoysvLWbNmDQAvvvgiCxYs4J133mH8+FZfNgAmTZpEUlISo0aNAho7yNetW4dG0/al6cKFC3nqqaewsrLi8OHDWFlZXfX4V199RWZmJosWLWpeFh8fz/r16/nNb37D8uXL0Wq1PPTQQwwZMqRdr0t0dDRvvvkmYWFhvPzyy+3aRlF6uuzkYo58fYG89HLsXHpx/1Mh+Axx6hLVAkR7mwS6qoiICHntzHRJSUkEBATc1v5K86rR6yUObtY3X/k22NjYNH/abykqKoqVK1cSERFhkOPe7e7kPaEohpKXXs6Rr9PITi7Bxt6C4dN88I90xURj+LJ6Qog4KeVNTziqLIeiKEonK8qp5Oi2C6SfKsTSxox75voRNNYdU7OuV/FBJYlO1tpVBDR2dBvT7NmzSU9Pv2rZihUrmDx5spEiUpSeJz+znNjtGaSfKsTcUsOI6T4MmeCFuWXXPRV33ciUTnVlqK6iKB3vUloZsdvTuZhQjEUvU4ZP9SZ0vBeW1l1/IiyVJBRFUQxASkl2UglxP2SQc64USxszImcNIGScJ+ZW3efU230iVRRF6QYatDpSjudxancWRTlV9LI15565fgTe497pN8J1BJUkFEVROkBNRT1n9+dwZm82NRVaHD1sGP+rAAYNd0Fj1n0nAVVJ4lrGH5bcbNu2bSQmJrJ06VK++uorBg0aRGBgIKCGzCpKV5GfWc7ZfTmcP5aHrkFP/2BHhtznhedg+y5xn8OdUkmiC5sxYwYzZswAGm9smzZtWnOSUBTFeOprG0g5nkfCT7kUXKzA1NyEwaNcGTLey2D3WBlLz08S3y2Fy2favbqNVtc46dCNZqhzDYH737zhfj788EM+/PBDoLG8hbe3Ny+//DLLli2jrq4OX19f1q5di42NDd7e3ixYsID//ve/aLVaNm/ejL+/PzExMcTGxvLLX/6Sbdu2sW/fPpYvX84XX3wBwObNm3n66acpLS3l448/5t57723381QU5dYVZleQsD+Xc8cuo63V4ehhzdiHBjFopCsW3agz+lb0zGfVBTz11FM89dRTaLVaxo8fz+OPP87y5cvZtWsX1tbWrFixgnfeeYdXXnkFACcnJ06cOMH777/PypUr+eijj5r3NXr0aGbMmMG0adOaZ3iDxsmGjh07xvbt23nttdfYtWtXpz9PRenpqsrqSDmeR/KRyxRlV6IxNWFgRF+C7vXAdUCfHtGkdCM9P0nc5BP/tSrzq9HrOq4sx+9+9zvGjx+Pvb09iYmJjBkzBoD6+vrmOktAc8nw8PBwvvzyy3btu+U2t1r+W1GUtjXU60g/VUjykctkJRYhJfT17sPYhwbhN9ylW9zf0FF6fpIwopiYGDIzM3n33Xf59ttvmThxIhs2bGh13Svlsm+lVPbtbKMoSut0Wj1ZScWknsgnPb6A+lodNvYWDJvcn8GRrti79qy+hvZSScJA4uLiWLlyJT/99BMmJiZERkayZMkSUlNTGThwINXV1WRnZzNo0KB27U+V2FaUjqfT6rmYVExaXD7ppxoTg0UvUwYMdWZwpBsefnYIk57dnHQzhp50yBLYT+M0pabAFinlMiGED/A54ACcAB6VUtY3zXX9KRAOFAHzpJQZhozRUN59912Ki4uJjo4GICIigpiYGObPn988k9vy5cvbnSQeeughFi1axKpVq5rno1AU5dbV1zaQnVRCWnw+GacKf04Mw/riO9QZrwAHNKbd976GjmbQUuFNM89ZSykrhRBmwAHgd8DvgS+llJ8LIT4ETkkpPxBCPA2ESimfEkI8BMyWUs670TE6vFR4B/dJKF2DKhV+d6soriXjdCEZZwrJPleCvkE2JoYwZ3zD++I52P6uSwxdolR40yRCV8qemjV9SWA88Mum5Z8ArwIfADObfgbYArwrhBCyu096oShKp5J6SV5GORlnCsk4XURRTuNpyNbZipBxnniHOuE20BZNJ8zb0N0ZvE9CCKEB4oCBwHtAGlAqpbzS05oNeDT97AFkAUgpG4QQZYAjUHjNPhcDiwH69etn6KegKEo3UFFcS1ZSMdlJxWQll1BbqUUIcBtox+gHB+Id6njXdj7fCYMnCSmlDggTQtgBW4HWrvmvXCm01kN03VWElHI1sBoam5s6KFRFUbqR+toGclNKyUosJiupmJLL1QBY9TGnf5AjXoEO9A92vKuGqxpCp41uklKWCiH2ApGAnRDCtOlqwhPIbVotG/ACsoUQpoAtUNxZMf4cbKcfUVGUm9Bp9eRllpN7vpSspGIuXyhDr5NozEzw8LMjYIw7/QIdcHC37vE3uHUmQ49ucga0TQnCCrgPWAHsAebQOMJpAfB10ybbmn4/3PT4j6o/QlHuTto6HZfTy8g9X0puSil56eXoGvQAOHnZMGSCF16BDrj52nbJaT97CkNfSbgBnzT1S5gAm6SU3wghEoHPhRDLgZPAx03rfwx8JoRIpfEK4iEDx3cd9flDUYyjrqaBS6mlXEotJed8KQWZFej1EiHAyas3weM8cPezw32gHZY2qgmpsxh6dNNpYGgryy8AI1pZXgvMNWRMiqJ0DTWV9VxKKSM3pZTc1FIKsyqQEkw0gr79+xA2sR/ug+xwG2DbrWZy62nUK38Xi4mJYdKkSbi7u9/ytgsXLryu4GBHx/PEE0/w+9//XpVH7yGqSusaE0JTUijOrQJAY2aC64A+RDzgjbufHS4DbDG7URVmpVP1+CSx4tgKkouT271+g1YPUmJ6gzepv4M/L414qSPC6xANDQ2Ymt76nzImJobg4ODbShKGcG08LSvhKt1PeWENuamlzX0KZQU1AJhZaHAbaMugES64D7Sjb/8+3Xrmtp5O/WUMaNasWYSHhxMUFMTq1asBsLGx4Q9/+APDhg1jwoQJFBQUAI0zzT333HOMHj2a4OBgjh07BsCxY8cYPXo0Q4cOZfTo0Zw7dw5oPKHOnTuX6dOnM2nSJADefvtthg8fTmhoKMuWLQMgIyODgIAAFi1aRFBQEJMmTaKmpoYtW7YQGxvLww8/TFhYGDU1Na0+h7i4OMaNG0d4eDiTJ0/m0qVL7VonKSmJESN+blHMyMggNDQUgL/85S8MHz6c4OBgFi9ejJSy1XiioqK4cjf9hg0bCAkJITg4mJde+jlB29jY8D//8z8MGTKEyMhI8vLybv8Pptw2KSWledUk/JTDzrUJfPLHg3z2p8PsjkniQnwB9m7WjJkzkLkvR/DEO/cy/Zkwwqd44zbQTiWIrk5K2a2/wsPD5bUSExOvW9ZepXlVsiin4ra3b6moqEhKKWV1dbUMCgqShYWFEpDr1q2TUkr52muvySVLlkgppRw3bpx84oknpJRS7tu3TwYFBUkppSwrK5NarVZKKeXOnTvlgw8+KKWUcu3atdLDw6P5GD/88INctGiR1Ov1UqfTyalTp8p9+/bJ9PR0qdFo5MmTJ6WUUs6dO1d+9tlnzcc8fvx4m/HX19fLUaNGyfz8fCmllJ9//rl87LHHpJRSLliwQG7evPmG6wwZMkSmpaVJKaV888035V//+terXhcppXzkkUfktm3bWo3nyu85OTnSy8tL5ufnS61WK6Ojo+XWrVullFICzdu/8MILzce41p28J5Tr6fV6WZRTKc/szZLf//uMXPPiT/LdJ3fLd5/cLT9+fr/87l9n5Kkfs2RhdoXU6/TGDldpBRAr23GO7fHNTca0atUqtm7dCkBWVhYpKSmYmJgwb15jOapHHnmkeU4IgPnz5wMwduxYysvLKS0tpaKiggULFpCSkoIQAq1W27z+xIkTcXBwAGDHjh3s2LGDoUMbxwlUVlaSkpJCv3798PHxISwsDLi1uSfOnTvH2bNnmThxIgA6nQ43N7d2r/OLX/yCTZs2sXTpUjZu3MjGjRsB2LNnD2+99RbV1dUUFxcTFBTE9OnT24zj+PHjREVF4ezsDMDDDz/M/v37mTVrFubm5kybNq35ue3cubNdz025NVIvKcqtIjelpLH5KLWUmorG96K1rTkeg+xx97PDY5Addi691H0KPYhKEq3oiBsz9u7dy65duzh8+DC9evUiKiqK2tra69Zr+c907T+WEII///nPREdHs3XrVjIyMoiKimp+3Nr65xIDUkpefvllnnzyyav2kZGR0TzvBDTOPdFW09K1pJQEBQVx+PDh21pn3rx5zJ07lwcffBAhBH5+ftTW1vL0008TGxuLl5cXr776aquvy7XHaIuZmVnz66bm1ehYlSWNZS6ykkrITi5uTgq9HSzpF+TYnBT6OFmppNCDqcZAAykrK8Pe3p5evXqRnJzMkSNHANDr9c2lvv/zn/9wzz33NG9z5ZP2gQMHsLW1xdbWlrKyMjw8GktbxcTEtHm8yZMns2bNGiorGwuZ5eTkkJ+ff8MYbzZHxeDBgykoKGhOAFqtloSEhHav4+vri0aj4a9//Wvz1dOVhODk5ERlZeVVZc/bimfkyJHs27ePwsJCdDodGzZsYNy4cTd8bsqt09bryDhTyE8bz/OfV4/wycuH+PHTZLLPleAV6MCEBQE8+rdR/Or10dy3MJDAMe7YOqurhp5OXUlcq4Pe71OmTOHDDz8kNDSUwYMHExkZCTR++k9ISCA8PBxbW9vmxABgb2/P6NGjKS8vZ82aNQC8+OKLLFiwgHfeeYfx48e3ebxJkyaRlJTUPCWqjY0N69atQ6Npe5TWwoULeeqpp7CysuLw4cNYWVld9bi5uTlbtmzh2WefpaysjIaGBp577jmCgoLavc68efN44YUXSE9PB8DOzo5FixYREhKCt7c3w4cPbzOeK9zc3HjjjTeIjo5GSskDDzzAzJkzb/wHUNqlqqyOzDNFpJ8uJCupGJ1Wj6mZCe5NZS68Ahxw9FBlLu5mBp1PojN09HwSZQXVNGj1OLrbdER417GxsWn+tN9SVFQUK1euJCLipuXdldug5pP4WXlhDSmxeaSfKiQvvRxobELyGeKEd4gTbn6qzMXdoEvMJ6EoStdQU1FPalw+KcfzuJRWBkDf/r0ZOcMH71BndbWgtEkliU7W2lUENHZ0G9Ps2bObm4SuWLFiBZMnTzZSRMqdknpJdnIJZ3/KIf1UIVIvcXC3JnLWAPwiXOjjZHXznSh3PZUkFIDmobpK91df20DigVzO7suhrKAGSxszhkzwwj/SFUcPwzSjKj2XShKK0kPUVmk5szebUz9mUVfVgJuvLcOn+TBwWF91V7Ny2ww9n4QX8CngCuiB1VLKfwohXgUWAQVNq/5RSrm9aZuXgV8DOuBZKeUPhoxRUbq7Bq2OU7uziPs+E22tDu9QJ8Lv74+rj62xQ1N6AENfSTQAf5BSnhBC9AbihBBXbon9u5RyZcuVhRCBNM4hEQS4A7uEEINk4xSoiqJcI+1kPge3pFJRVIt3qBMjZwzAyVM1KSkdx6DXoFLKS1LKE00/VwBJgMcNNpkJfC6lrJNSpgOptDLvRHfm7e1NYWFhh+7z9ddfv+r30aNHd+j+la6npqKe71ef5ft/ncXc0pQZz4Ux9elQlSCUDtdpDZVCCG8aJyA62rTot0KI00KINUII+6ZlHkBWi82yaSWpCCEWCyFihRCxV6qodmCkHbw/w7s2SRw6dMhIkSidISupmA1/OUr66QIiZw3gF3+MwMvfwdhhKT1Up3RcCyFsgC+A56SU5UKID4C/0lgm6a/A/wKP0/oZ+rq7/aSUq4HV0Hgz3Y2Offn116lLurX5JKSUVN5gPgmLAH9c//jHm+5r3bp1rFq1ivr6ekaOHMn7779/08dXr15Neno6b731FtBYiiMuLo7/+7//Y9asWWRlZVFbW8vvfvc7Fi9ezNKlS6mpqSEsLIygoCDWr1/ffMOelJIXX3yR7777DiEEf/rTn5g3bx579+7l1VdfxcnJibNnzxIeHs66devUOPkuTkrJyR0XOfJVGvZu1sx8bqgaraQYnMGvJIQQZjQmiPVSyi8BpJR5UkqdlFIP/Jufm5SyAa8Wm3sCuYaO0RCSkpLYuHEjBw8eJD4+Ho1Gw/r162/6+Jw5c/jyyy+b19u4cWNz3aM1a9YQFxdHbGwsq1atoqioiDfffBMrKyvi4+Ov2j/Al19+SXx8PKdOnWLXrl288MILzfNBnDx5kn/84x8kJiZy4cIFDh482AmvinK79HrJ3vXnOLw1Dd9hffl/L4arBKF0CkOPbhLAx0CSlPKdFsvdpJRXZq+ZDZxt+nkb8B8hxDs0dlz7AcfuJIb2fOJvqayghgat7o7LcuzevZu4uLjm2kQ1NTX07dv3po87OzszYMAAjhw5gp+fH+fOnWPMmDFA66XHHR0d24zhwIEDzJ8/H41Gg4uLC+PGjeP48eP06dOHESNG4OnpCUBYWBgZGRlXFRtUug6dTs+utYmkxuYTPqU/I2cOUFd9SqcxdHPTGOBR4IwQIr5p2R+B+UKIMBqbkjKAJwGklAlCiE1AIo0jo5YYZWRTB5SzklKyYMEC3njjjauWX6nk2tbj0FgUb9OmTfj7+zN79myEEO0uPX5tDG25tny4KrHdNUkp2ftZMqmx+Yx60Jdhk/obOyTlLmPo0U0HpJRCShkqpQxr+toupXxUShnStHxGi6sKpJR/k1L6SikHSym/M2R8hjRhwgS2bNnSXK67uLiYzMzMdj3+4IMP8tVXX7Fhw4bmpqa2So9D45wKLScjumLs2LFs3LgRnU5HQUEB+/fvv2pKUaXrO/L1BZKPXGbEdB+VIBSjULdhGkhgYCDLly9n0qRJhIaGMnHixKvmh77R4/b29gQGBpKZmdl8Up8yZQoNDQ2Ehoby5z//ubn0OMDixYsJDQ3l4YcfviqG2bNnExoaypAhQxg/fjxvvfUWrq6unfDslY6QGpfPie8zCbzXnYgHvI0djnKXUqXCr1FWUENDvU51CvYw3a1UeGleNZveOI6DmzWz/zAMjan6PKd0rPaWClfvPEXpYvR6ya6YREw0gsmLglWCUIxKvfsUpYs5uy+HvPRy7v3FIHo7WBo7HOUup5KEonQhlSW1HPkqjX6BDgwa4WLscBRFJQlF6UqO/jcdnU7P2PmD1b0QSpegkoSidBHFuVWcO3yJkHGe2DqrWeOUrkElCUXpIo5uu4CZhYbw+9X9EErXoZJEF/DEE0+QmJjYIftqWRE2IyOD4ODgDtmvYlgll6u4EF9A6HgvrGzMjR2OojRTSaIL+OijjwgMDOyQfV1bNlzpHuJ3XkRjZkJotKexQ1GUq/T4Oa5/2nSewqzKdq9/pVS42Q1KhTt52XDvLwbddF/XlvZ2dXXllVdeARoL+tXX15Oenk5UVBQrV64kIiICGxsblixZwq5du7C3t+f111/nxRdf5OLFi/zjH/9gxowZxMTEEBsby7vvvgvAtGnTeP755/n++++vKhv+t7/9DZ1Ox6JFizh06BAeHh58/fXXWFmp9u6upKqsjuSjlwkY7Y5Vb3UVoXQt6krCgK4t7T1mzBji4+OJj49nyJAhPP/889dtU1VVRVRUFHFxcfTu3Zs//elP7Ny5k61btzYnmLa0VjY8JSWFJUuWkJCQgJ2dHV988YVBnqty+5IPX0LfIAmb4HXzlRWlk/X4K4n2fOJvqbywBm1dx5TlaKu091tvvYWVlRVLliy5bhtzc3OmTJkCQEhICBYWFpiZmRESEkJGRsYtx+Dj40NYWBgA4eHht7UPxXCkXpJ48BIeg+ywc+ll7HAU5To9PkkYS1ulvXfv3s3mzZvZv39/q9uZmZk1j483MTFpLultYmLSXM7b1NQUvV7fvM2NSoZfWxK8pqbmjp+b0nFyUkopL6hhxDQfY4eiKK0yaHOTEMJLCLFHCJEkhEgQQvyuabmDEGKnECKl6bt903IhhFglhEhtmv96mCHjM6TWSntnZmby9NNPs2nTpjvqF/D29iY+Ph69Xk9WVhbHjv08L1NbZcOVrinxQC4WvUzxHeps7FAUpVWG7pNoAP4gpQwAIoElQohAYCmwW0rpB+xu+h3gfhpno/MDFgMfGDg+g2mttHdGRgZFRUXMnj2bsLAwHnjggdva95gxY/Dx8SEkJITnn3+eYcN+zqVtlQ1Xup762gYuxBfgN9wF0xsMlFAUY+rUUuFCiK+Bd5u+oqSUl4QQbsBeKeVgIcS/mn7e0LT+uSvrtbXPji7xUx5OAAAgAElEQVQV3pF9EkrX0RVLhafE5rHjowRm/2EY7n52xg5Huct0uVLhQghvYChwFHC5cuJv+n5l8mcPIKvFZtlNy67d12IhRKwQIragoMCQYSuKwaSdKMCqjzmuvrbGDkVR2tQpSUIIYQN8ATwnpSy/0aqtLLvuUkdKuVpKGSGljHB2Vm25SvfTUK8jM6GIAWHOmJioQn5K12XwJCGEMKMxQayXUn7ZtDivqZmJpu/5TcuzgZaDxT2BXEPHeK3uPVef0h1cTCymoU6nOqyVLs/Qo5sE8DGQJKV8p8VD24AFTT8vAL5usfxXTaOcIoGyG/VHKEp3lXYyH0trM9wHqb4IpWsz9H0SY4BHgTNCiPimZX8E3gQ2CSF+DVwE5jY9th14AEgFqoHHDByfonQ6nVZPxqlCfIf1RaNRRQ+Urs2gSUJKeYDW+xkAJrSyvgSuvw25s6n2JsWAspKLqa/VMUA1NSndgPoY002Vlpby/vvvN/+em5vLnDlzOuXYmzdvJiAggOjo6FYfj4+PZ/v27c2/b9u2jTfffLNTYusOLpwswNxSg5e/g7FDUZSbUkmim7o2Sbi7u7Nly5ZOOfbHH3/M+++/z549e1p9/NokMWPGDJYuXdrquncbvU5P+qlCvEOd0Jipfz+l6+vxtZv2xKwmP/NCu9fXafXo9RIzi7bvgO3bfwDRCxffdF/r1q1j1apV1NfXM3LkSN5//30+/fRT3njjDdzc3Bg0aBAWFha8++67LFy4kGnTpjVfDdjY2FBZWUllZSUzZ86kpKQErVbL8uXLmTlzJkuXLiUtLY2wsDAmTpzIkiVLmDZtGmfPnqW2tpbf/OY3xMbGYmpqyjvvvEN0dDQxMTFs27aN6upq0tLSmD17Nm+99Vab8W/YsIHXX38dKSVTp05lxYoV/OUvf+HAgQOkp6czY8YM3n777au2qa+v55VXXqGmpoYDBw7w8ssvU1NT01zafOHChVhZWZGcnExmZiZr167lk08+4fDhw4wcOZKYmBgAduzYwbJly6irq8PX15e1a9diY9P9b3DMSSmltkqL77C+N19ZUbqAHp8kjCUpKYmNGzdy8OBBzMzMePrpp1m3bh3Lli0jLi4OW1tboqOjGTp06A33Y2lpydatW+nTpw+FhYVERkYyY8YM3nzzTc6ePUt8fON4gJbVXd977z0Azpw5Q3JyMpMmTeL8+fNA46f8kydPYmFhweDBg3nmmWfw8rq+RHVubi4vvfQScXFx2NvbM2nSJL766iteeeUVfvzxx+b5L65lbm7OX/7yl6vmu7hy4r+ipKSEH3/8kW3btjF9+nQOHjzIRx99xPDhw4mPj8fT05Ply5eza9curK2tWbFiBe+8885NS6V3BxdOFGBqoaFfoGpqUrqHHp8k2vOJv6Xywhrqa3U4ed7Zp9bdu3cTFxfH8OHDgcZJhg4dOkRUVBRXbgCcN29e88m7LVJK/vjHP7J//35MTEzIyckhLy/vhtscOHCAZ555BgB/f3/69+/ffJwJEyZga9t4h29gYCCZmZmtJonjx49fFevDDz/M/v37mTVr1i28Cq2bPn06QghCQkJwcXEhJCQEgKCgIDIyMsjOziYxMZExY8YAjVcno0aNuuPjGpteL0mLL6B/kKOq1aR0Gz0+SRiLlJIFCxbwxhtvNC/76quvmueXuFbL8t9SSurr6wFYv349BQUFxMXFYWZmhre39w1Lg1/Zvi3Xlg6/Un78VvZxp1qWP28Zz5Vy6BqNhokTJ7JhwwaDxWAMl9PKqCmvx3eYGtWkdB+q58xAJkyYwJYtW8jPb7yZvLi4mKFDh7J3716KiorQarVs3ry5eX1vb2/i4uIA+Prrr5vLfZeVldG3b1/MzMzYs2cPmZmZAPTu3ZuKiopWjz127NjmmenOnz/PxYsXGTx48C3FP3LkSPbt20dhYSE6nY4NGzYwbty4dm17o9jaIzIykoMHD5KamgpAdXX1Ta+4uoO0k/loTE3oH+xo7FAUpd1UkjCQwMBAli9fzqRJkwgNDWXixIlcunSJV199lVGjRnHfffddVeJ70aJF7Nu3jxEjRnD06FGsra2Bxmae2NhYIiIiWL9+Pf7+/gA4OjoyZswYgoODeeGFF6469tNPP41OpyMkJIR58+YRExNz1Sf29nBzc+ONN94gOjqaIUOGMGzYMGbOnNmubaOjo0lMTCQsLIyNGzfe0nEBnJ2diYmJYf78+YSGhhIZGUlycvIt76crkXpJ2okC+gU5YG6pLuCV7qNTS4UbgiFKhXdEn0R7xMTEXNXBqxiOsUuFX0or48u347jvsUAGj3Q1WhyKckWXKxWuKHeztLjGpiafUCdjh6Iot0Rd9xrRwoULWbhwobHDYOTIkdTV1V217LPPPmseddSWH374gZdeeumqZT4+Pm12zt+tpF6SdjIfr0AHzK3Uv5zSvah3rMLRo0dva7vJkyczefLkDo6m58nLKKeypI7IWb7GDkVRblmPbW7q7n0tSscx9nvh/PE8TEwF3qqpSemGDD2fxBohRL4Q4myLZa8KIXKEEPFNXw+0eOxlIUSqEOKcEOK2P6JaWlpSVFR0eycHNUlYjyKlpKioCEtLS6McX9egJ+VYHgOGOGOhmpqUbsjQ79oY4F3g02uW/11KubLlAiFEIPAQEAS4A7uEEIOklLpbPainpyfZ2dnczvzXtZVaGrQ6CiqMc1JROp6lpSWenp5GOXbGmUJqq7T4j3IzyvEV5U4Zej6J/UII73auPhP4XEpZB6QLIVKBEcDhWz2umZkZPj4+t7oZAD9+lsTFs+UsXHHjmkqK0h7Jhy9jbWuOl6rVpHRTxuqT+K0Q4nRTc5R90zIPIKvFOtlNyzqVam1SOkpVWR2ZZ4sYHOmKiYl6ZyndkzGSxAeALxAGXAL+t2l5a/9FrXYqCCEWCyFihRCxt9OkdDOqy1vpCAk/5SL1koDR7sYORVFuW6cnCSllnpRSJ6XUA/+msUkJGq8cWpYj9QRy29jHaillhJQy4kqVUkXpSnRaPWf3ZdM/xBE7l17GDkdRblunJwkhRMsevNnAlZFP24CHhBAWQggfwA841tnxIVSzgHLnUmLzqKnQMmT89WXYFaU7MWjHtRBiAxAFOAkhsoFlQJQQIozGVp0M4EkAKWWCEGITkAg0AEtuZ2RTh1DtTcod0OslJ37IxMHdGk9/+5tvoChdmKFHN81vZfHHN1j/b8DfDBeRohheyvE8Si5XM3lRMEJdmSrdXI+94/q2qf9p5Q7odXqOf5uOo4cNvkNVf5nS/akk0QrV2qTcrjP7cijLr2HkDB+EGvaq9AAqSShKB6kqq+PYtgv0C3RQdZqUHkMliWuoz37K7ZBScmBzCg1aPffOG6T6IpQeQyWJ1qgKssotOn8sj9TYfIZP9Vb3RSg9yl1blrIwN5PUfetBW4NsqIWGWkRDLflZ4WhrvDjywVMgTJAmpghzG4SVLZpedphZ22PR2wHbvl44uXljZn5rc0crPU9pXjX7N5zDzdeWYVO8jR2O0gNJvZ6ykiKKLmdQW5xLXUUh2soirJwHEBo9x6DHvouTRCqR51Y0/14vTanDjPLavlToXQi5vBUNOkzRYSr0re5DLwWFwpYSU2cqLFypsx2Aad9B9PEMxNU3FFt71S7d09VVa/n2/dOYmJpw32OBqkaTctv0DQ1cvniegswEqnOTEUUp9K7KxLY+Hwd9EXaiDrtrtontPR5UkjAM39B7KPU+h4WVNeYWVpibmmIOmP/nHOJkPtav5TWvW1tTRWVZMdXlRdSUF1NbXkhdSQ660hw0lblY1uThVHMBt8pDmOXqIL5xuwLsybXyo8YxBMv+EbgHjsLZrT/CRLXy9QQNWh3frz5LeUENM/+/MPo4WRk7JKWb0NbVcPFcHEXnj6HPjce+PIn+2nTchZYrlb7KseaSqSeXrQeRbe2K6OOOub0H5nbu9LJzobe9M6H2hh9mfdcmCTNzC+ycXNu1rqWVNZZW1uB64xIL2vo6sjKTKcpIoPZSEpqiczhVJBOcdRxN9sdwEAqxI8s6mDqPSBwDo/AJGompmXlHPCWlEzVodXz34Rmyz5UwYUEA7n7qzmqlbeVF+aSd2EVN6gEci+Lw0abgK3T4AhVYkWXhx0mHBxF9A7H1CsBtQAh2Tm706QIDIO7aJNGWO/mTmJlb4OU3BC+/IVctr64s42LiMUrTjqO5fBL38lN4nD8A51dSudWKC1bBVLmNwGnI/fiGjMZEo7mzJ6EYVF1NAz+sPkNWUgnRj/rjH6kmFFKuVlZazIVj22lI2Y1zURze+kyGAvVSwwXzwcS6PoRZv2G4+Y/Cw8efQJOu+z+vkkQrOnpwUy8bW/xHTIQRE5uX5eekczF+F7oLB3ApOUFo+nuQ/h7FX/XhQp8RMPA+fEZOx9HFODOqKa0rL6zh2/dPU3q5mvG/CiBgtEoQCki9jgtnD1NwYjt9cvbhV5/EUKGjSlqQZhnMYdf7sfMfx4Cwsfhb2Rg73FuikoSR9PXwoa/HImARAEV52aQf/S+k7mZA+VEcTuyCE0tJ1fhS4DEBlxFz8AkcrvozjOjCyQJ+XJcEEqY/OwRPfzXb3N2suqqCpIPb0Cd+g2/pQXwpwxdI0wzghOcj9A6ajG/4eEItundflUoS1zJSE6CjiyeOM34D/Aa9TkfKmUMUxX+HXc4eRmb+G5OLq8kRLmT1HY/tsNkMCp+AxlT9+TpDXbWWQ1+kknjwEs79ejPp10HqXoi7VFlJASk/bcYkeTv+VccIF3WU04tzvSNJ8b2PAZHT8XXth6+xA+1A6izTGiPfS2ei0eAXdi9+YfcCUHg5iwsHNmOR9h3DLm/G/LsNFH1nS5rDWHoNnUvAqKkqYRiAlJJzRy9z6ItUaiu1DJvcnxHTfdCYqqu5u0lxwSXS9q2nV8o3DK49RYTQk48DZ52nYj1kJoNHTmG4uaWxwzQYg59ZhBBrgGlAvpQyuGmZA7AR8KZxTolfSClLRGMtg38CDwDVwEIp5QlDx9jVObl64TTn98DvqSgr5syBLyH5W4KKdmK9+78U7rYj1XkidiPnM3hYtGqSukNSSrKSijm6LZ38jHJcfPow/ZkwnPv1NnZoSieprCgl8cfPMUv6kuCaWIYLHVnCjWPuj+AQPptBQ8fR9y4ZYNIZHz9jgHeBT1ssWwrsllK+KYRY2vT7S8D9NM5I5weMpHE+7JGdEGMz4w84u7Hetg6ET30Cpj5BTVUFcfu2YJKwhaH5X2HxzWZyv3Uh020KrmMexieoU1+6bu9KcojdnsGl1DJs7C2IftSfgFFuqqLrXaC+toakn76g4dRmAisOMULUkyecOOE+H8fIX+IbMgqvu/ADmMGThJRyvxDC+5rFM2mcsQ7gE2AvjUliJvCplFICR4QQdkIINynlJUPH2ZI0dntTO1lZ9yb8gcfggccoLy3i9N4NWCRvZXjOZ5hu/oT0L7zJ830Qv/sex9FFTaPZFm2djnNHL3P6xyxKLldjbWvO2IcGETjGHY3Z3XdSuJtIvZ6k2D1UHv2EgKKdDKGaYvpwynkqdsPnM3j4fbh04eGpncFYDdkuV078UspLQoi+Tcs9gKwW62U3LbsqSQghFgOLAfr162f4aLuBPnaODJ/1W+C3FOVlk7p3PXYpXxCZ8g7a8//kpPVI5JBfEhw1F3OLntt+2l5SSvLSyzl35DIpsXnUVTfg3K839y0MYGC4i0oOPVx+Tgapuz7CI2MrgTKbGmnOGdtxmIU9RPC9M4hUN7g262q9na1d01/3sV5KuRpYDRAREdGxH/u7wB2Od8rRxRPHeS8BL5GZFEfu/jUMvPQtzod/S8nhlznhPAXnex/DN3SMsUPtdGUFNaQcz+Pc0cuU5lVjamaCT5gzIeM8cPW1VSW+e7D62moS9nyO5tR/CKqJpa+QJJkFERe4mMETHmVEHzWkuTXGShJ5V5qRhBBuQH7T8mygZbuIJ5Db6dF1j9amdukfEE7/gHAatH/n1IGvaYhbx7D8rZh/uZm0r30oGDiPgMlPYOvQM6falFJSmF3JhfgC0uMLKcqpBMBjkB3DJvfDd2hfzK262mclpSNlnD1Mwb7VDCrYwVAqycORY54L8Yr+NQEDQ4wdXpdnrP+ObcAC4M2m71+3WP5bIcTnNHZYl3V2f0RPZWpmzpDouRA9l7KiPE7uisHh/CYiz71JTfI7HLUbj/ekJbgE3tvtr6Ya6nXkppZy8WwxF04VUFFUCwLcfG0ZM2cgA4Y608exe9/gpNxEfRVnfliLxalPGNRwHldpxpk+YzEPf4Tge2bgooaMt1tnDIHdQGMntZMQIhtYRmNy2CSE+DVwEZjbtPp2Goe/ptI4BPYxQ8d3fcCdfsROZ+vowsim5qi00wco3reaoMIfsNn8PbkWA6gJeRTv8Y+h6dU9itZJKSnOreJiYjFZScXkppSi0+oxMRV4+TsQcb833qFO9Oqj2pl7upL0k+T9+AFe2f8lRFaTLrw4MuhFBk9exHDHvjffgXKdzhjdNL+Nhya0sq4Elhg2IqUl39B78A29h3OZuZz9YQ3+OV8QFPsatbFvkOA0kb5RT+Ea1PWuLiqKa8k9X0L2uRKyEoupKqsHwN61F8H3euAV6ID7IDvMzO/ukSl3A31dNSl7PsP0ZAy+dYn0kmYcNB9DWdCjTLl/Fj4W6qrhTqhXTwFgcH93Bi/+E7Xalzlw+Ed0x9cSXrAbmy3fcvFrH8qDFzB40q8xs+rT6bFJKSnLryE3tZTc86XkppRSUVwLgEUvU7wCHPAKdMArwIHeDmrk1t2iOOMMWTvfY0DONgZTRQbu7PJ6ln4TnmC8d39jh9djqCRxja71ebnzWZppuGfsRBg7kbyCQo7+8DFeaRsIPvkqlSff4pTrdDwmPoObr+E6/KReUny5qjEhpDYmheqmKwWr3ma4D7RjyH1eeAyyw8HdptvMBqfT66itLKW+rppe1vaYW1mr0VS3SN+g5dz+TYjj/8a/5iQ2UsPxXvdA+ONEjJuGt5k6pXU09Yq2oqNLhXdXLs5OuDzyEjrdC8Qe+oGGI/8i/NIWzD7byBnLCMTIxQSN/X8IzZ29jfR6SVF2Jbkppc1ftVVaAKztLPAYZI+7nx3ufnbYu/bqFifWurpqkg5sI2/fThrOpWCbVUKfsgbMdD+vU20hKHG2pK6/C+bDhjA4ejZeASO6xfPrbOVFl0ne/h79LmwgQBZwCUd+9PgNAyY9xZj+3sYOr0cTspufESMiImRsbGyH7e+njedJPnKZRX8f22H77Eku5WRw4fv38MvaQl+KyRUu5Po9zOD7f0Nv+/Z1DOp0egoyK5oTwqXUUuprG8+efZwsmxOCu589fZwsu81Js7KimJNb/03l9z/QN+ESvepADxQ6m1HZzwnh1hdhZ4vGzJyG6mpkYSEmOfk4ZZZhU9P4f5jvbEb5yAB8Z/4S/zHTMbkLy0C0lHn2MAW7VxFcvBNLoeWM2RCqw37NkPsewtLCwtjhdWtCiDgpZcRN11NJ4mo/bTpP8mGVJG6mrq6W+J3rsT61hmDtWWqlGYmOk3G+7xm8AiOvWrehXkdeenlz09HlC2U01OuBxo5mNz87PPzscBto1+36FPQ6HfHffUruF5/jHncRq3ooszGhaGh/7MaNx3/SXBz73rh9XK/Xk376J1J3fwX7j+BxvhSNhEJHM8rviyD00WfwGji0k56R8em09STs+gyLkx8xuD6RamnBKccpOI//LQODRxg7vB5DJYnb9NOm8yQfusSif4zrsH32dMmnjlC85z3CSn6gl6gj0XQo+f2exNwmgstp5eRllKPXSRDg5GmD+8DGKwW3gXbddljq5ZzznFyzkt7bD+NY0kC1BeQN98Fl1lxCJ/8SM7Pb/5Rbkn+RuC0fotu+m36p5eiBLH97es+ZzYi5v8Wsm09i05by/GzOf7cK7/RNOFFCNq5k+P6SoAeext6xZ97saUwqSdwmlSRuXU1lPZdSy0g7nUvWqQxqqqwBEwQN9OpdRf+hvviEuOPma4tFLzNjh3vbpJSc3b2JzE9W0+9ELmY6uOjbG7MHpzFy3jNY23T8fSXZyXGc/uyf2O46gUOZjtLeJhRNDmfYEy/h7h3U4cczhozTByj98R8ElvyIudBx0iICXfgiwsbPwVTd9GYw7U0S6i9wDXHXj2+6ufraBnJTSslOLiE7uaS51IXGzARXHy9cBvShquIEzpn/YkhDLNXnLDhTPBWNxe/wGhRm5OhvXX1lObFr36bhi29wvlyLm4Uge0IQAY8/y+QwwzZLevqH4/m3T9G+Wsfxr/5F+YaN+G45TtGXczgxzA3PBYsInTCv2/VdNGi1nN6zEavYDwioP0ultOKY04O4TvgtQwO733ukJ1NXEtc4sCmFxEO5LFZXEs0atDryLpSTfa6E7ORi8jIqkHqJxtQEt4G2eAy2x8PPjr79+1xXPTX1zGGKd68irGQH5qKBU1YjMR29hMAx07v85EilF9M48cHf6PP9Uaxr9GR5mKOfNYl7frUUG1tHo8WVlXiMM/9+m757ErCulVx2s0TMm07kr17AslfXnhiporyU0998QL+UT/CSl7iEM+kDHyXggd9i72C81/RupJqbbpNKEj8XxctqKnNxKa0MnVaPENDXuw+eg+3x9LfH1dcWU7P23dFccPkiKd/+H4OzNuJIGRkm/SgMfpyQ+xdhYWVj4Gd0ay4e28P5D1bicvQCJhLOh9jj9NhCRk16HNM7HO7bkSrLijj86Qo0X/yA2+V6yq1NKHpgOBFP/Q99PfyMHd5VcrMukP7t3wm+9AW2oorzZv5UhT9FyISHMVVluY1CJYnbdGBzCokH774kUVfTQFZiMRcTishMKGq+ec3B3RpPf3s8/R1w97PD4g4rptbWVHPq+49xOvMxvvp0SujDea+5+E19DgdX480NInU6zm5dQ2HMWlxTS6g2h/SxAwhY/AeCQscbLa720Ov1xH6zhqJP1uKdUEy9BjJHezPoyd/jHzHRqLElnzxI6Z5/MKxsN6boOd37Xqyjfodf+IQuV+rlbqOSxG06sDmFxAO5LP5nz04SUkqKcirJPFvExYTGqwWpl5hbNZa56B/sQL8gR6xtDTMWXer1nDn0LQ0H3yes+jANmHDWfgKO9z1H/+DOm+eivqKM2I9XwKZvsS+up8DOhMJpIxi16E+4ufh2Whwd5UL8fpL/9b94/HQe8wZI97fDYeECRsxY3Gn9FjqdjpM/bsb8+AeE1sdTjQUJLjPpd//vcfEO6JQYlJtTSeI29eQkoddLLqWWNs+tcKX+kZOXDf2CHOkf5IjrgD6YaDq3ryAj5Qw5P/yTsIL/Yi1qSbIIpWHkEoLGzcXEQJPNl2akcPK95djuOI5VnSS9vwVy3jSifvkS1pZdu12/PYovZ3Lsw+XYfXMI20o9l13M0f9iKqMfW4pVL8PU36qsrODUt//CI3kt3jKbfOFI5sBHCJj6LDZ2TgY5pnL7VJK4TQe2pJD4U89JEg1aHdlJJY2J4XQhtZVaNKYmeAXY4xPmTP9gw10t3KqS4kISv/k/fC+sw5VCLpp4cCng14Q+8CRW1h3Tb3Hx8C5SPvhfXI9nIIHkoQ70fezXjJmwAE0PnMu4vqaKQ5+9jX7D17hdqqXCWlB8/wginv4zTu4dc6V0Kfciqd/8g8DczThSzgVTXyqGPUXwxIVoVH9Dl9Xlk4QQIgOoAHRAg5QyQgjhAGwEvIEM4BdSypIb7ccQSSLhp1ye7MZJQqfVk5lQREpsHhlnimio02FuqaF/iBMDwpzpF+SAuWXX6YC9Vn1dHad2xGAf/y8G6tIopg/Jnr/Ab/pzOLt43XwH15A6HWe+/JjitWtxuVBKpSWkjxtI0OLnCQrqvn/nW6HX64n7/hMK1nyMz9ki6jWQM8aXgKdfxOc2h/Emnz5G8e6/E166Ewuh5Yz1aCzHPoPfiPtVf0M30F2SRISUsrDFsreAYinlm0KIpYC9lPKlG+1HJYlGer0k51wJKcfzSDtZQH1NA1a9zRgQ5syAMGc8BtujMe3aQ06vJfV6ko58h/bAKoZUH6FOmnHCfgrOk37PwMBhN92+vrKc2DVvwcb/Yl9UT769CQUzIhnzxJ9wc/bphGfQNaXE7yX5w7fpd+AC5g1wMdCRvr9exJD7H71pv4VOp+fkvq1ojrzP0PpYaqQ5iX2n4nH/73EdENpJz0DpCN01SZwDolrMfb1XSjn4Rvvp6CRxcEsKZ7tJkpBSUnCxguQjl0mNy6emvB4zSw0DwpwZNNwFT3/7Tu9fMJSclHguff8OIYXbsRBaTliORIx6hiH3TL3uOZblZBD73l+w/a7x/ob0fhbIh6YT9csXe0R/Q0fJz03l2Ad/xWn7cWyrGu+3MJ0/m5G/eh5zy15XrVtVVcXJ7R/hlvgxvjKTIuxI930Y/6nPYuPgaqRnoNyJ7pAk0oESQAL/klKuFkKUSintWqxTIqW8rtaBEGIxsBigX79+4ZmZmR0W18EtKZzdn8OTq6I6bJ8draaynvNH80g6lEtRThUaUxP6hzgyaLgL/YMdMe3Bs7GVF+Ry/tu/45vxOfaUc97El4KQJxg25XEKz8WS9O4K3A6lNt/f4Pj4Y4yZ9HiP7G/oKFVVpRyMeRPTTd/hlldPmY0JlTPuZcRv/ky1Ds5/808Csz/HiTIyNN6UhS0maNLjmPbQGlJ3i+6QJNzl/9/emYdHVd4L+P1lsocsJCQhISELEECWgCCLgAuIIKgo7tb1alGqVr1t79NWe29vW9vetm617kVbq9ZdUEsBEVDZCoSdhCUBsieE7HsmM9/945zUEQJkQmbmhHzv8+TJnHXemTnn+53zfd/5fUqViEgc8DnwEPBJV4KEKx65k7BgkHA6FYU5VeRsKOXI7gqc7Yq4lHBGTktk2MS4Xp0TqTvYWxrZt+IVYkgT88UAAB8MSURBVHa9SnVFDaUHIxlcAC0BcHhGGiPu+yGjMq39fIPVaHe0s3nZS9S+8Sbp+2tp9YeS4e1kptdQEzuOgOkPkTH1Kt3ecI5g+dxNSqkS8/8xEfkYmASUi0iCS3XTMa+LWewEaG5oI3t9CXu/KqahqpXgsADGXJTEyGkJxAyy1pPK3kQhNBxroGZdEPGlkYT3gwOT2xid1EJM/CACQhJ8rdjrsIkfYQNSCJkSSfuwoxzNiyR9vz81+wZQlNlO0vktKPTojX0NnwQJEQkD/JRS9ebry4FfAJ8AdwK/Nf8v84UfFugVfCy/jj3riji09RiOdidJI/oz7bphpI0dcFJ+pL5EbUUxW176JeGfrCem3kHJwACOPnQVF931U2rzczj++ZNMr/wA3vqAjaEXETDjYSZMubTXDHHqC1paWti54jWid/+Z8c48qokg97zbmflfj1DbVE3Wi78icdVu/B78H9Yl/47g79zIpFsfxhZoja7TGs/ik+omEUkHPjYn/YG3lVJPiEgM8B4wGCgAblBKVZ1uXz1e3fRhLnvXFXHfc5f02D67itPhJHf7MfasLaLscB3+QTZGTBnImIuTiE4M87qPlSjcv409z/+axHU5BNnh8PAIIu+8nSkL7j8pn1Jd2REO/+NJhhV+QBjN7LSNpiZzEZPn3EpIUN+qljsd1ZXH2PfpHxl29C3iqaLAL4njY+5l9NxFBIZ8+3irq6tg/ZJfE/rhauKPt1MTaaPpmkuZvOhx+sXE++gTaM4Gy7dJ9BQ9HSQ2fpjLHi8HCXurg+wNJexaXUh9VQuRcSGMuSSJEVMTzjpXUm9nz5r3KXz1BVJ2lOH0gyNTkhly3yOMmjzvjNvaG6s5sPx54nP+QqyzgiMkcij9DsbNv5+4mJ4f+6G3UJC7j6IVT5FZ8Slh0sq+4PHI1AcYOWMhcoYGfnt7Gxs/fI7Gv71DWm4DLQFQPnM0mQ8+TsKwTC99Ak1PoINEN9n4YS671xVxvxeCRHNDG3vWFrFnXTEtjXYShkQyfk4KqaNjkD5cPWK3t7Lp3WdoefN9ko820hgslF2eyfnfe4zE1NFu70+1t5H35dv4/+t5UtsOUqX6kRW7kJS53ydjqLWypXoK5XSSs3U1zV/+kXGN63Hgx57+sxkw+1FSRk058w5O3J9S7N64jPxXniNtawl+CgrHJ5Jy30OMvPgaD3wCTU+jg0Q38UaQaKhuYfuqAnLWl9Bud5I6dgDnXz6YhKFRZ974HKampoxNrzxBv4/XMqDawfFof5qvm8WF9/6sZ8ZvUIrS3V9Qs+YZhtespx0/NodeSvBFDzJx8iXnZLtFe1sru1a9QfjOV8hoP0gtYRwYdD1DrvwBMQmnH3u7q+Tn7WDni79m0Oq9hLVASWo/wm67hQtuelCn5bAwOkh0k40f5bJ7rWeCRGNNK1kr89n3dTE4IWNyPONnp/T59ob8g1nseuk3JK7ZR1gLFKf2I/T2m5l040MeG2ugvvgAR5c/yZDiZYTSwg7baGrPoXaLhtpKsj/9Iym5fyOeSgolkZKRdzN23v2E9PNMgr+amjI2vmoE+dgqB1X9/WlZeBlTFj1GWKRO8Gc1dJDoJhs/ymX3miLu/9MlPbbPxtpWdqwsYO/XxTgdipFTBzLhilQiBvTdh5GcTic71rxDyWuvkLajHIDCCYNI/e6DXq2usDdUcfCffyIu56/EOo+TTwKH0m5j7JWLiYvpfSOllRzOpnjFU4w69gmhtLI3MJO2CxYzbuaNHsuoeyJt9hY2vvssLW+9T8qRRpqChGOXjWX89x5n4BD3qws1nkEHiW7Sk0Giub6N7Svz2ftlMQ6HYviUgUy8IpXI2L4bHNpamtjw1h9of2cpSYXNNAYLx2Zncv79j/m0AFHtbeR99Xds/3qBtNb91Kgwtg9YQNKcR8jIOG1mGJ+jnE72bV6Jff1zZDZupB0/dkTOImrmIwwf572xOU7yUoqdX75P0Z9fJC2rDIDCC5JJu+9hRkyb7zMvjYEOEt1k08e57PyikMV/urTb+7C3Odi1upDtq/Jpb3WQMXkgE+elEhUXeuaNz1Eqy47yr5d+SdTyzfSvc1IxIAD79XOY+h8/JTTCQj2NlKJ075dUffEMI6rX4cSPLaEzCJz+IBOmXmapdouW5kZ2rXid/ntfJ8ORSw39yB50A0PmPUL8oFRf632Lowe2svvF3zJobTahrVA0JILwO77DBdctxubf+6v3eiM6SHSTswkSToeT/ZvK2PLpYRpr20jLHMDUa4fQf2DfbXPI3bGWnJefJHl9HkHtkD8iiv533MEFCxZ5rfqju9SV5HJk+VMMLfqIMJrZ6zeC2tF3cP7cuwgJ9d1vWlFwkMMrniOj5GP6U0+BXxLl593NmHn3ERxq7QSG1ZUlbHr1V0Qs/YqYGgcVMf40LZzJ1HseIzIqztd6fQodJLpJd4KEUoqjeyrZ9HEe1aWNDEyP4MKFQ/tsb6X2djubl71E3d/eJm1/DW02KJyWzvBF/8mwibN8rec29qYa9i9/gejsvzHIWUI14RxKuJrk2d8jId07VWQOh4M9Xy+FLa8ytnEzCtgVNo3AC+9j1IVXIl4amrSnaGtrZvM7z9L21gcMym+kIRiKLx3JmO/+iCHnTfW1Xp9AB4lusunjPHauLmDx810LEsfy69jwQS4lh2qIig9lyjXppI+LRSyWA8obVFUWs+n13xK2bB3xFe3UhvtRPX8qFyz6SY+NguZLlNPBgc3LadzwMpkNG/AXJ3uCxtOSeRdjZt5EcHDPtzUdL87j0KpXSC5YSpIqo4oIchIXknL5AySlZvT4+/mC7HUfUfDqCyRtL0YUHBkdTfTttzHpyu+e9DS9pufQQaKbdDVINNW1sXlZHjkbSwnpF8CkK9MYOT0R2zkyfoM75GxfzYElzzL461xC2qAsOYzAWxYy6dZHCfBAwWkFSgoPc/TzlxhS8AHxVFJDP3IHXEa/iTeTccHlZ1WVVl9bRfa69wjNfpdRLTvwE0V20Fhax9zGqMvuIPAc/U4rjuaw/ZX/o//KrYQ3OimPDaB5wSVMvecnRPXXCRt7Gh0kusmZgoTD4WTP2iK2fnaE9jYnY2cmMXF+Wp9Ln9Fmb2HT+8/R+M4HpB2sw26D4kkppN/zABnTruwzd1LOdjvZ6z+mOesdRtWtJ1RaKWMABQNmEDxyDqkT5xARGX36nShF6ZFs8rNWEHp4JSOasgiUdsoklrzEqxk8816Sh5znnQ9kAVqbG9j61tO0vbuMhMJGmoKgeEYG5y36AUPHdm+oVc3J6CDRTTYtzWPnqgIWv3BykCjIrmT9e4eoLmti8Khopt8wrM81SpcXH2Lbkv8jcvlmYmocVEfaaJg/jYn3/pjoxL47JChAfV01B758j4CcpQxtzCJMWrErG0W2JCr7DaM9IhlCohH/QBytzdgaSgmqz2dQ8wFiMYZyL5Z48mMvJWbidWRMnHXGXErnOjlfLePIay+QvKUAfyccGRFF2C03MHXhAwQG6Cy0Z4MOEt2ksyBRW9HE+vdzObr7OBGxIcy4YRgpY2L6ztWycrJ17TuUvrGEtK0lBDqgcFgkEbfczMTrF+OvU0afRHNzM3nbV9O4fw1BFXtJaM1jgLMKm3xzvjWpIMpsCVT3G0pr4mSSxs0iOWN8r2uE9gbHi3LJeuU3RP5zM5H1Tir6+1E5bxIT7v4hg5NG+VqvV9Krg4SIzAWeBWzAn5VSvz3Vup4MEm0t7WStyGfn6gJsNj8mzkslc2ZynxnPoaKqiE1v/oHAT9eSUthmZPy8eCTD73mYtPHWHwPcatjb22muq6KtrZXQsH6E9ouy3CBXVqe9tYWs956n8Z33Scirpc0fcsfHMeCWW5k6524CbTpXVFfptUFCRGzAQWA2UARsBW5RSmV3tn5PB4nNS/PYvqqAWXeOZNNHuTTWtjF88kCmXjuEsKhz/4rZ7rCzef27lP/9TVI35RPWCpVxwbBwLhPu+iEhUb0vVYXm3KRox3r2v/ZHor/aS0irojjOn/orppB5+/dJTxrjaz3L05uDxFTg50qpOeb0TwCUUr/pbH1PBImsFfkAxKWEM+OmDAamR/bY/q1KXlk2G//2e6JWbmNoUTt2G5RPTiftjvsYcvFVfaZqTdP7aKuvY8fbf6Tlw0+IK6in1R+yx/cnaOFVzL/qUYL9g32taEl6c5C4HpirlLrXnL4dmKyUetBlnUXAIoDBgwdPyM/Pd/t9mupqyd+z86T5ZXk15GZVkJY5gEHDouCENAxuF5VuFq7uF8bu7v/b06WH91C6diWROUUE2aExMgjb5PEMvngegR3ZQt39DO5+S25/ZPcDlvtOZ/e9dmELdzdw/9hwe3UPf0fu7b1b71F5aC/Fqz6l357D+NsVVdH+OKdOYPjl1xEWfvJdsMePVTx7TodFRRGTNNjN/Zvv0ouDxA3AnBOCxCSl1EOdrd/dO4mSgzn8/Wc/OitXjUaj8SUjp1/CvId+2K1tuxokrNi5vwhIdplOAkp6+k1iU9O5++mX3NrG7YDqdvx1bwN3fJxOB9lfLaPy838Sv6+cQAdU9fendcb5jFxwB9GxyZ1u5/5H8Nxn6C7u/26ederWZ3b3e3X3l3P7K3L3M7j/md1/i5M3cConBXs2ULryU6L3FhNkh7pwG40ThzNk7g3EJrsxMmG3fjd3V3dvg9AIz6f+seKdhD9Gw/UsoBij4fpWpdS+ztbv6TaJc42j+zax783niVqzk+haB01BQtmFQxl8012Muuga/HR3S00fobmhhi3vPkfzsn+QfLAWP6BgWCQhC+Yz4aYHCAs/w0OP5xi9troJQETmAc9gdIF9TSn1xKnW1UHiZKqPFbD93RdwLP+C5CMNOAUKRvYn7OqruODGBwgJ9czIZBpNb6Hw0Hb2vPEc4au3MaC6neZAI2NA/DU3Mn7OdwjoAw/q9eog4Q46SBg01B5n2wcv0vTPVSRnH8ffCcdiA2iaPZlxtz9MQpoeEUyjORGH08GOVW9R/v7fSdxylGA71IT7UXnhCFKvu43zpi84Z++2dZDoA7Q2N7Bt2avUfPYZg3aVEGSH6gg/qqedR+r1tzNy6pXn7AGu0fQ0zQ01ZH30MnX/WE7S3mMEOOB4dACNF49jyPV3MvT8medUV3AdJM5RGuur2PnpX6j+fCUDdxQS1qKoDxWOTRlK4oIbybzsZmw6vbJGc1ZUVxSR9f4L2FeuYfDBWvwUlMYH0jRtDOlX3cyIyfN6/QWYDhLnEFVl+exatoTmL9aRmF1BUDs0hAjl45OJvepaxs+7g8Cgvjs0qkbjScoLD7D7vZdRazYw6HAdfgqqomzUTB5O0vzrGHXJdb0yf5kOEr0Yp9PJoazVHFn1EbJpB4Py6rApqI6wUTV5GAPnXk3m7FsICNRPkmo03qSy5DC7li6hZc2XDMqpJNAB9aFC+bhkwmdcxNgrbiN6YIqvNbuEDhK9jOpjBexZ8RZ1X35FzO4CouqdAJTHB9I4ZTSpV93IyAuv6vW3uBrNuUJtTTk7P3mdhi/WEL+riLAWhRMoGRxK+wWjSZp1JSOmX2XZizkdJCxORXEuB9Z+TM2WTYTsO8LAkhb8FDQGC+XnxRMy/UJGzL2ZxHSdqEyjsTqOdjvZ6z+h6PNP8N+2l8SCJvwUNAVB2bBo/MaNJmnabIZPnUdgsDWqhnWQsBCtzQ3kZn1BadZ6WvZlE55TSFyF3VjmD2VpkTgzRzBo5jzOm361Za88NBpN1zhefpTsle/SsP5r+mUXEHv8m/O9PDWCttFDiMqcQPoFl5GQPsYnNQQ6SPgAh6Od4kM7KMneSu2hHOx5hwnJKyGutAV/o/aIxmChYmg0Mm4UCRfOYuTUKwkMscaVhUaj8QzHig6yf91SarZsIiw7n/jiZmxm0VsfKlQNjqJ92GBCh2UQO3wcKaMvJCJmoEeddJA4A/a2FpobaggOizhjz6B2exv1NeXUV5bRWH2MhspyGkryaSkpxnHsGH4V1YRU1BN9vJXA9m+2awgRKgdH4shIIXzMeFInzSQpY4JuV9Bo+jjNDTXkbvuC0u0baM3OJvRwGXFlrf++mASo6+dHTVwobbGR+MUNICg+gdDEZMISkgmPHURU7CAiYxKx+Qd0y0EHiTOwe90HBNz/MwAcAnZ/sAcIDpvg51T4OcHPqbA5Ich+6v00hAj1UYG0xvTDmZxAcPoQooePIXn0FKIHpumAoNFouoS9pZnCA9sozdlGXW4O7UcLCC6pIrSqiYg6x7cCSAeHJyYy/80vuvV+vTkLrFeITRnB7jsvxdncgrO1FdXWCq1t0O4Afxv4+YHNhvj5IWFh2MLDCYiMIiiyPyFRA+ifPJS4lBF9LimYRqPxDAHBIaRnziA9c8ZJyxyOdo6XHqai4ACNpYU0V1XQVl1JdJobWWy7SZ+9k9BoNJq+TFfvJHRdiEaj0WhOiQ4SGo1GozklXg8SIvJzESkWkZ3m3zyXZT8RkVwROSAic7ztptFoNJpv46uG66eVUn9wnSEi5wE3A6OARGC1iGQopRy+ENRoNBqNtaqbFgDvKKValVJHgFxgko+dNBqNpk/jqyDxoIjsFpHXRKS/OW8QUOiyTpE57yREZJGIbBORbRUVFZ521Wg0mj6LR4KEiKwWkb2d/C0AXgSGAOOAUuDJjs062VWn/XOVUq8opSYqpSbGxsZ64iNoNBqNBg+1SSilLuvKeiLyKvCZOVkEJLssTgJKzrSPrKys4yKSb04OAI67oeottJd7aC/30F7uob0MujTwhdcbrkUkQSlVak5eC+w1X38CvC0iT2E0XA8Dtpxpf0qpf99KiMi2rjwc4m20l3toL/fQXu6hvdzDF72bfici4zCqko4C9wEopfaJyHtANtAOPKB7Nmk0Go1v8XqQUErdfpplTwBPeFFHo9FoNKfBSl1ge4JXfC1wCrSXe2gv99Be7qG93KDXJ/jTaDQajec41+4kNBqNRtOD6CCh0Wg0mlPSq4KEiMT42qEzRCTM1w6dISIpIhLla48TERH9BKQbiEiX+rN7GxEZJSI/9rVHZ7hkcrAUVi3DTkevCBIi0k9EngQ+MTPFTvG1E/zb62ngTRG5zSons+n1FPAPjGdOLIHL77hCRJ4QkWm+doJ/e71gxczDIjIdOCIic33t0oHL7/gZcJOIDPG1UwciEioiz2McYw+JyHhzvk/LOquWYV3B8kFCRCYAn5qTi4EI4EbfGRmYJ+/XQDPwGjADuMWnUoCIXABsAKKB8UqpbB8rASAi/sDzGN2u78B4TmaWT6UAEUkEXsU4pm614F1OBFAFLBaRzlLXeBURuRz4EmgFbgP2A3U+lfo2/wnEAHcCwcDLAEqpTkaI9g5WLcO6imWDhIjYzJdlGA/W/UAptRsIBLb70CvAfFkNvKCU+qlS6lNgJ8bBiS9OZhevViAPIx27XUTGiUiqWUh7HZf3jQfSlFIPK6X2AaFAji+cTqAJeAYjJUx/YK7LsecTRMTP5RhSGIVxBPCIudyXwSIXuNk87jcAozGDva+8Ou4SzGMtAHhbKbVfKfV7oNy8q/b63YRVyzB3sVyQEJFhIvIy8JiIpCulipVS2SISKSLPYkTgy83qgXgveo0QkdeB/xWRFLOg+4tLIViMmQtFebFfcSdeu4G1wKMisg54Dnga40l3r9WHunj9wvQqBlpE5K8i8i/gSuBqEVkqIgle9BomIv/VMa2UqgF2K6Wagb8CN9HFnDae8jrhqnc8kAZ8D7hXRCaZ0173MjmilDokIoHm9JsYY8B49bh3cXsZeNwsK9oxLtQmu6y2GLhDRJK8dTdh1TKsu1gqSIjIo8AyjNQcccATItJxQtQBS5VSyRipPOzA/3rJKwZ4HSPPlA2j4LtRKWU3D0wwstru84bPabyeEJH5wBKMrLpLlVIzML4nG3CPj7x+bVZTzAH+AuxXSmUA9wL5wH97yetWYA3wIxFZZM6zmQECpdSHQA1ws8udmU+8+CYbwiFgl1IqB+O73Ixx9e4rLz8ApVSbOR3QMc+bV+qdlBW/FaOTxrMYVXMxpmcRRiD7ro+8LFGGnQ2+GpnuVFQBdyul/mWepF9gVAMcMa9S1gIopZpFZDsw2EteI4AmpdST5olwKzBLRPabV+4ACcBHACIyC9irlCr3gde1wAHgXqWUHUAptVNE6oFKD/uczut6jIAQDrSZXq0ish4vFXoYmYbvxGhHel5E3lRKNZnVJGJeaT6DceJ+aF61b1FKHfCFl7ksA/i5iNiBXUAQkOVhn9N6mb+pMs/JrzEuSh7zcr1/Z2XF+UqpNSLyIcaQBB31/gcxCmRERDx8x2PVMqzbWOpOAiMCbxGRILOAy8Woi/0WYiQIXIRR9+4NtgNBIjLBPBE2YJxA15g+gtGLaLiILMdomPXGCdOZVz5wS0eAMP3GApdijN/hDTrzKsb4vjYC80RkjohMxGhoLPaGlFLqK+BrpdQmYA/wc3ORX0cBp5TahvHbbcO4+rR3sitveYHRg2gPsFgpdRPwLvB9Tzudwcu1oN0D7BXv99bprKwIN5d9DxgoIv8tIjdi/I4t4JUqMauWYd1HKeX1P4yeN6Eu03KK9bYDmS7TccBbwFbgJg94JQAZncyPAH4K/Mpl3gLg90AkxtWeE1iH0ahnBa/+GPWz72OkXPe11zXAk+br/zB/xx3e/B1PWGc0RgE3xmWeDWOc9UKMQOtzL4wA5rrczwpeLvPTMXoQJfW0V1fdzPVOLCtGYtwFrQK+YyEvj5ZhHvkNvP6G8DhGfd1rwI/NebYT1hHgPGC5y/Rw8/U8D3nZMAY5WgIM7nhfl+WXmcvmm9MZ5g8dZE4/ajGvYHP6Rgt5bXPxCvWW12nW/QWwxHw92TzOBlnQy+ayrNMLKh95BXjiu3LHrZOywg+jurPHg+lZenm8DPPUnzcbmvqLyIsYPSEux2hgekhEwpXLuBEudYaRQL55u7gLmGMuW+4hxRjgGEYX0gkiEqiUUi6NcVsx+of/UkSSgbEYVU5RAEqppy3mFWF6vWchr0Kgn+nV1Mk+PeJ1mnV/CSwQkQZgHkZh7Kmqr+56zcWlWtg8N6zgdQV4vNvrad1OUVbsxOiSazUvb5RhHsHjDdci0k8p1QA0Ar9XSh0253c8YBIG1Hes73ISXIHRAyAaeFAZ9aOe8OqgHmN0PICZGAfbEfVNXXUt8IYYT1X/GqM30yLVw43TPeh1zKJePTo8Y1e9TthGMIaKXGIue1QptV57+c6rO25WKyu87eU1PHhLFoPxhO2HGE8iJ5jzAzD6oudj9EBYD1yhvrlV7EhfPhu4x8NeN2NWL2DcPr9hvn4e48GlW4Eoc55rVUqPV5VoL+94uWwfCMzVXr716iE3S5UVnvby9p9HqptEZDZGV69yjEaaS80vGWW0+G9VSqUopRZjNHo9ZS5zKvPbVUp9rpRa4mGvmRgBC4xGuaPm62jgD8A0ZTxsRYeX+bpHq0q0l/e8zO1FKdWmlFqhvXzn1RNuYL2ywpNevsBT1U3HMXqy/BX+ncUy2nwtyqxyMnkXuE1E4lQPV5G44wVMAq4TkSsx6vTfBPaLSLRSqkp7nTtergFMe/nU66zd+qCX1/FIkFBK7RCRgyLS0f+8Asg0l6mOhh0RGYbRXfOQFwLEqbzGmcvWichXwGdKqX+IkSjverzwLIn20l590cvKblb18gVnNXypiIQC04Av1TeP6Xe23u+AcqXUky7bfRcjTcTLSqnnuy3Rg16eRntpr77oZWU3q3pZiW5HPhG5FyNN8BTMbo3mfD+X1zazZ0QcxlO3HamGFbAcmOKBANFdr9kiEo2H0F7aqy96WdnNql5Wo1vVTWL0e78GuFwptd9lfkf+G0RkuFLqgBh9iGOBcSLyS4xeTVuUUofOXr9HvY5iPOzV42gv7dUXvazsZlUvS6K63h3M3+X1WOA18/UwjERaaeb0IOAdjG5jMRj1eE6Mx+Ov7er7aS/tpb16n5eV3azqZfW/Ln2xGF28ngVmm/OmASsxRmPbCryEkQ3yFoyW/8dP2MfDnvjBtZf20l7W8LKym1W9esvfmb5cAV7A6OL1HWA1RjZKMOrn1vJNIrLLMLqD+blsH+ihg1F7aS/tZREvK7tZ1as3/Z3pC47ASO0cbk7PMb/w2cD5QANm0ipz+TtASseP48EDUntpL+1lES8ru1nVqzf9nbZ3k1KqDqOR5i5z1gaMW7OFSqntwBvAPSIyX0T+gnFbV2xuq06377NBe2kv7WUdLyu7WdWrN9GVLrAfY7TqJygjydUuwC4iQ4EfYwxveAuQp5S6Xn0znKen0V7aS3tZx8vKblb16h2c6VYDY3CN3wE/cZm3HpjuMu3RvPLaS3tpL2t7WdnNql695e+MdxJKqVJgKXCFiNwgIqkYQwG2uazj8SEetZf20l7W9bKym1W9eg1uROMrMEaT24+RG93nEU57aS/tZS0vK7tZ1cvqf27lbhKRACOuWKvOTnu5h/ZyD+3lPlZ1s6qXlTmrBH8ajUajObc5J1PbajQajaZn0EFCo9FoNKdEBwmNRqPRnBIdJDQajUZzSnSQ0Gg0Gs0p0UFCo9FoNKdEBwmNRqPRnJL/B76LMuprDrg9AAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "times = pd.date_range(start=datetime.datetime(2014,6,24), end=datetime.datetime(2014,6,25), freq='1Min', tz=tus.tz)\n", "\n", "solpos = pvlib.solarposition.get_solarposition(times, tus.latitude, tus.longitude)\n", "print(solpos.head())\n", "solpos.plot()" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(0, 100)" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAD8CAYAAABthzNFAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xd4VHX2+PH3uTOTQkLvbQmsFKmhCAoKCAuy6iIoItgQxbq2tayyuiK7ukX96tp2XRXFXkBXVNQfyCqsuitNpCNIUYoQIIQUkkz5/P64d4YAAQJk7txJzut55snMnTtzz0ySOfNp54oxBqWUUupgVqIDUEop5U2aIJRSSpVLE4RSSqlyaYJQSilVLk0QSimlyqUJQimlVLniliBE5AUR2SEiy8tsqycis0VkrfOzrrNdROQJEVknIktFpEe84lJKKVUx8WxBTAWGHbTtbmCOMaYtMMe5DfBLoK1zuQb4RxzjUkopVQFxSxDGmHnA7oM2nwe85Fx/CRhRZvvLxvY/oI6INI1XbEoppY7O7TGIxsaYbQDOz0bO9ubAj2X22+xsU0oplSD+RAfgkHK2lVsDRESuwe6GIiMjo2eHDh2O+6DBbdsI79lD2sknH/dzKKWUK4yheMVK/I0b4W/Y8ISeatGiRTuNMUd9ErcTxHYRaWqM2eZ0Ie1wtm8GWpbZrwWwtbwnMMY8CzwL0KtXL7Nw4cLjD+bPf2bPO+/SfuGC434OpZRyQ6S0lDVdu9HwN7+hwbXXnNBzicimiuzndhfT+8A45/o4YEaZ7Zc7s5lOBfKiXVHxJRCJxP8wSil1oqKfVVJeh0t8xK0FISJvAAOBBiKyGZgE/AV4W0SuAn4ALnR2/wg4G1gHFAHj4xXXASyr/H4spZTyGqfytlhVIEEYY8Ye5q7B5exrgF/HK5bDsrQFoZRKEgloQVTrldQimiCUUskhduoece9ju1onCMQq864rpZSHGefLrItdTNU7QVgWekY9pVRScHo7RLuYXCJoF5NSKjlEv8xqF5M7xNIuJqVUcoj1dmgLwiXOGIR2MymlPC/6OaVjEC6JZmJNEEopr9Npri6zNEEopZJEbKGcjkG4IvZG60C1UsrjTEQHqV1mtyB0DEIp5XnRdRDu9TBV8wQRbUFoglBKeZ12MbkrVvRKu5iUUl6ng9QuE00QSqnkoLWY3Oa80drDpJTyPK3F5LLYNFdtQSilPE5rMblLtItJKZUsYiuptYvJHaKzmJRSySG2DsLFea7VO0FY0TEITRBKKa/TWkzuir7P2sWklPI6HYNwl+hCOaVUstAxCJc5mdhoC0Ip5XFai8ltOkitlEoWWovJZVruWymVLLQWk7t0HYRSKlkYrcXkMu1iUkolC63F5DJdB6GUShY6BuEyXQehlEoW0XUQOgbhDl0HoZRKGroOwmXRct/aglBKeZzWYnKb6DRXpVSyqCa1mETkNyKyQkSWi8gbIpImIq1F5GsRWSsib4lIStzj0HUQSqlkUR1qMYlIc+BmoJcxpjPgA8YAfwUeM8a0BXKBq+IeTLQvT7uYlFIeF+sKrwZjEH4gXUT8QA1gGzAImO7c/xIwIv5hRGsxaQtCKeVx1SFBGGO2AI8AP2AnhjxgEbDHGBNydtsMNI97MLG+PE0QSilvM+EwAOLzuXbMRHQx1QXOA1oDzYAM4Jfl7Frup7aIXCMiC0VkYU5OzonFEs3EzhuvlFKeFake01x/AWwwxuQYY4LAu0BfoI7T5QTQAtha3oONMc8aY3oZY3o1bNjwxCKJrqTWLiallNdFnBZEFU8QPwCnikgNsYfjBwMrgc+AUc4+44AZcY8kNkitLQillLeZcHQMogp3MRljvsYejF4MLHNieBa4C7hNRNYB9YEp8Y4l2pdntItJKeV1Ti0m8bn3se0/+i6VzxgzCZh00Ob1QG9XA9FprkqpJBH7IlvFu5g8Q1sQSqmkUR2muXpKrAWhg9RKKW+rFtNcvUR0kFoplSyqyTRX77CiXUw6BqGU8rhqMs3VM2KzAbQFoZTyuNgXWe1icom2IJRSyUJbEO7SFoRSKlnEqrlqC8IlPm1BKKWSRDRBVOXzQXhK7IxymiCUUt6m01xdpgvllFJJI6wL5dylpTaUUskiWotJE4Q7tAWhlEoWOs3VbdGyuTpIrZTyOp3m6q7oNFej01yVUh4XO7GZtiBcosX6lFLJItqC0Gmu7tBifUqpZBEbK9UWhEt0oZxSKlnoNFd3aQtCKZUsjImAZWkXk2u0BaGUShbhiKutB6jmCUJbEEqppBEJuzrFFap5gohmY6MrqZVSHmfCEVcHqKG6JwifLpRTSiWJSERbEG4SSxfKKaWSg4noGIS7tAWhlEoWYR2DcJWIgIi2IJRSnmciYR2DcJ3Ppy0IpZT3RYx2MblNRPSMckopzzPhEOL3u3rMap8g8Pl0oZxSyvtCYVdPNwqaIBC/HxMKJjoMpZQ6IhMOg18ThKvE74dQKNFhKKXUEZlwCPFpF5O7An5MUBOEUsrjqksXk4jUEZHpIrJaRFaJyGkiUk9EZovIWudnXVdi8Qcw2oJQSnmc3cVUPVoQjwOfGGM6AN2AVcDdwBxjTFtgjnM77uwxCE0QSilvM6Fg1W9BiEgtoD8wBcAYU2qM2QOcB7zk7PYSMMKVeHSQWimVDELhajHNtQ2QA7woIt+IyPMikgE0NsZsA3B+NirvwSJyjYgsFJGFOTk5JxyMDlIrpZJBdZnF5Ad6AP8wxnQHCjmG7iRjzLPGmF7GmF4NGzY88Wh0kFoplQQ8OYtJRPo53/ARkUtF5FERaXUCx9wMbDbGfO3cno6dMLaLSFPnOE2BHSdwjArTQWqlVFLw6CymfwBFItIN+C2wCXj5eA9ojPkJ+FFE2jubBgMrgfeBcc62ccCM4z3GsdBBaqVUMkhEF1NF2ishY4wRkfOAx40xU0Rk3FEfdWQ3Aa+JSAqwHhiPnazeFpGrgB+AC0/wGBWig9RKqWSQiC6mihwtX0QmApcC/UXEBwRO5KDGmCVAr3LuGnwiz3s8xO/HFJa4fVillDo2wRDiwUHqi4AS4Cqne6g58HBco3JTQLuYlFLel4iFchVqQWB3LYVFpB3QAXgjvmG5RweplVLJwJOzmIB5QKqINMde4TwemBrPoNykg9RKqaTg0VlMYowpAs4HnjTGjAQ6xTcs9+ggtVIqGXh1oZyIyGnAJcBMZ5u7UcaR+P2gC+WUUh7n1S6mW4CJwL+MMStEpA3wWXzDcpEOUiulkkEo7PospqOmI2PMPOxxiOjt9cDN8QzKTToGoZRKBiYUAq+tgxCRhtgrqDsBadHtxphBcYzLNTqLSSmVDEzYm4PUrwGrgdbAZGAjsCCOMblKWxBKqaQQDCKBE1qjfMwqkiDqG2OmAEFjzFxjzJXAqXGOyzUS8ENQZzEppbzLGINJQIKoSIdW9NNzm4icA2wFWsQvJJdpC0Ip5XXOZ5SkeC9BPCAitYHbgSeBWsBv4hqVi8TvB2MS0r+nlFIVYZxeDs+1IIwxHzpX84Az4xuO+8Rvv+EmFNIEoZTyJM8mCBFpjV2eO6vs/saY4fELyz3Rc7yaYAhSUxMcjVJKHcqzCQJ4D5gCfABE4huO+yTgvAVabkMp5VFeThDFxpgn4h5JokRbEDpQrZTyKC8niMdFZBIwC/u8EAAYYxbHLSoX7e9i0haEUsqbvJwgugCXAYPY38VknNtJr+wgtVJKeVHsC6wHE8RIoI0xpjTewSRCdAzCaEVXpZRHJaoFUZGV1N8CdeIdSKLEuph0kFop5VFe7mJqDKwWkQUcOAZRJaa5xs7xql1MSimPMqVOgvB7L0FMinsUCSQ6i0kp5XGebEGIiA/4vTHmFy7F4zoJpACaIJRS3uXJMQhjTBgocmoxVUkHrKRWSikPiiUIDxbrKwaWichsoDC60RhTJc4qF33DdR2EUsqrTKk9/Gu5XA6oIglipnOpksR5w01JcYIjUUqp8pkSO0GI1xKEMeYlNwJJFCvNPotqpFgThFLKmyJeSxAi8rYxZrSILMNeOV2WMcZ0i29o7pBUO0EYTRBKKY8yJfY6ZUnxSIIAbnF+nltmm2CfTe53cYvIZVa6tiCUUt62fwwixdXjHjZBGGO2OT83iUg2cDEwGtgAvONOePG3vwVRcpQ9lVIqMSIlJWBZ+xf2uuRIXUztgDHAWGAX8BYgxphKOaucs8ZiIbDFGHOuc2KiN4F6wGLgMjfqP1lpdpMtooPUSimPMiWlSGoqIuLqcY+0DmI1MBj4lTHmdGPMk0C4Eo99C7CqzO2/Ao8ZY9oCucBVlXiswxK/H/x+bUEopTzLlJRgpbjbvQRHThAXAD8Bn4nIcyIyGHsM4oSJSAvgHOB557Zglw+f7uzyEjCiMo5VEVZqKpHifW4dTimljokpLXF9BhMcIUEYY/5ljLkI6AB8DvwGaCwi/xCRoSd43L8Bv2X/+SXqA3uMMdHlzJuB5uU9UESuEZGFIrIwJyfnBMNwnjM9XVsQSinPipR4LEFEGWMKjTGvGWPOxZ7BtAS4+3gPKCLnAjuMMYvKbi7v0IeJ51ljTC9jTK+GDRsebxgHsFJTdaGcUsqzTEkpkoAupmMaEjfG7Ab+6VyOVz9guIicDaQBtbBbFHVExO+0IloAW0/gGMdE0tKIaAtCKeVRkX1FWDVquH7cipwwqFIZYyYaY1oYY7KwZ0n92xhzCfAZMMrZbRwww62YrNRUXSinlPIsU7QPKz3d9eO6niCO4C7gNhFZhz0mMcWtA9stCE0QSilvihQlpgXh7qqLgxhjPsceAMcYsx7onYg4rLRUIoVFiTi0UkodlZ0gqncLImEkLT1WDEsppbwmsm8fUh3GIDxh1/ewbDoYe6KUlaZjEEop74oUFWGl16CoNMTnS9dT8uoYWPdp3I9bPRPE6pnwzlUMePBDdhaUIKlp2oJQSnmSMYbIvn1YNWrww+4i7nvjc1LXfQwFO+J+7OqZIDIbAyCF29ldWIqkpWL26UpqpZT3mGAQQiGs9HT2FAVpQJ59R0blrAM7kmqaIBoB0JA8dhWUYmkLQinlUabInkBj1ajBnqJS6ste+46MBnE/dkJnMSWM04JoKHvYXVjKSTVqYPbtw0QiiFX1cmYwGGTz5s0U6ziLqsLS0tJo0aIFgUAg0aFUqkgsQaSTWxSkgbjXgqjmCSKP3YUlWJmZAEQKC/HVrJnIyOJi8+bN1KxZk6ysLNfLBSvlBmMMu3btYvPmzbRu3TrR4VSqiNP9bdWoQW5RqXYxxV16XYzlp6HsYVdhKVZNJ0Hk5yc4sPgoLi6mfv36mhxUlSUi1K9fv0q2kqMtCHHGIJr59kB6PfDHv3hf9UwQloVkNKK5fy+7C0tjrYZwQUGCA4sfTQ6qqquqf+ORomgLIoOdBSU09++Fmk1cOXb1TBAAmY1o6ttrtyAynBZEFU4QibRx40Y6d+58ws/zpz/9qRKiSYx58+bRo0cP/H4/06dPP+C+u+66i86dO9O5c2feeuut2PY5c+bQo0cPsrOzOf3001m3bh0AmzZtYvDgwXTt2pWBAweyefNmV17DwoULufnmmwH4/PPP+eqrr2L3XXHFFYe8LlU5IkWFgD0GsauglCbWHk0QcZfZ2B6DKCjFV1MTRDJI5gTxs5/9jKlTp3LxxRcfsH3mzJksXryYJUuW8PXXX/Pwww+zd689S+X666/ntddeY8mSJVx88cU88MADANxxxx1cfvnlLF26lPvuu4+JEye68hp69erFE088ARyaIFT8mDJjELsKS2hgdkPNpq4cuxoniEbUM7nsLizFinYxVdExCC9Zv3493bt3Z8GCBWzcuJEzzjiDHj160KNHj9gHzrZt2+jfvz/Z2dl07tyZ//znP9x9993s27eP7OxsLrnkEgBeffVVevfuTXZ2Ntdeey3hsH1G3MzMTO655x66devGqaeeyvbt28nPz6d169YEg0EA9u7dS1ZWFsFgkCVLlnDqqafStWtXRo4cSW5uLgADBw5k4cKFAOzcuZOsrCwApk6dyvnnn8+wYcNo27Ytv/3tb2Ovb8qUKbRr146BAwdy9dVXc+ONNwKQlZVF165dsQ6aJbdy5UoGDBiA3+8nIyODbt268cknnwB2l0k0WeTl5dGsWbPYYwYPHgzAmWeeyYwZRy58/P7775OdnU12djbt27ePDeIuWrSIAQMG0LNnT8466yy2bdsWe9133XUXvXv3pl27dvznP/8B7KRw7rnnsnHjRp555hkee+wxsrOzY/fPmzePvn370qZNG21NVKJwYbQFUYPc/GJqh3NjE23irXrOYgJKMhqQKwXsKijaP4spv+q3ICZ/sIKVW/dW6nN2bFaLSb/qdNT91qxZw5gxY3jxxRfJzs6mqKiI2bNnk5aWxtq1axk7diwLFy7k9ddf56yzzuKee+4hHA5TVFTEGWecwVNPPcWSJUsAWLVqFW+99RZffvklgUCAG264gddee43LL7+cwsJCTj31VB588EF++9vf8txzz3HvvfcycOBAZs6cyYgRI3jzzTe54IILCAQCXH755Tz55JMMGDCA++67j8mTJ/O3v/3tiK9lyZIlfPPNN6SmptK+fXtuuukmfD4ff/zjH1m8eDE1a9Zk0KBBdOvW7YjP061bNyZPnsxtt91GUVERn332GR07dgTg+eef5+yzzyY9PZ1atWrxv//9L/aYd955h1tuuYV//etf5Ofns2vXLurXr1/uMYYPH87w4cMBGD16NAMGDCAYDHLTTTcxY8YMGjZsyFtvvcU999zDCy+8AEAoFGL+/Pl89NFHTJ48mU8/3V/WISsri+uuu47MzEzuuOMOwE6M27Zt44svvmD16tUMHz6cUaNGHRqMOmYR50uC1KxJuHgLb9dO5/TUGrR04djVNkHMCO3ijy0bk/b9D0iN0wGIFFb9BJEoOTk5nHfeebzzzjt06mQnk2AwyI033siSJUvw+Xx89913AJxyyilceeWVBINBRowYQXZ29iHPN2fOHBYtWsQpp5wCwL59+2jUyF4AmZKSwrnnngtAz549mT17NgATJkzgoYceYsSIEbz44os899xz5OXlsWfPHgYMGADAuHHjuPDCC4/6egYPHkzt2rUB6NixI5s2bWLnzp0MGDCAevXqAXDhhRfGXtPhDB06lAULFtC3b18aNmzIaaedht9v/1s+9thjfPTRR/Tp04eHH36Y2267jeeff55HHnmEG2+8kalTp9K/f3+aN28ee8yRPPTQQ6Snp/PrX/+a5cuXs3z5coYMGQJAOBymadP93Rbnn39+7P3buHHjUZ8bYMSIEViWRceOHdm+fXuFHqOOLpy3F/x+inwpZFqb+FODejwqQU0Q8dSkdivYDrV9P5JvBcCyqkUXU0W+6cdD7dq1admyJV9++WUsQTz22GM0btyYb7/9lkgkQlpaGgD9+/dn3rx5zJw5k8suu4w777yTyy+//IDnM8Ywbtw4/vznPx9yrEAgEJvR4vP5CIXsU53369ePjRs3MnfuXMLhMJ07dyYvL++wMfv9fiIR+7TpB0+fTC1zfuDoMYwp9yy5R3XPPfdwzz33AHDxxRfTtm1bcnJy+Pbbb+nTpw8AF110EcOGDQOgWbNmvPvuuwAUFBTwzjvvxJLV4cyZM4dp06Yxb948wH7/OnXqxH//+99y94++vrLv39GUfU+O971Qhwrn5eGrVYvdhUEyA9vZDjSq485aj2o7BtGkXjsAaqVsY3dRECszk0hBYYKjqrpSUlJ47733ePnll3n99dcBu1+9adOmWJbFK6+8EhtD2LRpE40aNeLqq6/mqquuYvHixYD9wR8dQxg8eDDTp09nxw67YNnu3bvZtGnTUeO4/PLLGTt2LOPHjwfsxFW3bt1YP/orr7wSa01kZWWxaJF96vSK9Kn37t2buXPnkpubSygU4p133jnqY8LhMLt27QJg6dKlLF26lKFDh1K3bl3y8vJiLZDZs2dz8sknA/Z4SDRx/fnPf+bKK6+MPV+HDh0OOcamTZu44YYbePvtt0l3zkrWvn17cnJyYgkiGAyyYsWKo8YbVbNmTfKrwRcqLwjvzcNXuzY7C0pID9h/K43rt3Pl2NW2BdG0cVcA0gI72bG3hIaZmVV2oZxXZGRk8OGHHzJkyBAyMjK44YYbuOCCC5g2bRpnnnkmGRkZgD0Y+vDDDxMIBMjMzOTll18G4JprrqFr16706NGD1157jQceeIChQ4cSiUQIBAI8/fTTtGrV6ogxXHLJJdx7772MHTs2tu2ll17iuuuuo6ioiDZt2vDiiy8C9myh0aNH88orrzBo0KCjvr7mzZvzu9/9jj59+tCsWTM6duwY+2a/YMGC2AD4Bx98wKRJk1ixYgXBYJAzzjgDgFq1avHqq6/Guouee+45LrjgAizLom7durHxgc8//5yJEyciIvTv35+nn34asBNHed/cp06dyq5duxg5ciRgt0A++ugjpk+fzs0330xeXh6hUIhbb7011ro7ml/96leMGjWKGTNm8OSTT1boMer4RJwWxI78Enz+PYgx1K9/6BeBeJBkbgr26tXLRGeZHDNjOHVqZ07Ka8gFA16l2/03EmjRgpZPP1W5QXrAqlWrYt8+q7vp06czY8YMXnnllbg8f0FBAZmZmYRCIUaOHMmVV14Z+2COtw8//JD169fH1ipUR1Xxb33DBaPwNajPp5fdzaffXMyOjL18emXFW3vlEZFFxpheR9uv2rYgEKGJsQgFCti6p5gedeoQdqY3qqrppptu4uOPP+ajjz6K2zHuv/9+Pv30U4qLixk6dCgjRoyI27EOFh2YV1VLOC+PlDZt2L63mOLAPppJ/EtsRFXfBAE08aWzNVDMT3nF+OrVpWTV6kSHpOLIja6QRx55JO7HUNVLeO9efLVrs31vMXmBEO0D5U9njodqO0gN0DSlNrn+MNvyivHXrUdo9+5Eh6SUUjEmHCaydy++WrX4KW8vu3zQLD3+VVyjqnWCaJLeiD0+ix15O/HVq0dk71777E1KKeUB0Ykzvtq12Je/lrAIzWv9zLXjV+sE0bRmCwCChd/hr28vbgrpOIRSyiOivRpW3bpQuhaAZnVPcu341TpBNK/fHoBAeB2R2nUAdKBaKeUZoZydAOyrWZc0awsAzRqeeGXkiqrWCaJVM3uVambqFvJT7XpMYWfRkqo8Wu67apT7doOWFD9QKCcHgB2BTAKBnYgxNGl85PpelalaJ4j6DU4mIxLBStnJDp+9wjS0W1sQXpXMCaIqlPt2g5YUP1A0QWyxamBS8mgUgZS0I5dVqUzVOkGIZfEz46ckkM+PESdB7MxJcFRVm5b7dr/cN9jJplevXnTq1IlJkybFtmdlZcVKe/fu3TvWSrniiiu47rrrOOOMM2jXrh0ffvjhEZ9fS4rHR2hnDpKSwqYSH/kphWRZaa4ev1qvgwBoFajNN5HdrCux6JSaSuinKl6F8uO74adllfucTbrAL/9y1N203Peh3Cj3DfDggw9Sr149wuEwgwcPZunSpXTtapebqVWrFvPnz+fll1/m1ltvjSWDaGHD77//njPPPJN169bFCioeTEuKx0coJwd/gwb8uLuQ7SlhTk1v5Orxq3ULAiAroyk5ftiwYxeBJk0I/vRTokOqkqLlvl999dVY+e5gMMjVV19Nly5duPDCC1m5ciVgl/t+8cUXuf/++1m2bBk1nRM6lVW23Hd2djZz5sxh/fr1wKHlvqPlqidMmBCrs/Tiiy8yfvz4cst9RyueHkm03HdaWlqs3Pf8+fNj5b4DgUCFyoYPHTqUs88+m759+zJ27Nhyy31v3ryZ8ePHc9tttwH2Yry5c+fSvXt35s6dW6Fy32+//TY9evSge/furFixIvZeA7G6VGPHjj2guuvo0aOxLIu2bdvSpk0bVq8++kLSsiXF16xZEyspnp2dzQMPPHDAeImWFD+6UE4O/oYN2b5zNUWWReta7lRxjdIWRN2TiOSvZG/uN/ibNiXkNIGrrAp8048HLfd9ePEu971hwwYeeeQRFixYQN26dbniiisOeD3R9+pI18u7fTAtKV75QjtySG3ThqK9iyANWrs4gwkS0IIQkZYi8pmIrBKRFSJyi7O9nojMFpG1zs+6bsTTurH9bTZUvAy/tiDiRst9l8+Nct979+4lIyOD2rVrs337dj7++OMD7o/OnHrrrbc47bTTYtunTZtGJBLh+++/Z/369bRv3/6wr0NLilc+YwzBLVswTZoi5nsA2rTo62oMiWhBhIDbjTGLRaQmsEhEZgNXAHOMMX8RkbuBu4G74h3Mz7MGYc2fTGrKJkrqtieUk4MJhZAKnKFLHRst952Yct/dunWje/fudOrUiTZt2tCvX78D7i8pKaFPnz5EIhHeeOON2Pb27dszYMAAtm/fzjPPPENaWhpbt25lwoQJhxQ81JLilS+8cyemuJg9tRpgRXLIiERo2KiLu0EYYxJ6AWYAQ4A1QFNnW1NgzdEe27NnT1MZfvV8ZzP6qZ5mwZNTzMr2HUzJj5sr5Xm9YuXKlYkOwTOmTZtmLr300rg9f35+vjHGmGAwaM4991zz7rvvxu1YB/vggw/M448/fkyPadWqlcnJyTlk+7hx48y0adMqKzTXVKW/9cLFi83K9h3MJ89NNyP/0dWMmdKt0p4bWGgq8Pmc0K/JIpIFdAe+BhobY7YBGGO2iYhrw/VtA3X4xuxmI/XoBAR/2ERKi+ZuHV65RMt9q2QS3GyvnF5nZbAlNcS56W6chfpACUsQIpIJvAPcaozZe7QBsDKPuwa4BuzFR5WhY502zMrdw9J9QToBJRs3ktHX3b4+FX9a7vtQh5s9NHXqVFfjUIcKbrFnfH1fuoGidIvODdw/n3xCprmKSAA7ObxmjHnX2bxdRJo69zcFdpT3WGPMs8aYXsaYXg0bVk7Z2/aNewDwY/ESJD2d0gpOuVNKqXgp3bABf+PG5BbZEyU6tuzvegyJmMUkwBRglTHm0TJ3vQ+Mc66Pwx6bcEXHn9vTB4PhpQRatdIEoZRKuJK165DWbQha60kxhjatf+F6DIloQfQDLgMGicgS53I28BdgiIisxR60dm3Cfr0G7flZyBDNVsy8AAAfUklEQVRM20ZxkxaUrvverUMrpdQhTCRCyfr15DZsTmHaLk6K+AmkZLgeh+tjEMaYL4DDDTgMdjOWsrqkNOI/sp2tDVrS4vNPCeXm4q/rylIMpZQ6QHDLFkxxMetr1OWH1DAXZbRJSBzVvtRGVO8m3dnrs1juLwCguEwpApXcrrzySho1anRIyfFvv/2W0047jS5duvCrX/0qVhivtLSU8ePH06VLF7p168bnn39+yHMOHz68UkqYV9SECRNi5THKVrWtrFLqyluKnbIm62QjpZbQu8XpCYlDE4Sj+8/PBuC71DUAFK/QBFFVXHHFFbEKqWVNmDCBv/zlLyxbtoyRI0fy8MMPA/YCNYBly5Yxe/Zsbr/99tjKZYB3332XzMxMd4J3PP/887Eifslc9lxVTPGy5eD3sz7jO8QYenS8KCFxaIJwZGWdSaNwhNxa65FmzSk+hpIA6sh+//vf8/jjj8du33PPPTzxxBMYY7jzzjvp3LkzXbp0iZV8iJZ8jrrxxhtj0y6zsrKYNGkSPXr0oEuXLrECcjk5OQwZMoQePXpw7bXX0qpVK3butM/G1b9/f+rVq3dIXGvWrKF/f3tmyJAhQ2KlMcqW027UqBF16tSJlf0uKCjg0Ucf5d57763Qaz/77LNjZbBr167NSy+9RDgc5s477+SUU06ha9eu/POf/4y97oEDBzJq1Cg6dOjAJZdcElsZHS09Xl7Z83A4zNVXX02nTp0YOnQo+/btq1Bsyrv2LVtKJOvn5Gb+RNuwRa06R64QEC9aT8IhlsWpgSZ8Kj+R26wr/uXLEx1SXPx1/l9ZvfvoVTmPRYd6Hbir9+Grolx11VWcf/753HLLLUQiEd58803mz5/Pu+++y5IlS/j222/ZuXMnp5xySuwD+0gaNGjA4sWL+fvf/84jjzzC888/z+TJkxk0aBATJ07kk08+4dlnnz3q83Tu3Jn333+f8847j2nTpvHjjz8CdmmKGTNmMGbMGH788UcWLVrEjz/+SO/evfn973/P7bffTo0aNSr03kQX5S1atIjx48czYsQIpkyZQu3atVmwYAElJSX069ePoUOHAvDNN9+wYsUKmjVrRr9+/fjyyy85/fT93Qt/+ctfDih7vnHjRtauXcsbb7zBc889x+jRo3nnnXe49NJLKxSf8h4TDlO8fAVbu/Rkffo6Lk5PzPgDaAviAINaD6bIslhVq4Tgli2UVqFTOSZSVlYW9evX55tvvmHWrFl0796d+vXr88UXXzB27Fh8Ph+NGzdmwIABLFiw4KjPV16Z6C+++IIxY8YAMGzYMOpWYILBCy+8wNNPP03Pnj3Jz88nJSUFsMcsWrRoQa9evbj11lvp27cvfr+fJUuWsG7duli9oYrauXMnl112Ga+//jq1a9dm1qxZvPzyy2RnZ9OnTx927drF2rX2Cel79+5NixYtsCyL7OzsCpXBbt26dayE+rGUzlbeVLxyJZH8fFbW3EFYhKEdLkhYLNqCKOPUrlcQWPc6yxpt4jSg8MuvSLlodKLDqlRH+qYfTxMmTGDq1Kn89NNPseqj0e6Tg5Utsw2HL7Vdtkz04Z7rSDp06MCsWbMA+O6775g5c2bs+I899lhsv759+9K2bVvmzp3LokWLyMrKIhQKsWPHDgYOHFjuIHZUOBxmzJgx3HfffbHBZGMMTz75JGedddYB+37++efllhE/moMfo11Mya3wS/vMisua/0jdcIQuHcckLBZtQZSRUasZfSIZLG6xi3C9BhRW8/PhVqaRI0fyySefsGDBgtgHY//+/XnrrbcIh8Pk5OQwb948evfuTatWrVi5ciUlJSXk5eUxZ86coz7/6aefzttvvw3ArFmzYqcNPZJoqfBIJMIDDzzAddddB0BRURGFhYWAXWbb7/fTsWNHrr/+erZu3crGjRv54osvaNeuXSw5PPXUUzz11FOHHOPuu++ma9eusdYNwFlnncU//vGPWOny7777Lna8iihb9lxVPYX//S/7WmaxvGEhp/kb4gukHv1BcaIJ4iC/ajOM3QEfG5rXovB//8PoP2KlSElJ4cwzz2T06NH4fD7AThpdu3alW7duDBo0iIceeogmTZrQsmVLRo8eTdeuXbnkkkvo3r37UZ9/0qRJzJo1ix49evDxxx/TtGnT2JnoomdqW7NmDS1atGDKlCkAvPHGG7Rr144OHTrQrFmz2DkiduzYQY8ePTj55JP561//yiuvvHLU469evbrcU34+8sgjzJo1KzZQ/f777zNhwgQ6duxIjx496Ny5M9dee22FT5gD+8ueRwepVdURLihk3+LFfNcoQqFlMTJBs5diKlLy1auXyir3Xda+ghzT54WOZtJtPc3K9h1M/ty5lX4Mt3mhBHI4HDbdunUz3333XVyev7i42ASDQWOMMV999ZXp1q3ySiNXxDnnnGNKSkpcPaY6lBf+1k/EnvfeMyvbdzC/uberGfx8RxMOlcblOFSw3Le2IA6SltGAIf7mzOy4j1CNGuQ5J3BXx2/lypWcdNJJDB48mLZt28blGD/88AOnnHIK3bp14+abb46tZXDLhx9+GBvkVup45X04k3116/BpmxDDMttj+QIJjUcHqcsx/rTbee/LO1iV5afbp3MIFxTiy3S/DkpV0bFjR9avXx/XY7Rt25ZvvvkmrsdQKp5COTkUfvUVS7ql4hdh3MD7Ex2StiDK0+aks+gTSmNGj3wiRUXk/etfiQ5JKVXF5b7xBiYS4d2exZxp6tGwcddEh6QJ4nBu7nMXS1v62No0jd0vvYQJhxMdklKqiooUF5P7xpusywqwpb5wkwdaD6AJ4rC6dh7F6aEMpp9aQnDzZvLec+30FEqpaib31VcJ5+byWt8Qv5SmZLUelOiQAE0QR3T3oIdY0FbY3MRix+OPEykqSnRISqkqJpSbS84z/2RNa4tNLeC2s59IdEgxmiCOoFXr/lye2YN/DhHCO3aw47G/JTqkKqVsCesT5XZ11YNt3bqVUaNGAbBkyZJYDSaA+++/P+nOVa3cs/2BBwkVFfDsYOHXDYfQoOHJiQ4pRhPEUfx65PMEGvuZ1VPIfeUVCv/3v0SHVGWULWGd7Jo1a8b06dOBQxOEUoez95NP2DtzJu/2E1rUqcGlZz969Ae5SBPEUfgCqTx27ivMPAN+qgcbbr6ZUqfqp6qYwsJCzjnnHLp160bnzp1jZb2jJazBbgHcdddd9OzZk1/84hfMnz+fgQMH0qZNG95//30Apk6dynnnncewYcNo3749kydPLvd4Dz/8cKyU9qRJk44Y28KFC2OrnLt06YJ9ynT4/vvvGTZsGD179uSMM86IlRW/4ooruPnmm+nbty9t2rSJJYXoiXtKS0u57777eOutt8jOzo691pUrV8ZezxNPeKcLQSVO8Zo1/HDX3WxoAl/3Eh6+4G3E8tZHsq6DqIAmTbrw4Kn38cfCP3DfK/msGTeOk99+G3+DBokO7Zj99Kc/UbKqcst9p57cgSa/+91h7//kk09o1qxZrBheXl7eIfsUFhYycOBA/vrXvzJy5EjuvfdeZs+ezcqVKxk3bhzDhw8HYP78+SxfvpwaNWpwyimncM4559CrV6/Y88yaNYu1a9cyf/58jDEMHz6cefPmHbaMeK9evWKls++8806GDRsG2OUsnnnmGdq2bcvXX3/NDTfcwL///W8Atm3bxhdffMHq1asZPnx4rGsJ7JIif/jDH1i4cGGsNtP999/P6tWr+eyzz8jPz6d9+/Zcf/31BAKJXQSlEqd08xa+G38FhYES/j5S+L9Bf6N23axEh3UITRAVdEr2GG4tzufxwse4fdo2ll8wgk5vvE2gWbNEh+Z5Xbp04Y477uCuu+7i3HPP5Ywzzjhkn5SUlNiHc5cuXUhNTSUQCNClS5cDylcPGTIkVvPo/PPP54svvjgkQURLioN9gp+1a9ce9TwTb7/9NosXL2bWrFkUFBTw1VdfceGFF8buLykpiV0fMWIElmXRsWNHtm/fXqH34JxzziE1NZXU1FQaNWrE9u3badGiRYUeq6qWkvXrWXXxRYT2FfDEGGHSwAfo0HZIosMqlyaIYzDo1KvJSKvDk2YSN7yzi+XnDuPnTz5DrX59Ex1ahR3pm368tGvXjkWLFvHRRx8xceJEhg4dyn333XfAPoFAINa9Y1lWrIS1ZVkHFLKL7nO428YYJk6cyLXXXlvh+FasWMGkSZOYN28ePp+PSCRCnTp1Yi2Lg5Utr20qWGb8eMp4q6pn1yf/j82/vZ1if5h/XmTxuxGP0uXks47+wATxVodXEuiTfSG/v/IVXrgkQE6NID9OuIrFd9ykU2CPYOvWrdSoUYNLL72UO+64g8WLFx/3c82ePZvdu3ezb98+3nvvPfr163fA/WeddRYvvPACBQUFAGzZsiVW1nvw4MFs2bLlgP3z8vIYM2YML7/8Mg0bNgSgVq1atG7dmmnTpgF2Evj2228rHGPNmjXJz88/7teoqp7wnj0suG48O269lR/qhXn1kjQeuuYDTycH0ARxXLJa9uSJ2xey5KquzO0qpH/4KQsH9mHVs09j9JvhIZYtW0bv3r3Jzs7mwQcfrPD5nMtz+umnc9lll5Gdnc0FF1xwQPcSwNChQ7n44os57bTT6NKlC6NGjSI/P59IJMK6desOOTf1e++9x6ZNm7j66qtjg9UAr732GlOmTKFbt2506tSJGTMqvlDyzDPPZOXKlQcMUqvqyZSWsvSJh1kyqB/pc//Hx6fAD9eewZO/WUTDBok7lWhFSUWbyF7Uq1cvE50FkyhrN87n9Rd/Tc+5Bfz8J8it60OG/4pTbroHf4Ln5ketWrWKk0/2ztzq4zV16tQDBn+PxfLly3nhhRd49FFvTSNUlcsrf+uh3Fz++/gfSJk5i1r5EdY0h6WD6zN+3D9o2bxLosNDRBYZY3odbT8dgzhBbbN6M2nyAhYs/YCZrz7AyfP30ual91j++nts69iQ5iPG0WXU5YjOWEmozp07a3JQcRUpKWHJm1P46f1XabE6lwZhWNMc5p1Vj6GX3M+Ik705EH0k2oKoZBt+XMInr08mc/4aOq411CiFolTY1iaTQLcedBl1Lc0693A1Jq98q1Iq3tz8WzfGsPWb+Sx591lY+g3NN+wjNQiFqbCsg8W+Ph0Zecmfadr4JFfiORbagkiQ1i2zuf4uuzz4f+dPZ/H7z5O54kda/VBA3VXzyHtzHhtrwe7GaQR/1pza3XrT9ZeX0KDlzxMcuVLqSHauX83SWW+R9+3X+DdvpdHWEmoVQhsgNxOWnexjX6cseo+8kfFdhyU63EqhCSKOTus9itN624uotuVs4KuZz7B30VfU3LSbBjuKabz2e5jzPTmPvsH3GbCnro+ieumEGjUgrdXPaXJyL9r3+gU1Gzc/ZDrnsTLGnPBzKOVlldEbsnPLetYu+JRdq7+h5IcNWDt2USO3iDq5EWoVQVPsy446sKmlj4JW9anTqx+Dht9C39qNT/j4XqNdTAkSDO5j8cKP2fjVRwTXfkcgJ5eauSHq50KtfQfuWxyAvJpQWNNHcc0AoZo1sGrVIVCvEWmNmlOraWsaZLWjxUldqZFZ+5BjbdiwgZo1a1K/fn1NEqpKMsawa9cu8vPzad269QHbi/bksXX9MnZsXMGerevZt2sLoT27kNw8UvKLScsPUqMgTK0CSAse+Lx7MiC3NhTUTaGkcV1STmrHSWecS3aPs/H5kvf7dUW7mDRBeEw4HOS75V+wfuFn5G1cRThnO/68fNLyg2Tk23/ENYvAOsyvrSgVilOgJBWCqUIwxUe4bh0anj+BlEYtwPKBCFiCiCBigWUhYmFZPsTyIVb0uqUJRbnOGEMkEsZEIs5P+7qJRDAmApGI3Vooe4lEKN3+I7um/xMrN4+0YkP6PsjcB4HDnOsrZEFeJuRnQFGmj9KaqZjaNfE3bEJmm5Np3XMAbTuchj+QWv4TJLGkTBAiMgx4HPABzxtj/nKk/atigqiIgqJ8Nq1eyE8bllLw0wZKdu4gvCcXyS/AKizGV1KKryRMoCRCSqkhtQTSSiG9BPyRYztWSQBK/VAagLAPwpb9M+ITIj6IWELEZzm37YvxWUR8FsZvYXwWxu/D+HzgtxMQPh+Wz4/4/LHr+AJY/gDiT8Hy+RCfD/wBLJ8Pyx8Anx+f397H8vvtn4EAPl8A8QUIBPz4fAF7H18Af8CPJQEk4MdvBbB8KViBFHzixwoEsMRCfBZi+bAsC8vvw8K3f5vPAvHh8zv3Y2H5LAQnqcJxF1YzxjgX53okQsSEiYRDGBMhYgyRSMj+cDQRTCRMxESIhEIYEyYSdO4Lh+zHROxtoVCQcCSMCYUJh+zbJhIiEo4QCpYQCZYSDpcSDgUJh0oJh0NEQqVEwkEioRCRSBATDts/gyFMKIgJh+y1PaGQfT0YhnDIPsNiOAThCISi18MQDiORMIQjSCQCYYNEIlgRgxU2+ELOJQz+kMEfAn8Y/CEIONcP94F+OKU+KEmx/0ZLU6A0RShNsyhN9xNKDxDOSIOMTPx161OjfjPqNv85jVt3pmW7HqSlpR/X7zDZJV2CEBEf8B0wBNgMLADGGmMOe8KA6pogjlc4FCJv11Z2bt1IXs5mCnNzKNy7i2DhHoKFewnvK8CUFGNKSqG0FKs0iFUawgqGsYJhJBxBwhGssPPPHjZYEZzr2JeI8w/vXHzhY09KXhQBzEGNqQNuHvRvVB1XoIZl/5eIiOV8mXCuRywI+SHsE/unXwj7hLBfiPiFiN8i4vdh/H4I+CElBQmkYKWmYaWm40vPwJdRk0CtuqTVqkdmvcbUatiEeo1bUK/+z5K6uycRknEWU29gnTFmPYCIvAmcB1TOGWUUPr+feo1/Rr3GP3P1uJFwmJKCfIqL8ykqLKC0uJDifYWUFhcSChYTKi0iVLKPUKiYcGkJ4VCQiHMxoZD9zdn5JhsJhzChMDjfoAmHiUQi+7shnC6IiIkgxrntdEFErxtjkOgXI6eLQkz0OvZ9zu3YdgwSiSaJ/dmg7LVDvmpJ9H6DxPYV+7qIk2DEyTSyP+GI7L8lsv8+EYwIxhJ7e/Q6TpehZR1wn1iW/dNpDYnPj1iW/WFq+e1WnOVzWmwBLJ8fyxfA5/PbLblACr6UFKxAKr5AKr7UVAKp6QRS0wikpZOamk4grQZpGZlk1KhDanpNLMt3wn8vyju8lCCaA2VPtLAZ6JOgWFQlsnw+0mvXIb12HeomOhilVIV5KUGUNxp66JcykWuAa5ybBSKy5gjP2QDYWQmxuUFjjY9kiTVZ4gSNNV7cjLVVRXbyUoLYDLQsc7sFsPXgnYwxzwLPVuQJRWRhRfrZvEBjjY9kiTVZ4gSNNV68GKuXxtIWAG1FpLWIpABjgPcTHJNSSlVbnmlBGGNCInIj8P+wp7m+YIxZkeCwlFKq2vJMggAwxnwEfFSJT1mhriiP0FjjI1liTZY4QWONF8/F6pl1EEoppbzFS2MQSimlPCSpE4SItBKROomOoyJEpIOIXJXoOCpCRGomOoaqSETqJzqGihKRNiKSFDXoRaSTiNyd6DiORkSSbhlQUiYIEckUkUeBmUCzRMdzJCKS4cT6OXC+lz8knPf1MeAVERkpIlkJDumwnFj/LiLePus7sVj/D3hfRCaKyKmJjulwyvxvvYc9L9+zyryvHwIXeTWhOZ8BjwHvisgEEenubPf856/nAzyYiJwCfAnUA7ofqVZToonIBcBCIB+4GvAZY3YlNqryObWwnnRuPgycBPwpcREdnog0A54DRgMXi0jDBId0WCLSA/jAuXk9UAs7bs8RkQ7AHCDLGNPVGPN1omM6HBEZCswFSoBLgdXA3oQGdXh/BtKA64BS4O8AxhjPVylLmgQhItGTOpcA3wOPGWOCIpItIlki4pkZWWViXQNcYIyZZIz5AGjlVKxFPFJHu8z71hToYoz5jTHmS+ATIFtEbkpcdIdVBPwNe2FlXWCYk+C8aBvwa2PM7caYpUAKsDjBMR3OPmA+8DqAiPQSkc4e7XJcB4wxxvzO+XvtDAwGT/1vWU63UhtgojFmDfAu0ExE/hDdJ5ExHo2ng4NY3/2LwGQRaeX8k30G/EZEPsf+1vsY8FCiu28OjhVYYYxZKSIpIlIDuyncCMAkePpYmVj/4Lyvm4EdInKvs0sa9uLFy0WkUcICBUSkrYj8NnrbGLMHWGqM2Qe8BFxEBUsHxFs5sW5z/gZqisjj2K2HoU73WEJPQVZOrJuA/wC/FJHF2P9XE4GXRaRJgsIEDo0V2GCMWessqgV4FegEif3fKhunMSZijMkFMoDoGEkLYDZwrYg083orwtMJwvnAfxFYjr147kEROQeYgl276T1jzBnAZOf+hA0ClxPrZOBCAGNMqTGmCPvbbnNn/4R94z3M+3oucAtwpYhMAaYCb2N/YGQnKFRE5GLg38CdTh0uRMTnJAeMMe8Ae4AxZVpuCXG4WJ27C7D/XlsC1wJB7L+RhCgvVsf/wy5787Lzv3UjsAm499BnccdhYrXA/t9ybgei2xL1rfwI7+k1QF8ReQX4F/b/1mvAOa4Heaz2n7zEexegHzDHuW5h9zU+j90/Hjho3weAqzwW6z+BbmX2GQ4s8uj7OgU7edUAzgTqOffPBLITGGt/YBBwGnbXTA1nuwCWc72XE2d74DKgvcdi9ZWz7zjg9157X5376h6071XA77wWq/O3G13LNQD7dAEJibEC72lNoAdQ37n9NPAL57okMu4jXTzdgsB+k1NFpKexm2JfYn+bGWuMiZ09VkS6Yn+obUtMmED5sW7GPqdF1CZghYi4e0KGQx3ufb3KGFNkjPnMGLNbRHpjfxAnbPDPGDMP+I8x5r/AMuB+5y7LiR1jzELsc/osxJ4MECznqeLuCLEe0OUhItnY3yq/dzXAMg4Xq4hYxu4WwbndHbgS2JCIOOGI76sY5xPW2b48kTPEjvCe+owx+caYxcaYXc7/1clArvM4z65W9kSCEJGmItKunLsC2KU3RgIYYzYAS4FaIlJXROqLyDTsVsWTxi7V4bVYa4pIPWcfA9TG7m6Iu2OM9Vsn1joiUltE/ordqnjVOCdxSlCslPky8DB2/3gXY0zYeZxPRMZgd4NdY4zpH+94jyPWiPO4RiLyGvYMrCeMMa/HM84TjLWOiEQrJz9hjHnDg7GWPTlpHWA79peyuDqBv9VMEZmIPW42xRizKN6xnrBEN2Gw+8C3Yn8Y/ezgJhfwC+e+c5zb7bAHT9Oc26OTINb0JHxf+yUy1iPs+wfsfy6wTyglQPMkiLW38/OXSRBrH+fnkGSIlYO6m70ap/PzZLdirYyLF1oQ9YEd2NNXe4pIijHGlBloWoA93/mPItIS6Ir9LaEWgDHm7SSINcPFGE801uj7+mUiYz3Cvn8EzhORAuBs7P79LS7EGHW8sZ7rdDV87EaQjuON9Rwn1tluBOk43lh/Ca5ObT3uv1UR8RtjVrkRZGVxfe2AiGQaY8p2seSz/7wPg4Al2FPYov3LedjT7FphL9yKdiXsSKJY436WqKr4vh70GMFe2TvFue83xpgvNFaNVeOMIxebZvWxR+7fwT4ZUPNo0wt7Sh3O/bcCFwN1nG1lu0VqaKzVL9Yyj08BhmmsGqvG6c7FlS4mERmCvbhtO/b830HYi5vAHu3f6FyvBzyC3f+9Bw4c4Tf2WgKNtZrF6jxejL2e5BONVWPVON3hVhfTTuD/jDEvgV2FFfvNBegNXCD2Qq3N2CsiV4tIPWPMbpfi01g9HmvZhKaxaqwapztcSRDGmG9E5Dux51hHgByc1bnGmM9FZB7woTFmptjF+EaRoCm4GqvGqrFW71iTJU43VOoZ5cSuN9QPmGv2L4Evb7+HgO3GmP+rtIMfI401PjTW+NBYK1+yxJlIlZb1RGQCdsndU4HMMtutMtd9zgh/I+zVu4jIENm/kMwVGqvGqrFW71iTJc5Eq5QuJrHn0Y8AhhpjVpfZLmb/ysz2xpg1Ys8bbohdSvqP2AM+CysjDo1VY9VYNdaqEqcXHHeCEHvRR8i5WRfYYYxZLSJtge7AAmPMBhFpDvwfEBC7wmFL7MUtAeDvxph/ndhL0Fg1Vo1VY60acXqOOfb5wX7sqV2P4yzFx+7H+3/AGdgrdJ/BLhM9FnvU/96DnuOWYz3u8Vw0Vo1VY63esSZLnF69HOubLdiny3sVuAT4FLjeue9L7LnDXZzbv8CeBmaVeXyKay9MY9VYNdZqHWuyxOnly7G+4bWAr4Cazu2znF/AEOxa5wWUqcUPvAm0iv6yXH1hGqvGqrFW61iTJU4vX45pFpMxZi/2IM0VzqYvsZto5xtjFgMvA1eJyDkiMhW7ebfFeayri0c0Vo1VY63esSZLnF52PNNc/4U9ot/U2EWsvgWCInIS9nlX/43dl/e9MWaU2T8wlAgaq8aqsVbvWJMlTm861iYH0BR4CJhYZtsXwOllbrtWn11j1Vg1Vo012eP06uWYWxDGmG3Ae9hnS7pQRLKAYqC0zD4JOeXjwTTW+NBY40NjrXzJEqdXHXepDRH5JXAh0Bd4yhjzVGUGVpk01vjQWONDY618yRKn15xQLSYRCWCP53i+305jjQ+NNT401sqXLHF6SaUW61NKKVV1VMkStUoppU6cJgillFLl0gShlFKqXJoglFJKlUsThFJKqXJpglBKKVUuTRBKKaXKpQlCKaVUuf4/Q4iXLrWEuDAAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "pvlib.atmosphere.get_relative_airmass(solpos['zenith']).plot(label='kastenyoung1989, zenith')\n", "pvlib.atmosphere.get_relative_airmass(solpos['apparent_zenith']).plot(label='kastenyoung1989, app. zenith')\n", "pvlib.atmosphere.get_relative_airmass(solpos['zenith'], model='young1994').plot(label='young1994, zenith')\n", "pvlib.atmosphere.get_relative_airmass(solpos['zenith'], model='simple').plot(label='simple, zenith')\n", "plt.legend()\n", "plt.ylabel('Airmass')\n", "plt.ylim(0,100)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEKCAYAAAAFJbKyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XmcFPWd//HXp3suhlsuEYxgRDGoXOOtBAUNRgOe0SS6aMySzaUx8cxugiYxayKbrNnkR0LUQDbGI54Yj6h4kVVREPAADKgIIyDIfcwMc3x+f1QN04xzdDPdXTPd7+fj0Y/q+tbRny6K/sz3W1Xfr7k7IiIiyYpFHYCIiHQsShwiIpISJQ4REUmJEoeIiKREiUNERFKixCEiIinJWOIwszvNbL2ZvZVQtp+ZPW1my8Npz7DczOzXZrbCzN4ws1GZiktERNomkzWOmcCERmXXA3PcfQgwJ5wHOAMYEr6mANMzGJeIiLRBxhKHu78IbGpUPAmYFb6fBZydUP4nD7wC9DCz/pmKTURE9l1Blj+vn7uvBXD3tWbWNywfAKxOWK88LFvbeAdmNoWgVkLnzp1HDx06NLMRi4i0werNu9iyq5qBPTvRs7Qo6nAAWLBgwcfu3mdft8924miONVHWZF8o7j4DmAFQVlbm8+fPz2RcIiJt8p27F/Lo4jX88qIRTBoxIOpwADCzD9qyfbbvqvqovgkqnK4Py8uBAxPWGwisyXJsIiJpV1NbB0BBLHduYs32N5kNTA7fTwYeSSj/l/DuquOArfVNWiIiHVl1bdB4UhBvqmGlY8pYU5WZ3Q2MBXqbWTkwFbgFuM/MLgdWAReEqz8OfB5YAewCLstUXCIi2VRTF9Q4iuK5U+PIWOJw9y81s2hcE+s68K10fG51dTXl5eVUVlamY3c5o6SkhIEDB1JYWBh1KCJ5pUY1jvavvLycrl27MmjQIMxy5x+qLdydjRs3Ul5ezuDBg6MORySvVOsaR/tXWVlJr169lDQSmBm9evVSLUwkAvWJozCHahw5lzgAJY0m6JiIRKOmrr6pKnd+bnPnm4iItEN77qqK5c4fb0ocWfK1r32NJUuWpGVfXbp0Sct+RCTzavY0VeXOz23OXRxvr26//faoQxCRCDQ0VanGIS3YuXMnZ555JsOHD+eII47g3nvvZezYsdR3j9KlSxeuu+46Ro8ezfjx43n11VcZO3YsBx98MLNnzwZg5syZTJo0iQkTJnDYYYdx0003NflZt956K0cffTRHHXUUU6dOzdp3FJHk1F8c13McHcRNj77NkjXb0rrPzxzQjalfGNbiOk8++SQHHHAAjz32GABbt25l+vSGnuJ37tzJ2LFj+fnPf84555zDf/zHf/D000+zZMkSJk+ezMSJEwF49dVXeeuttygtLeXoo4/mzDPPpKysbM9+nnrqKZYvX86rr76KuzNx4kRefPFFxowZk9bvLCL7Lhef48idFNiOHHnkkTzzzDNcd911zJ07l+7du++1vKioiAkTJuxZ97Of/SyFhYUceeSRrFy5cs96p512Gr169aJTp06ce+65/OMf/9hrP0899RRPPfUUI0eOZNSoUSxbtozly5dn/PuJSPLqnxzPpec4crrG0VrNIFMOPfRQFixYwOOPP84NN9zA6aefvtfywsLCPbfHxmIxiouL97yvqanZs17jW2gbz7s7N9xwA1//+tcz8TVEJA121+g5DknCmjVrKC0t5eKLL+bqq6/m9ddf36f9PP3002zatImKigoefvhhTjzxxL2Wf+5zn+POO+9kx44dAHz44YesX7++qV2JSERy8TmOnK5xROXNN9/kmmuuIRaLUVhYyPTp07n66qtT3s9JJ53EJZdcwooVK/jyl7+81/UNgNNPP52lS5dy/PHHA8FF9z//+c/07du3qd2JSARqcvA5Dgv6F+yYmhrIaenSpRx++OERRZQ+M2fOZP78+fzmN79J2z5z5diIdCSDb3gMd3j3Z58n3k6Sh5ktcPey1tdsWu7UnURE2qH6v83bS9JIBzVVtVOXXnopl156adRhiIh8gmocIiKSEiUOERFJiRKHiIikRIlDRERSosSRBTfeeCPTpk1rdvmGDRs49thjGTlyJHPnzk15/zNnzuTb3/42AA8//HDaum8XEWmKEkc7MGfOHIYOHcrChQs5+eST27QvJQ4RyTQljgy5+eabOeywwxg/fjzvvPMOAO+++y4TJkxg9OjRnHzyySxbtoxFixZx7bXX8vjjjzNixAgqKir4xje+QVlZGcOGDdurq/RBgwbx8ccfAzB//nzGjh2712e+9NJLzJ49m2uuuYYRI0bw7rvvZu37ikj+yO3nOJ64Hta9md597n8knHFLi6ssWLCAe+65h4ULF1JTU8OoUaMYPXo0U6ZM4Xe/+x1Dhgxh3rx5fPOb3+TZZ5/lxz/+8V5Pid98883st99+1NbWMm7cON544w2OOuqoVkM74YQTmDhxImeddRbnn39+Wr6uiEhjuZ04IjJ37lzOOeccSktLAZg4cSKVlZW89NJLXHDBBXvWq6qqanL7++67jxkzZlBTU8PatWtZsmRJUolDRCQbcjtxtFIzyKTGXaDX1dXRo0cPFi1a1OJ277//PtOmTeO1116jZ8+eXHrppVRWVgJQUFBAXdi3f32ZiEi26RpHBowZM4aHHnqIiooKtm/fzqOPPkppaSmDBw/mr3/9KxCMpbF48eJPbLtt2zY6d+5M9+7d+eijj3jiiSf2LBs0aBALFiwA4IEHHmjys7t27cr27dsz8K1ERAJKHBkwatQoLrzwQkaMGMF55523506pu+66izvuuIPhw4czbNgwHnnkkU9sO3z4cEaOHMmwYcP46le/utcYHFOnTuXKK6/k5JNPJh6PN/nZF110EbfeeisjR47UxXERyQh1q55HdGxEsm/Q9Y8BsPKWMyOOpIG6VRcRkaxS4hARkZTkZOLoyM1vmaJjIiLpknOJo6SkhI0bN+qHMoG7s3HjRkpKSqIORURyQM49xzFw4EDKy8vZsGFD1KG0KyUlJQwcODDqMEQkB+Rc4igsLGTw4MFRhyEikrNyrqlKREQyK5LEYWZXmdnbZvaWmd1tZiVmNtjM5pnZcjO718yKoohNRERalvXEYWYDgCuAMnc/AogDFwE/B37l7kOAzcDl2Y5NRERaF1VTVQHQycwKgFJgLXAqcH+4fBZwdkSxiYhIC7KeONz9Q2AasIogYWwFFgBb3L0mXK0cGNDU9mY2xczmm9l83TklIpJ9UTRV9QQmAYOBA4DOwBlNrNrkgxjuPsPdy9y9rE+fPpkLVEREmhRFU9V44H133+Du1cCDwAlAj7DpCmAgsCaC2EREpBVRJI5VwHFmVmrBaEfjgCXAc0D9eKeTgU/2OS4iIpGL4hrHPIKL4K8Db4YxzACuA75nZiuAXsAd2Y5NRERaF8mT4+4+FZjaqPg94JgIwhERkRToyXEREUmJEoeIiKREiUNERFKixCEiIilR4hARkZQocYiISEqUOEREJCVKHCIiGeLeZJd7HZ4Sh4hIhlTXBomjMG4RR5JeShwiIhlSWVMLQHFBPOJI0kuJQ0QkQyqrg8RRUphbP7W59W1ERNqRquo6QDUOERFJUlV9U5VqHCIikozKsMZRohqHiIgko77GoWscIiKSlD01jkLVOEREJAkNd1UpcYiISBIaahy59VObW99GRKQdqdIDgCIikgrVOEREJCX11zhU4xARkaRU1eiuKhERSUFDjSO3fmpz69uIiLQjlTW6HVdERFJQpYvjIiKSiirVOEREJBWVe7pVz62f2tz6NiIi7Yi6HBERkZRUhImjkxKHiIgkY1dVmDiKlDhERCQJO3fXANC5qCDiSNJLiUNEJEN2VgWJo0uJEoeIiCRhR9hU1blYTVUiIpKE+hpH1+LCiCNJr0gSh5n1MLP7zWyZmS01s+PNbD8ze9rMlofTnlHEJiKSLhXVGnM8nW4DnnT3ocBwYClwPTDH3YcAc8J5EZEOz8yiDiGtsp44zKwbMAa4A8Ddd7v7FmASMCtcbRZwdrZjExGR1kVR4zgY2AD80cwWmtntZtYZ6OfuawHCad+mNjazKWY238zmb9iwIXtRi4gIEE3iKABGAdPdfSSwkxSapdx9hruXuXtZnz59MhWjiIg0I4rEUQ6Uu/u8cP5+gkTykZn1Bwin6yOITUREWpH1xOHu64DVZnZYWDQOWALMBiaHZZOBR7Idm4iItC6qxxm/A9xlZkXAe8BlBEnsPjO7HFgFXBBRbCIi0oJWE4eZnQgscvedZnYxQbPSbe7+wb5+qLsvAsqaWDRuX/cpIiLZkUxT1XRgl5kNB64FPgD+lNGoRESk3UomcdS4uxM8Z3Gbu98GdM1sWCIi0l4lc41ju5ndAFwMjDGzOJBbHa+IiEjSkqlxXAhUAZeHd0QNAG7NaFQiItJuJVXjIGiiqjWzQ4GhwN2ZDUtERNqrZGocLwLFZjaAoPPBy4CZmQxKRETar2QSh7n7LuBc4H/c/RxgWGbDEhGR9iqpxGFmxwNfAR4Ly3JrOCsREUlaMonjSuAG4CF3f9vMDgaey2xYIiLSXrV6cdzdXyS4zlE//x5wRSaDEhGR9iuZLkf6EDwxPgwoqS9391MzGJeISIdWVRMMG1sQy63R/yC5pqq7gGXAYOAmYCXwWgZjEhHp8LZWVAPQozT3npdOJnH0cvc7gGp3f8Hdvwocl+G4REQ6tC276hNHUcSRpF8yDwBWh9O1ZnYmsAYYmLmQREQ6vk07dwPQMwdrHMkkjp+aWXfg+8D/AN2AqzIalYhIB5fXNQ53/1v4ditwSmbDERHJDVt25XGNw8wGE4zYNyhxfXefmLmwREQ6ts1hjaNnPtY4gIeBO4BHgbrMhiMikhvqaxx52VQFVLr7rzMeiYhIDsn3i+O3mdlU4CmCcTkAcPfXMxaViEgHtyV8jqNn5/yscRwJXAKcSkNTlYfzIiLShIaL4/mZOM4BDnb33ZkORkQkV2zeld9Pji8GemQ6EBGRXJLvNY5+wDIze429r3HodlwRkWbUXxzPxRpHMoljasajEBHJMXUeTAvjyTTsdCwtJg4ziwM/dPfxWYpHRETauRZTobvXArvCvqpERESSewAQeNPMngZ21he6u0YBFBHJQ8kkjsfCl4iISFK9487KRiAiItIxNJs4zOw+d/+imb1J8KR4Inf34ZkNTURE2qOWahxXhtOzEsqMYPS/H2QsIhERadeaTRzuvjacfmBmI4AvA18E3gceyE54IiLS3rTUVHUocBHwJWAjcC9g7q5RAEVE8lhLTVXLgLnAF9x9BYCZaaxxEZFW1NY1viycW1p6APA8YB3wnJn9wczGEVzjSAszi5vZQjP7Wzg/2MzmmdlyM7vXzHKvZzARyQsbdwTd+vXuUhxxJJnRbOJw94fc/UJgKPA8cBXQz8ymm9npafjsK4GlCfM/B37l7kOAzcDlafgMEZGs+2hbkDj6dcuzxFHP3Xe6+13ufhbBHVWLgOvb8qFmNhA4E7g9nDeCgaHuD1eZBZzdls8QEYnK+u2VAPTtmqeJI5G7b3L337t7W0f/+2/gWhpGFOwFbHH3mnC+HBjQ1IZmNsXM5pvZ/A0bNrQxDBGR9GuocZREHElmZL2/XzM7C1jv7gsSi5tYtcmrS+4+w93L3L2sT58+GYlRRKQt1m0Laxw5mjiS6asq3U4EJprZ54ESoBtBDaSHmRWEtY6BwJoIYhMRabPyzbsAGNijU8SRZEbWaxzufoO7D3T3QQTPiTzr7l8BngPOD1ebDDyS7dhERNKhfFMFAAP3U+LItOuA75nZCoJrHndEHI+IyD5ZHdY4PrVfacSRZEYUTVV7uPvzBLf64u7vAcdEGY+ISDqs3Rpc4+jfXTUOERFJQTyWtmem2xUlDhGRbHrhF7DqlaijaBMlDhGRbKmthuduhjs/F3UkbaLEISKSLatfDaYDRkcbRxspcYiIZMuKZ4Lpp9va+Ua0lDhERLLl3TnB9JDx0cbRRkocIiJpVNfSWBxrFwfTAWXZCSZDlDhERNLowy3BU+MDWupuJB7pI3RtpsQhIpJGy9dvB+DTfbtEHEnmKHGIiKTR8o92ADCkceKo2R1BNJmhxCEikkbL1zeTOFa9HEw7+PUNUOIQEUmrFWHiOKRx4vjn34PpoR374T9Q4hARSatmE8eyR4Pp0LOyHFH6KXGIiKTRjqpgBOwepUV7L9iyKpj2PTzLEaWfEoeISDZZx+8xV4lDRCTTvIWHAjsgJQ4RkUxbvySY9jgo2jjSRIlDRCRNqmvrml7w1oPB9IhzsxdMBilxiIikSf3Df5+4o+qtB4LpMCUOERFJsHD1ZgBGHNhj7wWb3w+m+x+Z5YgyQ4lDRCRNFqwMEkfZQT2bXiEH7qgCJQ4RkbSZ/0GYOAYlJI7qyoiiyRwlDhGRNFm1aRcAB/dOuMbxzmPBtIOP+pdIiUNEJM1isYQmqQUzg+noyyKJJROUOEREMun9F4PpYWdEG0caKXGIiKSBt/Z0eLwwO4FkgRKHiEgarNwYXN/Yv1tJQ+HunRFFk1lKHCIiaTBn6UcAnHp434bCN+4Npp+ZFEFEmaPEISKSBnOWrgdgfGLiePm3wfTYb0QQUeYocYiIpMHL720E4IRP924o3LgimH7quAgiyhwlDhGRNCopjAdvEi+W58gT4/WUOEREMuH9F4Jpn44/4l9jShwiIm20YXsVAEXxhJ/U534WTE+5IYKIMkuJQ0SkjR5Z9CEAk0Yc0FC4el4wHfqFCCLKrKwnDjM70MyeM7OlZva2mV0Zlu9nZk+b2fJw2kz3kiIi7cu9r60G4KJjDgwK6mobFsZy7+/zKL5RDfB9dz8cOA74lpl9BrgemOPuQ4A54byISLu3fH0wgNOoT4V/7y6+J5gOHhNRRJmV9cTh7mvd/fXw/XZgKTAAmATMClebBZyd7dhERNrC6u+eeuLaYHr6zdEFk0GR1qHMbBAwEpgH9HP3tRAkF6BvM9tMMbP5ZjZ/w4YN2QpVRKRJa7ZUABBP7BF3d1ADof9REUSUeZElDjPrAjwAfNfdtyW7nbvPcPcydy/r06dP5gIUEUnCrJdWAvC1kwYHBeveDKYWjyagLIgkcZhZIUHSuMvdHwyLPzKz/uHy/sD6KGITEUnF7198D4BLTxwUFDzyrWB67oxoAsqCKO6qMuAOYKm7/zJh0Wxgcvh+MvBItmMTEdlX/bt3Ct6sXRxMjzgvumAyrCCCzzwRuAR408wWhWU/AG4B7jOzy4FVwAURxCYikrSPtjUaT3ztGw3vc6ybkURZTxzu/g+guSM6LpuxiIi0xS+efAeA744fEhT85YvB9MK7IoooO3LvyRQRkSx54PVyAL4+5tNBp4bb1wYLDj8rwqgyT4lDRGQfVNfW7XnfqSgOr0wPZg4YGVFE2aPEISKyD6Y//y4AE4btHxT8PezM8KK7I4ooe5Q4RET2wS+f/icAPz3nCFizqGFBt/4RRZQ9ShwiIinasmv3nve9uxTDjM8GMxc/EFFE2aXEISKSou/cvRCAK8YNgU3vNSw4ZHxEEWVXFM9xiIh0WO7O3OUfA3DluCHw4x7BgnNvjzCq7FKNQ0QkBfVdjAzp24X46pcbFhyVP88sq8YhIpKCW55YBsA9U46DaWEn3pc9GWFE2acah4hIkn73QnALbvdOhfSa94ug0GJw0PERRpV9ShwiIkmorfM9tY3nvj4U5k4LFlzzboRRRUOJQ0QkCef/7iUAjj+4F/v97sigcNyPoHS/CKOKhhKHiEgrFq/ewsJVWwD4S5dfNSw4+fsRRRQtJQ4RkRZU19Yx6bf/B8DfT1qO/TO8EP6DtRFGFS0lDhGRFgz59ycAuHj/1Rw2f2pQ+I2Xoag0wqiipcQhItKML/4+eE5jlP2Tn265Lig8/4/Q7zMRRhU9JQ4RkSZc89fFvPr+Jk6IvcWDxTcGhWfcCkecG2lc7YEeABQRaeTf/ncBT769jsvjj/PDwj8HhZ+fBsf8a7SBtRNKHCIiIXen7KfPsHFnFX8r+neOiK0MFlx0Nwz9fKSxtSdKHCIiwJotFZxwy7McaqtZUHJdw4JvL4Deh0QXWDukxCEiee8/n1jKH15Ywb1FP+XYWPB0OJ37wveWQlw/k43piIhI3nrlvY1cNOMVvhl/mPdK7mtYcMEsGHZ2dIG1c0ocIpJ33vpwK2f9z1yuKriflSUPNSwYPAYueRhi8eiC6wCUOEQkb8xevIYb736B24umsbJkRcOCrv3hGy/lZb9T+0KJQ0Ry2potFVzx53mcsu4OvlUwm4klCQsPGAWXPASdekQWX0ekxCEiOeedddu55f4X+Ny6P3BRwfPcD3v/2p36w6CDQrOIIuzYlDhEpMP7eEcVf3p2EbXzbueKggc5zGr4I+z1C+fDzsXO/C81R6WBEoeIdCjVtXU8++Zqlr5wL8M//hunxBfTG/geQOHe6/rJV2MnfReKu6K6RfoocYhIu1Rb5yx6fx3vzHuSwhV/Z2ztS/SxrRQCnwtfNLr5qeLIS+g09iro9WkAJYsMUeIQkcis3byDd96cz7blc+m0dj5HVL9Bf9sEBDlhdPgCPpEFNvQ9kZ4nTKbgiLOhoBiATtkKPM8pcYhIWtXU1rF63QbWrljEjtVvEN+whB7bl3NQ7Qf0tm17rds/fO3RRBXhw24jKD58Ar2P+eKemgRAn4xEL8lQ4hCRT6itrWPj5k1sXPMeOz5aSeXGVdRtLadg+xq6VHxIn5q1HGAbm9y2ABgcvvbSTLvRusID2dp7FF2GnEj/o04j1mvwXnc7DUjHF5K0UuIQ6YDcnYqq3WzbupmdWz9m19aPqdy2keodG6nbuRF2bSResZHCyk2UVG+mtGYL3eq20NO3ETdvdf9xoG/4alIrFw+qKGJd0UFs734osf2H0f2gEfQ7ZAQF3Q/4xC2w+4cv6TiUOEQILsTurq5l9+4qdlfuYnflTqqrdlFTuYOaql3UVFVQt3sndVW7qKuuoK5qJ747eE91BVZdgdUEr1hNBbHaKgpqKyioq6SwroqiukqKvYpOVNKZSmJJ/Hi3xIDS8JXyhkmqoJiN8T5sL+pHZWl/vNsAinodSNf9D6H3gYfRufenmu0AsBg4KNXYpMNQ4kjRB+8sYtMHbwEOdbW4O3gdXlcH1EFdHXgNXufgteDgdbW414KHZXV1e7bDaxPe128fvGzPNnU4Hs7XgTtWvz51WF1QHvNaCNcLXnXBi2BfsXCbGMG6wbKg3MJy89pw2njeiVEbvNyDfex51RIP10/mr9n2KE5wYTWjF1czcIvPdkrZaV3YFe9KZWF3qgu7UVPck9pOvbDOvSns2ofi7n3p3LMfXXv1p9t+/YgXFie1707AwPSHLDmgXSUOM5sA3Ebw//h2d78l4pD2snnDWg66+7P6S6opOXLfYzVxdlNIFcVUWTG7rYTqWBHVsRKqYyXUxIqpjZdQG++EF3TCC4OXFZRAYSnx4uBVWNyZwk6dKSoppahTFzqVdqW4tCudunSjsLhL2rrq7hq+RLKp3SQOM4sDvwVOA8qB18xstrsviTayBj169WPe4TdQ9MHzQAy3GI6BxXBiuCW+j4VtubGgLBaD+nUtKKtfl1jDvBFst2edWBwL5y1hu2BZuH4sWA+LY7F4sG2sAIvFMYvhCcssXEYshsUKsHBbixUEy+LB1GJxYvE4sXgBsVgwjcfjxAsKiMcLKYjHgmlhUF5QUEhBPE5BPE48Xh9rx8smheGrc9SBiLRj7SZxAMcAK9z9PQAzuweYBLSbxGGxGMdeeD1wfdShiIhEpj0ljgHA6oT5cuDYxiuZ2RRgSjhbZWZvZSG2jqA38HHUQbQTOhYNdCwa6Fg0OKwtG7enxNFUu8YnrrS6+wxgBoCZzXf3skwH1hHoWDTQsWigY9FAx6KBmc1vy/axdAWSBuXAgQnzA4E1EcUiIiLNaE+J4zVgiJkNNrMi4CJgdsQxiYhII+2mqcrda8zs28DfCW7HvdPd325lsxmZj6zD0LFooGPRQMeigY5FgzYdC3PvmA9siYhINNpTU5WIiHQAShwiIpKSDps4zGyCmb1jZivMLG+eyDOzA83sOTNbamZvm9mVYfl+Zva0mS0Ppz2jjjVbzCxuZgvN7G/h/GAzmxcei3vDmy1ynpn1MLP7zWxZeH4cn6/nhZldFf7/eMvM7jazknw6L8zsTjNbn/icW3PnggV+Hf6WvmFmo1rbf4dMHAndk5wBfAb4kpl9JtqosqYG+L67Hw4cB3wr/O7XA3PcfQgwh/x6vP1KYGnC/M+BX4XHYjNweSRRZd9twJPuPhQYTnBM8u68MLMBwBVAmbsfQXCzzUXk13kxE5jQqKy5c+EMYEj4mgJMb23nHTJxkNA9ibvvBuq7J8l57r7W3V8P328n+HEYQPD9Z4WrzQLOjibC7DKzgcCZwO3hvAGnAveHq+TFsTCzbsAY4A4Ad9/t7lvI0/OC4I7RTmZWQND7/Fry6Lxw9xeBTY2KmzsXJgF/8sArQA8z608LOmriaKp7krwbKMzMBgEjgXlAP3dfC0FyoYUxeHLMfwPXAnXhfC9gi7vXhPP5cm4cDGwA/hg2291uZp3Jw/PC3T8EpgGrCBLGVmAB+XleJGruXEj597SjJo6kuifJZWbWBXgA+K67b2tt/VxkZmcB6919QWJxE6vmw7lRAIwCprv7SGAnedAs1ZSw7X4Swei1BxB0dnxGE6vmw3mRjJT/z3TUxJHX3ZOYWSFB0rjL3R8Miz+qr16G0/VRxZdFJwITzWwlQXPlqQQ1kB5hEwXkz7lRDpS7+7xw/n6CRJKP58V44H133+Du1cCDwAnk53mRqLlzIeXf046aOPK2e5KwDf8OYKm7/zJh0Wxgcvh+MvBItmPLNne/wd0HuvsggnPgWXf/CvAccH64Wr4ci3XAajOr7/V0HMGQBHl3XhA0UR1nZqXh/5f6Y5F350UjzZ0Ls4F/Ce+uOg7YWt+k1ZwO++S4mX2e4K/L+u5Jbo44pKwws5OAucCbNLTr/4DgOsd9wKcI/uNc4O6NL47lLDMbC1zt7meZ2cEENZD9gIXAxe5eFWV82WB3MGCSAAAEbUlEQVRmIwhuEigC3gMuI/jjMO/OCzO7CbiQ4C7EhcDXCNrt8+K8MLO7gbEEXcl/BEwFHqaJcyFMrr8huAtrF3CZu7fYe26HTRwiIhKNjtpUJSIiEVHiEBGRlChxiIhISpQ4REQkJUocIiKSEiUOadfM7BwzW9ToVWdmTT0JnMz+fmxm48P33zWz0oRlO9IVd5Kx3GhmV6e4zdlm9qNmlu1z/GY2zcxO3dftJb/odlzpUMxsCvAV4BR3r2tt/Vb2tZKgB9WPw/kd7t6l7VEm/fk3AjvcfVoK27wETKyPudGyfY7fzA4C/uDup+/L9pJfVOOQDsPMDgV+BFxSnzTM7Bozey0cR+CmsGxQOB7FH8IxGZ4ys07hsplmdr6ZXUHQj9FzZvZcwmfcbGaLzewVM+vXRAzHmNlLYUeCL9U/qW1ml5rZg2b2ZDjewS8StrnczP5pZs+HMf2mif1+Otx2gZnNNbOhzXz/qoREN9jMXg6//08arfuJ4xKW/9CC8TqetmCciqsB3P0DoJeZ7Z/8v4jkKyUO6RDC/rn+QvB0+Kqw7HSCMQSOAUYAo81sTLjJEOC37j4M2AKcl7g/d/81QX88p7j7KWFxZ+AVdx8OvAj8axOhLAPGhB0J/gj4WcKyEQRPKx8JXGjBoFsHAD8kGDvlNOATCSE0A/iOu48Grgb+XxPrnAi8njB/G0GnhkcD6+oLmzsuZlYWHoeRwLlAWaP9vx5+hkiLClpfRaRd+Anwtrvfk1B2evhaGM53IfjBXEXQyd2isHwBMCiJz9gN/C1hm9OaWKc7MMvMhhD0IFqYsGyOu28FMLMlwEEEXT68UN/Nh5n9FTg0cYcW9HR8AvDXoPcHAIqb+Oz+BF2n1zuRhoT4vwQDFUHzx6Ur8Ii7V4Sf+2ij/a8nqIWJtEiJQ9q9sB+q8wh6e91rEfCf7v77RusPAhL7IKoFOiXxUdXecNGvlqb/f/wEeM7dzwk/5/mEZY0/s4Cmu6xuLEYwVsSIVtarIEhciZq6SNnccbmqlf2XhJ8h0iI1VUm7ZsHYCn8E/iUc8TDR34Gvhn+xY2YDzCyVgYq2E/wVnoruwIfh+0uTWP9V4LNm1tOCLr3Pa7xCOJ7K+2Z2AewZA3p4E/taChySMP9/BL0CQ3DDQL3mjss/gC9YMP52F4KRExMdCryFSCuUOKS9+zeCkcqmN7ol90J3f4rgusfLZvYmwRgUqSSCGcATiRfHk/AL4D/N7P8IemZuUTga3c8Iei9+hqB7761NrPoV4HIzWwy8TdNDIb8IjLSG9qwrCcacf42Emkhzx8XdXyPoQnsxwRgV8+tjCa8hHRKWibRIt+OKZJiZdXH3HWGN4yGCYQAe2sd93QY86u7PtDGWUoJENMXdXzezc4BR7v7Dfdmv5BfVOEQy70YzW0TQDPQ+wbgI++pnQGmrazVvRhjL68AD7l5/l1YB8F9t2K/kEdU4REQkJapxiIhISpQ4REQkJUocIiKSEiUOERFJiRKHiIik5P8DOEkfaD75UBwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(solpos['zenith'], pvlib.atmosphere.get_relative_airmass(solpos['zenith'], model='simple'), label='simple')\n", "plt.plot(solpos['zenith'], pvlib.atmosphere.get_relative_airmass(solpos['apparent_zenith']), label='default')\n", "plt.xlim(0,100)\n", "plt.ylim(0,100)\n", "plt.xlabel('Zenith angle (deg)')\n", "plt.ylabel('Airmass')\n", "plt.legend()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.5" } }, "nbformat": 4, "nbformat_minor": 1 }