Last active
November 7, 2018 19:17
-
-
Save rbiswas4/ac8278bd7cd93aa85718973c284990e8 to your computer and use it in GitHub Desktop.
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
{ | |
"cells": [ | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "from astropy.cosmology import Planck15", | |
"execution_count": 1, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "import numpy as np", | |
"execution_count": 2, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "import matplotlib.pyplot as plt", | |
"execution_count": 3, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "import seaborn as sns\nsns.set_style('whitegrid')\nsns.set_context('talk')", | |
"execution_count": 4, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "#def m(t, t0, m0):\n# return m0 - 2.5 * n * np.log10(t - t0)", | |
"execution_count": 5, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "z = np.arange(1.0e-8, 1., 1.0e-3)\nmu = Planck15.distmod(z)\nM = np.arange(-20, -12, 0.1)", | |
"execution_count": 12, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "# Take m5 ZTF (median) as 20.5\nm5_ZTF = 20.5\nmu_ZTF = m5_ZTF - M ", | |
"execution_count": 13, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "m5_LSST = 24.\nm5_spec = 23.", | |
"execution_count": 14, | |
"outputs": [] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "z_ZTF = np.interp(mu_ZTF, mu, z)\nfig, ax = plt.subplots(2, sharex=True, squeeze=True)\nax[0].plot(M, z_ZTF, '.')\nax[1].plot(M, Planck15.comoving_volume(z_ZTF)* 1600/40000,'.')\nax[0].set_ylabel('z')\nax[1].set_ylabel('Comoving Vol')\nax[1].set_xlabel('M')", | |
"execution_count": 15, | |
"outputs": [ | |
{ | |
"output_type": "execute_result", | |
"execution_count": 15, | |
"data": { | |
"text/plain": "<matplotlib.text.Text at 0x1118ddeb8>" | |
}, | |
"metadata": {} | |
}, | |
{ | |
"output_type": "display_data", | |
"data": { | |
"text/plain": "<matplotlib.figure.Figure at 0x111671e48>", | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAEiCAYAAAAmtt8/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlYVOXiB/DvzACyKqDiAoRrYIYi2OKSK4ZUbqRBVppLmuk1ta5bZWZlmRqZmje9Vy+WpmlqIiqu6U3q8UkjKZVfRQmImooawzbb+f1BMzEwwAycc2aG+X6ep+fpnDPnnJdXnS/vct6jEARBABERkYyU9i4AERG5HoYPERHJjuFDRESyY/gQEZHsGD5ERCQ7hg8REcmO4UNERLJj+BARkewYPkREJDuGDxERyY7hQ0REsmP4EBGR7Bg+REQkO4YPERHJzs3eBbCnM2fOAABUKpWdS0JE5Dz0ej0AICYmpt7XYMunAfR6vekPgazDOrMd68x2rDPbyV1nLt3yMbZ4oqKi6nX+N9//hKtqLfpERSDIz1PMojVa2dnZAIDw8HA7l8R5sM5sxzqznS11lpmZ2eD7uXT4NMSus/mY/8UluCkVMBy6gqWjIpEQHWLvYhEROQV2u9XDH0VlWLg7Cxq9gBKtAWVaAxbuzsIfRWX2LhoRkVNg+NRDXmEp3JXmVeeuVCKvsNROJSIici4Mn3oIDfSC1mAw26c1GODlocSZS7fYAiIiqgPHfOohyM8TS0dFYv4X5yrGfAAk9AhGwkcZcFcqoTUYOAZERFQLhk89JUSHoI3iDq6qtQjv1B4JH2WgTGtAGSpaRAt3Z6Fv5xacBUdEZAG73Rog0NsN9wR5oVRj4BgQEZENGD4isDQGpNHrcatEw/EfIiILGD4iMI4Bebor4dfEDW5KwCAAs7dlot97x7HrbL69i0hE5FAYPiJJiA7BybkDkZwUBZVSCa1eQFG5js8AERFZwPARUZCfJwK8PeChqj7+cy7/DqdhExH9hbPdRGZp/KdUq8P0LWfhoeI0bCIigC0f0VUd/2nipoBCoUC5zsBuOCKiv9i95bNv3z6sW7cOeXl5CA4OxtSpUzFy5MgaP3/9+nWsWrUKp06dwu3bt9G+fXs899xziI+Pl7HUtUuIDkHfzi2QV1iKWyUazN6WCa1eZzpu7IYL8C5FaKAXnwUiIpdj1/DZv38/Xn75ZYwfPx59+/bFkSNHMG/ePHh6emLo0KHVPq/RaDB58mQUFRVh5syZCAoKQnp6OmbNmgW9Xo/HHnvMDj+FZUF+ngjy88QfRWXshiMiqsKu4ZOcnIz4+HgsWLAAAPDQQw/hzp07WLVqlcXwOXnyJC5evIgdO3agW7duAIA+ffqgoKAAGzZscKjwMTJ2wy3cnQV3pRIavR4GASjXGVCu42oIROSa7Dbmk5eXh9zcXDz88MNm++Pi4pCTk4O8vLxq5/j4+CAxMRGRkZFm+zt06IDc3FxJy9sQxmnY/514P9Y+FQNPN/PXdnM1BCJyNXZr+eTk5AAA2rdvb7Y/LCwMAPDbb78hNDTU7FivXr3Qq1cvs31arRYnTpxA586dJSxtw9XWDVd5RWyOARGRK7Bb+BQVFQEAfH19zfb7+PgAANRqtVXXWb58OX7//XesXbu2XuXQ6/Wm18faqqSkBABsPv8fD7bAqozrcFMqoDMIGNTeByPXnDJtv9i7JWI7Na1XmRxdfevMlbHObMc6s50tdabX66FSqer8XG3sFj6CINR6XKmsvUdQEAQsX74cKSkpmDRpEmJjY8UsnqRiOzVFdFtvXFVr4alS4MW0y9DoBWj0FXWyKuM6ott6I9Db7pMRiYgkYbdvNz8/PwBAcXGx2X5ji8d43BKNRoP58+cjLS0NkyZNwty5c+tdDpVKhfDw8Hqda/wNob7nA8CZS7fQxO0qNJWmYjdxU6HYswU8PD0aXTecGHXmalhntmOd2c6WOsvMzGzw/ewWPsaxntzcXLMf9tKlS2bHq1Kr1Zg6dSrOnj2LhQsXYvz48dIXVkJcEYGIXJHdZruFhYUhJCQEBw8eNNt/6NAhtGvXDm3btq12jl6vx7Rp0/DDDz8gOTnZ6YMH4IoIROSa7DqoMH36dCxYsADNmjXDgAEDcPToURw4cADJyckAgMLCQuTm5qJTp07w9fXFtm3bcPr0aSQmJqJ169ZmTT+FQoHu3bvb60dpEK6IQESuxq7hk5CQAI1Gg40bN2LHjh0IDQ3FsmXL8MgjjwAAvvrqKyxYsACbN2/GAw88gPT0dADA9u3bsX37drNrqVQqnD9/XvafQSxcEYGIXIndp1MlJSUhKSnJ4rGEhAQkJCSYtjdv3ixXsezGmhURFuw6h6Ze7ugW0oytICJySnYPH6qurm64cp2AmZ99D4MgsBVERE6Jr1RwUEF+nogJC0C3kGbVuuEAoESj52QEInJaDB8HV3k2nLd79SeK+ZZUInJG7HZzAsZuuHP5dzB9y1nT2A/AyQhE5JzY8nESQX6eiO3SCu8k8JkgInJ+bPk4GT4TRESNAcPHCfGZICJydux2c2LWLs1z/gonJBCRY2HLx8nV1Q0nGASMXJOBJm5sCRGR42D4NAK1dcOV6wUAAjT6iv0Ld2ehb+cWHAsiIrtit1sjUrUbzkOlgKe7+R8xnwsiIkfAlk8jU7kbzstDiYSPMsyOc0ICETkCtnwaIePSPPe0acYJCUTkkNjyaeQ4IYGIHBHDxwXYOiEhoo0fSjUGPqRKRJJht5sLsWZCgrEl9OzG0+j33nHsOptvp9ISUWPGlo+LqWtCQtWWEF9cR0RSYMvHBdU0IcFSS8j44jq2gohITDaFz549e5CXl1fj8czMTCxYsKDBhSL5JESH4OTcgfjvxPuxZ0Yfi5+p/OI6zowjIjHYFD7z58/H6NGjceLECYvHc3NzsWfPHlEKRvKx1BKy9OI6jgcRkVjq1e02bdo0rFmzRuyykAMwtoQ+HNsDTdyqdMHpK8aC+IwQETWUzeEzf/58JCYmYs2aNZg6dSr+/PNPKcpFdmTpxXWcGUdEYrJ5tpubmxtef/113HPPPViyZAkef/xxrF69GhEREVAoFFKUkezE1plxfEaIiKxV79luY8aMQUpKCkpLS5GUlIS9e/eiSZMmYpaNHIAtM+PYEiIiazXoOZ/o6Gh88cUXeOGFFzBv3jz07NlTrHKRA2JLiIjE0uCHTFu1aoXPPvsMr776Kvbu3cuut0bOuFQPACwdFYmFu7PgrlSiXKeHUqlAmfbv5XssrRvX1cdeJSciR2JT+GzevBkdO3astt/DwwPvvfceunfvjh9//FG0wpFjq09LKDm+Lcr0AgLalrElROTCbAqf+++/v9bjTz31VIMKQ87H1pbQzH35cFcpYTh0hStoE7kwru1GorGuJQTTytocEyJyXQwfElVDxoTmD41AZIg/g4jIBTB8SDK2jgktTj0PnyYq6A0Cu+SIGjmuak2SqvqckIdKAW93pcXnhACguJyLmBK5ArZ8SDYJ0SFoo7iDq2otwju1r9YSqszSNG1jK4rdckTOj+FDsgr0dkOgtxvC/2oJLdydBZVCgWKN3uxzVbvk5u78ASqlEh6qv8OI3XJEzovhQ3ZTeUwoK/823j14scbJCToDoDMYUK7jTDmixoDhQ3ZlnB0XExaAR7q1qXFyQlXsliNybgwfchg1TdPW6PUwCID2r+eEAOu65RhGRI6L4UMOqXKXXGigF77++UatzwxV7ZbjGBGRY2P4kMOq3BKq65mhqjhGROTYGD7kNGzplquKqykQORaGDzml2rrlrBkjqrqaAseHiOTF8CGnVVO3nDVjREDFagoAJysQ2QPDhxqN+o4RWTNZgWFEJC6GDzValsaILK2mUBXDiEh6DB9yCTWtpmDNZAWGEZH4GD7kMiytpmDNZIWqGEZEDcfwIZdk7WQFMcKI07qJqmP4EEHaMOK0bqLqGD5EFogZRgCndRNVxfAhsoJYYVSf8aLCEh2uqrUIaFvGcKJGg+FDVA/WhJEY07oTegRj55l8uCkVMBy6wpYSNRoMHyIR1BRGDZ3WvfV0HgBA89c5dU1m+KOojMFEToHhQyQBsaZ1V1XbZIaEHsHY9f1luCs5nkSOj+FDJDGxJy9UZZzMYGwllcH68SS2lMheGD5EMrM2jKqO+egFwaZwsmY8qa6WEsOJpMLwIbKz2sIoyM8Tw9qrcFWtRZ+oCJsmM1RV03hSTS2lquHE8SUSE8OHyMFUDiMACPR2Q6C3G4L8PGuczFA1LMQYT6oaTraOLzGcqDYMHyInU9NkhiA/T8wacreo40lVWTu+xHCiujB8iJxY1VaSreNJDWkpVdbQLr2q4QSAQdXIMXyIGrG6xpNqailVDoj6jC9VZUs4lWp1UCgU8HRT1bsVxVUhHB/Dh8iF2NJSqhxOYo8vVVU1nAABWr0OgO2tKGtWhagcVgBbWfbA8CEik5rCyZbxJSnCqbK6WlF1rQpRuWxitLLq2ibLGD5EZJX6tprkDqeq6gqrhrayatu2dSzLlYLN7uGzb98+rFu3Dnl5eQgODsbUqVMxcuTIGj9fXFyMFStW4NChQygpKUHPnj3xyiuvoF27dvIVmoiqESucamuNSB1UVVnbyqpp25axLHsHm9wUgiDI86dowf79+zFnzhyMHz8effv2xZEjR7Bt2zasWrUKQ4cOtXjOlClTkJWVhblz58LHxwdr1qzB7du3kZaWBj8/P5vun5mZCQCIioqqV/mzs7MBAOHh4fU63xWxzmznKnVm7ThMba2ohq4K4UzclBAt2JaOikRXn2IA1v09a+h3J2Dn8BkyZAjuvfdeJCcnm/bNmjUL2dnZOHDgQLXPf/fdd3jqqaewYcMG9OvXDwBQWFiIwYMHY9q0aZgyZYpN92f4yI91ZjvWWXV1/Rb/zfc/VVsVwtIXryO1suzJ012JTQl3IdDbTbbwsVu3W15eHnJzczFnzhyz/XFxcThw4ADy8vIQGhpqduzUqVPw8fFBnz59TPsCAwNx33334eTJkzaHDxE5p9q6+ICaV4Ww1OUHoF5jVXVtO1N4uSuVuKrWItBbvkiwW/jk5OQAANq3b2+2PywsDADw22+/VQufnJwchIWFQaVSme2/6667LLaUiIiAusOqvmNVdW1bO5Zl72DTGgxo7esuyrWsZbfwKSoqAgD4+vqa7ffx8QEAqNXqaueo1epqnzeeY+nz1tDr9aZuDVuVlJQAQL3Pd0WsM9uxzmwndp35ArhVANyycburD7Ap4S5cVWtNX+7G/w/0dsOw9nfVe/tsQQlWZVyHm1KBcp0BCgXgoVJCZxAQ29EXR35Vw02psGr7Hw+2gCc0KCnRWFVner2+WiPAVnYLn7qGmpRKpU3nWPo8EZG9GbsAK2/Xdsza7dhOTRHd1rvGYHumR3Obto2BLRe7hY9xZlpxcbHZfmMLxtLMNV9fX+Tn51fbX1xcbLFFVBe9vmLJkNLSUpvPBQCDwdCg810R68x2rDPbuUqdeSmA9n4AoAXw9/+XlmrNjlmzbWudGb8/68tuzQXjWE9ubq7Z/kuXLpkdr3pOXl5etRbQpUuXLH5eaiqVqsFNT1fDOrMd68x2rDPbyV1ndmv5hIWFISQkBAcPHsSQIUNM+w8dOoR27dqhbdu21c7p27cv/vWvfyEjI8M0462wsBDfffcdpk6danMZYmJi6v8DEBFRvakWL1682F439/Pzw7p163Dr1i0oFAps3LgRe/bsweuvv47OnTujsLAQ2dnZ8PX1hYeHB4KDg3H69Gls3boV/v7+KCgowMKFCyEIApYuXQpPT+deboKIyFXY9SFTANi2bRs2btyIK1euIDQ0FFOmTDEtr7Nr1y4sWLAAmzdvxgMPPAAAuHPnDt59910cOXIEBoMBMTExmD9/Pjp06GDPH4OIiGxg9/AhIiLXw/nJFly4cAFdu3bF1atXbT734sWLmDBhAnr27Il+/fph8eLF9X4GiYiosbL7qtaO5tdff8XUqVOh0+lsPvf69esYP348QkNDsWLFCty6dQvvvfceCgoKsH79eglKS0TknBg+f9HpdNi+fTtWrlwJd/f6LTNx7Ngx3L59G7t37zbN1tPr9XjllVdw+fJlBAcHi1lkIiKnxW63v5w5cwYrVqzAxIkT8fLLL1v8zOHDh5GQkIDIyEj07dsXy5Ytg0ajMR0vLy8H8PcSQQDg7+8PALh9+7aEpScici4Mn7907NgRR44cwYwZMyw+aJWamooZM2agc+fOWLt2LZ5//nls374dL730kukz8fHxaNmyJd5++23cvHkTOTk5WLt2Le6++25ERETI+eMQETk0drv9pUWLFjUeEwQBK1aswMCBA7Fs2TLT/tatW2P69Ok4c+YMYmJi0LJlSyxevBhz5szBl19+CQBo27YtPv30Uz5tTURUCVs+VsjJycHVq1cxaNAg6HQ6038PPfQQ3N3dkZGRAeDv1tGQIUOwceNGfPjhh2jatCkmTJiAGzdu2PmnICJyHGz5WME4XvPaa6/htddeq3b8jz/+AACsWbMG9913H1auXGk6dt9992HIkCH4z3/+g3nz5slTYCIiB8fwsYJxhe0FCxZYXA8uICAAAHD58mUMHTrU7FhgYCA6duyIn3/+WfqCEhE5CXa7WaFjx44IDAzE5cuXERkZafovICAAK1aswK+//gqgYtXtM2fOmJ17584d/Pbbb5xmTURUiegtn7i4uHqdl56eLnJJxKNSqTBr1iy88cYbUCqV6NevH27fvo3Vq1ejqKgI99xzDwDgxRdfxIwZM/DSSy8hISEBarUaH3/8MQwGAyZOnGjnn4KIyHGIHj6BgYFQKBRiX9buEhMT4evri3//+9/YunUrfH19cd9992HOnDlo2bIlACA2Nhbr1q3DunXr8Pzzz6NZs2bo0aMHkpOTERYWZuefgIjIcXBhUSIikp1sEw6KioqQkZGBgoICuLu7o02bNnjwwQfNVgOQm3F8hs/gEBFZz/gK7Ya8kFOW8Pn888/x7rvvorS01OwV2F5eXpg3bx6SkpLkKIbojH8ADC/rsc5sxzqzHevMdnLXmeThc+zYMSxatAgRERGYOHEiOnbsCIPBgJycHKSkpOCNN95AmzZt0L9/f6mLUo2xkqOioup1fnZ2NgAgPDxctDI1dqwz27HObMc6s50tdZaZmdng+0kePuvXr0eXLl2wfft2eHh4mPZHRkYiPj4eSUlJ2LBhg13Cp6EKS3S4qtYioG0Zgvz4Cm8iImtJ/pzPxYsXMWLECLPgMfLw8MDw4cNx4cIFqYshul1n8zF+5yW8cugK+r13HLvO5tu7SERETkPy8HF3dze9asCS0tJSp5ua/UdRGRbuzoJGL6BEa0CZ1oCFu7PwR1GZvYtGROQUJA+fmJgYbN261eLCmtevX8dnn33WoBkT9pBXWAp3pXnVuSuVyCsstVOJiIici+RjPrNmzUJiYiLi4+MxatQotGvXDkDFStFffvklNBoNZs6cKXUxRBUa6AWtwWC2T2swIDTQy04lIiJyLpKHT0REBFJSUrBkyRJs3rzZ7FiXLl3w6quvomvXrlIXQ1RBfp5YOioS8784BzelAgYAS0dFctIBEZGVZHnOJyoqCrt27cK1a9dw+fJlCIKAkJAQtGrVSo7bSyIhOgRtFHdwVa1Fn6gIBg8RkQ1ED5+EhASMGDECjz76aLW3g7Zq1cqpA6eqQG83BHq7MXiIiGwk+oSDa9eu4Z133sGAAQMwefJk7N27F6WlHIgnIqK/id7y+d///odvv/0W+/fvx+HDh/H111/Dy8sLsbGxGD58OPr06QOlkq8RIiJyZaKHj1KpRO/evdG7d28sXrwYp06dQlpaGo4dO4Z9+/ahefPmeOSRRzB8+HDce++9Yt/erv4oKkNeYSlCA73YFUdEVAtJJxy4ubmhf//+6N+/PzQaDU6cOIG0tDTs2LEDn3zyCdq1a4dhw4Zh+PDhCAkJkbIoktt1Nh8Ld2fBXamE1mDA0lGRSIh27p+JiEgqsr1SwcPDA0OGDMGQIUNQWlqKkydP4vjx49i0aRPWrFmD8+fPy1UU0RlXPCjTGlCGiud/Fu7OQt/OLdgCIiKywC6DL/n5+fjll1/w66+/Qq1Ww9vb2x7FEA1XPCAiso1sLZ+LFy8iPT0d6enp+O2336BSqdCvXz8kJydj0KBBchVDElzxgIjINpKGz/nz55Geno6DBw8iNzcXABAdHY1x48YhPj4ezZo1k/L2sjGueFB1zIddbkRElokePllZWaYWTn5+PgRBQKdOnTBr1iwMGzYMbdu2FfuWDiEhOgR9O7fgbDciIiuIHj5jxowBAAQFBeHZZ5/F8OHD0aVLF7Fv45CC/DwZOkREVpBkeZ1hw4bhwQcfdLr39IiJz/wQEdVM9PBZunSp2Jd0Onzmh4iodlznRmSVn/kpKtfxLadERBYwfETGZ36IiOrG8BEZn/khIqobw0dkxmd+PN2V8GviBk93JZ/5ISKqQrYVDlwJn/khIqqd5OEzaNCgWqdcKxQKeHh4oHnz5ujWrRsmTJhQ7Q2ozojP/BAR1UzybrdevXpBrVbj8uXLaNKkCbp06YKoqCj4+/ujoKAAN27cQEBAAG7fvo2NGzdi5MiRKCgokLpYsvujqAxnLt3irDciIsjQ8rnnnnuQmpqKjz76qNoCopmZmZg4cSJGjhyJMWPGIDs7G5MmTcKqVauwbNkyqYsmGz73Q0RkTvKWz6ZNmzBu3DiLK1dHRUXhmWeewfr16wEA4eHhePLJJ3Hq1CmpiyUbPvdDRFSd5OFz8+ZNtGrVqsbjzZs3x7Vr10zbQUFBUKvVUhdLNnzuh4ioOsnDp1OnTti9ezc0Gk21YxqNBnv27EGHDh1M+3766adGtfI1n/shIqpO8vCZMWMGsrOzMWLECKSkpOCrr75CRkYGtmzZgqSkJFy4cAHTpk0DACxevBg7d+7EiBEjpC6WbPjcDxFRdZJPOOjfvz/WrFmDpUuX4p133jFNuxYEAW3atMEHH3yAuLg4FBYWYufOnRg2bBgmTpwodbFkxed+iIjMyfKQ6cCBAzFw4EBcvHgRubm50Ol0CAkJQWRkpCmM/P398f3338Pd3V2OIsmu6nM/fOUCEbkyWVc4iIiIQEREhMVjSqUSSqVrrPbDqddE5OpkCZ+TJ08iNTUVN27cgF6vr3ZcoVAgJSVFjqLYXeWp12WomIiwcHcW+nZuwRYQEbkMycNny5YteOuttwBUTKv28PAQ9foXLlzA6NGjcfToUbRu3VrUa0vBOPXaGDzA31OvGT5E5CokD5/NmzcjIiICGzZsEH3Ntl9//RVTp06FTqcT9bpS4tRrIiIZplpfuXIFiYmJogaPTqfDli1bMGbMGJSXl4t2XTlw6jURkQwtn7vuugs3btwQ9ZpnzpzBihUrMGnSJLRq1QqvvvqqqNeXGqdeE5Grkzx8pkyZgrfffhtxcXHo3LmzKNfs2LEjjhw5gubNm2PXrl0NupZer0d2dna9zi0pKQGAep/vC+BWAXALQGGJDlfVWrT2dUegd+N9zVJD68wVsc5sxzqznS11ptfroVKpGnQ/yb/lzpw5Ax8fH4wYMQLt27dHYGBgtff72DrbrTG876eyI7/8iVUZ1+GmVEBnEPBi75aI7dTU3sUiIpKM5OHzv//9DwDQunVrlJaW4vLly1Lf0iYqlQrh4eH1Otf4G0J9zwcqpl6v/vQ4NHoBGr0AAFj97Q2M6dc4x4HEqDNXwzqzHevMdrbUWWZmZoPvJ3n4HDt2TOpbODVOvSYiV+QaSwo4ME69JiJXJHrLZ9y4cZg2bRp69epl2q6LK61wUJVx6nXV5XYA4MylW5wNR0SNkujhk5+fb5o1Ydym2lWdev31zzfQ773jXPuNiBot0cOn6hgPx3ysY1z1mmu/EZErkHzM56OPPkJeXp5k109ISEB2drZTrOtmDb52m4hcgeTh8+GHH+Lhhx9GUlIStmzZgsLCQqlv6dQ4AYGIXIHk4XP06FHMnj0bGo0Gb775Jvr164fnnnsOqampKC3lb/NV1bT2G1AxAeGPojI7l5CIqOEkf84nODgYU6ZMwZQpU/D7778jLS0NBw8exD//+U94eXlh0KBBGD58OPr37y91UZwGJyAQUWMn63M+7dq1w/Tp05GamorU1FQ88MADSEtLw/PPPy9nMZxCkJ8nYsICAMA0AaGoXIcyrQELd2exBURETk3WFSz1ej2++eYbHDhwAMePH0dhYSGCg4Px2GOPyVkMp8IVEIioMZI8fAwGAzIyMnDw4EEcOXIEt2/fhr+/P+Lj4zFs2DBER0dLXQSnxgkIRNQYSR4+vXr1wp9//glPT0/T+E6fPn3g5tZ4XxsgJksrIMwfGmGaes3WDxE5I8kToFu3bhg2bBhiY2Ph7e0t9e0apcoTELLyb+Pdgxc5+YCInJrk4bNhwwbT/9+4cQMFBQVwd3dHq1atEBgYKPXtGw1jC+epf3/L1Q+IyOnJ0vf1448/YsmSJcjKyjLb3717d7zyyiuIjIyUoxhOj5MPiKixkDx8srOz8cwzzwAAnnjiCXTs2BEGgwE5OTlITU3FuHHj8Pnnn4v2iu3GrKbJB14eSq6ATURORfLw+eCDD+Dj44Pt27cjODjY7NgLL7yA0aNHY82aNVi1apXURXF6liYfJPQIRsJHGRwDIiKnInn4fPfdd5gwYUK14AEqXq395JNP4pNPPpG6GI1G5ckHXh5KJHyUwTEgInI6kq9woNFo4OPjU+NxX19flJXxaX1bGFc/KNUYLK6AfS7/DteBIyKHJnn4dOnSBfv27YNOp6t2TKvVIjU1FXfffbfUxWiULI0BlWp1mL7lLJ7deBr93juOXWf5Mj8icjySh8/kyZORlZWFp59+Gunp6cjOzkZ2djYOHDiAp59+Gj/99BMmTpwodTEapaorYDdxU0ChUKBcx3XgiMixST7mExsbi9deew0rVqzArFmzTPsFQUCTJk0wb948DB06VOpiNFqVx4BulWgwe1smtPq/W5nGbrgA71LOhiMihyHLcz5PPfUUHn30UWRkZODy5csQBAEhISHo3bs3/P395ShCo1b5FdyrO47MAAATMElEQVQ1dcN5qDgbjogch2wLrPn7++ORRx6R63YuqepUbI1eD4MAlOsMKNdxNhwROQ5ZwmfPnj04deoUrl+/DkOV38wBQKFQICUlRY6iNHrshiMiZyB5+CQnJ+Pjjz+Gu7s7mjdvDqVS1vfXuSR2wxGRo5M8fHbv3o2+ffti9erV8PLiO2jkxG44InJUkoePWq1GXFwcg8dO2A1HRI5I8vB56KGH8O2332LMmDFS34pqwG44InI0kofPa6+9hgkTJuCll15CbGwsmjdvDoVCUe1z9913n9RFcXnWdsNFtPFDqcbAlhARSUby8CkoKEBRURHS0tKwf//+ascFQYBCocCFCxekLgqh7m44wSBg5JoMNHFjS4iIpCN5+CxZsgR//vknJk2ahHbt2sHNTbZHi6gGtXXDlesFAAI0+or9C3adQ1Mvd3QLacZWEBGJRvIk+PnnnzFjxgw899xzUt+KbFS1G65cp4dSqUCZ9u9AKtcJmPnZ9zAIAltBRCQaycOndevWfLbHgVl6P1BVJRo9AI4HEZF4ZFnVOiUlBb/88ovUt6J6Mr4f6J42zUyrZHu7q6p9zjgeVPl1DX8UlfHdQURkM8lbPhcvXoRCocDw4cMRGhqKFi1aQKUy/2Lj8jqOw9gSOpd/B9O3nDXNggOqjwfN3fkDVEql2TRtYyuKLSMiqo3k4XP8+HGoVCq0bt0aWq0WV65ckfqW1EBBfp6I7eKJdxJqHw/SGQCd4e9p2pbCiGNERGSJ5OFz7NgxqW9BErFmPKiyqmHEMSIiqols8571ej1+/PFHXL58GR4eHmjTpg26du0q1+2pnozTsgFYfEBVqxdqPNfSM0NtFDpcVWsR0LaMYUTkwmQJn+PHj+ONN97AtWvXIAgVX1YKhQJBQUF4/fXXMWjQIDmKQQ1UuSUUGuiFr3++UWsYWRojUioUcFMqYDh0hWNERC5M8vD57rvv8I9//APNmzfH7Nmz0bFjRwiCgJycHGzduhUzZ87E5s2bER0dLXVRSASVW0K1hVFNY0QVYVQRUFXHiOYPjUBkiL8piP4oKmMwETVSkofP6tWrERwcjJ07d8LPz8/s2NixY/H4449j3bp12LBhg9RFIQnUFEb1GSNanHoePk1U0BsEJPQIxq7vL8NdyZl0RI2R5OFz7tw5TJ8+vVrwAICvry9Gjx7N4GlEGjJGBADF5RUPtG49nQcAKIP107rZUiJyHnZfaE2hUECr1dp0zr59+7Bu3Trk5eUhODgYU6dOxciRIyUqIdWXpW65+V+cg5tSAb0gWBVGRnVN67ampVQ5nAAwqIjsSPLw6d69O3bu3ImxY8fC29vb7JharcaOHTsQGRlp9fX279+Pl19+GePHj0ffvn1x5MgRzJs3D56enhg6dKjYxacGqtot10ZxB1fVWvSJijCNEakUChT/tYSPtaqGUV0tpcrhVKrVQaFQwNNNZXUrytZtIqqd5OEzY8YMjBs3Do899hiefvpptGvXDgBMEw6uXbuGN954w+rrJScnIz4+HgsWLABQ8bK6O3fuYNWqVQwfJxDo7YZAbzcE+XmatYyy8m/j3YMXTS2XymFhbZddZXWFEyCYXiVRVyvK1u2qYQagQUFWWGI+Pb22FhxDkpyF5OHTs2dPrF69GkuWLMF7771nepGcIAho2bIl3n//fTz44INWXSsvLw+5ubmYM2eO2f64uDgcOHAAeXl5CA0NFf1nIOkYW0YxYQF4pFsbsy/GWUPutnpad0PUFVS2blcOs6qtrPoE284z+abp6bW14BoaknWFFdCwoLNlu6H3qi2wnfnnkvLaclMIxgdvJKbX6/HTTz8hPz8fABAcHIyuXbva9H6fEydOYMqUKfjyyy8RERFh2n/+/HmMGjUKGzZsQL9+/ay+XmZmJgAgKirK6nMqy87OBgCEh4fX63xX1JA6q/yPpXIYidFSckVuSljd4hMjRK3dFjWwgUb1c0l17aWjItHVpxiAdf82G/rdCcgYPpXdvHkT/v7+1RYYrcu+ffvw0ksv4ejRowgJ+XvNsEuXLuHhhx9GcnIyHnnkEauvl5mZCb1eD19fX5vKYVRSUgIA1cayqGZi1pnxt9vWvu4I9HYz2z5bUIJVGdfhplRAZxAQ29EXR35Vw02pQLnOAIUCFV+6egMEGJ9BInJNHioF/vVoEAK8VFb921Sr1VCpVA0KH8m63T799FNs374du3fvrta6Wbp0Kb755htMmTIFzz77rNXXrCsn+d4g12IcP7K0HdupKaLbepuF0zM9mpu2AVgVVLZuM8zIGbkpFbhWrEeAl20NggbdU+wLCoKAefPmYe/evWjWrBkKCgpw1113mX0mJCQESqUSy5Ytw7lz5/D+++9bdW3js0LFxcVm+9VqtdlxW6hUqnp3m7HbzXaOWGe9egBj+onX9165S1DKLqSGXJtdk1SZAUC7Fr7w9nazqdutIUQPnx07dmDv3r0YO3Ys5s2bhyZNmlT7zOzZs/HCCy/g9ddfx5dffol+/fpZ9ZxO+/btAQC5ublmFXTp0iWz40S2qjwlvKHbVZ9vAswHgitPpLBme1h7lWl6etXjDbl2beNmzj42wjEf28d8Ar3Nf6mXmuhjPmPGjIGnpyc++eSTOj9rMBjw+OOPo0mTJti2bZtV1x88eDCioqKwcuVK075Zs2bhwoULSE9Pt6msZ86cAQCbx56M9Hp9g853Rawz20lZZwZBgCAACgWgVChq3QZg9Wcbut3Qe+n0eggCoFIpG9XPJeW1bfl7ZvxsTExMvf/uid7y+eWXX/Diiy9a9VmlUom4uDh8/PHHVl9/+vTpWLBgAZo1a4YBAwbg6NGjOHDgAJKTk+tb5HrjF6jtWGe2k7LOlAoFoLB+25bPNnS7Iee6VamzxvJzSXltuf9tih4+KpUKHh4eVn8+ICDApokCCQkJ0Gg02LhxI3bs2IHQ0FAsW7bMplluRg1JbSIiqj/RwycsLAw//vij1Z/PyspC27ZtbbpHUlISkpKSbC0aERE5CNHnJj/66KNITU3Fzz//XOdnf/75Z6Smptr0YCgRETk/0cMnMTERbdu2xTPPPIO9e/eaBqYqMxgMSE1NxYQJE+Dj44Px48eLXQwiInJgkqxwkJOTgxdeeAGXLl2Ct7c3unbtipYtW8JgMODmzZv46aefUFJSgjZt2mDt2rXo0qWL2EUgIiIHJtnyOhqNBlu2bEFaWhouXrwIna5iBWF3d3dERUXh4YcfRmJiok2TE4iIqHGQbW23wsJCqFQqNGvWTI7bERGRA7PLwqJEROTauBInERHJjuFDRESyY/gQEZHsGD5ERCQ7hk89XL9+Ha+++ioGDhyIHj16ICEhAQcOHDD7jE6nwwcffID+/fuje/fuGDt2LM6dO2enEjuWZcuWWXyJYHFxMZYuXYqBAwciOjoazzzzDOvsLzXVGQBs27YN8fHxiIyMRFxcHDZv3ixv4RxUbXVmpFarMWjQILzyyivyFMrB1VRnarUay5YtQ2xsLKKiojBs2DBs3bq1zhd81obhYyONRoPJkycjIyMDM2fOxJo1a3Dvvfdi1qxZ2Ldvn+lzb7/9Nv773//iueeeQ3JyMlQqFZ599lnk5eXZsfT298knn2Djxo0Wjy1atAg7duzA5MmT8eGHH8LDwwPjx49nndVSZ5s2bcIbb7yBuLg4rF+/HnFxcXj77bfx2WefyVxKx1JbnVX2zjvv4PLlyzKUyPHVVmezZ8/Grl27MGHCBKxbtw4DBgzAm2++ifXr19f/hgLZ5PDhw8Ldd98t/PDDD2b7J02aJAwfPlwQBEHIy8sTunTpImzdutV0vLy8XBgwYICwaNEiWcvrKK5evSrMmTNHiIiIEGJiYoTx48ebHS8tLRW6dOkirF692rRPrVYL3bp1M9vnSuqqM7VaLURFRQnJyclm++fMmSNMnz5dxpI6jrrqrLKvvvpK6NGjhxATEyMsXLhQvkI6mLrq7Pz588Ldd98t7N+/32z/okWLhJiYmHrfly0fG/n4+CAxMRGRkZFm+zt06IDc3FwAwLfffgu9Xo+4uDjTcQ8PDwwYMAAnT56UtbyOIjk5GefPn8emTZssLqek1WphMBjg6+tr2uft7Y0mTZrg9u3bchbVYdRVZ19//TVKSkowduxYs/0rV67EmjVr5CqmQ6mrzozu3LmDV199Ff/85z/RtGlTGUvoeOqqM0EQkJiYiF69epnt79ChA4qKinDr1q163ZfhY6NevXphyZIlUCj+fguTVqvFiRMn0LlzZwAVa9s1a9YMgYGBZueGhYWhoKAAZWVlspbZEUyePBlpaWl48MEHLR738/PDqFGjkJKSgnPnzuHOnTtYuXIliouL6/WupsagrjrLzs6Gv78/rly5gqSkJNx7773o37+/S4/51FVnRm+++SY6duzIV7Og7jq75557sGTJEvj7+5vtP3LkCFq2bFltv7VEf5+PM9PpdNixY0eNx4OCgjB48OBq+5cvX47ff/8da9euBVAxOFf5N3gjHx8fABUD656eniKV2r6srbNOnTrVea3Zs2djypQpGDNmDABAoVDgrbfeQnR0tGjldQRi1VlhYSG0Wi2mTZuGyZMn48UXX8Thw4fx9ttvw9fXFwkJCWIX3W7E/Ht2+PBhHD16FPv27TP7JbKxEbPOqkpJScHp06excOHCetchw6eS8vJyLF68uMbj999/v1n4CIKA5cuXIyUlBZMmTUJsbKxpf20a0194W+usJjdv3sQTTzwBDw8PrFy5Es2bN0d6ejoWLVoEb2/vRtX6EavOtFotiouLMWfOHDz99NMAKlrmBQUFWL16daMKH7HqrLCwEIsWLcLcuXMRHBwsYgkdj1h1VtWnn36Kd955B/Hx8Rg3bly9y8fwqcTHxwfZ2dlWfVaj0WD+/PlIS0vDpEmTMHfuXNMxX19fFBcXVztHrVabjjcWttRZbXbs2IGrV6/i8OHDCA0NBVDxRVpUVIQ333wT8fHxjSa0xaozY0u6f//+ZvsfeughHD9+HEVFRfDz82vwfRyBWHW2ePFidOrUCaNHjzattA9U/MKo0+ng5tZ4vhLFqjMjg8GA5cuXY+PGjXjsscewbNmyBv2bbDw1LSO1Wo2pU6fi7NmzWLhwYbWX4XXo0AG3b9/GnTt3zFbxvnTpEkJCQvgaCQsKCgrQsmVLU/AY9ezZE/v27cPNmzfRokULO5XOMYWFhQGo+EWoMq1WC6BxtbDFkp6eDgC49957zfZ/8cUX+OKLL3D06FGEhITYo2gOTavV4qWXXkJ6ejomTpyIuXPnNvjvF8PHRnq9HtOmTcMPP/yA5ORkDB06tNpnevfuDaDiL/oTTzwBoOIL4quvvsJDDz0ka3mdRfv27bFjxw5cunTJ9KUKAJmZmfD19eWrOCww/l1KS0vDzJkzTfuPHz+O8PDwRtXCFsvOnTur7Zs2bRq6deuGadOmISgoyA6lcnwLFy7EoUOHsGDBgjof3LUWw8dG27Ztw+nTp5GYmIjWrVsjMzPTdEyhUKB79+4IDg7GqFGj8NZbb6GkpARhYWHYtGkT/vzzT0yePNmOpXdco0ePxieffIIpU6ZgxowZaN68OY4dO4Y9e/bg5Zdfhru7u72L6HDuuusuPPnkk/j444/h5uaGqKgopKWl4dtvvzVNfiFzVR+RACoegwgICLB4jICvvvoKe/fuxaBBgxAVFWX2nQdUzIarT28Ow8dGxmb79u3bsX37drNjKpUK58+fBwAsWbIETZs2xfr161FSUoKuXbti06ZNZr/V09/8/Pzw2WefYfny5Xjrrbeg0WjQoUMHvP/++3j00UftXTyHtWjRIrRp0waff/451q1bh/bt22P16tWmyS9EDWX8zjt27BiOHTtW7fiJEyfQunVrm6/Ll8kREZHs+JApERHJjuFDRESyY/gQEZHsGD5ERCQ7hg8REcmO4UNERLJj+BA5gPnz5yM8PBxdunRBYWFhjZ8bMWIEwsPDMX/+fBlLRyQ+hg+RAzEYDDh+/LjFY3l5ebh48aLMJSKSBsOHyIGEhITg6NGjFo8dOXKk2gsKiZwVw4fIgQwePBgZGRkW33Z7+PBhDBo0yA6lIhIfw4fIgcTGxqK0tBQZGRlm+2/evInvv/8eDz/8sJ1KRiQuhg+RA4mJiUFAQEC1rrejR4/Cy8sLvXr1slPJiMTF8CFyICqVCgMHDsTx48dhMBhM+w8fPowBAwbwRYTUaDB8iBzM4MGDcfPmTdN7U9RqNb755hu+JoEaFYYPkYPp06cPPD09Te9OOXHiBJRKJfr372/nkhGJh+FD5GC8vLzQp08f07jP4cOH0bt3b/j4+Ni5ZETiYfgQOaDY2Fjk5OTg//7v/3Dy5EkMGTLE3kUiEhXDh8gBDRw4ECqVCsuWLUNZWRmf76FGh+FD5IACAgIQExODr7/+Gj179kRAQIC9i0QkKoYPkYMaPHgwAPDBUmqUFIIgCPYuBBERuRa2fIiISHYMHyIikh3Dh4iIZMfwISIi2TF8iIhIdgwfIiKSHcOHiIhkx/AhIiLZMXyIiEh2DB8iIpLd/wPdRBj4TXsAygAAAABJRU5ErkJggg==\n" | |
}, | |
"metadata": {} | |
} | |
] | |
}, | |
{ | |
"metadata": {}, | |
"cell_type": "markdown", | |
"source": "As LSST and the spec survey have higher m5 (24, 23) values and go deeper, their volume for a similar sky area to the ZTF high cadence field (~ 20.5, 1600 sq Deg) is deeper. How big must the sky area in these surveys be to have the same volume as the ZTF high cadence fields for objects with different intrinsic brightnesses?" | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "mu_spec = m5_spec - M\nzspec = np.interp(mu_spec, mu, z)\nmu_LSST = m5_LSST - M\nzLSST = np.interp(mu_LSST, mu, z)\nfig, ax = plt.subplots()\nax.plot(M, Planck15.comoving_volume(z_ZTF) * 1600 / Planck15.comoving_volume(zspec), '.', label='Spec Depth')\nax.plot(M, Planck15.comoving_volume(z_ZTF) * 1600 / Planck15.comoving_volume(zLSST), 'o', label='LSST Depth')\nax.set_ylim(ymax=100)\nax.set_xlabel('M')\nax.set_ylabel('Equivalent Area at improved m5 (sq Deg)')\nplt.legend(loc='best')", | |
"execution_count": 16, | |
"outputs": [ | |
{ | |
"output_type": "execute_result", | |
"execution_count": 16, | |
"data": { | |
"text/plain": "<matplotlib.legend.Legend at 0x111a88940>" | |
}, | |
"metadata": {} | |
}, | |
{ | |
"output_type": "display_data", | |
"data": { | |
"text/plain": "<matplotlib.figure.Figure at 0x111a31da0>", | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAaQAAAFYCAYAAAAcHIhXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlYlFX7wPHvzLAKmID7rryl4m5KiuIubmhJuSS5QKaYmqJlZuavRcslNXczt1Tc00pzN60w802x0lwqMxR3UQFlneX3x7wzMawzMAMD3J/r4sJ5nmdmbp5obs459zlHodPpdAghhBBFTFnUAQghhBAgCUkIIYSdkIQkhBDCLkhCEkIIYRckIQkhhLALkpCEEELYBUlIQggh7IIkJCGEEHbBwdwL7969y/Hjxzl16hSxsbE8ePAAlUqFt7c3VapUwd/fH39/f8qWLWvLeIUQQpRQirxWavjll1/4/PPPOXToEGq1GicnJ6pWrYq7uztarZaHDx8SFxdHamoqDg4OBAYGEhYWRqNGjQrrZxBCCFEC5JiQ4uLi+PDDD9m3bx/PPPMM3bt3x8/Pj9q1a6NUmvb0qdVqLl68yKlTp9i7dy+//fYbgYGBTJ8+nfLlyxfKDyKEEKJ4yzEhtW7dmsDAQF555RVq1Khh0Yv+/fffrFy5kqNHj3Ly5EmrBCqEEKJkyzEhXb16lZo1axboxf/55x9q165doNcQQghROuQ5hiSEEEIUBrOr7KKjo3M9r1AocHJywtvbm8qVKxc4MCGEEKWL2S2k+vXro1AozHrR8uXLM2nSJJ577rkCBSeEEKL0MDshHTp0iGnTpuHi4sKLL75I3bp1cXZ25sqVK3zxxRdcu3aNsWPHotPp2L9/P+fPn2fZsmV06tTJ1j+DEEKIEsDshDRlyhTOnj3Ltm3bcHNzMzmXmprKoEGDaNSoER988AFarZZhw4ah0WjYtGmTTQIXQghRspi9dNDhw4fp379/lmQE4OzsTL9+/di/f7/+RZVKevTowR9//GG9SIUQQpRoZickBwcHEhIScjz/4MEDMje2VCpV/iMTQghRqpidkNq2bcu6des4depUlnPnzp1j/fr1tG7dGoC0tDR2797NU089Zb1IhRBClGhmjyHdvn2bF198kZs3b9KgQQNq166Nk5MT//zzD7/99hsVKlQgMjKSqlWr4u/vT0JCAp999hlt27a19c8ghBCiBLBoYuyjR49Ys2YNhw8fJiYmBrVaTfXq1enWrRsjRoygXLly3L9/n48++ohnn32Wdu3a2TJ2IYQQJYis1CCEEMIumL1Sg8Fff/3FsWPHuHXrFoMHD6ZMmTJcvnxZuuaEEEIUiEUJ6cMPP2TDhg3odDoUCgVdunQhMTGR1157jS5durBgwQKcnJxsFasQQogSzOwqu40bN7J+/XrCwsLYvHmzscTbz8+Pl156iSNHjvDZZ5/ZLFAhhBAlm9kJafPmzQQGBvLGG2+YbClRrlw5pk2bRlBQELt377ZFjEIIIUoBsxPS1atX8ff3z/G8n58fN2/etEpQQgghSh+zE1K5cuW4fft2jucvX75M2bJlrRKUEEKI0sfshNS1a1c2bdrE5cuXjccM21FERUWxZcsWWdlbCCFEvpk9Dyk+Pp4XX3yR69evU79+fX777TdatWpFUlISv//+O1WqVGHbtm2UL1/e1jELIYQogSyaGPv48WNWrlzJoUOHuHr1KjqdjqpVq9K5c2fCw8Px9PS0ZaxCCCFKMFmpQQghhF0we2JsfHw8ly5d4uHDhygUCjw9Palfvz7u7u62jE8IIUQpkWdCOnz4MCtXruTcuXNZ9jtSKBS0atWKESNGEBAQYLMghRBClHy5dtnNmDGDyMhI3N3d6dChA0899RRly5ZFrVbz8OFDfv/9d06cOEFqaiphYWG88cYbhRm7EEKIEiTHFtKePXvYuHEjvXv35r333suxa+7Ro0fMmDGDNWvW0KxZM7p162azYIUQQpRcObaQhgwZQmpqKlu3bjXON8pN//79KVu2LKtXr7Z6kEIIIUq+HCfG/vHHHwQGBpqVjAC6devG+fPnrRaYEEKI0iXHhPTo0SO8vb3NfqEKFSqQkJBglaCEEEKUPjkmJI1Gg4OD+dslqVQqtFqtVYISQghR+pi9lp0QQghhS7k2gY4cOcL169fNeqELFy5YJSAhhBClU45VdvXr17f8xRQKSUxCCCHyJccW0tq1awszDiGEEKWcLK4qhBDCLkhRgxBCCLsgCUkIIYRdkIQkhBDCLkhCEkIIYRckIQkhhLALZq8NdPv27Xy9QaVKlfL1PCGEEKWL2WXf9evXN3vl74xkoqwQQghzmN1CmjlzJsuXLycuLo6+ffvi4+ODs7MzMTExfPnllyQlJfHiiy/i7Oxsy3iFEEKUUGYnpCtXrqDRaNi/f3+WbrjRo0czaNAg0tPTefPNN60epBBCiJLP7KKGL774gpCQkGzHhDw8PBgwYAC7d++2anBCCCFKD7MTklqtJi0tLcfz9+/fL1AgFy5coGHDhty6dcvkeFRUFM8//zxNmzalc+fOrFmzJstzz549y5AhQ2jevDnt2rVj/vz5pKenFygeIYQQhcvshPTMM8+wdu1afvnllyznTpw4wfr16+ncuXO+grh8+TKjRo1CrVabHI+OjiY8PJy6deuyePFi+vTpw5w5c1i9erXxmpiYGIYPH46zszOffPIJYWFhrF27lo8++ihfsQghhCgaZlfZXb16lYEDB/Lw4UN8fX2pWbMmOp2Of/75h0uXLlGnTh0iIyPx9PQ0+83VajVbt25l3rx5ODo68vDhQ7777jsqV64MwPDhw0lKSmLbtm3G58ydO5dt27Zx/PhxnJycePvttzl+/DgHDx7EyckJgE2bNjFjxgyOHj0qZedCCFFMmN1CqlmzJrt37yYkJITExEQOHTrE4cOHUavVjB49mu3bt1uUjABOnz7Nxx9/TFhYGK+//rrJudTUVE6dOkVgYKDJ8e7du5OQkEB0dDQAx48fp1OnTsZkBNCjRw80Gg1RUVEWxSOEEKLomF1lB1C+fHmmTZtmtTf38fHh8OHDeHt7s3PnTpNz165dIz09nTp16pgcr1WrFqCv+mvatCk3b97Mco2Xlxfu7u5cuXLF4phOnz4NgEqlsvi5QghRWmk0GgCefvrpfL+GRQkJ4MaNG1StWhXQd+Nt374dlUpFcHAwNWvWtOi1ypcvn+O5xMREANzd3U2Ou7m5AfDo0aMcrzFc9+jRI4viychwc+2dDkhK02Y57uqob/wqFGD5dGYhhCh8Fi0dNGLECBwcHNi1axdxcXH079+f+Ph4ADZs2EBkZGS+tj7PTl5DW0ql0qxrLGVoGTVr1szi5166dAmAevXqWfzcgtgZHcvUXWdxVCpJ12oJbl6NnWeuGx9/2K8xwS2qF2pM5iqqe1acyT2znNyz/LHkvmVX8GYpsz+xFyxYQGxsLAMGDABg+/btxMfHM2/ePA4ePEiFChVYuHBhgQMy8PDwAODx48cmxw2tHg8PD2PLKPM1husMr1HSBbeozveTO7EuzI+dr/qz88x1UtK1JKaqSUnXMnXXWe4kphR1mEIIkSuzE1JUVBTDhg3jxRdfBODIkSNUrlyZ3r17U7NmTQYOHGgsNLCGmjVrolKpuHr1qslxw+M6derg5uZGpUqViImJMbkmLi6Ox48fZxlbKskqerjwdC1PktO0OGZqGToqlVy7n1xEkQkhhHnMTkgJCQnUqFED0E+CPXfuHAEBAcbzZcqUyXXirKWcnZ1p2bIlBw8eNOmaO3DgAB4eHjRq1AiAtm3bcvToUZP3PnDgACqVCj8/P6vFU1zU8HIlXWs6ppSu1VLDy7WIIhJCCPOYnZCqVq3Kn3/+CcDevXsB6NSpk/H88ePHqV7duuMUo0ePJjo6moiICL777js++eQTVq9ezahRo3B11X/Ajhgxgrt37zJy5EiOHj1qnBQ7YMAAY/FFaVLRw4UP+zXGxVGJh7MDLo5KPuzXGIDTMQ+k604IYbfMLmro3bs3n376KbGxsRw/fpzKlSsTEBDAtWvXmD17NkeOHLH6wqpt2rRh8eLFLFq0iDFjxlCpUiUmT55MWFiY8RofHx/WrFnDnDlzeO211/D09CQ0NJRx48ZZNZbiJLhFddo9WZ5r95Op4eVK1J/3aD/naLEochBClF5mJ6Rx48ahUCjYvXs3TZo04a233sLR0ZGEhASOHz/OyJEjGTZsWL4DCQ4OJjg4OMvxbt260a1bt1yf27JlS5PVHIS+pVTRw4U7iSlM3XWWlHQtKei78qbuOku7J8tT0cOliKMUtqLT6bh37x4pKSlotVmnBZQ2KSn6noHM480id6mpqSgUCh4/fmyccmNLFs1DGjt2LGPHjjU5Vr9+fU6ePGmyUoKwH9fuJ+OoVBqTEfxb5CAJqWTS6XRcv36dxMREnJ2dZZI3lJqKW2tzdXUlPT2dq1ev4uHhQdWqVfM1ncZcOSakffv20bNnzzxfQKVS5fgL/80339C7d+/8RycKLKciB1cnJadjHlDDy1USUwlz7949EhMTqVSpEl5eXkUdjl0wtJBcXOR33RKGFvbjx4+5d+8e8fHxFi8RZ4kcU91nn31GcHAwBw8etGgrh5SUFL788kv69OnDqlWrrBKkyL/sihyCm1cjeNmPDF/zX9rPOcrO6NiiDlNYUUpKCs7OzpKMhFUolUrKly+Pk5NTgVa/MUeOLaQdO3awYcMGpk6dCkBAQAAdOnSgfv36VK9eHTc3N7RaLfHx8dy8eZNff/2V06dPc/ToUVQqFaNGjSI0NNSmwQvzZCxycHVSErzsRxlTKsG0Wq100wmrUigUODg42Hw8MseEpFQqGTZsGMHBwWzfvp3Nmzezb98+FIrsV0bT6XTUrVuXcePG0b9//2zXlxNFx1DkcDrmgYwpCSHsUp5FDR4eHoSFhREWFsbly5eJjo7m2rVrPHz4EIVCQfny5alatSqtW7emWrVqhRGzKAAZUxJC2CuLqux8fHzw8fGxVSyiEBjGlDIvxhq87EeZpySEKFIWbz8hij8ZUxKiZNPpdDkOr9gz2xWUC7uW22KsKoWCw+fvyDJDwi5dunSJiIgI2rZtS6NGjWjXrh0TJkzg4sWLRR2a0c6dO6lXr57JV+PGjenatSszZszg/v37NnnftLQ0Zs2axe7du43HpkyZkufiAvZCElIpl92Y0uM0DTP3npeScGF3Ll68yKBBg0hISOCdd95hzZo1TJ482bg1jjX25LGm5cuXs3XrVrZs2cKnn37KwIED+frrrxk0aJBNktL9+/dZu3YtarXa6q9dGCQhlXIZ5ym5Of1bKvw4VSN7KYks7iSmFOkivZ9//jne3t6sXLmSHj164OfnR9++ffn8888pV64cy5YtK5K4cuLr60uzZs1o3rw5/v7+vPLKK3z++efcvHmTjz/+uKjDszuSkIRxg7+3e/vi5mw6f0X2UhIGO6NjaT/naJFOqI6Li0On02WZD+Pm5sbUqVNNVpcZMmQIb7/9NgsXLqRDhw60a9eOiRMnZmmZ/Pzzz4SEhNC0aVOeeeYZpk2bRkJCgsk1f//9N2PGjKFVq1b4+fnx6quvZtmrzVwNGjSge/fufP311yQn//v/Vl5x7Ny5E19fX06dOkVQUBBNmjQhODiY48ePAxAbG0uHDh0AeOutt+jcubPJ+27fvp3AwEAaN27Ms88+S1RUVL7it6V8JaS///6bLVu2sGbNGg4ePGhyU0XxVNHDha6+FdFoTbeFz1gSLi2l0ivjIr1FuRNx+/btiY2NZdCgQURGRnL58mXjuR49etCvXz+T6w8cOMD+/ft55513eP311zl+/DgjRowwJrSff/6Z0NBQ3NzcWLhwIZMnT+bYsWO8/PLLxm6v27dvM3DgQK5du8b777/PrFmziI2NZfjw4SQlJeXr5/D39yc9PZ2zZ8+aHQfoJz2PHTuWvn37snjxYsqVK8eoUaM4d+4cFStWZPny5YB+654lS5YYnxcbG8vq1auZMGECixcvRqfTMXbsWB48eJCv+G0lxyq7du3a8X//938mg2FarZb333+fbdu2mfyF8sQTT/D+++/TvXt320YrbEpKwkVO7GWR3pCQEO7evcvatWt5//33AfDy8qJdu3YMGTKEJk2amFyfkpLCmjVrjOuvVapUifDwcL7//ns6duzIvHnz8PHxYcWKFcZFQ319fenXrx979+6lb9++rFu3DrVazdq1a/H29gagbt26hIaGcv78eVq2bGnxz2F4nXv37gGYFQfoq+fCwsIYOXIkoN+ip2vXrqxcuZJFixbh6+sL6HfcNvwb9J/dK1asoHbt2oB+A9Thw4fz22+/GVtV9iDHFpJh6fqMli5dypYtWwgMDGTt2rV89dVXzJ8/n2rVqhEREcGPP/5o84CFbRm679aF+bHzVX92nrle5H8Vi6JnLzsRKxQKIiIi+OGHH5g3bx4vvPACbm5ufP311wwYMIDIyEiT659++mmqVKlifNyxY0ecnJw4deoUycnJ/Prrr3Ts2BGtVotarUatVvPkk09StWpV4+fZ6dOnadGihTGJANSuXZujR4/mKxllZm4cBn369DH+28nJiY4dO3L69Olc36NChQrGZAQYN1PN3DVZ1Cyah7Rx40YCAwP55JNPjMfq1atH165dGThwIEuXLsXf39/qQYrCJcsMicyyaz1/2K9xkf0ePPHEEwQFBREUFATA+fPneeONN5g9ezZBQUE88cQT+rgrVjR5nkKhwMvLi4SEBBISEowthxUrVmR5jzt37gDw8OFDatWqZdX4Da9dsWJFs+MwqFChgsljLy8v4uPjc30/ww7bBoY5Sva2V5bZCSk5OZn4+Phst6RwcnLiueeeM0lUoviTZYZERpl3Ii7s//a3bt3ihRdeYPz48fTv39/knK+vLxEREYwZM4bY2FhjQnr48KHJdTqdjri4OLy8vHBzc0OhUBAWFpbt55phQzp3d/dsS7SjoqLw8fExaYGZ66effsLV1ZWGDRui0WjMisMgPj7epLUWFxdn8rg4y7WoIeNMX1dXV8qXL5/jVhSPHz+WFYZLGNm6QmRmmFBdFH+IVKhQAZVKxaZNm0hNTc1y/u+//8bV1ZWaNWsaj0VHR5t0S3377bekp6fTunVr3N3d8fX15Z9//qFx48bGrzp16vDJJ5/w66+/AvpuvzNnzpgkt+vXrzNixAhOnjxp8c9x6dIlDhw4wHPPPYerq6vZcRgcPXrU+O/U1FSOHj1K69atAWy6eV5hyLWFtGLFCk6ePGlcw65z586sX7+eHj16mOwQe/HiRdavX2+V/lRhX2SZIWEvVCoV06dPZ9y4cTz//POEhITg4+NDcnIyx48fJzIykokTJ5rsDvv48WNGjhxJaGgocXFxLFy4kLZt2xo/wMePH094eDhTpkyhV69epKWl8dlnn/HHH3/w5ptvAhAaGspXX33Fyy+/THh4OAqFgiVLllC3bl0CAwNzjfn8+fPcunULgKSkJM6dO8eaNWuoWbMmERERxuvMicNg9uzZpKamUr16ddauXcvjx48JDw8H9K05hULBiRMn8PHxoWnTpgW/8YUox4Q0bNgw/vrrL7777ju2b99uPK5QKFi5cqVxK/Pw8HCioqJwcXFh4sSJto9YFLq8xpR+i43Hs0zRdOOI0qVLly5s27aN1atXs2LFCuLi4nB2dsbX15dPPvkkyxI5fn5+NG/enLfffhsHBweCgoJ4/fXXjec7dOjAqlWrWLJkCePGjcPZ2ZnGjRuzfv16nnrqKQCqVq1KZGQkc+fOZfLkyTg7O+Pv78/kyZMpU6ZMrvGOHj3a+G8nJydq1KjBgAEDeOWVV0wSpzlxGLzzzjssW7aMmzdv0qRJEzZs2ECdOnUAKFOmDKGhoWzZsoXvv//eOEepuFDodDpdXhclJiby559/Gr8Mm/UBDB8+HE9PTyZOnEiNGjVsHrCtGZYeadasmcXPvXTpEqAv9CiJ7iSm0H7OUVLS/01IDkpQKZU4qfJXFl7S75kt5HXPYmJiAKw+EF/cDBkyBJVKxbp160rEFuY7d+7krbfe4rvvvqNy5cqF8p4Z71tev1cF+ew0MKuowcPDgxYtWtCiRYss59atW5fvNxfFS+ZKqzSNBq0OUtVaUtXShSeEKBiLt5/QarXExcXlWC5YqVKlAgcl7FfGMaUHSWlEbPmFdM2/M8mlLFwIkV9mJ6SbN28ybdo0Tp48iUajyfG6CxcuWCUwYb8MY0p3ElOkLFzYrQ0bNhR1CFYVHBxMcHBwUYdhU2YnpOnTp/PTTz/RvXt3qlevXuzLC0XByVJDQghrMjshRUdHM2LECJNSRSGkLFwIYS1mN3PKli1L+fLlbRmLKKZy233WUBYuq4ULIfJidkIaPHgwW7dutbvF+IT9yG6poeR0NWMio2VlByFEnszusnvllVc4efIkXbt25emnn8bb29tkaSHQT5o1LAkvSh8pCxdCFITZCWnt2rXGHQYzrqWUkSQkkVdZuEqh4PD5O3T1rShJSQhhwuyEtG7dOho1asTMmTOpU6eOLKQqcpRbWfjjNA0z957n/T2/82G/xjR0y+FFhBBWo9PpsvRo2SOzE1J8fDyvvvqqLPEizJaxC0+lUPA4TT9/7XGq/vvUXWdZG1wTrzIWz88WRSH5AZz/Ch7dAfeK4PssuHoWaggZlwPKzfXr11m2bBlRUVHExcXh7u5O48aNeeWVV/Dz8zO59v79+6xYsYJvv/2WW7duUaZMGRo0aMBLL71kXBvv5MmTDB06NM/4jhw5Ytz8LqMpU6awa9cu42OFQmFcmfzZZ59l6NChODhY//+D27dvM336dN555x1jXJ07d6ZNmzbMnDnT6u9XUGbfgaZNmxrX0BLCXIYuvMPn7zBz73ljMgJ9Bd4f91Io66LCs2qKdOHZK50Ojs6EHxeDQgnpyeDoCvveBP9x0OltsKO/vm/fvs2AAQOoWrUqkyZNonLlyty+fZsdO3YwbNgwFi5caFylOzk5mcGDBwMwatQoatasSWJiInv37mXs2LFMnTqVYcOG0bBhQ7Zu3Wp8j2PHjrF8+XKWL1+Ol5eX8XjmDQEzqly5MgsXLgT0K94kJCQQFRXFxx9/THR0NIsWLbL6/M6ffvqJY8eO8c4771j1dW3F7IQ0bdo0QkND8fDwoEuXLpQvXz7bbjtZOkhkVtHDha6+FXl/z+8mx5PT1cw8dhsHpQLtwZsyidZeHZ0JJ5aCOkPZfnqS/vuJpfrvnacVflw52L59O0lJSaxbt864uV1KSgqdO3dmyJAhJgnpwIEDXLlyhcOHD5ssDt21a1dSUlJYtGgRL730Eu7u7iaLhv7999+AfmNAcxc6dXJyyrLwaMeOHfHx8eHdd9/l66+/5rnnnivQz17cWVT2nZCQwKeffsrAgQPp0qULHTt2zPJlC5s3b6Znz540a9aMPn368PXXX5ucj4qK4vnnn6dp06Z07tyZNWvW2CQOkX+ZN/tzdlCgUChI0+hISteSkq5l6q6zMlfJ3iQ/0LeMDAkos/Qk/fnkh9mfLwJxcXFA1u25VSoVkyZNYuDAgcZj9+7dA/RjLJmFh4fz6quvkpaWZsNoYdCgQVSpUoVt27aZHN+2bRu9evWiUaNGdO7cmZUrV5rEOWXKFIYPH05kZCQBAQE0b96ckSNHcvXqVUC/OvjkyZMB/bYdU6ZMMT43PT2dWbNm4e/vT7NmzXj55Ze5du2aTX9Oc5jdQho8eHCRDIpt3bqVd999l7CwMAICAvjuu+944403cHR0pGfPnkRHRxMeHk7Pnj0ZP348p0+fZs6cOeh0Ol5++eVCj1fkzJyFWWVvJTtz/it9N11uFEo4/yU8PbxQQspL+/bt2bRpE/3792fAgAG0adOGWrVqoVQqadu2LW3btjVe265dO+bPn89LL73EwIEDadeuHb6+vjg6OtKkSROaNGli83gVCgWtW7dmz549qNVqHBwc+PTTT1mwYAHDhg0jICCAs2fPsmjRIu7fv2+SWM6dO8fly5d5/fXXcXR0ZN68eQwdOpS9e/fSsWNHxo0bx+LFi1myZInJ+P/u3bsJCAhg9uzZ3Lt3j48++ohJkyZlSYqFzeyEVFRLBu3atYtnnnnGuGuiv78/586dY9OmTfTs2ZNFixbh6+vL3LlzAf0vo1qtZsWKFQwZMsRkZ1tR9HKrwDNMos3v3krCBh7d0Y8Z5SY9WX+dnejUqRPTp09n/vz5zJ49G9BvoePn50dISIhJQqpfvz7z58/ngw8+YNGiRSxatAhXV1datmzJ888/T8+ePQslZm9vb9LT03n48CHOzs4sW7aMkJAQ3nrrLUCfOMuUKcPs2bMZOnQoVatWBfR71a1evdq4M6yPjw/PPvssu3btIiQkxNgN2aBBA5NiiypVqrB06VIcHR0B/R5ay5cvJzk5GVdX10L5mbNj9yukpqamGvuBDcqVK8fDhw9JTU3l1KlTWbYR7t69OwkJCURHRxdmqMIChi48J5WCMo5KYxdeqlpLYqpauvDshXtFfQFDbhxd9dfZkZCQEKKioliyZAkhISFUrFiRI0eOEBYWZvzj1aBHjx4cO3aMVatWERYWRt26dTl+/DgTJkwgIiIi2+48Wzpz5oxxzEutVhu/OnfujEaj4aeffjJeW6NGDZNtyuvVq0etWrU4depUru/RrFkzYzICjMmqqFfisft626FDh/LOO++wb98+AgICiIqK4tixY0RERHDt2jXS09ON2/caGHY0vHLlCq1bt7b4PTUaTb4qCpOS9P3sUo1onoZusKJ3RW4/1pCGA7O/v0O65t//+ZXAvp9+p6yLisrujlIe/j95/Z6lpKTg4eFh3O2zQHx64KybTG6d9TqdhlSfnmCN98uDVqtFoVCY9bMpFAoCAgIICAhAq9Vy9epV3nvvPVatWkXv3r2pW7euyfWtWrWiVatWANy5c4dZs2axd+9egoKCTFpVoB+DAf0fzObEotFo0Ol0OV5748YNnJyccHFx4e7duwCEhYXleG1KSgoajYby5ctneU1PT08ePHhASkpKtnHqdDqcnJxMnmfYUiglJcUCPaZAAAAgAElEQVTkuGEczvB+iYmJOf4MGo2mwPNT7f7/8N69e/PTTz8xYcIE47F+/foxYsQIzpw5A4C7u7vJcwwtqkePHhVeoCJfPF1VeLqqSMEJtdb0L9FUtdZYhafW6hjvX4Gu/ylbRJGWUi7lUPu9isN/V6BQZ+260zm6om4VDi5PFEFwWWk0GoKCgujTpw+vvvqqybmaNWvy5ptvMnDgQP7++2/q1q3LkCFDqF27Nh988IHJtRUrVuTdd9/lyJEjXL58OUtCsnbMp0+fpkmTJjg4OBg/z2bPnm1S+WdQoUIF47/j4+OznI+Li6Nx48Y2i9eW7D4hjR49mjNnzvDWW2/h6+vLr7/+yrJly3B3d6dXr165Pje/Nf0qlSpfE4ANf7HK5GHzGe5Z83r1mKV7Iss6eGkaHWn/azUt/uke/ds3LvXFDnn9nsXExADg4mKl+9Tt/8DBIes8JJ0WRZuxOHZ6G8dCKnhSKpUolcpcf7aKFSvy1VdfERoaiqenfuKu4a/6GzduoFAoaNiwIS4uLtSoUYPDhw/z2muvZfnwN9xnw7UZGbq7nJ2dzbrPKpUKhUKR7bVbtmzh9u3bTJkyBRcXF1q1aoWjoyMPHz40KQM/c+YMixcvZuLEidSsWROVSkVMTAy3bt2idu3aAFy8eJFr164RHh6Oi4uL8f0yxqlQKFCpVCax5PTzGO6bi4sLKpWKcuXKGXugMvvll1/yvA95seuEFB0dTVRUFB999JFxp0Q/Pz/Kli3L9OnTeeGFFwB4/PixyfMMLSMPD4/CDVgUiFTh2SmFQj/PqM2YTCs1PAeu5Qo9nJs3b2a7UoOvry9+fn68/fbbDBs2jODgYIYOHUqDBg1ISUnh9OnTbNy4kcGDBxu7+SMiIjh58iQvvPACQ4cOpXnz5iiVSs6ePcuaNWvo1KmT1VpHaWlpxg9trVZLfHw8P/74I5s2baJ3797GP7C9vLwICwtjwYIFPHr0iKeffpobN26wYMECPDw8ePLJJ42vqdVqGT16NBMmTECtVjNv3jz+85//0KdPH+Dfz8BDhw7Rvn17fHx8rPKz2EqOCen27dv5ekFrToy9ceMGAC1atDA53rJlS0C/XbpKpTLW3RsYHmceWxL2T6rw7Jirp12Udv/zzz989NFHWY4PHToUPz8/mjRpwq5du1i5ciUbN27k7t27KJVK/vOf/zB16lTjH7KgLwowXLt7924+++wzdDodtWrVYsSIEWYtF2SuW7dumcyBcnNz48knn2TatGkmx0GfKCtUqMCmTZtYsWIF5cqVIyAggIkTJ+Ls7GwS/+DBg3nvvfdIS0ujY8eOTJ061Vhd7OfnR9u2bZk3bx4nT55kxYoVVvt5bEGhy6GEpH79+vmad3ThwoUCB2Xw+++/ExwczIcffsjzzz9vPP7FF18wdepUtm3bxrx581Cr1URGRhrjnTt3Llu3buWHH36wuITR8BdM5hnV5pAuO8vlds92Rsdm6cLLWPTg4qhk56v+JKdpS1WLydwuu5y6VkqjjF1PJcWUKVM4ffo0hw4dstl7ZLxvef1eFeSz0yDHFtKoUaNMEpJarWbTpk04OTnRq1cv6tata6xc2bNnD0qlkvDw8HwHkp2GDRvStWtXZs6cSWJiIg0aNODcuXMsXbqU9u3b07RpU0aPHk1oaCgRERH069ePM2fOsHr1aiZNmlSk9fSi4PLqwtNpdTy35EecHaTFJERJkGNCyjwRdu7cuXh6erJt2za8vb1Nzr322msMHDjQmEGtacGCBSxZsoR169YRFxdHtWrVCAsLY+TIkQC0adOGxYsXs2jRIsaMGUOlSpWYPHlyjiWTonjJrQsvVaMDdKRpZPM/IUoCs4sadu7cyYgRI7IkI4CyZcsycOBAli9fbvVVZZ2cnJg4cSITJ07M8Zpu3boZl4kXJVPm3WhT1RqUSgUp6f8mKSl6EKXJrFmzijoEqzM7IanVauMkq+w8fGg/iyuKkiljF56rk5LgZT+anJeiByGKN7Mn6vj5+bF27Vr++OOPLOdOnDjB559/TocOHawanBCZVfRw4elanvhWeSLb1cMzLj301s7fOHzhdqlbfkipVBpn3gthLYZVMmzJ7BbS5MmTGTRoEP369aNp06bUrFmTlJQUYmJiuHjxItWqVTMudS5EYcir6CFVreO1zWfQ6nSlqrXk4OBAUlISWq3W6hu+idJJo9GQkpKS7ZCNNZmdkGrVqsXXX3/NypUr+eGHHzh79iwKhYIaNWowatQoRowYkWUJHyFsLbeiB4CktH+3S69fxaNUlIiXLVuWhw8fEhcXR/ny5Ytk2xhRcmg0GuO+UbZebMCilRoqVKjA22+/batYhMi3jEUPShQkpZt2WZWmEnE3Nzc8PDy4d+8eCQkJODjY9YIshcLQhVnQxT9Lm/T0dFJSUlCpVFSqVMnmU2ks/k3966+/OHbsGLdu3WLw4MGUKVPG5osPCmEOQxfeb7HxjImMJlX9b4uptJWIV61alfj4eB49epRl59TSKDExEdBvXSPMl5SUhEKhoHbt2oUyr9OihPThhx+yYcMGdDodCoWCLl26kJiYyGuvvUaXLl1YsGCBbIgnilRFDxe6NnDho+DSXSKuVCrx9PQ0Li5a2hlWHJDVKyxjWBWksBYZMHvEc+PGjaxfv56wsDA2b95s3LTKz8+Pl156iSNHjvDZZ5/ZLFAhLBHcojrfT+7EujA/vhybtfVuKBEfvua/tJ9zlJ3RsUUQpRAiI7MT0ubNmwkMDOSNN94wLnUO+ibwtGnTCAoKYvfu3baIUYh8saREfOqus5y/Gc/pmAelrkxcCHthdpfd1atXGTJkSI7n/fz8bLrInxAFIeviCWH/zE5I5cqVy3VLisuXL1O2rOzmKeyXpevilZYycSHshdlddl27dmXTpk1cvnzZeMwwvyEqKootW7bQqVMn60cohJUZSsQNXXhOKgUujqb/KxhaTDLGJEThMbuFNGHCBE6ePElwcLBxr6Rly5bx8ccf8/vvv1OlShVee+01W8YqhNXktS6etJiEKHxmt5CeeOIJtm/fzvDhw0lMTESlUnH69Gni4+MZOnQoX3zxBeXLl7dlrEJYVU5FD9JiEqJomN1CSktLw83NjYiIiCx7JQlR3EmLSYiiZ3ZCatu2LYGBgfTp04fWrVvbMiYhioSh6AHIc+8lqcoTwvrMTkgdOnRg//797Ny5kwoVKtC7d2+CgoJo2LChLeMTokhY2mJ6a+dvlHV1pEn1J6S1JEQ+mZ2QPv74Y9LS0jh27BjffPMNW7ZsYd26ddSuXZu+ffsSFBREjRo1bBmrEIXKkhZTad3qQghrsmgtOycnJwIDAwkMDCQ5OZlvv/2W/fv3s3LlShYtWkTTpk3ZsmWLrWIVosjk1WKC0rnVhRDWlO916V1dXWnVqhUJCQkkJCRw8uRJLly4YM3YhLAr2bWYzN3qwpDMJEEJkTOLE9KdO3c4cOAA+/fv58yZMwC0bNmSDz74gB49elg9QCHskSVbXUze8SsqpRInlRRACJEbsxPSxo0b2bdvH2fOnEGr1dKgQQMmTZpEUFAQlSpVsmWMQtglc7e6UGtBrdUak5Z06QmRPbMT0owZM6hWrRqvvPIKffv2xcfHx5ZxCVFsmDO+lJGUjAuRPbMT0qZNm2jRooUtYxGi2MqpIi9No0Grg3SNznitTLIVIntmJyRDMvrxxx85cuQI169fx8nJicqVK9OpUyfatGljsyCFKE4ytphqeLkS9ec9iyfZShGEKI3MTkharZbJkyfzzTffoNPpKFOmDDqdjuTkZDZs2ECvXr2YN2+eLWMVotjI2GKydJKtFEGI0srshLRmzRr27NnDgAEDGDNmjLGQ4fbt23z66ads3ryZJk2aMGzYMJsFK0RxZckk27yKIIQoqcxOSF988QWdO3fm/fffNzleqVIlpk+fzu3bt9mxY4ckJCHyUNAiiNAWnjxV3gXPqinSnSdKFLO3n4iNjaVdu3Y5nm/bti1Xr161SlBClHQ5bX3h7KDAUaUwuTZVo+/OS0xVk5KuZfnJOKYevGHcBuNOYgqnYx5wJzGliH4aIazDoi3Mr127luP5mJgYPDw8rBKUEKWJpUUQAMnp+nEnGW8SJYnZCalTp05ERkbSunVrOnToYHLuu+++Y8uWLfTt29fqAQpRGlhSBJGRTLoVJYnZCSkiIoITJ04QHh7OU089RZ06dQD4+++/+fPPP6lSpQrjx4+3WaBClCbZFUGoFAoep2lyfZ6UkIvizOyE5Onpyfbt2/n000/59ttvOXz4MABVq1Zl6NChhIeH4+npabNAhSitMraYzsY+5MO9F3BQKtDodHlOus2uS08SlLBXZiekL7/8kqeffpo333yTN99805YxCSEyMbSYnq7lSb0ySdx6lE7bZvXzHG/K3KUnCUrYM7MT0vvvv09YWBhjx461ZTxCiDx4lXHAq4wDFT1cLC4hlwQl7JnZCcnV1RVX16KZlPfzzz8zf/58zp8/j4eHB927d2fixIm4ubkBEBUVxYIFC/jrr7/w9vbmpZdeIiwsrEhiFaKwWbKOXmZ5JagpPerTuHo5SU6iUJidkKZPn87777+PSqWiZcuWeHt7o1RmncZk7a0ofvnlF0JDQ+ncuTPLly8nJiaG+fPnc//+fRYsWEB0dDTh4eH07NmT8ePHc/r0aebMmYNOp+Pll1+2aixC2LvcSsjzk6De3X0eN2cVGq1OWk/C5sxOSFOmTCEtLY3Zs2fnep21d439+OOPadasGQsXLkShUODv749Wq2Xt2rUkJyezaNEifH19mTt3LgDt27dHrVazYsUKhgwZgpOTk1XjEcLe5VRCnp8EBfA4VV/ZJ917wtbMTkhDhw5FoVDkfaEV3b9/n1OnTjFv3jyT9w4JCSEkJITU1FROnTrFhAkTTJ7XvXt3Vq1aRXR0NK1bty7UmIWwN9ZKUDL+JGzNonlIhe2PP/5Ap9PxxBNPMGHCBI4dO4ZKpSIoKIi33nqL2NhY0tPTjXOiDGrVqgXAlStX8pWQNBoNly5dsvh5SUlJAPl6bmkl98xy1rhn7sCDG9DQDdYG1+TWo3QquzsSfSOJhT/eRaWAZHXuLafMCeqN7b+iUipwUCpQa3WM969Ai6pljK/tVcbsjxurk9+z/LHkvmk0GlQqVYHez+LfkKtXr3L06FGuX7+OSqWidu3adOzY0SbbmN+/fx/Qdxd269aN5cuXc+nSJT755BNSU1MZOHAgAO7u7ibPMxQ7PHr0yOoxCVHSGKr2ALr+p6wxifx5N4VVp+/joFSQrtGiQ5+EcqLRgUajI+1/Laz5UXfsNkEJ+2T2b4ROp2PGjBls2bIFjcZ0trijoyNjxowhPDzcqsGlp6cD+s0B/+///g+ANm3aoNPpmD17NgMGDMj1+dkVXZhDpVJRr149i59n+CsiP88treSeWa4w79nwbin5Hn/KnKAWHL9TZBV88nuWP5bct19++aXA72d2Qlq5ciWRkZH06tWL4cOHU7t2bTQaDf/88w9r165l4cKFVKpUiX79+hU4KANDS6d9+/Ymx9u1a8esWbM4e/YsAI8fPzY5b2gZyWKvQhSMNQskpIJP5MXshLRt2za6devG/PnzTY57enrSvHlzRo0axZo1a6yakGrXrg1AWlqayXFDy6l69eqoVKos214YHmceWxJCFExRVvDdSUyRhFXCmZ2Q7t27l2uBQIcOHZgzZ45VgjLw8fGhWrVq7N27l8GDBxuPHz16FAcHB5o3b07Lli05ePAgw4YNM1biHThwAA8PDxo1amTVeIQQpgqrgi+4eTV2nrmOo1Iq+koysxNS48aN+eGHHwgJCcn2/JkzZ3jqqaesFhiAQqHg9ddfZ+LEibz++usEBwdz7tw5li9fzksvvYSXlxejR48mNDSUiIgI+vXrx5kzZ1i9ejWTJk0qspUlhCitzElQ5qxanjlBbfqvfi+2FMwvOc/YohLFg0Kn0+Xepv6fy5cvExoaytNPP83IkSOpW7cuDg4OxMbGsmHDBnbs2MGyZcvw8fExeZ41qu8OHz7M0qVLjUsDDRw4kFGjRhmLFg4dOsSiRYu4cuUKlSpVIiQkJN9LBxkG5po1a2bxc2Xg1HJyzyxXnO+ZIUmcjX3IrP0XLerey42DklxbVIZt39s2qy+tKQvkp6ghP5+dBmYnpCZNmqDValGr1cauMYVCgU6nw/AS2U2ctfbKDbYmCalwyT2zXEm5ZxlbMPkZf7KUq6MCHQoZn7JAYScks7vsQkNDC32lBiFEyWXu+FPmFk9+E1ZO276bMz4lCatw2PVKDUKI0iO3BFXRw4UJ3Z6ySovK0vEpSViFx+Kp0tevX+fu3btotdlP2W7RokWBgxJCiIwJKvNjc1tU5hRQZFbQhJV5wm/mZCXJK2dmJ6RLly4RERHBlStXsj2v0+lQKBTFbsxICFE8mduiymvbd0vllbAyTvjNnKzy09oqTQnM7IT03nvvcfPmTUaMGEH16tXzvSyPEELYQk4tqty2fbfW+FRmhgm/mZOVpa2tvBIYkGuyKm7JzOyEdOHCBUaPHs3IkSNtGY8QQlhdTtu+5zU+ZauEZZBXayu3BJacrq94dnFQZZusMv8cxaE1ZnZCqlixIo6OjraMRQghCoUl41PmJqz8jFdZKnMCAx3pGjVgmqwyJs6CtMYautn0x8nC7IQ0cuRIFi9eTKdOnYxrzAkhREmU34SVccJvXh/6tm5t5XU+r9bY1F1nWRtcs1C3CTH7nfr168e+ffsICgqibt26eHt7Z5mXpFAoWL16tdWDFEIIe5LbeFWvJlVybF1Z0j1o6wSWF0elkluP0u0zIS1YsICoqChUKhX3798nPj4+yzUycVYIUdrl1rrK/Div7sHcEljGMaTsklXG5ZTyk8zStVoquxfuMI3ZCWnHjh20bduWBQsWULZsWVvGJIQQpUZBEhhgdhGDJa0xw/O9ypjuNWdrZiek1NRUunfvLslICCGKUHYJC7JvbeV1PrfWWEUPF+NadoXF7ITUpk0b/vvf/+a5bbgQQoiikTlZ5XU+r8eFzeyENH78eF555RUmTZpE165dKV++PCqVKst1snSQEEKI/DA7IfXt2xeAb775hr1792Y5L0sHCSGEKAizE9IHH3xgyziKv+QHPPHXlzikxMEjX/B9Flw9izoqIYQoNsxOSP3797dlHMWXTgdHZ8KPi6mgA6UmFS66wr43wX8cdHobpBxeCCHylGNCio6Opk6dOnh6ehofm6PUjSEdnQknloI6BeOIWnqS/vuJpfrvnacVRWRCCFGs5JiQBg8ezNy5c+nTp4/xcW4TX0vlGFLyA/hxMahTsj+fngRRC8DVG5oNki48IYTIRY4J6YMPPjDZG13GkLJx/itQ5LENh1YNh6fDkXelC08IIXKRY0LKPGYkY0jZeHQH0pPzvk6Tpv8uXXhCCJEj2WWvINwrgqOr+dcbuvBOLNd39wkhhDCShFQQvs+CLvul3nNk6MKbVx++naGv0hNCCCEJqUBcPfXjQo5lLHueJk1fCHFiqb5KTwghhCSkAuv0NrQZAw4uaBUOWNTeSU/SV+klP7RVdEIIUWxYNSHdvHnTmi9XPCgU+iKFSRe522wcKLKu75crrQb2TZYxJSFEqWd2QmrUqFG2a9gZ7Nixg6CgIKsEVSy5evKw/ovE+Q6zrAtPmw7ndsqYkhCi1Mux7Pv27dv897//NT5Wq9WcOnUKjUaT5VqtVss333xjmwiLmbjGoyjv7a3vitNp/y35zo02Xf8lZeFCiFIsx4Tk6enJkiVLiImJAfTbk2/atIlNmzbl+GIhISHWj7C4MXThtRkDv2yGQ+/oK+vMISs7CCFKsRwTkpOTE6tWrSI2NhadTkdYWBgjR46kTZs2Wa5VKpV4eXnx5JNP2jTYYsXVE9q8Csn39S0fw/p2eZGVHYQQpVSuq33XqFGDGjVqAPqlg1q3bm18LMzU6W399x8X6wsYtOl5P0dWdhBClEJmFzX0798/12SUlJTEt99+a5WgSpQMVXg0Cgalo/nPlZUdhBCliNn7IT1+/Jj/+7//44cffiApKQldhmqwjIUOpWq1b0u4ekLP2foFWc1pJRlIF54QopQwu4U0d+5c9uzZQ40aNWjZsiVqtZru3bvTokULVCoVzs7OLFiwwJaxFn8FXdkh6hPYNEBaS0KIEsnshHT06FE6d+7Mjh07mDdvHgChoaGsX7+eLVu2oFAouHbtms0CLTEyrOyAysmy52rT4c+DMK+ezFkSQpQ4ZiekuLg4AgICAPDy8qJixYr8+uuvgH7S7AsvvMDu3bttE2UGY8eOpVu3bibHoqKieP7552natCmdO3dmzZo1No8j3zKOKXV9D5Rm95r+S50q6+AJIUocsxOSm5ubybhRrVq1uHTpkvFxvXr1uHHjhnWjy+Srr77i0KFDJseio6MJDw+nbt26LF68mD59+jBnzhxWr15t01gKzFAW3i7C8i48kIIHIUSJY3ZCaty4Mfv370er1W+34OPjw+nTp43n//nnHxwdLaggs9Dt27eZOXMmlStXNjm+aNEifH19mTt3Lu3btyciIoKXX36ZFStWkJZmxioJRa1AXXiylYUQouQwOyGFhoZy6tQpevbsSUJCAv369ePy5cuMGDGCGTNmsH79elq1amWzQKdNm0bbtm1NJuampqZy6tQpAgMDTa7t3r07CQkJREdH2yweqyloF55sZSGEKCHM/vRr27Yty5cv5/PPP8fDw4OmTZsyceJEli5dSlRUFA0bNmTq1Kk2CXL79u38/vvv7Nmzhzlz5hiPX7t2jfT0dOrUqWNyfa1atQC4cuUKrVu3tvj9NBqNSXekuZKS9Ksx5Oe5AHh1wbvBn3hd3IRSk2LZc9OT0P4wn7uPNCTU6YXWqWz+YihkBb5npZDcM8vJPcsfS+6bRqNBpbJwt4NMLPpzvEOHDnTo0MH4eOTIkQwdOpSkpCS8vLwKFEhOrl+/zkcffcRHH32U5T0SExMBcHd3Nznu5uYGwKNHj2wSky3FNR4FgNeFDSi06Vgy40ih01Dhl8VU+GUJ9xu8pH8tmbMkhCgm8lHiZcrFxQUXFxdrxJKFTqdj6tSpdOjQge7du2d7PjdKZf62e1KpVNSrV8/i5xn+isjPc03UXwDJ02HnSLh81OyJtApA8b9ry1/cSPlHl6DJAP1W63a6UKvV7lkpIvfMcnLP8seS+/bLL78U+P3sesfYyMhILl26xNSpU1Gr1ajVamMSUqvVeHh4APpVJDIytIwM54slV08YvA3aTchfwYMmDWKOw/4pUvQghCgWCtxCsqUDBw7w4MED2rVrl+Vcw4YNeffdd1GpVFy9etXknOFx5rGlYqcgW1kYpCfrv8tCrUIIO2fXCem9997L0vpZunQpFy5cYMmSJVSvXp19+/Zx8OBBhg0bhuJ/4yUHDhzAw8ODRo0aFUXY1pffrSwykr2WhBB2zq4TUt26dbMcK1euHE5OTjRu3BiA0aNHExoaSkREBP369ePMmTOsXr2aSZMm4erqWtgh21bGrSzM3Y02I1moVQhhx6w2hlRU20+0adOGxYsXc/nyZcaMGcPu3buZPHkyr7zySqHHYnPWWHbIMG/p+EJY1xtOr5OVHoQQdsHsT7RHjx7x7rvvFvn2E7NmzcpyrFu3blnWtyvRrNGFZyh6uBEN+96UFpMQosiZ3UL6+OOPZfsJe5Nx2SGHfHZPpidLi0kIYRdk+4niLGMXXs9ZUKut5eXhBlImLoQoYsVu+wmRDVdPeHo4DP8G2o7P37wlA0OLSTYDFEIUsmK1/YTIgzWKHgxkM0AhRCErNttPCAsUdK+ljNSpMr4khCgUxWb7CZEP1ih6ABlfEkIUCrMTkmH7iWrVqplsP/Hzzz+zceNGnnzySZttPyHyyZpFDyAVeUIIm7L77SeEFRiKHloM02/i9+NiQAHq5Py9nsxhEkLYgMWj3g8fPuTEiRPcvHmTHj164OLiQnx8vCSk4iDjYq3nv4LftkHsz5YvQWRgWLj1+EKI+dHut7oQQtg3i5YO+vzzz+nYsSMRERHMnTuXmJgYzpw5Q69evZgxY0ae+xMJO5G5TFzlXLDXkzEmIYQVmJ2Q9u7dy0cffUSnTp2YP3++Mfk0aNCATp06ERkZSWRkpM0CFTZgaDG9fgmeDARlAaskZYxJCFEAZiekVatW0aZNGxYsWECbNm2Mx6tWrcqyZcsICAhg69atNglS2FjmzQALUpEH0mISQuSL2Qnpr7/+omvXrjme79KlC7GxsVYJShQBa1fkgbSYhBAWsWilBsPW4Nm5ceMGLi4uVglKFKHsliGyVotp35sw9z+wthecWifJSQhhwqJ5SJs3b+b+/ftZzv35559s2rQJf39/qwYnipAtWkzqFP0mgTHHYe/rMK8+3r+tkO48IQRgQUKaNGkSarWaoKAg3n33XRQKBTt27GD8+PEEBwejVCoZP368LWMVRcEWLSbQr5WnTsHr/HqqHxkt3XlCCPMTUpUqVfjiiy/w9/fnhx9+QKfT8c0333Ds2DE6dOjA1q1bqVmzpi1jFUXJFi0mQKlTU+ZutBRACCHMnxh74sQJmjRpwscff4xWq+XevXtotVq8vb1lUdXSxNqrPuifLZNshRDmt5AmTJjAihUr9E9SKqlYsSKVK1eWZFRa2ajFZFIy/nE9iOwPx2ZLl54QpYDZLSStVkulSpVsGYsojmzQYgL+bTH9eVD/5VgG9k6Guh2gWkvwqCQtKCFKGLNbSOPHj+ezzz7j8OHDPHggf6mKTGzVYjJITwJNqj45HfsQ9r8lY05ClDBmt5C2bt1KQkIC48aNA0ChUKBUZs1n586ds150ovjJrsWkVeu/rCk9Sf9dxpyEKDHMTkj169enfv36toxFlCQZVxb//Us4ux2undRvq65Osd77ZNwKQ9nE6e0AABonSURBVLr0hCjWzE5Ic+fOzfOa7CbNilLO1RNahuq/kh9YZ9uL7GQec3JwhW8mQY1noPEAaCjJSQh7Z9H2E9lRq9UcOnSI0aNHm2zeJ0QW2Uyy1SidsckIkDo5y6oQMt4khH2zeIM+gwsXLrBz50727NnDw4cP0el0VK9e3ZqxiZIqQ3fe3W9X4vHPftzun7Nuiykjbbr+S8abhLBrFiWkBw8e8PXXX7Nr1y4uXboE6OckderUiUGDBhEQEGCTIEUJ5epJ/H+eI97nWerd2KEvgFAo/y1YsDYZbxLCruWZkDQaDceOHWPnzp18//33qNVqdDodNWrUIDY2lrlz59KrV6/CiFWUVJm3Vk+8DddPwZXvQKGyfoKSOU5C2KUcE9LFixfZtWsXe/bsIS4uDoBGjRoRGBhIYGAgDg4OdO3aFWfnAm5/LYSBYYzJwFAEkTFBoSz4pNvMDAlPEpQQRSrHhPTcc8/h4uKCn58f7du3p2vXrlSuXNl4/vr164USoCjFckpQtqjSy0gSlBBFIseE5OTkRFpaGvfu3ePq1av8/vvveHp6SotIFJ3sJt3acszJQBKUEIUix4T0448/cujQIXbv3k1kZCQbNmzAxcWFgIAAunXrxpNPPlmYcQrxr9zGnLQa668KkVnmBCVznoSwihwTkru7O/369aNfv37cu3ePPXv2sGfPHg4ePMihQ4dQKpUoFArOnTtHQEAATk5WXLdMCHNk16Vny1UhcmIY04o5Dtf+C/uk9SREfih0OstmCsbExLB79252795NTEwMCoUCDw8PgoKCCA4OplGjRraKtVD88ssvADRr1szi5xpK4evVq2fVmEoym92zwhpvyotjGX2rLWOCqtsR/j4Gj+6Ae0WLE5b8nllO7ln+WHLfCvLZaWBxQsro7Nmz7N69m71793Lv3j0UCgUXLlzIdzDZ0Wq1bN26lU2bNhEbG4u3tzddunRh3LhxuLu7AxAVFcWCBQv466+/8Pb25qWXXiIsLCxf7ycJqXDZ/J7pdIU73pQXpYO+S1HpqP/u6Jo1YeWRoOT3zHJyz/KnsBNSvldqAGjcuDGNGzfmrbfe4sSJE+zevbsgL5etVatW8cknn/Dyyy/Tpk0brly5wqJFi/jrr79YvXo10dHRhIeH07NnT8aPH8/p06eZM2cOOp2Ol19+2erxiGKmsOc45cUwvqVN13+XggkhjAqUkAwUCgX+/v74+/tb4+WMdDodq1atYuDAgUyaNAkAf39/PD09iYiI4MKFCyxatAhfX1/j4q/t27dHrVazYsUKhgwZImNbQs+cOU5FkaAyyytBOXtQ7u4dVOlJ8MhXkpUoUaySkGzl8ePH9O3bl549e5ocr1u3LgB//vknp06dYsKECSbnu3fvzqpVq4iOjqZ169aFFq8oRoprggIqGs5dlNaUKFnsOiG5u7szbdq0LMcPHz4MgK+vL+np6dSpU8fkfK1atQC4cuVKvhKSRqMx9p1aIilJ/+GRn+eWVnZ1z9zbgDtQ5TmUTRLwuPotDsn3cIn7nTK3f0anUKLU6Kv2FEUYpvG9/5esdP9LVlqVC4o9r/O4sh+p3g1Ru5bnceVWuN36GYeUONQu3iTW7IzWqWyRxV5U7Or3rBix5L5pNBpUKlWB3s+uE1J2fv31V1auXEnXrl1JTEwEMBY3GLi5uQHw6NGjQo9PlAxap7LE/+c542NlWtYEhU6DQqcp0uQE/yYo1f+SpfvN47jfPI5OoUKh06BTOKDQadCqnKl4aq4kLGG3rJqQ0tLSbDpmc/r0acLDw6levTozZszgypUruV6f3Rbr5lCpVPmqxpFKHssVq3vWuNW//84850nl+O+irUXMkKAUOs3/vusLKQwJy+PmcTxuHs+24q9y9DyT8SoA0h7lqzzdnhSr3zM7kp8qu4IwOyF16dKFqVOn0qVLl2zP79mzhw8++ICTJ08WOKjs7N27lylTplC7dm1WrVqFp6cn9+7dA/RjTRkZWkYeHh42iUWIbHfCtbfxp7zkVfGXUXbVf3U7Fmg+lRCZ5ZiQ7t+/z+XLl42Pr1+/ztmzZylbNmtzXqvVcujQIdLSbDP5cO3atcyePRs/Pz+WLl1qTDQ1a9ZEpVJx9epVk+sNjzOPLQlhE5YUSJi0StKLLGSLZU5W2c2nkoQlCijHhOTs7MykSZO4e/cuoC/t/vTTT/n000+zvV6n09lkX6Tt27cza9YsevXqxezZs026BJ2dnWnZsiUHDx5k2LBhKBT6zooDBw7g4eFR7FeNEMVUTgnK8MFctxP8fbT4tagyyqt1lVfCytwdWLejJC+Rc0Jyc3Nj+fLl/PHHH+h0OqZOncqAAQNo3rx5lmuVSiVeXl60adPGqsHFxcUxc+ZMqlWrRkhICOfPnzc5X7NmTUaPHk1oaCgRERH069ePM2fOsHr1aiZNmoSrq6tV4xEiXzInKCgeJecFYUl3oLS2xP/kOobUsGFDGjZsCMCNGzcIDAzkqaeeKpTAAH744QeSk5O5fv06ISEhWc7PmTOHZ599lsWLF7No0SLGjBlDpUqV+P/27j+sxvv/A/jzVI5+kSRD0RSShRQZLvo5acIyymK0ia2LRZk+1VrfDbF0tYZcLnFJIj/CdjVNrU4/zGUul024rsK1y6hJbKU4pc6pc3//aJ05Tqfuc7rP6XTO6/Hf/X7f55xXLzmv3vf9fr/v6OholbcOIkTjlHwwYedeX309u48zvR1tdVWwXj3mm9Ni4n6iV3vZ6SLay06zKGcsvHrJj2+OJ3//DUOREMNaH+jGaKq3Xr8v13nMMwL+nWEoLeJsNrvt6ViPCppW72V38uRJ/Pjjj/jnn3/Q3t4u18/j8aSLVgkhHHltBNXw75fEMEdH3ZxAoazXR1idx8x/z8V6fTFx16MvcQ/HLEZn/H/XRNK9MZWwLkhpaWlIS0uDhYUFxo4diwEDBqgzLkIIG72ZQKEvBas7ioqZomNFBe2V0ZhUby819nT8evHTgWLHuiCdPXsW7u7uOHz4MG1YSoi2YjuBoruC9cr9KtKDLkZjcn2s7o31NDrr+lKkFJt1YsoeT1rCXZ5YYl2Q6uvrsWHDBipGhPRnyhQsvjkAHtD6nEZX6qDs6Ky74sdVsXt1NHfxf7ByDEHd5E+4/9kVYF2Qxo8f3+NWPYSQfq6rgtWJLgf2H6oWu9dGc0PvZHccT0xVY7D/YV2QNm/ejMjISMycORMeHh7qjIkQoo3UdTmQipfWMmhvwdA7x4GX/weYDFH757EuSJmZmTA1NcWnn34KY2NjWFpaSndG6ESz7AjRY6pcDhS9oNGWlmNgAFT8oHjkzCHWBam1tRV2dnbSZw0RQohSurscCHA3W7C7dUjq+tl0mEF7a8e/gwawLkhZWVnqjIMQQv6j6uXBro755niqaDGxsjf/9ZDEcCAMzYf3fCIHlH4eUltbG27fvo3Hjx/D3d0dxsbGaG9vh4WFhTriI4QQeT0VrNeOu1xMzKaYsR2d9bgOqf8WMx4kwKT3ej6RA0oVpIsXLyIxMRF1dXUAgCNHjkAkEmHTpk3YuHEjwsLC1BIkIYRwRsliJnfcXUHrzb2x3qxDUhOJoTHqJ4ZgmAYmNABKFKTLly9jy5YtcHV1RVhYGL755hsAgK2tLSZMmICUlBRYW1tjyRLNL6YihBCNYVPQFPWpOjrrrvixWSemyjokRoL6f9chDeM+i11iXZD2798PZ2dnHDt2DI2NjdKC5ODggOzsbKxevRqZmZlUkAghRJHejs66Ou7EZbEzHw5Meg91VU9U/lFVwbogVVZWIjIyEgYGBvJvYmSEgIAAJCcncxocIYQQltRS7DRbkOSriwIDBgxAW5vi65YNDQ204SohhBCVsS5I7u7uOHv2LFpbW+X6nj59iuzsbLi5uXEaHCGEEP3B+pJdVFQUgoODsXjxYsybNw88Hg8CgQClpaX4/vvvIRKJEBERoc5YCSGE6DDWIyQHBwecOHECw4cPR1ZWFhiGwfHjx5GZmYkxY8bg6NGjcHJyUmeshBBCdJhKjzBvaGhAVVUVJBIJbGxsYG1trY7Y+sRvv/0GADA0NFT6tZ1P0VXltfqKcqY8ypnyKGeqUSZvnef25taN0js1AMCQIUMwZIhmFkr1J/TLrjzKmfIoZ8qjnKlG03lTOELy8fFR/s1ot29CCCEqUjhCGjVqlCbjIIQQoudUuodECCGEcI31LDs2KioquHw7QggheoT1pAaRSIS9e/fil19+QXNzMyQSibSvvb0dTU1NEAqFqKysVEughBBCdBvrEdKePXtw+PBhNDY2wsTEBI8ePcLIkSNhZGSE2tpaiMVifPHFF+qMlRBCiA5jXZDy8/Ph7u6O4uJiHDp0CACQkJCAgoICHDx4EG1tbbSXHSGEEJWxLkhPnjzB/PnzYWBggDfeeANWVla4ceMGAMDDwwOBgYE4c+aM2gIlhBCi21gXJGNjY5kR0JgxY3Dv3j3p8ZQpU1BdXc1tdIQQQvQG64Lk5OSES5cuSY/t7e2lIySgYwTF4/G4jY4QQojeYF2QQkJCIBAIEBISAqFQiIULF6KiogKxsbE4dOgQjh49ismTJ6szVkIIITqM9bRvf39/CIVCZGRkwMTEBLNnz8bKlStx4sQJAB07O8TGxqotUEIIIbqt1zs11NTUoLGxEQ4ODuDz+VzFRQghRM+wvmQXGhqK8+fPQygUyrSPGjUKTk5Oel2M/v77b8THx8PLywvTpk3D0qVLcfHiRZlz2tra8N1338HDwwNTp05FSEgIbt261UcRa5ekpCSEhobKtTc1NWHnzp3w8vKCq6srPvzwQ8rZvxTlDABOnToFf39/TJ48GX5+fjh27Jhmg9NS3eWsk1AohLe3N62pfIWivAmFQiQlJcHX1xcuLi5YtGgRsrOz0ZsxDuuCVFVVhbi4OMyZMwebNm1CUVERxGKxyh+sK0QiEcLCwnDlyhVEREQgLS0Nzs7O2Lx5My5cuCA9LzExEUePHsW6deuQmpoKQ0NDhIaG6v3MxKysLBw5cqTLvoSEBOTk5CAsLAx79+4Fn8/HmjVrKGfd5CwjIwNff/01/Pz8kJ6eDj8/PyQmJuLkyZMajlK7dJezV+3atQuPHj3SQET9Q3d5i4yMxPnz5/HRRx/hwIED8PT0xPbt25Genq76BzJKKC8vZ3bu3MnMmzePcXR0ZGbMmMHEx8czV69eVeZtdEphYSEzYcIE5ubNmzLta9euZRYvXswwDMNUV1czTk5OTHZ2trS/tbWV8fT0ZBISEjQar7aora1loqKimIkTJzJubm7MmjVrZPpfvnzJODk5Mfv27ZO2CYVCZsqUKTJt+qSnnAmFQsbFxYVJTU2VaY+KimI2bNigwUi1R085e1VpaSkzbdo0xs3NjYmLi9NckFqop7xVVFQwEyZMYH766SeZ9oSEBMbNzU3lz1Vqc9WpU6ciNjYWZWVlOH78OAICAlBSUoLQ0FB4eHhg9+7dqlfGfsrMzAzBwcFyMwzt7e1RVVUFALh69Sra29vh5+cn7efz+fD09JSZSq9PUlNTUVFRgYyMDDg5Ocn1i8ViSCQSmJubS9tMTU0xcOBANDQ0aDJUrdFTzi5fvozm5maEhITItKekpCAtLU1TYWqVnnLWqbGxEfHx8di6dSsGDx6swQi1U095YxgGwcHBmDVrlky7vb09Xrx4gWfPnqn0uSrv9j19+nQkJCTg8OHD8PDwwJMnT5CRkaHq2/Vbs2bNwrZt22TWYInFYpSVlWH8+PEAgPv378PCwgJDhw6Vea2dnR1qamrQ0tKi0Zi1QVhYGPLy8vD222932T9o0CAEBgYiMzMTt27dQmNjI1JSUtDU1IR3331Xw9Fqh55ydvfuXQwZMgSPHz/GihUr4OzsDA8PD72+h9RTzjpt374dDg4OWLFihYYi02495W3SpEnYtm2b3JPDi4qKYG1trfITxVV6hHlFRQXy8/ORn5+P6upqmJqaIjAwEIsWLVIpCG3V1taGnJwchf3Dhw/v8sm6ycnJePDgAfbv3w+g4+bfq3/pdzIzMwPQcfPe2NiYo6j7FtucjRs3rsf3ioyMxPr167F8+XIAHU8k3rFjB1xdXTmLVxtwlbP6+nqIxWKEh4cjLCwMmzZtQmFhIRITE2Fubo6lS5dyHXqf4fL3rLCwEAKBABcuXND5xf1c5u11mZmZuHbtGuLi4lTOI+uCVFlZiYsXL0qLkJGREebOnYvIyEh4e3tj4MCBKgWgzVpbW/HVV18p7Hd3d5cpSAzDIDk5GZmZmVi7di18fX2l7d3Rpf8EyuZMkbq6OgQFBYHP5yMlJQVWVlYoKChAQkICTE1NdWqUxFXOxGIxmpqaEBUVhVWrVgHoGMHX1NRg3759OlWQuMpZfX09EhISEB0dDRsbGw4j1E5c5e11x48fx65du+Dv74/Vq1erHB/rghQYGAgej4fp06cjLCwMCxYs0PlrrWZmZrh79y6rc0UiEWJiYpCXl4e1a9ciOjpa2mdubo6mpia513ROoe9q9NRfKZOz7uTk5KC2thaFhYUYPXo0gI4v1xcvXmD79u3w9/fXmULOVc46R9weHh4y7XPnzkVJSQlevHiBQYMG9fpztAFXOfvqq68wbtw4LFu2DG1tbdJ2hmHQ1tYGIyOVLiJpLa7y1kkikSA5ORlHjhxBQEAAkpKSevX/knW2t2zZgoCAAIwcOVLlD9NVQqEQn3zyCX7//XfExcVhzZo1Mv329vZoaGhAY2MjLCwspO0PHz6Era2tXq/hUqSmpgbW1tbSYtRp+vTpuHDhAurq6jBs2LA+ik472dnZAej44+hVncszdKWAc6mgoAAA4OzsLNN+7tw5nDt3DgKBALa2tn0RmtYTi8XYsmULCgoK8PHHHyM6OrrXv2OsC9K6det69UG6qr29HeHh4bh58yZSU1OxYMECuXNmz54NoOOXPygoCEDHl0ZpaSnmzp2r0Xj7i7FjxyInJwcPHz6UftECQHl5OczNzWUKO+nQ+buUl5eHiIgIaXtJSQkcHR11aiTOlbNnz8q1hYeHY8qUKQgPD8fw4cP7IKr+IS4uDj///DNiY2N7XHDMlsKC5OPjg7i4OOn1RDbXFXk8HoqKijgJrL84deoUrl27huDgYIwYMQLl5eXSPh6Ph6lTp8LGxgaBgYHYsWMHmpubYWdnh4yMDDx//hxhYWF9GL32WrZsGbKysrB+/Xps3LgRVlZWKC4uxg8//IDPP/+cHgbZhTFjxuCDDz7AwYMHYWRkBBcXF+Tl5eHq1avSCTZEVlcbQvP5fFhaWtJm0d0oLS1Fbm4uvL294eLiIvO9B3TMwlPlyo/CgjRq1CiYmprKHBN5nUP+06dP4/Tp0zJ9hoaGqKioAABs27YNgwcPRnp6Opqbm/HWW28hIyND5q9/8p9Bgwbh5MmTSE5Oxo4dOyASiWBvb49vv/0WCxcu7OvwtFZCQgJGjhyJM2fO4MCBAxg7diz27dsnnWBDCBc6v/eKi4tRXFws119WVoYRI0Yo/b693lyVEEII4YLKC2MJIYQQLrGe1MB2brk+rwonhBCiOtYF6a+//pJrk0gkePbsGVpbW2FjYyPdKocQQghRVq/vIbW3t0MgECA+Ph779+/HjBkzuIqNEEKIHuFsUkNycjKuX78uN9OMEEIIYYOzSQ1vvvkm7ty5w9XbEUII0TOcFCSRSITc3FxYWVlx8XaE6L2YmBg4OjrCyckJ9fX1Cs9bsmQJHB0dERMTo8HoCFGPXs+yE4lE+PPPP/H8+XN89tlnnAVGCOmYOFRSUoL3339frq+6upquShCd0qtZdkDHbgT29vYICAiQe1IlIaR3bG1tIRAIuixIRUVFGDp0aLcjKEL6E9YFqavtIQgh6uXj44MzZ86gpaVF7iGOhYWF8Pb27nKDUEL6I9qpgRAt5uvri5cvX+LKlSsy7XV1dbhx4wbmz5/fR5ERwj3Od2p4FY/HQ2ZmptKvI4R0cHNzg6WlJQQCAby9vaXtAoEAJiYmmDVrVh9GRwi3WBckPp+PP/74A7W1tbCwsMDo0aMxcOBAVFdX4+nTp+Dz+fTANEI4ZmhoCC8vL5SUlEAikcDAoOOiRmFhITw9PenhjkSnsC5Iq1atQkREBL788ksEBQXJPI8mPz8fMTEx2Lp1K/z9/dUSKCH6ysfHB+fPn0d5eTlcXV0hFArx66+/Yvfu3X0dGiGcYn0PKSUlBcuXL8fKlSvlHo62YMECrFq1Cnv27OE8QEL03Zw5c2BsbCydWFRWVgYDAwN4eHj0cWSEcIt1QaqqqoK9vb3Cfmtrazx+/JiToAgh/zExMcGcOXMgEAgAdFyumz17NszMzPo4MkK4xbog2dvbIzc3F2KxWK6vpaUF586dg5OTE6fBEUI6+Pr64v79+7h37x4uXbqEd955p69DIoRzrO8hrV+/HpGRkQgMDERQUBBGjx4NhmHw4MEDnDp1CjU1NcjIyFBnrIToLS8vLxgaGiIpKQktLS0yM+4I0RWsC5K/vz9aW1uRnJyMnTt3gsfjAQAYhsGYMWOQnp5Oj54gRE0sLS3h5uaGy5cvY+bMmbC0tOzrkAjhHOuCBADvvfceFi9ejNu3b6OmpgZAxy7fEydOlBYoQoh6+Pj44Nq1a7QYlugszp6HRAghhPSGwkkNGzduxPXr12XaJBIJ7ty5g5cvX8qdn5ubS5MaCCGEqExhQSoqKpKbxt3Y2IjAwECUl5erPTBCCCH6RenNVekKHyGEEHWg3b4JIYRoBSpIhBBCtAIVJEIIIVqBChIhhBCt0O3C2IaGBukCWKBjlh0A1NfXy7QDwLNnz9QQHiGEEH2hcGGsot0XGIbpdleGyspK7qIjhBCiNxSOkAIDAzUZByGEED1HWwcRQgjRCjSpgRBCiFaggkQIIUQrUEEihBCiFaggEUII0QpUkAghhGiF/weh5xvuvxWSNQAAAABJRU5ErkJggg==\n" | |
}, | |
"metadata": {} | |
} | |
] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "fig, ax = plt.subplots(1, 2, sharey=True, squeeze=True)\nax[0].plot(M, z_ZTF)\nax[0].plot(M, zspec)\nax[0].plot(M, zLSST)\nax[1].plot(Planck15.comoving_volume(zLSST)/ Planck15.comoving_volume(z_ZTF), zLSST)\nax[0].set_ylabel('z')\nax[0].set_xlabel('M')\nax[1].set_xlabel('V(LSST)/V(ZTF)')", | |
"execution_count": 26, | |
"outputs": [ | |
{ | |
"output_type": "execute_result", | |
"execution_count": 26, | |
"data": { | |
"text/plain": "<matplotlib.text.Text at 0x1126bec18>" | |
}, | |
"metadata": {} | |
}, | |
{ | |
"output_type": "display_data", | |
"data": { | |
"text/plain": "<matplotlib.figure.Figure at 0x112020b70>", | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ8AAAEiCAYAAAAmtt8/AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlcVWX+wPHPvRcu+yIgyA4q4i6K4l6apqLZQotajdXUtPmrprLMmpnSmhybJrO0daZ9slLLKVNDbdxyX3BBxQVZBWTft3u5vz+OUAyicrncc4Hv+/XixfGce5/z5XCuX57nPIvGZDKZEEIIIaxIq3YAQgghOh9JPkIIIaxOko8QQgirk+QjhBDC6iT5CCGEsDpJPkIIIaxOko8QQgirk+QjhBDC6iT5CCGEsDpJPkIIIaxOko8QQgirk+QjhBDC6iT5CCGEsDpJPkIIIazOTu0A1HTgwAEAdDqdypGIjspoNAIQHR1t1fPKvS3akiXua6n5NMNoNDZcYGF5cn07F/l9W1ZHuJ6duuZT/1dhVFRUk2NJSUkAREZGWjWmzqKzXN+EhARVznu5e1sNneX3bS1qX09L3NdS8xFCCGF1knyEEEJYnSQfIYQQVifJRwghhNVJ8hFCCGF1knyEEEJYnSSfZlQYKjhYeBBjXfvuSy/E/zqRVcKpnFK1wxCdXKce53M5/0z5J9vytqH11DKj9wy1wxHCIiprjEx/ewd6Oy3bnx2Pt6uD2iGJTkpqPs2w19gDsC1zm8qRCGE5Wi1otRoqaozsSylQOxzRiUnyaUYf9z4AHMo5JE1vosNwsNMRFeQJwL6UQpWjEZ2ZJJ9m9HFTkk9pbSlnis6oHI0QljM0rAsA+6XmI1QkyacZvg6+eOm9ADiQc0DlaISwnGFhyn197HwJFTUGlaMRnZUkn2ZoNJqG2s/BCwdVjkYIyxkS2gWNBox1Jg6lFakdjuikJPlcRn3yOZBzAJPJpHI0QliGh5M9kX5uANLpQKhGks9l9HXvC0BeZR7ppekqRyOE5dQ3vUnyEWqR5HMZQU5BeDh4ALA/Z7/K0QhhOcPCleRzMLWIWmOdytGIzkiSz2VoNVqifZVlYvdnS/IRHceIi8mnstbIkYxilaMRnZEknysY2m0oAPty9slzH9Fh+Lo70t3HBYDdyfkqRyM6I0k+VzCs2zAAssuzySzLVDkaISxneHdvQJKPUIcknyuI8IzATX+xZ1D2PpWjEcJyRnRXmt72pxTKcx9hdZJ8rkCn1RHtd/G5j3Q6EB3IiIs1H3nuI9QgyecqDPW7+NxHaj6iA/GT5z5CRZJ8rkL9c5+s8iwySjNUjkYIy6l/7rPrrCQfYV2SfK5CZJdI3PXuAOzN3qtyNEJYzuieSvLZl1JAVa3M3i6sR5LPVdBpdQ21nz1Ze1SORgjLGdXDB4BqQx0HU2WJBWE9knyuUky3GECp+ch4H9FReLno6euv1Op/OZuncjSiM5Hkc5WG+w8HlHnezhWfUzkaISxnTIRS+9lxRp77COuR5HOVunt0x9tRaR/fky1Nb6LjGNVDua+PZhRRXFmrcjSis5Dkc5U0Gk1D05s89xEdSUy4F/Y6DXUm6fUmrEeSTwvUN73tzd6LsU56BomOwVlvx5AQZWnt7adzVY5GdBaSfFpgRMAIAEprSjmef1zlaISwnGt6dQVg2+lc6VAjrEKSTwsEugYS4hYCwO6s3SpHI4TlXHsx+aQXVJKSX6FyNKIzUD35rF27lmnTpjFw4EBiY2NZs2bNZV9fUFDA/PnzGTNmDDExMTz00EOkpKRYJ1hghL9S+5HkIzqSvv7u+LjqAdh2SpreRNtTNfmsW7eOuXPnMmbMGJYvX05MTAzz5s1jw4YNl3y9yWRizpw5bNu2jblz5/Laa6+Rm5vL7NmzKS62zsSI9U1vhy4cotJQaZVzCtHWtFoNYyOU2s9WST7CCuzUPPmSJUuIjY1l/vz5AIwdO5bi4mKWLl3KlClTmrw+JSWFgwcPsnjxYm6++WYAevTowcSJE/n555+55ZZb2jzmmG4xaNBQW1fLwZyDjA4c3ebnFMIarunlw3eHMtl1Np9qgxEHO53aIYkOTLWaT3p6OmlpaUyaNKnR/smTJ5OcnEx6enqT91RXVwPg4uLSsM/DwwOAoqKiNoz2Vx4OHvTz7gfArvO7rHJOIayhvuZTWWtk3zmZake0LdWST3JyMgDh4eGN9oeGhgJw7lzTWQR69+7N8OHDWb58OWfPnqWgoIBXXnkFZ2dnJk6c2PZBXzQyYCQAv5z/xWrnFKKt+bg6MDBI+WPu55MXVI5GdHSqNbuVlpYC4Orq2mh/fa2mrKzsku976aWXeOCBB5g6dSoAer2e5cuXExwcbFYcRqORpKSkJvsrKpQeP5c6FmxUznWm6Aw7j+7EW+9t1rk7s8td347EaDSi06nTfNXcvX05g3x0HMmAn45mMKOX5eLuLL9va1H7elrivlat5nOlsQRabdPQzp49y4wZM+jSpQvLly/nX//6F+PHj+fxxx9n/37rrTLay7UXTjonAA4XHbbaeYVoazHBzgBkltSSUVyjcjSiI1Ot5uPm5gZAeXl5o/31NZ7647/1ySefAPDRRx81POsZPXo0d955J6+++irffvtti+PQ6XRERkY22V//F8WljgGMzBrJz+k/k1yXzCORj7T4vJ3dla5vR5GQkKDauZu7ty8nIsLEgv/mkldWzblqFyZEdrdILJ3l920tal9PS9zXqtV86p/1pKWlNdqfmpra6PhvnT9/nh49ejQkHlDmXIuOjubMmTNtGG1T9b3cdmXtkql2RIeh1WoYH6l0PPhvkjz3EW1HteQTGhpKUFBQkzE98fHxhIWFERAQ0OQ94eHhnD59mpKSkkb7Dx8+TGBgYJvG+79GBYwCoLi6mMT8RKueW4i2NKGPLwB7zxVQUiWzXIu2oeog0zlz5rB27VoWLlzItm3bePHFF1m/fj1PPPEEoMxmkJCQ0NAUd++996LT6fj9739PfHw827Zt4+mnn2bv3r089thjVo09yC2IMPcwAH7JlF5vouMYE9EVvU5LrdHE1iQZcCrahqrJJy4ujgULFrBjxw7mzJnDvn37WLx4cUNPti1btjBjxgwSE5WaRVBQECtWrKBr167Mnz+fp556iqysLD7++OOG91jTmMAxAGzP3G71cwvRVlwd7BjVU+nBufF4jsrRiI5K1RkOAGbOnMnMmTMveSwuLo64uLhG+3r06MG7775rjdCuaGzgWL448QXH8o5RUFWAl6OX2iEJYRGT+nZjS1Iu/z15gRpDHXo71aeBFB2M3FGtEN0tGic7J0yYpOlNdCgTLz73Ka02sOecLDAnLE+STys46BwY3k1ZYG5H5g6VoxHCcnzdHRkc4glAfKI0vQnLk+TTSmODxgKw8/xO6XItOpTr+/oBynOfujpZYE5YliSfVhobqCSfouoijuYdVTkaISxncr9uAGSXVHE4wzoT94rOQ5JPK/m7+hPRJQKALelb1A1GCAvq0dWVSD9lppH1x7JVjkZ0NJJ8LODaoGsB2JqxVeVIhLCsKf2V2s+6o1lXnI9RiJaQ5GMB9cnnTNEZMssyVY5GCMuZOsAfgIzCSo5lllzh1UJcPUk+FjDAZwBdHLoAsDVdaj+i4+jl50r3rsoyJ+uPZakcjehIJPlYgE6ra+j1Jk1voiPRaDRM7a/UfqTpTViSJB8LqW9625u9l7KaSy+EJ0R7VN/0lpJfQeJ5aXoTliHJx0JGB47GXmuPoc7AjvMy4FR0HH383ehxsent+8PnVY5GdBSSfCzExd6F4f7KbAf/TfuvytEIYTkajYbpg5QlTtYePi8DToVFSPKxoPHB4wFlluvaOlkHRXQcNwxUks/54ioOphWqHI3oCCT5WNC44HEAlNaUciDngLrBCGFBPX1d6evvDsAP0vQmLECSjwX5OvsywGcAAJtTN6scjRCWdWPUxaa3I1nUGutUjka0d5J8LOy6kOsA+Dn9Z+pM8gEVHceNF5/75JfXsON0nsrRiPZOko+FTQyZCMCFigscyzumcjRCWE6ApxMjuisLJn57SGbyEK0jycfCwjzC6OnZE4BNaZtUjkYIy4obHARAfGI2pVXSqUaYT5JPG5gQMgFQnvvIiHDRkcQO6IaDnZZqQ53MdC1aRZJPG5gYqjS9pZWmcarwlMrRCGE5bo72DYvMfXswQ+VoRHsmyacNRHaJJMhVaZ7YmLpR5WiEsKxbhyj39u7kAtILKlSORrRXknzagEaj4fqw6wGIT42XpjfRoYyN8MHP3QGAlQek9iPMI8mnjUwOnQzAueJznCk6o3I0QliOnU7bUPtZfSBDptsRZpHk00b6evcl0DUQUGo/QnQkt0UrySezqJKdZ/NVjka0R5J8muF+9nvCfpwBWUfMer9Go2FS6CQA4lMk+YiOpXtXV4aFKQsofr0/XeVoRHskyacZLtl7cCg5BzveMLuMSWFK8kkuTuZ04WlLhSaETbh9aDAAPx3LprC8RuVoRHsjyacZ5f4jlY2kDVBt3uJw/bz7NfR6W39uvaVCE8Im3DDQHzcHO2qMdayWbteihST5NKMsaBx1WnswVMKpDWaVodFoiA2PBWBDygbp9SY6FGe9HTcPVp5rfrk3Te5v0SKSfJpRp3elPGCU8o9jq80upz75pJemczz/uCVCE8Jm3Dk8BIDk3HL2nCtQORrRnqiefNauXcu0adMYOHAgsbGxrFmz5rKvr6ur491332XChAkMHDiQ6dOn8+OPP7ZJbKUhyjMbTm+ESvMW0IroEtEw19u6c+ssFZoQNqGPvzuDQzwB+HJPmsrRiPZE1eSzbt065s6dy5gxY1i+fDkxMTHMmzePDRuab+Z69dVXeeedd7j77rt5//33GTRoEE8//TRbt261eHxlgWPA3hnqauH492aXMyVsCgAbzm3AWGe0VHhC2IQ7Y5Taz4Zj2eSWVqscjWgvVE0+S5YsITY2lvnz5zN27FgWLFhAbGwsS5cuveTr09LS+Pe//81f/vIX7rvvPkaOHMkrr7zC0KFD2b59u8XjM9k5Qe9pyj+OrjS7nKnhUwG4UHmB/Tn7LRGaEDZj+qAAPJ3tqTHW8fU+qf2Iq6Na8klPTyctLY1JkyY12j958mSSk5NJT286dmDTpk04Ojpy8803N9r/xRdf8Kc//altAh1wh/I9ZQcUm7eGSbB7MAO7DgRgbfJaS0UmhE1wtNcx42K36y92p2GQVU7FVVAt+SQnJwMQHh7eaH9oaCgA586da/KepKQkwsPD2blzJzfeeCN9+/Zl0qRJrFvXhs9SeowHZx/ABMdWmV3MDd1vAGBT6iaqDFUWCk4I23D3iFA0GsguqWLj8Ry1wxHtgJ1aJy4tLQXA1dW10X4XFxcAysqajq0pKCggKyuL559/nieeeIKgoCBWrlzJk08+iZeXFyNGjGhxHEajkaSkpCb7KyqU2XqTziTjGzieLqdXUrXvM1J9prT4HAA9DD3QoqWstowVe1cwynuUWeV0FA3X9xLXviMxGo3odDrVzm3N6xsT5Mye9Are3Xyc7vqSRsc6y+/bWtS+npa4r1Wr+VxpTIBW2zS02tpaCgoKeOWVV7jjjjsYNWoU//jHP+jduzfLli1rq1ApudhhwLHoDA5mrs/jYe9BlGcUANvytlksNiFsxY19PAA4kl3F2XzpeCAuT7Waj5ubGwDl5eWN9tfXeOqP/5aLiws6nY7Ro0c37NNqtYwaNYpVq8xrEtPpdERGRjbZX/8XRWRkJPTqBQdehYKzhBXvhhHTzTrXLP0sDm47SEJRAj4hPng7eZtVTkfQ6Pp2YAkJCaqdu7l7u6306mXi44RSzlwoY3OGiamjfj13Z/l9W4va19MS97VqNZ/6Zz1paY17x6SmpjY6/luhoaHU1dVhMBga7a+trUWj0bRRpIBGA1GzlO0j34DRcPnXN2N8yHjc7N0wmAwy3Y7ocDQaDfePUT63Pxw+z4VSebYpmqda8gkNDSUoKKjJmJ74+HjCwsIICAho8p6xY8diMplYv/7X/7gNBgPbt28nOjq6bQMeOEP5Xn4Bzm42qwgHnQOTw5V1fr4/a/64ISFs1S2DA/Fy0VNjrOOLXalqhyNsmKrjfObMmcPatWtZuHAh27Zt48UXX2T9+vU88cQTgNLBICEhoaEpbuTIkVx77bW88sorfP7552zfvp3HHnuMzMxM5syZ07bBeoZA2FhlO+FLs4u5qcdNAJwoOMEpM58fCWGrHO113HVxyp0v9qRRWSODqsWlqZp84uLiWLBgATt27GDOnDns27ePxYsXM3WqMihzy5YtzJgxg8TExIb3vPXWW8ycOZMPPviAOXPmUFhYyEcffUT//v3bPuCoO5XvSeugwrx5rAZ1HUSou9KdfM2Zy08lJER79LuRoeh1WgrKa1h5QNb6EZem+txuM2fOJD4+nqNHj7Ju3bpGA0jj4uJISkpi+PDhDfscHR2ZN28e27dv58iRI3z11VfExMRYJ9i+N4HeDYw1Zs94oNFouLmn8jOuPbuWWmOtJSMUQnW+bo7EDVFmu/5we7IMOhWXpHryaVf0LtA/Ttk+9LnZxUzvPh2tRkthdSFbMyw/J50Qanvwmu5oNJBeUMm6Y9lqhyNskCSflhr8O+V79lHIOmxWEX4ufowOULqLf3fmO0tFJoTN6N7Vlcl9uwHw3pazstaPaEKST0sFDQWfi33rD3xqdjG3RNwCwI7MHWSXy1+GouN5eFwPAI5nlbA/s0LlaIStkeTTUhoNDJmtbB9dCTXll399M8YFjcPL0Ys6U53UfkSHFBXsyageykDqLw8XSu1HNCLJxxyDZoFOD9UlkGhejzV7nT039VS6XX93+jtZ50d0SI9dFwHA8QtVHM6uVDkaYUsk+ZjDxRv63KhsH/jE7GJujbgVgKzyLHae32mBwISwLSO6ezE0tAsAKw6btxqw6Jgk+Zgr+h7le8ZeyEm8/GubEeoeSkw3pZv4qlPmL9cghK3SaDQ8NkGp/SRkVbI/xbzxcaLjkeRjrrCx4KU8UGX/R2YXc3uv2wHYkrFFOh6IDumaCB96+TgAsGSTzOohFJJ8zKXRwLD7le3DX0N1qVnFTAiZ0NDxYPXp1RYMUAjboNFomD3YC4BfzuSzOzlf5YiELZDk0xqDZoGdI9SUKrNdm8FeZ9/w7Gf1qdXU1smMB6LjGRroTJ+ujgC8sfGU9HwTknxaxdkL+iuJg/0fgZkfqNt63YYGDbmVufyc9rMFAxTCNmg0GmYPUWo/e88VsPOs1H46O0k+rVXf9JZzDNJ2mVVEgGsA1wZdC8CKkyssFZkQNmWwvxMxYUoCeu2nJKn9dHKSfForMFr5AtjzntnFzOqtLFZ3IOcASQWyzr3oeDQaDc9MUWYHOZxexE+JOSpHJNQkyccShj+sfD+xFoozzCpiZMBIwj2UVSCl9iM6qmFhXkzo7QvA3386KTNed2KSfCyh783g4gsmo9ndrjUaTUPtZ23yWoqqiiwZoRA2Y+7kSDQaOJtbzuqD5v2xJto/ST6WYKeHofcp2/s/hlrzphG5sceNuNq7Um2sZtVpGXQqOqY+/u7cHKWs9/PGxlNU1BhUjkiooUXJZ82aNaSnN78yYUJCAvPnz291UO3S0PuV+d4qC+DwV2YV4WLv0jDb9YqTK6Tbteiwnp7UC72dlpySav65/Zza4QgVtCj5PPfcc9x2221s3XrpBdDS0tJYs6aTLg3t5gcDlNkK2P2u2d2u7+pzF1qNlgsVF4hPibdggELYjqAuzvx+tPKM872tZ7lQWqVyRMLazGp2e+SRR1i2bJmlY2n/RjyifM9LgjObzCoi0DWQCSETAPjs+GfSHVV0WI+O74GXi56KGiNvxMu0O51Ni5PPc889x4wZM1i2bBkPPfQQJSUlbRFX+9RtAIRfo2zvfMvsYmb3VdYLOp5/nP05+y0RmRA2x93Rnj9OVCYd/Xp/Oscyi1WOSFhTi5OPnZ0dL774Ii+//DI7d+7k1ltv5eTJk4DSY6vTG/WE8v3cNjifYFYRUb5RRHWNAuDjYx9bKjIhbM6dMSH08nPFZIIFPyRKTb8TMbu32+23386nn35KZWUlM2fO5Pvvv8fBwcGSsbVPPSeAbz9luxW1n3v73wvA9sztnC48bYHAhLA9djotf7lB+bzsSynkhyNZKkckrKVVXa2HDBnC6tWr6dGjB/PmzePzzz+3VFztl0YDox9XthO/g8IUs4oZHzyeMPcwAD5J/MQioQlhi8ZE+DCprx8Ai9adkK7XnUSrx/n4+fmxYsUKpk+fzr59+ywRU/vX/1ZwDwJTHfxiXu1Hq9Fyb797AViXvI7zZectGKAQtuVP0/qit9OSVVzFW5vPqB2OsIIWJZ/PPvuMUaNGNdmv1+t57bXX+POf/8zNN99sseDaLZ09jHpM2T70BZSaN4fV9B7T8XXyxWAySO1HdGgh3s48cq2yOOM/tydz5oJ562OJ9qNFyScmJgZvb+9mj991110sWrSo1UF1CENmg7M3GKth93KzitDr9Mzup/R8+/b0t+RV5lkyQiFsyiPjehDi5YyhzsSf10jng45OptdpK3rnX8f97PsXVJi3dv3tvW7Hw8GDamM1nx+XZ2qi43K017HgJqXzwa7kfNYkZKockWhLknza0rA/gIM71JQpsx6Ywdnembv73A3AVye/kglHRYc2PtKXKf26AfDy2hMUlteoHJFoK5J82pKT56/LLex5DyrNSxx39rkTN3s3KgwVfH5Caj+iY3vpxn64OdhRUF7Dq+tOqB2OaCOSfNraiEdA7wrVJWYvNueud+euvncB8OWJLymulpHgouPq5uHIs7G9AVh5IIOdZ+RZZ0ekevJZu3Yt06ZNY+DAgcTGxrZoYtKsrCyio6N555132jDCVnL2gpgHle1d70BloVnF3N3nblzsXSirLeOz459ZMEAhbM9dMSEMCfEE4Llvj8rYnw5I1eSzbt065s6dy5gxY1i+fDkxMTHMmzePDRs2XPG9JpOJ559/nrKyMitE2kqjHgO9G1QXKwnIDB4OHg3Pfr44/oU8+xEdmlarYfGtA9HrtKQVVPD3n2Rp+Y5G1eSzZMkSYmNjmT9/PmPHjmXBggXExsaydOnSK773yy+/JDk52QpRWoCz168933a/a3bPt9n9Zjc8+/k4UeZ8Ex1bhJ8bT1ycePSTnSnsTzHvcyNsk2rJJz09nbS0NCZNmtRo/+TJk0lOTr7sonXp6em8/vrrvPzyy20dpuWMfBQcPKCmFHYsMasId717w7ifFSdXyLgf0eE9eE13+ge6YzLBM6uOSPNbB2Kn1onray3h4eGN9oeGhgJw7tw5goODm7yvrq6O5557jtjYWK655ppWx2E0GklKalqlr6ioALjkMXN59ZpJ16PvU7fnfc51nYTB2bfFZQy3G85ndp9Raihl8dbFPBD+gMXis6a2uL62yGg0otPpVDu3rVzf1vy+5wz14LGsEs7llfPcit3MGdHV0uG1O2p/fixxX6tW8yktVabPcHV1bbTfxcUFoNlnOZ9++ikZGRntcrnuwshZGBy80Bqr8U78l1llONs5ExcYB8CmC5vIqTJv6h4h2ovuXg7MHqzMrPL9iWIOnq9QOSJhCarVfK40dYZW2zQvnj17ljfffJO33noLNzc3i8Sh0+mIjIxssr/+L4pLHWuV8vmw/hk8z36P55Q/gU/PFhcR1jOMn/J+Irs8m3Ul6/jboL9ZNkYraLPra2MSEsxb08kSmru31dDa3/cLESYScndyMK2It3YXsOGJ/ng421syxHZF7c+PJe5r1Wo+9cmjvLy80f76Gs//Jhej0cj8+fOZMmUKo0ePxmAwYDAo7b91dXUN2zYv+l7wDAWTETa/ZFYRDjoHHh30KAA/Jv/I8fzjlotPCBuk02p4444onOx1ZBVX8fx3R2Xut3ZOteRT/6wnLS2t0f7U1NRGx+tlZWVx+PBh1qxZQ79+/Rq+AN5+++2GbZtnp4cJf1G2T/wAaXvMKubGHjfS01OpNb2x/w35IIoOL8zHhRen9wXgx6NZrDqQoXJEojVUSz6hoaEEBQU1GdMTHx9PWFgYAQEBjfb7+vqyatWqJl8As2bNathuF/rFQcBgZTv+T2BG4tBpdTwV/RQAe7L3sCNzhyUjFMImzRgW3DD324vfJ5KSV36Fdwhbpeo4nzlz5rB27VoWLlzItm3bePHFF1m/fj1PPPEEAAUFBSQkJFBWVoZer2fAgAFNvkBJTPXb7YJWC5NeUbYz9iornpphTOAYhvsPB+D1/a9jqGsnTY9CmEmj0bAobgDd3B2pqDHyfysOUm0wqh2WMIOqyScuLo4FCxawY8cO5syZw759+1i8eDFTp04FYMuWLcyYMYPExEQ1w2wbYWOg9w3K9sYXobayxUVoNBqeGfoMGjQkFyez8tRKCwcphO3p4qLnzZlRaDVwLLOERetOqh2SMIPqc7vNnDmT+Ph4jh49yrp16xqthBoXF0dSUhLDhw9v9v1JSUk8+uij1gjV8q5fCFp7KE6DXcvMKiLSK5K4CKXr9fKE5TLpqOgURnT35vEJv85+8FNitsoRiZZSPfl0at49YMTFJRe2vwHF5i2e9X+D/w8XexeKq4tZnmDeqqlCtDePXRfBiO5eAMxdeZjUfHn+055I8lHbNc+Ciy/UVsDGP5tVhI+TD48MUuaO+zrpa5IKbGNUuxBtSafV8NbMwfi4OlBaZeDhLw5SVSvPf9oLST5qc3SHiS8p28dWQ4p5vdbu7H0nYe5h1JnqWLR3kXS9Fp2Cr7sjy+4cjE6r4URWCX9ec0zu/XZCko8tGDQLgoYp2z8+DYaWLx1sr7Nnfowy5dCBnAOsTV5ryQiFsFkjunvzzGRlpP/KAxn8e0/aFd4hbIEkH1ug1cK0f4BGC7knYbd5z21GBY7i+tDrAaXrtXQ+EJ3FQ9d0J7a/Mv5nwQ+JsvxCOyDJx1b4D/p1xdMti6Ew1axi5g2bh7OdMwVVBSw9eOV1kYToCDQaDX+/fRARvq7UGk088u+DZBW3fPiCsB5JPrZk/Avg2g0MlUrzmxlt134ufsyJmgPAylMZzdljAAAgAElEQVQrOXThkKWjFMImuTrY8cHsobg52pFbWs2Dnx2gskY6INgqST62xNEdpr6mbJ/ZqHRAMMOdfe6kr7cyB9ZLO1+ixtjyZ0hCtEfhPi68PWswWg0czSzmmVWHpQOCjZLkY2v63AiR05Tt9fOgPL/FRdhp7Xhp5EvoNDqSi5P559F/WjhIIWzXuEhfXpim/PG19kgWb246rXJE4lIk+dgajQam/h0c3KEiDzbMM6uYPt59mN1XWXL7wyMfytgf0an8fnQYM4cpKyEv3Xya7w7JDNi2RpKPLfIIhEkvK9tHV8LJdWYV82jUo4S5h2EwGfjzL3+WiUdFp6HRaHj55v6M6qGsgPrsqiPsSW55K4JoO5J8bNWQeyD8WmV77R+houVdRx3tHFk4eiEaNJwoOCHNb6JTsddpeffuaHpe7AH3h8/2czqnVO2wxEWSfGyVRgM3vg16NyjLgR+fMquYwb6DuavPXQC8f/h9WfVUdCoeTvZ8fO8wfFwdKKkycO/H+8gpqVI7LIEkH9vWJRSmvKpsJ34HR81bMO+JIU8Q7hGOwWTghR0vUG2stmCQQti2YC9nPr53GM56HZlFldz78T5KqmrVDqvTk+Rj6wb/DnpNUbbXPgVF6S0uwtHOkVfHvIpOo+NM0RnePPCmhYMUwrYNCPLgnbuGYHdxDrgHPt0vk5CqTJKPratvfnP2gepi+O5hqGv5h6a/T38eHqQs3/DFiS/4JfMXS0cqhE0bF+nL328fCMDecwU8tuIQBmOdylF1XpJ82gNXX7jp4nxvqTtgxxtmFfPAgAeI6hoFwAs7XiCvMs9SEQrRLtwyOIg/36CMAdp4PIdnVx2hrk4GoapBkk97ETkFhv1B2f7vq5C6q8VF2Gnt+Ns1f8PN3o38qnye3/48dSb5y090LvePCef/xvcE4NtDmbz0Q6LMgqACST7tyaRXwG8AmOpg9f1mzX4Q6BrIi6NeBGBX1i7+dfRflo5SCJv39KRe3DsqDIDPdqXytw0nJQFZmSSf9sTeEW7/GOxdoCQTvv0D1LW85jI5bDJ39LoDgGUJy9ibtdfSkQph0zQaDX+5oS+3RQcB8P7WZJZsPKVyVJ2LJJ/2xicCpl9cKuHsZtj2d7OKeTbmWfp49aHOVMez257lQsUFCwYphO3TajUsvnUgN0UFAPDWz2dYKvPAWY0kn/Zo4O0w7AFle8siOBXf4iIcdA78Y9w/Gp7/PL3laWqNMvZBdC46rYZ/3D6IqQOUheiWbDrFGxtPSROcFUjyaa8mv3px6W0TrH4A8s+2uIhgt2AWjV0EQEJuAov3LbZwkELYPjudlqUzBzckoLc2n+b1+CRJQG1Mkk97ZecAd3wGLl2V8T9f3QlVJS0u5trga3lk0CMAfJ30NStPrbR0pELYPPuLCWjaQH8Alv/3LIvWSyeEtiTJpxnnS2r5z/EiSm15Gg73ALjjc9DaQ+7Jix0QWj4A9eFBDzM+eDwAr+5+lf3Z+y0dqRA2z16nZemMqIZnQB9sS+al7xNlHFAbkeTTjE8P5vPOnjz+8Nl+agw2PBYmdCTccHHQ6akNsPEvLS5Cq9GyaOwienr2xGAy8OSWJ0krSbNwoELYPjudljfuiOL2i73gPt2VyjOrjshMCG1Akk8zrgl3BWB3cgHPf3fUtqvfQ2bDiDnK9q5lsP+jFhfhYu/C29e9jZejF0XVRczZPIfi6mILByqE7dNd7AX3uxGhAKw+mMGj/z4oc8FZmCSfZowOdWX2YC8AVh3I4K3NZ1SO6AomvQy9YpXtH+fCqZ9aXESQWxBLxy9Fr9WTUpLC4z8/LjNgi05Jq9Ww8KZ+DTMhxB/P4Z6P9sps2BYkyecy7hzUhVuHKNXvJZtO8c2+ls8obTVaHdz6T/AfBCYjrLwXMg+0uJgo3yj+OuavABy8cFCm4BGdlkajYe7kSP40rQ8Ae84VMOP93VyQ9YAsQvXks3btWqZNm8bAgQOJjY1lzZo1l319bm4uf/rTnxg/fjyDBw8mLi6O9evXt0lsGo2Gv906gLERPgDM/+4oG4/ntMm5LMLBFe5cCZ4hUFsB/74d8lo+aG5K+BTmDp0LQHxqPIv2LLLtZkch2tADY7vzxh2DGpZjuOWdnZy5UKZ2WO2eqsln3bp1zJ07lzFjxrB8+XJiYmKYN28eGzZsuOTra2pqeOCBB9i5cyePP/44y5Yto3///vzxj39k7dq1bRJj/VK8/QPdMdaZmPPlQXbb8lrwbn5w97fg7A0V+fDZzVCc0eJiZvedzd197gbgq6SveOfwO5aOVIh2I25IEB/eMxQne2VButve28n+lJYvbS9+pWryWbJkCbGxscyfP5+xY8eyYMECYmNjWbp06SVfv23bNk6ePMmbb77JLbfcwujRo1m4cCFjx47lww8/bLM4XR3s+OS+GLr7uFBjqOP+T/aRkF7UZudrNZ8IuGsV6F2hJAM+uwnKWjZ9jkaj4Zlhz3BD9xsAeO/we3ya+GlbRCtEuzA+0pevHxqBj6ueoopa7vznHr4/fF7tsNot1ZJPeno6aWlpTJo0qdH+yZMnk5ycTHp60+crLi4uzJgxgwEDBjTa3717d9LS2rZrsI+rA5/dH0OAhyPlNUZm/2sPiedtuDdY4BCY+SXoHCD/jFIDqmjZX2pajZaFoxcyLmgcAK/vf52vTn7VBsEK0T4MDPJk9SOj6N5V+UP08RWHWPbzaWmWNoNqySc5ORmA8PDwRvtDQ5XujefOnWvynpEjR7Jw4UI0Gk3DvtraWrZu3UpEREQbRqsI6uLMFw8Mx8fVgZIqA3f9cw/Hz7d8VgGr6X4tzLg4CPVCInx6Y4sTkL3WntfHvc5I/5EA/HXPX/km6Zu2iFaIdiHU24VvHxnF8HClN+zr8ad46pvD0hW7hezUOnFpaSkArq6ujfa7uLgAUFZ2dQ/0/v73v5OSksLy5cvNisNoNJKUlNRkf0VFBcAlj/11oi/PbjhPUUUtM9//hb9NDqSHt4NZ5297YbiOeoWAX55Hk3OUqg8mkXHdcowOni0qZU7QHErKS0gsSeTl3S+TnZPNZL/JZkd1uevbkRiNRnQ6nWrntpXr2xF/338e24Wlulo2ninlu0OZnMzI5y/XdcPLue3/W1X7elrivlat5nOlaqpWe/nQTCYTr732Gp9++in3338/EydOtGR4lxXWxYHXpgTg4aijpLqOZzdkkpRru90vy4LHc370XzFpdDgWnSZ488PoKlvWacJR58j8yPn0dVOWIP7w3If8mPVjW4QrRLtgr9Pw9Bhf7h/qjQY4kVvFYz+kcyrPdv8vsCWq1Xzc3NwAKC8vb7S/vsZTf/xSampqeO655/jxxx+5//77efbZZ82OQ6fTERkZ2WR//V8UlzoGEInSZHjnh3vIK6vm+Y3Z/POeoYzo7m12LG0qMhICg2DlfTgUJ9Nz2//B7DVKt+wW+DjiYx7/+XH2ZO/h49SPceniwoMDH2zUFHo1rnR9O4qEhATVzt3cva2Gjvz7/nNvGNk3hz9+nUBehYGn15/nrzf35/ahwW12TrWvpyXua9VqPvXPev63o0Bqamqj4/+rrKyM++67j/Xr1/P888+3KvG0Vi8/N755aAT+Ho6UVRu456O9bLLlcUB9piudEOwcoeAsfDQFcltWbXe2d2bZhGWMCRwDKCuh/n3/32UgqujUJvb147tHRxHm7UyNoY5nVh3hT2uO2va8kCpTLfmEhoYSFBTUZExPfHw8YWFhBAQENHmP0WjkkUce4fDhwyxZsoR77rnHWuE2q3tXV1Y+PJJwHxeqDXU89MUBvt5nw5Ny9poEd68GvZuyFPe/JkHanhYV4WjnyFvj32JK2BQAPj/+Oc/veF4WoxOdWoSfG//5vzFM6O0LwBe705jxwS6yiitVjsw26V566aWX1Dq5m5sb7777LoWFhWg0Gj766CPWrFnDiy++SEREBAUFBSQlJeHq6oper2fFihV8+eWX3H777QwbNozs7OyGr5ycHLp169ai82dnZwNc8n35+cozER8fnyuW4+5kz7QB/uxKziOnpJpNJy5gMsGI7l4tbo6yCs8Q6HEdnPwRKgvg6ErwiYSuV1+F12l1TAiZQGF1IYn5iZwuPE1CbgLjgsfhoLty54uWXN/27HL3WEc8b3M6y+/b0V7H9IEBaDTKdDxZxVWsOZRJvwAPQrydLXYeta+nJe4vVQeZxsXFsWDBAnbs2MGcOXPYt28fixcvZurUqQBs2bKFGTNmkJiYCMBPPymTZX799dfMmDGj0desWbNU+zkAuro58NWDIxum4lm6+TRPf3OYaoONdr8MiIL746FLOBiq4JvZsPNtaMF4BZ1WxwvDX+CxwY8BsCdrD7PXzSazLLOtohbC5mm1Gv44sRcf3zsMDyd78str+N1He3h782mMsjZQA42pE4+Oqn9oFhUV1eSYuQ/0ao11vPDdUb7Zr0xpEx3ahfd/F42Pq412xS7PU1ZBTb/Y9Db4bpi2BOz0LSrmh7M/8Jedf8FQZ8DL0Ysl45YwxG9Is69X+4GptVzuHuuI521OZ/l9/6/0ggrmfHmQIxnKgPTRPb1ZckcUvu6OrSpX7etpiftL9YlFOxp7nZbFtw7kmcnKTXEgtZCblv3CsUwbnQ3BxQdmfw/9b1X+fegL+HQ6lLas48T0HtP54PoP8HDwoKCqgPvj72fVqVVtELAQ7UewlzMrHx7J7JHK4PlfzuQTu3Q7W5JaNt1VRyTJpw1oNBrmjO/J+7+LxlmvTER467s7WX2g5RN8WoW9I9z6Lxj/gvLv9N3wwbWQvrdFxQzrNowvp35Jd4/uGOoMLNi1gAW7FsiaQKJTc7DTsfCm/rx3dzTujnbkl9dw78f7eHXdCdttlrcCST5taHK/bnz36GhCvZ2pNtTx9MrDzP/2qG1Ow6HRwLXPwswVSk+40iz4OBZ2v9ei50Ah7iH8e+q/uS74OgBWnVrF7PWzySi10cQrhJVM6d+NdU+MJTq0CwAfbEvm5uU7OZVTqnJk6pDk08Yiu7nx/Zxfu1+u2JvGLe/s5Gyuja4H0nsqPPhf6Nob6gywYR58fTdUFl51Ea56V5aMX8Ljgx9Hq9FyPP84d/xwBxtTN7Zh4ELYvqAuznz94Agev64nWg2cyCph+ts7+GxXSqebnFSSjxV4ONvz4eyhPDslEt3FBammv72Dr/el2eYN5xMBf/gZBs5Q/n1yLbw3FlJ+ueoitBotfxj4B96//n28HL0orS3lqS1PsXDXQipqK9oocCFsn51Oy1OTIln58EiCvZyoNtTxl/8kcs/H+zrVmCBJPlai1Wp4dFxPvnpwBAEejlTUGJm3+igPf3GA/DIbfCaid4Fb3oebloO9MxSnwyfTYOOLYLj6eEf4j2D1jasZ4T8CgJWnVjJj7QxOl7V8hVUhOpLoUC/WPT6WuMGBAGw7lcukJdtYuT/dNv8otTBJPlY2LMyLdU+MZdoAfwB+Ssxh0pJtbDiWpXJkl6DRKF2vH9wK/oMAE/zyJnwwHrIOX3UxPk4+vH/9+zwV/RR2WjtSSlJ44dgLrEhfQY2xpu3iF8LGuTna88aMKN67ewjeLnpKqww8s+oI93+6n5ySjj1BqSQfFXg661l252D+cfsg3C72fnn4i4PM+fdBckttsBbUtRfcvwnGzgWNVlkb6IPxsHkh1F7dB0Sr0XJf//v4cuqX9PTsSR11rM5czR0/3MGR3CNt/AMIYdum9Pcn/slrmDpAmTHg55MXuP6NrXx7MKPD1oIk+ahEo9Fwa3QQ8U9ewzW9ugLw49EsJvxjCyv2plFnayOh7fQw4c9KEuraG0xG2P4PeHcUJG+56mL6ePfh6xu+5paAW9Ci5WzxWe5edzd/3f1XSms6Z68fIQC8XR14565o3p41mC7O9pRUGXjqm8P8/pN9nC/qeM+CVJ3bTW2WmtutNdwc7bk5KoBgL2f2phRQXGlg84kLbD+dS78Aj1aPhLY49wAYMlvpfp2xFyry4fBXkHcGgmLAofmlMOrptDoCDAEM8RxCem06eVV5HMs/xn/O/gcfJx8iPCNsc048M8jcbgq15yJrTyK7uXHrkCDSCso5m1tOSn4FX+1Nw1mvY2CQJ1qNRvXraYn7S5IP6iYfUGpBfQPcuS06iNzSak5ml5JVXMWKfWnklFQzONgTZ71qSy81pbVTlujuPR1yjimzY184Dgc+ATQQMBh0l483Pz8fL70XD458EFd7VxJyEyipKWFz2mZ2Z+0m0iuSrs5drfLjtCVJPgq1/7Nsb1wc7LhhoD+R3dzYcy6f4koDW0/lsiXpAgMCPdDVKEM1JPm0U7aSfOo56+2Y0r8bw8O9OZpZRH5ZDccyi1mxJw2NRsOAQA/sdDbUUuraFaLuAs9QZW64qiI4txWOfgMuvkrzXDM1mPrr69vVlyjfKKb3mE52eTbJxclkl2ez+tRqMsoy6OfdD1e96yXLaA8k+Sgk+bScRqMhws+NGUNDKK6s4VhmCTkl1Xy1L43SKiN9fZ3w91PnDzRJPq1ka8mnXrCXM7NiQvBy0XMovYiSKgO/nMlj9cEM3Bzt6N3NDa3WRpqlNBrwHwhD7gFjLZw/pAxIPfE9nPoJuoRCl7AmSeh/r6+r3pXJYZOJ9ovmRMEJ8qvySSpM4uuTX1NWW0Zfr7442tlYE+RVkOSjkORjPkd7HRP7+DG6pw8JaUXkl9dwMq+ajWdK8O/iTC8/V6s3U0vyaSVbTT6gjAuKCunCrJhgagwmjp0vprjSwKYTF/jhSBaeTvZE+LmhtZVnI/aO0HMC9L9NaYbLOwVl2XDka0jZriSh3yzZ3dz1DXIL4raI2/Bz8SMxP5HS2lISchP45tQ3VBmqiPSKbFdJSJKPQu3PU0cQ6OnEzGEhOOvt2J+ST1lNHeuPZbMvpYBBwZ54ubRsJvrWkOTTSracfOo52uu4NrIrtwwOpKiilqScUgoratmQmM3aw+dx1uvo5eeGzlZqQs5eygzZPSZAwTkoTlMGqCZ8Cee2g7s/dAm77PXVarT09e7LHb3uwMnOiRP5JyirLeNAzgG+SfqGitoKIrpE4GTnZO2frsUk+Shs5fPU3um0GoaFeTHIo5qcMgPpxbWkF1ayYm8a5TUGBod0QW/X9k3zknxaqT0kn3oeTvZM6d+NaQP8KSiv4UxuGYUVtWw8nsO3B5XF23r5uVnlxrsqHoEQdSeEjICCZKU2VJyu1IRObaBG40iNeyg+XX2bLcJeZ0+0XzS39boNnUZHUmES5bXlHLxwkK9OfkVeZR7hHuG4O7hb8QdrGUk+Clv7PLV3VaWFjOvuxnVR3TmUVkRhRS0HUgtZfTADL2c9vbu5tWlTnCSfVmpPyaeet6sD0wb6M3WAPyVVtZzOKaWkysC2U7l8vjuVwvIawn1c8HCyVztU5TmPVzgM/p2ShEoyoSgNyrJxS/8Zj5R16DQmZQlv++ab0hztHBnuP5zbe92OXqcnqTCJCkMFx/KOseLkCk4XnsbHyQd/F3+b66ItyUdhq5+n9qr+eg7rHcqsmBDsdVoOpRVSXGkg/ngOm07kEOjpRKi3c5t8JiT5tFJ7TD71vF0diO3vz82DAzHWmTiVU0pFjZGDaUV8sjOFY5nFuDnZE+LlrP5zofokFHUndL8OKgsx5Z1GV1sKZ3+GPe9DUSq4+StfzXC0cySmWwwze8/E08GTs0VnKastI7k4mTVn1rA5bTMaNIR5hKHXWa/9+3Ik+Shs/fPU3vz2etrptIzo7s0tgwO5UFrNqZwyckurWZNwnkPpRQwK8qCLhZ8HSfJppfacfOp5Ouu5rrcvdw0Pxd3JnuTccsqqDSTnlfOfhPOsPpBBWZWBYC9n3B1toDbkEQT9b+Wc61BAg1NpKtRWKHPFHfgEkjYosyd4dW+2NqTX6YnyjWJWn1n08OxBTnkOORU55Fflsy1jG1+e+JLMskw8HTzxc/ZTtTYkyUfRXj5P7cWlrqe7kz1TB/gzLtKXtIIK0gsrSc2vYMXedMpqDAwI9MDRXmeR80vyaaWOkHzqOel1xIR7cc+oMPr4u5NfXk1GYSWlVQZ2J+fz8S/n2HuuABMQ7OWEg51lbkJzXSitpTxgFD6xz4FbgPI8qCJP6SF3Oh52v6sMXLVzVHrKaZvGq9PoiOgSwa29buWaoGuoraslpSSFKmMVJwpO8N2Z71h/bj3lteV0c+mGu976z4Yk+Sja2+fJ1l3uenbzcCRuSCD9Atw5mFZEUaXyPOjfu9OoqjUSFdwF+1aOF5Tk00odKfnU02mVgWm3RQdz46AAHPU6UvMrqKgxkl5YSfzxHD7acY7EzBLQQFAXp1bfiOZouL7dAiEoGoY9AN3HASbIPwuGKsg9AcdWwb5/Kp0W7J3APQi0TeP1dfZlQsgEZkTOwM/Zj7zKPPIq8yiqLmJP9h6+OPEFO8/vpKK2Aj9nP6sNXJXko2ivnydbdaXrqdFo6OHryp0xITjYaTmWWUxptYE95wr44fB59HZaIru5YXeJz9LVsMT9pTF11ClTr0JCQgIAUVFRTY4lJSUBEBkZadWY2kKtsY7tp3NZfSCTTSdyqDbUNRxzsNNyba+uTOrXjQm9fS3eNtycy17f6lI4/j0cXgEpO4Df3KJOXspqq71vUJKVffPdrU8WnOT7s9+z/tx68irzGh0b2HUgE0ImMD54POEe4a3/gZpxuXusLal13uZ0pM+TLWjp9SyqqOHD7cl8sC2ZWqPyeerp68rzU3szPtK3xU3Tlri/JPnQ8ZPPb5VWKd2z1x7JYvvp3IYbEUCrgejQLoyL9GVcZFf6+ru32fOSq76+Jech8TvlK2Nf42N2TkoCipgIPa9XmucuwVhnZG/2Xn5K+YlNaZsori5udDzUPZRrgq5hTOAYov2icdA5mPlTNSXJR9FRP09qMfd6nrlQyntbk/nuUCbGizPnDw7x5LkpvRne3fuqy5Hk00qdMfn8VnFlLZtP5BCfmMPWU7lU1hobHfdxdWBshA+jengzuqcPAZ6WG9Rp1vUtzoATP8DJHyH1FzDVNT7u3RO6j1cmPQ0bA05dmhRRW1fL/uz9bErdxJb0LVyovNDouKPOkehu0Yz0H8lw/+H06tILrcb8ZklJPorO8HmyptZez8TzxSxad5IdZ35tERgX2ZU/TuxFVLDnFd8vyaeVOnvy+a2qWiO7kvPZfCKHLUm5ZBQ2XT8kxMuZ4eFeDAv3YliYF2GtGEPQ6utbUQBnNsOpDXB2szKfXCMa8OsPYaMhZKQyzsitcft0namOEwUn2Jq+lR2ZOziWdwwTjT8O7np3hvgNYajfUAb7DqaPVx/sdVffa1CSj6KzfZ7amqWu557kfBatP0lCelHDvqkDunHL4CAm9mm+OU6STytJ8rk0k8nE2dxydpzOZceZfPYk51NabWjyOm8XPYNDPBkc0oVBQZ4MCPTAw/nq/mO26PWtM8L5BGXMUPIWZZ2hSy3P7RmirDkUNAwCh0C3AY2eGRVUFbAnaw87z+9kb9Zezpefb1KEg86Bft79GNR1EP19+tPfp/9lB7dK8lF05s9TW7Dk9TSZTMQfz2HpptMczypp2N/X3503Z0bRy6/pGl2SfFpJks/VMRjrSDxfwu7kfPalFLAvpZDiytpLvjbU25n+AR70DXCnb4A7fbq54+fu0OQ/5za9vjUVSgJK2QGpuyBzv9J77n9pdODbB7oNVGbm7jYA/PqBUxdMJhMZZRnsy97HgZwDHMg5QGZZ5iVP18WhC09GP8ktEbc0OSbJRyGfJ8tqi+tprDPx7cEMnln167L2TvY61j0xlnAfl0avtcT9ZUMrlAlbZafTMijYk0HBnjx0bQ/q6kyczS3jQGohB9MKOZxezKkLpZhMkJpfQWp+BT8ezWp4v6ezPb183Yjwc6WXnxs9fV2h3IC3cxuNNdI7Kx0Ruo9T/m2ogewjkL5XSUqZB5RpfkxGZTG8nGNw+Mtf3+8WgMa3D8G+fQjuGklcyI0w5GlyMZKQm8DhC4c5mneU4/nHqTJWUVhdyMbUjZdMPkK0FzqthtuHBhM7wJ81hzJZ/t8zZBVXcfx8SZPkYwmSfESLaS+OJYrwc2NmjLJMQnm1geNZJRzJKCbxfDHHz5dw5kIZhjoTRRW17E0pYG9KQaNyHO00dO+aS7iPC2E+zoR6uRDs5UyItzPd3B0tN1O3nR6ChipfPKrsK89TmuqyEpTZFbKPQuE55VjpeeXr7OZGxXR18uJ6755c79UdvPthCIglWa8nRVPH4JBrLROrECpzdbDj7hGhTBvgz6H0Qq6JaJsF61RPPmvXruXdd98lPT2dwMBAHnroIW6++eZmX19eXs7rr79OfHw8FRUVDB06lBdeeIGwsDDrBS2acHGwY1iY0hGhXrXBSHJuOSezSziVU8ap7FLO5JaRVlCByQRVBhPHs0oatTPXs9Nq8Pd0JNDTiQBPJwI8nPD3dMTfwxE/d0e6uTvSxVlv/qJ6Lj5KF+2Iib/uqy6FCyeVmlDuSbhwQlmXqPRiLa6yQKk5ZexVYgR6XfxizJMw8SXzYhHCBnVx0XNdb782K1/V5LNu3Trmzp3LPffcw5gxY9i0aRPz5s3D0dGRKVOmXPI9Tz75JEePHuXZZ5/FxcWFZcuWMXv2bH788Ufc3Jo+GBPqcbDT0cffnT7+jae1qao1smV/IpkltVQ7eJCaV8G5/HLSCyrIKlaezRjqTKQXVJJe0LTXXT17nYaurg50dXdUvrvp8XF1wNtFj/fF716uero46/F0tr/ylEIObhA8TPlqFHAJ5J9WZl7IP6vMtlCQrNSUKpSR5hSlt/j6CNGZqZp8lixZQmxsLPPnzwdg7NixFBcXs3Tp0ksmn/3797N161Y+/PBDrrnmGgCGDh3KhAkTWLFiBQ8++KBV4xfmcbTXEe7lQLiXA5GRPRsdq6o1kllUSUZhJRmFFZwvqlXsAcgAABgdSURBVCSzsJLzxVVkFVeSU1xNjVEZ31NrNHG+uIrzxZfoTHAJLnodns56PJzs8XS2x8NJ+XJ3ssfd0Q43R3vcLn53dbDDzdEOVwc7XBwccPONwiFgSNNebdWlUHZBWSpcCHHVVEs+6enppKWl8dRTTzXaP3nyZNavX096ejrBwcGNjv3yyy+4uLgwevTohn1eXl4MGzaMbdu2SfLpABztdfTo6kqPrpeee81kMlFYUUt2cRUXSqu4UFJNblk1uaXK9/yyavLKaigor6Gwoobf9uUsrzFSXlNJZlHztanL0WrARW+Hs4MOF70dTnodTvY6XB3tuG+0O9f2apu2cSE6ItWST3JyMgDh4Y3n1QoNVaZIOXfuXJPkk5ycTGhoKDpd4+aTkJAQ1q9f34bRCluh0WjwctHj5aKnL5efpdpgrKO4spbCihoKK2opLK+hqKKW4spaiiprKK6spbjSQEllLaVVtZRUGSitqqW0ykBFjbFJeXUmKK02XBzzVN3oWK2xTpKPEC2gWvIpLS0FwNW18V+4Li5Kl76ysrIm7ykrK2vy+vr3XOr1V8NoNDb0mf+tiooKgEseE61n7evrDrjbQaj7xX+gAfQXv5oy1pmorK2j4uJXZa2JSkMdlbV1VBlMF78r28Y6E9eGO1/yZzEajU3+WLKW5u5tNcjnybLUvp6WuK9VSz5XGtuqvcRU35d7z6VeL4S5dFoNrg46XB3UXfdIiI5KteRT3zOtvLy80f76Gsyleq65urqSkZHRZH95efkla0RXQ6fTXXKUsIzIblud5frWjwRXQ3P3tho6y+/bWtS+npa4r1WrLtQ/60lLS2u0PzU1tdHx/31Penp6kxpQamrqJV8vhBDCNqmWfEJDQwkKCmLDhg2N9sfHxxMWFkZAQECT94wZM4aSkhJ27tzZsK+goID9+/czatSoNo9ZCCGEZag6zmfOnDnMnz8fDw8Pxo0bx+bNm1m/fj1LliwBlMSSlpZGz549cXV1ZdiwYcTExPDUU08xd+5cPD09efvtt3Fzc2PWrFlq/ihCCCFaQNXkExcXR01NDR999BErV64kODiYxYsXM3XqVAC2bNnC/Pnz+eyzzxg+fDgAy5Yt429/+xuvvfYadXV1REdH8+abb+Lh4aHmjyKEEKIFVJ/bbebMmcycOfOSx+Li4oiLi2u0z8PDg0WLFrFo0SJrhCeEEKINdOr1fA4cOABwyf7qRqOx2WOi9TrL9a3/OaOjo6163svd22roLL9va1H7elrivla95mOr5EPStuT6di7y+7asjnA9O3XNRwghhDpkWgAhhBBWJ8lHCCGE1UnyEUIIYXX/3969R0VVvX0A/8KkgqAsQYUABcEGuSQhIIxaDqIoC6HUpVTGJcHUvGRZpmb2M7UUL1hqQRiYKGqkRip3BQwvIGCQKyNTxDHBEFIc7jOc9w/eOYvjGe54QH0+a/GHe2/OefbGmWfOPnvOpuRDCCFEcJR8CCGECI6SDyGEEMFR8iGEECI4Sj6EEEIER8mHEEKI4Cj5EEIIERwlH0IIIYKj5POIsrIyrF27Fm5ubnBwcMDMmTORkJDAaaNQKLBz505MnDgR9vb2ePPNN1FQUNBDET+5tmzZgsDAQF55Tk4OrKyseD8LFiwQPkjSbgqFAqNHj+b93RwcHNg2mZmZmDVrFuzt7TFp0iRERkb2YMS919WrV2Fra4vS0lJOeXvG7/fff4efnx8cHBwwYcIE7NixAw0NDUKF3m70VOtm6uvrERwcjIcPH2LZsmUYOnQokpKSsHz5ciiVSkyfPh0AsGnTJhw/fhwffvghjI2NERUVhcDAQMTFxWHYsGE93IsnQ3R0NCIjIyGRSHh1hYWF6N+/P6KiojjlAwcOFCo80glFRUWoq6vDli1bYG5uzpZrajZ9xs3Ly8PChQvh6emJ9957D7m5uQgJCQHDMAgKCuqhqHuf69evY8GCBVAoFJzy9oxfcXExAgMD4eDggJ07d+L69esIDQ2FXC7HunXreqI7LWMIKyUlhRGLxUx+fj6nPCgoiPHx8WEYhmFkMhljbW3NxMTEsPV1dXWMVCpl1q1bJ2i8T6LS0lLmgw8+YEaNGsU4OjoyAQEBvDZr165lZs+eLXxwpEt++eUXZtSoUUx1dbXa+oCAAN7fNSQkhHFycmLq6uqECLFXa2hoYA4cOMA4ODgwY8eOZcRiMVNSUsLWt2f81qxZw0ycOJEzngcPHmSsra2Z0tJSYTrSTjTt1oyOjg58fX3x4osvcsotLCxw69YtAMDFixehVCoxdepUtr5v376QSqU4e/asoPE+iUJDQ/HHH38gKioK1tbWattcvXoVVlZWAkdGuurq1asYPnw4tLW1eXV1dXXIycmBh4cHp3zq1KmorKxEXl6eUGH2Wrm5udi2bRvmzZuHDz/8kFPX3vE7d+4c3Nzc0LdvX7bNtGnToFQqkZmZ+fg70QGUfJqRSCT4/PPPoaGhwZY1NDQgIyMDL7zwAgDgxo0b0NPTg76+Pud3zczMcOfOHdTW1goa85MmODgYp06dgqurq9p6pVKJa9euobS0FDNmzICdnR2kUikiIyPB0NZTvVphYSH69u2LoKAgODg4wNnZGevWrYNcLodMJkNDQwNGjBjB+R0zMzMATVN2zzpLS0ukpqZiyZIlvM3i2jN+NTU1KCkp4bXR19eHrq5urxvjZ+aej0KhQGxsbIv1Q4cOhbu7O69869atuHnzJvbs2QMAkMvl0NXV5bXT0dEBAFRVVUFLS6ubon5ytHd8R44c2epxbt68idraWhQVFeGDDz7AoEGDcPr0aYSEhEAul2PZsmXdHTrpJn/++Sfkcjlmz56NhQsX4sqVK9i1axf7twTAe+2oXjdyuVzweHubwYMHt1j38OFDAK2PX0ttVO162xg/M8mnrq4O//vf/1qsHzt2LCf5MAyDrVu34ocffkBQUBAmT57Mlrem+VXTs6Sj49sSQ0NDREREwNraGkOGDAHQdEVaW1uLiIgIzJs3T+2Li/S80NBQ6OnpsVOmzs7OMDAwwEcffYRz5861+ruqRQlEvbbedzQ1NdvVpjd5ZpKPjo4OCgsL29W2vr4eq1atwqlTpxAUFISVK1eydbq6uqiqquL9jupTxbP6xtiR8W2Nrq4uXnnlFV65VCpFbGwsioqKePfkSO8wduxYXplUKuX8+9HXjup1M2DAgMcW19NANT6tjZ/qvael96feNsa9KxX2AnK5HG+//TYSEhKwZs0aTuIBmhYf3L9/Hw8ePOCUFxcXw9TUlHOjj3RcYWEhYmJieN9LUN1LGzRoUE+ERdpQXl6O2NhYyGQyTrnq72ZgYACRSMQu3FFR/fvR+xSEa/jw4W2On46ODgwNDVFcXMxpU15ejqqqql43xpR8mlEqlVi0aBHy8/MRGhqKgIAAXptx48YBAJKSktiy+vp6pKens3Wk84qLi7F+/XreysH4+HiYmprCxMSkhyIjrdHQ0MC6detw4MABTnl8fDxEIhHGjRsHJycnJCcnc6aHkpKSMGDAANjZ2Qkd8hOlX79+7Rq/8ePHIy0tDfX19Zw2IpFI7ZVpT3pmpt3a4/Dhw8jOzoavry+MjIzw22+/sXUaGhqwt7eHiYkJZsyYgY0bN6K6uhpmZmaIiopCZWUlgoODezD6p4NUKoWtrS0+/fRTVFRUwMjICCdOnMCZM2ewa9euZ/aeWm+nr6+PuXPnIjo6Grq6unByckJubi7CwsIwd+5cmJmZYdGiRXj77bfx/vvvY8aMGbh8+TK+//57rFixQu3ybMLVnvFTrSZ95513EBAQgJs3b2LHjh2YM2cOjI2Ne7gHXBoMrV9l+fv7IysrS22dSCTCH3/8AaDpSmfbtm04efIkqqurYWtri5UrV8Le3l7IcJ94fn5+EIlE2LdvH6e8oqICO3fuREZGBioqKvDCCy/g3XffZRd9kN6poaEB+/btw9GjR/HPP//A0NAQc+bMQXBwMHuzOyUlBV9//TWKiopgaGiIuXPnYt68eT0cee9z7NgxrF69GhkZGTAyMmLL2zN+OTk5CAkJwdWrVzFo0CC89tprWLp0Kfr06SN0N1pFyYcQQojg6J4PIYQQwVHyIYQQIjhKPoQQQgRHyYcQQojgKPkQQggRHCUfQgghgqPkQ9RatWoVrKysYG1tjYqKihbbvfrqq7CyssKqVasEjI50N39/f1hbW6OsrKzFNg8ePICdnR2WLFnClt29excuLi7sY3VU/29u377d5jkbGxtx6NAhzJkzB2PGjMFLL70ELy8vhIaGsk9oflR8fDz8/Pzg7OyM0aNHw8PDAxs3buTEnZWVpXYbdnU/QNP3YqRSKaqrq9Wes6qqCnZ2dti9ezesrKwQERHRar+WLVsGOzs7/Pfff5zyJUuW4NtvvwUA7Nq1q83Y/Pz8OtyX27dvt9lu06ZN7er340ZPOCCtamxsRFpaGmbNmsWrk8lk+PPPP3sgKtLdvL29kZWVheTkZMydO1dtm+TkZDQ0NMDHx4ct27RpE7y8vDq1ffzKlSsRHx8PT09PeHt7Q1NTE1euXMHevXuRmJiIw4cPc57lFxoairCwMEilUixevBhaWlq4du0afvrpJ8THx+PIkSMYNmwYLC0tERISwjuXhYUFFi5cyIvDyckJI0eOxO7du3nPcgSAS5cuoU+fPggMDERERASSkpIwf/58tX2Sy+XIyMjAhAkTOLGnp6cjLy+PjWvKlCkYPny42mOcPHkSZ8+ehbm5eYf70rxPc+bMUVtnaWnZrn4/dj20gyrp5T7++GNGLBYzkyZNYhYtWqS2TWRkJOPq6sqIxWLm448/FjhC0p0qKysZOzs7xs/Pr8U2gYGBnC2bs7OzGRsbG+bOnTtsG9X/G5lM1ur5cnNzGbFYzERGRvLq4uPjGbFYzISEhLBld+7cYaytrZkNGzbw2ufl5THW1tbMe++91+L5xGIx89Zbb7VYf+nSJcbGxoa5desWr27Tpk1McHAwwzAMs3z5ckYsFjO3b99We5y4uDhGLBYzp06dYsuUSiUzZcoUJiwsrMXzq1y5coWxtbVlpFIpU1lZ2eG+yGSyDr0eW+v340bTbqRV7u7uOH/+vNodWlNSUjBp0qQeiIp0twEDBsDNzQ05OTkoLy/n1ZeXlyMrKwvTpk1jn9y+b98+ODo64vnnn+/w+S5fvgyg6UGYj/L09MTQoUM5z1bMz8+HUqlU297BwQGjR4/mtO8oJycnDB8+nPdgVAC4cOEC+9Bg1VVf8wcLNxcfHw9dXV3O3lVnzpxBcXExvL29W42hpqYGK1asgFKpREhIiCBbILTW78eNkg9p1eTJk1FTU4Pz589zysvLy3H58mXenvLkyeXt7Q2lUomUlBReXWJiIpRKJfvmW1JSgrS0tE4/b0+1A+ePP/6IxsZGXv3p06dx8OBBXvvjx49zntissn//fqSnp3cqFhUPDw8cPXqU80GrrKwMf/31FyQSCQCw02nqkk9lZSUyMzPh4eGBfv36seUxMTGwsbFp88GemzdvRlFREebPnw9nZ+cu9aUj1PVbCJR8SKscHR3ZraybO336NLS1tdkXJXnyTZw4EXp6ekhMTOTVxcfHw8TEBE5OTgCAX3/9FUqlkrdZXHt5eHhAT08P0dHRmDx5MrZs2YKzZ8+yN78f3RfLxcUFpqamSEpKgpubG9avX4/U1FR2X63u2EfLxcUFDx8+RF5eHlt2/vx5GBgYsDf0+/TpA09PT+Tn56OkpITz++ruidXU1CA7OxsTJ05s9dzp6ek4fPgw7OzssHTp0i73pb6+HhUVFbwfdRvNqeu3ECj5kFaJRCK4ubkhLS2N8wk1JSUFUqmUNs97ivTt2xdTp07FpUuXOCu17t69i9zcXEyfPp3d0iI3Nxf9+/fv1EIDoGkLhoiICAwbNgz//PMPIiMjMX/+fIwdOxYLFy5EQUEBL7a9e/fCxsYG9+7dQ0xMDBYvXgyJRAJ/f39kZmZ2vuP/TywWA2haBaZy4cIFSCQSzlYe3t7eYBiGd/WTkJAAIyMjuLi4sGX5+floaGhgk5c6FRUV+OSTT6CtrY1t27Z1y9OnT506BYlEwvvZsGEDr626fguBkg9pk7u7O8rLy9k5dblcjgsXLtAWB08hHx8fKBQKpKamsmUJCQlgGIbziV4mk8HExKRL+yvZ29sjMTER4eHh8PX1hampKRoaGpCWlgZfX1+cOHGC037EiBE4duwY9u/fD39/f1haWkKpVCIrKwtBQUH47rvvOh0LAAwePBja2tqcZeLnz5/nXd2PGTMGw4YN4ySfiooKXLx4EV5eXuz2EQDYJeimpqYtnnfNmjW4d+8eVq9e3W27jU6YMAFRUVG8H3V7jqnrtxBoqTVp0/jx46GlpYUzZ85gzJgxyMjIgKamZptTCeTJ4+TkBGNjYyQmJmL27NkAmj5F29jYYOTIkWy7+/fvY+DAgV0+33PPPQepVMpO3924cQMxMTGIjo7Gxo0bMWXKFGhpabHtNTQ04OLiwl5d3LlzB0ePHkV4eDi++uorvPrqqzA0NOx0PLq6uuxV3/Xr13H37l21OxRPnz4dYWFhuHv3LgwNDZGYmAiFQsFJ0EDTOKmOq86RI0eQlpYGd3d3+Pr6djruRw0ZMqRDOys377dQ6MqHtElbWxvjx49n7/ukpKRg3Lhx7E1g8vTQ0NCAl5cXLl68iPv370Mmk6GgoID3pqqpqal2oUB77d69G7GxsbxyCwsLrF27Fq+//jru37+Pv//+GwAQHR2NvXv38tobGxtj6dKlWLZsGRQKRZdWvAFN32sTiUQAmq56zM3N1S4UUE29JScnA2i6OhSLxRg1ahSnneoqSN1Y3bx5E5s3b8aQIUOwcePGLsXdVc37LRRKPqRdJk+ejBs3buCvv/7C2bNnMWXKlJ4OiTwmqqm3tLQ0JCcnQyQSwcvLi9PGwMCA/VTfGT///DPCw8PBtLCXpeo+hGp76NTUVOzZs6fFb+Or2je/SuqMBw8ewMDAAID6KTcVS0tL2NraIjk5GeXl5cjJyeElaADssR4dK4VCgY8++gg1NTX48ssvoa+v36W4u6p5v4VCyYe0i5ubG0QiEbZs2YLa2lr6fs9TTCwWw8rKCunp6UhLS4OrqyuGDh3KaWNsbIx///0XSqWyU+fw9vaGTCZDeHg4r66urg4///wzzM3NYWFhwbavrq7G5s2beVcRjY2NiI2NxcCBA7u0RLmsrAwKhQLPP/88FAoFsrOzW5268vHxQV5eHntvSt33eExMTAAApaWlnPJvvvkGBQUF8Pf3x8svv9zpmLtD834Lie75kHYZNGgQHB0dkZmZCRcXF86jQ8jTx8fHB9988w1qa2vxxRdf8OpdXV1x7NgxXLt2jTfVBDQ9CkfdtKynpyckEgkWLFiArKwshIaGIj09He7u7tDX10dJSQlOnDiB0tJSREZGsgsaZs6ciV9//RVHjhzB5cuXMW3aNBgZGaG8vBwJCQkoLCzE9u3b0b9//073OT8/HwAgkUhQUFCA6upquLq6ttjey8sLISEh2LNnD5ydnWFkZMRrY29vj/79+yM/P5+9erxx4wbCwsKgo6MDsViMuLg4tcfX0dERZFFP834LiZIPaTd3d3dkZ2fTF0ufAdOnT8f27dvRr18/tVOsL7/8MjQ1NZGTk6M2+Zw8eVLtcS0sLCCRSKClpYX9+/fj0KFDSEhIwN69e1FVVQV9fX2MGzcOCxYs4Kz80tTUxM6dOxEXF4e4uDgcOHAADx8+hJ6eHhwdHbF+/XqMHj26S33Ozc2Fnp4eXnrpJXz77bewtbVtdVHFkCFD4OrqinPnzqmdcgOaloi7uLhwljGXlZVBqVSiqqoKn3zySYvHNzExEST5NO+3kDSYliZdCSGkFYsXL0ZFRQUOHTrU06F0WWNjI9zc3DBt2jSsXr26W4+dmpqKxYsXIzk5GWZmZt167K56nP1uC93zIYR0yrx585CXl4dbt271dChdlpWVhXv37iEgIKDbj+3u7g5zc3McP36824/dVY+z322h5EMI6RRHR0e4ubl1+cudvUF4eDjeeOONNp+/1hkaGhpYsWIFDh06BLlc3u3H74rH2e+2UPIhhHTaZ599huTk5Cf66ic7OxtFRUVYvnz5YzuHh4cHxowZg3379j22c3SUEP1uDd3zIYQQIji68iGEECI4Sj6EEEIER8mHEEKI4Cj5EEIIERwlH0IIIYKj5EMIIURw/wfH2u/QMxfuSwAAAABJRU5ErkJggg==\n" | |
}, | |
"metadata": {} | |
} | |
] | |
}, | |
{ | |
"metadata": { | |
"trusted": true | |
}, | |
"cell_type": "code", | |
"source": "", | |
"execution_count": null, | |
"outputs": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"name": "python3", | |
"display_name": "Python 3", | |
"language": "python" | |
}, | |
"language_info": { | |
"name": "python", | |
"version": "3.6.2", | |
"mimetype": "text/x-python", | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 3 | |
}, | |
"pygments_lexer": "ipython3", | |
"nbconvert_exporter": "python", | |
"file_extension": ".py" | |
}, | |
"gist_id": "ac8278bd7cd93aa85718973c284990e8" | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 2 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment