Skip to content

Instantly share code, notes, and snippets.

@eteq
Created August 16, 2019 15:55
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save eteq/004e01c82928386eaf600655bbd35081 to your computer and use it in GitHub Desktop.
Save eteq/004e01c82928386eaf600655bbd35081 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 31,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"\n",
"%matplotlib inline\n",
"from matplotlib import pyplot as plt\n",
"\n",
"from astropy import units as u\n",
"from astropy.coordinates import SkyCoord"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Start with a basic plot of the plane using the Astropy transformation machinery:"
]
},
{
"cell_type": "code",
"execution_count": 27,
"metadata": {},
"outputs": [],
"source": [
"plane = SkyCoord(b=0*u.deg, l=np.linspace(0, 360,100)*u.deg, frame='galactic')"
]
},
{
"cell_type": "code",
"execution_count": 37,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.PathCollection at 0x117f55940>"
]
},
"execution_count": 37,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlMAAAFACAYAAACY6/lAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3Xd8VFX6P/DPmUmFJCQkmTSQviEESEIvioBLXxEUUVdUQGTdVVRAEXb1J7qroNjr6iKKX0V0lSYgsKAUFUVK6CAgNW3SQ/qU5/dHigHTM5k75fN+veYFuXPn3mdS7jz3nOeco0QERERERNQ4Oq0DICIiInJmTKaIiIiImoDJFBEREVETMJkiIiIiagImU0RERERNwGSKiIiIqAmYTBE1kVIqWimVWOWRp5R6RCm1UCmVVGX72CqvWaKU2quUul7L2ImIqOkU55kish2llB5AEoD+AKYByBeRF6/apyuAewE8BeBDEZls90CJiMhm2DJFZFs3ADgjIudr2UcPwApAACi7REVERM2GyRSRbd0O4NMqXz+olDqklFqmlAoCABE5CqAFgO8AvKNBjEREZEPs5iOyEaWUF4BkALEikqaUCgOQgbIWqH8CiBCR6VrGSEREtseWKSLbGQNgv4ikAYCIpImIRUSsAP4DoJ+m0RERUbNgMkVkO3egShefUiqiynMTARyxe0RERNTs2M1HZANKqRYALgLoKCK55dv+D0A8yrr5zgH4i4ikaBYkERE1CyZTRERERE3Abj4iIiKiJmAyRURERNQETKaIiIiImsBD6wCIyH0ppTwB+Jc//AC0RNlNXnUPfS3PKZTNKl/1Yalm29WPEgD5AC6XP4qEhaRE1EAsQCeielNKeeO3xMcfVyZC/l5eXq38/PyCvb29W+v1+iClVCAAf6vV6m+xWPxExFuv1+t0Op1ep9Ppvby84O/vLwEBARIQEKD8/Px0np6eSqfTQa/XV/5b/n/l4eEBnU5X+XXFcwBgtVrFYrGI1WoVs9ksVqsVFV9bLBaxWCwV+8BqtYrVakVhYaE1NzdX8vLy5PLly6qwsFBntVotVqvVYrFYrCJi1uv1BXq9/rJSKs9qteZZrdYck8mUVVhYmFVQUJCF3xKx/Gr+X8DkjMj1MZkicmNKKQ8AYQAiAER6eHhEBgcHd/Hy8uoIoI3JZArV6XSeOp1Or9frdd7e3vDz85OAgAC0atVKBQYG6lu1aqUPCgryDAwM9A4ICND7+/uj4uHn53fF/318fKCU8yxHaDabkZ+fj8uXL1c+rvrampOTU5qTk2PKzs425+bmWnJzc625ubnIz89HYWGhKk/sLBaLxaLT6XI9PDySzWbz+by8vFP5+fkXUDZrfjKAFBEp0PgtE1EjMJkickHVJElRwcHBnb28vDqKSJTZbDbo9Xpvb29vfUREhLRt21bfvn1773bt2rWMiopSERERiIyMhMFggIcHqwFsQUSQm5uL5ORkpKSkIDk5GRcvXiw+e/Zs0fnz501JSUnq8uXLYjKZTDqdLqeGpCsFQDKTLiLHwmSKGq28C2cpgO4om5hyOoCTAD4D0B5lE1VOFpFspZQOwIcAOgO4r3yxX2oEpVQIgA4AImpKknx8fPTh4eHStm1bXfv27X3at2/fMjIyUkVGRiIiIoJJkgOrSLoqEq6qSdeFCxdMly5dUvn5+dbS0lKzTqfL0ev1yRaL5Xxubu7pgoKC8yhLui4BOC8iFo3fDpFbYDJFjaaUWg5gl4gsLV/ktwWAvwPIEpHFSqn5AIJE5HGl1GgAnQB8AWARF/ytXXltUicA0a1ateoZEBDQx2w2R+t0ugCDwaCio6N/lyRFRkYiNDSUSZKbuDrpSklJwYULF4rPnTtXeOrUKfO5c+dQXFxc4uHhccFsNh9MS0v7WUROAjgpIllax0/kSphMUaMopQIAHETZ8ilSZftJAENFJKV8bbrtIhKtlBoH4BqUJVMviMg0TQJ3IKqseCgCQLSnp2eMwWDoo5TqYTabI/38/Dy7dOkiCQkJLXr06OEXHR2NLl26wM/PT+uwyYlYLBacP38eJ0+exLFjx0oTExMvHz582JKRkWG1Wq25er3+ZF5e3t68vLzDKGtVPiMipVrHTeRsmExRoyil4gG8B+AYgDgA+wA8DCBJRAKr7JctIkHlNTwfA2gHYKaIHNYgbE2Ur9v3BwDRwcHBCS1atOhlMpk6e3h4tIyKikJcXJxnz549A2JiYvTR0dGIjIx0qiJtck6XL1/GL7/8gpMnT+LQoUP5Bw8eLDx16pQqKCgo9fDwSLZarYeNRuMes9l8AmWJVhpHJhJVj8kUNYpSqg+AHwEMFpGflFKvAcgDMKu6ZEqrOO1JKRUOIM7X17dbcHBwXxGJtVqtoa1atfKIiYlR8fHxLbt37+4bHR2NTp06wcfHR+uQiX5HRJCUlISTJ0/i+PHj5sTExLzDhw+bk5OTYTab8z09PU8VFhbuz8zMTETZzdRx1maRu2MyRY1Snjj8KCLty7++DsB8lBWY/66bT7tIm0d58X2foKCg61q2bHmD2Wzu1KZNG92AAQO84+LiArp27aqio6MREhLCViZyGUVFRTh9+jROnjyJI0eOFP3888/5R44cEZPJlAVgd0pKyv8A/Iyy7kJ+uJDbYDJFjaaU2gVghoicVEotRNns1QCQWaUAvbWIzNMsSBso76ZL8PPzGxQUFDTCZDJ1Cw4O9rr22ms9r7/++sC+ffuiU6dOTJrIbWVnZ2Pfvn344YcfirZv3375l19+ERFJMZvNO41G47cAfhaRJK3jJGouTKao0crrppYC8ALwK4BpKFva43OUFZtfAHCrM40cKl/epIe3t3f/0NDQkWazOcHf379F//799UOHDg3s16+fLiYmhiPmiOqQkpKCvXv34rvvvsvfsWNH4fnz58XDw+PXgoKCb7Ozs3eiLMFymmsDUW2YTJHbKp/7Klqv1/cLDw8fYbVa+3l7e7dKSEjQDRs2LLB///4ecXFx8Pb21jpUIqcnIjh79iz27t2LnTt35u7ataskLS3N7OnpeSInJ2drfn7+9wD2i0i+1rESNRSTKXIL5dMQtAPQNzw8/AadTjdYp9MZunXrhmHDhrUaOHCgd69eveDv7691qERuw2Kx4OTJk9izZ4/s2LEj+4cffrDk5uYWe3p6HszIyNhSXFz8I4BDIlKidaxEtWEyRS5LKRWh1+v/GBERcbvZbO7VuXNn3bBhw/wGDRrUok+fPggJCdE6RCK6SmlpKQ4dOoQ9e/aYt2/fnrNnzx4xmUzJhYWFq3Nycr4CkCgiVq3jJKqKyRS5DKWUP4AhERERt4rIsKioKJ+JEyf6jx492jc+Ph56vV7rEImoEZKTk7F161br6tWrs3766SerXq8/lpGR8VlxcfFmETmrdXxETKbIaZUXi/cNDg4e7+XlNd7Pzy9k7Nix3n/6058CBg8eDF9fX61DJCIbExEcP34cmzdvLv3yyy9zT58+bfLw8PghJSXlM6vV+q2IZGodI7kfJlPkNMrrnrq2bNlyTGBg4K1KqU7XXXedbsKECa2HDx+u2G1H5H5MJhP27NmDDRs2FKxbt64wMzMzT0Q2paWlrQbwvYgUax0juT4mU+TQqtQ93WY2m3t3795df/PNNweNGjXKo2PHjlqHR0QOJi8vDzt27MC6detytm7daiotLU0qr7daD9ZbUTNhMkUOpaLuKTIycpLVah0eGRnpM2HCBP8xY8b4JiQksO6JiBokKSmpst5qz549Vr1efzQjI+Nz1luRLTGZIs0ppf4QFBR0p4+Pz61+fn4hY8aMqax7atGihdbhEZGLEBEcO3YMmzdvLl21alXOmTNnTDqd7rvk5OT3Ubb0lUnrGMk5MZkiTSilugUHB9/t6el5W5cuXVree++9rceOHasPDQ3VOjQichMmkwnff/89Pv7445yvv/7apJTanZSU9C6ArSJSqnV85DyYTJFdlBeP9wgNDb1Hp9PdGhsb6zN9+vTgG2+8URcQEKB1eETk5qxWK3788Ud8/PHHuevWrSsF8HN5YrWFRexUFyZT1GzKE6heYWFhU5VSE+Li4nymT5/eeuzYsTo/Pz+twyMiqpaIYO/evfj444/zVq9eXWK1Wg+mpKT822q1fi0ihVrHR46HyRTZVHkC1S8sLGw6gD/17dvXa9q0acGjR49WrH8iImcjIkhMTMQnn3yS/8UXXxSZzeajqamp71oslvVcR5AqMJmiJitfMHhQeHj4dABjBg0a5Dl16tTgESNGwMfHR+vwiIhsQkRw9OhRrFixouCzzz4rLCkp+cVoNL5rMpnWiUiu1vGRdphMUaMppboZDIYHlFI3Dx061Gvq1Kmthw8fDi8vL61DIyJqdsePH8fKlSsLP/nkk4KioqLjqampL5V3BXJUoJthMkUNopQK9vPzu9vPz+9vf/jDHwJnzZrV+sYbb9R5e3trHRoRkWYOHjyId999N3fNmjUlIrIxNTX1DQAHhB+yboHJVDNRSukB7AWQJCJ/Ukp9COB6ABVNwVNFJLG8i+xDAJ0B3CciR7WItzZKKS+dTjc2PDx8bsuWLaPvv//+gClTpngbDAatQyMicihmsxmbN2/GG2+8kXHw4MH84uLi93Nyct4XkRStY6Pmw2SqmSil5gDoAyCgSjK1XkS+uGq/0QA6AfgCwCIRmW73YKtRXkjeOyIiYhaA0bfccovPzJkzA3r06KF1aERETiE7OxsrV640vfXWWznZ2dnn09PTXzaZTGtEpEjr2Mi2mEw1A6VUGwDLATwLYE4dydQ4ANegLJl6QUSm2Tveq+Jp5e/vP61FixYP9e7d23/WrFkhf/zjH+Hh4aFlWERETu3UqVP4z3/+k79ixYoii8WyOTU1dbEj9kRQ4zCZagZKqS8ALALgD+DRKsnUQAAlALYBmC8iJUopDwAfA2gHYKaIHNYo5oTw8PAFXl5ew/72t7/5zZgxwyc4OFiLUIiIXJbFYsGmTZvk+eefzzh9+nRaVlbW4pKSki9EpETr2KjxmEzZmFLqTwDGisjflFJD8VsyFQEgFYAXgPcAnBGRZzQMFUopX19f3ztatWo1r1u3bsHz588PueGGG6DT6bQMi4jILZw/fx5vvPFG/ooVKwotFsvnRqPxJRE5p3Vc1HBMpmxMKbUIwF0AzAB8AAQAWCUiU6rsMxTlSZZGMXYJCwubp9frx0+dOrXl3/72t5ZRUVFahEJE5PZKS0uxatUqywsvvJCVlpb2a2pq6r/Kp1iwaB0b1Q+TqWZ0dcuUiKSUF3a/AqBYRObbMRadh4fHjaGhoU+1bdu27eOPP956/PjxOtZCERE5jmPHjuHll1/OXb9+faHJZFqalZX1qohkaR0X1Y7JVDO6Kpn6BkAoAAUgEcD99liKQCnl6+fnN71ly5aPjx071n/BggWBXbp0ae7TEmlqzYEkLNl8Esk5RYgM9MVjo6IxISGqzufq8zyRPRQUFODDDz8sXbJkSV5xcfHGtLS0p0XkV63jouoxmXJRSqmQ4ODgR728vKbdd999fg899FALFpSTK6kp6VlzIAkLVh1Gkem3HhJfTz0W3Vw2rUdNz9X12ornmWiRPVmtVqxbt866cOHCLKPReCglJeXvIvKT1nHRlZhMuRilVJfw8PCFvr6+IxcsWNDqrrvu8uT6eOSsGpMwLdl8Ekk5v5/GJyrQFwBqfO77+cMxePE3NT7/2KjoWhOt2uIlsoUff/wRTz31VObhw4dT0tPTnzCbzV+JiFXruIjJlMtQSg2KjIxcFBER0W3hwoXBY8eOVRyVR86ssQlTck4RqruqqfJ/a3ru7OJx6DB/Q43PRwb61pqI1dWqRWQrv/76K5577rnc9evX5xUUFDyfn5+/jBOBaovJlBNTSimdTjfSYDC81Lt37/Cnn346uHfv3lqHRdQgNbXm1NZKVFvCVFvSAzS+Zaq2c55dPK7W134/f3it75WoMbKysvDGG28Uvvvuu5dLS0v/nZmZ+ZKIXNY6LnfEpgsnpJRSnp6efwoPDz9+6623frpr167Y9evXM5Eip1PRmpNUnqgk5RRhwarDWHMgCcnVJCYAKhOR6lQkKL6e+iu2+3rq8dio6FqfA1Dr87WdsyKumuKt670SNUbr1q3x1FNPtThz5kzYP/7xj8ejoqJOh4SEPKuUaqV1bO6GyZQTKU+ibgoLC/tl8uTJ//fDDz9Er1y5Mqhz585ah0ZUqzUHkjB48TfoMH8DBi/+pjKBWLL55BXdYgBQZLJgyeaTjU6YJiREYdHNPRAV6AuFspahiq622p4DUOvzdSVidSVbtb3Xur5PRLXx9fXF7NmzfU6fPm1YuHDh3KioqFOhoaFLlFKBWsfmLtjN5wSUUsrDw2NscHDwKyNHjgz95z//GdiuXTutwyKql9pqiWZ/llhj19krt8U73Mi6uqZcqC3e2uqxzi4ex5orspnS0lIsXbq09LnnnsstLi5empmZuVhE8rSOy5UxmXJg5TVRI0JDQ18bNmxY2KJFi4Lat2+vdVhEDVJbLRFQew2Ts9UY1RZvXTVV9am5ImqIkpISvPfeeyXPP/98blFR0b+zsrKW2GN+Q3fEZMpBKaX6hIWFLR08ePA1L7zwQlCnTp20DomoVjUlErW1yNTV+uRKmtpyVXEMZ0ouyTEUFxfj7bffLnnxxRdzi4qKFuXk5LwpImat43IlrJlyMEqp9hERERsHDRq0acuWLXFffvklEylyeLUVV9dWS1RXDZMrqeu91lVzxQJ2aiwfHx/MmTPH++TJk4aZM2f+Myws7LS3t/eE8uXNyAbYMuUglFKBBoNhUatWrSa98cYbwaNGjeIvOTmNpk52SXW3XLEbkGwlKSkJc+fOzdm+ffv5tLS0GSKyV+uYnB1XudWYUsorKCjokaioqLnPPPNM0D333OOp1+vrfiGRndXWxVTbtAAV+7B7qnZ1fZ/qmnoBYDcg1U9UVBRWrlwZeOjQocC//vWvmyIiIvakpqb+TUTOaR2bs2LLlEaUUsrX1/e2gICAJffff3/refPmtWjRooXWYRFVi60m2qvre8zRgNRYW7ZskVmzZmXm5OR8aTQaF4hIttYxORvWTGlAKdXTYDAcmjx58jsHDx5ss3DhQiZS5NDqmiOprjmYqOnq+h7XZx4rouqMHDlSHTt2LGTx4sXTo6KiTgQGBj6olGJ+0AD8ZtmRUso/LCzsPz179vxmy5Yt3T/88MPAsLAwrcMiqlNdXUzuVEiulbq+x/XpBiSqiV6vx7Rp0zxPnDhhuPvuu581GAxHlFIJWsflLNjNZwflXXq3t2rV6uWnn3669YwZM7y4CDE5qurqbmpbWJjdeI6BawOSLR05cgR33313ZlJS0ldGo/EREcnVOiZHxk/0ZqaU6hIWFvbThAkT3j58+HD4zJkzmUiRw6pp+P2wrqHsxnNwdXUDcmoFaoju3btj3759wYsWLfpzRETEiZYtW07hVAo1Y8tUM1FK+YaEhPyzdevWUz/88MPggQMHah0SUZ3qmuKArRqOrSkzsBPVJCsrC4888kjuli1bTqelpd0pIizEuwqTqWbg5eX1x6CgoA8ee+yxkEceecTHw4MzUJBzqM8s3OScOMM6NdVPP/2EqVOnZmZlZX1UPuqvROuYHAX7m2xIKeUXHh7+8cCBAz/fv39/m0cffZSJFDmkNQeSMHjxN+gwfwMGL/6msqunrlm4yXlxhnVqqv79++PIkSPBs2fPvt9gMBxVSsVrHZOjYDJlIx4eHkPCwsKOP/PMM7du3749KCqKd3PkmGr70OQUB66LUyuQLej1esyfP993+/btnbp3774tNDT0eaWUp9ZxaY3JVBMppXzDwsKW9u3bd/WePXvazJw504s1euTIavvQ5BQHrotTK5AtxcTE4MCBA60ffPDBB8unUYjVOiYtsWaqCZRSAwwGw2dPPPGE4cEHH/RhEkXOgHVRVB0WqFNjHTp0CLfddltmRkbGWxkZGc+IiKXuV7kWtkw1glLK22AwvNmrV68Nu3fvvmbWrFlMpMjhsC6KGqI+Xbw1/U6Re+vZsycOHToUPGPGjDkGg+GgUsrt6gIcvmVKKeUDYCcAb5QtzPyFiDyllOoAYCWA1gD2A7hLREqVUn4AVgDwAzBFRJJtHE8Xg8Gw4bHHHmszZ84cX84ZRY6otnXaAHANN6pWbaP5uPYf1cf+/fsxefLkjIyMjCdycnLe1Toee3GGZEoBaCki+eVFbt8BeBjAHACrRGSlUurfAA6KyDtKqfsBXARwAcCdIjLfVrH4+/vfGRIS8uqqVatCEhI4yz45rvosissh8NQQ7Aak+iooKMC0adNydu7cuaN8XqoCrWNqbg4/bl/Ksr388i89yx8CYDiAP5dvXw5gIYB3AOgBWMsfNul7U0r5GAyG94cMGTL2008/DQwICLDFYYmaTX3W0mPyRA3BAnWqr5YtW+Lzzz8PfP/998f8/e9/P6yUulFEjmodV3Ny+GQKAJRSegD7AHQG8BaAMwByRMRcvsslABWfDJ8A+BSAD4C7bHDuLgaDYcOTTz7Z9oEHHmBtFDmUmlqYIgN9q21FYF0UNVZ9fqfY4klV3XvvvV4DBgzocNNNN+0IDAz8hyt3+zlFwY+IWEQkHkAbAP0AxFS3W/m+OSIyRkSGicilppzX399/Svv27X/YtGlTF47WI0fD+aLInrj2HzVGbGwsDh48GDxy5MjFYWFha5RSLbWOqTk4RTJVQURyAGwHMABAoFKqomWtDQCbFZorpbzDw8NXDBky5I2DBw+yPoocEueLInuq63eKk35STSq6/Z577rkx5XNSddM6JltzhgL0UAAmEclRSvkC2ALgeQD3APiySgH6IRF52wbnM4SGhm5dsGBB50ceecSXrVHkqDhfFDkS/j5SfRw9ehTjxo3LMBqN0woLC9drHY+tOEPLVASAb5VShwD8DOB/IrIewOMA5iilTgMIBvB+U0+klOoRFha2/+OPP46dPXs2EylyaJwvihxJfdb+4xxVFBsbi3379oV07959eUhIyN+Vi3zQOnzLlL34+fndZDAYlm7cuDGka9euWodDdIXqCnsBzhdFjoNzm1FDlJaWYvr06bn/+9//NhuNxrtEpFTrmJrC7ZMppZQKCQn5R8eOHed8/fXXQa1bt9Y6JKIr1PUhxdFT5ChqGs3HOaqoOiKCF198sXDJkiUn09PTR4hIptYxNZZbJ1NKKS+DwbBi5MiRf1y2bFkrT0+3X/iaHBA/iMjZsZ6KarNhwwbL9OnTk41G4ygROa51PI3hDDVTzUIpFRQaGvrT448/Pvajjz5iIkUOi5MlkrNjfR/VZty4cfpvvvmm7TXXXLPTx8dnlNbxNIZTTNppa0qp8NDQ0F1Lly5tP378eLf8HpBjqq6bhBNwkrN7bFR0tV3VVeeoYne1e4uNjcXevXtDhg0btiIgIOBveXl5n2kdU0O4XTefUqp9WFjYjpUrV7YZOnSo27bMkeOpqTbqlt5R+HJfEot3yanVlDBxAWWqKjc3FzfccEP2mTNnFmRnZzvNjOlulUwppWLCw8O3ffXVVxF9+vTROhyiK9RWG/XYqGjeuZNLYk0gXa2oqAhjxozJPnr06JL09PRFWsdTH27TxaWU6h0VFbVxy5Ythm7dXG7yVXIBtdVGcWFiclWsCaSr+fr6YsuWLUGTJk2aZzAYWqenp88TB2/5cYtuLm9v7+vbtWu3adeuXUykyGGxSJfcEX/vqTpeXl5YvXp14OjRo2eGhoa+r5Ry6HzFoYOzhRYtWoxt167dl7t37w7p0KGD1uEQAah+NmguTkzuqD4LKHPmdPek1+uxfPnygNtvv32SwWD4wpETKpeumfLy8hravn37L3fv3t06ODhY63CIAHASTqKrsTidaiMimDdvXv5HH320zmg0TnHELj+XTaaUUv3at2+/8YcffgiOiIjQOhyiSiy4Jaof/q1QBRHB/ffff3n16tUfp6enP+BoCZVLFqArpXq0adNm/fbt25lIkcNhwS1R/fBvhSoopfDOO+/45+bm/vmbb77JAfB3rWOqymH7HxtLKdUlMjJyy7Zt20LbtWundThEv8OCW6L64d8KVaXT6fDJJ5+06tu3719DQkIe1zqeqlwqmVJKtQ0PD/9248aN4X/4wx+0DoeIheZETcDidLqaXq/H6tWrA2NjY+cHBwffr3U8FVymZkopZTAYDHvWrl3bbsCAAVqHQ8RCcyIbYHE6VaeoqAhDhgzJPnny5AN5eXmfah2PSyRTSimf0NDQnz/++ONuI0eOdKnWNnJeLJ4laj78+6K8vDz069cv8+TJk+NE5CctY3H6xEOVNUmtfOKJJzozkSJHwuJZoubDvy8KCAjA5s2bgyMiItYopdpoGYvTJx8hISH/b+zYscMeeughH61jIaqKxbNEzYd/XwQA7dq1w+effx4WGhq6VSnVQqs4nDqZ8vPzu6lTp04Pv/feewFax0LujYXmRPbFvy+qcO2116pFixZ1MBgMa5RSSosYnLZmSinVvX379t/u3bs3hLObk5ZYaE6kjdqK0/l3534efPDBy59//vm/jUbjPHuf2ymTKaVUSFhY2IFvv/22TUxMjNbhkJtjISyR4+AoP/dlsVgwfPjw7AMHDth9hJ/TdfMppXQGg+Gr999/P4KJFDkCFsISOY4lm09ekUgBQJHJgiWbT2oUEdmLXq/HV199FRQaGvq6Usqu/b1Ol0y1bt368VtvvTV23Lhx+rr3Jmp+LIQlchy8uXFvAQEB+OKLL0IMBsM6pZSXvc7rVMmUUireYDDMfemll/y1joXc19XF5sO6hrIQlshB8OaGEhISMGfOnLYGg+E1e53TaWqmlFItDAbD0R07drTv2rWr1uGQm6qpHuOW3lH49kQ6C16JNMYBIQQAVqsV1113XfbPP/98W2lp6f+a+3xOk0yFhYV98vTTT99y//33e2sdC7kvFpsTOb7qRvMBYGG6m0lJSUHv3r1TUlJSeopIRnOeyymSqZYtW9587bXXLt20aVOQRlNIEAEAOszfgOr+YhSAs4vH2TscIqon3gi5p6+++soyY8aM74xG4zBpxoTH4WumlFL8lUvMAAAgAElEQVRhQUFB76xYsYKJFGmO9RhEzomF6e7pxhtv1P/pT3+KDwoKeqA5z9PgZEoptUwpZVRKHamybaFSKkkplVj+GFvluSVKqb1KqesbE2B4ePh/3n777WBOzEmOgLMuEzkn3gi5r9dff72Vn5/fk0opQ3OdozEtUx8CGF3N9ldEJL78sREAlFIVleJDADQ4K/Tw8Li+a9eug8ePH89pEEgTV4/cA4BFN/dAVKAvFMq6CFhzQeT4eCPkvlq2bInXXnstOCws7L3mOodHQ18gIjuVUu3rubsegBWAoKyspN6UUp4Gg+HDZcuWtW5YhES2cfWooKScIixYdRiLbu7BGgsiJ1Nxw1NdYfrgxd9whJ+LmzBhgv7VV1+9Til1nYjssvXxG5xM1eJBpdTdAPYCmCsi2SJytHwV5+8APNaQgwUHBy/461//aujQoYMNQySqv9pmUubFlsj5TEiIuuJvt6Ybpop9yXUopfDBBx+0HjRo0HKlVLSImGx5fFsVoL8DoBOAeAApAF6qeEJEZolIbxH5pr4HU0q1DQgImLVgwYIWNoqPqMFYsErk2rj0jHvp2LEj7rvvPkNwcPDjtj62TZIpEUkTEYuIWAH8B0C/phwvLCxs+bvvvtva25tTSpF2WLBK5Np4w+R+/vGPf7T08/N7WCnVxpbHtUkypZSKqPLlRABHatq3Hse6Li4uLm7EiBEOP20DuZ6qBecFJWZ46q8s9WPBKpHr4A2T+/Hx8cFbb73VOjw8/E1bHrfBk3YqpT4FMBRACIA0AE+Vfx2PskLzcwD+IiIpjQkoIiLiwLZt2+K7devWmJcTNVp1y1B46hT8fDyQU2hicSqRi+HyUO5JRBAXF5dx+PDhgSJy2hbHdKgZ0JVS148bN27V+vXrOYKP7I4zJBO5n6uXnhnWNRRf7kvisjMubsuWLTJ16tSvkpOTb7LF8RwqmYqIiEjctm1bHFulSAtcKoaIeFPlHkQE8fHxGYcOHRokIqeaejyHqUtSSg3t06dPWyZSpBXWTxARi9Ldg1IKL774YnBERMTLtjieQyRTSikVHh7+2pIlS9i9R5rhDMlExJsq9/HHP/5RhYWFDVBKNfki7xDJFIDr+/fv36Zr165170lkI1wqhoiuxpsq96GUwksvvRQSERHxUt1713EsR6iZioqK2rFhw4Yh8fHxWodCbqKmUTxMnojo6qL0ikTq6m28Vjg/EUHXrl3Tf/nll+4iYmzscTRPppRSEd26dTt49OjRUE0DIbfCIlMiqi/efLm2t99+2/TEE0/8Kysr65nGHkPzbr7g4OC/zZ49O0jrOMi9sMiUiOqLy864trvuusvTy8vrPqVUo3MiTZMppZTOw8Nj+h133GHLBZeJ6sQiUyKqL958uTZ/f3/ccMMNLQEMa+wxtG6ZGjF69OgWLVu21DgMcjcsMiWi+uLNl+ubPXt2UGRk5ILGvl7TZCoyMvIfs2fPDtQyBnIfVUfvLdl8Erf0juLIPSKqE2++XF+fPn3g5+fXUyllaMzrNStAV0qFdu3a9ejx48dZeE7NjgWkRNQU1S07w/X7XMtbb71leuKJJxZmZ2c/19DXatYy5eXldeOUKVP8tTo/uRcWkBJRU0xIiML384fj7OJxeGxUNL7cl4SknCIIgKScIixYdRhrDiRpHSY1waRJkzx9fX3vbMxrNUumDAbDPRMnTvTR6vzkXlhASkS2wpsz1xQWFoZWrVqFKqUavBqLJsmUUsoLQExMTIwWpyc3xAJSIrIV3py5rsmTJ7fU6/VjGvo6rVqmrh81apSnUkqj05M7qFpwXlBihqf+yt83FpASUWPw5sx1TZw4sUV4ePjUhr5Ok2QqMjJyym233cZRfNRsKgrOK2oacopMgABBLTw5eo+ImoSj+1xXXFwcRCROKeXZkNfZfTSfUkqFhYVdOn/+fKS3t7ddz03ug8vFEFFz4ug+13X33Xdn/d///d8tIrK9vq/RYubxjtHR0Z5MpKg5saaBiJrThISoymTp6qlXKkb3VexHzmXSpEmtt2zZcguA7fV9jRbdfH2GDRvmp8F5yY2wpoGI7IWj+1xL37594eHhcW1DXmP3ZCoiImL4gAED+IlGzYo1DURkL2wJdy0REREQkUjVgFFydu/mU0oN7N27t71PS26iah1DK19P+HjqkFNoYg0DETWbyEDfams02RLuvDp06IDk5ORrAJyvz/52bZlSZcJCQ+23gsyJEycwcOBAeHt748UXX7ziufbt26NHjx6Ij49Hnz59KrdnZWVhxIgR6NKlC0aMGIHs7GwAgNVqxd13341Bgwbh6NGjdov14sWLGDZsGGJiYhAbG4vXXntN81hrkp2djYkTJ6Jnz57o168fjhw5Uvncpk2bEB0djc6dO2Px4sWV248ePYqBAwfinnvugdVqbfS5qxvBV2yy4pXb4vH9/OENSqS2b9+O+Ph4xMbG4vrrr7fbe7C1tWvXomfPnpW/4999913lc8uXL0eXLl3QpUsXLF++vHL79u3b0adPH8ybN0+LkOulup+Do/4MqlNcXIx+/fohLi4OsbGxeOqppwAAZ8+eRf/+/dGlSxfcdtttKC0tBQDk5+dj/PjxGD58OJKTk7UM/XdycnIwadIkdO3aFTExMdi9e7fDXZfq8tprr6F79+6IjY3Fq6++CqBh11Z7toRPnz4dBoMB3bt3r9z22GOPoWvXrujZsycmTpyInJycyucWLVqEzp07Izo6Gps3b67cvnLlSvTq1avy/dojzoULFyIqKgrx8fGIj4/Hxo0bNY+zJkOGDGkJoP4tPyJitweAjkOHDk0TO0pLS5M9e/bI3//+d1myZMkVz7Vr107S09N/95rHHntMFi1aJCIiixYtknnz5omIyNdffy1vvvmmpKamyrRp0+wWa3Jysuzbt09ERPLy8qRLly5y9OhRTWOtyaOPPioLFy4UEZHjx4/L8OHDRUTEbDZLx44d5cyZM1JSUiI9e/asfA/Tp08Xo9Eor7/+unz99deNPvegRduk3ePrf/cYtGhbg46TnZ0tMTExcv78eREp+7nY6z3Y2uXLl8VqtYqIyMGDByU6OlpERDIzM6VDhw6SmZkpWVlZ0qFDB8nKyhIRkcmTJ0thYaHMmTNHjh8/rlnsNanp5+CoP4PqWK1WuXz5soiIlJaWSr9+/WT37t1y6623yqeffioiIn/5y1/k7bffFhGRd955R9avXy+HDh2Sxx9/XLO4q3P33XfLf/7zHxERKSkpkezsbIe7LtXm8OHDEhsbKwUFBWIymeSGG26QX375pcHvYfX+SzJo0TZpX37N+cfqQ1d8vXr/JZvEu2PHDtm3b5/ExsZWbtu8ebOYTCYREZk3b15lrEePHpWePXtKcXGx/Prrr9KxY0cxm80iInLTTTeJ2WyW2267rfJ30Zaqi/Opp5763eew1nHWZN26dRIWFva61DO/sXfNVO8hQ4bYtfjcYDCgb9++8PSs/5QRa9euxT333AMAuOeee7BmzRoAgMVigU6ng06nq0gO7RJrREQEevXqBQDw9/dHTEwMkpKSNI21JseOHcMNN9wAAOjatSvOnTuHtLQ07NmzB507d0bHjh3h5eWF22+/HWvXrq2MVSnV5FhtVbewYsUK3HzzzbjmmmsAlP1cANjlPdian58fKrr9CwoKKv+/efNmjBgxAq1bt0ZQUBBGjBiBTZs2ASi783bE91Khpp+Do/4MqqOUgp9f2aXQZDLBZDJBKYVvvvkGkyZNAuAYf891ycvLw86dO3HvvfcCALy8vBAYGOhw16XaHD9+HAMGDECLFi3g4eGB66+/HqtXr27we7DX2n1DhgxB69ZXrnYycuRIeHiUVe0MGDAAly5dAlD2+XD77bfD29sbHTp0QOfOnbFnzx4AqIxdKdUsP4vq4qyJlnHWpHfv3vDw8Bhc3/3tmkyFhob2jY+Pb2HPc9ZGKYWRI0eid+/eeO+99yq3p6WlISIiAkBZImM0GgEAo0aNwo4dOzB+/HjMmTNHk5jPnTuHAwcOoH///g4Za1xcHFatWgWg7EPv/PnzuHTpEpKSktC2bdvK/dq0aVOZED788MMYN24cdu/ejZEjRzb63LYawffLL78gOzsbQ4cORe/evfHRRx8BgF3eQ3NYvXo1unbtinHjxmHZsmUAan8vM2bMwKBBg2C1WuGISz7VFLsj/wyqY7FYEB8fD4PBgBEjRqBTp04IDAys/FCs+jO588478frrr+PBBx/ErFmztAz7Cr/++itCQ0Mxbdo0JCQkYMaMGSgoKHC461Jtunfvjp07dyIzMxOFhYXYuHEjLl682KT3oOXovmXLlmHMmLLVUGr7O7/55pvRp08f9OnTB/7+/s0eV4U333wTPXv2xPTp0yu7Th0xzoiICFit1oj67m/XAnQfH59OUVGOUwD8/fffIzIyEkajESNGjEDXrl0xZMiQGvf38PDAypUr7RjhlfLz83HLLbfg1VdfRUBAQK37ahXr/Pnz8fDDDyM+Ph49evRAQkICPDw8qr2jqGglSUhIwE8//dTkcz82KvqKuV6AxtUtmM1m7Nu3D9u2bUNRUREGDhyIAQMG2OU9NIeJEydi4sSJ2LlzJ5588kls3bq11vcyatQojBo1yt5h1ltNsTvyz6A6er0eiYmJyMnJwcSJE3H8+PHf7VPxMwkMDMTXX39t7xDrZDabsX//frzxxhvo378/Hn744StqCa+m9TW0OjExMXj88ccxYsQI+Pn5IS4urjKhrU593oNWo/ueffZZeHh44M477wRQ898KUNbaVtHyZi9//etf8eSTT0IphSeffBJz587FsmXLHC7OivN7eXl5KKV0IlJnEaZdW6ZEpE1kZGSzn+ett96qLHCrrVizIhaDwYCJEydWNiuGhYUhJSUFAJCSklLZzaNlrCaTCbfccgvuvPNO3HzzzZXb7RlrTaq+h/z8fHzwwQdITEzERx99hPT0dHTo0AFt2rTBxYsXK19z6dIl2Op3oWINvtmfJcLbQ9eoJWOqvofIyEiMHj0aLVu2REhICIYMGYKDBw8263uwpZp+p4YMGYIzZ84gIyPDad5LdZw59uoEBgZi6NCh+PHHH5GTkwOz2QzAOd5XmzZt0KZNm8qW8kmTJmH//v0OcV1qiHvvvRf79+/Hzp070bp1a3Tp0qVJ70GLee6WL1+O9evX45NPPqlMRBztbyUsLAx6vR46nQ733Xdf5Weuo8VZISQkRAAE12dfuyZTZrPZEBYW1uzneeCBB5CYmIjExMQafyAFBQW4fPly5f+3bNlSOepg/PjxlaObli9fjptuuknTWEUE9957L2JiYn7XrGzPWGtS9T20aNGicgTS0qVLMWTIEAQEBKBv3744deoUzp49i9LSUqxcuRLjx49v8rltNYKv6nuYOHEidu3aBbPZjMLCQvz000+IiYlptvdga1XfS2FhYeVd3/79+1FaWorg4GCMGjUKW7ZsQXZ2NrKzs7FlyxaHbo2qyll+DrVJT0+vHHFVVFSErVu3IiYmBsOGDcMXX3wBQLu/54YIDw9H27ZtcfJkWffVtm3b0K1bN4e4LjVERRfehQsXsGrVKtxxxx1Neg/2nudu06ZNeP7557Fu3Tq0aPFbJc348eOxcuVKlJSU4OzZszh16hT69evXLDHUR0VyCpSVH1T9zHWkOCu0adNGAahfV199K9Vt8YiMjEyxTZ19/aWkpEhUVJT4+/tLq1atJCoqSnJzc+XMmTPSs2dP6dmzp3Tr1k3+9a9/Vb4mIyNDhg8fLp07d5bhw4dLZmamprHu2rVLAEiPHj0kLi5O4uLiZMOGDZrGWpMffvhBOnfuLNHR0TJx4sTKEWIiIhs2bJAuXbpIx44dr/h+N4WtRvBd7YUXXpCYmBiJjY2VV155pXJ7c7yH5rR48WLp1q2bxMXFyYABA2TXrl2Vz73//vvSqVMn6dSpkyxbtkzDKBvO2X4OVzt48KDEx8dLjx49JDY2Vp5++mkRETlz5oz07dtXOnXqJJMmTZLi4mKNI63bgQMHpHfv3tKjRw+56aabJCsry+GuS3W59tprJSYmRnr27Clbt24VkaZfW6uO7otbuFnin95sk5F9t99+u4SHh4uHh4dERUXJ0qVLpVOnTtKmTZvKz4e//OUvlfv/61//ko4dO8of/vAH2bhxY6PPa4s4p0yZIt27d5cePXrIjTfeKMnJyZrHWZuHH344F8AoqUd+Y7eFjpVS+vbt26eePXs2xC4nJLfQYf4GVPcbrACcXTzO3uEQEV3h6nX7gLJWqvqWIJB2Xn31VfPs2bNnisgHde1br24+pdRopdRJpdRppdT88m2xSqndSqnlSqn6HCc0LCzMsWfSI6fDNfiIyJFx3T7nFRUV5REcHNyxPvvWmQQppfQA3gIwBkA3AHcopboBmANgPIC9AOozDjmwvnNOENUX1+AjIkfGdfucV+vWreHt7R1en33rMzVCPwCnReRXAFBKrQRwEwA9AAFgRVmvSl08vby86r1oIFFtuAYfETkDrtvnvDw8PKDT6eo143d9uueiAFys8vWl8m2vAdgAYCCALfWJi8kU2YIt1+AjImpObD13XuWrkdgsmaouARIROSAi/UVkiohYqtnnah61TYRGVF+sQSAiZzEhIQqLbu6BqPKWKL1SldcrWywvQ81Hr9cDNkymLgFoW+XrNgAas2y5pWIyOqKmYA0CETmTCQlRlS1UlvIR9LZcr4+ah8ViAYB6JS71SaZ+BtBFKdVBKeUF4HYA6xoRl8lsNjvGypbk1DiCj4icDVvUnU95A5CpPvvWmUyJiBnAgwA2AzgO4HMROdqYuEymesVEVCvWIBCRs2GLuvMxm80QkdL67FuvIiYR2QhgY5OiAi5nZ2ezZYoajSP4iMhZcVSf88nNzUVpaWlWffa1Z0V4WmpqKkfzUaNcPYtwTpEJvp56vHJbPJMoInJ4j42KrnYmdLaoO67k5GRrRkbG6frsa7eFjkXEVFxcXJ9Rf0S/w3oDInJmVUf1KQBRgb5cUsbBnTt3rkBE6jXgzq5zFZjN5lKLxVIx3JCo3lhvQETObkJCVGXyVFG2MPuzRJYqOKhz586VAEipz752a5kCAA8Pj4z09HR7npJcBEfwEZGruHriYU6T4JguXrxohSMmUzqd7lJycmOmqCJ3xxF8ROQqWLbgHIxGowJgrM++dk2mSktLzzGZooaqaA4vMlmgV2VjGFhvQETOimULzqGkpMRSzxVe7JtMZWRk7Dt69GiJPc9Jzq1qczgAWEQqW6SYSBGRM2LZguPLzMyEiNRrWgTAzsmUxWLZt2PHjjx7npOcG5vDicjVsGzB8e3fvx9KqR/ru79dkykAJ44cOWLnU5IzY3M4EbmaimkSAn1/W0PXx9PeH8dUm59++qk4JSVlW333t+tPT0QsZrM5Kzc3156nJSfG5nAiclUlZmvl/7MLTRzR50B27NhxWUT21Xd/u6fCOp3u5wMHDtj7tOSk2BxORK6IJQyO7cSJEwBwqr772z2ZSk1N3bZnzx4WoVOdOIqPiFwVSxgcV1ZWFqxWa4aIWOveu4zdkymLxbJ3+/btLEKnWnEUHxG5MpYwOK59+/YBQL2LzwENkikAxw4cOCAWC5fpo5qxCZyIXBlLGBzXunXr8pKTk1c15DV2T6ZExKrX63fv2bPH3qcmJ8ImcCJyZVUXPgYAvVKVN4wsQtfW+vXrSwDUeyQfoE3LFJKSkpZ/8cUX+Vqcm5wDm8CJyNVNSIiqbKGyiADgOn1aO336NEwm068i0qA7d60mtvjfmjVr2MRANWITOBG5A5Y0OJa1a9eWZGdnf9TQ12mSTIlIfnFx8cXz589rcXpyEt4ev/16BrXw5Cg+InI5LGlwLCtWrMgrLCxc29DXaTblam5u7sdr1641aXV+clwVI/lyin779Sg21XuEKhGR02BJg+PIzc1FcnLyZRFpcB+rZslUQUHBqo8++ihHq/OT42KzNxG5C5Y0OI4NGzZYLRZLg0bxVdAsmRKR88nJyclnz57VKgRyUGz2JiJ3MSEhCrf0jqqcmFivFG7pHcWSBg289NJLmenp6e805rWarqyYkZGx+M033+SoProCm72JyF2sOZCEL/clVY7ms4jgy31JHM1nZ6dOnUJKSsoFEfm1Ma/XNJkymUyrP/3000KTiaVT9Bs2exORu2BZg2N444038tLT0xc39vWaJlMiUmK1WletW7eO1cV0BY7kIyJ3wLIG7ZWUlOC///1vkdlsbvAovgqaJlMAkJaW9vKSJUsytY6DHANH8hGRO2FZg/ZWr15ttVgsn4tIo7vJNE+mROTUhQsXWIhOANjkTUTuhWUN2nvxxRez0tPTX2nKMTRPpgAgMzPz2UWLFuVqHQdpj03eROROqq7RpwBEBfqyrMGOEhMTkZycfFpEmtSi4xDJVGlp6Zfr1q3LSk5O1joU0hibvInI3VSs0RcZ6IvknCIudmxHc+fOzUxJSXmkqcdxiGRKRKw5OTnzn3rqKbZOublhXUOhrtrGJm8icmUVtaJJOUUQcLFjezlw4ACOHTt2WkR+auqxHCKZAoCSkpIvvvrqq6ykJP7yuKuK+VakyjYFcAI7InJprBXVxty5czNTU1MfssWxHCaZEhFrbm7u42ydcl/VXVAEwLcn0rUJiIjIDlgran/79+/HiRMnTonIHlscr0HJlFKqq1Jqt1KqRCn16FXPnVNKHVZKJSql9lbZHqmU+kYptVYp5Vfb8YuLi79cv3595qVLlxr2Lsgl8IJCRO6ItaL2V14rZZNWKaDhLVNZAB4C8GINzw8TkXgR6VNl20MAZgFYCmBKbQcXEWt2dva8J554gq1TbogXFCJyR5wewb5+/vlnnDhx4hcR+dlWx2xQMiUixvKTN2RiKz0Aa/nj6tri3yktLV319ddfX0xMTGxIaOQCeEEhInfl4/nbx3GgL1d9aC5WqxXTp0/PTE1NnWnL49qyZkoAbFFK7VNKVQ3yTQDvArgfwMd1HkREjEbjnffcc0+m1cqZr90Nl5EhIndSMZIvu/C3NooSMz/7mst7771XajQavxCRI7Y8ri2TqcEi0gvAGAAPKKWGAICInBeRISJyo4hcrs+BRORQWlramvfff7/UhvGRA+MyMkTkjjiSz34yMjLw9NNPZxqNxkfr3rth6kymlFIPlBeVJyqlImvaT0SSy/81AlgNoF9TAktLS5vz//7f/8vKzOSyfe6AFxQickcceGM/s2bNysnNzZ0tIvm2PnadyZSIvFVeVB5fkTBdTSnVUinlX/F/ACMBNKkJTUTyLl++PPeRRx7JacpxyDnwgkJE7ogDb+xj9+7d+Pbbb38pKir6vDmO39CpEcKVUpcAzAHwhFLqklIqAEAYgO+UUgcB7AGwQUQ2NTW4goKCT7ds2XL6xx9/bOqhyMHxgkJE7ogDb5qf2WzG1KlTM9LS0qaIiNT9ioZr6Gi+VBFpIyIBIhJY/v88EflVROLKH7Ei8qwtgqsoRr/rrrsyiouLbXFIclC8oBCRO5qQEIVbekdBr8oGu+uV4qoPNrZo0aKC7OzsD0XkVHOdw2FmQK+JiPySlZW1+KGHHsrTOhZqPrygEJE7qlhGy1LeYGIRwZf7krgun438/PPPePPNN8+np6f/vTnP4/DJFABkZWW9vHbt2kMbNmyw1L03OSNeUIjIHXHwTfPJz8/H5MmTM4xG43gRacj8mA3mFMlUeXffzTNnzkxPS0vTOhxqBrygEJE74uCb5jNjxoycrKys+SJyprnP5RTJFACISHpmZua0W2+9NZuTeboeXlCIyB1x8E3zWLlypenbb7/9MS8vb5k9zuc0yRQAFBcXbzp58uSXr7zyCj9hXQwvKETkjjj4xvYuXLiA2bNnG41G4x3NNXrvak6VTAGA0Wh88IUXXkg6cOCA1qGQDfGCQkTuiuvy2Y7JZMLEiROz0tPTbxcRu81T6XTJlIiUGI3GMePHj2f9lAuZkBCFXte0umJbr2ta8YJCRC6L6/LZ3syZM/POnz//mtls/s6e53W6ZAoAROR0Wlran0ePHp1VUlKidThkA0+sOYzvz2Rdse37M1l4Ys1hjSIiImpeHHhjW6+//nrxxo0bv8nMzPynvc/tlMkUAJSWlm69cOHCs1OmTMm1U5coNaNPf7rYoO1ERM6OA29sZ+vWrdZnn332lNFovN1edVJVOW0yBQCZmZkv79y5c/3ixYsLtY6FmsZSw+9+TduJiJwdB97YxqlTp3DXXXelGo3GESKiSXeVUydTAGA0Gqe98sorx9evX88JPZ1Yxczn9d1OROTsOPCm6XJycjB69OjM1NTUsSKiWSG10ydTImJKT08fde+99yYfOXJE63Coke7o37ZB24mInNmaA0mVNVMVN41Rgb4cydcAZrMZN954Y7bRaPyriBzUMhanT6YAQEQyjUbjyNGjRxvPnj2rdTjUCP+a0ANTBlyDqg1RLTx16NOutXZBERE1g4pRfEnltVEWkcoWKSZS9WO1WnHHHXfknjx58vXLly//V+t4XCKZAgAROZGUlDR22LBhGcnJyVqHQ43Qp11r+Hj81uRdaLJiwarDXJ+PiFwKR/E1jYjgvvvuy9uxY8dyo9G4UOt4ABdKpgBARPZdunRp4vXXX5+ZkZGhdTjUQLzAEJE74Ci+ppk7d27++vXr16Snpz+idSwVXCqZAgCz2fzdxYsXp1x//fVZubm5WodDDcALDBG5A47ia7xnnnmm4JNPPtlqNBqnaTEFQk1cLpkCytbwu3Dhwv3Dhg3LLigo0DocqideYIjIHXAUX+O8+uqrRW+99daPRqPxVhFxqKniXTKZAoDLly//9+zZs4+NGjUqm7OkO4fqLjAKwLCuodoERERkYxzF1zjLli0rfe655w4ajcZxImLWOp6ruWwyBQDZ2dnvHz9+/JkbbriBLTTMa6gAAB2mSURBVFROYEJCFG7pHYWqM0sJgC/3JbEInYicHkfxNc57771XMn/+/CPp6ek3aDUpZ11cOpkCgMzMzFePHj0677rrrsvOybHbAtLUSN+eSMfVneAsQiciV8BBNg33/PPPFz3xxBM/p6enXysiDrvaicsnUwCQnZ299NSpUzMHDhyYZTQatQ6HasEidCJyVby+1Z+IYN68efkvv/zyt+np6cNFxKG/SW6RTAHA5cuXvzh37twdAwYMyLxw4YLW4VANWIRORK6K17f6sVqtmDlzZt7y5cvXGY3G8SJi0jqmurhNMgUARUVFWy5cuPCnwYMHp//yyy9ah0PV4CgXInJVvL7VzWw2Y/Lkybnr1q37yGg0ThERp1h310PrAOzNbDb/qJS6YdiwYVs2bNgQHh8fr3VIVEVFEebTXx1FdmHZzYi3h1vl/ETkYipG8CXnFCGwhSe8PXTILTIhMtCXxedVFBcXY/z48TmJiYmvG43Gp7SOpyHcLpkCABE5rJS6btSoUd8sW7Yscty4cfq6X0X2VGz6bQqRnCITFqw6DAC86BCRU6kYwVdReJ5daIKvpx6v3BbP61kVRqMRo0ePzr5w4cLCjIyM17WOp6Hc9pZfRE4bjcaEadOmHXrhhRcKHWgiVbfHES9E5Cp4Pavb4cOH0bdv3/Rjx47d4YyJFODGyRQAiEhmenr6gJdeemn9XXfdlVtaWqp1SASOeCEi18HrWe2++uory4gRIy5cuHBhSHFx8Wat42kst06mAEBESo1G4+2bNm16YfDgwdmZmZlah+T2OOKFiFwFr2fVExEsXry4aMaMGYlpaWkJInJC65iawqGTKaXUTUqpQ0qpRKXUXqXUtVWeu0cpdar8cU+V7UPL932hvucREcnIyHju6NGjd/fp0yfj+PHjtn4r1ADVjXgBgIISM2dCJyKnwhF8v1daWoo///nPua+88spao9E4UESytI6pqZQj1woppfwAFIiIKKV6AvhcRLoqpVoD2AugD8pWHNkHoLeIZCulPgMwFcC/APynodmuUio2LCxs0wcffBA5ZswYh042XdmaA0lXjOir4Oup5xpWROQUKkbxJeUUQa8ULCKIcvMRfOnp6RgzZkz2uXPnns/MzHxBHDkJaQCHThZEJL/KN7olULnSyCgA/xORLBHJBvA/AKPLn9OV72cFrljmrb7nPJqWltZr6tSp++fPn59vsTjFFBcuZ0JCFFp4/X6wKQs3icgZcB2+39u1a5f06tXLeOzYsTszMjKed5VECnDwZAoAlFITlVInAGwAML18cxSAi1V2u1S+DQCWAvgBgE5EGtVfJyLpRqNxwAcffPDOwIEDs1NSUhoZPTUFCzeJyFlxFN9vrFYr/vnPfxZOmjTpyKVLl/oUFhZ+rXVMtubwyZSIrBaRrgAmAPhn+ebqWpykfP/NItJLROY28byWtLS0eYmJibf26dMn9X//+5+17leRLbFwk4icFW8Gy2RkZGDo0KHZb7755odGo7G3iFys+1XOx+GSKaXUA+UF54lKqciK7SKyE0AnpVQIylqi2lZ5WRsAyc0RT2lp6bbk5OSEKVOmHGC3n31VV7ipAAzrGqpNQERE9cSbwbJuvYSEBOPevXvvTktLe8AZ1thrLEcvQO8M4Ex5AXovAF+hLHEKQlnRea/yXfejrAC92UYEKKX0oaGhi9q3b3/v2rVrW0dERDTXqaiKJ9Ycxic/XkDV31IWoRORI7p62Zj8YjNM1t+uXu5y7bJarXj22WcL33zzzTNGo3Gcq7ZGVeVwLVNXuQXAEaVUIoC3ANwmZbJQ1uX3c/njmeYeWikiFqPROC8xMXFynz59Ujdu3MhuPzv49kQ6rk733bXugIgcV9WCc0HZsjFQQKCvJxSAqEBft0ikUlJS3KJb72oO3TLlqJRSYWFhYZ8NHz487t///ndgQECA1iG5rA7zN/wumQLKuvvOLh5n73CIiKo1ePE3lSP3qooK9MX384drEJF9iQhWrFhhevTRR9MzMzOnlZaWbtE6Jnty9JYphyQiaWlpacPWr18/p3v37satW7eylaqZsO6AiJyBOxecG41GjBkzJnvu3LnrU1NTu7lbIgUwmWo0EZG8vLwPLl682OvOO+/8cfr06bn5+flah+VyOBs6ETkDd73x++9//2uOj49P3blz59TU1NSbRSRX65i0wGSqiUQkyWg0Xrt69eoF3bt3T9+5cyf7TW1oQkIUFt3cA0EtPK/YnlNkwoJVh5lQEZFDcLdlY7KysnDTTTdlz5o1a0tKSkpsYWHhOq1j0hJrpmxIKdXOYDB8ecstt3R56aWXAnx9XfuOxJ7cvR6BiBzT1SP4RIDcIhMiXXjZmPXr11v+8pe/ZOTk5MwqKCj4r9bxOAImUzamlNIFBQU9HBAQsOCdd94J5vp+tsFCdCJyNBUj+KrOdO7K0x8kJydj5syZ2Xv37t2XlpZ2p4gYtY7JUfCD3sZExJqVlfXK+fP/v727D2+yPvcA/r2fvNCmSZNCi20R5E0q1AKFooIORcZQ8YUBioBu0+PGdiETZROP87hzxnQcB5vs6HU5tyo6cGM67AU6BVSUKiKM0YqApV0tL22aviV9SZomeZ77/JHgKhR5adMnSe/PdT1X0qctvVVMv/n97t/vdzT/nnvueWfmzJnu48f7xMrQqOqr/QhCiNjVV46MCYVCWL16tX/ChAnHt27duqi2tnaGBKmvkjAVJcxcXVtbO3PHjh13FhQUHF+1alV7MJiwm79GnTSiCyFiTV9Ywbdr1y7k5uY2PPXUU//ncrlGBYPBhDtXrydImIqyQCCwra6u7tLf/OY3a8eMGdPw4Ycf6l1SXJJGdCFErEnkEfPGxkYsWrTIM2fOnD1HjhyZUldX9zAz+/WuK1ZJmOoFzNxRV1f3nxUVFVfOmzdv9x133OGpr6/Xu6y4Mzt/ECxm42n3E3FYXQgRm4r2V+PqVe9h2CNvwhcIwaTQVz4f7yv4NE3D888/H8jLy3O+/vrrP3S5XFcxc7nedcU6CVO9iJkrXS7XlDfffPP7Y8eOrVmzZo0/EAjoXVZc6QvD6kKI2JToR8Z8/PHHGDduXOPjjz/+stPpzPH5fBtZVqmdE1nNpxMisqSnpz+WkpJy369//ev+8+bNMxDR2b+xj5MtEoQQeknU15+Kigrcf//9TaWlpYdqa2t/wMyH9a4p3sjIlE6Y2VdfX//o0aNH85YuXfrXcePGNe7atUvvsmJeV43oBKDa046rV70nvVNCiKhJtJHxxsZGLF68uOWaa645sn379rlOp/MbEqQujIQpnTGzq7a2duGBAweumTt37oczZ85sqqio0LusmHWyEX1QpMGTgC/3n6r2tEszuhAiahKl4dzv9+OJJ55oz83NrVm/fv39LpdrtKqq7+tdVzyTab4YYzQar01PT3/+1ltvzXzyySdT09PT9S4pZiXqkLsQInacusN5mz+EoPbv35vxtEmnpmnYsGFD6NFHH3X7fL5nmpqanpIVej1DwlQMIiLFYrEstNls//ujH/3IsXz5covVatW7rJgju6ILIaKpqx3OTQZCitkYV0fGMDO2bt3KDz74YKPb7d7scrkeZuZGvetKJBKmYhgRJTkcjvuTkpJ+snTp0tRly5YlWywWvcuKGWcamTIQYc0d42L+BU4IEdviffSbmfHuu+/y8uXLG+vq6j6qra19iJkr9a4rEUnPVAxjZr/b7V5dW1s7bPXq1b8cOXJk3erVq/3t7fHZ7NjTzrQrusosvVNCiG6L54bznTt3YuLEiQ133333W59++unVTqdztgSp6JEwFQeYub2pqelJp9M5/Fe/+tWqkSNHup566im/1+vVuzRdnWxGN3SxpYRs5CmE6K54azhnZuzYsQMTJ05smD9//rb9+/dPdTqds5j5iN61JTqZ5otDRJSSlpb2QFJS0o+XLFmS+sADDyT35Z4q6Z0SQvSEzs3m2Y5kTLssA3/bV/2VnqlYbDhnZmzbto1XrFjRWFdXt9vpdP6UmT/Xu66+RMJUHCOi5EhP1UP33nuvbdmyZSkZGRl6l9Xr4r2vQQihv66azZNNBsydOAg7Pq//MmDFUsO5qqooKirSHn/88Sa32/1hJETJ3jo6kDCVAIion9Vq/W5KSsqjM2fOtP/sZz9zjBo1Su+yek1XL4JA+IiH/741N2Ze+IQQsSue3pR5vV4UFhYG1qxZ09zR0bHZ5XKtZOajetfVl0mYSiBEpCiKcsNFF1305GWXXXbxypUrB0yZMgV94Ziaov3V+J8tB8NnZXUSi0PyQojYEw/tAi6XC2vWrGl7+eWX24LB4HNNTU1rmdmjd11CwlTCIqIJ2dnZT6SlpRX8/Oc/7z9nzhzFYDh95Vsiiad3lkII/XXukVKIoHbx+zAWXj8OHz6MlStXut977z1PS0vLyvb29g3MHNC1KPEVspovCohoEhGpRDSv0z2ViEoi1+ZO93OJ6GMieomIeuy/BzP/s7q6+saDBw9OXLp06QvDhw+vX7t2bUcirwCM52XMQojedbI9oNrTDga6DFLJJgN+OjOn94tDuKn8gw8+wNSpUxumT5++b+PGjXe6XK4RPp/vRQlSsUdGpnoYERkAbAfgB/ACM78Wud/GzKctuSOiQgCPALgTQDkzvx2luhz9+/d/wGQyLb755pstDz30kH3MmDHR+FG6kU08hRDn6uteLzRm3ZrNW1pa8NJLLwXWrl3b4vP5djudzseYubRXixDnTcJUDyOiZQCCACYBeOMcwtQ6AD8BMB9AJTO/FeX6DIqi3JiZmfnYwIEDhz/88MMD5s6dq5jN5mj+2F5xpkZ0QHqnhBBfFWs9UiUlJVizZo17+/bt3o6Ojt97PJ7nmLmh1wsRF0Sm+XoQEQ0C8G0Az3Xx6SQi+gcR7Sai2Z3urwXwJoDJALZFu0ZmVlVVfaO6uvqqkpKSK5YtW/bs0KFD65YvX95WVVUV7R8fVbKJpxDiXMXChpx+vx/r1q1TL7/88oZZs2YVr1+//k6Xy3WJ2+3+pQSp+CIjUz2IiF4FsIaZd0dGnDqPTGUzcw0RDQfwHoDpzPwvHcv9EhGZ+/XrN69///6PjBgxInPFihXpN954I8Vrw/qZ3nECQFWMrMoRQvS+zg3nDosJbf4Qgtq/Xy16awS7vLwca9eubXnttdd8mqatr6+vX8vMJ6L6Q0VUGfUuIN4R0RIA3498aAfwl8hWBOkAbiKiEDMXMXMNADBzJRG9DyAfQEyEqUgz4ysAXiGi3IqKihUGg+GG73znO5b77rsvZfjw4XqXeF6yHcld9kIQwi+mMtUnRN9zahuA2xeEyUBwJJvQ3B6Meo+U1+vFpk2b1N/97ndN1dXVx+vr658MhUKbmTl49u8WsU5GpqKk88gUEaUB8DFzBxGlA/gYwG3MfEjXIr8GEaWYzea5AwYMeCAjI2PwkiVL0ubPn2+02+16l3ZWRfur8eDGki5Hp2JhmbMQovfpsXWKpmnYuXMnnn322cbi4uIOVVVfaWhoeC5WZiVEz5EwFSWnhKkpAH4PQEO4T+1pZi7Us77zQUQXp6Wl/YfZbL6noKAgZenSpenf/OY3EcvTgEMfefOMn5OpPiH6hs7Temf6TReNhvOKigr84Q9/aNuwYUO7pmnFTqfzaQAfMbPWoz9IxAwJU+KcUXj+siArK+vHAL41b968pMWLF6fm5ubqXdppzvQulAD8dv54meoTIsF93ereznpqZMrj8WDjxo3BZ555xtPU1HS0vr7+t8FgsIiZfd3+w0XMkzAlLggRmRVFmZWZmfmQzWbL+d73vmebP39+0rBhw/QuDYBM9QnR153pDVVn3W0493q9+Pvf/64VFhY2lpaW+jo6OgrdbnfhyR5Z0XdImBLdRkTpycnJc9LS0n5gtVqHLFq0KGXBggWWSy+9VNe6zjTVF0tnbQkhouPrVvUScMEN562trdiyZYv24osvNn722Wd+VVU31dfXrwNQyvILtc+SMCV6FBGlmc3m2enp6YuTk5OHL1iwwLJw4cKU0aNH93ot+b/YdtrBxwCQZjFh/+Pf6vV6hBDRd7JP6kyjUhcyMu3xeLB582a1sLCw6ciRI75gMLixsbHxZQCHJEAJQMKUiCIisptMplsGDhy42Gw259x+++3JixYtsubl5YG62Fizp435r7fgC57e7+lINqHk5xKmhEg0Z+uTOp9pvcbGRhQVFYVeeOEFd2VlZWsgEPhzU1PTemb+vKfrFvFPwpToFURkNRgMszIzMxcbDIbL586dm3TXXXfZ8vPzoxKsivZXY9nGkq5rgUzzCZGIvq5PatA5TOvV1dVh06ZNocLCQveJEyeaOzo6/uR2uzfIVgbibCRMiV5HRBZFUW7Iysr6IRGNnz17dr+77rorddKkSVCU7p1wdLYhfkAa0IVIRBf6BsrpdOLVV18Nrlu3zlNbW9vk8/leam5u/jMzV0WvWpFoJEwJXRFREoAZgwYNWqyqasGECROUOXPm9J8xY4ZhyJAh5/VnnetS6KdlawQhEsrZ/t/v/AbK6/WiuLgYW7Zsad66dWugvb29trW19aXW1taNcqSLuFASpkTMICIFwOV2u32W1WqdoyjKJddff73xtttuS5s2bRocDsdp39N5Uz6FCOpZ/j5Lv5QQiefrpveSDMC9ORp8lfvai4qK2lwul5eI3nE6nX8DUMzM3t6tViQiCVMiZhFRPwBXZWRkzDYajTfZ7fa0m2++OWnWrFm2yZMn461DDec0EnVSbx1iKoToXZ23QWBmhNw18FeVIHCshJM8VQ1JJsMup9P5V1VV32Vml67FioQkYUrEDSJyALguOzv7Dk3TvsGOizO0i8f3Mw4eB1PGJQgPbHXtXJpPhRDxadJjf0PVp5/AVHsg5Dt+KGA2m7/wehpf9LW4twAol+0LRLRJmBJxi4iGmPol3dA/a8h97T7fGHPmyCQeNNaQNDQfxtQMADIaJUQi8vl8KC4uxubNm5u3bt0aaGvvcAdg3OauqVoPYB8zh/SuUfQtEqZEQoj0W+WlZmTfriTb71MD7WlpI8YpC+fcYlxx79wu+62EEPFBVVXs27cPb731VntRUVFbbW2tl4jejfQ97ZS+J6E3CVMiIUX6rSYPHDhwtsFguMlutzumT59unjp1qn3SpEkYOnRor2wcKoQ4f83Nzdi3bx92797tf+edd1oOHTqkmUymvU6nc2Ok76lW7xqF6EzClOgTiCgNwBUDBgy4Njk5eZqqqkOHDBmiXHvttZZrrrnGOmnSJGRmZupdphB9jt/vR0lJCXbv3h3csWOHp6SkBMFg0E1Eu51O53Zm3gPpexIxTsKU6LOIaBCASRdddNH1RqNxKjNn5eTkYNq0aamTJ09OKigokOlBIXpQKBTCwYMHsWfPHnXHjh2evXv3am1tbV6j0fhPl8u1LRgMfgLgIDOffqimEDFMwpS4IER0HYCnAZgANDDztZH7NwBYC8AA4I/MvCpyPxfAHwEcAXAPM59+aJ7OKDzvN4KIJmVlZc1g5slGo3HA2LFjMW3aNPtVV11lzs/Ph8Vi0btUIWKepmmoqKjAnj17eOfOnZ6PPvoo2NTUFDAajZ+53e53vF7vLgAlzHzm4wqEiBMSpsR5i2xRsAvADcx8jIgGMnMdERkQDkszAJwAsBfAAmY+RESFAB4BcCfCQ/Zv61X/+SAiI4DRRqPxiszMzG+FQqECi8ViKygoUK677jrHlVdeacjLy4PJZNK7VCF0w8w4ceIE9u7di+Li4pYPPvjA73Q6VaPRWN7W1vaex+MpRniVXbPetQoRDUa9CxBxaSGATcx8DACYuS5y/woAFcxcCQBE9BcAtwE4hPBIFQPQED4qKy5EllgfiFyFQPgInMrKyrFvvPHG5AEDBswIBoN5DocjefLkycYpU6Y4LrvsMsrJyUF6ero0uYuE4/f7UV5ejrKyMpSUlHjff/99b2VlJSuKciwQCHxQX1//PoC9nV4XhEh4MjIlzhsRnZzeywVgA7CWmV8monkIj1bdF/m6uwFcycz3E1E+gOcAlAP4LjOf27blcYKIbAAmWCyWif379y9g5lxN0zJSU1ONY8aMofz8fGtubm5STk4ORo4ciX79+uldshBnxMyorq5GWVkZDh8+rJaUlLQcOHAgWFNTA1VVW41GY4XP59vf2Nj4DwD/AHBMGsRFXyZhSpw3InoGQAGA6QCSAXwMYBaAcQBmnhKmrmDmpXrVqjciSgFwKYCc9PT0CcnJyfnBYHCkyWRKyc7Oxrhx40zjx4+35+TkKDk5OcjOzpbRLNFrWltbceTIEZSVleHTTz9tKy0t9ZWXl5PX6w0YjcYaTdMO1NfX7w0Gg4cBlAFwSWgS4nQyzSfOCREtAfD9yId/BfB2ZKM8LxHtRDhInQAwuNO3XQygplcLjTGRf0clkWvjyftERMePH8/65JNPckwm0+iMjIxJiqJcHgqFsq1Wq2nUqFEYP368JS8vLyUnJwejRo1CSkqKbv8cIn6pqoqjR4+irKwMhw4dCkZGmdSGhgbWNM1jNBrLWltb/+nxeEoRDkz/YuaA3nULEU9kZEqcNyIaDeAZADMBmAHsQbix/HOEG9CnA6hGuAF9ITMf1KnUuBTZcHQEgByHwzHOZrNNDIVCOYqipA4cOJDy8vKMubm51sGDB5uzs7ORnZ2NrKwsWK1WvUsXOgiFQnC5XKipqTl5aeXl5W2lpaX+yspK8vv9fpPJdCwYDJa6XK69zFwGoIyZm/SuXYhEIWFKXBAi+imAexBuKP8jMz8duX8TwlsmGAC8wMxP6Fdl4iGidAA5AIY7HI5hKSkpIxRFGRIKhbIAWE0mk8nhcGDw4MEYMmSIediwYZbBgwebs7KyJHTFma5CUlVVlffo0aMdx44d02praykQCIRCoVCHyWSqI6ITfr//i8bGxnJN044jPMpUlWj9iULEIglTQiSQyF5ZdgBZALIBZNnt9mFWq3WkhK7YcKaQVFVV1XH8+PEvQ5Kqqh1Go/HUkFSD8NR5DYB6OdBXiNggYUqIPigSulIRDlxfhi6bzTaCiC7pHLosFgtZrVZOTU2Fw+Egu91ucDgcBofDYbLb7f3sdrvBZrPBarXCZrPh1OeJsnJRVVW0tbWhtbX1y8fOz1taWri5uTng8XgCHo8n5PF41ObmZq25uZlbW1vJ6/V+XUiqBuBEOCTVyWiSEPFFwpQQ4owioasfwltgWCOPX3melJTkSElJ6W82mwcYDAaHoih2Zk7VNM2mqqoVgFlRFIPBYFAURTGYzWbYbDZOTU1FamoqORwOxeFwGB0Oh8lms5mMRqNiMBhIURQyGAxkMBigKMppV+f7RARN06CqKjRNO+06eV9VVdY0jU8++ny+UHNzc9Dtdn8ZflpaWtDa2ko+n49UVdU0TVMjV8hgMPgMBkMrEbUyc4umae5gMOj2+XyNXq+3iZlbAbQCaIs8dv7YG4s7/wshuk/ClBCiVxGRGV2EssjzFABKp8tARCeTmMFoNBoVRTEQkVFRFEPkuQGAwswqM6uapoUil6qG09DJxxDCPX6drwC+Gno6ByGfbAMghDgXEqaEEEIIIbpB0bsAIYQQQoh4JmFKCCGEEKIbJEwJIYQQQnSDHCcjRDdENi9dFPnQCGA0gAxmbiKiKoQbmVUAIWYuiHxPNoD1kc8tYua2Xi9cCCFEj5EGdCF6CBHdAuBBZr4+8nEVgAJmbjjl61YB+BOA4QAGMfNzvV2rEEKIniPTfEL0nAUA/nwOX2fAv5fmU1QrEkIIEXUyMiVEDyAiC4ATAEaePECWiL4A4AbAAH7PzM9H7l+C8MhUM8IHQbfqU7UQQoieID1TQvSMWwB8dDJIRVzNzDVENBDAdiL6nJl3MvNRAFP1KVMIIURPk2k+Ic4TES0hopLIlR25fSdOmeJj5prIYx2A1wFc0buVCiGE6A0yzSdENxGRHcAXAAYzszdyLwXhI05aI8+3A/gFM7+tY6lCCCGiQKb5hOi+bwPYdjJIRVwE4PXwOcEwAnhFgpQQQiQmGZkSQgghhOgG6ZkSQgghhOgGCVNCCCGEEN0gYUoIIYQQohskTAkhhBBCdIOEKSGEEEKIbpAwJYQQQgjRDRKmhBBCCCG64f8B9SY2jE3LgW8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 720x432 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(10, 6))\n",
"plt.subplot(projection='hammer')\n",
"plt.scatter(plane.icrs.ra.wrap_at(180*u.deg).rad, plane.icrs.dec.rad)"
]
},
{
"cell_type": "code",
"execution_count": 76,
"metadata": {},
"outputs": [],
"source": [
"width = 20*u.deg\n",
"\n",
"upper = SkyCoord(b=width/2, l=np.linspace(0, 360,100)*u.deg, frame='galactic')\n",
"lower = SkyCoord(b=-width/2, l=np.linspace(0, 360,100)*u.deg, frame='galactic')\n",
"\n",
"combined_b = np.concatenate([upper.b.deg[::-1], lower.b.deg])*u.deg\n",
"combined_l = np.concatenate([upper.l.deg[::-1], lower.l.deg])*u.deg\n",
"combined = SkyCoord(b=combined_b, l=combined_l, frame='galactic')"
]
},
{
"cell_type": "code",
"execution_count": 80,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.patches.Polygon at 0x119975a90>]"
]
},
"execution_count": 80,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlMAAAFACAYAAACY6/lAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4wLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvqOYd8AAAIABJREFUeJzs3Xlc1HX+B/DXZ2a4hxtmuFQQEfACFMEjFSM8MFG8ytRSs922Y63217m11V7du5Vra7tlaWZY5lWeeWGkSKaIHKKioJzDfcNcn98f4KzmxTHwmRnez8djHsLwnfm+xoGZ93xOxjkHIYQQQgjpGonoAIQQQggh5oyKKUIIIYSQbqBiihBCCCGkG6iYIoQQQgjpBiqmCCGEEEK6gYopQgghhJBuoGKKkG5ijAUzxtKvudQxxp5ijL3GGCu65vr4a27zDmPsBGNsksjshBBCuo/ROlOEGA9jTAqgCEA0gGUAGjjn7/7qmBAADwN4FcDnnPMFvR6UEEKI0VDLFCHGFQsgj3NecJtjpAD0ADgA1iupCCGE9BgqpggxrvsBfHXN908wxjIYY2sZY64AwDnPAmAPIAXAvwVkJIQQYkTUzUeIkTDGrAEUAxjKOS9jjCkBVKCtBeovALw558tFZiSEEGJ81DJFiPFMB3CSc14GAJzzMs65jnOuB/BfAFFC0xFCCOkRVEwRYjwLcU0XH2PM+5qfJQLI7PVEhBBCehx18xFiBIwxewBXAAzknNe2X/cFgHC0dfPlA/gt57xEWEhCCCE9goopQgghhJBuoG4+QgghhJBuoGKKEEIIIaQbqJgihBBCCOkGmegAhJC+izFmBcCx/SIH4IC2D3k3u0hv8zOGtlXlr73obnLdry+tABoA1LdfmjkNJCWEdBINQCeEdBhjzAb/K3wccX0h5Ghtbe0sl8vdbWxs3KRSqStjzAWAo16vd9TpdHLOuY1UKpVIJBKpRCKRWltbw9HRkTs5OXEnJycml8slVlZWTCKRQCqVGv5t/5rJZDJIJBLD91d/BgB6vZ7rdDqu1+u5Vqvler0eV7/X6XRcp9NdPQZ6vZ7r9Xo0NTXpa2treV1dHa+vr2dNTU0SvV6v0+v1Op1Op+eca6VSaaNUKq1njNXp9fo6vV5fo9FoqpqamqoaGxur8L9CrOEmXzdScUaI5aNiipA+jDEmA6AE4A3ARyaT+bi7uwdZW1sPBOCn0Wg8JRKJlUQikUqlUomNjQ3kcjl3cnKCs7Mzc3FxkTo7O0tdXV2tXFxcbJycnKSOjo64epHL5dd9bWtrC8bMZztCrVaLhoYG1NfXGy6/+l5fU1Ojrqmp0VRXV2tra2t1tbW1+traWjQ0NKCpqYm1F3Y6nU6nk0gktTKZrFir1RbU1dWdb2houIy2VfOLAZRwzhsFP2RCSBdQMUWIBbpJkeTr7u4+yNraeiDn3Fer1SqkUqmNjY2N1Nvbm/fr10/q7+9vM2DAAAdfX1/m7e0NHx8fKBQKyGQ0GsAYOOeora1FcXExSkpKUFxcjCtXrrRcunSpuaCgQFNUVMTq6+u5RqPRSCSSmlsUXSUAiqnoIsS0UDFFuqy9C+cTAMPQtjDlcgC5ADYB8EfbQpULOOfVjDEJgM8BDALwSPtmv6QLGGMeAAIAeN+qSLK1tZV6eXnxfv36Sfz9/W39/f0dfHx8mI+PD7y9valIMmFXi66rBde1Rdfly5c1hYWFrKGhQa9Wq7USiaRGKpUW63S6gtra2guNjY0FaCu6CgEUcM51gh8OIX0CFVOkyxhj6wD8yDn/pH2TX3sALwGo4py/yRh7AYAr5/x5xtg0AIEANgN4gzb8vb32sUmBAIKdnZ1HODk5RWq12mCJROKkUChYcHDwDUWSj48PPD09qUjqI35ddJWUlODy5cst+fn5TefPn9fm5+ejpaWlVSaTXdZqtafLysp+5pznAsjlnFeJzk+IJaFiinQJY8wJwGm0bZ/Cr7k+F0AM57ykfW+6w5zzYMbYDAD90VZMvc05XyYkuAlhbYOHvAEEW1lZhSoUikjG2HCtVusjl8utgoKCeEREhP3w4cPlwcHBCAoKglwuFx2bmBGdToeCggLk5uYiOztbnZ6eXn/mzBldRUWFXq/X10ql0ty6uroTdXV1Z9DWqpzHOVeLzk2IuaFiinQJYywcwH8AZAMIA/ALgJUAijjnLtccV805d20fw7MBwAAAv+GcnxEQW4j2ffsGAwh2d3ePsLe3H6nRaAbJZDIHX19fhIWFWY0YMcIpNDRUGhwcDB8fH7MapE3MU319Pc6dO4fc3FxkZGQ0nD59uun8+fOssbFRLZPJivV6/RmVSpWm1WrPoq3QKqOZiYTcHBVTpEsYY5EAUgGM55wfZ4x9AKAOwJM3K6ZE5exNjDEvAGF2dnZD3N3dR3POh+r1ek9nZ2dZaGgoCw8Pdxg2bJhdcHAwAgMDYWtrKzoyITfgnKOoqAi5ubnIycnRpqen1505c0ZbXFwMrVbbYGVldb6pqelkZWVlOto+TOXQ2CzS11ExRbqkvXBI5Zz7t38/AcALaBtgfkM3n7ikPaN98H2kq6vrBAcHh1itVhvo5+cnGTNmjE1YWJhTSEgICw4OhoeHB7UyEYvR3NyMCxcuIDc3F5mZmc0///xzQ2ZmJtdoNFUAjpWUlPwA4Ge0dRfSmwvpM6iYIl3GGPsRwArOeS5j7DW0rV4NAJXXDEB345w/JyykEbR300XI5fJxrq6ucRqNZoi7u7v1XXfdZTVp0iSX0aNHIzAwkIom0mdVV1fjl19+wdGjR5sPHz5cf+7cOc45L9FqtUdUKtUhAD9zzotE5ySkp1AxRbqsfdzUJwCsAVwEsAxtW3t8jbbB5pcBzDenmUPt25sMt7Gxifb09Jyi1WojHB0d7aOjo6UxMTEuUVFRktDQUJoxR8gdlJSU4MSJE0hJSWlITk5uKigo4DKZ7GJjY+Oh6urqI2grsMzmtYGQ26FiivRZ7WtfBUul0igvL684vV4fZWNj4xwRESGZPHmyS3R0tCwsLAw2NjaioxJi9jjnuHTpEk6cOIEjR47U/vjjj61lZWVaKyurszU1NfsbGhp+AnCSc94gOishnUXFFOkT2pchGABgtJeXV6xEIhkvkUgUQ4YMweTJk53Hjh1rM3LkSDg6OoqOSkifodPpkJubi7S0NJ6cnFx99OhRXW1tbYuVldXpioqKfS0tLakAMjjnraKzEnI7VEwRi8UY85ZKpfd4e3vfr9VqRw4aNEgyefJk+bhx4+wjIyPh4eEhOiIh5FfUajUyMjKQlpamPXz4cE1aWhrXaDTFTU1NW2tqar4DkM4514vOSci1qJgiFoMx5ghgore393zO+WRfX1/bxMREx2nTptmFh4dDKpWKjkgI6YLi4mLs379fv3Xr1qrjx4/rpVJpdkVFxaaWlpa9nPNLovMRQsUUMVvtg8VHu7u7J1hbWyfI5XKP+Ph4m3vvvddp/PjxsLOzEx2REGJknHPk5ORg79696m+//bb2woULGplMdrSkpGSTXq8/xDmvFJ2R9D1UTBGz0T7uKcTBwWG6i4vLfMZY4IQJEySzZ892u/vuuxl12xHS92g0GqSlpWHnzp2NO3bsaKqsrKzjnO8pKyvbCuAnznmL6IzE8lExRUzaNeOe7tNqtaOGDRsmnTNnjuvUqVNlAwcOFB2PEGJi6urqkJycjB07dtTs379fo1ari9rHW30PGm9FeggVU8SkXB335OPjM0+v19/t4+NjO3v2bMfp06fbRURE0LgnQkinFBUVGcZbpaWl6aVSaVZFRcXXNN6KGBMVU0Q4xthgV1fXRba2tvPlcrnH9OnTDeOe7O3tRccjhFgIzjmys7Oxd+9e9ZYtW2ry8vI0Eokkpbi4+FO0bX2lEZ2RmCcqpogQjLEh7u7uD1pZWd0XFBTk8PDDD7vFx8dLPT09RUcjhPQRGo0GP/30EzZs2FCze/duDWPsWFFR0ccA9nPO1aLzEfNBxRTpFe2Dx4d7eno+JJFI5g8dOtR2+fLl7jNnzpQ4OTmJjkcI6eP0ej1SU1OxYcOG2h07dqgB/NxeWO2jQezkTqiYIj2mvYAaqVQqlzLGZoeFhdkuX77cLT4+XiKXy0XHIxZEp9OhsrIS1dXVUKvVhktraytaW1vR0tKC1tZWSKVS2NjYGC7W1tbXXVxdXeHm5kZj8/o4zjlOnDiBDRs21G3durVVr9efLikpWaPX63dzzptE5yOmh4opYlTtBVSUUqlcDuDe0aNHWy9btsx92rRpjMY/ke5qbm5GeXk5KisrUVlZiYqKClRUVKCqqgp6vXEmaTHG4OzsDDc3N7i5ucHV1dVw8fT0pE2u+xjOOdLT0/Hll182bN68uVmr1WaVlpZ+rNPpvqd9BMlVVEyRbmvfMHicl5fXcgDTx40bZ7V06VL3uLg42Nraio5HzFhraysuX76M/Px8XLp0CSUlJRD5miWTyeDl5QU/Pz/4+vrC19cXbm5uwvKQ3sU5R1ZWFjZu3Ni4adOmptbW1nMqlepjjUazg3NeKzofEYeKKdJljLEhCoXiccbYnJiYGOulS5e63X333bC2thYdjZip1tZW5OXl4fLly7h8+TJKS0uN1uLUU+zs7ODj4wNvb2/4+fkhICAANjY2omORXpCTk4OkpKSmL7/8srG5uTmntLT0vfauQJoV2MdQMUU6hTHmLpfLH5TL5Y8NHjzY5cknn3SbOXOmhN48SFe1trbi3LlzyM7Oxvnz56HVakVH6haJRAI/Pz8EBgZi4MCB8PX1hUQiER2L9LDTp0/j448/rt22bVsr53xXaWnpKgCnOL3J9glUTPUQxpgUwAkARZzzexljnwOYBOBqU/BSznl6exfZ5wAGAXiEc54lIu/tMMasJRJJvJeX1x8cHByCH330UafFixfbKBQK0dGImWpoaEBubi7Onj2LS5cumX0BdTs2Njbw9/dHQEAAAgICoFAo0Da0kFgirVaLvXv3YtWqVRWnT59uaGlp+bSmpuZTznmJ6Gyk51Ax1UMYY88AiATgdE0x9T3nfPOvjpsGIBDAZgBvcM6X93rYm2gfSD7K29v7SQDT5s6da/ub3/zGafjw4aKjETOl0WiQlZWFkydP4sqVK0LHPonk6OiI0NBQhISEwN/fn1qtLFh1dTWSkpI0q1evrqmuri4oLy//h0aj2cY5bxadjRgXFVM9gDHmB2AdgL8BeOYOxdQMAP3RVky9zTlf1tt5f5XH2dHRcZm9vf3vR40a5fjkk0963HPPPTSDiXRZcXExTp06hTNnzqClhZbruZadnR0GDx6MkJAQBAYG0nhDC3b+/Hn897//bdi4cWOzTqfbW1pa+qYp9kSQrqFiqgcwxjYDeAOAI4D/u6aYGgugFcABAC9wzlsZYzIAGwAMAPAbzvkZQZkjvLy8XrS2tp782GOPyVesWGHr7u4uIgqxAPX19Thz5gzS09OhUqlExzELMpkMAQEBCA4ORkhICGgtNsuk0+mwZ88e/tZbb1VcuHChrKqq6s3W1tbNnPNW0dlI11ExZWSMsXsBxHPOH2OMxeB/xZQ3gFIA1gD+AyCPc/5ngVHBGLOzs7Nb6Ozs/NyQIUPcX3jhBY/Y2FjqdiBddvnyZRw9ehS5ubl9thvPGCQSCQICAjBixAiEhoZSi5WFKigowKpVqxo2btzYpNPpvlapVO9xzvNF5yKdR8WUkTHG3gCwBIAWgC0AJwBbOOeLrzkmBu1FlqCMQUql8jmpVJqwdOlSh8cee8zB19dXRBRiAXQ6HbKyspCamori4mLRcSyOlZUVgoODMXz4cAwaNIhWZ7dAarUaW7Zs0b399ttVZWVlF0tLS//avsSCTnQ20jFUTPWgX7dMcc5L2gd2/xNAC+f8hV7MIpHJZDM9PT1f7devX7/nn3/eLSEhQUJjoUhX1dfX45dffsGJEyfQ0EALQfcGOzs7DBkyBMOGDYO/vz/NCrRA2dnZ+Mc//lH7/fffN2k0mk+qqqre55xXic5Fbo+KqR70q2LqIABPAAxAOoBHe2MrAsaYnVwuX+7g4PB8fHy844svvugSFBTU06clFqy0tBQpKSnIycmBTkcfnEVxcnJCREQERo4cCWdnZ9FxiJE1Njbi888/V7/zzjt1LS0tu8rKyl7nnF8UnYvcHBVTFoox5uHu7v5/1tbWyx555BH573//e3saUE66o6ioCEeOHMG5c+doPJQJYYwhMDAQI0eORHBwMHUDWhi9Xo8dO3boX3vttSqVSpVRUlLyEuf8uOhc5HpUTFkYxliQl5fXa3Z2dlNefPFF5yVLlljR/nikO/Lz8/Hjjz8iLy9PdBRyBw4ODggLC0NERAQ8PT1FxyFGlpqaildffbXyzJkzJeXl5S9rtdrvOOemvd9SH0HFlIVgjI3z8fF5w9vbe8hrr73mHh8fz2hWHukqzjlyc3ORkpKCwsJC0XFIF/Tr1w+jRo3CsGHDaJ04C3Px4kX8/e9/r/3+++/rGhsb32poaFhLC4GKRcWUGWOMMYlEMkWhULw3atQor9dff9191KhRomMRM3fu3DkcOHAAZWVloqMQI7C3t8eoUaMwevRoODk5iY5DjKiqqgqrVq1q+vjjj+vVavWaysrK9zjn9aJz9UVUTJkhxhiTyWQzPDw83p00aZLir3/9q+ugQYNExyJmLj8/HwcOHMCVK1dERyE9QCKRIDQ0FNHR0ejfv7/oOMSImpubsWbNmpb33nuvrqWl5ZPKysq3Oee1d74lMRYqpsxIexGV4O7u/m5sbKzHX//6V5eAgADRsYiZKyoqwsGDB2lMVB/i7e2NqKgoDB8+nLoALUhLSws++eST1jfffLOutbV1XUVFxd845zWic/UFVEyZgfYiKt7d3f2fU6ZM8fzLX/7iMmDAANGxiJkrLS3FoUOHkJubKzoKEcTe3h6RkZGIjo6Gg4OD6DjESNRqNT755BP13//+99r2lqo3Oed1onNZMiqmTFj7mKg4T0/PDyZPnqx84403XP39/UXHImauvr4eBw4cwOnTp2mJAwKgbV/A8PBwjBs3Dm5ubqLjECNpbW3Ff/7zn9a33nqrtrm5eU1VVdU7vbG+YV9ExZSJYoxFKpXKT8aPH9//7bffdg0MDBQdiZg5jUaDo0eP4qeffoJarRYdh5ggiUSCIUOGYPz48fD29hYdhxhJS0sLPvroo9Z33323trm5+Y2ampp/cc61onNZEiqmTAxjzN/Ly+ujgQMHRv373/92HzFihOhIxMxxzpGRkYEDBw6gro5a+knHBAYGYvz48Rg4cKDoKMRI6uvr8de//rVh3bp1ldXV1U+p1ertnIoAo6BiykQwxlwUCsUbzs7O81atWuU+depU2nSLdFt+fj727dtHGxCTLvPx8cGECRMQEhJCewFaiKKiIvzhD3+oOXz4cEFZWdkKzvkJ0ZnMHRVTgjHGrF1dXZ+yt7f/w5///GfXhx56yIq2gyDdVVtbiz179iAnJ0d0FGIhvLy8EBMTg5CQENFRiJFkZGTgd7/7XeXFixfTSktLH+Oc54vOZK6omBKEMcbs7Ozuc3JyeufRRx91e+655+zt7e1FxyJmTqfTITU1FcnJyTQuivQIb29vxMTEIDg4WHQUYiT79u3jTz75ZGVNTc23KpXqRc55tehM5oaKKQEYYyMUCsWX06dP93vrrbdclEql6EjEAuTn52Pnzp0oLy8XHYX0AT4+PoiJicHgwYNFRyFGoNPpsH79es0rr7xS3dDQ8Jfa2tqPaN+/jqNiqhcxxhwVCsU/vLy8EtevX+8eFhYmOhKxAPX19di3bx/OnDkjOgrpg/z8/BATEwPahcEyNDQ04KWXXqrbtGlTkUqlWsQ5PyU6kzmgYqoXtHfp3e/s7PyP119/3W3FihXWtAkx6S7OOdLS0nDw4EG0traKjkP6OH9/f8TFxcHX11d0FGIEmZmZePDBByuLioq+U6lUT9H2NLdHxVQPY4wFKZXKL+++++6gDz/80MXDw0N0JGIBVCoVduzYgcLCQtFRzIJEIoFcLoe1tTVkMhlkMhmsrKyu+1er1UKj0UCr1V73tUajQUNDA3Q6neiHYRaGDh2K2NhYWvzTAnDO8dlnn6lffvnlqtra2mebmpq+pKUUbo6KqR7CGLPz8PD4i5ub29LPP//cfezYsaIjEQug0+mQkpKCI0eO0Jv7NaRSKTw8PKBUKuHi4gJHR0c4OjrCyckJjo6OkMvl3ZrWzzlHU1MT6uvrUV9fj4aGBtTX16Ourg7l5eUoLy9HU1OTER+ReZNKpRg1ahQmTZpE29RYgKqqKjz11FO1+/btu1BWVraIc057UP0KFVM9wNra+h5XV9fPnn32WY+nnnrKljYSJcZQVFSEHTt2oKysTHQUoZycnODr6wuFQmG4uLu7Q3TXeWNjI8rLy1FRUYHy8nKoVCoUFxf36S5YGxsb3HXXXRgzZgysrKxExyHddPz4cSxdurSyqqpqffusv777y/0rVEwZEWNMrlQq1wQHB8dv3LjRlcYOEGPQaDQ4dOgQUlNTodf3vck1bm5uGDBgAPr3748BAwaYVfcR5xxlZWUoLCzElStXUFhYiMrKStGxep2TkxNiY2MxYsQIWvjTzOl0OrzzzjvN//znP4tVKtU8znm66EymgIopI5HJZBM9PDy+/POf/6x45JFHrOkFgxhDYWEhtm7d2qfegO3s7DB48GAMHjwYAwYMgFwuFx3JqJqamlBYWIi8vDxcuHChTz23fn5+mD59Og1StwA5OTlYsGBBVWlp6ScVFRUvc841ojOJRMVUNzHG7BQKxaqBAwcmbtq0ya1///6iIxELoNPpkJycjJSUlD7RGuXm5obg4GAEBwejf//+wrvselN1dbWhsLp06ZLFdwsyxhAeHo577rmHxlOZOa1Wi7/97W9NH330UaFKpZrDOc8SnUkUKqa6gTE2RqFQbHr55ZcVTzzxhC21RhFjKC8vx5YtW1BSUiI6So9yd3dHeHg4QkNDQbNc2+j1ehQWFiInJwfZ2dmorbXc2ei2traIiYlBVFRUnyqeLVFGRgbuu+++yoqKitUVFRV/5pz3udkxVEx1AWPMxtPT871+/fot/Oabb9xoV3ViDJxzHD9+HPv374dWqxUdp0fY29tj6NChCAsLg5+fn+g4Jo1zjsLCQmRnZ1t0YaVQKBAfHw9/f3/RUUg3aDQa/OlPf2pYu3ZtgUqlmtvXZvyZfDHFGLMFcASADQAZgM2c81cZYwEAkgC4ATgJYAnnXM0YkwPYCEAOYDHnvNjIeYIUCsXOZ5991u+ZZ56xo09UxBjq6uqwbds2XLx4UXQUo5NKpQgODsaIESMQFBQE2si78zjnKCoqQnZ2NjIzM1FXVyc6ktENHz4cU6dOtbgxcn3NyZMnsWDBgoqKioqXa2pqPhadp7eYQzHFADhwzhsYY1YAUgCsBPAMgC2c8yTG2BoApznn/2aMPQrgCoDLABZxzl8wVhZHR8dFHh4e72/ZssUjIiLCWHdL+rjc3Fxs27YNzc3NoqMYlYODA0aNGoXRo0fD0dFRdByLodfrkZeXh1OnTiE3N9ei1huzs7PDPffcg5EjR9KsPzPW2NiIZcuW1Rw5ciS5fV2qRtGZeprJF1PXYozZo62Y+h2AnQC8OOdaxthYAK9xzqcyxh4HcBFtxdSDnPPnjXBeW4VC8WlkZGT8V1995eLk5NTduyQEOp0O+/btw/Hjx0VHMSpfX19ERUVh6NChoDXWelZTUxMyMjJw6tQpi1p/bMCAAZg5cyaNpTNzn376qfqll14qUqlUMy19cLpZFFOMMSmAXwAMArAawDsAUjnng9p/3g/Abs75MMaYC4CvANiireuvW/ttXO3We+WVV/o9/vjjNMicGEVVVRU2b96M4mKj9kILI5FIMHToUIwdOxY+Pj6i4/RJxcXFOHnyJE6fPg2NxvxnqctkMkyYMAF33XUXdQ2bsaysLMyaNauyoqLij5bc7WcWxdRV7YXSVgB/AvDZr4qpXZzz4cY8n6Oj42IPD49/UrceMabMzEx89913FjEFXiqVYsSIEZgwYYJZLaZpyZqbm3Hy5EmkpaVZxKB1hUKBhIQEmrBgxq52+yUnJyerVCqL7PYzq2IKABhjrwJoAvA8btLNZ6Rz2CiVys9GjRo1nbr1iLFotVrs2bMHJ06cEB2l26ysrBAREYHx48fD2dlZdBxyE3q9HmfPnsXx48dRUFAgOk63MMYwduxY3H333dR1bMbau/2KVSrVDM55tug8xmTyxRRjzBOAhnNewxizA7APwFsAHgLw7TUD0DM45x8Z4XwKT0/P/S+++OKgp556yo669Ygx1NTU4Ouvvzb7bj0rKytER0djzJgxNOvKjJSUlCA1NRVnzpwx60VgPTw8kJiYSCuom7GsrCzMmDGjQqVSLWtqavpedB5jMYdiagSAdQCkACQAvuac/5kxNhD/WxrhFNqWQehWvwljbLhSqdy9fv167ylTptCaB8Qo8vLy8O2336KpqUl0lC6TSCSIiIhATEwMzcwzYzU1NUhJSUF6errZrmUmkUhw1113YdKkSTSWykxVVlZi+vTpVRcvXnyvsrLyDW7qhUgHmHwx1VvkcvkshULxya5duzxCQkJExyEWIiUlBQcOHIA5/52FhoYiNjaWZlZZkPr6evz000/45ZdfzHawure3NxITE6FQKERHIV2gVquxfPny2h9++GGvSqVawjlXi87UHX2+mGKMMQ8Pjz8OHDjwmd27d7vSIFpiDGq1Gtu2bUN2tvkOCwgICEBsbCwN/LVgTU1NOHbsGH7++We0tLSIjtNpMpkMsbGxGDNmDK1LZYY453j33Xeb3nnnndzy8vI4zrnZ7vrdp4spxpi1QqHYOGXKlHvWrl3rbGVlJToSsQBVVVVISkqCSqUSHaVLXFxcMG3aNFALbd/R3NyMlJQUpKWlmWVLVUBAAObMmUNd0GZq586duuXLlxerVKqpnPMc0Xm6os8WU4wxV09Pz4MvvPBC8NNPP00DzYlRXLx4Ed98841ZrmYuk8kwfvx43HXXXaAPFn1TXV0dkpOTcerUKbMbqG5vb4/Zs2dj8ODBoqOQLsjKykJ8fHxwNr4WAAAgAElEQVRFWVnZ4paWlr2i83RWnyymGGNenp6eP37yySf+CQkJNM+WGMWJEyewa9cus3sTAoCQkBBMnToVrq6uoqMQE1BZWYmDBw8iOzvb7Mb7jRkzBnFxcTQ43QyVl5dj8uTJVZcvX36srq5uk+g8ndHniinGmL9SqUxOSkryi4mJoRl7pNs459i7dy9SU1NFR+k0d3d3xMfHIzAwUHQUYoKKi4uxf/9+s9uA28fHB/PmzaOFZM1QbW0tYmNjq/Py8l6srq42mxXT+1QxxRgL9fLyOvDdd995R0ZGio5DLIBarcbmzZtx7tw50VE6RSKRYOzYsZg8eTItgkjuKCcnB/v27UN1dbXoKB1mY2ODhIQEDB06VHQU0knNzc2YPn16dVZW1jvl5eVviM7TEX2mmGKMjfL19d21b98+xZAhQ0THIRagrq4OGzduRGlpqegoneLl5YWEhATaQ490ilarxdGjR5GSkgK12nxmsUdFRWHq1KnU7Wdm1Go15s2bV5OamvpJeXn5c6a+FlWfKKZsbGwmeXt7bz506JBHQECA6DjEApSVlWHDhg2or68XHaXDZDIZJk2ahPHjx0MioR5u0jV1dXXYt28fMjMzRUfpMD8/PyxYsAC0NZh50el0WLZsWd2ePXu+LS8vX8E5N9kBqRZfTNnb28f7+fmtT05Odvf29hYdh1iAS5cuISkpyaw2Ku7Xrx9mzZpFC28SoykoKMCuXbtQVlYmOkqHyOVyzJ8/HwMGDBAdhXQC5xwrV66s37Rp036VSjXPVAsqiy6mrK2tY/z9/b89duyYm7u7u+g4xAJkZWVhy5Yt0Ol0oqN0iEQiwcSJEzFx4kRqjSJGp9fr8dNPPyE5OdkstqeRSCSYNm0aoqKiREchncA5x3PPPdewfv36HSqVarEpdvlZbDHFGIvy9/ffdfToUWqRIkaRlpaG3bt3m81UcTc3N8yZM4dWMCc9rrKyEt9//z0uXbokOkqHREREYMaMGTT5woxwzvHoo4/Wb926dUN5efnjplZQWWQxxRgb7ufndyAlJcWTmnSJMRw8eBBHjhwRHaPDRo4ciWnTpsHa2lp0FNKHnDx5Ej/88INZLFrr5+eH+++/H3K5XHQU0kF6vR4PPPBA7cGDBz9SqVQvic5zLYsrphhjQT4+PkcOHTrkRSvhku7inGPnzp04ceKE6CgdYmdnh1mzZtFWMESYhoYG7Nq1yyz2pXRycsLChQtBvRfmQ6fTISEhoeb48eNvVlRUvCU6z1UWVUwxxvp5eXkd27Nnj29YWJjoOMTM6fV6bNu2DRkZGaKjdIivry8WLFgAZ2dn0VEIQWZmJnbu3GnyrVTW1taYO3cugoODRUchHaRWqxEXF1eTmZn5YmVl5RrReQALKqYYYwqFQpG2ffv2AWPGjBEdh5g5nU6Hb7/91iw+XQNAdHQ0pkyZQmvpEJNSV1eH7du3Iy8vT3SU22KMYdq0aYiOjhYdhXRQc3MzJk6cWJ2bm/t4XV3dV6LzWEQxxRiz9fT0/HnDhg1DpkyZQlOWSLdotVp8/fXXZrGqOa3yTEwd5xxpaWnYv38/NBqN6Di3FR0djWnTpoE2vjcPdXV1iIqKqszNzZ3BOT8uMovZF1OsrUlq6x//+Mepv//9721F5yHmTaPRICkpyeQ/SQOAUqnEggULQMt+EHNQXl6OLVu2oKSkRHSU2woNDcXcuXNppp+ZKCgowNixY0tLSkpGc84LReUw+2LK09Pz1XvvvfeZzz77jJa2Jd2i0WiwceNGs5jePWTIECQmJsLKykp0FEI6TKfT4cCBAzh69KjoKLfVr18/PPDAA7CzsxMdhXRASkoKnzNnzrny8vKRnPMmERnMupiSy+Wzhg0b9tmPP/7oSm8qpDvMpZBijGHSpEmYNGkSdUUQs5Wbm4tt27aZ9OB0T09PLFmyhLagMROffvqp+qWXXkpWqVRTRaxBZbbFFGNsmL+//6ETJ054UDcH6Q6tVouvvvrK5Lv2rKysMHv2bBofRSxCTU0NvvnmGxQVFYmOcksuLi5YsmQJdaWbiSeeeKL+66+/XqNSqZ7r7XObZTHFGPNQKpWnDh065BcaGio6DjFjOp0OmzZtMvnB5rQeDrFEOp0Oe/fuRVpamugotySXy7FkyRIolUrRUcgd6HQ63H333dWnTp3q9Rl+ZldMMcYkCoXip7Vr146eMWMGzQMnXabX6/HNN98gJydHdJTb8vX1xcKFC2mlZmKxMjMzsWPHDqjVatFRbsrOzg6LFy+Gr6+v6CjkDurq6hAREVFx8eLFuzjnub11XrNbRsDNze35+fPnD6VCinQH5xzbtm0z+UJq8ODBWLp0KRVSxKINGzYMK1asgKurq+goN9Xc3Iz169ejoKBAdBRyB05OTti8ebOHQqHYwRjrtf20pK+99lpvnavbGGPhAwYMWLN161YXmrZKumPXrl04deqU6Bi3NXLkSMyZM4emaJM+wcHBAWFhYSguLkZNTY3oODfQ6XTIyspCv379TLboI228vb2h1+vtMzIyvJ999tmdvXFOs+nmY4zZKxSKrOTkZH/ad4x0x6FDh5CcnCw6xm1NmjQJkydPFh2DkF6n1+uxe/du/Pzzz6Kj3JSVlRUeeOABBAQEiI5CbkOv12PChAnVP//8831qtfqHnj6f2XTzKRSK/77++uveVEiR7khLSzPpQkoikeDee++lQor0WRKJBDNmzEB8fDwkEtN7i7q6jEp+fr7oKOQ2JBIJNm/e7Orh4bGOMebR4+fr6RMYg4ODw5zw8PDpv/3tb21EZyHmKysrC7t37xYd45akUinmz5+PyMhI0VEIES4qKgqLFi2Cra3pbWxxtaC6fPmy6CjkNry9vfHxxx8rFArFZtbDC/OZfDHFGFO6urr+e+PGja60SCHpqvz8fGzduhWm2q0tk8lw//33g5b6IOR/AgMDsWzZMpNcOFOtVuPLL7806XWyCDBz5kzpvffeG+7q6vp4T56n08UUY2wtY0zFGMu85rrXGGNFjLH09kv8NT97hzF2gjE2qSsBvby8/vvRRx+506JppKvKy8uRlJQErVYrOspNWVtbY9GiRQgKChIdhRCTo1Qq8fDDD8PT01N0lBu0trZiw4YNKCsrEx2F3MaHH37oLJfLX2GMKXrqHF1pmfocwLSbXP9Pznl4+2UXADDGrg5wmgig01WhTCabFBISMj4hIYGWQSBd0tDQgC+//BItLS2io9yUjY0NFi9eTINZCbkNZ2dnLF++HP369RMd5QbNzc344osvUFVVJToKuQUHBwd88MEH7kql8j89dY5OF1Oc8yMAOvpbIwWgB8ABdKqPjjFm5e7u/vnatWvdOhmREABt4xq++uork5xmDQC2trZYsmQJ+vfvLzoKISbPzs4ODz74oEm24DY0NOCLL75AQ0OD6CjkFmbPni0NDg6ewBib0BP3b8wxU08wxjLauwFdAYBzngXAHkAKgH935s7c3d1f/N3vfqegT+ykK64uymmq4xlsbGzw4IMPws/PT3QUQsyGlZUVFi5ciBEjRoiOcoPq6mps2LABra2toqOQm2CM4bPPPnNTKpXrGGNWxr5/YxVT/wYQCCAcQAmA967+gHP+JOd8FOf8YEfvjDHWz8nJ6ckXX3zR3kj5SB+TnJyMrKws0TFuytraGosXL4aPj4/oKISYHYlEgsTERIwcOVJ0lBuUlpYiKSkJOp1OdBRyEwMHDsQjjzyicHd3f97Y922UYopzXsY513HO9QD+CyCqO/enVCrXffzxx242NrQSAum8rKwsk11LSiaTYeHChSY59oMQc8EYw8yZM01yGZFLly5h+/btomOQW/jjH//oIJfLVzLGjNotYJRiijF27Vb2iQAyb3VsB+5rQlhYWFhcXJzJL9tATE9paSm2bdtmkksgSKVS3HfffTTYnBAjYIxhxowZGD16tOgoN8jIyMChQ4dExyA3YWtri9WrV7t5eXn9y5j32+lNvxhjXwGIAeDBGCsE8CqAGMZYONoGmucD+G1XA3l5eX34z3/+kwadk05rbm5GUlISNBqN6Cg3kEgkmDNnjkkOniXEXF0tqACY3PYzycnJcHNzQ1hYmOgo5Ffi4+Mlnp6e4xljgzjnF4xxn12ZzbeQc+7NObfinPtxzj/lnC/hnA/nnI/gnCdwzku6EoYxNmnUqFH9hwwZ0pWbkz6Mc47Nmzeb7My9GTNmYOjQoaJjEGKR4uPjTbLL77vvvsOVK1dExyC/whjDu+++6+7t7f3enY/uGJPqSvPy8vrg7bffplYp0mmHDh1CXl6e6Bg3FRMTg1GjRomOQYjFutpCFR4eLjrKdbRaLTZt2oS6ujrRUcivxMXFMU9Pz3GMMaN0F5hMMcUYi4mMjOxHrVKks86dO4cff/xRdIybioiIQExMjOgYhFg8xhgSEhIQEhJy54N7UUNDg0nvwNBXXdM69Q9j3J9JFFOMMebl5fXBO++8Q61SpFNqampMds+9QYMGYebMmaJjENJnSCQSzJs3D/7+/qKjXKe4uBjfffed6BjkV+655x6mVCrHMMaCu3tfJlFMAZgUHR3tZ2qfKIhp0+l02Lx5M5qbm0VHuYGXlxfmz58PicRU/sQI6Ruubhru5eUlOsp1Tp8+bXKD5Ps6xhjee+89D2OMnTKJV3ofH5/XX3vtNWqVIp1y8OBBFBYWio5xA0dHRzzwwAOgddIIEcPW1haLFi2Ci4uL6CjX2bNnj8nuytBXTZ48GY6OjlHd3QRZeDHFGPN2cXEJNbWBg8S05eXl4ejRo6Jj3ODqp2InJyfRUQjp0xwdHbF48WLY2dmJjmKg0+nwzTffmGRrel/FGMPKlStdXF1dH+3O/Qgvptzd3R97+umnXUXnIOajsbHRZMdJzZw5E76+vqJjEEIAeHh4YMGCBSbV3V5TU0MrpJuYJUuWWFlbWz/CGOvyL4rQ3zDGmEQmky1fuHBhpxcPJX3Xjh07THJ39rFjx9ICfYSYmICAAMTHx4uOcZ2zZ88iLS1NdAzSztHREbGxsQ4AJnf1PkSX63HTpk2zd3BwEByDmItffvkFubm5omPcICAgAHFxcaJjEEJuIjIyElFR3doy1uj27duHsrIy0TFIu6efftrVx8fnxa7eXmgx5ePj88enn37atEYIEpNVXV2NvXv3io5xA2dnZ5q5R4iJmzZtmkktmaDVarFlyxZaf8pEREZGQi6Xj+jqQHRhr/6MMU8nJ6cQ6hYhHcE5x7Zt26BWq0VHuY5UKsX8+fNhb28vOgoh5DYkEgkWLFgAZ2dn0VEMysrKcPDgQdExSLvf//73Li4uLiu6clthxZS1tfXMxYsXO4o6PzEvx48fR0FBgegYN4iNjYWfn5/oGISQDrC3t8eCBQsglUpFRzE4duwYLl++LDoGATBv3jwrOzu7RV25rbBiSqFQPJSYmGgr6vzEfFRXV+PAgQOiY9xg8ODBGDt2rOgYhJBO8PX1xZQpU0THMLja6q7RaERH6fOUSiWcnZ09GWOdXvdSSDHFGLMGEBoaGiri9MTM7Nixw+ReaBwdHTF79mwwxkRHIYR0UnR0NEzp/aeqqoq6+0zEggULHKRS6fTO3k5Uy9SkqVOnWtEbEbmTU6dO4dKlS6JjXIcxhjlz5tA4KULMWEJCgkmNn0pNTaXV0U1AYmKivZeX19LO3k5IMeXj47P4vvvuo1l85LYaGxuxb98+0TFuMG7cOAQEBIiOQQjpBjs7O8yZM8dkZuFyzvHdd99Br9eLjtKnhYWFgXMexhiz6sztev23iDHGdDrdPRMnTuztUxMzs3fvXpPbdkGhUGDy5C6v60YIMSEDBgzAXXfdJTqGQWlpKVJTU0XH6NMYY4iNjZUCGN+Z24koyQcGBwdb0Saw5Hby8/ORkZEhOsZ1JBIJEhMTIZPRgv2EWIpJkybBy8tLdAyDw4cPo66uTnSMPm3evHluSqVybmduI6KYipw8ebJcwHmJmdDr9di1a5foGDeYMGECvL29RccghBiRVCrF7NmzTWa5BLVabZLDG/qS0aNHQyaTdarJsteLKW9v77vHjBljOtt4E5OTlpYGlUolOsZ1lEolqGuaEMvk5eVlUt19mZmZyM/PFx2jz/L29gbn3Id1YpaciDFTY0eNGtXbpyVmorGxEYcPHxYd4zoSiQSzZs0ymU+uhBDjmzhxIjw9PUXHMNi7dy8456Jj9Fntk4z6d/T4Xi2mWBtlb/7Cnj17FmPHjoWNjQ3efffd637m7++P4cOHIzw8HJGRkYbrq6qqEBcXh6CgIMTFxaG6uhpAW/fTgw8+iHHjxiErK6vXsl65cgWTJ09GaGgohg4dig8++EB41luprq5GYmIiRowYgaioKGRmZhp+tmfPHgQHB2PQoEF48803DddnZWVh7NixeOihh3DgwAG0tLT0Wt5byc/Px5o1a/DRRx/hm2++gY+PD4COPQZTmo2zfft2jBgxwvA7npKSYvjZunXrEBQUhKCgIKxbt85w/eHDhxEZGYnnnntOROQOudnzYKrPwc20tLQgKioKYWFhGDp0KF599VUAwKVLlxAdHY2goCDcd999hu2TGhoakJCQgLvvvhvFxcUio9+gpqYG8+bNQ0hICEJDQ3Hs2DGTe126kw8++ABhYWH44IMPDAPAm5ub8cUXX2DVqlX44osvDJNhOOfYunUrPv300x5tQS8pKcHp06dvuH758uVQKBQYNmyY4bpnn30WISEhGDFiBBITE1FTU2P42RtvvIFBgwYhODj4ur1Nk5KSMHLkSLz//vs9kv9mOV977TX4+voiPDwc4eHh1w3nEJXzViZOnOgAoMMtP73dMhUQFBTUqyd0c3PDhx9+iP/7v/+76c8PHTqE9PR0nDhxwnDdm2++idjYWJw/fx6xsbGGF+t9+/YhOjoaW7duxXvvvddrWWUyGd577z3k5OQgNTUVq1evRnZ2ttCst/L3v/8d4eHhyMjIwPr167Fy5UoAgE6nw+OPP47du3cjOzsbX331leEx/OMf/8COHTsQHByMb775ptey3kpLSwt27tyJhQsX4oUXXjD8wXfkMURGRprUeIfY2FicPn0a6enpWLt2LVasaNt2qqqqCq+//jqOHz+OtLQ0vP7664Y3vH//+9/48ccfodPpcPbsWZHxb+pWz4OpPgc3Y2Njg4MHDxqemz179iA1NRXPP/88nn76aZw/fx6urq749NNPAQAbNmzAb3/7W3zwwQf48MMPBae/3sqVKzFt2jScPXsWp0+fRmhoqMm9Lt1OZmYm/vvf/yItLQ3Z2dkoKSlBZWUlUlJSEBAQgCeffBIBAQGGDyJ5eXnw9fXFfffdh2PHjvVotoMHD96wEfLSpUuxZ8+e666Li4tDZmYmMjIyMHjwYLzxxhsAgOzsbCQlJSErKwt79uzBY489Bp1OB6CtSPn555+RmpqKhoYGo2e/WU4AePrpp5Geno709HTEx8cLz3krY8eOdVAqlTEdPb63i6lREydO7NXB5wqFAqNHj4aVVceXjNi+fTseeughAMBDDz2Ebdu2AWh7EZdIJJBIJD3S/HqrrN7e3hg5ciSAtpW3Q0NDDYu7icp6K9nZ2YiNjQUAhISEID8/H2VlZUhLS8OgQYMwcOBAWFtb4/7778f27dsNWRljyM3NNYkWhTNnziA0NBTOzs6YOnUq+vXrBwAdegy9/f99J3K53LBKe2Njo+HrvXv3Ii4uDm5ubnB1dUVcXJzhhU+v15vkY7nqVs+DqT4HN8MYg1ze9lKo0Wig0WjAGMPBgwcxb948AKbx93wndXV1OHLkCB5++GEAgLW1NVxcXEzudel2cnJyMGbMGNjb20Mmk2HWrFnIy8tDbm4uwsLCALStPZSbmwvgf38fvbHodF1dHY4fP37ddRMnToSb2/W7nUyZMsUwy3jMmDEoLCwE0Pb+cP/998PGxgYBAQEYNGgQ0tLSAMDw/88Y65Hn4mY5b0VkzlsZNWoUZDJZh5dH6NViytPTc3R4eLjJLBvNGMOUKVMwatQo/Oc//zFcX1ZWZpi15e3tbWjKnTp1KpKTk5GQkIBnnnlGSOb8/HycOnUK0dHRJpk1LCwMW7ZsAdD2pldQUIDCwkIUFRUZihIA8PPzMxSEK1euxD333INTp04hMDCw17LeSmVlJZqbm7Fp0yY8+OCDWL9+PQDc8THMmDEDx44dM6l9vwBg69atCAkJwYwZM7B27VoAt38sK1aswLhx46DX601qy42rbpXdlJ+Dm9HpdAgPD4dCoUBcXBwCAwPh4uJieFO89jlZtGgRPvzwQzzxxBN48sknRca+zsWLF+Hp6Ylly5YhIiICK1asQGNjo8m9Lt3OsGHDcOTIEVRWVqKpqQk//PADHBwc0NDQAEdHRwBtH2IbGxsBAIMGDUJBQQGSkpIwZsyYHs+XkpLSqaEPa9euxfTpbbuh3O7vfM6cOYiMjERkZKThcfaGf/3rXxgxYgSWL19uaA03xZze3t7Q6/Udnr7dqwvm2NraBvr6+vbmKW/rp59+go+PD1QqFeLi4hASEnLbGVsymQxJSUm9mPB6DQ0NmDt3Lt5//304OTnd9lhRWV944QWsXLkS4eHhGD58OCIiIiCTyW76ieLqJ7uIiAg8/vjjJrOVgl6vR2lpKVJSUuDg4ICxY8dizJgxd3wMv/4EaSoSExORmJiII0eO4JVXXsH+/ftv+1imTp2KqVOn9nbMDrtVdlN+Dm5GKpUiPT0dNTU1SExMRE5Ozg3HXH1OXFxcsHv37t6OeEdarRYnT57EqlWrEB0djZUrV143lvDXRL+G3kxoaCief/55xMXFQS6XIywsDLa2trdcGV0ikRhaD3tDc3MzUlNTERMTc8dj//a3v0Emk2HRokUAbv23ArS1GF5tPewtv/vd7/DKK6+AMYZXXnkFf/jDH7B27VqTy3n1/NbW1jLGmIRzfscuk15tmeKc+10dyNuTVq9ebRjgdrvBmlezKBQKJCYmGpoVlUolSkpKALQNAlQoFMKzajQazJ07F4sWLcKcOXMM1/dm1lu59jE0NDTgs88+Q3p6OtavX4/y8nIEBATAz88PV65cMdymsLDQ8P+fk5MjvJBKS0vDmjVrsGbNGjg6OmLSpEkYMGAAPDw8MHHiRJw+ffq2j8GU3Op3auLEicjLy0NFRYXZPJabMefsN+Pi4oKYmBikpqaipqbGMEbGHB6Xn58f/Pz8DC3l8+bNw8mTJ03idakzHn74YZw8eRJHjhyBm5sbBg8eDC8vL9TX1wMA6uvr4eDgICxfamrqHVun1q1bh++//x5ffvmloRAxtb8VpVIJqVQKiUSCRx55xPCea2o5r/Lw8OAA3DtybK8WU1qtVqFUKnv8PI8//rhhgNutnpDGxkbDH8rVPeCuzjpISEgwzG5at24dZs2aJTQr5xwPP/wwQkNDb2ga782st3LtY7C3tzfMQPrkk08wceJEODk5YfTo0Th//jwuXboEtVqNpKQkJCQkgHOOQ4cO9XrmX4uKisKjjz6KRx99FOHh4SgqKoJWq0VTUxOOHz+O0NDQWz4GU3Pt89HU1GT41Hfy5Emo1Wq4u7tj6tSp2LdvH6qrq1FdXY19+/aZdGvUtczlebid8vJyw4yr5uZm7N+/H6GhoZg8eTI2b94MQNzfc2d4eXmhX79+hvFEBw4cwJAhQ0zidakzrnZDXr58GVu2bMHChQsxd+5cwxv86dOnERwcLCxfS0vLbVtd9+zZg7feegs7duy4bgP2hIQEJCUlobW1FZcuXcL58+cRFRXVG5Fv6mqBDbQNP7j2PdeUcl7l5+fHAHSoq69Xu/kkEoltZwaCG0NpaSkiIyNRV1cHiUSC999/H9nZ2aioqEBiYiKAtqbqBx54ANOmTQPQ1lW1YMECfPrpp+jfv3+vzTC7VdaMjAx88cUXhmUcgLZZc/Hx8cKy3kpOTg4efPBBSKVSDBkyxDAbSSaT4V//+hemTp0KnU6H5cuXY+jQoThz5ozJLdA5d+5ceHl5YcSIEZBIJFixYoXhj/5mj8GUffvtt1i/fj2srKxgZ2eHTZs2gTEGNzc3vPLKKxg9ejQA4E9/+lOHB4uKdqvfJXNSUlKChx56CDqdDnq9HgsWLMC9996LIUOG4P7778fLL7+MiIgIw8BuU7Zq1SosWrQIarUaAwcOxGeffWZ4TKbyunQnc+fORWVlJaysrLB69Wq4urrihRdewKxZs3DgwAE4Oztj/vz5QjMeP37csPTH4cOHDS3Mr7/+Ot544w20trYiLi4OQNsg9DVr1mDo0KFYsGABhgwZAplMhtWrV/faenkLFy68Iefhw4eRnp4Oxhj8/f3x8ccfA4DQnLfj7+9vg7Zi6o57m7HeGh3PGJP6+/uXXrp0yaNXTkhMHuccq1evRkVFhegoBi4uLnjiiSdo/z1CCADgq6++MrS8iTZt2rReGfRO2rz//vvap59++jec88/udGyHuvkYY9MYY7mMsQuMsRfarxvKGDvGGFvHGOvI/XgqlUrx896JycjKyjKpQgoAYmJiqJAihBh0ZOB3bzl27JhJLB/TV/j6+src3d0HduTYOxZBjDEpgNUApgMYAmAhY2wIgGcAJAA4AaAj85BdzKUbgfSOa1fjNgWenp6GdWUIIQRomyIvcrzUtWpraw0LBZOe5+bmBhsbG6+OHNuRFqUoABc45xc552oASQBmAZAC4AD0ADqyepmVtbV1z69yRszC+fPnUVpaKjrGdSZNmtQrC/ERQszLhAkTREcwMKflP8ydTCaDRCLp0EDvjhRTvgCuXPN9Yft1HwDYCWAsgI7s3SCjYopcZYqtUuY2iJkQ0jv8/Pzg7+8vOgaAtr1ar50VR3pO+4Q5oxVTNyuAOOf8FOc8mnO+mHOu68D9yGgsCgHaVrstKCgQHeM6EyZMoFYpQsgtjRs3TnQEg2v3kiU9p31GoWmTOUYAACAASURBVNGKqUIA/a753g9AV7Yt1/16w0bSN/X05qCd5erqet3O5oQQ8mtBQUEms3zImTNnDOv5kZ7TvtlyhwqXjhRTPwMIYowFMMasAdwPYEcXcmm0Wq1p7GxJhKmrqzO5AZTjxo275dYRhBACtG0vYgoLSQKAWq1GZmam6BgWr70BSNORY+/4DsI51wJ4AsBeADkAvuacZ3Ull0bToUzEgp04ccKkpvba29sbFkIlhJDbCQ8PR28vPH0r6enpoiNYPK1Wi/aJd3fUoY/jnPNdnPPBnPNAzvnfupirvrq6mlqm+jCdToeTJ0+KjnGdyMhIk3lxJISYNltbW5OZqHL58mVUVVWJjmHRamtroVarO/Sf3Jt9G2WlpaU0wrcPy8nJQUNDg+gYBhKJBJGRkaJjEELMSEREhOgIBtTV17OKi4v1FRUVFzpybK8VU5xzTUtLS0dm/RELZWozUEJDQ+Hk5CQ6BiHEjAwYMACurq6iYwBo20WC9Jz8/PxGznmHJtz16qhbrVarbh8dT/qYqqoq5Ofni45xHWqVIoR0xfDhw0VHAACUlZWZ3JZcliQ/P78VQIcW9erVYkomk1WUl5f35imJiTh16pToCNdxd3dHQECA6BiEEDNkKsUU0DZ8gvSMK1eu6GGKxZREIiksLu7KElXEnHHOkZGRITrGdUaOHCk6AiHETHl6esLT01N0DADA2bNnRUewWCqVigFQdeTYXi2m1Gp1PhVTfU9BQQFqa2tFxzCQSCQYMWKE6BiEEDMWGhoqOgIAoLi42KQm9liS1tZWXQd3eOndYqqiouKXrKys1t48JxHv9OnToiNcZ+DAgXB0dBQdgxBixkJCQkRHANDW8n/hQocmnJFOqKysBOe8w2tP9GoxpdPpfklOTq7rzXMSsXQ6ncn16VOrFCGku7y9vSGXy0XHAABcvHhRdASLc/LkSTDGUjt6fG/voXGW1sXoWy5cuICWlhbRMQysrKxM5hMlIcR8McYQGBgoOgYAKqZ6wvHjx1tKSkoOdPT4Xi2mOOc6rVZbZUrjZ0jPMrV9+AYNGgRra2vRMQghFmDgwIGiIwAAGhoaQDPljSs5Obmec/5LR4/v9d1dJRLJz6Y2TZ70DL1ej3PnzomOcZ0hQ4aIjkAIsRCmtLyKqa3jZ+7aZ0me7+jxvV5MlZaWHkhLS6NB6H1AQUEBmpubRccwkEqlCAoKEh2DEGIhnJycTGY19CtXroiOYDGqqqqg1+srOOf6jt6m14spnU534vDhwzQIvQ8wtfVPBgwYAFtbW9ExCCEWpF+/fqIjAAAKCwtFR7AYv/zyCwB0ePA5IKCYApB96tQpTtvKWD5T6+IbPHiw6AiEEAvj5+cnOgKAttYUU+oJMGc7duyoKy4u3tKZ2/R6McU510ul0mNpaWm9fWrSiyorK1FdXS06xnUGDRokOgIhxML4+PiIjmBQWloqOoJF+P7771sBdHgmHyCmZQpFRUXrNm/eTEu2WrC8vDzREa7j7OwMDw8P0TEIIRZGqVSCMSY6BgCgpKRD28iR27hw4QI0Gs1FznmnmvmEFFMAfti2bRu1R1owU1v3xJRm3RBCLIeVlRXc3NxExwAAlJWViY5g9rZv395aXV29vrO3E1JMcc4bWlparhQUFIg4Pelher3e5Kbp+vv7i45ACLFQprLpMa011X0bN26sa2pq2t7Z24lqmUJtbe2G7du3a0Sdn/SckpISk1r1HGibyUcIIT3B3d1ddAQAQEVFhegIZq22thbFxcX1nPOizt5WWDHV2Ni4Zf369TWizk96jqm1ODo6OprMWjCEEMtjKt18arUajY2NomOYrf9v787jo6rP/YF/vrNlmWQy2RcWQdDIvhiUUGWRooKKCoiy2KpFbEtRBK/1er21t/5qrRW3yqvaAi4tVwVFRBQFFIXKLoQ1QCKQhWTmTGZJJjNJZnt+f2SSm7Bmmcw5M3ner9d5nWGWcz4oSZ58z/c8388//zzg9/vbdRdfE9mKKSIqqaioqDh9+rRcEVgXKS0tlTtCK0q5dZkxFp2MRqPcEZo5HDxG0VFLly61WiyWv3Xks7IVUwBQVVX1whtvvMF39UUZpTWPU9Kty4yx6JOUlCR3hGZKa0kTKYqKilBZWVlKRB26e0rWYsrr9X7y/vvvu71enjoVLaqrq1Fbq6z6mIspxlhXSkxMlDtCs+rqarkjRKS//vWvNRaL5YWOfl7WYoqIGgKBwNr169e3ef0bpmxnz7Z73l6Xy8rKkjsCYyyKxcTEQKPRyB0DAOB0OuWOEHEaGhqwZs2aOp/P1+67+JrIWkwBgNlsfvkvf/mLVe4cLDSUVkzp9Xro9Xq5YzDGolxcXJzcEQCAJ6B3wCeffBLw+/2riajDl8lkL6aIqKi0tJQnokcJpXXgzcjIkDsCY6wbUMoi6rw+X/u99NJLNovF8kpnjiF7MQUAVqv1j3/605/4Qm8UUFoHXqX0f2GMRbeYmBi5IwBovGTF2q6goAAVFRXFRNSpER1FFFMej+fj9evX2yoqKuSOwjrB7XYrboiZiynGWDgoZc6U0homK92SJUuslZWVizp7HEUUU0QUcDgcTz377LM8OhXBlNh9l5t1MsbCQa1Wyx0BAI9MtceBAwdw7NixYiLa3dljKaKYAoCGhoaPPvvsM5vSJjCztlPiulBKaqbHGIteQgi5IwAAfD6f3BEixpIlS6wmk+nRUBxLMcUUEQWqq6t/y6NTkUuJI1MGg0HuCIyxboCI5I4AAPD7/XJHiAj79+/H8ePHi4hoTyiO165iSghxjRBipxCiQQjxxDmvnRFCHBZCFAgh9rV4PkcI8Y0Q4lMhRMKljl9fX//xhg0brErroM3axmazyR2hFY1Gg/j4eLljMMa6gUCA2yVGkuBcqZCMSgHtH5myAXgUwEsXeX0CEQ0norwWzz0KYCGA5QDmXurgRBSw2+1PPvPMMzw6FYGUtowBF1KMsXBRyuU1pVxuVLK9e/fi+PHjJ4lob6iO2a5iioik4Mnb09hKDSAQ3C77f9nj8azduHFjWUFBQXuiMQWoqamRO0Ir3KyTMRYuSlkWjYupSwsEAnjooYesJpNpfiiPG8o5UwRgkxDiByFEy5BvAHgLwC8B/OuyByEiSZLm/PznP7fysGnk8Hg8irslVykdiRlj0U8pd9FxMXVpf//73z2SJH1EREdCedxQFlM/IaKRACYDWCCEGAsARFRCRGOJ6A4iatOiQUR0yGw2r1uxYoUnhPlYF1Li4ppKaaLHGIt+Suk8rpR+V0pUVVWF//mf/7FKkvTE5d/dPpctpoQQC4KTyguEEDkXex8RVQT3EoBPAFzXmWBms3nx7373O5vVysv2RQKlXeIDAJ1OJ3cExlg34PP5FDMyz9/3Lm7hwoWO6urqx4moNtTHvmwxRUTLgpPKhzcVTOcSQuiFEIlNjwHcDKBTQ2hEVON0OpcsWrTI0ZnjsPBQWudzANBqtXJHYIx1A7W1If/Z3GFKWSNQaXbu3ImtW7eerKurW90Vx29va4QsIUQ5gMUAnhFClAshDAAyAfxbCHEQwB4AnxPRl50N53K53t+0aVPxrl27Onso1sWUMsTdkkqlmDZqjLEopqSReb6L+Xw+nw8PPPBAldlsnktd1BCsXRdXicgEoOcFXqoBMCwkiVqfj4QQc+6///7vDx8+nMYVt3IpZYi7JS6mGGPhoKQ5o3zjzfn+9Kc/uex2+ztEVNRV51D8TxsiOmmz2V549NFHlVP6s/MocWSK72phjIWDkhoW88hUa3v37sUbb7xRYrFYnu7K8yi+mAIAm8328qeffnro888/5z75CqXEkSmlLO/AGItuSrpRKiHhkguNdCu1tbWYOXNmlSRJU4moSxuBRUQxFew9NW3+/PkWs9ksdxx2AUrpsdISr1HFGAsHSZLkjtCM1yP9P/PmzXPYbLaniOjHrj5XRBRTAEBEFqvV+uA999xj52aeyqOUpRRaUmImxlh0CQQCilrkPTk5We4IivDBBx94t27duqumpmZlOM4XMcUUANTX13954sSJj1955RXlTdDp5pRY4Ho83POVMda1LBaLon5xS0lJkTuC7EpLS/H4449LkiTN6qq7984VUcUUAEiS9JsXX3zx7IEDB+SOwlpQYjGlxEuPjLHoUlFxwfaLsoiNje32E9C9Xi/uvvtum8ViuY+IwtanMuKKKSJqkCRp8tSpU3n+lIIo6TezJkqcFM8Yiy7l5eVyR2jGo1LA/Pnza0pKSl7z+Xz/Dud5I66YAgAiKjabzbNvvfVWG48+KIMSR6aU2K6BMRZdSktL5Y7QLD09Xe4Isnr99dfrv/jii2+sVutz4T53RBZTAODxeLaUlpb+ce7cudV8Czy7ECUt8cAYiz4ulwsWi0XuGM0yMzPljiCbLVu2BP74xz8WSZJ0X7jmSbUUscUUAFit1pe3bdu24YUXXnDLnaW7U6vVckc4T0NDA09CZ4x1mVOnTskdoZXs7Gy5I8iiqKgI999/v0mSpElEJMvlqogupgBAkqQHX3nllcINGzZwUyEZKXVRYSUt88AYiy4//tjl7YvaTAiBnJwcuWOEncPhwK233mo1mUxTiEi2idQRX0wRkddisdzyi1/8ouLIkSNyx+m2lLoOnsMRtps5GGPdCBGhuLhY7hjNUlNTERMTI3eMsPL5fLjjjjvskiT9iogOyplFmT8B24mIrJIk3XzrrbdKp0+fljtOt8TFFGOsO6moqFDUvMzuNioVCAQwa9as6hMnTrzudDrXyJ1HmT8BO4CIjp89e3bKhAkTqpTU96O7UOplPiWtmcUYix6FhYVyR2ilV69eckcIGyLCww8/XPPdd9+9K0nS7+XOA0RRMQUARPRDeXn53ePGjbMqqb1/d6DU4WUl3WnDGIseR48elTtCK3369JE7QtgsWbKkdsOGDessFssiubM0iapiCgB8Pt+/y8rK5o4bN87Gk4/DR6ldd5W0ACljLDqUlZXBbrfLHaOZwWDoNj2m/vCHP7hWrVq1RZKkB+VogXAxUVdMAY1r+JWWlv5ywoQJdpfLJXecbiEuLk7uCBfkdDoVNa+BMRb5Dh06JHeEVvr16yd3hLB49dVX65YtW7ZLkqR7iEhRnaKjspgCAKfTueb06dP/ccstt9i5S3rX0+l0cke4KJPJJHcExliU8Pl8UNqd43379pU7QpdbuXKl5/nnnz8oSdJtRKS49cuitpgCALvdvqKwsPAPEydO5BGqLqbUkSlAWQuRMsYiW2FhoaKWqhJC4Morr5Q7Rpf6+9//3vDUU08dsVgsE+Vqynk5UV1MAYDVan316NGjT9544412vk2+6yQkJMgd4aKUtBApYyyy7du3T+4IreTk5Cj6+29n/fnPf6575pln9loslhuISLGrnUR9MQUAdrt9eVFR0fz8/HwbT0juGklJSXJHuKiysjIoaJ4iYyxCmUwmlJSUyB2jlQEDBsgdoUsQEZ588snal19+eavFYrmJiJQzHHgB3aKYAgCn0/nRmTNnZo0ePdqqpFW+o4XBYJA7wkXV1dXxvCnGWKft3LlT7gjnGTRokNwRQi4QCGD+/Pk177777npJkqYSkVfuTJfTbYopAKirq9tUWlp6+09+8hPLyZMn5Y4TVVQqlaLnTXFnfMZYZzidTsVNPM/OzkZycrLcMULK5/Nh5syZ1evXr39PkqS5RBQR6+52q2IKAHw+367y8vKJEyZMMBUUFMgdJ6okJibKHeGilLa6O2MssuzcuRN+v7J+rkfbJb76+npMmTLFsW3bttfMZvNCJfWRupxuV0wBABEdrqiouPGWW24p+/zzz5X11RHBjEaj3BEuqqSkBD6f4u6mZYxFALfbrbiJ5wAwePBguSOEjCRJGDNmjH3//v3PSpL0rNx52qtbFlMAQETFkiSNePDBBw+9+OKL7ggqgBUrLS1N7ggX5fV6+VIfY6xDduzYAY/HI3eMVnr16oWUlBS5Y4TE4cOHMWrUKMuxY8dmVVVVvS53no7otsUUABCR1WKxjF66dOmG+++/v1ppXyyRRulf2EpbmJQxpnxOpxO7d++WO8Z5rr32WrkjhMRnn33mnzRpUmlpaenY+vr6r+TO01HdupgCACLySJJ035dffvniT37yE7vVapU7UsRS8pwpADh+/DgCAUWtQMAYU7hvv/0WXq+ybiaLiYmJ+Lv4iAgvvPBC3bx58wrMZvMIIjoud6bOUHQxJYS4UwhxSAhRIITYJ4S4ocVrPxdCFAW3n7d4fnzwvS+29TxERFVVVc8fPXr0Z3l5eVU8gtExSr+rxO1286U+xlibVVVV4cCBA3LHOM+gQYOg1WrljtFhHo8Hs2fPrn7llVc+lSQpn4hscmfqLEUXUwC+BjCMiIYDeAjAcgAQQqQAeBbA9QCuA/CsEKLpJ/mvANwIQC2EuKY9J3O73RvOnDkzfsKECeUbN27kIYx2Sk1NhUql7H9SR48elTsCYyxCbN68WZGj2SNGjJA7QodZLBaMGTPGvnnz5j9JkjQ7EnpItYWif/IRUW2LWyP1AJoe3wJgMxHZiMgOYDOAW4OvqYLvCwAQHTjnUbPZPPKBBx7Y/9RTT9Uq7VZYJVOr1YqehA4Ax44d47v6GGOXVVxcjBMnTsgd4zzZ2dno1auX3DE6ZPv27TRy5Ejp2LFjc6qqqv4cSa0PLkfRxRQACCHuFkIcB/A5GkenAKAHgLIWbysPPgc0jl7tAKAiog5dryMiiyRJo99+++2/5efn2ysrKzuYvvtRejFVX1/PE9EZY5fk8/mwceNGuWNc0OjRo+WO0G6BQADPPfece8aMGUfKy8vz3G63Mv/jdoLiiyki+oSIrgFwF4Dngk9faMSJgu//iohGEtGSTp7XbzabnywoKLgnLy/PtHnzZuWN9SqQkntNNVFivxjGmHJs374dSrwZyWAwRFxvqaqqKowfP97+xhtvvCNJ0rVEVHb5T0UexRVTQogFwQnnBUKInKbniWgbgH5CiDQ0jkS1HOfsCaCiK/J4PJ6vKyoqRsydO/cAX/a7vPT0dLkjXFZJSQl4wWvG2IVIkoR///vfcse4oGuvvRZqtVruGG22fft2GjFihLRv376fmc3mBdEyP+pCFFdMEdEyIhoenHQeL4QQACCEGAlAB8AK4CsANwshkoMTz28OPtdVmUySJF2/cuXKv+Xn59v4st/FZWVlyR2hTX744Qe5IzDGFIaIsH79esUtGwMAGo0GeXl5csdokwtc1tsgd6auprhi6hzTARwRQhQAWAbgXmpkQ+Mlv73B7Q9dfWslEfklSXqyoKBgZl5enumLL77gy34XkJGRERG/ORUUFKChoUHuGIwxBdm1axfKy8vljnFBw4YNg16vlzvGZVVWVnaLy3rnUnQxRUR/JqJBwZGqfCL6d4vXVhJR/+D2drgyBS/7DX/ooYe2z54921FTUxOuU0cEtVqNpKQkuWNcVkNDA/bv3y93DMaYQlgsFnz99ddyx7ggtVqNG2644fJvlBERYdWqVd6RI0dW7Nq1675ov6x3LkUXU0pFRGaz2Txhw4YNiwcPHixt2bKFR6laSE1NlTtCmyhxFXjGWPj5fD58/PHHim2bMnz4cEU3RZYkCZMnT7YvWbJkg8lkGujxeDbJnSncuJjqICKimpqat8vKykbOmTNn10MPPVRdW1srdyxFyMzMlDtCm9TU1CiyuzFjLLw2bdoEk8kkd4wLEkJgzJgxcse4qDVr1viGDx9u2rZt2wMmk2kaEVXLnUkOXEx1EhGdlSTphk8++eQ/Bw8ebNm2bVvUNCHrqB49elz+TQrx/fffK7LDMWMsPAoLC7Fnzx65Y1zUgAEDFDnab7PZcOedd9oXLly4qbKycpDb7V4vdyY5cTEVAkREdrv9byUlJaPuueee/b/+9a9r6urq5I4lm0hoj9DEbrfj0KFDcsdgjMnA4XDg008/lTvGRQkhMG7cOLljnGfDhg3+IUOGmLds2fKIyWS6LRrW1ussLqZCiIhKJEm67oMPPvj9gAEDLN11fb9ImIDe0nfffcdzpxjrZnw+H1avXo36+nq5o1zU0KFDFTVtoqKiArfffrt93rx5WysqKoa6XK41cmdSCo3cAaINEQUAvCKEWP3ggw+uHDZs2Kjly5cnR+paSh2h1WqRlJSE6urIuHRut9uxb98+XH/99XJHiXiBQAA+nw8+nw9+vx9+v7/V46YtEAi02hMRAoHAeRsRNb/W9LhpOa+Wjy8l2KoOQojzNpVK1erxhTa1Wt28b/lYo9Gct9doNIpf7Js12rBhAyoquqTXc0io1WrcdNNNcscA0Fh4vvrqq/UvvfSSxWq1PuL1eqNuOZjO4mKqixDRWQC36HS6m/Py8pY//vjjaUuWLInTarVyRwuLnJyciCmmAODbb7/F0KFDERcXJ3eULkVE8Hg8aGhoaN48Hs8FN6/XC6/X2+px0+bz+Vo9btp4/hlaFVsajQZarbZ537RpNBrodDpotdpW+3O3mJiY5n3TYy7WOm/Hjh0oKCiQO8YlDRs2TBGj/Dt27MCDDz5YZbfb37ZYLL8jIuUO5clIRNGizYolhIhJT0//fVJS0ry33347Ten9QkJhy5Ytil2S4WKuv/56TJ48We4Yl0RE8Hq9qKura7XV19c378/dGhoamvfcqDTyNRVXsbGxzftzt7i4uOZ90xYfHw+tVts8UtddFRYWYvXq1W0a1ZSLWq3GwoULZV3r1Gq14tFHH3V8/fXXJ81m81wiKpItTATgkakwIKIGAP8phPjHjBkzVo0dO/aaZcuWGSNponZ7xcfHyx2h3fbt24frrrsurHfOEBHq6+vhcrkuuLnd7lZbXV2dYnvhsPBoGjl0Op3t/qxarW4urOLj46HX61s9btoSEhKg1+sRGxsbVcWXyWTCJ598ouhCCgBGjRolWyEVCASwfPlyz+9//3urw+F4vK6ubjUp/T+YAvDIVJgJIUR8fPx0g8Hw2hNPPJGycOHCWJ1OJ3eskDt27BhWr14td4x2u+qqqzBnzpxOHycQCMDlcsHpdJ631dbWNm8ul4snvzPFUqvVSEhIOG9LTEyEwWBAYmIiEhMTodfrFV901dTUYPny5VD6qhXx8fFYuHChLFMOdu7cifnz51stFssnZrN5MRG1v2LvpriYkokQIj4tLe0ZvV4/7y9/+UvKjBkz1Er/ZtQeFosFy5YtkztGh8ycORMDBw686OtEhLq6OlRXVzdvDocDNTU1zZvT6eT5Q6zbUKlUMBgMzVtSUlKrzWg0IjY2VrZ89fX1WLlyJSRJki1DW918881hb9JZXFyM3/zmN7aDBw8eM5lM84moMKwBogAXUzITQmRmZma+kpGRcfObb76ZquROt+3hdDqxdOlSuWN0SEJCAh577DHY7XbYbDbY7XY4HI7mvcPhgMfjkTsmYxElJiYGycnJMBqNSE5Obt5SUlJgNBq7bIF0j8eDf/7znygrU/56u3q9Ho899hjCdbXCarXi6aefrvn0009NFovlEb/f/21YThyFeM6UzIjIDGC2EOKa6dOn/2Po0KEDly1bltK/f3+5o3VKJNxxpFarkZKSgtTU1FZbcnIyCgsLsXbtWrkjMhY1GhoaYDKZLrhsixCiuchq+jpMS0tDamoqkpKSOvz9xOv14v3334+IQgoAxo4dG5ZCqr6+HkuXLq3761//anc6nU+53e5VwbY+rIO4mFIIIjoO4EaNRjPuhhtu+PvUqVOznn/+eUNaWprc0TpEziH9c2k0GqSnpyMjIwPp6elIS0tDeno6kpOTL/pNesiQISgoKMCpU6fCnJax7oeIYLfbYbfbz/ua02g0SEtLa/66bfpaTklJuWSR5ff7sXr1apw+fbqr44dEeno6Ro0a1aXnCAQCWLVqle/pp5+2u93uN2w224vc6iA0+DKfAgkhVPHx8bMTExP//Ktf/cq4ZMmS+ISEBLljtdvzzz8f9sthRqMRmZmZyMrKQlZWFjIyMi5ZNF2K3W7H3/72N76kx5gCNf2SlJmZ2eprPi4uDoFAAB9//DGOHj0qd8w2mzNnDq666qouOTYR4auvvqLHH3/carfb15vN5ieJyNolJ+umuJhSMCFErNFo/E1sbOwTCxcuNCxatCguUloO1NfX44UXXuiy4wshkJycjJycHGRnZzfvQz0iduDAAUWv3cUYay0lJQW//vWv8d133+Hs2bOoqKhQ9JIxAJCVlYVf/vKXIT8uEeHrr7+mJUuWWCVJ+t5kMi0mIh5u7wJcTEUAIURccnLy47GxsY8tXrzYsGDBglild+q22Wx4/fXXQ3a82NhY9OzZs3nr0aNHWG4dJiKsXbsWhw8f7vJzMcY6R6VSYfr06Rg0aFDzc0QEm82G8vJylJWVoby8HJIkKepu21mzZiE3Nzekx9y2bRsWLVpUVVlZucdkMj1ORCdDegLWChdTEUQIoU9JSXkiNjb2V4899ljSggULYvV6vdyxLqiyshJvvfVWhz9vMBhwxRVXoHfv3ujduzcyMjJk62Pj8Xjw1ltvwWrlUXHGlEqr1WLmzJltulTW0NCAsrIylJaWoqSkBGfPnpWtGW6fPn3wwAMPhORYRIRvv/0WTzzxRFVFRcV+k8m0iNschAcXUxFICKFPTk5+LDY29tEFCxYYHnvssTilzan6/vvvsXnz5ja/X6/Xo2/fvs1bSkpKF6Zrv/LycqxcuVJRv80yxhrFxMRg7ty56OiC8j6fD6WlpTh9+jROnTqFioqKsHVJf/jhh9GjR49OHYOIsGnTJvrtb39rlSRpV2Vl5X8Eb2piYcLFVAQTQsQF51QtfuihhxIXLVqkV8oSNe+88w7OnDlz0ddVKhV69eqF/v37o1+/fsjOzlZ8B+X2FoiMsa4XHx+PuXPnIicnJ2THrKurw48//oiioiIUFxfD5XKF7Ngt9ezZE/Pmzevw5/1+G4jeJAAAFq1JREFUP9atWxf43e9+Z7Pb7f8OFlHFIYzI2oiLqSgghIhJSEj4uV6vf/qWW25J+q//+i/j1VdfLVueuro6vPzyy/B6va2ej4mJQf/+/ZGbm4urrrpKluUSOoOIsHr1ahQW8qg5Y0pgNBoxd+5cdGULGSJCeXk5jh8/juPHj4f0cv99992Ha665pt2fc7lcWLFihWfp0qXVDQ0N681m83NEVBKyYKzduJiKIkIIlUqlujUzM/P5a665pudzzz2XOmbMmLCP+Hz66ac4cOAAgMYCKjc3FwMHDkT//v2h0UR2azOPx4OVK1desPEgYyx8cnJyMGvWLCQmJob1vGazGUePHsWRI0dgs9k6fJycnBw8/PDD7fr+bDabsXTp0tr33nuv1uv1vmmz2V4jIkeHQ7CQ4WIqSgkhRubk5PwxOTk579lnn02ZNm2aqquWa2ipoqICK1asQL9+/TB06FDk5uZCq9V2+XnDqaamBm+++SbcbrfcURjrlgYOHIi77rorbMuuXMzZs2dRUFCAI0eOoK6url2fvfPOOzFixIg2vbewsBDPPfec/ZtvvnHU1NQ8V1dXt4qIuAGegnAx1QWEEKMA7AJwLxF9FHzOD6Dp/vpSIpoafH4QgOUATgJ4MNQt/YUQvTMzM/87JibmzsWLFxvmzZsX01V3APp8PhQWFqJv375Q2oT4UPvxxx/xr3/9K2yTVBljjcaNG4fx48crao6l3+9HYWEh9u3bd8m5ok3UajWeeOKJS051ICJs27YN//3f/11VXFxcYjabnw4EApuJv+koEhdTISaEUAPYDKAewMoWxVQtEZ1XYQghVgB4CsB9AIqI6MsuymVMSUl5TKvVPnL77bfHL168OGngwIGdPi4Robi4GOXl5cjPz1fUMjJdbdeuXfjyyy7538UYO4darcYdd9yB4cOHyx3lkiwWC3bv3o2CgoKLtlsYPHgwZsyYccHXampq8O6773pee+21GrfbvauysvIZIjrYlZlZ53ExFWJCiEUAvABGAdjQhmLqHQBPALgXwCki2tjF+dQqlWpyVlbWMxkZGVc++eSTqdOnT1e1d7jc7/fj0KFDKCgowIABA3DttddG3eW8ywku0YBdu3bJHYWxqKbX63Hvvfeid+/eckdpM5fLhT179mD37t3ndWCfPXs2zr1JqKCgAEuXLrVv3rzZ1dDQ8JbD4XiTiKrCmZl1HBdTISSE6AHgfwHcBGAFWhdTPgAFAHwAXiCidcHnRwB4E0ARgJ8TkT+MeftkZGQsVqvV986aNSt+4cKFCX369LnkZ/x+PwoKCnD48GH06dMHY8aMkX3egpyICOvWrcPBg/yLI2NdoWfPnpg5cyYMBoPcUTqkvr4eO3bswO7du9HQ0AAAWLRoEYxGI+rr6/HBBx/4X3rpJbvdbi+sqKj4fwC2hHq6B+t6XEyFkBBiDYClRLQrOOLUspjKIaIKIcSVAL4BMJGIfpQxbjMhhC4mJmZGSkrKU/369cv67W9/mzZ58mRx7oT12tpalJeXIyUlBenp6YqasyAnv9+P9957DyUlfGcyY6E0evRoTJo0CeG4eaaruVwufP311zh58iTuvPNOvPbaazUfffSROxAI/MtisbxGROVyZ2Qdx8VUJwkhFgB4OPjHJABNFUYaADeA+U2jUC0+8w5aFFpKIoQYlJWV9Vu1Wn3rz372s/h58+bpr7zySrljKZ7b7caKFSt4yRnGQkCn02HatGkd6sGkVC6XC2vXrvW/+uqr9srKylKLxfK8z+dbT0Tey3+aKR0XU12kZcEkhEgG4CaiBiFEGoCdAO4komOyhrwEIYRep9NNT01NfSw9Pb3XggULku+9915NUlKS3NEUy+FwYOXKlaipqZE7CmMRKz09HTNnzoRSVnPojEAggG3btmHZsmXW7du3N/j9/v+tqqp6UylXJVjocDHVRc4ppsYAeAtAAIAKwKtEtELOfO0hhOiZnJz8C51O92BeXp5+4cKFaT/96U+jYug91KxWK9555x04nU65ozAWcYYPH44pU6ZE/DzM4uJi/OMf/6hdtWpVXSAQ2F5ZWfkqgO95LlT04mKKtZlonCSVl52d/SiAm2fMmBH7yCOPGAYNGiR3NEWx2WxYvnw5N/VkrI1iY2Nx2223YciQIXJH6TCHw4EPP/zQ+8YbbzhsNluJxWJ5xev1riMi/kbQDXAxxTpECKFTqVS3ZWVlLU5MTMx94IEHEu+9997Yvn37yh1NEcrKyvDuu+9etM8MY6xR//79MXXq1Ii8W8/lcuGLL74IrFixwnrw4EF3Q0PDCrvdvoKIKuTOxsKLiynWaUKItLi4uGnJycnzExISes+ZM0c/a9as+KuuukruaLIqKirCBx98AL8/bN0uGIsYKpUKP/3pT5Gfnx9RdwY7nU589tlngbffftt65MiRer/fv9ZisbwD4CB3J+++uJhiISWESNbpdHelpaU9EhcXd+WsWbPiZ8+erR8wYIDc0WRRVFSEDz/8kEeoGGshNTUV06dPR05OjtxR2sThcGD9+vX+FStW2E6ePOn2er0fWq3W9wAc4wKKAVxMsS4khEjSarV3ZGRkPKLT6XLvueeeuDlz5iQMGTIkon4T7aySkhL885//5IIqSKvVIjY2FjqdDjExMdDpdNDpdNBqta02jUYDtVrdvFepVM17lUoFIQRUKhUAtPr31PQ9LRAIgIgQCASaN7/f37z5fD74fD54vd7mvcfjad43NDQ0bx4PrykbCkIIXH/99Zg4caLiV0ywWq1Yt26db+XKlfZTp045PR7P+zab7V9EdFzubEx5uJhiYSGESFCr1bdlZWU9olarB0+fPj127ty5iSNGjOgWhdXRo0fx0UcfRd3CyFqtFomJidDr9UhISIBer4der0d8fDzi4+MRFxfXvMXGxiImJiYi7wINBAJoaGhAfX096urqmje32w2Xy9W8r62thcvlgtPp5ALsHBkZGbjjjjvQq1cvuaNclCRJWLt2rW/FihX28vLy6oaGhn/a7fZV3MqAXQ4XUyzshBDxKpXq1uzs7F8KIYbfddddMXPnzjWMGjWqeaQhGu3btw8bNmyQO0abqVQqGI1GGI1GJCUlwWg0wmAwICkpCYmJiTAYDNDpdN2iGO4Ij8cDp9OJmpoa1NTUoLq6GtXV1XA4HM1bd5hPJ4TA2LFjMXbsWEUW0pWVlVizZo33nXfecZhMJpvb7X63urr6fSI6I3c2Fjm4mGKyEkLEApjUo0ePR/x+f97IkSNV06ZNS5k0aZI6khY1bav9+/fjs88+U8wIlRACRqMRaWlpSE1NRVpaGlJSUpCSkgKDwRDVxa3ciAhOpxM2mw1Wq7V5q6qqgs1mU8y/kc5ITU3FXXfdpajRKJfLhe3bt+Ozzz6r/uqrrzx1dXUmp9P5rtPp/JCXdGEdxcUUUwwhhArA4KSkpNsSEhKmqVSqK2666SbNnXfemTxhwgQYjUa5I4bEoUOHsHbt2rCfNzY2FtnZ2cjKykJmZiYyMjKQnp6u+Lkr3ZHP50NVVRUkSYLZbIbJZEJlZWXE9C5TqVTIz8/H+PHjZf/35fP58MMPP2Djxo1169atqzWbzS4hxJbKysqPAWwnIpesAVlU4GKKKZYQIgbA6PT09Ls0Gs2UpKSk5Ntvvz32tttuS8zPz0dMTIzcETts48aN2L17d5cdX61WIycnBz179kSPHj3Qo0cPGI1GviQXwYgINTU1qKiowNmzZ1FeXo6zZ8/C61XW0m5XXHEFpkyZgszMTFnOT0QoKirCpk2bvGvXrnUUFhb6NBrNnsrKytV+v/9rIjLLEoxFNS6mWMQQQhgBjM/JyZkZCARuvOKKK2Luvvtuw+TJk2MGDx4cUZekAoEA1qxZg8LCwpAcT61Wo1evXujbty/69OmDHj16QKPRhOTYTLn8fj8qKytRUlKC06dPo6SkRLbiSqfTYcqUKRg2bFjYi3ZJkrBlyxZat26dbceOHQEAJ+12+2q32/0lgCJuX8C6GhdTLGIJIXrrdLqbMzIy7vX5fENHjhwppk2blnLzzTerlTRH42L8fj/WrFmD48c7dqe10WjE1VdfjauuugpXXHFFxK9nxjrP5/OhrKwMRUVFOHHiBKxWa1jO27dvX0ydOhXJyclhOZ/b7cb27duxfv36pnlPUn19/ac2m209gB+IiPuQsLDiYopFheB8qyEGg+G2xMTEaWq1uveECRM0d911V/L48eMVO9/K6XTi9ddfb/NogtFoxODBgzFo0CBkZWXxZTt2SRaLBceOHcORI0dgsVhCfvy4uDhMmjQJXd3ixO/3t5r3ZDKZXEKIr4PznrbxvCcmNy6mWFQKzrfKz8jIuEutVk9JSkoyTpw4UTd27NikUaNGoU+fPoopRHbu3Imvvvrqoq9rNBoMHjwYI0aMQO/evRWTm0UOIoLJZMKBAwdw8OBBNDQ0dOp4QgiMGjUK48ePR3x8fIhS/p/q6mr88MMP2LVrV/2WLVtqjh07FtBqtXsrKys/DM57MoX8pIx1AhdTrFsQQiQDuC41NXVcXFzcBL/f36d3796qcePGxd9www0Jo0aNQlZWlizZiAirVq1CcXFxq+cTEhIwevRojBw5skt+YLHuyePx4ODBg9i5cydsNlu7P5+QkIDp06cjVIua19fXo6CgALt27fJu3brVUVBQAK/XaxdC7KqsrNxMRHvA856YwnExxbotIUQPAKMyMzNv0mg0Y4koOzc3FxMmTDDk5+fH5uXlhe3y4NatW/Hdd98BAOLj43HjjTciLy9P9tvKWfQKBAI4cuQIvvnmGzgcjjZ9ZuDAgbj99ts7XNz7fD4cPXoUe/bs8W/dutWxd+/eQG1trUuj0ew3m82bvF7vbgBHiUhZtygydhlcTLEOEUKMB/AqAC2AKiIaF3z+VgCvAVADWE5ELwSfHwRgOYCTAB4kooAcuS9FNF4/6yeEGJWdnT2JiPI1Gk3q0KFDMWHChKTRo0frRowY0SWjRAUFBVi3bh1GjRqFiRMnIjY2NuTnYOxCfD4fduzYgW3btl10/UiDwYDJkyejPQuWBwIBFBcXY8+ePbRt2zbH999/77XZbB6NRnPEbrdvcblcOwAUEFFdiP4qjMmGiynWbsEWBTsA3EpEpUKIDCKShBBqNBZLkwCUA9gLYBYRHRNCrADwFID70Dhk/6Vc+dtDCKEBMECj0VyXlZV1s8/ny4uPj0/My8tTjR8/3nj99derhwwZ0ukRpIMHD8JgMITs0glj7WWxWPDRRx/BbG7dhunKK6/Efffdd8m7RYkI5eXl2Lt3L7Zv317z3Xff1VdWVvo1Gk1RbW3tNw6HYzsa77Kr7uK/BmOy4GKKtZsQ4tcAcojomXOezwfweyK6Jfjn/wQAIvqTEOIdAE8AuBfAKSLaGN7UoRNcAmdofHx8fmpq6iSv1zvEaDTG5efna8aMGWO85pprRG5uLtLS0to0Wbyqqgo6nQ4Gg6HrwzN2CR6PBx9//DFOnDgBABg6dCimTp3aqmdZfX19c+uFgoIC17fffus6deoUqVSqUo/H853FYvkWwF4ikuT5WzAWflxMsXYTQjRd3hsEIBHAa0T0nhBiBhpHq+YF33c/gOuJ6DdCiBEA3gRQBODnRBRVK7wKIRIBjIyPj782JSUlj4gGBQKBdIPBoBk4cKAYMWJEwqBBg2Jzc3PRv3//Vt3b/X6/IheAZd1TU/8zvV6PuLg4HD9+3F9QUFBz+PBhb0VFBfx+v1Oj0RS73e4DVqt1H4B9AEp5gjjrzriYYu0mhHgDQB6AiQDiAOwEcBuAYQBuOaeYuo6IFsqVVW5CCD2AqwDkpqWljYyLixvh9Xr7a7VafU5ODoYNG6YdPnx4Um5urio3Nxc5OTnc+oCFjdPpxMmTJ3HixAkcOnSo9uDBg+6ioiLhcrk8Go2mIhAIHLZYLHu9Xm8hgBMAzFw0MXY+Xm+CtYkQYgGAh4N/XA3gy2CjPJcQYhsaC6lyAC1bj/cEUBHWoAoT/G9UENw+bHpeCCHKysqyd+/enavVagekp6ePUqlUg30+X05CQoL26quvxvDhw+OHDBmiz83NxdVXXw29Xi/b34NFLr/fj5KSEpw4cQLHjh3zBkeZ/FVVVRQIBBwajeaE0+nc73A4DqKxYPqRiDxy52YskvDIFGs3IcQAAG8AuAWADsAeNE4sP47GCegTAZxF4wT02UR0VKaoESnYcLQfgFyj0TgsMTHxWp/Pl6tSqQwZGRliyJAhmkGDBiX06tVLl5OTg5ycHGRnZyMhIUHu6EwGPp8PZrMZFRUVTVugqKio9uDBg/WnTp0S9fX19VqtttTr9R40m817iegEgBNE1P4mU4yxC+JiinWIEOI/ADwIIIDGFgivBp+fgsaWCWoAK4noj/KljD5CiDQAuQCuNBqNffV6fT+VStXb5/NlA0jQarVao9GIXr16oXfv3rq+ffvG9+rVS5ednc1FV4S5UJF05swZV0lJSUNpaWnAZDIJj8fj8/l8DVqtVhJClNfX15+2Wq1FgUCgDI2jTGeibX4iY0rExRRjUSTYKysJQDaAHADZSUlJfRMSEvpz0aUMFyuSzpw501BWVtZcJPn9/gaNRnNukVSBxkvnFQAsvKAvY8rAxRRj3VCw6DKgseBqLroSExP7CSGuaFl0xcfHi4SEBDIYDDAajSIpKUltNBrVRqNRm5SUFJOUlKROTExEQkICEhMTce7jlncuRjK/34/a2lo4nc7mfcvHNTU1VF1d7XE4HB6Hw+FzOBz+6urqQHV1NTmdTuFyuS5VJJ0FUInGIkni0STGIgsXU4yxiwoWXTFobIGRENy3ehwbG2vU6/UpOp0uVa1WG1UqVRIRGQKBQKLf708AoFOpVGq1Wq1SqVRqnU6HxMREMhgMMBgMwmg0qoxGo8ZoNGoTExO1Go1GpVarhUqlEmq1WqjVaqhUqvO2ls8LIRAIBOD3+xEIBM7bmp73+/0UCASoae92u33V1dVeu93eXPzU1NTA6XQKt9st/H5/IBAI+IObT61Wu9VqtVMI4SSimkAgYPd6vXa32211uVw2InICcAKoDe5b/tmlxM7/jLHO42KKMRZWQggdLlCUBR/rAahabGohRFMlptZoNBqVSqUWQmhUKpU6+FgNQEVEfiLyBwIBX3Dz+xuroaa9D41z/FpuHrQueloWQm5uA8AYawsuphhjjDHGOkEldwDGGGOMsUjGxRRjjDHGWCdwMcUYY4wx1gm8nAxjnRBsXjon+EcNgAEA0onIJoQ4g8aJzH4APiLKC34mB8C/gq/NIaLasAdnjDEWMjwBnbEQEULcAeBxIrop+OczAPKIqOqc970A4J8ArgTQg4jeDHdWxhhjocOX+RgLnVkA3m/D+9T4v1vzRZcmYowx1uV4ZIqxEBBCxAMoB9C/aQFZIcRpAHYABOAtIvp78Pkr0DgyVY3GhaCd8qRmjDEWCjxnirHQuAPA902FVNBPiKhCCJEBYLMQ4jgRbSOiEgBj5YnJGGMs1PgyH2PtJIRYIIQoCG45wafvwzmX+IioIriXAHwC4LrwJmWMMRYOfJmPsU4SQiQBOA2gFxG5gs/p0bjEiTP4eDOAPxDRlzJGZYwx1gX4Mh9jnXc3gE1NhVRQJoBPGtcJhgbA/3IhxRhj0YlHphhjjDHGOoHnTDHGGGOMdQIXU4wxxhhjncDFFGOMMcZYJ3AxxRhjjDHWCVxMMcYYY4x1AhdTjDHGGGOdwMUUY4wxxlgn/H+dJryTUwDHMAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 720x432 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(10, 6))\n",
"plt.subplot(projection='hammer')\n",
"\n",
"plt.fill(combined.icrs.ra.wrap_at(180*u.deg).rad, combined.icrs.dec.rad, c='k', alpha=.5, lw=0)"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python (astro36)",
"language": "python",
"name": "astro36"
},
"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.7"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment