{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.backends.backend_pdf import PdfPages\n",
    "import fiona\n",
    "from shapely.geometry import Point, shape\n",
    "from gisutils import df2shp, project\n",
    "from pydrograph import Nwis"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### instantiate an NWIS object using a polygon shapefile of the study area\n",
    "* `bbox.shp` is in UTM 83 zone 15 N\n",
    "* pydrograph will automatically reproject the extent to GCS NAD83 by default (EPSG: 4269), which is the typical GCS used by NWIS. Alternatively, the NAD27 datum can be used (EPSG: 4267) by specifying `datum=NAD27` to the `NWIS()` constructor."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "reading extent from ../data/bbox.shp...\n",
      "reprojecting extent from +init=epsg:26915 to +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/aleaf/Documents/GitHub/pydrograph/pydrograph/nwis.py:163: FionaDeprecationWarning: Collection.__next__() is buggy and will be removed in Fiona 2.0. Switch to `next(iter(collection))`.\n",
      "  g = shape(shp.next()['geometry'])\n",
      "/Users/aleaf/opt/anaconda3/envs/gis/lib/python3.8/site-packages/pyproj/crs/crs.py:280: FutureWarning: '+init=<authority>:<code>' syntax is deprecated. '<authority>:<code>' is the preferred initialization method. When making the change, be mindful of axis order changes: https://pyproj4.github.io/pyproj/stable/gotchas.html#axis-order-changes-in-proj-6\n",
      "  projstring = _prepare_from_string(projparams)\n"
     ]
    }
   ],
   "source": [
    "nwis = Nwis(extent='../data/bbox.shp')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### field measurements inventory table\n",
    "* fetches inventory tables of all misc field measurements within the polygon"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "getting site inventory for field_measurements...\n",
      "url: http://nwis.waterdata.usgs.gov/usa/nwis/measurements?nw_longitude_va=-92.700&nw_latitude_va=46.800&se_longitude_va=-92.600&se_latitude_va=46.700&coordinate_format=decimal_degrees&group_key=NONE&format=sitefile_output&sitefile_output_format=rdb&column_name=site_no&column_name=station_nm&column_name=site_tp_cd&column_name=dec_lat_va&column_name=dec_long_va&column_name=coord_meth_cd&column_name=coord_acy_cd&column_name=coord_datum_cd&column_name=dec_coord_datum_cd&column_name=district_cd&column_name=state_cd&column_name=county_cd&column_name=country_cd&column_name=land_net_ds&column_name=map_nm&column_name=map_scale_fc&column_name=alt_va&column_name=alt_meth_cd&column_name=alt_acy_va&column_name=alt_datum_cd&column_name=huc_cd&column_name=basin_cd&column_name=topo_cd&column_name=inventory_dt&column_name=drain_area_va&column_name=contrib_drain_area_va&column_name=tz_cd&column_name=local_time_fg&column_name=reliability_cd&column_name=project_no&column_name=rt_bol&column_name=peak_begin_date&column_name=peak_end_date&column_name=peak_count_nu&column_name=qw_begin_date&column_name=qw_end_date&column_name=qw_count_nu&column_name=sv_begin_date&column_name=sv_end_date&column_name=sv_count_nu&date_format=YYYY-MM-DD&rdb_compression=file&list_of_search_criteria=lat_long_bounding_box\n",
      "reading data with pandas...\n",
      "finished in 1.27s\n",
      "\n",
      "culling 1 sites to those within extent...\n",
      "finished inventory in 1.31s\n",
      "\n"
     ]
    }
   ],
   "source": [
    "field_sites = nwis.get_siteinfo('field_measurements')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### daily values sites inventory table"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "getting site inventory for daily_values...\n",
      "url: http://nwis.waterdata.usgs.gov/usa/nwis/dv?referred_module=sw&site_tp_cd=ST&nw_longitude_va=-92.700&nw_latitude_va=46.800&se_longitude_va=-92.600&se_latitude_va=46.700&coordinate_format=decimal_degrees&group_key=NONE&format=sitefile_output&sitefile_output_format=rdb&column_name=site_no&column_name=station_nm&column_name=site_tp_cd&column_name=dec_lat_va&column_name=dec_long_va&column_name=coord_meth_cd&column_name=coord_acy_cd&column_name=coord_datum_cd&column_name=dec_coord_datum_cd&column_name=district_cd&column_name=state_cd&column_name=county_cd&column_name=country_cd&column_name=land_net_ds&column_name=map_nm&column_name=map_scale_fc&column_name=alt_va&column_name=alt_meth_cd&column_name=alt_acy_va&column_name=alt_datum_cd&column_name=huc_cd&column_name=basin_cd&column_name=topo_cd&column_name=inventory_dt&column_name=drain_area_va&column_name=contrib_drain_area_va&column_name=tz_cd&column_name=local_time_fg&column_name=reliability_cd&column_name=project_no&column_name=rt_bol&column_name=peak_begin_date&column_name=peak_end_date&column_name=peak_count_nu&column_name=qw_begin_date&column_name=qw_end_date&column_name=qw_count_nu&column_name=sv_begin_date&column_name=sv_end_date&column_name=sv_count_nu&range_selection=days&period=365&begin_date=1880-01-01&end_date=2020-09-22&date_format=YYYY-MM-DD&rdb_compression=file&list_of_search_criteria=lat_long_bounding_box\n",
      "reading data with pandas...\n",
      "finished in 1.27s\n",
      "\n",
      "culling 1 sites to those within extent...\n",
      "finished inventory in 1.32s\n",
      "\n"
     ]
    }
   ],
   "source": [
    "dv_sites = nwis.get_siteinfo('daily_values')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### write shapefiles of the inventory tables\n",
    "* `shapely Points` were created from the lat/lon coordinates when the tables were assembled; these are used to write the feature geometries for the shapefiles"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "writing dv_sites.shp...\n",
      "writing misc_sites.shp...\n"
     ]
    }
   ],
   "source": [
    "df2shp(dv_sites, 'dv_sites.shp', epsg=4269)\n",
    "df2shp(field_sites, 'misc_sites.shp', epsg=4269)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### get daily values for a list of sites"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['04021520']"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sites = dv_sites.site_no.tolist()[0:2]\n",
    "sites"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "http://waterservices.usgs.gov/nwis/dv/?format=rdb&sites=04021520&startDT=1990-01-01&parameterCd=00060\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "dict_keys(['04021520'])"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dvs = nwis.get_all_dvs(sites, start_date='1990-01-01')\n",
    "dvs.keys()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "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>agency_cd</th>\n",
       "      <th>site_no</th>\n",
       "      <th>datetime</th>\n",
       "      <th>72604_00060_00003</th>\n",
       "      <th>72604_00060_00003_cd</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>datetime</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2005-05-26</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>2005-05-26</td>\n",
       "      <td>199</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-27</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>2005-05-27</td>\n",
       "      <td>200</td>\n",
       "      <td>A:e</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-28</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>2005-05-28</td>\n",
       "      <td>200</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-29</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>2005-05-29</td>\n",
       "      <td>187</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-30</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>2005-05-30</td>\n",
       "      <td>169</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           agency_cd   site_no    datetime 72604_00060_00003  \\\n",
       "datetime                                                       \n",
       "2005-05-26      USGS  04021520  2005-05-26               199   \n",
       "2005-05-27      USGS  04021520  2005-05-27               200   \n",
       "2005-05-28      USGS  04021520  2005-05-28               200   \n",
       "2005-05-29      USGS  04021520  2005-05-29               187   \n",
       "2005-05-30      USGS  04021520  2005-05-30               169   \n",
       "\n",
       "           72604_00060_00003_cd  \n",
       "datetime                         \n",
       "2005-05-26                    A  \n",
       "2005-05-27                  A:e  \n",
       "2005-05-28                    A  \n",
       "2005-05-29                    A  \n",
       "2005-05-30                    A  "
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dvs['04021520'].head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD8CAYAAACb4nSYAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAABDAUlEQVR4nO2dd5gb1dX/P0dbve5lXbANptjYlNAcAoQaTDDhDSUJid8k/ICXhFASIKQAKW9IITiNVCDhJbRAANMJYEwzHYw7xtjGva69677r9Vbd3x8zI41GI61WM7vSSufzPPusNBrNHE35zrnnnnuuGGNQFEVRioNIrg1QFEVRug8VfUVRlCJCRV9RFKWIUNFXFEUpIlT0FUVRiggVfUVRlCKiNNcGdMSQIUPMmDFjcm2GoihKj2HIkCHMmDFjhjFmsvezvBf9MWPGMGfOnFyboSiK0qMQkSF+yzW8oyiKUkSo6CuKohQRKvqKoihFhIq+oihKEaGiryiKUkSo6CuKTUNzG+u3N+baDEXpUlT0FcVmyp3vctJvZ+baDEXpUlT0FcXmw427c22ConQ5KvqKoihFhIq+oihKEaGirygedApRpZBR0VcUD6r5SiGjoq8oHlTzlUJGRV9RPGh4RylkVPQVxYNKvlLIqOgrioeoevpKAaOirygeVPOVQkZFX1EUpYhQ0VcUD+rpK4WMir6ieDDalasUMCr6iuJBPX2lkFHRVxQPqvlKIdOh6IvI3SJSKyIfupYNEpGXRGS5/X+g67MbRWSFiCwTkTNdy48RkUX2Z38REQn/5yhKcDRlUylkMvH07wUme5bdALxijBkLvGK/R0QOAaYAh9rfuV1ESuzv3AFcBoy1/7zbVJS8QDVfKWQ6FH1jzBvAds/ic4H77Nf3Aee5lj9sjGk2xqwGVgDHisgIoJ8x5l1jjXG/3/UdRckvVPSVAibbmP4wY0wNgP1/qL18JLDetd4Ge9lI+7V3uaLkHZq9oxQyYXfk+sXpTZrl/hsRuUxE5ojInLq6utCMU5RM0PCOUshkK/pb7JAN9v9ae/kGYLRrvVHAJnv5KJ/lvhhj7jTGTDTGTKyurs7SREXJDtV8pZDJVvSfAS6yX18EPO1aPkVEKkRkf6wO2/ftEFC9iBxnZ+38P9d3FCWv0OwdpZAp7WgFEXkIOBUYIiIbgJ8BU4FpInIpsA64AMAYs1hEpgEfAW3AVcaYdntTV2BlAvUCptt/ipJ3qOYrhUyHom+M+e8UH52eYv2bgZt9ls8BDuuUdYrSjYhYgq8duUohoyNyFcUmlm2gmq8UMCr6imLjDBJXzVcKGRV9RfGgMX2lkFHRVxSbiB3f0ewdpZBR0VcUG0HDO0rho6KvKA62p2/U01cKGBV9RbFxsndU85VCRkVfURSliFDRVxQbiYV3cmuHonQlKvqKYhOJ5emr6iuFi4q+otg4Mf2oar5SwKjoK4pNbESuxneUAkZFX1E8qOQrhYyKvqLYaMqmUgyo6CuKg5bZVIoAFX1FsYll76jmKwWMir6i2Eis4Fpu7VCUrkRFX1E8aJ6+Usio6CuKjXbkKsWAir6i2IjG9JUiQEVfUWxinr6Gd5QCRkVfUWy04JpSDKjoK0oM6XgVRenhBBJ9EfmuiCwWkQ9F5CERqRSRQSLykogst/8PdK1/o4isEJFlInJmcPMVJTxE58hVioCsRV9ERgJXAxONMYcBJcAU4AbgFWPMWOAV+z0icoj9+aHAZOB2ESkJZr6ihIdm7yjFQNDwTinQS0RKgSpgE3AucJ/9+X3Aefbrc4GHjTHNxpjVwArg2ID7V5TQiMX0c2uGonQpWYu+MWYj8HtgHVAD7DLGvAgMM8bU2OvUAEPtr4wE1rs2scFeloSIXCYic0RkTl1dXbYmKkqnELS0slL4BAnvDMTy3vcH9gF6i8jX033FZ5nv3WWMudMYM9EYM7G6ujpbExUlK1TylUImSHhnErDaGFNnjGkFngBOALaIyAgA+3+tvf4GYLTr+6OwwkGKkhdENGVTKQKCiP464DgRqRJrKOPpwBLgGeAie52LgKft188AU0SkQkT2B8YC7wfYv6KEis6cpRQDpdl+0RgzS0QeA+YBbcB84E6gDzBNRC7FejBcYK+/WESmAR/Z619ljGkPaL+ihI5KvlLIZC36AMaYnwE/8yxuxvL6/da/Gbg5yD4VpavQEblKMaAjchXFg4Z3lEJGRV9RbDRPXykGVPQVxSaep59jQxSlC1HRVxSbWMqm+vpKAaOiryg2OomKUgyo6CuKjRZcU4oBFX1F8aDhHaWQUdFXFAfN01eKABV9RbGJz5GrKIWLir6i2ES09o5SBKjoK4qNlmFQigEVfUXxoB25SiGjoq8oNjoiVykGVPQVxUbDO0oxoKKvKB5U85VCRkVfUWx05iylGFDRVxQbJ08/qpqvFDAq+opi48T0NcCjFDIq+opiox25SjGgoq8oNrGUzRzboShdiYq+otiop68UAyr6iuJBR+QqhYyKvqLY6MxZSjEQSPRFZICIPCYiS0VkiYgcLyKDROQlEVlu/x/oWv9GEVkhIstE5Mzg5itKeMRTNlX1lcIlqKf/Z+AFY8x44AhgCXAD8IoxZizwiv0eETkEmAIcCkwGbheRkoD7V5TQiKdsKkrhkrXoi0g/4GTgnwDGmBZjzE7gXOA+e7X7gPPs1+cCDxtjmo0xq4EVwLHZ7l9RwkbnyFWKgSCe/gFAHXCPiMwXkbtEpDcwzBhTA2D/H2qvPxJY7/r+BntZEiJymYjMEZE5dXV1AUxUlM6jHblKIRNE9EuBo4E7jDFHAXuwQzkp8Gs8+95dxpg7jTETjTETq6urA5ioKJmjHblKMRBE9DcAG4wxs+z3j2E9BLaIyAgA+3+ta/3Rru+PAjYF2L+ihEpE8/SVIiBr0TfGbAbWi8jB9qLTgY+AZ4CL7GUXAU/br58BpohIhYjsD4wF3s92/4oSNs6IXM3eUQqZ0oDf/w7woIiUA6uAS7AeJNNE5FJgHXABgDFmsYhMw3owtAFXGWPaA+5fUcLD8fRza4WidCmBRN8YswCY6PPR6SnWvxm4Ocg+FaXLUdVXChgdkasoNrGUTVV9pYBR0VcUmzAKrhlj2L6nJRyDFKULUNFXFJswSis/Mns9R//yJZZu3h2OUYoSMir6imIThqf/xnJrMOGK2oYQLFKU8FHRVxQbR/SDpGzGWgvaLaDkKSr6imITysxZmvap5Dkq+opiIyFUXItvQmVfyU9U9BXFQxC5Fq3PrOQ5KvqKYhNGwTUtz6zkOyr6imKjoRmlGFDRVxSgdncTu/a2AhAN4unHOnLDeXCs29bIddMW0NoeDWV7iqKiryjAsb9+hQXrdwIBY/r2/7AaC9c//gFPzNvI7NXbw9lgSLy/eruOReihBK2yqSgFR5DwTtgTsURstyxI66Mr+PI/3gVgzdSzc2yJ0lnU01eUEIkXbQtre9YW27WfQQkJFX1F8RBIX2OlHMIRac0AVcJGRV9RPGhpZaWQUdFXFA/B8vRDKOXgw0V368yiSjio6Pdw1m7bw71vr9bc8hAJI2UzLNVvz7ceXKXHo9k7PZyr/j2PDzfu5r+O2IchfSpybU5BECS8E/bsW++s3BbKdhTFQT39Hs6W3c2AeoRhEii8E0JN/p7ECx9u5qn5G3NthtIJ1NPv4TieZZAa8PlAU2s7dfXNjB5UlWtTAiEUV7rN5Q/MBeC8o0bm2BIlU9TT7+FEbNeyrb1ni/510xZw0m9n5kW5gTD6R3r22cgtD7y3lsN/NqPg+6kuvXc23390YbfvV0W/h+OEE3p6eGfmUmuawXSi//D76xhzw3M0t7V3qS0a3sktP336Q+qb2/JuFHLYvLK0lsfmbuj2/QYWfREpEZH5IvKs/X6QiLwkIsvt/wNd694oIitEZJmInBl034rL0+/hd0hpifU7WtO0WP7yynIAau1+jK4inwquFSOFErLMV8Lw9K8Blrje3wC8YowZC7xiv0dEDgGmAIcCk4HbRaQkhP0XNYXi6ZdGnDBVak8/Eume+WeDCbZlYw8/HTnFcWR6+jWdrwQSfREZBZwN3OVafC5wn/36PuA81/KHjTHNxpjVwArg2CD7V9yefu5j4UFwfke6+7y7vOggD5WSiLMNFaxsKYYQ2eZdTTnbd1BP/0/ADwG34gwzxtQA2P+H2stHAutd622wlyUhIpeJyBwRmVNXVxfQxMKmJFIYXpHjxadr0sdGu3a5p589JeqlBsYJ8RVykTknnJkLshZ9EfkvoNYYMzfTr/gs8z2rxpg7jTETjTETq6urszWxKHC8op4e089ELO3nQtfHekMordzDT0desGV37rzhriaXz7MgefqfBs4Rkc8BlUA/EXkA2CIiI4wxNSIyAqi1198AjHZ9fxSwKcD+FQon/hnJoG8iLqhdHN4J8F2n5RXt4ecjl+w3uIq12xrzIn23q8hl+C9rT98Yc6MxZpQxZgxWB+2rxpivA88AF9mrXQQ8bb9+BpgiIhUisj8wFtAqUgFxxLKn5+lHMghTxTutw933zsaWhPdBHiqxcFsBhya6mkJxZNKRy1/WFSNypwLTRORSYB1wAYAxZrGITAM+AtqAq4wxXZtwXQQUyg2SiVg68cGwO629uwyi14VyPnJJLIxXuI5+TtNRQxF9Y8xrwGv2623A6SnWuxm4OYx9KhZSaNk7mYR3Qv6pEc9MJcHCO/Y21NPPGscB6OnXdDpy6RPoiNweTiax8J5AJh3Ssd8asqCK5y4Iw9Pv4acjp3TUWrr24flMX1TTnSaFTi77fFT0eziFMiK3JINO2rgYhOsBetPKApVWtm18c3ldj38Q54qORP+pBZu44sF53WlSQaGi38MJy9NfvGkXh980g9ocpclFOhG6Cbsj13vkgnhhzsNr9pod3D5zRQCripeIrUphPzRvm7mCmctqO16xG8hlTF9Fv4cjIXn6/3p3LfVNbby0ZEsYZnUaySB0012dpEEeKiWuO2r11j2BbelTUXzVz2PnOWRh/N2MZVxyz+xQt5ktGtNXsibu6Qdzf8tstcpV6mcmgt5VdYa82hJW+CgMK92/NZ87h8OMURdKyDId6ukrWRNWPf1cl3MoyaQMQxd15HrVOYjYuE0LQ6TdxyOf+wjCPCexazEEB2T6ohqWbt4deDth01NH5Cp5QFghD0dQc+WBxPafNnunazpyvQQ5libF62xxn4+2qKE0T+vStkcNZSHZFmaWltPhu2bq2YG3FSY9ckSukh+EVXvH6YDM1bUoGcRxnSyb8DtyE/eZrqZ/h9sy/q+zpT1qKLdDb/lcXz5M2/J5gJsxhm/cN4c3Pg5WCFJj+krWhHWDZFLlsispiXn6aVbqIjEIM6bvfoC8ubyOs/78ZtY1ZIwxRA2UleSvCDqEaVs+x/Tbo4aXl2zh4nuCVZDRmL6SNU78M2hxqnh4J6hF2ZFJxkZXDUTzbi2smP6OxlaW1Oymtj67mb4cM8pKbU8/jweohmmbk7KZj0XrHIuCmqair2RNWBktkQwGR3UlGZVhsP93dTGzsGL6DuUl2d1mjh1OZlU+F3EL07bu9vTnrdvBtDnrO16R8O4P7chVsiasG8TxonPVwZRJR7Jk8GDIBu9vDnQsfez3lPbJGOdYlPWAiXK6IrzT1R32Dl+4/R0AvjxxdNr1fvvCUlraQkrnVdFXsiWskEf8RgtqUdD9dxzeCdsD9G4tyLH0+2q2D6mY6Jfmf0dumM5CPH04tE2Gwu2vrQxtWxreUQIT3NPPcXjHieNmMF1iV8d6A8X0fQI82W7OG945+y9vZW1XVxNWeGdvS7srSyvPVD9E3Nd5d7euVfQzoLU9yq7G1lyb4YtzuQS9QSKxlM1cj8hNs5KtBj98/INQz4f3J7cFcDH9Dl+2D1LnlFbYnv7Whuw6hLuDMMI7LW1RJvzvC7yytDa0bXaGzjoTb6/YmvW+3Hvq7t+pop8B1zw8nyN+8WKuzfDFuV6Cevo9KXsHYMbizaHt2+udB/P0k8lW9J1jURnWqKcuJAyn3JuB1t0pmy2dfNiv3daY9b7czlV333Ma00/DY3M3sLOxhecXWQJjjIl1JuYLzsUTdMh6rCM3RxO5ZTIiV1xFkFu7sOkfKHvHN6YfzA7H089nuiKzqLs94Jb2aLc9YN0/rbt/p4p+Gr7/6MKE9y3tUSrybBy840UG9/RzO/lHJn0KEZf2hVoYzhveCT2mH6wjN99E3y8EGIZweY9Td6eodjYzJ5LG/zv0f1/gu2eM4xsnHeD7udu56e7fmV9XU54TVrpWmDheZE+vvZNZ9o7L0w8xtcO7y0D9I2HG9GOiH66jsXlXE88s3JT198Pst0jcRuL7MAqudYbOXlPpGv17Wtr51XNLUn6eENPv5t+pnn4GlEaEtqjJT9EPydOPCWrOPH3rf6biEWaT2LvPIK0I/5h+dttyfmNlWbi+2dfueo+VdXs4Y8IwepV3/oHi93NC8fQ92/C7prNNNMjEvs7e35I051rmJBbS615dUU8/A+IZMvmXJ+1cO8Gzd6z/uff0M1s/zE6+pLBCoJh+iOGdWPZOuJ7+xp17gez7b7orvON33LK9PNMJujMuoNNOXQrNz+TB5F5Fwzt5gnNjQPyCzscCUGF7+jmL6TtT5GUwIhfCDe8kpWyG3ZEbOHsn3Ns0GnMUshR9322GH94J835LJ+hOQbvOZu+k8vMzMTuX8yRkfTWJyGgRmSkiS0RksYhcYy8fJCIvichy+/9A13duFJEVIrJMRM4M4wd0Fa/5zKWZj56+c/EEj+nndnCWZDBOwP1ZmB25SeGdQFU2fbYfNHvHlVESyjiKgKLvd42EcW94f5vfNrPdS3N7e8rPnMFvnQ7viDBzWW1SWCqTc+T+SnfPVhfEhWgDvmeMmQAcB1wlIocANwCvGGPGAq/Y77E/mwIcCkwGbheR/EqFcTFzaXK97Hz09B2Lwqu9E8ye7PffcUeu27YwUzbD7EAMdXCWT/ZOmNdgttvqqo5cbyvPV/Sz3E86QS/PUvR//p/FXHLPbO59Z03C8kwOq+mJnr4xpsYYM89+XQ8sAUYC5wL32avdB5xnvz4XeNgY02yMWQ2sAI7Ndv8Z2Jf1d5vb2nln5Vb69ypLWJ6Pw8JjTfWA3oLTVM1dTN/6n+4GiHabp5/9tsOMQ/uKfgi/24nlh+nph6FbyVlU4Z3jdBPjxDz9ToZ36pvaAFi7bU/C8kz6Snp8TF9ExgBHAbOAYcaYGrAeDMBQe7WRgLt+6QZ7WehEo4Zzb3ub66YtyOr7s1fvoLGlncmHDk9Yno+e/sL1O4EQPP1I986cFY0a7n93DU2tVrO7JIPwklsEwhSETMIKQcj2po5n78QbxGG0cBxzsr1m/L4WxjGbu3ZHh9vMdi9pY/qlweak8J7fTE53j4zpO4hIH+Bx4FpjTLoZiP36PXx/rYhcJiJzRGROXV3npyWLGsMHG3bxxLyNnf4uwMxltZSXRjht/NCE5WF4Wbe+uCzrh1E6grZCnJj6v95by+qtezpYu2Oue2QBP//P4pSfP71wI//79GL+9uoKwD0ZTGaeflfm6QfryO2K7J34bbp7b3g1hzJpHRpjuOmZxSyp2Z2wzEvQInit7VGufmh+wjL/lM3stp++Ize78I6D91LM5Hz31Jg+IlKGJfgPGmOesBdvEZER9ucjAKdHdAPgLlg9CvAdIWKMudMYM9EYM7G6urrTdpWWRJg0YSiHjOjX6e+CJfrHHTCYQb3LE5aH8UT+y6srsn4YpSPwiFzX6y/c/nYwY4An5m/knrfXpPx8T7Pl4W9vbAGscwbpbzz38Q9X9MPz9J1vnv2JEfFlAbN33B25YVTajPcDpT+GyzbX87W7ZnHvO2u48J+zYst9Pf0ATcQNOxr5yyvLk7cZ0JFxH/eWdB25dupYc5ain9yR2znbeoynL5Zr+E9giTHmVtdHzwAX2a8vAp52LZ8iIhUisj8wFgg20WQaSiORrA7m2m17WFW3h9MOrk5KlevuQRSZUGUPrmluDXiDuF7v6IaKoiWRxGwdJxsznZi7HaKwwlDb97TQ2hZi9o6BgVVl/PycQ2PLgg7Ocnv6u0Lw9J2HXEce6Q8eW8g7K7cBiYLo9xAL8hD+zkPz+avd4nPjt8nOjC1wm5nuuDnhnWw9fa/DtX5Hx4XYEmrvdHNMP8iI3E8DFwKLRGSBvexHwFRgmohcCqwDLgAwxiwWkWnAR1iZP1cZY1I/fgNSUiJZxT9fW2aFk047eGhSx053N8MywRGGhua2YBvq7s4kz0Tozu7TiYfbowqj7l171HD0L19i/PC+CcuDjci1ivKVu4Q6aIdp2IOzMo3pu6d5dAui39eCjFbf2+IvA36efmcuU/dD7X/unZNyPSe8ky60mA6vnZ//a8etsURPv3udyaxF3xjzFqnHJ5ye4js3Azdnu8/OUBYRVtXt4ZqH5/PnKUdl/L2Zy2rZf0hvxgzpzTpP6dR8zNN3bNrTElD0XQysKut4pYBIilLK6W68xPBO8HPRaB+zpZvrE5a3RU3WFVWNsW4Kt2BmLfqxPP2uGUO5tKae8cNTh0BLXBXFmhNEP/n3ZBsaAVIe56Ahy0y/HY/pp/ZB082x4LUzk2vT/ZUwruXOULAjckvsON3TCzIvLNXU2s67K7dx6sFWP0JyeCcPRd++ARuagom++5cNrCpPuV5YJGfrWP+9Oc9u3GLT3NbO3LXbueKBuVmLajqPPttTbbBaIWUu0b8jy2n2YjH9Lqqyee0jC9J+XpZiQnc/bzuI6KeqVtlRi6ujvpJUH9fWNyW8z2REbn2a+yubjnp3mMrJYOsuClb0nRPZGd5dtY3mtiinHmxl7VR6ilEF9fQ/3hL3KIPMzuQQjZrYhd3Q3MaVD87lwVlrs9qW+7rtjkdbbHpE+5hmct+4j39Ta5QrHpjH9A83J93EmZLuZs02rm9tUhK85Nlrtme1LceE0khubtOSFGrsJ7aZhnfO/OMbCZ3CkDpU53cO3LvuKK8+Vfz/+sc+SHjvPNymf5h6Yp6nF8STL47ed0CinRm2Tt24F6d7oHQFBSv6qS5YgPdWbeOPL32ctPy1pbVUlkX41P6DAKj0xFKDevqf/eMbsddNIVTsdDzBshKhuS3K84s28+MnP8xqW+4bOcwpE1NtKz5TlrOe9f/YMYNSbssdCmpqbac0NoG2yeohmu50Zv+AN0kiFjT/O4XD3eWUpriH0sX0N+7cy7/eXZNym8u21PPm8sRpBiMpVL/FR0zdQt5R6yLVZbzX41k7v3P+up0pt7Xv4CoArjz1QJ648tMJn6VzHlL1tbnvi8D9cZ2kYEU/VdMUYMqd7/FnT4qYMYaZy+r49IFDYoNhykokoekZZodLGGWaHWHqXRG8Qrb7sl2zrZH6puRsB79lHZHKG3Meyt50t+YOOnLPOWIfTjxoCE2t7bFY8Lx1Oznox9N5a3nn5ixN93ALUqLAq5XfOuXArLblHJtUotjVlKZoLTsi18uVSuqc5ysfnMdPn17MJlfBwo5I9etaPffIMb98iX++uTr2vqOwSKrTm86nSTXeoE+F1c/16YOGJH2W7lpJJejuB4V7DER3ULCin87T92PV1j2s297Iqa4BWSKSOIiim6r+ZYpz4fQJQ/Q9P+0/C2sS3n+8pZ7Db3qRx+Zu6NR2d+/1v+hLPLV2HA9ub5oO6agtqJVlEZrborFz/N4qK63w9Y+Ti+SlI12qXLZlLaLGJNVZ97YYM8U5Np29lsOiNIXjFBN9V/jTuZ7r7dRIrzedjtQduYn3yLY9LfzB1ULvKE05VXhn1up4uO3ON1ayeFNcdFOFWhxb/Fo/6VqFqfra3Jfe/e9mF5LNloIV/br65g7XcT/VZy61BOPUcakHg4WZvROG6DsPoVBE3/O+3NN5uMzOcHGOU6Zcet9s3+VO2QdHeJ2boDFF+h5Yxz8SEXqVl7Knuc1VJC7zfgE36U5ntuUOjEmOUXdGAN044pozTz9lTN/637cyft0517PzgOpM2muqn+cX3nHTnCbbBtJfD8418+vnl1KzK94ntDtFa9b5PX6tn3QPn1St41zmhBSs6GfSEecOPby2rI6xQ/swelBV6m2GmFqVboQgWHHgjh4M0TBF33OH9PZ0YmerOx9s2OW73PH0X/poC1sb4g/oVDnbYIlgiQj9e5WyZlsjm3dbN2u2Ubd0pQOC1Jr3HqpsszO6WvSr+1ak/TxVB7Jzqax1pTQ795LTOuhMR3iq3+cN73hp6tDTT02q/oBUg7ic68HvmKRLl67PILzT3RSs6LtPzg88E5w7OBfqnuY23l+9PanWjpcwR+R2dMGe9vvXGPeT6WnXcS7EPpWZif47K7ayora+4xVJ9ridkEW2sy15cd87t81cEdtqR55+SUSYvXoHED+GsdZCJ23w3nePXn48v/vSJ4DsH/CWp58oYtmKvuOThB3e2X9IbwCOO2Bw2vVSd+Rax8bd6e44KO7O9UxJmbLput/8+l82dtBvkK7PJpV3nipC4HTG+52LPWk6YlOHdxJt6860zYIVffchfTRFHNq5sd9ZuY2W9mgsP9/NKa5wT6gx/Q4yOjbs6LgjrD1DT//pBRsZc8NzfPWuWUy69Q3fdbz3R6PnIoxNoRvSIXDHvdtdqad7W9tT3qxRY4V3+nsGj3Um7dONN6YfkXjzPXtPP/l72YZ33Nk7F58wJqtt+G7X/m3pBiO98GENj8xZ7/uZI/pTjo2X0nJCLbHwTgfHz51tlSqm7x605Hc+vvWvuWmFPZ0J3tCQMyZn3Xb/EgrO/v0SRBqaUx/HVH0Ejm3nHLEPEG4hvY4oXNHPQAGcC2/mslp6l5cwcb/kdMHxI+JD9MOM6X/h9nd44L1gHTiOKPTtwNN/8L11nd52Y8pUs05vypfEiaG9Xo//A7E9aoV3BnjmOci0lkw6G8ASH2dQX9YPeJ+YftbhHVf2zo2fG5+dPX7bNY7op3Y8Ln9gXprvW//9RuyWZhjTdzsV77s6Vh+/4oTY645KPwBsbWhJvZM0JjS1RhM0Ysxgq/WTKtumNU2nutvT94YMG5r9xdzZ9eEj+wPh1FTKlAIW/Y7XaWm3TvxrS2s5ceyQpM5LSExLC7v2zk+eyi6n3iGWslmeXvRP8WnBePF6qN4wS2cmTu/sdHEtbdGE/TemiJE64R1vy0ak8975sx9s4vzbEquJRkRiojXp1tcz3pYbg19HbrYle+NCU1FawrlHWl5h0HEUztfTtTZThXbsLQDW8fqR/TBK7shN/5sbbe94vcez7t+rlAOqLQF2h3dSXXfpQjzONeX3W5rb2hM0orKshH6VpdTu9h/o154me8fdknMcsevOGIdI6vCOs95Au5Jvqg7krqBgRT+dOLkzDD7e0sCmXU2cdrB/PN8t+vlWe8e5JzqK6WfSEeg+XL3KSpKEN5JUNiE1mRynhLK3bdEEryxVXN9K2ZSkJrbz6/7ViZbTNQ8vYHeT9zcGj58bn5TNsDpyxw2zWp1BSh6AO7yTZgrBNKUfnNMbEeGykw9k/PC+tLRF+dury2MToXQUvnQ8au/1VBqJ8Or3TuXrx+2bEN5JddltSFPR0q+wnNMqbmqNJrg5G3fupbpvBXUNqWL6idk7V54aH3vhDOa01rN+d1lJhD4VpUnXmIPTInDqXKVKbe4KilP0JT5Tzkx7AvRTU4i+e9aifKu943hCHcX0M0kPdX7Z2zd8hqrykiThLelEB10mxynZ04+TKgZuefrJaXPpBnSlwq9MR0SEvgEzofw8/cCibx975zyn6zjszHb9UiLfXbmNcT+Z7vvgdc59/GFkLa8ojdDUFuX3L34cO/cdFRFznArvA7LMftiUlUQSsndS3c/pUrP9vuE8QJva2hO2WVffbIl+iu15s3e+cHR80j+3U+WEJnuVRehbUZoyXORsz/H0NbwTAj/63ISU1SLdszTNXFrLhBH9GN6/0nddd1pavs2R61y0HcX0mzydVn4PAef6H1RVTlVFGtG316vZtZfP/vF1Fm9KTsnsbAjohcWbEy76VJ5+u92R6/X0vTH+TPB7UIqQ0EncUR64H06VTTd+aahfuuMdLrs/dblfcGXvSKLoBx227zxw/a6DO99YmdJJ2L6nxbbL2oATVuvXqyxJtOas3c767Y18477ZvuE6ZxId73gI52FcXhJJ+CzVNZXuAeg/l6+1rLk1mtB6OGL0AIb2raQ2hei3xVJSLfv6VsavE7eT4jzgK8tK6FtZljK849jhFDfU8E4I7De4Nw9+4zjfzxwBe/GjzcxZu4PT0sS803U45prpi6wCUe7BJX4s8uTK+3kVTvxTBKrKSpNu1Linb138s1Zt5+MtDZz9l7f4z8LESqbeQWN19c0s35KYKuo9lK+6Bn2lvFGihohIUmptupIbqfALkUREEsZp7EmTlZEKy9O3jtXXPrUv4N9ymbN2By9+tCXtttpjnr713gnjBS3Q5VzTfmKcrqSHU9jO6dsqtycfGVhVzq7GxA7Vf7y+ilumL+HlJbUJ59bBEWtvy7HK7p8qK4lYLcBYJ318ncqyCA9c+ikg/bFwbl13iPa6M8YBUNfQlNCPdPS+A9J6+s45dLbldrTc95Nb9LfUN/HCYv8ibs45cByWXd0wcZFDwYo+pK5D7lxof3p5Oe1RkzK0A3DoPvF64/kW03/Srvy3orYhYbm7E62tPcpbKxJr0viKvuunLdtSz4zFWxI8XcfbdG54t7f5Hc/cpk4Jg4G9rQv67L+8yRl/TEwVTVcCIdVYAmdw1injqrnxrHg2izt+nGknp59YREToV1nGtZPGAtmFUaw6/Nbrm88/nAuP2y9rkW534sO26vcNEN55ZuGmWI0XR3DWbmtMGBgH6UOFjiA6cWsn1DGgqsx3trX45CTJD1hnQJP3M0dUe1eUEjVxsXWf12gUThxr1cD5xxurUtrriPr1kw/m4hPG8MFNn2XKJ60Hcc2upoRr/urPjKW6bwWNLe2+LamG5nbLIbIHLbofJE4LCOLhncqyCDvtY/Leqm1856H5CZ3Wzs+uLCuhsiyinn5YpKpD7vW8Rg7slXIb59gZExBssoO5a/3L63q9cMj84eIIsNfTdXuxfgLvt8z5jjseXbs7Lgji6chNF4N0BN0pUlXrEQtIPxr2Lz5T5xljrI5cu8XR3xXScYtgusFdHeHEqA+2O0yzCaN4wzsDq8rY3dSa8pymi/d7Ow97BwjvXP3QfM7685tEoybh2G/2tBLTib6zboursxJgQFV50vUwvF9lXPTbkn+704ryZsSVeM6v08GZ0AfkeVCkDMPZ3yktiXDTOYfSr7KMXuUlVJRaguyI/g1njWdg73KG2iOU/bz9hqY2epeXxu4D99iC3U2tMUfrc395E7DmNX7t+6cC8Ovnl/CfhZu4/vF4SeeoqxXX3yc81pUUuOj7F7o6wB6R6FCWJmPDvY0gMX13USc3G3cmZx/8+vklsdfpPNd4poCnY9Ml+jt9LqYde5Jzm52Jqd0da24P1fGanNDNzsaWhNmh3DgCN9gzsbx7AIqfCI6xy9ee63rQOsRywyVRAAG2uXK1t6XL2+4A50Z2jt+Li9OHX/yIGpOQATSwdznG+IcAIPHB6sUrrk54J0hM/7CbZrC7qY1xw/oAya2GdJlgzgxj3vCOX59Kdd+KmFfsH9O3lqUKmTrhE6eInt9oeGfy+VTn3Nmyt2N9UO9ytu9piQmv83F1GtHf09xG7wp/PfGeX7BaZ2OG9GbShKGxUiTuB6x7DEa/yjLN3gkLd3jH7WWOtS94h1TVBL0EiemnEki/1sP0RfEKl3vSeK6OiJaWRNhvcDwW7fZ8djYm3xAdDV+PfXdv/LvOs8epOb5rbysDe8dDIe7j6wj6voMT6xi5H0B+nWyO6N7z9pqkz+I569b7Yf3iHe9rtu2Jvc5kQhV3CCZx/9Z/R1Q7kwLq0Gb3Ozg4HXU7XOfB/fDbnCIvHJJbck54J1UaYCrcjoPTEhrc2xI4r1Pg5yjN/vEkjtlvIB/Zjkur52HkhPHcbG1ojsXnV9TFw4/OyFcnNTJVPr+z7esfX8TelnbWbE12jk63+3a+dtespM8gfs16U5ar+1aweVdT7KHgfO6Ivt811NDSlra/wxveco7R9ZPjYUj3eYuNthZRTz9M3CVt3bFLr96kqhvu8OcpRwLZT4YBqfsXdnhE2RjDJnfVvzQXg/NAOOmgIdx83uGx5avr4iLoxBWdpmtFacQ3t3lIH0ucykqExy4/HkjsXHKL9Kade5m3bieDelfExNed9eCI1Riv6Pts75YvuOzeGrf7ZU8npzd9cYQr28pdsmKVaxupaLPLPlzy6TFceeqBsVo0Tj2W/z7Wivt+8ZiRKbeRivaoSbieBtjZQO6Hr/tBviWN6Le2RxPGDjjpfdsbWjDG8PyimqSYvB9+zsogJ1XQdU62NTTz0PuJo7c/M34o1X0rOHh4X5ZtqccYkxTT9+vw3ry7iXL7ODzw3joWrN9Jc1t7LOa9yr5GUzlSx+w3MPZ6wv++wJf/8W7SOk5t+9UpzrnXk3cYOaAXm3c3xT+3VxjRzwrzrt+e7BTtaW5LCn299v1Tufn8wwB47gPLURtYVUZFaSRW0mXssL7859snMnpQL7Y2NLPWdlBinn5EGNqvgppdmc8/EJSCFn33AJMtrma09zor62A6unOPtG7+B95bl/IC69CWEv+m4XZPqOUZTyZMJh08J42r5sSxQ2LFwr7q8nwcD+TRy49nzdSzGTWwl+9vOHxkfz4xqj8iEuvjcHsf7gflCVNfZUVtAwN6lbGfne2yxrVNJywxrF9lrIkPcMdrK2ITnTjPT++guKW/nMyEEf34xv1zeHVpXPhjnr59hw7t518hclVdx+fHSUkc3q+SH04eHxsK74hASUSo7luRVT2U1vZorJQDuDz9PfFtuUMqHYm+u7+mrCTCgKoytjY089qyOq58cB4Tf/Vyh53XfimYjnc+Z+121tlpyb+bsSyh9szwfpXcffEnAThkRD927W1l4869sfCXE94567DhjB/e1715jIGPauId8ufd9nbCJDcfbdqFMSYm+t//7Diev/qk2OcdVQAF6/r62qf2pW9FqW8fUarwzuA+5ayobYh9x2lh9q8qY8zgKt5cXpe0rT3NbQnXMsCYIb35wlGjAPjjy/E6/1+eODoh5n/4qP48+q0TKCsRbptp9Ve1tJtY6/+g6j6s2daYkOG2bHM9T87fEMq0ql4KWvQBnv3OiUBiPM17k3RmFObdb63Oyo5UpZS98XVvKeJMYn2OB+Juxjq/caXdtB5gi8/4Ef1izXQ3e1vbYwPRnE60G55YxCx7ghK/mGpVeQkHDbVCZe4MImfd8pJIQsfqy0tq+fo/Z/HastqEjiw3lWUl3P61owH4waMfxLzKdpcggxWGuHbSWEZ5OuFX1SVmMvnREhMta+dTv3g4t375CA6zxR+guk9FUidnJuxuaqOfKy7uF95xH5N0LZNWlzA4DOlTwdaG5oS48/43Ps8rS1L3P/iJfj87z3zanA2c/LuZSTb+5OwJ/OPCY2LvnWPz4cZdsTCFs43BfSp44dqTefOHpyXs42WPTZfeZ41L+MaJ+7NpVxOPzo2L2snjqjnElSkHiZ31DgdU946VfgArv76+uc23z8zEPP3E+9vx5I/8xUv253EumDiad1Zui03M49DQ3O7byd2rvCR2HbW2R9nR2OqrJ8P7V3LJp/fn0bkbmLm0lplLa2PfGzHAuobdGW4vL9nCdx9ZmDbLLVsKXvSd8IPbo/LGkzsziXq2oyH3tlgX9x8uOIKJrqart2CU4106A8u8LQEHYwy9ykqYNME/3fRvdgbMHa+tBOLx4KNGD2DTriaenJ9YeXRvazSWhuZOR/vKne/x2xeW8s375ybt48NNu2Ie2c+eWRwTXCdbo6wkwifHDEz63t1vr4nVHYqIMOcnkzjrsOHM/+kZgFX697avHs22PS3c8PgiolETG+Dkzsi6dtI4Xr7uFA4e1pdJE4ZxxiHDYg+5dHg7SKvKS/nC0aMS1hk9qBfrU1Q6TeVZX/XgPBau35ngnTsetVtQ/21PXt+rrIQPN/rPNwDWg9gbFhzcu5yaXU1JOf6X3jeHf7y+kkfnrE+yz/HMLzp+Py47+QAmTRgaC2E57G5qjYV8AL5x0gEcMXpA7L3jyV/+wDxufXEZQFK1U2eMg1fzpl9zUsL7688az8T9BvLTpz6MPQjc15zD8T6ln1/93qlcdnK8BMKZhwynT0UpNz//UdLvdvyUiMegq08/KOG9++OLTxjDmMFVXP/4Bwkd7rv3tiYMyHIz1Q5ROg7biyly86+dNJaDh/Xlkntns2xLfazv6CsT49VKnbDm72ZYx7ijKEQ2dLvoi8hkEVkmIitE5Iau3t/g3uWURiTBg565LLH5lqq0q5tfn2+d2GEpRu52hJPBMOmQYfzkvw4BYMKIfqzdbnl6d7+1mjE3PBdrXj9mVxu8/IG5vvONLqmpZ29rO2ceOjy2zB3y+INn4nfnwv/KJ0cztG8F331kIWfc+jpLanbT1NpOU0t77MbzHo/bX1vpm22zZXczIsIPJx8MwGf+8DpTpy+NtWpKS4RHLz+BNVPPjsXNAd74OH78K8tKGNKngju+fkwsZg1WZsZ3J43j8XkbOOBHz8duggkjEr3ByrISZnz3ZO66aCIHDe3Duu2NtLVHWbh+J1OnL+W5D2pYVdfA2yu2xkTB6+n7se+gKtZvb4yFAFrbo8xcWsvn/vwm+9/4PB9s2Jmwflt7lOfsDviJrgddn4pSBvUuZ86aHbH9P7XACuEdf+BgltTsThqx29TazqxV21hZ2xBroTmUl0ZYsH5nzIteM/VszjrMugZumb6UHzz2AV+58z22uSemscXrqH0H8qPPTeCuiz6ZNFnQS4u3JIi+l8qyEv70lSMBK2Q4ckAv307fN394Gm/f8Blevu6U2LIJI/qxZurZPHnlCfzzoomUlUS45QuHJ3SM+k1e9PNzD01pj0P/KiuZ4L1V25PGozgj0b3n+YhRA2Idys7vcehdUcqvzz+ctdsa+e4jC9jT3EZbe5Qtu5tSjto/dB+rFeSI/TdOOsB3varyUu6++JOxeL+TBOF+KM1avS1h/g/vAysMgk+51AlEpAS4DTgD2ADMFpFnjDEfddU+IxGhLWp4fN4GJh82nOkf1nT8JR+++ql9+dGTi7jjtZVcceqBNDS1sc+AXkSjhsfmbuCMQ4bx5oqtjB7YixH9eyVcIHua27jv3TWA5dEcOXoAa6aezY+fXMRzi2owxvCLZ61DMGv1dq6dNJZ9XTfBhf+cRXXfCkYOqKK6bwWTJgyNCczJrnr/JxyYPGlzv8rSWJ8EWMPHn736RC67fy7Lt9Rz1p/fjH2WbryCw4kHDYndXHN/MgmAK089iKP3HciUO9/j76+v5H77tzrlaoGUfQnp8sKvmTSW2vomHpy1jsfmbqC8NMKRLu/Ty8HD+tLabvj+owtjwpq0zdPHZjQPwX6De9PcFmXjzr2MHlTFH178mL+/vjL2+Tl/e5szDx3GZScfwDH7DWKNHRe/8azxCROhiwgnjx3CUws2ccPji/jmyXFBOP+okby6tJYXFtdw2sFDufie2fTvVcbrroeiN331ouPH8KYdG7/0xP0BuO2rR/NRzW6qykt4Yt5G/jZzBWf+6U0u+fQYDtmnX6xjvdLTahjRvzI2mnv6hzW8vMRKj0xVf+i8o0Zy0NA+LKnZzecOH+G7jlu8n77q03zsilMftW/8YTh2WF/m/HgSc9ftYEifioQaVw7D+lWyZurZfPeRBTw5fyOzfnS67z4vmDiah95fx4X/fJ9JE4Zy2Mj+vL96O++stEI03jIlpSURTh8/LHYPjRyQeN2fcNAQrp88nt/OWMpby7cyuE85bVHDIR6Hw+HA6t5UlEZ41u7MPSLNNbrPgF7ce8mxScv/9tWj+Pa/5/N/ronfndHDYSNBy7R2amcixwM3GWPOtN/fCGCMuSXVdyZOnGjmzElfo6QjxtzwXNrP10w9O6PtnHHr6yyv7Th8AFZ4pldZCf2rylm6eXesI9S9r9/NWMptM1cmffeeSz7JaQcPxRjDK0tqufHJRb65w2d/YgS3ffXohGVXPjiX5+3yDCMH9GLjzr385OwJvt7H5l1NXP7AXBas3wlY6WVX2NUD97a0s2D9TvYbXMUJU18FLO/35etOoba+iVEDkz2zrQ3N/OrZj3huUQ2njKvmros+GftswfqdXHzP+4wd2ofZa3YA8PJ1p8T6BNLR2NLG1Q8t4Jj9Bsbs86O+qZXDb3oxYZmIf4XGkogw/ZqTYpUrvayorWfSrW9wwTGjOGlcNdc8PB9j4O9fP4Zj9hvIL579KFZ+YtTAXtTVN9PcFmX6NScltUb2NLfxuxnLuPedNbFlR+87gEe+dTyn/f61pAlz+vcqo7IswqQJw7jytIOSROmB99YypE85kw/zF97nF9XwjzdWsdA+rw7TvnU8x7oqQra1W4XufvrUhzw8Oz5hSqb3Q3cRjRpqdjclHQc32xqa+eurK3j2AyujyX3e5/xkEkP6JHYMb9jRyIm/mcnnDh/ObV892re1P3ftdh6bu5GVdQ0M7VvB7y84wvfhBHDubW/HjvfsH0/KqCPay3WPLOCJ+dYo+wuOGcXPzjk00FSoIjLXGDMxaXk3i/6XgMnGmG/Y7y8EPmWM+Xaq74Qh+uu3N/L8ohpumb40tuyPXzmCXY2t9K8q4/yjRqX5diJ3vLaS37ywNGn5OUfswzMLN1ESESYfOpzBfcppbGlna0Mzq+r28F+fGMGXjhnFAdVxkZu7djtfvCOeivabLx5ORIQvHj3Kt1nXHjXsaGxhzpodrNu+h/OOHMnQfolNzrr6Zp6Yt4GnFsSH3X9w02djnW5+tEcNG3fsZfSgXr4X/ytLtjC8f2WsGdsR0aiVB++3ra0NzTwyez2XnXxAVjVzOuL211bw2xeWcerB1dxz8ScTbGhrj9LQ3MaabY2MHtiLwX3S35g3PvEBD71viWF5aYRp3zo+oaXR2NLGtNnrmbN2B5V2C+7rx+2XcnuzVm3jNy8sZd66nTx82XEcd8BgtjU088ic9dTVNzOwqpyrTjsotOkRd+xpYXltAytqG9i1t5VvnrS/75iUddsaY525L193MgcN9X8Q9gSiUcO2PS0M6VNuz5dBSqEOkxc+3MzlD8zl7MNHcNvXju74Cz58vKWeL97+DhNG9GOanTYdhHwR/QuAMz2if6wx5jue9S4DLgPYd999j1m7tvODZFLR3NZOeUkkozh+R0SjhnZjAomXMzFI2LRHDUtqdjO0b0XSg0HJnA07GtnZ2MqwfpVZeW+KkitSiX63xvSx4vijXe9HAUnBV2PMncCdYHn6YRqQqjRDNkQiQiRp6Efn6ArBd7brTkFUsmPUwCpGJScgKUqPpbuzd2YDY0VkfxEpB6YAz3SzDYqiKEVLt3r6xpg2Efk2MAMoAe42xizuThsURVGKme4O72CMeR54vrv3qyiKohTBiFxFURQljoq+oihKEaGiryiKUkR0a55+NohIHRBeon56hgBbO1wrd6h9wVD7gqH2BaM77dsKYIyZ7P0g70W/OxGROX6DGfIFtS8Yal8w1L5g5It9Gt5RFEUpIlT0FUVRiggV/UTuzLUBHaD2BUPtC4baF4y8sE9j+oqiKEWEevqKoihFhIq+oihKEaGir3QKCWMigi4in21TCp+ecv0VneiLyEQRGZprO1IhIv1dr/PxIko9BVfuyevrWUSG2P+7fiqnLNB7IzD5fG/EyOubJExE5FAReQf4GTAgx+YkISKfEpGngbtE5H9EpMLkUS+7iBwvIo8CvxeRQ/JJuETkWBF5ALhFRA4Xkby5rsWiSkQeAp4GMMa059isBPTeCEY+3xt+5M3N0Q1cAzxpjPm8MeZjyB9vQUQ+AdwGPAY8CnwGOCinRrmwvb+/YZXE3op1LP/H/ixnx1BEIiLyM+AuYDpWqfCrgCNyZZMXY9Fovx0iIleAZXsOzfKi90aW5Ou9kY58uvC6DLtZbbBODiJyvoiMAnrZ73N9co4FVhhj/gW8BFQC65wP88C+w4CPjTH3AH8AngDOFZFxxhiTK/uMMVGsukwXG2MeBG4G9sOaoCcvsD39EcAW4FLgChEZYIyJ5lr4RaRERAaR3/fGMeT3vXEEeXhvpKMg8/RF5L+B8cAcY8x/RKQXMBf4AfDfWIWPNgMtxpjLcmjfPGPM0yIyDOtC/j1wEdZcwkuApcaY3+TAvlOAJmPMLPv9YOA9rEntV9lCcTXQ1xjzvRzbVgm0AGXGmGYRmQb8yxjzn+60y88+EYnYDyZE5Cng28D1wB7g/4wxK3Npn/2+EpgPfJ/8uDe89g0F1pM/98Z5wCHAQmPMcyJSDbwDTDbGrMzlvZExxpiC+QMEuBzrIr4EWAZ80/7sWixhvdB+PxLrZJ2VY/susz/bH/gt8P/s96cA/wGO70b7+mJ5KtuBu4GBrs9+BfzJfh0BTgTuAEbkyLZBzjF1rVNmn9NxObj20h27ccCt9utzgN3APKAC62GVa/t+CKzJ8b2Rzr7xwG9yfG9UA08Bb9j3cC3wJfuzqbm8Nzr7V1DhHWMd9eOBqcZqbl0FnCYik7EupFJgqL3uRuAtIJpj+04RkbOMMauxYpUb7NXnYl1Yzd1lH5bH/CrwdWATcIHrs8eB8SJyurG8121Y4rArR7Z9CWLH1GECsMUY87GI9BWRY7vJNj/73MduEzBORJ4Bfge8Dqw1xjQbY1rzwL7bscI51ZCbeyOdfcaYpVjCv95elIt740DgbWPMycaYvwPfA75rf/YQub03OkWPF30R+X8icordrAKr6TdSREqNMS8DHwCnYl1U3wEuEpEj7Q61SVgeTs7ts5uxM4Cf2XHAKcChWBdQd9g3wBjTjNUp+jLwMTBRRA62V10IPAz8SUQOAk7HarmU59C2cfZ6zlzPg4BGEbkYy1M9vCtjqpnah+XFbgJWAccYYz4PjBaRY7rKts7YZ4xpILf3RkfHD+BF4KYc3BunikgV1oPmfnt5CfCR/QewCOve+HN33RtB6JExffvEDwf+jeWNrAR6A1cAnwcOB+41xiwVkTHAn4BfGGPmichXsTomDwV+ZIxZnAf2/RH4uTFmgYj8G0u8SoGrjTEf+eyiq+y7xhiz1V5nLFYMtdkY80vXd38IHGz/fdMYsySHtjUZY37l+u4tWDHze7Ga2x+EaVsW9sWOnYj0N8bscm0n4X2O7PMev69gdUzm4t7o6Pj1wipYNhSroz4n94aIlBhj2kXk68A5xpgvu777Q6xQ3ni64N4IjVzHlzr7B5TY/8cBD9ivS7GaqPdhxXXvBi4E+tuf3wv8yrUNyTP77gN+ab8uA4bnwL6/Ao971j3ftvsgrAs/Yi8vzzPbquxlJwBfycNj1wuosJdH8tC+3th9Czm6N9LZN9Z1fktzeG884VnnfuDL9uvhrm10yb0R5p/TLM577Cb8L4ASEXke6Ae0Axhj2kTk20ANVs/6v4HzgFHALVhP7HedbRn77OSRfe1Y2TEYK8a7OQf2XQ1sEpFTjDGv28ufFJEJwAtAH+A0YIkxpiXfbBOR04wx74RpV5j2ET92ocfJQ7YvF/dGOvumEz+/S8iTewNoAFaLyC+AL4jIZGPMhrDvjS4h10+dDJ/Ap2DFlO8AvonVgz4ZKxvnWNd6VwEv2K8/ATwHzAKeBHqrfR3adwUw0/X+Auz0QmBosdmm9ql9fvZhhZeasMaI/Amo7ir7uuQ359qADE/MSdjpZPb72+2TcDEw114WwYrFPQqMsZcNAEaqfZ2ybxqwv+t7JxWrbWqf2udj335YmTx/Ao7uavu64q+nZO/MBaZJvKbF28C+xph7sZpk3zFWs3kU0GaMWQNgjNlprPQztS9z+9qNlT6KMeZNY8ybRWyb2qf2ue2LGmPWGmNWGmOuNcbM6wb7QqdHiL4xptFYOc1OoaozgDr79SXABBF5FitftttPhNpXmLapfWqfx765kBelHwLRYzpyIZYfa4BhwDP24nrgR1hpmKu7yXP2Re0rTNtA7QtKIdln7JhPT6VHePouolgpjVuBT9hP4J9iNbveyuVFY6P2FaZtoPYFRe3LF3LRkRDkDzgO6wS9BVyaa3vUvuKwTe1T+wrlr8eNyBWr7OuFWAWsurP2RkaofdmTz7aB2hcUtS8/6HGiryiKomRPT4vpK4qiKAFQ0VcURSkiVPQVRVGKCBV9RVGUIkJFXyk6ROQmEfl+ms/PE5FDMthOwnoi8gsRmRSWnYrSFajoK0oy52GVwO7UesaY/zXWbGiKkreo6CtFgYj8WESWicjLWLN+ISLfFJHZIrJQRB4XkSoROQFr8vLficgCETnQ/ntBROaKyJsiMj7FeveKyJfsba8RkV+LyLsiMkdEjhaRGSKyUkQud9n1A9uGD0Tk5zk4NEqR0aNq7yhKNog1F+0U4Cisa34eVvGsJ4wx/2ev8yusUZh/FWsC82eNMY/Zn70CXG6MWS4inwJuN8Z8xmc9767XG2OOF5E/Ys3e9mmgElgM/F1EPos1M9SxWHOqPiMiJxtj3uiyg6EUPSr6SjFwEvCkMaYRwBZrgMNssR+ANXvUDO8XRaQP1jSMj7pEvSLD/Tr7WQT0McbUA/Ui0iQiA4DP2n/z7fX6YD0EVPSVLkNFXykW/Iae3wucZ4xZKCIXA6f6rBMBdhpjjsxin85Q/qjrtfO+FMu7v8UY848stq0oWaExfaUYeAM4X0R6iUhf4PP28r5AjYiUAV9zrV9vf4YxZjfWXKgXgFVLXUSO8K6XJTOA/7FbE4jISBEZGmB7itIhKvpKwWOsGY4eARYAjwPOjEw/xZqj+CVgqesrDwM/EJH5InIg1gPhUhFZiBWPPzfFep2160Xg38C7IrIIeIxgDxFF6RAtuKYoilJEqKevKIpSRKjoK4qiFBEq+oqiKEWEir6iKEoRoaKvKIpSRKjoK4qiFBEq+oqiKEWEir6iKEoR8f8BXWGeRAqeKuYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "dvs['04021520']['Q'] = pd.to_numeric(dvs['04021520']['72604_00060_00003'], errors='coerce')\n",
    "ax = dvs['04021520']['Q'].plot()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### get a single site"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "http://waterservices.usgs.gov/nwis/dv/?format=rdb&sites=04015475&startDT=1880-01-01&parameterCd=00060\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>agency_cd</th>\n",
       "      <th>site_no</th>\n",
       "      <th>datetime</th>\n",
       "      <th>72578_00060_00003</th>\n",
       "      <th>72578_00060_00003_cd</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>datetime</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1978-09-19</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04015475</td>\n",
       "      <td>1978-09-19</td>\n",
       "      <td>143.0</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1978-09-20</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04015475</td>\n",
       "      <td>1978-09-20</td>\n",
       "      <td>125.0</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1978-09-21</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04015475</td>\n",
       "      <td>1978-09-21</td>\n",
       "      <td>100.0</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1978-09-22</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04015475</td>\n",
       "      <td>1978-09-22</td>\n",
       "      <td>80.0</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1978-09-23</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04015475</td>\n",
       "      <td>1978-09-23</td>\n",
       "      <td>83.0</td>\n",
       "      <td>A</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           agency_cd   site_no    datetime  72578_00060_00003  \\\n",
       "datetime                                                        \n",
       "1978-09-19      USGS  04015475  1978-09-19              143.0   \n",
       "1978-09-20      USGS  04015475  1978-09-20              125.0   \n",
       "1978-09-21      USGS  04015475  1978-09-21              100.0   \n",
       "1978-09-22      USGS  04015475  1978-09-22               80.0   \n",
       "1978-09-23      USGS  04015475  1978-09-23               83.0   \n",
       "\n",
       "           72578_00060_00003_cd  \n",
       "datetime                         \n",
       "1978-09-19                    A  \n",
       "1978-09-20                    A  \n",
       "1978-09-21                    A  \n",
       "1978-09-22                    A  \n",
       "1978-09-23                    A  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = nwis.get_dvs('04015475')\n",
    "df.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### just get the url for dvs for a site\n",
    "(`start_date` is 1880-01-01 by default; `end_date` is `None` by default)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "http://waterservices.usgs.gov/nwis/dv/?format=rdb&sites=04015475&startDT=1880-01-01&parameterCd=00060\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "'http://waterservices.usgs.gov/nwis/dv/?format=rdb&sites=04015475&startDT=1880-01-01&parameterCd=00060'"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "url = nwis.make_dv_url('04015475')\n",
    "url"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Fetch some miscellaneous measurements"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['04021520']"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sites = field_sites.site_no.tolist()[:5]\n",
    "sites"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "04021520\n",
      "http://nwis.waterdata.usgs.gov/nwis/measurements?site_no=04021520&agency_cd=USGS&format=rdb\n"
     ]
    }
   ],
   "source": [
    "fm = nwis.get_all_measurements(sites)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "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></th>\n",
       "      <th>agency_cd</th>\n",
       "      <th>site_no</th>\n",
       "      <th>measurement_nu</th>\n",
       "      <th>measurement_dt</th>\n",
       "      <th>tz_cd</th>\n",
       "      <th>q_meas_used_fg</th>\n",
       "      <th>party_nm</th>\n",
       "      <th>site_visit_coll_agency_cd</th>\n",
       "      <th>gage_height_va</th>\n",
       "      <th>discharge_va</th>\n",
       "      <th>measured_rating_diff</th>\n",
       "      <th>gage_va_change</th>\n",
       "      <th>gage_va_time</th>\n",
       "      <th>control_type_cd</th>\n",
       "      <th>discharge_cd</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>site_no</th>\n",
       "      <th>datetime</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"5\" valign=\"top\">04021520</th>\n",
       "      <th>2005-04-27 12:11:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>1</td>\n",
       "      <td>2005-04-27 12:11:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>KGG/GRM</td>\n",
       "      <td>USGS</td>\n",
       "      <td>4.19</td>\n",
       "      <td>129.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>0.00</td>\n",
       "      <td>1.1</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-23 12:20:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>2</td>\n",
       "      <td>2005-05-23 12:20:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>GRM/WWL</td>\n",
       "      <td>USGS</td>\n",
       "      <td>4.22</td>\n",
       "      <td>135.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>-0.01</td>\n",
       "      <td>0.9</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-26 14:51:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>3</td>\n",
       "      <td>2005-05-26 14:51:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>DLR/GRM</td>\n",
       "      <td>USGS</td>\n",
       "      <td>5.01</td>\n",
       "      <td>199.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>0.02</td>\n",
       "      <td>1.1</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-06-07 14:45:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>4</td>\n",
       "      <td>2005-06-07 14:45:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>DLR</td>\n",
       "      <td>USGS</td>\n",
       "      <td>3.91</td>\n",
       "      <td>113.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.7</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-07-22 12:50:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>5</td>\n",
       "      <td>2005-07-22 12:50:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>GRM</td>\n",
       "      <td>USGS</td>\n",
       "      <td>1.61</td>\n",
       "      <td>5.82</td>\n",
       "      <td>Fair</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.5</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                             agency_cd   site_no  measurement_nu  \\\n",
       "site_no  datetime                                                  \n",
       "04021520 2005-04-27 12:11:00      USGS  04021520               1   \n",
       "         2005-05-23 12:20:00      USGS  04021520               2   \n",
       "         2005-05-26 14:51:00      USGS  04021520               3   \n",
       "         2005-06-07 14:45:00      USGS  04021520               4   \n",
       "         2005-07-22 12:50:00      USGS  04021520               5   \n",
       "\n",
       "                                  measurement_dt tz_cd q_meas_used_fg  \\\n",
       "site_no  datetime                                                       \n",
       "04021520 2005-04-27 12:11:00 2005-04-27 12:11:00   CDT            Yes   \n",
       "         2005-05-23 12:20:00 2005-05-23 12:20:00   CDT            Yes   \n",
       "         2005-05-26 14:51:00 2005-05-26 14:51:00   CDT            Yes   \n",
       "         2005-06-07 14:45:00 2005-06-07 14:45:00   CDT            Yes   \n",
       "         2005-07-22 12:50:00 2005-07-22 12:50:00   CDT            Yes   \n",
       "\n",
       "                             party_nm site_visit_coll_agency_cd  \\\n",
       "site_no  datetime                                                 \n",
       "04021520 2005-04-27 12:11:00  KGG/GRM                      USGS   \n",
       "         2005-05-23 12:20:00  GRM/WWL                      USGS   \n",
       "         2005-05-26 14:51:00  DLR/GRM                      USGS   \n",
       "         2005-06-07 14:45:00      DLR                      USGS   \n",
       "         2005-07-22 12:50:00      GRM                      USGS   \n",
       "\n",
       "                              gage_height_va  discharge_va  \\\n",
       "site_no  datetime                                            \n",
       "04021520 2005-04-27 12:11:00            4.19        129.00   \n",
       "         2005-05-23 12:20:00            4.22        135.00   \n",
       "         2005-05-26 14:51:00            5.01        199.00   \n",
       "         2005-06-07 14:45:00            3.91        113.00   \n",
       "         2005-07-22 12:50:00            1.61          5.82   \n",
       "\n",
       "                             measured_rating_diff  gage_va_change  \\\n",
       "site_no  datetime                                                   \n",
       "04021520 2005-04-27 12:11:00                 Good            0.00   \n",
       "         2005-05-23 12:20:00                 Good           -0.01   \n",
       "         2005-05-26 14:51:00                 Good            0.02   \n",
       "         2005-06-07 14:45:00                 Good            0.00   \n",
       "         2005-07-22 12:50:00                 Fair            0.00   \n",
       "\n",
       "                              gage_va_time control_type_cd discharge_cd  \n",
       "site_no  datetime                                                        \n",
       "04021520 2005-04-27 12:11:00           1.1           Clear         NONE  \n",
       "         2005-05-23 12:20:00           0.9           Clear         NONE  \n",
       "         2005-05-26 14:51:00           1.1           Clear         NONE  \n",
       "         2005-06-07 14:45:00           0.7           Clear         NONE  \n",
       "         2005-07-22 12:50:00           0.5           Clear         NONE  "
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "fm.head()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:xlabel='datetime'>"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD8CAYAAACb4nSYAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy86wFpkAAAACXBIWXMAAAsTAAALEwEAmpwYAABT00lEQVR4nO2dd3wcxfXAv+9OXe62bGxjycbYGAw4gDEtVFNMd0IAgzCm/HDoDoGEohAIQQmkASGBxAFCsegQMB0MIfRi01wBF1k2Nu5VXXfv98fu3u3d7Z1Od6diab6fjz66nZvde7s78+bNmzczoqoYDAaDoWvga28BDAaDwdB2GKVvMBgMXQij9A0Gg6ELYZS+wWAwdCGM0jcYDIYuhFH6BoPB0IXIam8BmqNfv346dOjQ9hbDYDAYdijmzJmzXlWLotObVfoi8gBwErBWVfe00/4InAw0AEuA81V1s/3d9cCFQAC4UlVfs9P3Ax4E8oGXgWmaxCSBoUOHMnv27CRu0WAwGAwOIrLcKz0Z986DwISotDeAPVV1b+Ab4Hr7R/YAJgGj7XPuERG/fc69wFRghP0XfU2DwWAwtDLNKn1VfQfYGJX2uqo22YcfATvbn08FHlfVelVdBiwGxonIQKCHqn5oW/cPAxMzdA8Gg8FgSJJMDOReALxifx4MrHB9t9JOG2x/jk73RESmishsEZm9bt26DIhoMBgMBkhT6YtIGdAEVDhJHtk0QbonqjpdVceq6tiiophxCIOhdaiogKFDweez/ldUNHeGwbDDkXL0johMwRrgHe8akF0JDHFl2xlYZafv7JFuMHQMKipg6lSoqbGOly+3jgFKS9tPLoMhw6Rk6YvIBOBa4BRVrXF9NROYJCK5IjIMa8D2E1VdDWwTkQNFRIBzgefTlN1gyBxlZWGF71BTY6UbDJ2IZEI2HwOOAPqJyErgJqxonVzgDUuH85GqXqyq80XkSWABltvnMlUN2Je6hHDI5iuExwEMhvanqqpl6QbDDop09PX0x44dqyZO39DqDB1quXSiKSmBysq2lsZgSBsRmaOqY6PTzTIMBgNAeTkUFESmFRRY6QZDJ8IofYMBrMHa6dMty17E+j99uhnENXQ6OvzaOwZDm1FaapS8odNjLH2DwWDoQhilbzAYDF0Io/QNBhevzvuel75a3d5iGAythvHpGwwuLp4xB4AT9z6xnSUxGFoHY+kbDAZDF8IofYPBYOhCGKVvMBgMXQij9A0Gg6ELYZS+wWAwdCGM0jcYDIYuhFH6BoPB0IUwSt9gMBi6EEbpGwwGQxfCKH2DwWDoQhilbzAYDF0Io/QNBoOhC2GUvsFgMHQhjNI3GAyGLkSzSl9EHhCRtSIyz5XWR0TeEJFv7f+9Xd9dLyKLReRrETnOlb6fiMy1v/uriEjmb8dgMBgMiUjG0n8QmBCVdh3wpqqOAN60jxGRPYBJwGj7nHtExG+fcy8wFRhh/0Vf02AwGAytTLNKX1XfATZGJZ8KPGR/fgiY6Ep/XFXrVXUZsBgYJyIDgR6q+qGqKvCw6xyDwWAwtBGp+vQHqOpqAPt/fzt9MLDClW+lnTbY/hyd7omITBWR2SIye926dSmKaDAYDIZoMj2Q6+Wn1wTpnqjqdFUdq6pji4qKMiacwWAwdHVSVfprbJcN9v+1dvpKYIgr387AKjt9Z490g8FgMLQhqSr9mcAU+/MU4HlX+iQRyRWRYVgDtp/YLqBtInKgHbVzruscg8FgMLQRWc1lEJHHgCOAfiKyErgJuA14UkQuBKqA0wFUdb6IPAksAJqAy1Q1YF/qEqxIoHzgFfvPYDAYDG1Is0pfVc+K89X4OPnLgXKP9NnAni2SzmAwGAwZxczINRgMhi6EUfoGg8HQhTBK32AwGLoQRukbDAZDF8IofYPBYOhCGKVvMBgMXQij9A0Gg6ELYZS+wWAwdCGM0jcYDIYuhFH6BoPB0IUwSt9gMBi6EEbpGwwGQxfCKH2DwWDoQhilbzAYDF0Io/QNBoOhC2GUvsFgMHQhjNI3GAyGLoRR+gaDB6ra3iIYDK2CUfoGgwdG5xs6K0bpGwweGJ1v6KwYpW8weGDcO4bOSlpKX0SuEpH5IjJPRB4TkTwR6SMib4jIt/b/3q7814vIYhH5WkSOS198gyGDVFSEPgZHjIg4Nhg6CykrfREZDFwJjFXVPQE/MAm4DnhTVUcAb9rHiMge9vejgQnAPSLiT098gyFDVFTA1KmhQ62qso6N4jd0MtJ172QB+SKSBRQAq4BTgYfs7x8CJtqfTwUeV9V6VV0GLAbGpfn7BkNmKCuDmprQoSLWcVlZOwplMGSelJW+qn4H/AmoAlYDW1T1dWCAqq6286wG+tunDAZWuC6x0k4zGNqfqqrIYxHvdINhBycd905vLOt9GDAIKBSRcxKd4pHmOVomIlNFZLaIzF63bl2qIhoMyVNcHHGoTnGNSjcYdnTSce8cDSxT1XWq2gg8CxwMrBGRgQD2/7V2/pXAENf5O2O5g2JQ1emqOlZVxxYVFaUhosGQJOXlUFAQOgyKWMfl5e0olMGQedJR+lXAgSJSICICjAcWAjOBKXaeKcDz9ueZwCQRyRWRYcAI4JM0ft9gyBylpTB9euhQi4ut49LSdhTKYMg8WameqKofi8jTwGdAE/A5MB3oBjwpIhdiNQyn2/nni8iTwAI7/2WqGkhTfoMhc5SWwnUvAaDz5kFedjsLZDBknpSVPoCq3gTcFJVcj2X1e+UvB0x/2dDhMVOzDJ0VMyPXYPDATMg1dFaM0jcYPDDLMBg6K0bpGwweGJ1v6KwYpW8weGB0vqGzYpS+weCBce8YOitG6RsMHhiVb+isGKVvMHgQNJa+oZNilL7B4IXR+YZOilH6BoMHRucbOitG6RsMHhjvjqGzYpS+weCBGlvf0EkxSt9g8CBodL6hk2KUvsHggYnTN3RWjNI3GDwwOt/QWTFK32AwGLoQRukbDB4YS9/QWTFK32DwwMzINXRWjNI3GDwwKt/QWTFK32DwwETvGDorRukbDB4YlW/orBilbzB4YAx9Q2fFKH2DwQPj3jF0VtJS+iLSS0SeFpFFIrJQRA4SkT4i8oaIfGv/7+3Kf72ILBaRr0XkuPTFNxhaB6PyDZ2VdC39u4BXVXUUMAZYCFwHvKmqI4A37WNEZA9gEjAamADcIyL+NH/fYGgVjKFv6KykrPRFpAdwGHA/gKo2qOpm4FTgITvbQ8BE+/OpwOOqWq+qy4DFwLhUf99gaE3MKpuGzko6lv4uwDrg3yLyuYjcJyKFwABVXQ1g/+9v5x8MrHCdv9JOi0FEporIbBGZvW7dujRENBhSw1j6hs5KOko/C9gXuFdV9wGqsV05cRCPNM+qparTVXWsqo4tKipKQ0SDITWM0jd0VtJR+iuBlar6sX38NFYjsEZEBgLY/9e68g9xnb8zsCqN3zcYWg2zDIOhs5Ky0lfV74EVIrKbnTQeWADMBKbYaVOA5+3PM4FJIpIrIsOAEcAnqf6+wWAwGFpOVprnXwFUiEgOsBQ4H6sheVJELgSqgNMBVHW+iDyJ1TA0AZepaiDN3zcYWgVj6Bs6K2kpfVX9Ahjr8dX4OPnLgfJ0ftNgaAtM9I6hs2Jm5BoMHhhL39BZMUrfYPAgnYHcFRtr+HjphgxKYzBkjnR9+gZDpyQdQ//QP/wXgMrbTsyMMAZDBjGWvsHggXHvGDorRukbDJ4YrW/onBilbzB4YCx9Q2fFKH2DwYOgUfqGTopR+gaDB2YTFUNnxSh9g8EDo/INnRWj9A0GDzqSob9sfTVL1m1vbzEMnQSj9A0Gm31ueT30uSMtw3Dkn95m/J//195iRBAMKpuqG9pbDEMKGKVvMNhsqmkMfe5Iln5H5I5Z37DPb99g/fb69hbF0EKM0jcYiB24NUo/Ma/N/x7AKP0dEKP0DQZilXxHcu90RHxibYQXDLazIIYWY5S+wUBstI6x9BMjttI3jeOOh1H6nYGKChg6FHw+639FRXtLtMMRvapmJlRZpmP9P1i8PqPXSwefveO1aRx3PIzS39GpqKDppxdzX9EPqPNlwfLlMHWqUfwtJFp5ZWKP3EwrxLPv+7j5TG1EyL1jtP4Oh1H6OzplZbxUvC+3jr+Iuw4520qrqYGysvaVq6Uk01tpxR5NjJsiA7osEwpx7da69AVpBWydb5ar2AExSn9Hp6qK7TkFAGzO7x6RvsNQUWH1TpYvt8xjr95KMnnSoDUGcjOhEFdurk3/Iq2A49PfWF3PtU9/RW2D2e66RbSjS9Yo/R2d4mLEVlDqmF92+g5DWRnU1LAtJ595/Xex0qJ7K3aeCDLYo4lR+h3E0s/ySfOZ2gFHqrvfWswTs1fw2Cc7kJHR3rSyAdMcRunv6JSXIzk5gMsjUVAA5TvQ/vN2r+TC037NSef/laCjUty9FftzEKEuKycmPV2iLftMKP1MXEPomErfaYv8tqFR35S52M1l66sZet1L/HfR2oxds0NhGzDj/+9eTjn3L1ZaG7pk01b6IuIXkc9F5EX7uI+IvCEi39r/e7vyXi8ii0XkaxE5Lt3fNgClpcjkyYBt6ZeUwPTpUFrazoK1ALtX8vngUQA0+v0R6e7PNx57CaOufjasojPUo2mNgdxMXMPdeTtxr4FpXy9TOO6dsG8/c879L1ZsAuD5L77L2DU7FLahsqTvEL4aODImvbXJhKU/DVjoOr4OeFNVRwBv2seIyB7AJGA0MAG4R0T8Gfj9Lo/vwAMBCJ5/AVRW7lgKH6xeSUEBWQHLL9zoz47trdh5Hh9j2Qq12bkZ7dHExOln4JqZjmzJze44HXPH0nd6IsEMjuiGI4MydsmORTxDpY1csmmVIhHZGTgRuM+VfCrwkP35IWCiK/1xVa1X1WXAYmBcOr9vsGgNa6tNKS2F6dPJslVtoNijt2Ln8amTZ2hGezTBJ56IONb/vZP2NQMZ1lod6fWKRLqdWkNB77DluTnKywkUFkamtaFLNl3T4U7gl4DboTdAVVcD2P/72+mDgRWufCvtNEOaOJbRDl1HSkuRHlb0UeCLL72VeWlpaPwi8NVXmevRVFSg034WmfaPe1MbWHOdE9hn37QH59yKryNt7BIaX7b/BzIoW3i2byeltJStf58ePm5jl2zKSl9ETgLWquqcZE/xSPN8ryIyVURmi8jsdevWpSpil8Hv6xwTZZz7aEqwoItjYGbUii4rg9rI0Eitb2j5wJoTlWET+G5V2lEZ7vvsSG/Xces4lTqTDZLjKvps+aaMXbOj4Tv99PBBG7tk07H0DwFOEZFK4HHgKBGZAawRkYEA9n9nCH4lMMR1/s7AKq8Lq+p0VR2rqmOLiorSELFr0FkmyoQarwSBIK2i9KuqIsNdgaBIywfWosJKm3z+tKMy3A15R3q/vijNkUmDo1+3XCBcHjoj7blmUcpKX1WvV9WdVXUo1gDtW6p6DjATmGJnmwI8b3+eCUwSkVwRGQaMAD5JWXJDiPCKhx1IK6SA061PaOnbtmUm3QkUF7OoaGhEkoq0fGAtqpEI+Pye6S0h4HoUHakn54uK3glkcLXN/BxLLRV1z83cRTsY7VlVWyMc4DbgGBH5FjjGPkZV5wNPAguAV4HLVNVM48sgHUkppII/ieV6HSXTFMjgvZaX8/Eu+0QkaW5uywfWohqJJp9H6GkLiejRdKDXGwrZxBlPyrxPf0c3YhLRnnU1KxMXUdW3gbftzxuA8XHylQM70KyhHQOn+OzwSt/XvBXvdPgzeq+lpcia3LAjEtCLpkLpGS27Tnl5pE/f50s7KiPSvdNx3m+04yWT7jbnNjuxzm/Xd9lxAn8NKeNYWTt6JXH8xIEEpr4v5ALK8M2OGRNxqIcc0vJr2GGlDk2Dh6QdleFWph1J6Yfi9FthPMkpz/EakhUba/j6+22Z+8F2oD1fpVH6nQBHGXSkkL5UcBR6Qv+wo2RauYVL+VG6FHzTy6+kHZXh7vVkOu4/HXwxcfoZtPSbueahf/gvx92Z/jyK9sQofUNaOIZxukohGFT+9c5Sttc3ZUCqlhNW+vHvo7Us/Wh3RcrRFa7wzJPufo/6R9KM0w92TKUvMZZ+a7h3Mnu/yzdUs+j7rRm9ZqoY944hLYIZcu98sGQD5S8v5Mbn5mVAqpbjS0KBtErIJjB2aJ+I45QGiqPi9AG2XXVNxuL0+3brSNEskRMCM6v0E7t3UuXwP77NhDvfzeg1U8UofUNaZKriOdbt2m3ts3FHSyz9TCuEbH+krZ/S9T2Wf9baugzG6XccSz96u8SM+vSd/+leM3rN+g5Ee77KjETvGNqXsE8/vev4kgiZbE1CCj2J6J2MxukTGw2ZkvvIIx5fxTs9WSLi9DuQeyfGp59B2ZzynNY7tntdq/wF1PYexPDlyzMkXWYwlr4hLYIZsvTbe99TSWKQtq1iuJtSmW1kx+P7gu7pJylM8nIRMZDbcXR+q87IJRPl2e51HXzpg4y/6J+ZkSuDdLbJWYa2pKKC4A03ABD84MO0/MfRXfa2JhSnn1DpW/8zPZAbfc8pXd9e/tmv4QYjUFCYVpy+49/2Scey9MOLomU+XDgUvZNOj7MV16Z/4L1lfLhkQ1rXaM9IO6P0d2TsLqxushamCjY0pLXAl6+dF25Lxr3ja6WQzehKmJJP347Tdzs+HvzlnUwfcmDKcjlyZPt9CZenaGuio50y+T6cV5HWuE2m16Z3jQ/c8uICzvrXR2ldzlj6htSwu7BBsV5jUCStBb7a270TVujx8zjT/jM+OSuKlK9fWoo/Ly90+M/1efzu5UWpXauigsDV1wCQU7OdQNXK1K7TCviilmHIbJy+03tI45p2rysjRO9p605PEWPpG1LD2TfWUda28k+1axsOmUxbstR+P4llGBwZW30gNw0HekZWh7QVTXDjRgByGhsIfvVVm22e3RwhV6Cz8U0mXodtTQePmwBAMGq56xYRNTuakpK4WVdsrOGuWd/GV8QeUVmh9BQxlr4hNewubFjpS0R6S2l/S7+ZQdpLL0VWWA1a4ORT4dJLM/bb0becaCmI5siI0rcVTcBuyLODTVZD10abZzdHdKRX2mXGZU2HfPrbtqfXyLlmQweWLoub7ZKKOdwx6xsqN3godggZUXcfdCZnTSqPSU+F9twcxyj9HRm7C6u2YlDSW+CrtWLgk/994v/+pZfCvfcitpYJANx7b8YUf/QM3HTcR9Ex/ylhK5SQ0g80WT25Nto8u1lCA+rW+0hbcbmsaWdvgwCSsUauoSl+I+706moa4sxEt42oPx82mQ9LxsSkp0J77pNglH6yRE/06AjdbLsLG+xtzSYN5OVlZNu19nI3SqKBXLurHorTd5YtdnfhM0h6Sj8D1crpxdn3mR1o5N1h+3LbSVekf+0MEB50t47TNhTsxmxdQS9Wd+8HQNCXuUaurjH+Ku45Wdb7aozno4ozPtDw23SissKf27pnbZR+Mthdz4+C3anJyrEGdNLcBi9jlJaGQzZH75mWws/IAFoahENGPX4/YFVanx0OeenE69mQ3yOUnjbRIZtpOKmzMmHp24rGsfRzApYV+o89jk3/2hkgOooqbWvVbuT2v2IGN0ywGrZgM3McWhIxVJ/A0nca6cZ4czOixwds/rbT/t75kzAQ3UW8rXvWRukngZaVscaXx6Szb+PqE39uJaa5DV4mCa+9k+YyDK+8av2fO7ddejMJV9n0WxavuG5xYf9hofR0iX5y6fj0s6NnLqWC04vr2xeA3Ew0JBkkeuG7tN07HtZ0wOdP6KpsyWB+IveO445rTJDHy5has8VjuZLoSJ84BmJ7Lq9hlH4zLF23nWGT/s67Q38AwNwBw8NfdhD/qlNm0io7FRXo738H2APC7dCbSRinby9k5nNNfGry+WMWOEuV6J9sTMP6yoilD1BaSuAGy7DI2W+fZjK3LeGF7+wxlnStVQ9rOpiTk7Dn2pLfbHju+fBBlEHjWPoNLZyFLV6v2R6b+L5bX6qz7dBdDwPRreiNpd/B+GipFTI3c/fDAWjyu5YryvQEkBQJd7HTm7audfXWdZzQzzbuzYQmh3lVgnvugUsuiUhqOulkKz0DRA/kBtJx72TC0nfksN+p43fuOFjvasm6aiBDg5FRCj5UDuPQEmVZf8tvwwdRBk2Ovxmffhw8lb5tCB542UP8ePKfYtIdgsa903FxrLZgVjbg2vc0zW3wMolTZtKKXa+qCkVNqHu+ZRv2ZhJG7wDccw8yarfQYdPUn2bst6N/srUnfyWL0wDmZGJw2MUpf3uPRz9O8d1WVOCzFWYgEwZHHJorzy15Rw3RA7kug6ZZn35cPLS+yxD8umioZ7pFWPa2LmtG6SeiooLAL68FoCE7B7CVfklJRqJkMkXIp5/mtHVH2avbhGnD3kxSq2y6ZGvI4Apk0UornSUPMlmFHT2UaUv/q5VbuOE/c1t+ou2z9m3bEpEcXLU6Q5KFaa4daZF7x58dm2gbNNlZqSl9T0vfK9LHw0A0ln5HxC7c14+zFPva3B4AaLduUFnZYRQ+uDadSHfaur18QGiSVxv3ZpJZV8ddz9IZbI0meiAyHesrk5NtnHeakTDQTBBndmpw6dI2FyUZZelMlGvI8lD6tkGTvbwSgMbJU1oUwOA5B89rJrCHgegu421t6Zv19OMRVbg35XcHIFhb314SxSXk3knH8i0tRav9sNT2pZaUWAq/DRu3ZCx9Ny31wSYiuv1I61lmkGBQ8UlmlX5alqVtHWuUayPQ0JiOSCmRzH1k+4VAUGko6Bb5hWPQVFSQ8/4cGD2eRn9W2N8PzZZ98XLvOOdd95L1ubLSM0uEpd/GZS3lkiQiQ0TkvyKyUETmi8g0O72PiLwhIt/a/3u7zrleRBaLyNciclwmbiAh6UyoqqpibWFIdPZasxiAYEY775nBcU2kazHosVYMeHDosHbpzTS7DAORXep0YumjiXXvpGPppytNmIAqfp9kZmkHm5b7rl2Elv6IVB3qWmSurUjGOAhF5lzmmtTmtr7Lysiut9b4afTZNnCSAQzbFn3LLr98gf/tsl+L9Yu7N9jWq6emYz40AVer6u7AgcBlIrIHcB3wpqqOAN60j7G/mwSMBiYA94hIZoKsvUgyXjYuxcW8NTw8+cIpENqKIqeKo5/SVvqh0M92mpxll8ZkLdFMVpbM+vQz2QNRfCJkZVDpp1VOQkt/RO2cVTI0PaFSIBkL2RkArz/okHCi26CpqiLbnvjW4I7MSyKA4bnNOQR9Pv520Jkt1i87pE9fVVer6mf2523AQmAwcCrwkJ3tIWCi/flU4HFVrVfVZcBiYFyqv98sZWVsbYJbj7yQdQW9rLSWhCCWlzNrt4MYvGUtQzZ/T40dcxvMyUlPrlZYzuG9xeuAFHd7cuEo+7Ysguu2hd1l0VP7vXB3qTPq3om6VDoVMZN1OGAr/QhLP82yk1Y5cSaNde8eKWfffqlf01Un9tqY/LaGyVj6TvRd3MlZxcUupZ8dkZ4soYi+FugX3dGjd0RkKLAP8DEwQFVXg9UwAP3tbIOBFa7TVtppXtebKiKzRWT2unXrUhOqqoon9z6G+8b9iHeG7RuRngx1Z0zi/eFjGb9mIfmN9VQX2j79dGaA2r2P7avX8vcDfkJT1YqMTICa991WIH2LwTl79ZY6Ntc0pHUtKir4YNyxvDts37gK6vOqTexfPov/fG6tEx/eOSs5pZRuI+cmxtJ/7/2U30sme0pBtZ5L1uJv3T+Q1uS5pJWMrYxX9yhi9tijwr9VWkrwjDMjsqZ8z64e+cWnXsfcPiVJn5pMOXHcO/WushIha3k5ubYWbHQs/RYGMES4upLULzukpe8gIt2AZ4CfqerWRFk90jzvVlWnq+pYVR1bVFSUmmDFxRy15FMAstx7libZgn+4ZAO1Koy/7Rfk77M3NTsPdWRLTR4IDQ7ffvgU/nj4FMt9lMEJUJly7wAcWvZc6takXZHPPmoak8/8bVwF9c2abQChredC/tdE0+FdpOWbjiL47nsRx00NTSkr1UxW4aBaA7n+D94PpVX2Gmh9SLHsJDUWYr/D17N34qBLH+QnR18d8Tyii1rKRc+uE58O3oNXdzsk9vsEzz+Z15/jUaYiZC0txX/SiQA0+rNTCscOuCfjJalfgrqDWvoiko2l8CtU9Vk7eY2IDLS/HwistdNXAkNcp+8MrErn9xNSXk5WruWKafS3fELVrIVrKMjxc8CwPuRm+6m2l11Ny4izrYDqHCuOd1tuYUR6qhw20moYGwLB9BqlWbNCH7flFqZuTXqF9XkoKIlaa8fxW8eNv6+oIPDNN6HDxi++aplcCdAnn4w4bvL5U2+QFU7aeyB/On1M2nIFgtZAbtaWzaG0/w4fG86QQtlJqrG03+HU024MJanreUSPW6TcAFdVsbJHEaef8wfv7xOUv2TGXZz5DZFKP1J23XtvABqvuz6lAIboCZuJVvQM/eaOuAyDWDX2fmChqv7F9dVMYIr9eQrwvCt9kojkisgwYATwSaq/3yylpWT9zlLwTS1swVWVtxat5dAR/cjL9pOf7aemwXqRacXCO3HBASu8LdSdTGcCVEUFvrffBqwGqXbGoylfSu+/LzYxFcUXTxFFpTs+/OiGytPSr6iA888n2BSuUPUvvpy2a2z1llqGXvcS73SLfAch6y0FpapYDVomJlSFone6FYbSmnzpLQWSlGXpcd8N/qxwyGbUJZLqnXmNZxUXe0+cckhQ/pJRlo5P//ZXw1tWRp/n3Etza+/0y/L+PhTibOuXq5/8slm53O3VDqP0gUOAycBRIvKF/XcCcBtwjIh8CxxjH6Oq84EngQXAq8BlqpqhdXG9yfrJjwG4bsIVfP/lokiFn2BAdf6qrazeUsf43QcAkJftC72YtCx9O/LBcTc1+vzpTYCyu+CB+vBg6Paf/zJ1Rbg2zvhJSxVfPEUUle6EnjsNqfNsPa3GadOgsTGiK13ny7LS0+DLFZsBeGav8RHpIestBaWqakWxZ2LphFD0zuGHh9Kc5ZZTKjsVFTQdNT7i2BOP+67Jzg+HbEZVhERLF4d+xyua7oQTkObCPeOUv2QaGq9Y+ug67Bw2d731Td7vs2mPPSJ6CB8t3eD6LW+FEene2UFCNlX1PVUVVd1bVX9g/72sqhtUdbyqjrD/b3SdU66qw1V1N1V9JTO3EB/3Erd3vRl2CzQXzvnWorWIwFGjrDHo/OwMhWnakQ/ZhZZ7p6lf//SWcwhtqReWb1sw9d2GdMAA7y9aqviSnIoe3p4xMuv97y2LveYGqyIFXINm9Vk5zMkr4q5Z38bmT5LcLO932+TLSrlBtix9yHVZ+tvqUpu85Lh3/GP2Cqf5U1wKxC73jd+vDafFc594vMPqnr1DzyP6nTWr9OO5/F5+Gd+ttyY+N075q2tsXll6hc9+bY8lhTNZeZIdS4om2lL3uSKt4j0X9xnJuIMySQeZ2906+F1L3LrXbGnO5/zmwjX8YEgv+nXLBSAvU0of2HbaGXx7wmkANN54Y3oToJwt9VyN2337T+QPJYfHOyMhesEFsYluxZdsuGnUVPRgyVBPBeWPWlUzmfj2oOtea7NzOe2cP3HHrG9SXncoNzuO9Zabm3KDrGpFLbjdO1tqU1T6Ghun33TLralNnrPLfZMrAi1QW+dtJHgsdVxdflvoN6Mt/YamZhSXXVZfGHUoQ699kdqs3FC6b+KpcU8LFBbGNLxOY5qc7zw2beLf34/MY/9/7ovv4l5n7dbw2vkDeuRGfBc9MO4OQ44no7sHsK0uzjaNrUTnVfoVFWSP3iN06PvWZQ1WVaHALUf9H4v6lUSkr91ax5crtzB+VP9QciaV/oUPzub9xZbVmnacucfsyMd+cDz3HHRGSpfTo46KTHBbky2d7OZSSHXffOupoKL35HXqwdiS3jF5sTcTcfdq6rLCcyaagpqSxRRvKn1g2C4pN8iKpajdSyekOns46Fj6rsYu5TWHbMX7+ogDQ0kN/uz47ruo+68/7vjQ55hxGMclV1HBw0efy6qe/SMNA7us3vHDswH4rkdRKD2Ry7Tx3n/GyOEo/WZ7FyTnjg27FZVVm2s98zguyD0G9uDjG46O/C6BwRFPobtP2WqUfga49FKYPJmsyrCbQN5+O6IAbszvwQP7T6TUvbt9cTFvLbK6vo4/HyA/J3NK/5PKkLcr/ZBDuwve5PPTvW57mpLFVpDFn8wLVzjbSlxb2IumFq63X9vgrYyj19pxft5zmd677oKcHALi4+zPX2Hv1d9QlxP2Bc/8chWjbnyVhasTRQ3HEs/nms4SD8EgIJGWfks36HAIOHH6Lks/5WACW/GOWB9W8g3+rKTddxERMMFI91V9UxAqKth6+TR+PfZMzj6zPNIwsMtqbpPV46nLygn1IhP18Jpc8wGueOxzTr77PXJtIyyZRj6ZJZ/dv9+ckp5ycOw8grgbqie8Xvg3P1++qVkZM0nnU/oVFfCPf4BqRHy+r7EhrKDKy8nKs7po0RMy3ly0lsG98hm1U3jG4bzvIpeRzRRprx3jzI7My6N7g7eF0hKipblp5rzwQVUVW3MKGHf5DG4ZPzUivTnWbPVepC7aveNQ49VIlJbCAw8QzMrCr0Hy/D7q9tw79PXr878HWv6u4inQdAfXBIkYyI3X8DWHs+Cae0ZuyraCrXgLGsOuivruPeKOWySKqgqqUuAyhlSxDQOrHK7ubvXMQoaBXVbz7FPqdy4O9SITeebcWxi+8OUq5n63hTzbJVfXnEspAZXrrc1f/v3+Mp6eszKUvjXO2Iuz5IPXBjmJLPV4Yznue3728/hupdag8yn9srKQyeruuIvjkgCrAN55JxC5Pn7dGZN479v1jN+9f8QYwHuL12dGtihXSNO8eXEy0iL/edNee9N9xC5pixddyevdA2XFxVTn5gPw+sgDI9Kb44S/vut5D/GidzyVPljbB/bug//yyyg89CC29ewb+ioY1VtIlngKJ50wOlVFoiz92hQH60Jx+n630k9R6ztGQv9wL7bhD3+K68aKGax1z2oFeuRHhVq6NuIJ+PwR6c7v5x9grWdV+/hTod9NNLfEqze8cbs1W7zuxpubrR/xLv3lys0A/OaFBRFGyeYabyXdaD/zlm6FGc/Sb6/1raAzKv04luem/B5WSIUzo/C0nwDQlG8PFJWV8eHuB1LbGGD86vkR5xZkwqfv+MRdNL3+hndhtfM2rPiOeUXDmvWfB4JKt7z0V8mOLoZD+riiN8rLIS/fzpfEevvRsnrcg+PeefvrdSxYtTXUzY6r9AmHMNavXMX8VWFXTmDld9430Qzxuv/pjLcolsGRmwmlbw/k+lxGSFpjQaWl6D//GTqsP+mUuFmjn43b6lZVCnMiy5wWF4fGl+JtK+r0DtzPI1H76uUWq7bLR9226mbHl+K5jgI/u8ozf7wB96YEln4ittV7Xy+62KW1AVIL6XxKP47l+fzoI1jQb2jIxeMU6MZAIDRAOWv4OAoaajnglz+NKBA/3nfn9OWyfeKHLZ0TSmoKqLdP3M57y/iLOOn8v7KyR1Gzk1S6Ryl9L0v1mzXbeHXeah76oNJbxqhTIlZ2LC1F/nB7+Li5kEEvWaPuwd2bevzTcGOdyEcaUMX/9UI+qI60MgPz5sc5IzHRFteXvz6WCw4Zlqalb9kX7oHcujTcO36fMLhXfiitPg23BkS6hxKFKUZb2W4F7KwJ5Kbp1nKCzYTqOkERke84kaUf/7tn9hzf7PhSvNfYtHFTjBEG4aVBYvJn2NKPblC3JyjzmabzKX2vGHGbr4tKQj0BpzAoAjU1KPDW8P05dNnn5G3bElGALj1yePpy2b+b2xReyGx7bj5b1myIm/ezwaMA2JLXLSI9moCH1RWtGL5ds41j73iHi2d8xk0z53t2iaOtosJ/3hORT35sTXbTQYOaDxlMYlZuzFR0+7C2MRC3+xsIKr5336FH1MC12pWypcsaRxuS4rM23kh3aWVBMuPesWfkHj4yvAZVquGfDm6FE0/pNzQF2ePXr0WmnT05VB6cNYEivj/9TPQvrsn5HoZBvscgbEKffoIBjKV9d+bF3Q8LJ3iUubjlyFlmw8Z5V/HGXsKWfjNKP6pObftotmc2555/epjllt0Sx63UGnQ+pe/EF3ushhkUX6gnEF0YFvQfxuoeRYxfYq8M4SpABdkZ2GDM/t3cQPjlPr3XMYyZ9njcvE5XWTQyPZpAUGMskPo9x0QUwK3PvxR5kkeX2P1IetdsoSYnLyKf8wtJGcFJzMp1X6YxEN6eRjX+xJugKv7Nm+lVG2mROaGcSQ9yVlSwdcQo1p89JSLZt+ee+BcuoK4xmLJydSz9TPn0fSL4fMLfz7ZWi+2em2DZgiSIUPpxHphXb6LBtbNU8LtViEhEeG19U5CgbRgAnoaBEwlXXR++vnsg9crxI0LzY6D5CLcF/YeFDzzKXLyiGojaBMYnMLyokPXbvYMOnKUrshLNsrbdst3rqxmz6msAtr78uqcbyXkHA3pYEWjpNuQtofMpfbAK2kMPxVj8gby8uDMK39z1AESDoZU53QUoLycDjykUspbEksVRG1X4NJjQfx4IKv7Kyohr169eE6HUA//6V8Q5dVk5MV3iCKVfu41qew8BJ19YKSeh9b1kjboHTWBxVsfp7gaCir9Hd3rUV0ek+7DOj1dpI7Ar55k/vIzrJ1wR8ZWvajl1r78JwJjfvN78tTxQ7LV3EkXvJDlQH7QtfYAT9x4YMzEoFZKx9L3cW6E1cmpq0MWL8Qncf97+nLiXtepnXWOg2VmtjqVfXR9+v9PfCe+vm+0THrlwHMP6WWsNNTd+sbZbH+tDvPqhkWMrDtH10C9C/+55rN0WR+nbjU92Iks/NEPex/4rF9CnZgvbfDnebif7+fYptOaabDVKPwN4zChsOPc87xmFBQW8OXwcP1j1Df1qtsQUoEysoRJagqEgP+armArmRFnkWBVcBg1M6D8PBBX/Z3PIa3TNBIxS6vXrI2OBN+fZIalud4v9v2ftNgobaqnJcclaVRVqFJKKE4+W1aOr7/agRK8Q6tXNVlWCCr4jj6S7K+wQYGNhLwD+8sY3MefFYFfOhQNiI558qhlYotqK3sl+/LFQWt3vbw8r9hZMdAsEFb9r7KNbbhbb69Pz/zbn01+8drtnmK17c/FgXT0+EXrmZ3P0HtZExuuencsPb/9vwt92xnHi+bCzs3zsPrAHvz11T6B5S39dYe+E40tBVbrnxfaMQvMEbKobAvTvkcuarXUxeSFcR50G+LWfhd1KIwdEul/rs3LIbWqge30123ILPN1OTh3qbSt9Y+lniqhCULfPfqHPbqW/9u//4stBIxm/5FPPAhSxhEOa8vjPPism2TOWt7SU4K67AtD0yqsJ/eeBoJK1bQv5TeGKWu/MVrULXN3AyP1qNjvjBB5d4oonfkVBY13Y0rfzOc/MvU3du/94grpddm3WYg0sXRZzD27VakXvhPGK4HF0sX/M3uTsvlvkl0OGxOSPS4K5Bb5gkILGNOY8VFSg6zcg//gHcu7kUHLt9tqwYrcbned3P4yabNtyjzcQGYzcF7h7XnbcWPJkcZd9r1mtR//lfxx35zsx6Vtzw0oymJ8fkiv3ww8AeOeb2AX7vt9SF9GAOw17dZyGy/GZO+7KxmZ6DnWHH5lwfEmBfI+een3RgBijcEAPy9L36sk2Rrl3SvqGn0Wo4SwuptHnJ+Dzk99UT/f6mvDy6f36RdSRkKVfYNXTzUbpZ47/XnNE6HN9RMiZK8/oHwIw/pnprb4huM+jAYnXyjsiJprEpaqs3VbPur6DyGt0uXccq8xW6nWlkyPO25zfPb67JS8v0tIPzZy05bEL7NIHHmNyZTfOPfD/WFvQK8Zidd+ql0/b3cOp3FDD+675EF7WrNva6jG8JOK7/Qtd129u85cEcwt8GqRnqrObbQteg0FrXoi755KdG1bsVVV823cI0075Jb844Wfh8z0aI7d7B6B7XvqWvkYofQ+3UxzWFYZdKYFhw63wxYoKCu++M+45B/7+TSbf/3Ho2HmHbp++Gye6J7ShThxLf4+BPYDmd09TJaKn1KvAqhc1194QU8/7d8+loSnI1trY5+s0Pk6v3700S+h9lJdT26OX9X1TA7lNDVT2HmQJsWFDRK8u8InlRu5daMljLP0M4vgGIdJl4LZ2Zn65KmYWbqtQUYE+8khMclylb4uYKJJk5SbLKp21897kBcPXqffnRCj1X2/oFXHe5iG7xPRoQk/kV79C8vOZt9OuEYulOdaJ0zX97p/WVsifFO/FuMvt+7IVWzCoqBIalHtq9gomTf8wQobosLX128ONlpc16+T3ifDbhoUR3zW95XIrNLcuUIIIL58qkxa4rtWS3cNsC96K07dknfG4Zb1XOw1oVRUUF4f2XJ43wBUZ5tEYNQUjlX633Cy2p7hWixMx49aTb3/tss495pK4WdO9b6gn3NC3H9lZPigro8+mtZ75HYU827XMgNPOx2u4nDj+0KzbOAPgt5w6mr0G9+TLlVsSKn5FGTkgXK9fuPyH9OuWw/INVuROUXerfE4bPyL0ee22WBePM8ZUkBsbILKpptGqG6Wl1N3xV0v+xnoW9h9GZZ9BfN+tL3MHDLeWUgfL7//Sy4D1PrN8YpR+a+G2atxu20+XbWL/ob0z58bxwq5QM0YcFvPVzC+8NxBzlFwiS99tCbkXj6sbPCRCqW+KCgnbcufdse4W+2fkhON5c7C1lO/LM993zZy05bIf3sZ4g6ZVVaHeQA97/sBvXljAR0s3JtzByM1f34xdLjlk6X/xOb0uvYhr3nk49F29L6oyJloXqLSU7fdO9/xK+vahd/12fvrxM+Q21rds9zBngxGRUMTVD5d/yW7rKsNjKMXFUF5OtW0RVkf1pqJpCgYjxpRS9enPuvcJRt34Kkv6DSF4ww2h9P+4lwDwWn3WxZpDx4d6wo2BIDl+gaqqmEgqBy+XhfPO483FcJR+r4LEvu68H53KXHvJjf+NOihu4xwMQre8LCpvO5HK205kSJ8CduqZxyZ7D+iCHD9H7z6Anx09IhRJ4zWY6/RMokOjwSqXTjz+ff33seRrqucXdvm8++AzOfm8u7hl/EWhc3TzZsBahrlXQbZR+pnGsZTcYYDRYWvNraRZ3MfbMkyaBBXqPo/14+saAyFrJNFiXfEiaer//VBCN5XXdPO5320GIgdY3QN6oUbIVr6bB0VapiFJXP7/7lFT9TfXhq15rwiRwUHrnoe9/GxMRXZ6GL6XX4KaGgoawhbZpvyeoc+hCTsJfPfVE0/zTJdu3aChgSafn/rsXOr82cnvHuYKtfVp+CH2rtnKpoIeYcVeWsrWadcAsD2nIOFAZGOTRkzy6p6X3fKleCsqeONly8Vy39hTuW3vU0NfuSd9NbeO0oqNNaHy1hSw5Soupk+t9yJ3az0Gg51ysT2Oe2fxWsu11tMuN59Wbgrdg89VMHNWVnH5B1a489acgoSNc/RKqr0LcthYbZVDVctlJiL0f9uK1loz8cyYsrfdnlkbb+b7RrsRcSKRcgJNnD/nBfb8fjEV+5wAwKsjDw7lD/SxlhDxi9Aj3yj9jOMso+D2K0cry+Zm2rnDvlKaMp2gQu3ickE5vPtt2L+dqJIHX3gx9Plnc8Ox+PURKyLGyrt6S2wX9u//XQLA0vXbQ6Fk7rXBo6+y+fiTAbj27X8DsDW3MKTYnIZhyNbIrv/mJ/8Tluu9yHXNAXI3WZPVvu07JKYiO/fh32itVNq3Jry42oqe4fVkNhbYDUAC3319vA047Pd0//4TAZi5x+ER6Qmx3UZB8eF3lL4IfWq3srFH3wjFvnmstX5RXXZuwnGkxkDQcqPYdLN9+sEZSa7NBFBWRvF6Kxb+sR8cz3Z7QHbf9UtDPm4A+vTxPH3f9Us5a1wxW+uaQoZIYyBoKf3ycgri7Enw/fOx+yQ5Df3mBd+GZXdxgh3+WWhb/E/PWWn1bMrKIhZQzGlq5ILZMwG46mSrAfVqnJ1IKjd9CnNC6+MH7V3OqKig/88vB2BtYe+Ysre9rgmR+EuybKyObOCqiqzAgss/eCKUtsEplwUFBI6bAID/ySfouXAuW59/uWWuxDToEkrfqTSbqsNWZrQebG5NjauPDUeLNKYyWzOBAtpp/mcxL9vd8se1Aioq+OLe8HkTP3gu9Nm9LrhXo1G1MbbXMbzIanxGD+rJE1MtpdT3zj+GKmdw5sxQXlVl/S6j6O5TBuVYD3PNyL1Cis2J8Bn22nMRv7Hp1nDoYuDpZwAYvzi8VfKmfGuAbu7AEZYid1XkkHunl1V5Bm4LN4yVvQdFfm5m16uGgKVAJozeiaGuSAznPf3xpTsAwvMBkll+2A61bfJnWQqqpAQeeYTeU85m06CSCMXudtEkWvKhIRAk22WQ9LaV9OZpVxNcXmW5tZpzQVVVWQPLUfTZvC7c46uogK1bWdW9X0y+0/beibPGWUrMWbq6wVH6paXIdG9XWeUjT4U+P7rP8TB0KHXvWg39ysI+NIgle29XT2Hn3ta7cLta97zpNcadeEtEyGh+U32EWyl0d1GNs0LMzOHiPgV8v7WOpkDQnkhn7TbXbctGCutrWN3DfgausretvoluOVkRu2L9+7z9OXFvq5FyjKj+9rjAmeefACUlTFj8EWUfW6G7QZ+f9SNHW+HYY8YA4Lv0EvpsXs+6wl4tcyWmQZdQ+qUHWBXW8eNBrD85uxlLf8KeO/GrE3cH4LX5a1ouRILBw42SE/Oyr3kqvLlyXKVfVsYvj74EgDGrvmbo5tU8+pjlry1/eWFIWW98/GkAfnvqaJ699GCO3WMAlRuqYy539O4DyPIJu/bvxq5vzMQfDLCltiEUdaC/vimU9+IZc3jko+X069ONIZNPB6CqNmhVkoqKkEuq/9b1dKsPNzCfDBzF1gumWo2IbbHfPCu8ANimgp48/uh1ABx68f0s77WTVRmGDg01Er5TToaCAga5ehFuhbC4qLjZXa+cntDEfQbz9i+ODH9hv6eDl38F2PMZkt020Q7HbBIfvu7dQq6cPoU5bKppiOhxuaNXNiSYUGb5zsPV1BlsXO/L5ac/uoHdrnmOt3YZaymoc87xthaLi8MhvC76BOrYvG6TNcdhyhRobOTgSx+MyPP4i7/jrMt/wsgB3fEJLHCUfpOrMSotpag2dknrF4eHV2O94bjLWbmphtWrrXce8PlZ3NdqSALiY9cNKzitV0NEPTy3T/i5rHWWanZkr9mCD+XK9y2Fusxp9KMaZ2sXs8i63X/RVwQV1vfuj373HbJsSaixGLWuktmDd49sRCoq2PbIY3Rbtzri+R45qj9/OcNS3pc/+jlguYpO3GsgReedbfXggkEu+u8M3rjKGsv798gjYfJkGv9iGRVZ27ez6/oqFvUfxpbcwlBDs6m6geUbqltlNc4uofR/fsxIDh3RL8K6ijbWE06vtnH8jVc+9nnLhSgtRf85PcIv6bAxv0eEVRE9MSnubD2XVXPTm5a1tfta1/iAraw33WotlLZznwL2Le7NiAHdWLa+OmYwraYhEFq4TX5VRs+67fzzgNMsxQsE6sONptPw5S6Yx65X/B8A3/QtDlkrjbaCzgk0htwJAH8+bDJ7X/Ukc+uyQstM+IOR93vginlMf+a3AFx10tVW4vLlBH/5Syv/gQfA9OnslAP7fhcZxQOw+KwLmw27dQaUHbfdUxcfxBVH7Rqy1vv2s3oca3ceFtmAxJtJaw/U161cRZM/i8L14RnRvQtyCGpkRJL72S/36HU5NAYil9gosqOh1hX25o2RBwFwwek3s8dVT1lLCyxfbin/fv3CspWX05CXjy/qOQ/esIrq7Dz2vaKCT3eKnPcwZtU3HFz5Jfsv/BifT8jL9rNLToCFjzwLPh/bN26l27LwYPvHd0/mremRkT+fDd494viHl/ybD0rGcOiyz/AHAzy830mAtbz5EUtm8+ffnoM8+mjoee7/6D1xn0u2fS8nLnoPgE+HjPZsnIPR7p2KCno8YM1OP/DSh1hV2Affe++HXFunLvgf83falY+G2HsS9+kDU6eyvTFIt4aaGGvcvcfy1rpGmoLqaUCOmDWTiYve5V+7jWd+0VBWBXPIbazHh1JUsxkgvKHT8uU8csaVHP7HtwnsMjzjln+XUPoiwi4bv2P1irVhV8XLL0fkSTi92kYXL3ZfNLJiJUG9Wnu8jl05n/ueviWUvrGgJ0EkpMS/i9qybUtto7eyKS7m+K+t7vK+9lofuYGwYr79sCkAXHuUVRmdwbF9hvRGFWb961nUuWa/ftQ+PIP871dZ11++nNymBhr92Rz+0/u494DT+P0R58fc06L+w+hRX03P2m384YjzWNGjP9TU0HiNpaCzA00Ub1odc95dh5zFdcdfCYBfgzxZcS0AT9j/j138MZd98ASfDd6d6eN+BECDrSizP/kISkvJWreWZycM5KVXfx+67u55gZAFGfHM+vWLmCBT/4q1mJizPs7+Q/uEXXilpeQtXUz/7rlUnXdJeLvIfv0sheqeSeso2GnTaKhrYI+fW72qQVvXhRry3l9YMdkbh+wSenf/dC098NXKLZHyikBWFohQv3krOd+GZxkX2YON6wojt5SsycnnrwdP4pnRR7Ehv4cVFz55snWtsjLq9voBBYEGDls6hx4NNRxTOYc+28OhlL8+5uKI6x3/9fs8+kQZ/iE7h2Q78NNZzBq8N3cddCY12bn0nvlsqPz7ioewy6ZVHP2tNWB8SOUXoWtd999/R1z7lAX/4/8++Q+PjzmOm8dPpSYnn5xAY6RPvqyMkjXLiebFB6dRebvVWHDJJYwoFIZvWME9B55BdWPQ6rE4dbNfP4Krv0ceeyxcT8vKOPzryNDhOvFBXR0UFHD63FkM2rqWW8ZfRFNOLmzaBDU1bMnrFp6/ETV2cN+5YwFrUcPlG2p44avY8k5ZGb964x/0qd3CieffzVN7H0O9PTHv7C+ssY95O+3KnEGjQITpI8cDWLv/Zdjl0+ZKX0QmiMjXIrJYRK5rkx+tqGCn555kW04BW+yR/k3lt0dkaXZstqKCG269IDJtwwar0ouE//z+iEobUjYi1P70UgBOWPR+eGE3rLXHv+tRxJbcQn7wsyd47KyrAPjjzpYFWPFxFR/dfEekspk8GV2+nC8GjgwpfiBigta9B1lul6251uzboUf/ECoqOHRkP0blBbhyRSGnHXo58/rvQnDDRmoDSn5jnXV9EVY7+5gCtx9xPm8PHxv38Zz7uTWIfOglD/CXH55No21eZQeaeGf6RSz+wyns+X240Zzl2qc1p6mRcSvnU/mnUzng5MMsXzhw1XsVnLjoXX535IUMvfZFHhtjDX4Nu+N34UpQWsroL94LheSNHDWEJcvWQL9+fHfJVbyYP4QVPfqzZXsdW7bXhZ5fw21/ACD3rVlx72lo30KqFlaGlH3tlu18udMIJpx/N8dc+PfwjOUNG2DDBuYM3p2gHTq6j90Is3w5O995G2AvDrZ8ORsuujT0G4N65vHZu19ENiiABgIsKBpGdXYevV98zrrfigqKfm6tFfQze/DyiCWzeWqG1cDe9cOzufqkn3Pi+X/l1REHhWNsly+nbuHX5NfX8vAnD/DVfg3868mbqXfW0oGYJSnqs7IjLeeyMn76wZMA3HHoOdY9Lp8XsRsdBQXc8eKf+MPLd3L/078JXSvg87Hs9pN49pGrOfvzVzh28cdc9f6jjFq7jAfHWuv577bOVvBO77WqKrLXahNSvCUlcM89+Mpv5TfvPsjy3gN59AfHgz1W47yThqxscjdvDCvOqip61W1n9zXhRndR0TCoroYpU8gfvBM3vnU/CwfswtUTprHN3sB9Tbc+9KveHBbE1cveY5DVK3z78l8DcNLST2KVdFUV/Wq28O+nbmaXDSsAQvW2wLV8yhsjDuBPPyyN6B0nHT2WJBlYPjJ5RMQP/B04BlgJfCoiM1V1Qav+8LRpDBxgddf2vfJR3rjvEiafeWtElqXrm5mFWVZGz+1bOOuLV3nsBxN4eJ8TmTj/v+QGGsgNNDGv/y5kBQP0qK/m66ISRq9ZSlH1JmSDFY0SRJhtL5Vc0FiHAJW3n8R7JWM4Z1I5K3oN4Nk9j2JzXjfu338ivmCAn1x7HjdPe4LqoDDptN9w2txZ+FQpqt7EgVVzEZTVPYq46j27gJWU4NseeR+1WbmMWreMoupN9Pl2AUydSi7w8OO3c8ZhV/DZ4N056by7Qvn7OgVblcFb1vCdKyrGwZ3uWOg/f3cGu2xYyVUnX8NfDzmbZ/a0LJVBW63JP1kaZMiWNczbadeY6/WuswfkgkG4x+7SDx1K1vLl3DXzj2zNLeTdYfty70Gn07N2G2OWzQ1vwRfFyNVLeL4xn1v3nsh9dg/BTUFDLWVv3c8CW8nll98ChQHPa5VsWc1b6+vQDRsQ4PJTr+XNXceFvj/kkgc497OXKP3iFQZs38iyPpZf+dbX/s6wTfbcC7+fvSrn0a2+hstPvY6sQHl4i07gwHWLmRXsxZbtdXw5dB/+Ne5HqAjvDd0nlGe3Vd9C2bsAdN+ygd3WVfJ10VAAbnzrXwzf+B1v3HcJXxeVUJuVyx8PO5eLf1zGKQve5sRF77HHmqXUZOeR75530KcPu9rKx+HBfU8KfR677btIt1ZVFUNU+fjv5/LomAkM3rqWg6u+Ck+7tvN1LyvjjLlWQzr3jtMpP+r/OH3uLASrN7rv5hVw/hSYPp0XH5zGByVjqOq1Eyctsu4v5JMvLiZ7+XIqbz+Jt4fty3ln3MLpX73BkC1rYhqjQ5Yv54eVn1N+1IX8d5exjF/yMUv6DOH1kQeyNa+b5ZZxFGex5YKc8M0HoYYu1Li8/DJUVnI8cPUJF3PHnifw9tB9GL5xJUv7DmHigv+FH5Zr7GDgi8/SrxpeGTAagBM+ew1e+GPEc3F+d/d1lbx53yU0+rLICTZB377QrRt/eukvXHPiz/nHgaeHrnvunHBkXlLRY0kibbltl4gcBNysqsfZx9cDqOrv450zduxYnT3be03qFvwwcwaP4rRz/hTzVX5TPbVZudw6cU/OObDE42Qbnw9UmTtgOCe7lCTAAVVz+bh4r5hTCutrLNdH3XY253fnezsy4l/P3MIxdsTK8l47cfhP74s5t2/1Zub87RyCJUP5sjGXK0/+BSts37qbfts38b/pF1E4aIA1cFRRwX/+/AhXHWtZhPuvmM+nQ0Zz/Nfvc+9z9mMuKbEKkSob8nsw7ZRfhJTMIZVfUPHErwCo82dz/7gf8aN5b0UM8H3yt8ms7t6PPdcsCYcm2nzTr5jfHvV/vDtsX8as+ob/PHI1PntY7PNBu/GjyX9mz+8Xh5T/w0/cyGGVn4flqqy0PjuzQ2tqCCKs6d6HiyfewPFff8DFnzxjKRuP8ZHKPffniJNujkn3IqepgU//NpmeO/UL/66Lp46cxC8OmMw17zxMQUMdtxxtuckeeOpmVITyIy9gad8hZAcaGbZxFYv77kz3+ho+/+vZ1j0XFITmZiwsGsrPTr4mpKwB7nnu9/Sqr+bsM28lp6kxNBjds3Yb23MLCIpwxQePc+X7j5PlDC2q8n23vlz8o+u548W/hBsXFxvye3DL+Km8PuJAal2bx49Z9Q3PP/Jz66BvX6itZbXfitM/6NKHQvluef1ezt04L/KZ2C6/GNzvzOHSS+Hee2Pzdutm7V/tuMvs9xuioCDc0ER9P7//MEasX0FO755w111hZWrXyzp/Nncceg4v7H4oq3r0JyvQFNq964mKazlg5XyrzDzyCJxzDtty8tnrqqfIbazn87vPtqxtd5ny+fh8pxE8su9JVPYeSI+6an7/2t0M3LYhUk772Vwwdgpv2QbBG/ddwogNK+KWZ8/7HTqUy8acyUv2/gBHf/sxf37pL/R0ose8nnMziMgcVY3tnqtqm/0BPwHucx1PBv6W6Jz99ttP0wY0CPrknuP1/NN+rSXXvqgl176oZ51ZrouKSvTyRz/TzdUNia9RUuKspqJLeg8KXSP674Qpd2rJtS/q5NN/ozeNn6pXnzBNL/zxjXrsBX/T835ykz6297EatK+joI3iizj/zoMn6V8POlPXFPa28oiolpRo0L6Hel+Wfte9n765y1i9+6Az9NPBu4fz2TQ8MkP/eewFevaZt4au+9Yu+4V+07mmuuQIIDq3/y5am5UTTi8psfKCvlcyRt8rGRNxTsSfnU/9flXQ2qwcrfdlRV5rxgzVkhJdV9BTbz/sXN2akx/+vqDA+t7NjBmh68X8lZR4vycRvfqEaVpy7Yt6zfHTYs7bnp2nS3oP0g+H7KlrCnvFPDs3tdm5esrkP4ee4X6XP6Ib87pHXK+y10762yMv1J9OvF6nnXS1vj1s35j7dfLW+bP03/ueFLqe86w/GbyH3jR+qv7qmIv1xd0OiX+/Ue8s9Ne3r/UXlV6blaOfDxypFWOO09sPO1ffdb8/EUs++/iGYy/Vkmtf1MtP/oXW+rNjn8mMGdY7cv+G1ztzuOSS8Lvz+63jaJzn45RHr/ef6HvVmGcSBK3qOUBrsnK13pelW3IKIp+hauhZNYovtozGuW7oz++PlUNE3y/eW0uufVGPueDvkc842fsR0QVFQ3X4Nc/psRf8LfnnnABgtnrpYa/E1voDTvdQ+nd75JsKzAZmFxcXt/hmY4iqEPW+rLByi6c8opkxQzU7O6YQNIlP6/xZEYo86T9bUQYQbRJf/MruVeGSUYIlJbq09yD9pu+Qll/TKWjxCr9LwUcU4GSVQzIVuiXXc91zi99DvDJgN7ZLeg/SuQOG6+bcwpZfr7nnHK9R87rf5p6F+32JJCdjvOfl9UySfWdtSTJ1w+s5NVemWlLuWvIM49GSRiZJOorSPwh4zXV8PXB9onMyYunPmKHq81CqOTkte6AzZqgWJlnx4/05lbGkxLJ+klG8zm/Hq9DxCmOySkIkbC1GV+iWKt3o62ZCObTkevGUQGGh97tLdC/xrtW3r/XuPKzruA1cvLzxyoBTXqPvtyWNZTLlJZX329Fw36vTiMYrz9HnJHqOrWWYtNY1ougoSj8LWAoMA3KAL4HRic7JiNJXja14ffum/kDjVWJ3ZU2l8DWX1+ucZJRgugq4I1p4iUgkb0vvJRMuiObytvbzzeQ9GLzpgPUsntJv04FcABE5AbgT8AMPqGrCqY4ZGcg1GAyGLka8gdw2DdkEUNWXgZebzWgwGAyGjNMlZuQaDAaDwcIofYPBYOhCGKVvMBgMXYg2H8htKSKyDvCYCthq9APWN5ur/TDypUdHlq8jywZGvnRpa/lKVLUoOrHDK/22RkRme414dxSMfOnRkeXryLKBkS9dOop8xr1jMBgMXQij9A0Gg6ELYZR+LN4bfnYcjHzp0ZHl68iygZEvXTqEfManbzAYDF0IY+kbDAZDF8IofYPBYOhCGKVvaDEizh55HY+OLJuhc7OjlL0uqfRFZKyI9G9vOeIhIj1dnztiQcpuPku70aHLtIj0s//721sWL0zdSIuOXC9CdOgKkmlEZLSIfADcBPRqZ3FiEJEDROR54D4RuUBEcrUDjbSLyEEi8hTwJxHZoyMpLhEZJyIzgN+LyF4i0mHKtlgUiMhjwPMAqhpoZ7EiMHUjdTpyvfCiw1SMNmIa8B9VPVlVv4GOYy2IyN7A34GngaeAo4Bd21UoF7b19zesZbHXYz3LC+zv2u0ZiohPRG4C7gNewVou/DJgTHvJFI29p4WzI3Y/EbkELNnbUaxoTN1IgY5aLxLRkQpdq2J3qxXrBSEiPxKRnYF8+7i9X9A4YLGqPgK8AeQBVc6XHUC+PYFvVPXfwJ+BZ4FTRWSkqmp7yaeqQay1mc5T1QqgHCjB2qSnQ2Bb+gOBNcCFwCUi0ktVg+2t+EXELyJ96Nh1Yz86bt0YQwesF4notHH6InIWMApry7AXRCQfmAP8AjgLa/Gj74EGVZ3ajvJ9pqrPi8gArIL8J2AKsBJYCCxS1dvbQb7DgTpV/dg+7gt8BBynqkttRXEl0F1Vr25n2fKABiBbVetF5EngEVV9oS3l8pJPRHx2w4SIPAdcDlwLVAP/UtUl7SmffZwHfA5cQ8eoG9Hy9QdW0AHqhohMBPYAvlTVl0SkCPgAmKCqS9qzXiSN1x6KO/IfIMDFWIX4fOBr4CL7u59hKdbJ9vFgrBd2fDvLN9X+bhjwB+Bc+/hw4AXgoDaUrzuWtbIReADo7fruVuBO+7MP+CFwLzCwnWTr4zxTV55s+52ObIeyl+jZjQT+Yn8+BdgKfAbkYjVW7S3fL4HKdq4bieQbBdzeXnUDKAKeA96x6+9a4Cf2d7e1Z71o6V+nc++o9eQPAm5Tq8t1GXCkiEzAKkhZQH8773fAe0CwneU7XESOV9VlWL7KlXb2OViFq76t5MOymN8CzgFWAae7vnsGGCUi49WyXjdgKYct7STbTyD0TB12B9ao6jci0l1ExrWRbF7yuZ/dKmCkiMwE/gj8D1iuqvWq2tgB5LsHy51TBO1TNxLJp6qLsBT/CjuprevGcOB9VT1MVf8BXA1cZX/3GO1bL1pEp1D6InKuiBxud63A6voNFpEsVZ0FfAUcgVWorgCmiMgP7AG1o7EsnHaXz+7GvgbcZPsCJwGjsQpRW8jXS1XrsQZFZwHfAGNFZDc765fA48CdIrIrMB6r55LTjrKNtPM5+z33AWpE5DwsS3Wv1vSrJisflhW7ClgK7KeqJwNDRGS/1pKtJfKp6nbat2409/wAXgdubqu6Yct2hIgUYDUyD9vpfmCB/QcwF6te3NVW9SIddlifvv3idwIexbJGlgCFwCXAycBewIOqukhEhgJ3Areo6mcicjbWwORo4AZVnd8B5LsD+I2qfiEij2IpryzgSlVd4PETrSXfNFVdb+cZgeVDrVfV37rO/SWwm/13kaoubEfZ6lT1Vte5v8fymT+I1eX+KpOypSBf6NmJSE9V3eK6TsRxO8kX/fzOxBqcbI+60dzzy8datKw/1kB9xutGc7KJiF9VAyJyDnCKqp7hOveXWG68UbRCvcgY7e1fSuUP8Nv/RwIz7M9ZWF3Uh7D8ug8Ak4Ge9vcPAre6riEdTL6HgN/an7OBndpBvruBZ6Ly/siWe1eswu+z03M6mGwFdtrBwJkd8NnlA7l2uq8DyleIPbbQTnUjkXwjXO83q7XqRjOyPRuV52HgDPvzTq5rtEq9yOSf0yXeIbC78LcAfhF5GegBBABUtUlELgdWY42uPwpMBHYGfo/Van/oXEvtN9SB5AtgRceglo/3+3aQ70pglYgcrqr/s9P/IyK7A68C3YAjgYWq2tDRZBORI1X1g0zKlUn5CD+7jPvJMyxfe9SNRPK9Qvj9LiTDdSMV2YDtwDIRuQX4sYhMUNWVma4XrUJ7tzotaIUPx/Ip3wtchDWKPgErGmecK99lwKv2572Bl4CPgf8AhUa+ZuW7BPiv6/h07PBCoH9Xk83I17nlS0U2LNdSHdb8kDuBotZ6dq1yz+0tQAtezqHY4WT28T32izgPmGOn+bD8cU8BQ+20XsBgI1+L5HsSGOY679CuKpuRr3PLl4JsJViRPHcC+7b2s2uNvx0pemcO8KSE17V4HyhW1QexumVXqNVt3hloUtVKAFXdrFb4mZEvefkCaoWPoqrvquq7XVg2I1/nlq8lsgVVdbmqLlHVn6nqZ60sW6uwwyh9Va1RK6bZWajqGGCd/fl8YHcReRErZrbNX4aRr3PKZuTr3PK1ULY50O7LPqTNDjWQC6EYWQUGADPt5G3ADVhhmMvayHL2xMjXOWUDI1+6dGT5WiKb2j6fHZUdxtJ3EcQKaVwP7G23wjdidb3ea89CbWPk65yygZEvXTqyfB1ZtszSHgMJ6f4BB2K9pPeAC9tbHiNf15DNyNe55evIsmXyb4eckSvWsq+TsRawast1aZLCyJc6HVk2MPKlS0eWryPLlkl2SKVvMBgMhtTYEX36BoPBYEgRo/QNBoOhC2GUvsFgMHQhjNI3GAyGLoRR+oYuh4jcLCLXJPh+oojskcR1IvKJyC0icnSm5DQYWgOj9A2GWCZiLX/donyq+mu1dkIzGDosRukbugQiUiYiX4vILKwdvxCRi0TkUxH5UkSeEZECETkYa+PyP4rIFyIy3P57VUTmiMi7IjIqTr4HReQn9rUrReR3IvKhiMwWkX1F5DURWSIiF7vk+oUtw1ci8pt2eDSGLsYOt/aOwdBSxNqHdhKwD1aZ/wxr8axnVfVfdp5bsWZh3i3W5uUvqurT9ndvAher6rcicgBwj6oe5ZEv+qdXqOpBInIH1s5thwB5wHzgHyJyLNauUOOw9lSdKSKHqeo7rfYwDF0eo/QNXYFDgf+oag2ArawB9rSVfS+snaNeiz5RRLphbcH4lEup5yb5u87vzAW6qeo2YJuI1IlIL+BY++9zO183rEbAKH1Dq2GUvqGr4DX1/EFgoqp+KSLnAUd45PEBm1X1Byn8pjOVP+j67BxnYVn3v1fVf6ZwbYMhJYxP39AVeAf4kYjki0h34GQ7vTuwWkSygVJX/m32d6jqVqy9UE8Hay11ERkTnS9FXgMusHsTiMhgEemfxvUMhmYxSt/Q6VFrh6MngC+AZwBnN6YbsfYnfgNY5DrlceAXIvK5iAzHahAuFJEvsfzxp8bJ11K5XgceBT4UkbnA06TXiBgMzWIWXDMYDIYuhLH0DQaDoQthlL7BYDB0IYzSNxgMhi6EUfoGg8HQhTBK32AwGLoQRukbDAZDF8IofYPBYOhCGKVvMBgMXYj/B8C/Z72r6QHMAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ax = fm.loc['04021520', 'discharge_va'].plot(lw=0, marker='o', color='r')\n",
    "dvs['04021520']['Q'].plot(ax=ax)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "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></th>\n",
       "      <th>agency_cd</th>\n",
       "      <th>site_no</th>\n",
       "      <th>measurement_nu</th>\n",
       "      <th>measurement_dt</th>\n",
       "      <th>tz_cd</th>\n",
       "      <th>q_meas_used_fg</th>\n",
       "      <th>party_nm</th>\n",
       "      <th>site_visit_coll_agency_cd</th>\n",
       "      <th>gage_height_va</th>\n",
       "      <th>discharge_va</th>\n",
       "      <th>measured_rating_diff</th>\n",
       "      <th>gage_va_change</th>\n",
       "      <th>gage_va_time</th>\n",
       "      <th>control_type_cd</th>\n",
       "      <th>discharge_cd</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>site_no</th>\n",
       "      <th>datetime</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"11\" valign=\"top\">04021520</th>\n",
       "      <th>2005-04-27 12:11:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>1</td>\n",
       "      <td>2005-04-27 12:11:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>KGG/GRM</td>\n",
       "      <td>USGS</td>\n",
       "      <td>4.19</td>\n",
       "      <td>129.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>0.00</td>\n",
       "      <td>1.1</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-23 12:20:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>2</td>\n",
       "      <td>2005-05-23 12:20:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>GRM/WWL</td>\n",
       "      <td>USGS</td>\n",
       "      <td>4.22</td>\n",
       "      <td>135.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>-0.01</td>\n",
       "      <td>0.9</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-05-26 14:51:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>3</td>\n",
       "      <td>2005-05-26 14:51:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>DLR/GRM</td>\n",
       "      <td>USGS</td>\n",
       "      <td>5.01</td>\n",
       "      <td>199.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>0.02</td>\n",
       "      <td>1.1</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-06-07 14:45:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>4</td>\n",
       "      <td>2005-06-07 14:45:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>DLR</td>\n",
       "      <td>USGS</td>\n",
       "      <td>3.91</td>\n",
       "      <td>113.00</td>\n",
       "      <td>Good</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.7</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2005-07-22 12:50:00</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>5</td>\n",
       "      <td>2005-07-22 12:50:00</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>GRM</td>\n",
       "      <td>USGS</td>\n",
       "      <td>1.61</td>\n",
       "      <td>5.82</td>\n",
       "      <td>Fair</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.5</td>\n",
       "      <td>Clear</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020-02-19 13:08:30</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>154</td>\n",
       "      <td>2020-02-19 13:08:30</td>\n",
       "      <td>CST</td>\n",
       "      <td>Yes</td>\n",
       "      <td>JDL/RBB</td>\n",
       "      <td>USGS</td>\n",
       "      <td>3.44</td>\n",
       "      <td>24.00</td>\n",
       "      <td>Poor</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.3</td>\n",
       "      <td>IceCover</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020-04-08 08:08:47</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>155</td>\n",
       "      <td>2020-04-08 08:08:47</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>JDL</td>\n",
       "      <td>USGS</td>\n",
       "      <td>7.67</td>\n",
       "      <td>425.00</td>\n",
       "      <td>Fair</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.3</td>\n",
       "      <td>DebrisLight</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020-05-20 14:36:09</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>156</td>\n",
       "      <td>2020-05-20 14:36:09</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>JDL</td>\n",
       "      <td>USGS</td>\n",
       "      <td>3.27</td>\n",
       "      <td>37.20</td>\n",
       "      <td>Fair</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.3</td>\n",
       "      <td>DebrisModerate</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020-07-07 09:04:30</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>157</td>\n",
       "      <td>2020-07-07 09:04:30</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>JDL,GGS</td>\n",
       "      <td>USGS</td>\n",
       "      <td>2.62</td>\n",
       "      <td>1.73</td>\n",
       "      <td>Poor</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.3</td>\n",
       "      <td>DebrisHeavy</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2020-08-20 08:21:22</th>\n",
       "      <td>USGS</td>\n",
       "      <td>04021520</td>\n",
       "      <td>158</td>\n",
       "      <td>2020-08-20 08:21:22</td>\n",
       "      <td>CDT</td>\n",
       "      <td>Yes</td>\n",
       "      <td>JDL</td>\n",
       "      <td>USGS</td>\n",
       "      <td>3.42</td>\n",
       "      <td>17.90</td>\n",
       "      <td>Poor</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.3</td>\n",
       "      <td>DebrisHeavy</td>\n",
       "      <td>NONE</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>158 rows × 15 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                             agency_cd   site_no  measurement_nu  \\\n",
       "site_no  datetime                                                  \n",
       "04021520 2005-04-27 12:11:00      USGS  04021520               1   \n",
       "         2005-05-23 12:20:00      USGS  04021520               2   \n",
       "         2005-05-26 14:51:00      USGS  04021520               3   \n",
       "         2005-06-07 14:45:00      USGS  04021520               4   \n",
       "         2005-07-22 12:50:00      USGS  04021520               5   \n",
       "...                                ...       ...             ...   \n",
       "         2020-02-19 13:08:30      USGS  04021520             154   \n",
       "         2020-04-08 08:08:47      USGS  04021520             155   \n",
       "         2020-05-20 14:36:09      USGS  04021520             156   \n",
       "         2020-07-07 09:04:30      USGS  04021520             157   \n",
       "         2020-08-20 08:21:22      USGS  04021520             158   \n",
       "\n",
       "                                  measurement_dt tz_cd q_meas_used_fg  \\\n",
       "site_no  datetime                                                       \n",
       "04021520 2005-04-27 12:11:00 2005-04-27 12:11:00   CDT            Yes   \n",
       "         2005-05-23 12:20:00 2005-05-23 12:20:00   CDT            Yes   \n",
       "         2005-05-26 14:51:00 2005-05-26 14:51:00   CDT            Yes   \n",
       "         2005-06-07 14:45:00 2005-06-07 14:45:00   CDT            Yes   \n",
       "         2005-07-22 12:50:00 2005-07-22 12:50:00   CDT            Yes   \n",
       "...                                          ...   ...            ...   \n",
       "         2020-02-19 13:08:30 2020-02-19 13:08:30   CST            Yes   \n",
       "         2020-04-08 08:08:47 2020-04-08 08:08:47   CDT            Yes   \n",
       "         2020-05-20 14:36:09 2020-05-20 14:36:09   CDT            Yes   \n",
       "         2020-07-07 09:04:30 2020-07-07 09:04:30   CDT            Yes   \n",
       "         2020-08-20 08:21:22 2020-08-20 08:21:22   CDT            Yes   \n",
       "\n",
       "                             party_nm site_visit_coll_agency_cd  \\\n",
       "site_no  datetime                                                 \n",
       "04021520 2005-04-27 12:11:00  KGG/GRM                      USGS   \n",
       "         2005-05-23 12:20:00  GRM/WWL                      USGS   \n",
       "         2005-05-26 14:51:00  DLR/GRM                      USGS   \n",
       "         2005-06-07 14:45:00      DLR                      USGS   \n",
       "         2005-07-22 12:50:00      GRM                      USGS   \n",
       "...                               ...                       ...   \n",
       "         2020-02-19 13:08:30  JDL/RBB                      USGS   \n",
       "         2020-04-08 08:08:47      JDL                      USGS   \n",
       "         2020-05-20 14:36:09      JDL                      USGS   \n",
       "         2020-07-07 09:04:30  JDL,GGS                      USGS   \n",
       "         2020-08-20 08:21:22      JDL                      USGS   \n",
       "\n",
       "                              gage_height_va  discharge_va  \\\n",
       "site_no  datetime                                            \n",
       "04021520 2005-04-27 12:11:00            4.19        129.00   \n",
       "         2005-05-23 12:20:00            4.22        135.00   \n",
       "         2005-05-26 14:51:00            5.01        199.00   \n",
       "         2005-06-07 14:45:00            3.91        113.00   \n",
       "         2005-07-22 12:50:00            1.61          5.82   \n",
       "...                                      ...           ...   \n",
       "         2020-02-19 13:08:30            3.44         24.00   \n",
       "         2020-04-08 08:08:47            7.67        425.00   \n",
       "         2020-05-20 14:36:09            3.27         37.20   \n",
       "         2020-07-07 09:04:30            2.62          1.73   \n",
       "         2020-08-20 08:21:22            3.42         17.90   \n",
       "\n",
       "                             measured_rating_diff  gage_va_change  \\\n",
       "site_no  datetime                                                   \n",
       "04021520 2005-04-27 12:11:00                 Good            0.00   \n",
       "         2005-05-23 12:20:00                 Good           -0.01   \n",
       "         2005-05-26 14:51:00                 Good            0.02   \n",
       "         2005-06-07 14:45:00                 Good            0.00   \n",
       "         2005-07-22 12:50:00                 Fair            0.00   \n",
       "...                                           ...             ...   \n",
       "         2020-02-19 13:08:30                 Poor            0.00   \n",
       "         2020-04-08 08:08:47                 Fair            0.00   \n",
       "         2020-05-20 14:36:09                 Fair            0.00   \n",
       "         2020-07-07 09:04:30                 Poor            0.00   \n",
       "         2020-08-20 08:21:22                 Poor            0.00   \n",
       "\n",
       "                              gage_va_time control_type_cd discharge_cd  \n",
       "site_no  datetime                                                        \n",
       "04021520 2005-04-27 12:11:00           1.1           Clear         NONE  \n",
       "         2005-05-23 12:20:00           0.9           Clear         NONE  \n",
       "         2005-05-26 14:51:00           1.1           Clear         NONE  \n",
       "         2005-06-07 14:45:00           0.7           Clear         NONE  \n",
       "         2005-07-22 12:50:00           0.5           Clear         NONE  \n",
       "...                                    ...             ...          ...  \n",
       "         2020-02-19 13:08:30           0.3        IceCover         NONE  \n",
       "         2020-04-08 08:08:47           0.3     DebrisLight         NONE  \n",
       "         2020-05-20 14:36:09           0.3  DebrisModerate         NONE  \n",
       "         2020-07-07 09:04:30           0.3     DebrisHeavy         NONE  \n",
       "         2020-08-20 08:21:22           0.3     DebrisHeavy         NONE  \n",
       "\n",
       "[158 rows x 15 columns]"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "fm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "gis",
   "language": "python",
   "name": "gis"
  },
  "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.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}