00_data_source_exploration/weather_station.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# NOAA Historical Weather Station Data\n",
    "https://registry.opendata.aws/noaa-isd/\n",
    "\n",
    "The Integrated Surface Database (ISD) from the United States National Oceanic and Atmospheric Administration (NOAA) consists of global hourly and synoptic observations compiled from numerous sources into a gzipped fixed width format. The data format is described [here](https://www.ncei.noaa.gov/data/global-hourly/doc/isd-format-document.pdf). They also have conveniently converted all of this gzipped, fixed-width data into csv files and host it on an open S3 bucket (`arn:aws:s3:::noaa-global-hourly-pds`). You can [browse the bucket here](https://noaa-global-hourly-pds.s3.amazonaws.com/index.html). The data is also on this [ftp site](https://www1.ncdc.noaa.gov/pub/data/noaa/). This site contains other useful files as well like a [csv file with all the station IDs, locations, and record date ranges](https://www1.ncdc.noaa.gov/pub/data/noaa/isd-history.csv).\n",
    "\n",
    "In this notebook I will explore the data source by extracting the temperature and dew point data recorded at Centennial Airport in Colorado during 2021."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "First, download the csv with station metadata."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "--2022-06-29 03:26:20--  https://noaa-isd-pds.s3.amazonaws.com/isd-history.csv\n",
      "Resolving noaa-isd-pds.s3.amazonaws.com (noaa-isd-pds.s3.amazonaws.com)... 52.216.179.139\n",
      "Connecting to noaa-isd-pds.s3.amazonaws.com (noaa-isd-pds.s3.amazonaws.com)|52.216.179.139|:443... connected.\n",
      "HTTP request sent, awaiting response... 200 OK\n",
      "Length: 2902758 (2.8M) [binary/octet-stream]\n",
      "Saving to: ‘isd-stations.csv’\n",
      "\n",
      "isd-stations.csv    100%[===================>]   2.77M  13.4MB/s    in 0.2s    \n",
      "\n",
      "2022-06-29 03:26:21 (13.4 MB/s) - ‘isd-stations.csv’ saved [2902758/2902758]\n",
      "\n"
     ]
    }
   ],
   "source": [
    "!wget https://noaa-isd-pds.s3.amazonaws.com/isd-history.csv -O isd-stations.csv"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Read in the data from the csv, while parsing the correct data types, and merging the `USAF` and `WBAN` together for a single, combined ID"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<class 'pandas.core.frame.DataFrame'>\n",
      "RangeIndex: 29615 entries, 0 to 29614\n",
      "Data columns (total 12 columns):\n",
      " #   Column        Non-Null Count  Dtype         \n",
      "---  ------        --------------  -----         \n",
      " 0   USAF          29615 non-null  object        \n",
      " 1   WBAN          29615 non-null  object        \n",
      " 2   STATION NAME  28825 non-null  object        \n",
      " 3   CTRY          28720 non-null  object        \n",
      " 4   STATE         6622 non-null   object        \n",
      " 5   ICAO          10803 non-null  object        \n",
      " 6   LAT           28429 non-null  float64       \n",
      " 7   LON           28428 non-null  float64       \n",
      " 8   ELEV(M)       28325 non-null  float64       \n",
      " 9   BEGIN         29615 non-null  datetime64[ns]\n",
      " 10  END           29615 non-null  datetime64[ns]\n",
      " 11  station_id    29615 non-null  object        \n",
      "dtypes: datetime64[ns](2), float64(3), object(7)\n",
      "memory usage: 2.7+ MB\n"
     ]
    },
    {
     "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>USAF</th>\n",
       "      <th>WBAN</th>\n",
       "      <th>STATION NAME</th>\n",
       "      <th>CTRY</th>\n",
       "      <th>STATE</th>\n",
       "      <th>ICAO</th>\n",
       "      <th>LAT</th>\n",
       "      <th>LON</th>\n",
       "      <th>ELEV(M)</th>\n",
       "      <th>BEGIN</th>\n",
       "      <th>END</th>\n",
       "      <th>station_id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>007018</td>\n",
       "      <td>99999</td>\n",
       "      <td>WXPOD 7018</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>7018.0</td>\n",
       "      <td>2011-03-09</td>\n",
       "      <td>2013-07-30</td>\n",
       "      <td>00701899999</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>007026</td>\n",
       "      <td>99999</td>\n",
       "      <td>WXPOD 7026</td>\n",
       "      <td>AF</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>7026.0</td>\n",
       "      <td>2012-07-13</td>\n",
       "      <td>2017-08-22</td>\n",
       "      <td>00702699999</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>007070</td>\n",
       "      <td>99999</td>\n",
       "      <td>WXPOD 7070</td>\n",
       "      <td>AF</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>7070.0</td>\n",
       "      <td>2014-09-23</td>\n",
       "      <td>2015-09-26</td>\n",
       "      <td>00707099999</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>008260</td>\n",
       "      <td>99999</td>\n",
       "      <td>WXPOD8270</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2005-01-01</td>\n",
       "      <td>2012-07-31</td>\n",
       "      <td>00826099999</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>008268</td>\n",
       "      <td>99999</td>\n",
       "      <td>WXPOD8278</td>\n",
       "      <td>AF</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>32.95</td>\n",
       "      <td>65.567</td>\n",
       "      <td>1156.7</td>\n",
       "      <td>2010-05-19</td>\n",
       "      <td>2012-03-23</td>\n",
       "      <td>00826899999</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "     USAF   WBAN STATION NAME CTRY STATE ICAO    LAT     LON  ELEV(M)  \\\n",
       "0  007018  99999   WXPOD 7018  NaN   NaN  NaN   0.00   0.000   7018.0   \n",
       "1  007026  99999   WXPOD 7026   AF   NaN  NaN   0.00   0.000   7026.0   \n",
       "2  007070  99999   WXPOD 7070   AF   NaN  NaN   0.00   0.000   7070.0   \n",
       "3  008260  99999    WXPOD8270  NaN   NaN  NaN   0.00   0.000      0.0   \n",
       "4  008268  99999    WXPOD8278   AF   NaN  NaN  32.95  65.567   1156.7   \n",
       "\n",
       "       BEGIN        END   station_id  \n",
       "0 2011-03-09 2013-07-30  00701899999  \n",
       "1 2012-07-13 2017-08-22  00702699999  \n",
       "2 2014-09-23 2015-09-26  00707099999  \n",
       "3 2005-01-01 2012-07-31  00826099999  \n",
       "4 2010-05-19 2012-03-23  00826899999  "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "isd_stations = (pd.read_csv('isd-stations.csv',\n",
    "                            parse_dates=['BEGIN', 'END'], \n",
    "                            dtype={'USAF': str, 'WBAN': str}\n",
    "                            )\n",
    "                .assign(station_id=lambda df_: df_['USAF'] + df_['WBAN'])\n",
    ")\n",
    "\n",
    "isd_stations.info()\n",
    "isd_stations.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Filter the data to all stations in Colorado that have data from 2022 and have the word `AIRPORT` in the name."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "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>USAF</th>\n",
       "      <th>WBAN</th>\n",
       "      <th>STATION NAME</th>\n",
       "      <th>CTRY</th>\n",
       "      <th>STATE</th>\n",
       "      <th>ICAO</th>\n",
       "      <th>LAT</th>\n",
       "      <th>LON</th>\n",
       "      <th>ELEV(M)</th>\n",
       "      <th>BEGIN</th>\n",
       "      <th>END</th>\n",
       "      <th>station_id</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>19569</th>\n",
       "      <td>724620</td>\n",
       "      <td>23061</td>\n",
       "      <td>SAN LUIS VALLEY REGIONAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KALS</td>\n",
       "      <td>37.439</td>\n",
       "      <td>-105.862</td>\n",
       "      <td>2297.4</td>\n",
       "      <td>1973-01-01</td>\n",
       "      <td>2022-06-28</td>\n",
       "      <td>72462023061</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17465</th>\n",
       "      <td>720539</td>\n",
       "      <td>00165</td>\n",
       "      <td>STEVENS FIELD AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KPSO</td>\n",
       "      <td>37.283</td>\n",
       "      <td>-107.050</td>\n",
       "      <td>2335.1</td>\n",
       "      <td>2009-10-01</td>\n",
       "      <td>2022-06-28</td>\n",
       "      <td>72053900165</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17999</th>\n",
       "      <td>722101</td>\n",
       "      <td>03039</td>\n",
       "      <td>WOLF CREEK PASS AWOS-3 AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KCPW</td>\n",
       "      <td>37.450</td>\n",
       "      <td>-106.800</td>\n",
       "      <td>3593.9</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-28</td>\n",
       "      <td>72210103039</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17456</th>\n",
       "      <td>720531</td>\n",
       "      <td>00158</td>\n",
       "      <td>BLAKE FIELD AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KAJZ</td>\n",
       "      <td>38.783</td>\n",
       "      <td>-108.067</td>\n",
       "      <td>1583.1</td>\n",
       "      <td>2012-02-09</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72053100158</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17458</th>\n",
       "      <td>720533</td>\n",
       "      <td>00160</td>\n",
       "      <td>BOULDER MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KBDU</td>\n",
       "      <td>40.033</td>\n",
       "      <td>-105.217</td>\n",
       "      <td>1612.1</td>\n",
       "      <td>2010-02-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72053300160</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19607</th>\n",
       "      <td>724680</td>\n",
       "      <td>94015</td>\n",
       "      <td>BUTTS AAF AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KFCS</td>\n",
       "      <td>38.678</td>\n",
       "      <td>-104.757</td>\n",
       "      <td>1779.4</td>\n",
       "      <td>1966-09-15</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72468094015</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19590</th>\n",
       "      <td>724666</td>\n",
       "      <td>93067</td>\n",
       "      <td>CENTENNIAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KAPA</td>\n",
       "      <td>39.560</td>\n",
       "      <td>-104.848</td>\n",
       "      <td>1789.2</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72466693067</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19655</th>\n",
       "      <td>724767</td>\n",
       "      <td>93069</td>\n",
       "      <td>CORTEZ MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KCEZ</td>\n",
       "      <td>37.307</td>\n",
       "      <td>-108.626</td>\n",
       "      <td>1799.2</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72476793069</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20276</th>\n",
       "      <td>725700</td>\n",
       "      <td>24046</td>\n",
       "      <td>CRAIG-MOFFAT AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KCAG</td>\n",
       "      <td>40.493</td>\n",
       "      <td>-107.524</td>\n",
       "      <td>1884.9</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72570024046</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20265</th>\n",
       "      <td>725650</td>\n",
       "      <td>03017</td>\n",
       "      <td>DENVER INTERNATIONAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KDEN</td>\n",
       "      <td>39.847</td>\n",
       "      <td>-104.656</td>\n",
       "      <td>1647.2</td>\n",
       "      <td>1994-07-18</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72565003017</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19570</th>\n",
       "      <td>724625</td>\n",
       "      <td>93005</td>\n",
       "      <td>DURA-LA PLATA CO AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KDRO</td>\n",
       "      <td>37.143</td>\n",
       "      <td>-107.760</td>\n",
       "      <td>2020.6</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72462593005</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19598</th>\n",
       "      <td>724675</td>\n",
       "      <td>23063</td>\n",
       "      <td>EAGLE CO REGIONAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KEGE</td>\n",
       "      <td>39.650</td>\n",
       "      <td>-106.917</td>\n",
       "      <td>1980.3</td>\n",
       "      <td>1973-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72467523063</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17459</th>\n",
       "      <td>720534</td>\n",
       "      <td>00161</td>\n",
       "      <td>ERIE MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KEIK</td>\n",
       "      <td>40.017</td>\n",
       "      <td>-105.050</td>\n",
       "      <td>1564.2</td>\n",
       "      <td>2014-07-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72053400161</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17461</th>\n",
       "      <td>720535</td>\n",
       "      <td>00162</td>\n",
       "      <td>FORT MORGAN MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KFMM</td>\n",
       "      <td>40.333</td>\n",
       "      <td>-103.800</td>\n",
       "      <td>1393.2</td>\n",
       "      <td>2012-07-26</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72053500162</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19651</th>\n",
       "      <td>724760</td>\n",
       "      <td>23066</td>\n",
       "      <td>GRAND JUNCTION REGIONAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KGJT</td>\n",
       "      <td>39.134</td>\n",
       "      <td>-108.541</td>\n",
       "      <td>1470.4</td>\n",
       "      <td>1973-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72476023066</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17457</th>\n",
       "      <td>720532</td>\n",
       "      <td>00159</td>\n",
       "      <td>HARRIET ALEXANDER FIELD AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KANK</td>\n",
       "      <td>38.533</td>\n",
       "      <td>-106.050</td>\n",
       "      <td>2294.2</td>\n",
       "      <td>2009-11-21</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72053200159</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17463</th>\n",
       "      <td>720537</td>\n",
       "      <td>00163</td>\n",
       "      <td>HOLYOKE AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KHEQ</td>\n",
       "      <td>40.567</td>\n",
       "      <td>-102.267</td>\n",
       "      <td>1137.2</td>\n",
       "      <td>2010-11-29</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72053700163</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19612</th>\n",
       "      <td>724689</td>\n",
       "      <td>03026</td>\n",
       "      <td>KIT CARSON COUNTY AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KITR</td>\n",
       "      <td>39.241</td>\n",
       "      <td>-102.282</td>\n",
       "      <td>1278.2</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72468903026</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19577</th>\n",
       "      <td>724635</td>\n",
       "      <td>23067</td>\n",
       "      <td>LA JUNTA MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KLHX</td>\n",
       "      <td>38.049</td>\n",
       "      <td>-103.513</td>\n",
       "      <td>1278.7</td>\n",
       "      <td>1942-11-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72463523067</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19578</th>\n",
       "      <td>724636</td>\n",
       "      <td>03013</td>\n",
       "      <td>LAMAR MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KLAA</td>\n",
       "      <td>38.072</td>\n",
       "      <td>-102.687</td>\n",
       "      <td>1124.0</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72463603013</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19588</th>\n",
       "      <td>724665</td>\n",
       "      <td>93010</td>\n",
       "      <td>LIMON MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KLIC</td>\n",
       "      <td>39.273</td>\n",
       "      <td>-103.667</td>\n",
       "      <td>1633.1</td>\n",
       "      <td>1973-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72466593010</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29577</th>\n",
       "      <td>A00005</td>\n",
       "      <td>94076</td>\n",
       "      <td>MC ELORY AIRFIELD AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>K20V</td>\n",
       "      <td>40.054</td>\n",
       "      <td>-106.369</td>\n",
       "      <td>2258.9</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>A0000594076</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17675</th>\n",
       "      <td>720852</td>\n",
       "      <td>00280</td>\n",
       "      <td>MEADOW LAKE AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KFLY</td>\n",
       "      <td>38.946</td>\n",
       "      <td>-104.570</td>\n",
       "      <td>2095.2</td>\n",
       "      <td>2010-12-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72085200280</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19596</th>\n",
       "      <td>724674</td>\n",
       "      <td>94050</td>\n",
       "      <td>MEEKER AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KEEO</td>\n",
       "      <td>40.044</td>\n",
       "      <td>-107.888</td>\n",
       "      <td>1940.8</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72467494050</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19652</th>\n",
       "      <td>724765</td>\n",
       "      <td>93013</td>\n",
       "      <td>MONTROSE REGIONAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KMTJ</td>\n",
       "      <td>38.505</td>\n",
       "      <td>-107.898</td>\n",
       "      <td>1739.5</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72476593013</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19581</th>\n",
       "      <td>724645</td>\n",
       "      <td>23070</td>\n",
       "      <td>PERRY STOKES AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KTAD</td>\n",
       "      <td>37.262</td>\n",
       "      <td>-104.332</td>\n",
       "      <td>1749.3</td>\n",
       "      <td>1973-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72464523070</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19580</th>\n",
       "      <td>724640</td>\n",
       "      <td>93058</td>\n",
       "      <td>PUEBLO MEMORIAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KPUB</td>\n",
       "      <td>38.289</td>\n",
       "      <td>-104.506</td>\n",
       "      <td>1424.5</td>\n",
       "      <td>1942-10-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72464093058</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18336</th>\n",
       "      <td>722333</td>\n",
       "      <td>03069</td>\n",
       "      <td>SAGUACHE MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>K04V</td>\n",
       "      <td>38.097</td>\n",
       "      <td>-106.169</td>\n",
       "      <td>2385.4</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72233303069</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19604</th>\n",
       "      <td>724678</td>\n",
       "      <td>94033</td>\n",
       "      <td>STEAMBOAT SPRINGS/BOB ADAMS FIELD AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KSBS</td>\n",
       "      <td>40.517</td>\n",
       "      <td>-106.867</td>\n",
       "      <td>2096.7</td>\n",
       "      <td>2011-02-14</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72467894033</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17472</th>\n",
       "      <td>720544</td>\n",
       "      <td>00168</td>\n",
       "      <td>STERLING MUNICIPAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KSTK</td>\n",
       "      <td>40.617</td>\n",
       "      <td>-103.267</td>\n",
       "      <td>1231.1</td>\n",
       "      <td>2010-11-29</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72054400168</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19573</th>\n",
       "      <td>724627</td>\n",
       "      <td>03011</td>\n",
       "      <td>TELLURIDE REGIONAL AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KTEX</td>\n",
       "      <td>37.950</td>\n",
       "      <td>-107.900</td>\n",
       "      <td>2770.9</td>\n",
       "      <td>2006-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72462703011</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21677</th>\n",
       "      <td>745310</td>\n",
       "      <td>93065</td>\n",
       "      <td>USAF ACADEMY AIRSTRIP AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KAFF</td>\n",
       "      <td>38.967</td>\n",
       "      <td>-104.817</td>\n",
       "      <td>2003.2</td>\n",
       "      <td>1967-11-27</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>74531093065</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17464</th>\n",
       "      <td>720538</td>\n",
       "      <td>00164</td>\n",
       "      <td>VANCE BRAND AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KLMO</td>\n",
       "      <td>40.167</td>\n",
       "      <td>-105.167</td>\n",
       "      <td>1541.1</td>\n",
       "      <td>2010-12-02</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72053800164</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20283</th>\n",
       "      <td>725710</td>\n",
       "      <td>94025</td>\n",
       "      <td>YAMPA VALLEY AIRPORT</td>\n",
       "      <td>US</td>\n",
       "      <td>CO</td>\n",
       "      <td>KHDN</td>\n",
       "      <td>40.481</td>\n",
       "      <td>-107.218</td>\n",
       "      <td>2011.7</td>\n",
       "      <td>2005-01-01</td>\n",
       "      <td>2022-06-27</td>\n",
       "      <td>72571094025</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "         USAF   WBAN                               STATION NAME CTRY STATE  \\\n",
       "19569  724620  23061           SAN LUIS VALLEY REGIONAL AIRPORT   US    CO   \n",
       "17465  720539  00165                      STEVENS FIELD AIRPORT   US    CO   \n",
       "17999  722101  03039             WOLF CREEK PASS AWOS-3 AIRPORT   US    CO   \n",
       "17456  720531  00158                        BLAKE FIELD AIRPORT   US    CO   \n",
       "17458  720533  00160                  BOULDER MUNICIPAL AIRPORT   US    CO   \n",
       "19607  724680  94015                          BUTTS AAF AIRPORT   US    CO   \n",
       "19590  724666  93067                         CENTENNIAL AIRPORT   US    CO   \n",
       "19655  724767  93069                   CORTEZ MUNICIPAL AIRPORT   US    CO   \n",
       "20276  725700  24046                       CRAIG-MOFFAT AIRPORT   US    CO   \n",
       "20265  725650  03017               DENVER INTERNATIONAL AIRPORT   US    CO   \n",
       "19570  724625  93005                   DURA-LA PLATA CO AIRPORT   US    CO   \n",
       "19598  724675  23063                  EAGLE CO REGIONAL AIRPORT   US    CO   \n",
       "17459  720534  00161                     ERIE MUNICIPAL AIRPORT   US    CO   \n",
       "17461  720535  00162              FORT MORGAN MUNICIPAL AIRPORT   US    CO   \n",
       "19651  724760  23066            GRAND JUNCTION REGIONAL AIRPORT   US    CO   \n",
       "17457  720532  00159            HARRIET ALEXANDER FIELD AIRPORT   US    CO   \n",
       "17463  720537  00163                            HOLYOKE AIRPORT   US    CO   \n",
       "19612  724689  03026                  KIT CARSON COUNTY AIRPORT   US    CO   \n",
       "19577  724635  23067                 LA JUNTA MUNICIPAL AIRPORT   US    CO   \n",
       "19578  724636  03013                    LAMAR MUNICIPAL AIRPORT   US    CO   \n",
       "19588  724665  93010                    LIMON MUNICIPAL AIRPORT   US    CO   \n",
       "29577  A00005  94076                  MC ELORY AIRFIELD AIRPORT   US    CO   \n",
       "17675  720852  00280                        MEADOW LAKE AIRPORT   US    CO   \n",
       "19596  724674  94050                             MEEKER AIRPORT   US    CO   \n",
       "19652  724765  93013                  MONTROSE REGIONAL AIRPORT   US    CO   \n",
       "19581  724645  23070                       PERRY STOKES AIRPORT   US    CO   \n",
       "19580  724640  93058                    PUEBLO MEMORIAL AIRPORT   US    CO   \n",
       "18336  722333  03069                 SAGUACHE MUNICIPAL AIRPORT   US    CO   \n",
       "19604  724678  94033  STEAMBOAT SPRINGS/BOB ADAMS FIELD AIRPORT   US    CO   \n",
       "17472  720544  00168                 STERLING MUNICIPAL AIRPORT   US    CO   \n",
       "19573  724627  03011                 TELLURIDE REGIONAL AIRPORT   US    CO   \n",
       "21677  745310  93065              USAF ACADEMY AIRSTRIP AIRPORT   US    CO   \n",
       "17464  720538  00164                        VANCE BRAND AIRPORT   US    CO   \n",
       "20283  725710  94025                       YAMPA VALLEY AIRPORT   US    CO   \n",
       "\n",
       "       ICAO     LAT      LON  ELEV(M)      BEGIN        END   station_id  \n",
       "19569  KALS  37.439 -105.862   2297.4 1973-01-01 2022-06-28  72462023061  \n",
       "17465  KPSO  37.283 -107.050   2335.1 2009-10-01 2022-06-28  72053900165  \n",
       "17999  KCPW  37.450 -106.800   3593.9 2006-01-01 2022-06-28  72210103039  \n",
       "17456  KAJZ  38.783 -108.067   1583.1 2012-02-09 2022-06-27  72053100158  \n",
       "17458  KBDU  40.033 -105.217   1612.1 2010-02-01 2022-06-27  72053300160  \n",
       "19607  KFCS  38.678 -104.757   1779.4 1966-09-15 2022-06-27  72468094015  \n",
       "19590  KAPA  39.560 -104.848   1789.2 2006-01-01 2022-06-27  72466693067  \n",
       "19655  KCEZ  37.307 -108.626   1799.2 2006-01-01 2022-06-27  72476793069  \n",
       "20276  KCAG  40.493 -107.524   1884.9 2006-01-01 2022-06-27  72570024046  \n",
       "20265  KDEN  39.847 -104.656   1647.2 1994-07-18 2022-06-27  72565003017  \n",
       "19570  KDRO  37.143 -107.760   2020.6 2006-01-01 2022-06-27  72462593005  \n",
       "19598  KEGE  39.650 -106.917   1980.3 1973-01-01 2022-06-27  72467523063  \n",
       "17459  KEIK  40.017 -105.050   1564.2 2014-07-01 2022-06-27  72053400161  \n",
       "17461  KFMM  40.333 -103.800   1393.2 2012-07-26 2022-06-27  72053500162  \n",
       "19651  KGJT  39.134 -108.541   1470.4 1973-01-01 2022-06-27  72476023066  \n",
       "17457  KANK  38.533 -106.050   2294.2 2009-11-21 2022-06-27  72053200159  \n",
       "17463  KHEQ  40.567 -102.267   1137.2 2010-11-29 2022-06-27  72053700163  \n",
       "19612  KITR  39.241 -102.282   1278.2 2006-01-01 2022-06-27  72468903026  \n",
       "19577  KLHX  38.049 -103.513   1278.7 1942-11-01 2022-06-27  72463523067  \n",
       "19578  KLAA  38.072 -102.687   1124.0 2006-01-01 2022-06-27  72463603013  \n",
       "19588  KLIC  39.273 -103.667   1633.1 1973-01-01 2022-06-27  72466593010  \n",
       "29577  K20V  40.054 -106.369   2258.9 2006-01-01 2022-06-27  A0000594076  \n",
       "17675  KFLY  38.946 -104.570   2095.2 2010-12-01 2022-06-27  72085200280  \n",
       "19596  KEEO  40.044 -107.888   1940.8 2006-01-01 2022-06-27  72467494050  \n",
       "19652  KMTJ  38.505 -107.898   1739.5 2006-01-01 2022-06-27  72476593013  \n",
       "19581  KTAD  37.262 -104.332   1749.3 1973-01-01 2022-06-27  72464523070  \n",
       "19580  KPUB  38.289 -104.506   1424.5 1942-10-01 2022-06-27  72464093058  \n",
       "18336  K04V  38.097 -106.169   2385.4 2006-01-01 2022-06-27  72233303069  \n",
       "19604  KSBS  40.517 -106.867   2096.7 2011-02-14 2022-06-27  72467894033  \n",
       "17472  KSTK  40.617 -103.267   1231.1 2010-11-29 2022-06-27  72054400168  \n",
       "19573  KTEX  37.950 -107.900   2770.9 2006-01-01 2022-06-27  72462703011  \n",
       "21677  KAFF  38.967 -104.817   2003.2 1967-11-27 2022-06-27  74531093065  \n",
       "17464  KLMO  40.167 -105.167   1541.1 2010-12-02 2022-06-27  72053800164  \n",
       "20283  KHDN  40.481 -107.218   2011.7 2005-01-01 2022-06-27  72571094025  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "(isd_stations[ (isd_stations['STATE']=='CO') & \n",
    "               (isd_stations['END'] > '2022') & \n",
    "               (isd_stations['STATION NAME'].str.contains('AIRPORT'))\n",
    "            ]\n",
    "            .sort_values(by=['END', 'STATION NAME'], ascending=[False, True])\n",
    "            .head(100)\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I want to get the data for the Centennial Airport, so I will use `station_id = 72466693067`"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/tmp/ipykernel_4377/3032895316.py:6: DtypeWarning: Columns (42,43,54,94) have mixed types. Specify dtype option on import or set low_memory=False.\n",
      "  station_data = pd.read_csv(csv_path)\n"
     ]
    },
    {
     "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>STATION</th>\n",
       "      <th>DATE</th>\n",
       "      <th>SOURCE</th>\n",
       "      <th>LATITUDE</th>\n",
       "      <th>LONGITUDE</th>\n",
       "      <th>ELEVATION</th>\n",
       "      <th>NAME</th>\n",
       "      <th>REPORT_TYPE</th>\n",
       "      <th>CALL_SIGN</th>\n",
       "      <th>QUALITY_CONTROL</th>\n",
       "      <th>...</th>\n",
       "      <th>OD1</th>\n",
       "      <th>OE1</th>\n",
       "      <th>OE2</th>\n",
       "      <th>OE3</th>\n",
       "      <th>RH1</th>\n",
       "      <th>RH2</th>\n",
       "      <th>RH3</th>\n",
       "      <th>WA1</th>\n",
       "      <th>REM</th>\n",
       "      <th>EQD</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T00:00:00</td>\n",
       "      <td>4</td>\n",
       "      <td>39.84657</td>\n",
       "      <td>-104.65623</td>\n",
       "      <td>1647.2</td>\n",
       "      <td>DENVER INTERNATIONAL AIRPORT, CO US</td>\n",
       "      <td>FM-12</td>\n",
       "      <td>99999</td>\n",
       "      <td>V020</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>SYN08072565 32866 63404 10000 21056 38317 4016...</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T00:53:00</td>\n",
       "      <td>7</td>\n",
       "      <td>39.84657</td>\n",
       "      <td>-104.65623</td>\n",
       "      <td>1647.2</td>\n",
       "      <td>DENVER INTERNATIONAL AIRPORT, CO US</td>\n",
       "      <td>FM-15</td>\n",
       "      <td>KDEN</td>\n",
       "      <td>V030</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>MET10712/31/20 17:53:03 METAR KDEN 010053Z 340...</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T01:53:00</td>\n",
       "      <td>7</td>\n",
       "      <td>39.84657</td>\n",
       "      <td>-104.65623</td>\n",
       "      <td>1647.2</td>\n",
       "      <td>DENVER INTERNATIONAL AIRPORT, CO US</td>\n",
       "      <td>FM-15</td>\n",
       "      <td>KDEN</td>\n",
       "      <td>V030</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>MET10812/31/20 18:53:03 METAR KDEN 010153Z 340...</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T02:53:00</td>\n",
       "      <td>7</td>\n",
       "      <td>39.84657</td>\n",
       "      <td>-104.65623</td>\n",
       "      <td>1647.2</td>\n",
       "      <td>DENVER INTERNATIONAL AIRPORT, CO US</td>\n",
       "      <td>FM-15</td>\n",
       "      <td>KDEN</td>\n",
       "      <td>V030</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>MET10612/31/20 19:53:03 METAR KDEN 010253Z 000...</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T03:00:00</td>\n",
       "      <td>4</td>\n",
       "      <td>39.84657</td>\n",
       "      <td>-104.65623</td>\n",
       "      <td>1647.2</td>\n",
       "      <td>DENVER INTERNATIONAL AIRPORT, CO US</td>\n",
       "      <td>FM-12</td>\n",
       "      <td>99999</td>\n",
       "      <td>V020</td>\n",
       "      <td>...</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>SYN06472565 32966 60000 10000 21078 38343 4018...</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>5 rows × 97 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       STATION                 DATE  SOURCE  LATITUDE  LONGITUDE  ELEVATION  \\\n",
       "0  72565003017  2021-01-01T00:00:00       4  39.84657 -104.65623     1647.2   \n",
       "1  72565003017  2021-01-01T00:53:00       7  39.84657 -104.65623     1647.2   \n",
       "2  72565003017  2021-01-01T01:53:00       7  39.84657 -104.65623     1647.2   \n",
       "3  72565003017  2021-01-01T02:53:00       7  39.84657 -104.65623     1647.2   \n",
       "4  72565003017  2021-01-01T03:00:00       4  39.84657 -104.65623     1647.2   \n",
       "\n",
       "                                  NAME REPORT_TYPE CALL_SIGN QUALITY_CONTROL  \\\n",
       "0  DENVER INTERNATIONAL AIRPORT, CO US       FM-12     99999            V020   \n",
       "1  DENVER INTERNATIONAL AIRPORT, CO US       FM-15     KDEN             V030   \n",
       "2  DENVER INTERNATIONAL AIRPORT, CO US       FM-15     KDEN             V030   \n",
       "3  DENVER INTERNATIONAL AIRPORT, CO US       FM-15     KDEN             V030   \n",
       "4  DENVER INTERNATIONAL AIRPORT, CO US       FM-12     99999            V020   \n",
       "\n",
       "   ...  OD1  OE1  OE2  OE3  RH1  RH2  RH3  WA1  \\\n",
       "0  ...  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN   \n",
       "1  ...  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN   \n",
       "2  ...  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN   \n",
       "3  ...  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN   \n",
       "4  ...  NaN  NaN  NaN  NaN  NaN  NaN  NaN  NaN   \n",
       "\n",
       "                                                 REM  EQD  \n",
       "0  SYN08072565 32866 63404 10000 21056 38317 4016...  NaN  \n",
       "1  MET10712/31/20 17:53:03 METAR KDEN 010053Z 340...  NaN  \n",
       "2  MET10812/31/20 18:53:03 METAR KDEN 010153Z 340...  NaN  \n",
       "3  MET10612/31/20 19:53:03 METAR KDEN 010253Z 000...  NaN  \n",
       "4  SYN06472565 32966 60000 10000 21078 38343 4018...  NaN  \n",
       "\n",
       "[5 rows x 97 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "URL_BASE = 'https://noaa-global-hourly-pds.s3.amazonaws.com'\n",
    "year = 2021\n",
    "station_id = '72565003017'\n",
    "csv_path = f'{URL_BASE}/{year}/{station_id}.csv'\n",
    "\n",
    "station_data = pd.read_csv(csv_path)\n",
    "station_data.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There are a lot of columns. I'm mostly interested in temperature and dew point, but they are in a strange format. According to the data format document, the values include the value in degrees celsius multiplied by 10 followed by a data quality flag separated by a comma. We will need to parse this out. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "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>DATE</th>\n",
       "      <th>TMP</th>\n",
       "      <th>DEW</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2021-01-01T00:00:00</td>\n",
       "      <td>+0000,1</td>\n",
       "      <td>-0056,1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2021-01-01T00:53:00</td>\n",
       "      <td>+0000,5</td>\n",
       "      <td>-0067,5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2021-01-01T01:53:00</td>\n",
       "      <td>-0006,5</td>\n",
       "      <td>-0061,5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2021-01-01T02:53:00</td>\n",
       "      <td>+0000,5</td>\n",
       "      <td>-0078,5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2021-01-01T03:00:00</td>\n",
       "      <td>+0000,1</td>\n",
       "      <td>-0078,1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                  DATE      TMP      DEW\n",
       "0  2021-01-01T00:00:00  +0000,1  -0056,1\n",
       "1  2021-01-01T00:53:00  +0000,5  -0067,5\n",
       "2  2021-01-01T01:53:00  -0006,5  -0061,5\n",
       "3  2021-01-01T02:53:00  +0000,5  -0078,5\n",
       "4  2021-01-01T03:00:00  +0000,1  -0078,1"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "station_data[['DATE', 'TMP', 'DEW']].head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "station_data[['temperature_degC', 'temperature_QC']] = station_data['TMP'].str.split(',', expand=True)\n",
    "station_data[['dew_point_degC', 'dew_point_QC']] = station_data['DEW'].str.split(',', expand=True)\n",
    "\n",
    "station_data = (station_data\n",
    "                  .astype({'temperature_degC': float, 'dew_point_degC': float})\n",
    "                  .assign(temperature_degC=lambda df_: (df_['temperature_degC'] / 10).replace(999.9, np.nan),\n",
    "                          dew_point_degC=lambda df_: (df_['dew_point_degC'] / 10).replace(999.9, np.nan),\n",
    "                  )\n",
    "            )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "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>STATION</th>\n",
       "      <th>DATE</th>\n",
       "      <th>TMP</th>\n",
       "      <th>DEW</th>\n",
       "      <th>temperature_degC</th>\n",
       "      <th>dew_point_degC</th>\n",
       "      <th>temperature_QC</th>\n",
       "      <th>dew_point_QC</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T00:00:00</td>\n",
       "      <td>+0000,1</td>\n",
       "      <td>-0056,1</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-5.6</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T00:53:00</td>\n",
       "      <td>+0000,5</td>\n",
       "      <td>-0067,5</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-6.7</td>\n",
       "      <td>5</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T01:53:00</td>\n",
       "      <td>-0006,5</td>\n",
       "      <td>-0061,5</td>\n",
       "      <td>-0.6</td>\n",
       "      <td>-6.1</td>\n",
       "      <td>5</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T02:53:00</td>\n",
       "      <td>+0000,5</td>\n",
       "      <td>-0078,5</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-7.8</td>\n",
       "      <td>5</td>\n",
       "      <td>5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>72565003017</td>\n",
       "      <td>2021-01-01T03:00:00</td>\n",
       "      <td>+0000,1</td>\n",
       "      <td>-0078,1</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-7.8</td>\n",
       "      <td>1</td>\n",
       "      <td>1</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       STATION                 DATE      TMP      DEW  temperature_degC  \\\n",
       "0  72565003017  2021-01-01T00:00:00  +0000,1  -0056,1               0.0   \n",
       "1  72565003017  2021-01-01T00:53:00  +0000,5  -0067,5               0.0   \n",
       "2  72565003017  2021-01-01T01:53:00  -0006,5  -0061,5              -0.6   \n",
       "3  72565003017  2021-01-01T02:53:00  +0000,5  -0078,5               0.0   \n",
       "4  72565003017  2021-01-01T03:00:00  +0000,1  -0078,1               0.0   \n",
       "\n",
       "   dew_point_degC temperature_QC dew_point_QC  \n",
       "0            -5.6              1            1  \n",
       "1            -6.7              5            5  \n",
       "2            -6.1              5            5  \n",
       "3            -7.8              5            5  \n",
       "4            -7.8              1            1  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "station_data[['STATION', 'DATE', 'TMP', 'DEW', 'temperature_degC', 'dew_point_degC', 'temperature_QC', 'dew_point_QC']].head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAFQCAYAAABaoA7DAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAACBEElEQVR4nO2dd7wU1fXAv2f3NXoTEAQEFAugoGJF7B17r0FTTIwmGlv0Z2yJLWpsUWM09l4Sey/YG6ioVKkqVYp0Hrxyf3/cmd2Z2ZnZ2b77uN/P58FOPzt758y55557jiilMBgMBkPLJFZqAQwGg8FQOIySNxgMhhaMUfIGg8HQgjFK3mAwGFowRskbDAZDC6aq1AI42WCDDVTfvn1LLYbBYDBUFF988cUipVRXv21lpeT79u3L2LFjSy2GwWAwVBQi8n3Qtry5a0QkLiJfichL1nJnEXlTRKZa/3fK17UMBoPBEI18+uTPBiY5li8C3lZKDQDetpYNBoPBUETyouRFpBcwEviPY/VhwIPW5weBw/NxLYPBYDBEJ1+W/C3AhUCzY113pdQ8AOv/bn4HisjpIjJWRMYuXLgwT+IYDAaDAfKg5EXkYOAnpdQX2RyvlLpbKTVMKTWsa1ffwWGDwWAwZEk+omuGA4eKyEFAHdBeRB4BFohID6XUPBHpAfyUh2sZDAaDIQNytuSVUhcrpXoppfoCxwPvKKVOBl4ARlm7jQKez/VaBoPBYMiMQs54vQ7YV0SmAvtaywbDesW9H85k2k8rSy2GYT0mr5OhlFLvAu9anxcDe+fz/AZDsWhqVqyob6Bj65qsz9HcrPjbSxO5ubaK8Vfun0fpDIbomNw1BoMPf39tMkP/+ibL6xuyPoeI/n/l2sY8SWUwZI5R8gaDh09nLObu92cAsGx1LkpecpKjvqGJtyctyOkcBoNR8gaDh+Pv/jQv55m/rD6n4695ZRK/enAsX/3wc17kMayfGCVvMBSI1euSbppPpi/O+PjvF68GYGkOvYl88f3iVexy7dt8M3tpqUUxZIhR8gaDg7WNTXk7l9Nd88ingUkCA4lZhytU6H71DU3MW7bGd9vqdY05uZxs3vtuIXOX1fOvd6fnfC5DcTFK3mBw8MsHxuTtXEollXMslrl/3n5JNDeH73fJs+PZ+dp3fLft/Y/3GPLXNzK+dhCvjp9PfUP+XoSGwmOUvMHg4KNpwW6V61+bTN+LXs7qvPEMdXxDUzPvTNaTxMPtePjvl7MBuPLFCXz5w8889tkPPDnmBwDm5Tgu4Me6pjRvHUNZUVZFQwyGcsNhjHNnDq6KTCNtPp+5JPG5WaVT85r7P5rF/R/NSiwft32fjK4ZhgR8NpQ/xpI3rPfsdM3bnPPEVwW9xqKVa0O31zc00feil3nY8t073wnKR8l/8f3P9L3o5ZTZtM7jbnrzu8jyKaXoe9HL3PWe/4vs0ucnOK5h1HwlYZS8Yb1n/vJ6nhs3t6DX6NCqOmXdTyvqWbm2kfqGJibNWw7Apc+Np6lZIQ57eeK8FSmTsl78Wsv7/nfB6blve3tqyjqlFLMWrUpZX9+gXTA3vZH+xWBL1tjUzI9LVqfd31BajJI3GELwi2zxs6yzYYer32bkbR/w+0e/5Ig7P06sv/XtqS6L/La3pzLsqrdcx8bsQdkMZXnmi9nsceO7fDx9kWv9sjX6JRLF327Lds0rkxlx/Wh+Wp5/v78hfxglbzBEYMmqdYnPzRH1apTdvl+8OjHAavP5zMXc8PoU17p1jW7lmwiv9Fwknc6fMFf3GE685zPWNTbzm4fGMnHucm4fra3+Tq2TPY5lqxs45d7PWOBR4nYv48NpuhfxcxnE8RuCMQOvhvWapjQa21aaj32WjHNvVop4jsOPYb2BmAhffL8kcDskXyBNGVryMUcX4ZvZS3lz4gLmLVvD+DnLU7Y/8+VsPpi6KMVPb/dubGWfLo7fUFqMJW9Yb7lj9DS2u+rNxPIFT3+dso+tvqriyUcl6MWwvL6Bvhe9zKOfpU588g5WXvDMN4FyfZxmdmxzs+LeD2cCcN2rk0P3TTnW8VKIW92BxqbkOuc3W7JKDxYHvUfsr3TALR9w4K0fZCSHoXgYJW9Yb7nh9SmulAFPfzE7ZR/b4naq6IYAv/UCKybdGcYYxDM+14pKpta7k8c+/yFl3eT5KxKfnT2M2qo4AK1q4q79/S5vDxwbyg+j5A2GEPyM9iAVa1u2mQ6GZkq2p1+8cq3Lt/+eT2SO36njASGTJpSyMjA+ecN6RUNTM4Muez3yrE0/37kKPFQrvRkLV7FgeX2KMr77/en8tHwtfzl4YAYS+8iUJx+4HTbpZOnqBpRSLgW+ZPU61z4KnSPfWO+VgbHkDesV3y9eldG0fOX5X3/2V7LO9DS/eWgss39OxpALOuTwPx/OZM263HK/ZGvJxz35c4IGf58bN8d1nTYed02zUnw4NTg+38nUBStMiGWJMUresF4RNfwxuX/qAUHnEFfkyjJOvd8/2dkdo6dlJoSHbJW8eCKCNuna1ne/yfO0j76mKuZ7veZm5YrCCWPfm99np2vfzlBSQz4xSt6wXpEuZNLLa+Pnp6zLdTLUl1kWAalvaOKmN7/LOh2y94XVo2Od735eX7v320a9hx9PW2RdN5p8hsKQs5IXkToR+VxEvhaRCSJypbW+s4i8KSJTrf875S6uwZAbmSr5W95KTQ0QaMlHPGe274j/fDCD296eygMfz8rqeO9lGwO+iNcd5X2pBN3DRo8b7MT/fJaZgIaCkA9Lfi2wl1JqCDAUOEBEdgIuAt5WSg0A3raWDYaSko/Il2CffLCaf+HrZG6cT2ZkXiUKYJXly/cbMI2C97sH9kg8qx/51B12GRTCec6T47KSy1BYclbySmOnwqu2/hRwGPCgtf5B4PBcr2Uw5IozhW9UFnszSKaZHFQobCWdRf0RwM+3Hr7/uB/93UpNzf6vuZe+mZf4nM8KW4bcyItPXkTiIjIO+Al4Uyn1GdBdKTUPwPq/W8Cxp4vIWBEZu3BhtBF7gyFbrnp5UsbHHHHnxy4FWSofc3OOF/aq5iCL3F77+oQFvtubmlVal9M1WdxnQ2HIi5JXSjUppYYCvYAdRGRwBsferZQappQa1rVr13yIYyhjvvh+CU94Zl3WNzRx7SuTXIWvy4kfPOl0g9w1z341p6By2AOi2ar61GRm/me6+/0ZoYPLT4z5kTcn+r8AbGYuNimIy4W8ToZSSi0VkXeBA4AFItJDKTVPRHqgrXzDes5R//oEgON3SFYteviT7/n3+zOorYpx7n6bl0q0yAQZ1JkU6cgG20uzem12L8Pla9zZIr+ZvSxwX2eqAy9Rinln61Iy5J98RNd0FZGO1udWwD7AZOAFYJS12yjg+VyvZWiZ2FEeaxvLt3ao03rPVz75jLEU54OfpCZAi8K+N7/vWg4rZ5hrwrGocfSGwpMPd00PYLSIfAOMQfvkXwKuA/YVkanAvtaywZCCneAxk8iXr374uWRViWYsdFdWmrdsDWNnZT6gmyneyUyZMHVBsGVeCibPX85/PpiRUvHKkH9ydtcopb4BtvFZvxjYO9fzG1o+ySpH0Y+xKynNum5kIUQK5Rf3fe667p43vpt1WGMm5GIce634UnPALbqn8N53C3n4VzuWWJqWjZnxaigJTpeHZFnKLoib3pjCwMtey8u5olAMBV9JDL/unZRqV0F4C5FHYcmqdfS96GXemJA6G9mQilHyhpLgtNqDStlly23vTGO1TxKwlVkOWAI0NZXID+8gJV6/TJmzdE3KOmf4pzOGPpsX++T5OvulXTjFEI5R8oaS4LTkbXdNpikHMuWMR77I+th/FDhyJgpPjc2+0EipcRYrOfPRrxKfM0gImiCXsYn1EaPkDSk8P25OzpkS0+Gy5GPZu2u++D56sq8Ppi7K+PxBjJ5iIoIzwTnT+K1JyRj7RVn0TuyxidL3rSoDo+QNKZz9xDhueH1KQa/R7LLk7XWZn+fEez4N3FbIUMfTAtIIGzKjbW3msR8JO95o+UgYJW9w8cin2cVg+3HALe/T96KXedAna6JT/9rd72yUcljESanC2Q3RyeVFnK5C1phZS+h70ctZDe62JIySN7i49PnxeTuXPWvy5rdS/dl+lnw2z3uYf9bo+OJjD4p6CXoZZ/MbJdI7pDnYTjORbdbPloJR8gYXdVXx9DtlSIPPTFbn85kYeM2zJV/ogtqGVOz49zBynQwb9Xg7oieoEPn6glHyhgQr6htY01C4FLHzliVD61yWvD3wmoVT3vn4Tpm/gr4XvZxYLoaOz7bK0/rKz6vWuX6X1euaOOaujwP3HztrCb95aGxKQRJI3wuwdfv/PfttFpK2HPKaoMxQ2fiVussn/35vRuKzcjyzudhZzlJ1h9z+oWtbOp9tPrj8+QkFv0ZL4ukvfkxZN2bWzyilUsoOAvzlufFMnr+CBSvWslHHVkCyvZQsh1CFYSx5Q4JCJwhzPsP5cqU41UJt3N2ci6EDcplgtT6yZJV/rhq/TtysRasS4zruGdL6/y9/WMoKk/smLUbJr8eMnbXENTW8voCuGoClq5MPZBT929jUzB2jp4XmmV/tkLm2uvhKfj1392ZMXbW/yvF76X89e2nic9Bv+dOKypgFXEqMkl+POfquTzj94eQs0EKlh7WfT+cD7vdQe9c8N24uN7w+hZtDZps6Z8nWeC35IrhrTErdzAgad/FrD37uG2uLz6fw/dZn145R8oYEf31pYkHO6/d8RRlktQeBV3ny0DznqcA06r7PefmbecxdVu++huMSQaF9uWKKY0TDvk01Vf4qJ50Ods2rcNzzoJesUorHHakU1mMdb5S8oXg4C0c7n7kJc7UC/naOu1KR/fjWe5T8OU+OS3yuqYrx3ncLOfOxL1Ou57Tebngt/zN43/tuISvrjU8+H6Qbo1Eolq1p4ENPagoR+Gb2UpauXuda7x1f8p590cq1TJgbXBkr5fpK8cHUhRXZIzBK3lA0nO4T50N930c6m6B3ZqJtpIUNCPfoUBdyvSQNBUh+Nuq+z1N6D4bs8Pt5nApVKbjqpYmcfO9nLPOUMTz09o849t+fuNZ5k915lfMBt7zPyNvc0VhhPP75j5xy7+c8P25u5GPKBaPkDUXD+ZxF0bm2j/3lb+dlZUE5wzQr0QJriQT9DOkteXh/6kIA1joG2+0Zz98tcBsI69Ja8uuIyvPj5vDsVzoDqF8a5XLHKHkDUFglaFvwzitE8ck7B95m/5z5w+XsOaw1hT1KyuJVWqkG9XxWrw2P7Gr2WPU2QWOz3t5f0EvEG7nV3KxSegpnPzGOMbP0pLefV0V/OZQLRskbgOIMTDkftEzj5IP2DzuN8z3SsXV1Rtcz5BeltAvFORjqZKdr3057fCKRXZTr4XXX+O93yr2fu5ZvfGMKQ658g2Wr/ePv/1OBhUqMkjcAMHPxqvQ75UoEd02mPYofQop5GxdN+aBQTP0pl2LiKplH3vGzBv3EUX96bz2Cl76ZB8DSNZVnsQdhlLwBgL3/8V7Br+Fy1wQ8hUHKP5tqQOWg4of07lhqEcqC5uZoycsCj3cqdscv60xqN+7HpY59vMeXQ2soDTkreRHpLSKjRWSSiEwQkbOt9Z1F5E0RmWr93yl3cddfHvn0e/b6x7vsdeO7vGxZG4Vm0cq19L3oZT6alltFJfv5ckdLRJ8Uky1l8WCXgwxlQJR0v8Ove8dlWZ/9xLjEZ+2u0Zz1mLN8YPL+OsMove0r6s9gn68hpC5hpfUQ82HJNwLnKaW2BHYCzhSRgcBFwNtKqQHA29ayIUv+8tx4ZixcxYxFq7jq5cJMWvLypfXA/evd6SzMcPq430Oy0FHqLegZcj60Ttt9YTZFrJXvx6JSWeqgtMxZuibQZ6/wT2AWdawm6u9gR89EHeupBHJW8kqpeUqpL63PK4BJwEbAYcCD1m4PAofnei2DJszKyCeNVmv+cNoi9rzx3YyOPfepr13Ly9Y08NG0pDWX6UDqUf/6mJ9WZBaTXg7PYoUZfSUnMMQyglEQdquj9OpmLkqOS4VXHKusHzWvPnkR6QtsA3wGdFdKzQP9IgC6BRxzuoiMFZGxCxcuzKc4hhxxPhiZZlt88Wv3pJHlnrC0HwMGTMMKh3w7exljZi0J3O7FKf+bExeE7OnP6bv1z/gYL1v16hB533Z1qZm/Lz9kYM4ytASUY+DVSVOzf3dt7PfuduKaoxFgis/+2d0m353yE1MXpA4WfzR9MfMraBJc3pS8iLQF/guco5SKnChEKXW3UmqYUmpY165d8yWOIQ9kM9hpE3ckdVHAOk/vw06M5rWKwiyuXz04lmPu+iRwu5fXrfz42eatsfOX50KXNjWR9+3VqXXKOjNwq1HK37puDFDYf3rS3ZN0vgAe/cy/jrGz6dU3NHPq/WPY9+b3U/Ybdd/nHHhr6vpyJS9KXkSq0Qr+UaXU/6zVC0Skh7W9B/BTPq61PjJ+TvQcG/mkKp6Dknc+kSoskZR7OZvqUEFc8aIeu8gmv8ypu/SlOp7745FJD2L7vqmxCW1rg+v67NC3c1YyVSJBL3+nJb88JLe8fXxzs+LSgEIvUaK/bH4OiKNfbAUrvD6hsAV4MiEf0TUC3AtMUkrd5Nj0AjDK+jwKeD7Xa5UTPywOjs/ONy9+k/98GVH8+jm5HiV0MYH3YSrEoJY3j0kULjpwi/wLkob/O2jLlHVhsnt7Ry2BoPTQSsGPS1JnPTvbTzwkJai9V9g9cxUmybIXO8UqcnL/R+UzaSoflvxw4BRgLxEZZ/0dBFwH7CsiU4F9reUWwesT5rPbDaOz8vNmRUq7zz2/7WUFLlvnlTDIkvfqsGwUcjqyOWdddZzxGWQpzJWaqhh11XG6tat1re8c4u7JJtVDpRL0W7gGXkMjYvTGJSFpCZzRO9nWIrDPUU4ROPmIrvlQKSVKqa2VUkOtv1eUUouVUnsrpQZY/0cfMStzfmv5k72pcYtH7i3o3SlRvGf5a6lB0QpeS74QkQtRz3hq/DV+GX81sVxMJfrFX/YBYPT5eyTWbbdxJ7q3D86yuSibsNJyJ+DHmrXIf0Z21Be43az8XDrj5yzjdw9/4XIVZmMYvDFhfjK82ccN+YfHv2JsBoED+cLMeM2FIoVSlZFREBm3S15Fn4xSwvC0K6of4rLqhxPLuw3YIGWfQhSC2qFvZ9rV6dw6bWqrEr727u21Vf/HvQfk/6IVRpCXJWjg1YttmXdpU5uy7coXJ/DahPlMX5jMZJnNRLrTH/4iURvB2xNYsnodL349l4v+923G580Vo+RzoHTqqPLKEQWnMSi8Tz6M8bW/5PHqq3y3+fl4+3ROjYAB6CU/sY1M9d3WKU1ytF+N6OdaPm14XyAZG77n5u6os1N36eta/uiivQDYY3P/6LRZdScyq+5Eakm6Ks6vepJZdSeGyuXl3H03y2j/fBLYfjK05P3cMMvXpA7MRxnueHKM/8Qtfbz7OnZLCnMXFYr1TsmvWdfE3e9Pz4vvNx9G57gfl/LO5FTfvlKKu96bHhhPXgxe+Dr7AV9x1dcMfiF6fwbnQ5vL9aPSVurZOR59BvHv99jEd/2HtefwbO3l7BH7KmXba+fsFnrOuKd7YPt17V7NNn06sZ1M4fCYLnLhLaG3Yfs67j9te247YRvX+l03dfdE2pNsS2dV6TiIQTIrsS7dy6g2oHRfIYjqtos6GSqZWiN1W9r4+wCufDF6u/nQSg2yZNU61jaGp1XON+udkr/pzSlc88rklMk6peLwOz7ilw+MTVk/Z+karnt1Mpc9P75kM+xe+Tb7MDDvgxM8wzU4Tv697/IzOS7K7etB+twqACvShGM+UHNDyrWdfvV+G7RJOSbmeQp37NeZ9nVVnOF4oTxZezW31NwJwPHb96ZdbRXn7rsZG3dpTUxgz8270b7OraRHbt2DDdomB26V48W7VGk5rqn+T2Jd2CAvhEev5Av7p0oJrQ34ETN114S1BefAaxR3zaFDeka6Nrjz8LzybXFyT9msd0refkjXNBT3bZopV700CYDpC4uQArgIBFcEci8XIrrGyz5bpk6+rpHkgNz5VU8GHltbHc/oWt6XnXNg9bfxF/lT1dMp96ZTmxq+uWJ/tu3TCT7+J/xnH6rQ7bYrP9O/a1u+vXJ//rj3AN67YM+UnC77xcZwe/WtVMWEsX/ZN7HeeZmOotvVZjI7sS6dwizGb2PjvVIUd1+YMWSLHtXXHuW79g5w3UF4r6LYdYGDZ1oY0pJtmJVNUGECgNcKPJmi0DaZ9/xBD2CqJZ9/WWYuWulZk/rtO5Ocvm67MvzuUV0eXRYXVz8OwFvq73rF7LGwaCr0HArdrJj5N/7iOua0qteBk0PPe3fNzQA0bLORa73f94lLakRJLesYKN/zlXIP+BZzvMSrjCfN889FH9mStydDhSh55/2JEgAQZRBeKcVnMz0RNYUYvQ9hvbPky4mT7/0s7T4ihQniKfTz6o45Dr5eULe8MY8TfbwzHDffsG3KPhtLtDkP2/TpmLJOyE3WxC34z97w3O/gzp0C95UMfjnvjF0/OWuqky6eY4f1ZrPubbmq6j6erb2cXuJ2lxUzdbP3Ut7iHjYZD7xG/AqRylNGMJUe/ewHjr/7U39hisR6q+S99/mpMT/ynU8yIpt7P5yZ9/QCUePsKzKE0rPspyDGz1nG/76a47tfIUMpN2ibGka3ScwzRrNsTso+AJt2a8fw2LccEUsWwKjG7fp75eBG3j9wMXdU30KH5qUwfTS7x5K5VDbv3o4NHX76TvM/jix7KzKMjW9Ixvp3sFw0/zhmSHK7JFXAH/balFfP3o1jqnRelg64XYVKKaZefWBm18+Q2T+v5oGPZkbuJftNhvJ7Tu1toZa8o9Hm6pqyLxMU319M1jslH9RTuvC/37CfTzIim7+9NJEzHv3Cta4YL2Qp0HWKHYTp9x0O/ueH/O0ld4SC/XAV8t5+7aggZNMdj6X45EkctHUP3yRlj9Zcy801/0osV+P2sQ586xf0Gf0HRsY/565298LDh/Ngzd8T23t3buUa5Bz2/qmRZW8kfEzgCI+Lho//mfh4ftXTAOy5hWNMwqHkRcQ1uBrzWP5Nzam9g3wzZtbPXPHixLQD3EmZUhvK7x/9MmVdlIFXJ1F6LUX2umTNeqfkc+HHJWu411HItxgWtl+hhCg0NSsOu+MjRk8ufl64tyYuYIUjNXFTs+LaVydHOjaKxZUvnO6LA+Pugs4s/I5u7er46KK9cAWVjL0/5TwXhAzUdprzrnvF7C+4YvavqcOTqna5pydxRQf956EhzTDazUd6cu5881Ti4wHxMbzV5QZ3qKQKdjVV0cTtJ27D2dZkLLt3VQzlFvUSja4QSv3Zr+7vjEWr2OOG0b7FZ+xnzHnNq16elDcZvRS7Z26UfIZ4Lc9cyDYiLcpDtmxNA1//uJQ/PTUuu4v4EDWU86zHUy2p9yOGQ2ZqyW8ic+gd0Z9uE7NufJVDybcTTwqDBqub/f3HtMax7aVzUs63dWxG9Iu/fjG9GmbRbek49/qvH490+Ha7Hgg/joHVAdPjX7/YvdzgVnibrvoKcSp2r5Kfk+ytdpefOXBwj4R1b//+xejBRvWWPPtVMjooTK7T7h/DrMWrfUtnTpqnZ6k6Dap8pbQoB1erUfI5kGtjj2KlC7lF8eTT6Ir6fXO5L4mUsBFP8nbtBXxQ+6eMrlFtBabHSRNGu3oJ3H8gt1X9M3S375p7Rb/4j3qwfd/Gd6Mf42CHvp3g3n3goUNTN675Gcbel/4kXz7IFFvmPjvTuibO1nZxk3v2Sux2V80txGPCCCu9w26bFa/eQ1SDYsws/wHZUlMOyt2mxSn5pavXceWLE9JWMso1/DEXmpsVN70xJdrgjsD9H80quExR8FO8v3pgDDe/+Z1rXS531r5Gvn4dvyiduJUn3ztgmsIPOipir/g43822z3pOVRol33Mbpjf3cK3aRDwW5awPw89hs9YKDpjvkwNlnc8gX8zHh7/0Rzbf6SD9ecPBTPzjprzQ/T8w7xvfS27TpxOzrhvJ9lnkr5913ciMj4HsBt6jHBL28sh00mGUHnU5uO1bnJJ/+JPvuf+jWTz75eyAPTK/7avX5XfywqczFnPbO9Mi7euXVyMKxcrm+Pbkn7j17anpd4yI/d4LcRVnxOT5qRFTZ+yuZ5KmteR/nhm6+aCYtsq7VAfPdwAgVsUarIieviMA+KB5K/c+098JP4fN+9cHbxvn4/LZYHP/fVdZ7jPVDEumw4RnYUX5FLooVEz+sjXBv1WmJS6jhFCWg0Xf4pS8PeC3ep3/A7zWmuka9mN7GXjZ677rs+0NRJ3AkQv2FbIduPUjqgsll15Swicf4RwdWoXnWgFoVZNqydozFavSxbe//n+uRe8Yyu012o1zakPwwCu9d4TZYxgcm6WXO/UF4A9Vz4VfO4jFIcZBo0/dUT9LHrRSB1i5EB47Vn9e5pNwqyn6c5JPsqkQFuWI58YFpzO58Y3vArdlLkw5qHdNi1PyNuN8wuQAJlqDLLVVmU1P96UYIZRZ6mi7jeWzuxh1MCo/Pvn0+zpf1F1Z6tq2EQsZItN4zhOHf9yw3onPTkt+arMn9HDTfV2Lm8sPXHHooBQZ9ouNCRdy7jj38vTR4fvnm+Y0vZWVjkHrH32+S3129XGd3HzckPQ7eZj2k3eWcnqKOVkLsn82py5YmaggVQxarJJ/dbx/19Puvuezlmgh8Wu3qzPsVuaLvf/xXqT9cnnY7PS6UdxNTsV+bfU9rm0f1Z3N87WX8U+PW2zjDbQVv1P/zmzYNvmib/a+Dqe96Vp8vfYidu2Z+lTbKQQA6LlNynaakiF7h/deBcuD3Ij5wOee+fm9qlIngwGw1KfA9dpUJf/b3fpnJJXf5LN0nPZAmpdnBfPwp9+z/y3FKwTeYpV8Ou56b3rO51DoWqnDr3uHp8b+mLtQEWkfwU1RqQRZ8pdWPcwd1be41rWSpAIdERtPH1nArLoTGSJJxe7NmR4HuKIDT8w9gOcbz8hItrhKYxV3CB+AveWwzJRjKFU+FaO8L8Y1S/2VfI+hzoOSH1ulFhL3eyFcvN8miRz1Uci2XmrGFNlu+/f7M9jsL6+m3zGA/3yQQehtDqy3Sn6xI3l/LoOUq9c1MWfpGq6OMHkicb2MrpC6dyZRBKWYlZfLs9aUiMV2n+VXVa8yMv45QjPd0THiVQ53SxMx9o3pGO+j4h/gRxWN7NrWP11BlJwwrevCU/Gy16Xh2/PpTui/h98F3IsiASPYjv2cMg083GdXn+OX+9/DIIqQoRgofsTcwhVrWdfYzIyFqa6lKIO4USZc5YP1Vsk7yfbZU45qGPlMqOW+hs+6CI35HWum66KVxa9Ek4suU2lCKP9U9Qyf1Z3Fhix2KebWspZLqx8BYGDMx+0A/KP6Lga95BNfTrSxi3a1acZxqlNTILhozqObLVL4UYCSd/1AToXv01Px8+nfuXOEazukqJT5/1myl48bs5xShOdFyYvIfSLyk4iMd6zrLCJvishU63+fvmB+eeCjmdz9fnAXKDDdbQ7XbE5YzNEbciY9B789Fyxfy0dWpZkgvp69NHDb2sYm5i5LRmJ8+UP5TCix35VOv34HkpbSMXHty9xAlgXmcdlY/MdjDosHJwKLYsnX3bpF+A6xNJm7PdE6GdHe4wpqrIefZ8F/fw2N1ovc265E/JW0U/GvdQwA+r2EvIr/51nQGG0A/pFf7ZgQw8uR3hw7eaCUAS3TfkodSL3mleJY6unIlyX/AHCAZ91FwNtKqQHA29ZyQbkiTTmuoEaQi7smEaqY9RnSnD9AtpP+kz5NcdD5Pp7mroJ09L+iZ0EsdJWqJmvk1XmZ31S9nPjcQ7Srpom4y13jJvNfo0OrPJRWkDSW/tzUdA+R6eLx58dr4MVz4NunYaY9iOf5bZQKsOQd6+Z+lfzsp+SbPce/eHYkcX85vB+7bNIFgJiPlk9XhSobShlK4az8ZBNmcBaTvCh5pdT7gDeZxmHAg9bnB4HD83GtXAhqBNkG2ijlsDjzqOVveD2ZzCtMtKdDBnuDdPEp937OX54b71qXyfcfdtVbrKgvXOx0Q1NqWgO/XO+tqQ9U8t5JTvYAYbMK/pG6ts2D0gmKSc8HPTxhiNPeghmj3det9pYWjKDknfha/Z513kRqFk6/e7d2tVx2yMBEjqCOs99hVt2JibEUcJdEbAlMmJt5qGmUeR75oJA++e5KqXkA1v+pNdcAETldRMaKyNiFC/NT0zMIPyt01qJVWYf8KbKLR093tTtGJyN/wkQLixAK+k4fTlvEnKXZJ19avGodC5b7TLrJgJFb9wjc1uhjye8QS81gWS1NgUremwfdJiYhNzMfPRRH2l6OeSD38znZzNtRdmC7iTp7rH2loJuPiylobMBvMpVX8ffc1vfQMHdlt+8eA9yJ3H65a7/A/bOljOYfRcK+ZfOX1ed9Vr2Tkg+8KqXuVkoNU0oN69o1vwmQvErd2wYmzl3OHje+m1O3SmXhk8/m/P7bwo6Lfo22tZm5KsLOHcW/v9VGHQK32Za8fY3uLKG7LE3ZL0ZzsCXvo8z3j33us6eTPCv5QUfkfj4nte3TX9drdatm6GKV8Gvv8IEHKfm3/5a6znvOb57wPfTgkBe37cZy5qePxyT0ZZ8NxZ4MFYUoz+hO177Ncf/+NHjHHCmkkl8gIj0ArP+Lntjc64bw3vApC3QXa8ysZDdybgZW7oLl9Y70AdlImJ6wRhLWqKsyKO6Q6cPx99emBG6bESGqIOxWNdpK3rqzA2L+4XrVNEUu2Qfw75pbIu+blgNv8F8vBXycQl1B1u/ndcOsXZ7c5pTNqeTbODrYTT5Vp9LNmLW44eikO+mnFe7zKEv2uCeNxM3HDuWuk7eLdP4oBLXifL9MCkHUKnHZUEgl/wIwyvo8Cni+gNfyJdWSdy/bY0oNjvDH20dHSxwG0LF1dTK6JiPBMtk5mLBMfbtZ6WGjkKmSf2uSVq5+vYx05zp2WC/qG4LD/2x3jf2CPjz+ke9+tazjwupkzpj/Nu0aet20ZHIPary+b4tCKvmwQd1WVnZIr0Ke8R7UW8ojSMmnC8VMNwHMoiakwLkkLHn3Pa6pirHfwO7sumn0thpKwG9o56sqBaXMdmuTrxDKx4FPgM1FZLaI/Aq4DthXRKYC+1rLRcV7e71t4Lyndd3NtY3NgfuEnl85fPIZmPKZTNkOU5rewAcnZ/iUQAsi+3kCfuvCT3bCDn3o3Tk4nvyy5yfQ96KX+X6x7hH8pDr67nd3zc30kmQY6XIVoHijkibjJEMcszuDfuu4NZC2Uf6s0wSRLHmPMvv6cfjgH/qzU8k7k46lU/JhjSwq1rW95QRBF3C58ICATJkZEhRA8NYkfydCIcI4vYQ9DsvWNPD5zIDiL3kkX9E1JyileiilqpVSvZRS9yqlFiul9lZKDbD+L/y38eBVkM7Fv4zcMvG5fZ1zlDu6xmt2RNcUKoQyLPIlX+GMWSv5LM4VE/GtEzpq541dy3aCuc+boymAKHHuORFzyBxksVfVwqgX4aRnsrvGr9+GbqlJ0EKvCcmb7rXkf/gk+dm57dU/Jz+vSfNYprHkvxr6Ep8elWbCncdd028D9ws5X2mkMg0IiBVhKm66nu3rEwqf3rnkA6+FxHt/nV2nAwZvmPjs/CEys+SVw5LPSsS0hFryeXo4su1S+smWTiYR/669q7g0yQRyZ1c9G0kWXyXfoY9O9ZsPFjjnYAicqItis80psN9V0G93vdxvN2ideXENAHoNgw4B1mVoXpwAn7wTZxrhiC4YIK1PvtPkx9hwVfAYDYDE3O4a77OSrwHTdnWZhSQWI91Cuu9WjLHi9UvJK//P9mAfwBNjfuTMx6K5OpodSr5QtnxYI8imeo4fucwTSD2Xe2X39u4EVzER+k9/iMurHnStb7N6Hi/XXEwXtA/Z/m7bxKwxksFH+wthWbijqt5M3bbDr+FXb6T7GtFw+rElBpvtB1csg8Nuh13+AKNeyM91/FIjjPxHcOZISP4Q+aq04jp3hBfCu9fqZGgBVFtpvWPSbP3vflZs11yuZGpoxYtiyYdvL4bPvkUr+RR3jePzfEfXzqss/Yr9+p8fZlkNtBSWvJ0gqZDXCMOve1xX7fYdb9SxFafGX6MrydDKAV9ezWlV7kIsPac8wKDY9xwe12XwUgqr7PR7fyHCFFs+B0Kd58rnj+3134+8GXa7wL0u3e9j34OIkTAZ4Tzn2pAc75OSL7mUF7s1VmG7a7x3787RwfM9MmFMhv5tv5m4+earH5aGbjeWfI6k+uSTy8fclfRXZptbvlkpfnGfjr8unE8+XLbnx2WWEdCPbBvaiOtHp6yr8lhHvddM5orqhxhTd6a+luNe9xLH5Le4jtW3666m/CY1rTMXMK9KXvw/58pGw9zLbbrAXn/x33f7XwecpEiW/KqQyYqvXJj4+OcD3BOwYlX6t7XnNRTKXfP25J/4YfHqyPsXQ8lPmhc+E7YYsf0tWsnv+vfRjLrv80TaTzv0z0u25fi8Cq0QpBNt3rL0g02ZZsg8f7/NMtrfife56emZHtHUnIzs2FK+5zfxlwBFu0XjABgW0/7dD6Z6ErBlpbAtYeKZF61Iwfkw5vry2Po4x7kitKHEwE/AdX+w8hhl4muPSnMzjL5WJ0ILS8DmSFp25Lbu8YPqpbOApCX/yK/d4yT5VHNLVkfPuloMd006mopQvKhFK/llaxp477uFiaLef3rya9/9sn2bOn+fXI2ChgKlKgZ4+dto7ieA7TbuRKuaPCTrsvhD/H+uZdWQnChzT81NXFL9GIPke9ov0D2ifeJfAT4FuNMl//LDVop+k3wyxTnlvyH7tBAaCfichs6b+K9/7c8w6aVwd0q2qCZ47zqdCG1N+GzmITLNN0xSZr4LQA9ZzBkjetOtnTtvjbMkYzEphiWfDmPJ54mVa8MtnGzfpvWOSRa5Vr858Z7CTWteGzL5yMkem3flv2fsQjyLrzJrkT02kTy4U+tqWrdxT8dvbkxNbvZybYQUvLEY9No+db2nFqsLW8lvdWz686ejyWEhOjM3ZsOGWyU/2/drv6tDDrDaZ/eB+n+/4iRPngQf/zN7mfwqTYHbJ5+mqPfztZfxf1WPBm7/XdVL/Hns7inrN+wQnqzMGc9+0o59Qvf9eHp4Cm4nZWDIG0s+ExavzN5ay/ZGv+QYoM3VKBgz62e+nb2M+RHcL06iGAK2tTBvWTQLNJs8PD8sWW3JkxQoJgK7npNY3j32NRtMfiTjc1tSwXGPwh4Xe1aHyGpvOzQH5WfjVHCxHLMHOgeREy6Y0AkR+v9+u8Hp7wYPQtsuk3MnwzYnp5fDVQYw6D465KpKn6nz11VpyuGp5pSXhf38xWliz9hXeO9FrWMw32l97xEbRxXuPDyZTC4qD3dN4a/RYpT80Y6B1EzJR5cpH83lkNs/ZMT17+ThTG7sd9jw66KdO5u2X2WZ/84Xpoi43kIP1vydjcfdGPmcroFZgHbdYfc/u9eF+cftbdV1Oo497lBS1RkO5C5zpHXOJqXwflc5jveROeqgac9t0hcnabchbLJX+nM1OgyjoJel89lYHD3lRyjv/d21aLeZ38Vf5P6aG9gr5u4pOdujLebOsQk8UHM9Z1e53YHOcOh0lEPFKuOuyYCZi7KPtc1HlylfDaYhg0YaxIeeQcuZi1Zy6O0fpp+oZP8f4buMneW2mGqsWazOQWz9cGb+ff4Q/x99ZAHnVT2Vfuf53wZvcyrOUS/ApY6XRkDK3EikU7J+7PIHHVd/hScRlW3Vhj7snm2hLxnRmjDK4LAdk9+5P4FmitNNtSz3SC59ntmuRVvR9RE9SN9VkvfoikMG+p7imC30YHpfTxWwD9NUTHOSQQ6/gmHcNXki3YSDYtzoYnLyve6qUfd8MJNvZqfPcmcr9ygDUt6ek5310m3Jk1V85nnVz/Bg9XUc4ZeczCtbWFHp70MqXh349+BtTjbYLNUqrusQ7dgDrGscF+Kisr/PuhAjJaWsX8hja78AogxU24nWlAq25Fc4Bu3fujz9Ob34xe7bOX7WLIXGtWw573mGb7DaUWUt+X1nLV7tEs1bbtNvpnMt62iLdh/uvllw+vJ4GVjy+ZrQGMZ6oeTTkY8bXar2UogZc9l8F/uBaXJZ8kK2AXL9Yt5w1yzOs3ZF8LYNB/uv77OLe/mwO+CUZ6GNQ1n0Tx089GWn32nLfctDUrfZvQG77bXuEu2cEP4DZVIs3HZZqebgF8esD5OfO2QRBfO1X/55S/6/bwxXdWPrL//Coyt/jcJfcTuDGhLNy5LXT8mPrj2X8XW/pnObGv5+1NaBopWDu6YpDz33dBglD/y4JHVAcnyG+Z3TtZf5y+r5bMbi8J3S0Kl16oCfrSPqG5p4bXz0UMkwsvHJ2xaWreQ3lvls2ziutOV6spkc1L6n//oR5yc/ZxPO6eX8qXDelGRse6h7JYt7GEWB9bXSMytl5Z53MOol/b/zfmx5cOZyTPcZBwroYe07cMOUdc1Kub6KHSlmJxfzpi8G6GnVAV6yal1oWy6LgVfH81Eoj8J6oeSz0TMH//PD9Ds5SBdCOfK2Dzju7tzCJP96WKr1aUezXPfqZH73yJc5v0gAduyXgVVp0eRR8vdV38AdTX/N3yxMZ3GLHX4b7RhnqKJN6w1giwBltUHIJDCnHzwftVxbd9YDpPb9ybdVuSqCb7qjFY7o9xu1sXK8O8cf0oRQ+uJXUnDxVF/3VG21bZ0n8T67y9boMYImlbovkHDT2IRZ6+Wg5J0zux/4eFZBrrFeKPli8OPP4dOpF6/SjTOX9MCHDOnJQVu5rZ0Gq5HM/ln3RpatST6IM2pP4qz4s5HPbzf5vhu04fbqW3mxJkL8uoVSilvfmspfX9LZGjeJ5adXkaC2bfLzQddHO8ZPyV84HY4PiOVe9F2wss3EnZIJ9ssrLHNlNm0mikK2Y+N9X8TWfXAOvNqff/UmnPzf8HNf0UH/TX7Jf7uPW8nvznvdkXYB8F6dtavJO/mqDe6XSpQIW4CDY58wq+5E2hH8HF9zhE97ypBfxl9lVt2JiCW3M1Ah19rJQRglnyeiPoeZplCoZR0bkHQdVXnC77wNY4n1MmnPSmKiOL/66YyuZ3Nw/DO2is2KvH9TM9z81nc+G6JPMw+kNuJAp5eBh+V+bRtnzdY2vjXps2PEeXDYnTDoyPydE9y5fo74d/LzUfcmP9tx714lv+UhSQ3Y6KPkY3GyChp25unxGZCNNafOdWlWySvt0K8zl4/cnPv2amT7jTuCjxStxH2OsCAC58DrmVXPaRHFvwd0w9Fbc/z2uc/MvbhKFzWvspS800WTr/oQXoySLzKZxsU+XnMVY+vOSCxXeaaj/u/LOa6UCBf9T4cUflN3eg5SZk6gP/G5M/zXZ8LOZ2Z3XLY5Zlwx9Nb9diqLCJOCIlNVA9uclMZdk8XD71SizlnBWzlSNtuTurxKvro1/pa81TuIVYXf2x8+819//GPJz++mFoprPUkXW3Fa50op1lmDk7ts0oWa/53GXh//Ann1InsP1zk6kxxs37Rb21CfvPMFYCfGawxQiccM652XIiPeCCK3ks/59L4YJV9kMvkhq2hk25h7Asqy1clu+LVV93Ba/FWe/TI/8csp+WIywLZCalnH7dW3pj9g85GZnD07obIl3cBqrjNeo3DBDNh2lP6cSaOxZY+SrMxW1N6UBhJLvnTGO6pc2ZOnYtXhSn6JT+rgPru4j/n836n7WNSSfLE0NSsunHQ0z9VcSsdW1Un3z4q5AHRuXc1b5yajnexYe4BX/jgi1Ce/95bJHpn9YomhqKsunFr0RhA5B14L1cqNki8ymYygbyUzE5+Hb6p9wm9P/inhzzuhajSHxT+mc5v8WJZzlmafeMturLvHvubgeIAl52TEedFPHqZQNh4edmD4eX/1ptvK7TtCW67NafzZ8Txa8kG06QK17TI/zh4odVryQYquQy/Y9U+p/vX9rsb33rks+ZB765cHp1WntD2r5rY9AFisku65pmbo1DCfobHp2pL2vGDb1cbYtFvbxDPxu6oXE9tqqmKhlrwzZ04H0QPBzQidWxf+97WVfLOx5Fsembhr6kk2tt/trjMQntZxHDPrTrZyfMDQ2HTiTYUZsAE4OZ5acenFmv/jyZq/utbZbbU5apPKKD98yJPqdD+kHJZGyffeAU52WKobD4eJz7sjQpwDvjb5iK6JQjut9BKRLlGwlXyrTu71fhO4YlWwzxXQ1RNV1KaL/72zX37pfPJBKSPS/B7NVpbNW2ruZFadLpzetV0yTbRSQP893KcEWPkTE2t/yZGx91lFeLIz97FJeZYq/TsrhPatCt9Ts5W8c4yuUFWiCq7kReQAEZkiItNE5KL0R5QPH2UwRdpm9OSf+Gb20sDtzRlEFDpjgH9arrvKl9fryJJzq5IDqhvMepFMuLf6BnaL+add9nJ6PDU6YqvYLHaMTXatsy2S1kRMFJdJrHmYcghznWQaluicPHXsQ3Do7dBty9T9ijWJZqcz4Oj7YMgJ0Y+xX1CDj0quE4FTX4FfeEoU+r2swnLvz7JmIMfTKMEF49PL6UOVZ+C1ikZ+ObxvYrlZKZgz1rXPBm2qYNlsWsk6Tq16neebPJPZQhg/NxnQ0Koq+Ztu00e/ILfaKLsB//51yzk2Ptp3W63oqCL72XaOp1WkJS8iceAO4EBgIHCCiPgnoygg2Y5an/SfCG4HD6c9MIZDb/eZjm+RiSX/R0fypV03dYfYOSNf1lFFVI9eL1nI3vGveKgm2rT+BqLlafl0po7Pv7Dab4ajD5lYw2Hd/AEhqYYzRUSX3us7QkfmbHtK/s6dDbG4VtZ+L5X+e+gUyl6lbLtKRKDaSluA6Bm+3pm6mebgaVgV7bh3/uaz0pM6wW9MZvYY1+Lv4y+4fOpKkZLT3jtBsNqZlbK5mdYhtRGc+a7qapO9ZtvFs0Hb7Nw2/+Jarq++J1Gv2A+7gEoh60jYFNqS3wGYppSaoZRaBzwB5DGuTZNuduqNb3zHM1/MDt2nWGSSQsEZs9t97feB+y2odzfko2Lvh5w1wvWbkg9Ka4lmmdsDwh2ImCjOqbh3/J1728h/wA6n++/rpV3qLMmsmfOlLr13akBsdznxi+fhqHvgUnflrYQihvTRRX7bbaUa2ntK45MPvmDy45SX0+59bvUzLheGb8t1PE9bx2ZyabVjDkTTWuIx4bKqh/iu9hQOirknIzpDKH/oMoJVqpbpaqNE1E22aQ+6sBTwn41rYw/0frcgWeilUkMoNwIcOVqZba1LICKni8hYERm7cGFIDckQgsr6OTn/6WjuiUJht5dMLPmYOPZ9x0pV21uXTvtf066JTVXS7Orq/aPmrsBzrlUR/I1rkhkm2xBtMNYeUL65McBH7s0D77LkPQ9Tc5PbWvRTRofcpl0Q4HZn7HUpHPMg7OtnTabhh5CEZpVIOiUfapGnUfLZxMlnoTSd7dpXCS77MXjilxX++cuq16iRJu6suc212TnlpEbVU2P1AhLvOZ9Tfnrx3mllth9bFXKPqnwqaFVqdI3/JDbnglJ3K6WGKaWGde0anDEujK17ZTlZJiLP1VzK2Nrf0bom+wE3+0bYaYAn1p7GQ9XXpjnGcavsaAlrFuc6lXxAX/rye6q/e4lZdSeyEakvyoerr2FW3Ym0ZXW0R9ORZ7yRaN/5uwUr2Dc2lsuqH/bfwRspInG3a8FJp776L7GvTzPdbhT0tSJrjnC81HY7HwYdDsP/GEnujMhnYfBiYN/WoNQSvi6ziJb8igxnNHcfTDYvhhpHPuAtZvtM7Fs8Daa95X9wY/hEPGecfL+f3qJamthUZtOuTj9bfrdgww51fH3ZfmmkVo5/k9Q5xqv8yiRWpE8ebbk7p4n1Aubm+yLZ5FqJypT5Kxgam84GspzV65JhaSfH36SX/BRypBu7QZ37lO5RtJa17BYPyYWOpyHY08Atq+X4qncTm26ruZ27am4BYGAs1a0zIq4HwsbX/ZrNYz+mbE/BMQGmFreV9NsR/VzL1TRyQdUT/DR7hit8LQVvC47F4Y9fwW88CaxOeAIG7Afb/ya5rlgDncPPCd9+znj4jf+AWlnjfTnZvvpsLXmJwWKfWPgwdj0n49+xobYznRzhwVvMCUilEBTy2hCcomBU/HXarEm+qKqadI91WOw7+m1gR1T5yPvN01TNDC++45cZE9yuzCpS5zEUqoBIoZX8GGCAiPQTkRrgeOCFNMdkTCETDe1/S6p/ux2ruar6fh5JY4k7yaZosMunZzfkbJJEOXi4JjnTcJOubVzbttvYCrtzKPk2Hp/8ycO6O5YU28W+48yqFzi/+km2i00NvrA3J7vEdYbDjbZzr9/8QK0MnH3pYlnQ6eLSO2wEG+VQbKRQHHyz//rD/wXdBqV+r8Pv1KUE/eL9o/jkq1unj7Bx0nNbK3d9Zs9A9eBDXct1DQFjb0HplQNy9HflZ66sfpAdPkmOBdXX6jDV66r/E56F9X+/ps3T4TWDgw5fS/KexSTVkq9IJa+UagTOAl4HJgFPKaUm5Ps6merPfQd2T79TCHXWjDx7AoXNHaODy6NlY4zG/Nw1TdnXsvVy6nC3VX7QVlZcdki+mVhDcqBopz7teOJ3IwDYOTYx/GKtu8A+VzpOFOKTT6FIlnyluWNshv0y+XmrY5KftxgJv/841S0z6HAY9WKaRin+n3f5g34B21WlomDf10wfAk9SuPoqnzkL4AoUcNFYDw2ehGU0J3zv1Q0rHOuTz5ptNNpjAK2ph3/vBnPHRRLbPpf32zp99H6W/JylFZqgTCn1ilJqM6XUJkqpsJL0WZOuwst2MsWVgjTbijC1Vfp22WFa3kHJG16fEnhsdpa8422/wipzlibhV1BX0Y+NOronjsQEHSv++d2Jdc39dne5pWKNye+smhuwm3JPWcKDjSHhjCLuKBo/hdp1C/ey7U4o1uSjYl2nEBzzgJ7MdWDEDJ2B+Fnyjja1xJqF7XSnpWOdbRhk+Ax4eq3LV66GHX1yIXndNV2tuQ0/TYSf3DZlDY30ien23GqNu3SgjR1Vo4B2dVWct8XPMO/rjCtjicfvfu6+myY+x3188u9/l13gSToq1HRxE+auac9K/lt7JbdW35FY99oE/x836nVsRbqQjpGPdYsYTRHH7IFF0Dm4Ie1gUibs2moWv99jk8TyDv06Q/0y+CpZri428z0+rD0nuex4oCbPXuzqKs9TYWMjogtqJ07kUKi2QtnmZPch2/7Cvb3Q5KMYSKkYdASc9kp4yuKMCLjntusnXqXnE0RhoTVxzu939Na8deLJVFmtGvxnSjtfBlcsS85veP5MqG2fsvvjNam2ZpXjAbU/zl9Wz7dX7M+vdrIKp4RNFHNQa02s8oZQbtAm6a7xU/KFokUoeW88a3+Zy6/iOrzOdqnsHf+K7WQKHVnBrLoT6WrFsmbCleikSp1Ed/MysZqdMjqPO2WnjQOP2WrY7qkrjwxO7ARwRfWDqSu9g5sWP/74I3e+mxxAG9SzA6xe4ruvTUwlH6iLqh5PWFENKh4aF5zygPspVK9PMlFQo0jNtFLdNfmkIc08B+fLOdPxoYzvr7s9dJIV/pFC9svgl2/o/23jAFLy6NThNpJGxV8H3EaY3evu1clySdkuUkf20c3lB1+Jx1yyT0LJp5YxTDJQUgMkBnQLcEflSItq1bWso461PFNzBZdWP0I1jewgSRfKf2uv5PXaPwPwUe0fMj7/MTGtLK+qvh9IlhmLglPHOX/s0EEevwbtzUmy16WuRV+ZuvpMzQe++HEZPVnErLoT2T/2uV65bqXvvjZVS5LjDidUjYYv9EslTnP4S8/7gEdxjSS+fwRLfvBRuh5rLpRBzc+yIcq9iJLp0n3SzHZXza4Xf1up938mZlnBEXYbs3sb7XrifVGcGHcbPFdaRpHTtWJLmbgFC6zxJscL4/VanaHl6iMGE6M5YTSKgDS7UxfYOJe7ytKUr2HK/0Xgq9rfMrnutER1mBjNnFTljqHtbt3cGsm0gSZpTeYDJE6fvNPXHjoDNkrpvBHn0dBpQPg+AcnANurQikFWeoSj49aDss4RdtbGPW9h59gEOr3s8cVO0KkXYqIy64I6LXn7HngVSyaW/NH3pbp7MsYo+SQR7kWmBWEyfYkqlXoNv2diqWVVOyOy6jrqmrSe56tW/GRWrvO6xPz2GXg/eJxjxKZduaTNC4yp+z3trHkoYoVuuiYzAj0cWS+nq9RawusKlOKgRSl5ewq+nQ5aUCn52L1sKrPZRDLLx95ewkv9+eFsOM43ekOjorcsYDvxGbT1NuhF05IzXx0nrj70pvQCHJHq5hnasxX7xr5wr3T6zVclB4L2in3JQbHwXD5xn7Awp5wuXJZ8spSCi8TqAjfTrY+3rmOUfILAe+FYn3E4bxaWfINnxnWYUeQ1HFbMw2vJHxZLzSsVQ7nO67rEJMfcjxnvpV5S4LQ2nwDQTX5GRFCJ+QfJE73yxxEMcUzanKZcE/85OPYJbZvC07NkS4YZiiqDRHcLxYdNg9g1nhxhf79pK9ckpLdqLwSgb/1jRCXU9xx0TIBPXqH4oPZPAGxS/zBNzhmmXiV/uyemPHHyCDHLPnnX66b8j2Or3nPL5FOWDeC+mhvTXiLuExaWxOuTdyznw5LPhdqQyS/rLXm05O3Bz4xfogqme8aTwnq3TsNh7TKt9D0vhX6x1BQoMY+rMWlyCLR1lHpclTr5USTZM+8uP+u7Zj1DTj0xsGd7WJ50hW4uP/KD0qHc3VnC7TX/RA8XHBH8/bKkRVnySfTNjaEY3TzUteXz5mSYnp2zOlO6efxpG7ZPn8Pa6Xt3NqjRU5LWcitvmt4o7pqo+/lMfIlPSs5L26e/5dLJ2M+axBX7642aCPXBB1jydqRIXWqERF5pbeVr98u5vr7RxXL9RVHIUZX8b+0JhZ5z2gOlQahm+PKh1HU2dR3d25zGQOf+VptLb5DFPO6axOkEGHuf7zFjmjez9hFXRJNI8vm2lf+txw+1ZE/Kck9NsvedSQBHNrQYJR9zDZyoxP/VodZlfojyPIjLJ+/sGiY/u9KkQnQl32NI+n3SzVBcMkP/vyxitk5nlkiLtmHJzMKs8SBLfq9LdSKyLQ6OJlO2jDgPDr7FPZFofSWRgTOgUTt/o4aIY1O2gZHyoKRRbkrBDE8aCecz4W3TTnfNigUw/r+BPVPXYSjEERoca1hNe1bq6wfMpv2qWb8MY42rYb72DFTT5CpEknoH/b9vofuPLUbJX1rlSIxlNYQdYpNTZpZV5zDgCrCNhEzdDyHIku/kKDWW8kKKquSrIsTvpi1ZZwn4QsSoo8Gp2SZPrnpbf+i8Scq2UJfS2uX6f2+ukeo6nYis0L7yqhoYdpp74G59xU7dHOWe+7gvfAlS8ml9+j5K0Y506dzfNWYEuHuLdijop3fq/3vtEHiVGM001XZMLO/+5kF8U3c6uy4PTjl9epVOldzj9v6JdXGaXBrbNjwTBl7AeELoWFYeaBmtevk8Tqt6PWX1sfH32MKTkKuG8IbVWxZwdPw9fhF/nU0l1ardNDYHOiRzrn192X6R3sSCsCGLmVV3YiLOHtwj7imWvNeKCSLMSrdDLtO9CESgfnm064F+yILwK+RcFfKSsUPTfCauGEqFo1UPP9t/l5qINWiD2qft7rlwJpw1Fi6YrguYH3kPtOrs31OwJ1YlXEBOkX3Ume1uCelJVtFMzXIrbn3H3yVmwm655ovAY/yopskTKu1Nb+AZH5B51vbCKvmWMfC6wj+x5QHxMSnraryK1OKiA7fgulcnJwZBbbwDso0qrnNYW3RoHWHQE23Jf1ynrWTnDNId1yQba7U0utuBc2TfByWxwBfMuOb+DO3dKem3jFIF6Nnfpd/HJtuZlTufBZ/c7llpu2tahs3RInBqK3vMIoh4jeXaCDCgatr4r7eVfOvO7va09bHw5uXwjU+VsdWLoG13/2RyYe0nZEzInhBlnSTxaduVYcV3UvEGHthnSgRdeHrmo2vPo2/9Y1kFcmRCy3iqMihjlmItW/xu901oV5v+PE0++dWjVJDZo/lT3/VnLkqGRL5bex6z6k7kmZorog0KVwUP+FbRrO+LPZAqoi2kMPws8CCyzfOy/9XBU9lNCGN5EvR8rbN6pGd+Dic+mbr9imVw+dJgSz4syZm3LTgnAQYp87A2GfICOL/akac+qovUhyqa3BMdE+4a+9z+ytwo+ShkkHMkSMkDvH7Orinrqjz795Roxb3XNboby3FN0YttD4t9F23HEBdMFU2w7SgY4nhZhL0MV8yDZZnNF8gbhaqWYIjOb96BMwIqY6V7+a75OViJeo/95Rsw7Fd6oLvvbsHn9J5vz0v0/3UdQ67l0AMd+4SfL5BobXGh6pDiZqmiScfJ994JSCrvIHeNTTJopDBGTstQ8hlY8nW+M96A5mZ6fntnyuq/Vd3vWr4oYqHqy553V6wvSEKikIRJ1dIE25wEQ45LrtwgZGasak5aZkUnIITSUDw22g66D0ouO5XzBpslP/faPvXYqtpUpRqUvKzPjnDwTREGuj1twXb51C8NUfLOusGebJVRe54RLfmusowz48+71lWJtuTlR91rT/jk01jyiUCMAvVkW4iST/8Dqr11mtBd+wblpF6XOpsUKz9LGMvncv26q1ypjMGdNnTZmgYaGvMfyikhETMd8MlBY5UOzDve+PI/pyZfCqXYKYUNmdFzm+TnoT6pIyQOXTaBP03QA6ejXoKTnsntml5954zOClKGzvYz6HDP+SKqunXRZ7O73DzABVVPEh//VPKSXuPFGw1kkTAACzQmtX4o+ao6pHM/ALrWBijbgMHbtLx7HcObv+DQ+Ceu1asbktf535ezfWs65kzQgBaFKwrsy+896Q5adczs+L0vgx1+C1sfl35fQ3FottprK8eAaF1Ht/V9jJXxtOvm+v8OvaBNF+g3wp0eIxtSEto5loPce053TYoLN6KVnGk+HgddZTl1LyZ7EDGvJf+JfwK9hE/e+TLNIy1Eyadx1zTWJy2BoKK/t2V5g7/UDX1DWQzoWbQXVT3OBm2TrpSYSG7umi0P8V8fD/7ehR7MAZJRF+175HaeVp3goOujxfsbioOtmboPSj5f3qIugw7Xg6sFcTN4cx052voyd1i0/+FZqragUoJZEFQhykviWd31T+E7ZknLUPJRBl4zqUmZBfUq6TrxFrQWydEnH2i5BP98aa8XpQDCPlf4r7envp/xEZz2qmtT8+YjkwsH35L+GobypHVnOPm/cNwjOq3EKc/CCY8X7/o/z3QvRxl3cypo77MRVen7KfkNNo92rIekJR+u5uvsdCY5RPaEy9ESaIwwvbrASl55bqX3Z80pP4U3R0fipMEvt9bVaX7atDNg8Z+5CnDILfr/dhvCxru4NsWmvJxc6NTXOk/IxClD+bLpPknX2yZ75bHqVBaEGXJbHKwHfp3ypRSpidjb8EuDcOTdqeuctOnmuzoRQplY468Dfl9l5ZByVGTLJy1jMlS6H3Dzg6JlasxFBJr5Y/x/vtuWrW7IzScf9P2q/fPEA9Q0pomU2WBTmPtV+D6rF/uv77Nz+HE2Ju7dkC/8xt1iVdryPv7R1G0plntUJe9jydvGShBp0jskHoOA3sjeces5XJemKleW5GTJi8gxIjJBRJpFZJhn28UiMk1EpojI/rmJmYZO/cK3x+KZW/KjgvNW+F4CxbnVyYgCp36rqcqxwxSkLI/ymdy0/zXh5zr0nzpW+fC70l/364Bw0ahdXzOD1ZAv7Jw6Tn77QXB7T3HXRFXyPrN2s4z6SpkMtdkB4QfUL83qOunlyI3xwJGAa/6viAwEjgcGAQcAd4oUsEpyuh8wVpWZkh90BPQbwfw90udQt6kVd+NwZaMTaM72Vle3CVaW7R3VZVp11pWcNtkr/Hzb/kLHKnfbInw/gMFH+q8Pu99ON5Att5nsZMgVR76oBN0Hws5n+u/vfWbmf+tePj5gfMHPXZPBPBwnW26ow7Ulai9ik72zuk46clLySqlJSimfkkYcBjyhlFqrlJoJTAOC08Dlg4tDUuRK3B0Klo6jdFKjus4+DSsAb7bLlPoX2U706dw/mkV8wTQ4bwp0DC4MnjGDsihgsMfFyc/Gkjfkgy6b6rGBS6PNNgdS295KR7GQA/4OWxwEAw9PPc7PXeO0TzOQYc81b7KlfE+saQ2Mvja1yhVwS7Ujj1M6Ay1LCuWT3whwJmuZba1LQUROB04H6NOnj98u0QjJ40KsKrNSZVZMbrw2eqyvX81GJ81Zz+ZURPIn2l3KKAOqaREYcjy07pL5oc4eU+JBM5a8IQfs8bRMeuNeJR+vScbA97bsTU8NYyBAyTuevwxk2Kv+TfaqfZPvJ3wPE+/yrftweNyZSqJEhbxF5C0RGe/zd1jYYT7rfL+BUupupdQwpdSwrl19bnpUwrpUsXh0v5ojeVZ1bUgCJQ/Hxt91Le+1LjlTVhCasu00qebMLOJ8WM8jzoMj7srOF+n8HYy7xpAP/DJOpsP7HPTfI/l5o22t8/rMfp+dmrk213QbVY3WgKqPJe++TGE82mk1glJqH6XUYJ+/50MOmw04fR29gCynlEYkzE8ci2flV6uqCY5e8bJDzO21unC1u7h21pa8Usnv1t63M+SmVBEt/fe0BXCsNNE1hhzZ969w7EPp9/PiVfL2s7PzWcl1u10Y7Vwhkw6jUF2/RH8ISJkwX1kZNvumJkjMB4Vymr4AHC8itSLSDxgAfF6ga6Vn6EnBSj4oFhyoqnFY8l23zOiS30pyAoUIbCpZvuOclnzQIJOTbJR8rSf3TNg5Ntzaf71f8jPjkzfkyvCzs5tR7W17U616ss5EahkYcbnQ7Qdr7shy/3HDDeVn/aEcE5SJyBEiMhvYGXhZRF4HUEpNAJ4CJgKvAWcqlUOF6Fy4fCn02SlYyZ/2SvCxzu7cb9/L6LLfNXbl3g+Ts/aWqCy6nECKTz4oF3sudOjlWRFyvd99EH4uV2mcRPq9bCUzGLIjFoNN900ur7Ri2TM1PLqEZG6tEHKNrnlWKdVLKVWrlOqulNrfse1qpdQmSqnNlVKvhp2noNiKJkjJhw3YOgcxM86rIvztpYmJpa6SRjlvHNBVU0pb8D2GRk/gtfnI4BAxP0ZGDxUNZJc/6ARLztqvxpI3lBLnIOmuf9Ltc1BAWHAQdprlrY5NX3SnTGkZM16jEDSIGOar9yr2DbeG+d9EulwfWeBa/rJ5AMPjE4IPOPxOuNXHFaKaoGPvzHoSJzyWfh8nntQEWdGxD5z+rntdYuA199MbDBnjfOY790ttn1GwnyW/iYf5ZIuDC3bq9cfUClLm8Wo4++uAbV7rPbq26iFLXMs/KCu/hbMA8Zlj4IxP4PT3oNPGqSl7oTSRKfnyDdrnMeOvhnxx3ne64HcUnD3JfPYq24WHS4cT8DBMzmyGfSYYJR+r0rkp+u8JbT1Tp72j6hko3DUq+YIQEQSlR9G7bJrcqetmetZez6F6Ocos1KKQLyVvNy+j5Q0Zssne/om/2nWHNmkKiydwtLt8FqTZNIdJS1nOns2F9chdE6TkrR//F8+lP4czFei+f4U3Lwvc9Y3m7RKfO6z+nqOr3tULIUnFfEkZFK0gjE/ekC2n+Cf7ywhXEEAe22IuEw5jVf75cQpIy1XyAw+H4X9MLufjTe6y5MOt07WO/PK9Fn3oOCwDq/bo+92TOErNb97JbH8zGcpQSpyKPcuc8L7kktE2Slr0PNNyTa2eQ3VxYpu8+JkdyiqD82V96cFHBufwDopXzwdege0ars77GYWlVgWfZT/kLpPBkCkuJZ8mFDKTqJvVGeTQSaH4Bk/LVfKFwGmRprFO41KYKi8AnDsZfvla4c7v5Zxv4fxpmR+Xbd1cgyEf2Ep+78uCLa2DrPDhTKbxtPNMzvI+G/v+Lfz4C2eGb88zLUvJH35XMsF/QVwEjnOu+Tl0z4IU7rZp3yO0iHfu+FjybbPIK5RIdWAwlAKrHXuVspO21uBuJqX3nON7Va1SC9fXtQ8/vsgVtlqWkh96gn/60HzhbAh2trqd/FMN/KHqOX2IUtQ0rCycTIUgXyGU9mSUgPJoBkNBiTQmJBH2CTgv+I/11fgkPvPSbWD06+VIy1LyLgpgybuUvNW9S5O8aG1jM41RimaXGmf8ftgs4Eyw0zvnJf2xwZAhtrESZqUnXgQZWPKu+Ps4rp7vNifDxsMzO0eBaXlKvpBZGF0+eUvJp4l7XbOuyT8Op+c2eRMra9o7wjO7D05+7jUsdd9ssF1KA8OyUhsMBSKKkrfnrfQdEb3H6VTQWxzk3rZmaVpXLgCbFbYiqpOWG0JZaJ/8oqn6f2/x3VOehYd1RaX7qq9n4rwdSbHjL5wJ1QG56jv0hmU/5kXatJz1ebKQgrPbmY80B6B9j+dN8S/OYDAUmihWerct4NxJ2m//+sXB+zmpcvRM97vKvW3Wh7Augnt2z0ugoR4+vSPaNXOg5VnyYfHruQ4EOl8c09/W/9spTEEn5+q3B7NrtXWwV3yc/7umdedgJd/NSmm8/W9ykzUKNW2gVafk8r5/dbtt8kG7DfM729BgiEpUV0z7ntrqr+uol3c4PXx/5zictydf2x66Rpi5Hou7azQXkBao5G18tOvmB6Wuy4TlPiGBtpKsqoOj74VYjIVtkxMvTr73M2KZhGfZb4UB+4bvVwiGn+1bosxgqEgy9bfb7p3dLwrfz5mH3utbb7dhZmUKi0DLU/KhFaKsr9tvd93N+nWGMzhthf6L52EXazZtn51TdvMW7W6MZTPwavK9GAw5YZfTi+q6ddYk/u0H8MevIhwjbp2zyx90L/24R9MfG8V3nwdanpLf5mRdgHqrY302Wj9G5376x+iV4QzOQ27VAzUbD4dhv9S+Zp8c771XjXctr1iXwci9n4/cYDBkzg6n62d0y4hpfJ2Wf4+toXP/9Md4c1HZVv6WB6ev/+C0+I++L5qMWdDyBl4794cLZ/hvs29qdZYTibY4KDma3rkfXDAN6pfrZYe10Cju7tr7Uxezd9QenN3QbP+gwWDIjq6b6Wc0KrXtYNXCzK4RiyfDqcHtow8ad7OxS48OPkr/FYiWp+TD2Pp4WDITRpybv3PaXbWmtYlVPdZMTXweKtMyi9g//E746pFkRXmDwVAcTv4fTHw+OQs2Kk53kHOOyYjz4YsHwg60/i+sa7bluWvCqKqBfS7Xb+yodO6vR8yDaE4dVH2vxy8TnzeP/Zjiow+lfU/Y/cLCxvsbDIZUOveDXc/J/Djn4KvTBVObZuZrj6H6/wLPI1m/LPlsOOuL8O0+I/fNkvSnV9OYmZI3GAzlz2U/J5/9mHcGrP05jQ3ddTO4bEnBx99ysuRF5AYRmSwi34jIsyLS0bHtYhGZJiJTRKR407vyTSzm/hG92Gl4t/1F8hCH4l+oOhglbzC0NGIx/5QmsQyUvHf/ApGrJf8mcLFSqlFE/g5cDPxZRAYCxwODgJ7AWyKymVKZBIxXCLE4XLHMtapr2+RtXU4burCi2FLlnYaGBmbPnk19ffGLHhhKT11dHb169aK6urxiwMuOQtWVzYGclLxSyjHdk0+Bo63PhwFPKKXWAjNFZBqwA/BJLterFLbs3gassddaGpASFArIN7Nnz6Zdu3b07dsXMeMF6xVKKRYvXszs2bPp169fqcUpb8pQyedTil8Cr1qfNwKcCVhmW+tSEJHTRWSsiIxduDDD8KUyRew0xGiffEugvr6eLl26GAW/HiIidOnSxfTiIlGgurI5kFYKEXlLRMb7/B3m2OcSoBGwp3n5aQJfc1YpdbdSaphSaljXri0kkZXDJ39PzU2FLSBSRIyCX38xv30aNtgsdZ1TyWdSXjDPpHXXKKX2CdsuIqOAg4G9lUoEjM4Gejt26wWsP7XgPGGVrVgbsKPBYGgRHP4vGH01dNkkuc6p5NNViyoguUbXHAD8GThUKbXasekF4HgRqRWRfsAA4PNcrlVRNLtdNNW0vPFmg8HgoNcwnWbcGSfvqiBVugHrXJ1GtwPtgDdFZJyI3AWglJoAPAVMBF4DzmyRkTVBeGLna2VdiQRpOSxdupQ777yz1GKk5ZZbbmH16tXpd8wjbdtGKDcXwo033sgWW2zB4MGDGTJkCA899FCeJFvPcbq4SujuyjW6ZtOQbVcDV+dy/orFq+TRZfBUvKZFRMxf+eIEJs5dntdzDuzZnssPGRS43Vbyv//97/N63UxRSqGUIhYwd+KWW27h5JNPpnXr1r7b/WhsbKSqqjTzEu+66y7efPNNPv/8c9q3b8+yZct47rnnSiJLi8a26lsVt4g3rG9pDYrFDqe7Zr7ZSn7pLz8ulUQVz0UXXcT06dMZOnQoF1xwATfccAPbb789W2+9NZdffjkAs2bNYosttuDXv/41gwcP5qSTTuKtt95i+PDhDBgwgM8/1x7DK664glNOOYW99tqLAQMGcM899ySuE3TeLbfckt///vdsu+22/Pjjj5xxxhkMGzaMQYMGJfa77bbbmDt3LnvuuSd77qkL1Dit7GeeeYZTTz0VgFNPPZVzzz2XPffckz//+c9Mnz6dAw44gO22244RI0YwefLkwHsxc+ZMdt55Z7bffnsuvfRS1zY/+QH+9re/scUWW7DvvvtywgkncOONNwJwzTXXcOedd9K+vfYZd+jQgVGjRmX+AxnCsfPhHHBt8a9tWybl8LfddtupFkPDWqUub6/U5e3V6L+MUOry9urnOdNKLVXWTJw4saTXnzlzpho0aJBSSqnXX39d/eY3v1HNzc2qqalJjRw5Ur333ntq5syZKh6Pq2+++UY1NTWpbbfdVp122mmqublZPffcc+qwww5TSil1+eWXq6233lqtXr1aLVy4UPXq1UvNmTMn9Lwioj755JOEPIsXL1ZKKdXY2Kh233139fXXXyullNp4443VwoULE/u1adMm8fnpp59Wo0aNUkopNWrUKDVy5EjV2NiolFJqr732Ut99951SSqlPP/1U7bnnnoH34pBDDlEPPvigUkqp22+/PXGNIPnHjBmjhgwZolavXq2WL1+uNt10U3XDDTeo5cuXq44dO0b+DUrdBioSSweoxgalvnpUqaamglwGGKsC9KrJXVMoHAMwe8S/1h9MGFpeeOONN3jjjTfYZhtdDH3lypVMnTqVPn360K9fP7baaisABg0axN57742IsNVWWzFr1qzEOQ477DBatWpFq1at2HPPPfn888/58MMPA8+78cYbs9NOOyWOf+qpp7j77rtpbGxk3rx5TJw4ka233jqj73HMMccQj8dZuXIlH3/8Mcccc0xi29q1wRFZH330Ef/9738BOOWUU/jzn/8cel9WrFiR+L4AhxxyCKANPBMaWSTiVTD0xJJc2ij5QiECG24N879JrgrLgWOIjFKKiy++mN/+9reu9bNmzaK2NlmFKxaLJZZjsRiNjcmoJ69yE5HQ87Zpk6xBMHPmTG688UbGjBlDp06dOPXUUwMnCjmv493HPmdzczMdO3Zk3Lhx6b56oPwQfF9uvvlm33O0b9+eNm3aMGPGDPr3j1Agw1CRGK1TSE5/17Uo5nZnTbt27VixQucA2n///bnvvvtYuXIlAHPmzOGnn37K6HzPP/889fX1LF68mHfffZftt98+8nmXL19OmzZt6NChAwsWLODVV19NbHPKCdC9e3cmTZpEc3Mzzz77rK8s7du3p1+/fjz99NOAVtZff/11oOzDhw/niSeeAODRR5Nl5oLk33XXXXnxxRepr69n5cqVvPzyy4ljLr74Ys4880yWL1+e+G533313+M0zRGefK+HXb5dUBGPJFxJTwi9vdOnSheHDhzN48GAOPPBATjzxRHbeWdfXbdu2LY888gjxePT7vcMOOzBy5Eh++OEHLr30Unr27EnPnj2ZNGlS2vMOGTKEbbbZhkGDBtG/f3+GDx+e2Hb66adz4IEH0qNHD0aPHs11113HwQcfTO/evRk8eHBCAXt59NFHOeOMM7jqqqtoaGjg+OOPZ8gQ/6Lqt956KyeeeCK33norRx2VrCi03377+cq//fbbc+ihhzJkyBA23nhjhg0bRocOOnvqGWecwcqVK9l+++2prq6murqa8847L/J9NKQhm/z0eUZU1CK3RWDYsGFq7NixpRYjv1zRIfFx+ZkTaN+1VwmFyZ5Jkyax5ZZbllqMvHDFFVfQtm1bzj///FKLUjRWrlxJ27ZtWb16Nbvttht33303226bWfWxltQGWhoi8oVSapjfNmPJF5rOm8CS6QBIzAxyGUrD6aefzsSJE6mvr2fUqFEZK3hD5WKUfKHpPjCp5FvEVKjK54orrii1CJG4+uqrE356m2OOOYZLLrkk43M99thj+RLLUGEYJV9wkopdxWtD9jMY3FxyySVZKXSDwYkJ9yg0k19KfFQmJtlgMBQZo+QLjU+hb4PBYCgWRskXmkP/mfhYRoFMBoNhPcEo+UKz7S9KLYHBYFiPMUq+iHRoVVNqEVoUV1xxRSKbYqmYO3cuRx99dNr9rrnmmozOm+t3++677zjooIPYdNNN2XLLLTn22GNZsGBB1uczVC4musaQOa9eBPO/ze85N9wKDrwuv+csAj179uSZZ55Ju98111zD//3f/xVBIp0jZ+TIkdx0002JZGSjR49m4cKFdO/evSgyGMoHY8kXExNdkzNXX301m2++Ofvssw9TpkwB8M3F3tTURP/+/VFKsXTpUmKxGO+//z4AI0aMYNq0ab7nD8o1r5TiggsuYPDgwWy11VY8+eSTgE5eNnjwYAAeeOABjjzySA444AAGDBjAhRdeCOhc+GvWrGHo0KGcdNJJGX23oO9nr99pp53YfvvtueyyyxK56x977DF23nnnhIIH2HPPPRNyGtYzgnIQl+KvReWTd2LnlK5fUWpJsqYccomPHTtWDR48WK1atUotW7ZMbbLJJuqGG24IzMW+//77q/Hjx6sXX3xRDRs2TF111VWqvr5e9e3bN/AaQbnmn3nmGbXPPvuoxsZGNX/+fNW7d281d+5cV577+++/X/Xr108tXbpUrVmzRvXp00f98MMPSil3XvlMvptSwbnmR44cqR577DGllFL/+te/Etf405/+pG655Zas7nEY5dAGDP5g8smXmJq2sG4lYMJrcuGDDz7giCOOSJTWO/TQQ6mvrw/MxT5ixAjef/99Zs6cycUXX8w999zD7rvvzvbbbx96naBc8yeccALxeJzu3buz++67M2bMmJQc8nvvvXci+dfAgQP5/vvv6d27d1bfDQjNNf/JJ58kSvWdeOKJ61UuHkN0jJIvCpabxsRQ5ow3j3pYLvYRI0Zw1113MXfuXP76179yww038O6777LbbrtldA0713wUnPns4/G4K4d9OvxyxGeTa37QoEG89957kfc3tGyMT74YJB5eo+RzYbfdduPZZ59lzZo1rFixghdffJHWrVsH5mLfcccd+fjjj4nFYtTV1TF06FD+/e9/M2LEiNDr+OWa32233XjyySdpampi4cKFvP/+++ywww6RZa+urqahoSGj7wbhueZ32mmnRIUoO788aKv+448/duWNf+211/j22zwPlhsqgpyUvIj8TUS+EZFxIvKGiPR0bLtYRKaJyBQR2T93USuYjhtbH8zAay5su+22HHfccQwdOpSjjjoqoawfffRR7r33XoYMGcKgQYN4/vnnAW1V9+7dO1G2b8SIEaxYsSJRHjAIO9f8TjvtlMg1f8QRR7D11lszZMgQ9tprL66//no23HDDyLKffvrpbL311oEDr0HfLez73XLLLdx0003ssMMOzJs3L+EmatWqFS+99BL//Oc/GTBgAAMHDuSBBx6gW7dukeU1tBxyyicvIu2VUsutz38EBiqlficiA4HHgR2AnsBbwGZKqaaw87XIfPIAK3+C7z+GQYeXWpKsWV9yiVdSrvnVq1fTqlUrRIQnnniCxx9/PPECKATrSxuoRAqWT95W8BZtSPojDgOeUEqtBWaKyDS0wv8kl+tVLG27VbSCN5QnX3zxBWeddRZKKTp27Mh9991XapEMZUjOA68icjXwC2AZsKe1eiPgU8dus611fsefDpwO0KdPn1zFMRgicf/993Prrbe61g0fPpw77rijoNddvHgxe++9d8r6t99+my5dumR0rhEjRoTWgjUYIIK7RkTeAvycj5copZ537HcxUKeUulxE7gA+UUo9Ym27F3hFKfXfsGu1WHdNC2DSpElsscUWvhEghpaPUorJkycbd02ZkpO7Rim1T8TrPAa8DFyOttydwcG9gLkRz2MoQ+rq6li8eDFdunQxin49QynF4sWLqaurK7UohizIyV0jIgOUUlOtxUOBydbnF4DHROQm9MDrAODzXK5lKC29evVi9uzZLFy4sNSiGEpAXV0dvXpVZhH69Z1cffLXicjmQDPwPfA7AKXUBBF5CpgINAJnpousMZQ31dXV9OvXr9RiGAyGDMk1uuaokG1XA1fncn6DwWAw5IaZ8WowGAwtGKPkDQaDoQWT04zXfCMiC9G+/WzZAFiUJ3GKhZG5eFSi3Ebm4lGJctsyb6yU6uq3Q1kp+VwRkbFBsaLlipG5eFSi3Ebm4lGJckeR2bhrDAaDoQVjlLzBYDC0YFqakr+71AJkgZG5eFSi3Ebm4lGJcqeVuUX55A0Gg8HgpqVZ8gaDwWBwYJS8wWAwtGCMki8QUoGpGitRZqhMuY3MxaMS5c6nzEbJ5xER6S0i3QFUhQx2VKLMUJlyG5mLRyXKXSiZjZLPEyJyGPAucJOIjBaRjdMcUnIqUWaoTLmNzMWjEuUuqMxKKfOX4x/QHRgN7Gwt34TOqb9TqWVrSTJXqtxGZiN3KWU2lnx+WAT8DLQDUEqdC3wEXCIiG0FZ+gUrUWaoTLmNzMWjEuUuqMxGyecBpQuifAwMFpGu1rq/AzOAB63lsvILVqLMUJlyG5mLRyXKXWiZjZLPEhHZS0T+LiLHiMiWwBPAnsARItINQCl1NjBPRNqXUlabSpQZKlNuI3PxqES5iymzUfJZICJ7Afeiu1jbAP8CNgPOAQ4GThORQ0XkJGAYEC+RqAkqUWaoTLmNzMWjEuUuusylHnSoxD9gFHCZ9bkNcAAwHv0m7gmcCTwJvAoMKbW8lSpzpcptZDZyl5PMJf/ClfgHnAC84Fl3IPAWMMCxrl2pZa1kmStVbiOzkbucZDbumixQSj0OtBWRRxyrPwK+BQY49ltRbNmCqESZoTLlNjIXj0qUu9gyGyWfASLi9I0dDHQUkUcBlFLLrfVlVVmmEmWGypTbyFw8KlHuUslsUg2nQURiAEqpZmu5G3CCUupWEWkFvAisBn4E9gdGKqWmlEpeS8Zq9G+7zlruDhynlLqtXGWGypS7QmWuuDYN5l5nS1U+T9bSEJHDgeOARhF5TCn1KtAEzAdQSq0B9hGRkUAr4LYyaFQHAb8AWonIXZbMzcBCKE+ZoTLlrlCZD6fC2jSYe50TpR6EKNc/YCtgIrAfcCJ69Pu3QFWpZQuReSQwDj1afyrwDtCq1HK1RLkrVOaKa9PmXuf+Z3zywXQGZiql3lBKPQb8DjgMOMXeQURGisg5JZLPhYh0QI/Q/0kp9RrwNdpq+I2I7OnYr2xkhsqUuxJltqioNg3mXueFUr/xyvUP2AC4HxiO9fYFRgCTgP2s5a5A/zKQNW793936vwPwHnCn1aimA4eVk8yVKnclyuyQvWLatLnX+fszA68ORGQboBb4SSk1Q0QuQU9OuBn4XinVICKnAtsB5yidc6KkiMiBQD/gX8r6MUVkQ6C3UmqMtfwL9Gj+ycoatCo1lSh3hcpccW0azL3OJ2bg1UJEDgH+DkwA1ojIDKXUFSJyC/An4L9oX2AdepZac6lktRGRA4DrgT8ox9taKTUfa3DHohV6CnVDcSX0pxLlrlCZK65Ng7nXeafU3Zpy+AOqgeeAI6zlzdH5nG+1li8AHgDeRE9YGFoGMm+D7q4eby13RFs+fT37nQyMBQaVWuZKlbtCZa64Nm3udWH+jCUPKN2NmoZlESilpljdqvtE5BKl1NWiM8FtAcxRSs0pobg2tegGs1JEdgb+BqwC1onIJKXUZSKyC3AocKpSakIJZXVSiXJXnMwV2qbB3OuCCGj+9Nv2F8AcYDNrWYCt0YMnPUotX4DMe6FTlE4CTkd3A3cEnkZbE3GgY6nlbAlyV6jMFdemzb3O/996b8mLSFwp1aSUesga2HlCRI5TSk0FvhGdxL87MK+0kiYRkSqlVKNS6h0RaQA2VkrZeTA+s2bZtVV6YGdpyQT1UIlyV5rMljxSaW0azL0uFOu9kldKNVk/Dkqp60VX2XrXil/tBfQBFpdOQo2IxJQ1NVop1Sh6SrcopT4QkY8d+x2OlnlBaSQNplLkrsR7LSL9lVIzVHL6vC1z2bZpL+ZeF4b1TslbEyg2RMeuPmytPgmdO+Ip64eaDQxEdw1PVkr9WBppNSKyH7CbiDQppS63Vo8CZgJPW5YNInImunt7klJqdmmkTSIiu6Lvda1S6lFr9anosmZlKXcl3msr3PBGETkWmKi0v8Apc9m1aajY9lF597qUvqJi/wF7AHOBc4HPgXuAro7t4tk/VgYy74P2TR4HzAIuDtl3X2CLUstsyXIQeir3peiBtJvLXe5KvNfAEGAasHfA9ljYsmkf68G9LrUARf6RrgPOtz5XAw+hZ8/1cOxzAGUyKGXJ+CTwG2t5P+Ai4Gis2YDlJrMlz6bAGGBXa7k38Dzaaos59isLudGDZJV6r3cCbnLc5zMsI2ZjkjMty03mimoflXyvlVoPcteI5SSz+AIYICIbKqUagN+gp0pf6dhnG/SEhZJhy2zJ+CWwmYgcCjyC9vFdANwhIm2tQ0ouM6Tc6+uVUh+KzqG9CugGdFNKNTv2244ykFtpGtBx12V/rz33eTkwzAo3fALogi4jdz7JAhTbUgb3uRLbR6XeayctPq2BpdDni0gV+qH9G/Aw8L5SarWI1AAfAtcppf5XSlltbJmtzzugY4K3Bn5QSp0lOq/2W8DLSqnrSyiqC4/cMfuBVUopEXkYuFIpNU1EtlVKfVlicVOw/PHDgaHAj+V6rx1tuk4pVS8i16DH15qVUhdZbf0BYKpS6srQkxWRSmwflXqvnbRoS150nuZnReQ/wF+BFcDjwB/Qg2s9lc55MRpoLJ2kSRwy3ysiV6MnT/wFuAWYJyJtLKvzBfQEEQk5XdHwyH0t0N/aZFfD6Qi0Fl2B/kkR6V5q2UXkIEseAJTOGHg5cBMwtxzvteM+3wNcLiI9gaeAzYAdRWQrpVQjupxcXETiZSZ3JbWPirzXXlqsJS8ifdG5Ik5DK/DdgUOAw4HtgWPRvti56OiavVXpiwz0xS3zbsCR6IkWq4HbgVfQaUxPAo5SSk0qibAOfOQeARwF/MKWT0TuRcvdFfitKvFMRRHZCXgfLe9vVTLSChHZEe2Pfx3dJS+Lex3QPg4DjgHaAr8G+qKjU45GVxmaWApZnVRo++hLBd5rP1pkCKWItAFWAm8qpd6zVn8keoLF/9A/1tfo0fKhwF5Kqe9KIatNiMzrgHvRhRPuQ8s8gDJQOpBW7vtF5CSl1HRgEfpFe7BSanKJxAVARDqiLcmR6ARXT1nugwcBlFKficgL6AG1srjXIfe5ET1gfCTwZ3R73hKdvXFGKWR1UqHtoyLvdSClHvnN9x869ejtQA/gY+D/HNti6JCtSymT8KaIMl8GXOJYFy+2jFnKfQlwubW8P2WQ7xvdm/sH2kLvaq3bDW2RneaU3/q/5FWTMmjTJZe1BbSPirzXod+p1ALk+QfaDx2Nsou1vAm6BNcfHfvsA/y71LJmIfPdpZa10uV2yLyzY53tstzDUvQj0a6869E9XSmFrFnc57tKfX9bUPuoqHud9nuVWoA8/kD7o6cP/9ezfhgwxXr7Cnp22jtAeyNz8eQuA2WZIjPJykP2/5uh/a9zgcGVeJ9LLXNLah+VcK+j/LWIgVcR2R890ekuYDBQr5S6wLG9H/AftKU2HDhOKfVtKWR1yHQAcC0VJDMkpnVfQwXJHdY+xJGnRkR2R2c63EOVeBCtEu8zVGa7rkSZM6LUb5kc376C9pPdSLKLNQwdt3q9Z982QGugS5nIfS3JGX/bl7vMDpn+BgyvhHudpn383bFf3Po7FitVbKn/gKuBHSuwfVxbKe3D0UauBUZU2r2O+lfRlryItFNKrRArXbC1rgo9ceiPwAKl1J+t9Yl9SomI9EeX/vpe2a2szGWGhLWzkVLqXse6spY7k/ZRLojI9ujUtDOU1Zso9/sMICJboMcwZiqlVlnrylpuEdkU7Saa6lhXA2yFnktTdjJnQ8VOhnJMVNjD8QDHlJ6c8DVwG9BFRO4AnVK4dNJqLJlfQHf93hadopRylhkSM0GvRiftSlDOcmfaPsoB0XVCHwJOAe7ytI9vKMP7DAnX0v/QFvEHYqWAKPP2cQR6HsQdIvKiiBwvIhsoPTmyLGXOlopU8iIyBB0zPg04x/KlopLTpJuAccDd1v7dSiWrjeiSZf8AfqWU2gf4AbjV3l6OMgOIyAj0BKyTlVJvi0hHEdlQ9HR/W+6vKSO5K7R99EAXfD5eKXUc8D16VmVHEWnlUPRlIzOAiOyGbse/VkodAsxGx48DifZRVnKLSCvgeHTq4v2AF9GVp04QkS6Ol1PZyJwLFemuEV0v8UjgZev/Q4AblDVxwXqQlejkR1VKqbWlkzaR5Gh3dHa6x611PdAZ7U7w7FsWMttY1uRk9OSP+9HWz0p0Lc5bgdeVLrwSA6rLQe5Kax+WTG3QSa8eQd/jyeg47Tg6d86DSqnl5XSfAURkKNBO6SIfG6EV+pvW5ocoz/ZRDbwEPKuUustadySwK/CBUupZa13ZyJwTpR4UyPSP5IvJDnvrjC4o8DI6KgJgI6Cm1LL6yO7MXb8R2prsZC13tL9buf1Zss4H1gGnW+vORVv4bUstXwtqH6cAbwOfYk1+Q0+jf4oymCiURvYqdMZOW+4/Aa9RpqGG6JDJf2ENyFvrzkInoiu5fPn8qxhL3vKxHoEuljtaKfWuY9sGaIttL2AJOi/1KcoaACoVIrIPOlfOUuAlpdSn1vpqdGzwq0qpHUTkFLQV8UdVBlaDiGwNNClH/hCr53GsUupWx7pX0Pn5S56zowW0jxeUUp+LTrlwBfCeSlqUz6Mn8L1SGkndeOR+RSn1sbW+Rmmftr3fy8CflVLjSyGnE4/ML6MnPV2ITgv8suM7vA6crUqcWiGfVIRPXnS63ZuAd9EFcZ8RXX4LAKXUIqXU3ehwqCOBv5bBAzwSLfNkYBlwrYhsAjpPvFJqMTBBRC5Gj+TfUSYK/kB0D+MMEdnWXq+UmudR8Mehp34vLLqQHlpI+/i7iAxQSi1Fu2d2FJH9ROQwoB+6ilLJ8ZH7ake7dir4Y9Ht46dSyOnER+Zr0D3n/6IT//1SRP4oOgNmP8qw/m0uVEqCsu7AZ8qq3C4i04FbRaRZKfWMtW5/YBdgn1JbDtZAzaloy/xdEWmHrobTwdpuv1z3A/ZGy1zSBGmQGJDaHvg/tKzHigjKkdvb8mOfgM47crRSquRKnpbTPtpZu8xAD16ei65YdbJS6ocSiOoiXbu29mmNdjFdiO75lVTJB8g8AB0O/KWILEMnqzsbnYr8uDJp03mjUpT8D0CDiPRSSs1WSr0pImcDj4nIPKXUR8BnwG5KqZmlFRXQGfX+gy5xhtKx2jGs3BgqOcPy78Ab5aDgAZRSa0TkIaXULOvhuAw4xgo9HGvt0yQiS4DDVYlTMztoKe1jf3T7mAhMtAZja5RSP5dOVBeh7drapwld6enIMmkffjILyXs9B5gjIp9a2xtKJmmBqAh3DTABqAH+T0SqrOiIN9GzL3cCUEotLYcH2JKtWSn1ulJqldWgQIfENVr7HCoimyqlbis3359Sapb1/0/o+9sOOFpEeojIKSKyjVLqlTJ5gG1aSvtosPY5VHRBilXlouCjtmugr1LqmXJoHxnc68GWC7XFKXioACUvyQksv0Z3s/6J9puBVkAbl0o2P1TqSLZ9j38EFojIwbhrypYlomf4LUArykbgMbRfc13ogUXCfmArqX1Y9zRK+7gCWFNM2dKRQbsum0lDGdzr+mLKVWzK3l2j9ASWGqXUWmsA5R/AZSLSAf1QnxB+huIiydqVMcuKsBt9HLgDmIT2sU4rnZQasepWeta55FZKLRCR1cAgYPdSR9JYg8FzrBdQRbQPEdkTmK2UmmpZlypN+zilHNqHk0pq114q7V7nHVUGcZz2H7AF0M+zzg7z3Bc92CfolLD7oLuGpZZ5b3T92IuwYplJxmjviZ6EA7oa1WRgQKlltuTZA/jEuuf2PbYLZeyJlaAJbQ1fA2xTBjIfAnyBftnY62yZ9yvj9rEGPcGpvbXOvt/l3D72s9r1ubZMWIUyylXuNPqjLGUuyn0ptQCOH+NgdOKuq4EtPNsGAZ+jR75LLqtH5nHoSSB/B14FNvDIfLS1XAX0KrXMDtmPREcTXA9s7ngYXHLbspeBvL2Bb0lmOBSHzIPLtH0ciB6Q/D16dvCWjm1l2z6sF9M4dG3bP6Bjy+0sjVujBzHLSu6I+qOsZC7WX1lMhrLCmq5A+8baoEfEn1bW4I3o/CnrlK69KaoMhBY9hftu4Fql1IeiJ9xcby1PFZFdgUal1KciUq3KbFDHim2+xFpsBfwOHa43FFhpyZ3ItV5qRGQA8E+l1AEi0hM9o7Ireup8jKTM5dI+BgJ3osvHfSwij6GnyB9jbR8BNJRj+xCRS4BVSqlbrOXH0HnUD0UbBj2UUh+Vi9yiE6JdSTT9URYyF5NyGXhdgy6pdQk6cX9/dETHQACl1AfWD1RVDg+wxVK0krdDsxah44X3sJY/tB7geDk1KtHE0A/Ez8AZ6ERez6InE010KMuyUPAW09CDZfugc+h8j7Y2L0K7bMpGwVvMRSej+9ha/gPQVkT2hUSbLtewvbVAL0c0ynjgQ6zJQ0qHpJaT3KuJrj/KReaiUVIlLyKbiMjGOHI6Kx1SeAO6vuKRItJeRI4SkT5KR1GUFEvmvkqHtz2v9IBftbV5GrrBISL7i0gPVSYpSkWkv4j0AWqVHoSag57S3R54A23BL8GKNCgHZWnLLCKtLXlmogdSZyilbrcszceBsywLrRxk3kR0JaFVSqnp1rpqdOz4JGA7a50En6X4WHL3sRYfQGdlfEBEnkYXtzkJna1x8xKJGIZCP3tlrz9KQcmUvOh8zk+hrbJrRORUe5tSahL6h+oCPIm2mFuXQEwXDpnvE5EbbZkd1sFCYJnoqejXoZVoybHkfhr98P5NRH5lbZoA3Az8GzgNbbFdK7pwQknxkfkIdL7yJmBbK2IFdEbMJSUR0oOjfdyLo00rHYNdj27LZ4nIruXwQrJxyP2AiNyMnvW8L/p7PIbO+QJ6EH6DUsjoRUQOF5FnIGGQxK31Uq76o2SUYiAAbT1+ip5mviF6EPBp4BzPfteisx8OKoWcmcqMHoCdg3bhlFzmELmfQceVd0MPZh7r2L9Hmcr8P/T09Gr04No/0QpoHDCkTGX2a9NXo8dCSj6YHSL3s8DvPfv9Cm0U9CuFnB5ZtgWmAtPRidzs9Xb0jz3WWDb6o5R/pYqTb0Qrw7lKqfmiM78tAs4UkYVKqUetOOcuwIHKkQ2xhITJvEjpvCnz0d3yE5WjpFiJCZL7bGABsK1SqsEekFJKzSulsBZ+Mi9G+7UXAZeiLcpBwDSl1I8lkzRJ2jZt7fcOMEmVj+sgTO5l1rM4ED3WdJIqg1nDWLOblVJPi8jLIvKhUmpXpVSj5XdvFF1ToCvloz9KRkncNUqp1cB3aLdHO6UzAn4FPAcMFp0EawVwllLqq1LI6CWNzIOs3V4D9isjBR8m9zPoiAm7i1s2A1IBMn+JHvgbYe3zk1JqdJko+Cht2q6k9bZSam7pJHWTTm5rt2nAGUqpcSUR0oPSA9ZvWZ9HAktF5ENruVFEeiullqNlLgv9UUqKruQdA05/QXe1/2k1rhXAB8AO6OIazcqRurSURJFZRDZSSi1UVu6XciCC3NsDnZXVty0HIsg8DG2hlQ0R23TnEokXSMR23UMptU4ptbJUcvqhlPrZihJDKXUwWtG/LiKjgL+KSJtyMlxKSdGVvK1QlI46uRk9WPmqiGyGLurQGivhUbkQUeaS54L3ElHusoj+sTEyF4+IcpetolRWzV7r88HoGa83ATerEtcLKCeKNhnKihdvcn4Wkb5ot8wf0LGtfdADVeOKIlQaKlFmqEy5jczFoxLl9shs59HpCyxXSi0RXaz9HuCI9d0Hn0KuI7dhf+gZcrc4luOOz3ug83lsZm9Dx3CXdCS6EmWuVLmNzEbuHGTeHT0GNtBa3pkyiPwpx79C/kA7oFN6/gw85lhfDbRFh20dWeobUOkyV6rcRmYjdx5kPqLUclbCX8HcNSJyANBKKfWsiHwFTFZKneDY3klZgyeqTKbPV6LMUJlyG5mLRyXKnYnM6OGFsgkeKDcK6pMXkQ2Vjr0VYCw6pvk4a1sPVR4x2S4qUWaoTLmNzMWjEuWuRJnLkbwqeRHZA12ooZVS6jZrXY1Sap0V+/45Ot75dWA34ELlKVpRbCpRZqhMuY3MxaMS5a5EmSuCfPl9gIOAiejc2TOAOx3bqh2fl6NnL25VSj9VpcpcqXIbmY3cLU3mSvnL1w/UB/gY2Nta7oCeTJEoRmGt3wOdSbDkuSQqUeZKldvIbORuaTJX0l++ctesBa5SSr0tOoPhanTKWu9sylbAvqo8aipWosxQmXIbmYtHJcpdiTJXDDnNeBWd67sa+Fkp9QqA0lOgG9BdrmZrv52sba+W+geqRJkteSpObiNz8ahEuStR5kokayUvIiOBV9Alzh4WkS2s9XYu8g5AaxE5AXhERHrkKmyuVKLMUJlyG5mLRyXKXYkyVyyZ+ncAIVlUeQ+gO3AeutzZIMd+/wDeBN6jxD60SpS5UuU2Mhu5W5rMlf6X7Q8VR1db2YhkGOYf0XmpN7eWL0DX4dwiH4Lm/EUrUOZKldvIbORuaTJX8l9GcfIisinQCSvECfhCKXW9Y/uF6NzqvwGGAPNVifN9V6LMUJlyG5mLRyXKXYkytwgyePseDHyD7j7djk4eNAu42LFPX+CeUr+5KlnmSpXbyGzkbmkyt5S/SCGUIrILcCNwglLqKxG5G51AaBfgU2s22hPArsA2ItJZKVXS4sqVKDNUptxG5uJRiXJXoswtiohv4V2AUx3LXYGXrc/9gfvQ3a+xlMlMtEqUuVLlNjIbuVuazC3pL+qPFAfaOz73QteB7GGt2xioAjqU+gtVssyVKreR2cjd0mRuSX+R4uSVUk1KF8YFHQK1FFiilJonIicD/4fOL7EsyvmKQSXKDJUpt5G5eFSi3JUoc0si6yyUIvIAMA/YD90V+zaPchWESpQZKlNuI3PxqES5K1HmSiVjJS8igq7OMsn6f2+l1NQCyJY3KlFmqEy5jczFoxLlrkSZK51cLPlTgTGqgormVqLMUJlyG5mLRyXKXYkyVyq5KHlR2R5cIipRZqhMuY3MxaMS5a5EmSuVgpb/MxgMBkNpySnVsMFgMBjKG6PkDQaDoQVjlLzBYDC0YIySNxgMhhaMUfKG9RYRaRKRcSIyQUS+FpFzRSTm2ed5EfnE+ry/tf84EVkpIlOszw+JyB4issyxfZyI7FOab2YwJDHRNYb1FhFZqZRqa33uBjwGfKSUutxa1xFdwWglcJBSaqbj2HeB85VSY63lPazlg4v4FQyGtBhL3mAAlFI/AacDZ1mzMgGOAl5Ep8E9vlSyGQy5YJS8wWChlJqBfia6WatOAB63/k6IcIoRHnfNJgUS1WCITKSiIQbDeoQAiEh3YFPgQ6WUEpFGERmslBofcuwHxl1jKDeMJW8wWIhIf6AJ+Ak4Dl2PdKaIzEKXpjMuG0PFYZS8wQCISFfgLuB2K6fKCcABSqm+Sqm+wHYYJW+oQIy7xrA+00pExqFT3jYCDwM3iUhfoA/wqb2jUmqmiCwXkR2VUp8FnG+EdT6bq5RSzxREcoMhIiaE0mAwGFowxl1jMBgMLRij5A0Gg6EFY5S8wWAwtGCMkjcYDIYWjFHyBoPB0IIxSt5gMBhaMEbJGwwGQwvm/wEKAEyZwcfT7gAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "station_data.set_index('DATE')[['temperature_degC', 'dew_point_degC']].plot()\n",
    "plt.xticks(rotation=45);"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Open Weather API\n",
    "\n",
    "The NOAA data has historical weather for the weather stations, but it updates with a lag of a few days. I'd like to have more up to date weather data I can pull hourly. For that, I can use the [Open Weather API](https://openweathermap.org/).\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "from dotenv import load_dotenv\n",
    "import os\n",
    "import requests\n",
    "load_dotenv()\n",
    "OWM_KEY = os.getenv('OWM_KEY')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<Response [200]>"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "url = 'https://api.openweathermap.org/data/2.5/weather?'\n",
    "\n",
    "params = {'appid': OWM_KEY,\n",
    "          'lat': 39.847,\n",
    "          'lon': -104.656,\n",
    "          'units': 'metric'}\n",
    "\n",
    "r = requests.get(url, params)\n",
    "r"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'coord': {'lon': -104.656, 'lat': 39.847},\n",
       " 'weather': [{'id': 801,\n",
       "   'main': 'Clouds',\n",
       "   'description': 'few clouds',\n",
       "   'icon': '02n'}],\n",
       " 'base': 'stations',\n",
       " 'main': {'temp': 25.49,\n",
       "  'feels_like': 24.67,\n",
       "  'temp_min': 23.12,\n",
       "  'temp_max': 27.46,\n",
       "  'pressure': 1022,\n",
       "  'humidity': 22},\n",
       " 'visibility': 10000,\n",
       " 'wind': {'speed': 7.2, 'deg': 190},\n",
       " 'clouds': {'all': 20},\n",
       " 'dt': 1656473325,\n",
       " 'sys': {'type': 1,\n",
       "  'id': 3801,\n",
       "  'country': 'US',\n",
       "  'sunrise': 1656415956,\n",
       "  'sunset': 1656469859},\n",
       " 'timezone': -21600,\n",
       " 'id': 5412777,\n",
       " 'name': 'Barr Lake',\n",
       " 'cod': 200}"
      ]
     },
     "execution_count": 19,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r.json()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Index(['temp', 'feels_like', 'temp_min', 'temp_max', 'pressure', 'humidity',\n",
       "       'timestamp', 'lon', 'lat'],\n",
       "      dtype='object')"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = (pd.DataFrame(r.json()['main'], index=[0])\n",
    "    .assign(timestamp=pd.to_datetime(r.json()['dt'], unit='s').tz_localize('UTC'),\n",
    "           **r.json()['coord']\n",
    "           )\n",
    ")\n",
    "\n",
    "df.columns"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Timestamp('2022-06-29 03:28:45+0000', tz='UTC')"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pd.to_datetime(r.json()['dt'], unit='s').tz_localize('UTC')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "8ffe812cc978bbbceedabb2b9debbb786aca253c20b4f14ddd31576b2ba796e5"
  },
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}