00_Presentation.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h2><center> Evaluating the effect on marine biodiversity of Marine Protected Areas (MPAs) worldwide: </center></h2><h3><center> Presentation on Ocean Biodiversity Information System (OBIS) MPA Project </center></h3>\n",
    "<center>Abby Benson, Rene Francolini, Natalie Posdaljian, Natalie Low, Kyle Oliveira, Camille Ross, Sam McNeely, Nick Record, Camille Ross, Angel Segura, Ben Tupper </center>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "![MPA](MPA.png)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<b> OBIS </b>\n",
    "OBIS is a global open-access data and information clearing-house on marine biodiversity for science, conservation and sustainable development\n",
    "\n",
    "<b> MPAs </b>\n",
    "Conserved waters for a number of reasons including economic resources, biodiversity conservation, and species protection\n",
    "\n",
    "<b> Standardized Species Diversity: ES50 </b>\n",
    "Hurlbert's index, or ES50, is the expected number of distinct species in a random sample of 50 observations; it is a bias-indepenent indicator on marine biodiversity richness. It has the ability to somewhat correct for sampling bias."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Linking to GEOS 3.8.1, GDAL 3.1.2, PROJ 7.1.0\n",
      "\n",
      "── \u001b[1mAttaching packages\u001b[22m ─────────────────────────────────────── tidyverse 1.3.0 ──\n",
      "\n",
      "\u001b[32m✔\u001b[39m \u001b[34mggplot2\u001b[39m 3.3.2     \u001b[32m✔\u001b[39m \u001b[34mpurrr  \u001b[39m 0.3.4\n",
      "\u001b[32m✔\u001b[39m \u001b[34mtibble \u001b[39m 3.0.3     \u001b[32m✔\u001b[39m \u001b[34mdplyr  \u001b[39m 1.0.1\n",
      "\u001b[32m✔\u001b[39m \u001b[34mtidyr  \u001b[39m 1.1.1     \u001b[32m✔\u001b[39m \u001b[34mstringr\u001b[39m 1.4.0\n",
      "\u001b[32m✔\u001b[39m \u001b[34mreadr  \u001b[39m 1.3.1     \u001b[32m✔\u001b[39m \u001b[34mforcats\u001b[39m 0.5.0\n",
      "\n",
      "── \u001b[1mConflicts\u001b[22m ────────────────────────────────────────── tidyverse_conflicts() ──\n",
      "\u001b[31m✖\u001b[39m \u001b[34mdplyr\u001b[39m::\u001b[32mfilter()\u001b[39m masks \u001b[34mstats\u001b[39m::filter()\n",
      "\u001b[31m✖\u001b[39m \u001b[34mdplyr\u001b[39m::\u001b[32mlag()\u001b[39m    masks \u001b[34mstats\u001b[39m::lag()\n",
      "\n",
      "Loading required package: permute\n",
      "\n",
      "Loading required package: lattice\n",
      "\n",
      "This is vegan 2.5-6\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# install all libraries needed\n",
    "\n",
    "library(sf)\n",
    "library(robis)\n",
    "library(tidyverse)\n",
    "library(ggplot2)\n",
    "library(rnaturalearth)\n",
    "library(RColorBrewer)\n",
    "library(rnaturalearthdata)\n",
    "library(vegan)\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Download MPA polygons from NOAA"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "How to download NOAA MPAs"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "#### How to download NOAA MPAs, not to be run during presentation ####\n",
    "## Refer to mpa_poly.Rmd & get_mpa_poly.R"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Choose MPA polygon(s) of interest"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Given the MPA of interest, extract data from OBIS within the MPA"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# OBIS request function occurrence() in robis requires a text-formatted \n",
    "# geometry input this is a simple helper function to convert a standard \n",
    "# polygon shapefile into a suitable text-formatted geometry\n",
    "\n",
    "extract_polygon_geometry <- function(polygon) {\n",
    "  polygon %>%\n",
    "    st_geometry(.) %>%\n",
    "    st_as_text(.) %>%\n",
    "    return(.)\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reading layer `ca_all_mpas' from data source `/home/jovyan/ohw20-proj-species-marine-protected-areas/shapefiles/california_mpas/ca_all_mpas.geojson' using driver `GeoJSON'\n",
      "Simple feature collection with 155 features and 14 fields\n",
      "geometry type:  MULTIPOLYGON\n",
      "dimension:      XY\n",
      "bbox:           xmin: -124.5326 ymin: 32.53284 xmax: -117.1238 ymax: 42.00004\n",
      "geographic CRS: WGS 84\n"
     ]
    }
   ],
   "source": [
    "# Load sample polygons - California MPA file\n",
    "# Can replace this with global polygons when those are available\n",
    "\n",
    "ca_mpa_polygons <- st_read(\"shapefiles/california_mpas/ca_all_mpas.geojson\")\n",
    "\n",
    "\n",
    "# subset a couple of MPA polygons\n",
    "\n",
    "# Soquel Canyon is a little more offshore\n",
    "soquel_canyon_smca <- ca_mpa_polygons %>% filter(MPA_NAME == \"Soquel Canyon SMCA\")\n",
    "\n",
    "# Point Lobos is more coastal but an older MPA\n",
    "point_lobos_smca <-ca_mpa_polygons %>% filter(MPA_NAME == \"Point Lobos SMCA\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "![map](CalMap.jpg)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Get all species in MPA(s) from OBIS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Retrieved 5000 records of approximately 6477 (77%)\n",
      "Retrieved 6477 records of approximately 6477 (100%)\n",
      "Retrieved 1600 records of approximately 1600 (100%)\n"
     ]
    }
   ],
   "source": [
    "## pull all OBIS records for two California MPAs ##\n",
    "\n",
    "soquel_canyon_records <- occurrence(geometry = extract_polygon_geometry(soquel_canyon_smca))\n",
    "point_lobos_records <- occurrence(geometry = extract_polygon_geometry(point_lobos_smca))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A tibble: 6 × 148</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>brackish</th><th scope=col>date_year</th><th scope=col>references</th><th scope=col>scientificNameID</th><th scope=col>scientificName</th><th scope=col>coordinateUncertaintyInMeters</th><th scope=col>absence</th><th scope=col>individualCount</th><th scope=col>dropped</th><th scope=col>aphiaID</th><th scope=col>⋯</th><th scope=col>parvorder</th><th scope=col>parvorderid</th><th scope=col>organismID</th><th scope=col>category</th><th scope=col>organismQuantity</th><th scope=col>organismQuantityType</th><th scope=col>infrakingdomid</th><th scope=col>subkingdom</th><th scope=col>infrakingdom</th><th scope=col>subkingdomid</th></tr>\n",
       "\t<tr><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>⋯</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td> TRUE</td><td>2003</td><td>https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]</td><td>urn:lsid:marinespecies.org:taxname:558   </td><td>Porifera        </td><td>500m</td><td>FALSE</td><td>1   </td><td>FALSE</td><td>   558</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>FALSE</td><td>2000</td><td>NA                                                                                                                                                                    </td><td>urn:lsid:marinespecies.org:taxname:124713</td><td>Psolus squamatus</td><td>NA  </td><td>FALSE</td><td>NA  </td><td>FALSE</td><td>124713</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>   NA</td><td>1991</td><td>http://calcofi.org/publications/calcofireports/v57/Vol57-Sakuma_pages.163-183.pdf                                                                                     </td><td>urn:lsid:marinespecies.org:taxname:1128  </td><td>Euphausiacea    </td><td>NA  </td><td>FALSE</td><td>7150</td><td>FALSE</td><td>  1128</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>   NA</td><td>2007</td><td>https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]</td><td>urn:lsid:marinespecies.org:taxname:125267</td><td>Holaxonia       </td><td>50m </td><td>FALSE</td><td>3   </td><td>FALSE</td><td>125267</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>FALSE</td><td>2002</td><td>https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]</td><td>urn:lsid:marinespecies.org:taxname:132046</td><td>Polymastia      </td><td>500m</td><td>FALSE</td><td>2   </td><td>FALSE</td><td>132046</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>   NA</td><td>2000</td><td>NA                                                                                                                                                                    </td><td>urn:lsid:marinespecies.org:taxname:135392</td><td>Nanomia         </td><td>NA  </td><td>FALSE</td><td>NA  </td><td>FALSE</td><td>135392</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A tibble: 6 × 148\n",
       "\\begin{tabular}{lllllllllllllllllllll}\n",
       " brackish & date\\_year & references & scientificNameID & scientificName & coordinateUncertaintyInMeters & absence & individualCount & dropped & aphiaID & ⋯ & parvorder & parvorderid & organismID & category & organismQuantity & organismQuantityType & infrakingdomid & subkingdom & infrakingdom & subkingdomid\\\\\n",
       " <lgl> & <int> & <chr> & <chr> & <chr> & <chr> & <lgl> & <chr> & <lgl> & <int> & ⋯ & <chr> & <int> & <chr> & <chr> & <chr> & <chr> & <int> & <chr> & <chr> & <int>\\\\\n",
       "\\hline\n",
       "\t  TRUE & 2003 & https://www.ncei.noaa.gov/erddap/tabledap/deep\\_sea\\_corals.csv?ShallowFlag\\%2CDatasetID\\%2CCatalogNumber\\%2CSampleID\\%2CCitation\\%2CRepository\\%2CScientificName\\%2CVern {[}...{]} & urn:lsid:marinespecies.org:taxname:558    & Porifera         & 500m & FALSE & 1    & FALSE &    558 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t FALSE & 2000 & NA                                                                                                                                                                     & urn:lsid:marinespecies.org:taxname:124713 & Psolus squamatus & NA   & FALSE & NA   & FALSE & 124713 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t    NA & 1991 & http://calcofi.org/publications/calcofireports/v57/Vol57-Sakuma\\_pages.163-183.pdf                                                                                      & urn:lsid:marinespecies.org:taxname:1128   & Euphausiacea     & NA   & FALSE & 7150 & FALSE &   1128 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t    NA & 2007 & https://www.ncei.noaa.gov/erddap/tabledap/deep\\_sea\\_corals.csv?ShallowFlag\\%2CDatasetID\\%2CCatalogNumber\\%2CSampleID\\%2CCitation\\%2CRepository\\%2CScientificName\\%2CVern {[}...{]} & urn:lsid:marinespecies.org:taxname:125267 & Holaxonia        & 50m  & FALSE & 3    & FALSE & 125267 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t FALSE & 2002 & https://www.ncei.noaa.gov/erddap/tabledap/deep\\_sea\\_corals.csv?ShallowFlag\\%2CDatasetID\\%2CCatalogNumber\\%2CSampleID\\%2CCitation\\%2CRepository\\%2CScientificName\\%2CVern {[}...{]} & urn:lsid:marinespecies.org:taxname:132046 & Polymastia       & 500m & FALSE & 2    & FALSE & 132046 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t    NA & 2000 & NA                                                                                                                                                                     & urn:lsid:marinespecies.org:taxname:135392 & Nanomia          & NA   & FALSE & NA   & FALSE & 135392 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A tibble: 6 × 148\n",
       "\n",
       "| brackish &lt;lgl&gt; | date_year &lt;int&gt; | references &lt;chr&gt; | scientificNameID &lt;chr&gt; | scientificName &lt;chr&gt; | coordinateUncertaintyInMeters &lt;chr&gt; | absence &lt;lgl&gt; | individualCount &lt;chr&gt; | dropped &lt;lgl&gt; | aphiaID &lt;int&gt; | ⋯ ⋯ | parvorder &lt;chr&gt; | parvorderid &lt;int&gt; | organismID &lt;chr&gt; | category &lt;chr&gt; | organismQuantity &lt;chr&gt; | organismQuantityType &lt;chr&gt; | infrakingdomid &lt;int&gt; | subkingdom &lt;chr&gt; | infrakingdom &lt;chr&gt; | subkingdomid &lt;int&gt; |\n",
       "|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|\n",
       "|  TRUE | 2003 | https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...] | urn:lsid:marinespecies.org:taxname:558    | Porifera         | 500m | FALSE | 1    | FALSE |    558 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "| FALSE | 2000 | NA                                                                                                                                                                     | urn:lsid:marinespecies.org:taxname:124713 | Psolus squamatus | NA   | FALSE | NA   | FALSE | 124713 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "|    NA | 1991 | http://calcofi.org/publications/calcofireports/v57/Vol57-Sakuma_pages.163-183.pdf                                                                                      | urn:lsid:marinespecies.org:taxname:1128   | Euphausiacea     | NA   | FALSE | 7150 | FALSE |   1128 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "|    NA | 2007 | https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...] | urn:lsid:marinespecies.org:taxname:125267 | Holaxonia        | 50m  | FALSE | 3    | FALSE | 125267 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "| FALSE | 2002 | https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...] | urn:lsid:marinespecies.org:taxname:132046 | Polymastia       | 500m | FALSE | 2    | FALSE | 132046 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "|    NA | 2000 | NA                                                                                                                                                                     | urn:lsid:marinespecies.org:taxname:135392 | Nanomia          | NA   | FALSE | NA   | FALSE | 135392 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "\n"
      ],
      "text/plain": [
       "  brackish date_year\n",
       "1  TRUE    2003     \n",
       "2 FALSE    2000     \n",
       "3    NA    1991     \n",
       "4    NA    2007     \n",
       "5 FALSE    2002     \n",
       "6    NA    2000     \n",
       "  references                                                                                                                                                            \n",
       "1 https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]\n",
       "2 NA                                                                                                                                                                    \n",
       "3 http://calcofi.org/publications/calcofireports/v57/Vol57-Sakuma_pages.163-183.pdf                                                                                     \n",
       "4 https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]\n",
       "5 https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]\n",
       "6 NA                                                                                                                                                                    \n",
       "  scientificNameID                          scientificName  \n",
       "1 urn:lsid:marinespecies.org:taxname:558    Porifera        \n",
       "2 urn:lsid:marinespecies.org:taxname:124713 Psolus squamatus\n",
       "3 urn:lsid:marinespecies.org:taxname:1128   Euphausiacea    \n",
       "4 urn:lsid:marinespecies.org:taxname:125267 Holaxonia       \n",
       "5 urn:lsid:marinespecies.org:taxname:132046 Polymastia      \n",
       "6 urn:lsid:marinespecies.org:taxname:135392 Nanomia         \n",
       "  coordinateUncertaintyInMeters absence individualCount dropped aphiaID ⋯\n",
       "1 500m                          FALSE   1               FALSE      558  ⋯\n",
       "2 NA                            FALSE   NA              FALSE   124713  ⋯\n",
       "3 NA                            FALSE   7150            FALSE     1128  ⋯\n",
       "4 50m                           FALSE   3               FALSE   125267  ⋯\n",
       "5 500m                          FALSE   2               FALSE   132046  ⋯\n",
       "6 NA                            FALSE   NA              FALSE   135392  ⋯\n",
       "  parvorder parvorderid organismID category organismQuantity\n",
       "1 NA        NA          NA         NA       NA              \n",
       "2 NA        NA          NA         NA       NA              \n",
       "3 NA        NA          NA         NA       NA              \n",
       "4 NA        NA          NA         NA       NA              \n",
       "5 NA        NA          NA         NA       NA              \n",
       "6 NA        NA          NA         NA       NA              \n",
       "  organismQuantityType infrakingdomid subkingdom infrakingdom subkingdomid\n",
       "1 NA                   NA             NA         NA           NA          \n",
       "2 NA                   NA             NA         NA           NA          \n",
       "3 NA                   NA             NA         NA           NA          \n",
       "4 NA                   NA             NA         NA           NA          \n",
       "5 NA                   NA             NA         NA           NA          \n",
       "6 NA                   NA             NA         NA           NA          "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>6477</li><li>148</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 6477\n",
       "\\item 148\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 6477\n",
       "2. 148\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 6477  148"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A tibble: 6 × 137</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>brackish</th><th scope=col>date_year</th><th scope=col>references</th><th scope=col>scientificNameID</th><th scope=col>scientificName</th><th scope=col>coordinateUncertaintyInMeters</th><th scope=col>absence</th><th scope=col>individualCount</th><th scope=col>dropped</th><th scope=col>aphiaID</th><th scope=col>⋯</th><th scope=col>island</th><th scope=col>verbatimLocality</th><th scope=col>collectionID</th><th scope=col>islandGroup</th><th scope=col>typeStatus</th><th scope=col>georeferenceVerificationStatus</th><th scope=col>verbatimDepth</th><th scope=col>accessRights</th><th scope=col>subgenus</th><th scope=col>subgenusid</th></tr>\n",
       "\t<tr><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;int&gt;</th><th scope=col>⋯</th><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;lgl&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;int&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td> TRUE</td><td>2007</td><td>https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]</td><td>urn:lsid:marinespecies.org:taxname:558   </td><td>Porifera         </td><td>50m</td><td>FALSE</td><td>9 </td><td>FALSE</td><td>   558</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>FALSE</td><td>  NA</td><td>NA                                                                                                                                                                    </td><td>urn:lsid:marinespecies.org:taxname:274838</td><td>Sebastes proriger</td><td>NA </td><td>FALSE</td><td>NA</td><td>FALSE</td><td>274838</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>   NA</td><td>2000</td><td>NA                                                                                                                                                                    </td><td>urn:lsid:marinespecies.org:taxname:123081</td><td>Crinoidea        </td><td>NA </td><td>FALSE</td><td>NA</td><td>FALSE</td><td>123081</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td> TRUE</td><td>2007</td><td>https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]</td><td>urn:lsid:marinespecies.org:taxname:558   </td><td>Porifera         </td><td>50m</td><td>FALSE</td><td>2 </td><td>FALSE</td><td>   558</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>FALSE</td><td>2007</td><td>https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]</td><td>urn:lsid:marinespecies.org:taxname:132016</td><td>Poecillastra     </td><td>50m</td><td>FALSE</td><td>1 </td><td>FALSE</td><td>132016</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "\t<tr><td>   NA</td><td>2000</td><td>NA                                                                                                                                                                    </td><td>urn:lsid:marinespecies.org:taxname:123081</td><td>Crinoidea        </td><td>NA </td><td>FALSE</td><td>NA</td><td>FALSE</td><td>123081</td><td>⋯</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td><td>NA</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A tibble: 6 × 137\n",
       "\\begin{tabular}{lllllllllllllllllllll}\n",
       " brackish & date\\_year & references & scientificNameID & scientificName & coordinateUncertaintyInMeters & absence & individualCount & dropped & aphiaID & ⋯ & island & verbatimLocality & collectionID & islandGroup & typeStatus & georeferenceVerificationStatus & verbatimDepth & accessRights & subgenus & subgenusid\\\\\n",
       " <lgl> & <int> & <chr> & <chr> & <chr> & <chr> & <lgl> & <chr> & <lgl> & <int> & ⋯ & <lgl> & <chr> & <chr> & <lgl> & <chr> & <chr> & <chr> & <chr> & <chr> & <int>\\\\\n",
       "\\hline\n",
       "\t  TRUE & 2007 & https://www.ncei.noaa.gov/erddap/tabledap/deep\\_sea\\_corals.csv?ShallowFlag\\%2CDatasetID\\%2CCatalogNumber\\%2CSampleID\\%2CCitation\\%2CRepository\\%2CScientificName\\%2CVern {[}...{]} & urn:lsid:marinespecies.org:taxname:558    & Porifera          & 50m & FALSE & 9  & FALSE &    558 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t FALSE &   NA & NA                                                                                                                                                                     & urn:lsid:marinespecies.org:taxname:274838 & Sebastes proriger & NA  & FALSE & NA & FALSE & 274838 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t    NA & 2000 & NA                                                                                                                                                                     & urn:lsid:marinespecies.org:taxname:123081 & Crinoidea         & NA  & FALSE & NA & FALSE & 123081 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t  TRUE & 2007 & https://www.ncei.noaa.gov/erddap/tabledap/deep\\_sea\\_corals.csv?ShallowFlag\\%2CDatasetID\\%2CCatalogNumber\\%2CSampleID\\%2CCitation\\%2CRepository\\%2CScientificName\\%2CVern {[}...{]} & urn:lsid:marinespecies.org:taxname:558    & Porifera          & 50m & FALSE & 2  & FALSE &    558 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t FALSE & 2007 & https://www.ncei.noaa.gov/erddap/tabledap/deep\\_sea\\_corals.csv?ShallowFlag\\%2CDatasetID\\%2CCatalogNumber\\%2CSampleID\\%2CCitation\\%2CRepository\\%2CScientificName\\%2CVern {[}...{]} & urn:lsid:marinespecies.org:taxname:132016 & Poecillastra      & 50m & FALSE & 1  & FALSE & 132016 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\t    NA & 2000 & NA                                                                                                                                                                     & urn:lsid:marinespecies.org:taxname:123081 & Crinoidea         & NA  & FALSE & NA & FALSE & 123081 & ⋯ & NA & NA & NA & NA & NA & NA & NA & NA & NA & NA\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A tibble: 6 × 137\n",
       "\n",
       "| brackish &lt;lgl&gt; | date_year &lt;int&gt; | references &lt;chr&gt; | scientificNameID &lt;chr&gt; | scientificName &lt;chr&gt; | coordinateUncertaintyInMeters &lt;chr&gt; | absence &lt;lgl&gt; | individualCount &lt;chr&gt; | dropped &lt;lgl&gt; | aphiaID &lt;int&gt; | ⋯ ⋯ | island &lt;lgl&gt; | verbatimLocality &lt;chr&gt; | collectionID &lt;chr&gt; | islandGroup &lt;lgl&gt; | typeStatus &lt;chr&gt; | georeferenceVerificationStatus &lt;chr&gt; | verbatimDepth &lt;chr&gt; | accessRights &lt;chr&gt; | subgenus &lt;chr&gt; | subgenusid &lt;int&gt; |\n",
       "|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|\n",
       "|  TRUE | 2007 | https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...] | urn:lsid:marinespecies.org:taxname:558    | Porifera          | 50m | FALSE | 9  | FALSE |    558 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "| FALSE |   NA | NA                                                                                                                                                                     | urn:lsid:marinespecies.org:taxname:274838 | Sebastes proriger | NA  | FALSE | NA | FALSE | 274838 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "|    NA | 2000 | NA                                                                                                                                                                     | urn:lsid:marinespecies.org:taxname:123081 | Crinoidea         | NA  | FALSE | NA | FALSE | 123081 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "|  TRUE | 2007 | https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...] | urn:lsid:marinespecies.org:taxname:558    | Porifera          | 50m | FALSE | 2  | FALSE |    558 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "| FALSE | 2007 | https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...] | urn:lsid:marinespecies.org:taxname:132016 | Poecillastra      | 50m | FALSE | 1  | FALSE | 132016 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "|    NA | 2000 | NA                                                                                                                                                                     | urn:lsid:marinespecies.org:taxname:123081 | Crinoidea         | NA  | FALSE | NA | FALSE | 123081 | ⋯ | NA | NA | NA | NA | NA | NA | NA | NA | NA | NA |\n",
       "\n"
      ],
      "text/plain": [
       "  brackish date_year\n",
       "1  TRUE    2007     \n",
       "2 FALSE      NA     \n",
       "3    NA    2000     \n",
       "4  TRUE    2007     \n",
       "5 FALSE    2007     \n",
       "6    NA    2000     \n",
       "  references                                                                                                                                                            \n",
       "1 https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]\n",
       "2 NA                                                                                                                                                                    \n",
       "3 NA                                                                                                                                                                    \n",
       "4 https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]\n",
       "5 https://www.ncei.noaa.gov/erddap/tabledap/deep_sea_corals.csv?ShallowFlag%2CDatasetID%2CCatalogNumber%2CSampleID%2CCitation%2CRepository%2CScientificName%2CVern [...]\n",
       "6 NA                                                                                                                                                                    \n",
       "  scientificNameID                          scientificName   \n",
       "1 urn:lsid:marinespecies.org:taxname:558    Porifera         \n",
       "2 urn:lsid:marinespecies.org:taxname:274838 Sebastes proriger\n",
       "3 urn:lsid:marinespecies.org:taxname:123081 Crinoidea        \n",
       "4 urn:lsid:marinespecies.org:taxname:558    Porifera         \n",
       "5 urn:lsid:marinespecies.org:taxname:132016 Poecillastra     \n",
       "6 urn:lsid:marinespecies.org:taxname:123081 Crinoidea        \n",
       "  coordinateUncertaintyInMeters absence individualCount dropped aphiaID ⋯\n",
       "1 50m                           FALSE   9               FALSE      558  ⋯\n",
       "2 NA                            FALSE   NA              FALSE   274838  ⋯\n",
       "3 NA                            FALSE   NA              FALSE   123081  ⋯\n",
       "4 50m                           FALSE   2               FALSE      558  ⋯\n",
       "5 50m                           FALSE   1               FALSE   132016  ⋯\n",
       "6 NA                            FALSE   NA              FALSE   123081  ⋯\n",
       "  island verbatimLocality collectionID islandGroup typeStatus\n",
       "1 NA     NA               NA           NA          NA        \n",
       "2 NA     NA               NA           NA          NA        \n",
       "3 NA     NA               NA           NA          NA        \n",
       "4 NA     NA               NA           NA          NA        \n",
       "5 NA     NA               NA           NA          NA        \n",
       "6 NA     NA               NA           NA          NA        \n",
       "  georeferenceVerificationStatus verbatimDepth accessRights subgenus subgenusid\n",
       "1 NA                             NA            NA           NA       NA        \n",
       "2 NA                             NA            NA           NA       NA        \n",
       "3 NA                             NA            NA           NA       NA        \n",
       "4 NA                             NA            NA           NA       NA        \n",
       "5 NA                             NA            NA           NA       NA        \n",
       "6 NA                             NA            NA           NA       NA        "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>1600</li><li>137</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 1600\n",
       "\\item 137\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 1600\n",
       "2. 137\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 1600  137"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# take a quick look at both MPAs\n",
    "head(soquel_canyon_records)\n",
    "dim(soquel_canyon_records)\n",
    "\n",
    "head(point_lobos_records)\n",
    "dim(point_lobos_records)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Compute ES50"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Warning message in eval(expr, envir, enclos):\n",
      "“NAs introduced by coercion”\n",
      "Warning message in eval(expr, envir, enclos):\n",
      "“NAs introduced by coercion”\n"
     ]
    }
   ],
   "source": [
    "#convert data to allow ES50 calculation\n",
    "\n",
    "#convert individual counts from character to numeric\n",
    "#soquel_canyon_records\n",
    "soquel_canyon_records$individualCount <- as.numeric(soquel_canyon_records$individualCount)\n",
    "#convert NANs to 1\n",
    "soquel_canyon_records$individualCount[is.na(soquel_canyon_records$individualCount)] <- 1 \n",
    "soquel_canyon_records$Count <- 1 * soquel_canyon_records$individualCount\n",
    "\n",
    "#point_lobos_records\n",
    "point_lobos_records$individualCount <- as.numeric(point_lobos_records$individualCount)\n",
    "#convert NANs to 1\n",
    "point_lobos_records$individualCount[is.na(point_lobos_records$individualCount)] <- 1 \n",
    "point_lobos_records$Count <- 1 * point_lobos_records$individualCount\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 225 × 2</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>Category</th><th scope=col>x</th></tr>\n",
       "\t<tr><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td>Acanthoptilum              </td><td>268</td></tr>\n",
       "\t<tr><td>Actiniaria                 </td><td> 10</td></tr>\n",
       "\t<tr><td>Actinopterygii             </td><td> 11</td></tr>\n",
       "\t<tr><td>Aegina                     </td><td> 41</td></tr>\n",
       "\t<tr><td>Aglaophenia struthionides  </td><td>  3</td></tr>\n",
       "\t<tr><td>Agonidae                   </td><td>  4</td></tr>\n",
       "\t<tr><td>Alcidae                    </td><td>  4</td></tr>\n",
       "\t<tr><td>Alcyonacea                 </td><td> 20</td></tr>\n",
       "\t<tr><td>Anarrhichthys ocellatus    </td><td>  2</td></tr>\n",
       "\t<tr><td>Animalia                   </td><td>  1</td></tr>\n",
       "\t<tr><td>Anthomastus                </td><td>154</td></tr>\n",
       "\t<tr><td>Antipatharia               </td><td>  1</td></tr>\n",
       "\t<tr><td>Appendicularia             </td><td> 26</td></tr>\n",
       "\t<tr><td>Argentina sialis           </td><td>104</td></tr>\n",
       "\t<tr><td>Asbestopluma (Asbestopluma)</td><td> 11</td></tr>\n",
       "\t<tr><td>Asteroidea                 </td><td>  1</td></tr>\n",
       "\t<tr><td>Astropecten californicus   </td><td> 32</td></tr>\n",
       "\t<tr><td>Aurelia                    </td><td> 46</td></tr>\n",
       "\t<tr><td>Balaenoptera musculus      </td><td> 25</td></tr>\n",
       "\t<tr><td>Balanophyllia              </td><td>  1</td></tr>\n",
       "\t<tr><td>Bathochordaeus             </td><td>  4</td></tr>\n",
       "\t<tr><td>Bathylagidae               </td><td>  4</td></tr>\n",
       "\t<tr><td>Beroe                      </td><td>127</td></tr>\n",
       "\t<tr><td>Bolinopsis                 </td><td> 69</td></tr>\n",
       "\t<tr><td>Caniformia                 </td><td>  1</td></tr>\n",
       "\t<tr><td>Caridea                    </td><td>935</td></tr>\n",
       "\t<tr><td>Caryophylliidae            </td><td> 41</td></tr>\n",
       "\t<tr><td>Cerorhinca monocerata      </td><td> 45</td></tr>\n",
       "\t<tr><td>Cetacea                    </td><td>  3</td></tr>\n",
       "\t<tr><td>Chaetognatha               </td><td>  5</td></tr>\n",
       "\t<tr><td>⋮</td><td>⋮</td></tr>\n",
       "\t<tr><td>Sebastes pinniger     </td><td>  3</td></tr>\n",
       "\t<tr><td>Sebastes rufus        </td><td> 14</td></tr>\n",
       "\t<tr><td>Sebastes saxicola     </td><td>244</td></tr>\n",
       "\t<tr><td>Sebastes semicinctus  </td><td> 78</td></tr>\n",
       "\t<tr><td>Sebastes wilsoni      </td><td> 31</td></tr>\n",
       "\t<tr><td>Sergestidae           </td><td>477</td></tr>\n",
       "\t<tr><td>Siphonophorae         </td><td>  5</td></tr>\n",
       "\t<tr><td>Solmissus             </td><td>  7</td></tr>\n",
       "\t<tr><td>Stercorarius pomarinus</td><td>  2</td></tr>\n",
       "\t<tr><td>Sterna                </td><td>  1</td></tr>\n",
       "\t<tr><td>Sterna paradisaea     </td><td>  2</td></tr>\n",
       "\t<tr><td>Stylaster             </td><td>  1</td></tr>\n",
       "\t<tr><td>Stylaster californicus</td><td>  1</td></tr>\n",
       "\t<tr><td>Swiftia               </td><td>  4</td></tr>\n",
       "\t<tr><td>Swiftia pacifica      </td><td>  1</td></tr>\n",
       "\t<tr><td>Swiftia simplex       </td><td>  4</td></tr>\n",
       "\t<tr><td>Synodus lucioceps     </td><td>  4</td></tr>\n",
       "\t<tr><td>Tetronarce californica</td><td> 24</td></tr>\n",
       "\t<tr><td>Teuthida              </td><td>542</td></tr>\n",
       "\t<tr><td>Thenea muricata       </td><td>  5</td></tr>\n",
       "\t<tr><td>Thetys vagina         </td><td>147</td></tr>\n",
       "\t<tr><td>Trachipterus altivelis</td><td>  7</td></tr>\n",
       "\t<tr><td>Trachurus symmetricus </td><td>  1</td></tr>\n",
       "\t<tr><td>Tunicata              </td><td>  1</td></tr>\n",
       "\t<tr><td>Uria aalge            </td><td>158</td></tr>\n",
       "\t<tr><td>Virgulariidae         </td><td> 28</td></tr>\n",
       "\t<tr><td>Zalophus californianus</td><td> 14</td></tr>\n",
       "\t<tr><td>Zaniolepis            </td><td>139</td></tr>\n",
       "\t<tr><td>Zaprora silenus       </td><td>  1</td></tr>\n",
       "\t<tr><td>Zoarcidae             </td><td>  4</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 225 × 2\n",
       "\\begin{tabular}{ll}\n",
       " Category & x\\\\\n",
       " <chr> & <dbl>\\\\\n",
       "\\hline\n",
       "\t Acanthoptilum               & 268\\\\\n",
       "\t Actiniaria                  &  10\\\\\n",
       "\t Actinopterygii              &  11\\\\\n",
       "\t Aegina                      &  41\\\\\n",
       "\t Aglaophenia struthionides   &   3\\\\\n",
       "\t Agonidae                    &   4\\\\\n",
       "\t Alcidae                     &   4\\\\\n",
       "\t Alcyonacea                  &  20\\\\\n",
       "\t Anarrhichthys ocellatus     &   2\\\\\n",
       "\t Animalia                    &   1\\\\\n",
       "\t Anthomastus                 & 154\\\\\n",
       "\t Antipatharia                &   1\\\\\n",
       "\t Appendicularia              &  26\\\\\n",
       "\t Argentina sialis            & 104\\\\\n",
       "\t Asbestopluma (Asbestopluma) &  11\\\\\n",
       "\t Asteroidea                  &   1\\\\\n",
       "\t Astropecten californicus    &  32\\\\\n",
       "\t Aurelia                     &  46\\\\\n",
       "\t Balaenoptera musculus       &  25\\\\\n",
       "\t Balanophyllia               &   1\\\\\n",
       "\t Bathochordaeus              &   4\\\\\n",
       "\t Bathylagidae                &   4\\\\\n",
       "\t Beroe                       & 127\\\\\n",
       "\t Bolinopsis                  &  69\\\\\n",
       "\t Caniformia                  &   1\\\\\n",
       "\t Caridea                     & 935\\\\\n",
       "\t Caryophylliidae             &  41\\\\\n",
       "\t Cerorhinca monocerata       &  45\\\\\n",
       "\t Cetacea                     &   3\\\\\n",
       "\t Chaetognatha                &   5\\\\\n",
       "\t ⋮ & ⋮\\\\\n",
       "\t Sebastes pinniger      &   3\\\\\n",
       "\t Sebastes rufus         &  14\\\\\n",
       "\t Sebastes saxicola      & 244\\\\\n",
       "\t Sebastes semicinctus   &  78\\\\\n",
       "\t Sebastes wilsoni       &  31\\\\\n",
       "\t Sergestidae            & 477\\\\\n",
       "\t Siphonophorae          &   5\\\\\n",
       "\t Solmissus              &   7\\\\\n",
       "\t Stercorarius pomarinus &   2\\\\\n",
       "\t Sterna                 &   1\\\\\n",
       "\t Sterna paradisaea      &   2\\\\\n",
       "\t Stylaster              &   1\\\\\n",
       "\t Stylaster californicus &   1\\\\\n",
       "\t Swiftia                &   4\\\\\n",
       "\t Swiftia pacifica       &   1\\\\\n",
       "\t Swiftia simplex        &   4\\\\\n",
       "\t Synodus lucioceps      &   4\\\\\n",
       "\t Tetronarce californica &  24\\\\\n",
       "\t Teuthida               & 542\\\\\n",
       "\t Thenea muricata        &   5\\\\\n",
       "\t Thetys vagina          & 147\\\\\n",
       "\t Trachipterus altivelis &   7\\\\\n",
       "\t Trachurus symmetricus  &   1\\\\\n",
       "\t Tunicata               &   1\\\\\n",
       "\t Uria aalge             & 158\\\\\n",
       "\t Virgulariidae          &  28\\\\\n",
       "\t Zalophus californianus &  14\\\\\n",
       "\t Zaniolepis             & 139\\\\\n",
       "\t Zaprora silenus        &   1\\\\\n",
       "\t Zoarcidae              &   4\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 225 × 2\n",
       "\n",
       "| Category &lt;chr&gt; | x &lt;dbl&gt; |\n",
       "|---|---|\n",
       "| Acanthoptilum               | 268 |\n",
       "| Actiniaria                  |  10 |\n",
       "| Actinopterygii              |  11 |\n",
       "| Aegina                      |  41 |\n",
       "| Aglaophenia struthionides   |   3 |\n",
       "| Agonidae                    |   4 |\n",
       "| Alcidae                     |   4 |\n",
       "| Alcyonacea                  |  20 |\n",
       "| Anarrhichthys ocellatus     |   2 |\n",
       "| Animalia                    |   1 |\n",
       "| Anthomastus                 | 154 |\n",
       "| Antipatharia                |   1 |\n",
       "| Appendicularia              |  26 |\n",
       "| Argentina sialis            | 104 |\n",
       "| Asbestopluma (Asbestopluma) |  11 |\n",
       "| Asteroidea                  |   1 |\n",
       "| Astropecten californicus    |  32 |\n",
       "| Aurelia                     |  46 |\n",
       "| Balaenoptera musculus       |  25 |\n",
       "| Balanophyllia               |   1 |\n",
       "| Bathochordaeus              |   4 |\n",
       "| Bathylagidae                |   4 |\n",
       "| Beroe                       | 127 |\n",
       "| Bolinopsis                  |  69 |\n",
       "| Caniformia                  |   1 |\n",
       "| Caridea                     | 935 |\n",
       "| Caryophylliidae             |  41 |\n",
       "| Cerorhinca monocerata       |  45 |\n",
       "| Cetacea                     |   3 |\n",
       "| Chaetognatha                |   5 |\n",
       "| ⋮ | ⋮ |\n",
       "| Sebastes pinniger      |   3 |\n",
       "| Sebastes rufus         |  14 |\n",
       "| Sebastes saxicola      | 244 |\n",
       "| Sebastes semicinctus   |  78 |\n",
       "| Sebastes wilsoni       |  31 |\n",
       "| Sergestidae            | 477 |\n",
       "| Siphonophorae          |   5 |\n",
       "| Solmissus              |   7 |\n",
       "| Stercorarius pomarinus |   2 |\n",
       "| Sterna                 |   1 |\n",
       "| Sterna paradisaea      |   2 |\n",
       "| Stylaster              |   1 |\n",
       "| Stylaster californicus |   1 |\n",
       "| Swiftia                |   4 |\n",
       "| Swiftia pacifica       |   1 |\n",
       "| Swiftia simplex        |   4 |\n",
       "| Synodus lucioceps      |   4 |\n",
       "| Tetronarce californica |  24 |\n",
       "| Teuthida               | 542 |\n",
       "| Thenea muricata        |   5 |\n",
       "| Thetys vagina          | 147 |\n",
       "| Trachipterus altivelis |   7 |\n",
       "| Trachurus symmetricus  |   1 |\n",
       "| Tunicata               |   1 |\n",
       "| Uria aalge             | 158 |\n",
       "| Virgulariidae          |  28 |\n",
       "| Zalophus californianus |  14 |\n",
       "| Zaniolepis             | 139 |\n",
       "| Zaprora silenus        |   1 |\n",
       "| Zoarcidae              |   4 |\n",
       "\n"
      ],
      "text/plain": [
       "    Category                    x  \n",
       "1   Acanthoptilum               268\n",
       "2   Actiniaria                   10\n",
       "3   Actinopterygii               11\n",
       "4   Aegina                       41\n",
       "5   Aglaophenia struthionides     3\n",
       "6   Agonidae                      4\n",
       "7   Alcidae                       4\n",
       "8   Alcyonacea                   20\n",
       "9   Anarrhichthys ocellatus       2\n",
       "10  Animalia                      1\n",
       "11  Anthomastus                 154\n",
       "12  Antipatharia                  1\n",
       "13  Appendicularia               26\n",
       "14  Argentina sialis            104\n",
       "15  Asbestopluma (Asbestopluma)  11\n",
       "16  Asteroidea                    1\n",
       "17  Astropecten californicus     32\n",
       "18  Aurelia                      46\n",
       "19  Balaenoptera musculus        25\n",
       "20  Balanophyllia                 1\n",
       "21  Bathochordaeus                4\n",
       "22  Bathylagidae                  4\n",
       "23  Beroe                       127\n",
       "24  Bolinopsis                   69\n",
       "25  Caniformia                    1\n",
       "26  Caridea                     935\n",
       "27  Caryophylliidae              41\n",
       "28  Cerorhinca monocerata        45\n",
       "29  Cetacea                       3\n",
       "30  Chaetognatha                  5\n",
       "⋮   ⋮                           ⋮  \n",
       "196 Sebastes pinniger             3\n",
       "197 Sebastes rufus               14\n",
       "198 Sebastes saxicola           244\n",
       "199 Sebastes semicinctus         78\n",
       "200 Sebastes wilsoni             31\n",
       "201 Sergestidae                 477\n",
       "202 Siphonophorae                 5\n",
       "203 Solmissus                     7\n",
       "204 Stercorarius pomarinus        2\n",
       "205 Sterna                        1\n",
       "206 Sterna paradisaea             2\n",
       "207 Stylaster                     1\n",
       "208 Stylaster californicus        1\n",
       "209 Swiftia                       4\n",
       "210 Swiftia pacifica              1\n",
       "211 Swiftia simplex               4\n",
       "212 Synodus lucioceps             4\n",
       "213 Tetronarce californica       24\n",
       "214 Teuthida                    542\n",
       "215 Thenea muricata               5\n",
       "216 Thetys vagina               147\n",
       "217 Trachipterus altivelis        7\n",
       "218 Trachurus symmetricus         1\n",
       "219 Tunicata                      1\n",
       "220 Uria aalge                  158\n",
       "221 Virgulariidae                28\n",
       "222 Zalophus californianus       14\n",
       "223 Zaniolepis                  139\n",
       "224 Zaprora silenus               1\n",
       "225 Zoarcidae                     4"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 69 × 2</caption>\n",
       "<thead>\n",
       "\t<tr><th scope=col>Category</th><th scope=col>x</th></tr>\n",
       "\t<tr><th scope=col>&lt;chr&gt;</th><th scope=col>&lt;dbl&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><td>Actiniaria                   </td><td>  5</td></tr>\n",
       "\t<tr><td>Actinopterygii               </td><td>  2</td></tr>\n",
       "\t<tr><td>Aegina                       </td><td>  1</td></tr>\n",
       "\t<tr><td>Alcyonacea                   </td><td>  4</td></tr>\n",
       "\t<tr><td>Amphiophiura superba         </td><td>  5</td></tr>\n",
       "\t<tr><td>Amphiura (Amphiura) diomedeae</td><td>  2</td></tr>\n",
       "\t<tr><td>Anelassorhynchus abyssalis   </td><td>  2</td></tr>\n",
       "\t<tr><td>Asteroidea                   </td><td>  4</td></tr>\n",
       "\t<tr><td>Asteronyx loveni             </td><td>  1</td></tr>\n",
       "\t<tr><td>Balaenoptera musculus        </td><td> 22</td></tr>\n",
       "\t<tr><td>Bathyraja interrupta         </td><td>  1</td></tr>\n",
       "\t<tr><td>Bothrocara molle             </td><td>  1</td></tr>\n",
       "\t<tr><td>Calappidae                   </td><td>  1</td></tr>\n",
       "\t<tr><td>Chimaera                     </td><td>  1</td></tr>\n",
       "\t<tr><td>Chordata                     </td><td> 12</td></tr>\n",
       "\t<tr><td>Cnidaria                     </td><td>  3</td></tr>\n",
       "\t<tr><td>Crinoidea                    </td><td> 77</td></tr>\n",
       "\t<tr><td>Decapoda                     </td><td>  4</td></tr>\n",
       "\t<tr><td>Delphinidae                  </td><td>  7</td></tr>\n",
       "\t<tr><td>Desmophyllum pertusum        </td><td>  8</td></tr>\n",
       "\t<tr><td>Diphyes dispar               </td><td>  5</td></tr>\n",
       "\t<tr><td>Gorgonocephalidae            </td><td>  1</td></tr>\n",
       "\t<tr><td>Hedingia californica         </td><td>  1</td></tr>\n",
       "\t<tr><td>Heterochone calyx            </td><td>  1</td></tr>\n",
       "\t<tr><td>Holaxonia                    </td><td>105</td></tr>\n",
       "\t<tr><td>Hydrolagus colliei           </td><td>  3</td></tr>\n",
       "\t<tr><td>Icelinus filamentosus        </td><td>  1</td></tr>\n",
       "\t<tr><td>Labridae                     </td><td>  1</td></tr>\n",
       "\t<tr><td>Larus californicus           </td><td>  4</td></tr>\n",
       "\t<tr><td>Larus occidentalis           </td><td>  5</td></tr>\n",
       "\t<tr><td>⋮</td><td>⋮</td></tr>\n",
       "\t<tr><td>Parapagurus benedicti       </td><td>   1</td></tr>\n",
       "\t<tr><td>Pelecanus occidentalis      </td><td>   2</td></tr>\n",
       "\t<tr><td>Pennatulacea                </td><td> 143</td></tr>\n",
       "\t<tr><td>Pennatulidae                </td><td>   5</td></tr>\n",
       "\t<tr><td>Phalacrocorax penicillatus  </td><td>   2</td></tr>\n",
       "\t<tr><td>Phalaropus lobatus          </td><td>  19</td></tr>\n",
       "\t<tr><td>Phocoenoides dalli          </td><td>   4</td></tr>\n",
       "\t<tr><td>Pleuronectiformes           </td><td>  11</td></tr>\n",
       "\t<tr><td>Plexauridae                 </td><td>   9</td></tr>\n",
       "\t<tr><td>Poecillastra                </td><td>   2</td></tr>\n",
       "\t<tr><td>Porifera                    </td><td>1310</td></tr>\n",
       "\t<tr><td>Ptilosarcus                 </td><td>   3</td></tr>\n",
       "\t<tr><td>Puffinus griseus            </td><td>   2</td></tr>\n",
       "\t<tr><td>Rathbunaster californicus   </td><td>  11</td></tr>\n",
       "\t<tr><td>Rhabdocalyptus dawsoni      </td><td>  53</td></tr>\n",
       "\t<tr><td>Rhinobatidae                </td><td>   1</td></tr>\n",
       "\t<tr><td>Scleractinia                </td><td>  31</td></tr>\n",
       "\t<tr><td>Scorpaenidae                </td><td>   5</td></tr>\n",
       "\t<tr><td>Sebastes                    </td><td>  59</td></tr>\n",
       "\t<tr><td>Sebastes chlorostictus      </td><td>   1</td></tr>\n",
       "\t<tr><td>Sebastes proriger           </td><td>   1</td></tr>\n",
       "\t<tr><td>Sebastes ruberrimus         </td><td>   2</td></tr>\n",
       "\t<tr><td>Sebastes saxicola           </td><td>   1</td></tr>\n",
       "\t<tr><td>Sebastolobus                </td><td>   1</td></tr>\n",
       "\t<tr><td>Staurocalyptus              </td><td>   1</td></tr>\n",
       "\t<tr><td>Stylaster californicus      </td><td>   9</td></tr>\n",
       "\t<tr><td>Swiftia                     </td><td>  15</td></tr>\n",
       "\t<tr><td>Thrissacanthias penicillatus</td><td>   4</td></tr>\n",
       "\t<tr><td>Uria aalge                  </td><td>   1</td></tr>\n",
       "\t<tr><td>Xestospongia                </td><td>   6</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 69 × 2\n",
       "\\begin{tabular}{ll}\n",
       " Category & x\\\\\n",
       " <chr> & <dbl>\\\\\n",
       "\\hline\n",
       "\t Actiniaria                    &   5\\\\\n",
       "\t Actinopterygii                &   2\\\\\n",
       "\t Aegina                        &   1\\\\\n",
       "\t Alcyonacea                    &   4\\\\\n",
       "\t Amphiophiura superba          &   5\\\\\n",
       "\t Amphiura (Amphiura) diomedeae &   2\\\\\n",
       "\t Anelassorhynchus abyssalis    &   2\\\\\n",
       "\t Asteroidea                    &   4\\\\\n",
       "\t Asteronyx loveni              &   1\\\\\n",
       "\t Balaenoptera musculus         &  22\\\\\n",
       "\t Bathyraja interrupta          &   1\\\\\n",
       "\t Bothrocara molle              &   1\\\\\n",
       "\t Calappidae                    &   1\\\\\n",
       "\t Chimaera                      &   1\\\\\n",
       "\t Chordata                      &  12\\\\\n",
       "\t Cnidaria                      &   3\\\\\n",
       "\t Crinoidea                     &  77\\\\\n",
       "\t Decapoda                      &   4\\\\\n",
       "\t Delphinidae                   &   7\\\\\n",
       "\t Desmophyllum pertusum         &   8\\\\\n",
       "\t Diphyes dispar                &   5\\\\\n",
       "\t Gorgonocephalidae             &   1\\\\\n",
       "\t Hedingia californica          &   1\\\\\n",
       "\t Heterochone calyx             &   1\\\\\n",
       "\t Holaxonia                     & 105\\\\\n",
       "\t Hydrolagus colliei            &   3\\\\\n",
       "\t Icelinus filamentosus         &   1\\\\\n",
       "\t Labridae                      &   1\\\\\n",
       "\t Larus californicus            &   4\\\\\n",
       "\t Larus occidentalis            &   5\\\\\n",
       "\t ⋮ & ⋮\\\\\n",
       "\t Parapagurus benedicti        &    1\\\\\n",
       "\t Pelecanus occidentalis       &    2\\\\\n",
       "\t Pennatulacea                 &  143\\\\\n",
       "\t Pennatulidae                 &    5\\\\\n",
       "\t Phalacrocorax penicillatus   &    2\\\\\n",
       "\t Phalaropus lobatus           &   19\\\\\n",
       "\t Phocoenoides dalli           &    4\\\\\n",
       "\t Pleuronectiformes            &   11\\\\\n",
       "\t Plexauridae                  &    9\\\\\n",
       "\t Poecillastra                 &    2\\\\\n",
       "\t Porifera                     & 1310\\\\\n",
       "\t Ptilosarcus                  &    3\\\\\n",
       "\t Puffinus griseus             &    2\\\\\n",
       "\t Rathbunaster californicus    &   11\\\\\n",
       "\t Rhabdocalyptus dawsoni       &   53\\\\\n",
       "\t Rhinobatidae                 &    1\\\\\n",
       "\t Scleractinia                 &   31\\\\\n",
       "\t Scorpaenidae                 &    5\\\\\n",
       "\t Sebastes                     &   59\\\\\n",
       "\t Sebastes chlorostictus       &    1\\\\\n",
       "\t Sebastes proriger            &    1\\\\\n",
       "\t Sebastes ruberrimus          &    2\\\\\n",
       "\t Sebastes saxicola            &    1\\\\\n",
       "\t Sebastolobus                 &    1\\\\\n",
       "\t Staurocalyptus               &    1\\\\\n",
       "\t Stylaster californicus       &    9\\\\\n",
       "\t Swiftia                      &   15\\\\\n",
       "\t Thrissacanthias penicillatus &    4\\\\\n",
       "\t Uria aalge                   &    1\\\\\n",
       "\t Xestospongia                 &    6\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 69 × 2\n",
       "\n",
       "| Category &lt;chr&gt; | x &lt;dbl&gt; |\n",
       "|---|---|\n",
       "| Actiniaria                    |   5 |\n",
       "| Actinopterygii                |   2 |\n",
       "| Aegina                        |   1 |\n",
       "| Alcyonacea                    |   4 |\n",
       "| Amphiophiura superba          |   5 |\n",
       "| Amphiura (Amphiura) diomedeae |   2 |\n",
       "| Anelassorhynchus abyssalis    |   2 |\n",
       "| Asteroidea                    |   4 |\n",
       "| Asteronyx loveni              |   1 |\n",
       "| Balaenoptera musculus         |  22 |\n",
       "| Bathyraja interrupta          |   1 |\n",
       "| Bothrocara molle              |   1 |\n",
       "| Calappidae                    |   1 |\n",
       "| Chimaera                      |   1 |\n",
       "| Chordata                      |  12 |\n",
       "| Cnidaria                      |   3 |\n",
       "| Crinoidea                     |  77 |\n",
       "| Decapoda                      |   4 |\n",
       "| Delphinidae                   |   7 |\n",
       "| Desmophyllum pertusum         |   8 |\n",
       "| Diphyes dispar                |   5 |\n",
       "| Gorgonocephalidae             |   1 |\n",
       "| Hedingia californica          |   1 |\n",
       "| Heterochone calyx             |   1 |\n",
       "| Holaxonia                     | 105 |\n",
       "| Hydrolagus colliei            |   3 |\n",
       "| Icelinus filamentosus         |   1 |\n",
       "| Labridae                      |   1 |\n",
       "| Larus californicus            |   4 |\n",
       "| Larus occidentalis            |   5 |\n",
       "| ⋮ | ⋮ |\n",
       "| Parapagurus benedicti        |    1 |\n",
       "| Pelecanus occidentalis       |    2 |\n",
       "| Pennatulacea                 |  143 |\n",
       "| Pennatulidae                 |    5 |\n",
       "| Phalacrocorax penicillatus   |    2 |\n",
       "| Phalaropus lobatus           |   19 |\n",
       "| Phocoenoides dalli           |    4 |\n",
       "| Pleuronectiformes            |   11 |\n",
       "| Plexauridae                  |    9 |\n",
       "| Poecillastra                 |    2 |\n",
       "| Porifera                     | 1310 |\n",
       "| Ptilosarcus                  |    3 |\n",
       "| Puffinus griseus             |    2 |\n",
       "| Rathbunaster californicus    |   11 |\n",
       "| Rhabdocalyptus dawsoni       |   53 |\n",
       "| Rhinobatidae                 |    1 |\n",
       "| Scleractinia                 |   31 |\n",
       "| Scorpaenidae                 |    5 |\n",
       "| Sebastes                     |   59 |\n",
       "| Sebastes chlorostictus       |    1 |\n",
       "| Sebastes proriger            |    1 |\n",
       "| Sebastes ruberrimus          |    2 |\n",
       "| Sebastes saxicola            |    1 |\n",
       "| Sebastolobus                 |    1 |\n",
       "| Staurocalyptus               |    1 |\n",
       "| Stylaster californicus       |    9 |\n",
       "| Swiftia                      |   15 |\n",
       "| Thrissacanthias penicillatus |    4 |\n",
       "| Uria aalge                   |    1 |\n",
       "| Xestospongia                 |    6 |\n",
       "\n"
      ],
      "text/plain": [
       "   Category                      x   \n",
       "1  Actiniaria                      5 \n",
       "2  Actinopterygii                  2 \n",
       "3  Aegina                          1 \n",
       "4  Alcyonacea                      4 \n",
       "5  Amphiophiura superba            5 \n",
       "6  Amphiura (Amphiura) diomedeae   2 \n",
       "7  Anelassorhynchus abyssalis      2 \n",
       "8  Asteroidea                      4 \n",
       "9  Asteronyx loveni                1 \n",
       "10 Balaenoptera musculus          22 \n",
       "11 Bathyraja interrupta            1 \n",
       "12 Bothrocara molle                1 \n",
       "13 Calappidae                      1 \n",
       "14 Chimaera                        1 \n",
       "15 Chordata                       12 \n",
       "16 Cnidaria                        3 \n",
       "17 Crinoidea                      77 \n",
       "18 Decapoda                        4 \n",
       "19 Delphinidae                     7 \n",
       "20 Desmophyllum pertusum           8 \n",
       "21 Diphyes dispar                  5 \n",
       "22 Gorgonocephalidae               1 \n",
       "23 Hedingia californica            1 \n",
       "24 Heterochone calyx               1 \n",
       "25 Holaxonia                     105 \n",
       "26 Hydrolagus colliei              3 \n",
       "27 Icelinus filamentosus           1 \n",
       "28 Labridae                        1 \n",
       "29 Larus californicus              4 \n",
       "30 Larus occidentalis              5 \n",
       "⋮  ⋮                             ⋮   \n",
       "40 Parapagurus benedicti            1\n",
       "41 Pelecanus occidentalis           2\n",
       "42 Pennatulacea                   143\n",
       "43 Pennatulidae                     5\n",
       "44 Phalacrocorax penicillatus       2\n",
       "45 Phalaropus lobatus              19\n",
       "46 Phocoenoides dalli               4\n",
       "47 Pleuronectiformes               11\n",
       "48 Plexauridae                      9\n",
       "49 Poecillastra                     2\n",
       "50 Porifera                      1310\n",
       "51 Ptilosarcus                      3\n",
       "52 Puffinus griseus                 2\n",
       "53 Rathbunaster californicus       11\n",
       "54 Rhabdocalyptus dawsoni          53\n",
       "55 Rhinobatidae                     1\n",
       "56 Scleractinia                    31\n",
       "57 Scorpaenidae                     5\n",
       "58 Sebastes                        59\n",
       "59 Sebastes chlorostictus           1\n",
       "60 Sebastes proriger                1\n",
       "61 Sebastes ruberrimus              2\n",
       "62 Sebastes saxicola                1\n",
       "63 Sebastolobus                     1\n",
       "64 Staurocalyptus                   1\n",
       "65 Stylaster californicus           9\n",
       "66 Swiftia                         15\n",
       "67 Thrissacanthias penicillatus     4\n",
       "68 Uria aalge                       1\n",
       "69 Xestospongia                     6"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#calculate the number of unique species in each MPA\n",
    "#soquel_canyon_records\n",
    "SpeciesCount_SoquelCanyon <- aggregate(soquel_canyon_records$Count, \n",
    "                                       by=list(Category=soquel_canyon_records$scientificName),\n",
    "                                       FUN=sum)\n",
    "\n",
    "#point_lobos_records\n",
    "SpeciesCount_PointLobos <- aggregate(point_lobos_records$Count, \n",
    "                                     by=list(Category=point_lobos_records$scientificName),\n",
    "                                     FUN=sum)\n",
    "#visualize species counts\n",
    "SpeciesCount_SoquelCanyon\n",
    "SpeciesCount_PointLobos\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "2.87814552026609"
      ],
      "text/latex": [
       "2.87814552026609"
      ],
      "text/markdown": [
       "2.87814552026609"
      ],
      "text/plain": [
       "[1] 2.878146\n",
       "attr(,\"Subsample\")\n",
       "[1] 50"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "12.2153017829868"
      ],
      "text/latex": [
       "12.2153017829868"
      ],
      "text/markdown": [
       "12.2153017829868"
      ],
      "text/plain": [
       "[1] 12.2153\n",
       "attr(,\"Subsample\")\n",
       "[1] 50"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#calculate the ES50 value\n",
    "#Hurlbert's index, or ES50, is the expected number of distinct species \n",
    "#in a random sample of 50 observations; it is a bias-indepenent indicator \n",
    "#on marine biodiversity richness. It has the ability to somewhat correct for sampling bias.\n",
    "#0-20: low, 30-50: high\n",
    "\n",
    "ES50_SoquelCanyon = rarefy(SpeciesCount_SoquelCanyon$x,50) \n",
    "ES50_PointLobos = rarefy(SpeciesCount_PointLobos$x,50) \n",
    "\n",
    "ES50_SoquelCanyon\n",
    "ES50_PointLobos"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Retrieve taxonomic group  & Environmental Data of interest"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We are going to focus on the Soquel Canyon MPA and the genus Sebastes (rockfish)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "#input: polygon, taxa, metric\n",
    "#return: df with the metric value, lat, and long of each sighting within polygon\n",
    "\n",
    "EnvironmentalMetricPresent <- function(polygon, taxa, metric) {\n",
    "    \n",
    "    #obtain occurrences of taxa of interest within the polygon of interest\n",
    "    taxaPresence <- occurrence(geometry = polygon, taxa)\n",
    "    \n",
    "    #create dataframe of the environmental data as well as the lat/lon of the sighting\n",
    "    environmentalMetricDF <- data.frame(EnvironmentalData = taxaPresence[metric], \n",
    "                                   longitude = taxaPresence$decimalLongitude, \n",
    "                                   latitude = taxaPresence$decimalLatitude, \n",
    "                                   date = taxaPresence$eventDate)\n",
    "    \n",
    "    return(environmentalMetricDF)\n",
    "}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Retrieved 183 records of approximately 183 (100%)\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<table>\n",
       "<caption>A data.frame: 6 × 4</caption>\n",
       "<thead>\n",
       "\t<tr><th></th><th scope=col>sst</th><th scope=col>longitude</th><th scope=col>latitude</th><th scope=col>date</th></tr>\n",
       "\t<tr><th></th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;dbl&gt;</th><th scope=col>&lt;fct&gt;</th></tr>\n",
       "</thead>\n",
       "<tbody>\n",
       "\t<tr><th scope=row>1</th><td>12.78</td><td>-121.9688</td><td>36.83320</td><td>2009-05-21T05:12:46Z</td></tr>\n",
       "\t<tr><th scope=row>2</th><td>12.78</td><td>-121.9765</td><td>36.83096</td><td>2000-02-10T21:24:00Z</td></tr>\n",
       "\t<tr><th scope=row>3</th><td>12.75</td><td>-121.9790</td><td>36.84930</td><td>2011-05-21T10:06:00Z</td></tr>\n",
       "\t<tr><th scope=row>4</th><td>12.78</td><td>-121.9467</td><td>36.81833</td><td>NA                  </td></tr>\n",
       "\t<tr><th scope=row>5</th><td>12.75</td><td>-121.9817</td><td>36.84500</td><td>2000-05-12T04:02:00Z</td></tr>\n",
       "\t<tr><th scope=row>6</th><td>12.75</td><td>-121.9853</td><td>36.84600</td><td>2013-05-11T04:04:53Z</td></tr>\n",
       "</tbody>\n",
       "</table>\n"
      ],
      "text/latex": [
       "A data.frame: 6 × 4\n",
       "\\begin{tabular}{r|llll}\n",
       "  & sst & longitude & latitude & date\\\\\n",
       "  & <dbl> & <dbl> & <dbl> & <fct>\\\\\n",
       "\\hline\n",
       "\t1 & 12.78 & -121.9688 & 36.83320 & 2009-05-21T05:12:46Z\\\\\n",
       "\t2 & 12.78 & -121.9765 & 36.83096 & 2000-02-10T21:24:00Z\\\\\n",
       "\t3 & 12.75 & -121.9790 & 36.84930 & 2011-05-21T10:06:00Z\\\\\n",
       "\t4 & 12.78 & -121.9467 & 36.81833 & NA                  \\\\\n",
       "\t5 & 12.75 & -121.9817 & 36.84500 & 2000-05-12T04:02:00Z\\\\\n",
       "\t6 & 12.75 & -121.9853 & 36.84600 & 2013-05-11T04:04:53Z\\\\\n",
       "\\end{tabular}\n"
      ],
      "text/markdown": [
       "\n",
       "A data.frame: 6 × 4\n",
       "\n",
       "| <!--/--> | sst &lt;dbl&gt; | longitude &lt;dbl&gt; | latitude &lt;dbl&gt; | date &lt;fct&gt; |\n",
       "|---|---|---|---|---|\n",
       "| 1 | 12.78 | -121.9688 | 36.83320 | 2009-05-21T05:12:46Z |\n",
       "| 2 | 12.78 | -121.9765 | 36.83096 | 2000-02-10T21:24:00Z |\n",
       "| 3 | 12.75 | -121.9790 | 36.84930 | 2011-05-21T10:06:00Z |\n",
       "| 4 | 12.78 | -121.9467 | 36.81833 | NA                   |\n",
       "| 5 | 12.75 | -121.9817 | 36.84500 | 2000-05-12T04:02:00Z |\n",
       "| 6 | 12.75 | -121.9853 | 36.84600 | 2013-05-11T04:04:53Z |\n",
       "\n"
      ],
      "text/plain": [
       "  sst   longitude latitude date                \n",
       "1 12.78 -121.9688 36.83320 2009-05-21T05:12:46Z\n",
       "2 12.78 -121.9765 36.83096 2000-02-10T21:24:00Z\n",
       "3 12.75 -121.9790 36.84930 2011-05-21T10:06:00Z\n",
       "4 12.78 -121.9467 36.81833 NA                  \n",
       "5 12.75 -121.9817 36.84500 2000-05-12T04:02:00Z\n",
       "6 12.75 -121.9853 36.84600 2013-05-11T04:04:53Z"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>183</li><li>4</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 183\n",
       "\\item 4\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 183\n",
       "2. 4\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 183   4"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>'sst'</li><li>'longitude'</li><li>'latitude'</li><li>'date'</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 'sst'\n",
       "\\item 'longitude'\n",
       "\\item 'latitude'\n",
       "\\item 'date'\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 'sst'\n",
       "2. 'longitude'\n",
       "3. 'latitude'\n",
       "4. 'date'\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] \"sst\"       \"longitude\" \"latitude\"  \"date\"     "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "##Run above function to get Sebastes data within the Soquel Canyon\n",
    "\n",
    "Sebastes_SoquelCanyon <- EnvironmentalMetricPresent(extract_polygon_geometry(soquel_canyon_smca), \n",
    "                                           \"Sebastes\", \n",
    "                                           \"sst\")\n",
    "head(Sebastes_SoquelCanyon)\n",
    "dim(Sebastes_SoquelCanyon)\n",
    "names(Sebastes_SoquelCanyon)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       ".list-inline {list-style: none; margin:0; padding: 0}\n",
       ".list-inline>li {display: inline-block}\n",
       ".list-inline>li:not(:last-child)::after {content: \"\\00b7\"; padding: 0 .5ex}\n",
       "</style>\n",
       "<ol class=list-inline><li>12.73</li><li>12.78</li></ol>\n"
      ],
      "text/latex": [
       "\\begin{enumerate*}\n",
       "\\item 12.73\n",
       "\\item 12.78\n",
       "\\end{enumerate*}\n"
      ],
      "text/markdown": [
       "1. 12.73\n",
       "2. 12.78\n",
       "\n",
       "\n"
      ],
      "text/plain": [
       "[1] 12.73 12.78"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "## Get range of sst values for Sebastes sightings\n",
    "\n",
    "rangeValue <- range(Sebastes_SoquelCanyon$sst)\n",
    "rangeValue"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Visualize Environmental Data from OBIS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "options(repr.plot.width=10, repr.plot.height=8)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABLAAAAPACAMAAAASenCyAAAC91BMVEUAAAABAQECAgIDAwME\nBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUW\nFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJyco\nKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY4ODg5OTk6Ojo7\nOzs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExN\nTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5f\nX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBx\ncXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKD\ng4OEhISFhYWGhoaHh4eIiIiKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OVlZWWlpaX\nl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKip\nqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7\nu7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzN\nzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f\n39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx\n8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7///+EC70FAAAACXBI\nWXMAABJ0AAASdAHeZh94AAAgAElEQVR4nO3de5zcdX3v8R/hJoRorWLRivRiK1pUKhWqR61y\nqrUlBIuiJhLuiGCrgi3YeIGqoNSDIlpRa2utlqKnWsUiCp7qOVSlFVSEokCRqzDZJCRAgFw+\nf5yZ2d1k3ruTz34/4ZP9zXfzev2xOzs7l+++d/N87G4iNkZEVElN2wcgIioNsIiomgCLiKoJ\nsIiomgCLiKoJsIiomgCLiKoJsOZKpzcf3aaPv6z58Fbca6tPdfVzd55/y9bdNd7WfWzUQoA1\n+l3dbO5ft3ir9+3zmcwn/fSXplwx8Ie6e6BFk5e/3DTLnLuGTjV4x/2aF7/lni3d8F0Dizym\n/PF/9Ib9f2mnPX7rpOumvqMErGmDUBsB1uh3dTN/0WTfn60n3eu4KVcIWDvsdPfE5VfuMA2s\naXfdiud8sHn0ui3f8HO9KQ5s9uy9Wlz88J/ZuXn20Se/8knNrl+d8p4SsLb6o6LMAGv0u7r5\n9Vl/zpsaF6z9mw+MX1yx635TwZp+1614zpXNPjPd/J+bg0MPv2r+jp/vvd5wVvOUKRgWgLXV\nHxWlBlij3xSw3t586drD9tz12Z81e1HzL/2rvty8ZPy3RX/RfOmCJ3V/SHr4g7+zx66/fsod\nevPe5StftMfjj1q98byn7fb0szd2r9r48YP2eNS+b79Pb3l474et55vde8a+j9rlqW+916aA\nddye+41f/Ejzlj5Ymx9m4q4TJxn/HdaGC35n/h4v+bfNH8PFL37szk/8g6/Kxcnn7LWo/9Pe\nTwc/jskPbVOTYG0+47fn/cba7hWrn7Tz9+Xo432zOWji0lGn3qkf+bLmgiteuGD+87+ujzfk\ncAN3GvwgaLYCrNFvClh/2Zz5mJe+5Y+a5nL7WHNU/6ojm0+N03Bmc+rurz3BNry82fdN73h5\n88Rb5Obdy2/7hVef8uTmyNP3PmnpLs2nu3d9XfPE0/7ioGb/1XLLrxzdHHTeP9nDL2gOeOub\nntY8d/0UsI49qbmqf/HAX/qnPlibH2birhMnGQfrlc0zTl6yR//5+l3Y7Pn6dxz7izt8evDi\nxB37/cvZzWPPO2/V4Mcx8YCbmwBr8IynNe/oXvOnzZl69Xjfb35zw+COAx/5suaEXf/g1EU7\n7PR/9I7TDzdwp4EPgmYtwBr9poB1TrPL33dfvbWL1YpdHvtw9+KDj9lt9TgNZzeP+Zr1/iz9\n7oPW+5bpCLl59/Ku3zT72Y477ztm9onmELOLmgO6f/Y2vrE5Q295cf8noC80B3X/2D60b+8b\nOQHrqO81J/cuXd+cenEPrMGHGb/rxEn6p/pc8/Luw/zX7vPXTDzAM5sbuy9vW3CQXLx42o+E\ngx/HxANubgKswTOuffou19v3d3zOOr16vHX7Ni8d+BXg4JGXNfN6v1I/t/c91OAdpx1u8E4D\n76RZC7BGv6ub+YdP9CHrudL/uem7ze+aLWwu7V78YvOaCRrOGf+x5/lN/weVe3fZ5QG5+TnN\ny3qX928u6L78efN0s//ZXNa7ZuXOT9QHHv/z+d//u/991OnNe6aCZc98bI+SM5pr+2ANPsz4\nXSdO0j/VS5tv9y6fd9pNEw+w9w539V49pBeng6UfxxQZJsAaPKN9b6cXrT9w12unXj3e9c9o\nmn2WXPDD8bcGj7xs/LHXPmqHMbnjtMMN3mngnTRrAdboN/DPGpZYj4LTetf+pNm/973L8d2L\ni5tLNoH1lu7bGx/VrOrf85ndH9wGb35Oc3rv8gv7gDzQ/IrZHk3/1zH2283P5Jab8Vh9111n\n9UyaAtYHm8+ZbXjygdYHa/BhJsHqnWT8VPObB/QjOrnZ95N3Tbs4DawpH8db9EEGfuk+ecbe\nN2Mva943/erx1l90xF7dEfd+T8/awSMva97av7x/c6XccdrhBu808E6atQBr9Jv2I+G5vVc/\nbZ5tdv8ej19vaxc8Yd0msM7uvmt1s8v4TV/S/KvcfOLyi5rruy/Xdk14YLOFV8otJ/D45+c/\nqv/O6WAt3/WlZl/rPmcPLHmYSbB6J+mf6r7mUVM+oodfv3PTPOPPb9aL08Ca8nGcrQ8yCdbA\nGbsP91vNc9ZPv3pzN336iPnNgQ/pkZdN/KXnS3ryD9xx6uHkTgPvpFkLsEY/Byxb0lxuX2j+\n1DaB1Xvfmmbn8Zv+XvcnRhestc0O75ro5iFgfaxZ8ObPXvKvrx8Clh0x71ZbvNuqPljyMJNg\n9R+sd6oHmh03Tv2Y7vzrRQuaXS6Si9PAGvZxbG4CrMEzmt366OYXbp9+tfTz/bqHkiMvaz7Y\nf8/BXRj1jno4uZN8EDRLAdbo54F1SXOyvab/N3YDYNnuzcr+TX+r+U8fLHtMc8/QBx7H45eb\n/j9FeOcwsL7WvHf1bkts/EfCwYeZBpYtaDpDPqy1H93pFx4cvDj9d1hDPo7NTYA1eEaz3593\nTvMH06/ut3ySzY80x+qRlzVv77/ev/nu9DsOHm7wTlM/CJqNAGv088Ba9/i9H5i/b+/tQbBe\n2Hyl92psp90enAGslzb9f0xpYzYErAebPXpXbDxwGFgbnvLcT/f+qUQfrMGHmQ7Wi3u/7zI7\n++D/N/EAt9xp48e8bvDidLCGfBybGwdLzmh/3fyJHdtcOPXqfs9rJv/J1Jm9v+IcPPKy5vd7\nFx/cbd69csdphxu808A7adYCrNHPA8ve0Lx3/C/CBsH6u+Z5vb+8Oq33K3kfrIua/XrfM3xr\np1fqLb/cHNa9+IvNrd0/u2fu1ZwyHSx717yDf3XjBFiDDzN+10Gw/rY54D6z/37s7ivG739N\n85Le8VY/ecflAxcn7jjeOFhDPo7NTXyHNXjGWxY8ZY2teMKCW/Tq8d7T7Nn/K76Nn9+9+aYe\neVkzr2fph3uPOHDH6YcbuNPgO2nWAqzR7+pm9z+abNlUsL7dLNjhv3sXBsHauKjZ78/+4uDm\nN5fPBJa9tvnlU991+M4Lvqu3vHGHnY89yU5tfuM97znoaV9rHve+26aBdcu85iybAGvwYcbv\nOgjWhkOafU46ckHz8ckHWNz82invPHmf5k1ycfyO442DNeTj2NwEWANnvLX/S3P7XPPijXL0\n8duvX9I0v/qKJX+4dzPvvaZHPr054dFHnfWaHXf9jjzebdMPN3CngXfSrAVYo9/gf63h4Klg\nbfyV5gX9Ww2CZes++Jzdd933bb3fAM0A1oaP/+6CnZ689Popt7T3PX7X59jaZb++694nL7ej\n5+/1w2lg2e/P+5lNgjXwMON3HQTL1n3gWbvNf+EVmz6iDRc87/E7PuYFf7NRLo7fcbyJ/y3h\n9I9jcxNgDZzxpOa1/ff8YXO+HH3yHl9b8tTd5j36WaeMXzFw5Dc3F3/zhXvMf9G39fF+OP1w\nA3cafCfNVoBFRNUEWERUTYBFRNUEWERUTYBFRNUEWERUTYBFRNUEWERUTYBFRNUEWERUTYBF\nRNUEWERUTYBFRNUEWERUTY8QrPPeMGOvP+GkmW/UQie2fYChnTiaa530+rZPMLQTR/OT+PrR\n/ByeMJprlZzqHVlgnXjV7TN1449unvE2LXTrdW2fYGg3/OhnbR9hWDf/pO0TDO3aa9s+wdB+\nMppf8T+6oe0jDO26W2e8yS2HpoF164w3eaAzkv+R/g0r2z7B0FZ31s98o9nvwfvaPsHQxsba\nPsHQ7hvNr/jO6raPMLSVG2a8ycOABVihACsSYEUCrKIAKxJgRQKsSIBVFGBFAqxIgBUJsIoC\nrEiAFQmwIgFWUYAVCbAiAVYkwCoKsCIBViTAigRYRQFWJMCKBFiRAKsowIoEWJEAKxJgFQVY\nkQArEmBFAqyiACsSYEUCrEiAVRRgRQKsSIAVCbCKAqxIgBUJsCIBVlGAFQmwIgFWJMAqCrAi\nAVYkwIoEWEUBViTAigRYkQCrKMCKBFiRACsSYBUFWJEAKxJgRQKsogArEmBFAqxIgFUUYEUC\nrEiAFQmwigKsSIAVCbAiAVZRgBUJsCIBViTAKgqwIgFWJMCKBFhFAVYkwIoEWJEAqyjAigRY\nkQArEmAVBViRACsSYEUCrKIAKxJgRQKsSIBVFGBFAqxIgBUJsIoCrEiAFQmwIgFWUYAVCbAi\nAVYkwCoKsCIBViTAigRYRQFWJMCKBFiRAKsowIoEWJEAKxJgFQVYkQArEmBFAqyiACsSYEUC\nrEiAVRRgRQKsSIAVCbCKAqxIgBUJsCIBVlGAFQmwIgFWJMAqqhKwXkuSrgVYgQDLACs9wPLT\ntQArEGAZYKUHWH66FmAFAiwDrPQAy0/XAqxAgGWAlR5g+elagBUIsAyw0gMsP10LsAIBlgFW\neoDlp2sBViDAMsBKD7D8dC3ACgRYBljpAZafrgVYgQDLACs9wPLTtQArEGAZYKUHWH66FmAF\nAiwDrPQAy0/XAqxAgGWAlR5g+elagBUIsAyw0gMsP10LsAIBlgFWeoDlp2sBViDAMsBKD7D8\ndC3ACgRYBljpAZafrgVYgQDLACs9wPLTtQArEGAZYKUHWH66FmAFAiwDrPQAy0/XAqxAgGWA\nlR5g+elagBUIsAyw0gMsP10LsAIBlgFWeoDlp2sBViDAMsBKD7D8dC3ACgRYBljpAZafrgVY\ngQDLACs9wPLTtQArEGAZYKUHWH66FmAFAiwDrPQAy0/XAqxAgGWAlR5g+elagBUIsAyw0gMs\nP10LsAIBlgFWeoDlp2sBViDAMsBKD7D8dC3ACgRYBljpAZafrgVYgQDLACs9wPLTtQArEGAZ\nYKUHWH66FmAFAiwDrPQAy0/XAqxAgGWAlR5g+elagBUIsAyw0gMsP10LsAIBlgFWeoDlp2sB\nViDAMsBKD7D8dC3ACgRYBljpAZafrgVYgQDLACs9wPLTtQArEGAZYKUHWH66FmAFAiwDrPQA\ny0/XAqxAgGWAlR5g+elagBUIsAyw0gMsP10LsAIBlgFWeoDlp2sBViDAMsBKD7D8dC3ACgRY\nBljpAZafrgVYgQDLACs9wPLTtQArEGAZYKUHWH66FmAFAiwDrPQAy0/XAqxAgGWAlR5g+ela\ngBUIsAyw0gMsP10LsAIBlgFWeoDlp2sBViDAMsBKD7D8dC3ACgRYBljpAZafrgVYgQDLACs9\nwPLTtQArEGAZYKUHWH66FmAFAiwDrPQAy0/XAqxAgGWAlR5g+elagBUIsAyw0gMsP10LsAIB\nlgFWeoDlp2sBViDAMsBKD7D8dC3ACgRYBljpAZafrgVYgQDLACs9wPLTtQArEGAZYKUHWH66\nFmAFAiwDrPQAy0/XAqxAgGWAlR5g+elagBUIsAyw0gMsP10LsAIBlgFWeoDlp2sBViDAMsBK\nD7D8dC3ACgRYBljpAZafrgVYgQDLACs9wPLTtQArEGAZYKUHWH66FmAFAiwDrPQAy0/XAqxA\ngGWAlR5g+elagBUIsAyw0gMsP10LsAIBlgFWeoDlp2sBViDAMsBKD7D8dC3ACgRYBljpAZaf\nrgVYgQDLACs9wPLTtQArEGAZYKUHWH66FmAFAiwDrPQAy0/XAqxAgGWAlR5g+elagBUIsAyw\n0gMsP10LsAIBlgFWeoDlp2sBViDAMsBKD7D8dC3ACgRYBljpAZafrgVYgQCr23E/WjFTyztj\nM96mjTptH2BoU9dqG4hRS9fqjOgncUS/4pe3fYKhFXwO7z4kCyy+w0qO77D8dC2+wwrEd1gG\nWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqA\nFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A\n8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGW\nAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroW\nYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoP\nsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECA\nZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+u\nBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDS\nAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQ\nYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWn\nawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs\n9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAK\nBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD5\n6VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsA\nKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuw\nAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdY\nfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAy\nwEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cC\nrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkB\nlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiw\nDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1\nACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6\ngOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUC\nLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0\nLcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCV\nHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViB\nAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/\nXQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlg\npQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAm0nYN1+2qLeqzUf\nOGrxWXdvfj0RYCUHWH66FmAF2j7A+tbS8/pgvfv0m+8495QNm15PBFjJAZafrgVYgbYPsC6/\n5997YHUOvan73dVh10y+nnw3YCUHWH66FmAF2j7AMuuDdeXhG7sv33jR5Ovui3X3djvhZxtn\n6v7O2hlv00LrV7Z9gqGt7qyTt9sGYtTStcbGZvFTU959o/kV31nd9hGGtnL9jDd5KAzWpUf3\nLi67cPJ198UVB3Q74gcd2pa1DcSo1fbng9rozkPCYB0zAdYxm8C65g3dFv941UyNdVbMeJsW\nWrm87RMMbayzUt5uG4hRS9fqdGbxU1Pe2Eh+xa/qjLV9gqEtXznjTTphsL4z/qPgxZOvJ9/J\n77CS43dYfroWv8MKtF39Dmvs0J+a3bvo2snXk+8ErOQAy0/XAqxA2wdYKzqXLep01to5b775\n9jNP3bjp9USAlRxg+elagBVo+wDruIW9vmT3n7d0ydkrbNPriQArOcDy07UAK9D2AdYMAVZy\ngOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUC\nLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0\nLcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCV\nHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViB\nAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/\nXQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlg\npQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFW\nIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADL\nT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgG\nWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqA\nFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A\n8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGW\nAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroW\nYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoP\nsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECA\nZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+u\nBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDS\nAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQ\nYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWn\nawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs\n9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAK\nBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD5\n6VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsA\nKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuw\nAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdY\nfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAy\nwEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cC\nrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkB\nlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECA1e2EGx+cqdWdNTPe\npoXWrmj7BENb1XlA3m4biFFL11q+fBY/NeXdO5pf8Z1VbR9haCvWzniTNQuzwDr+hvtmalXn\n3hlv00Jrxto+wdBWdFbL220DMWrpWsuXz+KnpryVo/kV31nR9hGGNrZmxpusTAOLHwmT40dC\nP12LHwkD8SOhAVZ6gOWnawFWoDkI1gHXjb/+/NNLnwywkgMsP10LsALNQbCaq/qv1p21S+mT\nAVZygOWnawFWoDkHVrO555Q+GWAlB1h+uhZgBZpzYF3zoWbRcb2Of+dtpU8GWMkBlp+uBViB\n5hxYZi/7SfTJACs5wPLTtQAr0BwEKx5gJQdYfroWYAWag2DdfdST5o3/Eqv0yQArOcDy07UA\nK9AcBOtVOx18VP+3WMeVPhlgJQdYfroWYAWag2A97ovRJwOs5ADLT9cCrEBzEKzd74k+GWAl\nB1h+uhZgBZqDYL3gm9EnA6zkAMtP1wKsQHMQrP848MrgkwFWcoDlp2sBVqA5CNbz925236df\n6ZMBVnKA5adrAVagOQjWCw6erPTJACs5wPLTtQAr0BwEKx5gJQdYfroWYAUCLAOs9ADLT9cC\nrEBzEKzHTQLBrU8AABxzSURBVLag9MkAKznA8tO1ACvQHARrUb8Dd9vvlNInA6zkAMtP1wKs\nQHMQrInueuFXSp8MsJIDLD9dC7ACzV2w7KoDSp8MsJIDLD9dC7ACzWGw7tqt9MkAKznA8tO1\nACvQ3AVr43ufXPpkgJUcYPnpWoAVaA6C9ex++z2+eWvpkwFWcoDlp2sBVqA5C9Zvv+RDD5U+\nGWAlB1h+uhZgBZqDYMUDrOQAy0/XAqxAcxKs5V+58BOXBj4uwEoOsPx0LcAKNAfB2nDazr3/\noPv89xc/GWAlB1h+uhZgBZqDYL2/ecUnv/qVj72s+bvSJwOs5ADLT9cCrEBzEKynnzr++kT+\nn5/bCrD8dC3ACjQHwdr18vHXl/APR9sKsPx0LcAKNAfBmv/l8ddf3KP0yQArOcDy07UAK9Ac\nBOt/vLj/D7DWvvT3Sp8MsJIDLD9dC7ACzUGwLtnhKSe9+y9PeNK8r5c+GWAlB1h+uhZgBZqD\nYNk/79v7Zw3PvKT4yQArOcDy07UAK9BcBMvsju9d9fPAkwFWcoDlp2sBVqC5CNZd53df3HPW\n3cVPBljJAZafrgVYgeYgWP+11y7dl7c0e91U+mSAlRxg+elagBVoDoJ12FO/13t13VP/uPTJ\nACs5wPLTtQAr0BwEa8+/GX/9Mf5fc9oKsPx0LcAKNAfB2u0z46//YffSJwOs5ADLT9cCrEBz\nEKznvaz/p2X1c59f+mSAlRxg+elagBVoDoJ16Q6/dsqZ7zxmz3mXlj4ZYCUHWH66FmAFmoNg\n2WUH9P7h6LP4h6OtBVh+uhZgBZqLYJkt/+GPIx8WYCUHWH66FmAFmptgBQOs5ADLT9cCrECA\nZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+u\nBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDS\nAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQ\nYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWn\nawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs\n9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAK\nBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD5\n6VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsA\nKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuw\nAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdY\nfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAy\nwEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cC\nrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkB\nlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiw\nDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1\nACsQYBlgpQdYfroWYAUCLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6\ngOWnawFWIMAywEoPsPx0LcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAUC\nLAOs9ADLT9cCrECAZYCVHmD56VqAFQiwDLDSAyw/XQuwAgGWAVZ6gOWnawFWIMAywEoPsPx0\nLcAKBFgGWOkBlp+uBViBAMsAKz3A8tO1ACsQYBlgpQdYfroWYAXavsD6k4XdXmW25gNHLT7r\n7k1XA1ZygOWnawFWoO0LrGO+3Ol0ul8e7z795jvOPWXTUwBWcoDlp2sBVqDtC6xXXtV/1Tn0\npu53WYddM3k1YCUHWH66FmAF2q7Aenjh+W869uzb7crDN3bfeuNFk9cDVnKA5adrAVag7Qqs\nVUf+rxtuOPPI+y49uvfWsgu7L246v9vSG+6bqVWde2e8TQutGWv7BENb0Vktb7cNxKilay1f\nPoufmvJWjuZXfGdF20cY2tiaGW+ycmEUrH4PvOqyS4/ZBNYVB3Q74gcd2pa1DcSo1fbng9ro\nzkO2Ciw7+bPfGf+R8OLui7Hvdjvm5odnak3n/hlv00IPrWj7BENb1XlQ3m4biFFL11q+fBY/\nNeWtHs2v+M6qto8wtBUPzXiT+8PfYd3y4XVma191xdihPzW7d9G1k9fzO6zk+B2Wn67F77AC\nbVe/w1q9+Ly7bj/7mAftnDfffPuZp26cvB6wkgMsP10LsAJtV2DZTW9/9eve/XOz+89buuTs\nFZuuBqzkAMtP1wKsQNsXWFsIsJIDLD9dC7ACAZYBVnqA5adrAVYgwDLASg+w/HQtwAoEWAZY\n6QGWn64FWIEAywArPcDy07UAKxBgGWClB1h+uhZgBQIsA6z0AMtP1wKsQIBlgJUeYPnpWoAV\nCLAMsNIDLD9dC7ACAZYBVnqA5adrAVYgwDLASg+w/HQtwAoEWAZY6QGWn64FWIEAywArPcDy\n07UAKxBgGWClB1h+uhZgBQIsA6z0AMtP1wKsQIBlgJUeYPnpWoAVCLAMsNIDLD9dC7ACAZYB\nVnqA5adrAVYgwDLASg+w/HQtwAoEWAZY6QGWn64FWIEAywArPcDy07UAKxBgGWClB1h+uhZg\nBQIsA6z0AMtP1wKsQIBlgJUeYPnpWoAVCLAMsNIDLD9dC7ACAZYBVnqA5adrAVYgwDLASg+w\n/HQtwAoEWAZY6QGWn64FWIEAywArPcDy07UAKxBgGWClB1h+uhZgBQIsA6z0AMtP1wKsQIBl\ngJUeYPnpWoAVCLAMsNIDLD9dC7ACAZYBVnqA5adrAVYgwDLASg+w/HQtwAoEWAZY6QGWn64F\nWIEAywArPcDy07UAKxBgGWClB1h+uhZgBQIsA6z0AMtP1wKsQIBlgJUeYPnpWoAVCLAMsNID\nLD9dC7ACAZYBVnqA5adrAVYgwDLASg+w/HQtwAoEWAZY6QGWn64FWIEAywArPcDy07UAKxBg\nGWClB1h+uhZgBQIsA6z0AMtP1wKsQIBlgJUeYPnpWoAVCLAMsNIDLD9dC7ACAZYBVnqA5adr\nAVYgwDLASg+w/HQtwAoEWAZY6QGWn64FWIEAywArPcDy07UAKxBgGWClB1h+uhZgBQIsA6z0\nAMtP1wKsQIBlgJUeYPnpWoAVCLAMsNIDLD9dC7ACAZYBVnqA5adrAVYgwDLASg+w/HQtwAoE\nWAZY6QGWn64FWIEAywArPcDy07UAKxBgGWClB1h+uhZgBQIsA6z0AMtP1wKsQIBlgJUeYPnp\nWoAVCLAMsNIDLD9dC7ACAZYBVnqA5adrAVYgwDLASg+w/HQtwAoEWAZY6QGWn64FWIEAywAr\nPcDy07UAKxBgGWClB1h+uhZgBQIsA6z0AMtP1wKsQIBlgJUeYPnpWoAVCLAMsNIDLD9dC7AC\nAZYBVnqA5adrAVYgwDLASg+w/HStSsBqe7RRS9cBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQd\nwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCq\nTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq\n03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKsKtN1AKso\nwKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcCrCrTdQCrKMCqMl0L\nsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcB\nrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvK\ndC3AqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoy\nXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoC\nrCrTtQCrynQdwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UA\nq8p0HcAqCrCqTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3A\nKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM\n1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcCrCrT\ndQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq03UAqyjA\nqjJdC7CqTNeZVbCO+9GKmVreGZvxNm3UafsAQ5u6VttfW6OWrtUZ0U/ilK/4tkcbtXSdgs/h\n3YdkgXXizzbO1P2dtTPepoXWr2z7BENb3Vknb7f9tTVq6VpjY7P4qSnvvilf8W2PNmrpOivX\nzzjoQ/xIyI+EdaZr8SNhlek6/A6rKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvK\ndB3AKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoK\nsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcC\nrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq03UA\nqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKsKtN1AKsowKoy\nXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcCrCrTdQCrKMCqMl0LsKpM\n1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIA\nq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3A\nqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQew\nigKsKtO1AKvKdB3AKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrT\ntQCrynQdwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0\nHcAqCrCqTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqw\nqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKs\nKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcCrCrTdQCr\nKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq03UAqyjAqjJd\nC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzX\nAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCqTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCr\nynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCq\nMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CK\nAqwq07UAq8p0HcAqCrCqTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1\nAKvKdB3AKgqwqkzXAqwq03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQd\nwCoKsKpM1wKsKtN1AKsowKoyXQuwqkzXAayiAKvKdC3AqjJdB7CKAqwq07UAq8p0HcAqCrCq\nTNcCrCrTdQCrKMCqMl0LsKpM1wGsogCrynQtwKoyXQewigKsKtO1AKvKdB3AKgqwqkzXAqwq\n03UAqyjAqjJdC7CqTNcBrKIAq8p0LcCqMl0HsIoCrCrTtQCrynQdwCoKsKpM1wKsKtN1Rh6s\ntucasXQcwPLTtQCrynQdwKorHQew/HQtwKoyXQew6krHASw/XQuwqkzXAay60nEAy0/XAqwq\n03UAq650HMDy07UAq8p0HcCqKx0HsPx0LcCqMl0HsOpKxwEsP10LsKpM1wGsutJxAMtP1wKs\nKtN1AKuudBzA8tO1AKvKdB3AqisdB7D8dC3AqjJdB7DqSscBLD9dC7CqTNcBrLrScQDLT9cC\nrCrTdQCrrnQcwPLTtQCrynQdwKorHQew/HQtwKoyXQew6krHASw/XQuwqkzXAay60nEAy0/X\nAqwq03UAq650HMDy07UAq8p0HcCqKx0HsPx0LcCqMl0HsOpKxwEsP10LsKpM1wGsutJxAMtP\n1wKsKtN1AKuudBzA8tO1AKvKdB3AqisdB7D8dC3AqjJdB7DqSscBLD9dC7CqTNcBrLrScQDL\nT9cCrCrTdQCrrnQcwPLTtQCrynQdwKorHQew/HQtwKoyXQew6krHASw/XQuwqkzXAay60nEA\ny0/XAqwq03UAq650HMDy07UAq8p0nW0L1poPHLX4rLs3vQlYjzgdB7D8dC3AqjJdZ9uC9e7T\nb77j3FM2PQVgPeJ0HMDy07UAq8p0nW0KVufQm7rfZR12zeTbgPWI03EAy0/XAqwq03W2KVhX\nHr6x+/KNF02+DViPOB0HsPx0LcCqMl1nm4J16dG9l8su7L644oBuR/ygE67tuUYsxooU/3Ib\ngdoebdQKD3jnIVsN1jGbwLrmDd0W/3jVTI11Vsx4mxZaubztEwxtrLOy7SMMa8VY2ycYWqfT\n9gmGNjaSX/GrOqP5SVw+81d8Z6vB+s74j4QXT769NT8SjkYbVrZ9gqFN/ZFwRHrwvrZPMLRK\nfiQcjTZ0Vrd9hKFt0x8Jxw79qdm9i66dfBuwkgOsSIAVaLsEy8558823n3nqxsk3ASs5wIoE\nWIG2T7DuP2/pkrNXbHoTsJIDrEiAFWj7BGtKgJUcYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKs\nQIBlgJUeYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKsQIBlgJUeYEUCrECAZYCVHmBFAqxAgGWA\nlR5gRQKsQIBlgJUeYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKsQIBlgJUeYEUCrECAZYCVHmBF\nAqxAgGWAlR5gRQKsQIBlgJUeYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKsQIBlgJUeYEUCrECA\nZYCVHmBFAqxAgGWAlR5gRQKsQIBlgJUeYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKsQIBlgJUe\nYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKsQIBlgJUeYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKs\nQIBlgJUeYEUCrECAZYCVHmBFAqxAgGWAlR5gRQKsQIDV7cQP/e1MffyCT8x4mxb61EfbPsHQ\nPnbB37R9hGF94sK2TzC0j3yk7RMM7cLR/Iq/4GNtH2FoH/3UjDf5ZBpY//aFGfurY8+f+Uaz\n3+cvavsEQ1t27N+3fYRhXXxx2ycY2sknt32CoY3mWn9/7LK2jzC0iz4/822+lgVWQf94wFe3\n+XPMnd52wB1tH6GiXv7ytk9QUXce8La2j5AQYI1WgBUJsAIBVlmAFQmwIgFWIMAqC7AiAVYk\nwAoEWEREsxpgEVE1ARYRVRNgEVE1bQOwbj9tUe/V2LmvO+KMG/rX/HBhv68MXkfjOWt1+8bC\nf2/xbKOXt9Ylx7/ijd9r9XSjlrPWbX+55NVn/LjV021d+WB9a+l5/ZnecvpNd/7VkrW9iw93\nuv34VbcOXkf9vLXMVh55OGAN5K31jaVX3f3FE+5v94AjlbPWxhM+fP+Dn3nVaP6PC93ywbr8\nnn/vzbT67O6fuHsW/mTT9W//7PTryFmr2zmfPBKwBvLWOuHy9s41mjlrrVp4vdmKhRX+sLMt\nfofVn6nf9YtWTF781nHrpl1H5q515fFrAUvb4lrLF17+J6887fp2TjWqbflr68/OW732s8c/\n1M6xHknbFKzVJ39q8roNJ3192nXUa8trrVl6tQGWtsW1blj4tttWX/iaVe0ca0Tb8tfW2CkL\nFy69sZVDPbK2JVi3nfiRjZPXfevo9dOuo15bXuuDHzTAmtIW17ph4TVm61/7jbYONpJtca11\nb/rwqvsvXlLhzzrbEKxrFn9583VnXTj9Ouq1xbWuXroasKa2xbU6C3/afXnKxa2calTb4lrf\nP7T3O/hj/6WVUz2ith1YP37tf2y+6r7Drpl2HfXb4lrvP3zx4sWHHnF2S+cazba41oal3T+U\nDx3xrZbONZptca3/XNj769SlgNVtReeyRZ3O2odO+Fzvr1DX2mW9Va5ZeHf35abraDJnrdW9\na1532b1tH3GEctayi5dc3Tl/KV9bm3PWun/ph9c89IXD72z7iPHywTqu/0/TvnTN5D9Re//b\nu1d+89De33pdM/BPIqmfs1Y/fiQczFtrw98d+YozZv6/GNiO8ta65awlr/nzH7Z8wK2J/2kO\nEVUTYBFRNQEWEVUTYBFRNQEWEVUTYBFRNQEWEVUTYBFRNQEWEVUTYFFLXc3XHoXji4Za6ny+\n9igcXzQ0K915/FN2/aU/vn7zhZc1TXNA26ei2gIsmpUO2usTV/zDM59w/6YLP1nUXHVd26ei\n2gIsmo3ubc7ovrzx7Ds2XbDj+NqjcHzR0Gz08OP2+cYGuQBYtBXxRUOz0v/91eZxh//DuoEL\ngEXx+KKh2Wn95W99RvM7D2y+AFgUjy8amr0+0vzt5guARfH4oqHZ6D9e3ftPid/YnLvpgh3f\nrJvpXkRTAiyaje5a8KxPfv0fn/foGzddsHc2Z32+7WNRbQEWzUo/eMUTdn7SK74/cOG23975\naW2fimoLsIiomgCLiKoJsIiomgCLiKoJsIiomgCLiKoJsIiomgCLiKoJsIiomgCLiKoJsIio\nmgCLiKrp/wO08nU3haQ6mQAAAABJRU5ErkJggg==",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 480,
       "width": 600
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#input: environmental dataframe, taxa, and metric of interest\n",
    "#output: histogram of metric values\n",
    "\n",
    "GraphEnvironmentalMetric <- function(inputDF, taxa, metric) {    \n",
    "    \n",
    "EnvironmentalHistogram <- \n",
    "    ggplot(inputDF, aes_string(x = metric)) +\n",
    "    geom_bar(stat = \"count\") +\n",
    "    theme_bw() +\n",
    "    labs(title = paste(\"Environmental Metric\", metric, \"for Taxa\", taxa)) +\n",
    "    theme(plot.title = element_text(hjust = 0.5))\n",
    "    \n",
    "   return(EnvironmentalHistogram)\n",
    "\n",
    "}\n",
    "\n",
    "SebastesSST <- GraphEnvironmentalMetric(Sebastes_SoquelCanyon, \"Sebastes\", \"sst\")\n",
    "\n",
    "SebastesSST"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Map taxa Sightings with associated Environmental Factor"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABLAAAAPACAIAAACqxhfXAAAACXBIWXMAABJ0AAASdAHeZh94\nAAAgAElEQVR4nOzdaWAT1d7H8TNJk6ZpS2lZaxdABMqOQAHBpQKigNgLD3tl33pZFEFWZQfh\nogi4AIKARbgsoiC7FxCqCIooKHoRFxYBqUAtJW1JmzTzvJhrrGmgobSdNPP9vGpOZk7+kzNp\n8+vMnJFkWRYAAAAAAO3RqV0AAAAAAEAdBEIAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAA\nGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAA\nADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAA\nAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIA\nAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiE\nAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYR\nCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACN\nIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAA\nGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAA\nADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAA\nAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIA\nAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiE\nAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACNIhACAAAAgEYR\nCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAAGkUgBAAAAACN\nIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAA\nGkUgBAAAAACNIhACAAAAgEYRCAFowsSJEyVJWrZsmdqFeKkXX3xRkqQ33nhD7UJuqRSN4IkT\nJ5o1a2Y0GoOCgs6fP692OaWV9++TAOAbCIQAVHDixAmpIHv27CnCVwwNDa1SpUpwcHAR9um1\n3n333W3bthVtn84h+8c//nGrZXbs2KEs8+KLL95R554U7CUj6Empffr0+fLLLx988MGhQ4ea\nzeaieunp06cX+KkpW7ZsUb1c4Xz33XfDhw+///77K1eubDAYgoOD69Wr989//vPUqVPqFnZ7\nxfGRAYDSwk/tAgBoV2BgYNu2bW/1bKVKlYrwtSZMmDBhwoQi7NCbjR8/vmPHjk899VSR9yxJ\n0s6dO69cuVKxYsX8zyYlJUmSJMvynXbrScFeMoIFlpqdnf3dd9+VKVPmP//5j59fUf6RjYmJ\niY+Pdz68fPny0aNHK1So0LJlS2djYGBgEb7inVq3bt2AAQNsNlvDhg3bt29vNpuvXLly+PDh\nZcuWrV69esuWLe3bt1exvNsovo8MAHg/AiEA1VSuXHnr1q1qV+Frzpw5k5KSUkydN2zY8MSJ\nE2vXrh0zZozLU2lpadu3b69bt+533313R30Wa8FFy5NSb968KYQIDQ0t2jQohOjZs2fPnj2d\nD7du3dq5c+cGDRp4yYcoPT192LBhDodj8+bN//d//+dsdzgcs2fPnjZtWmJi4i+//FLkb8vd\nK0V7IAAUB04ZBeDtpkyZIknStm3bvv/++86dO1esWNFkMjVq1Gj9+vXKAnFxcZIkbd++3WVF\n5QzGNm3aiHxXoL3wwgtKn0uWLImIiHCeaGez2RYvXhwbGxscHGwyme67776RI0f+9ttvnheT\nd5kjR47ExcUFBwdXqFChf//+FotFluVFixbFxMSYzeY6derMnTs37/E0WZbffvvtBx54IDg4\nOCAgoHbt2lOmTMnMzPT81bt27Vq9enUhxMqVKyVJevDBB5X2GzduTJo0qXbt2gEBAf7+/jVq\n1Bg3btyNGzfudCyaNGlSoUKF1atX539qw4YN2dnZjz32mEv77TfKbcFuRyf/NYQOh2PJkiWx\nsbFBQUHBwcFt2rT55JNPbl//5s2bW7duHRYWZjQa77nnnvbt2+/evdvDBW713ub1j3/8IzQ0\nVAhx/vx55RzOn3/+WXiwX91qhyyE24/1oUOH9Hp9zZo1rVarcxWLxRIREWE0Go8fP+5JJ24d\nP348MzMzNjY2bxoUQuh0uqlTp/br169r165Xr151the4tztXP3DgwCOPPFKmTJmgoKAHH3xw\n3759nm+vohDD6kl5Be5OAFA6yABQ4pTvndWrV/dk4ZkzZwohpk+fHhIS0q5du+eee65jx47K\nb7D9+/fLsvzWW28JIfr16+eyYp8+fYQQq1evlmVZOdtw6dKlylPTp08XQowZM8ZsNvfq1WvI\nkCGyLOfm5iqntMXExDz77LNTpkxRHoaHh587d87DYpzLTJo0qWzZsj169BgxYkRkZKQQok+f\nPhMmTIiKikpMTOzbt6/RaBRCrFmzxlnw008/rbzc2LFjJ0+e3KJFCyFEo0aNbty44eGr79ix\no3///kKIFi1aLFy4cNOmTbIs5+TkPPTQQ0KIJk2aPP/8888++2ytWrWEELGxsXa7Xen5hRde\nEEK8/vrrtx+ygQMHJiYmCiG+/PJLlwWaNWtWqVKlTZs2CSFeeOEFDzfKbcFuR8dlBGVZ7tq1\nqxCiTp06w4cPT0hICAoKcnk/XSxfvlwIUaFChWHDhk2ZMmXgwIFhYWGSJDlXuf0Cbkt1sW3b\ntpdeekkIERoaunDhwoULF16/ft2T/crtJt/eli1bhBBt2rTJ2+jJWI8dO1YIMWXKFOdazzzz\njLJfed5Jfl9//bUQombNmrm5uQUWL3uwtyv75JAhQ/z9/Z944okxY8bEx8dLkuTn53fw4EHP\nSy3csBZYXoG7EwCUFgRCACq4o0A4d+5cIYTRaHz33Xedjc8//7wzBP7xxx9GozE0NDQnJ8e5\ngNVqDQkJCQgIUL7AucQJ5Vt7SEjIRx995FxF+Yb3wAMPWK1WZ6MyP0r37t09LMa5jL+//4ED\nB5SW8+fP6/V6g8EQExOTmpqqNL799ttCiCeffFJ5uHHjRuV7rfMbp8PhGDlypBBi4sSJnr/6\ne++9J4QYNGiQc4H3339f+b7r/IqcnZ0dExMjhNi2bZvS4mEg7Nev39GjR4UQw4cPz/usMmXI\nmDFjlFd3BkJPNip/wW5Hx2UElYOi7du3d27UDz/8YDabAwMDlSOx+dWvX18I8fPPPztbLly4\nEBwc3KJFCw8XyF9qfmlpaUKIKlWqOFs82a/cbvLtuQ2Enoz1zZs3a9eubTQaT506Jcvy119/\nrdfrGzdubLPZPO8kP5vNpizTrl27r7/++vbFe7JjKPukTqf78MMPnSu+/PLLQohWrVp5Xmoh\nhtWT8grsFgBKC04ZBaCalJSUrrfw2muvuSwcGxur/M9e0a1bNyHEjz/+KIQIDQ19/PHH09LS\nPv74Y+cCe/bsSU9Pj4+PdzsvpSRJQojatWu3a9fO2ZiUlCSEmDJlir+/v7Nx3LhxRqNx69at\nyrVhBRbjFBcXFxcXp/wcHR1dv359m802atSosLAwpfHJJ58UQvzyyy/KwxUrVggh5s6d6yxY\nkqRZs2YZDAalsDt69bwaN278wQcfvP7663q9XmkxGo3K9CTffvvtrda6ldjY2Pr1669fvz47\nO9vZqFQ4cOBAl4U936i83I6OC+W01cmTJzs3qlatWnPmzElMTLxy5YrbVa5fvy5JUt5pVyIj\nI69du3bkyBEPFygcT/YrTzbZE56MtclkSkpKcjgciYmJubm5iYmJfn5+a9ascV7dV7gdxs/P\nb8uWLXXq1PnPf/7TuHHjqlWrPv3000uWLDl58mT+hT3fMZo1a5Z3rpeRI0eaTKbDhw//8ccf\nHpZaiGH1pLxi2lsAoOQRCAGoJjMz8/1bUA5D5aWcsuWkXKnlDGm9e/cWQmzevNm5gHLuonLW\n6K088MADzp9lWf7qq6+EEHnnbBRClClTplatWjk5Od9//72HxSgaNWrk0o8QokGDBi4tzrU+\n//zz/K9etmzZevXqXb58+ddff72jV8+ratWqnTt3btq0qRDCYrGkpKSkpKQot0O4zVq3MWjQ\noLS0NOUIlRDC4XCsXbu2WbNmdevWdVnS843KL+/o5PfZZ58JIZo0aZK3cfTo0a+88sq9997r\ndpVOnTrJsvzoo4+uWrXKOYmIcuKuhwsUwh3tV7ffZE94ONaxsbETJ05MTk7u2LHj0aNHZ8yY\nkXfsCr3DxMTEfPvttxs3buzevXt2dva6detGjBjRoEGD6OjoOXPm5P0Pguc7hsu1miaTKSYm\nRpbl06dPe1hqIYbVk/KKY28BAFV43WRfALSjevXqypQbnqhcuXLeh8oRFfnPGVmeeuqpoKCg\nrVu3Llu2TK/XW63W7du3V6xY8fbHWypUqOD8OSMjw2q1Go3GkJAQt4tdu3bNw2IU5cuXz79M\n3sa8a928eTMjI0MIoVwIl9+lS5eio6M9f3UXW7dufeWVV7766qu8U4kU2tNPPz1hwoTVq1cr\nk17u27fv4sWLygl+ed3RRuWXd3RcZGZmZmZmmkymgIAAz8tetGhRbm7uqlWrBg0aJISoU6fO\nk08+mZiYWK1aNQ8XKIQ72q9us8me83Csp06dumXLlo8++qhx48bKWceF6CQ/vV7fvXv37t27\nCyHOnDnz2Wef7dixY+fOnS+++OK2bds+/fRTo9F4RztGeHi4y7PKMXbl1FxPSr3TYfWwvOLY\nWwBAFQRCAL7AbDbHx8evW7cuOTm5devWu3btslgsAwYMuP0c9waDwfnzbWKVw+FwLlBMlM4l\nSZo6darbBVxC4B1Zvnz5sGHDgoODExMTmzVrFhISotPptm7dqkzGUwjlypWLj4/fvHnzhQsX\noqKikpKSAgICevXq5bLYXW5U3tFxodPphBDKNW+ej4vBYFi2bNm0adO2bdu2e/fujz/+eP78\n+YsWLXr33XeVAFPgAoVwR/vVbTbZQ56PdUpKyoULF8SfN12IiIgoRCe3d++999577719+vT5\n/fff27Zte/To0VWrViUmJt7RjuE8F9RJWV3ZBzwp9U6H1cPyimNvAQBVEAgB+IjevXuvW7fu\n/fffb926tTJLxO3PF3URFBRkNpuzsrKuX7/uMum/Mld+kRy9uRWTyRQSEpKenj5ixIgifyFl\nbtIdO3Y8/PDDzsa7vNJp0KBBmzZtevfdd0eNGrVly5YuXbrkPwJWfBsVEBAQHBxssVhSU1Nd\njsQWKDw8fNiwYcOGDbNare+8886oUaOGDRsWHx/vvMCvwAXuSAnvV56P9aBBgzIyMubOnTtp\n0qTBgwfnvV9CoXeY1NRUZaZNl/ZKlSoNHz58+PDhX375ZWJi4h3tGHmPoDpfRfx5nNDzUj0f\n1jsqr2j3FgBQBdcQAvAR7dq1K1++/Pbt22/evLl9+/aYmBjlyiLPKcsrF6c5/fHHH6dPnw4I\nCMh/gVzRat68uRAi/530lMkzCi07O/vSpUtBQUF5vzHLsrxnz5676bZt27bR0dFbt25VpkXJ\nP52Mopg2Svw5WC63pJs7d27btm0PHz7sdpXz589fvnzZ+dBkMiUmJrZs2fL69etnzpzxZIG7\nKbUE9ivPx/qtt97au3fviBEjJk6cOHDgwD179ijTqNxRJy5atWpVvnz5Wy2mzPRjMpmUh57v\nGF988YXLNp4+fVqn08XExHhYaiGG1ZPyimlvAYCSRyAE4CP8/Py6det24cKFhQsXZmZm5p2H\n00PKtUAvvfRSTk6Os/Gll16y2+0JCQnF/S9/5dWnT5+e9+bdn376aaVKlZR5RD2kfOdWjqII\nIfz9/cPCwjIyMpTzA4UQsizPnDlTmRjj+vXrhatWp9MNGDDgq6++SkpKqlat2qOPPup2MU82\nyqVgD/Xr108I8corrzjvFX7u3LmXX375yJEjtWvXzr/8N998o8x7mXdwLRbLmTNn9Hp9xYoV\nC1yg0KWW2H7l4VifP39+3Lhx0dHRyr0uXnnllYoVK44dO/b8+fOed5Jfhw4dhBD9+vXbu3dv\n3nZZlt9///358+eLP2fEFXeyt+/fvz9vwl+xYsXNmzcfffTRMmXKeFJq4Ya1wPI86RYASgtO\nGQWgmsuXLyu3XnCrUaNGs2fPvqMOe/fuvXTp0nnz5kmSlJCQcKf19OnT54MPPvjwww+bNGnS\nvn17g8HwxRdf7N+/v2bNmvPmzbvT3u5U9+7dt27dun79+vvvv79Hjx7BwcHffffdtm3bAgIC\nxo0b53k/tWvXliRp586dgwYNMhqNS5cu7d+//6uvvtqmTRslRO3YsSMtLS0pKenxxx/fsGFD\nVFRUId4rIcSAAQNmzZq1f//+GTNm3OpCPk82Kn/Bnrx6nz59Nm/evGPHjrp167Zv3z4zM3Pr\n1q0Wi2XFihXKtKsuGjZs2Lt373//+9+1a9du3759uXLlrl27tnPnzosXLz777LPlypUrV67c\n7Re4m1JLbL8qcKx79+49cOBAi8WyYcMGZdKU0NDQxYsX9+rVa8CAAfv375ckyZMdJjIy0uWl\nJ06ceOrUqXXr1rVr165atWqNGjUym81paWknT568cOGCTqebM2eO80YsnuwYdrtdCDFo0KD2\n7dt37tz53nvvPXXq1Hvvvefv7z9nzhwPtzchIaEQw1pgeQXuTkU7rABQvErwnocA8D/KXc5v\nz3nHbeVu7C+//HLeHn766SchRMOGDfM2OhyOqlWrCiEeeughl1d0ua252z5lWbbZbIsWLWrc\nuLHZbPb394+JiZk0aVJaWppzAU+KcbvMI488IoRQbgWuUKbFz3sH89zc3BUrVjzwwAPBwcF+\nfn6RkZF9+/bNu4qHb8W8efPKly/v7+/fuHFj5YVeeOGF6tWr+/v7R0VFDR8+/Nq1a7Is9+/f\nPzAwsHLlyt9++63nN6bP2/jYY4/pdLrz5887W1xuTO/JRuUv2O1muoygLMs2m23BggUNGjQI\nCAgIDAx8+OGHP/7441vVr1Ty5ptvtmzZsnz58nq9PiQk5KGHHlq1apXD4fBwgfyl5pf/xvRy\nYfer23N7Y/oCxzoxMVEI0atXL5felON7r732miedfPvtt7eq6qOPPkpISLjvvvsCAgJ0Ol2Z\nMmUaNGgwYsSI/KsUuGOMHj1aCPHee+8dOHDg4YcfDgoKCgwMfOSRRz799FPPt/fbb78t3LAW\nWJ4n3QJAqSDJt52pHAAAAADgq7iGEAAAAAA0ikAIAAAAABpFIAQAAAAAjSIQAgAAAIBGEQgB\nAAAAQKMIhAAAAACgUQRCAAAAANAoAiEAAAAAaBSBEAAAAAA0ikAIAAAAABpFIAQAAAAAjSIQ\nAgAAAIBGEQgBAAAAQKP81C6gNLFYLHa7Xe0qvIjRaAwICMjKyrLZbGrXgiKmDO7NmzdzcnLU\nrgVFjMH1YQyuD2NwfRiD65afn19wcLDaVWgCgfAOOByO3NxctavwIrIs63Q6WZZ5W3wPg+vD\nGFwfxuD6MAbXhzG4bul0nMlYQnijAQAAAECjCIQAAAAAoFEEQgAAAADQKAIhAAAAAGiU708q\nI8vymjVrPvvss8jIyJEjR4aFhaldEQAAAAB4BXUC4YULF5KSkk6dOiXLcrVq1fr06RMTE6M8\ntWvXri1btqSmpkZERPTt2zc2NtZl3YsXL65evfr06dN2u11Zt06dOkKIZ5555ty5c87FTCbT\npk2bhBDJycmXLl1aunTpwYMH33rrrUmTJgkhxowZ8+uvv7722mv33HOPc5WRI0d27Nixffv2\nxbz1AAAAAOAVVAiEdrt9ypQpDRs2nD9/vk6n27hx44wZM1atWhUQELB///6NGzeOGjUqOjr6\nyJEjK1asqFu3rtlsdq4ry/LMmTMbNGiwfPlyvV6/efPm6dOnr1y5Mjg4OCMjY+jQoS1atFCW\ndM5U+9NPP7Vq1Uqv18fFxa1bt87Zlb+//5tvvjlnzpyS3HYAAAAA8B4qXEOYmZkZHx+fmJgY\nERERHh7erVu3zMzMy5cvCyE2btzYr1+/pk2bVqxYMT4+fvny5XnToBDixo0bKSkpbdu2NZvN\n/v7+HTp0sFqtyroWi6Vy5crl/+Q8NbRatWqHDx+22+3JycnVq1d3dhUfH3/u3Ll9+/aV4KYD\nAAAAgBdR4QhhSEhI586dlZ8tFsu2bdsiIyOjoqJSU1NTUlKEEM8888zly5erVKkyePBg56mk\nznVjYmL27NkTERFhMBj27NlTqVKlqlWr2my27OzsI0eOrF271mKx3HfffX379o2IiBBCtG7d\n+uzZs4mJiVFRUSNHjnR2FRgYOGDAgFWrVsXGxoaEhLgt9fPPP587d67z4YwZM+rXr1/kb0jp\nJUmSECIwMNAlt8MHKINrNpsDAgLUrgVFjMH1YQyuD2NwfRiD61Zubq7aJWiFapPKOByObt26\n2Wy2evXqzZo1y2AwpKamCiH27ds3fvz4kJCQDRs2zJgxY9myZS5pbeLEiVOnTk1ISBBChIaG\nTp061Wg0pqenly1b1m63Dx8+XAixfv36SZMmLV26NDAwUKfTDRkyZMiQIflraNu27cGDB5cv\nXz5u3Di3RdrtdovF4nyYm5vrPBMVTpIkKb/I4HsYXB/G4PowBteHMbg+jMF14XA41C5BK1QL\nhDqdbvHixWlpaTt37pw8efKCBQuU9h49ekRGRgohBg4ceODAgWPHjrVp08a5lt1unzlzZkxM\nzJw5cwwGw65du6ZNm/b666+HhoauWbPGudj48eP79et3+PDhxx577PZljBgxYtSoUceOHWva\ntGn+Zx988MGPP/7Y+TA9PV1JrVCYTKagoKCMjIzs7Gy1a0ERUwY3MzPTarWqXQuKGIPrwxhc\nH8bg+jAG1y2DwXCrk/hQtNQ83hUZGVm/fv3x48enp6cnJycrV/0FBgYqz+r1+rCwsLS0tLyr\nnDx58uzZs4MHDw4JCTGbzV27dvX39z906JBLzwEBARUqVLh27VqBNYSHh/fs2XPp0qVWq5X/\nygAAAADQFBUC4fHjx4cOHeo8piRJkp+fnxAiLCwsNDT0hx9+UNpzcnKuXr1aqVKlvOvKsizL\nct4jyHa7XQhx/vz5N954Q/lZCGG1Wq9evVq5cmVP6uncubPZbH733Xf1ev1dbxwAAAAAlBoq\nBMIaNWpYrdZFixZduHAhJSXl7bfftlqtTZo00el0nTp12rBhw4kTJ65du/bWW2+ZTCblPoR7\n9+7dvn27ECImJiY0NHTVqlUZGRk5OTkffPBBZmZm06ZNw8LCjhw58sYbb6SkpFy6dGnhwoVB\nQUEtW7b0pB69Xj9q1Khdu3ZxOigAAAAATVHhGsKgoKBZs2atXr167NixkiRFR0dPmTJFOZrX\npUuXrKysV199NSMjo1atWrNnzzaZTEKIEydO3Lhxo1OnTmazeebMmUlJSYmJibm5udHR0dOm\nTQsPDxdCKH2OHj3aYDDUqlVr7ty5/v7+HpZUs2bNDh06KJkTAAAAADRCkmVZ7RpKjfT0dJvN\npnYVXkS5BtpisTCpjO9xzhjEBe6+h8H1YQyuD2NwfRiD6xaTypQYbqIAAAAAABpFIAQAAAAA\njSIQAgAAAIBGqXZjekA7+s1/78Lv15XLdf389O/PSijj8YxHuFNX0zKnr9n3y29/yEIODw2e\n/HRczcgKahdVWlmt4qkpb9scspCFEKJCaPCmKT3VLgoAABQljhACxeupF9/9NeW648/Jm+z2\n3PhJa1StyJdlZOUkzN343/NXsm32HFvu+SvX/7lw648Xr6pdV2nVYfIKW+7/0qAQ4mqapf3E\nd9QsCAAAFDUCIVCMUlNvWLKsQgjp7+2D//WeKvX4vGlJ+2z23LwtDllMe2efWvWUakMXfJh/\nEmprjv10ynUVqgEAAMWDQAgUo1e3HXHbfu7KjRKuRCPOXE7N33g1LavkK/EBv6Rcc9cs/2vt\ngZIuBQAAFBsCIVCM9JL7j5jMJ6+YSFL+NlnibquFoRPu3zeTTl/ClQAAgOLD11KgGD3f62HX\ns0WFEEI0qnpPideiCfWrVsrfGFmB29oWRuMa7vfSlwa2K+FKAABA8SEQAsWojL9/dLm/0ohy\nwEWSpAUj2qtVkm97MeFRs8mYt0Wv080b/IRa9ZRq/xraQZfv3xmhgaayZU1qlAMAAIoFgRAo\nXkmTuz/UoKokhJCEThJhwYEfLxisdlE+y2j0e396wiMNq5UrExBi9r//vns2Tu0VXi5Y7bpK\nq/0LhoQE+juPct9fM/yDWX1UrQgAABQxSZa5usZT6enpNptN7Sq8iMlkCgoKslgs2dnZateC\nIqYMbkZGhtVqVbsWFDEG14cxuD6MwfVhDK5bBoMhJISLPkoCRwgBAAAAQKP81C4AADTt6KkL\nL2/65A/LTUmSIsuHTOvbplp4qNpFAQAAreAIIQCo5vjPv018+6Nr6VkOh5yb6zj/e9rQVz9I\nvcGNEwEAQAkhEAKAauasPeByIbc91zEjab9a9QAAAK0hEAKAatIsN/M3/vp7WslXAgAAtIlA\nCACqkfLf6U8IPz9+MwMAgBLC1w4AUE2Vim7mj2lcI7LkKwEAANpEIAQA1cwf9oTx78cDQ4MC\nxvd6WK16AACA1hAIAUA15cqY35/Z5+EG1SqWDYwoVya+ZZ3N05/20/GbGQAAlCv7hNAAACAA\nSURBVBDuQwgAagoyGWf0b6t2FQAAQKP4PzQAAAAAaBRHCOGDZv/74/3HzgghCyEkSfRu23Bw\n+2ZqFwUAAAB4HY4Qwtcs2/b5/mO/KGlQCCHLYt3eb/Z99Yu6VQEAAABeiEAIX7Mp+WT+xnnr\nD5R8JQAAAICXIxDC18juGnMdbpsBAAAATSMQQhMkSe0KAAAAAO9DIISvMejdTJUUHGgq+UoA\nAAAAL0cghK9JeqGHS4skxOYpfVQpBgAAAPBmBEL4mvAQ885/9a8QEqTXSXqdFFUx9ONXhxgM\napcFAAAAeB/uQwgfZDYYNk3rpXYVAAAAgLfjCCEAAAAAaBSBEAAAAAA0ikAIAAAAABpFIAQA\nAAAAjSIQAgAAAIBGEQgBAAAAQKMIhAAAAACgUQRCAAAAANAoAiEAAAAAaBSBEAAAAAA0ikAI\nAAAAABpFIAQAAAAAjSIQAgAAAIBGEQgBAAAAQKMIhAAAAACgUQRCAAAAANAoAiEAAAAAaBSB\nEAAAAAA0ikAIAAAAABpFIAQAAAAAjSIQAgAAAIBGEQgBAAAAQKP81C4AKB0yrTl7jv3s76dr\n0/DegACj2uVo2s1s21c/XjL66RtVDzcavfGX2Gc/XJm5erfssE8b1KlVTEW1y3F1/vf0n4+f\nrRkZWiHYpHYtAABAZd74XQrwNhOW7z56+qKQhRBiwaZPH2l47/R+bdQuSqMWvn9o++EfZFkW\nQuh10tNt7+//RBO1i/qb1uNWyrkO5ecXl38oCfHxq0PULcnpwpUbY5ftuHo9U3lYOSzo1X8+\nGV4uWN2qAACAijhlFCjAkm1fHP3hf2lQkfzNmbX7jqtXkXa9l3xy22enlDQohMh1yEn/+frg\nibPqVpVX6zGrnGlQIQvRevxKterJy+EQz7zxoTMNCiFS/sgY9fo2FUsCAACqIxACBdh66Pv8\njWv3nSj5SuDubZeX7ziqQim3IIvcfC1CtjvcLlzCDhz/+XqG1aUx9UbWoZPn1CgHAAB4BQIh\nUAC73fUrvhAix+amEcUt05qTr026nnlThVI8JqldgNPp3665bT91/koJVwIAALwHgRAogE7n\n5mPip+ezowKjn17O12j2N6hQyi15TwB0FRFaxn071xACAKBhfKkFCtC4Rnj+xpZ1q5R8JWjd\n6N78eSu+VW0VSrkFnds86B0h8fFmNY0GvUujv9GvXWwtVeoBAADegEAIFGD+sA7lygTmbQkv\nV4ZZRlXxfI+Ha0VVcD6UhGgeE9XnscYqluRi/4LB0t/jnyTEgQVeMcuoyeg3vd9jBv1fmdDo\np5s94DE/P/4QAACgXZJzvj4UKD093WazqV2FFzGZTEFBQRaLJTs7W+1ait3er3/a/fmPOkl6\n6sHaD9evpnY5xU4Z3IyMDKvVdRoS1R07ffHAN2f9/XRtGteoW9Xr7vInhOj84js3suyyJMwB\n+h2zB6hdzt9Yc+wffn7619/Tq1YsE/9AjHfeyBGF5s2fXNwlBteHMbhuGQyGkJAQtavQBALh\nHSAQutBUINQa/jj5MAbXhzG4PozB9WEMrlsEwhLDmUIAAAAAoFEEQgAAAADQKAIhAAAAAGgU\ngRAAAAAANIr55QBozsefn3nl/QM24ahbpeKikfFqlwMAAKAajhAC0JZ2Y9+etWn/zVyHPVd8\nc+bKo2Pe/s1iUbsoAAAAdRAIAWjIPxd/YHO9146cMG2DOtUAAACojUAIQENO//pH/kap5OsA\nAADwDgRCAFrienhQCCFkIXZ9+V3J1wIAAKA6AiEALXF3NFASokNsvRIvBQAAQH3MMgqgiH10\n+Meluz+35cqt6lWd3PsRtcv5m+oRFX6+eNWlUZY4aRQAAGgURwgBFKX2E9+Ztzk5PTM7y5qz\n99iPrce+bbWqXVMeK8b8w0//t997kpD+/XwPteoBAABQF4EQQJEZtXirNccmxF/X6cmy3OHF\nt1UsKb+9Lw+a0r2N0V8v6aSYKuU+fnVweHiw2kUBAACog1NGARSZ735Vzsb82xmYssPNPC7q\nat3i3tYt7lW7CgAAAPVxhBBAkZG8LvoBAADgdgiEAIoMeRAAAKB0IRACKDIt6kTlb9Tp+T0D\nAADgpfiiBqDIzB38RFCAKW+LJKT9Lw9Sqx4AAADcHpPKAChK2+f0OZ1yfW7Sfyw3bU+3ub/z\nQ3XUrggAAAC3RCAEUMRqVS77zoTualcBAACAgnHKKAAAAABoFIEQAAAAADSKQAgAAAAAGuX7\n1xDKsrxmzZrPPvssMjJy5MiRYWFhalcEAAAAAF6heAPhhQsXkpKSTp06JctytWrV+vTpExMT\nozy1a9euLVu2pKamRkRE9O3bNzY21mXdixcvrl69+vTp03a7XVm3Tp2/TVe4f//+xYsXT548\nuUWLFkKIZ5555ty5c85nTSbTpk2bhBDJycmXLl1aunTpwYMH33rrrUmTJgkhxowZ8+uvv772\n2mv33HOPc5WRI0d27Nixffv2xfJeAMXmRpZ1w4GTl66mV78nrHtcA5PRx//Rc+jkuU9OnvPT\nSW2b1Ghc456CVyjI4e/OJ397TpJE28bVm9aKvPsOn1+y/ZszV2RZqntv5cUjOtxNV3+k31z8\nwaFLqZaqlco+1+3BQJPxLmv7I+PmkrUHL15Nj6wY8mznB4LN/nfZYdEau3zPtz//ppfklvWq\nTu3TRu1yAADwfcX4xdFut0+ZMqVhw4bz58/X6XQbN26cMWPGqlWrAgIC9u/fv3HjxlGjRkVH\nRx85cmTFihV169Y1m83OdWVZnjlzZoMGDZYvX67X6zdv3jx9+vSVK1cGBwcrC1y/fj0pKclo\n/Ou7UUZGxtChQ5VwKITQ6f53NuxPP/3UqlUrvV4fFxe3bt065/L+/v5vvvnmnDlziu8dAErA\nvq9/nvfv5FyHQwjxybdn1+07MW/oE/ffVwQxyTsNeuX9M7/9ofy8++iP9993z6vDO95Nh4kL\nt56+cFX5+aMvf6xXrdLro566mw5bj31blmXl529/ufTo2BUHFgwpXFfr9n/z9s6jys+//Jb6\n8YlfJvaKa9e0RqFr25z8/ZIPj8hCFkKcvnBl/1c/je76UHzLmEJ3WLQeHbtCyEIIYRfiwPEz\nySfO7l8wWO2iAADwccV4DWFmZmZ8fHxiYmJERER4eHi3bt0yMzMvX74shNi4cWO/fv2aNm1a\nsWLF+Pj45cuX502DQogbN26kpKS0bdvWbDb7+/t36NDBarUq6yqWLVsWFxeXdy2LxVK5cuXy\nf3KeGlqtWrXDhw/b7fbk5OTq1as7l4+Pjz937ty+ffuK7x0AiltGVs6//kyDihx77uSVH+Vp\n8Ckv/fugMw0qjv/826rdXxa6w1ffO+RMg4rvzv6+9MMvCt1h+4nvONPg/8ii3fhVhejqj4yb\nK3cd/VtPsvjX+uSc3NzC1WbJyl6y7bCSBp0Wbz5082ZO4TosWo+PX/X30oRDljtPWatSOQAA\naEUxHiEMCQnp3Lmz8rPFYtm2bVtkZGRUVFRqampKSooQ4plnnrl8+XKVKlUGDx7sPJXUuW5M\nTMyePXsiIiIMBsOePXsqVapUtWpV5dkjR4788ssvo0ePPnjwoNJis9mys7OPHDmydu1ai8Vy\n33339e3bNyIiQgjRunXrs2fPJiYmRkVFjRw50vkSgYGBAwYMWLVqVWxsbEhIiNtNOHPmzM6d\nO50PO3TokPcUU+j1eiGEv7+/n5+Pn6PotTYmf2fPF/6s2fbPf7z0WJOad9OzMrhGo1H5wUsc\n+u58/sZdX/w4qmtc4To8+M2Z/I3/+eqn53u3LlyH2TZb/kabPTcwMPBOu1qw+TOXaCmEcMjy\nh4d+6N+hWSFqW7bjaP4OZSGvT/5u1P89VIgOi1aO3c2/MdIzrYV467TMOz+5KBIMrg9jcKGu\nYv8e73A4unXrZrPZ6tWrN2vWLIPBkJqaKoTYt2/f+PHjQ0JCNmzYMGPGjGXLlrmksokTJ06d\nOjUhIUEIERoaOnXqVOUE0YyMjGXLlj333HMmk8m5cFZWVtmyZe12+/Dhw4UQ69evnzRp0tKl\nSwMDA3U63ZAhQ4YMcXPKVtu2bQ8ePLh8+fJx48a5Lf78+fNJSUnOh82bN897jBGKvCfuooRd\nS7/ptv3K9ayAgIC7799LBnfvlz+s2fNl1s3s7Bw3cetmjr3QG5tts+dvtObYCt1h/sSlKESH\nV9Mz3bZfuGYpXHkpae47TPkjo0j2lrvm9r2TvaO2UsZLPrkoDgyuD2NwXdjtbv5GozgUeyDU\n6XSLFy9OS0vbuXPn5MmTFyxYoLT36NEjMjJSCDFw4MADBw4cO3asTZu/5g+w2+0zZ86MiYmZ\nM2eOwWDYtWvXtGnTXn/99dDQ0JUrVzZu3LhRo0Z5XyUkJGTNmjXOh+PHj+/Xr9/hw4cfe+yx\n25c3YsSIUaNGHTt2rGnTpvmfbdiw4ZIlS5wPo6Ki0tPT7/w98FlGozEgICArK8vm7qgISkC1\nSmXctsdEhN3lvqoM7s2bN3NyVD6fcPD89374+1mdLsKCTIXe2OAA/1RblktjSGDhOxSS5DYU\nFqLDWveU++J7N0dEm9xXuXDl1a5SIfn4z/nb61Wr6B2/2aT8mVAu1FunZd7zyUWRY3B9GIPr\nll6vDwoKUrsKTSiJM/0iIyMjIyPr1q3bu3fv5OTk5s2bCyGcZwHp9fqwsLC0tLS8q5w8efLs\n2bPz5s1TDgN27dp19+7dhw4dioqK+vrrr994443bv2JAQECFChWuXbtWYG3h4eE9e/ZcunTp\nm2++KUmSy7NhYWHNmv11alZ6ejrJJy/lxIbc3FzeFrU8EVtr+fYv0rOy8zZGlCsTE1XuLgfF\nSwZ37b7jt0+DQohhT7UodJED2jddsPETZwqRhZCEGPJks0J3WC283Nnfrin9OEVUDC5Eh33b\nNVl/4ITt7ydSmv0NjzSsWrjyejxcb+WOL1w69Dfon3ogxhs+wtGVQ35Nue7SGFvjHm+orRTx\nkk8uigOD68MYXKirGCeVOX78+NChQ7Oz//dVVZIk5UqzsLCw0NDQH374QWnPycm5evVqpUqV\n8q4ry7Isy448F0cpR4337t2bmZmZmJiYkJCQkJCQnp6+cOHCuXPnnj9//o033nAeWbZarVev\nXq1cubIndXbu3NlsNr/77ruct41SR6cTrz/7VMWyf11kVaVS6GujOqlYUtH64JPvbvOsn043\nuGPsg/WqFLr/js1r9WrdUKf7X3zT66R+7Rq3ub/wZ4aver5zoMmQNw0GmvzXTuxZiK70erFo\neHzem4gEBRhXjO1c6NqEEEue/YfZ3+B8aDYZ3hrd5W46LEJJ47uZzca//i8niXIhAfP/eVdT\nyAIAgAIV4xHCGjVqWK3WRYsW9e7d22AwbN++3Wq1NmnSRKfTderUacOGDcqRw/Xr15tMJuU+\nhHv37rVarZ06dYqJiQkNDV21alX//v2NRuOOHTsyMzObNm0aFxc3YMAA50s899xzffv2bd68\nuSRJR44csdvtPXv2zM3NXbNmTVBQUMuWLT2pU6/Xjxo1asKECRyVRmkUVaHsxqm9r6Zl/nDh\nWt2qFcLKmAtep/Sw5rifUXPe4CfMJr+6VcN1d/1PrSFPNhvQoenJMyk6IdW/twg63PFSf6tV\nzF23167LHdcrrmyeq53vVJ2q5XfPG3D+8vUvf7zUsm7kPeXdT3/lufsiyu2c2z8lPevEj1ea\nxIRXCPaumxDunN3vRnb2ok2HAgy6Zzo/6u9d1QEA4JuKMRAGBQXNmjVr9erVY8eOlSQpOjp6\nypQpylG7Ll26ZGVlvfrqqxkZGbVq1Zo9e7ZyauiJEydu3LjRqVMns9k8c+bMpKSkxMTE3Nzc\n6OjoadOmhYeHCyGctyIUQkiSFBwcXKZMGSGE8lqjR482GAy1atWaO3euv8ffJmrWrNmhQ4ft\n27cX/bsAlIgKoYEVQn1wMsbQMgE3r7k5f6Z5nagifBU/na5o79xoMokZgwq4gNlzVcLLVgkv\nW1S9CSGqVgqrVz06IyPDarUWYbdFooy/P/ejBwCgJEmut8zCrXENoQuTyRQUFGSxWJwnBsNn\nKIOremb47/nfR762zeW3VPPaUfOGPKFSRb7ASwYXReXHi1ffS/4uNT2rZlT54fGtwsLKMrg+\niU+uD2Nw3TIYDLe6MxyKFrePA+C96lSpNPTJFit2HlWuKJaEqB1diTQIOC398ItNn5xU5rY9\n/vNvWw59v3H6gMqh3KsDAOApAiEAr9bz0fo9H63/w69XL1y7/kjDe41M/gT86ftzV5xpUJFj\nyx38r/U75g1UsSoAQOlCIARQCsREV4iJrqB2FYB32XTw2/z3vUy9kXkp9Xq5wMLPZgQA0JRi\nvO0EAAAoPulZ7i83SknNKOFKAAClF4EQAIBS6b57yuVvlISoHV2x5IsBAJRSBEIAAEqlwR1i\nTUaDS2Nc4xpmk1GVegAApRGBEACAUslk9Fs8slOFkEAhCSGEJKQH61dd9ExXtesCAJQmTCoD\nAEBpVTOy3KZpvbNy7L+mpMVEVzCZTDr+0wsAuBP83QAAoHQzG/2YhhcAUDgcISx97HZ7t8Sp\nZ3+9KEm6RrWrrV40Te2KAAAAAJRKHCEsZf64caP+Y31O/XzGmpNzM9t65MSpOm0S7Ha72nUB\nAAAAKH0IhKXM472fk/9+G2JHbm77vmPUqgcAAABA6UUgLGUsGVn5Gy9dvlrylQAAAAAo7QiE\npY7srslNIwAAAADcHoGwlNHrpPyNRiP3IAYAAABwxwiEpczIgW7uOPzy5OElXwkAAACA0o5A\nWMr88+n/6/WPtpL0v+OEkqR7dnC3do+0ULcqAAAAAKUR9yEsfaaNHjxt9OAMq9Xk5+fnxwgC\nAAAAKCTiRGkVZDKpXQIAAACA0o1TRgEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFJPK\nAAXLyrI9/a8N1zOskiTKlTG/M6m72WBQuyigMMa9tevYj5ckIWQhHq5Tdcagx9SuCAAAqIkj\nhEABrl7P7PjiO2kWqywLh0NcvZ7VccI76elZatcF3LEnJ685dvqSkIUsCyGLT74/99QLa9Qu\nCgAAqIlACBTg6Zc25W/sPW9jyVcC3I09x05nWnNcGi03sz///qIq9QAAAG9AIAQKkGO3u7TI\nQmRluzYCXm7ph0eEkPO3/2vjwRKvBQAAeAsCIVAAyV2L5KYZ8Go5djdpUAiRbbeVcCUAAMB7\nEAiBArjNfjodgRClzIP1qrptf7zxfSVbCAAA8CIEQqAA/3yqef7DhBN7xalRC1B4LyQ86uaA\ntySe7fqQGuUAAACvQCAECtA1rn6PuHrO79GSTgzqGNu2SXU1awIKZcvUp/U6nTMV6iX9ljl9\nVK0IAACoTJJl91eVIL/09HSbjYtt/mIymYKCgiwWS3Z2ttq1oIgpg5uRkWG1WtWuBUWMwfVh\nDK4PY3B9GIPrlsFgCAkJUbsKTeAIIQAAAABolJ/aBQAASrecnNznlu765bershDlQ4JeHvr4\nPeX5ny4AAKUDRwgBAIWXk5v75Ivv/Pd8SrYtN8eW+9u19ISX3jt96YradQEAAI8QCAEAhTf6\nzR02u+PvbfLYN3apUw0AALhDBEIAQOGduZSavzErh/m3AAAoHQiEAIDCcztRNdNXAwBQWhAI\nAQCFV7lscP5Go0Ff8pUAAIBCIBACAApv4fAnhST9rUkSE7o/olI5AADgzhAIAQCFFxYSsOr5\nzqHBZkkSkpDM/oYpvVu3blJd7boAAIBHuA8hAOCuVAsv98GMBLWrAAAAhcERQgAAAADQKAIh\nAAAAAGgUgRAAAAAANIpACAAAAAAaRSAEAAAAAI0iEAIAAACARhEIAQAAAECjCIQAAAAAoFEE\nQgAAAADQKAIhAAAAAGgUgRAAAAAANIpACAAAAAAaRSAEAAAAAI3yU7sAAF7nza1HPjr2Y5bV\nbg4wPvVA7cEdmqpdEbxX9+nrrt3IkoUQQoQFm9+fkaByQRq2ZPvnez7/MSvbFmgytG9eK7FT\nc7UrAgCUAhwhBPA345bt3vzJd5asnFyHw5JpXbfv+JTVe9UuCl6q3YTVV/9Mg0KIPyxZjz2/\nUs2CNGzCij3vHThpuZmd63DcyMreeODbyW9/pHZRAIBSgEAI4C8//Hr12I8XXRoPnTx34ep1\nVeqBN9v5xQ82m92l0e5wrNpzTJV6tOznS6lHT11waTzy31/PXE5TpR4AQClCIATwl4PfnHHf\nfsJ9O7Rs3d5v3bbvOHK6hCvBwW/Oum8/8XMJVwIAKHUIhAD+4qfXu23X6/hdAVcGP8l9u569\npaQZ/dy/5wY/959oAACc+LMN4C+PN60pJNdv+ZKQHo+tqUo98GaJ8c3ctvdv36SEK8FjTWrk\nT+eSJLVrWkOFagAApQqBEMBfoiqW6dC8lhBCCOdEIaLbI/XLlTGrVRK81gO1q4QGmfI0yEKI\noABDe/59UOLCywU/1bKOS2N8yzqVQoNVqQcAUIrop0+frnYNpUZ2drbD4VC7Ci/i5+dnNBpz\ncnJyc3PVrgVFplXdKveUC75wNV3IonK54NFdHuwWV0/tolCUnJ9cu911Spg71ePRhqcuXrl0\n7YYkhCRJ999XOWlijyIpEneqRZ3oyAohF35Pl4UIDwse2aVlr9YN1S4KRakIP7nwNgyuW3q9\n3mQyFbwc7poky3LBS0EIIUR6errNZlO7Ci9iMpmCgoIsFkt2drbataCIKYObkZFhtVrVrgVF\njMH1YQyuD2NwfRiD65bBYAgJCVG7Ck3glFEAAAAA0CgCIQAAAABoFIEQAAAAADSKQAgAAAAA\nGkUgBAAAAACNIhACAAAAgEYRCAEAAABAowiEAAAAAKBRBEIAAAAA0CgCIQAAAABoFIEQAAB4\nNYdD2B0OtasAAN/kp3YBAAAA7h09dWHe+uS0TKuQ5WCz/7OdW7Zpcp/aRQGATyEQAgAAb3Tm\nctrEt/fI8v8eWrKyZ//7YEiQqWmtSFXrAgCfwimjAADAG81Z+7EzDf6PLP9rQ7I61QCAjyIQ\nAgAAb/R7miV/43WLteQrAQAfRiAEAADeyE+vz9cm6/WSCqUAgO8iEAIAAG/UrHZUvjapbrXK\nKpQCAL6LQAgAALzRxJ5xFcoG5m0JNhvnDGynVj0A4JOYZRQAAHgjnU5smtr7nT1fHfnvr7kO\nR+NaEYkdW+j4VzYAFCkCIQAA8F79n2jS/4kmalcBAD6L/7MBAAAAgEYRCAEAAABAowiEAAAA\nAKBRvn8NoSzLa9as+eyzzyIjI0eOHBkWFqZ2RQAAAADgFdQ5QnjhwoXZs2cnJCT07t37hRde\n+OGHH5xP7dq1a8iQIV26dBk1atSXX36Zf92LFy/OmjXr6aef7tmz56RJk/773/+6LLB///6n\nnnrq888/Vx4mJydfunRp6dKlrVq1euutt5TGMWPGdO3a9bfffsu74siRI3fv3l2U24k7lJGT\n8+8D32Tk5KhdCAAAAKAJKgRCu90+ZcqUwMDA+fPnL1iwoEKFCjNmzLh586YQYv/+/Rs3bhw2\nbNiyZcvatm27YsWKrKysvOvKsjxz5szQ0NDly5cnJSXVq1dv+vTpFovFucD169eTkpKMRqOz\n5aeffmrVqpVer4+Li/vpp5+c7f7+/m+++Wbxby48cvri763Hruw0MWnF9qOdJia1HvP26Yu/\nq10UAAAA4ONUCISZmZnx8fGJiYkRERHh4eHdunXLzMy8fPmyEGLjxo39+vVr2rRpxYoV4+Pj\nly9fbjab865748aNlJSUtm3bms1mf3//Dh06WK1WZV3FsmXL4uLi8q5VrVq1w4cP2+325OTk\n6tWrO9vj4+PPnTu3b9++4t9iFCxx4TZZdjgfykL+56vbVKwHAAAA0AIVriEMCQnp3Lmz8rPF\nYtm2bVtkZGRUVFRqampKSooQ4plnnrl8+XKVKlUGDx4cExPjsm5MTMyePXsiIiIMBsOePXsq\nVapUtWpV5dkjR4788ssvo0ePPnjwoHOV1q1bnz17NjExMSoqauTIkc72wMDAAQMGrFq1KjY2\nNiQkxG2pFovl4sWLzoflypXLe+wROp1OCKHX6/387mpHemnNPiG7NspCzFmbPK1/m7vpGYWm\nDK5Op7vLwYUXYnB9GIPrwxhcH8bguqW8LSgBqu12DoejW7duNputXr16s2bNMhgMqampQoh9\n+/aNHz8+JCRkw4YNM2bMWLZsmUtamzhx4tSpUxMSEoQQoaGhU6dOVUJaRkbGsmXLnnvuOZPJ\nlHd5nU43ZMiQIUOG5K+hbdu2Bw8eXL58+bhx49wWeezYsbxPLVmypFmzZne96b7GbDa7HMi9\nU1//ctlt+/FfLpctW/ZuesZduvvBhddicH0Yg+vDGFwfxuC6sNvtapegFaoFQp1Ot3jx4rS0\ntJ07d06ePHnBggVKe48ePSIjI4UQAwcOPHDgwLFjx9q0+esYkd1unzlzZkxMzJw5cwwGw65d\nu6ZNm/b666+HhoauXLmycePGjRo1uqMyRowYMWrUqGPHjjVt2jT/sxEREV26dHE+DA0NtVqt\nhdlaH6XX6w0Gg81my83NvZt+QoJMv6dluGs38IarpagGF16IwfVhDK4PY3B9GIN7KxwyLRlq\nvsuRkZGRkZF169bt3bt3cnJy8+bNhRCBgYHKs3q9PiwsLC0tLe8qJ0+ePHv27Lx585TDgF27\ndt29e/ehQ4eioqK+/vrrN954405rCA8P79mz59KlS998801JklyerVmz5uTJk50P09PTMzLc\n5BbNMplMBoPBarVmZ2ffTT8v9Hus3+wNwuW0UUlMG/A4b7halMHNzs4mk/seBteHMbg+jMH1\nYQyuWwaDweW8PxQTFc7NPX78+NChQ50RQpIkJf2HhYWFhoY6b0GRk5Nz9erVSpUq5V1XlmVZ\nlh2Ov2YfUY4m7927NzMzMzExMSEhISEhIT09feHChXPnzvWkns6dO5vNIC/W4gAAIABJREFU\n5nfffVev1xfJBuJORYcFNagW7tLYoOo90WFBqtQDAAAAaIQKRwhr1KhhtVoXLVrUu3dvg8Gw\nfft2q9XapEkTnU7XqVOnDRs2KEcO169fbzKZYmNjhRB79+61Wq2dOnWKiYkJDQ1dtWpV//79\njUbjjh07MjMzmzZtGhcXN2DAAOdLPPfcc3379lUOORZIr9ePGjVqwoQJQUHED9UsHtXxyKnz\nc9YcvJljCzAaXugb90DtKmoXBQAAAPg4FQJhUFDQrFmzVq9ePXbsWEmSoqOjp0yZUrlyZSFE\nly5dsrKyXn311YyMjFq1as2ePVs5UnzixIkbN2506tTJbDbPnDkzKSkpMTExNzc3Ojp62rRp\n4eHhQojg4GDnS0iSFBwcXKZMGQ9LqlmzZocOHbZv314MmwtPPVC7yo65/dSuAgAAANAQSZbz\nzfePW0hPT7fZbGpXURg5dvuxE6fCyobE3BddhN2aTKagoCCLxXKX1xDCCymDm5GRwfUMvofB\n9WEMrg9jcH0Yg+uWwWC41Z3hULSYusf3TZi7ZNt/PlWSv7+/YcZzg//xxCNqFwUAAABAfdzw\n0ce9vOzfH370ifM4cHa2bfL8Zf/98ZyqRQEAAADwCgRCH7f2/d0uLQ6HPHPxSlWKAQAAAOBV\nCIQ+Ltvu5qLHy79fK/lKAAAAAHgbAqGP00tuhjg4KLDkKwEAAADgbQiEPq5Jg5j8jf26diz5\nSgAAAAB4GwKhj1s+f2LlCmF5GqT2cS26PfmoagUBAAAA8BrcdsLHmYzGg+8t+WD3wU+/OBEY\nGPB/HVrfX7eG2kUBAAAA8AoEQk3o0j6uS/s4tasAAAAA4F04ZRQAAAAANIpACAAAAAAaRSAE\nAAAAAI0iEAIAAACARhEIAQAAAECjCIQAAAAAoFEEQgAAAADQKAIhAAAAAGgUgRAAAAAANIpA\nCAAAAAAaRSAEAAAAAI0iEAIAAACARhEIAQAAAECjCIQAAAAAoFF+ahcAAL7m6KnfNiUfl2Wp\nx6P1m8VEqV0OAADALREIAaAo9Zix/kp6hvLz1z9dCg0O+GDG0+qWBAAAcCucMgoARWbC8l3O\nNKhIs9x89s0datUDAABwewRCACgyx07/lr/x5JmUkq8EAADAEwRCACgyDlnO3yi7awQAAPAG\nBEIAKDI6neSm0U0bAACAVyAQAkCRebBelf9n7/7jmrzu/o+fEAKBoBhEASFYf6PV2h/Y6tpN\nq7i7SilV662zLdT6o66CWq1O3TAF7a1rZ+vPqqgwtLXQ1rrV6VjBFb5z6rzV6npXsTpFQWEV\nUCBoIIHr+0e2jGJAhIRArtfzwR9cJ9d18knOVn17ruucuxvDBwS3fSUAAADNQSAEALtJfGVs\nSLfO/zlWiKCuPr+ePd55FQEAADSFbScAwJ52L5ty4XpZetbXZqnuZ2OGhum6O7siAACARhEI\nAcDO+vXwS4gd4+wqAAAA7o1bRgEAAABApgiEAAAAACBTBEIAAAAAkCkCIQAAAADIFIEQAAAA\nAGSKQAgAAAAAMkUgBAAAAACZIhACAAAAgEwRCAEAAACgUadPn1YoFM6uwlEIhAAAAADQqL/8\n5S/OLsGBCIQAAAAAZKGoqGjWrFk9e/ZUq9WBgYGTJk3Ky8tr+qVnnnlm3rx5QgiFQhEeHu7M\n6h3D3dkFAAAAAEBbmDhxYn5+/qpVq3r37l1UVLRmzZqRI0devnzZ29u7sZc2bty4ePHi3//+\n9//7v/+r0Wic/Qnsj0AIAAAAwPVVVFQcO3Zs6dKlM2bMsLQ88cQTn3zyya1bt8xmc2Mv9evX\nz9/fXwjhktODgltGAQAAAMiBl5dX165dP/7440OHDtXV1Qkh+vTps2zZsh49ejTxkrOrdjgC\nIQAAAADXp1Kpfv/737u5uUVERHTv3v2FF17Ys2eP2Wxu+iWXRyAEAAAAIAtPPvnkhQsXDh06\nNH369HPnzr344osjRoy4c+dO0y+5NgIhAAAAALlQKpWjR49+9913v/322w8++ODEiROffPLJ\nPV9yYQRCAAAAAK7v5MmTU6dO/f77760tP/3pT4UQN27caOIlIYRlV3pXvYOUVUYBAAAAuL7g\n4OCDBw+eO3du/vz5oaGhpaWlGzZs6Ny584QJEzQaTWMvCSEsS8v8z//8z4MPPjhp0iRnfw47\nY4YQAAAAgOsLDAw8fPiwZfnQ8ePHL1y4MCAgICcnp0+fPk28JISYNWvWI488smrVql/+8pfO\n/hD2p5Akydk1dBjl5eUmk8nZVbQjarXax8ensrKyurra2bXAziyDazAYjEajs2uBnTG4LozB\ndWEMrgtjcG1SqVS+vr7OrkIWmCEEAAAAAJkiEAIAAACATBEIAQAAAECmCIQAAAAAIFMEQgAA\nAACQKQIhAAAAAMgUgRAAAAAAZIpACAAAAAAyRSAEAAAAAJlyd3YBAAAAAOBw392sq6m1c59K\nhRjYtWPPsREIAQAAALi+TSdr/lkl2bdPb5Uidbzavn22MQIhAAAAANcnSXWSZOdAKEkK+3bY\n9giEAAAAAGRAEsLegVDYu7+2RyAEAAAA4PokSXLADKF9+3MCAiEAAAAAGZDqhFRn7z65ZRQA\nAAAA2j1mCG0iEAIAAACQAUniGcK7EQgBAAAAuD5JkiR73zLKKqMAAAAA0CEwQ2gDgRAAAACA\n6+MZQpsIhAAAAABkgFVGbSEQAgAAAHB9zBDaRCAEAAAA4Pokqc4Bi8q42bfDtkcgBAAAACAD\nkgNm9Dr+FCGBEAAAAIAMOOSWUQIhAAAAAHQADlhUpuMjEAIAAABwfY5ZVIYZQgAAAABo/yRH\nbExPIAQAAACA9k+SHLAPoZ37a3sEQgAAAACuzxG3jDJDCAAAAAAdAbeM2kIgBAAAAOD6WFTG\nJgIhAAAAABmQHLDtRIfPgwRCAAAAADIgCWYIbXBzdgEOJ0lSWlra7Nmzk5KSysrKnF0OAAAA\nACeQJEmS6uz+c8/3PX/+/PDhw93dfzAVd/369WnTpgUEBHTu3HnkyJHHjx9vcFVOTo7iLps2\nbWqsvcVfi2NnCAsKCtLS0s6dOydJUq9evV5++eWwsDDLSwcPHty3b19paWlwcHBMTMywYcMa\nXFtYWJiamnr+/Hmz2Wy5dtCgQU30OW/evPz8fOvlarX6k08+EULk5uZeu3Zty5YtOTk527Zt\nW7ZsmRBi4cKFV69e3bBhQ48ePayXxMXFRUZGjhs3zqHfCQAAAAAncMiiMvd4PSMj44033hg7\nduyJEyfqt0dHR6vV6szMzE6dOun1+sjIyPz8fI1GYz1hxIgRBQUF1sP8/Pxx48aNHj26T58+\nNttb/AkcGAjNZnNCQsLQoUPfeecdNze3jIyMxMTElJQULy+vQ4cOZWRkxMfHh4aGHj16dPv2\n7Q8++KC3t7f1WkmSkpKSHnrooeTkZKVS+dlnn7311ls7d+708vJqrE+DwTB79uzhw4dbenBz\n+9fk54ULF5588kmlUjlq1KiPPvrI+haenp6bN29+++23HfcNAAAAAGgnmjmhd599Kpo+obq6\n+tixY6dOnaqfRMrKykJDQ1etWjVw4EAhxOrVq/fs2XP27Nn6k2Senp4hISHWw+nTpy9atMgy\nQ9ZYe8s48JbRqqqq6OjoOXPmBAcHBwUFTZ48uaqqqqioSAiRkZERGxsbHh7evXv36Ojo5OTk\n+mlQCFFRUVFcXBwREeHt7e3p6Tl+/Hij0VhUVNREn5WVlYGBgf7/5ufnZ+mqV69eR44cMZvN\nubm5ffr0sb5FdHR0fn5+dna2474BAAAAAO2FZYbQ7j9NiomJCQ0NbdDo5+e3d+9eSxoUQly7\ndk2pVOp0usY6SU9Pv3jx4vLly5vZfl8cOEPo6+s7YcIEy++VlZVffPFFSEiITqcrLS0tLi4W\nQsybN6+oqKhnz54zZ8603kpqvTYsLCwzMzM4OFilUmVmZgYEBDzwwAMeHh42+zSZTNXV1UeP\nHv3www8rKyv79u0bExMTHBwshBg9evTly5fnzJmj0+ni4uKsb6HRaKZPn56SkjJs2DBfX1+b\nH+HYsWOrV6+2HiYmJg4ZMsSuX1LHplAohBAajaZBnocLsAyut7e3l5eXs2uBnTG4LozBdWEM\nrgtjcG2qra21e58Th/q7u/1nQs9QXZtytOh+O5k2LMBfo7Ie1rb6FtSysrIZM2YsWrQoMDDQ\n5gm1tbV6vT4hIcHDw6M57ffL4auM1tXVTZ482WQyDR48eOXKlSqVqrS0VAiRnZ29ZMkSX1/f\n9PT0xMTErVu3NkhlS5cuXbFixYsvviiE0Gq1K1assH7Uu/ssLy/v0qWL2Wx+/fXXhRAff/zx\nsmXLtmzZotFo3NzcZs2aNWvWrLtri4iIyMnJSU5OXrx4sc3izWZzZWWl9bC2ttZ6JyqsLE+y\nOrsKOASD68IYXBfG4LowBteFMbgN1NXZe38IIT7/+vuiippWdvLR8eL6hxoP5bODu7a4t7y8\nvKioqLFjx65Zs6axcz799NOqqqqYmJhmtt+v+wiERqPxm2++KSws/PGPf+zv7282mxsslWOT\nm5vb+vXrb968eeDAgeXLl69du9bSPmXKFMvNr6+++upXX3114sSJMWPGWK8ym81JSUlhYWFv\nv/22SqU6ePCgXq/fuHGjVqu12aevr++uXbusly9ZsiQ2NvbIkSNjx45tury5c+fGx8efOHEi\nPDz87lefeuqpP//5z9bD8vJyS5qFhVqt9vHxMRgM1dXVzq4FdmYZ3KqqKqPR6OxaYGcMrgtj\ncF0Yg+vCGFybVCpVYzfxtYLkgH0IWx7jDx06NGXKFL1eHx8f38Rpu3fvnjRp0t3Jq7H2+9Xc\n+a61a9d279798ccfnzhx4sWLF4UQer1++vTpZrP5nteGhIQMGTJkyZIl5eXlubm5lqf7rEvo\nKJVKPz+/mzdv1r/km2++uXz58syZM319fb29vV944QVPT8/Dhw831meDd/Ty8urWrVtJSck9\nawsKCpo6deqWLVuMRiP/KgMAAAC4KskxWlbM4cOHJ0+evHv37qbT4K1bt7KysqKioprZ3gLN\nCoTbt29/8803n3766a1bt1obBwwY8OGHH77//vuNXfX111/Pnj3bOnekUCgs+dXPz0+r1ebl\n5Vnaa2pqbty4ERAQUP9ay5dbf6bYkjwb6/PKlSubNm2yplOj0Xjjxo3GbsNtYMKECd7e3rt3\n71Yqlc05HwAAAEDH45AVZe4RCIuLiwsLCy23GRYWFhYWFhoMhjt37sTGxi5YsGDIkCGF/1ZV\nVSWE2Llz5/r1662Xnzx50mQy9evXr0G3jbW3QLMC4aZNm+bMmfP73/8+NjbW2hgTE7N48eId\nO3Y0dlW/fv2MRuO6desKCgqKi4t37NhhNBofe+wxNze3qKio9PT006dPl5SUbNu2Ta1WW5ZY\nzcrK2r9/vxAiLCxMq9WmpKQYDIaamprPP/+8qqoqPDy8sT79/PyOHj26adOm4uLia9euvf/+\n+z4+Pj/60Y+a8+mUSmV8fPzBgwe5HRQAAABwVY7Ylf6e+1gMHz5cp9PNnDmztrZWp9PpdLod\nO3YcOXLk0qVLer1eV09qaqqoF4gsioqKFApFUFBQg24ba28BRXNmOb28vPbv3x8REWE0Gr28\nvI4ePWrZ7u/LL7989tlna2oafTTzypUrqampZ8+eVSgUoaGhL7300tChQ4UQdXV1H374YXZ2\ntsFgGDBgwOuvv25ZZfXdd9+tqKhYuXKl5dq0tLTvvvuutrbWcq1lhc/G+rx06VJqauqFCxdU\nKtWAAQNmzZrVYNaxvoULF44ZMyYyMtLasn379v379//85z9vYmP68vJyk8l0z69LPiy3vFdW\nVvIMoeuxPiDK8wyuh8F1YQyuC2NwXRiDa5MjniGMTTldVG7nv7VqPJX7XrexFkkH0qxAGBAQ\nsHPnzmeffbZBIPz0009nzpxZXl7u+DrbBQJhAwRCF8YfTi6MwXVhDK4LY3BdGINrkyMCYczO\nrx0RCH83d9i9z2vHmnXL6EMPPfSb3/zmzp079RvLysqSkpIsyRAAAAAA2jPLKiVtfMto+9es\nVUp/+ctfRkREPPTQQ5Z7LLdv375169Z9+/bduXOn/jIzAAAAANBO/WsZGHv32cE1a4Zw1KhR\nf/rTnzp16mRZ8SYlJSUtLS0sLCwrK+vJJ590cIUAAAAA0Hp1QrL7T4cPhM3dx3DMmDGnTp36\n/vvvr1+/LoTo2bOnZY94AAAAAGj/WrNtYBN92rfDtnd/G9t37969e/fuDioFAAAAAByFW0Zt\naSoQ+vj43PN6k8nECpMAAAAA2jlJsv+EnovPED777LPW30+fPn3p0qXw8PAePXrU1tbm5+ef\nOXPm0UcfHTFihOOLBAAAAIDWsTz1Z/c+O7imAmF6errll88+++zbb7+9cuVKUFCQ9dXz588/\n//zzP/3pTx1bIAAAAAC0Gs8Q2tSsVUYTExNXrFhRPw0KIQYMGDB//vyEhATHFAYAAAAA9mN5\nhtDuPx1csxaV+e677/z8/O5u9/f3z8vLs3dJAAAAAGB3zBDa0KwZQn9//9TU1AaNkiR99tln\nNoMiAAAAALQv9t+E0AEPJba5Zs0Qzpo1KzEx8e9///vTTz/drVs3IURxcfGf//znc+fOLV26\n1MEVAgAAAEBrscqoTc0KhHq93tvbe926dRs2bLA2+vv7JyQk6PV6h9UGAAAAAHbCKqO2NCsQ\nKhSKJUuWLF68uKCgoLi4WJKkbt26PfDAA25uzbrjFAAAAACcyyGrjAp5zBBaKBSK0NDQ0NBQ\nx1UDAAAAAA7hiEVBZXLLqL+/f2Mv1dTUVFRU2K8eAAAAALA/SaqT7H2Hpws8Q9isez6fukuf\nPn1u374dFBQUExPj6BIBAAAAoLUkx2xFeC/nz58fPny4u/sPpuKuX78+bdq0gICAzp07jxw5\n8vjx4w2uysnJUdxl06ZNllc/+OCD3r17e3p6PvTQQ3/4wx9a8600a4bwd7/73d2NxcXFU6ZM\nGTduXGveHgAAAADagOSAfQjvGQgzMjLeeOONsWPHnjhxon57dHS0Wq3OzMzs1KmTXq+PjIzM\nz8/XaDTWE0aMGFFQUGA9zM/PHzdu3OjRo4UQv/3tb1euXLljx47Bgwd//vnnCxYs+MlPftK5\nc+eWfYKWrwoTGBi4du1aVhkFAAAA0AE4YnrwXoGwurr62LFjEyZMqN9YVlYWGhqanJz8yCOP\n9O3bd/Xq1SUlJWfPnq1/jqenZ0g9iYmJixYtGjRokBBi1apVa9asiYyM7Nmz5xtvvHHx4sUW\np0HRmkAohAgJCWlQNwAAAAC0Q5ZnCO3+0/SbxsTE3L0qp5+f3969ewcOHGg5vHbtmlKp1Ol0\njXWSnp5+8eLF5cuXW07+xz/+IYQYOnSoRqN54oknjh492pqvpeWBUJKklJSUrl27tubtAQAA\nAKBNSI75aZWysrIZM2YsWrQoMDDQ5gm1tbV6vT4hIcHDw0MIUVhYKIRITU3NyMgoKCgYMWLE\n+PHjb9y40eICmvUM4cMPP3x3WcXFxSUlJW+++WaL3xsAAAAA2sbUp/q519tHvdJo+uCPf7/f\nTl4dM6hbZy/rYW1dqwJhXl5eVFTU2LFj16xZ09g5n376aVVVVYO1PH/1q1+FhYUJId59991d\nu3YdOHDglVdeaVkN97EPYX0qleqhhx6Kjo6eM2dOy3oAAAAAgDbz8f87f72sqpWdpGR/W//Q\nR62aOLxPy7o6dOjQlClT9Hp9fHx8E6ft3r170qRJ1kVKe/ToIYTo0qWL5VClUvXo0aOoqKhl\nNYhmBsLTp0+3+A0AAAAAwPkkB+wj39IODx8+PHny5I8++qjpXRtu3bqVlZX1xhtvWFt69OgR\nFBR09OjR8PBwIcSdO3euXr3aq1evlpUhmvkMYXh4+Llz5+5u37t3r2WhGwAA0KHdMjq7AgBw\nNKlOqrP/T9PvWVxcXFhYWFpaKoQoLCwsLCw0GAx37tyJjY1dsGDBkCFDCv+tqqpKCLFz5871\n69dbLz958qTJZOrXr5+1RalUzps3LykpKSsrq6CgID4+3sfHJyoqqsXfSrNmCE+ePGmprz6z\n2fztt99alrgBAAAdkdEolpz2qa0TQgiFQkhCLHnE2FNtdnZdAGB/zVkUtAV9Nn3C8OHDr1y5\nYvndso7o+++/P2TIkEuXLun1+vp7+G3cuDEuLi4rK6ukpGT+/PmWxqKiIoVCERQUVL/PxYsX\nV1RUvPzyyzdv3hw+fPihQ4fqb2B4vxRNb86oUCiavv7RRx89efJki9++YykvLzeZTM6uoh1R\nq9U+Pj6VlZXV1dXOrgV2Zhlcg8FgNDJr4GoYXBfWgsGde0QjxA//rFeIzSMM9i8OrcP/c10Y\ng2uTSqXy9fW1b5+TVv/uWqmd//vmo1Zlr5pi3z7b2D1mCE+fPp2bmzt//vzo6Gh/f//6LykU\nih49esyaNcuR5QEAAEf5c4FomAaFEJJYdUr9q0f5iykAl9OMfeRb0mcHd49AOHTo0KFDhx48\nePDdd9+tf+sqAADo6L4q9rbZ/r2xhYuQA0B7Jkmi6bsjW9ZnR9es/+JnZmY6ug4AANDGNKra\nMpON5eUUCjs/YwMA7QIzhLY0FQjDwsJiY2OXLVtm2fSwMXl5efauCgAAONycsOpfnlLd1SxF\nht5xQjUA4GBOWVSm/WsqEHbp0sXLy0vU2/cQAAC4jC5q0c2z9ka1sn6jWin9NLjD/4M3ANjg\niBlC0eH/g9lUIDx27FiDXwAAgCt567E7X15TfFHgLdUJhRAPd62bOYDpQQCuSZIkBzxD2OED\nIRvTAwAgaz8NljYNr9r8o6pNP6oiDQJwadK/Jgnt+9PBsTE9AAAAANcnSRLPEN7tHoHQujH9\nsGHDbJ7w6KOP2rkiAAAAALA7Vhm1hY3pAQAAALg+Vhm1qVUb0xsMhqKiIofVBgAAAAB2wgyh\nLa3amP5vf/vbf//3f5eWltq1JAAAAACwM1YZtalZgVAIceDAgY8//vjq1at1df+aFa2trf32\n2289PT0dVhsAAAAA2AkzhLY0KxCmp6f/7Gc/c3d3DwwMLCws7NGjR1lZmdFofPrpp998801H\nlwgAAAAArccM4d2atQ/hb37zm2eeeaasrKygoECpVP7pT3+qrKzcsGGDJEk//vGPHV0iAAAA\nALRanQN+5BEIv/vuu7i4uE6dOlkOJUlyd3ePj49/+OGHly1b5sjyAAAAAMAeHLErvUxmCE0m\nk1KptPyu0Whu3bpl+X3SpEn79u1zVGkAAAAAYCeSYzj7Y7VWswLhwIEDd+7cWVNTI4TQ6XR/\n+tOfLO1lZWXl5eUOrA4AAAAA7IIZQluaFQgXLlz42WefjR8/XggxceLENWvWzJ49Oykp6fXX\nXx86dKiDKwQAAACA1nLWDOH58+eHDx/u7v6D5TyvX78+bdq0gICAzp07jxw58vjx4w2uysnJ\nUdxl06ZNQoihQ4fWb/Tx8WnN19KsVUZfeukld3f3/Px8IcTSpUuPHTu2fft2IYROp1u/fn1r\n3h4AAAAA2oIkCanO3n3eo8OMjIw33nhj7NixJ06cqN8eHR2tVqszMzM7deqk1+sjIyPz8/M1\nGo31hBEjRhQUFFgP8/Pzx40bN3r0aCFEWVnZhg0bJkyYYHnJza1Zk3yNae4+hFOnTrX84u3t\n/eWXX168eNFkMvXt21elUrXm7QEAAACgDTjikb97dlhdXX3s2LFTp0599NFH1saysrLQ0NBV\nq1YNHDhQCLF69eo9e/acPXt22LBh1nM8PT1DQkKsh9OnT1+0aNGgQYMsl/fp06f+q63RwjTZ\nt2/fgQMHkgYBAAAAdAxSnUN+mhQTExMaGtqg0c/Pb+/evZY0KIS4du2aUqnU6XSNdZKenn7x\n4sXly5cLIaqrq2/fvv35558/+uijPXv2nDRp0nfffdeab6WpGcKwsLDmdJGXl9eaCgAAAADA\n0XoGajt5e1oPzbV13xV8f7+d9An291T9J0N5ejb3jsvGlJWVzZgxY9GiRYGBgTZPqK2t1ev1\nCQkJHh4eQoiKioqAgICampqtW7dKkpSYmPiTn/wkLy+vS5cuLSugqQ/g7+/fsk4BAAAAoF3p\nGegX6NfZenin2vTd1X/ebyd9enTtrPGyHiqVitaUlJeXFxUVNXbs2DVr1jR2zqefflpVVRUT\nE2M57NatW3FxsfXVjIyMoKCgvXv3zpgxo2U1NBUIDx8+3LJOAQAAAKBdyTn1XeH3t1rZyZ/+\ndq7+oY+3569in2lZV4cOHZoyZYper4+Pj2/itN27d0+aNKnBIqVWnTp1Cg0Nrb/8zP1q1Yo0\nAAAAANAxSI7ZirBFDh8+PHny5N27dzedBm/dupWVlRUVFWVt+b//+79Zs2ZZtogXQhgMhqtX\nr/bp06dlZYjmrzIKAAAAAB2XI1YZvWcgLC4uNpvNpaWlQojCwkIhRJcuXZRKZWxs7IIFC4YM\nGWJpFEJotVqNRrNz506DwTB//nxL48mTJ00mU79+/awdBgUF7du3r6amRq/Xm0ym5cuX+/n5\nTZo0qcWfgBlCAAAAAHJQ55ifpgwfPlyn082cObO2tlan0+l0uh07dhw5cuTSpUt6vV5XT2pq\nqhAiKytr//791suLiooUCkVQUJC1pWvXrtnZ2deuXXv00Ud//OMfm83m3Nxcb2/vFn8pzBAC\nAAAAkIFW3OHZeJ/3eD0/P7+RWmxfmZ6eXv/wpZdeeumllxqc8/DDD2dnZzezwHsiEAIAAACQ\nAWdsTN/+EQgBAAAAyIEjZggJhAAAAADQ7jliURlmCAEAAACgI5AcMKFHIAQAAACA9k+S6iTp\nHouCtqBP+3bY9giEAAAAAGTAIauMMkMIAAAAAO0eM4Q2EQgBAAAAyAAzhLYQCAEAAAC4Pkmy\n/6KgrDIKAAAAAB2CJITd8xuBEAAAe6uuE2kX1Neq3DTu4qe66ocbdqetAAAgAElEQVS1tc6u\nCADgArhl1AYCIQCgfblapXzn716SJAmhKBFi+zmvAb618x684+y6AAAdGxvT20QgBAC0L+99\n4yVJQgiFteV8ufLY96rh3U3OKwq4h+SzXmcrlF1U5rgBRn8fZ1cDwCZHLCrDLaMAANhRoUFp\nsrWC9x8LCYRopy7dcX/va7Xlr4Q3qt31f/fpoqp9exhz2kC745htJwiEAAC0jiSJv36vOl3q\nfsPodqNaobB1jsFksxlwvve+9qyrP6MtxC2TMjNfPPOAsyoC0Ai2nbCFQAgAcKaaOsWms+p/\nVCgth43FPrWyzSoC7sO5Mg9J2PhXjP1FPs88YHBCQQAaJwkh2fsOT7t32PYIhAAAp7licPvg\nnLehGbeCjggwO74c4L59Vepms50ZbaA9qqsTdXa+ZdT+HbY5AiEAwDlq6hQ7z6ubkwaFEI93\n4wFCtEdPd6379oaN9g4/ZQC4JlYZtcH2P2sBAOBo35UrS6ub9ceQv2ddd3WH/ydYuKSBfjU2\n099/9eB+UaAdkhzzcw/nz58fPny4u/sPpuKuX78+bdq0gICAzp07jxw58vjx4w2uysnJUdxl\n06ZN9c/57W9/q1Aofve737XmSyEQAgCco5lzg0qFNGdQtYNrAVpu0SPVP7hDVCH8VLXP9XRa\nPQAaZVlUxu4/TcrIyHj66acHDBjQoD06OrqgoCAzM/PUqVMhISGRkZFVVVX1TxgxYkRBPX/5\ny198fHxGjx5tPeGf//zn0qVLvby8WvmtEAgBAM4R4HXvf1XVaeqe7mFKzvOcf1Sz8mvv/1es\nquvw9+bA1fT2Mm8eYXhQW6tyE36q2mVDDSvZcwJolyTHaPpNq6urjx07NmHChPqNZWVloaGh\nycnJjzzySN++fVevXl1SUnL27Nn653h6eobUk5iYuGjRokGDBllPmDt37osvvti5c+dWfi08\nQwgAcI4HOtUO6mI+e+sHfxK5KySlQgRrpLHBNQ92Mf/+qmf2NQ/LS8V3FBmXPMtrFFGhNc6o\nF2jK6wMJgUC7Jzlgl4h7dRgTEyOEOHXqVP1GPz+/vXv3Wg+vXbumVCp1Ol1jnaSnp1+8ePHA\ngQPWls8///zUqVNpaWkfffRRCyv/NwIhAMA5FEK80r/608vixA13SQilQjwVaHq+Z42H27/+\ncC0xuh26rmpw1ZfXPJ4KMGk9mSgEANyfqeOfqr8GsOH2nc17Dt5vJ69Oiuim9a3X0No/j8rK\nymbMmLFo0aLAwECbJ9TW1ur1+oSEBA+Pf/0L6c2bN+Pi4tLS0jQaTSvfXRAIAQBOpHGXXuln\n/FlvRVm1wl8tqdx+8Mdq4W0bzzXUSeLabaXWk10oAAD359DRM9+XllsPzbW1LZgw/GPuCVW9\n5WG81B4/nzquxSXl5eVFRUWNHTt2zZo1jZ3z6aefVlVVWWYaLRYuXPhf//VfY8eObfH71kcg\nBAA4madSCvK28UeyRyPPuTfIjQAANMc/S29eLbK1Ucz9uH6jtP5hJ03L13Q5dOjQlClT9Hp9\nfHx8E6ft3r170qRJ1kVKs7KyMjMzv/322xa/bwMsKgMAaKd6d6rVuDfMfp1VUq9ObEEBALh/\nzlhltDGHDx+ePHny7t27m06Dt27dysrKioqKsrakpKTcunWrf//+/v7+/v7+33//fUxMzKRJ\nk1pWhmCGEADQbqmVUky/6h3nPU11/3rkw8NNiuln9GCGEADQAq3Ib0312aTi4mKz2VxaWiqE\nKCwsFEJ06dJFqVTGxsYuWLBgyJAhlkYhhFar1Wg0O3fuNBgM8+fPtzSePHnSZDL169fP2uHm\nzZvfffdd6+Gjjz66evXq6OjoFn8CAiEAoP0arDUnPFJ39J/uZdVuXdV1P+rOcjIAgBaSxL13\nibjvPu/V4fDhw69cuWL53bKO6Pvvvz9kyJBLly7p9Xq9Xm89c+PGjXFxcVlZWSUlJdZAWFRU\npFAogoKCrKf5+fn5+flZD93c3Lp27erv79/ij0AgBAC0a109655lnwkAQOtJkpDs/dDBvQJh\nfn5+I9fZvjA9Pb3+4UsvvfTSSy810X9xcXHTBdwTgRAAAACA62vOPvIt6NO+HbY9AiEAAAAA\nGXDGM4Ttn+sHQkmSdu3a9de//jUkJCQuLq7+HbcAAAAAZMIhM4St3pje6ZwTCAsKCtLS0s6d\nOydJUq9evV5++eWwsDDLSwcPHty3b19paWlwcHBMTMywYcMaXFtYWJiamnr+/Hmz2Wy5dtCg\nQU30mZube+3atS1btuTk5Gzbtm3ZsmVCiIULF169enXDhg09evSw9hwXFxcZGTluXMt3lgQA\nAADQXknC/vmtwwdCJ+xDaDabExISNBrNO++8s3bt2m7duiUmJt65c0cIcejQoYyMjNdee23r\n1q0RERHbt2+/fft2/WslSUpKStJqtcnJyWlpaYMHD37rrbcqKyub6PPChQtPPvmkUqkcNWrU\nhQsXrF15enpu3ry5jT87AAAAAOeQHMLZn6q1nBAIq6qqoqOj58yZExwcHBQUNHny5KqqqqKi\nIiFERkZGbGxseHh49+7do6Ojk5OTvb29619bUVFRXFwcERHh7e3t6ek5fvx4o9FYVFTURJ+9\nevU6cuSI2WzOzc3t06ePtavo6Oj8/Pzs7Ow2/vgAAAAAnEByzN70HZwTbhn19fWdMGGC5ffK\nysovvvgiJCREp9OVlpZaVk2dN29eUVFRz549Z86cab2V1HptWFhYZmZmcHCwSqXKzMwMCAh4\n4IEHPDw8bPYphBg9evTly5fnzJmj0+ni4uKsXWk0munTp6ekpAwbNszX19dmqZcuXTpw4ID1\ncPz48fVvMYVSqRRCeHp6uru7/sOocmMZXA8PD8svcCUMrgtjcF0Yg+vCGNw2xKIyNjjt7/F1\ndXWTJ082mUyDBw9euXKlSqUqLS0VQmRnZy9ZssTX1zc9PT0xMXHr1q0N0trSpUtXrFjx4osv\nCiG0Wu2KFSs8PDwa61MI4ebmNmvWrFmzZt1dQ0RERE5OTnJy8uLFi20WeeXKlbS0NOvhE088\nUX+OERbW7x+uh8F1YQyuC2NwXRiD68IY3AbMZrPd+3TIHZ4EwhZzc3Nbv379zZs3Dxw4sHz5\n8rVr11rap0yZEhISIoR49dVXv/rqqxMnTowZM8Z6ldlsTkpKCgsLe/vtt1Uq1cGDB/V6/caN\nG7Varc0+fXx8mi5j7ty58fHxJ06cCA8Pv/vVoUOHfvDBB9ZDnU5XXl7e+s/uMjw8PLy8vG7f\nvm0ymZxdC+zMMrh37typqWFDcFfD4LowBteFMbgujMG1SalU3vNv8vdNqnPAxvT27rDNOfNO\nv5CQkJCQkAcffHDatGm5ublPPPGEEEKj0VheVSqVfn5+N2/erH/JN998c/ny5TVr1qjVaiHE\nCy+88Mc//vHw4cNRUVE2+4yMjGy6hqCgoKlTp27ZsmXz5s0KhaLBq35+fo8//rj1sLy8nORT\nn+XGhtraWr4W18PgujAG14UxuC6MwXVhDG6bsTzzZ/c+7dth23PCojJff/317Nmzq6urLYcK\nhcLyBJqfn59Wq83Ly7O019TU3LhxIyAgoP61lnneurr/BHHLbHJjfTbHhAkTvL29d+/ezX3b\nAAAAgMtyxIoyHf+WUScEwn79+hmNxnXr1hUUFBQXF+/YscNoND722GNubm5RUVHp6emnT58u\nKSnZtm2bWq227EOYlZW1f/9+IURYWJhWq01JSTEYDDU1NZ9//nlVVVV4eHhjfTanHqVSGR8f\nf/DgQctDjAAAAABcj2N2nejwgdAJt4z6+PisXLkyNTV10aJFCoUiNDQ0ISEhMDBQCDFx4sTb\nt2+/9957BoNhwIABq1atstwaevr06YqKiqioKG9v76SkpLS0tDlz5tTW1oaGhur1+qCgICFE\nY302R//+/cePH2/JnAAAAABckOSANWA6fiBUuECobTM8Q9iAWq328fGprKy03qwLl2EZXIPB\nYDQanV0L7IzBdWEMrgtjcF0Yg2uTSqVqbGe4FnviudlXrhXbt8/OPt7f5abbt882xvZxAAAA\nAFyfJNVJ9l4U1AVm1wiEAAAAAOSAjeltIBACAAAAcH1sO2ETgRAAAACADDhil4iOP0PohG0n\nAAAAAKCtSY7ZivBezp8/P3z48AbbpF+/fn3atGkBAQGdO3ceOXLk8ePHG1yVk5OjuMumTZuE\nEOfOnXvuuee6du3q5+c3evToo0ePtuZbIRACAAAAcH1O2YcwIyPj6aefHjBgQIP26OjogoKC\nzMzMU6dOhYSEREZGVlVV1T9hxIgRBfX85S9/8fHxGT16dE1NTURERJcuXY4cOXL8+HGdTjd+\n/PjKysoWfy0EQgAAAAByIDnmpynV1dXHjh2bMGFC/caysrLQ0NDk5ORHHnmkb9++q1evLikp\nOXv2bP1zPD09Q+pJTExctGjRoEGDysvL33jjjc2bNw8YMKBv377Lly+/devWP/7xjxZ/KTxD\nCAAAAMD1NWdCrwV9Nn1CTEyMEOLUqVP1G/38/Pbu3Ws9vHbtmlKp1Ol0jXWSnp5+8eLFAwcO\nCCG6dev25ptvWtrLysrWr18fFhY2cODAFn8EAiEAAAAA17dq8Ryl8j83SFbdNu767MD9dvLf\nz43t2qWz9bCurrUJs6ysbMaMGYsWLQoMDLR5Qm1trV6vT0hI8PDwqN+o0Wiqq6tHjhyZnZ3t\n6enZ4gIIhAAAAABc37jRP2rQMmHcKGcU8h95eXlRUVFjx45ds2ZNY+d8+umnVVVVlplGK6VS\nefr06eLi4k2bNo0aNer48eNarbZlNfAMIQAAAAC0tUOHDj311FPz5s374IMPFApFY6ft3r17\n0qRJDRYpFUKEhYWNGjUqIyPjxo0bH330UYvLIBACAAAAQJs6fPjw5MmTd+/eHR8f38Rpt27d\nysrKioqKsrZ8+eWXffv2vX37tuXQzc1NpVI1kSfviVtGAQAAAMAhiouLzWZzaWmpEKKwsFAI\n0aVLF6VSGRsbu2DBgiFDhlgahRBarVaj0ezcudNgMMyfP9/SePLkSZPJ1K9fP2uHw4YNMxgM\nr7zyyltvvaVWqzds2FBVVfXMM8+0uEJmCAEAAADAIYYPH67T6WbOnFlbW6vT6XQ63Y4dO44c\nOXLp0iW9Xq+rJzU1VQiRlZW1f/9+6+VFRUUKhSIoKMjaotVqs7OzKysrH3/88Ycffvhvf/vb\n/v37+/Tp0+IKFXZfetWFlZeXm0wmZ1fRjqjVah8fn8rKyurqamfXAjuzDK7BYDAajc6uBXbG\n4LowBteFMbgujMG1SaVS+fr6OrsKWWCGEAAAAABkikAIAAAAADJFIAQAAAAAmSIQAgAAAIBM\nEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJFIAQAAAAA\nmSIQAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAA\nADJFIAQAAAAAmSIQAgAAAIBMuTu7AAAAIEe/Oulzq0ZIklAopDGB5gm9qp1dEQDIETOEAACg\nrcUf8blZLSRJCCEkSZFdpHrv72pnFwUAckQgBAAAbeqDs551dzVeMnDXEgA4AYEQAAC0qbxy\nG9lPEuJ3Vz3avhgAkDkCIQAAaBeMdQpnlwAAskMgBAAAbSrIW7q7UaEQE3uwrgwAtDUCIQAA\naFPLht6+u9FfLTy4YxQA2hyBEAAAtLVV4QaPf98fqhBS7061bz1icGpFACBTrOgFAADamtZD\nvD+CBAgAzscMIQAAAADIFIEQAAAAAGSKQAgAAAAAMkUgBAAAAACZIhACAAAAgEwRCAEAAABA\npgiEAAAAACBTBEIAAAAAkCkCIQAAAADIFIEQAAAAAGSKQAgAAAAAMkUgBAAAAACZIhACAAAA\ngEwRCAEAAABApgiEAAAAACBTBEIAAAAAkCkCIQAAAADIFIEQAAAAAGSKQAgAAAAAMkUgBAAA\nAACZIhACAAAAgEwRCAEAAABApgiEAAAAACBTBEIAAAAAkCkCIQAAAADIFIEQAAAAAGSKQAgA\nAAAAMkUgBAAAAACZIhACAAAAgEwRCAEAAABApgiEAAAAACBTBEIAAAAAkCkCIQAAAADIFIEQ\nAAAAAGSKQAgAAAAAMkUgBAAAAACZIhACAAAAgEwRCAEAAABApgiEAAAAACBTBEIAAAAAkCkC\nIQAAAADIFIEQAAAAAGSKQAgAAAAAMuXu7AIcTpKkXbt2/fWvfw0JCYmLi/Pz83N2RQAAAADQ\nLjg2EBYUFKSlpZ07d06SpF69er388sthYWGWlw4ePLhv377S0tLg4OCYmJhhw4Y1uLawsDA1\nNfX8+fNms9ly7aBBg4QQZWVlKSkpZ86cqamp6d279/Tp0/v37y+EmDdvXn5+vvVytVr9ySef\nCCFyc3OvXbu2ZcuWnJycbdu2LVu2TAixcOHCq1evbtiwoUePHtZL4uLiIiMjx40b59DvBAAA\nAADaCQfeMmo2mxMSEjQazTvvvLN27dpu3bolJibeuXNHCHHo0KGMjIzXXntt69atERER27dv\nv337dv1rJUlKSkrSarXJyclpaWmDBw9+6623KisrhRCrVq0qKSlJTExct26dv79/UlKS0WgU\nQhgMhtmzZ6f829atWy1dXbhw4cknn1QqlaNGjbpw4YL1LTw9PTdv3uy4jw8AAAAA7ZwDA2FV\nVVV0dPScOXOCg4ODgoImT55cVVVVVFQkhMjIyIiNjQ0PD+/evXt0dHRycrK3t3f9aysqKoqL\niyMiIry9vT09PcePH280GouKiiorK7t16zZ37tzevXsHBQXFxMRUVFQUFBQIISorKwMDA/3/\nzXpraK9evY4cOWI2m3Nzc/v06WN9i+jo6Pz8/OzsbMd9AwAAAADQnjnwllFfX98JEyZYfq+s\nrPziiy9CQkJ0Ol1paWlxcbEQYt68eUVFRT179pw5c6b1VlLrtWFhYZmZmcHBwSqVKjMzMyAg\n4IEHHvDw8LDc82lRWlrq5ubm7+9vMpmqq6uPHj364YcfVlZW9u3bNyYmJjg4WAgxevToy5cv\nz5kzR6fTxcXFWa/VaDTTp09PSUkZNmyYr6+vzY9QWVlZWFhoPezatauHh4f9vqEOz83NTQih\nVCrd3V3/YVS5sQyum5sbg+t6GFwXxuC6MAbXhTG4Nlm+FrQBh//Prq6ubvLkySaTafDgwStX\nrlSpVKWlpUKI7OzsJUuW+Pr6pqenJyYmbt26tUEqW7p06YoVK1588UUhhFarXbFiRYMwVllZ\nuXHjxueff16r1ZaXl3fp0sVsNr/++utCiI8//njZsmVbtmzRaDRubm6zZs2aNWvW3bVFRETk\n5OQkJycvXrzYZvEnTpyo/9IHH3zw+OOPt/orcTXe3t4NJnjhMhhcF8bgujAG14UxuC6MwW3A\nbDY7uwS5cHggdHNzW79+/c2bNw8cOLB8+fK1a9da2qdMmRISEiKEePXVV7/66qsTJ06MGTPG\nepXZbE5KSgoLC3v77bdVKtXBgwf1ev3GjRu1Wq3lhMLCwpUrVz788MOxsbFCCF9f3127dlkv\nX7JkSWxs7JEjR8aOHdt0eXPnzo2Pjz9x4kR4ePjdrwYHB0+cONF6qNVqLc8rwkKpVKpUKpPJ\nVFtb6+xaYGcMrgtjcF0Yg+vCGFwXxuA2hinTttEW33JISEhISMiDDz44bdq03NzcJ554Qgih\n0WgsryqVSj8/v5s3b9a/5Jtvvrl8+fKaNWvUarUQ4oUXXvjjH/94+PDhqKgoIcSZM2feeeed\nn/3sZ88++6zNd/Ty8urWrVtJSck9awsKCpo6deqWLVs2b96sUCgavNq/f//ly5dbD8vLyw0G\nw318clenVqtVKpXRaKyurnZ2LbAzy+BWV1fzjyCuh8F1YQyuC2NwXRiDa5NKpbIEATiaA+/N\n/frrr2fPnm2NCgqFwpLy/fz8tFptXl6epb2mpubGjRsBAQH1r5UkSZKkuro6a4t11vjs2bO/\n/vWvFy5cWD8NXrlyZdOmTdZzjEbjjRs3AgMDm1PnhAkTvL29d+/erVQqW/hRAQAAAKADcuAM\nYb9+/YxG47p166ZNm6ZSqfbv3280Gh977DE3N7eoqKj09HTLzOHHH3+sVqst+xBmZWUZjcao\nqKiwsDCtVpuSkvLKK694eHj84Q9/qKqqCg8Pr6mpWbdu3XPPPdezZ0/rBKCPj4+fn9/Ro0fN\nZvPUqVNra2t37drl4+Pzox/9qDl1KpXK+Pj4X/ziFz4+Po77NgAAAACgvXFgIPTx8Vm5cmVq\nauqiRYsUCkVoaGhCQoJl1m7ixIm3b99+7733DAbDgAEDVq1aZZkRPn36dEVFRVRUlLe3d1JS\nUlpa2pw5c2pra0NDQ/V6fVBQ0JkzZ4qLi/fs2bNnzx7rG7322muRkZGW91qwYIFKpRowYMDq\n1as9PT2bWWr//v3Hjx+/f/9+R3wPAAAAANA+KSRJcnYNHUZ5ebnJZHJ2Fe2IWq328fGprKzk\nGULXYxlcg8HA8wyuh8F1YQyuC2NwXRiDa5NKpWpsZzjYF/t7AAAAAIBMEQgBAAAAQKYIhAAA\nAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJFIAQAAAAAmSIQAgAAAIBMEQgB\nAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJFIAQAAAAAmSIQ\nAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJF\nIAQAAAAAmSIQAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABk\nikAIAAAAADJFIAQAAAAAmSIQAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAA\nyBSBEAAAAABkikAIAAAAADJFIAQAAAAAmSIQAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAA\nAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJFIAQAAAAAmSIQAgAAAIBMEQgBAAAAQKYIhAAA\nAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJFIAQAAAAAmSIQAgAAAIBMEQgB\nAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJFIAQAAAAAmSIQ\nAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABkikAIAAAAADJF\nIAQAAAAAmSIQAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAAyBSBEAAAAABk\nikAIAAAAADJFIAQAAAAAmSIQAgAAAIBMEQgBAAAAQKYIhAAAAAAgUwRCAAAAAJApAiEAAAAA\nyBSBEAAAAABkikAIAAAAADJFIAQAAAAAmXJ3dgEOJ0nSrl27/vrXv4aEhMTFxfn5+Tm7IgAA\nAABoF5wTCAsKCtLS0s6dOydJUq9evV5++eWwsDDLSwcPHty3b19paWlwcHBMTMywYcMaXFtY\nWJiamnr+/Hmz2Wy5dtCgQUKIsrKylJSUM2fO1NTU9O7de/r06f379xdC5ObmXrt2bcuWLTk5\nOdu2bVu2bJkQYuHChVevXt2wYUOPHj2sPcfFxUVGRo4bN66NvgUAAAAAcCon3DJqNpsTEhI0\nGs0777yzdu3abt26JSYm3rlzRwhx6NChjIyM1157bevWrREREdu3b799+3b9ayVJSkpK0mq1\nycnJaWlpgwcPfuuttyorK4UQq1atKikpSUxMXLdunb+/f1JSktFoFEJcuHDhySefVCqVo0aN\nunDhgrUrT0/PzZs3t+1HBwAAAIB2xAmBsKqqKjo6es6cOcHBwUFBQZMnT66qqioqKhJCZGRk\nxMbGhoeHd+/ePTo6Ojk52dvbu/61FRUVxcXFERER3t7enp6e48ePNxqNRUVFlZWV3bp1mzt3\nbu/evYOCgmJiYioqKgoKCoQQvXr1OnLkiNlszs3N7dOnj7Wr6Ojo/Pz87OzsNv74AAAAANBO\nOOGWUV9f3wkTJlh+r6ys/OKLL0JCQnQ6XWlpaXFxsRBi3rx5RUVFPXv2nDlzpvVWUuu1YWFh\nmZmZwcHBKpUqMzMzICDggQce8PDwsNwLalFaWurm5ubv7y+EGD169OXLl+fMmaPT6eLi4qzn\naDSa6dOnp6SkDBs2zNfX12apx44dW716tfUwMTFxyJAh9vsmOjyFQiGE0Gg0DXI7XIBlcL29\nvb28vJxdC+yMwXVhDK4LY3BdGINrU21trbNLkAunLSpTV1c3efJkk8k0ePDglStXqlSq0tJS\nIUR2dvaSJUt8fX3T09MTExO3bt3aIK0tXbp0xYoVL774ohBCq9WuWLHCw8Oj/gmVlZUbN258\n/vnntVqtEMLNzW3WrFmzZs26u4aIiIicnJzk5OTFixfbLNJsNlvuR7Wora11c2Nd1oYUCoXl\nP2RwPQyuC2NwXRiD68IYXBfG4DZQV1fn7BLkwmmB0M3Nbf369Tdv3jxw4MDy5cvXrl1raZ8y\nZUpISIgQ4tVXX/3qq69OnDgxZswY61VmszkpKSksLOztt99WqVQHDx7U6/UbN260ZD8hRGFh\n4cqVKx9++OHY2NjmlDF37tz4+PgTJ06Eh4ff/epTTz315z//2XpYXl5uSa2wUKvVPj4+BoOh\nurra2bXAziyDW1VVZXkWF66EwXVhDK4LY3BdGINrk0qlauwmPtiXM+e7QkJChgwZsmTJkvLy\n8tzcXMuGEBqNxvKqUqn08/O7efNm/Uu++eaby5cvz5w509fX19vb+4UXXvD09Dx8+LDl1TNn\nzvziF7+Iior6+c9/3sx/YgkKCpo6deqWLVuMRiP/KgMAAIgDTa4AABfySURBVABAVpwQCL/+\n+uvZs2db55QUCoW7u7sQws/PT6vV5uXlWdprampu3LgREBBQ/1pJkiRJqj+DbDabLb+cPXv2\n17/+9cKFC5999tn7qmfChAne3t67d+9WKpUt/lAAAAAA0OE4IRD269fPaDSuW7euoKCguLh4\nx44dRqPxsccec3Nzi4qKSk9PP336dElJybZt29RqtWUfwqysrP379wshwsLCtFptSkqKwWCo\nqan5/PPPq6qqwsPDa2pq1q1b99xzz/Xs2bPk35o57a5UKuPj4w8ePMjtoAAAAABkxQnPEPr4\n+KxcuTI1NXXRokUKhSI0NDQhISEwMFAIMXHixNu3b7/33nsGg2HAgAGrVq1Sq9VCiNOnT1dU\nVERFRXl7eyclJaWlpc2ZM6e2tjY0NFSv1wcFBZ05c6a4uHjPnj179uyxvtFrr70WGRnZnJL6\n9+8/fvx4S+YEAAAAAJlQSJLk7Bo6jPLycpPJ5Owq2hHLM9CVlZUsKuN6rCsG8YC762FwXRiD\n68IYXBfG4NrEojJthk0UAAAAAECmCIQAAAAAIFMEQgAAAACQKQIhAAAAAMgUgRAAAAAAZIpA\nCAAAAAAyRSAEAAAAAJkiEAIAAACATBEIAQAAAECmCIQAAAAAIFMEQgAAAACQKQIhAAAAAMgU\ngRAAAAAAZIpACAAAAAAyRSAEAAAAAJkiEAIAAACATBEIAQAAAECmCIQAAAAAIFMEQgAAAACQ\nKQIhAAAAAMgUgRAAAAAAZIpACAAAAAAy5e7sAgAAAAC0R9/fcTt2Q1VWrfD3rHsywKT1lJxd\nEeyPQAgAAACgoa9L3X/7nadZUlgOD11X/Xygsb9vrXOrgt1xyygAAACAHzCYxEf/+E8aFELU\n1CnSLqjrt8A1EAgBAAAA/MDFCuUdc8Psd6tGcdVAfHA1jCgAAACAHzDV2Z4JrOGOUZdDIAQA\nAADwAzpNnRCiwRoySjcR4lPnlHrgOARCAAAAAD8Q6F03MrCmwSxhZEi1jzsLjboaVhkFAAAA\n0NCkB2q6eUmHi1Wl1Ypu6rrRPczDu5mcXRTsj0AIAAAAoCGlm3g6yPR0ECHQxXHLKAAAAADI\nFIEQAAAAAGSKQAgAAAAAMkUgBAAAAACZIhACAAAAgEwRCAEAAABApgiEAAAAACBTBEIAAAAA\nkCkCIQAAAADIFIEQAAAAAGSKQAgAAAAAMkUgBAAAAACZIhACAAAAgEwRCAEAAABApgiEAAAA\nACBTBEIAAAAAkCkCIQAAAADIFIEQAAAAAGSKQAgAAAAAMkUgBAAAAACZIhACAAAAgEwRCAEA\nAABApgiEAAAAACBTBEIAAAAAkCkCIQAAAADIFIEQAAAAAGSKQAgAAAAAMkUgBPD/27vzqKjK\nx4/jD+uMrDGCQI6hWII2FqAouVZidTRcKhW11ELcKistl/p6XOqUSyfpVByzwiKPZie3suSY\nHu1gmrkgYiaLoiE4LAaJ0CAzzPeP++v+OH4d0pgled6vv2aee3nmgc8Z4DP33hkAAABIikII\nAAAAAJKiEAIAAACApCiEAAAAACApCiEAAAAASMrT1Qu4lWg0Gk9PfmL/78yZM0eOHElISIiI\niHD1WmBnZ8+ePXz4cN++fTt37uzqtcDOlHD79OnTpUsXV68FdlZcXPzzzz8TbpukhBsfHx8Z\nGenqtcDOzp07d+jQIcK9hoeHh6uXIAvqzU3QarWuXsK/S35+fnp6emhoaI8ePVy9FthZQUFB\nenp6SEjI3Xff7eq1wM4KCwvT09ODg4MNBoOr1wI7KyoqSk9P1+l0hNv2nD17Nj09/dVXX+3Z\ns6er1wI7U8JduHAh4cIlOGUUAAAAACRFIQQAAAAASVEIAQAAAEBSblar1dVrwK2qrq6upqYm\nKCjIx8fH1WuBnRFuG0a4bVh9fX11dTXhtklKuLfddpuvr6+r1wI7I1y4FoUQAAAAACTFKaMA\nAAAAICkKIQAAAABIis8hhH1YrdbMzMwff/xRr9c/99xzOp3O1SuCkxB9G0a4bRjhtmGE24YR\nLhyBawghSktLV69eXVRUtG3bNnXw999/z8jIyM3NvXr1amRk5NNPP92tW7cWxvft23fgwIH5\n8+fv27fv559/Xrhw4Zw5cyIiIl544QV1zmnTpoWGhr7++uvqyPz58zt06DB37lwnfrttzU3F\nZ2v/5oxG47p1606dOtXQ0NCrV68ZM2YEBga2ME70juO0cC9cuLBu3br8/Hyz2dylS5ennnqq\nR48egnAdyTnh5uXlvfbaa9fsOX369OHDhxOu4zjtmSuE+O6777Zu3Xrp0qWOHTtOmjQpPj5e\n8Mx1JKeFO3v27HPnzqm7abXaL7/8UhAuHINTRmWXnZ396quv6vX6a8bfeOONqqqqpUuXpqWl\nBQcHL1u2zGQytTBeWFjYv39/Dw+P+++/v7CwUAgRHx9/7Ngx9RUHo9FYU1Pz66+/NjQ0KCP1\n9fUFBQXKXy/8Mzcbn639VY2NjUuWLLFYLMuXL1+1atWVK1dWrFjRwrggeodxWrhWq3XZsmVB\nQUFr16797LPPDAbDkiVLamtrBeE6jNPCjY6Ozmhm+fLlWq32nnvuEYTrME4LVwixZ8+eTZs2\nTZ8+fc2aNYmJiR999FF9fb0gXIdxZrhXrlyZNm2a+uRds2aNMk64cAQKoewaGxvffvvthISE\n5oO1tbUhISHPPvtsZGRkeHj4pEmTLl++XFJSYmtcCNGlS5cDBw6YzeYffviha9euQoj4+Pjq\n6uri4mJlzqNHj/bo0SMkJCQvL08ZycnJsVqtcXFxzv2O25Sbis/W/s0VFxeXlZXNnDmzY8eO\nyiuOJ0+ePH/+vK1xQfQO47RwL1++bDQaExMTfXx8NBrNsGHDTCbTxYsXBeE6jNPC9fLyCm5m\n48aNo0eP7tSpkyBch3FauEKITZs2TZ48uXfv3h06dBg5cuTatWuVjxshXAdxZri1tbVhYWHq\nk1c9NZRw4QgUQtk9+OCDISEh1wz6+/svXLhQ+adBCHHp0iV3d/fg4GBb48o8wcHBM2bMyM7O\nnjFjhhCia9euOp3u6NGjys5Hjx41GAwGg0EdycnJ6d69u5+fnxO+zbbqpuKztX9zjY2NQghv\nb2/lblBQkIeHR1FRka1xQfQO47RwAwMDo6Ojs7KyamtrTSZTVlZWaGho586dBeE6jNPCbb5P\ndnb2xYsXx4wZo66BcB3BaeFeunTJaDQKIWbPnj1mzJiXX3759OnT6hoI1xGc+Te3oaHh4MGD\nL774YkpKyltvvVVaWqqugXBhdxRC/I3a2tr33ntv1KhRQUFBLYy7u7unpqZ+/PHHixcvbt++\nvRDCzc2td+/eyu8js9l88uTJmJiYmJiYY8eOKTMcO3aMExgczVZ8tkRGRgYEBGzYsMFsNpvN\n5k2bNimT2BoXRO869gpXCLFgwYKioqKJEyeOHTs2KytrwYIFyj8ohOsqdgxX0dTUtGHDhuTk\nZE/P/3szOcJ1FXuFe+nSJSHE7t27582bl5GRERUVtXTp0j/++EMQruvYK9z6+vrbbrvNbDbP\nmjVr/vz5V69eXbhwYV1dnSBcOAaFEC25cOHCyy+/bDAYJk+efCPj14iPjz99+nRdXd0vv/yi\n0Wi6du167733lpeXG43G3377raqqit9QDnWDMTXXrl27BQsWHDt2bMyYMU8++aQQokOHDh4e\nHrbGbc1D9I5mx3DNZvOyZcuio6M///zzL774IikpafHixdXV1bbmIVxHs2O46g4//vijyWR6\n4IEHWp6HcB3N7uGOGzdOr9f7+/s/88wzbm5uR44csTUP4TqaHcMNDAzMzMx86aWXunXr1q1b\nt3nz5plMpgMHDtiah3DRSnzsBGzKzc1duXLl+PHjH3300RsZ/18xMTGenp7Hjx/Pz8+PiYlx\nc3Pz8/O76667cnJyGhoawsLC1FMsYHc3HtM1DAbDhx9+WFdXp9FohBBfffWVcsaLrfHrInqH\nsm+4eXl5xcXFytuNCCGeeOKJnTt37t+/Pykp6bqTEK5D2f2Zq9i7d2+/fv1aeBFHQbgOZd9w\nlYvKfH19lX08PDx0Ol0LL+UQrkM56JmraNeuXUhISFVVla1JCBetxBFCXN+pU6dWrFgxZ86c\na3612Rq/Lo1GYzAYTpw4cfLkSfVq5tjY2BMnTuTm5vJ6lePcVEzNWSyW7Ozs6upqX19fT09P\n5Ur0Hj162Bq3NQ/RO47dw7VarVartampSd3TbDa3MA/hOo7dw1W21tXV5eTk9OnT52/nIVzH\nsXu4Op0uKChIvW7w6tWrlZWVoaGhtuYhXMexe7jnz59///331V/FJpOpsrIyLCzM1jyEi1bi\nCKHsqqurLRaLcqmJ8uKTn5+fu7t7WlraiBEjIiIi1FekWhhXDixcV3x8/Pbt2ysrK2NjY5WR\nuLi4rKysq1evjhw50rHfmwRuKj6tVnvd/bVa7ffff28ymZKSkjw8PDZv3rx///7U1NTy8vIP\nPvjgoYceCggIEELYGreF6FvJaeFGR0cHBQVlZGRMmTLF29t7x44ddXV1vXv3bmFthNtKznzm\nCiGKioosFkt4ePiNrI1wW8mZ4SYlJX3xxRd6vV6v12/cuFGr1bb8rz/htpLTwnVzczt48KDZ\nbE5OTrZYLJmZmX5+fv369WthbYSL1uCD6WU3derUioqKa0YiIiIWLVp0zZ7Tp0/X6/XXHR8+\nfLit+SsrK1NSUiIjI9PS0pSRpqamiRMnWiyWDRs2qG9vgH/mpuIbPnz4dfcfMWLEqlWrLl++\nrHyIbVlZ2QcffFBQUKDVagcPHjxlyhQlJlvjthB9Kzkz3PPnz3/22WcFBQUWi+WOO+548skn\ne/bs2cLaCLeVnBmuEGLfvn2rV6/evHnzjURDuK3kzHCbmprWr1+/e/fuK1euREVFzZo1q+Uz\nAwm3lZwZ7tmzZ9etW1dYWOjl5RUVFZWamtrC4V9BuGgdCiEAAAAASIprCAEAAABAUhRCAAAA\nAJAUhRAAAAAAJEUhBAAAAABJUQgBAAAAQFIUQgAAAACQFIUQAAAAACRFIQQAuMCAAQOio6Md\n/SgJCQn2epTk5GQ/Pz+7TAUAwL+Hp6sXAACAoyQnJ//555/K7ePHj8fGxlqtVtcuCQCAfxUK\nIQCgzXrxxRfV29nZ2S5cCQAA/06cMgoAcL2dO3cOGjTI39+/Xbt2BoPhnXfeUQ/lDRo0aODA\ngTk5OUOGDAkICOjQocP48eMrKiqUrU1NTUuWLOnUqZNWq+3Vq9f333///PPPe3t7K1vVU0Yf\neeSR2bNnCyHc3Nx69+4thIiJiYmJiWm+hlGjRgUHByu3rVbrsmXLlGl79uz51VdfXbPgH374\nYejQoQEBAT4+PnFxcRkZGY760QAA4EgcIQQAuNi2bdsee+yxhx9+eP369X5+ft99993cuXON\nRuPKlSuFEN7e3gUFBdOnT3/zzTfvueee7OzscePGaTSaTz/9VAixfPnypUuXjh07NiUlpaSk\nZPLkyZ06dVILoeq999575ZVXtm/ffvjwYV9f379d0qpVqxYvXjxx4sQpU6b8/vvvS5cubWxs\nVLfu2bPn4Ycf7t+//4YNGzQazZYtW1JSUqqrq+fOnWvPnwsAAI7nxtUUAADnGzBgQFVV1enT\np4UQ3bt3r6+vLywsVIvc6NGjv/3224sXL7Zv3z4xMXHPnj379+/v37+/sjUxMfHXX38tLS21\nWq3h4eEhISEnTpxwc3MTQhw6dCghIcHX1/fKlStCiISEhJqaGuVRpk6d+sknn6h/9ZTDg8eP\nH1eXNGrUqP3791dVVVmtVr1er9Pp8vLylE0XL16MiIjw9vZWpo2Li6utrc3NzfXx8VF2GDly\n5N69eysqKrRarcN/dgAA2A+njAIAXKmsrOz06dPDhg1rflgvKSmpsbHxp59+Uu76+PiobVAI\nodfrjUajEMJoNJaXlw8dOlRpg0KIvn37GgyGVi6ppKSkrKzswQcfVEfCw8OVE02FEBUVFTk5\nOcOHD3d3dzf9ZdiwYbW1tWqBBADgVkEhBAC4UmlpqRCiY8eOzQfDw8OFEGVlZcrdkJCQ5ls9\nPT2bmpqEEOXl5erOqqioqFYuSWmb1zzo7bffrtxQVvXuu++2a2bGjBlCiAsXLrTyoQEAcDKu\nIQQAuJJycE8peCrlxE5397951bKhoeF/d1OPFv5j172YwmKxNL/7zDPPpKamXrPPnXfe2cqH\nBgDAySiEAABX0uv14q/jhCrlrrKpBTqdTvx1nFCVn59/I4/r7u7e/H1ixF8HBsVfxwbVu4pz\n584pN+644w4hhMViSUhIuJEHAgDg34xTRgEArhQWFmYwGHbs2GEymdTBLVu2+Pj43HfffS1/\nbZcuXQIDA3fu3KmOHD582NaFfMqRQ7PZrNwNCgoyGo3qwcCKiooTJ04otzt37hwcHJyVlaUe\ntywoKMjNzVVu63S6Pn36bNu2raamRp08MzPzP//5jzo5AAC3CgohAMDFVqxYYTQaR44c+fXX\nX2dlZc2aNSsrK2vRokUBAQEtf6Gnp2dKSsrJkyeffvrpXbt2rV27duzYsc3ffqY55SLAN998\nc/PmzUKIESNGVFVVrVixory8PCcnJzk5OTIyUtnT3d195syZZ86cGTNmzJYtW9asWfPQQw/F\nxcWpU61cubK+vn7w4MGZmZm7du1atGjR1KlTS0tLPT057wYAcKuxAgDgdP3794+KilLv7tq1\na8CAAb6+vhqNJjY2NiMjQ900ZMiQiIiI5l+bkpKi/v0ymUzPP/98cHCwr6/vwIEDDx06NGHC\nBD8/P2Vr37591UcpKSmJjY318vJSRhoaGubMmdOxY0eNRnPvvfd+8803zz77rL+/v7Kz2Wxe\nsGBBWFiYt7d3z549t27d+txzz3l7e6tryM7OHjp0qL+/v5eXV7du3VauXNnY2Gj3nxIAAI7G\n5xACANqUxMTEU6dOqe9QCgAAWsApowCAW1haWtrjjz+uXrxXU1Nz5MgR5UPnAQDA3+JqBwDA\nLax9+/ZbtmwZPXp0amqqyWRKS0u7fPny3LlzXb0uAABuDRRCAMAt7KmnnhJCrF69esKECVar\nNSYmZseOHUOGDHH1ugAAuDVwDSEAAAAASIprCAEAAABAUhRCAAAAAJAUhRAAAAAAJEUhBAAA\nAABJUQgBAAAAQFIUQgAAAACQFIUQAAAAACRFIQQAAAAASVEIAQAAAEBSFEIAAAAAkBSFEAAA\nAAAkRSEEAAAAAElRCAEAAABAUhRCAAAAAJAUhRAAAAAAJEUhBAAAAABJUQgBAAAAQFIUQgAA\nAACQFIUQAAAAACRFIQQAAAAASVEIAQAAAEBSFEIAAAAAkBSFEAAAAAAkRSEEAAAAAElRCAEA\nAABAUhRCAAAAAJAUhRAAAAAAJEUhBAAAAABJUQgBAAAAQFIUQgAAAACQFIUQAAAAACRFIQQA\nAAAASVEIAQAAAEBSFEIAAAAAkBSFEAAAAAAkRSEEAAAAAElRCAEAAABAUhRCAAAAAJAUhRAA\nAAAAJEUhBAAAAABJUQgBAAAAQFIUQgAAAACQFIUQAAAAACRFIQQAAAAASVEIAQAAAEBSFEIA\nAAAAkBSFEAAAAAAkRSEEAAAAAElRCAEAAABAUhRCAAAAAJAUhRAAAAAAJEUhBAAAAABJUQgB\nAAAAQFIUQgAAAACQFIUQAAAAACRFIQQAAAAASVEIAQAAAEBSFEIAAAAAkBSFEAAAAAAkRSEE\nAAAAAElRCAEAAABAUhRCAAAAAJAUhRAAAAAAJEUhBAAAAABJUQgBAAAAQFIUQgAAAACQFIUQ\nAAAAACRFIQQAAAAASVEIAQAAAEBSFEIAAAAAkBSFEAAAAAAkRSEEAAAAAElRCAEAAABAUhRC\nAAAAAJAUhRAAAAAAJEUhBAAAAABJUQgBAAAAQFIUQgAAAACQFIUQAAAAACRFIQQAAAAASVEI\nAQAAAEBSFEIAAAAAkBSFEAAAAAAkRSEEAAAAAElRCAEAAABAUhRCAAAAAJAUhRAAAAAAJEUh\nBAAAAABJUQgBAAAAQFIUQgAAAACQFIUQAAAAACRFIQQAAAAASVEIAQAAAEBSFEIAAAAAkBSF\nEAAAAAAkRSEEAAAAAElRCAEAAABAUhRCAAAAAJAUhRAAAAAAJEUhBAAAAABJUQgBAAAAQFIU\nQgAAAACQFIUQAAAAACT1XzDzimL0HB3ZAAAAAElFTkSuQmCC",
      "text/plain": [
       "plot without title"
      ]
     },
     "metadata": {
      "image/png": {
       "height": 480,
       "width": 600
      }
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#input: environmental dataframe, taxa, and metric of interest\n",
    "#output: map of location and metric values\n",
    "\n",
    "MapEnvironmentalMetric <- function(inputdataframe, taxa, metric) {\n",
    "\n",
    "minLon <- round(min(inputdataframe$longitude), 3)\n",
    "minLat <- round(min(inputdataframe$latitude), 3) \n",
    "maxLon <- round(max(inputdataframe$longitude), 3) \n",
    "maxLat <- round(max(inputdataframe$latitude), 3) \n",
    "    \n",
    "world <- ne_countries(scale = \"medium\", returnclass = \"sf\")\n",
    "\n",
    "MetricMap <- \n",
    "    ggplot(data=world)+\n",
    "    geom_sf() +\n",
    "    geom_point(data = inputdataframe, aes_string(x = \"longitude\", y = \"latitude\", color = metric)) +\n",
    "    coord_sf(xlim = c(minLon, maxLon), ylim = c(minLat, maxLat), expand = TRUE) +\n",
    "    labs(title = paste(\"Environmental Metric\", metric, \"for Taxa\", taxa)) +\n",
    "    theme(plot.title = element_text(hjust = 0.5))\n",
    "    \n",
    "    return(MetricMap)\n",
    "}\n",
    "\n",
    "SebastesMAP <- MapEnvironmentalMetric(Sebastes_SoquelCanyon, \"Sebastes\", \"sst\")\n",
    "\n",
    "SebastesMAP"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "R",
   "language": "R",
   "name": "ir"
  },
  "language_info": {
   "codemirror_mode": "r",
   "file_extension": ".r",
   "mimetype": "text/x-r-source",
   "name": "R",
   "pygments_lexer": "r",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}