{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Stabilization"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import PFAS_SAT_ProcessModels as pspd\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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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                    \n",
       "1   log L/kg OC                    \n",
       "2   log L/kg OC                    \n",
       "3   log L/kg OC                    \n",
       "4   log L/kg OC                    \n",
       "5   log L/kg OC                    \n",
       "6   log L/kg OC                    \n",
       "7   log L/kg OC                    \n",
       "8   log L/kg OC                    \n",
       "9   log L/kg OC                    \n",
       "10  m                              \n",
       "11  kg/m2                          \n",
       "12  kg/m3                          \n",
       "13  kg add./kg soil                \n",
       "14  1:TRUE,0:FALSE                 \n",
       "15  fraction wet weight            \n",
       "16  kg/m3                          \n",
       "17  log L/kg                       \n",
       "18  log L/kg                       \n",
       "19  log L/kg                       \n",
       "20  log L/kg                       \n",
       "21  log L/kg                       \n",
       "22  log L/kg                       \n",
       "23  log L/kg                       \n",
       "24  log L/kg                       \n",
       "25  log L/kg                       \n",
       "26  log L/kg                       \n",
       "27  m/yr                           \n",
       "28  frac                           \n",
       "29  frac                           \n",
       "30  fraction                       "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Stabilization = pspd.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.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   1000     \n",
       "1   Solids flow    kg   850.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": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncominWaste = pspd.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": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAApoElEQVR4nO3deVxUZf//8dcww47iAneamKAkyjpsbqWJuWAqqWlKZir5q7Rssczu29KsNC0rcymzvLPMwNTSsnIjcUlNQFAE98QUzXBnX4br9wff5o5wwWQ4I3yejwePh3POda75nHnI2+OZc12XTimFEEKImmejdQFCCFFXSQALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmjEcCON3dzclKenp4VKEUKI2ik5OfmsUsr979tvKIA9PT1JSkqqvqqEEKIO0Ol0x6+0XW5BCCGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAIsKXFxcrKbvhIQEtm/fbqFqhNCeBLCwWhLAoraTABbXdfToUSIjIwkNDaVz584cOHAAgO+++4727dsTHBxM9+7dOXPmDAC5ubmMGjWKgIAAAgMDWblypbmvSZMmERQURIcOHcztr9RPZmYmCxYs4L333sNoNLJ161ays7N54IEHCA8PJzw8nJ9//rnmPwwhqpNSqso/oaGhStRuzs7OlbZ169ZNHTp0SCml1M6dO1VERIRSSqnz58+rsrIypZRSH3/8sRo/frxSSqkXX3xRPfPMM+bjz58/r5RSClDffvutUkqpCRMmqNdff/2a/UyZMkW9/fbb5n6io6PV1q1blVJKHT9+XLVp06Z6TloICwOS1BUy9YYW5RR1T25uLtu3b2fw4MHmbUVFRQCcPHmSIUOGcPr0aYqLi/Hy8gJg48aNxMXFmds3bNgQADs7O/r27QtAaGgoGzZsuGY/f7dx40YyMjLMry9fvkxubq5F71sLYUkSwOKaysrKaNCgAampqZX2jRs3jvHjxxMVFUVCQgKvvvrqNfuytbVFp9MBoNfrKS0tvaF+ysrK2LlzJw4ODjdzSkJYDbkHLK6pfv36eHl5sXz5cqD8ltWePXsAuHTpEs2aNQPgs88+Mx/To0cP5s+fb3594cKFa77H1fqpV68eOTk55tc9e/Zk7ty55tdX+kdBiFuJBLCoID8/Hw8PD/PPu+++y9KlS1m0aBFBQUH4+fmxevVqAF599VUGDx5MaGgobm5u5j5efvllLly4gL+/P0FBQWzatOma73m1fvr168c333xj/hJuzpw5JCUlERgYiK+vLwsWLLDMhyBEDdGV3x+umrCwMJWUlGTBcoQQovbR6XTJSqmwv2+XK2AhhNCIBLAQQmhEAljcsLVr1+Lj44O3tzczZsyotL+oqIghQ4bg7e1N+/btyczMBGDXrl0YjUaMRiNBQUF88803AJw4cYKIiAh8fX3x8/Pj/fffr8nTEUI7V3o4+Go/MhBDlJaWqpYtW6qjR4+qoqIiFRgYqNLT0yu0mT9/vnr88ceVUkrFxsaqBx98UCmlVF5eniopKVFKKXXq1Cnl7u6uSkpK1KlTp1RycrJSSqnLly+rO++8s1KfQtzKuMpADLkCFjdk165deHt707JlS+zs7Bg6dKj5qYg/rV69mhEjRgAwaNAg4uPjUUrh5OSEwVD+6HlhYaH5meCmTZsSEhIClD961rZtW7KysmrwrITQhgSwuCFZWVk0b97c/NrDw6NSWP61jcFgwNXVlXPnzgHwyy+/4OfnR0BAAAsWLDAH8p8yMzNJSUmhffv2Fj4TIbQnASxqVPv27UlPTycxMZE333yTwsJC877c3FweeOABZs+eTf369TWsUoiaIQEsbkizZs04ceKE+fXJkyfNo9iu1Ka0tJRLly7RuHHjCm3atm2Li4sL+/btA6CkpIQHHniAYcOGMXDgQAufhRDWQQJY3JDw8HAOHz7MsWPHKC4uJi4ujqioqAptoqKizEOKV6xYQbdu3dDpdBw7dsw8/8Px48c5cOAAnp6eKKV49NFHadu2LePHj6/xcxJCKzIZj7ghBoOBefPm0atXL0wmEzExMfj5+TF58mTCwsKIiori0UcfZfjw4Xh7e9OoUSPzzGjbtm1jxowZ2NraYmNjwwcffICbmxvbtm1jyZIlBAQEYDQaAZg+fTr33XefhmcqhOXJUGQhhLAwGYoshBBWRgJYCCE0IgEshBAakQAWQgiNSAALi1u/fj1ff/2N1mUIYXXkMTRhMUop3pz5FtPfegcbgx0/79jB2zNnYGMj/+4LARLAwkKUUgwfGUNcbCz/GvYWNg4uzJn/DOkZB/hu1UpsbW21LlEIzUkAC4vQ6XQ4ODrRslUrLqRvpMS+Ph5Nm1Cvfv1KE/AIUVfJ/wWFxXyyYD7PP/s0qDIoK2PkIw+zPPYL8zSUQtR1EsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFhZVUFCAsjFgU1aMg4OD1uUIYVUkgIVFJWzdjnJtij7/HO7u7lqXI4RVkQAWFvP777+zbt1aHLxCyft1N1FRUVqXJIRVkQAWFmEymRj9xFhc/LtReGQnffr0xc3NTeuyhLAqEsCi2h06dIjefe9ne8Zv2AVEUrR7NS8+/6zWZQlhdWRxLg0kJiUTGxt70/289+471VDN1bVs2Yr7+/evcvv8giJS9uxlX/o+HP174tilLznfTuP1Ka8QHh5uuUJrgby8PEpKSmjQoIHWpYgaJAGsgeTdu5k9ezb12w9CZ+/4j/sxuN5G6aUz1VhZRadMzny6+1yV2+tsDOj/dRcNfB+iMOMnLi7/DzOnv85TTz5psRpri6iBgzl37hw7t22Wp0XqEJ1SqsqNw8LCVFJSkgXLqRuUUoyIGc2aXzKo1/ff6Gz0WpdULY7P7FvhdVBIGC1beWtUza1ly+YEckugX88IvopdKguX1jI6nS5ZKRVWabsEsDZKSkroHNGdIyZ3XO4ZpXU51eL8hg8pK8rH3sMXncFe63JuKTobPY6twrn89WSeGz2MKa+8rHVJohpJAFuh8+fPExgSRpFvP5wDe2pdjrACpTnnuBj3Ip8umMfgwYOu2Eav1xMQEEBpaSleXl4sWbLEoveOR48ezfjx4/H19b1qmwULFuDk5MQjjzxisTpuZRLAVurgwYO063gXjpHP43BHoNblCCtQ9PsRclZNZXP8BkJDQyvtd3FxITc3F4ARI0bQunVrJk2aVNNlihtwtQCWx9A05uPjw4plseT++A4lF05pXY6wAvZNvHGKGEOvPv3Iysq6ZtuOHTua2xw9epTIyEhCQ0Pp3LkzBw4cAGDkyJGMGTOGDh060LJlSxISEoiJiaFt27aMHDnS3NeYMWMICwvDz8+PKVOmmLd37dqVPy+8XFxcmDRpEkFBQXTo0IEzZ8q/BH711VeZNWuWuf3EiRNp164drVu3ZuvWrQDk5+fz4IMP4uvry4ABA2jfvj11/YJOAtgK9OjRg5nTXiP322mUFeZqXY6wAk4+naBtT3r07kt+fv4V25hMJuLj480jDB977DHmzp1LcnIys2bNYuzYsea2Fy5cYMeOHbz33ntERUXx3HPPkZ6eTlpaGqmpqQBMmzaNpKQk9u7dy+bNm9m7d2+l98zLy6NDhw7s2bOHLl268PHHH1+xttLSUnbt2sXs2bOZOnUqAB988AENGzYkIyOD119/neTk5Jv5iGoFCWAr8eTYsUQP6EfuD7NQZSatyxFWwLndIM7YuPHgQw9TVlZm3l5QUIDRaKRJkyacOXOGHj16kJuby/bt2xk8eDBGo5HHH3+c06dPm4/p168fOp2OgIAAbrvtNgICArCxscHPz4/MzEwAvvrqK0JCQggODiY9PZ2MjIxKNdnZ2dG3b/nTLqGhoeZj/27gwIGV2mzbto2hQ4cC4O/vT2Cg3HKTALYi8+a8R2DzhuQlfKJ1KcIK6HQ6XLo/yc97D/PvSa+Ytzs6OpKamsrx48dRSjF//nzKyspo0KABqamp5p/9+/ebj7G3L38qxcbGxvznP1+XlpZy7NgxZs2aRXx8PHv37qVPnz4UFhZWqsnW1tb8iJxer6e0tPSKtf/5HtdqIySArYrBYGD118upd+EQeSnfa12OsAI6gy0ufV7ig08Ws2TJFxX2OTk5MWfOHN555x2cnJzw8vJi+fLlQPmz5nv27Kny+1y+fBlnZ2dcXV05c+YMP/74Y7WeB8Bdd93FV199BUBGRgZpaWnV/h63GglgK+Pq6kr8uh8oTV5BwbHdWpcjrIDeuQH1ov7DmHHPsGPHjgr7goODCQwMJDY2lqVLl7Jo0SKCgoLw8/Nj9erVVX6PoKAggoODadOmDQ899BB33XVXdZ8GY8eOJTs7G19fX15++WX8/PxwdXWt9ve5lchjaFZq69at9O7XH9dB07B1a651OcIKFBxNpCThQ1KSdtGiRQuty7lhJpOJkpISHBwcOHr0KN27d+fgwYPY2dlpXZrFyWNot5jOnTsz9713yPluGqb8S1qXI6yAY6tw9MYo7u11Hzk5OVqXc8Py8/O5++67CQoKYsCAAXzwwQd1InyvRa6Ardz4Fyaw+JsN1Bv4Kjq9rdblCI0ppciL/wBjY1j3w3fo9bVjHpHaTq6Ab1Gz3ppJuE9z8n5awI38YylqJ51Oh3PE4+z+9XeefX5CtfbdqVOnau1PXJ8EsJWzsbFh5VexNCrIIj9pldblCCug0xtwvu9FPotbwcKrDIT4J7Zv315tfYmqkQC+Bbi4uBC/7gfK0taQf/gXrcsRVkDvWI96US8zfsJLJCQkVEufLi4uACQkJNC1a1cGDRpEmzZtGDZsmPzvy0IkgG8RzZs358c135IfP5/iP37VuhxhBWwbNcM5cjz9HxjMkSNHqrXvlJQUZs+eTUZGBr/++is///xztfYvykkA30LatWvHxx/OJ+fb6ZhyL2hdjrACjp5GbNsN5d5e93Hx4sVq67ddu3Z4eHhgY2OD0Wi86pBjcXNkSaK/+etUfzVp8uTJdOnShe7du1+zXXT0UPZlZPDBFzOo/8Dr6AyVH+M5+WEMNnaOoNOhd26IW5/x6F0a/m+7Tfm/u416jMXBoy35h3aQ/c00bh/9IbaNy585VqqMC/EfU3h8D6BDZ7DD7f6J2DZoUu3nLm6Oc1AkuedP0Lf/AyRsXIfBcPO/1n8drizDiS1HAthKvPbaa1Vu+8ZrU9mXsZ9tG+biEjn+isvX3BY9Hb2TKxc2f8alnV/RqPvjFbb/Vd7+Ldh7+JKXsYUGnYcBkL9/K6acczSNmYdOZ0Pp5bPo7GStMmvlfE8MGd9O44mxT/HJwgValyOqSG5BXMW1vohITEykU6dOBAUF0a5dO3JycigsLGTUqFEEBAQQHBzMpk2bAFi8eDH9+/enR48eeHp6Mm/ePN59912Cg4Pp0KED58+fB8rnbF2xYgUAnp6eTJkyhZCQEAICAszzumZnZ9OjRw/8/f1pWM+ZnP0/k7P1s2ueh0Nz/2vOM1xWXEDRyXQa936GvANbzNtNuefRuzRCpyv/K2Ko74beweUffprC0nQ2elx6P8/yNeuZ/f5crcsRVSQBfA1X+iKiuLiYIUOG8P7777Nnzx42btyIo6Mj8+fPR6fTkZaWRmxsLCNGjDDPJrVv3z6+/vprEhMTmTRpEk5OTqSkpNCxY0c+//zzK763m5sbu3fvZsyYMeaJrqdOnUq3bt1IT09n6NChlJQUY/PrNvL2b73qORQc3YWdu6f59ZnY/3Dq03Gc/nw8APmHd+LgFYpto2boHepR9Hv5lzlObTqTf2QXpz4dx/mfPqH4zNHq+EiFBdnYO+MSNYmXX32NtWvX3vDxf95669q1K2vWrDFvnzdvXoWJ20X1kQC+hit9EXHw4EGaNm1KeHg4APXr18dgMLBt2zYefvhhANq0aUOLFi04dOgQABEREdSrVw93d3dcXV3p168fAAEBAf94PtXIyEgaNmzI8tgvKdi8kKLThyoc/2fQlhXl49phsHn7bdHTuX3UXJo+8i4A+fu34Ny2CwBObbuQn7EZKL/ibfb/PqJhlxHodDaciZtEQWbqzXycogbYNmiCS58JPBg9rMJ0lMI6yT3ga6iuLyL+Pv/qX+dmvdn5VAMCAljy6SKGj34C/ZC3MNR3A658r/fvTAU5FB7fS3F2JqADVQboaBARg06nQ2ewxbFVGI6twrBxakDB4Z04ehqrfN5CGw4efpR2Gk6P3n04mSmPLFozuQK+QT4+Ppw+fZrExEQAcnJyKC0tpXPnzixduhSAQ4cO8dtvv+Hj41Ot7/3X+VTXr1/PhQvlj6INGDCAl55/ltzvplNWXHkS7avJP/gzzn4ReIz5FI8x/8Vj7GIMDW6j6GQ6Rb8foTTnHFD+RERJ9jH09f9VrecjLCj7V1q18ta6CnEdEsA3yM7OjmXLljFu3DiCgoLo0aMHhYWFjB07lrKyMgICAhgyZAiLFy+ucOVbHaZMmcL69evx9/dn+fLlNGnShHr16gEw6d8vEdmlHXnrZ1e5v7z9m3Fq3bHCNqfWncjL2ExZ/iWyV77GqUVjOf3fp8BGT/3QvtV5OsJC8lK+p975g3z79fJrtps2bRp+fn4EBgZiNBr55Zerj7IcPXq0eYkiT09Pzp49W60111UyG9otpKioCL1ej8FgYMeOHYwZM8a8oOKf+++6pxuZhua43D1cu0KFZgqO7aZo41ySd+2gVatWV223Y8cOxo8fT0JCAvb29pw9e5bi4mJuv/32676Hp6cnSUlJuLm5VWfptZrMhlYL/Pbbb4SHhxMUFMTTTz9daUVae3t7fvxuFbbHd5CX/pNGVQqtlJw9Qd662Xy3auU1wxfg9OnTuLm5mf+X5ubmxu233058fDzBwcEEBAQQExNDUVERUHFpelF9JIBvIXfeeScpKSns2bOHxMRE85MYf+Xu7s7XX8VxYe1cVGmJBlUKLZgKLpPz3TTmvDuLzp07X7d9z549OXHiBK1bt2bs2LFs3ryZwsJCRo4cybJly0hLS6O0tJQPP/ywBqqvuySAa5nS0lImTpqMa3AkOoNM4F4XKFMJuWtm8ujDQ4iJGVWlY1xcXEhOTmbhwoW4u7szZMgQPvroI7y8vGjdujUAI0aMYMuWLdfpSdwMeQzNSsXExLBmzRr+9a9/sW/fviof9+TTz7Av6yL17h9nweqEtVBKkfvTAtr5NOedt9+6oWP1ej1du3ala9euBAQEMH/+fAtVKa5GroCt1MiRI294NNO8+fOJ++Z7XHq/gM5GlqqpC/KTVtG4IIuVX8ViY1P1X+eDBw9y+PBh8+vU1FRatWpFZmameWrLJUuWcM8991R7zeJ/5ArYSnXp0uWGpgDcsGEDL708BdcHZ2AjczbUCfmHf6EsbQ0bk3aZJ1OvqtzcXMaNG8fFixcxGAx4e3uzcOFCoqOjGTx4MKWlpYSHh/PEE09YqHoB8hiaVcvMzKRv377XvQVx8OBB2nW8C8feE3Bo7l9D1QktFf9xjMsrJxO//kfat2+vdTniOuQxtFrq3Llz3NurN3Ydhkn41hGmvAvkfDedhR/Ol/C9xUkA38KKi4u5L2oAhbeH4BzYU+tyRA1QpcXkrpnB00+M5qGHorUuR9wkCeBblFKKRx97giMXTDjf/YjW5YgaoJQid8NcuoS05Y3XpmpdjqgGEsBWKjo6mo4dO3Lw4EE8PDxYtGhRhf2z3nmX7zZuxTnyOXnioY7I27mMprpLxH3x+RVXQfm7tWvX4uPjg7e3NzNmzKi0/7nnnsNoNGI0GmndujUNGjQAYNOmTebtRqMRBwcHVq1aBZQ/nePl5WXe9+dQ+EuXLtGvXz+CgoLw8/Pj008/ra7TrtXkS7hb0Jo1axj6SAwNhs7EIDOU1Ql5B7bBL5+zd3cSTZs2vW57k8lE69at2bBhAx4eHoSHhxMbG4uvr+8V28+dO5eUlBT++9//Vth+/vx5vL29OXnyJE5OTowcOZK+ffsyaNCgCu2mT5/OpUuXmDlzJtnZ2fj4+PD7779jZ1d5zcK6SL6EqyXS0tKIHj6Cen0mSvjWEUWnD1OQ8BHrf1hTpfAF2LVrF97e3rRs2RI7OzuGDh3K6tWrr9o+NjaW6OjK95RXrFhB7969cXJyuub76XQ6cnJyym+T5ObSqFGjalkctLaTAL6F/PHHH/To3QeHzjHYN2ujdTmiBpRePkvOmjdZ8ukigoODq3xcVlYWzZs3N7/28PAgKyvrim2PHz/OsWPH6NatW6V9cXFxlYJ50qRJBAYG8txzz5kn63nqqafYv38/t99+OwEBAbz//vs3NDCkrpJP6BZRWFhIz/v6Uep1N86+XbUuR9SAsuJCcr+bzkvPP8uAAQMs9j5xcXEMGjQIvb7idwmnT58mLS2NXr16mbe9+eabHDhwgMTERM6fP8/MmTMBWLduHUajkVOnTpGamspTTz3F5cuXLVZzbSEBfAtQSjF8ZAwnihxw7iSPHtUFSpWRt342vTqHM+nfL93w8c2aNePEiRPm1ydPnqRZs2ZXbHulq1yAr776igEDBmBr+79JnZo2bYpOp8Pe3p5Ro0axa9cuAD799FMGDhyITqfD29sbLy8v82re4uokgG8Br78xjY07UnDu+bR5mXhRu+X9vBRPZxOff7qoSk88/F14eDiHDx/m2LFjFBcXExcXR1RUVKV2Bw4c4MKFC3Ts2LHSvivdFz59+jRQflGwatUq/P3LB//ccccdxMfHA3DmzBkOHjxIy5Ytb7juukbuklu55ctX8Pb783Ad8hY2tg5alyNqQF76T9ge38Ha3Un/eFkrg8HAvHnz6NWrFyaTiZiYGPz8/Jg8eTJhYWHmMI6Li2Po0KGVQj4zM5MTJ05Umoxn2LBhZGdno5TCaDSyYMECAF555RVGjhxJQEAASilmzpwpK2ZUgTyGZsWSk5O5594e1Os/BfsmssBiXVB4MoP872ewY9sW/Pz8tC5HVJOrPYYmV8BWKisri159+uEUMUbC9yadnP8Iptzz1dpn48CIav0fibLRY2h7L7nfv8XyL7+Q8K0jJICtUH5+Pj1694W2PXHy6aR1Obc8p9adKDpdPvdt8emD1dLn1MceqNZVr9965z2Oxq7n7Zlv0bt372rrV1g3CWArU1ZWxoMPPcwZvTsu7QZd/wBxXY16XHtOW1PuBfLT4ynZt47mzZry7szpFR69qgk7k1O5++67ee7Zp2v0fYW25B6wlZn470ksXPYd9Qe+Jmu61TBVZqLgyC8UbP2U/n0imff+e7i6umpdlqgFZCjyLWDJki/44JPFuPSZKOGrAZ2NHqfWnWj48Pv8uO8Md7bxJTk5WeuyRC0mAWwlduzYwZhxz1Av6j/onRtoXU6dZmPvhEv3sdBhJBHde7Jt2zatSxK1lASwFTh+/Dh9ovrj3GMcdu6eWpcj/o+TTyccez7Lff3uZ/fu3VqXI2ohCWCN5eTkcG+v+9Ab78exVbjW5Yi/cfQKwf6ex+hz/wByc3O1LkfUMhLAGjKZTAwYNISLLp44hVQeJiqsg3Obuylx9+G55ydoXYqoZSSANfTs8xPY/esZnCMe+0fj/UXNcezyKLErviYhIUHrUkQtIgGskYUff8xncStwvm8COr08jm3t9A4uOHQazsSXp2hdiqhF5DdfA5s3b+bJp8bh2vMpSi+epvTi6X/UT/Efxyg6ub+aq/sfx1ZhGFyrvupGwZFEDA2bYtvYw2I1VQcbOycMDZve8Fp6zm3uZt/Hi9m/fz9t27a1UHWiLpEA1sCJrNO0at0GMjdC5j/v5+C+PdVW05Xk7duIj39Qldv//rd6buTYmqKU4vLlS5zLzcfevydOxvvQO1VtsIVOb4uDf3dmz53PRx/Ms3Cloi6QkXCiWp05c4Z+/fqRmJgIlAeeNdq7dy9vvzubVd+vpX7/KVW+ai85d5Ki717j7O9XXt5HiCuRkXCiRtx2223s2rWLffv2ARASEqJxRVcWGBjIksX/5Z1pU7m86jXKCqv2iJmhUTNyc3PIzs62cIWiLpAAFhbh5+fH559/TkpKCvn5+VqXc1WPPfb/6NG1M3lpG6rUXqfTUd+jtQxRFtVCAlhYzPDhwwGIiYnRuJJre+G5pylNX1/l2yWmBneQlpZm4apEXSABLCyqcePGLFu2TOsyrqljx46Y8i9TVlC1VXxNBkcuyYq/ohpIAAuL8vCw7kfSoPy2gpNLPcqKC6rWXm+gsLDIwlWJukACWFiUjc0t8lfsRkYiyqhFUU1ukd8OIYSofSSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLCo07KysnjvvffIz72x6SWTk5P4+OOPLVSVqCskgEWdlrYvnedfeAHcW2Oo51alY5xad2LbL8m8PGUqBQVVm8JSiCuRABZ1WmSvnvzw/ffYnDtK/v7N121fVlJE4c5Y2vq0Zs/uJBwdHWugSlFbSQCLOi8yMpKdP2/Fdu/X5G79DKXKrtiuNPc8OStf5h6ff/HLz1to0qRJDVcqahsJYCEAX19f9u5OwtOURe6amZQVF1bYX3zmVy7FvchTIx5k5VdxcuUrqoUEsBD/x83NjZ83/0TP4JZcXv5vSi+fBSD/8E4ufzOFj+fN5vWpr6KTFTFENTFoXYAQ1sTe3p4vl3zGmzPfYvrbL2LrfRc2x7bz0/q1tGvXTuvyRC0jV8BC/I1Op+M/L03ki0Uf0fG2MlKTEyV8hUXIFbAQV9G/f3/69++vdRmiFpMrYCGE0IgEsBBCaEQCWAghNCIBLGo1vV6P0WjE39+ffv36cfHixX/c16VLl3jkkUfw9vamVatWPPLII1y6dMm8f8KECfj5+TFhwgSys7Np3749wcHBbN26tRrORNRGEsCiVnN0dCQ1NZV9+/bRqFEj5s+f/4/7evTRR2nZsiVHjhzh6NGjeHl5MXr0aPP+hQsXsnfvXt5++23i4+MJCAggJSWFzp07V8epiFpIAljUGR07diQrKwuArl27kpSUBMDZs2fx9PQEYPHixQwcOJDIyEjuvPNOXnzxRQCOHDlCcnIyr7zyirm/yZMnk5SUxNGjR4mKiiI3N5fQ0FBmzpzJiy++yOrVqzEajTJhj7gqeQxN1Akmk4n4+HgeffTR67ZNTU0lJSUFe3t7fHx8GDduHBkZGRiNRvR6vbndn7c30tPT+fbbb3FxcSE1NRWA2267jaSkJObNm2epUxK1gFwBi1qtoKAAo9FIkyZNOHPmDD169LjuMffeey+urq44ODjg6+vL8ePHa6BSURdJAIta7c97wMePH0cpZb4HbDAYKCsrn/WssLDixDv29vbmP+v1ekpLS/H19SU1NdV8DEBZWRmpqan4+vrWwJmI2kgCWNQJTk5OzJkzh3feeYfS0lI8PT1JTk4GYMWKFdc93tvbm+DgYN544w3ztjfeeIOQkBC8vb0tVreo3SSARZ0RHBxMYGAgsbGxvPDCC3z44YcEBwdz9uzZKh2/aNEiDh06RKtWrWjVqhWHDh1i0aJFFq5a1GY6pVSVG4eFhak/vzkWoipCQkJISUnhRv6eCVHb6HS6ZKVU2N+3yxWwEEJoRAJYCCE0IgEsarW1a9fi4+ODt7c3M2bMqLR/y5YthISEYDAYKn0ZN3HiRPz9/fH392fZsmXm7cOGDcPHxwd/f39iYmIoKSkBYOnSpQQGBhIQEECnTp3Ys2ePZU9O3PIkgEWtZTKZePLJJ/nxxx/JyMggNjaWjIyMCm3uuOMOFi9ezEMPPVRh+/fff8/u3btJTU3ll19+YdasWVy+fBkoD+ADBw6QlpZGQUEBn3zyCQBeXl5s3ryZtLQ0XnnlFR577LGaOVFxy5IAFrXWrl278Pb2pmXLltjZ2TF06FBWr15doY2npyeBgYHY2FT8VcjIyKBLly4YDAacnZ0JDAxk7dq1ANx3333odDp0Oh3t2rXj5MmTAHTq1ImGDRsC0KFDB/N2Ia5GAljUWllZWTRv3tz82sPDwzwXxPUEBQWxdu1a8vPzOXv2LJs2beLEiRMV2pSUlLBkyRIiIyMrHb9o0SJ69+59cycgaj2ZC0JYTPTwUaSkpADQOaI7G35cg4ODg8ZVVU3Pnj1JTEykU6dOuLu707FjxwrzQACMHTuWLl26VJrtbNOmTSxatIht27bVZMniFiRXwMIilFK41nPGRm/AxtYBGxsbTCZTjdbQrFmzCletJ0+epFmzZlU+ftKkSaSmprJhwwaUUrRu3dq8b+rUqWRnZ/Puu+9WOGbv3r2MHj2a1atX07hx45s/CVGryRWwsAidTseCD+YRbAwk87eTvP7qZAyGmv3rFh4ezuHDhzl27BjNmjUjLi6OL7/8skrHmkwmLl68SOPGjdm7dy979+6lZ8+eAHzyySesW7eO+Pj4CveOf/vtNwYOHMiSJUsqhLUQVyMj4USt9sMPP/Dss89iMpmIiYlh0qRJTJ48mbCwMKKiokhMTGTAgAFcuHABBwcHmjRpQnp6OoWFhYSEhABQv359FixYgNFoBMon8mnRogX16tUDYODAgUyePJnRo0ezcuVKWrRoYW4nvy8Crj4STgJYCCEsTIYiCyGElZEAFkIIjUgA13ERERGsW7euwrbZs2czZsyYK7b39PS87vSN06dPr/C6U6dOAGRmZuLv7w9AUlISTz/99D8t22zkyJFVms9XCGskAVzHRUdHExcXV2FbXFwc0dHR/7jPvwfw9u3bK7UJCwtjzpw5//g9hKgNJIDruEGDBvH9999TXFwMlF+lnjp1iqysLAICAvD392fixIlXPLZ///6Ehobi5+fHwoULAXjppZfM67ANGzYMABcXl0rHJiQk0LdvX6B8aK/RaMRoNOLq6spnn32GyWRiwoQJhIeHExgYyEcffQSUP1/81FNP4ePjQ/fu3fnjjz+q/TMRosYopar8ExoaqkTt06dPH7Vq1SqllFJvvvmmGjVqlGrevLn6448/VElJiYqIiFDffPONUkqpFi1aqOzsbKWUUufOnVNKKZWfn6/8/PzU2bNnlVJKOTs7V+j/z9fHjh1Tfn5+SimlNm3apPr06VOhXVJSkgoICFAXL15UH330kXr99deVUkoVFhaq0NBQ9euvv6qVK1eq7t27q9LSUpWVlaVcXV3V8uXLLfCpCFF9gCR1hUyVK2BR4TZEXFwcLVq0oGvXrri7u2MwGBg2bBhbtmypdNycOXMICgqiQ4cOnDhxgsOHD//jGs6ePcvw4cP58ssvcXV1Zf369Xz++ecYjUbat2/PuXPnOHz4MFu2bCE6Ohq9Xs/tt99Ot27d/vF7CqE1CWDB/fffT3x8PLt37yY/P9884OBaEhIS2LhxIzt27GDPnj0EBwdXWl24qkwmE0OHDmXy5MnmL+mUUsydO5fU1FRSU1M5duyYeSSaELWFBLDAxcWFiIgIYmJiiI6Opl27dmzevJmzZ89iMpmIjY3lnnvuqXDMpUuXaNiwIU5OThw4cICdO3ea99na2ponKa+Kl156icDAQIYOHWre1qtXLz788ENzP4cOHSIvL48uXbqwbNkyTCYTp0+fZtOmTTd59kJoR+aCEED5bYgBAwYQFxdH06ZNmTFjBhERESil6NOnD/fff3+F9pGRkSxYsIC2bdvi4+NDhw4dzPsee+wxAgMDCQkJYenSpdd971mzZuHn52e+8n7ttdcYPXo0mZmZhISEoJTC3d2dVatWMWDAAH766Sd8fX2544476NixY7V+DkLUJBmKLIQQFiZDkYUQwspIAAshhEYkgIVFXW9VYiHqMglgYTFVWZVYiLpMAlhYTFVWJRaiLpMAFhZzM6sSC1EXSAALIYRGJICFxdzsqsRC1HYSwMJi/roqcXFxMXFxcURFRWldlhBWQ4YiC4sxGAzMmzePXr16mVcl9vPz07osIayGDEUWQggLk6HIQghhZSSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiM3tCSRTqfLBo5brhwhhKiVWiil3P++8YYCWAghRPWRWxBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoZH/D4PlpPoKtMiZAAAAAElFTkSuQmCC\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": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABM3ElEQVR4nO3deVhV5frG8e8LCCg4a1laoZETCjibqUfNKQfMcmzS1F+lZR0tzXNsPJVm2WSaZlnapJadcp5T07QUkyytHHLIMTEn5un9/bGRg6KICCz25v5cF5fuvddwrw08PHutd61lrLWIiIiIiLgDL6cDiIiIiIjklJpXEREREXEbal5FRERExG2oeRURERERt6HmVURERETchppXEREREXEbal7zgTHGGmOCnc6RE8aY640xMcYYb6eziIg4QTVbxL0UqeY1/Rf+7FeaMSY+0+O7LzJPK2PMgTzMsNoYMyivlnelrLX7rbWB1trUvF62MWa6MSYp/f392xiz3BhTM/2154wxyed9T0ZmmreLMWajMSbWGHPcGPOpMabKBdbxXPofniaXma2XMWa9MSbOGLM6B9PfZYzZl57na2NMuUyvlTPGfJX+2j5jzF2ZXutvjEk9bztjjDHXXmbei64j/fWKxpjPjDGnjDEnjDGfXs7yRQoj1eys3LVmG2N8jTGvGWMOpM+71xjz5mVke8EY87MxJsUY89wlpjXGmHHpOY6n/99kej3cGLM5vf5vNsaEX+Q9OPv1U05z5mQd6a/XN8Z8m778o8aYxy53HUVZkWpe03/hA621gcB+oGum5/THPn+8kv5+VwH+AqZnem125u+JtfYVAGNMD+Az4E2gAhACJALrjDFlz86cXozuA/5O//dy/J2+/JcvNaExJgR4F7gXuBqIA97JNMkkICn9tbuByenznLXhvO0MtNYeusB6nsumKF9qHf8FjgDXA1cB4y+1XSKFnWq2I/KrZv8LaAg0BkoCrYAfLyPXLmAksDAH0z4A3A6EAaFAV+DB9Ky+wFzgE6AsMAOYm/78Wa+ct51hF1pJ+gebVhd4Ptt1GGMqAEtw/V0pDwQDy3KwXZKuSDWvF2OM8TPGvGmMOZT+9Wb6cwHAYuDazHvMjDGNjTEbjDEnjTGHjTETz/vBz+l6W6V/Cn3cGPNX+rLuz/R68fRPqvvS96itM8YUT38twhizLT3DamNMrUzz7TXGjDDGbE3/FDzNGHO1MWaxMeaMMWbF2YJijAkyrj2XPumPV6d/wv0ufdpl6b9oZ5d9X3qe48aYp9PX1fZS22qtjcNV3Opc4j0xwGvAi9baz6y18dbaI8AgIAYYlmnyFsA1wKNAn8v5HlhrV1hrPweyNJEXcDcw31r7rbU2BngauMMYUzL9Z+RO4GlrbYy1dh0wD1ejmycutQ5jTHvgOmCEtfaUtTbZWrslm+Wdc4g0fU/Di5kej0z/WTxkjBl0/vQiTlPNdsua3Qj4ylp7yLrstdZ+lKM33pVnhrV2MXAmB5P3A16z1h6w1h5Mz9c//bVWgA/wprU20Vo7ATBAm5xmyYFLrWM4sNRa+2n662estb9eaEHmAkcSMn8P03/mZhjXEbdf0+t3nh15KKzUvLqMBpoC4bg+qTUGnrLWxgK3AYfO22OWiusXsgJwM3ArMCSX664ElAYqAwOBSeZ/n1THAw2AZkA5XJ8604wx1YGZwD+BisAiYP55xfhOoB1QHdenzsXAv9On98LV8F3MXcD9uPbg+QJPABhjauPa43g3rqbxbO5LMsYEps930aYqXQ1cew+/yPyktTYN+DJ9m87qB8wHPk9/3DXT+u4yxmzNSbYcCAEyDhtZa3fj2gtaPf0rxVq7I9P0P6XPk1cutY6mwO/AjPQ/UJuMMf/IzYqMMR1xFda2uPYGtMp1apH8o5p9Lneo2d8Dw40xQ4wxddOb3szre8cYk/mI1pU4p2Zzbr0MAbZaa22m17eStzX7UutoCvxtXEPX/jLGzDfGXJ/LdT0LBAHVcL3X9+RyOW5FzavL3cB/rLV/WWuPAc+TzZ4za+1ma+331toUa+1eXLv+c9UsAMnp60621i7C9Um1hjHGCxgAPGatPWitTbXWrrfWJgK9gYXW2uXW2mRcBbM4roJ51tvW2qPpnzrXAj9Ya7dYaxOAr4B62WT60Fq7w1obj6sxDE9/vgeuPZDrrLVJwDOAvcgyznrCGHMS1yGfQP736RegV/peiLNf1+L64wJw+ALLOnz2dWNMCaAn8Fn6ezCHTEMH0vcAhF4iW04FAqfOe+4UrkNfgcDpi7x2VtPztnN3Ltaf3TqqAO2BVbj+sL6G6xBVBS5fL1zf/23pe16ey8UyRPKbava5Cn3NBsYC43B97yKBg8aYfmcntNYOsdbm9gPF+c6v2aeAwPSGObt6ftYT523njCtc//nrqIJr58tjuBr/Pbg+3ORGL2CMtfaEtfYAMCGXy3Eral5drgX2ZXq8L/25CzLGVDfGLDDGHDHGnAbG8L9f0Mt13FqbkulxHK4f/AqAP3ChRuecvOmfcP/k3E/URzP9P/4CjwOzyXTkAnnOrvfPTOuNA45nsxyA8dbaMtbaStbaiPS9lmd9nv7a2a9DQHT6a9dcYFnXZHq9O5CCaw8GwKfAbcaYipfIkxsxQKnzniuF6/BVdq+d9f1523nj2RfSf45Opv+xGAWMylQwF+Rg/eD6fu611k5L/4M6C9f36ZZcbOs53+Pz/i9SWKhmn6vQ1+z0Zn6StfYWoAzwEvBB5uETeej8mlkKiEnfE5qTmj3+vO3MaLIzN7VAc2BBpudGXWT9568jHtcQik3pH06eB5oZY0rnYluLZM1W8+pyCLgh0+Pr+d9YyAt9Sp0M/AbcZK0thevQjrnAdFciGkgAbrzAa+fkTf80eR1wMI8znO8wrk+MZ9dbHNdg87z0O3AA117VDOl7Ne4EVqY/1Q9Xgd5vjDmC65BVMVyHz/LaNlyHJs9mqQb4ATvSv3yMMTdlmj4sfZ5LstZ2OVsgcZ089nKmgtklfbJLrWMrWX9OL7V3pXim/2cumOd8j3H9XIkUNqrZOVOYanYG6xoXOwk4AdTO4zxwXs3m3Hq5DQg9b9hCKDmv2WUy1ex1QJdMz509AfhS6zi/Zl9OvQbVbDWv6WYCTxnX5YYq4Dq08kn6a0eB8ud9IiqJ6zBujHFdRmRwXgdK/2T+AfC6cZ1w4G2MudkY44frsFBnY8ytxphiwOO4zuxcn9c5zjMH6GqMaZY+Vus58vgPQPon4ydwfT/uMsb4G2MqAe/j+uT6hjGmMq4xa11wHR4Lx1WcxpHDqw6kv5/+uAbVe6Wvp9hFJv8U13a3MK4TQv4D/Ne6BtnH4jrT/z/GmABjzC1AN+Dj3Gz/heRgHV8BZY0x/dK3qweuYvZdNovtnz5tOK7xrSXTt/9z4H5jTK30oRlP59V2iOQh1eycKRQ1G8AY80/jOvmouDHGJ33IQEkuPaaW9PmLpddsL1wf5v3Nxa91+xGu8bWV04c2PM7/rpqwGtcY6EeN6yS/R9Kf/+YyNz07l1rHh0B347qcVjFcdXadtfb8oQZnlTPG3G6M8TLGDMC157p0+vZ/DvzLGFM2/W/jIxdZhkdR8+ryIq4xOFuBn3FdvuNFAGvtb7gK5R+Zxvg8gWsP3xngPWB2PuV6Ij3PJlyXdhoHeFlrf8c1KPttXJ/2u+K6hExSPuUAwFq7DRgKzML1aS8G16VUEvN4PbNxjV8bhusQ13ZcnzxvsdYeT38tylq7zFp75OwXrrE+ocaYOsaYu40x2X2SvhfXoZvJuK5aEI/rewlkXF+yRabtfghXE/sXroKbeWzWkPR8f+H6WRmcPs9ZN5us13ltdJlvy0XXYa39G4jA9fNyCtfwg27W2uiLLAugBK7v4Xu4/vDfB7SxrrN5J+AaP7sL10kWkMffY5ErpJqdA4WoZoNrOMNruIY4RAMPA3daa/8AMMZMMcZMyWY17+Gq031xnbAXz/+uuNLCGBOTadp3cZ3M+zPwC67La72bnjUJ12W07gNO4hqnfPt534uR59Xr7Grphd6PbNdhrf0G197/hbi+H8Fkf9TwCK6xrceADrhq9BRce2D/g2vP9x5gBa4PLB5fr421l9pbLXJhxnU26klch+L2OBxHcsgYY3F9z3blYNpauIq/33nj/ETEzahmux/juo7sJ9baLDfpucj0g4E+1trcnpDoFrTnVS6LMaarMaZE+uHz8bg+2e51NpXkJWNM9/RDXWVx7Tmar8ZVxD2pZns2Y8w1xphb0ocU1MA1ROIrp3PlNzWvcrm64Tr54BBwE65PeNp971kexHUoazeucVt5Pj5QRAqMarZn88U1JOIMrjG1czn3DpAeScMGRERERMRtaM+riIiIiLgNNa8iIiIi4jZ8LmfiChUq2KCgoHyKIiKSfzZv3hxtrc2PO7AVWqrZIuKusqvZl9W8BgUFERkZmTepREQKkDFm36Wn8iyq2SLirrKr2Ro2ICIiIiJuQ82riIiIiLgNNa8iIiIi4jbUvIqIiIiI21DzKiIiIiJuQ82riIiIiLgNNa8iIiIi4jbUvIqIiIiI21DzKiIiIiJuQ82riIiIiLgNNa8iIiIi4jbUvIqIiIiI21DzKiIiIiJuQ82riIiIiLgNNa9uLDAwsNAse/Xq1axfvz6f0oiIuD/VbJG8oeZV8oQKoYiI+1DNFnem5tXD7N69m44dO9KgQQNatGjBb7/9BsD8+fNp0qQJ9erVo23bthw9ehSAmJgY7r//furWrUtoaChffvllxrJGjx5NWFgYTZs2zZj+QsvZu3cvU6ZM4Y033iA8PJy1a9dy7Ngx7rzzTho1akSjRo347rvvCv7NEBEp5FSzRXLBWpvjrwYNGlgpPAICArI816ZNG7tjxw5rrbXff/+9bd26tbXW2r///tumpaVZa61977337PDhw6211o4cOdI+9thjGfP//fff1lprATtv3jxrrbUjRoywL7zwQrbLefbZZ+2rr76asZy+ffvatWvXWmut3bdvn61Zs2bebLRILgGR9jLqnSd8qWYXLqrZIjmXXc32cbZ1lrwUExPD+vXr6dmzZ8ZziYmJABw4cIDevXtz+PBhkpKSqFq1KgArVqxg1qxZGdOXLVsWAF9fX7p06QJAgwYNWL58ebbLOd+KFSvYvn17xuPTp08TExOTr2O+RETciWq2SO6oefUgaWlplClThqioqCyvDR06lOHDhxMREcHq1at57rnnsl1WsWLFMMYA4O3tTUpKymUtJy0tje+//x5/f/8r2SQREY+lmi2SOxrz6kFKlSpF1apV+eKLLwDXkJCffvoJgFOnTlG5cmUAZsyYkTFPu3btmDRpUsbjEydOZLuOiy2nZMmSnDlzJuNx+/btefvttzMeX6g4i4gUZarZIrmj5tWNxcXFUaVKlYyv119/nU8//ZRp06YRFhZGSEgIc+fOBeC5556jZ8+eNGjQgAoVKmQs46mnnuLEiRPUqVOHsLAwVq1ale06L7acrl278tVXX2UM/p8wYQKRkZGEhoZSu3ZtpkyZkj9vgoiIm1DNFskbxjUmNmcaNmxoIyMj8zGOiEj+MMZsttY2dDpHQVLNFhF3lV3N1p5XEREREXEbal5FRERExG2oefVAS5YsoUaNGgQHB/Pyyy9neT0xMZHevXsTHBxMkyZN2Lt3b8ZrY8eOJTg4mBo1arB06dJLLnPixIkEBwdjjCE6Ojpft0tExNMFBQVRt25dwsPDadjQdcS0f//+VK1alfDwcMLDw2nWrNkF592zZw9NmjQhODiY3r17k5SUBMD06dOpWLFixvzvv/9+gW2PSH5Q8+phUlNTefjhh1m8eDHbt29n5syZ51y7D2DatGmULVuWXbt2MWzYMJ588kkAtm/fzqxZs9i2bRtLlixhyJAhpKamZrvMW265hRUrVnDDDTcU+LaKiHiiVatWERUVRebxyq+++ipRUVFERUVd9LauTz75JMOGDWPXrl2ULVuWadOmZbzWu3fvjPkHDRqU79sgkp/UvHqYjRs3EhwcTLVq1fD19aVPnz4ZZ6+eNXfuXPr16wdAjx49WLlyJdZa5s6dS58+ffDz86Nq1aoEBwezcePGbJdZr149goKCCnozRUQkE2st33zzDT169ACgX79+fP31186GEsknal49zMGDB7nuuusyHlepUoWDBw9edBofHx9Kly7N8ePHLzpvTpYpIiJXzhhD+/btadCgAVOnTs14fsSIERmH/e++++4s8x0/fpwyZcrg4+O699D5dfrLL78kNDSUHj168Oeff+b/hojkI91hS0REpJBYt24dlStX5q+//qJdu3bUrFkTcA0bOLtX9XJ17dqVvn374ufnx7vvvku/fv345ptv8jK2SIHSnlcPU7ly5XM+VR84cCDj7ioXmiYlJYVTp05Rvnz5i86bk2WKiMiVO1tbr7rqKrp3787GjRsvOm2HDh0IDw9n0KBBlC9fnpMnT2bcFjZznS5fvjx+fn4ADBo0iM2bN+fzVojkLzWvHqZRo0bs3LmTPXv2kJSUxKxZs4iIiDhnmoiIiIzbBM6ZM4c2bdpgjCEiIoJZs2aRmJjInj172LlzJ40bN87RMkVE5MrExsZm3LI1NjaWZcuWUadOnYtOv3TpUqKionj//fcxxtC6dWvmzJkDuG4F261bNwAOHz6cMc+8efOoVatWPm6FSP5T8+phfHx8mDhxIh06dKBWrVr06tWLkJAQnnnmGebNmwfAwIEDOX78OMHBwbz++usZl74KCQmhV69e1K5dm44dOzJp0iS8vb0vukyACRMmUKVKFQ4cOEBoaKjOYhURyaWjR4/SvHlzwsLCaNy4MZ07d6Zjx47AuWNew8PDMy6Dldm4ceN4/fXXCQ4O5vjx4wwcOBBw1emQkBDCwsKYMGEC06dPL8jNEslzuj2siBQJuj2siIj70O1hRURERMQjqHkVEREREbeh5lVERERE3IaaVxERERFxG2peRURERMRtqHkVEREREbeh5lVERERE3IaaVxERERFxG2peRURERMRtqHmVDElJSbz33ntER0c7HUVERHJg/vz5/PLLL07HEClQal4FgOPHj9O81a0Mf2YsYQ0a8euvvzodSURELiI1NZWHH/0ndw14iJubt2TevHlORxIpMGpehTNnzhBavyE/bv2FFL/SHIuHBo2b6NO8iEgh1b1nHyZPmkhK4NUkl76ObrffzvSPPnE6lkiB8HE6gDgvMDCQh/5vIAkJCYx56SV69epN8E3BXHPNNU5HExGRC+hxewS1a9zEuJfHAvDkqH8RUqumw6lECoax1uZ44oYNG9rIyMh8jCNOM8awefNm6tev73QUkTxljNlsrW3odI6CpJrt+YwxREREMHfuXKejiOSp7Gq2hg2IiIiIiNtQ8yoiIiIibkPNq4iIiIi4DTWvIiIiIuI21LyKiIiIiNtQ8yoiIiIibkPNq4iIiIi4DTWvIiIiIuI21LyKiIiIiNtQ8yoiIiIibkPNq4iIiIi4DTWvIiIiIuI21LyKiIiIiNtQ8yoiIiIibkPNq4iIiIi4DTWvIiIiIuI21LyKiIiIiNtQ8yoiIiIibkPNq2Th5+fndAQREckhX19fpyOIFCg1r5IhOTkZgKCgIGeDiIhIjlWtWtXpCCIFSs2rZBg+fDgAAQEBDicREZFL2blzJwADBgxwOIlIwVLzKhkmTpzILbfc4nQMERHJgZYtWwJQs2ZNh5OIFCw1rwLAl19+CcDcuXMdTiIiIpdy/Phxjhw5wuuvv+50FJECp+ZViIqKokePHjRv3pzy5cs7HcfjxMXFsW/fPqdjiIiHiI+Pp0KFCgD885//dDaMB0pNTc0YkiGFk5rXIu7tt9+mXr16lCpVirVr1zodxyM9/58XCKvfUA2siFyxTZs2UaJECQDOnDmDMcbhRJ5n8eLF1KxVi+XLlzsdRS7Cx+kAeSk1NRVrbb4tPzExkRMnTuRo2rMF5VL/Xuq18/n6+uLllfvPHCdOnOC3337j6aef5ocffgCgdevWLFiwgLi4uFwvVy4uOTmZOHy5tUMntmz6npIlSzodSaRQSE1NJS0tLd+Wb63Ncc3O7PwafLmPMytWrBjFihW77AxnJSYmsmPHDj777DPeeuutjOdPnDiBj48PCQkJuV62XFhSUhJ+pSrQo3dfNm74jho1ajgdSc7jMc1rWloala+/gb8OH8KY/NmhnJaWmi/LvVy+/sVzPW9SQvw5j728vFm3fgNlymm4QL6xULrjY5w8uJXuPXqzdNF8vL29nU4l4rgeve/i6y8/x8s7f/4UpaWm5MtyL5ePb+6vnZ2SlJh1ecV8qXjV1VcSSS6hXMMupJWpwq0dbuOnzZs0pK6Q8Zjm1cvLi4cHD+a1qR9TqudLeBXzdzqSyDls9ZvZ/NVz/PPxEbz9pk6yEHnskcEsX7GC0j3HUKx8FafjiGQRc/IgnSK6s3bVCt0MohDxqDGvT43+N+1vaUDssglYm3+HokRyw3j7ENj5SWbM/IKp773ndBwRx7Vq1YrXXhnLmfkvkRp/xuk4IlkENL+PXSdSGfjAQ/k6LFEuj0c1r8YYPp7+Adf7JxK7fqbTcUSy8C5ekpIRTzF8xCjWrFnjdBwRxz34wAPc1/tOYhe9gi0kh/lFzjJe3gR0+CfzV6xl/Gs6YlZYeFTzCuDn58fShfPw2bOO2O2rnI4jkkWx8lUI6Dicbnf0YNeuXU7HEXHcW6+Pp361SsSsmqq9W1LoePmVILDrv3l+zMssWLDA6TiCBzavAFdddRUrliwiYe2HJB781ek4IlkUDwqnWKPe3NqhEydPnnQ6joijvL29+WrObMqc2UPcj/OcjiOShU/pqyjZeRR97+3Hzz//7HScIs8jm1eAOnXqMOuTjziz8BVSTv3ldByRLALCbyO2Qm263H4nKSk6XCpFW8mSJVm5dBEpW74mfvcmp+OIZOFXuSb+LQbQ7rbO/PWX+goneWzzCtC5c2eeGz2KM/NeJC1R1zCVwifgHwPYfvgMgx8e6nQUEccFBQWxcN7XxC5/m6Rje52OI5JFQO1WpFRtTvtOXXWNXQd5dPMK8PjwYdzevhUxS1/HFpLrtIqcZby8Cez0BJ/PX8pbEyY6HUfEcc2aNWPy229xZt4YUmNPOh1HJIuAZn35M9Gfe/rdrzHaDvH45tUYw/tTp1CjnC+xaz9yOo5IFl5+AQRGjGb0s8+zdOlSp+OIOO7ee+9h8MB+xCwch01JdjqOyDmM8SKg/aOs+D6K/7z4ktNxiiSPb17BdXu+hfO+ovjhLcRuXeZ0HJEsipWpRGCnJ+jZ5y5+/VUnGYq8POZFbgkNJmbFJO3dkkLHq5g/gV3+xatvTuSLL+Y4HafIKRLNK0C5cuVYuWwxSd9/SsL+rU7HEcnC/7o6+N3Sj7YdOhEdHe10HBFHeXl58flnn3B1WjSxG9UcSOHjU7I8Jbv+i/v/7wE2b97sdJwipcg0rwDVq1fny89ncWbReJJPHHI6jkgWAXVuJfG6RtzW9XaSkpKcjiPiqBIlSrB88QL4dRlxv693Oo5IFn6VginRejAdOnfl4MGDTscpMopU8wrQtm1bXhnzAjHzXiQtIcbpOCJZBDS/lz9OW+4f9IAOl0qRV7lyZZYunE/cqskkHd3tdByRLErUaAa12tPuti7ExenKRgWhyDWvAA8PGULf7hHELBqvKxBIoWOMFwEdhrHwm+94+ZVXnY4j4rgGDRrwwdR3OTN/DClnjjsdRySLgMY9OOpVgV533UNaWprTcTxekWxeASZOeIPQ68oSu/p9p6OIZOHlW5zAiNG8NG48c+fOdTqOiON69erJ448+TMyCsaQl6/qaUrgYYwhs+zDf/bSTf41+2uk4Hq/INq8+Pj7M/e8XlDyxg9gfda9iKXx8SlWkZJdR3NN/AD/99JPTcUQc9+zTT9G2aTixyyZgrfZuSeFifIoR2GUU77w/nY8//sTpOB6tyDavAKVLl+abZYtJ+fFL4v/QmYJS+PhdWwP/loNod1tnjhw54nQcEUcZY/hkxodc55dA7IZZTscRycI7oAwlI/7N4KGPsWHDBqfjeKwi3bwCVK1alflf/5fYZW+RHP2n03FEsgio1ZK04H/odoQigL+/P0sXzsN797fEbl/jdByRLHwrBhHQbiidI25n3759TsfxSEW+eQVo3rw5b7/xGmfmv0Rq3Cmn44hkEXBzXw6mBHLXvf0u+woE3t7ehIeHU6dOHbp27crJkyfzJ2S6QYMGsX379mynmTJlCh99pDveSe5cffXVrFiyiIS100g8+JvTcUSyKH5jI7zDu9Gm/W2cOXPmsuZVzb40czl/CBs2bGgjIyPzMY6zhj8xgulfLafkHc9hvIs5HUfkHGnJiZz58imG9u/Nf557NsfzBQYGEhPjuixcv379qF69OqNHj86vmIWWMWaztbah0zkKkqfX7Pnz59O33wDK9HkFn1JXOR1H5BzWWmJWTqZeecvSRfPx9vbO0Xyq2S7Z1Wztec1k/CvjaFzzemJXTtH1NaXQ8SrmR0CXf/HGxCnMmjU7V8u4+eabMy6kvXv3bjp27EiDBg1o0aIFv/3m2oPVv39/Bg8eTNOmTalWrRqrV69mwIAB1KpVi/79+2csa/DgwTRs2JCQkBCeffZ/zXSrVq042zAFBgYyevRowsLCaNq0KUePHgXgueeeY/z48RnTP/nkkzRu3Jjq1auzdu1aAOLi4ujVqxe1a9eme/fuNGnSBE9uxOTyde3alWf+9SQx88aQlqjra0rhYowhsPUDbN59hH8+PiJXy1DNvjA1r5l4eXkxZ/ZnlE88RFzkV07HEcnCJ7AcJbv+m0EPDWbTpk2XNW9qaiorV64kIiICgAceeIC3336bzZs3M378eIYMGZIx7YkTJ9iwYQNvvPEGERERDBs2jG3btvHzzz8TFRUFwEsvvURkZCRbt25lzZo1bN2a9bbLsbGxNG3alJ9++omWLVvy3nvvXTBbSkoKGzdu5M033+T5558H4J133qFs2bJs376dF154QbdflAsa8cTjdG3bgtilb+q63VLoGG8fAjuPZMasOUy9SP27GNXsi1Pzep7AwEBWLFlI2s8Lidv5g9NxRLLwvboaJdoMoWOXCA4cOHDJ6ePj4wkPD6dSpUocPXqUdu3aERMTw/r16+nZsyfh4eE8+OCDHD58OGOerl27Yoyhbt26XH311dStWxcvLy9CQkLYu3cvAJ9//jn169enXr16bNu27YJjpnx9fenSpQvgutD82XnPd8cdd2SZZt26dfTp0weAOnXqEBoamtO3SIoQYwzTpk7hxjJexK772Ok4Ill4Fy9JyYinGD5iFKtXr77k9KrZl6bm9QKuu+46Fi+YR9zKSST99YfTcUSyKFH9ZkzIbbTt2JnY2Nhspy1evDhRUVHs27cPay2TJk0iLS2NMmXKEBUVlfH166+/Zszj5+cHuI5GnP3/2ccpKSns2bOH8ePHs3LlSrZu3Urnzp0veCWEYsWKYYwBXCchpKSkXDDj2XVkN43Ixfj6+rJ4/tf4H9xM7M/LnY4jkkWxcpUJ6Dic2+/sya5du7KdVjX70tS8XkTjxo15b/IkzswbQ2rMCafjiGQR0OgOon2vpmefu3N0O8ISJUowYcIEXnvtNUqUKEHVqlX54osvANeJBZdzI4TTp08TEBBA6dKlOXr0KIsXL871dlzMLbfcwueffw7A9u3b+fnnn/N8HeI5ypcvz4qli0jc8AkJ+/WzIoVP8aBwijXuw63tb8vRFQRUsy9OzWs2+vbtw2NDHki/HWGi03FEzmGMIaDNEDZs38PIUf/O0Tz16tUjNDSUmTNn8umnnzJt2jTCwsIICQm5rNvQhoWFUa9ePWrWrMldd93FLbfcktvNuKghQ4Zw7NgxateuzVNPPUVISAilS5fO8/WI56hZsyZfzp5JzOLxJJ84fOkZRApYQFhHYiuG0OX2O3O0x1I1+8J0qaxLsNZyR8/efLszmsCOwzN2p4sUFqlxpzg1ayQTXnmJ/v37OR0nz6SmppKcnIy/vz+7d++mbdu2/P777/j6+uZqebpUVtExcdIkRr/4KqX7vIKXX4DTcUTOYdNSOTPvJe5oWY/3p05xOk6eKciarT2vl2CM4bOPZ1DZ+wyxP+Tu8kQi+cm7RGkCI0bzyD+H89133zkdJ8/ExcXRvHlzwsLC6N69O++8806ui6AULY88/DC9b+9EzKJXdQUCKXSMlzeBtz3OFwuW8eZbbzsdJ88UZM3WntccOnLkCKH1G2Ib3UNArRZOxxHJIv6PzSR9M5EtkRsJCgpyOk6hoz2vRUtKSgpt2t/G9pgSBLZ5wOk4IlkknzzCqdmjmDPzYzp27Oh0nEJHe17zQKVKlVi2aAHxa6aSeHiH03FEsiherQHe9e6gTfvbOH36tNNxRBzl4+PD/K/mUPLv34jZstDpOCJZFCtTicDOI+jV9+5zrhwgl1bge14z3/asID3zzDO0bNmStm3bXtFyvv76a+648068y1yL8fbBO6AsFToPxzuwLAcmD8DLtzh4uT4TlGs3BP8qtUg6to+/V7xLasxxsGkEhLShdLM+54yfPfThUIqVq0LFbk9eUT4p2qy1xH4zhdAyKSxfsjDHtyPMjWbNmrF+/fp8W35e057X3HPnur17924aNG7KmYQUvEuUBmOuuG6nxp7g+OIJpJw+BmmpeJe+iqt7Pp9Xmy1FTOwvK/Hd+hU//biJChUq5Nt6PKlmF5nmNS+VLVuW1BLlKdN7HKc2zMYmJ1Cu7YMcmDyAa/q94SqQ6dKSEzn8wcOUaz+E4lXrk5acwLGvxlIiuBEl67suBJwc/SfH5r5MWsIZrv2/qXj5+ju1aeIBbGoKZ+b+hz7tm/POxLecjlNoqHnNPXev22vXruUf/2jFVXePw79yLU6smXFFdfv4kokUq3AdpRp2AyDprz34XlXVqc0TDxCzdgZV0w7x3ZpvNLY/XaEcNrB69WpatWpFjx49qFmzJnfffTdnG+lNmzbRrFkzwsLCaNy4MWfOnCEhIYH777+funXrUq9ePVatWgXA9OnTuf3222nXrh1BQUFMnDiR119/nXr16tG0aVP+/vtvwHXv3zlz5gAQFBTEs88+S/369albt27G/YGPHTtGu3btCAkJYdCgQdxwww1ER0dnyV66dGnaNWtA7LI38atSm+QThy66nbHb1+BXuRbFq9YHwKuYP+XaPcSp7+f8b5pf1xBQpzX+QfWJ3/V9Hry7UpQZbx8COo3kkzlfM3lK/p3JGhgYCGT/uyyexV3rdosWLShXriwxS14nNf40/tfVuaK6nRr7N94l/7eHTI2rXKmA5vfyx2noP/D/8q1+elLNdnTM65YtW3jzzTfZvn07f/zxB9999x1JSUn07t2bt956i59++okVK1ZQvHhxJk2ahDGGn3/+mZkzZ9KvX7+Mu0P88ssv/Pe//2XTpk2MHj2aEiVKsGXLFm6++WY++uijC667QoUK/PjjjwwePJjx48cD8Pzzz9OmTRu2bdtGjx492L9//0WzT3jzdaoGpnF63af4VgzKeP7ozH9z6MOhHP5oOADJ0fvwrRR8zrzFyl6DTYonLTEOgNjf1hJQqyUBtVsSu31Nrt9PkbO8/QMpGTGaEaNGZzQM+elCv8vimdy1bgcGBnJXj9uJWTiOuJ3fX1HdLlmvM8cXT+DIzH9xav1sUs4cv9K3VYo4Y7wI6PBPFq3ewNhxr+T7+ty9ZjvavDZu3JgqVarg5eVFeHg4e/fu5ffff+eaa66hUaNGAJQqVQofHx/WrVvHPffcA7guRH3DDTewY4frxKnWrVtTsmRJKlasSOnSpenatSsAdevWzfV9eTt27EjZsmUvmr1jx44kxJzC6/RhfEpflfH81X3HcO39b3PNfa/n6D1IPLwT7+Kl8Cl1Ff43hJF09A9S48/kaF6R7BQrV5nA257g9jt7Zvyu5JcL/S6LZ3Lnur1h3VpS/tpD0v6fKNWkR8bzl1u3i1drQOUH36dkWAeS/z7A4emPkRp3KkfzilyMl29xAruOZsyrr1/WDQhyw91rtqPNa+b7717J/XHPv49v5nv85td9eVetWsUvv/xC5KaNJG+cTcKBbRecrliF60k6cu59jJNPHsH4FsfLrwSxv64h+fgBDkwewMF3/4+0pDjifnevT0BSePnfEEqxJn1p26ETJ07k322O8+p3WQo/d6/bhw8d5Iary5GwbeVFp7tU3QbwLl6SgNqtqNDlcfyuuYmEP3+57Dwi5/MpVYGSXUZxT/8Bl3X718vl7jW70F0qq0aNGhw+fJhNmzYBcObMGVJSUmjRogWffvopADt27GD//v3UqFEjT9ed+b68y5Yty9Ef+9q1azP7s0+IWfgqXOBi2AG1W5F4YDvxe6MA14kAJ1a8S6kmd2JtGnG/reOaAROpMvgDqgz+gKvueJrYX7/N0+2Soi0wrCNxV4fSKaI7ycnJTscRD+ROdTswMJAVSxaStnU+cbt+uOA02dVtgPh9P5GW7Br+kJYYR/KJw/iUqphHWyRFnd811fFv+X+0u60zR44ccTpOoVTomldfX19mz57N0KFDCQsLo127diQkJDBkyBDS0tKoW7cuvXv3Zvr06ed8csgLzz77LMuWLaNOnTp88cUXVKpUiZIlS15yvttuu40Xnn2atPhTpCXFn/OaVzE/Kt75NKfWz+Lgew9y+INH8K10EyXrdyHxz214B5bHp2T5jOn9rgsh+fh+UmL+ztNtk6ItoGV/dkQn8sBDQ9xuYL4Ufu5Wt6+//noWzZ9L3PKJkJZ1j1N2dRsg6cgujswYxqEPHuHIJ09QMqw9ftdUz9PtkqItoFYL0m5qRftOXYmPj7/0DEWM7rCVSWJiIt7e3vj4+LBhwwYGDx5MVFRUjua11jLogYf4am0UJSNGY7zy7/qaIrmRlhjHqc9H8ewTQxk+7J9OxylwulSWZ7qSuv3ZZzN58LHHKdPnFbwDLj5WVsQJ1lpilrzGP6pfxZefzzrn2vBFQaG8VFZhtH//fho1akRYWBiPPvoo7733Xo7nNcYw5Z2J1Lo6gLh1Fz5TVsRJXn4lKBkxmmdeGMOKFSucjiOSJ66kbt91V18eeWAgMfNfxl5g2JeIk4wxBLQdyqpNv/Cfl8Y4HadQ8XE6QGFy0003sWXLllzPX6xYMULr1mHHmqi8CyWSh7z8S+LlF8DBQxe/xqWIO7nSut2gXhjJk9+FtDTQETMpZIyPL16lr2b3H3udjlKoaM9rHnp36lQ+nfM1/m2GOB1FJAublkrMkte4s3M77rv3XqfjiDjuxx9/pP+gB1xDvXyKOR1HJIvYjXOoxEmm6G6J5yjSzeuAAQO46qqrqFOnzhUva9WqVTw+8l+UjHgK7+KXPslLpKDFfjudWhWLM3XypFyNnXrppZcICQkhNDSU8PBwfvjhwmdqAwwaNIjt27cDrjsjXehOdSKXKy9r9qFDh2jfqQvFWz2I33k3JBApDGJ//w7z23KWL15AiRIlLnt+T67ZRXrYQP/+/XnkkUe47777rmg5O3fu5PY7exJw2+MUK1c5j9KJ5J3Yn5ZQ4uhWFvy4iWLFLn8P04YNG1iwYAE//vgjfn5+REdHk5SUdNHp33///SuJK3JBeVWz4+LiaHdbF6jZjoCazfMonUjeSTyyi/hVU/j2mxVce+21lz2/p9fsIr3ntWXLlpQrV+6KlnHixAnaduiEb5O+FL8hLI+SieSd+H0/kbxxFiuXLc727kPZOXz4MBUqVMi4zFGFChW49tprWblyJfXq1aNu3boMGDCAxMREAFq1aoWnn+UuBS8vanZaWhq9776XI6YcAU165lEykbyTciaamPljmTHtPerXr5+rZXh6zS7SzeuVSk5OpsvtdxJ7VV0Cwjo6HUcki+S/DxK7+DW+mvM5N910U66X0759e/7880+qV6/OkCFDWLNmDQkJCfTv35/Zs2fz888/k5KSwuTJk/MwvUje+/dTz7Au6ncC2z5c5C49JIVfWnICMfPHMnLYUO68885cL8fTa7aa11yy1vLQkEf49WgsAS37Ox1HJIvUhBjOzHuJ8ePG0Lp16ytaVmBgIJs3b2bq1KlUrFiR3r178+6771K1alWqV3ddnL1fv358+63uDieF1yeffMqk9z4gsMsojI+v03FEzmFtGrFL36JD84Y8NfrfV7QsT6/ZRXrM65V48623mbNoOaV7jdMNCaTQsakpxC4cx729uvPQgw/myTK9vb1p1aoVrVq1om7dukyaNClPlitSEDZs2MCDDw+ldI8XdEMCKZRi18/khhLJfPThtDw5KuDJNVt7XnNh8eLFPP38CwR2HY2X3+WfASiSn6y1xK6eSr2qVzPhjdfyZJm///47O3fuzHgcFRXFjTfeyN69e9m1axcAH3/8Mf/4xz/yZH0ieWnfvn10jridgHZD8a0Y5HQckSxit6/CZ886li6clye3UPb0ml2km9e+ffty88038/vvv1OlShWmTZt2yXm2b99O77vuIbDzSIqVqVQAKUUuT+yP8yl9eg9fzZmNt3feHBWIiYmhX79+1K5dm9DQULZv387LL7/Mhx9+SM+ePalbty5eXl489NBDebI+kQvJTc0+c+YMbTt2xissghLBjQsgpcjlSTz4KwlrP2TFkkVUrFgxT5bp6TXbWGtzPHFRuE92dqKjowmt35Dk0DsIqHOr03FEsojfHUnSqklsidxIUFCQ03EKlezuk+2pinrNTk1NpWPnCH48lqYTtKRQSjn1FydnP8nsjz+kc+fOTscpVLKr2UV6z+vlSExMpGOXbiRd10SNqxRKScf2Ebt8AgvmfqXGVQQY9sRIIncdJLDNg2pcpdBJS4zjzLwXeW70KDWul0nNaw5Ya+k/8P/YE+NFQPN7nI4jkkVq3CnOzH+Jdya8yS233OJ0HBHHTX3vPaZ/9jmBnZ/EeOvWr1K4nL1d9+0dWvP48GFOx3E7al5zYMzL41j87Q8EtP8nxugtk8LFpiQTs+BlHrr/Pu67794CXXde3q5TJK+sWbOG4SNGUTJiNN7FSzkdRySL2LUfUaOCP++/O7nAjwp4Qt1WJ3YJX331FWPHv0Fgl3/j5evvdByRc1hriVk5iWZ1bmTc2JcKfP39+/dnyZIlBb5ekYvZvXs33e7oQUDH4RQrf53TcUSyiN26jOKHt7Bw7n9zdbvuK+UJdVvNazaioqK49/6BlOwyCp9SFZyOI5JF7Kb/clXKMT6f+QleXgX/65wXt+sUySsnT57k1va3UaxRL4oHhTsdRySLhP1bSfr+U1YuW+xY7fSEuq2bFFzEkSNHaN+pC8X/8QB+11R3Oo5IFnE7NsC2xaz4MZKAgACn44g4KiUlha7dexBTvhaB4Z2cjiOSRfKJQ8Qsfo2vv5iVcZcryR3teb2A+Ph42t3WhbSbWhNQq4XTcUSySDq6m7hv3mHJwvlUrlzZ6Tgijhv88FC2HTpNQKuBTkcRySI1IYaYeS/yypgXaNu2rdNx3J6a1/NYa+l7bz8OpZUioGlvp+OIZJES8zdn5o/hg6nv0rBhkbpsqcgFvTVhIp/PX0Jgpyd0u24pdGxqCrGLXuWuO7oxZPBgp+N4BDWv53n62edYE7mNgLaP6LqAUuikJScSO38sw4cOoVevnk7HEXHcsmXLGP3sc+m369bwGSl8YtdMI/T68rz91utOR/EYal4zmTlzFm+9M5XALv/Cq9iV31tYJC9Zm0bs8gnc2jSM55552uk4QO5u1ymSV3777Td69O5LYKcRFCt7jdNxRLKI/XEBJU/sYO6Xn+PjUzhOM/KEul043slCYOPGjfzf4CGUuuN5vAPLOh1HJIvYDbOoUiyOT2Z8WGiOCsycOdPpCFJEHT9+nFs73IZvs3vxv859r1cpniv+j82k/Pgl32z6gdKlSzsdJ4Mn1G01r8Cff/7JbV27UeLWR/C9qprTcUSyiN2+Bu/d37Lsx0j8/XW9YSnakpKS6NilG4mVGxJYt53TcUSySIreT+yyt1iyYC5Vq1Z1Oo7HKfLNa0xMDG07dsarTidK3NTE6TgiWSQe+p2EtdP4bs0qrr76aqfjiDjKWsv9gx5g9ylLyS4Fe0c5kZxIjTtFzLyXmPjm6zRv3tzpOB6pSI95TUtLo0fvuzjudy0lGnZ3Oo5IFimn/+LMgpf5dMaHhIaGOh1HxHHjXh3PwlXfEdhxmK4sIIWOTU0mZuE4/q/fXfTv38/pOB6rSDevT4x8ko2/7Sfg1ocKzRhCkbPSkuKJmT+Gp/81goiICKfjiDhu3rx5vPjyqwR2/TdevsWdjiNyDmstsSun0KTmDbw67mWn43i0Ijts4MMPp/P+RzMp3ecVjHfB31tYJDs2LZWYJW/QpU1zRj7xhNNxRBy3detW7r6vPyUjnsKn1FVOxxHJIi7yK8onHmLO7PWO3K67KCmSzeu6desYOuxxSvV4Ee8ShecMQJGzYtd9THAZLz54710dFZAi7+jRo7S7rTP+LQfhd20Np+OIZBG38wfSflnEys2bCAwMdDqOxytyzesff/xB19vvIKD9Y/hWuN7pOCJZxP68Ar+DkSz+MRJfX1+n44g4KiEhgfadupJ6Y0sCa//D6TgiWSQd/YO4lRNZtXwpVapUcTpOkVCkmtdTp05xa4dO+NS/k+LVGjgdRySLhD9/IXHDx6xdv47y5cs7HUfEUdZa7r6vPweSSxB4cx+n44hkkRpzgjPzx/D+lMk0atTI6ThFRpEZlJGSkkK3O3pypmx1Aup3cTqOSBbJJw4Ts+hV5sz6jJo1azodR8Rxz/3nBVb+sJWAdo9iTJH5cyVuIi05kZgFY3lsyAP06dPb6ThFSpGpBkMfG87WP/8moNUgp6OIZJGWGEvM/JcY+8LztG/f3uk4Io6bPftzXn/7HQK76nbdUvhYa4ldMZFWDUN44fnnnI5T5BSJ5vWdyZP57L9zCew0QtcFlELHpqUSs2g8vbvdxtBHHnY6jojjIiMjGfjgQ5Ts+m98Ass5HUcki9jvZ1PZ+wyffTxDJ9U6wOPHvK5YsYKR/36a0r3G4uWvMwCl8IldM42Qa0vxztsTnI4i4riDBw/SsXNXSrQZjO/VNzodRySL2F/X4rVzFct+jKR4cV1v2Ake3bz+/vvv3NmrD4G3PU6xstc6HUcki9ioRQRG/8r8zRvx8fHoX0eRS4qNjaVtx84Q0pES1Zs5HUcki8RDvxO/ZirfrVlFpUqVnI5TZHnsX8u///6bth064dv0bvyv1201pfCJ37OF5E2fs3Lj95QpU8bpOCKOSktLo2efu/nLpyKBje50Oo5IFimnj2XcrjssLMzpOEWaR455TU5OpnNEd+KvqUdAqE5+kcIn+fifxC59g3lffcmNN+rQqMiT/xrNhm1/EHjrwxpDKIVOWlICMfPH8K8Rw+nWrZvTcYo8j2terbUMfOAhfv87iYAW9zkdRySL1PjTnJn3Em+MH0fLli2djiPiuBkzPmLKhx8R2GUUxke365bCxdo0Ype+QadWN/PvUU86HUfwwOb1tdffYN7yNQR2GK4rC0ihY1OTiVk4jv539eL/BumybSLr16/n4ceGUbLraN2uWwql2HWfULWkZfq093RUoJDwqDGvCxcu5LmXXqZM73F4+ZVwOo7IOay1xK6aSqObqvDG+FecjiPiuL1799KlW3cC2j2Kb8UbnI4jkkXsLyvx/fMHluh23YWKxzSvaWlp9BswEGu8SFr5Nkn5sI5jO7fkw1Iv31W1Gud63r9+3Zjxf//SFSh1bbW8iHRFMmcqKBWq1SnQ9Xk3uZuUv3ZTNnY//139Pd7eOiog8ujwJ4iJjaPMT3OJ+2luni8/Ke40pw7uzvPlXq6KNRrmet5jv0fm2bLyyvmZCkLF4PACXZ/3Tbdgy11P4nczWPPdWipUqFCg65fseUzz6uXlxcplSzlx4kS+rWP79u1MnTo122nOHlLI6b8Xe+5iBg0aRJUqVS4z+f/MmTOH6dOnk5qaSsKpaBJORTNgwAB69uyZ62VeqSlTphAXF5dny1u+fHmW50qVKpVx56qGDRsW6FjTSZPeYc63n1Ii4RgrIzcSGKjrDYsATHh9PHsffSTflp+YmMizzz57yenOr72X+zg77dq145Zbbsnx9Ofbvn07Y8eO5ciRI4CrcQwKCuK9997L9TKv1Oeff86OHTvybHmbNm264N+Adu3aAXD99ddz330Fdw7L+vXrGfP2+6TFneKLzz6hdu3aBbZuySFrbY6/GjRoYMVz7Nmzx1aqVMkCduTIkU7HyRe///679ff3t4B94403HMkw8slR1qeYr92wYYMj6xcXINJeRr3zhC/VbM8SHx9v7733XgtYf39/p+Pki9OnT9vOnTtbwNarV8+RDPPnz7eAff2NtxxZv7hkV7M97oQtybmgoCAOHz7MmDFjeOWVV5gwwfPu8FS9enXi4+MZPHgww4YNY86cOQWe4f8GDWTpkiU0bdq0wNctIp7D39+fjz76iC1btpCQkEDVqlWdjpTnSpYsyYIFC1i2bBlbtmyhS5euBZ6hWbNmfDZzJv98bGiBr1tyxria25xp2LChjYws+LEukv/uv/9+pk+fTnJyssfe6alR48ZEbtpEWlqazhgtgowxm621zg8YLECq2Z5r3bp1tGjRguXLl9O2bVun4+SLWbNm0bdvX3bv3k21as6fnyEFK7uarT2vAsAHH3wAwMiRIx1Okn/WrV0LwPvvv+9wEhGRK9O8eXPgf+NCPVGfPn0AaN26tcNJpLBR8yqA6wSEJk2a8MYbbzgdJd/4+fkRGBjIAw884HQUEZErtnjxYgBOnjzpbJB8NG7cOPbv309qaqrTUaQQUfMqGebPnw9ASkqKw0nyz6JFi5yOICKSJzp27AjgkecrnDVixAjgf426CKh5lUwqVqwIwP79+x1Okn9CQ0MBz27QRaRoycvLVhU2Z89P8ORtlMun5lWyOH36tNMR8k3x4sUBSE5OdjiJiEjeiImJcTpCvisK2yg5p+ZVRERERNyGmlcRERERcRtqXkVERETEbah5FRERERG3oeZVRERERNyGmlcRERERcRtqXkVERETEbah5FRERERG3oeZVRERERNyGmlcRERERcRtqXkVERETEbah5FRERERG3oeZVRERERNyGmlcRERERcRtqXkVERETEbah5FRERERG3oeZVRERERNyGmlcRERERcRs+TgcQKQhpaWksX76c5ORkAJYvX46fnx/t2rXDy0uf4URECpsNGzZw4sQJAH744QcWLVpEkyZNKF++vMPJxGlqXqVI+OOPP7itUyf8S1UAoE///yPhVDS//for1atXdzidiIhkZq2lZ+++HDl6FIBv1m9i0aLOTP/oI/rde6/D6cRp2uUkRUJwcDCLFi7EhxRKNeiCj01hwfz5alxFRAohYwzrvl3NDUFVKVO3Fd4Gxr3yKvfdc4/T0aQQUPMqRUbHjh35/ru1tAkqzoZ139KpUyenI4mIyEUEBQWxZdP33NmsFjNnTGPkiCcwxjgdSwoBDRuQIqV27dp8Nedzp2OIiEgOlCpVivenvON0DClktOdVRERERNyGmlcRERERcRtqXkVERETEbah5FRERERG3oeZVRERERNyGmlcRERERcRtqXkVERETEbah5FRERERG3oeZVRERE8pW3tzfh4eHUqVOHrl27cvLkyVwv69SpU9x3330EBwdz4403ct9993Hq1KmM10eMGEFISAgjRozg2LFjNGnShHr16rF27do82BIpDNS8ittRERQRcS/FixcnKiqKX375hXLlyjFp0qRcL2vgwIFUq1aNXbt2sXv3bqpWrcqgQYMyXp86dSpbt27l1VdfZeXKldStW5ctW7bQokWLvNgUKQTUvIrbUREUEXFfN998MwcPHgSgVatWREZGAhAdHU1QUBAA06dP54477qBjx47cdNNNjBw5EoBdu3axefNmnn766YzlPfPMM0RGRrJ7924iIiKIiYmhQYMGjBs3jpEjRzJ37lzCw8OJj48v2A2VfOPjdACRK3HzzTezdetWwFUEx48fT8OGDYmOjqZhw4bs3buX6dOnM2/ePOLi4ti9ezfdu3fnlVdeySiCs2fPzljeM888Q3BwMLt372bYsGEZRbBv375MmjSJ+Ph4IiMj2bBhA8WLF3dqs0VE3FJqaiorV65k4MCBl5w2KiqKLVu24OfnR40aNRg6dCjbt28nPDwcb2/vjOnOHo3btm0b8+bNIzAwkKioKACuvvpqIiMjmThxYn5tkjhAzau4LRVBERH3EB8fT3h4OAcPHqRWrVq0a9fukvPceuutlC5dGoDatWuzb9++/I4pbkLDBsTtnC2ClSpV4ujRo5dVBP39/VUERUQK2NnhXvv27cNamzHcy8fHh7S0NAASEhLOmcfPzy/j/97e3qSkpFC7dm2ioqIy5gFIS0sjKiqK2rVrF8CWSGGg5lXcjoqgiIh7KlGiBBMmTOC1114jJSWFoKAgNm/eDMCcOXMuOX9wcDD16tXjxRdfzHjuxRdfpH79+gQHB+dbbilc1LyK21IRFBFxP/Xq1SM0NJSZM2fyxBNPMHnyZOrVq0d0dHSO5p82bRo7duzgxhtv5MYbb2THjh1MmzYtn1NLYaIxr+LWzi+CvXr1YurUqXTu3DlH80+bNo2hQ4dy4403Aq4TwFQERUTyVkxMzDmP58+fn/H/syfdAhk7E/r370///v0znl+wYEHG/8uWLcsnn3ySo3WdvxzxDMZam+OJGzZsaM9e0kI8kzGGLVu2EB4e7nQUkTxljNlsrW3odI6CpJrt+YwxdOvWja+//trpKCJ5KruarWEDIiIiIuI21LyKiIiIiNtQ8ypuacmSJdSoUYPg4GBefvnlLK8nJibSu3dvgoODadKkCXv37s14bezYsQQHB1OjRg2WLl16yWW2aNGC8PBwwsPDufbaa7n99tvzc9NERDxOQdbsgQMHEhYWRmhoKD169Mgy3lY8gLU2x18NGjSw4tkAu2XLFqdjZCslJcVWq1bN7t692yYmJtrQ0FC7bdu2c6aZNGmSffDBB6211s6cOdP26tXLWmvttm3bbGhoqE1ISLB//PGHrVatmk1JScnRMq219o477rAzZszI/42UPAdE2suod57wpZrt+QDbrVs3p2Nkq6Br9qlTpzKWO2zYMDt27NgC2lLJS9nVbO15FbezceNGgoODqVatGr6+vvTp04e5c+eeM83cuXPp168fAD169GDlypVYa5k7dy59+vTBz8+PqlWrEhwczMaNG3O0zNOnT/PNN99oz6uIyGUo6JpdqlQpwLVzLj4+HmNMwW6w5Ds1r+J2Dh48yHXXXZfxuEqVKhw8ePCi0/j4+FC6dGmOHz9+0Xlzssyvv/6aW2+9NaMwiojIpTlRs++//34qVarEb7/9xtChQ/Nr08Qhal5FcmjmzJn07dvX6RgiInIJH374IYcOHaJWrVrMnj3b6TiSx9S8itupXLkyf/75Z8bjAwcOULly5YtOk5KSwqlTpyhfvvxF573UMqOjo9m4cWOOb34gIiIuTtRscN0KvE+fPnz55Zf5sVniIDWv4nYaNWrEzp072bNnD0lJScyaNYuIiIhzpomIiGDGjBmA61axbdq0wRhDREQEs2bNIjExkT179rBz504aN258yWXOmTOHLl264O/vX6DbKiLi7gqyZltr2bVrF+Aa8zpv3jxq1qxZ4Nss+Uu3hxW34+Pjw8SJE+nQoQOpqakMGDCAkJAQnnnmGRo2bEhERAQDBw7k3nvvJTg4mHLlyjFr1iwAQkJC6NWrF7Vr18bHx4dJkybh7e0NcMFlnjVr1ixGjRrlyPaKiLizgqzZaWlp9OvXj9OnT2OtJSwsjMmTJzu5+ZIPdHtYOYduDyueSreHFU+k28OKp9LtYUVERETEI6h5FRERERG3oeZVRERERNyGmlchJSWFjl260eLW9gD0vuseWt7agf379zucTERELuTRx0fSrNWtgOvuVM1a3cp/v557iblEPIOuNiAkJycTGxvLD5E/UqbVAPZvX0VS8j68vPTZRkSkMAos7s/369ZQqmlPUk8eJmrzJsqVKe10LJECoe5EKF68OKuWL+H+e/sSs/4zGtcJJiryB6pUqeJ0NBERuYAxL/6H6R98SNLPSygTd5AfN22kVatWTscSKRC6VJacIyoqijp16uDjo53y4ll0qSzxRHv27KFs2bKUKVPG6SgieSq7mq0ORc6h67uKiLiPqlWrOh1BpMBp2ICIiIiIuA01ryIiIiLiNtS8ioiIiIjbUPMqIiIiIm5DzauIiIiIuA01ryIiIiLiNtS8ioiIiIjbUPMqIiIiIm5Dzasbad26NUuXLj3nuTfffJPBgwdfcPqgoCCio6OzXeaYMWPOedysWTMA9u7dS506dQCIjIzk0UcfzW3sDP3792fOnDlXvBwREXegmi2SP9S8upG+ffsya9asc56bNWsWffv2zfUyzy+E69evzzJNw4YNmTBhQq7XISJSFKlmi+QPNa9upEePHixcuJCkpCTA9Un70KFDHDx4kLp161KnTh2efPLJC857++2306BBA0JCQpg6dSoAo0aNIj4+nvDwcO6++24AAgMDs8y7evVqunTpAkCnTp0IDw8nPDyc0qVLM2PGDFJTUxkxYgSNGjUiNDSUd999FwBrLY888gg1atSgbdu2/PXXX3n+noiIFFaq2SL5xFqb468GDRpYcVbnzp3t119/ba21duzYsfb++++31113nf3rr79scnKybd26tf3qq6+stdbecMMN9tixY9Zaa48fP26ttTYuLs6GhITY6Ohoa621AQEB5yz/7OM9e/bYkJAQa621q1atsp07dz5nusjISFu3bl178uRJ++6779oXXnjBWmttQkKCbdCggf3jjz/sl19+adu2bWtTUlLswYMHbenSpe0XX3yRD++KyKUBkfYy6p0nfKlmO081WyR3sqvZ2vPqZjIfhpo1axY33HADrVq1omLFivj4+HD33Xfz7bffZplvwoQJhIWF0bRpU/7880927tyZ6wzR0dHce++9fPbZZ5QuXZply5bx0UcfER4eTpMmTTh+/Dg7d+7k22+/pW/fvnh7e3PttdfSpk2bXK9TRMQdqWaL5D01r26mW7durFy5kh9//JG4uDjCw8MvOc/q1atZsWIFGzZs4KeffqJevXokJCTkav2pqan06dOHZ555JuPkAGstb7/9NlFRUURFRbFnzx7at2+fq+WLiHgS1WyRvKfm1c0EBgbSunVrBgwYQN++fWncuDFr1qwhOjqa1NRUZs6cyT/+8Y9z5jl16hRly5alRIkS/Pbbb3z//fcZrxUrVozk5OQcr3/UqFGEhobSp0+fjOc6dOjA5MmTM5azY8cOYmNjadmyJbNnzyY1NZXDhw+zatWqK9x6ERH3opotkvd8nA4gl69v3750796dWbNmcc011/Dyyy/TunVrrLV07tyZbt26nTN9x44dmTJlCrVq1aJGjRo0bdo047UHHniA0NBQ6tevz6effnrJdY8fP56QkJCMvQf/+c9/GDRoEHv37qV+/fpYa6lYsSJff/013bt355tvvqF27dpcf/313HzzzXn6PoiIuAPVbJG8ZVxjYnOmYcOGNjIyMh/jiIjkD2PMZmttQ6dzFCTVbBFxV9nVbA0bEBERERG3oeZVRERERNyGmlcRERERcRtqXiXDkiVLqFGjBsHBwbz88stOxxERkWyoZktRpeZVANe1AB9++GEWL17M9u3bmTlzJtu3b3c6loiIXIBqthRlal4FgI0bNxIcHEy1atXw9fWlT58+zJ071+lYIiJyAarZUpSpeRUADh48yHXXXZfxuEqVKhw8eNDBRCIicjGq2VKUqXkVEREREbeh5lUAqFy5Mn/++WfG4wMHDlC5cmUHE4mIyMWoZktRpuZVAGjUqBE7d+5kz549JCUlMWvWLCIiIpyOJSIiF6CaLUWZj9MBpHDw8fFh4sSJdOjQgdTUVAYMGEBISIjTsURE5AJUs6UoM9baHE+s+2SLiLvK7j7Znko1W0TcVXY1W8MGRERERMRtqHkVEREREbeh5lVERERE3IaaVxERERFxG2peRURERMRtqHkVEREREbeh5lVERERE3IaaVxERERFxG2peRURERMRtqHkVEREREbeh5lVERERE3IaaVxERERFxG2peRURERMRtqHkVEREREbeh5lVERERE3Iax1uZ8YmOOAfvyL46ISL65wVpb0ekQBUk1W0Tc2EVr9mU1ryIiIiIiTtKwARERERFxG2peRURERMRtqHkVEREREbeh5lVERERE3IaaVxERERFxG2peRURERMRtqHkVEREREbeh5lVERERE3IaaVxERERFxG/8PUj5UJsd6VKAAAAAASUVORK5CYII=\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": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEGCAYAAACZ0MnKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAA1LklEQVR4nO3dd3xV9f3H8dc7CSEkJKyEnRDAsERERBD3Llpnndhaaa3UWmfXT3/ddml/HWprh6ta66xSixM3TpQlGxHCCnuPICPJ5/fHOcFrTG5ukntzMz7Px+M87rnnfs85nxsxn5zvlJnhnHPOxSIl2QE455xrPjxpOOeci5knDeecczHzpOGccy5mnjScc87FLC3ZASRabm6uFRYWJjsM55xrNmbMmLHJzPKq+6zFJ43CwkKmT5+e7DCcc67ZkLSips+8eso551zMPGk455yLmScN55xzMfOk4ZxzLmaeNJxzzsWsySQNSfdL2iBpXg2fS9KdkpZImiNpRGPH6JxzrV2TSRrAA8DYKJ+fDhSF2wTgr40Qk3POuQhNJmmY2ZvAlihFzgH+aYGpQEdJPRIRS3mF8dDUFbw4b10iLu+cc81Wcxrc1wtYFfG+JDy2tmpBSRMInkYoKCio841SBI+8v5K0FDF2aPf6Reuccy1Qk3nSiCczu9vMRprZyLy8akfCRyWJS47IZ+7q7cxbvT0BETrnXPPUnJLGaiA/4n3v8FhCnDu8F+lpKTwxfVXthZ1zrpVoTkljEvDVsBfVkcB2M/tc1VS8dMhswxlDu/OfWavZs788UbdxzrlmpckkDUmPAu8BAyWVSLpC0lWSrgqLPA8UA0uAe4CrEx3TxUcUsHNPGc/PTVhucs65ZqXJNISb2bhaPjfg240UDgBH9utMYZdMHpu2ii+N6N2Yt3bOuSapyTxpNEWSuOiIfD5YtoXijbuSHY5zziWdJ41aXDCiN6kp4nFvEHfOOU8atemak8FJg7ry1IwS9pdXJDsc55xLKk8aMRg3Kp9Nu/bx6sINyQ7FOeeSypNGDI4ryqN7TgaPT1uZ7FCccy6pPGnEIC01hQtH9mbK4o2s2fZJssNxzrmk8aQRo4tG5lNh8OSMkmSH4pxzSeNJI0b5nTM55qBcHp+2iooKS3Y4zjmXFJ406uDiI/JZve0T3lm6KdmhOOdcUnjSqIPTDu5Gx8w2PDbNx2w451onTxp10DYtlS8d1puX5q9jS+m+ZIfjnHONzpNGHV18RD77y42JM71B3DnX+njSqKOB3bM5rKAjj09bRTCHonPOtR6eNOrhkiPy+XjDLmau3JbsUJxzrlF50qiHM4f1JCs91UeIO+daHU8a9ZDVNo2zDu3JM7PXsnPP/mSH45xzjcaTRj1dfEQ+n+wv55nZvqqfc6718KRRT8PzOzKoe7ZXUTnnWhVPGvUkiYuPyGd2yXYWrNmR7HCcc65ReNJogPMO60V6WgpP+Kp+zrlWwpNGA3TMTGfswd2ZOLOEPfvLkx2Oc84lXK1JQ9JtsRxrrS4dXcCOPWU+ZbpzrlWI5Unj1GqOnR7vQJqr0X07c2h+R+55q5hynzLdOdfC1Zg0JH1L0lxgoKQ5EdsyYE7jhdi0SeJbx/djxebdvDDPu98651q2aE8ajwBnAZPC18rtcDP7SrwDkTRW0keSlki6qZrP+0h6NUxcb0jqHe8Y6uvUId3pl5vF36Ys9fmonHMtWo1Jw8y2m9lyMxtnZiuATwAD2ksqiGcQklKBuwiqvYYA4yQNqVLsd8A/zWwYcAvwm3jG0BCpKWLCcf2Yt3oH7yzZnOxwnHMuYWJpCD9L0sfAMmAKsBx4Ic5xjAKWmFmxme0DHgPOqVJmCPBauP96NZ8n1XkjetE1uy1/m7I02aE451zCxNIQ/kvgSGCxmfUFTgamxjmOXkDkYIeS8Fik2cCXwv3zgGxJXaq7mKQJkqZLmr5x48Y4h1q9tmmpfP2Yvry9ZBNzS7Y3yj2dc66xxZI09pvZZiBFUoqZvQ6MTHBc1fkecLykWcDxwGqg2sERZna3mY00s5F5eXmNFuClowvIbpvmTxvOuRYrLYYy2yS1B94EHpa0ASiNcxyrgfyI973DYweY2RrCJ40wnvPNbFuc42iQnIw2fGVMH/4+ZSnLN5VSmJuV7JCccy6uYnnSOIegEfxG4EVgKUEvqniaBhRJ6ispHbiEoNfWAZJyJVXGezNwf5xjiIuvHV1IWmoKd79VnOxQnHMu7mpNGmZWamblZlZmZg+a2Z1hdVXcmFkZcA0wGVgIPGFm8yXdIunssNgJwEeSFgPdgF/FM4Z46ZqdwfkjevPkjBI27NyT7HCccy6uYuk99SVJH0vaLmmHpJ2S4j6tq5k9b2YDzKy/mf0qPPYTM5sU7j9pZkVhmW+Y2d54xxAvE47rx/7yCh54Z3myQ3HOubiKpXrqt8DZZtbBzHLMLNvMchIdWHPWNzeL04d256GpK3xlP+dcixJL0lhvZgsTHkkLc9Xx/dm5p4xH3vdFmpxzLUeNvackVY6JmC7pceBp4ECVkJlNTGxozduw3h05qn8X7nt7GeOPLqRtWmqyQ3LOuQaL9qRROddUDrAbOC3i2JmJD635u+r4/mzYuZenZ62uvbBzzjUDNT5pmNnXACQdbWbvRH4m6ejqzgknEbwEOBboSdBVdx7wHPCCmVXEKe5m4diiXA7umcPf3yzmwsPzSUlRskNyzrkGiaVN40+xHJP0D4KxE/uA24BxwNXAK8BY4G1Jx9U/1OZHElcd35/ijaW8tGB9ssNxzrkGi9amMQY4CsiT9J2Ij3KA6irof29m86o5Pg+YGA7ai+vsuM3B6UO7U9A5k79OWcoXDu6G5E8bzrnmK9qTRjrQniCxZEdsO4ALqhauLmFI6i/pkPDzfWa2JB5BNydpqSlceVw/Zq/axtTiLckOxznnGiRam8YUYIqkB8L1NOpE0v8CBwEVktqa2WUNiLNZu/Dw3tzxymL+NmUpY/pXOzGvc841C7FMIxJTwpB0XbiYUqVDzezrZvYN4ND6BtgSZLRJZfxRhUxZvJEFa+I+mN455xpNLA3hsdoMvBgxV9RLkl6U9BLBnFKt2mVHFpKVnsqfX/842aE451y9xS1pmNnDBGM4hkmaBMwgmMr8QjP7frzu01x1yGzDFcf24/m565i5cmuyw3HOuXqpMWlIulJSUbgvSf8IJyycI2lEDaf1B54AJgDfBu4A2sU76Obqm8f1I7d9W3793ELMLNnhOOdcnUV70rieYD1wCMZcDAP6At8hSAafIekB4Abgp8B3zOxK4C/APZJ+EreIm7GstmnceGoR01dsZfJ8H7fhnGt+oiWNMjOrnKL1TOCfZrbZzF4BqluS7jAzu9LMvgycCmBms8zsLIL1vR1w8ch8DurantteXMT+8lY1QN451wJESxoVknpIygBOJhjZXam6KqcXJU2W9BrwSOQHZvbfhofaMqSlpnDT2EEs21TKox/4DLjOueYl2hrhPwGmE4z+nmRm8wEkHQ98bi1TM/sfSTlAhZntSkSwLcXJg7syum9n7njlY847rBfZGW2SHZJzzsWkxicNM3sW6AMMDtsnKk0jmJTwMyR9BdhVU8IIR4cf08B4WwRJ/PCLg9lcuo+/TVma7HCccy5m0Z40Ktfu3gpBDyrgJOBSgjaOblWKdwFmSZpB0N12I5BBMCr8eGATcFM8g2/OhvXuyNmH9uTet5bxlSP70KODdzJzzjV9sawRfqSkO4EVwH+BN4FBVcuZ2R3ACOBRII+gHWQEsBq4zMzONzMf2Rbh+18YiBn8bvLiZIfinHMxiTbL7a+BC4GVBIng58B0M3uwpnPMrBx4OdxcLfI7ZzL+6ELueauYK47py5CevvS6c65pi/ak8Q1gPfBX4CEz2wz4iLQ4+/YJB5GT0YbfvODLsDvnmr5oSaMH8EuCqUGWSnoIaCcpajuIq5sOmW249qSDeOvjTUxZvDHZ4TjnXFTRek+Vm9mLZnY5wfQgTwPvAKslPVLTeQ0haaykjyQtkfS5RnNJBZJelzQrnM7kjETE0dguG9OH/M7t+M3zCymv8Ic551zTFbVNw8z+N3x7nJk9BTwVjsU4t5ry36l6LJKZ/SHa5+G06ncRjCYvAaZJmmRmCyKK/Qh4wsz+KmkI8DxQGO26zUHbtFR+8IVBXPvoLJ6aWcJFI/OTHZJzzlUrWvXU2Ij92yp3zGyHmf2zmvLZtWy1GQUsMbNiM9sHPAacU6WMESw3C9ABWBPDdZuFM4f14ND8jvz+pY/4ZF95ssNxzrlqxa19wsx+3sBL9AJWRbwvAUZXKfMzgnU6riWY/+qUBt6zyZDED88YzEV/f4/73i7mmpOKkh2Sc859TrSk0TWsclLE/gFVq5sk/cDMfivpT1TTy8rMrotDvOOAB8zs95LGAA9JGmpmn5n5T9IEgunZKSgoiMNtG8eovp05dUg3/jalmEtGFZDbvm2yQ3LOuc+IVj11D0G1UvuI/WjVTZV9RqcTjAivutVmNRBZmd87PBbpCoL1OjCz9whGnOdWvZCZ3W1mI81sZF5eXgy3bjpuOn0Qn+wv545XfBykc67pqfFJo7K6SVKumW2q7UJm9kz4+mB4XvvwfayTF04DiiT1JUgWlxBMWRJpJcFI8wckDSZIGi2qn2r/vPZcOqqARz5YyaWjCxjcwwf8Oeeajmgr950paSMwR1KJpKNiuaCkoZJmAfOBBZJmSDq4tvPCea6uIVhPfCFBL6n5km6JWHf8u8CVkmYTjFIfby1wCbwbTx1Ax3Zt+J+n5lDma24455oQ1fQ7V9Ic4CIzWyRpNPBbMzu+1gtK7wI/NLPXw/cnAL82s5iSTryNHDnSpk+fnoxbN8gzs9dw7aOz+N8zBjHhuP7JDsc514pImmFmI6v7rLaV+xYBmNn7xNZtFiCrMmGE575B9Sv9uSjOHNaDU4d04/cvLWb5ptJkh+Occ0D0pNFV0ncqt2re16RY0o8lFYbbj6hm0SYXnSR+ee5Q0tNS+J+n5lDhI8Wdc01ALL2nKreq72vydYKp0SeGW154zNVRt5wMfnjGYN5ftoVHp/nSsM655Ku191RdmdlW4DoASZ2AbS2xsbqxXHxEPpNmr+E3zy/ipEFdfbEm51xS1boIU6wk/UTSoHC/raTXgCXAekktZuR2Y5PErV8aRllFBT/8zzw8/zrnkiluSQO4GPgo3L88vHZXgqVefx3H+7Q6BV0y+d5pA3lt0QYmzW4x020555qheCaNfRHVUF8AHg2nV19IHOe4aq2+dnRfhud35GeT5rN5195kh+Oca6ViWSP8O9VsV0gaXqXo3nBgXx5wIvBSxGeZcYy5VUpNEb+9YBi79pbx82cW1H6Cc84lQCxPGiOBqwhmoe0FfJNg2vR7JP0gotz1wJPAIuCPZrYMIFwoaVY8g26tBnTL5poTi5g0ew2vLFif7HCcc61QjSPCDxSQ3gTOqJxDKpxT6jmCxDHDzIYkPMoGaK4jwmuyr6yCs//8Nlt37+Pl7xxPTkabZIfknGth6jsivFJXILISfT/Qzcw+qXLcNYL0tBRuO38YG3fu5TfPL0p2OM65ViaWpPEw8L6kn0r6KcE64Y9IygK8cj0JDs3vyJXH9uPRD1by7tJaJyB2zrm4qTVpmNkvCNoxtoXbVWZ2i5mVmtmXExueq8kNpwygsEsmNz0115eHdc41mli73M4E/g38B9ggqcbl8CRlhnNP3RO+L5J0ZsNDdZHapady6/nDWLllN7e96NVUzrnGEUuX22uB9cDLwLMEjeDPRjnlHwRtHWPC96uBXzYsTFedI/t14WtHF/LAu8t5Ye7aZIfjnGsFYhl0dz0w0Mw2x3jN/mZ2saRxAGa2W5LqHaGL6ubTBzNr5TZ+8OQcBvfIoTDXZ6F3ziVOLNVTq4DtdbjmPkntAAOQ1B/vZZUw6Wkp3PXlEaSmim89PJM9+719wzmXOLEkjWLgDUk3x7iexk+BF4F8SQ8DrwI/iFLeNVCvju3448XDWbh2Bz+bND/Z4TjnWrBYqqdWhlt6uEVlZi9LmgkcCQi43sy8X2iCnTiwK98+sT93vb6UIwo7c/7hvZMdknOuBao1adRzXY0MYGt4/SGSMLM363EdVwc3njKAGSu28sOn5zK0VwcGdo91hV7nnItNjUlD0u1mdoOkZwjbJyKZ2dk1nHcbwTTp84GKyuKAJ40ES0tN4c5xh/HFO9/mWw/PYNI1x9C+rU8w7JyLn2i/UR4KX39Xx2ueS9Dbyhu/k6BrdgZ3XnIYX753KjdPnMudlwzHO6855+Il2nKvM8LXKXW8ZjHQBu8xlTRj+nfhu6cN5P8mf8Sowk5cNqYw2SE551qIWusuwtHcvwD6hOUFmJnlVCn3J4JqqN3Ah5JeJSJxmNl1cYzb1eJbx/dn+vIt/OLZhRya35FhvTsmOyTnXAsQS5fb2wmWb+1iZjlmll01YYSmAzOASQRJ5t3w/Yzws6gkjZX0kaQlkm6q5vM/Svow3BZL2hZD7K1WSor4w0XDyctuy9UPz2T77v3JDsk51wLEOrhvntWy8IaZPWhmDwIdK/cjjnWKdq6kVOAu4HRgCDBO0mfW6TCzG81suJkNB/4ETIwh9latU1Y6f770MNbv2MN3//0hFRXR105xzrnaxJI0fgA8X4fBfZdXc2x8LfcYBSwxs2Iz2wc8BpwTpfw44NFarumAwwo68b9nDOaVhRu4+63iZIfjnGvmYumP+StgF8HYixoH94VzTV0K9JU0KeKjbGBLLffoRfBEU6kEGF3DffoAfYHXosQyAZgAUFBQ44S8rcb4owqZvnwr/zf5IwZ2y+bEQV2THZJzrpmKJWn0NLOhMZR7F1gL5AK/jzi+E5hTj9hqcgnwpJnVOMmSmd0N3A3Bcq9xvHezJInbLhjGii2lXP3wTB7/5pHeMO6cq5dYqqeel3RabYXMbIWZvWFmY8xsSsQ208zKajl9NZAf8b53eKw6l+BVU3XWvm0a948/gi7t0/n6A9NYuXl3skNyzjVDsSSNbwEvStojaWe47YhzHNOAIkl9JaUTJIZJVQtJGkTQqP5enO/fKnTNzuDBr4+irMK4/B8fsHmXD6VxztVNLMu9ZptZipllhPs1dbmtt/BJ5BpgMrAQeMLM5ku6RVLkdCWXAI/V1pPL1ax/Xnvuu3wka7Z9whUPTvelYp1zdaJYfv+Gv7iPC9++YWbRVu5rUkaOHGnTp9c6TKTVmTx/Hd/61wxOGtSVv33lcNJSY1351znX0kmaYWYjq/ssluVebyVYvW9BuF0v6TdRyh8t6eVwAF6xpGWSvK9nE/OFg7vz87MP5pWFG/jxf+fjD2/OuVjE0nvqDGC4mVUASHoQmAXcXEP5+4AbCUaCe91HE3bZmELWbt/DX95YSs8OGVx7clGyQ3LONXGxzpvdkU/HWnSopex2M3uh3hG5RvX9Lwxk3fY9/P7lxXTvkMGFI/NrP8k512rFkjR+A8yS9DrBZIXHAZ+bGyrC65L+j2Caj8gJC2c2JFCXGJK49fxhbNi5l5smziUvuy0nDPTBf8656sXaEN4DOCJ8+4GZrYtS9vVqDpuZnVS/EBvGG8Jjs3PPfi7++1SWby7l8QljOKR3bQ+UzrmWKlpDeK1JQ9J5wGtmtj183xE4wcyejnOcCeFJI3YbduzhvL+8y96yCiZ+6ygKumQmOyTnXBI0qPcU8NPKhAFgZtuAn1Zzk6+Er9+pbqtn7K4Rdc2pHPxXwbh7pvqocefc58SSNKorU11bSFb4ml3D5pqBg7q2519XjKZ0XxkX/f09ijfuSnZIzrkmJJbqqfuBbQTrXQB8G+hsZuMTGlmcePVU/Sxcu4Ov3Ps+KSnikW+Mpqib533nWouGVk9dC+wDHg+3vQSJw7Vgg3vk8NiEIwG45O6pLFwb7+nGnHPNUSxzT5Wa2U1mNjLcbjaz0sYIziVXUbdsnvjmGNLTUhh3z1Tmlmyv/STnXIsWyzQiAyTdLeklSa9Vbo0RnEu+vrlZPPHNMWSlp3HpvVOZuXJrskNyziVRLNVT/yaYNuRHwPcjtmpJul5SjgL3SZoZy3ocrunK75zJE1eNoXNWOpfd+z4fLKttIUbnXEsVS9IoM7O/mtkHZjajcotS/utmtgM4jWDti8uAW+MRrEueXh3b8cQ3x9C9QwaX3/8B7y7ZlOyQnHNJEEvSeEbS1ZJ6SOpcuUUpr/D1DOAhM5sfccw1Y91yMnhswhgKOmfytQemMWXxxmSH5JxrZLEkjcsJqqPeJZi5dgYQrQ/rDEkvESSNyZKygYqGBuqahrzstjw64Uj657Xnygen88qC9ckOyTnXiGLpPdW3mq1flFOuIJjQ8Agz2w2kA1+LU7yuCeiclc6jVx7J4B7ZXPWvGTz2wcpkh+ScayQ1znIr6UvRTjSziTUcr5C0HhgiKdap110z0yGzDf/6xmi+/cgsbpo4l6Ubd3HT6YNJTfGaSOdasmi/1M+K8pkRTH3+OZJuAy4mWOWvPKL8m/UJ0DVd2RltuP/ykfzi2QXc89YyijeWcse4w2jf1v9WcK6limlq9DpdUPoIGGZme2st3Ah8GpHG8c/3lvPzZxZQ1LU9940/gl4d2yU7JOdcPTV0GpG6KgbaJOC6rgn76phC7h9/BKu3fsI5f36HWT4I0LkWKRFJYzfwoaS/S7qzckvAfVwTc/yAPCZefRTt0lO45O6pPDN7TbJDcs7FWY1JQ9KF4WvfOl5zEvALPttFN9pgQNeCFHXL5umrj2ZY7w5c++gs7njlY+JdBeqcS54a2zQkzTSzEZWvjRxX3HibRnLsLSvn5olzmThzNWcf2pPfXjCMjDapyQ7LOReDaG0a0bq5bA4H6fWVNKnqh2Z2dpWbPGFmF0maS9Bbqmr5YbUEORa4A0gF7jWzz009Iuki4Gfh9Web2aXRrumSp21aKr+/8FAO6tqe3774Eau27ubvlx1O1+yMZIfmnGuAaE8a6cAI4CHgG1U/N7MpVcr3MLO1kvpUdz0zW1FjEFIqsBg4FSgBpgHjzGxBRJki4AngJDPbKqmrmW2o5fv5k0YT8OK8tdzw+IdkZ7Th9ouHc/RBuckOyTkXRb2eNMxsHzBV0lFmtlFS+/B4tet/mtna8LXG5BDFKGCJmRWHAT8GnEMw1qPSlcBdZrY1vE+tCcM1DWOH9uA/XbK45pGZfOW+97n6hP7ceMoA0lIT0Q/DOZdIsfxf203SLGA+sEDSDElDayos6UhJ0yTtkrRPUrmk2pZ96wWsinhfEh6LNAAYIOkdSVPD6qyaYpggabqk6Rs3+qR6TcHgHjk8c+0xXHh4b+56fSkX3z2Vkq27kx2Wc66OYkkadwPfMbM+ZlYAfDc8VpM/A+OAj4F2BFVbd0UpH6s0oAg4Ibz+PZI6VlfQzO6uXGkwLy8vDrd28ZCZnsZvLziUOy4ZzkfrdnLGHW/x4ry1yQ7LOVcHsSSNLDN7vfKNmb0BZEU7wcyWAKlmVm5m/wBqfCoIrQbyI973Do9FKgEmmdl+M1tG0AZSFEP8rok5Z3gvnrvuGApzs7jqXzP58dPz2LO/vPYTnXNJF0vSKJb0Y0mF4fYjglHfNdkdNqJ/KOm3km6M4T7TgCJJfcNzLyEY7xHpaYKnDCTlElRXRYvDNWF9umTx5FVHceWxfXlo6grOvesdlmzYmeywnHO1iCVpfB3II5ig8CkgNzxWk8sIus1eA5QSPEGcH+0GZlYWlp8MLASeMLP5km6RVNm1dzJBN+AFwOvA981scwzxuyYqPS2FH35xCP8YfwQbdu7lrD+9wxPTV/lgQOeasLhPWNjUeJfb5mH9jj3c8NiHvFe8mbMP7ckt5xxMx8z0ZIflXKvUqBMWSjpT0ixJWyTtkLQzht5TrpXrlpPBv74xmu+eOoDn5q7llD9M4bk5a/2pw7kmJhEd5W8nWCK2i5nlmFm2meUk4D6uhUlNEdeeXMSka46me4cMvv3ITCY8NIN12/ckOzTnXKjWpCHp6FiORVgFzDP/E9HV08E9O/D01Udz8+mDeHPxRk79wxQefn8FFRX+T8q5ZKu1TaO6CQujTWIo6QiCWW6nAAcWYjKzPzQ83LrzNo3mbfmmUm6eOJf3ijczqm9nbv3SIfTLa5/ssJxr0eo1jYikMcBRQJ6k70R8lEPQO6omvwJ2ARmAt2S6BinMzeKRK0fzxPRV/PK5hYy94y1uOKWIK4/tRxufhsS5Rhdtltt0oH1YJjvi+A7ggijn9TSzGqcZca6uJHHxEQWcOLArP/nvfH774kc8O3stt50/jEN6d0h2eM61KrFUT/WpyySEkn4LvGJmLzU0uHjw6qmW58V5a/nxf+ezpXQfXx3ThxtOHkCHTF9h2Ll4iVY9FUvSGAB8Dygk4snEzE6qofxOgmlG9oWbguLJ6UHlSaNl2v7Jfm59YRGPTVtJh3ZtuPGUAVw6usCrrJyLg4YmjdnA3wiWbD0wQZCZNYslXD1ptGwL1uzgl88t4N2lm+mfl8WPzhzCiQO7Jjss55q1hiaNGWZ2eB1uJuDLQF8z+4WkfKCHmX1Ql6DjxZNGy2dmvLJwA79+fiHLNpVy3IA8fvTFwQzoll37yc65z2noiPBnJF0tqYekzpVblPJ/AcYAlUux7iI+U6M7Vy1JnDqkG5NvOI4fnzmED1du5fQ73uJHT89l8669tV/AORezWJ40llVz2MysXw3lZ5rZCEmzzOyw8NhsMzu04eHWnT9ptD5bS/dx+yuL+df7K8lMT+Xakw7i8qMKaZsWrae4c65Sg540zKxvNVu1CSO0P1zz28Kb5wEV9YrcuXrolJXOz88ZyuQbjmVkn078+vlFnPqHN3lqRgll5f5P0bmGiGUakUxJP5J0d/i+SNKZUU65E/gP0FXSr4C3gV/HJVrn6uCgrtn842ujePDro8jOSOO7/57NqX98k//MKqHcpyRxrl5iqZ56nKDn1FfNbKikTOBdMxteTdkU4EhgC3AyQXfbV81sYbwDj5VXTzkIGstfWrCe21/5mIVrd9AvN4vrTi7irEN7kpqiZIfnXJPS0N5T081sZKxtFJHlmgJPGi5SRUVl8ljMonU76ZeXxfUnF3HmME8ezlVqaO+pfZLa8WkbRX8iJiKsxquSzg+73jrXpKSkiLFDu/P8dcfy1y+PoE1KCtc/9iFfuP1NJs1e49VWztUilieNU4EfAUOAl4CjgfFm9kYN5StHhJcBe/AR4a4Jq6gwXpi3jjteXczi9bso6tqeq0/sz5nDevroctdqNah6KrxAF4K2CgFTzWxTfENMHE8aLhYVFcbz89Zy56sfs3j9LrrnZDD+6ELGjSqgQzuf18q1Lg1t0zgPeM3MtofvOwInmNnTcY4zITxpuLqoqDCmLN7IPW8V8+7SzWSmp3LRyHyuOKYv+Z0zkx2ec42ioUnjw6o9pZpaY3c0njRcfc1fs5373lrGpNlrqDBj7NDufOPYfowo6JTs0JxLqIYmjTlmNqzKsblmdkgcY0wYTxquodZt38OD7y3n4akr2LGnjBEFHbny2H6cdnB373HlWqSGJo37gW18On/Ut4HOZjY+jjEmjCcNFy+le8v49/RV3P/OclZu2U1+53aMG1XAhYfnk5fdNtnhORc3DU0aWcCPgVMIut2+DPzKzErjHORY4A6CpWTvNbNbq3w+Hvg/YHV46M9mdm9t1/Wk4eKtvMJ4ecE6Hnh3OVOLt5CWIk47uBuXjurDUf27kOJPH66Zq9ca4eGJqcCzZnZiQiL77H3uAk4FSoBpkiaZ2YIqRR83s2sSGYtztUlNEWOH9mDs0B4s3biLR99fyZMzS3h+7jr6dMlk3KgCLji8N7nt/enDtTxRO6KbWTlQISnRCzGPApaYWbGZ7QMeA85J8D2da7D+ee350ZlDmHrzydx+8XC65WRw6wuLGPObV/n2IzN5d8kmKnzAoGtBoj5phHYBcyW9DByokjKz6+IYRy9gVcT7EmB0NeXOl3QcsBi40cxWVVMGSROACQAFBQVxDNO56mW0SeXcw3px7mG9WLJhJ4+8v4qnZpbw3Jy1FHbJ5MKR+ZwzvCe9O3m3Xde8xdKmcXl1x83swbgFIV0AjDWzb4TvLwNGR1ZFhQMMd5nZXknfBC6uaZ3ySN6m4ZJlz/5yXpi3lkffX8UHy7cAMLpvZ740ohdjh/bwQYOuyYrHiPB2QIGZfRTv4MLrjwF+ZmZfCN/fDGBmv6mhfCqwxcxqrTbzpOGagpWbd/P0h6t5etZqijeVkp6WwqmDu3HuYb04fkAe6Wk+ZYlrOhrae+os4HdAupn1lTQcuMXMzo5jgGkEVU4nE/SOmgZcambzI8r0MLO14f55wP+Y2ZG1XduThmtKzIzZJdt5etZqnpm9hs2l++iU2YYzh/XkvBG9OCy/Iz7Xp0u2hiaNGcBJwBsRU6PPM7OhcQ7yDOB2gi6395vZryTdAkw3s0mSfgOcTTAR4hbgW2a2qLbretJwTdX+8gre+ngjE2eu5uUF69lbVkGfLpmccUgPTh/anUN6dfAE4pKioUljqpkdWWU9jc+NEm+qPGm45mDnnv28OG8dk2av4d2lmymvMHp1bMcZh3Rn7NAeHJbf0cd/uEZT73EaofmSLgVSJRUB1wHvxjNA51q77Iw2XDgynwtH5rNt9z5eXrCeF+YFAwjveWsZ3XMyGDu0O2cc0oPD+3Ty6Utc0sTypJEJ/BA4LTw0Gfilme1JcGxx4U8arjnbsWc/ry5czwtz1/HG4o3sK6sgt31bxg7txmlDujO6X2fapqUmO0zXwtSrekpSBnAVcBAwF7jPzMoSFmWCeNJwLcWuvWW8vmgDL85bx2uLNvDJ/nIy01M5tiiXkwd144RBeXTNzkh2mK4FqG/11IPAfuAt4HRgMHBD3KNzzsWkfds0zjq0J2cd2pM9+8t5d+kmXl24gdcWbWDy/PUAHNq7AycN6sbJg7tycM8cb0h3cRftSePA9Odhl9gPzGxEYwYXD/6k4Vo6M2Ph2p28tmg9ry7awIertmEG3XLactKgrpw4sCtj+nchO8MHE7rY1PdJY3/ljpmV+V8szjVNkhjSM4chPXO45qQiNu3ayxsfbeS1Ret5ZvZaHv1gFWkp4rCCjhxblMcxRbkM69WBNF8D3dVDtCeNcj6da0pAO2B3uG9mltMoETaQP2m41mxfWQXTV2zh7Y838dbHm5i3ZjtmkJORxtEH5XJMUS7HFeX5UrbuMxo8jUhz5knDuU9tKd3HO0s2hUlkI2u2B50g+3TJ5NiiXI7un8uovp3p4tO6t2qeNDxpOPc5ZkbxplLeWryRt5ds4r2lmyndVw7AgG7tObJfF47s14VRfTv72iCtjCcNTxrO1Wp/eQVzV2/n/eItTC3ezLTlW9gdJpGirkESGd2vM6P7dvHlbVs4TxqeNJyrs/3lFcxbvZ2pxVt4f9lmpi3bcuBJpH9eFiP7dObwwk6M7NOJvrlZ3r23BfGk4UnDuQYrK69g3podTC3ezAfLtjBjxVa2fxJ0suyclc6Igk6MLOzE4X06cUivDmS08ZHqzZUnDU8azsVdRYVRvGkX05dvZfqKrcxcsZXiTUGHy/TUFIb2yuHwPp0YUdCJQ/M70qNDhj+NNBOeNDxpONcoNu/ay4wVW5mxciszlm9lTsl29pVXAJCX3Zbh+R0PbIf07kCODzhskho6y61zzsWkS/u2nHZwd047uDsAe8vKWbh2J7NXbWP2qm18uGobLy9Yf6B8/7wshud3Ynh+Bw7N78jA7tk+AWMT50nDOZcwbdNSDzxZVNq+ez9zVm/jw5XbmF2yjSmLN/DUzBIA0lLEgG7ZHNKrA0N75XBwrw4M6ZHj7SNNiFdPOeeSysxYve0T5pRsZ97q7cxdHbxu3R00sqemiIPy2jM0TCSH9OrAoB45tG/rf/MmildPOeeaLEn07pRJ707BUrcQJJI12/cwL0wg81ZvZ8rijQeeSCAYxT64ew6De+QwuEc2Q3rm0KtjO29sTzBPGs65JkcSvTq2o1fHdnwhbB8BWL9jD3NLtrNw7Q4WrtvBwrU7mbxgHZUVJjkZaQzqkcOQMJEM7pHDQV3bk5nuv+rixX+Szrlmo1tOBt2GZHDKkG4HjpXuLWPRup1BIgm3J6avOjCaXYKCzpkM6JbNwG7ZDOgevPbNzSI9zWf6rStPGs65Zi2rbRqH9wkGFVaqqDBWbtnNwrU7WLx+F4vX7+Sj9Tt5bdEGyiuCx5K0FNEvL+tAMinq1p6DuranT5cs2vi08TXypOGca3FSUkRhbhaFuVmcfsinx/eWlVO8sTRIIut2snj9TmaXbOPZOWsPlEkLzy3qGiSRyq1/XnvvxYUnDedcK9I2LTVsOP/sckCle8tYunEXSzbs4uMNwetH63Yyef46wgcTJOjdqR0H5bWnX157+uVl0S+3Pf3zssjLbttqGuA9aTjnWr2stmkM692RYb07fub43rJylm/azccbdrIkTCZLNuzi3aWb2VtWcaBcdts0+uZl0S836zMJpW9uFu3SW9bTSZNKGpLGAncAqcC9ZnZrDeXOB54EjjAzH4ThnEuItmmpDOyezcDu2Z85XlFhrNn+CcUbSyneuIviTaUUbyzlg2VbePrDNZ8p2z0ng8LcTAq7BNVlhV2y6JubRZ8umc2yuqvJJA1JqcBdwKlACTBN0iQzW1ClXDZwPfB+40fpnHNBm0nl2JLjBuR95rPd+8pYFiaR5ZtKWba5lBWbd/PygvVsLt33mbI9OmQcSCYFnTPp0yXzwGt2E52Xq8kkDWAUsMTMigEkPQacAyyoUu4XwG3A9xs3POecq11mehoH9+zAwT07fO6zHXv2s3xTKcs37w5eN5WyfHMpk+evY0uVhNI5K/1AAunTOZOCLsHTSX6nTLpmtyUlJTltKE0pafQCVkW8LwFGRxaQNALIN7PnJNWYNCRNACYAFBQUJCBU55yru5yMNtW2nQDs3LOflVt2s3LzblZs2c2KzbtZuaWUGSu28szsNQca5AHS01Lo3akd+Z0yye9c+Zp54H2Hdm0S1jDflJJGVJJSgD8A42sra2Z3A3dDMPdUYiNzzrmGy85oU+MTyr6yCtZs+4QVW3azcstuSrbsZtXW3aza8gmzS7axLZyn68C12qbRu3MmE791VNwb4ptS0lgN5Ee87x0eq5QNDAXeCDNod2CSpLO9Mdw515Klp6UcGHdSnR179rNqS5BESrbuZtWW3WzctTchPbeaUtKYBhRJ6kuQLC4BLq380My2A7mV7yW9AXzPE4ZzrrXLifKUEm9NZqy8mZUB1wCTgYXAE2Y2X9Itks5ObnTOOeegaT1pYGbPA89XOfaTGsqe0BgxOeec+1STedJwzjnX9HnScM45FzNPGs4552LmScM551zMPGk455yLmScN55xzMZNZy55lQ9JGYEU9T88FNsUxnObAv3Pr0Nq+c2v7vtCw79zHzPKq+6DFJ42GkDTdzEYmO47G5N+5dWht37m1fV9I3Hf26innnHMx86ThnHMuZp40ors72QEkgX/n1qG1fefW9n0hQd/Z2zScc87FzJ80nHPOxcyThnPOuZh50qiGpLGSPpK0RNJNyY4n0STlS3pd0gJJ8yVdn+yYGoukVEmzJD2b7Fgag6SOkp6UtEjSQkljkh1Tokm6Mfx3PU/So5Iykh1TvEm6X9IGSfMijnWW9LKkj8PXTvG4lyeNKiSlAncBpwNDgHGShiQ3qoQrA75rZkOAI4Fvt4LvXOl6gkW/Wos7gBfNbBBwKC38u0vqBVwHjDSzoUAqwaqgLc0DwNgqx24CXjWzIuDV8H2DedL4vFHAEjMrNrN9wGPAOUmOKaHMbK2ZzQz3dxL8IumV3KgST1Jv4IvAvcmOpTFI6gAcB9wHYGb7zGxbUoNqHGlAO0lpQCawJsnxxJ2ZvQlsqXL4HODBcP9B4Nx43MuTxuf1AlZFvC+hFfwCrSSpEDgMeD/JoTSG24EfABVJjqOx9AU2Av8Iq+TulZSV7KASycxWA78DVgJrge1m9lJyo2o03cxsbbi/DugWj4t60nAHSGoPPAXcYGY7kh1PIkk6E9hgZjOSHUsjSgNGAH81s8OAUuJUZdFUhfX45xAkzJ5AlqSvJDeqxmfB2Iq4jK/wpPF5q4H8iPe9w2MtmqQ2BAnjYTObmOx4GsHRwNmSlhNUQZ4k6V/JDSnhSoASM6t8inySIIm0ZKcAy8xso5ntByYCRyU5psayXlIPgPB1Qzwu6knj86YBRZL6SkonaDSblOSYEkqSCOq5F5rZH5IdT2Mws5vNrLeZFRL8N37NzFr0X6Bmtg5YJWlgeOhkYEESQ2oMK4EjJWWG/85PpoU3/keYBFwe7l8O/DceF02Lx0VaEjMrk3QNMJmgp8X9ZjY/yWEl2tHAZcBcSR+Gx/7XzJ5PXkguQa4FHg7/ICoGvpbkeBLKzN6X9CQwk6CX4Cxa4JQikh4FTgByJZUAPwVuBZ6QdAXB8hAXxeVePo2Ic865WHn1lHPOuZh50nDOORczTxrOOedi5knDOedczDxpOOeci5knDZdwks6VZJIGJeHeyyXlNtK9eobdO2sr9274Wijp0hjKf6acpJGS7mxYtHUjabyknvU47wZJX61D+XRJb4bzRLkmyJOGawzjgLfD12ajrr+4zGyNmV0QQ7nKEcmFQK1Jo2o5M5tuZtfVJbZYhDM812Q8wTQcdbleGvB14JFYzwknCX0VuLgu93KNx5OGS6hwPqtjgCuImJJa0gmS3ohY2+HhcMRu5dPBzyXNlDS38glF0s8kfS/iGvPCCRaR9LSkGeG6CRNiiGuXpD+G5V+VlBcef0PS7ZKmA9dLOlzSlPDakyOmZThI0iuSZodx9g+fCOaFn4+X9N/weh9L+mnkvcPdW4FjJX0YrvlQKOmt8HozJR1VQ7kTFK7/oWDNhKclzZE0VdKwiJ/V/eH9iyVVm2TCn8PvJc0Gxkj6iaRp4c/2bgUuAEYSDAr8UFK7mn4uVZwEzDSzsoif7R8lTVewlscRkiaGP59fRpz3NPDl2v4buiQxM998S9hG8D//feH+u8Dh4f4JwHaCub1SgPeAY8LPlgPXhvtXA/eG+z8Dvhdx7XlAYbjfOXxtFx7vEnGt3GriMuDL4f5PgD+H+28Afwn324Qx54XvLyaYIQCCWYDPC/czCKbcLgTmhcfGE8yq2iUippHhZ7sifgbPRsSUCWSE+0XA9BrKHXgP/An4abh/EvBhxM/qXaAtkAtsBtrU8HO4KOJ954j9h4CzIn4uI2v7uVS59s8r/ztGXOO2cP96ginKe4QxlkT8N0sFNib7365v1W9eb+gSbRzBwj8QTAw4DqicWfYDMysBUDB9SSFBNRYEE8sRlv1SDPe5TtJ54X4+wS/dzVHKVwCPh/v/irgfEccHAkOBl8OHoFRgraRsoJeZ/QfAzPaE36HqPV42s83hZxMJnrimR4mpDfBnScOBcmBAlLKVjgHOD+N4TVIXSTnhZ8+Z2V5gr6QNBFNjl1Q5v5xgospKJ0r6AUEC6wzMB56pck61P5dqYuvB5+d5qpzHbS4w38KpuyUVE/x322xm5ZL2Scq2YH0X14R40nAJI6kzwV+/h0gygl8uJun7YZG9EcXL+ey/x73VHC/js1WqGeF9TiCYzXSMme2W9EblZ3UQOZ9OaeVXIPjF9pklUcOkUddrVve+qhuB9QQr6qUAe2K8T02i/Xwr7TGzcgAFy6D+heCJYpWkn1H9z7Han0s1Pqnm/MqYKqrEV1ElvrY0/Pu7BPA2DZdIFwAPmVkfMys0s3xgGXBsPa+3nHAqb0kjCNZIAOgAbA0TxiCCJWtrkxLGB0Ej89vVlPkIyFO4jrakNpIODv/6LZF0bni8raTMas4/NWxzaEewato7VT7fCUQmoA7AWjOrIJhAMrWGcpHeIqz/D5PnJqv/WiiVv+A3KWiLimzUj4yh2p9LNddbCBxU1yAkdSH4Hvvreq5LPE8aLpHGAf+pcuwp6t+L6imgs6T5wDXA4vD4i0CapIUEjcZTY7hWKTAqbLg+CbilagELevJcANwWNhR/yKdrMVxGUCU2h6B+v3s19/ggjHkO8JSZVa2amgOUh43pNxL8lX95eK9BfPrEU7VcpJ8Bh4dx3MqnU2HXmQVLv95D0P4ymWCZgEoPAH8LqxFTqfnnEukFguVl6+pE4Ll6nOcagc9y61olSbvMrH0Crz+eoJrnmkTdozmQ9B/gB2b2cR3OmQjcZGaLay3sGp0/aTjnEukmggbxmChY5+NpTxhNlz9pOOeci5k/aTjnnIuZJw3nnHMx86ThnHMuZp40nHPOxcyThnPOuZj9P03Fy6saF2alAAAAAElFTkSuQmCC\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 = pspd.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": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAi30lEQVR4nO3deVxU5eIG8OcwgAiDqGkuYKKSGMMyCCjUVcFAURa1XCBviuSt1KybpXaz0tLSrlZel1y6lNbPwC33XJJAKTUFRc0p3BXITHJjX2be3x/mXBdUFGbeAZ7v58Pn4xzOeeeZER6P75xFEUKAiIjMz0p2ACKi+ooFTEQkCQuYiEgSFjARkSQsYCIiSVjARESSWN/Pys2aNROurq4mikJEVDdlZGTkCSGa37r8vgrY1dUV6enpNZeKiKgeUBTlTGXLOQVBRCQJC5iISBIWMBGRJCxgIiJJWMBERJKwgImIJGEBExFJwgImIpKEBUxEJAkLmIhIEhYwEZEkLGAiIklYwEREkrCAiYgkYQETEUnCAqZKqdVqixk7NTUVu3btMlEaInlYwGTxWMBUV7GAqcpOnDiB8PBw+Pn5oVu3bvj1118BABs2bEDXrl3h6+uL0NBQnD9/HgBQUFCAESNGwMvLC97e3li9erVxrEmTJsHHxweBgYHG9Ssb5/Tp01i4cCE++eQTaLVapKWl4cKFC3j66acREBCAgIAA/Pjjj+Z/M4hqghCiyl9+fn6C6gcHB4fblvXs2VMcPXpUCCHEnj17REhIiBBCiIsXLwqDwSCEEOKzzz4T48aNE0IIMWHCBPHKK68Yt7948aIQQggAYv369UIIIcaPHy+mTp1613EmT54sZs6caRwnNjZWpKWlCSGEOHPmjOjUqVPNvGgiEwGQLirp1Pu6JxzVXwUFBdi1axcGDRpkXFZaWgoAyMnJwZAhQ3Du3DmUlZWhXbt2AIDt27cjKSnJuH6TJk0AALa2toiMjAQA+Pn54bvvvrvrOLfavn07dDqd8fHVq1dRUFBg0nlrIlNgAVOVGAwGNG7cGJmZmbd9b+zYsRg3bhyio6ORmpqKKVOm3HUsGxsbKIoCAFCpVKioqLivcQwGA/bs2QM7O7vqvCQi6TgHTFXSqFEjtGvXDitXrgRwberq4MGDAIArV67A2dkZALB06VLjNmFhYZg/f77x8aVLl+76HHcax9HREfn5+cbHvXr1wty5c42PK/tHgag2YAFTpYqKiuDi4mL8+vjjj7Fs2TIkJCTAx8cHGo0G69atAwBMmTIFgwYNgp+fH5o1a2Yc46233sKlS5fg6ekJHx8fpKSk3PU57zROVFQU1qxZY/wQbs6cOUhPT4e3tzc8PDywcOFC07wJRCamXJsfrhp/f3+Rnp5uwjhERHWPoigZQgj/W5dzD5iISBIWMBGRJCxgMostW7bA3d0dbm5umDFjhuw4RBaBBUwmp9frMWbMGGzevBk6nQ6JiYk3HcdLVF+xgMnk9u7dCzc3N7Rv3x62traIiYkxHkFBVJ+xgMnkcnNz0aZNG+NjFxcX5ObmSkxEZBlYwEREkrCAyeScnZ2RnZ1tfJyTk2M8442oPmMBk8kFBATg2LFjOHXqFMrKypCUlITo6GjZsYik48V4yOSsra0xb9489O7dG3q9HvHx8dBoNLJjEUnHU5GJiEyMpyITEVkYFjARkSQsYCIiSVjARESSsIDJbDIyMlBWViY7BpHFYAGTyVVUVOD5UaMR9EQ3dAsJxcWLF2VHIrIILGAyuSVLl+KzhQtQXlqMvbvS8Mq412VHIrIILGAyuWdiY/Hzzz8DABISEvDJrH9LTkRkGVjAZHL29vbGM986dux40w03ieozFjARkSQsYCIiSVjARESSsICJiCRhARMRScICJiKShAVMRCQJC5iISBIWMBGRJCxgIiJJWMBERJKwgImIJGEBk1kpiiI7ApHFYAGTWQghAACtWrWSnITIcrCAySxmzpwJAHB1dZUbhMiCsIDJLCZOnIi2bdvCyoo/ckTX8beBTK5nz54AgNTUVLlBJCooKEBWVpbsGGRhWMBkMsnJyVAUBSkpKVi1alW9nn6YPXs2fP0D8Ouvv8qOQhZEuf7hSFX4+/uL9PR0E8apP/R6PQbGDMXlK1eqPVbqd1tqINHdBYeFV3ndW/M0b9ESGm9tDSeqXU6fPIHf8ivQzE7g0P50PPTQQ7IjkRkpipIhhPC/dbm1jDAElJeXIysrC2cKreCgjajWWA07/InSc0drKNnt7N264ohjQJXXb9TFHiXZP8PRty+s7BwBAEdMFa628AlAi0e8UPTTCoRH9sOPO76Hra2t7FQkGfeAJcrNzYWPXwCsgkbA3v1x2XHIDIRBj/yNH6KP/6P4v6Vf8LjoeuJOe8CcA5bI2dkZWzdtQFHKApT+flx2HDIDxUoFdfir2JTyIz6cOUt2HJKMBSyZn58fPl+8CAUbp6Mi/0/ZccgMrGwbQh31JqbNmIn169fLjkMSsYAtwODBg/Day2NQsHE6DOUlsuOQGVg3ehiOkW9g6LA4HDp0SHYckoQFbCEmv/0WwoJ8UbhtDoQwyI5DZtCgtTvsuo9EWJ8InD9/XnYckoAFbCEURcFXSz5HmwYlKNyVKDsOmYmDRw/oO3RHr75RKCnh/37qGxawBbGzs8PWTeuhOpmGQl2q7DhkJg5BMcgpt8czw+JwP0clUe3HArYwLVq0wPYt36Ik7XOU5vKsqfpAUazgEPYyvv/pEKa8N/WO66lUKmi1Wnh6eiIqKgqXL182aa6RI0dCp9PddZ2FCxfiyy+/NGmOuozHAVuoDRs2IHZ4PBrH/BvWjR6WHYfMoKLgIq4kjcd/58/BkCGDb/u+Wq1GQUEBAGD48OHo2LEjJk2aZO6Y9AB4HHAtExUVhXf+NREF6z+AobRIdhwyA2t1UzhGvYnnXngR99rRCQoKQm5uLgDgxIkTCA8Ph5+fH7p162a83kRcXBxGjRqFwMBAtG/fHqmpqYiPj8djjz2GuLg441ijRo2Cv78/NBoNJk+ebFweHBxszKFWqzFp0iT4+PggMDDQ+KHhlClTMGvWLOP6EydORJcuXdCxY0ekpaUBAIqKijB48GB4eHhgwIAB6Nq16z1fX33BArZg419/DVGh3VCw9RMIg152HDID2xYdYN9zNMIjopCTk1PpOnq9HsnJyYiOjgYAPP/885g7dy4yMjIwa9YsjB492rjupUuXsHv3bnzyySeIjo7Gq6++iiNHjuDw4cPIzMwEALz//vtIT0/HoUOHsGPHjkoPiyssLERgYCAOHjyI7t2747PPPqs0W0VFBfbu3YvZs2fj3XffBQB8+umnaNKkCXQ6HaZOnYqMjIzqvEV1CgvYgimKgoTFC+HWWIXCH76SHYfMxL5jEKDpg9DwCBQWFhqXFxcXQ6vVomXLljh//jzCwsJQUFCAXbt2YdCgQdBqtXjhhRdw7tw54zZRUVFQFAVeXl5o0aIFvLy8YGVlBY1Gg9OnTwMAVqxYgc6dO8PX1xdHjhypdN7X1tYWkZGRAK6dPHR921s99dRTt63zww8/ICYmBgDg6ekJb2/v6r5FdQYL2MLZ2tpi84a1sMvNQOHh72THITNxCHgKebYtMChmKAyGa8eFN2zYEJmZmThz5gyEEJg/fz4MBgMaN26MzMxM49cvv/xiHKdBgwYAACsrK+Ofrz+uqKjAqVOnMGvWLCQnJ+PQoUOIiIio9HA4Gxsb43UrVCoVKioqKs19/Tnutg79Dwu4FnjooYewfeu3KN39fyg5e1h2HDIDRVHg0HM0dh85iYn/uvmDNnt7e8yZMwcfffQR7O3t0a5dO6xcuRLAtXvvHTx4sMrPc/XqVTg4OMDJyQnnz5/H5s2ba/R1AMATTzyBFStWAAB0Oh0OH+bP8HUs4FqiU6dOWL08EQWbZ6H80rl7b0C1nmJtA3XkG1j4xZdYuvTmQ718fX3h7e2NxMRELFu2DAkJCfDx8YFGo8G6deuq/Bw+Pj7w9fVFp06d8Mwzz+CJJ56o6ZeB0aNH48KFC/Dw8MBbb70FjUYDJyenGn+e2oiHodUy8+bPx6RpM+E05ENY2allxyEzKLtwBldXv42tm9abpCBNTa/Xo7y8HHZ2djhx4gRCQ0ORlZVVr66HzAuy1xEvjRmDw0d0WLl5Fhz7vQ3FSiU7EpmYbfO2cOj1CiL7DcCB9L217tZORUVFCAkJQXl5OYQQ+PTTT+tV+d4N94BroYqKCvTs1Qe6Anuoez4vOw6ZSUHGBjQ6swOZ6T+hUaNGsuPQfeCJGHWItbU1NqxZBceLv6LgwCbZcchMHDpH4qpTB/R/ejD0eh4XXhdwCuIObjzt05zeeecddO/eHaGhoXddz8nJCdu3fgu/LoEobtwaDdv53rZOzoJ4WNk2BBQFKocmaBYxDip1k/8tt7r272/TsNGwc3kMRUd348Ka99F65ALYPNQGACCEAZeSP0PJmYMAFCjWtmjWbyJsGres8ddOd6coChyC/4HMde/h5X++hvlzZ9fIuI8//jh27dpVI2PR/WEBW5j33nuvyut26NABG9Z+gz5R/WE96H1jad6oRewHUNk74dKOpbiyZwWahr5w0/IbFf6yEw1cPFCo24nG3YYCAIp+SYM+/0+0ip8HRbFCxdU8KLZ21XiFVB2KyhoOfSfgq6QJ8PJ8DC++8EK1x2T5ysMpiHtITU1FcHAwBg4ciE6dOmHo0KHGSwbu27cPjz/+OHx8fNClSxfk5+ejpKQEI0aMgJeXF3x9fZGSkgIAWLJkCfr374+wsDC4urpi3rx5+Pjjj+Hr64vAwEBcvHgRwLXz91etWgUAcHV1xeTJk9G5c2d4eXkZz/G/cOECwsLCoNFosHTpUjRQAVfWvgd98dU7vg67Np4ov/TbHb9vKCtGac4RPNTnFRT+utO4XF9wESp1UyjKtR8V60bNoOLRF1Kp7NRwjJ6E1ye+afz5qg61+trf591+1sk0WMBVcODAAcyePRs6nQ4nT57Ejz/+iLKyMgwZMgT/+c9/cPDgQWzfvh0NGzbE/PnzoSgKDh8+jMTERAwfPtx4ZtHPP/+Mb775Bvv27cOkSZNgb2+PAwcOICgo6I6X9GvWrBn279+PUaNGGS968u6776Jnz544cuQIBg4ciIsXL+KZgf1RsOlDCH15peMUn9gL2+auxsfnE9/Eb1+MxbkvxwEAio7tgV07P9g0dYbKztF4k1D7Tt1QdHwvfvtiLC5+/1+UnT9RU28rVYNNU2eo+7yO/k8PwtGjR2ts3Mp+1sl0WMBV0KVLF7i4uMDKygparRanT59GVlYWWrVqhYCAAABAo0aNYG1tjR9++AF///vfAVw7eaJt27bGX5CQkBA4OjqiefPmcHJyQlRUFADAy8vrgc+tDw8PR5MmTTB1ymQEPOqCwu8X3bTXcr1oDaVFcAocZFzeIvYDtB4xF62GfQwAKPplJxwe6w4AsH+sO4p0OwBc2+N1/sciNOk+HIpihfNJk1B8OrO6bynVALu23rDpGovQ3n1x6dKlGhmzsp91Mh3OAVfBjefQV+cc91vPxb/xPP3qnluvUqnwzcokdO4ShEsZ6+Dg3x9A5XO9t9IX56PkzCGUXTgNQAGEAYCCxiHxUBQFirUNGnbwR8MO/rCyb4ziY3vQ0FVb1ZdNJqT2CUfBpVz0jR6And9/Bxsbm2qNV1M/61Q13AN+QO7u7jh37hz27dsHAMjPz0dFRQW6deuGZcuWAQCOHj2Ks2fPwt3dvUaf+8Zz67dt22bc+1Gr1di+ZRMMhzag6PhPVR6vKOtHOGhC4DLqC7iM+hwuo5fAunELlOYcQenvx1GR/yeAa0dElF84BRUvEG9RHLrH4WheKZ5/cTTnbGsZFvADsrW1xfLlyzF27Fj4+PggLCwMJSUlGD16NAwGA7y8vDBkyBAsWbLkpr2KmjB58mRs27YNnp6eWLlyJVq2bAlHR0cAwCOPPIJvN6xD0XfzAEPV9l4Kf9lx7RKIN7Dv+DgKdTtgKLqCC6vfw28Jo3Hu85cAKxUa+UXW6Ouh6lGsVHAIH4dvtnyPT2bPkR2H7gPPhKuFSktLoVKpYG1tjd27d2PUqFHGi2tf9/XXiXjhldfQOObfUDk0kROUzKriynlcTpqIFcuWom/fvrLj0A14LYg65OzZsxg8eDAMBgNsbW0rvTvBM8/E4medDgsTZ6LRoA8kpCRzs3ZqAXXERMQMfRaHDmTUumtG1Ecs4Fro0UcfxYEDB+65XkWFHvjrItpUT4hrpygr/HuvFTgHbMHi4+Px8MMPw9PT8763/fLLr7AgYQkc+ow3QTKyROWXf0fBtzOxInEZ2rZte9d133//fWg0Gnh7e0Or1eKnn+78oe2Nt6d3dXVFXl5ejeauz7gHbMHi4uLw0ksvYdiwYfe13a5duzBq7CtwGjgVKofGpglHFsVQWoiC9dMwbcpkhIeH33Xd3bt3Y+PGjdi/fz8aNGiAvLw8lJWV3XH9//73vzUdl/7CPWAL1r17dzRt2vS+tjlz5gwi+w2AutfLN535RnWXMOhR8O0sDIrsjX++Mvae6587dw7NmjUzHp3TrFkztG7dGsnJyfD19YWXlxfi4+NRWloK4Obb01PNYgHXIfn5+ejZqw9U2n5o2CFAdhwyk8Idn8OjlSMWfjqvSuv36tUL2dnZ6NixI0aPHo0dO3agpKQEcXFxWL58OQ4fPoyKigosWLDAxMmJBVxH6PV69H96MC47toN952jZcchMCjM3Q52nw8a1q2FtXbUZRbVajYyMDCxevBjNmzfHkCFDsGjRIrRr1w4dO3YEAAwfPhw7d+68x0hUXZwDriNeGfc69p88j0YDJvMT8Hqi+HQmyvctx/afdqNx48b3ta1KpUJwcDCCg4Ph5eWF+fPnmyYk3RX3gOuARYsX48vlq6GOmABFxX9T64PyP3NQsOVjrF29Em5ubve1bVZWFo4dO2Z8nJmZiQ4dOuD06dM4fvzaVfC++uor9OjRo0Yz0+1YwBYsNjYWQUFByMrKgouLCxISEm5bJzU1Fa9N+Bcco9+CqqGjhJRkbvriq8hfPw2fzJyB4ODg+96+oKAAw4cPh4eHB7y9vaHT6TBjxgx88cUXGDRoELy8vGBlZYUXX3yx5sPTTXgqci12/Phx+HcNQoOwf/LqZPWE0Jfj6pp38ffIEMyd/bHsOFRFvClnHXPp0iU82asPbLrEsHzrCSEEClMWw69DK8z+aKbsOFQDOGFYC5WXlyOy/9MobK6B2ufuB91T3VG0fx2aFJ7FmpTdUKlUsuNQDWAB10KjxozFL+cL4Rj9T9lRyEyKTuyDPnM9tqfvNV56lGo/TkHUMrP/MxcrN26Dus9rUKy4F1QflF04jcJtc7Bp/dp7XuOhOqpz7RF6MCzgWmTLli14a/K7UEdPglUDB9lxyAz0hZeRv/4DLJw3B0FBQffeoBri4uKwZcsWkz4H3YwFXEvodDoMjh0KdcR42DRuKTsOmYGoKEPBphkYPTIOzz77d5M/34Nce4SqhwVcC+Tl5SE0vC8aPDEcdm3438P6QAiBgu3z8YT3o5j+/lTZcchE+CGchSsrK0N4ZD+UtekCteeTsuOQmRTuXYWW4iJWfL0eVlbcT6qrWMAWTAiBuOf+gVP5CtSRz8qOQ2ZSmPUj8Ms2fLc/Hfb29rLjkAnxn1YLNv3Df+Pb1N1w6P1PKAr/quqD0t+PozhlIbZ9uxGtW7eWHYdMjL/VFmrdunX44N8fQR31JqxsG8qOQ2ZQkZ+Hgg3TsTThM3Tu3Nnsz1+Va49QzeIUhAXKzMzE0OEj4Bj9FqwbNZcdh8zAUF6Cgg3TMf6fL+Hpp5+WkiExMVHK89Zn3AO2ML///jt69Y1Ewx7/QIPW7rLjkBkIYUDhtjno9YQf3n5rkuw4ZEYsYAtSXFyMXn2jYHg0GA6PdZcdh8ykcFciHrErxVdLPufF9OsZTkFYCCEEhg6LQ65eDXVgjOw4ZCaFuhRYn/oBW/enG2+SSfUHC9hCvDPlXaTsOwzHp6dxL6ieKM39BSVpX2B32g48/PDDsuOQBJyCsABJScsxe/4iOET+C1Y23AuqDyqu/IH8TR8i8aulvPhNPcYClmzfvn0Y+eIoOEa9CWs1z8OvDwylRSjY8D4mv/kGIiMjZcchiVjAEuXk5CA8Mhr2T46BbYv2suOQGQiDHgVbP0Z0WA+8/to42XFIMs4BS1JSUoLQ8AiUOLSGgzCg6OiuBx4rf/+3KD2XVYPpbmb/aCDsO1b9UoglZw+j5OxhOHaOhMq+kclyAdfuDnxlz8oH2laUFRv/rNxwsouVnRpNn/xHtbPdyq6tD4r2rEDHprZIWLyQc/3EApbF1tYWXl5ecLl8BSg+VK2xks9k1kyoOyg8koLANlW//nByxkYAwMWt8wAAzVu0hLfvbfcjrBHJOzfWyDg3lrG+rBje1fw7udWp40eR88MyPGyvwrf798HGxqZGx6daSghR5S8/Pz9BVFU7d+4UAAQAkZSUJDvOHeXm5gpPT08BQPTp08ckzzF16lRh79hIZGVlmWR8smwA0kUlnco5YDKZbt26QQiB3r17IyYmBidPnpQdqVKtW7fG4cOHsWrVKmzevBnTp0+v8ed49dVXcXB/Bjp27FjjY1PtpVwr56rx9/cX6enpJoxDdZWiKHBxcUF2drbsKHfVo0cP7Ny5E/fze0F0L4qiZAghbpuH4x4wmcXMmTORk5MDg8EgO8pdrVmzBgCQnJwsOQnVByxgMovXXnsNAHD69Gm5Qe7h+j3RFi5cKDkJ1QcsYDKL64dcnTt3TnKSqqktOal2YwGTWXFuleh/WMBERJKwgImIJGEBExFJwgImIpKEBUxEJAkLmIhIEhYwEZEkLGAiIklYwEREkrCAiYgkYQETEUnCAiYikoQFTHSDw4cPAwCys7Mt9g4eVHewgIn+kp2dDW9vbwDA2bNn0a1HMK/eRibFAib6S5s2bfD114mwd3RC81bO+HbjBt46nkyKt6UnukFsbAwCA7vC0dERzZo1kx2H6jgWMNEt2rVrJzsC1ROcgiAikoQFTEQkCQuYapWQkBBs3br1pmWzZ8/GqFGjKl3f1dUVeXl5dx3zgw8+uOnx448/DuDaHZw9PT0BAOnp6Xj55ZcfNLZRXFwcVq1aVe1xqG5gAVOtEhsbi6SkpJuWJSUlITY29oHHvLWAd+3adds6/v7+mDNnzgM/B1FlWMBUqwwcOBCbNm1CWVkZgGt7qb/99htyc3Ph5eUFT09PTJw4sdJt+/fvDz8/P2g0GixevBgA8MYbb6C4uBharRZDhw4FAKjV6tu2TU1NRWRkJACgb9++0Gq10Gq1cHJywtKlS6HX6zF+/HgEBATA29sbixYtAnDtLtAvvfQS3N3dERoaij/++KPG3xOqxYQQVf7y8/MTRA8KgEhLS6v2OBEREWLt2rVCCCGmT58uRowYIdq0aSP++OMPUV5eLkJCQsSaNWuEEEK0bdtWXLhwQQghxJ9//imEEKKoqEhoNBqRl5cnhBDCwcHhpvGvPz516pTQaDRCCCFSUlJERETETeulp6cLLy8vcfnyZbFo0SIxdepUIYQQJSUlws/PT5w8eVKsXr1ahIaGioqKCpGbmyucnJzEypUrq/0eUO0CIF1U0qncA6Za58ZpiKSkJLRt2xbBwcFo3rw5rK2tMXToUOzcufO27ebMmQMfHx8EBgYiOzsbx44de+AMeXl5ePbZZ/H111/DyckJ27Ztw5dffgmtVouuXbvizz//xLFjx7Bz507ExsZCpVKhdevW6Nmz5wM/J9U9LGCqdfr164fk5GTs378fRUVF0Gq199wmNTUV27dvx+7du3Hw4EH4+vqipKTkgZ5fr9cjJiYG77zzjvFDOiEE5s6di8zMTGRmZuLUqVPo1avXA41P9QcLmGodtVqNkJAQxMfHIzY2Fl26dMGOHTuQl5cHvV6PxMRE9OjR46Ztrly5giZNmsDe3h6//vor9uzZY/yejY0NysvLq/z8b7zxBry9vRETE2Nc1rt3byxYsMA4ztGjR1FYWIju3btj+fLl0Ov1OHfuHFJSUqr56qku4ZlwVCvFxsZiwIABSEpKQqtWrTBjxgyEhIRACIGIiAj069fvpvXDw8OxcOFCPPbYY3B3d0dgYKDxe88//zy8vb3RuXNnLFu27J7PPWvWLGg0GuOe93vvvYeRI0fi9OnT6Ny5M4QQaN68OdauXYsBAwbg+++/h4eHBx555BEEBQXV6PtAtZsi7uNqT/7+/iI9Pd2EcaguUxQFaWlp+Nvf/iY7CpFZKYqSIYTwv3U5pyCIiCRhARMRScICJiKShAVMdIMtW7bA3d0dbm5umDFjhuw4VMexgIn+otfrMWbMGGzevBk6nQ6JiYnQ6XSyY1EdxgIm+svevXvh5uaG9u3bw9bWFjExMVi3bp3sWFSHsYCJ/pKbm4s2bdoYH7u4uCA3N1diIqrrWMBERJKwgIn+4uzsjOzsbOPjnJwcODs7S0xEdR0LmOgvAQEBOHbsGE6dOoWysjIkJSUhOjpadiyqw3gtCKK/WFtbY968eejduzf0ej3i4+Oh0Whkx6I6jAVMJnf94tMAYDAYIISAoiiSU1Wub9++6Nu3r+wYVE9wCoJM7vMlS6BSqQAAPXr0QFz8SMmJiCwDC5hMbuBTT6FbSCgcHm4Dl7bt8eYbE2RHIrIInIIgk3NycsL32zYjMTERERERaNq0qexIRBaBBUxmYW1tjWeffVZ2DCKLwikIIiJJWMBERJKwgImIJGEB031RqVTQarXw9PREVFQULl++/MBjXblyBcOGDYObmxs6dOiAYcOG4cqVK8bvjx8/HhqNBuPHj8eFCxfQtWtX+Pr6Ii0trQZeCZF8LGC6Lw0bNkRmZiZ+/vlnNG3aFPPnz3/gsZ577jm0b98ex48fx4kTJ9CuXTuMHPm/Y4QXL16MQ4cOYebMmUhOToaXlxcOHDiAbt261cRLIZKOBUwPLCgoyHi5xuDgYFy/Y3ZeXh5cXV0BAEuWLMFTTz2F8PBwPProo5gw4doxwMePH0dGRgbefvtt43jvvPMO0tPTceLECURHR6OgoAB+fn748MMPMWHCBKxbtw5arRbFxcXmfaFEJsLD0OiB6PV6JCcn47nnnrvnupmZmThw4AAaNGgAd3d3jB07FjqdDlqt1niGHPC/6Y0jR45g/fr1UKvVyMzMBAC0aNEC6enpmDdvnqleEpHZcQ+Y7ktxcTG0Wi1atmyJ8+fPIyws7J7bPPnkk3BycoKdnR08PDxw5swZMyQlsnwsYLov1+eAz5w5AyGEcQ7Y2toaBoMBAFBSUnLTNg0aNDD+WaVSoaKiAh4eHsjMzDRuA1y7UE9mZiY8PDzM8EqI5GMB0wOxt7fHnDlz8NFHH6GiogKurq7IyMgAAKxateqe27u5ucHX1xfTpk0zLps2bRo6d+4MNzc3k+UmsiQsYHpgvr6+8Pb2RmJiIl5//XUsWLAAvr6+yMvLq9L2CQkJOHr0KDp06IAOHTrg6NGjSEhIMHFqIsuhXL9Oa1X4+/uL6590ExFR1SiKkiGE8L91OfeAiYgkYQETEUnCAiaz2LJlC9zd3eHm5oYZM2bIjkNkEVjAZHJ6vR5jxozB5s2bodPpkJiYCJ1OJzsWkXQsYDK5vXv3ws3NDe3bt4etrS1iYmKwbt062bGIpGMBk8nl5uaiTZs2xscuLi7Ga0gQ1WcsYCIiSVjAZHLOzs7Izs42Ps7JyYGzs7PERESWgQVMJhcQEIBjx47h1KlTKCsrQ1JSEqKjo2XHIpKOl6Mkk7O2tsa8efPQu3dv6PV6xMfHQ6PRyI5FJB1PRSYiMjGeikxEZGFYwEREkrCAiYgkYQETEUnCAiYikoQFTEQkCQuYiEgSFjARkSQsYCIiSVjARESSsICJiCRhARMRScICJiKShAVMRCQJC5iISBIWMBGRJCxgIiJJ7uuOGIqiXABwxnRxiIjqpLZCiOa3LryvAiYioprDKQgiIklYwEREkrCAiYgkYQETEUnCAiYikoQFTEQkCQuYiEgSFjARkSQsYCIiSf4f4E203kIcYBIAAAAASUVORK5CYII=\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.7.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
