Created
February 20, 2018 21:56
-
-
Save pllim/1c51018df5bc3c3285c4aa258f3530a9 to your computer and use it in GitHub Desktop.
synphot: Flat spectrum is not really flat in other units
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": 16, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"from synphot import SourceSpectrum, SpectralElement, Observation, units, models\n", | |
"from astropy import units as u" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 2, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"sp = SourceSpectrum(models.ConstFlux1D, amplitude=0.004*units.FLAM)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 3, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"bigbox = SpectralElement(models.Box1D, x_0=5000*u.AA, width=1000*u.AA)\n", | |
"smallbox = SpectralElement(models.Box1D, x_0=5000*u.AA, width=10*u.AA)" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"This is to show that flat spectrum is not really flat in other units." | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 6, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAEKCAYAAAAiizNaAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3X+UXlV97/H3hwkEVuWHwqiYgEkl\nLBusoH0a8VflgkpAVmKv3GWoWlQwbYWr1Vokq729QntvSXuveG0Bi0Ab8UeSC7WOiCLc+BMhyQQC\nEjQ6DSAjtEQDUaoNnfC5f5w95DA8M89MMs8ZMnxeaz0r5+yz9z57+2C+2efsZ2/ZJiIioin7THUD\nIiLimSWBJyIiGpXAExERjUrgiYiIRiXwREREoxJ4IiKiUQk8ERHRqASeiIhoVAJPREQ0asZUN+Dp\n6LDDDvOcOXOmuhkREXuVDRs2/MR2b6d8CTxtzJkzh/7+/qluRkTEXkXSfePJl0dtERHRqASeiIho\nVAJPREQ0KoEnIiIalcATERGN6mrgkbRQ0mZJA5LOb3N9pqRV5fpaSXNq15aV9M2STh5RrkfS7ZKu\nq6V9puS9S9JVkvYt6ZL08VLXnZJe3r0eR0REJ10LPJJ6gEuAU4D5wBmS5o/IdhbwsO2jgIuB5aXs\nfGAJcAywELi01Dfs/cD3RtT1GeDFwK8DBwBnl/RTgHnlsxS4bDL6FxERu6ebI54FwIDtLbYfA1YC\ni0fkWQysKMfXACdJUklfaXuH7XuAgVIfkmYDbwKuqFdk+3oXwDpgdu0enyqXbgUOkXT4ZHc2IiLG\np5uBZxZwf+18sKS1zWN7CNgOHNqh7MeA84DH2920PGJ7B/CVCbQDSUsl9Uvq37p1a6e+RUTEbupm\n4FGbNI8zT9t0SacBD9neMMZ9LwW+aftbE2gHti+33bLd6u3tuOJDRETspm4GnkHgiNr5bOCB0fJI\nmgEcDGwbo+yrgUWS7qV6dHeipE8PZ5L034Fe4IMTbEdERDSkm4FnPTBP0lxJ+1FNFugbkacPOLMc\nnw6sKe9o+oAlZdbbXKqJAetsL7M92/acUt8a228HkHQ2cDJwhu3HR9zjd8vstuOB7bYf7EqPIyKi\no64tEmp7SNK5wA1AD3CV7U2SLgT6bfcBVwJXSxqgGuksKWU3SVoN3A0MAefY3tnhlp8A7gNuqeYn\n8I+2LwSuB06lmqDwC+Bdk9zViIiYAFUDjKhrtVrO6tQRERMjaYPtVqd8WbkgIiIalcATERGNSuCJ\niIhGJfBERESjEngiIqJRCTwREdGoBJ6IiGhUAk9ERDQqgSciIhqVwBMREY1K4ImIiEYl8ERERKMS\neCIiolEJPBER0agEnoiIaFRXA4+khZI2SxqQdH6b6zMlrSrX10qaU7u2rKRvlnTyiHI9km6XdF0t\n7dyS35IOq6WfIGm7pI3l82fd6W1ERIxH13YgldQDXAK8ARgE1kvqs313LdtZwMO2j5K0BFgOvFXS\nfKrdSI8BXgDcJOno2i6k7we+BxxUq+tm4Drg622a8y3bp01e7yIiYnd1c8SzABiwvcX2Y8BKYPGI\nPIuBFeX4GuAkVftWLwZW2t5h+x6qbasXAEiaDbwJuKJeke3bbd/brc5ERMTk6GbgmQXcXzsfLGlt\n89geArYDh3Yo+zHgPODxCbTllZLukPRlScdMoFxEREyybgYetUnzOPO0TZd0GvCQ7Q0TaMdtwAtt\nHwv8DfBPbRsrLZXUL6l/69atE6g+IiImopuBZxA4onY+G3hgtDySZgAHA9vGKPtqYJGke6ke3Z0o\n6dNjNcL2z2w/Wo6vB/atTz6o5bvcdst2q7e3d9ydjIiIielm4FkPzJM0V9J+VJMF+kbk6QPOLMen\nA2tsu6QvKbPe5gLzgHW2l9mebXtOqW+N7beP1QhJzy/vjZC0gKrPP52cLkZExER1bVab7SFJ5wI3\nAD3AVbY3SboQ6LfdB1wJXC1pgGqks6SU3SRpNXA3MAScU5vR1pak91G9+3k+cKek622fTRXQ/kDS\nEPBLYEkJbhERMQWUv4OfqtVqub+/f6qbERGxV5G0wXarU76sXBAREY1K4ImIiEYl8ERERKMSeCIi\nolEJPBER0agEnoiIaFQCT0RENCqBJyIiGpXAExERjUrgiYiIRiXwREREoxJ4IiKiUQk8ERHRqASe\niIhoVAJPREQ0KoEnIiIa1dXAI2mhpM2SBiSd3+b6TEmryvW1kubUri0r6ZslnTyiXI+k2yVdV0s7\nt+S3pMNq6ZL08XLtTkkv705vIyJiPLoWeCT1AJcApwDzgTMkzR+R7SzgYdtHARcDy0vZ+VTbYB8D\nLAQuLfUNez/wvRF13Qy8HrhvRPopwLzyWQpctmc9i4iIPdHNEc8CYMD2FtuPASuBxSPyLAZWlONr\ngJMkqaSvtL3D9j3AQKkPSbOBNwFX1Cuyfbvte9u0YzHwKVduBQ6RdPik9DAiIiasm4FnFnB/7Xyw\npLXNY3sI2A4c2qHsx4DzgMcnsR0REdGQbgYetUnzOPO0TZd0GvCQ7Q2T3A4kLZXUL6l/69atE6g+\nIiImopuBZxA4onY+G3hgtDySZgAHA9vGKPtqYJGke6ke3Z0o6dOT0A5sX267ZbvV29vbocqIiNhd\n3Qw864F5kuZK2o9qskDfiDx9wJnl+HRgjW2X9CVl1ttcqokB62wvsz3b9pxS3xrbb+/Qjj7gd8vs\ntuOB7bYfnJQeRkTEhHUt8JR3NucCN1DNQFtte5OkCyUtKtmuBA6VNAB8EDi/lN0ErAbuBr4CnGN7\n51j3k/Q+SYNUI5o7JQ1PPrge2EI1QeGTwHsnsZsRETFBqgYYUddqtdzf3z/VzYiI2KtI2mC71Slf\nVi6IiIhGJfBERESjEngiIqJRCTwREdGoBJ6IiGhUAk9ERDQqgSciIhqVwBMREY1K4ImIiEYl8ERE\nRKMSeCIiolEJPBER0agEnoiIaFQCT0RENCqBJyIiGpXAExERjepq4JG0UNJmSQOSzm9zfaakVeX6\nWklzateWlfTNkk4eUa5H0u2SrqulzS11/LDUuV9Jf6ekrZI2ls/Z3etxRER0MqNTBkn7A6cBrwVe\nAPwSuAv4UtmierRyPcAlwBuAQWC9pD7bd9eynQU8bPsoSUuA5cBbJc0HlgDHlHveJOno2vbX76fa\nTvugWl3LgYttr5T0iVL3ZeXaKtvnduprRER035gjHkkfAW4GXgmsBf4OWA0MARdJulHSS0cpvgAY\nsL3F9mPASmDxiDyLgRXl+BrgJEkq6Stt77B9DzBQ6kPSbOBNwBW1dgo4sdRBqfPNY3c9IiKmQqcR\nz3rbHxnl2kclPRc4cpTrs4D7a+eDwCtGy2N7SNJ24NCSfuuIsrPK8ceA84ADa9cPBR6xPdQmP8Bb\nJP0W8APgA7br7QJA0lJgKcCRR47WpYiI2FNjjnhsf6nD9Yds949yWe2KjDNP23RJpwEP2d4wgXt9\nEZhj+6XATewaYT05s3257ZbtVm9vb7ssERExCcYc8UjqG+u67UVjXB4EjqidzwYeGCXPoKQZwMHA\ntjHKLgIWSToV2B84SNKngXcAh0iaUUY9T9zL9k9r9XyS6l1QRERMkU6P2l5J9Sjsc1TveNqNLEaz\nHpgnaS7wY6rJAr8zIk8fcCZwC3A6sMa2S8D7rKSPUk0umAess30LsAxA0gnAh2y/vZx/rdSxstT5\nhZJ+uO0Hy/0WUU1KiIiIKdIp8DyfalbaGVRB40vA58aazTasvLM5F7gB6AGusr1J0oVAv+0+4Erg\nakkDVCOdJaXsJkmrgbupJjKcU5vRNpoPAysl/QVwe6kb4H2SFpV6tgHv7NT2iIjoHtkjX7uMklGa\nSRWA/hq40PbfdLNhU6nVarm/f7RXVxER0Y6kDbZbnfKN53c8M6mmL58BzAE+DvzjnjYwIiKemTpN\nLlgBvAT4MnCB7bsaaVVERExbnUY87wD+DTia6l3JcLoA2z5otIIRERHtjBl4bGcR0YiImFQTDiyS\nfkXS2ySN+ePSiIiIdsYVeCTtJ+nNZYrzg8DrgU90tWURETEtdZpcMPwbnpOBrwFXAwtsv6uBtkVE\nxDTUaXLBDcC3gNeUVaKR9H+63qqIiJi2OgWe36BaTeAmSVuolqPp6XqrIiJi2uq0OvXttj9s+0XA\nR4CXAftJ+nLZRiAiImJCxj2rzfbNZRfPWVR74pzetVZFRMS0NeHp1LYft30D8OIutCciIqa5/EA0\nIiIalcATERGN6vQ7nr/hqdtVQ7VW2yFdaVFERExrnUY8/cCGNp9+4L92qlzSQkmbJQ1IOr/N9ZmS\nVpXrayXNqV1bVtI3Szp5RLkeSbdLuq6WNrfU8cNS536d7hEREc3r9Duez9ge2p2KJfUAl1DtYDoI\nrJfUZ/vuWrazgIdtHyVpCbAceKuk+VS/HzqGauvrmyQdXduF9P1UW1jXV8deDlxse6WkT5S6Lxvt\nHrvTp4iI2HOdRjzrhg/KY7eJWAAM2N5i+zGqH58uHpFnMbCiHF8DnKRq74XFwErbO8qKCQOlPiTN\nptqY7opa2wScWOqg1PnmDveIiIgp0GnEU/8L+tUTrHsWcH/tfBB4xWh5bA9J2g4cWtJvHVF2Vjn+\nGHAecGDt+qHAI7XRWT3/aPf4yQT709GG+7ZxxbfumexqIyIa88Zjnsdvv2x2V+/RKfC0m1gwXu1G\nFSPrGy1P23RJpwEP2d4g6YRx3ms87aCsxLAU4Mgjj2xTpLNHd+zkn7c+ultlIyKeDn766LO7fo9O\ngefFku6k+sv7ReUYdu1A+tIxyg4CR9TOZwMPjJJnUNIM4GBg2xhlFwGLJJ0K7A8cJOnTVDulHiJp\nRhn11O812j2exPblwOUArVZrtwLu647u5XVHv253ikZEPGN0Cjy/tgd1rwfmSZoL/JhqssDvjMjT\nB5wJ3EK1BM8a25bUB3xW0kepJhfMA9bZvgVYBlBGPB+y/fZy/rVSx8pS5xfGusce9CsiIvZAp8Dz\no05/SUtSuzzlfcq5VFsr9ABX2d4k6UKg33YfcCVwtaQBqlHIklJ2U9l07m5gCDinNqNtNB8GVkr6\nC+D2Ujej3SMiIqaGxoorkr4OXAt8wfaPaun7Aa+hGkl8zfY/dLeZzWq1Wu7v75/qZkRE7FUkbbDd\n6pSv04hnIfBu4HPlkdkjVO9WeoCvUv1uZuOeNjYiIp45xgw8tv8duBS4VNK+wGHAL20/0kTjIiJi\n+uk04nmC7f8AHuxiWyIi4hkgq1NHRESjEngiIqJR4wo8ZdHOkWknTHprIiJi2hvviGe1pA+rckBZ\nMPQvu9mwiIiYnsYbeF5BtezMd6hWJHiAiS8aGhERMe7A8x/AL4EDqH7Hc4/tx7vWqoiImLbGG3jW\nUwWe36RaseAMSdeMXSQiIuKpxvs7nrNsD68h8y/AYknv6FKbIiJiGhtv4HlI0shNar4x2Y2JiIjp\nb7yB50vs2qBtf2AusBk4pkvtioiIaWpcgcf2r9fPJb0c+L2utCgiIqa13Vq5wPZtVBMNIiIiJmRc\nIx5JH6yd7gO8HNjalRZFRMS0Nt4Rz4G1z0yqdz6LOxWStFDSZkkDks5vc32mpFXl+lpJc2rXlpX0\nzZJOLmn7S1on6Q5JmyRdUMt/oqTbJN0laYWkGSX9BEnbJW0snz8bZ58jIqILxvuO54LOuZ5MUg9w\nCfAGYBBYL6nP9t21bGcBD9s+StISYDnw1rI23BKqyQsvAG6SdDSwAzjR9qNlf6BvS/oysA5YAZxk\n+wdle+0z2bX99bdsnzbRPkRExOQbM/BI+iLVbLa2bC8ao/gCYMD2llLXSqpRUj3wLAY+Uo6vAf5W\nkkr6Sts7gHskDQALbN8CPFry71s+Bg4Fdtj+Qbl2I7CMXYEnIiKeJjqNeP7XHtQ9C7i/dj5IteZb\n2zy2hyRtpwois4BbR5SdBU+MpDYARwGX2F5bgtW+klrlh66nU60tN+yVku6gWmPuQ7Y37UG/IiJi\nD3QKPPfY/tFu1q02aSNHT6PlGbWs7Z3AcZIOAT4v6SW27yqP6i6WNBP4KjBUyt0GvLA8njsV+Cdg\n3lMaKy0FlgIceeTI38pGRMRk6TS54J+GDyRdO8G6B3nyqGM21YijbZ4yGeBgYNt4ytp+BPg6sLCc\n32L7tbYXAN8EfljSf2b70XJ8PdXI6LCRjbV9ue2W7VZvb+8EuxoREePVKfDURx6/OsG61wPzJM2V\ntB/VZIG+EXn6qCYBQPV4bI1tl/QlZdbbXKoRyjpJvWWkg6QDgNcD3y/nzy1/zgQ+DHyinD+/PIpD\n0oLS559OsC8RETFJOj1q8yjHHZV3NucCNwA9wFW2N5UZZ/22+6he/l9dJg9sowpOlHyrqSYiDAHn\n2N4p6XBgRXnPsw+w2vZ15ZZ/LOm0kn6Z7TUl/XTgDyQNUa2wvaQEt4iImAIa6+9gSTuBf6Ma+RwA\n/GL4EmDbB3W9hVOg1Wq5v7+/c8aIiHiCpA22W53yjTnisd0zeU2KiIjYzbXaIiIidlcCT0RENCqB\nJyIiGpXAExERjUrgiYiIRiXwREREoxJ4IiKiUQk8ERHRqASeiIhoVAJPREQ0KoEnIiIalcATERGN\nSuCJiIhGJfBERESjEngiIqJRXQ08khZK2ixpQNL5ba7PlLSqXF8raU7t2rKSvlnSySVtf0nrJN0h\naZOkC2r5T5R0m6S7JK2QNKOkS9LHS113Snp5N/scERFj61rgKdtTXwKcAswHzpA0f0S2s4CHbR8F\nXAwsL2XnU22DfQywELi01LcDONH2scBxwEJJx0vaB1hBta31S4D7gDPLPU4B5pXPUuCyLnU5IiLG\noZsjngXAgO0tth8DVgKLR+RZTBUwAK4BTpKkkr7S9g7b9wADwAJXHi359y0fA4cCO2z/oFy7EXhL\n7R6fKmVvBQ6RdPik9zYiIsalm4FnFnB/7XywpLXNY3sI2E4VREYtK6lH0kbgIeBG22uBnwD7Shre\n6/t04IgJtANJSyX1S+rfunXrBLsaERHj1c3AozZpHmeeUcva3mn7OGA2sEDSS2yb6tHcxZLWAT8H\nhibQDmxfbrtlu9Xb29u2QxERsee6GXgG2TXqgCpQPDBanjIZ4GBg23jK2n4E+DrVOyBs32L7tbYX\nAN8EfjiBdkREREO6GXjWA/MkzZW0H9WIpG9Enj52TQI4HVhTRi99wJIy620u1cSAdZJ6JR0CIOkA\n4PXA98v5c8ufM4EPA5+o3eN3y+y244Htth/sTpcjIqKTGd2q2PaQpHOBG4Ae4CrbmyRdCPTb7gOu\nBK6WNEA10llSym6StBq4m+qR2Tm2d5ZJASvKDLd9gNW2ryu3/GNJp5X0y2yvKenXA6dSTVD4BfCu\nbvU5IiI6UzXAiLpWq+X+/v6pbkZExF5F0gbbrU75snJBREQ0KoEnIiIalcATERGNSuCJiIhGJfBE\nRESjEngiIqJRCTwREdGoBJ6IiGhUAk9ERDQqgSciIhqVwBMREY1K4ImIiEYl8ERERKMSeCIiolEJ\nPBER0agEnoiIaFRXA4+khZI2SxqQdH6b6zMlrSrX10qaU7u2rKRvlnRySdtf0jpJd0jaJOmCWv6T\nJN0maaOkb0s6qqS/U9LWkr5R0tnd7HNERIyta4GnbE99CXAKMB84Q9L8EdnOAh62fRRwMbC8lJ1P\ntQ32McBC4NJS3w7gRNvHAscBCyUdX+q6DHib7eOAzwJ/WrvPKtvHlc8VXehuRESMUzdHPAuAAdtb\nbD8GrAQWj8izGFhRjq8BTpKkkr7S9g7b9wADwAJXHi359y2f4b27DRxUjg8GHuhGpyIiYs90M/DM\nAu6vnQ+WtLZ5bA8B24FDxyorqUfSRuAh4Ebba0ues4HrJQ0C7wAuqpV/i6Q7JV0j6Yh2jZW0VFK/\npP6tW7dOvLcRETEu3Qw8apPmceYZtaztneVx2mxggaSXlOsfAE61PRv4e+CjJf2LwBzbLwVuYtcI\n68mV25fbbtlu9fb2jtGtiIjYE90MPINAfXQxm6c+/noij6QZVI/Ito2nrO1HgK9TvefpBY6tjX5W\nAa8q+X5qe0dJ/yTwG3vUq4iI2CPdDDzrgXmS5kraj2qyQN+IPH3AmeX4dGCNbZf0JWXW21xgHrBO\nUq+kQwAkHQC8Hvg+8DBwsKSjS11vAL5X8h1eu9+i4fSIiJgaM7pVse0hSecCNwA9wFW2N0m6EOi3\n3QdcCVwtaYBqpLOklN0kaTVwNzAEnGN7ZwkiK8oMt32A1bavA5D0HuBaSY9TBaJ3l6a8T9KiUs82\n4J3d6nNERHSmaoARda1Wy/39/VPdjIiIvYqkDbZbnfJl5YKIiGhUAk9ERDQqgSciIhqVwBMREY1K\n4ImIiEYl8ERERKMSeCIiolEJPBER0agEnoiIaFQCT0RENCqBJyIiGpXAExERjUrgiYiIRiXwRERE\noxJ4IiKiUV0NPJIWStosaUDS+W2uz5S0qlxfK2lO7dqykr5Z0sklbX9J6yTdIWmTpAtq+U+SdJuk\njZK+LemoTveIiIjmdS3wlF1CLwFOAeYDZ0iaPyLbWcDDto8CLgaWl7LzqXYjPQZYCFxa6tsBnGj7\nWOA4YKGk40tdlwFvs30c8FngT8e6R0RETI1ujngWAAO2t9h+DFgJLB6RZzGwohxfA5wkSSV9pe0d\ntu8BBoAFrjxa8u9bPsNbqBo4qBwfDDzQ4R4RETEFZnSx7lnA/bXzQeAVo+WxPSRpO3BoSb91RNlZ\n8MRIagNwFHCJ7bUlz9nA9ZJ+CfwMGB4JjXaPn0xCHyMiYoK6OeJpN6rwOPOMWtb2zvI4bTawQNJL\nyvUPAKfang38PfDRCbQDSUsl9Uvq37p1a5siERExGboZeAaBI2rns9n1+OspeSTNoHpEtm08ZW0/\nAnyd6j1PL3BsbfSzCnhVh3s8ie3Lbbdst3p7eyfU0YiIGL9uBp71wDxJcyXtRzVZoG9Enj7gzHJ8\nOrDGtkv6kjIjbS4wD1gnqVfSIQCSDgBeD3wfeBg4WNLRpa43AN/rcI+IiJgCXXvHU96nnAvcAPQA\nV9neJOlCoN92H3AlcLWkAapRyJJSdpOk1cDdwBBwju2dkg4HVpT3PPsAq21fByDpPcC1kh6nCkTv\nLk1pe4+IiJgayj/+n6rVarm/v3+qmxERsVeRtMF2q1O+rFwQERGNSuCJiIhGJfBERESjEngiIqJR\nCTwREdGozGprQ9JW4L6pbsduOIxn3lJA6fP090zrL+y9fX6h7Y6/wE/gmUYk9Y9nKuN0kj5Pf8+0\n/sL073MetUVERKMSeCIiolEJPNPL5VPdgCmQPk9/z7T+wjTvc97xREREozLiiYiIRiXw7AUk9Ui6\nXdLwStyS9D8k/UDS9yS9r5b+cUkDku6U9PJaHWdK+mH5nDnavZ4u2vT5JEm3Sdoo6duSjirpMyWt\nKn1eK2lOrY5lJX2zpJOnpifjI+leSd8t/esvac+RdGP5zm6U9OySPi2+51H6/NeSvl/69fnhbVDK\ntbbfp6SFJW1A0vlT0Zfxatfn2rUPSbKkw8r5tPie27Kdz9P8A3wQ+CxwXTl/F/ApYJ9y/tzy56nA\nl6l2XT0eWFvSnwNsKX8+uxw/e6r7NcE+/wD4tXL8XuAfasefKMdLgFXleD5wBzATmAv8M9Az1f0a\no7/3AoeNSPsr4PxyfD6wfDp9z6P0+Y3AjHK8vNbntt9n+fwz8KvAfiXP/Knu20T6XNKPoNpC5r7h\n69Ple273yYjnaU7SbOBNwBW15D8ALrT9OIDth0r6YuBTrtwKHFL2MDoZuNH2NtsPAzcCCxvrxASN\n0mcDB5Xjg9m1I+1iYEU5vgY4SZJK+krbO2zfAwwAC7rd9klW79sK4M219L3+e27H9ldtD5XTW6l2\nH4bRv88FwIDtLbYfA1aWvHubi4HzqP47HzZtv+cEnqe/j1H9B/l4Le1FwFsl9Uv6sqR5JX0WcH8t\n32BJGy396apdn88Grpc0CLwDuKikP9G38hfWduBQ9r4+G/iqpA2Slpa059l+EKD8+dySPl2+53Z9\nrns31b/4YRr3WdIi4Me27xiRd7r0+Sm6tgNp7DlJpwEP2d4g6YTapZnAv9tuSfrPwFXAa6mG5CN5\njPSnnTH6/AHgVNtrJf0x8FGqYLTX97l4te0HJD0XuFHS98fIO237bPubAJL+hGr34c+UvKP1rd0/\nnveqPgN/QvWIcaTp8j0/RUY8T2+vBhZJupfqEcKJkj5N9S+ca0uezwMvLceDVM+Kh82meiQ1WvrT\nUbs+fwk41vbakmcV8Kpy/ETfJM2gegy3jb2rz9h+oPz5ENV3ugD41/JohfLn8CPV6fA9j9Znysvy\n04C3ubzUYPr2+XVU76zuKP/NzwZuk/R8pkmf25rql0z5jO8DnMCuF+0XAe+upa8vx2/iyS8j15X0\n5wD3UL2IfHY5fs5U92m8faYamf8EOLqknwVcW47P4cmTC1aX42N48svoLTxNJxcAvwIcWDv+DtUz\n+7/myZML/mq6fM9j9HkhcDfQOyJ/2++z/LexpaQNTy44Zqr7N5E+j8hzL7smF+z13/Nonzxq2ztd\nBHxG0geAR6keOQFcTzUTZgD4BdXsN2xvk/TnwPqS70Lb25pt8u6zPSTpPcC1kh4HHqZ6/g9wJXC1\npAGqkc6SUmaTpNVUf4kNAefY3tl868flecDnqzkRzAA+a/srktYDqyWdBfwI+C8l/3T4nkfr8wBV\ncLmxXLvV9u+P9X1KOpdqRlgPcJXtTc13Z1za9nmM/NPhe24rKxdERESj8o4nIiIalcATERGNSuCJ\niIhGJfBERESjEngiIqJRCTyx15N0saQ/rJ3fIOmK2vn/lvTBSb7no5NZX6nzOEmn1s4/IulD4ygn\nSWskHVRL++2y0vGLu9DOOZJ+Z7LrrdV/rqR3dav+mHoJPDEdfIeykoGkfYDDqH5wOOxVwM1T0K6J\nOo7qdxsTdSpwh+2f1dLOAL5N+V3TJJsDtA08ZfWIPXUV8L5JqCeephJ4Yjq4mV1L6BwD3AX8XNKz\nJc0Efg24XdKzJP0/Vfv6fFfSYgBJyyW9d7iyMtL4o3L8x5LWl/1QLmh383Z5yqjge5I+KWmTpK9K\nOqBc+82S9xZV+8/cJWk/4EIHnAq2AAAEUklEQVSqxV83SnprqX6+pK9L2qKy71IbbwO+UGvPs6iW\nHjqLWuCRdEKp6xpVe958pqzkjaRTS9q3Ve0BM7wP0utKezaq2h/pQKofML+2pH1A0jsl/V9JX6Ra\nAFO1fn13uC/l/t+QtFrVXlIXSXqbpHUl34sAbP8CuFfS3raaeIzXVC+dkE8+k/GhWmrkSOD3gN8H\n/pxqJPBq4JslzwzgoHJ8GNUvwgW8DPhGra67S11vBC4vefahWr7nt0qeR8ufbfNQjQqGgONKvtXA\n28vxXcCryvFFwF3l+J3A39ba8RGq0dzM0t6fAvu26ft9lKVYyvnbgSvL8XeAl5fjE6hW755d2noL\n8Bpgf6rVjueWfJ9j1/JMX6Ra2BLgWeV/wxOGr9faPUhZtgV4C9VS/T1Uv9b/EXB4KfdIOZ4J/Bi4\noJR5P/CxWp1/AvzRVP93lU93PhnxxHQxPOp5FdVfqLfUzr9T8gj4n5LuBG6iWkr+ebZvB54r6QWS\njgUetv0jqqDyRuB24DbgxcA8nmysPPfY3liONwBzVO2oeaDt4TZ9tkO/vuRqD5qfUC0S+rw2eZ5j\n++e18zOoFlil/HlG7do624Ou9nLaSBUgXwxscbXPDVSBZ9jNwEfLaOsQ79orZ6QbvWvZltcAn7O9\n0/a/At8AfrNcW2/7Qds7qDZw+2pJ/25py7CHgBeMcq/Yy2Wttpguht/z/DrViOJ+4I+An1G9M4Dq\nkVQv8Bu2/6OsBrx/uXYNcDrwfHb9pS3gL23/3Rj3bZtH1RbcO2pJO4EDaL+k/VhG1tHu/7NDkvax\n/bikQ4ETgZdIMtWow5LOG6O+Udtk+yJVq4OfCtwq6fWjZP232vFYfazf//Ha+eM8uW/7A78co57Y\ni2XEE9PFzVRL6W8r/9LeBhwCvJJq9APVlgkPlaDzn4AX1sqvpHofcjpVEIJq4cl3l3cmSJqlah+V\nuvHkeYKrHSN/Lun4klR/+f9z4MCJdLrYTLX1M6X9n7L9QttzbB9BtXrxa8Yo/33gV0uwBBh+v4Sk\nF9n+ru3lQD/V6KhTO79J9a6qR1Iv1aPHdRPs09FU/4CIaSiBJ6aL71K9B7l1RNr28pgKqk3FWpL6\nqUY/T2y25mpF4wOpdoIc3vXzq1SPwm6R9F2qgPSkv3DHk6eNs4DLJd1CNTrYXtK/RjWZoD65YDy+\nRPX+BKrHap8fcf1aRpmFVvrwS+C9wFckfRv411qb/rBMEriDagTyZeBOqlHWHapWSB/p8yXPHcAa\n4Dzb/zKB/kD1bu6mCZaJvURWp45omKRn2X60HJ8PHG77/XtQ3+FUo5w37Gmbyiy3S4Af2r54d+vb\nE5JeBnzQ9jum4v7RfRnxRDTvTWVUcxfVluV/sSeVlRHaJ1X7AelueI+kjcAmqkeSY73X6rbDgP82\nhfePLsuIJyIiGpURT0RENCqBJyIiGpXAExERjUrgiYiIRiXwREREoxJ4IiKiUf8fa67b4nWXFVoA\nAAAASUVORK5CYII=\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe44658550>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"sp.plot(bigbox.waveset, flux_unit=units.FLAM)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 9, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAERCAYAAACU1LsdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xd8VfX9x/HXhxX2HrLCHrKVgAOr\niKj8XLQWV+umom3VOupeiNZZq92W1la0VmWLKFYcqChDUCBh7w1hQyCEjM/vj3OI1zQhAXJzb5L3\n8/G4j5x77jfnfL7ccD/3+/2e8/2auyMiIgJQIdYBiIhI/FBSEBGRXEoKIiKSS0lBRERyKSmIiEgu\nJQUREclVKpOCmf3TzFLNLKWYjveBme02s8kFvP5HM0srjnOJiMSzUpkUgFeBQcV4vOeBa/J7wcyS\ngLrFeC4RkbhVKpOCu38O7IzcZ2btwm/8c83sCzPrfBTH+xjYl3e/mVUkSBj3Hm/MIiKlQaVYB1CM\nRgK3uPtyMzsF+Asw4DiPeSswyd03m9lxBygiEu/KRFIws5rA6cCYiA/vhPC1S4ER+fzaRnc//wjH\nbAZcBvQv1mBFROJYmUgKBN1gu929V94X3H08MP4YjnkS0B5YESaa6ma2wt3bH1ekIiJxrFSOKeTl\n7nuB1WZ2GYAFeh7nMd9z9xPcvbW7twYOKCGISFlXKpOCmb0JzAA6mdkGMxsK/BQYambzgYXA4KM4\n3hfAGOCc8HgFdiuJiJRlpqmzRUTksFLZUhARkegodQPNDRs29NatW8c6DBGRUmXu3Lnb3b1RYeVK\nXVJo3bo1c+bMiXUYIiKlipmtLUo5dR+JiEguJQUREcmlpCAiIrmUFEREJJeSgoiI5FJSEBGRXEoK\nIiKSS0lBRCTOHczM5un3F7Nxd3rUz1Xqbl4TESlPlmzZy6/enMfSrftoUb8615zaKqrnU1IQEYlD\n7s6rX63h6SlLqF21Mq/e0If+nRpH/bxRTwrhOsdzCFY6uyjPawnAa0BvYAdwhbuviXZMIiLxLHXf\nQe4Zs4DPlm3jnM6NeXZIDxrWTCiRc5dES+FXwGKgdj6vDQV2uXt7M7sSeBa4ogRiEhGJSx8v3sq9\nYxeQlpHFE4O7cvWprSjJNeKjOtBsZi2AC4F/FFBkMDAq3B5LsMhNydVeRCROpB/K5uGJyQwdNYdG\ntRJ497YzuOa01iWaECD6LYWXgHuBWgW83hxYD+DuWWa2B2gAbI9yXCIicSNl4x5+9da3rNy2n5+d\n0YZ7BnUioVLFmMQStaRgZhcBqe4+18z6F1Qsn33/sxScmQ0DhgEkJiYWW4wiIrGUneOM/HwVv5u6\nlAY1EnjjZ6fQr33DmMYUzZZCP+ASM7sAqArUNrN/u/vVEWU2AC2BDWZWCagD7Mx7IHcfCYwESEpK\n0vqhIlLqbdydzl1vz2PW6p1c0P0EnvpRd+pWrxLrsKKXFNz9AeABgLCl8Os8CQFgEnAdMAMYAnzi\nWjRaRMq4d+Zt5OGJKeTkOM8P6cGQ3i1KfOygICV+n4KZjQDmuPsk4BXgdTNbQdBCuLKk4xERKSl7\nD2by6MQUJs7bxMmJdXnxil60alAj1mF9T4kkBXefBkwLtx+N2H8QuKwkYhARiaXZq3dy59vz2LL3\nIHcM7MCtZ7enUsX4m2lIdzSLiERRZnYOL320jL9OW0mLetUZc8tpnJxYL9ZhFUhJQUQkSlZtS+OO\nt+exYMMeLk9qwaMXd6VmQnx/7MZ3dCIipZC78+bs9TwxeREJlSvw8tUnM6hb01iHVSRKCiIixWhH\nWgb3jUvmo8VbOaN9Q357WU9OqFM11mEVmZKCiEgx+XRpKveMWcDe9EwevvBEbuzXhgoV4uNS06JS\nUhAROU6HF8EZNWMtnZrU4vWhfTmxaX5zgMY/JQURkeOwcNMe7nhrHstT07ihX2vuG9SZqpVjM29R\ncVBSEBE5Bjk5zj+mr+L5/y6lXvUqvHZjX87s2CjWYR03JQURkaO0aXc6d4+ez4xVOzi/axOevrQH\n9WvEft6i4qCkICJyFCYv2MSD45PJynGe/XF3Lk9qGTfzFhUHJQURkSLYdzCT4ZMWMe6bDfRsWZeX\nruhFm4bxNW9RcVBSEBEpxKxVO7hr9Hw270nn9gHtue2cDlSOw3mLioOSgohIATKysvnd1GWM/HwV\nifWrM+aW0+ndKn7nLSoOSgoiIvlYsmUvd7w1jyVb9nFV30QevvBEasT5vEXFoezXUETkKOTkOK9M\nX83z/11K7WqVeOW6JM45sUmswyoxSgoiIqGNu9O5e/Q8Zq7ayXldmvD0pd1pUDMh1mGVKCUFESn3\n3J2J8zby6MSF5Ljz3JAeXBZHS2SWJCUFESnXdu0/xMMTU3gveTN9Wtfjd5f3omX96rEOK2aUFESk\n3Pps2TbuGTOfXQcOce+gTtx8ZjsqlrJZTYubkoKIlDvph7J5espiXpuxlg6Na/LP6/vQrXmdWIcV\nF5QURKRcmb9+N3e+PY9V2/cz9Iw23HN+p1I9q2lxi1pSMLOqwOdAQniese7+WJ4yicAooC5QEbjf\n3d+PVkwiUn5lZefwl2kr+cPHy2lUK4E3fnYK/do3jHVYcSeaLYUMYIC7p5lZZWC6mU1x95kRZR4G\nRrv7X82sC/A+0DqKMYlIObR6+37ufHse89bv5oe9mvH44G7UqVY51mHFpaglBXd3IC18Wjl8eN5i\nwOHlieoAm6IVj4iUP+7Of2av48nJi6lc0fjjVSdxcc9msQ4rrkV1TMHMKgJzgfbAn919Vp4iw4EP\nzew2oAYwsIDjDAOGASQmJkYtXhEpO1L3HeS+sQv4dOk2ftChIc8P6ckJdarGOqy4F9Vp/tw92917\nAS2AvmbWLU+Rq4BX3b0FcAHwupn9T0zuPtLdk9w9qVGj0r+ykYhE1wcpWzj/xc/5auUOhl/chVE3\n9FVCKKISufrI3Xeb2TRgEJAS8dLQcB/uPiMcnG4IpJZEXCJStuw7mMmIdxcxZu4GujWvzUtX9KJ9\n41qxDqtUiebVR42AzDAhVCPoGno2T7F1wDnAq2Z2IlAV2BatmESk7Jq9eid3jZ7Hpt3p3DagPbcN\n6ECVSmVzzYNoimZLoSkwKhxXqEBwldFkMxsBzHH3ScDdwN/N7E6CQefrwwFqEZEiOZiZzYtTlzHy\ni/Kz5kE0RfPqowXASfnsfzRiexHQL1oxiEjZlrJxD3eNnseyrWnlas2DaNK/noiUOlnZOfx12kp+\n//Fy6teowr+u78PZnRvHOqwyQUlBREqVldvSuHv0fOat383FPZvxxOCu1K1eJdZhlRlKCiJSKuTk\nOK/NWMMzHyyhauWKuhEtSpQURCTubdydzj1j5vPVyh2c3akRz/64B41r676DaFBSEJG45e6M+2Yj\nj08KVkR7+tLuXNmnZblcEa2kKCmISFzanpbBA+OTmbpoK31b1+eFy3uW6xXRSoqSgojEnQ9StvDQ\nhGT2ZWTx0AUncuMZbcr9imglRUlBROLGnvRMHp+0kPHfbqRb89q8eXkvOjbRNBUlSUlBROLC9OXb\nuWfsfFL3ZXD7OR24bUB7KlfUNBUlTUlBRGIq/VA2z0xZzKgZa2nXqAbjf346PVvWjXVY5ZaSgojE\nzDfrdnH36Pms3r6fG/u14d5BWi851pQURKTEHcrK4fcfL+Ov01bStE41/nPTKZzeTuslxwMlBREp\nUYs37+Wu0fNZvHkvlye14JGLulCrqtZLjhdKCiJSIrJznJGfr+J3U5dSp1oV/nFtEgO7NIl1WJKH\nkoKIRN2a7fu5e8x85q7dxf91O4Hf/Kg79WtoErt4pKQgIlGTk+O8MWstT72/hMoVjZeu6MXgXs00\nTUUcU1IQkajYsOsA941bwJcrdvCDDg15bkgPmtapFuuwpBBKCiJSrNydt79ez5PvLcY1iV2po6Qg\nIsVm85507h+XzGfLtnFa2wY8N6SHJrErZZQUROS4uTvjv9nI8HcXkpXtjBjclatPaUUFTWJX6igp\niMhxSd13kAfHJ/PR4lT6tK7H80N60rphjViHJccoaknBzKoCnwMJ4XnGuvtj+ZS7HBgOODDf3X8S\nrZhEpPi4O5Pmb+KxSQtJP5TNwxeeyA39NMV1aRfNlkIGMMDd08ysMjDdzKa4+8zDBcysA/AA0M/d\nd5lZ4yjGIyLFZEdaBg9PTGFKyhZ6tazLC5f3pF2jmrEOS4pB1JKCuzuQFj6tHD48T7GbgD+7+67w\nd1KjFY+IFI8pyZt5eGIK+w5mcd+gztz0gzZU0hTXZUZUxxTMrCIwF2hP8OE/K0+RjmG5L4GKwHB3\n/yCf4wwDhgEkJiZGM2QRKcCu/Yd4bNJCJs3fRPfmdXjh8p5aAKcMimpScPdsoJeZ1QUmmFk3d0/J\nc/4OQH+gBfBFWGZ3nuOMBEYCJCUl5W1tiEiUfbRoKw9MSGbX/kPcdW5Hft6/nRbAKaNK5Oojd99t\nZtOAQUBkUtgAzHT3TGC1mS0lSBJfl0RcInJke9IzGfHuIsZ9s4HOJ9Ti1Rv60LVZnViHJVEUzauP\nGgGZYUKoBgwEns1TbCJwFfCqmTUk6E5aFa2YRKTopi1N5f5xyWxLy+D2Ae25dUAHqlRS66Csi2ZL\noSkwKhxXqACMdvfJZjYCmOPuk4D/AueZ2SIgG7jH3XdEMSYRKcS+g5k89f5i3py9ng6NazLy2t70\naKHlMcsLCy4SKj2SkpJ8zpw5sQ5DpEz6csV27h27gM170hl2ZjvuGNhBy2OWEWY2192TCiunO5pF\nhP0ZWTwzZQmvz1xL24Y1GHPL6fRuVS/WYUkMKCmIlHOzV+/k12Pms37XAW7s14Z7zu9EtSpqHZRX\nSgoi5dSBQ1k8/9+lvPrVGlrWq85bN53KKW0bxDosiTElBZFyaNaqHdw7bgFrdxzg2tNacd+gztRI\n0MeBFDEpmFkS8AOgGZBOcK/BR+6+M4qxiUgxO3Aoi+c+CFoHifWr8+ZNp3JaO7UO5DtHTApmdj1w\nO7CaYLqKpUBV4AzgPjNLAR5x93VRjlNEjtOMlTu4b9wC1u08wPWnt+beQZ2oXkWtA/m+wv4iahDM\nYJqe34tm1ovgDmQlBZE4FXllUasG1Xl7mMYOpGBHTAru/udCXp9XvOGISHH6asV27h23gI2703Vl\nkRRJYd1Hf8izy4HtwKfuPj1qUYnIcUnLyOKZKYv598x1tGlYg9E3n0af1vVjHZaUAoV1H83NZ199\n4Hkze9vdX4pCTCJyHKYv38594xawaU86PzujDXefp9aBFF1h3Uej8ttvZi8DXwFKCiJxIpizaAlv\nzl5H24Y1GHvLafRupdaBHJ1juvTA3dPNtA6rSLz4fNk27h+3gC17D3LzmW2589yOmrNIjslRJwUz\nqwRcQ7AWgojE0N6DmTz13mLe+no97RrVYOzPT+fkRM1ZJMeusIHmffzvusrpwGfAzdEKSkQKN21p\nKg+MT2br3oPccpZmNJXiUdiYghZgFYkze9IzeXLyIsbM3UCHxjX56y/60aul1juQ4lFYSyHxSK/r\nTmaRkvXpkqB1kLrvIL/o347bz1HrQIpXYWMK7+Wzz4FGQGNAf40iJWDPgUxGTA7WSu7YpCZ/u6Yf\nPdU6kCgorPuoe+RzM2sN3Eew3vJTUYtKRHJ9vHgrD4xPZsf+Q9x6dntuO6c9CZX0fUyio6izpHYA\nHgJOAV4Abnf3zGgGJlLe7T5wiBHvLmL8txvpfEItXrmuD91b1Il1WFLGFTam0I0gGXQFngOGunt2\nSQQmUp5NSd7MI+8sZNeBQ9w+oD23DuhAlUoVYh2WlAOFtRTmA+sJxhb6An0jb1pz99sL+kUzqwp8\nDiSE5xnr7o8VUHYIMAbo4+5zjqYCImXJtn0ZPDYphfeTt9C1WW1G3diHrs3UOpCSU1hSGMr/3qdQ\nVBnAAHdPM7PKwHQzm+LuMyMLmVktgjUbZh3jeURKPXdnwrcbGTF5EQcysrnn/E4MO7MtlSuqdSAl\nq7Ck8G93zzqWA7u7A2nh08rhI78E8wRB19Svj+U8IqXdpt3pPDghmWlLt3FyYl2eG9KD9o11i5DE\nRmFfQ2Yf3jCzPx7twc2sopnNA1KBqe4+K8/rJwEt3X1yIccZZmZzzGzOtm3bjjYMkbiUk+O8MWst\n5734ObNW7eTRi7ow5pbTlRAkpgprKUTOetfvaA8eDkr3MrO6wAQz6+buKQBmVgF4Ebi+CMcZCYwE\nSEpKOtbuLJG4sXbHfu4bt4CZq3ZyersGPHNpDxIbVI91WCKFJoVi+QB2991mNg0YBKSEu2sB3YBp\n4eD1CcAkM7tEg81SVmXnOP/6cjW//XAplStU4JlLu3NFn5Zo1mGJF4Ulhc5mtoCgxdAu3CZ87u7e\no6BfNLNGQGaYEKoR3PD27OHX3X0P0DCi/DTg10oIUlYt37qPe8ct4Nt1uzmnc2Oe/FE3mtapFuuw\nRL6nsKRw4nEcuykwyswqEoxdjHb3yWY2Apjj7pOO49gipUZmdg4vT1vJHz9ZQY2Eivz+yl5c0rOZ\nWgcSlwpLCuvCq4gKZGaWXxl3XwCclM/+R/M7jrv3LyQWkVInZeMe7h27gEWb93JRj6YMv6QrDWsm\nxDoskQIVlhQ+NbNxwDuRM6KaWRXgDOA64FPg1ahFKFIKHczM5g8fL+dvn6+ifo0q/O2a3pzf9YRY\nhyVSqMKSwiDgRuBNM2sD7AaqEsyO+iHworvPi26IIqXL3LU7uXfsAlZu289lvVvw8IVdqFO9cqzD\nEimSwmZJPQj8BfhLeFdyQyDd3XeXRHAipcmBQ1k898FSRs1YQ7M61Xjtxr6c2bFRrMMSOSpFXqM5\nnBV1cxRjESm1vlyxnfvHL2D9znSuPa0V9w7qTM2Eo14CXSTm9Fcrchz2HszkqfcW89bX62nTsAaj\nbz6Nvm3qxzoskWOmpCByjD5evJUHJySzbV8GN5/VljsHdtTSmFLqFXWRnS7uvijPvv7uPi0qUYnE\nse1pGTz+7iLenb+JzifU4u/XJtGjhZbGlLKhqC2F0Wb2OsFsplXDn0nAadEKTCTeuDvjv9nIE+8F\n01vfObAjP+/fTovfSJlS1KRwCsEUFV8RzFn0BscwQZ5IabV+5wEenJDMF8u307tVPZ65tDsdmmg2\nUyl7ipoUMoF0oBpBS2G1u+dELSqROJGd47z61Rp++9+lVDAYMbgrV5/SigoVNEWFlE1FTQpfA+8A\nfYAGwN/MbIi7D4laZCIxtmTLXu4bl8z89bsZ0LkxT/6wG83qagI7KduKmhSGRsxeugUYbGbXRCkm\nkZg6mJnNnz9dwV+nraROtcr84aqTuLhHU01gJ+VCUZNCqpkl5tn3WXEHIxJrX6/Zyf3jgikqLj25\nOY9c2IV6NarEOiyRElPUpPAewYI7RjCm0AZYCnSNUlwiJWrfwUye+2Apr89cS4t6mqJCyq8iJQV3\n7x753MxOBm6OSkQiJeyjRVt5eGIKqfsOMvSMNtx1bkdqaIoKKaeO6S/f3b8xsz7FHYxISdq2L4PH\n313I5AWb6dSkFi9f05teLXUTmpRvRb2j+a6IpxWAk4FtUYlIJMrcnXHfbOSJyYtIP5TNr8/ryLAz\ndROaCBS9pRB5l04WwRjDuOIPRyS6Im9C69O6Hk9f2oP2jWvGOiyRuFHUMYXHox2ISDRl5zj/+nI1\nL3y4jIoVjCd+2I2f9k3UTWgieRwxKZjZuwRXHeXL3S8p9ohEitnizXu5f9wC5m/Yw8ATG/PED7vR\ntI5uQhPJT2Ethd+WSBQiUXAwM5s/frKcv322irrVK/Onn5zEhd11E5rIkRSWFFa7+7pjObCZVQU+\nBxLC84x198fylLkL+BnBOMU24EZ3X3ss5xOJNHv1Tu4fv4BV2/YzpHcLHr7wROpW101oIoUpLClM\nJLjSCDMb5+4/PopjZwAD3D0tXN95uplNcfeZEWW+BZLc/YCZ/ZxgSu4rjuIcIt+zJz2TZ6Ys4c3Z\n62hZvxr/HnoKZ3RoGOuwREqNwpJCZDu77dEc2N0dSAufVg4fnqfMpxFPZwJXH805RA5zd95P3sLw\ndxeyc/8hhp3ZljsGdqB6Fd2EJnI0Cvsf4wVsF4mZVQTmAu2BP7v7rCMUHwpMKeA4w4BhAImJeadg\nkvJu4+50Hp2YwsdLUunevA7/ur4P3ZrXiXVYIqVSYUmhp5ntJWgxVAu3CZ+7u9c+0i+7ezbQy8zq\nAhPMrJu7p+QtZ2ZXE6zkdlYBxxkJjARISko66uQkZdPhtQ5e+HApAI9c1IXrTmtFpYq6CU3kWB0x\nKbh7saxC7u67zWwaMAj4XlIws4HAQ8BZ7p5RHOeTsi9l4x4eGJ9M8sY9DOjcmBGDu9KiXvVYhyVS\n6kWtw9XMGgGZYUKoBgwkWNIzssxJwN+AQe6eGq1YpOw4cCiLlz5azivTV1OvehVdZipSzKI5CtcU\nGBWOK1QARrv7ZDMbAcxx90nA80BNYEz4n3qdboiTgkxbmsrDE1PYsCudq/omcv+gztSpXjnWYYmU\nKVFLCu6+ADgpn/2PRmwPjNb5pezYti+DJyYvYtL8TbRvXJMxt5xGn9b1Yx2WSJmk6/Ukbrk7o+es\n56n3l5B+KJs7B3bklv5tSahULENdIpIPJQWJSyu3pfHg+GRmrd5J3zb1eepH3TWbqUgJUFKQuHIo\nK4eXP1vJnz5ZQdXKFXj2x925rHdLzWYqUkKUFCRuzFmzkwfGJ7M8NY2LezbjkYtOpHGtqrEOS6Rc\nUVKQmNuTnsmzHyzhP7PW0bxuNf51fR/O7tw41mGJlEtKChIz7s6UlC08NmkhO9Iy+NkZbbjz3I7U\nSNCfpUis6H+fxETkfEXdmtfmn9f1oXsLzVckEmtKClKisrJzeG3GWn774VLc4eELT+T601trviKR\nOKGkICVmwYbdPDghmZSNe+nfqRFPDO5Gy/qar0gknigpSNTtO5jJCx8u47UZa2hQM0HzFYnEMSUF\niZrDA8mPv7uQ1H0ZXHtqK+4+vxO1q2q+IpF4paQgUbF+5wEefSeFT5duo2uz2oy8JomeLevGOiwR\nKYSSghSrzOwc/vHFan7/8TIqmGnhG5FSRklBis3ctTt5cHwKS7fu47wuTRh+SVea1a0W67BE5Cgo\nKchx233gEM9+sIQ3Z6+nWZ2q/P3aJM7t0iTWYYnIMVBSkGPm7kyct5EnJy9md3omN/2gDXcM1B3J\nIqWZ/vfKMVm1LY1H3knhyxU76NWyLq//qDtdmtWOdVgicpyUFOSoZGRl8/K0Vfx52goSKlXgiR92\n4yd9E6moqa1FygQlBSmyr1Zu5+GJKazatl9TW4uUUUoKUqgdaRn85v3FjP9mI4n1qzPqxr6c1bFR\nrMMSkShQUpAC5eQ4Y+au5+kpS9ifkcWtZ7fn1gHtqVpZaySLlFVRSwpmVhX4HEgIzzPW3R/LUyYB\neA3oDewArnD3NdGKSYpu2dZ9PDQhma/X7KJv6/r85kfd6NCkVqzDEpEoi2ZLIQMY4O5pZlYZmG5m\nU9x9ZkSZocAud29vZlcCzwJXRDEmKUT6oWz++MlyRn6+ippVK/Hcj3swpHcLrZEsUk5ELSm4uwNp\n4dPK4cPzFBsMDA+3xwJ/MjMLf1dK2CdLtvLoOwvZsCudH5/cggcv6EyDmgmxDktESlBUxxTMrCIw\nF2gP/NndZ+Up0hxYD+DuWWa2B2gAbM9znGHAMIDExMRohlwubdydzuOTFvLhoq20b1yTN286ldPa\nNYh1WCISA1FNCu6eDfQys7rABDPr5u4pEUXy65P4n1aCu48ERgIkJSWpFVFMMrNzeGX6an7/0XIA\n7hvUmaFntKFKJU1eJ1JelcjVR+6+28ymAYOAyKSwAWgJbDCzSkAdYGdJxFTezVq1g4cnprA8NY1z\nuzThsYu70KKeVkETKe+iefVRIyAzTAjVgIEEA8mRJgHXATOAIcAnGk+Iru1pGTwV3nPQol41/nFt\nEgM1eZ2IhKLZUmgKjArHFSoAo919spmNAOa4+yTgFeB1M1tB0EK4MorxlGvZOc5/Zq/j+Q+WkJ6Z\nzS/PbsetZ3egWhXdcyAi34nm1UcLgJPy2f9oxPZB4LJoxSCB5A17eHhiMvM37OH0dg0YMbgb7RvX\njHVYIhKHdEdzGbYnPZPffbiU12eupX6NBH5/ZS8u6dkMM91zICL5U1Iog9ydd+Zt4sn3FrNzfwbX\nntaau87rSO2qlWMdmojEOSWFMmZFahqPTExhxqod9GxRh39d34fuLerEOiwRKSWUFMqI9EPZ/OnT\nYHqKapUr8uQPu3GV1jkQkaOkpFAGfLRoK49NWsjG3cH0FA9c0JmGmp5CRI6BkkIptmHXAYZPWsRH\ni7fSsUlN3h52Kqe01fQUInLslBRKoUNZOfxj+ir+8PFyDOOB/+vMjWe0oXJFTU8hIsdHSaGUmbFy\nB4+8k8KK1DTO79qERy/uSvO61WIdloiUEUoKpUTq3oM89f5iJs7bRMv61fjn9UkM6KzpKUSkeCkp\nxLnM7BxGfbWGlz5azqGsHG4b0J5f9G+v6SlEJCqUFOLYrFU7ePSdhSzduo/+nRox/OKutG5YI9Zh\niUgZpqQQhyK7iprXrcbIa3pzbpcmmp5CRKJOSSGOZGXnMGrGWl6cuoxDWTncPqA9P1dXkYiUICWF\nOKGuIhGJB0oKMZa69yBPT1nChG830rxuNf52TW/OU1eRiMSIkkKM5O0q0lVFIhIPlBRiYNaqHTw2\naSFLtuzjrI6NGH5JV9qoq0hE4oCSQglK3XeQp99XV5GIxC8lhRJwuKvopanLyFBXkYjEMSWFKJu9\neiePvpOiriIRKRWilhTMrCXwGnACkAOMdPff5ylTB/g3kBjG8lt3/1e0YipJ6ioSkdIomi2FLOBu\nd//GzGoBc81sqrsviijzS2CRu19sZo2ApWb2hrsfimJcUZWVncNr4VVFGVk53Hp2e355trqKRKR0\niFpScPfNwOZwe5+ZLQaaA5FJwYFaFnx9rgnsJEgmpVJkV9GZHRvxuLqKRKSUKZExBTNrDZwEzMrz\n0p+AScAmoBZwhbvn5PP7w4BhAImJidEM9Zhs2XOQp6cs5p1wrqKXr+7N+V3VVSQipU/Uk4KZ1QTG\nAXe4+948L58PzAMGAO2AqWZKVouGAAAMlUlEQVT2Rd5y7j4SGAmQlJTk0Y65qDKysvnn9DX88ZPl\nZOU4tw9ozy3921G9isbvRaR0iuqnl5lVJkgIb7j7+HyK3AA84+4OrDCz1UBnYHY04yoOny5NZcS7\ni1i9fT/ndmnCIxd2IbFB9ViHJSJyXKJ59ZEBrwCL3f13BRRbB5wDfGFmTYBOwKpoxVQc1u7YzxOT\nF/HR4lTaNqzBqzf0oX+nxrEOS0SkWESzpdAPuAZINrN54b4HCS4/xd1fBp4AXjWzZMCA+9x9exRj\nOmYHDmXxl09XMvKLVVSuYDzwf525oV8bqlSqEOvQRESKTTSvPppO8EF/pDKbgPOiFUNxcHcmL9jM\nU+8vZvOeg/zopObc/3+daVK7aqxDExEpdhoRPYIlW/YyfNJCZq7aSZemtfnjVSeR1Lp+rMMSEYka\nJYV87EnP5MWpy3h95lpqVa3Ekz/sxlV9E6lYQZeYikjZpqQQISfHGT1nPc/9dym7DxziJ6ckcve5\nnahXo0qsQxMRKRFKCqFv1+1i+KSFzN+wh6RW9Rh+SV+6Na8T67BEREpUuU8K2/Zl8NwHSxgzdwON\nayXw0hW9GNyrme5GFpFyqdwmhcxw4rqXpi7jYFY2N5/VltsGdKBmQrn9JxERKZ9J4asV23ls0kKW\np6ZxVsdGPHpxF9o1qhnrsEREYq5cJYWNu9P5zXuLeD95Cy3rV+Pv1yYx8MTG6ioSEQmVm6Qw+uv1\nPDopBYC7z+3ITWe2pWplrXEgIhKp3CSFVg2qM6BzYx66sAvN61aLdTgiInGp3CSFU9o24JS2DWId\nhohIXNNsbiIikktJQUREcikpiIhILiUFERHJpaQgIiK5lBRERCSXkoKIiORSUhARkVzm7rGO4aiY\n2TZgbazjOAYNge2xDqKElbc6l7f6gupcmrRy90aFFSp1SaG0MrM57p4U6zhKUnmrc3mrL6jOZZG6\nj0REJJeSgoiI5FJSKDkjYx1ADJS3Ope3+oLqXOZoTEFERHKppSAiIrmUFEREJJeSwnEys4pm9q2Z\nTQ6fm5n9xsyWmdliM7s9Yv8fzGyFmS0ws5MjjnGdmS0PH9fFqi5FlU+dzzGzb8xsnplNN7P24f4E\nM3s7rPMsM2sdcYwHwv1Lzez82NSkaMxsjZklh/WbE+6rb2ZTw/dsqpnVC/eXife5gDo/b2ZLwnpN\nMLO6EeXzfT/NbFC4b4WZ3R+LuhRFfvWNeO3XZuZm1jB8Xibe4wK5ux7H8QDuAv4DTA6f3wC8BlQI\nnzcOf14ATAEMOBWYFe6vD6wKf9YLt+vFul5HWedlwInh9i+AVyO2Xw63rwTeDre7APOBBKANsBKo\nGOt6HaG+a4CGefY9B9wfbt8PPFuW3ucC6nweUCncfjaizvm+n+FjJdAWqBKW6RLruhW1vuH+lsB/\nCW6YbViW3uOCHmopHAczawFcCPwjYvfPgRHungPg7qnh/sHAax6YCdQ1s6bA+cBUd9/p7ruAqcCg\nEqvEUSqgzg7UDrfrAJvC7cHAqHB7LHCOmVm4/y13z3D31cAKoG+0Yy9mkXUbBfwwYn+pf5/z4+4f\nuntW+HQm0CLcLuj97AuscPdV7n4IeCssW5q8CNxL8Dd+WJl9j0HdR8frJYI/mJyIfe2AK8xsjplN\nMbMO4f7mwPqIchvCfQXtj1f51flnwPtmtgG4Bngm3J9bt/DDZA/QgNJXZwc+NLO5ZjYs3NfE3TcD\nhD8bh/vLyvucX50j3UjwbRnKRp3/p75mdgmw0d3n5ylbFupboEqxDqC0MrOLgFR3n2tm/SNeSgAO\nunuSmV0K/BP4AUFTMy8/wv64c4Q63wlc4O6zzOwe4HcEiaLU1znUz903mVljYKqZLTlC2TJbZ3f/\nHMDMHgKygDfCsgXVLb8vnfFa5/ze44cIuszyKivvcb7UUjh2/YBLzGwNQbN4gJn9m+DbwbiwzASg\nR7i9gaB/8rAWBN0sBe2PR/nV+T2gp7vPCsu8DZwebufWzcwqEXQt7aR01Rl33xT+TCV4T/sCW8Mu\nA8Kfh7sJy8L7XFCdCQdPLwJ+6mFHOmWgzvnU9yyC8ZH54d97C+AbMzuBMlDfI4r1oEZZeAD9+W7Q\n9Rngxoj9X4fbF/L9wanZ4f76wGqCgal64Xb9WNepqHUmaG1uBzqG+4cC48LtX/L9gebR4XZXvj8w\nuYo4HWgGagC1Ira/Iugnfp7vDzQ/V1be5yPUeRCwCGiUp3y+72f4t7Eq3Hd4oLlrrOtX1PrmKbOG\n7waaS/17fKSHuo+K3zPAG2Z2J5BG0I0C8D7BVQsrgAMEVynh7jvN7Ang67DcCHffWbIhHzt3zzKz\nm4BxZpYD7CLobwZ4BXjdzFYQtBCuDH9noZmNJviAyQJ+6e7ZJR99kTQBJgTj41QC/uPuH5jZ18Bo\nMxsKrAMuC8uXhfe5oDqvIPjgnxq+NtPdbznS+2lmtxJcvVMR+Ke7Lyz56hQq3/oeoXxZeI8LpGku\nREQkl8YUREQkl5KCiIjkUlIQEZFcSgoiIpJLSUFERHIpKUhUmdmLZnZHxPP/mtk/Ip6/YGZ3FfM5\n04rzeOExe5nZBRHPh5vZr4vwe2Zmn5hZ7Yh9Pwpn3ewchThbm9lPivu4Ece/1cxuiNbxJfaUFCTa\nviK8w9nMKgANCW52Oux04MsYxHW0ehFcm360LgDmu/veiH1XAdMJ79soZq2BfJNCeFf58foncHsx\nHEfilJKCRNuXfDftRVcgBdhnZvXMLAE4EfjWzGqa2ccWrMuQbGaDAczsWTP7xeGDhd/Q7w637zGz\nr8M57R/P7+T5lQm/TS82s7+b2UIz+9DMqoWv9QnLzrBg/YAUM6sCjCCY6HCemV0RHr6LmU0zs1UW\nrpuRj58C70TEU5NgupChRCQFM+sfHmusBWsWvBHOKIuZXRDum27BPP6H17E4K4xnngXrW9QiuHny\nB+G+O83sejMbY2bvEkz4ZhH1Sj5cl/D8n5nZaAvWAnnGzH5qZrPDcu0A3P0AsMbMStustlJUsb6l\nWo+y/yCYIiARuBm4BXiC4Bt0P+DzsEwloHa43ZDgblEDTgI+izjWovBY5xEsoG4EX24mA2eGZdLC\nn/mWIfg2nQX0CsuNBq4Ot1OA08PtZ4CUcPt64E8RcQwnaAUlhPHuACrnU/e1hFMohM+vBl4Jt78C\nTg63+xPMItsijHUGcAZQlWDmzTZhuTf5bkqVdwkmcgOoGf4b9j/8ekTcGwinWwB+TDClc0WCO3nX\nAU3D39sdbicAG4HHw9/5FfBSxDEfAu6O9d+VHtF5qKUgJeFwa+F0gg+7GRHPvwrLGPCUmS0APiKY\ncriJu38LNDazZmbWE9jl7usIPvDPA74FvgE6Ax34viOVWe3u88LtuUBrC1YSq+Xuh2P6TyH1es+D\nNQS2E0yI1ySfMvXdfV/E86sIJhMk/HlVxGuz3X2DB2txzCNIXp2BVR6sUwBBUjjsS+B3YSulrn+3\n1kFeU/276RbOAN5092x33wp8BvQJX/va3Te7ewbB4jgfhvuTw1gOSwWaFXAuKeU095GUhMPjCt0J\nvomvB+4G9hL0UUPQzdII6O3umeHMlFXD18YCQ4AT+O4D1YCn3f1vRzhvvmUsWBY0I2JXNlCN/Kc+\nPpK8x8jv/1OWmVVw9xwzawAMALqZmRN8W3czu/cIxyswJnd/xoJZai8AZprZwAKK7o/YPlIdI8+f\nE/E8h+/XrSqQfoTjSCmmloKUhC8JplveGX5D3QnUBU4jaDVAMK12apgQzgZaRfz+WwT970MIEgQE\nk6zdGPbRY2bNLZgLP1JRyuTyYLWsfWZ2argrciB4H1DraCodWkqwHCVh/K+5eyt3b+3uLQlm0jzj\nCL+/BGhr361vfXg8AzNr5+7J7v4sMIegVVFYnJ8TjI1UNLNGBN1ps4+yTh0JkruUQUoKUhKSCfrd\nZ+bZtyfseoFgwZYkCxZN/ynBhyEQzKpK8EG30b9b7exDgu6dGWaWTJAsvvdhWJQy+RgKjDSzGQTf\nqveE+z8lGFiOHGguivcI+ush6CqakOf1cRRwtVBYh3SCta4/MLPpwNaImO4IB4znE3xznwIsIGid\nzLdgpt68JoRl5gOfAPe6+5ajqA8EY0EfHeXvSCmhWVJFIphZTXdPC7fvB5q6+6+O43hNCVoH5x5v\nTOHVSH8Glrv7i8d6vONhZicBd7n7NbE4v0SfWgoi33dh2BpIIVhG9cnjOVjYsvm7Rdy8dgxuMrN5\nwEKCbrYjjaNEW0PgkRieX6JMLQUREcmlloKIiORSUhARkVxKCiIikktJQUREcikpiIhIrv8HhYju\nCCXc9U4AAAAASUVORK5CYII=\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe3fbe25c0>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"sp.plot(bigbox.waveset, flux_unit=units.FNU)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 7, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAERCAYAAABRpiGMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3XecFfX1//HXm9770nsTaSIuYIsg\nMYpoBMREjbEnaKLf/BKjApbYFYzRaGzBBJXESBKKYMWGYkFlUdhdkF6XImXpS9lyfn/MrF4322D3\n7t29e56Px33s3JnPzD0fL+7Zmc/M+cjMcM45545VlVgH4JxzrmLzROKcc65EPJE455wrEU8kzjnn\nSsQTiXPOuRLxROKcc65EKk0ikTRZ0jZJqcVo20HSe5KSJX0gqW1ZxOiccxVRpUkkwAvAsGK2fQSY\nYmZ9gXuBh6IVlHPOVXSVJpGY2TwgPXKdpC6S3pK0UNJHknqEm3oC74XLc4ERZRiqc85VKJUmkRRg\nEvB/ZnYScDPwdLh+MTA6XB4F1JfUNAbxOedcuVct1gHEiqR6wKnAfyXlrq4Z/rwZeFLSVcA8YBOQ\nVdYxOudcRVBpEwnB2dhuM+uXd4OZbQYuhG8Tzmgz21PG8TnnXIVQaS9tmdleYK2knwAocEK43ExS\n7n+b8cDkGIXpnHPlXqVJJJJeBuYDx0lKk3QtcBlwraTFwBK+G1QfAiyXtAJoATwQg5Cdc65CkJeR\nd845VxKV5ozEOedcdFSKwfZmzZpZx44dYx2Gc85VKAsXLtxhZglFtasUiaRjx44kJSXFOgznnKtQ\nJK0vTju/tOWcc65EPJE455wrEU8kzjnnSsQTiXPOuRLxROKcc65EPJE455wrEU8kzjnnSsQTiXPO\nxaFNuw9yz6tLyMrOifpnVYoHEp1zrrLIyTFe+nw9E95cRo7BqBPb0Ldto6h+picS55yLE2u272fc\n9BS+WJfOD7o148FRfWjXpE7UP9cTiXPOVXBZ2Tk899FaHnt3BbWqVeGPF/XlopPaEjH7a1R5InHO\nuQpsyeY9jJ2eTOqmvQzr1ZJ7R/aief1aZRpD1AbbJU2WtE1SagHbe0iaL+mwpJvzbBsmabmkVZLG\nRazvJOlzSSsl/VtSjWjF75xz5dmhzGz+OGcZFzz5CVv3HOaZy/rz7OUnlXkSgejetfUCMKyQ7enA\nb4BHIldKqgo8BZwL9AQuldQz3DwReMzMugG7gGtLOWbnnCv3Fq5P57wnPuKpuasZ2a8N7950Buf2\naRWzeKKWSMxsHkGyKGj7NjNbAGTm2TQQWGVma8zsCDAVGKHgYt9QYFrY7kVgZOlH7pxz5dOBw1nc\nPXsJFz07n0OZObx4zUD+9NMTaFQnthdnyuMYSRtgY8T7NGAQ0BTYbWZZEevblHFszjkXE/NWbGf8\njBQ27znIlad05OZzjqNezfLxK7x8RPF9+d1mYIWsz/8g0hhgDED79u1LJzLnnCtjuzOOcP/rXzNt\nYRqdE+ry3+tOIbFjk1iH9T3lMZGkAe0i3rcFNgM7gEaSqoVnJbnr82Vmk4BJAImJiQUmHOecK6/e\nTNnCnbOWsCvjCDec2YX/G9qNWtWrxjqs/1EeE8kCoJukTsAm4BLgZ2ZmkuYCFxGMm1wJzIpdmM45\nFx3b9h3irllLeDN1K71aN+DFawbQq3XDWIdVoKglEkkvA0OAZpLSgLuA6gBm9qyklkAS0ADIkfRb\noKeZ7ZV0IzAHqApMNrMl4WHHAlMl3Q98Bfw9WvE751xZMzOmLUzj/te/5mBmNrcOO45f/qAz1auW\n77KIUUskZnZpEdu3Elyeym/bG8Ab+axfQ3BXl3POxZWN6RncNjOFj1buYEDHxkwY3ZcuCfViHVax\nlMdLW845V2nk5BhT5q/j4TnLEXDfiF5cNqgDVaqUTXmT0uCJxDnnYmTVtn2MnZ7CwvW7GNw9gQdG\n9aZt4+gXWSxtnkicc66MZWbnMGneGh5/dyV1albl0Z+ewKgT25RZkcXS5onEOefKUOqmPdw6LZml\nW/ZyXp9W3H1BLxLq14x1WCXiicQ558rAocxsHn9vJZPmraFJ3Ro8+/OTGNa7ZazDKhWeSJxzLsoW\nrEtn7LRk1uw4wMWJ7bht+PE0rFM91mGVGk8kzjkXJfsPZ/HwW8uYMn89bRvX5p/XDuL0bs1iHVap\n80TinHNRMHf5Nm6fkcKWvYe45rRO3HxOd+rUiM9fufHZK+eci5FdB45w32tLmfHVJro2r8e060/l\npA6NYx1WVHkicc65UmBmvJGylbtmp7I7I5PfDO3KDUO7UrNa+SuyWNo8kTjnXAlt23uIO15J5e2l\n39CnTUOmXDOInq0bxDqsMuOJxDnnjpGZ8d+kNO57fSlHsnIYf24Prj29E9XKeZHF0uaJxDnnjsGG\nnRmMn5nMJ6t2MrBTEyaO7kunZnVjHVZMeCJxzrmjkJ1jvPDpOh6Zs5yqVcT9I3vzs4HtK1SRxdLm\nicQ554pp5Tf7uHV6Ml9t2M2ZxyXwwKg+tG5UO9ZhxZwnEuecK8KRrBye/XA1T76/iro1q/Lni/sx\nol/rCltksbRFbURI0mRJ2ySlFrBdkp6QtEpSsqT+4fozJS2KeB2SNDLc9oKktRHb+kUrfuecA0hO\n280FT37Mo++s4JzeLXnnpsGMrMCVeqMhmmckLwBPAlMK2H4u0C18DQKeAQaZ2VygH4CkJsAq4O2I\n/W4xs2lRitk55wA4eCSbP7+7guc+WkNC/Zo8d0UiP+rZItZhlUvRnGp3nqSOhTQZAUwxMwM+k9RI\nUisz2xLR5iLgTTPLiFaczjmX12drdjJuejLrdmZw6cB2jB9+PA1qxU+RxdIWy5ud2wAbI96nhesi\nXQK8nGfdA+GlsMckFVjEX9IYSUmSkrZv3146ETvn4tq+Q5ncPjOFSyZ9Ro7Bv34xiIcu7OtJpAix\nHGzP7wKjfbtRagX0AeZEbB8PbAVqAJOAscC9+R3czCaFbUhMTLT82jjnXK73l33D7TNT+WbvIX5x\neid+f/Zx1K4R/+VNSkMsE0ka0C7ifVtgc8T7nwIzzSwzd0XEZa/Dkp4Hbo56lM65uJZ+4Aj3vrqE\nVxZtpnuLejx92amc2D6+iyyWtlgmktnAjZKmEgy278kzPnIpwRnIt3LHUBTcLjESyPeOMOecK4qZ\n8WryFu6evYR9hzL57Vnd+PWQrtSoVrnKm5SGqCUSSS8DQ4BmktKAu4DqAGb2LPAGMJzgrqwM4OqI\nfTsSnK18mOewL0lKILgstgi4PlrxO+fi19Y9h7jjlRTe/XobJ7RrxMOj+3Jcy/qxDqvCiuZdW5cW\nsd2AGwrYto7/HXjHzIaWSnDOuUrJzJi6YCMPvv41mTk53HHe8Vx9WieqVuLyJqXBn2x3zlUK63ce\nYNz0FOav2ckpnZsyYXQfOjStnEUWS5snEudcXMvOMZ7/ZC2PvL2c6lWq8NCFfbhkQDt/Mr0UeSJx\nzsWt5VuDIouLN+7mrOObc//IPrRsWCvWYcUdTyTOubhzJCuHp+au4ukPVlG/VnWeuPREfty3lZ+F\nRIknEudcXFm0cTe3TlvMim/2M7Jfa/7w4140qVsj1mHFNU8kzrm4cPBINn96ezmTP1lLiwa1mHxV\nIkN7eJHFsuCJxDlX4X26egfjpqewIT2Dywa1Z9y5Pajv9bHKjCcS51yFtfdQJg+98TUvf7GRjk3r\nMHXMyZzcuWmsw6p0PJE45yqkd5Z+wx2vpLB932GuO6Mzvz2ruxdZjBFPJM65CmXH/sPcPXsJryVv\noUfL+jx3RSJ92zaKdViVmicS51yFYGbMWrSZe15dwoHD2fz+R925bnAXL7JYDngicc6Ve5t3H+SO\nV1J5f9k2TmwfFFns1sKLLJYXnkicc+VWTo7xry82MOHNZWTnGH84vydXntrRiyyWM55InHPl0tod\nBxg3PZnP16ZzWtemPDSqL+2b1ol1WC4fnkicc+VKVnYOf/94LY++s4Ia1arw8Oi+/CSxrZc3Kcc8\nkTjnyo2lm/cydnoyKZv2cHbPFtw3sjctGniRxfIumjMkTgbOB7aZWe98tgt4nGCWxAzgKjP7MtyW\nDaSETTeY2QXh+k7AVKAJ8CVwuZkdiVYfnHNl43BWNk++v4pnPlhNozrVeepn/Rnep6WfhVQQ0bxv\n7gVgWCHbzwW6ha8xwDMR2w6aWb/wdUHE+onAY2bWDdgFXFu6ITvnytrC9bs474mP+cv7q7igX2ve\n+d1gzvNKvRVKNKfanRfOvV6QEcCUcMrdzyQ1ktTKzLbk1zg8gxkK/Cxc9SJwN99PQM65CiLjSBZ/\nnLOcFz5dR6sGtXj+6gGceVzzWIfljkEsx0jaABsj3qeF67YAtSQlAVnABDN7BWgK7DazrDzt8yVp\nDMGZDu3bty/96J1zx+zjlTsYNyOZtF0HueKUDtw6rAf1avqQbUUVy28uv/NWC3+2N7PNkjoD70tK\nAfYW0v5/N5hNAiYBJCYmFtjOOVd29mRk8sAbS/lPUhqdm9XlP9edwsBOTWIdliuhWCaSNKBdxPu2\nwGYAM8v9uUbSB8CJwHSgkaRq4VnJt+2dc+XfW6lbuXNWKukHjvCrIV34fz/sRq3qXmQxHsSySM1s\n4AoFTgb2mNkWSY0l1QSQ1Aw4DVgajqXMBS4K978SmBWLwJ1zxbd932FueOlLrv/nQhLq1WTWDacx\ndlgPTyJxpNhnJJIaA62Bg8A6M8spov3LwBCgmaQ04C6gOoCZPQu8QXDr7yqC23+vDnc9HvirpByC\nRDfBzJaG28YCUyXdD3wF/L248TvnypaZMePLTdz72lIOHsnmlnOOY8wZnale1YssxptCE4mkhsAN\nwKVADWA7UAtoIekz4Gkzm5vfvmZ2aWHHDs8wbshn/adAnwL2WQMMLOy4zrnY27T7ILfNSOHDFds5\nqUNjJo7uS9fm9WIdlouSos5IpgFTgB+Y2e7IDZJOAi6X1NnM/MzAOUdOjvHPz9cz8c1lGHD3j3ty\nxSkdqeJFFuNaoYnEzH5UyLaFwMJSj8g5VyGt3r6fcdOTWbBuFz/o1owHR/WhXRMvslgZFHVpq39h\n23NLmjjnKq/M7Bye+2gNf353JbWrV+WRn5zA6P5t/Mn0SqSoS1tJwBKCsRH4/rMfRvCkuXOukkrd\ntIex05NZsnkv5/ZuyT0jetG8vhdZrGyKSiS/B0YT3Kk1FZhpZvujHpVzrlw7lJnNX95fybMfrqFx\nnRo8c1l/zu3TKtZhuRgpaozkMeCxsOrupcB7ktYDD5rZorII0DlXviStS+fW6cms2X6Ai05qyx3n\nHU+jOjViHZaLoWI9R2JmayXNAmoDlwPdAU8kzlUiBw4HRRZfnL+O1g1rM+WagZzRPSHWYblyoKjB\n9s7AJQSVejcSXN56wMwOlUFszrly4sMV27ltRgqb9xzkylM6css5x1HXiyy6UFH/ElYByQSlSPYC\n7YFf596NYWaPRjU651xM7c44wn2vfc30L9PoklCX/153Cokdvcii+76iEsm9fFdhN+9jqV5R17k4\n9mbKFu6ctYRdGUe48cyu3Di0q9fHcvkqarD97oK2SRpQ6tE452Ju295D/GHWEt5aspVerRvw4jUD\n6NW6YazDcuXYUV3klNSTYMzkUmAPkBiNoJxzZc/MmLYwjfteW8qhrBzGDuvBL3/QiWpeZNEVochE\nIqkDQeK4lGDGwg5Aopmti25ozrmysjE9g9tmpvDRyh0M7NiEh0b3oUuCF1l0xVPUXVufAg0J7ta6\nyMxWSlrrScS5+JCdY0yZv44/zlmOgPtG9OKyQR28yKI7KkWdkWwnmImwBZAArMQH2Z2LC6u27WPs\n9BQWrt/F4O4JPHhhH9o0qh3rsFwFVOjFTzMbQTA3yJfAPZLWAo0lFTkniKTJkrZJSi1guyQ9IWmV\npOTcApGS+kmaL2lJuP7iiH1ekLRW0qLw1e9oOuucC4osPvn+SoY//jGrt+/n0Z+ewAtXD/Ak4o5Z\nkWMkZrYHmAxMltQcuBj4s6R2ZtaukF1fAJ4kmM8kP+cC3cLXIOCZ8GcGcEV4Ga01sFDSnIj5UG4x\ns2lFd805l1fqpj3cMi2Zr7fs5by+rbj7x71IqF8z1mG5Cu6o7toys23AX4C/SJpeRNt5kjoW0mQE\nMCWcKfEzSY0ktTKzFRHH2CxpG8Fltd0FHcg5V7hDmdn8+d2VPPfRGprUrcFfLz+Jc3q1jHVYLk6U\npMZBSZ8jaUNQdiVXWrhuS+6K8BJaDWB1RLsHJP0BeA8YZ2aH8zu4pDHAGID27duXMFTnKq7P1+xk\n3IwU1u44wMWJ7bht+PE0rFM91mG5OFKSG8RLeltHfvt/O5AvqRXwD+BqM8sJV48HehAksSbA2IIO\nbmaTzCzRzBITErywnKt89h3K5M5XUrl40mdk5eTw0i8GMfGivp5EXKkr6vbfgorqiJInkjQgcoyl\nLbA5/NwGwOvAHWb2WW4DM8s9Wzks6Xng5hLG4Fxcmrt8G7fPSGHL3kNcc1onbj6nO3VqeJFFFx1F\n/ctaSHCWkF/SOFLCz54N3ChpKsEg+x4z2yKpBjCTYPzkv5E7hGMoWxRUjRwJ5HtHmHOV1a4DR7jv\ntaXM+GoT3ZrXY/qvTqV/+8axDsvFuaISyRAzW38sB5b0MjAEaCYpDbgLqA5gZs8CbwDDCSoMZwBX\nh7v+FDgDaCrpqnDdVeFEWi9JSiBIbIuA648lNufijZnxesoW7pq1hD0HM/nN0K7cMLQrNat5kUUX\nfQpumipgo/SlmfUvw3iiIjEx0ZKSkmIdhnNR8c3eQ9zxSirvLP2GPm0a8vBFfTm+VYNYh+XigKSF\nZlZkTcWizki8ToJz5ZSZ8Z+kjdz/+tccycrhtuE9uOY0L7Loyl5RiaSNpCcK2mhmvynleJxzxbBh\nZwbjZiTz6eqdDOrUhImj+9KxWd1Yh+UqqaISyUGCAXfnXDmQnWO88Ok6HpmznKpVxAOjenPpgPZe\nZNHFVFGJZKeZvVgmkTjnCrXim33cOi2ZRRt3M7RHcx4Y1ZtWDb0+lou9ohJJSW/xdc6V0JGsHJ75\nYDVPzl1JvZrVePySflxwQmuCu+Cdi72ipto9OXyu4zKgF8EzJUuBfxVUmsQ5V3oWb9zN2OnJLNu6\njx+f0Jq7f9yTpvW8yKIrX4p6sv144FXgE4KxEhE8G3K7pAvMbGnUI3SuEjp4JJvH3l3B3z5aQ0L9\nmjx3RSI/6tki1mE5l6+iLm09CfzKzN6JXCnpLOAp4MxoBeZcZTV/9U7Gz0hm3c4MLh3YnvHDe9Cg\nltfHcuVXkbf/5k0iAGb2rqS/RCkm5yqlvYcymfDmMv71+QY6NK3Dv345iFO7NIt1WM4VqahEUkVS\nzbzjIZJqFWNf51wxvb/sG26bkcq2fYf45Q86cdOPjqN2DS9v4iqGopLBFGC6pBvNbB1AOFnVEwQl\n3p1zJbBz/2HufW0psxZt5rgW9Xn28pPo165RrMNy7qgUddfW/ZJuBOZJqkMw2L4feMTM/NKWc8fI\nzJi9eDP3vLqUfYcy+e1Z3fj1kK7UqOblTVzFU5w5258EnpRUP3y/L+pRORfHtuw5yB0zU3lv2TZO\naNeIh0f35biW9WMdlnPHrKjbf2/KZ923y2b2aBRici4u5eQYUxds5KE3viYzJ4c7zjueq0/rRFUv\nb+IquKLOSCL/TLoO+GsUY3Eubq3bcYBxM5L5bE06p3RuyoTRfejQ1IssuvhQ1BjJPbnLkkZGvnfO\nFS0rO4fnP1nHn95ZTvUqVZhwYR8uHtDOy5u4uHI0I3sFz4BVAEmTJW2TlO+UuAo8IWmVpGRJ/SO2\nXSlpZfi6MmL9SZJSwn2ekP8f6cqpZVv3MvqZT3ngja85vWsz3rlpMJcMbO9JxMWdaD8L8gLB0/FT\nCth+LtAtfA0CngEGSWpCMDVvIkECWyhptpntCtuMAT4jmK53GPBmFPvg3FE5nJXNU3NX8/TcVTSs\nXZ2/XHoi5/dt5QnExa2iBttTCH6RC+giKTl3E2Bm1rew/c1sXvjcSUFGAFMsmO/3M0mNJLUiqOf1\njpmlh3G8AwyT9AHQwMzmh+unACPxROLKia827GLs9GRWfLOfUSe24c7ze9Kkbo1Yh+VcVBV1RnJ+\nlD+/DbAx4n1auK6w9Wn5rP8fksYQnLnQvn370ovYuXxkHMniT2+vYPIna2nZoBaTr0pkaA8vsugq\nh6IG29dLGgl0BVLMbE4pf35+5/p2DOv/d6XZJGASQGJi4lGP7zhXXJ+u2sG4GSlsSM/g5ye3Z+yw\nHtT3IouuEinq0tbTBPOQfArcJ2mgmd1Xip+fBrSLeN8W2ByuH5Jn/Qfh+rb5tHeuzO05mMlDb3zN\n1AUb6di0DlPHnMzJnZvGOiznylxRl7bOAE4ws+ywRMpHQGkmktnAjZKmEgy27zGzLZLmAA9Kahy2\nOxsYb2bpkvZJOhn4HLgC8FItrsy9vWQrd7ySyo79h7lucGd+d1Z3alX3Iouucipyql0zywYws4yj\nvdVW0ssEZxbNJKUR3IlVPTzeswR3XQ0HVgEZwNXhtnRJ9wELwkPdmzvwDvyK4G6w2gSD7D7Q7srM\njv2HuXv2El5L3kKPlvX525WJ9G3rRRZd5abghqkCNkoZBL/kIbxzK3xfrLu2yovExERLSkqKdRiu\nAjMzXlm0iXteXUrG4Wz+b2hXrh/ShepVvciii1+SFppZYlHtijojOb6U4nGuwtq8+yC3z0xh7vLt\nnNg+KLLYrYUXWXQuV1GJZIMVdspC8HR6UW2cq4hycoyXvtjAxDeXkZ1j/OH8nlx5akcvsuhcHkUl\nkrmSpgOzzGxD7kpJNYDTgSuBuQRjFs7FjTXb9zNuegpfrEvn9K7NeOjCPrRrUifWYTlXLhWVSIYB\n1wAvS+oE7AZqAVWBt4HHzGxRdEN0ruxkZefwt4/X8tg7K6hZrQoPX9SXn5zU1subOFeIoh5IPAQ8\nDTwtqTrQDDhoZrvLIjjnytLSzXu5dfpiUjft5ZxeLbhvRG+aN6gV67CcK/eKXbTRzDKBLVGMxbmY\nOJyVzZPvr+KZD1bTqE51nr6sP+f2bulnIc4VU7Sr/zpXri1cHxRZXLVtPxf2b8Od5/WksRdZdO6o\neCJxldKBw1k88vZyXvh0Ha0b1uaFqwcw5LjmsQ7LuQqpWIlEUk8zW5pn3RAz+yAqUTkXRR+t3M74\nGSmk7TrIFad04NZhPahX0/+mcu5YFff/nv9I+gfwMMFdWw8TTDp1SrQCc6607cnI5P7Xl/LfhWl0\nblaX/1x3CgM7NYl1WM5VeMVNJIOAiQRVgOsDLwGnRSso50rbW6lbuXNWKukHjvDrIV34zQ+7eZFF\n50pJcRNJJnCQoFBiLWCtmeVELSrnSsm2fYe4e/YS3kjZSs9WDXj+qgH0btMw1mE5F1eKm0gWALOA\nAUBT4K+SLjKzi6IWmXMlYGZM/3IT9722lIOZ2dxyznGMOaOzF1l0LgqKm0iuNbPc8rlbgRGSLo9S\nTM6VSNquDG6bmcq8Fds5qUNjJo7uS9fm9WIdlnNxq7iJZJukvBOff1jawThXEjk5xj8+W8/Et5YB\ncM8Fvbj85A5U8SKLzkVVcRPJ63w3Z3otoBOwnGAa3gJJGgY8TlCb629mNiHP9g7AZCABSAd+bmZp\nks4EHoto2gO4xMxekfQCMBjYE267yut9udXb9zN2WjJJ63dxRvcEHhzVm7aNvciic2WhWInEzPpE\nvpfUH7iusH0kVQWeAn5EMNf6Akmz8zyP8ggwxcxelDQUeAi43MzmAv3C4zQhmEzr7Yj9bjGzacWJ\n3cW3zOwcJs1bw+PvraR29ao88pMTGN2/jZc3ca4MHdNTWGb2paQBRTQbCKwyszUA4bzsI4DIRNIT\n+F24PBd4JZ/jXAS8aWYZxxKri1+pm/YwdnoySzbvZXifltx9QS+a1/cii86VteI+2X5TxNsqQH9g\nexG7tQE2RrxPI3geJdJiYDTB5a9RQH1JTc1sZ0SbS4BH8+z3gKQ/AO8B48zscD4xjwHGALRvn3d4\nx1VkhzKzeeK9lfx13hoa16nBsz/vz7DerWIdlnOVVnHvhawf8apJMGYyooh98ru2kHcmxZuBwZK+\nIhj32ARkfXsAqRXQB5gTsc94gjGTAUATYGx+H25mk8ws0cwSExISigjVVRQL1qUz/PGPePqD1Vx4\nYhveu2mwJxHnYqy4YyT3HMOx04B2Ee/bApvzHHczcCGApHrAaDPbE9Hkp8DMsIR97j65pewPS3qe\nIBm5OLf/cBYPv7WMKfPX06ZRbaZcM5AzuvsfCM6VB4UmEkmv8r9nEd8yswsK2X0B0C2cWXETwSWq\nn+U5fjMgPXxKfjzBHVyRLg3XR+7Tysy2KBhNHQmkFtYHV/F9uGI7t81IYfOeg1x1akduOec46nqR\nRefKjaL+b3zkWA9sZlmSbiS4LFUVmGxmSyTdCySZ2WxgCPCQJAPmATfk7i+pI8EZTd7nVV6SlEBw\n6WwRcP2xxujKt90ZR7j3taXM+HITXRLqMu36UzipgxdZdK68kVmBJxxIam9mG8ownqhITEy0pKSk\nohu6cuONlC38YVYquzMyuX5wF24c2tWLLDpXxiQtNLPEotoVdUbyCsEdWkiabmajSyM45wqybe8h\n7pyVypwl39C7TQNevGYgvVp7kUXnyrOiEknknVedoxmIq9zMjP8uTOP+15ZyKCuHscN68MsfdKKa\nF1l0rtwrKpFYAcvOlZqN6RmMn5HCx6t2MLBjEyaM7kPnBC+y6FxFUVQiOUHSXoIzk9rhMuF7M7MG\nUY3OxbXsHGPK/HU8/NZyqgjuG9mbywa29yKLzlUwhSYSM/PRTRcVq7bt49ZpyXy5YTdDjkvggVF9\naNOodqzDcs4dA78Z35WpzOwcnv1gNX95fxV1alblsYtPYGQ/L7LoXEXmicSVmZS0PdwybTHLtu7j\nvL6tuOeCXjSrVzPWYTnnSsgTiYu6Q5nZPPbuCp6bt4Zm9Wry18tP4pxeLWMdlnOulHgicVH1+Zqd\njJuRwtodB7hkQDvGDz+ehrWrxzos51wp8kTiomLfoUwmvrWMf362gXZNavPSLwZxWtdmsQ7LORcF\nnkhcqZu7bBu3z0xhy95DXHv4/EGhAAAVE0lEQVR6J35/dnfq1PB/as7FK/+/25Wa9ANHuO+1pcz8\nahPdmtdj+q9OpX/7xrEOyzkXZZ5IXImZGa8lb+Hu2UvYczCT3/ywGzec2YWa1fwxJOcqA08krkS+\n2XuI22em8u7X39C3bUP++YtBHN/KCx44V5l4InHHxMz494KNPPDG1xzJyuG24T245jQvsuhcZeSJ\nxB21DTszGDcjmU9X72RQpyZMHN2Xjs3qxjos51yMRDWRSBoGPE4wQ+LfzGxCnu0dCKbXTQDSgZ+b\nWVq4LRtICZtuyJ3WN5y6dyrQBPgSuNzMjkSzHy6QnWM8/8laHnl7OdWqVOHBUX24ZEA7L7LoXCUX\ntUQiqSrwFPAjIA1YIGm2mS2NaPYIMMXMXpQ0FHgIuDzcdtDM+uVz6InAY2Y2VdKzwLXAM9Hqhwss\n37qPW6cns3jjbob2aM4Do3rTqqEXWXTOQTQvaA8EVpnZmvCMYSowIk+bnsB74fLcfLZ/j4LKfkOB\naeGqF4GRpRax+x9HsnL487srOP8vH7ExPYPHL+nH369M9CTinPtWNC9ttQE2RrxPAwblabMYGE1w\n+WsUUF9SUzPbCdSSlARkARPM7BWgKbDbzLIijtkmvw+XNAYYA9C+ffvS6VEls3jjbm6dlszyb/Zx\nwQmtuevHPWnqRRadc3lEM5Hkd+E87yyLNwNPSroKmAdsIkgcAO3NbLOkzsD7klKAvfyvfGduNLNJ\nwCSAxMREn93xKBw8ks2j7yzn7x+vpXn9WvztikTO6tki1mE558qpaCaSNKBdxPu2wObIBma2GbgQ\nQFI9YLSZ7YnYhpmtkfQBcCIwHWgkqVp4VvI/x3QlM3/1TsbNSGb9zgx+Nqg9487tQYNaXmTROVew\naI6RLAC6SeokqQZwCTA7soGkZpJyYxhPcAcXkhpLqpnbBjgNWGpmRjCWclG4z5XArCj2odLYeyiT\n8TNSuPS5zwD41y8H8eCoPp5EnHNFitoZiZllSboRmENw++9kM1si6V4gycxmA0OAhyQZwaWtG8Ld\njwf+KimHINlNiLjbaywwVdL9wFfA36PVh8ri3aXfcPsrKWzfd5gxZ3Tmd2d1p3YNL2/inCseBX/k\nx7fExERLSkqKdRjlzs79h7nn1aXMXryZ41rUZ+JFfenXrlGsw3LOlROSFppZYlHt/Mn2SsjMmL14\nM3fPXsL+w1n87qzu/GpIF2pU8/Imzrmj54mkktmy5yB3zEzlvWXb6NeuEQ9f1JfuLerHOiznXAXm\niaSSyMkxXl6wgYfeWEZWTg53nHc8V5/Wiape3sQ5V0KeSCqBtTsOMG56Mp+vTefULk2ZcGFf2jet\nE+uwnHNxwhNJHMvKzmHyJ2v509srqFG1ChMu7MPFA9oRVJpxzrnS4YkkTn29ZS9jpyeTnLaHs45v\nwf0je9OyYa1Yh+Wci0OeSOLM4axsnpq7mqfnrqJh7eo8+bMTOa9PKz8Lcc5FjSeSOPLlhl2MnZbM\nym37GXViG/5wfk8a160R67Ccc3HOE0kcyDiSxZ/eXsHkT9bSskEtnr9qAGf2aB7rsJxzlYQnkgru\nk1U7GDcjmY3pB/n5ye0ZO6wH9b0+lnOuDHkiqaD2HMzkwde/5t9JG+nUrC7/HnMygzo3jXVYzrlK\nyBNJBfT2kq3c8UoqO/Yf5rrBQZHFWtW9yKJzLjY8kVQg2/cd5u5Xl/B68hZ6tKzP365MpG9bL7Lo\nnIstTyQVgJnxyqJN3PPqUjIOZ3Pz2d25bnAXqlf1IovOudjzRFLObdp9kNtnpvDB8u30bx8UWeza\n3IssOufKj6j+SStpmKTlklZJGpfP9g6S3pOULOkDSW3D9f0kzZe0JNx2ccQ+L0haK2lR+OoXzT7E\nSk6O8Y/56zj70Q/5fE06d/24J/+9/lRPIs65cidqZySSqgJPAT8imL99gaTZETMdAjwCTDGzFyUN\nBR4CLgcygCvMbKWk1sBCSXPMbHe43y1mNi1ascfamu37GTc9hS/WpXN612Y8dGEf2jXxIovOufIp\nmpe2BgKrzGwNgKSpwAggMpH0BH4XLs8FXgEwsxW5Dcxss6RtQAKwmziWlZ3Dcx+t5bF3V1CrWhUe\nvqgvPzmprZc3cc6Va9G8tNUG2BjxPi1cF2kxMDpcHgXUl/S9hyEkDQRqAKsjVj8QXvJ6TFLN/D5c\n0hhJSZKStm/fXpJ+lImlm/cy8ulPmPjWMs48LoF3bxrMTxO9Uq9zrvyLZiLJ7zdg3gnibwYGS/oK\nGAxsArK+PYDUCvgHcLWZ5YSrxwM9gAFAE2Bsfh9uZpPMLNHMEhMSEkrUkWg6lJnNI3OWc8GTH7N1\nz2Geuaw/f708keYNvFKvc65iiOalrTSgXcT7tsDmyAZmthm4EEBSPWC0me0J3zcAXgfuMLPPIvbZ\nEi4elvQ8QTKqkBauT+fWacms3n6A0f3bcuf5x9OojhdZdM5VLNFMJAuAbpI6EZxpXAL8LLKBpGZA\neni2MR6YHK6vAcwkGIj/b559WpnZFgXXfEYCqVHsQ1QcOJzFH+cs58X562jdsDYvXjOQwd3L71mT\nc84VJmqJxMyyJN0IzAGqApPNbImke4EkM5sNDAEekmTAPOCGcPefAmcATSVdFa67yswWAS9JSiC4\ndLYIuD5afYiGeSu2M35GCpt2H+TKUzpwy7Ae1Kvpj/M45youmeUdtog/iYmJlpSUFNMY9mRkct/r\nS5m2MI3OCXWZOLovAzo2iWlMzjlXGEkLzSyxqHb+p3AZeCt1C3fOWkL6gSP8ekgXfvPDbl5k0TkX\nNzyRRNG2fYe4a9YS3kzdSs9WDXj+qgH0btMw1mE551yp8kQSBWbGtIVp3P/61xzMzOaWc45jzBmd\nvciicy4ueSIpZRvTM7htZgofrdxBYofGTBjdl67N68U6LOecixpPJKUkJ8eYMn8dD89ZDsA9F/Ti\n8pM7UKWKP5nunItvnkhKwapt+xk3PZmk9bs4o3sCD47qTdvGXmTROVc5eCIpgczsHCbNW8Pj766k\ndo2q/OknJ3Bh/zZeH8s5V6l4IjlGqZv2cOu0ZJZu2cvwPi2554LeJNTPt36kc87FNU8kR+lQZjaP\nv7eSSfPW0KRuDZ79eX+G9W4V67Cccy5mPJEchQXr0hk7LZk1Ow7wk5Pacsd5PWlYp3qsw3LOuZjy\nRFIM+w9n8fBby5gyfz1tG9fmH9cO5AfdvMiic86BJ5IifbB8G7fPTGXznoNcfVpHbj77OOp6kUXn\nnPuW/0YsxPgZKbz8xQa6Nq/HtOtP5aQOjWMdknPOlTueSArRsWkd/m9oV24c2pWa1bzIonPO5ccT\nSSGuG9wl1iE451y551UEnXPOlUhUE4mkYZKWS1olaVw+2ztIek9SsqQPJLWN2HalpJXh68qI9SdJ\nSgmP+YT8MXLnnIupqCUSSVWBp4BzgZ7ApZJ65mn2CMG87H2Be4GHwn2bAHcBg4CBwF2Scke6nwHG\nAN3C17Bo9cE551zRonlGMhBYZWZrzOwIMBUYkadNT+C9cHluxPZzgHfMLN3MdgHvAMMktQIamNl8\nC+YIngKMjGIfnHPOFSGaiaQNsDHifVq4LtJiYHS4PAqoL6lpIfu2CZcLOyYAksZISpKUtH379mPu\nhHPOucJFM5HkN3Zhed7fDAyW9BUwGNgEZBWyb3GOGaw0m2RmiWaWmJDgT6E751y0RPP23zSgXcT7\ntsDmyAZmthm4EEBSPWC0me2RlAYMybPvB+Ex2+ZZ/71jOuecK1vRPCNZAHST1ElSDeASYHZkA0nN\nJOXGMB6YHC7PAc6W1DgcZD8bmGNmW4B9kk4O79a6ApgVxT4455wrQtTOSMwsS9KNBEmhKjDZzJZI\nuhdIMrPZBGcdD0kyYB5wQ7hvuqT7CJIRwL1mlh4u/wp4AagNvBm+CrVw4cIdktaXWufKRjNgR6yD\nKGPe58rB+1xxdChOIwU3P7nyRlKSmSXGOo6y5H2uHLzP8cefbHfOOVcinkicc86ViCeS8mtSrAOI\nAe9z5eB9jjM+RuKcc65E/IzEOedciXgicc45VyKeSMqYpKqSvpL0Wvhekh6QtELS15J+E7H+ibBc\nfrKk/hHHyLfEfnmVT59/KOlLSYskfSypa7i+pqR/h33+XFLHiGOMD9cvl3RObHpSPJLWhVMdLJKU\nFK5rIumd8Dt7J7eadbx8zwX0+Y+SloX9mimpUUT7fL/PoqaeKE/y63PEtpslmaRm4fu4+J4LZGb+\nKsMXcBPwL+C18P3VBFWMq4Tvm4c/hxM8bCngZODzcH0TYE34s3G43DjW/TrKPq8Ajg+Xfw28ELH8\nbLh8CfDvcLknQYHPmkAnYDVQNdb9KqS/64BmedY9DIwLl8cBE+Ppey6gz2cD1cLliRF9zvf7DF+r\ngc5AjbBNz1j37Wj6HK5vR/Ag9vrc7fHyPRf08jOSMqRg4q7zgL9FrP4VwZP7OQBmti1cP4JgrhYz\ns8+ARmEZ/XxL7JdZJ45SAX02oEG43JDv6qWNAF4Ml6cBPwxL4YwApprZYTNbC6wimKagIons24t8\nN/1BXHzP+TGzt80sK3z7Gd/VySvo+yzO1BMVwWPArXy/oGzcfs/gl7bK2p8J/oHlRKzrAlysoOT9\nm5K6hesLK6VfVHn+8iS/Pv8CeCMsznk5MCFc/23fwl9Ae4DCphUorwx4W9JCSWPCdS0sqBVH+LN5\nuD5evuf8+hzpGr4rZxS3fZZ0AbDJzBbnaRsvfc5XNKv/ugiSzge2mdlCSUMiNtUEDplZoqQLCQpX\n/oBSKKUfa4X0+XfAcDP7XNItwKMEyaXC9zl0mpltltQceEfSskLaxm2fzWwegKTbCaaHeClsW1Df\n8vvDtkL1Gbid4JJeXvHyPefLz0jKzmnABZLWEZyyD5X0T4K/QKaHbWYCfcPlgsrwF1mevxzJr8+v\nAyeY2edhm38Dp4bL3/ZNUjWCy17pVKw+Y8H0CLmXKWcSXLL5JryUQfgz9xJmPHzPBfWZcPD4fOAy\nCwcFiN8+DyYY81kc/ptvC3wpqSVx0ucCxXqQpjK+CKoe5w48TwCuiVi/IFw+j+8Pzn0Rrm8CrCUY\nmGscLjeJdZ+K22eCs+AdQPdw/bXA9HD5Br4/2P6fcLkX3x+cXUM5HWwH6gL1I5Y/Jbjm/Ue+P9j+\ncLx8z4X0eRiwFEjI0z7f7zP8t7EmXJc72N4r1v07mj7nabOO7wbbK/z3XNjLL23F3gTgJUm/A/YT\nXOIBeIPgTo9VQAbB3V1Y4SX2yz0Lphf4JTBdUg6wi+D6OcDfgX9IWkVwJnJJuM8SSf8h+KWUBdxg\nZtllH32xtABmBvcIUA34l5m9JWkB8B9J1wIbgJ+E7ePhey6oz6sIksU74bbPzOz6wr5P5TP1RNl3\np1jy7XMh7ePhey6Ql0hxzjlXIj5G4pxzrkQ8kTjnnCsRTyTOOedKxBOJc865EvFE4pxzrkQ8kbhy\nR9Jjkn4b8X6OpL9FvP+TpJtK+TP3l+bxwmP2kzQ84v3dkm4uxn6S9L6kBhHrRoXVZHtEIc6Okn5W\n2seNOP6Nkq6O1vFd7HkiceXRp4RPu0uqAjQjeIgt16nAJzGI62j1I3h24GgNBxab2d6IdZcCHxM+\nW1PKOgL5JpKwwkBJTQZ+UwrHceWUJxJXHn3Cd2VTegGpwD5JjSXVBI4HvpJUT9J7CuY2SZE0AkDS\nREm/zj1YeCbw+3D5FkkLwjkh7snvw/NrE/7V/rWk5yQtkfS2pNrhtgFh2/kK5uBIlVQDuJegIOci\nSReHh+8p6QNJaxTOPZOPy4BZEfHUIyg3cy0RiUTSkPBY0xTM+/FSWC0ZScPDdR8rmAcjdy6YwWE8\nixTMEVOf4KHYH4TrfifpKkn/lfQqQVFCRfQrJbcv4ed/KOk/CubTmSDpMklfhO26AJhZBrBOUkWr\n2OyKK9aP1vvLX/m9CMpLtAeuA64H7iP4S/00YF7YphrQIFxuRvDUsIATgQ8jjrU0PNbZwKSwTRWC\nki1nhG32hz/zbUPwV3sW0C9s9x/g5+FyKnBquDwBSA2XrwKejIjjboKzrZphvDuB6vn0fT1h+Y3w\n/c+Bv4fLnwL9w+UhBBWS24axzgdOB2oRVJTtFLZ7me9K8rxKUGwQoF7433BI7vaIuNMIS3UAownK\nm1cleKJ7A9Aq3G93uFwT2ATcE+7z/4A/RxzzduD3sf535a/ovPyMxJVXuWclpxL8gpwf8f7TsI2A\nByUlA+8SlN9uYWZfAc0ltZZ0ArDLzDYQJImzga+AL4EeQDe+r7A2a81sUbi8EOioYNa/+maWG9O/\niujX6xbMw7GDoHBji3zaNDGzfRHvLyUoekn489KIbV+YWZoF89ksIkh4PYA1Fsz1AUEiyfUJ8Gh4\nNtTIvpsvJK937LtSHacDL5tZtpl9A3wIDAi3LTCzLWZ2mGBSqrfD9SlhLLm2Aa0L+CxXwXmtLVde\n5Y6T9CH4i38j8HtgL8E1dwguASUAJ5lZZlhxtVa4bRpwEdCS734JC3jIzP5ayOfm20bBtL+HI1Zl\nA7XJvwx4YfIeI7//B7MkVTGzHElNgaFAb0lGcFZgkm4t5HgFxmRmExRUYB4OfCbprAKaHohYLqyP\nkZ+fE/E+h+/3rRZwsJDjuArMz0hcefUJQfnx9PAv4XSgEXAKwdkJBGXmt4VJ5EygQ8T+UwnGEy4i\nSCoQFAO8JhxzQFIbBXNJRCpOm29ZMKvdPkknh6siB8P3AfWPptOh5QTTzRLGP8XMOphZRzNrR1Ah\n9vRC9l8GdNZ3c97njs8gqYuZpZjZRCCJ4OylqDjnEYz1VJWUQHCp74uj7FN3gj8IXBzyROLKqxSC\ncYTP8qzbE14WgmCipERJSQRnJ99OIGVB1dj6BLPV5c5M+DbBpaf5klIIEsz3foEWp00+rgUmSZpP\n8Nf7nnD9XILB9cjB9uJ4nWD8AYLLWDPzbJ9OAXdZhX04CPwaeEvSx8A3ETH9Nhw0X0xwhvAmkExw\nFrRYQRXqvGaGbRYD7wO3mtnWo+gPBGNb7x7lPq6C8Oq/zpWQpHpmtj9cHge0MrP/V4LjtSI4C/lR\nSWMK7+J6ClhpZo8d6/FKQtKJwE1mdnksPt9Fn5+ROFdy54VnHakE0yTfX5KDhWdQzynigcRj8EtJ\ni4AlBJcACxsXirZmwJ0x/HwXZX5G4pxzrkT8jMQ551yJeCJxzjlXIp5InHPOlYgnEueccyXiicQ5\n51yJ/H/wZuxhv6aTlAAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe7068cdd8>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"sp.plot(bigbox.waveset) # PHOTLAM" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 8, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZIAAAEKCAYAAAA4t9PUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3XecVNX5x/HPA0uTKh0p0qsUcQV7\nL1giFoxYEo1ETaKJXTGosUQFNbaAPzXRRI2NgChYwIYFBWWld1bq0mHpy8KW5/fHvcRx3d3ZZXd2\nZne/79drXtx75tw7z3GQZ849955j7o6IiMiBqhLvAEREpHxTIhERkRJRIhERkRJRIhERkRJRIhER\nkRJRIhERkRJRIhERkRJRIhERkRJRIhERkRJJincAZaFx48betm3beIchIlKufP/995vdvUm0epUi\nkbRt25aUlJR4hyEiUq6Y2cqi1NOlLRERKRElEhERKRElEhERKRElEhERKRElEhERKRElEhERKREl\nEhERKRElEhGRCih14y4enbiIslhOXYlERKSC+XDuOgaOnMJb01ezdntmzD+vUjzZLiJSGWTn5PLY\npMU8/+Uy+rRuwP9d0ZcW9WvF/HOVSEREKoDNu/byx9dnMnXZFn511KHcfW43aiRVLZPPjumlLTMb\nYGaLzSzVzIbm834NM3srfP9bM2sbljcys8lmtsvMRuY55iEzW21mu2IZu4hIefH9ynTOfWYKM1Zt\n5W8X9+bB8w8rsyQCMUwkZlYVGAWcBXQHLjWz7nmqDQG2untH4ElgRFieCdwD3JbPqScA/WIStIhI\nOeLuvDRlOZc8P40a1arw9h+O4aIjWpV5HLG8tNUPSHX3ZQBm9iYwEFgQUWcgcF+4PQYYaWbm7ruB\nKWbWMe9J3X1aeL4Yhi4ikth27c3mzrFzeH/OOk7v3ozHL+5N/VrV4hJLLBNJS2B1xH4a0L+gOu6e\nbWbbgUbA5hjGJSJSri3ZsJPf/ed7VmzezdCzunLdCe3j+uM6lokkv1blvaG5KHUO7MPNrgWuBWjT\npk1pnFJEJO7enbWGoWPnUrtGEq/99iiO7tAo3iHFNJGkAa0j9lsBawuok2ZmSUB9IL00PtzdXwBe\nAEhOTo79EzkiIjG0NzuHh95fyCtTV9KvbUNGXnY4TevVjHdYQGwTyXSgk5m1A9YAg4HL8tQZD1wJ\nTAUGAZ95WTyGKSJSjqzZtoc/vDaD2au3ce0J7bn9zC5Uq5o4z5PHLJGEYx43AJOAqsBL7j7fzB4A\nUtx9PPAi8KqZpRL0RAbvP97MVgD1gOpmdj5whrsvMLNHCRLSQWaWBvzT3e+LVTtEROLpyyWbuPHN\nmWTlOM9d0ZcBh7WId0g/Y5WhA5CcnOxas11EypPcXOfvn6Xy1KdL6NKsLs9e3pf2TeqUaQxm9r27\nJ0erpyfbRUQSzNbd+7jprVl8sWQTF/ZtyUPn96RW9bJ7wLC4lEhERBLI7NXb+MNrM9i0cy8PX9CT\nS/u1Tvjn5pRIREQSgLvzn29X8eCEBTSpW4Mxvz+aXq0axDusIlEiERGJs4x92QwbN49xM9dwUpcm\nPPnLPhxcu3q8wyoyJRIRkThK3biL61+bwZKNO7nl9M7ccHJHqlRJ7EtZeSmRiIjEyTsz1/DncXOp\nWa0qr1zdj+M7NYl3SAdEiUREpIxlZuVw/4QFvPHdKvq1bcgzlx5O8/qJ8ZT6gVAiEREpQ8s37+YP\nr81g4bod/P6kDtx6emeSEugp9QOhRCIiUkbem7OWoWPnklTV+NdVR3Jy16bxDqlUKJGIiMRY5ISL\nfds04O+X9aVlg9ivpV5WlEhERGJo1ZYMrn99BnPXbOea49txx4CuCTXhYmlQIhERiZGJ89Zz+5jZ\nGPDCr47gjB7N4x1STCiRiIiUsn3ZuQz/cBEvfb2c3q3qM/KyvrRueFC8w4oZJRIRkVKUtjWDG16f\nyazV27jqmLbcdXZXaiQl7oSLpUGJRESklHy6cAO3jJ5NTq4z6rK+nNMr8dYOiQUlEhGREsrKyeXx\nSYt5/stldG9Rj2cv70vbxrXjHVaZUSIRESmBddv38MfXZ5KyciuX9W/Dved2p2a1in0pKy8lEhGR\nA/T54o3cMno2mVk5PD24DwP7tIx3SHGhRCIiUkyRl7K6Nq/LyMv60rFp2S6Dm0iUSEREimF1egZ/\nfCO4K+vy/m24pxJeyspLiUREpIg+nLuOO8bOAadS3ZUVjRKJiEgUmVnBXFmvTltJ71b1+fulfWnT\nqOI+YFhcSiQiIoVYtmkX178+k4XrdnDN8e24/cyuVE+qWHNllZQSiYhIAd6ekcbd78yjRlIVXroq\nmVO6Not3SAlJiUREJI/de7O59935jJ2RRr+2DXn60j60qF9xpn0vbUokIiIRFq7bwQ2vz2DZ5t38\n6ZSO/OnUTuV+BcNYUyIREQHcnde+XcUD7y2gfq1qvDakP8d0bBzvsMoFJRIRqfR2ZGZx19i5vD93\nHcd3asyTl/ShcZ0a8Q6r3FAiEZFKbfbqbdzwxgzWbsvkzgFdue6E9lSpYvEOq1xRIhGRSsndeXHK\nckZMXETTujUZfd1RHHFow3iHVS4pkYhIpbNl115uHzOHzxZt5IzuzXh0UC8aHFQ93mGVWzG9FcHM\nBpjZYjNLNbOh+bxfw8zeCt//1szahuWNzGyyme0ys5F5jjnCzOaGxzxjZuqDikiRfbV0EwOe/oop\nqZu57xfdef5XRyiJlFDMEomZVQVGAWcB3YFLzax7nmpDgK3u3hF4EhgRlmcC9wC35XPq/wOuBTqF\nrwGlH72IVDT7snN55IOF/OrF76hfqxrvXn8sVx3bDv0WLblY9kj6Aanuvszd9wFvAgPz1BkIvBxu\njwFONTNz993uPoUgofyPmbUA6rn7VHd34BXg/Bi2QUQqgBWbdzPouW94/stlXNqvDRNuOI5uLerF\nO6wKI5ZjJC2B1RH7aUD/guq4e7aZbQcaAZsLOWdannPmu5KMmV1L0HOhTZs2xY1dRCqIt2ekcc87\n80iqWoXnrujLgMM0Y29pi2Uiya+/6AdQ54Dqu/sLwAsAycnJhZ1TRCqgnZlZ3PPOPN6ZtZZ+7Rry\n1CV9OKSBpjmJhVgmkjSgdcR+K2BtAXXSzCwJqA+kRzlnqyjnFJFKbuaqrdz45izStmZwy+mduf7k\njlTVsyExE8sxkulAJzNrZ2bVgcHA+Dx1xgNXhtuDgM/CsY98ufs6YKeZHRXerfVr4N3SD11EyqOc\nXGfU5FQufm4qObnO6OuO5k+ndlISibGY9UjCMY8bgElAVeAld59vZg8AKe4+HngReNXMUgl6IoP3\nH29mK4B6QHUzOx84w90XAL8H/g3UAj4MXyJSya3fnsnNb81i6rItnNOrBQ9f0JP6tarFO6xKwQrp\nAFQYycnJnpKSEu8wRCRGPlmwgdvHzCYzK5f7z+vBxcmtdFtvKTCz7909OVo9PdkuIuVWZlYOD3+w\nkFemrqR7i3r8/bLD6dCkTrzDqnSUSESkXFqyYSd/fH0mizfsZMhx7bhjQBdqJFWNd1iVkhKJiJQr\n+9cNefC9BdSpkcS/fnMkJ3dpGu+wKjUlEhEpN7bs2sudY+fyycINHN+pMX/7ZW+a1q0Z77AqPSUS\nESkXPl+8kdvHzGF7RhZ3n9ONq49tp3VDEoQSiYgktMysHIZ/uIh/f7OCzs3q8MrV/TRPVoJRIhGR\nhLVw3Q5ufHMmSzbs4qpj2jL0rK7UrKYB9USjRCIiCSc313np6+U8OnEx9Q+qxstX9+PEzk3iHZYU\nQIlERBLKhh2Z3Dp6NlNSN3Nat2aMuKgnjerUiHdYUgglEhFJGBPnrWPo23PZm5XLwxf05NJ+rfWE\nejmgRCIicbd7bzb3T5jP6JQ0erasz9OD+9BeT6iXG8VKJGZWG8h095wYxSMilczMVVu56a1ZrErP\n4A8ndeCm0zpTPSmWE5NLaSs0kZhZFYIZeS8HjgT2AjXMbBPwAfCCuy+NeZQiUuFk5+Ty7Oc/8PSn\nS2leryZvXnMU/ds3indYcgCi9UgmA58AdwHz3D0XwMwaAicDw81snLv/J7ZhikhFsjo9g5vfmkXK\nyq0M7HMIDww8TFO+l2PREslp7p6Vt9Dd04GxwFgz07cvIkXi7oybuYZ7352PAU9d0ofzD28Z77Ck\nhKIlkrp57phwYFvkKob5JRoRkby2Zezj7nfm8d6cdRzZ9mCe+GUfWjc8KN5hSSmIlki+J0gekdmk\njpnNBn7r7itiFZiIVBxfLNnEHWNms2XXPm4/swu/O7GDlr+tQApNJO7eLr9yM7sQeA4YEIugRKRi\nyNiXzSMfLOLVaSvp1LQOL155JIe1rB/vsKSUHdBzJO7+tpndXdrBiEjFMXPVVm4ZPZvlm3cz5Lh2\n3H5mF82TVUEdUCIxszqAbvQWkZ/Jysnl75+lMmpyKs3q1uD13/bnmI6N4x2WxFC050huyaf4YOA8\nYGRMIhKRcit14y5uGT2LOWnbufDwlvzlvB66rbcSiHrXVp59B9YDV7j73NiEJCLlTW6u8/LUFQz/\ncBEHVa/K/13el7N6toh3WFJGog22319WgYhI+bRu+x5u/+8cpqRu5uQuTRhxUS+a1tPyt5VJtEtb\n/yLoheTH3X1I6YckIuXFu7PWcM8788jKcR664DAu69dGs/VWQtEubb2XT1kb4CZAt1+IVFKRDxf2\nbdOAJ37Zh7aNa8c7LImTaJe2xu7fNrP2wJ+BE4DhwIuxDU1EElHehwuvO6E9SVV1E2dlFvX2XzPr\nBgwDDgceA37n7tmxDkxEEoseLpSCRBsj+S+QDDwO3AzkAPX2XwMNJ28UkQpuxqqt3DZ6Nsv0cKHk\nI1qP5EiCwfbbgFvDsv0jaQ60j1FcIpIA9mbn8NQnS3n+ix9oXq+mHi6UfEVLJCe6+8oyiUREEsq8\nNdu5dfRsFm/YySXJrbn73G7UramHC+Xnoo2QjSvJyc1sgJktNrNUMxuaz/s1zOyt8P1vzaxtxHt3\nheWLzezMiPIbzWyemc03s5tKEp+I/FxWTi5PfbKE80d9zdaMfbx0VTIjBvVSEpECReuRHPAN4WZW\nFRgFnA6kAdPNbLy7L4ioNgTY6u4dzWwwMAK4xMy6Eyzx2wM4BPjEzDoD3YBrgH7APmCimb2v5X5F\nSsfi9Tu59b+zmLdmB+f3OYT7zutBg4OqxzssSXDREklLM3umoDfd/U+FHNsPSHX3ZQBm9iYwEIhM\nJAOB+8LtMcBIC0byBwJvuvteYLmZpYbnawVMc/eM8JxfABcAj0Zph4gUIifXeeHLZTz58RLq1kzi\nuSv6MuAwTXEiRRMtkewhWNzqQLQEVkfspwH9C6rj7tlmth1oFJZPy3NsS2Ae8JCZNQpjOxtIye/D\nzexa4FqANm3aHGATRCq+ZZt2cet/ZzNz1TbOOqw5fz3/MBrVqRHvsKQciZZItrj7ywd47vwui+Wd\nbqWgOvmWu/tCMxsBfAzsAmYD+T7T4u4vAC8AJCcnFzTNi0illZvr/PubFYyYuIia1ary9OA+nNf7\nEE1xIsUWLZHsK8G504DWEfutgLUF1EkzsySgPpBe2LHu/iLhU/Vm9nBYV0SKYXV6Brf9dzbfLk/n\nlK5NeeTCnjTTRItygKJNkXKUmVUHLicY+HaCMY7Xw/GLwkwHOplZO2ANweD5ZXnqjAeuBKYCg4DP\n3N3NbDzwupk9QTDY3gn4DsDMmrr7RjNrA1wIHF3k1opUcu7O69+t4qH3F1LFjEcH9eLiI1qpFyIl\nEu3J9m7ABOBrgrESA04ChpnZeXnuwPqJcMzjBmASwQSPL7n7fDN7AEhx9/EEPYtXw8H0dIJkQ1hv\nNEHSygaud/ec8NRjwzGSrLB86wG2XaRSWbttD3eOncNXSzdzXMfGjBjUi5YNasU7LKkAzL3g4QMz\n+xQY7u4f5yk/DRjm7ifHOL5SkZyc7Ckp+Y7Ji1R47s7YGWu4f8J8snOcP5/TjSv6a7p3ic7Mvnf3\n5Gj1ot7+mzeJALj7J2b29wOOTkTKxIYdmQwbN5dPFm6kX9uGPHZxLw5tpOnepXRFSyRVzKxG3vEQ\nM6tZhGNFJE7290IemDCfvdm53H1ON35zbDuqVlEvREpftGTwCsGYxA3uvgIgnMbkGeDVmEYmIgdk\n3fY9/PntuUxevIkj2x7MiIt60b5JnXiHJRVYtLu2/hoOmH9pZgcRDLbvAh53d13aEkkg7s7olNX8\n9b2FZOXm8pdfdOfKo9tSRb0QibGol6fcfSTB1CV1w/2dMY9KRIplzbY9DA3vyOrfriGPDtJYiJSd\naLf/3pJP2f+23f2JGMQkIkXk7rzx3Woe/mAhue48OLAHl/c/VL0QKVPReiR1I7avA56PYSwiUgyr\n0zMY+vYcvk7dwjEdGjHiol60bnhQvMOSSijaGMn9+7fN7PzIfRGJj9xc57XvVvHIBwsx4KELDuOy\nfnouROKnOLfwauJDkThbtSWDO8bOZtqydI7v1JhHLuxJq4PVC5H40rMgIuVAbq7zytQVjJi4mKQq\nxoiLevLL5NbqhUhCiDbYPpcfp3XvYGZz9r9FMK17rxjHJ1LpLd+8mzvHzOG7Femc1KUJj1zYkxb1\nNUeWJI5oPZJzyyQKEfmZnHC9kMcmLaJa1So8NqgXgzRTrySgaIPtK83sfKAjMNfdJ5VNWCKV29IN\nO7lj7BxmrtrGqV2b8rDWC5EEFu3S1rME65B8AzxoZv3c/cEyiUykEtqXncv/ff4DIycvpU6NJJ66\npA8D+2jVQkls0S5tnQD0dveccIqUrwAlEpEYmLV6G3eOmcPiDTs5r/ch/OUX3bV2upQLUZfa3b+g\nlLtnmH4WiZS6jH3ZPPHREl76ejlN69bkxSuTObVbs3iHJVJk0RJJ1zx3au2/c0t3bYmUgq9TNzP0\n7TmsTt/D5f3bMPSsrtStWS3eYYkUS7RE0q1MohCpZLZnZPHQBwsYnZJGu8a1efPaoziqfaN4hyVy\nQKIlklVe2Fq8gJlZtDoi8qOJ89Zzz7vzSN+9j9+d2IGbTutEzWpV4x2WyAGLlkgmm9lY4F13X7W/\n0MyqA8cBVwKTgX/HLEKRCmLjzkzuGz+fD+aup3uLevzrqiM5rGX9eIclUmLREskA4GrgDTNrB2wD\nagJVgY+AJ919VmxDFCnf3J0x36fx1/cXsicrh9vP7MK1J7SnWtUq8Q5NpFREeyAxE3gWeNbMqgGN\ngT3uvq0sghMp71anZ/DncXP5aulmjmx7MMMv6kUHLXsrFUyRJ2109yxgXQxjEakwcnKdl79ZwWOT\nFlPF0IJTUqFp9l+RUrZo/Q7uensuM1dt46QuTXjogp60bKBJFqXiUiIRKSWZWTk88+lSXvhyGfVq\nVePJS3pzfp+Wmt5EKrwiJRIz6+7uC/KUneTun8ckKpFy5uvUzfx53FxWbsngor6tGHZONxrWrh7v\nsETKRFF7JKPN7FXgUYK7th4FkoGjYxWYSHmQvnsff31/AW/PWEPbRgfx2m/7c2zHxvEOS6RMFTWR\n9AdGEMwCXBd4DTg2VkGJJDp3Z9zMNTz43gJ2ZmZz/ckd+OMperBQKqeiJpIsYA9Qi6BHstzdc2MW\nlUgCW7llN8PGzWNK6mYOb9OA4Rf2okvzuvEOSyRuippIpgPvAkcCjYDnzWyQuw+KWWQiCSYrJ5d/\nfrWcpz5ZQrWqVXRLr0ioqI/WDnH3e909y93Xu/tAgsRSKDMbYGaLzSzVzIbm834NM3srfP9bM2sb\n8d5dYfliMzszovxmM5tvZvPM7A0z07JxEnMzV23lF3+fwoiJizi5S1M+ueVEfnV0WyUREYreI9lo\nZm3ylH1R2AFmVhUYBZwOpAHTzWx8nru/hgBb3b2jmQ0mGIe5xMy6A4MJVmc8BPjEzDoDzYE/Ad3d\nfY+ZjQ7r/buI7RApll17s3l80mJenrqCZnVr8vyvjuDMHs3jHZZIQilqInkfcIJ1SGoC7YDFBP/Q\nF6QfkOruywDM7E1gIBCZSAYC94XbY4CR4eJZA4E33X0vsNzMUsPzrQpjrmVmWcBBwNoitkGkWD5e\nsIF7353H+h2Z/PqoQ7ntzC5aK0QkH0VKJO7eM3LfzPoC10U5rCWwOmI/jeDur3zruHu2mW0nGINp\nCUzLc2xLd59qZo8TJJQ9wEfu/lF+H25m1wLXArRpk7czJVKwDTsy+cu785k4fz1dmtVl1OV96dvm\n4HiHJZKwDujJdnefYWZHRqmW38XjvOuWFFQn33IzO5igt7J/JuL/mtkV7v6ffGJ8AXgBIDk5Weul\nSFQ5uc7r367k0YmL2ZeTyx0DunDN8ZqlVySaoj7ZfkvEbhWgL7ApymFpQOuI/Vb8/DLU/jppZpYE\n1AfSCzn2NIJbjzeFcb0NHAP8LJGIFMe8NdsZNm4us9O2c2zHRjx0fk/aNq4d77BEyoWi9kgib5LP\nJhgzGRvlmOlAp3AdkzUEg+KX5akznmBxrKnAIOAzd3czGw+8bmZPEAy2dwK+A3KBo8zsIIJLW6cC\nKUVsg8jP7MzM4omPl/DyNytoWLs6Tw/uw3m9D9H8WCLFUNQxkvuLe+JwzOMGYBLBQlgvuft8M3sA\nSHH38cCLwKvhYHo6QbIhrDeaYGA+G7je3XOAb81sDDAjLJ9JePlKpDjcnQ/nref+CfPZuHMvV/QP\nBtPr19JgukhxWWHLrZvZBH4+rvE/7n5eLIIqbcnJyZ6Soo6LBFZtyeDe8fP4fPEmureox8MX9qRP\n6wbxDksk4ZjZ9+6eHK1etB7J46UUj0jc7cvO5R9fLeOZT5eSVMW499zu/ProQ0nSYLpIiURLJMvd\nfVWZRCISQ9OWbeHud+aRunEXZ/dszr3n9qB5fU2KIFIaoiWSdwju0MLMxrr7RbEPSaT0bNm1l4c/\nWMTYGWm0bliLf111JCd3bRrvsEQqlGiJJPLWlfaxDESkNOXmOqNTVvPIh4vI2BdM837DyZ2oVV3T\nvIuUtmiJxAvYFklYC9ftYNi4ucxYtY1+7Rry0PmH0amZpnkXiZVoiaS3me0g6JnUCrcJ993d68U0\nOpFi2L03m6c/XcqLU5ZTv1Y1Hr+4Nxf11ZrpIrFWaCJxd10HkITn7kyav4EHJsxn7fZMBh/ZmjsH\ndOVgrZkuUiYOaK4tkUSxYvNu7pswn88Xb6JLs7qMufRwkts2jHdYIpWKEomUS5lZOTw7OZXnvlhG\n9aQq3BM+E6IJFkXKnhKJlDufLNjAfRPmk7Z1DwP7HMKws7vRtJ6eCRGJFyUSKTdWbcng/gnz+XTR\nRjo1rcMb1xzF0R0axTsskUpPiUQSXmZWDs9/sYxnP0+lahXjz2d35TfHttNlLJEEoUQiCW3yoo38\nZfx8VqVncG6vFgw7pxst6teKd1giEkGJRBLS6vQMHnhvAR8v2ECHJrV57bf9ObZj43iHJSL5UCKR\nhLI3O4d/fLmMkZNTMYw7B3RlyHHtqJ6ky1giiUqJRBLGF0s2cd/4+SzfvJuzezbn7nO6c0gDXcYS\nSXRKJBJ3a7bt4a/vLeDDeetp17g2L1/djxM7N4l3WCJSREokEjeZWTm8OGU5Iz9LxXFuP7MLvz2+\nHTWSNDOPSHmiRCJlzt35ZOFGHnxvAavSMzizRzPuPqc7rRseFO/QROQAKJFImfph0y4emLCAL5Zs\nomPTOvxnSH+O66S7sUTKMyUSKRM7M7MY+VkqL329nJpJVTU3lkgFokQiMZWb67wzaw2PfLiITTv3\ncvERrbhjQFea1K0R79BEpJQokUjMzFuznXvfnceMVdvo3ao+//h1Mn1aN4h3WCJSypRIpNRt2bWX\nxz9awpvTV9GodnUeHdSLQX1bUaWKVioUqYiUSKTUZOfk8p9pK3ni4yXs3pfD1ce248bTOlGvZrV4\nhyYiMaREIqVi6g9buH/CfBat38mxHRtx3y960KlZ3XiHJSJlQIlESmTttj089MFC3p+zjpYNavHc\nFX05s0dzzHQZS6SyUCKRA5KZFUyuOOrzVNzh5tM6c92J7alZTU+li1Q2SiRSLO7O+3PX8cgHi1iz\nbQ9nHdacYed0o9XBeipdpLJSIpEim7dmOw9MWMB3K9Lp1qIej1/cW0vdiggxfazYzAaY2WIzSzWz\nofm8X8PM3grf/9bM2ka8d1dYvtjMzgzLupjZrIjXDjO7KZZtENi4M5Pb/zubX4ycwg+bdvHwBT15\n74/HKYmICBDDHomZVQVGAacDacB0Mxvv7gsiqg0Btrp7RzMbDIwALjGz7sBgoAdwCPCJmXV298VA\nn4jzrwHGxaoNlV1mVg4vfb2cUZ+lsi8nl2uOb88Np3TU7bwi8hOxvLTVD0h192UAZvYmMBCITCQD\ngfvC7THASAtu9xkIvOnue4HlZpYanm9qxLGnAj+4+8oYtqFScncmzV/PQx8sZHX6Hk7r1oxh53Sj\nXePa8Q5NRBJQLBNJS2B1xH4a0L+gOu6ebWbbgUZh+bQ8x7bMc+xg4I3SDFhgwdodPPDefKYtS6dz\nszq8OqQfx3fSIlMiUrBYJpL8HiTwItYp9Fgzqw6cB9xV4IebXQtcC9CmTZtosVZ6m3ft5W8fLeGt\n6auoX6saDw7swaX92pCk2XlFJIpYJpI0oHXEfitgbQF10swsCagPpBfh2LOAGe6+oaAPd/cXgBcA\nkpOT8yYwCe3LzuXlb1bwzKdL2ZOVw5XHtOWmUztT/yCNg4hI0cQykUwHOplZO4JB8cHAZXnqjAeu\nJBj7GAR85u5uZuOB183sCYLB9k7AdxHHXYoua5XI/lUKH3p/ASu2ZHBylyYMO6c7HZvWiXdoIlLO\nxCyRhGMeNwCTgKrAS+4+38weAFLcfTzwIvBqOJieTpBsCOuNJhiYzwaud/ccADM7iOBOsOtiFXtF\nt3j9Th58bwFTUjfToUlt/v2bIzmpS9N4hyUi5ZS5V/yrPsnJyZ6SkhLvMOJu0869PPFxMA5Sp0YS\nN5/emSuO0iqFIpI/M/ve3ZOj1dOT7ZVAZlYOL05ZzrOTU9mbncuvj27Ljad24uDa1eMdmohUAEok\nFVhurvPu7DU8NnExa7dnckb3Zgw9qyvtm2gcRERKjxJJBfXtsi089MFC5qRt57CW9Xjikj4c1V5T\nmohI6VMiqWCWb97N8A8XMmnZioAuAAAOq0lEQVT+BlrUr8kTv+zN+X1aaplbEYkZJZIKYlvGPp7+\ndCmvTl1JjaQq3HZGZ4Yc155a1bU+iIjElhJJObcvO5dXpgYPFO7am80lR7bm5tM707RuzXiHJiKV\nhBJJOeXuTJy3nuETF7FySwYndG7CsLO70aW51kkXkbKlRFIOzVq9jYfeX8D0FVvp0qwuL1/djxM7\na2JFEYkPJZJyJG1rBo9NWsy7s9bSuE4NHrmwJxcf0UoTK4pIXCmRlAPbM7IY9Xkq//5mBQb88ZSO\nXHdiB+rU0NcnIvGnf4kSWGZWDq9MXcHIz1LZuTebi/q24pbTO3NIg1rxDk1E5H+USBJQbq7zzqw1\n/O2jJazZtoeTujThzgFd6daiXrxDExH5GSWSBPPlkk0M/3ARC9btoGfL+jw2qBfHdGwc77BERAqk\nRJIg5q3ZzoiJi/hq6WZaN6zF04P78Iteh+iJdBFJeEokcZa2NYO/fbSEcTPX0OCgatxzbneuOKoN\nNZL0RLqIlA9KJHGyLWMfoyan8vI3KzGD35/Ugd+d2IH6tbTErYiUL0okZSwzK4eXv1nBqMnBnViD\n+rbiljM606K+7sQSkfJJiaSM5L0T6+QuTbjzrK50ba47sUSkfFMiiTF358ulmxn+4SIW7r8T6+Je\nHNNBd2KJSMWgRBJDM1Zt5dGJi5i2LJ3WDWvxzKWHc27PFroTS0QqFCWSGFi6YSePTVrMRws20LhO\nde4/rweX9mtD9STNiSUiFY8SSSlK25rBU58s5e0ZadSunsStp3fm6uPaUVtzYolIBaZ/4UrBll17\nGTX5B/4zbSUYDDmuHb8/qSMNa1ePd2giIjGnRFICu/Zm88+vlvGPL5exJyuHi49ozY2nddKkiiJS\nqSiRHIC92Tm8Nm0VIyenkr57H2cd1pxbz+hMx6ZanVBEKh8lkmLIyXXGzVzDkx8Hz4Ic06ERdwzo\nSp/WDeIdmohI3CiRFIG78/GCDTw2aTFLN+6iZ8v6jLioF8d10rMgIiJKJFFMW7aFERMXMXPVNto3\nrs2zl/flrMOaY6ZnQUREQImkQDm5zm9fns7kxZtoXq8mwy/sySCtjy4i8jNKJAWoWsVo36QOR3do\nxK+PbkvNaprWXUQkPzH9eW1mA8xssZmlmtnQfN6vYWZvhe9/a2ZtI967KyxfbGZnRpQ3MLMxZrbI\nzBaa2dGxiv+ec7tz7QkdlERERAoRs0RiZlWBUcBZQHfgUjPrnqfaEGCru3cEngRGhMd2BwYDPYAB\nwLPh+QCeBia6e1egN7AwVm0QEZHoYtkj6Qekuvsyd98HvAkMzFNnIPByuD0GONWCUeyBwJvuvtfd\nlwOpQD8zqwecALwI4O773H1bDNsgIiJRxDKRtARWR+ynhWX51nH3bGA70KiQY9sDm4B/mdlMM/un\nmdWOTfgiIlIUsUwk+d0f60WsU1B5EtAX+D93PxzYDfxs7AXAzK41sxQzS9m0aVPRoxYRkWKJZSJJ\nA1pH7LcC1hZUx8ySgPpAeiHHpgFp7v5tWD6GILH8jLu/4O7J7p7cpEmTEjZFREQKEstEMh3oZGbt\nzKw6weD5+Dx1xgNXhtuDgM/c3cPyweFdXe2ATsB37r4eWG1mXcJjTgUWxLANIiISRcyeI3H3bDO7\nAZgEVAVecvf5ZvYAkOLu4wkGzV81s1SCnsjg8Nj5ZjaaIElkA9e7e0546j8Cr4XJaRnwm1i1QURE\norOgA1CxJScne0pKSrzDEBEpV8zse3dPjlqvMiQSM9sErIx3HMXUGNgc7yDKmNpcOajN5ceh7h51\nkLlSJJLyyMxSivJLoCJRmysHtbni0QyEIiJSIkokIiJSIkokieuFeAcQB2pz5aA2VzAaIxERkRJR\nj0REREpEiaSMmVnVcMLJ98J9M7OHzGxJuL7KnyLKnwnXZJljZn0jznGlmS0NX1cW9FmJIp82n2pm\nM8xslplNMbOOYXmx16dJRGa2wszmhu1LCcsamtnH4Xf2sZkdHJZXiO+5gDY/Fq4bNMfMxplZg4j6\nBa03VOgaRokkvzZHvHebmbmZNQ73K8T3XCB316sMX8AtwOvAe+H+b4BXgCrhftPwz7OBDwkmsDwK\n+DYsb0jwRH9D4OBw++B4t6uYbV4CdAu3/wD8O2L7uXB7MPBWuN0dmA3UANoBPwBV492uQtq7Amic\np+xRYGi4PRQYUZG+5wLafAaQFG6PiGhzvt9n+PqBYJbv6mGd7vFuW3HaHJa3JpjRY+X+9yvK91zQ\nSz2SMmRmrYBzgH9GFP8eeMDdcwHcfWNYPhB4xQPTgAZm1gI4E/jY3dPdfSvwMcHiXwmpgDY7UC/c\nrs+Pk3kWa32aWMdeyiLb9jJwfkR5uf+e8+PuH3mwPATANILJV6Hg77MoaxiVB08Cd/DT2c4r7PcM\nurRV1p4i+AuWG1HWAbgknPL+QzPrFJYXtCZLUdZ5SST5tfm3wAdmlgb8Chgelhd3fZpE5cBHZva9\nmV0bljVz93UA4Z9Nw/KK8j3n1+ZIVxP8IocK3GYzOw9Y4+6z89StKG3OV8wmbZSfMrNzgY3u/r2Z\nnRTxVg0g092TzexC4CXgeIq/VkvCKaTNNwNnu/u3ZnY78ARBcin3bQ4d6+5rzawp8LGZLSqkboVt\ns7t/CWBmwwgmX30trFtQ2/L7YVuu2gwMI7ikl1dF+Z7zpR5J2TkWOM/MVhB02U8xs/8Q/AIZG9YZ\nB/QKtwtbkyXaOi+JIr82vw/09h/XlHkLOCbcLu76NAnJ3deGf24k+E77ARvCSxmEf+6/hFkRvueC\n2kw4eHwucLmHgwJU3DafSDDmMzv8O98KmGFmzakgbS5QvAdpKuMLOIkfB56HA1dHlE8Pt8/hp4Nz\n34XlDYHlBANzB4fbDePdpqK2maAXvBnoHJYPAcaG29fz08H20eF2D346OLuMBB1sB2oDdSO2vyG4\n5v0YPx1sf7SifM+FtHkAwVIQTfLUz/f7DP9uLAvL9g+294h3+4rT5jx1VvDjYHu5/54Le+nSVvwN\nJ1hf5WZgF8ElHoAPCO70SAUyCNddcfd0M3uQYOEwCAbq08s25APnwTo11wBjzSwX2Epw/RwObH2a\nRNMMGBfcI0AS8Lq7TzSz6cBoMxsCrAIuDutXhO+5oDanEiSLj8P3prn77wr7Pi2fNYzKvjlFkm+b\nC6lfEb7nAunJdhERKRGNkYiISIkokYiISIkokYiISIkokYiISIkokYiISIkokUhCMbMnzeymiP1J\nZvbPiP2/mdktpfyZu0rzfOE5+5jZ2RH795nZbUU4zszsMzOrF1F2QTiTbNcYxNnWzC4r7fNGnP8G\nM/tNrM4viUGJRBLNN4RPuptZFaAxwQNs+x0DfB2HuIqrD8FzA8V1NjDb3XdElF0KTCF8rqaUtQXy\nTSTh7AIl9RLwp1I4jyQwJRJJNF/z45QpPYB5wE4zO9jMagDdgJlmVsfMPrVgXZO5ZjYQwMxGmNkf\n9p8s7AncGm7fbmbTw/Ug7s/vw/OrE/5qX2hm/zCz+Wb2kZnVCt87Mqw71YL1N+aZWXXgAYLJOGeZ\n2SXh6bub2edmtszCdWfycTnwbkQ8dQimmhlCRCIxs5PCc42xYM2P18KZkjGzs8OyKRasgbF/HZgT\nw3hmWbA+TF2CB2KPD8tuNrOrzOy/ZjaBYEJCi2jX3P1tCT//CzMbbcFaOsPN7HIz+y6s1wHA3TOA\nFWZW3mZrluKI96P1eumV90UwtUQb4Drgd8CDBL/UjwW+DOskAfXC7cYETwwbcDjwRcS5FoTnOoNg\n3Wwj+AH1HnBCWGdX+Ge+dQh+tWcDfcJ6o4Erwu15wDHh9nBgXrh9FTAyIo77CHpbNcJ4twDV8mn7\nSsKpN8L9K4AXw+1vgL7h9kkEsyO3CmOdChwH1CSYTbZdWO8NfpyOZwLBRIMAdcL/hiftfz8i7jTC\naTqAiwimNq9K8DT3KqBFeNy2cLsGsAa4PzzmRuCpiHMOA26N998rvWL3Uo9EEtH+XskxBP9ATo3Y\n/yasY8DDZjYH+IRg6u1m7j4TaGpmh5hZb2Cru68iSBJnADOBGUBXoBM/VVid5e4+K9z+HmhrwYp/\ndd19f0yvR2nX+x6swbGZYNLGZvnUaejuOyP2LyWY8JLwz0sj3vvO3dM8WMtmFkHC6wos82CdDwgS\nyX5fA0+EvaEG/uNaIXl97D9O03Ec8Ia757j7BuAL4Mjwvenuvs7d9xIsSPVRWD43jGW/jcAhBXyW\nVACaa0sS0f5xkp4Ev/hXA7cCOwiuuUNwCagJcIS7Z4WzrdYM3xsDDAKa8+M/wgY84u7PF/K5+dax\nYMnfvRFFOUAt8p8CvDB5z5Hf/3/ZZlbF3XPNrBFwCnCYmTlBr8DN7I5CzldgTO4+3ILZl88GppnZ\naQVU3R2xXVgbIz8/N2I/l5+2rSawp5DzSDmnHokkoq8Jph5PD38JpwMNgKMJeicQTDG/MUwiJwOH\nRhz/JsF4wiCCpALBRIBXh2MOmFlLC9aRiFSUOv/jwYp2O83sqLAocjB8J1C3OI0OLSZYapYw/lfc\n/VB3b+vurQlmhz2ukOMXAe3tx/Xu94/PYGYd3H2uu48AUgh6L9Hi/JJgrKeqmTUhuNT3XTHb1Jng\nB4FUUEokkojmEowjTMtTtj28LATBIknJZpZC0Dv53+JRHswYW5dgpbr9qxJ+RHDpaaqZzSVIMD/5\nB7QodfIxBHjBzKYS/HrfHpZPJhhcjxxsL4r3CcYfILiMNS7P+2Mp4C6rsA17gD8AE81sCrAhIqab\nwkHz2QQ9hA+BOQS9oNkWzECd17iwzmzgM+AOd19fjPZAMLb1STGPkXJEs/+KlICZ1XH3XeH2UKCF\nu99YgvO1IOiFnF7SmMK7uEYBS939yQM9X0mY2eHALe7+q3h8vpQN9UhESuacsNcxj2CJ5L+W5GRh\nD+ofFvFA4gG4xsxmAfMJLgEWNi4Ua42Be+L4+VIG1CMREZESUY9ERERKRIlERERKRIlERERKRIlE\nRERKRIlERERKRIlERERK5P8Bphz/C+p3oQoAAAAASUVORK5CYII=\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe3fbe2470>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"sp.plot(bigbox.waveset, flux_unit=units.PHOTNU)" | |
] | |
}, | |
{ | |
"cell_type": "markdown", | |
"metadata": {}, | |
"source": [ | |
"This is to mimic reported calculation." | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 10, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [ | |
"input_spectrum = sp.normalize(18 * u.ABmag, band=bigbox)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 13, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAERCAYAAABowZDXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAG0RJREFUeJzt3X2UZVV95vHvI69qg2C6QQXaVkYD\nqBG0Q4ytI2oCxNckY1ZEQwySkBdfQInRmDUrviQTDDNoMiZDiDBoBiFEIL6gQGtQAgLSDQ0NtJhO\nQ5A0kSatAY2JNv2bP84uuZS36p5u6lZVF9/PWnfVufvsu+/e3KKfOmffs0+qCkmSRnnUXHdAkrRj\nMDAkSb0YGJKkXgwMSVIvBoYkqRcDQ5LUy4ILjCRnJbknyc0z1N4lSb6V5DOTypPkD5N8Lcm6JG+d\nifeTpPlqwQUGcDZw9Ay2dypw7JDyXwEOAA6qqoOB82bwPSVp3llwgVFVVwCbB8uSHNiOFFYn+fsk\nB21De18A7h+y6zeB91XV1lbvnofTb0ma7xZcYEzhDOAtVfVc4LeBP5+BNg8EfjHJqiSfS/K0GWhT\nkuatnee6A+OWZBHwfOBvkkwU79b2/TzwviEv++eqOmpE07sB/1FVy1s7ZwEvnJleS9L8s+ADg+4o\n6ltVdejkHVV1IXDhdrZ7F3BB274I+L/b2Y4k7RAW/CmpqroPuD3JL8APvt307Blo+m+Bl7TtFwFf\nm4E2JWneykJbrTbJucARwGLgG8DvA38H/B/gicAuwHlVNexU1LD2/h44CFgE/CtwfFVdmmQv4Bxg\nKfBt4Deq6saZHY0kzR8LLjAkSeOx4E9JSZJmxoKa9F68eHEtW7ZsrrshSTuM1atX31tVS/rUXVCB\nsWzZMlatWjXX3ZCkHUaSf+pb11NSkqRexhYYSQ5IcnlbmO+WJCcOqbN3kouS3JTkK0meObDv6CS3\nJVmf5F3j6qckqZ9xHmFsAU5uC/M9D3hTkkMm1Xk3sKaqfgz4ZeBPAJLsBPwZ8DPAIcAxQ14rSZpF\nYwuMqrq7qq5v2/cD64D9JlU7BPhCq/NVYFmSfYHDgfVVtaGqvke3Euyrx9VXSdJoszKHkWQZcBhw\n7aRdNwI/3+ocDjwZ2J8uWL4+UO8ufjhsJto+oS0AuGrTpk0z23FJ0g+MPTDa4n8XACe1ZToGnQLs\nnWQN8BbgBrpTWeGHDb3CsKrOqKrlVbV8yZJe3wyTJG2HsX6tNskudGFxTlvo7yFagBzX6ga4vT0e\nQ3dzogn7AxvH2VdJ0vTG+S2pAGcC66rqtCnq7JVk1/b0V4ErWohcBzwtyVPa/tcCnxpXXyVJo43z\nCGMF3a1N17ZTTtB9K2opQFWdDhwMfCzJA8CtwPFt35YkbwYuBXYCzqqqW8bYV0nSCGMLjKq6kuFz\nEYN1rgaG3qmuqj4LfHYMXZMkbQev9JYk9WJgSJJ6MTAkSb0YGJKkXgwMSVIvBoYkqRcDQ5LUi4Eh\nSerFwJAk9WJgSJJ6MTAkSb0YGJKkXgwMSVIvBoYkqRcDQ5LUi4EhSerFwJAk9WJgSJJ6MTAkSb0Y\nGJKkXgwMSVIvBoYkqRcDQ5LUi4EhSerFwJAk9WJgSJJ6MTAkSb2MLTCSHJDk8iTrktyS5MQhdR6X\n5NNJbmx1jhvY90CSNe3xqXH1U5LUz85jbHsLcHJVXZ9kD2B1kpVVdetAnTcBt1bVK5MsAW5Lck5V\nfQ/4blUdOsb+SZK2wdiOMKrq7qq6vm3fD6wD9ptcDdgjSYBFwGa6oJEkzTOzMoeRZBlwGHDtpF0f\nBg4GNgJrgROramvbt3uSVUmuSfKz07R9Qqu3atOmTTPfeUkSMAuBkWQRcAFwUlXdN2n3UcAa4EnA\nocCHk+zZ9i2tquXA64APJTlwWPtVdUZVLa+q5UuWLBnPICRJ4w2MJLvQhcU5VXXhkCrHARdWZz1w\nO3AQQFVtbD83AF+kO0KRJM2RcX5LKsCZwLqqOm2KancCL2319wV+FNiQZO8ku7XyxcAK4NYp2pAk\nzYJxfktqBXAssDbJmlb2bmApQFWdDrwfODvJWiDAO6vq3iTPB/4iyVa6UDtl0rerJEmzbGyBUVVX\n0oXAdHU2AkcOKf8y8KwxdU2StB280luS1IuBIUnqxcCQJPViYEiSejEwJEm9GBiSpF4MDElSLwaG\nJKkXA0OS1IuBIUnqxcCQJPViYEiSejEwJEm9GBiSpF4MDElSLwaGJKkXA0OS1IuBIUnqxcCQJPVi\nYEiSejEwJEm9GBiSpF4MDElSLwaGJKkXA0OS1IuBIUnqxcCQJPUytsBIckCSy5OsS3JLkhOH1Hlc\nkk8nubHVOW5g3xuS/EN7vGFc/ZQk9bPzGNveApxcVdcn2QNYnWRlVd06UOdNwK1V9cokS4DbkpwD\nLAJ+H1gOVHvtp6rqm2PsryRpGmM7wqiqu6vq+rZ9P7AO2G9yNWCPJKELic10QXMUsLKqNreQWAkc\nPa6+SpJGm5U5jCTLgMOAayft+jBwMLARWAucWFVb6YLl6wP17uKHw2ai7ROSrEqyatOmTTPcc0nS\nhLEHRpJFwAXASVV136TdRwFrgCcBhwIfTrInkCFN1bD2q+qMqlpeVcuXLFkygz2XJA0aa2Ak2YUu\nLM6pqguHVDkOuLA664HbgYPojigOGKi3P91RiCRpjozzW1IBzgTWVdVpU1S7E3hpq78v8KPABuBS\n4MgkeyfZGziylUmS5sjIb0kl2R14BfBCulNH3wVuBi6uqlumeekK4FhgbZI1rezdwFKAqjodeD9w\ndpK1dKeh3llV97b3fT9wXXvd+6pq8zaOTZI0g6YNjCTvAV4JfJFuwvoeYHfg6cApLUxOrqqbJr+2\nqq5k+FzEYJ2NdEcPw/adBZw1cgSSpFkx6gjjuqp6zxT7TkuyD+2IQZK0sE0bGFV18Yj999AddUiS\nFrhRp6Q+Nd3+qnrVzHZHkjRfjTol9ZN0F9CdSzeHMe2chCRp4RoVGE8Afho4BngdcDFw7ohvR0mS\nFqBpr8Ooqgeq6pKqegPwPGA98MUkb5mV3kmS5o0+12HsBryc7ihjGfCnwLCrtiVJC9ioSe+PAs8E\nPge8t6punpVeSZLmnVFHGMcC36G7UO+t3WofQDf5XVW15xj7JkmaR0Zdh+EtXCVJwHYsPpjksUle\nn2Tai/okSQtLr8BIsmuSn01yPnA38FPA6WPtmSRpXhk16T1xDcZRwOXAXwGHV9Vxs9A3SdI8MmrS\n+1Lg74EXVNXtAEn+ZOy9kiTNO6MC47nAa4HPJ9kAnAfsNPZeSZLmnVFXet9QVe+sqgOB9wCHAbsm\n+VySE2ajg5Kk+aH3t6Sq6qqqejOwH/Ah4DVj65Ukad7Z5q/VVtXWqroUOGgM/ZEkzVNemCdJ6sXA\nkCT1Muo6jP8N1LBdwF5j6ZEkaV4a9bXaVdu5T5K0wIwKjHOqasus9ESSNK+NmsP4ysRGOz0lSXqE\nGnWEkYHtFePsyFw68bwb+N6WrXPdDUnaLnvuvgsfeM2Pjf19RgXGsAnvBef2e7/Df3z/gbnuhiRt\nl70es+usvM+owDgoyU10RxoHtm148I57U0ZakgOAjwFPALYCZ1TVn0yq8w7g9QN9ORhYUlWbk9wB\n3A88AGypquXbNLJt8Kk3v2BcTUvSgjEqMA5+GG1vAU6uquuT7AGsTrKyqm6dqFBVpwKnAiR5JfC2\nqto80MaLq+reh9EHSdIMGRUYd1bVtKelkmRYnaq6m+5mS1TV/UnW0a1Ddevkus0xwLmjuyxJmguj\nviV1eZK3JFk6WNjuwPeSJB8F3jDqTZIso1vp9top9j8GOBq4YKC4gMuSrHZlXEmae6OOMI4G3gic\nm+QpwLeA3enuiXEZ8MGqWjNdA0kW0QXBSVV13xTVXglcNel01Iqq2phkH2Blkq9W1RVD2j8BOAFg\n6dKlk3dLkmZIRpxxerBisguwGPhuVX1rG17zGeDSqjptmnoXAX9TVR+fYv97gG9X1f+c7v2WL19e\nq1Z5Abok9ZVkdd8vFW3L/TC+X1V3b0NYBDgTWDciLB4HvAj45EDZY9tEOUkeCxwJ3Ny3r5KkmTfq\nlNTDsQI4FlibZOK01buBpQBVdXor+zngsqr6zsBr9wUu6jKHnYGPV9UlY+yrJGmEsQVGVV3JQ68U\nn6re2cDZk8o2AM8eS8ckSdul1ympJIcMKTtixnsjSZq3+s5hnJ/knek8ui1E+Efj7JgkaX7pGxg/\nARwAfBm4DtjIAl6MUJL0w/oGxveB7wKPprsO4/aqcnlXSXoE6RsY19EFxo8DLwCOSfKJsfVKkjTv\n9P2W1PFVNXFF3L8Ar05y7Jj6JEmah/oGxj2T15MCvjTTnZEkzV99A+NiusUAQzeH8RTgNuAZY+qX\nJGme6RUYVfWswedJngP8+lh6JEmal3qvJTWoqq6nmwCXJD1C9DrCSPL2gaePAp4DbBpLjyRJ81Lf\nOYw9Bra30M1pXDBFXUnSAtR3DuO94+6IJGl+mzYwknya7ttRQ1XVq2a8R5KkeWnUEca0d7iTJD1y\njAqM26vqzlnpiSRpXhv1tdq/ndhI4iS3JD2CjQqMwTvmPXWcHZEkzW+jAqOm2JYkPcKMmsN4dpL7\n6I40Ht22ac+rqvYca+8kSfPGtIFRVTvNVkckSfPbdq0lJUl65DEwJEm9GBiSpF4MDElSLwaGJKkX\nA0OS1MvYAiPJAUkuT7IuyS1JThxS5x1J1rTHzUkeSPL4tu/oJLclWZ/kXePqpySpn3EeYWwBTq6q\ng4HnAW9Kcshghao6taoOrapDgd8FvlRVm5PsBPwZ8DPAIcAxk18rSZpdYwuMqrq73fubqrofWAfs\nN81LjgHObduHA+urakNVfQ84D3j1uPoqSRptVuYwkiwDDgOunWL/Y4CjefC2r/sBXx+ochdThE2S\nE5KsSrJq0yZvMy5J4zL2wEiyiC4ITqqq+6ao9krgqqraPPGyIXWGLn5YVWdU1fKqWr5kyZKH32FJ\n0lBjDYwku9CFxTlVdeE0VV/Lg6ejoDuiOGDg+f7AxpnvoSSpr3F+SyrAmcC6qjptmnqPA14EfHKg\n+DrgaUmekmRXukD51Lj6KkkabdTy5g/HCuBYYG2SNa3s3cBSgKo6vZX9HHBZVX1n4oVVtSXJm4FL\ngZ2As6rqljH2VZI0wtgCo6quZPhcxOR6ZwNnDyn/LPDZGe+YJGm7eKW3JKkXA0OS1IuBIUnqxcCQ\nJPViYEiSejEwJEm9GBiSpF4MDElSLwaGJKkXA0OS1IuBIUnqxcCQJPViYEiSejEwJEm9GBiSpF4M\nDElSLwaGJKkXA0OS1IuBIUnqxcCQJPViYEiSejEwJEm9GBiSpF4MDElSLwaGJKkXA0OS1IuBIUnq\nZWyBkeSAJJcnWZfkliQnTlHviCRrWp0vDZTfkWRt27dqXP2UJPWz8xjb3gKcXFXXJ9kDWJ1kZVXd\nOlEhyV7AnwNHV9WdSfaZ1MaLq+reMfZRktTT2I4wquruqrq+bd8PrAP2m1TtdcCFVXVnq3fPuPoj\nSXp4ZmUOI8ky4DDg2km7ng7sneSLSVYn+eWBfQVc1spPmKbtE5KsSrJq06ZNM911SVIzzlNSACRZ\nBFwAnFRV9w15/+cCLwUeDVyd5Jqq+hqwoqo2ttNUK5N8taqumNx+VZ0BnAGwfPnyGudYJOmRbKxH\nGEl2oQuLc6rqwiFV7gIuqarvtLmKK4BnA1TVxvbzHuAi4PBx9lWSNL1xfksqwJnAuqo6bYpqnwRe\nmGTnJI8BfgJYl+SxbaKcJI8FjgRuHldfJUmjjfOU1ArgWGBtkjWt7N3AUoCqOr2q1iW5BLgJ2Ap8\npKpuTvJU4KIuc9gZ+HhVXTLGvkqSRhhbYFTVlUB61DsVOHVS2QbaqSlJ0vzgld6SpF4MDElSLwaG\nJKkXA0OS1IuBIUnqxcCQJPViYEiSejEwJEm9GBiSpF4MDElSLwaGJKkXA0OS1IuBIUnqxcCQJPVi\nYEiSejEwJEm9GBiSpF4MDElSLwaGJKkXA0OS1IuBIUnqxcCQJPViYEiSejEwJEm9GBiSpF4MDElS\nLwaGJKmXVNVc92HGJNkE/NNc92MbLQbunetOzDLH/MjgmHcMT66qJX0qLqjA2BElWVVVy+e6H7PJ\nMT8yOOaFx1NSkqReDAxJUi8Gxtw7Y647MAcc8yODY15gnMOQJPXiEYYkqRcDQ5LUi4ExJkl2SnJD\nks+050nyh0m+lmRdkrcOlP9pkvVJbkrynIE23pDkH9rjDXM1lr6GjPmlSa5PsibJlUn+SyvfLclf\ntzFfm2TZQBu/28pvS3LU3IyknyR3JFnbxreqlT0+ycr2ma1MsncrXxCf8xRjPjXJV9u4Lkqy10D9\noZ9nkqNb2fok75qLsfQ1bMwD+347SSVZ3J4viM95SlXlYwwP4O3Ax4HPtOfHAR8DHtWe79N+vgz4\nHBDgecC1rfzxwIb2c++2vfdcj2sbx/w14OC2/VvA2QPbp7ft1wJ/3bYPAW4EdgOeAvwjsNNcj2ua\n8d4BLJ5U9sfAu9r2u4APLKTPeYoxHwns3LY/MDDmoZ9ne/wj8FRg11bnkLke27aMuZUfAFxKd7Hw\n4oX0OU/18AhjDJLsD7wc+MhA8W8C76uqrQBVdU8rfzXwsepcA+yV5InAUcDKqtpcVd8EVgJHz9og\nttEUYy5gz7b9OGBj23418NG2/QngpUnSys+rqv+sqtuB9cDh4+77DBsc20eBnx0o3+E/52Gq6rKq\n2tKeXgPs37an+jwPB9ZX1Yaq+h5wXqu7o/kg8Dt0v+cTFuznDJ6SGpcP0f0ibR0oOxD4xSSrknwu\nydNa+X7A1wfq3dXKpiqfr4aN+VeBzya5CzgWOKWV/2Bs7R+afwN+hB1vzAVclmR1khNa2b5VdTdA\n+7lPK18on/OwMQ96I91f2LCAx5zkVcA/V9WNk+oulDEPtfNcd2ChSfIK4J6qWp3kiIFduwH/UVXL\nk/w8cBbwQrpD18lqmvJ5Z5oxvw14WVVdm+QdwGl0IbLDj7lZUVUbk+wDrEzy1WnqLtgxV9UVAEl+\nD9gCnNPqTjW2YX+o7lBjBn6P7lTcZAvlcx7KI4yZtwJ4VZI76A61X5Lk/9H9RXFBq3MR8GNt+y66\nc6ET9qc7dTNV+Xw0bMwXA8+uqmtbnb8Gnt+2fzC2JDvTna7azI41ZqpqY/t5D91nejjwjXYKgvZz\n4tTjQvicpxozbRL3FcDrq520Z+GO+UV0czI3tt/5/YHrkzyBBTLmKc31JMpCfgBH8OAE8CnAGwfK\nr2vbL+ehk2RfaeWPB26nmyDbu20/fq7H1HfMdEev9wJPb+XHAxe07Tfx0Env89v2M3joJOkG5umk\nN/BYYI+B7S/TnZM+lYdOev/xQvmcpxnz0cCtwJJJ9Yd+nu13Y0Mrm5j0fsZcj29bxjypzh08OOm9\nw3/O0z08JTV7TgHOSfI24Nt0p2YAPkv3zYr1wL/TfZuKqtqc5P3Ada3e+6pq8+x2eftV1ZYkvwZc\nkGQr8E2689sAZwJ/lWQ93ZHFa9trbklyPt0/PluAN1XVA7Pf+172BS7q5urZGfh4VV2S5Drg/CTH\nA3cCv9DqL4TPeaoxr6cLhZVt3zVV9RvTfZ5J3kz3DaOdgLOq6pbZH04vQ8c8Tf2F8DlPyaVBJEm9\nOIchSerFwJAk9WJgSJJ6MTAkSb0YGJKkXgwMzYkkH0xy0sDzS5N8ZOD5/0ry9hl+z2/PZHutzUOT\nvGzg+XuS/HaP1yXJ3yXZc6Ds59rKpweNoZ/LkrxuptsdaP/NSY4bV/uaHwwMzZUv0678TvIoYDHd\nhV4Tng9cNQf92laH0n3vflu9DLixqu4bKDsGuJJ2XcoMWwYMDYx2tf3DdRbw1hloR/OYgaG5chUP\nLhXyDOBm4P4keyfZDTgYuCHJoiRfSHdfjbVJXg2Q5ANJfmuisfaX/clt+x1Jrmv3I3jvsDcfVqf9\nFb4uyV8muSXJZUke3fb9eKt7dbr7P9ycZFfgfXSLSq5J8out+UOSfDHJhrT7ngzxeuCTA/1ZRLfE\nyvEMBEaSI1pbn0h3z4lz2sq+JHlZK7sy3T0YJu5D8qLWnzXp7k+yB92Foy9sZW9L8itJ/ibJp+kW\n1svAuNZOjKW9/5eSnJ/uXi6nJHl9kq+0egcCVNW/A3ck2dFWF9a2mOtLzX08ch90SyosBX4d+A3g\n/XR/ea8Armh1dgb2bNuL6a6gDXAY8KWBtm5tbR0JnNHqPIpumZL/2up8u/0cWofur/AtwKGt3vnA\nL7Xtm4Hnt+1TgJvb9q8AHx7ox3vojp52a/39V2CXIWP/J9qSE+35LwFntu0vA89p20fQrea7f+vr\n1cALgN3pVj99Sqt3Lg8uQ/NpugXzABa1/4ZHTOwf6PddtOUpgP9Gt+T2TnRXN98JPLG97lttezfg\nn4H3ttecCHxooM3fA06e698rH+N7eIShuTRxlPF8un8Irx54/uVWJ8D/SHIT8Hm6JaH3raobgH2S\nPCnJs4FvVtWddGFwJHADcD1wEPA0Hmq6OrdX1Zq2vRpYlu4OcntU1USfPj5iXBdXdw+Ie+kWH9x3\nSJ3HV9X9A8+PoVu4kfbzmIF9X6mqu6q7l8oaumA7CNhQ3X0moAuMCVcBp7Wjm73qwXtVTLayHlye\n4gXAuVX1QFV9A/gS8ONt33VVdXdV/SfdjY8ua+VrW18m3AM8aYr30gLgWlKaSxPzGM+i+wv+68DJ\nwH1058ShO3WzBHhuVX2/rQ66e9v3CeA1wBN48B/bAH9UVX8xzfsOrZPuVrH/OVD0APBohi9NPZ3J\nbQz7/2xLkkdV1dYkPwK8BHhmkqL7K7+S/M407U3Zp6o6Jd1qwS8DrknyU1NU/c7A9nRjHHz/rQPP\nt/LQse0OfHeadrSD8whDc+kquiWxN7e/bDcDewE/SXe0Ad3S5/e0sHgx8OSB159Hd77/NXThAd2C\ndm9scwIk2S/dfQwG9anzA9XdIe3+JM9rRYOT0vcDe2zLoJvb6G5RSuv/x6rqyVW1rKoOoFvN9AXT\nvP6rwFPz4P3QJ+ZPSHJgVa2tqg8Aq+iORkb18wq6uZidkiyhO0X3lW0c09Ppgl8LlIGhubSW7jz/\nNZPK/q2dzoHuZjzLk6yiO9r4wU2KqlvhdA+6O59N3OXuMrpTRlcnWUsXJA/5h7JPnSGOB85IcjXd\nX+P/1sovp5vkHpz07uNiuvkB6E4/XTRp/wVM8a2mNobv0t0b/ZIkVwLfGOjTSW3y+ka6v/g/B9xE\nd1RzY7oVkye7qNW5Efg74Heq6l+2YTzQzT19fhtfox2Iq9VKPSRZVFXfbtvvAp5YVSc+jPaeSHdU\n8dMPt0/tW1N/BvxDVX1we9t7OJIcBry9qo6di/fX7PAIQ+rn5e0o4ma6W+v+wcNprB0R/WUGLtzb\nDr+WZA1wC92pu+nmbcZtMfDf5/D9NQs8wpAk9eIRhiSpFwNDktSLgSFJ6sXAkCT1YmBIknr5/xRf\nUbf/hwNmAAAAAElFTkSuQmCC\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe3fc11c18>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"# Still flat in FLAM space.\n", | |
"input_spectrum.plot(bigbox.waveset, flux_unit=units.FLAM)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 14, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEKCAYAAAAB0GKPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xd8VFX6x/HPk4TeS0CQ3qUjEREh\noKAgotgRRVdBURSluZbfqrurrp2qi9hYl1WxgYp0sBB6lRIQpHcFBaLSy/n9cS86YhInkMlNJt/3\n6zWv3Dlz585zdlifOeWeY845RERE/kxM0AGIiEjOoIQhIiJhUcIQEZGwKGGIiEhYlDBERCQsShgi\nIhIWJQwREQmLEoaIiIRFCUNERMISF3QAmal06dKuSpUqQYchIpJjLF68+AfnXHw450ZVwqhSpQqL\nFi0KOgwRkRzDzDaHe666pEREJCxKGCIiEhYlDBERCYsShoiIhEUJQ0REwqKEISIiYVHCEBGRsChh\nAC99vpbk7SlBhyEikq3l+oSx78AR3l2whauHz2bEjPWcOKE9zkVEUpPrE0bxgnmZ1KcV7c4py7OT\nVtPtzfnsTDkYdFgiItlOrk8Y4CWN4Tefy/PXNmTp1n10GDKTyck7gw5LRCRbUcLwmRk3nFeRCfe3\nonKpgtz99hIe+mg5+w8fCzo0EZFsQQnjFFVLF2JMrxbc06Y6HyzeSqeXZrF8276gwxIRCZwSRiry\nxMbwYIc6jL6zOYeOHuea4XMY/tU6jmtAXERyMSWMdDSvVorJfRJpX+8snp+8hpvfmMeOfRoQF5Hc\nSQnjTxQrmIeXb2rCC9c1ZPm2FC4bOpOJKzQgLiK5T8QShpmNNLNdZpYcUtbYzOaZ2VIzW2RmzVJ5\nX2Mzm2tmK81suZl1iVSM4TIzrk+oyMT7W1GldCHueWcJf/1wmQbERSRXiWQL4y2gwyllzwP/dM41\nBh73n5/qAHCrc66e//4hZlY8gnGGrUrpQnx09wX0vqgGHy3ZxuXDZrJ0qwbERSR3iFjCcM4lAXtO\nLQaK+sfFgB2pvO9b59xa/3gHsAsIa7/ZrJAnNoYH2tfmvTubc+TYCa57ZQ7//lID4iIS/bJ6DKMv\n8IKZbQVeBB5J72S/yyovsD4LYsuQ86uVYlKfRNrXP4sXpqyh6+vz2K4BcRGJYlmdMHoB/ZxzFYF+\nwJtpnWhm5YD/Abc7506kc15Pfzxk0e7duzM94PQUK5iHl7s24cXrG7FyewqXDUli/PI/NJpERKKC\nORe5rhQzqwKMd87V95+nAMWdc87MDEhxzhVN5X1Fga+AZ5xzH4b7eQkJCW7RokWZEXqGbf5xP33e\nW8rSrfu49twK/LNzPQrniwskFhGRcJnZYudcQjjnZnULYwfQ2j++GFh76glmlhf4GBiVkWQRtMql\nCvHh3Rdw/8U1+PjrbXQcOpOvt+wNOiwRkUwTyWm1o4G5QG0z22ZmPYA7gYFmtgx4Gujpn5tgZm/4\nb70BSARu86ffLjWzxpGKMzPliY2h/6W1ea/nBRw/4bhuxFxe+nytBsRFJCpEtEsqqwXZJXWqlINH\neeyTZMYt20FC5RIM7tKYiiULBh2WiMjvZOcuqVyjWIE8DL2xMYO7NGLNdz9z2dCZfLR4G9GUoEUk\nd1HCiCAz4+omFZjUtxV1yxflgQ+Xcc87S9i7/0jQoYmIZJgSRhaoUKIgo+9szkMd6jD9m+9pPySJ\npG+zdgqwiMiZUsLIIrExRq821fn4ngspWiAPt45cwD/GreTQ0eNBhyYiEhYljCxW/+xijL+vJbe1\nqMJbczZxxUuzWLkjJeiwRET+lBJGAPLnieUfV9bjv92bkXLwKFf9ezYjZqzX9FsRydaUMALUulY8\nU/om0rZOWZ6dtJqbXp/Htr0Hgg5LRCRVShgBK1EoL690O5cXrmtI8vYULhsyk0++3q7ptyKS7Shh\nZAMnN2ia1CeR2mcVoe/7S7lv9NekHDgadGgiIr9SwshGKpUqyPt3XcBf29dmcvJ3dBiaxJx1PwQd\nlogIoISR7cTGGPdeVIOx97SgQN5YbnpjPk+NX6XptyISOCWMbKphheJMuK8VtzSvzBuzNnLVv2ez\n+rufgg5LRHIxJYxsrEDeWJ68qj7/ue08fvjlMFe+NJs3Zm7ghKbfikgAlDBygIvqlGFK30Ra147n\nqQnf0O3N+exM0XawIpK1lDByiFKF8/HaLU159poGLN26j/aDk/hsmbaDFZGso4SRg5gZNzarxMT7\nW1EtvjD3jf6avu99TcpBTb8VkchTwsiBqpQuxEd3X0DfdjX5bPlOOgxJYtZaTb8VkchSwsih4mJj\n6NuuFh/f04KCeWPp9uZ8/jFuJQePaPqtiESGEkYO17BCcSbc3+rX1W8vf2kmy7buCzosEYlCShhR\n4OTqt+/ccT4HjxznmlfmMHjatxw9fiLo0EQkiihhRJELa5Rmct9ErmxUnqGfr+XaV+awbtcvQYcl\nIlFCCSPKFCuQh8FdGjP85nPZuucAlw+bychZG3Wzn4icMSWMKNWxQTmm9E2kRfVSPDF+FbeMnM+O\nfbrZT0ROnxJGFCtTND8jbzuPZ65pwNdb9tF+SBJjl2zTXhsiclqUMKKcmdG1WSUm9WlF7bJF6P/B\nMu55Zwl79h8JOjQRyWGUMHKJyqUK8f5dF/BQhzpM/+Z7Lh2cxBervw86LBHJQZQwcpHYGKNXm+p8\nem9LShfOS/e3FvHI2OXsP3ws6NBEJAdQwsiF6pYvyqe9L+Su1tV4b+FWLhs6k4Wb9gQdlohkc0oY\nuVS+uFgeuewc3u95AQ7HDa/O5dlJqzl8TEuLiEjqlDByuWZVSzKpTyJdEioyYsZ6Or88m292amc/\nEfkjJQyhcL44nr22IW/cmsAPvxym88uzGTFjPcd1s5+IhFDCkF+1q1uWKX0TubhOGZ6dtJobX5vL\nlh8PBB2WiGQTShjyO6UK5+OVbucy6IZGrN75Mx2GJvH2vM262U9ElDDkj8yMa86twJR+iTStXIJH\nP0nm1pELtLSISC6nhCFpKl+8AKO6N+Opq+qzePNe2g9O4qPFWlpEJLdSwpB0mRndmldmUp9WnFOu\nKA98uIw7Ry1m18+Hgg5NRLJYxBKGmY00s11mlhxS1tjM5pnZUjNbZGbN0njvZDPbZ2bjIxWfZEzl\nUoUY3bM5j15+Dklrd9N+cBITlu8MOiwRyUKRbGG8BXQ4pex54J/OucbA4/7z1LwA3BK50OR0xMYY\nd7SqxsT7W1KpZEHufXcJvd9dwl4tZCiSK0QsYTjnkoBT15twQFH/uBiwI433fg78HKnY5MzUKFOE\nMb1a8MCltZiy8jsuGZzE9FVayFAk2mX1GEZf4AUz2wq8CDySxZ8vmSQuNobeF9f8dSHDO0Yt4oEP\nl/HToaNBhyYiEZLVCaMX0M85VxHoB7x5phc0s57+eMii3bt3n3GAkjF1yxdlXO+W9L6oBmOXbKPD\n4CRmrf0h6LBEJAKyOmH8BRjrH38IpDronRHOudeccwnOuYT4+PgzvZychrxxMTzQvjZjerUgf95Y\nur05n8c+Sday6SJRJqsTxg6gtX98MbA2iz9fIqhJpRJMvL8VPVpW5e35m7ls6EwWbNSy6SLRIpLT\nakcDc4HaZrbNzHoAdwIDzWwZ8DTQ0z83wczeCHnvTLwWSFv/ve0jFadkrvx5YnmsU13eu7M5DkeX\n1+byrwmrOHRUy6aL5HQWTXftJiQkuEWLFgUdhvj2Hz7GM5O+4e15W6geX4hBNzSmUcXiQYclIiHM\nbLFzLiGcc3Wnt0RMoXxxPHVVA0Z1b8aBI8e55pU5DJy6hiPHTgQdmoichj9tYZhZAtAKKA8cBJKB\n6c65bNc5rRZG9pVy8ChPfLaKMUu2cU65ogy6oRHnlCv6528UkYjKlBaGmd1mZkvw7pUoAKwBdgEt\ngWlm9l8zq5QZAUv0K1YgDwNvaMRrtzRl98+HuPLlWfz7y3UcO67WhkhOEZfOa4WAC51zqa5pbWaN\ngZrAlkgEJtHp0npnkVClJI99mswLU9YwdeV3vHh9I2qWLRJ0aCLyJzToLYEZv3yHf7/GcfpeUpOe\nraoRF6thNZGslJEuqfRaGJhZfqATfxzDmOCcW3mmgUru1qlheZpXK8VjnyTz/OQ1TFn5PS9e11Ct\nDZFsKr0xjH8As4ELgPnAq8AHwDHgWTObZmYNsyJIiV6lC+dj+M3n8lLXJmz5cT+XD5vFK1+t19iG\nSDaUZpeUmV3unJuQ5hvNygCVnHPZpg9IXVI52+6fD/PYJ8lMXvkdjSoWZ+D1DalRRq0NkUjKlFlS\n6SUL//Vd2SlZSM4XXyQfr3Q7l2F+a6PjsFmMmLGe4yeiZ5xNJCdLr0uqtJn93czuN7PCZvaKmSWb\n2admViMrg5Tcw8y4slF5pvZrzUW143l20mqufWUO63b9EnRoIrleelNS3gXy4U2dXQBsAK4DxgNv\npPM+kTMWXyQfI7o1ZVjXJmz6cT8dh83kVbU2RAKV3hjGMudcIzMzYLNzrlLIa0v9bVazFY1hRKfd\nPx/m0U9WMGXl9zSuWJwXr29EjTKFgw5LJCpk1lpSxwGcl1FO3RFHU1gky5xsbQy9sbFaGyIBSu8+\njGpmNg6wkGP851UjHplICDOjc+OzuaB6KR79OJlnJq1m8srveOE6tTZEskp6XVKtU33B55ybEZGI\nzoC6pHIH5xzjlu3g7+NWcuDIcR64tBY9WlYjNsaCDk0kx8mUO73TSghmVhG4Ech2CUNyh9DWxt8+\nTubpiauZnPwdL1zfiOrxam2IREpYC/f4U2x7mVkS8BVQNqJRiYShTJH8vHZLU4Z0acz63fvpOHQm\nrydt0NiGSISkdx9GETO71cwm402rrQFUc85Vd849kGURiqTDzLiqydlM659IYq14/jXxG254dS7r\nd+u+DZHMll4LYxfQA/gXUN05NwA4kiVRiWRQaGtj3a5f1NoQiYD0Esb/AfmBV4BHzKx61oQkcnp+\nbW30S6RVTa+1cf0I3SUuklnSW0tqsHPufOBKvKm0nwDlzewhM6uVVQGKZFSZovl5/damDO7SiA0/\nePdtaAVckTP3p4PezrkNzrl/OecaAOcBxYBJEY9M5AyYGVc3qcC0fq1pW6cMz01ezdXD57D6u5+C\nDk0kx0pv0PsPk9qdcyucc//nnKue1jki2Ym3Am5Tht98LjtTDnLFS7MYMv1bjhxTa0Mko9JrYXxp\nZveZWaXQQjPLa2YXm9l/gb9ENjyRzNGxQTmm9mvN5Q3KMWT6Wq58eRYrtqUEHZZIjpJewuiAt57U\naDPbYWarzGwjsBboCgx2zr2VBTGKZIqShfIy5MYmvHFrAnsPHOGq4bN5fvJqDh09HnRoIjlCmkuD\n/O4kszxAaeCgc25fxKM6TVoaRMKVcvAoT0/4hvcXbaV6fCGev64RTSuXCDoskSyXWavVnrxYSaAI\ncBiIMbOSfgIRybGKFcjDc9c1ZFT3Zhw6eoLrRszhyfGrOHhErQ2RtISzNMgSYDfwLV531G5go5kt\nMbOmkQxOJNISa8UzpV8i3c6vzJuzNtJhaBJz1/8YdFgi2VI4CWMy0NE5V9o5Vwq4DPgAuAcYHsng\nRLJC4XxxPHlVfd7r2RyArq/P49FPVvDL4WMBRyaSvYSTMBKcc1NOPnHOTQUSnXPz8LZwFYkKzauV\nYnKfRO5oWZV35m+h/eAkkr7dHXRYItlGOAljj393d2X/8SCw18xi0c57EmUK5I3l0U51GdOrBQXy\nxnLryAU8+NEyUg4eDTo0kcCFkzBuAirgLQ3yKVDJL4sFbohcaCLBObdSCcbf15J7L6rOmCXbuWTQ\nDKat+j7osEQCFda02pxC02olEpK3p/DAh8tY/d3PdG5cnr9fUY+ShfIGHZZIpsjsabXxZvaCmU00\nsy9OPs48TJGcof7ZxRjXuyX92tVi4oqdXDJoBhOW7ww6LJEsF06X1DvAaqAq8E9gE7AwgjGJZDt5\n42Lo064mn93XkrNLFODed5fQ6+3F7Pr5UNChiWSZcBJGKefcm8BR59wM51x3oHmE4xLJluqcVZSx\nvVrw8GV1+Hz1Li4ZlMSYxduIpq5dkbSEkzBOTg/ZaWaXm1kTvEFwkVwpLjaGu1tXZ1KfVtQsU5gB\nHy7j1pEL2LrnQNChiURUOAnjKTMrBgwAHgDeAPqFc3EzG2lmu8wsOaSssZnNM7OlZrbIzJql8d6/\nmNla/6FVcSXbqR5fmA/uuoAnO9djyea9tB+SxH9mb9S2sBK1IjpLyswSgV+AUc65+n7ZVLyVbieZ\nWUfgQedcm1PeVxJYBCQADlgMNHXO7U3v8zRLSoKyfd9BHv14BV+u2U2TSsV5/tqG1CxbJOiwRP5U\nZs+Sqmpmg8xsrJmNO/kI5+LOuSRgz6nFQFH/uBiwI5W3tgemOef2+EliGt5y6yLZ0tnFCzDytvMY\nemNjNvnbwg6dvlYbNUlUiQvjnE+AN4HPyJw7u/sCU8zsRbyE1SKVc84GtoY83+aX/YGZ9QR6AlSq\nVCm1U0SyhJnRufHZtKxRmifGr2Lw9G+ZuGInz13XkMYViwcdnsgZC2cM45Bzbphz7kt/ltQM59yM\nM/jMXkA/51xFvLGQN1M5J7WtX1PtO3POveacS3DOJcTHx59BWCKZo1ThfAy9sQkjb0vgp0NHuWb4\nbJ4cv4oDR7SYoeRs4SSMoWb2dzO7wMzOPfk4g8/8CzDWP/4QSG3QextQMeR5BVLvuhLJti6uU5ap\n/RK56fxKvDlrI+2HJDFr7Q9BhyVy2sJJGA2AO4FngYH+48Uz+MwdQGv/+GK8PTZONQW41MxKmFkJ\n4FK/TCRHKZI/D09d1YAP7rqAPDExdHtzPn/9cBkpB7SYoeQ84YxhXA1Uc84dyejFzWw00AYobWbb\ngL/jJZ+hZhYHHMIffzCzBOBu59wdzrk9ZvYkv91R/oRz7tTBc5Eco1nVkkzs04phn6/l1aQNfLlm\nN092rsdlDcoFHZpI2P50Wq2ZvQ/c55zblTUhnT5Nq5WcYOWOFB4as5zk7T/Rvl5ZnuxcnzJF8wcd\nluRSmTqtFigLrDazKRmdVisif1SvfDE+uedCHr6sDl+t2U3bQTN4f+EWLS8i2V44LYzWqZWf4Uyp\niFALQ3KajT/s5+Exy5m/cQ8tqpfimWsaULlUoaDDklwkIy2MNBOGmZn7k2wSzjlZSQlDcqITJxzv\nLdzKMxO/4eiJEwy4pDa3X1iFuNhwOgBEzkxmdUl9aWb3mdnv7oYzs7xmdrGZ/RdviqyInIGYGOOm\n8ysxrX9rWtaI518Tv+GaV+bwzc6fgg5N5HfSSxgdgOPAaDPbYWarzGwD3jTYrnjrQb2VBTGK5Apn\nFcvP67c25eWbmrBj30GueGkWA6eu4fCx40GHJgKEufigmeUBSgMHnXP7Ih7VaVKXlESLvfuP8OSE\nVYxdsp3q8YV49tqGnFelZNBhSRTK7FlSOOeOOud2ZudkIRJNShTKy6AbGvPf7s04dPQE14+Yy98+\nXsFPh3TDnwRHo2oi2VjrWvFM65/IHS2rMnrBFtoNnMHkZO0nLsFQwhDJ5grmjePRTnX55N4LKV04\nH3e/vYSeoxbxXYr2E5esFc5+GHVTKWsTkWhEJE0NKxTn094X8shldUhau5t2g2bwv7mbOKEd/iSL\nhNPC+MDMHjJPATN7CXgm0oGJyB/liY3hrtbVmdI3kcYVi/PYpyu5/tW5fPv9z0GHJrlAOAnjfLyl\nxufgLQa4A7gwkkGJSPoqlyrE/3o0Y+D1jdiw+xcuHzaTQVPXcOiopuBK5ISTMI4CB4ECQH5go3NO\n+06KBMzMuLZpBab3b02nhuUZ9sU6Og6byfwNPwYdmkSpcBLGQryEcR7QEuhqZh9FNCoRCVupwvkY\n3MWbgnvk2Am6vDaPR8YuJ+WgpuBK5gonYfRwzj3u34vxnXOuM/BppAMTkYxpXSueqf0S6ZlYjfcX\nbqXdoBlMXLFTq+BKpgknYewys0qhDyDbrVQrIt4U3P/reA7jerekTJF83PPOEu4ctZgd+w4GHZpE\ngXCWN18BOMDwxjCqAmucc/UiH17GaGkQkd8cO36C/8zexMBpa4g148EOdejWvDKxMRZ0aJKNZOrS\nIM65Bs65hv7fmkAzYNaZBikikRUXG8OdidWY1q8151Yuwd/HreS6EXNY/Z1WwZXTk+E7vZ1zS/AG\nwEUkB6hYsiCjujdjcJdGbP7xAJ2GzeLFKZqCKxkX92cnmFn/kKcxwLnA7ohFJCKZzsy4ukkFWtcq\nw1MTVvHyl+uYuGInT1/TgObVSgUdnuQQ4bQwioQ88gETgM6RDEpEIqOkvwru/3o049gJx42vzeOh\nj5aTckBTcOXPhbUfRk6hQW+R8B08cpwhn3/LGzM3UqJgXh6/oi5XNCyHmQbFc5PM2tP7M7zZUaly\nzl15euFFjhKGSMYlb0/hkbErWLE9hda14nnqqvpULFkw6LAki2RWwmid3hudc9nuXgwlDJHTc/yE\nY9TcTbw4ZQ3HnaNP21rc0aoqeWK1A0K0y6yEUck5tyVTI4swJQyRM7Mz5SD/GLeSKSu/p3bZIjx9\nTQOaVi4RdFgSQZl1H8YnIRccc8ZRiUi2V65YAV69JYHXbmnKT4eOct2IOfzt4xVal0qA9BNG6MhX\ntUgHIiLZx6X1zmJa/9Z0v9DfGnbQDD5btkPrUuVy6SUMl8axiOQChfPF8Vinuozr3ZKziubnvtFf\nc9t/FrJ1z4GgQ5OApDeGcRzYj9fSKACc/FdigHPOFc2SCDNAYxgikXHqoHjfdrXo0VKD4tEgU8Yw\nnHOxzrmizrkizrk4//jk82yXLEQkcmJjjNsvrMr0Aa1pXSueZyet5oqXZrF4896gQ5MspJ8HIhK2\n0EHxlIPeoPijn2hQPLdQwhCRDDs5KH57i6q8O98bFB+/XIPi0U4JQ0ROS+F8cTx+RV0+vbclZYvm\no/e7X3P7WxoUj2ZKGCJyRhpUKMYn91zI453qsnDjHi4ZPIMRM9Zz9PiJoEOTTKaEISJnLC42hu4t\nqzKtf2sSa/42KL5kiwbFo4kShohkmvLFC/DarQm86g+KX/uKNyj+0yENikeDiCUMMxtpZrvMLDmk\n7H0zW+o/NpnZ0jTe28fMks1spZn1jVSMIhIZ7U8ZFG87cAYTlu/UoHgOF8kWxltAh9AC51wX51xj\n51xjYAww9tQ3mVl94E68vcMbAZ3MrGYE4xSRCDh1UPzed5fQXYPiOVrEEoZzLgnYk9pr5u3QcgMw\nOpWXzwHmOecOOOeOATOAqyMVp4hE1slB8cc61WXBxj20GzSDf3+5jiPHNCie0wQ1htEK+N45tzaV\n15KBRDMrZWYFgY5AxSyNTkQyVVxsDD1aeneKtz2nDC9MWUPHYTOZu/7HoEOTDAgqYXQl9dYFzrlv\ngOeAacBkYBlwLK0LmVlPM1tkZot2794diVhFJJOUK1aA4Tc35T+3ncfhY8fp+vo8+n+wlB9+ORx0\naBKGiO7pbWZVgPHOufohZXHAdqCpc25bGNd4GtjmnBv+Z+dq8UGRnOPgkeP8+8t1vJq0noJ543io\nQx1uPK8iMTHaUzwrZdYGSpHSDlidXrIwszL+30rANaTRGhGRnKtA3lgeaF+bSX1acU65Ivzfxyu4\ndsQcVu34KejQJA2RnFY7GpgL1DazbWbWw3/pRk5JAGZW3swmhhSNMbNVwGfAvc453f0jEqVqlCnC\n6DubM+iGRmz58QBXvDyLJ8ev4pfDafZES0Ai2iWV1dQlJZKzpRw4ynNTVjN6wRbKFsnP36+oS4f6\nZ+FNrJRIyO5dUiIiqSpWMA9PX92AMb1aULJQXnq9o3s3shMlDBHJds6tVIJxvXXvRnajhCEi2dLJ\nezc+H9BG925kE0oYIpKtnVUsv3fvxu26dyNoShgikiNcVLsM0/q15r6La/DZsh20HTiDd+Zv5sSJ\n6Jm4k90pYYhIjpE/TywDLq3NpD6JnFOuCH/7OJlrR8xh5Y6UoEPLFZQwRCTHqVGmMKPvbM7gLv69\nGy/p3o2soIQhIjmSmXF1kwp8MaANXZtVYuTsjbQbOINJK7TvRqQoYYhIjlasYB7+dXUDxobcu3H7\nWwvZ/OP+oEOLOkoYIhIVmoTcu7Fw4x4uGZzEkOnfcujo8aBDixpKGCISNU7eu/HFA21oX+8shkxf\ny6WDk/hy9a6gQ4sKShgiEnXKFs3PS12b8M4d55Mn1rj9rYX0HLWIbXu1xMiZUMIQkah1YY3STOqT\nyIMdajNz7Q9aYuQMKWGISFTLGxfDPW1qMH1Aa1rXiueFKWvoMDSJ2et+CDq0HEcJQ0RyhbOLF+DV\nWxL4z+3ncfyE4+Y35tP73SV8l3Io6NByDCUMEclVLqpdhil9E+nbriZTV31P24Ff8cbMDRw9rm6q\nP6OEISK5Tv48sfRtV4tp/RJpVrUkT034hk7DZrFg456gQ8vWlDBEJNeqXKoQI287j1dvacovh49x\nw6tz6f/BUnb/rJVwU6OEISK5mpnRvt5ZTOufyD1tqvPZsh1cPPArRs3dxHGthPs7ShgiIkDBvHE8\n2KEOk/ok0rBCMR7/dCVXvjyLJVv2Bh1atqGEISISokaZwrzd43xe6tqEH345zDXD5/DwmOXs3X8k\n6NACp4QhInIKM+OKRuX5fEAb7mhZlQ8Xb+OigV8xesGWXL1hkxKGiEgaCueL49FOdZlwf0tqlSnC\nI2NXcM0rc0jenjs3bFLCEBH5E3XOKsr7dzVn4PWN2Lb3AFe+PIvHP00m5eDRoEPLUkoYIiJhMDOu\nbVqBzwe0oVvzyrw9bzNtB37FR4u35ZpuKiUMEZEMKFYgD090rs+43i2pUKIgD3y4jOtfnZsruqmU\nMERETkP9s4sxtlcLnr+2IRt/2M+VL8/isU+S2XcgemdTKWGIiJymmBjjhvMq8uWANtzSvDLvzN/M\nxQNn8F6UzqZSwhAROUPFCubhn53rM/6+VlSPL8TDY1dw9fDZLNu6L+jQMpUShohIJqlbvigf3HUB\ng7s0YkfKIa4aPpuHxyxnT5Tc9KeEISKSicyMq5tU4IsBrelxoX/T34tf8b+5OX9tKiUMEZEIKJI/\nD492qsukPq2oW64oj326kiuLebcBAAANN0lEQVRemsXizTl3CXUlDBGRCKpVtgjv3nk+L9/UhD37\nj3DtKzl3CXUlDBGRCDMzOjUsz+cDWtPr5BLqL37FyFkbOZaDdvpTwhARySKF8sXxUIc6TO6bSONK\nxXli/CouHzaLeRt+DDq0sChhiIhkserxhRnVvRkjunk7/d342jzuH/0136UcCjq0dEUsYZjZSDPb\nZWbJIWXvm9lS/7HJzJam8d5+ZrbSzJLNbLSZ5Y9UnCIiQTAzOtQ/i+n9W3N/25pMXvkdbQd+xasz\n1nPkWPbspopkC+MtoENogXOui3OusXOuMTAGGHvqm8zsbOB+IME5Vx+IBW6MYJwiIoEpkDeW/pfU\nYlq/RC6oXopnJq2mw9AkZq7dHXRofxCxhOGcSwJSnT9mZgbcAIxO4+1xQAEziwMKAjsiEqSISDZR\nuVQh3vjLeYy8LYHjJxy3vLmAXm8vZvu+g0GH9qugxjBaAd8759ae+oJzbjvwIrAF2AmkOOemZnF8\nIiKBuLhOWab0TWTAJbX4cs0u2g78ipe/WMvhY8eDDi2whNGVNFoXZlYC6AxUBcoDhcysW1oXMrOe\nZrbIzBbt3p39mnAiIhmVP08s97WtyfT+rWlTqwwvTv2WSwcn8cXq7wONK8sTht/NdA3wfhqntAM2\nOud2O+eO4o1ztEjres6515xzCc65hPj4+MwPWEQkIBVKFGTELU0Z1b0ZsTFG97cW0f2thWz6YX8g\n8QTRwmgHrHbObUvj9S1AczMr6I91tAW+ybLoRESymcRa8Uzuk8j/dazD/A0/cungJJ6fvJr9h49l\naRyRnFY7GpgL1DazbWbWw3/pRk7pjjKz8mY2EcA5Nx/4CFgCrPBjfC1ScYqI5AR542LomVidLx9o\nQ6eG5Rj+1XraDpzBuGU7cC5rFjW0rPqgrJCQkOAWLVoUdBgiIhG3ePMeHv90JSt3/MT5VUvy1u3N\nKJA3NsPXMbPFzrmEcM6Ny/DVRUQkcE0rl2Rc75a8v3Ary7buO61kkVFKGCIiOVRsjHHT+ZW46fxK\nWfJ5WktKRETCooQhIiJhUcIQEZGwKGGIiEhYlDBERCQsShgiIhIWJQwREQmLEoaIiIQlqpYGMbPd\nwOag48ig0sAPQQeRxVTn3EF1zhkqO+fCWuo7qhJGTmRmi8JdxyVaqM65g+ocfdQlJSIiYVHCEBGR\nsChhBC837vWhOucOqnOU0RiGiIiERS0MEREJixJGhJhZrJl9bWbj/edmZv8ys2/N7Bszuz+kfJiZ\nrTOz5WZ2bsg1/mJma/3HX4KqS7hSqXNbM1tiZkvNbJaZ1fDL85nZ+36d55tZlZBrPOKXrzGz9sHU\nJDxmtsnMVvj1W+SXlTSzaf53Ns3MSvjlUfE9p1HnF8xstV+vj82seMj5qX6fZtbBL1tnZg8HUZdw\npVbnkNceMDNnZqX951HxPafJOadHBB5Af+BdYLz//HZgFBDjPy/j/+0ITAIMaA7M98tLAhv8vyX8\n4xJB1yuDdf4WOMc/vgd4K+R4hH98I/C+f1wXWAbkA6oC64HYoOuVTn03AaVPKXseeNg/fhh4Lpq+\n5zTqfCkQ5x8/F1LnVL9P/7EeqAbk9c+pG3TdMlJnv7wiMAXv3q/S0fQ9p/VQCyMCzKwCcDnwRkhx\nL+AJ59wJAOfcLr+8MzDKeeYBxc2sHNAemOac2+Oc2wtMAzpkWSUyKI06O6Cof1wM2OEfdwb+6x9/\nBLQ1M/PL33POHXbObQTWAc0iHXsmC63bf4GrQspz/PecGufcVOfcMf/pPKCCf5zW99kMWOec2+Cc\nOwK855+b0wwGHsT7d35S1H7PoC6pSBmC9w/pREhZdaCLmS0ys0lmVtMvPxvYGnLeNr8srfLsKrU6\n3wFMNLNtwC3As375r3Xz/0OTApQi59XZAVPNbLGZ9fTLyjrndgL4f8v45dHyPadW51Dd8X5hQxTX\n2cyuBLY755adcm601DlV2tM7k5lZJ2CXc26xmbUJeSkfcMg5l2Bm1wAjgVZ4TddTuXTKs5106twP\n6Oicm29mfwUG4SWRHF9n34XOuR1mVgaYZmar0zk3auvsnEsCMLO/AceAd/xz06pbaj9Uc1Sdgb/h\ndcWdKlq+51SphZH5LgSuNLNNeE3ti83sbbxfFGP8cz4GGvrH2/D6Qk+qgNd1k1Z5dpRanScAjZxz\n8/1z3gda+Me/1s3M4vC6q/aQs+qMc26H/3cX3nfaDPje74LA/3uy6zEavue06ow/iNsJuNn5nfZE\nb51b443JLPP/zVcAlpjZWURJndMU9CBKND+ANvw2APws0D2kfKF/fDm/HyRb4JeXBDbiDZCV8I9L\nBl2ncOuM13r9Aajll/cAxvjH9/L7Qe8P/ON6/H6QdAPZdNAbKAQUCTmeg9cn/QK/H/R+Plq+53Tq\n3AFYBcSfcn6q36f/b2ODX3Zy0Lte0PXLSJ1POWcTvw165/jvOb2HuqSyzrPAO2bWD/gFr2sGYCLe\nzIp1wAG82VQ45/aY2ZPAQv+8J5xze7I25NPnnDtmZncCY8zsBLAXr38b4E3gf2a2Dq9lcaP/npVm\n9gHef3yOAfc6545nffRhKQt87I3VEwe865ybbGYLgQ/MrAewBbjePz8avue06rwOLylM81+b55y7\nO73v08x6480wigVGOudWZn11wpJqndM5Pxq+5zTpTm8REQmLxjBERCQsShgiIhIWJQwREQmLEoaI\niIRFCUNERMKihCGBMLPBZtY35PkUM3sj5PlAM+ufyZ/5S2Zez79mYzPrGPL8H2b2QBjvMzP7wsyK\nhpRd7a98WicCcVYxs5sy+7oh1+9tZrdH6vqSPShhSFDm4N/5bWYxQGm8G71OagHMDiCujGqMN+8+\nozoCy5xzP4WUdQVm4d+XksmqAKkmDP9u+zM1Erg/E64j2ZgShgRlNr8tFVIPSAZ+NrMSZpYPOAf4\n2swKm9nn5u2rscLMOgOY2XNmds/Ji/m/7Af4x381s4X+fgT/TO3DUzvH/xX+jZm9bmYrzWyqmRXw\nXzvPP3euefs/JJtZXuAJvEUll5pZF//ydc3sKzPbYP6+J6m4Gfg0JJ7CeEus9CAkYZhZG/9aH5m3\n58Q7/sq+mFlHv2yWeXswnNyHpLUfz1Lz9icpgnfjaCu/rJ+Z3WZmH5rZZ3gL61lIvVacrIv/+TPM\n7APz9nJ51sxuNrMF/nnVAZxzB4BNZpbTVheWjAj6VnM9cu8Db0mFSsBdwN3Ak3i/vC8Ekvxz4oCi\n/nFpvDtoDWgCzAi51ir/Wpfi7atseD+IxgOJ/jm/+H9TPQfvV/gxoLF/3gdAN/84GWjhHz8LJPvH\ntwEvh8TxD7zWUz4/3h+BPKnUfTP+khP+827Am/7xHOBc/7gN3mq+FfxY5wItgfx4q59W9c8bzW/L\n0HyGt2AeQGH/f8M2J18PiXsb/vIUwLV4S27H4t3dvAUo579vn3+cD9gO/NN/Tx9gSMg1/wYMCPrf\nlR6Re6iFIUE62cpogfcfwrkhz+f45xjwtJktB6bjLQld1jn3NVDGzMqbWSNgr3NuC14yuBT4GlgC\n1AFq8nvpnbPRObfUP14MVDFvB7kizrmTMb37J/Wa4Lw9IH7AW3ywbCrnlHTO/RzyvCvewo34f7uG\nvLbAObfNeXupLMVLbHWADc7bZwK8hHHSbGCQ37op7n7bq+JU09xvy1O0BEY75447574HZgDn+a8t\ndM7tdM4dxtv4aKpfvsKP5aRdQPk0PkuigNaSkiCdHMdogPcLfiswAPgJr08cvK6beKCpc+6ovzpo\nfv+1j4DrgLP47T+2BjzjnHs1nc9N9Rzztoo9HFJ0HChA6ktTp+fUa6T2/7NjZhbjnDthZqWAi4H6\nZubwfuU7M3swneulGZNz7lnzVgvuCMwzs3ZpnLo/5Di9OoZ+/omQ5yf4fd3yAwfTuY7kcGphSJBm\n4y2Jvcf/ZbsHKA5cgNfaAG/p811+srgIqBzy/vfw+vuvw0se4C1o190fE8DMzjZvH4NQ4ZzzK+ft\nkPazmTX3i0IHpX8GimSk0r41eFuU4sc/yjlX2TlXxTlXEW8105bpvH81UM1+2w/95PgJZlbdObfC\nOfccsAivNfJncSbhjcXEmlk8XhfdggzWqRZe4pcopYQhQVqB188/75SyFL87B7zNeBLMbBFea+PX\nTYqct8JpEbydz07ucjcVr8torpmtwEskv/sPZTjnpKIH8JqZzcX7NZ7il3+JN8gdOugdjgl44wPg\ndT99fMrrY0hjVpNfh4N4e6NPNrNZwPchMfX1B6+X4f3inwQsx2vVLDNvxeRTfeyfswz4AnjQOfdd\nBuoD3tjT9Ay+R3IQrVYrEgYzK+yc+8U/fhgo55zrcwbXK4fXqrjkTGPyZ039G1jrnBt8utc7E2bW\nBOjvnLsliM+XrKEWhkh4LvdbEcl4W+s+dSYX81tEr1vIjXun4U4zWwqsxOu6S2/cJtJKA48F+PmS\nBdTCEBGRsKiFISIiYVHCEBGRsChhiIhIWJQwREQkLEoYIiISFiUMEREJy/8DHIHguL65kYsAAAAA\nSUVORK5CYII=\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe3fc06518>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"# But not in ABMAG (or FNU) space.\n", | |
"input_spectrum.plot(bigbox.waveset, flux_unit=u.ABmag)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 17, | |
"metadata": {}, | |
"outputs": [], | |
"source": [ | |
"obs = Observation(input_spectrum, smallbox, binset=smallbox.waveset)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 18, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"text/latex": [ | |
"$17.996364 \\; \\mathrm{mag}$$\\mathrm{\\left( \\mathrm{AB} \\right)}$" | |
], | |
"text/plain": [ | |
"<Magnitude 17.99636401 mag(AB)>" | |
] | |
}, | |
"execution_count": 18, | |
"metadata": {}, | |
"output_type": "execute_result" | |
} | |
], | |
"source": [ | |
"obs.effstim('ABmag')" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 24, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAERCAYAAACU1LsdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAGo9JREFUeJzt3XmYJVWd5vHvC8Umi6CUqEhZyqho\n28pSrogybjiura2jtqIoDm27oeI2znS7TY862Kg9ts3Q4oIiPijQreBGK4sgIlVsBZSgIzagKCCN\nFLbaAr/+IyKDJM28N7Mq42Zm3e/nefKpuBHnRpxDkve9cSLOiVQVkiQBbLbQFZAkLR6GgiSpYyhI\nkjqGgiSpYyhIkjqGgiSpsyRDIcknk1yX5JJ52t/Xk9yU5OQp65Pkr5NckWRdkjfMx/EkabFakqEA\nfBp42jzu73DgwGnWHwTsBuxRVQ8GvjCPx5SkRWdJhkJVnQncOHldkt3bb/xrknwnyR5z2N+3gPXT\nbPoL4L1VdXtb7rqNqbckLXZLMhRmcBTw+qraB3gL8PF52OfuwAuTrE7ytSQPmId9StKitWyhKzAf\nkmwHPBb4YpKJ1Vu1254HvHeat/20qg4YsuutgN9W1ap2P58E9pufWkvS4rNJhALNGc9NVbXn1A1V\ndSJw4gbu9xrghHb5JOBTG7gfSVoSNonuo6q6GbgyyQugu2vo4fOw638EntguPwG4Yh72KUmLVvqa\nJTXJ1sCZNF0wy4AvVdW7ppTZCjgG2Af4JfDCqvrJLPZ9HLA/sDPwC+BdwLeBvwfuBWwBfKGqpus2\nmm5/3wH2ALZr63FwVX0jyY7AscAK4Bbg1VV10Wz2KUlLUZ+hEGDbqrolyRbAWcChVfW9SWVeAzys\nql6d5EXAc6vqhb1USJI0VG/dR9W4pX25RfszNYGeA3ymXf4S8KRMulIsSRqtXi80J9kcWAP8J+Dv\nqurcKUV2Ba4GqKpbk/wKuDtww5T9HAIcArDtttvus8cesx6CIEkC1qxZc0NVLR9WrtdQqKrbgD3b\nvvmTkjy0qiZPTTHdWcEf9GdV1VE04xBYtWpVrV69upf6StKmKsm/zKbcSO4+qqqbgNP5w6kprqGZ\nRoIky4C7MmWksiRpdHoLhSTL2zMEkmwDPBn4wZRiXwZe3i4/H/h2+dBoSVowfXYf3Qv4THtdYTPg\n+Ko6Ocl7gdVV9WXgaOCzSX5Ec4bwoh7rI0kaordQqKqLgb2mWf9Xk5Z/C7ygrzpIkuZmkxjRLEma\nH4aCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaC\nJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKlj\nKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOoaCJKljKEiSOr2FQpLdkpyWZF2SS5McOk2Z\nnZKclOTiJN9P8tC+6iNJGq7PM4VbgcOq6sHAo4HXJnnIlDLvBC6sqocBLwM+2mN9JElD9BYKVXVt\nVZ3fLq8H1gG7Tin2EOBbbZkfACuT7NJXnSRJg43kmkKSlcBewLlTNl0EPK8t80jgvsB9pnn/IUlW\nJ1l9/fXX91tZSRpjvYdCku2AE4A3VtXNUzZ/ANgpyYXA64ELaLqd7qSqjqqqVVW1avny5X1XWZLG\n1rI+d55kC5pAOLaqTpy6vQ2JV7RlA1zZ/kiSFkCfdx8FOBpYV1VHzFBmxyRbti9fBZw5zdmEJGlE\n+jxT2Bc4EFjbdg9Bc7fRCoCqOhJ4MHBMktuAy4CDe6yPJGmI3kKhqs4CMqTMOcAD+qqDJGluHNEs\nSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoY\nCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKk\njqEgSeoYCpKkjqEgSeoYCpKkjqEgSeoYCpKkzrJhBZJsDTwT2A+4N/Ab4BLglKq6tN/qSZJGaWAo\nJHk38CzgdOBc4Dpga+CBwAfawDisqi7ut5qSpFEYdqZwXlW9e4ZtRyS5B7Biuo1JdgOOAe4J3A4c\nVVUfnVLmrsDn2n0sAz5UVZ+affUlSfNpYChU1SlDtl9Hc/YwnVtpziLOT7I9sCbJqVV12aQyrwUu\nq6pnJVkOXJ7k2Kr69zm0QZI0T4Z1H3150PaqevaAbdcC17bL65OsA3YFJodCAdsnCbAdcCNNmEiS\nFsCw7qPHAFcDx9FcU8iGHCTJSmCvdh+TfQz4MvAzYHvghVV1+zTvPwQ4BGDFiml7qyRJ82DYLan3\nBN4JPBT4KPAU4IaqOqOqzpjNAZJsB5wAvLGqbp6y+QDgQpq7mvYEPpZkh6n7qKqjqmpVVa1avnz5\nbA4rSdoAA0Ohqm6rqq9X1cuBRwM/Ak5P8vrZ7DzJFjSBcGxVnThNkVcAJ1bjR8CVwB5zaoEkad7M\nZpzCVsAzgBcDK4G/Bab7gJ/6vgBHA+uq6ogZil0FPAn4TpJdgAcBP55VzSVJ827YhebP0HQdfQ14\nT1VdMod97wscCKxNcmG77p20t7BW1ZHA+4BPJ1lLc73i7VV1w9yaIEmaL8POFA4Efk0zWO0NzZd/\noPkAr6r6g/7/CVV1FkMuTFfVz4Cnzrq2kqReDRun4NxIkjRG5vyhn2TbJC9JMnBgmyRp6ZlVKCTZ\nMsmfJDmeZkDak4Eje62ZJGnkhl1ofgrNXUcHAKcBnwUeWVWvGEHdJEkjNuxC8zeA7wCPq6orAZJ8\ndPBbJElL1bBQ2Ad4EfDPSX4MfAHYvPdaSZIWxLARzRdU1duranfg3TTzF22Z5GvtfESSpE3IrO8+\nqqqzq+p1NDOdfgR4fm+1kiQtiDnfklpVt1fVN3COIkna5Dg4TZLUMRQkSZ1h4xT+L83T0f5gE7Bj\nLzWSJC2YYbekrt7AbYvS5T9fz9k/chJWSUvPkx+8CyvufpfejzMsFI6tqk3mmcnv/9o6Tr/8+oWu\nhiTN2Q+vW8/7n/ew3o8zLBS+D+wNTVdSVc3qiWuL1a23FQ+/z1055pWPWuiqSNKsHfCRM/n9bdP1\n5M+/YaEw+XkI+/ZZkVFZtvlm3PUuWyx0NSRp1jbfbOCjaebVsLuPRhNNkqRFYdiZwh5JLqY5Y9i9\nXYY7nrzWfwfXPCozTpIGGhYKDx5JLUZodCdhkrT0DAuFq6pq4NfrJBlWRpK0NAy7pnBaktcnWTF5\nZfskticm+Qzw8v6qJ0kapWFnCk8DXgkcl+R+wE3A1jTPVPgm8OGqurDfKs4fz2ckabCBoVBVvwU+\nDnw8yRbAzsBvquqmUVSuD/GigiTNaNiZQqeqfg9c22NdJEkLbKxmSbX7SNJSNarPr7EKBUnSYLMK\nhSQPmWbd/vNemxGIIxUkaUazPVM4Psnb09imfc7C+/usmCRp9GYbCo8CdgO+C5wH/IwlOEGe01xI\n0mCzDYXfA78BtqEZp3BlVd3eW60kSQtitqFwHk0oPAJ4HPDiJF/qrVZ98pKCJM1otuMUDq6qicdv\n/hx4TpIDe6qTJGmBzDYUrps6/xFwxnxXpm+OU5CkwWYbCqfQPHAnNNcU7gdcDvxRT/WSJC2AWYVC\nVf3x5NdJ9gb+vJca9cxLCpI0sw0a0VxV59NcdJ5Rkt2SnJZkXZJLkxw6TZm3Jrmw/bkkyW1J7rYh\ndZKkTdmobqmf1ZlCkjdPerkZsDdw/ZC33QocVlXnJ9keWJPk1Kq6bKJAVR0OHN4e41nAm6rqxrk0\nYC4m+r8kSdOb7TWF7Sct30pzjeGEQW+oqmtpZ1WtqvVJ1gG7ApfN8JYXA8fNsj6SNDZGOeX/bK8p\nvGdjDpJkJbAXcO4M2+9C80Cf182w/RDgEIAVK6beBDXXumzU2yVpkzYwFJJ8BWbuyKqqZw87QJLt\naM4q3lhVN89Q7FnA2TN1HVXVUcBRAKtWrfLGUknqybAzhQ9tzM7bp7WdABxbVScOKPoiRtF15EUF\nSRpoWChcWVVXbciOkwQ4GlhXVUcMKHdX4AnASzfkOHOul6kgSTMadkvqP04sJBl4YXka+wIHAk+c\ndNvp05O8OsmrJ5V7LvDNqvr1HPcvSZpnw84UJn+tvv9cdlxVZzGLzpqq+jTw6bnsW5LUj2FnCjXD\n8pLk8xQkabBhZwoPT3IzzTf+bdpl2tdVVTv0WrseeEuqJM1sYChU1eajqogkaYARdXRs0NxHS5VT\nZ0vSYGMVCpK0FI2y23vsQsFrCpI0s7ELBUnSzMYqFLykIEmDjVUoSJIGG7tQcO4jSZrZ2IWCJGlm\nYxUK5UAFSRporEJBkjTY2IWC4xQkLUWj6ucYu1CQJM1srELBKwqSlqJR3jU5VqEgSRrMUJAkdQwF\nSVJnrELBYQqSNNhYhQJAvCdVkmY0dqEgSZqZoSBJ6oxVKHhJQZIGG6tQAJw4W5IGGLtQkKSlaFSz\nPBsKkrTIjfKmyfEKBQcqSNJA4xUKOHW2JA0ydqEgSZrZWIWCnUeSNNhYhYIkabCxCwUvKUjSzMYu\nFCRJM+stFJLsluS0JOuSXJrk0BnK7Z/kwrbMGX3VB7wjVZKGWdbjvm8FDquq85NsD6xJcmpVXTZR\nIMmOwMeBp1XVVUnu0WN9JElD9HamUFXXVtX57fJ6YB2w65RifwacWFVXteWu66s+E3yegqSlaFQd\nHSO5ppBkJbAXcO6UTQ8EdkpyepI1SV42w/sPSbI6yerrr7++38pK0iIzyq+yvYdCku2AE4A3VtXN\nUzYvA/YBngEcAPxlkgdO3UdVHVVVq6pq1fLlyze4LuVIBUkaqM9rCiTZgiYQjq2qE6cpcg1wQ1X9\nGvh1kjOBhwNX9FanvnYsSZuAPu8+CnA0sK6qjpih2D8B+yVZluQuwKNorj1IkhZAn2cK+wIHAmuT\nXNiueyewAqCqjqyqdUm+DlwM3A58oqou6bFOkqQBeguFqjqLWfTWVNXhwOF91ePOxxrFUSRp6Rq7\nEc3ekSpJMxu7UJAkzcxQkCR1xioUvKYgSYONVSg0vKggaekZ1ZfaMQwFSVpaRjlnm6EgSeqMVSh4\nSUGSBhurUADHKUjSIGMXCpKkmRkKkqTOWIVCOVBBkgYaq1AARylI0iBjFwqSpJkZCpKkjqEgSUvA\nqK6Ijl0oOE5B0lIzyo+tsQsFSdLMxioUvCNVkgYbq1AAiDelStKMxi4UJEkzMxQkSZ2xCoVy8mxJ\nGmisQgG8JVWSBhm7UJAkzcxQkCR1xioUHKcgaaka1dT/YxUK4DUFSUvQCD+3xi4UJEkzMxQkSZ2x\nCgUvKUjSYGMVCuDcR5I0yNiFgiRpZoaCJKkzVqEwqvt8JWmp6i0UkuyW5LQk65JcmuTQacrsn+RX\nSS5sf/6qr/rccdDejyBJS9ayHvd9K3BYVZ2fZHtgTZJTq+qyKeW+U1XP7LEekqRZ6u1Moaqurarz\n2+X1wDpg176OJ0mbslF1fmcU/exJVgJnAg+tqpsnrd8fOAG4BvgZ8JaqunSa9x8CHNK+fBBweb81\n7sXOwA0LXYkRG7c2j1t7wTYvJfetquXDCvUeCkm2A84A/rqqTpyybQfg9qq6JcnTgY9W1QN6rdAC\nSbK6qlYtdD1GadzaPG7tBdu8Ker17qMkW9CcCRw7NRAAqurmqrqlXf4qsEWSnfuskyRpZn3efRTg\naGBdVR0xQ5l7tuVI8si2Pr/sq06SpMH6vPtoX+BAYG2SC9t17wRWAFTVkcDzgb9IcivwG+BFtekO\nJjhqoSuwAMatzePWXrDNm5yRXGiWJC0NYzWiWZI0mKEgSeoYChspyeZJLkhycvv6iUnOT3JJks8k\nWdau3ynJSUkuTvL9JA+dtI8dk3wpyQ/aaUEes1DtmY15avOb2ulPLklyXJKtF6o9wyT5SZK17VQs\nq9t1d0tyapIftv/u1K5Pkr9N8qO23XtP2s/L2/I/TPLyhWrPMPPR3iR7Jjmn/R1fnOSFC9mmYebr\nd9xu3yHJT5N8bCHastGqyp+N+AHeDHweOJkmZK8GHthuey9wcLt8OPCudnkP4FuT9vEZ4FXt8pbA\njgvdrj7bTDOy/Upgm/b18cBBC92uAe39CbDzlHX/B3hHu/wO4IPt8tOBr9HMsvVo4Nx2/d2AH7f/\n7tQu77TQbeuxvQ8EHtAu3xu4djH/fz0fbZ70vo+2fx8fW+h2bciPZwobIcl9gGcAn2hX3R34XVVd\n0b4+FfjTdvkhwLcAquoHwMoku7QD+B5Pc/suVfXvVXXTiJowZ/PR5nbbMmCb9qziLjQj2peS59CE\nOe2/fzJp/THV+B6wY5J7AQcAp1bVjVX1rzT/nZ426kpvhDm1t6quqKofAlTVz4DrgKGjaReZuf6O\nSbIPsAvwzVFXdr4YChvnI8DbgNvb1zfQDMCbGO34fGC3dvki4HnQjcm4L3Af4P7A9cCn2i6ZTyTZ\ndkT13xAb3eaq+inwIeAqmm+Qv6qqxfxHVMA3k6xJM+UKwC5VdS0083wB92jX70pz5jThmnbdTOsX\no/lob6f93W8J/P9ea71xNrrNSTYD/gZ464jq3AtDYQMleSZwXVWtmVhXzbnji4APJ/k+sJ5mtliA\nDwA7pRmz8XrggnbbMmBv4O+rai/g1zSnqovOfLW57Zt9DnA/mq6FbZO8dHQtmbN9q2pv4L8Ar03y\n+AFlp5ucvQasX4zmo73NxuYb9GeBV1TV7dOUXSzmo82vAb5aVVdPs33J6HPw2qZuX+DZaeZs2hrY\nIcnnquqlwH4ASZ5K07dKNRMBvqJdH5o+9Stpuk6uqapz2/1+iUUaCsxfmw8Arqyq69ttJwKPBT43\n2ubMTtv9QVVdl+Qk4JHAL9pukmvbD77r2uLXcMeZEjRngz9r1+8/Zf3pPVd9g8xTeyfmNjsF+J9t\nN8uiNU9tfgywX5LXANsBWya5paoW69/ztDxT2EBV9d+r6j5VtZLmm/K3q+qlSe4BkGQr4O3Ake3r\nHZNs2b79VcCZ1cz99HPg6iQParc9CZj6zIlFYb7aTNNt9Ogkd2nD4kk0U6svOkm2TfM8ENpuvacC\nlwBfBibuIHo58E/t8peBl7V3qDyapmvsWuAbwFPT3JG1U7ufb4ywKbMyX+1tf+8n0fS9f3GkjZij\n+WpzVb2kqla0fx9voWn7kgoE8EyhD29tu1k2o+kS+na7/sHAMUluo/nQP3jSe14PHNv+If2Y9tv1\nEjKnNlfVuUm+BJxP09V0AYt36oBdgJOa7GIZ8Pmq+nqS84DjkxxME3IvaMt/lebulB8B/0b7u6yq\nG5O8DzivLffeqrpxdM2YtXlpL/BfaW6guHuSg9p1B1XVxJQ3i8l8tXmT4DQXkqSO3UeSpI6hIEnq\nGAqSpI6hIEnqGAqSpI6hoF4l+XCSN056/Y0kn5j0+m+SvHmej3nLfO6v3eee7aC9idfvTvKWWbwv\nSb7dDuSaWPfcJJVkjx7quTLJn833fift/3VJNqlbMHVnhoL69l2a0cq0c8PsDPzRpO2PBc5egHrN\n1Z4096bP1dOBi9pBexNeDJxFMwBwvq0Epg2FtFOab6RPAm+Yh/1okTIU1LezaUOBJgwuAda3I3u3\nohngdkGS7ZJ8K81zGdYmeQ5Akg+20wbQvn53ksPa5bcmOS/NnPbvme7g05Vpv02vS/IPaeb7/2aS\nbdptj2jLnpPk8DTPe9iSZkrwF6aZb3/i2QAPSXJ6kh8nmemD8iXcMRKWJNvRTBdyMJNCIcn+7b4m\nnqtxbDvamyRPb9edlWYe/4nnWDyhrc+FaSZT3J5mvqn92nVvSnJQki8m+QrNhG+Z1K61E21pj39G\nkuOTXJHkA0lekuY5GGuT7A5QVf8G/CTNJHfaFC303N3+bPo/NHPVrwD+HHg18D6ab9D70kx9Ac1I\n0h3a5Z1pRosG2As4Y9K+Lmv39VSaUdCh+XJzMvD4tswt7b/TlqH5Nn0rsGdb7njgpe3yJcBj2+UP\nAJe0ywcxaX584N00Z0FbtfX9JbDFNG3/F2D7Sa9fChzdLn8X2Ltd3h/4Fc08OpsB5wCPo5lj6mrg\nfm2544CT2+Wv0EzkBs1cO8va/Zw86XgH0czVc7f29Z/STNu9Oc1I3quAe7Xvu6ld3gr4KfCe9j2H\nAh+ZtM//ARy20P9f+dPPj2cKGoWJs4XH0nzYnTPp9XfbMgH+d5KLgX+mmZ54l6q6ALhHknsneTjw\nr1V1Fc0H/lNppsg4n+YhPg+YctxBZa6sO6ZcWEPzrIcdaT7AJ+r0+SHtOqWqfldVN9BMlrbLNGXu\nVlXrJ71+MfCFdvkL7esJ36+qa6qZTfRCmvDaA/hxVV3ZljluUvmzgSPas5Qdq+pWpndq3TGlxuOA\n46rqtqr6BXAG8Ih223nVzOHzO5ppriemM1/b1mXCdTSz22oT5NxHGoWJ6wp/TPNN/GrgMOBmmj5q\naLpZlgP7VNXvk/yE5lsyNDPHPh+4J3d8oAZ4f1X9vwHHnbZMkpXA7yatug3YhumnRB5k6j6m+3u6\nNclmVXV7krsDTwQemqRovq1XkrcN2N+MdaqqDyQ5heas63tJnjxD0V9PWh7UxsnHv33S69u5c9u2\nBn4zYD9awjxT0CicDTwTuLH9hnojsCPNVMPntGXuSvOsht8n+c80D+SZ8AWa/vfn0wQENDOMvrLt\noyfJrmlna51kNmU61TwRbX2amS/hzheC1wPbz6XRrctpHqREW/9jquq+VbWyqnajmUr8cQPe/wPg\n/m2QAXTPOk6ye1WtraoPAqtpziqG1fNMmmsjmydZTtOd9v05tumBNOGuTZChoFFYS9Pv/r0p637V\ndr0AHAusSvPQ9JfQfBgCUFWX0nzQ/bTueBLWN2m6d85JspYmLO70YTibMtM4GDgqyTk036p/1a4/\njebC8uQLzbNxCnc8R+HFNNNJT3YCM9wt1LbhNzQPb/l6krOAX0yq0xvbC8YX0Xxz/xpwMc3ZyUVJ\n3jTNLk9qy1wEfBt4WzXTt8/FvjRdfNoEOUuqNEmS7arqlnb5HcC9qurQjdjfvWjODp6ysXVq70b6\nO+CHVfXhDd3fxkiyF/DmqjpwIY6v/nmmIN3ZM9qzgUtonib3vzZmZ+2ZzT9k0uC1DfDf0jzS9FKa\nbrZB11H6tjPwlwt4fPXMMwVJUsczBUlSx1CQJHUMBUlSx1CQJHUMBUlS5z8AHOlMb/MZ2IgAAAAA\nSUVORK5CYII=\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe3f7cf208>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"obs.plot(flux_unit=units.FLAM, bottom=2.5e-16, top=3e-16)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": 28, | |
"metadata": {}, | |
"outputs": [ | |
{ | |
"data": { | |
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAERCAYAAABowZDXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3XmYXVWd7vHvWxkAk3AZEhBDyiBi\nRxwgWIJtFNFWBOzHYKvtgBEUbuynHYiiotiDQj8KrY22t1VMA7eljXDRJDYCKmkE6chgBgMhKSaB\nJiHRBIIktIgk9bt/rFVk5+ScOruqzj418H6e5zy1h7X3WSuVWr+919prbUUEZmZmzXQMdQbMzGxk\ncMAwM7NSHDDMzKwUBwwzMyvFAcPMzEpxwDAzs1JGXcCQdKmkTZLubMG5jpR0i6Q1ku6Q9K7Cvv+S\ntCp/Nkj64WC/z8xsONNoG4ch6VjgCeCyiHjpIM/1IiAi4l5JzwNWAC+OiN/VpFsI/EdEXDaY7zMz\nG85G3R1GRNwEbCluk3SopJ9IWpHvDGaUPNc9EXFvXt4AbAKm1Jx7EvAGwHcYZjaqjR3qDLTJfOCv\n8p3CMcA3SZV8aZKOBsYDv67Z9Tbg+ojY2pKcmpkNU6M+YEiaCLwa+L6k3s175H1/AZxb57CHI+LN\nhXMcBPw7cGpE9NSkfQ9wcavzbWY23Iy6PgwASdOBqyPipZL2Bu6OiIMGeK69gRuBL0XE92v27Q/c\nA0yNiD8MKtNmZsPcqOvDqJWbih6Q9E4AJUeUOVbSeGAxqQP9+3WSvJMUmBwszGzUqyxgSJom6QZJ\n3fmx1DP7SPtKSTskvaOw7VRJ9+bPqf343suBW4A/kbRe0unAKcDpkm4H1gCzS57uL4FjgdMKj9Ae\nWdj/buDysnkzMxvJKmuSyu3+B0XEyvwk0Qrg5IhYW5NuDLAE+ANwaUT8QNJ+wHKgC4h87Csi4rFK\nMmtmZk1VdocRERsjYmVe3gZ0A1PrJP0osJD0yGqvNwNLImJLDhJLgBOqyquZmTXXlqekcif0TOC2\nmu1TSY+lvgF4ZWHXVGBdYX099YMNkuYCcwEmTJjwihkzSg2xMDMzYMWKFY9ExJTmKdsQMPJjrQuB\neXXGKnwNODsidhQeeQUQu6vbdhYR80njLOjq6orly5cPPtNmZs8Skv67bNpKA4akcaRgsSAiFtVJ\n0gVckYPFZOAkSdtJdxTHFdIdTHq01czMhkhlAUMpClwCdEfEhfXSRMQhhfT/RnpE9Ye50/uLkvbN\nu48HPltVXs3MrLkq7zBmAXOA1ZJW5W3nAJ0AEXFRowMjYouk84BledO5EbGlUXozM6teZQEjIpZS\nvy+iUfrTatYvBS5tcbbMzGyARv1IbzMzaw0HDDMzK8UBw8zMSnHAMDOzUhwwzMysFAcMMzMrxQHD\nzMxKccAwM7NSHDDMzKwUBwwzMyvFAcPMzEpxwDAzs1IcMMzMrBQHDDMzK8UBw8zMSnHAMDOzUhww\nzMysFAcMMzMrpbKAIWmapBskdUtaI+nMOmlmS7pD0ipJyyW9prBvR96+StJVVeXTzMzKqeyd3sB2\n4KyIWClpErBC0pKIWFtIcz1wVUSEpJcDVwIz8r4nI+LICvNnZmb9UNkdRkRsjIiVeXkb0A1MrUnz\nREREXp0ABGZmNiy1pQ9D0nRgJnBbnX1vk3QXcA3wwcKuPXMz1a2STm5HPs3MrLHKA4akicBCYF5E\nbK3dHxGLI2IGcDJwXmFXZ0R0Ae8Fvibp0Abnn5sDy/LNmzdXUAIzM4OKA4akcaRgsSAiFvWVNiJu\nAg6VNDmvb8g/7wduJN2h1DtufkR0RUTXlClTWpl9MzMrqPIpKQGXAN0RcWGDNC/M6ZB0FDAeeFTS\nvpL2yNsnA7OAtfXOYWZm7VHlU1KzgDnAakmr8rZzgE6AiLgIeDvwfklPA08C78pPTL0Y+LakHlJQ\nO7/m6SozM2uzygJGRCwF1CTNBcAFdbbfDLysoqyZmdkAeKS3mZmV4oBhZmalOGCYmVkpDhhmZlaK\nA4aZmZXigGFmZqU4YJiZWSkOGGZmVooDhpmZleKAYWZmpThgmJlZKQ4YZmZWigOGmZmV4oBhZmal\nOGCYmVkpDhhmZlaKA4aZmZXigGFmZqU4YJiZWSmVvdNb0jTgMuC5QA8wPyL+uSbNbOC8vH87MC+/\nCxxJpwJ/k5P+Q0R8p6q8mpkBRATbe4IdPcHTO3rY0dN4ffuOYHtPzy7rO3qCp3t62LGjkK6n55l9\naVsPT9esb3/mfI3Xi9/Ze+z2HcELD5jI59/6krb8+1QWMEgB4KyIWClpErBC0pKIWFtIcz1wVUSE\npJcDVwIzJO0H/D3QBUQ+9qqIeKzC/JpZQUTUrSCLlVVxvW6l2hNs39HT5/qOvL5rxdvD0/XWm1Sq\nOyvmZvmsXxH3xND9e4/pEGM6xNjCz7FjOvpc/+3Wp1h63yMjP2BExEZgY17eJqkbmAqsLaR5onDI\nBFJwAHgzsCQitgBIWgKcAFxeVX7N+hIR9ASNK6JCRVVbce16dVniyrVepZrPU3vVWpum9sp1tyvZ\nmnP0VsC93/nM9+c0Q6VDMLajY2cFOkaM7ShUlmPSz3G9aWrW9xg3dmcl29GRjy9fEY/p6GDcmF3X\nxxa+pzdv42rWx44pfGdtvgp5qF0fI9HRoX7/O1245B6+fv29FfwG6qvyDuMZkqYDM4Hb6ux7G/Al\n4ADgLXnzVGBdIdn6vK3euecCcwE6OztblWVroqenpuJ55qqt8Xrdq706t/XlbvPrX7nu6OmrUq1f\nye5WiTYoy1CqrZh2q6ieqZxqK0IxftyYXdZ3qYjrVZZjChVx7XqhEh3b0bHL+pi8vvNcu6/3lffi\n+kAqT6te5QFD0kRgIal/Ymvt/ohYDCyWdCypP+ONQL3/LXX/YiNiPjAfoKura0j+qovtnrWV4G7r\nA6hUt9fePte9Aq25Ut2RK9ld8rD7FfHTTa+Q61eqQ3nrvuuV4M4Kr/GV6M79e4zr4DmF/eNqKtny\nFXGzK1Qxbkxtxdz8SrReJdshkFyB2tCrNGBIGkcKFgsiYlFfaSPiJkmHSppMuqM4rrD7YODGqvL5\n0KO/53M/XM3v/7ijeUVcqKx714fy1r1Mu2f9yitVTHuMG9v0Frz26rG4vvNqU4wZ06wiHkglu3sa\nV55mQ6PKp6QEXAJ0R8SFDdK8EPh17vQ+ChgPPAr8FPiipH1z0uOBz1aV15UPPcZ/3fsIr3j+vuw7\nYXzpds96lWyZds/dbsFLtnvW5mGg7Z5mZgNR5R3GLGAOsFrSqrztHKATICIuAt4OvF/S08CTwLsi\nIoAtks4DluXjzu3tAK/SV955BIdMnlD115iZjUhVPiW1lPp9EcU0FwAXNNh3KXBpBVkzM7MB8Ehv\nMzMrxQEDiPoPYJmZWYEDRoG7j83MGnPAMDOzUhwwzMysFAcMINyFYWYjULub0R0wCjyA2MysMQcM\nMzMrxQHDzMxKccDAfRhmZmU4YBTIIzHMzBpywDAzs1IcMMzMrBQHDBq8ys/MzHbhgGFmZqU4YBR4\n4J6ZjUTRpkc9HTDMzKwUBwzaF53NzFqp3a0ilQUMSdMk3SCpW9IaSWfWSXOKpDvy52ZJRxT2PShp\ntaRVkpZXlU8zMyunsnd6A9uBsyJipaRJwApJSyJibSHNA8DrIuIxSScC84FjCvtfHxGPVJhHMzMr\nqbKAEREbgY15eZukbmAqsLaQ5ubCIbcCB1eVHzMzG5y29GFImg7MBG7rI9npwI8L6wFcJ2mFpLl9\nnHuupOWSlm/evHlA+XMPhplZc1U2SQEgaSKwEJgXEVsbpHk9KWC8prB5VkRskHQAsETSXRFxU+2x\nETGf1JRFV1eX634zs4pUeochaRwpWCyIiEUN0rwcuBiYHRGP9m6PiA355yZgMXB0lXlNean6G8zM\nRq5SdxiSuoDXAs8DngTuBP4zIrb0cYyAS4DuiLiwQZpOYBEwJyLuKWyfAHTkvo8JwPHAueWKZGZm\nVegzYEg6DfgY6WmmFcDdwJ6kpqOzJd0J/G1EPFTn8FnAHGC1pFV52zlAJ0BEXAT8HbA/8M0UX9ge\nEV3AgcDivG0s8L2I+MnAi9mEG7LMzJpqdocxgdSX8GS9nZKOBA4DdgsYEbGUJu8oj4gzgDPqbL8f\nOGL3I8zMbKj0GTAi4htN9q/qa/9II3dimNkIFNGePthmTVJfr9kUwCPADfkOwszMhki73xLarElq\nRZ1t+wFflvT/IuJrFeSp7cKdGGZmTTVrkvpOve2SLgJuBkZFwDAzs+YGNA6jUSf4SOceDDOzxvo9\n0lvSWNLjsutbnx0zMxuumnV6b2P3UQpPAj8HPlRVptrNr8MwM2uuWR/GpHZlxMzMhrdmdxidfe1v\nMMJ7xPIwDDOzxpr1YVxTZ1sAU4ADgDEtz5GZmQ1LzZqkXlZcz++1OBt4I/DFynJlZmbDTqnHaiUd\nJunfSC84WgEcHhH/p8qMtZP7vM3MmmvWh/FS4HPAS4B/BE6PiB3tyNhQaPcwezOzVmjXRW+zPozb\ngXWkvoyjgaOLE/RFxMeqy5qZmfWl3Q/qNAsYp+MWGzMzo3nA+G5EbG9LToaQB+6ZmTXXrNP7l70L\nkkZNJ3cjHodhZtZYs4BRrEJnVZkRMzMb3poFjAE31kiaJukGSd2S1kg6s06aUyTdkT83SzqisO8E\nSXdLuk/SZwaaDzMza41mfRgzJN1ButM4NC+T1yMiXt7HsduBsyJipaRJwApJSyJibSHNA8DrIuIx\nSScC84FjJI0BvgG8iTQr7jJJV9Uc2zJ+gZKZWXPNAsaLB3riiNgIbMzL2yR1A1OBtYU0NxcOuRU4\nOC8fDdwXEfcDSLoCmF08tgruwjAza6xZwHgoou9niCSpRJrpwEzgtj6SnU4aSQ4psKwr7FsPHNPg\n3HOBuQCdnX3OlWhmZoPQrA/jBkkfrZ21VtJ4SW+Q9B3g1L5OIGkisBCYFxFbG6R5PSlgnN27qU6y\nukEpIuZHRFdEdE2ZMqVJcczMbKCa3WGcAHwQuFzSIcDvgD1Js9ReB3w1IlY1OljSOFKwWBARixqk\neTlwMXBiRDyaN68HphWSHQxsaF6cgfE4DDMbyVIjT/WN6s1mq/0D8E3gm7nynww8GRG/a3ZipTlE\nLgG6I+LCBmk6gUXAnIi4p7BrGXBYDlIPA+8G3luiPIPjTgwzs4ZKv9M7Ip4md2KXNIv07u/Vknrv\nQs4BOvP5LgL+DtifFJAAtufmpe2SPgL8lHQ3c2lErOnHd5uZjXrtvsYtHTD6KyKW0qQ8EXEGcEaD\nfdcC11aQNTMzG4BS78MY7dyFYWbWXNkXKB1eZ9txLc/NEPP7MMzMGit7h3GlpLOV7JUnIvxSlRkz\nM7PhpWzAOIb0mOvNpCeYNuDJCM3MnlXKBoyngSeBvUjjMB6IiJ7KctVuHohhZtZU2YCxjBQwXgm8\nBniPpB9Ulqsh4vdhmJk1Vvax2tMjYnle/g0wW9KcivJkZmbDUNmAsal2Ping563OjJmZDV9lA8Y1\npOEKIvVhHALcDbykony1lXswzGwka1cdVipgRMTLiuuSjgI+VEmOzMxsWBrQSO+IWEnqAB9V3Odt\nZiNJux/UKXWHIekThdUO4ChgcyU5MjOzYalsH8akwvJ2Up/GwtZnZ2h4GIaZWXNl+zC+UHVGzMxs\neOszYEj6EX10wEfEW1ueoyEkj9wzM2uo2R3GV9qSCzMzG/aaBYwHIuKhtuRkCIU7MczMmmr2WO0P\nexckjZpObjMz679mAaPYqP+C/pxY0jRJN0jqlrRG0pl10syQdIukpyR9smbfg5JWS1olaXntsVVw\nD4aZWWPNmqSiwXIZ24GzImKlpEnACklLImJtIc0W4GPAyQ3O8fqIeKSf32tmZhVoFjCOkLSVdPG9\nV14mr0dE7N3owIjYCGzMy9skdQNTgbWFNJtIExu+ZRBlGDT3YJjZSNaubtg+A0ZEjGnFl0iaDswE\nbuvHYQFcJymAb0fE/AbnngvMBejsrJ1Q18xs9Gr3UIABzSXVH5ImkkaFz4uIrc3SF8yKiKOAE4EP\nSzq2XqKImB8RXRHRNWXKlEHmdVCHm5mNapUGDEnjSMFiQUQs6s+xEbEh/9wELAaObn0OzcysrMoC\nhtK90iVAd0Rc2M9jJ+SOciRNAI4H7mx9Ls3MrKyykw8OxCxgDrBa0qq87RygEyAiLpL0XGA5sDfQ\nI2kecDgwGVic2+fGAt+LiJ9UlVGP2zMza66ygBERS2kytCEifgMcXGfXVuCIKvLVF3kkhplZQ5V3\nepuZ2ejggGFmZqU4YOCBe2ZmZThgFLkLw8ysIQcMMzMrxQHDzGyEizY1rDtg4BcomZmV4YBR4Lmk\nzMwac8AwM7NSHDDMzKwUBwwzMyvFAaPAXRhmZo05YJiZWSkOGGZmVooDBn4fhplZGQ4YBe1+obqZ\n2UjigGFmZqU4YJiZjXDtalavLGBImibpBkndktZIOrNOmhmSbpH0lKRP1uw7QdLdku6T9Jmq8gnt\nm7jLzKyV2t2KXtk7vYHtwFkRsVLSJGCFpCURsbaQZgvwMeDk4oGSxgDfAN4ErAeWSbqq5tiWcw+G\nmVljld1hRMTGiFiZl7cB3cDUmjSbImIZ8HTN4UcD90XE/RHxR+AKYHZVeTUzs+ba0ochaTowE7it\n5CFTgXWF9fXUBJvCuedKWi5p+ebNmweTTTMz60PlAUPSRGAhMC8itpY9rM62uh0NETE/IroiomvK\nlCkDyqPHYZiZNVdpwJA0jhQsFkTEon4cuh6YVlg/GNjQyrzV42EYZmaNVfmUlIBLgO6IuLCfhy8D\nDpN0iKTxwLuBq1qdRzMzK6/Kp6RmAXOA1ZJW5W3nAJ0AEXGRpOcCy4G9gR5J84DDI2KrpI8APwXG\nAJdGxJoK82pmZk1UFjAiYilNnlSNiN+Qmpvq7bsWuLaCrO3+Xe34EjOzEc4jvQvkkRhmZg05YJiZ\nWSkOGGZmVooDBh6HYWYjU7ub0R0wzMysFAeMAg/cMzNrzAHDzMxKccDA78MwMyvDAcPMzEpxwDAz\ns1IcMMzMrBQHDDMzK8UBAw/cMzMrwwGjwOMwzGwkatdFrwOGmdkI1e6LXAcMMzMrxQHDzMxKccAo\n8AuUzMwaqyxgSJom6QZJ3ZLWSDqzThpJ+rqk+yTdIemowr4dklblz1VV5dPMzMqp7J3ewHbgrIhY\nKWkSsELSkohYW0hzInBY/hwDfCv/BHgyIo6sMH9mZtYPld1hRMTGiFiZl7cB3cDUmmSzgcsiuRXY\nR9JBVeWpkfBADDOzptrShyFpOjATuK1m11RgXWF9PTuDyp6Slku6VdLJlWcSj8MwM+tLlU1SAEia\nCCwE5kXE1trddQ7pvdzvjIgNkl4A/EzS6oj4dZ3zzwXmAnR2drYw52ZmVlTpHYakcaRgsSAiFtVJ\nsh6YVlg/GNgAEBG9P+8HbiTdoewmIuZHRFdEdE2ZMqWFuTczs6Iqn5IScAnQHREXNkh2FfD+/LTU\nq4DHI2KjpH0l7ZHPMxmYBaxtcI5BcxeGmVlzVTZJzQLmAKslrcrbzgE6ASLiIuBa4CTgPuD3wAdy\nuhcD35bUQwpq59c8XVUJd2GY2UjUrreGVhYwImIpTergSI8nfbjO9puBl1WUNTOzUaHdF7ke6W1m\nZqU4YECbbubMzEY2B4wCeSCGmVlDDhhmZlaKA4aZmZXigIHHYZiZleGAUeAeDDOzxhwwzMysFAcM\nMzMrxQGD9g2rNzMbyRwwCjwMw8xGonY9uOOAYWY2QrX7ItcBw8zMSnHAwOMwzMzKcMAo8FxSZmaN\nOWCYmVkpDhhmZlaKAwZ+H4aZWRkOGGZmVkplAUPSNEk3SOqWtEbSmXXSSNLXJd0n6Q5JRxX2nSrp\n3vw5tap8mplZOWMrPPd24KyIWClpErBC0pKIWFtIcyJwWP4cA3wLOEbSfsDfA12kFqMVkq6KiMcq\nzK+ZmfWhsjuMiNgYESvz8jagG5hak2w2cFkktwL7SDoIeDOwJCK25CCxBDihqryamY1k7eqHrfIO\n4xmSpgMzgdtqdk0F1hXW1+dtjbbXO/dcYG5efULS3QPO5wUDPXJQJgOPDMk3Dx2X+dnh2VbmISvv\npMHVXc8vm7DygCFpIrAQmBcRW2t31zkk+ti++8aI+cD8QWVyCElaHhFdQ52PdnKZnx2ebWV+NpS3\n0qekJI0jBYsFEbGoTpL1wLTC+sHAhj62m5nZEKnyKSkBlwDdEXFhg2RXAe/PT0u9Cng8IjYCPwWO\nl7SvpH2B4/M2MzMbIlU2Sc0C5gCrJa3K284BOgEi4iLgWuAk4D7g98AH8r4tks4DluXjzo2ILRXm\ndSiN2Oa0QXCZnx2ebWUe9eVVeKpWMzMrwSO9zcysFAcMMzMrxQGjIpLGSPqVpKvz+hskrZR0p6Tv\nSBqbt+8raXGeGuWXkl5aOMc+kn4g6a48xcqfDlV5ymhRmT+ep5K5U9LlkvYcqvI0I+lBSaslrZK0\nPG/bT9KSPKXNkvzQxqiZBqcVZZZ0pKRb8u/5DknvGsoy9aVVv+O8f29JD0v6l6EoS0tEhD8VfIBP\nAN8DriYF5nXAi/K+c4HT8/KXgb/PyzOA6wvn+A5wRl4eD+wz1OWqssykwZkPAHvl9SuB04a6XH2U\n90Fgcs22fwQ+k5c/A1yQl08CfkwaY/Qq4La8fT/g/vxz37y871CXreIyvwg4LC8/D9g4XP9vt6K8\nheP+Of99/MtQl2ugH99hVEDSwcBbgIvzpv2BpyLinry+BHh7Xj4cuB4gIu4Cpks6UNLewLGkR5OJ\niD9GxO/aVIR+a0WZ876xwF75buQ5jLzxN7NJgZ788+TC9tE6DU6/yhwR90TEvQARsQHYBExpd6YH\nob+/YyS9AjgQuK7dmW0lB4xqfA34NNCT1x8BxknqHQX6DnYOTLwd+AsASUeThukfDLwA2Az839zM\nc7GkCW3K/0AMuswR8TDwFeAh0lXn4xExnP/AArhO0gqlKWoADow0loj884C8fdDT4AwTrSjzM/Lv\nfzzw60pzPXCDLq+kDuCfgE+1Kc+VccBoMUl/DmyKiBW92yLdj74b+KqkXwLbSLP5ApwP7Ks0VuWj\nwK/yvrHAUcC3ImIm8D+k299hp1Vlzm3Bs4FDSE0VEyS9r30l6bdZEXEUadblD0s6to+0g54GZ5ho\nRZnTznT1/e/AByKip07a4aAV5f1r4NqIWFdn/4jSlskHn2VmAW+VdBKwJ7C3pO9GxPuA1wJIOp7U\njkuk+bU+kLeL1Ib/AKk5Zn1E9E7Y+AOGacCgdWV+M/BARGzO+xYBrwa+297ilJObU4iITZIWA0cD\nv83NLhtzhbgpJ+9rGpzjarbfWHHWB6xFZSY3uV4D/E1uvhmWWlTePwVeK+mvgYnAeElPRMRw/Xtu\nyHcYLRYRn42IgyNiOukK+2cR8T5JBwBI2gM4G7gor+8jaXw+/AzgpojYGhG/AdZJ+pO878+A4rtE\nho1WlZnUFPUqSc/JgeTPSNPiDzuSJii954XcVHg8cCdpupveJ51OBf4jL4/4aXBaVeb8u19Mau//\nflsL0Q+tKm9EnBIRnfnv45Okco+4YAG+w2inT+Wmmw5SM9PP8vYXA5dJ2kEKCKcXjvkosCD/gd1P\nviofQfpV5oi4TdIPgJWk5qtfMXynWzgQWJziGmOB70XETyQtA66UdDopAL4zpx8N0+C0pMzAX5Ie\n6Nhf0ml522kR0TuF0HDRqvKOGp4axMzMSnGTlJmZleKAYWZmpThgmJlZKQ4YZmZWigOGmZmV4oBh\nQ0LSVyXNK6z/VNLFhfV/kvSJFn/nE608Xz7nkXnAYu/65yV9ssRxkvSzPICtd9vbJIWkGRXkc7qk\n97b6vIXzf0TSqHuM1HblgGFD5WbSKG7yXDuTgZcU9r8a+MUQ5Ku/jiQ9e99fJwG35wGLvd4DLCUN\nfmy16UDdgKE87fwgXQp8rAXnsWHMAcOGyi/IAYMUKO4EtuURz3uQBvf9StJESdcrvVdjtaTZAJIu\nyFMtkNc/L+msvPwpScuU3knwhXpfXi9NvgrvlvSvSu9quE7SXnnfK3PaWyR9Wel9HeNJ07a/S+l9\nCb3vdThc0o2S7pfUqBI9hZ0jhJE0kTTFyukUAoak4/K5et+LsiCPgkfSSXnbUqX3MPS+h+R1OT+r\nlCaunESav+u1edvHJZ0m6fuSfkSaXE+Fcq3uLUv+/p9LulLSPZLOl3SK0ntMVks6FCAifg88qDSZ\noI1WQz2/uj/P3g/pXQOdwIeAvwLOI115zyJNFwJphO3eeXkyaRStgJnAzwvnWpvPdTxpdLhIF0RX\nA8fmNE/kn3XTkK7CtwNH5nRXAu/Ly3cCr87L5wN35uXTKLzfAPg86e5pj5zfR4Fxdcr+38Ckwvr7\ngEvy8s3AUXn5OOBx0rxEHcAtwGtIc3atAw7J6S4Hrs7LPyJNmgdp7qKx+TxXF77vNNLcR/vl9beT\nplYfQxrh/BBwUD7ud3l5D+Bh4Av5mDOBrxXO+TngrKH+f+VPdR/fYdhQ6r3LeDWpIrylsH5zTiPg\ni5LuAP6TNIX0gRHxK+AASc+TdATwWEQ8RAoGx5OmFVlJekHTYTXf21eaB2LnFBUrSO/q2IdUuffm\n6XtNynVNRDwVEY+QJqY7sE6a/SJiW2H9PcAVefmKvN7rlxGxPtKMrqtIgW0GcH9EPJDTXF5I/wvg\nwnx3s09EbKe+JbFzGpLXAJdHxI6I+C3wc+CVed+ySHMiPUWahrx3yvnVOS+9NpFmGbZRynNJ2VDq\n7cd4GekKfh1wFrCV1CYOqelmCvCKiHha0oOkq2tIM/i+A3guOytbAV+KiG/38b1100iaDjxV2LQD\n2Iv601b3pfYc9f7OtkvqiIgeSfsDbwBeKilIV/kh6dN9nK9hniLifEnXkO7WbpX0xgZJ/6ew3FcZ\ni9/fU1jvYdey7Qk82cd5bITzHYYNpV8Afw5syVe2W4B9SNNB35LT/C/SuzaelvR60suWel1Bau9/\nByl4QJrp9YO5TwBJU5VnzS0ok+YZkd6Et01pBlLYtVN6GzCpP4XO7ia9JIuc/8si4vkRMT0ippGm\ne39NH8ffBbwgBzmAZ96LLek2XE+xAAABQklEQVTQiFgdERcAy0l3I83yeROpL2aMpCmkJrpf9rNM\nLyIFfhulHDBsKK0mtfPfWrPt8dycA7AA6JK0nHS3cVdvwohYQ6oEH46db0C7jtRkdIuk1aRAsktF\nWSZNHacD8yXdQroafzxvv4HUyV3s9C7jGna+B+M9pOm+ixbS4KmmXIYnSS/m+YmkpcBvC3malzuv\nbydd8f8YuIN0V3O7pI/XOeXinOZ24GfApyNNsd8fs0jNhjZKebZasxIkTYyIJ/LyZ4CDIuLMQZzv\nINJdxZsGm6f81NQ3gHsj4qsDPd9gSJoJfCIi5gzF91t7+A7DrJy35LuIO0lvEfyHwZws3xH9qwoD\n9wbgfyu95nYNqemur36bqk0G/nYIv9/awHcYZmZWiu8wzMysFAcMMzMrxQHDzMxKccAwM7NSHDDM\nzKyU/w+EyEaAyLwn9QAAAABJRU5ErkJggg==\n", | |
"text/plain": [ | |
"<matplotlib.figure.Figure at 0x7fbe3f7b43c8>" | |
] | |
}, | |
"metadata": {}, | |
"output_type": "display_data" | |
} | |
], | |
"source": [ | |
"# ABmag gave me weird plot, so using FNU instead.\n", | |
"obs.plot(flux_unit=units.FNU, bottom=2e-27, top=2.4e-27)" | |
] | |
}, | |
{ | |
"cell_type": "code", | |
"execution_count": null, | |
"metadata": { | |
"collapsed": true | |
}, | |
"outputs": [], | |
"source": [] | |
} | |
], | |
"metadata": { | |
"kernelspec": { | |
"display_name": "Python 3", | |
"language": "python", | |
"name": "python3" | |
}, | |
"language_info": { | |
"codemirror_mode": { | |
"name": "ipython", | |
"version": 3 | |
}, | |
"file_extension": ".py", | |
"mimetype": "text/x-python", | |
"name": "python", | |
"nbconvert_exporter": "python", | |
"pygments_lexer": "ipython3", | |
"version": "3.6.3" | |
} | |
}, | |
"nbformat": 4, | |
"nbformat_minor": 2 | |
} |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Don't leave comment here. I won't get notified.