{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Get and analyse data from AURN\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get and merge data by a bounds of longitude and latitude\n",
    "\n",
    "Users can give a bounds of longitude and latitude to download the station data."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Load the ObsAQ package and Define the range"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import obsaq\n",
    "\n",
    "lon_min = -9\n",
    "lon_max = 1.8\n",
    "lat_min = 49\n",
    "lat_max = 61\n",
    "bounds=[lon_min, lon_max, lat_min, lat_max]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Get the station metadata\n",
    "\n",
    "NOTE: This is the all station information and the next step will choose the stations by bounds."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "meta = obsaq.meta()\n",
    "site_table = meta.get_metadata('aurn')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Choose the stations by range"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Site is selected by bounds: [-9, 1.8, 49, 61]\n"
     ]
    },
    {
     "data": {
      "application/vnd.microsoft.datawrangler.viewer.v0+json": {
       "columns": [
        {
         "name": "index",
         "rawType": "int64",
         "type": "integer"
        },
        {
         "name": "site_id",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "site_name",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "location_type",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "latitude",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "longitude",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "parameter",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "Parameter_name",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "start_date",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "end_date",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "ratified_to",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "zone",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "agglomeration",
         "rawType": "object",
         "type": "unknown"
        },
        {
         "name": "local_authority",
         "rawType": "object",
         "type": "string"
        }
       ],
       "ref": "700ecc9a-102b-4301-8f9b-3de9d1c44355",
       "rows": [
        [
         "0",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "O3",
         "Ozone",
         "2003-08-01",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "1",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "NO",
         "Nitric oxide",
         "1999-09-18",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "2",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "NO2",
         "Nitrogen dioxide",
         "1999-09-18",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "3",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "NOXasNO2",
         "Nitrogen oxides as nitrogen dioxide",
         "1999-09-18",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "4",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "SO2",
         "Sulphur dioxide",
         "2001-01-01",
         "2007-09-30",
         "2007-09-30",
         "North East Scotland",
         null,
         "Aberdeen City"
        ]
       ],
       "shape": {
        "columns": 13,
        "rows": 5
       }
      },
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>site_id</th>\n",
       "      <th>site_name</th>\n",
       "      <th>location_type</th>\n",
       "      <th>latitude</th>\n",
       "      <th>longitude</th>\n",
       "      <th>parameter</th>\n",
       "      <th>Parameter_name</th>\n",
       "      <th>start_date</th>\n",
       "      <th>end_date</th>\n",
       "      <th>ratified_to</th>\n",
       "      <th>zone</th>\n",
       "      <th>agglomeration</th>\n",
       "      <th>local_authority</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>O3</td>\n",
       "      <td>Ozone</td>\n",
       "      <td>2003-08-01</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>NO</td>\n",
       "      <td>Nitric oxide</td>\n",
       "      <td>1999-09-18</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>NO2</td>\n",
       "      <td>Nitrogen dioxide</td>\n",
       "      <td>1999-09-18</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>NOXasNO2</td>\n",
       "      <td>Nitrogen oxides as nitrogen dioxide</td>\n",
       "      <td>1999-09-18</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>SO2</td>\n",
       "      <td>Sulphur dioxide</td>\n",
       "      <td>2001-01-01</td>\n",
       "      <td>2007-09-30</td>\n",
       "      <td>2007-09-30</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  site_id site_name     location_type  latitude  longitude parameter  \\\n",
       "0     ABD  Aberdeen  Urban Background  57.15736  -2.094278        O3   \n",
       "1     ABD  Aberdeen  Urban Background  57.15736  -2.094278        NO   \n",
       "2     ABD  Aberdeen  Urban Background  57.15736  -2.094278       NO2   \n",
       "3     ABD  Aberdeen  Urban Background  57.15736  -2.094278  NOXasNO2   \n",
       "4     ABD  Aberdeen  Urban Background  57.15736  -2.094278       SO2   \n",
       "\n",
       "                        Parameter_name  start_date    end_date ratified_to  \\\n",
       "0                                Ozone  2003-08-01  2021-09-20  2021-09-20   \n",
       "1                         Nitric oxide  1999-09-18  2021-09-20  2021-09-20   \n",
       "2                     Nitrogen dioxide  1999-09-18  2021-09-20  2021-09-20   \n",
       "3  Nitrogen oxides as nitrogen dioxide  1999-09-18  2021-09-20  2021-09-20   \n",
       "4                      Sulphur dioxide  2001-01-01  2007-09-30  2007-09-30   \n",
       "\n",
       "                  zone agglomeration local_authority  \n",
       "0  North East Scotland           NaN   Aberdeen City  \n",
       "1  North East Scotland           NaN   Aberdeen City  \n",
       "2  North East Scotland           NaN   Aberdeen City  \n",
       "3  North East Scotland           NaN   Aberdeen City  \n",
       "4  North East Scotland           NaN   Aberdeen City  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "final_sites = meta.get_site(bounds=bounds)\n",
    "final_sites.head(5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Start to download the merged targeted station data\n",
    "Download the final merged file for all targeted station data using memory.\n",
    "\n",
    "- `pollutant`: See names of pollutants. Define one pollutant, diverse pollutants or all pollutants to download the data for them.\n",
    "               Pollutants can be defined as \"PM2.5\",\"PM10\",\"O3\",\"NO\",\"NO2\",\"NOXasNO2\" or \"SO2\".\n",
    "- `start`: the start date of data to be downloaded.\n",
    "- `end`: the end date of data to be downloaded.\n",
    "- `year`: the year of data to be downloaded. Defaults to 2010.\n",
    "- `output_dir`: the directory to save the downloaded data. \n",
    "- `download_mode`: \"Stream\" for saving final and intermediate files while \"memory\" for only the final file.\n",
    "- `save_per_site`: whether save files for every station individually.\n",
    "- `save_merged`: whether save the merged file for all selected data.\n",
    "- `add_site_id`: whether include site id in the downloaded file.\n",
    "\n",
    "Warnings are normal for processing data with different format. Ingoring them can be fine. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "     site_id        site_name     location_type   latitude  longitude  \\\n",
      "0        ABD         Aberdeen  Urban Background  57.157360  -2.094278   \n",
      "1        ABD         Aberdeen  Urban Background  57.157360  -2.094278   \n",
      "2        ABD         Aberdeen  Urban Background  57.157360  -2.094278   \n",
      "3        ABD         Aberdeen  Urban Background  57.157360  -2.094278   \n",
      "6        ABD         Aberdeen  Urban Background  57.157360  -2.094278   \n",
      "...      ...              ...               ...        ...        ...   \n",
      "3057    YK11  York Fishergate     Urban Traffic  53.951889  -1.075861   \n",
      "3058    YK11  York Fishergate     Urban Traffic  53.951889  -1.075861   \n",
      "3059    YK11  York Fishergate     Urban Traffic  53.951889  -1.075861   \n",
      "3060    YK11  York Fishergate     Urban Traffic  53.951889  -1.075861   \n",
      "3061    YK11  York Fishergate     Urban Traffic  53.951889  -1.075861   \n",
      "\n",
      "     parameter                             Parameter_name  start_date  \\\n",
      "0           O3                                      Ozone  2003-08-01   \n",
      "1           NO                               Nitric oxide  1999-09-18   \n",
      "2          NO2                           Nitrogen dioxide  1999-09-18   \n",
      "3     NOXasNO2        Nitrogen oxides as nitrogen dioxide  1999-09-18   \n",
      "6         PM10  PM10 particulate matter (Hourly measured)  1999-09-18   \n",
      "...        ...                                        ...         ...   \n",
      "3057     NV2.5       Non-volatile PM2.5 (Hourly measured)  2011-08-31   \n",
      "3058      V2.5           Volatile PM2.5 (Hourly measured)  2011-08-31   \n",
      "3059        wd                    Modelled Wind Direction  2010-08-01   \n",
      "3060        ws                        Modelled Wind Speed  2010-08-01   \n",
      "3061      temp                       Modelled Temperature  2010-08-01   \n",
      "\n",
      "        end_date ratified_to                    zone agglomeration  \\\n",
      "0     2021-09-20  2021-09-20     North East Scotland           NaN   \n",
      "1     2021-09-20  2021-09-20     North East Scotland           NaN   \n",
      "2     2021-09-20  2021-09-20     North East Scotland           NaN   \n",
      "3     2021-09-20  2021-09-20     North East Scotland           NaN   \n",
      "6     2021-09-20  2021-09-20     North East Scotland           NaN   \n",
      "...          ...         ...                     ...           ...   \n",
      "3057  2018-10-10  2018-10-10  Yorkshire & Humberside           NaN   \n",
      "3058  2018-10-10  2018-10-10  Yorkshire & Humberside           NaN   \n",
      "3059     ongoing       Never  Yorkshire & Humberside           NaN   \n",
      "3060     ongoing       Never  Yorkshire & Humberside           NaN   \n",
      "3061     ongoing       Never  Yorkshire & Humberside           NaN   \n",
      "\n",
      "     local_authority  \n",
      "0      Aberdeen City  \n",
      "1      Aberdeen City  \n",
      "2      Aberdeen City  \n",
      "3      Aberdeen City  \n",
      "6      Aberdeen City  \n",
      "...              ...  \n",
      "3057            York  \n",
      "3058            York  \n",
      "3059            York  \n",
      "3060            York  \n",
      "3061            York  \n",
      "\n",
      "[1625 rows x 13 columns]\n",
      "Saved merged file: data/test_pm25_final/MERGED_aurn_PM2.5_20171201_20181130.csv\n"
     ]
    }
   ],
   "source": [
    "merged_df = meta.download_sites(\n",
    "    port=\"aurn\",\n",
    "    pollutant=\"PM2.5\",\n",
    "    start=\"2017-12-01\",\n",
    "    end=\"2018-11-30\",\n",
    "    output_dir=\"data/test_pm25_final\",\n",
    "    download_mode=\"memory\",   \n",
    "    save_per_site=False,    \n",
    "    save_merged=True,\n",
    "    add_site_id=True\n",
    ")\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Merge data by site_id"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Load the ObsAQ package and check the site information"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.microsoft.datawrangler.viewer.v0+json": {
       "columns": [
        {
         "name": "index",
         "rawType": "int64",
         "type": "integer"
        },
        {
         "name": "site_id",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "site_name",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "location_type",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "latitude",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "longitude",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "parameter",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "Parameter_name",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "start_date",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "end_date",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "ratified_to",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "zone",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "agglomeration",
         "rawType": "object",
         "type": "unknown"
        },
        {
         "name": "local_authority",
         "rawType": "object",
         "type": "string"
        }
       ],
       "ref": "fbf595a4-859d-4ecc-9a6a-b305c1ae7b06",
       "rows": [
        [
         "0",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "O3",
         "Ozone",
         "2003-08-01",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "1",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "NO",
         "Nitric oxide",
         "1999-09-18",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "2",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "NO2",
         "Nitrogen dioxide",
         "1999-09-18",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "3",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "NOXasNO2",
         "Nitrogen oxides as nitrogen dioxide",
         "1999-09-18",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ],
        [
         "4",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "SO2",
         "Sulphur dioxide",
         "2001-01-01",
         "2007-09-30",
         "2007-09-30",
         "North East Scotland",
         null,
         "Aberdeen City"
        ]
       ],
       "shape": {
        "columns": 13,
        "rows": 5
       }
      },
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>site_id</th>\n",
       "      <th>site_name</th>\n",
       "      <th>location_type</th>\n",
       "      <th>latitude</th>\n",
       "      <th>longitude</th>\n",
       "      <th>parameter</th>\n",
       "      <th>Parameter_name</th>\n",
       "      <th>start_date</th>\n",
       "      <th>end_date</th>\n",
       "      <th>ratified_to</th>\n",
       "      <th>zone</th>\n",
       "      <th>agglomeration</th>\n",
       "      <th>local_authority</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>O3</td>\n",
       "      <td>Ozone</td>\n",
       "      <td>2003-08-01</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>NO</td>\n",
       "      <td>Nitric oxide</td>\n",
       "      <td>1999-09-18</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>NO2</td>\n",
       "      <td>Nitrogen dioxide</td>\n",
       "      <td>1999-09-18</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>NOXasNO2</td>\n",
       "      <td>Nitrogen oxides as nitrogen dioxide</td>\n",
       "      <td>1999-09-18</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>SO2</td>\n",
       "      <td>Sulphur dioxide</td>\n",
       "      <td>2001-01-01</td>\n",
       "      <td>2007-09-30</td>\n",
       "      <td>2007-09-30</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  site_id site_name     location_type  latitude  longitude parameter  \\\n",
       "0     ABD  Aberdeen  Urban Background  57.15736  -2.094278        O3   \n",
       "1     ABD  Aberdeen  Urban Background  57.15736  -2.094278        NO   \n",
       "2     ABD  Aberdeen  Urban Background  57.15736  -2.094278       NO2   \n",
       "3     ABD  Aberdeen  Urban Background  57.15736  -2.094278  NOXasNO2   \n",
       "4     ABD  Aberdeen  Urban Background  57.15736  -2.094278       SO2   \n",
       "\n",
       "                        Parameter_name  start_date    end_date ratified_to  \\\n",
       "0                                Ozone  2003-08-01  2021-09-20  2021-09-20   \n",
       "1                         Nitric oxide  1999-09-18  2021-09-20  2021-09-20   \n",
       "2                     Nitrogen dioxide  1999-09-18  2021-09-20  2021-09-20   \n",
       "3  Nitrogen oxides as nitrogen dioxide  1999-09-18  2021-09-20  2021-09-20   \n",
       "4                      Sulphur dioxide  2001-01-01  2007-09-30  2007-09-30   \n",
       "\n",
       "                  zone agglomeration local_authority  \n",
       "0  North East Scotland           NaN   Aberdeen City  \n",
       "1  North East Scotland           NaN   Aberdeen City  \n",
       "2  North East Scotland           NaN   Aberdeen City  \n",
       "3  North East Scotland           NaN   Aberdeen City  \n",
       "4  North East Scotland           NaN   Aberdeen City  "
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import obsaq\n",
    "\n",
    "meta = obsaq.meta()\n",
    "site_table = meta.get_metadata('aurn')\n",
    "\n",
    "site_table.head(5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Choose the stations by site_id"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Site is selected by site_id: ABD\n"
     ]
    },
    {
     "data": {
      "application/vnd.microsoft.datawrangler.viewer.v0+json": {
       "columns": [
        {
         "name": "index",
         "rawType": "int64",
         "type": "integer"
        },
        {
         "name": "site_id",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "site_name",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "location_type",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "latitude",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "longitude",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "parameter",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "Parameter_name",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "start_date",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "end_date",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "ratified_to",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "zone",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "agglomeration",
         "rawType": "object",
         "type": "unknown"
        },
        {
         "name": "local_authority",
         "rawType": "object",
         "type": "string"
        }
       ],
       "ref": "d96b94db-dcc2-4375-84e6-4b988ec79665",
       "rows": [
        [
         "0",
         "ABD",
         "Aberdeen",
         "Urban Background",
         "57.15736",
         "-2.094278",
         "O3",
         "Ozone",
         "2003-08-01",
         "2021-09-20",
         "2021-09-20",
         "North East Scotland",
         null,
         "Aberdeen City"
        ]
       ],
       "shape": {
        "columns": 13,
        "rows": 1
       }
      },
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>site_id</th>\n",
       "      <th>site_name</th>\n",
       "      <th>location_type</th>\n",
       "      <th>latitude</th>\n",
       "      <th>longitude</th>\n",
       "      <th>parameter</th>\n",
       "      <th>Parameter_name</th>\n",
       "      <th>start_date</th>\n",
       "      <th>end_date</th>\n",
       "      <th>ratified_to</th>\n",
       "      <th>zone</th>\n",
       "      <th>agglomeration</th>\n",
       "      <th>local_authority</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ABD</td>\n",
       "      <td>Aberdeen</td>\n",
       "      <td>Urban Background</td>\n",
       "      <td>57.15736</td>\n",
       "      <td>-2.094278</td>\n",
       "      <td>O3</td>\n",
       "      <td>Ozone</td>\n",
       "      <td>2003-08-01</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>2021-09-20</td>\n",
       "      <td>North East Scotland</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Aberdeen City</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  site_id site_name     location_type  latitude  longitude parameter  \\\n",
       "0     ABD  Aberdeen  Urban Background  57.15736  -2.094278        O3   \n",
       "\n",
       "  Parameter_name  start_date    end_date ratified_to                 zone  \\\n",
       "0          Ozone  2003-08-01  2021-09-20  2021-09-20  North East Scotland   \n",
       "\n",
       "  agglomeration local_authority  \n",
       "0           NaN   Aberdeen City  "
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "final_sites = meta.get_site(site_id='ABD')\n",
    "final_sites.drop_duplicates(subset='site_id')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Start to merge the selected station data\n",
    "- Use \"start\" and \"end\" to define a period of time.\n",
    "- Other parameters are the same as \"Download data by a bounds of longitude and latitude\" part"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "   site_id site_name     location_type  latitude  longitude parameter  \\\n",
      "0      ABD  Aberdeen  Urban Background  57.15736  -2.094278        O3   \n",
      "1      ABD  Aberdeen  Urban Background  57.15736  -2.094278        NO   \n",
      "2      ABD  Aberdeen  Urban Background  57.15736  -2.094278       NO2   \n",
      "3      ABD  Aberdeen  Urban Background  57.15736  -2.094278  NOXasNO2   \n",
      "6      ABD  Aberdeen  Urban Background  57.15736  -2.094278      PM10   \n",
      "7      ABD  Aberdeen  Urban Background  57.15736  -2.094278      NV10   \n",
      "8      ABD  Aberdeen  Urban Background  57.15736  -2.094278       V10   \n",
      "9      ABD  Aberdeen  Urban Background  57.15736  -2.094278     PM2.5   \n",
      "10     ABD  Aberdeen  Urban Background  57.15736  -2.094278     NV2.5   \n",
      "11     ABD  Aberdeen  Urban Background  57.15736  -2.094278      V2.5   \n",
      "12     ABD  Aberdeen  Urban Background  57.15736  -2.094278        wd   \n",
      "13     ABD  Aberdeen  Urban Background  57.15736  -2.094278        ws   \n",
      "14     ABD  Aberdeen  Urban Background  57.15736  -2.094278      temp   \n",
      "\n",
      "                                Parameter_name  start_date    end_date  \\\n",
      "0                                        Ozone  2003-08-01  2021-09-20   \n",
      "1                                 Nitric oxide  1999-09-18  2021-09-20   \n",
      "2                             Nitrogen dioxide  1999-09-18  2021-09-20   \n",
      "3          Nitrogen oxides as nitrogen dioxide  1999-09-18  2021-09-20   \n",
      "6    PM10 particulate matter (Hourly measured)  1999-09-18  2021-09-20   \n",
      "7          Non-volatile PM10 (Hourly measured)  2009-02-20  2019-11-07   \n",
      "8              Volatile PM10 (Hourly measured)  2009-02-20  2019-11-07   \n",
      "9   PM2.5 particulate matter (Hourly measured)  2009-02-20  2021-09-20   \n",
      "10        Non-volatile PM2.5 (Hourly measured)  2009-02-20  2019-11-07   \n",
      "11            Volatile PM2.5 (Hourly measured)  2009-02-20  2019-11-07   \n",
      "12                     Modelled Wind Direction  2010-08-01  2021-09-20   \n",
      "13                         Modelled Wind Speed  2010-08-01  2021-09-20   \n",
      "14                        Modelled Temperature  2010-08-01  2021-09-20   \n",
      "\n",
      "   ratified_to                 zone agglomeration local_authority  \n",
      "0   2021-09-20  North East Scotland           NaN   Aberdeen City  \n",
      "1   2021-09-20  North East Scotland           NaN   Aberdeen City  \n",
      "2   2021-09-20  North East Scotland           NaN   Aberdeen City  \n",
      "3   2021-09-20  North East Scotland           NaN   Aberdeen City  \n",
      "6   2021-09-20  North East Scotland           NaN   Aberdeen City  \n",
      "7   2019-11-07  North East Scotland           NaN   Aberdeen City  \n",
      "8   2019-11-07  North East Scotland           NaN   Aberdeen City  \n",
      "9   2021-09-20  North East Scotland           NaN   Aberdeen City  \n",
      "10  2019-11-07  North East Scotland           NaN   Aberdeen City  \n",
      "11  2019-11-07  North East Scotland           NaN   Aberdeen City  \n",
      "12       Never  North East Scotland           NaN   Aberdeen City  \n",
      "13       Never  North East Scotland           NaN   Aberdeen City  \n",
      "14       Never  North East Scotland           NaN   Aberdeen City  \n",
      "Saved merged file: data/test_pm25_siteid/MERGED_aurn_PM2.5_20171201_20181130.csv\n"
     ]
    },
    {
     "data": {
      "application/vnd.microsoft.datawrangler.viewer.v0+json": {
       "columns": [
        {
         "name": "index",
         "rawType": "int64",
         "type": "integer"
        },
        {
         "name": "site_id",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "Date",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "time",
         "rawType": "object",
         "type": "string"
        },
        {
         "name": "PM<sub>2.5</sub> particulate matter (Hourly measured)",
         "rawType": "float64",
         "type": "float"
        },
        {
         "name": "status.7",
         "rawType": "object",
         "type": "unknown"
        },
        {
         "name": "unit.7",
         "rawType": "object",
         "type": "unknown"
        }
       ],
       "ref": "e3452f15-4185-4d58-aacb-af4953cb02a1",
       "rows": [
        [
         "0",
         "ABD",
         "30-11-2017",
         "24:00",
         "2.2",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "1",
         "ABD",
         "01-12-2017",
         "01:00",
         "2.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "2",
         "ABD",
         "01-12-2017",
         "02:00",
         "3.2",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "3",
         "ABD",
         "01-12-2017",
         "03:00",
         "4.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "4",
         "ABD",
         "01-12-2017",
         "04:00",
         "2.4",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "5",
         "ABD",
         "01-12-2017",
         "05:00",
         "3.0",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "6",
         "ABD",
         "01-12-2017",
         "06:00",
         "3.0",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "7",
         "ABD",
         "01-12-2017",
         "07:00",
         "6.4",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "8",
         "ABD",
         "01-12-2017",
         "08:00",
         "3.3",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "9",
         "ABD",
         "01-12-2017",
         "09:00",
         "3.7",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "10",
         "ABD",
         "01-12-2017",
         "10:00",
         "2.8",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "11",
         "ABD",
         "01-12-2017",
         "11:00",
         "4.9",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "12",
         "ABD",
         "01-12-2017",
         "12:00",
         "4.3",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "13",
         "ABD",
         "01-12-2017",
         "13:00",
         "4.2",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "14",
         "ABD",
         "01-12-2017",
         "14:00",
         "1.4",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "15",
         "ABD",
         "01-12-2017",
         "15:00",
         "4.9",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "16",
         "ABD",
         "01-12-2017",
         "16:00",
         "4.3",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "17",
         "ABD",
         "01-12-2017",
         "17:00",
         "3.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "18",
         "ABD",
         "01-12-2017",
         "18:00",
         "4.5",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "19",
         "ABD",
         "01-12-2017",
         "19:00",
         "4.4",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "20",
         "ABD",
         "01-12-2017",
         "20:00",
         "3.5",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "21",
         "ABD",
         "01-12-2017",
         "21:00",
         "2.4",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "22",
         "ABD",
         "01-12-2017",
         "22:00",
         "0.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "23",
         "ABD",
         "01-12-2017",
         "23:00",
         "-0.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "24",
         "ABD",
         "01-12-2017",
         "24:00",
         "6.2",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "25",
         "ABD",
         "02-12-2017",
         "01:00",
         "2.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "26",
         "ABD",
         "02-12-2017",
         "02:00",
         "1.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "27",
         "ABD",
         "02-12-2017",
         "03:00",
         "2.3",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "28",
         "ABD",
         "02-12-2017",
         "04:00",
         "1.2",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "29",
         "ABD",
         "02-12-2017",
         "05:00",
         "0.7",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "30",
         "ABD",
         "02-12-2017",
         "06:00",
         "1.8",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "31",
         "ABD",
         "02-12-2017",
         "07:00",
         "1.4",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "32",
         "ABD",
         "02-12-2017",
         "08:00",
         "-0.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "33",
         "ABD",
         "02-12-2017",
         "09:00",
         "3.7",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "34",
         "ABD",
         "02-12-2017",
         "10:00",
         "2.6",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "35",
         "ABD",
         "02-12-2017",
         "11:00",
         "5.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "36",
         "ABD",
         "02-12-2017",
         "12:00",
         "3.2",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "37",
         "ABD",
         "02-12-2017",
         "13:00",
         "0.9",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "38",
         "ABD",
         "02-12-2017",
         "14:00",
         "2.0",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "39",
         "ABD",
         "02-12-2017",
         "15:00",
         "3.9",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "40",
         "ABD",
         "02-12-2017",
         "16:00",
         "2.9",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "41",
         "ABD",
         "02-12-2017",
         "17:00",
         "3.9",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "42",
         "ABD",
         "02-12-2017",
         "18:00",
         "5.0",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "43",
         "ABD",
         "02-12-2017",
         "19:00",
         "2.8",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "44",
         "ABD",
         "02-12-2017",
         "20:00",
         "5.9",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "45",
         "ABD",
         "02-12-2017",
         "21:00",
         "2.2",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "46",
         "ABD",
         "02-12-2017",
         "22:00",
         "0.1",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "47",
         "ABD",
         "02-12-2017",
         "23:00",
         "3.3",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "48",
         "ABD",
         "02-12-2017",
         "24:00",
         "2.7",
         "R",
         "ugm-3 (TEOM FDMS)"
        ],
        [
         "49",
         "ABD",
         "03-12-2017",
         "01:00",
         "0.7",
         "R",
         "ugm-3 (TEOM FDMS)"
        ]
       ],
       "shape": {
        "columns": 6,
        "rows": 8760
       }
      },
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>site_id</th>\n",
       "      <th>Date</th>\n",
       "      <th>time</th>\n",
       "      <th>PM&lt;sub&gt;2.5&lt;/sub&gt; particulate matter (Hourly measured)</th>\n",
       "      <th>status.7</th>\n",
       "      <th>unit.7</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>ABD</td>\n",
       "      <td>30-11-2017</td>\n",
       "      <td>24:00</td>\n",
       "      <td>2.2</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>ABD</td>\n",
       "      <td>01-12-2017</td>\n",
       "      <td>01:00</td>\n",
       "      <td>2.1</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>ABD</td>\n",
       "      <td>01-12-2017</td>\n",
       "      <td>02:00</td>\n",
       "      <td>3.2</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>ABD</td>\n",
       "      <td>01-12-2017</td>\n",
       "      <td>03:00</td>\n",
       "      <td>4.1</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>ABD</td>\n",
       "      <td>01-12-2017</td>\n",
       "      <td>04:00</td>\n",
       "      <td>2.4</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8755</th>\n",
       "      <td>ABD</td>\n",
       "      <td>30-11-2018</td>\n",
       "      <td>19:00</td>\n",
       "      <td>3.0</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8756</th>\n",
       "      <td>ABD</td>\n",
       "      <td>30-11-2018</td>\n",
       "      <td>20:00</td>\n",
       "      <td>5.6</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8757</th>\n",
       "      <td>ABD</td>\n",
       "      <td>30-11-2018</td>\n",
       "      <td>21:00</td>\n",
       "      <td>3.0</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8758</th>\n",
       "      <td>ABD</td>\n",
       "      <td>30-11-2018</td>\n",
       "      <td>22:00</td>\n",
       "      <td>0.3</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8759</th>\n",
       "      <td>ABD</td>\n",
       "      <td>30-11-2018</td>\n",
       "      <td>23:00</td>\n",
       "      <td>5.7</td>\n",
       "      <td>R</td>\n",
       "      <td>ugm-3 (TEOM FDMS)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>8760 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "     site_id        Date   time  \\\n",
       "0        ABD  30-11-2017  24:00   \n",
       "1        ABD  01-12-2017  01:00   \n",
       "2        ABD  01-12-2017  02:00   \n",
       "3        ABD  01-12-2017  03:00   \n",
       "4        ABD  01-12-2017  04:00   \n",
       "...      ...         ...    ...   \n",
       "8755     ABD  30-11-2018  19:00   \n",
       "8756     ABD  30-11-2018  20:00   \n",
       "8757     ABD  30-11-2018  21:00   \n",
       "8758     ABD  30-11-2018  22:00   \n",
       "8759     ABD  30-11-2018  23:00   \n",
       "\n",
       "      PM<sub>2.5</sub> particulate matter (Hourly measured) status.7  \\\n",
       "0                                                   2.2            R   \n",
       "1                                                   2.1            R   \n",
       "2                                                   3.2            R   \n",
       "3                                                   4.1            R   \n",
       "4                                                   2.4            R   \n",
       "...                                                 ...          ...   \n",
       "8755                                                3.0            R   \n",
       "8756                                                5.6            R   \n",
       "8757                                                3.0            R   \n",
       "8758                                                0.3            R   \n",
       "8759                                                5.7            R   \n",
       "\n",
       "                 unit.7  \n",
       "0     ugm-3 (TEOM FDMS)  \n",
       "1     ugm-3 (TEOM FDMS)  \n",
       "2     ugm-3 (TEOM FDMS)  \n",
       "3     ugm-3 (TEOM FDMS)  \n",
       "4     ugm-3 (TEOM FDMS)  \n",
       "...                 ...  \n",
       "8755  ugm-3 (TEOM FDMS)  \n",
       "8756  ugm-3 (TEOM FDMS)  \n",
       "8757  ugm-3 (TEOM FDMS)  \n",
       "8758  ugm-3 (TEOM FDMS)  \n",
       "8759  ugm-3 (TEOM FDMS)  \n",
       "\n",
       "[8760 rows x 6 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "meta.download_sites(\n",
    "    port=\"aurn\",\n",
    "    pollutant=\"PM2.5\",\n",
    "    start=\"2017-12-01\",\n",
    "    end=\"2018-11-30\",\n",
    "    output_dir=\"data/test_pm25_siteid\",\n",
    "    download_mode=\"memory\",     \n",
    "    save_per_site=False,\n",
    "    save_merged=True,\n",
    "    add_site_id=True\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Analyse the data and draw plots"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Temporal trend of the air pollutant concentration.\n",
    "After data filtering and preprocessing, we can do various analysis based on the data. Here, we use a trend analysis to examine how PM2.5 concentration varies over time."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABW4AAAHqCAYAAACUWtfDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAADJaklEQVR4nOzdBZQcVdbA8TtxEpLggSwBQnBncXcCARYJy+JBFltgP8hiwQkSnGAhuAeCBIknxN3d3d0mLjP9nVsz1VNd01LVXdVV3fP/nTOZTE9Pd3Xpq/vuu68gEolEBAAAAAAAAAAQGpWCXgAAAAAAAAAAQCwCtwAAAAAAAAAQMgRuAQAAAAAAACBkCNwCAAAAAAAAQMgQuAUAAAAAAACAkCFwCwAAAAAAAAAhQ+AWAAAAAAAAAEKGwC0AAAAAAAAAhAyBWwAAAAAAAAAIGQK3AAAACJUDDjhAbrvtNgmLuXPnSkFBgXz11Vd5/TnjmTFjhlx88cVSt25dYx38/vvvQS8SAABAhUHgFgAAIIs0+OXkq2/fvkEvKjwwePBgef7552Xt2rWSi5o1ayYTJkyQl19+Wb799ls58cQTkwa3za/KlSvLfvvtJ1dffbWMHTs25rnmc/7973/Hfa2nnnoq+pyVK1dGH+/QoYP861//kgMPPFBq1qwphx56qPzvf/9zvG41SB7vWDvssMNcrRMAAIBsqZK1dwIAAIAR/LL65ptvpGfPnuUeP/zww7O8ZEhk//33l82bN0vVqlXTCty+8MILRtBwl112ifndtGnTpFKl8OZR6GceMmSIEUh94IEHHP3NDTfcIE2aNJGioiKZMmWKfPTRR9K1a1cZOnSoHHfccdHn1ahRQ3799Vdp06aNVKtWLeY1fvjhB+P3W7ZsiXn87rvvlvr168vNN99sBIU1oPzBBx9Ily5dZPTo0bLTTjulXL7q1avLZ599FvOYZhMDAACEEYFbAACALNKgk5UGtDRwa388XxQXF8u2bduMQFyu2bFjh7H8Glj0Y/k1iBhmK1asML7bA87J/P3vf4/Zl8844wz5xz/+YQRwP/744+jjl1xyifz5559GUPfKK6+MCXTPmTNHmjZtagR2rX755Rc599xzYx474YQTjKzg77//PmEGr1WVKlXy9lgDAAD5J7xd/AAAABWUBgtbt24tRx55pBEwrFevntxzzz2yZs2acjVSL7/8cqOsgg5h14zDo48+OlpmQYeW68/6GhrgGjNmTMzfaxbozjvvLLNnz5bGjRtLrVq1jIzGli1bSiQSiXnuxo0bjWHpDRo0MAKOOkz9zTffLPc8HXqu2ZkaSNPl1+d269bN+J0+//TTT5fdd9/dWFZdJg3GubV9+3bZbbfd5Pbbby/3u8LCQuPzPvLII8bPGjR+9tlnjffSzEr9jGeddZb06dMn7lB/XUZd940aNTKWffLkyXFr3I4fP95YfzpsX99v7733ljvuuENWrVoVfY6WSHj00UeN/zds2DA6NF9fL1GNW90W//znP43Pp+UATj31VOncuXPMc3T76uv89NNPRgmDfffd11iGCy64QGbOnOloHeq+cOmll0qdOnWMfUD/VjsRrMuumcZKP4O+ny6vW+eff77xXYOxVn/729/k7LPPlnbt2sU8rvuN7rNHHXVUudeyB22VlmJQmt3rlGYD634CAAAQdmTcAgAAhIwGaTVIqIHJ//73v0bQS4eEa7Bt0KBBMUP2NVB34403Gn+jmYQaeLziiiukbdu28uSTT8p//vMf43mtWrWS6667rtzwfA1iafajBghff/11I8j63HPPGdmmGsBVGpzVrEkNdt55553GkPfu3bsbAb1FixbJO++8E7P8vXv3NoKKGsDdY489ogG/d99913idm266yQio/vjjj0aQslOnTnLZZZc5Xj/6+TVgp4FpzeK0DrXXybO2bt0q119/vfGzBuh0aLwO4b/rrrtk/fr18vnnnxuB6uHDh8cM31dffvmlMURfh+Vr4FYDqBpIt9MsaQ2y6jbSoO2kSZPkk08+Mb5rAFQDnddcc41Mnz7dGPqv60jXhdpzzz3jfq5ly5YZge1NmzYZ210D3F9//bWxzjTAbQYpTa+++qqxLTVIvW7dOmP76bodNmxY0vWny6jBaw3aPvbYY8b61PWogdF+/frJKaecYiy7Zto+/PDD0fIHGuB1a9asWcZ3/Sx2ut/+3//9n2zYsMF4bd3nfv75Z2nevHm5MgmJLF261PhurttUdN3q59bvu+66q/HZXnvttbQ+GwAAgO8iAAAACMz999+vKavRnwcMGGD8/P3338c8r1u3buUe33///Y3HBg8eHH2se/fuxmM77bRTZN68edHHP/74Y+PxPn36RB9r1qyZ8diDDz4Yfay4uDhy2WWXRapVqxZZsWKF8djvv/9uPO+ll16KWaZrr702UlBQEJk5c2b0MX1epUqVIpMmTSr3WTdt2hTz87Zt2yJHHXVU5Pzzz495XD+XLlsy5ufs2LFjzONNmjSJHHjggdGfd+zYEdm6dWvMc9asWROpV69e5I477og+NmfOHOP16tSpE1m+fHnM883fffnllwk/i/rhhx+M5/Xv3z/62BtvvGE8pq9hZ/+cDz30kPFc3QdM69evjzRs2DBywAEHRIqKiozHdBvq8w4//PCYz/buu+8aj0+YMCHpurvqqquM7Ttr1qzoY4sXL47Url07cvbZZ5f73PoZUjGf+8ILLxj7zdKlSyN9+/aNHH/88cbjv/76a/S5+rPu96tXrzaW49tvvzUe79y5s7E/zZ07N/Lcc88ZzzP3wUTuvPPOSOXKlSPTp09PuYxPPPFE5PHHH4+0b9/e2Fbm/n/GGWdEtm/fnvLvAQAAso1SCQAAACGiGYc6pP+iiy6SlStXRr90qL9mBdqH+B9xxBFy2mmnRX/WbElziLpO4GR/XLNE7awTT5mlDjQj9q+//jIe08mfKleubGSBWmnpBI3DaZ1Sq3POOcdYLjvr5FFa9kGzRDXzUyeWcks/n2ZZtm/fPuY1NRP2X//6V/QxXW4zI1czZ1evXm1kdmppiXjvq7VVE2XEJvosmh2q20izllU6n8dczyeffLKceeaZ0cd0m2v2r5ZX0LINVprta8021nWZaBtbM6x79OghV111lVHmwbTPPvsYGbADBw7MqIyAZmvr+tMsZM3g1YxbzWjVDF47zXjVbG/NSFZaNkEzjs0SDano8zV7WvfDgw8+OOXzNetcs5Q181wzsjWrXUtNaBZ7OiU7AAAA/EbgFgAAIERmzJhhBDT32msvIwBm/dIh5cuXL495vjU4qzToq7QWbbzH7XVydai9NYCnDjnkEOO7WYt13rx5Ru3b2rVrxzzv8MMPj/7eSuu5xqMlETS4qfVYtQSBfiadtEo/r1s6yZQGWf/44w+jNILS0gla/9YauFVabuCYY44x3leH7Ov7at3YeO+baNntNACsw/y1/rAGcfU1zb9N5/OY61FrB9slWs/2ba+B0Hjb2D7hmJYJSPQ+GtxesGCBpEuDzBo879Wrl4waNcrYX7UcQyIaLNbnz58/3yhzoT87MWDAAKNsh5a80OBrurQUhB4DZicFAABAmFDjFgAAIEQ0cKZBW52kKR57NqhmlMaT6HH7ZGJ+sGajWgNtWqtVJ6Rq06aNkeGptVW1pqx9giqnNGtSa7Nqxq9mkGpd3cMOO0yOPfbY6HO+++47YwIw/b3W5NV1q+tGsy/N+quplj0ezdocPHiw8ZpaJ1czY3XbaQZpvJq4fghyGyeima8XXnih4+frPqG1hJs1a2YE4HW9pjJu3Djj73QCM82U1SB+unR7azBfA/EAAABhQ+AWAAAgRBo1amRk/51xxhmOg4iZ0CCjDq03s2yVTqilzEnFdOi6LpNO7GXNup06dWr096n8+uuvRsarTmqmgTqTBm7TpUFgDQBruQQtL6CToj311FMxz9HAnmYUazauloGwDulPl2a0akbpCy+8IM8++2xMtrSd9T1T0fWok8fZuVnPqWjgv2bNmgnfR7NP7dnaftJ9XIPqGmC/9NJLU04ypsF2DY5rAF5LS2Q6qZju01rmwkl5DAAAgGyjVAIAAECIaMah1iF98cUXy/1Oa7OuXbvW8/f84IMPYrI19WfNhr3ggguMx5o0aWIsk/V56p133jECkxpwc5Idqs/V1zFpKQYdHp8uDTJee+210rFjR/n222+N9WMvk2BmpVqzUIcNGyZDhgxJ+33jvaZq3bp1uefWqlXL+O5ku+l6Hj58eMyybdy4UT755BMjiB6vbnA6y37xxRcbJSbMUhhq2bJlRuazBsDr1Kkj2fTII48YgfRnnnkm6fOWLl1qLLtud+0ASBZs1QCvNaNa6xBrkNZOjzPdjhoMBgAACBsybgEAAEJEJ/a65557jKH8Y8eONQJVGkTVbE6duOzdd981gpVe0SzYbt26GUPVdQIzLTug9V+ffPLJaGDsiiuukPPOO8/IZtVgn5Yi0AmuNPj30EMPGVnCqVx22WXy9ttvGwEyrWOqtU8//PBDOeigg2T8+PFpL78Gat9//30j8Hf00UdH68GaLr/8ciPb9uqrrzaWYc6cOdK2bVsjCKo1g9OhgU3N9n399deNmrp/+9vfjPWhr22nk8opXXda2kG3pa5PM6Br9cQTTxgTdWkgXCeC0zrAWp9XX1czljVg6YWXXnrJqCurQdr//Oc/RqkBLTmhpQr0M2Wb7k/W8haJ6L6j2eFaM1cnUdMvk9Ya1gn9TGangxmc1qDv8ccfLzfccINRTkNp8FezdvV1r7zySh8+GQAAQGYI3AIAAISMBhY14KfBNA2gamBNMy5vvvlmo4SClzQDUwO39913n1GvVUshaBDUWgJAA4Z//vmn8ZiWJdDyBro8b7zxhvzvf/9z9D7nn3++fP755/Lqq68awV6dyOu1114zAmuZBG5PP/10Y2i/Tqhlz7ZVWt9Wg3a6LjVQpwFbHZavQfC+ffum/b6anfrggw8awWfN2NQAuwa9dRI3q5NOOsnI6tRtqutZS1NoIDZe4FaDj1o39/HHHzeC0ZolqpOqaUaxBp29cuSRRxo1h1u0aGF0EOgyadBe14t+DyutbaviBZe1w8MauLXbZZddjCC+Bqw1GK6Z39pp8MorrxgZv14FxQEAALxUEAly9gIAAAAERoOaWgM23cxTAAAAAP6haxkAAAAAAAAAQobALQAAAAAAAACEDIFbAAAAAAAAAAgZatwCAAAAAAAAQMiQcQsAAAAAAAAAIUPgFgAAAAAAAABCpkrQCxAGxcXFsnjxYqldu7YUFBQEvTgAAAAAAAAA8pBWrV2/fr3Ur19fKlVKkVMbCVCbNm0iRx99dKR27drG16mnnhrp0qVL9PfnnHOO1t+N+brnnntiXmPevHmRJk2aRHbaaafInnvuGXnkkUci27dvd7UcCxYsKPc+fPHFF1988cUXX3zxxRdffPHFF1988cUXX+LDl8YjUwk043bfffeVV199VQ4++GAj2vz111/LlVdeKWPGjJEjjzzSeM5dd90lLVu2jP5NzZo1o/8vKiqSyy67TPbee28ZPHiwLFmyRG699VapWrWqvPLKK46XQzNt1YIFC6ROnTqefkYAAAAAAAAAUIWFhdKgQYNoPDKZAo3eSojstttu8sYbb8idd94p5557rhx33HHSunXruM/t2rWrXH755UaZg3r16hmPtW3bVh5//HFZsWKFVKtWzfEKq1u3rqxbt47ALQAAAAAAAABfuIlDhmZyMs2e/fHHH2Xjxo1y2mmnRR///vvvZY899pCjjjpKWrRoIZs2bYr+bsiQIXL00UdHg7aqcePGxgqYNGlS1j8DAAAAAAAAAOTF5GQTJkwwArVbtmyRnXfeWX777Tc54ogjjN/deOONsv/++xvFesePH29k0k6bNk06dOhg/H7p0qUxQVtl/qy/S2Tr1q3Gl0kDvQAAAAAAAAAQFoEHbg899FAZO3askR78yy+/SLNmzaRfv35G8Pbuu++OPk8za/fZZx+54IILZNasWdKoUaO037NVq1bywgsvePQJAAAAAAAAAMBbgZdK0Dq0Bx10kJxwwglGQPXYY4+Vd999N+5zTznlFOP7zJkzje86KdmyZctinmP+rL9LREsuaKDY/NJJyQAAAAAAAAAgLAIP3NoVFxfHlDGw0sxcpZm3SkssaKmF5cuXR5/Ts2dPo7CvWW4hnurVqxvPsX4BAAAAAAAAQFgEWipBM18vvfRS2W+//WT9+vXSrl076du3r3Tv3t0oh6A/N2nSRHbffXejxu3DDz8sZ599thxzzDHG31988cVGgPaWW26R119/3ahr+/TTT8v9999vBGcBAAAAAAAAIBcFGrjVTNlbb71VlixZInXr1jUCshq0veiii4zyBX/99Ze0bt1aNm7cKA0aNJCmTZsagVlT5cqVpVOnTnLfffcZ2be1atUyauS2bNkyyI8FAAAAAAAAABkpiEQiEangCgsLjcCx1rulbAIAAAAAAACAoOOQoatxCwAAAAAAAAAVHYFbAAAAAAAAAAgZArcAAAAAAAAAEDIEbgEAAAAAAAAgZAjcAgAAAAAAAEDIELgFAAAAAAAAgJAhcAsAAAAAAAAAIUPgFgAAAAAAAHDppU6T5azXewe9GMhjVYJeAAAAAAAAACDXfDZwTtCLgDxHxi0AAAAAAAAAhAyBWwAAAAAAAAAIGQK3AAAAAAAAABAyBG4BAAAAAAAAIGQI3AIAAAAAAABAyBC4BQAAAAAAAICQIXALAAAAAAAAACFD4BYAAAAAAAAAQobALQAAAAAAAACEDIFbAAAAAAAAAAgZArcAAAAAAAAAEDIEbgEAAAAAAAAgZAjcAgAAAAAAAEDIELgFAAAAAAAAgJAhcAsAAAAAAAAAIUPgFgAAAAAAAABChsAtAAAAAAAAAIQMgVsAAAAAAAAACBkCtwAAAAAAAAAQMgRuAQAAAAAAACBkCNwCAAAAAAAAQMgQuAUAAAAAAACAkCFwCwAAAAAAAAAhQ+AWAAAAAAAAAEKGwC0AAAAAAAAAhAyBWwAAAAAAAAAIGQK3AAAAAAAAABAyBG4BAAAAAAAAIGQI3AIAAAAAAABAyBC4BQAAAAAAAICQIXALAAAAAAAAACFD4BYAAAAAAAAAQobALQAAAAAAAACEDIFbAAAAAAAAAAiZQAO3H330kRxzzDFSp04d4+u0006Trl27Rn+/ZcsWuf/++2X33XeXnXfeWZo2bSrLli2LeY358+fLZZddJjVr1pS99tpLHn30UdmxY0cAnwYAAAAAAAAA8iBwu++++8qrr74qo0aNkpEjR8r5558vV155pUyaNMn4/cMPPywdO3aUn3/+Wfr16yeLFy+Wa665Jvr3RUVFRtB227ZtMnjwYPn666/lq6++kmeffTbATwUAAAAAAAAAmSmIRCIRCZHddttN3njjDbn22mtlzz33lHbt2hn/V1OnTpXDDz9chgwZIqeeeqqRnXv55ZcbAd169eoZz2nbtq08/vjjsmLFCqlWrZqj9ywsLJS6devKunXrjMxfAAAAAAAAIJkDnuhsfJ/76mVBLwpyiJs4ZGhq3Gr27I8//igbN240SiZoFu727dvlwgsvjD7nsMMOk/32288I3Cr9fvTRR0eDtqpx48bGCjCzduPZunWr8RzrFwAAAAAAAACEReCB2wkTJhj1a6tXry733nuv/Pbbb3LEEUfI0qVLjYzZXXbZJeb5GqTV3yn9bg3amr83f5dIq1atjMi2+dWgQQNfPhsAAAAAAAAA5GTg9tBDD5WxY8fKsGHD5L777pNmzZrJ5MmTfX3PFi1aGOnI5teCBQt8fT8AAAAAAAAAcKOKBEyzag866CDj/yeccIKMGDFC3n33XfnXv/5lTDq2du3amKzbZcuWyd577238X78PHz485vX09+bvEtHsXv0CAAAAAAAAgDAKPOPWrri42KhBq0HcqlWrSq9evaK/mzZtmsyfP9+ogav0u5ZaWL58efQ5PXv2NAr7arkFAAAAAAAAAMhFgWbcasmCSy+91JhwbP369dKuXTvp27evdO/e3ag9e+edd0rz5s1lt912M4KxDz74oBGsPfXUU42/v/jii40A7S233CKvv/66Udf26aeflvvvv5+MWgAAAAAAAAA5K9DArWbK3nrrrbJkyRIjUHvMMccYQduLLrrI+P0777wjlSpVkqZNmxpZuI0bN5Y2bdpE/75y5crSqVMnozauBnRr1apl1Mht2bJlgJ8KAAAAAAAA7UfMl8d/nSBTX7xEalStHPTiADkn0MDt559/nvT3NWrUkA8//ND4SmT//feXLl26+LB0AAAAAAAASFe3iUuN71u3FxO4BfKhxi0AAAAAAAAAVHQEbgEAAAAAAAAgZAjcAgAAAAAAAEDIELgFAAAAAAAAgJAhcAsAAAAAAAAAIUPgFgAAAAAAAABChsAtAATk+T8nycf9ZgW9GAAAAAAAIISqBL0AAFBRfTV4rvH9nnMaBb0oAAAAAAAgZMi4BQAAAAAAAICQIXALAAAAAAAAACFD4BYAAAAAAAAAQobALQAAAAAAAACEDIFbAAAAAAAAAAgZArcAAAAAAAAAEDIEbgEAAAAAAOC5SNALAOQ4ArcAAAAAAAAAEDIEbgEAAAAAAOC5gqAXAMhxBG4BAAAAAACANEUiFIWAPwjcAgAAAAAAAEDIELgFAAAAAAAA0kTCLfxC4BYAAAAAAAAAQobALQAAAAAAAJAmEm7hFwK3AAAAAAAA8E2E0CaQFgK3AAAAAAAAQJrmrtoY9CIgTxG4BQAAAAAAgG8KpEDyWeHm7UEvAvIUgVsAAAAAAAAACBkCtwAAAAAAAAAQMgRuAQAAAAAAACBkCNwCAAAAAAAAaYoEvQDIWwRuAQAAAAAAgDRFIoRu4Q8CtwAAAAAAAPBNhJxUIC0EbgEAAAAAAIA0kXALvxC4BQAAAAAAgG8KpEDyWTGBW/iEwC0AAAAAAAAAhAyBWwAAAAAAACBNTE4GvxC4BQAAAAAAANJEqQT4hcAtAAAAAAAAkCYybuEXArcAAAAAAABAmgjbwi8EbgEAAAAAAOCbSJ6HNkm4hV8I3AIAAAAAAABpyvfANIJD4BYAAAAAAAC+KZACyWdMTga/ELgFAAAAAAAAgJAhcAsAAAAAAACkKUKRW+Rj4LZVq1Zy0kknSe3atWWvvfaSq666SqZNmxbznHPPPVcKCgpivu69996Y58yfP18uu+wyqVmzpvE6jz76qOzYsSPLnwYAAAAAAAAVDXFb+KWKBKhfv35y//33G8FbDbQ++eSTcvHFF8vkyZOlVq1a0efddddd0rJly+jPGqA1FRUVGUHbvffeWwYPHixLliyRW2+9VapWrSqvvPJK1j8TAAAAAAAAKg4mJ0NeBm67desW8/NXX31lZMyOGjVKzj777JhArQZm4+nRo4cR6P3rr7+kXr16ctxxx8mLL74ojz/+uDz//PNSrVo13z8HAAAAAAAAKqbi4qCXAPkqVDVu161bZ3zfbbfdYh7//vvvZY899pCjjjpKWrRoIZs2bYr+bsiQIXL00UcbQVtT48aNpbCwUCZNmhT3fbZu3Wr83voFAAAAAAAA7+V7Rmp+fzpU2Ixbq+LiYnnooYfkjDPOMAK0phtvvFH2339/qV+/vowfP97IpNU6uB06dDB+v3Tp0pigrTJ/1t8lqq37wgsv+Pp5AAAAAAAAkP+YnAx5H7jVWrcTJ06UgQMHxjx+9913R/+vmbX77LOPXHDBBTJr1ixp1KhRWu+lWbvNmzeP/qwZtw0aNMhg6QEAAAAAABBPgRRIPiNsi7wulfDAAw9Ip06dpE+fPrLvvvsmfe4pp5xifJ85c6bxXWvfLlu2LOY55s+J6uJWr15d6tSpE/MFAAAAAAAAuEXCLfIycKup5Bq0/e2336R3797SsGHDlH8zduxY47tm3qrTTjtNJkyYIMuXL48+p2fPnkYw9ogjjvBx6QEAAAAAAAAit8jDUglaHqFdu3byxx9/SO3ataM1aevWrSs77bSTUQ5Bf9+kSRPZfffdjRq3Dz/8sJx99tlyzDHHGM+9+OKLjQDtLbfcIq+//rrxGk8//bTx2ppZCwAAAAAAAPiFjFvkZcbtRx99JOvWrZNzzz3XyKA1v9q3b2/8vlq1avLXX38ZwdnDDjtM/ve//0nTpk2lY8eO0deoXLmyUWZBv2v27c033yy33nqrtGzZMsBPBgAAAAAAgIqgmMAt8jHjNtWsezphWL9+/VK+zv777y9dunTxcMkAAAAAAACA1CKUSkA+T04GAAAAAACA/JTvgU1KJcAvBG4BAAAAAACANBUTuYVPCNwCAAAAAADANwVSEPQiADmJwC0AAAAAAAAAhAyBWwAAAAAAACBNVEqAXwjcAgAAAAAAAGnK98nXEBwCtwAAAAAAAECaiouDXgLkKwK3AAAAAAAAQJqKqZUAnxC4BQAAAAAAgG/yvZRAfn86BInALQAAAAAAADxXUQKaETJu4RMCtwAAAAAAAPBNgRRIPismbgufELgFAAAAAACA5/I7XFuGhFv4hcAtAAAAAAAAkCYmJ4NfCNwCAAAAAAAAaaLGLfxC4BYAAAAAAABIUxFFbuETArcAAAAAAABAmoqI28InBG4BAAAAAADgm4jkd2STUgnwC4FbAAAAAAAAeK6ihDOZnAx+IXALAAAAAAAApIkSt/ALgVsAAAAAAAB4rkAqBjJu4RcCtwAAAAAAAECaiNvCLwRuAQAAAITGwBkrpZgxpwCAHFLEdQs+IXALAAAAIBTGLVgrN38+TL4fNi/oRQEAwDECt/ALgVsAAABARA54orN82n920ItRoa3bvN34vnz91qAXBQAAxyLUSoBPCNwCAAAApT4bSOAWAACv5Xtck4Rb+IXALQAAAAAAAJCm4nyPTCMwBG4BAAAAhIJ521sQ8HIAAOBGEYFb+ITALQAAAAAAAJAu4rbwCYFbAAAAAAAAIE2USoBfCNwCAAAACJcCiiUAAHIHk5PBLwRuAQAAAAAAgDQVEbmFTwjcAgAAAAAAAGmKUCoBPiFwCwAAACAUuPGF6bFfxsmpr/QKejEAeCTfz+75/vkQnCoBvjcAAAAAlLvxpcItfhq5MOhFAADHKJUAv5BxCwAAAJQqIGQIAABcKmbECHxC4BYAAAAoFWGwIwAAcKm4OOglQL5yVSqhuLhY+vXrJwMGDJB58+bJpk2bZM8995Tjjz9eLrzwQmnQoIF/SwoAAACgQigg8RkAkEPIuEWgGbebN2+Wl156yQjMNmnSRLp27Spr166VypUry8yZM+W5556Thg0bGr8bOnSobwsLAAAAAAAqhk/7z5aFazYFvRhASpS4RaAZt4cccoicdtpp8umnn8pFF10kVatWLfcczcBt166dXH/99fLUU0/JXXfd5cfyAgAAAL6hxm3AuPEFYPFylynSacIS+eP+M4JeFCApMm4RaOC2R48ecvjhhyd9zv777y8tWrSQRx55RObPn+/V8gEAAACoYAigAzAVUTwUOYDALQItlZAqaGul2biNGjXKZJkAAAAAAACQJyJ5HtikVAICDdyqQYMGGbVsR40aJddcc43079/ft4UCAAAAghBhrD4AhEqex/uQJ8i4ReCB2xdeeEHeeecdmTJlilxyySXy0EMPZfzmrVq1kpNOOklq164te+21l1x11VUybdq0mOds2bJF7r//ftl9991l5513lqZNm8qyZctinqOlGS677DKpWbOm8TqPPvqo7NixI+PlAwAAAJA9BM4BALlo2w5KeiDgwG3lypXlhBNOkJtvvlnuvvtuqVWrVsZv3q9fPyMoO3ToUOnZs6ds375dLr74Ytm4cWP0OQ8//LB07NhRfv75Z+P5ixcvNjJ+TUVFRUbQdtu2bTJ48GD5+uuv5auvvpJnn3024+UDAABAxUJt1XAoYDMAAHJIz8mxCYZAVicnU7vuuqs0b948+nOxBwXCu3XrFvOzBlw1Y1bLMZx99tmybt06+fzzz6Vdu3Zy/vnnG8/58ssvjZq7Guw99dRTjYnTJk+eLH/99ZfUq1dPjjvuOHnxxRfl8ccfl+eff16qVauW8XICAAAAAAAAQCgzbl9//XU58cQTjf9v3bpVHnvsMc8XRgO1arfddjO+awBXs3AvvPDC6HMOO+ww2W+//WTIkCHGz/r96KOPNoK2psaNG0thYaFMmjTJ82UEAAAA4A9KBAIAAKSRcbvvvvtG/1+9enW58sorxUuawat1c8844ww56qijjMeWLl1qZMzusssuMc/VIK3+znyONWhr/t78XTwaeNYvkwZ5AQAAAIQDlRIAAABcBG5NkUhEfvnlF+nTp48sX768XMmEDh06pLUgWut24sSJMnDgQPGbToqmk60BAAAAAIDwIhMfQEXmuFSCSbNib7nlFpkzZ47svPPOUrdu3ZivdDzwwAPSqVMnIxhszezde++9jUnH1q5dG/P8ZcuWGb8zn6M/239v/i6eFi1aGGUZzK8FCxaktdwAAAAAAABIjvg7kKWM22+//dbIqm3SpIlkSrN3H3zwQfntt9+kb9++0rBhw5jfn3DCCVK1alXp1auXNG3a1Hhs2rRpMn/+fDnttNOMn/X7yy+/bGT/6sRmqmfPnlKnTh054ogj4r6vlnrQLwAAAADhQWYdAABABoFbzao98MADxQtaHqFdu3byxx9/SO3ataM1afU9dtppJ+P7nXfeKc2bNzcmLNNgrAZ6NVh76qmnGs+9+OKLjQCtZgHrBGr6Gk8//bTx2gRnAQAAgNxTQJFbAMgL9McBWS6V8Pzzzxv1YTdv3pzhW4t89NFHRqmCc889V/bZZ5/oV/v27aPPeeedd+Tyyy83Mm7PPvtso/yBtY5u5cqVjTIL+l0DujfffLPceuut0rJly4yXDwAAAAAAAAByIuP2uuuukx9++MEoS3DAAQcYpQysRo8e7apUQio1atSQDz/80PhKZP/995cuXbo4fl8AAAAA4UNmFgDkFwZQAFkO3DZr1kxGjRplZLbWq1dPChjHBAAAAAAAfECHDoCKzHXgtnPnztK9e3c588wz/VkiAAAAAAAAAKjgXNe4bdCggTFJGAAAAAB4ySylxqg+APDWAU90lvd7zQh6MQD4Hbh966235LHHHpO5c+e6/VMAAAAAAABP5sSBO236zgrsvdmcQJZKJWht202bNkmjRo2kZs2a5SYnW716dZqLAgAAAAAAAD9EqBgM5H/g9p133mHoEgAAAPISzVwAQL4KIuuVUDGQpcBt79695ZxzzpHbbrstw7cEAAAAwomhnAAQLpRM8A5rEsjjGrf//ve/Zc8995Qbb7xR2rdvL4WFhf4uGQAAAIAKGVQg8xkA8gOncyBLgdvZs2dL37595YgjjjAmKKtXr55cdNFF8v7778v8+fMzXAwAAAAAAAAAgOvArTrmmGPk6aefluHDh8usWbOkadOm0rVrVzn00EPluOOOk2effVZGjhzp5iUBAACA0CDTEwDCgQoJPmCdAvkduLWqX7++3HvvvdKlSxdZuXKlEdCdO3euXHLJJfLKK694u5QAAAAAAABIW4TILZC/k5MlU6tWLbn22muNr6KiIlm9erUXLwsAAACgAiHDDgDyE0FjIEuB2/feey/u4wUFBVKjRg05+OCD5ayzzkpzcQAAAABUdAVMZwMAAOA+cPvOO+/IihUrZNOmTbLrrrsaj61Zs0Zq1qwpO++8syxfvlwOPPBA6dOnjzRo0MCPZQYAAAAAAEDIRzWQZwtkucat1q896aSTZMaMGbJq1Srja/r06XLKKafIu+++K/Pnz5e9995bHn74YX+WGAAAAECe4hYfAPzCGRaoABm3OgnZr7/+Ko0aNYo+dtBBB8mbb74pTZs2ldmzZ8vrr79u/B8AAAAA3CqgUgJQ4RFkzA+czoEsZ9wuWbJEduzYUe5xfWzp0qXG/+vXry/r16/PcNEAAAAAAAAAoGJyHbg977zz5J577pExY8ZEH9P/33fffXL++ecbP0+YMEEaNmzo7ZICAAAAAAAgLZEgitwCyG7g9vPPP5fddttNTjjhBKlevbrxdeKJJxqPffbZZ8ZzdJKyt956K7MlAwAAAFChEFMAYMd5wTuBrkq2I5CdGrc68VjPnj1l6tSpxqRk6tBDDzW+rFm5AAAAAJAOaiICAACkEbg1HXbYYcaXtfbtt99+K4899phXywYAAAAAAAAAFZLrwO0dd9wR9/F58+bJ8OHDCdwCAAAASAsjaQGYqMfqvSBWKVsRyHLgds2aNTE/FxUVyezZs2XKlCnSpk2bDBcHAAAAQEVXQK0EAAAA94Hb3377Le7jL7/8svz+++9yzz33eLFcAAAAAAAAyGH0wwGZqSQeueGGG6Rv375evRwAAACQddxgBouR0QDsIgy2B1CBeRa4HTdunBx//PFevRwAAACQdYQHwqGAEDoAAID7UgnNmzcv99iyZcvkjz/+kMsuuyzm92+//XbmSwgAAAAAAICcRccokKXA7ZgxY+I+ftJJJ8ny5cuNL1XAjAIAAAAAACANBPoAII3AbZ8+ffxZEgAAACBgpB4Ei1qWAJCb3u81Q848eA85fr9dg14UIK94VuMWAAAAALzA4D0AyC1v9Zwu//p4aNCLAVTMwO0ll1wiQ4emPgDXr18vr732mnz44YdeLBsAAAAAAKjAIiTiA6jAHJVK+Oc//ylNmzaVunXryhVXXCEnnnii1K9fX2rUqCFr1qyRyZMny8CBA6VLly7GBGVvvPGG/0sOAAAAeIz4QLAI0ABA7qLcDRBQ4PbOO++Um2++WX7++Wdp3769fPLJJ7Ju3broJGRHHHGENG7cWEaMGCGHH364D4sJAAAAAAAAABWH48nJqlevbgRv9Utp4Hbz5s2y++67S9WqVf1cRgAAACArKK0KAOFABj4AuAjc2mnZBP0CAAAAAAAAEiEQD/g4ORkAAAAAAAAAIHsI3AIAAAAIBTMhS+fRAADkFrJqAe8RuAWAHLJ20zbpNnFJ0IsBAAAAZAWxQAAVGYFbAMgh//tpnNz73eigFwMAAF8sL9wS9CIAADxE4B0IOHC7ffv2TF8CAODQyg1bg14EAAB8saxwi7zUeUrQiwEgJCKE/ADAeeD2p59+km3btkV//uCDD2T//feXGjVqyB577CEtW7b0axkBADYRCkgBAPLMmk1l9xoAgPxAxXIgM1WcPvGGG26QJUuWyF577SVffvmlPProo/LYY4/JKaecImPGjJFWrVpJ/fr15d///neGiwQAAACgIuNGHwByD6klQICBW2t2V9u2bY0MWw3eqiZNmshuu+0mbdq0IXALAFmgp2Qm3AYA7zCSIXgFhGsBIG9R+gLIQo3bgtIowezZs+Xiiy+O+Z3+PHPmTFdv3r9/f7niiiuMTF197d9//z3m97fddpvxuPXrkksuiXnO6tWr5aabbpI6derILrvsInfeeads2LDB1XIAAAAACA86JwGY6FgDUJE5zrhV3bp1k7p16xp1bTdt2hTzuy1btkQDu05t3LhRjj32WLnjjjvkmmuuifscDdRqaQZT9erVY36vQVst4dCzZ09jorTbb79d7r77bmnXrp2rZQGAXELzFQAAAACA/OYqcNusWbPo/3v37i2nnXZa9OehQ4dKo0aNXL35pZdeanwlo4HavffeO+7vpkyZYgSTR4wYISeeeKLx2Pvvv2+UbnjzzTeNTF6govln28FSo2pl+fbOU4JeFAAAAABIC4m2AOCiVEJxcXHM11NPPRXz+3r16hkTlHmtb9++xoRohx56qNx3332yatWq6O+GDBlilEcwg7bqwgsvlEqVKsmwYcMSvubWrVulsLAw5gvIFyPmrpEBM1YGvRjwGUPGAMBbnFaDR3kEAMi/exQur0jl68Fz5YAnOge9GPmRcZvM5ZdfLl7TMglaQqFhw4Yya9YsefLJJ40MXQ3YVq5cWZYuXWoEda2qVKliTJSmv0tEA8wvvPCC58sLAACA3Oa29BcAAADS12HMoqAXoWIEbv1w/fXXR/9/9NFHyzHHHGOUY9As3AsuuCDt123RooU0b948+rNm3DZo0CDj5QWAbKHnGgCQzwifA0B+4HwOZKlUQiqHH364kQXrpwMPPFD22GMPmTlzpvGz1r5dvnx5zHN27Nghq1evTlgX16ybW6dOnZgvAAAAgFI0ABAunJUBVGSeZdxq+YF169aJnxYuXGjUuN1nn32Mn3VytLVr18qoUaPkhBNOiE6apjV4TzmFiZkA5C/iCgCAfENWFgDkL+5fgIADt1dddZXrv9mwYUM0e1bNmTNHxo4da9So1S+tQ9u0aVMje1Zr3D722GNy0EEHSePGjaNZvloH96677pK2bdvK9u3b5YEHHjBKLNSvX9+rjwYAoUF7BwD8wfk1XKg1DAAA4GGphHSMHDlSjj/+eONLad1Z/f+zzz5rlF0YP368/OMf/5BDDjlE7rzzTiOrdsCAAUapA9P3338vhx12mFHztkmTJnLmmWfKJ598EuCnAgD/RQgxAIAvCBgCAAAgJzNuu3TpIh06dDCyYe+44w4jYGpas2aNkR2rpQqcOvfcc5PWEevevXvK19BladeuneP3BAAAABA+xMwBILeRWgIEmHGrwVHNfl26dKkMGTLEyIzVbFfTtm3bpF+/fj4sIgDAjhpRAIB8RhAXAPIDty1AljJu33jjDXn77bflv//9r/HzTz/9ZGTdbtmyxShjAAAAAOS6ZKPBkD0L12wOehEAhAWn5ZzBJRSZtsEoWZVB4HbGjBlyxRVXRH++7rrrZM899zSycHVSsKuvvtrpSwEAAAChQsA2XH4auUCebHJ40IsBAMgQYTggS4HbOnXqyLJly6Rhw4bRx8477zzp1KmTXH755bJw4cIMFwUAAAAIFpkeQWLdA0Au4xIKBFjj9uSTT5auXbuWe/ycc86Rjh07SuvWrb1eNgBAAiSGAQAAIJ/R3s2vbcbmBHwO3D788MNSo0aNuL8799xzjeDtrbfemuZiAAAAAEAJkrYAAKhY6KzJsFSCZtbqVyJaNkG/AAD+i9BnDQAAAABAXnOccQsAAADkK7rDwlUbkVrDAEycnwFUZI4zbitXruzoeUVFRZksDwDAAYaRAADyGWFbAMgP3LYAWQrcRiIR2X///aVZs2Zy/PHHZ/i2AJB/tmwvksLN22WvOvHrgXuBgC0AoCJYtXFb0IsAAACyiFvdDAO3w4cPl88//1zeffddadiwodxxxx1y0003ya677ur0JQAgr9333SjpM22FzH31Mt/fi4saAAAA8hlzOuQHRlAAWapxe+KJJ8pHH30kS5YskebNm8tvv/0m++67r1x//fXSs2fPDBcDAHLfoJmrgl4EAAByFjf3AAAAGU5OVqNGDbn55pulV69eMnHiRFm+fLlccsklsnr1arcvBQD5JYt3nFq+BgDgHU6rAAD4h/sXwOdSCVYLFy6Ur776yvjatGmTPProo1KnTp00FwGAH2at2CCN9tw56MWoUCqRKgQAAAB4ioAfUJGOdW6q08643bZtm7Rv314uvvhiOfjgg2X06NHSunVrWbBggbz66qtSpUpaMWAAPuk9ZXnQi1DhVCrI3kWG5isAAAAAAPnNcbR1n332kdq1a0uzZs2kTZs2stdeexmPb9y4MeZ5ZN4C4UAx//wO3AIAAAAAgPzmOHC7Zs0a4+vFF1+Ul156KW5Kc0FBgRQVFXm9jACQE7IZt2XEGAB4iw5PAAgX2rv5gc0IZClw26dPnwzfCgDyGxm3AACkT5NAAAAAkEbg9pxzznH6VACokLJ6v0nXNQAgzxC2BdIzbel6Y9TAYXtTthDhw7kdTnGLm+HkZMXFxfLaa6/JGWecISeddJI88cQTsnnzZqd/DgB5j4xbAACA8Jm7cqMc8ERnWbw2P+9fG7fuL5e0HiD5imAOgIrMceD25ZdflieffFJ23nln+dvf/ibvvvuu3H///f4uHQDkkErZrHFLExYAAMCRvtOWG9+HzVkV9KIAFRY1iwGfA7fffPONtGnTRrp37y6///67dOzYUb7//nsjExdA+HBhzM/afARsAcAfXDeB/EX9ZABA3gdu58+fL02aNIn+fOGFFxoXwMWLF/u1bACQU8xbgkgW7v4JMAAA8g2xNfiN9lNuYXMBFQvn6AwDtzt27JAaNWrEPFa1alXZvn2705cAgApR47aomCsOAABAWNApAADIVVWcPlEzyG677TapXr169LEtW7bIvffeK7Vq1Yo+1qFDB++XEgByqMZtNuK2hIYBAADcIZsLyD4OOyBLgdtmzZqVe+zmm2/O8O0B+IULZHD104q5KwAAAFk2afE6mbdqkzQ5ep+gFyW85awCXg6kh6Y1gIrMceD2yy+/9HdJAHiKBk5ww/Cyse6zUUcXAHLVwjWb5MzX+shfzc+Wg/aqHfTiwKGCaHgN6bjsvYHG97mvXhb0ooQPtRIAIPSYiDvDGrcAcgsnPQBARTVq3hrj+8AZK4NeFAAAKjS6TYDMELgF8hQJmfkdNGfzAgAAuCyVQAM5p7C9AIDALYCQ+GrQHOkweqHkA9qYAIBcVVQcIViCvMWeDQDINQRuAYTC8x0nS/Ofxkle1LjNwntxTw0A8EOjJ7sEdj2mzBP8QolboLxnfp8oT/w6PujFAKK4x42PwC2Qp8iWAQDAOS6bZX4bsyjoRQD8wXEORH07dJ78OGJB0IsBIAUCt0Ce4gY0v4PmZCUBAAA4U1Ba5Zb2U25iuwH5jUERyRG4BfIUzZv8XPcE5AEAANyhVAIQHG5fkAr7SHIEboE8RYAvz9c92xcAkGdouwCw4pQAAARugbzFkKIAseqBrNqyvUjWbd4uYfXZgNmyYv3WoBcDKXDdBPIfnQNA9gWV8D51aaHMX7UpoHcHvEPgFshTNEzzG5sXKPOPDwbKsS/0kDBav2W7vNR5irTowKzNuYIh1UD+MQ9r2k9Axbk/vaT1ADn7jT7ZfVPABwRuAcAjTHwBBGP6sg0SVsWlp4OtO4qDXhQ4RMcnkH/okAEA5CoCt0Ce4r4zv2/6CSwAAHJdkdmzUYpLG/xG+yk/t5uWSyq2nU8AIF8QuAXyVISWaWBY8wAAJLdg9SZp9GQXGTBjRdCLggo0Kgr5ScslvdFjWtCLASBDhDDiI3AL5ClOevmNcgxAjig9VAsYpxt6XDeza+6qjcb3EXNWB70oqAhKT8G0n/L3vNx3Gp1AAPITgVsA5TJ1lxduCXoxcpIZlyHbGQByFzH27CADEoCXaH+HF1sGqdAiCHHgtn///nLFFVdI/fr1jUyU33//vdzJ99lnn5V99tlHdtppJ7nwwgtlxowZMc9ZvXq13HTTTVKnTh3ZZZdd5M4775QNG8I7SQmQLR/0mSkP/jDG9d99MWiunPxKL6NWFMLXOKFNCgDIVwRe4Dd2MQBwZ/XGbbJ+i7+xAU7NIQ7cbty4UY499lj58MMP4/7+9ddfl/fee0/atm0rw4YNk1q1aknjxo1ly5aybEAN2k6aNEl69uwpnTp1MoLBd999dxY/BRBeazdtc/03YxesNb4XErhNu6eQyckAmMxhuWQS5A7Or9nF6kZW22gBLwdQEdEGym1/f7GnnPtG36y8F+Vs4qsiAbr00kuNr0Q97q1bt5ann35arrzySuOxb775RurVq2dk5l5//fUyZcoU6datm4wYMUJOPPFE4znvv/++NGnSRN58800jkxeoyNK5+axepaQ/Z8v2Iu8XCAAAgJIUyDLqjPtv4qJ1slutalJ/l50Cuaeh0y38whKU00DkgXvUkl/uOz3t1ygqjhgB6UqVKsa5ZdVG9wlhqAA1bufMmSNLly41yiOY6tatK6eccooMGTLE+Fm/a3kEM2ir9PmVKlUyMnQBuGdeerYVFQe8JLkrG42ScDR7ACB/cF4FKgCie765/P2Bcs4bfYJeDMDR0P+R89Zk9BqNnuwi93w3yrNlAnIycKtBW6UZtlb6s/k7/b7XXnvF/L5KlSqy2267RZ8Tz9atW6WwsDDmC0CJ7aUBW9q1GWDdAShlnktJ9sodbKvssrY3uHx6Y9HazUEvQuhwWJfXY9JSOeCJzrJh6w7PXnN7kcdHcST3sjkRbjs8TE7qOXmZZ68F5GTg1k+tWrUysnfNrwYNGgS9SIAv0mnAbC+m0ZPpMDw/1+DkJSUdTUzgAgAA4vlt9MKgFyG0HTL51Hpat2m7LFi9Ke2/7zxhifF9TZ4MgaZpDCe27mBUaZhxHOdY4Hbvvfc2vi9bFtuLoT+bv9Pvy5cvj/n9jh07ZPXq1dHnxNOiRQtZt25d9GvBggW+fAaEQ++py2TwrJVBL0bOMAOCnDTTx7oDAABBoR1SMVzbdrCc9TqlCQA3KjGsBjkotIHbhg0bGsHXXr16RR/TkgZau/a0004zftbva9eulVGjymqL9O7dW4qLi41auIlUr15d6tSpE/OF/HXHVyPlxk+peewUjf30zVm5MWvvxXZCPg1Ze+q3CbJi/VbJZ9wm5A7Or9k9Jqyjg1j38Fs+7WMzlm8IehFCJY82rS+CHK0Xpm1D3Ba5KNDA7YYNG2Ts2LHGlzkhmf5//vz5xpDjhx56SF566SX5888/ZcKECXLrrbdK/fr15aqrrjKef/jhh8sll1wid911lwwfPlwGDRokDzzwgFx//fXG84BsXADbDZsvW3cUSRhlcn2mTlT6WHf+6T99hczkRiWvTF26Xr4fNl9a/zU96EVBBUcJmiyz3DwvL9wi//l+lGxjCKsn2JOTlErIo+O8SgWZzd6pfNq2CJ/l67cYNaH7Tosd8e23iYvWydeD52b1PRE+VYJ885EjR8p5550X/bl58+bG92bNmslXX30ljz32mGzcuFHuvvtuI7P2zDPPlG7dukmNGjWif/P9998bwdoLLrhAKlWqJE2bNpX33nsvkM+Dikdno3zytwmydvM2+c+5B0k+oe2TPtadf279Yrjxfe6rlwW9KIAjnA5yT6psHO08mrVigzQ+MnFZLrjz6YDZ0mXCUjnr4D2DXpS8QDukvALGPeQtArbhlw9H36I1JZM+9p++Us49dK+sve/l7w80vjc7/YCsvSfCJ9DA7bnnnpv0RKtZty1btjS+Etltt92kXbt2Pi0h/KD1Zk9puLtUzoNeYjMzZPO2fMy4RbpYdwASTV6I3HfxO/1E5/GkA8m7torZXtnBBKmh8euohVKzWmW59Oh9gl4UVGBuRrFx9gg/YuxIhV0kx2rcIj/NXL7eqDf72YDZQS8KUqD3OtzYPAAQDGKL/mVAtuw4KZBlQXn/+3mc3Pf9aMm7UglBLwj8w8ZFFlASD0EgcAvfa8HMtUzYVLhlh/F9ybotAS4VkiEgmDmC3gBMnA9yB1sq+HW/vYitAH9xSs4/bFNnWE8IM8alJUfgFr4689U+cu6bfcs9ni8jRsN+AcykRzDkH00q+n5Bby+QW/Lksgd4Tq9mxWFvUOUY2gjwy+wVG6SIYQfIYVxuwonNkhyBW/hqW1Hs7MCcKP23d52yyfsywbYCAAB+sXbi0+ZAtrCrpW/9lu1y/lv95L1eM0I5eoVti2xgokN/MVItPgK3QAbCmDmsE0l4g5NmqDNu2TxATuBQBVJfz8i4RbYmiCQokL4t20sScqYsKZQwYYuGH9sIyAyBWyADYW/7ZbJ8ExeFq1EGWM1bVVY7G8gFYezog3sbtpbU6kfmCnKoPZVrWJ/lcQrOzRIcbvZlgvLJsXZKUEoGuYjALZBnIh5dzJ77k5mdw9wgqIhNjvYj5kf/v2bT9kCXBXCK+8j82la/jFyQjUWpcLiRBpApziLhlU8dJ1yvEAQCt6gwtWEOeKKzfD14blbf88AWnbNeB8ra28xlJX8DNbd9OVwqGmtnwmtdpwa6LIBTNPBzD9nR2T9G0r1uaptn87Yirxcp53HWCfa4Li6OyMI1myTs3B53Ya3tSaZt7ojkQbkV+CuX9xE/EbhFhToUvxg0J6vvp5Ouvt87+wX8kf97+bxV4b8h8Ku2mhoye1WgywK4R4M/VxADyL50J6n/avBcOfzZbrJlO8FbOAs6ZuP4btt/lpz5Wh8p3BK+0UGvd5sqbfvNknzEuRtOsJ8gFxG4RSDMDqul67bIxEXrJFc56XhL92YkXVyLAMBbGhQqSvdkzkk5r5BxE66MuX7TVxjfybqF6jFpqVz14aCgF0PGLyi5t9kSwv2yTd9Z8mrpqKVcOp05OUMwwiU5MpOB3EXgFlllv16c/1Zfufz9gZ6+x7pN26VV1ylxb7K9bp84uf4FOVPy8DmrZeTc1YG9f0VFw8hb//fjGBnBflzhDJ29Sk555S/ZUVSWaR2Uw57pJvd9NyroxQDyMwgeSb+tZLbruOraVNB2iJZUGrtgbdLnENxLX1jXXTiXCoAbOdSPFAgCtwj0wNzkQ0/0e71nyMf9ZsvYBWvKv28AXcvZbjvb3+/atkOyuwCo0A1IDVrPXL7B09f8Y+xiefTncZ6+JsLvy0FzZFnhVlm1cZuEQY/Jy4zvnw+cYwSVnarI54Oc42Bj5VKGWthZ12W6bSUyoOEUu0puxvcjefZ5EDx2k3BvF47j+AjcIu8kO9jTHuqahm07gs8SQzAq8gXnu2Hz5cK3+8m8VRs9fd25FbCmb0VXu0ZV4/v6kNUIfLHTZLn+k6FJn7O8cIsxIebclRuj5wOCBrmDbZV96TbPohm3FfnCG0dFXRvJdoOyfSVbS4NsY9uGVz5sGq4zCBKBW2RVJOAbrvmr/Qn+xHvLjVt3SBDCOoypYqm422DmsvXG97WbvAm20UhC0JUS0tkHR8wtGfHx15SSLF3kD2K63otk0HYx23zZnk8AuYcOmfxdJzQVkU0VbX8LWwJFRUXgFlmVzYzXoE+qFeycnhGCc0iEm/GKyzwtBFknXG0vcv/+1iAUnWm5h0tSdlhjRlvTHqVU8iocZ3CKPSV9BXRd5TT2/czuOytqaZ5s7TcVc+06R+AWBq1JOXuFt3Up4zFvwP0871X0gz7TG84gbljzLTjHTb93dhRTcqSiMgMxQQdurbPVpzOSgvNBnqmgN25+6zx+SWabg+MsBuedJEF+1k2U23URZAeJk2UlESS88qmsTUVrBlTK0gfO/T3DXwRuYdCalOe/1c/398nXGEy8E00+XJiyJd/WVX59moqTpY+QiYTjujF1aWH0///6JP3JHitYOz8nkbWZe20A87jiUuHfvtx7KiVfwsrcyhsyLM+W681wzt1AHt2nczjHReAWWVWWcevfLaz50kEf80G9fy42vnJwkfNuG3jF64++g7vxCsvc8kFn3G7cVnJDXLtGFZm4qCyI6xR7cO5J1kQhAO8dL5qCZTVuOdL8snjtlqAXASnoRJiZ3TNF8iJoFN5PgTBhP3GH9RUOBG6RVdluWBdnKejDjRwQy6u+maIk9UWzdXwj2OtF0AEZc4hYnRpV0/r7fBtRUJFpVpu1dAa8kckhYtbcDPo8kc9yac0uLdyS00FKtwo8anO5PXzCOlSc04Dz9TMnzWA//KNlK4Oa3DxdOiLjgCc6y7a069TDKQK3yKpsNqz1rWjI54Z820z5dFMQtGQZt+s2M8tpPtmyvTiUk5N5tRxhudHVDo8t28uCj9uLiuWzAbNlRxGN7nd7zUj6+6Oe6y4vd5mSteVBatFgHJdd/+TJyi2rsSl5I5JhDcp8WRf58jmy6bw3+wa9CDnLr/1Ny1be8+0oyaXP+8PwBcb3tZu2Zfw+IWkmhxaBW2RVNmoVWsswFAV4Jc/VRkQQSYz5FujsPXV50IuQN5IFy9KfhRxh9OvohQlKJUhOZ/6G7VrQosMEOeyZbtFAbafxi+WlzlOk4/jFUtF9OWhu0IuABFnrLTtOlgWrN5X7nRmwqgj10B//Zbz8PLLkJjnXzjthkM8zwmeccSv5Id/uJ/JJmLaMm/Ojdm6bpUiycQYZt2CthE6S9VWjauW4yReZvA3HcXwEbpFV2Qyk9pi0VBau2ZyV9wrT6SXTYbn5MqxXh2385/tgei27T2Iij2xk3DIsJ/9Yzz9moDTogIz5/ssKtzr+G/NjhDFQ0L408PP1kHnG983bSo6j7UnKklRUGiT7dgjBXP+UHB+p9ry1m7bLF4PmyDN/TEz0EoFn5mfr2H30l/GOnuvl2sj/NZu7Cjya9d1t2z+Iw81JMKcCnAaQZa26TJVz3+ybvbZo+JqNSY+96lVKwolbdlBGym8EbpFVZsPAq3OSDvd8qdPkmGGf5mt/NnCOXNc2/RnAc1Wml5Ug2jzWhpaXgeMuE5ZKEEJ4zS2nw+iFMmTWKgm7ZDVutxXRSMg31oaxeS4IOiDjVVvdrMUZFtq5qTaVTr7WpwKMFHB7fdEg2TN/TPJteeDMkNkl16o1mxKXx/HjnvqfbQfLLZ8Pk4ou34Jh2UhQyFbGWFmphPT+fnPp/dNPI2NHvAD5wItjfeKidcb3HT4PG94a4sBnstVonnqCbqtXBARuK6iP+82SS1r3z/r7el1C789xi40A7e9jFsX9/aqNmddbcSJeeylX0/yDzrjNh9GO6Tags6n5T+Pkhk+HStglayhRKiG/M6zNIK7bU9LIuauNjHszIBn8KAYJJXP9ajaj6joxmI6ubNHhjg1bdJFfRhGgyDWdxy8xvscrlVBWt9T7A23E3DUyYMZKz18XwciBplna0h3dUVg6V4Dbjjsnb6dJNdbEmnS5ObJDerkNjSDvTaPn6iy+Z8LLgouFqFx6U6e3I34ue/P24ySskn3uaFkzD27J8vkc7QUCtxVUq65TZerS9Vl/32hvjEdHpjmrfKZDhPwweXFhIO+b6b1L0IFTv3s0syGM+2O2eH3vnGxoEoHb/AyumcxN73Z4Wrth843vS9clnl3cjSBrpfvJ/Fwf9JkpFYEZoG4/omT/SGbB6uyUWaronEwspgFZ88Z587aihAGrfD1OwyDoDv1clK0RFua7pPtuZWV9xHOnvNLLmNAxm9hVYeXF7mBef/y+Px00q6STMNfuIL2cSDha45bjOC4Ct8gqr9Poi+M1OEJyxvsiRyc5CbpUQtD1LHM1cKsBrzDVfPXqpoUatxXLDktpDLNjTrPtzIm0gspOz0RYRl+s27xdvh1aUtfWq4kkcom5HcJWsqIi0IDrFwPnpBUA1NOAeeMcr4ZeWcZtxouZV7xcH5EK1FGQyv3fj5bR89dI2IQxYUCvOcnacF7i+E9Pz8n5PSeHF51OlUqvP37fn4bvCHa2HunYyx4CtwgkcOvVjZP5ekE1WJLdkAc1XD7TIEEQJ2DrMufFBDkBbPvz3uyb1cyGwTNXGl9+S9ZQ0rpTr3eb6vsyIHu2WzIazAw6nZDn8Ge7yfiFa7N+/On5MNMOAj8zmtx4sdNkeeb3ssmdvBjCmlMi2b8Wzl+1ySjbMWFhSY28imbS4nXS+q/p0qbvTGnZabKMSWO2bL0GmIHbeJuk15SSwEOyzaXLodth9ooNEkbrNm2Xbh6XKglLh1EoA7cZvEbnCUvkyQ4TJGxS3QfpCJTn/5wU7RBNez8J/W4V+gUMlbu+GSkVkZv9vko049Z58DKd9pX5CmGc1NZRqQQCuL4jcIusMAtu25Omzjho94xukKKTnVnOcWHJpgnHUuSgPDjvW7d932nLZc7Kjb6/58I1m2VbFrMSb/xsmPEVZOD2pc5TpE3fWb4vAwLKuLVseu3Que5jZ5NNdhi9yLPG76Y4Q7NzlX2YuZ6XKlKmhLk/Ob1h8yK5RgOGamjp5FoVTbMvhkvrv2ZEb2Ktx3dszcPEK1tvBisnOZY3lu7XyV5j+JzVxvcJpZPMhM2jv4yTe78bFdoRR/lymsjm/UG2JydL5aXOk+WrwXNl0drYMjBe7nJTlhQGvh/my77ql3xeP9qeeb/XDFljmd8m4mGniJ6fCxyu13wr5ZaqnJEK6eUrrxC4RVZ0GrekXP1CdeWxf8voQmKeJILKuDUbgTPjZHEElgWccY3bADJuLW+Zbz12t305Qi54q69UNF7t/uaMx6gYrIEde1DR+qMOVTWDHDoZWTyZnEv0BnTj1h2yfP1WyVRozmhxjsmlhd7UAc4FyQIpGljUofxu95+KFPjORKJOFCedK0bGbeWCvJ712pxIN6z7UziXKn0hXc2+KsvmSzAS0mmjLcHTCrdsl0vfHSBBCWOGuR7P//56pDE6DP6fR+at2iRv9Zwur3SZUvaeEXejY+wdG6pyabRMM24jPn7YMCd8JfvkXmbcBjGBXS4hcIusMA/ARWtKTohm+8D8nu6kEuaNu1l/xvqa2dRlwtK4tZ1yURANWutbhvXGyg37R6AXMn3WnvNEwnqzi8xKJdjPBeaPOrP8NW0Gy8f9S7Ktu0+KP8TYPiTUDb0BffCHMbLawf6XzsiQsAjLCJWsZtzG2S0+HzjHGMrvPnDr2eLlJS+y3rV9mGnpKXM5KtT28rLGbRZXnA7p1wCKmSWdq0G+bE9OluozlV2HYpfLqzmjt4akZnqYDnHNuvxryjIj2xnZ2/bWmEKi48J+Sus3Y4Xxfcis8qNjqlQqCZcVFUUSngv7TFse83Mm97J+tBW15NdvYxamfy6PpNe28vBtQOAWqejwPu1FzZQ928H8btYtS+cEN2remuiNVlD1ZJPJ5rB1L096QQdO8yHIGcae/2x/9unL1nvyekvWpc4IrFA34xUo49Y+bNjct9Zv2WF8X7B6c9LtPy7DuqKTFxemX8an9HuYLk0FFXyCP3NbRhyuB0fDcqXi0LI/l78/IKMZ7xMdT0mHYRZnPkTai0mpkB2zSkew6f7muTCdkD1i7tKp9m3z9/b7Ja8OCT87Jt0ct350MuiIjAGlgT3kFi92BzM5bEdxccLj5c6vY2sFh+1S83G/WfJw+3EycVF65UycfB6SaPxH4BZJXf/JUHn4x7EZv05pZ1W5XhmznIAlycoxnWk8+vohnJzMXsstVwQROLWe7PPhxO/FR1i4ZlNOrovtO0qW+WnLJEiZWLZ+S+g7G7KtybsD5L7vRkk+spbTsV8X3G7mR34el/H53Yuak2HeO+MNC8xX5v4T77xqTj5i5WTbV6Rzz8udp7i+6Ys2zRJMCuWk5abZU45uGpM8K5pxK+Hm5fKF/bPm8lDZMB72kTQnc3Y9IiQSzpi4n9ukVZcpcsvnw9MegRPG/SUIZR14klPMSj3aJnC67Oncv5nXKT+OJXMU8JbSOYe8RI3b7CFwC0eTHnk1ZMhsWJtDWM1erHRufmInJIv//yAFdUOX9iiI6F1t8pmP/S+VIDkv04+wYv1WOfO1PvLloLlS0a3akLqhnA/7jBuTlxRKV49nIA+SdfZd66y99nNoohl9/dr8+vZe7FtlQ1GDuzrNW7VRFsS5lv8yaqFUFMmujfFuZtyWShgzf03CkQYVdRSGfZ9Pp32ibZOMM24trwX3srraCvLj82TrmI+3b8fbz7tPWhbz/PI1bn1cyBxnjvzaEdBISngvklbGbcTxX6bTdvTz+pTp8W1fNB0Z0WH0wpjfeVnjFvERuEVWGh/RE0b04C75XimDGrcxwdoQtjjMySZygU7KcuRz3Y3hoom2t04AdGzLHsakQH5a70FpjqBlevEyy5OEZUIDnY3eackUrzssdLh6tt8T2WWtLWq9MXK6Xb3e/NbGc77sW+e80VfGLVhb7vENW3P/fOuUuS23xKnF+GGfknrJVhu2lpTkSMZ6vby6zWC5+J3+GS9nPrG3/dKRuMBFvOclX46wH815crrxJtHDx3URycdSCUk6Oa2jB1ZsiJ1w0/V6Dt/tVui3bQhvUSsMt/v32k3b4tYqVm5GYGUShM0kprFqw9akHQzpl7iNxLSNLnq7nzT/aVzM77yscUsna3wEbpEV//fjWGOGbvMwjNa4NYevpdGJaT2vmbVy7Y8HSbMmg+H+ZPfTyAWyaVuRkXWU6Lo0rTSTaObykvpjXrKen72YDCho+ZYBeum7/eWY53vElHFIfdH1Lrs0Fa7vuU1nAjZtt5SYGT2/fKAxG8zGecSzwG14d1CtFV9RRFycU9RprXqnfM6VHwySWz4fluGS5S/7/AbpJAJ4MTS1LBgY3mPRa15+1mxmjJcF2X14zyxu/iBHWNhLtVn3hcveG2j7ncvlDfAQSvbWYT6yK9BpJ3Tcnkde6jwlYWlGDeBat6VOlOvlPDdeJKGd8NJf8uyfkzx7bU3csY96Oeq57jH3ucnKUMFbBG6RNcPmrLLUQYnEnEjSybi11mmKLZUQkshtDpm6dH100p9EgQpfz8eW117jUzmGrPJoZQ2YuVLCYK4lsJaqUyKIDMV0zh8Ip3T2H69v8K3ZSp7UuDVfIoSXppUOSpHkimWFW5KO2PDjpkKvnQNmhOM8HUZmGy96GEUSTRoWSXp8O9l0yQ5Vs28/7JeKsJbUyOZ6i54m47zn8sIt3kx8muEHCuN2sn4knUQp5nfJ/s6j9w/LyMewH+MIx/6QTnugZFRqmeVJ7oXijezJloFJ2iRuPrcGbQ9/tpt8OmB2wvOEtpHNCYPzLWkpjAjcIiWvAqH6OvY6KGambFrBnoLyk58hs0aLkRWd4u/9bpq17Fi+pzDXZHrxKgg8azt9ZaWSs3cFz5fh7EhcxzYZ6+Y366d7sT/p8aejEZAbTnmll1z67oCEv5+QYemZpMHFBL8LSSwjcGUZt7Ecnbojzp6nN5m5WiohVyft8UOyieTOer0P5UgSiiTJuE3yVy5r3CZ6KT9PdW7ak2HK+gtiUY54tlvSUSB+L1OY1r8rDpa73/QVtqcl/ht754kTtrk80xbvWE6natHm0vknXukyVdqPiN8Wfu7PiTKwNMlIj9PxC9fKAU90Tnt0Lk2m5Ah3ISXPAjAFZScMe43bdG62rQHlsPT0hkEmW0svSJGA97F0hpcgfb2mLDMustYJolJJdryZjbZstt3SKbWCcCpy1dgtv5Nlcv6Ys3KjvNVjWkzwuMuEzCeBy9HbmLybTHVMhqU3kjVTNiYIGubqPaxXzEtFJh0qTv/yofZjPKubeuWHg+SH4fMdvnP+C6RdGGdjmbUm3fy91mM0S3Bls0M521m51tW1vVzGbeJliY6AlNyWswFDj2nZu2yNAmnYorO832tG0udMXVpo3GPMX7XJ9yMi3iTbXr6n3ifFTgKY+LnLC7eGa58syOxPRsxZHfc55qSHSi/zI+eucVWSKpEQrblQIXALXx24Z634WRel382ZGtOanMxyRlmZg5mJ2fTzyAWOJ3lKdaHx42RqfctMgvDNfxorD/6Q+OYtW8I4jC6Rn0eWzApauNl5iYpk+4h5f57NNUDGbW6zHvJO4q4aYFXLCsuf9910QNh90HumvN97ZtK6ZVZOGuX0KYZH1crlN8bclRtl3qqS/SmVZGUztrsIKFXIwG20Bl7850VSdihHMiv7UeDuWqET+bXoMEFymZeXxayWSvC4rMW5b/SRv7/YM+axdF86VIEYm0iSc1Wu7gvpCNPi5dK9QLr7wgd9Zjoatj92YTBzFiQ6ZtPZMtuLYkslJHuNW78YLme+1ieNd8m83ehVs9OaxJAoKSJmxFsk4qj8UTL5fcRkjsAtstZA0WM5Wiqh9Pg3JydLJxnDMh+ZPPrL+Ap7ozx45kppbBk6Fm+76fpp8l7iIaTWxl6izZ6t9Wrdrm51GL1IOo5bLEFzO0ImzDcDbtk7Z7L5nhWVZhSd+VpvXyYOzLYek5Ym3Xd0GFavqctiMgCsz7dnPk5Y6Hx4/LrN22I6M1JJdt1ympWB7LFur0GlQ/vOfbOvnPNGX4d/n3hDJrqpqWjtkUSZrvbJyT4bMFsebj/W0U1au+Hz5SeHx2SqORHCfih6GmBL8++6TVwi17UdYnut7K25Ao/XZbyAfj6ek2Nq3LoolRA9L7o8WdnL6IVmlYZmQfJzP3PLmoxTlrUZ8fe9LLuml++kgdvo69om7Ipn0drEI4ByYZ8stMwZkGgi29i2biTnM/fDjsAtfL3gWA/gNZvKGk9mI95sTGdaKqEie637NJm2bL0ngTLtXQticrJIgknnsumhH8fI8ARDQdxyu6rs6zbspT+Sl0qI/Z7KM79PlCd/yyyzqaIXxJ+ypNAYIv5jjg7tte5NP49aKJ1KZ/CN5x8fDJK1pRMYmsNBI7Y63YkanqksWF3SyP5q8FxPOwzMoEe4j+r8t2Rd2U3Uy3Fmjs4k41YnLYmnot+4J8q41Zm7fxuzKOHNoNV7KYbiOlqOFBtEb7A7jM4sOJwvNNN4+Fxv2kKZ8PPQSTcQnc7xHMQpwBpgKlmGxEvR+Mh6xvcrjtnH1XvYXzMbnfVJa/X6/u75we9OmGSTgSUqKRQKDu671m3enqV74QLP79H6Tl1R+h7OP4CT2Iy1Haz/jdYo54D0BYFbZM3D7cdFL+zRGreV0p+5O9E5tqLdHBd4OERd60sGMTmZtcEXVOD297GLjVILXnDbgM329S2mh9TFu2t299cpglrmRdxpUOvbofOk3bDMAo65nLH8yM/j5Lcx3gQMQh7vd9zIXJuibEfvqctjrhvWzW8P3K7c4LyMjnaA5VOmt9aV+2LgHMk1P41YIKe88pfnr2utVzx7pfvs9GQlndzU3qyQ7ZMEbbwXO002vvt9KCWbnGztpm1yxqu9pflP4zyZ3DDXh1ZH8qRUgp+fw8nrZDu5xLrvuCmVULNaFeN7nRpVXb9P7OPh4G/A39/n56NMRlH6G0h2+ESLQTNXlXXEF+iE695uYbcTBSb6ey/at2agOuV7Wv6vp52ya21664aJOpMjcIuEIj7ckJuvGa1xGy2VkE7GLawXxXSC3/GK2nOuzJzuzjrccNjsVQ6fn921nu6ucuNnw+S5PyfF/V3XCUvk0/6zy8qhZPEj5XLG7S+jFhqdWpkw13nYM7Wt2UDnvdlXxi5IUPPM4fGw3TYcNF4ALZ3MSqecHrbWbEytwZvNkhZuJ1jqO60kKB6k17tPjVu/2Etbthd7Ogni1h1FcW9Sc+SQ9I15TiorlZC5N7pPlSveHyiPW0pkpV6Oku/xArOfDpjtuNMoV+q6ZiLeMlx2tLtszMyYZS1CsDJsIrkyOZm9VIKTvw/h+o6KVIyO/PxOLgpou0Q8zmSPuAtsZnNRo/cCXiyMw/aL9Zpq1Li1LYtbHL05HLh9/vnnjUaf9euwww6L/n7Lli1y//33y+677y4777yzNG3aVJYtK5vdDpnx4+JnvqSZuVK5UgaB24p+R5ToxiiDzfZq16mBNHry7USt2+Le70bLvz4ZGsrPbz3evMoKue/70fJylymBNM52uC0qnGdybVboVRu2GROMfdQ3/qQWEZeTjGnWdqLArT0D10tO9/VHfx4fPVf/76dxcuHb/SRbdAkXrnE22Zq67csRUlE4nZQsVR1bddl7A43vVSrHNqvzMYbgpo1QbkSQByvkwz6zZMKiddJ+5ALPa9zuKCr2pBM8DNK+cY7zh2ZbPacybpO8QDaPy2y9l/V9ymfcRjxbTscZjFkW9Pvn4qicbIh3rz5lSUn25+9jFvleB9ZJh4Se952WE7S+2oI1/ix7umdbR3ubi13SyWqx7uLGecZWKqHP1OVGOTdUgMCtOvLII2XJkiXRr4EDSxrH6uGHH5aOHTvKzz//LP369ZPFixfLNddcE+jy5qNMrjv2k495youWSig9wLWm1qZt7m6wE8ZtUwR0tUGjdTVXrPc3oydbzE/bdWLiupDJPPZLWbaf09nU/dzHzjho90CWwb4c2WysZbttN2a+s9ldtTdZh1vbJTvCgrjn3RzmullZpI3KT/rPctwIzSadDMoMINpvzFdvLD95zHiHk4rZb0i32gK+m1wEgN1yuq9bh6iNX5TdmZU1uzfdmY2D402QqMfk5B35V7cZ7Or1Rjio+1k16HGhYWObYToScLZfqmutnkPdtkW9FAlDHVcJVjaOoHQ/Y5izOa3b26z/XvY775WfnCwc6yYcS1Hi9W5Tg16E0LAeOg/+MMb4/lD7sXL7l8Oz9r5mDEATB6z76/u94ycRpJrAu3oVb8NomZ5eoveeSU6iEY+TYqwlpIy4re33t381Qi59N/Xk6MiTwG2VKlVk7733jn7tsccexuPr1q2Tzz//XN5++205//zz5YQTTpAvv/xSBg8eLEOHOstyQ/YvfuaNbtnkZCU/j5i7Rv4cu9jVa6WbKTh1aaGRofV2z+kSplnhh8xKPqxebyb0eXZmEMQcbu22YWmdrVkzihJOTibZUX+XnaSiyXaDd9pSZ72f2riJZ9icsn213/QVtuzIxJ9lzcZtvvSuxzsuKhLzkNXshVe6TJXOE9LrxPHTTZ8NMwKI2hGQati0/rqawwbxH7brxhZbxq2f99lk0zij27xFB/fD2jPVOckkd2ZtUzc0qJeKPeM2TEPgk1m3abvjjls3n6XAQX3gbIhOmOLgxnzx2i2SbebqeavHtArfEenVaDpfa9y6em72OyV22EslJIm/pOpUKfd8CbcwnWu7ltZVD3qRtCzV8vVbsrJ+7C9vHs5aQzye9Vuy235/v/cMo1SX9X3nu0hash7P3SeV1c33dPRc2jVunT1v5vL1cuJLf8n6FJP3xitHlqwUWGyN2/RESy0EftSEU+gDtzNmzJD69evLgQceKDfddJPMn19Sr23UqFGyfft2ufDCC6PP1TIK++23nwwZMiTpa27dulUKCwtjvlCeFyd3+8mnbMb5SLkMlic6THA1MUS6k5OZQ76sAU7NUBszP/Xsxn7RWm03fJq8w+Hkl3vJUc9197VkhJ6AU20CP0pUWE/QQU8O4oVIhsea31kdVW3BhWSB1ng0OGhq9sVw+cDSW20fojd41srosXXOG32MSWCcStWoMNlvdHX9dZmwJJSZp34wjx9zJulEM9yHhXlTmWw3d7rtlhVuSZpx66dkN8PxFGRxwpqwZYb9MNz5sPaCLH1+t5caJ9cmr2veZcvlHwyQs17v42Mpp9IHAtoty+ruJV8AvX4FWXrny0Fz5ZshyScA9cLkxYWJM4vjrKJsViZzuq3cMutQl7x2eq/h5s+yPTlZsk7FZEEQz1ZzVo7tJG8SrkueISwV/a76cJBc+cGgQN476FVg3y0mLS6J+Wy11bl30jFrHSGm56dZHs9XkGkQ22nZNE3Y0ol7py1NPmGZk3b4jnI1btMvgYkcD9yecsop8tVXX0m3bt3ko48+kjlz5shZZ50l69evl6VLl0q1atVkl112ifmbevXqGb9LplWrVlK3bt3oV4MGDXz+JLnJn96O0lIJxSKzVmyICQA5PXHe//1oo6co0YjEVBfKyqVPsAaYPugz0xg2udwWBMiWBQ5qECbKKixfjiJ9ukoCueG3vOU2Bz18YefHbMheDFVv2KKzcfPiNHD7+K/OsuTWWDLXrJtP3+vGT4dFhyQXumyUOG3EbLedN0bPXyv/+X60PN9xkixZ5yzDV4/9HpOW5nR9w0pmp5TkisRL6rQRbX+aPePWT2FumIZ40XL282dybgjLDXwiC1Zv9il4VaIsyz6YHTNa4zbF22tmcBDXAOv+4dXbJ/usTd4bIM3jTIqppVXindcSddhrlrY1IOqFTDO2TPaP8cSvE8p+l4X9MMiMsfLJAO7/xvV7Sjj4ud4L0vyLMJz+l6wL5v42aE7uaXX7/DZmkaPX03kKzDZqtkaSpDOq8KVOk6MjHK2jiyIu7vPNZJBstonDch4Jq1AHbi+99FL55z//Kcccc4w0btxYunTpImvXrpWffvopo9dt0aKFUWrB/FqwwHkmSL7zusZpuRpIllIJ8eoaOjnodRiw9hSlm/lpBjesJ9x5q0o+96YcHKLm9Q2hdRP0iTPDuN+BXTeTEoWV2xu/dBrZbn09eK7xums3bZf9d69Z9l5JLpNOt4V1WJ41K02DxZmId46IZ7stWGdmEn03dL6c1spZhu/Jr/SSu78dJe/+FZ4SKk5NKK0HW5BjLZ/Ek5xEyg31TPgatg+bzYzb9CbVFKnoQeVsrCPrp3++42SjFECm7BPfuTn/5/DmyEhZqYjkWfZ+r5+yYKCTjNtgN1a2An7W2ttq/qpNxsSJG+O0gxMdkpqlrSPGnNKM9GzNpdBrSmyN6+FzVnseDPbibyYuWicjHdTPdvM+9n0o2SKUJcOX/O+Rn8fJvz5OPHI10T2AF8ewOZw/0TI6EcZzbQgXKS0vdppslHNxY8yC7Nb1T7budaRB2eOxW8XpSDXz+qDXFKftVLfs+3C80bbJ2n2aJPPZwDnyZIeSzqrTLSMdrZ87VVvLSakEu2hSXb7s9CET6sCtnWbXHnLIITJz5kyj3u22bduMQK7VsmXLjN8lU716dalTp07MF0osttSf9OPiVxa4LR9sUdnIdDAzL6wBprDUYUtHueFYGX4Ea/D69i9HyOwVJUNBCuLUCnY6lD0V6yJ7kb0RRMaMtTHrOnBrb2T7sB/Gu0DXrFY56d84rTNq3WbWYaZuhwrqerO+ltPArX1oayZDFM1hVEHQz/tq16muy4W06loycmFNaXAqV2pDRZI87nS4sv1QeanzFAlrBov1GPS7Ayw39oD4kh2/GvRxcnzY1++zf07MeLmcTlqVD52PyaSz76a6JkayNvw+9XLm8qiLTCxLEDRLZZSLMmNXtxmUsiSHefxneor8bti8uGXS1BYHSRoayE5U+sTL8/fl7w+Ua9smL/HnRMx1P+J8ee2/+mXUQhk2J/NAslv9p68wSsKNX5heoC9X2j1B25jBnBCfD5zjeCIvU4fR8TNZszW60/o2H/ef5VkpE31dnRcmTPpMXWF8NydfN+/Rrff11utbqk2QTsat2c6Nlzxgv8eLJwzZ6WGWU4HbDRs2yKxZs2SfffYxJiOrWrWq9OrVK/r7adOmGTVwTzvttECXM5f5nXFqXlj1gI53wvMi0yHVichsu2lHkt4EvtKl7Eb/qd/KhlIFIZ0LWSXbUZzpGlxqKxcxd9XGuK97SesBcs+3o9J6D73x1vqpZg+n9WM7zWpKJogaddbPkHHGrfjLXLySGUALMr6E/m6ZIMq66t1mzv376xFy6NPdYhryTth7hcM+LDmRd3pOl7b9ZsnUFHWnUmWohL3/KTpsOsmCOs1kcNuw9NIqhx0LVuZHnr7M29po+ZRxmywgqkGfdxxkxdtPwU47gbwYrhh/qLnkDTd7Vlm9u9K/DWi3LLuZTP68ihS4tZ9/rdlobvZfN9t09or4E57Ge69MA3H25apiCdyudpCBf/YbfYzaoDGvaVsm3VcOebqr9I0zOq3sb7IjNuM2VrJd2vxMTrdjwqz5DD/p9NIMcHMEZLr8XN9uXzuM5/2+DtvVft/HmvtkNldRrepV4u6/m7cXuR7F223i0tDNJ6EJVVZbbHV81YM/jHH8essKSwLATqW69dZSl9Z7vKQqxmU4vwK3jzzyiPTr10/mzp0rgwcPlquvvloqV64sN9xwg1Gb9s4775TmzZtLnz59jMnKbr/9diNoe+qppwa96DnLeuH1N+M2/lBYN5lmZuZsukE7XYYWHSbIJ/1nGydtNXR29nuZrdJZ5256DdMJDCcLMOgQr3QsX79Vjn+xZ9yh9MtdXijiSffGK5MeYOtfus3ctj/bz5vbXlOWy3WWYXDJGtvpNDpjMm5dvkCfabENSh3qk+3gnS6y7j/JbmIT0YnYDniic8zIBTeinUou9t/r4mTqhL294+TzOd2mXnT0pKvIwbWmIEG2ywsdJ4mfcjlum+i0oZlgqt2wkklqXZXQ8GA/cdqxHcmz7eFpqYTAlsOs/x1Jmf2YSeB26botoZsYMBH7UtaoWilUE21lXHPV9vfWjNt4I/7imbMyeaB547YdRvDmi0FzQzU5mV2i/f6S1v2l9V8zvHmPsOz2YVkO2xwK2jaMV34u24JsMwXRuWzd909vtHvc53SdmHxupEQd934lDug9cjrsqzReW8rNBGhuJ3iOWEdMxPl9t0lLc/HwDZVQB24XLlxoBGkPPfRQue6662T33XeXoUOHyp577mn8/p133pHLL79cmjZtKmeffbZRIqFDhw5BL3ZOi3ce9fIgMtvD+j1ecMtNwCvRzd3gWascLUPvKctlYGngMOgGh/n+6SyGfT0ku2lI54bk0/6zbe9X9obprjd7RoX1wmpuk0ykU5cn0FIJLmYAzpS14ajvoxN5JXL18X9z/fr77lozWoYh0QSCXrN3AmX2tgVGZ45O3qI34W70n76yXD0pNyqXps9vd5ExPjfD7JQgG+xa2iHRsHKnNxe/lgbzgrBtRyTtc5/fWX1BX9Pi0RtXNzcCD/0YmxnyZ2lmf6pMZ+0Atn/+TG9WNRtMzwtOHPN8D8ln6exbKUslRPwd9mstlaAzaeu+qEPh40l35JdmYJ3aqpc0L5285pnfJ8pprXrJssItaZeA0kxxXdZhs5O3a+NJFUC2/9r+sTc77ajw+Fzj1eRkdhq4zTSobv9zMwBcNU5jx+t2nJ7XkpVhufe7shFw5T5mgkWxju7JdGm9+rTJXifZ5gvjNc/cK+aX1nXuZ0tOCEJYskSzNrLB8jbWt7QfnwW2Y03r+eq5N5mwD85IdUykrHHr8gMabftoZ62rP0U+BG5//PFHWbx4sWzdutUI4urPjRo1iv6+Ro0a8uGHH8rq1atl48aNRtA2VX1blGfOOuhHo6NcUDE6JCcizduPzehEnuh8MyZJIMoaNDCzbK2P+eXwZ7pJt4lLxG+pbg5u/2qEPPeHu1p/yW6S16dZK8leb87r1Z9sYoVk0p3wTll33VTBCS3JoUNGEg5r87FdVWAbRmPWQkp0I+pWoz13Nr6fffCeaWecmEPmnPK6ztSs0rrOmk3jRqbnkSqVSydOzLA1GPYG04Vv9ze+j12wVrpMWBJ3eJbTQEe6mQleeLj9WPnP96NCuY3CWiphi4ObxgJLCRZroGW4ZQKfZAEY7fiw/3pchhOkuJ1o0W2AyG1d61xhXlPLSiV4+zmveH+gq+Wwdg6PTVBL00kmfTyTl5SM0jBnJ/926DyjDvYpr/QyZiLXQLEGAvRrsMP9aV5pqaruk2In2jJpIE9fb0iKhAUn7Nedf348OPp/r+YzcMKrduEQW7B7z9rVy97DozCj2SFkzea182qff7PHNDnsmW6Oaq67mZys7EnBnoMyaX/nSq3bTPYFLybY1DrtmdS4TUeizWq2URa7TJBwK5Jg/SfLTtVl03q+uab/jBXlrkmpO/Ai0nXCkrj3ra4zbiORaNmYMB+HuSzUgVtkxxJL4DZ28hQf3qz0NbXgf7yZa10FbtO8xsf7XD0mx28Up2PUvNVywVt9Yx7TIHGbvmVF0ZU2tluX1uqzDyd0w1oyQodAJHuFATNWytdDYidsCITPPXJuaoR6xXqRSrUbfz9svnS2BKzK17j1fsWkE0Q1swTcMLPm9TOke4zGC+YlU65xkUH7P5N7h0wDrubNn2Zouc32tQq6wXT9J0Ok0/iyusfJvNm9/AzFm7ft8HQ435eD/GmA63m9ywTnnRslW9cMZPm7jXpNDX5IZrrHiPUG3uxgtQc22w2fn3TEhdfHgNvN5eZc8Fq3qXLgk10cPVdnl/56cPlh2dnkZt2aW9Kv/X12iqHs9uXQNpbZzkp0qp+0KPMJKtdu2lZuJFg/y031Nx61w1Zu2Bp3Ii4Vcbkd7dtoomU93PblCMmUzinhjrf7TKrJWN2YVVqr15y8zOx0zZSeNxLdB7jpSC/XpnSwKp/5Y5IsWZdemaeS9/DvmubmpUPaZ2lId9FmLl8vx7bskXGphSOe7S5v90xdI97Lz1g5QaP6qyxdx6z7Q8k5LvVWyNV+1PEL18W970+2TsYsWCv3fT86bmzA7egTfbrZJo53HDq5v4peq129c8VB4Ba23qjsvFe8k4vbmx17jVsNFjhaBg8/ZLyXeq/XzGijLpX2IxbEvl4ay2BdDRtc1K5xwxie5eF6s9fACcsJOqMat5HU+7FmzzT7YnicP078Wl5LVBt6Uho1XZPVPdZVkG4Q1G2jyd64cDscrfdUbzpuUtXDS7XvmZOnPNBujDHkNl0LVm82Xu+N7lOjN/bZpLXCW/zqbKLHeNkWP49a6HiYrhMvdJyc8jmrsrSe7DPufj9snjzQriz73iuZZpiG5RxrPt1+jA+MczNi7cjx+hwaSeNc0HHcYkczKn9k69RNNUGa2eEbJC1jMdKSAe10f0+8Hr3ZYGYd5ETLoYth7huJrk9OS2Ikc1zLnuVn205Q62/4nNWZjxxIY/XZj5Fk7e9EI/NKXsfZm1tHuTnhyzwbHr/e/5WWc6linyXY+p4u3rTRk13kpc5T4r+Oi+WKpNnZ8n7vmalfOxLOgGlY7iPSmRTRaRLFpDTnFQlSpQTZ6K93K99pn9UJ/FJk3OaLW+Pdc1psLZ3ATOsw2yUqJ5TOBInwBoFbxIit/+K9KaVDyfzIWHM6sVgYe9IyGRpWYLupc/MaToe/pTvZkuMJS/LoIplsP/6g90zpF2dG13LD2kp/3KmqNxkiup3NG8VEN6taRsMLvUsz/aybdLda1VL+3Zj5a8p+iESMrFM3pRV0YrAvSoc2fezyxvuOr0bG/Jzu7mh+9nToLpNsuGUif9tlp3KPte03Sxau2Swf9pll1OkKQiTDSaCcTgTlVKpri5MJr7xg7ltmWeanfpsoncZ7V0ZHjyMdCul3/bjZKzak9R5u/8YsV2KfdDTZjZVm3Hp947Vojbtr4NVtBsfM3tx+5AIjKON4RuUkUq3C5eu3yDdD/Mtm0lWrZSz+93NJLVdHgYvSzef3pf6tHtOSdliW5FtFknZiZloCKh7NGE907dXJQuc5uEFOVk5j0KzUpRfGL1xrjPQys0TdBG694PTylo3mYLrvYf87c/LedK7difw62l3ddifbzemmzYemeCSHJsyrCIKuqZtoNGSyNkI+BW4TjfpUOgmz+VnjBdidTCaWcL3Z1qGOuCiLdfjdhZu/CNwiYW0vL4Jp9ppFqYawu0nLT7ceUr6dkK3rYdJid72xWkvSCa2FGn0/yZz1JipfWHerRDVXE33eRKUSdq5RJVpfsdeU9LNCrVkUsx1mg2cuEr1Zr1P6OVIFOkyzVm6UPi6CoD+NXCg3fjpMWpYGKevXrSHpCmoeaA1MmRm3jv+mqDhpJlRYO6r8yM5KJVW94lSTWnp3rva3s0qPo+Y/jfW1M0xr7Z3/Vj/5sM/M6CRMXyUoR2FfjjW2IeTxWC/tWhs0XtsgWek17fTx+tN/4bLchnaiWnl53k3Vhnni1wny7B+Tkk5k5AUnu5g5uqZsmeP/kVcZ9madT+1sNkd/lC5IdDmiGbfirYjLdZWoZmW852qb4IPS4y2etXFey/46HUYvina6BNE2trfZdf+8v93ohKNCPM+a9zExJVXwVCfEy7QOcUEmE946XpmR0JdiSsWPa599QmXHf2fbal6uu0d+Hic3fjrUs9fLV+VLJZT/fz60m3UiS7ftay2PYrat3N5/xBNTQ1gkZg4VazmfPAvDZBWBW9iGeFseD2BZ3NxopHuK0aGpYZskJN1GQaa1/ZyeqLWgvSmdJdQG629jFsatN1fyXULB6WIMmLFCJtjKfaRa73oBSzSMs9ywNtsDN302TO78OjYrNN2bimkuJ/5Kl34Gs95tosymRDqPXyLPuJxEzyrTyQ68mCPDzDI4+KkujoY2a5C7cpLhlvFsSnK+NHvP3dzE6HM1AGc93r2QTiaX18tgzv6dyKqN3pZKKBuOHbszmavCzyCJluoxjz0//Flav/id0v1aJ2F6PkE5Cvu2f6uHu2H+Q0snGNphpig7mJBQb2jDck3xpSsoxWfbXrpuvJ60sfxiRJKOgtCRHvZSCYloBq9XNAh8+qu95XLLxGXWCVHNXdKryZCiUmRw2a+DMxMEUK2sy6jtDj+lu7s4PdTsn3/UvDXGtd5eszk6ibGPdyFpZ9wmWCb7iAC7xq37yw1ZDLKlalMm/LsMVnmm51z7fUHc90j6/pnvL3r/F+8+1MxST1SW6IWOkxx1Pnl5XdL7iWx0OAdBa/lqNqgXrKu8pMZ5+cftQfZfE9yrhdlLaYyu0wkbzWuzF1fDLwaVncv1Za0TvCUKoKeavByxCNwi5gCKd7OXycXQ7YnATcZtup1DPRLMzpvOcGevzyvprGrr9jFv2JwMmXOzvo36p5ZGt9t9QhusD7cfF7feXBi4/Ty3fD5crvggdjbrVKtSZ6B3+v7R3CSP1k+1KpUCOa9EM/gL0uuECIJX9/Jm/Uf9HO/1mhE3oGVtHKaTcZvsRsFJ4+eT/rNktKVEhR7nGoDTOt2pOiF0yK0GZ/yYbE55XSoh1fG5xoMZm53IRmdVybHn3+uvK82adfIZ7G0KJ0F8LfNhPw/YAyP945ScMZ10wG4hKr8T8WX7ahAvUQmj6qXne7N2nVvWLJl4nKzai9/pL7d/OcIyOZnzv83UD3EmrrPm+5qzqlf1aEIp+3vEo5//DdtEjL5nRKfY9+zHiNvOJLfHmNs2uy81bj1+zX+esK/xfb/dann2mp7slWl+zp4OJmo212GPybHDqDMNtJv3JZnKZCme+n2CHPZMt4Sl+JbHGaV45YeD5MtBc+W93uXbeV5M+KtCcznzaVk/GxBb3kyvHU3eG+DR8lizbOM/bj8/PffnpIQTTYZVOhP63njyfgnP40//7myeikTzfOjrWq9xnSw1/3MxozksCNwiprG21DKjaBAXinKzwyeRbrZEqkkgUvWcp5udm2xpM1nX1viWBm6dvNY/PhhkBFOcZsLNXF6WGfJOz+kxvc/p3CDbJycLmhf7eqr1kKxunv0vzYudeQOeytSlyXumq1bO/qle90vzcA6q/ECQbvxsWNKs0d/GLIqpEVmScetuTbXoMCGjLP5XukyVaywlKszzQaps1z/GLnI1y7WZNemGl5OTOTk+UwWrvFBWKMHnhqsxAZN/b+Dm2mu/nKbaxxNdk+wZt8loBnDYb3S1raOdH25mijdFSjsP//HBQBk2e5UU2mrVmx118SZCc0KzA+Oxb/ZU69ja4Zt6cjLvRJJ1VkYiMe0X7Xz6cXjsJLF+0KGp9qGsbs5xXu7P1iB2JiMj3M7NYB82nvJ1bY+bpVm8kd4KtX/W/XaraXyvZusE+PfXI6V7aZJIsvUT7xhNdH51s8SJ5k1IZZWL4db23SXTfTRZ1r2boHAmy/FHisz/ZPeH8Y5n+6YckGRSTSc8HyUQEonOPdqO0TJxXrVn9GUSJcas3BB/3/9nW2cTnwdJy0N1TiNBIiZmYNu3vhua+bwP1rkjrCMhk3UShrzpFjgCt5CFq8uCtZ8OsNRx82s4WRJuhvb5uVgakNUhXKkuFtaTT6rnJvptplkX1veNV+cskf98P1qWOpwAyqijWfo++v8nfysLGOlkK/Z6finZMm7DcpNtLoduf2smoqMGRorXThavsF/EzMDZTtVST06mDf9LWifvmS4IIFi2bUdRdF3FO4ckGpIWtv0gE4//mjiwag+O6g2B20kAk40OMAORbo4tcyhrqswrc9Zrp6+dzmrsmkZAK1Wdbjcdg2Hft+K/dsSScRv8SVVnKbbXKtbtGm+CxlSjRnRYtRvBf/rkZ18zk/irwWVtroEOb+rNbHS90fzXJ0Pl4R/Hxu2oS3dSmES18soCas7Wrp5Ogsi4nbqkUJr+fV85aK+dyx40j7uYju6IXPvRYFfBqmTcfjYnI55SXbvTaZ8nWk7XgVtxvz84WY5Er2bPWHYrYlv3idp4bpgjCuyr7i+HcxLEa3Pp8WetQ9y8/Vh5pUvJNTcdWrJkdY5kDWbCi1NLqrJeyZoQ8e4BnY6CzDVedwwfundt6T5pabl7gj/HLTbKxPVN0mZIv8Zt7PMSzRcxw5K4FFY6wiUdun+a62HeKm/nP4kkm+8lPA20nEPgFjJ+0dq4WUdBFKB3k+pf4OMy/DxqgTFpQqreUWtDN50TkU6kocNyJpROopHOa1iX4ZP+s11tNaeBInuAqMuEsqCKnpjnuJx4xZ4NaN/XzKGM2WJfZ1pEXTMRdUIJp7WuUm27pA1C299u3FrSeNmpamVPSgo4vbdbbMm4z5QGDZLVTtKhSPGGpAVtaeFW2by9ZP97zzKpm9WajduMbLlUNbg6WoYGWXvGNUPO3ruvN4CfWepBZUqH76m5DmYrT1XCRBvp1qx7U++py+SQp7umDEKGYdj62W/0kZs/HxboMuhaMI8JP2rcWjPV/IzbOj2fnPxKL7n/+9HlHk+WOZfoRjdZjUEN6to7QsKwz5WIJF2H1kziD/qkHmqbqrSEld+xeyeruOycEome9+/8aoSvy/XzqIVGuYTosVYckVmlgbCSWQDKQo5e1gH2u9SA04zVstd3/nst5fT1kLn+Lr89cFu6HRJ9Kn8O4ZIX7TtthdHGm7ncXc3/SIJjzE1Ncc0gvPfbUSWvl+DPdPJHU4cxi4y2vavltLzu2a/3katK2wNeSPRJw3LG9fPaZmbcamDd3u75esg86TZxiaMszkxocDPopAfrR9d7ZmsyTzoKt+yQe74dVe6eQNvaxu+TXP/TPWcla1PkmnQ/i94/mtetVNnmbiVrfzi5VoZ9EsSgELhF3HR4PVmaDdpMboDclB1QU5c4b0T5lQmsJxQzczVVANHpDJXKydKmc6KyZm24qRGcaqi13dQ4Qcx0paq/eeRz3SWb7Pu4eSMcr7ZRwu2YMnCb5P0TBOOdDJ1PNemSmxu+uZb6RJnSxlxZxm3532tPehhNXrwumqlm1qm10qHJ75bWrO1YOkmTm0ylU17pJZe2HlCu7q3XGZLmLj1lifMJHsomEpJyHTUXvt0vdpZ2nZRjYcm6ShX8CEsMzaxTF6g0bvRdvnTpBEzev7527Oq50toJlap9MCROmYxkJWDidURpJsjTv0+Me7Oi79/0o8HypO1aFoKE46Sd0eZqs16zUzWXEq3rxBNO+rMS3OxaZiDe3B9XbtgqvVzMI5AJM7CiAUnrKAFz+TVzWTPxveJ2bfs9eVxi5TvMtXSPmxFbZa/iolSCwzZ7tA64B/uv/ZplN89Fx2ayfSz50N/Y3/33hzHSbdJSXyeptL5stuu357LUGbcR4zqogfW2/WeV+/2935XvqPRiPVn/RNuQx7fsGf15XZa2b7L2qnVIfDrHrZP7GG8ybsuOuVe7pp/Fng3ZyNbWUkF+HbbpTjKYn8VAvEPgFnEd/2JPufubkh7hTExf5m6IwWvdpjp+brzra7Z7IWMKnXtwQss04zaT4ZGpfG+5MGfKbMC7bfj7JWLbp6KBZRcX01SNFQ1yJXx/e4Zj9MYquxf3LwZ5l/GpnR7RwG0WL8WDZmZWQ0wDR+Zyx1u3OjT5q9JZsLUBqMO+Gz3ZRX4e6bxOYrwhWWEY2m7etNjPT3NWlpzHl6+PX1qlTZ+Z8v2weQlf16+b01xjzGgc/b/3r19WRzTiuhMv1fVJg6QnvfyXMUtwQYb7bbUkNbfjHXPnvNE37nOPfaFHNKPJOilGmLI1/s9WxsAMXprbytq5nayje+m6LdJhdEltaafCcNhNXFTSceSyD98T5q5pHS1gXSXNfyqrMe4Ft+s7nclkktFSME7af16VqEpUizbx8yPuSiV4sP/azwt2bk9fiSZ0S3dd+nWIevG6CxJMpOV3gDbZyycNBnmwWKn6FrS2uF4H3SYaWWXa1NNrsbUE0bEte0i2ed2me6LD+LiPe72naWk5zbZXsxyOFK1To4r4ZXuSc7aZZT98jn8JBzrqNt4x5cX9SLLX4J4gfQRukdDYBSUlFMJ6eMULBnkx9PrVrlMdf2bricnJiUiHc2pdLS9jNF7XpfFavCCVmUhallkRDmVF2m0/W5jZlnaptmmyG6p0J5JQ253UyXMYNx000/0kUonoEP1oADSLd+03fRZ/OPzr3aYa5Q2cMLP+UpWh0GybwbNWOp6NOZkwBW5/HOFush4tKfHUb+UzIk0h+GihYa1D6zVrECGT14/3p2Ydc71+WbOSfhq50PXrJ5w4Ko2OR/MmupwQ73M6u7K5eNZSCclON3d9MzJmMsN47RaTva6s18o6H5y/wWQXmf9eMY8B6+Sc+lg6q8VaczQRt6+8aesOeaFj2ezlmWradohRuqZseeK36b0aIuz28ybcXTyYjCuR2NEB4rlojdskB1vyQKT793Ry3HkRXD3r9T6u7i9CfMr1LOPWqsfkpVlv6/lRGz8dXrdfErW1te6tGjKr7N5kybrNCedVibffW89TOiGvW7ectr94ST+L3ovoPcP5b/WVU16J34YxS1fqROJ+0REn8XapR5K0NZxKtqtGktWKzvid8xuB2wruoCe7pKw/N3vFRvlxuHfZll5xOQG7Y8sKnc8uHjs5WersCp2sSOtqJcy4Fff8qKHkVrK2zt3fjCz//NJbS3M1BB20sm+OZBmiMxJkkWfSUC5fUzTBL9IulRAMM9CwYcuOQIfV6f71mXXixVK3fTlcHvxhTMKsP/t+aa8Jrb/+uF9Jr7hZpzoRbaQlCxz7fQxop9HbPaYlfZ9Ml0EzvlZtKH/+pHe9TCTNdZIo2zlRJl0mN3jJls1+bkyntl2yEjBOPqcTyYKcQZuyRIcnmhm3kZgyLerWL4bL0bZyQTrBUDKfDyxfA9PvrOOwH9XmMRDTUWApleCGteaojiaxBhLUb2MWRq8FTj3fcbJ8OchdXVk765E0rjQwm4gZuNXJftRyy5wW6fAqczfR62r5hkzZTzUZZxknLJWQbtax8wUyS7dl83Jq3l/oaEZtv9gnlrRmNF9gOUa85uVn1kkzEwXD7EklySQqszJ1afJOKjfXfl3vOieCk5F72ZStWzYdaWK/zz2tVW9jUsls7TteVyy44dOhxvfPBsyWBas3JyxlYo7s85NmbsdrJ3hx7k22nzuqcRv2BkZACNxWcPGGU7aPk23VOY0en3SDAE5rXy4uPaEHyU3gVmvEmkOOD3qqLCsik2CW2RAPM3MG7GSTkwVdaN9+4bIvn9XwOHVPS56byfun/1r2WdvDyL6M2d7eGoiNV1dTh0zFm0AsUfkJneQj0TlOJxpMpmeK7Ay/A7dfDJxjZMYOmLEi4fvNiDNRy5UfDJQ3e0x39B5aR/KEl/6yTGi23qi/FnTHjNVLnSZn5X2iE+9YAgf6SNmESc5fS4NEJ7/cS0bNW+0wkFI24UQ6Um2uTAOCybKaFq9N77puX6JUE4sGqf3IBXGH8ZtZR/2nr5D1W3cYE2qZtQNXlU7Skoj+rTkBl1mKyOuBDubs2tGRMuE5rBO2EbtOWGIEakwlx0V6C26Ornj81/HRm29T5/HpZd8lG6ZubQ/qTXy8tkeyT5K97ePsjcqXhPJfpZjIbSad6xFp/dd0WWI7P5nHsdPzrX1kTsRF+98MtDt5J6/Wrbn6zOy/X0ctjHltrZEZFk6vSzp64T9xJs0sUXruzGAFvtZ1atISZm7KGDV5b4Ac83yPmDI7YRB0m86cc+W6j4dIsy+GRx+Pt1SZLqkmJPjB2gyyTgxv9VaPacaExn7Rjif7pnQ6MtGvUgnUuE2OwC3SnlRAG5Z6gC9cs8nT2ptauN8J++Q+fkg1asba42qeiDKZM83tBWb8wvAHbpMxz90rfGyM2GcbT7o89lmOXWyQTG6SEmVhOHnJ9XGyWb04FuNNzJYue/afZo7oLLLZosemmXVlDsXWmZ0T2afuTsb3hnvUStrgdhMcSzWk3I/JqqxWb9we8/njBc93qlq53GNuMjw04GR14dv9jYa135/Njc8GelfHOR3W4KpTM0sDcjr6Jelrl54xNGDVfVL6pTt0v169cZvxZTJndV6wZpP0mRq7nb0cLRMvA96J8SHJRHLMwebXDDbNaNbzpZPzvH0yRb8ybqMDQkKfc6vn3QXyu2W2bD2F/zDcXSkY042flpThWVA6eal5DM9ftcmoneiVN7pPM4ap67KHWbJTmE4KpvcHTiZOKn868G6/snYSxQskOD0Pb9xWJK3/miGP/Rpbi9N8zZjZ6m2f2foOay1lKnT9WLOknV4SnJVKEM/u/fTzbCxtR9eqXjnmtatW9j7M8lD7sUa7wS+rk7RtzcBoJsdzn2krpG3f8hOXmdyUKrFf88Nyxo03oqfsWPDjHSNG9m3faWUTW2pHpdZ/tWeBl/vLDBfIrzaj9dyUqOTT+71nSteJyTsF+2Qw2ad2EFvXj5elOJLdH63ZtE3aj5hvXDvhDoFblBOvMR7vRsvMaBlsGzKWr/4Yu8i4kdWJiP7z/ahoFpvJi9Od9TyngZXz3+wrkxcXhr5nKtl1Md7kENGM1kji3sZknzuetv3iN5SOeDZ2yGk8H5U2sux1ltxs00xuYu3rz83szk5qQsab2T1VoPc4y6y1mbLX4c32bNraoDRvmHQ2Z816MYeLxnPMvnWN72cctIerQGUmrLUu/WBOPGeer61BaLPhlmkWRbyGms54H5a6bEHT1WOuCfsq0W3w2C/jYjIw3Z7nvVrNupx/f7Gn8WW69N0B0QBpovOJBiLML6/qCOarRNeLeCMArAH0ZOynVfvh+NRvEzzJpjHr06W6Hw7DVo74WLbFPN7OfqOPLxneWlIj2bVHg0vm8PlsadVliowo7SCITqIaid9eVjNXrE+5z+u8AVq3Mvo8Dy8Xml1p3V7pvnRZlnn8TnZr5+T0OCNXYv7G8v9fRy/0qaSQN6+lNW518qurPywZml69igZuI0mTeryQ6YRMek+hgb6/Ji8rN8mak+vPwtLOmXQtTZIlmSyom0pRgmEUfk5gFU+8fdUcWffOX85GaLnV9KPBctuXI5KW5ghRjkBKXh065iRm6bKuMy+TLJK1+/XapqUjtVxdueXxbAnyE4FbOAogxLvQRR+qAEeZNlR0dugTX+ppBDe6TFhabgi413UctZ7m7JUb5bM4tev8bjR5QQMQt8c5KVvrJGrdX12f2riKN0TIDevkLG5pJoUyC97bA8uOeLj5X+4yxfH7OymVkGqSLaspSwo9HwZlf73tO7J70tD3N5eh+8Sl8u3QeUmf36p0X0rUSLZOwOYVP+vAXvH+wKQBfy1x4EXZjUTrY/edq2X0uvnC2MIJhtZqRpdmZT//Z+LJirRDxyx14ec+5MXrJJo8xBDey1bWJOq8ipdxHK/zMx57QN2+Gb8vLbtg0qDfTZ8NjSklkGgm+ZjXzeFGX8Y1TmPKYwW3HrQO8jmv95WL3+mf5FneL9/H/WfLbaVDk5N9/LJLvrPJwfS67GZk3b+/HmEEBnSovrUGqF/ntVHz1sQ9dUUnJ4skbu/EBDrtf2/tQHW4LE6e59Wuqfc99raBtdPX/Dw6aWU2RBz+UjMYT23VS/79zUgje/0JS6Z09NY1yUpyO0lmorIy8ZgdH26Y63ziovgJLfaa236LF+AzE5nc1vpOxrqnxZvo2t1rhMsqj+anyfR6bD0vennvlywIbCZLacZvqnMuYhG4haNh1Wavvg5/MicvM7Nw127e5nnD4bL3BhiZIfEmuskG+7KX1bGKfXzXmtVcfd6Uz7H8flJpxmmH0YmLhIclbhtvOd7uOc0YMpTs+ToM9P9+HOPqgqwBYbNgvdfMxlqyyckSyeR6l8mx4uQYcVML2cys85I2Bqx1orYVZbfGrbUxokN0nE5sqI1vs3GfbukX5WTzuql75pZ94jS92f3BMuFkh9KsHycTM6aTZXxovdoSNhMWrovOWJwtui+ZjWx7ECHZTaR5vtQOnVs+j98ZVvIa4cn+npsk2Gg/v2q20FelGeEVxaatRa6CdE5otu5QS/A21Q3dsDmrZNDMVTFD8lPV0lXmqSqsN8TJjgmncygk0rBFl+jwdr8/f6rAcLKsPjfnQbfMRA5z6eItpbno9hrfsc+Jn7lpfd7bPeNn7/01ZbkRsNXAdbJzolW1yund8mqtdmumn5XZLDA/iwaT7Z1v9o7zggRtE6eBZevTfhkVvwRTpDRTP15SRCZ0M1mvDzWrVTG+66TLYfajZe4Wc/9N1uTKtD1mlhZKVp/VjVSj1EbNX1NutI61Q85r8XIa9Fj1anLR9JXfbkF2sCUru5FqP9inbo2kv9e4zJj5azJu91l3dS9XVbJDaPKSwoSjbQsyLJuV7wjcohzNgrQbVjoM418fDzHqb2mgzWxovdJlqvSeuszTXiAzaKnDVbRx9k7P6YEOt03UC/X3/XaJ/t+LCTt0nel76efN5KSVbNbuMLBmcHcavyS6vZ248O1+Rk96vCHF6bAPHdXtmGxyMjsNgGkgLKNSCQn+1slrull3QbI2POOdY/xk7fnVUgFOJ3jQycvMWb/TyZIwOTknDJiencmUdN++pHX/mBti81zeprRTzosgUpgmsojnig8Gyj3flpS8yRY915n7gn2V2G+ENKhsZtE47Ujy6gbl51GZ19ZMduNqvzxpPcPnO04O/AYrm/w6B17/ydDoZCbOZ7ov+7+TpkMubSd7MMxp2Qmnr+3nuohkOmdCxNl50LUCe5s3knq+AIedUyXPc7YY5v6daq6Hw/Yu6Tg069zHLmdqyerxlwVsJVoD3B6QsXb66rVWR1eYYksWxX/t8stc9vgvSc7VD7QbbWScekk7Laz7jJMMfS94eZSZHR47koyocjoh1Sulo+PsdB/QRCev9EuQBGPtUNd7I6t0sqCd1vbVc1+8ToGWHbMz+asTWhJN909txwdJR0ppUpCuGzfHS6qRks1/GitXtxmc8bFhPc84Hb3ghJlVC28RuEU5my2NikTOeLV3TIPsjq/KNw68uFfXl9DedK2B9duYxJmnXjNn755VWhg+Xk94yczd1p8l4yCbvoYOD9DP60Si9wxjoMQqnXsRDdxr76LJ3khJR7wJqnTdlQ2lSv0aLTpMMLI+MpuczN3jVvN8aDgPmul9ENFaHsHveq6phr0lGm6WdMibz4ucrC6YLr9XWeYaBLSfHsxgRrJhS5nwM5s410Qs6/xjS11u+02kBpXtM9dHn1tUbNRAL19v0ZtlNEtnZGJTknZEohqDx1tq6ua7RWv9C3iYdax1/9DzxjAXNc4TbZvvLOVlohm3Kfa3MJRy8qP2rEk/v5+dkDraKp3RP35zUiXN3Dfc1LPWodbxJs1MxAyIptoPzWXoM215WufIytbPYPs8ZlvbabasfW1Y/87+Ek7a8YnWr2bO+zH/yOj5a8tlemaa8e20FIydTmzUefyS6M9OtoBm+JnXpmTr12l5sWT1RTOZNMrOWgvZz3bXCw4Dr7ru4nUKaOe0l9xM5GY/hrS+dePW/WWji9EtftDzmtaK1nkmPhvgvIzEthRBdHPiukUZ1mO2rrMPeqdO3vCaTn7YzlbGCYkRuIWrk721pzhVPUSvMhHMYWn/+3mcZMuIuSWNETOrM16jTFeT9fH5HgXQUjUAtcfu7R7TjExRv0oG+C2d+7mvh8w1ehe99FaP6fH3f9tQQCcy2dsT/a2TQ2iED5MSJJu4K11bLeURsh3I+3ZI8pq2yWgPtGYzLAnwWNNJhTTL3E/xhrit96j33U2N5WyzzuqdrrgZZwk7YyLlailbbyIHOug00fV52DPdYs5fegPsR+Bz4Zr0rmubkrQPEp3/3UzKmOsyqcnulO5RV7cZJP/6JLYDQCcESZRRFm/b6JDbp3+fGGd4dniP62xZnyQb09c60Q7oyLhkNbHtE4mlEh0BYLaPkmx+8zxnHaJuf7q1FubIuWukZafJ8vnAOY7vHZyuH3Of3rK9OK2RUYlGsOn10WzK/Dxqofw0YkHKQPvjllqr9jqO9ra/tpPiZUBan5boXPqM5Xj1U6YjITV797w3+6b1t03bDpb7242OBkl/tpR8qVG1fHhDRyJo7VtTssBtsmxcp4Lot7J2fKSTwKNzXAQ9J0O8NkGqdakJZ/FqomqQflnA5Ru0Lri5ttzMI5GqPIa5TtKt/RtvW7rpOPPKM39MNMomrt1UkkBCqyI5ArcVULpDxX601ERUz/6ReBIVrw4+80DONrPhaAZu410Axy1cG3PCW71pm7G8mZxE9dXiXQ+tFyQttP9eaa/YBw6GNmdDvPqfyRqwTjNxrEEDpw0Kp6YuLYzW2bFyE1S0DrnPpKPC/Fu9EXnuj7IGt5NXzHRCqWzZur04sIxwbTila8D0FUZ5gGx2HNlnCtabQpM5DDpdiRrczb4oX8PvmwwC3lZeDr/y+mbzn22HZPx6bs4ZiZ46bVniemf206V5M28dhZLWkGcH0h3ZsDJO7bIwZWIGbaVHE5Okuq6YHU7W2tQ6dFQzaM1rtLVDPl4Gn33/fq3b1ITHlY5iyaVSCpnQ69jXg0tK6fglkxJMmm2sI9a0HfXbmIUpJ6BKxRwBUGAvlRDnueYuY62fbN8trJPpmYEyHd4c8XB0oHV5rcvlReBWrx3W62mrrlPibi/r57YfSxpMjj7P9nf63FRDz7VGdZCCHN1nzca8/asR8qalI7N6lcrlnm+/1ifPuM08cOtmZJfZEaCZmZmYtWJDRiPb7Ht6+xHzoyNQred1s23qh3SuHxr4a/pR/LIBG3zsXHPi4L12jp4D3Hy0VKM5vBqNEXE4Usov5j1NBWk2ZIzAbQV002fD0vq7JzpMcPzcvtOWSxcPhkxMs9WKSjf7xy2zppU5qU+867teDKyZZHoh05vcTGqvJrpg6QUpzKyNz2Q0sK03kE6Hv1z1Ydnn1hnXE3Ez/MR0Sev4k3AVFVlLJUQc1/DJ5KKjN1d6E6I3019bgmX6/taaSE7qI7XpO1Matoit2xsG1t7jlzuHpw5WJvXt/Ka1P630uDnllcwyb79KEGiwZ+Dp0KXda5VNvmhOCJMPta4OfLJL6G5ek02wY8+CMDPNslEr2um53U4nAUzEaY3pfOZlrdVEFqzenLA2te5DZnCiraVkR7yYutMRIZqJq6M1sj3DeVA6jFnke+e5F6NTznytjzzcfpynWXA6TD1Zm8esK6/nRs1y1HkwrCM47Fm+5mvpRKpap9kJc9LkVKwBjnQyRK2BW+sIDa1jag3crtlUloGbTnai/XnaFo1XXs18WhjOo06vfdYkHJ0TQjtz09kWZgamBijNEljxasnGu8ewj4pMtuzbLKW9spUxf+Onw+ScN9LLPjZd9t5Az9olS9Ztlsd/1YBoSTvUuns6PfbSYZYodGN2acA63rEWdOlALbc2afG6ckkYdvaSIYlOG9pWdzPhdCqLLclmQSRZmHMomejWT47AbQWUbLZnr+gMrG4CvU5nzdUGaBDiNTAe+2V8TLBDgxOZZtFoFnMmGRZB2eKwoP1/fxxr3EDqJHdOaMNUC+yn6v1OVZNRy0r87rBGcu9py6I3r6myG3VYn7XHP12apa1Z1VVsmR164dYMa5P1/4m83m2a8Xf6mtpgMRs0QVtlOTbGpTGTdUVlDdod+0IP396nwW41y2Uw1KpeMmO0yU0dqlwpa6sdCo0yDOTGC66YD3mRXWo//5k3hFUrF3iWHeS195PUSnNTcxXpe6h94kxK3S/j7bfxbnK1Bl2iOpramW5mPurP1pu/TIf6h92rCSYm8tJgH+rNe+GjfrOiLdWybLJIWXkxy36kdUV1HgxtMyfK8k1nklWn2WHWU3D5jNcio6PB3imuwUVTsjq99vupeBmT1slRk4nYTuPbi4sTTKgWCc1532lQ+vy3ykZu6JwQt30xXB617A+pWOe36DhusVxgeb3L3x/oeeeZF0HxVEPd7cxEIWvgNBM6wanbfcRsr+h724Oz1m3tb/vO/YvrRObW65fWB3c74ZqfnLSd7/nWWYk6nQPnqg8HJR2llW5bzU1NcgSDwG0FlEvHZSbZq15nVsQLKlsbgbVtQY50dJ7gbWH3bHE6ZG1Jac+efdbdZK5pM1heTzA0041EwzvtNDNlTOkEDNrbnM191d6LmsmFWYMjGpCyNpiDdN93o4JehJxkBkP81jvORBqZ3Bg6vVkNg0wzMjQzykqHi5ujQxKMsnXllS6x567R89bG3FzemiRbN4gJQpxksetEQcpJDU54743u0+Ie39aMLXMfa/Je/NEpZmf6zZ8PMzLqzNfLlU6bTFlLTPjFOvw7THRb31HaWW2eH74dOs8YdaadxfECesk6wtOZI8JaDzzZ6Chr51mRrXapBpC1BMQFb/czgs1mwDm2vnvi17Z3QmvSit1fk8tPghuPPWlDh7onG1URhuPMaUa4Bkm1jr4Z6NZJztxMuPVAuzExEwJ7QffbRFm/XmSUZlpO4LRWvTP6+15Tl8vBT3WVGxNMcprqva/7eGjM8RVJowawJs1kY4SK2Yljv36poCcnUxEPj6VkneKZspbayDat//3uXzM8zSbOR5lHmpBz/IrbagN/l52qSpXK+dcf8GKn+EO7rTc+1atW9m1W9rD7z/ejZcLzF0vtGlV9mchjeOlkcXYj5652nHXhZnKpfpZ6gEqHZGmgwUmpgnRps8jL3s50hzjnegAy3+jEQEHJJHBrDmOsCDQzyqTnCOvkfn7MCt+yU0l9+SP2qVOuVqQf/Mj0nrNio+yy05qkJSLgLyc1EE98qWwyn0R0tMhxLcsmxsuhPhukqWyCuhJaOsGc90LLpMS7dHhddsg64un7YfPl5lP3j/s86xk40SVNr1c6A70GiEY9c1FMsDfTchVtHJYKsr+PZgrGmzzTfNr2EFxj3WzTk192X+ZJs7V3qhZbr9arTP57vxtt3LfOffWy6LwXpj/HLZZmpx9gtFv/tstO4jUdOn9k/brRLG1riQTNAk9UVzkdg12WrjFHeVmDpdqBYL09sXcm+y3exGNWydZXGOYBSbX8anYaJSK8ZmYuB6HH5GUyhqBtSvkXYUNgWQLawL/ig0GevqZOrhBmAy3LF4Z6U0HqP922rQpSl77IdOb3a9sOKVe7L5kuaWY0P9R+jNzz7aiUJRkyZa9xlAkd6o7cd3Wb4OpbD5+TurGJ2Ixdez1a66RQTqWqM2Z2yqTqKAsz7RRYEOBNAsS3CRdzsdwTMrNiQ2zbzp7Zmo1AcqIMR2vAKdlIEK1FaXYwW0eRJZoPwWuntYoNbD74Q1mWqZU5IZ7bofh+SHZ982KSQq2PfNRz3cUP5j2bZhnqNdy6nbV0x7VtB8sZr2aW9ZrI/d+PTrgOD3umm5EpG1SAUd/fTgPZ2eqQ6znZ/VwKyYL5ZHDmjuFPXhD0IoQegVt4asoS93Wqwi5Z48OPILjTGX7DRoc55EJmsA7dub9dWaPJiS4Tlko22GtcZWsYEhCPm6GMFdkjP4+LXifajyibRV1pje5Mbuq8qC8YRq26TpVqpTV6kV/S7SBF7vqXZVi19tpnO+lCgzP2yTxN1pFMRQ6CnQe26CwH19tZss062bEan2AuAC1zYtaKDTM3lydz8iarg5/ydhLRRLRebrxybBMXFcYdgeeFuas2Ge13nYvCbUZsELQ+dbyArh+clrVD/vFiMs58R+AWntu0Lb/KBWR7gg0dGpSr68laL6rz+PDevLlZtmwFR3M4BgNUaBqk0NqOOoHKZ5bhuyqd6idOgx4auNVyMbnKbUcVckP3Sc5qeiI/b7j/+8OYwEojxWuvWScMm+1gVJN+lJ2qxg7PDxu/MkG9TG453cUyxqtNag9k++nj/rMT/q6Zi1F96WzHrhOzkxiSK1688sigFwEBsXdGeZG1n28KIqwVKSwslLp168q6deukTp2SmnH5zO9A1Je3nSS3l05akA905u5sNiBy2aQXGkvfaStk1cat0XpncOb4/XaJTooGAKlUr1Ip6eQ1AAAAQNjtXadGTElFs/50vit0EYdkcjJ4Lt9qvRK0dTfE5Zsh84JejJxE0BaAGwRtAQAAkOvSnQenIqFUAjxHoLPiImgLAAAAAADgDQK3FcyQLBRBz7eMWwAAAAAAACDbCNxWMFt2FPn+Hm/3nO77ewAAAAAAAAD5jMBtBZzMBAAAAAAAAEC4EcWrYAjcAgAAAAAAAOFHFK+CWbF+W9CLAAAAAAAAACAFArcVzL3fjQp6EQAAAAAAAACkQOAWAAAAAAAAAEKGwC0AAAAAAACAQBUVR4JehNAhcAsAAAAAAAAgUOs2bw96EUKHwC0AAAAAAACAQO1as2rQixA6eRO4/fDDD+WAAw6QGjVqyCmnnCLDhw8PepEAAAAAAAAAOLCtqDjoRQidvAjctm/fXpo3by7PPfecjB49Wo499lhp3LixLF++POhFAwAAAAAAAJBC9SqVg16E0MmLwO3bb78td911l9x+++1yxBFHSNu2baVmzZryxRdfBL1oAAAAAAAAAFKIRJicLO8Ct9u2bZNRo0bJhRdeGH2sUqVKxs9DhgwJdNkAAAAAAAAApEbctrwqkuNWrlwpRUVFUq9evZjH9eepU6fG/ZutW7caX6bCwkLflxMAcl2lApHi0gvpofVqx/wuIpG4F9wZyzdkbflyxU5VK0u1KpWYMdWmRtWSvuQt28vXtdL1ZbDuZgUi23aUPLda5Url9sHtRbT6cvX40HNNQUGBbmIpMP9fILJ2k7Nj5m+77CSV9UVK6f+N15Gy1zXPTfrcujtVLX0fkYmLStqExzbYpWQ5RGT0/LXGYycfsJvxwJTFhbJnneoye8VGOeOg3WXQzFVJl2e/3WrKbrWqydgFJa9jd+Hh9eSvKcui79toz1rSYfSics+75vi/SZXKuvwln+fHEQskXxxZv47MW7VJGu21s5FpU1QckUmLS7aFrg/dhnpdKY5Eot9XbtgmG7buiG5HfWzJui0BfxIAAJCJSpY2HPIkcJuOVq1ayQsvvBD0YgCB2HfXnaRSQYFxQ7p47Raj+LfeKO+/e03jhlYfN38f+3PJjaJ5A7y8cCtBuQrGDNqqacvWB7koOU2DC1t3FCX8/S41q0YDVHvsXM2IUK7cUNbZmKsuPHwvqVmtivw5brHx88kNd5Phc1Yb/z9x/13lwD1rGecZezDqllP3N845+jtVEoAzz1FlQT0zmFUSbCsw9tGek0uCYU481eRw4+/nrNwo3w+bX+73z19xhGzdUSy1a1SVJ3+bILmgy3/Piq67zwfOlhMP2E3+vt+u8mn/2dJ+5AK59bT95c4zG8otnw+X+as3pXy9RxsfKucdupc0eW+A58s6peUlslM1ZzXNDniic8zPc1+9zPX7tek7U17vNk0GPHZezA3Cd0PnyQG715IzD94j+ljbfrNkt5rV5LqTGjherm/vPNlYr6r/Y+eV+/0Pd50qN3w61Pj/Z81OlCbvDpDJSwrl+3+fIjtXr2Lsh2NKA8Zq1NMXyu47V495v0wDt3rMafDZL5cds490Hr/EOKeNffbi6Od/+rLDpdnpB8jBT3U1Omymvnhp3L/XAG7h5h1SN8ns0hq4rVGlklSpXNK5s2bjNnnnr+nyzOVHyJvdp8nH/WfLjafsJ4fvXVsWrt0sH/ebHd1nJi5aJ5e/P1BOO3B3GTK7JPje9O/7ylvXHWv8/+N+s+T0RnvIFR8MjAbVx1mC7xpI7zBmkRGY12383B8T5esh8+Iu5+xXmkT3M3M9zHqlidGmWrB6k5z1ep9op2iy6+u0ly6RQ5/uFv352hP2lV9GLRSvnH3InvLs5UcY635p4ZboPgwg93V68Mxom6D9iAXy1eC50nCPWsb1xg967v3nCfvK1W0GGz+3v/tU+WzgHBkxd7Xjjtiwe/OfJdcLPWfqbdLW7UXyzB+Tsvb+2nlcv+5Osk/dGvJe75nRtprJWK7S+zf9pv/XJW3914yUr/2fcxuV+zvzgY3bdsh3Q+cbbcjPB87x5bMhzwK3e+yxh1SuXFmWLYu9OdOf995777h/06JFC2MyM2vGbYMGyRvj8NdetavL8vXuAxN71q5unEjiBTVq16hinFiM7AwR2bQtcaCkIlm4ZnO5xzZvL5KpSwnEubF7rWpywv67ysCZKzPat244uYH8MNz/rCnNRlq0tmzb33b6Acbx88Pw+TH7xHmH7il9pq2I+xqvX3tMNGMtGjQrDZZZs+Mq2TLl7vl2lOPlrFq5QL6545SyjDtrllxpbOWaNoPllIa7yQPnH2TcVH5w4/Fy4B47y82fD5PVG7cZz+n58NnG31z4dr9y7/HH/WfIlR8Oivv+w5+6ILr8J7z0V8zvLjlyb+k/Y0V0e09u2dh47mHPdDPW2zv/Os4I6pmfvWGLLuVe/9h968ofD5wZ85gGhV7tWjJC5PYzDpDnrjhSzn+zr8xeuVFGPn1R3ECV3fH77RIT5PHCDSfvJ62uOToa3EiHBjTMgOBnzU4yvr93w/HR3z/4wxjpOG6xESjZf/dacYNRL151VNqfIdV6SxT0a37RIdLsy+Gy587Vo8fDbWc0jP7ei8Dtq9ccLU90cPc61atUMoLH8QJw711/fLntdET9OtH/v35tyc2F0u364AUHyb671jR+rr9LjYSB2x4Pny1vdJ9mBMH19bSTzykzCNXu36fIjZ8NKxeAWrZuq1z2/gDp88i5joO26rC9axvXLA32mhmXbt13TiO544yG5bI6bj51/3LPvfeckpsXNzTIp0HJ3+8/I/rYY5ccagSL1WmNdo95fosmh8lzf06SWqXrQW8GL3ir7PxVvWri9fPcFUfICx0nG/+f06qJce77vx/HyB9jF8uYZy6S41/sGQ3K9Z9esj+Pe+5iI9PYyTGSLj1nauD2gfMOinbe/DVluexas5pUKV3vugyJ6OdIFrRVGuS22rVWNWl5Zck54z/nHWS0/x5rfGg0sGsGbtVRf6trHP+aud9+xHzjZvvlq8vON/eUbvf7zm1kXCv12vFSp8lG4OGFfxxpnCM1cFuzdJu9cOVR0qLJ4cY1Qdugw5+8UI58rpscv9+ucbOHzKzwBruVHVPdHz5bZixbLxe90z/62FkH7yEDZqyMji7Q47dF6blDt709cGu/3ifzrxMbyOJ1m42AtWZ9v3jlUcY6VAfXq23sJ8e+0EOyTc8767dsl7EL1smLnUr2bSs99g9/tiSArct41zcjjY5AXT/pzEZ+8F47J01GuOLY+sa1yg9Njt5bukxYmjDQlur6q9tQ2wA6iufSo/Yx7nvKvkrug4qLS79HIvKPD0raP0fsU8foLEpGOwZ0RICRBV/6em37zkr5d17TZd1Zj6nSzl6kJ96+lG7QVtvrqUY2jZ63RsZa2qatuk412vd6vI2Yu0aCsGvNqrKmNGisiQH16lSXN3tMl/9ecLC81yt1MDPeNUjP5Waik7b/W155pDybIHh79N/qyv3nNZJ7vxud9HWdnnPe+udxxogeff+HLjwkeu+Vij5XO6qf/n1i3N+bbYlkXrrq6Ggnrd6bacBY24tWd53VUD4d4C6we+AeJfcDyLPAbbVq1eSEE06QXr16yVVXXWU8VlxcbPz8wAMPxP2b6tWrG18VkTaiu00qaxycXnrjoCccs1FYv24NOXDPnWVZ4RZPMiq1Mfjr6LJGpWaQaaPdfG0NIJkBoC8GOTuw/++Cg2OyrjSYMm7hOiNzVBsu/7voEGNorXkCVXpxea1bbPkMbeA9c/nhnvaM6TodPKskc+PL206S1n9NN5bN6ps7TpZbv3CXxVC7ehUjMPS/n8eVG2L9V/NzogEqDQYcXG9n2XPnGsZn1xvB3lOXG7/TrJLW1x9nBJYkJsBWGojTnwrKB9/iZbppsPfo51M35jv/90y57L2yhoIGC1es3yo3n7qf0VNn3txvLy6W3v87Vxav3Synv9rbeLzDf06X3lOWy4d9ZyatdTP8yQtkl5rVjEa+ZiTFuxH97T+nR3t849H9/o1/His32YILpi9vP0lu/3KE8X8N2P563+kxv7/vu1FGJtpjv46PPnb+YXsZvZXXti2rt33uoXvK581Oku6TlsoZB+1h3LS+cvXR0SDfr/edJk0/SlyfW7PDqletJHvVriFDZ68y1p1unzo7VYkJFGq2kd78aRBMGyO63p/5faLRuP/uzlOkVunNri7fig1bjQCVbpe96tSQXlOWyZ1fjyz33tedmF4HV5ub/m5kmJ5/aD2ZsGidEWA1g5jWY0OP3ykvXpLy9UY+faFxDtHzljXoNvqZi+RfHw+RYXNWGzecVuYNimZ6aUBFA0VzV22Up3+bGHODq+vV9PO9p8khe9U2Agfbi4qNQIMeKxooGr9wrZE9al8eqwN2r2ncaFStXMkIqmoG3hmNyrL4rEGhzwbMNob9alaf+uK2k6TvtJLj1gwSH/Fs9+jPQ1tcYDQ2dVlmLt8gh+5dWz7pP9sIqukQa/N8rucGPabs5xvdh0fNW2MEir8cNNfYrzQg+PXguUaDU8+XZnkCDW6Mf/5iOab0eNebTA3ymUPKNUPQPG4OqbezTF9Wdt3o9tBZMcttZ66ykrNNbFDOzEjLhGYCDJm1KpplooF+7WxJRc8jnR48y/j/4JkrZeS82BuMCc9fHPf898u9p0mdnaoa+6F5Y2A/f2iD9Me7TzWOtUSB2zMP2iO6nI9cfIhxU6G6/N9ZUqNqZTmj9BwZFSnZTnrN0UxRzWRMRoNIZtBWNb/oULnu4yHlAkUNdttJDqlXWz699UQjgK/vYc0UtJ5rzWNMM0P1tWat2GhsfzNwcNTf6kT3GXPG4P12rykTnm8sbmkwVPdXDfa6Cfha6bGs69JLeqPz5G8TjfNqnRpVy2WSNj5y72jg1sxGalea3X3WwXsa10BToz13jmasmkH7eNc1/RzHNdglGrg1b7Q08Pv8FUcaQTgzYPvs5YfL690qSY/Jy5IGTD+55QS5u7TDTdtsem6yMjNN7TSgeOlRe8tTv000bjpNZiBNbzBbdposlx69d3Q5Kzu4wUyXfsYnmxwe89hP95wWbROZ9Fx3y2kHGF/xPH7JYcaXevryI4yvRJOnmNtJr8u6b85uVb5jqE6NKlK4JXGHw0F77Rzz85btZZ3Dut40YKzZatqREm/mbb12/TpqobzVc3r0OqBBWb1OmOchbUefeuDuMR0IVx3/t7jr8J1/HSsPtx8XbWPqOebM10oyhD+66e9y6dH7GP/XDGdty/04fEG0nTrz5UvloBTno3jnPd3/9Rxzwv67yej5a6LHgdKRArpuBz9xvrG+dRl1u6ql67bIqa16RZ97zzkHyt51akSPj3jMNt2hT3eN2zFmvuc71x1rrNOP+s5K+jms53tte+2yU9Wk91Mf3PB3OXBC+Y5evW7pOVezFP/1SUl2/nUn7is/jYwN1L9w5ZGuzmVmYP+3+0+Xu78ZJf1KO3MSZYjbnX3wHvJipylGm0hHbljvf7RtYE0G+OzWE+XURrsb7Ut930OeTr0v6Huv3bxd7vhqRLS8jF77TA/9OMYIdPf63znyQe+ZMctgmv7SpcZxrW3PX0ctStjZqveB2u7UT6qfV6+tZpvwiUsPi3aqp6Jt0eXrt8jJL5fte2rqi5fIxe/0N+7f9BhMte9cdES9mJFCHR84UzpNWGx0OF1w2F7Guj22Qd2YfcDsuP/HsfWNEU0nHbCr0RmkbUFrNm3bm/9eFsjX78Vlwf1Hfi45vp3S+2wzkK/nwZ9GLpBlhVuNc5vus+Z9sCZZ6HO0xM/4heuMexbzPe2BW41R6GsqvX+Zuyr1KKB0aJvc1HH8YmM96D22tsPdBKVN937nPDlF6X1QqqCtsWwOO4qs5ztTvNGy5v/1uxloto6ktNO4gvkaemxYR+SW/Fw2Olev45r5q/eldloCya3CLfmRje21nA/cKs2ebdasmZx44oly8sknS+vWrWXjxo1y++23B71ooXP5sftEA7cjnrrQaFCkYr0xe7LJYfJKl7KLmN7sm8O7rLTBro13HW6oDUttIPabvlw6jl8SbVxNWrzOyJKz3nTpcLq3e06Xt687Vi4+cm85vVUvo3Gr2Rof9JmZMDvqwQsOdvT57z3nQGPYnV4YRs5dbbyH0psSbYhqhoQuszaAvho0R55P0tAzfXzLCdGMwnZ3nWIMVzj3zb5GYGXvujWMQJ0G1A7fp45MWVJoBIz0QpqINpDnrtwk79p6/Sa80DiaWXDem32NwPq/zzpQLj9mH6ORr41TvXBff/J+MX/3WtNj5Pavhsu3d5wSzaLwws6VCoybg39agpLxHLZ3WcaXOZxw7aZtRqD14iP2NoJK1mXWBrZmS+6xc3VjeXVY7yOND02aFaQBEGUdRmpmceuNgAbM4vUaauZft4lLjAagBjo1kGr66vaT5LYvRxi9h/eXZgt1f+hsady6v9HosPvo5hOM70fvW9cI8miQVPdtHZqs++zVbQYZ20dfSy+YTUpvdJQu28DHz5Pda1U3/kaze7Rxpzfe2pgysxY1AGvNzNGbLrsrj6sv715fltH4atNjjO+btu2QyYsLjSGQ1mNO39sMVprr8dxD9zI+tzbAzE4N6024W9bPqseZ1ivUzFQz40obzuccuqdcn2IYskn3jWQ3TNaAmdLGhK6Tb4fMM44Zpcegfv1y32kyY9kGY9iWflmdpPUs4zT0tHddh8+mWp6+j5YMj163abu80WOqnNJw94Q3Qrr/6LHU9IR9jZ8P2KOW3LZHWYanBon13NHsi5Ih7XqzqttO16Nmc5m956p29apG0E0b1npu0C9r4FdvAk85cHeZuXy9HLRXbSPD157dWq1K7HJqEEpvhMbMX2P87aoNW41G3fP/ONI4bvTmRA8LY59/qacRSDIzbPU6kajjRY8H/Tx6rjQ9eP7B0qLDeBmfRkDPTo9rK22M2zMfUzn9oD2MLytd75rprcFvzUjSzik9tvR4V4OeON9Y3/o7cxid/l6DVrqtzWNN6TWi77QVxn6kgVcdlq0dbP94f6Acv/+ucu0JDaKBW70+WTsIzOvPM1eUfM7zDtsreq3Wa51TWrZixsuXRgO+3955itERZg0GmEFbZe+k0ONOb9q1Lqnuq70sAUhd5ZWlIBoIn75svWxMM0vWpMvlddDVC3pMavAyEQ3gWGmnnX4lojeUZpkO6znIZB77qt+j58bUd9Xnm9d8DSQPmbXSCNbrdU+DKPFoe+KDG/9u/P+yo/eRHpOXGs/X6/htXw43brw1EKP7mRm41VEj2vnUafwSo6NegyUf3vT3mBvPmqXbSjvBzDIE5nX60UvKhnRmg+7r+uUFs12hAXn7Y/ZjxEqzcu0BHO2gqVypUvQ1dL/oOnGJcS7X9qM9yKHthETB4/q77GS0i83ArXkdULqNdDub5wonrj5+32jgVs9pmr2s2eR6DJpBW+MzlA7PffjCQ4ys41TrQbd/30fPNWoZf9hnpjxy8aHRjEDr331449/ljjPWGO1i3ZfNto9+Tju9llg7G1pcergR+E4WuNWMMG2XJgramuV8jBFEZx9obDttW35089+NwKdZZsOknfgm7fRquHst4zr3YufJMRmjf99vF7ng8Hox7YJr/v43I+CmAXHznKvXXOu20KCdJr3857xG0nXCUtfnwjeuPUa+Hz7fCPpr55A1cKvbbuGaTUnrS2r73TyOX7v2GCMbXQOj2omr60iDLm91n2aUDdH2sraZzMx4bQ+b7Xnt5NNOJDN7/N9nNjTan/remkGo2eSacKHtcavW1x8vra8v+b/e++i61WtQvPWgn1E7yDSTUkt/2GnwN2Zfu+nvctWHg+SYfesa21qDpdpG1nOXudwaALUG3vS9lT5PP98hT3U1Oqv03KnLZJbKOb7BLsZ9oGYwz1250Qii6TbWfV87vLUz/u6zDzT2mRNe7CnPXnGEcU+ho0o1cHvMvrvI57edFO1sNQN2v/2nbFSHdUSTHivaxnPCSeDWTLbQjj3z/tn0v4vdn8fNfUjXq7aD2t5Sch9l6j11mdzx1chyHRmaEGUO87dmc359x8lGG9nMatUAaSbllEz2ziAdQaQJIPZM9rKgeMn/9ZyibUO9Jt579oHSZeJSeeKSw+Rvu5bUY9d7/WFzVsk3peV19J66JJhecl4fv2id0XGgx46OAtP2o9K25wPtxhj/1/sobWPGC8ZHl6X0Ne3LqZcObbNp3EUTTcxAuSZPaLteDwvtcy15bsnfFZn/L3285OfYzgB7SSFNfHEySlCPOV0efR1dBsQRyRPvv/9+ZL/99otUq1YtcvLJJ0eGDh3q+G/XrVtnlO7Q7/lu6/aiyP6PdzK+lhducfQ3ExaujYyZvyayZuNW4+c7vhweea3rlMj8VRuNnxes3mi8XsdxiyKrN2yNvNJlcmRHUbFnyzxq3upIkYev58bXg+dE19c/2w6OFG7eFpmyZF3k1Ff+ij6+Yv2WyNQlhcb6TPS5dR1u3rYj5rGJi9Yafzd63mpj/ZrrU63duC0yb+VG4/2uaTMo5ndq244iY1uGgS7fZwNmR2YtXx8ZMWdVZOm6zcby3fjpkMjY+WuM57QfPj/S7Ith5daB2/1AX2Pd5m2RN7pNja5/fR+7VRu2RtZv2R7pPWWZ8XyTrrM/xi6KvPfX9MigmSuij4+cuypSXFwc3Va3fTHM2Od037bve6e36hXpPXVZyuX9fui86DGjNm3dERm/YK3jz6vPN+myfTFwdsr1p/uU9e8yNXflhshvoxca7+vlMW1d3sbv9IusXO/sXJSOLdt3RLbvCMexki26Po99oXtkzooNMY/r8anHYkW2ZO3myNs9phnH1LgFayIDppedB7z21+SlkQ1btsc8Zj+XLyvcbBxfM5evN66fiejymuco9fmA2ZEvB842/u/lMf/VoDnGMjllnofhbh0//+dER8+dsWx95M6vhvty/jXpNVKvVXquTPU+1n1t8MyVkUd/Hpv0+bpvv9hxUmTj1tjjIJ/oZ7Qem0qPiU/7z8r4tbUdc8KLPSKTF69Leqxd9HbfaLvofz+VbZPn/pgYeanTpIgX3u81PXLX1yMcP//DPjMi93wz0vi/uez2r//+MDru+tR2WKb09bV9Yf1Zv85/s0/khT8nRVp2nBR5oN1o4zE9L2vbVf9/3Avdje/nvdHHuF58M3hO5PqPy9qZeoyc+VovY/83/eODgdHXf/b3CcZjx7fsUW576f3DWa/1jrQbNq/ctSDVufTlzpON84buE8c83z0yaZE396zm59avv7fsYdx3uDV8zqpy57SP+s40XrPL+MXlnt9v2vLIBW/1jR43um/pc+Pd69nXU7r0HNRz0lLj/9rG1/ebtrTQ1WvofZ55LjPXWbxjwrx/cKr7xCWRf7w/wLh3ikfbBrrNZ9vadF5ef/W403ajbkvz2F20ZlP0Hl+/dPn0ntVry9aVtIPi0Xu1eJ/TfEzbcfr9xJd6GtvGfNxsM+n/tS2WKT0XPPzjGOP19Dj2it6f6DVer/V2eh7U99PzuG4L/b95TtXj9ueRCyJe+rjfzGhcwwt6vNnbwMgsDlmg/0gFpzVu69atK+vWrZM6dcqyA/PVvFUb5b1eM40eV2bsS04PDx1Gbh0+jfCYvUJ7qGtEh/0DALLnrR7TjOw362gFAN4zs/0yyRwLUqMnuxjZVJrFqBPGmZnk1trhXtJSDToyJVU26g4tg1S5kmzeVmTUzNWSL9Zh+enSskM6SZ5OyueEjvjRclHWkQ3Zcl3bIXLYPrWjtaG9oCN6Lny7v1HOIl5mdK7TLMUek5ZKt4fODmwZBsxYIfvU3alcaRWv6QhZHV10a4JSMn7S0V0638SfD5xhZBzHOx+ax7DftGyW1up3M5FqJjR79b3eM+T2MxoapRx0QtVr/r5vXh5PFVmhizgkgdsKGLgFAAAAkBt0CLCWlurZ/BzJ1U6e93vPNAItd38z0qj52vt/5ziaRCcb9HZYy5q9dNVRcScoBBAeN302VFZt2BZo4BzwAoFblwjcAgAAAAAAAAhTHNL/vHIAAAAAAAAAgCsEbgEAAAAAAAAgZAjcAgAAAAAAAEDIELgFAAAAAAAAgJAhcAsAAAAAAAAAIUPgFgAAAAAAAABChsAtAAAAAAAAAIQMgVsAAAAAAAAACBkCtwAAAAAAAAAQMgRuAQAAAAAAACBkCNwCAAAAAAAAQMgQuAUAAAAAAACAkCFwCwAAAAAAAAAhQ+AWAAAAAAAAAEKGwC0AAAAAAAAAhAyBWwAAAAAAAAAImSpBL0AYRCIR43thYWHQiwIAAAAAAAAgT5nxRzMemQyBWxFZv3698b1BgwZBLwoAAAAAAACAChCPrFu3btLnFESchHfzXHFxsSxevFhq164tBQUFQS8OKmhvi3YcLFiwQOrUqRP04gBZxzEAJMbxAcTHsQE4w7ECJMbxgSBoKFaDtvXr15dKlZJXsSXjVgv9Vqok++67b9CLARgXCi4WqMg4BoDEOD6A+Dg2AGc4VoDEOD6QbakybU1MTgYAAAAAAAAAIUPgFgAAAAAAAABChsAtEALVq1eX5557zvgOVEQcA0BiHB9AfBwbgDMcK0BiHB8IOyYnAwAAAAAAAICQ+f/27gemqioO4PgPFTSw8E8vxRTUFG1ppOVQs5VUQjPFyLByC82ZYq50uQrXStcf0wrMamUtyf5Mgk2ltD/GP1MjA4YFZogGuRaIuaCMAIPbfmd7bw8FBQJ9f76f7Qq8d+85597t5733d889hx63AAAAAAAAAOBiSNwCAAAAAAAAgIshcQsAAAAAAAAALobELdCKNWvWyIQJE+TSSy+VK664QmbNmiUlJSXN1qmrq5OHH35Y+vfvL71795a7775bjh8/3mydRx55RK6//noz2Pl1113XYl1ffvmlTJw40dRls9lMOeXl5edtY1pamowePVp69eolY8eOlc8++6zZ91u3bpVp06aZ9vn4+MiBAwc6dCzgnTwhBlatWmW+DwgIkL59+8ptt90m+/fv79DxADwtPubNm2fODc5LVFRUh44H4GnxcWZs2JeXXnqpQ8cE8NRY0bbo+WTQoEHi7+9vziOlpaUdOh6Au8TGwYMHzXpDhw4154b169eftc7XX38tM2bMMLGh62zfvr1DxwIgcQu0Yvfu3eZE8O2338pXX30lp0+fNknQv//+27HO8uXL5dNPPzUXNLr+b7/9JjExMWeV9eCDD8qcOXNarKesrEyio6MlIiLCJFb1xPH777+3WI6zb775Ru677z5ZsGCBFBYWmpOZLsXFxY51tK1TpkyRtWvX/q9jAe/kCTEQGhoqr7/+uhQVFcnevXvNxZXuw4kTJ/7XsQE8IT6U3mBXVFQ4li1btnT4mACeFB/OcaHLpk2bzI233qgDncXdY0XnOde/f/75Z0lPTzfrhISEmAflzvsAeFps1NbWyvDhw+XFF1+UgQMHtriOtjUsLEzeeOONdu8/0IwFoE2qqqosDZndu3ebv6urqy1fX18rLS3Nsc6hQ4fMOrm5uWdt/8wzz1hhYWFnfa7b9+jRw2psbHR89sknn1g+Pj5WQ0NDq+2JjY21pk+f3uyz8PBwa9GiRWetW1ZWZtpVWFjYjj0GPCcG7Gpqakz7MjIy2rDHgGfHR1xcnBUdHd2BvQU8Pz7OpLESERHRhr0FvCdWSkpKTFuKi4sd32sdNpvNeuedd9q9/4C7xIazkJAQKykp6ZzraLu2bdvWpvKAM9HjFmijmpoa87Nfv37mZ0FBgXnyp0+U7fQ1ouDgYMnNzW1zufrqRrdu3SQ5OVkaGxtNPR988IEp19fXt9XttA7nulVkZGS76ga8KQYaGhrk7bfflsDAQPP0G+hM7hofOTk55hXEUaNGSXx8vJw8ebLNbQM8PT7s9NXbnTt3ml6HQFdyt1ipr683P3UYBTutR19L1zedAE+NDeBCInELtEFTU5MsW7ZMbrzxRhkzZoz5rLKyUvz8/KRPnz7N1h0wYID5rq2GDRsmu3btkpUrV5qLHC3v119/ldTU1HNup3VoXf+nbsAbYmDHjh1m3Cu9qUhKSjKvW11++eVtbh/gqfGhwyS8//77kpmZaYbU0dcM77jjDnPjAnh7fDjbvHmzGfvwfK/OAt4WK/ZEWUJCgvzxxx/mIbmeT7RsHWIE8NTYAC4kErdAG+j4OjqWU0pKSqeXrSeWhQsXSlxcnOTl5ZkbZz0JzZ4924wbdezYMZN0si8vvPBCp7cB8OQYmDp1qhmzSsdp00RVbGysVFVVdfp+wHu5a3zce++9MnPmTDPZjI5RqA85tA7thQt4e3w40/Ft586d26xXIdDZ3DFWtEeiToZ8+PBh0xNSJyfLzs42DwG1FyPgrbEBdKYenVoa4IGWLl1qbmZ1VsjBgwc7PtdByPWpcnV1dbMnffo6XWsDlLdEByvXV7fXrVvn+OzDDz+UIUOGyP79++WGG24wSSc7++shWseZs2a2t27AG2IgICBARowYYRadMXbkyJHy7rvvmt4hgLfHhzOdZEN7ox85ckRuvfXWNrcR8OT42LNnj5nJ/OOPP27HngPeEyv6qrluq6+Ya1ttNpuEh4ebMgFPjQ3gQuIxGNAKfcKmJ4pt27ZJVlaWeY3CmV6k6FNmfcXUTi/s9ancpEmT2lyPzkh55hPp7t27O14L6dGjhyPppIv9ZKF1ONet9BXw9tQNeGMMaJn2MdmAjvLE+NBXA3WM26CgoDa3D/D0+NAHfdpexkZHV/CkWNHklyZtS0tLJT8/X6Kjo9vcPsDdYgO4oM6argyAER8fbwUGBlo5OTlWRUWFY6mtrXWss3jxYis4ONjKysqy8vPzrUmTJpnFWWlpqVVYWGhmXw0NDTW/61JfX2++z8zMNLNWrl692jp8+LBVUFBgRUZGmtkpnes60759+8wMmC+//LKZQVNnytSZNYuKihzrnDx50tS1c+dOM5NlSkqK+Vv3A/D0GDh16pSVkJBgZpYtLy837Zs/f77Vs2fPZrMfA94YH3/99Ze1YsUKEx9lZWVWRkaGNX78eGvkyJFWXV1dlx03eAd3jw+7mpoay9/f33rzzTc7/RgBnhIrqampVnZ2tnX06FFr+/btpsyYmJguOV7wHq4eG7q9vaygoCBzTaW/a312eq1lX0fvxRMTE83vv/zyS5ccM3guErdAK/Q/15aW5ORkxzr//POPtWTJEqtv377mwv6uu+46Kyl68803t1iO3ijbbdmyxRo3bpwVEBBg2Ww2a+bMmebi6Hz0QklPQH5+ftY111xjErTOtK0t1a0XXYCnx4C2TdszaNAg871eVGm53333XacdI3gvd48PvRmZNm2aKU9vwvUGZeHChVZlZWWnHSN4L3ePD7uNGzdal1xyiVVdXf2/jwngqbHy6quvWoMHDzbnEk2iPfXUU46kGOCpsaHbt1Su1menDzRaWicuLq5TjxU8n4/+c2H7+AIAAAAAAAAAzoUxbgEAAAAAAADAxZC4BQAAAAAAAAAXQ+IWAAAAAAAAAFwMiVsAAAAAAAAAcDEkbgEAAAAAAADAxZC4BQAAAAAAAAAXQ+IWAAAAAAAAAFwMiVsAAAAAAAAAcDEkbgEAAIAuMHToUFm/fv3FbgYAAADcFIlbAAAAuLV58+aJj4+PWXx9fWXAgAFy++23y6ZNm6SpqanN5bz33nvSp0+fdtff2nZ5eXny0EMPtbs8AAAAQJG4BQAAgNuLioqSiooKKS8vl88//1ymTp0qjz76qNx5553y77//XpQ22Ww28ff3vyh1AwAAwP2RuAUAAIDb69mzpwwcOFCuvPJKGT9+vKxcuVLS09NNEld7xKrExEQZO3asBAQEyJAhQ2TJkiVy6tQp811OTo7Mnz9fampqHL13V61aZb6rr6+XFStWmLJ12/DwcLP++bY7c6gE/W7jxo0mmawJ3auvvlpyc3PlyJEjcsstt5iyJ0+eLEePHm22b7ofuk+9evWS4cOHy+rVqy9aMhoAAAAXDolbAAAAeKSIiAgJCwuTrVu3mr+7desmGzZskIMHD8rmzZslKytLHn/8cfOdJkw1yXrZZZeZnru6aLJWLV261CRYU1JS5IcffpB77rnH9PAtLS0953YtefbZZ+WBBx6QAwcOyOjRo+X++++XRYsWSUJCguTn54tlWaY+uz179pj1tffwjz/+aBK/moh+/vnnu/z4AQAA4OIicQsAAACPpclRHT5BLVu2zAyhoD1hNan73HPPSWpqqvnOz89PAgMDTa9Y7bmrS+/eveXYsWOSnJwsaWlpctNNN8lVV11lErNTpkwxn7e2XWu0d25sbKyEhobKE088Ydo2d+5ciYyMND1wNUFr782rtHftk08+KXFxcaa3rY7dq8lfTeACAADAs/W42A0AAAAAuor2YNWkqsrIyJA1a9bITz/9JH/++acZbqCurk5qa2tbHYu2qKhIGhsbTaLVmQ6f0L9//3a359prr3X8rpOoKR2+wfkzbZO2T3vxfv/997Jv375mPWy1PedrNwAAANwfiVsAAAB4rEOHDsmwYcNMz1YdWzY+Pt4kQfv16yd79+6VBQsWSENDQ6sJUB0Dt3v37lJQUGB+OjtXz9rW+Pr6On63J5Rb+qypqclRv/a6jYmJOassHfMWAAAAnovELQAAADySjmGrPWaXL19uEq+aDH3llVfMWLfKPkyCnQ57oL1ZnY0bN858VlVVZYZKaElL23UWnZSspKRERowY0SXlAwAAwHWRuAUAAIDb06ELKisrTQL1+PHj8sUXX5hhEbSXrU7uVVxcLKdPn5bXXntNZsyYYYYfeOutt5qVoWPfag/XzMxMM6mZ9sLVIRJ0DFotQ5O+msg9ceKEWUeHPZg+fXqL23XWEAZPP/202Yfg4GCZPXu2STrr8Am6PzpGLwAAADwXk5MBAADA7WmiNigoyCRRo6KiJDs7WzZs2CDp6elmiANNqCYmJsratWtlzJgx8tFHH5nErrPJkyfL4sWLZc6cOWKz2WTdunXmc52ETBO3jz32mIwaNUpmzZoleXl5Jpl6ru06g05atmPHDtm1a5dMmDBBJk6cKElJSRISEtJpdQAAAMA1+Vg6YwMAAAAAAAAAwGXQ4xYAAAAAAAAAXAyJWwAAAAAAAABwMSRuAQAAAAAAAMDFkLgFAAAAAAAAABdD4hYAAAAAAAAAXAyJWwAAAAAAAABwMSRuAQAAAAAAAMDFkLgFAAAAAAAAABdD4hYAAAAAAAAAXAyJWwAAAAAAAABwMSRuAQAAAAAAAMDFkLgFAAAAAAAAAHEt/wFtmePGKJSqUQAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 1400x500 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# Preprocessing\n",
    "date_str = merged_df[\"Date\"].astype(str).str.strip()\n",
    "time_str = merged_df[\"time\"].astype(str).str.strip()\n",
    "\n",
    "# Handle non-standard \"24:00\" timestamps.\n",
    "is_2400 = time_str.str.startswith(\"24:00\")\n",
    "dates = pd.to_datetime(date_str, dayfirst=True, errors=\"coerce\")\n",
    "time_fixed = time_str.where(~is_2400, \"00:00\")\n",
    "dt = pd.to_datetime(\n",
    "    dates.dt.strftime(\"%Y-%m-%d\") + \" \" + time_fixed,\n",
    "    errors=\"coerce\"\n",
    ")\n",
    "dt = dt + pd.to_timedelta(is_2400.astype(int), unit=\"D\")\n",
    "\n",
    "# Write back to the dataframe\n",
    "merged_df[\"datetime\"] = dt\n",
    "pm_col = \"PM<sub>2.5</sub> particulate matter (Hourly measured)\"\n",
    "\n",
    "# Draw the plot\n",
    "plt.figure(figsize=(14, 5))\n",
    "plt.plot(\n",
    "    merged_df[\"datetime\"],\n",
    "    merged_df[pm_col],\n",
    "    linewidth=1\n",
    ")\n",
    "plt.xlabel(\"Datetime\")\n",
    "plt.ylabel(\"PM2.5 (µg/m³)\")\n",
    "plt.title(\"Temporal variation of PM2.5\")\n",
    "plt.tight_layout()\n",
    "plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "obsaq39v2",
   "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.20"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}