eteq/check_sdss_cmd.ipynb

## check_sdss_cmd.ipynb
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/erik/miniconda3/lib/python3.5/site-packages/astroquery/sdss/__init__.py:29: UserWarning: Experimental: SDSS has not yet been refactored to have its API match the rest of astroquery (but it's nearly there).\n",
      "  warnings.warn(\"Experimental: SDSS has not yet been refactored to have its API \"\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "\n",
    "from astropy import units as u\n",
    "from astropy.coordinates import SkyCoord, Angle\n",
    "from astropy.table import Table, QTable\n",
    "from astropy.io import fits\n",
    "\n",
    "from astroquery import sdss\n",
    "\n",
    "first = lambda x:next(iter(x))\n",
    "\n",
    "%matplotlib inline\n",
    "from matplotlib import style, pyplot as plt\n",
    "from mpl_toolkits.mplot3d import Axes3D"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "qry_template= \"\"\"SELECT  p.objId  as objID,\n",
    "p.ra, p.dec, p.type, p.flags, p.specObjID, dbo.fPhotoTypeN(p.type) as phot_sg,\n",
    "p.modelMag_u as u, p.modelMag_g as g, p.modelMag_r as r,p.modelMag_i as i,p.modelMag_z as z,\n",
    "p.modelMagErr_u as u_err, p.modelMagErr_g as g_err, p.modelMagErr_r as r_err,p.modelMagErr_i as i_err,p.modelMagErr_z as z_err,\n",
    "p.psfMag_u as psf_u, p.psfMag_g as psf_g, p.psfMag_r as psf_r, p.psfMag_i as psf_i, p.psfMag_z as psf_z,\n",
    "p.extinction_u as Au, p.extinction_g as Ag, p.extinction_r as Ar, p.extinction_i as Ai, p.extinction_z as Az,\n",
    "ISNULL(s.z, -1) as spec_z, ISNULL(s.zErr, -1) as spec_z_err, ISNULL(s.zWarning, -1) as spec_z_warn, s.class as spec_class, s.subclass as spec_subclass\n",
    "\n",
    "\n",
    "\n",
    "FROM dbo.fGetNearbyObjEq({ra.deg}, {dec.deg}, {rad.arcmin}) n, PhotoPrimary p\n",
    "LEFT JOIN SpecObj s ON p.specObjID = s.specObjID\n",
    "WHERE n.objID = p.objID\n",
    "\"\"\"\n",
    "\n",
    "radius = Angle(30*u.arcmin)\n",
    "def get_qry(coord, radius=30*u.arcmin):\n",
    "    rad = Angle(radius)\n",
    "    return qry_template.format(ra=obj1800235p127.ra, dec=obj1800235p127.dec, rad=rad)\n",
    "\n",
    "def do_qry(coord, radius=30*u.arcmin):\n",
    "    tab = QTable(sdss.SDSS.query_sql(get_qry(coord, radius)))\n",
    "    tab['coords'] = SkyCoord.guess_from_table(tab, unit=u.deg)\n",
    "    tab['sep'] = coord.separation(tab['coords']).to(u.arcmin)\n",
    "    return tab"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "obj1800235p127 = SkyCoord(178.33*u.deg, 4.96*u.deg)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "&lt;QTable length=24800&gt;\n",
       "<table id=\"table4529391432\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>objID</th><th>ra</th><th>dec</th><th>type</th><th>flags</th><th>specObjID</th><th>phot_sg</th><th>u</th><th>g</th><th>r</th><th>i</th><th>z</th><th>u_err</th><th>g_err</th><th>r_err</th><th>i_err</th><th>z_err</th><th>psf_u</th><th>psf_g</th><th>psf_r</th><th>psf_i</th><th>psf_z</th><th>Au</th><th>Ag</th><th>Ar</th><th>Ai</th><th>Az</th><th>spec_z</th><th>spec_z_err</th><th>spec_z_warn</th><th>spec_class</th><th>spec_subclass</th><th>coords</th><th>sep</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th>deg,deg</th><th>arcmin</th></tr></thead>\n",
       "<thead><tr><th>int64</th><th>float64</th><th>float64</th><th>int64</th><th>int64</th><th>int64</th><th>bytes6</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>float64</th><th>int64</th><th>bytes6</th><th>bytes19</th><th>object</th><th>float64</th></tr></thead>\n",
       "<tr><td>1237671142020219643</td><td>178.32879731</td><td>4.96462952148</td><td>3</td><td>68987912192</td><td>0</td><td>GALAXY</td><td>22.65771</td><td>21.9395</td><td>21.12205</td><td>20.8422</td><td>20.54399</td><td>0.2825266</td><td>0.08393663</td><td>0.0612877</td><td>0.06515449</td><td>0.1543406</td><td>23.03028</td><td>22.73921</td><td>21.96251</td><td>21.72451</td><td>21.20847</td><td>0.1049995</td><td>0.07725767</td><td>0.0560337</td><td>0.04248866</td><td>0.03012499</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.328797309864,4.96462952148232</td><td>0.28692367495620213</td></tr>\n",
       "<tr><td>1237671142020219091</td><td>178.325317321</td><td>4.96270838587</td><td>3</td><td>316728336843024</td><td>0</td><td>GALAXY</td><td>24.87228</td><td>24.01616</td><td>22.38811</td><td>23.12548</td><td>22.34377</td><td>0.7724417</td><td>0.4432389</td><td>0.164077</td><td>0.4081114</td><td>0.5367754</td><td>24.80732</td><td>24.07956</td><td>22.99819</td><td>23.35466</td><td>22.45049</td><td>0.1051598</td><td>0.07737558</td><td>0.05611922</td><td>0.04255351</td><td>0.03017097</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.325317321334,4.96270838587141</td><td>0.32365997771680666</td></tr>\n",
       "<tr><td>1237671142020219092</td><td>178.323746</td><td>4.96402411303</td><td>3</td><td>193589476982800</td><td>0</td><td>GALAXY</td><td>21.9166</td><td>20.21304</td><td>18.80011</td><td>18.2708</td><td>17.86793</td><td>0.2013644</td><td>0.02702857</td><td>0.01253898</td><td>0.01146572</td><td>0.0226814</td><td>22.73883</td><td>21.26597</td><td>19.87917</td><td>19.43731</td><td>18.87523</td><td>0.1051809</td><td>0.07739112</td><td>0.05613049</td><td>0.04256206</td><td>0.03017703</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.32374599952,4.96402411303364</td><td>0.4450260723158584</td></tr>\n",
       "<tr><td>1237671142020219616</td><td>178.326423263</td><td>4.95492414353</td><td>3</td><td>68988044032</td><td>0</td><td>GALAXY</td><td>22.90097</td><td>22.55076</td><td>21.35245</td><td>21.00606</td><td>20.68048</td><td>0.2817276</td><td>0.1147992</td><td>0.0579804</td><td>0.06061608</td><td>0.1530485</td><td>23.29821</td><td>22.95922</td><td>21.82498</td><td>21.55735</td><td>21.09591</td><td>0.1052939</td><td>0.07747427</td><td>0.0561908</td><td>0.04260778</td><td>0.03020945</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.326423262661,4.95492414352573</td><td>0.372105631901309</td></tr>\n",
       "<tr><td>1237671142020219994</td><td>178.331466288</td><td>4.95221141823</td><td>3</td><td>281543964623104</td><td>0</td><td>GALAXY</td><td>24.80282</td><td>23.58797</td><td>22.56713</td><td>22.44884</td><td>21.65725</td><td>0.5354075</td><td>0.2248108</td><td>0.1367196</td><td>0.1655365</td><td>0.260006</td><td>24.72056</td><td>23.79318</td><td>22.76958</td><td>22.62515</td><td>21.80913</td><td>0.1051695</td><td>0.07738271</td><td>0.05612439</td><td>0.04255743</td><td>0.03017375</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.331466287681,4.95221141823446</td><td>0.47546340651109925</td></tr>\n",
       "<tr><td>1237671142020220012</td><td>178.335168491</td><td>4.96292771944</td><td>3</td><td>281543964758272</td><td>0</td><td>GALAXY</td><td>24.16623</td><td>25.29726</td><td>22.725</td><td>23.52194</td><td>22.79156</td><td>0.5560204</td><td>0.5022258</td><td>0.1684994</td><td>0.3984331</td><td>0.4509773</td><td>24.62042</td><td>24.96176</td><td>23.02436</td><td>23.07375</td><td>22.47514</td><td>0.1048441</td><td>0.07714327</td><td>0.05595073</td><td>0.04242574</td><td>0.03008038</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.335168491305,4.96292771944105</td><td>0.35539572629502586</td></tr>\n",
       "<tr><td>1237671142020219615</td><td>178.317003361</td><td>4.95822916419</td><td>6</td><td>123214323908624</td><td>0</td><td>STAR</td><td>23.32777</td><td>22.98223</td><td>21.81454</td><td>21.5239</td><td>21.10142</td><td>0.3153538</td><td>0.1296772</td><td>0.06708524</td><td>0.07051217</td><td>0.1579347</td><td>23.06081</td><td>23.00879</td><td>21.86376</td><td>21.58507</td><td>21.15269</td><td>0.1055713</td><td>0.07767835</td><td>0.05633881</td><td>0.04272002</td><td>0.03028902</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.317003360724,4.95822916418773</td><td>0.7841112951278447</td></tr>\n",
       "<tr><td>1237671142020219981</td><td>178.32244716</td><td>4.94928679774</td><td>6</td><td>281543964754176</td><td>0</td><td>STAR</td><td>24.26053</td><td>24.22479</td><td>22.80974</td><td>22.63634</td><td>22.53944</td><td>0.4671116</td><td>0.3091245</td><td>0.1454135</td><td>0.1667313</td><td>0.3482596</td><td>24.16979</td><td>24.21153</td><td>22.78297</td><td>22.60784</td><td>22.50568</td><td>0.1055523</td><td>0.07766435</td><td>0.05632866</td><td>0.04271232</td><td>0.03028357</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.322447159862,4.94928679774194</td><td>0.7855019319828696</td></tr>\n",
       "<tr><td>1237671142020219991</td><td>178.320466426</td><td>4.95769301608</td><td>3</td><td>281543964754176</td><td>0</td><td>GALAXY</td><td>24.22022</td><td>23.96807</td><td>22.65197</td><td>22.00484</td><td>21.83047</td><td>0.6756395</td><td>0.3776841</td><td>0.1836063</td><td>0.1495008</td><td>0.3621331</td><td>24.27066</td><td>24.37602</td><td>23.07774</td><td>22.49855</td><td>22.51784</td><td>0.1054515</td><td>0.07759023</td><td>0.05627491</td><td>0.04267156</td><td>0.03025467</td><td>-1.0</td><td>-1.0</td><td>-1</td><td></td><td></td><td>178.320466426481,4.95769301607523</td><td>0.5864430956914526</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>1237671142557221005</td><td>178.081853462</td><td>5.38133697946</td><td>3</td><td>35255507619856</td><td>945883779803670528</td><td>GALAXY</td><td>18.21699</td><td>16.85933</td><td>16.20369</td><td>15.84113</td><td>15.57382</td><td>0.03288162</td><td>0.005410521</td><td>0.004512412</td><td>0.004550104</td><td>0.01044649</td><td>21.49107</td><td>20.34904</td><td>19.4799</td><td>19.1081</td><td>18.9646</td><td>0.0665713</td><td>0.04898253</td><td>0.03552621</td><td>0.02693845</td><td>0.0190997</td><td>0.0759038</td><td>1.768862e-05</td><td>0</td><td>GALAXY</td><td>STARFORMING</td><td>178.081853461525,5.3813369794599</td><td>29.30808935468947</td></tr>\n",
       "<tr><td>1237671142557155465</td><td>178.042175303</td><td>5.26205585697</td><td>3</td><td>105658611531864</td><td>945886803460646912</td><td>GALAXY</td><td>19.38679</td><td>18.26323</td><td>17.74543</td><td>17.35562</td><td>17.24966</td><td>0.04905115</td><td>0.01246719</td><td>0.009075034</td><td>0.008557597</td><td>0.02494653</td><td>21.50059</td><td>20.75753</td><td>20.07365</td><td>19.68711</td><td>19.74026</td><td>0.07829967</td><td>0.05761215</td><td>0.04178514</td><td>0.0316844</td><td>0.02246464</td><td>0.1179374</td><td>7.379925e-06</td><td>0</td><td>GALAXY</td><td>STARFORMING</td><td>178.042175303131,5.26205585697271</td><td>24.986461827541287</td></tr>\n",
       "<tr><td>1237671142557155481</td><td>178.028095756</td><td>5.29936368886</td><td>6</td><td>68987912704</td><td>3655870471926292480</td><td>STAR</td><td>19.36934</td><td>18.00363</td><td>17.41057</td><td>17.16948</td><td>17.03721</td><td>0.02601778</td><td>0.006407063</td><td>0.005529598</td><td>0.005722572</td><td>0.01017298</td><td>19.36505</td><td>18.00363</td><td>17.39515</td><td>17.17069</td><td>17.0165</td><td>0.07393174</td><td>0.05439827</td><td>0.03945416</td><td>0.0299169</td><td>0.02121145</td><td>0.0002282232</td><td>8.827133e-06</td><td>0</td><td>STAR</td><td>F9</td><td>178.028095756444,5.29936368886114</td><td>27.204887455707308</td></tr>\n",
       "<tr><td>1237674599997243524</td><td>178.196246238</td><td>5.40427591448</td><td>6</td><td>68987912704</td><td>3655872396071641088</td><td>STAR</td><td>19.6337</td><td>18.03608</td><td>17.42512</td><td>17.17045</td><td>17.08224</td><td>0.03471343</td><td>0.005930416</td><td>0.00516722</td><td>0.005612459</td><td>0.01158681</td><td>19.61094</td><td>18.06244</td><td>17.40708</td><td>17.14714</td><td>17.08634</td><td>0.07856219</td><td>0.0578053</td><td>0.04192523</td><td>0.03179063</td><td>0.02253995</td><td>0.0002329271</td><td>8.470722e-06</td><td>0</td><td>STAR</td><td>F9</td><td>178.196246238395,5.40427591448054</td><td>27.828949343294976</td></tr>\n",
       "<tr><td>1237671142557024333</td><td>177.863617084</td><td>4.97750808348</td><td>3</td><td>35253360136464</td><td>5365108589491580928</td><td>GALAXY</td><td>24.45719</td><td>21.82636</td><td>20.2683</td><td>19.48768</td><td>19.11165</td><td>1.087922</td><td>0.08223589</td><td>0.0315064</td><td>0.02298855</td><td>0.0510853</td><td>24.55219</td><td>22.5255</td><td>21.0977</td><td>20.29562</td><td>19.86438</td><td>0.1089869</td><td>0.08019155</td><td>0.0581616</td><td>0.04410218</td><td>0.031269</td><td>0.4195722</td><td>1.570988e-05</td><td>0</td><td>GALAXY</td><td>AGN BROADLINE</td><td>177.863617083733,4.97750808348156</td><td>27.897601541129657</td></tr>\n",
       "<tr><td>1237671142557024575</td><td>177.893712396</td><td>4.95316554845</td><td>3</td><td>68987912448</td><td>5365109414125301760</td><td>GALAXY</td><td>21.16026</td><td>19.80788</td><td>18.32243</td><td>17.77203</td><td>17.40221</td><td>0.2264779</td><td>0.03268445</td><td>0.01398877</td><td>0.01241305</td><td>0.02713966</td><td>23.29521</td><td>21.79142</td><td>20.1517</td><td>19.55021</td><td>19.28202</td><td>0.110814</td><td>0.08153591</td><td>0.05913664</td><td>0.04484152</td><td>0.0317932</td><td>0.3002202</td><td>7.056786e-05</td><td>0</td><td>GALAXY</td><td></td><td>177.893712395721,4.95316554845204</td><td>26.082588489241488</td></tr>\n",
       "<tr><td>1237671142557155672</td><td>178.069885236</td><td>5.30598820756</td><td>6</td><td>68987912704</td><td>5365118759974137856</td><td>STAR</td><td>22.46594</td><td>21.68303</td><td>21.51815</td><td>21.514</td><td>21.16424</td><td>0.2131306</td><td>0.04262485</td><td>0.0526799</td><td>0.0704174</td><td>0.1544516</td><td>22.3522</td><td>21.70647</td><td>21.51176</td><td>21.52892</td><td>21.20439</td><td>0.07079662</td><td>0.05209148</td><td>0.03778109</td><td>0.02864826</td><td>0.02031197</td><td>2.287982</td><td>0.001234605</td><td>0</td><td>QSO</td><td></td><td>178.069885236071,5.30598820755819</td><td>25.934005836537974</td></tr>\n",
       "<tr><td>1237671142020350606</td><td>178.399101342</td><td>5.2060694252</td><td>6</td><td>68987912448</td><td>5457349098158616576</td><td>STAR</td><td>24.21147</td><td>21.45792</td><td>21.2584</td><td>21.5177</td><td>21.1074</td><td>0.5453206</td><td>0.03804974</td><td>0.04327755</td><td>0.06922841</td><td>0.1379617</td><td>24.19868</td><td>21.4653</td><td>21.28599</td><td>21.57843</td><td>21.2114</td><td>0.08113739</td><td>0.05970012</td><td>0.0432995</td><td>0.0328327</td><td>0.0232788</td><td>3.066346</td><td>0.0001556588</td><td>0</td><td>QSO</td><td>BROADLINE</td><td>178.399101341775,5.20606942519762</td><td>15.3308708567223</td></tr>\n",
       "<tr><td>1237654605870530735</td><td>178.223507451</td><td>5.18908594275</td><td>3</td><td>69591892296</td><td>5457360093274894336</td><td>GALAXY</td><td>22.47614</td><td>19.7799</td><td>18.26422</td><td>17.68018</td><td>17.32185</td><td>0.5374746</td><td>0.02487411</td><td>0.01088206</td><td>0.009298238</td><td>0.02455113</td><td>23.24504</td><td>20.91237</td><td>19.50935</td><td>18.89311</td><td>18.52195</td><td>0.08873758</td><td>0.06529226</td><td>0.04735539</td><td>0.03590816</td><td>0.02545933</td><td>0.2762799</td><td>4.641956e-05</td><td>0</td><td>GALAXY</td><td></td><td>178.223507450825,5.18908594274838</td><td>15.147153090894598</td></tr>\n",
       "<tr><td>1237671142557221579</td><td>178.148570927</td><td>5.41030232073</td><td>3</td><td>35253360001296</td><td>5457362567176056832</td><td>GALAXY</td><td>25.71203</td><td>22.76553</td><td>20.94545</td><td>19.82676</td><td>19.43449</td><td>0.7300709</td><td>0.1923033</td><td>0.05425442</td><td>0.02974834</td><td>0.06576987</td><td>25.29397</td><td>23.5197</td><td>21.7002</td><td>20.5754</td><td>20.27556</td><td>0.06945611</td><td>0.05110514</td><td>0.03706571</td><td>0.02810581</td><td>0.01992737</td><td>0.5517851</td><td>0.0001713383</td><td>0</td><td>GALAXY</td><td></td><td>178.148570926859,5.41030232072854</td><td>29.112038895915756</td></tr>\n",
       "</table>"
      ],
      "text/plain": [
       "<QTable length=24800>\n",
       "       objID              ra      ...         sep        \n",
       "                                  ...        arcmin      \n",
       "       int64           float64    ...       float64      \n",
       "------------------- ------------- ... -------------------\n",
       "1237671142020219643  178.32879731 ... 0.28692367495620213\n",
       "1237671142020219091 178.325317321 ... 0.32365997771680666\n",
       "1237671142020219092    178.323746 ...  0.4450260723158584\n",
       "1237671142020219616 178.326423263 ...   0.372105631901309\n",
       "1237671142020219994 178.331466288 ... 0.47546340651109925\n",
       "1237671142020220012 178.335168491 ... 0.35539572629502586\n",
       "1237671142020219615 178.317003361 ...  0.7841112951278447\n",
       "1237671142020219981  178.32244716 ...  0.7855019319828696\n",
       "1237671142020219991 178.320466426 ...  0.5864430956914526\n",
       "                ...           ... ...                 ...\n",
       "1237671142557221005 178.081853462 ...   29.30808935468947\n",
       "1237671142557155465 178.042175303 ...  24.986461827541287\n",
       "1237671142557155481 178.028095756 ...  27.204887455707308\n",
       "1237674599997243524 178.196246238 ...  27.828949343294976\n",
       "1237671142557024333 177.863617084 ...  27.897601541129657\n",
       "1237671142557024575 177.893712396 ...  26.082588489241488\n",
       "1237671142557155672 178.069885236 ...  25.934005836537974\n",
       "1237671142020350606 178.399101342 ...    15.3308708567223\n",
       "1237654605870530735 178.223507451 ...  15.147153090894598\n",
       "1237671142557221579 178.148570927 ...  29.112038895915756"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tab = do_qry(obj1800235p127)\n",
    "tab"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "cmd_mags = {band:tab['psf_' + band] - tab['A' + band] for band in 'ugriz'}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$24.083189 \\; \\mathrm{{}^{\\prime}}$"
      ],
      "text/plain": [
       "<Quantity 24.08318915758459 arcmin>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "central_radius = 4*u.arcmin\n",
    "\n",
    "outer_size_factor = 20 #20x stars in the outer annulus\n",
    "outer_radius = 30*u.arcmin\n",
    "inner_radius = np.sqrt(outer_radius**2 - outer_size_factor*central_radius**2)\n",
    "inner_radius"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x1165f27f0>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEOCAYAAABiodtuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGWlJREFUeJzt3X+QXXV5x/H3hwLKj5QUWUOIobFOUBGNyC3jiFFHrXUo\nGdRWK2MRR6apVVtwGB2KjrY6dlprcbTTjk0LFdoMZRiwisXa6GiBGUV304RfQaIOaEgkUUyIE1Ro\nnv5xz47L3r17f+z9nvM953xeM5ncvT+fPXfvec73+T7nexURmJmZzXVE1QGYmVl+nBzMzKyHk4OZ\nmfVwcjAzsx5ODmZm1sPJwczMeiRNDpKulrRX0t1zrrte0rbi3wOStqWMwczMRqeU5zlIehnwU+Da\niDhjgdv/FjgQER9OFoSZmY3syJRPHhG3Slqz0G2SBLwJeGXKGMzMbHRJk8MA64GHI2JnvztI2ghs\nBDjuuOPOes5znlNWbGZmjTAzM/OjiJga9XFVJocLgOsWu0NEbAI2AXQ6nZieni4jLjOzxpD04DiP\nqyQ5SDoSeANwVhWvb2Zmi6uqlfXVwH0Rsaui17cG2L3/Ma7/1vfZvf+xqkMxa5zUrazXAV8Hni1p\nl6SLi5vezICSktkgt+3cx83b93Dbzn1Vh2LWOKm7lS7oc/3bUr6utcP6tVNP+t/MJqfKCWmzJTll\n+TH8/m+eWnUYZo3k5TPMzKyHk4OZmfVwcjCbIHdQWVM4OZhNkDuorCk8IW02Qe6gsqZwcjCbIHdQ\nWVO4rGRmZj2cHKwx2jwZ3Obf3dJwcrDGaPNkcJt/d0vDcw7WGG2eDG7z725pJP2a0Eny9zmYmY1O\n0kxEdEZ9nMtKZmbWw8nBzMx6ODmYmVkPJwdrBLdymk2Wk4M1gls5LYU2H3S4ldUawa2clsLsQQfQ\numVRkiYHSVcD5wF7I+KM4roXAp8Gngo8AbwzIr6ZMg5rPq9pZCm0+aAjdVnpM8Br5133MeAvIuKF\nwAeLn81ap80li7qYPeg4ZfkxVYdSuqTJISJuBR6ZfzXwq8XlE4DdKWMwy5XnSSxnVcw5XAp8SdLH\n6Sanl1QQg1nl2lyysPxV0a30x8B7ImI18B7gqn53lLRR0rSk6X37fHRlzdLmkoXlr4rkcBFwU3H5\nBuDsfneMiE0R0YmIztSUj67MzMpSRXLYDby8uPxKYGcFMZiZ2SJSt7JeB7wCOEnSLuBDwB8Cn5R0\nJPAzYGPKGMzMbHRJk0NEXNDnprNSvq6ZmS2Nl88wM7MeTg7Wej4ZzayXk4O1nk9GM+vlhfes9Xwy\nmlkvJwdrPS/aZ9bLZSUzM+vh5GBmZj2cHMzMrIeTg5mZ9XByMLPa8Dkp5XFyMLPa8Dkp5XErq5nV\nhs9JKY+Tg5nVhs9JKY/LSmZm1sPJwczMejg5mJlZDycHM7M+2tw66+RgZhPXlJ1qm1tn3a1kZhM3\nu1MFat1d1ObW2aTJQdLVwHnA3og4o7huHfBp4HjgAeAtEfFoyjjMrFxN2am2uXU2dVnpM8Br5133\nz8DlEfF84LPAexPHYGYlm92pnrL8mKpDsTElTQ4RcSvwyLyrTwNuLS5vAX43ZQxmZja6Kiak7wHO\nLy6/EVjd746SNkqaljS9b1/7JoRs6ZoyMWpWtiqSw9uBd0qaAZYBv+h3x4jYFBGdiOhMTdW7dmnV\naHO3idlSlN6tFBH3Aa8BkHQa8Dtlx2Dt0ZSJUbOylZ4cJD09IvZKOgL4AN3OJbMk2txtYrYUSctK\nkq4Dvg48W9IuSRcDF0i6H7gP2A38S8oYzMxsdElHDhFxQZ+bPpnydc3MbGm8fIaZWcVy7KpzcjAz\nq1iOXXVeW8nMrGI5dtU5OZiZVSzHrjqXlczMrIeTg5mZ9XBysAXl2D1hZuVxcrAF5dY94WRlVi4n\nB1vQ+rVTbFi3MpvuiVGTlZNJfvyepJNi27pbyRaUW/fEqK1+Tfmayibxe5JOim3r5GC1MGqyyrFv\nPDe79z/GbTv3sX7tVCnf2Ob3JJ0U21YRMbEnS6nT6cT09HTVYVgiZe+oDK7/1ve5efseNqxb6SP5\nBpM0ExGdUR/nkYNlwSWH8vlI3hbj5GBZ8I6qfLnNK1lenBwsC95RmeXFraxmZtbDycGexL3oZgZO\nDjZPbmdGm1k1ks45SFoNXAusAALYFBGflHQicD2wBngAeFNE/CRlLDYcTwybGaQfOTwBXBYRpwMv\nBt4l6XTgcuArEbEW+Erxs2VgdmLY5xqYtVvS5BAReyJia3H5ILADWAWcD1xT3O0a4HUp4zAzs9GU\nNucgaQ1wJnAHsCIi9hQ3/ZBu2Wmhx2yUNC1pet8+18DNzMpSSnKQdDxwI3BpRDw697bort+x4Boe\nEbEpIjoR0Zmacg3czKwsyZODpKPoJobNEXFTcfXDklYWt68E9qaOw8zMhpc0OUgScBWwIyKunHPT\n54GLissXAZ9LGYeZWROlPC8p9cjhHOBC4JWSthX/zgX+CvgtSTuBVxc/m5nZCFKel5T0PIeIuB1Q\nn5tflfK1zcyaLuV5ST5D2ibCy26YlS/leUlODjYRZSy74QRkVh4v2W0TUcayG/5CILPyODnYRJTx\nfQxe98nqrk5fh+uyktWG132yUeRYhqzTqsdODtZIOe4YrFw57ojXr51iw7qVExv9pvw7d1nJGsnz\nE5ZjGXLc8mu/clTKv3MnB2ukHHcMVq4mfS95vySQ8u/cycEap06TfmbD6JcEUiZAzzlY4+RYazZb\niiqaMTxysMZxScls6ZwcrHGaVGs2q4rLStZabnc168/JITPeYZXHcxNm/bmslBn355dnlLkJd0BZ\n2zg5ZMaTqeUZZW7CSdtmteVAwckhM55MzdOkknZbdixN1pYDBScHsyFMKmm3ZcfSZG0Z3SedkJa0\nWtJXJd0r6R5JlxTXv7H4+bCkTsoYzHIy6YXX6qjuTRdtWR049cjhCeCyiNgqaRkwI2kLcDfwBuAf\nE7++WVZcNvToqS4GJgdJRwCXR8RfjvrkEbEH2FNcPihpB7AqIrYUzz3qU5pZzbWlLFN3A8tKEXEY\nOG+pLyRpDXAmcMcIj9koaVrS9L597kU3a4K2lGXqbtg5hzslfagYRYxM0vHAjcClEfHosI+LiE0R\n0YmIztSUjzJsNKPWtuteCzebpGF39icCbwZ2S/qcpI9IeuMwD5R0FN3EsDkibhozTstEnXago54B\n7TOmzX5pqAnpiHgTgKSnAM8Dng+cDdyw2OPUnVS4CtgREVcuLVTLQZ0mE0etbbsWbvZLioh0Ty69\nFLgNuAs4XFx9BfAU4O+AKWA/sC0ifnux5+p0OjE9PZ0sVhuOT+IyqxdJMxEx8ikDSVtZI+J2oF9L\n0mdTvral4VZMs3ykPFjzqqxmZjWVcp7My2eY9eESmuUu5TyZRw7WOJPqqHL3kuUu5TkjHjlY40yq\no8rdS9ZmTg7WOJPaqXvy3drMZSWrtYVKSF6ewWzpnBys1jwvMHl1Ogt+FE39vVJxWclqzfMCk1en\ns+BH0dTfKxUnB6s1zwtMXlMTblN/r1SSLp8xSV4+w8xsdOMun+E5B2uFSdWbXbe2tnBysFaY1MS1\nJ8CtLTznYK0wqXrzJJ7Hy3JYHXjkYCOrY2llUuc+TOJ5PPqwOvDIwUbmlsClcdeM1YGTg43MO7el\nqXv7rcti7eCyko2syctT1LFkVjaXxdrBIwezOVwy+6V+IwSPHNsh6chB0mpJX5V0r6R7JF1SXP83\nku6TdKekz0panjIOK1edj77Xr51iw7qVQ+/46vy7DtJvhFDVyLHJ2zpHqctKTwCXRcTpwIuBd0k6\nHdgCnBERLwDuB/4scRxWojqXHUbd8dX5dx1k1ESZWpO3dY6SlpUiYg+wp7h8UNIOYFVE/Pecu30D\n+L2UcVi52lR2aPLvmtvEed23dd0m8ktbW0nSGuBWuiOGR+dcfzNwfUT82wKP2QhsBDj11FPPevDB\nB0uJ1Sw3dduxWK/rv/V9bt6+hw3rVpaadLNeW0nS8cCNwKXzEsP76ZaeNi/0uIjYFBGdiOhMTdXz\naMFsEiZRUnHNvlq5lekGSd6tJOkouolhc0TcNOf6twHnAa+KuiwNa1aRSZRU3IlVrdzKdIMkTQ6S\nBFwF7IiIK+dc/1rgfcDLI+JQyhjMmmASO5a61+ytXKlHDucAFwJ3SdpWXHcF8CngKcCWbv7gGxHx\njsSxmLVaFUeuniupr9TdSrcDWuCmW1K+rpnlwaWs+vIZ0maWTFWlLI9Yls5rK5lZMlWdTe0T5pbO\nIwczaxxPvi+dk4OZNU7d2kZz5LKSWU35pDZLycnBrKZcV7eUXFYyGyDXzhfX1S0ljxzMBsj1CL3J\n38jXRrmVCT1yMBvAR+hWhtxOGHRyMBvAnS82aQuVKnM7CHFZyawGcis51EWu222hUmWKMuHu/Y9x\nxLEnPG2cxzo52MRU9UHMZQeQMo5c5z1yl+t2K+u7HW7buY8jnrrsxHEe67KSTUxVNdNcarVz41i/\ndmqiHU65lRwg3y6uuXLcblBeqXL92ikO/+zgI+M81snBJqaqD2IuO4C5cUw6YeU47/GF7bu5YWYX\nBw49zsaXP2vkx5eRXHLcbmU6ZfkxHD504MfjPNbJwSamqg9iLjuAuXGMk7DqcCQ+VxBEBMF4X+RY\n9Yivbtu7bE4OZgmMk7Cq3lmOasO6VSw/9uixR2xVj/jqtr3L5uRgjVW3I8Oqd5ajWuqIreoRX922\nd9ncrWSNlWunSj/jtDLm0qlVRz7DfHFJRw6SVgPXAiuAADZFxCclfQQ4HzgM7AXeFhG7U8Zi7dOG\nI0OXRiwVRYw3mTTUk0srgZURsVXSMmAGeB2wKyIeLe7zp8DpEfGOxZ6r0+nE9PR0slgtjbqVdurG\n29cGkTQTEZ1RH5e0rBQReyJia3H5ILADWDWbGArHwZjtDpa9upV2UkpRAkp1Vq1LVVbanIOkNcCZ\nwB3Fzx+V9APgLcAH+zxmo6RpSdP79nnnUkdlnQlapoV2nsPsUOuSKOsSJziRpVRKcpB0PHAjcOns\nqCEi3h8Rq4HNwLsXelxEbIqITkR0pqaas3NpkyZO+i208xxmh7pYosxhJzcbw2krltUmoeeWyHJ4\nHycleSurpKPoJobNEXHTAnfZDNwCfCh1LGaTsNBE9/q1U+w/9AsOHHqc3fsfWzAZLta6mcPE8mwM\nG9aVF8NS50xyazrI4X2clNTdSgKuAnZExJVzrl8bETuLH88H7ksZhzVPlROxC+3kT1l+DMuPPZqb\nt+/hhGOPGnnHkMNOrooYxtmZzn/vc9oJ5/A+TkrqkcM5wIXAXZK2FdddAVws6dl0W1kfBBbtVDKb\na/f+x7hyy7d58EeHgHyO0JayY8hhJzcbw2xppIzEO842y/noPIf3cVKSJoeIuB3QAjfdkvJ1rdlu\n27mPB398iF8/6djSjtCGGak0ZcdQ5s53nG02akJxu+94vHyG1c7cnUNZH/acj1YnLffSyKgJpU3v\n3SQ5OViWFjvaq+II/bQVyzj5hEc4bcWyUl+3Ck0ZAc3KPdnlymsrWZZya1G8/+GD/PDAz7n/4YNV\nh2IjamI79XwpWmg9crAs5Xa0Nyge17WtSilKZ04OlqXcShuD4nFd26qU4mDKycGWxEfMXbmNdOrG\nf0dLk+JgynMOtiS5zQ1UpQ117ZT8d5QfjxzsSUY9gvMRs496J8F/R/nxyMGeZNQjOB8x+6h3Evx3\nlB+PHOxJfAQ3uhTbbPf+x/jC9t0EwYZ1q7zTrLk6ji6dHOxJcusSqoMU2+y2nfu4YWYXEcHyY4/2\ne0I9d7Cz6tjN5uRgQL0/eE20fu0UBw49ThAexRXquIOdVccRuZODAfX+4DXRKcuPYePLn1V1GFmp\n4w52Vh1H5E4OBtT7g2d5SD36rOMOts7crWSAu0Vs6dy11SweOZjZRHj02SxODmY2ES77NIvLSmZW\nqRTLTdvSJU0OklZL+qqkeyXdI+mSebdfJikknZQyDjPLl+cq8pS6rPQEcFlEbJW0DJiRtCUi7pW0\nGngN8P3EMZhZxjxXkaekI4eI2BMRW4vLB4EdwKri5k8A7wMiZQxmljd3yuWptDkHSWuAM4E7JJ0P\nPBQR28t6fTMzG14p3UqSjgduBC6lW2q6gm5JadDjNgIbAU491V0QZmZlST5ykHQU3cSwOSJuAp4F\nPBPYLukB4BnAVkknz39sRGyKiE5EdKamXI+sI3eimNVT6m4lAVcBOyLiSoCIuCsinh4RayJiDbAL\neFFE/DBlLG2S0w7ZnShm9ZR65HAOcCHwSknbin/nJn7N1stph7x+7RQb1q10J4otKKcDGXuypHMO\nEXE7oAH3WZMyhjbKqTUwh7NmvRx5vrwacL68fEYD5bBDzol3QPnK6UDGnszJwRqvjTuguoyWfCCT\nr1qvrVTnemWdY6+bNp5kldO8k9VTrZPDuB+AHHbM/vBaSm4EsKWqdVlp3HJBDjXoNpY6mqIOJRuX\na2ypap0cxv0A5LBj9oe3vr6wfTc3zOziwKHH/T3P1li1LiuNq401aJucIIgIosFrRuZQerVq1Xrk\nYDZXWeWeDetWsfzYo7MoCab6nXMovVq1WjlysGYqa5I/p5HnoN95sRHAYrd5Qts8chhTHSYl2yaH\nuaSyDfqdFxsBLHab58TMyWFMHnbnp207tGEOUBZLHm1MpjY8J4cx+YNlVRvmAGWxhNm2ZGqjcXIY\nkz9YVjUfoFhKnpA2K9mk2kRzmhi35nFysIlzj/ziRumq8ra0qrisZBPnyfrFjVIO8ra0qjg52MS5\nFr64UearvC2tKoqoxxIAnU4npqenqw7DzKxWJM1ERGfUx3nOwZbMdXGz5kmaHCStlvRVSfdKukfS\nJcX1fy7pIUnbin/npozD0vJ3U5g1T+o5hyeAyyJiq6RlwIykLcVtn4iIjyd+fSuB6+JmzZM0OUTE\nHmBPcfmgpB3AqpSvaWZmS1dat5KkNcCZwB3AOcCfSHorME13dPGTBR6zEdhY/PhzSXeXE20SJwE/\nqjqIJegb/xHHnvC0I5667MTDPzv4yOFDB35cclzDaOy2rwnHX61nj/OgUrqVJB0P/A/w0Yi4SdIK\nuhs7gI8AKyPi7QOeY3qcGfdcOP7q1Dl2cPxVa2v8ybuVJB0F3AhsjoibACLi4Yj4v4g4DPwTcHbq\nOMzMbHipu5UEXAXsiIgr51y/cs7dXg/UuVxkZtY4qecczgEuBO6StK247grgAkkvpFtWegD4oyGe\na1OSCMvj+KtT59jB8VetlfHX5gxpMzMrj8+QNjOzHk4OZmbWI9vkIOmNxZIbhyX1bcOS9FpJ35b0\nHUmXlxnjYiSdKGmLpJ3F/7/W534PSLqrWEak0pUFB21LdX2quP1OSS+qIs5+hoj/FZIOzFm25YNV\nxLkQSVdL2tvvXJ4abPtB8We77aH/Uj/z7pPtezBk/KO9BxGR5T/guXRP3vga0Olzn18Bvgv8BnA0\nsB04verYi9g+BlxeXL4c+Os+93sAOCmDeAduS+Bc4IuAgBcDd1Qd94jxvwL4QtWx9on/ZcCLgLv7\n3J7tth8y/my3fRHfSuBFxeVlwP01+/sfJv6R3oNsRw4RsSMivj3gbmcD34mI70XEL4B/B85PH91Q\nzgeuKS5fA7yuwliGMcy2PB+4Nrq+ASyf15ZcpZz/FgaKiFuBRxa5S87bfpj4sxYReyJia3H5ILDQ\nUj/ZvgdDxj+SbJPDkFYBP5jz8y7yWbtpRXTXlgL4IbCiz/0C+LKkmWK5kKoMsy1z3t7DxvaSoiTw\nRUnPKye0ich52w+rFtt+3lI/c9XiPVgkfhjhPaj0m+AkfRk4eYGb3h8Rnys7nlEtFv/cHyIiJPXr\nGX5pRDwk6enAFkn3FUdhNnlbgVMj4qfFMvH/AaytOKa2qMW2L5b6uRG4NCIerTqeUQ2If6T3oNLk\nEBGvXuJTPASsnvPzM4rrSrFY/JIelrQyIvYUQ8+9fZ7joeL/vZI+S7c8UkVyGGZbVrq9BxgY29wP\nS0TcIukfJJ0UEXVYVC3nbT9QHbb9Qkv9zJP1ezAo/lHfg7qXlb4FrJX0TElHA28GPl9xTLM+D1xU\nXL4I6BkJSTpO3e+5QNJxwGuobimRYbbl54G3Fl0bLwYOzCmdVW1g/JJOlqTi8tl0//5zXEV2ITlv\n+4Fy3/ZFbD1L/cyT7XswTPwjvwdVz7IvMvv+ero1vZ8DDwNfKq4/Bbhlzv3OpTsz/1265ajKYy/i\nehrwFWAn8GXgxPnx0+2s2V78u6fq+BfalsA7gHcUlwX8fXH7XfTpIss4/ncX23k78A3gJVXHPCf2\n6+h+98njxd/9xTXb9oPiz3bbF/G9lO78353AtuLfuXV5D4aMf6T3wMtnmJlZj7qXlczMLAEnBzMz\n6+HkYGZmPZwczMysh5ODmZn1cHIwM7MeTg5mZtbDycGsBJJeIunDVcdhNiyfBGdmZj08cjBbhKTn\nSrq1WOb4vZK+M+bz3CBp/aTjM0vFycGsD0lHApuBSyLiBXTXwhp3YcQz6K57Y1YLlS7ZbZa5NwDb\nI+J/i5/vZd7S68N8J4mkpwJHR8SBlMGaTZKTg1l/L6C7uuWsM4D/mnuHGO47SZ5HN7GY1YbLSmb9\n/Rg4DUDSC4E/oLvc8aiej0tKVjMeOZj196/Af0q6C/ga8EBEfG+M53k+8M1JBmaWmltZzfqQdHxE\n/LS4/F7ghIj4QMVhmZXCZSWz/t4j6R5J24A1wEcqjsesNB45mJlZD48czMysh5ODmZn1cHIwM7Me\nTg5mZtbDycHMzHo4OZiZWQ8nBzMz6/H/qkNnysUwXj0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x107565208>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "mag = cmd_mags['r']\n",
    "color = cmd_mags['g'] - cmd_mags['i']\n",
    "\n",
    "msk = (tab['phot_sg']=='STAR')&(tab['sep'] < central_radius) #arcmi\n",
    "plt.scatter(color[msk], mag[msk], s=2, alpha=.5)\n",
    "\n",
    "plt.xlim(-1, 2.5)\n",
    "plt.ylim(25, 17)\n",
    "plt.xlabel('$g-i$')\n",
    "plt.ylabel('$r$')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.text.Text at 0x11772cc50>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEOCAYAAABiodtuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzsvX10m9d5J/i7FAkSIEHwU6AoihStkE4oJ7RjSo6TpZyk\ndpxzIp3stLHzodNtZqZR0tndbjed7mRaN3HSzOx0J/HZtHP2JMq0nW1HzknktNvK8cRR0kZiMrYs\nyjajiIoIyRQ/RBIEQQIgCZAAybt/gPflfS/u+4kXIEi+v3N8aAF473s/n+e5zyehlMKFCxcuXLjg\nUbbdHXDhwoULF6UHlzm4cOHChYscuMzBhQsXLlzkwGUOLly4cOEiBy5zcOHChQsXOXCZgwsXLly4\nyEFBmQMh5C8JIbOEkF9yn32XEPLm5n93CSFvFrIPLly4cOHCOkgh4xwIIScALAH4a0rpA5Lvvw4g\nTin9SsE64cKFCxcuLKO8kI1TSi8TQg7LviOEEABPA/hgIfvgwoULFy6so6DMwQD9AMKU0pDWDwgh\nZwCcAYDq6uqH3/72txerby5KBJn1DSytrKGmqhwV+3K1oAvLacRSGdR5K1Bf7XGkzUI/78KFCDv7\n2CyuXbs2RylttvrcdjKHTwL4jt4PKKVnAZwFgL6+Pjo4OFiMfrnYQZiKpTAQiqC/qxmtdV5Tz3z3\n6jguDE3jVO8BfPxYu+V3fuvSbZwfnMRTfW347GNvs/z8boGduXchh9Fc5jPXhJAxO33aFuZACCkH\n8OsAHt6O97vYOTA6FK11XssEvr+rWfXXKggICCEgILae3y0YCEVwYWgaAGwxWRdbMNrH2zHX23Vz\neBzAryilk9v0fhc7BIU4FHYYCo+Tva0I+CpsM5fdgnyZLA/3FqIPJ+faLArtyvodAK8AuJ8QMkkI\n+ZebX30CBiolFy6A7GE41XvA8FBMxVL47tVxTMVSBe8TYy6lRsSKOQeAs/PAhICBUMSBnhUWxZ5n\nO4zTiT4W2lvpkxqff7qQ73Wxe2BWyndVHPpzUOqS+XZIxnZR7L1m5338M3axnQZpFy4cg0hcSp0Y\nOgFxjHoEdjuYp5U1yFfVV0wUm5HZeZ8TfSxoEJyTcL2VXFgB80jq72pEnc+D7qAfI+FFQ0K1k5iK\nFa+r7RhXvl5hLpwBIeQapbTP6nPuzcHFrgSTmOLJDC4MTaMlMI+Z+CoAfcl5J6mnrEiH2yGZ7yRV\nkYtcuMzBxY6CWQmYEcOpWAoBX4Xq5qAHLYL2xvgCzl0Zw5M9LZhPpnPePxVL4cLQPRAQnOxtBYCC\n+a2LYyxVFLJ/O+WGt1P6KYPLHFzsKFiV7HkC9VB7vaXf8zh3ZQwXh8MYnVuGt6I85/0DoQjOD06C\nEIKArwIAdPtpdRw7mcgUAjvlhrdT+imDyxxc7Chsl6ri9CMdAKC6OYj9iiXTICCq77T6aXUcO5nI\nFAI7RWW1U/opg2uQduFiB8C9ObiwC9cg7cLFLkap2xdc7D64KSVduHDhwkUOXObgwgWKnxLBhYtS\nh8scXOwpiEyA/fvC0L2C5fbZyYxnJ/e9kNgL8+IyBxd7CmKCN/ZvAmIqwZ8T79xJ2Ml9LyT2wry4\nBmkXewqiayH/t1BeQDvZnXEn972QMDsvO9nLzL05uChpmL2+m/2dmGY6n7TTdt+5k2Cn73tB5WIW\nO/mG4TIHFyUNs4dL63eFJFT8O12CuIWdQBDzXS+zYzRbj6QU4aqVXGwbzFy5zV7ftX5nJbLYqgqA\nf6cbwbyFfFRRxVLD5Lte/V3NiCcziCXTmIqlNPu6k+NTXObgYttg5oCaPVxav7NCqPLJ22SWWOwF\n5EMQi8Vk87WltNZ5EfBV4MLQNOp8nh3LAPTgMgcX24ZiGPWsECqrBEPs124nFsVAsQzgTkj0Tve1\n1IzXha4h/ZeEkFlCyC+5zx4khLxKCHmTEDJICDleyD64KF2YNXYWS4ct9sdILy32q5j65d1q49hu\n472VeXW6r6Vmqym0Qfq/APiw8Nn/BeDLlNIHAXxx898uXGjCCaJrh5gaHVaxX4UgbFr9LjVCsltg\ndl4LwZyt7PNiCAcFVStRSi8TQg6LHwOo3fz/AICpQvbBxc6HEyoAO7psI7VBoY2NU7EUnrt4C2Nz\nSQDqfrvxB4WB2Xk12k92VERW9lMxbDPbYXP4PQAvE0K+huzN5b3b0AcXewx2iGm+xJ8RiAafBy8P\nz+D0Ix2mCg4xDIQiGIsm0dHky+n3TvaCKWXw86pF4KdiKcSSaZzoatLcT4Um3sUQDraDOfwOgP+d\nUvp9QsjTAP4CwOOyHxJCzgA4AwDt7e5BcGEf20FMGYFIZdYQCi9hdG4Zf/7Jd5uWJIsRvV1sWJGo\nnTbQWm1Pi8APhCIYCEVxqvdA3i7YdlGM/bwdQXC/BeBvN///PABNgzSl9CyltI9S2tfc7F6fXZQe\n9HS/TIf8uRNH0BWsQWadWrIRbLdxthCwYitx2q5i1Z7QHfRLbQBmbAO7Ye224+YwBeAxAD8F8EEA\noW3og4tdjlIIpuKlu6MHA0p/9jKsSNROS99W7QmneuUqob2i0itomVBCyHcAvB9AE4AwgC8BuAXg\nG8gyphUA/4pSes2oLbdMqAseMuLPf3Zh6B7OD07iqb42fPaxt+k+Z+dd+bS3W7Fb5qIUx5FPn0qy\nTCil9JMaXz1cyPe62P2QSez8ZwQEhBAQEMPn7LyLYa9IkWZQSCNsMQl2Ka7pQCiCF65N4urdeXz+\nifuLwrTcCGkXOxIyFYH4WcBXIdUXi89pgRGkrO7ZdRs1QiGNsHs9d1V/VzOu3p3H2FwSA6FIUeag\noGolJ+GqlVwUG9kKcdM41XtgTxKkQmM7PZd2IuzOgV21kpuy24ULDWh5pey21BXbNR4r3ki7wfvH\nLtj6ACjqHLjMwcWehBFB1JPSSjV1hV0iv10pI2TMd7cxXiewXfvNtTm42JMw0mHrfb9dqSuM1Ap2\n9fJOpIxwKl3EXrctyGBmfQqhdnOZg4sdAb1UBnYOhdGB0/veKW8WpyJ2GewyLbPj0WvfKaLu5ozK\nhdH66OXgygcuc3CxI6CXysAqUTJDlIvhzmi179udCFCvfaeIejGSGe42w7ZeDq584DIHFyUFrcOr\nRXzsEKVSUV2wPncH/fju1XFDglWK/vcMZvu23cS5VNbeSRQqB5drkHZRUtAyvml5q+h5sWgZN7uD\nfrQEKtEd9Ds/AAtgfR8JL5akgbsQ2G5jPm8Ed8L4vR0GdPGdhfLkcm8OLkoKTuqctaTEkfAiZuKr\nGAkvWkqhbRdG0vJe0rOXgjGf7QUWxwLYv0W8ODSF89cmEU9mcOaxI471Vw/Fuv24zMFF0aFHLJ1U\nnTipihJhRT2idZhlBGu3w8z6FkL1JFsDJ/YBBQWlFBTFCyYuFoN11Uouio7tqglt9DmDGVWBbAxa\nz2kF0zk1D7stNsDqvJgZv2wNrNYMl+F4ZyMebK/D8c5G08/YfVex4d4cXBQdxVYt5OMy2t/VbNpA\nrnVD0JKWnZqH3WZk5efFzNqZGb+ZG4udebSrosxnzVy1kotdi2J73eTjMmqF4Fsl9k7Nw26zWfDz\nomUX4JmGU+O3047dd9t9biqWQjyZQX9Xo+mgOLtwE++52BHIpw5Dd9CPkfCiLR32drte7nVozf9e\nTYpodtz87z5xvKP06jm4cOEU8qnDoFXRywxKObZgL6C1zitV7WlJ3tvFzIv1XrM3DiduUy5zcFE0\n5HOAnLryuzeB0kC+3l5aTNuMZ5jT616o9BUymBVWnBBqXG8lF0VDPt45dgJ9ZM9sdxBWKaOYHjRW\n1kHL28vKbwu57mbTV9iZX0fWZF9FhZ3H3JuDi6KhmIZT0SCXj/Fyr9w2iun1ZGUdtKRg2boU2jNM\nBrPpK5wuUWsW+7z+BjvPFZQ5EEL+EsBJALOU0gc2P+sF8E0ANQDuAjhNKU0Ush8uSgPF1N/zhyqe\nzOD518Zx6VYEz5zswcePtSsSmSzL64WheyAgONnbitY6Ly4M3cP5wUnEkml89rG3FaX/24FiMm8n\n9oIVwlnIvSdrW8a4iukN5YS3UqHVSv8FwIeFz/4zgC9QSt8J4O8A/EGB++BiD4ER/Wzd56x6gYJi\neXUNw9MJRa2gpWYYCEVwfnAS569NKt8REBBCQECKPp5iYrurrVlVoVhRNxUbsv3llGrU6vvXU4vz\nlh7eREFvDpTSy4SQw8LH3QAub/7/RQAvA/jjQvbDxe7AVCyFF4emQEFxqvcgWuu8ORIak/Kf6mtT\npPxTvQcBZIm8KImJxupYMo0PH21BrbdC+e5kbysCvoqSJEK7CVZvaNuVisMMtjv2RPX+9UzGThvb\nYXO4AeCjAP4/AE8BOKT1Q0LIGQBnAKC93XUn3OsYCEVw/tokKKWo83nw8WPtOaqFRGoNiZU1JFJr\nynOtdV6F2PDqJFkVsoFQNMeH3HVnLQ4KcUPbrujxYuyZQuco2w7m8C8A/Bkh5I8B/AOAtNYPKaVn\nAZwFskFwxemei1JFf1cz4skMKKjmDaDWW47aqnLUenO3tpHL4XZLe3sdVm5oZm8Eu3lNRcbn9C2p\n6MyBUvorAB8CAEJIN4CPFLsPLkoPZquznextlepxGU71HkSdzyMlBkYuh6V4QyiWWqQUPLKszL/Z\nG4GZEpvbPW4rfdFLG+L0LanozIEQsp9SOksIKQPwDLKeSy72OMxsbDPBRmZLWRYiEMpJIsPaiycz\nuByaA1BYtchOS963G5MWmumL+Bv+d07fkgrtyvodAO8H0EQImQTwJQA1hJD/efMnfwvgrwrZBxc7\nA2Y2NpP8m/0exJMZTMVStjw/CoGBUAQvXJvE1bvz+PwT9+fNIJhx9smjwaJ45BRL/eIUE91tSQuZ\nM0RvWwCxZFpzb+v11+n9XWhvpU9qfPWNQr7Xxc6DmY3NDkQsmcbl0BwCvgrTh4ERpQafBy8Pz+D0\nIx2qFMv5Eq3+rmZcvTuPsbkkBkKRvA8pM84GvJ6iSLTFUqmVkqQO6I+70Convn3mDNESqMTQ5Kri\ncGGlv07DjZB2sWPADsZULKVpV9ACI0qpzBpuzy4BgIo55Eu0Wuu8+PwT9+cdeMSwW91nrUrq22kT\nsLsnzPZZrBsCQJVBmG9PDMwsBlzm4GLHQStTp1bKDPZ7AKqbAw8n1AtmI2XttrUbYHVc23nT0NsT\neutqts+iDYz9ViwaxAIzCSGmbstOMVSXObhwHMWQ9mQHkP8MgOp7/vA9cbQlp71CEWMzVeVcaGM7\nbQJ6e0KPAZjts9k919/VjPHoMkbnkugO+g1/7xRDdZmDC8dRDGlPdgD1PtNDIZkZ36dS07fvBBSC\naTux3t1BP1oC81Ji7XSfW+u8aG+sxtBkwlQ5UqcYqsscXDiOfDan2YMrO4DiZ8XMfKkFvk+l4hlj\nhFLy/S8EnFhvvdrRhXBrjiXTONHVlGOLkL3HKebkMgcXjiOfzZlvsRY9u4MWikW0d4odoRRvOE4S\nXCfWW68Np+ePT+siq03i1HtEuMzBRUlB69DpHQTRJVDL7iD7fWud1zTRLoZEXQpSuxHx3I4+OkkI\nnWDSZoMtnYBWe4UWalzm4KKkYKdYi8wlUM/uYJfQyJ4TCWW+hLMUpHaj4jrFitrmsVNUcoB6/pxg\npFrrUeibqMscXOwIWEmLYWR34H3KZQV/tCB7TiTmVoi7UwVhjOCUpM/G1t/VWPQ6CnyMi5U12044\nXVvajmrVdWV1sWchOwhG0jwjNN+9Oo4LQ9OIJzNKwJkZIzh7Dsgl5mYZjxbhKIQ06NRtRGTC+cIO\nESvkzcrpYDM+0aNVQUSrPTNjd11ZXexZaNkYZPEOetI8n47DymHq78qmDo8l0znPyBiIrE2jDLFO\nwqnbiNOMyw4RK6R6yWqwmRGcdmPm29NjrK4rq4s9B3YgYsk0BkJRAObiHbQOi910HK11XgR8Fbgw\nNK2ZA8fogDotheuhVLykRIJmh4gVciz9Xc2IJdOqioFGMFtwx260tVZ7esKHmEHALgilO6OGTl9f\nHx0cHNzubrgoArQOCzsQJ7qaTKmBtqOPLrTB1k+stLeT4cSY7LRhtP/4Nj9xvOMapbTPar/cm4OL\nkoORGqgUCHKpSONOwgmGp9eGXlTxTkWhYya0YLT/nOiXyxxclByM1EAyiMFvLw5NgYLiVO/BHCLl\nSv1yOKEX12tDL6p4p6LQMRNmINvPqjb3VVTYaddlDi5KDnYOixj8dv7aJCilUpvAdsUSlDpTKrQU\nvJNiFXYSjPZzWaWv1k67LnNwUTLQc0EF9FNhMGNiPJlBZ1M1uvbXoFPwBGLtdQf9ONVb/FiCUghw\n00OhpeDdqIorBRgx3Y3VZMJOu2X2u+TCRS5YkNJULGX5WUY8B0KRnH+L38kQml3CxZthvDw8g8TK\nGtobq6W5aEbCi0oab6fGZaZ//V3NRQ8e220ws79kv8lnXzrRt3zfn9fz65mMnXcWlDkQQg4RQv6J\nEDJMCLlBCPnfNj9vIIRcJISENv/uDgWkC1NEUgs88ZyKpTARTaK2ah+6g/5NY2YluoN+6UEZCEUw\nNpdER6MPpx/pkBLhfIiz0bjMtM0kZ5kNpNCEa7fAzP6S/Ub8zMk5Z229ODSl2Tez50KrX3rP53Pm\n9FBotdIagN+nlL5OCPEDuEYIuQjg0wB+Qin9D4SQLwD4AoB/U+C+uCgC8tEri37cP7wxA0IIRsKL\nAKAYM0fCiyr1DEtp/ERPUIlslRk881FrGI2rEJlojWAmOny3wcz+MhP74qSKbyAUwQuDk2j2e9Df\n1ZiXvcWMp57dtq2iqHEOhJC/B/CfNv97P6V0mhByAMBPKaX36z3rxjnsLfCpDI51NuDq6LzifQSo\n7Q8yP/GdRCTt9lU27t0US1DINdRr2+p7lVQo0SQ+9nCb455H+YIQUtpxDoSQwwAeAnAFQJBSylxL\nZgAENZ45A+AMALS37+yN7kIbWq54n33sbQCyBO9yaE6Vz95Iqip14y8Pu7cOcdxTsRTiyYym9FpM\nOEHkilWEyep7Zbm6Pv/E/XlHJBv1q9goikGaEFID4PsAfo9SqrKc0+zVRXp9oZSepZT2UUr7mptd\nI95uRb76fD5dANPVlpLxt1DGSNGGMRCK4HJoDnU+j+LttR22DCZJv3BtMi89uJNraGUujN4r26/8\nWmi9q1DrYfS+ko1zIIRUIMsYzlFK/3bz4zAh5ACnVpotdD9clC7s6vP1EvAVUwIzkpLzlYDNPs8S\nAo5Hl3H20h0AKHrdBYBzDtBwJTZ7m3ByDcWaH3r9yDf6WGu9CnUTMnpfScY5EEIIgL8AcJNS+hz3\n1T8A+C0A/2Hz798Xsh8uSht2iQDb/FvJ0nLVKcWwPRgd+nwNhlrPy9QbAV8Fvjc4AUIIPvbwQUck\nb6tzqJXmZDtVfXyf8u2HXeZhZR9YmXOj952zGedQUIM0IeR/ADAA4DqAjc2P/xBZu8P3ALQDGAPw\nNKV0Xq8t1yDtQsRULIUXh6bw5sQCIktpqTGwEMn6nK7+ZhdahnirNQmsJHErlLE1nzm0Y0CWeXmZ\nmbdCr7WYedjMnPPBnSPhxZy+laRBmlL6MwBE4+tfK+S7Xex+MEk5spjWrIvAPrNas0EP26nC4iGT\nGHlDvlkU+ubD982uEVgGu2VLZf0wW8uh0Lcf1v6JribVrU+PKbFnWgLzGJtL4urdeXz+ifvzZl5u\n+gwXJYs3xhdw7soYTj/SgYfa65WbAp9QT6bCkElSAEzVbDAjGRYrR5BRX5xiSnrjKdatyM6cMqLY\n39WIE11NiCXTmIqlbPWTpV8xquVQ6LVnmWuPdTaoYnX0mBLrS3fQj3NXxjAWTWIgFMl7b7jMwUXJ\n4tyVMVwcDgMAHmqvz0p3QkI9sbBJa51XJUnNxFcBZA+UmcNiRjIs1k0hH2naCjE3K9EbGXKd7oNR\nO1p2BCYE8M+YYbTsxsXXqWZzwNt17Ky9KOhoQStzrR5T4vsUrK1yxKUWcJmDC4fhpKR5+pEO1V/m\njUNB0d/VnKOfBbJElJek+JuDGdiVDAshYecjTQPOqD1EAvzC4KSh2sKpPhi1I6u0Fk9m8vYU4n8L\nwJGxnLsyhh/dCGN0bhl//sl3m7qVintKy2NPvE07Jbi4zMGFo3CSOD3UXq+SnlrrvDjz2BHl37yx\nWVQrmJVyReTrOQU4p4u205d81R4yDyieAF+9O69SW8iYYj594Nuz0g7r51QspTge8DBqy+i9+Uri\npx/pwOjcMjJrVOV+Le5Pfr7PXrqD89cmEU9mVPueB3+bBtS3JnZbIR6vz06fXeagg52UgqFUkK/k\nreVxofWuWDINAKCgyu2hzufRNVIarauddS+WHcII+UqOekxOFgks+72TeaastqP1bqM+ieozs88x\nGO2Zh9rr8eeffLdmXI4MFBSUUlB5jLDSV3abJiCqNpladp8vYGtT7lnmsBvy75ci8pW8RTuBDPza\n1fk8Ku+OeDKDc6+OY3VtHa11XmlJSiM9up115yVXpq/eiQKFEZMT19dpprhdTDbfOAgrtqqpWArj\n0SRqvft0S6ae6j1o6ETB36bFWxNTx15Pxm2Fqe9Z5mBmMUtFGtyNEJmzkZ2Av1mcuzKGsbkkgNyA\nq6lYCi8Pz2AkvIyFZEZaklKPELAMr71tAV3vFy3hYicJFDLffqvM3Wnj/Ha5BcvsFzK9vxas0IqB\nUAQ/Gg6DUqpbMtVMZgAxzkYUUL721IP4+tOppGGnJNizzMHMYpZSEqzdBr1YAdlhUflyR7dSM4hr\n1FrnxTMf6cG3B+6gs7EG3UF/jiSvRQjYewZCUbQEKjE0uapIbmZvFztJoDDr278TYVclLD7H7FqA\nPrO3on6KJdP4UE8Qtd5yw30iY+Di3jP6t13sWebgEv7thVUiKrtZiHEN7Dcj4UU885GjaK3z4luX\nbuP84CRiybQ0OExLTcK/R3bYtPpvdV9tp12L2WyMfPu3G2bnSMy19cI1Y88qEeJaO83smfBhNtpc\nxsDFPjG7A7vpOtXnPcscXGwvrBr6uoN+xJMZvDYaVVz2ZBGygNr1kICAEAKiGaiv3S92g5EdNjse\nUTKju579w46R3grsRFNvB8xKwuJcXr2bjRi2EhAmrrXTtiQ7QpHIwGW35YCvAheGppX4HycEX5c5\nuLANJ6VeIx1+S2AeQxNxVQCcLEL2eGcjTvVuHb6Tva1S10Yr0GJk7P3xZMZU3iaZ0V20f/BxBFaM\n9KWOfPaKWYIq2p/s1FgwWmvAuvebmfb1fm+GgRdCnekyBxe24aTxVdYWX7zmeGcjuvfPKwFwgNyw\nzBgHg9PqQ5k/vJi3SYtYaKnGtOIIjIz0Own57BWzayiTqJ1aez3iWwpOCIVQk7vMoYSw0+IqnJRW\nWE4Z3rVvIBRRKsA91F6vpAZgEAkrr3dlz+czl7L1kPnhT8VSKpdDLWJhZHQXpV2j3+8k2NkrpXIe\nzKTxiCXTiCczinebEwbxcGLFMOVGIefIZQ4lhFKQQKzAqQAhQJ5TRmQYWukbWPssEI4FwRlFlxrB\nrCHaSd//Qt50zBAPszmArMLOuErlPJhJ48HibZjR2G7feUN6cnUdP7+TtaVprYXWjdsJhuEyhxLC\nTnKDtAK7MSUj4UWMzSVx7soYgrVVOWoXJrWz4u5PvCOopDm+MHRPGl1qxchrhhHIUEqecFaJlJjs\n0CqclGStnAerUe9OFNOR/Ya5TncH/Srbl1nwhvTjnfXwVQaVYDYzfVPOw2YckOvKuktQSkTFSdiN\nKenvasalkVkM3p3HhaF7ONV7EEF/FZLpNdVtgpWkPNbZgJHwIsKJFRAQPNXXhlO9B1VtWjHy2l2P\nUlGHANYDusRkhyIKWRJVL6+TEYzem08sgF4/xAC0rcJI9glz1/4adO/3myrWJHrNDYQiqjigfOAy\nBxcFhx0iyw7dfY01CIWXQUByIkuDtVWIJzN4vGc/TvUezCH8p3oPqA6XaOC2auQ1S/S3Ux2iR2DN\nBHSJyQ5FGI3N6u1Xrw640disvFfmxKD3e6trHUumUefz2L4x8O2xOAjZDYfPwsr2MJ+VmNne9PIx\nmYXLHFyUJNihO9HVhN/u71QOG5+ymzdYm0nBwf8+WFuFkfCirT4B+kQ/HwKZ701Dr49OqC2N2rAi\nCIgqEK22tVKzW3mvVU8mq2vNUoXnc2Pg2+OdG5id7fQjHTh3ZQyDdxdQXkYQml3CTHxVySvWHfTn\n2N5KVq1ECDkE4K8BBAFQAGcppd8ghDwF4FkA7wBwnFLqFocucRRTVcJSDJzoalKu1uz9/FVbNFiz\nA3/xxgy+NziBBp9H6kqqJaVaKS6jB6s3JSdvGjKvL7v9koFvQ6sWs9l9wqtAGGFjf2UeYmLpzELC\nKiM91tmgiqexe15kzg3MzsZyivUcqEXvoYDq9sun+nBqngp9c1gD8PuU0tcJIX4A1wghFwH8EsCv\nA/hWgd/vwiEUU1VyYegezg9O4qm+thwCwa7v/V3NmlWzvnn5Dt6ciOGbl+/giaMtqrazEt493NdU\ng5ZAZY7rrBHDKMTYnXREuDqaDRbs3j9vaFDOl+HL5svKPpEx69qqWYyEl1TpTvjfWYlEz0eQMbvW\nrN/ijcFJ2wtzb5Y5UWhVi3NCgDNkDoSQMgBfoJT+e6uNU0qnAUxv/v8iIeQmgIOU0oubbVttck9j\nu/Pw8H8LhalYCkMTcaxk1jE0EVcV7wHUlb54NdK3Lt1WkpN97sQRfPPyHXzuhNqFdSC0VRjlwfa6\nHMYiG6Neeguz49F6RjRmOgEzNQAY8mX4svmysk9kCRAnokmEZpdV6U628yZmBK3x5nNeZC7bZuJd\nnBZeDJkDpXSDEHISgGXmwIMQchjAQwCuWHjmDIAzANDevvu8eKxiOw2dxfKkGghFEFlcRTBQhcjS\nqspYyR80VuSHeYjwycn6u5rxoZ4gRueWVSm3eWOdzCCt5THF/tqZf71n8l1PGeMxUwOAwaonkwjZ\nfNndJ3wOo0ON+XnaFEuQAewXF9IDUyWFwkt47uItS4kDnYRZtdIvCCFfAvAnlNINqy8hhNQA+D6A\n36OUJsxlmSo1AAAgAElEQVQ+Ryk9C+AsAPT19eVvft/hKOamB7bnpiIzKjMiynIYiQa3/i51crIL\nQ/fwnwdGUV1ZrgQlyWwWMilMy9uH93SyMv96RtaJaBK1VfoFX/TAB0wxAmKFKGl5Mtm5IemhkLmH\nSgFOVxZkqiRmrGcxPcWGWebQAOAxAL9DCLkC4BcAfkEpPW/0ICGkAlnGcI5S+re2e+qiaAdHlu20\n2JszWFulUvfEkxm8ObmAyGJaFewGZOeFubICAAFBdWU5eg7Uqrw+zBA/Zu8QU3wPhCK4OBxGR1O2\nHK/ZA6+XyO2HN2ZACNEt+KIHxVhpgoBYMbY7eUN1MihrO95npWIkbw/jf2t3Prv3+9G1v8ayMFjU\nCGlK6dMAQAipBHAUwDsBHAegyxxI1qjwFwBuUkqfs91LF0UF28z9XY1F8w4R3w2oaxIHfBWIJNLo\naPKppH9WDGVoIo7IUjao7VhnA0ZmF/FkTwvOvXoXo3NJ/LOHDipj0TusidQaEitrSKTWVJ/zXiOM\nCeWrr2e3ne6gH2cv3QEFVdKRm4GYi4mfD2Z/AWCK0cv0/06su5NBWU69z0qUvJXoft4elm86FdFN\n2wpeHJrKO3UMYNFbiVK6CuD1zf/M4H0AfhPAdULIm5uf/SGASgB/DqAZwA8IIW9SSp+00hcXhYPT\nXg923y1+zogpj4FQBM9fGcfS6hoeua9BIf4z8VW8PDyDS7ciWFpdg69yH7721IO67wCAWm85aqvK\nUetVHw0ZIdZqwwxa67ZSMX/36rhiKLfqmy7eTAZCEXznyjiSmaz2l+X5t8LonbyhFnsvmXmflSh5\nM4SdVz3K0sPbmU+j9+rdDqw4JeihoK6slNKfAZpVVv6ukO92YR/bqffVM/CJyc2ArG2itqoc6xtA\nZ2ONIhGe6s1+F/RXYnQuqUoHoTc+PYOu+JyVOTKK7uWD+/JBf1c25cjwVCInvfl2GDXF9A7bnU6E\nj6Fh6VbMEH7+ea2xOHlu8gnSy9nD+yoq7PSBULoz7Lx9fX10cNCNldvL4FUm7GDHkxm8eH0adIOi\nye/B0uo6Hn/HfqnuN9+grXywlXPHXHlIBjv9K6XcToD9sRfiXez7/q5G6R7RAm+Hu3gzjI5Gn+IE\nsB3zbeWd+6rr7q4vxzqtvsNNn+FiRyE0u4SxaBIjs9mMrc21HrTVVWF4ehGLKxn0dTaAgEgrtNkJ\n2nLq4NvV49sxZlqRds18ny+K6WVn9C4j+4AWeDtcR6NP5QSwHS7mVm4pdC2dsvMOlznscZSalKmH\ngdBWBtbTj3Tg2wN3MDyVwMl3tuLBQ/WKQRfI6l3fnMh6NwHaxeKNiIlTB99qziG2Jk4QVqMxFJq4\nGY2dTyZnxSBv511G9gEt8OvAivAwF2RxjcyeqWKdPVLusdW4yxz2OLYzsM4qRP35g231CIWXUeut\nyPHKqPN5EFlKq7xWrARt8R4tLMum2WjnfA+66G2S77qYlabtGD/NQq+NgVAEz782juXVrIcYM9QX\nqvAQYK+WM3/b5KPrZU4BZs6UnbNnZy02VpOmY8t4uMxhj8OMZKq1IUW3yXyJotHGFw/hyd5WRfoT\nn7ViiNUrB8rnzNFKeW3Xt15rvPFUOlt2MpU21U6+yMf4aRZ6bfR3NePSrQiGpxMqb7RSKjwk9pf9\nlb3D7G3Prour5bVYz2RMv4CDyxz2OMxIUFobciAUUaWtyDcFhNWNz3vCMEOhLB+NEWTvtaKCsuvL\nL76XzQVAUOfzIOD1mG7LynuswgnVlqwNfu2fOdmT4yZ8+pEOJFfX0OKvUqVBMQuz49YTfmSfa0WW\n83E5ZvevVTtXMe03LnNwYQgzsQeyzWo2upSlgDj9SIcpFY74/LlXx0GxgabqypxoYTPGWObayEo8\nst+aJaR23UXFeeWNnnwNi3yRL0FxwkVT1oYYtS7iofZ6PHb/flwYmsYhGykkzI5bT/gxGwAn7h0z\n+1frNwMhea10oLhu5vueffbZorwoX5w9e/bZM2fObHc39iT8VRV44GAA/qotd2m2sZ/oacH7375f\n9R3DS9encWFoGlUVZXjgYEDadkO1BzenExgJL2EqlsInjrWjtc5r6ln2/IWhewjNLmO/vxKHGn34\nwP37cSDgNdWHl65P48WhaaxTipXMOi4Oz+r+VtaW1vy8dH0aDdUe5XPxM/G5hmoPqirK8ERPCzoa\nqzEQiqCh2oPFlTW8dH0amfUNXB6J4Ma9OL70D79Eg88Dr2dfzntEyPpXCmDj5W+Nq5l13I4sKePh\nf2O1/2bHrfUOM+9m77g8ElHtDTP7V+s3DdUe3JxJYCyaRE1lue7+N4Mvf/nL088+++xZq8+5zMGF\nLZgl/OxwMQInEjF2uEbCWdfUmqpyNFR7cH0yjvZGL57oadEkuosraxgIRVBeVobx+SRa66qQSK2j\nqcaj9El2wPk2OhqrlYN4pLkG7zoU0CQGZogFa/sXk7EcRiObs6lYCs9fGcPrYwvobvHj0SNN8FdV\nqH47Pp/EhaFpjM4t4fWxGH5yM4yR2SVMxlLwV5WbYqKlCEb8meG/qcaTw6CLwdi03mGW6QO5e8PM\nXmG/6Q76cXkkkiM01FSWW2KKWn2zyxxctZILWzBzZTfSzfK/+/wT9+P5V8fww19O44fXZ7CUXsPH\nHm7LuZLzV30AKjVMIrWGWm+5qk+86yK79ovqArGYiggrdRdY27JqXLI507LbyH6r9I8C37s2gc+d\nOIImf6Vm5TdZ/wvhbZQPRMP/VCxlOuW4033MJ8meqO4xo/5hv8nHbiHrG2tHsWHZjJB2mYMLW7C6\neY2YSWudFzOLK3j1rXn4PPvwgbfvl/5WbIelnSAgGJqMayYqk+m3eT1xf1dz1uMomutxZMWwOR5d\nRm1VOY51NuR42MjmTMtuo5Wqg7XZ3eLHuStjCPorpdXwtMZu12lAa27yhVapVzvQIo5mmYUVG8N4\nNIkXrt3LO7kdb/PK18akZcMqq/TV2mnPZQ4uTCMfyUzv0LN2n+xpQXJ1Dfc11eBT7+kw5c7KEsv1\nttXmlP3kIRqNRYkNgBJgJzO883+1MBCK4OUbYVNpuPm55NONy76XzQNz83zfka3i8rxB1E7/9cal\nNTf5Qlbq1e4+0yKOgDmGZuU2/LUf/grxVAaJlC0vUQUXhu7h+SvjONpaq2TRtQvxbLBxnHPjHFwU\nGmYOm52DzatiHrt/v6VnmeRJQKQStJ5KiL9BXB2dx/HO+rzcR/lbgB6x5sfMYDWtB0skyALEzl66\no5mmOV8PF6veWFb2gIwgizU17NbOsMoUzURys374N7P3+r35kVACglRmA8PTiyovOyfAxvMJN87B\nRaFh5rDZUWGw9mLJtOVnmeQZ9Fepbg7sIMeSaQyEotI2+RvE5dAcWgKVGJrMZjPl9claY5IRrTqf\nB93BrMpHVMPwUdd8VblwYiXHbmA01w+11+Oh9nrFlhJPpfNK02w126hRxLPZdZS1TUBACFEC4uyq\nxcwyRbPMh49ct1KOVQ8ne1sVtWgxYheswGUOLkzDzGGzo8LgjcZmD5yY3iKWTOO1u8msHr62Stcw\nrNVnZvAVk7JpjUkkWuzfLYF5qRqG/34mvqrYR8R0DPycGIEf52dO3JeX6shOIjr+9+KaGN2etMBH\nvgPae8opI7nZsfN1EqzexvQC6viKg6WU68xlDi4chdWMoFrPGuUx+uoPbuCN8RiaairxlY8+gGBt\nVTZj62YQnJYqxEx/jnU2IOCryDFYi6lCtIiXVoUx2ffi5+x9AEzNm1MqHytMXauetuh5pOehpgct\nzx/e48zqjU4P4ti1cjrlc1soZL4loDBMxWUOLgoKu5td77mBUATD9xJYWM5gfimNc1fG8Pkn7kfX\n/hp07/erjM5m2s2V6HOJG4Acl1M9ZiYzRvO/57/XMpCbmTerEqzWvFpNNyIrYanHLL97ddxUWU4r\nfddiaFp1wGWQEVVZTqd8056bZb6y31lxsQWc8yZzg+BcFBR2I1z1nmuo9sBTXoaOxmoEA5X49Hs7\n8cb4Ai4Oz6K90YvZxVUlEEgMDGqo9mA1s46VzDpa67yqgKUP3L8fTTUeJSiJBWb1dzWjqmIfxuaX\n0dfRgGOdDaqgJUA7KFArMMlozB2N1bYjg+3Oa75tiEFjYvQwC+LTC9jTmy/xvVrBa6+PLWB4ehFH\nWwN4+HCD7lhk69bsr8Ryeg2nH+mwFGnPvm+o9uSMQS+Yjx9za50XDxwMqKLiz16+g1duR1FTlY2W\nfmN8AV//0S00+yuV/umta0kGwRFCDgH4awBBABTAWUrpNwgh/xHAKQBpAHcA/HNKaayQfXGxPbDr\nKaMnlcv0tFdH59Hf1agU+gHUdgAgK40NhCKgoBgIRZV6zaJEv1VNTH27SKTWUestzxqbNzOwsjb5\n1N58v6zEB2jFNjiJQuVJ0oOWOk0GPQlYS70kQrRZmOkb/1tm7Df6ndb3Ysp1I8nf6DY7Fk2i2e9B\nPJnBVCyl3GySq2sq7z6n90uh1UprAH6fUvo6IcQP4Boh5CKAiwD+LaV0jRDypwD+LYB/U+C+lDxK\nyRhVShAPj6gT5tUc/V3NmsZMmZHaTMplXsdOQJQMrN1Bv0L8P/Zwm1T9Vaj4AC3YqYFQ6H3HEy7m\nLJCP+kUrStnOWPLNoCr7njdcA7n1OcR+ysYsMtRYMo3LoTkEfBWKG3OLv8qUa3mZL9BoOEAJCsoc\nKKXTAKY3/3+REHITwEFK6Y+4n70K4GOF7MdOQSH0hjsRvNfLSHgR3UE/HmxL4tKtWTT4PPjm5TsI\nhZcAZKU80ShrFN3KEw5ZfIDYhhbzMSL+WsbiQhJwOzUQCrnvxH4bvcuKR5ys1KeWF5VWpblCMEbR\ncC0yC1k/tcYcrK1SXJZ5RvhQez3eGF/AzOIKGnweTZvOQCiCsiq/vm5NA0UzSBNCDgN4CMAV4at/\nAeC7Gs+cAXAGANrbdz+xtOMGulNg5RDKDMRvzS3hn25FMLGQBCFl6ArWKBKUHkGRHUJeLcEfXJFo\ny6rB8e8y8hRiHk2il1MhCfjpRzqQTK8haKEGgsxbyiyhtKoysbvHZcWcXhyayvGYkrU/EIrg/LVJ\nUEoVVaJW/5yAuB9lXk5G8yD2i1epnb10R4mNmImv4uXhGczEV5Xzwo+lv6sZGyuL83bGURTmQAip\nAfB9AL9HKU1wn/8Rsqqnc7LnKKVnAZwFgL6+PsvRPTtNTVMIvaEWij03Vg6hTEc9EU1icCyGhw41\n4OjBWttuimJfjnc2IjS7hOOdjfizH4fws9sRTC2k8NGHDioBdLztgcFKvQnRy4mPbjab91/mOirD\nQ+31eKw7WwOh3WTErcxbyornkt4z/Pzns+dkMSWXQ3M40dWkumX2d8kj4VkOLpEBOimQmSkOpPcZ\nDzHnFAPP6J7qa8Op3gNo8Hnw8vAMnuxpwejcMmLJtCIYtNZ5bZcJLbi3EiGkAsDfA/g7Sulfcp9/\nGsCnAHyUUmoY3m3HW8lsTYC9iELPjcxLyKyXDPv9G+MLykE73FSNg/Ve/MbDbUpaa/YOVudA5uEi\n8xJhHkuRxCp+8qswpmIrqCwvw2t3o5hfzqDKU4bJhRUcbqxWUniLKcet1Jso30dwtDWAd7XV4fJI\nBD2tAfzGw4dwIGCubsVL16dxcXgWvYfq8OiRJsO5F+daXAuzaafNQO8ZkVgajdWKpxL7N0vxrecJ\n5a+qwMOHG9B3uEHxnGJ1Izoaq9HRWI0Xh6YwODaveLDZgZNn6vJIBK+PxVTp59k8VJQR9LTW4lTv\nQTx6pAmDY/N4fSyG9kYfVjLrOD84ifJ9BH2bnlp/8qdfT33pj77wDat9KChzIIQQAH8FYIJS+iXu\n8w8D+AqAJymlcTNt2WEOTrjt7VbYmRsrbpniQdFy5dNq8/krY/ibV8aUTc7n/hcJtBkXSf59A6EI\nVjLruPCLaYTjK+hprUVTdSVG55Jorffidz/YjfZGH57oaZHWV3jgYACZ9Q2Mzi2pCgvJxuKvqkDf\n4QY8zBEmvp9GhJz9ZjWzjjuzSzg/OInWOq/yTq3x8dKr2HcrRYuMoPeM+B6jseoRVy032X1lBKNz\nS3iypwXtjT7D/SwyFVYv47++Oo4bUwk0+Stz3q21R8XPZXtC73m97xuqPVjJrGE1s4EDHMPiGZ1Y\nJKq/qxmh8FKOG++zX/ri+LNf/ONvaq+iHIVWK70PwG8CuE4IeXPzsz8E8GcAKgFczPIPvEop/ZzT\nLy+Gmmanqa4Y7MyNHdWQ0XVdq00xv47st1oukqJBWzTQXRiaxoNtAXTtr0FTjQdt9T4lKtrIi4ap\nJWLJtJLygnng6OVxEudEKyGgbD5a67ygoDh/bQKra9mKdX/5z49rzqdYYlLmucX/NQO9OdWC+B6Z\noT9fewTLrTWfTFtSoYmpWpjaSXRFHghl65NfDs0p/dTqP59hlvfKGghlS+FeujWLBw/VK9lX+e9l\na17n8+DC0LSijjSjtpK68ZZi4j1K6c8A7nRv4aVCvreY2EseRlYOr1nmo9Xmsc4GjMwuorOpWpXC\ngv1WLyJZNGgDagMdkCUGiZU1AMArb2XtdSd7W6WHbyqWUozKr41GMRCKqtxhreRxEvvI+qTn1ghs\nMUtKKe5Gl3D20h3FwC3O59W7Wd/4gU2bgywdBf9OM947enOqBaM9YMQ8zMCuzUB8F5/JVkza2N/V\nKHV91mO64vpevTuPq6PzCM0uI+DLSvyy2ul6YzObLsSNkC4RFFo9U0qQqXasQEvtIivFePbyHYzN\nJTE+v4wf35xFRRnB+9++X/mtngpCjHjm14a970CdF1UVZSgjwO3ZJRyq9+HizbASicpHuQ6EIvib\nV8aU6/q7DgVUdZ5ZJPW72uoQXU6jo7HaUM3S39WcUz+5qqIMjx5pkqppDtR5kUqvI7K4gvQ6xe3I\nMirKiKreMj8+oxKTWvMnfs7WjI1RNqd24UQJUKM27Jw1NgfM5qSlXhTXiu8LTxdY1HNleRmOtm7t\nnZvTCYzNJREMVOHjx9oN7WVatEbvLEzFUvi///w/pe3YHNzcSnmi0OqZUkM+fWfXa17lofU7FmjW\n4q/CSHhJ5SMuVs+SRVDLchjx4NUL7Y3V2SCj20uo2JetxcCPszvox6EGL0CJqsIbH0nd36VdSY5B\ndH3lbzOxZFqJgJUlGhwIRXCwzovWeh+a/R482FYPCmo7T5KWxK0lrfIeW2bdb+3A6ZKm4n410764\nTlrf6YHfX9+6dFspYcvf9lh5WrO3Hq111evTQCiCMm+tsReDBC5z2AY46T5XaMj8ywF76ZMVlcfc\nlspD63f830ONuamvsy6m2cRvjEjLImaN+sgf4tDsEsaiyZyMqQOhCCYXVkApVaXVFn/DAuK6g37F\nF50PutI63K11XhAQnL82CQA5KRdYIrknjwbxsYfblP4zvblRNle7mXFla1FoWBE+tMbFJ9471XsQ\nQC6z09sreszVqjDIXJkTK2uorSpXxVk4pQIyDPxcz6zYaXdHMoedagRmKGY8Q77QCsYx+p0MrXVe\nU9KSUY4h0TjMgtRkEbNafWQFeVjAm9g3UbIXDZbiHuwO+tHR5MPpRzowEl7E86+NI5leBwExVWM4\nnkpnbw+pdE5fmb0h4PVIx8XsILFkWklbzuZAye80p77RmCXChdireudXz64kQistxdBEDGsb2SAx\nMR+T3l4xShMvfscHTYppQXgD/lN9bcrNQWvfF4qmtdZ5sbYYnbLz7I5kDsVQy+S7WDudgTGYlRyt\nXrftQqYyYCqPqVhKN+Ea38fnLt7KiVLWk+xFAi+mheY9VbqDfvgr9yGztp6jDtPaEwGvB3U+j1Km\nVJxP2bhEQ3g8mVHUcbykPBZNornWo1JbMYbHB0zpzbUTKh7RYwzQTq4HICcoT2xLTEvBGOHkwgr6\nDterajIb7RWjJIky9dQX//46QuFlJQGeGKT3wrVJdDT6dFWoYvvxZEbXa84s+LkqSW+lQqEYV918\nGdBOtivwMEvM9TxgnIToBcL/ZX049+pdjM4l8ZkT90nrJgDqGswyGBFG0dVWVDFVecrRXVOJoYk4\n3hhfUBIEau0J0QWR7+tULAUACCdWVO6kfMW1kfBijjsus80cP9yA0bklXLwZVtWiAIAXrt3TvN04\nsYdlbdgVOMS2xLQUjBF2t9TkEGTxtsm8jXhBw2yeLPauheU0KvYR3NdUI/VYMqNCFdtnpXJZLXJZ\nPigGvT3Kq89KMvFeoVAMtUy+DGgn2xXswqyrXT7g0wrIGNJAKILvDU5iaXUNqcwaGmsqpQnuZGmZ\n+b4axSzoEfPuoB8djT5El1bx8ztz8FXuy0kQKEMsmcaLQ1M5LqqiG6lWUSIxzQezzdR69+HGVAJH\nW2tV75YlhONhdg9rra/oOMBgReDQs3/ofa9lR9AqqCQ+K3Ny4Pdaf1ezQsDZeoljCvqrkEyv5aTA\n0BsrsyPFkxnNfFAMesybjSeezJR+4r2dhnwZ0E62K1iBnt+3URCRHfDqG5k03t/VjPHoMkbnkkil\n11WqIzNMSi9mQcufnI+DONnbqvTxwbYAGmsq8WRPi6Lv1gIzXPI5mFjbjMAe62xQbgyi0TyWTGM8\nmsS3Lt3G8c5G1c1ieCqBVCaGzqYalZEaAJ7uO6RSv/Aw2sNmgsR4xwE9yNZGzwtNr798YkUz3kbi\nerK+G2V81asyNxCK4EfD4RwnBiPwYwCga+fSY958O79Tyon3dgN2iw2BwSifvFnI1Dzid3wQUb6Q\nXd/5v611XvzBh98BALh4YwaJ1Qye7GnJ6asZTymtq7r4PCPsaxsUI7OLqqAm3puKDyATbzzdQT+6\nW2rQXFOlsgMwAnuiq0kVmSyqywDg/LVJeCvKEJpdUt0svnXpNmqrylHrLVf12Szh1gK/vr1tAVwa\nmUV30C/15uJh1g6ht15mouBlz4k3AK0brhYTkY1HNiZm0xGjrvnf82nEWb95Rsgb2VlUPnM8YOo0\nM4GIG8l4VPdHGnCZg0mUmg0hX2YljsfumERdu1Y6BL0+2knNwKAnTc4n0/BWlGN+k9iayWwqa4/v\nn1YkayyZxtBEXHGF5SX/8egy1jc2cG8hiYfbG9DflVstbCS8iERqHT7PGgZCURBkbxANPg9aApWY\nXEjie4NRxT1TXHsCAp9nH3oO+BWPKdZPKymjjfYVPxfsRnOytxXPXbyFn9+Owucplxr49WIPtPqj\nR4xFVZtRW7Jn+Wdkn+sZy7XGxBNuPS+1gZA6jTgAzRggmeOBFXffPWVz2A6Umg0hX2blxHiMrrlm\nVWt6B13rt3q/YX3jmcFAKJviub+rUepyqPVv/p2y1N1snJ997G2Ka2N30I+pWArPvzqGS6FZLCyl\nEV9Zw9oGxcF6H1rrvDmuq8yWwtIuvzm5gEgijY4mH2biq6jYR5BYWUMitZYzB2+ML+DNyQV85J0H\ncPo9WQP7SHhR1T+ze8RofsW1YjcPozTkzHuHpYvo71IH/cnsR3pODmyPaeXV0hMw+Gf10ncbqbrE\nWzMTECKLxulFZDcL3oAt2jfYM2yP6nnl8Qz83JUxNwiu0NgOG4JZf3ArzzHYGY/YrkhIzOqozR50\nGbTy3ItgzOBU7wEAUKRcAKo+i2PIx7uGt4dcHZ3H//vKXaTXN9AdrMEH39GCuaUVhYiKrquvjUZx\ndXQeydU19B6qQ2QxjZrKckSXsrYLCorRuSWMRpfweE9QdYP59sAd/OOvZvHBt0OlxuLHwGBkcGdE\nTitam72T1RBg68Ab+GUuqPFkBs1+D8bmsjcrMakcWzOzqiV+r8nyavHvVlw6saW6EY3Tsv1r1B+R\naNf5PIgsrZoqCytzj+bjbPRuMUbnjD1bWzWL4elFYGNt77iy7hUY6U7zla7t9oddnbXULHrPy67O\nWgedBzvk49FlvDkeQ9f+Gl0jH0/E/tfvvI7MGsXp92RVXRRUIX4i4ZcxH7NMj6mAuoN+vDYaRcBb\njvpqD77y0Xfm9FX0diIgSGU28IvJOACCx9+xH0MTMfz8ThTein148FA9OpuqMXwvgddGoypjaGdj\nDWqrYuhsrFGNRSYBszV8sC2A2qp9GL6XwNmkOokfixTny6GK0jwLrJMZW8U5ZIy6ty0An2cF3UE/\ngrVVqj6y/48l05iIJpXEgjLGbEVo4s8CoBYMtBghuwEe62hQ1lPWvp63VDixgucu3pIGyMlgxdhs\npIZlz0xEk9nUM5nVZelLDeAyhxKGXdUPu7LqBTjl0x+m87RSIY153TTXVJr2/eaxJQ2V56Ty1sPf\nvXEPofASOhp8iCWzKhxRauX7IXpDidBTO7UEKpVnmZGRgCiEUIyozc7jPZzqPYiTva2goFm1xNIq\n6nwefObEEYAAK+kNXBwOo9lfiYryfTljP/2eDqWWAZBrdOUJZHfQj1rvLG5MJxAKLyKemkNjTaUy\nFwOhrL9/TVU5Lo3MYjyaxNBktuQKLz3zwXbinIyEFzE2l8S5K2MI1lYp/Yol0xia3JpbmcooNLuE\nwbvzqNi3T7o+/HzLAsZEO4fMzsTPk5i2BNiqwz06twxvRTmujs6riLDWvuW/44Msjx1uMHTxtnIL\nl6lhxblkt6ZDjT78JOV6K+062FVltdZ5EfBV4MLQtKaPdD79YTpPWc1hPWMfS38MQFNtoQU91ZPo\nSspUXudeHcfq2hoeaK1FZUUZfnxzVmWYlfVf7ztxfOxAisForA2WyoIROUZ0kqvrSGXWcHV0AYHq\nCmWNmN3i2wN3MB5Nor+rGY91789G2nKpOURhQbZPtNQfA6EIRmaWsLZBUe+rAKUUPQf8OYLIpZFZ\n/Px2Vu0kepo1+DzYoBRP9rQoaiy+hkR/VzMu3Ypg8O4CXhyawpnHjqh8+GVSPc+Yelpr8WDbVlyI\nllQdS6bxwmC2TkLvobqcYDFRtSgTWvhYD/Ye5t32ZE8L5pNpJTCN9dMMeBuM7Jakt0ZG4Pfoa6NR\nRZpYNR8AACAASURBVMgQ1YRsX3xiL0VI72aY8fc2A7u3Dr1+MMgCikRCKb6XXeGZ1Mv84s2mCuCl\nYZY1lRHjeDKTEyPQHfQjnkojupyGZ18ZpmJr6GyuRjyZQTiRzUP24+EZvHwjjPHoMgJej+JW2N/V\njK/+4AaGpxJKegw2J7FkGg+2BRBPZnDu1bt4+UYYT/W14bOPvQ1TsZRiCH5xaAqDdxfQc2Ar8Oz0\nIx1IpteQSq/h9bEYVjLr6PRVqxjRSHgRIzNLCIWXVbcBNkdm/OX52Ag2VzxhZevA4iZkknd30A+f\npxxP9rTgrbklnHt1DASA31uOoYk4bs8u4eXhGRw9mJ2L5lqPqoZE76EARsKLqtQWouFflOrFsTLI\nGDLPgAfvLmAkvJQjCIlMUUbgeU+uLceDA/jaUw+q5lP09jKCGGQpYyp6Kio98L8dCS9qxuXkC5c5\nlBis+ntrIR8DulayNj29qBmPHqbOYXEPZiUy0YjaEphHKLyEWDKNgNeDj7yrBU/1tYGAKH1iRs8m\nfyUeOBDAT34VxsjMEqLLGYzMZlVHtVXlWNuguDwyh4VkGqtrG8o7X70zDwqKV+9ElRsJu/20BCrx\n5mRceZ6lyDh3ZUyZMwqK8jKC3kMBhdAFa6vg85Tj3sIK3t1RB29FOT5z4j7lgMeTGcRTaRyq96G9\ncUsNxqTub126rbodyeaIrQWLYWBts3ZEz6rTj3RIvbQYcfvu1XElqygA1FaV48mjQfgqgzj9SIci\nnZ/oalLZUI53NmJoIobF1JrStnhLYFK9TCXEg7dhiO18/on78eLQFOKpLdsBoDY8sz3K/+X3pZXf\n2MHFGzP45uU7+NyJI3jiaIvtdmUCmxZDdQIucygx8Itt5F9vBVZuH6JOmf9cSy9q1oDGB4eZlchE\nP+/uoB9f/cEwostpNPkrFXWCVhBTOLGCn47MIrK0ipa6KkU9w1z9RmaWUO+rwEIyg5/+KoKFZBqe\n8jI011RidG4Zk7EVJeYgnkwj6PcoEcuMIZy7MqbMGbvuP913CMc6G5T03QQEY3NJdAVzc/+0BOYR\nT6XxvasTmE9m0JsK4E5kS+U0ENqKoAa2blzhxEpWt++vUmwDotqBT1/B5ujlX87glbeiiC6u4sPv\nzDLqHw/PKqohNu/dQT+ePBrE8PQi2ut9aK2vUqlvmO2EuWSyz0fCiwjNLmMkvKRKuc73zYzXGWuL\nzbFYNY15/bBbrFhljc2zFiHWinqXQfR84oPY9Jj1Ny/fwZsTMXzz8h0VczDjRs1DJiRqxV6Y8fwz\nQkGZAyHkEIC/BhAEQAGcpZR+gxDyJwA+CmADwCyAT1NKbaWVLTXkG5zGL7ZWzhw70Lt96N0IZF4Q\nsg1nxoB2oiujHF6jZ3iI7x0IRfD0w4fwvWsTePrhQ4bufwOhiHKLeOYjPYpUPBVLoXu/H137a3C8\nsxHfvvwW/vFXYaTXKbr2V+Mr/+MDeG00qtxIvvriMH4VXsw+15I9fDyjYX+/PXAHw/cS+OQj7Ur6\n7uXVNTzd14aP9bUp4+BVSTPxVXTv9yOzsYF1CtyLpfD5D6kDCflkbMxr7EfDYYzMLKK/e4txym4Q\nok2ooixbaYwQojDejiafohoCoOy99sZqDE0m0HOwNocQXRi6h6GJGCKLaZVap7+rOYdp2DX893dt\nlT4VAwxlwgAAw6R3LA7l4s0ZLK9uqNJ+y8BnbWWp0QfvLqC8jEjtevxe/NyJI8rNQes3Wm7UZjIZ\nyAzbvLHaLgp9c1gD8PuU0tcJIX4A1wghFwH8R0rpHwMAIeR3AXwRwOcK3JeiwEk30nztBmbbMusp\nIftcZgzWerdMjWQ1JmMrFUVl1pNkTO2xIRsjL62y4u/MaM1UGw+116P3UAD/eCuMDQpk1jfw2mhU\nJRV2NvlQW1WOd7UFsJhaw/nByRx7w9XReQxPJZDMbChM5dKtCIanE6p6DKK9hv0NeCvw//z0Nv7V\n+9+mqJMYE+FtHyw529o6RXeLH5/pP6KogZhxWBZsJs4Fi1fobKoGBUX3/i3jNPN4u6+pJselk+2b\n71wZx+LqOh7pzDUga+0HvbXi2+b3yulHOpQAQ4YLQ/fwndfGcWlkFs985KgyXxeG7iHor0LX/hrd\n9r9zdRzzyxn4PGWayQf537OsrQQEY9Ekeg7UovdQQLnl8+eAMcdYMo3+rmZ8/3feZzgHMq2BGI+i\nF5jIt8ELKyUZIU0pnQYwvfn/i4SQmwAOUkqHuZ9VAwYrs4PgBEHnCabTnkYymO2z7Br83MVbGByd\nR0X5PlVAk+zdMjWSFjPVunKzWIJjHQ24OjavGEzjyYyqDVEFAABXR+dVCeLEcZ/qPYiRmUX85OYs\nFlfWcX5wUiUVnn7PYbQ3VqPB58H/+d+GMR1fxb2FrHGb1Xd475FG9ByoxX1NNQpxfOZkj9IXUcUm\nGow/frwd/d1ZCZARHL5uBJvP/q5shbiPvPOA8h5VsFl0K9jshWuTSs6nq6PzKkI4OreMmfgqXh6e\nUUU8A1A83liuJlHK7+9qxqWRWQxPJdB7qE557vlXx/C9a5OYiCbxrz/8duleMlLzdAf9ONG15Y4t\nu2kQECTTGxi+l1BuCLz67bf7O3WZ0w+vz+DNiQU8eKgexzsbNZP18XuE/RUdKZhthneKMPIYFOdA\nJgQZGZqz6VmSqPXuQ4PPk2M3+u7V8dLPykoIOQzgIQBXNv/97wD8TwDiAD6g8cwZAGcAoL19+/MZ\nmUG+xiugOKmv7bQpdT2MJtFzUO16qAXZ3MgYk8wgPhCK4IXBSWyAomzT42kmvor5zYC8C0PTil6f\n/WWMAIA0AaDscJIygg1ki7kf2V+tSKq8AffclTGMR1OgAOaWssyB1XeYW0ojsbKOQ40+laGVXfnH\no8sqD6ezl+7kVDNj88zSMbCqZuJaMNdgpvtOpDJ4+UYYTx4N4rH79ys2iQ1KMRJewrkrYxiaiINS\nqrjZNtd6VBlf+XVgMRGgUIiPWP70mY8czWHCw9MJLCTTGJ5OWN6zbJ07mnzo2l+DgVBUJVDw7znZ\n24p4KoPR6JKyTrz6TbzJiILG7z7epYzZ6MYv7hXxVid7r1VBUU+FpDWHA6Gt7K+MwfNj6O9qxkYp\nZ2UlhNQA+D6A36OUJgCAUvpHAP6IEPJvAfwvAL4kPkcpPQvgLAD09fUV/HZRCOJrB1rGOidVVlbb\nFF0P2b8ff8d+lepFK12BGZURe5avaMbqMcdTaTT7K3EvlkJXsFql6786Oo/+rkYQEKU859MPt+VI\nXGI/+EjYq2NZNQzzTw/6KzE0mVAk1XNXxvCjG2Hcmkng3e0NOHa4Hl5POR7r3o9/ff5NPNnTgt/u\n79Q0BA6EstHhK+l1FbFPpDKIpzJIpDIqqfm+pjj++r+PoXwfQZO/Esc61cIf2yOJ1JpCHLqDfqxt\nUIzOJfF4T4syl2Ug6N6csxZ/GMPTCfgq9qHZX4nI4ioCvgppfQvmVruQzACgmIimsLi6pjAXZuwV\n901bvReV5WVoq/fixaEpPP/aOC7dmsUzJ4+q9oRsf/D2he79/hw7imi4bW/0YWgyrqwT88YS516m\n2+dteXZu/GK7/Hv1aInWd1o2M610KNlbwzLa6r3oOeDH4z0teRuheRScORBCKpBlDOcopX8r+ck5\nAC9BwhyKjUIQXzvQMtY5aYPQa1MrCph3PWQS71MPt2lucACW5pM929tWi47G3HrMTz/chsfub1Zd\nm/l+nOxtxZsTC7gxlYDfW67SQcsiqr99+Q7+6VYEr9yZw9LqOpKr63js/mbFW6edI2JP9rTg6ug8\nZpfS+MnNWQQDlfg/fq1LCWwDgK899aAq1oGfT6buEUtYUlCsrW/gxnRcCWQ61QvcmE5gJrGCgDfr\nsfTty2+hs8mHgNejqhnR4q/CoXovKKV4rLsZEwtJDN2L4dsDd5BIraO/qxFP9ARBQRGsrcKhRh9e\nvD6N6/fiONHVjJXMOoanEopNgmdQsWQaHz7aghvTcSVyu7mmEk1+dUyDuF/a6n04EMj+paBYXFnD\nq2/N48LQPeWWB0ApaDMeTSoxHa11ubW8ReLIGO3Vu1m1osweorfHZXvezo1fZiNgfRZpiV40NNsj\nE9Ekaqv2qcbCx6zIBI6Xb4SxtkHRWONBsLYqh8EPhCKlqVYiWb+7vwBwk1L6HPd5F6U0tPnPjwL4\nVSH7YRaFIL48zN5MtPrhhMpKhFF0rZZkpVVFjP9dOLFi2l2RfzaezGBoMqHoxzubfHgrsqwQfB5i\nPx48VK9Ey7KxfOe1cSTT2RgG3iPlvqYaDI7F8J7OJpCy7G1BTLvM1mk+mcZ+fxWaaz2YmE9iZGZR\nUTMBWxGxbO7Go0mEF1cUBsdiAR5sg2rOAl4PyveVYXQuiUfva1LqIrzncANmYil0B/2YW07jF5Mx\n/Pfbc1inWfXRpzazr8aSabw1t4xEKoO55VXMLqaRTK+huaYKj3XXqohRIpXB6NwyDgSqMB1fwS+n\nYrg9u4zbs0sYn1/Go/c1gYKCxZIwG8Sn3tOx6ZUUx2QsBa9nHx69r0kJ3hOL/Yh5o4Ym4vjFZBxD\nEzF85sQRxRX1wtA9UEoxGl3Ca6PzihGdeYgxNPg8SGXW0ODbqq/NvJGYKkVMcaEHp86RaCPgs86K\nAX78mZJpBgZCEby8eQN8bTSqOAt88/IdJS+YOC6myhqaiGNsLokXh6aUzAVsLkpZrfQ+AL8J4Doh\n5M3Nz/4QwL8khNyPrCvrGErEU6kQxJeH2ZtJvgFs+arGzEhWshoBDOPRJL46cgOdjTXSG5BWH3nD\ndcBXoXhqnOhqwocfOCBNacH349yrYzm3GcVoei+Rw8g+9Z4OxQefSc3tm6oY0U9eZHqMMTB1zBvj\nC4p66VTvASX9BMuyyuoeMJ06U8sc62zAm5MLaK6pymZejS7hZ6EobvgSqCgvw43pRfgqyvCutjrc\niSxiLJrC8HRCpZoZmojjymgUc4sZgFKUEwJ/ZXkOY2exDe+5rwG/8/4jGJ5K4F5sBesbFL+YiGNy\nYQVPbarjxKRuzCOLtwWx2I6salHblvOZE/fhqz8YxuTCiuKKyoLVWCwIa+vclTEM3p1XMfOXh2eU\niOwnjraobhfM2yqeSmtWGyyUDU/qbrsZjyEa9/n9w2x1LPcU21/M9Zcg616cyqwhFF5CV1DudcWv\nC/NskqW+L8liP5TSnwHSDGkvWW1rYdnZJHLbgULfTIAtrxneuwUwfxB4tYLe77UY2EAoW8QksZLG\nB9+em5eH/caM8Y/3bmLqBd5NkxEv1sbrYwuYW1zF62MLqrZkRlPZfPDvHZnNBl4xlVQilT24LyYz\nONnbqkqvACBHvcTST7T4q1SxBoyIhMJLeO7iLXTv9ytFfn58cxY1leU4EMgGlr2j1Y+eAwHUestx\nqvegwpToBsXXXr6F8egy/uDD78AzJ3sUyf718QVkNigSKxkVI/34sXZMRJP45VQCRw8EFKnztx49\nDILsbYaprNh8sFgQvh1GlGPJNEZmllBRTnBfUw3mN6O5ZRgJL6KMEHS3VKskaTYvD7XXK5lLu4N+\nJNNrKmYu3s74tTp76Q6GJuL4UE9Q06tHPHe8bevizbC0wI4ZiPuYzY1RABrPSJhqjgXzsf5RUITC\nS0ivreNzJ47o9o33XtNyLrCDfc8++2xeDRQLX/nanz37zl/7dTxwMOBIe1OxFF66Po2Gag/8VRXG\nDzgAf1UFHjgYMP0+O318fWwBQ5NxbGxkDzd77qXr07gwNI2qijLNOWSS4Su3oxhfSOL1sZjm77X6\n1lDtwfoGRX11BT7TfwSP97Sovp+KpXB9Mo72Ri+eEL4T222t86rmq6Hag2ub45uYT+L6vYSqfz+4\nPoWJhRS8njL82juCynOLK2sYn0+io7Ea/qoKZZyXbkXw5ngMkwtJ7CsjuDwSUb23pqocq5kNnL82\nidfHY7gxlcBbc8toqvGo5mQqllLUXp9+bycOBLw4EPDiQ0dbcLipGiuZNaxmNnBgM1L3gYMBvDmx\ngKHxGNobfTje2YB3tgZwdWweSyvrONxYjaX0Op482oLffbwbfYcb4K+qUNr8+zenFDfVj7yrFYsr\na5hdXMWHjrZgOraC5fQ6qivLMXh3ATdnEnjgYACLK2u4G13Gw+31+I2+bODgxeFZHO9swJnHjuB9\nb2vGw4cbsLiyplpXcd/4qyrQUO1BKLyEZGYNqfQG5pZXdfdKQ7UHNZXl+MD9+/HG+AIaqj3oaJTP\ny4FA1pa0v7ZK2R9s3AcC2dvd81fG8PrYAg7UeXErnMCNqQTe3VGH33pvp/SciOeOjam90Yt1SjE2\nl8TaxgbuRJZU+1m2x/nPqir2YXRuCR+4fz8OBLzK3LwxvpAjVPHz+OiRJmV/9Xc15/TZX1WBO5El\n/LdfziCzTtHTWosHDgZ06cFL16dxcXgW7zoUQE9rQLXfv/zlL08/++yzZ3MmxgA7Jn1GnVe78pEd\nFMr47KTHk50+nuxtzUq9UXV0qJlbC7vu6mUAFfvG13Zgv/8DDd929pyYU0f8ns/wKaqdWEK3ziYf\n2hurVcbATx3vQGJlDZl1qhq7GCT03MVbGAkvoWIfwfB0YjPVw5K0Ch3zEmI3h4A3V5U2EIpgaDKu\nSME8WP+/NzgBCorPPvY2tNZ58WBbPULhZSUw7rtXxxXpWjb3/L76zIn74Kvcl2PnONULPHOyB+de\nHcPwdBw1VeWqiGfxBsPmQ7aubB7Y97xKj60hy6dkJCkzKZ8ZltmeAbIJGMX4GKNoez6egFcraqWO\nED+TuYnKSm/Kzp/ocCGqTcXyr0DWK+7SrQh627aSMBqpjnk1k8xuwdaB0RlRbeUEbdsxzKF+U6Jz\nCk6peIy8FJzoo1b6aBlEbw/+c6OoZH5OmDeQUd/YoTJT4lN8h9b3l0ZmMXg3690iuiXyxIAPQGIB\nQx/qCap8zcV3DoSyka7dwRpVUNjxzsYcgsIIrl5aBTNjYjEQvJFVNNqKf0WPJ1FdKHqxsaI1ABBe\nXMHr4zG8722N+PADbdK54GMv+LkU+yESdvE3vArKCOKesZNJlI2VrbE8en5rL8aSaZwfnER3Sw0S\nqXUAW0yPjZ23c2ntG57xsKjzY50NOTnPmINEIpVV67E8YCMzi/BVtpimYbyaCcj1WhL3Az8PTtG2\nHcMcnIZTxmczkpaZDaGnA9cr+yiD0di0XOms3HhE47Gslq8srYZR33ipWuZ+qhdIx5LHdTT5NJ/h\n12ckvJijY2d9l7kVitCLZOe/ExmBbBwyIgdsrY+MwfDP8gWMxFoCIhHkbwAySVS2PiJDMLMPtWJe\nGMRkfUbgAxFlzEhc2/6uZrw4NAVCCDoba1Qp0MUzIGOW/GfMI+tUb26xKB5MeOGNwyzFCW8zsQpm\noznR1YQXh6bw6ltRVbyMlQSCZrFnmYNTMCNpGS0Un9RL9nunDdmy9uzeeLQkUPH6b4WYHOts0C2g\nLusDI3ysxCWvVpL1V5zvqVhKiTQmIIpbIVMZaMV+aM2ZTGgwy3xl6yNjMFrP8Lc+WYyAVhEgs4kZ\nzUIv5sVOMSre6C9jDq11XiXxYWdjDdClnjc+UFNWIU5PjSRG1wO5Bm6eMDNnChakyeJM7IIXgF64\ndg/ptXUc62xQ4mUKoSbfMQbps2fPPnvmzJnt7kYOtIzMWYNVGZjByciY9MrtKDqafPjEsXapgcqK\nIdtOn8X+6oEfy+LKWtaIfSeKmspyxSDZUO1BMr2GMkIUY6JRm8wYHgxUKQRLnC+9eWTjqqmUG/oY\nZPP90vVp/NXP7+LKW/MI1lbCW7FvM1jtoNQwazRn7LvuoB+XRyK4PhnHxeFZ5Xkz4+A/N9oDZvch\n/2/e4J9Z38Do3BLe2RrA4Ni8ql96zgxa45iKpfCLyRgON1bjiZ4WdDRWq/rB3seMuWacL5r9lVhO\nr+H0Ix2az3z9R7fw8o0ZjISXMJNI4dEjTXj0SJPK0M6M8YcbqxFeXFGel60n++yJnhY8eqQJ/qoK\nTQM3Pz/sNwcCXtyOLOHFoWnFOUBm3Gb9y6xvKI4R4vo3VHswEl5EW70PfYfr8en3buWO0tuLX/7q\nv5979ot//E3ppOrAZQ4FgpkNxNBQ7UFNVTk+cSw30MUJmDl4Yn+NmBkby/h8UsrY/FUVmFxIYfDu\nAsYXkobM7aXr0/jprQjKyoB/9lAb3hhfkM7X81fG8DevZNNK9B3ODfw0w0hl891Q7cHN6QTiqQyq\nK8uRWFnD8c4GPHqkSflePHx672LfXR6JKJ4xvYfqlOfNeo9pEV6zXmxiH7X6fHkkgtfHYggvruCV\n21EVIdMjPFrj4L1nGFHl33t5JKJ4xT1wMKBIvnrzIXotMWGipmpLKGn2V2IhmUZTdSVmE6vKd3w/\n+7uasZpZxyt3orh4M4zK8jI8vOkRZpYpi15LN+7FUUaAzqYaKYO+OZPAWDSJtfUtryh+zOPzSVwY\nmsbo3JKm5xeb0+OdDTmeWWydZALVn/zp11Nf+qMvfEM6qTpwmUORoHfAzBC0fFxv7RAi9sxqZh23\nBRc/fiwdjdWoqVK7KfK/Y4eCv1XI3tsd9GMqlkIitaYUpZfN1+tjCxieXsTR1gAeljAHO4yQfZZ1\noazEh462oKnGY5oR6EEmeYpStd4NR4vwGq2nFsT54ee/vAyoKCtDemMdU7EVhZB1NFYrfWd4Y3wB\nX//RLbyzNaDo8q3cRMW9obXerH+Z9Q3811fu4tyr42it8+KN8QW8cidXKDkQ8OIj72rFsc4Glauo\neGO6HVnCxZthLK+u4W37qzG5kJLuGa399NL1abwwOImbMwmsZNbxs9tR3J5dQpO/Eg3VHpW7reIa\nXZl1jWY3SH7M7Gb1gfv35+w92ZyKLsd6++LZL31x3L05mMROiHEQYZYgyMYmXuHNtM824kpmXdnM\nDdUevHR9WkUs2ME7e/lOjmrJX1WBZn8l/v/2rjw2ruO8/z7epHiJpESKoi47pBJJtnzIshGHjgNE\nduLacF1YqBM3DXrYDdCgOZwUhmukdoP+0SJwarcNEKUpkKTKYTluYzlODSVNJCawZVG2GFmyJeoi\nRVISL/HSUuI1/WN3VrOzM+/NvPeWu2vPDxC42n1v5pvzm/nOd89NoJAIaxuWpN1MuLilobIED92y\nOrlRyD4OHCtqy9FQWaLdVMNsnKIoICp/FHHc+XO/7xtDR/dI8lStg26DVYksxc3Ii+4f7e/Ff73e\ni+ICws1r65L91VAZPwHv6x7Grevqce2ySnSdHUNnz8WU064ovtlz9AIYgC9vW58mOtWNn9wvfuP9\nyuFzeOFgH/YeH8JvTwyje3AKl+fmk3PloYR8Xe5/vxtT3ZISlBYV4KY1tSgrLsKeo4PKg5BuPonM\n7dpllWhZWo6ZuQWUFBF6R2LY1dmHrr5x9IxewqaVNVf9OGrLlWI9k7kntkEn5hR9R3gZT//91/re\n034OUSJXAuzZwFQprWqbX9Ytr/K3rqtPS8Auxx7i5qFyWlFe9/ELU+gevIRV9RWeij4/KwsTxWjU\nynsTmM4n/pypCaeuP+TvufJ/boHh+OCkp7evHItK1V98rIcmZ5LJbXZ3xXNcc+WuymvZdl2ZWNW0\nt8a9iY+fn8L1LTUoLy5K5r32su7ymyuyRzIP1+IXU0wsW/RH6eiOh84+ezGG7VtasH1LSzzm0cjV\nmEfiHH+r9yKe2XMMd29owunhS8kw6LwfdXSLJrWiGS3/XgxVb2PUocL7kjlkYwMJC1PzNN1il7/z\nK1+085cXSu9IDC8c7E86+siWMjI9cpx7mSZTPQunSdyk5HejMuOzgel8CtJm0/rFAGyypZa4Ucox\nsWT/F7lMAiWtx+TNUz5otLdezX4mhroR6z/SP47nftWNG1fX4nN3fiCQ3w536tM5tQH6MDIyZHNb\nOV6YaH307b0nkvlCePRc3ndtjfF0s+saKtJymquYDre8Oj18CRPTc5idn0f34FQyh4X4rEgnjzcm\npw6WD1tinUGRV8whKu/jbGwgiwW5bUH7TF5wYjm7Y/0pp0+v/uTBwXR0yrF7dBB9DyamU6OARumV\nHgSm88nruTBt4H0sb3Qc8kapM1WVo6vK9vxirCsdHSoTVfFG8XznWbw9EPdKX7lUPTe8+kXl6PbH\nt6Snz/TyCxEh33a8+oabZm9X5Ao5fmESE5fnsLp+SZpPz1u9F9E9OJXiM8NvXPzmcKjvInqGY2is\nKlOGIOee17oYUvLBQ1WnLfKKOQQVB2V78wiCqGgO47+gOpkA6RFZZVptaLcRyXDfg3hI7hUpJ0mv\nMjIx/qJfhJjsKCiiEHXqmI/fRslDliyrLsHmlhqcHYnh23tPYOu6+qToIky4efG7uooSPPXSESyw\n9Ox2InQMi5fFo7HqNj/u3+DnjMrL0z0nnsjv3tiI08Mx3LKuTnlrktvN4SXWbagqxbaNTcn29o5c\nSkSlncOTf7AxSQsDw+zcPE6PTOHh29b43pjFOlFYHEjRmVcK6ce/9HljW3wRYZSU2YJoDRFGkW3j\nvyBCVrKKdvtyYDG5f23625S+uiUlKC4gbGiuxn2bV6ZYz/iVkYnx534Rr50cwanhqZQghyJUtuyi\npRBXJF/fUptmpRKV4YSsxJdpOtw3jsnLsxiamkFRIeF3J0dw6OwYRi7NJM19TeejnznotcsrcdOa\npRibnsUnNq3wNZBYXV+OtfVLcHl2HoUFlJx7b5wexf+9O4Ta8mKl1ZpsSqxSNquek/1Q2hqrUFxA\nmJyew+nhSxiamkFjdVnaPOLlqKyIVPOTK/Qvzczhro1NyTX20qEBnBqJYSI2i+VCPc215egZjWFg\n7HKK0Ydufol+Jf/y7HPvfVPWx/7mr0OZE9pukNmEiRmoF/ikUZkhmkAVjVNeRHVLSvCj/b0YnLiM\n9Y1V2LahKbnZ6KKuyjC14qoqi28CPEKpTRszMf7cL6J/7HLaQhYh+4TITPQHr/Xg6LlJrEkoDy7u\nZQAAG2RJREFU62VzySBMTWYqch+L5R7uG8cP3+hFVVkx7mhrwCc3rUBteTEWGMPA+NWNaHZ+AXu7\nhzA8OYOa8mJjelQMjvtUyNFt5f4tS+Sx+M3xQbzVM5YSKfjy7DyODExgQ3N1ir+LWN+JwSm8+GY/\nPthUidLighQnRFV9V2bncXl2PiVGWkNlCa7MLuAH+3sxPHkZm1pqlI6qHF4OcSJjvq65BgxIOvTx\ndzvPXERbYxW2rqsDQEmrI17G3PxCijWSaHL+8uGBpMUgj/xaVVYc2JQ1r8RKQZGPOgZdAD1TyGIK\nE7GA+Izumix+39Edz93AGMMjd1yTlAl7RV3VwVb0owo5oisjE+PfXFuOJ+/dkBQt6cbIT8yiUtbr\n3jXtIz8RlVju7q5+XLoyh1PDl3D3piY0VpehpqIYj9xxbTJHN8/2VgBCa2NqJFw5KZAJLSbtEsPQ\n6CIFq/QfYn0HzozidyeHUVHaiIdvW5v2vFyvqC/hYzMem8X49AxKCuPBP7mllEnwSlH0yAM7coVy\ne2s9bllbh8bqsmRZdRUlKe2U4zfJMbS4mG08NovXTo6ge3AS17fUps+l+dlZ7WTxwPuCOWQTYRWN\nQTc1eQGayLT9lHMiTTxKZDwKKpIWKjaWYGLf2MrcVeazmTZRlsdSNIfUQR5DOdiel0LWS+8TJvKt\nWC43n+QMSrZS4wrg9tZ6PLilxVcxDPgfMnTtEjOh8YCI8vs8Wq1ubYjPc32EbPqq609ZqVtbURLX\nxVSV4vpVtRianEnqDXRjIcYa6x2JYdfBPlQUFyRDwnMzZjFEOAC80NmHBTAUgFIYoMhAd3f1Y3J6\nLi2f9PHBSXQPTuLK3AI2r6pJ2Wfe6r2IotqmQBH/HHPIMKJWoqu+9zp5cZhs2n7KObldPD8AgBQL\nFdMge+LJ34apyM9zGm3LsEW2/WNMx8fmUCEzKNkcVdXPMi1yLgeTQ4aqXTy8thisUWzLjr0n03Il\ncKgi5IoBCFXPyr4CqjVz4MwoekZi2PahxqT59MDYNHpHYqguvxq1V3XQqS4vREVJATY0V6fcevjN\nQ/RD4NkB+c3My5dly9qlKd/3DMdwfUst1iV8UkQz4p37e0AlFbWeA6BBRpkDEa0C8H0AjQAYgB2M\nsWeF3x8D8A0Ayxhjw5mkJSxMbgCqkMJBNyzdRuQVOVJ+VoTJhiFe5f3K8xKX+IEzhuPnp9DWVKlc\nDH5QPZ9p8WGmmY8fbMYnTB2yOapXdFsT0YoqdDuHuFEfOD2Kuzc2KhMqAakOfPJ6lKPNelmSqXx4\nZHp4uaJol5fzkwO9+Pnhc4jNzOOGltG0m4SKcTbXlqekQpVFcap6RHAxl+zLIot5ZVHUw7euwbMz\nsTHlYPsg0zeHOQCPMcbeJKIqAAeJaA9j7GiCcdwFoDfDNEQCkw1YFVI46IZlIvP3ezYoTMrzEpcA\n3sy0ozueca6tqTLNmzfTZsfZEvMFgc5cNtNMSr496G6oXqIV0XvZK3Q7L4OLp+7bvELbx6IJtZdI\nSNSFic9zur36Ty5XNebtrcuw99ggjgxMKD3NxXfEm4vcVrHf/Px+dL4sYl2qdt24einmxs73KDvU\nBxllDoyxcwDOJT5PEtE7AFYCOArgmwD+FsDPMklDVDBZkKqQAkGh24i8TsymDmVB67aBFzPViSpU\n7/lt5rabfZSioUwzMnmTkzeSTNHRXFsOBoZdnX0AkBQ7eW3IOvATr5/SXZUsyktUKtct/yam2DTZ\n8GV6ZFplep68d6N2k1aB6+g2t9QAYGhbXmXEnFT1e0kHdLe5IFg0nQMRrQVwI4D9RHQ/gH7GWBeR\ntwejCmEWQ9B3TTbMG1cvNUqXmClkWyYuwmvjsFmcfm2ybXOUp27Tum0txcQbgpxHOAwdNlA5znlt\nyLp2qE68Oksy1Slb1yavOcQNBXS6BVOIYmJuPcTpsT1AcR1dU01p8oakEx+pbm2q/lD1aZRzYVGY\nAxFVAvgpgC8iLmp6AnGRkt97jwJ4FABWr77aUF3cFJNF6NV5uexJHSToXJD22Hg7e/0W9PZhq0g3\n2exlOsOeukXTQ1WoAxlBLMWA9DzCOrQ1VqGpZjRUqAQZ925uTopM+CZlMqY2ujJetp/eIgiSBg/D\nMTy4pcWYeXd0D+GFg304cGYUsStz+N3JeJyjL29br6XHZm2qFPgiVDof8X3xvY7uuGf7gTOjSfFs\nlIefjDMHIipGnDHsZIy9SETXAVgHgN8aWgC8SURbGWPnxXcZYzsA7ACALVu2sGSZ0qlGDkoF2ClS\nOXLh5K2baEGUzkHaw9/hkVe9+nQx+stvQwqzYfn95ldeU00pekZi2Lm/B43VZcZM2/YZPwu1A6dH\n0XV2HG3LRyO5ucoRPk3TvHq1Q/e9id7ChFbVeuH+EX4bpaxMPnBmFD3DMdy6rg4VpUW4e0OT5+Zv\nuzb9xkjVVzrdB7em4grqKPVimbZWIgDfBfAOY+wZAGCMHQawXHjmDIAtftZKs/MLyc9yPl1dUCoT\nE08RmVb0mUA30YLQFuYdboftFV46F/qLwyvWkRedYfqorbEKO/f3KCOhirCxFFPBz0JNDsMdFrxs\nVd7koNCJoMZiM2l2+0FoBfTrxe9GKD8rWw55WYcNjKnzUS8GZFqjRqZvDrcD+AyAw0R0KPHdE4yx\nV2wLmro8l/wsTzS+OKrLi0KdnBfbGkUF3WYVhLYw74h22LrFZXO6U23cYQL2ydApb/3oDNpHXARx\n94YmzwBwIqKKkKv6zitaqi09Kr2Pqb7EZs1xOXx7a71Vv5h48jfXlqOtsQrP7DmWYlqugjwHxP/H\nfRouobqsSDnGHd1DaREBwoqnbYwyojZGEZFpa6XfAt4xcxlja03KKi0u1DZejhLKkUsnW1MEZVBR\n6Eu8ZPNhoNu45UUQRkwlxtO30T8EhSha8kukJL8DhIuQK34XlqGK4hTV2IuOZ/dubvYVeZquOfHE\nTaA02bkp3To/DEBtWs7rNu2zju4hvHrkAohIOcY6XUAYcauKOQfRWYVF3gTee+a5bz0Vu+ZjysBZ\nuuBtYVNz5hN0Qdq8InvKv0UVvVQV0XX6yhwKCpBMEA8gJdLrvuNxqxJd/lw/7Ds+hMP9E9i6rs4z\n9aZpPm2vtJsDY9PJ4IKf3LTCmGa/HMA6WlTP8e95mlWbMRMjjnLa+eYil9PZM5oMcDc4eUX5jJyf\n2WTNvXL4HPYcHcTmVbXYtqHJKsgkT3t7XXMNOntGtX24rKoUl2bmUoLb8bptogYXFRI2NtcoA0mq\n9hjTQI+6sVWlNPUrz+uZp59++tx7Ok1oZVmRtfwzl62PTGEqerFV+KkC10V101LZla+qr8ChvvGU\n0xc/oXKZrspj1RSmtJs819E95Ou4xUUJNubLsnOYSQ4KOZ+BSIOfXsDPuEGV5U8uR3UrN7kV+Pmm\niHJ6W9k5z1Xw6tHzac5kcv3f2H5D2vs21l3cFNcGslhKbJfYL1xXqgoHYotMiJfyhjkUFxZYbxxe\n1+Zcg6mVkq11hxfTkAPXRSVKMpGR+z1vy9h1CzKIiW17q1201CDw26D8Nn+dwtVE/m9qDaP63sTX\nxMQ3RZbT28w90SBAZRbqV79JTnUOOfqsLgqtyfoFkEKXqSGBX3t8xzxgsp+8YQ5BIC6CxTC7tNnQ\nTEMRyAvZdmPSWYmMx2bx8Q3LPRXEQeElIzd9PgxjDzvWfqfFsGaWgHqD0ilaVW3W0WAi/w9zCNA5\n63n9lRGWucpmoeJpGYCv9ZBN/Vd1S/GQF/wvN/W2Xb/iZ52u1IZeVQBL2YmuoLSi2rehCuSNzmHH\njh1PPfroo1bviLK7xUj4YyPLlJ/V0SfLH8PqUfhkerNnDDetWZoin8+ljHlif+jk4Sbv+vWTl04m\nLLz6U04uI+t8xGQtNnQGkf+HbZNujur0KlHrAuXkRbsO9mFjcw3u/OBy5fM29fP+/Nj65WioLMHH\n1i/H9JU57Dl6Hm8PTCQTH5msX5u1LI63PI7ibx3dQ3jtxAjWNFQks/adGJpK0UO5ZD8+EM0PMyVa\n4uKI8dhsiuu77lnx72KZ0Xo5B+WSdZfYH1HclnSwvWXY3K5E8Yd4suXvq5LL6Gz+VeaNO1/vwemR\nKTxwQwtGYzNJmvxED2HmfpBTN6c5aqjCY+zu6tdGbg0ClfPazv09ODMSQ1tTVSTr1y8Mhrxvvdw1\ngB++0Yu9x4bwyB3XpOTZ4M8DV+ddULxvmAOQmckqD6ycqUmHsMwgCpt5E5lppukzeS4TjJP7XkxM\n2zkw2cwhfiARr/0AlGIfbvOvi7kj6yi4efDE5RlcmLiC8uIiT5qimvtBdANhFdg6qBTrXpFbo/Kr\nEQNsRrFWVGMjiode7hpIMUpgYLh0ZQ5Hz03g+IVJrQiXGz0EFSvlNXOwHdxMnIzlgTWpI4qNOEqb\neQCRWk7Y0Be0HVE4GskpTk1gO74qxT9/3+Z2JOso2luXoXcklnZz0IEzl7qKEuUtJhM3aVNGEnQO\n1FWUYHp2DnUVJco65T6Nyq8m6gCbqhtmc+3VGEuyUYKcvc+v3J1XYhNB6Mpr5mA7uJkQLfmJh1Sn\nFdmENIp6wyLqEAym9AVth40Xqa5eXcRTMRonT9DCyzXZ8HRKYU6XiVWQil7xb3NtOb76iQ/6tpND\nZf4JwGj9BA2+aIqgc+DVo+dxYnAKrx49j20bm9J+l/vUVnG+WJBP+oA6ban4vIl5LS/3ofdjDukg\ngxu1aMlvUatOK6qTZNT12sLUcgKIVhQUtB1+p0KTenU3pO/sO4lfHxtC7MocPrp+ufV8EUUCNu/5\n0RumHPF0Kpt/+o25V99GsZ6Cts0vf4oqfLhpWBUvZMp/Kio9pByuJijymjkE6bxMnhZszfyy5W9h\nG5BQRpQM1i/iqO47v1NhGFzTUInOnjFc01CZUq7pptBcW46J6Vk8f7APk9Nz+IrFCT8qePWXKBIx\nGT8vnwzbfteNrRx7y6Sv/cQ7mVKIZ6pcna8Or9N0zxBFprWCyM0WBYHfzANw++eBsenkd3wAMrEx\n80nT0T2krS+T9ZtCRacN4puFfx4DGarxUNFi+p0IuV9VdelokZ/99G1r8Nhdbfj0bWtSyrXpt1PD\nU5iYnsXe7kEtDZmEF61+fSND1HfIsJ3PurHddbAPuzr7kt/v7urHd/adwu6ufqNyZYhe2FyWH9U4\ntLcuCxSt1qbfxX6yXa/trcuw/eYWbE9YMQVFXt8c/BCEw3uFfvZDJpKuZAJhT9l+Dly6PlPpCVTm\nmyr6bGk28dLVea6qTnBtjVVWoZkfuLEFRwYmMDsfr8v2hBlWdMFFW70jl7Bj70ncu7lZaSZpQleU\ntzLd2Mr6H13OFpvTM/fClrO4ecHEoimouKej+2oiIb8Ag15rwA9cZMppdx7SCgTVSehCP/vhjdMj\nONQ7htbllVlNF+qHqOTX/K+pkl2lJ1CZb8qbcxBm7Tf2potP9pC9b/MKAPCNXzMam8Hy6jKsqbfX\nLSX7czi40QK3dnm+82xajCjbdRGlfktXlhweXs7ZEoahHekfT7Nq4vCLVBClCKm99WoiIb8Dg9xP\nQZnR+9aU1Q9BdRImeXs5xM1rYnouLe9uWGQ7eKDJyamjewjHL0yhuJA8b03yeyY3raDMWie/1bXB\ni6GNxWYwMT2H1uVxHcQPX+/B8wf7cKR/AhtXVivj7YTRLe3u6kfnmVFsaK5OYcC284DTLps82q4L\nU71QUDpVG3BYnZL4/jN7jmmtmvxM0aO8NTXXZjY5jww+/j8uKAy0z7+nmUMQeFmxqCBuXtu3tOAv\n29eFUs6pyg8iGouKoZjUL56ITIKZcZgEQLNl1n5taG+1N2XuHpxCz0gMD97cgubacpwansL49Ax+\n9e55vHNuAmsaKtKig/ptwF5jRCAUFxbihpalgUVBnAbbiKIqqOrW0ROETtGSKuqENYC3VZOfhVDU\nVoG68sKsWd27zbVxp1wqeY/fHGbnFzIycYBwAyNuXl5iDz8XeRNRjAmdUflRcJic7oOeiExOZbbM\n2qsenhlMJa7xiqopmx4/cse1uDB5BdMz81jTUIGHb13jmTReBK9HzM0tMiyOB29eiXs3N6e1wasO\n23nMb73j0zOoKS9J0UuIsNEBBdG76ez8RYQR73hZNXlt/gNj09jd1Q8C4ZZ1dcporFGho3vIKuGR\n/K6ub9pbl2Hh8uRoEJryhjlMXZ7LWJyWMBPPdPNS1WG6OfLn/XIA8Hqi8KPgMA1vHOSEFfWpzK+e\nnxzo1caV0s0BlXjoxtVL8a+fuillIza9LfF6xNzcsnJcDmfN2+B361F5uYsbnLz581svjzCqC/mi\nGifd2NmEw+ZQxUiSEaV4xxS7u/rxHx2nsaS0CMcHJ5W5I6JCe+sy7D02iM4zF/Fy14DVgcjrANlc\nW46F2PhIEJryhjnokv0EOfXL75hOvLA3DLkOkwXvV4bXM1GccIIuyjAyaT+7dyBY2AevvrE98ZqO\nnSofAN8Adb4w8mcOv0OMysu9o1ufuKi9NR6C4+i5cWxYURN4/utCjZuCt8sr2ZPpQSJKkSqBsKS0\nCBtWVFvdDoOgubYcm1fV4ui5CRzqu+gbuFN+1+YAaYqMMgciWgXg+wAaATAAOxhjzxLRUwAeAcAN\nd59gjL3iVZYu2Y9qwfhNEPkd04kX9oZhSr9tGbbP2CCM2V5QmbRKCb27qx+7Oq+eck2Tn+isoGTo\nTrxe9KpEATrrl6vWTunl2CjHxb9yXbqsbbrERc215VhdX4GuvnFUlxelMDHV2vEaU9FEMwrH1LDB\n+GQag0C0mLK5HdpAbOd9m1fG9VsG1kw6RHnDyvTNYQ7AY4yxN4moCsBBItqT+O2bjLFvhK1A1Rl+\nEyRoB2biarsYt5ZsQLfgvcJSi+962b379Zk8/mKsJC/Ft6pMr7qSiviRGF7uGkBNRXGKLkG0flGF\nrbAdU85E3uq9iGf2HENjVRm6+saTdenEP1wpLSbFkW8t47HZFCbGyzTtI1MTTa92iQi6yUe5RhdD\n7Cm3M6w1U5Q0E2PRBFozqozoZwD+DcDtAKZsmMOWLVtYZ2en0bNRbaS5tCFfzbO8ImMTVhaBRN3u\nMG2wGQv52a/sOoQ9Ry9g24bGtJzCtm2WGY2sYL6jtSHNZl8Hk/5QtZu35/YP1OOjbcuNRW1e9YUd\n+2yFfM9WeVFhMegiooOMsS3W7y0WcyCitQD2AdgE4MsA/gzAOIBOxG8XFxXvPAqAp3/bBODtxaCV\no6Cipr6grKpu4fLkaFCljoAGAMOB3y4sLi4oraheuBKbQMAoi37g7cXC3AwKikqkdgejX6QbQKbb\noKFhZVFVfdF8bHyIzUzHxJ982iyWUVxQWlFdUFpRQ8XlVWwmNjY3dr5H/t2qbQbvFFTU1APUBLDz\nnC4qKa8orKhZJrbHaK6a0JiZeRZu7keAkGtZTf8irMmIsJ4xZh22YVGYAxFVAtgL4B8ZYy8SUSPi\nnc0AfB3ACsbYn/uU0RmE++UKHP3ZQz7TDjj6s433K/0ZD7xHRMUAfgpgJ2PsRQBgjF1gjM0zxhYA\nfAfA1kzT4eDg4OBgjowyByIiAN8F8A5j7Bnh+xXCYw9gkcVFDg4ODg7eyLS10u0APgPgMBEdSnz3\nBIBPEdENiIuVzgD4K4OydmSEwsWDoz97yGfaAUd/tvG+pH9RrZUcHBwcHPID7+lkPw4ODg4OweCY\ng4ODg4NDGnKWORDRdiI6QkQLRKQ1wyKiTxDRMSI6QUSPLyaNXiCiOiLaQ0Tdib9K11wiOkNEh4no\nEBGZefllCH59SXE8l/j990R0Uzbo1MGA/juJaDzR14eI6GvZoFMFIvpPIhokIqVxRh70vR/9Odv3\nQDzUDxH9moiOJvadLyieydkxMKTfbgwYYzn5D8CHAKwH8BsAWzTPFAI4CeAaACUAugBsyDbtCdr+\nGcDjic+PA/gnzXNnADTkAL2+fQngHgC/AEAAbgOwP9t0W9J/J4CXs02rhv47ANwE4G3N7znb94b0\n52zfJ+hbAeCmxOcqAMfzbP6b0G81Bjl7c2CMvcMYO+bz2FYAJxhjpxhjMwB+DOD+zFNnhPsBfC/x\n+XsA/jCLtJjApC/vB/B9FsfrAGols+RsIpfngi8YY/sAeMXdz+W+N6E/p8EYO8cYezPxeRLAOwBW\nSo/l7BgY0m+FnGUOhlgJ4Kzw/z6E7JAI0cgY40H6zyMemVYFBuCXRHQwES4kWzDpy1zub1PaPpwQ\nCfyCiDYuDmmRIJf73hR50feJUD83Atgv/ZQXY+BBP2AxBlnN50BEvwTQpPjp7xhjP1tsemzhRb/4\nH8YYIyKdzfBHGGP9RLQcwB4iejdxCnOIHm8CWM0YmyKiewD8D4DWLNP0fkFe9H0i1M9PAXyRMTaR\nbXps4UO/1RhklTkwxj4esoh+AKuE/7ckvlsUeNFPRBeIaAVj7Fzi6jmoKaM/8XeQiP4bcfFINpiD\nSV9mtb994EubuFgYY68Q0beIqIExltWgcIbI5b73RT70vSrUj4ScHgM/+m3HIN/FSgcAtBLROiIq\nAfAQgJeyTBPHSwA+m/j8WQBpNyEiWkLxPBcgoiUA7kL2QomY9OVLAP40YbVxG4BxQXSWbfjST0RN\nRESJz1sRn/9ho+0uFnK5732R632foC0t1I+EnB0DE/qtxyDbWnYP7fsDiMv0rgC4AODVxPfNAF4R\nnrsHcc38ScTFUVmnPUFXPYBfAegG8EsAdTL9iFvWdCX+Hck2/aq+BPA5AJ9LfCYA/574/TA0VmQ5\nTP/nE/3cBeB1AB/ONs0C7T8CcA7AbGLe/0We9b0f/Tnb9wn6PoK4/u/3AA4l/t2TL2NgSL/VGLjw\nGQ4ODg4Oach3sZKDg4ODQwbgmIODg4ODQxocc3BwcHBwSINjDg4ODg4OaXDMwcHBwcEhDY45ODg4\nODikwTEHBwcHB4c0OObg4LAIIKIPE9E/ZJsOBwdTOCc4BwcHB4c0uJuDg4MHiOhDRLQvEeb4q0R0\nImA5u4ioPWr6HBwyBcccHBw0IKIiADsBfIExdj3isbCCBkbchHjcGweHvEBWQ3Y7OOQ4/ghAF2Ps\nrcT/j0IKvW6Sk4SIygCUMMbGM0msg0OUcMzBwUGP6xGPbsmxCcD/ig8ws5wkGxFnLA4OeQMnVnJw\n0GMEQBsAENENAP4E8XDHtrgOTqTkkGdwNwcHBz1+AODnRHQYwG8AnGGMnQpQznUA3oiSMAeHTMOZ\nsjo4aEBElYyxqcTnrwKoYYw9mWWyHBwWBU6s5OCgx5eI6AgRHQKwFsDXs0yPg8Oiwd0cHBwcHBzS\n4G4ODg4ODg5pcMzBwcHBwSENjjk4ODg4OKTBMQcHBwcHhzQ45uDg4ODgkAbHHBwcHBwc0uCYg4OD\ng4NDGv4fsQOuVDVaaecAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1170ae080>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "mag = cmd_mags['r']\n",
    "color = cmd_mags['g'] - cmd_mags['i']\n",
    "\n",
    "msk = (tab['phot_sg']=='STAR')&(tab['sep'] > inner_radius) #arcmi\n",
    "plt.scatter(color[msk], mag[msk], s=2, alpha=.5)\n",
    "\n",
    "plt.xlim(-1, 2.5)\n",
    "plt.ylim(25, 17)\n",
    "plt.xlabel('$g-i$')\n",
    "plt.ylabel('$r$')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [default]",
   "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.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}