{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Stabilization"
   ]
  },
  {
   "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)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Model document\n",
    "\n",
    "The process includes mixing an additive into the soil to which the PFAS will preferentially sorb. The PFAS in the contaminated soil will either be released to surface or ground- water, and the remainder will sorbed to the soil/additive mixture. The flow of PFAS from soil stabilization 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 through the soil. Model predictions are based on achievement of equilibrium. It is further assumed that the soil and additive are well mixed. The partition coefficient is used to estimate the concentration of PFAS in the liquid and solids. By default, it is assumed that no volatilization occurs, but a user can enter a fraction of PFAS that is volatilized.\n",
    "\n",
    "The concentration in the liquid changes throughout the year as PFAS 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 location 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; however, this is outside the scope of the current effort and will depend on how the land is used.\n",
    "\n",
    "<img src=\"../Images/ProcessModels/Stabilization.png\" alt=\"Drawing\" style=\"width: 400px;\"/>\n",
    "\n",
    "### Assumptions and Limitations:\n",
    "\n",
    "1.\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",
    "2.\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",
    "3.\tFuture work is required to implement a dynamic (i.e., non-equilibrium) 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",
    "4.\tDefault soil properties reflect representative values for loam and silty clay loams, but the values for different situations are provided and may be input by the user.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for Stabilization model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": true
   },
   "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>Reference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log Koc (soils)</td>\n",
       "      <td>2.19</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log Koc (soils)</td>\n",
       "      <td>3.04</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log Koc  (soils)</td>\n",
       "      <td>1.88</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log Koc (soils)</td>\n",
       "      <td>1.37</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log Koc (soils)</td>\n",
       "      <td>1.77</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log Koc (soils)</td>\n",
       "      <td>1.97</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log Koc  (soils)</td>\n",
       "      <td>2.63</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log Koc  (soils)</td>\n",
       "      <td>3.24</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log Koc (soils)</td>\n",
       "      <td>1.51</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Soil Log partition coefficient</td>\n",
       "      <td>SoilLogPartCoef</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log Koc (soils)</td>\n",
       "      <td>2.79</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>Stabilization</td>\n",
       "      <td>Stabil</td>\n",
       "      <td>depth_mix</td>\n",
       "      <td>Depth of contamination to be treated</td>\n",
       "      <td>1.00</td>\n",
       "      <td>m</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Stabilization</td>\n",
       "      <td>Stabil</td>\n",
       "      <td>appl_dens</td>\n",
       "      <td>Mass applied per area</td>\n",
       "      <td>20.00</td>\n",
       "      <td>kg/m2</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Stabilization</td>\n",
       "      <td>Stabil</td>\n",
       "      <td>bulk_dens</td>\n",
       "      <td>Bulk density of contaminated soil to be treated</td>\n",
       "      <td>1175.00</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>Stabilization</td>\n",
       "      <td>Stabil</td>\n",
       "      <td>add_mass_ratio</td>\n",
       "      <td>Mass of additive per mass of soil</td>\n",
       "      <td>0.05</td>\n",
       "      <td>kg add./kg soil</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Stabilization</td>\n",
       "      <td>Stabil</td>\n",
       "      <td>is_stay_inplace</td>\n",
       "      <td>Is the stabilized soil stay in place</td>\n",
       "      <td>0.00</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Additive Properties</td>\n",
       "      <td>AddProp</td>\n",
       "      <td>ts_cont</td>\n",
       "      <td>Additive total solids content -wet</td>\n",
       "      <td>1.00</td>\n",
       "      <td>fraction wet weight</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Additive Properties</td>\n",
       "      <td>AddProp</td>\n",
       "      <td>bulk_dens</td>\n",
       "      <td>Additive bulk density - wet</td>\n",
       "      <td>1000.00</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log K (Additive)</td>\n",
       "      <td>3.90</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log K (Additive)</td>\n",
       "      <td>4.30</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log K  (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log K (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log K (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log K (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log K  (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log K  (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log K (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>Additive Log partition coefficient</td>\n",
       "      <td>AddLogPartCoef</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log K (Additive)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>Precipitation Data</td>\n",
       "      <td>Precip</td>\n",
       "      <td>ann_precip</td>\n",
       "      <td>Annual precipitation</td>\n",
       "      <td>1.13</td>\n",
       "      <td>m/yr</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>Precipitation Data</td>\n",
       "      <td>Precip</td>\n",
       "      <td>frac_runoff</td>\n",
       "      <td>Fraction of incoming precipitation to run-off</td>\n",
       "      <td>0.30</td>\n",
       "      <td>frac</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>Precipitation Data</td>\n",
       "      <td>Precip</td>\n",
       "      <td>frac_ET</td>\n",
       "      <td>Fraction of incoming precipitation to ET</td>\n",
       "      <td>0.65</td>\n",
       "      <td>frac</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</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.00</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                              Category   Dictonary_Name   Parameter Name  \\\n",
       "0   Soil Log partition coefficient      SoilLogPartCoef  PFOA              \n",
       "1   Soil Log partition coefficient      SoilLogPartCoef  PFOS              \n",
       "2   Soil Log partition coefficient      SoilLogPartCoef  PFBA              \n",
       "3   Soil Log partition coefficient      SoilLogPartCoef  PFPeA             \n",
       "4   Soil Log partition coefficient      SoilLogPartCoef  PFHxA             \n",
       "5   Soil Log partition coefficient      SoilLogPartCoef  PFHpA             \n",
       "6   Soil Log partition coefficient      SoilLogPartCoef  PFNA              \n",
       "7   Soil Log partition coefficient      SoilLogPartCoef  PFDA              \n",
       "8   Soil Log partition coefficient      SoilLogPartCoef  PFBS              \n",
       "9   Soil Log partition coefficient      SoilLogPartCoef  PFHxS             \n",
       "10  Stabilization                       Stabil           depth_mix         \n",
       "11  Stabilization                       Stabil           appl_dens         \n",
       "12  Stabilization                       Stabil           bulk_dens         \n",
       "13  Stabilization                       Stabil           add_mass_ratio    \n",
       "14  Stabilization                       Stabil           is_stay_inplace   \n",
       "15  Additive Properties                 AddProp          ts_cont           \n",
       "16  Additive Properties                 AddProp          bulk_dens         \n",
       "17  Additive Log partition coefficient  AddLogPartCoef   PFOA              \n",
       "18  Additive Log partition coefficient  AddLogPartCoef   PFOS              \n",
       "19  Additive Log partition coefficient  AddLogPartCoef   PFBA              \n",
       "20  Additive Log partition coefficient  AddLogPartCoef   PFPeA             \n",
       "21  Additive Log partition coefficient  AddLogPartCoef   PFHxA             \n",
       "22  Additive Log partition coefficient  AddLogPartCoef   PFHpA             \n",
       "23  Additive Log partition coefficient  AddLogPartCoef   PFNA              \n",
       "24  Additive Log partition coefficient  AddLogPartCoef   PFDA              \n",
       "25  Additive Log partition coefficient  AddLogPartCoef   PFBS              \n",
       "26  Additive Log partition coefficient  AddLogPartCoef   PFHxS             \n",
       "27  Precipitation Data                  Precip           ann_precip        \n",
       "28  Precipitation Data                  Precip           frac_runoff       \n",
       "29  Precipitation Data                  Precip           frac_ET           \n",
       "30  Volatilization                      Volatilization   frac_vol_loss     \n",
       "\n",
       "                              Parameter Description   amount  \\\n",
       "0   PFOA Log Koc (soils)                             2.19      \n",
       "1   PFOS Log Koc (soils)                             3.04      \n",
       "2   PFBA Log Koc  (soils)                            1.88      \n",
       "3   PFPeA Log Koc (soils)                            1.37      \n",
       "4   PFHxA Log Koc (soils)                            1.77      \n",
       "5   PFHpA Log Koc (soils)                            1.97      \n",
       "6   PFNA Log Koc  (soils)                            2.63      \n",
       "7   PFDA Log Koc  (soils)                            3.24      \n",
       "8   PFBS Log Koc (soils)                             1.51      \n",
       "9   PFHxS Log Koc (soils)                            2.79      \n",
       "10  Depth of contamination to be treated             1.00      \n",
       "11  Mass applied per area                            20.00     \n",
       "12  Bulk density of contaminated soil to be treated  1175.00   \n",
       "13  Mass of additive per mass of soil                0.05      \n",
       "14  Is the stabilized soil stay in place             0.00      \n",
       "15  Additive total solids content -wet               1.00      \n",
       "16  Additive bulk density - wet                      1000.00   \n",
       "17  PFOA Log K (Additive)                            3.90      \n",
       "18  PFOS Log K (Additive)                            4.30      \n",
       "19  PFBA Log K  (Additive)                           0.00      \n",
       "20  PFPeA Log K (Additive)                           0.00      \n",
       "21  PFHxA Log K (Additive)                           0.00      \n",
       "22  PFHpA Log K (Additive)                           0.00      \n",
       "23  PFNA Log K  (Additive)                           0.00      \n",
       "24  PFDA Log K  (Additive)                           0.00      \n",
       "25  PFBS Log K (Additive)                            0.00      \n",
       "26  PFHxS Log K (Additive)                           0.00      \n",
       "27  Annual precipitation                             1.13      \n",
       "28  Fraction of incoming precipitation to run-off    0.30      \n",
       "29  Fraction of incoming precipitation to ET         0.65      \n",
       "30  Fraction of PFAS lost to volatilization          0.00      \n",
       "\n",
       "                   unit  Reference  \n",
       "0   log L/kg OC         NaN         \n",
       "1   log L/kg OC         NaN         \n",
       "2   log L/kg OC         NaN         \n",
       "3   log L/kg OC         NaN         \n",
       "4   log L/kg OC         NaN         \n",
       "5   log L/kg OC         NaN         \n",
       "6   log L/kg OC         NaN         \n",
       "7   log L/kg OC         NaN         \n",
       "8   log L/kg OC         NaN         \n",
       "9   log L/kg OC         NaN         \n",
       "10  m                   NaN         \n",
       "11  kg/m2               NaN         \n",
       "12  kg/m3               NaN         \n",
       "13  kg add./kg soil     NaN         \n",
       "14  1:TRUE,0:FALSE      NaN         \n",
       "15  fraction wet weight NaN         \n",
       "16  kg/m3               NaN         \n",
       "17  log L/kg            NaN         \n",
       "18  log L/kg            NaN         \n",
       "19  log L/kg            NaN         \n",
       "20  log L/kg            NaN         \n",
       "21  log L/kg            NaN         \n",
       "22  log L/kg            NaN         \n",
       "23  log L/kg            NaN         \n",
       "24  log L/kg            NaN         \n",
       "25  log L/kg            NaN         \n",
       "26  log L/kg            NaN         \n",
       "27  m/yr                NaN         \n",
       "28  frac                NaN         \n",
       "29  frac                NaN         \n",
       "30  fraction            NaN         "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Stabilization = ps.Stab()\n",
    "Stabilization.InputData.Data[['Category','Dictonary_Name','Parameter Name', 'Parameter Description', 'amount', 'unit','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming Contaminated Soil to Stabilization"
   ]
  },
  {
   "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>850</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>150</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>8.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1e+06</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   850  \n",
       "2   Moisture flow  kg   150  \n",
       "3   Carbon flow    kg   8.5  \n",
       "4   PFOA           μg   1e+06\n",
       "5   PFOS           μg   1e+06\n",
       "6   PFBA           μg   1e+06\n",
       "7   PFPeA          μg   1e+06\n",
       "8   PFHxA          μg   1e+06\n",
       "9   PFHpA          μg   1e+06\n",
       "10  PFNA           μg   1e+06\n",
       "11  PFDA           μg   1e+06\n",
       "12  PFBS           μg   1e+06\n",
       "13  PFHxS          μg   1e+06"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncominWaste = ps.IncomFlow()\n",
    "IncominWaste.set_flow('Contaminated Soil', 1000)\n",
    "IncominWaste.calc()\n",
    "ContaminatedSoil = IncominWaste.Inc_flow\n",
    "ContaminatedSoil.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in Stabilization (Stabilized soil stay in place)"
   ]
  },
  {
   "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>99.92</td>\n",
       "      <td>0.01</td>\n",
       "      <td>0.07</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>99.97</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>67.14</td>\n",
       "      <td>4.69</td>\n",
       "      <td>28.17</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>43.04</td>\n",
       "      <td>8.14</td>\n",
       "      <td>48.82</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>61.85</td>\n",
       "      <td>5.45</td>\n",
       "      <td>32.70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>71.26</td>\n",
       "      <td>4.11</td>\n",
       "      <td>24.64</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>91.58</td>\n",
       "      <td>1.20</td>\n",
       "      <td>7.22</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>97.78</td>\n",
       "      <td>0.32</td>\n",
       "      <td>1.91</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>49.24</td>\n",
       "      <td>7.25</td>\n",
       "      <td>43.51</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>94.00</td>\n",
       "      <td>0.86</td>\n",
       "      <td>5.14</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Volatilized  Remaining in Soil  Leachate  RunOff\n",
       "PFOA   0.0          99.92              0.01      0.07  \n",
       "PFOS   0.0          99.97              0.00      0.03  \n",
       "PFBA   0.0          67.14              4.69      28.17 \n",
       "PFPeA  0.0          43.04              8.14      48.82 \n",
       "PFHxA  0.0          61.85              5.45      32.70 \n",
       "PFHpA  0.0          71.26              4.11      24.64 \n",
       "PFNA   0.0          91.58              1.20      7.22  \n",
       "PFDA   0.0          97.78              0.32      1.91  \n",
       "PFBS   0.0          49.24              7.25      43.51 \n",
       "PFHxS  0.0          94.00              0.86      5.14  "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Stabilization.InputData.Stabil['is_stay_inplace']['amount'] = True\n",
    "Stabilization.calc(Inc_flow=ContaminatedSoil)\n",
    "Stabilization.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deVwV9f7H8ddZFFkUU9RUzA0VBQ6L4lVc0grUVHJNyZ9p6L25l6mVdUvtXtcsTbEsW1xSoVJTSQ2XXNMQARURpRQXpHJFdg6H+f3B7RSKCslhDvB5Ph4+HvfM+c7MZ7j6bpj5LhpFURBCCFH2tGoXIIQQlZUEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESvQlaezk5KQ0adLEQqUIIUTFdOzYsWuKotS5c3uJArhJkyZERUWVXlVCCFEJaDSaC0Vtl0cQQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJICFEEIlEsBCCKESCWAhhFCJBLAQQqhEAlgIIVQiASyEECqRABZCCJVIAAshhEokgIUQQiUSwKIQBwcHix27SZMmXLt2rdjtY2Nj2bZtm8XqEUJtEsDCakkAi4pOAlg80NWrVxk4cCC+vr74+vpy6NAhACIjI/Hz88Pb2xs/Pz/OnDkDgMlkYurUqXh4eGAwGFi6dKn5WEuXLsXHxwcPDw8SEhLueZzc3FzefvttwsLC8PLyIiwsjIyMDIKDg/H19cXb25vNmzeX/Q9DiNKkKEqx/7Rt21YRFZu9vf1d24KCgpQDBw4oiqIoFy5cUFxdXRVFUZTU1FTFaDQqiqIoO3fuVAYMGKAoiqJ8+OGHyoABA8zfXb9+XVEURWncuLGyZMkSRVEUZdmyZcqoUaPue5wvvvhCGT9+vLmO6dOnK2vWrFEURVFu3ryptGjRQklPTy/FqxfCMoAopYhMLdGinKJy2rVrF/Hx8ebPt2/fJi0tjdTUVEaMGEFiYiIajQaj0WhuP2bMGPT6gr9etWrVMu87YMAAANq2bcvGjRsB7nmcO0VERLBlyxYWLlwIQHZ2NhcvXqR169alf9FClAEJYPFA+fn5HD58GFtb20LbJ06cSPfu3dm0aRNJSUl069YNKPitSqPRFHksGxsbAHQ6HXl5eQC89dZbRR7nToqisGHDBlq1alU6FyaEyuQZsHiggIAAQkJCzJ9jY2OBgjvXhg0bArBy5cpC7ZcvX24O2Bs3btz3+Pc6TvXq1UlLSzN/7tGjB0uXLqXgNzqIiYn5+xclhBWQABaFZGZm4uzsbP7z/vvvs2TJEqKiojAYDLRp04bly5cD8OqrrzJ9+nQ6deqEyWQyH2P06NE89thjGAwGPD09Wbdu3X3Pea/jdO/enfj4ePNLuLfeeguj0YjBYMDd3Z233nrLMj8EIcqI5o+7ieJo166dEhUVZcFyhBCi4tFoNMcURWl353a5AxZCCJVIAAshhEokgEWJ7dixg1atWuHi4sK8efPu+l5RFCZNmoSLiwsGg4Ho6GigoNtY+/bt8fT0xM3NjRkzZpj3mTZtGq6urhgMBvr378+tW7fK7HqEUIsEsCgRk8nE+PHj2b59O/Hx8axfv75QH2GA7du3k5iYSGJiIp988gljx44FCrqg7dmzh+PHjxMbG8uOHTs4cuQIAP7+/sTFxXHixAlatmzJ3Llzy/zahChrEsCiRCIjI3FxcaFZs2ZUrVqVoUOH3jUkePPmzTz//PNoNBo6dOjArVu3SElJQaPRmCf7MRqNGI1Gc3/hgIAA88CNDh06cPny5bK9MCFUIAEsSiQ5OZlGjRqZPzs7O5OcnFzsNiaTCS8vL+rWrYu/vz//+Mc/7jrH559/Tq9evSx0BUJYDwlgUSJFdVu8c9Tb/drodDpiY2O5fPkykZGRxMXFFWo3e/Zs9Ho9w4YNK8WqhbBOEsCiRJydnbl06ZL58+XLl2nQoEGJ29SsWZNu3bqxY8cO87ZVq1YRHh7O2rVr7zmUWYiKRAJYlIivry+JiYmcP3+e3NxcQkNDCQwMLNQmMDCQ1atXoygKR44cwdHRkfr163P16lVz74asrCx27dqFq6srUNCzYv78+WzZsgU7O7syvy4h1CCT8YgS0ev1hISE0KNHD0wmE8HBwbi5uZmHJ48ZM4ann36abdu24eLigp2dHV988QUAKSkpjBgxApPJRH5+Ps8++yx9+vQBYMKECeTk5ODv7w8UvIj745hCVFQyFFkIISxMhiILIYSVkQAWQgiVSAALIYRKJICFEEIlEsDC4iIiIti4cZPaZQhhdaQbmrAYRVGYO38Bcxa8h1ZflUOHD/Pu/HlotfLffSFAAlhYiKIoDB8ZTOj69dQdtgBtNQeWLHuJU/EJbP12A1WqVFG7RCFUJwEsLEKj0VDN1o5mzZtz89QujDY1cK7/KNVr1DDPeiZEZSe/CwqL+XT5Mqa8PAmUfMjPZ+Tz/8fX67+UeR6E+B8JYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALi8rKykLR6tHm51KtWjW1yxHCqkgAC4vae+BHFMf66DKvU6dOHbXLEcKqSAALi/n111/5/vsdVGvaloxz0QQGBqpdkhBWRQJYWITJZGL0mHE4uD9B9s9H6N27D05OTmqXJYRVkQAWpe7s2bP06vMMP8ZfpKpHT3KiN/PqlJfVLksIqyOLc6ngaNQx1q9f/9DHWfT+e6VQzb01a9acZ/r1K3b7zKwcYo6fIO5UHLbuAdh27UPaltn8Z8Zb+Pr6WrDS8i8jIwOj0UjNmjXVLkWUIQlgFRyLjmbx4sXU+McgNDa2f/s4esd65KX+VoqVFXbFZM8X0deL3V6j1aOr24mabZ4jO34Pt75+g/lz/sOE8eMtVmNFEThgMNevX+fIwX3SW6QS0SiKUuzG7dq1U6KioixYTuWgKAojgkcT/lM81ftMR6PVqV1Sqbgwv0+hz54+7WjW3EWlasqX/fv2km6EvgHd+Wr9Wlm4tILRaDTHFEVpd9d2CWB1GI1GunR/ip9NdXB4/AW1yykVN3Z+RH5OJjbObdDobdQup1zRaHXYNvfl9sa3mTx6GDPe+rfaJYlSJAFshW7cuIHBpx05bfpibwhQuxxhBfLSrnMr9FW+WB7C4MGDimyj0+nw8PAgLy+Ppk2bsmbNGos+O46KimL16tUsWbLkvu38/Pz48ccfLVZHeSYBbKXOnDlD+46dsO05hWqPGdQuR1iBnF9/Ju3bWezbvZO2bdve9b2DgwPp6ekAjBgxgpYtW/Lmm2+WdZmiBO4VwNINTWWtWrXim7D1pG9/D+PNK2qXI6yAzaMu2HUfS4/efUlOTr5v244dOxZq8+677+Lr64vBYGDGjBkAJCUl4erqyujRo3F3d2fYsGHs2rWLTp060aJFCyIjIwGIjIzEz88Pb29v/Pz8OHPmDAB79+6lT5+C5/szZ84kODiYbt260axZs0J3xQ4ODub23bp1Y9CgQbi6ujJs2DD+uNHbtm0brq6udO7cmUmTJpmPW1lJAFsBf39/5s9+h/Qts8nPTle7HGEF7Fr5QesA/Hv1ITMzs8g2JpOJ3bt3m0cYRkREkJiYSGRkJLGxsRw7doz9+/cD8PPPP/PSSy9x4sQJEhISWLduHQcPHmThwoXMmTMHAFdXV/bv309MTAzvvPMOb7zxRpHnTUhI4PvvvycyMpJZs2ZhNBrvahMTE8PixYuJj4/n3LlzHDp0iOzsbF588UW2b9/OwYMHuXr1amn8qMo1CWArMX7cOIL69yV920KUfJPa5QgrYN9+EL9pnXj2uf8jPz/fvD0rKwsvLy9q167NjRs38Pf3BwoCOCIiAm9vb3x8fEhISCAxMRGApk2b4uHhgVarxc3NjSeffBKNRoOHhwdJSUkApKamMnjwYNzd3Zk8eTKnTp0qsq7evXtjY2ODk5MTdevW5bff7u4K2b59e5ydndFqtXh5eZGUlERCQgLNmjWjadOmAAQFBZXmj6tckgC2IiFLFmFo9AgZez9VuxRhBTQaDQ5PjefQiUSmv/mWebutrS2xsbFcuHCB3Nxcli1bBhR0b5w+fTqxsbHExsby888/M2rUKABsbP7slaLVas2ftVoteXl5ALz11lt0796duLg4tm7dSnZ2dpF1/fVYOp3OvP+D2pTkfVNlIQFsRfR6PZs3fk31m2fJiPlO7XKEFdDoq+DQ+3U+/HQla9Z8Weg7R0dHlixZwsKFCzEajfTo0YPPP//c/IIuOTmZ33//vdjnSk1NpWHDhgCsXLmy1K7hD66urpw7d858xx0WFlbq5yhvJICtjKOjI7u/30besW/IOh+tdjnCCujsa1I98A3GTnyJw4cPF/rO29sbT09PQkNDCQgI4LnnnqNjx454eHgwaNAg0tLSin2eV199lenTp9OpUydMptJ/DGZra8uHH35Iz5496dy5M/Xq1cPR0bHUz1OeSDc0K3XgwAF69e2H46DZVHFqpHY5wgpk/XIU496PiImKpHHjxmqX87ekp6fj4OCAoiiMHz+eFi1aMHnyZLXLsjjphlbOdOnShaWL3iNt62xMmalqlyOsgG1zX3RegTzZ4+kS3dlakxUrVuDl5YWbmxupqam8+OKLapekKrkDtnKvTJ3Gyk07qT5gJhpdFbXLESpTFIWM3R/iVRu+37YVna5izCNS0ckdcDm1cMF8fFs1ImPPcnmLLNBoNNh3f5Hoc7/y8pRppXpsPz+/Uj2eeDAJYCun1WrZ8NV6amUlkxn1rdrlCCug0emxf/pVVoV+wycrVpTacWUeh7InAVwOODg4sPv7beSfDCcz8Se1yxFWQGdbneqB/+aVaa+zd+/eUjlmcYYSi9IlAVxONGrUiO3hW8jcvYzc38+pXY6wAlVqNcS+5yv0GziYn3/+uVSPXdRQYlH6JIDLkfbt27Pio2WkbZmDKf2m2uUIK2DbxIsq7YfyZI+nuXXrVqkdt6ihxKL0yZJEd/jrVH9l6e2336Zr16489dRT920XFDSUuPh4PvxyHjUG/geNvupdbS5/FIy2qi1oNOjsH8Gp9yvoHB75c7u24L+7tfzHUc25NQC3j37LzX2raDTxS7Q29gDkG7O5vmMpxt+TAAWtjQN1n51VcAxhNew9e5J+4xJ9+g1k767v0esf/p91cYYbi4cnAWwl3nnnnWK3/e87s4iLP83BnUtx6PlKkcvX1Auag87OkZv7VpF65CtqPfVioe13yji9D5v6Lcg8exgHj4L/CKRFbUFnV5M6owrmGjBev4xGK39lrJH948HEb5nNmHET+PST5WqXI4pJHkHcw/1eRBw9ehQ/Pz88PT1p3749aWlpZGdn88ILL+Dh4YG3tzc//PADUDCmvl+/fvTt25emTZsSEhLC+++/j7e3Nx06dODGjRsAjBw5km+++QaAJk2aMGPGDHx8fPDw8CAhIQGAq1ev4u/vT9u2bXGqWYO004dIO7DqvtdRrZH7A+cZNt5MQcnNpmaX4WTE7zNvN2XcRF+9tvlzldrOaPTSF9kaabQ6HHpN4evwCBZ/sFTtckQxSQDfR1EvInJzcxkyZAgffPABx48fZ9euXdja2ppnpDp58iTr169nxIgR5tmk4uLiWLduHZGRkbz55pvY2dkRExNDx44dWb16dZHndnJyIjo6mrFjx7Jw4UIAZs2axRNPPEF0dDSDBw/GaMxFe+4gGacP3PMasn6JpGqdJubPv61/gytfTCRl9SvmbRmn92HXuis2jdww3kjGlFHwLNHBw5/UnzaQsmYKN/evwXjj/pODC3VpbexxCHyTf898hx07dpR4/z8evXXr1o3w8HDz9pCQEEaOHFlaZYq/kAC+j6JeRJw5c4b69evj6+sLQI0aNdDr9Rw8eJDhw4cDBbM+NW7cmLNnzwLQvXt3qlevTp06dXB0dKRv374AheZivdOAAQMAaNu2rbnNwYMHGTp0KAA9e/bkkUce4ev168ja9wk5KWcL7f9H0ObnZOLYYbB5e72gOTR4YSn1n3/fvC3z9H7sW3dFo9Fi17IjmWcOAlC1XjMavvgpju0Hkp+dRsrqVzBeu/QwP1JhYVVqPopD72k8GzSM06dPq12OeAB5oHcf95rTtKhnrvfrJ1mcuVjvtc9fX4AUdQ4PDw/WfPEZw0ePQTdkAfoaTsC9n/XeKff38xhvXuG3sP/NN5ufh97xUar7FCwVo61qi10rP+xa+aHRaMg6d1QmB7Jy1ZzdyPMbjn+v3lxOki6L1kzugEvI1dWVK1eucPToUQDS0tLIy8uja9eurF27FoCzZ89y8eJFWrVqVarn7ty5M1999RVQsPrBzZsFXdH69+/P61NeJn3rHPJzi55E+14yTu+jZqfncB77ecGf8avJS79OXurvZF+Ox/S/JZIUk5Hca5fQ1ahbqtckLOTqOZo3d1G7CvEAEsAlVLVqVcLCwpg4cSKenp74+/uTnZ3NuHHjMJlMeHh4MGTIEFauXFnozrc0zJgxg4iICHx8fNi+fTv169enevXqALw5/XV6dm1PRsTiEh0z4/QB7Fp2LLTNrkUHMk7vJ+9WCr+te50rn40n5YuXChaLbNWp1K5HWEZGzHdUv3GGLRu/vm+72bNn4+bmhsFgwMvLi59+uvcoy9GjRxMfHw8UvCS+du1aqdZcWclsaOVITk4OOp0OvV7P4cOHGTt2LLGxsYW+7/T4EyTpG+HQebiKlQq1ZJ2PJmfXUo5FHqZ58+b3bHf48GFeeeUV9u7di42NDdeuXSM3N5cGDRo88BxNmjQhKioKJyen0iy9QpPZ0CqAixcv4uvri6enJ5MmTWLFHROx2NjYsH3rt1S5cJiMU3tUqlKoxXjtEhnfL2brtxvuG74AKSkpODk5mX9Lc3JyokGDBuzevRtvb288PDwIDg4mJycHKOgZITdfpU8CuBxp0aIFMTExHD9+nKNHj5p7YvxVnTp12PhVKDd3LEXJu3u5cFExmbJuk7Z1NkveX0iXLl0e2D4gIIBLly7RsmVLxo0bx759+8jOzmbkyJGEhYVx8uRJ8vLy+Oijj8qg+spLAriCycvL47U338bRu6cMmqgkFJOR9PD5jPq/IQQHv1CsfRwcHDh27BiffPIJderUYciQIXz88cc0bdqUli1bAjBixAj2799vydIrPemGZqWCg4MJDw+nbt26xMXFFXu/8ZNeIi75FtWfmWjB6oS1UBSF9D3Lad+qEe+9u6BE++p0Orp160a3bt3w8PBg1ar7j6oUpU/ugK3UyJEjSzyaKWTZMkI3fYdDr6lotLJUTWWQGfUttbOS2fDVerTa4v9zPnPmDImJiebPsbGx1KtXj6SkJPPUlmvWrOHxxx8v9ZrFn+QO2Ep17dq1RFMA7ty5k9f/PQPHZ+ehreZgucKE1chM/In8k+Hsioo0T6ZeXOnp6UycOJFbt26h1+txcXHhk08+ISgoiMGDB5OXl4evry9jxoyxUPUCpBuaVUtKSqJPnz4PfARx5swZ2nfshG2vaVRr5F5G1Qk15f5+ntsb3mZ3xHb+8Y9/qF2OeADphlZBXb9+nSd79KJqh2ESvpWEKeMmaVvn8MlHyyR8yzkJ4HIsNzeXpwP7k93AB3tDgNrliDKg5OWSHj6PSWNG89xzQWqXIx6SBHA5pSgKo/41hp9vmrDv/Lza5YgyoCgK6TuX0tWnNf99Z5ba5YhSIAFspYKCgujYsSNnzpzB2dmZzz77rND3C997n627DmDfc7L0eKgkMo6EUV+TSuiXq4ucke+vduzYQatWrXBxcWHevHl3ff/uu+/i5eWFl5cX7u7u6HQ6bty4wZkzZ8zbvby8qFGjBosXF8wvMnPmTBo2bGj+btu2bQAYjUZGjBiBh4cHrVu3Zu7cuaV/8RWUvIQrh8LDwxn6fDA1h85HL7OTVQoZCQfhp9WciI6ifv36921rMplo2bIlO3fuxNnZGV9fX9avX0+bNm2KbL9161YWLVrEnj177jpOw4YN+emnn2jcuDEzZ87EwcGBqVOnFmq3bt06tmzZQmhoKJmZmbRp04a9e/fSpEmTh7rmikRewlUQJ0+eJGj4CKr3fk3Ct5LISUkka+/HRGwLf2D4AkRGRuLi4kKzZs2oWrUqQ4cOZfPmzfdsv379eoKC7n6evHv3bpo3b07jxo3vez6NRkNGRgZ5eXlkZWVRtWpVatSo8eALExLA5cnvv/+Of6/eVOsSjE1DV7XLEWUg7/Y10sLnsuaLz/D29i7WPsnJyTRq9Oek+c7OziQnF72cVGZmJjt27GDgwIF3fRcaGnpXMIeEhGAwGAgODjbPRz1o0CDs7e2pX78+jz32GFOnTqVWrVrFvcRKTQK4nMjOzibg6b7kNe2MfZtuapcjykB+bjbpW+fw+pSX6d+/f7H3K+qx4r2eGW/dupVOnTrdFZi5ubls2bKFwYP/XM5q7Nix/PLLL8TGxlK/fn2mTJkCFNxx63Q6rly5wvnz53nvvfc4d05W4igOCeByQFEUho8M5lJONez9pOtRZaAo+WRELKZHF1/enP56ifZ1dnbm0qU/1+67fPnyPef5LeouF2D79u34+PhQr14987Z69eqh0+nQarX885//JDIyEih4BtyzZ0+qVKlC3bp16dSpk0xdWUwSwOXAf/47m12HY7APmIRGI/+XVQYZh9bSxN7E6i8+e2CPhzv5+vqSmJjI+fPnyc3NJTQ0lMDAwLvapaamsm/fPp555pm7vivquXBKSor5f2/atAl394KBP4899hh79uxBURQyMjI4cuQIrq7yiKw4ZC4IK/f119/w7gchOA5ZgLZKNbXLEWUg49Qeqlw4zI7oqL+1rJVeryckJIQePXpgMpkIDg7Gzc2N5cuXA5jnd9i0aRMBAQHY29sX2j8zM5OdO3fy8ccfF9r+6quvEhsbi0ajoUmTJubvx48fzwsvvIC7uzuKovDCCy9gMBj+zqVXOtINzYodO3aMx5/0p3q/Gdg8KgssVgbZl+PJ/G4ehw/ux83NTe1yRCm5Vzc0uQO2UsnJyfTo3Re77mMlfB/S5WXPY0q/UarHrG3oXqq/kShaHfrWT5L+3QK+XvelhG8lIQFshTIzM/Hv1QdaB2DXyk/tcso9u5Z+5KQUzH2bm3KmVI45618DS3XV6wXvLeKX9RG8O38BvXr1KrXjCusmAWxl8vPzefa5/+M3XR0c2g9Su5wKoZb//ee0NaXfJPPUboxx39OoYX3enz+HHj16lFF1BY4ci6Vz585MfnlSmZ5XqEueAVuZ16a/ySdhW6kx4B1Z062MKfkmsn7+iawDX9Cvd09CPliEo6Oj2mWJCkCGIpcDa9Z8yYefrsSh92sSvirQaHXYtfTjkf/7gO1xv9HCtQ3Hjh1TuyxRgUkAW4nDhw8zduJLVA98A519TbXLqdS0NnY4PDUOOoyk+1MBHDx4UO2SRAUlAWwFLly4QO/Aftj7T6RqnSZqlyP+x66VH7YBL/N032eIjo5WuxxRAUkAqywtLY0nezyNzusZbJv7ql2OuINtUx9sHv8XvZ/pT3p6utrliApGAlhFJpOJ/oOGcMuhCXY+dw8VFdbB3rUzxjqtmDxlmtqliApGAlhFL0+ZRvS537Dv/q8Sj/cXZcu26yjWf7ORvXv3ql2KqEAkgFXyyYoVrAr9Bvunp6HRSXdsa6er5kA1v+G89u8ZapciKhD5l6+Cffv2MX7CRBwDJpB3K4W8WykP3qkIub+fJ+fy6VKu7k+2zduhdyz+qhtZPx9F/0h9qtR2tlhNpUFb1Q79I/VLvJaevWtn4las5PTp07Ru3dpC1YnKRAJYBZeSU2je0hWSdkHS3z/OmbjjpVZTUTLidtHK3bPY7X+9o56S7FtWFEXh9u1UrqdnYuMegJ3X0+jsijfYQqOrQjX3p1i8dBkffxhi4UpFZSAj4USp+u233+jbty9Hjx4Fil6dwRqcOHGCd99fzLff7aBGvxnFvms3Xr9MztZ3uPZr0Uv8CFEUGQknykS9evWIjIwkLi4OAB8fH5UrKprBYGDNys95b/Ysbn/7DvnZxetipq/VkPT0NK5evWrhCkVlIAEsLMLNzY3Vq1cTExNDZmam2uXc07/+9U/8u3Uh4+TOYrXXaDTUcG4pQ5RFqZAAFhYzfPhwAIKDg1Wu5P6mTp5E3qmIYj8uMdV8jJMnT1q4KlEZSAALi6pduzZhYWFql3FfHTt2xJR5m/ys28Vqb9Lbknq7eG2FuB8JYGFRzs7W3SUNCh4r2DlUJz83q3jtdXqys3MsXJWoDCSAhUVpteXkr1hJRiLKqEVRSsrJvw4hhKh4JICFEEIlEsBCCKESCWAhhFCJBLAQQqhEAlgIIVQiASyEECqRABZCCJVIAAshhEokgIUQQiUSwEIIoRIJYCGEUIkEsKjUkpOTWbRoEZnpJZte8tixKFasWGGhqkRlIQEsKrWTcaeYMnUq1GmJvrpTsfaxa+nHwZ+O8e8Zs8jKKt4UlkIURQJYVGo9ewSw7bvv0F7/hczT+x7YPt+YQ/aR9bRu1ZLj0VHY2tqWQZWiopIAFpVez549OXLoAFVObCT9wCoUJb/IdnnpN0jb8G8eb1WXnw7t59FHHy3jSkVFIwEsBNCmTRtOREfRxJRMevh88nOzC32f+9s5UkNfZcKIZ9nwVajc+YpSIQEsxP84OTlxaN8eArybcfvr6eTdvgZAZuIRbm+awYqQxfxn1kw0siKGKCV6tQsQwprY2Niwbs0q5s5fwJx3X6WKSye0539kT8QO2rdvr3Z5ooKRO2Ah7qDRaHjj9df48rOP6Vgvn9hjRyV8hUXIHbAQ99CvXz/69eundhmiApM7YCGEUIkEsBBCqEQCWAghVCIBLCo0nU6Hl5cX7u7u9O3bl1u3bj3U8b799lsMBgOurq54eHjw7bffmr9LSEjAy8sLb29vfvnlF5YsWULr1q0ZNmzYw16GqKAkgEWFZmtrS2xsLHFxcdSqVYtly5b97WMdP36cqVOnsnnzZhISEtiyZQtTp07lxIkTQEE4P/PMM8TExNC8eXM+/LuKB6oAAAl0SURBVPBDtm3bxtq1a0vrckQFIwEsKo2OHTuSnJwMwN69e+nTp4/5uwkTJrBy5UoAmjRpwowZM/Dx8cHDw4OEhAQAFi5cyBtvvEHTpk0BaNq0KdOnT+fdd99l27ZtLF68mE8//ZTu3bszZswYzp07R2BgIIsWLSrbCxXlhgSwqBRMJhO7d+8mMDCwWO2dnJyIjo5m7NixLFy4EIBTp07Rtm3bQu3atWvHqVOnePrppxkzZgyTJ0/mhx9+YPny5TRo0IAffviByZMnl/r1iIpBAlhUaFlZWXh5eVG7dm1u3LiBv79/sfYbMGAAAG3btiUpKQkARVHuGoZc1DYhiksCWFRofzwDvnDhArm5ueZnwHq9nvz8P2c9y84uPPmOjY0NUPASLy8vDwA3NzeioqIKtYuOjqZNmzaWvARRgUkAi0rB0dGRJUuWsHDhQoxGI40bNyY+Pp6cnBxSU1PZvXv3A48xdepU5s6da74jTkpKYs6cOUyZMsXC1YuKSoYii0rD29sbT09PQkNDGT58OM8++ywGg4EWLVrg7e39wP29vLyYP38+ffv2xWg0UqVKFRYsWICXl1cZVC8qIo2iKMVu3K5dO+XOX8GEuB8fHx9iYmIoyd8zISoajUZzTFGUdndul0cQQgihEglgIYRQiQSwqNB27NhBq1atcHFxYd68eXd9n5CQQMeOHbGxsTH39/3DBx98gLu7O25ubixevNi8fdq0abi6umIwGOjfv795ePPOnTtp27YtHh4etG3blj179lj24kS5JwEsKiyTycT48ePZvn078fHxrF+/nvj4+EJtatWqxZIlS5g6dWqh7XFxcaxYsYLIyEiOHz9OeHg4iYmJAPj7+xMXF8eJEydo2bIlc+fOBQoGb2zdupWTJ0+yatUqhg8fXjYXKsotCWBRYUVGRuLi4kKzZs2oWrUqQ4cOZfPmzYXa1K1bF19fX6pUqVJo++nTp+nQoQN2dnbo9Xoef/xxNm3aBEBAQAB6fUEHog4dOnD58mWgoJdFgwYNgII+w9nZ2eTk5Fj6MkU5JgEsKqzk5GQaNWpk/uzs7GyeC+JB3N3d2b9/P9evXyczM5Nt27Zx6dKlu9p9/vnn9OrV667tGzZswNvb2zygQ4iiSD9gYTFBw18gJiYGgC7dn2Ln9nCqVatWZucvqutbcYcNt27dmtdeew1/f38cHBzw9PQ03/X+Yfbs2ej1+rummzx16hSvvfYaERERf794USnIHbCwCEVRcKxuj1anR1ulGlqtFpPJVKY1ODs7F7prvXz5svkRQXGMGjWK6Oho9u/fT61atWjRooX5u1WrVhEeHs7atWsLhfrly5fp378/q1evpnnz5qVzIaLCkjtgYREajYblH4bg7WUg6eJl/jPz7bvuIC3N19eXxMREzp8/T8OGDQkNDWXdunXF3v/333+nbt26XLx4kY0bN3L48GGgoGfF/Pnz2bdvH3Z2dub2t27donfv3sydO5dOnTqV+vWIikdGwokKbdu2bbz88suYTCaCg4N58803Wb58OQBjxozh119/pV27dty+fRutVouDgwPx8fHUqFGDLl26cP36dapUqcL777/Pk08+CYCLiws5OTnUrl0bKHgRt3z5cv773/8yd+7cQnfKERER1K1bt+wvXFiVe42EkwAWQggLk6HIQghhZSSAhRBCJRLAlVi3bt34/vvvC21bvHgx48aNK7J9UlIS7u7u9z1mUlJSoRddUVFRTJo0CYCVK1cyYcIEAJYvX87q1asfpnygYP22a9euPfRxhFCDBHAlFhQURGhoaKFtoaGhBAUF/e1j3hnA7dq1Y8mSJXe1GzNmDM8///zfPo8QFYEEcCU2aNAgwsPDzcNlk5KSuHLlCp07d2batGm4u7vj4eFBWFjYXfsmJSXRpUsXfHx88PHx4ccffwTg9ddf58CBA3h5ebFo0aK7Vh/+w8yZM1m4cCFXrlzBy8vL/Een03HhwgWuXr3KwIED8fX1xdfXl0OHDgFw/fp1AgIC8Pb25sUXX5R5hkW5JgFcidWuXZv27duzY8cOoODud8iQIWzcuJHY2FiOHz/Orl27mDZtGikpKYX2rVu3Ljt37iQ6OpqwsDDzY4Z58+bRpUsXYmNji7UacIMGDYiNjSU2NpZ//vOfDBw4kMaNG/PSSy8xefJkjh49yoYNGxg9ejQAs2bNonPnzsTExBAYGMjFixdL+aciRNmRgRiV3B+PIZ555hlCQ0P5/PPPWbNmDUFBQeh0OurVq8fjjz/O0aNHMRgM5v2MRiMTJkwgNjYWnU7H2bNnH6qOQ4cO8emnn3LgwAEAdu3aVWjmstu3b5OWlsb+/fvZuHEjAL179+aRRx55qPMKoSYJ4EquX79+vPLKK0RHR5OVlYWPj0+xXo4tWrSIevXqcfz4cfLz8x9qjoeUlBRGjRrFli1bcHBwACA/P5/Dhw9ja2t7V3tZBl5UFPIIopJzcHCgW7duBAcHm1++de3albCwMEwmE1evXmX//v20b9++0H6pqanUr18frVbLmjVrzPM8VK9enbS0tGKf32g08uyzzzJ//nxatmxp3h4QEEBISIj5c2xsrLm2tWvXArB9+3Zu3rz59y5cCCsgASwICgri+PHjDB06FID+/ftjMBjw9PTkiSeeYMGCBTz66KOF9hk3bhyrVq2iQ4cOnD17Fnt7ewAMBgN6vR5PT08WLVr0wHP/+OOPHD16lBkzZphfxF25coUlS5YQFRWFwWCgTZs25uHDM2bMYP/+/fj4+BAREcFjjz1Wyj8NIcqODEUWQggLk6HIQghhZSSAhRBCJRLAwqIetCqxEJWZBLCwmOKsSixEZSYBLCymOKsSC1GZSQALi3mYVYmFqAwkgIXFPMyqxEJUBhLAwmIedlViISo6CWBhMX9dlTg3N5fQ0FACAwPVLksIqyGT8QiL0ev1hISE0KNHD/OqxG5ubmqXJYTVkKHIQghhYTIUWQghrIwEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJICFEEIlEsBCCKESCWAhhFCJBLAQQqikREsSaTSaq8AFy5UjhBAVUmNFUercubFEASyEEKL0yCMIIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJICFEEIlEsBCCKESCWAhhFDJ/wNbasOS7XgMPQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "Stabilization.plot_sankey()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdeZyNdf/H8dd3ZswYM4OskYqixjbG1iKJGPsSWXOLUNGiFEolrbbcLVKJFlos3frJvichZZ0kkQhly24Ws39/f5xjGgYzw8xcc868n4/HPDjnXMv7Omfmcz7Xbqy1iIiIiIh4Ah+nA4iIiIiIZJaaVxERERHxGGpeRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpes5kxpqAxxhpjyjmdJTOMMTcZY046nUNExCmq2yKeJV80r8aY6DQ/KcaYM2ked89g3ObGmD+yMcuPxpj/ZNf0rpS19ndrbdGcmLYxZroxJt79Ph83xiwyxlRyvzbKGJN43mczIM247Y0xG4wxMcaYo8aYz4wxZS4wj+buL50B57+WQbYgY8zXxpi97vFvy2D4ksaYue48fxpjOp73egtjzO/u15el/RI87304+7MuK3nd06lrjIk0xsQaY9YZY6qd9/pN7vc4yhhzxBjzalbnIZJXqG5fnKfWbfdKwjhjzH73uLuNMaOzkG20MWarMSbZGPNsBsP6GGPeMsaccGd57bzXbzTGrHLX01+NMQ3SvNbPGJN03nJGG2OKZzarezoZfW8EG2MmGmOOGWNOGWOWZmX6+Vm+aF6ttcFnf4B9QJs0z33pdD4v96r7fb8OiAImpXltStrPxlo7DsD9xfQpMAYoDoQBvsD3xpjC502/J3Dc/W9WWGAlcB9wIhPDT3QPVwroA3ySpqCXAWYAg4ESwDbgi/PGf/W8Zb3lQjNxf0mma6SNMYHAbHeOq4D/AbOMMX5pXl8GzHNnvA74KhPLJZInqW47Kqfq9nCgMlALCAGaAFuykGsH8DSQmSbvcSACqOKeXxdjTK80r88EVgHFgNeAb4wxaVcIvjtvOYOttcfOn4m7ob9YI33R7w23yUAAcJM7xyUbcvlXvmheM2KMCTTGvGeMOWiM+dsY84YxpoB7LWsWcEPaNS9jzB3GmJ/ca0oH3Gt3fpcx3+bGmD+MMc+5t5TtT7tFwbi2Do4zxvzlntfKNM3KvcaYbcaYk8a1pa9SmvEOGWOecq9NRhtjPjDGlDHGLDXGnHavSRd2DxtqjElKM+6Pxpjh7n9PG2MWGGOuSvN6X2PMPnfeIe551c9oWa210cB0oNqlhnMv31jgRWvtV9baOGvtAeB+9yCPpRm2CNAO6AfUMOdticwgT6y1dpy1dg2QkkGmq4A2wDBrbYy19ltgMXD2s+oEbLDWzrbWngFeBOoZY8pnNk8mRABx1tr3rbXxwH9xFf+z7/2DwA5r7Xhr7Rn3zy8XWZ5zPnP3c6lblowxfu7fu2PGmF3GmAHnDy/iNNVtj6zbdYGvrbWHrcvurKyIWGs/sdYuBqIzMXhPYIy19qC1dh/wNtDLnTcMV8P4qjvrNGAXcE9ms2Qko+8Nd4bGwCPW2mPW2mRr7cZLTO+cz8y4muaP0jy+rM/YU6l5dXkZ11pidaA20BAY4l7Lag/sPm/NKxHXH2Mx4E5cv6B9L3Pe1wMGKOue5gRjTLD7tXFAKK4/+GLAC4B1N2mTgUdwrdGtBOacV4jbA3fhWuvsimur3VNAaSAY6H+JTPfh+gMrAxQFngAwxoQDbwKdgXLunxKZWUh30e0GbM5g0GrA1bi2LKay1iYD/4eriTurM3AE1xr0Cv4tlGfnucMY0yEz+TIQCkRZa/emee5noKr7/1Xdj89mPYlrS1FVss/580gBtqaZx23AX+4vuqPGmOXGmMqXOa/HcP3uVANuATpeenARR6hun8sT6vaPwDPGtVs+XX00xnxsjHkzM9kyoQppaibpa/bv7o0NF3o9O2T0vXEbroZ5tLtm/2yMaXs5M7qSz9hTqXl16Q4Mt9YetdYexrULocfFBrbWrrPWrnevKe0CPsJVcC5HLDDSWptorZ2Fa3d2RWNMAVzN2OPW2kPuea1yF4NuwCxr7XfW2gRgBK5f1Dpppvu2e3n2AT8Aa6y1v7j/WGcDNS+RaZK1dpe1NgZXYxjufr4zrrXmH91b/14g49+h543rxIIduHYhpf2y6OHeAnH2p5h7OSxw+ALTOsi5f5A9gWnWWgtMBf5jjPE9+6K19mZr7f9lkC8zgoFT5z13CteWz8y8Du73Ic3Ph9mcoRyu3+NRuL5QV+DaDeZL1nUG3nRvsTiGazegSF6jun0uT6jbL+PaAtoT2OTeYt7t7IDW2j7W2qcyyJYh9+cQwLk1M6s1+67zlvPXLMbITM2uDezHVbMHAdOMMRWzOB+4vM/Yo3n1wmWGMcbgWmNMu3a0F7jmEuNUMcYsNMYcNsacxrWb+HLXco64t6KdFYvrl74M4AfsvsA4ZdPmdRfG/edlTltEzlzgcTAXd+gCec7O96808z1N+j/O871urS1qrS1jrW1/3lro5+7Xzv4cB47i2qJR+gLTKuN+HWPMjcAdwNldTl/jOhY04gLjXalo4PxjbQvjOhYsM6/Dv+/D2Z+HAYwxAWkLJK4vsqVpnnsyk/M4A3xrrV3u/mIciet4tcsphOd8zuf9X8RxqtsXlOfrtrvZf8daezuuev0m8Jm7nmcba20iEM+5NTOrNXvlecuZulXWvYfrbM1+EhiepmbPzOQ8zuD6nEZbaxOstUtxrbA0voxFvpzP2KPl++bVvdXuEK7dQGddh6uogGtt8nyTgE3AjdbawsAruP5ws9NBIAm44QKvHSBNXvfWtWv4N3NOOYhrbfHsfAsDRbJ5HltxFexOaZ90L2N7YLn7qbMnaC01xhwCfsf1pXHOoQPZZDtQ2BhzXZrnagBn18R/dT8+m7UIrs8nwzV1a2182gIJbAAi0jz39kXm4YNrV93ZeWzh3N9Vy6V/J33cWyfOSnuiwjmfM3BtRsshkptUt7OcKa/U7VTWdd7Bm7iazNBszgOuE2drpHl8fs2+yRhT8CKvX5K1NiJNzX4beDlNzT57mFVG3xtna/aFflcvJm3ei9bsHPqM85R837y6TcO15lTcGFMKeJ5/zxY/DJRKczwTuDb7n7LWRruP23kwuwO51xw/A94xxpQ2xvgaY+q7i8EMoL0xpoG7AXkWOIar8clJXwH3Gtclm/xxFf9LnuyUVdbaJOAZ4FVjTCf3lsmywBRczel491aXHsBzuHaNnf25D7jHpL8iwQW5p322GPifV8jSZjqB6yz+V4wxhYwxDYHm/LvVdyZQ1xjTxj2Nl4EfrLV7srj4l7IUCHQfKxYADARigNXu1z8DGhpj7nL/jjwD/Alc7HJBBuhtXJeTaQ3cDIS4j7/7ChhojLnauE5+GZSNyyGSXVS3MydP1G0AY8zTxpg7jeuSWQWMMQ/hOizh54tOPA33OAVx9S5+7ulcrI/5DBjsrmPX4tpCOtmddwuwExjmztoZ116q2Zf3DqSXie+NZbiulDPI/XvSCNdxsMsuMdn7jTH+7kML2uGq2QXIhc84r1Hz6vIirrW0X4FIYA3/Huf3MzAH2Jvm+J6BQF9jTDTwHq6ilBMG4DqgezOuIvcqYNx/eH2AD3GdsNQYaOcuIDnGWrsZ1+WgZuHaWnAQ166J+GyezxRcy/csrsuMnD1rvr51nQx1F64THj5wH1d2yFp7CFcTuR/X8T8Y15ny915iVntx7bopjuvkiTPGmKvd475sjJmVZtgHce1iPIqrAPax1u505z0AdMG1C+w4ri2i518Tcpg593qBf2fxPTnDv1dWOInrZI57zn7m1nVlgd7ubCdwXYKmvXvX5IUk4zrZ5QjwKK7jBV8HKuD6ovkB19/EelwFOFs/Y5FsoLqdCXmobuOe5zhcKxf/AA/gqmN/AxhjJhtj3ubiPsdVs9vjel/P8G+9b2KMOZpm2HG4tvj+huv343/W2slpXu8ENMBVT1/CVS/TXjaxoUl/ndfqmX0/3C71vRGP66TBTsBp4F2gq/t47IuJx/W9Ndu9/BFA59z6jPMS49r7IpJ1xnUpkONAWWvtQafzSOYYY0KBrdbaTF0myBjTHhhlrb05Z5OJSE5T3fZM7sPjOlprV2diWK//jLXlVbLEGNPWuK6vGIxrS+NP3vrHkV8ZY0KMMU3du7Kuw3Xm6qyMxhORvEl12/vlt89YzatkVSdcJ0r8jetkg0veplE8kg+uS26dwnXYwCZchxWIiGdS3fZ++eoz1mEDIiIiIuIxtOVVRERERDyGmlcRERER8RiZOtv4rBIlStjy5cvnUBQRkZyzcePGo9bakk7nyE2q2SLiqS5Vs7PUvJYvX54NG3L6esoiItnPGLM346G8i2q2iHiqS9VsHTYgIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1rx4sODg4x6Zdvnx5jh49munhIyMjWbBgQY7lERHxdKrZItlDzatkCxVCERHPoZotnkzNq5c5cuQI9957L3Xr1qVu3bqsWbMGgHXr1lGvXj1q1qxJvXr12LFjBwDJyckMGjSI6tWrExYWxrvvvps6rXfffZdatWpRvXp1tm/fftHpJCQk8OKLLzJjxgzCw8OZMWMGMTEx9O7dm7p161KzZk1mz56d+2+GiEgep5otchmstZn+qV27tpW8IygoKN1z3bp1s6tWrbLWWrt3714bGhpqrbX21KlTNjEx0Vpr7dKlS22HDh2stda+//77tkOHDqmvHTt2zFpr7fXXX2/HjRtnrbX2vffes3369LnkdD799FP76KOPpuYYOnSo/fzzz6211p44ccJWqlTJRkdHZ+PSi2QNsMFmod55w49qdt6imi2SeZeq2X5ON8+SvZYtW8a2bdtSH58+fZqoqChOnTpFz5492blzJ8YYEhMTU4fv168ffn6uX4VixYqljtuhQwcAateuzf/93/8BXHQ651uyZAlz5sxh7NixAMTFxbFv3z4qV66c/QstIuKhVLNFsk7Nq5dJSUlh7dq1BAYGnvP8448/TqNGjZg1axZ79uyhYcOGgGvLuzHmgtMKCAgAwNfXl6SkJACGDRt2wemcz1rL119/zc0335w9CyYi4oVUs0WyTse8epmmTZsyfvz41MeRkZGAa+37mmuuAWDy5MnnDD9hwoTUQnf8+PFLTv9i0wkJCSEqKir1cbNmzXj33XdxbfmHzZs3X/5CiYh4KdVskaxT8+rBYmNjKVeuXOrPm2++ybhx49iwYQNhYWFUqVKFCRMmADBkyBCGDh3KHXfcQXJycuo0+vbty3XXXUdYWBg1atRg6tSpl5znxabTqFEjtm3blnrw/7Bhw0hMTCQsLIxq1aoxbNiwnHkTREQ8hGq2SPYwZ9eyMqNOnTp2w4YNORhHRCRnGGM2WmvrOJ0jN6lmi4inulTN1pZXEREREfEYal5FRERExGOoefVCixYt4uabb6ZixYqMGjUq3evWWgYMGEDFihUJCwtj06ZNGY57/PhxIiIiqFSpEhEREZw4cQKAY8eO0ahRI4KDg3nsscdyfuFERLxY+fLlqV69OuHh4dSp49pj2qtXLypUqEB4eDjh4eHUq1fvguP++eef3HrrrVSqVIkuXbqQkJAAwHfffUeRIkVSx3/llVdybXlEcoKaVy+TnJzMo48+ysKFC9m2bRvTpk075xqCAAsXLmTnzp3s3LmTiRMn0r9//wzHHTVqFI0bN2bnzp00btw4tbEtWLAgr776auq1AUVE5MqsWLGCyMhI0h6v/MYbbxAZGUlkZCQ//PDDBcd75plnGDhwIDt37uSqq67i448/Tn3tzjvvTB3/xRdfzPFlEMlJal69zLp166hYsSI33HAD/v7+dO3aNd1t/mbPns3999+PMYbbbruNkydPcvDgwUuOO3v2bHr27AlAz549+eabbwAICgqifv36FCxYMHcXVEREUllr+fbbb+nYsSNwbp0W8TZqXr3M/v37ufbaa1MflytXjv3792dqmEuNe/jwYcqUKQNAmTJl+Oeff3JyMURE8iVjDE2bNqV27dpMnDgx9fnBgwen7vbv3r17uvGOHTtG0aJFU++8dX7tX7t2LTVq1KBFixb8+uuvOb8gIjlId9jyMhe69Nn5d2O52DCZGVdERHLOmjVrKFu2LP/88w8RERGEhoYCrsMGzm5VvZBL1e9atWqxd+9egoODWbBgAffccw87d+7MmQUQyQXa8uplypUrx19//ZX6+O+//6Zs2bKZGuZS45YuXZqDBw8CcPDgQUqVKpWTiyEiki+drbmlSpWiffv2rFu37qLDNmvWjPDwcPr27UuJEiU4efJk6p230tbvwoULExwcDEDLli1JTEzk6NGjObwkIjlHzauXqVu3Ljt37uTPP/8kISGB6dOn07Zt23OGadu2LZ999hnWWn788UeKFClCmTJlLjlu27ZtmTJlCgBTpkyhXbt2ub5sIiLeLCYmJvWWrTExMSxZsoRq1apddPjFixcTGRnJRx99hDGGRo0aMXPmTODcOn3o0KHULbPr1q0jJSWF4sWL5/DSiOQcHTbgZfz8/Bg/fjzNmjUjOTmZ3r17U7Vq1dRbDvbr14+WLVuyYMECKlasSKFChfj0008vOS7As88+S+fOnfn444+57rrr+N///pc6z/Lly3P69GkSEhL45ptvWLJkCVWqVMn9hRcR8WCHDx+mffv2ACQlJXHffffRvHlzpk+fzuDBg3nttddSh123bh3+/v7njD969Gi6du3KCy+8QM2aNenTpw8AM2fO5IMPPsDPz4/AwECmT5+uQ8LEo+n2sCKSL+j2sCIinkO3hxURERERr6DmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVVIlJCQwadIkjh496nQUERHJhLlz57J161anY4jkKjWvAsCxY8eo37AxT704khq16/Lbb785HUlERC4iOTmZRwc8yX29+3F7/QbMmTPH6UgiuUbNqxAVFUVYrTps2rKVpIAiHDkDtW+5VWvzIiJ5VPtOXfngvfEkBZcmsci1tLvnHiZ/9oXTsURyhZ/TAcR5wcHB9HuwD3FxcYx4/XU6d+5CxUoVKVOmjNPRRETkAjre05YqN1di9KiRADzz7FCqVg51OJVI7jDW2kwPXKdOHbthw4YcjCNOM8awceNGatWq5XQUkWxljNlora3jdI7cpJrt/YwxtG3bltmzZzsdRSRbXapm67ABEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeZV0AgICnI4gIiKZ5O/v73QEkVyl5lVSJSYmAlC+fHlng4iISKZVqFDB6QgiuUrNq6R66qmnAAgKCnI4iYiIZGTnzp0A9O7d2+EkIrlLzaukGj9+PHfccYfTMUREJBMaNGgAQGhoqMNJRHKXmlcB4OuvvwZg9uzZDicREZGMHDt2jEOHDvHmm286HUUk16l5FSIjI+nYsSP169enePHiTsfxOrGxsezdu9fpGCLiJc6cOUOJEiUAePLJJx1O432Sk5NTD8mQvEnNaz737rvvUrNmTQoXLsyqVaucjuOVXn7lVWrUqqMGVkSu2Pr16ylUqBAAUVFRGGMcTuR9Fi5cSGjlyixdutTpKHIRfk4HyE7JyclYa3Ns+vHx8Zw4cSJTw54tKBn9m9Fr5/P398fH5/LXOU6cOMH27dsZNmwYP/30EwCNGjVi3rx5xMbGXvZ05eISExOJxZ/GzVqyef2PhISEOB1JJE9ITk4mJSUlx6Zvrc10zU7r/Bqc1cdpFShQgAIFCmQ5w1nx8fH8/vvvTJ06lXfeeSf1+RMnTuDn50dcXNxlT1suLCEhgYDCJejYpRvr1q7h5ptvdjqSnMdrmteUlBSuue56/jl4AGNyZoNySkpyjkw3q/wLBl72uAlxZ8557OPjy+of1lK0mA4XyDEWijR/gpP7t9C+YxcWL5iLr6+v06lEHNexy3188/VX+PjmzFdRSnJSjkw3q/z8L//a2UkJ8emnV8CfkqVKX0kkyUCxOq1JKVqOxs1a8PPG9TqkLo/xmubVx8eHR/v3578TP6dwp9fxKVDQ6Ugi57A33c7GWS/x5NODefdtnWQh8sRj/Vm6bBlFOo2gQPFyTscRSSf65H5atm3PqhXLdDOIPMSrjnl94fnnaHpHbWKWjMPanNsVJXI5jK8fwa2eYcq0/zFx0iSn44g4rmHDhvx3zEii5r5O8pkop+OIpBNU/37+OJFMn4f65ehhiZI1XtW8GmP4fPInXFcwnpgfpjkdRyQd38AQQtq+wFODn2XlypVOxxFx3MMPPcT9Xe4lZsEYbB7ZzS9ylvHxJajZk8xdtoqx/9Ues7zCq5pXgICAABbPn4Pfn6uJ2bbC6Tgi6RQoXo6g5k/RrkNH/vjjD6fjiDjunTfHUuuGq4leMVFbtyTP8QkoRHCb53h5xCjmzZvndBzBC5tXgFKlSrFs0QLiVn1K/P7fnI4jkk5g+XAK1O1C42YtOXnypNNxRBzl6+vLrJkzKBr1J7Gb5jgdRyQdvyKlCGn1LN169OSXX35xOk6+55XNK0C1atWY/sVnRM0fQ9Kpf5yOI5JOUHgLYkpUofU995KUpN2lkr+FhISwfPECkjZ/w5ld652OI5JOwDWhFLyzNxEtWvHPP+ornOS1zStAq1ateOn5Z4ma8xop8bqGqeQ9QXf1ZtvBKPo/+rjTUUQcV758eebP+YaYpe+ScGSP03FE0gmq0pCkCvVp2rKNrrHrIK9uXgGefmog9zRtSPTiN7F55DqtImcZH1+CWw7iq7mLeWfceKfjiDiuXr16fPDuO0TNGUFyjA6pkbwnqF43/oovyH96PqBjtB3i9c2rMYaPJk7g5mL+xKz6zOk4Iun4BAQR3PZ5nh/+MosXL3Y6jojjevT4D/379CR6/mhsUqLTcUTOYYwPQU0HsOzHSF557XWn4+RLXt+8guv2fPPnzCLw4GZitixxOo5IOgWKXk1wy0F06nofv/2mkwxFRo14jTvCKhK97D1t3ZI8x6dAQYJbD+WNt8fzv//NdDpOvpMvmleAYsWKsXzJQhJ+/JK4fVucjiOSTsFrqxFwR0+aNGvJ0aNHnY4j4igfHx++mvoFpVOOErNOzYHkPX4hxQlpM5QHHnyIjRs3Oh0nX8k3zSvATTfdxNdfTSdqwVgSTxxwOo5IOkHVGhN/bV1atLmHhIQEp+OIOKpQoUIsXTgPfltC7I4fnI4jkk7A1RUp1Kg/zVq1Yf/+/U7HyTfyVfMK0KRJE8aMeJXoOa+REhftdByRdILq92D3acsDfR/S7lLJ96655hoWz59L7IoPSDi8y+k4IukUurkeVG5KRIvWxMbqyka5Id81rwCPPvII3dq3JXrBWF2BQPIcY3wIajaQ+d+uYdSYN5yOI+K42rVr88nED4maO4KkqGNOxxFJJ+iWjhz2KUHn+/5DSkqK03G8Xr5sXgHGj3uLsGuvIua7j5yOIpKOj38gwW2f5/XRY5k9e7bTcUQc17lzJ54e8CjR80aSkqjra0reYowhuMmjrPl5J0OfH+Z0HK+Xb5tXPz8/Zv/f/wg58Tsxm3SvYsl7/AqXJKT1s/ynV29+/vlnp+OIOG74sBdocls4MUvGYa22bkneYvwKENz6Wd7/aDKff/6F03G8Wr5tXgGKFCnCt0sWkrTpa87s1pmCkvcElL2Zgg36EtGiFYcOHXI6joijjDF8MeVTrg2II2btdKfjiKTjG1SUkLbP0f/xJ1i7dq3TcbxWvm5eASpUqMDcb/6PmCXvkHj0L6fjiKQTVLkBKRXv0u0IRYCCBQuyeP4cfHd9T8y2lU7HEUnHv2R5giIep1Xbe9i7d6/TcbxSvm9eAerXr8+7b/2XqLmvkxx7yuk4IukE3d6N/UnB3NejZ5avQODr60t4eDjVqlWjTZs2nDyZs7fc3LBhAwMGDMhwuHr16uVoDvFepUuXZtmiBcSt+pj4/dudjiOSTuCNdfENb8fdTVsQFRWVpXFVszOm5tXtgQd60bdHV2Lmj8Em63aEkrcYYwiKeJwV639h+MuvZGncwMBAIiMj2bp1K8WKFeO9997LoZQuderUYdy4cRkO98MPum6nXL7q1asz9bPJRM0fRdLpf5yOI5JOoVptORlSgfYdu5CcnPkrG6lmZ0zNaxpjx4zmltDriFk+QdfXlDzHp0AAQa2H8tb4CUyfPuOypnH77befcyHtN954g7p16xIWFsbw4cMB2LNnD6GhofTt25dq1arRvXt3li1bxh133EGlSpVYt24dAOvWraNevXrUrFmTevXqsWPHDgC+++47WrduDcBLL71E7969adiwITfccMM5BTI4ODh1+IYNG9KxY0dCQ0Pp3r176t/fggULCA0NpX79+gwYMCB1uiIAbdq04cWhzxA9ZwQp8bq+puQtxhiCGz3Exl2HePLpwZc1DdXsC1PzmoaPjw8zZ0ylePwBYjfMcjqOSDp+wcUIafMcffv1Z/369VkaNzk5meXLl9O2bVsAlixZws6dO1m3bh2RkZFs3LiR77//HoA//viDJ554gi1btrB9+3amTp3K6tWrGTt2LCNGjAAgNDSU77//ns2bN/PKK6/w3HPPXXC+27dvZ/Hixaxbt46XX36ZxMT0ezY2b97M22+/zbZt29i9ezdr1qwhLi6Ohx9+mIULF7J69WqOHDmSpeWV/GHwoKdp0+ROYha/ret2S55jfP0IbjWEKdNnMnHSpDS/DXkAACAASURBVCyNq5p9cWpezxMcHMyyRfNJ+WU+sTt/cjqOSDr+pW+g0N2P0Lx1W/7+++8Mhz9z5gzh4eEUL16c48ePExERAbgK4ZIlS6hZsya1atVi+/bt7Ny5E3CdyFi9enV8fHyoWrUqjRs3xhhD9erV2bNnDwCnTp2iU6dOVKtWjYEDB/Lrr79ecP6tWrUiICCAEiVKUKpUKQ4fPpxumFtuuYVy5crh4+NDeHg4e/bsYfv27dxwww1UqFABgG7dul3O2yVezhjDxxMncGNRH2JWf+50HJF0fANDCGn7Ak8Nfpbvvvsuw+FVszOm5vUCrr32WhbOm0Ps8vdI+Ge303FE0il00+2Yqi1o0rwVMTExlxz27PFTe/fuJSEhIfX4KWstQ4cOJTIyksjISP744w/69OkDQEBAQOr4Pj4+qY99fHxISkoCYNiwYTRq1IitW7cyd+7ci14JIe20fH19U8fPaBgduiOZ5e/vz8K531Bw/0ZiflnqdByRdAoUu4ag5k9xz72d+OOPPy45rGp2xtS8XsQtt9zCpA/eI2rOCJKjTzgdRySdoLodOOpfmk5du2fqdoRFihRh3LhxjB07lsTERJo1a8Ynn3xCdHQ0APv37+effzJ/4supU6e45pprAJg8efJlLcOlhIaGsnv37tStBjNmXN5xvpI/FC9enGWLFxC/9gvi9v3idByRdALLh1Pglq40btoiU1cQUM2+ODWvl9CtW1eeeOQh9+0I452OI3IOYwxBdz/C2m1/MuTZCx+7dL6aNWtSo0YNpk+fTtOmTbnvvvu4/fbbqV69Oh07dszSJV2GDBnC0KFDueOOO7J0Jm1mBQYG8v7779O8eXPq169P6dKlKVKkSLbPR7xHaGgoX8+YRvTCsSSeOOh0HJF0gmo0J6ZkVVrfc+8Ft2ieTzX7wkxWNvPWqVPHbtiwIUeC5FXWWjp06sL3O48S3PwpjDFORxI5R3LsKU5NH8K4Ma/Tq1dPp+Nkq+joaIKDg7HW8uijj1KpUiUGDhx4WdMyxmy01tbJ5oh5Wn6s2QDj33uP5197gyJdx+ATEOR0HJFz2JRkoua8TocGNflo4gSn42Sr3KrZ2vKaAWMMUz+fwjW+UcT8pN2Wkvf4FipCcNvneezJp1izZo3TcbLVpEmTCA8Pp2rVqpw6dYqHH37Y6UjiAR579FG63NOS6AVv6AoEkucYH1+CWzzN/+Yt4e133nU6TrbKrZqtLa+ZdOjQIcJq1cHW/Q9Ble90Oo5IOmd2byTh2/Fs3rCO8uXLOx0nz9GW1/wlKSmJu5u2YFt0IYLvfsjpOCLpJJ48xKkZzzJz2uc0b97c6Th5jra8ZoOrr76aJQvmcWblROIP/u50HJF0Am+ojW/NDtzdtAWnT592Oo6Io/z8/Jg7ayYhx7cTvXm+03FE0ilQ9GqCWw2mc7fu/Pbbb07H8Si5vuU1ODg49Uy53PTiiy/SoEEDmjRpckXT+eabb+hw7734Fi2L8fXDN+gqSrR6Ct/gq/j7g974+AeCj2udoFjEIxQsV5mEI3s5vuxDkqOOApagqndTpF7Xc46f/efrV0mOPUmZHv+9onySv1lrifl2AmFFk1i6aD6+vr45Nq969erlqdsFZkRbXi+fJ9ftXbt2UfuW24iKS8K3UBEw5orrdnLMCY4tHEfS6SOQkoxvkVKU7vRyNi215DcxW5fjv2UWP29aT4kSJXJsPt5Us/NN85qdrrrqKpILFadol9GcWjsDmxhHsSYP8/cHvSnT8y1XgXRLSYzn4CePUqzpIwRWqEVKYhxHZo2kUMW6hNRy3TYtJS6aA588hilQkFKdXqJA0audWjTxAjY5iajZr9C1aX3eH/+O03HyDDWvl8/T6/aqVau4666GlOo+moLXVObEyilXVLePLRpPgRLXUrhOOwAS/vkT/1IVnFo88QLRq6ZQIeUAa1Z+i7+/v9Nx8oQ8edjApe6Nu379eurVq0eNGjW45ZZbiIqKIi4ujgceeIDq1atTs2ZNVqxYAbiuVXbPPffQpk0bKlSowPjx43nzzTepWbMmt912G8ePHwegV69ezJw5E4Dy5cszfPhwatWqRfXq1dm+fTsAR44cISIiglq1avHwww9z/fXXc/To0XTZixQpQkS92sQseZuAclVIPHHgossZs20lAddUJrBCLQB8ChSkWEQ/Tv04899hdvxAYMVbCKrcgNjfvs+Gd1fyM+PrR1DLIXwx8xs+mJBzZ7Jm5j7X4l08tW7feeedFCt2FdGL3iT5zGkKXlvtiup2csxxfEP+3UKmxlWuVFD9Huw+Db36PJhj9dObarajx7xe6N64CQkJdOnShXfeeYeff/6ZZcuWERgYmHqHiV9++YVp06bRs2fP1LtDbN26lalTp7Ju3Tqef/55ChUqxObNm7n99tv57LPPLjjvEiVKsGnTJvr378/YsWMBePnll7n77rvZtGkT7du3Z9++fRfNPu7tN6kQnMLp1V/iX7J86vOHpz3HgU8f5+BnTwGQeHQv/ldXPGfcAleVwSacISU+FoDY31YSVLkBQVXuIkbNq2QD34LBhLR9nsHPPp/aMOSkC/0ti3fy1LodHBzMfR3vIXr+aGJ3/nhFdTukZiuOLRzHoWlDOfXDDJKijl3ReypijA9BzZ5kwXdrGTl6TI7Pz9NrtqPN64Xujbtjxw7KlClD3bp1AShcuDB+fn6sXr2aHj16AK4LUV9//fX8/rvrxKlGjRoREhJCyZIlKVKkCG3atAE4556+5+vQoQMAtWvXTh1m9erVdO3aFYDmzZtz1VVXXTR78+bNiYs+hc/pg/gVKZX6fOluIyj7wLuUuf/NNENf/NqwyTEnSDxxkIByVSlQ7BqMjy8JRy6cWSQrChS7huAWg7jn3k6pfys55UJ/y+KdPLlur129iqR//iRh388UvrVj6vNZrduBN9Tmmoc/IqRGMxKP/83ByU+QHHvq4m+aSCb4+AcS3OZ5RrzxJrNnz87ReXl6zXa0eb3YvXEvdCOAS23Szsw9fS82Ttr79mZls/mKFSvYunUrG9avI3HdDOL+/vWCwxUocR0Jh3ae81ziyUMY/0B8AgoR89sqUuKi2T+hD39/0JukU4eJ+W1VpnOIXErB68MocGs3mjRryYkTOXeb48zcC1u8g6fX7YMH9nN96WLE/br8osNlVLcBfANDCKrSkBKtnyagTCXi/tqa6RwiF+NXuAQhrZ/lP7168/PPP+fYfDy9Zue5S2WFhoZy4MAB1q9fD0BUVBRJSUk0aNCAL7/8EoDff/+dffv2cfPNN2frvOvXr89XX30FwJIlSzL1ZV+lShVmTP2C6PlvwAUuhh1UpSHxf2/jzJ5IwHUiwIllH1L41nsBiPltJaU7v0y5/p9Qrv8nXN3zbR33KtkquEZzYkuH0bJtexITE52OI17Ik+p2cHAwyxbNJ2XLXGL/+OmCw2RUt8/s/ZmURNfhDynxsSSeOIhf4ZLZtUiSzwWUuYmCDR4kokUrDh065HScPCnPNa/+/v7MmDGDxx9/nBo1ahAREUFcXByPPPIIycnJVK9enS5dujB58uRz1hyyw/Dhw1myZAm1atVi4cKFlClThpCQkAzHa9GiBa8OH0bKmVOkJJw55zWfAgGUvHcYp36Yzv5JD3Pwk8fwv7oSIbVak3TqMMmnj+BfNjR1+AJFr8YnoBDxB3Zk67JJ/hbUoBe/H43noX6PeNyB+ZL3eVrdvu6661gwdzaxS8dDSvotTpeq2wAJh/7g0JSBHPjkMQ59MYiQGk0JKHNTti6X5G9Ble8kpVJDmrZsw5kzZzIeIZ/RHbbSiI+Px9fXFz8/P9auXUv//v2JjIzM1LjWWvo+1I9ZqyIJafs8xifnrq8pcjlS4mM59dWzDB/0OE8NfNLpOLlOl8ryTldSt6dOncbDTzxN0a5j8A26+LGyIk6w1hK96L/cdVMpvv5q+gUPzfFmefJSWXnRvn37qFu3LjVq1GDAgAFMmjQp0+MaY5jw/ngqlw4idvWFz5QVcZJPQCFC2j7Pi6+OYNmyZU7HEckWV1K377uvG4891IfouaOwFzjsS8RJxhiCmjzOivVbeeX1EU7HyVP8nA6Ql1SqVInNmzdf9vgFChQgrHo1fl+ZubV+kdzmUzAEn4Ag9h+4+DUuRTzJldbt2jVrkPjBh5CSAtpjJnmM8fPHp0hpdu3e43SUPEVbXrPRhxMn8uXMbyh49yNORxFJx6YkE73ov9zbKoL73ZcvEsnPNm3aRK++D7kO9fIr4HQckXRi1s3kak4yQXdLPEe+bl579+5NqVKlqFat2hVPa8WKFTw9ZCghbV/ANzDjk7xEclvM95OpXDKQiR+8d1nHTr3++utUrVqVsLAwwsPD+emnC5+pDdC3b1+2bdsGuO6MdKE71YlkVXbW7AMHDtC0ZWsCGz5MwHk3JBDJC2J2rMFsX8rShfMoVKhQlsf35pqdrw8b6NWrF4899hj333//FU1n586d3HNvJ4JaPE2BYtdkUzqR7BPz8yIKHd7CvE3rKVAg61uY1q5dy7x589i0aRMBAQEcPXqUhISEiw7/0UcfXUlckQvKrpodGxtLRIvWEBpBUGj9bEonkn3iD/3BmRUT+P7bZZQtWzbL43t7zc7XW14bNGhAsWLFrmgaJ06coEmzlvjf2o3A62tkUzKR7HNm788krpvO8iULL3n3oUs5ePAgJUqUSL3MUYkSJShbtizLly+nZs2aVK9end69exMfHw9Aw4YN8faz3CX3ZUfNTklJoUv3HhwyxQi6tVM2JRPJPklRR4meO5IpH0+iVq1alzUNb6/Z+bp5vVKJiYm0vudeYkpVJ6hGc6fjiKSTeHw/MQv/y6yZX1GpUqXLnk7Tpk3566+/uOmmm3jkkUdYuXIlcXFx9OrVixkzZvDLL7+QlJTEBx98kI3pRbLfcy+8yOrIHQQ3eTTfXXpI8r6UxDii545kyMDHuffeey97Ot5es9W8XiZrLf0eeYzfDscQ1KCX03FE0kmOiyZqzuuMHT2CRo0aXdG0goOD2bhxIxMnTqRkyZJ06dKFDz/8kAoVKnDTTa6Ls/fs2ZPvv9fd4STv+uKLL3lv0icEt34W4+fvdByRc1ibQszid2hWvw4vPP/cFU3L22t2vj7m9Uq8/c67zFywlCKdR+uGBJLn2OQkYuaPpkfn9vR7+OFsmaavry8NGzakYcOGVK9enSlTpmTLdEVyw9q1a3n40ccp0vFV3ZBA8qSYH6ZxfaFEPvv042zZK+DNNVtbXi/DwoULGfbyqwS3eR6fgKyfASiSk6y1xHw3kZoVSjPurf9myzR37NjBzp07Ux9HRkZSunRp9uzZwx9//AHA559/zl133ZUt8xPJTnv37qVV23sIingc/5LlnY4jkk7MthX4/bmaxfPnZMstlL29Zufr5rVbt27cfvvt7Nixg3LlyvHxxx9nOM62bdvoct9/CG41hAJFr86FlCJZE7NpLkVO/8msmTPw9c2evQLR0dH07NmTKlWqEBYWxrZt2xg1ahSffvopnTp1onr16vj4+NCvX79smZ/IhVxOzY6KiqJJ81b41GhLoYq35EJKkayJ3/8bcas+ZdmiBZQsWTJbpuntNdtYazM9cH64T/alHD16lLBadUgM60BQtcZOxxFJ58yuDSSseI/NG9ZRvnx5p+PkKZe6T7a3yu81Ozk5meat2rLpSIpO0JI8KenUP5yc8QwzPv+UVq1aOR0nT7lUzc7XW16zIj4+nuat25Fw7a1qXCVPSjiyl5il45g3e5YaVxFg4KAhbPhjP8F3P6zGVfKclPhYoua8xkvPP6vGNYvUvGaCtZZefR7kz2gfgur/x+k4Iukkx54iau7rvD/ube644w6n44g4buKkSUye+hXBrZ7B+OrWr5K3nL1d9z3NGvH0UwOdjuNx1LxmwohRo1n4/U8ENX0SY/SWSd5ikxKJnjeKfg/cz/3398jVeWfn7TpFssvKlSt5avCzhLR9Ht/Awk7HEUknZtVn3FyiIB99+EGu7xXwhrqtTiwDs2bNYuTYtwhu/Rw+/gWdjiNyDmst0cvfo161Gxk98vVcn3+vXr1YtGhRrs9X5GJ27dpFuw4dCWr+FAWKX+t0HJF0YrYsIfDgZubP/r/Lul33lfKGuq3m9RIiIyPp8UAfQlo/i1/hEk7HEUknZv3/USrpCF9N+wIfn9z/c86O23WKZJeTJ0/SuGkLCtTtTGD5cKfjiKQTt28LCT9+yfIlCx2rnd5Qt3WTgos4dOgQTVu2JvCuhwgoc5PTcUTSif19Lfy6kGWbNhAUFOR0HBFHJSUl0aZ9R6KLVyY4vKXTcUTSSTxxgOiF/+Wb/01PvcuVXB5teb2AM2fOENGiNSmVGhFU+U6n44ikk3B4F7Hfvs+i+XO55pprnI4j4rj+jz7OrwdOE9Swj9NRRNJJjosmes5rjBnxKk2aNHE6jsdT83oeay3devTkQEphgm7r4nQckXSSoo8TNXcEn0z8kDp18tVlS0Uu6J1x4/lq7iKCWw7S7bolz7HJScQseIP7OrTjkf79nY7jFdS8nmfY8JdYueFXgpo8pusCSp6TkhhPzNyRPPX4I3Tu3MnpOCKOW7JkCc8Pf8l9u24dPiN5T8zKjwm7rjjvvvOm01G8hprXNKZNm847708kuPVQfApc+b2FRbKTtSnELB1H49tq8NKLw5yOA1ze7TpFssv27dvp2KUbwS0HU+CqMk7HEUknZtM8Qk78zuyvv8LPL2+cZuQNdTtvvJN5wLp163iw/yMU7vAyvsFXOR1HJJ2YtdMpVyCWL6Z8mmf2CkybNs3pCJJPHTt2jMbNWuBfrwcFr/Xc61WK9zqzeyNJm77m2/U/UaRIEafjpPKGuq3mFfjrr79o0aYdhRo/hn+pG5yOI5JOzLaV+O76niWbNlCwoK43LPlbQkICzVu3I/6aOgRXj3A6jkg6CUf3EbPkHRbNm02FChWcjuN18n3zGh0dTZPmrfCp1pJClW51Oo5IOvEHdhC36mPWrFxB6dKlnY4j4ihrLQ/0fYhdpywhrXP3jnIimZEce4roOa8z/u03qV+/vtNxvFK+PuY1JSWFjl3u41hAWQrVae90HJF0kk7/Q9S8UXw55VPCwsKcjiPiuNFvjGX+ijUENx+oKwtInmOTE4meP5oHe95Hr149nY7jtfJ18zpoyDOs276PoMb98swxhCJnpSScIXruCIYNHUzbtm2djiPiuDlz5vDaqDcIbvMcPv6BTscROYe1lpjlE7g19HreGD3K6TheLd8eNvDpp5P56LNpFOk6BuOb+/cWFrkUm5JM9KK3aH13fYYMGuR0HBHHbdmyhe739yKk7Qv4FS7ldByRdGI3zKJ4/AFmzvjBkdt15yf5snldvXo1jw98msIdX8O3UN45A1DkrJjVn1OxqA+fTPpQewUk3zt8+DARLVpRsEFfAsre7HQckXRid/5EytYFLN+4nuDgYKfjeL1817zu3r2bNvd0IKjpE/iXuM7pOCLpxPyyjID9G1i4aQP+/v5OxxFxVFxcHE1btiH5xgYEV7nL6Tgi6SQc3k3s8vGsWLqYcuXKOR0nX8hXzeupU6do3KwlfrXuJfCG2k7HEUkn7q+txK/9nFU/rKZ48eJOxxFxlLWW7vf34u/EQgTf3tXpOCLpJEefIGruCD6a8AF169Z1Ok6+kW8OykhKSqJdh05EXXUTQbVaOx1HJJ3EEweJXvAGM6dPJTQ01Ok4Io576ZVXWf7TFoIiBmBMvvm6Eg+RkhhP9LyRPPHIQ3Tt2sXpOPlKvqkGjz/xFFv+Ok5Qw75ORxFJJyU+hui5rzPy1Zdp2rSp03FEHDdjxle8+e77BLfR7bol77HWErNsPA3rVOXVl19yOk6+ky+a1/c/+ICp/zeb4JaDdV1AyXNsSjLRC8bSpV0LHn/sUafjiDhuw4YN9Hm4HyFtnsMvuJjTcUTSiflxBtf4RjH18yk6qdYBXn/M67Jlyxjy3DCKdB6JT0GdASh5T8zKj6latjDvvzvO6Sgijtu/fz/NW7Wh0N398S99o9NxRNKJ+W0VPjtXsGTTBgIDdb1hJ3h187pjxw7u7dyV4BZPU+Cqsk7HEUknJnIBwUd/Y+7Gdfj5efWfo0iGYmJiaNK8FVRtTqGb6jkdRySd+AM7OLNyImtWruDqq692Ok6+5bXflsePH6dJs5b439adgtfptpqS95z5czOJ679i+bofKVq0qNNxRByVkpJCp67d+cevJMF173U6jkg6SaePpN6uu0aNGk7Hyde88pjXxMREWrVtz5kyNQkK08kvkvckHvuLmMVvMWfW19x4o3aNijwz9HnW/rqb4MaP6hhCyXNSEuKInjuCoYOfol27dk7Hyfe8rnm11tLnoX7sOJ5A0J33Ox1HJJ3kM6eJmvM6b40dTYMGDZyOI+K4KVM+Y8KnnxHc+lmMn27XLXmLtSnELH6Llg1v57lnn3E6juCFzet/33yLOUtXEtzsKV1ZQPIcm5xI9PzR9LqvMw/21WXbRH744QcefWIgIW2e1+26JU+KWf0FFUIskz+epL0CeYRXHfM6f/58Xnp9FEW7jMYnoJDTcUTOYa0lZsVE6lYqx1tjxzgdR8Rxe/bsoXW79gRFDMC/5PVOxxFJJ2brcvz/+olFul13nuI1zWtKSgo9e/fBGh8Slr9LQg7M48jOzTkw1awrVfmWyx73n9/Wpf6/YJESFC57Q3ZEuiJpM+WWEjdUy9X5+d7anaR/dnFVzD7+77sf8fXVXgGRAU8NIjomlqI/zyb259nZPv2E2NOc2r8r26ebVSVvrnPZ4x7ZsSHbppVdzs+UG0pWDM/V+flWugNb7Dri10xh5ZpVlChRIlfnL5fmNc2rj48Py5cs5sSJEzk2j23btjFx4sRLDnN2l0Jm/73YcxfTt29fypUrl8Xk/5o5cyaTJ08mOTmZuFNHiTt1lN69e9OpU6fLnuaVmjBhArGxsdk2vaVLl6Z7rnDhwql3rqpTp06uHmv63nvvM/P7LykUd4TlG9YRHKzrDYsAjHtzLHsGPJZj04+Pj2f48OEZDnd+7c3q40uJiIjgjjvuyPTw59u2bRsjR47k0KFDgKtxLF++PJMmTbrsaV6pr776it9//z3bprd+/foLfgdEREQAcN1113H//bl3DssPP/zAiHc/IiX2FP+b+gVVqlTJtXlLJllrM/1Tu3ZtK97jzz//tFdffbUF7JAhQ5yOkyN27NhhCxYsaAH71ltvOZJhyDPPWr8C/nbt2rWOzF9cgA02C/XOG35Us73LmTNnbI8ePSxgCxYs6HScHHH69GnbqlUrC9iaNWs6kmHu3LkWsG++9Y4j8xeXS9VsrzthSzKvfPnyHDx4kBEjRjBmzBjGjfO+OzzddNNNnDlzhv79+zNw4EBmzpyZ6xke7NuHxYsWcdttt+X6vEXEexQsWJDPPvuMzZs3ExcXR4UKFZyOlO1CQkKYN28eS5YsYfPmzbRu3SbXM9SrV4+p06bx5BOP5/q8JXOMq7nNnDp16tgNG3L/WBfJeQ888ACTJ08mMTHRa+/0VPeWW9iwfj0pKSk6YzQfMsZstNY6f8BgLlLN9l6rV6/mzjvvZOnSpTRp0sTpODli+vTpdOvWjV27dnHDDc6fnyG561I1W1teBYBPPvkEgCFDhjicJOesXrUKgI8++sjhJCIiV6Z+/frAv8eFeqOuXbsC0KhRI4eTSF6j5lUA1wkIt956K2+99ZbTUXJMQEAAwcHBPPTQQ05HERG5YgsXLgTg5MmTDifJOaNHj2bfvn0kJyc7HUXyEDWvkmru3LkAJCUlOZwk5yxYsMDpCCIi2aJ58+YAXnm+wlmDBw8G/m3URUDNq6RRsmRJAPbt2+dwkpwTFhYGeHeDLiL5S3ZetiqvOXt+gjcvo2SdmldJ5/Tp005HyDGBgYEAJCYmOpxERCR7REdHOx0hx+WHZZTMU/MqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw/pwOI5IaUlBSWLl1KYmIiAEuXLiUgIICIiAh8fLQOJyKS16xdu5YTJ04A8NNPP7FgwQJuvfVWihcv7nAycZqaV8kXdu/eTYuWLSlYuAQAXXs9SNypo2z/7Tduuukmh9OJiEha1lo6denGocOHAfj2h/UsWNCKyZ99Rs8ePRxOJ07TJifJFypWrMiC+fPxI4nCtVvjZ5OYN3euGlcRkTzIGMPq77/j+vIVKFq9Ib4GRo95g/v/8x+no0keoOZV8o3mzZvz45pV3F0+kLWrv6dly5ZORxIRkYsoX748m9f/yL31KjNtyscMGTwIY4zTsSQP0GEDkq9UqVKFWTO/cjqGiIhkQuHChflowvtOx5A8RlteRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpeRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpeRUREJMf4+voSHh5OtWrVaNOmDSdPnryi6X3zzTeEhYURGhpK9erV+eabb1Jf2759O+Hh4dSsWZNdu3Yxbtw4KleuTPfu3a90MSQPUfMqHkeFUETEcwQGBhIZGcnWrVspVqwY77333mVP6+eff2bQoEHMnj2b7du3M2fOHAYNGsSWLVsAVz1v164dmzdv5sYbb+T9999nwYIFfPnll9m1OJIHqHkVj6NCKCLimW6//Xb2798PwHfffUfr1q1TX3vssceYPHkyAOXLl2f48OHUqlWL6tWrs337dgDGjh3Lc889R4UKFQCoUKECQ4cO5Y033mDBggW8jDvJ0gAADF1JREFU/fbbfPTRRzRq1Ih+/fqxe/du2rZty1tvvZW7Cyo5Ss2reDQVQhERz5CcnMzy5ctp27ZtpoYvUaIEmzZton///owdOxaAX3/9ldq1a58zXJ06dfj1119p2bIl/fr1Y+DAgaxYsYIJEyZQtmxZVqxYwcCBA7N9ecQ5al7FY6kQiojkfWfOnCE8PJzixYtz/PhxIiIiMjVehw4dAKhduzZ79uwBwFqLMeac4S70nHg3Na/icVQIRUQ8x9lDvfbu3UtCQkLqoV5+fn6kpKSkDhcXF3fOeAEBAYDrPIekpCQAqlatyoYNG84ZbtOmTVSpUiUnF0HyGDWv4nFUCEVEPE+RIkUYN24cY8eOJTExkeuvv55t27YRHx/PqVOnWL58eYbTGDRoECNHjkzdALFnzx5GjBjB008/ncPpJS9R8yoeS4VQRMSz1KxZkxo1ajB9+nSuvfZaOnfuTFhYGN27d6dmzZoZjh8eHs7o0aNp06YNoaGhtGnThjFjxhAeHp4L6SWv8HM6gMiVSFsIe/TokVoIK1WqlOVCmJiYSIECBVQIRUSyUXR09DmP586dm/r/MWPGMGbMmHTjnN2gAK7zEL777rvUxx06dEg9DOx8L7300kWnI97DWGszPXCdOnXs+btYxbsYY9i8ebOaN/E6xpiN1to6TufITarZ3s8YQ7t27c65PrWIN7hUzdZhAyIiIiLiMdS8ioiIiIjHUPMqHmnRokXcfPPNVKxYkVGjRqV73VrLgAEDqFixImFhYWzatCnDcc9edqtSpUpERERw4sQJAL788kvCw8NTf3x8fIiMjMz5hRQR8RK5WbPXrVuXWq9r1KjBrFmzcn4BJXdZazP9U7t2bSveDbCbN292OsYlJSUl2RtuuMHu2rXLxsfH27CwMPvrr7+eM8z8+fNt8+bNbUpKil27dq295ZZbMhx38ODBduTIkdZaa0eOHGmHDBmSbt5btmyxFSpUyOEllJwA/H979x/T9J3HcfxV6IwgO/WcMkj9MdK4KBWVACOLKJkJXDZTNHpKL3eaY1tmnHFxFy/+53a53NAsYTHE8MemopmWmJF1WSb+SLZDRaEMS05J1D/4oaXZ0Ojmeduo5Xt/mGvCiR0qUD7l+Uj6x7f9fL7f97d/vPNq++3302o9Rr9LhAc9O/FJssrKyuJdRkxj3bPv3btnhcNhy7Isq7e315o5c2Z0G+aI1bP55hXGaWlpkdPpVFZWliZNmqTy8nL5fL5BY3w+nzZu3CibzabCwkLduXNHoVAo5lyfz6dNmzZJkjZt2jTkHyCOHj0qj8cz+icJAAlirHt2amqq7PYHN1P6+eefWXQmARFeYZxgMKjZs2dHtx0Oh4LB4LDGxJr73XffKSMjQ5KUkZGh77///qFj19XVEV4B4DHEo2c3NzcrOztbixYtUk1NTTTMIjEQXmEca4jbuw21xOtQY4Yz91Gam5uVmpoql8s1zEoBAPHo2S+99JIuX74sv9+vDz744KEVF2E2wiuM43A4dP369ej2jRs3lJmZOawxseamp6crFApJkkKhkGbNmjVon16vl29dAeAxxatnS9KCBQs0ZcoUXbp0aUTPCfFFeIVx8vPzde3aNXV2dqq/v19er1dut3vQGLfbrUOHDsmyLF24cEFTp05VRkZGzLlut1u1tbWSpNraWpWVlUX3NzAwoGPHjqm8vHzsThQAEsBY9+zOzk7dv39fktTd3a0rV65o3rx5Y3fCGHVcBALj2O12VVdXq7S0VJFIRBUVFcrOzlZNTY0kafPmzXr11Vf11Vdfyel0KjU1VQcOHIg5V5J27typ9evX65NPPtGcOXN07Nix6DEbGxvlcDiUlZU19icMAAYb65599uxZVVZW6plnnlFSUpL27dun5557Lj4nj1HB8rAYhOVhkahYHhaJiOVhkahYHhYAAAAJgfAKAAAAYxBeAQAAYAzCK3T//n39blWZilaWSJI2/OGPWr6yVD09PXGuDAAwlG1/+ateLl4p6cFKUy8Xr1T9575fmQUkBu42AIXDYd27d0/NrW2aVlyhno6v1R/uVlISn20AYDxKS5msC2f/qd8U/l6ROyEFvvXrt9OmxrssYEyQTqCUlBR9fapBf/6TR/9uOqICl1OB1mY5HI54lwYAGMI//v43Hdx/QP3/atC0/wTV5m9RcXFxvMsCxgS3ysIggUBALpeLdaCRcLhVFhJRZ2enpk+frmnTpsW7FGBExerZJBQMwv1dAcAcL7zwQrxLAMYclw0AAADAGIRXAAAAGIPwCgAAAGMQXgEAAGAMwisAAACMQXgFAACAMQivAAAAMAbhFQAAAMYgvBqkuLhYJ06cGPTcRx99pC1btgw5vqurSy6XK+Y+u7q6dOTIkeh2a2urtm3bJkk6ePCgtm7dKkmqqanRoUOHnqZ8SdK8efN08+bNp94PAIx39GxgdBBeDeLxeOT1egc95/V65fF4nnif/98I8/LytHfv3ofGbd68WRs3bnzi4wDAREPPBkYH4dUg69at05dffqlffvlF0oMm1tvbq2XLlmnHjh1yuVxatGiR6urqHprb1dWloqIi5ebmKjc3V01NTZKknTt36syZM1qyZImqqqr0zTffaNWqVQ/Nf++99/Thhx+qt7dXS5YsiT6Sk5PV3d2tvr4+rV27Vvn5+crPz9e5c+ckSbdu3VJJSYmWLl2qt956S5ZljeI7BADjBz0bGB2EV4PMmDFDBQUFamhokPTgE/yGDRtUX1+vQCCg9vZ2nT59Wjt27FAoFBo0d9asWTp16pTa2tpUV1cX/ZmpsrJSRUVFCgQC2r59+6/WkJmZqUAgoEAgoDfffFNr167V3Llz9c4772j79u3y+/367LPP9MYbb0iS3n//fS1btkwXL16U2+1WT0/PCL8rADA+0bOB0WGPdwF4PP/7GaqsrExer1f79+/X4cOH5fF4lJycrPT0dK1YsUJ+v185OTnReeFwWFu3blUgEFBycrKuXr36VHWcO3dOH3/8sc6cOSNJOn36tDo6OqKv//jjj7p7964aGxtVX18vSXrttdc0ffr0pzouAJiEng2MPMKrYVavXq13331XbW1t+umnn5Sbmzusi/KrqqqUnp6u9vZ2DQwMaPLkyU9cQygU0uuvv64vvvhCaWlpkqSBgQGdP39eKSkpD4232WxPfCwAMBk9Gxh5XDZgmLS0NBUXF6uioiJ60f/y5ctVV1enSCSivr4+NTY2qqCgYNC8H374QRkZGUpKStLhw4cViUQkSc8++6zu3r077OOHw2GtX79eu3fv1vz586PPl5SUqLq6OrodCASitX366aeSpOPHj+v27dtPduIAYCB6NjDyCK8G8ng8am9vV3l5uSRpzZo1ysnJ0eLFi/XKK69oz549ev755wfN2bJli2pra1VYWKirV69qypQpkqScnBzZ7XYtXrxYVVVVv3rspqYm+f1+7dq1K/oHgN7eXu3du1etra3KycnRwoULVVNTI0natWuXGhsblZubq5MnT2rOnDkj/G4AwPhGzwZGlu1x/kmYl5dntba2jmI5ADA6bDbbt5Zl5cW7jrFEzwZgqlg9m29eAQAAYAzCKwAAAIxBeAUAAIAxCK+Iamho0Isvviin06nKysp4lwMAiIGejYmK8ApJUiQS0dtvv63jx4+ro6NDR48eHXQDawDA+EHPxkRGeIUkqaWlRU6nU1lZWZo0aZLKy8vl8/niXRYAYAj0bExkhFdIkoLBoGbPnh3ddjgcCgaDcawIAPAo9GxMZIRXSJKGut8vSwQCwPhEz8ZERniFpAef2q9fvx7dvnHjhjIzM+NYEQDgUejZmMgIr5Ak5efn69q1a+rs7FR/f7+8Xq/cbne8ywIADIGejYnMHu8CMD7Y7XZVV1ertLRUkUhEFRUVys7OjndZAIAh0LMxkdmGum7mUVgnG4CpYq2Tnajo2QBMFatnc9kAAAAAjEF4BQAAgDEIrwAAADAG4RUAAADGILwCAADAGIRXAAAAGIPwCgAAAGMQXgEAAGAMwisAAACMQXgFAACAMQivAAAAMAbhFQAAAMYgvAIAAMAYhFcAAAAYg/AKAAAAY9gsyxr+YJutT1L36JUDAKNmrmVZM+NdxFiiZwMw2CN79mOFVwAAACCeuGwAAAAAxiC8AgAAwBiEVwAAABiD8AoAAABjEF4BAABgDMIrAAAAjEF4BQAAgDEIrwAAADAG4RUAAADG+C/6lk7bd3X3kAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "Stabilization.plot_sankey_report(margin = 0.5)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Sensitivity analysis to precipitation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, 'Percent of Incoming PFAS that \\n remains in the Soil (%)')"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEGCAYAAACZ0MnKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3hUdbrA8e+bRkhIb7QUCKEXgYCg2Bv2xY7utS9id/fu6urdvdetctfdu7prW9Z1LWvBrmsDFRWsQJDeDSWhJSH0mvLeP84JDiGZTJKZnJT38zznmTlnfnPOOyPmnfOroqoYY4wxgQjzOgBjjDFthyUNY4wxAbOkYYwxJmCWNIwxxgTMkoYxxpiARXgdQKilpqZqTk6O12EYY0ybUlBQUKaqabWPt/ukkZOTw7x587wOwxhj2hQRWV/XcaueMsYYEzBLGsYYYwJmScMYY0zALGkYY4wJmCUNY4wxAWs1SUNEnhKREhFZUs/rIiJ/EZE1IrJIREa0dIzGGNPRtZqkATwNjPfz+tlAnrtNAh5vgZiMMcb4aDVJQ1VnAeV+ilwIPKuOr4FEEekWqnjeWrCRl+ZsCNXpjTGmTWo1SSMAPYAin/1i99hRRGSSiMwTkXmlpaVNutg7izbzxxmrqKiqbtL7jTGmPWpLSUPqOFbnClKqOlVV81U1Py3tqFHwAbliVCZlew4yc0VJk95vjDHtUVtKGsVAps9+T2BTqC52Ut80MuI7MW1uUcOFjTGmg2hLSeNt4Gq3F9UYYKeqbg7VxSLCw7h0ZCafrixh8879obqMMca0Ka0maYjIi8BXQD8RKRaRG0RksohMdou8BxQCa4C/A7eEOqbL8jOpVnhlXnGoL2WMMW1Cq5nlVlUnNvC6Are2UDgAZKXEcHyfFKbNLeK2U/oQFlZXs4oxxnQcreZOo7W6fFQWG3fs54vvyrwOxRhjPGdJowFnDswgMSaSl6xB3BhjLGk0JDoynAnDezBj6RbK9x7yOhxjjPGUJY0AXDEqi4oq5fX51iBujOnYLGkEoF/XOIZnJTJtbhFOe7wxxnRMljQCdMWoTFaX7GH+hh1eh2KMMZ6xpBGg84Z2JzYqnGlzbRJDY0zHZUkjQLGdIjh/WHf+vXAzuw9UeB2OMcZ4wpJGI1w+KpP9FVW8syhks5cYY0yrZkmjEY7JTKRfRpyN2TDGdFiWNBpBRLh8VCYLi3awfPMur8MxxpgWZ0mjkSYM70FUeJhNmW6M6ZAsaTRSUmwUZw3uyhvfbuRARZXX4RhjTIuypNEEV4zKZOf+CqYv3eJ1KMYY06IsaTTB2N4pZCZ3tioqY0yHY0mjCcLChMvzM/nyu22s37bX63CMMabFWNJooktGZhIm2N2GMaZDsaTRRF0Tojm1fzqvFBRTWVXtdTjGGNMiLGk0w+WjsijdfZBPVpZ6HYoxxrQISxrNcEq/NNLjOtkkhsaYDqPBpCEivQI51hFFhIdxyciezFxRwpadB7wOxxhjQi6QO43X6jj2arADaauuGJUFwNNfrvM2EGOMaQER9b0gIv2BQUCCiFzk81I8EB3qwNqKrJQYzh7Sjee/Xs8tp+QSHx3pdUjGGBMy/u40+gHnAYnA+T7bCOBHoQ+t7bj5pFx2H6zkhW+sbcMY077Ve6ehqm8Bb4nIWFX9KtSBiMh44GEgHHhSVafUej0JeArIBQ4A16vqklDHFYjBPRIY1yeVpz5fy3XH59ApItzrkIwxJiQCadP4VkRuFZHHROSpmi2YQYhIOPAocDYwEJgoIgNrFbsPWKCqQ4GrcRJMqzH5pFxKdh/kjfkbvQ7FGGNCJpCk8RzQFTgL+AzoCewOchyjgTWqWqiqh4CXgAtrlRkIfAygqiuAHBHJCHIcTXZ8nxQG94hn6qxCqqrV63CMMSYkAkkafVT1l8BeVX0GOBcYEuQ4egC+83EUu8d8LQQuAhCR0UA2TgI7iohMEpF5IjKvtLRlBt6JCJNPyqWwbC8zbPZbY0w7FUjSqHAfd4jIYCAByAlyHFLHsdo/16cASSKyALgd+BaorOtkqjpVVfNVNT8tLS24kfpx9uBuZKfE8MRn36FqdxvGmPYnkKQx1W2E/iXwNrAM+EOQ4ygGMn32ewKbfAuo6i5VvU5Vj8Fp00gD1gY5jmYJDxMmndibhcU7+apwm9fhGGNM0DWYNFT1SVXdrqqfqWpvVU1X1SeCHMdcIE9EeolIFHAFToI6TEQS3dcAbgRmqWqrW6j74hE9Se3SiSc+K/Q6FGOMCbp6u9zWEJFOwMU4VVKHy6vqr4MVhKpWishtwHScLrdPqepSEZnsvv4EMAB4VkSqcO52bgjW9YMpOjKc647P4cHpK1m6aSeDuid4HZIxxgSNNFT3LiIfADuBAuDwotiq+qfQhhYc+fn5Om/evBa95s79FRw/ZSan9k/nLxOHt+i1jTEmGESkQFXzax9v8E4D6Kmq40MQU7uV0DmSK4/N4snZhfz0zH5kpcR4HZIxxgRFIA3hX4pIsLvYtnvXH9+L8DDh77OtbcMY037UmzREZLGILALGAfNFZKWILPI5bvzomhDNhOE9eHleEWV7DnodjjHGBIW/6qnzmnNiEYkFDqhqVYOF26lJJ+bySkExz3y5jv88s5/X4RhjTLPVe6ehqutVdT3w25rnvsdqlxeRMBG5UkTeFZESYAWwWUSWisiDIpIXuo/ROvVJ78KZAzN49qv17D1Y5zhEY4xpUwJp0xjku+NOLjiyjnKf4MxAey/QVVUzVTUdOAH4GpgiIj9sZrxtzuSTctm5v4IX59i06caYts/fIkz34sws21lEagbRCXAImFrHW05X1YraB1W1HGf1v9dEpMOtUDQ8K4ljeyXz5Oy1XD02h6gIW5bdGNN2+aueekBV44AHVTXe3eJUNUVV762j/BEJQ0SiReRGEbldRFLqKtNRTD45ly27DvDWAps23RjTtgUyjchRCSJANQsqHQDebOI52oWT+6bRv2scf5tVSLVNm26MacOCVlciIi+ISK7PoWTgeeBFIClY12mLaqZNX1Oyh49XlHgdjjHGNFkwK9h/AfxGRP4oIgnAH3EmHZwB3B/E67RJ5w3tRo/Ezjwyc7VNm26MabP8De6L8W24FpF+IvJjEbmorvLuqntX4lRFTcNZje8MVT1OVV8NduBtTUR4GHeensfC4p28s2iz1+EYY0yT+LvT+AB3sSUR6QN8BfQGbhWRB2oXFpEkEbkVZ1nWy3AmOZwuIs0aJNieXDyiJ/27xvGH6Ss4WNlhxzwaY9owf0kjSVVXu8+vAV5U1duBs6l7tPibwEEgGnhOVZ8FzgdGisjbdZTvcMLDhHvPGUBR+X6e+2q91+EYY0yj+UsavhXvpwIfAqjqIaC6jvIpwAs4jd893LL7VfVXwE1BibYdOKlvGifkpfLXmWvYua9D9kA2xrRh/pLGIrdR+8dAH5wGbUQksZ7y/4OTWN4Afu77gqpaJb6Pe88ewK4DFTz66RqvQzHGmEbxlzR+BJThtGucqar73OMDcXpGHUFVX1PV41X1RFX9KOiRtiMDu8dz0fCePP3FOorK9zX8BmOMaSX8jQjfr6pTVPVOVV3o81IR0LV2eRGZKiKD6zqXiMSKyPUiclXzQ24ffnpWX0TgjzNWeh2KMcYELKBxGiKSKiI3i8gs4FMgo45ijwH/LSLLReQVEXlMRJ4SkdnAl0Ac0OG73tboltCZG8b14q0Fm1hUvMPrcIwxJiD+JiyMAyYAVwJ9cdoqeqtqz7rKq+oC4DIR6QLkA92A/cByVbWf03W4+eRcps0t4nfvLuelSWMQEa9DMsYYv/zdaZQANwC/A3JV9T9xZrj1S1X3qOqnqvqiqr5pCaN+cdGR3Hl6Ht+sLefj5Ta9iDGm9fOXNO7DGXPxOHBvrXmlTJBMHJ1F79RYHnh/OZVVdfVkNsaY1sNfQ/ifVfVY4AKcdTTeBLqLyD0i0relAmzvIsPDuHt8f74r3cu0eUVeh2OMMX4FMjV6oar+TlWHAKOABOD9YAciIuNFZKWIrBGRn9fxeoKI/FtEFrpLyF4X7Bi8ctagDPKzk/jzh6vZY8vCGmNaMX8TFo6pfUxVF6vqfap6VFWV+wf97fo2f0G4S8g+ijNFyUBgoogMrFXsVmCZqg4DTgb+JCJRDX7CNkBEuO/cAZTtOcjUWYVeh2OMMfWqt/cUThfaEQAi8pWqjm3gXEcN+GuE0cAaVS10r/cScCGwzKeMAnHidDHqApQD7eZn+YisJM4d0o2/zyrkqmOzyIiP9jokY4w5ir+k4dv/s8G/YKr6WTPi6IEzaLBGMXBsrTKP4KzPsQlnzMflqtquWo7vHt+PGcu28OcPVzHl4qFeh2OMMUfx16YR5k53nuLzPLlmq11YRF52HxeLyCKfbbGILGogjroGKNReqegsYAHQHTgGeERE4us8mcgkEZknIvNKS0sbuHTrkZ0Syw/HZPPyvCJWbd3tdTjGGHMUf3caCUAB3/9Bn+/zmuKsreHrTvexKetnFAOZPvs9ce4ofF0HTFFn2bs1IrIW6A/MqX0yVZ0KTAXIz89vU8vk3XFqHq8WFPPAe8v553WjvQ7HGGOO4K/LbY6q9lbVXnVstRPG4ZlsVXU9cAAY4m773WP+zAXyRKSX27h9BU5VlK8NwGkAIpIB9APaXatxUmwUt53Sh09WlvLJShvwZ4xpXfz1nkoXkYdE5B0R+X19VUF1vO8ynF//l+Ks4PeNiFzi7z2qWgncBkwHlgMvq+pSEZksIpPdYr8BjhORxcDHwD2qWhZITG3Ntcfn0Ce9C794Y4l1wTXGtCri1PbU8YLIBzjVU7NwqpziVPXaBk8oshBnbfASdz8N+MjtKtvi8vPzdd68eV5culkK1m/nkie+5D/GZPPrC+ucPNgYY0JGRApUNb/2cX8N4V1V9b9Udbq7zGug3XnCahKGa1sD1zF1GJmdxLXH5fDsV+uZs7bc63CMMQbw/8dcfHtMAeG19uvzgYhMF5FrReRa4F3gvWAG3VH89Mx+9EzqzM9fW8SBiiqvwzHGGL9Jo6b3VM0Wj9ODqgCot75HVX8G/A3nzmQYMFVV7wlWwB1JbKcIHrhoCIVle3n449Veh2OMMfV3uVXVnKaeVFVfdxdsOgGn15NpohPy0rgsvydTZxVy7pBuDO6R4HVIxpgOLGhtDW4vq8Hu827AYuB64FkRuStY1+mI/uucgSTHRnH3q4uosOnTjTEeCmYDdS9VXeI+vw74UFXPB8bgJA/TRAkxkfzmwsEs27zLJjQ0xngqmEmjwuf5abiN36q6G7Cfx800fnBXzh3SjYc/Ws2akj1eh2OM6aAaTBq+8035bJF1FC0SkdtFZALO7LgfuO/vDNRV3jTS/RcMonNUOPe8tojq6jY1O4oxpp0I5E5jPlAKrAJWu8/Xish8ERnpU+4GYBBwLc4MtDvc42OAfwYt4g4sLa4T/33eQArWb+fZr9Z5HY4xpgMKJGl8AJyjqqmqmoKzUNLLwC04a24AoKolqjpZVS9U1Rk+xz9R1eastWF8XDSiByf1TeMP01dSVL7P63CMMR1MIEkjX1Wn1+y4CeFEVf0a6BSyyEydRITfXzQEAe57YzH1TQNjjDGhEEjSKBeRe0Qk293uBra7S7RaA7cHeiR25udn92f26jJeLSj2OhxjTAcSSNK4Emd9izeBt4As91g4ziy2xgNXHZvNqJwkfvPOMkp2H/A6HGNMB9Fg0lDVMlW9XVWHq+oxqnqbqpaq6iFVXVO7vIj0FZGPRWSJuz9URH4RiuA7srAwYcrFQzlQWc19ry+xaipjTIsIpMttXxGZKiIzRGRmzebnLX8H7sUdt6Gqi3AWVTJBlpvWhXvG9+ej5Vv5x+drvQ7HGNMB+FvutcYrwBPAk0AgU63GqOockSOW/baVhELk+uNzmLu2nCnvr2B4ViIjs/1NQGyMMc0TSJtGpao+rqpzVLWgZvNTvkxEcnHWEcddtW9zMII1RxMR/nDpUHokdebW579l256DXodkjGnHAkka/xaRW0Skm++ocD/lb8WZGr2/iGwE7gJuDkawpm7x0ZE8dtUIyvcd4q5pC6iy0eLGmBAJJGlcA/wM+JLv19bwt55GoaqeDqQB/VV1nKquC0Ksxo9B3RP41QWDmL26jEdmHtU/wRhjgqLBNg1V7dWYE4pIJ+BiIAeIqGnbUNVfNyE+0whXjMpk7tpyHvp4FSOzkxiXl+p1SMaYdqbepCEip6rqTBG5qK7XVfX1et76FrAT547EKthbkIjw2wmDWbJpJ3e+9C3v3nECXROivQ7LGNOO+LvTOAmYCZxfx2sK1Jc0eqrq+OYGZpomJiqCx64awQWPfMHtL87nhR+NITI8mDPgG2M6Mn/Lvf6P+3hdI8/5pYgMUdXFzYrMNFmf9DgeuGgId760gD9OX8m95wzwOiRjTDvRYJuGiCQCV+O2UdQcV9U7apVbjHMHEgFcJyKFONVT4hTXocEL2zTkwmN6MGdtOX+bVUh+TjJnDMzwOiRjTDsQyOC+94Cvcdb89jdB4XnNCURExgMP48xp9aSqTqn1+s+Aq9zdCGAAkKaq5c25bnv2y/MGsrB4B//58gLeveMEMpNjvA7JGNPGSUNzFonIfFUdEfAJRZ5T1f9o6Fit18NxFnk6AygG5gITVXVZPeXPB36sqqc2FE9+fr7Om1dvD+F2b8O2fZz719nkpMTyyuSxREeGex2SMaYNEJECVc2vfTyQFtLnRORHjRjcN6jWhcOBkfWUrTEaWOOO8TgEvARc6Kf8RODFAGLv8LJSYvjTpcNYvHEnv323zhxsjDEBCyRpHAIeBL7Cz+A+EblXRHYDQ0Vkl7vtBkpwuuH60wMo8tkvdo8dRURigPHAa/WdTEQmicg8EZlXWlrawKXbvzMHdWXSib3519cbeOGbDV6HY4xpwwJp0/gJ0EdVy/wVUtUHgAdE5AFVvbeRcUgdx+qrNzsf+MJfW4aqTgWmglM91chY2qW7z+rHqq27+cWbi8mI78RpA6xh3BjTeIHcaSwFAl6MugkJA5w7i0yf/Z7ApnrKXoFVTTVaRHgYj145gkHdE7jthW9ZULTD65CMMW1QIEmjClggIn8Tkb/UbEGOYy6QJyK9RCQKJzG8XbuQiCTgDDpsqLrL1CG2UwRPXTuKtLhOXP/0XNaV7fU6JGNMGxNI0ngT+B1HTljob2r0RlPVSuA2YDqwHHhZVZeKyGQRmexTdAIwQ1Xtr10TpcV14pnrRwNwzT/nUGZTqRtjGqHBLrcA7q//vu7uSlWtaKB8OJDBkYMBPWmB7ehdbuvz7YbtTPz71/TLiOPFSWOIiQqkecsY01E0ucutiJwMrAYeBR4DVonIiX7K3w5sBT4E3nW3d5oWtgmV4VlJPDJxBIs37uTW5+dTWeVv3KYxxjgCqZ76E3Cmqp6kqicCZwF/9lP+TqCfqg5S1SHuZlOItEKnD8zgNz8YzCcrS/mvN5YQyF2nMaZjC6ROIlJVV9bsqOoqEYn0U74IZ2p00wZcdWw2W3Ye4K8z19AtMZq7Tu/b8JuMMR1WIEljnoj8A3jO3b8K/w3hhcCnIvIuPutpqOr/NTlKE1I/OaMvm3ce4KGPVtMtIZrLR2V5HZIxppUKJGncjLPu9x04g/Bm4bRt1GeDu0W5m2nlRIQHLhpCye6D3PfGEtLjojmlf7rXYRljWqFAJiyMBQ6oapW7Hw50UtWAB/x5yXpPBW7vwUoun/oV35Xs5aVJYxiWmeh1SMYYjzRnwsKPgc4++52Bj+q4wEPu479F5O3aW1MDNy2nZvBfalwUVz81h0XFNmrcGHOkQKqnolV1T82Oqu5xJw2srabN449Bicx4Ij0umhduHMPEv3/NVX//hqevH83I7CSvwzLGtBKB3GnsFZHD62mIyEhgf+1CqlrgPn5W1xa8kE2oZSbH8PJNY0npEsXV//iGbwq3eR2SMaaVCCRp3AW8IiKzRWQ2MA1nyg/TjnVP7MzLN42la0I01/xzDl+s8TvJsTGmg2gwaajqXKA/Ti+qW4ABNXcVpn1Lj49m2k1jyUmJ5bqn5/LJyhKvQzLGeCyQOw2AUcBQYDgwUUSuDuRNIhImIvFNDc54L7VLJ1780Rjy0rtw07MFzFi6xeuQjDEeCmTuqedwGrfH4SSPUcBR3bB8yr8gIvFuV91lwEoR+VmQ4jUeSIqN4oUbxzCgezy3PD+fdxdt9jokY4xHAuk9lQ8M1MAnJhqoqrtE5CrgPeAenBHkDzYxRtMKJMRE8q8bRnPdP+dy+4vzqag6hh8Mr3NFXmNMOxZI9dQSoGsjzhnpzk31A+Atdxp1mwmvHYiLjuSZ60dzbK8UfvzyAl6eV9Twm4wx7UogSSMVWCYi0wMcrPc3YB0QC8wSkWxgV/NDNa1BzQDAcX1SufvVRTz/zXqvQzLGtKBAphE5qa7jjRl7ISIR7up8Lc6mEQmNAxVV3PL8fGauKOGu0/O487Q8RMTrsIwxQVLfNCINtmk0dmCeiHQCLgZyap3/1405j2ndoiPDeeKHI7n39cU89NFqCkv38odLhhIdGe51aMaYEKo3aYjIbupuixBAVbW+rrRv4aynUYDP1Oim/YmKCOOPlw4lNz2WP3ywkg3l+5h69UjS46K9Ds0YEyL1Jg1VjWviOXuq6vgmvte0MSLCLSf3oXdqLHdNW8CER7/kyWvyGdDNhucY0x4FOrivMb4UkSEhOK9pxcYP7sYrNx1HZXU1lzz+JTNXbPU6JGNMCIQiaYwDCkRkpYgsEpHFIrIoBNcxrcyQngm8des4eqXFcuMz83hydqGtO25MO+OvTaOTqjalTeLsZsRj2riuCdG8fNNYfjJtIb99dznfle7l1xcOIjI8FL9PjDEtzd//yV/B4WlEGuQzx9TuejbTQcRERfDYVSO45eRcXpyzgWuemsPOfRVeh2WMCQJ/XW6jROQa4DgRuaj2i6r6eq1DLwDn4fSaUpxeVoeLA739BSIi44GHgXDgSVWdUkeZk4GHgEigTFXrHENivBcWJtw9vj+5aV34+euLmPDYFzx5TT6907p4HZoxphnqHdwnIuOAq4DLgNojwFVVrw9aEM6646uAM4BiYC4wUVWX+ZRJBL4ExqvqBhFJV9UG5+q2wX3em7uunJueK+BQZTW/v2gIFwzr7nVIxpgGNHpwn6p+DnwuIvNU9R+NvFgSkAcc7rCvqrP8vGU0sEZVC933vwRciDNLbo0rgddVdYN7PlvcoY0YlZPM27cdzx0vfssdL37LF6vL+J8LBhITFch8mcaY1iSQ1snnROQOEXnV3W53JySsk4jcCMwCpgO/ch/vb+AaPQDf2e+K3WO++gJJIvKpiBT4W9NDRCaJyDwRmVdaWtrApU1L6JkUw7SbxnLrKbm8XFDEBY98wYotNiWZMW1NIEnjMWCk+/gYMAJ43E/5O3HW3FivqqfgLNzU0F/uuiYtql1vFuHGcS5wFvBLEelb18lUdaqq5qtqflpaWgOXNi0lMjyMn53Vn+euP5ad+yu44JEv+NfX661brjFtSCBJY5SqXqOqM93tOpykUJ8DqnoADnfbXQH0a+AaxUCmz35PYFMdZT5Q1b2qWoZzNzMsgPhNKzMuL5X37zyBMb1T+MWbS7jl+fnWu8qYNiKQpFElIrk1OyLSG6jyU77YbbR+E/hQRN7i6ARQ21wgT0R6iUgUcAVHN76/BZwgIhEiEgMcCywPIH7TCqV26cTT147i3rP78+GyrZzzl9kUrN/udVjGmAYE0hL5M+ATESnEqUbKBq6rr7CqTnCf3i8inwAJwAf+LqCqlSJyG077RzjwlKouFZHJ7utPqOpyEfkAWARU43TLXRJA/KaVCgsTbjopl9G9krn9xW+57G9f8ZMz+nLzSbmEhdk068a0Rg2upwGHpzvvh5M0VjQ0UtztPZWJT1JS1fnNC7VprMtt27DrQAX3vr6YdxdtZlyfVP5wyVC6J3b2OixjOqz6utwGlDQaeaHfANcChTh3BOCM6zg1qBcKkCWNtkNVmTa3iF/9exnhYcI94/tx1bHZdtdhjAdaMmmsBIao6qGgnriJLGm0PRu27eO+Nxbz+ZoyRuUk8cBFQ+mTbiPJjWlJ9SWNUMwitwRIDMF5TQeRlRLDczeM5sFLhrJq6x7OeXg2f/14NYcqqxt+szEmpBpMGiLycSDHfDwAfCsi00Xk7ZqtOUGajkdEuDQ/kw9/ciJnDMzgTx+u4oJHPmdh0Q6vQzOmQ/M3NXo0EAOkug3bNRXL8YC/yYOeAf4XWMz3bRrGNEl6XDSPXjWCC5du4ZdvLWHCY19w/fG9+MmZfW0aEmM84O//upuAu3ASRAHfJ41dwKN+3lemqn8JTnjGOM4c1JUxuSlMeX8FT36+lunLtvDAhKGMy0v1OjRjOpQGG8JF5HZV/WvAJxT5P+AgzuC8w11zrcutCZavC7dx7+uLWVu2l4uG9+Du8f3pmhDd8BuNMQFrVu8pETkOyOHIcRfP1lP2kzoOW5dbE1QHKqr468zV/H3WWsLDhMkn5TLpxN50jgr3OjRj2oUmJw135b5cYAHfTx+iqnpH0KMMAUsa7duGbfuY8sFy3lu8hW4J0dw9vh8XDuthYzuMaabmJI3lwEANcECHiGQAvwe6q+rZIjIQGNvYNTmCxZJGxzBnbTm/eWcZizfuZFhmIv993gBGZid7HZYxbVZzxmksAbo24lpP48whVdPDahVOg7oxITO6VzJv3Xo8f7p0GFt27ufix7/ithfmU7x9n9ehGdOuBNJnMRVYJiJzOLJh+4L6yqvqyyJyr1uuUkT8zYprTFCEhQkXj+zJ2UO68sRnhUyd9R0zlm3lxnG9uOWUPnTpZF10jWmuQP4vur+R59wrIim4iyiJyBhgZyPPYUyTxURF8JMz+nLFqEwenL6Sxz79jpfnFXPrKblMHJ1FdKQ1lhvTVIH2nsoG8lT1I3cti3BV3V1P2RHAX4HBOFVbacAlqrooeGEHzto0zIKiHUx5fzlfF5aTEd+Jm0/K5QpLHsb41ZyG8B8Bk4BkVc0VkTzgCVU9rY6yYcAYYA7fT6W+UlU9W5bNkoap8dV32/jzR6uYs7acrvEwyhYAABT+SURBVPHR3HJKLpflZ1ryMKYOzUkaC4DRwDeqOtw9tlhVh9RT/itVHRuEmIPCkobxpap8VbiNhz5czZx1TvK49ZRcLhuVSacISx7G1GhO76mDvtOci0gEbntFPWaIyMUiYh3lTasjIhyXm8q0m8bw/I3H0jOpM798aymnPPgp//p6PQcrrc+GMf4EcqfxB2AHcDVwO3ALsExV/6ue8ruBWKASOIBTRaWqGh/EuANmdxrGH1Xl8zVl/PnDVczfsIPuCdHcfHIul4zMtNHlpkNrTvVUGHADcCZOApiOsz53cFdvChFLGiYQqsrs1WU89JGTPBJjIvnhsdlcPTab9Hib18p0PM1JGrHAAVWtcvfDgU6q2iZGTVnSMI2hqsxbv50nZxcyY9lWIsPCuOCY7tx4Qi/6d/XkZtkYT9SXNAIZp/ExcDqwx93vDMwAjgteeMa0DiLCqJxkRuUks65sL099sZZX5hXzakExJ+SlcuMJvTkxLxVrsjMdVUC9p1T1mIaOtVZ2p2Gaa8e+Qzz/zQae+XIdJbsP0i8jjhtO6MWFx3S3Hlem3WpO76m97oC9mhONBPYHMzhjWrPEmChuPaUPn99zKn+6dBgicPerizh+yic8OH0FReVtoqbWmKAI5E4jH5gGbHIPdQMuV9WCoAYiMh54GAjHaWifUuv1k4G3gLXuoddV9dcNndfuNEywqSpfrNnG01+uZeaKEhQ4MS+NiaOzOG1AOpHhgfwWM6Z1a1KbhttzKgroz/cjvFcEe4S327j+KHAGUAzMFZG3VXVZraKzVfW8YF7bmMYSEcblpTIuL5VNO/YzbW4R0+YWMflfBaTHdeLyUZlcPiqTnkkxXodqTNAFcqcR8hHeIjIWuF9Vz3L3a2bIfcCnzMnATxubNOxOw7SEyqpqPllZyotzNvDJyhIATuqbxpWjszi1fzoRdvdh2pjm9J6aISIX41QHhWpsRg+gyGe/GDi2jnJjRWQhTlXZT1V1aV0nE5FJOPNlkZWVFeRQjTlaRHgYZwzM4IyBGWzcsZ9pczYwbV4Rk54rICO+E5eM7MmE4T3pk97F61CNaZZA7jRqRnhX4TSAB32Et4hcCpylqje6+/8BjFbV233KxAPVqrpHRM4BHlbVvIbObXcaxiuVVdXMXFHCC3M2MGtVKdUKQ3okMGF4D84f1p20uE5eh2hMvZo8uK8lBFI9Vcd71gH5qlrm79yWNExrULL7AG8v2MQb325k6aZdhIcJJ+SlMmF4D84c2NWmLDGtTnNGhAtwFdBLVX8jIplAN1WdE8TgInCWhT0N2AjMBa70rX4Ska7AVlVVERkNvApkN1RlZknDtDartu7mzW838ua3G9m08wCxUeGMH9yNCcN7MDY3hfAwGzhovNecpPE4UA2cqqoDRCQJmKGqo4Ic4DnAQzhdbp9S1d+JyGQAVX1CRG4DbsaZCHE/8BNV/bKh81rSMK1VdbUyZ105b8zfyHuLN7P7YCUZ8Z04e3A3zh7clfycZEsgxjPNSRrzVXWEiHzrs57GQlUdFqJYg8qShmkLDlRU8fHyEt5asJHPVpVysLKa1C6dOGtQBucM6caxvZKtB5ZpUc3pPVXhjqOoWfM7DefOwxgTJNGR4Zw7tBvnDu3G3oOVfLKyhPcXb+H1+Rt5/psNJMVEcubArpw9pCvH5aYSFWEJxHgjkKTxF+ANIF1EfgdcAvwipFEZ04HFdorgvKHdOW9od/YfquKzVaW8v2Qz7y7ezLR5RcRHR3D6wAzOGtSVE/JSiYkK5H9jY4IjoN5TItIfp5FagI9VdXmoAwsWq54y7cWBiiq+WFPGe4u38OGyLew6UElURBhje6dw2oB0TumXTmayjUI3wdHoNg0RiQYmA32AxcA/VLUypFGGgCUN0x5VVFUzd105M5eXMHNFCYVlewHolxHHqQPSOa1/OsOzkqwh3TRZU5LGNKACmA2cDaxT1btCGmUIWNIwHUFh6R5mrijh4+UlzF1XTmW1khgTySn90jmlfzrj+qSSHBvldZimDWlK0lisqkPc5xHAHFUdUWfhVsyShulodu6vYPbqUmYuL+GTlSVs31eBiDMafVyfVE7IS2NEdqKtBWL8akrSmO+bJGrvtxWWNExHVlWtLCjaweery/h8TSnzN+ygqlrpHBnOmN7JjMtL48S8VPqkd7HVCM0RmpI0qoC9Nbs4y7zuIwRzT4WSJQ1jvrf7QAVfF5Yze3Upn68uO9wW0jU+2pnuvU8qx/ZOpltCZ48jNV5r1XNPhZIlDWPqV7x9H5+vLmP2mjK+WFPGjn3OUjnZKTGM6ZXCmNxkju2VQvdESyIdjSUNY4xfVdXKii27+LqwnG8Kt/HN2nJ27neSSFZyDGN6OwlkTG4KPSyJtHuWNIwxjVJdrazYspuvC7fxzVonidTcifRM6syonGRGZicxMjuJvhlx1r23nbGkYYxplupqZeVWJ4nMWVvOvPXbKd19EIC4ThEMz05iZFYS+TlJHJOZSGwnG6nellnSMMYElapSVL6feevLKVi/nYL121m5dTeqECYwoFs8+dlJjMhOYljPRLJTYqyHVhtiScMYE3I791ewoGgHBeucO5EFRTvYd6gKgMSYSIb1TGRYZiLDMxMZ2jOBlC62emFr1ZxZbo0xJiAJnSM5qW8aJ/VNA5wlb1dt3cPC4h0sLNrBgqIdPDJzNdXub9Ws5BiGZSYyrGcCx2QmMrB7vE3A2MrZnYYxpkXtPVjJ4o07WVi0g4XFO1iwYQebdh4AnGqt3LQuDO6R4Gzd4xnUI4Eu1j7S4uxOwxjTKsR2imBM7xTG9E45fKxk1wEWFu9kycadLN20ky+/K+ONbzcCIAK9UmIZ1COBIT3iGdw9gUHdE0iIifTqI3RoljSMMZ5Lj4/mjIHRnDEw4/Cxkt0HWLpxF0s27mTxxp3MX7+dfy/cdPj17gnRDOgWz4Bu8Qzs7jxmJ8cQZl1/Q8qShjGmVUqPiya9fzSn9E8/fGzbnoMs2bSL5Zu/3z5dVUqV20gSExVOv65x3yeTbnH0zYgjLtruSoLF2jSMMW3agYoqVm/dw/LNu1jmk0x2Hfh++Z8eiZ3pm9GFvl3j6JfhJJI+6V2IjrSZfutjbRrGmHYpOjKcIT0TGNIz4fAxVWXTzgMs37SLVSW7WbVlNyu37uGLNds4VFUNOI3u2Smx9M3oQr+MOPLcRNIrNdaSiR+WNIwx7Y6I0COxMz0SO3O6TztJZVU167btY9XW3azcstt53LqbD5dtPdwNOEycrsB90rvQJ91JJHnpXchN72K9uLCkYYzpQCLCw9xk0IVzhnQ7fPxARRWFpXtZU7qHNSV7WFOymzUle/hsVSkVVd9X4XdLiKZPehd6p8bSO60LvdNiyU3rQtf46A7TAG9JwxjT4UVHhjOwu9MLy1dFVTUbyve5ieT77dWCYva6I90BOkeGk5Ma6yQRn4TSKzW23TXCt5qkISLjgYeBcOBJVZ1ST7lRwNfA5ar6aguGaIzpYCLDw8hN60JuWhfOGvT9cVWlZPdBvivdQ2HpXmcr28Pi4p28v3jz4aougNQuUeSkxJKT6iSR7JSYw/ttsbqrVUQsIuHAo8AZQDEwV0TeVtVldZT7X2B6y0dpjDEOESEjPpqM+GiOy0094rWDlVWs37aPwtI9rC3bx7qyvazdtpfZq0t5taD4iLJpcZ3IcZNIdkoMWSmxZCfHkJ0SQ2JMVEt+pIC1iqQBjAbWqGohgIi8BFwILKtV7nbgNWBUy4ZnjDGB6RQRTl+3W29t+w5Vsq5sH+u27XW2sr2sK9vHZ6tKKXGnma8RHx1BdkosWSkxZCXHkJ0cc/h5t4TOnq1f0lqSRg+gyGe/GDjWt4CI9AAmAKfSQNIQkUnAJICsrKygBmqMMU0VExVRZ9sJwP5DVWwo38f6bXvdx31sKN/Hsk27mLF0yxEN8hFhQvfEzmQmdyYrOYaeSTFkJseQmdSZzOQYUmKjQjYNfWtJGnV9utqjDh8C7lHVqoa+DFWdCkwFZ3BfUCI0xpgQ6uyOZu/X9eg7lKpqZdOO/Wwo30dRuZNMirbvp6h8Hx8u20rZnkNHlI+JCqdnUmeevHoUWSkxQY2ztSSNYiDTZ78nsKlWmXzgJTdhpALniEilqr7ZMiEaY4w3wsPEuZNIrjsB7D1YSbGbRIq276OofD9F2/eFZFLH1pI05gJ5ItIL2AhcAVzpW0BVe9U8F5GngXcsYRhjjDNzcH13KcHWKpKGqlaKyG04vaLCgadUdamITHZff8LTAI0xxgCtJGkAqOp7wHu1jtWZLFT12paIyRhjzJHCvA7AGGNM22FJwxhjTMAsaRhjjAmYJQ1jjDEBs6RhjDEmYJY0jDHGBKzdrxEuIqXA+ia+PRUoC2I4bYF95o7BPnP719zPm62qabUPtvuk0RwiMq+uhdXbM/vMHYN95vYvVJ/XqqeMMcYEzJKGMcaYgFnS8G+q1wF4wD5zx2Cfuf0Lyee1Ng1jjDEBszsNY4wxAbOkYYwxJmCWNOogIuNFZKWIrBGRn3sdT6iJSKaIfCIiy0VkqYjc6XVMLUVEwkXkWxF5x+tYWoKIJIrIqyKywv3vPdbrmEJNRH7s/rteIiIviki01zEFm4g8JSIlIrLE51iyiHwoIqvdx6RgXMuSRi0iEg48CpwNDAQmishAb6MKuUrgP1V1ADAGuLUDfOYadwLLvQ6iBT0MfKCq/YFhtPPPLiI9gDuAfFUdjLPI2xXeRhUSTwPjax37OfCxquYBH7v7zWZJ42ijgTWqWqiqh4CXgAs9jimkVHWzqs53n+/G+UPSw9uoQk9EegLnAk96HUtLEJF44ETgHwCqekhVd3gbVYuIADqLSAQQA2zyOJ6gU9VZQHmtwxcCz7jPnwF+EIxrWdI4Wg+gyGe/mA7wB7SGiOQAw4FvvI2kRTwE3A1Uex1IC+kNlAL/dKvknhSRWK+DCiVV3Qj8EdgAbAZ2quoMb6NqMRmquhmcH4ZAejBOaknjaFLHsQ7RL1lEugCvAXep6i6v4wklETkPKFHVAq9jaUERwAjgcVUdDuwlSFUWrZVbj38h0AvoDsSKyA+9japts6RxtGIg02e/J+3wdrY2EYnESRjPq+rrXsfTAo4HLhCRdThVkKeKyL+8DSnkioFiVa25i3wVJ4m0Z6cDa1W1VFUrgNeB4zyOqaVsFZFuAO5jSTBOaknjaHOBPBHpJSJROI1mb3scU0iJiODUcy9X1f/zOp6WoKr3qmpPVc3B+W88U1Xb9S9QVd0CFIlIP/fQacAyD0NqCRuAMSIS4/47P4123vjv423gGvf5NcBbwThpRDBO0p6oaqWI3AZMx+lp8ZSqLvU4rFA7HvgPYLGILHCP3aeq73kYkwmN24Hn3R9EhcB1HscTUqr6jYi8CszH6SX4Le1wOhEReRE4GUgVkWLgf4ApwMsicgNO8rw0KNeyaUSMMcYEyqqnjDHGBMyShjHGmIBZ0jDGGBMwSxrGGGMCZknDGGNMwCxpmJATkQkioiLS34NrrxOR1Ba83nsikthAmSdrJoQUkfsCPO99tfa/bHqUjSciJ4tIowfFichwEWnU3F4i8pKI5DX2WqZlWJdbE3Ii8jLQDWfGzftb+NrrcGY4LWvCeyNUtTL4UR1xjT2q2iVY5ZoZS72fV0TuB/ao6h8bec5XgN+q6sJGvOck4Ieq+qPGXMu0DLvTMCHlzmd1PHADPlNSu79cP/VZ2+F5d8Ruzd3Br0RkvogsrrlDEZH7ReSnPudY4k6wiIi8KSIF7roJkwKIa4+I/Mm9xscikuYe/1REfi8inwF3ikiaiLwmInPd7fiazyUi/3TjWyQiF/vEnioiOe7nesZ9/VURifG5Rr6ITMGZfXWBiDxf3+eop9we91FE5EH3u1gsIpc39P3W+h5qf97zReQbcSY0/EhEMtzveDLwYzeGE+r7XmqdOw4YWpMw3P9+z4jIDPd7ukhE/uDG/YE4U9kAzAZOF2dWWtPaqKpttoVsA34I/MN9/iUwwn1+MrATZ26vMOArYJz72jrgdvf5LcCT7vP7gZ/6nHsJkOM+T3YfO7vHU3zOlVpHXApc5T7/b+AR9/mnwGM+5V7wiSsLZ6oVgP8FHvIpl+R7PSDHvcbx7vGnamJ3r5HvPt9TK676Pkftcnvcx4uBD3FmL8jAGfnbzd/3W+s8tT9vEt/XQNwI/Kme777O76XWuU8BXvPZvx/4HIjEWctjH3C2+9obwA98yn4IjPT6369tR2+WyU2oTcSZghyciQEn4kzpADBHVYsBxJm+JAfnjwo4E8sBFAAXBXCdO0Rkgvs8E8gDtvkpXw1Mc5//y+d6+BwHZ8K7gT4/0uPdX9Cn43PnpKrb67hGkap+4XONO3Cm6Q7m5xgHvKiqVTgT1H0GjAJ24f/79eX7eXsC08SZ4C4KWFvPdev8XtRZj6VGN5yp2H29r6oVIrIYJ9F94B5f7MZXowRnVtqONAtxm2BJw4SMiKQApwKDRURx/kioiNztFjnoU7yKI/89HqzjeCVHVqlGu9c5GeeP2FhV3Scin9a81gi+jXt7fZ6Huefd71vYreppqEGw9ut+yzfxc9Q1lX8Nf9+vL9/P+1fg/1T1bTee++t5T53fSy37OTr+gwCqWi0iFereVuAkcd/4ot33m1bG2jRMKF0CPKuq2aqao6qZOL9cxzXxfOtwp/IWkRE4ayQAJADb3T+0/XGWrG1ImBsfwJXU/QscYAZwW82OiBxTz/G61l/Oku/X4J5YzzUqfOry/X0O33K+ZgGXi7PWeRrOynxz6vksgUgANrrPr/E5vhuI89mv73vxtRzo08Q4+gLtfaLQNsmShgmliTh11b5ew/kj3RSvAcluVcvNwCr3+AdAhIgsAn4DfB3AufYCg0SkAOdu6Nf1lLsDyHcbs5fhNAgD/BZIchugF+LU39e2HLjGjSsZeLyOMlOBRW4Dt7/P4VvO1xvAImAhMBO4W50p0JvqfuAVEZkN+PY4+zcwoaYhnPq/l8NUdQWQ4FbnBUxEMoD96q46Z1oX63JrOiQJcRdWt8fRO6o6OFTXaAtE5MfAblUNeKyG+55dqvqP0EVmmsruNIwxofQ4R7atBGIH8EwIYjFBYHcaxhhjAmZ3GsYYYwJmScMYY0zALGkYY4wJmCUNY4wxAbOkYYwxJmD/DxxJplU3smc3AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "precip = np.linspace(0,10,30)\n",
    "frac_remain = []\n",
    "Stabilization.InputData.Stabil['is_stay_inplace']['amount'] = True\n",
    "for i in precip:\n",
    "    Stabilization.InputData.Precip['ann_precip']['amount'] = i\n",
    "    Stabilization.calc(Inc_flow=ContaminatedSoil)\n",
    "    frac_remain.append(sum(Stabilization.report()['Remaining in Soil'])/sum(ContaminatedSoil.PFAS))\n",
    "\n",
    "plt.plot(precip,frac_remain)\n",
    "plt.xlabel('Annual precipitation rate (m)')\n",
    "plt.ylabel('Percent of Incoming PFAS that \\n remains in the Soil (%)')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in Stabilization (Stabilized soil is transported to landfill)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "scrolled": true
   },
   "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>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>100.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Volatilized  Remaining in Soil  Leachate  RunOff\n",
       "PFOA   0.0          100.0              0.0       0.0   \n",
       "PFOS   0.0          100.0              0.0       0.0   \n",
       "PFBA   0.0          100.0              0.0       0.0   \n",
       "PFPeA  0.0          100.0              0.0       0.0   \n",
       "PFHxA  0.0          100.0              0.0       0.0   \n",
       "PFHpA  0.0          100.0              0.0       0.0   \n",
       "PFNA   0.0          100.0              0.0       0.0   \n",
       "PFDA   0.0          100.0              0.0       0.0   \n",
       "PFBS   0.0          100.0              0.0       0.0   \n",
       "PFHxS  0.0          100.0              0.0       0.0   "
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Stabilization = ps.Stab()\n",
    "Stabilization.InputData.Stabil['is_stay_inplace']['amount'] = False\n",
    "Stabilization.calc(Inc_flow=ContaminatedSoil)\n",
    "Stabilization.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deVhU9eIG8HcYBIFBSlESMVBRSZhhEbziipa4AZUrZKahXdc00xYtI1tckp96UYs2c7kKdNMUTZA0Fbcb67gRqCmyaCUuCMgyDOf3BzU3FBWVme8A7+d5eJ7mcM73vDPB6+E755yRSZIEIiIyPBPRAYiImioWMBGRICxgIiJBWMBERIKwgImIBGEBExEJYvogK9va2kpOTk56ikJE1DilpqYWSJLU+vblD1TATk5OSElJqb9URERNgEwmu1jbck5BEBEJwgImIhKEBUxEJAgLmIhIEBYwEZEgLGAiIkFYwEREgrCAiYgEYQETEQnCAiYiEoQFTEQkCAuYiEgQFjARkSAsYCIiQVjARESCsICpVgqFQm9jOzk5oaCgoM7rq9Vq7N69W295iERhAZPRYwFTY8UCpjq7cuUKRo4cCR8fH/j4+ODIkSMAgKSkJPTq1Quenp7o1asXsrKyAABarRbz5s2DUqmESqXC6tWrdWOtXr0aXl5eUCqVyMzMvOs4FRUVeO+99xATEwMPDw/ExMSgpKQEoaGh8PHxgaenJ3bs2GH4F4OoPkiSVOev7t27S9Q0WFlZ3bEsJCREOnTokCRJknTx4kXJxcVFkiRJKiwslDQajSRJkvTjjz9KI0aMkCRJkj799FNpxIgRuu9dvXpVkiRJcnR0lCIiIiRJkqS1a9dKkyZNuuc433zzjTRjxgxdjvnz50ubNm2SJEmSrl+/LnXu3FkqLi6ux2dPVL8ApEi1dOoDfSYcNW179+5FRkaG7vHNmzdRVFSEwsJCTJgwAWfPnoVMJoNGo9GtP3XqVJiaVv+YtWzZUrftiBEjAADdu3fHtm3bAOCu49wuISEBsbGxCA8PBwCUlZUhJycHTz31VP0/aSI9YgFTnVVVVeHYsWOwsLCosfzVV1/FgAED8P333yM7Oxt+fn4Aqv+6kslktY5lbm4OAJDL5aisrAQALFy4sNZxbidJErZu3YquXbvWzxMjEoRzwFRn/v7+WLNmje6xWq0GUH3k2q5dOwDA+vXra6wfGRmpK9hr167dc/y7jWNtbY2ioiLd48GDB2P16tWo/ssOSE9Pf/gnRSQQC5hqdevWLTg4OOi+VqxYgYiICKSkpEClUqFbt26IjIwEALz55puYP38+evfuDa1Wqxtj8uTJePLJJ6FSqeDu7o4tW7bcc593G2fAgAHIyMjQvQm3cOFCaDQaqFQquLm5YeHChfp5EYj0TPbXUURdeHt7SykpKXqMQ0TU+MhkslRJkrxvX84jYCIiQVjARESCsIDJIOLj49G1a1c4Oztj6dKlouMQGQUWMOmdVqvFjBkzEBcXh4yMDERFRdU4n5ioqWIBk94lJSXB2dkZHTt2hJmZGYKDg3n5MBFYwGQA+fn5aN++ve6xg4MD8vPzBSYiMg4sYNK72k51vNsVckRNCQuY9M7BwQG5ubm6x3l5ebC3txeYiMg4sIBJ73x8fHD27FlcuHABFRUViI6ORlBQkOhYRMLxZjykd6amplizZg0GDx4MrVaL0NBQuLq6io5FJBwvRSYi0jNeikxEZGRYwEREgrCAiYgEYQETEQnCAiaDSU1NRUVFhegYREaDBUx6V1lZiX9Omw7f3n3Rd8Az9/1oIqKmggVMerd+wwZ8GfkZNOWlSDp6CLNfnyc6EpFRYAGT3r0QEoJTp04BAL7++musDP9EcCIi48ACJr2ztLTUXfnWpUsX2NraCk5EZBxYwEREgrCAiYgEYQETEQnCAiYiEoQFTEQkCAuYiEgQFjARkSAsYCIiQVjARESCsICJiARhARMRCcICJiIShAVMBiWTyURHIDIaLGAyCEmSAABt27YVnITIeLCAySCWL18OAHBychIbhMiIsIDJIN566y04OjrCxIQ/ckR/4W8D6d3AgQMBAAcOHBAbRKDi4mJkZWWJjkFGhgVMerNv3z7IZDLs378f3333XZOefli1ahU8vX2QmZkpOgoZEdlfb47Uhbe3t5SSkqLHOE2HVqvFqOBxuFFY+MhjHfgxvh4S3ZvfoCF1Xvf2PK3tnoCryqO+IzUo2ed/xaWiStg2l3AiLQWtWrUSHYkMSCaTpUqS5H37clMRYQjQaDTIysrCxRITWHkMf6SxLDpdRfnlM/WU7E6Wzv/AaWufOq/fooclynJPwdpzGEyaWwMATusrXEPh7gO7J5W49fO3GBLwLI4c/AlmZmaiU5FgPAIWKD8/H+7dfWDi+zIsu/YSHYcMQKrSomjXMgz17ox/b/iG50U3EXc7AuYcsEDt2rXDnh924tb+z1D+2znRccgAZCZyKIbMwQ/7j2DZ8nDRcUgwFrBg3bt3x7ovPkfxriWoLLoqOg4ZgImZBRSBC/DR0uWIjY0VHYcEYgEbgTFjRmPurBko3rUEVZoy0XHIAExbtIF1wNsY99JEnDhxQnQcEoQFbCTCFr6LQb6eKEmIgCRViY5DBmBu3xXN+03GoKHD8fvvv4uOQwKwgI2ETCbDpvXr0N68DCVHo0THIQOx6tYf2k794D8sEGVl/OunqWEBG5HmzZtjzw+xkJ8/hJKMA6LjkIFY+QYjT2OJF16aiAc5K4kaPhawkbGzs8Pe+N0oO7QO5fm8aqopkMlMYDVoFn76+QTe/+DDu64nl8vh4eEBNzc3BAYG4saNG3rNlZKSglmzZt13vV69eArlw2IBGyGlUoktG9ej6IelqLz5h+g4ZAAmzcyhCJyPFas/RUzMt7WuY2FhAbVajVOnTqFly5ZYu3atXjN5e3sjIiLivusdPXpUrzkaMxawkQoMDMR7899CcexiVJXfEh2HDMBU0RLWgQswacpU3O+CJ19fX+Tn5+seL1++HD4+PlCpVAgLCwMAZGdnw8XFBZMnT4abmxvGjRuHvXv3onfv3ujcuTOSkpIAAElJSejVqxc8PT3Rq1cv3U2DDhw4gICAAADA+++/j9DQUPj5+aFjx441ilmhUOjW9/Pzw6hRo+Di4oJx48bpplR2794NFxcX9OnTB7NmzdKN29SxgI3YG/PmIvCZvijesxJSlVZ0HDIAM7tOsBw4HUOGByIvL6/WdbRaLfbt24egoCAAQEJCAs6ePYukpCSo1WqkpqYiMTERAHDu3DnMnj0bJ06cQGZmJrZs2YLDhw8jPDwcixcvBgC4uLggMTER6enp+OCDD7BgwYJa95uZmYk9e/YgKSkJixYtgkajuWOd9PR0rFq1ChkZGTh//jyOHDmCsrIyTJkyBXFxcTh8+DCuXLlSHy9Vo8ACNmIymQxffxEJ58fkKDm8SXQcMhDLLr6A61A8M2Q4SkpKdMtLS0vh4eGBVq1a4dq1axg0aBCA6gJOSEiAp6cnvLy8kJmZibNnzwIAOnToAKVSCRMTE7i6uuLpp5+GTCaDUqlEdnY2AKCwsBCjR4+Gm5sb5syZg9Ona79zx/Dhw2Fubg5bW1u0adOm1lPnevToAQcHB5iYmMDDwwPZ2dnIzMxEx44d0aFDBwBASEhIfb5cDRoL2MiZmZkhbud2NM9PRcnJH0XHIQOx8hmBAjM7jA4eh6qq6vPC/5oDvnjxIioqKnRzwJIkYf78+VCr1VCr1Th37hwmTZoEADA3N9eNaWJiontsYmKCyspKAMDChQsxYMAAnDp1Cjt37rzr6XB/H0sul+u2v986PLPj7ljADUCrVq2wd89ulB/7N8pyToqOQwYgk8lgNXA6jp0+j7fmv1PjezY2NoiIiEB4eDg0Gg0GDx6MdevWobi4GED1TZ7++KPub94WFhaiXbt2AID169fX23P4i4uLC86fP6874o6Jian3fTRULOAGwsXFBVtjolAcFw7N9cui45AByEybQRHwNiK/2YgNGzbW+J6npyfc3d0RHR0Nf39/vPDCC/D19YVSqcSoUaNQVFRU5/28+eabmD9/Pnr37g2ttv7fa7CwsMCnn36KIUOGoE+fPrCzs4ONjU2976ch4u0oG5g1a9finY+Ww2bsMpg0V4iOQwZQceUibm5diD0/xKJ3796i4zyU4uJiKBQKSJKEGTNmoHPnzpgzZ47oWAbD21E2EjNnzEDw88NRHBfOMyOaCLPWjrDyn42AZ5/X/Rnf0Hz55Zfw8PCAq6srCgsLMWXKFNGRjAKPgBugyspKDPQfioxiSygG/lN0HDKQ4tSdaHHxINQpP6NFixai49AD4BFwI2Jqaoqd338H62uZKE7/QXQcMhArrwDctOmE50aO0ctcLRkePxOuFgqFQveOsiG999576NevH5555pn7rmtjY4O9e3aje4+eKH3MHhYdPGt8P++zUJiYWQAyGeRWj8N2+OuQKx7/33KT6n97Ww6ajuYOTwEAbiZvx/WDG9D+1X/DxNwKAFClKcPV+NXQ/JENQIKJuQJtxiyqHoMMSiaTwcrvFah3fIBZr83F2tWr6mXcXr168XJiQVjARuSDDz54oPU7deqEndu3YWjgczAd/TGatWpf4/t2IYsht7TB9YMbUPjfb9HymSk1lt+u5JeDMG/bGbfOHINCWf2PQFFKLOSWj6H1pOpzTjVX8yAz4Y+NKDK5KayGvYlN0W9C6fYUptbDXCrLVxxOQdzDva5tT05ORq9eveDu7o4ePXqgqKgIZWVlePnll6FUKuHp6Yn9+/cDqD638rnnnkNgYCA6dOiANWvWYMWKFfD09ETPnj1x7do1AMDEiRPx3XffAQCcnJwQFhYGLy8vKJVKZGZW3xntypUrGDRoELy8vDBlyhS8+OKL+GjReyiK/Rja0pu1Po/m7d2guX7pns9Vc/0ypIoyPNZ3PEoyDuqWa0uuw9T6fx+h3qyVA2SmzR7yFaX6IG+ugHXQO5j31gLdz9ijqMu9HEg/WMD3Udu17RUVFRg7diz+9a9/4fjx49i7dy8sLCx0VyadPHkSUVFRmDBhgu6qolOnTmHLli1ISkrCO++8A0tLS6Snp8PX1xcbN26sdd+2trZIS0vDtGnTEB5e/QGOixYtwsCBA5GWlobnn38eOTk5eHHcOLw8bgyKf1gGSXvn9fmlvybBrLWT7vHvUQtw6ZtXcXnj67plJb8chOVT/WDe3hWaa/nQllTf6lChHITCn7fi8qa5uJ64CZpr+bcPTwI0a9kOiqHz8NzI0Thz5ky9jVvbzzvpDwv4Pmq7tj0rKwtt27aFj48PAKBFixYwNTXF4cOHMX78eADVF044OjrqfjkGDBgAa2trtG7dGjY2NggMDASAGtfk327EiBEAqj+48691Dh8+jODgYADAkCFD8PjjjwMAViz/BD6dHVDy0+e6o5a/iraq/BZseo7WjWsXshj2L69G25dW6Jbd+iURVk/1g0xmAssuvriVdRgAYGbXEe2mfAWbHiNRVVaEyxtfh6Yg95FfV3p0zR1VaPaPEDwzeBiuX79eL2PW9vNO+sPJvPu427XtMpnsjnXv9edaXa7Jv9s2f7/u/m77kMvl2PafaHj18MX11B0A7j7Xe7uKPy5Ac/0Sfo9ZWL2gqhKmNk/A2qv6loEmZhaw7NoLll17QSaTofR8MprZtr/HiGQoCvchKL6ej2FBzyPxpx/RrNmjTQ/V5X4PVH94BPwQXFxccOnSJSQnJwMAioqKUFlZiX79+mHz5s0AgDNnziAnJwddu3at13336dMH335bfcPuhISEGkc+CoUCe+N/QNWJnZAqy+s8ZskvB/FY7xfgMG1d9deMjagsvorKwj9QlpcBbVn1GSGSVoOKglzIW7Sp1+dEj8aq30ScKSjHP6dO55xtA8MCfghmZmaIiYnBq6++Cnd3dwwaNAhlZWWYPn06tFotlEolxo4di/Xr19c4oqgPYWFhSEhIgJeXF+Li4tC2bVtYW1vrvv/kk09i984dkEqLoCnIqdOYJb8cqr4F4t9Ydu6Jkl8SUXnjMn7f8jYufT0Dl7+ZDfMnnGHZtWFeDttYyUzksBryOrbF/4SVq+7/CRZkPHglXANTXl4OuVwOU1NTHDt2DNOmTYNarb5jvS1bojBl9lw8FvwJ5FaPC0hKhlZZ+DtuRL+FbzdvwLBhw0THob+525VwnANuYHJycjBmzBhUVVXBzMwMX375Za3rvfBCCE5lZCAyajlajF5s4JQkgqmNHRTD30LwuPE4kZ4KJycn0ZHoPljADUznzp2Rnp5ep3UrK7VALW8WUiMmVV+iXNubxGR8OAdsxEJDQ9GmTRu4ubk98LYbN27CZ1+vh9XQN/SQjIyR5sZvKN69HN9GbYajo+M91/3444/h6uoKlUoFDw8P/Pzzz3ddd/LkycjIyABQfYFQQUFBveZuyngEbMQmTpyImTNn4qWXXnqg7Y4ePYppr86GzagPIbd6TE/pyJhUlZegOPYjfPR+GIYMGXLPdY8dO4Zdu3YhLS0N5ubmKCgoQEVFxV3X/+qrr+o7Lv2JR8BGrF+/fmjZsuUDbXPx4kUEPPs8FP6zalz9Ro2XVKVF8e5wjA4YjNdmv3rf9S9fvgxbW1vdGTq2trawt7fHvn374OnpCaVSidDQUJSXV5/K6OfnB775rh8s4EakqKgIA/2HQu7xLCw6+YiOQwZScnAdurW1RuSna+q0vr+/P3Jzc9GlSxdMnz4dBw8eRFlZGSZOnIiYmBicPHkSlZWV+Oyzz/ScnFjAjYRWq8VzI8fghnUHWHoFiY5DBlKijoOiIAO7tm+FqWndZhQVCgVSU1PxxRdfoHXr1hg7diw+//xzdOjQAV26dAEATJgwAYmJifqMTuAccKMx+/V5SDv/O1o8H8Z3wJuI0mw1NMkx2PvzMTz22IPN9cvlcvj5+cHPzw9KpRIbNmzQU0q6Fx4BNwKff/EFNsZshWL4m5DJ+W9qU6C5mofi+BXYvvU/cHZ2fqBts7KycPbsWd1jtVoNOzs7ZGdn49y5cwCATZs2oX///vWame7EAjZiISEh8PX1RVZWFhwcHPD111/fsc6BAwcw9835sA56F3IL61pGocZGW3oTRbEfYeXypfDz83vg7YuLizFhwgR069YNKpUKGRkZWLp0Kb755huMHj0aSqUSJiYmmDp1av2Hpxp4KXIDdu7cOXj/wxfmg16DhZOH6DhkAJJWg5vfL8KLAQOwetWK+29ARoEfytnIXL9+HU/7D0WzHsEs3yZCkiSU7P8C3Tu1xar/Wy46DtUDThg2QBqNBgHPjURJa1co3O990j01HrfSduDxkhx8v/8Y5HK56DhUD1jADdC0Ga/il99LYB30mugoZCC3fk2GVh2LvSlJNW4/Sg0bpyAamFX/Wo3/7EqAYuhcyEx4FNQUVFzJRklCBH6I3X7fezw8ike59wg9HBZwAxIfH493wxZBEfQOTMytRMchA9CW3EBR7GJEromAr6/v/Td4BBMnTkR8fLxe90E1sYAbiIyMDIwJGQfF8DfQ7LEnRMchA5AqK1D8w1JMnzwR48e/qPf9Pcy9R+jRsIAbgIKCAjwzZBjMe09A8/b887ApkCQJxXvXoreqM5Z8/KHoOKQnfBPOyFVUVGBIwLOoaN8DCrenRcchAylJ+g5PSNfw7ZZYmJjwOKmxYgEbMUmSMHHSK7hQJIMiYLzoOGQgJVlHgF8S8GNaCiwtLUXHIT3iP61GbMmyT7D7wDFYDX4NMhn/VzUF5b+dQ+n+SCTs3gV7e3vRcUjP+FttpHbs2IHFn/wfFIELYGJmIToOGUBlUQGKdy7Bhq+/hJeXl8H3X5d7j1D94hSEEVKr1Rg34WVYB70L0xatRcchA6jSlKF45xK88dpMjBw5UkiGqKgoIfttyngEbGR+++03+A8LgEX/V2Bu31V0HDIASapCSUIE/Ht3x8J33xEdhwyIBWxESktL4T8sEFWd/WD1VD/RcchASo5G4cnm5di0fh1vpt/EcArCSEiShHEvTUS+VgFFz2DRcchASjL2w/TCYexJS9F9SCY1HSxgI/He+4uwP/kkrEd+xKOgJqI8/xeUHfoGxw4dRJs2bUTHIQE4BWEEoqNjsGrt57AKmA+TZjwKagoqC/9A0Q/LELVpA29+04SxgAVLTk7G5KnTYB24AKYKXoffFFSV30Lxzo8RtuBtBAQEiI5DArGABcrLy8OQgCBYPj0DZnYdRcchA5CqtCjeswJBg/pj3tzXRcchwTgHLEhZWRmeGTIcZVb2sJKqcOvM0YceqyhtN8ovZ9VjuposO/eEZZe63wqxLOckynJOwtorAHLLFnrLBVR/OnDhf//zUNtKFaW6/5b97WIXk+YKtHz6lUfOdrvmju649d9v0aWlGb7+IpJz/cQCFsXMzAxKpRIONwqB0hOPNNa+i+p6SlW7ktP70bN93e8/vC91FwDg2p41AIDWdk9A5XnH5xHWi32Ju+plnL+XsbaiFKpH/H9yuwvnziDv8Ga0sZRjd1oymjVrVq/jUwMlSVKdv7p37y4R1VViYqIEQAIgRUdHi45zV/n5+ZKbm5sEQBo6dKhe9vHhhx9KltYtpKysLL2MT8YNQIpUS6dyDpj0pm/fvpAkCYMHD0ZwcDDOnz8vOlKt7O3tcfLkSXz33XeIi4vDkiVL6n0fc+bMwfG0VHTp0qXex6aGS1ZdznXj7e0tpaSk6DEONVYymQwODg7Izc0VHeWe+vfvj8TERDzI7wXR/chkslRJku6Yh+MRMBnE8uXLkZeXh6qqKtFR7un7778HAOzbt09wEmoKWMBkEHPnzgUAZGdniw1yH399JlpkZKTgJNQUsIDJIP465ery5cuCk9RNQ8lJDRsLmAyKc6tE/8MCJiIShAVMRCQIC5iISBAWMBGRICxgIiJBWMBERIKwgImIBGEBExEJwgImIhKEBUxEJAgLmIhIEBYwEZEgLGCivzl58iQAIDc312g/wYMaDxYw0Z9yc3OhUqkAADk5Oejb3493byO9YgET/al9+/bYsiUKltY2aN22HXbv2smPjie94sfSE/1NSEgwevb8B6ytrWFrays6DjVyLGCi23To0EF0BGoiOAVBRCQIC5iISBAWMDVKcrkcHh4ecHNzQ2BgIG7cuPFI423fvh0qlQouLi5QKpXYvn277nuZmZnw8PCAp6cnfv31V0REROCpp57CuHHjHvVpUCPHAqZGycLCAmq1GqdOnULLli2xdu3ahx7r+PHjmDdvHnbs2IHMzEzExsZi3rx5OHHiBIDqcn722WeRnp6OTp064dNPP8Xu3buxefPm+no61EixgKnR8/X1RX5+PgDgwIEDCAgI0H1v5syZWL9+PQDAyckJYWFh8PLyglKpRGZmJgAgPDwcCxYs0L0516FDB8yfPx/Lly/H7t27sWrVKnz11VcYMGAApk6divPnzyMoKAgrV6407BOlBocFTI2aVqvFvn37EBQUVKf1bW1tkZaWhmnTpiE8PBwAcPr0aXTv3r3Get7e3jh9+jSGDRuGqVOnYs6cOdi/fz8iIyNhb2+P/fv3Y86cOfX+fKhxYQFTo1RaWgoPDw+0atUK165dw6BBg+q03YgRIwAA3bt3R3Z2NgBAkqQ7LsiobRnRg2IBU6P01xzwxYsXUVFRoZsDNjU1RVVVlW69srKyGtuZm5sDqH4Tr7KyEgDg6uqKlJSUGuulpaWhW7du+nwK1ASwgKlRs7GxQUREBMLDw6HRaODo6IiMjAyUl5ejsLAQ+/btu+8Y8+bNw5IlS3RHxNnZ2Vi8eDHmzp2r5/TU2PFKOGr0PD094e7ujujoaIwfPx5jxoyBSqVC586d4enped/tPTw8sGzZMgQGBkKj0aBZs2b45JNP4OHhYYD01JjJHuRuT97e3tLtf4oR1ZVMJsOhQ4fQp08f0VGIDEomk6VKkuR9+3JOQRARCcICJiIShAVMRCQIC5job+Lj49G1a1c4Oztj6dKlouNQI8cCJvqTVqvFjBkzEBcXh4yMDERFRSEjI0N0LGrEWMBEf0pKSoKzszM6duwIMzMzBAcHY8eOHaJjUSPGAib6U35+Ptq3b6977ODgoLuJD5E+sICJ/lTbOfG83wPpEwuY6E8ODg7Izc3VPc7Ly4O9vb3ARNTYsYCJ/uTj44OzZ8/iwoULqKioQHR0dJ1vY0n0MHgvCKI/mZqaYs2aNRg8eDC0Wi1CQ0Ph6uoqOhY1Yixg0jtJknTzq1VVVUZ9L91hw4Zh2LBhomNQE8EpCNK7devXQy6XAwD69++PiaGTBSciMg4sYNK7USNGoO+AZ2DVpj0cHDtiwdtvio5EZBQ4BUF6Z2Njg58S4hAVFYXhw4ejZcuWoiMRGQUWMBmEqakpxo8fLzoGkVHhFAQRkSAsYCIiQVjARESCsIAJAODn54c9e/bUWLZq1SpMnz691vWzs7Ph5uZ2zzGzs7OxZcsW3eOUlBTMmjULALB+/XrMnDkTABAZGYmNGzc+SnwAgJOTEwoKCh55HCJDYQETACAkJATR0dE1lkVHRyMkJOShx7y9gL29vREREXHHelOnTsVLL7300PshaqhYwAQAGDVqFHbt2oXy8nIA1eV56dIl9OnTB2+88Qbc3NygVCoRExNzx7bZ2dno27cvvLy84OXlhaNHjwIA3n77bRw6dAgeHh5YuXIlDhw4gICAgDu2f//99xEeHo5Lly7Bw8ND9yWXy3Hx4kVcuXIFI0eOhI+PD3x8fHDkyBEAwNWrV+Hv7w9PT09MmTKl1ruZERkzFjABAFq1aoUePXogPj4eQPXR79ixY7Ft2zao1WocP34ce/fuxRtvvIHLly/X2LZNmzb48ccfkZaWhpiYGN00w9KlS9G3b1+o1WrMmTPnvhns7e2hVquhVqvxyiuvYOTIkXB0dMTs2bMxZ84cJCcnY+vWrZg8ufpKukWLFqFPnz5IT09HUFAQcnJy6vlVIdIvngdMOn9NQzz77LOIjo7GunXrsGnTJoSEhEAul8POzg79+/dHcnIyVCqVbjuNRoOZM2dCrVZDLpfjzJkzj5TjyJEj+Oqrr3Do0CEAwN69e2t8NNDNmzUHTucAAAL7SURBVDdRVFSExMREbNu2DQAwfPhwPP7444+0XyJDYwGTznPPPYfXX38daWlpKC0thZeXV53eHFu5ciXs7Oxw/PhxVFVVoXnz5g+d4fLly5g0aRJiY2OhUCgAVN/A59ixY7CwsLhjfWO9qQ9RXXAKgnQUCgX8/PwQGhqqe/OtX79+iImJgVarxZUrV5CYmIgePXrU2K6wsBBt27aFiYkJNm3aBK1WCwCwtrZGUVFRnfev0WgwZswYLFu2DF26dNEt9/f3x5o1a3SP1Wq1LtvmzZsBAHFxcbh+/frDPXEiQVjAVENISAiOHz+O4OBgAMDzzz8PlUoFd3d3DBw4EJ988gmeeOKJGttMnz4dGzZsQM+ePXHmzBlYWVkBAFQqFUxNTeHu7o6VK1fed99Hjx5FcnIywsLCdG/EXbp0CREREUhJSYFKpUK3bt0QGRkJAAgLC0NiYiK8vLyQkJCAJ598sp5fDSL9kj3IO8fe3t5SSkqKHuMQETU+MpksVZIk79uX8wiYiEgQFjARkSAsYDKI+Ph4dO3aFc7Ozli6dKnoOERGgQVMeqfVajFjxgzExcUhIyMDUVFRNc7rJWqqWMCkd0lJSXB2dkbHjh1hZmaG4OBg7NixQ3QsIuFYwKR3+fn5aN++ve6xg4MD8vPzBSYiMg4sYNK72k515BVsRCxgMgAHBwfk5ubqHufl5cHe3l5gIiLjwAImvfPx8cHZs2dx4cIFVFRUIDo6GkFBQaJjEQnHm/GQ3pmammLNmjUYPHgwtFotQkND4erqKjoWkXC8FJmISM94KTIRkZFhARMRCcICJiIShAVMRCQIC5iISBAWMBGRICxgIiJBWMBERIKwgImIBGEBExEJwgImIhKEBUxEJAgLmIhIEBYwEZEgLGAiIkFYwEREgrCAiYgEeaBPxJDJZFcAXNRfHCKiRslRkqTWty98oAImIqL6wykIIiJBWMBERIKwgImIBGEBExEJwgImIhKEBUxEJAgLmIhIEBYwEZEgLGAiIkH+H7oD9kHAcjFqAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "Stabilization.plot_sankey()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>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>1050.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Solids flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>900.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>150.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>8.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1000000.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Parameter Unit     Amount\n",
       "0   Mass flow      kg   1050.0   \n",
       "1   Solids flow    kg   900.0    \n",
       "2   Moisture flow  kg   150.0    \n",
       "3   Carbon flow    kg   8.5      \n",
       "4   PFOA           μg   1000000.0\n",
       "5   PFOS           μg   1000000.0\n",
       "6   PFBA           μg   1000000.0\n",
       "7   PFPeA          μg   1000000.0\n",
       "8   PFHxA          μg   1000000.0\n",
       "9   PFHpA          μg   1000000.0\n",
       "10  PFNA           μg   1000000.0\n",
       "11  PFDA           μg   1000000.0\n",
       "12  PFBS           μg   1000000.0\n",
       "13  PFHxS          μg   1000000.0"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Stabilization.Stabilized.report()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "50.0"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Stabilization.Stabilized.Additive_mass"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Additive Log K for PFOA: 3.9 log L/kg\n",
      "Additive Log K for PFOS: 4.3 log L/kg\n",
      "Additive Log K for PFBA: 0.0 log L/kg\n",
      "Additive Log K for PFPeA: 0.0 log L/kg\n",
      "Additive Log K for PFHxA: 0.0 log L/kg\n",
      "Additive Log K for PFHpA: 0.0 log L/kg\n",
      "Additive Log K for PFNA: 0.0 log L/kg\n",
      "Additive Log K for PFDA: 0.0 log L/kg\n",
      "Additive Log K for PFBS: 0.0 log L/kg\n",
      "Additive Log K for PFHxS: 0.0 log L/kg\n"
     ]
    }
   ],
   "source": [
    "for key, val in Stabilization.Stabilized.AddLogPartCoef.items():\n",
    "    print('Additive Log K for {}: {} {}'.format(key, val['amount'], val['unit']))"
   ]
  }
 ],
 "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
}
