US Endangered Species Data.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Last updated: 27 Jul 2018\n",
    "\n",
    "# Exploring US Endangered Species Data with Python\n",
    "\n",
    "This notebook retrieves data from the US Fish & Wildlife Service on endangered and threatened species listed under the Endangered Species Act. It runs a simple analysis and returns a graph of species added under the ESA over time."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "from bs4 import BeautifulSoup \n",
    "import pandas as pd\n",
    "import re, requests\n",
    "import matplotlib.pyplot as plt\n",
    "plt.rcParams['font.size'] = 16\n",
    "plt.rcParams['figure.figsize'] = [16, 8]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Data source\n",
    "\n",
    "Data are retrieved from the U.S. Fish & Wildlife Service [Environmental Conservation Online System](https://ecos.fws.gov/ecp/), which reports on Threatened & Endangered species in the United States that are protected by the Endangered Species Act (ESA). \n",
    "\n",
    "There are a number of different reports are available here: https://ecos.fws.gov/ecp/species-reports. They are mostly provided as sortable tables, viewable online."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## What question are we asking?\n",
    "\n",
    "After looking at some of the available datasets, it seems like there are data available on how many requests were made to list species, and how many species were officially listed.\n",
    "\n",
    "For the purposes of this analysis, let's ask the following question:\n",
    "\n",
    "> *How has the rate of species being listed under the ESA changed over time?*"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Downloading listing years\n",
    "\n",
    "We can use the `requests` and `BeautifulSoup` libraries to scrape the [U.S. Federal Endangered and Threatened Species by Calendar Year](https://ecos.fws.gov/ecp0/reports/species-listings-count-by-year-report) report to get the years that species were listed under the ESA."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Years when species were listed: 1967, 1970, 1972, 1973, 1975, 1976, 1977, 1978, 1979, 1980, 1981, 1982, 1983, 1984, 1985, 1986, 1987, 1988, 1989, 1990, 1991, 1992, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010, 2011, 2012, 2013, 2014, 2015, 2016, 2017, 2018\n"
     ]
    }
   ],
   "source": [
    "url = 'https://ecos.fws.gov/ecp0/reports/species-listings-count-by-year-report'\n",
    "r = requests.get(url)\n",
    "soup = BeautifulSoup(r.text, \"html5lib\")\n",
    "rows = soup.table.find_all('tr')\n",
    "listing_years = set()\n",
    "for row in rows:\n",
    "    values = row.find_all('td')\n",
    "    if len(values) > 0:\n",
    "        listing_years.add(values[0].string)\n",
    "listing_years = sorted(list(listing_years))\n",
    "print('Years when species were listed: {}'.format(\", \".join(listing_years)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "## Downloading species listed for each year\n",
    "\n",
    "We can take the list of years we determined from the last step, and again use `requests` and `BeautifulSoup` to query the [Species Listed During Calendar Year](https://ecos.fws.gov/ecp0/reports/species-listings-by-year-report?year=2018) reports for each year in `listing_years`."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "3579 species records read from Species Listed During Calendar Year\n"
     ]
    }
   ],
   "source": [
    "url = 'https://ecos.fws.gov/ecp0/reports/species-listings-by-year-report'\n",
    "all_species = []\n",
    "for y in listing_years:\n",
    "    params = {'year': y}\n",
    "    r = requests.get(url, params)\n",
    "    soup = BeautifulSoup(r.text, \"html5lib\")\n",
    "    rows = soup.table.find_all('tr')\n",
    "    for row in rows:\n",
    "        species = {}\n",
    "        data = row.find_all('td')\n",
    "        if len(data) > 0:\n",
    "            splink = data[2].a['href']\n",
    "            spcode = re.findall('spcode=(\\w{4})', splink)[0]\n",
    "            species['spcode'] = spcode\n",
    "            species['sciname'] = data[2].string\n",
    "            species['status'] = data[3].string\n",
    "            species['listyear'] = y\n",
    "            all_species.append(species)\n",
    "print(f'{len(all_species)} species records read from Species Listed During Calendar Year')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Let's take the `all_species` records and put them into a `pandas` DataFrame. There are a number of duplicate values in the DataFrame that we can also drop."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": true,
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1811 species records remaining after dropping duplicates\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>spcode</th>\n",
       "      <th>sciname</th>\n",
       "      <th>status</th>\n",
       "      <th>listyear</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>A009</td>\n",
       "      <td>Antilocapra americana sonoriensis</td>\n",
       "      <td>Endangered</td>\n",
       "      <td>1967</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>A0A8</td>\n",
       "      <td>Arctocephalus townsendi</td>\n",
       "      <td>Threatened</td>\n",
       "      <td>1967</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>A00F</td>\n",
       "      <td>Canis rufus</td>\n",
       "      <td>Endangered</td>\n",
       "      <td>1967</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>A004</td>\n",
       "      <td>Mustela nigripes</td>\n",
       "      <td>Endangered</td>\n",
       "      <td>1967</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>A000</td>\n",
       "      <td>Myotis sodalis</td>\n",
       "      <td>Endangered</td>\n",
       "      <td>1967</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "  spcode                            sciname      status listyear\n",
       "0   A009  Antilocapra americana sonoriensis  Endangered     1967\n",
       "1   A0A8            Arctocephalus townsendi  Threatened     1967\n",
       "2   A00F                        Canis rufus  Endangered     1967\n",
       "3   A004                   Mustela nigripes  Endangered     1967\n",
       "4   A000                     Myotis sodalis  Endangered     1967"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_species = pd.DataFrame(all_species)\n",
    "df_species = df_species[['spcode', 'sciname', 'status', 'listyear']]\n",
    "df_species = df_species.drop_duplicates()\n",
    "print(f'{len(df_species)} species records remaining after dropping duplicates')\n",
    "df_species.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Let's look at the numbers of species listed in each year."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "slideshow": {
     "slide_type": "-"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6gAAAHlCAYAAAD4GNGzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3Xu4ZHdZJ/rvmzQgEO5pbgkhqCCico0RhdEIKAEUcESBeDAiTjzjBcTLgDrnMOOoA2c8chln5GQICqMNAl6Iiggi4BkZAiEEkhBuAkK4tgOIihcC7/yxVuOm6U5X7aq99y/Vn8/z1LOrVlW99a69115V37V+tVZ1dwAAAGCvnbDXDQAAAEAioAIAADAIARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIawb68bSJKTTz65Tz/99L1uAwAAgB3wpje96a+6e/+xHjdEQD399NNz8cUX73UbAAAA7ICq+stFHmeILwAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQ9u11AwBwrXagjv2Yc3rn+wCADWAPKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIRwzoFbVc6vqY1V1+ZZp/6mq3l5Vb62q362qm26576er6t1V9Y6qeuBONQ4AAMBmWWQP6q8nOfuwaa9M8tXdfdck70zy00lSVXdJ8qgkXzU/579W1Ylr6xYAAICNdcyA2t1/luTjh017RXdfPd98fZJT5+sPS/LC7v7H7n5vkncnOXON/QIAALCh1vEd1O9P8kfz9VOSfGDLfVfN0wAAAOAarRRQq+pnk1yd5DcPTTrCw/oozz2vqi6uqosPHjy4ShsAAABsgG0H1Ko6N8m3Jfme7j4UQq9KcrstDzs1yYeO9PzuPr+7z+juM/bv37/dNgAAANgQ2wqoVXV2kicleWh3f3rLXRcmeVRVXa+q7pDkjknesHqbAAAAbLp9x3pAVb0gyVlJTq6qq5I8JdNRe6+X5JVVlSSv7+7/s7uvqKoXJXlbpqG/P9zdn92p5gEAANgcxwyo3f3oI0y+4Boe/wtJfmGVpgAAADj+rOMovgAAALAyARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGMIxA2pVPbeqPlZVl2+ZdvOqemVVvWv+ebN5elXVs6rq3VX11qq65042DwAAwOZYZA/qryc5+7BpT07yqu6+Y5JXzbeT5EFJ7jhfzkvyq+tpEwAAgE13zIDa3X+W5OOHTX5YkufN15+X5OFbpj+/J69PctOqus26mgUAAGBzbfc7qLfq7g8nyfzzlvP0U5J8YMvjrpqnAQAAwDVa90GS6gjT+ogPrDqvqi6uqosPHjy45jYAAAC4ttluQP3ooaG788+PzdOvSnK7LY87NcmHjlSgu8/v7jO6+4z9+/dvsw0AAAA2xXYD6oVJzp2vn5vkpVumf+98NN97J/nrQ0OBAQAA4JrsO9YDquoFSc5KcnJVXZXkKUmemuRFVfW4JO9P8l3zw1+W5MFJ3p3k00keuwM9AwAAsIGOGVC7+9FHuev+R3hsJ/nhVZsCAADg+LPugyQBAADAtgioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIezb6wYAgDU7UMd+zDm9830AwJLsQQUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYgoAKAADAEPbtdQMAAADbcqCO/Zhzeuf7YG3sQQUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCHs2+sGAGAhB2qxx53TO9sHALBj7EEFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIbgKL4A7KxFjr7ryLsAQOxBBQAAYBACKgAAAENYaYhvVT0xyQ8k6SSXJXlsktskeWGSmye5JMljuvufli5uSBgAAMBxZdt7UKvqlCSPT3JGd391khOTPCrJ05I8vbvvmOQTSR63jkYBAADYbKsO8d2X5PpVtS/JDZJ8OMn9krxkvv95SR6+4msAAABwHNh2QO3uDyb5pSTvzxRM/zrJm5J8sruvnh92VZJTjvT8qjqvqi6uqosPHjy43TYAAADYEKsM8b1ZkocluUOS2ya5YZIHHeGhR/yiaHef391ndPcZ+/fv324bAAAAbIhVhvg+IMl7u/tgd38mye8k+YYkN52H/CbJqUk+tGKPAAAAHAdWCajvT3LvqrpBVVWS+yd5W5JXJ3nE/Jhzk7x0tRYBAAA4Hmz7NDPdfVFVvSTTqWSuTvLmJOcn+cMkL6yqn5+nXbCORgGAXbbIKd8Sp30DYG1WOg9qdz8lyVMOm/yeJGeuUhcAAIDjz6qnmQEAAIC1EFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYwr69bgAAdt2BOvZjzumd7wMArslx+H5lDyoAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhrBvrxsAAADYcwfq2I85p3e+j+OcPagAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQ9i31w0AXCsdqMUed07vbB8AABvEHlQAAACGIKACAAAwhJUCalXdtKpeUlVvr6orq+rrq+rmVfXKqnrX/PNm62oWAACAzbXqHtRnJnl5d985yd2SXJnkyUle1d13TPKq+TYAAABco20H1Kq6cZJvTHJBknT3P3X3J5M8LMnz5oc9L8nDV20SAACAzbfKHtQvTXIwya9V1Zur6jlVdcMkt+ruDyfJ/POWa+gTAACADbfKaWb2Jblnkh/t7ouq6plZYjhvVZ2X5LwkOe2001Zo4xicCgIA9t4i78feiwGOe6vsQb0qyVXdfdF8+yWZAutHq+o2STL//NiRntzd53f3Gd19xv79+1doAwAAgE2w7YDa3R9J8oGq+op50v2TvC3JhUnOnaedm+SlK3UIAADAcWGVIb5J8qNJfrOqrpvkPUkemyn0vqiqHpfk/Um+a8XXAI53huoDABwXVgqo3X1pkjOOcNf9V6kLAADA8WfV86ACAADAWgioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQ9u11AwAAwA45UIs97pze2T5gQfagAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYwr69bgDguHegjv2Yc3rn+wAA2GP2oAIAADAEARUAAIAhCKgAAAAMwXdQAQCA3bPIsRcSx184TtmDCgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGMK+vW4ArtEiJ3J2EmcAANgI9qACAAAwBAEVAACAIQioAAAADEFABQAAYAgCKgAAAEMQUAEAABiCgAoAAMAQBFQAAACGIKACAAAwhH173QAAwJ44UIs97pze2T4A+Dx7UAEAABjCygG1qk6sqjdX1R/Mt+9QVRdV1buq6req6rqrtwkAAMCmW8cQ3yckuTLJjefbT0vy9O5+YVU9O8njkvzqGl4HgN1i6CMAsAdW2oNaVacmeUiS58y3K8n9krxkfsjzkjx8ldcAAADg+LDqEN9nJPk3ST43375Fkk9299Xz7auSnLLiawAAAHAc2PYQ36r6tiQf6+43VdVZhyYf4aFHHP9VVeclOS9JTjvttO22AYtbZMii4YoAALBnVtmDep8kD62q9yV5Yaahvc9IctOqOhR8T03yoSM9ubvP7+4zuvuM/fv3r9AGAAAAm2DbAbW7f7q7T+3u05M8Ksmfdvf3JHl1kkfMDzs3yUtX7hIAAICNtxPnQX1Skh+vqndn+k7qBTvwGgAAAGyYdZxmJt39miSvma+/J8mZ66gLAADA8WMn9qACAADA0gRUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDWMt5UNkQB+rYjzmnd74P2EmWcwBgJy3yWSPxeeMo7EEFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCPv2ugFggx2oYz/mnN75PgAAuFawBxUAAIAhCKgAAAAMwRBfgE1iWDUAcC1mDyoAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABD2LfXDQAALOVAHfsx5/TO9wHA2tmDCgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQxBQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhCKgAAAAMQUAFAABgCAIqAAAAQ9h2QK2q21XVq6vqyqq6oqqeME+/eVW9sqreNf+82fraBQAAYFPtW+G5Vyf5ie6+pKpulORNVfXKJN+X5FXd/dSqenKSJyd50uqtcq1xoBZ73Dm9s30AAADXKtveg9rdH+7uS+brf5PkyiSnJHlYkufND3tekoev2iQAAACbby3fQa2q05PcI8lFSW7V3R9OphCb5JbreA0AAAA228oBtapOSvLbSX6suz+1xPPOq6qLq+rigwcPrtoGAAAA13IrBdSquk6mcPqb3f078+SPVtVt5vtvk+RjR3pud5/f3Wd09xn79+9fpQ0AAAA2wCpH8a0kFyS5srt/ectdFyY5d75+bpKXbr89AAAAjherHMX3Pkkek+Syqrp0nvYzSZ6a5EVV9bgk70/yXau1CAAA7LlFztTgLA2saNsBtbv/R5KjLaX3325dAAAAjk9rOYovAAAArGqVIb4AwLosMnQuMXwOgI1mDyoAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCE4zQwsy6kgAGDvLfJ+7L0YrnXsQQUAAGAIAioAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCEIqAAAAAxBQAUAAGAIAioAAABDEFABAAAYwr69bgAAgNmBWuxx5/TO9nG4Rfra7Z6AjWQPKgAAAEMQUAEAABiCgAoAAMAQfAf12m7U76qwuywHAAAczbXos6I9qAAAAAxBQAUAAGAIhvgCAKzKaVgWcy0aZgjsDXtQAQAAGIKACgAAwBAEVAAAAIYgoAIAADAEARUAAIAhOIrvXnG0PwAAgC9gDyoAAABDEFABAAAYgiG+wBcy/BwAgD1iDyoAAABDEFABAAAYgiG+AAAwGl+54ThlDyoAAABDEFABAAAYgoAKAADAEARUAAAAhiCgAgAAMAQBFQAAgCE4zQzsJYeQBwDg2mSHP7/agwoAAMAQBFQAAACGYIjvMgzHBAAA2DH2oAIAADAEARUAAIAhGOILALCJfDVp9/mdw8rsQQUAAGAIAioAAABD2LEhvlV1dpJnJjkxyXO6+6k79VoAALAtiwzLTQzNhV2yI3tQq+rEJP8lyYOS3CXJo6vqLjvxWgAAAGyGnRrie2aSd3f3e7r7n5K8MMnDdui1AAAA2ADVvf7hClX1iCRnd/cPzLcfk+TruvtHtjzmvCTnzTe/Isk7Fih9cpK/WkOL66qzzloj9rTOWnra/Vp62t06o9bS0+7X0tPu19LT7tfS0+7X0tPu19LTemvdvrv3H6vQTn0H9UiD+b8gCXf3+UnOX6po1cXdfcYqja2zzqb3tM5aetr9Wnq69va0zlp62v1aetr9Wnra/Vp62v1aetr9Wnram1o7NcT3qiS323L71CQf2qHXAgAAYAPsVEB9Y5I7VtUdquq6SR6V5MIdei0AAAA2wI4M8e3uq6vqR5L8cabTzDy3u69YQ+mlhgTvQp111hqxp3XW0tPu19LT7tYZtZaedr+Wnna/lp52v5aedr+Wnna/lp72oNaOHCQJAAAAlrVTQ3wBAABgKQIqAAAAQxBQAQAAGIKACgAAwBAEVADguFVVN6mqR1bVj1fVE+frN11j/W/ZxnNuXFVfdoTpd12yzq2r6tbz9f1V9S+r6quW7ecotX9xTXXuMPd15yWfd1pVfcl8varqsVX1n6vqX1fVUmepqKqHHqq1qqr6xqr6ivn6favqJ6vqIdusdVJVPWJeLn+0qs6uqm19dq+qO1fVk6rqWVX1zPn6V26n1lHqP3abPd2/qk46bPrZS9Y5s6q+dr5+l/l/+cHL9nOU2s9fU537zn1965LP+7qquvF8/fpV9e+r6ver6mlVdZMlaz2+qm63zHOOUue6VfW9VfWA+fY5VfUrVfXDVXWdVesnAx/Ft6p+Oclvd/efr6HWA5M8PMkpSTrJh5K8tLtfvs1apyZ5VXe/b8v07+/u5y5Yo5J819zLS5LcL8nDkrw9ybO7+3NL9HNyd//Vltv/R5Izk1ye5L/1En/gqvqOJK/t7o9X1f4k/2+SeyR5W5Kf6O6rFqyzzr/dzZP8SKa/2QVJfibJ1ye5Mskvdvcnlqj1zUm+M8ntklyd5F1JntPd795GX3fO9Dfbukxd2N1XLlvrGl7jsd39a0v2dEqSi7r7b7dMP3vZZb2qzkzS3f3GqrpLkrOTvL27X7ZMnSPUfX53f+8qNeY69828nHf3K5Z87tclubK7P1VV10/y5CT3zLSc/2J3//WCdR6f5He7+wPLdX/EWofOF/2h7v6TqjonyTdkWs7P7+7PLFHry5J8R75wOX/BovN1WK21rTuv4TX+7+7+uSV7WmkdPD9+bevhI9T+0+6+3zaeZ32+eK21rM+r6nuTPCXJK5J8cJ58apJvSfLvu3vlD6dV9f7uPm2Jx393kmck+ViS6yT5vu5+43zfJd19zwXr/GCm9VsleVqS70tyRZL7JPl/uvuCJXp61uGTkjwmyfOTpLsfv0St3+vuh8/XH5ZpXl+TaZ33H7v71xesc3mSM7v701X1tCRfluT3Mv0vp7u/f4me/j7J3yX5oyQvSPLH3f3ZRZ+/pc4zMv3P7st0qsX7zzW/Kcmbu/unlqj13Ul+Kslbknxzktdl2rH0NUm+p7svW6LWk5I8OskLkxz63z810/vOC7v7qYvWuobXWHY5f3ySH870/3/3JE/o7pfO9y2znD8lyYMy/c5fmeTrMi1PD8j0d/yFJXq68PBJmX73f5ok3f3QJWq9obvPnK//q0zz+rtJvjXJ7y/6O6+qK5LcbT6F5/lJPp3pPev+8/R/uURPf51pOf+LTMv5i7v74KLP31LnNzP9vm+Q5JNJTkryO3NP1d3nLlvzi15j4IB6MMlfJtmf5LcyfcB68zbqPCPJnTKtRLf+U35vknd19xOWqPWLSe6b5JIk357kGd39n+f7lvln+q9Jbpnkukk+leR6SX4/yYOTfHTJnj7/ulX1b5P8iyQHknxbkqu6+4lL1Hpbd99lvv5bSV6f5MWZ/sm/p7sX2gq8rr/dXOtlSS5LcuMkXzlff1GmDw936+6HLVjnqUluleRVmT5wvzfJO5P8UKYPRi9eoqcdX9HPr7Pwyn5dK/r58WtZ2Y+4op+fv5aV/bpW9HOttazs5+Xg25O8NtP65NIkn8gUWH+ou1+zRE9rW3ce43WWWc7Xsg6eH7+W9XBVvfXwSZl+b+9Iku5eeI+X9fmerM/fkeTruvuTh02/WaaNfXdasM7h67vP35Xkft19wyV6ujTJg7r7w/PGwucn+Znu/p2qenN332PBOpdlWn9fP9Pf8Mu7+yPzvL26u+++RE9XZXofeMU8T0nyS0l+Mkm6+3lL1Pr8PFTV6zItj++tqpMzbXi624J1ti7jb0rytYc2LFXVWxatc6inTMH2EZney78603vMC7r7tUvUuWJ+7vUzbfA4ZQ7Q18kUUL96iVpvTXLv+fknJ/nN7n5gTXvRn93d37BErXcm+arDN3bOG0ev6O47LtHTEe9Kcqfuvt4SPV2W5Ou7+2+r6vRM78P/vbufuY3l/O6Z1uEfSXLqlg3QFy27Ds60Ee85mTZeVqb390clyZLLwtbl/I1JHtzdB6vqhkle391fs2CdK7v7Kw/1t/V9rqouXfL/+M1J7pXpfeCRSR6a5E2Z5vF3uvtvFqzz1u6+a00jFT6Y5Lbd/dl5w+9blvmdH1V3D3nJ9I+cJHdM8n9l2ur39kxbOu+0RJ13HmV6ZfqQtUxPlyXZN1+/aZKXJXn61n4XrTP/vE6S/5XkuvPtfYfuW/b3NF+/JMkNt9RettY7tlx/02H3Xbrbf7utrzv/vT64Qk+Xbbm+L8mfz9dvlmkv3DI9vTPJdY4w/brbWKbeepTLZUn+ccll86T5+ulJLs4UUpdaNrfUOjFTWPpUkhvP06+f5K1L1LkkyW8kOSvT1uOzknx4vv5NS/a0dTl/Y5L98/UbbmM5v3JrjyssU2/OtDX7WzPtDTqY5OVJzk1yo2WXgy3L5keTnDjfriV/55dtee4Nkrxmvn7aNpaDda47P3WUy98kuXrJ+Vt5HXyo1vxzpfVwkgvn5fzOSW4///99YL5++xWWc+vzY/zttvy9Vl2f3+QI02+yzHKeaUPQQzKv37Zczsq0wWPpZXPL7dtk+hD5+By2zjpGnUu2XH/L0Za1BWvdKNOezgOZQleSvGeZGkfp6w3b7SvTHsr7zdd/+9D/W5JbHD6/y/Q03771/Pv+n0k+sESdy+efXzIvE9efb5+Y5G3LLgf55x1J1z9s/bDscv72I62P5vXUO5ao89FMYfD2h11OzzQCaJme3nbY7ZMyvYf+8nbWUUdafpapMz/+hCRPzLRx/u7ztO0u52+Z10m3SHLx0XpeoM6Lkzx2vv5rSc6Yr98pyRuX7Onw5fw6mULqC5IcXKLO5Zk+894s0/v4zefpX5Itn7FWuSw1Rn+XdZJ097uS/Ick/2HeavToTB9KvnzBOv9QVWd29xsOm/61Sf5hyZ72dffVc1+frKpvT3J+Vb040x9qUYdqfKaq3tjd/zTfvrqqlh1Scv2qukemf6oTu/vvttRettZrqurnkvzH+frDu/v35qFUywwPXNffLklOmLf23ijJSVV1ene/r6pukeV+55+rqpt398eT3DbTm0W6+xPzFp9lfG6u8ZeHTb/NfN8ybpXkgZneyLaqTMN5FnViz8N659/PWUleUlW3zz9v7V7U1T0Nbfp0Vf1Fd39qrvv3VbXM/J2R5AlJfjbJT3X3pVX1973EFsgtDi0HJ2R6wz449/R3VXX1krUur38ePv2Wqjqjuy+uqjslWXgo7fTy/blMexReMW8hf1Cm5fyXMu1xWtQJ85bsG2YKljdJ8vFMW4SX/T7HviSfnZ97o7nR99fy3wtZ57rzk5n2bnz08Duqapkh0utaBydrWg9390NrGk57fpJf6u4Lq+oz3X34+mER1ueLWef6/BeSXFJVr8i0YSGZNuh8S6b5XdTrk3z6SOu3eS/tMv6mqr6su/8iSXrak3pWpuGry3x/9HNVdZ2e9pp9/juQNX3XcqnvMfa0Z+XHqupeSX6jqv5w2Rpb3K2qPpXpvel6VXXrnvbsXjfz33JBP5Dk+VX17zIt05fOe4huluTHl+zpC5ab7v5Ikmcledb8PrqoP6yq/z/TB/XnJHlRVb0+08aKP1uyp5cleXlVvTbTe8uLk88PlV92Of+xJK+qqnflC5fzL8807H5Rf5BpY/ilh99RVa9ZsqePVNXdD9XqaU/qtyV5bqZhzIv6p6q6QXd/OtPewUP93CTrZN/KAAAK60lEQVRLfiab39OfPr+nPL2qPppsOyvdJNOGpUrSW5bzk7Lc3+8HkjxzHlXzV0n+5/y++YH5vmUcvpx/JtNG1gvnPc6LuiDTRo8TM33Ge3FVvSfJvTONLlzdOlLuTlyy5Na9a6hzryQXZdpl/4r5cuU87V5L1vqDHGHPT5KfT/K5Jer8Uea9XYdNv3UO25q4QK3XJHn1lstt5ulftMVmgVrXSfLvkrx/vnwu05aRA0lO2+2/3Vzr0Zm22H000/eN/iTTlq0PJjlviTqPzBQoXzHP20Pm6fuTHFiyp7OTvHv+O54/X14+Tzt7yVoXJLnvUe5buK9Mw2bvfti0fZmGhn12yZ4uSnKD+foJW6bfJEtsvd/yvFMzvbH+SpL3b3M5eF+S92QayveeJLeep5+U5beQ3iTJr2camntRplD6nkzDYu+2RJ2jLueZt5ovUeuJcw9/mWmr/auS/LdMW9CfskSdJ2TaA39+pjePQ1td9yf5syV7umfWt+78+UzfFTvSfU9bos5a1sHzc9a2Hp6fd8NMW/4vzDQcdzvL+atjfb5InbWtz+fn3SzTEL6fyDRk9VFJbrau+d5GP3dLcsej/E2/Z4k6p+XIo31OSfKAFfqrTF+z+I01z/dNMw35XPZ5X5np++PfmWlI8wnbqHHWGufj6zMNzU2m78X+ZJLv3mZfD56f/y1bpp2Q5HrbqHVCpgDxnZmGMt8784ibvbhk+mxw66Pcd58l6hzxd5Hk5CRfs2KPD8n0tYF1zvcNktxhG8+70bxuuFeSW23ztZca9XKMWrfNNLT30P/uI3KU9/ntXEb+DupJveVAL2uod+tMK+XK9OHhI9uocf1k2pN0hPtO6e4PfvGzlqp/w0xDuj62Sp251omZ/mk/vc3n3yTT3or/tY3nrvtvd2KmvWZXz+Pd755peNiHl6xz8yRfmuTdfdj3jbbR0wmZDoTw+WUq01CLpQ+qsA5VdWqmPZ9ftFxX1X16iQOcVNX1uvsfjzD95EwfmBc+MMNhz39Ipjedn9nO849S8waZVtTv3cZzb5RpediXaZ3wRXv3jvH8O3X3O5d93Wuod9sk6e4P1XQE0QdkCvSH78E8Vp2vyvSB7fLufvsa+lp53bkuO70OnuustB6uqrtl+pD97FV72VLT+vyL66xtfT7Xu1W2HAxs2fXBuutsek/rrKWn1azrf3yd64pN7mmdtTa5p5ED6nWTfKbnBudhSffMNGb9j5aoc9fuPtqXupftaS21RuxpnbU2vae53mlJPtXTMMPTMw1nvbK7r1hTrbd39+V7UWfTe1pnrU3vaa53RrYcKXWV0LuuWnra/Vqb2lNV3T3JszONrLgq04aYUzMNS/+h7r5kwTr3SPKrc52tRwNeqs4Rejq81r/uBQ9Sta46C9Radv6u6Xe1rvlb5+98meVguJ4WeJ2ljr6703U2vad11trknkYOqG/JNOTiE1X1U5mOQvmyTOP4L+7un16wzmczDQt8Qaajsb1thZ7WUmvEntZZ6zjo6clJfjDJP+afj2L455mGy1zQ3b+827X0ZP52oKdvynRqkk9mGlL055mGQn4myWN6idPrrKuWnszfDvR0aZIf7O6LDpt+7yT/Xy9+RNm11Nn0ntZZS09L1Trad3Iryc929813s86m97TOWpve01FfY+CAennPh+OuqouT/IueDtKyL9P34BY6hHFNX5h/TKbvvjwy02khXpDpdCDvW7KntdQasad11joOeroi016pG2T6buSX9j8fOvyiXu4w8muppSfztwM9vTnJt87Pv0OSX+7u76iqb8l00KuFTza+rlp6Mn870NO7+iin2Kiqd3f3QgeBWledTe9pnbX0tFStf0jynzIfHO4wT+zum+5mnU3vaZ21Nr2no+o1fvF3nZdMRzD96vn6yzMfsCDTkdEWPrx2vviQymdmOpDFB5K8bsme1lJrxJ42ff7W3NOhU4KcmOlE6lsPJLTsod/XUktP5m+netpSb+upIa7Yi1p6Mn870NOzkvxhpg2X3zBfHjlP+5XdrrPpPW36/I3Y01zrdTnKwe2y3Gl01lJn03va9PlbZ09Hu4y8B/WuSf57pvMIJcl9Mh1l866ZtpYeWLDOEU/2W1WV5Bt7myfdXaXWiD2ts9Zx0NOvZzolwg2TfDrTFqSXZzrJ9426+7t3u5aezN8O9PTcTAfleFWmI2R+sLt/vKYDU13S3Xfe7Vp6Mn/r7mmu96C5ztaD3l3Y3S/bizqb3tM6a+lp4TpfkeTjPZ+i7bD7btULHnhpXXU2vad11tr0no76GqMG1CSp6Wh/35rpZLT7Mv1j/nEvccS+qjpn0TC7W7VG7GmdtY6DnvYl+a5MH5Bekumw9o/OdLqD/9LzuQt3s5aezN8O9HSdJP8qyV0ybSh8bnd/tqYj6d6ylzjP57pq6cn8rbsnAMYzdEAFANgpNZ2C56cz7aW65Tz5Y0lemuSpi24QX1edTe9pnbX0tK1aD890vuBt1VpXnU3vaZ21Nr2nozlh1QI7papOqqqfq6orquqvq+pgVb2+qr5vL+psek/rrHUc93TuXtXS0+7XOo56unyN/zMr1dLT7tfa9J6SvCjJJ5J8c3ffortvkeSbMx0h+MV7UGfTe1pnLT0tX+usw2p9Ysla66qz6T2ts9am93REw+5BraqXJvndJH+S5LszfafqhUn+babvm/zMbtbZ9J7WWUtP5m/UntZZS0/mb9Se1lnrOOjpHd39Fcvet1N1Nr2nddbSk/kbtad11tr0no6q13CkpZ24JHnLYbffOP88IdNJ53e1zqb3tOnzN2JPmz5/I/a06fM3Yk+bPn8j9rTp87fmnl6R5N8kudWWabdK8qQkf7LbdTa9p02fvxF72vT5G7GnTZ+/dfZ0tMuwQ3yT/F1V3TdJqurbk3w8Sbr7c0lqD+psek/rrKWn3a+lp92vpafdr6Wn3a+16T09Msktkry2qj5RVR9P8pokN8+0d3a362x6T+uspafdr6Wn3a+16T0d2TpS7k5cMp1O5g2Zxtr/jyR3mqfvT/L43a6z6T1t+vyN2NOmz9+IPW36/I3Y06bP34g9bfr8rbOn+Xl3TvKAJCcdNv3svaiz6T1t+vyN2NOmz9+IPW36/K2zpyPWX0eR3b4keexIdTa9p02fvxF72vT5G7GnTZ+/EXva9PkbsadNn79l6yR5fJJ3JPm9JO9L8rAt912y23U2vadNn78Re9r0+Ruxp02fv3X2dNTXWEeR3b4kef9IdTa9p02fvxF72vT5G7GnTZ+/EXva9PkbsadNn79l6yS5LPMegCSnJ7k4yRPm22/e7Tqb3tOmz9+IPW36/I3Y06bP3zp7OtplXwZVVW892l2Zvoi7q3U2vad11tLT7tfS0+7X0tPu19LT7tfa9J6SnNjdf5sk3f2+qjoryUuq6vZZ7vus66qz6T2ts5aedr+Wnna/1qb3dETDBtRMbzIPzHROna0qyev2oM6m97TOWnra/Vp62v1aetr9Wnra/Vqb3tNHquru3X1pknT331bVtyV5bpKv2YM6m97TOmvpafdr6Wn3a216T0e2jt2wO3FJckGS+x7lvgO7XWfTe9r0+Ruxp02fvxF72vT5G7GnTZ+/EXva9Plbc0+nJrn1Ue67z27X2fSeNn3+Ruxp0+dvxJ42ff7W2dPRLjUXAwAAgD018nlQAQAAOI4IqAAAAAxBQAUAAGAIAioAAABDEFABAAAYwv8G6Mww9xAy2koAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1152x576 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "year_counts = df_species['listyear'].value_counts().sort_index()\n",
    "plt.rcParams['figure.figsize'] = [16, 8]\n",
    "year_counts.plot(kind='bar', color='orange');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "source": [
    "Let's also look at the percentages of each current conservation status for all species"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "slideshow": {
     "slide_type": "slide"
    }
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtkAAAHICAYAAACbL/Y9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3XmYZGV99vHv79TS+zr7DDPUMOygEIQIKohL3Mo9KEaNg1tMNCFu0TJxj0u5RF831BgF9xAMcatEcEcUREGgUREGKGafnr33rep5/zine7p7uqd7Zqr7qeX+XFddU+fUqeq7unuq73rqOeeYcw4RERERESmdwHcAEREREZFqo5ItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmIq2SIiIiIiJaaSLSIiIiJSYirZIiIiIiIlppItIiIiIlJiKtkiIiIiIiWmki0iIiIiUmJx3wFEpLKkMjkDLFq0aZfp64bz2XRx0UOKiIh4Zs453xlEZJGlMrkOYCmwBOiM/p3ressxfKki0APsjy77Jl2f6bIPeDifTe89xqcmIiJSFlSyRapQKpOLAeuBU4CTJl3WR5dWf+nm5SDwwCyXrRodFxGRcqeSLVLBUplcHDgbeARw+qTLyUDSY7SFNAzkCQv3JuBe4LfAXflsesRjLhERkQkq2SIVJJXJrQUeDVwY/fsooMFrqPIxAtwN3Ab8Jrr8UaPeIiLig0q2SJlKZXJNwAWEZXr8stprqMrTC9zBpOKdz6bzXhOJiEhNUMkWKROpTK4NeArwJOAi4Cwg5jVUdeoGfgL8H3BDPpve5TmPiIhUIZVsEY9SmdxZwDOANPBYdFjNxeaA3wE/iC635LPpMb+RRESkGqhkiyyiVCbXQDhS/YzocqLfRDLNQeBHRKU7n01v9ZxHREQqlEq2yAJLZXLrCUeq08ClQL3XQHI07iEs3P+dz6Zv9R1GREQqh0q2yAJIZXLLgRcDfw2c5zmOlMYm4OvA1/LZ9CbfYUREpLypZIuUSCqTqwOeDbwMeBqaX13Nfg18Fbg2n03v8R1GRETKj0q2yHFKZXKPISzWLwQ6PMeRxTUK3AB8DfhOPpse8pxHRETKhEq2yDFIZXIpwqkgLyM8u6JID/DfhIX7p/lsWi+uIiI1TCVbZJ5SmVwM+EvgtcAlgPlNJGVsE/AZ4Ev5bLrHdxgREVl8Ktkic0hlcs3Aq4B/BFJ+00iF6QO+Anwqn03f6zuMiIgsHpVskVmkMrk1wJXA3wDtnuNIZXOEx9/+WD6b/oHvMCIisvBUskWmSWVy5wJvAi4HEp7jSPXpAv4N+EY+mx71HUZERBaGSrYIkMrkjPCwe28Gnug5jtSGbcAngM9r3raISPVRyZaaFu3M+FLgLcCZnuNIbTpIOLL9sXw23e87jIiIlIZKttSsVCb3XOADwBm+s4gA3cD7gc/ls+kR32FEROT4qGRLzUllchcDWeAxvrOIzCAPvBv4aj6bLvqNIiIix0olW2pGKpM7G/gg8EzfWUTm4ffA2/PZ9Ld9BxERkaOnki1VL5XJrQPeS3iGxsBzHJGjdSvwtnw2/TPfQUREZP5UsqVqpTK5JcA/A68D6jzHETleNxKW7Tt8BxERkbmpZEvVSWVySeCNQAZo8xxHpJQccA3w5nw2vc9zFhEROQJ9dL4AzKxgZndOumSO8v55M1u6UPlKwcwuNbPv+84xXSqTexxwJ+HcaxVsqTYGvBy4N5XJvcx3GBERmZ1GsheAmfU555qP4/554Hzn3J7SpZrX1zXC34k5j2hgZpcCb3bOlcVOhKlMrg34MPBqwiIiUgt+CvxtPpu+z3cQERGZSiPZiygaoX6Pmd1hZl1mdnq0fomZ3WhmvzOzzzOpJJrZt83sdjP7vZn9zaT1fWb2fjO7y8xuNbMV0foN0fJvzOy9ZtY36T7/FK2/28zeE61Lmdkfzewq4A5grZk9xcxuiXJeZ2bN0bZPM7N7zexm4PmL8T2bj1Qm9wLgj8DfoIItteUJwN2pTO5dqUxO+x2IiJQRleyF0TBtusjlk27b45w7D/gs4Sm8Ad4F3Oyc+zPgu8C6Sdu/wjn3KOB84EozWxKtbwJudc6dA9xEOIIL4WmaP+GcuwDYPv4gZvYU4BTgz4FzgUeZ2SXRzacBX4m+fj/wduDJUc7fAm80s3rgC8CzgIuBlcfzDSqFVCa3NpXJfRf4L2CV7zwintQRHlf7rlQmd6nfKCIiMi7uO0CVGnTOnTvLbddH/97OodHgS8avO+dyZrZ/0vZXmtnzoutrCYvyXmAE+P6kx/qL6PpFwHOj698APhpdf0p0+V203Bw91mbgYefcrdH6CwlPL/7LcPYISeAW4HTgIefc/QBm9jXCkeNFl8rkAuBK4F8Jn4eIhG+Wf5rK5L5MuGPkok43ExGRqTSSvfiGo38LTH2Tc9jk+Gje85OBi6IR698B9dHNo+7QhPrpjzUTAz7onDs3upzsnPtidFv/tO1+OGm7M51zr5wt42JLZXLnEB43+OOoYIvMZCPhjpFX+A4iIlLLVLLLw03ASwDM7OlAR7S+DdjvnBuI5m9fOI/HuhX4y+j6iyatvwF4xaT51WvMbPks93+smZ0cbddoZqcC9wLrzWxDtN1fzfvZlUAqk4unMrn3EU5fuWAxv7ZIBVoCXJ3K5L6VyuQ65txaRERKTiV7YUyfk52dY/v3AJeY2R2EUzo2R+t/AMTN7G7CqRG3znL/yV5POIf6NsJ5ygcBnHM3Ek4fucXMuoBvAS3T7+yc2w1cAXwz+rq3Aqc754YIp4fkoh0fH55HlpKIztj4c+Bf0BQnkaPxl8Cd0aEtRURkEekQflXGzBoJ54Q7M3sR8FfOuef4znWsUpnc84Avcmh0X0SOXoHwjfr78tl0wXcYEZFaoJJdZczsYuDThHOrDxAenWST31RHLzoc2b8RnhJdRErjJuCl+Wx6i+8gIiLVTiVbyk4qkzsNuBY4x3cWkSq0D3hVPpv+H99BRESqmUq2lJVUJrcR+AzhccBFZOF8DnhDPpse8h1ERKQaqWRLWUhlcs3AVcBf+84iUkPuAf4qn03f4zuIiEi1UckW71KZ3LmE00NO9Z1FpAYNAq/MZ9Pf9B1ERKSaqGSLV6lM7nLgaqDBdxaRGvde4N35bFp/FERESkAlW7xIZXJGeHzwd/jOIiIT/gu4Ip9ND/oOIiJS6VSyZdGlMrlG4MvAZb6ziMhhbgOek8+md/oOIiJSyVSyZVGlMrk1wHeB83xnEZFZbQGenc+m7/QdRESkUum06rJooh0cb0MFW6TcrQVuTmVyFXu2WBER31SyZVGkMrmnEZ5tbrXvLCIyL03A9alM7i2+g4iIVCJNF5EFl8rkXgV8Foj7ziIix+Rq4DX5bHrUdxARkUqhki0LKpXJvQ/4F985ROS4/ZRwnnaf7yAiIpVAJVsWRCqTC4DPA6/ynUVESuZXwNPz2XSP7yAiIuVOJVtKLirYXwSu8BxFRErvN8BT89n0ft9BRETKmUq2lFRUsK8B/tpzFBFZOHcCf5HPpvf4DiIiUq5UsqVkUplcDPgK8GLfWURkwd0DPDmfTe/yHUREpBypZEtJpDK5OPB14IW+s4jIovkT8MR8Nr3ddxARkXKj42TLcUtlcgngP1HBFqk1pwE3pTK5db6DiIiUG41ky3GJCvZ/Ac/1nUVEvHmYcET7Qd9BRETKhUay5ZilMrkk8N+oYIvUuhOBn6cyuVN9BxERKRcayZZjksrk6oDrgWf4ziIiZWMb8Jh8Nr3ZdxAREd80ki1HLTpM3zdQwRaRqdYA/5fK5Dp8BxER8U0lW47Fx4Dn+w4hImXpTOC7qUyu3ncQERGfVLLlqKQyuTcA/+g7h4iUtccBX48+9RIRqUl6AZR5S2VyLwD+zXcOEakIzwc+6TuEiIgv2vFR5iWVyT0O+BFQ5zuLiFSUf85n0x/0HUJEZLGpZMucUpnc6cAvgU7fWUSkIl2Rz6a/7DuEiMhiUsmWI0plciuAW4GU5ygiUrnGgGfms+kbfAcREVksKtkyq1Qm1wT8HHiU7ywiUvH6gEvz2fTtvoOIiCwG7fgoM0plcjHgWlSwRaQ0moH/TWVya30HERFZDCrZMptPAmnfIUSkqiwHrktlcknfQUREFppKthwmlcm9BHit7xwiUpUeTXhCKxGRqqY52TJFKpM7A/gN0OQ7i4hUtZfks+lv+A4hIrJQVLJlQiqTayQs2Gf6ziIiVa8f+PN8Nv0H30FERBaCpovIZJ9FBVtEFkcTcH0qk2vxHUREZCGoZAsAqUzuVcDLfOcQkZpyGvBF3yFERBaCposIqUzuHMITztT7ziIiNekN+Wz6//kOISJSSirZNS6VybUCvwVO8Z1FRGrWKOGJan7lO4iISKlouoh8ERVsEfErAfxXKpNb7juIiEipqGTXsFQmdyVwme8cIiLAGuDrqUzOfAcRESkFlewalcrkzgM+4juHiMgkTwZe5zuEiEgpaE52DUplcgnC42Gf4zuLiMg0A8A5+Wx6k+8gIiLHQyPZtemtqGCLSHlqBL6cyuT090lEKppexGpMKpM7E3iH7xwiIkfwGOBNvkOIiBwPTRepIdHI0C+BC31nERGZwzBwnk67LiKVSiPZteVKVLBFpDLUAV/UtBERqVR68aoRqUzuJOD9vnOIiByFC4G/9x1CRORYqGTXji8Q7lAkIlJJPpDK5E70HUJE5GipZNeAVCb3auCJvnOIiByDJuDzvkOIiBwt7fhY5VKZ3Brg90Cb7ywiIsfhZfls+qu+Q4iIzJdGsqvfZ1HBFpHK95FUJtfiO4SIyHypZFexVCaXBp7lO4eISAmsADK+Q4iIzJemi1SpVCYXA+4GzvSdRUSkRIaA0/LZ9GbfQURE5qKR7Or1SlSwRaS61AMf9B1CRGQ+NJJdhVKZXDNwP7DSdxYRkRJzwIX5bPo230FERI5EI9nV6Z9QwRaR6mTAx3yHEBGZi0ayq0wqk1tNOIqtE8/MYnTvVnZ/90MTy2MHdtL+uJdS6NvLwKbbsFicePtKlj7j9QT1zYfdv+e336HvrhvAQfM5T6X1gucAsP9nVzP44O0kl69n6TPfBEDfPT+hONRL6/nPWZwnJ1I7XpDPpr/lO4SIyGw0kl193osK9hEllpzA6pd/itUv/xSrNv4/LFFH46kXUZ86l9Wv/AyrX/FpEp1rOHjrdYfdd2R3nr67bmDlyz7Gqld8isEHbmN03zaKw/0Mb/sjq1/xaZwrMrI7T3F0mP57fkTLn6U9PEuRqvehVCaX9B1CRGQ2KtlVJJXJPQJ4ue8clWTo4btItK8i3rachvXnYUEMgLrVpzHWu+ew7Uf3bqVu9ekEiXosiFG39mwG7r8FMFxhDOccbmwEC2L03HY9LY96NhaLL/KzEqkJJwH/4DuEiMhsVLKry4fRz/So9P/xJhrPuOSw9X13/5CGk84/bH1y6YkMbbmHwmAPxdEhBh/8LYWePQR1jTSe9hh2XHMl8bYVWF0TIzvuo/GUCxfjaYjUqrenMrklvkOIiMxEQ2xVIpXJ/QXwNN85KokrjDK46TY6Hr9xyvqDv7oWghhNZ1562H0SS9fS+ujL6L72HViinuTy9RCNfrc9+jLaHn0ZAHv/75O0X/xSeu+6gaGHfkdieYr2x7xowZ+TSI1pB96NRrRFpAxp1LMKpDI5Az7iO0elGXzwdpIrNhBr6phY19f1YwYeuI2lz3ozZjbj/VrOeQqrrvgEK1/yIYL6FhIdq6fcPrLrAQDiHWvov+cnLHtuhtHdDzO6b9vCPRmR2vU30Q7fIiJlRSW7OjwXOMd3iErT/4ef0zRpqsjgg7fT8+tvsfwv30mQqJ/1foX+AwCM9XQzcN8tNJ75+Cm3H/jF12h73EugOAauGK60ADc2XPonISJJ4I2+Q4iITKeSXR3e4jtApSmODjGUv5PG0x4zsW7fDz9HcWSQXde+ne1X/wN7b/g0AGO9e9l13bsmttv97Q+w/T/+ju5vvZfOv/hbYpMO8zdw3y0kV55CvGUJQX0zdatPZ/sXXwcGyeUnLd4TFKktr0llch1zbyYisnh0nOwKl8rkLgZu8p1DRMSzd+az6X/1HUJEZJxGsivfW30HEBEpA1emMjmdI0BEyoZKdgVLZXJnAc/wnUNEpAwsBV7pO4SIyDiV7Mr2T8DMh8AQEak9b05lcjo0rYiUBZXsCpXK5E4AXuw7h4hIGVmHXhdFpEyoZFeuNwAJ3yFERMrMW6JzB4iIeKWSXYFSmVw78GrfOUREytBZwLN8hxARUcmuTH8HtPgOISJSpt7mO4CIiI6TXWFSmVwdkAdWeo4iIlLOLsxn07/2HUJEapdGsivP5ahgi4jMRYfzExGvVLIrz8t9BxARqQCX6+Q0IuKTSnYFSWVy64HH+84hIlIBWoHLfIcQkdqlkl1ZNqKTz4iIzJemjIiIN9rxsUJEx319AFjvO4uISAU5JZ9Nb/IdQkRqj0ayK8fjUcEWETla2o9FRLxQya4cV/gOICJSgTamMjn9rRORRacXngqQyuSa0Q48IiLHYg3wVN8hRKT2qGRXhhcATb5DiIhUqFf4DiAitUcluzJc4TuAiEgFe3Yqk1viO4SI1BaV7DKXyuQ2ABf7ziEiUsGSwEt9hxCR2qKSXf5eho6NLSJyvLRfi4gsKpXs8ne57wAiIlXgIk0ZEZHFpJJdxlKZ3KnAab5ziIhUgRiQ9h1CRGqHSnZ5e47vACIiVeRZvgOISO1QyS5vKtkiIqXz1FQml/QdQkRqg0p2mUplcsuAi3znEBGpIi3Apb5DiEhtUMkuX89EPx8RkVLTlBERWRQqceVLO+iIiJSeSraILApzzvnOINOkMrk4sAdo851FRKQKnZPPpu/2HUJEqptGssvThahgi4gslGf7DiAi1U8luzw9zXcAEZEqpikjIrLgVLLLk0q2iMjCuSCVya3wHUJEqptKdpmJDt13nu8cIiJVzICLfYcQkeqmkl1+LiH8AyAiIgtH5yEQkQWlkl1+Hu07gIhIDXiM7wAiUt1UssuPSraIyMI7L5XJ1fkOISLVSyW7jKQyuRhwvu8cIiI1IAk8yncIEaleKtnl5RFAo+8QIiI1QvOyRWTBqGSXF00VERFZPCrZIrJgVLLLy4W+A4iI1BDt/CgiC0Ylu7xoJFtEZPGsSmVyJ/oOISLVSSW7TKQyuTbgdN85RERqjKaMiMiCUMkuHxegk9CIiCw2TRkRkQWhkl0+NB9bRGTxaSRbRBaESnb50HxsEZHFd1Yqk9OniCJScirZ5eNc3wFERGpQA7DWdwgRqT4q2WUglck1AGt85xARqVGn+g4gItVHJbs8bEA7PYqI+KKSLSIlp5JdHjb4DiAiUsNUskWk5FSyy8PJvgOIiNQwlWwRKTmV7PKgki0i4o9KtoiUnEp2eVDJFhHxJ5XK5BK+Q4hIdVHJLg8q2SIi/sTQvjEiUmIq2Z6lMrkkOkariIhvmjIiIiWlku3fesJRFBER8UclW0RKSiXbP00VERHxTyVbREpKJds/zQMUEfFvne8AIlJdVLL9W+87gIiIsMR3ABGpLirZ/umFXUTEP70Wi0hJqWT71+E7gIiIqGSLSGmpZPunki0i4l+rTkgjIqWkku2fSraISHno9B1ARKpHxZRsM+sws0f6zrEAVLJFRMqDpoyISMmUdck2s5+ZWauZdQJ3AVeb2cd85yoxlWwRkfKgki0iJVPWJRtoc871AM8HrnbOPQp4sudMJZPK5OqBet85REQEUMkWkRIq95IdN7NVwAuB7/sOswA0ii0iUj5UskWkZMq9ZL8XuAHY5Jz7jZmdBNzvOVMpqWSLiJQP7fgoIiUT9x3gSJxz1wHXTVp+EPhLf4lKTiVbRKR8aCRbREqmrEu2mV0NuOnrnXOv8BBnIahki4iUD41ki0jJlHXJZuo87HrgecB2T1kWQrvvACIiMiHpO4CIVI+yLtnOuf+evGxm3wR+5CnOQtALuohI+Yj5DiAi1aPcd3yc7hRgne8QJVT0HUBEKt/gg7ez7QuvYdvnX83BW6877PahLfew45p/5OEPP5v+e2+eWD+6dys7rvlHtn/p7xne9kcAXLHArv/8F4qjQ4uWv4yUrGSbWcrM7pm27t1m9uajfJyfmdn5x5Hj3Wa2zczuNLM/mNlfHcdjvdfMSnYYXTN7tplljmL7djN77aTlS83suI48ZmYvM7N7zOz30ffnqH4+R3jca8zssuj6f5jZmTNsc4WZfXqOx7nUzB5TikxzMbNC9Htyj5ldZ2aNc2z/z9OWfxX9mzKzF09af76ZfbJEGfNm1mVmd5nZjWa2shSPu1DKumSbWa+Z9YxfgO8Bb/Wdq4QOm28uInI0XLHAvh9+luUveA+rX3UV/X/4OSN7Nk/ZJt66jCXPeD1NZz5+yvreO/+P9sdfwbLnvo2e2/4nXPe7/6XprCcSJGryEP5l/enufJjZTG8UPu6cOxd4DvB5M0scy2M7597pnCvZp8nOue8657JHcZd24LVzbjVPZvZ04PXAU5xzZwHnAQdL9fjjnHOvcs794RjvfilQ8pJtZjP9rg865851zp0NjAB/O8fDTCnZzrnxnCngxZPW/9Y5d+VxxJ3uCc65c4DfTs8As/4f8KKsS7ZzrsU51zrpcur0KSQVTiPZInJcRnbcR7x9FYn2lVgsQdMZlzB4/61Ttom3rSC5fD3Y1Jd8i8VxYyO4sWEIYhSH+hjcdBtNZz9xMZ9COVm0P87RCPWHzOw2M7vPzC6O1jeY2X+a2d1mdi3QMOk+TzGzW8zsjmiksTlanzezd5rZzcALZvuazrn7gQGine7NbIOZ/cDMbjezX5jZ6WbWFj1eEG3TaGZbzCwxbXT2UWb28+i+N5jZKjNbbma3R7efY2bOzNZFyw9MHxmdPJIbPfYnzexXZvbg+NeZJgtsiEZbPxKtazazb5nZvWb2dTOz2fLN8HhvA97snNsefX+GnHNfiO7/ajP7TTRi+t/j2WfLaaFPR6PhOWD5tJ/1+dH1l0c/758Dj520zbPM7Ndm9jsz+5GZrTCzFGHRfUP0nC82s2VRnt9El4nHmPRY9WZ2dTTi+zsze8Kk7/d1ZvY94MYZvh+T/QI4ObrfS6Pf0zvN7PNmFjOzLNAQrft6tF3fpJ/TxdFtb7BJnziYWaeZfTv6/b7VzB4ZrX+3mX0p+l49aGbzKeU3TcrYZ+EnLb8GLjKzJ0XPvSt63Lpouwuin91d0XNqiZ7PR6Lv591m9ppo21VmdpMdGt2/ONr2mmi5y8zecKSAZf2u3cx+7Jx70lzrKphKtlQFS+zbaowWrGgEgDkIimDj1x2Yc5izaNmBCwgohts5w4Cg6MJ/nWE48/usKkP3zj81FpJ19Wt7d+0DSATW1LvrweTa3l37p2/bP9S7pLN/z+Dy3l0DAMtOuSD2px99YakrjNmpj3/Z3p0/+lzzhkc8YaCjr3t4sZ9HOXBMlITFEnfO/bmZPQN4F+EZjf8OGHDOPTIqIHcAmNlS4O3Ak51z/Wb2VuCNhOeTABhyzj3uSF/MzM4D7nfOdUer/h34W+fc/Wb2aOAq59wTzewu4PHAT4FnATc450aj/oqFI+GfAp7jnNttZpcD73fOvSIqeK3AxYQjjRdbWP67nXMDc3w/VgGPA04Hvgt8a9rtGeDsaFQeM7sU+DPgLMKDIvwSeGxUtA7LB0w/MtnZwO2zZLl+UuF+H/DK6DFny/k84DTgEcAK4A/AlyY/YFT03wM8inDE/KfA76KbbwYudM45M3sV8Bbn3JvM7HNAn3Puo9FjfIPwk4mbLXwDcwNwxrTsrwNwzj3CzE4HbjSzU6PbLgIe6ZzbN8vzHh/lfjrwAzM7A7gceGz0O3AV8BLnXMbM/n78ZzFNhvDNyzOjx7t00m3vAX7nnHuumT0R+Aow/hinA08AWoA/mdlnnXOjs+UEngl0RdebgHucc+80s3rCc6o8yTl3n5l9Bfi7KPu1wOXRuVdagUHCn+1B59wFURn/pZndSHi28Rucc++3cHS8Mcq6Jhrtx8yOeACLsizZ0TeoEVhqZh2Ef6sBWoHV3oKVnkq2VIX2DR/bPmZjfz7Xdi66zPcX35xz5igGjiLhv84cRYOiOVxQPLQc3eYCR9HCC9OvBw4XPQ5BMboPOCvigvC+0ddz2MR1CIoTt7kgfKxDX3t8ObydoBguT2w79b4W3Zdw2TFxPwdRjvF1k7clmHizMnEbQRFu2r1/6e97+5e8avNHHwoc3Lhr3/IHDw62vvbBcNnCiwXO8aEDW+vO37p571MLN+81F70JOi3xUOAS5B/8asPV23efeGXDloevuvnf1xecs1etW7bl5Mb6IXNY9IbJJi6Tl8EIH2/iNhzB+DLj24bbBdELejDptgAmbjdjYjmIJtUdWo62gfDxx2+b9m8wednm/6ltO7x8npvOabbpgJPXXx/9ezvhR+wAlwCfBHDO3W1md0frLwTOJCwAEO44f8ukx7r2CFneYGavBk4CngZg4Sj4Y4DrxsszUDfpsS4nLIEvAq6a9ninERbUH0b3jQE7ott+RThCewnwgejrGeHI6Fy+7ZwrAn8wsxXz2B7gNufc1ug53Un4fTxwhHzzdXZUrtuBZsIHFfT3AAAgAElEQVQye6SclwDfdM4VgO1m9pMZHvPRwM+cc7ujvNcC4+X3BODaqIgngYdmyfVk4MxJP7NWM2txzvVO2uZxRG8InHP3mtnDk77OD49QsBui7yGEP68vAn9D+KbgN9HXbAC6Z777vDyO6HwnzrmfmNkSM2uLbss554aBYTPrJnyzsnWGx/ipmRWAuwnfeAIUgPGZDqcBDznn7ouWv0z4xuPHwA7n3G+ir98D4SdEwCPt0KcnbYT7AP4G+FL0pvLbzrk7zexB4CQz+xSQY45PBMqyZAOvIZwntZrwxWf8t6kH+IyvUAtAc7Kl4sXo7W61pj37Sj+VEWdmzogVF/2oD5UziD6wqZHub/fx0ctiKwF2fx+ggQ88M7Zm+rZbB+PsPDex4hcXHP7Sv/mqfax462redPOBpY2XLCG5NMlbvtu9dO3frl3w57BYZnvTNrEM/bfO/TDztZfDz4XQydTyNP6JQYGpf49n+ttghAVpth0X+4+Q5ePOuY+a2fOBr5jZBsI3HgdmGYn8LvBBM+skLFjTC6MBv3fOXTTDfX9BOIp9IvAdwv2oHFMPyTubyZ+gzPc/4eT7jH8fj5Rvst8z8/MDuAZ4rnPuLjO7gnBu9Fw55/M3fbZtPgV8zDn33Wjk992zbBcAFznnBo/wNY70vTvS78ng9N8HC5v1l51zbzvC/Y7GTNnGvycz/Sxn8gTn3J5p64aiNzizfY3x9bP93/oH59wNh91gdgmQBr5qZh9xzn3FzM4BnkpY3F/I4Z+QTCjLOdnOuU8459YTftxwknNufXQ5xzl3xD1xK4xGsqXitTGwozPefqQXfFlADesbGN41zMjuEYpjRQ7++iAtf9ZyVI/Rf28/iY4EdSvrKI4UJ8aBiyPV9RLlzKwYWGwsZomxuNWNJKx+OGmNg3XWPFhvrf31VrK/ic65PmCHmT0JwrmohKO6Nx/xjuE805dE9zkbeGS0/lbCqRDjc1AbJ00BmG+m6wmncGyMRvEeMrMXRI9nUXkYz34b8Ang+5PKy7g/AcvM7KLovgkzO2tS/pcSTkspAvuAZxBO5ThevYRTCeZypHyTfRD4sEVHqDCzOjs0F7iF8OeXIPp5zOEm4EXRnN1VhNMepvs1cGk0eptg6vz5NmBbdH3jpPXTn/ONwN+PL5jZTG+SJv8OnUp4VLY/zeM5zOTHwGVmtjx6vE4zOzG6bdRm3on2SD+nydkuBfaMjyiX0L1Aavz/CvDXwM+j9avN7ILo67dEU2NuIJxOkojWn2pmTdHz7I6mDX0ROC+athVE+we+g3Bn2VmV60g2AM65T0UvMmcSnoxmfP1X/KUqqer6CyY1qd317uuIr4gz9rDvKDXJYsbql64m/9E8rujouLiD+jX17Lp+Fw3rG2j9s1YGHhxg86c2U+gv0HtnL93/080pHzgFAOcc3d/rZt1rw6OjdlzawdbPb8UVHKs3VtPsvHk50vzPY/Ey4DNm9m/R8nuccw/McZ/PAldH00TuJCy7RHOLrwC+Gc0bhfCj8vtmfJTZvRf4hpl9gbDsfNbM3g4kgP8E7oq2uxa4jqkjuERZRqKP1j8ZfdQfB/4f4ehxPppWcFO0+c3ACc65w/YROFrOub1m9ksLD434f4Qf18+03az5pm33v9F0jx9FI7aOQ/Oo30FYih8mnPc7V7n/H+CJ0bb3EZa66bl2mNm7Caf57CCcbz/+Kd27CafubCN8Q7U+Wv894Ftm9hzgH4ArCX+n7o6e100cfhSQq4DPmVkXMAZc4ZwbnjTFZN6cc3+Ifj9utHBn2FHCEdyHCef0321mdzjnJr8RuRsYs3Bu/zUcmnc+/jzHf78HmPqGoiScc0Nm9nLC72eccNrH56Lfi8uBT5lZA+F87CcD/0E4zeiO6PdgN/Bcwt/9fzKzUaCP8P/zmij/+BvyI47wm3PlO2PBzN5F+CTPBP6XcCL+zc65mfY6rjipTO75HJpDJFKRzi/e+b0THtvd8OP9P3kUdtjH4yKV5EtdG7te6TuEiFSHspwuMsllwJOAnc65lwPncGjHjGpQvu9wROapheGhZU17Wyk2zjVCJ1LuSj2SLSI1rNxL9mA0p2ssOtRKN+He0dVC00Wk4jUxVFxl2zuLxVXTd0QRqTQjvgOISPUo95L92+gYhF8gPMrIHUTz06pETZ63WKpLE8Mst+4VRXeqRgGl0qlki0jJlPuOj+OnT/2cmf0AaHXO3X2k+1SYvb4DiByPBgb74gHxBgZbjFOCcF8kkYpVkyfhEZGFUdYj2Wb24/Hrzrl8dGD+Hx/pPhVmt+8AIsejIbFrW10hVgQoxlfHcRz3EQREPNJrsoiUTFmWbAtPy9pJdMbH6LiMnWaWorrO+Kg5rFLRrH7n3gYXfiCWbC/GrdDwoOdIIsfjaM8KKCIyq7Is2YRnfLyd8Dz2t0eX3xKePapqTkaTz6b7CY/TKFKRgvrtA3UWCwDqW0YTNrRSI9lSyVSyRaRkyrJkTzrj4/uBc6PrVwMPEh7EvZro40mpWLFEdyHhwjMcLA32Ft3oSZrTKpVsp+8AIlI9yrJkT3KZc67HzB4H/AXhmYM+6zdSyWnKiFSs9tEhR6IQAKxke+DGTo/NdR+RMqaRbBEpmXIv2YXo3zThKTG/AyQ95lkIGsmWitVYqAuCYKQIsJpt9b3BSUmcjpojFam/a2NXr+8QIlI9yr1kbzOzzwMvBP7XzOoo/8xHSyPZUpHijI02FBoTxMJDC69ie8toe0NdTDs/SmXSVBERKalyL6wvBG4AnuacOwB0Av/kN1LJaSRbKtLKYNe2xGh7YLFRB7CCnUuK7cklNrSyx3c2kWOgqSIiUlJlXbKdcwPOueudc/dHyzucczf6zlViKtlSkRLJXd3NI0vGgmDEANo4sNS1JFbTt0HzsqUSaSRbREqqrEt2jdB0EalIVr+jt3mk0xEbAyDABUHM9dnghoTnaCLHQiPZIlJSKtn+dfsOIHIsivW7RhtH2gMLxiZeR+oZ2nuQVAynN49ScTSSLSIlpZLt3ybfAUSORbFub6x+rCVp0Ug2QAs9fSNN9cPxQv1DHqOJHAuNZItISalk+3c/UPQdQuRojSV6GxNjFregaOPrOtk74tqSQWxw5UGf2USOgUq2iJSUSrZn+Wx6GMj7ziFydJwbiQ13xgouFsQPrV3BTlfsSLZb/0nx2e8rUpY0XURESkoluzzc6zuAyNFYxb5do86tCIqFhMXcpPXbksXW5Fr6Tl7uMZ7IsdjqO4CIVBeV7PLwJ98BRI7Gsvj2nU1DQcExFrfgUMlezbYmkkH7/tG1reZ0eEqpGDu7NnZpZ10RKSmV7PKgkWypKFa380DLcMt+54bGLOYm5mSvZEc7wGiibk+iUKedH6VS3Ok7gIhUH5Xs8qCSLZWlftdg83BHj3ODI8Gkkr2U3SsAXFO8Jz64Qjs/SqW4y3cAEak+KtnlQdNFpKKM1e2mZbhzgOLgKMGhkp1ktB5X3F9sS7hY/3rt/CiVQiVbREpOJbsM5LPpXcB+3zlE5mu4rifZPLJkxLnBgk0q2QBJRncX2+ta6D95ma98IkdJ00VEpORUssuHRrOlYgzHB1uaRzqLzg06C9yU15Em+nuK7YnVe0dOPNmc0xlNpdwNAvf5DiEi1Uclu3yoZEtFaKXv4EDMLW0c6TTckLNYccpIdjv7BqmPrxwiOVKnnR+l/N3TtbGr4DuEiFQflezy8UffAUTm4xTbun0EVtePtSbMjYEVp7yOLKM7LCxx25IYWt7jJaTI/Gk+togsCJXs8nGH7wAi89GU3LEXs/pkoak+buaml+xVbI8BuIb4/njfer3GSLnTfGwRWRD6A1g+bgX0kaWUPVff3Q8Qd/VNccywwpTXkdVsawAotiYK1rdBZ36UcqeRbBFZECrZZSKfTfcCXb5ziMxlrG73GM45c8m2IMARHDaS3QrgOpIN3aPrN5hzu/wkFZmTQyVbRBaISnZ5+aXvACJzGa47ELQOsA8LOh0uCKw45fZl7FoKUGxLrhikrrG+UJf3EFNkPh7q2tjV6zuEiFQnlezyopItZW8g0d/YOpDch1kzQRELClOOLtJKbyfODbmm+FoHQ3WDyw74yioyB83HFpEFo5JdXlSypazVMTLUFx9rax5p7wFwuJiZO2y7GIVdmAUEtjnRf2Js0YOKzI+miojIglHJLiP5bHozsNV3DpHZnGQ7tvUFwbLm4c4BABcjCKaNZAM0MHgAwNXH9ljfhqWLnVNknm7xHUBEqpdKdvnRaLaUrZNt894CrGoeWTLsXGHUmZlZ8bCS3crBfgDXkhjZOXrSyYFzOxc/rcgRDQK/8B1CRKqXSnb5udl3AJHZxOp39WAWNI90juGG+1wQG5upZC9hzyhAsSNZ109Dc8NYMr/oYUWO7OddG7uGfIcQkeqlkl1+NJItZWukfs8wQONoJ84N9xMEhcDcYa8jK9lhAMX25FKA+qGlBxc3qcicbvAdQESqm0p2+bkb6PMdQmQmg8n9ANSPtcdxQwMuiBVmGslezbYEgGtJnOigkOw/Ua81Um5+4DuAiFQ3/eErM/lsukB49keRsjNQ15sASBaa65wbHHJBUDBzh5XsVWxrBiCwJMaWWN9JSxY5qsiRbO7a2HWv7xAiUt1UssvTT3wHEJnOKBZ7EyNNAHHX2EBxcJggcGbFw15HVrCzc2KhLta9dXTDhsC5HYsYV+RINFVERBacSnZ5yvkOIDLdCbZnx75Y0IFzzqhrdW5w1AUxh7nDjoPdyd7lOFcEKDbHB3ppamsaS+QXPbTIzFSyRWTBqWSXoXw2fTew2XcOkcnOsIe7B8xWNA9ywFm83RUHCy4I3EzTReIUEobbA1BsT8YBGoaW9Cx2ZpEZFIAf+w4hItVPJbt8aTRbysrq2NYeZ7aks5d9YJ24waILAsyKM57RMcnwXgDXXtcBUN+/7rAyLuLBr7s2dh3wHUJEqp9Kdvn6vu8AIpMV6ncPADQPt/ZgFnduiPBkNIcfwg+gmb5egGJrYi1ArO+kzpm2E1lkmioiIotCJbt8/QQY8B1CZNxg3f4CQPNIR3iISTeECwKDw+dkA3SyLzzRRyJodbBz8+jJG2LObV+svCKz0KH7RGRRqGSXqXw2PQTc6DuHyLi+ut4YQPNI5zCAc8OxaCQ7PtP2y9npJhaSwbYDtHQ0j8UfXpSwIjPbC/zWdwgRqQ0q2eXtet8BRMb1JIbqAJpGlowC4EbiLrCAWV5HVrF9YoTbNcX7AJqGOnXmR/HpR10bu4q+Q4hIbVDJLm/fA0Z9hxDp5ODe3QmaITylOoBzo0lnFtgMh/ADWM22xvHrxbakAdT3r9XOj+KTdigXkUWjkl3G8tn0AXSoKSkDpwdbdu4PgvAoIWPt4euGG6tzZnGYebrISna0jV8vdiRbAeLa+VH86Qf+x3cIEakdKtnlT1NGxLuzLH9g2Gw1QLLYlgjXFuqdWcyMGUv2MrqXjV8vtiXXADw0espJMee2LXxikcN8q2tjV5/vECJSO1Syy9+3CU+eIOJNc2LnAGbRKdWbGsK1xaaiMetIdiMDrTjXC0BdbJmD/ftoW9IyFtOJlsSHa3wHEJHaopJd5vLZ9G7gR75zSG0bqts3sW9A4OqbwmuuyUEcmHFONkCc0d2HFmwrQMtQh04EIostD/zcdwgRqS0q2ZXhP3wHkNrWW9/rAJoGXU8xSLY55xzQ7IwEzDxdBKCRwYlC7RrjBwAa+tfqdUcW21e6Nna5uTcTESkd/bGrDN8Bun2HkNp1MDmYAOjsZa+zWAeM9AHOQYLwMqM2DkycUKnYliwCJPrWty90XpFJHPBl3yFEpPaoZFeAfDY9CnzFdw6pTQ0MD3QnCg0AHf2xg5i1Uhzuj1liyOGSHKFkL6V7Yn+CYnuyCeDB0VNPijm3dcGDi4R+0bWx60HfIUSk9qhkVw5NGREvTrZt23bG400ATcOdfQDODfXHLTHocPHZji4CsIrtE8fFdu3JlQDddCxrHwu086Mslmt8BxCR2qSSXSHy2fSfgF/4ziG15+zgof09QbAUoGm0cwjAuaGhWJAYxgqzFmyA1WyvG7/uGmJrXHisYlqGOnTmR1kM/cB1vkOISG064h9IKTv/AVzsO4TUltMt3/99WA3QNDp+SvXB4YQlhrFi3ZHuu5LtLRMLZkbMNlNwZzT3r4GWfQsZ26uRvSNs+8I2xg6OgUHHpR0sfcpSDt52kO5vdzO8Y5gN79xAw/qGGe+/54Y97P/5fjCoP6GeNa9cQ5AM2PK5LQxtHaLl3BZWXrYSgO7vdFO/tp7W81oX8ylWiut1bGwR8UUj2ZXlOkCHP5NFVazbN4JZHKBhtNMBuOLgSDxIjpgVj/hGfQU7l0xedvWxfQCJvvUdC5W3HFjMWPmilZzywVM46R0nse/H+xjaNkTdCXWs+4d1NJ7aOOt9R/ePsveHe9nw7g2c8v5TcEXHwV8fZGjLEACnvO8UBu4boDBQYPTAKIMPDqpgz+4a3wFEpHapZFeQfDY9CHzDdw6pLb11PcXx6/VjneEcazc4FrfEiAXFWY+RDdDO/mU4Nza+7FoTYwCbRk89Me7clgWK7F2iPUFDKhyljjXEqFtdx9j+MepX11O36oiD/wC4oqM4UsQVHG7EEe+IQwzcqMMVHW7MQQDd13ez/PnLF/rpVKqHgZ/6DiEitUslu/J8wXcAqS37kwMTrxPJYlscwLnBIkFs1IJicfZ7QoALAooTh58sdtTVAWxn6aqOMauJnR9Hdo8w9PAQDRtmnhoyXaIjwdKnLeW+N93Hva+/l6AhoOXsFupX15PoTPDAux6g7YI2RnaNANBw4vwetwbp2Ngi4pVKdoXJZ9N3Arf7ziG1IUZhbE9yLDm+HHctdQCuOOhcEBsNrDg2+71DdQxNTL4uticnhl3bhtp7Sp233BSGCmz+9GZWvnglsYYjDvofuk9/gd7f9XLqR07l9I+fTnG4yIFfhbPEVr1kFSf/68ksffrScBT7ecvp/m43mz+zmX0/q9457segAHzJdwgRqW0q2ZXp330HkNqwznZt35yIT8xvCGgMJxO7IbBgzIK5S3YLvRM7nrmm+DoHowAt/att9ntVPjfm2PLpLbRf1E7b+W3zvl/f7/tILE0Qb41jcaP1/FYGNg1M2abnjh4a1jdQHC4yvG2Yda9bx4FfHaA4fMQPFmrJtV0bu/K+Q4hIbVPJrkxfRWeAlEVwpm3e3R2LtQA0DLveYlDXAuDccEAQFOYzkr2EPSMTC4HFCdgMkOhLzb95VhjnHNu+tI26VXUsfdrSo7pvYkmCwQcGKQ4Xcc7R/4f+KfO43Zhj7w/3svTpSymOFGH8rYoLbxMc8EHfIUREVLIrULQD5Md955Dq94jgwf6BwFYALOllXzFIhMXYDceKQWzUgsKcQ6cr2DllG1cX2w1w/9hp6+LOVeW87IH7BzjwqwP0/bGPTe/YxKZ3bKL3rl56bu/h3jfcy+ADg+Q/nif/0TwQHlEk/7HweuOGRlovaGXTuzax6e2bcEVHx6WHDsay98d7aX9sO0FdQP3aenBw/9vvp/GURmJN85uSUuW+17Wx6x7fIUREdJzsynUV8Fag3XcQqV4nBNvGimbLAZYcdAfBzgBwbjRBEPTPZyR7FdunnHbdtSSGGCyw2a1Yc/IYv9yVYN3CpPen6dQmzr7m7Blva33U4YfbS3QkSL0xNbG84nkrWPG8FTPef+lTD42Mmxlr/27t8YWtPh/wHUBEBDSSXbHy2XQP8CnfOaS69df1jI5fbx1s7sMsKsyjdS4IihYUC3M9xiq2NU1eLrYnJ0p3+1CbzvwopfTjro1dv/YdQkQEVLIr3SeITlMtshD21R3a4a5xdMngxIIr1LsgKAZWmLNkr2THlBPPFNuTEyeoae1bVdU7P8qi0yi2iJQNlewKls+m9wKf951Dqtee5MjElLKm0c7RQ7cUG4tB4II5jpMNsJTdU86W4loT61y4cxrJ/hN1qkIplVu7Nnb9xHcIEZFxKtmV76PAsO8QUn1WsK97czI2UbIbxpZMGrUuNjoL3Hx2fKxjpAHn9k+siAWNGNsA/jR2RtXu/CiLTkcUEZGyopJd4fLZ9A7gGt85pPqcFmzZuSV+6BjZ9YUpsz6aiwEusLnnZAMkGdkzdUWwE+Aht2rtslGq9vTqsmi6gO/5DiEiMplKdnX4EDDnUR5EjsYj7KGevbGgc3w5UeiIATg3NgwkXRDYfEayARrpn7KDY7EpMbEvQedQS9Wf+VEWXFanUBeRcqOSXQXy2fRDwDd955DqclaQHxsyWzm+HKc1PL26G+4FcBYwnznZAB3sH5q87NqTE689bf2rVI7keDwAXOs7hIjIdCrZ1eODRDuTiZRCY3xPEbOJHRMDmuoBnBsaACgGWGCFef3OLWPXlE9aih3JibM9JvvWteKcfnflWH24a2PXvKYtiYgsJpXsKpHPpv9IeLp1kZI4kByYGKWuG3H9xaChGYCoZDszC4LivMrxKrZPORVhsTUxcQaVewtnnJCAh0sSWmrNNrRPioiUKZXs6vI2dNxsKZHddSMTxbizl72FWF0LgCsODgMUjWC+c7JXsb1+yopkrMPBHoD73ZoTl4+6rSULLrXkX7o2do34DiEiMhOV7CqSz6a3E+4EKXJcmhjs3ZJk4kQxyw66HmexTgDnwpLtzILA5j2SPcO5xIOoWJst1c6PcvRuAb7iO4SIyGxUsqvPRwEdd1iOy6m2dfsDicTE68OS3lgfZuE8ajc4BlA04kEwvznZy9m1bPo61xifKNbt/SuOO7PUlCJwpY4oIiLlTCW7yuSz6UHgrb5zSGU7O3ho//Z4vHF8uXGkY+L86q4YlmxnxCyY12wRWunpxLkpRxgptiUmClJd/7oW7fwoR+FLXRu7fus7hIjIkahkV6F8Nv2fwK9855DKdbY9NHQwFiwZX24cW3LolOpu0AEUjdh8jy4CEKPQPXm52BHO8Qb4w9gZaxKQP57MUjMOAP/sO4SIyFxUsqvX69Eh/eQYnRxstVFYPb5cP7Zk4hB84wPSRVzc5nl0EYB6BvdPXi62JVaNX7/XrU2tHC1uO67QUive1bWxa7fvECIic1HJrlL5bPo3wNd855DKNJLoBbNJp1TvnCjTzg0TEBsumosFwfwPT9zKwalHvmmIr3LQA+AIgmVDzdr5UeZyD3CV7xAiIvOhkl3dMuiQfnIM9iRHbPJyothx6LXCDcdjQXyoSDER2PzmZAMsZc/oYStjtmX8akf/cn3yInO5smtj19jcm4mI+KeSXcV0SD85FgnGRrYlp5bnuGtLTCy4kXjckkPRdBGbfv/ZrGDHYetcY+zA+PX6/rU686Mcybe6Nnb91HcIEZH5UsmufjqknxyV9bZj26ZkYkrLDmiZmDri3GgyHiSGirikzfM42QCr2Z6cvq7YmpwY3f792JkrtfOjzGIAeJPvECIiR0Mlu8pFh/R7je8cUjnOtIf3PphMTIxcJ0fdYDHWOHE4P9xYXdwSIw6XCILCvEeyV7Gtefo6156ceNzfu9RJq0cL2vlRZvKhro1dGiwQkYqikl0D8tn0D4Av+s4hleERwUP9u2KxiULc0cveQqx+UkEuNiSCulGHqzM7mukiOzunryu2JyfOQlMkiK0YatLOjzJdHviw7xAiIkdLJbt2vBHYMudWUvPOsIfpC4Ll48vLetzBQizRdmiLYlM8SI44SAZHUbKXsGf59DnXril+goPh8eVO7fwoh3tD18auobk3ExEpLyrZNSKfTfcAr/SdQ8rf8mB3oggTI8zLDtIPwaRRaNcct+QYUHc0Oz7GKSQMt2fKSrMYgU1MA2jsX9OCc/M/ZIlUu691bez6tu8QIiLHQiW7huSz6R8C/+47h5S3fclhw2zitaFtoGlg/JjZzhULQGMiSBaABjM375INkGR4z/R1rj6YWHfX2Nkrko6Hjj29VJGHgNf5DiEicqxUsmvPm9ARHGRWzu1MTj3BTMNYx/Chm4f7AMzioxixo5mTDdBMX+9hX7ElMfH4XW79htVjhcOP9Se1pgC8tGtjl+boi0jFUsmuMflsuo9w2ojmvsphVrN356ZkYspJY6aeUj0s2USnejQrHtVrSCf7DptbW2yvmzg8YIFYfPVgg4qVvK9rY9evfIcQETkeKtk1KJ9N/wT4rO8cUn7OCB7uvn/aMbLriksOvSFzQwMAzoIxADN3VK8hy9l12HzrYkdy6eTlJQPL9Aawtt0C/KvvECIixyvuO4B48xbgacBJvoNI+XhE8FDPDxPxusnrEsVOGx/adm5oCKAYC8+nPt+R7NtuG+Cqz+xld0/mMSNWR9DUQsMznkvTi1+Ba06sc1AY+OMvYvt/9iU+3L/36YVCkbWvW0v9CfVs+dwWisNhNw8SAateuoru73Sz7u/Xsfn/t3fmYXZUZf7/vHWX3tPd6ewhpAHZExZBxRGRRVFs1xHEBWncxRnBdcwM81MQ0XYflcWVAXVURBgUoiijLIosWYBcAmEJuUnozp70frvvUuf3xznVt/r27c7St3PTnffzPPXcW6dOVb11qurU97z1nlPXbmDh5QvxKtRfMEXoAd6baE3kdptTURTlAEf0K8YHL81Llp4B3AfsVVytMrlIvbCCnX/5Efg+FQuOZ7D9afzBfvBz4EUw2UFMZgB8n3hEsmnfRDHYq8LYn6qKOvoHh4dTRyIRRHL+jBlRb+fOLL4P2SzU1nr09fk0zYjQ3eVz3fVz+filHcycGWFXTyQ7WDcrKnV15J57Bm/2XKZ94vNUvRBt33bbl+cbY6iSzCCz+3ON5zRV9z3dR/8L/VTMryCzNcP8D85n43UbaXp9E37Kx6vyaDy9sRzFqkwMrYnWxM/KbYSiKEopUPfPQUyyreUB4BvltkOZOIyfY+c9NzDrgquY+4Hv07f6PhrO/jBeNI5U1rGRUJYAACAASURBVCKeh1TW4VVOg2gFmZyJ4sH8j8wHIFof5ZDGBZ2V8VreftpHqamoGRAg6kVynucRjYqJxYTqao94hXDssXH6+nw+/JHppAd9KiqE39zSxbx5MT75qZmYXM7z6uoRiRA76VQq/+k1dH/3q/Q8dqcn8UpmX3Al0173L56/3frO+1/oZ9op0xCx7cBsV5ZMV4a6k+vofrybhlc1lKtoldLzaxXYiqJMJVRkK/8B3FtuI5SJIb3pWaINc4k1zCGzdR3R+pkMPP8I0cZ5VC08EYnGqT7spXgVVXhVdYhA3Ul1ZHdmiTXFyKVyBrxILBJlWnUTOZP1PE8AxPdtxEhnZ46qKo+BlGHmTBuBdsYZNfT0GFreVMf69Rmqqz22bslSEUekspLchnVUnvV6clu3EGmaSaazI2Yyabb9ro2t9/wkZgaNZLoyxOpj9D7ZS2pDCj/js/H6jcw5fw7bfr+NWW+eNSS+lUnPBuDSchuhKIpSSjRcRKF5ydJZwArgkHLbopSWvjV/Z2DdSprOu4y+NX+nZ+VdSDROtG4GROMMrn+CioUnMrj+CXKDfZi+ndQsrmVg/QAmZ8j15hARGmtn09m7Dd/kDCCeiB+JRj2RrN/UFPH6+3P09hpyOTj22Ap6enK0t2e5uLWRVU+kePHFDNGocMiCGE/EXpFJr3gkJtEoZnCA6d+/id5vfW2nDDLdTw8gnpDd8qxfcWild8hHD6Fijg0RH9wyyNbbttJ0bhMbb9hIxaEVCILEhNn/PHsonzLp8IGzEq2JB8ptiKIoSilRT7ZCsq1lK3ABkC63LUr5EJPzAZPtyjLtlGnUHl8LwPmnfXTrJ1q+zvymw5lV19gb9Tx8Y6SiooLaWsnNnhMlGvU48sg48Tg888wgHR1ZqquFW3/TyerVAzQ0RPj3f59FJm3IrFrhxRadSGTeArwZs+i88nN4c+ZK7aKzmdf6bea+71vEItGcROwwkxuu38DGH25k8683M+ufZ9H+03ZmvmUmsfoYuVSO2W+bzdY7tpaz6JTx8TUV2IqiTEVUZCsAJNtaHgY+VW47lNISrWsi271t6H+udwfR+tk2zc9h/OzQb9T4OQMSnxknNj1Gri+H5+H3ZwZq73z0RlrP/g+293TWGqt9JZVKceG7mnZu2ZxlxkyP7dtzzD8khjFw5JExBgYMtbUeZ5xRwwsvpJlW7/Hss4OYvj4vsqAZf+d26q/4Kv7O7cRfeXrFwIZVAGT7OhnMZGONJ9Vu23rHVma/bTaVh1SS7cyS7cwicSHWFMNkDCZrwAM/rV9in6QsA75YbiMURVEmAhXZyhDJtpbrAe14NIWIzz2K7K4OMp2bic06jGzXNiqPeBmZXR2k1j+ByaSHfrODfYKB7lXd1BxbQ9+aPowPg+lUVSwaZ0vnBjzPy+V8QzQSyYoIfX25aDwu9HQbGho8envsd2t6enxiMejq8lm8uJKqKo/P/9tmMhk7XMngH3+HVFaRfvIxyGSgtqZaYlV9L/7gA7RfdzEgbLlrx/Sex3swYuh6tIuKhRVsvXMrcy+ey5Zbt9D3fB+5nhwbrt3AjPNmjF0QyoHIduCdidZEZrc5FUVRJiEak60Mo3nJ0irsxyBOLLctSmlIrV3Gzr/8GIxPxfzjGOxYQ65nB3geEqvEH+iF3O51jicevhnpMZ42TRgchJoa2LnT4HlgDHiuCe95wsknV9LX5/O5f5vFD3svenTZ937/8mxyLfg+1RdeTN0H/oWKezetNH39L93626uY/c6recuCL/xx6eae8zp+3oFEhAWXLtC466lDBjgn0Zr4W7kNURRFmShUZCsjaF6y9HBgOaADEB9E3BX/j7+1Hpabl/bkiCDtw/cc8YdI7WVvDOYHOn+4AtN3StNh5/whWdn1xled/ot2zzPz92Y/v6D1gT/KW84oTI8/uu1+b1f6NcH8hxq+tfSWudta9vV4lAOaDydaEz8ptxGKoigTiYaLKCNItrW8AFwEaAvsIGKebK9NC/PCaRW56QXXQLYCwPds3SFCZK/3Q0dRd7RfHx82Hl9N39xqjNEv/009vqcCW1GUgwEV2UpRkm0tfwCuLrcdyv6jJ5KOI1IVTouapuGZTLYSwJfgK6Emurf7mUt7XbF0vzE+LTy/MrdoZoUxL+zt9pUDmj8Dny63EYqiKPsDFdnKWFwF3FpuI5SJp4L0QHvcGxFwHaWhwFPtVwP4MlR3xPZ2X7PZ3FQs3a+PDws7We4fffihmdzmvd2+csDyLHBhojWhbycURTkoUJGtjEqyrcUH3od+EXLKc4R0tD9TEesuTBepLxDRpgbAz4eJ7LUnu4FdMzEmO2JBRWSmgc5gNkVFdfNArHdvt68ckGwD3phoTXTuNqeiKMoUQUW2MibJtpZB4G3A4+W2RZk4jveSO9bE48OEbzRn0sarK4yfrgXwMYHI3mtPtofxPPxtRRdGZWN4dnZ/o/YLmPykgDcnWhNry22IoijK/kRFtrJbkm0t3cB5wLpy26JMDItlXWpdLDasPmjoZUcmVjskso1J9wMRQXI5MUHevfZkA1QwuLNYuqmODvN01vTNryrq9VYmCz7wnkRr4pFyG6IoirK/UZGt7BHJtpbNwOuxr32VKcYx3gY2RyO14bQZ3XRmo1U1QwlmsBcgIrFUjlwEjC+yb3VIHd09xdL9abFhceErcsc3VWrnx8nM5YnWxB3lNkJRFKUcqMhW9phkW8tzwBsBjZOdYhwqW2u6PW9mOG1Wp+n1vVh9MG/MYB9A1IulchhPxN/nL/VNZ8dgsXS/saImPP+of8zhzZnsln3dj1JWvp1oTVxbbiMURVHKhYpsZa9ItrUsB96B/WKbMkWoo6c+B3PDabM6GTTiTR9K8FMpgIjEBnPiRzxv30X2bDYXjbU29fE54fk+qmqbB2J9+7ofpWzcAny23EYoiqKUExXZyl6TbGv5M3AJ+rGaKYHg+1tieIgMG65vWn9VOjxutjEDgwAxL5728WMi/j4PxTaXjqIdJk11ZL6B/nDavL6Gkd9yVw5kbgEuSrQmSlI/iEiTiDzups0i0u7+d4rIU6XYR5F9niQib9x9znHv50wRuWui96MoSnlQka3sE8m2ll8Cnym3Hcr4WSDbOtbGY12F6ZW56enwvDH9aYCoxAZ9TGQ84SLzaK8tukBEiMiGcFJV/1zt/Dh5+CXw3kRromTnyxizwxhzkjHmJOAHwHfc/5OwHSvHRET2pXPuSdjQOEVRlH1GRbayzyTbWr4DfKXcdijj41hZv21NPD4iJCPuTx8uYPxUFiDqxbM+flzE32chNYeOhtGWmcrIjvD8ityiRu38OCn4BXDxfv7YTEREfiwiq0Xkz+LevIjIfSLyFRG5H7hcRGaKyG0issxNr3L5Xi4i/xCRx9zv0SISB74EXOg85heKSI2I3OjWfUxE3urWv0REbheRu0XkORH5emCYiJwrIg+JyEoRuVVEal36G0RkjYj8Hfjn/VhWiqLsZ1RkK+Mi2dZyBfCFctuh7DuLvXW9z8ZjI4RR1Ewf9rrfmFQOIOZVpH1MzPP2XWTPYNus0ZaZabFh233YP/aww7Tz44HOzUBrGb7meCRwnTHmeOyHjN4RWtZgjHmNMeZbwHexHvCXuTw/cXnWAGcYY07G1mNfMcak3f9bnAf9FuAK4K9u/bOAb4hI0En3JOBCYDFWmC8QkRnAfwKvNca8FFgOfFpEKoEfA28GXg0M64OgKMrUYp/GuFWUMMm2lqublyztA75VbluUved4Sfp/isXihekRmT6sEW7MgAGISjxnoGI8MdkVpKswphORER5tvyFeEdmUGprvoab+8IFI39OFn8VRDhT+G/hQojVRjtj5dcaY4ENZK4Dm0LJbQv9fCxwnIsH8NBGpA+qBm0XkSGwfk9E+rnQu8BYRCTpzVgKHuv9/McZ0AbgY8YVAA3Ac8KDbZxx4CDjG2fycy/8L4CN7ecyKokwS1JOtlIRkW8u3gY+jnSEnHYfLpopt0Uh9YbpH/fBGuG+Fb8yL5wwmPh5PNkCMdNEx1/2G+MzCtPl99dr58cDkJ8AHyySwAcJDQeYY7jgKh0B5wCuD2G5jzHxjTA9wNXCvMWYR1rtcOcp+BHhHaP1DjTFPj2GDAPeE8h9njPmgy6N1pKIcJKjIVkpGsq3lBuD92AeNMkmYKV2N/SLDwjc832T9SP0w77YxgwIQ8ypyQIVIblznuYa+EZ0tAUxt7FADwwR8Tf+cSozRYSMPLH4IfKRUo4hMMH8G/jWYEZGT3N96oN39vySUvweoC83/CfiEOLe0iJy8m/09DLxKRF7i8leLyFHY8JTDROQIl+/de38oiqJMFlRkKyUl2dZyM/AedBztSYMv6UYjMsx73NjL9ky8bngIiUlHAaJe3Acqx+vJbmBXqugCT2J4DBthZFlucWOVMWvHsz+lpFwPXDpJBDbAZcCpIrLKhXR8zKV/HfiqiDwIhIewvBcbXvK4iFyI9XjHgFUi8qSbHxVjzDasaP+ViKzCiu5jjDED2PCQpa7j4/qSHaGiKAccGpOtlJxkW8tvmpcsTQG3AhpJewDTRNf29fFYNzDMk93UTWcmWl0VTjNYkR2RqI8QFy83rhCBmWzNJTmi6DJTEdkqqdzhwfxD/nHNR2d+/eRTFRXHjGefSkn4fqI1cVk5dmyMuTL0PwksCs1/M/T/zIL1tmM7JxZu7yHgqFDS/3PpO4GXFWT/aJH1bwJuCs2/KfT/r0W2gTHmbmxstqIoUxz1ZCsTQrKt5U5sjGP/7vIq5eMYb+PmZ+OxnYXps7pMby5SMXwsa5OpADCelwPwxB+XyJ5LR2S0ZaY2NhCe76Su8ahURK+l8vNf5RLYiqIokw0V2cqEkWxruQd4Aza+UTkAWSTrutbE4yPCNmZ1MuB70cZhiSZnRbaIARBv30cXAZhH+2idzPAb4iNGeZjfP007P5YPA1yVaE18qtyGKIqiTBZUZCsTSrKt5W/YcWU3ldsWZSSLvHWZ50fqWWZ0SQZkuMgmVw3gi+cDeOPs+DiXjmmjLfMb402FabX9syqwYxgr+5cU8K5Ea+LKchuiKIoymVCRrUw4ybaWFcArgCfKbYsynCOlPboxGh3hUa7KNmQRKagf/GoA44kV2d74wkVmsWXGaMvMtNihpmCos0dyJ0yr0i8/7m/agVcnWhO/KbchiqIokw0V2cp+IdnWshE4Hbir3LYoeebKjmk7I970wvSYP73YyCF1AL6rNbxxdnycRncTxgwWXRjxqpGhodUA+Id/fPMRmUzRsbWVCWEZ8PJEa2JFuQ1RFEWZjKjIVvYbybaWXuCtwH+V2xbFUktqzqDI3ML0mCkIxza5NG6kGH8oJjs37uHbIuS2jrow7m0Oz+6gfsbRKSk+7J9San4NvCbRmugotyGKoiiTFR3CT9mvJNtafOBTzUuWrgG+z+ifMVYmmGoG+roi4iFSW7gswnQZlmAGe4HpAL5bMt7RRQAqSe3so25BsWV+TawvMjjc0b2gvy5HQ3Hnt1ISDPDFRGtizHGgFUVRlN2jnmylLCTbWn4InAOM7slUJpSXSHvH2nisePl7jcOG1zNmYOgT1TlXa5TCkz2N7r7RlpmG2Ij6qTo1Wzs/Thz9wAUqsBVFUUqDimylbLiRR04FNOazDCzy1u1cE4+P+LS555ucidQPf8NgBobCNIY82ePs+AjQxPZRvwzqN1TUF6Y9ml1cV61ffpwIXsR2cLyt3IYoiqJMFVRkK2XFdYh8NfA/5bblYGOxrBtYE4+NiL2o72NHOl4/LJTM+ANDH4fJh4uM35M9Z4yRHf362Igwkgf9Rc1HpLXzY4l5FHhZojWxstyGKIqiTCVUZCtlJ9nWkkq2tVwEXA5owO1+4mhvo7wQj4346mJTN7vSsdr4sESTGjovvhvaz/PG/22YubTHR10YjzQaGCaot9I487gBvUZKyC+xHRw37zanoiiKsleoyFYOGJJtLd8DXg6sLrctBwMLZFvNpki0ujB9VpfpzUarasJpxqSGwjpyYjwoTUz2PNpHdLocRsxrL0xa0F+rX34cP73ABxOtifcmWhMDu82tKIqi7DUqspUDimRbyyrgZcB15bZlqtNA78yuiDezMN1+Uj0+7GuMxk8Nfd3RFyMAnvjjFtmz2dw41nJTHekuTKvpnxUbdXxtZU94EDgx0Zq4sdyGKIqiTGVUZCsHHC585F+Bt1AQLqCUhgi5rEdudhZGjJE9q9NkjXjDP2tuUkPe45Ane9x2NLF9NsaMKtb9+viIZY/6i2prtPPjvpABrsCGh+iXMxVFUSYYFdnKAUuyreVO4ATgz+W2ZaqxULa0b4pFtiMyYpzyulRlDpFhYSTGDEUUmEBke54vhevuLVFyMcHsGG253xgfEU7y99ziBUekM6OuoxTlaeC0RGviK4nWxPhbR4qiKMpuUZGtHNAk21o2A28APo12iiwZx8n67c/H40XfEsSLfFLdmAEBiEh0IId1aksJwkUA4gyO+rbCr4/PK0zrYMbcRSlfr4U9w2A/+nSKjh6iKIqyf9EvPioHPMm2FgN8p3nJ0r8CvwKOLbNJk57F3gt9a+IjP/YCEDWNIzsWmkEPICKxgYzYz9F4Upr+h7X09gxSVXxhVXSugR6BunDywlRNjiJmKsPoAN6faE3omyBFUZQyoJ5sZdKQbGt5AjgF2ymyJF7Ug5XjZIN5Jh4vGjbgyfSRYSAmHQWIerGBnBspW7zcuMNFABrZOfboFhHZWJhU1z8zHo5hUUbwG2CxCmxFUZTyoSJbmVSEOkW+Cnii3PZMVpplU+W6WLTomyzxGkaMnW1MJg4QldhgltJ6smexZcwNmarIzsK0h/zFVdr5sShdwPsSrYkLE62JEeWmKIqi7D9UZCuTkmRby0NYr/ansWP+KnvBDOmevjUSrStMF2N8E2kcIbIx2QqAqBdP5/Dd6CJ+SeqPuXSMGbbmT4uPiBH/u794wZHa+bGQu4ETEq2JX5TbEEVRFEVFtjKJSba15JJtLd8BjgF+W257JhOVpOf3eTKrMH1aHzsG4yM92ZCrAohJPJ0TPwogJfJkz6N9lIBsi2mMj/hgzgYze/4JqWymWP6DkOeANydaE+clWhMbym2MoiiKYlGRrUx6km0t7cm2lguA8wANIdgNc9i5Je3h+TC7cNmMHjrT8WlFPMt+NUDUi2d8TARAZPxD+Fl7NjWMtdxviI9oDAAc1l9zsA9F1w18Djg+0Zq4q9zGKIqiKMNRka1MGZJtLXcDi4CrgXSZzTlgOcbbsCUZjXUgMkIkz+w0PZlYdRHPsqkBiHnxbA7ryfbclx/Hy0y2zBhruamJLjBFhm+sSc2IYUyqFDZMMnzgp8CRidbENxOtCfXoK4qiHICoyFamFMm2loFkW8sXgMXAX8ptz4HIIlnX/Vw8VrRT3KwuUrlI5bAPwBj7RcZagKhUZH1MDEoXk11Dfz3G9I2aQSSCx4gwiIf9xVW1B1/nxweBlydaEx9KtCa2ltsYRVEUZXRUZCtTkmRby7PJtpbXAu8Eni23PQcSi7xkdk1FvL/YslmdJut70YLwjXQvIAAxL+6bQGRLaUQ2QJTsmILRVEa2F6Y94C+ed2Q6fbCMoLEReE+iNXF6ojWxotzGKIqiKLtHRbYypUm2tdwKHAd8EEZ6Qw9GjpD22HOxWNFxxpu6xQeZPizRHxzyMse8Ct9g4gBCacJFAKro7xxruamLjQj/WWfmHXpSKjti5JEpRgr4EnBMojXxq3IboyiKouw5KrKVKY8bheRG4EjgE8DmMptUVubIroYNsWhFsWUVuWk+IsNGFzFmYEhkR60n24psMSWrP+rpLOpZD/AbKmLF0g/vr5rKnR9/gxXXX0y0JsYsH0VRFOXAQ0W2ctCQbGtJJ9targWOAD4PHCyhBsOoYWDOjkik6IgeUTN9xLh8xqSGvqwYlbgxUAmA+EWG+ts3ZrBtzM57fkO8aOfImtSMyBTr/GiA3wGnuQ/K6NsXRVGUSYqKbOWgI9nW0p9sa/k6cBj2VXxPmU3ab9TS3+2JaUqJzC223JPGkYlmYGhkj5gXF6ACSuvJnsOmMbdl6mKHGjuqxjAe9hdX1frm+VLZUUaywM+ARYnWxNsSrYlHym2QoiiKMj5UZCsHLcm2lu5kW8sXgcOBb2LjX6c0R8mLHV2edCFSXzSDNI6oE4yfGoqHjkrM4DzZIqZknuz5tBcNXxkiIpUIGwuTH/BPmHtUJr2rVHaUgRTwfeCIRGuiNdGaeKrcBimKoiilQUW2ctCTbGvZnmxr+RxWbH8Z2FZmkyaMRd66XS/EYpuKLjTGmMjIrz0a0z/UudCIGATX4bGUnuyOEZ94H0FFZEth0nNm/sJT+jOTMS67E7gGWJhoTVymYSGKoihTjyJfdlOUg5NkW8tm4P81L1l6DfAe4HLghPJaVVoWybrBZ+LxoiN5TOtn52DF9JEjhpjUUJhGzmNoVJJSerJns3n67vL4NdH+yEChnhY5LFVVmu+77x82Ad8BfpBoTRw0YUqKoigHIyqyFaWAZFvLAHAjcGPzkqVnYsX2W5gCb36O8TZGbo/HR3w9EaCph12DFfVFwkUGhoS1LxIW2SWrPxrZNRNjcoUjm4TxG+LRyI6RptempnuYrn5EqktlzwSwFvg6cHOiNVG0/BVFUZSpxaQXDYoykSTbWu5LtrW8HXgJ8G2gq8wmjYv5sr12bTxWdHzrmV2mJxOvGxkbbYYGF8EXwuNrl8yT7eFHPPwxw3RMY7xIr0x42F8cr/MP2C8/Pgi8Gzg60Zr4kQpsRVGUgwf1ZCvKHpBsa1kHfKZ5ydIvAq3AZcBR5bVq76mnb1Z7tK5oDPOsTlLZaNW0wnRjBoYa474X1tilCxcBqGBwR4rqOaMt96fFFxRLf8A/Ye7R6b9sXF5VWUpzxsN67EghP0u0JqbCyCeKoijKPqCebEXZC5JtLb3JtpbrgGOAc7FialLE1sbIpiP4czsjkaZiy2d1mkzOi48Q2Zj0kJj2RYa84CKlbaTX0j12Oca8aabIh4TWmAXNp6QyRb9guR/pw14LZwOHJVoTX1CBrSiKcnCjnmxF2QeSbS0GuAe4p3nJ0o8BbwLeC5wHxMtp22gcLh3tCM0ZmFds+cwuTFddZERIhjGZoa8tZsWExGzpYrIBmtgxuI1RHdmWuNdB2h+WyeB5R/bHfXbbdbLkGOB+4CbgtkRrone/W6AoiqIcsKjIVpRxkmxrSQG3Arc2L1naCJyPHZ3kDA6gt0XHyfrtWyKRKkSKKtnaVJwukSJD6WWG4rSNDDuektYfs9jCGo4fM4+pjvZIOj0ivWZgumB6+xCpKaVNo7AWuBn4eaI1kdwP+1MURVEmISqyFaWEJNtadgE/Bn7cvGTpIcC7sIL75LIaBiz21vU/H49theLu4og0Fh8Kz+SGgp1zMiwqo6T1xzzad7s9vyEuXudIkf2wvyg+zf/H2u5IZKKGXOwCfgvclGhN/H2C9qEoiqJMIVRkK8oEkWxreRH7JclvNi9ZeixwAfAG4OWUcGSOPeVYWc/qeHzUuGePooN3AP7Q0Hg5MeGRSUotsnfrhfYb4kU/WnO/f+Kso9P3dyyrKlmxGmAlcDfwJ+ChRGsiO/YqiqIoipJHRbai7AeSbS1PA18CvtS8ZGkDcA7wejcduj9sWOhtqbo9HusfNUOksejQfuAPid8cw0R2rEjmfWY2mxp2l8evj88vlr7aNB/+2VS6Y9n4RhjZCvwZK6r/lGhNTNkvfyqKoigTj4psRdnPJNtaOoHb3ETzkqVHkxfcZwIT8lGVJnqa1sWaRsZagP2kujfyk+qO2uBPTvKfUi/16CIz2TZrt5kqI7MMdAoME+Q+XuToVGxvRxjJAg+R91avTLQmyj1KiaIoijJFUJGtKGUm2dbyDPAM8L3mJUvjwOlYwf067GfdSxADYUyczPzN0UjRcJHaFJ3piukjYrKNyQwAlQAeXiYnvgEQ8TOU2JNdwWA1xnQhUj9mxqhsJGtGeL1rBhrBpHoRqS22mmM9VlDfDfwl0ZroHpfRiqIoijIKKrIV5QAi2daSBv7qps83L1laBZwEnBqajmEvRy2Zz/bNIszt9byi3uKmHnYOVBYJFzGDvTiRHfFi/Rnsd2xE/CwlFtkAMdJbM1SMKbJNdbRTujMj0h/2F0Wm+Y+s7Y5ETnRJncAKYDmwDFieaE2sL7XNiqIoilIMFdmKcgDjhgd8yE0ANC9ZWoMdrSQQ3adgvz45qvA+1lu/NQ1NuVFGFpnZZXrS8foRotmYwT5gBkBEYoMZyYU92VX7elyjUUN/Vycjv+wexp8W872RInvX/bkTO05P3X/3H2prvoYV1s9r+IeiKIpSLlRkK8okI9nW0gf83U0ANC9ZWge8FBte8hLgSPfbDMQWe+u6X4xFNyGysNg2Z3XSn4nVjOw16KdSwd+YxAaz5ASGPNklp4Fdqc5RRzkBYLM/vWIDL/b/CFgDrAZWJ9ta2u3iT/K1iTBMURRFUfYSFdmKMgVItrX0YL8+eH84vXnJ0iiwsJ6+hR3R6CHAsdjRTA4B5mO//lg1u9NkstUVMwq3a8zAQPA/6sUHsy5cxPMmRGSn57NxZ5LDHwHagY1AElgXTJvPOsnGlL9nAvauKIqiKCVERbaiTGGSbS1ZYC20rAVIFMmz+ObF04AZkVx6ph+pmA00uqkBsjls6EWNJ5F+Dy8DIOIHYSRewWSAQWCgYOoFekJTNzZmejOwBdhyztlrd55T+iJQFEVRlLIgxmjIoqIoiqIoiqKUkr0aoUBRFEVRFEVRlN2jIltRFEVRFEVRSoyKbEVRFEVRFEUpMSqyFUVRFEVRFKXEqMhWFEVRFEVRlBKjIltRFEVRFEVRSoyKbEVRFEVRFEUpMSqyFUVRFEVRFKXEqMhWFGVKIyJXiMhqEVklIo+Lm7dPVAAAGR5JREFUyCtc+idFpHoP1t+jfHuwnTNF5K6CtJtE5Py93E5SRGbsJs8FIvK0iNwrIqeKyPf2xeYi271ERK4txbYURVGmOiqyFUWZUhSI6meBdwEvBR4FPg5sdFk/CYwqnkPi9JPAB/ZWXIbFrYicCRy/1wdTHA/4k4g8JiKvLrLfmcAtwB+MMWcZY5YbYy4r0b4PWErVGCrY5hWuYfa4iORC/y/blwbSXuy35Mcyyn7uE5FTC9L+1x3j8yLSFTrmf5poexRlqjElRXZBZfi4iCyZ4P29ZT/s48w9qeRG8zS59G2uPJ4SkQ9PgI0jKuwieYY9PETkDyLSUGI7LhERX0ROCKU9KSLNJdh2TkTSIpJyU8dYXkIR+Y+C+X/s436/ICIXheb/LCLf3s06uz0fk43Qvf2kiNxZeO2IyCuBNwEvNcacAFwDbDLGDBpjPmSM+YcxpkNELgPmAfeKyL1u3RtEpEdE1orIVcaY5cDzLt/ngbe5fL2h/Z0vIje5/xc4u54QkWeAS0Li9kzgHcCrnf1pEUkAbwHe69ZPishVIrJSRBIicoxLb3Ln+zER+SFQBTxnjDkZWCgij7pt/lBEIsCNgAEuFZFvhD3oInKliNzoro0XXDkEx3JRkW0hIu8XkWdF5H7gVXt4niJ7kq/EjNlo2heMMdcYY04yxpwEpIL/xpg9ejMwjnIo+bHsKcaYt7vj/RDwt9Ax71PdpSgHM7sV2SJyiIj8TkSecw+f74pIfJS880Tkt3uwzX0WVu4h8dlR0ttF5HFAgJ3AfOA44Gr3gKopckxbXcv9xVGO9ycicptbtlVEForIx0TkYne8txljfg+cMdoxiX1lm5PhXoGL3bL7nCBMu6ldRN4cOqZuEWkHbgN+KCLHuWUrRWSL29ZzInJ7sAw4RkTmFTFlJzAbaAK+IiKzRWSGe7h/LLApZPfpImJCtqXcQx4RaXbzQ40Zt21E5GERGZDhDZ1g258EqkVkpog8AswFFu/mnF8nIjtFJBPaZ0JEtobyREVku4h81SW9CFwjIveLyDrgCOBuJyxOEpE3jnW9uuu8XUKv5sUK5BTQASwwxlQZY+YZYy5z53GbWO/pGhG51l0PV4rIWW5djDH/JCGRVrDP20RkUER6RWSdWE9Z1C1+D/DmUPajgNe5cnhKRL4sIhWhbb3CHXN4+wtF5Nfu/5Uicldw/EGZFOQvqadOdiOQ95BA6CzCXs//UrB8LrDdGDPo5m8DZokVie0i8hGX/hXsucwBORF5OfY6jALXA68RkQ8A52LP99eAO9y6ERF5REQeA74AVLr0a4HHgS1AFzDflXEz8DHy1/knsHXU64DfA7eLSNIt226MeSlwA/ba+QvwDPZcvhVIADOA14rI08C7gVc5UZTDCvYGZ8cO4L+CgnHX3bnAOwEfOA/4ortHbga+hxV2H3Hber+I/Ar4ATDolh8H1InI31wdtFJc41+smL9XRH7p7ERE7hCRFWLfLHwkbIuIXCO2QfKwiAR1x2yx9fETbgq2XbQBENpesUbTu9398aSIfI0iiG28LnN5fiQiUizfGJwhIv9w9cr5Y5TDaA2YG0RkuSufq8Y4lnNF5CFX3reKSK1LH61hViO2MbVM7LPvrS69SkR+7eqpW7ANtj3G7Ssorx+IJebO8ekuzzdCx/I6ydfXPxanHUTkRVcHPeZsOWovy11RJgfGmFEn7IPgUeD9bj4C/BT4RpG80bG2VaoJuBL47FjpQC8wAPQDM4EY8Bz2oXq0O6aHgQ8D9wNZYBmwEvtwvMxt5xisWNsM/A04xqXfBHwHuBf4FnAJcG1o2Q1u2QvAa4DfYB9qG0P2ngs8BPQA3cB/uPLtcnlXYx/Wf8d64zZjRUUG66lLA9uBx4C/YB+Mm7EP8Hbg1ILyuQQrAja76QXg5cBLXDmtcmVyQqg878Y+bJ8DPuzSzwTuApqBJ902P+CW3Qec6rbTCSx3x3GVW36Zszvh0n8BJIEZbvln3DafBD7t0pqBp4Fn3boDQC1WaCRDx/dG4EFgrTvW64F1wCbgfLfN093vDVgB8iRWPEVD1801wBNYUbHCldWM0H56wza7tHpXhhngKXfsm4ANrvxS7tw95tbvdTZ0u7zbgf8B/uiOcZs7zqw7noux10TWbesit6//dfv/lZvPufN6qEvfDPwa+IdLvwIYcMsKj+ES3DUcSrsJOL+E925v6P/NwBXj3MbHgOtD859z5yzlzu/1wOuBpViPdBp7v12C9fQG5ysJ/B82lKTPlX2Pu57ucst/Dqxy+xnA1iGPY6/j1S59Hfn75WfA/wJ3he6n60PzPcD7yNcXv3TbfdLt93x3Hdziro9/w9YpCWfjz7H3bMbNr8GK8UewdUinu75uI3/PGjd/BbZxcC323nrI7bODfD35DLbO+447lgbsPfhZrOiudMdxJLA8VDf0AYeFzsl091vljq3JzRvgze7/14H/dP9vAT4Zet7UA8cCdwIxl349cHGRayNJvi6Zh73/ZmIbTn8F3lZknemh/z8PbNrdtRe6P27FOquOA54vVg5j2R8qnwi2/jyhyLHMAB4Aatz854EvhPJ9wv3/OPAT9/8rwEXuf3DuaoBPAze69BOwdcqpoxzvmbjrtcj5FGy9c15oW09hn2srsM/camxo1hEuz/8A/+r+vwhcGnou/KBU9YxOOh1I09gL4RzggYK0adgHWDX2YXWrq0D+ihNeLl81VlyuwlacjwQ3c1CBkBdQP8Y+rP4MVLk8H8YK3yewD4Zql34luxfZOTftdDfzhcA3sBX7AFZADWAfRL5LN9jXwX2hNB/7wMkCZ7i0TW7bwbYM9uGdw4qnTGh7ObfuFrfdlFtnp0u/1+Xpw4q0p916WVfGK52tfmh/v3ZlbUJTF/ahnXF5feyDud+V63KsyMq6ZTvceitDx9KDrSTbsQ/nvtD2HnPbmgdsdf9NaF83u23vcmXgu99g28Ydc7f7vyO0vsF6yJJufpc75ow7zn6XnsFeCwb7sF7m8l1DvmGSdvu9w5V1cC6yzpZg24FNwf77sAIjOGc+VuA8E9puqqDMw1MutK+gTDpDx7+T/LVUuO5zWIHXV2DXxpB96dD6q0Ppy9z5yrmyDc5vluHXYXfBuUgXHP9G8vdMxp3jPuz18BDwJ7fOf42rstm9QF6GrS+uCqVf7NKewIqgXmAhtmHZib02D8Xeuz1YwRPFNvT6XJ5HsA3GHux1dl/ovF/vtrnBHePDLt+vXTn+ya2zGlsPHevS7yEveHdiPehdWEH6Zex1tpThIvubwINufhO2UfV7t+3TyNeLb8DWhRkg7uz7mNvvaW7+KWy9+1WssAmE7v+5Y7kf+Khb52ysyM4BX8IK5cOxjYQn3TF/D/iq28YGrDBbDqx3x/e4S78G+JE7FwmX3u/WOxO4t0i9/ISbuoDTXPogIO7/heTF4TagomAb/4q9Px930zPAlUWuryR5YfpW4GehZR8Evl1knXdgr48E9l5asifXr5u/CXhvaL6nWDmMZb87ryux1/g24F1FjuVN2Po8WP8p4KehfPPd/1cA/+f+L3fnNlhnA/bavQM4O2TbSvZOZF+AbWCucsf02dCyL2DrqKChcArw19Dy1wO/cf9fBGa7/68C7p5osaOTTuWYdhcucjy2VTqEMabb3bAvcUmvBFqNMWcXrPtxYJexcZFXuxuuGEcC1xljjsc+EN/h0m83xrzMGHMiVnx+cDe2DjMTW6l/CbjTGHML9mGUwnqbYtiH10uxHoJAbJyGrfzBPrh97KtgD/tgAajDVmxgH+YD5AVM4MH1XdozLj3qthMFKrDCMYutFD1ny9PGmGPdep3kPakRrOfoKrfd00Lb7cdWpJuxD96bsALCB/5ojKl25boW+zD3nd2B+G/EesC63X6iwHTyD8VNWI/FiVhP1DpsIyt83fjAWW59D1sxC9Zrkg3li7rjCMqwDStwDPBPWI/Vi8621VgvYNrtdwtWmHeS90Ce58r+bdgQlIyb7saKlKVYsYizLY0VCzHsA3WZ22YvtpEwy+WdifX8Tnd2Qv6tQXA8BvugMdhz8S63D5wNPViRFAj3erfMc/NgxRnYc/Aq7PUReKODV9ZBWMl3nQ0xrBC4waWfiBVPgr22Axs8lw/sdVSH9WQFjZVPYM/5c9jzL9gH8YArlxTwQ+w93ow9T0uBd0jxMKS9wr0qPwd7TSIi52LrgZcDJwGniMgZInI81ut6tqsHLsdeD6uw98EKrDf2e9i3FBF3nMud3SuA72PfrmzAerNXYcVG0Hg6GrjdldNOrMiNYAXBAPY67seGQm1ydkewb7mCcxXH3ssdWG/3f2IFU+Gr+G3APBE51u3vcHfMA27bs7De42+6ZZD3Pn/P2dvobJqDDQU5H3vu5onIIViPYr2z54shW4thQv9XAueLSHAf1GPrqvdir+dzsCErZ2DP0Rbs9XeqO/6A4LoOOnu+FnilO3+PkQ+tyRhjgv3nsPXDaAhws8nHBh9tjLlyjPzBOmNnEKnEXtvnG2MWY509lWOvNYLB0P/wPvsK0kfYLyKHYRs857jn5NJR9i/APaH1jzPGhJ+HgQ3hchTgHaF1DjXGPO2Whc/7HiO2L821wNudvTcW2LsI25AKrqHdnYNidivKlGJ3IlsofkOG0+8xxuwskud0rDcIY8yT2IdbMdYZYx53/1dgH44Ai1zcXwJb0e9tz3zD8Jv8LOxDI2gcdBljXiiyXp1b9263fuAVvCBkY+CtXOPS2rFl+TO3vrj5IM4sEONdbltZNwVl8gLwiIj8DftQrCEvygXbwLjCpc1x6wh5r3IMKw7fjw1PiQAvc3GjCeyD7n2h8liOFVIV2NCLKjctdGknY8VHg7OzE/uAv9j97nLbSbnpC9iK8vmCspxD3gtcSb7B1ocVp7hynOXWDbzC12EbOY+4PDvc7wqX51BXvhms8P+ms7PCbVuwD//gAeBjPXKfcNvvcXm2u7I+35VZ8DYiuEY2u/15Ll9wv2TJi6BVbr/B/RDBXgMXYD3/EppizmawDUcfe2790L4i7vcQ7LkB22iIumUnYGN5cds73207HN8sbt9hGkPLXod907QQK27nueOpwr596sCK94w7vndjw8TWYMt+X6kSG7u/A9uIuceln+umx7Bi7xhn19nAb40x2wFcPRO8oTgUK+7qsHWNYEMulmHL/iFs434B9rxdghWHr3fbAHsNHuvKYwDb4D4Mez57sedoHrbBWYlt0An2fsD9riffGJoLfEZEnsRez12hY+/BlvlFwH+77dZiz+FPsB7iFPBqt58Kt/67gD+4/UZdOQX3w3PAf2IbdDPdbwVWXD+M7VS5AyuMIR/WADbG/+8h+05z23rQbet27P32Lqw3+iG3zZWu3DcZY3xsvTJa5756rKOl38UKnzZKvjB/AS4F2xgTkWkubagBICLTRWRhkXV7yDeMH8HG1c9wjbp3Y737YYL6YbuLcZ6QkUIY3f5p2Pqqy8WknxdaJ3wsDwOvEpGXuPWr9yCG+U/AJ4IYcxE52aU/QL6j7SJsfbKnVOHekopIHXmHGCJyIfZ6PhO4zp23p4AjRSSoKy9i5DlQlCnN7kT2aqxYGcLdPAuw3lEY3mIflnUPbQh7AsIt2puw8VuLsV7cvfEw5LBel5OBp53H4gfY4/0rLq5XRGIuv+fsfdjtf5vzAHe4ZRHyxxu2N7w/sBVQBfbB9rjbVyCwwAqhHcDL3HaqyTcG3oN9rZjDCudu7APYc7YvcfZEsA9/sA/q87Ge1whWrP/KbfNoN/8yZ9NyrCj33XYD4TmXfAjFLmzlvhzrPdvi8gVe6iOxD/6gg12M4ecl8PQGHvMu8uI3XE6FDTfj9rUAK7ivwQrd5aFtNboyEpf3nW7fr8VW/hXOxnc4u2rd+gFN5BuHR5OP1TfY0KQglOIN5L31x7vfnVgvU2CLwXpIfayH7wzy5/iD2IZFEPIShP74ofUhL/7OduuGPeFZrFc/KKdPYc9Lzu336NB2gpClP7j5FLbcN5BvhPjY0KTgrUotViDuwt6/f8R6x7Pk77/Nbj7wiv6R8XucUsZ20FuIFWpBp0XBhioEXreXGGN+yuiNfLD3x2XYGFODFRUnYh/y/4N9xf5d7D2yxRgTx9ZnD2FFvGAbx5dhr6e1xph3unIJzsHVWA/2HGzfizuwginoMBnBhmesI98gvAFb7xQKzzuBt2MbkJ/D3stBY+dX2HL+gGtQXIK9fk7FXps3YxsPy7Ai5bPO9ve6t3SfxdY3v8G+2cIYc6axo6P0YEdaeRP2XAedJc8GvmRsB9IB7LVwudvva4wxp2DFbsyVcRorrC/HNmJaReRhrCNhtGfA3UBURFa5snx4lHxhLgfOcs6BFcDxxpincI0Jt617sPVWIT8C/igi9xpjNgH/jr3unwBWGmN+F85sjOnEeq8T2HO7bA/s22tGs98Y8wS2Ybka6xV+cJRj2Ya9Jn7l1n+Y0d9OBFyNPXerXKPvapd+A/bZtwob5//oXhzHDuy1+CT2TeMjAK7x8GVsv5012Ldg3zHG9GPrw9vd+RzElreiHDyMFUuCfRAtJ99JI4K9Sb7l5i8h1FmK4THZnwNucP+Pwz64isVkPxla/7PkY9W2YwVXDFsp3eTSr2TPYrKDkI0E9mEbdBqb4/4PYsX7neRjWR8nH6tZi32QGbf8KfLhAc+Rj18NOgn5WE+awYrZDPk42SDGOYcVTzFXrkEc7oCbZpEPPbk9tL0fY8NyAiH4JvIxtA9ivYwG61H7d/JhAVXYEVay2FfzHW6de7Fi4kryoj4Iw9iCfWX9I6xwCcpxNVa4PYGNsQ/y92Gvgywu7j50DJvIx0NvwDYIgpjgJPbVo49947Ee2ygYdPb8xNmXIR8S81W3vZ+69QOvZjf2Gklh+wgMurQ1IVuCsJS021/arduPfcsReNF/Rd4z3+6OL0s+ZjeIZQ7+p12ZBftY65YHMf+BwA62H8RK3+d+B8kL6MCm4PoMGj9LsF7HoMwvCh3X067cgv33OhtWu7zrXfq1IZu/7/a5zdn3dayH0mCFRhLrKQ6uy6+5++pu4PR9jU1jeEz2ya5MY1jv7CNArVs2H3svHI8N/wg6y013x/d74H0u7THy8ciXYxsYu7DXyxFYz3UfVggvw16fJ2IbRztc+dyIiwPGXmvbQsd8F/BCwXFciK0rVmGFYBBn/H5n7/3YOufaMcoi6fL/3M2f7tZ9ECtYgk50PwM+7v4vxNY9M9x0p7PhH8Ail+fLuI6Dbn4NcEhh+RfYch+jxOXqpJNOOum079PuM1iP252ucl/rHtDBA+kSRhfZNcBv3UPgZvcwOtItS7J7kX2pezDe5/Z5k0u/ktFFdrt7+D3n1ttOXrhswb7C/T93HAMuX/Ca3seK0sDjF4giHzeKhZtPkhdOCaxIe8BtLxiNIxDmTzsbgo5nA+Q7Bd6PfegHHs6wGLwb68m/siA9i/XIgBUTQXo/+djwHVgxncGKtL9ihfFOrPgYCNnzQexDOBB+j2I9lr0u/6aCfRusADmVvOBLk2+cXB0qp6AzaLD+18iL7CvcOQg6Zx6B9VL3YjvhdWCvmztwwoC8yE5jr53b3f9nyHdcCzo3fsbt5ynynQ19rNjqJu+J9F0Z/cRta6M7ztXYa+kb5BsJ213+TdjrqN3Nd2Cv26ARFz6fwTnZSL4xFdgTdGbsI9/BLLjmUm6doBE2iBXSK93+nyQvwINzH+TrctsOwpw2kO/sup28cO4mL/67yXu8s9hG4h1YAZshP6pOyUS2m7+TvFi+HHs/JbDe5mBEglZ3vE+QrwOasdf1Kqxn+dDQNoNr7DWhtMOc7U+4ayIYmeEmCkZPwYr+7uCYJ6zitXXa88BR+62yV5Gtk0466bRfp6B3d8lxcXAxY8yAiBzhHoZHGWPSE7LDcSAivcaYYNzRGHYkky4ROQvrIa03xvSNdUyT6Xj3BBG5EvtQ/ma5bVHKg9gP2XzHGDPiq4JTlYPxmBVFUZSJYSJ79FZjB9OPYcNOLj1QBWcgsB0zgOdFJIhXv8oYE8QcjnVMk+Z49wSz+977yhRG7BdML8V1kjoYOBiPWVEURZk4JsyTrSiKoiiKoigHK7v9rLqiKIqiKIqiKHuHimxFURRFURRFKTEqshVFURRFURSlxKjIVhRFURRFUZQSoyJbURRFURRFUUqMimxFURRFURRFKTEqshVFURRFURSlxKjIVhRFURRFUZQSoyJbURRFURRFUUqMimxFURRFURRFKTEqshVFURRFURSlxKjIVhRFURRFUZQSoyJbURRFURRFUUqMimxFURRFURRFKTEqshVFURRFURSlxKjIVhRFURRFUZQSoyJbURRFURRFUUqMimxFURRFURRFKTEqshVFURRFURSlxKjIVhRFURRFUZQSoyJbURRFURRFUUqMimxFURRFURRFKTEqshVFURRFURSlxKjIVhRFURRFUZQSoyJbURRFURRFUUqMimxFURRFURRFKTEqshVFURRFURSlxKjIVhRFURRFUZQSoyJbURRFURRFUUqMimxFURRFURRFKTEqshVFURRFURSlxPx/7W5kfX0PWc8AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 576x576 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "status_count = df_species['status'].value_counts().sort_index()\n",
    "plt.rcParams['figure.figsize'] = [8, 8]\n",
    "status_count.plot(kind='pie', autopct='%1.1f%%');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Downloading petition data\n",
    "\n",
    "One of the ways that species get listed under the ESA is by petition from the public. ECOS makes these petitions available in a nice JSON format for us."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "750 petition records downloaded\n"
     ]
    }
   ],
   "source": [
    "url = 'https://ecos.fws.gov/ecp/report/table/petitions-received.json'\n",
    "params = {'active': 'any'}\n",
    "r = requests.get(url, params).json()\n",
    "print('{} petition records downloaded'.format(len(r['data'])))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Fields available in petition data:\n",
      "Petition Title; Date Received by the FWS; Where the species is believed to or known to occur; Petitioner Name; Requested Action; Petition Finding(s); Active; Petitions Document(s)\n"
     ]
    }
   ],
   "source": [
    "columnheaders = [x['title'] for x in r['metadata']['columns']]\n",
    "print('Fields available in petition data:\\n{}'.format(\"; \".join(columnheaders)))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's put the petition data into a DataFrame as well."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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>Petition Title</th>\n",
       "      <th>Date Received by the FWS</th>\n",
       "      <th>Where the species is believed to or known to occur</th>\n",
       "      <th>Petitioner Name</th>\n",
       "      <th>Requested Action</th>\n",
       "      <th>Petition Finding(s)</th>\n",
       "      <th>Active</th>\n",
       "      <th>Petitions Document(s)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Petition to List the Dunes Sagebrush Lizard as...</td>\n",
       "      <td>06/01/2018</td>\n",
       "      <td>NM, TX</td>\n",
       "      <td>[Center for Biological Diversity, Defenders of...</td>\n",
       "      <td>[Listing: Threatened or Endangered, APA: Desig...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to List the Dunes Sagebrus...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Petition to List U.S. Populations of Lake Stur...</td>\n",
       "      <td>05/23/2018</td>\n",
       "      <td>AL, AR, GA, IA, IL, IN, KS, KY, LA, MI, MN, MO...</td>\n",
       "      <td>[Center for Biological Diversity]</td>\n",
       "      <td>[Listing: Threatened or Endangered, APA: Desig...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to List U.S. Populations o...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Petition to List the Siskiyou Mountains Salama...</td>\n",
       "      <td>03/12/2018</td>\n",
       "      <td>CA, OR</td>\n",
       "      <td>[Cascadia Wildlands, Center for Biological Div...</td>\n",
       "      <td>[Listing: Threatened or Endangered, APA: Desig...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to List the Siskiyou Mount...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Yellowstone Bison DPS</td>\n",
       "      <td>03/05/2018</td>\n",
       "      <td>CO, IA, MT, OK, SD, UT, WY, Canada, Mexico, Un...</td>\n",
       "      <td>[James Horsley]</td>\n",
       "      <td>[Listing: Threatened, APA: Designate Critical ...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Yellowstone Bison Petition Volume 2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>PETITION TO LIST THE REFUGIO MANZANITA (ARCTOS...</td>\n",
       "      <td>11/30/2017</td>\n",
       "      <td>CA</td>\n",
       "      <td>[Los Padres ForestWatch, The Chaparral Institu...</td>\n",
       "      <td>[Listing: Endangered, APA: Designate Critical ...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'PETITION TO LIST THE REFUGIO MANZAN...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                      Petition Title Date Received by the FWS  \\\n",
       "0  Petition to List the Dunes Sagebrush Lizard as...               06/01/2018   \n",
       "1  Petition to List U.S. Populations of Lake Stur...               05/23/2018   \n",
       "2  Petition to List the Siskiyou Mountains Salama...               03/12/2018   \n",
       "3                              Yellowstone Bison DPS               03/05/2018   \n",
       "4  PETITION TO LIST THE REFUGIO MANZANITA (ARCTOS...               11/30/2017   \n",
       "\n",
       "  Where the species is believed to or known to occur  \\\n",
       "0                                             NM, TX   \n",
       "1  AL, AR, GA, IA, IL, IN, KS, KY, LA, MI, MN, MO...   \n",
       "2                                             CA, OR   \n",
       "3  CO, IA, MT, OK, SD, UT, WY, Canada, Mexico, Un...   \n",
       "4                                                 CA   \n",
       "\n",
       "                                     Petitioner Name  \\\n",
       "0  [Center for Biological Diversity, Defenders of...   \n",
       "1                  [Center for Biological Diversity]   \n",
       "2  [Cascadia Wildlands, Center for Biological Div...   \n",
       "3                                    [James Horsley]   \n",
       "4  [Los Padres ForestWatch, The Chaparral Institu...   \n",
       "\n",
       "                                    Requested Action  \\\n",
       "0  [Listing: Threatened or Endangered, APA: Desig...   \n",
       "1  [Listing: Threatened or Endangered, APA: Desig...   \n",
       "2  [Listing: Threatened or Endangered, APA: Desig...   \n",
       "3  [Listing: Threatened, APA: Designate Critical ...   \n",
       "4  [Listing: Endangered, APA: Designate Critical ...   \n",
       "\n",
       "                                 Petition Finding(s) Active  \\\n",
       "0  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "1  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "2  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "3  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "4  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "\n",
       "                               Petitions Document(s)  \n",
       "0  [{'name': 'Petition to List the Dunes Sagebrus...  \n",
       "1  [{'name': 'Petition to List U.S. Populations o...  \n",
       "2  [{'name': 'Petition to List the Siskiyou Mount...  \n",
       "3  [{'name': 'Yellowstone Bison Petition Volume 2...  \n",
       "4  [{'name': 'PETITION TO LIST THE REFUGIO MANZAN...  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_petitions = pd.DataFrame(r['data'], columns=columnheaders)\n",
    "df_petitions.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The `Date Received by the FWS` field contains the date that the petition was submitted. We can add a column to the petition data that contains the year component of this date."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Petition Title</th>\n",
       "      <th>Date Received by the FWS</th>\n",
       "      <th>Where the species is believed to or known to occur</th>\n",
       "      <th>Petitioner Name</th>\n",
       "      <th>Requested Action</th>\n",
       "      <th>Petition Finding(s)</th>\n",
       "      <th>Active</th>\n",
       "      <th>Petitions Document(s)</th>\n",
       "      <th>year_received</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Petition to List the Dunes Sagebrush Lizard as...</td>\n",
       "      <td>06/01/2018</td>\n",
       "      <td>NM, TX</td>\n",
       "      <td>[Center for Biological Diversity, Defenders of...</td>\n",
       "      <td>[Listing: Threatened or Endangered, APA: Desig...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to List the Dunes Sagebrus...</td>\n",
       "      <td>2018</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Petition to List U.S. Populations of Lake Stur...</td>\n",
       "      <td>05/23/2018</td>\n",
       "      <td>AL, AR, GA, IA, IL, IN, KS, KY, LA, MI, MN, MO...</td>\n",
       "      <td>[Center for Biological Diversity]</td>\n",
       "      <td>[Listing: Threatened or Endangered, APA: Desig...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to List U.S. Populations o...</td>\n",
       "      <td>2018</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Petition to List the Siskiyou Mountains Salama...</td>\n",
       "      <td>03/12/2018</td>\n",
       "      <td>CA, OR</td>\n",
       "      <td>[Cascadia Wildlands, Center for Biological Div...</td>\n",
       "      <td>[Listing: Threatened or Endangered, APA: Desig...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to List the Siskiyou Mount...</td>\n",
       "      <td>2018</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Yellowstone Bison DPS</td>\n",
       "      <td>03/05/2018</td>\n",
       "      <td>CO, IA, MT, OK, SD, UT, WY, Canada, Mexico, Un...</td>\n",
       "      <td>[James Horsley]</td>\n",
       "      <td>[Listing: Threatened, APA: Designate Critical ...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Yellowstone Bison Petition Volume 2...</td>\n",
       "      <td>2018</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>PETITION TO LIST THE REFUGIO MANZANITA (ARCTOS...</td>\n",
       "      <td>11/30/2017</td>\n",
       "      <td>CA</td>\n",
       "      <td>[Los Padres ForestWatch, The Chaparral Institu...</td>\n",
       "      <td>[Listing: Endangered, APA: Designate Critical ...</td>\n",
       "      <td>[{'name': 'Petition findings not yet made', 'u...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'PETITION TO LIST THE REFUGIO MANZAN...</td>\n",
       "      <td>2017</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                      Petition Title Date Received by the FWS  \\\n",
       "0  Petition to List the Dunes Sagebrush Lizard as...               06/01/2018   \n",
       "1  Petition to List U.S. Populations of Lake Stur...               05/23/2018   \n",
       "2  Petition to List the Siskiyou Mountains Salama...               03/12/2018   \n",
       "3                              Yellowstone Bison DPS               03/05/2018   \n",
       "4  PETITION TO LIST THE REFUGIO MANZANITA (ARCTOS...               11/30/2017   \n",
       "\n",
       "  Where the species is believed to or known to occur  \\\n",
       "0                                             NM, TX   \n",
       "1  AL, AR, GA, IA, IL, IN, KS, KY, LA, MI, MN, MO...   \n",
       "2                                             CA, OR   \n",
       "3  CO, IA, MT, OK, SD, UT, WY, Canada, Mexico, Un...   \n",
       "4                                                 CA   \n",
       "\n",
       "                                     Petitioner Name  \\\n",
       "0  [Center for Biological Diversity, Defenders of...   \n",
       "1                  [Center for Biological Diversity]   \n",
       "2  [Cascadia Wildlands, Center for Biological Div...   \n",
       "3                                    [James Horsley]   \n",
       "4  [Los Padres ForestWatch, The Chaparral Institu...   \n",
       "\n",
       "                                    Requested Action  \\\n",
       "0  [Listing: Threatened or Endangered, APA: Desig...   \n",
       "1  [Listing: Threatened or Endangered, APA: Desig...   \n",
       "2  [Listing: Threatened or Endangered, APA: Desig...   \n",
       "3  [Listing: Threatened, APA: Designate Critical ...   \n",
       "4  [Listing: Endangered, APA: Designate Critical ...   \n",
       "\n",
       "                                 Petition Finding(s) Active  \\\n",
       "0  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "1  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "2  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "3  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "4  [{'name': 'Petition findings not yet made', 'u...    Yes   \n",
       "\n",
       "                               Petitions Document(s) year_received  \n",
       "0  [{'name': 'Petition to List the Dunes Sagebrus...          2018  \n",
       "1  [{'name': 'Petition to List U.S. Populations o...          2018  \n",
       "2  [{'name': 'Petition to List the Siskiyou Mount...          2018  \n",
       "3  [{'name': 'Yellowstone Bison Petition Volume 2...          2018  \n",
       "4  [{'name': 'PETITION TO LIST THE REFUGIO MANZAN...          2017  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_petitions['year_received'] = df_petitions.apply(lambda x: x['Date Received by the FWS'][-4:], axis=1)\n",
    "df_petitions.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Some petitions request that multiple species be listed under the ESA. We can find these exceptions if we look in the `Petition Title` field. We can use a function to check if the Petition Title refers to multiple species. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "def parseSpeciesNumber(petitionTitle):\n",
    "    nspecies = re.findall('\\d+', petitionTitle)\n",
    "    if len(nspecies) > 0:\n",
    "        return int(nspecies[0])\n",
    "    else:\n",
    "        return 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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>Petition Title</th>\n",
       "      <th>Date Received by the FWS</th>\n",
       "      <th>Where the species is believed to or known to occur</th>\n",
       "      <th>Petitioner Name</th>\n",
       "      <th>Requested Action</th>\n",
       "      <th>Petition Finding(s)</th>\n",
       "      <th>Active</th>\n",
       "      <th>Petitions Document(s)</th>\n",
       "      <th>year_received</th>\n",
       "      <th>n_species</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>171</th>\n",
       "      <td>475 Southwestern species; List T/E</td>\n",
       "      <td>06/25/2007</td>\n",
       "      <td>AL, AR, AZ, CA, CO, IL, KS, KY, LA, MT, NC, NM...</td>\n",
       "      <td>[Forest Guardians , Nicole J. Rosmarino, Conse...</td>\n",
       "      <td>[Listing: Threatened or Endangered]</td>\n",
       "      <td>[{'name': '12 Month Not Warranted on 09/27/200...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to list 475 Southwest spec...</td>\n",
       "      <td>2007</td>\n",
       "      <td>475</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>101</th>\n",
       "      <td>404 Southeast aquatic species; list w/CH</td>\n",
       "      <td>04/20/2010</td>\n",
       "      <td>AL, AR, CO, CT, DC, DE, FL, GA, IA, IL, IN, KS...</td>\n",
       "      <td>[Alabama Rivers Alliance, Center for Biologica...</td>\n",
       "      <td>[Listing: Threatened or Endangered, APA: Desig...</td>\n",
       "      <td>[{'name': '12 Month Warranted But Precluded on...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Withdrawal of Three crayfish', 'url...</td>\n",
       "      <td>2010</td>\n",
       "      <td>404</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>214</th>\n",
       "      <td>List 225 Candidate Species with Critical Habitat</td>\n",
       "      <td>05/11/2004</td>\n",
       "      <td>AL, American Samoa, AR, AZ, CA, CO, DE, FL, GA...</td>\n",
       "      <td>[Barbara Kingsolver, , Biodiversity Conservati...</td>\n",
       "      <td>[Listing: Threatened or Endangered, Listing: T...</td>\n",
       "      <td>[{'name': '12 Month Not Warranted on 05/11/200...</td>\n",
       "      <td>No</td>\n",
       "      <td>[{'name': 'Petition to list 225 species- cover...</td>\n",
       "      <td>2004</td>\n",
       "      <td>225</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>168</th>\n",
       "      <td>206 Mountain-Prairie Region species - List as T/E</td>\n",
       "      <td>07/30/2007</td>\n",
       "      <td>AR, AZ, CO, IA, ID, KS, MT, NE, NM, NV, OR, UT...</td>\n",
       "      <td>[Forest Guardians  Nicole J. Rosmarino  ]</td>\n",
       "      <td>[Listing: Threatened or Endangered, Listing: E...</td>\n",
       "      <td>[{'name': '90 Day Not Substantial on 01/06/200...</td>\n",
       "      <td>Yes</td>\n",
       "      <td>[{'name': 'Petition to list 206 Mountain-Prair...</td>\n",
       "      <td>2007</td>\n",
       "      <td>206</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>421</th>\n",
       "      <td>Northern goshawk (Accipiter gentilis) in the f...</td>\n",
       "      <td>09/26/1991</td>\n",
       "      <td>AK, AL, AR, AZ, CA, FL, GA, IA, ID, IL, IN, KS...</td>\n",
       "      <td>[Dr. Robin D. Silver  and  10 others ]</td>\n",
       "      <td>[Listing: Endangered, APA: Designate Critical ...</td>\n",
       "      <td>[{'name': '12 Month Not Warranted on 06/22/199...</td>\n",
       "      <td>No</td>\n",
       "      <td>[{'name': 'Northern goshawk petition', 'url': ...</td>\n",
       "      <td>1991</td>\n",
       "      <td>100</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                        Petition Title  \\\n",
       "171                 475 Southwestern species; List T/E   \n",
       "101           404 Southeast aquatic species; list w/CH   \n",
       "214   List 225 Candidate Species with Critical Habitat   \n",
       "168  206 Mountain-Prairie Region species - List as T/E   \n",
       "421  Northern goshawk (Accipiter gentilis) in the f...   \n",
       "\n",
       "    Date Received by the FWS  \\\n",
       "171               06/25/2007   \n",
       "101               04/20/2010   \n",
       "214               05/11/2004   \n",
       "168               07/30/2007   \n",
       "421               09/26/1991   \n",
       "\n",
       "    Where the species is believed to or known to occur  \\\n",
       "171  AL, AR, AZ, CA, CO, IL, KS, KY, LA, MT, NC, NM...   \n",
       "101  AL, AR, CO, CT, DC, DE, FL, GA, IA, IL, IN, KS...   \n",
       "214  AL, American Samoa, AR, AZ, CA, CO, DE, FL, GA...   \n",
       "168  AR, AZ, CO, IA, ID, KS, MT, NE, NM, NV, OR, UT...   \n",
       "421  AK, AL, AR, AZ, CA, FL, GA, IA, ID, IL, IN, KS...   \n",
       "\n",
       "                                       Petitioner Name  \\\n",
       "171  [Forest Guardians , Nicole J. Rosmarino, Conse...   \n",
       "101  [Alabama Rivers Alliance, Center for Biologica...   \n",
       "214  [Barbara Kingsolver, , Biodiversity Conservati...   \n",
       "168          [Forest Guardians  Nicole J. Rosmarino  ]   \n",
       "421             [Dr. Robin D. Silver  and  10 others ]   \n",
       "\n",
       "                                      Requested Action  \\\n",
       "171                [Listing: Threatened or Endangered]   \n",
       "101  [Listing: Threatened or Endangered, APA: Desig...   \n",
       "214  [Listing: Threatened or Endangered, Listing: T...   \n",
       "168  [Listing: Threatened or Endangered, Listing: E...   \n",
       "421  [Listing: Endangered, APA: Designate Critical ...   \n",
       "\n",
       "                                   Petition Finding(s) Active  \\\n",
       "171  [{'name': '12 Month Not Warranted on 09/27/200...    Yes   \n",
       "101  [{'name': '12 Month Warranted But Precluded on...    Yes   \n",
       "214  [{'name': '12 Month Not Warranted on 05/11/200...     No   \n",
       "168  [{'name': '90 Day Not Substantial on 01/06/200...    Yes   \n",
       "421  [{'name': '12 Month Not Warranted on 06/22/199...     No   \n",
       "\n",
       "                                 Petitions Document(s) year_received  \\\n",
       "171  [{'name': 'Petition to list 475 Southwest spec...          2007   \n",
       "101  [{'name': 'Withdrawal of Three crayfish', 'url...          2010   \n",
       "214  [{'name': 'Petition to list 225 species- cover...          2004   \n",
       "168  [{'name': 'Petition to list 206 Mountain-Prair...          2007   \n",
       "421  [{'name': 'Northern goshawk petition', 'url': ...          1991   \n",
       "\n",
       "     n_species  \n",
       "171        475  \n",
       "101        404  \n",
       "214        225  \n",
       "168        206  \n",
       "421        100  "
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_petitions['n_species'] = df_petitions.apply(lambda x: parseSpeciesNumber(x['Petition Title']), axis=1)\n",
    "df_petitions.at[576, 'n_species'] = 1\n",
    "df_petitions.at[577, 'n_species'] = 1\n",
    "df_petitions.at[582, 'n_species'] = 1\n",
    "df_petitions.at[598, 'n_species'] = 1\n",
    "df_petitions.at[715, 'n_species'] = 1\n",
    "df_petitions.sort_values(by='n_species', ascending=False).head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's look at a bar graph to see the number of species petitioned per year."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6gAAAH0CAYAAAAwkgCNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3X24rHVdL/73hwcBRUVgi8pGocJEiwdBxDRT8QG1hDpieFTQPGHqSbP0pFlHsietflmc05GDYcLJNNEMTnoUw8i0QJEnQTRQCbYo7FAwREzx+/vjvjcsF2vvvZhZw/7O2q/Xdc21Zu6513u+c6+ZWfOe+2GqtRYAAADY0rbZ0gMAAACAREEFAACgEwoqAAAAXVBQAQAA6IKCCgAAQBcUVAAAALqgoAIAANAFBRUAAIAuKKgAAAB0QUEFAACgC9tt6QEkye6779723nvvLT0MAAAAZuDTn/70v7XW1mxuvi4K6t57753zzz9/Sw8DAACAGaiqf13OfDbxBQAAoAsKKgAAAF1QUAEAAOhCF/ugAgAAbAnf+c53sm7dutx6661beiirwo477pi1a9dm++23n+j3FVQAAGCrtW7dutz73vfO3nvvnara0sOZa6213HDDDVm3bl322WefiTJs4gsAAGy1br311uy2227K6Qqoquy2225TrY1WUAEAgK2acrpypl2WCioAAABdsA8qAADAaKVXpra2snmzdNJJJ+We97xnjj322C02BgUVAACA/MIv/MKWHoJNfAEAALakq666Kvvtt19+/ud/Po94xCPy1Kc+Nd/61reWnPfEE0/Mwx/+8Oy///455phjkiQnnHBCXvCCF+RJT3pS9t1337ztbW+7ff4/+IM/yKMe9ajsv//+ecMb3nD79NNOOy37779/DjjggLzgBS+4PecP//APkyRf+MIXcsQRR+Tggw/Oj//4j+dzn/tckuT000/Pj/zIj+SAAw7I4x//+BVfFtagAgAAbGFXXHFF3vWud+Vtb3tbnvOc5+R973tfnv/8599pvje96U350pe+lB122CE33njj7dMvueSSnHvuufnmN7+Zgw46KM985jNz6aWX5oorrsgnP/nJtNbyrGc9Kx/72Mey22675Xd+53fyiU98Irvvvnu+9rWv3el2jj/++Jx00knZd999c9555+VlL3tZPvrRj+aNb3xjPvzhD2fPPff8vttfKQoqAADAFrbPPvvkwAMPTJIcfPDBueqqq5acb//998/znve8HHXUUTnqqKNun37kkUdmp512yk477ZQnPvGJ+eQnP5mPf/zjOeuss3LQQQclSW6++eZcccUVufjii/PsZz87u+++e5Jk1113/b7buPnmm/NP//RPOfroo2+f9u1vfztJ8tjHPjYvfOEL85znPCc/8zM/s2L3fwMFFQAAYAvbYYcdbj+/7bbbbnQT3w984AP52Mc+ljPPPDO/9Vu/lcsuuyzJnb/eparSWsvrXve6vOQlL/m+60488cRNfh3M9773veyyyy656KKL7nTdSSedlPPOOy8f+MAHcuCBB+aiiy7Kbrvttuz7uTn2QQUAAJgD3/ve93LNNdfkiU98Yn7/938/N954Y26++eYkyRlnnJFbb701N9xwQ84555w86lGPytOe9rS8/e1vv32eL3/5y7n++utz+OGH5z3veU9uuOGGJLnTJr73uc99ss8+++T0009PkrTWcvHFFycZ9k199KMfnTe+8Y3Zfffdc80116zofbQGFQAAYNTz18Lcdtttef7zn5+bbroprbW86lWvyi677JIkOfTQQ/PMZz4zV199dX7jN34jD3rQg/KgBz0ol19+eR7zmMckSXbeeef8xV/8RR7xiEfk9a9/fX7iJ34i2267bQ466KC84x3v+L7beuc735mXvvSl+e3f/u185zvfyTHHHJMDDjggr3nNa3LFFVektZbDDz88BxxwwIrex2od/AUOOeSQdv7552/pYQAAAFuZyy+/PPvtt9+WHsZUTjjhhOy888559atfvaWHkmTpZVpVn26tHbK5393sJr5V9cNVddGC0zeq6peqateq+khVXTH+vN84f1XViVV1ZVVdUlWPnPieAQAAsNXY7Ca+rbXPJzkwSapq2yRfTvL+JK9NcnZr7U1V9drx8q8meXqSfcfTo5O8dfwJAMBWYBPHXrmTDjbmgy69/OUvzyc+8Ynvm/bKV74yL3rRi+407wknnHA3jWr27uo+qIcn+UJr7V+r6sgkTxinn5rknAwF9cgkp7Vh2+Fzq2qXqnpga+0rKzRmAACAFdNa2+RRbbeEP/3TP93SQ5jItLuQ3tWj+B6T5F3j+T02lM7x5/3H6XsmWXgop3XjNAAAgK7suOOOueGGG6YuVgzl9IYbbsiOO+44ccay16BW1T2SPCvJ6zY36xLT7vTXrqrjkxyfJA9+8IOXOwwAAIAVs3bt2qxbty7r16/f0kNZFXbcccesXbt24t+/K5v4Pj3JBa2168bL123YdLeqHpjk+nH6uiR7Lfi9tUmuXRzWWjs5ycnJcBTfuzxyAACAKW2//fbZZ599tvQwGN2VTXyfmzs2702SM5McN54/LskZC6YfOx7N97AkN9n/FAAAgM1Z1hrUqrpnkqckecmCyW9K8p6qenGSq5McPU7/YJJnJLkyyS1J7nyYKQAAAFhkWQW1tXZLkt0WTbshw1F9F8/bkrx8RUYHAADAVuOuHsUXAAAAZkJBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAvLKqhVtUtVvbeqPldVl1fVY6pq16r6SFVdMf683zhvVdWJVXVlVV1SVY+c7V0AAABgNVjuGtQ/SfKh1trDkhyQ5PIkr01ydmtt3yRnj5eT5OlJ9h1Pxyd564qOGAAAgFVpswW1qu6T5PFJTkmS1tp/tNZuTHJkklPH2U5NctR4/sgkp7XBuUl2qaoHrvjIAQAAWFWWswb1B5KsT/LnVXVhVf1ZVd0ryR6tta8kyfjz/uP8eya5ZsHvrxunfZ+qOr6qzq+q89evXz/VnQAAAGD+LaegbpfkkUne2lo7KMk3c8fmvEupJaa1O01o7eTW2iGttUPWrFmzrMECAACwei2noK5Lsq61dt54+b0ZCut1GzbdHX9ev2D+vRb8/tok167McAEAAFitNltQW2tfTXJNVf3wOOnwJJ9NcmaS48ZpxyU5Yzx/ZpJjx6P5Hpbkpg2bAgMAAMDGbLfM+X4xyTur6h5JvpjkRRnK7Xuq6sVJrk5y9DjvB5M8I8mVSW4Z5wUAAIBNWlZBba1dlOSQJa46fIl5W5KXTzkuAAAAtjLL/R5UAAAAmCkFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0IVlFdSquqqqPlNVF1XV+eO0XavqI1V1xfjzfuP0qqoTq+rKqrqkqh45yzsAAADA6nBX1qA+sbV2YGvtkPHya5Oc3VrbN8nZ4+UkeXqSfcfT8UneulKDBQAAYPWaZhPfI5OcOp4/NclRC6af1gbnJtmlqh44xe0AAACwFVhuQW1JzqqqT1fV8eO0PVprX0mS8ef9x+l7Jrlmwe+uG6cBAADARm23zPke21q7tqrun+QjVfW5TcxbS0xrd5ppKLrHJ8mDH/zgZQ4DAACA1WpZa1Bba9eOP69P8v4khya5bsOmu+PP68fZ1yXZa8Gvr01y7RKZJ7fWDmmtHbJmzZrJ7wEAAACrwmYLalXdq6ruveF8kqcmuTTJmUmOG2c7LskZ4/kzkxw7Hs33sCQ3bdgUGAAAADZmOZv47pHk/VW1Yf6/bK19qKo+leQ9VfXiJFcnOXqc/4NJnpHkyiS3JHnRio8aAADmQC2189tGtDvtFAdbn80W1NbaF5McsMT0G5IcvsT0luTlKzI6AAAAthrTfM0MAAAArBgFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4su6BW1bZVdWFV/e14eZ+qOq+qrqiqv6qqe4zTdxgvXzlev/dshg4AAMBqclfWoL4yyeULLr85yVtaa/sm+XqSF4/TX5zk6621H0rylnE+AAAA2KRlFdSqWpvkmUn+bLxcSZ6U5L3jLKcmOWo8f+R4OeP1h4/zAwAAwEYtdw3qHyf5b0m+N17eLcmNrbXvjpfXJdlzPL9nkmuSZLz+pnF+AAAA2KjNFtSq+skk17fWPr1w8hKztmVctzD3+Ko6v6rOX79+/bIGCwAAwOq1nDWoj03yrKq6Ksm7M2za+8dJdqmq7cZ51ia5djy/LsleSTJef98kX1sc2lo7ubV2SGvtkDVr1kx1JwAAAJh/my2orbXXtdbWttb2TnJMko+21p6X5O+TPHuc7bgkZ4znzxwvZ7z+o621O61BBQAAgIWm+R7UX03yy1V1ZYZ9TE8Zp5+SZLdx+i8nee10QwQAAGBrsN3mZ7lDa+2cJOeM57+Y5NAl5rk1ydErMDYAAAC2ItOsQQUAAIAVo6ACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQhc0W1Krasao+WVUXV9VlVfWb4/R9quq8qrqiqv6qqu4xTt9hvHzleP3es70LAAAArAbLWYP67SRPaq0dkOTAJEdU1WFJ3pzkLa21fZN8PcmLx/lfnOTrrbUfSvKWcT4AAADYpM0W1Da4eby4/XhqSZ6U5L3j9FOTHDWeP3K8nPH6w6uqVmzEAAAArErL2ge1qratqouSXJ/kI0m+kOTG1tp3x1nWJdlzPL9nkmuSZLz+piS7reSgAQAAWH2WVVBba7e11g5MsjbJoUn2W2q28edSa0vb4glVdXxVnV9V569fv3654wUAAGCVuktH8W2t3ZjknCSHJdmlqrYbr1qb5Nrx/LokeyXJeP19k3xtiayTW2uHtNYOWbNmzWSjBwAAYNVYzlF811TVLuP5nZI8OcnlSf4+ybPH2Y5LcsZ4/szxcsbrP9pau9MaVAAAAFhou83PkgcmObWqts1QaN/TWvvbqvpskndX1W8nuTDJKeP8pyT5P1V1ZYY1p8fMYNwAAACsMpstqK21S5IctMT0L2bYH3Xx9FuTHL0iowMAAGCrcZf2QQUAAIBZUVABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKAL223pAQAArAZVy5+3tdmNA2CeWYMKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFxRUAAAAuqCgAgAA0AUFFQAAgC4oqAAAAHRBQQUAAKALCioAAABdUFABAADogoIKAABAFzZbUKtqr6r6+6q6vKouq6pXjtN3raqPVNUV48/7jdOrqk6sqiur6pKqeuSs7wQAAADzbzlrUL+b5Fdaa/slOSzJy6vq4Ulem+Ts1tq+Sc4eLyfJ05PsO56OT/LWFR81AAAAq85mC2pr7SuttQvG8/+e5PIkeyY5Msmp42ynJjlqPH9kktPa4Nwku1TVA1d85AAAAKwqd2kf1KraO8lBSc5Lskdr7SvJUGKT3H+cbc8k1yz4tXXjNAAAANioZRfUqto5yfuS/FJr7RubmnWJaW2JvOOr6vyqOn/9+vXLHQYAAACr1LIKalVtn6GcvrO19tfj5Os2bLo7/rx+nL4uyV4Lfn1tkmsXZ7bWTm6tHdJaO2TNmjWTjh8AAIBVYjlH8a0kpyS5vLX2RwuuOjPJceP545KcsWD6sePRfA9LctOGTYEBAABgY7ZbxjyPTfKCJJ+pqovGab+W5E1J3lNVL05ydZKjx+s+mOQZSa5MckuSF63oiAEAAFiVNltQW2sfz9L7lSbJ4UvM35K8fMpxAQAAsJW5S0fxBQAAgFlRUAEAAOiCggoAAEAXFFQAAAC6oKACAADQheV8zQwAc6I2dsz1jWhtNuMAAJiENagAAAB0QUEFAACgCwoqAAAAXVBQAQAA6IKCCgAAQBcUVAAAALqgoAIAANAFBRUAAIAuKKgAAAB0QUEFAACgCwoqAAAAXVBQAQAA6MJ2W3oAANOqumvztzabcQAAMB1rUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXFFQAAAC6oKACAADQBQUVAACALiioAAAAdEFBBQAAoAubLahV9faqur6qLl0wbdeq+khVXTH+vN84varqxKq6sqouqapHznLwAAAArB7LWYP6jiRHLJr22iRnt9b2TXL2eDlJnp5k3/F0fJK3rswwAQAAWO02W1Bbax9L8rVFk49Mcup4/tQkRy2YflobnJtkl6p64EoNFgAAgNVr0n1Q92itfSVJxp/3H6fvmeSaBfOtG6cBAADAJq30QZJqiWltyRmrjq+q86vq/PXr16/wMAAAAJg3kxbU6zZsujv+vH6cvi7JXgvmW5vk2qUCWmsnt9YOaa0dsmbNmgmHAQAAwGoxaUE9M8lx4/njkpyxYPqx49F8D0ty04ZNgQEAAGBTttvcDFX1riRPSLJ7Va1L8oYkb0rynqp6cZKrkxw9zv7BJM9IcmWSW5K8aAZjBgAAYBXabEFtrT13I1cdvsS8LcnLpx0UAAAAW5+VPkgSAAAATERBBQAAoAsKKgAAAF1QUAEAAOiCggoAAEAXNnsUXwAAWM2qlj9va7MbB2ANKgAAAJ1QUAEAAOiCggoAAEAX7IMKAAAksT8uW541qAAAAHRBQQUAAKALNvEFAGAu2PwUVj9rUAEAAOiCggoAAEAXFFQAAAC6YB9UAACAu8ld2Zc62fr2p7YGFQAAgC4oqAAAAHRBQQUAAKAL9kEF2ATfuQcAcPexBhUAAIAuWIMKsAVYMwsAcGfWoAIAANAFBRUAAIAuKKgAAAB0QUEFAACgCwoqAAAAXVBQAQAA6IKCCgAAQBcUVAAAALqgoAIAANAFBRUAAIAubLelBwBsPaqWP29rsxsHAAB9UlCZCUUEgGn4PwKwdVJQAWCVuytlL1H4eqOsw6Z5jqwuCioAAMAiiu+W4SBJAAAAdMEaVABgYtYwzC9/O6BHCioAy+LN7OxZxgBs7RRUAFYlBwZiKT4EAOibfVABAADogjWoAAAwZ2wNwGqloALAXeSNIQDMhk18AQAA6IKCCgAAQBds4svcsWkdAACsTgoqjBRf2DI89wBgZayG/6ldF9TVsIABAABYnq4LKgAAsDpY+cRyzOQgSVV1RFV9vqqurKrXzuI2plG1/NNqZ1kAAACbcnd2hhVfg1pV2yb50yRPSbIuyaeq6szW2mdX+rZ641Mh7m4ecwAArCazWIN6aJIrW2tfbK39R5J3JzlyBrezVbGmk7uTxxsAAFvCLPZB3TPJNQsur0vy6MUzVdXxSY4fL95cVZ9fZv7uSf7tznl3cZTLzJ633Flmz1vuCmXfrct4hayKZeFxMfvcFcq2LGafu2S2ZTG/ubPMnrfcWWbPW+4KZa+KZbxC2ZbF7HOXzN5E7kOWEziLgrrUkO60cWFr7eQkJ9/l8KrzW2uHTDKwLZU9b7mzzJY7++x5y51l9rzlzjJ73nJnmT1vubPMnrfcWWbPW+4ss+ctd5bZcmefPW+5s8yet9xZZc9iE991SfZacHltkmtncDsAAACsIrMoqJ9Ksm9V7VNV90hyTJIzZ3A7AAAArCIrvolva+27VfVfk3w4ybZJ3t5au2wFb+IubxbcQfa85c4yW+7ss+ctd5bZ85Y7y+x5y51l9rzlzjJ73nJnmT1vubPMnrfcWWbLnX32vOXOMnvecmeSXc13TwAAANCBWWziCwAAAHeZggoAAEAXFFQAAAC6oKACAADQBQUVALYSVXXfqvrZqvrlqnrVeH6XGd7eU6b8/ftU1Q8uMX3/aXLHjAdU1QPG82uq6meq6hHT5i5xO787g8x9xvE+bAWyHlxVO47nq6peVFX/o6peWlUTf9tDVT1rQ+5Kq6rHV9UPj+cfV1WvrqpnrlD2zlX17PH58YtVdURVTf1+uaoeVlW/WlUnVtWfjOf3W4kxb+T2XrQCGQ+rqsOraudF04+YMvfQqnrUeP7h4+vRM6bJ3MjtnLbSmWPu48YxP3XKnEdX1X3G8ztV1W9W1f+tqjdX1X2nyH1FVe01zdg2knuPqjq2qp48Xv7PVfU/q+rlVbX9it7WvB3Ft6o+2lp70gplPTHJf0qyV5LvJrkiyZ+11q6cMvdhSY5MsmeSluTaJGe21i6fbsQbvb0Xtdb+fMqMh2UY73mttZsXTD+itfahKXIPTdJaa5+qqocnOSLJ51prH5xmvEvczmmttWNXMnPMfVySQ5Nc2lo7a4qcRye5vLX2jaraKclrkzwyyWeT/G5r7aYJczd81/C1rbW/q6r/nOTHklye5OTW2ncmHfOY/4NJfjrf/xx516TjXZD7tCRH5fufI2dM81jbzO3999baG6fMeFqStUnObq1dtWD6z7XW3j5hZiU5OsMyeG+SJ2V47fhckpNaa9+bZsyLbmvq186q2r219m8LLj8/4/MjydvaFP9Qquqnk/xDa+1rVbUmyf+X5KAMz5Ffaa2tmzD3j5K8r7X2iUnHtpHcXZP81wyP3VOS/FqSx2R47v1ua+3rU+av+P+nqjo2yRuSnJXky+PktUmekuQ3W2sr/mauqq5urT14wt99TpI/TnJ9ku2TvLC19qnxugtaa4+cYlwvyfA6XEnenOSFSS5L8tgkv99aO2XC3BMXT0rygiSnJUlr7RUT5v5Na+2o8fyRGZbLORle73+vtfaOSXLHvEuTHNpau6Wq3pzkB5P8TYbXo7TWfm7C3G8l+WaS/5fkXUk+3Fq7bdJxLsj94wyvO9tl+ErDw8fb+IkkF7bWXjNF9nOSvCbJxUmemOSfMqzM+dEkz2utfWbC3F9N8twk706y4bVsbYb/3+9urb1p0jFv4jYnfu6Nv/+KJC/P8Jp2YJJXttbOGK+b+PlXVW9I8vQMf7+PJHl0hsfykzM8Rn5nwtwzF0/K8Df8aJK01p41Se6Y/cnW2qHj+Z/PsFzen+SpSf7vpH+/qrosyQHjV3SenOSWDO8FDh+n/8yEuTdleO59IcNz7/TW2vpJshblvjPD3+2eSW5MsnOSvx7HW62146a9jdtvq+eCWlWXLJ6U5KFJPp8krbWJP0Gtqjcl2SPJ2RneKH8pyb8keVmGNxenT5jrReiOXC9Cd+TO6kVoZi8W4+Pip5L8Q5JnJLkoydczFNaXtdbOmTD3jzM8j0/L9z9Hjk1yRWvtlZOOeRPXzKzKAAAUH0lEQVS3Oe1z5HeTPC7JBRmWyR+31v7HeN00z5H/leT+Se6R5BtJdkjyfzMs7+smXRazeu1ceF+r6teT/HiSv0zyk0nWtdZeNUnumPfZ1trDx/N/leTcJKdneL14XmttojVxVbU+yb8mWZPkrzJ8wHLhpONckPvBJJ9Jcp8k+43n35Oh7B3QWjtyiuxZ/X/6fJJHt9ZuXDT9fhk+nHzohLmLX5NvvyrJk1pr95ow96IkT2+tfWX8sPO0JL/WWvvrqrqwtXbQJLlj9mcy/E/aKcPj44daa18dl8Xft9YOnDB3XYb/c2dluP9J8odJXp0krbVTJ8y9/f5W1T9leE58qap2z/Ch2QGT5I55C597n07yqA0fjlXVxZNmV9WFGUruszO8B/qRDP9P39Va+4cpxnvZmLVThg9a9hzL9fYZCuqPTJF9SZLDxrzdk7yztfa0GtbYn9Ra+7EJc/8lySMWf2g8fsh8WWtt3ynGu+RVSR7aWtthktwx+zNJHtNau7mq9s7wnuX/tNb+ZJrn35h7YIb/d19NsnbBh/fnTfP/KcMHmn+W4UPfylDOjkmSKR9zC59/n0ryjNba+qq6V5JzW2s/OmHu5a21/TaMf+F7iaq6aIrXoQuTHJzh/+fPJnlWkk9nWB5/3Vr79wlzL2mt7V/DlhVfTvKg1tpt44ftF0/Ty+6ktdbtKcmZSf4iycOSPCTJ3kmuGc8/ZMrszyw4v12ST4zn75dhbdmkuf+SZPslpt8jw5vvSXMv2cjpM0m+Pe2ySLLzeH7vJOdnKKnJ8GI/Te62GcrTN5LcZ5y+U5JLpsi9YHxcPCHDJ6ZPSPKV8fxPTLksLlxw/lNJ1ozn77XwMTNB7uULx7/ououmeVyMP7dLcl2SbcfLNc0yXvj3G8/fM8k54/kHT/m4+JeNTK8pnyPf2Mjp35N8dwWWxXbj+V2SfDDJWxY/ZibJHX9un+SGJPdY8Pec5vE2k9fORc+PC5Lca8H4Jx7vmPH5Bec/vei6aZ4jF44/903yGxnWkn0uw5rEh06Re1G743H75ZUa78LHxYLHwkr+f7rvEtPvO+Vz7+tJnrnhNXjB6QkZPmiZejmMlx+Y4U3WKxa/jk6QfcGC8xcv9ZiZMPfeGdZu/mWG4pQkX5xmrEuM95MrNd7x9z+c4YOEJHnfhteIJLstXjaTjnm8/IDxb/fPSa6ZIvfS8eeO42Nvp/Hytkk+O+Wy+EzuWHmz06LXvGmee59b6rV3fE3+/BS512Uoew9ZdNo7w5ZV0yyLzy66vHOSDyX5o5V4TV58frw8Te42SV6VYYXIgeO0qZ97Y87F4+vvbknO39j9mSD39CQvGs//eZJDxvMPTfKpKXIXP/e2z1BS35Vk/RS5l2boM/fL8N5q13H6jlnwPnclThPvW3B3aK09a9zs6+Qkf9haO7OqvtNa+9cViP9eVe3aWvtakgdleGFLa+3r4ycBE+eOeYvH+MDxukntkeRpGV6MF6oMm6BMY9s2btbbWruqqp6Q5L1V9ZDc8QnwJL7bhs15bqmqL7TWvjHexreqapplcUiSVyZ5fZLXtNYuqqpvtSk+HVtgm/ET9G0y/JNanySttW9W1XenyL207tgU++KqOqS1dn5VPTTJNJvhbjN+AnuvDCXyvkm+luGTyZXYH2C7JLeNefdOktba1TXdvga3VtWhrbVPLpr+qCS3TpF7Y4ZP/q9bfEVVXTNFbjKU0+8mSWvtxqr6qSQnV9XpGV6sJ7Uh8ztV9anW2n+Ml79bVRNvCjfD186dquqgDM+PbVtr31ww/mk33Tunqt6Y5PfG80e11v5m3NR1mk3K2zjGK5L8VpLfGteGPDfDBw0/NGHuhteKeyfZuar2Hl8/d8t0j4lkdv+ffifJBVV1VoYPLJLhA6enZFg2kzo3yS1LvQaPa20n9e9V9YOttS8kSRvWpD4hw+an0+4r+r2q2r4Na7Ru33exhn0mJ97fsA1rJn6pqg5O8hdV9YFp8hY4oKq+keF/8g5V9YA2rPG9R8bHxxT+S5LTquqEDM+1i8Y1MPdL8stT5H7fY7W19tUkJyY5cXx/MakPVNU/ZnhT/GdJ3lNV52b4UORjU+Qmw2vCh6rqHzJsAXZ6cvsm/dM8934pydlVdUW+/7n3Qxl2FZjU32ZYwXDR4iuq6pwpcpPkq1V14IbsNqxJ/ckkb8+wyfOk/qOq7tlauyXDWr4kSQ37XE783rANa/3fMv5ffktVXZesWM+5b4YPxypJW/D82znTPS7+S5I/GbdI+rck/zy+X7lmvG5Si59738nwwfWZ45rqSZ2S4cOWbTO8Bz+9qr6Y5LAMW46unJVsu7M6ZXjz/UcZFu66Fcr82Qwl8qwkVyd55jh9TZK/nCL3iCRXZtgf4uTx9KFx2hFT5J6S5HEbuW7i8Y6//9GMnzYtmLZdhs2pbpsi97wk9xzPb7Ng+n0z5affY87aDP88/meSq1focXFVki9m2KTui0keME7fOdN9snffJO/IsD/AeRlK6RczbD57wBS5rxpz/jXDJ9NnJ3lbhk+B3zDlsnhlhrX0J2d4QdrwKd+aJB+bIveR4zL47Pj8OyvD5uXnJTl4itzfzrAf1VLXvXnKZfG3WWLt/Hib35si9/9l3Hph0fQHZNFakgnzV/S1M8nfLzo9cJx+p0+VJ8jePskJ4+vx1RneqPx7hjVRD54id6q1S5vIfW6GtRfXZdhX9O8yfHL/5STHT5k9k/9PY8b9Mmzy9isZNjs9Jsn9ZrGMVmAZH5Bk3408Vp43ZfaDs/TWTnsmefIKjb8y7CbyFzNcRrtk2AxzJbL2y7AP/H/KsPnzNlPmPWGG9/sxGTbFTYZ9Zl+d5DnTjnnMe8aY95QF07ZJssOUudtkeCP/nzJs9nxYxq2UejxleI/1gI1c99gpcpdcjkl2T/KjKzj+Z2bYJWKWy+ieSfZZgZx7j693ByfZYwXyJt46aBnZD8qwae+G159nZyPvvaY5db0P6mJVdUCGF+KTVihv1yQ/kOTKtmifnClzt8mwA/+eGf5Brcuwqn7qgwPMQlWtzbC286tLXPfYNuHBRapqh9bat5eYvnuGN7YTHWxgibxnZnix/LWVyNvIbdwzw4vGl6bMuXeGx9x2GQrDndb2TZD5oCRprV1bw9E4n5yhsC9eQzlJ9iMyvGm5tLX2uWnzFmU/IAueI0s9/nqx4RPH1tq3lrhuz9bal+/8W1Pd3r0ybD57/Qrlrehr5xL522Z403HLCuXdN8Na6xtWIGvntuDAbytpvN/VhjXe22XY1O7LrbWvrED2TP4/jdl7ZMEBylbidWgec2eZPW+5s8yet9xZZy9xWzN5jZrxa99cjdmymL/crgtqVe3fWtvYDuA9Zz84yTfasCng3hk2Sb28tXbZDHI/11q7dLoRzy573nJnmT1vuXdD9iFZcITSlSrAs8qdZbbc2WfPW+4ssqvqwCQnZdiiY12GD4fWZtg8/mWttQsmzD0oyVvH3IVHB542d+F4F+e+tE1xsKtZZW8md1bLeJbLYlZ/v+5yZ529iduc6iB+d3fuLLPnLXeW2Vt7btf7oCa5sKq+lGGn3ne11j7be3ZVvTbJS5J8u6o2HLnvE0l+s6pOaa39UU+58zhmy2L2uTMe809k+AqRGzNszvKJJPerqu8keUFrbaL9RWeVO49jnrfceRzzPC6LDLsZvKS1dt6i2zsswwE6Jj0S7J/PKPcdm8h9xxS5s8zeVO6slvE7pshNZjfmecudWXZVbWxf3sqwC9FEZpU7y+x5y51lttxN3Fbna1AvzPD9Yc/NsE/ONzMUyne3Bd9D2FN2DYc/PyTDdulXJfmBdsehqM9rEx7+fFa58zhmy2L2uTMe84VJnjpm7ZPkj1prP11VT8lw0KuJvvh6VrnzOOZ5y53HMc/psriibeTrLKrqytbaRAeMmrfcWWbPW+4ss+ctd5bZVXVrkj/IeGC8RV7VWtulp9xZZs9b7iyz5W5c72tQWxs2JXx9ktfX8D1oxyT5x6q6pk34fVQzzr6tDUep/Y8k38rw1RFpw1FgpxjuzHLnccyWxexzZ5m9bbvjC6OvznBI/LTWPlLDd6T2ljvLbLmzz5633Flm/78ajip7Wu44kuheGb6D+ENbUe4ss+ctd5bZ85Y7y+wLkvxNa+3Ti6+oqmmO1jqr3Flmz1vuLLPlbkT3a1DbEl8EXMO748e3FfrS3ZXMrqp3ZPiKgXsluSXDpwwfyvBl1fdurT2np9x5HLNlMfvcGY/57RkOPHF2hiNGfrm19ss1HIjqgtbaw3rKnccxz1vuPI55HpfFmP30MXPhQfzObK19cNLMecydZfa85c4ye95yZ5VdVT+c5GsLPnhaeN0ebcKDMM0qd5bZ85Y7y2y5m7itzgvqf26t/eU8ZddwJMejM7y5eG+GQ7U/N8On4H/axu8N7CV3HsdsWcw+d8Zj3j7Jzyd5eIYvv357a+22Go6Ue/824Xd1zip3Hsc8b7nzOOZ5XBYAMA+6LqgAwMqo4et7Xpdh7dD9x8nXJzkjyZvahF9nM2+58zhmy2L2uXfTmI/K8H3GXefO45gti/nNXco2KxU0C1W1c1W9saouq6qbqmp9VZ1bVS/sNXsTucf1mDuPY7YsZp97N4350hk991Y0dx7HPG+58zjmeVwWSd6T5OtJntha2621tluSJ2Y4WvDpW1HuPI7Zsph97t0x5icsyv16p7nzOGbLYn5z76TrNahVdUaS9yf5uyTPybAv3LuT/HqGfXJ+rbfsecudxzFbFrPPnccxWxbzmzuPY57TZfH51toP39XrVlvuLLPnLXeW2fOWO8vsecudZfa85c4yW+4mtNa6PSW5eNHlT40/t0nyuR6z5y13HsdsWVgWlsXqyp3HMc/psjgryX9LsseCaXsk+dUkf7e15M7jmC0Ly8KysCxWa+5Sp6438U3yzap6XJJU1U8l+VqStNa+l2Ta79CYVfa85c4ye95yZ5k9b7mzzJ633Flmy5199rzlzjL7Z5PsluQfqurrVfW1JOck2TXDmtqtJXeW2fOWO8vsecudZfa85c4ye95yZ5ktd2NWsu2u9CnJ/kk+mWHb/48neeg4fU2SV/SYPW+58zhmy8KysCxWV+48jnkel8WY8bAkT06y86LpR2xNufM4ZsvCsrAsLIvVmnun21nJsLvzlORF85Y9b7nzOGbLwrKwLFZX7jyOuddlkeQVST6f5G+SXJXkyAXXXbC15M7jmC0Ly8KysCxWa+6St7WSYXfnKcnV85Y9b7nzOGbLwrKwLFZX7jyOuddlkeQzGT/1TrJ3kvOTvHK8fOHWkjuPY7YsLAvLwrJYrblLnbZLx6rqko1dlWGn3O6y5y13ltnzljvL7HnLnWX2vOXOMlvu7LPnLXfG2du21m5OktbaVVX1hCTvraqHZLp9W+ctdx7HbFnMPncex2xZzD53Hsc8b7l30nVBzfCP+GkZvl9noUryT51mz1vuLLPnLXeW2fOWO8vsecudZbbc2WfPW+4ss79aVQe21i5KktbazVX1k0nenuRHt6LceRyzZTH73Hkcs2Ux+9x5HPO85d7ZSq6OXelTklOSPG4j1/1lj9nzljuPY7YsLAvLYnXlzuOY53RZrE3ygI1c99itJXcex2xZWBaWhWWxWnOXOtUYCgAAAFtU79+DCgAAwFZCQQUAAKALCioAAABdUFABYI5U1QerapcVyDmhql69EmMCgJXS+9fMAEDXqmrb1tpt086zXK21Z6xEDgD0yBpUALYaVfVbVfXKBZd/p6peUVWvqapPVdUlVfWbC67/m6r6dFVdVlXHL5h+c1W9sarOS/KYjdzWVVX136vq40mOrqofrKoPjXn/WFUPG+fbo6reX1UXj6cfG6c/v6o+WVUXVdX/rqptF+TuXlVvrqqXLbi9E6rqV8bzG7s/r6+qz1fV3yX54RVarACwYhRUALYmpyQ5LkmqapskxyS5Lsm+SQ5NcmCSg6vq8eP8P9daOzjJIUleUVW7jdPvleTS1tqjW2sf38Tt3dpae1xr7d1JTk7yi2Peq5P8r3GeE5P8Q2vtgCSPTHJZVe2X5GczfLfcgUluS/K8RdnvHufZ4DlJTq+qpy51f6rq4PH+HpTkZ5I8ahnLCwDuVjbxBWCr0Vq7qqpuqKqDkuyR5MIMRe2p4/kk2TlDwftYhlL60+P0vcbpN2QojO9bxk3+VZJU1c5JfixDgdxw3Q7jzyclOXYc321JbqqqFyQ5OMmnxvl3SnL9ovtyYVXdv6oelGRNkq+31q6uqlds5P7cO8n7W2u3jGM6cxnjB4C7lYIKwNbmz5K8MMkDkrw9yeFJfq+19r8XzlRVT0jy5CSPaa3dUlXnJNlxvPrWZe5T+s3x5zZJbhzXhi5HJTm1tfa6zcz33iTPznBf3r3gd5e6P7+UpC3z9gFgi7CJLwBbm/cnOSLDmtMPj6efG9dypqr2rKr7J7lvhrWSt4z7ix426Q221r6R5EtVdfR4G1VVB4xXn53kpeP0bavqPuO0Z4/jSFXtWlUPWSL63Rk22312hrKaTdyfjyX56araqaruneSnJr0/ADAr1qACsFVprf1HVf19hjWatyU5a9zn85/HzWlvTvL8JB9K8gtVdUmSzyc5d8qbfl6St1bVryfZPkO5vDjJK5OcXFUvzrDp8Etba/88znfWuK/sd5K8PMm/Lrovl41l88utta+M05a8P621C6rqr5JcNOb845T3BwBWXLVmax8Ath5j4bsgydGttSu29HgAgDvYxBeArUZVPTzJlUnOVk4BoD/WoALAFKrq/Un2WTT5V1trH94S4wGAeaagAgAA0AWb+AIAANAFBRUAAIAuKKgAAAB0QUEFAACgCwoqAAAAXfj/ATnO8LDfjq6JAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 1152x576 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "petitions_per_year = df_petitions[['year_received', 'n_species']].groupby('year_received').sum()\n",
    "plt.rcParams['figure.figsize'] = [16, 8]\n",
    "petitions_per_year.plot(kind='bar', color='blue');"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Combining species and petition data\n",
    "\n",
    "We can merge the two DataFrames that summarize species and petitions by year."
   ]
  },
  {
   "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>Species Petitioned</th>\n",
       "      <th>Species Added</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>1967</th>\n",
       "      <td>0.0</td>\n",
       "      <td>81.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1970</th>\n",
       "      <td>0.0</td>\n",
       "      <td>60.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1972</th>\n",
       "      <td>0.0</td>\n",
       "      <td>3.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1973</th>\n",
       "      <td>0.0</td>\n",
       "      <td>8.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1974</th>\n",
       "      <td>9.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "      Species Petitioned  Species Added\n",
       "1967                 0.0           81.0\n",
       "1970                 0.0           60.0\n",
       "1972                 0.0            3.0\n",
       "1973                 0.0            8.0\n",
       "1974                 9.0            0.0"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df_merged = pd.concat([petitions_per_year, year_counts], axis=1, sort=True).fillna(0)\n",
    "df_merged = df_merged.rename(index=str, columns={\"listyear\": \"Species Added\", \"n_species\": \"Species Petitioned\"})\n",
    "df_merged.head()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can plot the merged DataFrame to see how the species added and the species petitioned compare each year."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6gAAAHlCAYAAAD4GNGzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzs3Xu8XGV9P/rPQxLuV028ESW0QoVCiBgCVMRUaqqgQAWL1CqkWrQiVfy1B37Vc7Q/UOk5HmmpHBEvBX8ComALVYoKCliKSBLDRYKCNEgAIeEqAirkOX+slbiT7JCZZJKsTN7v12tee80za333M3vWXjOftZ5Zq9RaAwAAAOvbJuu7AwAAAJAIqAAAAHSEgAoAAEAnCKgAAAB0goAKAABAJwioAAAAdIKACgAAQCcIqAAAAHSCgAoAAEAnjF3fHUiS8ePH10mTJq3vbgAAALAWzJ49e1GtdcKq5utEQJ00aVJmzZq1vrsBAADAWlBKuauX+QzxBQAAoBMEVAAAADpBQAUAAKATOvEdVAAAYHj85je/yYIFC/LUU0+t766wjm2++eaZOHFixo0bt1rLC6gAAMBALViwINtss00mTZqUUsr67g7rSK01Dz74YBYsWJCdd955tWoY4gsAAAzUU089lec+97nC6UamlJLnPve5a3TkXEAFAAAGTjjdOK3p6y6gAgAA0Am+gwoAAKxVk07+xkDrzT/tkFXO89GPfjTnn39+xowZk0022SSf+cxnsu+++w6sDwcffHDOP//8bL/99qtd4yMf+Ug++9nPZsKECXn66afzsY99LIceeuhK57/qqquy6aab5g/+4A+SJGeddVa23HLLvP3tb88555yTGTNm5EUvelGS5J3vfGc+8IEPZPfdd1/t/q3KpEmTMmvWrIwfP35gNQVUAABgqFx33XX5+te/njlz5mSzzTbLokWL8utf/3qgv+Oyyy4bSJ0TTzwxf/M3f5N58+blVa96VR544IFsssnoA12vuuqqbL311ksD6rvf/e6lj51zzjnZY489lgbUz33ucwPp37pmiC8AADBU7rvvvowfPz6bbbZZkmT8+PFLg9ukSZNy0kknZdq0aZk2bVruuOOOJMnChQtzxBFHZJ999sk+++yTa6+9Nkny+OOPZ+bMmdlzzz0zefLkXHzxxUvrLFq0KEnypS99KdOmTcuUKVPyrne9K88880yeeeaZHHvssdljjz2y55575vTTT3/WPu+2224ZO3ZsFi1aNGpf5s+fn7POOiunn356pkyZku9973v5yEc+kk984hO56KKLMmvWrLz1rW/NlClT8uSTT2b69OmZNWtWkuSCCy7InnvumT322CMnnXTS0t+59dZb54Mf/GD22muv7Lfffrn//vuf9W/x4IMPZsaMGXn5y1+ed73rXam1DuT1GklABQAAhsqMGTNy9913Z9ddd8173vOeXH311cs8vu222+YHP/hB3vve9+b9739/kuR973tfTjzxxNxwww25+OKL8853vjNJcsopp2S77bbLzTffnJtuuimvec1rlqk1b968XHjhhbn22mszd+7cjBkzJuedd17mzp2be+65J7fccktuvvnmzJw581n7fP3112eTTTbJhAkTRu3LpEmT8u53vzsnnnhi5s6dm1e96lVLlz3yyCMzderUpb93iy22WPrYvffem5NOOinf+c53Mnfu3Nxwww35t3/7tyTJL3/5y+y333658cYbc+CBB+azn/3ss/4t/v7v/z4HHHBAfvjDH+bQQw/Nz372s9V5eZ6VIb4AAMBQ2XrrrTN79ux873vfy3e/+90cddRROe2003LssccmSY4++uilP0888cQkyRVXXJFbb711aY3HHnssv/jFL3LFFVfky1/+8tL2HXbYYZnfdeWVV2b27NnZZ599kiRPPvlknve85+WNb3xj7rzzzpxwwgk55JBDMmPGjFH7evrpp+dLX/pSttlmm1x44YUppay0L6vjhhtuyPTp0zNhwoQkyVvf+tZcc801Ofzww7PpppvmDW94Q5LkFa94Rb797W8/69/immuuyde+9rUkySGHHLLC32IQBFQAAGDojBkzJtOnT8/06dOz55575txzz10aUEdeCmXJ9OLFi3Pdddctc/QxSWqtz3rplFprjjnmmHz84x9f4bEbb7wx3/zmN3PmmWfmK1/5Sr7whS+sMM+S76COtLK+rI5nG4Y7bty4pc9tzJgxefrpp1f5+9f25YMM8QUAAIbKj3/849x+++1L78+dOzc77bTT0vsXXnjh0p/7779/kmZY8Kc+9alllhmt/eGHH17mdx100EG56KKL8sADDyRJHnroodx1111ZtGhRFi9enCOOOCKnnHJK5syZ03P/V9aXbbbZZqVHUlf22L777purr746ixYtyjPPPJMLLrggr371q1fr9x944IE577zzkiT/8R//scLfYhAcQQUAANaqXi4LM0iPP/54TjjhhDzyyCMZO3ZsXvrSl+bss89e+vivfvWr7Lvvvlm8eHEuuOCCJMkZZ5yR448/PpMnT87TTz+dAw88MGeddVY+9KEP5fjjj88ee+yRMWPG5MMf/nDe9KY3La21++6759RTT82MGTOyePHijBs3LmeeeWa22GKLzJw5M4sXL06SUY+wrszK+vLGN74xRx55ZC655JL88z//8zLLHHvssXn3u9+dLbbYItddd93S9he+8IX5+Mc/nj/8wz9MrTUHH3xwDjvssNX6/R/+8Idz9NFHZ++9986rX/3qvOQlL+n5OfWqrI0zL/Vr6tSpdckZpgAAgA3bvHnzsttuu63vboxqbVy7k2WN9vqXUmbXWqeuallDfAEAAOiEVQbUUsrvlVLmjrg9Vkp5fynlOaWUb5dSbm9/7tDOX0opZ5RS7iil3FRK2XvtPw0AAIBVmz9/vqOnHbbK76DWWn+cZEqSlFLGJLknyb8mOTnJlbXW00opJ7f3T0ry+iS7tLd9k3y6/QkAAJ036eRvLHN/XX9/EjZm/Q7xPSjJT2utdyU5LMm5bfu5SQ5vpw9L8sXa+H6S7UspLxxIbwEAABha/QbUtyS5oJ1+fq31viRpfz6vbd8xyd0jllnQti2jlHJcKWVWKWXWwoUL++wGAAAAw6bngFpK2TTJoUm+uqpZR2lb4VTBtdaza61Ta61TJ0yY0Gs3AAAAGFL9XAf19Unm1Frvb+/fX0p5Ya31vnYI7wNt+4IkLx6x3MQk9655VwEAgA3SR7YbcL1HVznLRz/60Zx//vkZM2ZMNtlkk3zmM5/JvvsO7tQ4Bx98cM4///xsv/32a1xrr732yu677770mqzLmz9/ft7whjfklltuWeGx6dOn5xOf+ESmTl3lFVySJFdddVU+8YlP5Otf//oa9Xlt6SegHp3fDu9NkkuTHJPktPbnJSPa31tK+XKakyM9umQoMAAAwNp23XXX5etf/3rmzJmTzTbbLIsWLcqvf/3rgf6Oyy67bCB15s2bl8WLF+eaa67JL3/5y2y11VYDqbuh6mmIbyllyySvTfK1Ec2nJXltKeX29rHT2vbLktyZ5I4kn03ynoH1FgAAYBXuu+++jB8/PptttlmSZPz48XnRi16UJJk0aVJOOumkTJs2LdOmTcsdd9yRJFm4cGGOOOKI7LPPPtlnn31y7bXXJkkef/zxzJw5M3vuuWcmT56ciy++eGmdRYsWJUm+9KUvZdq0aZkyZUre9a535ZlnnskzzzyTY489NnvssUf23HPPnH766aP29fzzz8/b3va2zJgxI5deeunS9tmzZ2evvfbK/vvvnzPPPHNp+5NPPpm3vOUtmTx5co466qg8+eSTSx/71re+lf333z9777133vzmN+fxxx9Pklx++eV52ctelgMOOCBf+9rISNc9PQXUWusTtdbn1lofHdH2YK31oFrrLu3Ph9r2Wms9vtb6u7XWPWuts9ZW5wEAAJY3Y8aM3H333dl1113znve8J1dfffUyj2+77bb5wQ9+kPe+9715//vfnyR53/velxNPPDE33HBDLr744rzzne9MkpxyyinZbrvtcvPNN+emm27Ka17zmmVqzZs3LxdeeGGuvfbazJ07N2PGjMl5552XuXPn5p577sktt9ySm2++OTNnzhy1rxdeeGGOOuqoHH300csM8Z05c2bOOOOMXHfddcvM/+lPfzpbbrllbrrppnzwgx/M7NmzkySLFi3KqaeemiuuuCJz5szJ1KlT88lPfjJPPfVU/vIv/zL//u//nu9973v5+c9/vmZ/3LWs37P4AgAAdNrWW2+d2bNn5+yzz86ECRNy1FFH5Zxzzln6+NFHH73055IAeMUVV+S9731vpkyZkkMPPTSPPfZYfvGLX+SKK67I8ccfv3TZHXbYYZnfdeWVV2b27NnZZ599MmXKlFx55ZW588478zu/8zu58847c8IJJ+Tyyy/Ptttuu0I/b7jhhkyYMCE77bRTDjrooMyZMycPP/xwHn300TzyyCN59atfnSR529vetnSZa665Jn/+53+eJJk8eXImT56cJPn+97+fW2+9Na985SszZcqUnHvuubnrrrty2223Zeedd84uu+ySUsrSZbuqn++gAgAAbBDGjBmT6dOnZ/r06dlzzz1z7rnn5thjj02SlPLbC48smV68eHGuu+66bLHFFsvUqbUuM//yaq055phj8vGPf3yFx2688cZ885vfzJlnnpmvfOUr+cIXvrDM4xdccEFuu+22TJo0KUny2GOP5eKLL86RRx75rL9ztMdqrXnta1+7womW5s6d+6y1usYRVAAAYKj8+Mc/zu233770/ty5c7PTTjstvX/hhRcu/bn//vsnaYYFf+pTn1pmmdHaH3744WV+10EHHZSLLrooDzzQXNTkoYceyl133ZVFixZl8eLFOeKII3LKKadkzpw5yyy3ePHifPWrX81NN92U+fPnZ/78+bnkkktywQUXZPvtt892222X//zP/0ySnHfeeUuXO/DAA5fev+WWW3LTTTclSfbbb79ce+21S79T+8QTT+QnP/lJXvayl+W///u/89Of/jRJVnqm4K5wBBUAAFi7ergszCA9/vjjOeGEE/LII49k7NixeelLX5qzzz576eO/+tWvsu+++2bx4sVLA9sZZ5yR448/PpMnT87TTz+dAw88MGeddVY+9KEP5fjjj88ee+yRMWPG5MMf/nDe9KY3La21++6759RTT82MGTOyePHijBs3LmeeeWa22GKLzJw5M4sXL06SFY6wXnPNNdlxxx2z4447Lm078MADc+utt+a+++7Lv/zLv+Qv/uIvsuWWW+aP//iPl87zV3/1V5k5c2YmT56cKVOmZNq0aUmSCRMm5JxzzsnRRx+dX/3qV0mSU089NbvuumvOPvvsHHLIIRk/fnwOOOCAUS9X0xWl1rq++5CpU6fWWbOcSwkAgPVv0snfWOb+/NMOWU892XDNmzcvu+222/ruxqgmTZqUWbNmZfz48eu7K0NrtNe/lDK71rrKi7Ua4gsAAEAnGOILAABsNObPn7++u8CzcAQVAAAYuC58lZB1b01fdwEVAAAYqM033zwPPvigkLqRqbXmwQcfzOabb77aNQzxBQAABmrixIlZsGBBFi5cuL67wjq2+eabZ+LEiau9vIAKAAAM1Lhx47Lzzjuv726wATLEFwAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBO6CmgllK2L6VcVEq5rZQyr5SyfynlOaWUb5dSbm9/7tDOW0opZ5RS7iil3FRK2XvtPgUAAACGQa9HUP8pyeW11pcl2SvJvCQnJ7my1rpLkivb+0ny+iS7tLfjknx6oD0GAABgKK0yoJZStk1yYJLPJ0mt9de11keSHJbk3Ha2c5Mc3k4fluSLtfH9JNuXUl448J4DAAAwVHo5gvo7SRYm+ZdSyg9LKZ8rpWyV5Pm11vuSpP35vHb+HZPcPWL5BW3bMkopx5VSZpVSZi1cuHCNngQAAAAbvl4C6tgkeyf5dK315Ul+md8O5x1NGaWtrtBQ69m11qm11qkTJkzoqbMAAAAMr14C6oIkC2qt17f3L0oTWO9fMnS3/fnAiPlfPGL5iUnuHUx3AQAAGFarDKi11p8nubuU8ntt00FJbk1yaZJj2rZjklzSTl+a5O3t2Xz3S/LokqHAAAAAsDJje5zvhCTnlVI2TXJnkplpwu1XSinvSPKzJG9u570sycFJ7kjyRDsvAAAAPKueAmqtdW6SqaM8dNAo89Ykx69hvwAAANjI9HodVAAAAFirBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATegqopZT5pZSbSylzSymz2rbnlFK+XUq5vf25Q9teSilnlFLuKKXcVErZe20+AQAAAIZDP0dQ/7DWOqXWOrW9f3KSK2utuyS5sr2fJK9Pskt7Oy7JpwfVWQAAAIbXmgzxPSzJue30uUkOH9H+xdr4fpLtSykvXIPfAwAAwEag14Bak3yrlDK7lHJc2/b8Wut9SdL+fF7bvmOSu0csu6BtAwAAgJUa2+N8r6y13ltKeV6Sb5dSbnuWecsobXWFmZqge1ySvOQlL+mxGwAAAAyrno6g1lrvbX8+kORfk0xLcv+Sobvtzwfa2RckefGIxScmuXeUmmfXWqfWWqdOmDBh9Z8BAAAAQ2GVAbWUslUpZZsl00lmJLklyaVJjmlnOybJJe30pUne3p7Nd78kjy4ZCgwAAAAr08sQ3+cn+ddSypL5z6+1Xl5KuSHJV0op70jysyRvbue/LMnBSe5I8kSSmQPvNQAA0JdJJ39jhbb5px2yHnoCK7fKgFprvTPJXqO0P5jkoFHaa5LjB9I7AAAANhprcpkZAAAAGBgBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE7oOaCWUsaUUn5YSvl6e3/nUsr1pZTbSykXllI2bds3a+/f0T4+ae10HQAAgGHSzxHU9yWZN+L+PyQ5vda6S5KHk7yjbX9HkodrrS9Ncno7HwAAADyrngJqKWVikkOSfK69X5K8JslF7SznJjm8nT6svZ/28YPa+QEAAGClej2C+o9J/o8ki9v7z03ySK316fb+giQ7ttM7Jrk7SdrHH23nBwAAgJVaZUAtpbwhyQO11tkjm0eZtfbw2Mi6x5VSZpVSZi1cuLCnzgIAADC8ejmC+sokh5ZS5if5cpqhvf+YZPtSyth2nolJ7m2nFyR5cZK0j2+X5KHli9Zaz661Tq21Tp0wYcIaPQkAAAA2fKsMqLXW/1lrnVhrnZTkLUm+U2t9a5LvJjmyne2YJJe005e299M+/p1a6wpHUAEAAGCkNbkO6klJPlBKuSPNd0w/37Z/Pslz2/YPJDl5zboIAADAxmDsqmf5rVrrVUmuaqfvTDJtlHmeSvLmAfQNAACAjciaHEEFAACAgRFQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6IRVBtRSyuallB+UUm4spfyolPL3bfvOpZTrSym3l1IuLKVs2rZv1t6/o3180tp9CgAAAAyDXo6g/irJa2qteyWZkuR1pZT9kvxDktNrrbskeTjJO9r535Hk4VrrS5Oc3s4HAAAAz2qVAbU2Hm/vjmtvNclrklzUtp+b5PB2+rD2ftrHDyqllIH1GAAAgKHU03dQSyljSilzkzyQ5NtJfprkkVrr0+0sC5Ls2E7vmOTuJGkffzTJcwfZaQAAAIZPTwG11vpMrXVKkolJpiXZbbTZ2p+jHS2tyzeUUo4rpcwqpcxauHBhr/0FAABgSPV1Ft9a6yNJrkqyX5LtSylj24cmJrm3nV6Q5MVJ0j6+XZKHRql1dq11aq116oQJE1av9wAAAAyNXs7iO6GUsn07vUWSP0oyL8l3kxzZznZMkkva6Uvb+2kf/06tdYUjqAAAADDS2FXPkhcmObeUMiZNoP1KrfXrpZRbk3y5lHJqkh8m+Xw7/+eT/O9Syh1pjpy+ZS30GwAAgCGzyoBaa70pyctHab8zzfdRl29/KsmbB9I7AAAANhp9fQcVAAAA1hYBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThi7vjsAAMCGZ9LJ31ihbf5ph6yHngDDxBFUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOmGVAbWU8uJSyndLKfNKKT8qpbyvbX9OKeXbpZTb2587tO2llHJGKeWOUspNpZS91/aTAAAAYMPXyxHUp5P8j1rrbkn2S3J8KWX3JCcnubLWukuSK9v7SfL6JLu0t+OSfHrgvQYAAGDorDKg1lrvq7XOaad/kWRekh2THJbk3Ha2c5Mc3k4fluSLtfH9JNuXUl448J4DAAAwVPr6DmopZVKSlye5Psnza633JU2ITfK8drYdk9w9YrEFbRsAAACsVM8BtZSydZKLk7y/1vrYs806Slsdpd5xpZRZpZRZCxcu7LUbAAAADKmeAmopZVyacHperfVrbfP9S4butj8faNsXJHnxiMUnJrl3+Zq11rNrrVNrrVMnTJiwuv0HAABgSPRyFt+S5PNJ5tVaPznioUuTHNNOH5PkkhHtb2/P5rtfkkeXDAUGAACAlRnbwzyvTPK2JDeXUua2bX+X5LQkXymlvCPJz5K8uX3ssiQHJ7kjyRNJZg60xwAAAAylVQbUWut/ZvTvlSbJQaPMX5Mcv4b9AgAAYCPT11l8AQAAYG0RUAEAAOgEARUAAIBOEFABAADoBAEVAACATujlMjMAwBqYdPI3lrk//7RD1lNPAKDbHEEFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6QUAFAACgEwRUAAAAOkFABQAAoBMEVAAAADpBQAUAAKATBFQAAAA6Yez67gAAJMmkk7+xQtv80w5ZDz0BANYXR1ABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4QUAEAAOgEARUAAIBOEFABAADoBAEVAACAThBQAQAA6AQBFQAAgE4Yu747AAAbnY9sN0rbo+u+HwDQMY6gAgAA0AkCKgAAAJ0goAIAANAJAioAAACdIKACAADQCQIqAAAAnSCgAgAA0AkCKgAAAJ0goAIAANAJAioAAACdIKACAADQCQIqAAAAnSCgAgAA0AkCKgAAAJ2wyoBaSvlCKeWBUsotI9qeU0r5dinl9vbnDm17KaWcUUq5o5RyUyll77XZeQAAAIZHL0dQz0nyuuXaTk5yZa11lyRXtveT5PVJdmlvxyX59GC6CQAAwLBbZUCttV6T5KHlmg9Lcm5+DQ4FAAAbGElEQVQ7fW6Sw0e0f7E2vp9k+1LKCwfVWQAAAIbX6n4H9fm11vuSpP35vLZ9xyR3j5hvQdu2glLKcaWUWaWUWQsXLlzNbgAAADAsBn2SpDJKWx1txlrr2bXWqbXWqRMmTBhwNwAAANjQrG5AvX/J0N325wNt+4IkLx4x38Qk965+9wAAANhYrG5AvTTJMe30MUkuGdH+9vZsvvsleXTJUGAAAAB4NmNXNUMp5YIk05OML6UsSPLhJKcl+Uop5R1Jfpbkze3slyU5OMkdSZ5IMnMt9BkAAIAhtMqAWms9eiUPHTTKvDXJ8WvaKQBgHfvIdsvdf3T99AOAjdqgT5IEAAAAq2WVR1ABAAA2OEaGbJAcQQUAAKATBFQAAAA6wRBfAIbf8sO8EkO9gA3GpJO/sULb/NMOWQ89gbVPQAWgu3x/CAA2KhtOQLX3GwAAYKhtOAEVAOg+O5QBWANOkgQAAEAnOIIKAACsF8ufAMrJn3AEFQAAgE4QUAEAAOgEQ3wBoFdOAAT0wfVLWSuG/L3IEVQAAAA6wRFUAACAlVn+iOUQHa3sIkdQAQAA6ARHUAEAADpqhUvxbL6eOrKOCKgA64shQwAAy9j4AuqQn/UKAIaC92uAjdLGF1CBjZMPuwAAnSegAjBURr3u4JB/XwcAhoWz+AIAANAJjqACwEpsbGdOBDZCvgJDxziCCgAAQCc4ggqwobP3GwAYEgIqAADQDXa6bvQM8QUAAKATHEFl42KvHAAAdJYjqAAAAHSCI6gADN7yoxWMVGAYWK8B1joBFVaHocIAADBwhvgCAADQCY6gsnYZDsUwsl4DAP0yAq8nAioAa2TSyd9YoW3+5uuhIwDABk9ABdbYqAHltEMGU9zexnXH3xo6bflt7cC2s8BSdrqufwIqAMByfEhdd9bqTk5ggyOgAgDDzffGgfVghVEPdnL1REAFAFhLHB1ctwyDhg2fy8wAAADQCY6gAsAGwtG4VfPdUYANm4AKAKw237Ha8BgGC3SZgMqGwQkuAABg6AmoAABAX1YcPfFnK87kgAKrQUAFgA3ZkIwwMVQYgERAhfVr+Q+WyQb74RIAOsf7LGxwXGYGAACATnAEFegUw/wAetTFo4NDMuQcWH8EVFbUxTc81p0h+XAh6MLKuVboeuZ9dtWG5L0IumT0bf9yJ7fqwP+agLoBcx0zAABgmAioXWAvIcPIEQIAYB0xMmR4rJWAWkp5XZJ/SjImyedqraetjd8D68qwDhe1MQeADZQDHKwno35+HOBIzoEH1FLKmCRnJnltkgVJbiilXFprvXXQvws2NGs16Hqj2mgM6w4TAGADNcDPoWvjCOq0JHfUWu9MklLKl5MclmT4AqpAwHriyCeDIOgCSXwlY13z+RGeVam1DrZgKUcmeV2t9Z3t/bcl2bfW+t7l5jsuyXHt3d9L8uNVlB6fZNEAujioOoOsNax1BllLnXVXS511V2tY6wyyljrrrtaw1hlkLXXWXa1hrTPIWuqsu1rDWmeQtXqps1OtdcIqK9VaB3pL8uY03ztdcv9tSf55AHVnDah/A6nTxT51rU4X+zSsdbrYp2Gt08U+da1OF/s0rHW62Keu1elin4a1Thf71LU6XezTsNbpYp+6Vqerfdpk5dF1tS1I8uIR9ycmuXct/B4AAACGyNoIqDck2aWUsnMpZdMkb0ly6Vr4PQAAAAyRgZ8kqdb6dCnlvUm+meYyM1+otf5oAKXPHkCNQdYZZK1hrTPIWuqsu1rqrLtaw1pnkLXUWXe1hrXOIGups+5qDWudQdZSZ93VGtY6g6w1sD4N/CRJAAAAsDrWxhBfAAAA6JuACgAAQCcIqAAAAHSCgAoAAEAnCKgAAKMopWxXSjmqlPKBUsqJ7fT2A6z/2j7n37aU8rujtE/us84LSikvaKcnlFLeVEr5/X5qrKTux9a0Rltn57ZPL+tzuZeUUjZvp0spZWYp5Z9LKX9VSun5yhWllEOX1FlTpZQDSym/104fUEr5m1LKIatZa+tSypHtunhCKeV1pZS+P8uXUl5WSjmplHJGKeWf2undVqdPK6k/czX6c1ApZevl2l/XZ51ppZR92und2//bg/up8Sy1vziAGge0fZrR53L7llK2bae3KKX8fSnl30sp/1BK2a6POn9dSnlxv/1eSa1NSylvL6X8UXv/z0opnyqlHF9KGbfG9bt4Ft9SyieTXFxrvXZA9f44yeFJdkxSk9yb5JJa6+WrUWdikitrrfNHtP9FrfULPdYoSd7c9uOiJK9JcliS25KcVWtd3E+flqv9nVrra/pcZnytddGI+3+eZFqSW5J8tvaxgpRS/iTJ1bXWh0opE5L8v0lenuTWJP+j1rqgxzoDef1LKc9J8t40r/fnk/xdkv2TzEvysVrrw33U+sMkRyR5cZKnk9ye5HO11jtWo18vS/Oaj1wfL621zuu31krqz6y1/kuf/dkxyfW11sdHtL+un/+RUsq0JLXWekMpZfckr0tyW631sj66P1rdL9Za374mNdo6B6Rdt2ut3+pjuX2TzKu1PlZK2SLJyUn2TrNef6zW+miPdf46yb/WWu/uv/fL1Flyfel7a61XlFL+LMkfpFmvz661/qaPWr+b5E+y7Hp9Qa/PablaA9nOPkv9/6vW+r/67E9nt9dt/fW2zR7m7XVbb4232aWUtyf5cJJvJbmnbZ6Y5LVJ/r7WOogPrD+rtb6kx3n/NMk/Jnkgybgkx9Zab2gfm1Nr3bvHOu9Ksx0rSf4hybFJfpTklUn+71rr53usc8byTUneluSLSVJr/ete6rS1/q3Weng7fVia53lVmm3bx2ut5/RY55Yk02qtT5RS/iHJ7yb5tzT/u6m1/kWPdZ5M8ssk/5HkgiTfrLU+0+vzGVHnH9P8f45Nc/nFg9qar07yw1rr3/ZR60+T/G2SG5P8YZL/SnOgac8kb6213txjnZOSHJ3ky0mW/J9PTPO+8uVa62m99ulZfkc/6/VfJzk+zf/6lCTvq7Ve0j7Wz3r94SSvT/O3/naSfdOsQ3+U5vX7aB/9v3T5pjR/8+8kSa310B7r/KDWOq2d/ss0z/Nfk8xI8u+9/q1LKT9Ksld7Kc+zkzyR5j3poLb9TT3WeTTNev3TNOv1V2utC3tZdpRa56X5W2+Z5JEkWyf5WtunUms9ZnXqLq3f0YC6MMldSSYkuTDNB6Yfrmatf0yya5oN5sh/xLcnub3W+r4e63wsyQFJ5iR5Y5J/rLX+c/tYP/9A/1+S5yXZNMljSTZL8u9JDk5yfx/9uWn5pjTP88dJUmvtaW/qyL6XUj6U5FVJzk/yhiQLaq0n9lKnXf7WWuvu7fSFSb6f5KtpNg5vrbX2tKd4UK9/KeWyJDcn2TbJbu30V9J8uNir1npYj3VOS/L8JFem+QD+30l+kuQ9aT44fbWPPnljePY6A3lTaGsN5RvDoN4U2tf+jUmuTrP9mZvk4TSB9T211qv66NNAtrOr+B39rNed2l63tTq1zR7W7XVbayDb7FLKj5PsW2t9ZLn2HdLs0Nu1xzrLb9eWPpTkNbXWrXqsMzfJ62ut97U7BL+Y5O9qrV8rpfyw1vryHuvcnGYbvUWa1+6ltdaft8/ru7XWKT3WWZBmO/+t9rkkySeS/E2S1FrP7aVOW2tp/0sp/5VmHfzvUsr4NDuZ9uqxzsj1enaSfZbsSCql3NhHnR+mCbVHpnl/3iPN+8cFtdar+3heP2qX3SLNTo4d2/A8Lk1A3aOPWjcl2a9dfnyS82qtf1yao+dn1Vr/oMc6P0ny+8vvzGx3fv6o1rpLH/0Z9aEku9ZaN+uxzs1J9q+1Pl5KmZTm/fV/11r/aTXW6ylpttM/TzJxxI7l63vdxra15qTZWfe5NDsoS5r37bckSa/rwHLr9Q1JDq61LiylbJXk+7XWPXusM6/WutuSvo18DyulzO3jf/aHSV6RZjt/VJJDk8xun9vXaq2/6KVOW+umWuvk0oxMuCfJi2qtz7Q7dm/s5+89qlpr525p/mmTZJck/2eaPXu3pdmTuWuftX6ykvaS5oNTr3VuTjK2nd4+yWVJTh/Z317rtD/HJXkwyabt/bFLHuuxzqVJvpTkZUl2SjIpyd3t9E79/q3b6TlJthrRv5770y7z4xHTs5d7bO66fv2X/M72tb5nDfpz84jpsUmubad3SHM0rq/1Mcm4Udo37XN9vGklt5uT/KrP9XrrdnpSkllpQmrf63WSMWlC02NJtm3bt0hyUx915rTr9fQ0e5inJ7mvnX51n3/rkev2DUkmtNNb9fm/Nm9k/9ZkvU6zt3tGmiNEC5NcnuSYJNv089qPWBfvTzJmxHrez9/65hHLbpnkqnb6Jf289kvW65W097udfWwlt18kebrP59aZ7XW7TKe22RnS7fXI123Ea7Va2+w02+vtRmnfrs/1+uEkh6Tdjo24TU+zo6Pv59Xef2GaD5d/neW2TauoM2fE9I0rW796qLNNmiOd56cJXklyZz+v1Ur69IM16NM304T+JLl4yf9Wkucu/1x77U97/wXt3/m6JHf3UeeW9ufm7XqwRXt/TJJb+/wb3ZzfHljaYrltQT/r9W2jbXPabdGP+6hzf5pAuNNyt0lpRvj0WufW5e5vnea98ZOrsy0abZ1ZjW3IJklOTLPDfUrb1ve6neZo9w7t+jdrZf3toc5Xk8xsp/8lydR2etckN/RRZ/n1elyakHpBkoV9Prdb0nx23SHNe/Rz2vbNM+Jz0+reeh6Pv47VJKm13p7klCSntHuIjk7zQeOlfdR6qpQyrdb6g+Xa90nyVB91xtZan2779Ugp5Y1Jzi6lfDXNC9SrJTV+U0q5odb66/b+06WUnoeP1FoPbYdonZ3kE7XWS0spv6m13tVHX5Jki1LKy9P8M46ptf5yRP/6Hc5yVSnlfyX5eDt9eK3139qhVv0MGRzU679Ju0d4myRbl1Im1Vrnl1Kem/5es8WllOfUWh9K8qI0byyptT7c7inqx+K2xvKv0wvbx3r1/CR/nOYNb6SSZthPr8bUdlhv+7eZnuSiUspO+e0e8V48XZvhT0+UUn5aa32srflkKaWf5zU1yfuSfDDJ39Za55ZSnqx97LEeYcnrv0maN/WFbZ9+WUp5uo86t5TfDpu+sZQytdY6q5Sya5Keh9M2v7ouTnO04VvtHvTXp1mvP5HmCFQvNmn3dG+VJlhul+ShNHuN+/3ex9gkz7TLbtN28mel/++PDGo7+0iaIx73L/9AKaWfodGd2l63y3Rtmz2s2+tkcNvsjyaZU0r5VpqdCUmzA+e1aZ5rr76f5InRtmPtUdpe/aKU8ru11p8mSW2OpE5PM4S1n++PLi6ljKvNEbSl34UszXcue/4+Y22Otry/lPKKJF8qpXyjn+WXs1cp5bE07zublVJeUJujupumff169M4kXyylfCTNejy3PWq0Q5IP9FFnmfWk1vrzJGckOaN9f+zVN0op30vzof1zSb5SSvl+mh0U1/RRJ2n+ny4vpVyd5r3jq8nS4fH9rNfvT3JlKeX2LLtevzTNMPtefT3NDu65yz9QSrmqjzo/L6VMWVKnNkdS35DkC2mGL/fq16WULWutT6Q5SrikL9ulv89Xad+rT2/fM04vpdyfrFZm2i7NTqSSpI5Yr7dOf6/ZO5P8UztqZlGS69r3xLvbx3q1/Hr9mzQ7Ty9tjzT34/NpdnaMSfOZ7aullDuT7JdmlOCaWdOEuzZu6XPv/SpqvSLJ9WkO1X+rvc1r217RR52vZ5QjOElOTbK4jzr/kfaI1XLtL8hyew17rLdVmr1Ml6YZ3tXv8lcl+e6I2wvb9hX29vRQa1ySjyT5WXtbnGavyvlJXrKuX/80H5Dub29HJLkizd6we5Ic10edo9IEym+1z+uQtn1CkvP77NPrktzRrgdnt7fL27bX9VHn80kOWMljPfcpzdDZKcu1jU0zdOyZPupcn2TLdnqTEe3bpY89+yOWm5jmzfdTSX62mq///CR3phned2eSF7TtW6e/vbLbJTknzdDc69OE0jvTDI3daxDrddq96j3WObH9/Xel2aN/ZZLPptm7/uE+6rwvzVH3s9O8ySzZOzshyTV9/q33zmC2s6em+f7YaI/9Qx91Orm9bpdd0232dzOAbXaGdHvd1hrkNnuHNMP6/keaoatvSbLDIJ7zavyN9kqyy0pey7f2UeclGX0kz45J/mg1+1bSfI3iSwN+ztunGf7Z73K7pfm++BFphjNv0ufy0wf4HPZPMzQ3ab4T+zdJ/rTfPrXLH9wu/9oRbZsk2azPOpukCRJHpBnGvF/aETXr+pbm/f4FK3nslX3UGfVvkGR8kj3XsI+HpPl6wKCe85ZJdl6N5bZptwOvSPL81Vi+r1GoPdR7UZqhvUv+V4/MSt7D+7119TuoW9cRJ2sZUM0XpNn4ljQfCn7e5/JbJM0RoVEe27HWes+KS/VVf6s0Q7UeWM3l90qzET9rTfoxot6YNP/sT6zm8tulOYrx4GosO7DXv30epTZHPMamGY5yT631vj7rPCfJ7yS5oy73faTV6NMmaU6asHR9TDNEo+8TMKypUsrENEc/V/h/KKW8svZ44pNSyma11l+N0j4+zQfonk7eMMryh6R5g/q71Vl+JTW3TLNh/+8+l9smzTowNs02ZIWjfKtYftda60/6WeZZar0oSWqt95bmjKJ/lCbIL38Ec1V1fj/NB7lbaq23DaBfa7SdHZSub6/bGp3ZZg/b9rqtNcht9vMz4uRf/f7vq7Ph9qlrdQZda5TaA/l/HtY6XexT1+oMqlZXA+qmSX5T2861Q472TjNO/T/6rDW51rqyL3KrM6A6XexT1+qMqPeSJI/VZujhpDTDWufVWn80gDq31VpvUaebfRrWOiPqTc2Is6aubugd1jpd7NOw1hlErVLKlCRnpRlFsSDNjpeJaYaiv6fWOqfHOi9P8um2zsizAfdbZ2R/lq/zV7XHk1Otozo9P6+21rP9jdZ5n7pWZ9C1nuV39Hwyuo2xThf71LU6g6rV1YB6Y5rhFQ+XUv42zVklL0szZn9WrfV/9lHrmTRD/C5Ic/a1W1ezT+psYH3qWp221slJ3pXkV/nt2Q6vTTO85vO11k+qs+Z1utinYa3T1np1msuUPJJm6NG1aYZG/ibJ22qPl9YZ1jpd7NOw1hlwn+YmeVet9frl2vdL8pna+xlh1dnA+tS1OgPu08q+i1uSfLDW+pyNuU4X+9S1OoOuNWqRjgbUW2p76u1Syqwkr6rNyVbGpvk+Wz+niv5hmutyHZ3meym/TBM0vlxHXBtPnTWr08U+da1OW+tHaY5SbZnmO5K/U397yvHra4+nnFdnw+vTsNZpa/0wyYx2+Z2TfLLW+iellNemOdlVTxclH9Y6XezTsNYZcJ9uryu55EYp5Y5aa08ngFJnw+tT1+oMuE9PJfl/0p4Ebjkn1lq335jrdLFPXasz6FqjqgP8suygbmnOQrpHO3152hMSpDkLWr+X9Vj+lMrT0pyg4u4k/6XOYOp0sU9dq9Muu+QSIWPSXGx95AmF+jlNvDobWJ+Gtc7IWiPqjbxkxI829jpd7NOw1hlwn85I8o00Oyb/oL0d1bZ9Sp3B1Olin7pWZ8B9+q+s5OR16e8SOkNZp4t96lqdQdca7dbVI6iTk/zvNNcPSpJXpjlj5uQ0e0LP76PWqBf5LaWUJAfW1bjYrjobRp+6Vqdd5pw0l0zYKskTafY8XZ7mguDb1Fr/VJ01r9PFPg1rnbbWF9KcsOPKNGfOvKfW+oHSnJBqTq31ZRtznS72aVjrrIVar29rjDyp3aW11st6raHOhtmnrtUZVK1Syu8leai2l15b7rHn1x5PujSsdbrYp67VGXStUet3MaAmSWnO5jcjzUVox6b5J/xm7fNsfKWUP+sn0Kqz/msNa5221tgkb07z4emiNKe/PzrNpRDOrO01DdVZszpd7NOw1mlrjUvyl0l2T7Nj8Qu11mdKczbd59Uer/U5rHW62KdhrTPoWgCse50NqAAA60tpLr/zP9McsXpe2/xAkkuSnNbrDnN1Nrw+da3OWurT4WmuDaxOx/vUtTqDrjWaTdZk4bWllLJ1KeV/lVJ+VEp5tJSysJTy/VLKseurljr/f3t3zxpFHEVh/LlqaWWQWAhaqY1gb4oEBC0UOy3FLxA7BfGraCmigtpKCgWxUPEtNumClVXSxNa/xWwRdjfIsnd3z1zPA0PYHfhxmhCG2c30b5Oa8w/rlp08R3FTVWfI+pH0u1bKUdxU1Um2ngK7wFprbam1tgSs0f134Gd20hzFTWrOLDatDjm7dmQ3qTnZ1kiSd1Aj4hXwAtgAbtB9R+oJ8IDuuyT3523Z6d8mNUdxU1VHcVNVR3GTmqO4qaqTvGmrtXZ20nN2JnMUN6k5ipuqOoqb1Jxsa2xtyv+yNIsD+Db0+uPg5yG6h8jP3bLTv01qjuKmqo7ipqqO4iY1R3FTVSd502vgLrC8771l4B6wYSfHUdyk5ihuquooblJzsq1xh+RHfIHfEbECEBHXgB2A1tofIBZk2enfJjVHcVNVR3FTVUdxk5qjuKmqk2ndBJaAtxGxGxE7wBvgGN2dWTs5juImNUdxU1VHcZOak22NNu0V7iwOusfJfKD7XP074Mzg/ePA+iIsO/3bpOYobqrqKG6q6ihuUnMUN1V1ZmCdAy4BR4fev2Inz1HcpOYobqrqKG5Sc7KtEXtaYN4HcFvNstO/TWqO4qaqjuKmqo7iJjVHcVNVZ1ILWAe2gJfANnB937nPdnIcxU1qjuKmqo7iJjUn2xrrTwvM+wB+qll2+rdJzVHcVNVR3FTVUdyk5ihuqupMagGbDO4MAKeBT8CdwesvdnIcxU1qjuKmqo7iJjUn2xp3HEGwiPh+0Cm6L+DO3bLTv01qjuKmqo7ipqqO4iY1R3FTVSfZOtxa2wNorW1HxCrwPCJOMdl3We30b5Oao7ipqqO4Sc3JtkaSvECl+wNyme5ZOvsL4P2CLDv926TmKG6q6ihuquooblJzFDdVdTKtXxFxobX2FaC1thcRV4FHwHk7aY7iJjVHcVNVR3GTmpNtjTbtLdhZHMBDYOWAc48XYdnp3yY1R3FTVUdxU1VHcZOao7ipqpO86SRw4oBzF+3kOIqb1BzFTVUdxU1qTrY17ogB5JxzzjnnnHPOLTTV56A655xzzjnnnPvP8gWqc84555xzzjmJfIHqnHPOOeecc04iX6A655xzzjnnnJPIF6jOOeecc8455yT6C9QgF1s7ShPPAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 1152x576 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "df_merged.plot.bar();"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Conclusions\n",
    "\n",
    "There seem to be a couple broad trends we can draw from all our analyses.\n",
    "\n",
    "- Looking at the numbers of petitions per year, it seems that petitioning to add species to the ESA has gotten more common over time. There are few petitions before 1990, and what seem to be a significant amount after 2000. \n",
    "\n",
    "\n",
    "- Looking at the species added per year, the data seem to be multimodal. Perhaps this could reflect the political climate of the time, and the executive branch's views on protecting endangered and threatened species?\n",
    "\n",
    "And there are more analyses we could pursue, if we wanted to delve more deeply:\n",
    "\n",
    "- Explore the actual `Petition Finding(s)` from the petition data. How successful are petitions?\n",
    "\n",
    "\n",
    "- Do petitions or FWS scientists have a greater impact on getting species listed under the ESA?\n",
    "\n",
    "\n",
    "- How do the types of species (plant, animal, invertebrate, etc.) affect the success of a petition, or the likelihood that a species will be listed?\n",
    "\n",
    "Feel free to clone this notebook and do some of these further analyses, if you are interested in learning more!"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}