admin_first_time_setup.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "910686c3-6e3b-42ac-9214-6095bd98cb02",
   "metadata": {
    "tags": []
   },
   "source": [
    "# initialize data for this repository\n",
    "Eric Gagliano (egagli@uw.edu)   \n",
    "Updated: February 20th, 2024"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f65f8deb-f078-4348-969f-99b8e2bdda79",
   "metadata": {},
   "source": [
    "**This notebook is not meant to be run by the user, but is meant to document workflow that is used to pull in the SNOTEL and CCSS data and create the CSVs. The auto-update script and github action will keep all station data and the station geojson updated. To reinitialize and start over with a fresh set of data:**\n",
    "1. remove old data\n",
    "    - `!rm -rf data/`\n",
    "    - `!rm all_stations.geojson`\n",
    "    - `!mkdir data/`\n",
    "    - `git add` / `git commit` / `!git push` deletions\n",
    "2. set up snotel_ccss_env for this notebook\n",
    "    - `!mamba env create -q -f environment.yml --force`\n",
    "    - switch kernels to snotel_ccss_env\n",
    "3. run this notebook and push updated data\n",
    "    - run section 1\n",
    "    - run section 2, and keep re-running until no more stations added\n",
    "        - there are some stations that there is truly no data for, but often stations will fail and a re-run will bring the data in\n",
    "        - I usually run ~4 times before no new stations are added\n",
    "    - run section 3\n",
    "    - `git add` / `git commit` / `git push` both `data/` folder and `all_stations.geojson`\n",
    "    - might have to use [git gc](https://git-scm.com/docs/git-gc) to clean repo\n",
    "4. check automatic workflow\n",
    "    - make sure github action is working and running daily \n",
    "    - check most recently updated CSVs and make sure start and end dates correct"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "22422f0b-7216-4131-874f-b7838d9a7c29",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# !rm -rf data/\n",
    "# !rm all_stations.geojson\n",
    "# !mkdir data/"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "3b8a4568-5cce-4fd8-bf09-0907b898781d",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "#!mamba env create -q -f environment.yml --force\n",
    "# now switch kernels"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4aae0cfe-b9fd-45e2-9dab-faaf1e00baa8",
   "metadata": {},
   "source": [
    "## 1. Create separate SNOTEL and CCSS station geodataframes"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "ed1373a9-1542-4166-8d7a-eb45755e0e3d",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import geopandas as gpd\n",
    "from shapely.geometry import Point\n",
    "import requests\n",
    "from urllib.parse import urlencode\n",
    "import matplotlib.pyplot as plt\n",
    "from datetime import datetime\n",
    "from io import StringIO\n",
    "import numpy as np\n",
    "import mgrs\n",
    "import datetime\n",
    "import os\n",
    "import time\n",
    "import glob\n",
    "import ee\n",
    "ee.Initialize()\n",
    "import sys\n",
    "sys.path.append('snotel_ccss_stations')\n",
    "import snotel_ccss_stations\n",
    "import shapely\n",
    "import ulmo\n",
    "import fsspec\n",
    "import tqdm\n",
    "\n",
    "wsdl_url = 'https://hydroportal.cuahsi.org/Snotel/cuahsi_1_1.asmx?WSDL'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "1771c2f0-4055-43f0-a816-c9a9a4b2802c",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/eric/repos/snotel_ccss_stations/snotel_ccss_stations/snotel_ccss_stations.py:62: UserWarning: Could not infer format, so each element will be parsed individually, falling back to `dateutil`. To ensure parsing is consistent and as-expected, please specify a format.\n",
      "  stations_df['endDate'] = pd.to_datetime(stations_df.endDate)\n",
      "/home/eric/repos/snotel_ccss_stations/snotel_ccss_stations/snotel_ccss_stations.py:63: UserWarning: Could not infer format, so each element will be parsed individually, falling back to `dateutil`. To ensure parsing is consistent and as-expected, please specify a format.\n",
      "  stations_df['beginDate'] = pd.to_datetime(stations_df.beginDate)\n",
      "/home/eric/repos/snotel_ccss_stations/snotel_ccss_stations/snotel_ccss_stations.py:70: FutureWarning: Downcasting behavior in `replace` is deprecated and will be removed in a future version. To retain the old behavior, explicitly call `result.infer_objects(copy=False)`. To opt-in to the future behavior, set `pd.set_option('future.no_silent_downcasting', True)`\n",
      "  stations_df['isActive'] = stations_df.isActive.replace({'True':True, 'False':False})\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 3.28 s, sys: 83 ms, total: 3.36 s\n",
      "Wall time: 5.54 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "snotel_stations_gdf = snotel_ccss_stations.all_snotel_stations(from_ulmo=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "5ca4b458-6606-4054-94ef-9ace1da69256",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "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>name</th>\n",
       "      <th>elevation_m</th>\n",
       "      <th>latitude</th>\n",
       "      <th>longitude</th>\n",
       "      <th>county</th>\n",
       "      <th>state</th>\n",
       "      <th>beginDate</th>\n",
       "      <th>endDate</th>\n",
       "      <th>HUC</th>\n",
       "      <th>isActive</th>\n",
       "      <th>mgrs</th>\n",
       "      <th>mountainRange</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>code</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>301_CA_SNTL</th>\n",
       "      <td>Adin Mtn</td>\n",
       "      <td>1886.712036</td>\n",
       "      <td>41.235828</td>\n",
       "      <td>-120.791924</td>\n",
       "      <td>Modoc</td>\n",
       "      <td>California</td>\n",
       "      <td>1983-10-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>180200021403</td>\n",
       "      <td>True</td>\n",
       "      <td>10TFL</td>\n",
       "      <td>Great Basin Ranges</td>\n",
       "      <td>POINT (-120.79192 41.23583)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>907_UT_SNTL</th>\n",
       "      <td>Agua Canyon</td>\n",
       "      <td>2712.719971</td>\n",
       "      <td>37.522171</td>\n",
       "      <td>-112.271179</td>\n",
       "      <td>Kane</td>\n",
       "      <td>Utah</td>\n",
       "      <td>1994-10-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>160300020301</td>\n",
       "      <td>True</td>\n",
       "      <td>12SUG</td>\n",
       "      <td>Colorado Plateau</td>\n",
       "      <td>POINT (-112.27118 37.52217)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>916_MT_SNTL</th>\n",
       "      <td>Albro Lake</td>\n",
       "      <td>2529.840088</td>\n",
       "      <td>45.597229</td>\n",
       "      <td>-111.959023</td>\n",
       "      <td>Madison</td>\n",
       "      <td>Montana</td>\n",
       "      <td>1996-09-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>100200050701</td>\n",
       "      <td>True</td>\n",
       "      <td>12TVR</td>\n",
       "      <td>Central Montana Rocky Mountains</td>\n",
       "      <td>POINT (-111.95902 45.59723)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1267_AK_SNTL</th>\n",
       "      <td>Alexander Lake</td>\n",
       "      <td>48.768002</td>\n",
       "      <td>61.749668</td>\n",
       "      <td>-150.889664</td>\n",
       "      <td>Matanuska-Susitna</td>\n",
       "      <td>Alaska</td>\n",
       "      <td>2014-08-28 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>190205051106</td>\n",
       "      <td>True</td>\n",
       "      <td>05VPJ</td>\n",
       "      <td>NaN</td>\n",
       "      <td>POINT (-150.88966 61.74967)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>908_WA_SNTL</th>\n",
       "      <td>Alpine Meadows</td>\n",
       "      <td>1066.800049</td>\n",
       "      <td>47.779572</td>\n",
       "      <td>-121.698471</td>\n",
       "      <td>King</td>\n",
       "      <td>Washington</td>\n",
       "      <td>1994-09-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>171100100501</td>\n",
       "      <td>True</td>\n",
       "      <td>10TET</td>\n",
       "      <td>Cascade Range</td>\n",
       "      <td>POINT (-121.69847 47.77957)</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>877_AZ_SNTL</th>\n",
       "      <td>Workman Creek</td>\n",
       "      <td>2103.120117</td>\n",
       "      <td>33.812420</td>\n",
       "      <td>-110.917732</td>\n",
       "      <td>Gila</td>\n",
       "      <td>Arizona</td>\n",
       "      <td>1938-12-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>150601030802</td>\n",
       "      <td>True</td>\n",
       "      <td>12SWC</td>\n",
       "      <td>Southwest Basins and Ranges</td>\n",
       "      <td>POINT (-110.91773 33.81242)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1228_UT_SNTL</th>\n",
       "      <td>Wrigley Creek</td>\n",
       "      <td>2842.869629</td>\n",
       "      <td>39.132332</td>\n",
       "      <td>-111.356850</td>\n",
       "      <td>Sanpete</td>\n",
       "      <td>Utah</td>\n",
       "      <td>2012-10-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>140600090303</td>\n",
       "      <td>True</td>\n",
       "      <td>12SVJ</td>\n",
       "      <td>Colorado Plateau</td>\n",
       "      <td>POINT (-111.35685 39.13233)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1197_UT_SNTL</th>\n",
       "      <td>Yankee Reservoir</td>\n",
       "      <td>2649.321533</td>\n",
       "      <td>37.747971</td>\n",
       "      <td>-112.774948</td>\n",
       "      <td>Iron</td>\n",
       "      <td>Utah</td>\n",
       "      <td>2012-10-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>160300060202</td>\n",
       "      <td>True</td>\n",
       "      <td>12SUG</td>\n",
       "      <td>Colorado Plateau</td>\n",
       "      <td>POINT (-112.77495 37.74797)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>878_WY_SNTL</th>\n",
       "      <td>Younts Peak</td>\n",
       "      <td>2545.080078</td>\n",
       "      <td>43.932251</td>\n",
       "      <td>-109.817749</td>\n",
       "      <td>Park</td>\n",
       "      <td>Wyoming</td>\n",
       "      <td>1979-10-01 00:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>100800130103</td>\n",
       "      <td>True</td>\n",
       "      <td>12TWP</td>\n",
       "      <td>Greater Yellowstone Rockies</td>\n",
       "      <td>POINT (-109.81775 43.93225)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1033_CO_SNTL</th>\n",
       "      <td>Zirkel</td>\n",
       "      <td>2846.832031</td>\n",
       "      <td>40.794880</td>\n",
       "      <td>-106.595352</td>\n",
       "      <td>Jackson</td>\n",
       "      <td>Colorado</td>\n",
       "      <td>2002-08-14 06:00:00</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>101800010202</td>\n",
       "      <td>True</td>\n",
       "      <td>13TCF</td>\n",
       "      <td>Southern Rocky Mountains</td>\n",
       "      <td>POINT (-106.59535 40.79488)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>930 rows × 13 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                          name  elevation_m   latitude   longitude  \\\n",
       "code                                                                 \n",
       "301_CA_SNTL           Adin Mtn  1886.712036  41.235828 -120.791924   \n",
       "907_UT_SNTL        Agua Canyon  2712.719971  37.522171 -112.271179   \n",
       "916_MT_SNTL         Albro Lake  2529.840088  45.597229 -111.959023   \n",
       "1267_AK_SNTL    Alexander Lake    48.768002  61.749668 -150.889664   \n",
       "908_WA_SNTL     Alpine Meadows  1066.800049  47.779572 -121.698471   \n",
       "...                        ...          ...        ...         ...   \n",
       "877_AZ_SNTL      Workman Creek  2103.120117  33.812420 -110.917732   \n",
       "1228_UT_SNTL     Wrigley Creek  2842.869629  39.132332 -111.356850   \n",
       "1197_UT_SNTL  Yankee Reservoir  2649.321533  37.747971 -112.774948   \n",
       "878_WY_SNTL        Younts Peak  2545.080078  43.932251 -109.817749   \n",
       "1033_CO_SNTL            Zirkel  2846.832031  40.794880 -106.595352   \n",
       "\n",
       "                         county       state           beginDate    endDate  \\\n",
       "code                                                                         \n",
       "301_CA_SNTL               Modoc  California 1983-10-01 00:00:00 2100-01-01   \n",
       "907_UT_SNTL                Kane        Utah 1994-10-01 00:00:00 2100-01-01   \n",
       "916_MT_SNTL             Madison     Montana 1996-09-01 00:00:00 2100-01-01   \n",
       "1267_AK_SNTL  Matanuska-Susitna      Alaska 2014-08-28 00:00:00 2100-01-01   \n",
       "908_WA_SNTL                King  Washington 1994-09-01 00:00:00 2100-01-01   \n",
       "...                         ...         ...                 ...        ...   \n",
       "877_AZ_SNTL                Gila     Arizona 1938-12-01 00:00:00 2100-01-01   \n",
       "1228_UT_SNTL            Sanpete        Utah 2012-10-01 00:00:00 2100-01-01   \n",
       "1197_UT_SNTL               Iron        Utah 2012-10-01 00:00:00 2100-01-01   \n",
       "878_WY_SNTL                Park     Wyoming 1979-10-01 00:00:00 2100-01-01   \n",
       "1033_CO_SNTL            Jackson    Colorado 2002-08-14 06:00:00 2100-01-01   \n",
       "\n",
       "                       HUC  isActive   mgrs                    mountainRange  \\\n",
       "code                                                                           \n",
       "301_CA_SNTL   180200021403      True  10TFL               Great Basin Ranges   \n",
       "907_UT_SNTL   160300020301      True  12SUG                 Colorado Plateau   \n",
       "916_MT_SNTL   100200050701      True  12TVR  Central Montana Rocky Mountains   \n",
       "1267_AK_SNTL  190205051106      True  05VPJ                              NaN   \n",
       "908_WA_SNTL   171100100501      True  10TET                    Cascade Range   \n",
       "...                    ...       ...    ...                              ...   \n",
       "877_AZ_SNTL   150601030802      True  12SWC      Southwest Basins and Ranges   \n",
       "1228_UT_SNTL  140600090303      True  12SVJ                 Colorado Plateau   \n",
       "1197_UT_SNTL  160300060202      True  12SUG                 Colorado Plateau   \n",
       "878_WY_SNTL   100800130103      True  12TWP      Greater Yellowstone Rockies   \n",
       "1033_CO_SNTL  101800010202      True  13TCF         Southern Rocky Mountains   \n",
       "\n",
       "                                 geometry  \n",
       "code                                       \n",
       "301_CA_SNTL   POINT (-120.79192 41.23583)  \n",
       "907_UT_SNTL   POINT (-112.27118 37.52217)  \n",
       "916_MT_SNTL   POINT (-111.95902 45.59723)  \n",
       "1267_AK_SNTL  POINT (-150.88966 61.74967)  \n",
       "908_WA_SNTL   POINT (-121.69847 47.77957)  \n",
       "...                                   ...  \n",
       "877_AZ_SNTL   POINT (-110.91773 33.81242)  \n",
       "1228_UT_SNTL  POINT (-111.35685 39.13233)  \n",
       "1197_UT_SNTL  POINT (-112.77495 37.74797)  \n",
       "878_WY_SNTL   POINT (-109.81775 43.93225)  \n",
       "1033_CO_SNTL  POINT (-106.59535 40.79488)  \n",
       "\n",
       "[930 rows x 13 columns]"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "snotel_stations_gdf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "0240ed52-6ecf-48ff-b3d6-0466cc6ee450",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 4.88 s, sys: 150 ms, total: 5.03 s\n",
      "Wall time: 1min 12s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "ccss_stations_gdf = snotel_ccss_stations.all_ccss_stations()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "388da4a2-2659-4918-9664-515a6c4c5cb2",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "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>name</th>\n",
       "      <th>elevation_m</th>\n",
       "      <th>latitude</th>\n",
       "      <th>longitude</th>\n",
       "      <th>county</th>\n",
       "      <th>state</th>\n",
       "      <th>HUC</th>\n",
       "      <th>mgrs</th>\n",
       "      <th>mountainRange</th>\n",
       "      <th>beginDate</th>\n",
       "      <th>endDate</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>code</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>FRW</th>\n",
       "      <td>Farewell Gap</td>\n",
       "      <td>2895.60</td>\n",
       "      <td>36.415211</td>\n",
       "      <td>-118.578979</td>\n",
       "      <td>Tulare</td>\n",
       "      <td>California</td>\n",
       "      <td>180300070202</td>\n",
       "      <td>11SLA</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2004-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-118.57898 36.41521)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>BNK</th>\n",
       "      <td>Bonanza King</td>\n",
       "      <td>1965.96</td>\n",
       "      <td>41.083118</td>\n",
       "      <td>-122.631271</td>\n",
       "      <td>Trinity</td>\n",
       "      <td>California</td>\n",
       "      <td>180102110304</td>\n",
       "      <td>10TEL</td>\n",
       "      <td>Klamath Mountains</td>\n",
       "      <td>2011-01-03</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-122.63127 41.08312)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CDP</th>\n",
       "      <td>Cedar Pass</td>\n",
       "      <td>2164.08</td>\n",
       "      <td>41.583000</td>\n",
       "      <td>-120.303000</td>\n",
       "      <td>Modoc</td>\n",
       "      <td>California</td>\n",
       "      <td>180200020603</td>\n",
       "      <td>10TGM</td>\n",
       "      <td>Great Basin Ranges</td>\n",
       "      <td>NaT</td>\n",
       "      <td>NaT</td>\n",
       "      <td>POINT (-120.30300 41.58300)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CRL</th>\n",
       "      <td>Charlotte Lake</td>\n",
       "      <td>3169.92</td>\n",
       "      <td>36.777523</td>\n",
       "      <td>-118.426010</td>\n",
       "      <td>Fresno</td>\n",
       "      <td>California</td>\n",
       "      <td>180300100206</td>\n",
       "      <td>11SLA</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2004-10-01</td>\n",
       "      <td>2023-12-11</td>\n",
       "      <td>POINT (-118.42601 36.77752)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CHM</th>\n",
       "      <td>Chilkoot Meadow</td>\n",
       "      <td>2179.32</td>\n",
       "      <td>37.408390</td>\n",
       "      <td>-119.492188</td>\n",
       "      <td>Madera</td>\n",
       "      <td>California</td>\n",
       "      <td>180400061101</td>\n",
       "      <td>11SKB</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2006-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-119.49219 37.40839)</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SLT</th>\n",
       "      <td>Slate Creek</td>\n",
       "      <td>1737.36</td>\n",
       "      <td>41.043980</td>\n",
       "      <td>-122.480103</td>\n",
       "      <td>Shasta</td>\n",
       "      <td>California</td>\n",
       "      <td>180200050304</td>\n",
       "      <td>10TEL</td>\n",
       "      <td>Klamath Mountains</td>\n",
       "      <td>2004-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-122.48010 41.04398)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SLI</th>\n",
       "      <td>Slide Canyon</td>\n",
       "      <td>2804.16</td>\n",
       "      <td>38.091234</td>\n",
       "      <td>-119.431881</td>\n",
       "      <td>Tuolumne</td>\n",
       "      <td>California</td>\n",
       "      <td>180400090501</td>\n",
       "      <td>11SKC</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2005-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-119.43188 38.09123)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SLK</th>\n",
       "      <td>South Lake</td>\n",
       "      <td>2926.08</td>\n",
       "      <td>37.175903</td>\n",
       "      <td>-118.562660</td>\n",
       "      <td>Inyo</td>\n",
       "      <td>California</td>\n",
       "      <td>180901020601</td>\n",
       "      <td>11SLB</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2004-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-118.56266 37.17590)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>STL</th>\n",
       "      <td>State Lakes</td>\n",
       "      <td>3169.92</td>\n",
       "      <td>36.926483</td>\n",
       "      <td>-118.573250</td>\n",
       "      <td>Fresno</td>\n",
       "      <td>California</td>\n",
       "      <td>180300100305</td>\n",
       "      <td>11SLA</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2018-11-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-118.57325 36.92648)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TMR</th>\n",
       "      <td>Tamarack Summit</td>\n",
       "      <td>2301.24</td>\n",
       "      <td>37.163750</td>\n",
       "      <td>-119.200531</td>\n",
       "      <td>Fresno</td>\n",
       "      <td>California</td>\n",
       "      <td>180400060903</td>\n",
       "      <td>11SLB</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2011-01-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-119.20053 37.16375)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>130 rows × 12 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                 name  elevation_m   latitude   longitude    county  \\\n",
       "code                                                                  \n",
       "FRW      Farewell Gap      2895.60  36.415211 -118.578979    Tulare   \n",
       "BNK      Bonanza King      1965.96  41.083118 -122.631271   Trinity   \n",
       "CDP        Cedar Pass      2164.08  41.583000 -120.303000     Modoc   \n",
       "CRL    Charlotte Lake      3169.92  36.777523 -118.426010    Fresno   \n",
       "CHM   Chilkoot Meadow      2179.32  37.408390 -119.492188    Madera   \n",
       "...               ...          ...        ...         ...       ...   \n",
       "SLT       Slate Creek      1737.36  41.043980 -122.480103    Shasta   \n",
       "SLI      Slide Canyon      2804.16  38.091234 -119.431881  Tuolumne   \n",
       "SLK        South Lake      2926.08  37.175903 -118.562660      Inyo   \n",
       "STL       State Lakes      3169.92  36.926483 -118.573250    Fresno   \n",
       "TMR   Tamarack Summit      2301.24  37.163750 -119.200531    Fresno   \n",
       "\n",
       "           state           HUC   mgrs       mountainRange  beginDate  \\\n",
       "code                                                                   \n",
       "FRW   California  180300070202  11SLA       Sierra Nevada 2004-10-01   \n",
       "BNK   California  180102110304  10TEL   Klamath Mountains 2011-01-03   \n",
       "CDP   California  180200020603  10TGM  Great Basin Ranges        NaT   \n",
       "CRL   California  180300100206  11SLA       Sierra Nevada 2004-10-01   \n",
       "CHM   California  180400061101  11SKB       Sierra Nevada 2006-10-01   \n",
       "...          ...           ...    ...                 ...        ...   \n",
       "SLT   California  180200050304  10TEL   Klamath Mountains 2004-10-01   \n",
       "SLI   California  180400090501  11SKC       Sierra Nevada 2005-10-01   \n",
       "SLK   California  180901020601  11SLB       Sierra Nevada 2004-10-01   \n",
       "STL   California  180300100305  11SLA       Sierra Nevada 2018-11-01   \n",
       "TMR   California  180400060903  11SLB       Sierra Nevada 2011-01-01   \n",
       "\n",
       "        endDate                     geometry  \n",
       "code                                          \n",
       "FRW  2024-03-15  POINT (-118.57898 36.41521)  \n",
       "BNK  2024-03-15  POINT (-122.63127 41.08312)  \n",
       "CDP         NaT  POINT (-120.30300 41.58300)  \n",
       "CRL  2023-12-11  POINT (-118.42601 36.77752)  \n",
       "CHM  2024-03-15  POINT (-119.49219 37.40839)  \n",
       "...         ...                          ...  \n",
       "SLT  2024-03-15  POINT (-122.48010 41.04398)  \n",
       "SLI  2024-03-15  POINT (-119.43188 38.09123)  \n",
       "SLK  2024-03-15  POINT (-118.56266 37.17590)  \n",
       "STL  2024-03-15  POINT (-118.57325 36.92648)  \n",
       "TMR  2024-03-15  POINT (-119.20053 37.16375)  \n",
       "\n",
       "[130 rows x 12 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ccss_stations_gdf"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "600b9202-ab70-4495-b7d6-390744f66078",
   "metadata": {},
   "source": [
    "## 2. Download data for all SNOTEL and CCSS stations (remember to re-run) and then save a compressed file of all CSV data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "id": "2b1781d8-db7c-4042-b9d9-ecca06135dae",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "#snotel_ccss_stations.download_snotel_data_csv(snotel_stations_gdf)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "id": "8c9dd490-91b5-4ee9-8c8a-9030d85431cc",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "#snotel_ccss_stations.download_ccss_data_csv(ccss_stations_gdf)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "id": "8adf28f4-54c8-4d2f-9e50-c58da2fd035a",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "#os.system('tar --lzma -cvf data/all_station_data.tar.lzma data/*.csv')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "cf4f9da2-4ecf-4f84-94d3-1227c2f4fb1e",
   "metadata": {},
   "source": [
    "## 3. Merge SNOTEL and CCSS station geodataframes, add new column \"csvData\" and set to True if we have data, save final geojson"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "f5808290-897e-47b9-a507-2d73b9201cf4",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "all_stations_gdf = pd.concat([snotel_stations_gdf,ccss_stations_gdf])\n",
    "all_stations_gdf\n",
    "condition = all_stations_gdf.index.str.contains(\"SNTL\")\n",
    "all_stations_gdf['network'] = np.where(condition,\"SNOTEL\",\"CCSS\")\n",
    "all_stations_gdf = all_stations_gdf[['name','network','elevation_m','latitude','longitude','state','HUC','mgrs','mountainRange','beginDate','endDate','geometry']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "e45d82fb-1ea3-4a8b-83d3-a3a632186249",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "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>name</th>\n",
       "      <th>network</th>\n",
       "      <th>elevation_m</th>\n",
       "      <th>latitude</th>\n",
       "      <th>longitude</th>\n",
       "      <th>state</th>\n",
       "      <th>HUC</th>\n",
       "      <th>mgrs</th>\n",
       "      <th>mountainRange</th>\n",
       "      <th>beginDate</th>\n",
       "      <th>endDate</th>\n",
       "      <th>geometry</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>code</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>301_CA_SNTL</th>\n",
       "      <td>Adin Mtn</td>\n",
       "      <td>SNOTEL</td>\n",
       "      <td>1886.712036</td>\n",
       "      <td>41.235828</td>\n",
       "      <td>-120.791924</td>\n",
       "      <td>California</td>\n",
       "      <td>180200021403</td>\n",
       "      <td>10TFL</td>\n",
       "      <td>Great Basin Ranges</td>\n",
       "      <td>1983-10-01</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>POINT (-120.79192 41.23583)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>907_UT_SNTL</th>\n",
       "      <td>Agua Canyon</td>\n",
       "      <td>SNOTEL</td>\n",
       "      <td>2712.719971</td>\n",
       "      <td>37.522171</td>\n",
       "      <td>-112.271179</td>\n",
       "      <td>Utah</td>\n",
       "      <td>160300020301</td>\n",
       "      <td>12SUG</td>\n",
       "      <td>Colorado Plateau</td>\n",
       "      <td>1994-10-01</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>POINT (-112.27118 37.52217)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>916_MT_SNTL</th>\n",
       "      <td>Albro Lake</td>\n",
       "      <td>SNOTEL</td>\n",
       "      <td>2529.840088</td>\n",
       "      <td>45.597229</td>\n",
       "      <td>-111.959023</td>\n",
       "      <td>Montana</td>\n",
       "      <td>100200050701</td>\n",
       "      <td>12TVR</td>\n",
       "      <td>Central Montana Rocky Mountains</td>\n",
       "      <td>1996-09-01</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>POINT (-111.95902 45.59723)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1267_AK_SNTL</th>\n",
       "      <td>Alexander Lake</td>\n",
       "      <td>SNOTEL</td>\n",
       "      <td>48.768002</td>\n",
       "      <td>61.749668</td>\n",
       "      <td>-150.889664</td>\n",
       "      <td>Alaska</td>\n",
       "      <td>190205051106</td>\n",
       "      <td>05VPJ</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2014-08-28</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>POINT (-150.88966 61.74967)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>908_WA_SNTL</th>\n",
       "      <td>Alpine Meadows</td>\n",
       "      <td>SNOTEL</td>\n",
       "      <td>1066.800049</td>\n",
       "      <td>47.779572</td>\n",
       "      <td>-121.698471</td>\n",
       "      <td>Washington</td>\n",
       "      <td>171100100501</td>\n",
       "      <td>10TET</td>\n",
       "      <td>Cascade Range</td>\n",
       "      <td>1994-09-01</td>\n",
       "      <td>2100-01-01</td>\n",
       "      <td>POINT (-121.69847 47.77957)</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",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SLT</th>\n",
       "      <td>Slate Creek</td>\n",
       "      <td>CCSS</td>\n",
       "      <td>1737.360000</td>\n",
       "      <td>41.043980</td>\n",
       "      <td>-122.480103</td>\n",
       "      <td>California</td>\n",
       "      <td>180200050304</td>\n",
       "      <td>10TEL</td>\n",
       "      <td>Klamath Mountains</td>\n",
       "      <td>2004-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-122.48010 41.04398)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SLI</th>\n",
       "      <td>Slide Canyon</td>\n",
       "      <td>CCSS</td>\n",
       "      <td>2804.160000</td>\n",
       "      <td>38.091234</td>\n",
       "      <td>-119.431881</td>\n",
       "      <td>California</td>\n",
       "      <td>180400090501</td>\n",
       "      <td>11SKC</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2005-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-119.43188 38.09123)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>SLK</th>\n",
       "      <td>South Lake</td>\n",
       "      <td>CCSS</td>\n",
       "      <td>2926.080000</td>\n",
       "      <td>37.175903</td>\n",
       "      <td>-118.562660</td>\n",
       "      <td>California</td>\n",
       "      <td>180901020601</td>\n",
       "      <td>11SLB</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2004-10-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-118.56266 37.17590)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>STL</th>\n",
       "      <td>State Lakes</td>\n",
       "      <td>CCSS</td>\n",
       "      <td>3169.920000</td>\n",
       "      <td>36.926483</td>\n",
       "      <td>-118.573250</td>\n",
       "      <td>California</td>\n",
       "      <td>180300100305</td>\n",
       "      <td>11SLA</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2018-11-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-118.57325 36.92648)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>TMR</th>\n",
       "      <td>Tamarack Summit</td>\n",
       "      <td>CCSS</td>\n",
       "      <td>2301.240000</td>\n",
       "      <td>37.163750</td>\n",
       "      <td>-119.200531</td>\n",
       "      <td>California</td>\n",
       "      <td>180400060903</td>\n",
       "      <td>11SLB</td>\n",
       "      <td>Sierra Nevada</td>\n",
       "      <td>2011-01-01</td>\n",
       "      <td>2024-03-15</td>\n",
       "      <td>POINT (-119.20053 37.16375)</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1060 rows × 12 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                         name network  elevation_m   latitude   longitude  \\\n",
       "code                                                                        \n",
       "301_CA_SNTL          Adin Mtn  SNOTEL  1886.712036  41.235828 -120.791924   \n",
       "907_UT_SNTL       Agua Canyon  SNOTEL  2712.719971  37.522171 -112.271179   \n",
       "916_MT_SNTL        Albro Lake  SNOTEL  2529.840088  45.597229 -111.959023   \n",
       "1267_AK_SNTL   Alexander Lake  SNOTEL    48.768002  61.749668 -150.889664   \n",
       "908_WA_SNTL    Alpine Meadows  SNOTEL  1066.800049  47.779572 -121.698471   \n",
       "...                       ...     ...          ...        ...         ...   \n",
       "SLT               Slate Creek    CCSS  1737.360000  41.043980 -122.480103   \n",
       "SLI              Slide Canyon    CCSS  2804.160000  38.091234 -119.431881   \n",
       "SLK                South Lake    CCSS  2926.080000  37.175903 -118.562660   \n",
       "STL               State Lakes    CCSS  3169.920000  36.926483 -118.573250   \n",
       "TMR           Tamarack Summit    CCSS  2301.240000  37.163750 -119.200531   \n",
       "\n",
       "                   state           HUC   mgrs  \\\n",
       "code                                            \n",
       "301_CA_SNTL   California  180200021403  10TFL   \n",
       "907_UT_SNTL         Utah  160300020301  12SUG   \n",
       "916_MT_SNTL      Montana  100200050701  12TVR   \n",
       "1267_AK_SNTL      Alaska  190205051106  05VPJ   \n",
       "908_WA_SNTL   Washington  171100100501  10TET   \n",
       "...                  ...           ...    ...   \n",
       "SLT           California  180200050304  10TEL   \n",
       "SLI           California  180400090501  11SKC   \n",
       "SLK           California  180901020601  11SLB   \n",
       "STL           California  180300100305  11SLA   \n",
       "TMR           California  180400060903  11SLB   \n",
       "\n",
       "                                mountainRange  beginDate    endDate  \\\n",
       "code                                                                  \n",
       "301_CA_SNTL                Great Basin Ranges 1983-10-01 2100-01-01   \n",
       "907_UT_SNTL                  Colorado Plateau 1994-10-01 2100-01-01   \n",
       "916_MT_SNTL   Central Montana Rocky Mountains 1996-09-01 2100-01-01   \n",
       "1267_AK_SNTL                              NaN 2014-08-28 2100-01-01   \n",
       "908_WA_SNTL                     Cascade Range 1994-09-01 2100-01-01   \n",
       "...                                       ...        ...        ...   \n",
       "SLT                         Klamath Mountains 2004-10-01 2024-03-15   \n",
       "SLI                             Sierra Nevada 2005-10-01 2024-03-15   \n",
       "SLK                             Sierra Nevada 2004-10-01 2024-03-15   \n",
       "STL                             Sierra Nevada 2018-11-01 2024-03-15   \n",
       "TMR                             Sierra Nevada 2011-01-01 2024-03-15   \n",
       "\n",
       "                                 geometry  \n",
       "code                                       \n",
       "301_CA_SNTL   POINT (-120.79192 41.23583)  \n",
       "907_UT_SNTL   POINT (-112.27118 37.52217)  \n",
       "916_MT_SNTL   POINT (-111.95902 45.59723)  \n",
       "1267_AK_SNTL  POINT (-150.88966 61.74967)  \n",
       "908_WA_SNTL   POINT (-121.69847 47.77957)  \n",
       "...                                   ...  \n",
       "SLT           POINT (-122.48010 41.04398)  \n",
       "SLI           POINT (-119.43188 38.09123)  \n",
       "SLK           POINT (-118.56266 37.17590)  \n",
       "STL           POINT (-118.57325 36.92648)  \n",
       "TMR           POINT (-119.20053 37.16375)  \n",
       "\n",
       "[1060 rows x 12 columns]"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "all_stations_gdf"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "f6cd9584-0d2c-4fff-9550-2542b78b9fbf",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0.5, 1.0, 'stations by network')"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAk0AAAJaCAYAAAAyKXN7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAACZoElEQVR4nOzdd3RUVdcG8GdKKoQACSm0UATpL1WaSA1VUIpKUQEVUUABRRQsgAooKGIBFEQEFUFFEAVpUhQBqdKlSK9BIIQSUmbu98f+JskkU+6dnvD81ppFcuuZm+jsnLPPPjpFURQQERERkUN6fzeAiIiIKD9g0ERERESkAoMmIiIiIhUYNBERERGpwKCJiIiISAUGTUREREQqMGgiIiIiUoFBExEREZEKDJqIiIiIVGDQRJSPzJ8/H1OnTnXrGtOnT8eXX36ZZ/uJEyeg0+ls7gsE5cqVw/333+/vZmh24MABjB07FidOnPB3U7JYftbvvfeev5tClK8waCLKR7wZNMXHx2Pz5s3o1KmTW9cnawcOHMC4ceMCKmgiItcY/d0AIgoMISEhaNSokb+bQV6kKApu377t72YQ5VvsaSIKEJcuXcLTTz+NMmXKICQkBCVKlEDTpk2xZs0aAECLFi2wbNkynDx5EjqdLutlMW7cODRs2BDFixdHkSJFULduXcyePRs51+QuV64c9u/fjw0bNmSdX65cOQD2h+c2btyI1q1bIyIiAuHh4WjSpAmWLVtmdcyXX34JnU6HdevW4dlnn0V0dDSioqLQrVs3nDt3zurYtWvXokWLFoiKikJYWBjKli2L7t2749atW6qe0+LFi1GrVi2EhoaiQoUK+Oijj7L23bhxA0WLFsXAgQPznHfixAkYDAZMnjzZ7rVzDltNmTIF5cuXR+HChdG4cWNs2bIlz/Hbt29Hly5dULx4cYSGhqJOnTr47rvvrJ7LQw89BABo2bJl1jP/8ssvMW3aNOj1eiQlJWUd//7770On02Hw4MFZ28xmM4oVK4YXX3wxa9uVK1cwaNAglCpVCsHBwahQoQJeffVVpKWlWbVPp9NhyJAh+PTTT1G1alWEhIRg7ty5Nt97RkYG+vbti8KFC+OXX36x+4yI7mgKEQWEdu3aKSVKlFBmzpyprF+/XlmyZInyxhtvKAsWLFAURVH279+vNG3aVImLi1M2b96c9bLo16+fMnv2bGX16tXK6tWrlbfeeksJCwtTxo0bl3XMzp07lQoVKih16tTJOn/nzp2KoijK8ePHFQDKnDlzso5fv369EhQUpNSrV09ZuHChsmTJEqVt27aKTqfLapeiKMqcOXMUAEqFChWU5557Tlm5cqXy+eefK8WKFVNatmyZddzx48eV0NBQJTExUVmyZImyfv165ZtvvlEee+wx5erVqw6fT0JCglKqVCmlbNmyyhdffKEsX75c6dOnjwJAmTx5ctZxw4cPVwoVKqQkJydbnf/SSy8poaGhyn///Wf3HpZnUK5cOaV9+/bKkiVLlCVLlig1a9ZUihUrZnXNtWvXKsHBwUqzZs2UhQsXKitWrFD69etn9QyTkpKUCRMmKACUadOmZT3zpKQk5Z9//lEAKPPnz8+6Zvv27ZWwsDClUqVKWdv++usvBYCyfPlyRVEUJTU1ValVq5ZSqFAh5b333lNWrVqlvP7664rRaFQ6duxo9X4AKKVKlVJq1aqlzJ8/X1m7dq2yb9++rPdpeW5Xr15VWrZsqcTFxSnbt293+HMgupMxaCIKEIULF1aGDRvm8JhOnTopCQkJTq9lMpmUjIwM5c0331SioqIUs9mcta969epK8+bN85xjK2hq1KiREhMTo1y/fj1rW2ZmplKjRg2ldOnSWde1BE2DBg2yuuakSZMUAMr58+cVRVGUH374QQGg/P33307fQ24JCQmKTqfLc25iYqJSpEgR5ebNm4qiKMq///6r6PV65YMPPsg6JjU1VYmKilL69+/v8B6WZ1CzZk0lMzMza/vWrVsVAMq3336bta1KlSpKnTp1lIyMDKtr3H///Up8fLxiMpkURVGU77//XgGgrFu3Ls/9SpcurTzxxBOKoihKWlqaUqhQIeXll19WACgnT55UFEVRxo8frwQFBSk3btxQFEVRPv30UwWA8t1331ld691331UAKKtWrcraBkCJjIxUrly5YvN9Tp48WTl+/LhSrVo1pVq1asqJEyccPh+iOx2H54gCxD333IMvv/wSb7/9NrZs2YKMjAxN569duxZt2rRBZGQkDAYDgoKC8MYbb+Dy5ctWQ0Bq3bx5E3/99Rd69OiBwoULZ203GAx47LHHcObMGRw6dMjqnC5dulh9X6tWLQDAyZMnAQC1a9dGcHAwnn76acydOxfHjh3T1Kbq1avjf//7n9W23r17IyUlBTt37gQAVKhQAffffz+mT5+eNTQ5f/58XL58GUOGDFF1n06dOsFgMNh9H0ePHsU///yDPn36AAAyMzOzXh07dsT58+fzPBtbWrdunTX8umnTJty6dQsvvPACoqOjsXr1agDAmjVr0LhxYxQqVAiA/JwLFSqEHj16WF2rX79+AIDffvvNanurVq1QrFgxm/ffuXMnGjVqhNjYWPz5559ISEhw2maiOxmDJqIAsXDhQvTt2xeff/45GjdujOLFi+Pxxx/HhQsXnJ67detWtG3bFgAwa9Ys/Pnnn9i2bRteffVVAEBqaqrm9ly9ehWKoiA+Pj7PvpIlSwIALl++bLU9KirK6vuQkBCr+1esWBFr1qxBTEwMBg8ejIoVK6JixYr48MMPVbUpLi7O7racbRk6dCiOHDmSFXhMmzYNjRs3Rt26dVXdx9n7uHjxIgBgxIgRCAoKsnoNGjQIAPDff/85vU+bNm1w6tQpHDlyBGvWrEGdOnUQExODVq1aYc2aNUhNTcWmTZvQpk2brHMuX76MuLg4q3w2AIiJiYHRaMzzM7H187NYvXo1Ll68iKeeegpFixZ12l6iOx1nzxEFiOjoaEydOhVTp07FqVOnsHTpUrzyyitISkrCihUrHJ67YMECBAUF4ZdffkFoaGjW9iVLlrjcnmLFikGv1+P8+fN59lmSu6OjozVft1mzZmjWrBlMJhO2b9+Ojz/+GMOGDUNsbCx69uzp8FxbAaRlW85Ap1WrVqhRowY++eQTFC5cGDt37sTXX3+tua32WN73qFGj0K1bN5vH3H333U6v07p1awDSm7R69WokJiZmbX/ttdfw+++/Iy0tzSpoioqKwl9//QVFUawCp6SkJGRmZub5meQOrnJ66aWX8O+//+Lxxx9HZmYmHn/8cadtJrqTsaeJKACVLVsWQ4YMQWJiYtawEyA9HrZ6jXQ6HYxGo9WQUmpqKr766qs8x9q7Rm6FChVCw4YN8eOPP1odbzab8fXXX6N06dKoXLmy1reWxWAwoGHDhpg2bRoAWL1Pe/bv34/du3dbbZs/fz4iIiLy9CI9//zzWLZsGUaNGoXY2NisWWyecPfdd6NSpUrYvXs36tevb/MVEREBIG8vVU7x8fGoVq0aFi1ahB07dmQFTYmJibh06RKmTJmCIkWKoEGDBlnntG7dGjdu3MgTEM+bNy9rv1p6vR6fffYZhg4din79+mHGjBmangPRnYY9TUQB4Nq1a2jZsiV69+6NKlWqICIiAtu2bcOKFSusejJq1qyJH3/8ETNmzEC9evWg1+tRv359dOrUCVOmTEHv3r3x9NNP4/Lly3jvvfeyPrBzqlmzJhYsWICFCxeiQoUKCA0NRc2aNW22a+LEiUhMTETLli0xYsQIBAcHY/r06di3bx++/fZbh70Ytnz66adYu3YtOnXqhLJly+L27dv44osvAMCqN8WekiVLokuXLhg7dizi4+Px9ddfY/Xq1Xj33XcRHh5udeyjjz6KUaNG4ffff8drr72G4OBgTW115rPPPkOHDh3Qrl079OvXD6VKlcKVK1dw8OBB7Ny5E99//z0AoEaNGgCAmTNnIiIiAqGhoShfvnxWz1jr1q3x8ccfIywsDE2bNgUAlC9fHuXLl8eqVavQpUsXGI3Z/6t+/PHHMW3aNPTt2xcnTpxAzZo1sXHjRkyYMAEdO3ZU9Rxze//99xEREYFBgwbhxo0beOmll9x9PEQFk58T0YlIUZTbt28rzzzzjFKrVi2lSJEiSlhYmHL33XcrY8aMyZoVpiiKcuXKFaVHjx5K0aJFFZ1Op+T8T/iLL75Q7r77biUkJESpUKGCMnHiRGX27NkKAOX48eNZx504cUJp27atEhERoQDImo1na/acoijKH3/8obRq1UopVKiQEhYWpjRq1Ej5+eefrY6xzJ7btm2b1fZ169ZZzRzbvHmz0rVrVyUhIUEJCQlRoqKilObNmytLly51+owSEhKUTp06KT/88INSvXp1JTg4WClXrpwyZcoUu+f069dPMRqNypkzZ5xeP+czyFnCwAKAMmbMGKttu3fvVh5++GElJiZGCQoKUuLi4pRWrVopn376qdVxU6dOVcqXL68YDIY8z/inn35SACiJiYlW5wwYMEABoHz00Ud52nL58mXlmWeeUeLj4xWj0agkJCQoo0aNUm7fvp2nzYMHD1b9PidPnqwAUN544w2bz4foTqdTlByV74iICoj09HSUK1cO9957r1XBSSIiV3F4jogKlEuXLuHQoUOYM2cOLl68iFdeecXfTSKiAoJBExEVKMuWLUP//v0RHx+P6dOnqy4zQETkDIfniIiIiFRgyQEiIiIiFRg0EREREanAoImIiIhIhYBLBDebzTh37hwiIiI0F84jIiIi0kpRFFy/fh0lS5aEXm+/PynggqZz586hTJky/m4GERER3WFOnz6N0qVL290fcEGTZb2m06dPo0iRIn5uDRERERV0KSkpKFOmTFYMYk/ABU2WIbkiRYowaCIiIiKfcZYWxERwIiIiIhUYNBERERGpwKCJiIiISIWAy2kiIiLKz0wmEzIyMvzdDMohKCgIBoPB7eswaCIiIvIARVFw4cIFJCcn+7spZEPRokURFxfnVg1IBk1EREQeYAmYYmJiEB4ezgLNAUJRFNy6dQtJSUkAgPj4eJevxaCJiIjITSaTKStgioqK8ndzKJewsDAAQFJSEmJiYlweqmMiOBERkZssOUzh4eF+bgnZY/nZuJNvxqCJiIjIQzgkF7g88bPh8BwREVGAOHsWWL0aSEkBihQBEhOBUqX83SqyYNBERETkZ9u3AxMmAD//DGRmZm83GoHOnYHRo4H69f3XPhIcniMiIvKjxYuBe++Vf3MGTIB8b9m/ZIn32nDhwgU899xzqFChAkJCQlCmTBl07twZv/32W9Yxu3btwkMPPYTY2FiEhoaicuXKGDBgAA4fPpx1zKJFi9CwYUNERkYiIiIC1atXx4svvpi132QyYeLEiahSpQrCwsJQvHhxNGrUCHPmzPHem/MgBk1ERER+sn070KsXkJbm+Li0NKBnTzne006cOIF69eph7dq1mDRpEvbu3YsVK1agZcuWGDx4MADgl19+QaNGjZCWloZvvvkGBw8exFdffYXIyEi8/vrrAIA1a9agZ8+e6NGjB7Zu3YodO3Zg/PjxSE9Pz7rX2LFjMXXqVLz11ls4cOAA1q1bhwEDBuDq1auef2NeoFMURfF3I3JKSUlBZGQkrl27hiJFivi7OURERE7dvn0bx48fR/ny5REaGqr6vG7dpCdJy/GLFrnQQAc6duyIPXv24NChQyhUqJDVvuTkZAQHByMhIQH33nsvFttobHJyMooWLYphw4Zh9+7dWLdund171a5dG127dsWYMWM8+yZUcPQzUht7sKeJiIjID86eBZYu1XbO0qVynqdcuXIFK1aswODBg/METIBU0V65ciX+++8/jBw50uY1ihYtCgCIi4vD/v37sW/fPrv3i4uLw9q1a3Hp0iWPtN/XGDQRERH5werVgMmk7ZzMTGDNGs+14ejRo1AUBVWqVLF7zJEjRwDA4TEA8Nxzz6FBgwaoWbMmypUrh549e+KLL75AWo6xxylTpuDSpUuIi4tDrVq18Mwzz+DXX3/1zJvxAQZNdGe5fh1YuVKmqJw54+/WENEdLCXFt+fZYsnQcVTDSG0WT6FChbBs2TIcPXoUr732GgoXLowXX3wR99xzD27dugUAqFatGvbt24ctW7agf//+uHjxIjp37oynnnrK/TfjAwya6M6QkQG8/DIQFwe0bw906QKULQt07QpcvOjv1hHRHcjVtF1PpvtWqlQJOp0OBw8etHtM5cqVAQD//POPqmtWrFgRTz31FD7//HPs3LkTBw4cwMKFC7P26/V6NGjQAMOHD8fixYvx5ZdfYvbs2Th+/Lh7b8YHGDRRwacowOOPA5MnA///107W9p9/Bpo2BfLJzA0iKjgSE6UOkxZGI9CmjefaULx4cbRr1w7Tpk3DzZs38+xPTk5G27ZtER0djUmTJtm8RnJyst3rlytXDuHh4TavbVGtWjUAcHhMoGDQRAXfn38CCxZIkJSbyQQcPw588onv20VEd7RSpaRwpRZduni+Qvj06dNhMplwzz33YNGiRThy5AgOHjyIjz76CI0bN0ahQoXw+eefY9myZejSpQvWrFmDEydOYPv27Rg5ciSeeeYZAFJOYOTIkVi/fj2OHz+OXbt24YknnkBGRgYSExMBAD169MAHH3yAv/76CydPnsT69esxePBgVK5c2WnOVCBg0EQF35w5jv+cM5uBmTN91x4iov83ejQQEqLu2NBQYNQoz7ehfPny2LlzJ1q2bIkXX3wRNWrUQGJiIn777TfMmDEDAPDAAw9g06ZNCAoKQu/evVGlShX06tUL165dw9tvvw0AaN68OY4dO4bHH38cVapUQYcOHXDhwgWsWrUKd999NwCgXbt2+Pnnn9G5c2dUrlwZffv2RZUqVbBq1SoYtXa7+YGmOk3lypXDyZMn82wfNGgQpk2bBkVRMG7cOMycORNXr15Fw4YNMW3aNFSvXl11g1iniTyuXTtg1SrHxxiNkvdEROQCV+s0AVLpu2dPxwUuQ0Kkw/zBB91q5h3N53Watm3bhvPnz2e9Vq9eDQB46KGHAACTJk3ClClT8Mknn2Dbtm2Ii4tDYmIirl+/rvW9EXlOXJzzxIHoaN+0hYgolwcfBDZulMKVuf9XZTTK9o0bGTAFAk19YSVKlLD6/p133kHFihXRvHlzKIqCqVOn4tVXX0W3bt0AAHPnzkVsbCzmz5+PgQMHeq7VRFo89hgwb579/QYD8MQTvmsPEVEu9etLpe+zZ6UOU0qKzJJr08bzOUzkOpcHENPT0/H111/jhRdegE6nw7Fjx3DhwgW0bds265iQkBA0b94cmzZtshs0paWlWRW+SvFkAQoiAGjdWsoMrFol+Us5GY1AiRLA88/7p21ERDmUKgX07evvVpA9LieCL1myBMnJyejXrx8AWSEZAGJjY62Oi42Nzdpny8SJExEZGZn1KlOmjKtNIrJNpwN+/BHo1y9v33ejRjK7LtfvLRERUW4uB02zZ89Ghw4dULJkSavtuauKKorisNLoqFGjcO3atazX6dOnXW0SkX1hYcDs2cDp08DXXwNffAHs2QP88QdQvry/W0dERPmAS8NzJ0+exJo1a/Djjz9mbYuLiwMgPU7x8fFZ25OSkvL0PuUUEhKCELXzLYnccesWcPgwEB8P1KrF5G8iItLEpZ6mOXPmICYmBp06dcraVr58ecTFxWXNqAMk72nDhg1o0qSJ+y0lclVmJvDGGxIsNW8uOU7x8ZI4wErgRESkkuaeJrPZjDlz5qBv375Whah0Oh2GDRuGCRMmoFKlSqhUqRImTJiA8PBw9O7d26ONJlLNsoRK7orgmZnAN98AO3cCmzcDhQv7r41ERJQvaA6a1qxZg1OnTuEJG1O0R44cidTUVAwaNCiruOWqVasQERHhkcYSafbHH8C339reZzIB+/cDn30GvPiib9tFRGTL2bPA6tXZNQcSE1lzIIBoqgjuC6wITh7Vr5/0KGVm2j+mUiXJdSIicpE7FcEBANu3AxMmyCLiOf9/ZTTKAnWjR0sxJ3KZzyuCE+U7x487DpgAmVFHROQvixcD994r/+b+/1VmZvb+JUu8cvukpCQMHDgQZcuWRUhICOLi4tCuXTts3rwZgCyhptPpsGXLFqvzhg0bhhYtWlhtu3LlCoYNG4Zy5cohODgY8fHx6N+/P06dOpV1jE6nc/iylDKyt3/BggUAgPXr10On0yE5Odkrz8WWwF8dj8gdMTFS8dtksn9MVJTv2kNElNP27UCvXo4XngNkf8+esp6Kh3ucunfvjoyMDMydOxcVKlTAxYsX8dtvv+HKlStZx4SGhuLll1/Ghg0b7F7nypUraNSoEYKDgzF9+nTUqFEDJ06cwGuvvYYGDRpg8+bNqFChAs6fP591zsKFC/HGG2/g0KFDWdvCwsKyvp4zZw7at29vdZ+iRYt64F27hkETFWyPPQb88IP9/QYD0L+/79pDRJTThAnOAyaLtDRg4kRZb8VDkpOTsXHjRqxfvx7NmzcHACQkJOCee+6xOm7gwIGYMWMGli9fjo4dO9q81quvvopz587h6NGjWWWIypYti5UrV6JSpUoYPHgwfv3116x9ABAZGQmdTme1LaeiRYva3ecPHJ6jgq1TJ6BpUwmOcjMapVbTkCG+bxcR0dmzwNKl2s5ZulTO85DChQujcOHCWLJkidWSZrmVK1cOzzzzDEaNGgVz7uWoIDPrFyxYgD59+uQJcsLCwjBo0CCsXLnSqvcqP2LQRAWbwQAsXw507SrLqQDZ/9auLV3dXEKFiPxh9WrHqQO2ZGbKir4eYjQa8eWXX2Lu3LkoWrQomjZtitGjR2PPnj15jn3ttddw/PhxfPPNN3n2Xbp0CcnJyahatarN+1StWhWKouDo0aOa2terV6+swM7yOnbsmKZreBKDJir4ihQBvv8e+PdfKS/wySfAtm3yuusuf7eOiO5Uri5Q7+GF7bt3745z585h6dKlaNeuHdavX4+6deviyy+/tDquRIkSGDFiBN544w2kp6druodlor6jZdVs+eCDD/D3339bvfy5Ri1zmqjguX0b+O47YNkyID0dqFMHeOopWWPu6af93ToiIuFqWR0vlOMJDQ1FYmIiEhMT8cYbb+Cpp57CmDFjsmayWbzwwguYPn06pk+fbrW9RIkSKFq0KA4cOGDz+v/88w90Oh0qVqyoqV1xcXG4K4D+uGVPExUshw5J3aW+fSUBfMkSYNw4ICEB+Oorf7eOiChbYqLkVmphNAJt2ninPTlUq1YNN2/ezLO9cOHCeP311zF+/Hik5Ojx0uv1ePjhhzF//nxcuHDB6pzU1FRMnz4d7dq1Q/Hixb3edm9i0EQFx+3b8j8hy3RWS7Ki2Sx5AH37An/+6b/2ERHlVKqUFK7UoksXj1YIv3z5Mlq1aoWvv/4ae/bswfHjx/H9999j0qRJeOCBB2ye8/TTTyMyMhLf5lptYfz48YiLi0NiYiJ+/fVXnD59Gr///jvatWuHjIwMTJs2TXP7kpOTceHCBatX7mBu7969eYbwvIXDc1RwfP+940KVBgPw3nsym46IKBCMHi2TVdSUHQgNBUaN8ujtCxcujIYNG+KDDz7Av//+i4yMDJQpUwYDBgzA6NGjbZ4TFBSEt956K8+6stHR0diyZQvefPNNDBw4EOfPn0dUVBTat2+Pr7/+GmXLltXcvv42SsJMnDgRr7zyStb39913X55jvLXYCZdRoYKjZ08JnGxMh80SFCT/c9KYjEhE5Ihby6gsWSL//3IUOIWEyMLjDz7oTjPvaFxGhSin9HTA2d8AmZnOjyEi8qUHH5TyJ9265c1xMhpl+8aNDJgCAIfnqOCoWxf46Sf7QZFeD1SvLv8SEQWS+vWl0vfZs1KHKSVFZsm1aePRHCZyD4MmKjjuvddxL5LZDAwd6rv2EBFpVaqUTFqhgMQ/ualgOHzYedf1ww9znTkiInIZe5q0MJul23TuXOlCLVMG6NcPaNWKicX+9sYbwM2bjofmPviAQ3NE5FUBNreKcvDEz4afIGqlpgIdOwLt2gELFwIbNshMhjZtgAceUL9KNXleSorkAmRmOj4uV00RIiJPCQoKAgDcunXLzy0heyw/G8vPyhXsaVJr6FBZXBHIXmDR8iG9bBnw0kvARx/5p213ukuXnAdMBgNw5oxv2kNEdxyDwYCiRYsiKSkJABAeHq55nTXyDkVRcOvWLSQlJaFo0aIwGAwuX4tBkxpJScCcOfbr/5jNwMyZwJtvAkWL+rRpBCAqSobdHNVnMpuB2FjftYmI7jhxcXEAkBU4UWApWrRo1s/IVQya1Fi/3nlPRloa8Mcf2kvik/uKFgXuv196/Cy9gLmZzUCu6rVERJ6k0+kQHx+PmJgYZGRk+Ls5lENQUJBbPUwWDJrUcBYwWfA/Ev956y0ZPk1Ptx04DR8OuFDCn4hIK4PB4JEPaAo8TARXo35958fodEC9et5vC9lWqxawbh1QqZL19rAwmVk3ebJ/2kVERAUGe5rUqFwZaN1aZszZ6nUyGoH27YGEBN+3jbI1bAgcOABs3gwcOgRERMhsx4gIf7eMiIgKAC7Yq9aZM0DTpvJvzoRjvR4oVw7480/AzQQzIiIi8j0u2OtppUsDO3cCY8ZIUcvQUOlZevNNYPt2Bky+dvAg8NxzspZcjRrACy8AR4/6u1VERFSAsaeJ8p9582Q5FL0+e7jUYJC8sgULgO7d/ds+IiLKV9jTRAXTnj0SMJnN1vllJpO8evYE/v3Xf+0jIqICi0ET5S8ff2x//ThFkdf06b5tExER3RE4e44Cy8GDwO+/y1BbixYyczGn1asd180ymbKXuyEiIvIgBk0UGM6fBx59FFi71np7u3bAV18BJUrI92pS8AIrTY+IiAoIDs+R/12/DjRvLnWwcvvtN+lxSk2VXqR69STp2x6DAWjVymtNJSKiOxd7msj/5s6VcgG2eogyM6VgZb9+srbf+fPOrzdokMebSERExKCJ/O/LL50f8913jvcbjTKjbu5c4O67PdIsIiKinDg8R/6XlOReHlJMDPDkk8Du3UCfPp5rlzuOHQNefBGoWBEoVQp48EFgzRp/t4qIiNzAnibyv4QE4OxZ6+Vp1DIYpCr4p596vl2uWrMG6NwZyMiQPCwAWLYM+Okn4KWXgHffldmBRESUr7CnifxvwADXAiZAgpJDhzzbHndcvSq9Smlp2QETkF0mYfJkYNEivzSNiIjcw6CJ/K9nT1kM2dGsOEciIz3bHnfMnQvcumV/uFGvBz74wLdtIiIij2DQRP4XHAysXAk89RQQEpK9PTRUShE4GsrS6wMnjwkANm50vN9sBrZscb1njYiI/IZBEwWGQoUkL+nCBckJ+u034OJFYOlSoHRpmR2Xm8EAREUBTz/t+/bao9Opy1diThMRUb7DoIkCS9GiQOvWUqCySBF5bdgAVKki+43G7ACqfHnZZ6kWHghatXI8E9BgcN57RkREAUmnKIG15kRKSgoiIyNx7do1FClSxN/NoUChKMD69cC6dfJ1s2ZAmzb2F+/1l+vXZTbgtWv2h+B++QXo1Mm37SIiIrvUxh4sOUD5g04HtGwpr0AWESHlBdq3B27cyA6cjEaZQTdhAgMmIqJ8ikETkac1bgwcPgx8/jmweDFw+zbQsKEs71Kvnr9bR0RELuLwHBEREd3R1MYeAZYQQkRERBSYGDQRERERqcCgiYiIiEgFJoITedr588AXXwC7dwNhYUCXLvIKCvJ3y4iIyA0Mmoi0unwZSE0F4uLyViqfO1eWg7GUGtDpgHnzgEqVgNWrpYYTERHlSxyeI1Lrl1+ARo2A6GigTBkgNhZ47TXg5k3Zv3490L+/1GMym+VlMsm+48eBtm1lHxER5UsMmojU+PRToHNnYNu27G1XrkixygoVgOeeA15+2f7yKJmZUrvpl198014iIvI4Bk1Ezpw7BwwZIl/nXhpFUYCkJGD6dGDrVvtLpwAylPfzz95rJxEReRWDJiJn5sxxvAgv4DhYsjCZgD//lFdg1ZQlIiIVGDTZk5wsPQd79mTnpdCd6eef1QVFzigKcOQIcO+9stjwtWvuX5OIiHxGc9B09uxZPProo4iKikJ4eDhq166NHTt2ZO3v168fdDqd1atRo0YebbRXXb4MPPGEJPk2bAj8739A+fKS08LegTvPwYPAX3957nqW4GvDBuCRRzx3XSIi8jpNJQeuXr2Kpk2bomXLlvj1118RExODf//9F0WLFrU6rn379pgzZ07W98HBwR5prNclJwNNmwJHj1r3Lp0+DTz7LHDmDPD2235rHnlReroEMleuABUrysK6Oh0wYwag13umpyknkwlYuRLYtQuoU8ez1yYiIq/QFDS9++67KFOmjFVAVK5cuTzHhYSEIC4uzu3G+dz77+cNmHIaPx7o21dq7lDB8dlnUjrgv/+yt1WvLtv/+MPzAZOF0QgsWeK9oOn2bWDRIimyGRoqs/8aNPDOvYiI7gCahueWLl2K+vXr46GHHkJMTAzq1KmDWbNm5Tlu/fr1iImJQeXKlTFgwAAkJSV5rMFeoygyBOcof8lgkErPVHBMnQo884x1wATIsFyrVkBamvfurdMBt25559qrVwOlSgGPPirvceJE4J57gBYtZAiaiIg00xQ0HTt2DDNmzEClSpWwcuVKPPPMM3j++ecxb968rGM6dOiAb775BmvXrsX777+Pbdu2oVWrVkiz8+GTlpaGlJQUq5df3L6d94MzN0UBjh3zTXvI+1JSgFGjbO+zFKa8dUuCZW/IyJAeLU/bvRu4/34Zbrbcx1JUc+NGoFMn7/WeEREVYDpFUZ/dHBwcjPr162PTpk1Z255//nls27YNmzdvtnnO+fPnkZCQgAULFqBbt2559o8dOxbjxo3Ls/3atWsoUqSI2qapYzYDq1YBv/4qOSz16wM9ewKFCsm+sDDZbo/RKEtkzJjh2XaRf8yZAzz5pPME/9BQ6XHy5EQAnQ6IiJB16sLDPXddAOjdG/j+e8fVx1etAhITPXtfIqJ8KiUlBZGRkU5jD009TfHx8ahWrZrVtqpVq+LUqVMOz0lISMCRI0ds7h81ahSuXbuW9Tp9+rSWJql36hRQsybQoYMUIpw9WwKgUqWANWtkfTBnCeuZmUCvXt5pH/ne+fPqepEmTZKAWu+hCh0Gg7y+/trzAZPZDPzwg+OAyWiUoIqIiDTRlAjetGlTHDp0yGrb4cOHkeBgEdLLly/j9OnTiI+Pt7k/JCQEISEhWpqhXVoa0Lo1cOKEfJ/zAyUlBWjf3nktJoMBaNkSaNbMa80kH4uLU1eDq0sX4OGHgc8/l5w2d4Zo9XoJ3F99Vdax87T0dBmOc8Rslt97IiLSRNOfzsOHD8eWLVswYcIEHD16FPPnz8fMmTMxePBgAMCNGzcwYsQIbN68GSdOnMD69evRuXNnREdHo2vXrl55A6osWiSz4mz99a0o6j44H3wQWLzY/tpilP907w44CtgNBgmSExKkbterr0o1b6OGvzUsvy/DhwOXLgHXr0uxTG/VLgsJAUqWdN6mu+/2zv2JiAowTUFTgwYNsHjxYnz77beoUaMG3nrrLUydOhV9+vQBABgMBuzduxcPPPAAKleujL59+6Jy5crYvHkzIiIivPIGVFm82L2hlYgIYOFCoHBhz7WJ/C8yEnjrLdv79Hp5vfuu9fa4OAmebNHpgKAgGfK1qFpVcqfefx+Ijvb8cJytNgwa5Pj3XVGkgCsREWmiKRHcF9QmY2nSoQOwYoV717hyBShWzDPtocChKMC0acAbbwBXr2Zvr1wZmDULuO8+2+e8/74EXDmHuWrWlFy5+vWlV0mvB6KifN87efOmtPvvv61nyVmKdE6eDIwY4ds2EREFMLWxx50RNL30ktSqcZQc64jRKB9E+aWyOWmXlgasXSvBcYUKMnzmLNhJTZVzUlIkyKpbNzCGb00mCegmTZI2WkRHAx98ILWbiIgoC4OmnA4fBqpUcW3KuNEouS8LFnimLUTepCjA448D33yT9/ddr5filitWyDAiEREB8FLJgXyrcmXgnXfk69y5Hjqd/AVuKwdEr5dk4NGjvd9GKjgURdaxe/55YMAAmYjgq79N1q6VUga27mc2y/5vvvFNW4iICpg7I2gCgJEj5cOrXr3sbaVKyfISR44APXrINr0+e3ZUbKwsqlqrlu/bS/nTP/8AZctKj87HH0uZgh49gCJFgKVLvX//mTMdz+7T62W5ICIi0uzOGJ7L7epVqWdTooR1D9PRo/LBlpoqSb0dO2qbXk53tvPnZSr/9ev2j/ntN1nTzltq15ZlVByJjpZEdSIiAqA+9rgzIwJ7s+Duugt44QXftoUKjg8/dBwwAUD//lJk1VsJ45ahZkdryxUt6p17ExEVcHfO8ByRt33+ufNjTp0C9uzxXht693YcMOn1kihORESaMWiiwHfiBLBzZ+APKaldmuTMGe+1oVcvmfhga1jZaJQh6Wee8d79iYgKMAZNFLhWrQIaNADKl5cE/rg4oFs3yT0LRHFx6o6LifFeG8LCgHXr5LkB2YsDA0ClSjKrr0QJ792fiKgAuzMTwSnwLVoEPPSQ5P7kHG4yGGQm2tatkoMWSKZMAV580fExpUoBp0/7pgjmtm2SeG4yAU2bAs2bB0bxTSKiAMPilpR/3b4ti84mJ9uuN2QwyMxGb0zhT08HLl+WwKxQIW3npqbKWnMnT9o/ZskS4IEH3GoiERF5FotbUv71009SFsJePG8yAb/8Aly44Ll7XrgADB4sMytLlpSg6cEHJZdKrbAwYO9e2+vVFS4siz4zYCIiyrfuzJIDFNj+/VeSlh2tFago0qOjNo/IkbNnZa25Cxey72k2S2D266+y7EjLls6vc/GizKAzGIDGjYHwcOl5SkxkzS8iogKA/xenwFO0qPQmqTnOE154QQpT5r6nySTB2aOPSqkAS0K1LatXS8/U7dvZOVh6PfDHH0CbNgyYiIgKAA7PUeDp1s1xgKLTATVqyNR6d126JEnn9oI0sxk4d056nOw5fVqG3XIGTJZzMzJkGZVDh9xvKxER+RWDJgo8cXHAkCH2Z3opCjB+vGdmgh096rxXy2AA9u+3v//TTyWB3FZRSUte1rRprreRiIgCAoMmCkzvvQc8/7wELHo9EBQk2wsXBubNA7p08cx91MyQM5sdH7dsmePAKzNT8qOIiChfY6IFBSaDAZg6FRg5EvjhB+DKFaBCBaB7d+2lABypUQMoV06qjjviaNabo4R1LccQEVFAY9BEga1kSelx8ha9HhgzRhbStbf/8ceBMmXsX+PeeyVnyV5gZDRKcUkiIsrXODxH+UN6OrBlC/D771J80pP69QPefdd6KNAy261HD2DGDMfnDx7sfHjuuec81lwiIvIPBk0U2MxmCWhKlpTaR82bA/Hx0vvjyeBp5EiZBTdhgvQ6vfQSsGePFKQMDXV8bs2awMcfy9c5SwtYvn7nHaBJE8+1lYiI/ILLqFBge/ZZmZ2Wm8EgC9D+9ZdU7w4Ef/4JfPABsGaNzJpr3hwYPlxdYUwiIvIbrj1H+d/OnUC9evb36/XAuHHAa6/5rk1ERFTgcO05yv+++MJxJW2zGZg503ftISKiOxqDJgpcJ044n6p/9qz9hX2JiIg8iEETBa6oKOdrtkVGeqYyOBERkRMMmihw9enjuKfJYABq1QI6dABatQJGjQKOH/dd+4iI6I7CoIkCV5s2MgPN1uK9er3kNG3YAKxcCaxbB0yeDNx1FzBnju/b6i9HjwLDhgFlywIxMUDHjrKsC4csiYg8jkETBS69Hvj5ZwkEdDp56XP8ylqG5SwBgskkgdSTTwKbN/u+vb62Zo3UiJo2TWpMXboErFoF3H+/LHjMwImIyKMYNFHg2rFD1pr7+WcJABQFqFoVeOwxCY7MZtvnGQzAlCm+bauvJScDDz4IpKVZD2FaKpNPnw7Mn++PlhERFVhce44C08aNMjyXO6fpn3/kpdPZ70nJzASWL5egKzlZhuwSErzeZJ/66ivg1i37z0Cvl8CxTx/ftouIqABjTxMFHkWRpUwyMvKu6WYyycvZ0NOtW0D9+hJ4lSsHJCZKsFVQbNrkeNag2SzFQZ2VbCAiItUYNFHg2bhREpztDb+5Yt06WbvuyBHPXdOf9HrnpRYseWBEROQRDJoo8HgjsDGZgOvXgdGjPX9tf2jdOm8vXE4GA9Csme2Zh0RE5BIGTRR4IiPVH5uz+KXeya+zyQQsXgxcueJauwJJz55AdLT9oMhkAl56ybdtIiIq4Bg0UeBp1855JXCL7t0lyAoPB0qXdt6zYjIB586530Z/Cw+X+lRFilgHi5bnNmGClB4gIiKP4ew5CjyFC0ul7507nR9rNgPp6cDt20Bqqro8qNBQ99voD2az5Gb9/TcQFiZB0ZEjUsxz8WJ5//XrA4MGAbVr+7u1REQFjk5RAqsCXkpKCiIjI3Ht2jUUKVLE380hf5k2TQo0OmM0Zs8Qc1SGIKe77pJK4iVLutdGX9q2TYbkjh2TniXL++zdG5g5U3qeiIjIJWpjDwZNFJhSUiSouXnT89c2GoFGjYA//vD8tffskQDHaJRyB6VKuX/NQ4eAevWkNy138rdeLxXTly7lTDkiIhepjT2Y00SBqUgR4LvvgKAgbcneamRmSlmDv/92/1oWx44BTZoA//sf8NRTQL9+sh7co4+6H/i9845U/rY1W85sBn75BfjrL/fuQURETjFoosDVsSOwfbsEHoULSwBVsyZw333qE8Xt0euB9etdO/fGDeD994FKlSS3KC5OgqWtW62PM5uBBQuAzp1drzllMslyKI6KVBqNwDffuHZ9IiJSjUETuWf/fmDECODhhyUBedMmzy4UW6uWJDpfvy4J33//LQGKJ7jSzqtXpUjmyJFSgPP2beDiRQmkbPUEmUySvL16tWttTE2V9+2IogCXL7t2fSIiUo1BE7nGbAaefx6oUQP48ENg0SJg1iygaVPggQckmPCWTp3cXx7EbAaaN9d+3vDhwMGD2nqODAZg3jzt9wKAQoWAYsWcH1eunGvXJyIi1Rg0kWveew/4+GP5OjNTgghLILNsmfQ6eUtiogzTuTNE16gRULeu9ORs3gz8+af0Zjly+bIMgzmqxG2LyQRcuOBaO3U6YMAAx/WnTCbgiSdcuz4REanGoIm0S0sD3n3X/n6zGZg71ztFJBVFhuiee869kgEDBsiSKnFxksB9771AbKxc99Yt2+fs2eNaD5fRKEnhrho5Us63FyS+8oqUUSAiIq9i0ETabd3qfCkSsxlYscKz9924UXKc6tUDnn4aOHUKCAnRfh29Hnj9dZmVlpKSvT01FZgxQyqS28ojcrVnKzPTvZ6gqCjpDXvoIes2xMdLb9+ECa5fm4iIVGPQRNqlpTk/RqfzbF7Tpk1Aq1bAgQPa25Kb2Sy9YLYSwU0mCc5szUarX1/buniAPIeHH5aeLHfExsosunPnpL7Utm0SNA4ZwvpMREQ+wqCJtKte3Xm9JEXx7FIeL7wgAY2rU/e10OuBzz7Luz0sDBg6VH2QEhIii+Z+/bXnApsSJSQAq1/f/bILRESkCYMm0i4+Huja1f6HtsEggVXjxp653+HDUrzRFwETIPc5ccL2vtdfl6VLgOz3byuA1OmkFyw5mcENEVEBwaCJXPPxx7JESO5ZXUajFKKcP99zvStnz3rmOlpERdnebjQCX30lQ3iPPQbcc4/tYM4y9DdzpiymS0RE+R6DJnJNfLxU637ppew6QuHhsoTIzp2SsO0pcXGeu5Yaej3Qt6/9/Tqd1KP64gugWTPH5QAMhuzSDERElK9xwV5yn6LIUFRIiPeSkuvUkSn/nhqi0+ulrblrLhmNEhDu3q2uqOR99zlf+LdwYec1oIiIyG+4YC/5jk4HhIZ6dxbXe+/J9T1xD4MBaNFCXoBc05KXVLcu8Pvv6gImQN37Dg52taVERBRAmKFK+UPr1sAvvwADB8pUewudTvsaciaTDCu2by9r561dKz1YTZvKrDQtOncG1qyxv99olGVliIgo3+PwHOUvZjOwYQNw8qQkaz/zjPbK46NHA+PHe6Y9164BFSvKLLncQ306nfRq7djh2RwvIiLyKK8Nz509exaPPvoooqKiEB4ejtq1a2PHjh1Z+xVFwdixY1GyZEmEhYWhRYsW2L9/v2vvgig3vR5o2RLo1096ee65R9v5993nuYAJkGKXq1cDxYtnt88y3BcUBHz3HQMmIqICQlPQdPXqVTRt2hRBQUH49ddfceDAAbz//vsoWrRo1jGTJk3ClClT8Mknn2Dbtm2Ii4tDYmIirjMRltQwmYAvv5RhsvBwIDpaepMOHbJ9fMuW6q9tNAJVqnikmVbq1AGOHwdmzQK6dQO6dJHA7PRpqWdFREQFgqbhuVdeeQV//vkn/rAzW0hRFJQsWRLDhg3Dyy+/DABIS0tDbGws3n33XQwcONDpPTg8dwfLzAR69AB++sk6V8nScxMbC9y8Kf+OHAn07y/r4Gkpovnbb7IcCxER0f/zyvDc0qVLUb9+fTz00EOIiYlBnTp1MGvWrKz9x48fx4ULF9C2bdusbSEhIWjevDk2bdpk85ppaWlISUmxetEd6qOPJGACrJO7FUV6oM6dkxyiw4elHlTx4sBddwHlyjmfwWYwyNCclp4pIiKiHDQFTceOHcOMGTNQqVIlrFy5Es888wyef/55zJs3DwBw4cIFAEBsbKzVebGxsVn7cps4cSIiIyOzXmXKlHHlfVB+pyjAhx9qO+faNaBMGWDiROcz6Dp0AH7+mYvb5ke7dkkh0W++AS5e9HdriOgOpiloMpvNqFu3LiZMmIA6depg4MCBGDBgAGbMmGF1nC7XB5OiKHm2WYwaNQrXrl3Lep0+fVrjW6AC4fJl61ICat2+LRXIFyyQxWxzCgkBHnxQygr8/DPA4d785fBhSfSvWxd48kng0UeB0qWBp5+WnzsRkY9pqtMUHx+PatWqWW2rWrUqFi1aBACI+//lLi5cuID4+PisY5KSkvL0PlmEhIQgJCREU6OpAAoKcv3c6dOBGzckCXvFCuDMGSAmBujYEQgL81wbyXfOnJG6WVevWm/PzARmzwYuXMjOfSMi8hFNPU1NmzbFoVyzmA4fPoyEhAQAQPny5REXF4fVq1dn7U9PT8eGDRvQpEkTDzSXCiRFkV4FV3uCbt2Sf4OCpAzBs88C3bszYAp0R44AQ4cCFSpID9JDD0kNLgB4/30JmHLXvgKkVtfPPwN28iSJiLxFU0/T8OHD0aRJE0yYMAEPP/wwtm7dipkzZ2LmzJkAZFhu2LBhmDBhAipVqoRKlSphwoQJCA8PR+/evb3yBiifW7xYqnP/+6/r1zCysH2+s3y5lGMwm6X3CACWLAF++AEYMwaYM8d2wGRhNALz5klvFBGRj2j6tGnQoAEWL16MUaNG4c0330T58uUxdepU9OnTJ+uYkSNHIjU1FYMGDcLVq1fRsGFDrFq1ChERER5vPOVzCxYAvXq5P8TCEgL5y8WL0hOYkWGdwG8JnsaNc36NzEwmhRORz3EZFfKP9HSgZElJAHfX+fPA/+fTUT4wYQLw+uvSy2SLwSB1uTIy7F/DaJR1CD/5xDttJKI7iteWUaE7nKJI3kmfPlK1u21bGUpJTdV2nV9/dT9g0ulkEV8GTPnLxo32AyZAhuXMZgme7MnMBJ54wvNtIyJygEETqWc2AwMGAC1ayJpqO3ZIhe0nngD+9z/g7Fn11zp9Wt2w3KRJUorgp59k6nmhQkCxYjIF/dIloFMnl98O+Yman3tQkCSH28pX0+lk7cG6dT3eNCIiRxg0kXpTp8p0byA7/8TSY3D8uCT2qh3tLVFC3bGJiVLAsksXCdJu3ACuXAE+/xyIitL8FigAtGrlOHAyGoHWrWV2XMeO1sdGRACvvSY/fyIiH2NOE6mTmQmULSv5Q45s2qRuLTjLGnI3b9rer9MBd98NHDjAWjwFzZUrQEKClIqwN0yXc43AU6eAPXukWGnTprKQMxGRBzGniTzryBHnAZPBAKxZo+56hQoBb7/t+JjJkxkwFUTFi0suWmioJHxbWIbipkyxnhFZtixw//3S68iAiYj8iAVuSB1HNXMsdDp1x1kMHSrnvP46cP26fK0oUs17+nT5oKSCqXlzCcRnzpRClWlpQJMmwODBkh9ny549wKxZcl7x4sAjj8jviKOEcSIiD+LwHKmTlibDadeuOT5u5UqZUafFrVtS7PDSJelVaNvWvWVVqGBRFCmA+v770huVmSmBkskE1Ksnv3PMbyMiN6iNPRg0kXqvvgq8847tPBSDAShfHjh0yHrIhchd06dLD5QtBoPM5lQ7LExEZANzmsjz3ngjO9ckZ2BkMABFi8oyGAyYyJNMJgnUHe3/7Tdg927ftYmI7lj8hCP1QkJkGG3OHClsWawYUK4cMHo0sG8fUL26v1tIBc3hw1LTyxGDQX4viYi8jIngpE1QkBQW7NfP3y2hO0F6uvNjdDp1xxERuYk9TUQUuCpVkvIUjmRmSs+nhcmUXXyViMiDGDQRUeAKDweeesp+WQGDQWZctm8PLFsmSeFBQfKqVw/45hv1VeqJiJxg0EREge3tt2WdudyFTo1GCaoWLQLee09qNm3cmB0k/f038OijwJAhDJyIyCMYNBFRYCtcGFi/XgKju+6SYKlYMeCZZyQwMhqBV16RY3MWV7WUxpg+nYniROQRrNNERPnbM8/IQtL28pgMBqBNG2DFCt+2i4jyDdZpIqI7w7ZtjhO/TSZgxw7ftYeICiyWHCDfun1b8k5u3gSqVZPZUUTuCAtzfkxoqPfbQUQFHnuayDcUBZg8GYiPl9XqH3wQqFwZaNlSChgSuerBBx1XojcagW7dfNYcIiq4GDSRb4waBYwcCSQnW2//4w+gcWPgxAl/tIoKgv79ZRkfW2UJ9HoJmoYM8XmziKjgYdBE3nfiBDBpku19JhOQkgK89ZZPm0QFSFQUsHo1ULy4fG8wSLCk08nQ3dKlHAYmIo9gThN537x58iGWczp4TpmZUoRw2jTmnpBr6taV4HzBAmDNGvlda9oUePxx6YUiIvIABk3kfWfPOg6aACAtDbh6VXKeiFwRHg488YS8iIi8gMNz5H2xsdmFBu0JCmKPABERBTQGTeR9jz7quJfJYAAeekjd1HEiIiI/YdBE3le5MjBwYN61wwAJmMLCgNdf9327iIiINGDQRL4xbZqUHMid6F21KvD770CVKv5pFxERkUpce458KzkZWLVKKoJXrw40aGC7B4qIiMhH1MYenD1HvlW0KPDww/5uBRERkWYcniMiIiJSgUETERERkQoMmoiIiIhUYNBEREREpAKDJiIiIiIVGDQRERERqcCgiYiIiEgFBk1EREREKjBoIiIiIlKBQRMRERGRCgyaiIiIiFRg0ERERESkAoMmIiIiIhUYNBERERGpwKCJiIiISAUGTUREREQqMGgiIiIiUoFBExEREZEKDJqIiIiIVGDQRERERKQCgyYiIiIiFRg0EREREanAoImIiIhIBQZNRERERCowaCIiIiJSQVPQNHbsWOh0OqtXXFxc1v5+/frl2d+oUSOPN5qIiIjI14xaT6hevTrWrFmT9b3BYLDa3759e8yZMyfr++DgYDeaR0RERBQYNAdNRqPRqncpt5CQEIf7iYiIiPIjzTlNR44cQcmSJVG+fHn07NkTx44ds9q/fv16xMTEoHLlyhgwYACSkpI81lgiIiIif9EpiqKoPfjXX3/FrVu3ULlyZVy8eBFvv/02/vnnH+zfvx9RUVFYuHAhChcujISEBBw/fhyvv/46MjMzsWPHDoSEhNi8ZlpaGtLS0rK+T0lJQZkyZXDt2jUUKVLE/XdIRERE5EBKSgoiIyOdxh6agqbcbt68iYoVK2LkyJF44YUX8uw/f/48EhISsGDBAnTr1s3mNcaOHYtx48bl2c6giYiIiHxBbdDkVsmBQoUKoWbNmjhy5IjN/fHx8UhISLC7HwBGjRqFa9euZb1Onz7tTpOIiIiIvEJzInhOaWlpOHjwIJo1a2Zz/+XLl3H69GnEx8fbvUZISIjdoTsiIiKiQKGpp2nEiBHYsGEDjh8/jr/++gs9evRASkoK+vbtixs3bmDEiBHYvHkzTpw4gfXr16Nz586Ijo5G165dvdV+IiIiIp/Q1NN05swZ9OrVC//99x9KlCiBRo0aYcuWLUhISEBqair27t2LefPmITk5GfHx8WjZsiUWLlyIiIgIb7WfiIiIyCfcSgT3BrXJWERERESe4JNEcCIiIqI7BYMmIiIiIhUYNBERERGpwKCJiIiISAUGTUREREQqMGgiIiIiUoFBExEREZEKDJqIiIiIVGDQRERERKQCgyYiIiIiFRg0EREREanAoImIiIhIBQZNRERERCowaCIiIiJSgUETERERkQoMmoiIiIhUYNBEREREpAKDJiIiIiIVGDQRERERqcCgiYiIiEgFBk1EREREKjBoIiIiIlKBQRMRERGRCgyaiIiIiFRg0ERERESkAoMmIiIiIhUYNBERERGpwKCJiIiISAUGTUREREQqMGgiIiIiUoFBExEREZEKDJqIiIiIVGDQRERERKQCgyYiIiIiFRg0EREREanAoImIiIhIBQZNRERERCowaCIiIiJSgUETERERkQoMmoiIiIhUYNBEREREpAKDJiIiIiIVGDQRERERqcCgiYiIiEgFBk1EREREKjBoIiIiIlKBQRMRERGRCgyaiIiIiFRg0ERERESkAoMmIiIiIhUYNBERERGpwKCJiIiISAUGTUREREQqMGgiIiIiUoFBExEREZEKmoKmsWPHQqfTWb3i4uKy9iuKgrFjx6JkyZIICwtDixYtsH//fo83moiIiMjXNPc0Va9eHefPn8967d27N2vfpEmTMGXKFHzyySfYtm0b4uLikJiYiOvXr3u00URERES+pjloMhqNiIuLy3qVKFECgPQyTZ06Fa+++iq6deuGGjVqYO7cubh16xbmz5/v8YYTERER+ZLmoOnIkSMoWbIkypcvj549e+LYsWMAgOPHj+PChQto27Zt1rEhISFo3rw5Nm3aZPd6aWlpSElJsXoRERERBRpNQVPDhg0xb948rFy5ErNmzcKFCxfQpEkTXL58GRcuXAAAxMbGWp0TGxubtc+WiRMnIjIyMutVpkwZF94GERERkXdpCpo6dOiA7t27o2bNmmjTpg2WLVsGAJg7d27WMTqdzuocRVHybMtp1KhRuHbtWtbr9OnTWppERERE5BNulRwoVKgQatasiSNHjmTNosvdq5SUlJSn9ymnkJAQFClSxOpFREREFGjcCprS0tJw8OBBxMfHo3z58oiLi8Pq1auz9qenp2PDhg1o0qSJ2w0lIiIi8iejloNHjBiBzp07o2zZskhKSsLbb7+NlJQU9O3bFzqdDsOGDcOECRNQqVIlVKpUCRMmTEB4eDh69+7trfYTERER+YSmoOnMmTPo1asX/vvvP5QoUQKNGjXCli1bkJCQAAAYOXIkUlNTMWjQIFy9ehUNGzbEqlWrEBER4ZXGExEREfmKTlEUxd+NyCklJQWRkZG4du0a85uIiIjI69TGHlx7joiIiEgFTcNzRERELlEUYNs2YOdOIDgYaNsWKF3atWulpwMGg7w8TVGAv/8G1q8Hrl8H6tcHWrcGQkI8fy/Kdxg0ERGRd/3+O/Dww8DFi9nb9Hqgd2/gs8+A8HDn1zCbgS++AD78ENi3D9DpgFatgJEjJQDT6vp1YPt2CZLq1AGKFQM2bgT69weOHrU+NjwcGD8eGDpU7nv1KjB1qhwfEgJ06QI8/ri690H5GnOaiIjIc37+GXjnHWDXLiAjQ7ZlZto+Vq8H2rQBVqyQYMQes1mCkm++keMsH1t6veyrWhVITZXAp3dv4Mkn5Wtbbt8GRo2SYC01VbaFhADt2wO//AKYTPbb8dZb0ss1frzcN6fQUOC774DOneV7kwlITgYKFZKeNT2zYQKZ2tiDQRMREbnvwgWgXTtgzx7t565dC7RsaX///PlAnz7qrqXXA7GxwIYNQKVK1vtMJuD++4FVq/IGPWrkDNjs3XvZMmDNGmDGDODWrex9d90FvPIK0K+fd4YVyS0MmoiIyDdSUoD//Q84ccK18yMigLFjgUGDpMcmMxM4dUqCkLJlgXvvBf76S32gYzBIwHTgQHYP1n//AS+/LEN83hQWJr1Z9j5aH3oI+PZbBk4BhrPniIjIN2bNcj1gAiS/6KWXZKjuzTeBkiWBihWB8uXltX27tp4hkwn45x/pwQKATz6Ra3o7YAJkyM9RX8T33wNz5ni/HeQVDJqIiCjb5cvAG29IkGEwADExMqyUM4k7t88/d/++ZjPw55/AmDHApUvZ20+dys6N0sJolCG6+fOB555z7RreoNMBH33k71aQizh7joiIxLlzQJMmwJkz2QnRly4B770HzJsnQU358nnPO3XKt+1US1EkAHSWi+RLiiKz/8xmJofnQ/yJERGRePpp4OzZvDPITCYJnvr2tX1eIH74Z2YC5coB//4bOAGThdHoeLYgBawA/E0nIiKfO3ECWL7cfnmAzEzgjz+A/fvz7ouN9WrTNDMapQxBtWr+bkleOh1wzz3+bgW5iEETERFJFWw1PTI7duTdlpjo8eZoYpmJptPJKyYG+OknSSb31Sy14sXVHacoMsxZuzZw8qRXm0Sex6CJiIikAKOrxz33nGfbopVlONFolDpMO3ZIyYGYGKBrV9nubVeuaDt+3z4ppXDzpnfaA8iyNY8/LsHj3XcDL7wgw5XkMgZNREQkH+BhYY6PMRplHbbcqlUDJk70Tru0yMiQ4pJ9+gBpabJtyhQgOto3gZMWZrMk3JctKxXULdXJPeX992UY8NtvgWPHgMOHZdZetWryjMglDJqIiAgoUgQYMsR+grJeL8uTlChhe/8rrwCLFgENG2ZvK1UKGDcOKFzY8+21x2yWxXYt0/rLlJE6T48/HpiL7l65Isu6REdLmz2RtL5+PTBihHydM0fNZJLAsnt3mSlJmjFoIiIiMX488Mgj8nXuXKDChaV2U84aSrl16wZs2SIVwvfskSDsyhWpkK2WZSaeO7PLzGbg44+zA5BSpYDZs6UGVWSk69f1plu3ZEHgqVPdv9YHH9jvWVMUCZw8UVvrDsRlVIiIKJuiAJs3AwMGZC9DkvNjQqcDKlQABg6UY4oWtT7fbJZep/fekwBIr9deWPKDD2Rh3MKFJch59FHX3ktKiizRAkguz7hxsuivK+vO+UqhQrKOnzu9c5GR8t4dadECWLfO9XsUMFxGhYiItNPpgIMHJWAC8g4XKYoEICNHynpzp09b7x8/Hpg8WY6zDAdpdeqUXH/QIKBXL1leRWstKJ0uO2l961agenXgq68CO2ACJDH8p5/cu4aaXjrWiXIJgyYiIsqmKNJLpOZD9dQpmZ1mceOGBE3uWrEi+2u9HliyRJK71ZYPMBiAtm0lhykzE+jYMTsxPD9ISnLv/JYtHSe+6/VyDGnGoImIiLJdvSqL3arN3NixA9i4Ub7+6CPPBCdXr1p/X6iQLONy6pT8O2uWTKO3FxhYhggBKdh5+bL7bbIwGKTSuDeTysPD3Tt/+HD7RUotPXBPPeXePe5QDJqIiCibK2muluTl775z//56ve317QBJRH/sMfnA37ABqFxZthuNEszo9UBQEPDll5KzA0ghSU+x5Gh99x3w11/azrXkVqnhbsXw++6TvDDAOrA0GCRg+vFHID7evXvcoQKscAUREfmVKzk/W7dK/tLu3Z65/8CBzo8rVUpm6P36qwzfpaYCtWoB/ftLUUsLT66LV7w4sGAB0KCBfF+3rlRSV/PMrl+XdqkZeqtQwa1mAgCGDZPg6ZNPJHC0FP589lnpKSOXcPYcERFl+/FHqeOjRXy8rF3n7pCVXg80bQqsWaO+Qrkz27ap67kJCnKetG40ylCf5bNpwwZJUrc3FJbbPfdIgOlI06bZw53kM5w9R0RE2qkNAHKqUkWCnJIlXb+vwSAz5X791XMBEyC9Qs56VoKDgaNHnVdEtySVWzRvLu2NjlbXlr//zluiIbePP1Z3LfILBk1ERJStfn3t09Etw2mPPeb6fU0mqaFUrpxUFvekjRvt1z2y5ChdvqxuKZM//7ROVG/TRobs1EhPB5KT7e/X690vN0BexaCJiIiyVagAdOigbnq/Xg+0a5c9nPfWW+7X//nvP6BHD6nT5CmlSgEXLwJPPCEz8YDstm/dCjzwAHDkiPrrLVli/X1IiGfWtjObWak7wDFoIiIia59/LgvJOkqijoiQ9c1++ik7YBg+3DNrpwFSIHPDBs9cC5Bp/LNnSy2pjAwZaluxQtama9s2e/kYNSwlDBRFgrtmzTxXNPPiRc9ch7yCs+eIiMhafLzUX5oxQ2oiXbwIxMYC/foB994rAVONGtb1hI4cAaZN82w7hg8Hdu707DWB7CAvNRVo1Qo4dEjb+RERkvi+YYMEd4DngqawMLmWJ2f9kcdw9hwREbknMxOoVk3bEJcaOp2soebOOmyOfPEF8OST2s7R67MDJKNRcrE8/TH6xBPS22dvqDM1VZaEWbdOhgaffhpo0sSzbbjDcPYcERH5xhdfeD5gAiQYuXHD89e1mD1b+zk5e5QyMz0fMAHyPO0V5Vy4UGbgDRwoCehz50qZgoSEvOsAkscxaCIiIve8+653rlu4MBAV5dlrWhYDTkgANm3y7LU9adasvNt++w3o2VNm4eV26hRw112SSJ/T1atSsb1rV+DBB6VSeO5lakg1Ds8REZHrFMU7+Tc6nVS1njLFc9fcsgVITJThLZPJc9f1hqiovAFQjRrA/v2Oz9PrpZeqUSNg/Xqgc2fg5k3rY8LDgaVLJZ+LAKiPPRg0ERGR665fz66Q7Sk6naw/t3Wr53qa0tJkptyVK4EfMFlcupRdOPPiRSAuTv25Dz8M/PyzvO/cSep6veRC/fOPzJIk5jQREZEPuJNz9PbbkvScczHbkBBgwABZENeTQ3M//CBBSH4JmABg/Hj5NyMDeOMNbed+9530qNma1Wc2yxDf9Onut/EOw6CJiIgkqfnNNyUvJj4eqFdPXXXq6Gjr0gNabNggydjXrgHHjgH79klg89ln6pcmUWvTJllfLj85cUKG1hISgJkzPXttkylvkU5yinWaiIjudCdPAjVrylCbxYULkjh8zz3A5s3285aCgmTK+9Sp2u+7di1w/LgMxZUv70rL1VNT4TzQXLoEtG8PnD/vnevbSignh9jTRER0p6tf3zpgymnrVqB3b8fnP/OMa/c1mYCPPnLtXK1at5Zhrvzkzz9l3TxvMBqBhg29c+0CjEETEdGd7Kef8s7Syu2HH2T4zp6lS12/v7fXWjObZar+wYOSI+WrStuNGvnmPq7KzASGDPF3K/IdBk1ERHeyzz5zfozJBCxbZnvf6NHuLa5744b3ijJ++aUUgmzTBhg1SmbO5U6Mtiyp0rixZ++9ZYtnr+eOnIGiZZhy7FgpikmaMGgiIrqTqc1rsTVL7ttvgYkT3W+DJ3ubDhyQatlhYUD//tbDjjkr7BQpAtSuDTz6qMzU27QpO4AqSEaNAlq0kNyzoCDgvvukFMGYMf5uWb5UAH9DiIhItWbNZPjKmZYtrb/PzJQFfD3BU4vyLl0KdO8uvUnOFtBNSZFK5m3bZm+77z5JTi8oPv44ewjOEjDaW8+OVGFPExHRnWzUKOczy6pXB0qWtN62bJnnZl8tWwY8+6z9ZHQ1kpKAHj0kmHMWMAHynr/5xnrbSy+5fv9AlHOZGJ2OAZMHMGgiIrqTBQc7Xrg2NBRYvjzvdk+u26YowKefyjDSrVuuXWPECG2z40wmyXHKqX377IKS9hgM+aeK9sKFUuCSPIZBExHRna5vX2DNGqBatextBoMMXR0/bjtI8MYstF27JHjS6vZtya/SQq+3XRtq9GjJcWrXDihUKHt7uXJSHuH2beDJJ7W30V16vQSwWgp0ms32E/jJJQyaiIhI6hjt3y+9PhkZMsy1cqX99c68sdirogAzZmg/b/FixyURbDGbgaeesr3vnnuAFSsk+T0zU9ZvO3RIgqwvv5ShQF+yBExGo/ZaUykp3mnTHYqJ4EREZE3NLLLq1b1z75MntZ+za5f2c55/HqhVy/lxBgMwb54M/+UezvOmu+4CSpSQtfgSE4EKFYBevbRfR817JNUYNBERkX1JSTLEc/Ys8L//AfffLwnFR4965345F+91RlEkt0rLkJ5OJwHQu++qO37OHFlU2NeOHpVXuXIyTOjKOnGlS8saguQxDJqIiCivtDRg0CAJGnLWNwoNBd5+W5K2PU2vl7pJaly9KuUF1q1Td7xOJ7lbY8fKArhqpKf7f0bdiRPAsGFA8eLaz42JAb77Tp5TQaxB5QfMaSIiImuKAjzyCPDFF9YBEyCJ0CNGAD/+KD0ZnlSoEDB8uLr2PfAA8Pvv6q+9aJEEgGoDJkCS4y9fVn+8N7kyNLhrF9Czp/QOpqV5vk13IAZNRERkbds2WZPOkYkTgeee035te4FW6dLAhg0yHOXMH3/Iy2RSd89ChYCuXVU3McvFi9rPCSSWgHflSuD11/3blgKC/XVERGRt3jznxyiK9Do1ayYBjFoNGsjMtHnzZLZeqVIS0LRtq66Mwf792hea1euBd94B/v0XOHxYhq169gS6dHE8hd/TPWn+9MEHQJkyQMWK8qw5XOcSPjUiIrJ27py64/75RwInLS5dkpl3ahOxAWDHDllYeP164MgRbfcDpNL4qFHW2374Qdqxfj0QHS3bzGbplfn6a+llKlsWiIoKnCE6d2RmyoxBQAqavvwyMG4cq4RrxKCJiCi/SksDzp+X4acSJbSfn5kpeS9ms8yysvQ+qK14XaiQ9jpARYqoP/bCBekR2rBBPtxz51e5a/9+oHNn4PvvgWPHJJDYskXKDJhM8jzs1X/K2R69Xt3SLYEiPR146y3p8du82fkyOpSFOU1ERPnN1asyoyo6WgouxsQATZtKL4kamZnAQw8BYWFSyLFRI5kV9+CD8oGqtuJ1oUJSP0gLtUHW559n5zkBng+YLLZskWGr5s3layA7VypnwFSsmPV5VarI8jJffQXce69cIz7eO230lm3bgJEj/d2KfMWtoGnixInQ6XQYNmxY1rZ+/fpBp9NZvRo1auRuO4mICJCAqUkT4JNPpGK1xZYtQIcOwNy5js83m2VY6ocfrIMCk0mSv+++W/ar6RHasgV4+GFt7Vcz1LViBTBggPpEb18wmWQob8kSYPt26aXq0EFKJGzYAJw65drQoYXR6Fpvobs++shzCy/fAVwOmrZt24aZM2eilo1qo+3bt8f58+ezXsttLfZIRETavf22fDjnDijMZumNGThQAit73n1XkqHtOXECGDNGXZHJq1elNEFsrKqmA7AO9Ox5+23vrG3nDksP2QMPyFBm7lwgRZHhxPr1Xbt+ZibQv797bXT1vpYeNnLKpd/KGzduoE+fPpg1axaK5e6yBBASEoK4uLisV3FXinIREZG19HQZtnLUA5OeDnzzjf39H3/s/D6ffio9To5yXQwGOSYsTGbEqVW0qOP9ycnAn38GZo6QvR6ZxYtluZK77pJeKC0sweHQocD777vXPldduJB32+HDkjzfpw/QowfQuLEMQd51lxQ93bfP9+0MAC4FTYMHD0anTp3Qpk0bm/vXr1+PmJgYVK5cGQMGDECSrxc3JCIqiC5dcp4TZDQCBw/a3//ff87vc/Uq8MwzjoMzkwlo1w74+29g7Vrn17S45x7H+1NT1V/Ll4xGWUYmt1mzgG7dZLjOFaVLA99+KzlT/hqOrFs3+2tFkSrod98NTJ4MzJ8vhUG3bAHOnJGyDTNmADVrAtOm+ae9fqR59tyCBQuwc+dObNu2zeb+Dh064KGHHkJCQgKOHz+O119/Ha1atcKOHTsQEhKS5/i0tDSk5ahUmsIVmYmIbCtUyPkxiuJ4aC0oCMjIcHwNo1ECgfvvl3Xn7CVha62XBDgviFmiRGBO83/4YUm4zyk5Ofv9uJqofuqUvNfPPnOreS4rVUp6jyymTAHee0++dhbEDRkiwdN993mvfQFGU0/T6dOnMXToUHz99dcIDQ21ecwjjzyCTp06oUaNGujcuTN+/fVXHD58GMuWLbN5/MSJExEZGZn1KlOmjPZ3QUR0JyhaFGjZ0vGwmWVmnD0tWzq/T7Nmco9Fi6SSdM40DHempzdoYLu3JiejEXj22cDLabp8Weo95TR/vmeSqF96SUpH+MMXX2R/nZ4uRUC1yF3/qqBTNFi8eLECQDEYDFkvAIpOp1MMBoOSmZlp87y77rpLeeedd2zuu337tnLt2rWs1+nTpxUAyrVr17Q0jYjozrBunaLodPKS/o3sl8GgKG3bOj7/6FFF0evznmt56XSKsm+f9Tm3byvK7t2K8tZb9s9z9ipSRFGOH1f3Hq9fV5QGDey38/77FaVOHdfb4uorKEhR3ntPUcxmaefw4bLN1+3w1GvUKOvn/uefrl1n+3Z1P1dHLl1SlI8/VpRXXlGUDz5QlPPn3b+mBteuXVMVe2ganmvdujX27t1rta1///6oUqUKXn75ZRhs/AVy+fJlnD59GvF26leEhITYHLYjIiIbWrSQHJgnnpD8H6NRProyM4HERGDhQsfnV6wI/PyzzALLXbjRYJBCj9WrW28PCZFE5/btXW/3nj22F8s9dkyGpn7/Hbh1S4aEjhyRXo8SJaSAZ0qK9Dx16iQ9G2XKqC/A6UkZGbJY8ZUrwPjxUjnd2VCnWt4o3unMqVPW37uaT9aokQxVqhk+zk1RZEbn66/Lz95SKHTECOCVV6QIZyBVLXc3OmvevLkydOhQRVEU5fr168qLL76obNq0STl+/Liybt06pXHjxkqpUqWUlJQUj0Z7RER3tJQURfnsM0UZOlRRXn1VUXbu1HZ+WpqijBmjKPXqyWv0aNnmiKMeKmevzZvzXm/+fEUxGqWHzNY5lt60tm2l98li7lxt965SxbM9NDqdokyf7v+eIndfBoOinDuX/VzPnXP9Z9yvn7bfP4uPP3Z83fHjXbuuRmpjD52iuBfatmjRArVr18bUqVORmpqKBx98ELt27UJycjLi4+PRsmVLvPXWW6pzlVJSUhAZGYlr166hiJZy+0RE5F1BQfaXFXGmUiVg06bsdd727AHq1FFXWkCnk+n43bsDX34JrFolZQnUqFcPqFZNyjB4soyBoyVW8pO5c4HHH8/+vnt3KXKqdSafTgfcvCklKNRKT5c8PUc9XGFhMmvUlV4sDdTGHm4HTZ7GoImIKEBVrw4cOOD6+UajLBr77rtShHPePPWBR9GiwLVr2Qniaj7UK1WS9j73nNS3KghBjqdVqCABrCUoOXdOhttOn9Z+rT59ZLFjtRYtkhpQzpQoIcOxjz8uBUDVFF7VSG3sEWDTE4iIKGDZmQWtWmYm8MEHEjitWKEtiElOlgEbk0ldwFS9uhRoNBqBrl0ZMNlz/Hh2iQEAKFlSCnQ++KD2ay1YoK03T+1aiZcuATt2yHqLtWsDZ89qb5uHMGgiIiJ1ypUDBg927xqKIhXHvb3e2bx52V+3aSPDdEbNpQkLPkUBpk+3DnZiYqTKec2a2q5lMskQrFpay1coigR5Wtc79CAGTUREpN4nn0i1cHfo9bJenTeDmJy9GHo9sHy55FABcl9v1oHS6bJzt7TwV22qpCTbleY/+ED7tdT2HgHWlcjVUhQJzP7+W/u5HsCgiYiItJkxQ0obFC7s2vk6nQyfeWvITK+XQCCnmBjgr79kyZdBg4BHH/Ve0KYowMiR6o/X6aRchD/X27NVsLp1a5n2r4WWEkIPPwwEB2u7vsWKFa6d5yYGTUREpN3DD0tuiSsVwk0moGlTYMIE+d7TwYvZLLWcctPppCL6hx/KrLFhw7xXA0htsBEWJrWnfv0VeOMN2eaPYcR//7W9feJEYMAA9depVk39sZGRsr6dK44ede08NzFoIiIi1xQpIvkwWhkMQO/eEiysXy9r3EVFySyphg3dD2T0epnJ5czEicCTT2a3yZMBlJpeowcfBC5ckEKZBgMwbpwsjNu7t/TE1a8vxR3HjpV8Mm/R64HvvrO/v3VrddeJiJCfpRbPPw/Mng3YKYBtV+XK2o73EGbFERGRNooiwc6nnwLr1mk//513snN+mjeXl8WsWcDWre61b+xYyZlyxmiU+7VtK0GWL4fHjEYpo5B7envDhvLKbcwY4MUXZUFdT9PppJyDPfffr64u1QMPuDbc9sQTUk7g/vuB1avV/RweeED7fTyAPU1ERKSeosjq9q1aSZ2dS5e0nT94sHz423P1qnsJ0WPGAK+9pv54RQFGj5YPal+WLVQU20OI9pjNsnyON5hMjgtMFiqkrrjkkSOut8FolOFJNQFTkybA3Xe7fi83MGgiIiL1Zs3KHpLTWjV6yhSZfedIhQrarwtIb8lLL0kvk5ZhtvXrJT/GlXu6w2wG+vZVd+zVq8A99wDnz3uvPbNnA7nWls2Sng7cuOH8Grdvu9eGJk1kuM6R2Fip7u4nDJqIiEgdRZHEXVdyf7p3B4YPd35c585A8eLa7mE0yofpiBHa27Vnj3+m+o8YIYsnO2M2Ax07Art2ebc9Oh0wdartfZMnOw8q9XqgQQP32zF1KjBzZt4crrAw6aH8+2/v5nc5waCJiIjUuXhRemVcGcaqX1/dcSEh0psF5A1mDAbrmWWWwKpOHWDjRikroFVYmO+n+o8cKUvJqLF6tSSHe7uNmZmy5lxu58+rG+40m6WUg7t0Opmtd+wYcOYMcPIkkJYG3Lollcvj4ty/hxuYCE5EROq4k/PTsaP6Y7t1kyn4o0cDO3fKNp0OaN9eFu4NCgLWrJEP+oYNpdq3Wikpsujv11/LsFe5cnJtX+Uz6XTSfrU9ad9957vFgdPS8m5LTFR3bnR0dvFQT9DpgFKlPHc9D2HQRERE6sTGAgkJ8te/Vvv3A7VqqT++XTt5HT0KXLkiC7bm7GV4+mntbTh+XGbqnTkj31uW5fBlArheL4viOpOSIkvB/PSTbwImgyFvhe5t2+Tnpoafe4B8hcNzRESkjl4PvPCCa+cuXeraeXfdJUnQ7n4oK4rURTp/Xr62BEq5c3VcKdaphU7nvCbRv/9KnabnnwcuX/ZueyxMJuC556y3rVihvkdMS0CcjzFoIiIi9YYMAZo1037eggWSdO0vf/wh93fUa1OsmOTTeLMid2am1CSyx2yWekWW4M7bLEHRU09Jsn5OGRnqg6Zu3TzbrgDFoImIiNTT6+WDXeuMM4MBGDgQmDRJFlz15ZAYIKUFnAVDV69Kb8uzz3qnDTqdPANHNYbWrAH++cd3JRB0OpmtNnNm3gCpQQN1CejBwUCnTt5pX4Bh0ERERNrs3q19NpfJJLPARo+Wdefq1pV8Il9RG6QpirbimDkNHQrcvAn075+3GGRYmLz3adMcX2PdOt+uPWc221//rWNHKcDprLfp7bdtL/hbADFoIiIibVxdmR7I7kHZt0+G+a5e9UybnLnvPucJ1VFRQKVKsgZeVJT6axcvLrP6PvgACA8HvvhCikHu3g3Mnw8sWSLlGt5+23nOlNnsvUWE7TlyxPZyOAYDsHixrClnr02jRrlWHyufYtBERESOpafLB//UqTJVv1Mn93tDMjNlFpmlJpMWKSnAxx9L7aeEBKBlS2D5cse9SS1aAFWr2h9W1Okk8To4WL4eNsx5O2Jjgd9+k/yjF17IG1jUqgX06iXrpEVE2L7Ghg3AQw9Jj07FisDBg5JL5GujR9veXq+eBH/PPy/BodEoAWW3bsCJE8CECfK+TSZJ9h8xQl5Ll/q+yroP6BTF1wPLjqWkpCAyMhLXrl1DkdwLGRIRkW99/73k+Fy+nJ3PFBwsQY8nCi4GB8vw0KBB6gKx48eBe++1PW3/nnuAP/+0fR1FkQBm4ULb161XD9i8WWooATLMVr++5BfZYjQCv/8ONG7svM32vPEG8NZb1nWYDAYJNizP2leCg2UZFFd6uQ4ckED6xIns55eRITWwli0DqlXzZEu9Qm3swZ4mIiKy7ZdfgEceyZ72bvkQT0/33DBSerr06jzyiPOeCUWRPBt7dY62bgU6dLC9b+FC+wETAOzYIcGbZY21QoWkTlG/fnmDsNq1gb/+ci9gWrZMAibAetjQZJLn6usq5enpwKOPak/Q/+8/6cU7fVq+z8jI7ik7eVKGYH1VNsEH2NNERER5KYrUCvrnH9/NdJszR4IUe9avl6E4Z44ckfpOOTVpIoGOs2CkQgW5T5ky2dvS0qRcwfXrQI0ari3Xklvr1jI05yhQ9HVvEyCV2Nu3V3/8O+8Ar77quJ01a0oPoL0hygDAniYiInLdvn2SX+PL5UU+/tjxMb/9pu5athae3b5dXQBy6pTk6+R83yEhMv2+VSvPBEyArJXnrGfN1wETALz0krbjFyxw3s69eyVo9UVlcy9j0ERERHn9959v76cokhvjyO3b6q519mzebZZcG2cyMyXA+usvdcc7k5ICHD4MXLpkvd3XM+TU2r9fW7CWkqLuuH37JBE/93PIZxg0ERFRXmXL+v6eISGO97dooe465crl3daxo/oZfwYDsHatumPtOXECeOwxWcj27rulhyoxUWpVATLM6O0lW1wpDaEo2gLGGjXUv4+jR6XHKTVVe7sCBIMmIiLKq2JFmaXm7Q/2nHIv45Fbhw7q8mJsLVPy4ovapsBrHZa8fl2Su7/8Evj5ZxnOW7DAunzAunWSGL1qlZQo8PaU/Lp1XevR0pK4/cwz2t7H0aOy9p6nevJ8jEETERHZ9uGH0luRO3Dy1tBSYqLj/Xq9FIt0pGtXoE6dvNsbNQLmzlUXBJpMEjCqkZEheUBRUbJmXP/+QJcuMryZO4fHZJJXv37S02SpxO2NwFSvl0T+3JXJ1ShfXv2xHTo4Tt63JSVFfta2hlEDHIMmIiKyLTYW6NEjb5DkbBjNFXq99MA4c//9MsOrWDHr7TqdfHg7Cqoee0zqPDVtav8Yg0GGnO67z/4xyckyxf72bbnme++pL0ipKFIMc8UKKQK5Y4cU6PQ0s1mKa27apG3WWsWKEmyppdMBs2fLz0UtRQFu3QJmzFB/ToBgyQEiIsrr0CHpbbl61evDSP8hCrMwAN/q+uCaoRiql7iEZ4aG4P6XqtpfF9hsluGuvXtlXbc6daTY5JUr0lPyyCOAvc+Q9HSp0r1iRXYxSUACt9hYuU7ukgWAzHgbN04W1QVkvTW1yek5GQzAu+/KkOHNm0DRop6fWabXy9It0dHSxs6ds9vt6JzNm6VIqFYpKUDDhvaLgdpSubL8ngUAtbEHgyYiIrKmKJIPs3evVwKm84jDN+iDcygJPcyYg35IRjGYoQeggwGZMMGIntX34OvdtRyPXmVmykK5M2ZIr4fBINtCQ4GPPgKeesr+ed9/D3z2mdR1KlpUeo2eflpyer77TgLGihWlkvjatbLcCeCZUgCffSb3OnPGuiaUpzz4oKwbl9PGjRKobd2a9/iyZYFFi6QKuquuXAGaN5eZcmqULp1dFNPPGDQREZFrtm6VXgMPUwCMwxi8jdegQAcDTMiApRRA3jwpHcyY1PcARnxZw/5Fhw6V+k72Psp++MF5grnFtm2Sk7R/vwRgluVNvPEx2bs38NVX0gtUtKjn15v75BNg8GDb+65cAZKSZIafySTDcbZmHLri1i1Jdt+50/FxBoMU0fzlF8/c100sbklERK7Zvt0ryd5TMQzjMBYmGGGGARkIhgRLtu+lQIep38bY7+y6cAGYNs1+UKPTSW9O164yo+6nn2wPg6WmSkHLe+6RgAmQa2ZkeK+45/z5wKefAuHhwMMPu78Acm5RUfb3FS8OVKkiQUunTp4LmAB5P3/8AQwc6Pg4k8l+UBfAGDQREZG14GD3g4WBA2W4p2RJQK9HGoLxJt7QeBEdzqbH4OQxO1HT4sWOh8oURXpVliwBvv5ahqzq1pVgK6enn847lOVtOh3w/vvSxrFjZZabp2bR6XTqZ/95Q3i4BIQjRsj3ORPTLMH4009rW64lQDBoIiIia4mJ7vU03X03Tgz/EG/u64aBjXbjDeMEfIteSEYx5+faYLcp166pDzQsQeC+fTLTy/L9iRMSUPmaogDHjkkAd9ddsjZbgwbuX1evl9yr0qXdv5a7Jk2SIcgaOYZXq1QBPv9cgqpArYrugIf7A4mIKN9LSJDZZ99/ry0R3GiE8nhfvFLsM0yuaoBeD+j10TCbR8JkZwjOMTPKBF9CQoXY7E0XLsiaeKGhMktO66wzRZFp/r//LknLS5e60C4vqF5dZq79+ae2XiKdTt6T5d+6dSXJPBDodMCjj8orJUXaV6RIvgyWLBg0ERFRXrNmSYCyfr31tPzc9HrJkZk2DejWDZPeN2LSK7LLUsvRXs6SczoM73UBen2stOX554Eff8xuS3S0DAXduqX90kuWSNB0/bqLbcula1epdH3oEBAXJ3lC48fbP16nk6AvLs56e+PG2QnoajRpIkndpUpJEvsjj3injpbFjRsSAEVHWy/Tcu6c/A7Mny/7q1YFnn1W2mM02i//kM8waCIiorwKFwZ++01eX30lH8xFi0pdoVWrpNZRZCTw1FO4OvAVLN0UjaQPgLfecnZhBWqDqIpBpzB0di2prt24sUzPzxm8WRYV1unkpaUUwJYtUuhSS+FHR2rVktwki8xMYN48CSZsBZyKItP/c/e6nD0r+U3Xrqm770svSc0pb9u2DXjzTWD5cnnO4eHAE08Ar78u7/Hee+V3w+LPP+W1YIEEumoXTA5wLDlARESapN3MxLzZGfhhWSi2btUhOVnb+XqYYIbzXKRQpOKWKRS6Ua9I0rSjocK6dZ1Pc7elXDnJa3JX9ep56xPt3SvLpVy9mh3QWXqRnnhCevNyJkmvWQO0bastCX/xYklw96bVq2WWndls/TPQ6aQye3Ky44D1jTekKGgAY50mIiLyuLFjZdTJnQLW0frL+M9cHGp6nNJTbiMooaQEHvYYjbIA7tChQLVq6ntpPCkyEjajx0uXJPF5wQIZCqxRAxg0CGjXzrqXKT1deve01ms6fz7vEJ8nZWRIUvmlS67PqAwJkXX/Zs+WBXujo6UERN++nuvpc5Pa2IPDc0REpEqPHlJFwF1XzZFQEzDF6S8iyBDhOGACpJfj5Ekpb+CvfoASJexvHzVKXo68/772gMlg8G7ABAA//yxDs+5ISwN69szOjTtxQmqBffABsGFDYMz0U4klB4iIyKnvv/dEwKQAUGBS8fe6AZl4ptkBmSUXGur4YL1eei8AGSazu2Cdl+j1wJNPuneNn37Sfo4lAPGmAwc89zwtQ3uKIq9Tp6SwZz7CoImIiBzKzLS/hJtaep0l58V5D5MBmbg7+DiGf1VXPrD79HFcMTszU6a1A1Jl2hNrw6llNMracc4qYDvjag9Z8+bAkCHA5MlSPiE11b125FaokPeeZ2amlFnYscM71/cCBk1EROTQu+/KLHJ3RBhToYfzmk8huI0n796IjQejUaRMpGwcNUpma9kqZKnXA507Z6+V17OnrDXnq1pAVavKsiHFXCvcmaVtW9fOO3VKpvqPHCmz6MLDJdApXVrWgBs7VmbkucrbSeZ6vQzR5RMMmoiIyK7MTEk9cVdE0G3o4bw35ezB6/jsnxYoViFHEFKxonywVqwo31sCIr1eepgWLszeZjBI0nXnzr4JnPbtk9IF7nr5Zc+199YtCZQ2bpRZa6VLS86UKzxVx6qAYCI4ERHZdfIkcPmye9fQ6YA0QzgyYb9WjwGZaBa5F1FV6tg+oHZt4J9/pFdn927Jc+rQwXYS8Zgxvqv0rSiSIX/xontBT+HCQL9+wJw5HmualREjZJjtpZe0neeJzH9HzGbgvvu8ew8PYk8TERHZ5YnOD0UB/rvhOJnbBCNGPuckH0enkw/Y554DBgywHTCdPg1MnOhGa11w6ZIUcnSXt4epRo6UGYa//KL+nOvXvddjZzTKsGr9+t65vhcwaCIiIrsSEpxPXlNDUXJ/8MpQnQEZ0MGMDzAMHRb0larf7vjmG+2zvSIjZdFcwDpA0HKd7du13dMWX9SXOn9ehi7V9mjdfbf3yjiULCnTMvMRBk1ERGSXwQD06uWNK+sAKGiCTTiIqhiGD2X6fJcu7n1Inz2rLthZvhz44Qepdp2UBBw5Avz9t0wTvOsuoFIloGxZ9fd1Nffn6lVg61bgiy+slyHxtoEDgStXnB/Xq5f31rJbuFBmHuYjDJqIiMihTz+VlBtviMAN3I3D8k1mJrBrlyQwuyo21vkUeYNB1krr3h1o0yZ74dn//Q+YOVMCqPXrZWaaWq1ba2vnxYtAq1ZAVJQMUT35JHD7trZruCMjQ9YUdKZIEeeFOV1lSezPRxg0ERGRQ8HBwL//eqf4tJK7bpPRCKxc6foF+/RxHDQZjUD79sC33wKffCK5SLZ6tubOVX/PiAigUSP1x1+6JD1Z69b5r4K5Xi+J9Wo88ojn728wSGmEfIZBExERORUTI4vZf/st0KABEB8PVKkCDBsmo1yvv649X1gPM5ojV/KzTufewnblywPPP2/nhnoJUpYvl+Gp55+XHqeaNfMutrt+vfqijlOmaMt/6tUrMKbyq1337e67ZU0/TzKZfDfD0YO4YC8REbnt+nWJV65eVRtrKDDAhPOIRwn8Z71ryRIp1Ogqs1nqE733ntQssihcWL7P3UCDQYahdu2SzPfTp+XNmJwX48SrrwJvv62+bZmZ0nUXCB+9f/0F3HOPumOXLAG6dvXcvXU6+fm88ILnrukGtbEHe5qIiMhtERHAqlVaZtrpoACYgFH4Hc1kLp3BAJQqBXTq5F5j9HoJmi5ckDpDX34JTJ8O3LhhO6IzmSTqsxSA/OgjdQFTYqK2gAkAjh71f8Ck00kuV4MG6s958EFg9mzPlR9QFOmuzGcYNBERkUfUrautsLUZRnyM59Ecv6Mh/sLFsHKycK2jdea0iIgAunUD+vaV9c2crV9nyWNasEDd9V2pBxUIeTyJiTKmqvoHZZZA8oUXPBvwlSvnuWv5CIMmIiLymHbttH2umv5/YYqdunpoV/YATLXreadh//3nPFcqJUUChP/+c3wcAAQFAfVcaGvp0t6bimiPJThKTJTgceVKqU2lhqIAzz4LDB3q+TpSI0d69no+4FbQNHHiROh0OgwbNixrm6IoGDt2LEqWLImwsDC0aNEC+/fvd7edRESUD9xzj7ZRHwuTYsDuA8FYscLzbQIgNZec9WDFxWUnizvjao+RTgcMH+7auVoZDFJjqUMH4LffZPy0bt28x/37L7Bli+3Con/8IWUYvGHTJu9c14tcDpq2bduGmTNnolatWlbbJ02ahClTpuCTTz7Btm3bEBcXh8TERFwPhJkCRETkVTqd5Ea7eq7XCkQ/8YTjniaDQWbUAeqKWlav7npbxo2TqtyeVqQI0KSJrDN39Ki839u3gWXLpCZUbmvXSpR7111A48ZSaDIxUdb2AyQh/oknPN9OC7WzEwOIS0HTjRs30KdPH8yaNQvFimWvRK0oCqZOnYpXX30V3bp1Q40aNTB37lzcunUL8+fP91ijiYgocLlao1FRsj+vPa52bRlmssVgkNlyQ4fK94MHO7/eiBG2tysKsHkzMH8+sGIFkJ6e9xidTmajeTqn5/p1eS8TJzovHLlsWfZwXU7r1kngtXy5rAn377+ebWNuN2549/oe5lLQNHjwYHTq1Alt2rSx2n78+HFcuHABbdu2zdoWEhKC5s2bY5Odbri0tDSkpKRYvYiIKP+KinL9XEsHiVd88okEFMWLZ28zGoGHH5ahIksnwJNPSl0iW7WXdDqgRQvbPUUbNkjxqiZNpMhmhw4yQ+zTT/Meu3mzLBvjSYoiw2mzZjk+LjNTFjxWlLy9PSYTkJYGPPoocPmyZ9tny/nz3r+HB2kOmhYsWICdO3dioo1ZAxcuXAAAxMbGWm2PjY3N2pfbxIkTERkZmfUqk8/WoSEiImvNmwPh4a6de+OGF3ub9HrglVekSufGjZLnc/as9AqVKJF9XOHCUtzy/vutzzcagf79gV9+kR6dnDZtkmn8R49ab79yRXq4pk7Ne3zua3jKoEFAs2YS/NiyerUEK/Zyt0wmKbilpuyCu9QmpAcITUHT6dOnMXToUHz99dcIdVCMQ5drGqOiKHm2WYwaNQrXrl3Lep0+fVpLk4iIyMf++w+YPFlm8z/0EPDZZ9ajLEFB2gpk55aa6n4bHQoJAZo2lTyfmBjrfadPA889B1SoIBWrixSRXqUvvpBga/ZsSQI/c0Zmf8XHA2FhEjCZTPbzdEaPtq4CbjB4t17Txo0yu89We44d81y9JXcEB+d9/gFOUzGMHTt2ICkpCfVyTLM0mUz4/fff8cknn+DQoUMApMcpPkfRqqSkpDy9TxYhISEI8dYKykRE5FHLlwM9ekgnhqLIZ++iRZL8vWKFpMEA7qw9q6Bvq9NoUOoc7mtqxqPv1UaROBe7rbQ6dEjWkLOUHgDk6+XLJffHsijv/v3Sk3Ptmvpk5tRU4McfpWYUIPlE3k6E3r9fAr8HH7TeXqyY/wtsAvIM8xlNfwu0bt0ae/fuxd9//531ql+/Pvr06YO///4bFSpUQFxcHFavXp11Tnp6OjZs2IAmTZp4vPFEROQ7//wjn7+3b8vnvSUlRlGA5GSgbdvsNJicaUPa6HAsowwWnmiIwd80Qql4E9ZM3uWZN+BIWpr0PiUn287zSUqS2XWKIj1P6teLEUajdf5OzZpAy5aeK+Rpz7vv5t12//3SO+Zvrv+S+I2moCkiIgI1atSwehUqVAhRUVGoUaNGVs2mCRMmYPHixdi3bx/69euH8PBw9O7d21vvgYiIfOCjjyRmsNVJYTJJx8ucOfL9gAHu3En3/y89bqIQ7h9ZFYd/9cAsrkOH5E28954kbed8Iw895DjxOTNTikJ+/TVw/Lj2e2dmAiVLWm9bsEASx73p1Clg1CgZPrz/fklKd6cuhCflsyRwQOPwnBojR45EamoqBg0ahKtXr6Jhw4ZYtWoVItSupkxERAFp8WLHM9vMZhkNGjFC0oJmzZIOGnco0CMNIXjvqX8w86yTafT2JCfLbLBlyyTZSqeTKK9KFSkMFR4O/PyzisYoktPkivDwvAveligBNGwI7Nvn2jXVOHdOEtBMJnnfy5cDY8ZIMnhmpqyd57Xpig7o9dbJ9/mETlECYWAzm9qVhomIyLeiomQymCMNGgBbt8rXx44BPXsC27a5f+9Q3EJqepBkmWthMgH33iuNyD0bzGCQRO9Bg4AJE9Tl+VSqBBw5oq0NAPDhh8Dzz1tv27JFikr6msEAREfLDygpCahcGcjI8H07Fi2S2QQBQG3swbXniIhIlQYNHM+SNxqlwLRFhQoSQP39t/Q65Sjhp9lthOHSjB+0n7h8uQQntqbPm0yS6P3bb+qm/4eEyOw6LUJDZXph7oAJAD7/3Ps5TbaYTMDFizI8WK6cFPX09Wy6Bg2ALl18e08PYNBERESqPPec49I9mZnAM8/k3f6//wFPPQUcOODO3XXYOXev9tPmz3ccEJlMMstMzRBVerr2aYFmsxTPtOXwYc8MjblS30GvR9ZCf+PH551h501160p+mD8CRjcxaCIiIlU6dpTACbD+nLbEJB98ANSoYf/8ixfdu3/qVRfqGPz3n/MijampMkznrLfFlWyW9HQplGlLVJR7Ba0sPv8cWLNGEr7ffRc4eND5dRUle0guOFiKb/oi97hqVSnfULSob4pnehiDJiIiUkWnk9SchQtlGE6vl86C1q2BVauAYcMcnx8d7d79g6Nd+FCvWNF5j0ZCggQegGeCmNzsVebu1cv9Wk1GI9C7t/wQJkyQgptVqgB16jh+LzqdDJFZvP8+cOuWe21xRq+XWQKdOklumtEI1KolyfX5ZPFeBk1ERKSaTiejTZs3y8hSerqMtCQmOj+3f3/37h3RsLr2k5580vEQmF4vY4oPPSSRX6NGrjfQnrp1bW/v2tVx15wa0dGSa5Xb0KH2AxGdTgKWJ5/M3jZ3rvd7fmrUkHHalSuz77V/v2x77LF8ETgxaCIiIpfodNryh4cO1T75LZsZdZ9rqv20Bg2kKKUtBoMUmXz2Wfm+TRvgzz9lPToHS4VpEhEhM+5sCQqSHiF3VKhge/ujjwJPPCFf58zpMholUPzmG8CyUoelyJa37dkjw4I5gzNLoDR/vtTACnAMmoiIyCdiYoDu3V09W4+Xux2BOdOF3ojp0yXXJ2ddoNBQCSo2bJC15HIqWVK60DzB0SL06eky7d4dgwbZ3q7TyZDjggVSCyo0VPK2evaUKY09emQfazD4fw04vV4KjwY41mkiIiKf2bXL/miVcwqeq74WH+1r7drpGRnA7t0SrFSvDkRG2j+2RAlJIneHwQAMHiyJYLacPg2ULev69WNjpeJ3cLDr17B47TVg4kT/DpEFBXkuWNWIdZqIiCjg1KkjnR2u0eHj/a1xeOlB104PCpIVhZs0cRwwAe40MpuiZA/92RLu5kLEmzerD5hmzpSaTEFBck6DBtLLZjF8uCTEe6Nek16vbrjT9bFbn2HQREREPvXtt+50jigY8vh1KGYvD5K8957zoCY4GChfXr7OOUPPaJTgY9Ysx2vLBQcDhQu71r7o6Ox7O6IokqU/cCBw8qQkxWdkANu3Ay1ayDIqf/0liwWePm2/rIJe73qQV6wY8Prrjo8xGoEHHnDt+j7EoImIiHxu2jRXz9Rh9bV78Na9Kz3ZnLxCQiSwsNcj9dprUkrg2DFg3TpZDLd4cQlmHn5YAhFLIrY9w4e7Ps1fbY/QhAlSw8me11+XGYM//5x3lmGFCnL+0qXyXvfvd62WU8uWkphesqTtQqM6nQRrL76o/do+xpwmIiLyi+bNgd9/d+1cAzJx5vfjiGtmZ2aap9y8KbO6vv1Wvq5XT3pt3J31duWK5CS5WhG8aVNg40bnx6lZMNAenU4qlg4dmr1t717gkUekgKYWQUESTK5fL7MTDQbJn7KUP/jmG+vkdB9TG3swaCIiIr8YP146bFyhgwmTWqzAiHWdPNsoX1m5Emjf3vXzP/sMePppx8fcvg2Ehbl+D51OhgCPHrXu2VIUCdiWLpVhTC3X69FD1pxbvlySvuvXlx45P8/eUxt75L+FX4iIqECwt7qIGjqYceZcPs4w2b/fvfO/+EJ6qo4ckZIJnTsDpUt7pm0WiiLDj6mp1vlMOh3QrJkshKz1et9/L1XLH33Us231kXz8G0dERPnZpUuun2uGEWGhATVQoo3W4a3c/vpLFtl9+WWp1VS2rFTWzrlkS2ioDM+5Q6fLntWmKFL8c9o0uefIkdqvZzAAX32l7tijR+U+xYpJG6pUkVpOWhdN9iAGTURE5BeVKtnOC1brh4PVkJHqYk6Qv3lqnTdLXSVFkd6nfv2s97/wguvX1uvlh9S1q6xtV6oUcO+9smrzjBmuXdNkApYtAy5ccHzc5s1A7doyAzE5WXK/Dh+WBQ4TE6X3yw8YNBERkV8MGODOcmc6HM0ohyVv7PRkk3zH3tIq7lAUqQC+b1/2thEjZFFce6Ki7C9obDZLoLJsGbB2LXD+fPZ93HHsmOQynTtne39GBtCtm/Qo5UyUVxR5bdoEvPmme21wEYMmIiLyi8REmVDlOjO+WxD4i7za5KwcgTssw183b0r59T178h4THi7J2Dt2OF7qxRsURWbQ9e5te//PP0tPlL2I2myWnq6cQ5E+wqCJiIj8QqeTmeb9+7t6BT3OpbhYHNKffvsNaNfOe9ffvVv+bd/efsL5rVsSeCQkSBK5NyqBO7NhgzyL3LZtc14d/No14MQJrzTLEQZNRETkN0ajpMi46uBtFVWxA8nMmUCbNsA//3jvHiEhkmjurI7Tzz8Dhw5Jm/xVfeidd/JuCw5W1x5PrLmnEYMmIiLyqxo1ZH1cV1xND/fXGq/aJCVJbaaBA71/r3btgEmT1B27bJlfZ6Nh7dq8iwR36OC46KdOB9x1l6yl52MMmoiIyK+CgoBOLteo1KF1a0+2xk2rVgGNG8ubMhhkGZYyZaSmkjvFLNUKCgIee0z90JXrmfiO6VWGF2azJH7n1LChPEN7UysVBRg1yi9DigyaiIjI75o1c/3cjRulc+LyZc+1xyVjx0ovz5Yt0lNiNgMpKcCZM567h7NA4a23ZH24hAR11/NWIBcSAlStqu7Y3M9HpwMWLwaqVZPvLcGTZZbf6NHuJMK5hUETERH5nbsz8FeuBNq2hf+G6tavB8aN8+49dDpJABs+HChaVLZZenTCwoCJE7MLTtqbmZZTbCxQs6YsleJpqanWa9Y58tFHebfFxsrMvh9+ALp3l2B00CBZ+278eP8kroNrzxERUQBQFAmc/v3XvevMnw/06uWZNmnStSuwZIl37zFuHPDGG/J1Whrwyy/AyZOSEPbAA0DOz0yzWcoNWGbS2bJsGdCxo1QW/+knz7e3USP1S61Mniw1pfyEC/YSEVG+sn490KqVexO52rcHfv3VY01SLzrae+ODhQsDEyYAQ4Zo62FJSpJkse3brbcHBQFz5gB9+sj3FStKwUlvCA5W3/1nCeL8QG3sweE5IiIKCC1aSEqOO9audTzxymvsVdX2hC++kGE5rUNSMTHA1q3AunWy/MjAgcAnn8iif5aACciu9O0NPXqoP3byZO+1w0PY00RERAFj4UKgZ0/3rlG+vKz1qnYCl0c8/bSsk+ZpkZFSHTs01PPXtoiLAy5e9M61T55Un5QOyEw6bwagdrCniYiI8p0mTdy/xvHjwKJF7l9Hk+ef906UNmGCdwMmQF3SuCsKF9a+REtg9ePkwaCJiIgCRpkysri9u95+2/1raFKjhmShe2pWV+HCMpQ2aJBnrufIc8/J7DtPS08Hrl61TlB3JCrK+fIpfsagiYiIAsry5VJqyB3eymt26JFHpJurRw8pCWCr58kSVA0YIMuYtG8vAYvRKG86MRGYN0+SuAcP9k27y5eXopyRkY6P0+mAYsWkZ0pNkJWeLhGw2gKaTZuqO86PmNNEREQB58oV6bxxNUc5Ls67+c2qmM3AV18BH34I/P23BEbt28vU+vvu83PjbEhNBb77TkoZnDsns+oef1ymNJrN2blG337rnSG9AwfUF8T0MJYcICKifO3gweyi0Fq1agX89ptn2+MWs1l6avxUlNGjBg2SpHdPTlOsXh3Yt89z19OIieBERJSvVa3qem7Sjh1+Kj1gj15fMAImwDvvo21bz1/TCxg0ERFRwHr1VVlJIyRE23nXrgF//umdNt3xWrTwfER6/bpnr+clDJqIiCigde8uS6RpNXBggPU2FRQPPgiULu3ZEgtqk8X9jEETEREFPFc+Uw8dAsaMcfGG6enAxx9LwvMzz8jFCoKLF4FJk4Ann5Qq4Rs3aq+NFBQkUxyLF/dcu2JjPXctL2IiOBERBbyePaVauFaFCklB7cKFNZw0e7YESrm7qWrVkgVovVHTyBc+/VRqMpnN2b1EmZlA8+ay2HDRotqud/ky0KCBlFlw16OPykxDP2EiOBERFRjdu7t23s2bwB9/aDhh8WLgqadsj+vt2SM1jQKrr0GdpUuBZ5+V92U2y7+W9/jHH9rWiLOIipJlTzxh8WLPXMfLGDQREVHA69jR9WLRt29rOHjIEMf7L14E3nnHtYb401tv2Z/1ZjZLfQZXMuc91et261a+CEYZNBERUcArVEhGzFzxv/+pPPDcOXk5M2mSaw3xlwsXgO3bnQclAweqv2ZyMvDKK8CNG241Lctdd+WLkgwMmoiIKF8YP1778iqlSwMVKqg8+PJldcclJwNpadoa4k+3bqk7bv9+4MQJ58fNnStDc+++67my675aMsZNDJqIiChfiIjQPjJWo4aGg8uUUX+spjE/PytVCggPd36cTiez4hz580+gXz8Z0vOUGjVc70b0MQZNRESUb8THazteU6mCokWBEiWcH1e8OJCfZneHhAD33+/8OL3eeTDoLOdLq8qVga1btVcv9RMGTURElG9Y1oxVq3ZtjTdQM+19+PB8kX9jZeJE58eYTI4TwBQF2L3bc20CgBkz8lUJBwZNRESUbzRurG0W3fDhGm/Qrh0wcqT9/ffcA7zwgsaLBoAKFYCHH7Yf7BkMQMWKQMuW9q+RmurZGW5Go/xA8xEGTURElG9ERwOPP67u2EKFtA/nAZAE55UrrbupoqKkvPi6derygwLRxx9L8GQwWG83GqW3Z8ECx0ujhIXlPdcdpUrlq14mgEETERHlMx9+qC4YeuQRN27Sti2wa5es/HvxorzGjs2/ARMAxMRI/tDIkRIEAkBoqEShO3cC9es7Pl+nAzp39lx7Tp0CkpI8dz0f4DIqRESU7+zeDdSp43i0aPduWfmEbFAUSfoOCdG28O65c0ClSurLGDizfz9QrZpnruUGLqNCREQF1v/+B3z5pXze50wONxqlQ2T27AIQMO3YAQwaBHToADz2GPDrr56b6q/TydCYloAJAEqWBPbuBapXd78Nej0QF+f+dXyIPU1ERJRv7d0LTJsGrF4tnSetW8useNVVwANJZiawbBnw44/A2rXAmTPZ+/R6CZiaNQN++cU7JQ9u3QI2bwbS0+UBlizp+PjDh6U7b+VKyYe6eVP9vfR64MEHgUWL3Gqyp6iNPRg0ERER+dv585JHtW+f4+O8EWxkZgJvvglMnQpcv559n65dJSKNjXV+jbQ0iVx37wZee8358RERkl9VpYpbTfcUBk1ERET5gdkM1KwJHDyobkq/TgccOSIlAtylKED//sC8eXnvbTAA5coB27YBxYqpv+YLLwAffGB/f6NGwKxZGsu1exdzmoiIiPKDefOAAwfU10BSFGDFCs/ce9s2WUvO1r1NJlmL7qOPtF3z/fel1yomJnubwQC0aiXDf5s3B1TApIWmoGnGjBmoVasWihQpgiJFiqBx48b49ddfs/b369cPOp3O6tWoUSOPN5qIiKjAeP117eekp3vm3nPmOC6zbjIBM2dqu6ZOBwwdCpw9C2zfLuvVJSUBv/0mvUz5mKaC9KVLl8Y777yDu+66CwAwd+5cPPDAA9i1axeq/38mffv27TFnzpysc4KDgz3YXCIiogLk6FHrhG+1PJW+cuqU5DQ5cv689ERpXTrGaATq1XO9bQFIU9DUOVdRq/Hjx2PGjBnYsmVLVtAUEhKCuHw2hZCIiMgvDh507bznn5eilFrWlLElJkaCG0eBU/Hi+W+tPS9xOafJZDJhwYIFuHnzJhrnWDtm/fr1iImJQeXKlTFgwAAk5bNqn0RERD5TqJBr5926JcWo3PXYY44DJoNBEsUJgAuz5/bu3YvGjRvj9u3bKFy4MObPn4+OHTsCABYuXIjChQsjISEBx48fx+uvv47MzEzs2LEDISEhNq+XlpaGtLS0rO9TUlJQpkwZzp4jIqKCLy1N1oS5elX7ufXrSyK3OxRFimeuXp23cKbBIL1Mf//tvGZTPue1kgPp6ek4deoUkpOTsWjRInz++efYsGEDqtkog37+/HkkJCRgwYIF6Natm83rjR07FuPGjcuznUETERHdEd57D3jpJe3nRUQAKSnu3//WLeDZZ4FvvpHEb4v69WVb5cru3yPA+axOU5s2bVCxYkV89tlnNvdXqlQJTz31FF5++WWb+9nTREREdzRFAUaNAiZNUl92wOLgQc8ViDx3Tnqc0tMlgbtuXc9cNx9QGzRpSgS3RVEUq6Anp8uXL+P06dOId7AcdUhIiN2hOyIiogJPpwPeeQd45hng44+BDz+07vFxZPNmzwVNJUsCfft65loFlKagafTo0ejQoQPKlCmD69evY8GCBVi/fj1WrFiBGzduYOzYsejevTvi4+Nx4sQJjB49GtHR0ejatau32k9ERFQwlCsnhSHLlweee07dOR99BBQuLEuruDuTjpzSFDRdvHgRjz32GM6fP4/IyEjUqlULK1asQGJiIlJTU7F3717MmzcPycnJiI+PR8uWLbFw4UJERER4q/1EREQFS9Om6o/9+2/g4YelbtOrrwIjRsi6ceQVXHuOiIgokJw9C5Qu7dq5Tz0lFbxZV0kTn+U0ERERkQeVKgXcfTdw6JD2cz//XOoqpacDixYB168DVasC/foBsbEeb+qdhj1NREREgWbqVGD4cO3nGY1SWykpKXtNObNZhuw++QQYONCjzSwo1MYeHPgkIiIKNM8953ghXXsyMyVgsnydmSlBU2amzM5r2VIW0SWXMGgiIiIKNAYD8MILnr/u+vVAgwaSM/XWW8CFC56/RwHGoImIiCgQvfsu0KuXd6599iwwZgxQsaIEUqQKgyYiIqJANX++JIQ/+KDnE7kVRZZQad8euHjRs9cuoBg0ERERBbLKlYHFi50XvIyMdO36aWlAo0Yy044cYtBEREQU6G7fBiZPdnzMjRuuX//ECSlLQA4xaCIiIgp0f/0FXLvm+BiTCahQwfV7/PgjcPSo6+ffARg0ERERBbq0NHXHPfCA6/fQ64Hly10/P7e0NGDVKuD772W5lwKAQRMREVGgq1FD3ZpyjzwCvPmma/dQFODqVfl6+3bg6aeBZs2ALl2Ar79WH7gpiiwkHB8PtGsna+PVqQPUrQvs2OFa2wIEK4ITERHlB926AUuXyjBcbgYDUL269OjodMA77wCjRmm/R1wc0LkzMGuWFNfMzJRgzWyWpV3WrgVKlnR8jYkTgdGjbbcxOFiGGmvW1N42L1IbezBoIiIiyg/OnwcaNwbOnLEOnIxGICIC+OMPCZwACXJq1QL279d2D0uAZIvRCNSuDWzdan9B4P/+k6AqI8P2foMB6NAB+Plnbe3yMi6jQkREVJDEx8uw2csvA9HRsq1wYVkeZdeu7IAJkODn118BrZ0P9gImQHqdtm8HNm2yf8yCBbZ7wixMJmDZMuDSJW3tChAMmoiIiPKL6Ghg/HgJOtLSgJQU4OOPgYSEvMeWKQPs25cdYHmC0QisXGl//7lz0pvkiKLk22KaDJqIiIjyo+Bg+8NkFmXKAAcOAM8/D4SHu39PnU56nOyJi3Pc02S5RkyM+23xAwZNREREBVmJEsCHHwLJycDhw0BYmOvXysgAGja0v79nT8ez/AwGmVHHoImIiIgCVlAQUKkSMHCguvIFtoSGAp062d8fEwO8+qrtfXq9tGHCBNfuHQAYNBEREd1J3nlHensA5/lHuaWnA1euOD5mzBhg0qS8Seh33w2sWyc1m/IplhwgIiK601hmsX3+ObBli7bZbJs3ywK/zty6BaxZI8nqlSoB99zjPAfLT9TGHkYftomIiIgCgcEglb67dJG6SStWqD83IkLdceHhcv0ChMNzREREd7KwMHU5Tjqd9BhVq+b9NgUoBk1ERER3si5dHBe1tFAUWdcuQIfYfIFBExER0Z3skUeA0qUdJ4UbjcAnn0hJgTsYgyYiIqI7WVgY8NtvQKlS8r3RKMN1Op0EUk89BVy4AAwe7N92BgAmghMREd3pKlcGjhwBfvxR1qxLTwcaNAD69gWiovzduoDBkgNERER0R1Mbe3B4joiIiEgFBk1EREREKjBoIiIiIlKBQRMRERGRCgyaiIiIiFRg0ERERESkAoMmIiIiIhUYNBERERGpwKCJiIiISAUGTUREREQqMGgiIiIiUoFBExEREZEKDJqIiIiIVGDQRERERKQCgyYiIiIiFRg0EREREanAoImIiIhIBaO/G5CboigAgJSUFD+3hIiIiO4ElpjDEoPYE3BB0/Xr1wEAZcqU8XNLiIiI6E5y/fp1REZG2t2vU5yFVT5mNptx7tw5REREQKfT+bs5XpeSkoIyZcrg9OnTKFKkiL+bky/xGXoGn6P7+Aw9g8/RfXyG2iiKguvXr6NkyZLQ6+1nLgVcT5Ner0fp0qX93QyfK1KkCH+x3cRn6Bl8ju7jM/QMPkf38Rmq56iHyYKJ4EREREQqMGgiIiIiUoFBk5+FhIRgzJgxCAkJ8XdT8i0+Q8/gc3Qfn6Fn8Dm6j8/QOwIuEZyIiIgoELGniYiIiEgFBk1EREREKjBoIiIiIlKBQRMRERGRCgyafGT8+PFo0qQJwsPDUbRoUbvHffnll6hVqxZCQ0MRFxeHIUOGWO3fu3cvmjdvjrCwMJQqVQpvvvmm07VyChI1z1Gn0+V5ffrpp1bH3MnPUe3vIgBcvnwZpUuXhk6nQ3JystW+O/kZAs6f4+XLl9G+fXuULFkSISEhKFOmDIYMGZJnXc07+Tk6e4a7d+9Gr169UKZMGYSFhaFq1ar48MMP8xzHZ+j4v+ehQ4eiXr16CAkJQe3atW0ecyc/Qy0CriJ4QZWeno6HHnoIjRs3xuzZs20eM2XKFLz//vuYPHkyGjZsiNu3b+PYsWNZ+1NSUpCYmIiWLVti27ZtOHz4MPr164dChQrhxRdf9NVb8Ss1zxEA5syZg/bt22d9n7PS653+HNU+QwB48sknUatWLZw9e9Zq+53+DAHnz1Gv1+OBBx7A22+/jRIlSuDo0aMYPHgwrly5gvnz5wPgc3T2DHfs2IESJUrg66+/RpkyZbBp0yY8/fTTMBgMWX9Q8hk6/+9ZURQ88cQT+Ouvv7Bnz548++/0Z6iJQj41Z84cJTIyMs/2K1euKGFhYcqaNWvsnjt9+nQlMjJSuX37dta2iRMnKiVLllTMZrM3mhuw7D1HRVEUAMrixYvtnsvnKBw9Q0WR59S8eXPlt99+UwAoV69etdrHZyicPcecPvzwQ6V06dJZ3/M5Ci3PcNCgQUrLli2zvuczFGqe4ZgxY5T//e9/ebbzGarH4bkAsXr1apjNZpw9exZVq1ZF6dKl8fDDD+P06dNZx2zevBnNmze3KlbWrl07nDt3DidOnPBDqwPXkCFDEB0djQYNGuDTTz+F2WzO2sfn6NyBAwfw5ptvYt68eTYXr+Qz1O7cuXP48ccf0bx586xtfI7aXbt2DcWLF8/6ns/QfXyG6jFoChDHjh2D2WzGhAkTMHXqVPzwww+4cuUKEhMTkZ6eDgC4cOECYmNjrc6zfH/hwgWftzlQvfXWW/j++++xZs0a9OzZEy+++CImTJiQtZ/P0bG0tDT06tULkydPRtmyZW0ew2eoXq9evRAeHo5SpUqhSJEi+Pzzz7P28Tlqs3nzZnz33XcYOHBg1jY+Q/fxGarHoMkNY8eOtZl0nPO1fft2Vdcym83IyMjARx99hHbt2qFRo0b49ttvceTIEaxbty7rOJ1OZ3We8v+Jerm35yeefI4A8Nprr6Fx48aoXbs2XnzxRbz55puYPHmy1TEF7Tl68hmOGjUKVatWxaOPPurwuIL2DAHP/y4CwAcffICdO3diyZIl+Pfff/HCCy9Y7S9oz9EbzxAA9u/fjwceeABvvPEGEhMTrfbxGbqvoD1Db2EiuBuGDBmCnj17OjymXLlyqq4VHx8PAKhWrVrWthIlSiA6OhqnTp0CAMTFxeWJ+pOSkgAgz18J+Yknn6MtjRo1QkpKCi5evIjY2NgC+Rw9+QzXrl2LvXv34ocffgCQ/T/P6OhovPrqqxg3blyBfIaAd34X4+LiEBcXhypVqiAqKgrNmjXD66+/jvj4+AL5HL3xDA8cOIBWrVphwIABeO2116z28Rm6ryA+Q29h0OSG6OhoREdHe+RaTZs2BQAcOnQIpUuXBgBcuXIF//33HxISEgAAjRs3xujRo5Geno7g4GAAwKpVq1CyZEmP/gfka558jrbs2rULoaGhWdNxC+Jz9OQzXLRoEVJTU7O+37ZtG5544gn88ccfqFixIoCC+QwB7/8uWgLQtLQ0AAXzOXr6Ge7fvx+tWrVC3759MX78+Dz7+QzdVxCfodf4Lwf9znLy5Ell165dyrhx45TChQsru3btUnbt2qVcv34965j/a+eOXY4L4zCOe6dHElJiMzFgsEryDygZZVAWm9XIH2AwU/4EgzrZnMmsHBulLMrEiOF6txM9es/9LM87nO+nznTufnWu7ruu5T71el35fF7r9VqO46hWqymXy+nxeEiSrterksmkms2mHMfRfD5XJBLRaDT6X5/167xyXCwWmkwmchxHh8NB0+lUkUhEvV7PneH3HE324ivbtr/dnvN7hpJ3jpZlaTabyXEcHY9HWZalfD6vcrnszvB7jl4Z7nY7JRIJtVotnc9n97lcLu4MMvQ+z/v9XpvNRt1uV9ls1l1zv98lkeFPUJp+SbvdViAQ+PbYtu2uud1u6nQ6isViisfjajQaOp1Ob3O2260qlYq+vr6USqU0HA59dSXUK8flcqlisahwOKxQKKRCoaDxeKzn8/k2x885muzFV59Kk+TvDCXvHFerlUqlkqLRqILBoDKZjPr9Pjm+8MpwMBh8fJ9Op9/mkOG/z3O1Wv245ng8umv8nOFP/JH45ScAAIAXbs8BAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAYoDQBAAAY+AuMdLGP7qVqwAAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 700x700 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "f,ax=plt.subplots(figsize=(7,7))\n",
    "all_stations_gdf.plot(ax=ax,column='network',cmap='bwr',legend=True)\n",
    "ax.set_title('stations by network')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "id": "f8863b10-c1e7-4af4-a024-400cd7f9d984",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "current_files = [os.path.splitext(os.path.basename(path))[0] for path in glob.glob('data/*.csv')]\n",
    "idx_csv = all_stations_gdf.index.isin(current_files)\n",
    "all_stations_gdf['csvData'] = False\n",
    "all_stations_gdf.loc[idx_csv,'csvData'] = True"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "37f7a420-6e0d-40e1-abc2-88025024d5dc",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "all_stations_gdf = all_stations_gdf[['name','network','elevation_m','latitude','longitude','state','HUC','mgrs','mountainRange','beginDate','endDate','csvData','geometry']]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "21c93281-4a9a-45fa-a95a-86bf7ead3bd1",
   "metadata": {},
   "outputs": [],
   "source": [
    "all_stations_gdf.to_file('all_stations.geojson')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9e342821-78d1-4233-822b-1e4d5d7c76a4",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "snotel_ccss_env",
   "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.18"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}