Skip to content

Instantly share code, notes, and snippets.

@smsharma
Created November 15, 2019 16:39
Show Gist options
  • Save smsharma/1abe6d8b3eadd265394df346b3b55df9 to your computer and use it in GitHub Desktop.
Save smsharma/1abe6d8b3eadd265394df346b3b55df9 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": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import healpy as hp\n",
"from tqdm import *\n",
"from astropy.io import fits\n",
"import urllib.request\n",
"from NPTFit import create_mask as cm \n",
"\n",
"%matplotlib inline\n"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"('../data/gll_psc_v18.fit', <http.client.HTTPMessage at 0x181d13aef0>)"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Download 4FGL catalog\n",
"url = 'https://fermi.gsfc.nasa.gov/ssc/data/access/lat/8yr_catalog/gll_psc_v18.fit' \n",
"urllib.request.urlretrieve(url, '../data/gll_psc_v18.fit')"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [],
"source": [
"# Open catalog and grab quantities of interest\n",
"\n",
"gll_psc_v18_fits = fits.open('../data/gll_psc_v18.fit')\n",
"\n",
"b_4fgl_ary = gll_psc_v18_fits[1].data['GLAT']\n",
"l_4fgl_ary = gll_psc_v18_fits[1].data['GLON']\n",
"TS_4fgl_ary = gll_psc_v18_fits[1].data['Signif_Avg '] ** 2"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"application/vnd.jupyter.widget-view+json": {
"model_id": "bf1235e7d08943bf8d12e40acf54f98b",
"version_major": 2,
"version_minor": 0
},
"text/plain": [
"HBox(children=(IntProgress(value=0, max=5099), HTML(value='')))"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"nside = 256\n",
"mask = np.zeros(hp.nside2npix(nside))\n",
"\n",
"for l, b, TS in tqdm_notebook(zip(l_4fgl_ary, b_4fgl_ary, TS_4fgl_ary), total=len(l_4fgl_ary)):\n",
" if 9 < TS <= 49:\n",
" mask_rad = 0.3\n",
" elif TS > 49: \n",
" mask_rad = 1.0\n",
" mask += ~cm.make_mask_total(nside=nside, mask_ring=True, inner=0, outer=mask_rad, ring_b=b, ring_l=l)\n"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAn8AAAFyCAYAAACEOLy4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJztvWvMdtt613WNZ60GgiBii1ltPAVZEaIhDXQn+KFq1d1NsCZNdFOCBj+gwAc0EqVRAyQkJPoBNdFggocPBrUpmzRUDulyNYVaJaIWF9VwWkZAcGVFK4ptsO2m7/TD88z3He94x+E6jjHmPf+/ZO/1vPc95xjXHMf/uK5rPk86joMAAAAAAMA9eFptAAAAAAAAmAfEHwAAAADAjYD4AwAAAAC4ERB/AAAAAAA3AuIPAAAAAOBGQPwBAB6SlNL7KaWU/ftrgur5hSml35pS+qaI8gEAwBuIPwDAVqSU/rGU0kcppT+XUnoqvvubU0r/e0rp96WUvtAp41uJ6C8Q0d+ZUvr5KaU/TET/utGub04p/ZHKV/8rESUi+vss5Wf1/J6U0j/tURYAANSA+AMAbMVxHN9PRL+XiP5HIvqW4utvJKJPiejfP47jv++U8V8S0Q+8/Px/EtHvcjDtR4jo36jU9YqI/jeH8k9+DxH9UcfyAADgLSD+AAA7chDRf05Ev6b4/OnlO85vp3/V+LlJSukpDxW/ZdBx/LXjOH6IU46G08t5HMefPI7jM2s5AADQAosEAGBXvo+IviWl9DOIiFJKHxDR5+VFKaXfkFL6p1JK/0xK6cu9AlNK76WU/r2U0l9IKX1TSum7Uko/mFL6upTS7yCi7yeir00p/ZKU0r+aUvqOlNK/9nLvv5hS+hNZWb86pfTrUkrfRkT/QPb5z0sp/baU0q9KKf3OlNLPqdjxHSmlH0sp/ZqXf/9uIvr9KaW/K6X0lZTSb2mV9fKcfz2l9I+nlP7llzD4N6aUvjWl9Ofp2TsKAABNIP4AAFtyHMdP0bMY+5UvH/3i4zj+bH5NSunbiejnHMfx+4/j+E+J6B9OKTXFz3EcP01Ev5mIfpyIfpiIficR/dXjOH6UnsXmryWi/4eI/kMi+l3HcXw3Ef28lNI3Hcfx7xLRV1/q/QIRffE4jv/4OI4/RET/S1bNv01Ef+Q4jt9HRH+ciH5DxY7vpmfP5o++fPTHiOg7juP4S/Qc9v1ZrbJenvP7iejT4zj+LSL6a0T0PxHRDxLRv3Icx59sPT8AABBB/AEA9iYP/dbCsd9BRH8q+/cnRNT1/h3H8TdervsCEX1IRCml9LOI6OuP4/grRPT3EtHfQkT/aErpS0T0l4novZfbv/ry33+SiP5EVuyPZT9/GxH97S/3/mx6I/BK/rP82Y7jOMv+anZNq6zvJaJvfxG6P0RE3/zyvz/We3YAACAien+1AQAA0OG/IqL/KKX0y+ht79pJuYY9ERHnV7p8LxF9Oz0LuO8joi/RG3H5ioh+/DiOj17+/dG7t7/OPazxioh+4DiOH2t8f/JDRPS7U0q/iIj+orCsP0hE30NEf46I/h0i+o1E9EPHcfy/gzoBAACePwDAlvyMlNLTy5u0f4CIfuNxHOcbtYneCLXvoWcP3skve7l+xPfRs/fuR4novyCi7ySiP/Py3Z8noqeU0tcSPecappR+SXH/9xJR/tnfRm+E6PdQ9pZySumLNQOO4zjoWcT99uM4/ruGndWyjuP4P+hZfP7cl1D4LybeSzAAAADPHwBgL14Ezm8ioi+klH47Pf/aly+8fPcv0LPQ+fUppR87juO7Ukq/OaX0a4noZxLRf3Mcxx9PKf2DRPTLiehfSin9m0T064joF6aU/sBxHH/qOI4fTyl9PxH9t8dx/HRK6UeO4/jTRM95gSmlX0VE35lS+mF6Dsl+d0rpnyOiX5pS+ieO4/iDKaVf9PL7+P5vIvo76DlM/INE9FuI6LemlH4uEf00EdV+N+DJ76VnD+T57N9ARP88Ef3dKaX/ZFDWd9Nz3iIR0X9NRC0BCQAAb5GeD58AALA/pzfw/HUsBxYwAAAQA/EHAAAAAHAjkPMHAAAAAHAjIP4AAAAAAG4ExB8AAAAAwI2A+AMAAAAAuBFX+VUveCsFAAAAAOBdan/9qMtVxB8A4CZ88an719mW89Fnn7zz2Ze+ofnnhJfw8auvrDYBALAxV/lVL5cwEgDQRirqaiLrZJXY6tl0spsQzDntP22ESATgIRB7/iD+AAAuWDx2HFF1MlNceYm9WUK2VU+rjtHzPX3wqdkmAEA4EH8AAB9mhV97gkUimkqvlhc7eiC1SMVh754aEIsALAHiDwAgY3aOnURMlGhEipdA++izT14L0quJvh5lu3GeUdqHEIUAhALxBwBo4yX0ogWQtPxo4SetYzeBuMp72asXghAANyD+ALgjPVFX8+xw0HrozvK9wrk9RuV4vJk7Q1jm9ewkGonGbzdHC128lALAEIg/AO5AS+xZhYpF8OVh0VH9+TUzPU+juloh0JIdBNqsNmzVmTNTCEMMAvAOEH8APCISzx7R+O1O6Ruq3OulXrsdRNTJriKvBVeIRbaz1qvMLZtTHsQgABB/AFwSbS4eZ4Oc5a15dFaJValgf9R+1YpziENwAyD+ALgC0W/YjsK3VxYIK8KdXLzz9iK8kVyBaPndizv3EcQgeEAg/gDYiZ3/VNmVwpyR4UUvNOF3zrX5PR7PbfEE9/qhJXy1z9q7x7stvvQN3whRCK4MxB8AK7GKvVlvfF5BTF2ZljixvOXMEVeWF3tmjAHN+F41ViEGwYWA+ANgNlzBh9w7QCQTgNG/yuZK3t+VQAiCzYH4AyAKicir/a67Wois9l2rrFGdvXJqdY6u15TxCAJ3tVeUm2vHFYoab2A+fj3CuBJmC1pLf0MUgk2A+APAE8tbuEQywSf9CxKc6zj1cpD+7j6ubRoiXyZY4QmLEFgz+pxroyY/L/KXgFvzMFvjD0IQLATiDwArEsHX+sXGrU1SkgdmFX3e9DbxWd6y2V6nlR5Lj/xP64sRXr+w29MObTmSvtUeds5rIATBZCD+AJDAEXqcMNxVwpoRHqHVz77SHq4oibZxRRuMBN0VUgCi2q0m3CEIQSAQfwCMePX5h69/zr12o5P9yMO3K9ocMk45q59/tSeUW7fHrzvxvr52H8f+Wj5r7V6P3zUZMd8489qrrlbdTx986l42uDUQfwC0yEVfSe93iO0oeiRYN+GeB4ezOffEhXVzjxTjWlFkKa91PfceDZpctxae4rd3v7UtWvZABIKLAvEHABFf6J14/dJYLqMXQrjiYOWv9eiVx/WueNgRQZRY2fkgIRHpOT3vea+MGbl7rXvK+6RjuVWHR18iPAwUQPyB+5Ln73kky3tjFQSSF0e8PWw1JOHk6PyvyNytlrjhvuSTf7+z+Buhed5RGed1Mw9f2pdHVqR8QAgCJhB/4D7s+KfTJN4tzVvAM0QdB2kIMPLXs0SXzxXtVxFxPWaIU64wXvXizqjuqLebOQckiEHQAOIPPD7cN3QjN62aCFu9ge3A7PD5LKLFpdYDpc1F7HnGrS+acO7p2Rh9UODY06tfe51nigZEICiA+AOPSfmGbo0Ir4U01JVfV9vUdwxHX50ZYiGqDk25mrdsNR7j0Vj1COWvEno9JDl+o+s826i19uClEUAQf+BROMWe9xuUHDSJ75INopdHJrXHG6sYsW5q0vskeZS1ujQbeKT48xwT0mtqaAWp5F4vPL3OkR7s6DULYvCWQPyBa8Px8EUg9ZZo3jBc+Waupr6yTmmen6QuS3uWgkmbd3l+3xJg0V7bVQJKIwB7fbHCg7eD+JQQae9ZNkTgrYD4A9ck9/TNOnXn9ZWfta6vXSOpS3oP9z4LnPD16HOt14pzfcvGWd6xmewoAGd7Q7VYw91WT7T1Xs39nLIhAm8BxB+4BudLGxJx5fUGnTSUZfXYjcrw2gCsQobrWePetwM72zaDGc9/xTxWad6k9oAiLUOKZk3EyyIPCcQf2JsvPn1Z7GWQki+Is8NQlk3DQ/RJ6i3v9RCM2nKigPjTe653I2IuS/J7PQ6BZzmjZ5mxbkEEPhQQf2Avyr+0Eel5sNbhGQKU2GS138MjannpwpPIt2qvKHgstHIhVxyGJON/dP0qPA9po3ST2Xz02ScID18biD+wDzNe3vAQC56LsLSsFRvAri+iROc/7SIoamPWW/Ra+8kjX3LmAYjDjuHpKM+sdjzBG3hZIP7AWl59/mHYSxtlmR517Cz8PDYr736YkaA+KnsnISeF8/LM+VlUXZb783KsL1j06twlQjDzcOZ9sLI+LzyBlwLiD8znFHyPyugNQUl+Xe0tWUl95fel0F4tzqwhaO5LJb1re/dxvU9egrvVHqOXA2rM9mxzX2DQvDxheeHCO0dW2x6cseQh3iI866N5VvYRhOD2QPyBeez4t3W90Wxc+b21Tassa3RNz57yusicuVadM8r13LBbbbSTF7jG6BDgaU95L3d8cdrWIho590mRHuBGdXuEcqPmsaTevP7zc4jAbYH4A3Mo39p9RM+fZkPO75UKt/way6YSFXJv1edVrnfZtfJLRgLFq16OcJMIVU4d5zVSz7QFz2foHQx2elmihnfIffbz9exHXuCWQPyBGFq/okXrcfLOj/Nk5HHQ1C3Z5Cx5URwbtc8RnZc2a4OLTvyXeoe8PIze4sIiBmupCaP6PFgRKo28vyzDqxzP1AyIwS2A+AP+1Lx8ltCGNg9LUocFqwdEen+ZZxMVkpOWV7MxZ9RvLS9PZPjOszwLESK81zdWId+yo9dnGlsivNK1Ok68x5RHmVq8PbOezwEBuByIP+DH6EUOzUIi3TykYaTyvt3CQTXOZ5YKr1Zy9vmZ1Ibe/dy+1gr7nYRfL2+NO/684fZPzV5uuZJUg9b1O7Crhzo6dUJz8ObcJxkjyAlcBsQfsBH15q51gWox8iK1BNJINHJtjYazWFs28FHZUma03yhfLvIZVnmCpLlxHLQHNmv4sDy4RMLxRM6wo1bniVV8Sa/TIk1ngBCcCsQf0LHTm7ucZPj8e+nmN1v4jcQEx56eXbXQ3Kjdanh6SFbm03nWMRJ/OxwQiMYe8h7WsDOXK3oRJUjaR7LWrPKmn+Vx1uDevQgJTwHiD/Ap//QakV+4jFsW955eKI4LxysWJfq49kRg8Va1wom9sma16RW8lDtR28g9UgakNtSY3f69Q8rs9tDkBxPJRGCEKCzhiFh4A8OA+AM8asKPKCafhZvLxvWA9UK93MU0OnzH8cSsCDVJ8vVyRuH1Vhm7iqqdEvmJ4l+E2EXsnPVxclxX9ceVPJTSFIWcFXm2EIBhQPyBPpGir1dOLTTJsYMbArbkEs72wElCv7PgbhQS8bf6mXrsJFJ3sGWnfttJmEtSL3axUSICV74Mg3CwKxB/oI41p097GvdcFKWLzA6bao8V9nFelvD0vHg844yNdaX4WT1OV3vaWoyEl+RFp5Yg0qxnNXs4n3u3bzkvrHX1Dqetcj1sgAh0AeIPvItW+Hl52DzFHzcncPWGOsLbPk1Yd6agstYZKVB2GCurBNgOzz6CK7i430XYUn7H9ZTvhvQZvHILIQDNQPyBN3h5+4jkLwiUrBAaUXV7eKJWe/086hy1g1d90W11tQ3am51Cvi24L3bNyEHleBx3z4VtPYNEUI9SeGrlj4AIVAPxd3ckgi8yvHcyS/RZ6pntddkppymKleH+GawMlUaEwTlCXprfGdVGmoMeR+xYQqWSMnbI8T3tKKm1jcTWkSdQ8uwQgiIg/u4MR/hFC4+rbdSz7eWI5FF42zOfB/DYpf1WH1S0OXeta/JrLdEFywtf0nI8bNkFbv+W382yCW8GsxGLv6cIK8B8esLvo88+YYsODmd5rYXv/N8utOzhPm+0LXn45LSpZlvvO27dQE7e3ivGzMiO3vz2qCvi+hLJ2NSuL15rU77+ea11nPaL6OeRN1ezP3jS+u0UwA48fxem96fYaqc268mXcxLcLbelxcwcwVb9klN2ZLtyvS6rwlWWej08uzuN6dG8zv/t9bxlfdrydlgLPNNTol7a8ijrivW3gCeQhdjz936EFSCe3omotSj1BOBosksShK8EVxTn1ISbNOzUup5TjvfCLOn7FcwI13P6cKV4zL1M3DQGTkivVl5NUFqIOCR6CcrypQWpB3I0NyxpGrMF2E6CLyfvo1effwgB6AQ8fxej5+07ORcOa34YES+BuobHCxjWcjzRblY75QRpxsMsr5ekz6Vtat1IV7wUYKmf47nz9IBxbbEcPr3tOD/rtcfqtWeXF0O4RK8V+b5GBC9gAXL+HpVXn3/IEn5Eb5/ga/Ry9qTCr8wLseTAeOS0WO5vtYuX8FudexjhwZGQt681f1GCRTCNxoKk/lm5sHluWzk3R0T2h8ZDHmUH0Xpxl9MSn721ejbSA3+EzXnfnXsi0AHP3wWw/r6+nJ644wq/nldRu7BKPAJRi3fLO9CzhVumJfTl6QldmbeXw813rIUnW9da7auN6x3y/SI84bNyXlvtJ1lrPGnNR45QWd33OTsJ1xpR0ZtWufi1MPhVLw+FRfRFLxqlC95a/mixkNQzCvFYyvaCU6d321o8lla88t6sZWjqtOSFRdiz2gYtXmKuF671yK+bJYx7NtS+bxFh0w4h59FYr7XJjUPBEH+Pwkj4jTxvJ5H5RZ6eGK5HclRH7dodFjIinXezdo1EAGjGQ7T4izogRLJLHqq3+NtlbkjxFoC9Omb0vTTdpuSKfehBrY9uKgAh/h6BUW5fbcHuCbHWwmINqXmFSbUhQY49VjSnz9b1UvtqdVuEMKderVjUPo+GFUJsB/EX2TfRZUQzIzwa3Q7a9UFyz6547Bs5NxSAEH9XJhd9FtEj2fQjwralHdq6d/H0tOrnnMS9xaHG67dSrESFpE488iRni1cL3OcaHaii+6Ws01LejBQB7zaRrH0eHsxZAt1zTR4dbqXPX9730Wef3EkEQvxdkdLT5+HtkoqX6IW6l4fHWbhmbbxWz1p5j1aISz2/2rLK66/gPRi1qWc6wpVpjWWLUJDkplna3DtH1Hu+SOquldk7LHvaEUGEwyBCAN4sFIxf9XJFtMLv/P78H6d86b0lH33W/rNuvXvK6/JNRBqKjMLa9py21OTzSMuQtlfres92P8cNt8zetdZQt7adrkb+3K2fpf0yS5CM+phjb219kTyr5R6ubXkdve8teM9vz+hQ/l1+ANXUUfNua8u6A/D8LaL8/USep2RNmREhWE5YlFuG1hMm9ZLldWk8JK17uO1reY5aXb3wc8vbELERRXhhRukO0g3W6v25iucmxxpqlN6juTcf+9EeQa5nUzIXudENTlkn1pBo+bm0TCkj21t9PLqWU1/Jg3oB4fm7Al7Cr+eF40yOlgdktGnORutNy7+TClHtYq3xqnA2FEl/ntdLn/0KYuVktEHMSA8Y9e3seSPxqp5Y2knrVcnrl3jwyp8l9XH7yeNgm1N6H73GhLackfCdwcg7z/HUjvq014/4xdDPwPM3kZHoi85R0Xq4pEJIGhaWkoslzilxZIPFgzGyMUcb6ql5uHon+Lwsbeh0ByHY6hftPOHeI7GtLNNLWHkyqz+lzz6rrbw9XLP7uJwHnmv/CqRe3xOPdaC87oF+OTRe+NgVifCrfe9BK0w2W/R4IjnJzlrwuG2gCW/Wrpee2Gu5Mbsxyg+qfS4pxzMMt2P71YgWLdL5Nqv9ItMYZovqs06rAJSGrVfQqlsSIuaU/yACEOJvR7749GV23kX5vSfek6b3vbTsiHu0nghvz5Cnp6FXZs/L90hCxdqvO3rnopEIIM1YkYZvZ3P1PvcSsFdbB2bxAAIQ4m8n8r/SsYvbXeptaoUPuSFKri2c6897vMVM5LNpN1CJF1PbV5y6a0J/9UZqrX+m/V51eZQjSYGQ1CXxPEN0vM3quVTS8jKu9BTPFKwXFoF44WMXNMIvKuFW6l0cTTYPcSSx5yy7LN9jE4wMCY2uKZ+JSy7yes+SXzuypfXvGYv+rHuvRG9saNvgnDPefTpbtORto51HlvkXhXd/a6gd9KL6d3TALZnRX6M/q/pIQPxNphWuq23AXgM9r3O0+I/qjBCuozbxaoezzLI9enbk93Lr0FzTE2CnbaO+Kz2i2kXbox2iOG2zPh/n/lXPOhL0nPu0WDf6KIHJQfv8Mw/dJZwxOLstZ/YfRwDu4BF9RBD2dWR0aigncm9xqCWma3OdetfXytfkAOZ4eOTysjzK9woVzsgHtOY/Su+VwLX96gu2R+jcMm57Za3KGe4RnfMXFVpujddd2tXC6mfQrmOjULNGEGvX/wuFgRH23ZnWol07aXFCet429U583Pq9hZ9HWZ45YpzPuYyEn2fZmnKl3ufVXsHd8NpodwgHWpkVOozkSu3dYgdvduv7XNSNIhzccjW23AWIPyekuQLlIO95TDS5bjXBYw2LaEPC3PLzcrQvL9RstYSiRuF4bc4exxaJvaNrNc9e2tl71kdbULnpANxyNDmwq8KnJZz8uNJOT7tb9WrblsMO7S6hdXhe+RyaA6jlkCoZn9x6HjkHEGFfI7XBYcnTsHqpuDZwQ8LlNREhWI+8Fq1IHF1rDWdZ+pPbLl5hdy41u1pj5WqbqAUPD3NrPqwMp3utSWVZOV5zSDrmIlIwQB1O32jbe8Y6u3kIWBz2fT/CijshPamMFnGvBYZ70vHwakioTVLvRVWywEjKk4amLd6w3INZ877tdKJvfReVc8gVxLPbyDsfdFS2dx7qjPbyiEJwxx/XntZc05QH2mjz9LhREksI2DvH8ArA86ek9/cBOblctcFaE2wrTtqj8iTJuz28Q0Mau7jCynIi1S56eX1R3r0IT275nUe4/UoCR4vFG22NEniUM6PeKFuv7uGLPCR6lCX13pbMiGZw69rUA4gXPmbQygPwyuWSJtpLKPPpyv/Nguu9GdlltVuSz6jJobEuWpF5Ta26Zt2XM+qHHXIKI+aIZu5F5nKNclwl946w5jW22k7TTzXhtAsje/K5kws1iaCOysfTlLOL+G7Z+Sh5gAj7Csn/VJtFCIy8T9Fo67GEOSX5fd4LTVR4NoJR+0R5/XrXaG0qx7bl2ayeOy/PhRde3k1NblRUXpS0zBNu2aPwXv7d7p5eLqN1s2wTi7c9qr2k4X5tGFdL7blzG8rvv/j05V09gGwQ9hWgVfxWweMdoo0MG3HLHOVY9Mrg2N+7xhpuOu/bIf+Ow8jO1veaUI1lrNbqWhECqtGzrzeuRnN/5hjSphOsDrWebVhry1abl6xKMYkmHz+WfhqNg6usdd6MBOjTB59OsmQIwr5RvPr8Q3aI8Lzm/N+qCcPNPZzpEZGEZLgh9HLRn3FitIaaOeFsb3oifadFXSP+o9F4jjnjcXXb79TvNfI2bHlnyjaMfiaNpzWas30s3uNeiknEON0t4lJj9My93P/dgfhjkHewdcDmm34tV6PEKzdOW3ZeB7c8ax6P1LZa/a3vNP3nFZ6T1jcSh5w+53gDRmVxvX5elPXVxtNIcEVw2lHrH01e1kp6c6T2ecTmb/VOtT47y5bmvHHmm5TauPVMaZH0I4eZh4HdDx4nI2F8VQGIsO+AsmM14dte/kIZPuSEK2eFNyRwRUbre4nA8ghv7NhOEW0g9aLtvCBzx8/MZxgJ5t7c5ojqVghcM8c433Ps4H7nxYy+nTVPokPorb3Cs56d1gjpc1rsH/Xd4hxAhH09KXP8RrlILSQelZaXoyxjtfBreT16NnBP45HPYfVKcr1tNQ9vaQfns1b5XCRem50W9RravpvtafMKQUakY3i24aw0ix3GZYQHz7PdWv3q7c3sOTKi6dlbrrfltTMO/Vd7Cxjir8GZ41fbPCW5Wr1QSqusUY6Q5HPviTojbJVP3h1CZBY44b5SPI826Jq4sISlIkJ6EWiecYYwGbWdd2iOyxXnPseGK7I64tDDIxVmJqPoSW2v1kSIent56/mvFAJG2LfC2YG9ztaenmuij5tTFRHu1dhhrZNTdq8OqeveY+HVhFlzRidmj/ynVjmaUIfXIqrBMzRZK2vnjfjEGh7U3l9rR85Gv3Polzv+W3XUPpeMoZUhY24agKZM6X2j8ix9VCuHY+do3ejd27NxwVvACPta4Qi/GpxQILeskpZ3weq16S1qHJu0tFzynnWMyo7GK/TngaZ+bYqD5R4Lmo3rql6kaLgevrPNPbwq3Lo0cO+ViI9Weo6lfgmSdZobRViNdU8YPWtr3+E4U7gRnPO7K3gA4fnLqMXsOSeN1vecMrRYTvYlNa+U1hMl8R7mgtN7EYr0Urbq4Tx7KaqsHr/I55KcoPP7RteOnjvSMzer/ax4eve1dc4SLbv2QUlvXEasYT0bTrwPyBrvWIQNEkHWsoEjIDXRGI6H+Pxu4ksg8PxpkQq/kqssYCWticY9WXrkAUULv13gLlC9Ns29vSvGnKTO3vNe2SPLHfMezzizjzXhXg67zkcprTm3akzP9ChGi3aOhzJiPGpSeCRl7fwSCMRfh94mG+HN4yKtu7U4edstCa9IwnQeAnNV+T1GoaOVm2YpkjzEpleoylpHdPkeSMKKPVaN7bz+/L8z18odntsbazrGqMycsvwZXtteKoFUqHmG+8vrW/dcyQl0+7BvT5l7iaVRuCAvn7NIjtzZtTIkLnJL3b17LCHOiPCoR/+eIRGPzW3lgcKKJezvXUckrTETGfbzmD/S+72eZ7QpR7Vb1FzihPxmh39r6RmR0RRNCDh67nLq7+2PtbJ613DsKe8NDgEj7OuFRjzUPuMKv1Z5tXt6pyJOuLCHNKw3KtNyOrUsFC2vScsT0lsIRnkhXt66qwm+E8kzW7yIq0LdLTxCyLvC8TBxvIq7p8to1uDatT2PVQSjeRTp9bR42SLt4rbFjHVk9bjmcGvxx4nH98K+58bfCx1ywsbndRxquQ8euUulnZY8QO31LbwmUmsRkITOOWVasC5Ms0Nd1vqulg/mFQLXIG2rWTZKBaClzzW5lp5jrNemV9jwteRjfseQp2UNr5UVMcd3y/+7bdjXoyOkIdBeiMPL1cwtY0btdRQWAAAgAElEQVSYYpe8yBLNKVVyqpTi3feWcmplcsbSI298VqTtU16vSbXg1CEpq/YMmnE7Y6x4z4VZ43v2HrCKFevFjJA4t9yg8C/CvhyihZ9msfY4bXDLmLVArPSSlHVawjJ5SN1DVHrdLw1hR7HrJlNDGtbjhDZHSNvHc+xy8fKOWcPAHswcjzt7rne17UrrhQRue+/iAbyd58/S8F6iqSxHenK4sqclysvV85Z4nPzL+jyx9mfr2a3jyUtYcpLFPds04mWFneab5vlqCe81dnrOk137gchf7Gq9qZq6VnC1vUsTTVvkARR7/t73rP1KaDbIE89wi+Z05j0IZxGxiHO8Ipz2aG2Op+eylQPp8RyeHsNRrmatrlEOT2SO1vm9R1ty+5lzzW5zp0RqY36tZw6exaayfk3e1u79JMVjb3mkNrnCXLwqtwr7nl4/y4sLXt6qqPtWJPxrkrAj0bYZ5z4vQWQJJ+4azqmx04sR+edcUUq012Zaewav8LSFO4zlGjuNjROPNo0YTxoht7p9Obm3+bUSe1eHf28j/mrCj9tRtdw16eQovUflSbx2GpfmJ0nwKFv6/N45gD2PVStvqqTniW31UYQnTNqWHtdI8chJjBIpUTlrvXHTqzMXZL0+j5zjXFZvsCet5+utlZx21tbLva+1rnE80B609g8PvMvprQncumYfGLzz1nP7VwrA24i/nFFnchZ1S715qNI7Nymva2QLl102hxNuu0kT5WuhMW1/c+7rheK45Zf1tNqG02YrvTCeqRSjMLYH1rkr8UL2bPCC+zy1MWcZN9w5UNrn9ewSr3/Npt5nEu9y7TupTZ7toz30SMsbfedxvRet9B8Nqz31RDcRf7m6Hok+TQhwJrVBEyEoOXZcAa4IrrVfxDNKBWkNy6ZQliPxhHt7qaLal+NJiPaeasS85HtJ7ubuaNcsiWDl4DUeORGO3vfa9ISZeAvAqyA5IHFZ5f17+Ld9JQ3bCztwvQvSeiTlrMxDihbDu+VYtezh2Ol1zXndKKTDEW21a6R2cuqq3Vuzv5cGIaWsZ5exxG1frt2t9pzBSEjXbGuN3dp88vbk5WVy29XTDiseh7KIZ/HcB2vlastrrXFWPNq2t4afc+dL3/CN1rd/8Xv+POhttrM8a6vh5mFwTrmceryu86Z8Lo5n0Dsc1/uM68FaEVrshaCtZXPqWYXU0yTpm9FBYTa9iMPsNINV3uZouDbPmAOtvXF1juFO878GNz1nJg8t/qTu1BkdYjllem+aLTjhMq+2koZtZizeu20QnNw+S5iTKyC5KRMSVi+CkaH9iHJH7dz6zvqcnkI9ilbZEsFc+3fUetArW9NOkXbm1NYez37daf3lRkQ86pgd/n3YsK9XQ5bu2RNrmMojfBZFK3QSWT8nzFp+F21Hq76oMJI21KwJdXlieV7PsNFuYV8OrXBp+ZmkrFp5O7ZDVD/V1mxN+oXH2lMTCZJxyrWhV48UySHSuw899zrP9JsoW2rPqwz/IuxLxBN+5amrdwKzhM1aZXKZfQqamT8yqnOmDS0s7T/aTKK4ine05rnUlOsdeppJz/Mpfa7ewcDSvj3bRteMhIT1EM0tV+uVjqBmS6utuDas8qB718l5jijPd+T60fNKl9/N8gA+pOcv4k+4WU4I2pO3pk6PezSnL4t3oXe61YjRli3c59LWySmbW6flgODppbbCGRdWr0Vr/ErLiYDjteUKFc9QYITHVXpfSURfeUcOtFGbFtI1WnK9x7OedWr3hJOat1NanqS+Gqu9f5L6FN4/eP40cE7Ep0LXDA7tgDrvm+11y59z1C5n21lszNt29uRr2eNZHocrtR+H3pzy9kJ5M9Ou3XK7pGVz7I/yPnLI7fMUvLOR2u7xrKujXkRz2jp6n40UvBYeTvy9+vxD0YDhhH6laF3Mnh4QDbW6ON6u3Rjlxe3KTPE7qz04z2N95lL4WoXwapHAESzSUHDZPty0Fw80EQjv+r1zyLhzqDYmV49LIn2qhLaPvNI5eiHzXvvUxnyrHG80/TYj9PtQYd9Xn39IRHzB0jrVeC0UXqFarltaExpolSsN10WLFuvpSXt/r00tISBLeRHtLgmJW8eFNPy+y0lZg7SvRikY0vXq7C9OOZJyNbTmoCU0bU2XmLWO9MrhPienn2tle6wXUWuOxLnQCj9znnc0b6K8c71nHK1vH332CT198Cm3KoR9e4wGsPcJcba3aaewHwfJSYsbvuFuJJJ6ewJmZUjEC8kzeHjypOPU0sa7hJet4d1Iz4ll3GnnsFf9VjETEeYelZtv+j0BIikzeu2ISCuxRGh6KSXaMst7PNaNkTfyrEt6rwcPJ/5Wix/LgPEMV0aEgL3R5LH0Tkk7bPK74t02nM0tqvz8Gm6/rxobMzbMaFptLNmAvfE8fEUIQE+vpeUQUF4/cx5I5ueJ5oBQficR1JzyI1npsHl/Sa0BPMfI296dUSN7d4CkPK7ok3q8OGGN1WLZWr9mMVu1Gc9o61YogbOxtD738k6MyMdrK8QTiXd9rVArdx2KzD9qjQ9rmkMuyspn6a3BI898ay2XhLHL+rwP6lF1a66v0eqbGWjq5ISxy2ul143qW8Vpw6vPP5SEfkU8jPgrWaX4IyaV5KQpyc1YLfw8kCyMq0TfrotthLiybKqtVAwOree1bvIeaNq552njPhN3MytFt9Y+ji2c8izeEMlc230O53iM41oZkQes0QGOu1eNDrKSvpE4REYH46vvnw8R9h29GbNTJ1nyD1ZvYi1WbrCrF2UOnvbNDKnl/5Z6Hjgb+KhMzxD/7mMkR5IDNIoGcDZYL6+wlKiwcaTH1LM8z7QXK1HritabJs3946wl0vazrl8W8rrPF1m9ufzbvrP/Hp4Wbn6H132tMnoLfVTZGqTlrj6NefSTpI7Rhq1ts96CVgvjWdCOHWk4R1qup0fI+owWD1iJdzjbUmat3azlRq1FvbqudKho4T3Wel4+Tp9L1tLRtRH9xB1n3o6JwS9+xtu+V0A7IB5hoeHSEyQ1z0/kqXg1PU/XrHYYhV6iaXm/PZ9XWhbX6yAtLwLJYUD7PYfICIY0dK1l1wjMo+DpnYtYDyPSClbwsDl/I7xOmWfIkxs60bj6PWmdhCyCNCL02uqf/POdQr6zQt+1ttCOZW7uS453m2vKWu3d9WzjsuyoPMxafdw1oDXnPMKzeZ2SOXTeW7aXtN163m9rSDsyROvpcbLa2cvt4/ZnOaYkY2zWHhlZ5wou7fnThny9N2lJzkyv7ty74WljtFdsVdlemxB4Q54bExk69C4rmiuPMWn+VEkvjGdpl1YOIzfXKrpPJO2hucablXnhHntMecg4y12Fpi0j+8A7xe3S4s8Da26J5ZrW9b2JpB1cM/NhPBgtJjsKP0vOzIrnGB1ERtfsQi0tYCbaNSR6TuZCvvYd57PyOw8vUQ/pgXkEd25J6+CKVInw5l6/8/otyRku79v5uWpI0nF2XEcv+8KHVQW3BpvnguwdGrbY1jtJXU0YnlzV7hLLKdeaVsAtU9PWM/onso6ybE9vxC5jV7vpRtjPFT4Sj3QvrFsrjzvWe2HOmo0S0avJOy1Z5Xks65eIHu85poXbF5znbY0FS1Sl8eIHXvjgEiGCap3cOzFLB7plQtROqpKJuTKk0OLKgi8C7YGgV87qNo5Ig7AQFSlYyci+vP1X90UtT9Cbnsf0tKH3uXRd5ZY/QhN65Xj/NbT2tpF90SlKHLg2jLznu+cM3lL8RS1emhObd5jDo/6rEenxeVRGz7jDIryLd8xiw65jSSPieodbDySibnSwnkXvACUREJ6C1kMASveHK4ZtR2ifp+x7zs8ruGTYtxXyHYVWOS59ziD2GOjc3Ijy2siwHsfFPWPArlpIIkKovTrKUNSs8FtEakNZlvccmTkeLCGZUZkeEQBPRnWP1qlZtufjaRTKze2Ltumshxt+lY5pTvt6pY20bNWGoUc2c/fq2s8cW6PXD+4ax9lbpfZVQr8I+1o9abPETV7fyhOA9PQc7cl4JI+kBOspM3IM1fqk5zlqXR9lmyRsyWG10JxdtyXtZNb6taN33/rsknE5EuFagTYbzjOXYrr8efT9LDTzppbGsmr/fxjxN9qIetdI8XDL5xN2xgml56GxlufFboJvN3u8kAj+UshZ+p5zX+1gpPFK1b4f/TsK7ho0c7xZ88p6Hq0VeIXko59jN5E2mmPSdr3qYWAFq+fMw4R9W6fYHRR2j9K+qJBcqw5uPTPCKjv01QwbrpIfw1mcRqGbWjkeoZKy3F7IuXYdp5+1Y6EXguKKJs/x0VtjvOrThh41c2HUlpJxU14vtafWtt7h8NlpAdHrkyQ0nF/DSVGYva56jl8Jtwz79v7Ica8xr7DZPgpc71BrEb5SGJxo77DiKqyeBM51nI1DU7cUrmcv38R6Y2b2+IgSmtxr8/bQ3C+ttwf3YNLzGp6CwKtdJWXVbBu1reZwpkVyoCqfe9QOXnZLx5KkbYl87OzpIC6X8/y9+vzDh908PbxOLa/D1TyiV2FWO0bW80hjQbLYcjyElvpbXrYWXt7F0ecR/d3y0vTq0IpfrSe3Vk5EdEVigzdebWO5z4JXnVKvoaWOk5l77EeffUJPH3yaf/T4nj8OXqdAzWk0Am79s0JIq9sjCq9nWzluIut9tH6P3NC4HuQo8aEpx2PszBgjtTbTtKNH/1/9sDSb2WtjWZfH+ObYX/O+WvFus0uJv5Gr03NQrdjoWsm3tRPFKmHxqOz8bKtDylcX/L0QnTej5PkyjBVlU7lJ5XW31pkorIfSljd3tE7OYnZ7nozmZS8UOVN8rSYqHGyZ67XypDZYQ7/vm+6ejNZlLb3XUu9Z94z8oha1TWWXhbIkwh0fTW1DjURax4q23KUfa+HOPPTT8sh5zdmWUGmF47g5ZpocyNLjkdtRWw/ynLVWmVwbenVx7dfYEBFem12/d8i8lQIwejYvr7SkDaWpEmVdtfEXRS1tRLLPWlJLzns/fqUqgogu5vmbyepNjAPnJFJb0O+MtB16YYOrtmn5TCPPReRBJrIN82cazYPeJhsVCuWGhnM7yoOHVViNruE+e8sT52lfzTZL/3A8wi3hdFVW72vccWgRRrX6Iryzq9vSymXEX+uvekSiWZS591nxHsicRTTiVF372bvs8rOWdyb/L8e+K07+VeLVM9fnqptwrw0k81kiyPL/XungYg2Jca/N+6P2mYddEczKV7Q8LzckKmF1CsxZ7i5rv0UXXSrsO6Jc5E7Xr3fId0WIq1dn6xlHE2W0wHHCSFo8wvI9wT0qTyrmpKEfbfvUQgkaGzhjVPNMvbqkWPrcMgdHbdo6OIyuPTeFlQKhFv460djmtcbV2q+3nrRC5T0816OVcPtJeliQeJrz+3prWXTouLSlVeeMeWfxjI9orfvlnPDkMp4/Lp4CRVKXV3mjUMRItLQ8keXk6NmeX7/Toki0xtNaTsDRyc8SavfOGfJcnDzLm01kWLcsb1TnDHoexvy/PTSHjd53nmPRWtbo2VZ7dqTepVVe9dE6OHMPWe2R83AyzVxfH078RTOjcyImzEjs9a6L8HRqQg6zhYdmMdlNHLXysCL6coYHIGJxzw86lv7rhXN74pNjn9Wr1RJg5XNbxkdEiLBXh3UstJ51tfCTwBmzkv70Tj3wQhPJmYFGoI8iKLPSMx4q7EskD7lJhU3uQfNcsHqhJm6drXpam9LInV8ra2ZYSyImOG0cbTsnlK0db15h5yjPVHRYjqg/vkchk55NVtFS1tcqr2aPdp1q9S+373t2lt+t2nh7z9e6xhrt0awrq+jtGdryuGFGzaF4NHa15WnsqZVXK4tjq3QPlVBbU7zG4UN6/jxPLzWlnp+kJAufdDLVwo3a8ktak4e7Yc4UgXm9V4Gz+VrKiriHU55nuZL24Hquy3t6J+3eGJZ47KK9R5pQXNTc1KYyWLyXXmvs6L4rrS8rsXrIZ3uVd+dsz9YaEvXsl/nzbhF/1m2kpq1qm+MJat0nPX1wbNhxcau1keZ5pZ4K7olKakvt+l7u1Ww8x/TIC6O1o+W543jYRoKO66Vr0bqvNZ48UiY4Hn3r+rCjxyvCBs58zttE4sVdjSTC5CXEJW3JLdeyF3vslZ5jLafm+eS2Z28NePkzb4/9593ymLj19JHD6WztiVdDpBdhJ7TiWFJ+Pk6kY8ajz1tlrDrJRnjwZlAumJx2bXnTW8/O9axz+y6/Lrq/R+W32qzcSDzX1Z4tnmVpIzit6/KfrV7LnfDuV0mUaAbW5/O2tRWWHoWrW/PU275LiT+icZ6KNETAUdoejEJTVq7uFrfmb5SipmwPS/u0cmusXL3PopF6XfIFsneAihI4Fq9DD45gH3lDWp/nB6OyvKiNO8KzEn2oiyjjilxBBEc4ibRlSOdQ1BpS41Jh3xreYZBaWV5ljurghiK55WrCTr1BPjvsri2zFY6UhoZr10rbR9OeMxcADaUXziNcXt5LxG/3k1rKQO3f5T09G7zxCpNJkI5lTcqFNk2DM2ckfdRKG/GKMozGV+szj7o8WL221NIiyvXZIyTd2gM8yvXGWv7Hr75C9OhhXyK+q/lqROYZRJXtwUzvgiTkWbuWG1orJ3PPllbZszwLkpNxuTi3vEb5NRryNjvt06ZmaLxls706Ht4JD5u9NshRWL52rbUNJOuIxkvaupcztzWUc8srhWDnvdG6F1w1pWWV3ZcRf/nA8AhJcBcb7/BHXq+H/R4LZ1QYbMQs0ac9WUk9T606R2KwJaI4dXIEqfbe2rWW+mrXetnHpdb2s7x8HK+fdq0pxbKk3vyakTAa9ZdHf3JTQHprs+RA0xNaUWgOXavW6SgivJBXE34nK+y+TNjX+2/7Rrtyo+vshTw4dWjCmqMyR2VzrvNcEFqb4OjZJaExTlh3tBlbw2+SDWFkt0XwziYPG436misQave1PJ5cGzn1ju716oveM7TGZa19Z8Jd12pjYPSc1rmnZeQA4DomVodxHxnOHukZ3teWqw37PtwveeawQvjtRC8kIy1j5/arbdwtNJubZBPXtNdoYbfaXNscdxd8OSNPjSZUVrZ57d8aGzXzpCVsLVypfzlwPY15O3L6RLpH7CTCdrJlBRHP/4htepmwryeP1olE+nAm9xpueKW24Fo9JVakZXLDP5J8sloeG9eW2rUem/ijCYEcaxqEd5lXZsZzaw45kvlpmUMz+723NtTC4VE5hyN2nwuW9IQaGg+t5VA6o31v6fkjmisALbk8LVplScKyozBo714OI2E4g9biyAnXckXyKKybU058jX1erA7nXYFywx2Nm9FiPhpTrfEU0U9aj3QUXs+pOXSOGK1dZVtKPEWrD8c1WuNCG23Q2m3ZM2r1ctbfvL7cESDZF2rXcVItRp97eiBv6fmbScuLJF2cvBYySSjsFKwSD9cVkExCSVkcgV87vXPs8BKr3DIjGbXBbCwHqRxOe3IW75rHyrud8jFTeqQ1XpHofqytRTuMnROrJ78FZw32RpMu0WOnfuLgfXDgotlvLXZexvP3iDF3Iv5zeT275HTNOe32PuuFL67Sn9YTPKfdON4/rmCoLQ6jcL30+TxypXqbeC/Prjd+OWkJ5XW1hb7lEdcstLWyPDeWWvn553mdHM8kp+z8e26YslV/bS7k/d8L10o889x1qVZPXo4Ebr0tRvZEUfNinj+PxkSkTR7XSiJb1v4efb6ay3j+EKLSn8zz+4n0nqIczb2rT+1cT2f52ZXxtF/qgdCM1ZqnvNYnGqGnRVuGNDzUWuN6YofbviOPwQzPHceO8pqyryXP26pT40UcCU7pmlqru/wsnwvWvc+rjNy+Gcz2ekbsiWX/7rKnXMrzd0WiRKt1ImqFmHXgrhz4K4Rd1ILp0Q9ReWQntXwZK2XuTS2douVN8qiTi9SGcrOv2dArq/zMEi7mepnLelcf0KXP2rNVI44tY6QlAD2weKxnMZpjq7xn3nVGlPfxK929l/w9f7u6UaORhLsi6tTUE9FXXJtmjJPaJmvtpzJEtnqsa8MZmn7ibFQjD5jUO1m7LyI0Wys3f26OV+y8p2d3S4T16irrtLRhxJzo2aid59y+X7nutfpRMt5H13LncdQ6FHVAXuGpXMFt/rxb73Rcu+78+aqeQyK/cK01bFyzaYTn6bVmuySHKQLvsETpLZu1aEnGxXmt1+bWEgwznt0SSo2qryY+OJtza23UeLy586olTjhifuQp7AnGmvdR2m+RYbgI4addD7j9zW1vMI/Itehy4o+LxJNgwRpOkno6or18nBP66vwTrVdiFrO8jdLvtPdY+77XHrVw5+zN5qxvRniMK55Kcd0TWRzPj1d7SrxMvetbIfvy56hDJnetWyl6POfBrDllmUOPIDC995vINrlMzt+JZJH27IjWAsoJn0jCIFwiQ6m9emoLtNQea7imV5+nV9PSxtEhkjJUumIB7YWDPcvlzHmLN5xTRi/UxhVrtfBoL0TbE0DStAAPAVgrQyuay3WE6y3swX0+6YE7v2+FQJHsOy1m2e69nnsw+yBJFB8q9+Cynj9Oo84SSKPrdnGbewgEzzBu/l+Pemubs8XeSK/hKMSiEdNezPKKtOqpiQDLHPLyNOf/LcvMba55tHrP2bK5Z8v5c+3frbYatUG5YbWesWdTC+149jwwj8poCXavjdw6Di1w++i0kTNWVj4Phx32Xa8oWUQ7X87zJyUi94JTZ81DU7Nr9uSpeSlye0564ZuVtnNZdUofUfN4SNA8U2sD13qkW+LMYqNGCHFo3avx/mg93dJ6at/3vKw1O8+NrxVa7XnsIp9zFBFp1dvyeOfX9srjiO+IKEev/B7W+iT3z1zPr+ARuwuX9PzNHjitxaSXs2Ipu4XHKUDqcRuVE3m60oju2s9aj+0jc4VnXT2u8vnWEgIjL6nn/G6tNzXbRpssx7MogeNpL8O7XmuHxzpcu3eVQLnC3NTwqM/VY2eRe0nxZyVfYDVhR23OiAauC94brSeGK1A9wrKtckcneo5NeTmebveZi8HuYZkeXNtH13HzpXrh1F7dLTtG3hSugGzZzG2b894ZITDOmhEduuXQa4uVod3Vc1XjvbekYjwyO4ScR1zy9/xZ0ZzuyvCKxn09qrcVxqihDUl5nWxHoURr+dz6oydYL/y9M5yxs4uXQyNsOKG+Vn1lOacNrTneql/6fSvU2QpJcteL/PvReK2tMTP6njPWpOOx15acNdA7DMsV/5y+tK5vM+c2dwxxwvJXWV898NpbHv73/L084DJme2yiyorwtM2mFYKLrGfGfbPYxT6ON8iy+Um8+rktEgEhtVHiEeAIxPJ7qfdG2r7asTOqR1Nu6Z0vPxvVs8s8qFE+x+62cq+bMZ/vhEUXXUb8rcIrn8cT6WZT+0wa1tKWcVW4Cw03P2s2PRFwtX7TjDXOPZ6pGTUPW+2z3iZWSzcY1duiF9ZcFarThltXzqVaf3HyJSFWeKxKa1qJ9HmixtLDv+1bg3MKLb0B+X3aUMkoVCS180TihRjZzj1Bt7C6/yV1SPvBGuJqCemz7F590WGY2rOV7TQrxGdFIhA0onDUDh5t5FXGKLRZerJq87YVUmzVVbO/Nb5r12rp1TmaPxLvE9fu3lrdGkOcdABuWaP7Ho1ROLz2maZdVqW79NJAuPd4cinPn8bFqQ0n5P9tfa9hpkfihDvYvQda9IlNk7MjCQVq6+LaMovReN6VqHbyEJG1e63t2xJlvfBfZNh2JH7ynzmb2C4enChPMlF7rW3dX65LmjbaMTrFgfOs0jSGKCLXIm0/WVPhLiX+erRCKq0Jt9NiVBKR87fqWbkbiHedHv1rDa+ftkiRbtARHs0Z9Ba+yPE6Gh+jsN7Ia7jLhsVlF0/nbpT9zTlE1w6NHuuI9J4r9Idknqx+ntX1R3CZt33PH1pv/UpDnzmWHJjepm85lUtCvhJvHrdtPAZ763m8w0WtOltYTvAaL2q5eXi2xaoQRiScZ5KG7jh92OuzHtrQU1m3dlyO2kJbF9cWyXprqXfnsa5t012eadZ+tZrRHK+lHnDnWau+iLlH9I7n73Hf9j1puTp3Gnze3jZpDonk+5o4jG7LyDokXpnZjEJ4gAe33SLa13MMaTbPmrd0ZJNXCkPNY8opy/sg6Y1HyJqTv9W7lyv8Z6ZCjGx5dDz28dp679F2Hr/95HKeP6K2969GT7l7nFxaSv6KpyILkd48LVav58hrJCl7ZB83bNwaazu1u4bRCZmz+Zbl9b6X2uZVFre8fOy15hZnzo3GmGTsrIiu1MrseWe8kfRV/n1UBGVky6gcz3Yqx+gV1qBWX1nWmlodUk0gbb+K+Ht8z5+G1ompdorO7+FwhQE/k51OhFL3fO3+ltfEo9+9vCZXHoP53Gzl4kUKuZ3G60lu08jrvyp1YgU1z/mOtlptKueEtb4ZHsOrrUGl51zKjuNOyi3EnwbtgLjqoLiq3SvQHhZaYbOrLZxE+40X7SapEdfWZ+9tPFpvtYaawOCIDmnqwhXH94lVJEjZbV71uGK/atbb3EnEEbuaVKzWOhA5Hi75e/4+fvUVUei3JDpE2eqwiFCIllpYydIuM1z/kaF0Tc5OeU/U80vD1Z602jzqObkhl1rqRt7+UYItnyuWNmjdyx2HZ/29sHCvHg2j1IRyfFwlFDjCYzM/0a5hvdwxTYTjEfqlFa6dFZXp2ZWXkc/V2nV5nRyN4PXXzh7O89drvFUDfufTXM02rQdzR/e59KSnPRWerOxrad3Sfo54NsmJehXc5/YKKXPapCZ+NTa2yh+V0Vo3LOVK2XnM1JD0Vf5sEi9Tb057t5N07SgPatZIWdS6a5kPHvNrFg8n/loLpTRUt7Ngs1CbcLsOTiL+RjNih9y5ndp55/FdztfVJ3lt7qjV7nzsc9ev0YY6ateRoNpt3Ow0p05mjScybzAAACAASURBVLWSnujTRDbycrn5uK17JfdcCa9nKft9xri+ZNiXSB/6lUwCrhu5tuDXPvMOHXkyCqNJkGyAUq/c1cMWV7M7OnzHKbs1JqWHGA8xZsn1kZTNzaWr5eu1yq/ZaA0Ve46LndbDHqP1rfUco/WLs7aN9poeFgHYssF6Xc1z6TUGWv0UPc64ZXPmXmmrV8iX6MLiz0Jkx7fi+2X9nBPTea3Vlp4d1jJm08qbWG1jdP2z8rpqrF4oOfMlKqQ1o1yLZ2SUR1R+FjGfZ60Ru61FRHLRHnV44JQTlbYhvbZ1oLQ6RSyiy0pUCkLkeL902NdTBVvQDlSu+NLiYZcm/8FSv6R8zuc7YQm7RNfRKmuXdj3HpVc4WJPy4e3luiLlBu6Rc1iW3btvt3ZrCb9oO1sRphkHQQ9hxjm0WMqfNU5mtPmJt9655C95zrG89SuBI2pmuKy198z0kEULwBV1rfYwcsff6BrJ9dq2ne2h8R4DkWOq5/nrpYx41R0lYGfO85yVnsCe189jzS6v8ZijM720J9q1SJJ2tUNEKHo+DMQffsnzDFphGukG2YMz4Dkn79y2GeLorKv3fUSdkUTXs8KbsZsHpYfGS3dlZoeQHrENZ9DyZs0UHDWP+Oj6GUjborxeksvaSoXZLUVgNy7v+SPay/s3q6xWsnsrIf4Rc3Jme/0i6pI8w6h9a3b2TsWceiP71Ot0bvWecMqfnf5Qm8faw+Uum6CHp6xVVu86TfmzGNnHOVTv9Gxe0QetF9W6hlrI+9K73xgh33t6/mbl/pWnrJ1yolYvcp45OSOv5uw2z/t85UI7qrsm8lqnYq86e3D7ySvfh7OBrvSyauoubd9pzZHiOXekwoLzuRSvvuA8S55jKRlPK8fKrLaRhIcjifLARumbhxB/q5EOKmsY69zUJR4d7cDfeaOZkWdkOcl623JFLBuw5JmtIc3edV7jbJboHnHVsZQzEl6jcWdtA68DBedAN0qj8dx/LHAjCpxDDFdIcdeR1p4ZwVXm10OEfU8iw78tL0qro1uu356Y4Eye2ilH6/Xz8BZaw1M9uyRlRYS5drFjFr1QhFUEW0LPHniGl6x1asrhhJKkNpXlRJffu15iB9fOsmzPMefdVhah45Xmww2btuqQ7F/n9zPCsC1bIvF6LoHX755h3ytQm7iWxagc2CuERnmKjvSQXD3k5UHks3uHkqS5ifl3kf08M2WDmyowsuNRvBW1OiLXrrLsiLDz1Q54PTyepdcuZf97jjnrulELqUd6dHfgNuKP05m9a1qdWeYBju4ZLXYrRZwHngLwChNoJhwRoVm0NP3vEZIdfael1gYe42n14SPfWL3zqWZ5Xlcf3nrrtaYsD6R5YrumBEiE3ww0hy/u3PI+pK5on4cK+xK1Q78aF3uU69wD75CS5rpot7pUiO/CzLHWQhIOG33eqjMibDYjnLYinO+9LmjmxigNZXS/1LZWv0rs9LJpN7gh1ryfaj/n1554riWcNaIlnjgpJFp7W+VpU59qeKXA5Pf1yiy/E77ogbBvrcE0qlo7IGfmLvU+4yLJNZm9KHM8E6s9CTWk3gVLuN+Tns1efZ9vDLX6PBdv8AbvtJMe3LHf8pyU9+/Wxx5eH0779KJNrc8i8+dOgSIR7qPvdxL2Gk84dxz0rnEQfioe8m/7eueg7UIuwrxCP616uKdSDxHARXJquhKSthiFUK35eD082lnjbfSwK79GM3es822G129nRoKh58nahZnrTHnIs64NXrZ4MjtszrnXImw193k7cyQ8XNiXiOjV5x++/rl2gpwxSaLCSdHXj8rphdG4uRW1kJN0cyjL5NTPJTqcLa1LYs/V2mMUipSOKen9O2BZl1qes9512nAsxwbNutezyys87e25PsuIEIQ7eMe0KTecfYJb1i5wxZhHisl5/dMHnzKtew3CvkTthtsxpDQzdMp9Hu6Gwrmv9p2kXVsh50chInTrXe5VFulVWMKA5byY7WWxHggj7Y0cd1q7a/PL2g4RbbjKnuh+m41X30jaRCH8VDyk5+9E+nv/NMmfrTI0eVxep1Hu954eFmk5kvusnhzN6VwrgDVww5ZSO6LK5SIdX73Q36x+19SRI33O0XyQzG3JdVZPqdYz1CpHcric5Z2UlKkt22PP4ZSn2Y/O+zzn0SN55SMiiYZcP3j+cixJk9oO1Zx8ok9LVq+dJ6smvCbHa6fFSePR29HT3avXu/4Z/eeRA+mRs3SWY8kN83iWHI5HlGNzvsl62+iNdQxL+lFSpuU+b1t2W1t79Gz1fo5Zf6b25KHFnwXuCd2rrvO/0hCCJR/B82SpLU9SrtVeabuW169KfI4Ir0Xnv94RqfDy8p5ZBN/5nVWASjyGq9I+PCIr579b8yf32O6SqiK15UriDOh56LDviTb86zkBuK5zawhGEz7lPm/rulZImGtLzaYRvYRwj/DhzLDvCC+RUCvvDOvk5XmF13KsaQURbb+rAPYK8a98Pmk4nFuOpizPcdRbQ7mh+lV4h293erYRM9aS2vzj1uXg9ROHfW8h/ohi/+6vJ9pBytnAvRZkL3IbpYtJS8Cc//a0z7tML1HlJf44341y8npYPcIRbb/rBk0UM09XbNZS8clZsyTl1crwEqGaQzZYx8w+kdblFO5Fzt8KagtUa2McbZi5y13rqm+FOmcsRFqbo0LGUpty7r5w94RgVEirTH2YkfskfRZNeoYWr+dfMZal4cPa9aN/j/rBa/3L7+WWw4307BQi3onoNrl7u9/G80d0He+fF7uGtaT0vBajZ+Qs/mf5Hp6KXp2lp3NkV69ebXtwyx/Z1gvv1u7XlFm7rnx2a4ja6oXt9eWM+ecVKbgiK7053vMsL7v83jNUXatndN/MMRLl9e7hnd4lKdfxJQ94/nrMfpsG2DknU++UZsnjq93POQ22vKsWO7h4ebB65Us80DXvjCU0x702qh294do5yxPRi0rc2RPSoiWe8vlhbbuW1zNC+F2FiPk58iR71TFZ+Km4lfi7Gx6LhwbPTWRk/0igRL4oMPpsVXne1NqwDNdzQnYnnGcb9ZslVLuibbnPvMoTt9N42w1NyDeK1QeEKzNL+F2FW4V9T2aGf1ckWrfguvGj3P21kJ01xFp+37vGk1F4YpS8bg39juwa1dsKqebs+IKANvymDX1xbDg/G72o0Lt+JhFhsCuEkSPnm6ZMzfonrYsz7ntjV5LSYp1PHikYM/CoM8DrJw77vu9twRX4+NVXwgVgLSy4ygtX+0yS80HkY7vl5CXNZ4pu656IOr/nlCHB67l23aC5bZa3g0dYLLJNS5F/1ulZd49aHT1vee9QsCr3a0be1+zyJGNgVFevrcq+7nnxW/8eIU3Z4KaTeNbrgVd9q8O9J7f0/J1ECkBPT0OtXE1yt8Senv3SBWv2prGruFnZJjnafrfUO/I0jjy62vu9kI6rlt2rvas5kvbnXhch1CK94p5o5pnkEL7DHMjxzEdc5fXTro0befxO8MKHhF0UuIboU4/VA8AJgUWwq/AjWmfbSICUi9oubVjznnuWrc0btLBL2xK9m7fJvd7KqB1HOXbafthlbGttX83Zfq121M4NSQi77HvPPo1u4930xi3DvjOwurY1p4yWd0KSG1XeV7OnVbe0fAmzQ0/c/DMPD+wMpGMgv36VN7UWMvVgZT8Q+XmEPO6dvXFy295TXJT3zvCg9coYpYxIyvLC4/ml91nXTs922UFcz+bWYd+cmSFgzULp6WL3Kq8nDi15I56Lc9R9msVoVYhmRkhlFOa05DeO8oS8x5rknsjctF6ItefNnTW2ynaQjJPyfmnd0vtr9502e62NXn3g2ZZe13Bs0NxflhNhoweWQ9kkjx/CvjsyCrtx2DHcFDUhy3BfL0ypLXeEJmwtKX+19ymafLHUjI1R6M+C1h7ts3iEi0tbvKiF0kbXe9g004PvZfMIyxgpy4mob0WqkHf53nVIPNGP5h2E+HshWp1bNhzt/T0bvARovrhabBzdq22LKJskXr7aPY+ykLSE+czn2zn/q2WbJO+NI840c8PrUCXB0u65RywqNHmledkbF5oDR7QHPbdJms8b1S9e3vrWM+yW55eDsG9BRPjXY0Msy/AO41psiqzfo45WyKdVpmf4xqMcC9p8PWm4KCc6nFeG63bZsHvjLP/MWqYn0nXEEjLfJVzX+yz/3DPXVZsmog2DzvJ2tuilBUWmTUiprSWSfu/ZP1n4IexrJbLDpKedXhnlCXpFGNErzFHCWSg8yvYqU1Ov9l7piV7qCbLYOcublnt/evVJQ5q1+3ei9SzWPtPco9kcPdszMoXDK0xvnfMeh92cmeKKK/A8I1tatI6UjYSfCnj+Gow8gJrTLPeU05r4I1G0iweEyOaZnHUa5Hj/rPV7P0ukh6Ism9sG3huVl4es1/aWftF4MSLquwK1dcCrzy229PqnpOxjbogzv3+Gh7W813MMPhLcvtaySPjB8+fFqAPzDYZ7wuOccmo5EbUyys9qOUK7eS52o+fBqH2Xtyu3faNPtpF9zPXweD3b6jw867Ury/RitP5I6c2v2bTmopcn3UtUWNu+tjftKPxWzINRnda14QoevxOIvw6SjpQKQE2ZGiGxaqPZZYGJRiIAPWh5mnYWFKvwSEuQeog49bUOFdZyPci9VFFjiuP1ijq8SttQamvL08l9lt79jzbHZ+8RHO9ryzNdlnN14UeEsC8LzksgXqG3vLwTjdgrB7GH6ARjj6x3CH7Uh7uf6CM9nq2yo+u/QlhNYkt0GEzaLrPTPqR1nWPPO3w9CpF75KJJ7ZkRsVgZ8pes173xsoHwQ9g3Ak7Heg9gz1yz/N+cE6QlLPLo5B6YmieW46HRtueqXCkLkWFTjocqynu0e7triXwuaarMTOHXu67mYY+eix7eaqs9ErGpDZfOFH61eauZy+U9Gwg/Ffjzbkw+fvWV0L8CUsNjYmhPj+V9O3k1dsCywEru4YQhrtYf3h6FVc/v7eWdyagPZqQqRMPx0Na+75WV3zO7/yVrMDc30dvJILk+OpLhvWfVyriq8CNC2FfMbAGopRX+PT/Lr9MsfmU5VyBysY4MkUQQZe+MduAu6tzQ8Kic3n3S+1cTmdLB6XtLOFrbR5rw/CgEK7FHgyUU31vnIsSVJIzaGiPWvvYQ+FI2E34I+0azWYe/pkwKroUny2tr37W4yubGwdv7EJ0gr2GUJF4u2JbQTavcSMocS+7hpHfdVZPqtXZr+qpVF6d+iWj36oeeaJux9s0eT9z6vOap1KvLqdd6EJvBrjpAAsSfghUdL1kQo7x0p6CctXCD+ZSh/tkiT7KBcATg6F6tDbUyvclFecumVmrGCK3ws8BdP1r3etYzEi2cMnrkfWZpN0k+MVdc1w5+rfK4Hkgvj6KH8Iv2+j2C8CNC2NfMzDCwNFxhuUZqV7kRe5Zfq69ktWeybNPV+WDSfqiNCc04WfHcXmNO4z30DK1pPCeRubiteR1pQ+Tz9OqLaLsTjmDzWoc5dZbXlte31rJWub3wL8cWiY3ccWipt8Xmog9h30fGy2VuOdG2OE+JowlpraPEejr3uOa0o3YPx2MjRSoSPEMtXGZ6fz3DRJK24ubKeoYwNd9p7aiNIYsNNZu8ytIwY4z2DuDe9Wu9ldLQbateSxmcOvK6VuT4PRrw/Dmx4ncB5uWeeJ7KW/dyvUTez+u5WXAXC+tpenSyl/bPqO255UThHSaOmjOR9WrbgON55Mz1md7BK2LxaHPvGY0fL097BNpIzkr7IyMOm3v8TuD5W4XkdwFGnPg430V4m3qTjptbwq2vrMOSlxOxSLTaIGoh3CG3MtqGVZugNW/Jq9+57bvDWJiB9Tkl7TnKt+zB7fuoPcGCNGpAtF64RnkfLyL8VED8LaJcXLwGa2SYaeRZaG14UnvO60cC0PqcnLKk9VlCLho0Hqp8jPRE/eh5a2VwhM9OG50Ey9zySMfQphDM3Iwjwpk5q8PBXvXnfd3yru0wT2o2lJ+tFn4lq6IFVwNh3yBeff4hEfHCfed13kmyLSG2Kox2otngotpj1v3akO7o+poQl7avNn9mVbh51WajTQHQpiT0DlLcg1lpi9QeyXqxY/qBdc6fRPX37DQOS3+27NlFbHnZcVFvH8K+u9GbQKVXycOrlSfG1mzZZZK2iLQv+iTd2jSinqk8MJyfeZTreV0kO3hHvOmNmdZ3lnFWi0JIxdOu/eBtV+nZ9PJ09vpbSstbJw1518bUKL1nJR5RlF3HcQQQf0E8ffDp659HYiefZB75Qr3T/6oBbnmmlYuLNG+xd40153Fkn+f10utyr8UVF1CJ3btsdiNqz1Q7MHiwU5t4P58kDKuZj61Dv9d6qVl3znFyxbms4XzOfN9+dBD2nUDrTWBvr815Wut5+HbLzxihDRdHPFtpS97eNXZp3108vruyOnwdPVZHdVjHbykYuHVpw+GjurzqyO/3XDd77W2dq7UxZfXiSu676jpz0VBvjjjsC/E3iVeffzgUZLNyO64kADk5jOd30YuWJDF7t3atCVVrbtIV4Ni9U66aF5r8xHx8WMK9Ee08q4967eYlzFrMSHmJyF2clQvLsYM7hnM7HsTbh5y/XXn64NNpIUOiergynxSRuWie5Hb22qC2MEvDHdrFefd2zEP+nuVdlZqQ6OXKju73sEcTYouww2NdiB4fM4Qft96InL+IfpXacbVcT8ka94DCTwU8fws43wQ+sXideteNFhRuvZZEcm49Vqwn9qt6tk4k4bYTa9jtCpRjcKWXjzM3vezReHpa13OjCL0yVhG1Bs0cR5o2tnr6JXZJ9y/vMZ7DHe8PKPoQ9r0KX3z6cncz0Oay9CaENLfHa8JyhYmmDkvIqrRNEi7YaYPjckXx5slqoTJT/Flt0pYzSs2YiTRHMDonblSet8DeKQcvKrXpRNIeD5DfVwPi74qUvxPQM+mYc2rfwSvmkYsmvVeDNW9pln079JPVFu9n2UH8Sbyzs2zzbOfaYSwnSrRw0z560ZDad70yZ681M+rU4LWuccvRjNcHFXw5yPm7Irm46y1itTw+btmt7zgTaMaCo815sQhFDZacqBXt6J33WF7furf2nbZsD1uJ5oqpK2zcvT7Tltcri3NYkORu9ern1De61pJDfCdGB5reOsEtp7yu1jeadebOQPxtwMevvvLO7wVsiQxtaNQ6ITwSwaPqkIjYOywM1ue0jBdNCGt0zejf52c9u0cHKy+uNL7KthilmvQYXSfNPeQiGatWMagtSwunf3bnXJvzfpqxj5x8/Oord/D6qUDYd0PyF0JGpyqPPJ1RPavxfE6i9ol+5zaQwvWClM/uEX7X5pjWyqpdX/teE47mXC8N/0ls59h5lTHp4enU5HDVsI5li02S8kZlrl6XPNIQolMZWm10Q8GHsO8jwHkTyRIeObGECC1I67J6snohgpZdVk/Eag+QZuO1iA2N11abtnB+5uEp5nh4tB4jq7frKsLPA43I6vXN6vnnwYxoCxer8It6llq5NxR+KuD525jWXwY58TpV9RbeiJNbtMexZXPPU+Bp01U8qkRzXjJoicqeZ0MrRGd5yzw9pI+A9bks99fG0Y5ev7zMyIgOt/5ePS1Pvncbe3Nj4Ye3fR+RkQgkkoUINCHBnqDaaXFtlV+rI7J+a1LzzJDfzHbw2jw0otWzTXfc+FbgMf/z9aX8OdIO7j2cddBKL5XB8zC6Y3ke3Fj0nSDs+4hwBrbErV5eOwqRnItk6x5tqHNkszV04x3Sk9ZtCUnNXFBn1rUyHCcR3hJmpkusppz35Weack406R1lqgZ3HeLWU3veE246idQWTjvkdUlTaDjlSmzlXhsFhJ8OeP4uRu9vBPdoeSrKz7l5SPmpOcoL0go9eNcTTa9Nd3iOWUnZ3uVHhuVq9A5MO/RjJCOBog1997DkTEa8nNEq3zoWRmvwKBWnV6dWrHGiQauB6HsLhH3vAicUXEMa8s2v44RJvBcHzqlYGupeuXBZw8GaOjQeL4s9KzYIie3cUHSN8toV6QqR9XGFBre9R/metbK5SMqQjm3uGPZeV0ZjaiT8tOHvWlnc71tlRo1biL4qEH93QisAT0YL3C6CiUh+Oh7ds+J5WrlDO3narIv2qjZuPXPteTSeldp1EUi8OhH1jYSwVpx5iYpeXVavVkS9EWW11m2tSGtdP2MdkR4UIfyaQPzdlVMIak64Uae51WKrtKHGDgIluk6vBG+ph5V7j4cdnmFJTf0je0rbvOrp1WvdfLVrScQY0Bz+ynt7h9xd1ieL164sg/ud9HA0Gwg+FhB/d6YmAInqp3nrxsIpY5dFRBvCnmWLpoxdQ6r5PbND2eXYXnn4WOU9lIRjo+yJTv0o0YSdud9H4rUWlPePDh7cuWSNAHDm6AgIPzYQf3cnDwWPFoYZAlB7TwSrBZ9HvtiMtuyJqKg6uTaN7FgZstXYUY4Dq3DawdvuyQyP7io8nmPkES3L95zLo5QKjj09IPxEQPyBN0h+SbQm5+NRFuEW0tASJ3fSWpbULk96uUYzbZAmxBOtH5/SHNXWNbvh3d6WeeQtbGplzFjzWqJqJHwttln6ceSZlZQHwacG4g+8y8gb2CMqhByNlyBptZcl/zEn0sPz6OK8xqOIv9X2cohoa634ixB+ZTle/ePh9ffMX7T2oyW/9QSizwzEH6jDyQfswcktkZZ1xQ3Oa1PxKlNS3xXa+1F5tL7QhrU969UKM83LDh65lJ6eOe+DbYn0YJu3k9Q2CD8XIP5An5YXkGhuqKC8P3Iz9HjJReol9QoBe3A1j9KIRxNSNbgpGbu9rOA91izhcumLIto2lYQ9teV75AvvBkSfKxB/gEeZD5jnlFgWvx1OtKN6tOFaiUdhdD3H++AZ2pkBNx/Pavvdxd/q5+eGX2f0s5fws8CJinDrXC3qZwDRFwLEH+Dz6vMP3/ksIkk6+l5pHSfanJZeiMNjY+Z4fa62OXgLlkfwZnok9te+j8ZT7Ejr0dhStqeHbSNGL4xcdcxa+OizT+jpg09Xm/GoQPwBHd5CcGc8E66jRMjqDV6CJB/0EUSblVHfcj1eu7WfxaaRWDqZ9cyeEZArzeUa1qgJBN8UIP6AjTuJQAs1Dw3RPdtK43155Hbi5utxrlnRThF1a4RR5NzyfMZZbyF7Ygmpc8uH6JsKxB/w44tPX3ZbvKJzc0Zv2UUsursu7ES6l09GSeXWNxlXbfTRlG875nA8WZ7PXvYRd+xHvYjE7fPV3k2Plz1mev49Uli087gH8vmWAfEH/LF4A3virrf5S+uzbDIWRqJ2Vm5Rre6aXa17e9e1rq8Jl1EZeTmj/r+CAJSO20jxx51rnPu92v4u3l4r1gMEtw7J9RIg+pYD8QdiaP2ewBOpOGvdqxGBrRBsfn30JlTzuNTsiKq3Vg9HUEns5Io/TjneYaZVzHjBR2tLDlf8XaHNHxmJ8NtB9BFB+G0CxB+YQ+1XxeRwRMSqjcYS4okKC3kky5dIQ38WpCJIYlOUJ4RTHyc0PdPLOwqp1+zVhOFb1/SuXUnLczZr3NfGpNc43K2tIfa2BOIPzKP2t4Otb4ZJ75uJt42zRIv1Wkn9EkHnIf4iwsXc/NRVwsKz7JOakLV6ErX2eOcV7wS3T3ddDyH8tgXiD8ynJgI1rFq8y7CxNPRprTuqvJnPEQ33ZZPWNZq6yvK4HsgdDwac8ktB7i3+IsP9VxjPo4PKznmvEH3bA/EH1uIlBInaLzNYPQNeosiy2ZfPphFtu4uMq5atqXelF7d1v8az2gsRe4vqFjPHvxfccVF+zylXI5g9vPAQfJcC4g/sQysvcFUYjVvnjCR+zYstrfs8Q0jSTd4jr0zLKiHA2ZCjhJL0xSrJQYdTPrduzZgf2bCjd887H1VaXn5NT4ByRT0E32WB+AN7kb8lXMsjan0Wyaz6LAn1retr90k3FM49nkjbwVN8Rt3bKi/CG5sjORRI7ylDkJIXWiR5mT2vGMfWHUTfSaTgb80TiSCUAOF3aSD+wN7UfmfgSXlS9dpMpfRCsvnnmvKk93pTe7ZV9mi9n5IyuYJ7dhu0Uhp6aASVZZxyy9jlBZnd4M4t7RzUHF7yPsFf4HgoIP7ANdj9z8hxN2KPMOlKAbaSnhdY2h5SzxdXmHh7dnqebm3ItiYKPYQzx46WPRGeuruLSQsQfQ8PxB+4JuWfkttxYe9trJZcqxWh752IDJdqPS+R3sGRcMs/93jxI8czP1RTvgZNHVq7HvEghnDuLYD4A9em97awdWH2fmPRIhg4ye+PtgnNwEu0eZRjGW+jxH7tuNol5YBIJuK80ix2P2B6AtF3KyD+wGPQ+gXSHnlq3iKQiJ+c3StLGr67irdwtjfFq10s+ZFWGzRvxmpfHprBipesuDmSvUPc6nYjkrcdRN8tgfgDj0kuBncTPR5vjXrnurU2txUb2279JcESSpbeI7GhVab3G8weWHIJverTvrQSDefFH4g9wADiDzw++a+POYlerL3r0rzZqfmVHZZrPXkkAbiqfTmiZsacmJmDJ8FDqM96A94anofgAwUQf+A+1LyB3uHg/DPPX/ki+fUe3uLjykJsdzxfFNEIhOi+1eQi7ipIZ6PN883bFaIPNID4A/fF+qflWhtbxMsZo7qsXj7L246a+7hEeFV2yAEty2uh7YfRvTM9a5xcup3RzqHetVFA7AEmEH8AENmEIPdtwRKrh0fj5bO8aNKzhWtDBDv+uo3Rr+fxFuH5vZrQ5Qw8X2rxeAHLK69v9RyA4AMKIP4AaGH1DOb0xECkF2mmh8uyOUe87To7D4sr+mvXz2CVSPHyKp/3Sp5jtzdxR4xedoHQA05A/AEwYiQCZ/5KD02ds/DMcVxxvxRLqDUyt46T2xdhg2edml9d41HvSnprA0QfcAbiDwAJnt7AKMqXTjiCZAfvo8ev7cBnygAAB11JREFU+NDmRnLLG3n3VngbPUL/rTpWpQPs+CtzkL8HHgiIPwCszBKEGgHDfQN09UZXEv1772aFs6W0npsTAr2y0NkBrvc0ImUDQg9MBuIPAE+ihaCXR6S818PrJql/db7bjuImypM4y9u3ksg0iKixA8EHFgLxB0AUs4VgKeR2e2lCm6gfYYuEWfmVoxCzxeN7NUGnFaMRaQyeQPCBTYD4A2A2K7yDOZK8tSiPEHeTjvZI7gBH6EYJ50ixxBGgM/pXMsa9bIDIA5sD8QfASlpCUPsmqVY07ORpa7GDN0sqhr1+t5+3+NOK71EuHKf+FYJvxpiH4AMXAuIPgJ2oiUEvAccNHe6e/7Xyd+itDF17CjZuXbVyeuJtt7ejW3iE8iH2wIWB+APgCniFinf8/YAjemK0JwSjQ4fRda5iJDJ3GUPevxi8dT9EHnhAIP4AuCJX+H2DPaQb9+rcvzuJv9loPNI51r8YUgKxB24AxB8Aj8jVxKE0XFhe0yuXm4vmyQqxGvW7EXcUtJr2PccCxB0AEH8A3IIriEHNnyhrXcupQ1sGt0zveqT1e/yqH005UUjsy699+uDTMJsAuCgQfwDckVeff9jcNK+42VvKt5Yn+dU6EUS94GEtcxal7RB7AAyB+AMAtHn1+YfN71YKgt3DkqvhCucd/j5xrd7aNRB1ALgB8QcAkGEJIbc8iz2Po5c3UvpnzkquJjJXeXE9fhUP8vIACAXiDwDgg1UU1lj1J812z3+r4eXFk4hGTR5e/h1EHgBLgPgDAKxhVkhZms+m9ZjNzptbFfru5UpCzAFwCSD+AAB70hKHXiFgKZrfRbezp1DC+VzIuwPgIYD4AwBcmzLczM05ixKAVwQeOwBuBcQfAAAAAMCNEIu/9yOsCED8YAAAAAAA4F2eVhsAAAAAAADmAfEHAAAAAHAjIP4AAAAAAG4ExB8AAAAAwI2A+AMAAAAAuBEQfwAAAAAANwLiDwAAAADgRkD8AQAAAADcCIg/AAAAAIAbAfEHAAAAAHAjIP4AAAAAAG4ExB8AAAAAwI2A+AMAAAAAuBEQfwAAAAAANwLiDwAAAADgRkD8AQAAAADcCIg/AAAAAIAbAfEHAAAAAHAj3l9tAACz+Lr09cdP0U8+/yMlSucX6fX/vSFl/65+X3yWih9ScXFxe7+MZ47WfYzyj0p5rXoopefru9dUymfWcTQ+735GwmcoPh/WOapbWm/1u0NtP5U9UhsKbw2dt69PFVvyIZqK78p78vLyYXp+nt///P3bn791f/ZMve/fLrsoL7393dv1F98Vz1HW+ebro1l++Sxvfs7LzspIRNmK8s7153f5FYmIfvhHfvKj4zh+BQEwGYg/cBt+in6Sfvn73/r8j/RE6SkRpSeip/QsBp+e3uxE589PiVJ2zTvfn/97/f2LMz2/Pv+eiI6Unn3u2fevPyN697Pz5zc7FR2vPz+vfymXiOiJXl9/fv/8X3pt2/nvt76nymeJXuqiF9uL7/L/UvlZeuvf737/7r1vfZ599qburDxG2e+UT+3733x+vHVNec/r798p88iuyf5dfJ/Ka5+vePnu7f++dDdROiil42VoZNecVWXfP73++Y2oOz97oreveR4OL9+9/v7l5+z787un2nevP3/11r+fv39F773+7Pn7915+Pu9//v787uW/5/f0/PNTevX6vvz+97L73kuvXl//fO9z3c/XvHq55rkuIqL3Xj4773mPjpd7zrKPN2W//vn5f8+2v1xDRO8lovcovZ7C76VET5RePksv/36ip5dOfS89/+u9r//06wiABSDsCwAAAABwIyD+AAAAAABuBMQfAAAAAMCNgPgDAAAAALgREH8AAAAAADcC4g8AAAAA4EZA/AEAAAAA3AiIPwAAAACAGwHxBwAAAABwIyD+AAAAAABuRDqOd/6qJwAPSUrpfyain1htBwAAvPAzj+P4+1cbAe4H/rYvuBM/cRzHN602AgAAiIhSSv/DahvAPUHYFwAAAADgRkD8AQAAAADcCIg/cCf+g9UGAABABtYksAS88AEAAAAAcCPg+QMAAAAAuBEQfwAAAAAANwLiDzw8KaV/NqX0j6SUvnO1LQAAkJNS+qWrbQD3A+IPPDQppV9BRF89juMHiOgnUkr/0GqbAAAgpfQ1KaVvI6LvWm0LuB8Qf+DR+RYi+uzl579IRN+8zhQAAHjmOI6vHsfxh+jN+gTANCD+wKPztUT0/738/NeJ6G9daAsAAACwHIg/8OgcRPTey8/vEcY8AACAm4ONEDw6/xcR/U0vP/9sIvqrC20BAAAAlgPxBx6dP0pEv+Dl519ARD+40BYAAABgORB/4KE5juMjIvqalNKvJKIfP44D4g8AsAUppV9NRH9PSunXp5S+drU94D7gz7sBAAAAC0gpfQ0R/Q16zkf+6QMbMpgExB8AAAAAwI1A2BcAAAAA4EZA/AEAAAAA3AiIPwAAAACAGwHxBwAAAABwIyD+AAAAAABuBMQfAAAAAMCN+P8BI7dtaYlVPnsAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 612x388.8 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"mask_4fgl = mask > 0\n",
"hp.mollview(mask_4fgl)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [],
"source": [
"np.save(\"/Users/smsharma/Desktop/mask_4fgl_Paddy.npy\", mask_4fgl)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment