Skip to content

Instantly share code, notes, and snippets.

@saimn
Last active April 21, 2017 19:39
Show Gist options
  • Save saimn/a381dca41c1be3df7e4facdd80d28f1a to your computer and use it in GitHub Desktop.
Save saimn/a381dca41c1be3df7e4facdd80d28f1a to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:25.012422Z",
"start_time": "2017-04-21T19:38:24.553594Z"
},
"collapsed": true
},
"outputs": [],
"source": [
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:25.035089Z",
"start_time": "2017-04-21T19:38:25.014429Z"
},
"collapsed": true
},
"outputs": [],
"source": [
"%load_ext autoreload\n",
"%autoreload 2"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:26.491684Z",
"start_time": "2017-04-21T19:38:25.037170Z"
}
},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import warnings\n",
"from astropy.io.fits.verify import VerifyWarning\n",
"\n",
"from mpdaf.obj import Spectrum\n",
"\n",
"warnings.simplefilter('ignore', VerifyWarning)"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:26.573478Z",
"start_time": "2017-04-21T19:38:26.494046Z"
}
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"WARNING: UnitsWarning: 'erg/cm2/s/A' contains multiple slashes, which is discouraged by the FITS standard [astropy.units.format.generic]\n",
"WARNING: FITSFixedWarning: The WCS transformation has more axes (2) than the image it is associated with (1) [astropy.wcs.wcs]\n",
"WARNING: FITSFixedWarning: 'cdfix' made the change 'Success'. [astropy.wcs.wcs]\n"
]
}
],
"source": [
"sp = Spectrum('../mpdaf/lib/mpdaf/data/obj/Spectrum_Variance.fits', ext=(0, 1))"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:27.069024Z",
"start_time": "2017-04-21T19:38:26.576376Z"
}
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEOCAYAAACXX1DeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsnXeYJFW5/7+nOofpnriT08Iu7KokASUpCGJAEBAEBUEE\nwXDN/BS9hnsV9XpV1Ovl6iqgYAAkSEZAguSwZNgFdtkJOzszO7Fz7j6/P6rf06eqq3t6dmd2NpzP\n88zT09VVp05Vdb/ve973Pe9hnHMoFAqFQlEJbak7oFAoFIqdG6UoFAqFQlEVpSgUCoVCURWlKBQK\nhUJRFaUoFAqFQlEVpSgUCoVCURWlKBQKhUJRFaUoFAqFQlEVpSgUCoVCURX7UndgIWhubuZ9fX1L\n3Q2FQqHYpXj22WenOOctc+23WyiKvr4+rF27dqm7oVAoFLsUjLGhWvZTrieFQqFQVEUpCoVCoVBU\nZUkVBWPsK4yxVxljrzDGrmWMuRljjYyx+xhjG4qvDUvZR4VCodjTWTJFwRjrBPBFAAdzzt8KwAbg\nTACXALifc74CwP3F9wqFQqFYIpba9WQH4GGM2QF4AYwC+DCAq4ufXw3g5CXqm0KhUCiwhIqCc74F\nwM8ADAMYAxDmnN8LoJVzPlbcbRxA6xJ1UaFQKBRYWtdTA/TRQz+ADgA+xtjZ8j5cX37Pcgk+xtiF\njLG1jLG1k5OTi95fhUKh2FNZStfTcQAGOOeTnPMsgJsBHA5gK2OsHQCKrxNWB3POf8c5P5hzfnBL\ny5zzRRQKhUKxjSylohgG8E7GmJcxxgAcC2A9gNsAnFvc51wAty5R/xQKhQXhcBiRSGSpu6HYgSzZ\nzGzO+VOMsRsBPAcgB+B5AL8D4AfwN8bY+QCGAHx0qfqoUCjKmZmZAQAEAoEl7oliR7GkJTw4598D\n8D3T5jT00YVCoVAodgKWOj1WoVAoFDs5SlEoFAqFoipKUSgUCoWiKkpRKBQKhaIqSlEoFAqFoipK\nUSgUCoWiKkpRKBQKhaIqSlEoFAqFoipKUSgUCoWiKkpRKBQKhaIqSlEoFAqFoipKUSgUCoWiKkpR\nKBSKeaGvJ6bYk1CKQqFQKBRVUYpCoVDMGzWq2LNQikKhUCgUVVGKQqFQKBRVUYpCoVBsE+l0eqm7\noNhBKEWhUCi2iUgkstRdUOwglKJQKBQKRVWUolAoFApFVZSiUCgU80alx+5ZLKmiYIzVM8ZuZIy9\nxhhbzxg7jDHWyBi7jzG2ofjasJR9VCgUij2dpR5R/ArAPzjn+wLYH8B6AJcAuJ9zvgLA/cX3CoVC\noVgilkxRMMaCAN4F4EoA4JxnOOchAB8GcHVxt6sBnLw0PVQoFAoFsLQjin4AkwD+wBh7njF2BWPM\nB6CVcz5W3GccQOuS9VChUCgUS6oo7AAOAvAbzvmBAOIwuZm4HjGzjJoxxi5kjK1ljK2dnJxc9M4q\nFArFnspSKooRACOc86eK72+Erji2MsbaAaD4OmF1MOf8d5zzgznnB7e0tOyQDisUCh2V9bRnsWSK\ngnM+DmAzY2yf4qZjAawDcBuAc4vbzgVw6xJ0T6FQKBRF7Et8/i8A+AtjzAlgE4DzoCuvvzHGzgcw\nBOCjS9g/hUJRAcbYUndBsYNYUkXBOX8BwMEWHx27o/uiUCgUCmuWeh6FQqFQKHZylKJQKBTzIpkt\n4Mnh2FJ3Q7EDUYpCoVDMi8eHovj2vVuWuhuKHYhSFAqFYl7kCyo1dk9DKQqFQjEvXptMLXUXFDsY\npSgUCkXNJDN53PLq7FJ3Q7GDUYpCoVDUDLeuqKPYzVGKQqFQKBRVUYpCoVDUjCrxtGeiFIVCoVAo\nqqIUhUKhUCiqohSFQqGoGeV52jNRikKhUCgUVVGKQqFQKBRVUYpCsUszMDCAbDa71N3YY1Ar2+2Z\nKEWh2OXJ5XJL3QWFYrdGKQrFLg3nXFm5CsUioxSFQqGoGaWS90yUolAoFNtENBq13B6LxbBli1qv\nYndCKQrFLo9yPe1cJBIJZDKZpe6GYgFRikKxS6OUxI6lltvNGFv8jih2KEpRKHZZCoXCUndBUQGl\nwHcvlKJQ7LKEQqGl7oJCsUew5IqCMWZjjD3PGLuj+L6RMXYfY2xD8bVhqfuoUCiKqIHCHsmSKwoA\nXwKwXnp/CYD7OecrANxffK9QWKJcHArF4rOkioIx1gXgBABXSJs/DODq4v9XAzh5R/dLoVAoFCWW\nekTxSwBfByBHJVs552PF/8cBtFodyBi7kDG2ljG2dnJycpG7qVAo5gPnHLFYbKm7oVgglkxRMMY+\nBGCCc/5spX247lew9C1wzn/HOT+Yc35wS0vLYnVToVBI8HkEKaamphaxJ4odiX0Jz30EgJMYYx8E\n4AYQYIz9GcBWxlg753yMMdYOYGIJ+6jYyVExCoVi8al5RMEY8zHGbAt1Ys75NznnXZzzPgBnAniA\nc342gNsAnFvc7VwAty7UORUKxY5DKfHdh4qKgjGmMcY+zhi7kzE2AeA1AGOMsXWMsZ8yxvZepD79\nF4D3MsY2ADiu+F6hUOwE1Cr7lZLYvajmenoQwD8BfBPAK5zzAqDPcwBwDICfMMb+zjn/8/Z2gnP+\nEICHiv9PAzh2e9tUKBRLh1IUuxfVFMVxnPOypcM45zMAbgJwE2PMsWg9UyhqQAmknRP1XHYvKrqe\nSEkwxvZijLmK/x/NGPsiY6xe3kehUOwZ1Cr+laLYvaglmH0TgHwxJvE7AN0A/rqovVIoFLs0auXB\n3YtaFEWBc54DcAqAX3PO/x+A9sXtlkJRO5UW0FEsDUpB7H7UoiiyjLGPQU9VvaO4TcUmFDsFnHMk\nk8ml7sYeg1ICeya1KIrzABwG4Iec8wHGWD+APy1utxQKxa6MUii7F3POzOacrwPwRen9AICfLGan\nFAqFQrHzsNRFARWKbYaW3FTW645D3ek9E6UoFLs0SknsfDy4YQbj0cxSd0OxgChFoVAoFpRv3P4m\nzrt5BJn8zrOmeTqdRi6XW+pu7LJUq/VkY4xdxBj7AWPsCNNn3178rikUip2N+QzgCjvRYG90dBRq\n3Zptp9qIYg2AdwOYBvA/jLHLpM9OXdReKRSKBWUprOnt9QrG43GDa3FmZgajo6Pb2JeddwJgPB5H\nJBJZlLbz+by49lQqhVQqhfHxcYyNjc19sES1rKdDOef7AQBj7H8B/B9j7GYAHwPAtrnnCoVih7N5\n82Z0dHTA5XItdVcE2WwWDkdpSlYikUA2m0UwGEQul8PExAQ0TUNvby8ACEG3LXDOkU6nAehuKJvN\nBrt94ZbjicVi8Hq90DTd9p6ZmUEymURnZ+ecx05PTyOfzyMQCCxYfwD9moeHhwEAfr9fTExljM1b\naVYbUTilE+Y45xcCeAHAAwD88+zzDmVgYACZzPyCaZlMRvkwFbstnHMUCpVjBpxz5HI55PN5TExM\nVBQk81nhjnOORCIhLNpwOGz4fGRkBJOTk+CcIxQKYXx8HNPT08hmsxgeHgbnHPl8HgMDA0gmk8jl\ncpbXIPe1UCgglUqJ6wEgZAHnHFu3bsXIyAgmJibEtqGhIczMzJT1r1YymQwmJyexefNmsS0WiwnF\nZEUulxNWPln9+Xx+XudNp9NV+0wKqFAoIBaLoVAobPPIqppKXcsYez/n/B+0gXP+fcbYKIDfzPtM\nO5hsNgun0zn3jkW2bNkCl8uFjo6OReyVYqHZWd0JBP0wydJcjPaBUqowACHwafSQTqeF4EwkEvB6\nvYjFYmCMwefzAQAGBwfBOUddXR1isRjq6+uhaRrsdjs458hkMjWNRgxCuyiYASAYDCIUCiEYDIrP\n8/k8otEo7HY7ZmdnhRLYvHmz4b7l83lkMhmDGyWXyyGVSsHr9WJ4eBgdHR1gjGFkZETcj0KhgNbW\nVqGMAF2A02x+Uj65XA6zs7Ow2+2G/tXKli1bRJ8452CMib5WYvPmzVi2bJnYj4S4FVNTU/B6vfB6\nveCci9Hh6OgoOOeYmZlBc3Mzpqam0NvbC03TwDlHJBJBJpOBpmlwOBzI5/PQNA2FQsEwkquFioqi\nuNqc1fYrAFwxr7PsYLZVeOzsQkdhZFd4XoODgwAgfsDbCgnIdDoNTdPgcrkQDocxMzMDh8OB9vZ2\nMMagaRqmpqaQTCbR398PQBdkhUIBU1NT4Jyju7tbWNTd3d0YGRkR7ZPlGY1GEYlE0N/fj3g8jsnJ\nSeECkvsEAJFIBB6PBw6Hw+Br51K/Q6EQCoUCBgYGUFdXh7q6OhQKBTDGkM1mhUUtKzz6vFAoYGZm\nRgjVgYEBBAIBw7m2bNkCxhgymYxwKdHoiNoCSgoEACYnJ4Vw31YymYwQ9PRns9kM57SCc46pqSnU\n19eLe2TVj0gkIv5aWlqQTCaRzWaFMiXhv3XrVmiahsHBQSxfvlx8Ru69QqGAXC4HxhhyuRzsdjvd\n35rCCHM66YrLn54AoE/en3N+WaVjloLJyUk0NzejxusuY1cQOgojj7w5i2c2zuK8Q1uXuisVIQEw\nNDSEnp4eIURkCoUCEokE/H6jRzefz4v9Y7GYEPSMMfT19WFmZgaFQgGZTAZDQ0NwuVzo7OwU50wm\nk5idnTUIMrJIqZ2RkREh0Egw5/N5zM7OQtM0RKNRcd7Z2VnA4TNc28DAAAqFAoLBIJqbmzEzMyNd\nPMR5NU1DNpsFY0wIPjpfPB5HNpsV10u/RXLF0HZ6T8fJwrVQKMBut4tzkFLRNM2wL7XBOYfT6TSc\ng5TWxMQEbDabsPjJ+t68eTOcTieamppgt9uxadMmwzOmERKN0mg7YwypVAqJRAKNjY3is3w+j+np\n6aqKQs7UotEZXS8ZDnKsxW63I5/PY2hoSNx7UibytdI+zc3NNRV4rSWaczuAFICXAew8idEmaLg8\n3yGVYtflj08M4/nRxE6tKACjZWyz2ZDJZMT3lAK4MzMziEQiaG1txczMjMj46e3tFdagLCDJVUTb\nSBgRnHOR2ULCgoQpCRlyUVA7dA56tdlsQlCRK8PmtxvOISuZeDxusKLzhQLy+dJIIZfLwWazif7K\nwpyEFykV+p/OQYqGFCeNJug4WVhnMhnRD5vNJtwtslVNriK6N6SMKM5ALp1wOIy99toLgO7Olv9k\n9xL1c3Z2FjMzM4aAcS6Xw9jYGDjnaGxsNNwj2e2UyWQM7nK6B0ApAC0/R3Kb0bMiRZBOp5HP58W1\nyoqWMQaHwyFGFQs2ogDQRdlPuyrZbBbpdLrMYlsq6Au/WH7rnREKbJK1tb2Mj4/j+dHEnPvlcjlo\nmjave53L5ZBOp7err+TSkH3PJJxGRkawbNkyhMNhpNNpoURSqRSGhobEvowxTE9PIxqNwmazGdwI\nNptNCD66NhImctCS/qf+0HEADAqDYgDyyAOAELIk3Cckq5YEEY0KYrGYISCbLxSQy5UUBV0jCUPq\nnyx0yfqVYxFAKc5C/XY4HGJ/WcHYbDZxDCkn8s/L8QDqv/wqu77I/Ub3I5lMij6Ru8ns2qFtshIG\njDEXABgaGjIoNurP9PQ0/H6/GCVOTk6KZ0yGBcWf5HtFbZBS2Lp1q2jD5XIZJhvKLkKn0wnGWE0/\njFoUxd2MseM55/fW0uDOgDyEy2Qy2LJlCwA9RWxgYADNzc2oq6sT+9aiVGU3wPYyMjKCQCCApqam\nBWmv1mtYKsLhsBh6k998e5FLi88VNPT5fFi2bJlheyqVgtvtrtjfcDiM5cuX19SXqakpRKNR9PX1\nIZlMIpPJCJeP2a2QTqeF4CHLj4Q3CVMSzpqmieArCQOKUchWOCkBt9uNgYEB4UKx2+3IZDJCgZjd\nD/Q5uWBIYNH3iYSUfJ8L0r2mc5ACyufzhkyffC4PrrEy4W8eLcjnJteTz+dDOp1GNltaRNN8XaQU\nZMHv8XjE/rFYTBwvn6O+vl4YA/L1ycFoefQ1NjZmUKRkkctWO90fek8uMFkxydcuKxZSUNQOKRLq\nk/wdoeOpryTH5BEWfe5yueB2u2G32xGPx+F2uw1KOpvNQqvRgqplpycB/J0xlmSMRRhjUcbY4swO\n2Q4qCYtEotzqpC9PMpkUwca5GB4etmxrW5C/PGYSiYTBb0ykUimEQqGy/ePxeM3XsFREo1HhSjH4\nsKswPj4uApHm7XIWC1C9UJ1slcqMjY2V3WPZEiXS6bRlP2j/UCgkLOmZmRmMj4+LayRLWXYtTE1N\niR+8WWjKlj+9ylYqWYkkHLPZrEF4yH8kxHK5nAgWu91ucYxsmVJ7svFEbZIgy2QyZYKbziHvJ48o\nsrmsaNfsoqH2vV4v8vk8nE4nOOdipAAATqcTLpcLwWBQKAAS8uReKhQK4lgaXdlsNthsNgSDQRE0\nLxQK8Pl84hrJ+pbjCXSv6JWuJ5VKifudzWbFM8jn8/B6veIeyPeBnp98jzjX5zXQ/ZddR3Qv5Psv\nKxIavXg8HjHfgpQ//Zn3d7lcwtUUDAYNWWuBQIAU1IKNKC6Dvh7Fy7ya6bYTMDIyYmmxGoSKySKS\nt891eZXmWUSjUSSTyTKrtRrmc9H7WCyGTCaDaDRqSNULhUJIJpOor683HCdfw84EWdlOp9PQx1Ao\nhIaGBqRSKfFFBlA2IkokEmCMlY0AabtMrd/K6elppNNptLe3i+c9Pj6OhoYGpNNpTE9Po6enR1iC\nqVRKPAurEQn5nWn/aDQqfuDkQy71kQtXAo0ogFKwV7ZE6Vgy9uQgrllhMMbg9XpFQJiOkd1KNHIi\nC1nTNHg8HuFOYYwhEAiIe1JXV4d4PG7wwZesdhiuiaxhGunI3+tCgQs3GaXH6m3ojQSDQTDGxCsp\nA4q1OBwO4XJxuVxC0NXX1yMcDgsFSFZy0ZVieE42m018d2Q0TYPT6RQpseQWk7OY6B7TKEdW6nRP\n6Nz0POV7pmkaRkdHxfOj+yM/b7p/nOsxIJ/PJ5QSBaYDgYBIZy7GFVBXV2dQyn6/H5zrCQw0N0Me\nLJjvC7nlao1R1KJNNgN4ZaGVBGOsmzH2IGNsHWPsVcbYl4rbGxlj9zHGNhRfG7bnPNW6TT8Mq5EC\nWQLyvjQZyNxmLBZDPB4XE3sqTQyS26IvJr0fHBzE6OiosLyt+m3eJme1mJVYJBJZsAmEZgve/Jl5\ncmMikUAkEhFWtNx3EsBjY2MYHh7GxMSEGBElk0nDqInuA/mKzdtlSGjL+5hn8dKPKBqNiueZSCQw\nOzsrLP3h4WHRB7lchPyjBCD6JFveNJwn689s6VM75hiCPEqQrVHalslkxEjA5XIZzgfoAjUQCBj8\n/bQP5dCTwKVz0jFer1coD/qchKvf7y/rr3wfzMHwQqFgmF1cKI4QCIfDAafTibq6OgQCgYpGQq0x\nJbofNpsNHo9nXq5hxpjB9ahpGrxeL3w+HzRNQzAYNIy+aAQjj4jo/tvtdjGqoc/l/80jNzmu4nK5\nhJDP5/MYHx83PH+/3y/65vf7DfdMznaikRQpmrlmeTPG6LtUU4JSLSOKTQAeYozdDUA4IPn2p8fm\nAHyNc/4cY6wOwLOMsfsAfBLA/Zzz/2KMXQLgEgDfqKVB85AdKGl4snqsoLQz+Z6Nj48jn8+jp6fH\n0MbIyIjBwiXhwznHli1b0N3dbekXpy8HPdxkMil89tPT0wAg+icLvHA4XCakCFmAR6NRNDSUdOrU\n1BQaGxvLRiDzYXp6GvX19YhGo2hsbCz7IRYK+oxPEgB0nVu3bjVcg+zDJZ8v3et4PC5eo9EoEomE\nGEnJz3FoaEj8oMpGFICYQBWPxw0uEs71FMKBgQEhfMiVRHMGaOKVLPjMSmd8fByZTAbt7e3iWBIk\nZHGSoiDLXQ4C03eInrEcXJWDqQAMrgryl5OB4vF4RJwjEAiIa6K4hnwdPp+v7JmRdUrPxeFwGILh\n8jOz2Wyor69HIpGAy+UC5xyzqahhHyq3Qc9RFvAOh0O4UwDUnBxQSykLn88n/O61MJPIodE7t7iT\nlSoJf7r3JLRpEht952mE5nK5EIlExIhQzlSi50ijRopzABCKjtqVvxfyaLAW6JnUQrHtmgYAtZx9\noPjnhFTWY3vhnI8BGCv+H2WMrQfQCeDDAI4u7nY1gIcwD0UhT6OnjBH60dNMRllwyQKBHtzMzIzw\nH9IEIdkilo+h4Tt9Jme7yFDNFWpLttAMk5SKbc3OzsLtdgslYv6RhUIhIVi2d7CXSqVgt9vFUJe+\nrJFIxHD/zFCpB/lzytqh66dJYLLQk1MH6VXOEZdTP0m40sjCfJz+v36+8fFxsY2ug85HfTIrAtmt\nQ/2Vj5menhZWPY2QqP+kjOgYs/Ut99HswyfhUygU4PF4YLfbxcjWbMF7vV6k02khmBwOh8HNRASD\nQUQiEeF7r5QqblbiJIQ0TbM0LGRhb/4ukLIJBAJlbtDFdIvK36m5eHE0ji/fPoQHL1o9r3PIgXH5\nXJWUE41EAP2+JJNJw0iTlAQp+WQyKZ4pTaCkoDU9/8VIUilwDm2e7dayFOp/bnOPaoQx1gfgQABP\nAWgtKhEAGAdQc5K8LERYcWKP/J4+J/cDuW3kCTnj4+MGP7j8o5cFiNW5GWNlLhBz38zCXRa28vAf\ngMEHTtspS2V6etpgAZt/lJzrs2FrGVGMjY3B4/Ggra1NKDTztcv9n5ycRFtbG5LJZJlCpP1kwQtA\n5G2bM2DMx5jdGfIojM5vfgYFrufrk8A2ty0rJznbhD6XBbisXGSfMl3P9PQ0mpubxXH03GQ3H233\ner0iOEpt0KQw2p/cQHReylKhUUo0GhU+aYJKOZghdwLnvKIFur2CR9Osj6fZ4jLbab+UMTCTQn+j\nLqQpNlMLqdz2T/+q9b7J+9EIr66uDpqmiZErISsiQFcWqVRKfI8XglS2ALdDQyydRybP4bAxnPTH\n1/GP8/eFy157yvicexbjBPXS+wbG2D3b2G+r9v0AbgLwZc65IZuK678Gy68bY+xCxthaxthamlgk\nCxh6L7/S/6lUSsw2NX9GQSlZQJJrif7S6bTwY8sKyOpcVr51WQhSRoVZKMvCUxZslOprvjZyJxBm\ngVkJ2oe+xObz0/8TExPCBSKPogxZLlIqomwRU7aIOQfcnBEiKz35HpGQlbcZsm9yJZ8vjQSpH7Lf\nn6w52k+OH8j9on7S9cr7y/dFDjZTG+Tu4VwPsno8HsO9yuVy8Pl84vqdTqdox+/3w+VyiYClpmnw\n+/1lQp8+s8LtdsPj8RgE1qvjCTz4Zql43PZYql5v7XNLtkdPXP74OC65e9iw7VM3bMJ0vDSKq7WW\nm8NWDPAXv1+h5I4p/knfU5pkSIrNHFz/9WPjyOR1I8Hn88HtdiMQCFgG4efDPzeE8YGrXsPvntqK\nE//4Oj7ypzdw9nUbAQDp3MJVjyVaOOciwsg5nwVQe3pPFRhjDuhK4i+c85uLm7cyxtqLn7cDsMxN\n5Jz/jnN+MOf8YJoWL31mcGVYCWGzIDS/B+RMD6MAjcVimJmZwcDAAEKhUJkCIME7ODiIoaEhw/Hy\nxB2zErFSGGbfNudc+IPN10dCj3zvdXV1GBwcrBjjACBKD5vPYW6XMkJoezgcFv0jpSkH4mUBTILX\nPLFKThGUhS3tZxX8kxUHkZO20fF0Dll5mNMRzcrAvC8vjuBkZcc5F0YGuVzy+Tw8Ho8YiVBQFChl\n13Beckva7XbU1dUZ0j4BCIEiY7Vtvlz2yBi+/0/dwGCMbdfE07n6MhGzdjfNJHJIZkvP7I9rJ3DK\n1a8DAPIFjuHZND7+1w3YEs5gMpbFjS/P4KnhGJ4dieG2daWU6tcn519m3O/UxdxHrnkDP394DKdc\n84b47KnhKB4diOBHD2wxzBGplUQmb7guGStDLRAIGOJGk7Esbn5lBpGU/hulwDgFpyvBOcdLY3GE\nU9ZKbzKexQ8f0J/5tS9Mi+10ntw8V5WqRVHkGWM99IYx1ovtMxaoHQbgSgDruTEwfhuAc4v/nwvg\n1lrbtBL2ZgFeaR/5Vf680iil2nEEnZNKJtNn5lm7slKodG75/LK7St5naGgImzdvxsTEhKWLh/pE\no4dsNmuoNVPJ3Ub9lBc7kS10GdkXb84Ioldqk4SwPJlJFtryfTGPCOTz5iUhbx6xmJVHpX0oXVAe\njVCmizz6MR9Lfn1y+VBtITk+ILu96L1cymKxMbsYzCMUzjn+874RzCRyOGbNOty3YdvKbQPAd+8t\nxQhvWBfFTx7UK6t+5E9v4JePlr4/Vz87hVAqjw9e9RpeHEvg3L+9ibFoFmdftxEf/csGsd/Fdw7j\nF4+M45VxPX7z7/eU2geA+zeG8T/Fdv/+ygzOvnYjfvavURyzZh2ueXYS174whYtuHgAAhFJ53Pma\nbtQcs2Yd8gWOS+7ejO/cO4L7NoSRzc9frH3m5gF88KrX8NXbBxFLlwyyE//wGmaSect4Riydx3BI\nN8roWl8cs56jFUvnDc/jvjdCeGwwitcnU/jSbUM45eo3sG5rAr99civ+/kpJoV5aNAwqkZ+nUqwl\nmP3vAB5ljP0L+oJFRwG4cF5nseYIAJ8A8DJj7IXitm8B+C8Af2OMnQ9gCMBHa2nMasQgxx7mGk3I\nglg+lraZ3Q5W7ct55NSenM9sFvpyCi59Jte0kfeVM2uszmN1DVSrPpvNiqyoWCyG6elpgx+UrsU8\ncpLdXnIZAUAfiZDwlRWRLJzlZwPoQViawUsWOY1S6D5Rtoh8vykQKMcPzOfkrFTKwO/3i0CifP/y\nef2Ha3YlArq/mNZOAHQ3EOWuO51OEcym/lHMRaZSPIhcJBQIXygmYlkMh9I4uMs4QvjUDW/ih+/r\nRnvAiZlEDusndHfqMWvW4arTl4Nz4PwbN+GWc1YiV+B4bCiGhzZF4Cta3j96YAvenErhnb1+HNDh\nwy8fGcNhvX68o6fOYCFOJXIwe0dki/8fG3Xh94839O/hvW+EsaLZjdPeVqpIkMwWhBKoxhduHRT/\nb5pOYTyaRU+DE5ferwvEv786Kz7fEtG/U39YW33pU1khAbqVPZ9lnfIFjs1h/VzPjyZw5l834I7z\n9gUAxDIFPDSYwFmN5aO3Xz8+jnvfCOO3p5bmfF16/xYcu7cxWymZLeAnD43i0cEoVja7ceurM4br\nBHSL/fMbhhIRAAAgAElEQVS3DIr3fQ0uHNjpwxtTSVTjFw+P4T+O76r5WuccUXB9PYqDAFwP4DoA\nb+ecb3eMgnP+KOeccc7345wfUPy7i3M+zTk/lnO+gnN+HOe8tqm80IXma5NJ3Lo+TOcw+LZpm1kQ\n0rFmC9rs6rByDZnPQ39ysFL+TB5NyNa22S9vZdHL/bSy+q38+ZxzMVs4l8uJLB65HXpPGVayi0h2\nyVjdK1lg05oCcts0K7S+vh5er9dQHsHr9Yp1D6pdp/l88jkBIJ83lnimfV0ul8H/73Q6hRuI/MW0\nPykfSi+UA8putxt+vx8ej8eQU19rAUqmaXC53WX9ng8vjMYxHc/i1GteR4FzDIfSOOMvG/D/7hzG\nTCKHn/2rNOdjYCaN0UgGmVwBjw1GDe089GYE371XX7Ph5GvewGl/3oBfPKJb5GRtA8D1L03j63cO\nI5kt4NZ1s7jkbqMlDwCfumUMlz08hs/+Xa+imqkhaHz541txy6vGn/RcAt3Mna+F8O/3bMYnrntz\nXseZmUkY3Ta/eWIrpor3OJTMYWg2jVxxlPHyeALvv3I9nhzW7+c1z07iuN+vNxwfzxQwHCrVpPrr\nyxFhTPzfE+MYj+pKhUYuX75tsGr/ntkcw6PF5/en5ybLlIQVX71jCDe+PI3UHDGIJ4ZjeN8Vr83Z\nHlFTci7nfArAHTW3ukQUCgVc/+I0Hh6I4bT9msuyYCoJ2M3hLD7992HcdU6/SA01p0mSIDa7c2S3\nAh1L2T1yfRa5cJj5/LJ1b66JI+9L+1m5rKwydGTrm+IMshKTZ5zKIxbz+cyK1uwiI0Ui983v91v6\n1x0OB1KplEHI+v1+MUGP2qXzORwOMMbKrHFZ6Go2G1wuJ3K5nFiO0uv1ivkMFIeheIBcAsHn84lR\nms/nM7iD5LgDXYu8Ylst5AtcCJT7zl9pOaK4fd0smn12HNarm+cFzrF+axIdQScaPPpP9Cu3D2Hv\nJjdmk3n85omtuPHlkrD9yJ90n/tZBzajPaCPXi6+cxjH7BXAg28aq+1c89xUTf0GgGyB44NXlYTJ\ny+MJbI0aYxC3ry8Jr+dHjbGzSvzq0fG5d6rCwwOLU0HoztdCQlk+MhDFZY+M4dvHduLYvYNYP5FE\nOsfxzbs349BuH57ebH2t515fUl40Cl4/kcQNL81gRZMbbXVODM3qbqe5hPkrW7etZNBvntg6907z\nZOEWjV1iSkLO+N7sMqJXWQiTZWFlwcrtmwWk+Rh5Up0stGVLV/aJyyUVrEYpssCW5wWYs3nkfehY\n2f1DCosm6FGfzErJrCysFK2siORtlL0lu3qshKLsriMouyedTovaP3IqLVBKNZSfY+nZlKcakiKy\n2WxlE7joumkfOod5hGB2LckuRLo/c/Ef940Y2oul85gIpbC8yY2XxxIYj2VxWdGiX3NqPzh0vzcA\nnLS6AQ6N4dPv0HNHNk7rbh1ZScjcvzGMsw9qEe/NSmJ7+eKtgwi4rQOsx6xZt6DnqsZUfPGzluiZ\n3L8hjGP3DhqEbyUlYcbGgKc3x/CNu/TsLZImm2YqL5FKZPIF3PBS6Tnfv9H6WX793R34b2k0CQAU\npz56eQAPbSodd+zeAXzusDZ89u+bMBGb3z3cLepcxzJ53L9Rz8KRRYjZErYSeLrQKbk7zMFS2c9t\nZeWbLXgrZWO1nXMuLFZZcFZqS871N7uvzP2RrXs5o4j6b3ZR0bGViszJio/SZMlnT5VHJycnDffG\nLEST2QIyef3a7Ha7pVAmi54EtBzw9fl8IlV0PJ7HsxOSEq/wveCcY3MobZlKStlJVsL+6mcnK2aT\n+P1+4TKjY632fXwwio1TKeE6IE784+s4/8ZNeGE0ji/eNogfPVAKOl5084BQEgAQTedx0ysz+Ort\nQxWu0Mjdr4fwQo1W/bZCWTOLzZWn1Va518xdn9JjBN842npJ45+e0IO7i/vUwhPDMdz9WnkxzlrQ\nGBNKAtANmkvvH6lyhD6a/P1TW/FZ6XtQjfeuCOL0/RotP/vOcZ344hFt+NZ7OgEA3z62C41eO646\nfa8ar6DEvBQFY8zHGPsEY+zOeZ9pERmLZPHjh8Z1YWuyyM3KoVxJ6BF6wBhPsHK5WAl7s8VN7cjC\nWd4X0AUglQOg/SlrRm6D81IZYXMcQnYhmV02chuykjH/bx6dWCkdm80mJnjJcwkKhYIo5iYrNysL\n/pK7h/G1O4bwkwd1y8eqtMRckALx+Xy46bUkLn2w5L6o5AZ6bksc51xf7sfeMJXEq1UCqH9cO4kn\nhmKWn1kpuZOvfgOXPz6OZ0diOOvaDRiaTePf79mMT9+0qeI5vlKD8KcRwbqJ6oFJYjSSrdjuL07s\nrXrsPi21lcKYLwGX/pxP2Lc80N9eVx7jeWubPjJc3qT355JjOvCD93VXbP84KQB81oHN8Dg0PHjR\nahy/MoirTl+OVn9xZKkBh/f68bY2L9yOktj71jHWCkXGbLHLOCpMQAQAsw0yNJuuODIAdDflS2MJ\n/PWF6YqjjjvO2wdH9Oky4fBeP+w2hs++03pOssYYTnlrI967Ioh7zi8pR59z/ssl1DLhzskYO4Ux\ndgP0khvvAfDbeZ9pByALSXoPADeviyCcKpVqMAtEhspZUPKEqUrHW7liKln8Y9Esrnx2VljLsgUu\nu7fMI5lKsQmy6uVt8uhBPk6eH2AeXdA55c8LhYKoyEnHy8qMfLDUZ1IUZgv+qeEY1k8k8cCbEeRr\nyN/O5AoifVFmPJrBJ//2Jl4eNwrOfN46iJoxpTtSe9//5xZcfKdRoI5HM4aJWFPFiV3/dssAHqrg\nwvnj2glsnNLdQc9tiePiO4cxGsnik3+zDrI+sHHb0063F6eNob+xck7P5ScvzDohZvoaXLjpEytx\n9kHNAICPHaBnPH3mnctwxWnLcfsn98ENZ6/Aneftg1Pf2ggNRul6QLsPR/bV4T/e24U7ztsH91xg\nHA0c0l2aAHhO8RyALiT7G9247qwVuOXclbjvglX44ft7RKrw5w9vxU8+2IP3rqzHneftg8++sxWr\nlxndl3Pxnr0DuPfTq2re/7oXp6t+PhkvV/QXv6u0UunlJ/fB57Th0vd14y2tHrx7uW6QMcbwlaPa\ncMPZK8S+PzBlNDlNKdJu+/yy7yoqCsbY8YyxP0Cv8/QRANcAmOGcn8c5v31eZ1kCSNA9uTmOPzwf\nwstbUxgJZ/DnF0MGRSDHIeg4q3iBjHnkUWkUIrdN+70ykcatr5esVbcpG2YkkoVms5ed8843Ykhn\nyzOw/vuxabw4nqrYb/lPTmeV3UlmxQFArAMg1xeSFRCdg1xP5OJijCFf4HhqOIonhqJlrpe7Xw+h\nwHlVhZEuCvg1T201+L7fmExhOJTBVtOkrsquJ+P7436/HlPxLMKpHFI5jkcGIvjUDW/izekUPvbX\njTj7uo0iI+XKZyaRL3C8ujWJlyxGH8esWYern50SowYKOlfjB/dXz23fFj66XyO+9q7Kyx4/eNFq\nXH/WCqxa5sHZB+qC1GljuOr05ThjP11oH78yCFvRMu4KOnH08gA+tMo4AqCRwVxQii1xSLcPjV67\nENAfP6AZK5vdOGP/ZnidNvhdNjT7HPA6bfjCEW04sNOLfYujmwcvWo3W4qjj3csD8DltcNo0/O2s\nkkCs99jxrWM68PdzVpYJQyLotpe5GE97WxMO7dZTV71OGz66fxMuP6UfDR4b1pzaj9+cMrfiPKqv\nevHCfVvmVjxy/amP/XVj2ecnrGrA+Ye0oNFrx+rWUsmS/z25H8evLD2jk1Y3otnnwEXvWIb371OP\nI/ur9+3K0/fC7+fh3qv27f4HgEcAHMk5HwAAxtivam55B0LfgXyhPEZx6b+mxP9PjyRw/SsRnLVf\nEJxzXP18CId0umCX3CDmEQUAPDwYx6NDcXzjyCaxj7yv+VgSxGRhk8tGt/yNQsc8J+Df7tInot9y\nZqdh+++fC2P/riC6fbqwdjgccLvdeHzzKDoCSaxuMo5QaL0Bp9NpcFFRm5R9ZHUNdXV1Za4hn88n\nChEaRhHFQ+ka7Q4nXhxLWKZTAsDPHx7DvzZFsHYkjitPW44rnp7AGfs3Yf8OH37+8Cj+7fA2MUOW\ngnmc68/VYbO2gu7aEEPzVo7T9is9n3yhfHIWoGfcRNP6s6U00Qtu1IV9PFPAHVIWz5Zijvxt62bw\nmXcuQyJTQL3HjtFIBmae3bK4sQEA2KvJhUvf1y0EyleOaseRfXVo9NqhMeCn/ypNaFve6BLui2VF\n9wtlVdk13dr+zGFuTCayOLy3NBmizmXD996rW6OdASfWPKV/H798VJuY3U386v3L8KNHZw2K+89n\n7i1mPt9wegeaG3Vh1uCx44azV8DvsmHNRyoLqHPf3oJPSAF5K4Ke0nfzkC7fgs5NufmcfQAA6VwB\nHQEHDur04Y71IQTdNtz0iZWYTebw5+em8OWjSsr5lnNW4qXxBL537wg4dJdQNJ1Hdp6zn81Q7OHs\ng1oMSQrVOPOA5rl3AtARmF9912qK4iAAZwL4J2NsE/Q5FAuzFugi8eunZw3mpSzs9UwXYyDupvVR\nZGDDu7uNE+Lof3IJPTWSwBMjKcsRBuccP3x4Gp87tAENbr16JE3UkoO69L/5K+10OrE1FMd3/7EZ\nP35Pk6FdSvUkmKaBsZIi+HrRfTIQyomMKM45onk7HtuUwAf3rdffF6uukp9f7hOgVx0Nh8NiYpnN\nZkOBc6RzHHe/HsJsIofzD10mlE4+n0d9fT0e2BjGZY+M4W9n9IgJeb96bKshH9+KtSO6UD2/KKDT\neY6ftXtxx/oQPry6EU0+49dyzZMTuP6laXz+MGtf7PWvxgDoo7TLLVIDz/jLG2jy6sLSPMIxs16K\nB9Ds1XwB+N/HtuL29bN48KLVOOvacstvR/Dd47rQVqf/wE/Ytx4nrS6VlX//PvW47oVpMQFszanL\nYTMZ2J6ib14WX985tuSi+PWH+1AnjRzOPKAZJ6xqQK7A0eCx48i+AI6/ojR3oK/BievOWoFwMoeT\nr3kD3zuuy+CXNyv2Zt/c804YY6hgD5TalUY/izWz3WXX8JeP6SOXrxSVgsYYmn0Og5IAgKDHjqP6\nA7jv06swGc+irc6JK5+ewJ+fr56KfMpbdUXwkw/24Pv/HEE8Y/RunHfwglRKWhAqup445y9wzi/h\nnO8F4HsADgDgYIzdzRhbiJnZC85zY6UAkBw7IGRBT7V4ZpJG5UELvbwxnRXB2q3FVLJMnmMsZpy9\nzDnHM6MpnHfLGCbiej2bJ0eSYnIXTeHnXJ+45XTpP/R/u2UAW6NZFDjwvYemMBTK4pbXS1YpuX0e\nHSwt2kMGijhvUeA+uyVuKCB2z8YYfvbwmEEhUKzF5/OJbB/aTpPMPB6PmF381+en8PG/bsCvHxsX\nX3hKXSV+cP8WxDMFkb2VyrM5lYQVz22J478f0gOGn75pE75lKgb32qQuvK2UgEylzydiOYMCqIZc\nLkH2F9N8gcn4wpfOvveCVSKIS3zn2E7x/z8/vQprTu1Hd1D/7hzY4cV5BxstTI0xfPc4Xei/q78O\ndlvl8iCVgv9vbfOit8EYx6hz2YRbrdJxQY8dD160GkfvFUC5KbTwMMawepkHx+w197oVRDQaLVvI\nqlY0xmoqy23TmFDka0dixWP153fNGaVMo5vPWYkHL1qNLx7RBgA4tNuPKyzcQM65NOYOpKasJ875\n45zzLwDoAvALAO9c1F5tI5F0QVhL5i81B5A3CVoAeGwwiqFQyZVQV1cHn8+Hr983gZGidfb6tC4c\n7nwjhs/dOWFof83akqsiksrjhD+8jh89rAtWWnXL6/Ph5Ou24M/PT8Fu03905Pv+/j9HMBjSFdHm\naKnPN748gxdH4/juvSOYThaDzQV9Ylg4lceTw+VZOST4C9KPlVJLza4ku90Ov9+PcIbhu/dsBmMM\n2QLH65NJrB2J4fnRBEIWqZCZXKHs3pLSiaS3PXWSyjwAwGumwm+V6uAsNmGL67/krmGLPa35rMUI\n6HrJvw4AvzqpFw4bw/ErdAV9+yf3wQMXrsJ79g7i4ne14+jlAdg0hpUtpYqwl53YhyYL63zvZjfu\n/tS++M/jK2cJASWDY77YTYLLShFVSQJaUC4/pR+fOqR2i5vSv3cUFJO54NBlsGkM3fUurDm1H184\nos0yntVW58QDF5YC4985tlPEjRYSzvWioubConNR0fXEGOvjnA+aTlIAcC+Ae4tF/To559UTg3cw\nsaL/OZrK4ZybS1276dUI3pwtrXomT6ai4OnzUwz7aBlRWfFPz03hxH1LVstLU6UJbfQjuWtD6YZX\nKCKJfw3o+1z1zCS+/u5SOp6cQw/oOdvE/0nWscujB7EKxYll1z8Twt2vG63nmUQOGvTU29FI6Qch\nF6iTaxGRG+rl8TAeKbpjHh2IVgy40jV/8/4p7NfmxheOaij73MJ1v9shpy3+7IQeXHynUXHccs5K\nnHzNG3hLqwcf3a+pbJZsnSkovF+7PrL90Kp6HLt3AF4pdfGEVQ04YdX8VgKWUz+tOOvA5qopndXQ\nGMM5BzWL2d1WI4yF9ATRaoH0W5NjfTs7FNR/nxRwXtniwcoqAW7GGO69YFXFWNx84ZyLFfccDodY\nZ6fSyLAa1WIUP2WMadCrtz4LYBKAG8De0FegOw66S2rJFYVsIb08of+QZ5J5ZCSDcCBktCYufaCU\nG/27Z3Vr9gf3bwEHsMyv35b7NoQNFtILoyXL1upmf/N+a58kuYj0fs1/Vmm+OPD76h1DOPktDZaW\n3C8fGcP3i/nm0aJln84VkC9wg/AxQ7Ncv3X3sMjPtuJrdwzhshP7sHEmA7tdwz9eN7qYbDYbvnyH\ndQB7d+XtXX788sRe3LZ+FrOJHJ4fTSDoseMH7+tGi8/6p+VxaLjq9OVo9jkMSoMxVvU5bS+hUAh1\ndXW44NDt83vPJWLmI+JSqRRcLpfl9zmdTiOZTKKurk5MuqRCj9uiKKhyQaFQQDqdhs1mE+1Q0Ukq\n1TJfZMORIO9FLcuvUjai0+mct5IoFAqWE0qBUhFKqrxcyzrklah4FZzz0xljqwGcBeBTANoBJACs\nB3AXgB9xzrfN6bcD0ExfaVmZnHK99QQa2kWe3m5VeviRoQRWNjmEP9KKv700g6P3bkBn0Gko6XvF\n05bLa1TlPCkv/5ZXZ3GkhUB/cyaFSCqPDVNJvFks9fDjB7dgIpbD/0mpfuPRDFp8DjGs/c2TusX7\nxHBMZMVY8fxowlD07ScPGe/h9i6yslB4HRoSlYZ2NeLQWM0ZK/t3+LB/hw/rtiYwXHRhys/nf07q\nwxdvG8ReTS68Oa0bMbRK245kLkuSPqPJl9tqtZsFJlXZLRQKiEQihlFtOp1GKpVCIBBAPB43pGaT\nkEulUsjlcoa6W2YymYyoCUaCk0rLu91upNNpoQhoxcB0Og2HwwHOuViCljFW03rdZsLhsIhnErXM\nFwJKriAAFe85XTc9w0wmI5ZOpSxFqjpN9wEoKWJyu1HpGsrEnA9V1R3nfB30MuO7HJ+5Y/sKjxEP\nWEy2+vkTc1dxvGLtNK5YO42DOn1lueXbi1Xmzmgkiw8XF4IhXhpLYNYUrP/YXzeizqXhtk/uKyph\nEnN9tW8pVq80l3Gwsqis+NCqetyx3jQS0fSsooXi+JVB0c9a2K/Ni5fGE/jFib0ieH37eftgJJzB\nBTduwlWnL8f37h0R2USVWN3qNeS5E29r9+LM/ZvwllYPvnPvjhl8UwxOrk0lW9ROp1NY0eSCjcfj\nInONBA61QdmAei2w6uc2e7VisZgo8W5WVpSpR3XCqL9yAUgquU9KhJQAWf75fB6JREJchzzHx2az\niQC23W4X7dNE0VgsZqgnti2QEKf+U7/n8u5RvMSccGPVfiQSMVQ+llfFpOsgaBEtUggOh0NkKtI9\nojpqxb7XNITZbYoC7qw8twPy6ytBSuIbdw3hQ6sacFRxEk40XcCXbh3E0aasESo3XQlayUuOgQB6\nldFKy++e+tZGfKGY3QFAKIprztgL3fV6hk21gnIBt61MMX3+8FZc/rh1htO5b28RiuLGs1fgtD9v\nwNvaPLj43R1waAzXvjCNJ4ejmE7kUODAh4opppQ6+sCFq8AYw15NJav/6jP2AmNsmwvfXfTOVowt\nUACHhCAJ/lQqBU3TRHYdVbjNZrPCeg+Hw7DZbEgmkyJmRVZsIBAQVjtZpYwxZLNZJBIJ1NfXI5fL\nIR6PIxgMGowJElYGSzpnjI/RfuTWCYfDog8k7EnIVqrMrGka0um0EHSxWAyBQACFgl4ZmJQHJXPQ\nvnJ7pBwACIGeTCYN9cTkxbOoj7LCJaVH5etDoZC4RqpQHAzqc7S+engzphOl0RO5f6h4Jd1/edla\nK4OLRkDyOuuapgll6nK5RH9pUq3b7Rbp7lTo0263GxSFzWZDJBJRikJR4unNcTy9OY7bP7mP2PbS\neMJyxnE1tlQQdr9+bLwsk+MbR3fg/ftYL+QDQCgJQPdrmw3Vw3r8eGI4hl9/uM9QuhkA/CZf/l9O\nbce165PYd5kH9R47/vKxvdHgscPj0NBT78TqZV70FM/3laPawNGGT1y3EZk8x3tXBPHeFfqP+xcn\n9hp+qL89tR99DSUf+sFdPjEHZL601Tnwh4/OvxgbQSvukVDw+/1CiAAly5vKppNvmgQKWe+apglh\nDYCEhdiHhAoFPlOplBgNpNNpw4iAMSYWx/J6vfo8nJQxDVlWADS7H9AFp7xErMulPx9ZqdA5zWX7\naaRAsQb5XLQvuXHIlSavPCgrDlpylNqMRCKiLbrnJPypAoFsoZOApn7T/XAB6JAGmPTsyNqn+0/K\nivpqM03+lRcpo76S4pDnRdGxtFokXRNdo7y0al1dnegPrzGyXS3rycE533H5ZIpFZ3uWuKx2vNmd\nBKCikjhz/6aymjf/vHAV7lgfwoEdXkTTebT69ZIOE7Eseupd+Mf5++L9V5bWRTAH/HxODV97Vymb\nTJ51uubU5Yb9adLjpe/rNsStGGM4oKNUNwgA9jFlqPz0hF4cs2YdvnNsZ5myqkYul0MsFkNP0Nr/\nTe4Tr9crhCGg+/DJAs1ms2I1QHJHyMKdziNDVqssYGRhJFvuNH+IBKosqElg6Uqj5CqRz59Op0vC\nSwNyhdJ55TbpWqhfJGjl/tG+ZkvYZrOJwK9ZUMrXJlvS9EpKRHYzkdKgNmhf6hspKfl8pEDo2ux2\nO3K5nGHBLNlVBkCsN0/nJiVM10T3MJFIGGJ98Xhc3CNS1qQA5Ocqr/xot9sRi8XgdDrhcrlQV1eH\nUCgkrjEYDKKhoUFkQxXm8n0VqTaieIIxNgK9lMc/zKmyil2P/3ls/nGbVr+jrLYS8f3ju0QZDBma\ncWrFkX11ZfMTNMYMs4wJmvzlsmu49uN7I+i244NXvYZDu/340pFtNS2AUylddFsDyteftUKUxKgV\n+hHLAp1+7ORDpxn0wWBQWLUAhPVoFkqysDJXDJDdJSRESVCQ4CdBRQKWhCv1k1KqzW6gjDQXQRbu\nVDDSZrPhS4c14eePTZcpCiuLXm6LXjnnaGxsxMzMjFBeNPqgY0hwk+Ck+9TZ2YktW0op3iToZSXR\n0dGBfD6PkZERw7opFAiWA75yLEP+nNw3dJ/kxcNkhULPmO4Bubbkfsn3kMhkMiKIb7PZxOiJ9qfX\nZcuWIZvNYnZ2VvSblKLf7xeKJxwOo6mpSQTrpWdRk6KoNjP7YABfLr79JWPsGcbYL5heLHA+S8sq\ndmGsZowSh/WUZzq9q7+ubNawzFvavPh6hbUCqtFW54THoeGeC/ZFncuGD+1TW3aKuaJwLeTzeSST\n1jO5rZQEuSXM55WXuZW3x+NxZDIZhMNhZDIZQ1YOjRRk4VcoFERpFfLt0z52u11Yj4wxUY+LhJrL\n5TJYzLIgYYyJ6sC03fwqu2YcDofBRUjCV351uVxiUiltNysL8zFOp9PwOaD77un8FBMgAUl9MSsK\nh8MBu92Onp4e9Pb2orOzEx6PR7Rvt9tRX18Pt9st7uOyZcsM1+l0OsW8A3k5XLmv8mfNzc3iXsrX\nSgouk8mIERHFReR+y/dINh7M+9E1t7W1oa2tTZzH5/Ohvr7eMOqiv8bGRjgcDjQ2NqK/v9+Q0dXb\n20uj1+1TFMUfwCDn/Lec85MBHA7gdujzJx5hO9maFLsS/Y0unPLWRrylteTa8FSwfLvr505TbPAs\nXv6932XDtR/fu2z7QZ2+spm6APCfx3eXTSqjVEXy324PTlspsEhUi8dFIhGR/iiTTCYrKhAS3rUS\njUYRi8VQKBQQCoUQjUbFtkgkIqxfoGR55nI5g0KQBaG8nfzNJMzcbrcQUuZXs9+d3DYkKLu6uoRg\nkhUObWtrazMoDbkth8MhlApB7csCUp4gR//39fVh2bJl4lnRMcuWLUN/fz+am5vR09MDxhjcbrc4\nngL0nZ2daG5uNghOl8sl+kMC0OfzGfricrnQ2dmJ/v5+dHV1oaenB42N+miX1oQJBALo7+9He3u7\ncP1R/2SlICsHuneapqGurg7Lli0zKEB5dJZIJAyBe3l0KAt1UpKU/kr31mazoaOjQ7Tr9Xrh9XrR\n0dFheBbNzc2iH26326B4rKD7FI/Ha/pR1hzMLsYrHij+gTHWWf0IRSVohalvS9VNVzS55x1cJswp\nsAsFjSasSg7Mx0qPx+OGYou1JFqQRUXCgiABXmOyBgAYzi23Y26bcv7l9in7JxgMGvalUQSNPDRN\nEz5lGsVQGiq1R7582UUhW6zkiigUSqsKAiW3jFlwk9ujsbERs7Ozwq1D7THG4PF44PV6xbwB2t7R\n0YHR0VHRh66urjKXFN07WZHJ9112HcmxAqaV2qBj/X4/JicnDQKY3CIk6Lu6uqBpGoaH9dnuHo9H\nuOb8fj/8fr0seD6fx+DgoKFPjDG0tZWy62TonDI2m82wv9vtRnt7OxKJBEKhENra2jA4OCgUm3x9\nsi5m9bEAAB+vSURBVLtKbk++XtllJwfjZcXPGBP3YHZ2FoVCAeFw2KDE29raxPumpiZxXpfLhfb2\nUnFCimX5fD6Mjo6itdW6gKZMT08PEolETdkZ25xAzDlf+OL6uyGH9/rLttGPWS789uMP6LOqP3+4\n8QEvRukcWpeAOGP/JsP77x3XhfsuWCVSRF12DQd1GgO9NGHvNimTygr5B8MYq3lUQZOxQqGQwXWT\nTCYruoWskH/QyWRSpDHSdmo/EokgFosJPzFjTOSrWymlTCYjRhG0D7l8zK4QunY5b5/OI+8rW/qy\nm4gEq5VriDGG+vp6dHR0CAUiW/l2ux3BYFBY9GSlOp1OLF++XFjq8qp9LS0t8Pl86O/vNwTCGWMA\nK4kMSu2ULeTe3l5RGNDv9xuEcWdnJzo6OtDR0YGWlnL3JF17Z2enWA+lr6+vbD+5P42NjeL6FwKy\n1jVNrwQdDAZRX18PTdPQ3d2N3t5etLa2or29Hf39/eIYWWlR0U3zCM/sOpRHILIbUHbLUYyiv79f\nFDIFYGlE0fbOzs6y9d6tqGUfYqddM5sx9n7G2OuMsY2MsUuWuj/byvelAm37t+v5cuFwGOl0WhQO\nA0pFxOSlHWUu+1AP/v0Y3YLob3Th8pP7cGZRwAdcxsd4+cl9CLhsOP8Q61jBoT0l5dURcJQtSXn0\nXoEyt9Kl0nKU3zi6A58uloKQ3Uw3fWKl+D8WiyEUCiEcDhsmFZFlTQKaSioAEOmctB/9yVU/5e0y\noVDIMGIwxwYoDZOsehLccvsUuKT+UF/of3nhJ0prlAU7UHLHWPnn6XN6T6/kqiChQv532f1A2zo7\nO8Wx3d3dQpCSH1+2fjs6OtDQYEwSMLss/H6/ELIkwOrq6tDa2mpwQ7W2tqKnp8cwivB6vejr60NT\nU5Oh/9S82+2G11vKESWFZ/X8ZEhgzrUfCXNyMy00TU1NaGxsRGtrK/r7+4Wy9vl8BiFNz7mpqQld\nXV1oaWlBIBAoe8byd4OujdaBl0eJspIxj/CWip1yHgVjzAbgcgDvhV5L6hnG2G3FmeJLwv7t3pqq\nmPqcmqGuvFwB8ktHtuFTN+hrMFSKIRWKQdGL3tGC5U0eXP64ntmzsp4h69Pbuvhd7VjdqpeEPrDN\niUIui189HcZ4VBeMq1u9uLVo6V/5zKRoe9Uyj6GcBwD84fS9kMpxpHIFuOwauoLlMRHOOTwOTWRA\nvWfvgIgVyJkaVNcmEokIN4Y5a4ZGFZqmIZFIiG20kBLl19MxRCgUMizLaoZcP+Y+UfYJjWQYK+Wg\nywtLOZ1Ow2xXum56T3EOeYKUWaCRoDbnvZO7gc4jW/uyNS+7mEhRcM7R1NSE6elpoVBoUSqzRdjX\n14eBgQERpK1kcZpfyW1Blqu5TVk5mO88Y0xY3qWsK/2zZDK5TSUxaqG1tRVbt25dcgFK9PYa1yRv\naGjA9LSe+dXb2wtN0zA9PW2YwwKUUnQZY+IYAHMqyR3NnIqCMXY7yr8fYQBrAazhi1Pv6VAAGznn\nm4p9uA7AhwEsqKLYq8mFXJ4byozLvHt5Hf61SXdV1JhFhpXNbjw/mrCcREYBa1nYnbS6AbetmxWK\nI5NJY82HWtHTpFstlxePNZQ+ID+lxrGyngFw4vBeP25+ZRYHFGf5kE963xY3OoJOPD4YxaHdRvcR\noK+l67QDZ+xfvjIW1YqhWjY//kC3UHSAruzkoDJQqvBpEC6SsiBhKCsOxphICWVMz9wh65TaZIyJ\nzB/aT4ZcP0DJl0/9kNNGzZaqnIEjf077p1Ipg+JKpVJCWcgTruRAruzHNgeX5Ywl+mttbYXH48Gm\nTZvEhLmGhgbhrzcL8JaWFktDgzE9XbJQKBgsefM+8xFA8nMEAJutssgwt2vOBFtIvF5vxZjEzoD8\n/aXXpqYm+Hy+skKITqcTfX19YIwhHo+LEW9n584TBq5lRLEJQAuAa4vvzwAQBbASwO8BfGIR+tUJ\nQC5FOgLgHQt+Fs7xkX29uOxJXVEcv5cP976pW6VrTu3HyhaPKNuQn2dBonsuWCUW3AF0a76tzol7\nLtgXqXisePrS2tGxaARHdrth43m0Sa4g2T0l3BfFUQpZ0GQdA8AvTuwDUMrq+fmHemHXGJJHFOBz\nGLMg+hqqD9dpqVQShubU0Hg8XiYcaJvZKifLWhbgcrEzEp5yTEC2yOk6zYKLtstBYRL6sotH3pes\neHMuvjxhyqCYpZGCuS/ySILcNjSi6u7uxtDQkCGDiGIKMzMz4jgS6g6HA8Fg0OBG6ezstAzEVspo\noYBvNbbHWqXRR/X2t6npeUGxgJ0ZxsqLDFqN8mhfAGhra8PMzIyII+0s1KIoDuecHyK9v50x9gzn\n/BDG2KuL1bG5YPoqexcCgLOtPH2zFjTTD+aIbg++eZw+hCShQfQE7Xh9OoNcjfpCQwFva9O/yBe+\nYxn663XB5HQ4kAJEuujRywOYTegC6+IjmgzlFzjnuOAAP2ZTpR/nr07qxfIml0FgAcC7er1wSa4I\n2k6lq80Lz3/usFa8tdi/aDQq1uCWrXf5FdBHRBe9YxmQzwE2p3CpvKvPB86ME7zMIwqn0ykmDcmz\nTOX+0h9licg+fmrHKvuGRgRUVoGCxXSc3LY881dWVNQvSoWkORBk/ZGioJgKuQry+bxhljQtfDU2\nNmbIrwdKriW32w1N09DQ0GAIUPb09JR9jxba907xhm3l4L5GXP3E0Bx7lRT9ns58lZndbheJBzsT\ntSgKP2Osh3M+DACMsR4AZLYs1lI1WwDIy3R1FbcJOOe/A/A7AHC1r7DM1bzg0GUVy3qftLoBfUE7\nbJoucP/j6Ga8rbX0o5RdKgd0eHHBQfU474AgPnqjsXDeod0+rG604Y8vRrB/q0tYoslkUlh3Hzug\nGaFQCPF43GCRMcZwYIcHKwLGBYUA3d2SSCSwT7PLsJ0WugnHYgZ/+somJw7saTC0Y5UCSpy+XynT\nSS7vbJ6fQIKOhOhpb60XJQIAXfh9892tcLvdBtePPCqQUzrtdrvIJaf9yFI3++blgDBgrIUjCyHz\nfAAS9PIoQ3YR0bXJlr6c3SOPSuSJTNRPqrKay+XQ0tIi6hxRWqfL5RICQh55kOLweDwis2cp2B5L\n/MT9O/CFa5+vad9qefx7AjtTjGF7qUVRfA3Ao4yxN6GbCv0APscY8wG4epH69QyAFYyxfugK4kwA\nH59vI6xgrH3zyQOC+MtLYezd5MJXjmpHKpXCQxv1OkUHtrsN+fb0d+GhzegN2GG3adBY+XDCYwM+\nvK8fH1jhB8Bx6cMzlftTHCnIFi5l3sjWLn3BSKBSv8iXTxkVsu/fXOeHrGqrhU3kfHn6nwQgtUev\n9BePG9fllt02yWQSbrfbUA9HtqJJUFLWiDxSoP5RRpGsHGRr3xxzkH+EsvvHZrOhpaUFY2NjQrjL\nyk6OjcijlUAggFAoZLhul8tlmONgnucgRm1FwUtuIsaMOf3y/ad7sVRKYkcy3xnxip2XqmNDpq9w\nFwKwAno5jy8B2IdzfifnPM45/+VidIpzngPwbwDugb5Q0t8453O6ub7zbmNAtpAvrt52VBO1i/8+\nfhkuPrIF+Xxed4VIWUkkjOSUzg/0O3FAu7uqX5cxBrdDg8uuVaxDLwtmOpeVa0T+o/0YY4b69eaJ\nUfRHwoxqyzidTkQiEYMLjVxe1A85ziG3J98T6gPNX6D+yMKSkIO0JPRlX34gEDBY8fI+sq9fHpVQ\nppD5Tz4nCX4rhQToQlzTNDGpy5x6SAFqcgn19/cLod/V1QWv14tAICB8xzQ/QcZqtTaaLSuPMBR7\nBmYX7K7MXAsXFRhjl3PODwTw4g7qE537Lugr6c0JpWoubzAF/YpS+9AuL4BpMKbvQ/V2OOfYr82N\nzx9SL37A5G+WhTFZovS+2WvDVKIUI5AFzmlvqcNB7W5RMZQmzMguJUqJk/31AMr85gAM1iydg9YX\noPe0H6Wdcs7hdrstJ4vJBepoJGL+Msv9kTN25HkJmUxG+Nfl65eFNm2jEgQ2m02sH2C320VK4cjI\niGHimlUsQs4gMVuqcgyDUk01TS9Jkc/nMTExAafTie7ubpGn3t7ejpGREbS3tyOTyYi4gdvtFnMP\nqCQE1djJZDKi/7KiqQZNXAOALVu2zGuS066KPHlvT6anp2e3cb/V8q29nzH2EQA38510LOktriDX\nbFqnmH7G8qiBftwklHxOG967V7nAA0ojDKpvzznH5Se0wmMDPnWbvnBOwZQCeEC7F/u36b7sfD6P\naDRaJlCpCqjZqjVn7lAf5MC67EIxKynqLwlrCromEgnkcjnU19cjmUxC0zSDC4vaoPLH5n6QsDS7\npswKJplMCuHKmD5rls5BkOKUJ4K1trZiZGRE+PHl9X3r6urgdrtF+QcrXC4XWlpaDAK/p6cHdrvd\nEIwnAdbf3y+ev9vtFnGj7u5uQ18bGhoM/aQRhlVAvRYWavbwzo7fX56GvSeyuygJoDZFcRGArwLI\nM8aS0OUv55wvzkya7YAxhqtPbsOvnwlj7ZakmG+gaRpWNTuxsslV9gM3C1xze+Z9uoN6ts/RfR6s\nHU3h4A6XweqXrXX6SyQSQuBS0FZ2tchCmVwnclaTnCJqpSjMgpu20Xs6J9VPIkUViUSEm0iebEZ+\neGrTfB55JCT7/+XYQnd3t6X1LM89IBwOB7q7u8UCNhSEp5La5OqiGIL8bOQ2ZEuWzu10OssmkdF1\nmRVDLT/s5uZmsZjNfAWBcjspdlXmVBSc8/Ja0jsZvFByy9R77EJBEIwx/OwDHQY3j5USkN8DpRGF\n+X8A+PI79SqUslKQ/yfhKwejZR88uWLI5SIvu2gOBMtBXLPgljN05PgHHU+jEU3TxAiD4he0mIp5\nRrDchtw+Xa/sVqL2qC3aVs3F4na7y4K5VNqAPqdtRH19PYLBIKamppDL5RZkxu+2uEjkgL5CsadQ\ny8xsBuAsAP2c8x8wxroBtHPOn1703s0Ts2vi7R1ufOaQBiG8zBO8zG6fsgXSLf63ClCZ3T4k0OW1\ngElJyQFgOk7e3/yZLLCpPVIgJLSpfVngy9Y9ubnk4xkrzZGgc8mjG8J8LkrxlUs8p1IpMYubMTbn\npKxtccFQP6mYnLLOFYodRy0h+f8DcBhK6akxQFSW2CkgJ4Y5wNjis+PEfYMGt46VNS6/J6VB7+Uq\nnnLKplV8wDwakAWwPNKg81DlTqs/OVOIjjPPF5ALsdH1m7fL/Zffm89jtaaBpmnwer3o6uoq6ztl\nEMmK02azoa6ublFLK2xLbEChUGwftSiKd3DOPw8gBQCc81kAO8/ccgkSIO11DjCUp51C2iYLTLOy\nsFIeVoKyknCn85gFr1kBUODVKs5A26napyy8zemglfpgVhyU1invK1+TWYmQMqORCmUMybEUuU2z\nolYoFLsHtQSzs0yv5soBgDHWAmB+hY92ECSszju4GSft4xfCjz4zjwCsXECy+4ksbZokR0JcnjQH\nWAeRqVqolUKh5QllwRsIBAzzG8wzkOX3cswD0DOJzCumycdU66ccg5A/B/Q5AHIbskIhaEEaeb0E\nhUKxe1GLovgfAH8HsIwx9kMApwH49qL2ahshoedz2OCoKw9UytY8/S/PT6jkfpKrYFqNACopIavR\nCe0jT8CSha9ZuFv1TRbsNHfAatRkVS2VXhsaGsRqa7SdAuptbW1lcYa+vj4MDg5a9g8o1VtarLLS\nip2PnTRbXrEIzOl64pz/BcDXAfwYwBiAkznnNyx2x+ZDg8eG89/eaBDGQHlqp9mKBqxnRcuf+3y+\nMiVg5Y6Sz2flnpJdRjL19fVlQpzOKe/b3t4u1sJljInZw5VGC11dXfB4POju7jb0uaurS9THp7Zo\n1jItflNpHYNAICAykeQCdrTYO60EtqP4/+3df2xdZ33H8fcndnLjJI4d2y0JjtuEkYHSdpQQMiZa\npq2FNhMjrUBbOjTY2FZtqyaxrVQt2Q/+oBKlmoYmBlLoNihDRIy1tBMqaoOmMSqVLJSkabqkpEq7\nNm1+UDWBltGu6Xd/nOc4x9fXJ3Zs32Of+3lJls99zrn3Pufx8fme5znneZ5NI8u46fK5N4CaWd1M\nqptoRBwADsxyXs5Zo3sBH7g4G1iv+aq61evmZqfmoa2LJ+ienp7RyUaKj58WayJ5c1DxOxqNxmhH\nvbwD2bJly0an4sw1j6BarAX09PSMnpgXLFjA8PAwzz77LKdPn2ZwMBuWJJ+gpqenh6NHj47eV1i4\ncCGrVq0a7STX19fHqVOnxsytIGn0UdW8l3RzB7miwcFBTp48Oa75qTiZTzt95tp148rTque7VPVT\nj4FIaD1WUfEqu1XtIt+muZkHzjy11DwAXfE9zTUTyK60izWN/ESf3+9o/qze3t5xaQMDAwwPD7Ny\n5cqznoCXLl1KV1fXuCknm/O6YsWKMZOx53nKh5wuBo4yzTWdYnq7x88fGBho6/fZ5Czsqs1pxZLa\n/EWbT9h5Wv66edrK/Cdvhsm1aoZqbtJqFXia85EHgLxTl5RNGbly5cox480X39toNGg0GuMmLWk+\nMU8UPKRs2IxiZ7YFCxaMNlE1B4GJTvplenp6WvaTmEyQmWlV1GLs7H553SCf3ry66mzYDKrNCGXN\nNYpikGi+YQ2MmWWsVW0kDyLF3szNHeJajeJalE9G39PTw7Fjx0aXW+W9q6urtP9B3tRz3nnnTTjF\nZD7V4mQ03xeZrInyuHr1aj/xZAA0FnaxYdjjPdVJrQJF84m++LtZXsMovh9oOS4QZDWEF154YXSc\no2KQOP/88zl+fPwESfmJs6enpzQI5AFlorwW89Q8plHRVJ44yp+UmqqJ3lPViKHnUiuy2eX+NPVT\nq0CR/y6exPP+DPm65cuXjw4lXjyYy6YglLJ5jnt7ezl16tTocNv59yxdupTBwcExzVvFRwfzG9MT\nqeJKvNFoTBgUzaaj+SLM5r/a3KMoag4Cxdd9fX2jJ+78RJk/pdRq1rGFCxeO6eOQX7U3/yMsX768\n7W30ZmbtUItAkTd7DAwMjLtnkF/Z508fdXd3j3sEtOzqp7u7e0yfgbJ7EsVtfEVlZnVRi0DR3Mmr\n1Ym6+JRS3gch19vby7Jlyyb9fZJKh5tetWrVWUdQNTObL2oRKODMlX5zEJjM0z1DQ0Nj+iGUWbBg\nAUNDQ+OG7iiqqgNaJ3I5m82+2gSKXG9vLyMjI+PSZ+oGW16bWLJkyZgOd2ZmdVWbp55yxR7QxSeg\nent7Z3R2ssWLF7NmzZoZ+zybOtcmqnHJqmXse+5FDwrYQWoXKIoajQb9/f1jpjA1s5lRVd8Za79a\nB4r+/v5J33uw+cdXtNVyoOgclTSwS7pd0gFJj0i6W1J/Yd0tkg5JOijpqml8h+8fdAA3P5nNvqrO\npA8AF0fELwCPA7cASFoPbAUuAq4GPqdsdr2zaj5hjIyMuAOc2axwTa7TVBIoIuL+iHg1vXwIyIea\n3ALsiIiXI+IwcAjYdC7fkXewMzOz6ZkLbTMfAe5Ly8PA04V1z6S0cSRdL2m3pN0nTpxgaGhoSp3m\nzMxscmbtslvSTqDVkKnbIuKetM024FXgK1P9/IjYDmwH2LhxY8zko69mZnbGrAWKiLiybL2k3wHe\nC1wRZx5fOQIUe8utTmlmLflm9tzkv0u9VPXU09XATcD7IuKnhVX3AlslNSStBdYBu6rIo5mZZaq6\n4/tZoAE8kK48HoqIP4yI/ZK+BjxG1iR1Q0S0ns7NzCrh2kLnqSRQRMQbS9bdCtzaxuyYmVmJufDU\nk9m09Pf3n30jmzHuEd953NnA5rWzTTNr1RgYGPDfpUYcKMxsxvX19VWdBZtBbnoysylZscSDAXYa\nBwqbt9xWbtYeDhRmZlbKgcLMzEo5UJiZWSkHCpvX3EvYbPY5UJiZWSkHCjMzK+VAYfOW50Q3aw//\np9m8tWjRIi644IKqs2FWex7Cw+a1rq6uqrPQca5Yt4JuPPp/J3GgMLMp2bx+iHeuXlR1NqyN3PRk\nZmalHCjMzKyUA4WZmZVyoDAzs1IOFGZmVsqBwszMSjlQmNk58cRRnaPSQCHpzyWFpKFC2i2SDkk6\nKOmqKvNnZq151N7OUlmHO0kjwHuA/ymkrQe2AhcBrwd2Svr5iHA3ULM5xsGic1RZo/hb4CagWH/d\nAuyIiJcj4jBwCNhURebMrJybnjpHJYFC0hbgSETsbVo1DDxdeP1MSmv1GddL2i1p94kTJ2Ypp2Zm\nNmtNT5J2AitbrNoGfJys2emcRcR2YDvAxo0bfWljZjZLZi1QRMSVrdIlXQKsBfamNs7VwMOSNgFH\ngJHC5qtTmpmZVaTtTU8RsS8izo+INRGxhqx5aUNEHAXuBbZKakhaC6wDdrU7j2ZmdsacGmY8IvZL\n+hrwGPAqcIOfeDIzq1blgSLVKoqvbwVurSY3ZmbWzD2zzWxKPFd55/Ff3MymZGBgoOosWJs5UJjZ\nlLhG0Xn8Fzczs1IOFGY2ZR7nqbM4UJiZWSkHCjMzK+VAYWZTJomurq6qs2FtUnmHOzObf9auXVt1\nFqyNXKMwM7NSDhRmZlbKgcLMzEo5UJiZWSkHCjMzK+VAYWZmpRwozMyslAOFmZmVUkRUnYdpk3QC\neKrqfMygIeBHVWeiYi4Dl0Gn7z/MfhlcGBHnnW2jWgSKupG0OyI2Vp2PKrkMXAadvv8wd8rATU9m\nZlbKgcLMzEo5UMxN26vOwBzgMnAZdPr+wxwpA9+jMDOzUq5RmJlZKQeKNpH0pKR9kvZI2p3SBiQ9\nIOmH6feKwva3SDok6aCkqwrpb0ufc0jS32keTV48QRl8QtKRlLZH0q8Vtq9jGfRL+rqkA5L+W9Iv\nddJxMMH+d8wxIOlNhf3cI+nHkj4654+BiPBPG36AJ4GhprRPAzen5ZuB29LyemAv0ADWAk8AXWnd\nLuAdgID7gM1V79s0y+ATwI0ttq1rGXwJ+P20vAjo76TjYIL976hjoLB/XcBR4MK5fgy4RlGtLWT/\nOKTf1xTSd0TEyxFxGDgEbJK0ClgeEQ9FdqTcWXhP3dSuDCT1Ae8C/gEgIl6JiJN0yHFQsv8TqdX+\nt3AF8EREPMUcPwYcKNongJ2Svi/p+pT2uoh4Li0fBV6XloeBpwvvfSalDafl5vT5olUZAPyJpEck\n/WOhyl3HMlgLnAD+SdIPJN0haSmdcxxMtP/QOcdA0Vbgq2l5Th8DDhTtc1lEXApsBm6Q9K7iynRV\nUPdH0FqVweeBNwCXAs8Bf1Nh/mZbN7AB+HxEvBV4iayZYVTNj4OJ9r+TjgEAJC0C3gf8S/O6uXgM\nOFC0SUQcSb+PA3cDm4BjqQpJ+n08bX4EGCm8fXVKO5KWm9PnhVZlEBHHIuJ0RLwGfIGsXKCeZfAM\n8ExEfC+9/jrZibNTjoOW+99hx0BuM/BwRBxLr+f0MeBA0QaSlkrqzZeB9wCPAvcCH06bfRi4Jy3f\nC2yV1JC0FlgH7EpV0x9Lekd6wuFDhffMaROVQf7PkVxLVi5QwzKIiKPA05LelJKuAB6jQ46Difa/\nk46Bgus40+wEc/0YqPrOfyf8kFWr96af/cC2lD4IfBv4IbATGCi8ZxvZEw4HKTzNAGwk+0d6Avgs\nqdPkXP8pKYMvA/uAR9I/xaq6lkHK+6XA7rS/3wBWdNhx0Gr/O+0YWAo8D/QV0ub0MeCe2WZmVspN\nT2ZmVsqBwszMSjlQmJlZKQcKMzMr5UBhZmalHCjMzKyUA4WZmZVyoLDak3SJpKck/VHJNj2S/kNS\nVyHtGkkh6c2zkKd+SX88059b+PxFkr4jqXu2vsM6hwOF1V5E7CMbqfNDJZt9BLgrIk4X0q4Dvpt+\nz7R+oGWgUGZa/5sR8QpZT9/fnM7nmIEDhXWO48BFJes/SGGsHEnLgMuA3yMLMnn6mjQz2xck7Zd0\nv6SetO4v0yxk35X0VUk3pjGuvilpr6RHJeUn7k8BP5dmObs9fe5BSXeSDcswIunP0nselfTRwvcf\nkPRFSY9L+oqkKyU9qGx2tHxAPciGyPjg9IvOOp0DhXWKTwENSRc2r0hDPr8hIp4sJG8BvhURjwPP\nS3pbYd064O8j4iLgJPB+SW8H3g+8hWxk0I1p26uBZyPiLRFxMfCtlH4z2aQ1l0bExwqf+7n0uUPA\n7wK/SDaL2R9Iemva7o1kQ3G/Of38FllQuxH4eCGfjwJvn2wBmU3EgcJqT9JmsoHYvknrWsUQ2Qm/\n6DpgR1rewdjmp8MRsSctfx9YA7wTuCcifhYRPwH+La3fB7xb0m2SLo+IUyVZfSoiHkrLlwF3R8RL\nEfEicBdweeH790U2LPd+4NuRDdq2L+UFgNSM9ko+aq/ZuXKgsFqTtBi4jex+wD7g4hab/S+wuPCe\nAeBXgTskPQl8DPiNNJwzwMuF954mm5CnpVQj2ZC++5OS/qokuy+dbX9afP9rhdevtchLA/jZJD/X\nrCUHCqu7vwDuTM1KLQNFRLwAdKWgAvAB4MsRcWFErImIEeAwZ67oW3kQ+HVJi9P9jfcCSHo98NOI\n+GfgdrKgAfAToOxK/z+BayQtSfN3XJvSJk3SIPCjiPi/qbzPrJkDhdVWmiDn3cBnUtJENQqA+8ma\neyBrZrq7af2/UvL0U0T8F9lcCo8A96XvOgVcAuyStAf4a+CTafvngQfTjerbW3zew8AXgV3A94A7\nIuIHJbvbyq+QNbeZTYvnozADJG0A/jQifnsan7EsIl6UtAT4DnB9OuFXQtJdwM2p+cvsnLkzjhnZ\nFbykf5fU1dSXYiq2S1pPdr/jSxUHiUXANxwkbCa4RmFmZqV8j8LMzEo5UJiZWSkHCjMzK+VAYWZm\npRwozMyslAOFmZmVcqAwM7NSDhRmZlbq/wE5UQU0odEVjAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fd8c0a602e8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"sp.plot(noise=True)"
]
},
{
"cell_type": "markdown",
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T13:17:31.986622Z",
"start_time": "2017-04-21T13:17:31.982647Z"
}
},
"source": [
"## With Spectrum1D"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:27.133184Z",
"start_time": "2017-04-21T19:38:27.070862Z"
},
"scrolled": false
},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"INFO:root:Successfully loaded reader \"wcs1d-fits-reader\".\n",
"INFO:root:Successfully loaded reader \"tabular-fits\".\n",
"INFO:astropy:uncertainty should have attribute uncertainty_type.\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"INFO: uncertainty should have attribute uncertainty_type. [astropy.nddata.nddata]\n"
]
}
],
"source": [
"s1d = sp.to_spectrum1d()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:27.157228Z",
"start_time": "2017-04-21T19:38:27.136117Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"Spectrum1D([ 8.8293047 , -2.09127831, -0.3879379 , ..., 15.92895603,\n",
" 16.99621773, 16.33918953])"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"s1d"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:27.183685Z",
"start_time": "2017-04-21T19:38:27.159803Z"
}
},
"outputs": [
{
"data": {
"text/latex": [
"$[4602.604,~4603.2345,~4603.8649,~\\dots, 7183.0286,~7183.659,~7184.2895] \\; \\mathrm{}$"
],
"text/plain": [
"<Quantity [ 4602.60402868, 4603.2344769 , 4603.86492512,...,\n",
" 7183.02859577, 7183.65904399, 7184.28949221]>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"s1d.dispersion"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:27.208580Z",
"start_time": "2017-04-21T19:38:27.185816Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"UnknownUncertainty([ 121.84641266, 123.27562714, 121.95872498, ...,\n",
" 85.33888245, 85.86727142, 84.69287872])"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"s1d.uncertainty"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:27.423205Z",
"start_time": "2017-04-21T19:38:27.210847Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7fd8ade536d8>]"
]
},
"execution_count": 10,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD8CAYAAAB0IB+mAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8FOX9B/DPN3dICEm4BAIkHHJDgHDfchsFqa3iVcAD\nbdHqT1sb6oVVLNV6X5VaUavYUvGqKHKLighBLrkPAwECCXc4EpLs8/tjZzezm9lNsrvZTWY+79cr\nr+zOzO48k8l+95nvPIcopUBEROYXFuoCEBFRcDDgExFZBAM+EZFFMOATEVkEAz4RkUUw4BMRWURA\nAr6IJIrIhyKyU0R2iMgAEUkWkaUiskf7nRSIfRERkW8CVcN/EcBipVRHAD0A7ACQBWC5Uqo9gOXa\ncyIiChHxt+OViDQAsAlAG6V7MxHZBWC4UipPRJoBWKWU6uDXzoiIyGcRAXiPNAAFAOaJSA8AGwDc\nC6CpUipP2+YogKZGLxaR6QCmA0BcXFzvjh07BqBIRETWsWHDhuNKqcaVbReIGn4GgLUABimlfhCR\nFwGcBXCPUipRt90ppZTXPH5GRobKzs72qzxERFYjIhuUUhmVbReIHP4hAIeUUj9ozz8E0AvAMS2V\nA+13fgD2RUREPvI74CuljgLIFRFHfn4kgO0APgMwRVs2BcCn/u6LiIh8F4gcPgDcA+B9EYkCsB/A\nNNi/TBaIyG0ADgC4LkD7IiIiHwQk4CulNgEwyh+NDMT7ExGR/9jTlojIIhjwiYgsggGfiMgiGPCJ\nLKC4tAz/zc4FpzS1tkC10iGiWuyl5Xvw6sp9iI+OwPhuzUJdHAoR1vCJLKCgsBgAcLaoJMQloVBi\nwCcisggGfCIii2DAJzK54tIyLMg+FOpiUC3AgE9kcmcuMG9Pdgz4REQWwYBPZHJseU8ODPhERBbB\ngE9kchLqAlCtwYBPRGQRDPhEJsccPjkw4BMRWQQDPpHJMYdPDgz4RCbHlA45MOATEVkEAz4RkUUw\n4BMRWQQDPhGRRTDgExFZBAM+kcnp5y3nHObWxoBPRGQRDPhEJidi/JishwGfyOSYxiEHBnwiIotg\nwCeyENb2rY0Bn4jIIhjwiSyEN22tjQGfiMgiGPCJiCyCAZ/IQnjT1toY8ImILCJgAV9EwkVko4h8\nrj1PFpGlIrJH+50UqH0RkW9409baAlnDvxfADt3zLADLlVLtASzXnhMRUYgEJOCLSAqATABv6hZP\nBPCO9vgdANcEYl9EROSbQNXwXwDwIACbbllTpVSe9vgogKZGLxSR6SKSLSLZBQUFASoOERnhTVtr\n8zvgi8hVAPKVUhs8baOUUgAM/9WUUnOVUhlKqYzGjRv7WxwicqOMP3pkQREBeI9BACaIyJUAYgAk\niMh7AI6JSDOlVJ6INAOQH4B9EZEfeNPW2vyu4SulZiqlUpRSqQAmA1ihlLoZwGcApmibTQHwqb/7\nIqLqEzDKk11NtsOfA2C0iOwBMEp7TkRBxpQOOQQipeOklFoFYJX2+ASAkYF8fyLyD2/aWht72hIR\nWQQDPpGF8KattTHgExFZBAM+EZFFMOATWQhv2lobAz4RkUUw4BMRWQQDPpHJ6dM4bKVjbQz4REQW\nwYBPZHL6Wj1v2lobAz4RkUUw4BOZHGv15MCAT2QhvGlrbQz4RCbHIE8ODPhEJqdP6Rild46eKcK2\nI2eCVyAKmYCOh09Edc/AOcthU8DUgamYNaFLqItDNYg1fCKLs2m1/rfX5IS0HFTzGPCJTE6fxWE+\n39oY8ImILIIBn8hC2Cbf2hjwiUxOMcqThgGfiMgiGPCJLIQ3ba2NAZ/I5JjRIQcGfCILYfC3NgZ8\nqnX25hfi9IVLoS4Gkekw4FOtM+q51ch86dtQF8M0WKsnBwZ8qpUOn74Y6iKYEm/aWhsDPtUqNhur\no4GmwL8p2THgU63y1nc/h7oIpsb0jrUx4FOtsjf/XKiLYDoM8uTAgE9kcoz35MCAT0RkEQz4RCan\nHzyNrXSsjQGfyEKYz7c2Bnwik2OMJwe/A76ItBSRlSKyXUS2ici92vJkEVkqInu030n+F5eIqqsu\n1OptNoXPtxxhP4waFogafimAB5RSnQH0BzBDRDoDyAKwXCnVHsBy7TkRUQX/3ZCLu+dvxHs/HAh1\nUUzN74CvlMpTSv2oPS4EsANACwATAbyjbfYOgGv83RcR+aL237Q9fs4+WN7RM0UhLknwnCsuxf3/\n2YRT54M3UGBEIN9MRFIB9ATwA4CmSqk8bdVRAE09vGY6gOkA0KpVq0AWh4jgmtKpifTO2aIShIkg\nPtoeThb/lIfTF0owuW/VP8+OLyIzZ3Q2555G04QYnCsuwfwfcrF2/wlszzuLxHpRePTqzkEpQ8AC\nvojEA1gI4D6l1FnRVSWUUkpEDE+lUmougLkAkJGRYeLTTVR9y3ccw9GzRbipX+tQF8Wj7rOWIDJc\nsGf2lQCAu977EQCqFfDDtHjhaEJqsymEhYXucmTRljzsLziHe0a2D8j7ldkUJr76neE6WxBvsgSk\nlY6IRMIe7N9XSn2kLT4mIs209c0A5AdiX0RWcts72Xjo45/8eg9/w4lSCn9dvBM/Hz/vXPbBuoPY\nV3AOWw+dAQCUlCnkny3C2aISn/bhiO1vrN6PF5ftQZs/feEM/j8dPoPdxwpRUFjs34FUw4z5P+LZ\npbsD9n4PLNjkcV1ZEC9r/K7hi70q/08AO5RSz+lWfQZgCoA52u9P/d0XERlbs+84eqQkIi66eh/p\nY2ddc+YHTpxH64ZxePf7HLRKrofhHZrgwIkLeH3VPry+ah+2/3ksfvfBJizbcQwA0PGy+s7X9n1q\nOeKiwp3PT1+4hMR6UQDsXxpKwaXWXlJmw5HTF9EiMRY/HT7rXP78MnugPVdciu/3ncD0f20AANSP\nicDWWWOrfGwrd+Xj54LzuLZXChrUi3S+57miUlzWIMbj6076mFPvM3sZUpJi8fFvB7ksLyopwyeb\njnh8XanN5tP+fBGIlM4gALcA2Coijq+xP8Ee6BeIyG0ADgC4LgD7IiKdCa98i7RGcfhUCyg5czIr\nbKPPGLjftP1ia57L82HPrMK2x8fi0U+3Od+vVFcD7f/UcpwtKnU+33m00OX15y+VOR/nnLiAdC3g\n3/TmD1iz7wR2Pzkelz/8pctrhl7eGKt3F1Qod7dZSzCyYxPn80Ldfqti2rz1AIA/f74dmx8dgwb1\nIpH50jc4cOKC4d9pz7FCnL9Uhvm6lkInz19CclyUy3YbDpzElkNnMG1QGnJPXsCY51fj5Rt6oqCw\nGAWFxZj+bjaWbD+GtEZxWPn74ZUOCFjd4/JHIFrpfKuUEqVUd6VUuvbzhVLqhFJqpFKqvVJqlFLq\nZCAKTNah6kIDcj8opVBUUuay7L21B/C+FnD00zymZi3C7mOFSM1ahNSsRQCADQdOYcuhM85gD9hT\nLfmFRTh06gKe/Hw7bDblMh7+mYuuKRejP3GXx75yPs45ft4l5XC2GsFpxY5jKCopQ/7ZIqzZdwIA\nKgR7AIbB3mH5zsBkgvs8tQwAcODEBQDG8y6Mfn41rnn1O5wrLj/G2Yt2uGxTZlO49vXv8fj/tqOg\nsBhvr8nBxZIy3P5utnObJdvtVz8/Hz9fpVZHS7XtgyGgrXSIAkmp2tuM0FfbjpxB28bxuHipDD2f\nWAoA+ObBEWiZXA8A8PAn9nz9Tf1aV8jd/+2rXc7H/1l/EH9cuLXC+8/8aCt6tExERJhgw4FTGN+t\nGeKiy9MsT32xE5P7tsLFS2VomhCDOV/u9Fre4X9bhelD2/h0rC+t2IuXVuz16bXe7D5WiFbJ9RAT\nWX5ce/MLkRATiSYJMbDZFLI+2uLymkulrmmTRVvzcHWP5s5Khb6RyRdbjzofX7jk+gX35/9tcz7u\nM3tZpWXt/5fleDizk9dtikttmPXZNozu3BSD2jWq9D39YcqAP/HV73C8sBjfZV1Rpe0vldpgU8rl\nH4hCQx/ga1v9/sKlUoSJ+PR/snb/CWw8eBp/XbwTv+jVAjvyylMhmS99g+jIcKx/aJTLvha5pVuW\n6GqCRsHeYXPuaSTE2D/ad723Ae/e2tdlffdZSwDY0zWXyirPH89dvb/SbYLluSW78NKKvZjUswXG\ndrkMQ9o3Qlx0BEY9txoAMG9qH7y0Yg82Hjxd4bUf/XjI+dgRyNNmfoERHRpj3rS+FbYHKt5QXfjj\n4WqX+Um3qwQjb6/JwdtrcgxTTYFkyoC/ObfiyfZmzPNfI8dDXo9Cx177qh1V/DKbQudH7akOo/+T\n7/edgE0plxpa3pmLaNYgFgAwee5a5/KP3ILG2aJSoKjUJc3g2JevHKmXgsJiU3VmclwxfLzxMD7e\neBgzRrTFH8Z2dK6f9vZ6j6+9f8Fml+cXtfsNK3d5Tim5tzoK5g3WmmDKgF9dOVpej2qX2lTDn1fJ\n1Is3/MMe0OdN7YNG8dF4YdluLN+Zj+eu64EJPZpXaR/b885WvpEPPAXBSa8ZtwuvS/YcO4fck9X/\n/CoFDHl6pe59Cg23W7u//NZj/tkiFJXUXMDv3ya5xt7bgQGfaq3adM/W/bL8ja/3YVC7Rth//Dx+\n98FG53L34Hr/gs0VapaeBPt4jdIedc2S7ccqNC2tine/P4Dj58rb9a/9ufI2JV9tO1rpNtXx4yOj\n0Uu7jwMAB4NQ8TTN8MgbD57CKyv2hLoY5IcNB045W3MAcGld4mCzqQotTbwpLCrB377ahZIq5KqN\n7Dx6FutzXINB7skL+MuXO3HVy9+6BHt/3fefwL1XbdeuSbzHdV2aJ+DVG3tV+b02a52/qsP9auqR\nT7x3blv8Ux4e+XSb120AYOFvBnpc98mMQcjs1gwAcEPflkiOi8KoTuXNTmOiav4eomkC/qTX1uBv\nSwLXM46qpqikDKt3FyA1axGyc/xreXvt62uczeYA4xrv88t2o8fjS6o84NSzS3bjlZV78cnGym+2\nvbBsN1KzFuHMxRJsPHgK+WeLMO6Fb/Crv3/vst3ba3KqtO/q2ldwvvKNaqFfD6j+sA9L7huK/U9d\niT+O61hh3Z3D2iIhtvYkH/bmn3MOF1GZLs0T8ML16Ybr0lsm4tWbemHe1D547OouAIBXdF9sjeOj\n/S9sJUwT8Ck0shZuwa/fWgcAeP+Hg5VuX1hU4tK+vLoWbbG3XDmhC/hr9h1HatYiwzyso2ZfVGrD\npVKby2W8uxeW2a8Qb3pzLSa9tgZ9n1puuJ2vVwu+StJ6iQZLQ7eORg6emhc6evf+YWwHtGkcBwAY\n2Lahx/d/+truCAsThIUJfjO8Lb68d4hzXc6cTEzo0dw5tk6HpvVxVfdmPh3Hm7/O8Ol17r7Z4/mm\nbsvkWKzRtQaMjggz/LLqk1o+HciIjk2cLb30Lb6mDUoNQGm9s0TAX7b9GF5bFfj2wFV1+PRFdHzk\nS483huqy9TmnnI+/2XO80u0znlyG9D8vdVm2ale+szORXlVz2p9rXwL6POzh0xdx4lwxIrSu/GVl\nNjzz1U5kPLmsQlO7E+eKkf7nJc7n+m7+Rt79Prhjtg+9vLHHdZP7tAz4/jY8MtqZenB4/voeuH1I\nG2x4eBT2PXUlpg5Mda4b1ck+EO7Q9o2x4M4B+OCO/ph/R3+P739NzxYuzzs1S0DOnEyX1k99UpNx\nY79WeOfWvnjlxl5YcOeAah/HqM6GA/RW2+P/2+7yPK1RnPPxwrsGonliLO4a1haAvT3/iA5N8Pz1\nPXD3iHbO7Rw1em9G6HoV1xTTB3ylFG5/NxtPL94Fm03hy615QZ9V58uteSgqsWH+usprwA55Zy5W\nCEw/HT6Duav3Bbp4ftG3mz9+rrjS3rHFWgeYia9+hxxtMK5/eQigRjl8vbe/+xl9Zy9zfjGcuXAJ\nDyzYjPPFpRg0ZwV6P7nMOXbL+pxT+Mc39pY2+hp6UUkZej+5DKcv+DboV02aOjAVr93UC7MndXMu\na6MLNgAw59ruNbJv9/+9ST1TAAAN46MRHiaYNaE8gPVunYScOZnoltIAjeKjMUCr3X9139AK7/v2\ntD6Iiqg87ERFhOGpSd2cY970TUsOarPpbi0aeFynz9M3SbCXL2t8R2f5RASTeqbg7ivKA35XL++X\nNb4jnv1VD0RH1HwOv/YkymrIs7q8/n+yczHzo62YPalrlYebPXOhBIXFJUhJqudzGcRLd9HlO47h\nTx9vxeJ7hyJJu5Q+cvoiBs5ZgXuuaIcHxnRwbnvVy98CAKYPbetzWQJh3c8ncd0b3xuu23m0EAuy\nc/HAmA4oKbUhNircsKPS5tzTGP63VciZk+mxN+2Jc5dQL9n1X3S/c8RGhVnOmpc9ODnu4SzUdbBx\n3HDVd2LSD0frfkM22DK7N8OM4e1w5UvfVFhXXGrDlW417bAwQc6cTMMrou+yrkDT+tFo91DF4Quq\nypG2mTYoFYu3HcXXfxiOcA/DFC/8zUBsPHjKcB0ANE0oz0n//eZeGNfVt9SM3v/uHoyG8VH448It\nFa4oE+tFOr+4P7yr+lcEet7uISTHRWHb45UP4hYTGV6lLynH1UEwmK6G796FWv/hz9M6oOSfLcaF\nS6VIzVpUacueUc9/jcF/Xel1m8o4Pi5Gld/b3snGsbPFzjw4AORrw8B+7WWMkVApLbPh6cWeu+OP\nf/EbzPsuB10f+wo9n1iK69/4Hj8dPoPUrEU4cvqi4ft56lw15OmVWKI1hVNKebwy+2BdrsfyGKVn\nVu4scAaq9VVojhdIc2/p7Xw8b1ofPPurHujcPMFw27FdylMSn8ywj8DoGHtn1tWdXfLCANAiMRYR\n4a4fafcv00W/G4y2jV2vEvRuH2IfRqFfm4bImZOJ1g3jPFZ2erdOcm5vRF+WQAR7AOiW0gDNE2Px\nUGYn9GiZ6Pwb/bJ3CjY9OgYvTk7H4vuGICPV3qb9579cibemVj+XHxXuPTTGRUdUe2TS2qDulbgS\nnR5d7PLcJcg6x82Asybw3lrvaRZvY3AvWJ+L5omxGNze+/gXjg+dPt2Re/KCy6BJO4+WB6Yw5/bl\n73HghHELjuPninHo1EWkt0z0WgZv1v18EnNX78fcW3pXmHQiNWsRrklvjhcm9wRgH3kw+4DnWp27\nzYfOOG/mrjAYCKvdQ1+iQaznm5LPLd2Njbmn8foq11SWoyu9L2bMt7e4WHDngBoZ68Wb9Fbl52lE\nh/Kc7S96tsBHWkuiRvFRyH54tMvrkuu53kidOigNUwelAQCe/mV3l15qOXMy8eXWPPRv0xBREWEu\ng6E1iI3E8geGY9wLq7HzaCEen9AFj31WeXNDX8TVYDPDjpcl4FPtS7C4tMwZoCemu94fEBFc0bEp\n1mRdgYslZbDZFEY/7/1/J7FeJB69ugtW7lpVYV2rZN+v9GsD0wV8b5MJONYIxHlZX9VJdZRSFVIz\nDy60D9C0+g8j0Kqh538Ex6sKzhWjqKQMMZHhLr387O+v397+Cn3qYdgzq5yPP1h3EDdoswlNePlb\nHDlT5FN+s7i0DI9+sg3/ybbXkE9duISGuqZhZ7QvxU82HXEGfF9uWH6g3bt4x0NzRm/t6nceLaww\nBG+gGF1x+CohJqLSkSSv7tEckWFhzu31nrs+3RnwHTdBq+q6jIo3bsdrqSBPV0Vf3jsEJ8/bz3dN\nBXxvqcxAqkruu3mifYgL96GKWzes59IUWP85WnzfEIx7wTXVFhNZt5Midbv0VXBU1wvPET9Fyh8f\nMRhn5IutefjpsGtnjvd/OOic3cfdmBe+xpq9x/HAgs1ebwh/sfUobn17vWGzRMerbvnnD1iqTS7h\n6f7nW9+Wd/M3Kn9VLdue7wz2Rn791g8uz32dGMJhTyXjggfb4p8C13PyjiFtXIJBo/hoPPPL8huq\nOXMy8fINPZEUF4U/jO3gTNEY0d8Q9Zf7FZsjCIuIy5d7Tbl3ZHtcEYTWJ1Wl/3Msf2CYy43l129y\n7ezV8bIEPHKV61yzjrGR6irT1fC9cbT6EHifR/K379sv+fXf9o5hax3LPttcPgZ5UYkNN75pD473\nXNEOqVpLirHPr8auY4W4Uze87Jp9Jyo0SwTsVyZr95/AN3uOO29GOcqYc9w1neNLG6O9+YX497pc\nPJTZCSVlCpfKbNie5/oF5l4jc+/B6D52e123OIBd5e8c1hYfbzqM/QXnsenR0UisFwWbTeEPH27B\nHUPSXLadoWuup7d39niEiRjO5VpXh4n+v9GXh7oILvQ3oFMbxiE8TPD2tD743+Y851WR3m2D0/DE\n5+XNMj11qqorTB3w3fPer66054FFAKOKuPvEyd7GtvDUpV7/wdyltbt/o4rDy+pHVNQb/rdVLs+N\nmj4apZz0pry1HodPX8S0wWkY8cwqw2FxHe8786MtGNvlsgrrgznZcm2y84lx6PjIYq/bREWE4f3b\n+2HN3hPOaf3CwgT7nrqyymlD9xuunvZTXXcOa4M3vq49QxyHUpjuM+II/sM7NMHwDp6vQj767UCc\nOn8Jg9o1qvNDqJs6paPPe+uJiGHwem3VXtzyz/I0xtBnfG+dc8JLj86q2nm00GX2HQejLvi7j3lP\nlziCeUmpzeMY6Ae1UQc/WJeLqfMqjrDoKaVlBhPTXUe07N26vAWM0Yfc0Svy2l4pWPHAMAD2y/1r\ne6e4bBceJgHJZackxeKeK9rhrSl9/H4vqp5erZIwslPTOh/sAZPU8KubW1ZKYeSzX1dYHsixeErK\nAlMbnjZvneHykjIbInU1wspqkY4vDverBb1Jr63Bvqeu9Lj+N+9XbTyRuuixq7sgPjoCBYXFWLL9\nGAT2duY7PAxZnDW+Iyb3aYUOukm8a5KIuPTJIN84Knp1vbWNr0xRwz9YzfGwgzHIWqDSH/qhC/SU\ncu1zMPr51XjwQ+NheItKyqo8F6nZ55H1JCYyDLMndcOjV9tv0o3s1BS9Wyfh5v72DnrulfSo8LCg\nBftAqqO3AgLGkcqtaprNbEwR8IN57lKzFmGdl846w55ZhSlvrXMOIVBT/rhwS4UJoRdkl3cyW7Ur\nH88usc+B6q0vgbvZX1Q+HZsZOZr2pSTVw4aHR+GuYa4dij51a1UTrCaHFFiOZttGN8atwBQBP9g8\nDSvg8PXuAvzq72tqtAwfVzLc79R56/Hyir04eqYID/y3ahNwAMC873L8LFntoB9zRp+mGtDGeBRH\nfeuNhvHRFQJ6l+YNnH0fqO5yTloe4nKEiikCfm2sbB0/51+bdV8tWJ/rkup5ZeUebD9SM1PnOWx6\ndDR2PTnOZVnOnEw8d12PgO+rjW5YgJv6tcITE43brDtas8y/vZ9LMH/3tr7Y/ufKx0FxFx4m+Msv\nulW+IdVqbRvH4+b+rfDGLYEZOrmuMUfAt+z3dUUPLtziMnvPe2sPGrb0qS6j3P7DmZ2w7fGxSKwX\n5dLbsWsL+9gw3mY1qo6RHZvg8qb291r2f8Ocy2dP6oab+rXGncPaOMeMf//2fvjwrgGYpA3Bm6Z9\nQdzSvzUaxEYiMjwM9aJ8b6vw2+GhHbjOJ9a8LWMoLEzw5DXdAva/WdeYopVObazhh5K33rO+Gv+i\naxfz+Xf0w8C2rmMIvXpjL8yY/yNitODfPcXz+D6tkutVerN94yOjse3IWfRunYRYt3FZHJNihIUJ\nZo7vhA05p5B94BSS6kWhc/ME9G6dhBv6tUJCjP2L4IlruuKJa7o6X3/vyPZ4cXn1p8R8cFxHPGgw\nS1NdYfRZqcqwEGQOpgj4VPPcx7NxD/YAnGOXl+p6tW2ZNQbXv7EWD4y+HNvzzuLWwWnYdPA0urZI\nQGxUODo87LlDU1JclOHAdFtnjUGsW5vol2/siX+vy0WnZvaWMyLiDPZGru2V4gz4dw7zPOKjFVzV\noznmV2G2Mqr7GPApYKK1vHk9XW08ISbSOYWdYwYifRAf1akplu04huqobxDImzWIrVY3/pSkWEwd\nmIob+7XC5U3rXvPKQLp3ZHvM/+GgT714qW4xRcBnSqd26NI8AVnjO+IXblPYefPmlAzszT+HUc99\njTd/nYET54vxx4Vba7CUdmFuszZZWUwQZlqi2sEUAZ+Ca/4d/QyXi4hPs/e0axLvHJRudS2c9IXI\nLExxDWfRzqEBNc5gsLROzSrOxJTeMtEwf091A1u0WZspAn5d8fz1gW+XHgizru6MYR0aV1geFRFW\noQu6t9mpAoHf3SHA7wDLYMAPkvrREZjUM6XyDYOsfZN4TB2UZtiipaTUhn9OdR2dMTKc0aGuqeqX\nKM+s+THgB4v2aXo4s1Noy6EzbVAqlt5v78g0vmvFlE5JmQ3D2rvW/JsmxASlbBR8vLoyPwb8IHHU\nnvp7GMslFB67uryVir0Dk2uHosb1oysMMvVwpuuUb1T3sZWbdZgi4NeFm7bVHV1x6f8NRXpLe0/V\n2MhwbH5sjOF2vxvZHoB9gmx3nQ1uujo0qV9xPtNoXTvsKQNa4+UbelbYxr3Ha6BZdXjmYPH2b8i4\nb36mCPh1geODVtWUSFqjOOeIj09e0xUNYiNx57A2+P0Y185FV3Vvhpw5mRWG83331r74QuvwBACj\nOzdFlq4G/3uDyTRu6Fc+GmTv1OSgTHJNRMFT4+3wRWQcgBcBhAN4Uyk1J9D7UHUg++ioPTU2qFnf\nN6o9fjx4Gu/e2hepWYsA2OfedD+qmePt+f8jZ4qcXeEdFWL3inF8jOupvap7M0xMb4G7hrWtMFuW\ng34ANH0mp0PT+s75ecl8WLO3jhqt4YtIOIBXAYwH0BnADSJSJ5LAgWqNcqNWa54yMNXjNveNuhzv\n3trXZZmIbuxut6Log3vrhvap2hwpGk9XEhPTy3u/GgV7B0caKSKsfBvmeInMoaZr+H0B7FVK7QcA\nEfk3gIkAttfwfv22Z/aVztp2dXx+z2AUFBZj2tv2ScCHtm+MpyYZj6N+57A2eH+t8aBVoqvhVwy4\n9jWzJ3V1TqzcJCEGGx4ehQaxkTh5/hKa+Nia5jLtdfopGts2jsfOo4Xo0dLz6JeBUvuv1YjqrprO\n4bcAoB+r95C2zElEpotItohkFxT41q0+UPf5Vv1+OHY/Od7Zzd8XXVs0wIiOTTBlgH0u1IQYz9+p\nM8d3wk9o27/DAAAKDElEQVSPu07GMfTy8maQdwxpg/oxERjczrVpZMM4e23evRNUw/hoRISH+Rzs\ngfKBzdJ0M0Y50lATDW4MU93H6RqtI+Rj6Sil5gKYCwAZGRk+he5ATRjeqH604YiB6S0TsSn3dLXe\na+aVndCrdRIGtK3YDPN/dw/G8XPG88zOvaU3Tp63z5bVtUUDbJ1VcXame0a2Q6vkesjs1qzScqx4\nYFil2+jd1K8VxnRpiib1y780bhuchrX7T2BCehACPqv4AXf7kDTMXb0fAPP1VlfTNfzDAFrqnqdo\nywKqzOZblOjdOsn5eNn9wxAfbfz998mMQcjsXjG43qs1iTQSExmOiektDGtP3VLsVwGeXtc8MdZr\nuaMjwnFdn5ZVqpm1aRyPNo2rPruPiLgEewBomVwPi+8bikZstVMnuZ9Psq6aDvjrAbQXkTQRiQIw\nGcBngd5JqQ8Bf8/s8fjP9P7O55VNeWY0td2tg9Pw2d2DKtxwJT+wChp0/JNbR42mdJRSpSJyN4Cv\nYG+W+ZZSalug91OdlE54mKDMpry2VHH45sERzs5I4e6jiMHedNHbNH7kA6Z0QoapfPOr8Ry+UuoL\nAF/U5D76pVV9uIJ1fxqJC5fKqrRty+R6XteH8RNCJsJOzuZnip62+tp3Ze3nG8ZHuwTyu4a1Rd+0\nZJ/2qw/4jvbw5J+60ImuTjP4eLDeYh0hb6UTKBmtk5B94BQezuyMxz6retYoy23AME+MbljqPyiL\n7x2KS2W2Ku+XqLZh4Dc/U9Tw9ZLjomrkfRvFR2PjI6PRKN74/WOjwmt8chAiIn+YJuA7EgGXNYjB\nsvur1/a8qpLiojBGNxUgc/hkBpz20DpMk9JxjjsD4yaW6x4aiYtVvFnrzeMTuuC2wWm4UFxm2EmL\n/MMbh0Q1xzQB38FTpTtQnU8iw8PQthodmYhqE9bmrc10Ad/dk9d05TysdQhr+EQ1x/QB/+b+rUNd\nBKJajU1hrcM0SWj+y5oDzyNRzTFPwHdGCqZviKqDeX3rMF1Kx3HT9u8390JivZppk09UVxk1amBK\nxzpMF/AdxnWtfKx4IiIrMU9KJ9QFoIBQbKZDVGNME/Bv0VrjpDaMq2RLqs0c0zPeNjgtxCWxHuby\nzc80KZ1f9k7BL3unhLoY5Kf0lon48K4BSA/ChOlEVmOagE/mkZHq23DVVDnW4a3NNCkdIiLyjgGf\nyOKiI8IBAL8f2yHEJaGaxpQOkcWFhwly5mSGuhgUBKzhExFZBAM+EZFFMOATWYhwljZLY8AnIrII\nBnwiIotgwCcisggGfCIii2DAJ7IQ3rK1NgZ8IiKLYMAnIrIIBnwiIotgwCcisggGfCILSIixj5PI\njrbWxoBPRGQRDPhEFsAxdAhgwCcisgwGfCIii/Ar4IvIMyKyU0S2iMjHIpKoWzdTRPaKyC4RGet/\nUYmIyB/+1vCXAuiqlOoOYDeAmQAgIp0BTAbQBcA4AK+JSLif+yIiIj/4FfCVUkuUUqXa07UAUrTH\nEwH8WylVrJT6GcBeAH392RcR+U+pUJeAQimQOfxbAXypPW4BIFe37pC2rAIRmS4i2SKSXVBQEMDi\nEJEDG+kQAERUtoGILANwmcGqh5RSn2rbPASgFMD71S2AUmougLkAkJGRwfoHEVENqTTgK6VGeVsv\nIlMBXAVgpFLOC8bDAFrqNkvRlhERUYj420pnHIAHAUxQSl3QrfoMwGQRiRaRNADtAazzZ19EROSf\nSmv4lXgFQDSApVpPvrVKqbuUUttEZAGA7bCnemYopcr83BcR+Yk5U2vzK+Arpdp5WTcbwGx/3p+I\nAoP3bAlgT1siIstgwCcisggGfCIii2DAJyKyCAZ8IgtRHFvB0hjwiSyAE6AQwIBPZAms2RPAgE9E\nZBkM+EQWwJQOAQz4RESWwYBPZCHM5FsbAz4RkUUw4BNZCDP51saAT2QhTOlYGwM+kQWwZk8AAz4R\nkWUw4BMRWQQDPhGRRTDgExFZBAM+kYVwDDVrY8AnsgAOpUMAAz6RJcRGhQNg4Le6iFAXgIhq3r9u\n7YdFW/PQKD461EWhEGINn8gCUhvFYcaIdqEuBoUYAz4RkUUw4BMRWQQDPhGRRTDgExFZBAM+EZFF\nMOATEVkEAz4RkUUw4BMRWYSoWjSakogUADigPW0E4HgIixNMPFZz4rGaU2081tZKqcaVbVSrAr6e\niGQrpTJCXY5g4LGaE4/VnOrysTKlQ0RkEQz4REQWUZsD/txQFyCIeKzmxGM1pzp7rLU2h09ERIFV\nm2v4REQUQAz4REQWEfSALyLhIrJRRD7Xns8SkcMiskn7uVK37UwR2Ssiu0RkrG55bxHZqq17SaT2\nTdwmIjlaGTeJSLa2LFlElorIHu13km57sx2rWc9rooh8KCI7RWSHiAww8Xk1OlbTnVcR6aA7nk0i\nclZE7jPleVVKBfUHwP0A5gP4XHs+C8DvDbbrDGAzgGgAaQD2AQjX1q0D0B+AAPgSwPhgH0cVjjMH\nQCO3ZU8DyNIeZwH4q4mP1azn9R0At2uPowAkmvi8Gh2rKc+r7jjCARwF0NqM5zWoNXwRSQGQCeDN\nKmw+EcC/lVLFSqmfAewF0FdEmgFIUEqtVfa/8LsArqmxQgfWRNg/RNB+X6NbbrZj9aTOHquINAAw\nFMA/AUApdUkpdRomPK9ejtWTOnusbkYC2KeUOgATntdgp3ReAPAgAJvb8ntEZIuIvKW7bGoBIFe3\nzSFtWQvtsfvy2kYBWCYiG0RkurasqVIqT3t8FEBT7bEZjxUw33lNA1AAYJ7Y05JvikgczHlePR0r\nYL7zqjcZwAfaY9Od16AFfBG5CkC+UmqD26rXAbQBkA4gD8CzwSpTDRuslEoHMB7ADBEZql+p1QDM\n0ibW6FjNeF4jAPQC8LpSqieA87Bf6juZ6Lx6OlYznlcAgIhEAZgA4L/u68xyXoNZwx8EYIKI5AD4\nN4ArROQ9pdQxpVSZUsoG4B8A+mrbHwbQUvf6FG3ZYe2x+/JaRSl1WPudD+Bj2I/rmHbZB+13vra5\n6Y7VpOf1EIBDSqkftOcfwh4UzXheDY/VpOfVYTyAH5VSx7TnpjuvQQv4SqmZSqkUpVQq7JdNK5RS\nNzv+oJpJAH7SHn8GYLKIRItIGoD2ANZpl1hnRaS/dgf81wA+DdZxVIWIxIlIfcdjAGNgP67PAEzR\nNpuC8nKb7ljNeF6VUkcB5IpIB23RSADbYcLz6ulYzXhedW5AeToHMOF5DdWd8OEob6XzLwBbAWyB\n/Q/ZTLfdQ7DfAd8F3d1uABmw/6PtA/AKtB7DteUH9kvezdrPNgAPacsbAlgOYA+AZQCSTXyspjuv\nWhnTAWRrx/UJgCQznlcvx2rW8xoH4ASABrplpjuvHFqBiMgi2NOWiMgiGPCJiCyCAZ+IyCIY8ImI\nLIIBn4jIIhjwiYgsggGfiMgi/h9x7b9aqiY4iAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fd8adedf828>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"plt.plot(s1d.dispersion, s1d.flux)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"metadata": {
"ExecuteTime": {
"end_time": "2017-04-21T19:38:27.445564Z",
"start_time": "2017-04-21T19:38:27.424888Z"
}
},
"outputs": [
{
"data": {
"text/plain": [
"True"
]
},
"execution_count": 11,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"np.shares_memory(sp._data, s1d.data)"
]
}
],
"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.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment