Created
August 29, 2023 00:18
-
-
Save christinahedges/432f29e472425365196148837fe88962 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": [ | |
{ | |
"cell_type": "code", | |
"execution_count": 1, | |
"id": "54eec79b-8dd8-4cf8-8f0d-6a4f0a437e52", | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"import pandorasat as ps\n", | |
"import matplotlib.pyplot as plt\n", | |
"import numpy as np\n", | |
"import astropy.units as u\n", | |
"from astropy.coordinates import SkyCoord\n", | |
"plt.style.use('seaborn-v0_8-white')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 15, | |
"id": "39281541-535b-4709-901f-5a30191e2c77", | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# Find a place to point\n", | |
"c = SkyCoord.from_name('GJ 436')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 16, | |
"id": "80006b0d-e134-480c-b838-127b4170baa5", | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"# Point the telescope\n", | |
"p = ps.PandoraSat(ra=c.ra, dec=c.dec, theta=10*u.deg)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 18, | |
"id": "d3ad3182-6361-4bf4-916a-c0f71b17d8c7", | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/plain": [ | |
"<matplotlib.legend.Legend at 0x35eb7e790>" | |
] | |
}, | |
"execution_count": 18, | |
"metadata": {}, | |
"output_type": "execute_result" | |
}, | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAi8AAAG7CAYAAAALy3WMAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/P9b71AAAACXBIWXMAAA9hAAAPYQGoP6dpAABp90lEQVR4nO3dd1jV5t8G8PuwpwKK2+IEFyqKIIoDRFAQ91Yc1bq11WqdrdZZq3Vr0bq3IooK4sKJA9xbq+IeiAOQvfL+4Y+8UoaMA8mB+3NdXHJyniTfgIGbJ0+eKARBEEBERESkItSkLoCIiIgoJxheiIiISKUwvBAREZFKYXghIiIilcLwQkRERCqF4YWIiIhUCsMLERERqRSGFyIiIlIpDC9ERESkUopEeElISEC7du0QFBSU43UPHDgADw+PTN9fu3YtHB0d81IeERER5UChDy/x8fEYN24cHj58mON1L168iN9++y3T91+8eIEVK1bkpTwiIiLKoUIdXh49eoTu3bvj+fPnOV53xYoV+OGHH1CxYsVM20yfPh01a9bMS4lERESUQ4U6vAQHB8PW1ha7du1K997ly5fRuXNn1K1bF+7u7jhy5Eia98+dO4d169bB2dk5w237+PggNjYWXbt2zZfaiYiIKGMaUheQn3r37p3h8rCwMAwdOhRjx45Fs2bNcP36dUyaNAklSpSAtbU1AGDHjh0AkOE4mY8fP2LhwoXYsGEDbt26lX8HQEREROkU6vCSmW3btqFJkybo27cvAMDMzAz37t3Dpk2bxPCSlblz56JTp06oXr06wwsREVEBK5LhJSQkBCdPnoSVlZW4LDExEZUrV/7mumfPnsX169cxe/bs/CyRiIiIMlEkw0tSUhLc3d0xbNiwNMs1NL795Th06BDevn0LOzs7cVuJiYmwsrLCP//8k62eGyIiIsq9IhleKleujGvXrsHMzExctn79eiQkJKQLNP81fvz4NG2OHj2KLVu2YMuWLShdunS+1UxERERfFOq7jTLTu3dv3L59G4sXL8bTp09x8OBBLFq0COXKlfvmuiVKlICZmZn4UaJECWhoaMDMzAw6OjoFUD0REVHRViR7XsqXLw9PT08sXLgQ69atQ+nSpTFp0iS0b99e6tKIiIjoGxSCIAhSF0FERESUXSrT87J3715Mnjw53XKFQoH79+9LUBERERFJQWXCi6urK5o1aya+TkpKQv/+/dGyZUvpiiIiIqICpzLhRUdHJ82A2NWrV0MQBIwfPz5dW2trayQkJMDU1LQgSyQiIqI8CAsLg5aWFi5fvpxlO5W82yg8PBz//PMPfv75Z2hpaaV7Pz4+HklJSRJURkRERLmVlJSE+Pj4b7ZTmZ6Xr+3YsQOlSpVCmzZtMny/VKlSAICAgICCLIuIiIjyoFWrVtlqp3I9L4IgwMvLS3wuERERERUtKhdebt26hdDQULi5uUldChEREUlA5cLL2bNnYW1tjeLFi0tdChEREUlA5cLLzZs30aBBA6nLICIiIomoXHh5+PAhqlWrJnUZREREJBGVCy/v379HsWLFpC6DiIiIJKJyt0rfvHlT6hKIiIhIQirX80JERERFG8MLERERqRSGFyIiIlIpDC9ERESkUhheVISFhQV+/vnndMv37t0LR0dH8bWHhweWL18OAHj58iUsLCzSfNSuXRv29vaYNWsWEhIS0mw/9aNGjRpo2LAhBg4cmOkA6b1798LCwgJeXl5KPlIiIqKsMbyoEF9fX1y4cCHH63l5eSEwMBCBgYE4cuQIxo0bh927d2PNmjVp2i1fvhyBgYE4ffo0tmzZglKlSqF///74999/023Tz88P3333Hfbv35/r4yEiIsoNhhcVUr58ecycOTNNj0l2mJiYwNTUFKampqhQoQI6d+6MDh064Pjx42naFS9eHKampihdujRq1aqFP/74A7Vq1cKiRYvStPvw4QMuXLiAkSNH4vLly3jx4kWej42ISC4EQcCDBw+QnJwsdSmUCYYXFfLTTz8hNDQU69aty/O2tLS0oK6unmUbhUKBbt26ITAwEHFxceLyw4cPw9DQEO3bt0epUqXY+0JEhUZKSgpGjRqFGjVq4Ndff5W6HMqEyk1Slx8EQUBMTEyB7lNPTw8KhSJH65QuXRpjxozB4sWL0a5dO1SsWDHH+xUEAcHBwTh48CB++OGHb7avVq0aEhMT8fTpU9SoUQPAl0tGLVu2hJqaGhwdHeHj44ORI0fm+HiIiOQkJSUFw4cPFy+ph4SESFwRZabIhxdBEGBvb4/z588X6H6bNm2Ks2fP5vgXvoeHB/bu3Ys5c+bA09MzW+u0a9dO3E9CQgJMTEzQr18/DBo06JvrGhoaAgCio6MBAG/evMHVq1cxcOBAAICzszN27NiBK1euwNraOkfHQkQkFykpKRgyZEianm1eNpIvXjYCVKrHQF1dHTNmzMCpU6fSjVnJzJo1a+Dj44NVq1bhu+++g7W1NYYNG/bNy0YAEBUVBQAwMDAA8KXXRVtbG/b29gAAGxsbFC9eHPv27cvlERERSSs5ORmDBg3CunXroKamBmdnZ3E5yVOR73lRKBQ4e/asSlw2StWgQQN06dIFc+bMweDBg7/Zvly5cqhQoQLMzMywevVqdOjQAfPnz8e0adO+ue6DBw+gqamJSpUqAfgSXuLi4tCwYUOxTXJyMg4fPoxff/0VOjo6uTomIiIpJCcnY+DAgdiyZQvU1dWxbds2REZG4ujRowwvMlbkwwvwJcDo6+tLXUaOjB8/Hm3atMnx4N3vvvsOo0ePxp9//gl3d3fUq1cvy/ap88hoa2vjyZMnuHv3LqZNmwZbW1uxzaNHjzB27FgcO3YM7u7uuToeIqKClpSUhP79+2P79u3Q0NDAjh070LVrV6xfvx4Ae17kjJeNVJSxsTHGjx+PV69e5Xjdfv36oWrVqpg5cyZSUlLE5REREQgLC0NoaCju3LmDn3/+GXfu3MFPP/0E4Euvi5GREXr06AFzc3Pxw9XVFdWqVYOPj4+Sjo6IKH8lJiaiT58+YnDZvXs3unbtCgBQU/vyq/Hrn48kLwwvKqxr166wsrLK8XoaGhqYNm0abt++DW9vb3H56NGjYW9vj5YtW2Lo0KFITEzEzp07UaVKFQBfwou7uzu0tLTSbbNXr144f/48QkNDc39AREQFIDExEb169cLu3buhqakJb29vdOrUSXw/dTwge17kSyEIgiB1EcrWqlUrAEBAQIDElRARkZwkJCSgZ8+e2LdvH7S0tODt7Y127dqlabN9+3b06dMHrVq1yvaNEaQc2f39zTEvRERUJCQkJKBbt244cOAAtLW1sW/fPrRt2zZdO/a8yB/DCxERFXqJiYno0aMHDhw4AB0dHfj4+MDFxSXDtgwv8sfwQkREhVrqGBcfHx9oa2tj//794lwuGWF4kT8O2CUiokIrKSkJffr0gbe3N7S0tODj45NlcAEYXlQBwwsRERVKSUlJ8PDwgJeXFzQ1NbF37160adPmm+sxvMgfwwsRERU6ycnJGDBgAHbu3AlNTU3s2bMHbm5u2VqX4UX+GF6IiKhQSU5Oxvfff49t27aJE9C1b98+2+tzkjr5Y3ghIqJCIyUlBT/88AM2b94MdXV17Ny5Ex07dszRNtjzIn8ML0REVCikpKRg6NCh2LBhA9TV1bF9+3Z06dIlx9theJE/hhcV0Lt3b/z8888ZvnfgwAE0atQIf/31Fzw8PMTliYmJWL58OVq1aoU6deqgZcuWmDdvHqKiosQ2kyZNgoWFhfhhaWkJV1dXbNq0CRlNvBwTE4P69eujd+/eyj9IIqI8SElJwYgRI7B27Vqoqalh69at6N69e662xfAif5znRQW4ublh8eLFSEhISPdcIX9/fzg7O6dbvnDhQpw/fx6zZ89GxYoV8eLFC8yZMwfPnj2Dp6en2K5t27aYOnUqgC/h5OLFi/jjjz8QGRmJ0aNHp9nmiRMnYGpqiqtXr+LFixeoWLFiPh0xEVH2CYKA0aNHY/Xq1VAoFNi8eTN69uyZ6+0xvMgfe15UQNu2bREbG4sLFy6kWR4VFYXAwMB0z+UAgH379uHHH3+EnZ0dKlSoADs7O8yYMQMnT57Eu3fvxHY6OjowNTWFqakpzMzM0KNHD0yZMgWrV69O95BFX19fODk5wdzcnE+QJiJZEAQBP/74I1atWgWFQoGNGzeiT58+edomw4v8MbyoABMTE9jZ2eHo0aNplh8/fhxGRkawtbVNt45CocDFixfTjJa3srKCn58fjI2Ns9xf+/btoampiTNnzojLIiIiEBgYCGtrazg4OMDHxyfDS0tERAVFEASMGzcOy5cvh0KhwLp169CvX788b5fhRf4YXv4nOjo604+4uLhst42Njc1W25xq164dAgIC0pxMhw8fhqurq3hb39f69euHLVu2wNHREdOnT8eRI0cQFxeHatWqQVNTM8t9aWtro0KFCnj06JG47OjRo1BXV0eTJk3QqlUrvHz5EpcvX87xcRARKYMgCJgwYQKWLFkCAFizZg0GDhyolG0zvMgfw8v/GBgYZPrx39HqpUqVyrTtf59QWqlSpQzb5ZSTkxNiYmJw6dIlAMDnz58RGBgId3f3DNuPHDkSCxYsQJkyZbB7926MGTMGzZo1g7e3d7a/Hl+HLD8/PzRp0gS6urqwtLREmTJlsG/fvhwfBxFRXgmCgMmTJ+Ovv/4CAHh6emLw4MFK2z7Di/wxvKgIAwMDtGzZUrx0dPz4cVSoUAF16tTJdJ327dtj586dOH/+PBYuXIjq1atj6tSpuH379jf3FxUVJYassLAwBAcHw8nJCcCXS1KtW7fG4cOH0/U0ERHlJ0EQMG3aNMyfPx8AsHLlSgwdOlSp++AkdfLH8PI/UVFRmX78t7fi3bt3mbb19/dP0/bp06cZtssNd3d3HD9+HIIgwN/fP8OBugBw//59/PHHH+JrY2NjuLu7Y8uWLShTpgwuXryY5X7i4+Px9OlTVK9eHcCXO5qSk5Px66+/olatWqhVqxa2bduG6OhoHDt2LFfHQkSUG9OnT8fcuXMBAMuWLcOIESOUvg/2vMgfb5X+H319fcnbfkuLFi0wefJkXLx4ERcuXMCUKVMybJecnIwNGzagffv2qFWrlrhcS0sLOjo6MDExyXI/Bw8ehEKhgIODAwDg0KFDsLOzS7e/kSNHwsfHJ0fTbhMR5dbMmTMxa9YsAMDixYvTTeegLAwv8sfwokK0tLTQunVrzJ8/H+bm5qhUqVKG7WrXro2WLVtixIgR+Pnnn2FlZYX3799j3759SEhISPM4+Li4OISFhQH4Ms/L2bNnsWjRIgwfPhwmJiZ4+fIlrl27hqVLl8Lc3DzNfnr06IG//voLoaGhKF26dL4dNxHRnDlzMH36dABf5rH66aef8m1fDC/yx8tGKqZdu3a4d+9epgN1Uy1ZsgQdOnTAihUr0LZtWwwdOhRRUVHYunVrmgHD/v7+sLe3h729PTp37gwfHx9MmzYNw4cPB/Cl18XY2BiOjo7p9tG5c2doaGhg//79yj1IIqKv/PHHH5g2bRoAYP78+ZnOOK4sDC/ypxAK4WQdrVq1AgAEBARIXAkREeXFggUL8MsvvwD40vuS2eVyZXrx4gW+++47aGlpIT4+Pt/3R/8vu7+/2fNCRESytHjxYjG4zJw5s0CCC8CeF1WgUuElISEBv//+Oxo1aoQmTZpg0aJFnOWViKgQWrp0KcaNGwfgyx1Gv/76a4Htm+FF/lRqwO7s2bMRFBSEdevWITo6GmPHjkW5cuXy9AAuIiKSl5UrV4oDcqdOnSoO1C0oqeEF+DLXS0azmJO0VOY7Eh4eDm9vb8yaNQt169aFnZ0dvv/+e9y4cUPq0oiISEk8PT0xatQoAMCkSZMwa9YsKBSKAq3h67DCierkSWV6Xq5cuQIDAwPY2NiIy4YMGSJhRUREpExr1qwR73ScMGEC5s6dW+DBBUjb85KcnAwNDZX5VVlkqEzPy4sXL1C+fHn4+PigTZs2aNWqFVauXMlUTERUCKxfv16c5n/cuHGYP3++JMEFSB9eSH5UJk7GxMTg2bNn2LlzJ+bNm4ewsDD89ttv0NXVxffffy91eURElEubNm0SH6z4448/YuHChZIFF4DhRRWoTHjR0NBAVFQU/vrrL5QvXx4A8Pr1a+zYsYPhhYhIRW3ZsgUDBw6EIAgYNWoUFi9eLGlwARheVIHKXDYyNTWFtra2GFwAoHLlynjz5o2EVRERUW5t374dAwYMgCAIGD58OJYtWyZ5cAEYXlSByoSXevXqIT4+Hk+ePBGXhYSEpAkzRESkGnbt2gUPDw+kpKTghx9+wIoVK2QRXIC0dxsxvMiTyoSXKlWqoGXLlpg8eTLu37+Ps2fPYs2aNejVq5fUpRERUQ54eXmhT58+SElJwaBBg+Dp6SmruVQUCoVYD8OLPKnMmBfgy5NEZ82ahV69ekFXVxd9+vSBh4eH1GUREVE27d27F7169UJycjIGDBiANWvWyCq4pFJXV0dKSgrDi0ypVHgxNDTEn3/+KXUZRESUC/v370ePHj2QnJwMDw8PrF27VpbBBfj/S0ecjkOe5Pm/hoiICpWDBw+iW7duSEpKQu/evbFhw4Y0A2Plhs83kjeGFyIiyleHDh1C165dkZiYiB49emDTpk2yDi4Aw4vcMbwQEVG+OXLkCDp37oyEhAR07doVW7duVYnp9hle5I3hhYiI8sWxY8fQoUMHxMfHo3Pnzti+fbtKBBeA4UXuGF6IiEjpTpw4gfbt2yM+Ph4dOnTAjh07oKmpKXVZ2cbwIm8ML0REpFSnTp1Cu3btEBcXh3bt2mH37t3Q0tKSuqwcYXiRN4YXIiJSmjNnzsDNzQ2xsbFwdXXFnj17VC64AAwvcsfwQkRESnHu3Dm4uroiJiYGLi4u8Pb2hra2ttRl5QrDi7wxvBARUZ5duHABbdq0QXR0NFq3bo19+/ZBR0dH6rJyjZPUyRvDCxER5UlQUBBcXFwQFRUFR0dH+Pj4QFdXV+qy8oQ9L/LG8EJERLl26dIlODs74/Pnz2jZsiUOHjwIPT09qcvKM4YXeWN4ISKiXLl69SqcnZ0RGRmJZs2awdfXt1AEF4DhRe4YXoiIKMeuXbsGJycnhIeHo2nTpjh06BD09fWlLktpGF7kjeGFiIhy5MaNG3BycsKnT59gZ2cHf39/GBgYSF2WUjG8yBvDCxERZdvt27fh5OSEjx8/wsbGBv7+/jA0NJS6LKVjeJE3hhciIsqWO3fuwNHREe/fv4e1tTWOHDmC4sWLS11WvmB4kTeGFyIi+qZ79+7B0dERYWFhaNCgAY4ePQojIyOpy8o3DC/yxvBCRERZevDgARwdHfHu3TvUr18fx44dg7GxsdRl5StOUidvDC9ERJSphw8fwsHBAW/fvkXdunVx/PhxmJiYSF1WvmPPi7wxvBARUYYePXoEBwcHvHnzBnXq1MHx48dRokQJqcsqEAwv8sbwQkRE6YSEhMDBwQGvXr1CrVq1EBAQAFNTU6nLKjAML/LG8EJERGk8ffoUDg4OePnyJWrUqIETJ06gVKlSUpdVoBhe5I3hhYiIRM+ePYODgwOeP38Oc3NznDhxAqVLl5a6rALH8CJvDC9ERAQAePHiBRwdHfH06VNUq1YNJ06cQNmyZaUuSxIML/LG8EJERHj16hUcHBwQEhKCqlWr4uTJkyhfvrzUZUmG4UXeGF6IiIq4169fw8HBAY8fP0blypVx8uRJVKhQQeqyJMXwIm8ML0RERdjbt2/h6OiIhw8fwszMDCdPnkTFihWlLktynKRO3hheiIiKqNDQUDg6OuLBgwf47rvvcPLkSZiZmUldliyw50XeGF6IiIqgd+/eoVWrVrh37x4qVKiAkydPonLlylKXJRsML/LG8EJEVMS8f/8eTk5OuHPnDsqVK4eTJ0+iSpUqUpclKwwv8sbwQkRUhHz48AFOTk64desWypYti5MnT6JatWpSlyU7DC/yxvBCRFREfPz4Ea1bt8aNGzdQunRpnDhxAubm5lKXJUsML/LG8EJEVASEh4fD2dkZ165dQ6lSpXDixAnUqFFD6rJki+FF3hheiIgKuYiICDg7O+PKlSsoWbIkAgICUKtWLanLkjWGF3ljeCEiKsQiIyPRpk0bXLp0CSVKlMCJEydQp04dqcuSPYYXeWN4ISIqpD5//oy2bdvi4sWLMDExQUBAACwtLaUuSyVwkjp5Y3ghIiqEoqKi4ObmhvPnz8PIyAjHjh1DvXr1pC5LZbDnRd4YXoiICpmoqCi4urri7NmzKF68OI4fP44GDRpIXZZKYXiRN5UKL8eOHYOFhUWajzFjxkhdFhGRbERHR6Ndu3Y4e/YsihUrhqNHj6Jhw4ZSl6VyGF7kTUPqAnLi0aNHcHBwwKxZs8Rl2traElZERCQfMTExaNeuHU6fPi0GFxsbG6nLUkkML/KmUuHl8ePHMDc3h6mpqdSlEBHJSkxMDNzd3XHq1CkYGhri8OHDsLW1lboslcXwIm8qddno8ePHqFSpktRlEBHJSmxsLNq3b48TJ07AwMAAhw8fhp2dndRlqTSGF3lTmfAiCAKePHmCwMBAuLi4wMnJCQsXLkRCQoLUpRERSSY2NhYdOnRAQEAA9PX1cfjwYTRp0kTqslQew4u8qcxlo9evXyM2NhZaWlpYsmQJXr58idmzZyMuLg7Tpk2TujwiogIXFxeHjh074tixY9DX14e/vz+aNm0qdVmFAsOLvKlMeClfvjyCgoJQvHhxKBQK1KxZEykpKZgwYQImT54s/kcjIioK4uLi0KlTJxw9ehR6eno4dOgQmjVrJnVZhQYnqZM3lblsBABGRkZQKBTi66pVqyI+Ph4RERESVkVEVLDi4+PRuXNnHD58WAwuzZs3l7qsQoU9L/KmMuHl7NmzsLW1RWxsrLjs3r17MDIygomJiYSVEREVnPj4eHTp0gX+/v7Q1dWFr68vWrRoIXVZhQ7Di7ypTHixsrKCtrY2pk2bhpCQEJw+fRp//vknBg8eLHVpREQFIj4+Hl27doWfnx90dHTg6+sLBwcHqcsqlBhe5E1lxrwYGBhg3bp1mDt3Lrp06QJ9fX307NmT4YWIioSEhAR0794dvr6+0NHRwcGDB+Ho6Ch1WYUWw4u8qUx4AYDq1atjw4YNUpdBRFSgUoPLgQMHoK2tjQMHDsDJyUnqsgo1hhd5U5nLRkRERVFiYiJ69uyJ/fv3Q1tbG/v370fr1q2lLqvQY3iRN4YXIiKZSkxMRK9evbBv3z5oaWnBx8cHLi4uUpdVJDC8yBvDCxGRDCUmJqJ3797w9vaGlpYW9u3bhzZt2khdVpHB8CJvDC9ERDKTlJSEvn37Ys+ePdDU1MTevXvh6uoqdVlFCiepkzeGFyIiGUlKSoKHhwd2794NTU1NeHt7w83NTeqyihz2vMgbwwsRkUwkJSWhX79+2LlzJzQ1NbFnzx64u7tLXVaRxPAibwwvREQykJycjAEDBmDHjh3Q0NDA7t270b59e6nLKrIYXuSN4YWISGLJyckYOHAgtm3bJgaXjh07Sl1WkcbwIm8ML0REEkpOTsb333+PLVu2QF1dHTt37kSnTp2kLqvIY3iRN5WaYZeIqDBJvVS0detWqKurY8eOHejSpYvUZREYXuSO4YWISAJJSUno378/tm/fLva4dO3aVeqy6H8YXuSN4YWIqICl3g69c+dOaGhoYOfOnexxkRmGF3ljeCEiKkBJSUno06cPdu/eLQ7O5RgX+eEkdfLG8EJEVEASExPRp08feHl5QVNTE15eXujQoYPUZVEG2PMibwwvREQFIPUhi97e3uLMuZyATr4YXuSN4YWIKJ8lJCSgZ8+e4tOhvb290a5dO6nLoiwwvMgbwwsRUT5KSEhAjx494OPjIz4dmg9ZlD+GF3ljeCEiyicJCQno1q0bDhw4AG1tbfj4+KBNmzZSl0XZwPAibwwvRET5ID4+Ht26dcPBgwehra2N/fv3w8XFReqyKJsYXuSN4YWISMni4+PRpUsX+Pn5QUdHB/v374ezs7PUZVEOMLzIG8MLEZESxcXFoXPnzvD394eOjg4OHjwIJycnqcuiHGJ4kTeGFyIiJYmLi0PHjh1x5MgR6OrqwtfXF46OjlKXRbnASerkLU/hpWbNmjleR6FQ4O7du3nZLRGR7MTGxqJjx444evQo9PT04Ofnh5YtW0pdFuUSe17kLU/hRRAELF++HMWLF89W+/DwcIwZMyYvuyQikp2YmBh06NABx48fh56eHg4dOoQWLVpIXRblAcOLvOUpvNjb26Nx48YwNDTMVvvIyEjY29vnZZdERLISExOD9u3bIyAgAPr6+jh06BCaN28udVmURwwv8qaWl5XXrl2baXDJ6JJSsWLFsHbt2rzskohINqKiouDm5oaAgAAYGBjg8OHDDC6FBMOLvOXbgF1BEPJr00REkouMjISrqyvOnTsHQ0ND+Pv7o2nTplKXRUrC8CJvOQovPj4+Wb7fsWNH8XOFQpGbeoiIZO/Tp09o06YNgoODYWRkhCNHjsDGxkbqskiJUsNLSkoKBEHg7zSZyVF4efDgQbplKSkpOHLkCMLCwtKEFyKiwujDhw9o3bo1rl27BhMTExw7dgwNGjSQuixSstTwAnz5Pff1a5JejsLLxIkTxc8TEhKwZ88ebN68GU2aNMEPP/yg9OKIiOTk3bt3cHJywq1bt2BqaoqAgABYWlpKXRblg6/DSnJyMsOLzOR4zEtMTAy2bduG3bt3w8HBARs3bkSZMmXyozYiItl48+YNWrVqhXv37qFMmTIICAhArVq1pC6L8knqJHUAJ6qToxyFl6VLl+LgwYNwd3fH7t27YWxsnF91ERHJxsuXL+Ho6IiHDx+ifPnyOHHiBMzNzaUui/LRf3teSF5yFF7+/vtv6OnpYceOHdi5cyeA/7+rSKFQ4MKFC8qvkIhIQs+ePYOjoyNCQkJgZmaGEydOoEqVKlKXRfmM4UXechRe7t+/n191EBHJzuPHj+Ho6Ijnz5+jSpUqOHHiBMzMzKQuiwoAw4u85XjMi5eXF7Zv344nT54AAKpUqYK+ffuic+fOSi+OiEgqDx48QKtWrfDq1SuYm5vjxIkTKF++vNRlUQFheJG3HIUXLy8v7NixA5MmTULt2rUhCALu3r2LP//8E4IgoEuXLvlVJxFRgbl79y4cHR0RGhqKWrVqISAggDcmFDFfD9hleJGfHIWXnTt3YtWqVWlOYltbW6xYsQKjRo3KMLycPn063TJ9fX1Uq1YNRkZGOa+YiCgf3bx5E05OTggLC0O9evVw7NgxmJqaSl0WSUBdXR3JyckMLzKUo/ASFxeX4V8fZcuWRVxcXIbrbN26FTdv3oStrS0AICgoCDVr1sSbN28wfPhwTmxHRLJx9epVtG7dGh8/fkTDhg1x9OhRmJiYSF0WSYThRb5yFF60tbUzfU9LSyvD5QkJCfD39xd/AHz8+BE///wzvLy80Lt371yFlyFDhsDExAR//PFHjtclIspIUFAQXFxcEBERgcaNG8Pf35+9w0Ucn28kXzm+28jOzi7dckEQEBUVleE6oaGhKF68uPi6WLFiePPmDYoVKwYNjZw/F9LPzw+nT59Gp06dcrwuEVFGAgMD4erqis+fP8Pe3h6HDh2CoaGh1GWRxFLHvXCSOvnJUXq4e/dujnfQpk0b9OrVC05OTgCAEydOoE2bNoiJiUHlypVztK3w8HD8+eefnI6biJTm1KlTaNeuHaKjo+Hg4ICDBw9CX19f6rJIBtjzIl857/rIoREjRsDR0RHXrl2DIAiYMmUK6tatCwBYvHhxjrY1f/58dOjQAe/evcuPUomoiDly5Ag6deqE2NhYODs7Y9++fdDT05O6LJIJhhf5yvfw0rFjRxw6dEgMLLl14cIFXL58GQcPHsSMGTOUUxwRFVn79u1Djx49kJiYiHbt2sHLyws6OjpSl0UywvAiX2rfbvLFo0ePsG7dOoSGhgIAjh07lq31atSogYcPH+auuv+Jj4/H9OnT8dtvv/GHCxHl2datW9GtWzckJiaie/fu2Lt3L3+2UDoML/KV7Z6XFStWYNiwYVi5ciV69+6Nc+fOoXXr1t9c7+HDh+jUqRMqVaoEHR0dCIIAhUKBPXv2ZLvIFStWoE6dOmjWrFm21yEiysiaNWswbNgwCIKAAQMGYO3atWlmUyVKxfAiX9kOL8WKFUONGjXw+++/Y+7cuXjw4EG21vP09Mx1can8/Pzw/v17WFlZAfhy+zXw5Xr1tWvX8rx9IioaFi9ejHHjxgEARo0ahaVLl6aZSZXoawwv8pXt8NKkSRMAX54ePWXKFKxbty5b65UvXx6nT5/G8+fP4eHhgbCwMERGRuaoyC1btiApKUl8vXDhQgDA+PHjc7QdIiqaBEHA7Nmz8dtvvwEAJk2ahLlz50KhUEhcGckZw4t8ZTu8tGnTRvxcoVBg8ODB2Vpv/vz5+PjxI27cuAEPDw+oq6tj0qRJ8PLyynaR/30YWuptjHy6KxF9iyAImDhxIhYsWAAAmD17NqZOnSpxVaQKGF7kK1d3Gz1//hxbt27Fq1ev0nxTM7pEdOHCBfj4+Igz6ZqYmCA+Pj531RIR5UBKSgpGjRqFv//+G8CXy0Y//fSTtEWRyuAkdfKVq/AycuRI9O3bF23btv3m9WINDQ2kpKSI3bPh4eF5vsbMxwIQ0bckJSVh0KBB2Lx5MxQKBdasWZPtHmMigD0vcpar8KKrq4sePXpkq62Hhwd++uknfPr0CStXroS/vz+GDRuWm90SEWVLQkIC+vTpgz179kBdXR2bN29G7969pS6LVAzDi3zlKrwMGTIECxYsQJMmTdI8kLFRo0bp2nbo0AF16tTBxYsXkZKSgiVLlqBatWq5r5iIKAuxsbHo2rUrDh06BC0tLezatYtPr6dcYXiRr1yFl8DAQAQHB+P58+fi5SCFQpFheAGAqlWromrVqrmvkogoGz5//owOHTrg5MmT0NXVhY+PD5ydnaUui1QUw4t85Sq8XLhwAYcPH+ZthkQkG+Hh4Wjbti0uXrwIQ0ND+Pn5cWJLyhOGF/nK1chZS0tLPH/+XNm1EBHlSlhYGBwcHHDx4kUYGxsjICCAwYXyjOFFvnLV8/Lw4UO0a9cOlStXhpaWVq6m/CciUobnz5/D2dkZDx48QKlSpXDs2LE8PwiWCGB4kbNchZdVq1Ypuw4iohx78OABWrdujRcvXuC7777DsWPHYG5uLnVZVEgwvMhXri4bbdu2DYaGhihfvjzKly8PAwMD7NixQ9m1ERFl6urVq2jWrBlevHgBCwsLBAYGMriQUnGSOvnKVXg5f/48ihUrJr4uXrw4zp07p7SiiIiycubMGTg4OCAsLAwNGjTA2bNnUbFiRanLokKGPS/ylavwkpKSgqioKPF1VFRUmgcnEhHlFz8/P7i4uCAyMhLNmzfHyZMnYWpqKnVZVAgxvMhXrsa89OvXD71794arqysEQYC/vz/69++v7NqIiNLYsWMH+vXrh6SkJLRr1w67d++Grq6u1GVRIcXwIl+5Ci9du3ZF3bp1ERQUBAD466+/UL169TRtBEHIe3VERP/z999/Y+TIkRAEAX369MGGDRugqakpdVlUiDG8yFeun5Bobm6OM2fOoGPHjumCCwDcv38/3bLIyEj88MMPud0lERVBgiBgzpw5GDFiBARBwMiRI7F582YGF8p3DC/ylauel1SBgYEICgpC8eLFs9U+IiICgYGBedklERUhgiBgwoQJ+OuvvwAAv/76K37//XfO7k0FguFFvvIUXgRBwKhRo3K0Dn/oEFF2JCUlYejQoVi/fj0AYNGiRRg7dqzEVVFRwvAiX3kKLxldGiIiyqv4+Hj07t0be/fuhZqaGtauXYuBAwdKXRYVMQwv8pWn8EJEpGxRUVHo3Lkzjh07Bi0tLezYsQOdO3eWuiwqgjhJnXwxvBCRbHz8+BFubm64ePEi9PX14ePjAycnJ6nLoiKKPS/yxfBCRLLw8uVLtGnTBnfu3IGxsTH8/f1ha2srdVlUhDG8yBfDCxFJ7t69e3BxccGLFy9Qrlw5HDlyBHXq1JG6LCriGF7kK9fzvPzXxIkTcebMGX6TiShHgoKCYG9vLz5g8fz58wwuJAsML/KltJ4XAwMDTJ06FYmJiXB2doarqytsbW15azQRZcrf3x9du3ZFTEwMbGxs4Ofnh5IlS0pdFhEAhhc5U1rPy6+//oozZ85g2bJl0NDQwPjx49GsWTPMmTMH169fV9ZuiKiQ2Lp1K9q3b4+YmBi4uLggICCAwYVkheFFvpQWXoAvE9DZ2Njgt99+w+HDh9G1a1fs3r0bvXr1QqtWrbB69WrEx8crc5dEpIIWLVoEDw8PJCUloU+fPjhw4AAMDAykLosoDYYX+VLqgN3o6GicPHkShw8fRmBgIEqXLo2BAwfC1dUVYWFhWLhwIYKDg7Fu3Tpl7paIVIQgCJg4cSIWLFgAABg7diwWLlwozqdBJCcML/KltPAyfPhwnD9/HsWKFUPbtm2xefNm1K1bV3zf3NwckZGRmDp1qrJ2SUQqJDExET/88AM2bdoEAJg/fz4mTJjAcXEkW5ykTr6UFl5KliyJ1atXZzlI19raGl5eXsraJRGpiOjoaHTv3h2HDh2Curo61q5diwEDBkhdFlGW2PMiX0rrq01KSkKdOnXSBZeIiAiMGTMGAGBqaoqqVasqa5dEpAI+fPgAJycnHDp0CLq6uvDx8WFwIZXA8CJfeep5uXbtGp49ewYA8PHxQe3atdMNugsJCUFgYGBedkNEKurFixdwcXHBvXv3YGxsDF9fXzRp0kTqsoiyheFFvvIUXnR1dbF8+XIIggBBELB27do0A+8UCgX09PQwfvz4PBdKRKrl7t27cHFxwcuXL1GhQgUcOXIEtWrVkrosomxjeJGvPIWXGjVqICAgAADg4eGBFStWoHjx4kopjIhU14ULF+Dm5oZPnz6hZs2aOHLkCCpWrCh1WUQ5wvAiX3kKL69fv0bZsmWhUCgwf/58REdHIzo6OsO25cqVy8uuiEhF7N+/H7169UJsbCwaN24MX19flChRQuqyiHKM4UW+8hReHB0dce7cOZQoUQKOjo5QKBQQBEF8P/W1QqHAvXv38lwsEcnb33//jVGjRiElJQWurq7YvXs39PX1pS6LKFcYXuQrT+ElICAAJiYm4udEVDQJgoApU6bgjz/+AAAMHjwYf//9NzQ0+OB6Ul0ML/KVp58s5cuXFz/39PSEm5sbH8ZIVMQkJCRg0KBB2Lp1KwBg5syZmDZtGn8OkMrjJHXypbQ/i2JiYjBy5Ejo6urCxcUFbdu2hbW1tbI2T0QyFBkZiS5duuD48eNQV1fHP//8g4EDB0pdFpFSsOdFvpQWXv766y8kJCQgMDAQx44dE4NM27Zt4erqCktLS2Xtiohk4PXr12jbti1u3rwJfX19eHt7w8XFReqyiJSG4UW+lHpBWktLC46OjnB0dERCQgI2btwIT09PbNy4kQN2iQqRu3fvok2bNnjx4gVKly6NQ4cOoUGDBlKXRaRUDC/ypdTwkpycjKCgIBw9ehTHjx9HSkoK3N3d4ebmppTtP3v2DDNnzsTVq1dRvHhx9O3bF4MHD1bKtokoe86cOYMOHTogPDwcFhYW8Pf3R+XKlaUui0jpGF7kS2nhZdKkSTh58iRSUlLg5OSEefPmoUmTJuI3P69SUlIwZMgQWFpaYt++fXj27BnGjRuH0qVLw93dXSn7IKKseXl5oW/fvkhISECTJk1w4MABzuFChRbDi3wpLbwkJCRgzpw5aN68ObS0tJS1WdH79+9Rs2ZNzJgxAwYGBqhUqRLs7Oxw5coVhheiArBkyRKMGzcOgiCgU6dO2LZtG3R1daUuiyjfMLzIl9Jm2E19ftH79+8zbJvXGXZLlSqFJUuWAPgyp8TVq1dx6dIlTJ8+PU/bJaKsJScnY8KECVi8eDEAYNSoUViyZInSelWJ5IrhRb7yZYbdr+d3yI8Zdh0dHfH69Ws4ODjw7gaifBQTE4O+ffti3759AID58+djwoQJnMOFigSGF/nK8wy7xsbG4ucFZdmyZXj//j1mzJiBefPmYdq0aQW2b6KiIjQ0FO3bt0dwcDC0tLSwadMm9OzZU+qyiAoMJ6mTL7W8rFy+fHnxm5v6ROny5cun+TAwMMD8+fOVUmwqS0tLODg4YPLkydi5cycSEhKUun2iou7u3bto3LgxgoODUaJECQQEBDC4UJHDnhf5ylPPy7Vr1/Ds2TMAgI+PD2rXrg0DA4M0bUJCQhAYGJiX3QD4Mpbm+vXrcHJyEpdVq1YNiYmJiIqKEp+xRER5c+LECXTu3BkRERGoVq0aDh06hOrVq0tdFlGBY3iRrzyFF11dXSxfvhyCIEAQBKxdu1bsiQG+PFVaT09PHMybFy9fvsSoUaNw+vRplC5dGgBw+/ZtmJiYMLgQKcmmTZswePBgJCUloWnTpvDx8UHJkiWlLotIEgwv8pWn8FKjRg1xrIuHh4d46Sg/WFpaonbt2pgyZQomT56MV69eYcGCBRg2bFi+7I+oKBEEAdOnT8esWbMAAD179sSGDRugo6MjcWVE0mF4kS+l3So9f/58REdHIzo6OsO2eb1VWl1dHatWrcKsWbPQo0cP6OrqwsPDA/369cvTdomKuvj4eAwaNAjbtm0DAEyZMgWzZs1K04tKVBQxvMhXvtwq/V/KulW6dOnSWLFiRZ63Q0RffPz4EZ06dcKZM2egrq6O1atXY9CgQVKXRSQLDC/yledbpVPHmxTkrdJElHePHz+Gm5sbHjx4gGLFimHPnj1o3bq11GURyQbDi3zlKbyUL18+3ecpKSlQU1PDu3fvcOXKFdSoUYMPbSOSmQsXLqB9+/Z4//49KlasiEOHDqFOnTpSl0UkKwwv8qW0i9pXrlxBs2bNEBwcjHfv3qFz58747bff4O7uDn9/f2XthojyaMeOHXBwcMD79+/RoEEDBAUFMbgQZYCT1MmX0sLL3Llz4erqinr16mH37t3Q1tbGuXPnMGvWLCxbtkxZuyGiXEq9o6h3796Ij49H+/btcfr0aZQtW1bq0ohkiT0v8qW08PLw4UP0798furq6OHHiBJydnaGlpQUbGxu8fv1aWbsholyIjY1Fr169MHPmTADAhAkTsHfv3nSTShLR/2N4ka88jXn5WsmSJfHo0SPExMTg7t27mDRpEgDg/Pnz/MuOSEJv375Fhw4dEBwcDA0NDaxevRrff/+91GURyR7Di3wpLbwMGDAAI0eOhJqaGiwtLWFjYwNPT0+sWLEC8+bNU9ZuiCgHbty4AXd3d7x48QImJibw9vZGy5YtpS6LSCUwvMiX0sJLv379YG1tjdevX8Pe3h4A0LhxY7Rs2RI1atRQ1m6IKJsOHjyIXr16ITo6GhYWFvD19UW1atWkLotIZTC8yJfSwgsA1KpVC7Vq1RJf169fX5mbJ6JsEAQBixYtwoQJEyAIAlq1agUvLy8YGxtLXRqRSmF4kS+lhZfIyEisX78et27dQlJSUrqZdjdv3qysXRFRJhISEjBixAisW7cOADB06FAsX74cmpqaEldGpHoYXuRLaeHll19+wa1bt+Du7s47GIgk8PHjR3Tp0gWnTp2CmpoaFi1ahDFjxkChUEhdGpFKYniRL6WFl/Pnz2Pr1q2oW7eusjZJRNn077//ol27dnj48CEMDQ2xc+dOuLq6Sl0WkUrjJHXypbR5XkqXLs2n0FKR9vnzZ4SFhRX4fo8cOQIbGxs8fPgQZmZmOHfuHIMLkRKw50W+lHrZaMaMGRgzZgzMzMzSXWMvV66csnZFJLnIyEhcvnwZwcHBuHTpEq5fv46QkBBoamriyJEjcHBwyPcaBEHAkiVLMH78eKSkpKBJkybYu3cvSpcune/7JioKGF7kS2nhZfTo0QCAIUOGAAAUCoU4aFehUODevXvK2hVRgUtKSgIAaGh8OWU8PT0xceLEdO0SExMxePBg3Lp1C3p6evlWT3x8PIYNG4aNGzcCAAYOHIi///4b2tra+bZPoqKG4UW+lBZeAgICxM/Dw8NhaGgIhULBS0mk0m7duoUNGzZg27ZtWL16NTp27AgAsLGxQaVKlWBjYwMbGxtYWVmhSpUqaNu2LYYOHZqvIeLt27fo3LkzLly4wIG5RPmI4UW+lBZeypUrB09PT2zcuBGfP3/GkSNHsHTpUujp6WHatGnK2g1Rvvv48SN27NiBDRs24MqVK+JyX19fMby0aNECT548SbfurVu3xN6Z/HD58mV07NgRr169gpGREXbt2gVnZ+d82x9RUZYaXgRBgCAI/ANBRpTWLbJy5UocOHAAf/zxB7S0tAAAnTp1wrlz5/Dnn38qazdE+SYmJgY9e/ZE2bJlMWrUKFy5cgWampro0qULfH194enpKbbN7IdYfgaXnTt3olmzZnj16hVq1KiB4OBgBheifJQaXgD2vsiN0sLLvn37MHPmTDg4OIg/2Js2bYr58+fD399fWbshyjd6enq4f/8+EhISUK9ePSxduhSvX7/Gnj174Obmlu1g8ujRI/z+++9YsmSJUupKSUnB1KlT0atXL8TFxcHV1RUXL15E9erVlbJ9IsoYw4t8Ke3PxA8fPqBUqVLplhcrVgwxMTHK2g2RUsXGxkJHR0cM3EuXLkWxYsVgZWWV620+f/4cM2bMQM2aNfHTTz/lqb7IyEj07dsXBw8eBPDlrr65c+em+aFKRPnj6z9YkpKSOCBeRpTW89K4cWNxSvJUUVFRWLRoEWxtbZW1GyKliYyMRKtWrTBp0iRxWYsWLfIUXACgZMmSAL4E+rx4/PgxmjRpgoMHD0JbWxtbtmzB/PnzGVyICsjXYSU+Pl7CSui/lBZeZsyYgbt376Jp06aIj4/HiBEj0KJFC7x69YoDdkl2Pn/+jLZt2+LChQv4559/8PbtW6Vtu0SJEgC+hJf/PuMruy5fvgwbGxvcuXMHZcuWxZkzZ9C3b1+l1UhE36ahoSH+sRAXFydxNfQ1pV02KlOmDPbs2YMLFy4gJCQESUlJqFy5Muzt7Xm7NMlKdHQ03NzccP78eRgZGeH48eMoU6aM0rafGl6Sk5MREREBIyOjHG/jjz/+wMePH2FtbY39+/dzkkciiejo6CA6OprhRWaUkirev38vDmays7ODlZUVkpOT8fHjR37DSVbi4+PRuXNnnD17FsWLF8exY8fQoEEDpe5DR0cH+vr6AHJ/6ejjx48AgHHjxjG4EElIR0cHAC8byU2ewkt0dDSGDRuGZs2a4enTpwCAvXv3omvXrtiyZQtWr14Nd3d3pXbJE+VWUlISevfujaNHj0JfXx/+/v6wtrbOl32l9r68f/8+V+unDnLPz1l6iejbUsML/xCXlzyFl+XLl+PVq1fYunUrqlSpgpiYGMyZMwd169bF0aNH4e/vD3t7eyxcuFBZ9RLl2unTp7F3715oaWlh//79sLOzy7d95XXQbmp4Se3BISJppA7aZXiRlzyNeTl69Cjmzp2Lhg0bAgACAwMRHR0NDw8P8cGMnTt3xtChQ/NeKVEetWrVCps3b0axYsXQqlWrfN3Xxo0boaGhATMzs1ytHx0dDYA9L0RSY8+LPOUpvISFheG7774TX58/fx7q6uqwt7cXl5UsWRKxsbF52Q1RnsTGxkJXVxcA4OHhUSD7tLS0zNP6vGxEJA8ML/KUp8tGpUuXxosXLwB8efbD6dOnUa9ePRQvXlxsc+3aNZQtWzZvVRLl0oIFC2BjY4M3b95IXUqOMLwQyQMH7MpTnsJLhw4dMGfOHAQEBGDu3Ll48+YNevfuLb5///59LFq0CG3atMlzoUQ5tXr1avzyyy+4ffs29u/fX6D7Dg4Oxu+//46dO3fman2GFyJ5YM+LPOUpvAwfPhx2dnaYMmUKDh48iDFjxqBdu3YAgPnz56Njx44wNzfH8OHDlVIsUXZt375d/H83efJkDBs2rED3HxwcjBkzZmDPnj05XjcxMRFJSUkAOGCXSGocsCtPeRrzoqGhgcmTJ2Py5Mnp3uvYsSPc3d1Rq1atvOyCKMcOHDiAfv36QRAEjBw5EnPmzCnwGr6eZTenvn4WGHteiKTFnhd5UtoMu/9lYWGRX5smylRAQAC6d++O5ORkeHh4YNmyZeJDFwtSXm6VTg0vampq0NLSUmpdRJQzDC/yxHn7qdBITk7Gjz/+iPj4eHTs2BHr16+X7NEUeZmk7uvbpKUIXkT0/zhgV54YXqjQUFdXh7+/P4YOHYqdO3emeZx9Qfu65yWnD2fkYF0i+WDPizwxvJDKS0hIED+vWLEiPD090zzKXgqpPS8JCQmIiorK0boML0TywQG78sTwQirt+fPnqFWrFry8vKQuJQ09PT3xh15Ox73w0QBE8sGeF3lieCGVFRoaCicnJzx+/Bi///47EhMTpS5JpFAocPr0ady9exfly5fP0brseSGSD4YXeZJuUEAuhIaGYs6cObh48SK0tbXh6uqKcePGSX6JgArep0+f4OzsjIcPH8LMzAyHDx8Wn6clF7a2trlaj+GFSD44YFeeVCa8CIKAMWPGoFixYti2bRsiIiIwZcoUqKmpYeLEiVKXRwUoKioKrq6uuHnzJsqUKYPjx4+jQoUKUpelNHwoI5F8sOdFnlTmslFISAiuX7+OefPmoXr16rC2tsaYMWPg6+srdWlUgOLi4tCxY0dcvHgRxsbGOHr0KKpVqyZ1WRk6fPgwZsyYgdOnT+doPfa8EMkHB+zKk8r0vJiammLt2rXiLaipcnonB6m2devWISAgAAYGBjh8+HCen96cn/bv3w9PT08AQIsWLbK9HgfsEskHe17kSWXCS7FixdCsWTPxdUpKCrZu3YrGjRtLWBUVtOHDh+Pp06dwc3ODjY2N1OVkKbcT1bHnhUg+GF7kSWXCy38tWLAAd+/ezdWD70i1CIKAlJQUqKurQ01NDQsWLJC6pGzJ7SMCGF6I5IMDduVJZca8fG3BggXYtGkTFixYAHNzc6nLoXw2efJk9OjRQ+V+eLDnhUj1ccyLPKlceJk1axY2bNiABQsWwMXFRepyKJ/NmzcP8+fPh7e3N44fPy51OTmS2ydL824jIvngZSN5UqnLRitWrMDOnTuxaNEitGnTRupyKJ+tWLECU6ZMAQAsXLgQbm5uEleUM3m9bMQBu0TSY3iRJ5UJL48fP8aqVaswZMgQNGzYEGFhYeJ7pqamElZG+WHz5s0YPXo0AODXX3/Fzz//LHFFOcfLRkSqj+FFnlQmvAQEBCA5ORl///03/v777zTvPXjwQKKqKD+cPXsW33//PQBgzJgx+P333yWuKHcqVKiAs2fPokSJEhAEAQqFIlvrMbwQyQcH7MqTyoSXIUOGYMiQIVKXQfksISEBffv2RXJyMnr27InFixdn+5e+3Ghra8Pe3j7H6zG8EMkHB+zKk8oN2KXCTUtLC9u2bUPr1q3xzz//QE2t6P0XZXghkg9eNpInlel5oaLD3t4eR48elboMpdixYwcePHiAXr16wcLCIlvr8G4jIvlIDS/JyclISkqChgZ/bcpB0fuzlmTp+PHjuHPnjtRlKJ2npyd+//133LhxI9vr8G4jIvlIDS8Ae1/khOGFJBcSEoJu3brBxsYGly5dkrocpcrNHUe8bEQkH6ljXgAO2pUThheSVFxcHLp164bw8HDUq1cP9erVk7okpcrNRHUML0Tyoa6uLl4qYs+LfDC8kKTGjBmDq1evomTJkti1axe0tLSkLkmpcjpRXVJSEhISEgAwvBDJBQftyg/DC0lm06ZN+Oeff6BQKLB9+3ZUrFhR6pKULqeXjWJjY8XPGV6I5IHhRX4YXkgSN2/exLBhwwAAM2bMQOvWrSWuKH/ktOcl9U4jhUKRZqAgEUmH4UV+GF5IEkuXLkVcXBxcXFwwbdo0qcvJNzntefl6vIuqTs5HVNhwll354Q3rJInVq1ejatWqGDJkSKGeiK5p06YIDAxEmTJlstWeg3WJ5Iez7MoPwwtJQkNDQ3xidGFmYmKCpk2bZrs9wwuR/PCykfwU3j95SXY+f/6MJUuWiL+gKT2GFyL5YXiRH4YXKjDLly/H2LFj0bZtW6lLKVArVqzAb7/9hvDw8G+2ZXghkh+GF/nhZSMqEBEREVi4cCEAYOjQoRJXU7CmT5+Ojx8/okePHjAyMsqyberdRnw0AJF8cMCu/LDnhQrE0qVL8enTJ9SsWRM9evSQupwClZPbpdnzQiQ/HLArPwwvlO/Cw8OxaNEiAF96IdTV1SWuqGDl5HZphhci+eFlI/lheKF8t3jxYkRERKB27dro1q2b1OUUOPa8EKk2hhf5YXihfPXx40csXrwYAPD7778X6jldMsOeFyLVxvAiPxywS/kqJiYGzs7OePz4MTp16iR1OZLIyZOlU8MLB+wSyQcH7MoPwwvlqwoVKmDPnj2IjY0tkr0uAC8bEak6DtiVH4YXKhC6urpSlyCZXr16oXnz5tl6anbqrdIML0TywctG8lM0/xSmfBcVFYVx48bh8ePHUpciOTMzMzRp0iRb4YU9L0Tyw/AiPwwvlC/Wr1+PxYsXw83NDYIgSF2OymB4IZIfhhf54WUjUrqkpCTxDqMff/wRCoVC4oqkFRkZiY0bNyIqKuqbD6PkgF0i+Ukd88IBu/LB8EJK5+3tjadPn6JkyZLo37+/1OVILjY2VgxxEydOzHKSPva8EMkPe17kh5eNSOmWLFkCABg5ciR/CQMwMTEBAAiCgE+fPmXZluGFSH4YXuSH4YWU6t69e7h48SLU1dUxfPhwqcuRBU1NTRQvXhzAtyeq491GRPLD8CI/DC+kVBs3bgQAuLm5oXTp0tIWIyPZnaiOPS9E8sPwIj8ML6RUBgYGKFGiBAYMGCB1KbKS3YnqGF6I5IcDduWH4YWU6tdff8Xr16/h7u4udSmykt3nG/FuIyL5Yc+L/PBuI1I6LS0tqUuQnexcNkpJSRF/OLLnhUg+GF7khz0vpBTh4eE4efIkUlJSpC5FliZNmoRz585leet4bGys+DnDC5F8MLzID8MLKYW3tzccHR15uSgTtWvXRpMmTVCqVKlM23z+/BkAoFAoxB+WRCQ9hhf5YXghpdi1axcAoEmTJhJXorpSLymZmJgU2SdwE8nR1wN2+bgTeeCYF8qzsLAwnDhxAgDQo0cPiauRp6dPn2L//v0wNDTE999/n2Gb1MG8qXcmEZE8pPa8pKSkICkpCZqamhJXRPzzjvLM29sbycnJaNCgAapVqyZ1ObL077//4qeffhJnH84IwwuRPH19GZeXjuSB4YXyLPWSEXtdMpd6t9HHjx8zbZMaXkxNTQukJiLKntTLRgDDi1wwvFCevHnzBqdPnwYAdO/eXeJq5Cs7t0qz54VIntTU1MQpIBhe5EElw0tCQgLatWuHoKAgqUsp8nx8fCAIAmxtbVGpUiWpy5Gt1IczxsXFiRPR/RfDC5F8cZZdeVG5Abvx8fH4+eef8fDhQ6lLIQBDhgxBvXr1kJycLHUpsmZoaAgNDQ0kJSXh48ePGc7jEhYWBoDhhUiOdHR08PnzZ/a8yIRKhZdHjx7h559/5q1qMqKurs7bo7NBoVDAxMQE7969w4cPH1ChQoV0bdjzQiRfnOtFXlTqslFwcDBsbW3FAaJEquRbg3YZXojki+FFXlSq56V3795Sl0BfmTJlCj58+IBRo0bB0tJS6nJkb+3atVAoFKhdu3aG7zO8EMkXw4u8qFR4IflISUnB+vXrERoaim7dujG8ZMO3Lq8xvBDJFwfsyotKXTYi+bhy5QpCQ0NhYGCA5s2bS12OyouNjUV0dDQAzvNCJEfseZEX9rxQrvj5+QEAnJ2dxfkPKGtXrlzB2bNnUbNmTbi4uKR5L3X+F01NTRgaGkpRHhFlgeFFXtjzQrni6+sLAHBzc5O4EtVx5MgRjB07NsMB51/fJq1QKAq6NCL6BoYXeWF4oRyLiIjA1atXAQBt2rSRuBrVkdXdRhzvQiRvqeElNjZW4koIYHihXLh8+TIEQYCZmRnKlSsndTkqI3WW3YweEcDwQiRvxsbGALJ+PhkVHJUd8/LgwQOpSyiyYmJiYGFhgfr160tdikphzwuR6kodSJ96iZekpbLhhaTj7u4Od3d3PhIgh7J6OCPDC5G8MbzICy8bUa6pq6tLXYJK+fqy0X8fccHwQiRvqeHl3bt3EldCAMML5VBiYiJ7XHIpteclKSkJUVFRad5LDS+c44VIntjzIi8ML5QjBw4cgJGREYYOHSp1KSpHV1cXvr6+OH/+vHjnQio+UZpI3hhe5IVjXihHgoODERUVxblIckGhUGQ6Lw4vGxHJW2p4ef/+PVJSUqCmxr/9pcSvPuVIUFAQAMDGxkbiSgoXhhcieUsNL8nJyQgPD5e2GGJ4oexLTk7G5cuXAQC2trYSV6Oajh07hiVLluDWrVviMkEQGF6IZE5bW1t8dAcvHUmP4YWy7e7du4iOjoaBgQFq1KghdTkqydPTE2PHjsXZs2fFZZ8/f0ZiYiKA/x/US0Tyw3Ev8sHwQtmW2utibW3N26RzKaNZdlN7XfT09KCnpydJXUT0bQwv8sHwQtl2+/ZtAEDdunUlrkR1ZTTLLi8ZEakGhhf5YHihbLO0tES7du3QtGlTqUtRWRn1vKT+IOQcL0TyxvAiH7xVmrJtwIABGDBggNRlqDT2vBCpLoYX+WDPC1EByqrnhYN1ieSN4UU+GF4oW8LDw/H27dt0z+ShnMno4YyPHz8GAFSqVEmKkogomxhe5IPhhbJl+/btKFu2LHr16iV1KSrN0tISvr6+2LFjh7js/v37AMDbz4lkrlSpUgAYXuSAY14oW1LvNGLvQN4YGxune0QAwwuRamDPi3yw54WyJTW81KlTR+JKCpeIiAi8ffsWAGBhYSFxNUSUla/DCy+hS4s9L/RNgiAwvCiRl5cXXr58id69e+PZs2cAgHLlyqFYsWISV0ZEWUkNLwkJCfj8+TPPWQkxvNA3vXnzBp8+fYKamhovbSjBtGnT8O+//6Jhw4Z4+vQpAPa6EKmC1FmwY2JiEBYWxvAiIV42om9K7XWpXr06dHR0JK5G9aXO5/L+/XuOdyFSMRz3Ig8ML/RNvGSkXAwvRKqL4UUeeNmIvsna2ho//fQTGjRoIHUphULqDz+GFyLVw/AiDwwv9E3NmzdH8+bNpS6j0EjteXn37h0ePXoEgOGFSFUwvMgDLxsRFbDU8PL06VMkJiZCV1cXFSpUkLgqIsoOhhd5YHihLIWGhuLs2bP4/Pmz1KUUGqnh5cWLFwC+3GmkpsZTkUgVMLzIA39iUpb8/PzQvHlzdOjQQepSCo1WrVrBz88PLVu2BMBLRkSqhOFFHhheKEtXr14FAA7WVaKKFSvC1dUV4eHhABheiFRJ6vON3rx5I3ElRRvDC2WJ4SX/8E4jItVTu3ZtAMCdO3cQFxcncTVFF8MLZSo5ORk3btwAwPCiTMnJyVizZg2CgoIAAPXq1ZO4IiLKrkqVKqF06dJITEzElStXpC6nyGJ4oUz9+++/iImJgb6+PqpXry51OYWGmpoaRowYgeTkZJiZmfHRAEQqRKFQwM7ODgBw4cIFiaspuhheKFOpl4zq168PdXV1iaspPBQKBbS0tAAAzZo1g0KhkLgiIsqJ1PBy8eJFiSspuhheKFMc75I/kpKSkJCQAACwsrKSuBoiyqnGjRsD+NLzIgiCxNUUTZxhlzI1YMAAVK5cGfXr15e6lELl3LlzSE5OBgCULl1a4mqIKKesra2hoaGB169f48WLF/juu++kLqnIYXihTFlaWsLS0lLqMgqd/fv3i59/+vRJwkqIKDf09PRQr149XLlyBRcuXGB4kQAvGxEVIEEQ0oSX9+/fS1gNEeXW15eOqOAxvFCGhg4din/++Uccm0HKERwcjJCQEHEANMMLkWriHUfSYnihdIKDg7FmzRoMHz4cL1++lLqcQiM5ORmjR48GALi5ucHPzw8//fSTtEURUa6khpdr165xsjoJMLxQOjNnzgQAeHh4oEqVKhJXU3j8888/uHTpEooVK4bVq1fD1dUV1apVk7osIsqFypUro2zZskhMTMTWrVulLqfIUanwEh8fjylTpsDa2hr29vZYv3691CUVOpcuXYKfnx/U1NQwZcoUqcspNN69e4fJkycDAGbPno0yZcpIXBER5YVCocCECRMAAJMmTcLHjx8lrqhoUanw8ueff+L27dvYtGkTpk+fjhUrVuDw4cNSl1VovHv3Dj179gQA9O7dm7PqKsnbt2/RoUMHhIeHw8rKCsOHD0dkZCQ2b94MT09PqcsjolwaNWoUateujQ8fPmDq1KlSl1OkqEx4iYmJgZeXF6ZOnYratWujdevWGDx4MLZt2yZ1aYVCTEwM2rdvj5CQEFSpUgWLFi2SuqRC4cqVK7CxscHFixdhZGSEdevWQUNDA+Hh4ejfvz9+/PFHTnJFpKI0NTWxcuVKAMDq1atx6tQpaQsqQlQmvNy/fx9JSUlpZiRt2LAhbty4gZSUFAkrUz2CIODZs2c4ceKEuOzq1au4fv06TExM4O/vD1NTUwkrVG2fPn3Cvn374OLiAmtra7x48QIWFhYIDg4W//+WLFkSAJCQkIDPnz9LWS4R5UGLFi3Qu3dvCIIAR0dHjB49GuHh4VKXVeipzCR1YWFhMDY2Fp8JA3z5BRAfH4/w8HCYmJjkew3x8fHYsmVLlv8xv56NNiQkBBEREeLr//6FXa9ePaipfcmPT58+zXLCstq1a0ND48u368WLF/jw4UOa97/edq1atcSv08uXL/HmzRu8evUKz549w9OnT/Hs2TPExsZCTU0N3t7eYtvZs2dDEAQ8evQIjx49yrDm/8rL+6q4LgAkJiYiJiYG0dHR4r/v3r3Dq1ev8ODBAzx8+FBsq1Ao0LVrV6xZswZGRkbicj09Pejp6SEmJga+vr4oU6YMSpQokeYJ06dOnco0mBsZGaV5bMPZs2eRmJiYYVtDQ0M0atRIfH3+/PlM747Q19eHra2t+PrixYuIiYnJsK2Ojg6aNGkivr506VKmQUxTUxPNmjUTX1+9ejXT80hNTQ0tW7YUX9+4cSPd//evOTo6ip/funULYWFhmbZt0aKFeJv63bt38fbt20zb2tvbi+fGgwcP8OrVq0zb2tnZQVdXFwDw6NEjPH/+PNO2tra20NfXBwA8efIET548ybSttbU1ihUrBgB49uwZHj9+nGlbKysrGBsbA/hy3v/777+Ztq1bt64YoN+8eYN79+5l2rZOnTooVaoUACA0NBR37tzJtG3NmjVRtmxZAF+mAbh582ambc3NzVGhQgUAXwL/tWvXMm1btWpVmJmZAQAiIyNx+fLlTNtWrlwZlStXBgBERUUhODg407bfffedOGg+NjY2y9uey5cvLz5ENSEhAYGBgeJ7PXr0wNu3b3HixAmsWLECq1atgpWVFWxsbFCmTBl8+vQJOjo60NDQSPOhUChgZGSESpUqidu6fv16pjUUK1YMVatWFZ+HltEf76nvGRgYpLn0f+vWLSQlJWX4LDVdXV3UqFFDfH3nzp1Mp8nQ0dFBzZo103ydJSGoiH379gktW7ZMs+z58+eCubm58ObNmzTLHR0dBUdHR6XX4OHhIQDgBz+++VGlShVhwoQJQkhISKb/n8zMzNKs4+rqmuZ9XV3dTLffokWLNG1LliyZadtGjRplud+vP2rVqpWmbc2aNTNta2ZmlqattbV1pm1NTU3TtG3RokWmbfX09NK0bdu2bZZf66917do1y7bR0dFi2/79+2fZ9t27d2Lb4cOHZ9n2yZMnYtvx48dn2fbOnTti2+nTp2fZNjg4WGw7f/78LNuePHlSbLtixYos2/r5+YltN2zYkGXb3bt3i2137dqVZdsNGzaIbX19fbNsu3LlSrHtyZMns2z7559/im2DgoKybDt9+nSx7e3bt7NsO2HCBLHtkydPsmw7YsQIsW1oaGiWbYvCh5qamvDq1StB2bL7+1tlel60tbXTJcHU1zo6OgVSg6urKw4dOoSoqKgM31coFDA3Nxdfv3r1CtHR0Zluz9zcXEzBr1+/zvLyQfXq1cWelzdv3mTZ+1O9enVoamoC+DJY9NOnT9DW1oaOjg50dXWhq6sLbW1tsdcntfbMfOupx3l5XxXX1dDQgL6+vth7oqenh5IlS6J8+fKoVKkSGjZsiBIlSmS5feDLHQqrVq0Se3pS/7JMVbt27Ux7SP77F0/NmjUz7bmrWrVqmtcWFhYwNDTMVtvq1atn+kTxcuXKpVs3s3r/2zNauXLlTHtTUnswUpmZmaFOnToZtv2v7777Lsu2X39fK1SokGXbr4+7XLlyWbZNPd8AoEyZMlm21dbWFj8vVapUlm2//lqULFkyy7apvTnAl693Vm2//v4bGRll2Ta15wcAihcvnmXbr3sXDQ0Ns2z79f8JfX39LNum9hIBX74mWbVN7SUCvnyts2r79V1/mpqaWbb9+v+7urp6lm0dHBzQpEkT3Lx5E69fv8bevXshCEK6D+DL16l8+fLiuvfv3890u/r6+mnaPnz4ME1P8def6+npiT1bAPD48WPxmWr/paOjg4oVK4qvnz59mmlPrpaWFszMzFC9evU035eCphAE1RgtePXqVfTt2xc3b94Uf4lfvHgRQ4cOxbVr19L8Im7VqhUAICAgQJJaiYiIKOey+/tbZQbs1qxZExoaGmmuB165cgWWlpZpggsREREVbirzW19XVxcdO3bEjBkzcPPmTRw/fhzr169Hv379pC6NiIiICpDKjHkBgMmTJ2PGjBno378/DAwMMHr0aDg7O0tdFhERERUglQovurq6mD9/PubPny91KURERCQRlblsRERERAQwvBAREZGKYXghIiIilcLwQkRERCqF4YWIiIhUCsMLERERqRSGFyIiIlIpDC9ERESkUhheiIiISKUwvBAREZFKUanHA2TXu3fvkJycLD5am4iIiOTvzZs3UFdX/2a7Qtnzoq2tDQ2NQpnLiIiICi0NDQ1oa2t/s51CEAShAOohIiIiUopC2fNCREREhRfDCxEREakUhhciIiJSKQwvMhMfH48pU6bA2toa9vb2WL9+faZtT506hQ4dOsDKygru7u4ICAhI8761tTUsLCzSfERHR+f3IeRITo53+PDh6Y7n5MmT4vsbN25Es2bNYGVlhSlTpiA2NrYgDiFHsnu8Hh4e6Y7VwsICkydPBgBERESke8/W1rYgDyXbEhIS0K5dOwQFBWXa5u7du+jWrRvq1auHLl264Pbt22ne9/X1hZOTE+rVq4eRI0fi48eP+V12rmXneAvDuQtk71gLw3mb6lvHWxjO29DQUIwZMwY2NjZo1qwZ5s2bh/j4+AzbSnreCiQrM2fOFNzd3YXbt28LR48eFaysrAR/f/907e7duyfUrl1b2LRpk/D06VNh69atQu3atYV79+4JgiAIb9++FczNzYXnz58L7969Ez9SUlIK+pCylN3jFQRBaN26tbB///40xxMfHy8IgiAcPnxYaNiwoXDixAnhxo0bgqurq/D7778X5KFkS3aP99OnT2mO89ixY0Lt2rWFmzdvCoIgCJcvXxZsbGzStHn//n1BH843xcXFCSNHjhTMzc2FixcvZtgmOjpaaNq0qfDHH38Ijx49EmbNmiU0adJEiI6OFgRBEG7cuCHUrVtX2Ldvn3Dv3j2hb9++wpAhQwryMLItO8dbWM7d7ByrIBSO81YQsne8qn7epqSkCN27dxcGDx4s/Pvvv8KlS5eE1q1bC3/88Ue6tlKftwwvMhIdHS1YWlqmOTFWrlwp9O3bN13bBQsWCIMGDUqz7PvvvxcWLVokCIIgnDt3TmjatGn+FpxHOTne+Ph4oWbNmkJISEiG2+rdu7ewbNky8fWlS5eEunXrCjExMcovPJdycrxfS0pKElxdXYXFixeLy3bv3i306NEjv0pViocPHwrt27cX3N3ds/yB7+XlJTg6Ooq/nFNSUoTWrVsL3t7egiAIwoQJE4SJEyeK7V+/fi1YWFgIz58/z/+DyIHsHm9hOHeze6yF4bwVhOwf79dU8bx99OiRYG5uLoSFhYnLDh48KNjb26drK/V5y8tGMnL//n0kJSXByspKXNawYUPcuHEDKSkpadp26tQJ48ePT7eNz58/AwAePXqEypUr52/BeZST4w0JCYFCoUDFihXTbSc5ORm3bt2CtbW1uKx+/fpITEzE/fv38+8Acignx/u1vXv3IiIiAj/88IO47NGjR6hUqVJ+lptnwcHBsLW1xa5du7Jsd+PGDTRs2BAKhQIAoFAo0KBBA1y/fl18/+vvbdmyZVGuXDncuHEj32rPjeweb2E4d7N7rIXhvAWyf7xfU8Xz1tTUFGvXrkXJkiXTLI+KikrXVurzljO5yUhYWBiMjY2hpaUlLitZsiTi4+MRHh4OExMTcXnVqlXTrPvw4UNcuHABPXv2BAA8fvwYsbGx8PDwwJMnT1CzZk1MmTJFVj8Uc3K8ISEhMDAwwC+//ILg4GCUKVMGo0ePRosWLRAZGYn4+HiUKlVKbK+hoQEjIyO8ffu2QI8pKzk53lSCIGDt2rXo168f9PX1xeWPHz9GUlISunbtitDQUFhbW2Py5MlpvgZS6927d7bahYWFoVq1ammWlShRAg8fPgTwZcbs/x5XiRIlZPW9BbJ/vIXh3M3usRaG8xbI/vGmUtXztlixYmjWrJn4OiUlBVu3bkXjxo3TtZX6vGXPi4zExsam+cUGQHydkJCQ6XofP37E6NGj0aBBA/GRCCEhIYiIiMDw4cOxatUq6OjoYMCAARkmaKnk5HhDQkIQFxcHe3t7rF27Fi1atMDw4cNx69YtxMXFpVn3621l9XUraLn5/gYFBeHt27fo3r17muUhISGIiorC5MmTsXjxYrx79w7Dhg1DcnJy/hSfjzL7uqR+TeLi4mT/vc0tVT13s6swnLe5UVjO2wULFuDu3bsYO3ZsuvekPm/Z8yIj2tra6b6xqa91dHQyXOf9+/cYOHAgBEHAsmXLoKb2JY+uW7cOiYmJYupfuHAhWrRogZMnT8Ld3T0fjyL7cnK8I0aMgIeHB4oXLw4AqFGjBu7cuYPdu3eLJ1ZG29LV1c2v8nMsN9/fI0eOoHnz5jAyMkqz3M/PDwqFQlxv2bJlsLe3x40bN9CgQQPlF5+PMvu6pB5bZu/L6XubG6p87mZXYThvc6MwnLcLFizApk2bsHjxYpibm6d7X+rzlj0vMlK6dGl8+vQJSUlJ4rKwsDDo6OigWLFi6dqHhoaiT58+SEhIwObNm9NcdtDS0krTXamtrY0KFSogNDQ0fw8iB3JyvGpqauIPwFRVqlRBaGgojIyMoK2tjffv34vvJSUlITw8HKampvl7EDmQ0+8vAJw9ezbDB4zq6uqmCTwlSpSAkZGRrL6/2VW6dOk03zvgyy/21C7nzN6X0/c2p1T93M2uwnDe5oaqn7ezZs3Chg0bsGDBAri4uGTYRurzluFFRmrWrAkNDQ1xwBMAXLlyBZaWluJfZaliYmIwePBgqKmpYevWrShdurT4niAIcHJywt69e9O0f/bsGapUqZLvx5FdOTneSZMmiXMlpLp//z6qVKkCNTU1WFpa4sqVK+J7169fh4aGBmrUqJGvx5ATOTle4MslhRcvXqBhw4ZplkdFRaFRo0a4ePGiuCw0NBSfPn2S1fc3u+rVq4dr165B+N9j1gRBwNWrV1GvXj3x/a+/t2/evMGbN2/E91VNYTh3s6swnLc5pern7YoVK7Bz504sWrQIbm5umbaT+rxleJERXV1ddOzYETNmzMDNmzdx/PhxrF+/Hv369QPw5a/01OvEq1evxvPnzzF//nzxvbCwMHz+/BkKhQItW7bE8uXLERQUhIcPH+KXX35BmTJl0KJFC8mO779ycryOjo44ePAgfHx88OzZM6xYsQJXrlxB3759AXwZULdu3TocP34cN2/exIwZM9C9e3dZdT/n5HiBLwM5U//q/pqBgQEaNmyIefPm4ebNm7hz5w7Gjh2LZs2awcLCokCPKbe+PtY2bdogMjISc+bMwaNHjzBnzhzExsaibdu2AIBevXph//798PLywv379/HLL7+gZcuWGd7BIleF7dzNSmE7b7+lMJ23jx8/xqpVq/DDDz+gYcOG4v/NsLAwADI7b5VywzUpTUxMjPDLL78I9evXF+zt7YUNGzaI75mbm4v30Lu4uAjm5ubpPlLvq4+LixPmzZsnNG3aVKhXr54wdOhQ4fXr11IcUpaye7yC8GWOBGdnZ6FOnTpCp06dhODg4DTbWr16tWBnZyc0bNhQmDx5shAXF1dQh5FtOTlePz+/TOf7CA8PFyZNmiTY2toKVlZWwvjx44Xw8PD8Lj/X/js3xn+P9caNG0LHjh0FS0tLoWvXrsKdO3fSrO/t7S20aNFCqF+/vjBy5Ejh48ePBVZ7bmR1vIXl3E31re9tYThvv/at41Xl83b16tUZ/t80NzcXBEFe561CEP7X50NERESkAnjZiIiIiFQKwwsRERGpFIYXIiIiUikML0RERKRSGF6IiIhIpTC8EBERkUpheCEiIiKVwvBCpCKGDBmSbqp1X19fWFhYYPny5WmWr1q1Ch06dMjXeiwsLBAUFJSv+wC+PMxt9+7d4msPD490x5sVCwsLWFhYYNCgQTne96RJkzBp0qQcr5dfrl+/Lh5PTr4GRIUNwwuRirC2tsatW7fSLAsKCkKpUqXShYjr16/DxsamIMvLN35+fvD09MzTNpYvX45FixbleL2pU6di6tSpedq3MtWuXRuBgYGwsrKSuhQiSTG8EKmIhg0b4vHjx4iOjhaXBQUFYdCgQbh+/Xqa56vcuHGj0IQXZUwCXrx48XRPN84OQ0NDGBoa5nn/yqKpqQlTU1NoampKXQqRpBheiFSEpaUlNDU1cefOHQDA27dv8fr1a3Tr1g2Ghoa4evUqAODJkyeIiIiAtbU1BEGAp6cnHB0dUadOHdjb22PFihUAgDNnzqBevXqIjY0V9xEYGIgGDRogLi4OgiBg5cqVsLe3h7W1NYYNG4bXr19nWFtCQgJmz54NW1tb2NraYvz48QgPDwcAvHz5EhYWFjh69CicnJxgaWmJoUOHiu+n7tfd3R1169bF4MGDMWvWLEyaNAlBQUGYPHkyXr16BQsLC7x8+RLAl6fxDh48GJaWlnBxccH58+ez9TVMreXUqVNwdHSElZUVZs+ejX///RedO3dG/fr1MXToUERFRQFIf9lo//79aNOmDerVq4eePXvi7t27adq1b98ednZ2ePr0KSIiIvDrr7+iSZMmaNiwISZMmICIiAgAX0Kno6Mjtm/fjmbNmqF+/fqYMGECEhISAACRkZEYPXo0rK2t0ahRI4wfP16siYgYXohUhpaWFurVq4ebN28CAC5evIg6depAX18fjRo1Ei8dXb9+HdWrV4exsTF8fHywadMmzJkzB4cPH8bIkSOxfPly3LlzB02aNIGuri7OnDkj7uPo0aNwdHSEjo4Otm7dioMHD+Kvv/7Crl27UKJECXz//fdITExMV9uiRYtw+/Zt/PPPP9i8eTOioqLw448/pmnj6emJRYsWYevWrbh16xY2bNgAAHjx4gWGDx+Otm3bwsfHB5aWlti2bRsAwMrKClOmTEGZMmUQGBiIsmXLAgB8fHzg6uoKPz8/1KlTB7/88kuOemjWrFmDVatWYdasWdiyZQtGjRqFn3/+GevWrcP169exZ8+edOucPXsWU6dORf/+/XHgwAHUqVMHQ4cOFQPH/v378dNPP2H16tWoVKkSRo0ahXv37sHT0xMbNmzA48eP0wShd+/e4ciRI1i7di2WL1+Oo0ePwsfHBwCwbNkyhIWFYceOHdi8eTPu37+PVatWZfv4iAo7DakLIKLss7a2FsNLUFAQbG1tAQA2Njbw9fUFkHa8S9myZTFv3jzY2dkB+PKY+pUrV+Lhw4eoXbs2nJ2dcfToUbi4uCA5ORnHjx/H7NmzAQBr167F9OnTxX3MnDkT9vb2OHv2LBwdHcWaYmNjsXXrVnh7e8PCwgIA8Oeff8LW1hYPHjyAvr4+AGDMmDGoW7cuAMDd3V0cv+Pl5YW6detixIgRAIAff/xR7EnR0tKCoaEh1NXVYWpqKu7TxcUFnTt3BgD88MMP8PX1xYcPH1CyZMlsfR1HjBiBGjVqoEaNGpg7dy7c3NzQtGlTAICdnR1CQkLSrbNr1y60a9cOvXr1AgD88ssv0NTUFHtTLC0txa/L/fv3ERwcjMOHD6Ny5coAgAULFsDV1VXcdmJiIqZNm4bq1avDwsICzZo1w61bt9C9e3e8evUK+vr6qFChAnR1dbF06dJsHRdRUcGeFyIV8vWg3aCgIDGk2NjY4Pbt20hISMD169fRqFEjAEDjxo1hbGyMv/76CyNGjICDgwPCwsKQkpICAHBzc8OpU6eQkJCAy5cvIzExEfb29oiOjsbbt28xduxYWFlZwcrKCtbW1ggPD8fTp0/T1PTixQskJiaiZ8+eYtsWLVogJSUlTVszMzPxcwMDA7EH58GDB7C0tEyzzfr162f5dahYsWKabQFAfHx8Nr+KadfX0dFB+fLl07xO7U352pMnT1C7dm3xtZaWFiZOnCiGqq+3ERISgmLFionBBQCqVq2K4sWLpwlG//2aJCUlAQD69euHq1evws7ODsOHD8etW7dQqVKlbB8fUWHHnhciFWJlZYV3797h1q1bePfuHRo0aAAAqF69OgwNDXHp0iU8evRIDDVeXl6YO3cuunXrBmdnZ0ycOBH9+vUTt9eoUSPo6enh/PnzOHv2LJycnKClpSUO/l26dGmaX8AA0g18TU5OBgBs374denp6ad4rUaKEOLYls0Gm6urq6S75fOsSkLq6erplObls9N/11dS+/XechkbWPy61tbXFz7W0tDJsk5ycLH69MmqXegx2dnY4ffo0AgICcOrUKfz2228IDAzEwoULv1knUVHAnhciFaKnp4eaNWti165dsLS0hK6uLgBAoVCgUaNG2Lt3LypVqgQTExMAwI4dOzBy5EhMmTIFHTt2hLGxMT58+CD+klRTU0ObNm1w6tQpBAQEwM3NDQBQrFgxlChRAmFhYTAzM4OZmRnKli2LBQsW4MmTJ2lqqlixItTV1REeHi62NTAwwLx58/Dhw4dvHlP16tXFQcipvn6tUChy/wVTIjMzM9y/f198nZycDEdHR1y5ciVd28qVKyMyMjJNL8ujR48QFRWVLgxmZOPGjbhz5w46deqEpUuXYt68eTh69KhyDoSoEGB4IVIxjRo1gp+fX7pboW1sbBAQECBeMgIAY2NjXLhwAU+ePMHt27cxduxYJCYmprks4ubmhv379yM+Ph6NGzcWlw8YMABLlizBiRMn8PTpU0ybNg1Xr15FlSpV0uzXwMAA3bp1w4wZMxAUFIRHjx7hl19+wbNnz1ChQoVvHk/37t1x/fp1rFmzBk+ePIGnpycuX74shhZdXV1ERETg6dOn4mUVKXh4eODAgQPYt28fnj17hnnz5kEQhDSXklJVrVoVzZs3x8SJE3Hz5k3cvHkTEydORKNGjWBubv7Nfb19+xYzZ87E9evX8fTpUxw5cgS1atXKj8MiUkkML0QqpmHDhoiJiREH0qaysbFBbGxsmlAzZcoUREVFoUOHDhg9ejQsLCzQunVr3Lt3T2xTv359GBsbw9nZOc2lkUGDBqFr16747bff0LFjR7x+/Rrr1q3LcL6USZMmwc7ODmPGjEH37t2hoaGBNWvWZHh557/Kly+PZcuWwdvbG+7u7rh27RpatWolXmZq3LgxzMzM4O7unqbugtaoUSNMnz4dK1euRPv27cU7iXR0dDJsP3/+fFSsWBEDBgzAoEGDUL16daxcuTJb+/rxxx/RoEEDDB8+HB06dEBMTAwWLFigzMMhUmkKQRkzQBER5dK///6LpKSkND0LQ4YMgaWlJUaPHp3n7VtYWGDz5s3pwp4q8/DwgI2NjVK+PkSqiD0vRCSp58+fY+DAgTh37hxevXoFLy8vXLhwAa1bt1baPiIiIsRbmlVZYmIiwsLCMpxrh6go4d1GRCQpJycnPHz4EFOnTsWHDx9QuXJlLF68GDVq1FDaPkaPHg17e3usW7dOaduUwp07d9CjRw8AEOelISqKeNmIiIiIVAovGxEREZFKYXghIiIilcLwQkRERCqF4YWIiIhUCsMLERERqRSGFyIiIlIpDC9ERESkUhheiIiISKUwvBAREZFK+T8iS53BFb6hoQAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<Figure size 640x480 with 1 Axes>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"# Plot the sensitivities of the instruments\n", | |
"w = np.arange(0.1, 2, 0.01) * u.micron\n", | |
"fig, ax = plt.subplots()\n", | |
"ax.plot(w, p.NIRDA.sensitivity(w), label='NIRDA', c='k', ls='-')\n", | |
"ax.plot(w, p.VISDA.sensitivity(w), label='VISDA', c='k', ls='--')\n", | |
"ax.set(xlabel='Wavelength [microns]', ylabel=f\"Sensitivity [{p.NIRDA.sensitivity(w).unit.to_string('latex')}]\")\n", | |
"plt.legend()" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"id": "72d47073-adee-492b-94bf-aed9b1d7dc6c", | |
"metadata": {}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 3 (ipykernel)", | |
"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.9.13" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 5 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment