misc/Jupyter-notebook-access/EPN_TAP_request.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "3bff6f22",
   "metadata": {},
   "source": [
    "## Query all EPN-TAP services"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "9f0667ce",
   "metadata": {},
   "source": [
    "We show here how to send the same EPN-TAP query to all existing EPN-TAP services. This will retrieve all results from many independent services / institutes, including those unknown to the user."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "2949a0d9",
   "metadata": {},
   "source": [
    "EPN-TAP services are first identified from the IVOA registry. Then the query is sent to all of them in sequence."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "a554f163",
   "metadata": {},
   "source": [
    "### Grab a list of all EPN-TAP services from the registry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "82c90a59",
   "metadata": {},
   "outputs": [],
   "source": [
    "import pyvo as vo\n",
    "from matplotlib import pyplot as plt\n",
    "from astropy.table import Table\n",
    "from astropy.table import vstack\n",
    "import urllib\n",
    "from astroquery.utils.tap.core import TapPlus  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "d01266b6",
   "metadata": {},
   "outputs": [],
   "source": [
    "query = \"SELECT table_name, access_url, ivoid FROM (SELECT ivoid, table_name, table_utype FROM rr.res_table WHERE table_utype LIKE 'ivo://vopdc.obspm/std/epncore%' OR table_utype LIKE 'ivo://ivoa.net/std/epntap#table-2.%' OR table_utype='ivo.//vopdc.obspm/std/epncore#schema-2.0' ) AS epntables NATURAL JOIN rr.capability NATURAL JOIN rr.interface WHERE standard_id = 'ivo://ivoa.net/std/tap' AND intf_role='std'\"\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "779d26d2",
   "metadata": {},
   "outputs": [],
   "source": [
    "registry_url = \"http://voparis-rr.obspm.fr/tap\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "662ef537",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><i>Table length=249</i>\n",
       "<table id=\"table4426342752\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>table_name</th><th>access_url</th><th>ivoid</th></tr></thead>\n",
       "<thead><tr><th>object</th><th>object</th><th>object</th></tr></thead>\n",
       "<tr><td>thmsm.epn_core</td><td>http://dorotka.ufa.cas.cz//tap</td><td>ivo://iap/tap</td></tr>\n",
       "<tr><td>cpstasm.epn_core</td><td>http://dorotka.ufa.cas.cz//tap</td><td>ivo://iap/tap</td></tr>\n",
       "<tr><td>omega_maps.epn_core</td><td>https://idoc-dachs.ias.u-psud.fr/tap</td><td>ivo://idoc/tap</td></tr>\n",
       "<tr><td>omega_cubes.epn_core</td><td>https://idoc-dachs.ias.u-psud.fr/tap</td><td>ivo://idoc/tap</td></tr>\n",
       "<tr><td>gaia_dem.epn_core</td><td>https://idoc-dachs.ias.u-psud.fr/tap</td><td>ivo://idoc/tap</td></tr>\n",
       "<tr><td>eit_syn.epn_core</td><td>https://idoc-dachs.ias.u-psud.fr/tap</td><td>ivo://idoc/tap</td></tr>\n",
       "<tr><td>vvex.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>vims_satellites.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>tnosarecool.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>titan_profiles.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>spectro_planets.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>spectro_asteroids.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>planets.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>moonsprop.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>meteor_showers.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>m4ast.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>iks.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>hst_planeto.epn_core</td><td>http://voparis-tap-planeto.obspm.fr/tap</td><td>ivo://padc.obspm.planeto/tap</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>mp3c_best_epn.epn_core</td><td>https://dachs.oca.eu/tap</td><td>ivo://oca/tap</td></tr>\n",
       "<tr><td>uset_sunspot_group.epn_core</td><td>http://vo-tap.oma.be/tap</td><td>ivo://ksb-orb.dc/tap</td></tr>\n",
       "<tr><td>uset_sunspot_drawing.epn_core</td><td>http://vo-tap.oma.be/tap</td><td>ivo://ksb-orb.dc/tap</td></tr>\n",
       "<tr><td>rob_spoca_ch.epn_core</td><td>http://vo-tap.oma.be/tap</td><td>ivo://ksb-orb.dc/tap</td></tr>\n",
       "<tr><td>pretharsis.epn_core</td><td>https://dachs-vo.ipsl.fr/tap</td><td>ivo://geops.ipsl/tap</td></tr>\n",
       "<tr><td>exotopo.epn_core</td><td>https://dachs-vo.ipsl.fr/tap</td><td>ivo://geops.ipsl/tap</td></tr>\n",
       "<tr><td>mpc.epn_core</td><td>http://dc.zah.uni-heidelberg.de/tap</td><td>ivo://org.gavo.dc/tap</td></tr>\n",
       "<tr><td>vir.epn_core</td><td>http://vo-node1.iaps.inaf.it/tap</td><td>ivo://ia2.inaf.it/iaps/tap</td></tr>\n",
       "<tr><td>ndcs.epn_core</td><td>http://vo-node1.iaps.inaf.it/tap</td><td>ivo://ia2.inaf.it/iaps/tap</td></tr>\n",
       "<tr><td>cemls.epn_core</td><td>http://vo-node1.iaps.inaf.it/tap</td><td>ivo://ia2.inaf.it/iaps/tap</td></tr>\n",
       "<tr><td>vir.epn_core</td><td>http://vo-node1.iaps.inaf.it/tap</td><td>ivo://ia2.inaf.it/tap</td></tr>\n",
       "<tr><td>ndcs.epn_core</td><td>http://vo-node1.iaps.inaf.it/tap</td><td>ivo://ia2.inaf.it/tap</td></tr>\n",
       "<tr><td>cemls.epn_core</td><td>http://vo-node1.iaps.inaf.it/tap</td><td>ivo://ia2.inaf.it/tap</td></tr>\n",
       "<tr><td>climso.epn_core</td><td>http://climso.irap.omp.eu/__system__/tap/run/tap</td><td>ivo://ov-gso/climso</td></tr>\n",
       "<tr><td>amdadb.epn_core</td><td>http://cdpp-epntap.irap.omp.eu/__system__/tap/run/tap</td><td>ivo://cdpp/amda</td></tr>\n",
       "<tr><td>illu67p.epn_core</td><td>http://cdpp-epntap.irap.omp.eu/__system__/tap/run/tap</td><td>ivo://cdpp/illu67p</td></tr>\n",
       "<tr><td>transplanet.epn_core</td><td>http://cdpp-epntap.irap.omp.eu/__system__/tap/run/tap</td><td>ivo://cdpp/transplanet</td></tr>\n",
       "<tr><td>spicam.epn_core</td><td>http://vo.projet.latmos.ipsl.fr/tap</td><td>ivo://latmos.ipsl/tap</td></tr>\n",
       "<tr><td>abs_cs.epn_core</td><td>http://spectrum.iaa.es/tap</td><td>ivo://spectrum.iaa/tap</td></tr>\n",
       "</table></div>"
      ],
      "text/plain": [
       "<Table length=249>\n",
       "          table_name                                access_url                                 ivoid            \n",
       "            object                                    object                                   object           \n",
       "----------------------------- ----------------------------------------------------- ----------------------------\n",
       "               thmsm.epn_core                        http://dorotka.ufa.cas.cz//tap                ivo://iap/tap\n",
       "             cpstasm.epn_core                        http://dorotka.ufa.cas.cz//tap                ivo://iap/tap\n",
       "          omega_maps.epn_core                  https://idoc-dachs.ias.u-psud.fr/tap               ivo://idoc/tap\n",
       "         omega_cubes.epn_core                  https://idoc-dachs.ias.u-psud.fr/tap               ivo://idoc/tap\n",
       "            gaia_dem.epn_core                  https://idoc-dachs.ias.u-psud.fr/tap               ivo://idoc/tap\n",
       "             eit_syn.epn_core                  https://idoc-dachs.ias.u-psud.fr/tap               ivo://idoc/tap\n",
       "                vvex.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "     vims_satellites.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "         tnosarecool.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "      titan_profiles.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "     spectro_planets.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "   spectro_asteroids.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "             planets.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "           moonsprop.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "      meteor_showers.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "               m4ast.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "                 iks.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "         hst_planeto.epn_core               http://voparis-tap-planeto.obspm.fr/tap ivo://padc.obspm.planeto/tap\n",
       "                          ...                                                   ...                          ...\n",
       "       mp3c_best_epn.epn_core                              https://dachs.oca.eu/tap                ivo://oca/tap\n",
       "  uset_sunspot_group.epn_core                              http://vo-tap.oma.be/tap         ivo://ksb-orb.dc/tap\n",
       "uset_sunspot_drawing.epn_core                              http://vo-tap.oma.be/tap         ivo://ksb-orb.dc/tap\n",
       "        rob_spoca_ch.epn_core                              http://vo-tap.oma.be/tap         ivo://ksb-orb.dc/tap\n",
       "          pretharsis.epn_core                          https://dachs-vo.ipsl.fr/tap         ivo://geops.ipsl/tap\n",
       "             exotopo.epn_core                          https://dachs-vo.ipsl.fr/tap         ivo://geops.ipsl/tap\n",
       "                 mpc.epn_core                   http://dc.zah.uni-heidelberg.de/tap        ivo://org.gavo.dc/tap\n",
       "                 vir.epn_core                      http://vo-node1.iaps.inaf.it/tap   ivo://ia2.inaf.it/iaps/tap\n",
       "                ndcs.epn_core                      http://vo-node1.iaps.inaf.it/tap   ivo://ia2.inaf.it/iaps/tap\n",
       "               cemls.epn_core                      http://vo-node1.iaps.inaf.it/tap   ivo://ia2.inaf.it/iaps/tap\n",
       "                 vir.epn_core                      http://vo-node1.iaps.inaf.it/tap        ivo://ia2.inaf.it/tap\n",
       "                ndcs.epn_core                      http://vo-node1.iaps.inaf.it/tap        ivo://ia2.inaf.it/tap\n",
       "               cemls.epn_core                      http://vo-node1.iaps.inaf.it/tap        ivo://ia2.inaf.it/tap\n",
       "              climso.epn_core      http://climso.irap.omp.eu/__system__/tap/run/tap          ivo://ov-gso/climso\n",
       "              amdadb.epn_core http://cdpp-epntap.irap.omp.eu/__system__/tap/run/tap              ivo://cdpp/amda\n",
       "             illu67p.epn_core http://cdpp-epntap.irap.omp.eu/__system__/tap/run/tap           ivo://cdpp/illu67p\n",
       "         transplanet.epn_core http://cdpp-epntap.irap.omp.eu/__system__/tap/run/tap       ivo://cdpp/transplanet\n",
       "              spicam.epn_core                   http://vo.projet.latmos.ipsl.fr/tap        ivo://latmos.ipsl/tap\n",
       "              abs_cs.epn_core                            http://spectrum.iaa.es/tap       ivo://spectrum.iaa/tap"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "services = TapPlus(registry_url)\n",
    "job = services.launch_job(query)\n",
    "r1 = job.get_results()\n",
    "# lists all published EPN-TAP services:\n",
    "r1\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8fb09a31",
   "metadata": {},
   "source": [
    "**This lists the available services; table name and server url are printed for the first entry**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "b1aeebb7",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "          table_name         \n",
      "-----------------------------\n",
      "               thmsm.epn_core\n",
      "             cpstasm.epn_core\n",
      "          omega_maps.epn_core\n",
      "         omega_cubes.epn_core\n",
      "            gaia_dem.epn_core\n",
      "             eit_syn.epn_core\n",
      "                vvex.epn_core\n",
      "     vims_satellites.epn_core\n",
      "         tnosarecool.epn_core\n",
      "      titan_profiles.epn_core\n",
      "     spectro_planets.epn_core\n",
      "   spectro_asteroids.epn_core\n",
      "             planets.epn_core\n",
      "           moonsprop.epn_core\n",
      "      meteor_showers.epn_core\n",
      "               m4ast.epn_core\n",
      "                 iks.epn_core\n",
      "         hst_planeto.epn_core\n",
      "           exoplanet.epn_core\n",
      "                          ...\n",
      "       mp3c_best_epn.epn_core\n",
      "  uset_sunspot_group.epn_core\n",
      "uset_sunspot_drawing.epn_core\n",
      "        rob_spoca_ch.epn_core\n",
      "          pretharsis.epn_core\n",
      "             exotopo.epn_core\n",
      "                 mpc.epn_core\n",
      "                 vir.epn_core\n",
      "                ndcs.epn_core\n",
      "               cemls.epn_core\n",
      "                 vir.epn_core\n",
      "                ndcs.epn_core\n",
      "               cemls.epn_core\n",
      "              climso.epn_core\n",
      "              amdadb.epn_core\n",
      "             illu67p.epn_core\n",
      "         transplanet.epn_core\n",
      "              spicam.epn_core\n",
      "              abs_cs.epn_core\n",
      "Length = 249 rows\n",
      "thmsm.epn_core\n",
      "http://dorotka.ufa.cas.cz//tap\n"
     ]
    }
   ],
   "source": [
    "# accessing elt: \n",
    "print(r1['table_name'])\n",
    "print(r1['table_name'][0])\n",
    "print(r1['access_url'][0])\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "decc550b",
   "metadata": {},
   "source": [
    "### Sending a simple query to all services"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7409fabf",
   "metadata": {},
   "source": [
    "Here we're looking for values of an optional parameter. This will just print the table name if the parameter is not provided, or the number of values otherwise."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "dd193450",
   "metadata": {},
   "source": [
    "**The query is a bit long to execute, due to the large number of EPN-TAP services published in the registry (in January 2024)**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "8c5b5cab",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "thmsm.epn_core\n",
      "cpstasm.epn_core\n",
      "omega_maps.epn_core\n",
      "omega_cubes.epn_core\n",
      "eit_syn.epn_core\n",
      "vvex.epn_core\n",
      "vims_satellites.epn_core\n",
      "tnosarecool.epn_core\n",
      "titan_profiles.epn_core\n",
      "spectro_planets.epn_core\n",
      "spectro_asteroids.epn_core\n",
      "planets.epn_core\n",
      "moonsprop.epn_core\n",
      "meteor_showers.epn_core\n",
      "m4ast.epn_core\n",
      "iks.epn_core\n",
      "hst_planeto.epn_core\n",
      "exoplanet.epn_core\n",
      "dynastvo.epn_core\n",
      "crism_speclib.epn_core\n",
      "bdip.epn_core\n",
      "basecom.epn_core\n",
      "apis.epn_core\n",
      "hfc1t3.epn_core\n",
      "bass2000.epn_core\n",
      "mdisc.epn_core\n",
      "jasmin.epn_core\n",
      "hrsc3nd.epn_core\n",
      "usgs_wms.epn_core\n",
      "planmap.epn_core\n",
      "planetserver_crism.epn_core\n",
      "mars_craters_lagain.epn_core\n",
      "mars_craters.epn_core\n",
      "sbnaf.epn_core\n",
      "vcd.epn_core\n",
      "mcd.epn_core\n",
      "mars_dust.epn_core\n",
      "sshade_spectra.epn_core\n",
      "mp3c_best_epn.epn_core\n",
      "pvol.epn_core\n",
      "rwcalerts.epn_core\n",
      "lofar_jupiter_raw.epn_core\n",
      "lofar_jupiter.epn_core\n",
      "ionosondes.epn_core\n",
      "cassini_rpws.epn_core\n",
      "mp3c_best_epn.epn_core\n",
      "uset_sunspot_group.epn_core\n",
      "uset_sunspot_drawing.epn_core\n",
      "rob_spoca_ch.epn_core\n",
      "pretharsis.epn_core\n",
      "exotopo.epn_core\n",
      "mpc.epn_core\n",
      "climso.epn_core\n",
      "amdadb.epn_core\n",
      "illu67p.epn_core\n",
      "transplanet.epn_core\n",
      "spicam.epn_core\n",
      "abs_cs.epn_core\n"
     ]
    }
   ],
   "source": [
    "for ii in range((r1['table_name']).shape[0]):\n",
    "\tStable_name = r1['table_name'][ii]\n",
    "\tService_url = r1['access_url'][ii]\n",
    "\tquery2 = \"SELECT distinct spatial_origin FROM \" + Stable_name\n",
    "\tanswers = TapPlus(Service_url)\n",
    "\t# remove one uncorrect entry\n",
    "\tif r1['access_url'][ii] == 'http://vo-node1.iaps.inaf.it/tap': \n",
    "\t\tcontinue\n",
    "\t# skip the largest server\n",
    "\tif r1['access_url'][ii] == 'https://vo-pds-ppi.igpp.ucla.edu/tap': \n",
    "\t\tcontinue\n",
    "\ttry:\n",
    "\t\tjob2 = answers.launch_job(query2)\n",
    "\t\tr2 = job2.get_results()\n",
    "\t\tStable_name, r2[0][0]\n",
    "\texcept:\n",
    "\t\tprint(Stable_name) # command required here\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "239cb42d",
   "metadata": {},
   "source": [
    "### Sending a more realistic query to all services"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "268c6fe5",
   "metadata": {},
   "outputs": [],
   "source": [
    "# first define a result table \n",
    "# actually doesn't work - all grouped in first column later\n",
    "d = {\n",
    "    'service_title': [''],\n",
    "    'granule_uid': [''],\n",
    "    'time_min': [0.],\n",
    "    'solar_longitude_min': [0.],\n",
    "    'solar_longitude_max': [0.]\n",
    "}\n",
    "import pandas as pd\n",
    "df = pd.DataFrame(d)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "01446a07",
   "metadata": {},
   "source": [
    "Then sends a more realistic query to all services and concatenate results - we're asking for any data during a specific season on Mars, whatever the year"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "a9b15fa1",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "thmsm.epn_core\n",
      "cpstasm.epn_core\n",
      "omega_maps.epn_core\n",
      "omega_cubes.epn_core\n",
      "gaia_dem.epn_core\n",
      "eit_syn.epn_core\n",
      "vvex.epn_core 0\n",
      "vims_satellites.epn_core\n",
      "tnosarecool.epn_core\n",
      "titan_profiles.epn_core\n",
      "spectro_planets.epn_core\n",
      "spectro_asteroids.epn_core\n",
      "planets.epn_core\n",
      "moonsprop.epn_core\n",
      "meteor_showers.epn_core 2\n",
      "m4ast.epn_core\n",
      "iks.epn_core\n",
      "hst_planeto.epn_core\n",
      "exoplanet.epn_core\n",
      "dynastvo.epn_core\n",
      "crism_speclib.epn_core\n",
      "bdip.epn_core\n",
      "basecom.epn_core\n",
      "apis.epn_core\n",
      "hfc1t3.epn_core\n",
      "hfc1ar_majtrim.epn_core\n",
      "hfc1ar.epn_core\n",
      "bass2000.epn_core\n",
      "mdisc.epn_core\n",
      "jasmin.epn_core\n",
      "hrsc3nd.epn_core 169\n",
      "usgs_wms.epn_core\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: In merged column 'granule_uid' the 'description' attribute does not match (Internal table row index, which must be unique within the table. Can be alphanumeric. != Internal table row index Unique ID in data service, also in v2. Can be alphanumeric.).  Using Internal table row index Unique ID in data service, also in v2. Can be alphanumeric. for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'ucd' types <class 'str'> and <class 'str'>, choosing ucd='time.start' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: In merged column 'time_min' the 'description' attribute does not match (Acquisition start time (in JD), as UTC at time_refposition != Acquisition start time (in JD)).  Using Acquisition start time (in JD) for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: In merged column 'solar_longitude_min' the 'description' attribute does not match (Min Solar longitude Ls (location on orbit / season) != Solar longitude (minimum)).  Using Solar longitude (minimum) for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: In merged column 'solar_longitude_max' the 'description' attribute does not match (Max Solar longitude Ls (location on orbit / season) != Solar longitude (maximum)).  Using Solar longitude (maximum) for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'utype' types <class 'str'> and <class 'str'>, choosing utype='ivo://vopdc.obspm/std/epncore#schema-2.0' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'description' types <class 'str'> and <class 'str'>, choosing description='EPN-TAP access to HRSC level 3 panchromatic image data' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "planmap.epn_core\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: UnitsWarning: The unit 'pix' has been deprecated in the VOUnit standard. [astropy.units.format.utils]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "planetserver_crism.epn_core 1244\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: In merged column 'granule_uid' the 'description' attribute does not match (Internal table row index Unique ID in data service, also in v2. Can be alphanumeric. != Internal table row index, which must be unique within the table. Can be alphanumeric.).  Using Internal table row index, which must be unique within the table. Can be alphanumeric. for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'ucd' types <class 'str'> and <class 'str'>, choosing ucd='time.start;obs' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: In merged column 'time_min' the 'description' attribute does not match (Acquisition start time (in JD) != Acquisition start time (in JD), as UTC at time_refposition).  Using Acquisition start time (in JD), as UTC at time_refposition for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'ucd' types <class 'str'> and <class 'str'>, choosing ucd='pos.posAng;pos.heliocentric;stat.min' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: In merged column 'solar_longitude_min' the 'description' attribute does not match (Solar longitude (minimum) != Min Solar longitude Ls (location on orbit / season)).  Using Min Solar longitude Ls (location on orbit / season) for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'ucd' types <class 'str'> and <class 'str'>, choosing ucd='pos.posAng;pos.heliocentric;stat.max' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: In merged column 'solar_longitude_max' the 'description' attribute does not match (Solar longitude (maximum) != Max Solar longitude Ls (location on orbit / season)).  Using Max Solar longitude Ls (location on orbit / season) for merged output [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'utype' types <class 'str'> and <class 'str'>, choosing utype='ivo://ivoa.net/std/epntap#table-2.0' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'description' types <class 'str'> and <class 'str'>, choosing description='PlanetServer2/CRISM data: http://access.planetserver.eu/' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "mars_craters_lagain.epn_core\n",
      "mars_craters.epn_core\n",
      "sbnaf.epn_core\n",
      "vcd.epn_core\n",
      "mcd.epn_core 0\n",
      "mars_dust.epn_core\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'ucd' types <class 'str'> and <class 'str'>, choosing ucd='pos.ecliptic.lon;pos.heliocentric;stat.min' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'ucd' types <class 'str'> and <class 'str'>, choosing ucd='pos.ecliptic.lon;pos.heliocentric;stat.max' [astropy.utils.metadata]\n",
      "WARNING: MergeConflictWarning: Cannot merge meta key 'description' types <class 'str'> and <class 'str'>, choosing description='The database contains some outputs of the MCD for all the scenarios\\navailable, at different positions of latitude and longitude. These MCD\\noutputs are provided as Votables containing profiles of temperatures,\\npressures, density and abundances of O2 and O3 for altitudes between 0\\nand 249.5 km from the surface.' [astropy.utils.metadata]\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "sshade_spectra.epn_core\n",
      "mp3c_best_epn.epn_core\n",
      "pvol.epn_core\n",
      "rwcalerts.epn_core\n",
      "lofar_jupiter_raw.epn_core\n",
      "lofar_jupiter.epn_core\n",
      "ionosondes.epn_core\n",
      "lars.epn_core\n",
      "gris.epn_core\n",
      "chrotel.epn_core\n",
      "cassini_rpws.epn_core\n",
      "soir.epn_core 0\n",
      "nomad.epn_core 371\n",
      "gem_mars.epn_core 2000\n",
      "mp3c_best_epn.epn_core\n",
      "uset_sunspot_group.epn_core\n",
      "uset_sunspot_drawing.epn_core\n",
      "rob_spoca_ch.epn_core\n",
      "pretharsis.epn_core\n",
      "exotopo.epn_core\n",
      "mpc.epn_core\n",
      "climso.epn_core\n",
      "amdadb.epn_core\n",
      "illu67p.epn_core\n",
      "transplanet.epn_core\n",
      "spicam.epn_core 118\n",
      "abs_cs.epn_core\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING: MergeConflictWarning: Cannot merge meta key 'description' types <class 'str'> and <class 'str'>, choosing description='The database contains profiles of atmospheric CO2 density and\\ntemperature derived thereof, and ozone concentration profiles, all\\nderived from the first year(s) of SPICAM UV stellar occultation\\nobservations (Forget et al. 2009, Lebonnois et al. 2006). The database\\nalso includes the full dataset of four Martian years of aerosol\\nextinction profiles from SPICAM UV solar occultations (Määttänen et\\nal. 2013).' [astropy.utils.metadata]\n"
     ]
    }
   ],
   "source": [
    "for ii in range((r1['table_name']).shape[0]):\n",
    "\tStable_name = r1['table_name'][ii]\n",
    "\tService_url = r1['access_url'][ii]\n",
    "\tquery2 = \"SELECT * FROM \" + Stable_name + \" WHERE target_name = 'Mars' and solar_longitude_min >= 30 and solar_longitude_max <= 40\"\n",
    "\tanswers = TapPlus(Service_url)\n",
    "\t# remove one uncorrect entry\n",
    "\tif r1['access_url'][ii] == 'http://vo-node1.iaps.inaf.it/tap': \n",
    "\t\tcontinue\n",
    "\t# skip the largest server\n",
    "\tif r1['access_url'][ii] == 'https://vo-pds-ppi.igpp.ucla.edu/tap': \n",
    "\t\tcontinue\n",
    "\ttry:\n",
    "\t\tjob2 = answers.launch_job(query2)\n",
    "\t\tr2 = job2.get_results()\n",
    "\t\tprint(Stable_name , r2['service_title'].shape[0] )\n",
    "\t# cumulative results\n",
    "\t\tdf = vstack([df, r2['service_title','granule_uid','time_min', 'solar_longitude_min','solar_longitude_max']], join_type='outer')\n",
    "\texcept:\n",
    "\t\tprint(Stable_name ) # a command is required"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "1b0894e8",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><i>Table length=3905</i>\n",
       "<table id=\"table4426514496\" class=\"table-striped table-bordered table-condensed\">\n",
       "<thead><tr><th>col0</th><th>service_title</th><th>granule_uid</th><th>time_min</th><th>solar_longitude_min</th><th>solar_longitude_max</th></tr></thead>\n",
       "<thead><tr><th></th><th></th><th></th><th>d</th><th>deg</th><th>deg</th></tr></thead>\n",
       "<thead><tr><th>object[5]</th><th>object</th><th>object</th><th>float64</th><th>float64</th><th>float64</th></tr></thead>\n",
       "<tr><td>.. 0.0</td><td>--</td><td>--</td><td>--</td><td>--</td><td>--</td></tr>\n",
       "<tr><td>-- .. --</td><td>meteor_showers</td><td>MAR_6WU_2018_MULT_0::v1.0</td><td>2458161.863194444</td><td>34.151</td><td>34.151</td></tr>\n",
       "<tr><td>-- .. --</td><td>meteor_showers</td><td>MAR_6WU_2021_MULT_0::v1.0</td><td>2459537.4125</td><td>35.389</td><td>35.389</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0376_0000.nd3.09</td><td>2453133.87744213</td><td>30.15999984741211</td><td>30.15999984741211</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0380_0000.nd3.12</td><td>2453135.06644676</td><td>30.700000762939453</td><td>30.700000762939453</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0380_0001.nd3.11</td><td>2453135.07356481</td><td>30.700000762939453</td><td>30.700000762939453</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0383_0000.nd3.09</td><td>2453135.95583333</td><td>31.110000610351562</td><td>31.110000610351562</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0389_0000.nd3.11</td><td>2453137.67203704</td><td>31.889999389648438</td><td>31.889999389648438</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0394_0000.nd3.07</td><td>2453139.07164352</td><td>32.529998779296875</td><td>32.529998779296875</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0394_0009.nd3.07</td><td>2453139.07917824</td><td>32.540000915527344</td><td>32.540000915527344</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0397_0000.nd3.08</td><td>2453139.91438657</td><td>32.91999816894531</td><td>32.91999816894531</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0401_0000.nd3.14</td><td>2453141.03354167</td><td>33.43000030517578</td><td>33.43000030517578</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0401_0001.nd3.14</td><td>2453141.04221065</td><td>33.43000030517578</td><td>33.43000030517578</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0405_0000.nd3.10</td><td>2453142.15415509</td><td>33.939998626708984</td><td>33.939998626708984</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0411_0000.nd3.09</td><td>2453143.83871528</td><td>34.70000076293945</td><td>34.70000076293945</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0416_0000.nd3.09</td><td>2453145.23615741</td><td>35.33000183105469</td><td>35.33000183105469</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0420_0000.nd3.08</td><td>2453146.358125</td><td>35.84000015258789</td><td>35.84000015258789</td></tr>\n",
       "<tr><td>-- .. --</td><td>hrsc3nd</td><td>h0420_0001.nd3.08</td><td>2453146.36568287</td><td>35.849998474121094</td><td>35.849998474121094</td></tr>\n",
       "<tr><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td><td>...</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0436A2_o3dens_SPICAMmeasurements</td><td>2453150.43312</td><td>37.89</td><td>37.89</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0409A6_o3dens_SPICAMmeasurements</td><td>2453142.86246</td><td>34.47</td><td>34.47</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0433A5_o3dens_SPICAMmeasurements</td><td>2453149.59333</td><td>37.51</td><td>37.51</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0395A1_o3dens_SPICAMmeasurements</td><td>2453138.93386</td><td>32.68</td><td>32.68</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0388A1_o3dens_SPICAMmeasurements</td><td>2453136.96869</td><td>31.78</td><td>31.78</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0426A1_o3dens_SPICAMmeasurements</td><td>2453147.45563</td><td>36.55</td><td>36.55</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0405A1_o3dens_SPICAMmeasurements</td><td>2453141.54687</td><td>33.87</td><td>33.87</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0392A3_temp_SPICAMmeasurements</td><td>2453138.10524</td><td>32.3</td><td>32.3</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0426A1_temp_SPICAMmeasurements</td><td>2453147.34952</td><td>36.5</td><td>36.5</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0398A1_temp_SPICAMmeasurements</td><td>2453139.64026</td><td>33.0</td><td>33.0</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit2888A1_temp_SPICAMmeasurements</td><td>2453836.75638</td><td>37.6</td><td>37.6</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0395A1_temp_SPICAMmeasurements</td><td>2453138.98209</td><td>32.7</td><td>32.7</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0409A6_temp_SPICAMmeasurements</td><td>2453142.93734</td><td>34.5</td><td>34.5</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0436A2_temp_SPICAMmeasurements</td><td>2453150.44868</td><td>37.9</td><td>37.9</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0388A1_temp_SPICAMmeasurements</td><td>2453137.01027</td><td>31.8</td><td>31.8</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0419A1_temp_SPICAMmeasurements</td><td>2453145.58252</td><td>35.7</td><td>35.7</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0405A1_temp_SPICAMmeasurements</td><td>2453141.61723</td><td>33.9</td><td>33.9</td></tr>\n",
       "<tr><td>-- .. --</td><td>spicam</td><td>orbit0433A5_temp_SPICAMmeasurements</td><td>2453149.56235</td><td>37.5</td><td>37.5</td></tr>\n",
       "</table></div>"
      ],
      "text/plain": [
       "<Table length=3905>\n",
       "   col0   service_title               granule_uid                   time_min     solar_longitude_min solar_longitude_max\n",
       "                                                                       d                 deg                 deg        \n",
       "object[5]     object                     object                     float64            float64             float64      \n",
       "--------- -------------- ------------------------------------- ----------------- ------------------- -------------------\n",
       "   .. 0.0             --                                    --                --                  --                  --\n",
       " -- .. -- meteor_showers             MAR_6WU_2018_MULT_0::v1.0 2458161.863194444              34.151              34.151\n",
       " -- .. -- meteor_showers             MAR_6WU_2021_MULT_0::v1.0      2459537.4125              35.389              35.389\n",
       " -- .. --        hrsc3nd                     h0376_0000.nd3.09  2453133.87744213   30.15999984741211   30.15999984741211\n",
       " -- .. --        hrsc3nd                     h0380_0000.nd3.12  2453135.06644676  30.700000762939453  30.700000762939453\n",
       " -- .. --        hrsc3nd                     h0380_0001.nd3.11  2453135.07356481  30.700000762939453  30.700000762939453\n",
       " -- .. --        hrsc3nd                     h0383_0000.nd3.09  2453135.95583333  31.110000610351562  31.110000610351562\n",
       " -- .. --        hrsc3nd                     h0389_0000.nd3.11  2453137.67203704  31.889999389648438  31.889999389648438\n",
       " -- .. --        hrsc3nd                     h0394_0000.nd3.07  2453139.07164352  32.529998779296875  32.529998779296875\n",
       " -- .. --        hrsc3nd                     h0394_0009.nd3.07  2453139.07917824  32.540000915527344  32.540000915527344\n",
       " -- .. --        hrsc3nd                     h0397_0000.nd3.08  2453139.91438657   32.91999816894531   32.91999816894531\n",
       " -- .. --        hrsc3nd                     h0401_0000.nd3.14  2453141.03354167   33.43000030517578   33.43000030517578\n",
       " -- .. --        hrsc3nd                     h0401_0001.nd3.14  2453141.04221065   33.43000030517578   33.43000030517578\n",
       " -- .. --        hrsc3nd                     h0405_0000.nd3.10  2453142.15415509  33.939998626708984  33.939998626708984\n",
       " -- .. --        hrsc3nd                     h0411_0000.nd3.09  2453143.83871528   34.70000076293945   34.70000076293945\n",
       " -- .. --        hrsc3nd                     h0416_0000.nd3.09  2453145.23615741   35.33000183105469   35.33000183105469\n",
       " -- .. --        hrsc3nd                     h0420_0000.nd3.08    2453146.358125   35.84000015258789   35.84000015258789\n",
       " -- .. --        hrsc3nd                     h0420_0001.nd3.08  2453146.36568287  35.849998474121094  35.849998474121094\n",
       "      ...            ...                                   ...               ...                 ...                 ...\n",
       " -- .. --         spicam orbit0436A2_o3dens_SPICAMmeasurements     2453150.43312               37.89               37.89\n",
       " -- .. --         spicam orbit0409A6_o3dens_SPICAMmeasurements     2453142.86246               34.47               34.47\n",
       " -- .. --         spicam orbit0433A5_o3dens_SPICAMmeasurements     2453149.59333               37.51               37.51\n",
       " -- .. --         spicam orbit0395A1_o3dens_SPICAMmeasurements     2453138.93386               32.68               32.68\n",
       " -- .. --         spicam orbit0388A1_o3dens_SPICAMmeasurements     2453136.96869               31.78               31.78\n",
       " -- .. --         spicam orbit0426A1_o3dens_SPICAMmeasurements     2453147.45563               36.55               36.55\n",
       " -- .. --         spicam orbit0405A1_o3dens_SPICAMmeasurements     2453141.54687               33.87               33.87\n",
       " -- .. --         spicam   orbit0392A3_temp_SPICAMmeasurements     2453138.10524                32.3                32.3\n",
       " -- .. --         spicam   orbit0426A1_temp_SPICAMmeasurements     2453147.34952                36.5                36.5\n",
       " -- .. --         spicam   orbit0398A1_temp_SPICAMmeasurements     2453139.64026                33.0                33.0\n",
       " -- .. --         spicam   orbit2888A1_temp_SPICAMmeasurements     2453836.75638                37.6                37.6\n",
       " -- .. --         spicam   orbit0395A1_temp_SPICAMmeasurements     2453138.98209                32.7                32.7\n",
       " -- .. --         spicam   orbit0409A6_temp_SPICAMmeasurements     2453142.93734                34.5                34.5\n",
       " -- .. --         spicam   orbit0436A2_temp_SPICAMmeasurements     2453150.44868                37.9                37.9\n",
       " -- .. --         spicam   orbit0388A1_temp_SPICAMmeasurements     2453137.01027                31.8                31.8\n",
       " -- .. --         spicam   orbit0419A1_temp_SPICAMmeasurements     2453145.58252                35.7                35.7\n",
       " -- .. --         spicam   orbit0405A1_temp_SPICAMmeasurements     2453141.61723                33.9                33.9\n",
       " -- .. --         spicam   orbit0433A5_temp_SPICAMmeasurements     2453149.56235                37.5                37.5"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8e2db595",
   "metadata": {},
   "source": [
    "Then plot cumulative results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "66f9747c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAvwAAAEmCAYAAAAeFURbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAABQxElEQVR4nO3de1wU9f4/8NcCgiIXBRVCSBTTvGeJRunxhmXeKLKOl8jSfmZWX+3ivVN28ZKpqfnVTlZWluLJwLKOnKwEsVJB40jqVw3hiHg7KnIRg1zm98e06MAyOwuzzM7M6/l48FiZfbt+lpHd937m/Xl/LIIgCCAiIiIiIkPy0HoARERERETkOkz4iYiIiIgMjAk/EREREZGBMeEnIiIiIjIwJvxERERERAbGhJ+IiIiIyMCY8BMRERERGRgTfiIiIiIiA/PSegCuVllZidOnT8Pf3x8Wi0Xr4RARERERqUIQBJSUlCAsLAweHrXP4xs+4T99+jQiIiK0HgYRERERkUvk5+cjPDy81vsNn/D7+/sDEH8QAQEBGo+GiIiIiEgdxcXFiIiIqMp3a2P4hN9WxhMQEMCEn4iIiIgMx1HZOhftEhEREREZGBN+IiIiIiIDY8JPpLXSUuCBB4Du3cXb0lKtR0REREQG4jYJ/6JFi2CxWDB9+vSqY4IgYP78+QgLC0OTJk0wYMAAHDp0SLtBEqmtd2/A3x/YuhXIzhZv/f3F40RERPVRUQGsWAE8+6x4W1Gh9YhII26R8GdkZOC9995D9+7dJceXLFmC5cuXY/Xq1cjIyEBoaCiGDBmCkpISjUZaR5zBJXt69wYyMuzfl5HBpJ+IiOpu5kygSRPgueeA1avF2yZNxONkOpon/KWlpRg/fjzWrVuH5s2bVx0XBAErVqzAvHnzEB8fj65du+Ljjz9GWVkZNm7cqOGIncQZXLKntLT2ZN8mI4MfDomIyHkzZwJvvQVUVkqPV1aKx5n0m47mCf/TTz+N4cOHIzY2VnI8NzcXZ8+exT333FN1zMfHB/3798dPP/1U6+OVl5ejuLhY8qUZzuBSbf76V3XjiIiIALFsZ9ky+Zhly1jeYzKaJvyJiYk4cOAAFi1aVOO+s2fPAgBCQkIkx0NCQqrus2fRokUIDAys+tJsl13O4JKcH35QN85IrFYgNRXYtEm8tVq1HhERkX68807Nmf3qKivFODINzRL+/Px8TJs2DZ9++ikaN25ca1z1jQQEQZDdXGDOnDkoKiqq+srPz1dtzE5JSFA3jozl99/VjTOKpCQgMhIYOBAYN068jYwUjxNRTfyATNXt3q1uHBmCZgn//v37cf78edxxxx3w8vKCl5cX0tLSsGrVKnh5eVXN7FefzT9//nyNWf8b+fj4VO2qq+nuujk56sYRGV1SEjB6NHDqlPR4QYF4nEk/kRQ/IJM9TZqoG0eGoFnCP3jwYGRnZyMrK6vqq1evXhg/fjyysrLQrl07hIaGYseOHVV/p6KiAmlpabjrrru0GrZykZHqxhEZmdUKTJsGCELN+2zHpk/n7CWRDT8gU218fdWNI0Pw0uof9vf3R9euXSXHmjZtiuDg4Krj06dPx8KFC3HLLbfglltuwcKFC+Hr64tx48ZpMWTnTJ4MbNumLI7I7NLTayYuNxIEID9fjBswoMGGReSWHH1AtljED8hxcYCnZ4MPjzR25Ii6cWQImiX8SsycORNXr17F1KlTUVhYiD59+uDbb7+Fv7+/1kNzTOleAXrbU4DIFc6cUTeOjMtqFT/4nTkD3HQT0K+f+ZJafkAmOTLrHOsUR4bgVgl/amqq5HuLxYL58+dj/vz5moynXm66Sd04IiPj7wspkZQkzmzfmOyGhwMrVwLx8dqNq6HxAzLJGToU+PFHZXFkGpr34Tesfv0APz/5GH9/MY7I7Pr1ExO32macLBYgIoK/L2bGmvXr+AGZ5Ozbp24cGQITflexWoGyMvmYK1e4CJEIEEsyVq4U/1w96bd9v2KF+Uo3SMRF3VL8gExycnPVjSNDYMLvKmvWKNv4Ys2ahhkPkbuLjwe2bAFat5YeDw8Xj5upZIOknKlZNwN+QCY5StuRa9W2nDTBhN9V2IefyHnx8UBeHrBzJ7Bxo3ibm8tk3+xYs15TfDzw4ouAR7W3cQ8P8Th/Z8zr7rvVjSNDcKtFu4ZSfZayvnFEZuHpyc4iJMWa9ZqSkoClS2uWOVmt4vE772TSb1bVPwTWN44MgWfbVX76Sd04IiKzYs26lNyaBhszrWkgqf/7P3XjyBCY8LtKXp66cUREZsWadSmuaSA5V66oG0eGwITfVaKi1I0jIjIzLuq+jmsaSI6jluDOxpEhsIbfVTZsEPvsK4kjIiLH4uOBuDjutNuqlbpxZCyjRgFbtyqLI9Ngwu8qfn5AdDSQkVF7THQ0P2ETETmDi7qJ5J0/r24cGQJLelxp377aS3aiorjLHREROY8JHcnZvFndODIEJvyulJQEnDhh/74TJ8y1FTwREamDbUpJDpuGkB1M+F2FbdOIiMgV2KaU5HCnXbKDCb+rsG0aERG5AtuUkpyJE9WNI0Ngwu8qbJtGRESuwjalVJuCAnXjyBDYpcdVWGNJVDelpUBCApCTIy5u37CB3ayI7GGbUrKntrWDdY0jQ2DC7yq2GsuCAvt1/BaLeD9rLImu691b2so2O1vczyI6ml2tiOxhm1Kq7upVdePIEFjS4yqssSRyTvVk/0YZGeL9REQkr3NndePIEJjwu1J8vLiTXfUZfkEQj7PGkkhUWiq/SR0g3l9a2jDjISLSK6UTiZxwNBUm/K40cybw5Zf27/vyS/F+M7JagdRUYNMm8ZatSSkhQd04IiKzatZM3TgyBCb8rlJRASxfLh+zfLkYZyZJSUBkJDBwIDBunHgbGclNyMwuJ0fdODIuThgQyattf4a6xpEhMOF3lTVrHL8RWa1inFkkJQGjR9fcn6CgQDzOpN+8oqLUjSNj4oQBkWP//a+6cWQITPhdhTOWUnI7D9uOcedh89qwQd04Mh5OGNjHKx5U3Z496saRITDhd5U2bdSN0zvuPExy/PzE1ptyoqPZj9+sOGFgH694kD3l5erGkSEw4aeGwZ2HyZF9+2pP+tmH39w4YVATr3hQbYKC1I0jQ2DC7yq5uerG6R13HiYl9u0DSkqA++8HunUTb0tKmOybHScMpBxd8RAEc17xIBEX7ZId3GnXVfgLJ8Wdh0kpPz8gOVnrUZA74YSBlKMrHsD1Kx7chdd8lO5Xwn1NTIUz/K7Sp4+6cXrHnYeJqK5sEwZyIiLMM2FQUKBuHBlLt27qxpEhMOF3lYgIdeOMID4e2LIFaN1aejw8XDxulp2HL13SegRE+uLpCYwdKx8zZox5JgzYdpHkPPqounFkCEz4XaVfPyA4WD7GTDNSNvHxQF4esHMnsHGjeJuba55kHwBiYrQeAZG+WK1i20k5iYnmqVlv2VLdODKWRo3UjSNDYA2/q3z5JXDxonyMWUtYPD3NXVd67JjWIyDSF9asS4WGqhtHxnL2rLpxZAiazvCvXbsW3bt3R0BAAAICAhATE4Pt27dX3V9aWopnnnkG4eHhaNKkCTp16oS1a9dqOGKFbB0U5AQHA3FxDTMe0q8PP9R6BETaY806kXIs+SI7NE34w8PDsXjxYmRmZiIzMxODBg1CXFwcDh06BAB47rnnkJKSgk8//RRHjhzBc889h2effRZffvmllsN2TMls1MWL5uoZTXXj7a31CIi0xwRG6vx5dePIWFjyRXZomvCPHDkSw4YNQ4cOHdChQwcsWLAAfn5+2PPnds8///wzJkyYgAEDBiAyMhKTJ09Gjx49kJmZqeWwHWPPaFKLWRIYIjlMYKRatVI3joyFJV9kh9ss2rVarUhMTMSVK1cQ8+eixr59++Krr75CQUEBBEHAzp07cezYMdx77721Pk55eTmKi4slXw2OPaNJLWZJYIjkVO/sVd84IiNTunjdLIvcCYAbJPzZ2dnw8/ODj48PpkyZguTkZHTu3BkAsGrVKnTu3Bnh4eHw9vbG0KFDsWbNGvTt27fWx1u0aBECAwOrviK0aHvZpYu6cWReTGCI2Ie/Opb0kJzUVHXjyBA0T/g7duyIrKws7NmzB0899RQmTJiAw4cPAxAT/j179uCrr77C/v37sWzZMkydOhXfffddrY83Z84cFBUVVX3l5+c31FO5btAgdePIvMySwBDJsW3cZ7HY37jPYjFX1zOW9JAcpWXP7l4eTarSvC2nt7c32rdvDwDo1asXMjIysHLlSqxYsQJz585FcnIyhg8fDgDo3r07srKysHTpUsTGxtp9PB8fH/j4+DTY+O06fVrdODIvsyQwRI7YNu6bNk3aFCE8XEz2zbSXB5EcX19148gQNE/4qxMEAeXl5fjjjz/wxx9/wMNDehHC09MTlZWVGo1OobAwZbuphoW5fixEREYRHy+2M05PF5se3HSTeBXMbB+MWdJDcmJixL2AlMSRaWia8M+dOxf33XcfIiIiUFJSgsTERKSmpiIlJQUBAQHo378/ZsyYgSZNmqBNmzZIS0vDJ598guXLl2s5bMfS0hzvsmuLM5vSUiAhAcjJAaKigA0bAD8/rUdFRHpgtTLZB9gYgoicpmnCf+7cOSQkJODMmTMIDAxE9+7dkZKSgiFDhgAAEhMTMWfOHIwfPx6XLl1CmzZtsGDBAkyZMkXLYTsWFASEhADnztUeExIixplJ795ARsb177OzAX9/IDoa2LdPu3ERkftLSrJfzrNypfnKeWyLmAsKAEGoeb/FIt7PNUDm9Gdrc9XiyBA0Tfg/+OAD2ftDQ0Oxfv36BhqNys6eBZo2BcrKat7n62u+La2rJ/s3ysgQ72fST0T2JCUBo0fXTG4LCsTjW7aYK+m3LWIePbrmfbZFzWZaxExSStcxar3ekRqU5l16DGvmTPvJPiAenzmzYcejpdLS2pN9m4wMMY6I6EZWqzizb28m23Zs+nRz9hS3d5U4KMh8H4BIatcuZXEnT7p2HORWmPC7QkUF8NZb8jFvvSXGmUFCgrpxRGQe6enSMp7qBAHIzxfjzMJ2xePixZr32TtG5lFRIa5xUaKkxLVjIbfChN8VFixQN07vcnLUjSMi81CavCiN0zu5Kx6AWNJj1iseJJZ6KfVnS3QyByb8ruBodt/ZOL2LilI3jojMgx1ppHjFg+Qoacdps2GD68ZBbsfpRbt5eXlIT09HXl4eysrK0LJlS/Ts2RMxMTFo3LixK8aoP0pLdcxS0rNhg9iNR0kcEdGN2JFGilc8SE5tV36qa9qULbFNRnHCv3HjRqxatQr79u1Dq1at0Lp1azRp0gSXLl1CTk4OGjdujPHjx2PWrFlo06aNK8fs/lq0kG/JeWOcGfj5ia035RbuRkfzxYeIarqxI43FIk1ozNiRplUrdePIWJS2+x440LXjILejqKTn9ttvx/Lly/HII48gLy8PZ8+exf79+7F7924cPnwYxcXF+PLLL1FZWYlevXrh888/d/W43VtWlrpxRrBvn5jU28M+/EQkJz5e7DzTurX0eHg4O9IQ3chDYaW20jgyDEUz/K+//jqGDx9e6/0+Pj4YMGAABgwYgDfeeAO5ubmqDVCXWrZUN84o9u3jTrtEVDfx8UBcHHfaPX9e3TgylqZN1Y0jw1CU8Msl+9W1aNECLcxSqlIbpYul0tOBAQNcOhS34+cHJCdrPQoi0iNPT/O9ZlbHRcwkR8l6OWfiyDB4TccVuKiKiIhcoV8/IDhYPiY42DyLmIlIEcUJ/549e9ClSxdERUUhMTHRlWPSPy6qIiIiooZ28KC6cWQYihP+qVOn4vXXX8d3332HJ554AuXl5a4cFxlVRYXYUePZZ8Vbs7QmJSJSQ3q64910L15kH36zUrrTMndkNh3FCf+lS5cQHh6OkJAQlJeXo6yszJXj0jcuqrJv5kzA1xd47jlg9Wrx1tdXPE5ERI6xZJTkNGqkbhwZhuI+/PPmzcMjjzyC5s2b47HHHkPz5s1dOS5946KqmmbOtL+zsNV6/fiSJQ07JiIiveH7C8lp1w44fFhZHJmKRRCUbssG5Ofno7i4GF26dHHlmFRVXFyMwMBAFBUVISAgoGH+UasViIx0vDNkbq45WspVVIgz+VZr7TGenkBZGeDt3XDj0optsyAllP96Ehmf1cq2nHx/ITlTpgB//7vjuCefBN591/XjIZdTmuc61aUnIiJCV8m+Zmw7QwI1kzsz7gy5Zo18sg+I969Z0zDjISL9SUoSE92BA4Fx48TbyEjxuJnw/YWI6oBtOV0lPh548cWau9l5eIjHzbQzZE6OunFkXFYrkJoKbNok3jr6oEjmkJQEjB4NnDolPV5QIB43W9LPnYepNoWF6saRYTDhd5WkJGDp0poJi9UqHjfTG1RUlLpxZEycwSV7rFZg2jT75Su2Y9Onm+/DYXw8kJcH7NwJbNwo3ubmMtk3OyX1+87EkWEw4XcFuTcoGzO9QU2d6vjysqenGEfmxBlcqk16es3/FzcSBCA/n20oiQBA6VrFhlrTSG6DCb8r8A1KytsbGDFCPmbECHMs2KWaOINLctiG0j5eESN7HL3XOhtHhsGE3xX4BiVltQL798vHHDjAhM6s+AGZ5HDn8pp4RYxq46Ww27rSODIMp8/4lStXsHjxYnz//fc4f/48KisrJfefOHFCtcHpFvskSzlK6IDrCd2AAQ0yJHIj/IBMpJyjK2IWi3hFLC6OnXrMKC9P3TgyDKcT/ieeeAJpaWlISEjATTfdBIszPcXNok8fdeP0jgkdyeEHZJLDnculnLkixgkU81GakzF3Mx2nE/7t27fjm2++wd133+2K8RjDihXK4+bMceVI3AMTOpLTr5/YTlAuiYmIEOPIfPj6IcUJFJLTpw/wv/+rLI5Mxeka/ubNmyMoKMgVYzGO995TN07v+vUDgoPlY4KDmdCZlacn4OMjH+PtzfIEs+LrhxTXNJCciAh148gwnE74X3/9dbz88ssoKytzxXiM4fff1Y0jMrLSUsebruXkiHFERFS7u+5S1gb7rrsaZjzkNpxO+JctW4Z//etfCAkJQbdu3XD77bdLvgjAHXeoG6d36enAxYvyMRcvsguLWSUkqBtHxsLXDymuaSA5P/3kuOOd1SrGkak4XcN///33u2AYBjNlCvDNN8rizIA1pyTH0ey+s3FkLHz9kOKaBpJTUKBuHBmG0wn/K6+84opxGEtJibpxesc3KJITFQVkZyuLI/Ph64eUbZF7QYH91pwWi3i/WdY0kNR//6tuHBkGN95yBX9/deP0jovuSM6GDerGkbHYElw5Zuri5OkJrFwp/rl6a0Xb9ytWcJG7WTl6r3U2jgxDUcIfFBSECxcuALjepae2LwK79BA5w88PiI6Wj4mOFuPIfDw9gbFj5WPGjDFXghsfD2zZArRuLT0eHi4ej4/XZlykPa7xoFooKul5++234f/nbPQKpT3mFVi7di3Wrl2LvD93fOvSpQtefvll3HfffVUxR44cwaxZs5CWlobKykp06dIF//jHP3DzzTerNg7Vcac7KWcW3XGjGHPatw9o395+nX5UlHg/mZPVCmzaJB+TmAgsWmS+pD8uTnzdPHNGLGnq189cPwOq6c/JWdXiyDAUJfwTJkyw++f6Cg8Px+LFi9G+fXsAwMcff4y4uDj88ssv6NKlC3JyctC3b19MmjQJr776KgIDA3HkyBE0btxYtTG4BGuSpbjojhxJSgJOnLB/34kT4v2ctTQnRzvLAubdWdbT03zPmeQ5+l1xNo4Mw+lFu2oaOXKk5PsFCxZg7dq12LNnD7p06YJ58+Zh2LBhWLJkSVVMu3btGnqYztuwQVl9vllqkrnojuRYrcC0afYXINpMny7OZnL20nw4YUCknNLqB3eukiCXcJtFu1arFYmJibhy5QpiYmJQWVmJb775Bh06dMC9996LVq1aoU+fPti6davs45SXl6O4uFjy1eBYkyzFRbskx9EMriBcn8El8+GEAZFygwapG0eGoXnCn52dDT8/P/j4+GDKlClITk5G586dcf78eZSWlmLx4sUYOnQovv32WzzwwAOIj49HWlparY+3aNEiBAYGVn1FaLV99OzZ9bvfaMrL63c/GRdncEmOrUtP9Y40NhaLubr0EMnp2VPdODIMiyDIXUd3vYqKCpw8eRKXL1/GF198gffffx9paWlo1qwZWrdujbFjx2Ljxo1V8aNGjULTpk2xqZZFXOXl5Si/IXksLi5GREQEioqKEBAQ4PLnA0AsUQgJkV+oGhwMnDtnjhKF778HYmMdx333HTB4sOvHo7XaEhd7tP31bBipqcDAgY7jdu5kvbJZJSUBo0eLf77xd8L2u8TONESivn2BH390HHf33cDu3a4fD7lccXExAgMDHea59Z7hLy4uxtatW3HkyJE6/X1vb2+0b98evXr1wqJFi9CjRw+sXLkSLVq0gJeXFzp37iyJ79SpE06ePFnr4/n4+CAgIEDy1eBSU5V1pUlNbYjRaE/p8zTLz4OkOINLjrANJZEyMvlRneLIMJxO+B9++GGsXr0aAHD16lX06tULDz/8MLp3744vvvii3gMSBAHl5eXw9vZGdHQ0jh49Krn/2LFjaNOmTb3/HZdigkukHDcSIiXi48VWxjt3Ahs3ire5ueZO9q1W8X1k0ybx1mrVekSkNaVlzFqVO5NmnE74d+3ahX5/zrQlJydDEARcvnwZq1atwhtvvOHUY82dOxfp6enIy8tDdnY25s2bh9TUVIwfPx4AMGPGDGzevBnr1q3Db7/9htWrV2Pbtm2YOnWqs8NuWJcvqxund0rLMFiuYV6cwSUlbG0ox44Vb838ITApCYiMFMvhxo0TbyMjxeNkXjNnqhtHhuF0wl9UVFS1o25KSgoefPBB+Pr6Yvjw4Th+/LhTj3Xu3DkkJCSgY8eOGDx4MPbu3YuUlBQMGTIEAPDAAw/g3XffxZIlS9CtWze8//77+OKLL9C3b19nh92wlLbbNEtbzgEDHHck8vNjwm92nMElUsa2pqF6d6uCAvE4k37zysxUN44Mw+k+/BEREfj5558RFBSElJQUJCYmAgAKCwud3hDrgw8+cBgzceJETJw40dlhaqu0VN04I/DxkX++Pj4NNxZyX9xIiEie3L4VgiCWwXHfCvOqrFQ3jgzD6Rn+6dOnY/z48QgPD0dYWBgG/PnmvGvXLnTr1k3t8emTt7e6cXqXnq5sETP7rBORnIoKcT3Hs8+KtxUVWo+o4XHfCpLTrJm6cWQYTs/wT506FX369MHJkycxZMgQeHiInxnatWuHBQsWqD5AXZoxA3jtNWVxZsA+60RUXzNnAsuXSxemvvgi8PzzwA27sRseX09JDtcQNiyrVfxwfeaMuPlfv35ue2WtTm0577jjDjzwwAPwu6Eu+9Zbb8VLL72k2sB0bd48deP0jjtlElF9zJwJvPVWzS40Vqt43EwLEFu1UjeOiOpGZwvnVdtpt7S0VHYHXFPx9hY33pITEmKekh72WSeiuqqoEGf25Sxfbs7yHqLqWNLTMHS4cF61hJ9uUFoq7qIr59w58yzaZZ91IqqrNWsc95e3WsU4Mzh/Xt04MhaW9Lieo4XzgLhw3s32xWDC7wrjxqkbZwTss05EdaG03bOTbaF1iyWSJMdDYVqnNI5q0unCeacX7ZICv/6qbpxRxMeLreJ0ssCFiNxAbaWAdY3TO1uJZEGB/RlGi0W8nyWS5nT33erGUU06XTivOOHv2bMnLDIvqGVlZaoMyBCaNxc3DFISZzbss05ySkuBhAQgJweIihI3p3O0aRsZW58+wP/+r7I4M7CVSI4eLSb3Nyb9LJGkf/9bedzQoa4di1Hp9Cqb4oT//vvvd+EwDObVV4GRI5XFEZGod28gI+P699nZgL8/EB0N7Nun3bhIW2Fh6sYZga1Ecto0aWlBeLiY7LNE0ry2bVMeN2uWa8diVDq9yqY44X/llVdcOQ5jue8+oFEj4I8/ao9p1EiMI6Kayf6NMjLE+5n0E13HEkkibej0KhtXbbiCpyeQmCgfk5jodv8ZiDRRWlp7sm+TkWGerlYkxa40tbOVSI4dK97yPYWUlumwnKd+dNiIRFHCP3ToUPz0008O40pKSvDmm2/if5XUWxpdfDzwxRdAaKj0+E03icfd8D8DNZAmTdSN07uEBHXjyFh0Wi9LpIm9e9WNo9rFxwN5ecDOncDGjeJtbq7b5neKSnoeeughPPzww/D398eoUaPQq1cvhIWFoXHjxigsLMThw4exe/du/POf/8SIESPw1ltvuXrc+sBLrjXpaBtql7l6Vd04vcvJUTeOjMVWLyvXBo8b9xGJ2CWwYemoEYmihH/SpElISEjAli1bsHnzZqxbtw6X/9y0wWKxoHPnzrj33nuxf/9+dOzY0ZXj1R8d/WdwuaQk+4vMVq5020/E1ACiosQFukriyHw8PcWSFbmJpDFjzDdxQGTPtWvqxpFhWATB3hJjx4qKinD16lUEBwejUaNGao9LNcXFxQgMDERRURECAgK0Ho552bahrv7fzbbAxU1r3lzCmX7hdfv11JfSUrEbjyMlJWzRaUZWKxAZ6XiGPzeXST/R8OHAP//pOG7YMOCbb1w/HnI5pXlunRftBgYGIjQ01K2TfbdgtQKpqcCmTeKtm2213CB0ug01NRA/P7H1ppzoaCb7ZuVoV0vALXe1bBB8f6HqlOZkzN1Mh116XCkpSZyZGjgQGDdOvI2MFI+biU63oaYGtG9f7Uk/+/Cbm053tXQ5vr+QPSEh6saRYTDhdxVbCUv1RLegQDxuphdlvmGTEvv2iWU7998PdOsm3paUMNk3O3bpqYnvL1QbtrGlWjDhdwWWsEjxDZuU8vMDkpOBgwfFW5bxkK1LT21rXywWc3Xp4fsLyaneCry+cWQYTPhdgSUsUv36AcHB8jHBweZ5wyYi5Wy7WgI1k3433tXSZfj+QnIiI9WNI8NQ1JbTnoqKCpw/fx6VlZWS4zfffHO9B6V7LGEhIlKPbVdLe219V6wwT4cvgO8vJE/BJqlOxZFhOJ3wHz9+HBMnTqyx864gCLBYLLDyMiJLWKpLTwcuXpSPuXhRjOOeBURkDzcyFLVqpW4cGcuJE+rGkWE4nfA/9thj8PLywtdff42bbroJFmd6iptFv36AhwdQ7eqHhIeHeUpYOCNFRGrgRoZE8pTuN8R9iUzH6YQ/KysL+/fvx6233uqK8RhDUZF8sg+I9xcVAUFBDTMmLfGKBxGROtiFheTcfbeycp2773b9WMitOL1ot3Pnzrhw4YIrxmIc/furG6d37LJBRKQOlvSQHG68RbVwOuF/8803MXPmTKSmpuLixYsoLi6WfBGA06fVjdM7dtkgIiJyvWbN1I0zkooKMdd49lnxtqJC6xE1KKdLemJjYwEAgwcPlhznot0bhIUBly4pizMLdtkgIqo/lvSQnMuX1Y0zipkzgeXLpftTvPgi8PzzwJIl2o2rATmd8O/cudMV4zCWtDTHfedtcWbCLhtERPXDkp6GZbXq6z3L3oZs9Ykzgpkzgbfeqnncar1+3ARJv9MJf3+z1J3XR1AQEBICnDtXe0xIiDkW7FbHLhtERKQHSUn2r0qvXOm+V6UdtcB2Nk7vKiqAZcvkY5YtA954A/D2bpgxaaROG29dvnwZH3zwAY4cOQKLxYLOnTtj4sSJCAwMVHt8+nX2rLh1tb2kPyREvJ+IiMgZLOlpGElJwOjRNWfCCwrE41u2uGfSLzfRWJc4vXvnHWVdE995B3jhhYYZk0acXrSbmZmJqKgovP3227h06RIuXLiA5cuXIyoqCgcOHHDFGPXr7FnxU3TXruJsfteu4vdM9omIqC7Y5tj1rFZxZt9e2Yvt2PTp0npwd9G4sbpxepeerm6cjjmd8D/33HMYNWoU8vLykJSUhOTkZOTm5mLEiBGYPn26U4+1du1adO/eHQEBAQgICEBMTAy2b99uN/bJJ5+ExWLBihUrnB2ytgIDxU+Oq1eLt7wKQkREdcU2x66Xni4t46lOEID8fPdMEk+eVDdO765cUTdOx+o0wz9r1ix4eV2vBvLy8sLMmTORmZnp1GOFh4dj8eLFyMzMRGZmJgYNGoS4uDgcOnRIErd161bs3bsXYXrravPZZ0DTpsDAgcC4ceJtZKR4qZCIiMhZbHPsenreHb6kRN04vevUSd04HXM64Q8ICMBJO58M8/Pz4e/v79RjjRw5EsOGDUOHDh3QoUMHLFiwAH5+ftizZ09VTEFBAZ555hl89tlnaKSnjSLatwceeQQoL5ceP3VKrP9j0k9E5ByrFUhNBTZtEm/dsaSiIcTHiy0FPaq9hXt4iMfdsbZcT/RcNtWunbpxerdtm7pxOuZ0wv/Xv/4VkyZNwubNm5Gfn49Tp04hMTERTzzxBMaOHVvngVitViQmJuLKlSuIiYkBAFRWViIhIQEzZsxAly5dFD1OeXm59puBRUcDOTm13y8I7lv/R0TkjpKSxCukvGIqPuelS2u+h1it4nEz/kzUpOeyqSefVDdO74qK1I3TMae79CxduhQWiwWPPvoorl27BgBo1KgRnnrqKSxevNjpAWRnZyMmJga///47/Pz8kJycjM6dOwMQd/X18vLC//zP/yh+vEWLFuHVV191ehyqKS0FlJQ22er/2KLSfPz9lV1OdfKKGZFh6bVjiivILSi1mT5d3POEZT11YyubGj1aTO5v/Fm7e9nUf/+rbpzecV+CKk7P8Ht7e2PlypUoLCxEVlYWfvnlF1y6dAlvv/02fHx8nB5Ax44dkZWVhT179uCpp57ChAkTcPjwYezfvx8rV67ERx99BEttn7LtmDNnDoqKiqq+8vPznR5TvThzlcMd6//I9f7+d3XjiIxMzx1TXEHPC0r1xLY7fOvW0uPh4e79AXPrVnXj9E7pxJkJJtjq1IcfAHx9fdGtWzcUFxfj22+/RceOHdGpDosevL290b59ewBAr169kJGRgZUrV6JTp044f/48br755qpYq9WKF154AStWrEBeXp7dx/Px8anTBw/VOLNw2R3r/4iI3IkzCa4ZrpjqeUGp3uhxd3h2pZHivgRVnE74H374YfzlL3/BM888g6tXr6JXr17Iy8uDIAhITEzEgw8+WK8BCYKA8vJyJCQkIDY2VnLfvffei4SEBDz++OP1+jdcSmlvWw8P96z/I9dTuuW90jgiI2OCK6XnBaV6pLfd4aOigB07lMWZgdIKEScqSfTK6ZKeXbt2od+fiWpycjIEQcDly5exatUqvPHGG0491ty5c5Geno68vDxkZ2dj3rx5SE1Nxfjx4xEcHIyuXbtKvho1aoTQ0FB07NjR2WE3nMmTlcWNGePeswRERO6ACa6UbUGpHHddUEquN2KEunF617atunE65nTCX1RUhKCgIABASkoKHnzwQfj6+mL48OE4fvy4U4917tw5JCQkoGPHjhg8eDD27t2LlJQUDBkyxNlhuQ+lWzOvX+/acZD7UrrlvdI4IiNjgivl6el4rRgnlMzrhrbmqsTp3Z13qhunY06X9ERERODnn39GUFAQUlJSkJiYCAAoLCxEYye3av7ggw+ciq+tbt+teHsDM2YAb71Ve8yMGWIcmRNLeuwrLQUSEsSWtlFRwIYNgJ+f1qMirdkSXLnXVDMluFaruA+BnMREYNEi8/xM6LrKSnXj9I47D1dxeoZ/+vTpGD9+PMLDwxEWFoYBf9a27dq1C926dVN7fPq0ZAkQEmL/vpAQ8X4iuq53b7FLwtatQHa2eOvvLx4nc1Oa4LJLz3Xs0mNewcHqxundiRPqxumY0wn/1KlTsWfPHnz44YfYvXs3PP7c6a9du3ZO1/AbVu/eta/4PneOSYzZnT2rbpze9e4NZGTYvy8jg78vZscEV4qLmElObZONdY3Tuz/+UDdOx+rUlvOOO+7AHXfcITk2fPhwVQake6WltScvNhkZYhzLFcyJG6Ncx98XcoQJrhQXMZOc0FB14/SusFDdOB2rU8J/6tQpfPXVVzh58iQqKiok9y1fvlyVgenWuHHK4776yrVjIfeUm6ssrmVL147DHSQkKI9LTnbtWMg9McGVsi1iLiiwvxmZxSLeb5ZFzERyGjcGysqUxRmc0wn/999/j1GjRqFt27Y4evQounbtWtWH//bbb3fFGPVF6cZbzmzQRcZhtQIbNyqLrb7DoxHl5KgbR8bDBFfK0xNYuRIYPVp87jf+TGy9xFes4IJds2LJqJSXwjRXaZyOOV3DP2fOHLzwwgv49ddf0bhxY3zxxRfIz89H//798dBDD7lijPqidJdfLXcD1orVCqSmigvwUlPNs8juRunpwIULjuNatjRHAqN08xezbBJDNdkSXKDm5jhmTXDj44EtW2pOCoSHi8fj47UZF2mPJaNSt96qbpyOOZ3wHzlyBBMmTAAAeHl54erVq/Dz88Nrr72GN998U/UB6s6TT6obZxRJSUBkJDBwoFjONHCg+H1SktYja1hK64zHjzdHArNhg7pxZExMcGuKjwfy8oCdO8Wrhjt3iuWCZvxZ0HVKS0HNUDIKsC3nDZxO+Js2bYry8nIAQFhYGHJuuNR+QcnMpdE9/7y6cUaQlCRefq7eaaOgQDxupqRfaW99s+yC6OcHREfLx0RHc8EuMcElUoL7vEj5+qobp2NOJ/x33nknfvzxRwBiZ54XXngBCxYswMSJE3GnCXYqc8i28ZYcM228ZbUC06bZr721HZs+3ZzlPSTat6/2pD86WryfCBCveg0YIG7ENWCAOa6C1YZXTcmerCx14/SuTRt143TM6VUKy5cvR2lpKQBg/vz5KC0txebNm9G+fXu8/fbbqg9Ql2wba9nbGXLGDHNtvOWoh7YgXO+h/ecmbobGBVX27dvHnXaJlLJdNa0+kWK7amrWMicCfvpJ3Ti9y8tTN07HnE7427VrV/VnX19frFmzRtUBGUZqqnPHjYo9tKW4oKp2fn5svUnkiNUKTJ5c+1VTi0W8ahoXZ+4rIGbFxiFSnGSr4nRJDwBcvnwZ77//PubMmYNLly4BAA4cOICCggJVB6db3Dn0OvbQluKCKiKqjwULgIsXa7//xqumZD5cpCpl74NxfeJ0zOkZ/oMHDyI2NhaBgYHIy8vD//t//w9BQUFITk7Gf/7zH3zyySeuGKd+cOdQKfbQluIuiERUV1br9RaljpjlqilJ5ecriysudu043IXSDbVMsPGW0zP8zz//PB577DEcP34cjW/4Ad13333YtWuXqoPTpbFj1Y3TO/bQJiJSR3o68OdVdYfMctWUrrNalX/Qq/5+bFR//KFunI45nfBnZGTgSTs95Fu3bo2zJqiBcog77dbEHtrXnT+vbhwRmYfSZC442DxXTem69HTlHe/uusu1Y3EXly8ri/uz3byROV3S07hxYxTbuRR09OhRtGTdMdUmPl5cRJaeLr5p3XST+IZktpl99kgmorpSOmv/P/9jvtdWEktnlerUyXXjcBcVFco/AJlgEbPTM/xxcXF47bXX8Meflz8sFgtOnjyJ2bNn48EHH1R9gLrzl7+oG2ck7KFNRHVltYpdzjZtEm/NuHeHbU2UXDlGcDAwb17DjYnch9LubhYLMHWqa8fiDt55R3lsx46uG4ebcDrhX7p0Kf773/+iVatWuHr1Kvr374/27dvD398fCxYscMUY9WXSJHXjyFhY0kPkPG4yJZJbE2Xz3nucTDGr4GBlcUOHmmPzz927lcf+85+uG4ebcLqkJyAgALt378YPP/yAAwcOoLKyErfffjtiY2NdMT796dNH3TgyFrYpJXION5mSsq2JmjZNuqlheLj4YcBMPwuSOndOWdygQa4dh7tQWqbTtCkQGOjasbgBpxN+m0GDBmGQWf7TOOORR5THbdvm2rGQ++nXT2zH+udu1Xb5+3PBHREglu1Mm8ZNpuyp/jMxQR9xcoBNQ6SU7jXQvbtrx+EmFJf07N27F9u3b5cc++STT9C2bVu0atUKkydPRrkJVjk7xF84kmO1AmVl8jFXrpizPpmouvR06Sx2dWbcZMp2xaP6As3Tp8XjZitzousOHVI3Tu+UduhRGqdzihP++fPn4+DBg1XfZ2dnY9KkSYiNjcXs2bOxbds2LFq0yCWD1BVu8kBy1qwBKivlYyorxTgis1PahtIsm0w5uuIBiFc8OGFgTv7+6sbpndI9K5TG6ZzihD8rKwuDBw+u+j4xMRF9+vTBunXr8Pzzz2PVqlX4xz/+4ZJB6goX7ZKcnBx144iMjG1spXjFg+S0aKFunN7xA5CE4hr+wsJChISEVH2flpaGoUOHVn0fHR2NfKVbOhsZZ/hrZ7WyD39kpLpxRGQevOJBctgUQkrpZrAm2TRW8Qx/SEgIcnNzAQAVFRU4cOAAYmJiqu4vKSlBo0aN1B+h3vznP+rGGQXb6om6dVM3jsjI2MZWigkdyWGCK/XnflGqxemc4oR/6NChmD17NtLT0zFnzhz4+vqi3w2dRA4ePIioqCiXDFJXlP4MzPSzsi0yq34p2tZWz0xJ/8WL6sYRGRkTXClHG29ZLEBEBLt8EQGAr6+6cTqnOOF/44034Onpif79+2PdunVYt24dvG/YuOHDDz/EPffc45JB6orS3evMsMsdwEVm1bHGkkg5JrhSchtv2b5fscJ8pZIkatpU3Ti9e+YZdeN0TnHC37JlS6Snp6OwsBCFhYV44IEHJPd//vnneOWVV1QfoO7I9VevS5zecZGZVHa2unFERsYEt6b4eODFFwGPam/fHh7icW68ZV5cpCrVpIm6cTqnOOG3CQwMhKedF9egoCDJjL9p9e+vbpzecZGZ1J/rYFSLIzI6286yrVtLj4eHm2+XXUAsgVy6tOZVUatVPG6mEkmS8lLYh0VpnN598426cTpnkrPegE6fVjdO71iDK9WmjbpxRGYQHy/upmv2Ll9yJZI2Zt15mNgFrjou2pVgwq+2sDBlmziEhbl+LO7AVoNbUGD/TcpiEe83Sw0uEdWNpycwYIDWo9CWMyWSZv9ZEfGKuoTTJT1qWrt2Lbp3746AgAAEBAQgJiYG27dvBwD88ccfmDVrFrp164amTZsiLCwMjz76KE67+8x4Wpq6cXrHGlypEyfUjSMi82CJJMlhgitVUqJunM5pmvCHh4dj8eLFyMzMRGZmJgYNGoS4uDgcOnQIZWVlOHDgAP72t7/hwIEDSEpKwrFjxzBq1Cgth+xYUBBwwwZldoWEiHFmwRrc6/iGXbvSUuCBB4Du3cVbsyxsJ1KKJZIkp6BA3Ti9u3ZN3Tid07SkZ+TIkZLvFyxYgLVr12LPnj2YNGkSduzYIbn/nXfeQe/evXHy5EncfPPNDTlU55w9C4SGAufO1bwvJMQ8m17ciDW4otBQdeOMondvICPj+vfZ2WIniehoYN8+7cZF5E5sJZJyZT1malNKUnJrO+oSp3fNmgEXLiiLMwFNZ/hvZLVakZiYiCtXrkh28L1RUVERLBYLmsmcnPLychQXF0u+NLFmTc1Z/KAg8bhZ2Wpwx44Vb82W7ANAx47qxhlB9WT/RhkZ4v1EJL5mjh0rHzNmjDlfWwn4/Xd14/QuOlrdOJ3TPOHPzs6Gn58ffHx8MGXKFCQnJ6Nz58414n7//XfMnj0b48aNQ0BAQK2Pt2jRIgQGBlZ9RUREuHL49iUlAQ8+WHPx7qVL4nG2TTOvqVMdvxl7eppnY7bS0tqTfZuMDJb3EAFil55Nm+RjEhPNs5EhSd1xh7pxeqe0OYpJmqhonvB37NgRWVlZ2LNnD5566ilMmDABhw8flsT88ccfGDNmDCorK7HGwQz5nDlzUFRUVPWVn5/vyuHXZLWKMyxyxozhC7JZeXsDI0bIx4wYIcaZQUKCunFERuaoSw9gro0MSSo2Vt04vVOa/zV0nqgRzdtyent7o3379gCAXr16ISMjAytXrsTf//53AGKy//DDDyM3Nxc//PCD7Ow+APj4+MDHx8fl467V9u2Oe7r+8YcY5yjxMxqrlTX8Viuwf798zIEDYpwZfjY5OerGERkZF/2TnH79xO53cjX6Fot51ngorfDQohJEA5rP8FcnCALKy8sBXE/2jx8/ju+++w7BwcEaj06BV15RN84okpLEzT4GDgTGjRNvIyPNV97EGTqpqCh144iMjF16SE56uuMFuYJgnveXQ4fUjdM5TRP+uXPnIj09HXl5ecjOzsa8efOQmpqK8ePH49q1axg9ejQyMzPx2WefwWq14uzZszh79iwqKiq0HLa8vDx144wgKQkYPbpmoltQIB43U9LPGTqpDRvUjSMysn79AEcTX8HB5pnBJanUVHXj9M7R5JqzcTqnaUnPuXPnkJCQgDNnziAwMBDdu3dHSkoKhgwZgry8PHz11VcAgNtuu03y93bu3IkB7rqLYECAsp12HZQmGYbcVvCCIF5eNNNW8Erbf5mkTRj8/MQOCXILd6OjxTgiIqpdZaW6cXrnqLza2Tid0zTh/+CDD2q9LzIyEoIee8VOngzMnasszgy4FbzUhx8qj7vvPteOxV3s21d7a0724Se6Lj0duHhRPubiRfO8npKU0rJnPZRHq6FpU3XjdM7tavh174UX1I3TO5awSCktXzJTmRMgJvUlJcD99wPduom3JSVM9oluxNdTktOypbpxenflirpxOseEX23e3sCMGfIxM2aYp+0iF5lJ8ZJr7fz8gORk4OBB8ZZlPERSrVqpG0fGcu6cunF6166dunE6x4TfFZYsEZN6i0V63GIRjy9Zos24tGDbCr76z8LGYjHXVvBKW8Zq2VqWiIj0JytL3Ti9i4tTN07nmPC7ypIl4vbVb78NPPOMePv77+ZK9gFxIe7KleKf7X0AAoAVK8yxYBcA1q1TN46IzOP8eXXjyFiKitSN07tu3dSN0znNN94yNG9vsQON2cXHA1u2iN16blzAGx4uJvvx8ZoNrcFxq28iqiuW9JCc3Fx14/Tu00+Vxw0b5tqxuAEm/NQw4uPFy2Zm32mXM3RERESul52tbpzOMeGnhuPpyVZxnKEjorrihAHJiYpStmusWXYuV7rfkUn2RWINvytZreKOdps2ibdWq9YjIiIivWLXM5Lz5JPqxundyJHqxukcE35XSUoCIiOBgQOBcePE28hI8/VXJ6mzZ9WNIyLzsHU9k2OmrmckVViobpzeVVSoG6dzTPhdISkJGD265g6zBQXicSb95sWEn4jqytMTGDtWPmbMGPOtjSIR31+k1q5VN07nmPCrzWoVu9EIQs37bMemT2d5j1kdPqxuHBGZh9UqlojKSUzk+4tZXbigbpzelZWpG6dzTPjVlp5ec2b/RoIA5OeLcWQ+69erG0dE5uHo/QXg+4uZOfq/4Wyc3jVurG6czjHhV9uZM+rGkbHYu/JTnzgiMo+CAnXjyFgcre9wNk7vOnRQN07nmPCrjV0UiIjIFf77X3XjyFhatlQ3Tu+UtCh1Jk7nmPCrzdZFwWKxf7/Fwi4KRETkPCZ0JKd5c3Xj9I5X1CWY8KvN0xNYuVL8c/Wk3/b9ihXsokBERM5p3VrdODKWr75SN07vWOIkwYTfFeLjgS1bgLAw6fHWrcXj8fHajIuIiPSLffhJzuXL6sbpXVqaunE6x4TfVfbsqbkw9/Rp8TgREZGz2Ief5Pz73+rG6Z2fn7pxOseE3xVmzgTeeguorJQer6wUj8+cqc24iIhIv9iHn+RUzznqG6d3b72lbpzOWQTB2KsViouLERgYiKKiIgQEBLj+H6yoAHx95V9wPT3FjR68vV0/HnIvtS3mtsfYv5pE5KzUVGDgQMdxO3cCAwa4ejTkbkJDgXPnHMeFhJhjt93mzZWVLzVrBhQWuno0LqM0z+UMv9rWrHE8u2K1inFERERKcZ8XkjNkiLpxeldUpG6czjHhV1tOjrpxREREAPd5IXm7dqkbp3dK17KYZM0LE361RUWpG0dkFhUVYsvaZ58VbysqtB4RkXvhPi8kR2mZjhnKeQCxxEnNOJ1jwq+2qVMdf1r09BTjiEg0c6a49uW554DVq8VbX18ucCe6Efd5ITkeClM6pXF616GDunE6Z5Kz3oC8vYHnn5ePef55LtglsrF1taq+9sVqZVcrouri44EXX6yZtHl4iMe5z4t5cSdmqXbt1I3TOSb8rrBkCTBjRs1ZFk9P8fiSJdqMi8jdVFQAy5fLxyxfzvIeIpukJGDpUvsfkJcuFe8nc7rlFnXj9I6L3CWY8LvKkiVi68233waeeUa8LStjsk90I3a1IlLOagWmTZNv2Tt9Ovvwm1WnTurG6d2hQ+rG6RwTflfy9ARuuw246y7xlnWVRFJHj6obR2Rk6enAqVO13y8IQH6+GEfm89pr6sbpXWmpunE6x4TfVZKSgMhIcZOUcePE28hIXm4luhG7ShApxxIFkjN3rrpxesc2thJM+F0hKQkYPbrmTExBgXicST+RiC/IRMq1aqVuHBnLzp3qxumdkl12nYnTOSb8apOrsbQdY40lkYht04hIDRZLzS+zUdr9zyxdArnTroSmCf/atWvRvXt3BAQEICAgADExMdi+fXvV/YIgYP78+QgLC0OTJk0wYMAAHHL3xRWssSRSjvtWECl3/ry6cUYhtxGZmTzyiLpxenflirpxOqdpwh8eHo7FixcjMzMTmZmZGDRoEOLi4qqS+iVLlmD58uVYvXo1MjIyEBoaiiFDhqCkpETLYctjjSWRcty3gkg5lsDV5CipN1PS/9xz6sbpXaNG6sbpnKYJ/8iRIzFs2DB06NABHTp0wIIFC+Dn54c9e/ZAEASsWLEC8+bNQ3x8PLp27YqPP/4YZWVl2Lhxo5bDlscXZCLncN8KImX69QPCw+VntCMixDgzUJrMmyXp9/YWXzPlzJhhngmUHj3UjdM5t6nht1qtSExMxJUrVxATE4Pc3FycPXsW99xzT1WMj48P+vfvj59++qnWxykvL0dxcbHkq0HxBZnIedy3gsgxT09g5Urxz9XfY2zfr1jBFtBmZptAqc5iMd8Eyrx56sbpnOYJf3Z2Nvz8/ODj44MpU6YgOTkZnTt3xtk/2/CFhIRI4kNCQqrus2fRokUIDAys+oqIiHDp+GvgCzJR3Xh7iwva33lHvDXLLBSRM+LjgS1bgNatpcfDw8Xj8fHajIvcx5IlQHm5dALl99/NlewD7MNfjZfWA+jYsSOysrJw+fJlfPHFF5gwYQLS0tKq7rdUS5oFQahx7EZz5szB8zfUBBcXFzd80m97QZ42TbqANzxcTPb5gmxezZopawHWrJmLB0JEuhUfD8TFic0fzpwRS0T79eNEEl1nm0AxM5ZYS1gEQW6P7oYXGxuLqKgozJo1C1FRUThw4AB69uxZdX9cXByaNWuGjz/+WNHjFRcXIzAwEEVFRQgICHDVsO2zWvmCTFJnzyp7cTlzBggNdf14iIj0ypnafPdKdaghWK3ihqcFBfbPv8UiTsTm5uo6N1Oa52pe0lOdIAgoLy9H27ZtERoaih07dlTdV1FRgbS0NNx1110ajtAJnp7AgAHA2LHirY7/Q5FKQkMBX1/5GF9fJvtERI4oTeKZ7JsTS6wlNC3pmTt3Lu677z5ERESgpKQEiYmJSE1NRUpKCiwWC6ZPn46FCxfilltuwS233IKFCxfC19cX48aN03LYRPVz5QrQtKm4KLU6X1/T9AQmIqo3QZCf6Weyb24ssa6iacJ/7tw5JCQk4MyZMwgMDET37t2RkpKCIUOGAABmzpyJq1evYurUqSgsLESfPn3w7bffwt/fX8thE9XflStiec9tt4k1/c2aAVlZnNknInJWbUk/k30CuOblT25Xw682TWv4iYiIiIhcRLc1/EREREREpB4m/EREREREBsaEn4iIiIjIwDTfeMvVbEsUiouLNR4JEREREZF6bPmtoyW5hk/4S0pKAKDhd9slIiIiImoAJSUlCAwMrPV+w3fpqaysxOnTp+Hv7w+LM7vyVVNcXIyIiAjk5+ez24/GeC7cC8+H++C5cB88F+6D58J98FyoTxAElJSUICwsDB4etVfqG36G38PDA+Hh4ao9XkBAAP+TugmeC/fC8+E+eC7cB8+F++C5cB88F+qSm9m34aJdIiIiIiIDY8JPRERERGRgTPgV8vHxwSuvvAIfHx+th2J6PBfuhefDffBcuA+eC/fBc+E+eC60Y/hFu0REREREZsYZfiIiIiIiA2PCT0RERERkYEz4iYiIiIgMjAk/EREREZGBmSLhX7RoEaKjo+Hv749WrVrh/vvvx9GjRxX//R9//BFeXl647bbbJMc/+ugjWCyWGl+///57VczatWvRvXv3qk0mYmJisH37drWemu5oeS6qj8NisWD69On1eDb6puW5mD9/fo37Q0ND1XpquqP170VBQQEeeeQRBAcHw9fXF7fddhv279+vxlPTHS3PRWRkpN2Yp59+Wq2npytanotr167hpZdeQtu2bdGkSRO0a9cOr732GiorK9V6erqi5bkoKSnB9OnT0aZNGzRp0gR33XUXMjIy1HpqpmGKhD8tLQ1PP/009uzZgx07duDatWu45557cOXKFYd/t6ioCI8++igGDx5s9/6AgACcOXNG8tW4ceOq+8PDw7F48WJkZmYiMzMTgwYNQlxcHA4dOqTa89MTLc+FTUZGBt577z1079693s9Hz7Q+F126dJHcn52drcrz0iMtz0VhYSHuvvtuNGrUCNu3b8fhw4exbNkyNGvWTK2npytanouMjAzJfTt27AAAPPTQQ+o8OZ3R8ly8+eabePfdd7F69WocOXIES5YswVtvvYV33nlHteenJ1qeiyeeeAI7duzAhg0bkJ2djXvuuQexsbEoKChQ7fmZgmBC58+fFwAIaWlpDmP/+te/Ci+99JLwyiuvCD169JDct379eiEwMNDpf7958+bC+++/7/TfM6KGPhclJSXCLbfcIuzYsUPo37+/MG3atLoN3IAa8lzY+3t0XUOei1mzZgl9+/atx2iNTcv3i2nTpglRUVFCZWWlU3/PqBryXAwfPlyYOHGi5Fh8fLzwyCOPODtsQ2qoc1FWViZ4enoKX3/9teR4jx49hHnz5tVl6KZlihn+6oqKigAAQUFBsnHr169HTk4OXnnllVpjSktL0aZNG4SHh2PEiBH45Zdfao21Wq1ITEzElStXEBMTU7fBG0xDn4unn34aw4cPR2xsbP0GbkANfS6OHz+OsLAwtG3bFmPGjMGJEyfq9wQMpCHPxVdffYVevXrhoYceQqtWrdCzZ0+sW7eu/k/CILR6v6ioqMCnn36KiRMnwmKx1G3wBtOQ56Jv3774/vvvcezYMQDAv//9b+zevRvDhg2r57MwhoY6F9euXYPVaq1xhbhJkybYvXt3PZ6BCWn9iaOhVVZWCiNHjnQ4o3Xs2DGhVatWwtGjRwVBsD8j+fPPPwsbNmwQsrKyhF27dgkPPvig0KRJE+HYsWOSuIMHDwpNmzYVPD09hcDAQOGbb75R9TnpVUOfi02bNgldu3YVrl69KgiCwBn+GzT0ufjnP/8pbNmyRTh48GDV1ZaQkBDhwoULqj83vWnoc+Hj4yP4+PgIc+bMEQ4cOCC8++67QuPGjYWPP/5Y9eemN1q8X9hs3rxZ8PT0FAoKClR5LnrX0OeisrJSmD17tmCxWAQvLy/BYrEICxcuVP156VFDn4uYmBihf//+QkFBgXDt2jVhw4YNgsViETp06KD6czMy0yX8U6dOFdq0aSPk5+fXGnPt2jWhV69ewtq1a6uOKSlBsFqtQo8ePYRnn31Wcry8vFw4fvy4kJGRIcyePVto0aKFcOjQoXo9DyNoyHNx8uRJoVWrVkJWVlZVDBP+67T4vbhRaWmpEBISIixbtszpsRtNQ5+LRo0aCTExMZK4Z599Vrjzzjvr9gQMRMvfi3vuuUcYMWJEncZtRA19LjZt2iSEh4cLmzZtEg4ePCh88sknQlBQkPDRRx/V+7noXUOfi99++034y1/+IgAQPD09hejoaGH8+PFCp06d6v1czMRUCf8zzzwjhIeHCydOnJCNKywsrPqPZfuyWCxVx77//vta/+4TTzwhDB06VPbxBw8eLEyePLlOz8EoGvpcJCcn13gcAILFYhE8PT2Fa9euqfr89MRdfi9iY2OFKVOm1Ok5GIUW5+Lmm28WJk2aJIlZs2aNEBYWVr8no3Na/l7k5eUJHh4ewtatW+v9PIxAi3MRHh4urF69WhLz+uuvCx07dqzfk9E5LX8vSktLhdOnTwuCIAgPP/ywMGzYsPo9GZPxavAaIg0IgoBnn30WycnJSE1NRdu2bWXjAwICanQMWbNmDX744Qds2bKl1r8vCAKysrLQrVs3h+MpLy937kkYhFbnYvDgwTUe5/HHH8ett96KWbNmwdPTsx7PSp/c6feivLwcR44cQb9+/Zx/Igag5bm4++67a7TXO3bsGNq0aVPHZ6Nv7vB7sX79erRq1QrDhw+v+xMxAC3PRVlZGTw8pMscPT09TduW0x1+L5o2bYqmTZuisLAQ//rXv7BkyZK6PyEz0uBDRoN76qmnhMDAQCE1NVU4c+ZM1VdZWVlVzOzZs4WEhIRaH8Pepaj58+cLKSkpQk5OjvDLL78Ijz/+uODl5SXs3bu3KmbOnDnCrl27hNzcXOHgwYPC3LlzBQ8PD+Hbb79V/XnqgZbnojqzl/RoeS5eeOEFITU1VThx4oSwZ88eYcSIEYK/v7+Ql5en+vPUAy3Pxb59+wQvLy9hwYIFwvHjx4XPPvtM8PX1FT799FPVn6ceaP0aZbVahZtvvlmYNWuWqs9Lj7Q8FxMmTBBat24tfP3110Jubq6QlJQktGjRQpg5c6bqz1MPtDwXKSkpwvbt24UTJ04I3377rdCjRw+hd+/eQkVFherP08hMMcO/du1aAMCAAQMkx9evX4/HHnsMAHDmzBmcPHnSqce9fPkyJk+ejLNnzyIwMBA9e/bErl270Lt376qYc+fOISEhAWfOnEFgYCC6d++OlJQUDBkypF7PSa+0PBckpeW5OHXqFMaOHYsLFy6gZcuWuPPOO7Fnzx7TzipreS6io6ORnJyMOXPm4LXXXkPbtm2xYsUKjB8/vl7PSa+0fo367rvvcPLkSUycOLHOz8EotDwX77zzDv72t79h6tSpOH/+PMLCwvDkk0/i5Zdfrtdz0istz0VRURHmzJmDU6dOISgoCA8++CAWLFiARo0a1es5mY1FEARB60EQEREREZFrmLIPPxERERGRWTDhJyIiIiIyMCb8REREREQGxoSfiIiIiMjAmPATERERERkYE34iIiIiIgNjwk9EREREZGBM+ImITGj+/Pm47bbbNPv3//a3v2Hy5Mn1fpzy8nLcfPPN2L9/vwqjIiJS165duzBy5EiEhYXBYrFg69atTj+GIAhYunQpOnToAB8fH0RERGDhwoVOPQYTfiIig7FYLLJfjz32GF588UV8//33mozv3LlzWLlyJebOnVt17LHHHsP9998v+d423kaNGiEkJARDhgzBhx9+iMrKyqo4Hx8fvPjii5g1a1ZDPgUiIkWuXLmCHj16YPXq1XV+jGnTpuH999/H0qVL8X//93/Ytm1bjV26HfGq879ORERu6cyZM1V/3rx5M15++WUcPXq06liTJk3g5+cHPz8/LYaHDz74ADExMYiMjJSNGzp0KNavXw+r1Ypz584hJSUF06ZNw5YtW/DVV1/By0t8Cxs/fjxmzJiBI0eOoFOnTg3wDIiIlLnvvvtw33331Xp/RUUFXnrpJXz22We4fPkyunbtijfffBMDBgwAABw5cgRr167Fr7/+io4dO9Z5HJzhJyIymNDQ0KqvwMBAWCyWGseql/TYZtgXLlyIkJAQNGvWDK+++iquXbuGGTNmICgoCOHh4fjwww8l/1ZBQQH++te/onnz5ggODkZcXBzy8vJkx5eYmIhRo0Y5fB4+Pj4IDQ1F69atcfvtt2Pu3Ln48ssvsX37dnz00UdVccHBwbjrrruwadMmZ35MRESae/zxx/Hjjz8iMTERBw8exEMPPYShQ4fi+PHjAIBt27ahXbt2+Prrr9G2bVtERkbiiSeewKVLl5z6d5jwExERAOCHH37A6dOnsWvXLixfvhzz58/HiBEj0Lx5c+zduxdTpkzBlClTkJ+fDwAoKyvDwIED4efnh127dmH37t3w8/PD0KFDUVFRYfffKCwsxK+//opevXrVaYyDBg1Cjx49kJSUJDneu3dvpKen1+kxiYi0kJOTg02bNuHzzz9Hv379EBUVhRdffBF9+/bF+vXrAQAnTpzAf/7zH3z++ef45JNP8NFHH2H//v0YPXq0U/8WS3qIiAgAEBQUhFWrVsHDwwMdO3bEkiVLUFZWVlVrP2fOHCxevBg//vgjxowZg8TERHh4eOD999+HxWIBAKxfvx7NmjVDamoq7rnnnhr/xn/+8x8IgoCwsLA6j/PWW2/FwYMHJcdat27t8MoCEZE7OXDgAARBQIcOHSTHy8vLERwcDACorKxEeXk5Pvnkk6q4Dz74AHfccQeOHj2quMyHCT8REQEAunTpAg+P6xd+Q0JC0LVr16rvPT09ERwcjPPnzwMA9u/fj99++w3+/v6Sx/n999+Rk5Nj99+4evUqAKBx48Z1HqcgCFUfMGyaNGmCsrKyOj8mEVFDq6yshKenJ/bv3w9PT0/JfbY1VjfddBO8vLwkHwpsa5VOnjzJhJ+IiJzTqFEjyfe2DjnVj9m65FRWVuKOO+7AZ599VuOxWrZsafffaNGiBQCxtKe2GEeOHDmCtm3bSo5dunSpzo9HRKSFnj17wmq14vz58+jXr5/dmLvvvhvXrl1DTk4OoqKiAADHjh0DALRp00bxv8WEn4iI6uT222/H5s2b0apVKwQEBCj6O1FRUQgICMDhw4drXMZW4ocffkB2djaee+45yfFff/0VPXv2dPrxiIhcqbS0FL/99lvV97m5ucjKykJQUBA6dOiA8ePH49FHH8WyZcvQs2dPXLhwAT/88AO6deuGYcOGITY2FrfffjsmTpyIFStWoLKyEk8//TSGDBni1GsoF+0SEVGdjB8/Hi1atEBcXBzS09ORm5uLtLQ0TJs2DadOnbL7dzw8PBAbG4vdu3c7fPzy8nKcPXsWBQUFOHDgABYuXIi4uDiMGDECjz76qCQ2PT3d7poBIiItZWZmomfPnlUTEs8//zx69uyJl19+GYC47unRRx/FCy+8gI4dO2LUqFHYu3cvIiIiAIivmdu2bUOLFi3wl7/8BcOHD0enTp2QmJjo1Dg4w09ERHXi6+uLXbt2YdasWYiPj0dJSQlat26NwYMHy874T548GZMmTcKSJUuq1gxUVlZW9dW3SUlJqapfbd68OXr06IFVq1ZhwoQJkrUGP//8M4qKipzuWkFE5GoDBgyAIAi13t+oUSO8+uqrePXVV2uNCQsLwxdffFGvcVgEuVEQERGpTBAE3HnnnZg+fTrGjh0LQNxkq3379nXajfKhhx5Cz549JTv3EhHRdSzpISKiBmWxWPDee+/h2rVrKCwsxDfffIPU1FTExsY6/Vjl5eXo0aNHjZp+IiK6jjP8RESkmQceeAAZGRmYMGEC3njjjRrtNomIqP6Y8BMRERERGRhLeoiIiIiIDIwJPxERERGRgTHhJyIiIiIyMCb8REREREQGxoSfiIiIiMjAmPATERERERkYE34iIiIiIgNjwk9EREREZGBM+ImIiIiIDOz/A1xMw+ZTJ8BuAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 900x300 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=(9, 3))\n",
    "plt.xlabel('Time (JD)')\n",
    "plt.ylabel('Season (Ls in °)')\n",
    "p=plt.plot(df['time_min'], df['solar_longitude_max'], 'ro')\n",
    "p3= plt.show()\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}