Skip to content

Instantly share code, notes, and snippets.

@dokeeffe
Created December 14, 2016 17:04
Show Gist options
  • Save dokeeffe/0d9cfd39020ad6ffd2603d650ba43f81 to your computer and use it in GitHub Desktop.
Save dokeeffe/0d9cfd39020ad6ffd2603d650ba43f81 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"import astropy.units as u\n",
"from astropy.coordinates import EarthLocation\n",
"from astropy.time import Time\n",
"from pytz import timezone\n",
"from astroplan import Observer\n",
"from astropy.coordinates import SkyCoord\n",
"from astroplan import FixedTarget\n",
"from astroplan.plots import plot_airmass \n",
"import matplotlib.pyplot as plt \n",
"%matplotlib inline\n",
"\n",
"longitude = '-8d21m00.000s'\n",
"latitude = '+52d15m00.600s'\n",
"elevation = 100 * u.m\n",
"location = EarthLocation.from_geodetic(longitude, latitude, elevation)\n",
"\n",
"observer = Observer(name='bobs',\n",
" location=location,\n",
" pressure=0.615 * u.bar,\n",
" relative_humidity=0.11,\n",
" temperature=0 * u.deg_C,\n",
" timezone=timezone('UTC'),\n",
" description=\"BOBS\")"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
" \n",
"coordinates1 = SkyCoord('18h57m11.16s', '+51d16m08.9s', frame='icrs')\n",
"coordinates2 = SkyCoord('00h38m04.07s', '+34d42m42.2s', frame='icrs')\n",
"coordinates3 = SkyCoord('03h09m28.54s', '+30d40m26.0s', frame='icrs')\n",
"\n",
"\n",
"\n",
"exo1 = FixedTarget(name='HAT-P-37 b', coord=coordinates1)\n",
"exo2 = FixedTarget(name='HAT-P-19 b', coord=coordinates2)\n",
"exo3 = FixedTarget(name='WASP-11 b', coord=coordinates3)"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhUAAAGLCAYAAABwYGDQAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzs3Xe8FNX9//HX595LVS5I70VRBEWaBsEK9v5LrCiKxood\nYozxG42axJYoFtQYW1QQo0YlsXdQ1KggiooUBRHpXWnCvef3x5nxLtfb2b0zs/t+Ph772NnZmeVz\nmL2znzltzDmHiIiIyNbKizoAERERyQ5KKkRERCQtlFSIiIhIWiipEBERkbRQUiEiIiJpoaRCRERE\n0kJJhYiIiKSFkgoRERFJCyUVIiIikhZKKkRERCQtYpFUmNk+ZvYfM/vOzIrN7Ogq7LO/mU02sw1m\nNtPMhtVGrCIiIlK2WCQVwDbAVOB8oNKbkZhZZ+A54HWgF3A7cL+ZHZS5EEVERKQiFrcbiplZMfD/\nnHP/qWCbm4DDnHO7pawbBzR2zh1eC2GKiIhIKXGpqaiuPYHXSq17GRgQQSwiIiICFEQdQA21BhaX\nWrcYKDSzes65jaV3MLNmwCHAXGBDxiMUERHJHvWBzsDLzrnl5W2U1KSiJg4BxkYdhIiISIKdAjxW\n3ptJTSoWAa1KrWsFrCmrliIwF2DMmDF07949g6FFb8SIEYwaNSrqMDJO5Yy3pUuXMvrh0dTrWo/6\njepXuv1bD73F/mfsn/nA0mDD9xvYOHsjFw67kBYtWlRr36Qez+pSObPL9OnTGTp0KAS/peVJalLx\nHnBYqXUHB+vLswGge/fu9O3bN1NxxULjxo2zvoygcsbdwoULadayGc16NaNRs0aVbv/hsx/SdY+u\ntRDZ1vt++fcsX7Oc3XbbjTZt2lRr36Qez+pSObNWhd0HYtFR08y2MbNeZtY7WLV98LpD8P4NZvZw\nyi5/D7a5ycy6mdn5wHHArbUcuoiIiARikVQAuwMfA5Px81TcAkwBrg3ebw10CDd2zs0FjgAOxM9v\nMQI40zlXekSIiIiI1JJYNH845yZQQYLjnDujjHUTgX6ZjEtERESqLhZJhaTXkCFDog6hVqic2WXX\nwbtW+H5xkbFpYx02bajDpo11yMsvpk69zRTU3USdepuxvHhN5FeeXDmeKmduUlKRhXLlS65yZpcu\nffZk9getWfpNC1YvacyaJY1Zs6SQ1Usbs35NA4o2VXy6Kqi7iUbNv6ew+RoKW6yhUfPvadJmJa22\nX0yrHRZTr+GPtVSSiuXK8VQ5c5OSChGpdZs2FvDtZx2Z83FnFs1qw6LZrflhReUjREIF9TZRXJRH\n8eb8n9Zt/rEOKxc0ZeWCpmXus13bFbTeYRHtesynS585tO66iLz8ZNRuzJs3j2XLlkUdhmSx5s2b\n07Fjx63+HCUVIpJxzsHir1oz+4OufD15e+ZN6/jzmgdzNO+wjJZdltCk9SoKW66mccs1FLZYzTZN\n1lGn/o/Uqb+JgrqbMfO7FBcZm38sYPOPddi4ri7fL2vEmqWFrFlWyJoljVkxvymLvmrN98sKf0o4\npr/dA4D6jdbTpfccuvSdQ7e9ZlDYYk0t/69Uzbx58+jevTvr1q2LOhTJYg0bNmT69OlbnVgoqRCR\njFm1og0z/jWIWe/3Zfm3zbd4r1HzNWzf72va95hP666LaNllMXUbbKrW5+flO+o22ETdBpto2Hgd\n27VZVeZ261Y3ZNHsViya3YZvPu3E3Kmd2fB9A6a/3YPpb/fghduPoGPPb9hl0Of02O8Ltm36Q43L\nnG7Lli1j3bp1OTFxn0QjnNhq2bJlSipEJF5WrYKHHoL772/OF1/c8NP6grqb2GGPr9i+31ds3+9r\nmnVY/lONQ6Y1bLyO7fvNYft+cxh44rsUF+WxYEZb5kzpwuwPujJvWqefHi+NPpTOveey+9Ef0W2v\nL8kvKK6dICuRCxP3SfIpqRCRtPjiC7jzTnj0UVi7FqAOeXmb6dxnBr0O+ZJue82ITWfJvPxi2veY\nT/se89ln6NusWVrI52/14PM3d+W76e2ZM2V75kzZnkbN17D70R/R94jJbNt0bdRhi8SekgoR2Spv\nvgnXXw+vpUw9t8sucOqpq5k650raHlS/StN0R6mwxRoGHP8+A45/n5ULm/Dxi32Y8lw/vl9WyJsP\nDmbCI/vR84Bp7HvqRJq2WxF1uCKxpaRCRGrk3XfhqqvgjTf867w8OOYYuOgi2H9/WLRoHdNvWou/\nY3JybNdmFYN//Sb7Dp3IFxN68OGzv2D+Fx345OXefPrqbvQ6ZCr9jnwRKPfuzyI5S0mFiFTL5Mk+\nmXjxRf+6Th045xy47DLo3DnS0NKqoG4Rux00jd0OmsZ309sx4eH9mPW/nZj6Yl8+faUXXbtPZOFp\neVTzfmIiWS0u9/4QkZhbtgzOOgt2390nFPn5/vWsWTB6dHYlFKW16/4dJ9/4GGfedT877D6b4qJ8\nZn42iH33bcFf/wo/xqOriOSwGTNmkJeXx9133x1pHEoqRKRCxcXwj39At27wwAN+3SmnwJdfwn33\nQadO0cZXm9r3mM/Qv47hxD+NpkWbWaxdm8fll0OvXvD661FHl0wPP/wweXl5TJkypcz3999/f3bb\nbbcy3ysuLqZt27bk5eXx8ssvb/HeV199RV5eXqWP/Px8FixYUG5899577xbbN2jQgJ133pkRI0ZU\naUKyNWvWcNVVV3HIIYew3XbbkZeXxxNPPFHu9qNGjaJ79+7Ur1+fDh06cMUVV7BhQ4V3G48VNX+I\nSLmmTIHzzoMPP/Svd9sN7r4b9tor2rii1q77HA47/nr6dBnF9dc34csv4cAD4YQT4I47oFWrqCNM\nFqtgbHFF773xxhssWrSILl26MHbsWA455JCf3mvdujVjxozZYvubb76ZJUuWcMstt+BcyWyqTZuW\nPQtragw33ngj7dq1Y/369UycOJE77riDl19+mU8++YQ6deqUu+/ChQv5y1/+QpcuXejTpw8TJkwo\nd9tLLrmEO++8k5NPPpkRI0Ywbdo0brnlFmbMmMEzzzxTYYxxoaRCRH6mqAhuvBGuuQY2b4bCQvjT\nn+D886FAZw0AzBwnnrie009vwtVXw113wRNP+NEw998PRx8ddYTZb8yYMfTr149hw4Zx5ZVXsn79\neho0aADANttsw8knn7zF9o8++ijr1q2r0f06jjjiCHr08LOxnnnmmTRq1Ih77rmHF154gWOOOabc\n/bp06cLixYtp0aIFkyZNYp999ilzu2+++Ya77rqLc889l3vuueen9Z07d+byyy/n9ddf54ADDqh2\n3LVNzR8isoU5c2C//eAPf/AJxXHH+aaOiy9WQlGWJk187cTkyb4mZ+lSPwrmrLPg+++jji57bdiw\ngWeeeYYhQ4Zw/PHHs27dOsaPH19r//7gwYNxzjFnzpwKt6tbty4tWrSo9PPeffddnHOceOKJW6w/\n6aSTcM7x+OOPVzm2m2++mY4dO9KwYUMOPPBAZsyYUeV9t5aSChEB/P05/vlP/8M4aRI0agSPPOKv\nvjXCoXK9e8MHH8Dll4OZ73/Sq5f/v5TKrV69muXLl2/xWLZsGZs2lT11+/jx41m7di0nnXQSrVq1\nYv/992fs2LG1Fu/s2bMBaNasWVo+b+PGjQA/1bSEGjZsCMDkyZOr9Dn33nsv999/P5dccglXXHEF\nkydPZvDgwaxcuTItcVZG1x0iwvr1flho2AS9zz4+ocjmER2ZUK8e3HQTHHEEnHZaSa3PLbf4mp5M\nT0vuHNTGfccaNkxvWZxzFVbt77rrrj9bN3bsWAYOHEjbtm0Bf0V/wQUXsHz58rT90KdatWoVy5cv\nZ/369UyYMIEbb7yRwsJCDjvssLR8frdu3XDOMWnSJPr37//T+okTJwLw3XffVelzvvnmG2bPnk3z\n5v5eO4MHD2bffffllltu4c9//nNaYq2IkgqRHDdvHvzyl75TZn6+7ztx+eV+WWpm333h009h+HB4\n7DG49FLfPHLvvVDqQjSt1q2DbbfN3OeHfvgBttkmfZ9nZtx9993suOOOP3tv5MiRFBdvef+VFStW\n8PLLL3P77bf/tO7YY4/lggsu4IknnmD48OHpCw6f9Oy9995bxLvDDjtw3333/fTjvbUGDBhA7969\n+dOf/kTLli3Zd999mTZtGhdddBF169Zl/fr1Vfqc448/fouY9t57b3r16sULL7ygpEJEMmviRN9n\nYulSaN4cnnzSz4YpW6+w0Nf8/OIX8Jvf+HuifPYZPPNMbg3Drao99tijzBumbbfddixfvuXspY8/\n/jibN2+md+/efPXVV4D/4e/fvz9jx46tdlJRXFzM0qVLt1jXrFkzCoJORGbG/fffT6dOnSgoKKBN\nmzZbJECV7V9V48eP54QTTmDYsGE456hTpw5XXHEF48ePZ8mSJVX6jK5du/5s3U477cSrr75arVhq\nSkmFSI666y5/Bb15M/Tpox+7TDCDSy7x/VROOAE+/thPHvbUU75ZJN0aNvS1CJkWNPNH5rHHHgNg\n4MCBW6wPh5/OnTuXztVou5s1axbdu3fHzHDOYWa89957/OIXv/hpm/79+/80+qMm+1dFhw4deO+9\n95g9ezaLFy+mW7duNG/enNGjR9OzZ89qfVZUlFSI5JjiYvjd7+Bvf/OvhwzxQyCj/qHIZoMG+eaP\nsJnp4INh7FhfS5ROZultloijuXPn8u6773LxxRez7777bvFecXExQ4cO5bHHHuPKK6+s8md27NiR\n11LviAflJhCZ2L+0rl27/lTjMGXKFFauXMlBBx1UpX1nzZr1s3UzZ86sVpK1NZRUiOSQTZvg7LPh\n4Yf96xtvLBmtIJnVsSO88w6ceir8+9++5uLuu/3kYlJ1Y8aMwcz47W9/S7t27X72/n333cfYsWOr\nlVQ0aNCAwYMH1zimrd2/PEVFRfzud7+jsLCQs846q0r7PPXUU1x//fW0bNkSgLfffptPPvmkWv8f\nW0NJhUiOWL/e/5A995zvhPnAAzBsWNRR5ZYGDeBf/4ILLvCdNocPh8WL4eqrczuxS53dsjJjx46l\nd+/eZSYUAEcffTQXXXQRU6dOpXfv3ukKcavcdtttrF27lm+++QaAp59++qcahZEjR/40jHT48OGY\nGb169WLDhg2MGTOGqVOn8vjjj9OqitO0du7cmb322ovzzjuPH374gdtuu402bdowcuTIzBSuFCUV\nIjlg1So46ih/pVy/vp974qijoo4qN+Xnwz33+Km8r7vOz1q6ZImfQCtXR9xUNBV36vsff/wxM2fO\n5Oqrry5326OOOoqLL76YMWPGlJlUVPZvZcJNN930U0dLM+PJJ5/kySefBODss8/+Kano168fo0eP\n5tFHH6WgoIA999yTCRMm/KzvSHnMjLPPPpt169Zxxx13sGzZMgYMGMDo0aMrnYo8XZRUiGS5FSvg\ngANg6lRo3Bj++18/D4VExwyuvRZatoSLLvLNIOvX+74teTk2JeGwYcMYVkGV2ZtvvvnTcp8+fSgq\nKqrw8zp27FjuNi+++GK14zv33HM599xzq71fqoULF1Zpu7POOqvKzRyldevWbYty//a3v63R52yt\nHPv6iuSW1avh0EN9QtGyJUyYoIQiTi64wM9jkZcHDz3kX1ejJUAkdpRUiGSptWv9zI4ffgjNmvlb\nc/fqFXVUUtpJJ/k5LMzg73+HESOUWEhyKakQyULr1/u7ZE6a5Js8XnkFypjpWGLi5JN9x1mA22+H\nK65QYiHJpKRCJMv8+KOf/+CNN/yUzS+9BGVMVCgxc8YZvgMnwM03+z4XIkmjjpoiWaS42A8TfeEF\nP3zxuedgzz2jjkqq6rzzYMMG3wRy7bX+7rB77BF1VCJVp5oKkSxyzTXw+ONQUADPPpuZqaAlsy69\ntKSW4oIL4N13o41HpDqUVIhkiUce8XcYBfjHP/xU0JJMV13lb51eVOSnVBdJCiUVIllg4kQIh7f/\n/ve+fV6Sywzuu8/fMXbduqijEak69akQSbhZs/yNqjZtguOPhz//OeqIJB3q1vX3COnTB+bNg+nT\np0cdkmSpdH63lFSIJNjKlX4uihUroH9/f6OwXJuRMZs1bQpjxjRn0KCGDB06NOpwJIs1bNiQ5s2b\nb/XnKKkQSajiYt/uPmsWdOoE48f7ER+SXfbZpyNPPDGdE09cxubNviPnqadGHZVkm+bNm9OxY8et\n/hwlFSIJdcstfshovXp+pEcVb2IoCfSrX3Xkjjs6cv75cOedcOyxsPfeUUcl8nOqKBVJoHfe8R0y\nwf/IxOQOz5JB550HQ4b4ESEnnujvbCoSN0oqRBJm6VL/o1JUBEOHloz6kOxm5ocKd+8OCxbAKaf4\n74BInCipEEmQMJFYsMD/uNxzj/+xkdyw7bbw1FPQsCG89lrJvCQicaGkQiRBrr/e3xysYUN48kn/\nIyO5pUcPuPdev3zddT65EIkLJRUiCfH++34aboC774Zddok0HInQ0KFwzjn+Tqannw6rVkUdkYin\npEIkAdat8zcKKy72PyjDhkUdkURt1CjYcUf47ju45JKooxHxlFSIJMAf/gAzZ0LbtnDHHVFHI3HQ\nsGHJZGePPOKHFYtETUmFSMy9/Tbcdptfvv9+2G67aOOR+BgwAC6/3C+fe64fGSQSJSUVIjG2dq1v\nM3cOzjwTDjss6ogkbq65Bnbd1c9bMXy4/66IREVJhUiM/e538PXX0KGDn0FTpLR69XzzR0GBvwHZ\n449HHZHkMiUVIjH1+utw111++cEHoXHjaOOR+OrTB666yi9fcAEsWhRtPJK7lFSIxNCGDX7IIPjp\nmQ88MNp4JP5+/3vo29ffufayy6KORnKVkgqRGLr5Zt/s0batXxapTJ06flIsMxg7FiZMiDoiyUVK\nKkRi5uuv4YYb/PKtt0KjRtHGI8mx++5+FAj4ZpBNm6KNR3KPkgqRmLn0Ut/8MXgwnHBC1NFI0vzl\nL9C8OXz+Odx+e9TRSK5RUiESI//9r38UFMDo0bpZmFRf06YlTWbXXAPz50cajuQYJRUiMbF+fcl0\nyyNH+ruQitTEsGEwcKCf52TkyKijkVyipEIkJm66CebMgXbtSoYHitREXp6/6Vxenr+b7SuvRB2R\n5AolFSIx8PXXcOONfnnUKN3SXLZer15w0UV++aKL1GlTaoeSCpEY+P3vYeNGPx/FccdFHY1ki2uv\nhRYt/M3oHnww6mgkFyipEInY5MnwxBO+U+Ytt6hzpqRP48Zw9dV++ZprfB8LkUxSUiESsSuu8M9D\nh8Juu0Ubi2Sfc86B7bf3U3ePGhV1NJLtlFSIROjVV+G116BuXbjuuqijkWxUt66fuwL8UFPdHl0y\nSUmFSESKi0tqKYYPh86dIw1HstgJJ/j7gnz/fUmCIZIJSipEIvLkkzBlip+G+//+L+poJJvl5fkh\ny+CHms6ZE208kr2UVIhEYNMm+MMf/PJll/ke+iKZdOCBcNBB/runeVAkU5RUiETg/vth9mxo2VIz\nHkrtCedCGTsWPv442lgkOympEKll69aVdMq86ipNdCW1p29fOPlkvxwONRVJJyUVIrXs/vv98L7O\nnf1wP5HadM01vo/Fc8/BJ59EHY1kGyUVIrXoxx/hr3/1y1dc4Yf7idSmHXf0o0EArr8+2lgk+8Qi\nqTCzfczsP2b2nZkVm9nRlWy/X7Bd6qPIzFrWVswiNfHoo/5W1G3a+DtJikThyiv985NPwowZ0cYi\n2SUWSQWwDTAVOB9wVdzHATsCrYNHG+fcksyEJ7L1iopKOsr95jdQv3608Uju6tkTjj4anCv5Toqk\nQyySCufcS865q51z44Hq3PlgqXNuSfjIVHwi6fDUU37ER9OmcO65UUcjuS6srRgzBr75JtpYJHvE\nIqmoIQOmmtkCM3vFzAZGHZBIeZwrab++5BKN+JDo9e8PBxwAmzf76btF0iGpScVC4FzgWOBXwLfA\nW2bWO9KoRMrxwgvw6ac+mbjwwqijEfHCmVwfeAAWLow2FskOBVEHUBPOuZnAzJRV75vZDsAIoMLu\nbyNGjKBx48ZbrBsyZAhDhgxJe5wi4GspwvstDB/umz9E4mD//WHAAHjvPbj11pKRSZLbxo0bx7hx\n47ZYt3r16irtm8ikohwfAHtVttGoUaPo27dvLYQj4k2c6E/a9erBiBFRRyNSwszXVhx5JNxzjx/m\n3KxZ1FFJ1Mq60J4yZQr9+vWrdN+kNn+UpTe+WUQkVm64wT//+td+KKlInBx+OPTuDWvXwj/+EXU0\nknSxSCrMbBsz65XSJ2L74HWH4P0bzOzhlO0vMbOjzWwHM9vFzG4DBgGjIwhfpFzTp8PLL/srwssu\nizoakZ8zK6lBu/tu33FTpKZikVQAuwMfA5Px80/cAkwBrg3ebw10SNm+brDNp8BbQE/gAOfcW7UT\nrkjV3HWXfz7qKNh++2hjESnPiSf6O+XOnw/PPht1NJJksUgqnHMTnHN5zrn8Uo9fB++f4ZwbnLL9\nX51zOzrntnHOtXDOHeCcmxhdCUR+bs0aeDioX7voomhjEalIvXol96G5885oY5Fki0VSIZKNHn4Y\nfvgBdt7ZzwcgEmfDh0N+vu9YrBuNSU0pqRDJgOJiGB308LnwQt9uLRJn7drBscf6ZdVWSE0pqRDJ\ngNdeg5kzoVEjOO20qKMRqZqwmW7sWFi+PNpYJJmUVIhkQHild/rpPrEQSYK99oI+fWDDBj/Lpkh1\nKakQSbOvv4bnn/fLF1wQbSwi1WFWUltx110aXirVp6RCJM3uvttPzX3wwdCtW9TRiFTPSSf5WTXn\nzYP//jfqaCRplFSIpNG6dSXVxhpGKknUoEHJ8NI77og2FkkeJRUiafTYY7BqFXTpAocdFnU0IjUz\nfDjk5cFbb8GMGVFHI0mipEIkjcJainDMv0gSdegAhx7ql//5z0hDkYRRUiGSJtOnw/vv+2Ti1FOj\njkZk65xxhn9+5BEoKoo2FkkOJRUiaRJOyX3YYdC6dbSxiGyto46Cpk1hwQJ45ZWoo5GkUFIhkgZF\nRfDoo345vMITSbJ69eCUU/zyQw9FG4skh5IKkTR45RV/RdesGRx5ZNTRiKRHmCCPHw8rVkQbiySD\nkgqRNAg7s518MtStG2koImnTpw/06gU//uhHNolURkmFyFZasQKefdYvq+lDsk34nVYTiFSFkgqR\nrfT44/5KbrfdoHfvqKMRSa9TToE6dWDKFPj006ijkbhTUiGylcKmjzPO0C3OJfs0b+5HgoBqK6Ry\nSipEtsLnn8OHH0JBge9PIZKNwiaQMWN8rZxIeZRUiGyFsJbiiCOgZctIQxHJmEMP9XOvLFtWcgde\nkbIoqRCpoc2bNTeF5IaCgpJZYjVtt1RESYVIDb32GixeDC1awOGHRx2NSGaFScVLL8Hq1dHGIvGl\npEKkhp54wj8ff7zvHS+SzXbdFXbe2fep+M9/oo5G4kpJhUgNbNpUMjfF8cdHG4tIbTCDE07wy2FC\nLVKakgqRGnj9dVi5Elq1gn32iToakdoRJhUvvwyrVkUbi8STkgqRGgiv1I491t/qXCQX7LKLf2za\npCYQKZuSCpFqUtOH5DI1gUhFlFSIVJOaPiSXhYn0K6/4vwORVEoqRKpJTR+Sy7p3h549fY3d+PFR\nRyNxo6RCpBrU9CGiJhApn5IKkWpQ04dISUL96quwYkW0sUi8KKkQqQY1fYhAt27Qq5efqj6suRMB\nJRUiVaamD5ESagKRsiipEKkiNX2IlAgT69deg+XLo41F4kNJhUgVqelDpMSOO0KfPlBUpCYQKaGk\nQqQKNm8uGT6npg8R75e/9M/PPRdtHBIfSipEquD9930v96ZNYe+9o45GJB6OPNI/v/oqbNgQbSwS\nD0oqRKogvBI79FAoKIg2FpG46N0b2raFtWthwoSoo5E4UFIhUgXPP++fjzgi2jhE4sSspLZCTSAC\nSipEKvXNN/DZZ5CX52sqRKREalLhXLSxSPSUVIhUIqylGDjQ96kQkRKDB0O9ejB3LkyfHnU0EjUl\nFSKVUNOHSPm22cYnFqAmEFFSIVKhdevgjTf8cljNKyJbUr8KCSmpEKnAG2/4oXIdO8Iuu0QdjUg8\nhbV4kybpBmO5TkmFSAVSmz7Moo1FJK46dYKePaG4GF56KepoJEpKKkTK4VxJUqGmD5GKqQlEQEmF\nSLmmTYNvv4UGDWDQoKijEYm3sAnkxRf9tPaSm5RUiJQjrKUYPNgnFiJSvj339EOuV62Cd9+NOhqJ\nSrWTCjM71Mz2Tnl9gZlNNbPHzGy79IYnEh01fYhUXX4+HH64X1YTSO6qSU3FX4FCADPrCdwCvAB0\nAW5NX2gi0Vm+HN57zy+HJ0oRqZj6VUhNkoouwBfB8rHAc865K4ELgMPSFZhIlF56yfdk79nTDycV\nkcodcoivsZg+3c+wKbmnJknFj0DDYPlA4JVgeQVBDYZI0r34on/WLJoiVdekCfTv75dffz3aWCQa\nNUkq3gFuNbOrgF8AQcszOwHz0xWYSFScKzkhHnJItLGIJE04ZbeSitxUk6TiQmAzcBww3Dn3XbD+\nMEDTnkjiffklLFoE9ev7Hu0iUnUHHOCf33hDdy3NRQXV3cE5Nw/4WX9459yItEQkErHwXh977eUT\nCxGpugED/BDsxYvhiy80vX2uqcmQ0r7BqI/w9TFm9qyZXW9mddMbnkjtC6ttw2pcEam6evVg72DS\nATWB5J6aNH/ci+8/gZltDzwOrAOOB25OX2gita+oCN56yy+H1bgiUj3h346SitxTk6RiJ2BqsHw8\nMNE5dzJwOn6IqUhiffIJrFwJjRpBv35RRyOSTGEt31tvacruXFOTpMJS9jsQP/EVwLdA83QEJRKV\nsD/FfvtBQbV7HIkIQN++fnjpmjUwZUrU0UhtqklS8RHwBzM7FdiPkiGlXYDF6QpMJAphda2aPkRq\nLj8f9t/fL6sJJLfUJKm4FOgLjAb+4pybHaw/DtBtZCSxfvwR3n7bL6uTpsjWUb+K3FSTIaWfAj3L\neOu3QNFWRyQSkQ8/hLVroXlz2HXXqKMRSbYwqZg0CTZs0PDsXJG2W5875zY45zal6/NEalvYn2LQ\nIMhL219phUVHAAAgAElEQVSGSG7aeWdo08YnFOHN+ST71WSeinwzu8zMPjCzRWa2IvWRiSBFaoP6\nU4ikj5mm7M5FNbke+yMwEvgX0Bh/u/OngWLgmrRFJlKL1q0ruZpSfwqR9FC/itxTk6TiFOBs59wt\n+HuAjHPOnQVcB+hOCZJI777rO2q2bw9du0YdjUh2CJOKDz/0w0sl+9UkqWgNTAuWf8DXVgA8B+hG\n0ZJIqU0fZtHGIpItOnb0SXpREUycGHU0UhtqklTMB9oEy18BBwfLewAbq/thZvb7oH/GGjNbbGbP\nmNlOVdjveDObbmbrzewTMzusuv+2SCjspKmmD5H0UhNIbqlJUvEMEHZluxP4k5nNAh4BHqzB5+0T\nfE5//AyddYBXzKxBeTuY2UDgMeA+oDcwHnjWzHrU4N+XHLd6NXz0kV8eNCjaWESyTZioh4m7ZLea\nzFNxRcryv8xsHjAAmOWc+28NPu/w1NdmdjqwBOgHvFPObhcDLzrnbg1eX21mBwEXAudXNwbJbW+/\nDcXFsOOO0KFD1NGIZJd99vHP06b5fhWFhdHGI5m11aPxnXPvOedurUlCUY4mgAMqGp46AHit1LqX\ng/Ui1TJpkn/eb79o4xDJRm3aQJcu4By8/37U0Uim1eiWSWbWFtgbaEmpxMQ5d0dNgzEzA24D3nHO\nfVHBpq35+X1GFgfrRarl3WBy+YEDo41DJFvttRfMmeP/1g4+uPLtJbmqnVQEzRP3Aj8Cy/G1CiEH\n1DipAO4GegB7bcVnVGjEiBE0btx4i3VDhgxhyJAhmfonJcY2bYIPPvDLSipEMmPgQBgzpqRWUOJt\n3LhxjBs3bot1q1evrtK+Namp+BN+ToobnHPFNdi/TGY2Gjgc2Mc5t7CSzRcBrUqtaxWsr9CoUaPo\n27dvzYKUrDN1qp9GuGlT2KnSMUciUhN7BZeJ77/vh5fm50cbj1SsrAvtKVOm0K9fv0r3rUmfiobA\n4xlIKI4BBjnn5lVhl/coGYESOihYL1JlqU0fmp9CJDN22cV30PzhB99hU7JXTZKKB4Dj0xWAmd2N\nn6XzZGCtmbUKHvVTtnnYzK5P2e124FAzG2lm3czsGvxokdHpiktyg/pTiGRefj7sGcy3HP7NSXaq\nSVLxe2A/M3vLzO40s1tTHzX4vPOAQuAtYEHK44SUbTqQ0gnTOfcePgk5B5gK/Ao4ppLOnSI/o6RC\npHaEf2PqV5HdatKn4vfAIcCM4HXpjprV4pyrNLFxzv1snkPn3L+Bf1f33xMJffstzJ/vr6L22CPq\naESyW5hUqKYiu9UkqfgN8Gvn3D/THItIrQpPbr17Q8OG0cYiku3694e8PJg7FxYsgLZto45IMqEm\nzR8bAVVgSeKp6UOk9hQWQs+eflm1FdmrJknF7cBF6Q5EpLYpqRCpXeHQUiUV2asmzR+/AAab2ZHA\n58Cm1Dedc79KR2AimbR2LXz8sV9WUiFSOwYOhLvvVmfNbFaTpGIV8HS6AxGpTR995CfhaddONxET\nqS1hTcWUKbB+PTQo917UklTVSiqCe3P8EVjqnFufmZBEMk+TXonUvk6d/A3GFi6EDz+EffeNOiJJ\nt+r2qTBgNtA+A7GI1Br1pxCpfWYaWprtqpVUBFNzzwKaZSYckcxzTkmFSFTCJhD1q8hONRn9cQXw\nVzPbNd3BiNSGmTNhxQqoX9/PUSEitSe1psJVe7pEibuadNR8BH9TsU/M7Edgi74Vzrmm6QhMJFPC\nWoo99oC6daONRSTX9OnjE/oVK2DGDNh556gjknSqSVJxadqjEKlFavoQiU7duj6hf/tt3wSipCK7\nVDupcM49nIlARGqLkgqRaA0Y4JOKDz6AM8+MOhpJpyolFWZW6JxbEy5XtG24nUgcrV4NXwT3sh0w\nINpYRHLV7rv758mTo41D0q+qNRUrzayNc24JfvKrsrrXWLA+P13BiaRbOItmp07QokW0sYjkqn79\n/PO0afDjj+rblE2qmlQMBlYEy4Mq2K7n1oUjklnhlVHfvtHGIZLLunSBJk1g1Sr47DP9PWaTKg0p\ndc5NcM5tTln+6QFMAboBf8XfbEwktqZM8c/hlZKI1D6zkkRCTSDZpSbzVABgZvua2cPAQuAy4A1g\nz3QFJpIJ4QlMSYVItMK/QSUV2aW69/5oDZwOnAkUAk8A9YD/55z7Iu3RiaTR99/7ia9A1a0iUVNS\nkZ2qXFNhZv8FZgC74eeqaOucuyhTgYmk29Spfga/9u2hZcuooxHJbWFS8emnvrOmZIfqNH8cBjwA\n/NE597xzrihDMYlkhDppisTHDjtA48Y+ofj886ijkXSpTlKxN9AImGxm/zOzC82seYbiEkk7ddIU\niQ911sxOVU4qnHPvO+fOBtoA9wInAQuCzzjIzBplJkSR9FBNhUi8hAl+mPBL8lV79Idzbq1z7kHn\n3N74eSluwd+5dImZ/SfdAYqkw9q18OWXflk1FSLxoM6a2afGQ0oBnHMznHOXA+2BIekJSST9PvkE\niouhdWto0ybqaEQESpKKTz6BTZuijUXSY6uSipBzrsg596xz7uh0fJ5Iuqk/hUj87LADFBbCxo0l\n9+SRZEtLUiESd+pPIRI/eXnqrJltlFRITlBNhUg8qV9FdlFSIVlv/fqScfBKKkTiRUlFdlFSIVnv\n00+hqMjf6rxdu6ijEZFUqZ01N2+ONhbZekoqJOulNn2YRRuLiGypa1do1Ag2bFBnzWygpEKynjpp\nisSXOmtmFyUVkvXUSVMk3tSvInsoqZCstnEjfPaZX1ZNhUg8KanIHkoqJKt99pmfqa9pU+jUKepo\nRKQs6qyZPZRUSFZL7U+hTpoi8bTjjr6z5vr1MH161NHI1lBSIVlN/SlE4i8vD3r29Mthc6Ukk5IK\nyWpTp/rnPn2ijUNEKqakIjsoqZCs5VzJTJrhCUtE4mnXXf3ztGnRxiFbR0mFZK158+CHH6BOHd9m\nKyLxFSYVqqlINiUVkrXCk9POO/vEQkTia5dd/POcOf5iQJJJSYVkrTCpCE9WIhJfLVpAq1Z+WdN1\nJ5eSCslaYVIRVquKSLyps2byKamQrBV20lRSIZIM6leRfEoqJCsVFZVUoSqpEEkGJRXJp6RCstJX\nX/n7fjRoAF26RB2NiFSFhpUmn5IKyUrhlU6PHn62PhGJvx49/POiRbBsWbSxSM3odCtZSf0pRJKn\nUSPo3Nkvh3/DkixKKiQraeSHSDJpBEiyKamQrKSkQiSZ1Fkz2ZRUSNbZuBFmzvTLSipEkkVJRbIp\nqZCsM3MmbN4MhYXQrl3U0YhIdaQmFc5FG4tUn5IKyTqpnTTNoo1FRKqnWzfIz4dVq+C776KORqpL\nSYVkHfWnEEmuevV8YgFqAkkiJRWSdZRUiCSb+lUkl5IKyTpKKkSSTUlFcimpkKyybh18/bVf1i3P\nRZJJSUVyKamQrDJ9uu8x3qIFtGwZdTQiUhNhUvHFF/7mgJIcSiokq6jpQyT5tt/e3wxw/XqYMyfq\naKQ6lFRIVlFSIZJ8+fklNxdTE0iyKKmQrKIbiYlkB90GPZmUVEhWCa9q1ElTJNnUWTOZlFRI1li9\nGr791i8rqRBJNtVUJJOSCskaYdNH+/bQpEm0sYjI1gn7VMye7e/lI8mgpEKyRphUqJZCJPnat4f6\n9WHTJvjmm6ijkapSUiFZY8YM/9y9e7RxiMjWy8uDHXf0yzNnRhuLVJ2SCskas2b55/BEJCLJttNO\n/llJRXIoqZCsEZ54whORiCSbkorkiTypMLPfm9kHZrbGzBab2TNmVuHPgpkNM7NiMysKnovNbF1t\nxSzxs3kzfPWVX1ZSIZIdlFQkT+RJBbAPcCfQHzgQqAO8YmYNKtlvNdA65dEpk0FKvM2b5zt01a/v\nO3iJSPIpqUiegqgDcM4dnvrazE4HlgD9gHcq3tUtzWBokiDhSadrV9/BS0SSL0wq5s3z9wFpUNml\npkQujqffJoADVlSy3bZmNtfM5pnZs2bWoxZik5hSfwqR7NOsGWy3nV+ePTvaWKRqYpVUmJkBtwHv\nOOe+qGDTGcCvgaOBU/DleNfM2mY+SokjjfwQyT5magJJmsibP0q5G+gB7FXRRs6594H3w9dm9h4w\nHTgX+GNF+44YMYLGjRtvsW7IkCEMGTKkhiFLHKimQiQ77bQT/O9/Sipq07hx4xg3btwW61avXl2l\nfWOTVJjZaOBwYB/n3MLq7Ouc22xmHwNdK9t21KhR9O3bt4ZRSlwpqRDJTqqpqH1lXWhPmTKFfv36\nVbpvLJo/goTiGGCQc25eDfbPA3oC1UpGJDts3Fgyja+aP0Syi5KKZIm8psLM7gaG4PtHrDWzVsFb\nq51zG4JtHga+c85dGby+Ct/8MRvfsfNyoCNwfy2HLzHw1VfgHBQWQsuWUUcjIumkpCJZIk8qgPPw\noz3eKrX+DOCRYLkDUJTy3nbAP/DzU6wEJgMDnHNfZjRSiaXUpg+zaGMRkfTqGjRqL1sGK1ZA06bR\nxiMVizypcM5V2gTjnBtc6vVIYGTGgpJECZMKNX2IZJ9tt4V27eC77/wor/79o45IKhKLPhUiWyMc\nTqpOmiLZSU0gyaGkQhJPIz9EspuSiuRQUiGJp+YPkewWJhUzZkQbh1ROSYUk2vffw6JFfllJhUh2\nUk1FciipkEQL+1O0bAlNmkQbi4hkRphUzJoFxcXRxiIVU1IhiaamD5Hs16UL5OfDunWwYEHU0UhF\nlFRIomnkh0j2q1MHtt/eL6sJJN6UVEiiqaZCJDeoX0UyKKmQRNNwUpHcoKQiGZRUSGI5p6RCJFco\nqUgGJRWSWMuXw6pVfnmHHaKNRUQyS0lFMiipkMQKTy4dOkDDhtHGIiKZFSYVX38NmzZFG4uUT0mF\nJJZGfojkjrZt/cVDURHMmRN1NFIeJRWSWBr5IZI78vJK/tbVBBJfSioksdRJUyS3qF9F/CmpkMRS\nUiGSW5RUxJ+SCkmk4mKYPdsvq/lDJDeEs2qqT0V8KamQRFq40N8HID/f3xdARLJf587+ee7cKKOQ\niiipkEQKr1Q6dvT3BRCR7BdeQMydq7uVxpWSCkmk8EolvHIRkezXvr0fBfLjj7BoUdTRSFmUVEgi\nKakQyT116vjJ7kBNIHGlpEISKWz+UFIhklvCv3l11ownJRWSSOFVijppiuSW1H4VEj9KKiSR1Pwh\nkptUUxFvSiokcYqKYN48v6ykQiS3qKYi3pRUSOIsWACbN0NBgb/JkIjkDtVUxJuSCkmc8AqlY0c/\n+ZWI5I6wpmLePF9rKfGipEISRyM/RHJX27Z+aOnmzfDdd1FHI6UpqZDE0cgPkdyVn+9rKUH9KuJI\nSYUkjkZ+iOQ29auILyUVkjhKKkRym0aAxJeSCkkcJRUiuU01FfGlpEISZfNm+PZbv6ykQiQ36Rbo\n8aWkQhIlnKOiTh1o0ybqaEQkCmr+iC8lFZIoYXVnp06ao0IkV4U1Fd9+C5s2RRqKlKKkQhJF/SlE\npHVrqFcPioth/vyoo5FUSiokUZRUiEhenq+tBHXWjBslFZIoSipEBNSvIq6UVEiiKKkQEdCw0rhS\nUiGJoqRCREA1FXGlpEISI3WOCt33QyS3qaYinpRUSGLMn+9vdVy3ru/9LSK5SzUV8aSkQhIjPHl0\n6uR7f4tI7gprKhYsgI0bIw1FUujULImh/hQiEmrRAho2BOdg3ryoo5GQkgpJDCUVIhIyU7+KOFJS\nIYmhpEJEUqlfRfwoqZDECK9GNPJDREA1FXGkpEISQzUVIpJKNRXxo6RCEmHTppIbBympEBEoORco\nqYgPJRWSCPPn+zsS1qsHrVpFHY2IxEFYU6Hmj/hQUiGJoDkqRKS0sKZi8WJYvz7SUCSg07MkgvpT\niEhp220HhYV+WU0g8aCkQhIhPGFo5IeIhFLnqlBSEQ9KKiQRwjZT1VSISKoOHfxz2JFboqWkQhJB\nzR8iUpZ27fyzkop4UFIhiRDe8rxjx2jjEJF4ad/eP3/3XbRxiKekQmKvuLjkhBFelYiIQMk5QUlF\nPCipkNhbtsxPfmUGbdpEHY2IxElYU6Hmj3hQUiGxF16BtGwJdetGG4uIxItqKuJFSYXEnpo+RKQ8\nYU3FypWwbl20sYiSCkkAJRUiUp7CQthmG7+s2oroKamQ2AvbSpVUiEhpZupXESdKKiT2wquP8MQh\nIpJK/SriQ0mFxJ6aP0SkIqqpiA8lFRJ7SipEpCKqqYgPJRVZaNy4cVGHkFblJRXZVs7y5Eo5p70+\nLeoQakWuHM/aLGeUNRW5cjyrKvKkwszOM7NPzGx18HjXzA6tZJ/jzWy6ma0P9j2stuJNgmz6kq9d\nC6tW+eXSfSqyqZwVyZVyfvbGZ1GHUCty5XjWZjmjrKnIleNZVZEnFcC3wO+AvkA/4A1gvJl1L2tj\nMxsIPAbcB/QGxgPPmlmP2glXalN4kth2Wz90TESkNPWpiI/Ikwrn3PPOuZecc18552Y75/4A/ADs\nWc4uFwMvOududc7NcM5dDUwBLqytmKX2qD+FiFQmPD8sWuSn9JfoRJ5UpDKzPDM7CWgIvFfOZgOA\n10qtezlYL1lGSYWIVKZlSygoAOd8YiHRKYg6AAAz2xWfRNQHvgd+6Zz7spzNWwOLS61bHKyvSH2A\n6dOnb0WkybB69WqmTJkSdRhp8cEH/rl+fShdpGwqZ0WSWs6lS5eyfMlyfvjkB+o3ql/p9mtXrmX2\nh7NrIbKtt+H7DWxcspFPP/2UhQsXVmvfpB7P6qrtcjZv7hOKN96Anj1r7Z/NmeOZ8ttZ4R+zOecy\nH00lzKwA6Ag0Bo4Dzgb2LSuxMLONwGnOuX+lrBsOXO2cK/celmZ2MjA23bGLiIjkkFOcc4+V92Ys\naiqcc5uBr4OXH5vZL4BLgOFlbL4IaFVqXatgfUVeBk4B5gIbahysiIhI7qkPdMb/lpYrFklFGfKA\neuW89x5wAHBHyrqDKL8PBgDOueX4USMiIiJSfe9WtkHkSYWZXQ+8CMwDGuFrE/YDDg7efwSY75y7\nMtjlduAtMxsJPA8MwQ9FPbuWQxcREZEUkScVQEvgYaANsBr4FDjYOfdG8H57YHO4sXPuvaB/xF+C\nxyzgGOfcF7UatYiIiGwhFh01RUREJPliNU+FiIiIJJeSChEREUmLrEoqzMyijkFEcpvOQ9kll45n\nOsqaFUmFmfUys22cOoiISERy5TxkZgPMrLIZjBMvV44npLesiU4qzOxwM5sDPIOfNOtPZtYk6rjS\nzcyOMbOngvINSlmfVRl0jpXzRTO7x8xOMLOGwfqsKifkRllz6Dx0iJnNxc/3Mzk4pjtGHFba5crx\nhMyUNbFJhZntBtwM3IOf2vt2/CycfzCz5lHGli5mVsfMbgH+CSwH9sHfFv73ZlYnWzLoHCpnnpld\niy/nF/iZYEcDt5tZPeecy5Yf21wpay6chwDMrBNwHfAgcCAwAjgEuN7MOkQZWzrlyvGEDJbVOZeo\nByXDYIcCC4HmKe9dCrwP/C7qONNU1k74E/JhKeuuBz4ALg1e50Udp8pZ5XK2BKYAp6asOwf4GLgp\nW8qZC2XNlfNQSjkPB9YCnVLeOwF4G7gr6jh1PONT1sTUVIRXNS4oOdAF+Ar4MWWze/A/TgebWS3e\npy5jWgNNgJUp627D/9iebmadnXPFWXDFl9XlTIm7Pn4ytxUpbz8ePM42s52TXE7I/rLmynnIzOrD\nFuVsii9TnZTNngFeAfYzs71rN8L0yJXjCbVX1tgnFWZ2oJm9ANxjZmeYWePgrXeB/kCHYLt859xG\n/EmrEZCoL7mZ/cLMDjazzikn2npAEdAo5QuxBP/H/D3B1OQpX5LYy6FyDjSzU8xsF6BBsLoJsABo\nb2Z5AM65NcB/gKnA/wXrElNOyI2y5tB56CAzewt41Mx+Z2ZdgrdmAD2BncJtnXObgBfwx/mY2o51\na+TK8YTaL2uskwozGwE8jc+mGgG3AHeYWXPn3OvA5wQnJ8ABOOdeATYBvYLPiPWVkJltZ2b/AV7D\nt2+9D1wB4JybiK9yPBnIT9ltEjAT2MXMmtVuxDWTQ+UsNLN/Ay/h2yffAu4EcM59CiwDjsQ3D4Rm\nA68DXVNO4rGXK2XNhfMQgJmdCTyFryGcD5wL/NPM2jjnPsT/CI00s23CfZxzk/HHuYsFIgi9WnLl\neEJEZY26faeCdp9C/Jf4wpR1Q/Ff+L8Hr08AioGDSu37IDAx6jJUsZznB2XaCd8McB3+xPuX4P0T\n8dVT+xC0hQXrz8PfhM2iiFvlLLecJ+HvX9Mt+CP+NTAHuD94f0DwnT0VqJey37HAd0CzqMugsm5R\nxlw5D9XH36DxLynr9g3K+ULweg98jeJZQN2U7W4Gvoi6DDqe8Shr5AWv4D+kOT5bTu3o1QC4AN/2\nvmew7jHga2AQ/ip3O/xV8HlRl6GK5fxfeICD19vge1ZvBnYK1r0RlKlHyna/w9/uvVHUZVA5tyjn\n88BTKa8t+BEtBg4P1j0GTAYOTdnuNPxVQ8uoy6CyblHGXDkPGf5qdmTKugLgiOB4HhmsGwV8C5wR\n/D80Bl4F/i/qMuh4xqOskRe8gv+QTsBHwG+AgpT12wPPAf8JXlvwpV4QPH8T7Ncl6jJUUj4LHo8D\nj5R6ryXwDvBy8LotvhngbeAyYBj+Si/2vZFzqJx5wfNdpGT4lPS0fgqYHCw3w3dw+xK4CRge/PHf\nQAJqZHKsrFl9HkopT9OgPHcD9VPWbwuMBaalrPs7sAif7M/D11btFHUZdDzjUdaoC92A4Aq0rBNM\ncHJ6C+iYsi4P3347GegdrGuOz7KuBM6I+mBWUub8Uq/vBF4Edk1ZF17xrQb6B+t2x18lTMLf7v3S\nqMuicpZZzpFB7IOC1+EP8O74dsoDg9dd8H1KXgCmxb2cuVbWUuXMuvNQOeUchU/o90xZZ8DR+Gat\n8Dg3APrim0FOiTruXD6exPA3NMr/jGvwVy/dw4IGzycDewXLPfBteOcDDVL2HQwsDf+jyvrPjPsD\n6BA8DwLmBgc5NZvcGd/2dVqp/ZqWPrnH+ZFt5cRP/FO3jPU7Bs+746sOR7FlP4KmwJvAH0vtV5+E\nzdeQTWUF2pZ6Hda2ZN15qJwfnVOBY4PlHfDV4DeQ0t8Fn0DMBwbEvZzAKWzZfJrNx/MaYvgbWuuj\nP8xsqJktww8T3An/o4LzY9a7AmOAQ82soXPuC3x13IX42dtC2+Kvbgn2dbUVf1WZ2fFm9oCZXWxm\n3VPWNzez54AbzM8W+SYwAT+j2QEpH7Ec2BVYF+wXDstb4ZwrqrWCVMLMfmVmD5nZdWY2MGV9sywr\n54lmNg8/KdfOYY9oM2tkZo8AY8zPFPkRvgqxP/4EF9oE7AgsCfYLy7nBOVdci0WpVPDdvc/MLkkd\nq25mjbOlrGY2xMw+Av5tZs+b2ZHgzyVZdh461sxmmFn3oGyW8l57fIfpo8ys0Dn3FXA/voynpXxM\nOOR7DcSznABm9iDwKHBSOEIl244nJOA3tBazqt2A6fjhR2cD7fBVp78plWWdypZV5HXw7T9f4KvQ\nf43vUHQPMbuSxVcVbov/Yi/Dt0V+js8If5my3V3AXwlmMgM64jspvoPvcb1t8P8wmZi2VeJ/NCbi\n+zzci79iXQnskrLNnUkvZxD3Ofi240vwfQTql3r/MvzVepvgdSvgIXzb5FFAG3wzz2ekVC3H7RGU\n7Un8LHv34KvC5wPDUrb5bZLLiq9FeSD43l6E73D4JLCKLWvQEnseCuLdFrg2OC4bgedLvR+eb48s\n9TdbANyBb18fh+9MPSs4lzWMulzllDU/eL4JP9phOTAwWBfWVCT6eAYxJ+I3tDb/Q64OTj7bpqz7\nFLg9WC6rSjn8QrQLTgAvAp8Q457GQD/8UMndUtY9jW/X+lXwulHKe+EXYWDwR7wmKOMPpAwFitMD\nX006HriPLatJ1wBnp7xumLKcuHIG8TbAdzQ8LywTvkqxdeo2Kcvhd3Z7fPK4PDgRfA+MiLo8lZT1\nOPwonXYp654Kvs/HlXFME1dW/NXa+8DuKeva4zsIn1bOPkk8D+2IT/bPwU+xvRk4Kngvjwqqu/EJ\nyXH4ROI94Mqoy1OF8m6Lb0Ztie9k+BTljC5K4vEM4k3Eb2im/xMGEFzRlFpfEDyPAj6oYP9OpV7X\nA+pEfXDLiTU8eKcFJ+HUjjE7BQfz30DT1O3L+JzewR90bIdQ4ocdnQy0T1nXNChjtyp+RuzLGcTZ\nBN/OvDt+TPf84EdpBX6e/CbBdmX2FQA6Bz9k29ZWzFtR1qeBfwfL2wbPp+GHFL4OtAjWlffdjX1Z\n8QnQaaRcoQU/RN+SMuy11D6dSr2O5Xmo9PkWf4GTj69BfQCYXcn+3ctYV5DOGDNUzrygnG8Ef68H\nBN/ZsGNps1L7J+J4llXW0scljr+hmfqPOCA4Ec8NTsL/oGQugryU7W7EXxm1KOMzrg2+JL2iPrAV\nlHPf4CSaWm16ErCYkg6KYbJxTlDWYSnb7oDvbbtLbcSb5nKmLh8ZHOdV+Mz5MaBflpSzN35a6fPw\nzT2n4tswb8JfGd2Qsm0H/ERAYYe2WFajVlDWm4AvS213A34W1ElsWQsV+7ICB+Gr8i8lGF1Uxjb5\n+KF33wJ9y3g/Ceeh1PPtt/h+EV1T3jdgF3wi/Nuyjhd+7oKPgIOjLk81ypn6u9IS39QTnnefxDfb\nzAGuT9LxLKOsifoNzcR/Rgd8ldl1wY/Jcfj2m38TZE0E1TT4arkfgcYp+4c/wkcE/1ld0x1jGsrY\nHHgYnw1PBTqnvFcP34fi8uB1mFE2wbf3/S1l3TH4KuO4VhWXW86Ubdri+5D8CX9vgMH4xGIswRVC\n0suJ7xezHngsZV0B8OfgmHYO1u2Nr6W6Neoy1fC7uz2+c+UEfN+Jd4MT2+Bg2+tStt0rrmXF9+v4\nLz65HxN8H1cBe6Rskzpr6y+Bz0t9RuzPQ0F8lZ5vg+3y8In99wS1paXKuQ++3T2WP7aVlTNY92zK\nttu4l4IAAAzYSURBVO/iO5e+QsoPbtyPZ1WPKTH+Dc3Ef8hB+J78O6Ss+1Vwovp7qW374zOxX6Y7\njgwe8AL8BD4v4avD1+LH4KdOW3s1vtNiWF0c9ie4Ffio1OcdHnWZalLOlC+uEXRcpKTD1Ijgj6Aw\n6eUMtglniXyo1L6/wv/opjYD7RV1mWpQ1tThoHvh+8pMxnfqCjvZPkLK7JlxLSu+z8s/8ZOtdUlZ\n/z/gwWA5r9Q+DwEPp7xuTcyb5VJirc75ti0+QX4oeN0Hn0zEaphvNcv5Nr4JoDu+ZuIh/I/sA/g+\nJd8APaOOP4PHNHa/oZkYUtoU31kr9cZQ44GXgX3MbL+U9cvwPVPzIRk3aXHObQamAHc5557AVxmP\nxHfeC92Pn3HuTjPb1pUMoysE5phZvpnlB5/3Qu1FX3WVldMF32jnbQh3C45hW3zv8aKUz0tkOYNt\n/o2vRtzHtrzFcx38ld/6lG0n1UbcNVFBWbunbDPJOXc2vlnjIufcMjNrif8B+hjAzArCbWu7DJVx\nzq3Dj3b4p3NuThgrfuKt8LtbDP6ujOZv8d0f+K+Z5ZnZX/Df3SNrP/oaqfL51jm3AD+3wVAzexqf\nOA6ERJx7yyvni/iatBPxQyS7AUc45850zp2Lv+q/yMzq1XK8WyPZv6EZyLJ2xZ9kjy61vjf+CunW\nUutn4k9ykWdY1SijlXodDqtMvTIfiB/Z8AS+n8Uw/FC9U2sz1gyV82dXcZTUUoT3dSizJ30cH1Up\nJ35WyJn4E/Fl+Ds4zgWuijr+dJaVLZsF6uNPWMPxyUgirvhI6YhGSS3hWOAfZZRxB/xtva8Ljud0\ngg5+SXhU53yLH8n0B3yt23vAflHHn4Zy9gWexfc52J1gdBIlTcyHUkZfmTg/kv4bmqn/lBfwQyi3\nLbX+keALEP6hNwJuI7hZTdIelDQFHI+f8Kf0nd6Oxg9tmhqcsM6JOuYMlbMv/pa6E/GZc2ISihqU\n8w58rcVnST2eVSxruyCZ+BDfH2ZI1DFvZXnfIegkXSqpODX4kV1IjIc2V1K2Ss+3wbn25uDH6uwo\n4sxgOccAzxDcZyjqOGvrmAavY/cbmqn/kF7Bieo8tmyb/jMwK+pCZ6jM7+JnF2wZvE7tHNQp6vgy\nWU78HUfvBa6OOr4MlrNlyntNoo4v08c0eB5CMLFOkh/4DqiL2HJEUphU7UrC7kVSRvmqdL4lYVfs\nNS1nNjySXNawvTGtnHOfmNlNwFXAJjN7HJ8t747PKrOGmRU431Z9Nn5SkZPMbAdggJmd75z7yDn3\nTbRRbr0Kyrk3MBS4IHg/0So5nsOdc5Odc6uijTI9KjqmZjbMOTcu2gi3jpmZ82fivYEfnHOTg/V/\nBNqY2R+cc5/ha50Sq6rnW+fclIhCTItc+l1JclnDHvyZ+XCzu/DDtebhe1SvBY53fj7yrGNmH+AP\n+jzgXOfcyxGHlBEqZ/bJ5rKa2Wj8uec1fNt7Q3zfplciDSzNcuV8myvlhGSWNdNJRX18z/K+wEbn\nXKwzrJoKru6exVezXuyceyDikDJC5cw+2V7W4Bw0Dd8h80f8nVNvijaqzMih821OlBOSWdaMJhW5\nwsw6A6cDNznn1le4cYKpnNknF8pqZq/i5zAY6UqGP4tIBiipEJGsZmb5zrmiyrcUka2lpEJERETS\nIhMzaoqIiEgOUlIhIiIiaaGkQkRERNJCSYWIiIikhZIKERERSQslFSIiIpIWSipEREQkLZRUiIiI\nSFooqRCpgJntZ2ZFZlYYwb/dysxeNbMfzGxFbf/7knnB96s4+I49HYM4iqOMQ5JPSYXkrJSTeXEZ\njyIzuxqYBLRxzq2JIMQRQCtgN2CnCP79LZhZgZndZGafBonOd2b2sJm1KbXddmY21sxWm9lKM7vf\nzLZJeb+emT0UfM6m8n7EzKyumf3FzOaa2QYz+9rMTq8kxn3M7D9BbMVmdnRNylDOZ99uZh8FsVR4\nG3Ez62pm31cxGXT44/tT2czsTTO7tYzPPd3MVgbLc8r57obf3wdT9htkZs+b2TIzW2tmn5nZ38ys\nbbDJJPxdMJ+oQrwi5VJSIbmsNdAmeL4UWI3/EQ/X/805t9k5tySi+HYAJjvnvnbOLStrAzMrqMV4\nGgK9gWuBPvhbMncDxpfa7jH8nRUPAI4A9gXuTXk/H1gH3A68WsG/9yQwCDgD/6M7BJhRSYzbAFOB\n8/E/1jUtQ1kc8ADweEUbBcfkMWBCFT4ztLSKiaujpFy747+rrYFjg/U7UvL9vSSI51z8//MC4Ff4\nY3MeUAiMBEj5nmflTeWkFjnn9NAj5x/AMGBFGev3A4qBwpTtVuJ/LL8E1uKv7hoE780BVuB/MC3l\nc+oCfwPmAz8A7wH7VRDPHKAo+LeLgAeD9cX4H4TxwedcnRLn/4AN+B+PG4C8lM97E7gDGBXEtwg4\nE/8j+yCwBn8nz0Or+f+2exBf++B19yDGPinbHAJsBlqXsf9DwNNlrD80iLPJVhzTYuDo6pahCtv/\nEZhSwfs3AQ+X950q4/tVFH6/Sh2vW6v5PS3rc9oF34m/lfPvl96+zOOhhx5VfaimQqRypa94GwIX\nASfgfzAHAc/gfwgPA4YC5wLHpexzF9A/2Kcn/ir8RTPboZx/c3fgZeBf+CvPS1Le+yPwNLAr8GBQ\nhf08PqnYDZ90nAn8odRnngYsBfbAJxh/D+KYhL9qfwV4xP5/O3cXYlUVhnH8/2hiljFiNQp+QBZU\nEEhFVGCFhheFgRdS5EVFQWRCEanRRSSCIVkRGHlhYgYaoV0VViQJgfShoSahY1Cilkqp4fiBSb5d\nvPvg9njOnHNmzjhmz+9mZq+999pr73Ox3r3Wu7Z0eU8Po8oI8vn8VWzfBRyJiC2lY9YXx9zZQr0P\nAZuBlyTtk9QlaXGLbWtW9T30mqQp5KjB7L7W1SYPA0OAxbV2xsBM69kl7EIOnZpdKi4DnomI3QCS\n1pKBRGdEnAR2StpABhtrJI0n58vHRcSBoo63JD1ADu1Xd/5ExCFJp4CTEfFH1e5VEbGysiFpIbAn\nIp4rinZJehVYBCwonbctIl4rzlkEvEwOuy8vyhYAs8jA5PtGD0HS0OIaqyPiWFE8Gjhnuigi/ily\nC0Y3qrNkAnAP+ZY9HbgGWAqMJAOmtqhzD72t62ryTX9mRByT1I4m9tUNwNGIODjQDbH/B49UmLXu\nRCWgKBwEdhcBRbmss/j/FjKPYFeRvNctqZvMNag3UtGTH6q2byKnU8o2AsMljS2V/Vj5JyLOAIeA\n7aWySsfTSQNF3sAa8g3/2aZb3rxB5PTFzIjYHBGfk/P/jxeJnpNKz/KopEdbvUC9e5C0rlT39roV\nnG8ZGfBtrFTVapv6gaidW2LWLzxSYda601XbUaesErQPJ3MKbiM7yrLevB0f78U50Fy7ocHLRqkz\nHgdMqXrDP0BVUCJpMDnCcIDm7Qd+q6p7B9lJjgU2ARNL+1p6E29wD0+ROTJQ+/nUMxmYJmlu5TLA\nIEl/A09HxPst1HUU6KhRPoJMKG7WLqBD0iiPVtiF4KDCrP9tIUcqRpXeYttpB5nVXzYJ6I6Ife28\nUKkzngBMjogjVYd8A4yQdGspr+J+soP9roVLbQRmSLoiIk4UZTeSQdm+iDgF/NIf9xAR+3tTL5lP\nMri0PR2YB9xNJs+2oguYWqP8djJQaNZacnpnHvBi9U5JHRHRSpBi1iMHFWaN9WkYOyJ+lrSaTIKc\nQwYZncAUMs/hsz62713geUlLgHfI6ZD5wJt9rPccRWf8MbkkcxowRNKoYvfhiDgdETslfQEskzSL\nXPWyBPiwlE+CpJuBoeQIxnBJEwEiYltxyGoy12SFpPnAtcDrwPIioKjXxivJPILKbzahqPtwROxt\n5h56qPt64CpyueawSpuBnyKXZHZVHX8HcCYidtSrswdLgdmS3iaXsZ4q2vtI8bdmE6sLImKfpBeA\nJZI6gA+A3eRoz2NANzC3+jyz3nJQYdZYO+aknyA7yTfIZX5/At8Cn/S1LRHxu6QHyQz/reRSzGXA\nwp7Oa6GsYgxnO7Stxd/KnP1k4OuibCYZ3KwnRxbWcu7qFYB1wPjS9painsHFPR2XNJUMSDaR+R8f\nAa/00D7IVTMbOPs9h0pgtRJ4soV7qOU9Mg+movIBrOuAPQ3a1ZKI+FXSveRv+CUZnO0EZkREvW97\n1PztImKppC5gDrlqaBgZWHwKnPeBLbO+UIRzeMzMBoKk+4CvgJEXwzSEpBVAR0RUT6eZNcWrP8zM\nBt5eSasG6uKV1TTkKJNZr3mkwsxsgBTfyRhTbB6LAfok/MXSDvvvc1BhZmZmbeHpDzMzM2sLBxVm\nZmbWFg4qzMzMrC0cVJiZmVlbOKgwMzOztnBQYWZmZm3hoMLMzMzawkGFmZmZtcW/89K+Barv2C0A\nAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7ffbab05a4d0>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# time = Time.now()\n",
"time = Time('2016-12-14 19:00:00')\n",
"start = Time('2016-12-14 18:26:00').datetime\n",
"end = Time('2016-12-14 21:16:00').datetime\n",
"plot_airmass(exo2, observer, time)\n",
"#plot_airmass(exo2, observer, time)\n",
"#plot_airmass(exo3, observer, time)\n",
"p = plt.axvspan(start, end, facecolor='g', alpha=0.5)\n",
"\n",
"plt.legend(loc=1, bbox_to_anchor=(1, 1)) \n",
"plt.show() "
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.12"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment