{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Getting Started" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "SORA version: 1.0dev\n" ] } ], "source": [ "## SORA package\n", "from sora import Occultation, Body, Star, LightCurve, Observer\n", "from sora.prediction import prediction\n", "from sora.extra import draw_ellipse\n", "\n", "## Other main packages\n", "from astropy.time import Time\n", "import astropy.units as u\n", "\n", "## Usual packages\n", "import numpy as np\n", "import matplotlib.pylab as pl\n", "import os" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Before analysing stellar occultations data, let's predict them.**\n", "\n", "To predict stellar occultation we needs the intended Solar System body ephemeris and a time window." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Obtaining data for Chariklo from SBDB\n", "###############################################################################\n", " 10199 Chariklo (1997 CU26) \n", "###############################################################################\n", "Object Orbital Class: Centaur\n", "Spectral Type:\n", " SMASS: D [Reference: EAR-A-5-DDR-TAXONOMY-V4.0]\n", " Relatively featureless spectrum with very steep red slope.\n", "Discovered 1997-Feb-15 by Spacewatch at Kitt Peak\n", "\n", "Physical parameters:\n", "Diameter:\n", " 302 +/- 30 km\n", " Reference: Earth, Moon, and Planets, v. 89, Issue 1, p. 117-134 (2002), \n", "Rotation:\n", " 7.004 +/- 0 h\n", " Reference: LCDB (Rev. 2021-June); Warner et al., 2009, [Result based on less than full coverage, so that the period may be wrong by 30 percent or so.] REFERENCE LIST:[Fornasier, S.; Lazzaro, D.; Alvarez-Candal, A.; Snodgrass, C.; et al. (2014) Astron. Astrophys. 568, L11.], [Leiva, R.; Sicardy, B.; Camargo, J.I.B.; Desmars, J.; et al. (2017) Astron. J. 154, A159.]\n", "Absolute Magnitude:\n", " 6.58 +/- 0 mag\n", " Reference: MPO647128, \n", "Albedo:\n", " 0.045 +/- 0.01 \n", " Reference: Earth, Moon, and Planets, v. 89, Issue 1, p. 117-134 (2002), \n", "\n", "Ellipsoid: 151.0 x 151.0 x 151.0\n", "\n", "----------- Ephemeris -----------\n", "\n", "EphemKernel: CHARIKLO/DE438_SMALL (SPKID=2010199)\n", "Ephem Error: RA*cosDEC: 0.000 arcsec; DEC: 0.000 arcsec\n", "Offset applied: RA*cosDEC: 0.0000 arcsec; DEC: 0.0000 arcsec\n", "\n", "\n" ] } ], "source": [ "# First, let's consider an Solar System Body\n", "\n", "chariklo = Body(name='Chariklo', \n", " ephem=['guidelines/input/bsp/Chariklo.bsp', 'guidelines/input/bsp/de438_small.bsp'])\n", "\n", "print(chariklo)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Ephemeris was split in 1 parts for better search of stars\n", "\n", "Searching occultations in part 1/1\n", "Generating Ephemeris between 2017-06-20 00:00:00.000 and 2017-06-26 23:59:00.000 ...\n", "Downloading stars ...\n", " 5 GaiaDR3 stars downloaded\n", "Identifying occultations ...\n", "\n", "2 occultations found.\n" ] }, { "data": { "text/html": [ "
| Epoch | ICRS Star Coord at Epoch | Geocentric Object Position | C/A | P/A | Vel | Dist | G | long | loct | M-G-T | S-G-T | GaiaDR3 Source ID |
|---|---|---|---|---|---|---|---|---|---|---|---|---|
| arcsec | deg | km / s | AU | mag | deg | hh:mm | deg | deg | ||||
| object | object | object | float64 | float64 | float64 | float64 | float64 | float64 | str5 | float64 | float64 | str19 |
| 2017-06-21 09:57:43.440 | 18 55 36.17454 -31 31 19.03261 | 18 55 36.17500 -31 31 19.60516 | 0.573 | 179.41 | -21.84 | 14.663 | 15.254 | 225 | 00:56 | 128 | 165 | 6760228702284187264 |
| 2017-06-22 21:18:48.200 | 18 55 15.65251 -31 31 21.67062 | 18 55 15.65249 -31 31 21.62190 | 0.049 | 359.72 | -22.00 | 14.659 | 14.224 | 53 | 00:50 | 149 | 166 | 6760223758801661440 |
"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Now, let's start instantiating the Occultation**\n",
"\n",
"An occultation is defined by the occulting body, the occulted star, and the time of the occultation "
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\u001b[0;31mInit signature:\u001b[0m\n",
"\u001b[0mOccultation\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\u001b[0m\n",
"\u001b[0;34m\u001b[0m \u001b[0mstar\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\n",
"\u001b[0;34m\u001b[0m \u001b[0mbody\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\n",
"\u001b[0;34m\u001b[0m \u001b[0mephem\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\n",
"\u001b[0;34m\u001b[0m \u001b[0mtime\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\n",
"\u001b[0;34m\u001b[0m \u001b[0mreference_center\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'geocenter'\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\n",
"\u001b[0;34m\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
"\u001b[0;31mDocstring:\u001b[0m \n",
"Instantiates the Occultation object and performs the reduction of the\n",
"occultation.\n",
"\n",
"Attributes\n",
"----------\n",
"star : `sora.Star`, `str`, required\n",
" the coordinate of the star in the same reference frame as the ephemeris.\n",
" It must be a Star object or a string with the coordinates of the object\n",
" to search on Vizier.\n",
"\n",
"body : `sora.Body`, `str`\n",
" Object that will occult the star. It must be a Body object or its name\n",
" to search in the Small Body Database.\n",
"\n",
"ephem : `sora.Ephem`, `list`\n",
" Object ephemeris. It must be an Ephemeris object or a list.\n",
"\n",
"time : `str`, `astropy.time.Time`, required\n",
" Reference time of the occultation. Time does not need to be exact, but\n",
" needs to be within approximately 50 minutes of the occultation closest\n",
" approach to calculate occultation parameters.\n",
"\n",
"reference_center : `str`, `sora.Observer`, `sora.Spacecraft`\n",
" A SORA observer object or a string 'geocenter'.\n",
" The occultation parameters will be calculated in respect\n",
" to this reference as center of projection.\n",
"\n",
"\n",
"Important\n",
"---------\n",
"When instantiating with \"body\" and \"ephem\", the user may define the\n",
"Occultation in 3 ways:\n",
"\n",
"1. With `body` and `ephem`.\n",
"\n",
"2. With only \"body\". In this case, the \"body\" parameter must be a Body\n",
"object and have an ephemeris associated (see Body documentation).\n",
"\n",
"3. With only `ephem`. In this case, the `ephem` parameter must be one of the\n",
"Ephem Classes and have a name (see Ephem documentation) to search for the\n",
"body in the Small Body Database.\n",
"\u001b[0;31mFile:\u001b[0m ~/miniconda3/envs/sora-develop-stats/lib/python3.9/site-packages/sora_astro-1.0.dev0-py3.9.egg/sora/occultation/core.py\n",
"\u001b[0;31mType:\u001b[0m type\n",
"\u001b[0;31mSubclasses:\u001b[0m \n"
]
}
],
"source": [
"Occultation?"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"1 GaiaDR3 star found band={'G': 14.223702}\n",
"star coordinate at J2016.0: RA=18h55m15.65210s +/- 0.0197 mas, DEC=-31d31m21.6676s +/- 0.018 mas\n",
"\n",
"Downloading star parameters from I/297/out\n",
"GaiaDR3 star Source ID: 6760223758801661440\n",
"ICRS star coordinate at J2016.0:\n",
"RA=18h55m15.65210s +/- 0.0197 mas, DEC=-31d31m21.6676s +/- 0.0180 mas\n",
"pmRA=3.556 +/- 0.025 mas/yr, pmDEC=-2.050 +/- 0.020 mas/yr\n",
"GaiaDR3 Proper motion corrected as suggested by Cantat-Gaudin & Brandt (2021) \n",
"Plx=0.2121 +/- 0.0228 mas, Rad. Vel.=-40.49 +/- 3.73 km/s \n",
"\n",
"Magnitudes: G: 14.224, B: 14.320, V: 13.530, R: 14.180, J: 12.395, H: 11.781,\n",
" K: 11.627\n",
"\n",
"Apparent diameter from Kervella et. al (2004):\n",
" V: 0.0216 mas, B: 0.0216 mas\n",
"Apparent diameter from van Belle (1999):\n",
" sg: B: 0.0238 mas, V: 0.0244 mas\n",
" ms: B: 0.0261 mas, V: 0.0198 mas\n",
" vs: B: 0.0350 mas, V: 0.0315 mas\n"
]
}
],
"source": [
"star_occ = Star(coord='18 55 15.65250 -31 31 21.67051')\n",
"#star_occ = Star(code='6760223758801661440')\n",
"\n",
"print(star_occ)"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Stellar occultation of star GaiaDR3 6760223758801661440 by 10199 Chariklo (1997 CU26).\n",
"\n",
"Geocentric Closest Approach: 0.049 arcsec\n",
"Instant of CA: 2017-06-22 21:18:48.200\n",
"Position Angle: 359.72 deg\n",
"Geocentric shadow velocity: -22.00 km / s\n",
"Sun-Geocenter-Target angle: 166.42 deg\n",
"Moon-Geocenter-Target angle: 149.11 deg\n",
"\n",
"\n",
"No observations reported\n",
"\n",
"###############################################################################\n",
" STAR \n",
"###############################################################################\n",
"GaiaDR3 star Source ID: 6760223758801661440\n",
"ICRS star coordinate at J2016.0:\n",
"RA=18h55m15.65210s +/- 0.0197 mas, DEC=-31d31m21.6676s +/- 0.0180 mas\n",
"pmRA=3.556 +/- 0.025 mas/yr, pmDEC=-2.050 +/- 0.020 mas/yr\n",
"GaiaDR3 Proper motion corrected as suggested by Cantat-Gaudin & Brandt (2021) \n",
"Plx=0.2121 +/- 0.0228 mas, Rad. Vel.=-40.49 +/- 3.73 km/s \n",
"\n",
"Magnitudes: G: 14.224, B: 14.320, V: 13.530, R: 14.180, J: 12.395, H: 11.781,\n",
" K: 11.627\n",
"\n",
"Apparent diameter from Kervella et. al (2004):\n",
" V: 0.0216 mas, B: 0.0216 mas\n",
"Apparent diameter from van Belle (1999):\n",
" sg: B: 0.0238 mas, V: 0.0244 mas\n",
" ms: B: 0.0261 mas, V: 0.0198 mas\n",
" vs: B: 0.0350 mas, V: 0.0315 mas\n",
"\n",
"Geocentric star coordinate at occultation Epoch (2017-06-22 21:18:48.200):\n",
"RA=18h55m15.65251s +/- 0.0323 mas, DEC=-31d31m21.6706s +/- 0.0341 mas\n",
"\n",
"###############################################################################\n",
" 10199 Chariklo (1997 CU26) \n",
"###############################################################################\n",
"Object Orbital Class: Centaur\n",
"Spectral Type:\n",
" SMASS: D [Reference: EAR-A-5-DDR-TAXONOMY-V4.0]\n",
" Relatively featureless spectrum with very steep red slope.\n",
"Discovered 1997-Feb-15 by Spacewatch at Kitt Peak\n",
"\n",
"Physical parameters:\n",
"Diameter:\n",
" 302 +/- 30 km\n",
" Reference: Earth, Moon, and Planets, v. 89, Issue 1, p. 117-134 (2002), \n",
"Rotation:\n",
" 7.004 +/- 0 h\n",
" Reference: LCDB (Rev. 2021-June); Warner et al., 2009, [Result based on less than full coverage, so that the period may be wrong by 30 percent or so.] REFERENCE LIST:[Fornasier, S.; Lazzaro, D.; Alvarez-Candal, A.; Snodgrass, C.; et al. (2014) Astron. Astrophys. 568, L11.], [Leiva, R.; Sicardy, B.; Camargo, J.I.B.; Desmars, J.; et al. (2017) Astron. J. 154, A159.]\n",
"Absolute Magnitude:\n",
" 6.58 +/- 0 mag\n",
" Reference: MPO647128, \n",
"Albedo:\n",
" 0.045 +/- 0.01 \n",
" Reference: Earth, Moon, and Planets, v. 89, Issue 1, p. 117-134 (2002), \n",
"\n",
"Ellipsoid: 151.0 x 151.0 x 151.0\n",
"\n",
"----------- Ephemeris -----------\n",
"\n",
"EphemKernel: CHARIKLO/DE438_SMALL (SPKID=2010199)\n",
"Ephem Error: RA*cosDEC: 0.000 arcsec; DEC: 0.000 arcsec\n",
"Offset applied: RA*cosDEC: 0.0000 arcsec; DEC: 0.0000 arcsec\n",
"\n",
"\n",
"\n"
]
}
],
"source": [
"occ = Occultation(star=star_occ, body=chariklo, time='2017-06-22 21:18')\n",
"\n",
"print(occ)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**After that, we instantiate the observers and their light curves**"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Observers**\n",
"\n",
"Now let's define our observers, they can be setted manually or from the MPC database"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Site: Tivoli\n",
"Geodetic coordinates: Lon: 18d01m01.24s, Lat: -23d27m40.19s, height: 1.344 km\n",
"\n",
"\n",
"Site: Observatorio Pico dos Dias\n",
"Geodetic coordinates: Lon: -45d34m57.54s, Lat: -22d32m07.74756091s, height: 1.811 km\n"
]
}
],
"source": [
"### User\n",
"\n",
"out = Observer(name='Outeniqua' ,lon='+16 49 17.710', lat='-21 17 58.170', height =1416)\n",
"ond = Observer(name='Onduruquea' ,lon='+15 59 33.750', lat='-21 36 26.040', height =1220)\n",
"tiv = Observer(name='Tivoli' ,lon='+18 01 01.240', lat='-23 27 40.190', height =1344)\n",
"whc = Observer(name='Windhoek' ,lon='+17 06 31.900', lat='-22 41 55.160', height =1902)\n",
"hak = Observer(name='Hakos' ,lon='+16 21 41.320', lat='-23 14 11.040', height =1843)\n",
"\n",
"print(tiv)\n",
"\n",
"print('\\n')\n",
"\n",
"### MPC Database Search\n",
"\n",
"opd = Observer(name='Observatorio Pico dos Dias',code='874')\n",
"\n",
"print(opd)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Light Curves**\n",
"\n",
"Now let's define our light curves, they can be instanciated from different way: \n",
"- **(i)** Manually with arrays containing the flux and the times; \n",
"- **(ii)** Read an ASCII file; \n",
"- **(iii)** Already obtained times."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Outeniqua (Namibia)**"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Light curve name: Outeniqua lc\n",
"Initial time: 2017-06-22 21:20:00.056 UTC\n",
"End time: 2017-06-22 21:23:19.958 UTC\n",
"Duration: 3.332 minutes\n",
"Time offset: 0.000 seconds\n",
"\n",
"Exposure time: 0.1000 seconds\n",
"Cycle time: 0.1002 seconds\n",
"Num. data points: 2000\n",
"\n",
"\n",
"There is no occultation associated with this light curve.\n",
"\n",
"Object LightCurve model was not fitted.\n",
"\n",
"Immersion and emersion times were not fitted or instantiated.\n",
"\n",
"\n"
]
}
],
"source": [
"out_lc = LightCurve(name='Outeniqua lc',file='guidelines/input/lightcurves/lc_example.dat',\n",
" exptime=0.100, usecols=[0,1])\n",
"\n",
"print(out_lc)\n"
]
},
{
"cell_type": "code",
"execution_count": 10,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZMAAAEPCAYAAACHuClZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABsrUlEQVR4nO2debgUxfX3v+eugIooIKuAe0RFVAQHBQeviOC+L4m4RFF+LsHXJaIYMCBEjYoJiqKicY0aN6KikisDCIPKjoiKyBowwkXBBe425/2jp5qanuptume5l/o8z31uT3d1dXV3dZ2qc06dImaGRqPRaDRBKMp3ATQajUbT8NHCRKPRaDSB0cJEo9FoNIHRwkSj0Wg0gdHCRKPRaDSBKcl3AfJFq1atuEuXLvkuhkaj0TQY5s+fv5mZW6uO7bLCpEuXLpg3b16+i6HRaDQNBiJaY3dMq7k0Go1GExgtTDQajUYTGC1MNBqNRhMYLUw0Go1GExgtTDQajUYTGC1MNBqNRhMYLUw8EI/HMW7cOMTj8XwXRaPRaAqSXXaeiVfi8Tj69euH2tpalJeXo7KyEpFIJN/F0mg0moJCj0xciMViqK6uRiKRQE1NDWKxWL6LpNFoNAWHFiYuRKNRc7usrCzlt0aj0WgMtDBxQVZpaRWXRqPRqNHCxAdakGg0Go0aLUw0Go2mgGio3qMF7c1FRPsCeA5AWwAJAJOY+RFLGgLwCIBBAH4FcAUzL8h1WTUajSYo8XgcJ510Empra1FWVtagVOuFPjKpA3ALMx8K4DgA1xNRV0uagQAOSv4NATAxt0XUaDSacIjFYtixYwfq6+sbnPdoQQsTZt4oRhnM/BOA5QA6WJKdBeA5NpgLoAURtctxUTUajSYwDdl7tKCFiQwRdQFwFIBPLIc6AFgn/V6PdIGj0Wg0BU9D9h5tEMKEiHYH8DqAYcy8zXpYcQrb5DOEiOYR0bxNmzaFXUyNRqMJjYYkSIAGIEyIqBSGIHmRmd9QJFkPYF/pd0cAG1R5MfMkZu7BzD1at1YuY6zRaDSaDChoYZL01HoawHJmfsgm2RQAg8ngOABbmXljzgoZMg3VLVCj0ezaFLRrMIDjAVwGYCkRLUruuxNAJwBg5scBvAfDLfgbGK7BV+a+mOEggkrW1dU1OLfAbBKPxxGLxRCNRvXz0GgKlIIWJsz8MdQ2ETkNA7g+NyXKLk899RSqq6sBwHQL3NUbz4bsd6/RNHT8dOQKWpjsanTp0sXcbmhugdlC+N0DWsBqNLlE1ZFzoqBtJrsa3bt3BwAccsghugeepCH73Ws0DRm/Eyi1MCkgmjVrBgBo3769FiRJGrLfvUbTkPHbkdPCpIAQwmT79u15LklhogWJRpM7/HbktDApIEpLSwEAq1ev1q7BWUS7XwdDP79dDy8dOW2ALyAWL14MAPjuu+9QUVGh1TpZQHuHBUM8v5qaGpSXl+vnpzHRI5MCYsGCnZHzG1rE0IbC9OnTG2xU1kJAGGUTiUTBPj89csoPemRSQBx99NHmdnFxccF7LjXEyYTHH3+8ua29w/xT6N518Xgcffr0QX19PZo2bapHTjlEC5MCwph/aWBEkilchLpjx44daNKkCT766KPAH20uhJNwvy4pKQnU0DREQRoGhe5dF4vFUF9fD2DXm5eU7zqphUkW8ftyZTVXXV1dQX8IsVjMnK1fXV0duKxvvvkmzjvvPBARSkpKcOWVV+Lyyy8P/f5FmXfbbbdAgkTbXQrTu67QR07ZohDqpLaZZAnxckeMGIGKigpP+tujjjrK3C4uLsbatWs9630//vhjjB07NiV9NnXH0WgUxcXFAIxeftCPdsKECWBmUxf/xBNPeH5ufhDCpLy8POM8Xn/99UZpd2kMtoZCGTnl+lkWxAqNzLxL/h1zzDHsFRjro3hOz8w8duxY87zi4mIeO3as6zlz5swxzykrK+OioiJu2rQpz5kzx9N5RGSmnzNnDpeXlzMAT3lkwmWXXcYA+Nprrw2c1+WXX27eu/W5ZfL87VixYgUD4I4dO3pKP2fOHB47dmzK83v00UfNMmXr2eaaOXPmcHFxcUodsiPM95ENRPlU7y4XzJkzh5s0aeLpWarI5PnKbUeYddJaFgDz2KZN1WquLBF0uF1TU2P+d1MhTZ8+HYDRMZB7JWGqoVR07NgRANCpUyfXtG4qv/bt2wMAOnfujDVr1gDIjppCPBMxp8cJEcW5trY2xQ32sMMOA2DEUnvppZcKUt3jl8Zoa+jbty/q6+vRpEmTnI5S8hFPrhBGZFqYZIlMXi5LBniBlwb1xBNPVKYvKipCIpFAaWlpXnXH8Xgc0WgUdXV1tnMTREN2zDHHYM2aNejQoQNee+21rNlMysrKXNPKdiG5URDOEfvuu2+Db3AB4/2sXbvW/N1YbA11dXUAci8c8223Ces+/arotM0kB3h9uSph4iSIhF42kUikpY9EIqYN5v77789roxeLxVBTU+M4N0Hcg7DDtG3bNitlnjdvHgCgtrbWNW2+G4VcIAT9k08+ae5rbA4FuX538rObNm1ag3yWol7Iv93QI5MCx0mQRKNR1NbWokmTJsr0LVq0AAB07drVNv9cuBN6aZTFyEQIk6Co7isej+Omm24CAKxatQrxeNzxnt1Glyrhn08yeZdC0Mtku/GbNm0aPvvsM/Tr1y/QtWKxGD766CMMHDjQMZ9MhGNY30XPnj0zPteJbH+31nrhZWSnhUkBoWqcmFk550R+2dbGQFS0rVu3Ol4vHo+joqICNTU1WXUn9KLyEyOTkhKjSgZpqO1sHbFYzByRMLMv1YecrhDnAGW6SmeuR1tTpkzBWWedBSIKZMsQ3pLMjL/+9a+O+WQiSMJys62vr/dkn/NbvmxPzLTWCy/1RKu5coDcMM6ZMyfNhVeVzmkfkN7bF8TjcZx44om48847zXkrdnnEYjFs3749p+6EdpXeOjIJIkyErcOqVotGo+aHXVRU1KjUVuKe/b7LXKtgpk6dCiDdWcQvsVjMrCNh190w3WxFvQ4TlbNE2FjrhZd6ooVJDonH4zjhhBNw1113eZ5DYdeoyi/3vffeM7fl3rdsS1FRSDYBq80kCHb3FYlEMG7cOADGAmQNUZdtRyG9SyfkkEFBymm93xYtWth20oLmHeRZehEmfuekFOq71mquHCBUVarelNyg+RmZyBx77LHmtlyxhDeXNQ9Z3yrIt9E1zJGJk1rt0EMPBQDsscceGecvKCSbSaauobm+h27dugEw3MpfffXVjOucfN748eNx7bXXAgCaNm0auIxhutmKem1n4xAdTGY21X65LF+YaGESAlOnTsWiRYtcjWFuPYpMP2z5PPn6hx12GJYuXZqSds6cOTjxxBNRX1+fMgrIZ4V8//33MX/+fADhCBMZ64f7z3/+M3CehWgzkfHzLt1Gr0FwMhKHuZpoVVWVuW21HwYlaBnr6+tN1XNdXV2anSgWi5nvIBOVVaEIEkALk8A8+eSTGDJkCIqKimznUIiGMay5J1asDUJJSQnq6uqw2267paX96KOPTP978T9Twmjw33zzTZx77rnm702bNoWWt4xwNhCrWP7000+h5t9QyZYwsXOCEIT5fq2dtEJaqbS+vj5F9WzVSARVWc2aNQt9+vQJq7i2jBs3zrVsoQsTIiph5mCtVJ7x43b37rvvAkCasdeth6HKN1M1lzVNcXEx6urqzCG2fLxv376u+fklSE9dGGQFGzdutE3r5srrhHA2EGzbts13HqJe7LXXXgAKS82VKdkwEAPuEz7DfHbWTlrv3r1Dyzso9fX1jgIjkw6mbFsZMGBATlRdI0aMEPHs0nuoSTwLEyKaBOAmZt7hkGY/AC8DOM5HOQsKVY9KIKSzPG/hhx9+AGA0qGVlZWjZsqXpoinw6vKbqTCx9i6FqkjVUPTq1UuZh50LcrYROnRBp06dMGfOHPO3/OEEWX0yGo2aNiQAaN68ua/zZddb4b7cGMjWyMSuAc12HQujUfXyTNw6nEQEZkZ9fb1ngeFnbpAgyOx+P51m0VkGYGts9PNVXA3gOCK6kJm/tB4kovMBPAnA31daYKh6VAIhnYWAEb7egOEdNHny5JTznZgxY4bnIa3bS7cKnKIiw0lPpcayE06ZCJMwepfWCZX77rtvSt5hfTiRSATHHnssPvnkEwDA7rvv7uv8f//73+Z7bQwjEkG2hIm1AQWMztg+++wDoLCfoVt0BGEDqa+vt1Vty8JEJgxhJ7cbmYZKEvdQW1vrea5KUoVoqx/2I0zuBTAcwDwiupGZnwEAIioDMB7AtQB+AHCOjzwLDqchqVWVJVcUp/UxxIcj97IHDhyYtqCU6gMTPWIAKS9dzkul5gJ2ChP5uJMwyZQgvU27sov9LVu2NI9lolOWVWNBPLh69OhhbpeWlnrqMDQEwhQmotPTsmXLFKM4AHMSoBjVbdy4MeXdvPvuu1i4cCEqKiryblR2M+LLNpDq6mqMGjUKo0aNUk5stQoTq3YjE+RzX3nllYzycrLjqOjatSueeuop9O7d+xe7NJ6FCTPfTUQxAC8AeIqITgLwCICnAHQDMBvAJcy83muehYjbkNSuQWvbtm3a+VbkXnZtba0n1+CZM2ea29YJeALR2xY4qbnsGo9EIhFaKBM/WMvz3XffATCexZw5c3DDDTeYx8aPH+9bpyyrxtxUj04I9+vmzZvj/vvvx3XXXVcwveuPP/44pZ648dFHH5mzvH/zm9+4pveiDpF7ulamTp1qRtEVz2zDhg046aSTTIeQ008/HQAwduzYvLu7uo1MrB3MadOmYdasWSnlLioqQn19fdo3eNddd4UaxViet2PF6b35NfwffPDB4YZTYeZKIjoSwPMALk3+JQCMATCKmbPnZ5gHVA9PNTIAgNatW9vmIz6gTIansqeGeOkfffRRSu9p9uzZKedYRyaqsnjdn+0YQNbrvvDCC+b2008/nfJhW3u7dtipxrwEd7RDPNNmzZrhiCOOyDifsInFYujXr5/n0aHwagOMke4bb7zhmH727Nno168fEomEGVrErhx2z1du8ORyincj1+VCCH8vl0fl9KHy1rSW225kokrrB/E9Cuw6h6KzUFNTo1TF+TX8e+k4ZWJJ/BnAJgCiVmwFMDNbgoSIJgM4HcD3zHy44ngUwNsAViV3vcHMf87kWtYXpSIWi2Hp0qX4wx/+kLLfy8OWX5gqvLoqDzlQnHjp1o9WzmfcuHFmBfYjTFSV0s29Mwys5ZE/vgMOOCDlmFcVl12vK4grtNNoL59MmzYNgPeRllXQfvzxx47pJ0yYYNa3HTt24LnnnlOmc3o3YlS3xx574KabbsK9994LwHim0WgUv/76q5m2EGZ0P/vss+a2V6cPa7nthIkqrVfE9yirWBOJhLLdEiFhAHcBfdxx7v5SoQuT5KjkFQAHAfgAwJsA/grgfSK6D8DdWRAqzwKYAEBdiw1mMfPpQS4iB3ez7pe56667TNdbGSfds+pFHHPMMb7LKCqDNRKpnNddd91lbgcdmdi5d4aJnVszM+Oggw5KOZZJUEa5IQgiTEQ55cahENRcTs9ENaq0rn3j9kz3228/c5uZ8cwzz9iWo7i4GPX19ejevTsWLVpkHhMN62677YZzzz3XFCa33347IpEIfv75ZwDGaD3fKq54PJ7yDXldWM5abidhkuk9WkdxAPDZZ5/hwgsvTBuZ+lFj1dfXp3goqoRTqMKEiK4H8EDynDuZ+b7k/ukwBMwdAKJEdAkzr7XPyR/MPJOIuoSVnx2vvPKKKckFQsBYyqMUHH4NmaqX8+mnn7qmmz17Nv7zn//YppG3VY2nk81EIBtSBV57U37VYtby9O/fH1OnTsUvv/ySsXHYLiJAEJuJKEt9fX1BzYAXa9bsscceKRMx33jjDZx33nlpk2ll1/Dx48c76twB4MADD0z5rapT1sZnv/32SxEmwrvQ+j4PPvhgADvfhRfh5oYqvpWsqnKrn/KMdGDn6MkNu3Kr6nCm9xiNRtM6suPHjweQXp+tHaoff/wRY8eOVYb9l4WJUINaHUy8fIt+RiZ/B7AWhpHdfGPMvIKIjgPwIIDrASwE0FKdRdaIENFiABsA3MrMy1SJiGgIgCFA+lKznTt3TksvDxVlysvL02bZOjVOXuaPWHtEqnTxeBx9+/ZNe7HW38ItMRNvLrFOSk1NTUqcIy+9KbkielWLWcsjJjGuXbsWX36Z5oHuCbuKH/bIpBAQ5dlzzz1ThMnbb78NwHgWsseR/AyGDRuGNm3aOOZvFZzWGeazZs0yQ8GLstgJ20QikXJMbIv3lamQtjpcWBEjCy/10yo4JkyY4Lvxj8fj5rMIs75EIhGcddZZeP3118198pwsO1auXInLLrvMNuz/nDlzMHfuXESjUbPNy2TU7Sdq8NsAjpIFiYCZa5j5RgDnpp+WdRYA6MzMR8IQeG/ZJWTmSczcg5l7WA3mYr6D3BPbe++9lfmojJCZ9nQFdrYaOV9rr8nu2mI9dVVFdrOZyENpeUjtZX6LmGUuN2BOkVDj8TheeeUV5TFmxrJly9LSe4muavcBhyFMgoagCcqjjz6Kfv36YdKkSQDsoy0ffvhO86LwOKqoqMDzzz9v7q+pqXFtjOQGvkuXLml1/6mnnkqJtmA9Ry6jVZisWLEC48aNw+TJk1PS+WX69OnmtsqtVwgIuX7axcGy1vMhQ4ak1DuxbUXujAk7I4CUEVoYtGvXTrnfSRCLKB2y8V9m4MCBuPvuu1FRUYGWLVsq8wpVzcXMrvNHmPktIprvNc8wYOZt0vZ7RPQYEbVi5s2Z5CcLkJtvvtns5ct4DYUiY1UFzJs3Dy+88II53LYbSsv5/vjjj8o01o9QxOQShk0/I5NoNGrec2lpqdlIOPnHi8onn5tIJPDhhx+aLpNW5JUiVRBRmtvqCSecgEQi4TrJKpvCRH7WzIzZs2ejsrIS/fv3z7quf9KkSaartKhLoicuVEkC8exat26NTZs2gZlRXV2d8v7LyspsIyII5Hw7deqUdo/du3dPM8o7CROZ0aNHp/z+9ddfMwqZI9uB5DorUMXB8upN+eqrr+Kiiy5CUVGRqQpS1S8xCdk6afmdd97B1Vdf7ek+vDgA2QkNVTslEB1lEaXDet+yXbSqqgpHHHEEFi9enJJmy5YtruUPfT0TZl4Xdp5OEFFbSj5hIuoJ4568+ZAqkB9adXV12kdqh5MwmTt3LioqKjBixAhz34UXXog777zTXNfEywd0//33e7q2XDmAnaqjeDyOW265xTGPSCRizlJ+8cUXzeNe1mCJRCJpLtJihGJFjIDsnlvHjh1TjL/AzsbILbpqUDXXtGnTMGbMmJR7VdlMfvrpJ0SjUYwcOdLz+jRBkNUb4rfdcsfiuYr5T4BxD3L4msrKSnTv3t3xmnLjpXpXwubSuXNnM63XkYkKt8ZUheyN9NRTT9mmi0Qi5vM466yzPOX96quvAjDKXltbi5qaGqUwGTBgAIYOHYqWLVumtBnvvfee53rRp08f3HnnnY5pnISJwDp6Erapww8/3LETJgSNKtTQp59+6noffgzwndxTGYRpgCeilwFEAbQiovUARgIoTV7ncQDnAxhKRHUAtgO4mD2MyeyixspLbCYSCaUwUakGnC45c+bMND2kaNjcvKS8DC+taf73v/+l/J44cSKOOOII3HjjjbYjAbkBFmvKy8ZZZjZdQ50En7zqo8hXuK/KyL1JFU2bNrUVCiUlJcpepejZyeu7yHgxwMfjcZxyyikAjAl01vT19fXm9rZt2zy/xzA477zz8OGHH6b8tlNziTLK77uoqChlrk4kEsE333zjeE23zlQy+B9at26NtWvXKs9RjerscBstqNbikd+lm3Bs3rw5vvvuO7z66qt48803Xd+ZGOEREUpLS21nx9fU1ODxxx9H06ZN0bZtW2zYsAHAzqjBMnajfDf7Sjwex2effaY8JgsTuxH/oYce6nivToImkUi4Cno/BvjVALwYBthnvs6ZMV/icnwCDNdhX6xYsUI5IrC+BNUHIOtondIJ+vTpY86IteLmJWWd3a7Cem2r00AikcBTTz3lGCbCq9B65plnMHjwYF/CT7XPTb0iG3StXHnllWnXtxpXVXgZmVjnYcjlsW7vscceKSrBbM+PkCdLlpaW4ogjjjCfkdyAyz3IVq1amdvl5eU4/vjjU/J0a8DcRiZi37p168xt6+RSMcJll/hvu+++u2NjN2nSJAwdOhTMbHZ4rOVyG/nI8dVqampcO0fChnrggQfixBNPdBz5AEadkctWWlqa4hUJpMb489r5UM0xcSpDJjiVxcsy137UXM/Z/E2B4eVFAGbAmB1f8DCzUtKK8OJOqNJs3rzZ1jjXq1cvDBw4UJmXW4U65xz3UGdiYSk7SktLzQCKdngx7ANGg+zUQ/HqiODWS/31119t0wwePDhtn9X4r8KLMJE/GLm3bydMZJVgtm0m1pnPsdjOtcDlez7ppJOwfPlyAKk2wClTpqTEGBP5OOE2MhG6dXk0bB0Zjxw5EoC7W7XTsXg8jqFDh5orh8r3KwvEJ5980rG8e+65p+NxO1asWOEqSACjcyh7hj7wwANpwtXJAcAOqy0mm6i0Nsccc0x44VSY+Qq7Y0RUBOBuANcBuNxrnvmEiJSS1ktl27w53ba/evXqtNDzAmZO6SHKuL0gL72MuXPnOh6/9957sWTJEsc0Xnt3Yc1QdusRb9iwwVYFo3pmcu/ProFUuXk7IZdRFtjysxK9UDHfI5vIz12o+kR9k2eR19TU4PPPPweQ+i6PPfZYX5NtreevW7cuTW++YMGCtHOsdV22mTh1Nn766Sdb+6Fq/odsxxI88sgjTreDww8/3PwWysvLlR2TILz44ot46aWXzN+/+c1vlBGq/X5H2R71CubMmaP0QPOyZEMoBnhmTjDzPTBUYX8JI89s07FjR2WltfOYkrGqCgBg69atqK6utm0knT5auxENkGrDscNtNv1vfvMb13hSXt0yw5qh7HY9ZsZXX33lOT+596cShvF4HL/8sjPgqZh1bUXuLcrvUm505KWQ7SbkZQP5uU+cODEleKX8sZeWluKwww4DkPosEomE6xwlK/LIZM2aNWnzOFT1ytohE3Y0ZnYdudr11oWnoGDChJ2abauHnRPCtbZPnz6YPn166KPJX375JW0Ua73GYYcdlhKW3wthl1Nl9x03bhzGjBmjTB+qa7BH5gAIV9RnCVmvKTNv3jzXc62qAgDmIlkq7GbNA6lrI6i46KKL8I9//MOxPF4CD7pFh5X95IVAzWTiUlhqLiJKi83lhDV8hFUlYG2k7Bww7HqAcn4333wzAONehSosF8JERswjEfVGDq9/3333mWpNuQF+6KGH0uIw+bGZAEh7rocccggAQ4Bs3boVQLo6cfLkybj00ksDG+DluiXXeT8TA0UePXv29NRAuzkoWLn66qtT8lXd8yGHHIINGzbg/PPP9+wt6uZJ5fdbtXoGAobHpl15tm7dKgSf7UqLYbsG7+10sUIik4ZS4PRRqIa0n3zySUpvVmbUqFGora21zdNNkADhTKQTAeP69u1rCpO77747LZ3cKHs1tttdz4l27dopoxIAUE5clD/gBx54IO0cayNlt7aJF8cCq4cUkHthIt6XCEooN/pdu3ZVvoexY8fi/PPPT8vHD9ZVJkXdE4IESFcDH3nkkQBSPeFUtGjRIiXsifyeVYEMBQ8++KDn8ov7ZWZPk2D9RmGora3Fpk2bzN/Lli1Tjj7EHBsvz//55593XYrYbwQOVYfYyell/vz5YmrDwXbXCW1kQkQnA7gIwOdh5ZlNgggTp56QrEoRnH766bb6epXbrF/chIlXT61YLJaSl6z7FaxduzZNry0asng8bqs+suL2EZWXl9s+Z7c1IeQAkXZumGJipyi3n3hicsOdL2GyePFi3HDDDeb7kuudPIHNquay2uDcZmirFneS5yqp3tHXX3+d8tvLpFlgpzpMXiJZhL23dgbk3yJwpBfE9f/73/96Cv3jddkDQWlpaYrN6NZbb1Xmade5VCFHMbbDqf6p1MVCwFuRl7a2uYatQdXzyISIPrL5m0lEq2BEES4HkFH491yTLWGiytfJ8BukHIIwRibz5s3zZOR7/PHHzQl6ctlFUEw3YTJ79mzcddddri7PzIxVq1bZHnPyhpGfx5133okTTjgBf/zjH5Vpn3vuOfTu3RsjRoxANBrFdddd51guAObyA8xsCpNcx+yaOnUqamtrzXdgVduphAmAtOiwQ4cOdbyOtWERXmICa8gbIH22tJyHdWa1iunTp5v2R7v5O26LytkhyrJu3TrX0CqA/0mUDz74oG0YJsHmzZt9xSNzi5/mhjxfSmBXXz10pmwftp+RSdQh8x9gCJO/MvNHPvLMG7KNIBaLmSv8ecFvw+EU6iAM3MojN3p2XHLJJfjoI2+vToTlFhCRbVBMK8ID6aGHHnJMt3nzZjzxxBO2x528YVQeS9boAT///DPGjRtn6o5Fo+J0TcFpp51mlj9fwuTf//53ym87Rw1rY3XHHXfgz382+nuqkOYC8V1Y36mIoyXw0sOWhcmQIUNc0/Xt2zel/NZ5GkEQ32GHDh3MfWGuoXLooYfaLiAmX0+en+TmsdmiRYtAZVIFfLWb/OjhWX9td8CPa3DooVfyyY4dOzB06FBMnjw5pYfnBaeGwxoKHDDmpdjFtglD0ISh5hLLCHshkUigZcuW5j2tWbPGc2Mqyuq26qHbCMeqlpD13l5GaitWrHA0ODohP898CRPr9dasWWNuO9Up2Q5l14AK21ldXV2ajcTqOeVl2d+FCxea207vXdyTPKG1vr4eN954Y0oeQRDll8PMhLmGihd156JFi8z389vf/tZ2jRihTpbdvjOhpKQkTWCplrsA3OcVAQi+BnxjY+PGjXj88cczOtep4VB5Ce25555ZFSZuDVkikUhrFKz4mcFNRFi4cKHpbJDJc7SLCOAVqyCR151RqV5UyF52PXr0wIIFCzw1BoUgTKxYG3k7NZfcqPTq1QudO3dOEUQAUmxnbipdLx53XjwkgdTZ6fI1RKiSTBk3bhxatmyJqqoqbNy4Me0aYbrdivk9TsjXlh0XrIhVHtetCxbu8NZbb8XYsWNT6oJwHbfiQe1m62C1ywqTIPhtOOzckIFwDLdWlYcVL8LKLayENb9nnnkG5513HoDMGtIwDdZWdYxKmNgZFsViQ2eccQbatm2Ld955x/V64nnK6sN8h6aX+eKLL8yZ+daepghHDgAff/yxcjKaPNfKSegzsydvJ7cFuAR1dXWIx+Npi78FRQRPJCKlw4Qc5uXll18OdK077rjDNVSQjHXxMZnq6mo899xzGQW/lBFOBLKXnd11PQgTWz2YrTAhooznizCz0xK7DR6vPS1Btnutjz76qONxL9Fa3QLkWamrq0tptEX+XkdZXp5Jhw4d8N///ld5TPbAso6orMv9AuoFzQDgd7/7nektY414bEchjkxkFi9ejJNPPhlAeuMgC8v+/fsr1U4PP/ywuX3LLbfYRqueM2eO0rhrxev3UlNTg969e6cIwEMPPTTN6J8psuur/A7llQZfe+21QNeora1VRsiww879HTC+2++++y5wx0vYAWfMmGGqzuwcITyoudShPOA8MnkW3gI7ylDynEYtTNzCl1jJdq/Vy2xyNzJRtckhWvbff3907tzZsxHfC6o5O4LevXunrBwn06VLl7T048ePx7XXXpu2v2PHjua212egEiaFNDLp3r27ZzuZKp0sGK3LAMjMmDHDkxD1q6KS6/MJJ5wQmjCREeouACm2IS82ICdKSkrSImc74WRDKioqQtu2bR3ddf0gRzAW8dKseFm3xA4nYXJlxrk2crzEqZEJakBzw80jxC1aK+Bf7WRtRJo0aZKyzG8YWMOqWxEB/6zrpag+ULv5ArJtwa8wyadrsBNdu3Y1vROd3ntpaalSoMjnON2X7HXlRBCboJ2XZdC1Y+QRr1CvxWKxjBwyZEaPHu1pXoh8bTtKS0sxePBgLFmyxNPyvF5o2bJl2nwymRkzZrhl4d81mJndp17voviZJAWkR1ENm+HDh+Oee+6xPe7XqJwJe++9t7mGQ1i4CRMA5oqOMqoPws25wI8w+eKLLwAYHmdiglohCRNZGIiyqpgyZQquuuqqtPcm38vKlSttz/djG/BDmzZtzG/Gzh4oO1xkgrz87Zw5c3Daaacpvdf8cvDBB9sGdVXhNDIZN24cIpGIp0jmXvnXv/7lGJTUQ1sR3DVYsxO/DUc255gAznpXcX23kUnQMu6xxx6OjgaZ4GVCl4pM9d5en8Ftt90GwJjpvW2bsWp0Iam5mNlU/TlNEvRiGLfOZpfJ1sqSBxxwgGsHLNM1OwTCQQEwIj2I9xf0PSYSCddJizLffvttRscyZdq0aY5RN5wcLoqLi1FfX2/rGuw4piOiyUR0pteCNjQybawKDbcPy0sjuWDBgkCNQ11dXag9KMCTMVCJ6n5VC5pZ8arqk9MJNczf//73rC/b65X33nvPdu6CTCKRcJ1ouv/++9semzlzpu+yecHLe/AyanVCriOy/SQobmH2rTiN/F5//XXE4/GctlPCcUNFUshkHOjxCgDd5R1ENDK5RG6Dp3379vkuQii4Bbrz4s01ZcqUtPDifqirqwu9d+4nKoGMSgj16dNHmVZ8+DNmzMD333+f0fUA4P3338/JOvBemDp1qqd0n376qdLDTcbLCqJh40WYDBs2LLRryGrSoA13IpHwZcR2Srtx40ZUVFQ4RiQPG3nEZoM6QioyjxrcOLr0jYQVK1Y4HvfSU5o/f75rw+JEVVVVwQgTVUiInj17KtOuX78eAPCf//wnsCea39XzsoXX8COzZ892feeTJk2yPdazZ0/Pc0j84EWYqDz2/GD3TbitSOrG8uXL8fHHH3tO7xarrKamxqyjucCDW7N67Qbs4jYTuzkMjQ1mdjTEAsH1s0uWLCkYI7QcAlxgt+SpHFsqSPmLiopCjfEUBK/CxEuP18lA/MADD2QlWrKXPJ3UQ16wEyZr164NlK/XUaFbOQTFxcVpEQqyiZOtqlWrVti8eXNmNhNN4yCRSLgu2xsUZvYUSiJf2AkTuWcYRGVz0UUXhRrjKRd4WdPciXvuuSe0mFkyXkbS7733XqBrZGvJAK8TXwVuarWLL744SHF846T6d5s/o4XJLgAz265fIBPUqFnI2AkTuVEJMmHt7LPPblCCBAinQc2Gp6KXcvldtMpKtjwsPdgcUlAtUiVzzjnnBCmOb5zC3bs9My9qri5EJM9O6gIARNQHNrYTZs6Om4cmo9mwXj+cU045xfcwvaFg57UkP095wSy/ZNv92w8N3UtRpaYMm2y9L7/5HnHEEbbh4AH7RayyhVP53bxGvQiTy5N/MgQgZlcej/lqMqB3796oq6vzFdLlH//4h6cZtEEX4Slk7EYm3bt3x4IFCwAE66mrVrPLFw1dmHid/Bok4nYQz70w83VTE+baDplNYTIT/uNzabJIy5Yt0alTJ1/CxOt8gIbeCDlh1/uTQ+MEESZuDg65xO9Ss0HI9sJvTuy9994Z32u2olL4Hdm7CZNcT4Z1epdu7YOjMGHmaEYl0mQNImrUjX62sFsmWMxgB4IJk0MPPTTjc8MmW5MJVRx22GF5c7ywW13SC36CMfohbMO+2yJyYeNUfrdOAxWSrjeXJIh2zRvXZI1ceLNkxwdJo/FGMTCfmZVeAwVv2yCiyQBOB/A9Mx+uOE4AHgEwCMCvAK5g5gWu+YZdUI0mB+h6qylUGoJr8LMATnU4PhDAQcm/IQAm5qBMGo1GsyuScWyuvJN0M3YKdnMWgOfYYC6AFkTUziG9RqPRaDLjELsDBa/m8kAHAOuk3+uT+xxDgWp1gaYhouutJs/YVsHGIExUN6c0rhPREBiqMG3I1IROtht67TGSW8Qa5BpvNAZhsh6AHOqzIwDlrCdmngRgEgCw9ubSaDSa0Ch4m4kHpgAYTAbHAdjKzOGtdlOAaFWHRqPJE7ah1n2PTIioG4BLARwKYDdmPjm5vwuAngCmMXNoq7kQ0csAogBaEdF6ACMBlAIAMz8O4D0YbsHfwHANvtJLvrkelsgCIOi1RV56aGVPrlUUuRLw+p3nDi91qLGrwnzdHzN7/gPwZwB1MEwOCQD10rH9AdQDuNFPnvn6Sz6j0P6OPvpox+OJRIKZmRcsWOApv5KSEuX+U089lZmZY7FYqOUP669r166e07Zr1y7w9c4++2wuLi5O2//444/n9L6ffPJJzjY7duzIqGyDBg3Ke70I+69du3a230gYf2+88Qaz0VAo/1q1asUA+Kabbsr7s3D6u/322wOdf8UVV6Tts2tTPau5iOhiACMATIOxlO84+TgzfwtgHoBGu2a8E7/8YrtmDICdcW28xtqxC2sg9i9fvtxH6dIpLy8PdL4dHlZqMwka0qK0tBRt27YtiEW5OAeRJDK9xo8//hhuQQqAjRs3ZjVuVVFRkeMSzPvtt1/Wri1o0aJF4DzOOOOMQOf7WXnSj83kJhiqpLOYeQkAVQjJ5TAmD+5y/Pzzz57SeYm1U1xcnBJ/S94WwmTevHk+S5jKAw88EOh8O/xETXUTwG7cd999GDx4cEGsw1LIwsTPmuSNjeLiYk+LnvXt2zfl99y5cx2XYP71118BZL60tBdKSuytEH4WcgvyfSiWDLathH6EyREAPmBmpzjEGwA03jjmDuyxxx6e0nkRJrfccgt69+5t/pYbkbACyRVCA+xnFKPi8MMPRyQSwYABA0IqUeYUsjBp1apVyCVpOIwePdrTomVdu3ZNGa0ff/zxjssfC83AG2+8EbyQNjRt2tT2mB/NQpBv/fnnn7fusl0D3o8wIbhPz2gDQL0KUSPHTZiIIbOXofn111+PXr16ATBGJXJlECodr+t82/H//t//C3R+riGitJ6a6J21bds2H0VKoZCFya4cZXr48OGevpXBgwfj4YcfNn8fe+yxqKqqch0ByO8kbNXx9u3bbY99+OGHnvOxEyZeRjcKFXJzVTrAnzBZAaC33UEiKgZwAoBlPvJsNLh9sM899xwAbyMTIsL+++8PwKjUJ554onlMjEzOPvtsNG3aNON1y3O9TkJQLrnkEvz73/9O2Sc+kiBrt4fFypUrs36NTIXJ7NmzQy5Jw8Lt2ywrK0MkEkG3bt3MfUVFRYhGo7YCoqysDMXFxSkNddgdCqeR+/HHH+85H7vv4+OPP3Y918+oxs9X+CqAo4noFpvjwwEcCOAlH3k2Gtwq7DPPPIN4PI6lS5d6yktUzKOPPjrlmBAmkUgElZWVuPzyyzMqb5C1IMKiQ4cOntOed955OO6441L2vfPOOwDUz/7bb78NVjifPPLII44GW0E8Hse4ceM8pbWSqYozF6OmQsats6HqlBQVFZnfmIqnn34ao0ePxgUXXGDuc3IEOfjgg0MdIfrJy04gOKn/hBC94oorPF/HjzAZD2AxgPuJ6BMY0XpBRH9N/r4HwFwkZ5jvari93Lq6OsRiMfz0k63KMQXRAHz//feYMWOGuX/r1q3mdiQSwaWXXppBaYEXX3zR9zluH6VVJeeGH08RVd4PPvggJk2apCzXgw8+6DnvMKivr3c02AKGIDnxxBNx5513oqKiwrdAyVQoFELHIZ+4fZviuJzOra4fffTRGD58OM4//3xzn5N3YuvWrdGkSRMvxUV5eXmogieTkbuoa506dUo7ZHsdH5lvB9APwPMAjoYxQZEA/D8AxwB4AcCpzNyw9Cch4WUoHY1GMXDgQNe0CxYsMHuhGzZsSOnxWN083fKya9z79Oljbjt5jch4ceW95pprXK8t8FPJv/rqK2U5X3/9deUSrGGveOdGcXExotGoY5pYLGaqOWtqalyFj5VMhckf//jHjM7zSqHbZLzWM+vIBDDemep8sa9Hj53rRNmNYgBgr732cjwuM336dBx44IGe0nohE2GiErBJvrK9jp8LMPNWZr4ChqF9IIDfATgDQDtmvpyZvXW7GyFuH1RlZSUikQgikYhK2qcg91g7dOiQ0ih79RoT/Pa3v1Xub9asmbntVZi8+uqrjseZOeXe/vznPzum91PJR44cqVzH/bzzzsPGjenRc7LRG3d6Tueee66r15AsbETnwg+ZCpMwGyYVhaxGi8fjGQkT8T0Lu4n1+xZp5XPc3r8XrzKRbrfdbJcN8Y1dp85pZCw8yRTtmq0/f0aWS2bewswfMPNLzPwuM2/KJJ/GhJswkSuS23A3EomYH2i7du1SGh2ru6DddXfffXcAxvBaxQsvvGBue1VN9ezZ0/G4MFoK+vXr55jej0qsrq4Os2bNStk3fPhwDBkyBB07dkxL//e//91z3l6RG8327dunHHvjjTdc1VZyHRCdi0yvn4vzColMRz8VFRWoqqrylFYWDJ9++imAnbZJeQQil8eroFqzZo2ndEGxfgvMbFtGp5GxsJn4ee5+ZsC/QkQDiSj/rjNZoqioCD179sR1113n+1yHYWFGeYkGgIhSBIJVffPFF18o8xA9G3m+isyNN95obmfSc1PRpk2blAbSTkUg8DrREzBGBVbhJPTV1oYdAIYNG+Y5b6/Iws86Wc2LzUTGryABdm1hMmjQoIzOq6mpwaZN3vq6S5YsMbfPPPNMs3MQiURw1FFHKc/x+u0sXboUkyZ5Nyd7cdRRobJD2pVR7viVlJSkpBOq9awIEwAXAHgHwH+J6AEiOsLHuQVPUVERysvLMX78eEyc6H/lX/HQ/Xgo2RGPx1OEifxCrQ3D4sWLHctj7VEJZDuM1xGCWzoxK1ggVAR2ldmPAf7ee+9NE4wiX1WFr66u9py3V2RnB6tQ92IzsTJz5kyMHTvWsyG+IQiTbE2G9TsxtbS0FMXFxSgrK/M8D2nBggXmttWmZfcM/TS2Tz/9tOe0shraD6rnbzcNQO7QXHLJJSmdNTGay5YwiQB4AkAZgFsALCKieUR0IxE1+Cm2Z555pqvqQdbDW/XQ4qFnqquXK0Hv3r1ThImMtRHr3r2753yBnT0Q2X/eawPg1guzTg4TKgI7j7N99tnHcZavTNeuXW3LowppkY1G7ZJLLjHfr/VZDB061NdoQ3h2jRgxwrNnV75GGEJl6oWTTz45o2u41QO/RuSrr74ao0ePRmVlJfbZZx/HtOIbk9W4Xm1aTuWyzlFRjaDtcIvLZVdfrPX+888/9zQy69ChQ0qkBFHXFKPt4GvAM/MnzPx/ANoBuBBG6PduAB6BMVp5g4jOJqIGueDWoEGDlI2B3JiPHTvW3LaOQES6PffcM+08L5x00kmmLeXYY4+1HZn873//S6lIhx9+uDI/u+u/9NJLGDNmTIpnSRBhMmTIEBx77LEA1DPRI5FIii++zLPPPovx48d7urbqfkS5//vf9CUWJkyY4ClfPxx33HH4/e9/DyBdfdi5c2dfeYmPlJk9e3bla2Ti1aUVgGvDbYdbPZDrqJdvq3Pnzhg+fDgikYjnb1GossrKyjzbtJyEiVV1dPvtt3sqB+A+MhGToK1Y7/X999/3fP8qRw3FTPuD7c73bf9g5hpm/hcznwFjrfVbYAR4PBvA67BZ5bDQ+eabb5T77733XlNwHHbYYeb+bdu2paQTv4UwcTJWq8IktGrVyvxo7YI8Asa8E7kna1dRxH5rQyL84+UPxaqeskP14Zx88sl46KGHPJXFSl1dnWfjqCoPUR6Vd9yQIUM85SvjxZW5S5cuANLjXfltsDPx7HK7hl2jH1SYePX2A4yecCbnV1VVOT5/+VjXrl1TXNvl/QI/c0as55SXl6cJErtn6JT3L7/8Yh4/7rjjfI1c3dzwn3nmGU/lmTp1qqfnz8wpxnshzBT3bSuZAhnTmXkTMz8M4CgAt8JY6yRY0Kg88fDDDyuHjsOHDzeH4PKLkicPAunzP+wMdvF4HGvXrlUeU9lGVI2ol55sTY0Rj3P+/Pm21xB4NYSrIhX/5z//cW2s7IRJSUmJZzuDkzARwv7UU081j2Uyw/yyyy5zPF5UVGTraOG3wZZn83vtBbtdQxGUD4B9R8krfoTJwoUL0/YNHDjQ9bxoNOrYgMrf3rZt25TP65BDDlGeK96V3X3U1NT4ciFW5W2HUEv7CScfj8ddI2rb2UGsAjmRSOD3v/89rrvuOtf6LdsZRbkVc29CiRqcBhEdQkRjAawB8ACMFRCD1dw84cUbR644VgPv3nvv7ek6dteQGwoiMl+mqrLKPVm7yix6/BdffHHK/iBxrKyuuYBhVBS9Ub8jk7Fjx3ruranKbQ2DIa8xUVFR4SlfGZWLsfVjCkuYyOm9PgO3iZjHHHOMcv9XX9nOM/NUH4LGPtuxIz32a3l5eUrDJ+xrhx56qDKPVatWmdvr1q1LCcoocBs9XHvttcrjtbW1qKioMJ1ZVPn4HZmI2F1CgKnOt/suYrGYrTARz81O8FrLU1ZWhsGDB2PixIn461//qjwHMELNy5E2xDuTjfLJvNfZ5eG7lhDRXkQ0lIjmAvgCwB0wIkk+DaAPM6u7BwWOkzeOqAjyi7Lq6b32PLz0xK1qLmulGz9+vNkAufWMrIElg7guq6KvMrOyN+rlmrLa0A2nkYn4v3LlSnNbjMy8lgUATjjhBMcyZEuYhHWO7I0k4zQB7ve//72pmrXD7lm6IZ7RAQcckHZs/PjxGD16dMq+SCSCgw9Wq+StnTBVz1xWmcrzOkSdECpKFTU1NeakWC/vRjxruU2wjoavueYajBkzxjYPuwnILVu2tH1nzIxrrrnGdja9rCEhopS2wjpqkcv7z3/+E2+99VZafrLTQPK52A5T/cwzOZ2IXoNhE5kAoAeA/8CYBd+WmYcwc4MNT3rzzTe79hCdGiIhXITtxK5CRiIR09Alr4LGzLYuwNbrDhs2zNVm0qRJE2UPRpW+R48enrzQrG7IRITy8nLTo2z9+vVK9ZJdGf14XDkJE3GsS5curj03pzDh1kCSquuJa1kXnFq9erXjuVas9UMVAHLWrFkp+9waublz55rb8n1aDavyCKZ58+au8eL8LAAlN67ieare3eGHH47hw4en7berK/K9Wa8jkKPgPvXUU+ZzEwu2yaMbgVxfxLIPKqwhewYPHpymGpMFXl1dHTp16mTacVTvbq+99lJea9iwYbbvur6+Hp06dbJtq7Zs2ZJiL5U7etbvTS5vfX298pqtWrVC06ZN5W8qlPVMpgA4D4ZKawSAzsw8IDkLvsGvYSJCvqtQjUyslV6oEuzmfciInuBtt92mPG6dtGjFi81E9PysPRhVfuecc07KENeOZ599NkXoXHvttaisrDR7UatXr1aql+waCD/qEyJKE1SLFi1KyWffffdFZWWl8r4F99xzj2d1nNWFWhYmVqcFuyjFdlGCZZVVPB5Hv379cPfdd5vOFc888wz69u2b4jrsJkxkpw+54TjllFPM7bKyshTd+SOPPOKYp1/kMoo5TqtXr854RCw3jKWlpSgqKkJJSYnS3Vy+tlBbx+NxvPvuuwCAJ554Iu2cK6+80qwvdnOy4vE4pk6dmrKvtrYWsVgs5b7kkXsikUDLli1tHWHke7NSU1OjVHOJBt1JuyHm1whEtHIg3WYUjUbN77m4uNicEyaf/9hjj1nbklDCqUwC0JuZf8PM45g53R+zAeNlaCsarWbNmuHaa69V9nK95CNelvxy5fOYWekaLPSwXmwmw4YNQzQaTevB2LnYetHb19fX45xzzjF/T5w4EZFIBBs2bHBUL9mV0W72vgoiShOgn3zyCYCd7yWRSCASiaR5q8n8+OOPtkJMhNAQWF2o5XfRuXPnlPe3ffv2NKERj8fRt29fZZRg+X1/8MEHqK6uRn19vdlReP311817EvtUDhDyvchzjuzWzSGilJneiUTCXJfDKaChdduK3CjJeQPGyMjqaeZVzSdG2OXl5ZgwYQLGjBmDmTNnpqzxY4WIzG8kFouZ5VDZnI466iizvoj7s5YtFoulaQpE/vIzkRfTKioqcl1cq3Xr1mavXyC+b5VrsGjQI5GIrYNJTU1Nip1DRCsH1A40wuvxoosuwvTp0zFmzJiUYK3C49LpmxL4mWdyHTPPdU/Z+LCOTEpKSjBx4kRMnz7drAii0nvxfBFp5Mptp+aS991+++0pFcp6XMZu9GKnLnLzfhKVXBYmgmhyprudesmujLfddptnrysR90v+8MQCQSJ/L5GCjz/+eFMoWznzzDNTfls/IFmYtG3b1pxzAhiz2a0TEGOxmKnbt74P+b7FPB1gp3PFQQcdZF5T7LM+K6FmFMiee7Iwsapf5PchGunRo0ejf//+ac/k4osvNr3l5Dkb1ncqjNtySHbxPvbff/+0eSR24UKs6sO3337brPNDhgxJmzsihGBxcXFKvCxhK5DrjBhpymWXVcZ2wsRav8WIXBZA1nTl5eWIRqOOI7Jt27Zh/PjxGDlypLlP3KtqEqdcH62jIrG91157pQhaueNpXSQtFouZ86PatWtndsQGDx5spiktLfXscdlo42yFiUrNBRgflxAMotKLmFBOPS9xjmxEtKZXjUwOOOAATz0EwH7ugp0wscbRsq4iN2LECFRWVio9hoQnjp16ye6DEqoCLxARIpEITjrpJHOfUOuIcs+aNctsGOyE1HHHHWeW1YrVQ8/J/sPM6Ny5c8ozk0cRQPpcElkNIseaEgug7b777mYjJdSuvXr1MvdZbTrXXnttSiN95ZVXKu85Go2m6L1/97vfmcfkRlrlYHHOOeeYaswbbrghZaQgI0LjyB5x8ndTVVWVUg+EilImHo+nPfOePXsq67zI65RTTjF703L+whgfiURw0UUXATAM4tOnT0ffvn3NdPL7squncv2OxWLmiFzcmyqdeGdOaq6vvvrK1CAIVNoE+fnI6eR3IOpTixYtTI+4Aw44IKXjKX87AGzVcPL1n332We8el3YHiOgjIqokoo7Sby9/3oL2Fxh2jb/8Ap16Gb169cLw4cOVnitWRCNt5ytup+ZSYXfMbu6CnTCRe1VNmzbFY489htGjR5tlve222xCJRBxXbbMTdGHOM5FHH+JDFs4P06dPN0cG1p6bfI4oq1s5nUZ2zJz2QQOpQtwaJVj2NpLVgSJOWvPmzdOenxDe48aNS6ujEydOTMlTHo04NXLykq/y9exGdiKvI4880sznvvvuS0kjjNsqNZfwlJRVXaowQLJKSuAWE6t169Zmb1oePcj1Sghm0fuWZ6LLaZ1UUnb121pnrOnk46rOSXV1dUoEb1F/VROJ5VFvJBLBk08+CcBwib/pppvMdLKrvFxea9mHDRtmer3ZPWe3SOEyTjqZKIwJKs2k315oECFKO3TooAzDAaS+9IqKCtNg7lTZ/BgYxSRB61BfpeaKx+PmkNePEc+uN+HUwFZWViIWi6X0ju6++24zjfzfD7L6QcyhYWbce++9nns9wgAvN/Dz58/HgAED0L59e9NpQfQ0ReNVU1OTEtTSqfzCuCtGKE4jO2ZGJBLBBx98YPZ0O3TogNdee015T9Z9ZWVlZiQEUT5VyJCNGzeiX79+qK2tVaoQxaijpqYGxcXFZl5FRUUpDbNYS0e+nhW3hlvYpCKRCMaNG5eSRujjVWuiy/Xr0ksvxerVq9PcwuPxOKLRaMpzAQy7mCrIo7VHbVd/gZ2zycV7tVsKwE7NFQbMbHZw5Pxl9Rywc6SkEibimCivEMiySkwOOW+t61ZhVl1d7bq8tZ9Yg7ZfFjMXMXMxM38t/fbyl52woSFjjSMlv+BYLGZ+FDU1NeYEnjCESTweN33U77jjDuX1mdnsMcyYMUMVH8f3dZ3Si3v1MroIIkyOPPJIjB492oxrJewCXhAfo8ycOXMAGBFlZVdo0ZhUVlamecw5lT+RSKSoiryMsuQYXe3bt3cUjrKK8O233065LpDaEIue/pdffonq6mokEgmlUV0edcjrdTs1iHbPwE7IqFQhVpWYqNOqkYm4nrwwHDOnddqA9NUKzznnHFd1o8Cu/orRoxAmcgMpp82GMLGOZOVrFxcXY8KECbj88stT1JDRaNQ0wMvfnXXEpbKZWM+RsX4/xcXFpibF7p79RD/QNhMFoockXq6ojKqP0K4xt3s5VmOonF7OS545Kz5KP5Xczmbg5hFkxeoIkIkwEec0b94cw4cPN38vX75cmV6lShMfo9w7F/GZVLpqsf8Pf/iDa96yekQ2PqqwNmJ2H7QKOWyNbHSXRxPxeBxDhw41bSHy7HW7XqJoRK+44grrnAAA6XXBrpwqNZesZpXrnzWmmjhX5aFo51Ivfwt2DiPV1dXK/X4afvH9CtWiXQOZqfuyVyKRCGKxGK677jpcd911mDVrFoYMGaKsv0IlOGTIEDzxxBOmZ6Gqs2J9BnZOEmIUK0bgEyZMSIkaoSKUkYkVIppMRGe6pDmdiCZ7vnoBYe3hqAzKTgvW2L1AK1ZBpbo+YEyKEg1DJiMiu7Dm559/ftp+LxPuwhiZiN6oGFGMHDlSWcZbb71VmY94LwJ5VGDXK7UKD1F+q5uumFXsdeKqU7gbO2RhItyaAaSsC9+vXz888cQTSnfWO++8MyU/67OT661smPca4t7NZiLXUatgF8/ZzmaiwuoYEI1Glb1nN3WjG6Kcc+fORTwet20gg4aNUaGyqUycODHFiC/2y/VX3Nfll1+e4sWmytuq1XC6dmVlpeleLQdDtXuOWREmAK4A0N0lzZEALveRpytEdCoRfUVE3xDRHYrjUSLaSkSLkn9/CuO6spFWrAfwf//3f07l9Jyv28Q6oY8X6UQAQ78Tn1Q9OuFBJTcu48ePd21sMhEm1jyFqkrch108NHkVSIGI/6VSSzhhTSN+y95rbrOKZazP248wkWdoy27IYoS6fv16VFdXp7xnufyjRo1KyU8lJES9tRr7vXjNuY1MrDYYuZctYsDJgkPMGlcFSF22bJmyRy4LDtF7dlI3ehEm69YZ4aQ+++wzVFRUpAhymWzaTDLFi/ONrNWQnXdU34edg0AYaq6w1x4pB6BWvGYAERUDeBRAfwDrAXxGRFOY2TrbbRYzn+4nb+sMUztDlOwVpNJZ271sawiKSZMmmT0B2RgqYGazUZk3bx5OO+00M52Yz+Cnktu5BgsPKrlxEY26U2PqVZhY9eCVlZUpH6nwgKqpqfHlvqyK/xVEmHgphwqnj89NsIgRGZDqzWWN7CzTu3dvM8CmtbG3GmRlZMO81/tT2UysDZWMXI/FaHL9+vXmcWHrE+fH43FzwuQtt9yCo48+Ou1bkLdnzpzp6kji5ZuQe9c1NTW20R6ybYDPJL3f8PmrV682JwN7OdetzmZrZAI4eGoRUTmAvgC8B/JxpyeAb5j5W2auAfBPAGeFkfFPP/2U8rBXrFihTGenlrJD9sKRueGGGxx7/1u2bDF7bRdccIFnl2S/rsGikkaj0RTVhF1jI9KL2EhuFVTu8YvGzuoBpbJvuN3TkUcembbPS2wvO2HipRwqVL10p3LL2C1LMGXKlJTfsnOIdaKZQGWQlcnk/uzUXF4abtH7X7BgQUpEAlFWILVjJs/MtsNL58ZLIz1gwIA0dZpTnmHip5xO5ztRVVVldjhWrVqFP/3pT57PFWR9ZEJE1u76zUR0pSJpMYDWMEYmj3u+ujsdkBryeD0AVTS2CBEthhGE8lZmXqbKjIiGABgCAG3atEF5eTl27NgBZrb1LJJdDnv06JES5wjY2ZubO3duyqSgdu3apbhnuvX+N23aZL58ux5nGK7BohzDhw9H//798f777+OOO+5Qpo/H46Yr4ymnnILKykrbNSMEqh6/eEayG6eXxq1Zs2Zo06YNVq1a5bhsrxN2NhM/5ZDJtMG59957bdcif+edd1J+y0EFVQ381Vdfjf333x9RiwusFb/3Z519Dhj369Zjj8fjeOONNwAYIWlKS0vNUXZ9fb35zETHzOtoKR6PhzIycXIblsmFzcQrXkYmYuLn5s2bce+995rfqnj2Xq7tlsZP+d2eXhEASv6xtG39qwWwFMB9ANTRCzNDdSfW2rMARtDJIwH8HcBbdpkx8yRm7sHMPTp27IjKykozUqhTeGqhZ7RO4InH46a64vTTT08ZTbRp0yZlLRERXsGOsrIyZewtILORiQpr3mINFqew31YB5/bBqXrEftUH4polJSWO78XLvYfdQIhrLly40NcCXCNGjMD111+vPGZVL7k9p/79+3uOhOCVeDxuhmCX+eqrr8w5D6pZ60D6ZEMRPFF0rvyMBq1qUrtnLN7DF1984ek92DloyOcWkppL4FTH5UXPEokESkpKUkI6+an7Ydyz49WYuQsz78fM+8Fo2B8Wvy1/BzJzL2a+k5m9rQHrjfUA9pV+d4RlWWBm3sbMPye33wNQSkSpa6raEIlEbCOFekHV2ErlMtfH6Nmzp/LjkSvy0qVLU6JzenEB9Mvw4cNT8hahuVeuXKlMb1XxRaNRTxXUr5HPiuoDUp3rpREJW3UhbGvz589Pa+yIyDZKMGAffNGql3Z7xmHekxxHTDUKGjlypOmefP311yvvS4xGZffq4cOHm50VubxOc5lEOQROjgNi9LZs2TLP3moCO4GV7XkmfvAyMjn55JNN9Z0cY82PmsutfH6eq59uWz8A//CRPgw+A3AQEe1HRGUALoYRCt+EiNpS8okQUU8Y9+RtYXGJTCpQNBpNmyyneoHC0GhF/lASiYRtdM6wRibymi3xeNw0RNo1fk6jjEzIRJhYK7vXnqvA6rmTaaMjkOf8WBu7n3/+GRUVFSmh5GXsjJnW2eTdunVzLFdQYaJ6hnYdhbq6OvPZ29k57EYcqqVf3RCOA272DWHzVL0HN+w6gdm0mWTjfOtzFy7EYv2aTA3wfr8xgZ+owTOYeY3X9GHAzHUAbgDwAYDlAF5l5mVEdB0RXZdMdj6Az5M2k78BuJh9SIZMXnZ1dbWpz83EiCsQPX9gZ3wsJzLxIJKxLuQjPnYnY6i1Jxlk0mIYwsRrz1WgipTqhNuHdPbZZ9s2dlu3bsWOHTtSQsnL2K0dYlUzqsKPywRV3ameYSQSSZm3IygtLbWNeSWjGnGId/bWW295bpC8flP9+vXzVC4Vqk4gkLqUgZ9ORzZx+76dnnumBniVE40XfLsGJ722joVhHFcuW8fMz/nN146k6uo9y77Hpe0JMFZ+DHod1zRivYvq6mrT7dXJyOmWZyQSwZgxY3D77bc75hPUI8SaD+DfGCoIOmnRbzmt54qeq9dyn3TSSSgvL3eMtyUjPiRrBGCBkzF3zz33xNq1awGoPeRUTgRA+mRYu7AmgqC9XbtnuM8++6SkO/DAA/Hcc8an7Ga8VrF582YAwKuvvoq3337bc4fLi+OAV6O6n3OFxyIzp3zfYZBt1+BMz1XVJZUTjRd8CRMiugrA/QD2sksCw0AemjDJNn4+TFll4uTjDzjPRJUR4aKd1pAPS81l1V1n8jHmS80l8FvuSCSCv/3tb+Z6G254+ZDsGjtZKKjuVY7HJSOvZwFkX5jYPUOrG6gcdTaTRlUEUpUFc5hOA5l44zmdK9xrAffv2ytBvbkyOT+T6AyqCCB+2wbPwoSITgXwFIBlAO4F8CAMz6lPYUQUPgXAa7CMIhoKXhq6iooKjB071vOEu7AMeWICpMpQnqkwATL7GHOt5hJYK7ufcsv5ufU4g/R4ZRKJRNqo5sEHH1SmtS5FkG1hAqifodWms2rVKkf3XDfkAKl+VVH5oF+/fr4nenql0EcmThNSveKnpLfAMGz3ZuaHk/sWMfNfmPlUANcAOBeA2jWoQPHzYQad6JbJ9ePxuDkP4ZFHHslIlyuuEUYjlG81VybILq1edMBuHkde8LNWi3VEkAthosIqTFauXOnbU0oQj8fN2fDyyoeFTFAbqNfoxn4IMjIJIkwywU/LcDSAfzPzT6rzmflpALMB3BW4VDlEzFT3EuwQ8N7QyDFynHB7mTEplpWqtxtkZJIJmeSRb2Eizw/KVQ/5/vvv99wYWQNb5kuYWCMBZ+IpJbBGPrDmXaj47Ui4OWxk05vLDiFMPv/885w6EvgRJrsBkGOE7ADQ3JJmHtQz1AuSeDyON998E4Cxal0YD1718r2MTOwaS6HDt/NaaUjCxMsa7dbzfvzxRwAw4w1lgjyXKIhR1U/9EO6ZXrCG4s+HMInH43jvvVQNdSaeUgK3ettY8OpduG7dOl/1J4iaS6jFFy1a5DqyzHQSrgo/Jf0ORsgUwUYA1tgae8IIrdIgkN1j7SLYBsFPb9ourdvQO1PX4FySaeNXV1dnzsq+8847M67sdosheSFTn/slS5Z4Tis6NIINGzbYpDTIhjCRR8CCTNU9QHCVUUNBeMbZCc1ly4zITmvXrs1IZZjpyESEVnEbWYrJnyKichCB4sebaxlShccsABcTUR9mnkVEhwO4MJmuQZCpe6xfgnpjBfFa8XutbJCpqko2THuJbGyHn8inVmLJVTfluSNeyvDyyy+bq2X6Zdu2bY7Hs/EexUhCrOoIwFyCIVPCqreFjJvDhli+AvDnIRZkZNKvXz9zyWq3dq20tDRN8GT6zvwIk6kAxhNRe2beAMNF+AIAMSLaAmBvGK7BYzIqSR4Iy3NHhVebiZw+E2TjelFRkaOKJN/CRKXmcuoJlZSUoKioKLCwDyJMMu1wLF26FF9++WXG13UiG+9RfAsTJ07E888/H3r+jRknoenFA9SJTN61n3bNj+Bxw48weQKG6+8PAMDMXxBRBYARAA6AYS8Zz8wfZFyaPBB270m8/CVLlphzR8KaJ+J0fps2bXDTTTeZq/GpXDrzLUw2bdqUVi5hqFUJ05KSEkydOjWwsPeydIAdmXY4mNnV9pEp2XqPkUgEVVVVWpiESJD6A2SumvbaroXZofYsTJi5FsD/LPvmAvC1KFVjR4TxXrhwoeMyvwLRMKxcuTKQT781KrFqPsXcuXPNddPDwGt5Fy9eDMDQz1rLFY1GzZ6RNVQ9EI6wDzIyybQMYm7F9u3bA13bLu9sIQuqIPVRs5NM6k+QSYt+CatDnR+LbCOmurraDEAnGkenCiE8eb755pvABjpVTJ14PG5WzAEDBgT22MjEIP3111+b21aDoGyotSJioAUlqDBxwq58J554ou2yzEHJZgMjIgQD/hwONNkhX9qETNDCJGTkAHReVgJUGej8YI23ZXXHtIuQmil+Ay0CwMCBAx09XoRvv0B4cO3YsSOUBi2bwsTu/lu1apW1Xn02Gxh5NB1GfdFkRjbWVMk2tmouIkrAYZleB5iZw15bvsEg6yD/97//4ZFHHnH8+E877TQ89NBDoRjo7PSfYRnYgMzWFverl33zzTdD8zABUoVJ2Kobu/t3W5EyCNkUJnL4e6flnDW5oSGNTJwa/ZnITJjs8ggd5Pjx4z2lDdOjzKr/zEb+meTnppeVRx+TJ09GWVkZ6urqQhGAcoDOsKPB2uVzwAEHhJK/ii+++AIVFRVZy1/QkBqyxkajGpkwczSH5dilCWIAy+U8lWzlB6R6ddXX1+Oaa65Bp06dQhGAsqqmuro61Oi1dio4v7P9/XDbbbehR48eWVGjCWcJYOc6N9oInz8akkDfZdVRuSAXvYuGVNmckL26xNKvYTViLVu2NLcTiUTK76DY2RTk9bnDpra2NmuN/OGHH25uN+YwKJrwyUiYENFuAA4GsDszz3JLv6uSS/e+hk42J5BWVVWZo56ioqJQgw5Go1HlPJlsTVgEjPqUrUZexBQ79NBD8fTTT+tRSZ5oVGouFUTUEcAjAM6AEYOLRR5EdAKASQD+j5lj4RazYZNNYdKYBFW2wm9YRz1hNsSRSATt27c3F4ISNGvWLG1t97C48MILs9bIixA2Xbt21YKkAGhI37efxbHaAfgEQBsAUwDsA0CubZ8k910EIBZeERsuWs1VGGRz1AMAu+22W9q+V155JWvvf4899shKvsDOiMVe3No12aMhajX8jExGwhAWJzNzjIhGQhImzFxLRLMAHB9yGRs8DalCNFayGXRQNSPdunpimBx99NFZy1sLk8KgIQoTP5MWBwGY4qLCWgugfaASNSIaot5T459cf/BHHHFE1vIWQtC6+qMmt9TU1AAA5s2bl+eSeMePMGkDYIVLmloYi2hpJLTNpHGT63eQzevpkUn+icfj+OGHHwAAF1xwQYMJaeNHmGwBsK9LmoNhLKKlQW5GJnr0k3+0MNGEiRwCSbiBNwT8CJPZAM4koraqg0R0EIBTAUwPo2CNiVw0NnqEkj9yvYJlLoTJ0qVLG0yPuLEhvA8b2pLHfr6CBwA0ATCDiAYCaAYYc06Sv/8NIAHgwdBL2UDJ5ahBj1DyR2MamYgIz2Es46rJjIa65LGf9Uw+IaIhAB4H8I50SKwxWgfgKmZuMMv25gptM9GESTbfuQiKGVaQTU1mNMQlj32Nz5n5GQCHA/gbgE8BrASwAMBjALox84uhl7ABo20muwaNSc111llnOS4XoNHY4dv/j5lXALjZ7jgRtWbmTXbH/UJEp8KYdV8M4Clm/ovlOCWPDwLwK4ArmHlBWNcPA20zadw0JjVXtid4ahovoTmTE9GeAP4I4AYAzUPKsxjAowD6A1gP4DMimsLMX0jJBgI4KPnXC8DE5P9dCj1CyR+NSZgADVPFosk/nsbnRNSZiM4lojOIqI3lWBMiGg7gWwB3eM3TIz0BfMPM3zJzDYB/AjjLkuYsAM+xwVwALZKhXzSanOClcQ9TAMirIWo0hYJrw09Ef4NhG3kNwFsAVhPR/yWPRQF8BWAMgKYw1E37h1i+DgDWSb/XJ/f5TQMAIKIhRDSPiObJy+Vmi1yOFrSaK394efb77LNPaNe79tprtZeVpuBwFCZEdDkMtRUALAfwJYBSAH9LHnsfRviUJwAcyMw3M/P3IZZP9ZVaW2gvaYydzJOYuQcz92jdunXgwnklFw29VnPlj1yPTMSiVRpNIeFmM7kCQA2AfswcBwAi6gtgGoCnYYwCzmDmbI271yN11n1HABsySJMXdAO/a+BFUFRXV6ekD1I3SktLtZeVpuBwU3N1A/CmECQAwMwzYai7CMa8kmwqcD8DcBAR7UdEZQAuhhH+XmYKgMFkcByArcy8MYtl8o325mrceHn2zZo1M7fFXI5MefLJJ7WBXFNwuAmTPQGo1h8VAR+zqrhl5joYarYPYKjZXmXmZUR0HRFdl0z2Hgzj/zcAngTwf9kskx/0yGTXwIsw2X333c3tv/zlLw4p3TnqqKMCna/RZAM3NVcRjEjAVmoBgJm3h14iC8z8HgyBIe97XNpmANdnuxxB0DaTxo2XSYtymm7dugW63uLFi3HYYYcFykOjCRsvbry6lcoQ3cDvGnjpLMjCJOhaIZ999lmg8zWabOClVo8iolGqA0RUr9jNzKxX1pHQNhONHNI9E2FSUlJiLlzVs2fP0Mql0YSFl5EJ+fzLbaCiAkaPTHYN/Kq5li9f7vsaV1xxhbn2ezaX7dVoMsXxK2Dmokz+clX4hkIuRg0//vijnsiWJ3766SfXNPLIZNGiRb6vse+++6JFixYA9ChUU5johj+LiJHJ7Nmzs9bQL1hgxLT88ccf9foTeSAej2PJkiWu6WQB0KuX/9BxiUTCFEi5jlKs0XhB18os8v33RjCAGTNmZK2hX7x4sbkt1p/Q5A6vz/uXX34xt2Vh4GeUoYWIppDRtTOL7LPPPuZs52w19AMHDmyQS3w2FqLRKMrLy13Tffnll+b2VVddZW7L6i8nYcHMZlqxtK5GU0hoYZJFTj755Kw39JFIBB999FGDW+KzsRCJRDB+/HjXdLIzRm1trTkikfeXl5fbjlS0MNEUOtqFN4vkaqEhvf5EfqmqqnJNI0aoRITS0lIQEerq6lBWVobt2425v5WVlRg1ahQ+/PDDtPPXrVunhYmmoNHCJMvohr7x42XE2b17dyxcuBCDBg3CXXfdBQBmJ6N3794AjLqy3377Kc9fvXq1FiaagkYLE4na2lqsX78eO3bsyHdR8kaTJk3QsWPHwMEIdyXkzkJ5eXlKhGCBsIfsv//+ZnprJyMej+Prr79WXmPr1q344YcfAAALFy5E9+7dwyi6RhMatKtOrOvRowfPmzcvZd+qVauwxx57oGXLlrukLz8zo6qqCj/99JNtD1mjRtSXOXPmmCMNmaKiIiQSCZSWlmLGjBmmIInH42b6pk2bonfv3qisrHS8Vnl5OaZPn65HvJqcQ0TzmbmH6pg2wEvs2LFjlxUkgNEgtmzZcpcemQXFroEXnbZEIpHi1ReLxcxRS01NDbZt2+Z6Db04lqYQ0cLEwq4qSAS7+v1nC2ZGUVFRmlefcC0WHn+CkpKSFLdhgSoPjaYQ0DYTjSYHFBUV4eSTT8aoUaNSRi+yx1/Lli0xdOhQ89g111yDxx83VlsYMGAAzj33XFRVVWXVM1CjyRQtTAqM9evX4/rrr8cXX3yBRCKB008/HQ888ABeeuklzJs3DxMmTMh3ETUOqKIciFGHVZAIhMffuHHjzH3MjE6dOpm/BwwYgCFDhmSn0BpNCGg1V0Di8TjGjRsXSqgUZsa5556Ls88+GytWrMDXX3+Nn3/+2XQl1RQ+zz33XNo+rxNKrSovWZUlq8A0mkJEj0xsGDZsmGt0161bt2LJkiVIJBIoKipCt27dsOeee9qm7969u+Ns6Y8++ghNmjTBlVdeCcDo0T788MPYb7/9MHr0aKxbtw6nnnoqVq1ahUsvvRQjR47EL7/8ggsvvBDr169HfX097r77blx00UWZ3LImSwwfPtxTOqdJrlqYaAodLUwCsHXrViQSCQCGl87WrVsdhYkby5YtwzHHHJOyr3nz5ujUqRPq6urw6aef4vPPP0ezZs1w7LHH4rTTTsOaNWvQvn17vPvuu2aZNPlj8ODBmDx5MmprazNaz8ZukuvatWvDKJ5GkzW0MLHBS7yleDyOiooK1NTUoKysDC+++GIgw6gIt2G3v3///mjZsiUA4Nxzz8XHH3+MQYMG4dZbb8Uf//hHnH766ejTp0/G19cEJxKJIBaLIRaL4c477wyUl6w6vf/++zFo0CBteNcULNpmEgChlggryOJhhx0G60TKbdu2mXGZrIKGiHDwwQdj/vz5OOKIIzB8+HD8+c9/DlQGTXAikYhn1ZYT8lyS+vp6PbdEU9BoYRIQ0XCE0WOsqKjAr7/+ahpx6+vrccstt+CKK65As2bNMG3aNGzZsgXbt2/HW2+9heOPPx4bNmxAs2bN8Lvf/Q633nqruViWpuFjNcDruSWaQkYLkwKCiPDmm2/itddew0EHHYSDDz4YTZo0wdixYwEAJ5xwAi677DJ0794d5513Hnr06IGlS5eiZ8+e6N69O+69916MGDECAPCnP/0JU6ZMyeftaAIid1D08gKaQkfH5pJYvnw5Dj300DyVqHDQz8E/qvVJ5HhdmQoCVb4aTb7Qsbk0mhwjG8+ztWSzRlNIaGGi0WSBWCxmjirCWLJZCyNNoaOFiUaTBaLRaOAlm/XoRtOQ0PNMNJosEMaSzWJ0w8zm6EYb4TWFSsEKEyLaG8ArALoAWA3gQmb+QZFuNYCfANQDqLMzDmk0uSboks1idCMmxWrXYE0hU7DCBMAdACqZ+S9EdEfy9x9t0vZj5s25K5pGk33CGN1oNLmikG0mZwH4R3L7HwDOzl9RcgcR4bLLLjN/19XVoXXr1jj99NN95dOlSxds3uwsX72k0eSXMCfFajTZpJCFSRtm3ggAyf/72KRjAB8S0XwiclzwgYiGENE8Ipq3adOmkIsbDrvtths+//xzbN++HQAwbdo0dOjQIc+l0mg0GmfyquYiov8AaKs45GcBj+OZeQMR7QNgGhF9ycwzVQmZeRKASYAxadEx12HDAJcQ9L7p3h3wEEBy4MCBePfdd3H++efj5ZdfxiWXXIJZs2YBALZs2YKrrroK3377LZo1a4ZJkyahW7duqKqqwiWXXIJNmzahZ8+eKZPcXnjhBfztb39DTU0NevXqhccee0y5JKxGo9FkSl5HJsx8MjMfrvh7G8D/iKgdACT/f2+Tx4bk/+8BvAmgZ67Kny0uvvhi/POf/8SOHTuwZMkS9OrVyzw2cuRIHHXUUViyZAnGjh2LwYMHAwDuuecenHDCCVi4cCHOPPNMM2T58uXL8corr2D27NlYtGgRiouL8eKLL+blvjQaTeOlkA3wUwBcDuAvyf9vWxMQ0W4Aipj5p+T2KQDCCZvrYQSRLbp164bVq1fj5ZdfxqBBg1KOffzxx3j99dcBACeddBKqqqqwdetWzJw5E2+88QYA4LTTTsNee+0FwIjpNH/+fBx77LEAgO3bt2Offew0hhqNRpMZhSxM/gLgVSL6PYC1AC4AACJqD+ApZh4EoA2AN5MzjUsAvMTM7+epvKFy5pln4tZbb0UsFkNVVZW5XxWjScy0tlsL5fLLL09ZX1yj0WjCpmAN8MxcxcwVzHxQ8v+W5P4NSUECZv6WmY9M/h3GzPfmt9ThcdVVV+FPf/oTjjjiiJT9ffv2NdVUsVgMrVq1QvPmzVP2T506FT/8YEzJqaiowL/+9S98/72hJdyyZQvWrFmTwzvZtdCz1DW7KgUrTHZ1OnbsiD/84Q9p+0eNGoV58+ahW7duuOOOO/CPfxje0yNHjsTMmTNx9NFH48MPP0SnTp0AAF27dsWYMWNwyimnoFu3bujfvz82btyY03tp7OiwJxqNDkGfsk+HXjfQz8Ef48aNw4gRI5BIJFBcXIzRo0eHstKiRlNo6BD0Gk0WiUajKC8vDxTUUaNp6BSyAV6jaRDosCcajRYmaTCz0itqV2FXVXsGJWhQR42moaPVXBJNmjRBVVXVLtugMjOqqqrQpEmTfBdFo9E0MPTIRKJjx45Yv349CjVuVy5o0qQJOnbsmO9iaDSaBoYWJhKlpaXYb7/98l0MjUajaXBoNZdGo9FoAqOFiUaj0WgCo4WJRqPRaAKzy86AJ6JNAAo5SFUrAI1tGcTGeE9A47yvxnhPgL6voHRm5taqA7usMCl0iGieXdiChkpjvCegcd5XY7wnQN9XNtFqLo1Go9EERgsTjUaj0QRGC5PCZVK+C5AFGuM9AY3zvhrjPQH6vrKGtploNBqNJjB6ZKLRaDSawGhhotFoNJrAaGGSRYjoECJaJP1tI6JhyWM3EtFXRLSMiO5P7islon8Q0VIiWk5Eacv1EdEUIvpc+l1ORK8Q0TdE9AkRdWko90REZUQ0iYi+JqIviei8fNxTFu7rkuT+JUT0PhG1ysd9ZXBPZUT0TLLsi4koKuV1THL/N0T0N0qu09BA3pXyvoioGRG9m6x7y4joL9I1Guy7kvLMbVvBzPovB38AigF8B6AzgH4A/gOgPHlsn+T/SwH8M7ndDMBqAF2kPM4F8BKAz6V9/wfg8eT2xQBeaSj3BOAeAGOS20UAWuX7noLeF4zgqd9L93I/gFH5vi+P93Q9gGfEPgDzARQlf38KIAKAAEwFMDDf9xT0vpLvrV9yfxmAWYVwX0HfVXJfztuKnL30Xf0PwCkAZie3XwVwsiLNJQD+nWyQWgL4GsDeyWO7A/gYQFdLBfkAQCS5XQJjFiw1kHtaB2A3xTl5u6eg9wWgFMCmZENAAB4HMCTf9+Xxnh4F8DvpdyWAngDaAfjScu9P5Puegt6XIt0jAK7J930Fvad8tRVazZU7LgbwcnL7YAB9kkPNGUR0bHL/vwD8AmAjgLUA/srMW5LHRgN4EMCvlnw7wGiUwcx1ALbCaNxyQcb3REQtksdHE9ECInqNiNok9+XznoAA98XMtQCGAlgKYAOMD/rp5DmF/q4WAziLiEqIaD8AxwDYN1nu9VJe65P7gIbxruzuyyRZH8+A0SgDDfddAXlqK7QwyQFEVAbgTACvJXeVANgLwHEAbgPwalIH3RNAPYD2APYDcAsR7U9E3QEcyMxvqrJX7Mu6v3fQe0qm7wijB3Y0gDiAv4rsFZfMiQ97CO+qFIYwOSp5bAkAYU8p9Hc1GYagmAdgPIA5AOrgXO6G8K7s7kvkUwKj8f4bM38rdisuWfDvKp9thRYmuWEggAXM/L/k7/UA3mCDTwEkYARquxTA+8xcy8zfA5gNoAcMXfUxRLQaxvD1YCKKSXntC5gfxZ4AxGimkO+pCkbPSVT61wAcLeWVj3sCgt9XdwBg5pVs6BNeBdBbyqtg3xUz1zHzzczcnZnPAtACwIpkenn5zY4wRl0ir4J+Vw73JZgEYAUzj5f2NdR3lbe2QguT3HAJdg5bAeAtACcBABEdDMP4txmGuuQkMtgNRm/kS2aeyMztmbkLgBMAfM3M0WReUwBcntw+H8BHyUYs2wS9J4Zhc4gmz68A8EVyO1/3BAS8LwD/BdCViERk1f4Alie3C/pdJb2bdkvu7w+gjpm/YOaNAH4iouOSveLBAN5O5lXw78ruvpK/x8BoVIdZ8m6o7yp/bUVYxhf9Z2tMawajF76ntK8MwAsAPgewAMBJvNNw9hqAZTAa1tsU+XVBqlGtSfKcb2B43OzfUO4JhpF6JgxVUCWATvm6p5Dv6zoYAmQJDIHZsoG8qy4AvkqW/T8wwo2Lc3ok068EMAE7o2c0hHelvC8YIyxO7l+U/Lu6ob8r6dwuyGFbocOpaDQajSYwWs2l0Wg0msBoYaLRaDSawGhhotFoNJrAaGGi0Wg0msBoYaLRaDSawGhhoik4iOgKImIiuiLfZfECET2bLK/4uyPfZSpEiChGRBm7jxLR55bnHA2vdJqglOS7AJrGTQaNx5VZKUhueATAjzBmHmvC5zEYEXKjAE7Mb1E0VrQw0WSbexT7hsGYdSwaX5lFAFYBmAsjiGJDYjwzr853IRorzPwYABDRKGhhUnBoYaLJKsw8yrovqb7aE86N79bslUqj0YSNtploCg47mwkRrU7+7U5EDxPROiLaTsbKdGcn05QQ0Z1EtIKIdhDRSiK6weFaA4joPSLaTETVyfQP0M4Q+WHcz/5krCj5TbK8W8hYIe9xIkoLAU7GSo3TieiH5D0sJ6IRRFRuk/9viGhy8tlUE9H3RDSLiIYq0laQsfrjlmTeXxPRX4hoT0XaWPI9yM+0Ovnc70tGuFWV52Iimp+81++J6Hkiam+TlojociKaQ0SbkmVaR0QfENFF7k9XUyjokYmmoVEKYBqMhajehhG76BIArxPRKTBWk+sFYzXAagAXAPg7EW1i5lfkjIjoTzDUcFsAvANjhcRuAG4FMIiIIsy8LUhhiagdgM8ANAfwHoDXYcRI2g/AZTDiXFVJ6Z8GcBWS0WJhqAGPg7FGRQUR9WdjLQqR/jQY8ZbKAbwPI0hgCwBHArgdwEQp7bXJ378kz/kehv3hjwDOIKLjmflHxW28BKAPjGe6DcCgZN77wGLjIqKbATyULPdzyf8DYIRIV40274URon8VjAjLW2EsxnUsjHf3iuIcTSGSi6Bs+k//yX8wlrhlSEsSW45fkTx+hc15/0ZyGdPk/j7J/VtgNNwtpGP7A6gBsNCSV7/kOXPk9JbrP+zxfp61ux8ANyaP/UFxbDcATRXXfUPenzw2ypoPjFD4W5P3d6Ii/47SdmcYwnUbgN9Y0j2WzHuSZX8suX8+kqtjSuX+BsZ6Lm2l/V2S19iC1OWmi2AIUTaanJRrVMEQnM0U5W9l87zFs4jmuy7rv51/Ws2laYgMY+Zq8YOZZ8Ho2e4F4I8s9a7ZWOxoNoAjiKhYyuOm5P9r2NIbZ+ZnYTgC/DbEMm+37mDmX5hZ3v8HGIs2XWXZDxgjkypLmS6HMeKZyMwzFPnLqyP+DsYobgIzf2lJeheAnwBcZqNK+yPvXPETzPwLgBdhCIkeUrrfJq/xd5ZsYcycgLGwU0KRNwDUwhBM1vJvtkmvKUC0mkvT0PiRmVcq9m+AoTqarzj2XwDFANomtwFjEaFaABcQ0QWKc8oAtCailsxcpTjulSkAxgJ4lIgGwFiHezaALzjZzQYAImoGQzW1GcAwY9mQNKoBHCr9Pi75f6qHcoiFxz6yHmDmH4hoIYC+AH4DY0lYmXmK/NYl/++luIZKsH1LROtgjJBkXoQxeltGRK8lz40zs3bAaGBoYaJpaNg1MnUAYNMICRtDqbSvJYz6P9LlertDsmn4hZnXEFFPGKqZUwGcmzy0joj+ysx/S/7eC8ayqq09lEnQIvn/v06JkggDu527tdjfwnrAOnJLIp6pPNoT1/gf1HyHdGFyM4w1Uq4CcEfyr46I3gNwCzN/Y5OXpsDQai7NrspWAD8wM7n8rQl6IWZezswXwRBgPWA0mEUAHiGi30vlAQzbjmOZpKx/TP7v4PF+AWN0pqKdJV0miHPb2BxPuzYz1zPzI8x8ZPK882As5XwmgPftPNg0hYcWJppdlbkA9iKiw3J1QTbW7Z7PzPfB8EADgLOTx36GsWrjYUS0t8cs5yb/D/SQdmHyf9R6IOkG3R3ADuxcYjgTFiT/p00oJKL9kVx/3A5m/p6Z32DmC2Go4w4AcHiA8mhyiBYmml2Vh5P/n1TNgSCi3YjoOOt+vxBRTyJS9dTFvl+lfQ/BsNVMVs1zIaK9iOhoadc/YHhnDSWivor0HaWfL8CwEd1IRAdako6GYch/QXZsyIAXpWt0kcpRBOABWNobIipPznshy/5SGK7fQOrz0RQw2mai2SVh5koyAjKOA7AiqaNfBcNG0hlG7/pjGHaOIFwK4HoimgHDnfYHGD3uM2AY1MdLZZpMRMfAmCuzkog+ALAWRsO6HwwD+TMw1pgHM28moksB/AvAdCKaCmPd+eYw5svsmzwPzLyaiIYBeBTAAiJ6FcCm5H1GAHwJY75JxiSvcQeABwEsJKJXYKi+BsCwxSxJlkvQFMb65auJ6BMAa2DMwekPw9FgCjMHGSlpcogWJppdFma+j4hmw3ATPgHAWTAav/8CmARjsl5QXoYxobA3DG+npsn8/wngQWb+3FKm65NC4ToAJ8NohLfAECoPwBhhyOnfJaIeMARBBYBTYAisL2EISjntY0T0DYxJmecBaAbDK+sBAGNtDO2+YOaHiGgjDFfgK2C4HH8AY5Kj9Xn+kix3PxjP5+xk+pUAhgKYHLQ8mtxBkneiRqPJACJ6Fsacj/1YB3rMOmQEehwJoB8zx/JbGo1A20w0mvBYRXo9k6xByfVM4N11WpNDtJpLownOWzBCvQj0eibZQaxnIlidp3JoFGg1l0aj0WgCo9VcGo1GowmMFiYajUajCYwWJhqNRqMJjBYmGo1GowmMFiYajUajCYwWJhqNRqMJzP8Hf67WZCbOaMAAAAAASUVORK5CYII=",
"text/plain": [
"
\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**Finally, the log contains all the details**"
]
},
{
"cell_type": "code",
"execution_count": 36,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Stellar occultation of star GaiaDR3 6760223758801661440 by 10199 Chariklo (1997 CU26).\n",
"\n",
"Geocentric Closest Approach: 0.049 arcsec\n",
"Instant of CA: 2017-06-22 21:18:48.200\n",
"Position Angle: 359.72 deg\n",
"Geocentric shadow velocity: -22.00 km / s\n",
"Sun-Geocenter-Target angle: 166.42 deg\n",
"Moon-Geocenter-Target angle: 149.11 deg\n",
"\n",
"\n",
"5 positive observations\n",
"1 negative observations\n",
"\n",
"###############################################################################\n",
" STAR \n",
"###############################################################################\n",
"GaiaDR3 star Source ID: 6760223758801661440\n",
"ICRS star coordinate at J2016.0:\n",
"RA=18h55m15.65210s +/- 0.0197 mas, DEC=-31d31m21.6676s +/- 0.0180 mas\n",
"pmRA=3.556 +/- 0.025 mas/yr, pmDEC=-2.050 +/- 0.020 mas/yr\n",
"GaiaDR3 Proper motion corrected as suggested by Cantat-Gaudin & Brandt (2021) \n",
"Plx=0.2121 +/- 0.0228 mas, Rad. Vel.=-40.49 +/- 3.73 km/s \n",
"\n",
"Magnitudes: G: 14.224, B: 14.320, V: 13.530, R: 14.180, J: 12.395, H: 11.781,\n",
" K: 11.627\n",
"\n",
"Apparent diameter from Kervella et. al (2004):\n",
" V: 0.0216 mas, B: 0.0216 mas\n",
"Apparent diameter from van Belle (1999):\n",
" sg: B: 0.0238 mas, V: 0.0244 mas\n",
" ms: B: 0.0261 mas, V: 0.0198 mas\n",
" vs: B: 0.0350 mas, V: 0.0315 mas\n",
"\n",
"Geocentric star coordinate at occultation Epoch (2017-06-22 21:18:48.200):\n",
"RA=18h55m15.65251s +/- 0.0323 mas, DEC=-31d31m21.6706s +/- 0.0341 mas\n",
"\n",
"###############################################################################\n",
" 10199 Chariklo (1997 CU26) \n",
"###############################################################################\n",
"Object Orbital Class: Centaur\n",
"Spectral Type:\n",
" SMASS: D [Reference: EAR-A-5-DDR-TAXONOMY-V4.0]\n",
" Relatively featureless spectrum with very steep red slope.\n",
"Discovered 1997-Feb-15 by Spacewatch at Kitt Peak\n",
"\n",
"Physical parameters:\n",
"Diameter:\n",
" 302 +/- 30 km\n",
" Reference: Earth, Moon, and Planets, v. 89, Issue 1, p. 117-134 (2002), \n",
"Rotation:\n",
" 7.004 +/- 0 h\n",
" Reference: LCDB (Rev. 2021-June); Warner et al., 2009, [Result based on less than full coverage, so that the period may be wrong by 30 percent or so.] REFERENCE LIST:[Fornasier, S.; Lazzaro, D.; Alvarez-Candal, A.; Snodgrass, C.; et al. (2014) Astron. Astrophys. 568, L11.], [Leiva, R.; Sicardy, B.; Camargo, J.I.B.; Desmars, J.; et al. (2017) Astron. J. 154, A159.]\n",
"Absolute Magnitude:\n",
" 6.58 +/- 0 mag\n",
" Reference: JPL Horizons, \n",
"Phase Slope:\n",
" 0.15 +/- 0 \n",
" Reference: JPL Horizons, \n",
"Albedo:\n",
" 0.045 +/- 0.01 \n",
" Reference: Earth, Moon, and Planets, v. 89, Issue 1, p. 117-134 (2002), \n",
"\n",
"Ellipsoid: 151.0 x 151.0 x 151.0\n",
"\n",
"----------- Ephemeris -----------\n",
"\n",
"EphemKernel: CHARIKLO/DE438_SMALL (SPKID=2010199)\n",
"Ephem Error: RA*cosDEC: 0.000 arcsec; DEC: 0.000 arcsec\n",
"Offset applied: RA*cosDEC: 0.0000 arcsec; DEC: 0.0000 arcsec\n",
"\n",
"\n",
"###############################################################################\n",
" POSITIVE OBSERVATIONS \n",
"###############################################################################\n",
"\n",
"-------------------------------------------------------------------------------\n",
"Site: Outeniqua\n",
"Geodetic coordinates: Lon: 16d49m17.71s, Lat: -21d17m58.17s, height: 1.416 km\n",
"Target altitude: 56.6 deg\n",
"Target azimuth: 115.3 deg\n",
"\n",
"Light curve name: Outeniqua lc\n",
"Initial time: 2017-06-22 21:20:00.056 UTC\n",
"End time: 2017-06-22 21:23:19.958 UTC\n",
"Duration: 3.332 minutes\n",
"Time offset: -0.150 seconds\n",
"\n",
"Exposure time: 0.1000 seconds\n",
"Cycle time: 0.1002 seconds\n",
"Num. data points: 2000\n",
"\n",
"Bandpass: 0.700 +/- 0.300 microns\n",
"Object Distance: 14.66 AU\n",
"Used shadow velocity: 22.004 km/s\n",
"Fresnel scale: 0.040 seconds or 0.87 km\n",
"Stellar size effect: 0.010 seconds or 0.23 km\n",
"Inst. response: 0.100 seconds or 2.20 km\n",
"Dead time effect: 0.000 seconds or 0.00 km\n",
"Model resolution: 0.004 seconds or 0.09 km\n",
"Modelled baseflux: 1.029\n",
"Modelled bottomflux: 0.109\n",
"Light curve sigma: 0.307\n",
"\n",
"Immersion time: 2017-06-22 21:21:20.175 UTC +/- 0.027 seconds\n",
"Emersion time: 2017-06-22 21:21:30.201 UTC +/- 0.029 seconds\n",
"\n",
"Monte Carlo chi square fit.\n",
"\n",
"Minimum chi-square: 192.774\n",
"Number of fitted points: 200\n",
"Number of fitted parameters: 2\n",
"Minimum chi-square per degree of freedom: 0.974\n",
"\n",
"immersion:\n",
" 1-sigma: 76880.325 +/- 0.027\n",
" 3-sigma: 76880.347 +/- 0.121\n",
"\n",
"emersion:\n",
" 1-sigma: 76890.351 +/- 0.029\n",
" 3-sigma: 76890.348 +/- 0.101\n",
"\n",
"\n",
"-------------------------------------------------------------------------------\n",
"Site: Onduruquea\n",
"Geodetic coordinates: Lon: 15d59m33.75s, Lat: -21d36m26.04s, height: 1.220 km\n",
"Target altitude: 56.1 deg\n",
"Target azimuth: 114.7 deg\n",
"\n",
"Light curve name: Onduruquea lc\n",
"Initial time: 2017-06-22 21:11:52.175 UTC\n",
"End time: 2017-06-22 21:25:13.389 UTC\n",
"Duration: 13.354 minutes\n",
"Time offset: -0.190 seconds\n",
"\n",
"Object LightCurve was not instantiated with time and flux.\n",
"\n",
"Bandpass: 0.700 +/- 0.300 microns\n",
"Object Distance: 14.66 AU\n",
"Used shadow velocity: 22.004 km/s\n",
"Fresnel scale: 0.040 seconds or 0.87 km\n",
"Stellar size effect: 0.010 seconds or 0.23 km\n",
"\n",
"Object LightCurve model was not fitted.\n",
"\n",
"Immersion time: 2017-06-22 21:21:22.023 UTC +/- 0.010 seconds\n",
"Emersion time: 2017-06-22 21:21:33.634 UTC +/- 0.011 seconds\n",
"\n",
"\n",
"-------------------------------------------------------------------------------\n",
"Site: Tivoli\n",
"Geodetic coordinates: Lon: 18d01m01.24s, Lat: -23d27m40.19s, height: 1.344 km\n",
"Target altitude: 58.5 deg\n",
"Target azimuth: 112.4 deg\n",
"\n",
"Light curve name: Tivoli lc\n",
"Initial time: 2017-06-22 21:16:00.094 UTC\n",
"End time: 2017-06-22 21:28:00.018 UTC\n",
"Duration: 11.999 minutes\n",
"Time offset: -0.150 seconds\n",
"\n",
"Object LightCurve was not instantiated with time and flux.\n",
"\n",
"Bandpass: 0.700 +/- 0.300 microns\n",
"Object Distance: 14.66 AU\n",
"Used shadow velocity: 22.004 km/s\n",
"Fresnel scale: 0.040 seconds or 0.87 km\n",
"Stellar size effect: 0.010 seconds or 0.23 km\n",
"\n",
"Object LightCurve model was not fitted.\n",
"\n",
"Immersion time: 2017-06-22 21:21:15.478 UTC +/- 0.011 seconds\n",
"Emersion time: 2017-06-22 21:21:19.838 UTC +/- 0.038 seconds\n",
"\n",
"\n",
"-------------------------------------------------------------------------------\n",
"Site: Windhoek\n",
"Geodetic coordinates: Lon: 17d06m31.9s, Lat: -22d41m55.16s, height: 1.902 km\n",
"Target altitude: 57.4 deg\n",
"Target azimuth: 113.4 deg\n",
"\n",
"Light curve name: Windhoek C14 lc\n",
"Initial time: 2017-06-22 21:12:48.250 UTC\n",
"End time: 2017-06-22 21:32:47.963 UTC\n",
"Duration: 19.995 minutes\n",
"Time offset: -0.375 seconds\n",
"\n",
"Object LightCurve was not instantiated with time and flux.\n",
"\n",
"Bandpass: 0.700 +/- 0.300 microns\n",
"Object Distance: 14.66 AU\n",
"Used shadow velocity: 22.004 km/s\n",
"Fresnel scale: 0.040 seconds or 0.87 km\n",
"Stellar size effect: 0.010 seconds or 0.23 km\n",
"\n",
"Object LightCurve model was not fitted.\n",
"\n",
"Immersion time: 2017-06-22 21:21:17.234 UTC +/- 0.024 seconds\n",
"Emersion time: 2017-06-22 21:21:27.189 UTC +/- 0.026 seconds\n",
"\n",
"\n",
"-------------------------------------------------------------------------------\n",
"Site: Windhoek\n",
"Geodetic coordinates: Lon: 17d06m31.9s, Lat: -22d41m55.16s, height: 1.902 km\n",
"Target altitude: 57.4 deg\n",
"Target azimuth: 113.4 deg\n",
"\n",
"Light curve name: Windhoek D16 lc\n",
"Initial time: 2017-06-22 21:20:01.884 UTC\n",
"End time: 2017-06-22 21:22:21.894 UTC\n",
"Duration: 2.333 minutes\n",
"Time offset: 0.000 seconds\n",
"\n",
"Object LightCurve was not instantiated with time and flux.\n",
"\n",
"Bandpass: 0.700 +/- 0.300 microns\n",
"Object Distance: 14.66 AU\n",
"Used shadow velocity: 22.004 km/s\n",
"Fresnel scale: 0.040 seconds or 0.87 km\n",
"Stellar size effect: 0.010 seconds or 0.23 km\n",
"\n",
"Object LightCurve model was not fitted.\n",
"\n",
"Immersion time: 2017-06-22 21:21:17.288 UTC +/- 0.028 seconds\n",
"Emersion time: 2017-06-22 21:21:27.228 UTC +/- 0.034 seconds\n",
"\n",
"\n",
"###############################################################################\n",
" NEGATIVE OBSERVATIONS \n",
"###############################################################################\n",
"\n",
"-------------------------------------------------------------------------------\n",
"Site: Hakos\n",
"Geodetic coordinates: Lon: 16d21m41.32s, Lat: -23d14m11.04s, height: 1.843 km\n",
"Target altitude: 56.8 deg\n",
"Target azimuth: 112.5 deg\n",
"\n",
"Light curve name: Hakos lc\n",
"Initial time: 2017-06-22 21:10:19.461 UTC\n",
"End time: 2017-06-22 21:30:19.345 UTC\n",
"Duration: 19.998 minutes\n",
"Time offset: -0.200 seconds\n",
"\n",
"Object LightCurve was not instantiated with time and flux.\n",
"\n",
"Bandpass: 0.700 +/- 0.300 microns\n",
"Object Distance: 14.66 AU\n",
"Used shadow velocity: 22.004 km/s\n",
"Fresnel scale: 0.040 seconds or 0.87 km\n",
"Stellar size effect: 0.010 seconds or 0.23 km\n",
"\n",
"Object LightCurve model was not fitted.\n",
"\n",
"Immersion and emersion times were not fitted or instantiated.\n",
"\n",
"\n",
"###############################################################################\n",
" RESULTS \n",
"###############################################################################\n",
"\n",
"Fitted Ellipse:\n",
"equatorial_radius: 138.657 +/- 0.373\n",
"center_f: -13.613 +/- 0.120\n",
"center_g: -2.094 +/- 0.499\n",
"oblateness: 0.086 +/- 0.003\n",
"position_angle: 123.956 +/- 1.496\n",
"polar_radius: 126.769 km \n",
"equivalent_radius: 132.579 km \n",
"geometric albedo (V): 0.060 (6.0%) \n",
"\n",
"Minimum chi-square: 12.130\n",
"Number of fitted points: 10\n",
"Number of fitted parameters: 5\n",
"Minimum chi-square per degree of freedom: 2.426\n",
"Radial dispersion: -0.232 +/- 0.797 km\n",
"Radial error: 0.532 +/- 0.222 km\n",
"\n",
"Ephemeris offset (km): X = -13.6 km +/- 0.1 km; Y = -2.1 km +/- 0.5 km\n",
"Ephemeris offset (mas): da_cos_dec = -1.280 +/- 0.011; d_dec = -0.197 +/- 0.047\n",
"\n",
"Astrometric object position at time 2017-06-22 21:18:48.200 for reference 'geocenter'\n",
"RA = 18 55 15.6523911 +/- 0.034 mas; DEC = -31 31 21.622094 +/- 0.058 mas\n"
]
}
],
"source": [
"print(occ)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"**You can find more information about each Class at their specific Jupyter-Notebook.**\n",
"\n",
"**The END**"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3.9.12 ('sora-develop-stats')",
"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.9.12"
},
"vscode": {
"interpreter": {
"hash": "e791c7fa47f20d05b023e3f24fcf5d35e048e2b324f60c72e6b67828804630a8"
}
}
},
"nbformat": 4,
"nbformat_minor": 4
}