{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Land Application"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import PFAS_SAT as ps\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from IPython.display import Image\n",
    "import pandas as pd\n",
    "pd.set_option('display.max_colwidth', 0)\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Model document\n",
    "\n",
    "Finished compost and dewatered, stabilized WWT solids are both suitable for land application. Once these materials are land applied, PFAS may volatize, be released to surface or ground- water, or sorbed to the soil/waste stream mixture. The flow of PFAS from land application is modeled using a liquid-solid partition coefficient normalized to the amount of organic carbon, combined with a water balance to track the flow of PFAS from the soil. Model predictions are based on achievement of equilibrium. By default, it is assumed that no volatilization occurs, but a user can enter a fraction of PFAS that volatilizes. It is further assumed that the land applied material is well mixed with the top layer of soil (thickness is a user input). The partition coefficient is used to estimate the concentration of PFAS in the liquid and solids. The concentration in the liquid changes throughout the year as PFAS runs off or is leached to the groundwater (i.e., it is assumed that annual precipitation is uniform throughout the year and continuously removes PFAS from the mixture). The user enters a run-off coefficient based on the soil type, land use, grade, and vegetation. The run-off is assumed to be released to surface water. Another fraction of the precipitation is removed via evapotranspiration (ET) based on the local climate and vegetation. The remaining precipitation is assumed to leach into groundwater. The PFAS remaining in the soil may be taken up by and bioaccumulate in plants. While PFAS uptake by plants was not modeled in this initial version of the SAT, it may be an important PFAS fate pathway as plants may enter the food chain depending on what is grown. \n",
    "\n",
    "<img src=\"../Images/ProcessModels/Land_Application_Diagram.png\" alt=\"Drawing\" style=\"width: 400px;\"/>\n",
    "\n",
    "### Assumptions and Limitations:\n",
    "1.\tThe organic carbon-normalized partition coefficient assumes that the organic carbon in the soil has the same PFAS sorption capacity as the organic carbon in either the compost or the dewatered, stabilized WWT solids in the land applied material.\n",
    "2.\tThe water balance model is averaged over a year and ignores potential effects from intense rains that may lead to substantial additional erosion and loss of solids and associated PFAS.\n",
    "3.\tApart from precipitation, the water balance does not consider other external water inputs such as irrigation. This could be included by adjusting the precipitation input value. \n",
    "4.\tVolatilization is assumed to be zero by default due to a lack of data. However, the user may assign a fraction of the PFAS to be volatilized.\n",
    "5.\tFuture work is also required to implement a dynamic model to account for changes in the organic C content of over time as land-applied materials decompose, and to account for episodic precipitation events. \n",
    "6.\tThe loading rate for the land application of dewatered WWT solids will vary based on solids properties among other factors. The loading rate default values given here reflect a typical annual loading rate for application to agricultural land growing corn. "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for Land Application model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Category</th>\n",
       "      <th>Dictonary_Name</th>\n",
       "      <th>Parameter Name</th>\n",
       "      <th>Parameter Description</th>\n",
       "      <th>amount</th>\n",
       "      <th>unit</th>\n",
       "      <th>minimum</th>\n",
       "      <th>maximum</th>\n",
       "      <th>Reference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log Koc (soils/sediments)</td>\n",
       "      <td>2.190</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>1.300</td>\n",
       "      <td>4.500</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log Koc  (soils/sediments)</td>\n",
       "      <td>3.040</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.400</td>\n",
       "      <td>4.700</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log Koc  (soils/sediments)</td>\n",
       "      <td>1.880</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>1.300</td>\n",
       "      <td>1.880</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log Koc (soils/sediments)</td>\n",
       "      <td>1.370</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log Koc (soils/sediments)</td>\n",
       "      <td>1.770</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>1.310</td>\n",
       "      <td>2.100</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log Koc (soils/sediments)</td>\n",
       "      <td>1.970</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2.190</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log Koc (soils/sediments)</td>\n",
       "      <td>2.630</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.300</td>\n",
       "      <td>3.180</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log Koc (soils/sediments)</td>\n",
       "      <td>3.240</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.650</td>\n",
       "      <td>3.780</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log Koc (soils/sediments)</td>\n",
       "      <td>1.510</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.790</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log Koc (soils/sediments)</td>\n",
       "      <td>2.790</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.050</td>\n",
       "      <td>2.875</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>Land application</td>\n",
       "      <td>LandApp</td>\n",
       "      <td>appl_dens</td>\n",
       "      <td>kg dry applied per area</td>\n",
       "      <td>2.000</td>\n",
       "      <td>kg TS/m2</td>\n",
       "      <td>0.800</td>\n",
       "      <td>18.000</td>\n",
       "      <td>[7]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Land application</td>\n",
       "      <td>LandApp</td>\n",
       "      <td>depth_mix</td>\n",
       "      <td>Depth of soil mixed</td>\n",
       "      <td>0.100</td>\n",
       "      <td>m</td>\n",
       "      <td>0.050</td>\n",
       "      <td>0.200</td>\n",
       "      <td>[8]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Soil Properties</td>\n",
       "      <td>SoilProp</td>\n",
       "      <td>bulk_dens</td>\n",
       "      <td>Soil bulk density - wet</td>\n",
       "      <td>1300.000</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td>1100.000</td>\n",
       "      <td>1800.000</td>\n",
       "      <td>[9]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>Soil Properties</td>\n",
       "      <td>SoilProp</td>\n",
       "      <td>ts_cont</td>\n",
       "      <td>Soil total solids content -wet</td>\n",
       "      <td>0.850</td>\n",
       "      <td>kg TS/kg</td>\n",
       "      <td>0.700</td>\n",
       "      <td>1.000</td>\n",
       "      <td>[10]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Soil Properties</td>\n",
       "      <td>SoilProp</td>\n",
       "      <td>C_cont</td>\n",
       "      <td>Soil organic C content -dry</td>\n",
       "      <td>0.015</td>\n",
       "      <td>fraction TS</td>\n",
       "      <td>0.005</td>\n",
       "      <td>0.050</td>\n",
       "      <td>[9]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Precipitation Data</td>\n",
       "      <td>Precip</td>\n",
       "      <td>ann_precip</td>\n",
       "      <td>Annual precipitation</td>\n",
       "      <td>1.000</td>\n",
       "      <td>m/yr</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>[11]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Precipitation Data</td>\n",
       "      <td>Precip</td>\n",
       "      <td>frac_runoff</td>\n",
       "      <td>Run-off fraction of precipitation</td>\n",
       "      <td>0.050</td>\n",
       "      <td>frac</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>[12]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Precipitation Data</td>\n",
       "      <td>Precip</td>\n",
       "      <td>frac_ET</td>\n",
       "      <td>Evapotranspiration fraction of precipitation</td>\n",
       "      <td>0.500</td>\n",
       "      <td>frac</td>\n",
       "      <td>0.033</td>\n",
       "      <td>1.000</td>\n",
       "      <td>[13]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>Volatilization</td>\n",
       "      <td>Volatilization</td>\n",
       "      <td>frac_vol_loss</td>\n",
       "      <td>Fraction of PFAS lost to volatilization</td>\n",
       "      <td>0.000</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                     Category  Dictonary_Name Parameter Name  \\\n",
       "0   Log partition coefficient  LogPartCoef     PFOA            \n",
       "1   Log partition coefficient  LogPartCoef     PFOS            \n",
       "2   Log partition coefficient  LogPartCoef     PFBA            \n",
       "3   Log partition coefficient  LogPartCoef     PFPeA           \n",
       "4   Log partition coefficient  LogPartCoef     PFHxA           \n",
       "5   Log partition coefficient  LogPartCoef     PFHpA           \n",
       "6   Log partition coefficient  LogPartCoef     PFNA            \n",
       "7   Log partition coefficient  LogPartCoef     PFDA            \n",
       "8   Log partition coefficient  LogPartCoef     PFBS            \n",
       "9   Log partition coefficient  LogPartCoef     PFHxS           \n",
       "10  Land application           LandApp         appl_dens       \n",
       "11  Land application           LandApp         depth_mix       \n",
       "12  Soil Properties            SoilProp        bulk_dens       \n",
       "13  Soil Properties            SoilProp        ts_cont         \n",
       "14  Soil Properties            SoilProp        C_cont          \n",
       "15  Precipitation Data         Precip          ann_precip      \n",
       "16  Precipitation Data         Precip          frac_runoff     \n",
       "17  Precipitation Data         Precip          frac_ET         \n",
       "18  Volatilization             Volatilization  frac_vol_loss   \n",
       "\n",
       "                           Parameter Description    amount         unit  \\\n",
       "0   PFOA Log Koc (soils/sediments)                2.190     log L/kg OC   \n",
       "1   PFOS Log Koc  (soils/sediments)               3.040     log L/kg OC   \n",
       "2   PFBA Log Koc  (soils/sediments)               1.880     log L/kg OC   \n",
       "3   PFPeA Log Koc (soils/sediments)               1.370     log L/kg OC   \n",
       "4   PFHxA Log Koc (soils/sediments)               1.770     log L/kg OC   \n",
       "5   PFHpA Log Koc (soils/sediments)               1.970     log L/kg OC   \n",
       "6   PFNA Log Koc (soils/sediments)                2.630     log L/kg OC   \n",
       "7   PFDA Log Koc (soils/sediments)                3.240     log L/kg OC   \n",
       "8   PFBS Log Koc (soils/sediments)                1.510     log L/kg OC   \n",
       "9   PFHxS Log Koc (soils/sediments)               2.790     log L/kg OC   \n",
       "10  kg dry applied per area                       2.000     kg TS/m2      \n",
       "11  Depth of soil mixed                           0.100     m             \n",
       "12  Soil bulk density - wet                       1300.000  kg/m3         \n",
       "13  Soil total solids content -wet                0.850     kg TS/kg      \n",
       "14  Soil organic C content -dry                   0.015     fraction TS   \n",
       "15  Annual precipitation                          1.000     m/yr          \n",
       "16  Run-off fraction of precipitation             0.050     frac          \n",
       "17  Evapotranspiration fraction of precipitation  0.500     frac          \n",
       "18  Fraction of PFAS lost to volatilization       0.000     fraction      \n",
       "\n",
       "     minimum   maximum  Reference  \n",
       "0   1.300     4.500     [2,3,4,5]  \n",
       "1   2.400     4.700     [2,3,4,5]  \n",
       "2   1.300     1.880     [2,3,4,5]  \n",
       "3  NaN       NaN        [2,3,4,5]  \n",
       "4   1.310     2.100     [2,3,4,5]  \n",
       "5  NaN        2.190     [2,3,4,5]  \n",
       "6   2.300     3.180     [2,3,4,5]  \n",
       "7   2.650     3.780     [2,3,4,5]  \n",
       "8  NaN        1.790     [2,3,4,5]  \n",
       "9   2.050     2.875     [2,3,4,5]  \n",
       "10  0.800     18.000    [7]        \n",
       "11  0.050     0.200     [8]        \n",
       "12  1100.000  1800.000  [9]        \n",
       "13  0.700     1.000     [10]       \n",
       "14  0.005     0.050     [9]        \n",
       "15 NaN       NaN        [11]       \n",
       "16 NaN       NaN        [12]       \n",
       "17  0.033     1.000     [13]       \n",
       "18 NaN       NaN        NaN        "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "LandApp = ps.LandApp()\n",
    "LandApp.InputData.Data[['Category','Dictonary_Name','Parameter Name', 'Parameter Description', 'amount', 'unit','minimum','maximum','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming Compost to land application"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Parameter</th>\n",
       "      <th>Unit</th>\n",
       "      <th>Amount</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Mass flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>1000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Solids flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>600</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>400</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>VS flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>360</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>168</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>5240</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>7375</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>3267</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>3424</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>16407</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>826</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1045</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>2710</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>7360</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>190</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Parameter Unit Amount\n",
       "0   Mass flow      kg   1000 \n",
       "1   Solids flow    kg   600  \n",
       "2   Moisture flow  kg   400  \n",
       "3   VS flow        kg   360  \n",
       "4   Carbon flow    kg   168  \n",
       "5   PFOA           μg   5240 \n",
       "6   PFOS           μg   7375 \n",
       "7   PFBA           μg   3267 \n",
       "8   PFPeA          μg   3424 \n",
       "9   PFHxA          μg   16407\n",
       "10  PFHpA          μg   826  \n",
       "11  PFNA           μg   1045 \n",
       "12  PFDA           μg   2710 \n",
       "13  PFBS           μg   7360 \n",
       "14  PFHxS          μg   190  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncominWaste = ps.IncomFlow()\n",
    "IncominWaste.set_flow('Compost', 1000)\n",
    "IncominWaste.calc()\n",
    "Compost = IncominWaste.Inc_flow\n",
    "Compost.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in land application"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Volatilized</th>\n",
       "      <th>Remaining in Soil</th>\n",
       "      <th>Leachate</th>\n",
       "      <th>Runoff</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>PFOA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>25.35</td>\n",
       "      <td>67.19</td>\n",
       "      <td>7.47</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>81.56</td>\n",
       "      <td>16.60</td>\n",
       "      <td>1.84</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>7.10</td>\n",
       "      <td>83.61</td>\n",
       "      <td>9.29</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.10</td>\n",
       "      <td>89.91</td>\n",
       "      <td>9.99</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>3.68</td>\n",
       "      <td>86.69</td>\n",
       "      <td>9.63</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>11.14</td>\n",
       "      <td>79.98</td>\n",
       "      <td>8.89</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>59.62</td>\n",
       "      <td>36.34</td>\n",
       "      <td>4.04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>87.89</td>\n",
       "      <td>10.90</td>\n",
       "      <td>1.21</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.45</td>\n",
       "      <td>89.59</td>\n",
       "      <td>9.95</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>69.75</td>\n",
       "      <td>27.22</td>\n",
       "      <td>3.02</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Volatilized  Remaining in Soil  Leachate  Runoff\n",
       "PFOA   0.0          25.35              67.19     7.47  \n",
       "PFOS   0.0          81.56              16.60     1.84  \n",
       "PFBA   0.0          7.10               83.61     9.29  \n",
       "PFPeA  0.0          0.10               89.91     9.99  \n",
       "PFHxA  0.0          3.68               86.69     9.63  \n",
       "PFHpA  0.0          11.14              79.98     8.89  \n",
       "PFNA   0.0          59.62              36.34     4.04  \n",
       "PFDA   0.0          87.89              10.90     1.21  \n",
       "PFBS   0.0          0.45               89.59     9.95  \n",
       "PFHxS  0.0          69.75              27.22     3.02  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "LandApp.calc(Inc_flow=Compost)\n",
    "LandApp.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deVxU9f7H8dfMIAgM4oILiGuoKDAMyJD7VrhDZipyzY0ol7LyXjPr3rI0zS3zupSamksqVHbTzH3LXUEZzUzUklQ0V5BFEBjm9wc/p8gNjOEM8Hk+HjwezVm+8zkTvjnzPef7PSqz2YwQQoiSp1a6ACGEKK8kgIUQQiESwEIIoRAJYCGEUIgEsBBCKEQCWAghFGJXlI3d3NzM9evXt1IpQghRNh05cuS62Wyu/tflRQrg+vXrExcXV3xVCSFEOaBSqX6733LpghBCCIVIAAshhEIkgIUQQiESwEIIoRAJYCGEUIgEsBBCKEQCWAghFCIBLIQQCpEAFkIIhUgACyGEQiSAhRBCIRLAQgihEAlgIYRQiASwEEIoRAJYCCEUIgEsikSj0aDX6/H19SU0NJSUlJRif487d+7w9NNPo9friYmJYc+ePfj4+KDX68nMzCz29xNCKRLAokgcHR0xGo2cOHGCqlWrMm/evGJ/j/j4eHJycjAajYSHh7Ny5UrGjBmD0WjE0dGx2N9PCKVIAIvH1rJlS5KSkgDYtWsXPXv2tKx75ZVXWLp0KZD/JJXx48cTGBiIn58fp06dAuDmzZv06tULnU5HixYtOH78OFevXuX555/HaDSi1+tZsGABX375JRMmTGDAgAElfoxCWJMEsHgsJpOJ7du3ExYWVqjt3dzcOHr0KCNGjGDGjBkAjB8/noCAAI4fP87kyZMZNGgQNWrUYNGiRbRt2xaj0ciwYcMICwtj+vTprFy50pqHJESJkwAWRZKZmYler6datWrcvHmTkJCQQu3Xu3dvAJo3b05iYiIAe/fuZeDAgQB06tSJGzducOvWLavULYQtkgAWRXK3D/i3334jOzvb0gdsZ2dHXl6eZbusrKwC+zk4OAD5F/Fyc3MBMJvN97SvUqmsVboQNkcCWDwWV1dXZs+ezYwZM8jJyaFevXqcPHmSO3fucOvWLbZv3/7INtq1a2fpVti1axdubm5UqlTJ2qULYTOK9Fh6If4sICAAf39/oqOjGThwIP369UOn09GoUSMCAgIeuf97773H0KFD0el0ODk5sWzZshKoWgjbobrf18AHCQoKMsfFxVmxHCGEKHtUKtURs9kc9Nfl0gUhhBAKkQAWQgiFSAALIYRCJIBFoW3atIkmTZrg5eXFlClT7llvNpt59dVX8fLyQqfTcfToUQASEhLQ6/WWn0qVKjFr1izLfnPmzKFJkyb4+PgwduxYAG7cuEHHjh3RarW88sorJXOAQpQwuQtCFIrJZOLll19m69ateHp6YjAYCAsLo1mzZpZtNm7cyJkzZzhz5gyHDh1ixIgRHDp0iCZNmmA0Gi3t1K5dm2effRaAnTt3snbtWo4fP46DgwNXr14FoGLFikycOJETJ05w4sSJkj9gIUqAnAGLQjl8+DBeXl40bNgQe3t7+vfvz9q1awtss3btWgYNGoRKpaJFixakpKRw+fLlAtts376dJ554gnr16gHw6aefMm7cOMtAjRo1agDg7OxMmzZtqFixYgkcnRDKkAAWhZKUlESdOnUsrz09PS0T8RRlm+joaCIiIiyvT58+zZ49e3jyySdp3749sbGxVjoCIWyPBLAolMIMG37UNtnZ2axbt46+fftaluXm5pKcnMzBgweZPn06/fr1u287QpRFEsCiUDw9Pblw4YLl9cWLF/Hw8CjSNhs3biQwMJCaNWsW2Kd3796oVCqCg4NRq9Vcv37dikcihO2QABaFYjAYOHPmDOfOnSM7O5vo6Oh7pqIMCwtj+fLlmM1mDh48iKurK+7u7pb1q1evLtD9ANCrVy927NgB5HdHZGdn4+bmZv0DEsIGyF0QolDs7OyYO3cuXbp0wWQyERkZiY+PD/Pnzwdg+PDhdO/enQ0bNuDl5YWTkxOff/65Zf/bt2+zdetWFixYUKDdyMhIIiMj8fX1xd7enmXLllm6LerXr09qairZ2dl8++23bNmypcBdF0KUdjIXhBBCWJnMBSGEEDZGAlgIIRQiASyEEAqRABZCCIVIAAuru337Nr2e60dTPz3nzp1TuhwhbIYEsLCqS5cuYWjZhj2/3ORaTQOBhifZt2+f0mUJYRMkgIXVHD16FF1gEFeq6tB2eR1tUC/sO71Clx5hLFu2XOnyhFCcBLCwiitXrtCqdRvS75hwaNLWMriiYl0darf6DBkymG+/XfuIVoQo2ySAhVXUqFGDFcuXoW/WmDsX/pjP13Q7GfO1X1m9ejVt2rRWsEIhlCcBLKxCpVLRt29f6tStd886B0cn+vfvL3M+iHJPAlgIIRQiASyEEAqRABZCCIVIAAshhEIkgIUQQiESwEIIoRAJYCGEUIgEsBBCKEQCWAghFCIBLIQQCpEAFkIIhUgACyGEQiSAhVXl5eUVXGBWpg4hbJEEsLCqGzdvorZ3srw259zBwcFBwYqEsB0SwMJqkpOTOXzwAA6ezSzLsq//RuMmTRSsSgjbIQEsrGb+/AU4PRGExrmyZZnp9wTatWqhYFVC2A4JYGEV33//PR9MnY6Doa9lmdmUw51Tu4noH65gZULYDglgUazi4uKIGDiYfgMGUyn0LSpUq2NZl3FiOzo/X7y9vRWsUAjbYad0AeXR6/96g8WLF//tdtJvJRdDNQ+nda1SqO3MZjPZWZlUdK1GBZ/OVB00G42Tq2V9buo1MvevZOHeH6xVqhCljgSwAnp07cyCBQuo/Ox47Kq4P3Y7LhnJ5N3JKMbKCrKr7I5KrSn09qoKDqjsHCxPQL7LlJVO+neTeXvcWPz8/Iq7TCFKLQlgBYSEhDBj2hTenjgN1/CpqCtqH6udP59h2qrsq+fI2DyTQX178fa4N5Uup1Qwm833/BETZZP0ASvk5ZEjiXg2lPSNMzDnmZQup9jl3rpK6u7PSftmPB+9/29mf/yRhEohGI1G6jV4gl9++UXpUkQJUJnNhR+aFBQUZI6Li7NiOeVLbm4uT3fpzom0img7DVO6nL8lLzuTy0tfI/fWFfj/PygVKlRgzJgxeHp6Klxd6ZGQkMC8TxfgWacuxiOHqVy58qN3EjZPpVIdMZvNQfcslwBW1q1bt9AHPUl6w044B/RQupzHlncngwuz8m8vc3DzxNnTG3WFigpXVTrZ1dWTc/44TZ0z2Ll1E3Z20lNY2kkA27Bff/2V5sEtsH9qFI4NApUuR9gAc56JtLUf0KdTEAs//UTpcsTf9KAAlj5gG9CwYUPW/W8NGZtnkXP9gtLlCBugUmvQdh9DzNqNzJk7T+lyhJVIANuItm3bMufjj0j7bhKmzFSlyxE2QO3gjDb037z1zni2bNmidDnCCiSAbcjQoUOIGtif9PVTMZtylC5H2IAKVdzRdn+DPv3/walTp5QuRxQzCWAbM2PaVJ70rkv6jvkUpX9elF0V6/hi33IgT3Xpxo0bN5QuRxQjCWAbo1ar+TpmFdUyk7gd963S5QgbofV7mju1g+gW2ovs7GylyxHFRALYBmm1WrZv3kDej+u5ffaQ0uUIG+HcZiBnU/KIfHGYfDsqIySAbVSdOnXYuH4dt7fNI/vqOaXLETZApdag7Tqa9Tv2MnX6DKXLEcVAAtiGBQcH89mn80j7bjKmDOvPfCZsn9reEW3o23wwZTrr1q1TuhzxN0kA27iIiP68NuJF0tdPwZwrfX8C7CrVwKXnOAYMGsLx48fRah9vMqfCqF+/PtevXy/09kajkQ0bNlitnrJGArgUmPj+e7QLbEr61jnS9ycAcPBoQsV2UYR062FTvxMSwEUjAVwKqFQqor9Yjoc6lYyDMUqXI2yEc7P2mJ5oR1ZWFllZWQXWXbt2jeeeew6DwYDBYGDfvn0AHD58mFatWhEQEECrVq1ISEgAwGQyMWbMGPz8/NDpdMyZM8fS1pw5cwgMDMTPz89yL/L92snOzubdd98lJiYGvV5PTEwMGRkZREZGYjAYCAgIYO3atSX06ZQOEsClhKOjI1s2fIf6zE4yTu1VuhxhI5xb9icPGDBoSIEz4ddee43Ro0cTGxvLmjVriIqKAsDb25vdu3cTHx/PhAkTePvttwFYuHAh586dIz4+nuPHjzNgwABLW25ubhw9epQRI0YwY8aMB7Zjb2/PhAkTCA8Px2g0Eh4ezqRJk+jUqROxsbHs3LmTN954g4wM6z1EoLSRaZZKEXd3d7ZsWE/bjk9h51oTB/dGSpckFKZSqVFp7Nl+6DjvTZjI++PfBWDbtm2cPHnSsl1qaippaWncunWLwYMHc+bMGVQqFTk5OZbthw8fbpl5rWrVqpZ9e/fuDUDz5s355ptvAB7Yzl9t2bKFdevWWYI7KyuL8+fP07Rp02L+JEonCeC/0Gq1pKenl/j7vvvuu7Rr146nn376odvp9XqWL1nEwKjhaMKnYVfJ7Z5tLn4aidreEVQqNM5VcOvxTzTaKn8sV+d/8akaMpKKnvn/EFJjvyX5h2XUGfUFagdnAPJysrixaQ45VxMBM2oHLTX6vZ/fhrAp2tC3mDnnDZp5exMe3o+8vDwOHDiAo2PB/1ejRo2iY8eO/O9//yMxMZEOHToAD38Kh4ODAwAajYbc3FwA3nnnnfu281dms5k1a9bQpEmT4jnQMka6IGzEhAkTHhm+dz377LOM+9frpH83mbzsrPtuUzNiMh6Rc7Gv5cWtg18WXD50Dh5D51jCFyDj5x9wcG/E7dMHLMvS4tahcaqMxwvz8HjhE6p1exWVWv5m2yI7bVVcQt/mhWHDiYuLo3PnzsydO9ey3mg0AvlnrrVr1wZg6dKllvWdO3dm/vz5loC9efPmQ9/vQe24uLiQlpZmed2lSxfmzPnj4nF8fPzjH2QZJAH8ALt27aJDhw706dMHb29vBgwYYPklio2NpVWrVvj7+xMcHExaWhpZWVkMHToUPz8/AgIC2LlzJ5D/y9mrVy9CQ0Np0KABc+fOZebMmQQEBNCiRQvLL/qQIUP4+uuvgfxbf8aPH3/PhY9r164REhJCYGAg5xPPkXMzibQN0zGb8x54HBXr+JKTfOmhx5qTfBlzdhaV2w4k4+QfTy02ZSRj51LN8rpCNU9UdhUe49MU1mTOucPFeYO5+vUEMu/k0r5jJ8aNG0dcXBw6nY5mzZoxf/58AMaOHctbb71F69atMZn+eBRWVFQUdevWRafT4e/vz6pVqx76ng9qp2PHjpw8edJyEe6dd94hJycHnU6Hr68v77zzjnU+hFJKTmceIj4+np9++gkPDw9at27Nvn37CA4OJjw8nJiYGAwGA6mpqTg6OvLf//4XgB9//JFTp07RuXNnTp8+DcCJEyeIj48nKysLLy8vpk6dSnx8PKNHj2b58uW8/vrr97z33Qsfn3zyCTNmzGDRokW8//77dOrUibfeeotNmzbx2Wef0cg5j4v7VqJtM/C+x5D5y2Hsq9e3vL6y+m1Qq1FpKuA+aCaQf/br1LQdDnV8yPk+CVNGChrnymj9Qrjy5TtkJOyjYj09Wt9OVKhau5g/ZfF31XvzuwKv0w9/TfiAQcQd3Iezs3OBdS1btrT8XgJMnDgRADs7O2bOnMnMmTMLbJ+YmGj576CgIHbt2vXQdqpWrUpsbGyBNhYsWPB4B1YOyBnwQwQHB+Pp6YlarUav15OYmEhCQgLu7u4YDAYAKlWqhJ2dHXv37mXgwPwQ9Pb2pl69epZf0I4dO+Li4kL16tVxdXUlNDQUAD8/vwK/4H/25wsfd7fZu3cv/fv3B6Br165UqVKFmFUrsD9/kIyfdhTY/8rqt7n0+Sjy7tzGtUVfy/K7XRB3wxfg9s+7cW7aDpVKjVPjltxOyL/Lwr5mQ2oPW4Rr8HPkZaVxefk/ZcL4UsDZ8BxX7arTt/8A8vIe/O1IKE/OgB/i7sUH+OMCxIMuVjzsZvg/t6NWqy2v1Wq1pc/tQfv8+cLH/d6jWrVqbNu0gRat26JxrUVFz2ZAftAW5rH12VfPkZN8iSsx///VMC8XO9dauAT2zK/R3hGnJq1watIKlUpF5q+xVHCr88h2hXJUKhUVg8PZvHgEly9ftvTVCtsjZ8BF5O3tzaVLlyxfs9LS0sjNzaVdu3asXLkSgNOnT3P+/Pliv/Lbpk0bvvwy/4Lali1bSE7Onx+iWbNmxKz6gvTvp+U/lbgIMn7+gcqt/4HniCX5Py8vJzf9Brm3rpJ18SSmrPw7QsymHLKvX0BTqUaxHpMofnnZmWRsmMYHkyZJ+No4CeAisre3JyYmhlGjRuHv709ISAhZWVmMHDkSk8mEn58f4eHhLF26tMCZb3EYPz7/0TSBgYFs3LgRd3d3XFxcAOjWrRsT3v0Paes+gCIMTc34eQ9OjVsWWObUqAUZP+8mN+UyV1aN49Lil7n8+Ws41PLCqUnrYj0mUbzM5jwyNs+ie4eW6HV+NGnSBC8vL6ZMmXLf7Xft2oVer8fHx4f27dtbln/88cf4+Pjg6+tLRETEPSPtZsyYgUqlsswTsXLlSvR6veVHrVZb7rwQDyZPRS5F7ty5g0ajwc7OjgMHDjBixIgCv+Rms5nIF4exdt9xXELfRqXWKFitUEL6nmU0yLvE7h1b8fX1ZevWrXh6emIwGFi9ejXNmjWzbJuSkkKrVq3YtGkTdevW5erVq9SoUYOkpCTatGnDyZMncXR0pF+/fnTv3p0hQ4YAcOHCBaKiojh16hRHjhzBza3gveg//vgjzzzzDL/++mtJHrpNk6cilwHnz5/HYDDg7+/Pq6++ymeffVZgvUqlYuGn82ha3ZGMPcsUqlIoJf3HbThciGXjd99iNBrx8vKiYcOG2Nvb079//3vmYVi1ahW9e/embt26ANSo8Uf3Um5uLpmZmeTm5nL79m08PDws60aPHs20adMeOHBj9erVREREWOEIyx4J4FKkUaNGxMfHc+zYMWJjYy13YvxZhQoVCO3ZnbzfExSoUCgl6+JPZO9fzrbNG3BzcyMpKYk6df64WOrp6UlSUlKBfU6fPk1ycjIdOnSgefPmLF++HIDatWszZswY6tati7u7O66urnTu3BmAdevWUbt2bfz9/R9YS0xMjARwIUkA26jIyEhq1KiBr69vkfbbtGkTH3w4Dccuo61UmbA1OSm/k/79dL6KXmWZY+F+XYt/PWPNzc3lyJEjfP/992zevJmJEydaQnnt2rWcO3eOS5cukZGRwRdffMHt27eZNGkSEyZMeGAthw4dwsnJqci/t+WVBLCNGjJkCJs2bSrSPidPnqRfxAC0Pd6gQuVaVqpM2JK8Oxmkr5vEB++9S5cuXSzLPT09uXDhj3u2L168WKAb4e42Xbt2xdnZGTc3N9q1a8exY8fYtm0bDRo0oHr16lSoUIHevXuzf/9+fvnlF86dO4e/vz/169fn4sWLBAYG8vvvv1vajI6OlrPfIpAAtlHt2rUrMCPVo1y/fp2nu3bHofVgKnr6WLEyYSvMeSbSN35E356def21UQXWGQwGzpw5w7lz58jOziY6OpqwsLAC2zzzzDPs2bPH0s976NAhmjZtSt26dTl48CC3b9/GbDazfft2mjZtip+fH1evXiUxMZHExEQ8PT05evQotWrl/7HPy8vjq6++sgwWEo8mAzHKgOzsbLr2fIbsOsFofZ9SuhxRQjJ+WEKzWlrmfzL3nnV2dnbMnTuXLl26YDKZiIyMxMfHxzInxPDhw2natCldu3ZFp9OhVquJioqydB306dOHwMBA7OzsCAgI4KWXXnpkPbt378bT05OGDRsW74GWYXIbmg1LTEykZ8+enDhx4oHbmM1mBgwawuajv6Dt+SYqlXypKQ8yjBtxPrMF45HDVK5cWelyxCM86DY0OQMu5T6cOo0Nuw5Qqe9kCd9yIjPRSE5sDNsOHZDwLeXkX2wptnbtWiZP+wht6NsySXo5kXPjIhmbZvLtmq/w8vJSuhzxN0kA26iIiAhatmxJQkICnp6eLF68uMB6o9HIgMFDcek5DrtK1RWqUpQkU2Yaad9NYub0KQ98AoUoXaQLwkatXr36get+//13OnfviWP7F3HwkEe9lAdmUy4ZG6YxuH8fXnrxRaXLEcVEzoBLmczMTDp3DyWvUQecm7ZTuhxRAsxmMxk7FxDYsBazPpqudDmiGMkZcCly946HJJMWbQu517K8uH10HVXSf+N/Ow+i0cgES2WJBHAp8u5777Mz9kdcnvvggROhiLIl85dYTMa1bIs7bJl6VJQdEsClRHR0DLPmLcC1/zTUFYp3nmFhm7KvJZKxdQ5bN66nXr16SpcjrED6gEuB2NhYooaPwCX0bey0hR+eLEovU0YKaesm8+mc/9KyZctH7yBKJTkDtnEXL16ka88wnJ56GfuaMsTz70g/voUbG2cXe7uutYr3GXkOTzxJ9qVTjIwawsCBzxdr28K2SADbsIyMDJ7u2gOVTzecGrVQupxSz7FxKyrfvkXOzUvkZaWRfeVXTKlXLes1Gg3Dhw+nUaNGhW7Ty8uLxo0bF1uNBw8eZNCgQfR45lk+nDSx2NoVtkkC2Ebl5eXRJ/wfXLd3x9nQW+lyygRNRS2uLfred11Oyu/cOb6Jz1esQh8QwIwPP1Dkq7+9vT1DX3iRubNnoVZLD2FZJ/+HbdTYcW9z8OdEnDsNlzseSkCFyrXQthuC24uLSXDyI6R7GKNe/yeZmZklWke9evVYsmghTk5OJfq+QhkSwDZo6dJlLFz6Bdoeb6Kyq6B0OeWKys4erS6EKgP/y+od8TRp5sehQ4eULkuUURLANmbv3r288vo/0Yb9G42Tq9LllFsaJ1e03ceQ6d+Xpzp3Y/369UqXJMogCWAbcu7cOUJ79ca582vYu9VVuhwBOHu3weWZ/9D/+cGsW7dO6XJEGSMBbCNSU1N5qkt37AJ749iwudLliD9x8GiCyzP/4R+DhrJnzx6lyxFliASwDTCZTPR6rh+plb1wCuipdDniPhzcG+McMoo+/f9Bamqq0uWIMkIC2AaMeu2fGH+7hnP7KLnjwYY5PmHA5O7HqNf/qXQpooyQAFbYp/Pn88XX3+LcfSwqjdyWbeuc2g7hm3Xr2bp1q9KliDJAAlhB27dv5423/oNL2L/RVNQqXY4oBLWDM47tX2LEq6MpygNthbgfCWCFnD59mt59w9F2+xcVqtZWuhxRBI5eBq4lp3LgwAGlSxGlnASwApKTk3m6S3fsW/yDinV1SpcjikilUmPn05mP/jtH6VJEKSedjgr419hxXL58maou8WRejH/sdq7+uLcYq7qX1r0BTm6FOzvPuZ1O8i9GACpWqUklz8JPaFOy1OQ5aFE3bEnF+v6oVI93DuLk+zQbFr/E9evXcXNzK+YaRXkhAayAt998g9DuXf92O59/XoWdO3cWQ0X3998P/kOVKlUKtW1GRgYDBw4EICv5ClnJV1iyZAmVK1e2Wn2Pw2QyceHCBeYtWMS1Hxbi+NTLVKzjW+R2NI4uuNRtyv79+wkLC7NCpaI8UBXlQkJQUJA5Li7OiuWIsuDChQvUrZs/ki8tLQ2t1vYuMJrNZtavX88/Bg3BuesYKtYreldQ6p4VjOjYiEkfyLSR4uFUKtURs9kc9Nfl0gcsil2dOnW4c+cOAO7u7gpXc38qlYrQ0FC+/fpL0jZ+hDk3u8ht2NV8gj0HDluhOlFeSAALq7C3t2fRokWkp6eTnV30cCspTz31FHq9PxkJ+4q8r30tL348brRCVaK8kAAWVvPCCy8AMHGibX9FH/3KCNRndxd5P41TZdJvpVihIlFeSAALq7P16waNGzfGlH6z6DtqKmDKzSEvL6/4ixLlggSwsLqSfqpEUdnb25Nnyi3yfiqVClQqGREnHpsEsBBCKEQCWAghFCIBLIQQCpEAFkIIhUgACyGEQiSAhRBCIRLAQgihEAlgIYRQiASwEEIoRAJYCCEUIgEshBAKkQAWQgiFSAALIYRCJIBFufbi8JH46wNQO1d9rP0rVqqGs0slpn80s5grE+WBBLAo1TQaDXq9Hl9fX0JDQ0lJKdoE6SNeisLFtTK4Ny3UtJJ3Lp/h5rYFAJhNOVSs50+deg0I79unwHatWrUqUh2ifJIAFqWao6MjRqOREydOULVqVebNm1ek/QMDAzkef4Rat34mffOsRz4bzsG9EVWfHoYpM5W0/71P81r2HD18wPIQ0rv2799f5GMR5Y8EsCgzWrZsSVJSkuX19OnTMRgM6HQ6xo8fD0BiYiLe3t5ERUXh6+vLgAEDOHnyJJWcHMg6e5Cbq/6FKSOZO5cS+H3FGC59/iq/rxhDzo2LAGSdP87vq8ZxK3osjapXxNO9BqGhoTRs2JDZs2db3vvuk6B37dpFhw4d6NOnD97e3gwYMMBypr1hwwa8vb1p06YNr776Kj179iypj0rYCAlgUSaYTCa2b99OWFgYAFu2bOHMmTMcPnwYo9HIkSNH2L07/7lvZ8+e5bXXXuP48eOcOnWKVatWsX//fr5YsQJ3bQVSosdiNudRc8BUPIbOxrXt86TsXg7AnctnybmcwEeTxtOtSxcSEhLYvHkzhw8f5v333ycnJ+ee2uLj45k1axYnT57k119/Zd++fWRlZTFs2DA2btzI3r17uXbtWsl9WMJmSACLUi0zMxO9Xk+1atW4efMmISEhQH4Ab9myhYCAAAIDAzl16hRnzpwBoEGDBvj5+aFWq/Hx8eGpp55CpVLh7++Pk5Mjn348nVvffsDvK8dyafFIkrd/Rvb186THf0923NcEBzXnxagoAHr06IGDgwNubm7UqFGDK1eu3FNjcHAwnp6eqNVq9Ho9iYmJnDp1ioYNG9KgQQMAIiIiSugTE7ZEAliUanf7gH/77Teys7MtfcBms5m33noLo9GI0Wjk7Nmzlqc0Ozg4WPZXq9WW12q1mtzcXJ5/fgCd2rVGcyuJSv6dqd773+Sl30B7dhsL539CtWrVLPv/uS2NRkNu7r3PlrvfNphtsscAAAp+SURBVPIcOQESwKKMcHV1Zfbs2cyYMYOcnBy6dOnCkiVLSE9PByApKYmrV68Wuj17e3tmzphOlYv7ubrin1RQmTEeOUzt2rWLpV5vb29+/fVXEhMTAYiJiSmWdkXpYqd0AUIUl4CAAPz9/YmOjmbgwIH8/PPPtGzZEsi/KPbFF1+g0WgK1dbYsWMZPHgw1apWxa9ZCEePHqVy5crFVqujoyOffPIJXbt2xc3NjeDg4GJrW5QeqqJ8FQoKCjLHxcVZsRxR1qhUKtq3b8+uXbuULsXmpKeno9VqMZvNvPzyyzRq1IjRo0crXZawApVKdcRsNgf9dbl0QQihkM8++wy9Xo+Pjw+3bt1i2LBhSpckSph0QQihkNGjR8sZbzknZ8BC3IcMJRYlQQJYiPuQocSiJEgAC3EfhRlKLMTfJQEsxCPcbyixEMVBAliIR7jfUGIhioMEsBCPUJjhxkI8DglgIYRQiASwEEIoRAZiCHEfdyfx6dChAx06dLAsnzt3rkIVibJIzoCFEEIhEsBCCKEQCWAhhFCIBLAQwKRJk/Dx8UGn06HX6zl06NADt42KiuLkyZMA1K9fn+vXr5dUmaKMkYtwotw7cOAA69ev5+jRozg4OHD9+nWysx/8ePpFixaVYHWiLJMzYGEVZrOZTZs2AXDs2DFiY2MVrujBLl++jJubm2XAhZubGx4eHmzfvp2AgAD8/PyIjIzkzp07QP6dEfJgAlEcJICFVVy9epXQsGcASElJYdDQFxSu6ME6d+7MhQsXaNy4MSNHjuSHH34gKyuLIUOGEBMTw48//khubi6ffvqp0qWKMkYCWFhFzZo12bVzB5WqVMPHP5AdWzcrXdIDabVajhw5wsKFC6levTrh4eEsWLCABg0a0LhxYwAGDx7M7t27Fa5UlDXSByyspnXr1ly68BsVKlTA3t5e6XIeSqPRWAZd+Pn5sWzZMqVLEuWAnAELq3J2drb58E1ISODMmTOW10ajkZo1a5KYmMjZs2cBWLFiBe3bt1eqRFFGyRmwKPfS09MZNWoUKSkp2NnZ4eXlxcKFC4mIiKBv377k5uZiMBgYPny40qWKMkYeSy+EEFYmj6UX9+jQoQObNxe8ODZr1ixGjhx53+0TExPx9fV9aJuJiYmsWrXK8jouLo5XX30VgKVLl/LKK68AMH/+fJYvX/53ygdkIIQo3SSAy7GIiAiio6MLLIuOjiYiIuKx2/xrAAcFBTF79ux7ths+fDiDBg167PcRoiyQAC7H+vTpw/r16y0DDBITE7l06RJt2rThjTfewNfXFz8/P2JiYu7ZNzExkbZt2xIYGEhgYKDlKcLjxo1jz5496PV6Pv74Y3bt2kXPnj3v2f+9995jxowZXLp0Cb1eb/nRaDT89ttvXLt2jeeeew6DwYDBYLA8h+3GjRt07tyZgIAAhg0bJg/IFKWaBHA5Vq1aNYKDgy0j1qKjowkPD+ebb77BaDRy7Ngxtm3bxhtvvMHly5cL7FujRg22bt3K0aNHiYmJsXQzTJkyhbZt22I0Ghk9evQja/Dw8MBoNGI0GnnxxRd57rnnqFevHq+99hqjR48mNjaWNWvWEBUVBcD7779PmzZtiI+PJywsjPPnzz/yPTZt2kSTJk3w8vJiypQp96xfuXIlOp0OnU5Hq1atOHbsGABZWVkEBwfj7++Pj48P48ePt+xz7NgxWrZsiZ+fH6GhoaSmpgKwdetWmjdvjp+fH82bN2fHjh2PrE+UX3IXRDl3txvimWeeITo6miVLlrBixQoiIiLQaDTUrFmT9u3bExsbi06ns+yXk5PDK6+8gtFoRKPRcPr06b9Vx759+1i0aBF79uwBYNu2bZYJbwBSU1NJS0tj9+7dfPPNNwD06NGDKlWqPLRdk8nEyy+/zNatW/H09MRgMBAWFkazZs0s2zRo0IAffviBKlWqsHHjRl566SUOHTqEg4MDO3bsQKvVkpOTQ5s2bejWrRstWrQgKiqKGTNm0L59e5YsWcL06dOZOHEibm5ufPfdd3h4eHDixAm6dOlCUlLS3/psRNklZ8DlXK9evdi+fTtHjx4lMzOTwMDAQn2t//jjj6lZsybHjh0jLi7uoZPXPMrly5d54YUXiImJQavVApCXl8eBAwcsZ8dJSUm4uLgAoFKpCt324cOH8fLyomHDhtjb29O/f3/Wrl1bYJtWrVpZgrxFixZcvHjR8j5368nJySEnJ8fy3gkJCbRr1w6AkJAQ1qxZA0BAQAAeHh4A+Pj4kJWVZeniEeKvJIDLOa1WS4cOHYiMjLRcfGvXrh0xMTGYTCauXbvG7t27CQ4OLrDfrVu3cHd3R61Ws2LFCkwmEwAuLi6kpaUV+v1zcnLo168fU6dOtQz7hfz5Gf78+B+j0WipbeXKlQBs3LiR5OTkh7aflJREnTp1LK89PT0feka6ePFiunXrZnltMpnQ6/XUqFGDkJAQnnzySQB8fX1Zt24dAF999RUXLly4p601a9YQEBBQ4KnKQvyZBLAgIiKCY8eO0b9/fwCeffZZdDod/v7+dOrUiWnTplGrVq0C+4wcOZJly5bRokULTp8+jbOzMwA6nQ47Ozv8/f35+OOPH/ne+/fvJzY2lvHjx1suxF26dInZs2cTFxeHTqejWbNmzJ8/H4Dx48eze/duAgMD2bJlC3Xr1n1o+/c7m3/QGfTOnTtZvHgxU6dOtSzTaDQYjUYuXrzI4cOHOXHiBABLlixh3rx5NG/enLS0tHtG+/3000+8+eabLFiw4JGfgSjHzGZzoX+aN29uFqI02b9/v7lz586W15MnTzZPnjz5nu2OHTtmbtiwoTkhIeGBbb333nvm6dOn37M8ISHBbDAYLK8vXLhgbtSokXnv3r1/s3pRVgBx5vtkqpwBizLNYDBw5swZzp07R3Z2NtHR0YSFhRXY5vz58/Tu3ZsVK1YU6Aa5du0aKSkpAGRmZrJt2za8vb2B/Ok2Ib+v+oMPPrAMU05JSaFHjx58+OGHtG7duiQOUZRiEsCiTLOzs2Pu3Ll06dKFpk2b0q9fP3x8fJg/f76lW2PChAncuHGDkSNHotfrCQrKHzF6+fJlOnbsiE6nw2AwEBISYrmnefXq1TRu3Bhvb288PDwYOnQokP/Y+rNnzzJx4kRLl8rdsBbir2QuCCGEsDKZC0KUuEcNgBCivJMAFlZxdwDExo0bOXnyJKtXry4wsEIIIQEsrKQwAyCEKO8kgIVVFHUAhBDlkQSwsIr7XdwtyhBiIcoDCWBhFZ6engWG5168eNEyR4IQIp8EsLCKwgyAEKK8k+kohVX8eQCEyWQiMjISHx8fpcsSwqbIQAwhhLAyGYghhBA2RgJYCCEUIgEshBAKkQAWQgiFSAALIYRCJICFEEIhEsBCCKEQCWAhhFCIBLAQQihEAlgIIRQiASyEEAqRABZCCIVIAAshhEIkgIUQQiESwEIIoRAJYCGEUIgEsBBCKKRIT8RQqVTXgN+sV44QQpRJ9cxmc/W/LixSAAshhCg+0gUhhBAKkQAWQgiFSAALIYRCJICFEEIhEsBCCKEQCWAhhFCIBLAQQihEAlgIIRQiASyEEAr5P/qP0I41Zv7wAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "LandApp.plot_sankey()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdeVxU1fvA8c8Z9l3NNPddQQTBQHNJcd+QctfMJbVcM8uybLPsp1mpuWdmavo1t8xccU0tMxdUNMWtlFTcNwRk5/7+mIFAAQGBy8Dzfr181cyce+9zZ4ZnnnvuuecqTdMQQgghhBDCHBj0DkAIIYQQQoiskuJVCCGEEEKYDSlehRBCCCGE2ZDiVQghhBBCmA0pXoUQQgghhNmQ4lUIIYQQQpgNKV5zmVLKVimlKaXK6x1LViilaiql7ukdhxBC6EXythDmpUgUr0qpyFT/kpRS0ake93nMsu2UUn/nYiz7lVIv59b6npSmaWc1TSuWF+tWSq1QSsWa3uc7SqktSqkaptcmK6XiH/psRqVatrNSKkgpFaWUuqWUWqKUKpPONtqZfnRGPfzaY2JzNS2XevtjM2hrr5RapJS6qJS6r5Q6rJRqnUHbz03rbZLB+5D872B24jWtx1cpFayUeqCUOqiUqpPqtX5KqbNKqXCl1HWl1PdKKYfsbkOIgkLydsbMNW+bDhJmKqXCTMueV0p9kcW4aj603UhTrh2RQft3lVIXTDk7TCn1pVLKItXrPkqpfaaceUkp9V6q14YqpRLS2d5T2Xwvn1ZKbTC9HxeUUt1SveallDqilLprep+3KqVqZmf9RVmRKF41TXNM/gdcBDqlem6Z3vEVcp+Z3veKQATwXarXfkj92WiaNhPA9MO0CPgSeArwBCyA35RSzg+tvz9wx/Tf7Ep8aPtfZtDOBvgHaAIUAyYCa5RSZVM3Ukq5Ae2B2+ms47OHtlU/vQ2ZfiSfS+d5O2AdMB8oDqwG1iqlLE1N9gANNU1zAaoDzsD4TPdeiAJM8rau8ipvjwfcgHqAE9AKOJ6VgEwFe+rvhA+QCKzNYJE1QF1N05yBukAjYEiq11cBgUAJUxxjlFJtUr2++6H9dNQ07ZHcbiro33v4eZP5wF2gFDAIWJh8IIDxO93ZtP1SwA7gf495G4RJkSheH0cpZaeUmqOUuqqUuqyU+kopZWU6yloLVE195KWUaqyUOmA6YruilPo6VRGRne22U0r9rZR6Xyl103R02CfV6w6mo9RLpm3tSd6OUqqrUipEKXVPKbUj1R8ESqlrSqm3lFInTTF/o5Qqo5TabjoK3ZKcTJSxBzIh1bL7lVLjTf+9r5TarJQqnur1wcrYA3lTKTXWtK2UXsaMaJoWCawA6mTWzrR/U4CPNU1bpWlajKZpV4B+piYjU7V1AV4AhgJ1VaqeyNykadpdTdP+T9O0i5qmJWma9jNwDfB+qOlcYAyQ8MhKnlxrIEbTtLmapsUCUzEm/yamGP9NlVgVkISxiH3Ew5+56bmUniWllKXpe3dbKfWPUmrUw+2F0JvkbbPM277AGk3TrmtG55/gQKQ/sN20nfRi/1vTtPupnkrJicrYA1sBWKZpWqKmaWeA/YB7DmN5hOn97wR8pGlalKZpvwJbgT6m+O6Y8nbybU4TySBnm9aX5jNTxqJ5QarHOfqMzZUUr0afYjxK9ACeBfyAsaZioDNw/qEjr3iMf4wlgOcxfkEH53DblTAWG2VN65ynlHI0vTYTcMX4B18C+BDQTEXaYmA4xiO2PcD6hxJxZ6AZUBvohbHX7i2gNOAIDMskppcw/oGVwdjT+AYYT3MA04AeQHnTv5JZ2UlT0u0NHH1M0zrAMxh7FlNompYI/IyxiEvWA7gJ/ATs4r9EmbzNM0qpLplsy8L0w3NJKfWdUqpEFvelPFAZCEn1XF/ghqZpO7OyjhxwB44lP9A0LQk4Qapkq5RqqZQKB8KBDhi/PzkxEuN3pw5QH+iWeXMhdCF5Oy1zyNv7gXeV8bT8I4WiMg53mpaFuAzAy8APj2n3ilIqAuPvRE1gQaq4ZgH9TQc87hi/Q7mZv12BCE3T/k313DHS5mwbZRy7HAt8BXyekw09yWdsrqR4NeoDjNc07ZamadeB/wP6ZtRY07SDmqYdMh2x/YPxD6JZDrf9APhc07R4TdPWAhpQXSllhbEYe13TtGumbf1u+qPrDazVNG23pmlxwCSMX1SfVOudbtqfi8A+4A9N0/7SNC0aY0J8uNcwte80TftH07QojIWhl+n5HhiPmvebev8+5PHfoQ9Mf5xnMJ5CSv1j0dfUA5H8r4RpPzTgejrrukraP8j+wHLTkeuPwMsq1ZgmTdNqmXpJ03MV46mrikADjD8Oix6zLyilbEzb+kbTtAum54phPB02JpNFP3hoX7993LYe4oixKE0tHGPvKwCapu00DRuoiDGRXczmNpL1AKZpmnbV9KOf0XAKIfQkeTstc8jbnwLTMebuI6Ye897JDTVNG6Rp2luPiQ2gJcbc90tmjTRNW6RpmhPGoQrfYSxik/2C8bOKxtgRMFvTtNRDGJo9tJ8nsxBXalnJ2bGacexyMeBNUnVQZFNOPmOzVqh3LiuUUgrjEWPqo6N/gXKZLFNbKRWojBfG3Ac+JudHOTdNvWjJHmD80pcBLIHz6SxTNnW8psQY9lDMqZNIdDqPHcnYtXTiSd7upVTbvc+jf5wPm6hpWjFN08pomtb5oaPQpabXkv/dAW5h7NEonc66ypheRylVDWgMJJ9yWoNxLGi6F1I9TNO0cE3Tjpp+XK4Ao4COSinbjJYx9ZCswDimNXWCnQjM1zTtciabnPjQvg4xrdMmdYLE+EO2PdVzo03LR2Icx5qaM8bxaA/v2yWMPdE5HT+V5nN+6P+F0J3k7XQV+LxtKvZnaJrWEGO+ngYsMeXz7OgPrNQ0LSYrjTVNO43xuoUZAEqp0sBmYBzGaxoqA12VUgNTLbbnof1M3WO6PVXOHg2MT5WzfzI1y07OjgDmAStTD/fIhpx8xmatyBevpl67axhPAyWriDGpgPFo8mHfAUeAappxMPgEjH+4uekqxrGTVdN57Qqp4jX1Npbjv5jzylWMpyOSt+sMuOTyNk5gTNjdUz9p2sfO/HdaJ/kCre1KqWvAWYw/GmmGDmSDhvEzTPdzNJ2mWgLYA71MPzzJWgLvmMYYXQOeBtalKjwz3qjpyDv5HxAEtE713HRT05MYLzpIHU8d0/PpsQQy+0EwmHqJkqW+cjnN54xxbJgQBYbk7WzHVFDydgpN0x5omjYN4ylz16xuSCnlZFpnpkMG0pE6J9bAeEp/hakD41+MQx46ZGVFmqa1TpWzpwOfpsrZycOsTgPOSqmKqRatS8Y524CxV/aZTDadunMlw5ydR59xgVLki1eT5RiPnJ5SSpUCPuC/XqvrQKlU45nA+AUL1zQt0jRW5tXcDkjTtHiMxdIMpVRppZSFUqqJKRmsBDorpZqaCpD3MPYGBuV2HA9ZhfHo1FcpZY0x+Sc9Zpls0TQtAXgX+Ewp1d3UM1kWY6KyBGabel36Au9jPDWW/O8l4EX16IwEj1BKNVRKVVdGpTAmoG2m03MPt1XA9xiTQ2fTaZnUmmAcd5ccx23gFdJeofuktgN2prFiNhhPMUUBe00x9jONxUUpVQXjZ5PZ+C0FDFRKGZRS/kAtwMnUu7wKeFMp9YwyXvzydi7uhxC5RfJ21hSIvA2glBqjlHpeGafMslJKvYZxWEJ2Tpd3B8I0TduXWSOl1KvKOFWVUkp5AGP5LyeeAhyUUt1Mr5fDOLY/p6ftH6Fp2l1gIzBBGadb9APaYTpbqJRqr5TyNH1HXDD+Bl0BzmWy2n5KKWulVHWMFys7mb5Lef4ZFzRSvBp9jPHim5NAMPAH/43zOwasB/5V/43veRMYrJSKBOZgTEp5YRTGUx1HMSa5zwBlGpczCPgW4xielsALpgSSZzRNOwq8g/FK3jCMR3vhGI+cc3M7P2Dcv/cwTjPyl+mlJpqm3cM4Tq0UxnGn15L/YRznFYZx/A/KeKV81ww2UxPj1CSRGD/ze6TqtVVKfaqUWpuq7QCMF2DcUP9dwdzVFO+th+JIAu6Yxp4l+0ilnS8wsyEG6b0n0fw3s8I9jBdzvJjqM/cEDiqlooDfMX5vh2eyykTTMjeBERjHC04EqmD8odmH8W/iEMYEnKufsRC5QPJ2FhSgvI1pmzMxHlzcwHiQ/2LykCul1GKl1HQy1590el2VUq2UUrdSPdUc4/cjEuN44TXAJ6ZYb2MsgsdhzKeHMea6r1It76cenefV47FvRFqvYhyacgvjxXqDNE1LLk5LYPzNuo+xYC0LtHvM9yEW49CTdRi/V62BHvn1GRckStPSO7sixOOZxubcAcpqmnZV73hE1iilXIETmqZlaZogpVRnYLKmabXyNjIhRF6TvG2eTEPSummatjcLbQv9Zyw9ryJblFIByji/oiPGwfYHCusfR1GllHJSSrUxnc6qiPHK1YwmAhdCFHCStwu/ovYZS/Eqsqs7xgslLmO82CDT2zQKs2QAJmM87XQI40Uu/6drREKIJyF5u/ArUp+xDBsQQgghhBBmQ3pehRBCCCGE2chW8aqU2pJXgQghRF4qivmrKO6zEKJwyCx/Zelq42TOzs5tfXx8ZJyBEMIc3dc7gPwmOVsIYcYyzNnZKl5r1KhBUFBez6cshBC5TymV2eTfhZLkbCGEucosZ8uYVyGEEEIIYTakeBVCCCGEEGZDilchhBBCCGE2pHgVQgghhBBmQ4pXIYQQQghhNqR4FUIIIYQQZkOKVyGEEEIIYTakeBVCCCGEEGZDilchhBBCCGE2pHgVQgghhBBmQ4pXIYQQQghhNqR4FUIIIYQQZkOKVyGEEEIIYTakeBVCCCGEEGZDilchhBBCCGE2pHgVQgghhBBmQ4rXQsbCwgIvLy/q1KlDp06duHfvXq5vIzY2llatWuHl5cXKlSv5/fffcXd3x8vLi+jo6FzfnhBCFEaSr4XIGSleCxk7OzuCg4M5ceIEJUqUYM6cObm+jaNHjxIfH09wcDA9e/Zk2bJlvP322wQHB2NnZ5fr2xNCiMJI8rUQOSPFayHWsGFDwsLCANi9ezf+/v4pr40cOZLFixcDULlyZcaPH0+9evXw8PDg9OnTANy5c4cXX3wRT09PnnvuOY4fP86NGzd4+eWXCQ4OxsvLi2+//ZZVq1YxYcIE+vTpk+/7KIQQhYHkayGyTorXQioxMZGdO3cSEBCQpfYlS5bkyJEjDBs2jClTpgAwfvx4vL29OX78OJMmTaJfv36UKlWKBQsW8PzzzxMcHMyQIUMICAjgq6++YtmyZXm5S0IIUShJvhYie6R4LWSio6Px8vLiqaee4s6dO7Ru3TpLy3Xp0gWAZ599ltDQUAD27t1L3759AWjRogW3b98mPDw8T+IWQoiiRvK1EDkjxWshkzyG6t9//yUuLi5lDJWlpSVJSUkp7WJiYtIsZ2NjAxgvIEhISABA07RH1q+UyqvQhRCiSJF8LUTOSPFaSLm4uDBz5kymTJlCfHw8lSpVIiQkhNjYWMLDw9m5c+dj19G0adOUU0u7d++mZMmSODs753XoQghRpEi+FiJ7LPUOQOQdb29v6taty4oVK+jbty89evTA09OTGjVq4O3t/djlP/nkE1555RU8PT2xt7fnhx9+yIeohRCi6JF8LUTWqfRONWTEx8dHCwoKysNwhBAibyilDmua5qN3HPlJcrYQwlxllrNl2IAQQgghhDAbUrwKIYQQQgizIcWrEEIIUYBs2bKFWrVqUb16dSZPnvzI65qmMWrUKKpXr46npydHjhwB4NKlSzRv3hw3Nzfc3d2ZMWNGfocuRL6Q4rWIyGkyPHPmDF5eXin/nJ2dmT59en6HL4QQRUJiYiIjRowgMDCQkJAQli9fTkhISJo2gYGBnDt3jnPnzjF//nyGDRsGGKfYmjp1KqdOnWL//v3MmTPnkWWFKAykeC0CniQZ1qpVi+DgYIKDgzl8+DD29vZ07txZj90QQohC7+DBg1SvXp2qVatibW1Nr169WLduXZo269ato1+/fiileO6557h37x5Xr16lTJky1KtXDwAnJyfc3NxSbjkrRGEixWsR8CTJMLWdO3dSrVo1KlWqlJ/hCyFEkREWFkaFChVSHpcvX/6RAjQrbUJDQzl69CgNGjTI24CF0IEUr0VAbiXDFStW0Lt377wNVgghirCs3CnrcW0iIyPp2rUr06dPlxsViEJJitciIDeSYVxcHOvXr6d79+65H6AQQgjA2HFw6dKllMeXL1+mbNmyWW4THx9P165d6dOnD126dMmfoIXIZ1K8FgFPmgzBOCa2Xr16lC5dOu8DFkKIIsrX15dz585x4cIF4uLiWLFiBQEBAWnaBAQEsGTJEjRNY//+/bi4uFCmTBk0TWPQoEG4ubnx1ltv6bQHQuQ9KV6LgCdJhsmWL18uQwaEECKPWVpaMnv2bNq2bYubmxs9evTA3d2defPmMW/ePAA6dOhA1apVqV69Oq+++ipz584F4I8//mDp0qX8+uuvKTPEbN68Wc/dESJPyO1hi4jNmzczevRoEhMTGThwIB988EFKIhw6dCiapjFy5Ei2bNmCvb09ixYtwsfHeFe2Bw8eUKFCBc6fP4+Li4ueuyFEjsntYYUQwnxklrOleBVCFAlSvAohhPnILGfLsAEhhBBCCGE2pHgVQgghhBBmQ4pXIYQQQghhNqR4FUIIIYQQZkOKVyGEEEIIYTakeBVCCCGEEGZDilchhBBCCGE2pHgVQgghhBBmQ4pXAUBiYiKfTfycdevW6R2KEEKITOzYsYOnnynL+E8nkJ0bDQlRWEjxKrh//z4t23Zg6sIV9Bn4Gp/930RJiEIIUQDN/eYbXuzeC57rz8yFK+jWszcxMTF6hyVEvpLitYi7cOECXj4N+CvcCpcuEyjW60umzl9Crz59iY2N1Ts8IYQQJm+PfY8RI0Zg4dYS++r1sfR+gfWbAmnQuKkUsKJIkeK1CPvjjz+o59uAiMrNcGwxFGVhiaVTSZy7TWLH8Ys0fN6Pmzdv6h2mEEIIoHKVKvTq0xfHyweI/GsHdzdNpVuXF2nQoD42NjZ6hydEvpHitYj65Zd1NPNrTtxTNbCp/hxKqZTXDNa2WHn5c/KfS3h4P8ulS5d0jFQIIQTAyGFDWL70B571rQ+ahrWtLcuX/sD8ubPT5HAhCjspXoso52LF6PnSy7gXTyImaM0jr0dt+ZoOzRvRqmUrnJ2ddYhQCCGEEOJRlnoHIPTRwq8ZLfya8c033/Dxos2PvK6AGdOmULFixfwPTgghhBAiA9LzKoQQQgghzIYUr0IIIYQQwmxI8SqEEEIIIcyGFK9CCCGEEMJsSPEqhBBCCCHMhhSvQgghhBDCbEjxKoQQQgghzIYUr0IIIYQQwmxI8SqEEEIIIcyGFK9CCCGEEMJsSPEqhBBCCCHMhhSvQgghhBDCbEjxKoQQQgghzIYUr0IIIYQQwmxI8SqEEEIIIcyGFK9CCCGEGdE0DdD0DkMI3UjxKoQQQpiRqKgoAGxsbHWORAh9SPFaxBmP4NN9JV/jEEII8XiapnHhwgVQBpxciukdjhC6kOK1iIuMjCTJwvqR5xMT4rGxsdEhIiGEEBk5dOgQt+9FQOQt2rdto3c4QujCUu8AhL6OnQhBsy+R5jktIY74mAe4uLjoFJUQQoj0TJ0+C0u3FsQe28yb3+/VOxwhdCE9r0VYeHg4a9euxcGtWZrn426GUqFSVWxtZTyVEEIUFLt27WLjpk1gbU+dOu64urrqHZIQupCe1yJs8Q8/YF+1HhaOxdM8H3vlDH71fXSKSgghxMN27dpFQJdu2D3/CrF/LuPTH3/QOyQhdCM9r0XU6dOn+fiTCVh5BaR5XtM0OLOL/n166xSZEEKIZIcOHaLXy/3wf7ErNo0HEH9kLePeeZO2bdvqHVqRduXKFeLj4/UOo8iSntci6NixY7Rs2x6rhn2xKVMzzWtxV89io8XSpo1cCCCEyDsXL14kNDQ0z9bfrFmzxzfKQ+06BmBjZ5ejZTVNIyrqARcunOfmnXtYubXEtkFPYv9cxrh33mTcu2NzOVqRXa6169C9e3cWzJ+HUkrvcIqcQlO8njt3rsAfBSUkJBAXF0dsbGya/8bFxWW6nLe3N9WqVXuibWuaxp9//snUGbMJ3LwJhxbDcHB7Pm2bpERi9i5mwrvvYDBIp7wQIu+MHP0WW7ZsoVi5GnqHkif2RxbHQpXM2cJJSSQ+uA1lSmJwvEXEgdXUqOnKJ/Nm0bZtWyIiInI3WJFtiYmJrFi7gTru7rw5epTe4RQ5KuN5Ph/l4+OjBQUF5WE4OXPz5k1Kly6NffFSWNva5/r67165kOvrzK469XxzvGxcbCx3bt0iDgss3dtgX6cVFnZOj7SLDPqFypEnObhvrxSvotBRSh3WNK1IDeYuqDkbIDQ0FG+f+lg3H4FdtSL1sTxW/J0wrnw3JM1zNnYOOkUj0mPlVBynTu9zf83HrP5xCe3bt9c7pEIns5xdKHpen376ad4e+x7frdqIU5cJKEurXF2/Q/R9Io9vy7SNMliiLCzBwgplYWl8bGllfJz8/yrnBeHNHC8JymCBtZ0T9sWeQWUQQ3RoMHGHf2bFoQNSuAoh8lzlypXZuG4tbTsGYOH8GdZPV9I7pALDqkQ5Kr27Ue8wRBY4dnyHni+9zP4/fqd27dp6h1NkFIqeV4CkpCQ6de7K/otROLYZJWNQsiH6/GGitk1n07q1uo8TEyKvSM9rwbRkyVJGvP0exXp9hYW9zC0tzE/UiZ1YHf+Z40eCKFkyh0NFxCMyy9mFpovNYDCw6sf/USrhBlGH1ugdjllIio8hcvf3xO+aw9ZNG6RwLWD27z/ApUuX9A5DiDzVr19fhr7Sj8iNk9ESCvZ1C0Kkx6FOS+IqNKCd/wvExsbqHU6RUGiKVwAHBwd2bNkEIVt5cHaf3uEUWEmxUUQc3sC9paNoWtGac6dDaNy4sd5hiYf0G/QqLdq0JzIyUu9QhMhTX3w+kYbuVYncOYfsnA0UoqBwaPIyFyINDBj0qnyH80GhGPOaWrly5Yy9iC1bY+lSGuvST3aVfmEQcXQz8XfCiL8ZSsL1f0iIieTp0qV5c9RImjZtysmTJ/UOUaQjKjKScM2OLt17EbhxHRYWFnqHJESeMBgMrF6xDJ/nGnPz0M841u+qd0hCZItSBhzajCZw9TgmTf6CD8a9p3dIhVqhGfP6sFWrVjNo+Ou49PwSS6en9A5HV/9+4Z/msVPpitg6ldApGpFVSZbW2LV8ncgtU+gf0JLp06bqHZJZkzGvBV9YWBh16/lgaDwI+5oN9Q5HiGxLuH+LeyvHsnTBPDp37qx3OGat0M82kJ4ePbpz8tQpZiz8HOdu/4fBylbvkHQjV62aN8eO77Jw2VjquNdm8KBBeocjRJ4pV64cWzZtwK9VGyxdSsmZM2F2LJ1L4uT/Hn1fGcTeKlXw8vLSO6RCqVCNeX3YJx9/RKvnvIjaPgtNS9I7HCFyxMLOGaeADxg9Zix79uzROxwh8pSPjw8L539LxIZJJETe0TscIbLNpkxN7Jq9RpsO/ly7dk3vcAqlQl28KqX43w+LqGAdTdSfK/QOR4gcs3qqAg5t3+SFLt34559/9A5HiDzVo0d33np9OJEbPicpXq7eFubHwe15kmo0p3V7f6Kjo/UOp9Ap1MUrgK2tLVs3rcfin9+ICpFeK2G+7Kp4Y+nTnVZtOxAeHq53OELkqU8+/ogWDTyJ2j5Trt4WZsnhuZ5cSXSid9/+8h3OZYW+eAUoXbo02wM3EfPbAmLDTusdjhA55ujdkYgSrvi/2JWEhAS9wxEizyil+HHJYspbPSDqz+V6hyNEtimlcGj9OnuCTvLR+E+ytayFhQVeXl7UqVOHTp06ce/evbwJ0iQoKIhRo0Y9tl2jRo3yNI6sKhLFK4CnpyfLliwmYtMXJNy/oXc4QuSYg98gTl65z/DXH59ohDBntra2bNu8Qc6cCbNlsLLB0X8cM+bOZ/nyrA9ftLOzIzg4mBMnTlCiRAnmzJmTh1Eax5rPnDnzse327SsYc+gXmeIVICAggI/GvUPkhkkkxckYFGGelMECxw5vs3JdILNm521CE0JvKWfOfv+e2Ctn9A5HiGyzcCyOU8D7vDpsOAcPHsz28g0bNiQsLCzl8VdffYWvry+enp6MHz8egNDQUFxdXRk8eDB16tShT58+7Nixg8aNG1OjRo2U7R48eJBGjRrh7e1No0aNOHPG+De1e/du/P2N02p+8sknDBw4ED8/P6pWrZqmqHV0dExp7+fnR7du3XB1daVPnz4pQyM2b96Mq6srTZo0YdSoUSnrzU1FqngFGPv22/i3aELUlq/RkhL1DkeIHDHYOODY6QPGfTSebdu26R2OEHnK09OTZT8sImLjZDlzJsySdamq2LccSftOL2Trtt+JiYns3LmTgIAAALZt28a5c+c4ePAgwcHBHD58mN9++w2Av//+mzfeeIPjx49z+vRpfvzxR/bu3cuUKVOYNGkSAK6urvz2228cPXqUCRMm8P7776e73dOnT7N161YOHjzIp59+Snz8o7duPnr0KNOnTyckJITz58/zxx9/EBMTw5AhQwgMDGTv3r3cvHkzu29VlhS54lUpxcLvvqVaMQNRe5fqHY4QOWZVvAyOHd6hW6+XOH1axnKLwi0gIIAP3n1bzpwJs2VfowGGOh1o1a7jY2/7HR0djZeXF0899RR37tyhdevWgLF43bZtG97e3tSrV4/Tp09z7tw5AKpUqYKHhwcGgwF3d3datmyJUgoPDw9CQ0MBCA8Pp3v37tSpU4c333wzwztsduzYERsbG0qWLEmpUqW4fv36I23q169P+fLlMRgMeHl5ERoayunTp6latSpVqlQBoHfv3jl9uzJV5IpXAGtrawI3/IJt2GGi/tqudzhC5JhthTrYNHyZVm07cPv2bb3DESJPvTf2HTo2b0zU1q9l7m5hlux9OnPbpizder5EUlLG3+HkMa///vsvcXFxKWNeNU1j3LhxBAcHExwczAqXS9MAACAASURBVN9//80g081rbGxsUpY3GAwpjw0GQ8oFvh999BHNmzfnxIkTbNiwgZiYmHS3n3pdFhYW6V4gnF6b/JpVoUgWrwBPPfUUO7ZuJnbf/4i5dELvcITIMQeP1sSUe5b2nV4kLi5O73CEyDNKKRYtmE81FwNRvy/ROxwhsk0phUPLoRw4/S9vj333se1dXFyYOXMmU6ZMIT4+nrZt27Jw4cKUntuwsDBu3Mj6UJrw8HDKlSsHwOLFi3O0D5lxdXXl/PnzKT29K1euzPVtQBEuXsH4Jv+08kciN39F/N2reocjRI45NOnLP/eSGPjqEJlPUBRqyWfObMKCiPprh97hCJFtysIKx47vsmDJchYuXPTY9t7e3tStW5cVK1bQpk0bXnrpJRo2bIiHhwfdunUjIiIiy9seO3Ys48aNo3HjxiQm5v51P3Z2dsydO5d27drRpEkTSpcujYuLS65vR2Xnh87Hx0cLCgrK9SD0Nmv2HD6cNAWXnl9gsHHQOxwhciQpLpr7q8fx/qghjH3nbb3DKXCUUoc1TfPRO478VFhzNhgvKKnfqAn27d/BtkIdvcMRItvibl3k/k8fsmXjOpo0aaJ3OLkmMjISR0dHNE1jxIgR1KhRgzfffDPb68ksZxfpntdkr48cQc8X2hO5+SuZgUCYLYO1HY6d3mfC51+yYcMGvcMRIk+5urry0wrTmbN7cv94YX6sS1bEoc0bdHqxC+fPn9c7nFzz3Xff4eXlhbu7O+Hh4QwZMiTXtyE9ryYJCQm0aNOekEh7HFu8pnc4QuRYbNhpIjdO4o89u/D09NQ7nAJDel4LpxkzZ/Px5Kly5kyYragjG3G8sIvgoAN5cordXEnPaxZYWlqyYe1PON05TeTRTXqHI0SO2ZRzxfb5QbRu3zHd6U2EKEzeGDWSHgHtiNw8Rc6cCbNk792RiGI1eKFLd7ntdxZJ8ZqKi4sLO7ZuJiFoNdEXjuodjhA55lC7GYnVmtK2Y0C6U6Ek3yUlL1SuXJlbt25luX1wcDCbN2/Os3hE4ffN7Fm4l3Umas/3eociRLYppXDwG8zxS3cYOSr7Y0Ozo1GjRnm6/vwixetDqlWrxvq1a4ja+jXxt7N+FwwhChqHhr24FGfHy/1fKdAzEEjxKp5U8pkzx1uniAqW75IwP8rCEscO77B87XrmzJ2bZ9vZt29fnq07P0nxmo6mTZsyfeqXRKyfSGL0fb3DESJHlDLg0Pp1duwP5tPP/u+x7W/evEnXrl3x9fXF19eXP/74A8j4XtiJiYm8/fbbeHh44OnpyaxZs1LWNWvWLOrVq4eHh0fK3b/SW09cXBwff/wxK1euxMvLi5UrVxIVFcXAgQPx9fXF29ubdevW5cG7IwqbYsWKsXNbIPGHVhEdGqx3OEJkm8HWEceAD3n3g4/Zvj1vbqCUfNZt9+7d+Pn50a1bN1xdXenTp0+B7uR4mBSvGRg8aBCv9OlB5KYv0BIfvaevEObAYGWLY6f3mTpzDqtWrc607RtvvMGbb77JoUOHWLNmDYMHDwYyvhf2/PnzuXDhAkePHuX48eP06dMnZV0lS5bkyJEjDBs2jClTpmS4HmtrayZMmEDPnj0JDg6mZ8+eTJw4kRYtWnDo0CF27drFO++8Q1RUVB69Q6IwqVatGut+/omoLdOIv31Z73CEyDar4mVxbD+Gbj17p3QU5JWjR48yffp0QkJCOH/+fEqHhTmw1DuAgmzaV19y6tQLHPn1WxxajUAppXdIQmSbpWMJnDq9z8DXhlC1ahV8fNK/4H7Hjh2EhISkPL5//z4RERGEh4fTv39/zp07h1KK+Pj4lPZDhw7F0tKYRkqUKJGybJcuXQB49tln+fnnnwEyXM/Dtm3bxvr161OK3piYGC5evIibm9sTvhOiKGjWrBlfT/mCMR98ikuvL7Gwc9Y7JCGyxbaiJ4nP9aFV2w4cO3IoTW7NTfXr16d8+fIAeHl5ERoaajbzzUrPayYsLCz4efUKSkRf5sFhOXUpzJd16WrYtxhGu46dCAsLS7dNUlISf/75Z8o9s8PCwnBycsrwXtiapmV4QJd8z+vU98TO6j21NU1jzZo1KXFI4Sqy69XBgxnwUg+iNn0pZ86EWXLwbEN0GS86BnTO8ED/SSXnaUibq82BFK+P4ejoyI4tm0g8voEHfx/UOxwhcsy+ZiNwb0erdh3Tfb1NmzbMnj075XFwsHHcYEb3wm7Tpg3z5s1LSXh37tzJdPsZrcfJySnN7Q3btm3LrFmzUsZfHT0qM3/oKS9npsjMxx9/zI4dOb/969dTvsSnRjmids3n8jcDufL9CK4sHMn1lR+RGHkX4L/nF73OlUWvE3P5FABxN//l2vL3CZv/GmHzX+XeH8tTvo+JUXe58dOnXFk4kisLhnF99fgn31kh0uHwfH/O3Ilj0GtDzWo8an7I12EDjo6OREZG5ucmAWMSbNq0Ka1atcrR8hUrVmTz+l9o1a4jd7aCwdYJlMLCoTglO76FhWNxLn8zEIO1HRiMxwMlWg/HtrwbcTf/5c6Ob0mMuAVoOLi3wKVRrzQ9VjfWfEbig3uU6Ts1N3ZXiAw5+HblxrYwoqKOpZwuAnjrrbeYOXMmI0aMwNPTk4SEBJo2bcq8efMYO3Ys/fv3Z9q0abRo0SJlmcGDB3P27Fk8PT2xsrLi1VdfZeTIkRluO6P1NG/enMmTJ+Pl5cW4ceP46KOPGD16NJ6enmiaRuXKldm4cWPevCGiwJowYcITLZ985szb9znuxj2gdP+vsbB34e6eHwjfv4oSrYx3/SndexIW9v9NDJ8UH8vNnz+jRJvh2FWpR1J8DDfXfk7k0U041fPn3u/LsK3shbPPCwDE3bjwRHEKkRFlsMCx7VusXz2OKVOn8c7bY/QOqcDI1zts6VW85pZly36kb79+PPPKLKxLVuTunh/Q4mMo0WoIl78ZSBlTckyWFB/L1YUjHkmC9tV9carnb2wTE8mVhSNRVraU6v4JVsWe0Wv3RBGhJcRz/+ePebWHP19OnqR3OPlG7rCVM8l5e/fu3XzyySeULFmSEydO8Oyzz/K///0PpRSHDh3ijTfeICoqChsbG3bu3ImVlRXDhg0jKCgIS0tLpk2bRvPmzVm8eDG//PILiYmJnDhxgjFjxhAXF8fSpUuxsbFh8+bNlChRggEDBuDv70+3bt2oXLky/fv3Z8OGDcTHx7N69WpcXV25efMmL730Erdv38bX15ctW7Zw+PBhSpYsmRL/xYsXqVKlCsU7vIWjux/R5w9z//B6Snf/NN28HXFsG7GX/qKk/3+FQvzdq1xfPo7ywxdzY80EHOq0xKFW4yd6X4XIqoTwG9xb+S4rlizE399f73DyTYG7w1ZmUzQcOnSIRo0aUbduXerXr09ERAQxMTG88soreHh44O3tza5duwDjqccXX3yRTp06UaVKFWbPns20adPw9vbmueeeSzmNOWDAAH766SfAOIH6+PHjH5nG5+bNm7Ru3Zp69eoxZMgQKlWq9MhE6336vISTkyNR22agJcRhW6EO8XevZLifUSF7sCnnhl2VeoDxyu8SrYcSvv+n/9qc2Ydd9fo4uDXlwanfcukdFiJjytIKx47vMmPmTI4cOaJ3OMKMpHd1clxcHD179mTGjBkcO3aMHTt2YGdnx5w5cwD466+/WL58Of37908Z53zixAl+/PFHDh48yAcffIC9vT1Hjx6lYcOGLFmyJN1tpzeDxaeffkqLFi04cuQInTt35uLFi48sV7FiRUqVKkX0ngXE3Qwl+p+DWD9dOeX168vf58qi17m65C0A4m/9i/Uz1dOsw6p4GbS4aJJiH+Dk3ZHbgTO5tnwc4ftWkhBx+4nfVyEyY+lSCqeOY+n9cr90v+NFkW5jXs0xCQIUL1aM5zxqELljDg/+frIkCPDg1B4c3JriULsZUVK8inwSfXInVapWo1atWnqHIsxI8tXJBoMh5erkM2fOUKZMGXx9fQFwdnbG0tKSvXv30rdvX8A4TVqlSpU4e/YsYBwq4uTkxNNPP42LiwudOnUCwMPDg9DQ0HS3nXoGi+Q2e/fupVevXgC0a9eO4sWLp7usjY0NJVwcuf7DaBKjwnF5rnvKa6V7T6LsK7Mo029aqiUynlnGruqzlBuyAKe6bYm/c5mri98g8UF45m+cEE8oKSYSS0tLrK2t9Q6lQNCteDXXJAhwOfQ8UWf2EXfl1BMlwcSou8TfvYpNeXesSpRDGSyIu5l+zELklgdn/4STgezYsgkHB4csLbNlyxZq1apF9erVmTx5crptdu/ejZeXF+7u7jRr1gyAM2fO4OXllfLP2dmZ6dOn59q+iPyV3tXJGc06kdmQtNTrMRgMKY8NBkOGVzynN4NFdoa9HT16lHfGvot17F2UZcYFgFXJisRdO5fmufh711DWdhhs7I0x2DnhUNuPkv5jsClTg5hLJ7IchxDZFXczlMhtM9m0/heeeSZ7QwsnTpyIu7s7np6eeHl5ceDAgQzbDh48OGW6xOze5ju/6Va8mnMS3LNnDxfO/4Ozisn0Ti6PS4JRp34nKSaSsHmDuPzNQBLCrxN16vcsxyFEdsVdP8+DX+cSuHF9mgu2MpOYmMiIESMIDAwkJCSE5cuXp5kPFuDevXsMHz6c9evXc/LkSVavNt4QoVatWilTXh0+fBh7e3s6d+6c6/sl9OPq6sqVK1c4dOgQABERESkX/C1btgyAs2fPcvHixVzv6W/SpAmrVq0CjPMD3717N9P2k/5vAk28ahG5Y06GOd+hth+xl0NScntSfCx3d3yLc4OuAET/e4ykeOOZv6TYB8TfvYql89O5tUtCpJEYdY+I9ZOYN3smjRo1ytayf/75Jxs3buTIkSMcP36cHTt2UKFChQzbL1iwgNq1az9pyPmiQE2VZU5JsEyZMmzbvJHo3d8Se/Vcum0elwSjTu2hdI9PKT9sIeWHLeSZ/tNl3KvIMwmRd4jYMJEF875JObuRFQcPHqR69epUrVoVa2trevXq9cgtW3/88Ue6dOlCxYoVAShVqtQj69m5cyfVqlWjUqVKT7YjokCxtrZm5cqVvP7669StW5fWrVsTExPD8OHDSUxMxMPDg549e7J48eI0nQ25Yfz48Wzbto169eoRGBhImTJlcHJyyrC9wWBg5bKlPKPdIergT+m3sbLh6a4fEb5vBWHfDeHqwpFYP1Mj5SLbuGt/c+2HN7mycCTX/vc2TnXbYFOmZq7ulxAAWkIckZsmM2xQf/r2fTnby1+9epWSJUum/N2VLFmSsmXLsnPnTry9vfHw8GDgwIHExsYC4Ofnx5Ne4JlfCtQdtlInwejoaOzs7NixYwfDhw9n6NCheHh4YGlpmWdJsHfv3qxcuZJmzZo9NgkCeHt7s2ThAvq/Ohy0xEdeT06Cd7bP4872byApCQf35jjV8ych/DqJ929iXdY1pb1VsWcw2NgTe+UMNmVlLKLIPUnxsURt/JzRI4bSq1fPbC0bFhaW5mi9fPnyj5x6Onv2LPHx8fj5+REREcEbb7xBv3790rRZsWIFvXv3zvlOCF0kzxDj5+eHn59fyvOp5wT29fVl//79jyybej7fZAMGDGDAgAEpj1MP70r9WuplU7fx8fFh9+7dALi4uLB161YsLS35888/2bVrV7q/DamXt7e3Z3vgRuo+60vxFoPSzDSQzPrpyjzzUvrDY1wadMXF1AEhRF7RNI3IHXNp7FmDyZP+L0fraNOmDRMmTKBmzZq0atWKnj170qBBAwYMGMDOnTupWbMm/fr145tvvmH06NG5vAd5K1+nyirIYmNjsbCwSEmCw4YNS5mk/XE+mziJqd8uwbn7RAxWtnkcqRDZo2kakYFTaVarFGtWrcj2bY5Xr17N1q1bWbBgAQBLly7l4MGDzJo1K6XNyJEjCQoKYufOnURHR9OwYUM2bdpEzZrGHqm4uDjKli3LyZMnKV26dO7tXDbIVFmFz7lz5+jRowdJSUlYW1szd+7cLJ9VOHLkCM/7taB4n2lYuujznRQiI5EHVlPq9jGC9v+Bvb19jteTmJjI77//zq5du/j2228ZN24ca9as4bffjGd5d+7cyZw5c/j555/x8/NjypQp+Pj4ULlyZYKCgtJMO5ffMsvZBarnVU8XL15MkwS/++67LC/74fvjCDp8mD/3/Q/7ZoPzMEohsi9q/wrKWUaybMnmbBeuYOxpvXTpUsrjy5cvU7Zs2UfalCxZEgcHBxwcHGjatCnHjh1LKV4DAwOpV6+eboWrKJxq1KiR4zuwnT17DoO1HUo6HEQB8+DMPji1je2HDz1R4QrGa3uSz5p4eHjwww8/5FKU+ipQY171lJwEjx07xqFDh7I1JvDGjRvsP3AQrZSMexIFS9Sp3zCc2822zRuws7PL0Tp8fX05d+4cFy5cIC4ujhUrVhAQEJCmzQsvvMDvv/9OQkICDx484MCBA7i5uaW8vnz5chkyIJ7IwIEDKVWqFHXq1HnidR06dIjBQ4fh1On9dIcNCKGX2Gt/82DXN2zdtCHldto5debMGc6d+++anODgYEqXLk1oaCh///03YDyTljw7jDkpksVrbibBmJgY2nYMILFaUxzcmuZCdELkjtgrZ4jZ8x3bAzdle3qV1CwtLZk9ezZt27bFzc2NHj164O7uzrx585g3bx4Abm5utGvXDk9PT+rXr8/gwYNT/r4ePHjA9u3bU6aoEyInBgwYwJYtW554PZcvX6adfwD2LUdgXbpqLkQmRO5IiLhNxIbPWTj/W5599tknXl9kZCT9+/endu3aeHp6EhISwuTJk1m0aBHdu3fHw8MDg8HA0KFDcyH6/FUkx7z+9ttvODo60q9fP06cyPn8fJqm0a1nb3afvoZj+zEoVSSPBUQBlHD/JvdWjGXZou944YUX9A6nQJAxr+YvNDQUf3//HOftqKgonm3QiFulnsWxvlx0JQqOpPgY7q/+gDcH92H8xx/pHU6BUOBuD6u3pk2bUqJEiSdezycTPuPXA8dxaD1KCldRYCTFRRO5YSLvjx0jhasQJklJSXTv1Ydb1qVx8JWzAKLg0LQkorbNpHWjenz80Yd6h2MWpOLKoVWrVjNt1lwcOo3DYJW703YJkVOalkTU1q/p4NeIce+O1TscIQqMse+9z58hF3BoMTxHFy4KkVei9i2nom0sSxcvlO9mFslsAzkQFBTEwNeG4Nz5Eywdn7wHV4jcErV3KVWdFYu//06SoBAmixf/wPzF/8Ol15coSyu9wxEiRVTIbizO/87Wo4extZWZL7JKel6zKSwsjHYdO2HfYhjWpavpHY4QKSL/2oH1pYMEbvgFa+uM790uRFGyd+9eRo5+C8eAD2RmAVGgxIadIua379mxZXO6dyV8Elu2bKFWrVpUr16dyZMfveFGeHg4nTp1om7duri7u7No0SIALl26RPPmzXFzc8Pd3Z0ZM2bkaly5RYrXbIiKiqJVu47g3g77mtm7x7AQeSnm0gni9i1hx5bNuk4qLURe6d27Nw0bNuTMmTOUL1+e77///rHLXLhwgU4vdsGhzRtYl6yYD1EKkTUJ4TeI2PQFPy79AQ8Pj1xdd2JiIiNGjCAwMJCQkBCWL19OSEhImjZz5syhdu3aHDt2jN27dzNmzBji4uKwtLRk6tSpnDp1iv379zNnzpxHli0IimTxmpMkmDzY/4bl0zj4ylWqouCIv3eNyM1fsWr5MmrXrq13OELkieXLl3P16lXi4+O5fPkygwYNyrT9/fv3adm2AxbeXbCr+uTTDgmRW5JiHxCxfiIfj3uXTp065fr6Dx48SPXq1alatSrW1tb06tWLdevWpWmjlCIiIsJ4B8bISEqUKIGlpSVlypShXr16ADg5OeHm5kZYWFiux/ikiuSY1+XLl2d7mffe/5A/T57HucsEGUsoCoyk2Cgi10/k/z4ZT7t27fQOR4gCITExkRe79uC+SzUc6vnrHY4QKbSkRCK3fk1A66a88/aYPNlGWFgYFSpUSHlcvnx5Dhw4kKbNyJEjCQgIoGzZskRERLBy5UoMhrT9maGhoRw9epQGDRrkSZxPokj2vGbXkiVLmbfwBxz935PB/qLA0JISiQycSnf/Nox+43W9wxGiwHj9jbcI/vcmDn6vSmeDKFCi9i6hRnFLvp8/L8++m+nN3//wtrZu3YqXlxdXrlwhODiYkSNHcv/+/ZTXIyMj6dq1K9OnT8fZ2TlP4nwSUrw+xr59+xg+ajSOnWSwvyhYovYsxK20A/PmztY7FCEKjG/mzeN/P/2CQ4exKIsieXJRFFBRx7dhe+UIm9evzdOLasuXL8+lS5dSHl++fJmyZcumabNo0SK6dOmCUorq1atTpUoVTp8+DUB8fDxdu3alT58+BfbOiFK8ZiI0NBT/Fzrj0HoU1k9X0jscIVJEBQficPMkm9b9jKWl/EALAfDrr7/yznsf4BTwARa2jnqHI0SKmIt/EfvnMnZs2cxTTz2Vp9vy9fXl3LlzXLhwgbi4OFasWEFAQECaNhUrVmTnzp0AXL9+nTNnzlC1alU0TWPQoEG4ubnx1ltv5WmcT0J+9TKQMtjf60XsqhWpO0qKAi46NJj4QyvZeeBPihUrpnc4QhQIZ8+epXO3Hji2fxurEuX0DkeIFPF3rxAZOIW1q5bj6uqa59uztLRk9uzZtG3blsTERAYOHIi7uzvz5s0DYOjQoXz00UcMGDAADw8PNE3jiy++oGTJkuzdu5elS5fi4eGBl5cXAJMmTaJDhw55Hnd2qPTGRmSksN0nOyOJiYm0ae/PsbsGHFoMkzFTosCIv32Z8NXvs/GXNfj5+ekdjlnJ7D7ZhVVRydl3796lbj1fHri2x7GuXLgoCo6kmEjCV77LxA/fYeSIEXqHY1Yyy9kybCAdb7z1NkcvXMfB7zUpXEWBkRgdQcSGiUz7arIUrkKYxMfH0yGgMw9Ke0rhKgoU40W1U+jVuaMUrrlMiteHfDt/PktWrpHB/qJA0RITiNr8Jf17deO1V1/VOxwhCgRN03ht6HDO3IrBoekAvcMRIo2o3QuoU64Yc2YWzLtUmTOpzlLZtWsXY8aOw6XH51jYOekdjhCA8Qc6ate31Kv6DNOnfqV3OEIUGNO+nsHPW37FpcdklMFC73CESBF1dBNOd86wYftBuag2D8g7anLu3Dle7Nodh/ZjZLC/Gbm6+A3irv+Tb9sr0Xooyto+X7ZlsLLFrmZDHhxZT/HIf1m7az8WFvIDLQTA5s2bGf/ZRIr1+gKDTf78TQqRFdEXjhAftJodB/fj4iJTbOYFKV4xDvZv2bY91g16Y1eprt7hiGzIz8IVwD3+LCWLlcqXbW3bGkj8ncuokC3sCDqIk5OcDRAC4OTJk/Tq0xfHju9h6VJa73CESBF/6xKRW6ezZcMvVKtWTe9wCq0iX7zGx8fj/0IXHpSSwf7mqNK7G9M81hLiibv1L3HXzmG4fZ74sBBc7K0Z88br9O/Xz6yOgmvU9uTfA6vYs+tXKlWSeYaFALh58yat2nXApskr2JavrXc4QqRIjL5PxIaJzJo2heeff17vcAq1Il28aprG0OEjOXXjAU4Bb+odjsgFytIKm2eqY/NMdcD4GcdeOsGEBT/z3vsf0qNHD8aMHoWHh4fOkT5eh/bt8PV+l4YNG+odihAFQmxsLG07BhBfqSGO7s31DkeIFFpiPJEbv2DQyz0ZOPAVvcMp9Ip08Tp9xix+2rwdlx5fyGD/QkophW1FD6jogW3kHTb8tY2fmrbg5d69mPrVZBwcHPQOMUMzpn6pdwhCFBiaptHvlUH8+8ASx+Z99A5HiBSaphH56zzq16rA1K8kb+eHIjtVVmBgIB99+hmOnT6Qwf5FhKVjCZwa9qJ4v1n8tO8UNd3q8Pvvv+sdlhAiC/5v0uds3XsIhzajUarI/nSJAuhB0C88FR3GmlXLMRjku5kfimTPa0hICD1fehnHju9hVewZvcMR+czCzgnHdm/y4NwB2r/QlZd79SzwvbBCFGU///wzX0ybQbFeX2KwttU7HJFFiVF3IZ2beCbFPQAtCXLxjKcyWGCwc8619WVle8rSmgfnDpD010Z2BB3E0dEx37Zf1BW54vXWrVvGwf6N+8tg/yLOvkYDbMrXZvWub9nbuCl7dm7jqaee0jssIUQqR48epd/AwTi98DGWTiX1Dkdkw+XZffN1ezZ2+XgW1WBF8c4f8mDHbHZuC6RixYr5t21RtIrX2NhY2vm/QFyFBjjWaal3OKIAsLBzwrH9GK7tXYLPc435Y8+vlC1bVu+whBDA1atXadPBHzu/IdiUqaF3OCKbHp4NJjVN09DiY4i7EUrc9b+xuHOB+CuncbKz5q1RI3hlwACKFSuWj9Fmj6NzMcLXf86Cb+bQoEEDvcMpcorM4AxN0xgw6FUuRBpwaPKy3uGIAkQphePz/Yko/xyNnvfj5s2beockRJEXHR1N6/b+aDVb4ODaRO9wRC5TSmGwtsO2vBvOz3bCofUoXPrNIb7hYCYtWk/ZCpXoP3AwN27c0DvUdNnZ2/HmqBG89FJvvUMpkopMz+ukyV8Q+NsBnLtNksH+Il2ODboTGR9D0xatOfDHbzg759/4KSHEf5KSkujZpy9XtWI4PtdT73BEPlFKYVuhDlSog23UXTYF/cI6tzp8O3cWPXsWrO/Bv+f/wc7OTu8wiqwiUcWtXbuWz6d8jaP/+zLYX2TKofHL3LAuw6Ahw/QORYgi68OPx/P70dM4th6JUkrvcIQOLByK49jsFew6vMtro9/F/4UuBeqsmL29vXw3dVToi9fg4GD6vjIIJ//3sHSWwf4ic0opHJoNZuvOPaxfv17vcIQocpYt+5FZ8xbg6P8eytJa73CEzmzKuVKsz9fsv2HA26cBoaGheockCoBCXbxeu3bNONi/2WvYlKmpdzjCTBisbbFvNZJXXh3C3bt3S/AN4AAAHD9JREFU9Q5HiCLjwIEDDBk+EqeAD7BwKK53OKKAUJbWODZ7hdhabfFt2Jh//vlH75CEzgpt8Zo82D+pRnMc3OQewyJ7bCt6oCrX57XhI/UORYgi4eLFi7Tv9AL2rUdiXaqK3uGIAsihnj/U7UwTvxaEhYXpHY7QUaEsXjVN46W+/bmS5IyDDPYXOWTfuB9bd+5h+/bteociRKEWGRlJq3YdMXh0xL66TDskMubg1Z64qk15sVtPkpKS9A5H6KRQFq8fjf+E3UEncWglg/1FzhmsbbH27cFnk7/SOxQhCq3ExES6dO/FbbtyOPh01jscYQYcG3Tnn+vhzJw1W+9QhE4KXfG6fPkKZsydj6P/OAxWNnqHI8ycvWsTDgcFyRgrIfLImHfGcujsZRxbDJXOBpElymCBXavX+fDj8ZKbi6hCVbwePHiQV4cNxyngfSwcZbC/eHIGKxvs3Fswa85cvUMRotBZ8P33LFy2CseOY1EWVnqHI8yIVYly2Pp2o+/AwXqHInRQaIrXS5cuGQf7txyJdamqeocjChEbj3YsXLSY6OhovUMRotD47bffeGPMWJw6fYCFvYve4Qgz5FCvE3+dCOH48eN6hyLyWaG4w1ZUVBSt2nUk7qkaONg5EnP5ZK6uP+7q39z99btcXWd22JeugkO56jle/uaR/y44etq7FRSqU3MGkixtoVR1HGo1zpN5Ia2Kl8H6meqsXbuWl156KdfXL0RR9PWsOURH3CNu6eg8WX9CXGyerDernGrUx7p4mRwtm/AgnPATuwGwK1cL+3KuuRhZAaAUCRZ2WJetiV2VeiiDRc5WY7DA2r0N02bMYvH3+v1Gi/ynNE3LcmMfHx8tKCgoD8PJmaioKFq27cCD6Jg8Wf/1K5e5ce1Knqw7qxYvXpzjZQcMGPDIc998802huLVdYmIi4eHh/LxhM0cOH8bOvSV2Pp1zvSfn/qF1dK6q+P67b3N1vSL/KKUOa5rmo3cc+amg5myAhIQE4uPj82z9Y8aMYd26dRm+rpTCxsYGa2vrdP9rZWWFjY1NynOpX7e2tsbSMvO+n2LFimEw5Ozk5q1bt/jss8/SPFemTBnee++9HK2voElMTOT2nTus3bCZS1euYeXeBgevDhhs7LO9roSI29xZ8jpXL1/ExUV68AuTzHJ2oSheRdZomsYXX3zB/7d3/1FR1/kex58DI0hiP/yFEfmDUEiZEZTxqLFGbUGhV9fSkrrbD3LV1OrYbUu9+aMf29WWXT1qrRp5Lcqgu9ZqFga6IoreaERQIxVv4g90Wy1/jYrAMPcPjrORomIzDMO8Hudwjt/h+/nMe/jjc15+5vP9fKZOnQrUHeIQEhLi4apcZ+/evcx+M42PV35O2+HTadXuFpf1XXlwJ+12fkxpyVaX9SlNS+FVvNGRI0e4//77KSkpoWvXri3uhCmr1cprb8wm76ti2g6ffk0nYdo+n8Mbz/w7Y8eOdUOF4imXG7NbzJpXuTKDwcCUKVOcazc7d+5MZaV7Zqs9ISIigvQli5g9axqnPn2FGtuPLus7IOQ2/m/Pt26dKRIR+bmbb76Z4uJiNm7cyP79+zGZTJ4uyaXi4uL424r/4YWnn+T0pzOxV9oa3Ycj1MT6/AI3VCfNlcKrD2rdujVnz54FYMCAAR6uxvUmPP00Tz/1GJUbl7qsT7+AIK67KYTS0lKX9SkicrXi4+PJzc1l586drFmzxtPluJTBYGD6y//J6BFDOZff+HE7oHMEhV/rGwZfovDqo4KCgpg6dSolJSU0ZumIt5g25SXOfleE3XbcZX0aQ25j61YtGxARz7jnnnsAuP/++z1ciXv84dVZnC37X+znTjWqXUDH7hzc/512hPEhCq8+7PXXXwdgzpw5Hq7E9W688UZGjhzJ2R05Luuz+rqOHDhwwGX9iYg0Vk5O3Zj244+uWxbVXHTo0IGhQ/+NszvXNqqdwdiKtiFd2LFjh5sqk+ZG4dWHXXgS9sIDXC3NY4+m4HfEdYOZwRjAGf3PXkQ86N577wVg7ty5Hq7EPcY8+Rh+B4sa3c7Y5kaOH3fdN23SvCm8+rjw8JZ7oEPHjh2xnzvtsv4MxgBsZxReRcTzysrKPF2CW3Tu3Pnaxm3/Vlo24EMUXn1cmzZtPF2C2wQEBOCotbuuQ4Mfta7sT0TkGtlsjX8q3xsEBgbisNc0vqHBH7td47OvUHgVEREREa+h8CoiIiIiXkPhVURERES8hsKriIiIiHgNhVcRERER8RoKryIiIiLiNRReRURERMRrKLyKiIiIiNdQeBURERERr6HwKiIiIiJeQ+FVRERERLyGwquIiIiIeA2FVxERERHxGgqvIiIiIuI1FF5FRERExGsovIqIiIiI11B4FRERERGvofAqLY7D4aC2tpba2lq39e1wOFzet4iIiFyZwqu0OFarFX9/f6Kjo/G/PsRl/RpvCCF98SL8/f357eNPuqxfERFft2vXLsIje3F3UjK0at3o9o5WrXl60rN073k7K1d95oYKpTlReJUWJyYmhiee+h3Xh4bTOmGsy/oNui2OdgMeJLRLN16eNsVl/YqI+LqOHTvSvl07Tpx30Cb5pUa3bz34Kc61CaWmqorboyLdUKE0Jwqv0uK0atWKpe8sZuYLz3Ii8yUqD5X+4j5rq85xevUcwvkH24usREVFuaBSEREBaN++PQUb/s7we36F7bM3qDl19Krb1p4/i23Nn4m+uQ0lRV/Ts2dPN1YqzYHCq7RIBoOB5yc/x8cfvsfZz2dz5pu/X3NfNaf+yan/mUqypSeb8tbRvn17F1YqIiIAAQEBfPj+Ml6cNIYTmS9y/vDuK7apOfk9Jz+ewr8NjCb/77m0a9euCSoVT1N4lRYtOTmZLZvyMRb/FdumDByOxj3Edf7wbk5kvsSUZ8aSsWwpAQEBbqpUREQMBgPTprzEh//9DrZVr3OmdEOD91Ye+pYTmS8x/fmJLHv3HVq1atWElYonKbxKi9e7d2+2F1npWn0A2+o3qa2qvKp2Z0o3YFv1OsuXpTPlxd9jMBjcXKmIiAAMHz6cTRvWY7Aux7Z5+UU7vJz5Zj1nPv8vsj5Yxgv/8bzGZx+j8Co+oWPHjmzOX8+9Md049ddp1Jw61uC9DocD2+bl+G39iIL8PIYNG9aElYqICECfPn3YXmQl1LYH25o/UVt9HoejFtumD/Df9jFbNm5gyJAhni5TPEDhVXxGYGAgH33wPi+Mf4ITWS9y/siei+6prT6PLTuNW86Usb3Iitls9kClIiIC0LlzZ74qyOfOnp04veJlbJ+/SZeqcrYXWYmOjvZ0eeIhCq/iUwwGAy9Pm0pG+iJOr3yNM7s2OX9XY/uRU399mYSousEyJMR1e8SKiMi1CQoKYsXHmTw/9reMHBzDlo15dOrUydNliQcZPV2AyOX4+/tjMpmoqamhe/fuZGRkcOONN/7ifkeMGMHGbt1ITB6K7fghAsItnP7sDUY/MJzWgQG0bn35TbIHDRrE5s2bf3EdIiJyZQaDgenTpnq6DGkmNPMqzVpQUBDFxcXs3LmTdu3a8dZbb7ms79jYWLYXWbn51Lf8mDWNd99ewDtLFrNgwYIrtlVwFRER8QyFV/EaAwcOpKKiwnn9xz/+EYvFgtlsZubMmQCUl5cTFRXFmDFjiI6O5tFHH2Xt2rXccccd9OjRg8LCQgAKCwsZNGgQycnJtAk08ve1OTz00Cjy8vIYOnQoALNmzSI1NZWEhATCw8OZP3++872Dg4MByMvLIyEhgZEjRxIVFcWjjz7qfCr2iy++ICoqivj4eJ599llnvyIiInLtFF7FK9jtdtatW+d88j8nJ4eysjIKCwspLi5m69at5OfnA7B3716ee+45tm/fzq5du1i+fDmbNm0iLS2NN954A4CoqCjy8/PZtm0br7/+On/6058u+b67du3iyy+/pLCwkFdeeYXq6uqL7tm2bRvz5s2jtLSU7777joKCAiorKxk3bhzZ2dls2rSJo0ev/rQYERERaZjWvEqzdu7cOWJiYigvL6dfv37ce++9QF14zcnJITY2FgCbzUZZWRldunShe/fumEwmoG6P11//+tcYDAZMJhPl5eUAnDx5kscff5yysjIMBsMlQynAkCFDCAwMJDAwkE6dOvH9998TFhZW757+/fs7X7tQa3BwMOHh4XTv3h2AlJQUlixZ4vK/j4iIiK/RzKs0axfWvO7fv5+qqirnmleHw8HUqVMpLi6muLiYvXv38tRTTwF1W2Jd4Ofn57z28/OjpqYGgOnTp3PXXXexc+dOPvvsMyorL31wwU/78vf3d7a/0j0/31BbREREXEPhVbzCDTfcwPz580lLS6O6upqkpCSWLl2KzWYDoKKign/+859X3d/Jkye55ZZbAFi2bJnL642KiuK7775zzvRmZWW5/D1ERER8kcKreI3Y2Fj69OlDZmYmiYmJPPLIIwwcOBCTycTIkSM5ffr0Vff14osvMnXqVO644w7sdrvLaw0KCuLtt9/mvvvuIz4+npCQEG644QaXv4+IiIivMTTm6824uDiH1Wp1YznS1MxmMzt27NDX3G5gs9kIDg7G4XAwceJEevToweTJkz1dls8yGAxbHQ5HnKfraEoas1smg8HAkCFDWL16tadLEXGby43ZmnkVcZN33nmHmJgYevfuzcmTJxk3bpynSxIRaRIXthN0h27dunHs2LGrvr+4uJgvvvjCbfVI09NuAyJuMnnyZM20ioh4WHFxMVarleTkZE+XIi6imVfxWYMGDfJ0CSIiPuPo0aM8+OCDWCwWLBYLBQUFwL8OjYmNjWXQoEHs3r0bqNvf+4UXXsBkMmE2m+udfrhgwQL69u2LyWRi165dDfZTVVXFjBkzyMrKIiYmhqysLM6cOUNqaioWi4XY2FhWrlzZ9H8M+UU08yo+S0e8iog0neeee47JkycTHx/PgQMHSEpK4ttvv3UeGmM0Glm7di3Tpk1jxYoVLFmyhH379rFt2zaMRiM//vijs68OHTpQVFTE22+/TVpaGunp6Q328+qrr2K1Wlm4cCEA06ZN4+6772bp0qWcOHGC/v37c88999CmTRtP/WmkkRRexWcFBwdjs9nIy8tj1qxZdOjQgZ07d9KvXz8++OADDAaDp0sUEWkx1q5dS2lpqfP61KlTnD59usFDY9auXcv48eMxGuuiSrt27ZxtH3jgAQD69evHJ598Alz94TM5OTmsWrWKtLQ0ACorKzlw4AC333676z+0uIXCqwh1R7x+8803hIaGcscdd1BQUEB8fLynyxIRaTFqa2vZsmULQUFB9V5/5plnuOuuu/j0008pLy8nISEBqDuMpqFJhAuHw/z08JgLh8/8vJ+fczgcrFixgsjISNd8MGlyWvMqwr+OePXz83Me8SoiIq6TmJjo/Ooe6h6kgoYPjUlMTGTRokXOcPrTZQOX0lA/bdu2rbcPeFJSEgsWLHBuEblt27Zr/1DiEQqvIlzdMbAiInJ1zp49S1hYmPPnz3/+M/Pnz8dqtWI2m+nVqxeLFi0CGj40ZsyYMXTp0gWz2UyfPn1Yvnz5Zd+zoX7uuusuSktLnQ9sTZ8+nerqasxmM9HR0UyfPt09fwRxGx1S4ON8+ZCCn655TUtLc274PWnSJOLi4njiiSc8W6C4lA4pkJZChxSIL9AhBSIiIiLSIuiBLfFZNpsNgISEhHoL+3+6JktERESaF828ioiIiIjXUHgVERERj1uzZg2RkZFEREQwe/bsi36/a9cuBg4cSGBgoHOPVoDdu3cTExPj/Ln++uuZN29eU5YuTUzLBkRERMSj7HY7EydOJDc3l7CwMCwWC8OGDaNXr17Oe9q1a8f8+fP529/+Vq9tZGSkc9stu93OLbfcwogRI5q0fmlamnkVn/GHP/yB3r17YzabiYmJ4auvvmrw3jFjxjhPgunWrRvHjh1rqjJFRHxOYWEhERERhIeHExAQwOjRo1m5cmW9ezp16oTFYqFVq1YN9rNu3Tpuu+02unbt6u6SxYM08yo+YcuWLaxevZqioiICAwM5duwYVVVVDd6fnp7ehNWJiPi2iooKbr31Vud1WFjYZScYGpKZmUlKSoorS5NmSDOv4hOOHDlChw4dnIcRdOjQgdDQUNatW0dsbCwmk4nU1FTOnz8P1O1AoP0xRUSaxqX2Gm/oaNiGVFVVsWrVKkaNGuWqsqSZUnj1Udu3b+e1115jx44dALz22mvk5uZ6uCr3SUxM5ODBg/Ts2ZMJEyawYcMGKisreeKJJ8jKymLHjh3U1NTwl7/8xdOliohc0tGjR9mzZw9Q9x/yPXv2cPbsWQ9X5RphYWEcPHjQeX3o0CFCQ0Mb1Ud2djZ9+/YlJCTE1eVJM6Pw6qO2FhUzY8YMAjr3IODmHsyYMYNVX3zp6bLcJjg4mK1bt7JkyRI6duzIww8/zOLFi+nevTs9e/YE4PHHHyc/P9/DlYqIXNroR39LZFQUAEVFRURGRpL+7n97uCrXsFgslJWVsW/fPqqqqsjMzGTYsGGN6uOjjz7SkgEfoTWvPurJJx4juM11PPm7seBnZMYb/8W0KS95uiy38vf3dx5IYDKZeO+99zxdkojIVVv6zmLuSUrmx6BbsB8rJzlhIOPH/c7TZbmE0Whk4cKFJCUlYbfbSU1NpXfv3ixatAiA8ePH849//IO4uDhOnTqFn58f8+bNo7S0lOuvv56zZ8+Sm5vL4sWLPfxJpCkovPqwUaNGEhUVyfHjxxk8eLCny3Gr3bt34+fnR48ePQAoLi4mJCSEkpIS9u7dS0REBBkZGdx5550erlRE5NK6du1K0df/y+8mTCLWnMiLL/xHo9eFNmfJyckkJyfXe238+PHOf3fu3JlDhw5dsu11113HDz/84Nb6pPlQePVxJpPJ0yU0CZvNxjPPPMOJEycwGo1ERESwZMkSUlJSGDVqFDU1NVgslnoDpYhIc9O2bVsyM/Stkfg2w6We8GtIXFycQ09gi4g3MhgMWx0OR5yn62hKGrNFxFtdbszWA1siIiIi4jUUXkVERETEayi8ioiIiIjXUHiVFm/NmjVERkYSERHB7NmzL/r9hx9+iNlsxmw2M2jQIEpKSpy/69atGyaTiZiYGOLifGq5pIiISLOk3QakRbPb7UycOJHc3FzCwsKwWCwMGzaMXr16Oe/p3r07GzZs4KabbiI7O5uxY8fWO1N7/fr1dOjQwRPli4iIyM9o5lVatMLCQiIiIggPDycgIIDRo0ezcuXKevcMGjSIm266CYABAwY0uI+giIiIeJ7Cq7RoFRUV3Hrrrc7rsLAwKioqGrz/3Xff5f7773deGwwGEhMT6devH0uWLHFrrSIiInJlCq/NWEJCAl9++WW91+bNm8eECRMueX95eTnR0dGX7bO8vJzly5c7r61WK88++ywAy5YtY9KkSQAsWrSI999//5eUD9StGT127Ngv7udaXWof44ZOpFm/fj3vvvsuc+bMcb5WUFBAUVER2dnZvPXWW+Tn57utVhHxfhq3RdxP4bUZS0lJITMzs95rmZmZpKSkXHOfPx8E4+LimD9//kX3jR8/nscee+ya36e5CAsL4+DBg87rQ4cOERoaetF927dvZ8yYMaxcuZL27ds7X79wb6dOnRgxYgSFhYXuL1pEvJbGbRH3U3htxkaOHMnq1as5f/48UDeAHT58mPj4eH7/+98THR2NyWQiKyvrorbl5eX86le/om/fvvTt25fNmzcDMGXKFDZu3EhMTAxz584lLy+PoUOHXtR+1qxZpKWlcfjwYWJiYpw//v7+7N+/n6NHj/Lggw9isViwWCwUFBQA8MMPP5CYmEhsbCzjxo275MxnU7JYLJSVlbFv3z6qqqrIzMxk2LBh9e45cOAADzzwABkZGfTs2dP5+pkzZzh9+rTz3zk5OVecIRER36ZxW8T9tNtAM9a+fXv69+/PmjVrGD58OJmZmTz88MN88sknFBcXU1JSwrFjx7BYLAwePLhe206dOpGbm0vr1q0pKysjJSUFq9XK7NmzSUtLY/Xq1QDk5eVdtobQ0FCKi4sBeOutt9iwYQNdu3blkUceYfLkycTHx3PgwAGSkpL49ttveeWVV4iPj2fGjBl8/vnnHl8najQaWbhwIUlJSdjtdlJTU+nduzeLFi0C6mYqXn31VX744Qfn13pGoxGr1cr333/PiBEjAKipqeGRRx7hvvvu89hnEZHmT+O2iPspvDZzF76CujAILl26lIyMDFJSUvD39yckJIQ777yTr7/+GrPZ7GxXXV3NpEmTKC4uxt/fnz179vyiOgoKCkhPT2fjxo0ArF27ltLSUufvT506xenTp8nPz+eTTz4BYMiQIc6n+D0pOTmZ5OTkeq+NHz/e+e/09HTS09MvahceHl5vz1cRkauhcVvEvRRem7nf/OY3PP/88xQVFXHu3Dn69u17VQvy586dS0hICCUlJdTW1tK6detrruHIkSM89dRTrFq1iuDgYABqa2vZsmULQUFBF93f0ANRIiK+QOO2iHtpzWszFxwcTEJCAqmpqc4F/4MHDyYrKwu73c7Ro0fJz8+nf//+9dqdPHmSm2++GT8/PzIyMrDb7QC0bdvWuY7zalRXV/PQQw8xZ86ceutBExMTWbhwofP6wldUgwcP5sMPPwQgOzub48ePX9sHFxHxUhq3RdxL4dULpKSkUFJSwujRowEYMWIEZrOZPn36cPfdd/Pmm2/SuXPnem0mTJjAe++9x4ABA9izZw9t2rQBwGw2YzQa6dOnD3Pnzr3ie2/evJmvv/6amTNnOhf/Hz58mPnz52O1WjGbzfTq1cu5hnTmzJnk5+fTt29fcnJy6NKli4v/GiIizZ/GbRH3MTTmqcK4uDiH1Wp1YzkiIu5hMBi2OhyOOE/X0ZQ0ZouIt7rcmK2ZVxERERHxGgqvIiIiIuI1FF593Jo1a4iMjCQiIoLZs2d7uhwRERGRy1J49WF2u52JEyeSnZ1NaWkpH330Ub09AEVEpHnRhIOIwqtPKywsJCIigvDwcAICAhg9ejQrV670dFkiInIJmnAQqaPw6sMqKiq49dZbnddhYWFUVFR4sCIREWmIJhxE6ii8+rBLbZOmU1ZERJonTTiI1FF49WFhYWEcPHjQeX3o0CFCQ0M9WJGIiDREEw4idRRefZjFYqGsrIx9+/ZRVVVFZmYmw4YN83RZIiJyCZpwEKmj8OrDjEYjCxcuJCkpidtvv52HHnqI3r17e7osERG5BE04iNQxeroA8azk5GSSk5M9XYaIiFzBTycc7HY7qampmnAQn6TwKiIi4iU04SCiZQMiIiIi4kUUXkVERETEayi8ioiIiIjXUHgVEREREa+h8CoiIiIiXkPhVURERES8hsKriIiIiHgNhVcRERER8RoKryIiIiLiNRReRURERMRrKLyKiIiIiNdQeBURERERr6HwKiIiIiJeQ+FVRERERLyGwquIiIiIeA2FVxERERHxGgqvIiIiIuI1DA6H4+pvNhiOAvvdV46IiNt0dTgcHT1dRFPSmC0iXqzBMbtR4VVERERExJO0bEBEREREvIbCq4iIiIh4DYVXEREREfEaCq8iIiIi4jUUXkVERETEayi8ioiIiIjXUHgVEREREa+h8CoiIiIiXkPhVURERES8xv8DtHElbmFo9tcAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "LandApp.plot_sankey_report(margin=0.5)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
