{ "cells": [ { "cell_type": "markdown", "id": "e7c67773", "metadata": {}, "source": [ "(bolometric)=\n", "\n", "# Apply Bolometric Scaling" ] }, { "cell_type": "markdown", "id": "394c5dde", "metadata": {}, "source": [ "This notebook shows how to obtain the bolometric scaling of a light curve to calculate the bolometric luminosity from some observed photometry. First, import the required functions." ] }, { "cell_type": "code", "execution_count": 1, "id": "95166051", "metadata": {}, "outputs": [], "source": [ "from slsne.lcurve import get_bolcorr\n", "from slsne.utils import calc_flux_lum, get_lc\n", "import numpy as np" ] }, { "cell_type": "markdown", "id": "0c0a8912", "metadata": {}, "source": [ "Lets import a test light curve from the reference data directory." ] }, { "cell_type": "code", "execution_count": 2, "id": "29f351cd", "metadata": {}, "outputs": [], "source": [ "# Get a light curve\n", "phot = get_lc('2018lfe')\n", "phot = phot[phot['UL'] == 'False']\n", "# Define the redshift and date of peak\n", "redshift = 0.35\n", "peak = 58468.55" ] }, { "cell_type": "markdown", "id": "55ccc6fd", "metadata": {}, "source": [ "The `get_bolcorr` function will calculate the corresponding bolometric correction factor for all photometry in the table." ] }, { "cell_type": "code", "execution_count": 3, "id": "ca03c446", "metadata": {}, "outputs": [], "source": [ "bol_scaling = get_bolcorr(phot, redshift, peak)" ] }, { "cell_type": "markdown", "id": "a4b57435", "metadata": {}, "source": [ "Before plotting the results, we'll use a couple of built in utilities to calcualte $L_\\lambda$ and import the bolometric light curve." ] }, { "cell_type": "code", "execution_count": 4, "id": "ecba7546", "metadata": {}, "outputs": [], "source": [ "# Calculate F_lambda and L_lambda\n", "F_lambda, L_lambda, lambda_AA = calc_flux_lum(phot, redshift, return_lambda=True)\n", "\n", "# Read in bolometric data from database\n", "bolo = get_lc('2018lfe', 'bol')\n", "bolo['Phase'] = (bolo['MJD'] - peak) / (1 + redshift)" ] }, { "cell_type": "markdown", "id": "dbc9d0be", "metadata": {}, "source": [ "And finally plot the data" ] }, { "cell_type": "code", "execution_count": 5, "id": "e3433f32", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAskAAAG4CAYAAAC3o3oUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC7sElEQVR4nOzdd3hb5fn/8ffRtmR573go0yN7h+mEFaDwIwFKC2WvlgZSdmi/ZRRKGaWMlBQoZZdRWkigtIVSsA2EkeWEJB5Z3vEesi1Zssb5/aFY8cq249i5X9eVy9bR0TmPHII/fnw/96OoqqoihBBCCCGECNIM9QCEEEIIIYQ41khIFkIIIYQQohcJyUIIIYQQQvQiIVkIIYQQQoheJCQLIYQQQgjRi4RkIYQQQgghepGQLIQQQgghRC8SkoUQQgghhOhFQrIQQgghhBC9SEgWQgghhBCiFwnJQgghhBBC9DJiQ7LT6SQtLY0777wzeMxmszFlyhSmTZvGggULhnB0QgghhBDiWKYb6gEMlocffph58+b1Of71118TGho6BCMSQgghhBDDxYicSd6+fTtFRUWcc845Qz0UIYQQQggxDB1zIfmLL77g/PPPJykpCUVRWLVqVZ9zVqxYgc1mw2QyMXfuXNasWdPj+TvvvJNHHnmkz+sURSE7O5vZs2fz5ptvDtZbEEIIIYQQw9wxV27hcDiYOnUq1157LRdeeGGf5//2t79x++238/zzzzN37lyefvppFi5cSHFxMXFxcXzwwQdMmDCBCRMm8PXXX/d47VdffcWoUaOorq7mjDPOYPLkyUyZMqXfcbjdbtxud/Cx3++nqamJ6OhoFEUZ2DcthBBCiEGhqiptbW0kJSWh0Rxzc4PiWKYewwB15cqVPY7NmTNHXbJkSfCxz+dTk5KS1EceeURVVVW955571OTkZDUtLU2Njo5Ww8LC1N/85jd9rn3nnXeqr7zyyj7vff/996uA/JE/8kf+yB/5I39GwJ+KiooBySbi+KGoqqpyjFIUhZUrV7Jo0SIAOjs7MZvN/OMf/wgeA7jqqqtoaWnhgw8+6PH6V199lS1btvDEE0/gcDjw+/1YrVba29vJzs7m+eefZ/bs2f3eu/dMst1uJzU1lYqKCsLCwgb8vQohhBBi4LW2tpKSkkJLSwvh4eFDPRwxjBxz5Rb709DQgM/nIz4+vsfx+Ph4ioqK9vva2tpaFi9eDIDP5+OGG27YZ0AGMBqNGI3GPsfDwsIkJAshhBDDjJRKikM1rELyobr66quDn48ZM4ZNmzYN3WCEEEIIIcSwMawq2GNiYtBqtdTW1vY4XltbS0JCwhCNSgghhBBCjDTDKiQbDAZmzpzJZ599Fjzm9/v57LPPOOGEE4ZwZEIIIYQQYiQ55sot2tvb2bFjR/BxSUkJGzduJCoqitTUVG6//XauuuoqZs2axZw5c3j66adxOBxcc801gzKeFStWsGLFCnw+36BcXwghhBBDz+fz4fF4hnoYYhDp9Xq0Wu1Bn3/MdbfIzc1lwYIFfY5fddVVvPrqqwA8++yz/P73v6empoZp06axfPly5s6dO6jjam1tJTw8HLvdLgv3hBBCiGHiQN+/VVWlpqaGlpaWoz84cdRFRESQkJBwUAs5j7mQfKySkCyEEEIMPwf6/l1dXU1LSwtxcXGYzWbpgjFCqaqK0+mkrq6OiIgIEhMTD/iaY67cQgghhBDiaPD5fMGAHB0dPdTDEYMsJCQEgLq6OuLi4g5YejGsFu4JIYQQQgyUrhpks9k8xCMRR0vX3/XB1J9LSBZCCCHEcU1KLI4fh/J3LSH5AFasWEFWVtZ+d+cTQgghhBAji4TkA1iyZAkFBQWsXbt2qIcihBBCCDFgrr76ahYtWjTUwzhkiqKwatWqQb+PLNwTQgghhOhm/fr1R/V+M2fOPOTXXH311bz22mvBx1FRUcyePZvHH3+cKVOmDOTwBl1X+9/m5mYiIiIOeH51dTWRkZGDPi6ZSRZCCCGEGIbOPvtsqqurqa6u5rPPPkOn03HeeecN9bAGTWdnJwAJCQkYjcZBv5+EZCGEEEKIYchoNJKQkEBCQgLTpk3jnnvuoaKigvr6egA2b97MaaedRkhICNHR0dx44420t7fv83put5ulS5cSFxeHyWTi5JNP7lFumpubi6IofPLJJ0yfPp2QkBBOO+006urq+M9//kNmZiZhYWFcdtllOJ3O4Ov8fj+PPPIIo0ePJiQkhKlTp/KPf/wDgNLS0uAmcpGRkSiKwtVXXw3A/Pnzufnmm7n11luJiYlh4cKFQN9yi8rKSi699FKioqKwWCzMmjWL77777oi/vlJuIYQQQggxzLW3t/PXv/6VcePGER0djcPhYOHChZxwwgmsXbuWuro6rr/+em6++ebgDsa93X333bz33nu89tprpKWl8fjjj7Nw4UJ27NhBVFRU8LwHHniAZ599FrPZzCWXXMIll1yC0Wjkrbfeor29ncWLF/PHP/6RZcuWAfDII4/w17/+leeff57x48fzxRdfcPnllxMbG8vJJ5/Me++9x0UXXURxcTFhYWHBfsYAr732GjfddBOrV6/e5/vOzs5m1KhRfPjhhyQkJLBhwwb8fv8Rf00lJAshhBBCDEMfffQRoaGhADgcDhITE/noo4/QaDS89dZbuFwuXn/9dSwWCwDPPvss559/Po899hjx8fE9ruVwOHjuued49dVXOeeccwB48cUX+fTTT3nppZe46667guf+9re/5aSTTgLguuuu45e//CU7d+5kzJgxAFx88cXk5OSwbNky3G43v/vd7/jf//7HCSecAMCYMWP46quveOGFF8jOzg4G8Li4uD41yePHj+fxxx/f59fgrbfeor6+nrVr1wavM27cuMP6evYmIfkAVqxYwYoVK/D5fEM9FCGEEEKIoAULFvDcc88B0NzczJ/+9CfOOecc1qxZQ2FhIVOnTg0GZICTTjoJv99PcXFxn5C8c+dOPB5PMPwC6PV65syZQ2FhYY9zuy8MjI+Px2w2BwNy17E1a9YAsGPHDpxOJ2eeeWaPa3R2djJ9+vQDvscDLWrcuHEj06dP7zHTPVAkJB/AkiVLWLJkSXDvdyGEEEKIY4HFYukxa/qXv/yF8PBwXnzxxUG9r16vD36uKEqPx13Husodumqg//WvfzFq1Kge5x3M4rvuIb8/3UszBpos3BNCCCGEGAEURUGj0dDR0UFmZiabNm3C4XAEn1+9ejUajYb09PQ+rx07diwGg6FH7a/H42Ht2rVkZWUd9piysrIwGo2Ul5czbty4Hn9SUlIAMBgMAIf1W/spU6awceNGmpqaDnuM+yIhWQghhBBiGHK73dTU1FBTU0NhYSG33HIL7e3tnH/++fzkJz/BZDJx1VVXsWXLFnJycrjlllu44oor+pRaQGDG9qabbuKuu+7i448/pqCggBtuuAGn08l111132GO0Wq3ceeed3Hbbbbz22mvs3LmTDRs28Mc//jHY5zktLQ1FUfjoo4+or6/fbweO3i699FISEhJYtGgRq1evZteuXbz33nt88803hz3mLlJuIYLa2tpwOp14vV4URcFgMGA2mwkJCZF97YUQQohjzMcff0xiYiIQCKMZGRn8/e9/Z/78+QB88skn/OIXv2D27NmYzWYuuuginnzyyX1e79FHH8Xv93PFFVfQ1tbGrFmz+OSTT454446HHnqI2NhYHnnkEXbt2kVERAQzZszgV7/6FQCjRo3iN7/5Dffccw/XXHMNV1555T47cPRmMBj473//yx133MG5556L1+slKyuLFStWHNGYARRVVdUjvspxoKsm2W63ExYWNtTDGVBNTU1UV1fjcrn6fV6n0xEREUFkZCRWq1UCsxBCiGFjf9+/XS4XJSUljB49GpPJNEQjFEfTofydy0zycczlcrFr1y46Ojr2e57X66WhoYGGhgYMBgNJSUlER0cfpVEKIYQQQhx9UpN8nGpvb6eoqOiAAbm3zs5OSktL2bZtW3B7SCGEEEKIkUZC8gGsWLGCrKwsZs+ePdRDGTAej4ddu3YdUe/ntrY2CgsLaWtrG8CRCSGEEEIcGyQkH8CSJUsoKCjosXf5cFdSUoLH4zni63i9XrZv305dXd0AjEoIIYQQ4tghIfk409zcPKCzv6qqUlFRQVVV1YBdUwghhBBiqElIPo6oqsru3bsH5do1NTVUVlYOyrWFEEIIIY42CcnHkaampn22eRsItbW1MqMshBBCiBFBQvJxwu/3U11dPej36dr5RwghhBBiOJOQfJyorq7G7XYflXtVVVUNyh7qQgghhBBHi4Tk44DP5zvqHShKS0txOBxH9Z5CCCGEEANFQvJxoKGhAb/ff1TvqaoqO3fulA1HhBBCjHh+P5SWwubNgY9H+VuuGCSyLfVxoL6+fkju6/F42LlzJxkZGSiKMiRjEEIIIQZTYSGsXAlFReBygckEGRmweDFkZg716MSRkJnkAxjuO+61trYeVi2yqqo0NDRQUlJCaWnpYe/O53Q6KS8vP6zXCiGEEMeywkJYvhzy8yEmBtLTAx/z8wPHCwsH795+v5/HH3+ccePGYTQaSU1N5eGHH+a0007j5ptv7nFufX09BoOBzz77DIA33niDWbNmYbVaSUhI4LLLLutRlpmbm4uiKHz22WfMmjULs9nMiSeeSHFx8eC9oWOQoqqqOtSDGA5aW1sJDw/HbrcTFhY21MM5aDt37qSlpeWgzvV6vaxZs4b//e9/rF69msbGxuBzJpOJ9PR0Jk6cyBlnnMHkyZMPaXZ49OjRREVFHerwhRBCiCOyv+/fLpeLkpISRo8ejclkOqTr+v3w6KOBQJyVBd2/JaoqFBTAjBmwbBloBmFKctmyZbz44os89dRTnHzyyVRXV1NUVITZbObmm2+muroao9EIwFNPPcXy5cvZtWsXiqLw8ssvk5iYSHp6OnV1ddx+++1ERETw73//GwiE5AULFjB37lwee+wxYmNj+dnPfobP52P16tUD/2aOokP5O5eQfJCGY0ju7Oxk8+bNBzzP7/fz1ltv8corr2C320lNTWX+/PlMnjyZiIgIPB4PxcXFFBYWsnHjRurq6pg4cSI33HADJ5100kGFZa1WS2ZmZvAfrBBCCHE0DFZILi2F++4LzBz3FwvsdmhshAcfBJvt8Mffn7a2NmJjY3n22We5/vrrezzncrlISkri+eef55JLLgFg6tSpXHjhhdx///39Xm/dunXMnj2btrY2QkNDgyH5f//7H6effjoA//73v/nBD35AR0fHIX+tjiWH8ncuNckj2MHUItfW1vLAAw+wdu1afvjDH7Jo0SImTJjQJ/jOmTMHCATqb775hldeeYVbb72VU045hbvuuoukpKT93sfn81FSUkJ6errUJwshhBj22toCNcgWS//PWyxQVRU4b6AVFhbidruDAbY7k8nEFVdcwcsvv8wll1zChg0b2LJlCx9++GHwnPXr1/PAAw+wadMmmpubg4v7y8vLycrKCp43ZcqU4OeJiYkA1NXVkZqaOvBv6hgkNckjWHNz836fX7NmDZdeeillZWU899xzLFu27IAhVqPRcNJJJ/Hiiy/y+OOPs23bNn74wx/y2muvHbCDhsPhOCobmgghhBCDzWoNLNLbV7dThyPwvNU68PcOCQnZ7/PXX389n376KZWVlbzyyiucdtpppKWl7RmXg4ULFxIWFsabb77J2rVrWblyJUCfjlR6vT74eVc2ONrdsoaShOQRyu1273fB3saNG7n99tvJysri7bffPuSFiYqicNppp/Huu+9y0UUX8eyzz3Lbbbdht9v3+7qamhra29sP6V5CCCHEsSY1NdDFoqIiUIPcnapCZWWgu8VgTLqOHz+ekJCQ4EK83iZPnsysWbN48cUXeeutt7j22muDzxUVFdHY2Mijjz7KKaecQkZGxlHfS2G4kJA8QrW2tu7zuZqaGu644w6ysrJ44oknCA8PP+z7WCwWbr/9dp555hm2bNnCFVdcQVFR0T7PV1WV0tLS4+onUSGEECOPRhNo8xYTE1ikZ7eD1xv4WFAQOL5o0eAs2jOZTCxbtoy7776b119/nZ07d/Ltt9/y0ksvBc+5/vrrefTRR1FVlcWLFwePp6amYjAY+OMf/8iuXbv48MMPeeihhwZ+kCOAhOQRal8h2ev18qtf/QqTycTvf//7ASu+P/HEE3njjTcIDw/nuuuuIzc3FwC/6me3czc7Wnew27kbv+rH7XZL2YUQQohhLzMTli6F6dMDi/S2bQt8nDEjcHww+yTfe++93HHHHdx3331kZmbyox/9qMeM8KWXXopOp+PSSy/t8b0+NjaWV199lb///e9kZWXx6KOP8sQTTwzeQIcx6W5xkIZTdwtVVdm0aVO/vY1feOEFXn75ZV588cUeBfkDxe12c++995KXl8ctD95C55hOSttLcfvcGLVGbKE2FiQsYEzYGLKysob1ClkhhBDHvsHqbtGd3w/l5YFFelZroMRiMGaQD0VpaSljx45l7dq1zJgxY2gHcwyR7hbHOafT2W9ALi0t5ZVXXuHqq68elIAMYDQa+d3vfsddj93F0988TZY3i6lpUwkJCcHpcbKxcSPF9mIuTLsQs8VMRnrGoIxDCCGEOFo0moFv83a4PB4PjY2N/PrXv2bevHkSkI+AhOQRaF+lFitWrCAuLq5HAf9g0Gg1TPrhJLZ/u52C3ALisuOISYuhtL2UFncLdo+dCkcFW1q2cKXuSuaNnTeo4xFCCCGOF6tXr2bBggVMmDCBf/zjH0M9nGFNQvII1F9ILioqIicnh/vuu29QN/Tw+2FTWSNbaiqZPeF0ilxryV2TSxJJGC1GLHoLBo0Bp8/JpqZNPLn6Se4PvZ+J8RMHbUxCCCHE8WL+/PlIJe3AkIV7B7BixQqysrIOuUXaUPH7/Tj6adr45ptvMmrUKM4999xBu3dJiYlXX03gtXcsFO3UULApnrCwcwgfl8ju+t3oOnUYNAYMWgOKopBgSqCpo4k31r2BX5VuF0IIIYQ4dkhIPoAlS5ZQUFDA2rVrh3ooB6Wtra3PT5B2u53PPvuMCy+8EJ1ucH55UFJi4u23YykuNhNpNhEWokMX0k5tiwt/SAQmJYJdJbtwu914/V60ihaD1kCcKY4t1VvYXrd9UMYlhBBCCHE4JCSPMP2VWnzyySf4fD7OO++8Qbmn3w85ORHY7XpGj+4g3pRIpG8cnaYqLBFOfH4/YaHj0Gi07Ni5g7bONiINkYTqQwnRhuD2u9lRvmNQxiaEEEIIcTgkJI8wbf1sEp+Tk8OcOXOIjo4elHvW1BgoLTURF+dGUUBBg811DiZ/FA5tJVqdD7fHx6iU0XTqOmmoaSDVnIqCQoevA4PGgN/l73fsQgghhBBDQULyCOL3+3G5XD2O2e12NmzYQHZ29sFfp58NQHQ6HRERESQkJBAXF0dkZCShoaEYDAacTi1ut0JIyN664kjveCY6riO+czY6jLj0tfg0fsYnjMexzcGuLbtQVZU6Vx2jQ0eTEJJAZWXlgH0thBBCCCGOhHS3GEE6Ojr61CN/9dVX+Hy+gw7JJW0l5NTkUOYoQ9WqhJpCmRg/kR9O+iFjY8f2+5rQUJVVq7woig5FcQTHEOkdT0T7LVjtMymyvElsbCMTolMIzQplw/cbUGNVxiaOZX7CfDSKBqeznea6fCLDdKCxgj4VFPk5TgghhBBHn4TkEaSjo6PPsby8PCZNmkRsbOwBX1/SVsLbJW/jxMmExAmEmcJweBxsrt9M9XfVLJ27lMzYvntspqUpTJ6sJz9fz4QJZlpb7bjd7sCTqgZt2ZksmBhLROI/KHOUEjk2khh7DIWfF/Kzm37GaOtoTEoJEdocaChH7TShaELAkAHWxWAcxH09hRBCCCH6ISF5BOldauFyufjmm2+47rrrDvhav+ontzYXr97LnFFzUBQFgDBjGFmxWRTUF7CqaBXpMeloes3uajSweDFUVMC2bTqSk6MxGDqorW2jtlrLeMNOFqe2EG85i/IkcPpdeJO93Hbtbfz54T9zwgsRxBreRa/YcXvjcHjiCQ1RwJ0P3gqIXCpBWQghxDHLr/opt5fT5m7DarSSGp7a53ulGH4kJI8gvUPy2rVr6ejoOKhSi1pXLU2aJsbGjg0G5C6KopAclkxhQyHl9nJsEbY+r8/MhKVLYeVKKCoClysEW8cuLut8i4nezVj/40T93EiMzUbLggW4Ro/joQcf4uc//xllW54mZZaFDv9oQKHd0YHZEodGkwWeAmhbBYZ0Kb0QQghxzCmsL2Rl0UqKGopweV2YdCYyYjJYnLG439++iuFDQvII0jsk5+bmkpqaiu0gNpS3RFrwVnux6C39P2+wUNVWRZt73x0oMjMhPR3Ky6FzUyEJf3+WUGMDrpgk7F4vitOJubgYfU0N9ZdeyqxZs7h1yYW01L1PTcPZhEcFwrnf76e9vZ0waxhok6GzEDzlYDjw+xBCCCGOlsL6QpZ/t5wGZwMp4SlY9BYcHgf51flU2Cv2WaY4ENra2vjZz37GqlWrCAsL4+677+aDDz5g2rRpPP3004Nyz+ONTM2NEKqq0tnZGXzs8/n44osvmD9/fp+Z4d7MZjO2JBsmnQmHp+9ufQCOTgcmnQmr3gKlpbB5c+Cjv+dOeRoN2FL9TNi6kjB3A5qJWZgTEoiKjYWwMDpGj0ZvtxOemwt+P1devpiYqFD+l/MNXq937/0cDnx+H2gs4HeBX9rDCSGEOHb4VT8ri1bS4GwgKzaLMGMYWo02WKbY4GxgVdGqQdtR9vbbb2f16tV8+OGHfPrpp3z55Zds2LBhUO51vJKZ5BHC7Xb36GxRUFBAc3PzAUstFEXBZrNhNBnJiMkgvzqfrNisHsFaVVUqWys5051M6vNvQ1ExuFxgMkFGRqAgObPbT8rl5YGai5QU2HMdo9FIVFQUTU1NuOPiCCkpwVBTA6PCyZw4g+JdX7J+/Xrmzp0bvGd7ezvhFkBjCnS7EEIIIY4R5fZyihqKSAlPOawyxSPR1tbGa6+9xltvvcXpp58OwCuvvEJSUtKA3ud4JzPJI0TvUotNmzZhNBqZOHHifl+XkJBASEgIGkXD4ozFxJhjKKgvwO6y4/V7sbvsFNQXMKlRy8WfVaPJ3wgxMZCejhodQ/tX+TTev5yq/xXunVRuawuEaEvP0g2DwRDY0MRiQensRON00qkmoDdPIvvEMXz//ffU1tYGz3c6Hfg95WDIxK9LprSllM21myltKR20n8yFEEKIg9HmbsPlde23TNHlde23TPFw7dq1C4/Hw5w5c4LHwsPDSU9PH/B7Hc9kJnmE6B2St27dSkZGBjrdvv+KLRYLiYmJwceZsZksnbs0uAChqq0Kk87EjPhpXLepnihHFWRlgaJQXw9FRWE01GeR0FxAxfZV7LgoncUXaci0WqnxtuDZXURKSs+QrtfrcXU20Eo7frMZ0NDiW0BSSjUnzKhn3ZrPWXjOhRh0HoxKHQ53InW6yfxt9eOyKEIIIcQxw2q0BssUw4xhfZ4Plika5Tehw5WE5BGid0jesmULp5122j7P12q1jBkzps+viDJjM0mPSe/ZyqbZj1L+APawFFx1Cg4HFBSA0wnh4Qr+kGRGtxWS+3U5y6tsLL05FePYNDrWrKY01BwsubBF2ChtLsFRWoxm2mw6ExJY37geVVU5OfYykkabqKt/h92luYwfPwmHL4OCmvH8o+lDml3NR31RhBBCCLEvqeGpByxTnJE4g9Tw1AG/95gxY9Dr9axdu5bU1MD17XY727Zt49RTTx3w+x2vJCSPEN1DcmNjI9XV1UyaNGmf58fGxmIwGPp9TqNoetRP7crbjHuji2IsdPqgqSmwXi8tDYxG8PothLdXMTG1jdwGWPWhhmVX3051dS0tG9fiToylTa/SVFNCaM1uYpOTibv2Clb6y1nXsI5ZMbNwqaPBehu72kz8/lev8NsHf8zocVP4su5jqlobmJV2aL2bhRBCiMHUVaZYYa+goL6A5LBkLAYLjk4Hla2VxJhjWJSxaFC+P1mtVq666iruuusuoqKiiIuL4/7770ej0Rxwsb44eJIsDmDFihVkZWUxe/bsoR7KfnUPyVu2bAHYZ0hWFOWgduADKCyEl961UttqIkLvIDQUfD7weqGqChwOMHQ68OpMeExWkpMDrym3ZDLq/x4jZM5JaJpaiCrdSWjDesKndhB3vZGakN9j9f6JReMzmBk9c8/dNCy6YCHX/iieztpnCedZ5kd/zhWjqolSanuMq/eiCCGEEOJo6ypTnJ44ncaORrY1bqOxo5EZiTMG/TedTz75JCeccALnnXceZ5xxBieddBKZmZmYTKZBu+fxRmaSD2DJkiUsWbKE1tZWwsPDh3o4/fJ4PPh8vuDjLVu2EB0dTUJCQr/nR0RE7HMWuTu/P7A5yHZ3KgvGZpBYk09JSBagEBoKHR3QUK9iM1dSnTgDe3gqFl8gPLe1AZMzGfPYC/z1vWuZ5l9Lh1GhcVQyxUoH+rZOpoYbibFuo6w1k8b2OExKCbGGd7nw3NG8/+EaOvPr0adqGRvaQgpfk+89hQbf3nDfu3ez7HgkhBDiaOu3TPEofP+xWq28+eabwccOh4Pf/OY33HjjjYN63+OJhOQRoL9Fe5MmTdrnr1wOdha5q5NbcqqGovjFhLdWkGwvoElNxuuzEKV1EN1cSXNYDJ9ZF1FdpiEyMtAZzmoFVD/flrxDXEwhPpOHZv0Y7K42wk0mOlUTsRGngqeA1KgNtHt+QIQvB71ixxsyifRMlW83fkdqXArlSig2Szvj9QU0+k5FJfC+ui+KkB2PhBBCDJXeZYpHQ35+PkVFRcyZMwe73c6DDz4IwAUXXHBUxzGSyTTbCNA9JPv9/mBI7o/JZMJqPbiVtt07uTXEZvLd3KU0pEwnQddIYts2Iv2NbDXN4LOspVRHZFJSAvn5gZbJqfGF7NrxU7Q1DzHDWsq4UC1jjDX4vbXYXXZUVErtZaBNRuMpZvyoRkzaMtz+OEBh0qRJRIZE0FzVTLvXQaPbTCQ1hGvtwN5FEZkxmTg6HSz/bjn51fnEmGNIj04nxhxDfnU+y79bTmF94RF/jYUQQohjzRNPPMHUqVM544wzcDgcfPnll8TExAz1sEYMmUkeAbqH5NLSUhwOxz5DcnR09EFf12oNzAo7HBAWFgjKq2PS8Y0qpWnXNhxAhWYCmTE2IvWBDfja2yFrbCHVlctoqNtBZ2sKIckKXr2KydPMiTEmytw6/FojNa1F6P0xjApRMWqawaylsT0EAI1GwwknnMC/c/9DeFw4dTiw6n2o3jbsbiW4KOL/pf8/Pij+ILjjkSzuE0IIcTyYPn0669evH+phjGgSkkeA3ov2FEUhM7NviYGiKIcUklNTAxvq5ecH2yMTHV5MxnkrCTUU0dLswuM10ebK4PuSxSxcmMmoUX4s6kp2l9rZbZ/NrNGRqPqdtHa6sRqSiDB4MOkc+BUverUVl2cbHaqJEN0ojHqwmMCx5+0kJyeTFptK/ff1TJ8/Fqe3ie9rduEzpTAjcQaLMhYRog8Zsh2PhBBCCDFySUgeATweT/DzLVu2MHr0aEJDQ/ucFx4ejl6vP+jrajSBHacrKgJ9kWdMKmT2uOUYtQ1U1KSg0VmYMtlBhCWf8/QVhCYvRaMNYdPnRTgdM4lUwkhOUWluDifS0ESoKQJ8LYRrKkEbBZpYFDx4MBHSuQv8DViNHXT6xuHxeAE44YQTePfdvxHV6sCcdDpnJ57DtMwZjI8bj0bRsLl28wF3POq+uE8IIYQQ4mBISB4Bene2GIhSiy6ZmbB0Kaxc6Wds5Er8ngZKarOIiVHIyIDYmDBQs8BTAI5VrCs+lxCNC4/fgqpCabmCbdQJ4PoOvHWgtoPqAUzgr8ei9YFGAd9u8Leg8TcSFeKn0ZeA128kIkzD2QtSWbexjAzDnWR5QrHk72J3Uj3NsVbsnjaMWqPseCSEEEKIASUheQToCskdHR3s2LGDiy++uM85Wq32gC3scktzUVDItmX3OJ6ZCS3Gd0lozyVUO4lMg0J4eHAjvcAn2mRqKgr5Pv8Uzj3JxMQpDkorwygtBYjFNmouuPKhczeoKqgOUN2gMQVmldGDYgRfLVp/AxFmPe1OL6pqICZ5ITn/9zaGvz6NzZZIXWsVDqOPiqRQNp0wmtpIF3WOOuYlzzuqOx4JIYQQYuSSkDwCdIXkwsJC/H5/vzPJERERwQC5r37CCgo5pTkAPYJyXmkeW3fncXFcJ7ERFuins1xppYX2piqmTLeSkJIB7nxsqYGeysGgnDQZvLtAkwwaA/hbQBu7N21rrOB3gzYWg2kybuU0Gpo9aMo6uMO8gdLN6/kuwUlnlJ5YxcyUKi9pn5TSeIqVb6wtfFv5LVmxWUdtxyMhhBBCjFwSkoe57qUWBQUFmEwmxowZ0+e8qKgogP32E+4Kxt2Dcl5pHjmlOZydmk2cyQF+B2j7ljVolXbi47zEjqsG5oCnHDwF2JKTAQtaxQG+CtCEgTENOneAtvt0NKB2gkYP+lTw1RCfNIGmdgcROc+TEBZG7ugQvG31TIyZSKcCjaOMxNY385MqLVWzkvADDc4GqtqqMOlMwcV90idZCCGEEIdKQvIw1z0k79ixg7Fjx6LT9fxr1Wg0hIaGUlhfyPLvltPgbCAlPAWL3oLD4yC/Op8KewVL5y7tEZS/KPsCn+pjgW0B89JOgcZd4M4HTVbPcOutIyXyC0APLS+BJgQ0UaAdBf5GbAlVgbIK/UmgnwCdBXvqkrstIlRVUFtBmwi6BOjcjkZ1MEanw1lWRkuMlShtNBUllTh9jaRMNGMKc+IfB/Guaq7VpPB1SBgXzfop4cZw2XFPCCHEiDV//nymTZvG008/PdRDGdEkJA9zXq83+PnOnTsZN25cn3NCQ0NBgZVFK/ffT7jgfdJTLyW7PYpNNa00x1nRanV7Sy+si8FbEVikp00GjQU8FeD6ClAh5BTQpQZmm30VoImGsCtBFx8opdCnQmcxNPw28NFvDBxXOwMBWTGDMQP8zkCo1lgJ8bahaDTUGLTojXqSbVGEj3FgDDfg7dCiosPi8hKvOpjpbSfC28qk1JOPxpdeCCGEAED1+7GXl+Nua8NotRKemoqiGbxJmvfff/+QulWJwyMheZjrPpNcUVHB6aef3uecsLAwyu3l++0nPKPFzNgP/oHTtx5nWxOnepppSoul4OR08tLyAkHZmAmRS6FtJXQWgacSvDtACYGQU0EXF7igNiww2+wpANdaiF4GXTO6xkyI/hXU3bYnKLsDJRbaxEBA1sQEXmecAfpUvrW/i03nJtSvR6vRMG5GJD5dK027XRgsBgzeTkLQ06CEYKWdkOZ8VPVsFJlBFkIIcRTUFxZStHIlDUVFeF0udCYTMRkZZCxeTGw/exYMhK4SSjG4JCQPc10h2eVy0d7eTmxsbJ9zwsLCqGut22c/4Ziyemb/czP+hgZKx8TSGOVnjGEC01og/dMKVvIuwN6gbEgP1By7i8D+Z9ClBeqLu9vT8YLOwsC5Btve50wTIe4paHwEfHWAbk+HC8OeWeoYsC4CRYMpoZqq2HpGlTmJjQ1DH9aM126g3dWO2WAmvENDe0IYddpOYjWJdNY0UFq8m7QJyQziD/FCCCEE9YWFfLd8Oc6GBsJTUtBbLHgcDqrz87FXVDB36dJBCcpSbnF0SIwY5rpCckNDA9C3F7JeryckJASr0YpJZ8LhcfR4XvGrZHxZhLG5jR3xeiqVVtKix5CSMhGyshjlNnLuVg+5uz4nrzRvz4s0gdCrHwXoQNN34xIgUI7hd4G/n408TBMh5v/APB/w7WkPt2cGOXJpIIw785gWbsd80UIqDB2MbvdgUsCvAaMWwls68Zr0fB+qYm8yU7ttIoUFnaxY3s6jj0Jh4eF/XYUQQoj9Uf1+ilauxNnQQGxWFsawMDRaLcawMGKzsnA2NFC0ahWq3z/UQxWHSULyMNcVkhsbGwGIiYnp8XxYWKATRWp4KhkxGVTYK1BVNfh8eJ2dmPJ6KsIhPCSC9Oj0vds3KwokJzOmxsVCQyYqgdfl5kJeHoF6Yo0pUIO8R2lZ4A8QOL6ntrhfxkyIvgcSXoDoO8GQCSHzggEZRw5YFpB51iOwdCmVoyLRd0KUx0c4esp1XtbFhVDTkYqhfjZxZgtmixFLWCj5+bB8uQRlIYQQg8NeXk5DURHhKf2XMYYlJ9NQWIi9vHyIRiiOlITkYa73TPK+QrJG0bA4YzEx5hgK6guwu+x4/V46mxtxtjWhCbVyQvIJjI4c3fMGFgu4XMwLn8h823wgkJ1zciDvm1QwZAQW6alqICCX7nndVyXwVX4g+Op7beSRlxdI2rB3VjriGghbHAjHDQ8FAzLmwKLB2WdcyVcXnkbJ+HTsWQkYTj6T7zVWKnZOQOsYhTG0mQhzNa2ecSimJLKyoKEBnnw/j8935R75F1oIIYToxt3WhtflQm/pW8YIYLBY8LpcuNv6+W2qGBYkJA9z3WeSdTpdn131ukIyQGZsJkvnLmV64nQaOxrZ1riNGo0Da3gs8yInE2vpW8+MwwEmE1j3zgZnZ8OCBZCTo2Ft0WLQxlC/u4DqKjs2mxdbsh18ZfBtO+SP2rtoDwIBOSenZwu5LuZsULSg+gIfzT03NPFpoDRhIm3hkXgMjczJOAFnk5UQbTOh5s3Uexy0Vk8krnYrkfZSlNQc1tTlUFfXz72EEEKII2C0WtGZTHgcjn6f73Q40JlMGK37+G2qOObJwr1hrvtMcnR0dI9f+ZjN5mDP5K5d9rx+Lz+e9GMAHJ0OrHoLqY1vo8nfCPFqr809VOo2VNKePoMxqT1ng7OzYeNGWPGXTOZOW8r4xJVMmlBEQmwV+E0wfyGEj4LVu8CQF3hBV0BesCDwuDdn3t6ArPoCj817NzRZYFtAti2br7a+xfbNr5EQ4WJCqo/OZicefxrqpm1klP2OSG84zTo7IdFeQsKuYqKln3sJIYQQRyA8NZWYjAyq8/OJzcrq8f1XVVVaKytJnDGD8F7fP8XwISF5mOvqk9zQ0LDPUov97bJnixoDiy+EikooKIDk5ECJhcMBlZV4ImL4t34Rk7/U9Mi1eXnQ3AxaLdS2ZNLQls4Zi8sDi/S6eiKfqQkE5Jwc+OIL8Pn2H5C7l1jsebyxZiM55c3BgAxw8sTLaNdE8MGXf6O8dhrtn+zkNn0pMWaoirGzy+DC5OlkSqWRcXxDVO1smCy77gkhhBg4ikZDxuLF2CsqqC8oICw5GYPFQqfDQWtlJeaYGDIWLRrUfslicMnf3AGsWLGCrKwsZs+ePdRD6Vf3covenS0sFktwl7386nxizDGkR6cTY44hvzqf5d8tp7C+EDIzYelSmD4dGhth27bAxxkzGPXoUiZfkhmoQd7T3KJrQjgyElJSAkHZ59OQ940NTJMDNcZdJRbZ2V0nBD4eTECGwEfLAiJ8+SxKjdy7ockeZ2eeS/aUH2CITmJhRye+2jIaLHNo0+rxaxScBhMV2mySDA0krlkFsrpYCCHEAIvNzGTu0qUkTp9OR2Mjjdu20dHYSOKMGYPW/k0cPTKTfABLlixhyZIltLa29qn3PRZ0L7eYOHFij+eMJiMrNx9gl72iVaTHpKPJzIT0dCgvh7a2QA1yaipoNGTv+TfefUI4MjIwk9w1MdwVnKFXDs7L2xuQfb7A495BWVV7BuQu5mxsSX7w1oJr894Z6j0B/IezLyLq+1w04Y1870hDay/EZPaCR4u7U8UYVkZs+iic36+lcv1/MY3P2P9W1ao/0NO5+2y4bEoihBBiP2IzM4lJTz+qO+7ldi1+F4NKQvIwt6+ZZEVRqOmo2e8ue8lhyRQ2FFJuLw+0fdNowGbr9z7Z2XsDclc3m+6VE10fewTl3jXI+0rSlvn9vzl3IXR8E9jdr+39QDs5Q0Zge2xjJlqtlrlZRlpT2ykydNBAPfGNU4jQjiUucTeeiK186/IQW9XMu188RX1lfLDMJDM2s++9unYS9Lv63EsIIYTYF0WjIWIf3z/F8CUheZjz+Xz4fD6ampp61CTr9XraO9v3ucsegMVgoaqtijb3gdvTbHw6l+SdCpVjs/H7AzPJwZyblweqSvb8+UBgYrjfRXr9JuluHLmBhYPm7EBobV4OvgbQpoDaQF1bGYpSTqy3IrjhSHhyMo3hDpJCKmne2UBIZxjTT4zGGm8it7yV9vp6rDorSYkT8Ifo+KrsKzbWbOTGGTdy2pjTArPKve+lswR6PLvzodu9hBBCCHH8kJA8zPl8Ppqbm/H7/T1CstFoxGA0BHfZCzOG9Xmto9OBSWfCatx/e5q8PNi1UWEBOdhOhbxTs4M1ytl0C8N0y725av+L9Loed9vQJEhRArXJqj8wg+xrAH0WeMvAW45PCaG4tYOJFBKrWxXYHjs1lcpUK3Eb3PiT09j69RpmzkymtLoaraIh022lbGw0JVYfhbUF1DvraaprYnvjdi7OupgLMxaRqflg7726Zty1YaDJCmyT3bbnXlJ6IYQQQhw3JCQPY36/H1VVg7vtdS+3MJlMJIcnkxGTQX51fo+aZAi0p6lsrWRG4gxSw/fdniY4IXxNNjaAnByyFwALstn1Sh5p5GC7pp8wvGdWGeg5Qww9z3Xm7alJnr/3+daV0FkIhomBgOwpBb2NRLMNt1LKjtYifOSSEHYpebvL2HzSGM5p9JHZ7KXFWs6WrXnEjTYwsU2PMyaWwpPG8231muAPBRGmCFrdrXxV/hVOZwF3TWgnwjy2b+9mRQFtcmAsnvLAgkQhhBBCHBckJA9j+9ttz2g0BnfZq7BXUFBfQHJYMhaDBUeng8rWSmLMMSzKWLTvhWwE8uveCeG95RLZ2i9Iw0fztAXY+utY0V3XDDH0XJzXvatFF3M2dO4Cx6egugPH9LbAHwjUTqs+Ghwb+GDtE1R7YliQfQmjZ8XQ9OdHOV8txtm8G6NXj3dyGPb5CZRQQnNNM37VT1NHEz7Vh8fnITokmoZ2O9WtbYSHTqbfLUc0FvBWBRbzCSGEEOK4ISF5GOu+aA96ziQbjUZg7y57XX2Sq9qqMOlMzEicwaKMRX0XsPXSfUIY6LGCzzZWi+3Wg9iooysYdw/K/bV9C56/AJqfAb87sIBuDUAJzIvE2bSbONWDHR0dfh2jN5WTjQonaNBf246/PIHvV7ezO8JB2uxo4q31nGlvYnuTluI2DUadEZ2qAxWaXc2UtWmpdjhI6qgm3JzSd+x+R2AMGtkxSQghhDieSEgexrpvJBIeHo5erw8+1xWSIRCU02PSKbeX0+Zuw2q07r8V2v4cTEu3/nQPys4vAjvq9ReQATwloIkGXw2oBvDsgNVV0KGBDB1tzQ7MRgO2jWVE5DezMX4jE8p24m2vhaQMUs9IJn/Dh/hbmnCosYRqOlgYr6XSFYOKgtPrxGq0khiaSLmjjs0tHmYklIOa3GfHQXyVYJwRaAcnhBBCiOOGhORhrHu5Re/d9rqHZACNogmUKhyJg23pti/mbGh5DfCDYWzfgOzMg4588LdA+FXg+gY61sHUcugEvlVQPCqe8RosOzWc1biLLzJ/wGf+dURWFRKbMBFTlA2/sY4kQyJ1tR34lBY0epV0q58EUyc723zotXpizDFoNBqMuhA+qm7nknQLEZ6CQA2yZk93C18laGPAukgW7QkhhBDHGQnJw9i+eiTrdDo0/TQxzy3NRUHps3sdQF5pHioq823z+7/ZobR0671Qr4szD7y7A5uDdD3uOseZBy2vBD6PuCZw3DUTOm8DbznM8uPy+XF97SNsfSK6qDA6MtyMi9hOe3kIjX4HqlbFptFgDbUyLul0dq8rQhPeiENR0SpeTIqHUEMYMeYYLAYLqqrS4e2gxh+JPeQqEo2FgT7J3qpAiYVxRiAgS/s3IYQQ4rgjIXkY6z6TPGrUqOBxk8nU7/kKCjmlgUDbPSjnleaRU5rDAtuCfl8H9F7Bt1d/Ld36W6jXFYJ1CRB6DrT/Bxoega4JcEcOmKaBafqex7mBxXraeLCcB74G2ua2oNvQSpg+HDoUOtKiiPTXkTFqNi2dadDgwhYJ9Q1mOjoSSAjbTvv2OlLmRaMoDtBYibfEY9AZcHs7MHrsjNNribEmYzRPhYgLZMc9IYQQQgASkoe17iF56tSpweO9Sy3w+6G8nOy2KIy6DD7Z9TkQCMrdA3J/M8xBfVbwddM7OPdeqAd9Z4m7jtXcCqape493X9Dnbwt0uNDbQJ9C7IZS0JeCVsEZ0gEVPhjlJzFRR0jbBFztNeR9oaKqChMm6ElNHcV7739HkiaEQlckWsNoOjqaCPfaOc2gMiZUh1FRiDb7sFb9m07N2RjCx+37fQohhBD9Uf0yyTICSUgexnw+X7BPcn+dLQAoLISVK6GoCFwu5plMxMUb+fekd/mi7At8qu/AAflwdA/Knj37WHcPyFG/APcm6NgA3rr+O150lgbKHvwO+K4JviqFk204M0x0fF9EyNZWzNpYnPPmgbKW8rI6YqwFtHUkM36shZamMBacEMuugjYK4qZi1oWwMCGe2Z4qTKoLB16M+hDSotLobN5Ki7OGiPSrJSgLIYQ4eO5CaFsZKNfzuwLftwwZYF18zJTrqaqKz+dDp5PYdyjkx5xhzOfz4XA4cLlcfXokA4GAvHw55OdDTAykp0NMDGNKW5j33hoiS2vQKtqBD8hdzNmgaAM/UfdeqOfMA10qhMwA1QnV1/VtCadPDfyP5qsN8FXJ3oDcXEPIlHTMk8LgOy/m7eFsqz+fmrpxtDQnYDI0Urd7G1ZLJ5GJZ/D7Z110bLYxxjiasZ01hKothIe0MjWygznRnURpt2HQ1eLtKKW98r+oqn9wvh5CCCFGFnchNC8Hdz5oYgK7s2piAo+blweeH6xbu90sXbqUuLg4TCYTJ598MmvXrgUgNzcXRVH4z3/+w8yZMzEajXz11Vfs3LmTCy64gPj4eEJDQ5k9ezb/+9//elzXZrPxu9/9jmuvvRar1Upqaip//vOfe5zz9ddfM23aNEwmE7NmzWLVqlUoisLGjRuD52zZsoVzzjmH0NBQ4uPjueKKK4L7OgwXEpKHMZ/Pt+8eyX5/YAa5oQGysiAsLNC2LSyM0iQzIXYHmau34fd5ySvNO7gbOnID4bY/zrzA872Pqb5AUFZ9e1/bfcY48eVAgO7cCd6KnkFa0QR+EteEwVwF5kaC6iMk3Io5tAlOzoQzr2LjRoVPvp5HhO18pmdfiiv0Qf6Rcy9rSx8matyjjM08j5V/XslFkacyLySGcWEWbJYwYkNtGA1JoJhR/DXolEo6m9fjdew+uK+HEEKI45fqD8wg+xpAnwXasMD3O21Y4LGvAdpWBc4bBHfffTfvvfcer732Ghs2bGDcuHEsXLiQpqam4Dn33HMPjz76KIWFhUyZMoX29nbOPfdcPvvsM/Lz8zn77LM5//zzKS8v73HtP/zhD8yaNYv8/Hx+/vOfc9NNN1FcXAxAa2sr559/PpMnT2bDhg089NBDLFu2rMfrW1paOO2005g+fTrr1q3j448/pra2lksuuWRQvhaDRebdhzGv17vP3fYoLw+UWKSk9Oj9W9pSSqm9jLHjJjLRFYLFkMnH/Szm69eh7JzXu3Si67FrI/iaex7X7dnEo3MnND0TKMUIvplM+MFjwV9lmUO7fpUV6DyRZ8ikZGMui6bnMe2swJjmJYDbDx9/Bv8v/HPuPeUUTv6knOVP1XLHRfUQolDnTsJsVoiOAYvZCGosGrUOr7sEv0d21xNCCHEAnvJAiYW25/dZIPBYmwydhYHzDLYBvbXD4eC5557j1Vdf5ZxzzgHgxRdf5NNPP+Wll15i9uzZADz44IOceeaZwddFRUX1WMP00EMPsXLlSj788ENuvvnm4PFzzz2Xn//85wAsW7aMp556ipycHNLT03nrrbdQFIUXX3wRk8lEVlYWVVVV3HDDDcHXP/vss0yfPp3f/e53wWMvv/wyKSkpbNu2jQkTJgzo12OwSEgexnw+X5+QrNVqAzVHbW3gcoHFEjy/tKWU0pZSbBE2UkKTYds25oVPxG1N7LfrRR8Hu3Nef8e6L9YzTe//tU3PgCu/Z2s4CARlQ3q/iyJUFaZNV5jWnAN5BBcRZmdD+MY8IvK/oG3WTNLS/sDuyvWYTW341HD0eoW2dnC5YdQosJgV/H4zGtrQYD/8vxQhhBDHB39boAZZZ+n/eY0l0FLUP/ATLzt37sTj8XDSSScFj+n1eubMmUNhYWEwJM+aNavH69rb23nggQf417/+RXV1NV6vl46Ojj4zyVOmTAl+rigKCQkJ1NXVAVBcXMyUKVN6dNKaM2dOj9dv2rSJnJwcQkND+x27hGQx6LrKLYxGI5Y9YThYj2y1gskEDkeg1GIPW4QtsKmI3R543mol2zYZABW19y36Opid81R1H9tNd2sX11+QjvrF3uPdz4dA6UU/P4nPnw/Mzw4E5G79mvPfeZrIko2kXnUNj35zKsnJjdhrP8Dl8aBoW/BofVjMUTic0NgAIckq3k4PJrMBnSnkwF8HIYQQxzeNde/icm1Y3+f9jsDzGuvRH9seFkvPAH/nnXfy6aef8sQTTzBu3DhCQkK4+OKL6ezs7HFe9x18IRCU/f6DLxtpb2/n/PPP57HHHuvzXGJi4iG8g6ElIXkY8/l8tLW1ERYWhrLnVz3BkJyaChkZgUV7WVmgKHt33FNVqKyEGTMC53EQpRbdmbP3BmRF2zcMW+bv/7UQqF8+UJA+FN03NvniCyIbd5Jjg5AQyC3NI3yqB2vxdMoc/yYhQkeIqkEhFEwOvP5WOlp0mMyhhMbGoWjD917Xmbcn9O/nPQkhhDj+dC0ud+eDJqtnyYWqBnZtNc4InDfAxo4di8FgYPXq1aSlpQHg8XhYu3Ytt9566z5ft3r1aq6++moWL14MBMJsaWnpId07PT2dv/71r7jd7mDm6Fow2GXGjBm899572Gy2Yd1RQxbuDWM+n6/Hf6TQLSRrNLB4caCrRUFBYObY6w18LCgIHF+0KHDeoeq+IK+rjnhf5/VezNfFMr9vQO5izj68UJqdHVic6PNhix7LmEXXkFeWQ71mE03hq2nKrOT7dgWXDzy+RHw0olPq0fh9GEweIuI7MUTO3Ps/tK5Z7d61ZkIIIUTX4nJtDHgKwGcH1Rv46CkIHLcuGpR+yRaLhZtuuom77rqLjz/+mIKCAm644QacTifXXXfdPl83fvx43n//fTZu3MimTZu47LLLDmmGGAi+5sYbb6SwsJBPPvmEJ554AiA4YbdkyRKampq49NJLWbt2LTt37uSTTz7hmmuuCe7xMBxISB6mVFXF7/fjdrt71AX16JGcmcmai+axyxYBjY2wbVvg44wZsHQpeSF15JbmHtqNu5dJxNwbqC9uW9k3KA9FwMzLA58vGJSzyyA7bQGtlnzKvWuosH7GdxWz2bnDittSQ4XBy25iqVJtWJPCMZjcoBsFKlD4Dqx9BxoywHTK0XsPQgghhg9jJkQuBeN08DdC57bAR+OMPccHr0/yo48+ykUXXcQVV1zBjBkz2LFjB5988gmRkZH7fM2TTz5JZGQkJ554Iueffz4LFy5kxowZh3TfsLAw/vnPf7Jx40amTZvG//3f/3HfffcBe3f8TUpKYvXq1fh8Ps466ywmT57MrbfeSkREBJrDmZwbIsN3Dvw41/WTmMvl6n8meY+OcTbe0OxioSGTeeETA7XKqanklX954K2oe9tXHTEEgnLX4/7OG2x5eYFSi66ts/c8viR7Ab8P11De7iRCZ6HVksrzuyq5MLKeOWNMYGgg3GAiJPps0I+Crd/Bv1+H7fXgjYRQB2TsCszKZx4bTeGFEEIcQ/azuHwwmUwmli9fzvLly/s8N3/+fNR+yhZtNhuff/55j2NLlizp8bi/8ovu/Y8BTjzxRDZt2hR8/Oabb6LX60lN3Vta0jVrPZxJSB6mukJy75lkg8HQ47yuWuOPS3NwWxPJtk0++K2oe9vXgryuoOzKh4aH+l/MN5h6B2QIfixf+QqXR8ILyhgaO3bRqduMVj+B3+S6mVmcxezoE7n9J2ehiUmFomJ48Y+ojQ5cCSZaQydgcvuwbtiApqICli6VoCyEEKKvfSwuH6lef/11xowZw6hRo9i0aRPLli3jkksuISRkZC18l5A8THUPyd1nj3uHZNgblHNKc45sK+r91QlH/WJvQO5vMd9gUtWeAXmPvDTYZYMFCdM4e+qtXPv3W/i+85/oYqpRNB6Kttu595wFjM1IC1zj3Sdx1TVRM8pPh7eRmromWrzhxIRGM7WyjbBVq1DHj8deWYm7rQ2j1Up4airKMPrVkRBCCHGkampquO+++6ipqSExMZEf/vCHPPzww0M9rAEnIXmY8nq9QKDcomsmWa/XB4vme8u2ZQcDcp+tqHNzA7XD2f0E27w93R3mz9//gPrbXe9IgrJjz5j6u0bvjhPz5+/dDXDP+cHZ8kXXMC0ONlY/zVmnRGPcmUW7y0GbtZVtO1fx74J4TjrpYaLrPsX1/RdsDVHZ2RKCzWLGZm6n3uOl1FGDW6Njat6nlDQ1UVVbi9flQmcyEZORQcbixcTKDLMQQojjxN13383dd9891MMYdBKSh6nuM8kRERFA/7PIXfJK84IB2af6yCvN2xuUFaVHj+G9L+pWxrA/+9pdDw4/KB/K7n79nK+iBmbL46C06hXy6yAqZDqXpKRS56ijKCYR99Zm3vn2Hc470cq5/laqTZ3UheqJdcXiRKHFpyfO0IJeE05ZRRt1Wytp9mkxT56M3mLB43BQnZ+PvaKCuUuXSlAWQgghRhAJycNU94V7XTPJWq0WAL/qp9xeTpu7DavRSklzCXllecESi65ZVthTitG9xzD0WPjWXxlDD/vbXe9IgvLB7u63j/Pn2+YHz2/WTqNVq9DsamZhwnQmKB1sU0PwN/n5+2t/Z9OmjWROyKA1AcLbjHTumY23+yPACwoq0VUddHjcMDYZ457NWYxhYcRmZVFfUEDRqlXEpKdL6YUQQggxQkhIHqb6626h0+korC9kZdFKihqKcHld2F12vH4vV029Kjhz3L1GOfi412Yc+HwHDshwcLvrHa6D2d3vIM5vdUCzK7BQ8QRbNs7ab5lQ+y03zTqJDfkb+OyTfE5KSaHNbWJ0kZeG0WqwdZ3dHwFtHcQ2OtkRHYLLYiS+2y0VRSEsOZmGwkLs5eVE2GyH/36FEEIIccyQaa9hqr/uFrvadrH8u+XkV+cTY44hPTqdMGMYPr+Pbyq/obC+MPj6bFs2C2wLem5F3W0zDrTaAwdkGJxNQXpfo6vO+WAWBPZzfrD0Ys8PB+b4eYTEz2OKSeXps2YzOdHPZ+sa+TwtDbtVT2xpPUaHC8Xnx+hwEVtWT7tGYWtqJEa9qc8tDRYLXpcLd1vbkb1XIYQQQhwzJCQPU71nkv2qn08rPqXB2UBWbBZhxjC0Gi0T4yaSbcumwdnAqqJV+NW9O+tk27IDZQldem3GQV7eUX5X/Wh6OrCrX/cFgV3629GvvwWEpfOhrGe4NsfPA0VDvH4WKdb5PPXMZ7TFJfPu6QlUj08gxN5BdEUjIfYOdo9L4PNRZrwRcYQb925ZraoqrpYWWsrK8Hu9GCyWQfsyCCGEEOLoknKLYap3C7iajhp22neSEp3Sp8OFoigkhyVT2FBIub0cW4St7wX3sRkHcHAzyoPBmQeujYHPzacGPnbVKHd93n0B3z4WEEYYYNV/95Sa7HkrztpvKSv1U1auZ/q0ycwcn0PV/6qIPT+BPyc5meS2EenR0az3skVvJ/VjH5nNoah+P56GKlzVu2mrrsPu9OBqbsYcG8vmt98m88ILZQGfEEIIMQKM2JDsdDrJzMzkhz/8YXBP8f0dH256l1s4vU46/Z1Y9P3PZloMFqraqmhz91MSsJ/NOIYsKHcF3ohrAo+7wq9lAbS8EjgWcc3e8ov9LCCcNj7wHlb9N5twfS4Tkosp2tTBtwXzmDhvHrPSv+GndWvI+yaPhSHXUxA2gW1N23B52zHpTEyPmcVpP72Gmt//CfWvfybK6cDQ7sbrh2aDjqr4JKxTplCzcSOtlZXS6UIIIY43fj+Ul0NbW3BnW2Qh97A3YkPyww8/zLx58w76+HDj8/nw+Xx4PB6MRiNmnZkQfQgOj4MwY1if8x2dDkw6E1ajte/F9rEZR/DxkSy+O1z9LQh05ATKKACM03o+d4AFhNOmqtg9sHNzMf6K/1DScg4T583b8xZP4LSFpxGm28CW3P9x+dJ/w1SC3UFSw1Mp+frfdDZuw+ztwOXy0e5XMOi1xKh+zO5m7NpOIqTThRBCHH8KC2HlSigqApcLTCbIyIDFi2WX1mFuRH4X3759O0VFRZxzzjkHdXw48vv9uN1uILB/e0JIAhOiJlBhr+izX7uqqlS2VpIZk0lqeGrfi82fv++Z4uzsA28kMhh6LwjsviDPMBaib93/+d3tWUCYnQ0NrVMpaTmHmJhqsmfvqW925pGQ4EATdRor/1PCv//1b+IMcUyOnxwoTfH7KX9tOTqPi9ZZGXgiLJAUSWdyNM22eIw+L/pNa0BVe3S6EEIIMcIVFsLy5ZCfDzExkJ4e+JifHzheWHjgawyAzs7Oo3Kf480xF5K/+OILzj//fJKSklAUhVWrVvU5Z8WKFdhsNkwmE3PnzmXNmjU9nr/zzjt55JFH+rxuX8eHI1VVgyHZaDSiUTRckH4BMeYYCuoLgq3f7C47BfUFxJhjWJSxCI1yzP2VH5z+FuQdorw8qHbMo6Dup+yqWUBpQU5gK+09ZRrZ0Zcz2ZPMU089RUlJSfB1u7d+S/SazYToQ1D8/kB7OIsJDDrQaHCEmrE2NuGpr5ROF0IIcbzw+wMzyA0NkJUFYWGBhe9hYYHHDQ2walXgvAE2f/58br75Zm699VZiYmJYuHDhgN9DHIMh2eFwMHXqVFasWNHv83/729+4/fbbuf/++9mwYQNTp05l4cKF1NXVAfDBBx8wYcIEJkyY0ON1+zq+L263m9bW1h5/jjUulwsg2AJuYtxEls5dyvTE6TR2NLKtcRuNHY3MSJzB0rlLyYwdpr/26V5vHHNv4KMj55CCcvey63vvhTGTstlVoqW0dG+rOEtoKHfMmEHSjh28/fbbNDQ0AKB+/jmRDQ5Ukwn0OtBqwOsLXtunN6D1+vA72+l0ONCZTBit/ZS1CCGEGDnKywMlFikpwd76QYoCycmBmeRB+s3ia6+9hsFgYPXq1Tz//PODco/j3TFXk3zOOefstxziySef5IYbbuCaawILup5//nn+9a9/8fLLL3PPPffw7bff8s477/D3v/+d9vZ2PB4PYWFhdHR09Hv8vvvu6/c+jzzyCL/5zW8G5T0OhO4zyV0hWafTkRmbSXpMeo8d91LDU4f3DPIR7ujX77rE2XmUFvjYVaIFfNjMeZCdzZjOTm7csIF3li9n27nnEvn994QXl1IaY8Fu1YPFCOFmaGwHvRYUBa2nE59OixJiwV5ZSeKMGYSn9lPWIoQQYuRoawvUIO+r/afFAlVVgfMGwfjx43n88ccH5doi4JgLyfvT2dnJ+vXr+eUvfxk8ptFoOOOMM/jmm2+AQLjtKql49dVX2bJlSzAI7+t4f375y19y++23Bx+3traSkpIy4O/pSHTNJHftuNe1LbVG0fTf5m04GoAd/fqsS9wTvG1ZCyhzZNPszsO2J3AbzzyT6bffTsktt1B78820Z2dj/eHlNPoqMX1fgGtCKEpqDDjc0OIEsx5Lm4PmqCjs9a2Y4+LIWLRIFu0JIcRIZ7UGFuk5HIESi94cjsDzg/SbxZkzZw7KdcVewyokNzQ04PP5iI+P73E8Pj6eoqKiAb2X0WgMhs9jVfeQrNFo+vRHHhH2t2PfgXbf26PHusNeM9OB4JwNToIz0xNuuIFJL73E5vx8srKyMJ54IqnGpVT9dhmmbbtwJ8aijk9At6saa10LLr2O1lFjSZw1i4xFi6T9mxBCHA9SUwNdLPLzAzXI3b8HqypUVsKMGYHzBoFFNrAadMMqJB+qq6+++pCODye9yy26ZpHFAexnZjq/Jh+lNZ9plZB90kkUFhSw5ptvsLz3HpUnJdN8w+WEf5yDftsONG4P/tgQmqeOJ/KcS1hw6nmEp6bKDLIQQhwvNJpAm7eKCigoCNQgWyyBGeTKykCXi0WLpF/yMDasQnJMTAxarZba2toex2tra0lISBiiUQ2d7jPJEpIP0n5mpls109m16hUiSjeStvhqvDEx/OW++0h86zl2VqWScclPOeWCu9i99VucTbWYo+JJmjgPjXZY/TMSQggxUDIzYenSvX2Sq6oCJRYzZgQCsvxmcVgbVt/dDQYDM2fO5LPPPmPRokVAoF/wZ599xs033zy0gzvKes8k63TD6q/ymJRdBmmlkGODMTaFn0/+OX/+75951rmV+7dGM61URbHpSJ5y8lAPVQghxLEiMzPQH1l23Btxjrlk1d7ezo4dO4KPS0pK2LhxI1FRUaSmpnL77bdz1VVXMWvWLObMmcPTTz+Nw+EIdrsYaCtWrGDFihXBbaCPJd37JCsahdKW0pHR0WKoqCq2xdcwJg1ySnP4QvmCmRfN5IOnP+DicXGkNTURNdRjFEIIcezRaMBmO2q3y83NPWr3Op4dcyF53bp1LFiwIPi4q8PEVVddxauvvsqPfvQj6uvrue+++6ipqWHatGl8/PHHfRbzDZQlS5awZMkSWltbCQ8PH5R7HC6Xy4VWq6Wio4J/FPyDWn8tLq8Lk85ERkwGizMWD9/eyENhzwq/bOCLsi/wqT6mTJ5Cc1ozd330EWOvvZaZnZ0YDIYhHaYQQgghBt8xF5Lnz5/fZ1vl3m6++ebjrryit65yC0OSgbdL3saluEhPTMeit+DwOMivzqfCXjG8NxEZInmlefhUH1pFC1pYfOti/nnxHXz3+OPEPfEEY8eO7fWCvMCCwKHYvlsIIYQQg0J+Hz+Mdbg6IBPsnXbSo9IJM4ah1WgJM4aRFZtFg7OBVUWr8KsDvyXmSJVXmkdOaQ4LbAu4N/teFtgW0BTexNSzp+H46CMa3nuPlpaWbi/Ys1PJSGy/J4QQQhzHJCQPU6qq0uhthBiIM8Wh6bVAQFEUksOSKWwopNw+OFtijjTdA3K2LdAiLtuWzeljTkedk0pxjIZtL7xA08qV+P3+YECujoyk9CA2NRFCCCHE8HHMlVuIg+fwOFAMCiG6kD4hGcBisFDVVkWbe3C2xBxpVNQeAblLti2btbF/x5GRyke5O0j/+9+J3rqV8NBQqiMj+aQkn4h4FRvzh2bgQgghhBhwEpIP4FjubuHr8KEL0dHh7eh3t72EdUWEdLZhPbOfLTGljraP+bb5+3xu9rk/pKm9icbyJjbk5ZGYmIinro5PojSU2mDBvOlHbZxCCCGEGHxSbnEAS5YsoaCggLVr1w71UHpQVRVtuxaz00ydq65PSFZVlXpnA6fu9JK6saTni6WO9pBl27JZeMlPmTvLRojPydaPP2bH92vwe6tZsOiaPrPPQgghhBjeZCZ5GOt0dxLvjifcEE5xczG2SBsWgwVHp4PK1kpiZmWSOv0ENLl5oGggO3tvQF6wIPBYHLTsMkiOy+KNsaW85qzm/5mNLPREkFkG2IZ6dEIIIYQYSBKShylVVXG5XEToIrh09KUUqAV8VfkVHp+HlPAUZiTOYFHGItJiMyEij9KVr6B+8Bqjw1IlIB+OPT9cmMfNwr15HZE7W8lpb+QkawreTz8N/EPKzqZ0TxmLTcpYhBBCiGFNyi2GMbfbjclkYrR1NPecfA/XTr+WjJgMzh1/LstOXhbsj5yXBrtay1D8ftBqJSAfDlUNLtLTT0wlKt1GTUcrf8vPYYvfj+r3U5qXR6mUsQghxHFHVf142itxtxThaa9EldarI4LMJA9jLpcLo9EIgF6n58eTfkxiaCI5pTl8WfYl2bZs8krz2LXqFRaEpWGLHgs+X2BWVILyISlVFHJL8im1wemLfkX4xmYeWfZTvqvZTOyuOOq+NGHw+7EtWIBNvrZCCHHc6LTvoL3yEzrtO/H73Gi0RgzhYwlNXoghfNxQD08cAQnJw1TXjntGoxGNRhNcuNe1gCynNIcvyr4gOX8nC0rBtvianjXJIEH5EGysDgTkrkV6aprKD3fu5LmH76ZQ/Ral0c75d9wnAVkIIY4jnfYdtBS/itfdhN6ciE5rRvU5cTVtxeuoJiL9agnKw5iUWxzAihUryMrKYvbs2UM9lD66yi10up4/62TbstEqWpLzdzJmY9negAyBjwsWBIJyXt4QjHp4ipg3rUcXC0VROOuGG5g5YTbl1S20W/T4kpOPyVaBQgghBp6q+mmv/ASvuwlD2Hg0eiuKRotGb8UQNh6vu4n2yv8OSunFRx99RERERPB7zsaNG1EUhXvuuSd4zvXXX8/ll18+4Pc+nkhIPoBjuQWcy+XCZDKh1Wp7PJdXmodP9aFFw65paeSl9XpxV1CWXeIO2nzb/D5t3pry85k9cyZhSiRFqwuo+OorKioqhmiEQgghjiavYzed9p3ozYl92rAqioIuJJFO+w68jt0Dfu9TTjmFtrY28vPzAcjLyyMmJobc3NzgOXl5ecyXReRHRELyMNZVbtE9JHffWvnKW19mzKJryCnNIa+016xxdrZsJHIEuhbpTVu0iEVPPsnahga+ev11tv33vzQ1NQ318IQQQgwyv7cdv8+NojX3+7xGF4Lf58bvbR/we4eHhzNt2rRgKM7NzeW2224jPz+f9vZ2qqqq2LFjB9lSAnhEJCQPU/3NJHcPyF2zntm2bBbYFvQflMVh6QrItgULGH/mmZx11lnMuvhi3tuwgdLcXDa8/z5ut3uohymEEGIQaXShaLRGVJ+z3+f93g40WiMaXeig3D87O5vc3FxUVeXLL7/kwgsvJDMzk6+++oq8vDySkpIYP378oNz7eCEheZjqvnCvKySrqD0CcpeuoKwi5RUDQlV7dLFITEzkjjvuoNFq5d8FBfh8PkpKSlClnEUIIUYsnSUJQ/hYPM7qPv+/V1UVb0c1hvBx6CxJg3L/+fPn89VXX7Fp0yb0ej0ZGRnMnz+f3Nxc8vLyZBZ5AEh3i2Gqa6ay+8K9+bb5+zxftk0eOL03ClEUhYkTJ3LXXXexbNkyZjc0EONwUFVVRXJy8tAMUgghxKBSFA2hyQvxOqrpbN2OLiQxUGLh7cDbUY3OGEVo8lkoyuDMR3bVJT/11FPBQDx//nweffRRmpubueOOOwblvscTmUkehrpKLYA+NcliaFitVi655BLOOussHn/8cerr66mtrZX6ZCGEGMEM4eOISL8aU9RE/J3NdLaV4O9sxhQ1adDbv0VGRjJlyhTefPPN4AK9U089lQ0bNrBt2zaZSR4AMpM8THXNJEtIPnYkJyfzq1/9iosuuojf/va3PP3005SVlRESEkJISMhQD08IIcQgMISPIzJsDF7HbvzedjS6UHSWpEGbQe4uOzubjRs3BkNyVFQUWVlZ1NbWkp6ePuj3H+kOKSR/+OGHh3yDM888UwLCIOhebiEh+dig0+mYPHky9957L7feeisffPABixYtYufOnWRkZPTpZy2EEGJkUBQN+tCjX1739NNP8/TTT/c4tnHjxqM+jpHqkL5rL1q06JAurigK27dvZ8yYMYf0umPJihUrWLFixTG1SUTvcgsJX8eOqKgozj33XHJycnjyySeZM2cOSUlJ7Nq1i/Hjx/fppRnk90N5ObS1gdUKqamgkWooIYQQYqgccrqqqakhLi7uoM61Wq2HPKBjzZIlS1iyZAmtra2Eh4cP9XCCZCb52JWamsodd9zBd999xwMPPMDzzz9PW1sb5eXlpKX13tkFKCyElSuhqAhcLjCZICMDFi+GzMyj/waEEEIIcWgL96666qpDKp24/PLLCQsLO+RBif3rav8GgZCskRnHY4rRaGT8+PE88MADbNiwgb/97W8ANDQ0UFtb2/PkwkJYvhzy8yEmBtLTAx/z8wPHCwuH4B0IIYQQ4pDS1SuvvHJIs8PPPfccMTExhzwocWBd5RbF7mK+rPiy33PySvPILc09iqMSXeLj4znllFP48Y9/zLPPPktpaSkAlZWVtLS0BE7y+wMzyA0NkJUFYWGg1QY+ZmUFjq9aFThPCCGEEEfVYU9BdnR04HTu3WWmrKyMp59+mv/+978DMjCxf10h2aA38GXFl3120+vafU9hHzWwYlApioLNZuPmm28mPj6e+++/H6/XC0BJSUng3055eaDEIiUFetcqKwokJwdmksvLh+AdCCGEEMe3ww7JF1xwAa+//joALS0tzJ07lz/84Q9ccMEFPPfccwM2QNFX94V7cxPmcmrqqT22ne5ve2px9JnNZsaMGcNvfvMbCgsLeeGFFwDw+/3s2LEDT1NToAbZYun/AhZL4Pm2tqM4aiGEEELAEYTkDRs2cMoppwDwj3/8g/j4eMrKynj99ddZvnz5gA1Q9M/tdqMoCnq9nlNTT2WBbQE5pTk8lPeQBORjSEJCAvPmzWPJkiW88sorfP755wB4PB7KmppQjUZwOPp/scMRWMQ3AhbACiGEEMPNYYdkp9MZrE/+73//y4UXXohGo2HevHmUlZUN2ABFX10L90wmU7ClWLYtG62ixaf60CpaCcjHiK6yi6uvvpozzjiDBx54gJKSEgDs4eE0xsWhlpeDqvZ8oapCZWWgu0Vq6hCMXAghhDi+HXZIHjduHKtWraKiooJPPvmEs846C4C6ujrpaHEUuFwujEYjEAhieaV5wYDsU319apTF0DGZTKSmpnLfffeRkJDAnXfeSXt7O2g01J5wAm0GA6xdCzt3BhbrtbRAQUGgy8WiRdIvWQghhBgCh/3d97777uPOO+/EZrMxd+5cTjjhBCAwqzx9+vQBG+BQW7FiBVlZWcyePXuoh9KD2+0OhuQvy78Mlljcm31vsPRCgvKxIyYmhlGjRvH73/+ehoYG7r//fvx7ulZ0qCrusjL44gv45z8hLy+waG/pUumTLIQQQgyRww7JF198MeXl5axbt46PP/44ePz000/nqaeeGpDBHQuWLFlCQUEBa9euHeqhBHUvt1jfuJ688rweNcjZtmwJysegtLQ0xo8fz29/+1vy8vL49x/+QOzbb2Oor6d5yhTaTz8dsrMhMRE6OoZ6uEIIIY5x99xzD0ajkcsuu2yohzIiHXJIvu+++1i/fj0QWJQ0ffr0HptZzJkzh4yMjIEboehXV7mFqqpkp2X3qUHuCsoq6j6uII42rVbL6NGjOfXUU7nx+utx/+1vdG7ahM9kwm+10qrT4YiPh9mzobERnnwS9iz0E0IIIXr75S9/yR/+8AfefvttduzY0ef5k046iRtvvHEIRjYyHHJIrqys5JxzziE5OZmbbrqJ//znP3R2dg7G2MQ+dJ9JnhUzi+y0/hfpZduymW+bf3QHJ/YrNDSUxMREfn7eeZwaG8ua0lLUkhIM1dUA2O12Ova092PNGqirG8LRCiGE2J8HHniAhx56qN/nHnroIR544IFBvX94eDjXXXcdGo2GzZs393jO7/ezadMmZsyYMahjGMkOOSS//PLL1NTU8Pbbb2O1Wrn11luJiYnhoosu4vXXX6epqWkwxil6cblcmEymoR6GOAyJiYmEazTMyMqiKjSUnJIStOXlwaDsLCjAW1UFUVEwceIQj1YIIcS+aLVa7rvvvj5B+aGHHuK+++5Dq9UO+hi8Xi9ms5ktW7b0OF5cXIzD4ZCQfAQOqyZZo9Fwyimn8Pjjj1NcXMx3333H3LlzeeGFF0hKSuLUU0/liSeeoKqqaqDHK/bo3d1CDC+jMjLQW62cm51NiaqSV16OvrKS0A0bMOzeTavZjCcxUXokCyHEMezee+/lwQcf7BGUuwLygw8+yL333jvoY/j1r39Ne3t7n5C8YcMGdDodU6ZMGfQxjFQD0lsqMzOTu+++m9WrV1NeXs5VV13Fl19+ydtvvz0Qlxe9dC+3EMOTfuxYQmbMIMrhYOFZZ7GlrY2SsjLw+1GVwGbiDTExdMTGDvVQhRBC7Ef3oGw0Go9qQF6/fj3PP/88P/jBD/oNyVlZWfvMCt999x1paWm8++67gz7O4WrAG7DGxcVx3XXX8cEHH3DnnXcO9OXFHjKTPMxpNJh+/GMMSUmkORxcnJVFa3MztZWV6JuawOOh6dRT2b5zJ263e6hHK4QQYj/uvfdeDAYDnZ2dGAyGoxKQ/X4/P/3pT7n55pu58sor2b59Ox6PJ/j8hg0b9ltqMXfuXJ5//nneeeedQR/rcKU73Bfefvvt/R5XFAWTycT48eP5f//v/xEVFXXYgxP965pJ7grJYpjKzMTyy1/iuf9+UsvK0I4axddVVYyPi2NcfDy6piZcHg/bt28nPT0dvV4/1CMWQgjRj4ceeigYkDs7O3nooYcGPSj/8Y9/pKGhgQcffJDy8nI8Hg9FRUVMnjwZgI0bN3LRRRft9xrx8fFs3bp1UMc5nB12SM7Pz2fDhg34fD7S09MB2LZtG1qtloyMDP70pz9x++2389VXX5GVlTVgAxYB3cstZCZ5GKurIyw9nYpZs+hMS2Pre+9xzd//zp+nT+fcdesAaJ85k23btpGeno5Od9j/ZIUQQgyC3jXIXY+BQQvKVVVV3Hvvvbz99ttYLBbGjx+P0Whky5YtTJ48mZ07d9LS0nLARXu//vWvqaiokIm3fTjs77gXXHABUVFRvPLKK8FtqO12O9dffz0nn3wyN9xwA5dddhm33XYbn3zyyYANWARmkruXW4hhTFXRnHEG0TNn0rhtG5fcfTdlisKNb77J69ddx2w10Ofa5XKxbds2JkyYIEFZCCGOEf0t0uv6OJhBeenSpZxzzjn84Ac/AECn05GZmRmsS96wYQMQ6L7RvVbZYDAwYcIEAFauXElkZCQzZsyguLhYFvj147C/2/7+97/n008/DQZkCPTre+CBBzjrrLP4xS9+wX333cdZZ501IAMVPXX/qU9mkoep3FxQFMjOJhQYNWoUlZWV3HnnnYwuL+eVl1+mcdIkTt5zekdHB9u3b2f8+PESlIUQ4hjg8/n6XaTX9djn8w34PT/66CM+//xzCgsLexyfPHlyn5A8b968HuecfPLJfPnllzidTh544AE++eQTHnjgAbZu3SohuR+H/Z3WbrdTV1fXp5Sivr6e1tZWACIiIob9RiMrVqxgxYoVg/If+pGQPskjgKJATk7g8+xs4uPjcbvddHz8MTdPnMgTLhfLli3jmWeeYdasWQA4nU6Ki4uZMGGC1CgLIcQQ299mIYNVanHeeefR3Nzc5/jrr78e/PyRRx7hkUce2ec1Hn74YW688UYSEhKYOHEiBQUFgzLW4e6wu1tccMEFXHvttaxcuZLKykoqKytZuXIl1113HYsWLQJgzZo1wWn94WrJkiUUFBSwdu3aoR5KkN/vl5nkkSA7GxYsCATlvDwAUnbtInbLFhxz5vCjP/2JadOmcfvtt/dYWOFyuSguLh72P4AKIYQ4+rZv305ubi433XQTABMnTpTFe/ugqOqeosdD1N7ezm233cbrr7+O1+sFAjUxV111FU899RQWi4WNGzcCMG3atIEa75BpbW0lPDwcu93eo8RkKDQ3NxMVFcUDDzzAeeedx4wZMyQoD2d5eYGgrNWCz4c/O5vihAScTicdHR0sWbKEsrIyXnjhBcaNGxd8mcFgYNy4cYSEhAzh4IUQ4ti2v+/fLpeLkpISRo8eLb+dPU4cyt/5Yc8kh4aG8uKLL9LY2Eh+fj75+fk0Njby5z//GYvFAgTC8UgIyMcah8MBIP+gR4rs7GBARqtFs2ABY8eORa/XExISwjPPPEN8fDxLliyhsrIy+LLOzk6Ki4tpa2sbwsELIYQQI9NhhWSPx8Ppp5/O9u3bCQ0NZcqUKUyZMoXQ0NCBHp/oh8vlApByi5EiLy8YkPH5IC8Pg8HA2LFj0Wg0WK1Wnn32WSwWCz//+c+pra0NvtTn87F9+3YaGxuH8A0IIYQQI89hhWS9Xs/3338/0GMRB6mjowOQmeQRoavUYsECuPfeHjXKFosFm80GQFRUFH/605/w+XwsWbKkx6INVVUpLS2lpqZmiN6EEEIIMfIcdrnF5ZdfzksvvTSQYxEHyel0AoGZZJlFHsa6B+Ts7MCxXov5IiMjGTVqFAAJCQn86U9/orW1lVtuuYX29vYel6uqqqK8vDzwwO+H0lLYvDnw0e8/eu9LCCGEGAEOuwWc1+vl5Zdf5n//+x8zZ84M1iF3efLJJ494cKJ/XTPJspnIMKeqPQNyl67He9bUJiQk0NnZSX19PWlpaaxYsYIbb7yRW2+9lWeeeabHv736+nooLCR53To0xcXgcoHJBBkZsHgxZGYerXcnhBDDxmH2MBDD0KH8XR92SN6yZUtwu8Nt27b1eE5mNwdX10yyyWSSr/VwNn/+vp/rFZxTUlLo7OzEbrczfvx4nnnmGW655RZ+9rOf8cwzzxAVFQWAqaQE09tv0+5wEJqRgcZqBYcD8vOhogKWLpWgLIQQe3T1m3c6ndIp6DjRlaEOZq+Bww7JOV2bIIijrvfCPTHyKYrCmDFj2LZtGw6HgylTpvDnP/+ZpUuXcv311/Pss8+SlJBARE4Oerud9tGjcXV2EqWq6MLCICsLCgpg1SpITwfNYVdaCSHEiKHVaomIiKCurg4As9ksk08jlKqqOJ1O6urqiIiIQKvVHvA1R7S37ZdffskLL7zArl27+Pvf/86oUaN44403GD16NCeffPKBLyAOS/eFe/KP+fih0WgYN24cxcXFuFwu0tPTeemll1iyZAnXXXcdL917L7bSUtxxcaAoeL1eGhoaiIyMDPxAlZwMhYVQXg57FgQKIcTxLiEhASAYlMXIFhEREfw7P5DDDsnvvfceV1xxBT/5yU/YsGEDbrcbCGxX/bvf/Y5///vfh3tpcQDdyy3E8UWn0zF+/HiKiorweDwkJyfz0ksvsXTpUh6+5x7eyszE1O0fv9/vp6mpifDwcMwWC1RVgfRVFkKIIEVRSExMJC4uDo/HM9TDEYNIr9cf1Axyl8MOyb/97W95/vnnufLKK3nnnXeCx0866SR++9vfHu5lxUHoKrcwGAwyk3wcMhgMjB8/nuLiYnw+HzExMbzwwgs8cfPNfJ2fT1ZoKAndtoNXVZWWlhZ8qorVZAKrdQhHL4QQxyatVntIAUqMfIddmFhcXMypp57a53h4eDgtLS1HMiZxAB0dHRgMBjRSV3rcCgkJYfz48cH/oVutVn753HM4U1PZlZdHwdatPV+gqnhLSmiIi8O3p6WcEEIIIfbtsFNWQkICO3bs6HP8q6++YsyYMUc0KLF/TqdTdtsTWCwWxo0bF/xhyWQ2M/eRRwgbM4bmr79mQ24ufo8HTXs7ISUleCIiqJ0zh505ObjXrZP+yUIIIcR+HHa5xQ033MAvfvELXn75ZRRFYffu3XzzzTfceeed3HvvvQM5RtGL1+s9qNYlYuQLDQ1l3Lhx7NixA7/fj3f8eEY99hisWEFtbi6FDQ1MmjkTR0YGHePGEZGXh6m0FGdnJ0p0NIYpU6R/shBCCNGPww7J99xzD36/n9NPPx2n08mpp56K0Wjkzjvv5JZbbhnIMQ6pFStWsGLFCnw+31APRYh+Wa1Wxo4dy86dOzGvWQOKQvTjj9P4v//xwCOP4Nu8mfvPPpt577+PzuHAkZ6OPyQEd0cHYatXY5b+yUIIIUQfh11uoSgK//d//0dTUxNbtmzh22+/pb6+noceemggxzfklixZQkFBAWvXrh3qofRLyi0EQFhYGOPHj0ej1WJdt47Q/HwmnHUW97z5Js1WK9XLluEvLsadmIg/NBS0WvyhobQkJdFaUoLv/fel9EIIIYTo5oj6JENgpX1WVtZAjEUcJNk+U/QnNDSUpEsvZTdgXbMGgKSZM1n1k5+gf+IJPqurI7SujmlJSXtfpCg4o6LwffsthsJCLBMnDs3ghRBCiGPMIc0kf//99/gPYbZp69ateL3eQx6UOHgykyy6s1gsJP/kJzjnzsW6bh0Jf/kL0Vu2EJeeTsjUqaxZu5acnJwe5UP+kBBUl4uyLVuoqakZwtELIYQQx45DCsnTp0+nsbHxoM8/4YQTKC8vP+RBCSEOX0hICMk/+QlagwHF50M1GvEkJDA7K4vTTzuNkpIS/vnPfwY3pdF0dKAaDPjMZqqqqtixY4f8cCuEEOK4d0jlFqqqcu+992I2mw/q/M7OzsMalDh4MpMs+mP67jv0kZE0trTgVVU0DgfGujrGjhlDWFgY//3vf1m1ahULzzqL5NZWHBkZdO7Zqc9ut1NQUMDo0aOxHsrGI35/YMvrtrbAhiWpqTAQvbwH67pCCCHEfhxSSD711FMpLi4+6PNPOOEEQkJCDnlQQogjkJcHOTlozziDqJNOouqtt9B99BF0dhJSUkJ8XByLzz+frz75hO0ffoh/7lx08+f3CJ4ej4dt27YRHx/PqFGjDvzDWGEhrFwJRUXgcoHJBBkZR95ebrCuK4QQQhzAIYXk3NzcQRqGOBTdF+7JTLLoYU9AZsECyM5GB6Rcfjm7tVr44AOUzk70djtRnZ2cM3cuH2zbxq+//pqz587lMputz39PtbW1tLa2YrPZ9v0bpMJCWL4cGhogJQUsFnA4ID8fjqS93GBdVwghhDgIR9zdQghxDFHVYEDuotFoGJWURFNmJg2JiXSOGoXG6cRvNjM9Lo7L77uP/z71FKWlpSxbtgydruf/Fjo6OigqKiIhIYHExMSeQdrvD8z0NjRAVhZ0PRcWFnhcUACrVkF6+qGVSAzWdYUQQoiDJCF5mJOZZNHD/Pn9HlY0GqIVBeLiKO3WAi5s/XquTksj4corufHNN6mqquJ3v/sdERERPV6vqirV1dW0tLQwevTovWVU5eWBUoiUlL1BNnhTBZKTAzPC5eVgsx38+xis6wohhBAHSaZghDgeZGfDggVEf/89E6qr0Wg0hK5fj3XdOtpmzWLq0qWsWLGC4uJirrzySoqKivq9TEdHB4WFhVRXVwfKftraArXCFkv/97VYAs+3tR3aeAfrukIIIcRBkpA8zMlMsjhoe4Kydd06Jr7/PhH5+bTNmkX7zJkAzJo1i7/+9a+Eh4dz3XXX8dFHH/V7GVVV2b17N0VFRbj0+sBiOoej/3s6HIHnD6VLBgTOH4zrCiGEEAdJQvIwJDvuicOWnQ1aLQatlpj4ePynnNLj6cTERP7yl7+wcOFCHnjgAR577DE8Hk+/l3I6nRS0t9M6ahRqeXmgHro7VYXKysDiutTUQxtnamqgi0VFxcBeVwghhDhIR1yTPHv27H5nM1VVRVEU1uzZHlcMDplJFockLw98PtBq0fp8pNfUUGaz0dTUFDzFaDRy7733MmnSJB5//HGKi4t57LHHiI2N7XM5VVGomDmTxJISwr7/Hp3NtrcLRWUlxMTAokWHvrhOowm0eauoCCzSS04emOsKIYQQB+mIQ/I//vGPgRiHEGKw9WoPR14empwcRms0mDMyqKqqCv6WQlEULrzwQiZMmMDdd9/N5Zdfzq9//WtO6TXzDKBrbMSemUlHfT1RVVWEQKAUYsYMGDUKamsPr1VbbS2ccEIgFBcVQVXVwFxXCCGEOAhHHJLT0tKCn9fU1LB27VoA5syZQ3x8/JFeXhyAzCSLg9I7IMPejzk5xAPmGTPYtWtXjy2pJ02axBtvvMGDDz7Ibbfdxvnnn88dd9xBaGjo3msrCsaqKtpmzqQpKYlQVSUpPR3j7t2B+3b7f8QhURTYtSswzssu27vjXknJkV1XCCGEOAgD9rvKt956i5NPPpl//etffPTRR5xyyim88847A3V5IcSR6Kd/MhBczIeqYrVayczMxNKro0R0dDRPP/00v/71r/n888/50Y9+xLfffht8vn3mTNpmzcK6fj2G6mqaRo2i7IsvsH/4Id5TT+17z4PVNba8PCgrg8mTAx/z8vp/L0IIIcQAGrA+yY899hhr164lMjISgObmZubPn8+Pf/zjgbqF2EN23BOHbB/9k4EeYdNgMJCenk55eTkNDQ3B44qisGjRIubNm8eDDz7IzTffzEUXXcQvfvELzGZzsEOGdd06QvPzUXw+ambNwhUdTVJ9fb/1zAel22w3X3wRqKeWgCyEEOIoGLCZZL/f3+NXsKGhofj9/oG6vBDiKFEUhbS0NNLS0tD0WhiXkJDAihUruOeee/jXv/7FpZdeyvr164HAjLKq1aL4fKhaLe0zZ+L1eikvLw+0i3O5Dm9AezpydC04lIAshBDiaBiwkHz55Zdz4okn8uijj/Loo49y0kknceWVVw7U5cU+yEyyGCwxMTGkp6djMpl6HFcUhYsvvph33nmHuLg4fvrTn/KHP/wBwzffBAOy4vMRuic8AzgcDgoKCqipqel7I78fSkth8+bAx94/XHfryIHPF3gshBBCDLIBK7dYtmwZp59+OqtXrwbgueeeY+aeX8EKIYYns9lMZmYmlZWV1NfX93guOTmZF154gbfffpvNzz7L//71L9Jvuon4H/4wuJsfECzFUFWVqqoq7HY7o0ePxmAwBLaWXrky0L3C5Qp0r8jICLR/y8zstyMHOTmBAciMshBCiEE0YCEZAjt2zZo1ayAvOeRWrFjBihUr8Pl8Qz2UfslMshhsGo2G1NRUwsPDKSsr67G5iEaj4acZGfgWLeL369Zx5+9/z+XV1fz0pz8F6BOUAdrb2ykoKCDN6STyjTegoQFSUvb2Qc7PD/RHPuGEQHeLfXTk6PFYCCGEGGBHXG7x/fffc9FFF3H99dezfft2zjzzTEJDQ5kyZUqPFfDD1ZIlSygoKAi2tjsWyI57YiiEh4eTlZVFREREzydUFe0ZZ/Czt9/m5z//OW+//TZXXHEFa81m2mbN6rtjHuDzeHC++SatJSV4J0yAsLBAOUVYGGRlBYLzl18GQvB+OnIIIYQQg+WIQ/INN9zAT37yE0466SROPPFEbrrpJhobG3n66ae55ZZbBmKMQohjhE6nY+zYsaSmpgYX9bXPmkX7zJnodDquvvpq3njjDfR6PVdffTVPrF1L89Spfa5jqKnBVFqKMyqK+oYG2tvb9/7wpyiBHfbcbhg9uv+BZGfvv2OHEEIIcYSOOCR3dnZy4YUXcs011xAVFcWFF16I0WjktNNO6/FrWTE4pNxCDIXY2FgyMjL6LOoDGDduHK+99hrXXXcdr7zyCldddRVFRUU9ztE4nShuN/6QEFRVpbW1lcbGxr0bmVgsgRrltraj8XaEEEKIPo44JHcPab03IejdPkoIMXKEhISQmZlJdHR0n+d0Oh033ngjr732GqqqctVVV/HHP/4x2AbObzajGo1oOjqCr+ns7KS+vp62tjbU9vbAIj6r9ai9HyGEEKK7I1649/333xMXF4eqqrS0tBAXFwcE6mbtdvsRD1Dsn8wki6Gk0Wiw2WxYrVbKy8v79EbPyMjgjTfe4LXXXuMvf/kLOTk53HfffUybMgWXzYa5uJgOiyVQYkHg/xttra34y8ownnACptTUIxug3w/l5Xu3tE5NBfnhXQghxEE44pAc/PXoHhUVFTz44IO8+OKLR3ppsQ+ycE8ca6Kjo7FYLJSVldHe3t7jOZ1Ox3XXXcfpp5/Ogw8+yA033MCPf/xj7jj3XFJqaggpKcEdF4c/JARNRwfGujrcERFUTp5MSGkpycnJgXZxh+pA7eWEEEKI/RjwKZWmpiZefvnlgb6s2AeZSRbHCpPJRHp6OqmpqWi12j7P22w2XnzxRW699Vbef/99Hr3tNraGheFMT0dvt2MqL0dvt+PIyMAxeTK6xkaam5vZunUrVVVVh9aGsbAQli8PtJOLiYH09MDH/PzA8cLCAXznQgghRqIB7ZMshBCxsbFERERQVlbWp+RKq9Xyk5/8hFNOOYV/3X03ea++yhfz53POTTdhUVX8ZjOG3buxrl9PZ2IiENjyvqamhvr6euLi4oiPj+83hAf5/YEZ5IaGQDu5rh8ku9rLFRTAqlWB4CylF0IIIfZBvkMMczKTLI5Fer2ecePGMXr0aHS6vj+Lp6am8tO33iL+Rz9C++WXPL5kCbkNDRiqq7GuX0/bnrZy3fl8Pqqrq9m8eTO7d+/uU/8cVF4eKLFISdkbkLt0tZcrLAycB6h+Py2lpdRu3kxLaSnqvq4rhBDiuCIzyUKIQRMVFUVYWBjl5eU0Nzf3eE6j0TDnrrtwTppE3J/+RMXNN7N27FjSf/YzNPvZ0r4rLDc1NZGSkkJ4eHjPE9raAjXIvbrtBFksUFUFbW3UFxZStHIlDUVFeF0udCYTMRkZZCxeTKzULQshxHHtkEPyhRdeuN/nW1paDncs4iB1X7gnM8niWKfT6RgzZgzNzc2Ul5f3WexrPuccflBZyfaiIlZ/+y0/f/BBbm9v57zzztvvf99ut5sdO3YQERFBcnIyRqMx8ITVGlik53AESix6czjAZKKxtpbv3nsPZ0MD4Skp6C0WPA4H1fn52CsqmLt0qQRlIYQ4jh1ySO4za9PP81deeeVhD0gIMTJFRkZitVopKyvr8cN06Pr1KH4/4zMzSU1NpXXbNn7zm9/w6aef8qtf/YqEhIT9XrelpQW73U5cXByJiYloU1MDXSzy83vWJENgK+vKStTp0yn49lucDQ3EZmUFw7gxLIzYrCzqCwooWrWKmPR0FKlbFkKI49Ihh+RXXnllMMYhDpPMJIvhpGtb68bGRioqKghZswbrunXBGuTQ9ev5hcnE7NmzWfLuu/zoRz/i1ltvZdGiRfv9b11VVWpra2lsbCQxMZHYRYtQKioCi/SSkwMlFg4HVFZCTAxtc+bQ8PrrhKek9LmuoiiEJSfTUFiIvbycCJttkL8qQgghjkUyRSKEOOqio6OZ1NhI7JYtPRbptc+cSdusWZzodvPJr37FmWeeycMPP8ySJUuoqqo64HW9Xi8VZWVsq6ig7eSTwWYLdLnYtg0aG2HGDFi6lI74eLwuF/p91C0bLBa8Lhdu2RZbCCGOW7JwbxiSmmQxEui0WqIuughlyhTcFRV4PB6AYGAOUVV+/etfc8YZZ/Dwww/zox/9iCuuuILLL78cyz7CramkhIicHELXr0f1eLCPGoVp9GiMp50GU6cGdtz78kss1dXoTCY8DgfGfuqWOx0OdCYTRtkWWwghjlsykyyEGBrz50N2NpGRkUycOJHo6OjgU+0zZ9I+axYA8+bN45133uGSSy7htddeY9GiRbzzzjt9FgCaSkqIffttzMXFeENDUfx+vE4nzq1bafnb33A1N8OXX0JODpb4eGIyMrBXVPTZwVJVVVorK4nJzCT8SLfFFkIIMWxJSB7mZCZZjARarRabzcbYsWPR6/V9nrdYLCxdupT333+fk08+mT/84Q9ceeWVbNmyJXCC309ETg56u52O0aNxjxmDOzkZQ3MzPpMJta6O1gcfpOX99/GcfDLKggVkLF6MOSaG+oICXHY7fq8Xl91OfUEB5pgYMhYtkkV7QghxHJPvAEKIY0ZERARZWVlERkb2+3xCQgL3338/r732GhqNhmuuuYbHHnsMz86dmEpLccfFBbtZdCYm4k5Kwlhdja6ujtAtW2hKSGBzVBSVlZVETZjA3KVLSZw+nY7GRhq3baOjsZHEGTOk/ZsQQgipSR7uZCZZjDRdfZWbmpooLy/H5/P1OScrK4tXX32Vd999l+eff56aTz/lzykphPVqF9eZmIihuhpFqwW/n44xY3p0wkhKSuKkZctorajA3daG0WolPDVVZpCFEELITPJw1LuGUoiRKCoqiqysLCIiIvp9XqfTcdlll/Huu++SlJ7Od5s388V//kNra2vwHEN1NYqqgs8HGg0hu3YFn6v59lvyV65ka0EBHouFuEmTiLDZUDQaSvPyKM3NHeR3KIQQ4lgmIXmYk5lkMZIZDAbGjh3LmDFj+q1VhkAJxl1//CNjzjkHS1MTf3/3XdauXYumogLj7t24ExPxxsXRPmkSpl27CF2/Hgj822lYt47qb76hvLycLVu2UF9fT0luLqU5OT03IRFCCHHckXILIcQxLzIykvDwcKqrq6mtre372xSNhqhrr2WcwUDi999TtmkTJVu2YB03jkSnE09UFPU//jG6pias69YFXrOn1VzDnsfRM2eSv3Ilzfn5pC9cSOoppxzNtyiEEOIYIyF5mJOZZHG80Gg0jBo1iujoaCorK7Hb7T2ed40eTdPll5OYk0Pq11+zu7SU9UVF1NbWEn/jjUy02VBGjw6c7Pdj2L2bUeHhmGw2KteupTE/H9XnI2bWLHwpKWzZsoWEhARiYmLQSI2yEEIcdyQkCyGGFZPJxLhx42hpaaGiooLOzk4AQtetA0Wh5uqrMZx9Nnqnk9DKSla8/TYhjz1G+vvvc8q99zIjKoqInBxM//wnittNgtFIWEMDlamptMXEEL1nhtnj8VBRUUFNTQ3x8fHExsZKWBZCiOOIhORhSHbcEyLQLi4sLCxYgoGiBEspunbts40bx+uhobTExPDM99/zuyuv5ImUFKanpOBJS8MfEkJnWRlR9fVY2tspzsqicf36YFCGQFiurKykpqaGuLg44uLi0Gq1Q/KehRBCHD0SkoUQw1ZXCUZERASlRiNt0CMoh65fj3X9epQLL+TW++/H/qtf0bF6NS9VVjKpvZ3MxEQ8zc2Ys7KIcblQvV7WrV0L0CMoA3i9Xnbv3k1NTQ2xsbHEx8fvczGhEEKI4U9CshBi2LNYLGRmZlIVFhYMyqH5+Sg+H22zZtE+cybm3bsZZzbTsWgRmp07Kd24EcfWrcRnZJCZkIDb6STabid59Ggquy3m683v91NbW0tdXR3R0dEkJCRgNBqP8jsWQggx2CQkD3NSbiFEgEajISUlBcdPfoJj82a8bjeqVhssvdA4nShuN5qEBGbOnMmYqCiKi4v5eutWtlRUMHfWLDI0GqLHjsUVE3PAfuSqqtLQ0EBDQwORkZEkJiYSEhJyNN6qEEKIo0BCshBiRLGsW4c5Kor2jg7a7XZC16+nfeZM/GYzqtGIpqMDf2gokaNHM2/0aCY0NfHdd9/x7eef0xwdjdLYSPzcuYd0z+bmZpqbm4mIiCAxMRGz2TxI704IIcTRIku1hyFZuCfEPuTlQU4OymmnYX38caIvvpjYLVsIXb+ezoQEXDYbxro66PZvKCoqinPOPpvzp05ls9fLBbfcwjPPPIPD4Tjk27e0tFBYWMj27dt77PwnhBBi+JGZZCHEyLAnILNgAWRnA2A480yiDAYsH39MtdFIy4IF6GtqCCkpwR0Xhz8kBE1HB8a6OnTjxrHw7rsp/+orXn75ZT7++GNuvfVWzjrrrEP+YbS1tZXW1lZCQkKIj48nKipKfqAVQohhZsTOJDudTtLS0rjzzjuBwAzPrFmzmDZtGpMmTeLFF18c4hEODPnGK8QeqtojIAdlZ2M8+2xsaWkkLFhA8xVX8P/bu/Owpu50D+DfkECAQBK2sCkgyqYg7shUWRRFUFurrcWl1d7WmXbUdmrpVKcd7XTa2ta5tvWWqbW31Y5L3VrcFxQBxxardcW1alUEEsIiBMIO5/7hkCuuqIEk8P08T56H5Jyc8+b3kJw3b35LVVAQrMvLYZubC+vycuiDg1GUlISmoCC88MIL2LhxI0JDQ/HWW2/h5ZdfxqVLlx4qpOrqaly5cgU5OTlQq9VoaGgwwgslIqL20GErye+//z4GDx5suO/o6Ij9+/fD3t4eer0eoaGhGD9+PFxcXEwYJREZTUzM3bf9J3F2ycyEk0KBgjffhPrkSYj0ejTZ26POwwMOx47BoaQElQMGwNPTE4sWLUJ2djYWLVqEyZMnIykpCTNmzICDg8MDh1ZfX2+YPs7Z2RkqlYqD/IiIzFyHrCRfuHAB586dQ0JCguExsVhsGExTW1sLQRDuO3rdErCSTPQARCJYZWWhy+XL8B82DFa9e6POywsOx47dmF/5lvdTZGQk1q5di5deegnff/89JkyYgE2bNj10RbipqQnFxcU4c+YMfv31V1y/fr1DfA4REXVEZpck79+/H2PHjoWXlxdEIhE2bdp02z4pKSnw8/ODra0tIiIicOjQoRbbk5OTsXDhwtueV1ZWhvDwcHTp0gVvvPEGXF1d2+pltCleVIkeUnT0jS4ZGRmwP3wYISEh6J6XB8XRo4b5lG9lY2OD559/Hhs3bsSAAQPw3nvvYdKkScjMzHyk92JFRQV+++03nDp1ChqNhl0xiIjMjNklyXq9HuHh4UhJSbnj9nXr1mHOnDlYsGABjh49ivDwcMTHx0Or1QIANm/ejMDAQAQGBt72XKVSiRMnTuDy5ctYs2bNjaVsLRwryUQPRhg6FBXBwahcuxbVc+ZAcfQo3CZOhMPo0fd8P3l4eOD999/HypUr4ebmhuTkZLz44os4ceLEI8VTV1eH/Px85OTk4OrVq6iqqnqk4xERkXGYXZ/khISEFt0kbrV48WLMmDEDzz//PABg6dKl2L59O7755hvMnTsXBw8exNq1a7FhwwZUVlaivr4ecrkc8+fPNxzD3d0d4eHh+Pe//42nnnrqjuepra1FbW2t4b45TefESjLRwyk6exbnUlNRfO4cAo4ehVgshp2bG2RJSfD29oarqyvy8vJQVlZ212OEhIQgJSUFP//8M5YsWYIXXngBMTExmDVrFvz8/B46tuauGMXFxZDJZHBzc4OTkxOsrMyulkFE1ClY1KdvXV0djhw5gri4OMNjVlZWiIuLQ3Z2NgBg4cKFuHbtGq5cuYJ//OMfmDFjBubPn4/CwkJUVFQAAMrLy7F//34EBQXd9VwLFy6EQqEw3Lp27dq2L+4hsZJM1DpFZ8/i5yVLoD52DN41NZA5O0Mik0FfUIALc+ei6OxZSKVSdO/eHYGBgZDJZHc9lkgkwuDBg7Fq1Sr8/e9/x6+//oqJEyfi/fffR1FR0SPHqtfrDbNi5OXltfjCTkRE7cOikuTi4mI0NjbC3d29xePu7u7QaDT3fO7Vq1cxdOhQhIeHY+jQoZg9ezbCwsLuuv+8efNQXl5uuF27ds0or4GI2p/Q1IRzqamoKi5GgL09VFotyvz9kT9yJOoGDoTiwgWoFy+G0NQE4MZsOMHBwejRo8c9Z6GwsrJCQkICNm7ciD/96U/Yt28fxo0bhy+//BLV1dWPHHdDQwMKCwtx6tQpXLhwAWVlZfwliYionZhddwtjmj59uuHvQYMG4fjx461+rlQqhVQqNX5QRsZKMtH9lefmovjcOfiJRHC6ehXX/fxQ/p+uEeXdusGuqgryQ4dQuX49HJOSDM9r/iWpuLgY+fn5dx1cZ2Njg8mTJ+Pxxx/H8uXL8e233yI1NRUzZ85EYmIixGLxI7+G5gVKbGxs4OrqCldXV1hbWz/ycYmI6M4sqpLs6uoKsVh824C7wsJCeHh4mCgqIjJ3tRUVaKipgUQqbZEgN9MHBUHj7Iz6uwyac3V1RWhoKNzd3e/5xdTBwQGzZ8/Ghg0b0LdvX/ztb3/DxIkTsWPHDqPNXlFXV4eCggLk5OTg0qVLKC8vN8pxiYioJYtKkm1sbNC/f3+kp6cbHmtqakJ6ejoiIyNNGFn7uvnnVlaSie5P6ugIia0titzcbkuQAaBOr0dZ166wGjbsrscQi8Xo0qULQkND4ezsfM/zeXt7Y+HChfj222/h4+OD+fPnY+LEidi+fbvRkmVBEFBWVoaLFy8iJycHBQUFqK+vN8qxiYjIDJPkyspKHD9+3NA14vLlyzh+/Dhyc3MBAHPmzMFXX32Fb7/9FmfPnsXLL78MvV5vmO3C2FJSUtCzZ08MHDiwTY5PRG1P4eMD1+BglF+7dlufXkEQoMvLg2tICBQ+Pvc9lo2NDbp164aQkBA4Ojrec99evXrhk08+wb/+9S/4+vpiwYIFePrpp7Ft2zajzotcV1cHtVqNnJwc/Pbbb4ZBykRE9PBEgpmNAsnMzERsbOxtj0+bNg0rVqwAAHz++edYtGgRNBoN+vTpgyVLliAiIqJN49LpdFAoFCgvL4dcLm/Tc93PK6+8gtWrVyMtLQ0BAQEmj4fIEjTPblFVXAx5ly6wkclQp9dDl5cHe1dXRLzyCtxCQh74uGVlZa2egeLcuXP46quvkJWVBXd3d0yaNAnjxo17qKWu78fOzg4qlQrOzs6cRo46NXO6fpNlMbsk2VyZ05uMSTLRw7l5nuSGmhpIbG3hGhKC4HHjHipBbiYIAoqKiqBWq1tVIb548SJWrVqFXbt2QSqVYvz48UhKSrpt5h5jkEgkcHFxgZubm0UMRiYyNnO6fpNlYZLcSub0Jrs5SQ4MDLzvT75E9P+EpiaU5+aitqICUkdHKHx8IDJSpbWxsREajQZarRZN/5lO7l60Wi3WrVuH77//HtXV1YiPj8ezzz6LgIAAo8RzK4VCAZVKZfLPMKL2ZE7Xb7IsHXoKuI6K32uIHtyVzExAJIJfdDSUtwzeu5KVBQgC/GJiHukcYrEY3t7ecHd3R2Fh4X2TZZVKhdmzZ+O//uu/sGnTJqxZswY7duzA4MGDMXXqVERERBh1cG7zvO+2traGrhjGmJ6OiKgjYkc1C8fZLYhaSSTClYyMGwnxTa5kZeFKRgZgxPeSRCKBt7c3wsLCoFKp7vs+lclkmDJlCjZv3oz33nsP169fx6xZszBlyhSjTh/XrKamBrm5ucjJyUFubi5qamqMenwioo6ASfJ9cHYLoo7BLzoafrGxLRLl5gTZLzYWftHRRj+nRCJB165dWzVtXPP+o0aNwqpVq/DFF1/A1dUV8+fPxxNPPIGVK1eisrLSqPE1NjaiqKgIp0+fxvnz51FSUtKqbiJERJ0B+yS3kjn1aZo9ezbWrFmDtLQ0BAUFtcnIeKKOqjkxFonFEBob2yxBvhO9Xo9r165Br9e3+jkXL17E6tWrsXPnTsMgv2eeeabNFlASi8VwdnaGq6sr7O3t2+QcRO3JnK7fZFlYSbZw7G5B9GD8oqMNCbJILG63BBm40a0iODgYfn5+kEhaNySkR48eWLBgAbZu3Yqnn34aqampeOKJJ/DXv/4Vv/76q9FjbK4unz17FmfPnkVxcTGry0TUKTFJtkAs/hM9vCtZWYYEWWhsvK2PcntwcXFBaGgo3NzcWv0cNzc3zJo1C9u3b8err76KY8eOYfLkyZg5cyYOHjzYJp8LVVVVuHr1Kk6ePInc3FxUV1cb/RxEROaKs1tYOFaSiVrv1j7IhkF7QLtWlIEb3Rp8fHzg4uKCq1evtjoBlclkmDx5MiZOnIj09HSsXLkSs2bNQkBAAKZOnYqRI0fC2traqLE2V5eLiorg6OgIlUoFpVJp1HMQEZkbVpKJqFO40yC9Ow3ma28ymQwhISHw9PR8oC+9EokE8fHxWLlyJZYuXQqVSoUFCxbgiSeewL/+9S+jD/JrVlFRgUuXLuHUqVMoLCxEY2Njm5yHiMjUWEm2cKwkE7WSINxxkJ7hvgm7MYlEInh5ecHJyQlXr159oIF9IpEIAwYMwIABA3Dp0iWsXr0aX3zxBb7++ms8+eSTSEpKapNBfrW1tcjLy0NBQQFcXFygUqlga2tr9PMQEZkKZ7doJXMaHTtr1ix89913SEtLQ8+ePWFnZ2fSeIjIuIqKipCfn//QVdri4mKsW7cOGzduRFVVFaKiojBhwgQMGjQIVkZaXfBO5HI5VCoVFApFm52D6EGZ0/WbLAu7W9yHOc6TfPP3GlaSiToeNzc39OrV66H7/bq6umLmzJnYvn075syZg9zcXMyaNQtPPvkk1q5d22YD8HQ6HS5evIjTp09Dq9WyKwYRWTRWklvJnL6Jzpw5E2vXrkVaWhp69erFnziJOrCysjJcu3YNdXV1D30MQRBw8uRJbNiwAXv27IGDgwMmTpyIpKSkNq36WllZwcXFBW5ubvzFi0zGnK7fZFlYSbZwrCQTdWxKpRI9e/aESqV66GOIRCKEh4fjvffeQ2pqKhISErBy5Uo8/vjj+OKLL6DT6YwY8f9rampCUVERzpw5g/Pnz6O0tJRTWBKRxWCSTERk5sRiMbp27YqQkJBHXgXPy8sLycnJ2LJlC5588kmsWrUKY8eOxWeffQaNRmOkiG9XWVmJy5cv4+TJk8jPz0dtbW2bnYuIyBiYJFs4VpKJOg97e3uEhITAx8cHYrH4kY7l7OyMP/3pT9i8eTPGjx9vWMnvrbfewpkzZ4wU8e0aGhqg0Whw6tQpXLhwAeXl5W12LiKiR8Ek2QLx50qizs3NzQ2hoaFwcXF55GO5urri1Vdfxfbt2/Haa6/h1KlTeO655zBjxgxkZma26eC75oF+zXMuNzQ0tNm5iIgeFJNkC8dKMlHnJJFI4Ofnh+DgYMhkskc+nkwmQ1JSEn744Qd8/PHHaGpqQnJyMiZMmIB169ahqqrKCFHfWfOcyzk5Obhy5UqbLYRCRPQgmCQTEVkwmUyG4OBg+Pn5GWU5arFYjGHDhuHrr7/GihUrEBISgsWLF2P06NFYsmQJCgsLjRD1nTU1NaGkpATnz5/HyZMnkZubi4qKCv56RkQmwRX37iMlJQUpKSlmM9/nrRcLVpKJCABcXFygVCqh0WhQWFholMQyNDQUCxcuhFqtxrp16/D9999j9erVGDFiBKZMmYKQkBAjRH5n9fX1KCoqQlFREaytraFUKuHs7AwHB4c2OycR0c04T3Irmcs8i4IgYObMmVi3bh3S0tIQHh4OiYTfdYjo/9XW1iI/Px/Xr1836nH1ej22bNmC7777DgUFBejXrx8mT56MqKioNl3J72Y2NjZwdnaGs7Mz516mVjGX6zdZHna3sDCsJBPR/UilUvj7+xutv3IzmUyGSZMm4YcffsBHH32EhoYGJCcnIykpCbt27WqXgXd1dXXQaDQ4c+YMTp8+jYKCgjZbQZCIOjcmyUREHVRzf2V/f39IpVKjHVcikWD48OH45ptv8PXXX8PT0xNvv/02nn76aWzYsKFNB/ndrKamBmq1GmfOnMGpU6eQn58PvV7fLucmoo6PSbKFYyWZiO7HyckJvXr1QpcuXR55fuVbhYeH47PPPsPKlSsRGBiIRYsWYfTo0fj0009x9epVo57rXmpra6HRaHDu3Dnk5OTg2rVrnCWDiB4JO7NamObuFkyOiehBiEQiuLu7w8XFBRqNBlqt1qizRoSEhOCjjz6CWq3G+vXrsWnTJqxatQrh4eF4/PHHERcXZ9SuH/dSV1cHrVYLrVYLiUQCpVIJJycnODo68rOTiFqNA/dayVw6/jc2NmL27NlYv349du/ejX79+vFDn4geWFsN7rv5+JmZmdiyZQsOHToEW1tbxMXFYezYsejbt69JPreYMHdO5nL9JsvDSrKF4XcaIjKG5sF9er0eeXl5Ru+aIJVKER8fj/j4eGg0Gmzbtg1bt27F1q1b0bVrV4wZMwZjxoyBu7u7Uc97Lw0NDSguLkZxcTGsra3h5OQEZ2fndqtwE5FlYZ9kIqJOTCaTISgoCN27d4etrW2bnMPDwwMvvvgiUlNTsXTpUvTu3RvLly/HmDFjMHv2bOzZswe1tbVtcu67qa+vh1arxblz53D69GloNBrU1dW1awxEZN5YSbZw/LmQiIxBqVRCoVCguLgYarUa9fX1Rj+HlZUVBgwYgAEDBuCNN97A3r17sWXLFsybNw9yuRzx8fF44oknEBQU1K6fbTU1NcjPz0d+fj7kcrlhYZb2mvuZiMwTk+T7MPcV94iIjEUkEsHNzc0wuK+wsBBNTU1tci4HBweMGzcO48aNw5UrV7BlyxZs374dGzZsQGBgIJ555hkkJiYaZantB6HT6aDT6WBlZQWFQgEnJycoFAomzESdEAfutZK5dPyvr6/Hq6++ivXr1yMtLQ39+vUzWSxE1LHV19cjPz8fJSUl7XK+hoYGHDx4EN9//z3+/e9/w83NDZMnT8a4cePg6OjYLjHciZWVFeRyOZRKJZRKpdGn0aO2ZS7Xb7I8rCRbGH6nIaL2Ym1tDT8/P7i5uSE3N7fNFwmRSCQYMmQIhgwZgt9++w0rV65ESkoKli5dipEjR2LChAno1atXu3cza2pqQllZGcrKyiASiSCXy+Hk5MSEmaiDY5JswdgfmYjag0wmQ0hICIqLi5Gfn98uy0/7+/tjwYIFmDlzJjZv3oxNmzZh69atCAwMxIQJEzBq1CiTzEohCALKy8tRXl4OkUhk6JLBPsxEHQ/f0URE1Cqurq4IDQ2Fu7t7u31Jd3V1xQsvvIBNmzbh008/hYeHBz766CMkJCTggw8+wLlz59oljjsRBAFlZWW4fPkyTpw4gd9++w3l5eX8xY+og2Al2cLc/OHLSjIRtTexWIwuXbrAzc0N165dQ3l5ebudt7krhkajMVSXf/jhB/Ts2ROjR4/G8OHD4erq2i7x3KqpqQnXr1/H9evXIZFIDHMwOzg4mCQeInp0HLjXSubS8b+2thavvfYa1q9fj71796JPnz4mi4WISKfTIS8vD9XV1e1+7oaGBhw4cACpqak4ePAgBEFAv379MGLECAwfPhxKpbLdY7qVjY2NYZU/JsymYS7Xb7I8rCRbGFaSicicyOVyQ3/lgoKCdumv3EwikSAmJgYxMTEoKytDZmYm0tLS8NFHH+Hjjz/GwIEDMWLECMTGxposOaqrq4NWq4VWq4W1tTUUCgWUSiXkcjk/w4nMHCvJrWQu30RramowZ84crF+/Hunp6QgPDzdZLEREN2tsbIRarYZWqzVpv9zS0lLs27cPaWlpOHbsGMRiMQYPHoyRI0ciKirKLCq6YrEYjo6OhkVcJBLWrNqKuVy/yfLwXWnBWIUgInNyc3/lvLw8lJWVmSQOZ2dnPPXUU3jqqadQVFSE9PR0pKWlYf78+bCxscFjjz2GESNGYOjQobCzszNJjI2NjYZp5YAbM4goFAooFArY29ubJCYiaolJsoVh4Z+IzJ1UKkX37t1RUVGBvLy8Np9f+V7c3NyQlJSEpKQkaDQa7N27F2lpafjLX/4CW1tbDB06FMOHD8djjz1msoQZAPR6PfR6PQoKCmBjY2NYvMTR0ZFTyxGZCJNkIiJqE46OjggJCUFpaSny8/NRV1dn0ng8PDwwdepUTJ06FXl5eYaEee7cuZBKpfjd736H4cOHY8iQISbtklFXV4fi4mIUFxfDysoKDg4OhiqzVCo1WVxEnQ2T5PtISUlBSkoKGhsbTR0KAA7cIyLL4+zsDKVSCa1WC41GYxafp126dMH06dMxffp05OXlYd++fUhPT8fbb78Na2trQx/m6Ohok3Z/aGpqgk6ng06nw7Vr12Bra2tImB0cHHgdIGpDHLjXSubS8b+qqgrJyclYv349MjIyEBYWZrJYiIgeVENDA9RqNYqKisyy+5hGo8G+ffuwd+9enDx5Era2toiOjkZCQgIGDx5sVgPsxGIx5HK5IWk2p9jMiblcv8ny8B1lYVhJJiJLJpFI0LVrV6hUKhQUFKC0tNTUIbXg4eGByZMnY/LkySgoKMDu3buxc+dO7N69GwqFAiNGjMCoUaPQu3dvk/cVbmxsNCxgAtwY/Nc8JzO7ZRA9OibJRETU7qRSKbp16wZ3d3fk5eWhoqLC1CHdxsvLC88//zymT5+OixcvGpLljRs3wtPTE/Hx8Rg+fDiCg4PNomjRPPgvPz8f9vb2cHJyglKphK2tralDI7JITJItmDl8KBMRPQp7e3sEBgaivLwceXl5qKmpMXVItxGJRAgICEBAQABmzZqF48ePY9euXUhNTcWKFSvg7u6OqKgoxMTEoH///mbR7aGqqgpVVVXIz89nP2aih2T6dzI9EHPsw0dE9KgUCgXkcrlJVu57EFZWVujXrx/69euHP//5zzh27BgyMzORlZWFDRs2wMHBAUOGDEF0dDQiIyPNYuGSmpoa1NTUoLCw0LCISXPSbG1tberwiMwWk2QLxmoAEXUkIpEIbm5ucHZ2hkajgVarRVNTk6nDuiuJRIKBAwdi4MCBSE5Oxvnz55GVlYWsrCzs2rUL1tbWGDhwIKKjoxEdHQ1XV1dTh3zbIib29vaGhFkmk5k2OCIzw9ktWslcRsdWVFTgzTffxPr167F//3707NnTZLEQEbWluro6FBQUoKSkxNShPLCCggJDwnzs2DE0NjYiNDQU0dHRiImJgZ+fn9kVOqytrQ2zZcjlcojFYlOHZBTmcv0my8MkuZXM5U12c5L873//GyEhISaLhYioPVRXV6OgoMBky1w/qvLychw4cABZWVnIzs5GdXU1fHx8DBXmsLAws0tIRSJRi24Zljxbhrlcv8nysLsFERGZNTs7O3Tv3h2VlZXIy8uDXq83dUgPRKFQYPTo0Rg9ejRqa2tx+PBhZGZmYvv27Vi5ciWcnZ0RGxuL+Ph49OnTx+RTywE3xr/cvIiJnZ0dFAoFlEolu2VQp8Ek2cJwnmQi6qwcHBwQHBxsNstcPwypVIohQ4ZgyJAhaGxsxKlTp5CZmYk9e/bg+++/h0qlQlxcHEaOHIlevXqZzed8dXU1qqurodFoYG1tbUiYHR0dzSKpJ2oLTJKJiMiimOMy1w9DLBYjPDwc4eHhmD17NnJycpCWlobdu3djzZo18Pb2xrBhwxATE4OwsDCzSUbr6+tRXFyM4uJiWFlZGfoxK5VKs5j+jshY+N9sYVhJJiK6MRWbh4cHXF1dUVBQgOLiYoueItPKysqQML/22ms4evQo0tLSsG3bNqxcuRIuLi6IiopCbGwsBgwYABsbG1OHDABoamoyzJZx9epVw6p/CoUCdnZ2pg6P6JEwSSYiIoslkUjg4+MDlUqF/Px8ix3cdzOJRIJBgwZh0KBBmDdvHnJycpCZmYmMjAykpqZCJpMhKioKcXFxiIyMNJuEGWi56p9UKjVUmLmICVkiJskWjB84REQ32NraWvTgvrsRi8Xo06cP+vTpg1dffRWXLl1CRkYG9u7di507d0Imk2HkyJGYMGECgoODTR1uC7W1tdBqtdBqtRCLxYaZMhQKhdnN5kF0J0yS7yMlJQUpKSlm0+fNkn9OJCJqa82D+65fv478/HzU1taaOiSjEYlE6NGjB3r06IEZM2bg0qVLSEtLw9atW5GamoqQkBDEx8cjNjYW3t7epg63hcbGRpSWlqK0tBQikQgODg6GbhmWPL0cdWycJ7mVzGWexbKyMvzlL3/B+vXrkZ2djYCAAJPFQkRkzgRBQFFREdRqtdkuc20MDQ0N+Omnn7Bp0yYcPHgQdXV1CAwMRExMDGJjY9GjRw+z/uWxraeXM5frN1keVpItDL/TEBG1jkgkgkqlgouLi0Usc/2wJBIJoqKiEBUVBb1ej+zsbGRkZGD16tVYtmwZvL29DQmzOS5ccqfp5ZpX/TOXGT2oc2KSbMHMuTJARGQuxGIxvL29oVKpDMtcd9SCg0wmQ1xcHOLi4lBfX49ffvkFGRkZ2LVrF1avXg1nZ2cMHz4ciYmJCA0NNbvryK3Tyzk6Ohq6ZVhbW5s6POpkmCQTEVGnYG1tDV9fX0Oy3BFmwrgXa2trREZGIjIyEnPnzkVOTg4yMjKwe/dubNiwAV27dkVCQgISExPRpUsXU4d7m6amJpSXl6O8vBwAYG9vb0iY7e3tTRwddQZMki0M50kmIno0Ny9znZ+fj8rKSlOH1OZunod59uzZOHLkCHbs2IFVq1Zh2bJl6N27NxITExEXFwelUmnqcO+oqqoKVVVVKCgogI2NDRQKBZycnDi9HLUZDtxrJXPp+F9aWoq3334b69evx88//4zu3bubLBYioo5Ap9OhoKCgw0wb9yBqamqQmZmJnTt34uDBgxCJRHjssceQkJCAoUOHWsTMExKJBEqlEs7OznB0dLxtu7lcv8nysJJsYfidhojIuORyOeRyOXQ6HdRqdaeoLDeztbXFqFGjMGrUKJSUlGDPnj3YsWMH5s6dCwcHB8TFxSExMRF9+vQx20F0DQ0Nhn7M1tbWcHJy4gImZBRMki0Y3/xERMbTnCxXVlZCrVZDp9OZOqR25eLigqSkJCQlJeHKlSvYsWMHdu7ciU2bNsHT0xOjRo1CYmIiunXrZupQ76q+vv62BUx4raSHxSTZwrCSTETUthwcHBAQEICqqipoNBqUlZV1us9ePz8//PGPf8RLL72EEydOYMeOHdi4cSOWL1+OkJAQJCQkIC4uDiqVytSh3lXzAiad6ZcBMi4myRaM346JiNqOvb09/P39UVtbi8LCQpSUlHTIeZbvxcrKCn379kXfvn2RnJyMH3/8ETt37sSSJUuwePFi9OrVyzAHs5+fn6nDJTIqJslERET3IJVK4ePjAy8vL8NP+Y2NjaYOq91JpVIMGzYMw4YNg06nw4EDB5CZmYmvv/4aKSkp8PX1RWxsLGJiYtCzZ0+z7cNM1FpMki0Mp4AjIjINiUQCLy8veHh4oKioCFqtFnV1daYOyyTkcjkSExORmJiImpoaHDp0CBkZGUhNTcWKFSvg5uaG6OhoxMTEYMCAAZBImG6Q5eF/LRER0QOwsrKCu7s7VCoVSktLUVhYiOrqalOHZTK2traGZbEbGhpw4sQJZGZmIjMzExs3boSDgwOGDh2KmJgYREZGciEQshhMki0MK8lEROZBJBLBxcUFLi4u0Ol00Gg0qKioMHVYJiWRSNC/f3/0798fc+bMwa+//mpImHfu3AkbGxtEREQgJiYGUVFRcHJyMnXIRHfFJJmIiOgRNU8fp9frDTNidHYikQhBQUEICgrCH/7wB+Tl5RkS5vfeew8ikQh9+vQxdMvw9vY2dchELTBJtmCsJBMRmReZTIbu3bujuroaBQUFTJZv0qVLF0ydOhVTp05FaWkp9u/fj8zMTHz++ef45JNPEBgYiJiYGMTExCAgIIDXODI5JskWprPN1UlEZIns7OzQvXt3VFVVoaCgAOXl5aYOyaw4Oztj3LhxGDduHPR6PbKzs5GRkYHVq1dj2bJl8Pb2RnR0NGJjY9G7d2+IxWJTh0ydEJNkC8Zv2URE5s3e3h49evSAXq9Hfn5+p++zfCcymQxxcXGIi4tDfX09fvnlF2RmZmL37t1Ys2YNnJycEBUVhZiYGAwaNAhSqdTUIVMnwSTZwrCSTERkeWQyGQIDA1FZWYmCggImy3dhbW2NyMhIREZG4s0338Tp06eRkZGBzMxMbN68GXZ2doiOjkZiYiIGDRrEqeWoTfG/y4KxkkxEZFkcHBwQGBiIiooKFBQUcMnke7CyskJYWBjCwsIwe/ZsXL58Gfv27cPu3buxa9cuuLi4YOTIkRg9ejSCgoJ4TSSjY5JMRETUzhwdHREUFASdTge1Ws1k+T5EIhH8/f3h7++PF154AefPn8f27duRlpaG7777Dt26dUNiYiJGjRoFT09PU4dLHQST5PtISUlBSkqK2SxBynmSiYg6juap43Q6HQoKCqDX600dktkTiUQIDg5GcHAwXn31VRw6dAg7duzA//7v/yIlJQX9+vVDYmIihg8fDkdHR1OHSxZMJLCTa6vodDooFAqUl5dDLpebLA61Wo2///3vWL9+PU6cOMF5JYmIOhAmyw9Pr9cjIyMDO3fuxOHDhyGRSBAVFYXY2Fi89dZbJr9+k+VhJdnCsJJMRNRxsbL88GQyGcaMGYMxY8agqKgIu3fvxo4dO/DWW2+ZOjSyUFamDoCIiIhaksvlCA4ORkBAABwcHEwdjsVxc3PD1KlTsWbNGnzzzTemDocsFCvJFoyVZCKijq25sszZMB6ev7+/qUMgC8VKsoVhF3Iios6neTaMwMBADkYjaiesJFswVpKJiDoXR0dHODo6orKyEmq1GjqdztQhEXVYTJItDCvJRETk4OCAgIAA6PV6qNVqlJeXmzokog6HSbIFYyWZiKhzk8lk6NGjB6qqqqDRaHD9+nVTh0TUYTBJJiIisnD29vbw9/dHTU0N1Go1rl+/zl8eiR4RB+5ZGH7oERHR3dja2qJbt27o1asXXF1d+Ysj0SNgkmzB+OFHRER3IpVK4evri9DQUKhUKlhZ8XJP9KD4rrEwrCQTEVFr2djYoGvXrggLC4OHhwfEYrGpQyKyGOyTbMFYSSYiotaQSCTw9vaGh4cHtFottFotGhoaTB0WkVljkmxhWEkmIqKHJRaL4enpCXd3dxQXF6OwsBB1dXWmDovILDFJtmCsJBMR0cOwsrKCSqWCm5sbSktLodFoUFNTY+qwiMwKk2QiIqJOSiQSwcXFBS4uLigrK4NarUZVVZWpwyIyC0ySLczN3S1YSSYiImNRKpVQKpWoqKiARqPhktfU6TFJJiIiIgNHR0c4OjpyFT/q9JgkWxhWkomIqD00r+JXW1uLwsJCFBcXc/A4dSpMkomIiOiupFIpfHx84OnpCa1Wi6KiIjQ2Npo6LKI2xyTZgrGSTERE7cXa2tow13JxcTG0Wi2nj6MOjUmyheFPXUREZEpisRju7u5QqVQoLS1FYWEhqqurTR0WkdExSbZgrCQTEZGp3Dx9nE6ng0ajQUVFhanDIjIaJskWhpVkIiIyN3K5HHK5HFVVVSgsLMT169d5vSKLZ2XqAOjhsZJMRETmxN7eHt26dUNoaChUKhWsrJhmkOViJZmIiIiMysbGBl27doWXlxeKioqg1WpRX19v6rCIHgiTZAvDeZKJiMhSiMVieHh4wN3dnYP8yOIwSSYiIqI2desgP61Wi/LyclOHRXRPTJItDCvJRERkyZoH+dXU1KCwsBClpaVoamoydVhEt2GSTERERO3O1tYWvr6+8Pb2RlFREYqKithvmcxKhx12WlVVBV9fXyQnJwMArl27hpiYGPTs2RO9e/fGhg0bTBzho2MlmYiILJ1EIoGnpyfCwsLg5+cHOzs7U4dEBKADV5Lff/99DB482HBfIpHg008/RZ8+faDRaNC/f38kJiZCJpOZMMoHx3kniYioI7q533JFRQUKCwvZb5lMqkMmyRcuXMC5c+cwduxYnDp1CgDg6ekJT09PAICHhwdcXV1RWlpqcUkyERFRR+fo6AhHR0fU1NRAq9WipKSE/Zap3Zldd4v9+/dj7Nix8PLygkgkwqZNm27bJyUlBX5+frC1tUVERAQOHTrUYntycjIWLlx413McOXIEjY2N6Nq1q7HDb3McuEdERJ2Fra0tfHx8EBYWBi8vL0gkHbK2R2bK7JJkvV6P8PBwpKSk3HH7unXrMGfOHCxYsABHjx5FeHg44uPjodVqAQCbN29GYGAgAgMD7/j80tJSPPfcc1i2bFmbvQYiIiIynpv7Lfv4+MDGxsbUIVEnYHZfyRISEpCQkHDX7YsXL8aMGTPw/PPPAwCWLl2K7du345tvvsHcuXNx8OBBrF27Fhs2bEBlZSXq6+shl8sxf/581NbWYty4cZg7dy5+97vf3TOO2tpa1NbWGu4394vS6XRGeJUPr7KyEnV1dWhqakJFRQUaGxtNGg8REVF7kkql8PHxwfXr11FYWIiampp77q/X6wFwTA89BMGMARBSU1MN92trawWxWNziMUEQhOeee054/PHHb3v+8uXLhddff10QBEFoamoSkpKShAULFrTq3AsWLBAA8MYbb7zxxhtvHeB26dKlh01HqJMyu0ryvRQXF6OxsRHu7u4tHnd3d8e5c+fu+dwff/wR69atQ+/evQ39nFeuXImwsLA77j9v3jzMmTPHcL+srAy+vr7Izc2FQqF4tBfSiel0OnTt2hXXrl2DXC43dTgWjW1pPGxL42A7Gg/b0njKy8vh4+MDZ2dnU4dCFsaikuQHNX36dMPfQ4YMeaCRsVKpFFKp9LbHFQoFP7CMoHnFJXp0bEvjYVsaB9vReNiWxmNlZXbDsMjMWdR/jKurK8RiMQoLC1s8XlhYCA8PDxNFRUREREQdjUUlyTY2Nujfvz/S09MNjzU1NSE9PR2RkZEmjIyIiIiIOhKz625RWVmJixcvGu5fvnwZx48fh7OzM3x8fDBnzhxMmzYNAwYMwKBBg/Dpp59Cr9cbZrtoK1KpFAsWLLhjFwxqPbaj8bAtjYdtaRxsR+NhWxoP25IelkgQzGtOlMzMTMTGxt72+LRp07BixQoAwOeff45FixZBo9GgT58+WLJkCSIiIto5UiIiIiLqqMwuSSYiIiIiMjWL6pNMRERERNQemCQTEREREd2CSTIRERER0S2YJLeCn58fRCJRi9uHH37YYp+TJ09i6NChsLW1RdeuXfHxxx+bKFrzlpKSAj8/P9ja2iIiIgKHDh0ydUhm7Z133rntfy84ONiwvaamBjNnzoSLiwscHBwwYcKE2+YR76z279+PsWPHwsvLCyKRyLDSZjNBEDB//nx4enrCzs4OcXFxuHDhQot9SktLMWXKFMjlciiVSrzwwguorKxsx1dhHu7XltOnT7/t/3TUqFEt9mFbAgsXLsTAgQPh6OgIlUqFcePG4fz58y32ac17Ojc3F6NHj4a9vT1UKhXeeOMNNDQ0tOdLMbnWtGVMTMxt/5cvvfRSi33YlnQvTJJb6d1334VarTbcZs+ebdim0+kwcuRI+Pr64siRI1i0aBHeeecdLFu2zIQRm59169Zhzpw5WLBgAY4ePYrw8HDEx8dDq9WaOjSz1qtXrxb/ewcOHDBse+2117B161Zs2LABWVlZKCgowPjx400YrfnQ6/UIDw9HSkrKHbd//PHHWLJkCZYuXYqff/4ZMpkM8fHxqKmpMewzZcoUnD59Gnv27MG2bduwf/9+/P73v2+vl2A27teWADBq1KgW/6ffffddi+1sSyArKwszZ87EwYMHsWfPHtTX12PkyJHQ6/WGfe73nm5sbMTo0aNRV1eHn376Cd9++y1WrFiB+fPnm+IlmUxr2hIAZsyY0eL/8uYCFtuS7kug+/L19RU++eSTu27/5z//KTg5OQm1tbWGx958800hKCioHaKzHIMGDRJmzpxpuN/Y2Ch4eXkJCxcuNGFU5m3BggVCeHj4HbeVlZUJ1tbWwoYNGwyPnT17VgAgZGdnt1OElgGAkJqaarjf1NQkeHh4CIsWLTI8VlZWJkilUuG7774TBEEQzpw5IwAQDh8+bNhn586dgkgkEvLz89stdnNza1sKgiBMmzZNeOKJJ+76HLblnWm1WgGAkJWVJQhC697TO3bsEKysrASNRmPY54svvhDkcnmLa1Bnc2tbCoIgREdHC6+++updn8O2pPthJbmVPvzwQ7i4uKBv375YtGhRi59jsrOzERUVBRsbG8Nj8fHxOH/+PK5fv26KcM1OXV0djhw5gri4OMNjVlZWiIuLQ3Z2tgkjM38XLlyAl5cX/P39MWXKFOTm5gIAjhw5gvr6+hZtGhwcDB8fH7bpfVy+fBkajaZF2ykUCkRERBjaLjs7G0qlEgMGDDDsExcXBysrK/z888/tHrO5y8zMhEqlQlBQEF5++WWUlJQYtrEt76y8vBwA4OzsDKB17+ns7GyEhYXB3d3dsE98fDx0Oh1Onz7djtGbl1vbstnq1avh6uqK0NBQzJs3D1VVVYZtbEu6H7Nbcc8cvfLKK+jXrx+cnZ3x008/Yd68eVCr1Vi8eDEAQKPRoFu3bi2e0/ym02g0cHJyaveYzU1xcTEaGxtbfBgBN9rp3LlzJorK/EVERGDFihUICgqCWq3G3/72NwwdOhSnTp2CRqOBjY0NlEpli+e4u7tDo9GYJmAL0dw+d/p/bN6m0WigUqlabJdIJHB2dmb73mLUqFEYP348unXrhkuXLuEvf/kLEhISkJ2dDbFYzLa8g6amJvzpT3/CY489htDQUABo1Xtao9Hc8f+2eVtndKe2BIDJkyfD19cXXl5eOHnyJN58802cP38eP/zwAwC2Jd1fp02S586di48++uie+5w9exbBwcGYM2eO4bHevXvDxsYGf/jDH7Bw4UIuc0ltKiEhwfB37969ERERAV9fX6xfvx52dnYmjIzo/yUlJRn+DgsLQ+/evdG9e3dkZmZi+PDhJozMfM2cOROnTp1qMcaAHs7d2vLmPu9hYWHw9PTE8OHDcenSJXTv3r29wyQL1Gm7W7z++us4e/bsPW/+/v53fG5ERAQaGhpw5coVAICHh8dto4+b73t4eLTp67AUrq6uEIvFd2wntlHrKZVKBAYG4uLFi/Dw8EBdXR3Kyspa7MM2vb/m9rnX/6OHh8dtg0obGhpQWlrK9r0Pf39/uLq64uLFiwDYlreaNWsWtm3bhoyMDHTp0sXweGve07zetHS3tryTiIgIAGjxf8m2pHvptEmym5sbgoOD73m7uY/xzY4fPw4rKyvDz4eRkZHYv38/6uvrDfvs2bMHQUFB7GrxHzY2Nujfvz/S09MNjzU1NSE9PR2RkZEmjMyyVFZW4tKlS/D09ET//v1hbW3dok3Pnz+P3Nxctul9dOvWDR4eHi3aTqfT4eeffza0XWRkJMrKynDkyBHDPvv27UNTU5PhYkt3lpeXh5KSEnh6egJgWzYTBAGzZs1Camoq9u3bd1s3vda8pyMjI5GTk9PiS8eePXsgl8vRs2fP9nkhZuB+bXknx48fB4AW/5dsS7onU48cNHc//fST8MknnwjHjx8XLl26JKxatUpwc3MTnnvuOcM+ZWVlgru7u/Dss88Kp06dEtauXSvY29sLX375pQkjNz9r164VpFKpsGLFCuHMmTPC73//e0GpVLYYWUwtvf7660JmZqZw+fJl4ccffxTi4uIEV1dXQavVCoIgCC+99JLg4+Mj7Nu3T/jll1+EyMhIITIy0sRRm4eKigrh2LFjwrFjxwQAwuLFi4Vjx44JV69eFQRBED788ENBqVQKmzdvFk6ePCk88cQTQrdu3YTq6mrDMUaNGiX07dtX+Pnnn4UDBw4IAQEBwqRJk0z1kkzmXm1ZUVEhJCcnC9nZ2cLly5eFvXv3Cv369RMCAgKEmpoawzHYloLw8ssvCwqFQsjMzBTUarXhVlVVZdjnfu/phoYGITQ0VBg5cqRw/PhxYdeuXYKbm5swb948U7wkk7lfW168eFF49913hV9++UW4fPmysHnzZsHf31+IiooyHINtSffDJPk+jhw5IkRERAgKhUKwtbUVQkJChA8++KDFh78gCMKJEyeEIUOGCFKpVPD29hY+/PBDE0Vs3v7nf/5H8PHxEWxsbIRBgwYJBw8eNHVIZu2ZZ54RPD09BRsbG8Hb21t45plnhIsXLxq2V1dXC3/84x8FJycnwd7eXnjyyScFtVptwojNR0ZGhgDgttu0adMEQbgxDdxf//pXwd3dXZBKpcLw4cOF8+fPtzhGSUmJMGnSJMHBwUGQy+XC888/L1RUVJjg1ZjWvdqyqqpKGDlypODm5iZYW1sLvr6+wowZM2778su2FO7YhgCE5cuXG/ZpzXv6ypUrQkJCgmBnZye4uroKr7/+ulBfX9/Or8a07teWubm5QlRUlODs7CxIpVKhR48ewhtvvCGUl5e3OA7bku5FJAiC0H51ayIiIiIi89dp+yQTEREREd0Nk2QiIiIiolswSSYiIiIiugWTZCIiIiKiWzBJJiIiIiK6BZNkIiIiIqJbMEkmIiIiIroFk2QiIiIiolswSSaidrdixQoolUpTh3EbkUgEkUhklrEBwPTp0w0xbtq0ydThEBF1aEySicjobk7mbGxs0KNHD7z77rtoaGgwdWj3tXz5cvz666+G+ytWrDC8FrFYDCcnJ0RERODdd99FeXl5u8b22WefQa1Wt+s5iYg6KybJRNQmRo0aBbVajQsXLuD111/HO++8g0WLFpk6rPtSKpVQqVQtHpPL5VCr1cjLy8NPP/2E3//+9/jXv/6FPn36oKCgoN1iUygU8PDwaLfzERF1ZkySiahNSKVSeHh4wNfXFy+//DLi4uKwZcuWFvvs3r0bISEhcHBwMCTVzQ4fPowRI0bA1dUVCoUC0dHROHr0qGG7IAh455134OPjA6lUCi8vL7zyyiuG7bW1tUhOToa3tzdkMhkiIiKQmZn5UK9FJBLBw8MDnp6eCAkJwQsvvICffvoJlZWV+POf/2zYb9euXRgyZAiUSiVcXFwwZswYXLp0ybB92LBhmDVrVotjFxUVwcbGBunp6QCAf/7znwgICICtrS3c3d3x1FNPPVTMRET0aJgkE1G7sLOzQ11dneF+VVUV/vGPf2DlypXYv38/cnNzkZycbNheUVGBadOm4cCBAzh48CACAgKQmJiIiooKAMD333+PTz75BF9++SUuXLiATZs2ISwszPD8WbNmITs7G2vXrsXJkyfx9NNPY9SoUbhw4YJRXo9KpcKUKVOwZcsWNDY2AgD0ej3mzJmDX375Benp6bCyssKTTz6JpqYmAMCLL76INWvWoLa21nCcVatWwdvbG8OGDcMvv/yCV155Be+++y7Onz+PXbt2ISoqyijxEhHRg5GYOgAi6tgEQUB6ejp2796N2bNnGx6vr6/H0qVL0b17dwA3ktp3333XsH3YsGEtjrNs2TIolUpkZWVhzJgxyM3NhYeHB+Li4mBtbQ0fHx8MGjQIAJCbm4vly5cjNzcXXl5eAIDk5GTs2rULy5cvxwcffGCU1xYcHIyKigqUlJRApVJhwoQJLbZ/8803cHNzw5kzZxAaGorx48dj1qxZ2Lx5MyZOnAjgRp/n5j7cubm5kMlkGDNmDBwdHeHr64u+ffsaJVYiInowrCQTUZvYtm0bHBwcYGtri4SEBDzzzDN45513DNvt7e0NCTIAeHp6QqvVGu4XFhZixowZCAgIgEKhgFwuR2VlJXJzcwEATz/9NKqrq+Hv748ZM2YgNTXVMDAwJycHjY2NCAwMhIODg+GWlZXVovvDoxIEAcCN7hgAcOHCBUyaNAn+/v6Qy+Xw8/MDAEPMtra2ePbZZ/HNN98AAI4ePYpTp05h+vTpAIARI0bA19cX/v7+ePbZZ7F69WpUVVUZLV4iImo9VpKJqE3Exsbiiy++gI2NDby8vCCRtPy4sba2bnFfJBIZkk4AmDZtGkpKSvDZZ5/B19cXUqkUkZGRhi4bXbt2xfnz57F3717s2bMHf/zjH7Fo0SJkZWWhsrISYrEYR44cgVgsbnEeBwcHo73Gs2fPQi6Xw8XFBQAwduxY+Pr64quvvoKXlxeampoQGhraopvJiy++iD59+iAvLw/Lly/HsGHD4OvrCwBwdHTE0aNHkZmZibS0NMyfPx/vvPMODh8+bLbT0hERdVRMkomoTchkMvTo0eOhn//jjz/in//8JxITEwEA165dQ3FxcYt97OzsMHbsWIwdOxYzZ85EcHAwcnJy0LdvXzQ2NkKr1WLo0KGP9DruRqvVYs2aNRg3bhysrKxQUlKC8+fP46uvvjKc88CBA7c9LywsDAMGDMBXX32FNWvW4PPPP2+xXSKRIC4uDnFxcViwYAGUSiX27duH8ePHt8nrICKiO2OSTERmKSAgACtXrsSAAQOg0+nwxhtvwM7OzrB9xYoVaGxsREREBOzt7bFq1SrY2dnB19cXLi4umDJlCp577jn893//N/r27YuioiKkp6ejd+/eGD169APFIggCNBoNBEFAWVkZsrOz8cEHH0ChUODDDz8EADg5OcHFxQXLli2Dp6cncnNzMXfu3Dse78UXX8SsWbMgk8nw5JNPGh7ftm0bfvvtN0RFRcHJyQk7duxAU1MTgoKCHqIFiYjoUbBPMhGZpa+//hrXr19Hv3798Oyzz+KVV15pMX+xUqnEV199hcceewy9e/fG3r17sXXrVkPXh+XLl+O5557D66+/jqCgIIwbNw6HDx+Gj4/PA8ei0+ng6ekJb29vREZG4ssvv8S0adNw7NgxeHp6AgCsrKywdu1aHDlyBKGhoXjttdfuOi/0pEmTIJFIMGnSJNja2rZ4TT/88AOGDRuGkJAQLF26FN999x169er1wDETEdGjEQk3dwIkIurERCIRUlNTMW7cuDY9z5UrV9C9e3ccPnwY/fr1e+Dnt1ecRESdGSvJREQ3mTRpErp06dImx66vr4dGo8Hbb7+NwYMHP3CC/NJLLxl14CEREd0dK8lERP9x8eJFAIBYLEa3bt2MfvzMzEzExsYiMDAQGzdubLH4SWtotVrodDoAN6bMk8lkRo+RiIhuYJJMRERERHQLdrcgIiIiIroFk2QiIiIiolswSSYiIiIiugWTZCIiIiKiWzBJJiIiIiK6BZNkIiIiIqJbMEkmIiIiIroFk2QiIiIiolv8H38BbWDLd5ZqAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "from slsne.utils import plot_colors\n", "import matplotlib.pyplot as plt\n", "\n", "bands = np.unique(phot['Filter'])\n", "for i in range(len(bands)):\n", " band = bands[i]\n", "\n", " use = phot['Filter'] == band\n", " plt.errorbar(phot[use]['Phase'], L_lambda[use].value / bol_scaling[use], fmt = 'o', color = plot_colors(band), alpha = 0.5, label = band)\n", "\n", " plt.errorbar(phot[use]['Phase'], L_lambda[use].value * lambda_AA[use], fmt = 'x', color = plot_colors(band), alpha = 0.5)\n", "\n", "# Plot Bolometric\n", "plt.fill_between(bolo['Phase'], bolo['L_low'], bolo['L_high'], color = 'k', linewidth = 0, alpha = 0.2, label = 'Bolometric')\n", "plt.plot(bolo['Phase'], bolo['L_med'], color = 'k', linewidth = 1)\n", "\n", "# Plot Parameters\n", "plt.errorbar([], [], fmt = 'x', color = 'k', label = r'$\\lambda L_{\\lambda}$')\n", "plt.legend(loc='upper left', bbox_to_anchor=(1, 1))\n", "plt.yscale('log')\n", "plt.xlabel('Phase [Days]')\n", "plt.ylabel(r'L$_{\\rm Bol}$ [erg/s]')\n", "plt.ylim(1e42, 1e45)\n", "plt.xlim(-50, 250)\n", "plt.show();" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.4" } }, "nbformat": 4, "nbformat_minor": 5 }