{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Wastewater Treatment (WWT)"
   ]
  },
  {
   "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",
    "The WWT process models in the SAT framework are designed so that by changing default parameters, a variety of WWT technology components and configurations can be represented. A relatively small amount of the incoming PFAS may be released in initial screen rejects and grit. Volatilization is another potential pathway, and its importance is currently poorly understood. Most PFAS entering conventional WWTPs will exit in the effluent or with the solids. In an increasing number of cases, some fraction of the PFAS are removed in tertiary treatment through reverse osmosis (RO), ion exchange (IX) or granular activated carbon (GAC) adsorption systems prior to release. These systems then produce new waste streams (i.e., RO concentrate and spent GAC, respectively) that must be managed.\n",
    "\n",
    "The mass flow of PFAS through the WWT processes is modeled using a water and mass balance through the system. The incoming PFAS-waste is diluted with the rest of the influent wastewater. When PFAS is partitioned between liquids and solids, it is assumed that equilibrium is achieved. The model does not consider transformations of PFAS. A relatively small fraction of the influent water leaves in the screen rejects and grit. The water lost with the material is assumed to have the same PFAS concentration as the influent. The remainder of the water and PFAS enters primary settling (if the WWTP has primary settling), and a fraction of the water at the influent PFAS concentration is again lost. The rest of the PFAS-containing water enters biological treatment, where PFAS partitions between the solids and liquid based on the organic carbon-normalized partition coefficient for each PFAS. The effluent from biological treatment enters secondary settling, which acts in a similar manner to primary settling. The effluent from secondary settling can then either be released to surface water or continue to tertiary treatment. The primary and secondary solids can then be thickened, dewatered, dried, stabilized, anaerobically digested or composted. The effluent can be treated through tertiary treatment processes designed to remove PFAS, such as GAC, RO, or IX.\n",
    "\n",
    "<img src=\"../Images/ProcessModels/WWTP.png\" alt=\"Drawing\" style=\"width: 700px;\"/>\n",
    "\n",
    "\n",
    "<img src=\"../Images/ProcessModels/WWT_Settling_Solids_Diagram_v3.png\" alt=\"Drawing\" style=\"width: 700px;\"/>\n",
    "\n",
    "\n",
    "\n",
    "### Assumptions and Limitations\n",
    "\n",
    "1.\tThe model assumes that liquids lost in screen rejects, grit, and settled solids have the same PFAS concentrations as the liquid entering the process.\n",
    "2.\tVolatilization is assumed to be zero. However, the user may assign a fraction of the PFAS that volatilizes/aerosolizes.\n",
    "3.\tFuture work and additional data are required to include PFAS transformations during WWT processes. \n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for WWT model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Category</th>\n",
       "      <th>Parameter Name</th>\n",
       "      <th>Parameter Description</th>\n",
       "      <th>amount</th>\n",
       "      <th>unit</th>\n",
       "      <th>minimum</th>\n",
       "      <th>maximum</th>\n",
       "      <th>Reference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>WWT Plant</td>\n",
       "      <td>Des_Cap</td>\n",
       "      <td>Design capacity of the WWT plant</td>\n",
       "      <td>7.6000</td>\n",
       "      <td>Million Liter/day</td>\n",
       "      <td>3.00</td>\n",
       "      <td>300.000</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Screen Rejects</td>\n",
       "      <td>frac_sr-grit</td>\n",
       "      <td>Mass fraction influent into screen rejects and grit</td>\n",
       "      <td>0.0001</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Screen Rejects</td>\n",
       "      <td>sol_cont_sr_grit</td>\n",
       "      <td>Solids content of the screen rejects and grit</td>\n",
       "      <td>0.6500</td>\n",
       "      <td>fraction</td>\n",
       "      <td>0.50</td>\n",
       "      <td>0.900</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Primary Settling</td>\n",
       "      <td>is_prim_set</td>\n",
       "      <td>Is there primary settling</td>\n",
       "      <td>1.0000</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Primary Settling</td>\n",
       "      <td>frac_prim_solids</td>\n",
       "      <td>Mass fraction influent into settling that exits in the solids</td>\n",
       "      <td>0.0040</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Primary Settling</td>\n",
       "      <td>sol_cont_prim_solids</td>\n",
       "      <td>Solids content of primary solids - wet</td>\n",
       "      <td>0.0600</td>\n",
       "      <td>fraction</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.080</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Biological Treatmet</td>\n",
       "      <td>sol_cont</td>\n",
       "      <td>Mixed liquor suspended solids (MLSS)</td>\n",
       "      <td>3000.0000</td>\n",
       "      <td>mg/L</td>\n",
       "      <td>2000.00</td>\n",
       "      <td>4000.000</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Biological Treatmet</td>\n",
       "      <td>C_cont</td>\n",
       "      <td>Carbon content of the MLSS</td>\n",
       "      <td>0.5000</td>\n",
       "      <td>kg TS/kg</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Biological Treatmet</td>\n",
       "      <td>VS_cont</td>\n",
       "      <td>VS content of the MLSS - dry</td>\n",
       "      <td>0.8000</td>\n",
       "      <td>kg VSS/kg TS</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Thickening</td>\n",
       "      <td>is_prim_thick</td>\n",
       "      <td>Does primary sludge go to thickening</td>\n",
       "      <td>1.0000</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>Thickening</td>\n",
       "      <td>is_sec_thick</td>\n",
       "      <td>Does secondary sludge go to thickening</td>\n",
       "      <td>1.0000</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Thickening</td>\n",
       "      <td>sol_cont_thick</td>\n",
       "      <td>Solids content after thickening</td>\n",
       "      <td>0.0500</td>\n",
       "      <td>fraction</td>\n",
       "      <td>0.03</td>\n",
       "      <td>0.100</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Dewatering</td>\n",
       "      <td>is_sol_dew</td>\n",
       "      <td>Are solids dewatered</td>\n",
       "      <td>1.0000</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>Dewatering</td>\n",
       "      <td>sol_cont_dewat</td>\n",
       "      <td>Solids content of dewatered solids</td>\n",
       "      <td>0.2200</td>\n",
       "      <td>kg TS/kg</td>\n",
       "      <td>0.15</td>\n",
       "      <td>0.300</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Secondary Settling</td>\n",
       "      <td>is_sec_set</td>\n",
       "      <td>Is there secondary settling</td>\n",
       "      <td>1.0000</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Secondary Settling</td>\n",
       "      <td>frac_sec_solids</td>\n",
       "      <td>Mass fraction influent into settling that exits in the solids</td>\n",
       "      <td>0.0160</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Drying</td>\n",
       "      <td>is_sol_dry</td>\n",
       "      <td>Are solids dried</td>\n",
       "      <td>1.0000</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Drying</td>\n",
       "      <td>sol_cont_dry</td>\n",
       "      <td>Solids content of dried solids</td>\n",
       "      <td>0.9000</td>\n",
       "      <td>kg TS/kg</td>\n",
       "      <td>0.65</td>\n",
       "      <td>0.950</td>\n",
       "      <td>[1]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>Drying</td>\n",
       "      <td>frac_PFAS_to_Vol</td>\n",
       "      <td>Fraction of PFAS lost to volatilization</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>Volatilization</td>\n",
       "      <td>frac_vol_loss</td>\n",
       "      <td>Fraction of PFAS lost to volatilization</td>\n",
       "      <td>0.0000</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log Koc (WWT)</td>\n",
       "      <td>2.1900</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>1.30</td>\n",
       "      <td>4.500</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log Koc (WWT)</td>\n",
       "      <td>3.0400</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.40</td>\n",
       "      <td>4.700</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log Koc (WWT)</td>\n",
       "      <td>1.8800</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>1.30</td>\n",
       "      <td>1.880</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log Koc (WWT)</td>\n",
       "      <td>1.3700</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log Koc (WWT)</td>\n",
       "      <td>1.7700</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>1.31</td>\n",
       "      <td>2.100</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log Koc (WWT)</td>\n",
       "      <td>1.9700</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>2.190</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log Koc (WWT)</td>\n",
       "      <td>2.6300</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.30</td>\n",
       "      <td>3.180</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log Koc (WWT)</td>\n",
       "      <td>3.2400</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.65</td>\n",
       "      <td>3.780</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log Koc (WWT)</td>\n",
       "      <td>1.5100</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1.790</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log Koc (WWT)</td>\n",
       "      <td>2.7900</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>2.05</td>\n",
       "      <td>2.875</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                     Category        Parameter Name  \\\n",
       "0   WWT Plant                  Des_Cap                \n",
       "1   Screen Rejects             frac_sr-grit           \n",
       "2   Screen Rejects             sol_cont_sr_grit       \n",
       "3   Primary Settling           is_prim_set            \n",
       "4   Primary Settling           frac_prim_solids       \n",
       "5   Primary Settling           sol_cont_prim_solids   \n",
       "6   Biological Treatmet        sol_cont               \n",
       "7   Biological Treatmet        C_cont                 \n",
       "8   Biological Treatmet        VS_cont                \n",
       "9   Thickening                 is_prim_thick          \n",
       "10  Thickening                 is_sec_thick           \n",
       "11  Thickening                 sol_cont_thick         \n",
       "12  Dewatering                 is_sol_dew             \n",
       "13  Dewatering                 sol_cont_dewat         \n",
       "14  Secondary Settling         is_sec_set             \n",
       "15  Secondary Settling         frac_sec_solids        \n",
       "16  Drying                     is_sol_dry             \n",
       "17  Drying                     sol_cont_dry           \n",
       "18  Drying                     frac_PFAS_to_Vol       \n",
       "19  Volatilization             frac_vol_loss          \n",
       "20  Log partition coefficient  PFOA                   \n",
       "21  Log partition coefficient  PFOS                   \n",
       "22  Log partition coefficient  PFBA                   \n",
       "23  Log partition coefficient  PFPeA                  \n",
       "24  Log partition coefficient  PFHxA                  \n",
       "25  Log partition coefficient  PFHpA                  \n",
       "26  Log partition coefficient  PFNA                   \n",
       "27  Log partition coefficient  PFDA                   \n",
       "28  Log partition coefficient  PFBS                   \n",
       "29  Log partition coefficient  PFHxS                  \n",
       "\n",
       "                                            Parameter Description     amount  \\\n",
       "0   Design capacity of the WWT plant                               7.6000      \n",
       "1   Mass fraction influent into screen rejects and grit            0.0001      \n",
       "2   Solids content of the screen rejects and grit                  0.6500      \n",
       "3   Is there primary settling                                      1.0000      \n",
       "4   Mass fraction influent into settling that exits in the solids  0.0040      \n",
       "5   Solids content of primary solids - wet                         0.0600      \n",
       "6   Mixed liquor suspended solids (MLSS)                           3000.0000   \n",
       "7   Carbon content of the MLSS                                     0.5000      \n",
       "8   VS content of the MLSS - dry                                   0.8000      \n",
       "9   Does primary sludge go to thickening                           1.0000      \n",
       "10  Does secondary sludge go to thickening                         1.0000      \n",
       "11  Solids content after thickening                                0.0500      \n",
       "12  Are solids dewatered                                           1.0000      \n",
       "13  Solids content of dewatered solids                             0.2200      \n",
       "14  Is there secondary settling                                    1.0000      \n",
       "15  Mass fraction influent into settling that exits in the solids  0.0160      \n",
       "16  Are solids dried                                               1.0000      \n",
       "17  Solids content of dried solids                                 0.9000      \n",
       "18  Fraction of PFAS lost to volatilization                        0.0000      \n",
       "19  Fraction of PFAS lost to volatilization                        0.0000      \n",
       "20  PFOA Log Koc (WWT)                                             2.1900      \n",
       "21  PFOS Log Koc (WWT)                                             3.0400      \n",
       "22  PFBA Log Koc (WWT)                                             1.8800      \n",
       "23  PFPeA Log Koc (WWT)                                            1.3700      \n",
       "24  PFHxA Log Koc (WWT)                                            1.7700      \n",
       "25  PFHpA Log Koc (WWT)                                            1.9700      \n",
       "26  PFNA Log Koc (WWT)                                             2.6300      \n",
       "27  PFDA Log Koc (WWT)                                             3.2400      \n",
       "28  PFBS Log Koc (WWT)                                             1.5100      \n",
       "29  PFHxS Log Koc (WWT)                                            2.7900      \n",
       "\n",
       "                 unit  minimum   maximum Reference  \n",
       "0   Million Liter/day  3.00     300.000             \n",
       "1   fraction          NaN      NaN                  \n",
       "2   fraction           0.50     0.900     [1]       \n",
       "3   1:TRUE,0:FALSE    NaN      NaN                  \n",
       "4   fraction          NaN      NaN        [1]       \n",
       "5   fraction           0.04     0.080     [1]       \n",
       "6   mg/L               2000.00  4000.000  [1]       \n",
       "7   kg TS/kg          NaN      NaN                  \n",
       "8   kg VSS/kg TS      NaN      NaN                  \n",
       "9   1:TRUE,0:FALSE    NaN      NaN                  \n",
       "10  1:TRUE,0:FALSE    NaN      NaN                  \n",
       "11  fraction           0.03     0.100     [1]       \n",
       "12  1:TRUE,0:FALSE    NaN      NaN                  \n",
       "13  kg TS/kg           0.15     0.300     [1]       \n",
       "14  1:TRUE,0:FALSE    NaN      NaN                  \n",
       "15  fraction          NaN      NaN        [1]       \n",
       "16  1:TRUE,0:FALSE    NaN      NaN                  \n",
       "17  kg TS/kg           0.65     0.950     [1]       \n",
       "18  fraction          NaN      NaN                  \n",
       "19  fraction          NaN      NaN                  \n",
       "20  log L/kg OC        1.30     4.500               \n",
       "21  log L/kg OC        2.40     4.700               \n",
       "22  log L/kg OC        1.30     1.880               \n",
       "23  log L/kg OC       NaN      NaN                  \n",
       "24  log L/kg OC        1.31     2.100               \n",
       "25  log L/kg OC       NaN       2.190               \n",
       "26  log L/kg OC        2.30     3.180               \n",
       "27  log L/kg OC        2.65     3.780               \n",
       "28  log L/kg OC       NaN       1.790               \n",
       "29  log L/kg OC        2.05     2.875               "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "WWT = pspd.WWT()\n",
    "WWT.InputData.Data[['Category','Parameter Name', 'Parameter Description', 'amount', 'unit','minimum','maximum','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming Landfill Leachte to WWT"
   ]
  },
  {
   "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>10.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>990.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Volume flow</td>\n",
       "      <td>L</td>\n",
       "      <td>1000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>5.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>5700.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>90.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>750.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>680.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1650.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>550.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>50.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>30.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>190.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>270.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   10.0  \n",
       "2   Moisture flow  kg   990.0 \n",
       "3   Volume flow    L    1000.0\n",
       "4   Carbon flow    kg   5.0   \n",
       "5   PFOA           μg   5700.0\n",
       "6   PFOS           μg   90.0  \n",
       "7   PFBA           μg   750.0 \n",
       "8   PFPeA          μg   680.0 \n",
       "9   PFHxA          μg   1650.0\n",
       "10  PFHpA          μg   550.0 \n",
       "11  PFNA           μg   50.0  \n",
       "12  PFDA           μg   30.0  \n",
       "13  PFBS           μg   190.0 \n",
       "14  PFHxS          μg   270.0 "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncomingWaste = pspd.IncomFlow()\n",
    "IncomingWaste.set_flow('LFLeachate', 1000)\n",
    "IncomingWaste.calc()\n",
    "LFLeachate = IncomingWaste.Inc_flow\n",
    "LFLeachate.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in WWT"
   ]
  },
  {
   "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>DryerExhaust</th>\n",
       "      <th>WWT Effluent</th>\n",
       "      <th>solids</th>\n",
       "      <th>Screen Rejects</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>PFOA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>80.21</td>\n",
       "      <td>19.78</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>37.37</td>\n",
       "      <td>62.62</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>88.74</td>\n",
       "      <td>11.25</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>95.48</td>\n",
       "      <td>4.51</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>90.82</td>\n",
       "      <td>9.17</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>86.70</td>\n",
       "      <td>13.29</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>60.27</td>\n",
       "      <td>39.72</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>27.40</td>\n",
       "      <td>72.59</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>94.26</td>\n",
       "      <td>5.73</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>51.35</td>\n",
       "      <td>48.64</td>\n",
       "      <td>0.01</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Volatilized  DryerExhaust  WWT Effluent  solids  Screen Rejects\n",
       "PFOA   0.0          0.0           80.21         19.78   0.01          \n",
       "PFOS   0.0          0.0           37.37         62.62   0.01          \n",
       "PFBA   0.0          0.0           88.74         11.25   0.01          \n",
       "PFPeA  0.0          0.0           95.48         4.51    0.01          \n",
       "PFHxA  0.0          0.0           90.82         9.17    0.01          \n",
       "PFHpA  0.0          0.0           86.70         13.29   0.01          \n",
       "PFNA   0.0          0.0           60.27         39.72   0.01          \n",
       "PFDA   0.0          0.0           27.40         72.59   0.01          \n",
       "PFBS   0.0          0.0           94.26         5.73    0.01          \n",
       "PFHxS  0.0          0.0           51.35         48.64   0.01          "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "WWT.calc(Inc_flow=LFLeachate)\n",
    "WWT.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAu0UlEQVR4nO3deVhV1f7H8ffhHFAZxDkVTDSUEDhM4ng1ccIcUHPOITWz1GwmbVJvZtLNyhQjNUtNA1NzzJzHmxQo4oQpKk5oJQ7IPK7fH1zPL1IQDNgM39fz8Dydfdbe+7tP8nk26+y1lk4phRBCiNJnpnUBQghRWUkACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoxFCUxnXq1FEODg4lVIoQQlRMhw8fjldK1f379iIFsIODA4cOHSq+qoQQohLQ6XQX77dduiCEEEIjEsBCCKERCWAhhNCIBLAoE/R6PR4eHri4uODu7s4nn3xCTk5OiZ/v7k9gYGCB7a2trUusFoC5c+eSkpJSoucQZU+RvoQToqRUq1aNqKgoAP7880+efvpp7ty5w7///e887bKysjAYHv6f7d39/3q+smDu3LmMGDECS0tLrUsRpUjugEWZU69ePRYtWkRQUBBKKZYuXYq/vz+dO3emS5cujBo1ivXr15vaDx8+nA0bNpCdnU1AQAA+Pj4YjUYWLlwIwN69e+nQoQP+/v60aNEi3/MmJCTg5OTE6dOnARg2bBiLFy82vf/OO+/g7u5OmzZt+OOPPwDYtGkTrVu3xtPTk65du5q2z5gxgzlz5pj2dXV15cKFCyQnJ9OrVy/c3d1xdXVl1apVzJs3j6tXr+Lr64uvr2+xfY6iHFBKFfrH29tbCVESrKys7tlma2urfv/9d/XNN98oOzs7dePGDaWUUnv37lV9+/ZVSil1+/Zt5eDgoDIzM9XChQvVzJkzlVJKpaWlKW9vb3X+/Hm1Z88eZWlpqc6fP286tpmZmXJ3dzf9hIaGKqWU2r59u2rTpo0KCQlRfn5+pvaA2rhxo1JKqYCAANN5bt68qXJycpRSSi1evFi99tprSimlpk+frj7++GPT/i4uLio2NlatWbNGjRs3zrT99u3bSimlGjdurK5fv/4PPkFRlgGH1H0yVbogRLnQrVs3atWqBcATTzzBxIkTuX79OmvXrmXAgAEYDAa2b9/OsWPHWLNmDZB7RxsTE4OFhQWtWrWiSZMmpuPl1wXRrVs3Vq9ezaRJkzh69Khpu4WFBb179wbA29ubHTt2AHDlyhWGDBnCtWvXyMjIyHOO+3Fzc+P1119nypQp9O7dmw4dOvyjz0WUb9IFIcqk8+fPo9frqVevHgBWVlZ53h81ahQrVqzgm2++YezYsUDuX3Pz588nKiqKqKgoYmNj6d69+333z09OTg6nTp3C0tKSW7dumbabm5uj0+mA3C/wsrKyAJg8eTIvvvgix48fZ+HChaSlpQFgMBjyfIl4d3vz5s2JjIzEzc2Nd999l/fff7/In42oOCSARZlz/fp1XnjhBV588UVT6P3d6NGjmTt3LoCpX9fPz4/g4GAyMzMBOHPmDMnJyUU692effYazszPfffcdY8aMMR0rPwkJCdjZ2QGwbNky03YHBwciIyMBiIyMJDY2FoCrV69iaWnJiBEjCAgIMLWxsbEhMTGxSLWK8k+6IESZkJqaioeHB5mZmRgMBkaOHMlrr72Wb/tHHnkEZ2dn+vXrZ9o2btw4Lly4gJeXF0op6tatm+fLuvud764ePXowZswYvvrqK8LDw7GxsaFjx4588MEH9zyJ8VczZsxg0KBB1KxZk86dO5uCdsCAASxfvhwXFxdat25N8+bNATh+/DgBAQGYmZlhbm5OcHAwAOPHj6dHjx40bNiQPXv2FPJTE+WdThVhTbiWLVsqmQtClAUpKSm4ubkRGRmJra2t1uUIUSCdTndYKdXy79ulC0KUOzt37sTZ2ZnJkydL+IpyTbogRLnTtWtXLl687+RSQpQrcgcsKr2tW7fi5OSEo6PjA4ckC1GcJIBFpZadnc2kSZP46aefiI6OJiQkhOjoaK3LEpWEBLCo1MLDw3F0dKRp06ZYWFgwdOhQNmzYoHVZopKQABaVWlxcHI0aNTK9tre3Jy4uTsOKRGUiASyEEBqRABaVmp2dHZcvXza9vnLlimlkmxAlTQJYVGo+Pj7ExMQQGxtLRkYGoaGh+Pv7a12WqCTkOWBRqRkMBoKCgvDz8yM7O5uxY8fi4uKidVmikpChyEIIUcJkKLIQQpQxEsCiUGbNmoWLiwtGoxEPDw9+/fVXrUsSotyTPmDxQGFhYWzevJnIyEiqVKlCfHw8GRkZD9zvny6gKURFJ3fA4oGuXbtGnTp1qFKlCgB16tShYcOGRERE0K5dO9zd3WnVqhWJiYn3LKCZnJzM2LFjadWqFZ6enqZRZgUtoNmpUycGDhzI448/zvDhwynK9xRClCv3Wyguvx9ZlLNySkxMVO7u7qpZs2ZqwoQJau/evSo9PV01adJEhYeHK6WUSkhIUJmZmfcsoPnWW2+pb7/9Viml1K1bt1SzZs1UUlJSgQtoVq9eXV2+fFllZ2erNm3aqAMHDpTo9aWnp6uvv/5a3bx5s0TPIyov8lmUU+6AxQNZW1tz+PBhFi1aRN26dRkyZAgLFy6kQYMG+Pj4AFC9enVTd8NfF9Dcvn07gYGBeHh40KlTJ9LS0rh06RLbt29n+fLleHh40Lp1a27cuEFMTAwArVq1wt7eHjMzMzw8PLhw4UKJXdvNmzfp4NuVSQHv4u7lY6pBiNIgHXSiUPR6PZ06daJTp064ubmxYMGCfNv+dQFMpRRr167FyckpTxv1vwU0/fz88mzfu3evqavj7nnvLoBZErb89BNHDv1KZkYGl29A8KLFfPrxf0rsfEL8ldwBiwc6ffp0njvDqKgonJ2duXbtGhEREQAkJibeNyj9/PyYP3++qR/3yJEjpu3/dAHN4jBi+HAy0tMBOHr0qISvKFVyByweKCkpicmTJ3P79m0MBgOOjo4sWrSIMWPGMHnyZFJTU6lWrRo7d+68Z9/33nuPV155BaPRSE5ODk2aNGHz5s1FWkBTiIpKRsIJAeh0Oo4ePYrRaNS6FFEByUg4IYQoYySAhRBCIxLAokB6vR4PDw9cXV3p06cPt2/fLrD9l19+yfLly4t8ntu3b/PFF188ZJVClE8SwKJA1apVIyoqihMnTlCrVq0CHz8DeOGFFxg1alSRzyMBLCojCWBRaG3btjWtl3bu3Dl69OiBt7c3HTp04LfffgNgxowZzJkzp8A2f/zxB/3798fd3R13d3cOHjzI1KlTOXfuHB4eHgQEBHDt2jU6duxouvs+cOCANhctRAmSx9BEoWRnZ7Nr1y6effZZAMaPH8+XX35Js2bN+PXXX5k4cSK7d+/Os09+bV566SWeeOIJ1q1bR3Z2NklJSQQGBnLixAmioqIA+OSTT/Dz8+Odd94hOzublJSU0r5kIUqcBLAoUGpqKh4eHsTFxeHs7Ey3bt1ISkri4MGDDBo0yNQu/X+DGe4qqM3u3btN/cR6vR5bW1tu3bqVZ38fHx/Gjh1LZmYm/fr1w8PDo4SuUAjtSBeEKNDdPuCLFy+ilGLBggXk5ORQo0YNoqKiTD+nTp3Ks19h2hSkY8eO7N+/Hzs7O0aPHv1QX+wJUdZJAItCsbS0ZN68eXzyySdYWlrSpEkTVq9eDeTO63D06NE87atXr55vmy5duhAcHAzkdm0kJCRgY2NDYmKiaf+LFy/yyCOP8NxzzzFu3DgiIyNL4zKFKFUSwKLQPD09MRqNhISEsHLlSpYsWYK7uzsuLi6meX4hd1QZkG+bzz//nD179uDm5oa3tzfR0dHUrl2b9u3b4+rqSkBAAHv37sXd3R1PT09WrVrFyy+/rMk1C1GSZCiyKFaTJ0/Gy8uLMWPGaF1KkchQZFGSZCiyKHHvvfcev/76K/7+/lqXIkS5IAEsis3MmTMJDw+ndu3aWpciRLkgASyEEBqRABYPtHXrVpycnHB0dCQwMPCe99PT0xkyZAiOjo60bt06zxJCs2fPxtHREScnJ7Zt2/bAYwYFBeHo6IhOpyM+Pr5Er0sIzd1vobj8fmRRzsonKytLNW3aVJ07d06lp6cro9GoTp48mafNggUL1PPPP6+UUiokJEQNHjxYKaXUyZMnldFoVGlpaer8+fOqadOmKisrq8BjRkZGqtjYWNW4cWN1/fr1UrtOQB09erTUzicqF2RRTvEwwsPDcXR0pGnTplhYWDB06NA8j5wBbNiwgWeeeQaAgQMHsmvXLpRSbNiwgaFDh1KlShWaNGmCo6Mj4eHhBR7T09MTBweH0r5MITQhASwKFBcXR6NGjUyv7e3tTRPy3K+NwWDA1taWGzdu5LtvYY4pRGUgASyEEBqRABYFsrOz4/Lly6bXV65cwc7OLt82WVlZJCQkULt27Xz3LcwxhagMJIBFgXx8fIiJiSE2NpaMjAxCQ0PvGWjh7+/PsmXLAFizZg2dO3dGp9Ph7+9PaGgo6enpxMbGEhMTQ6tWrQp1TCEqAwlgUSCDwUBQUBB+fn44OzszePBgXFxcmDZtGhs3bgTg2Wef5caNGzg6OvLpp5+aHitzdnamZq3a1K1XDz8/PxYsWIBer8/3mADz5s3D3t6eK1euYDQaGTdunGbXLkRJk7kgRIlQStG771Ns2Zwb0p4tW3Fgz06srKw0ruz+ZC4IUZLymwtCJmQXJUKn0/Fo40dB5QDwmONjVK1aVeOqhChbpAtClJjg+Z/TpEkTAFZ/twK9Xq9xRUKULRLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASxKVFGmOxWispEAFiXqxo0bWpdQaDVq1NC6BFHJSACLEpWYmEinTp20LqNAISEhQO7qzEKUJglgUWL27dsHwOrVqzWupGBPP/00AGZm8usgSpf8ixMlIjk52XTnW6dOHW2LKcDMmTMBiIiI0KyGO3fuSF95JSUBLIrd4sWLsba2BiAtLU3jau4vPj4ed3d3pk2bxttvv03Llvcs11Vq6j1Sn1mzAzU7v9COrAmngX37D7B8xcp/fJyvFy8shmryV79BQ3r27lPo9n+v55kxzzL55VeLu6x/ZMniheTk5JheW1pacv3GLca/MFGzmtLTUgmc8xkuzo/Tv39/zeoQpU8CWAOnT5/mmyVfYfuv4ZhVsXzo41g88hgZf5wrxsrySqzZnPUXdEXaR2ewoEaHEaA3Z+NlgJwH7VKqqjTxJvVcBJbN21HlUTcA1pzTtsa6T72L3roWI8c8ywEHBzw9PTWtR5QeWZZeI+PGv8AP+49g4/8OOjNZrFJA8qkDEP4txyIP0aBBA63LEcUov2XppQ9YI19+EUSL+tYk7/9G61JEGWHl3AHVvDPde/YhNTVV63JEKZAA1ojBYGDz+rVY/XmC5KNbtS5HlBFWbYZwNac6Q0eMkicjKgEJYA3VqFGDXdt/IjM8lNQLUVqXI8oAnU6HVdcX2X84mnenTde6HFHCJIA15ujoyPq1q0ne+imZN65oXY4oA8zMq2Dd+y3mBS/mu+9CtC5HlCAJ4DKgU6dOfPpxIImbZpGdmqh1OaIM0FvXxMb/HcZPmMSvv/6qdTmihEgAlxHjn3uOZ4YOJHnLf1DZWVqXI8oAi3pNsOz6Ij379OXy5ctalyNKgARwGTL3k4/xalqf5D2L5AsYAYBls9aYufWii19PkpKSHvo4r776KnPnzjW99vPzY9y4cabXr7/+OjqdjvXr15u2OTk58cEHH5heDxgwgFWrVuHh4YGHhwfW1tY4OTnh4eHBqFGjHrq2ykwCuAzR6/WsW7OKGkkXSIncqHU5ooywbNmfG1XtGDB4WJ5RfEXRvn17Dh48CEBOTg7x8fGcPHnS9P7BgweZOXOmqc2NGzewsrIiLCzM1CYsLIwnnniCqKgooqKiaNmyJStXriQqKorly5f/gyusvCSAyxgbGxt2bdtCdtQGUs9pN0GMKDt0Oh3WXV4g/PRlXg9486GO0a5dO1OYnjx5EldXV2xsbLh16xbp6emcOnWKjh07mgL44MGD9OnTh+vXr6OUIjY2lmrVqlG/fv1iuy4hAVwmNW7cmB83rid5x3wyrl/QuhxRBuj05lj3epMlK1bx9ddFH7zTsGFDDAYDly5d4uDBg7Rt25bWrVsTFhbGoUOHcHNzo3Xr1pw4cYKMjAxTGycnJ06dOsXBgwdp165dCVxZ5SYBXEa1bduW4Pmfk7jxQ7KTb2tdjigD9Ja22PR5h5dee4MDBw4Uef927dpx8OBBU7i2bdvW9Lp9+/ZUqVIFFxcXIiMj+eWXX2jduvU9bUTxkgAuw0aOHMHEcaNJ+jEQlZWpdTmiDDCv0wgrv1fw7z+A8+fPF2nfu/3Ax48fx9XVlTZt2hAWFpbn7rZ9+/bs37+fxMREatasSZs2bUwBLHfAxU8CuIybPWsm7Y3NSNq5QJ6MEABUa+KFwXsgXbo/SUJCQqH3a9euHZs3b6ZWrVro9Xpq1arF7du3CQsLM4Vru3btWLhwIe7u7gAYjUZ++eUXLl26hKura4lcT2UmAVzGmZmZ8f13K3gkJ57k8DValyPKCCvPXiTWcqJP/4FkZRXuuXE3Nzfi4+Np06ZNnm22tramVUvatWvH+fPnadu2LZA7Z0m9evVo2bKlLNlUAmQ6ynIiLi4Od28fzNqOwdJJ/hQUoHKySVw/k8Fd2/DlF/O1LkcUQKajLOfs7OzY9uMmUvYEk/77Wa3LEWWAzkyPdc83CFm3iQVffKF1OeIhSACXI97e3nyzeBGJm2aTlXhD63JEGWBW1Rpr/3eY8s40duzYoXU5oohkSaK/sba2/kdDPh/WtGnT6NixI127di2w3aBBAzl56hRzl8ym+sAPMDOvek+bK8FjMbOoBjodequa1On1Gnrrmv+//X99ebW6TaSqvTMpZ8K4vm4WDccFY167EQBK5XBr12LSLh4FdOgMFtTpOwXzGvIgflljXrMh1k++wcAhwwgP+xknJyetSxKFJAFcRrz//vuFbjv9vXc5GX2K3dvnYd3zDXS6e/+QeWTYh+gtbbm1bxkJv3xPra7P59n+V8mn9lPFvgXJ0fup0WE4ACmnDpCdeIMGY4PQ6czIuhOPzuLesBdlQ9VH3chuM5zHH38cV1dX9PrcZa6GDh3K1KlTOXDgAC+88ALm5uaEhYUxbdo0tmzZQs+ePbGyssLa2po33nijWGuaO3cu48ePx9Ly4dc9rOikCyIfe/fupVOnTgwcOJDHH3+c4cOHmx4Di4iIoF27dri7u9OqVSsSExNJS0tjzJgxuLm54enpyZ49ewBYunQp/fr1o1u3bjg4OBAUFMSnn36Kp6cnbdq04ebNmwCMHj2aNWtyn3JwcHBg+vTpeHl54ebmxm+//QbA9evX6datG66urlhXq0LimV9I3Le0wOuo2siVzFtX830/JyOV9Csnqf3kyyT/tt+0PTvpJnrrWqZwN1Svg76q9cN9mKJUWBm7g5keK9taREREEBUVxdSpUwFYuXIlb731FlFRUVSrVo1FixZx7NgxPv744xKrZ+7cuaSkpJTY8SsCCeACHDlyhLlz5xIdHc358+f5+eefycjIYMiQIXz++eccPXqUnTt3Uq1aNRYsWIBOp+P48eOEhITwzDPPkJaWBsCJEyf44YcfiIiI4J133sHS0pIjR47Qtm3bfCcxqVOnDpGRkUyYMIE5c+YA8O9//5vOnTtz8uRJhgwZQmZGBoaLv5AcvTffa0g9F45FXQfT6z9C3ubqN5O5tvw1AFJifqFqE2/Ma9mhr2pj+oLP8vEOpJwN5+o3k7m5+6sSXX1ZFB+d3pwzNzN4dvwLphuGr776iu+//5733nuP4cOH4+/vT1JSEt7e3qxatSrP/p06deLuk07x8fE4ODgAkJ2dTUBAAD4+PhiNRhYuXAjkf6Myb948rl69iq+vL76+vqX3AZQz0gVRgFatWmFvbw+Ah4cHFy5cwNbWlgYNGuDj4wNA9erVAfjvf//L5MmTAXj88cdp3LgxZ86cAcDX1xcbGxtsbGywtbWlT58+QO4zmMeOHbvvuZ966ikg94u3H374wXSOdevWAdCjRw9q1qzJmlUh9PLvh8H2EarYOZv2/yPkbTAzw6KuA7U6jDRt/3sXRMqp/dh4+wNg6dyRlOh9VKnviKF6HeyeW0jaxaOkXTrGH6HvUKfvVKo5ePyDT1SUNJWVQdL1a6xcGcnObVv57LNPGTduHP/973/p3bs3AwcOBHK/64iKigJgxowZDzzukiVLsLW1JSIigvT0dNq3b0/37t2B3BuVkydP0rBhQ9q3b8/PP//MSy+9xKeffsqePXtMzxiLe0kAF6BKlSqm/9br9YV+4L2g45iZmZlem5mZ5XvMu20edN4WLVrw3bfLGPbMWPRD/oPBth5w/77ev8tOTSTt4rH/TfijA5UD6KjhOxadTofOYE61x1pS7bGWmFnWIDXmFwngMk5nsKDhswtIu3yCP79/j0GDBhXLcbdv386xY8dM3WQJCQnExMRgYWFx3xuVf/3rX8Vy3opOuiCKyMnJiWvXrhERkTtVZGJiIllZWXTo0IGVK1cCcObMGS5dulTs30a3b9+e77//Hsj9hbh16xYAffr0YdpbU0jcOIuc9ML3uaWc/hkrF1/sJ3yD/YSvsZ+4FEONR0i/cpL038+aHnVTKofM67Hoq9cr1usRJSMnPYXUvYv46KP/FHn0msFgMM05fLcLDUApxfz5801zAcfGxprugIvrRqUykgAuIgsLC1atWsXkyZNxd3enW7dupKWlMXHiRHJycnBzc2PIkCEsXbo0zz/M4jB9+nS2b9+Oq6srq1evpn79+tjY2AAQ8Mbr+HfrSNK2zwp9vORT+7Bs3jbPNsvm7UiO3kdOSgLX177P1SUTufb1i2Cmp7p372K9HlEykrZ+wlM9OvPaqy8XeV8HBwcOHz4MYLrbhdwVNIKDg8nMzJ0U6syZMyQnJxd4LBsbGxITZY3DgshQ5HIkPT0dvV6PwWAgLCyMCRMmmPrxADIyMujg25Vzqh7WHcdoV6jQzMWPemNpZYXjY4+h0+no0aMHgYGBjB49+p4+4LvPu8+YMcP0GNpvv/3G4MGD0ev19OrVixUrVnDhwgVycnJ499132bRpE0op6taty/r16zly5Ahz5sxh8+bNALz44ou0bNmS0aNHM3/+fIKCgmjYsKHpqaDKKr+hyBLA5UhMTAyDBw8mJycHCwsLvvjiC9OXgXfduHEDd28f0lv0yX0sSVQayce2UTX6R45GRlCrVi2tyxF/kV8Ay5dw5UizZs04cuRIgW1q165NyLfL6NS5M5auXdCZ6UupOqGltIvHyPjlO37+5aCEbzkiAVzBpKWl8dJrAdRsM0jCt5LIvBlH0k9z2LD2e5o3b651OaII5Eu4Mmrs2LHUq1evSJNgK6UYOXosl9KqYNVuWAlWJ8qK7LQkEjfO4j+zP6BLly553tu6dStOTk44OjoSGBh4z76XLl3C19cXT09PjEYjW7ZsAWDHjh14e3vj5uaGt7c3u3fvNu0TEhKCm5sbRqORHj16EB8fD+T2I9vZ2ZmWrL97LFEwCeAyavTo0WzdurVI+8z8YBY7w45g1f2l+84PISoWlZ1F8paPGT6gLxMnTMjzXnZ2NpMmTeKnn34iOjqakJAQoqOj87T54IMPGDx4MEeOHCE0NJSJEycCuaMwN23axPHjx1m2bBkjR+YO5MnKyuLll19mz549HDt2DKPRSFBQkOl4r776qukxtZ49e5bw1VcM8ltaRnXs2LFIfXlr167lP3ODsOr91n1nSBMVi1KK5L1fYXy0NkHz7n30MDw8HEdHR5o2bYqFhQVDhw5lw4YNedrodDru3LkD5A6saNiwIQCenp6m/3ZxcSE1NZX09HSUUrnnTU5GKcWdO3dM7cTDkQCuACIjI3nm2eew6fMWBpvaWpcjSkHKkc1UTzjLxh9Wm2Y++6u4uDgaNWpkem1vb09cXFyeNjNmzGDFihXY29vTs2dP5s+/d1WNtWvX4uXlRZUqVTA3Nyc4OBg3NzcaNmxIdHQ0zz77rKltUFAQRqORsWPHmgYJiYJJAJdzV69epXvP3lj6TqBKfUetyxGlIPX8YbIif2DXti2muUgeRkhICKNHj+bKlSts2bKFkSNHmkbBAZw8eZIpU6aYJt7JzMwkODiYI0eOcPXqVYxGI7NnzwZgwoQJnDt3jqioKBo0aMDrr7/+zy6ykpAALsdSUlLo2qMXOHeXdeIqiYz4SyRv/5xN63+gSZMm+bazs7Pj8uXLptdXrlzBzs4uT5slS5YwePBgANq2bUtaWprpS7UrV67Qv39/li9fzmOPPQZgGvTz2P8GeQwePJiDBw8C8Mgjj6DX6zEzM+O5554jPDy82K65IpMALqdycnIY/PQI/jCrg1WrgVqXI0pBdkoCSRtnETT30wdOduPj40NMTAyxsbFkZGQQGhqKv79/njaPPvoou3btAuDUqVOkpaVRt25dbt++Ta9evQgMDKR9+/am9nZ2dkRHR3P9+nUg92kJZ+fcGfiuXbtmardu3TpZwr6QJIDLqGHDhtG2bVtOnz6Nvb09S5YsyfP+W++8x8/HYrDuOgmdTqdRlaK0qKxMkn78iOeeeZrRo595YHuDwUBQUBB+fn44OzszePBgXFxcmDZtGhs3bgTgk08+YfHixbi7uzNs2DCWLl2KTqcjKCiIs2fP8v7775seK/vzzz9p2LAh06dPp2PHjhiNRqKionj77bcBePPNN02Pp+3Zs4fPPiv8nCSVmQxFLodWrFjJhFffpMbQ/6C3qqF1OaKEKaVI2jGfVg2r8uPGdUWe4UxoT4YiVxBhYWG88OJLVB/wvoRvJZEc8QN1M35nzaqDEr4VjARwOXLx4kV6+ffDqtvkPMsMiYorJeYXOPkTOw9HYGVlpXU5ophJAJcTiYmJdPHrid7Dn2qP+Tx4B1HuZfxxnpRdC9izY5tpxQlRscjfM+VAdnY2/QcO4bZ1Yyy9+mpdjigF2Um3SNz0IV99GXzPlKOi4pAALgdefeNNIs//jpXv8/LEQyWQk5lO0uYPeWXS8wwdOkTrckQJkgAu4xYtXszSkNVY9XwTnV56jCo6pRTJO+fTqaUr78+YrnU5ooTJb3QZtm/fPl4LmIrt4Nnoq9loXY4oBcm/hGKnT+K7b7fIXzuVgARwGXX27Fn6PjUQqx6vYV7L7sE7iHIv+dR+zGL2sj3yENWqVdO6HFEKJIDLoNu3b9PFryfmPkOo5uChdTnlllI5XF30fJ5tWbev5dO68Gp2Gf+Pj/FXhup10VvXIm3fYv67bw/169cv1uOLsksCuIzJysqid78BJNdpgbXHk1qXU74pVSyB+3eDmhfvr82i4FlY1ajDymXf4O7uXqzHFmWbBHAZM2HSZKJ/T8LG/2WtSyn3dGZ6Gk/ZXGCb7OTbpF44gjq1A+78wWuvvsKbb7yOubl5KVUJy5cu5e0pb9C3rzxiWNnIXBBlyOfzgpgW+Cm2QwIxqyKjnkpbxp+xpB/8lrqGNL7/7luMRmOpnDcxMRFra2v50q0Cy28uCHkMrYzYtm0b70z/N9b+70j4asSiXhOs+75H/KOdaNuhE18EB5fKeW1sbCR8KynpgigDTp06xaChT2PdawrmNeQLGC3pdDqsjd2p0siNqdNncCcxialvBmhdlqig5A5YY/Hx8XT160mV9s9Q1d5F63LE/5jXbED1gbOY9fFnfPW3uZiFKC4SwBrKyMjgyT79SG/kg5VrF63LEX9jqF6H6v2m8cprb3D27FmtyxEVkASwRpRSjBk3nvN3wOpfI7UuR+TDvHYjqvoMZNjIZ/IsWClEcZAA1kjgfz7mxz0HsfJ7BZ1O/jeUZVZefTj7ewLLli3TuhRRwchjaBrYtGkT/Z8aQI1uE9FXr/vQx0m/dJzU2MPFWFle1sbuGGo2LHT7pCNbMK/TmCqNynZftlkVSyzqOqAzWBR6n9Rzh6gRvZbfThyVJxZEkcmSRGVIcloGXq3bwp0jcOfhjxMRtr/4irqPm7+fxad9x0K3jzhzEM4cNL0uyr6lJScnh9u3bnPtahzVXLpQzasveuuaD9yvalMvft+/mPDwcFq3bl0KlYrKQO6ARbG6c+cOw4cPZ/Pm3BFoOTk5ZfKOMSYmhk8++5yV36+her8ZmNdp9MB9En9dTb9mVVmy6MtSqFBUJDIQQ5SK6tWrs2nTJs6cOQNQaqPJiqpZs2Z8+UUQcz/6kDsb3icnPfmB+5jbufJz2K+lUJ2oLCSARYlo1qwZy5Yt48SJE6SmpmpdTr6efXYsnf7VjuQTux/Y1qJeE87H/EZGRkYpVCYqAwlgUWJGjRoFwPPPP/+AltoKePUlsk9u40HdcWYWVbGqXZ/Tp0+XUmWiopMAFiXK3Nycb7/9VusyCtShQwdSblxDZaY/sK2hmg137vyDb06F+AsJYFGiGjRooHUJD6TT6ahqaYXKTHtwWzODdEGIYiMBLEpUWXwC4v4KWWe5uR5RHkgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIrAknKrXTp0+zcuVK0lMfvCLGXd988w23bt3iqaeeKsHKRGUgd8CiUouJOcusD2djaNoKM0vbB7ZXzt35dsUKftj4YylUJyo6CWBRqfXu3YudO7Zj+P0kKce2Fdg24/oF0n5ezptT3mL514tLqUJRkUkXhKj0fH19OfRrGF26P0nSrTisOo5GZ6bP0yb1XATJO+YTPP9zRo4coVGloqKRO2AhyF1E9GhkBM0tbpG46UNy0lMAUEqRfHg9mXuD2fHTZglfUawkgIX4n5o1a7Jv13b6tjeSsGoKmTfjSNr1BdUv/cyRQ+G0bdtW6xJFBSMBLMRfmJub8/XihUwPeIk/vnkRj1qKqEO/0rhxY61LExWQ9AEL8Tc6nY7XXn2FIYMHUb9+ffR6/YN3EuIhSAALkQ87OzutSxAVnHRBCCGERiSAhRBCIxLAQvxNp06dOHToEAA9e/bk9u3b97SZMWMGc+bMKeXKREUjfcBCFGDLli1alyAqMLkDFpVCcnIyvXr1wt3dHVdXV1atWsWuXbvw9PTEzc2NsWPHkp6efs9+Dg4OxMfHAzBr1iyaN2/Ov/71L06fPm1qM2/ePFq0aIHRaGTo0KGldk2i/JM7YFEpbN26lYYNG/Ljj7mT6CQkJODq6squXbto3rw5o0aNIjg4mFdeeeW++x8+fJjQ0FCioqLIysrCy8sLb29vAAIDA4mNjaVKlSr37a4QIj9yBywqBTc3N3bs2MGUKVM4cOAAFy5coEmTJjRv3hyAZ555hv379+e7/4EDB+jfvz+WlpZUr14df39/03tGo5Hhw4ezYsUKDAa5pxGFJwEsKoXmzZsTGRmJm5sb7777LuvXry+2Y//4449MmjSJyMhIfHx8yMrKKrZji4pNAlhUClevXsXS0pIRI0YQEBBAWFgYFy5c4OzZswB8++23PPHEE/nu37FjR9avX09qaiqJiYls2rQJgJycHC5fvoyvry8fffQRCQkJJCUllco1ifJP/l4SlcLx48cJCAjAzMwMc3NzgoODSUhIYNCgQWRlZeHj48MLL7yQ7/5eXl4MGTIEd3d36tWrh4+PDwDZ2dmMGDGChIQElFK89NJL1KhRo5SuSpR3OqVUoRu3bNlS3X0+UojCcHBw4OLFixTl35kQFY1OpzuslGr59+3SBSGEEBqRABZCCI1IAAshhEYkgEWFtnXrVpycnHB0dCQwMPCe9/fv34+XlxcGg4E1a9aYtu/ZswcPDw/TT9WqVU2PrgUFBeHo6IhOpzONkvuriIiIe44nxP1IAIsKKzs7m0mTJvHTTz8RHR1NSEgI0dHRedo8+uijLF26lKeffjrPdl9fX6KiooiKimL37t1YWlrSvXt3ANq3b8/OnTvvu0pGdnY2U6ZMMbUVoiDyGJqosMLDw3F0dKRp06YADB06lA0bNtCiRQtTGwcHBwDMzPK/F1mzZg1PPvkklpaWAHh6eubbdv78+QwYMICIiIhiuAJR0ckdsKiw4uLiaNSokem1vb09cXFxRT5OaGgow4YNK9T51q1bx4QJE4p8DlE5SQCLEqGUwrdbDy5evAiAu3er+842VtZdu3aN48eP4+fn98C2r7zyCh999FGBd9NC/JV0QYgS4+XhzoH9e9GZGXjcqTk6na5Uz29nZ8fly5dNr69cuVLkdd6+//57+vfvj7m5+QPbHjp0yDQdZXx8PFu2bMFgMNCvX78inVNUHhLAokTodDo++fgjWnp5En/zJi9OnFDqAezj40NMTAyxsbHY2dkRGhrKd999V6RjhISEMHv27EK1jY2NNf336NGj6d27t4SvKJD8rSRK1LBhQ5k8aWKphy+AwWAgKCgIPz8/nJ2dGTx4MC4uLkybNo2NGzcCuY+M2dvbs3r1ap5//nlcXFxM+1+4cIHLly/fM0nPvHnzsLe358qVKxiNRsaNG1eq1yUqDpkLQgghSpjMBSGEEGWMBLAQQmhEAriS8/X1Zdu2bXm2zZ07N99nWf+6SGV+Pvzwwzyv27VrB+T2qbq6ugK5Twy89NJLD1u2yejRo2XIryi3JIAruWHDhhEaGppnW2EHHuTn7wF88ODBe9q0bNmSefPmPfQ5hKgIJIAruYEDB/Ljjz+SkZEB5N6lXr16lbi4ONzc3HB1dWXKlCn33bdfv354e3vj4uLCokWLAJg6dSqpqal4eHgwfPhwAKytre/Zd+/evfTu3RuAnj17mia9sbW1ZdmyZWRnZxMQEICPjw9Go5GFCxcCuQM8XnzxRZycnOjatSt//vlnsX8mQpQapVShf7y9vZWoeHr16qXWr1+vlFJq9uzZasyYMapRo0bqzz//VJmZmcrX11etW7dOKaVU48aN1fXr15VSSt24cUMppVRKSopycXFR8fHxSimlrKys8hz/7uvY2Fjl4uKilFJqz549qlevXnnaHTp0SLm5uanbt2+rhQsXqpkzZyqllEpLS1Pe3t7q/Pnzau3atapr164qKytLxcXFKVtbW7V69eoS+FSEKD7AIXWfTJU7YJGnGyI0NJTGjRvTqVMn6tati8FgYPjw4fddsn3evHm4u7vTpk0bLl++TExMzEPXEB8fz8iRI/nuu++wtbVl+/btLF++HA8PD1q3bs2NGzeIiYlh//79DBs2DL1eT8OGDencufNDn1MIrUkAC/r27cuuXbuIjIwkJSUFDw+PB+6zd+9edu7cSVhYGEePHsXT05O0tLSHOn92djZDhw5l2rRppi/plFLMnz/fNCVkbGysTPEoKhwJYIG1tTW+vr6MHTuWYcOG0apVK/bt20d8fDzZ2dmEhITcMxosISGBmjVrYmlpyW+//cYvv/xies/c3JzMzMxCn3/q1KkYjUbTPAoAfn5+BAcHm45z5swZkpOT6dixI6tWrSI7O5tr166xZ8+ef3j1QmhH5oIQQG43RP/+/QkNDaVBgwYEBgbi6+uLUopevXrRt2/fPO179OjBl19+ibOzM05OTrRp08b03vjx4zEajXh5ebFy5coHnnvOnDm4uLiY7rzff/99xo0bx4ULF/Dy8kIpRd26dVm/fj39+/dn9+7dtGjRgkcffZS2bdsW6+cgRGmSochCCFHCZCiyEEKUMRLAQgihEQlgUaIetCqxEJWZBLAoMYVZlViIykwCWJSYv65KbGFhYVqVWAiRSwJYlJjiWpVYiIpKAlgIITQiASxKTHGsSixERSYBLErMX1clzsjIIDQ0FH9/f63LEqLMkKHIosT8dVXi7Oxsxo4dm2fVYSEqOxmKLIQQJUyGIgshRBkjASyEEBqRABZCCI1IAAshhEYkgIUQQiNFegpCp9NdBy6WXDlCCFEhNVZK1f37xiIFsBBCiOIjXRBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoZH/AxXezKO8aR1WAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "WWT.plot_sankey()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABiUElEQVR4nO3dd3gUVdvH8e+kkkYAAZUiRaQFkhASqtIEgpSI0gUBEUVB8LEg+mBBBcTGgxQBEUUUKTY6SJEekBJCl95BaQJJCKnz/rFh3wRICJBkssnvc125dGfPzLlnE87ec86ZM4ZpmoiIiIiIOAInqwMQEREREcksJa8iIiIi4jCUvIqIiIiIw1DyKiIiIiIOQ8mriIiIiDgMJa8iIiIi4jCUvGYDwzBMwzAqWB1HZhiG8YBhGNGGYThbHYuIiBXUZos4lnyVvKb8g7/2k2wYRmyq113T2aeRYRgnsjCGlYZh9M6q490t0zSPmabpbZpmUlYf2zCMKYZhxKd8vhcMw1hqGEbllPeGGIaRcN3v5I1U+7Y2DGOjYRgxhmGcNwxjmmEYpW5Sx5CUL57adxHbtZ+bfhkYhjHhunJxhmFEpXq/iGEYv6XEetQwjKdSvdfTMIyk6/aPNgyjxG3Gm1EdjQ3D2GEYxsWUz+o3wzBK3s7xRXIjtdk3ctQ22zAMN8MwPjcM40TKvkcMwxh1G7HVSzl+lGEY2w3DeDiDsoZhGB+nxHE+5f+NVO8HGoaxxTCMKyn/DUznM7j2sy3zn2Km6mhsGMYKwzAuGYZx5HaPLfkseU35B+9tmqY3cAxok2rbNKvjy6M+Sfm8SwFngCmp3puZ+ndimuYnAIZhtAd+BEYBRQE/IA5YaxhG4Ws7pzRG3YELKf+9o9hS/dz0y8A0zReu+9uZDvyUqsg4IB64F+gKjDcMwy/V++uvq8fbNM1T19eT8uUwJJ1YM6pjNxBqmmYhoASwHxifuY9AJPdSm22J7Gqz3wKCgVqAD9AIiMhMQIZhFAHmAZ8ChYBPgHmpvw+u8zzQFggA/IE2QJ+UY7kBc4AfgMLAd8CclO1pPoNUPwHpxLXSMIxGN9l+qzpigG+AgZk4fbmJfJW8pscwDHfDMEYZhnEq5WdUyjYvYBFQInWPmWEYtQzDWJ/S03XaMIyx1/3hZ7beRilXoa8ZhnEm5VjPpHrfI+VK9WjKFdpawzA8Ut4LMwxjV0oMKw3DqJJqvyOGYQxMuTqNMQxjsmEY9xqGsSjlqnXZtX/0hmGUNWw9ly4pr1cahvGhYRjrUsouMQyjaKpjd0+J57xhGO+k1NX0VudqmuYVbI1btVt8JgbwOTDUNM0fTdOMNU3zb6A3EA28kqr4I8D9wACg8538Dm5Xyt9EO2yNUerX75imGW2a5lpgLvB0NtR50zpM0/znumQ4CUh3CNS4bog0padhaKrXb6T8LZ4yDKP39eVFrKY22yHb7BDgN9M0T5k2R0zTnJqpDx7qAX+bpvmTaZpJpmn+AJwFnkynfA/gc9M0T5imeTIlvp4p7zUCXIBRpmnGmaY5GjCAJpmMJTMyrMM0zY2maX4PHLrVgYybjCSk/h2m/M19ZxjGv4Zh7Elpv7Ns5CG3UvJqMxioAwRiu1KrBbxtmmYM8Bhw6roesyRs/yCLAnWBR4G+d1j3fYAvUBJ4Fhhn/P/V5GdATWz/cIsAbwDJhmFUxNb79x+gGLAQ21Vo6sa4HdAMqIjtqnMR8N+U8k7YEr70PAU8AxQH3IDXAQzDqAp8ia3n7/5Ucd+SYRjeKfttvUXRSsADpO3ZxDTNZOCXlHO6pge2q/FZKa/bpKrvKcMwtt+irr6GbWhsi2EY7W59FoDtcz0LrE55XRFINE1zX6oy27D1PGSVW9Zh2ObBXQRisf2+PrmTigzDaAG8CjTFlgA3urOQRbKV2uy0HKHN3gC8ahhGX8Mwqqckvanr+9IwjC8zCukmr9NLrP2wtZHXpG4v/YDtpmmaqd7fTta22TlRxzXvAWWB8tg+627ZUEeuo+TVpivwgWmaZ0zTPAu8TwY9Z6ZpbjFNc4NpmommaR4BJgIN77DuhJS6E0zTXIjtSrWSYRhOQC/gZdM0T6ZcbYabphkHdAIWmKa51DTNBGwNpge2BvOaMSk9cieBNcCfpmluNU3zKvAbUCODmL41TXOfaZqx2BLDwJTt7YF5pmmuNU0zHngXMNM5xjWvpyRVBwBv/v/qF6BjSi/EtZ8S2L5cAE7f5Finr71vGIYn0AH4MeUz+JlUUwdSegD8M4hrNPAQtsb+HWCKYRj1b3EuYEuYp6ZqlLyBy9eVuYRtWOyaOted58FM1JPaLetImQdXCNvn8zbw123WcU1HbL//XSk9L0Pu8Dgi2Ultdlq5vs0GPgI+xva72wycNAyjx7WCpmn2NU0zvQuK9dh607sYhuGast+DgGc65b2xtZHXXAK8UxLm69+79n7qNvv1687zu3TqSU9m6sgqHYHhpmn+a5rmCWzfbXmeklebEsDRVK+Ppmy7KcMwKhqGMd8wjL8Nw7gMDOf//4HervOmaSamen0F2x9+UaAAcLNEJ028KVe4x0l7Rf1Pqv+Pvclr7wxi+vsm8Vyr93iqeq8A5zM4DsBnpmkWMk3zPtM0w0zTTH0+s1Leu/ZzCjiX8t79NznW/anefwJIxNaDATANeMwwjGK3iOda7BGmaZ5P+TJbmLJ/ekNQgK13E1tPZOqhrmig4HVFCwJRqV5vuO48H0x1zPnXGkjgTeDNVA3m/Nuo49p5XeD/51e5ZHQ+6UjzO77u/0VyC7XZaeX6NjslmR9nmmZ9bPNWhwHfpJ4+kR7TNM8Dj2MbFfoHaAEsA9IbHr++zSwIRKd0OmSmPf3suvO0J9mpk1rgYWB+qm1vplP/zerIKvmyzVbyanMKKJPq9QMp2+DmV6njsfVsPWSaZkFsQzvXD2ncrXPAVWxXl9dLE2/K1WRp4GQWx3C909gm8V+r1wO4J4vr2IutQeqQemNKr0Y7YHnKph7YGuhjhmH8jW3IyhXb8NmdMLn17/BpYJ1pmqnnKe0DXAzDeCjVtgBgV6YqNc3W1xpIYAQwIlWD2foO63DB1qN8feOZmkeq//dN9f9pfsfY/q5Echu12ZmTm9psO9M2L3Yc8C9QNTOVmKa5yjTNENM0i2BriysDG9MpvgtbG3lN6vZyF+B/3bQFfzLfZhdK1WavBVqn2jYiK+q4CY/rXuf7NlvJq8104G3DMIqlTHR/F9tdgmC7yrvHMIzUfyw+2IZxow3bMiIvZnVAKVfm3wAjDdsNB86GYdQ1DMMd27BQK8MwHjUMwxV4DdudneFZHcd1fgbaGLYlS9ywDSln6RdAypXx69h+H08ZhlHAMIz7gK+xJWP/M2zLQD0KtMY2PBaIrXH6mEyuOmAYRnvDMLwNw3AyDKM5tnlCc2+xW3fS3nmLaZtj9yvwgWEYXilTDx4Hvs9MHJlxqzoMw3jSMIxKKedSDBgJbE3phU1Pz5S/qUBs81t9Uv6WZgHPGIZRJWVqxjtZdR4iWUhtdubkijYbwDCM/xi2m488DMNwSRn69+HWc2pJ2b9GypSBgtimXRw3TfP3dIpPxTa/tmTK1IbX+P+2eyW2OdADDNtNfi+lbP/jNk89IxnWkdJWF8DW4WKkfGYZ3UBYxDCMtin79cLWc+1r2JZ3nAW8ZRhG4ZTvxpcyOE6eoeTVZii2OTjbgR3Ylu8YCmCa5l/YGspDqeb4vI6thy8KmATMzKa4Xk+JZxO25aA+BpxM09yLLdkag+1qvw22JWTisykOAEzT3AX0B2Zgu9qLxraUSlwW1zMT25X1K9iGuHZju/KsnzJ89DQQaZrmEtM0/772g22uj79hGNUMw+hqGEZGV7kvY+v1uIht+ZXnTNNcCWkWAX/gWmHDMOpiu7r96cZD0TclvjPY/lZeTPmsrqlr3LjOa8htfiwZ1VESWIzt73EHkIxtWkVGPLH9Didh++LvDjQxTXMRts9xBbY5bxtSymfp71jkLqnNzoRc1GaDbTrD59imOJwD+gHtro1kGbb1tCdkUM0bKfsdxzYdwd7GGYbxiGEY0anKTsR2M+8OYCewIGUbKZ95W2xt3kVs85TbXve7eOO69voctyETdTTANhVkIbZRg1hgSQaH/Bvb3NazQCi2NnoCth7YD7D1fB/GNpXiZ/JBe22Y5s1GWERuzbDdjXoR21DcYYvDkUwyDMPE9js7kImyVbA1/u7XzfMTEQejNtvxGLZ1ZH8wTfOGh/SkU/5FoLNpmnd6Q6JDUM+r3BbDMNoYhuFp2NZT/Azble0Ra6OSrGQYxhMpQ12FsfUczVPiKuKY1GbnbYZh3G8YRv2UKQWVsE2R+M3quLKbkle5XY9ju/ngFLalpjqb6r7Pa/pgG1o8iG3eVpbPDxSRHKM2O29zwzYlIgrbnNo52Nb2zdM0bUBEREREHIZ6XkVERETEYdxW8moYxuLsCkREJDup/RIRyRtu6wk8BQsWDA0ODtY8AxFxRNc/YjfPK1q0qFm2bFmrwxARuW1btmw5Z5rmTZ+aeVvJ60MPPcTmzZuzJioRkRxkGMZ+q2PIaWXLllWbLSIOyTCMo+m9pzmvIiIiIuIwlLyKiIiIiMNQ8ioiIiIiDkPJq4iIiIg4DCWvIiIiIuIwlLxKtnB2diYwMBA/Pz8CAgL4/PPPSU5Ozvb6rv2MGDEiw/Le3t7ZFgvAqFGjuHLlSrbWISIikh/d1lJZIpnl4eFBZGQkAGfOnOGpp57i8uXLvP/++2nKJSYm4uJy53+G1/ZPXV9uMGrUKLp164anp6fVoYiIiOQp6nmVbFe8eHG++uorxo4di2maTJkyhbCwMJo0acKjjz5K9+7dmT17tr18165dmTNnDklJSQwcOJCQkBD8/f2ZOHEiACtXruSRRx4hLCyMqlWrplvvpUuXqFSpEnv37gWgS5cuTJo0yf7+4MGDCQgIoE6dOvzzzz8AzJs3j9q1a1OjRg2aNm1q3z5kyBA+++wz+77VqlXjyJEjxMTE0KpVKwICAqhWrRozZ85k9OjRnDp1isaNG9O4ceMs+xxFREREyavkkPLly5OUlMSZM2cAiIiI4Oeff2bVqlU8++yzTJkyBbAlnOHh4bRq1YrJkyfj6+vLpk2b2LRpE5MmTeLw4cP2/b/44gv27dsHQGxsbJppAzNnzsTX15exY8fSs2dPZsyYwb///stzzz0HQExMDHXq1GHbtm00aNDAntQ+/PDDbNiwga1bt9K5c2c++eSTDM9r8eLFlChRgm3btrFz505atGjBgAEDKFGiBCtWrGDFihXZ8XGKiNwVTe3S1C5HpmkDYolmzZpRpEgRABo2bEjfvn05e/Ysv/zyC+3atcPFxYUlS5awfft2fv75Z8CW2O7fvx83Nzdq1apFuXLl7MdLb9pAs2bN+Omnn+jXrx/btm2zb3dzc6N169YA1KxZk6VLlwJw4sQJOnXqxOnTp4mPj09Tx81Ur16d1157jUGDBtG6dWseeeSRu/pcRERygqZ2aWqXI1PPq+SIQ4cO4ezsTPHixQHw8vJK83737t354Ycf+Pbbb+nVqxcApmkyZswYIiMjiYyM5PDhwzRv3vym+6cnOTmZPXv24Onpyb///mvf7urqimEYgK1HIDExEYD+/fvz0ksvsWPHDiZOnMjVq1cBcHFxSdMrcW17xYoViYiIoHr16rz99tt88MEHt/3ZiIhYSVO7xNEoeZVsd/bsWV544QVeeukle8J4vZ49ezJq1CgAe2MXGhrK+PHjSUhIAGDfvn3ExMTcVt3/+9//qFKlCj/++CPPPPOM/VjpuXTpEiVLlgTgu+++s28vW7YsERERgG3KwrXpC6dOncLT05Nu3boxcOBAexkfHx+ioqJuK1YREatoapc4Ek0bkGxxraFKSEjAxcWFp59+mldffTXd8vfeey9VqlShbdu29m29e/fmyJEjBAUFYZomxYoVS3P1f7P6rmnRogXPPPMMX3/9NRs3bsTHx4cGDRowdOjQG4bFUhsyZAgdOnSgcOHCNGnSxN4Qt2vXjqlTp+Ln50ft2rWpWLEiADt27GDgwIE4OTnh6urK+PHjAXj++edp0aKFvYEUEXEkmtoluZppmpn+qVmzpimSHWJiYszy5cubFy9etDoUyaOAzeZttHd54UdttqTHy8srzeuDBw+aRYoUMZOTk81vv/3W7NevX5r3R4wYYY4cOdKsVauWuWvXLtM0TfPJJ580Fy9efMOxV6xYYbZq1SrD+q5JSkoyH374YbNUqVLm9u3bb1r+p59+Mnv06GGapmk2bNjQnDNnjr2ehg0bmqZpmh9++KH58ccf2/d58MEHzcOHD5umaZrnz583v//+e7NBgwbm+++/b5qmaZYpU8Y8e/bsTWOS3CGjNlvTBsRyy5Yto0qVKvTv3x9fX1+rwxERyVc0tUscjaYNiOWaNm3K0aNHrQ5DRCTf0NQuTe1yZIatZzZzgoODzc2bN2djOCIi2cMwjC2maQZbHUdOUpstWeXKlStUr16diIgIjZBJjsiozda0AcmTFi9eTKVKlahQocItF8MWEZH0aWqX5DaaNiB5TlJSEv369WPp0qWUKlWKkJCQW643KCIiN6epXZLbqOdV8pyNGzdSoUIFypcvj5ubG507d2bOnDlWhyUiIiJZQMmr5DknT56kdOnS9telSpXi5MmTFkYkIiKpaWqX3A0lryIiIpJjrk3tWrRoEbt372b69Ons3r3b6rDEgSh5lTynZMmSHD9+3P76xIkT9nUBRUTEWpraJXdLyavkOSEhIezfv5/Dhw8THx/PjBkzCAsLszosERFBU7vk7mm1AclzXFxcGDt2LKGhoSQlJdGrVy/8/PysDktERESygJJXyZNatmxJy5YtrQ5DRESuo6ldcrc0bUBERERyjKZ2yd1Sz6uIiIjkGE3tkrul5FVERERylKZ2yd3QtAERERERcRhKXkVERETEYSh5FRERERGHoeRVRERERByGktd8ZNiwYfj5+eHv709gYCB//vmn1SFlmyNHjrBt2zarwxAREZEsptUG8on169czf/58IiIicHd359y5c8THx99yv8TERFxcHOvPZPny5TzZoROJScl89vFwXnzhBatDEhERkSyintd84vTp0xQtWhR3d3cAihYtSokSJdi0aRP16tUjICCAWrVqERUVxZQpUwgLC6NJkyY8+uijxMTE0KtXL2rVqkWNGjWYM2cOAElJSQwcOJCQkBD8/f2ZOHEiACtXrqRRo0a0b9+eypUr07VrV0zTzJHzXLJ0Gc1DWxBruhJvuNH3xRf55POROVK3iMjdyk8jZCJ3SslrPtG8eXOOHz9OxYoV6du3L6tWrSI+Pp5OnTrxxRdfsG3bNpYtW4aHhwcAERER/Pzzz6xatYphw4bRpEkTNm7cyIoVKxg4cCAxMTFMnjwZX19fNm3axKZNm5g0aRKHDx8GYOvWrYwaNYrdu3dz6NAh1q1blyPnWaVyJYYPG0rCxb9JvPQPI0aM4NFGDXOkbhGRu5F6hGz79u0sW7aM0qVL33K/xMTEHIgu6yQkJNCjV2/uK/UAmzdvtjoccUBKXvMJb29vtmzZwldffUWxYsXo1KkTEydO5P777yckJASAggUL2qcINGvWjCJFigCwZMkSRowYQWBgII0aNeLq1ascO3aMJUuWMHXqVAIDA6lduzbnz59n//79ANSqVYtSpUrh5OREYGAgR44cyZHzLF26NIMGDaJGjRoADBo0iJo1a+ZI3SIidyM/jJBdvHiRRxo3ZdbcRUQXfoh6Dz/CTz/9nO31St7iWJMZ5a44OzvTqFEjGjVqRPXq1Rk3bly6Zb28vOz/b5omv/zyC5UqVUpTxjRNxowZQ2hoaJrtK1eutDe+1+p1tJ4BEZGc1rx5cz744AMqVqxI06ZN6dSpE3Xr1qVTp07MnDmTkJAQLl++nGaEbPv27RQpUoT//ve/NGnShG+++YaLFy9Sq1YtmjZtyrRp0+wjZHFxcdSvX5/mzZsDthGyXbt2UaJECerXr8+6det4+OGHs/UcExMTiY2N5er5U3D+FH41gvnrwMFsrVPyHvW85hN79+6194oCREZGUqVKFU6fPs2mTZsAiIqKummSGRoaypgxY+xX5Vu3brVvHz9+PAkJCQDs27ePmJiY7D4VEZE8KT+MkBUtWpRtWzYyePBgAHZGbOKdtwZle72St6jnNZ+Ijo6mf//+XLx4ERcXFypUqMBXX33FM888Q//+/YmNjcXDw4Nly5bdsO8777zDf/7zH/z9/UlOTqZcuXLMnz+f3r17c+TIEYKCgjBNk2LFijF79uycPzkRkTxCI2Qit6bkNZ+oWbMm4eHhN2wvWrQoGzZsSLOtZ8+e9OzZ0/7aw8PDPk8qNScnJ4YPH87w4cPTbL/W8F4zduzYuwteRCQf2Lt3L05OTjz00EPA/4+QLV68mE2bNhESEkJUVJR92kBq10bIxowZg2EYbN26lRo1athHyJo0aYKrqyv79u2jZMmSOX1qIllKyauIiEguoBEykcwxbufuwuDgYFPLWogjCAoKYuvWrTm2vqzkfoZhbDFNM9jqOHKS2mzJrd5++22GDRumNlrSlVGbrRu2RERERMRhKHnNB5ydnQkMDKRatWq0adOGixcvZlh+woQJTJ069bbruXjxIl9++eUdRikiIiJya0pe8wEPDw8iIyPZuXMnRYoUyfDuVYAXXniB7t2733Y9Sl5FREQkuyl5zWfq1q3LyZMnATh48CAtWrSgZs2aPPLII/z1118ADBkyhM8++yzDMv/88w9PPPEEAQEBBAQEEB4ezptvvsnBgwcJDAxk4MCBnD59mgYNGth7fdesWWPNSYuI5HIaIRPJPCWv+UhSUhLLly8nLCwMgOeff54xY8awZcsWPvvsM/r27XvDPumVGTBgAA0bNmTbtm1ERETg5+fHiBEjePDBB4mMjOTTTz/lxx9/JDQ0lMjISLZt20ZgYGBOnq6IiMPQCJlI5mmprHwgNjaWwMBATp48SZUqVWjWrBnR0dGEh4fToUMHe7m4uLg0+2VU5o8//rBf9Ts7O+Pr68u///6bZv+QkBB69epFQkICbdu2VfIqIpIJdevWZfv27YBt9Ktfv36cPXsWT09PJk2aROXKlRkyZAje3t68/vrr6Zb5559/eOGFFzh06BAA48ePZ/To0fYRsmbNmvHqq6/SqVMnLl++TGJiIuPHj+eRRx6x8vRFbknJaz5w7Yr+ypUrhIaGMm7cOHr27EmhQoWIjIxMd7/k5ORblslIgwYNWL16NQsWLKBnz568+uqrd9RTICKSX1wbIXv22WcB2+jXhAkTeOihh/jzzz/p27cvf/zxR5p90itzbYTst99+IykpiejoaEaMGMHOnTvt7frnn39OaGgogwcPJikpiStXruT0KYvcNk0byEc8PT0ZPXo0n3/+OZ6enpQrV46ffvoJsD1CcNu2bWnKFyxYMN0yjz76KOPHjwdsje2lS5fw8fEhKirKvv/Ro0e59957ee655+jduzcRERE5cZoiIg7n2gjZfffdxz///HPDCFlgYCB9+vTh9OnTafbLqMwff/zBiy++CPz/CNn1QkJC+PbbbxkyZAg7duzAx8cn+09W5C4pec1natSogb+/P9OnT2fatGlMnjyZgIAA/Pz8mDNnjr2cYRgA6Zb54osvWLFiBdWrV6dmzZrs3r2be+65h/r161OtWjUGDhzIypUrCQgIoEaNGsycOZOXX37ZknMWEcntro2QHT16FNM0GTduXJrRr2s/e/bsSbNfZspk5NoIWcmSJenZs+cd3QQmktM0bSAfiI6OTvN63rx59v9fvHjxDeXPnz9PmTJlAChXrtxNy9x7771pkt1rfvzxxzSve/TocUcxi4jkR9dGyNq2bUvfvn3to18dOnTANE22b99OQECAvXzqEbLry1wbIfvPf/5jnzZwsxGyUqVK8dxzzxEXF0dERISmd0mup55XSeOdd97hzz//tK9IICIiOUsjZCIZU8+rpPHhhx/y4YcfWh2GiEi+ohEykcxTz2s+sXjxYipVqkSFChUYMWLEDe/HxcXRqVMnKlSoQO3atTly5Ij9vY8++ogKFSpQqVIlfv/991sec+zYsVSoUAHDMDh37ly2npeISH6jETLJ75S85gNJSUn069ePRYsWsXv3bqZPn87u3bvTlJk8eTKFCxfmwIEDvPLKKwwaNAiA3bt3M2PGDHbt2sXixYvp27cvSUlJGR6zfv36LFu2zN4rICIiWefDDz9k48aN3HPPPVaHImIJJa/5wMaNG6lQoQLly5fHzc2Nzp073zCUNGfOHPvQUfv27Vm+fDmmaTJnzhw6d+6Mu7s75cqVo0KFCmzcuDHDY9aoUYOyZcvm9GmKiDg8jZKJ3JqS13zg5MmTlC5d2v66VKlSnDx5Mt0yLi4u+Pr6cv78+XT3zcwxRUQk8zRKJpI5Sl5FRERyAY2SiWSOktd8oGTJkhw/ftz++sSJE5QsWTLdMomJiVy6dIl77rkn3X0zc0wREck8jZKJZI6S13wgJCSE/fv3c/jwYeLj45kxY8YNd6mGhYXx3XffAfDzzz/TpEkTDMMgLCyMGTNmEBcXx+HDh9m/fz+1atXK1DFFREREspqS13zAxcWFsWPHEhoaSpUqVejYsSN+fn68++67zJ07F4Bnn32W8+fPU6FCBUaOHGmf1O/n50fHjh2pWrUqLVq0YNy4cTg7O6d7TIDRo0dTqlQpTpw4gb+/P71797bs3EVEHIVGyUQyyTTNTP/UrFnTFHEENWrUMG1/3iI2wGbzNtq7vPCjNtuxJCQkmOXKlTMPHTpkxsXFmf7+/ubOnTvTlBk7dqzZp08f0zRNc/r06WaHDh1M0zTNnTt3mv7+/ubVq1fNQ4cOmeXKlTMTExMzdcwyZcqYZ8+ezZmTTDF48GC10ZKhjNps9byKiIjkAholE8kcw5bcZk5wcLC5efPmbAxHJGsEBQWxdetWbufvW/I2wzC2mKYZbHUcOUlttuRWb7/9NsOGDVMbLenKqM1Wz6vc0qxZPzF0+AiSkpKsDkVERETyORerA5DcKzk5mcHvvMfYrybjUrAYq1av4defZuDj42N1aCIikoppmiQkJODm5mZ1KCLZTj2vkq7Zs+cwYvhQos+d5uKh7Sz7fSFvDn7H6rBERCSV+Ph4uj/zLD4FfZk16yerwxHJdup5lXS1bPkYkZGRPP744xw9epTIyEgtsSIikotER0fzaGhLtu7YjVuFunTu0oX9Bw4y+L9vWh2aSLZRz6ukq0CBAgQEBPDAAw8AEBAQQNGiRS2OSkRErnFzc6NIkSIU9nInZtcKGj7ajPikZKvDEslW6nkVERFxUG5ubiyaN5uff/6ZDh06sGLJIqtDEsl26nkVEREREYeh5FVEREREHIaSVxERERFxGEpeRURERMRhKHkVEREREYeh5FVEREREHIaSVxERERFxGEpeRURERMRhKHkVEREREYeh5FVEREREHIaSVxERERFxGEpeRURERMRhKHkVEREREYeh5FVEREREHIaSVxERERFxGEpe5Zbi4uKsDuG2OWLMIiIicmtKXuWWDhw4YHUIt80RYxYRuVMJCQlWh3BbHC1eyV2UvMotXbhwge7du1sdxm2Jj4+nb9++VochIpIjDh06ZHUIt+Xw4cNWhyAOTMmrZOi7774D4KuvvrI4ksz78ssvARg1apS1gYiI5JC3337b6hBuy08//UTx4sWtDkMclJJXSVdiYiI9e/akQIECuLu7Wx1OpiQkJNCvXz/uueceXF1drQ5HRCTbXbx4EYDff//d2kAyae/evQCsWbPG4kjEUSl5lZuKj4+3J39///23xdFkTlxcHG5ubgAcOXLE2mAy4cKFC7w1+B1iY2OtDkVEHFjJkiUBaN68ucWRZE7lypUBqFixosWRiKNS8ippxMfH8/zzz9t7Wo8ePYqvr6/FUWUsPj6e3r17U6BAAQBOnDiBt7e3xVHd2rFjxxgxfChdnu6BaZpWhyMiDigoKIgrV66wfv16q0O5JdM08fLyAmD//v0WR5P15sybz5hxX1odRr7gYnUA2Wnh4t/59bfZVodxg8lfTbjjfR9u0JCKlSpnYTQ230yaeMO2p57uwdDhI7K8rqwyaeL4G7Z16/EM7384zIJobt+5c+coeF8ZVm3exTvvDWHoB+9bHZKIOIDk5GTefvttPvroIwAWLVpEnTp1LI4qfQkJCbz88suMH29rszdu3EiFChUsjirrzZ03n28mTaTEfffSrl07q8PJ0/J08hoRsZVvv/2WgvW7YrjkjfmP2+OKsvuIkeXHLVC2BlePbMW1eHm8/ZsBsOg0QHKW15VV3B+oTtyxHbjd+yBe1ZsCsOAk5OaY0yqCx6Mv4VroPr74ciB+VarQpUtnq4MSscTqNWt5b+hwq8NIY+WSRXe8b/XAmtxTrFgWRgPJySarl6ed1+pfI5hPR43h01FjsrSurPD3qZPs3rEtzbbgOvV4850h1gSUzfbs2ol3QCg9nn2OcuXKERQUZHVIeVaeTl4HvzWITVu2sPbAYbxbvIphZH3SdycK1mxjdQg38AlqbXUIty03fo53yifsvzz3Yl8efLA8tWrVsjockRxnJiex+o+lFHz4aVyLlLQ6HAA8K0dz9ei2Wxe8idOlm/C3a4EsjghcCm/HrVg5vPwaAQangdNZXkvWSC5TBZfjZ3AvWRnPivUAOGpxTNkqpAZFygdz5UANQlu2ZlvEZkqUKGF1VHmScTtz7YKDg83NmzdnYzhZLzY2llr1HuG0b1W863SyOhyRdF3Z/ydJaycRuXkjpUuXtjqcPMcwjC2maQZbHUdOcrQ2e/yECQwa8hG+nT/BuUDun7cukp7oDbMofmE7mzesw9PT0+pwHFJGbXaev2HLw8ODJQvn4bTvD2L+Wmt1OCLp8nyoNk7VWtK0RSuio6OtDkckx734wgt069CW6AUfYyYlWh2OyB3zqt2Bf4x76PhUN5KTHWUqm+PI88krwP3338+ShfO5uuor4k7nvTscJe/wDH6C8+4laN/pKTV4ki+NGTWSwDLFiFn5lVbhEIdlGAbezfqxbts+/vv2u1aHk+fki+QVIDAwkKnffE3U/I9IvHzO6nBEbsowDLwefYE//zrK628MsjockRzn7OzMnF9/ouClQ8REzLc6HJE7Zri44d36TcZN+oYffphmdTh5Sr5JXgGeeOIJ3nztP0TPG05y/FWrwxG5KcPZFe9Wg/h66nS+/XaK1eGI5LiCBQuy/PeFJG39ldiDjjNnV+R6zl6F8QkbzAv9+jvEWryOIl8lrwCD33qTFg1qEbNkFKapYVnJnZw9ffEOG0z/V15j7VrN1Zb8p1y5csyf8xsxS0cTfzZP36MueZxbsbJ4NutPq7C2HD2qv+WskO+SV8Mw+O6brynnnUzMOnXjS+7lVvQBvJq/TJu2T3Lo0CGrwxHJcfXr12fcF/8jat4wkq5csjockTvmWaEWTgFhNG3RiqioKKvDcXj5LnkFcHd3Z/H8Obgd20DMrj+sDkckXR7la+IS1I5HQ1ty6ZK+vCX/6dGjO316Pk30/BGYiQlWhyNyx7xqPs4Fz9I82aEzSUlJVofj0PJl8gpQtGhRli1eSNzaKVw9sdvqcETS5VmjFVGFHuLxJzuQmKjlgyT/+WTEcOr4lSPmjy+1AoE4LMMw8G7Sh80HTvLK629YHY5Dy7fJK0DVqlWZ+eMPRC/4hMRL/1gdjshNGYaBV6PebD9+gf4vv2p1OCI5zsnJiZ9n/EjR+H+I2fSr1eGI3LFrN+RO+XEWk77+2upwHFa+Tl4BHnvsMT54922i5g4lOe6K1eGI3JTh7IJ3y4H8+Oscvhw/3upwRHKcl5cXyxYvgF2LuLJPd22L43L2KIhP2GBeeX0Qq1atsjoch5Tvk1eAV/4zgCcfa0r04s8xkzUPRXInpwLeeIe9zRv/fYdly5ZZHY5IjitVqhSL5s/lyh9fEv/PQavDEbljrveUxqvFqzz+ZHsOHtTf8u1S8optWPar8eOoUsyDmDVTrA5HJF2uhUvg/dhrtOvYmb1792bZcV955RVGjRplfx0aGkrv3r3tr1977TUMw2D27Nn2bZUqVWLo0KH21+3atWPmzJkEBgYSGBiIt7c3lSpVIjAwkO7du2dZrJK/hYSEMHniBKLmDScx+oLV4YjcMY+ygbiGdOTR5o9x8eLFbKunUaNGbN5sWy+5ZcuWN61ryJAhfPbZZ9kWQ1ZT8prC1dWV+XN+xeufHcRs/93qcETSVeABf9zqdKVpaEsuXMiaL+/69esTHh4OQHJyMufOnWPXrl3298PDw/nwww/tZc6fP4+Xl1eaRbfXr19Pw4YNiYyMJDIykuDgYKZNm0ZkZCRTp07NkjhFADp16sgrL71IzPyPSE6IszockTvmFdiS6Huq0OaJ9jlyQ+7ChQspVKhQtteT3ZS8plK4cGGW/b6Q+A0/cvXodqvDEUmXl39zYu8PpFXYEyQk3P3yQfXq1bMnort27aJatWr4+Pjw77//EhcXx549e2jQoIE9eQ0PD6dNmzacPXsW0zQ5fPgwHh4e3HfffXcdi0hmvP/euzQOqU7M0jFagUAcmlejZ9l16jJ9+w/I9D4xMTG0atWKgIAAqlWrxsyZM1m+fDk1atSgevXq9OrVi7i4Gy/sypYty7lz5wAYNmwYFStW5OGHH04zkjd69GiqVq2Kv78/nTt3vvsTzAZKXq/z0EMPMfuXn4he9BkJF05aHY5Iurwe6cHeC/E8+/wLd/3lXaJECVxcXDh27Bjh4eHUrVuX2rVrs379ejZv3kz16tWpXbs2O3fuJD4+3l6mUqVK7Nmzh/DwcOrVq5dFZyZya4ZhMG3qFEq6RBOzfrrV4YjcMcPJGe+WrzNzziLGjB2XqX0WL15MiRIl2LZtGzt37qRFixb07NmTmTNnsmPHDhITExmfwc29W7ZsYcaMGURGRrJw4UI2bdpkf2/EiBFs3bqV7du3M2HChLs+v+yg5PUmGjduzGcfDydq7jCSrkZbHY7ITRlOzniHvsrcpav4fOT/7vp49erVIzw83J6Y1q1b1/66fv36uLu74+fnR0REBBs2bKB27do3lBHJSR4eHixZOA+nA6uI2a27tsVxObl74d1mMG+98x5Lliy5Zfnq1auzdOlSBg0axJo1azhy5AjlypWjYsWKAPTo0YPVq1enu/+aNWt44okn8PT0pGDBgoSFhdnf8/f3p2vXrvzwww+4uLjc/cllAyWv6XihTx+e7vgEMQs/wUzSwvCSOzm5e+LdZjBDho1gwYIFd3Wsa/Ned+zYQbVq1ahTpw7r169P06tav359Vq9eTVRUFIULF6ZOnTr25FU9r2KF++67j2WLF3J1zWTiTmXdTYwiOc218P14txxI+85P8ddff2VYtmLFikRERFC9enXefvvtNDfT3q0FCxbQr18/IiIiCAkJyZUPx1HymoHR//ucGmWLE7NykuZUSa7l4lscn1Zv0Lnr0+zcufOOj1OvXj3mz59PkSJFcHZ2pkiRIly8eJH169fbE9N69eoxceJEAgICANsV+oYNGzh27BjVqlXLkvMRuV3+/v5M++5bouaPIPHyWavDEbljBUpXw73u0zwa+hjnz59Pt9ypU6fw9PSkW7duDBw4kPXr13PkyBEOHDgAwPfff0/Dhg3T3b9BgwbMnj2b2NhYoqKimDdvHmC7Yff48eM0btyYjz/+mEuXLhEdnftGoJW8ZsDZ2Znffp6J7+VDxETMtzockXS5l6xCgQa9aNqiJWfOnLmjY1SvXp1z585Rp06dNNt8fX0pWrQoYEteDx06RN26dQFwcXGhePHiBAcH4+Sk5kSsExYWxuBBrxM9bxjJ8bFWhyNyx7yqNyWuZDCPtWlLfHz8Tcvs2LGDWrVqERgYyPvvv8/QoUP59ttv6dChA9WrV8fJyYkXXngh3TqCgoLo1KkTAQEBPPbYY4SEhACQlJREt27dqF69OjVq1GDAgAG5cnUC43Z6FIODg81ra4XlJ0eOHKFGcC3cGvfD48Fgq8MRSVf02h8oHXeY9WtWUqBAAavDyVUMw9himma++gec39ps0zTp2r0nv0ccwLv1mxiGLqjEMZnJSUTN/5iWIRX5fso3GIZhdUg5LqM2W/+yM6Fs2bLMn/MbMUtHE3/2qNXhiKTLq/5THLvqTvdnntVUF8l3DMNgyuRJlC9oELP2e6vDEbljhpMz3i1eYd7SVXyndbJvoOQ1k+rXr8+Xo0cRNW8YSVcuWR2OyE0ZhhNezQewZN0Whg4bbnU4IjnOzc2NRfNm4358E9E79BhlcVxJUedIunKJhypUsDqUXEfJ623o3v1pXnimO9HzR2Am3v3C8CLZwcm1AN5t3uKT/43hl19+sTockRxXtGhR2wNnwqdy9fid38QoYpWk2MtEzR3GqM8/ydQyhIsXL6ZSpUpUqFCBESNG3LTMrFmzqFq1Kn5+fjz11FMAHD16lKCgIAIDA/Hz88u167peT3Neb1NycjJt2rbjzxMxeDUbkC/noYhjiPv7AFGz32f1H8sICgoCbDchVq9e3V6mc+fOvPnmm6xZs4YXXngBV1dX1q9fz7vvvsvChQtp2bIlXl5eeHt78/rrr2dpfKNGjeL555/H09MzS4+bHs15zX9+//132nXuim+nj3EtpKe/iWMwkxK4/NsQuoc9yhcjP79l+aSkJCpWrMjSpUspVaoUISEhTJ8+napVq9rL7N+/n44dO/LHH39QuHBhzpw5Q/HixYmPj8c0Tdzd3YmOjqZatWqEh4dTokSJ7DzFTMmozbZ89Vlvb29LlmF49913adCgAU2bNr2t/ZycnJg1/QeC69Tn7KZf8a7VjhPje+Hk5gGGgbNXYYq2ehVn78L/vz3lLuwizfpSoFQV4s8e5cKyiSRFnwczGS+/JvjW65wmET71bX9ci5Si2OODsvS8Jf9wv68CiY36ENqyNdsiNlOiRAk8PDyIjIy8oey0adN466236NatGwBfffUVFy5cwNnZmSFDhmRLfKNGjaJbt245lrzK3XO09jo0NJShQ97jvY+H4dtpBE7uXlnaXifF/Mv5RaNty3MlJ+HsW5x7O7yfHR+B5BOmaRKz4itCHirFyE8/ydQ+GzdupEKFCpQvXx6wdUrMmTMnTfI6adIk+vXrR+HChQEoXrw4YJtmc01cXBzJyclZdSrZyvLk1SoffPDBHe/r5eXFssULCAgK5koh29XJvV2G4+zpy7+rvuPShlkUadonzfZrkhPiOPvrhxRp3hePckEkJ1zl7G8fEb11AT5BrQFIOHcckpOJO7GL5PirOLnprnG5M16VHyb64imaPdaaTevX3rTM119/zaxZs/j9999ZtGgRUVFRREdHU7NmTd566600ZRs1asRnn31GcHAw586dIzg4mCNHjpCUlMSbb77JypUriYuLo1+/fvTp04eVK1cyZMgQihYtys6dO6lZsyY//PADY8aM4dSpUzRu3JiiRYuyYsWKnPg4xEHdTXv9n5f7s3PXLn5d9Dk+YYOBrGuvL66ZRoGygRQMfhyA+DOH7+IsRSBmyxwKxxzj15UbcHZ2ztQ+J0+epHTp0vbXpUqV4s8//0xTZt++fYDt/p2kpCSGDBlCixYtADh+/DitWrXiwIEDfPrpp7mi1/VWcs2c15UrV9KoUSPat29P5cqV6dq1q/1u6U2bNlGvXj0CAgKoVasWUVFRXL16lWeeeca+Ftm1L78pU6bQtm1bmjVrRtmyZRk7diwjR46kRo0a1KlThwsXLgDQs2dPfv75Z8C2msB7771HUFAQ1atXtz/Z4uzZszRr1gw/Pz969+5NmTJlOHfuHAAlS5Zk8YJ5XPnjS0j1BK4CpauR8O+pdM8zZvcq3EtWwaOcbRjXybUARZq9wKUNP/9/mT2r8KrWmAJlg4g9sCGrPmLJp7xqd+BvCtOn3wBiY2MJDAy0/8ycOZPevXsTFhbGp59+yrRp05g7d669h7ZTp06ZqmPy5Mn4+vqyadMmNm3axKRJkzh82PZFvnXrVkaNGsXu3bs5dOgQ69atY8CAAZQoUYIVK1YocXVAjtZeT/hyLFXv8yZm1TdpzuNu2+ukmAs4+xS1l3crXu6uP1vJv64c2Ejy9nks/30h3t7eWXrsxMRE9u/fz8qVK5k+fTrPPfccFy9eBKB06dJs376dAwcO8N133/HPP/9kad3ZIdckr3DzL7n4+Hg6derEF198wbZt21i2bBkeHh6MGzcOwzDYsWMH06dPp0ePHly9ehWAnTt38uuvv7Jp0yYGDx6Mp6cnW7dupW7dukxNZ8mJokWLEhERwYsvvshnn30GwPvvv0+TJk3YtWsX7du359ixY2n2CQ4O5puvJpIce5HEmIsAxB7ciFuxsvYy/0z/L6e+7c/pqa8CkHDuKG73pb1z0LXw/ZjxsSTHXQEg5q81eFVpgFfVBnpet9y9pASSo85Q9oHS9qT02k9mk9NbWbJkCVOnTiUwMJDatWtz/vx59u/fD0CtWrUoVaoUTk5OBAYGcuTIkSypU6zlSO21i4sL82f/gve53ZipHmBwt+21T41WnF80mr+nv8Wl8JkkRqX/RCSRjMSfOcyVpWNYOHc2DzzwwG3tW7JkSY4fP25/feLECUqWLJmmTKlSpQgLC8PV1ZVy5cpRsWJFext9TYkSJahWrRpr1qy58xPJIbkqeb3Zl9zevXu5//777U9/KFiwIC4uLqxdu9Y+P69y5cqUKVPG3i3euHFjfHx8KFasGL6+vrRp0wawPS0ovS/OJ598EoCaNWvay6xdu5bOnTsD0KJFC/tckdQ6duxAQZ+CnPn+VU598xLJcVfwrdPB/v69XYZT4pkx3N99ZKY+g7jT+3H2KIhLweIUKBNA/D+HSIqNytS+ItczTZPoZeN4OKAi77/3zl0fz8XFxT4n6lryca2eMWPG2JPiw4cP07x5cwDc3d3t5ZydnXPlc7Ll9jlae12oUCGW/b4QMz6G01Nf5dS3/e+6vfYoX5OSfb7GJyCUhAsnOD3lZS2lKLctKeZfouYNZ8K4MWmecJhZISEh7N+/n8OHDxMfH8+MGTMICwtLU6Zt27asXLkSgHPnzrFv3z7Kly/PiRMniI21XdD9+++/rF27lkqVKt31OWW3XJW8ZtWXXOrjODk52V87OTmle8xrZe6k3kKFfHnsscfwua8s97R6FacC6Xf3uxZ9gPi/D6TZlnDxbww3D5zcPYnZs4qE8yc4Mb4XJyc+R3L8Fa7sXXdb8YhcE7PxZ+4zLzBz2vdZ8vjWsmXLsmXLFgD7MC7YbowZP348CQm2JeT27dtHTExMhsfy8fEhKkoXZo7KEdvrChUqULxYMdyS4ygWNoiirV+7q/YawNnDB6+qjSja+jXc739IS3PJbTET44meP4J+z/WiW7eud3QMFxcXxo4dS2hoKFWqVKFjx474+fnx7rvvMnfuXMDWRt9zzz1UrVqVxo0b8+mnn3LPPfewZ88eateuTUBAAA0bNuT1119PsyJNbpWrktebqVSpEqdPn2bTpk0AREVFkZiYyCOPPMK0adMA2xflsWPHsvxqoX79+syaNQuwDYv++++/6ZadMG4MpVxjiFk/PcNjelVtRNyJ3cQeiQRsNwT8u2wiBWu3wzSTufLXWu7vNZZSL35DqRe/ofiT7xCzZ3WWnZPkHzF718GeJSxdNN9+R//1c17ffPPN2zrm66+/zvjx46lRo4Z9PiFA7969qVq1KkFBQVSrVo0+ffrcMql4/vnnadGiBY0bN779k5NcyRHa6wIFCjB0yDu2B87cYlQro/YaIPboNpITbCMQyXFXSPj3NC4Fi2XRGUleZx8ZC6zE8KF3flMiQMuWLdm3bx8HDx5k8GDbjYkffPCBvQfWMAxGjhzJ7t272bFjh32UolmzZmzfvp1t27axfft2nn/++bs7qRyS61cbcHNzY+bMmfTv35/Y2Fg8PDxYtmwZffv25cUXX6R69eq4uLgwZcqUNFfwWeG9996jS5cufP/999StW5f77rsPHx+fm5YtUKAASxbOwz8omJhCJfGq2vCm5Zxc3SnW7h0uLJ3AhaXjITkZL7/G+AS1Ju74Tpy978HF5x57effSfiScP0Zi9AVcvItk6flJ3hX39wFiV0xg9R/L0tw5mpSUdNPyU6ZMSfM69XJIqZfKqly5Mtu3b7e/Hjp0KGDrJRs+fDjDh6d9qlejRo1o1KiR/fXYsWPt/9+/f3/69++f6XOS3M9R2uvu3btz+NgJpi34BJ+272E43/yrMKP2GiD+7wP8u3QCODmDmYxPQHPc76+YpecleVfMnz+ljIzNzZKRsfxEDynIQFxcHM7Ozri4uLB+/XpefPHFm66Rmdr27dup37Ax3m0G414i988bkbwnMeo8l2a8wZSvxtGuXTurw8k19JCCvO122+ukpCRCW7Yh8jx4PdpXD5yRHBXz11rY8J19DW65Ua5+SEFuduzYMTp27EhycjJubm5MmjTplvv4+/sz7btveapnb5w7f4xLweI5EKmITXLCVaLnDeeNV/orcZV85Xbba2dnZ377eSZBtepyMWIOXjXb5kygku/Fnd5P7MqJN4yMSeapnzoDDz30EFu3bmXbtm1s2rTJfgftrbRp0wa/yhWJ1VxVyUGmmUzMktE0r1+Ttwf/N0uOeavnZR87dozGjRtTo0YN/P39WbhwYZbUK3K77qS99vHxYdj77xIbMTcHIhSBxKhzRM8fwXeTJ9kf2y23T8lrNhj20Qj2njiLV43WVoci+UhM+HQeKBDH91O+yZIh0KSkJPr168eiRYvYvXs306dPZ/fu3WnKDB06lI4dO7J161ZmzJhB375977pekZxy+PBh+vR9Ca9H+1kdiuQDtpGxjzQylgWUvKbo1asXxYsXp1q1and1nN9++40Rn4/Cu81/9VhXyTExu1fgcngtvy+Ym2U3wqR+Xrabm5v9edmpGYbB5cuXAbh06ZKGwCRHZEV7ffnyZR4NbYlL0JN4lK+ZhdGJ3Mg0k4n5/YssHRnLz5S8pujZsyeLFy++q2Ns3bqVp595Fp/Wb+GS6pGBItkp7uQerq75lmWLF1K8eNbNsb7Z87JPnjyZpsyQIUP44YcfKFWqFC1btmTMmDFZVr9Ieu62vU5MTCTsyQ5c8n0QT42QSQ6ICf+RBzzis2xkLL9T8pqiQYMGFCly50tRnT59muYtW+PRqA/u9z+UhZGJpC/x0hmiFnzCjB+m3vWowZ2YPn06PXv25MSJEyxcuJCnn37a/gQukexyt+11/5dfZfvRc3g3ek6JhGS7mF0rcDm8LktHxvI7Ja9ZIDY2lmaPtcas2ASvyg9bHY7kE8lxV4ieN4z3/juIVq1aZfnxM/O87MmTJ9OxY0cA6taty9WrV9M8vEAkt/ly/Him/TIHr1ZvpLu+q0hWsY2MfZPlI2P5nZLXu2SaJp26Ps1psxBedTpZHY7kE2ZyEtG/jySsWUNef+3VbKkjM8/LfuCBB1i+fDkAe/bs4erVqxQrpicMSe60fPly3njrbXzCBuOcwWNhRbKCfWRs2veWjIzlZUpe79Lb777Hmq1/4d3sJQ0/SY6JWTOVikXcmPzVhGz7u8vM87I///xzJk2aREBAAF26dGHKlCn6dyC50r59+3iyQye8H3sd1yIlb72DyF1IjrtC1Nyh2TYylt9pzOQu/PjjdEZP+JpCnT/BcHGzOhzJJ2K2L8Hj9FYWRmzC1dU1W+tq2bIlLVu2TLPtgw/+/xncVatWZd26ddkag8jdunDhAo82fwzX2l0oUMbf6nAkjzOTk4he/DltmzfKtpGx/E7Ja4ouXbqwcuVKzp07R6lSpXj//fd59tln0y3/559/8vyL/SjY7gOcvQrnYKSSn109tp34DdNYu37dXd2wIuLIbqe9TkhIoNXjTxJ7XwDeAS1yONL8LTH6AifHdU+zzb1oaTxLVMDZ3StTxzi3aX6WxOLk5kGRgEez5FjpMYECNZ/gauQCKt3jztfZODKW3yl5TTF9+vRMlz127BiPtXkcz2Yv4Va8XDZGJfL/Ev49RfSiz/lt1nQqVapkdTgilslse22aJs/1eZG9567i06Zn9gYlNzCcXXApUhIzMR4zMYHkK5eIO3ecuHP/fyNouXLlCA4Opnbt2nh4eNxwjH5ZlLw+0aYlTZo0yZJjpef7aT+yde4nFHWNZ0HE5mwfGcvPlLzepujoaJq2aIVT9VZ4VqhtdTiSTyRfjSZ67jA+HvYBzZo1szocEYcw8n+j+O33Ffh2HIHh5Gx1OPmOs0dBSj438YbtZnISCeePE//3QS6dO8gfm/eweOlyevbowYCX+lKhQgV7WUd6at+WyO1s3z6d5Zs3amQsmxmmaWa6cHBwsLl58+ZsDCd3S05O5rHWj7P5nwS8m+oGLckZZnISUbM/pGPTOkz4Ug8BuFOGYWwxTTPY6jhyUn5usxcuXEjHrj0o1PljXHzvtTocuYWEi38Tt30xsbuWExQUxMBXBtCqVSucnR3nouPQocNcvPgvQUFBVoeSJ2TUZmu1gdvw2sA32LTvON5NXlDiKjkmZuXX+JcuzNjR/7M6FBGHsHPnTjp1fRrvVm8ocXUQroXuw7tBT4o+N5k9XgE885//cn+pMvxv1GiSkpKsDi9Typcvp8Q1hyh5zaTJ33zD5B9m4tVqEIaz5rFIzojZugCff/cx59efcHHRLB+RWzlz5gxNW7SkwMPPUKBUVavDkdtkuLjhXa0JPh0/xjl0IB+O/YbgOvXZt2+f1aFJLqLkNRPWrFnDgFcH2ha29ihodTiST8QejiBxy88s/30hvr6+VocjkuvFxcXRovXjJJath5dfY6vDkbvkfl8FfNp9yPHCgdQIqc0nn37uML2wkr2UvN7CoUOHaPNEO7xD/4PrPaWtDkfyiYRzx4n5fRRzf/uF8uXLWx2OSK5nmibdn3mWo1dc8arf1epwJIsYhhPeQW0o1PlTRkyYql5YAZS8ZujSpUu2ha1rtsejnOaxSM5IunKJqHnDGPO/z3nkkUesDkfEIQwd/hG/r92EV/OXMQx9teU1roXvx6f9UI4XCiQopDbTp8+wOiSxkCbRpSMxMZE2T7TncpFK+NTQo90kZ5hJCcQs+ITeT3fmmWd6Wh2OiEP49ddf+XjkFxTq/AlObgWsDkeyiWE44V2zDW6l/Xiu3wAuXb7EC336WB2WWECXp+noN+Bldp68iHej3laHIvmEaZrELJ9ArcoP8NknH1sdjohD2Lp1K9179can9Vu4+BS1OhzJAW7Fy+PbfigD336fjz/9zOpwxAJKXm9i7LhxzPhtAd6Pva6FrSXHXNk8myJXT/LzzB9xctI/TZFbOXXqFM1btsajUR/c73/I6nAkB7kWKYlvh+EM/fR/jPvyS6vDkRymb8jrLF26lDfffg/vsLdxKuBtdTiST1zZ/yfJO+az/PeFeHvr707kVq5cuUKzx1pjVnoUr8oPWx2OWMClYDEKPvE+g955n2nTfrQ6HMlBmvOayt69e2nfqQvej72Oa+H7rQ5HbuHi2h+5tC57GqxCFWvh6pP9j/dzur8KTsXKcWX5OP5YsojSpbWihcitJCcn07lbd/6mMN61O1odjljItfD9FGz7Ln369ady5UrUrFnT6pAkByh5TXH+/HkeDX0MtzpdKfBAdavDkUxwLVYGt3sfJDkhjsR/T8JtPOr4Vt5+viOFChXKsuPdzLFjxxg57isAJo0fR61atbK1PpG8YvA777Em4i8Ktv9QTzsU3IqVpcAjz9DxqafZvX0r7u7uVock2UzJKxAfH0/LsCe4WiIIb//mVocjmeRVqT5eleqn2WYmxhN/9gjxBzcSt2sZVatW4Y1XBvD444/j6pq7nowWGRnJBx98wOB33qVLl85WhyPiEKZN+5GxX02mUOdPMFzcrA5Hcgmvqo24eGg977z3Pp+MGG51OJLN8v2cV9M06d3nRQ78m4TXw92tDkfukuHihvv9FfF5uBv39J7E4aJ16fPmUMpWqMjq1autDi+NypUr88svv/Lh+0OsDkXEIWzYsIE+fV/Cp81/cfYqbHU4kosYhoFn4xf5csJENm/ebHU4ks3yffL62ecjmbt0NV6h/9HKAnmM4eyKV5UG+LQfRkLw07Rs254X+vYnJibG6tAAKFCgAE8++YSGPUUy4dixY7Rs8ziezV7CrXg5q8ORXMjZuzAFHnmGTl27ExcXZ3U4ko3ydfI6f/583h8+Au82/8XJ3dPqcCQbeT5Um8JPj+bn8D1UrFKNiIgIq0MSkUyKjo6maYtWOPm3xrNCbavDkVzMq2ojLroU5p333rc6FMlG+TZ53bFjB12e7oFPqzdx8S1udTiSA5w9fPBu8QrxQZ1p9GgzVq1aZXVIInILSUlJPNmhM+c9SuIV/ITV4UguZxgGHo36MHbcOC5evGh1OJJN8mXyeubMGZo91ooCj/TCvWRlq8ORHOZV+WE8WrxGq8efUAIrksu9NvANNu0/gXeTFzTFRjLFxecePMvXZMp331kdimSTfJe8Xr16ldBWYSSWexivqo2sDkcs4lEmAK8Wr9Gm7ZOaQiCSS309eTLfTJuFd6tBGM65a7UQyd1cqrfg8y/GYmbhEoqSe+Sr5NU0TZ7u2YtjV93xqtfF6nDEYh5lA3F/5Fkef7IDV65csTocEUll9erVvPzaG/iEDcbZo6DV4YiDcS9ZlegEWLFihdWhSDbIV8nrh0OHsWz9VryaD8Aw8tWpSzq8qjbkSsEHeOOtwVaHIiIpDh48SNgT7fAO/Q+u9+ipc3L7DMPAqWpzPv3faKtDkWyQbzK4n376mU9GjcWr9Vs4uRawOhzJRTwbPceUqT8QHh5udSgi+d6lS5doGtoSl+AOeJQLsjoccWBeVRuxauUKTpw4YXUoksXyRfK6ZcsWnnnueXzavIWLzz1WhyO5jLOnL56Nnqdztx4kJSVZHY5IvpWYmEibtu24XKQS3jVaWR2OODgnd088K9bl119/tToUyWJ5Pnk9deoUoa3a4Nn4Rdzvq2B1OJJLeVaqRwzuLFiwwOpQRPKtvv0HsPPUJbwb9bY6FMkjzGIPsXbDJqvDkCzmYnUA2enKlSs0bdEKqjTHs1I9q8ORXM6oGsqn/xtNWFiY1aGI5DvTZ8xg0oTxFAxqzeV1P1odDgAXw2fe8b4+Qa1wLuCdhdGACVxKicmrakNcC92XpcfPaiZguLjj7FMUz4p1cXLzyPEY3O6rwKZ1E3K8XsleeTp5fWvwO+z7aw9FqhXjytIvrA7H7mzE0jve16vEQ3jel/WPRjwbscT+/4UeCsHVxxGeG26Q7FIAs0gZvKo0wMnt7uYye1V+mIhJ33LgwAEqVFAvvUhOql2rFh9++KHVYaTxzl0kr881qkShQoWyLhggOTmZISlT82N229aoDgoKom3btllaT1YxTZOomBgiInewftLXeFZpRIGabXP0wUBuxcpw8vgRYmJi8PLyyrF6JXsZt7MGWnBwsLl58+ZsDCdr/fXXX/z5559Wh3GDnj173vG+w4YNo0SJElkXTIq5c+fy22+/pdn2yiuv4O/vn+V1ZZXk5GQuX77M/MVLCQ9fh2fVxnjU6oCzp+8dHzN69RS61XmAUSM/z8JIJTcwDGOLaZrBVseRkxytzZbM27VrF7Vr1yYmJobXX3+dTz/91OqQMnT8+HFGjxnL+Enf4NPmLdxLVMqxuqNnvM7cHydTr55GYB1JRm12nk5e5c6sX7/e/o98/vz5tGqV+2+cOHr0KMM++pgZv83D54khdzycdvXELorsmMnubVuyOEKxmpJXyYs+/vhj3nzzTb744gsGDBhgdTi3NHfuXJ7q/gy+HYbjek+pHKkzevl43n26uUN8PvL/Mmqz8/wNW3L76tati2maBAcH07p1a3bs2GF1SLdUpkwZvprwJR8OfoOo34aQdOXSHR3HrXh5Du7bQ0JCQhZHKCKS9QYNGkTbtm15+eWXHeJpUmFhYbz/7tvErf02x+o07ynHmvUbc6w+yX5KXiVdmzbZ7tDMzVMHrvfygP60b/MYsZt/u3Xhm3By88CzyH3s2rUriyMTEckeP/30EwDjxo2zOJLMefGFPsSd3kfCxb9zpD4Xn6KcOp0zdUnOUPIqGVq7di0AZ86csTiSzBv81iCu7FpOckLcHe3vWvxBtmzRtAERcQwuLi4UK1aM/v37Wx1Kpnh6etKje3eu7rzzm5dvh+HiRuzV2BypS3KGklfJUP369QFy7d2sN/Pggw9Ss2ZNruy7sydmxXvfy8GDh7I4KhGR7DNo0CCrQ7gtjRs+gsulkzlSl+HixtWrV3OkLskZSl7llkqWLMn69eutDuO2PNa0Mcnnj93RvoaLOzGxukoXEcdRsmRJq0O4LQULFsRMyKGE0tmF+Lj4nKlLcoSSV7ml0qVLWx3CbfP19cU58c4SUMPFlZgrSl5FxHE4OTnW17lhGHmyLskZjvXXLpZwdna2OoTb5mgNuYiIiGSOvuFFRERExGEoeRURERERh6HkVUREREQchpJXEREREXEYSl5FRERExGEoeRURERERh6HkVUREREQchpJXEREREXEYSl5FRERExGEoeRURERERh6HkVUREREQchpJXEREREXEYSl5FRERExGEoeRURERERh6HkVUREREQchpJXERERyRHHjh3jnmL30rZde0x3nxyp03D34uihAxQqWpyuPZ7JkToleyl5lTzl77//5uuvv2bdunV3dZy9e/bw9ddf89dff2VRZCIiUrx4cZo8+iimmzdudbvmSJ2uhe6jYP2nSIiPp0fXLjlSp2QvJa+Sp2zctJnnnnuO2au24FLx4Ts6hkf5YLb/HcNzzz3Hr7PnZnGEIiL5V4ECBZg1fRpvDHiBSz8NJu7vA9lan5mcRPTqb/E4uIItGzfQvHnzbK1PcoaSV8lTwtq0ZtasnyDqH5JiLt/RMZJjL5N07ijDPhrBW4MGZnGEIiL5m2EYvPfO20yZOJao2e9zZW94ttSTHB9L9PyPedA4w7aITVSuXDlb6pGc52J1ACJZrUOH9pQvX47QVm2IvngSr1rtMQwjU/vG7F7F1dVfM+27b3n88cezOVIRkfyrffv2lC9fnuYtW992W30riZfPED13OG2aPsLkrybg5uaWJceV3EE9r5In1axZk+0Rmyl+fhvRv3+BmZiQYXnTNIkO/xFj84+sXbVCiauISA4ICgq6rbY6M+JO/sXFGYP478t9mPrtZCWueZCSV8mzSpQoweYN66hb1ofLv7xDUszFm5ZLTogjetHnlIzZz/aIzQQEBORsoCIi+ViatvrX9NvqzIjZvYro+cOZ/t1k3hj4epb15EruouRV8jRPT0/m/vozz3cO4+KMgcSfPZLm/cToC0T98jYNKxXnz3Wrue+++6wJVEQkH7O31Z3CuDjjjRva6lsxzWSiw3/Eact01q1aQZs2bbInUMkVlLxKnufk5MTHHw3jy5GfcPmXd4k9uAmA+H8OcWnGG/Tv2YlfZs2gQIECFkcqIpJ//X9b/TGXfn6HKylt9a0kJ1wletHnlLpygO0Rm/H398/mSMVqumFL8o2nn+5GhQoP0iqsLVcPbSLp4Hq+njCezp07WR2aiIikSN1Wx1w8iWfQ4+kO/ydGXyBm3kc8WieAH777Vp0Q+YR6XiVfqVu3Lls3b6T5Q76sWPq7ElcRkVzoWltd6Hg4Mcu/xEy68Uau+H8OcmnGQF5+tgs/zfhRiWs+ouRV8p0yZcrw04wfCQkJsToUERFJR5kyZYjYtIHAe+Dyb++TFBtlf+/KvnAu/zaEyV+OYci77+jGrHxGyauIiIjkSj4+Pvy+cB5Pt2nCpRkDSTh/guiNP5Mc/i0rly2hY8cOVocoFtCcVxEREcm1nJ2dGf2/z6lWtTL9B/yHBx+qyNItmyhZsqTVoYlFlLyKiIhIrvf8c88R2rw5xYoVw9PT0+pwxEJKXkVERMQhlClTxuoQJBfQnFcRERGxxCuvvMKoUaPsr0NDQ+ndu7f99WuvvYZhGMyePdu+rVKlSgwdOtT+ul27dsycOZPAwEACAwPx9vamUqVKBAYG0r1795w4DclhSl5FRETEEvXr1yc8PByA5ORkzp07x65du+zvh4eH8+GHH9rLnD9/Hi8vL9avX28vs379eho2bEhkZCSRkZEEBwczbdo0IiMjmTp1as6ekOQIJa8iIiJiiXr16tkT0V27dlGtWjV8fHz4999/iYuLY8+ePTRo0MCevIaHh9OmTRvOnj2LaZocPnwYDw8PPdo7n1HyKvlKo0aN2Lx5MwAtW7bk4sWLN5QZMmQIn332WQ5HJiKS/5QoUQIXFxeOHTtGeHg4devWpXbt2qxfv57NmzdTvXp1ateuzc6dO4mPj7eXqVSpEnv27CE8PJx69epZfRqSw3TDluRbCxcutDoEEZF8r169eoSHhxMeHs6rr77KyZMnCQ8Px9fXl/r16+Pu7o6fnx8RERFs2LCBN954g0OHDhEeHs7WrVupX7++1acgOUw9r+LwYmJiaNWqFQEBAVSrVo2ZM2eyfPlyatSoQfXq1enVqxdxcXE37Fe2bFnOnTsHwLBhw6hYsSIPP/wwe/futZcZPXo0VatWxd/fn86dO+fYOYmI5BfX5r3u2LGDatWqUadOHdavX5+mV7V+/fqsXr2aqKgoChcuTJ06dewJr3pe8x/1vIrDW7x4MSVKlGDBggUAXLp0iWrVqrF8+XIqVqxI9+7dGT9+PP/5z39uuv+WLVuYMWMGkZGRJCYmEhQURM2aNQEYMWIEhw8fxt3d/aZTDERE5O7Uq1ePzz77jPLly+Ps7EyRIkW4ePEiu3btYtKkSfYyr732Go0aNQLA39+fDRs28M8//1CtWjULoxcrqOdVHF716tVZunQpgwYNYs2aNRw5coRy5cpRsWJFAHr06MHq1avT3X/NmjU88cQTeHp6UrBgQcLCwuzv+fv707VrV3744QdcXHStJyKS1apXr865c+eoU6dOmm2+vr4ULVoUsCWvhw4dom7dugC4uLhQvHhxgoODcXJSKpPf6NtYHF7FihWJiIhg4cKFvP322zRp0iTLjr1gwQJWr17NvHnzGDZsGDt27FASKyKShZydnbl8+XKabVOmTEnzunjx4pimmWbbypUrb3q89LZL3qHLFXF4p06dwtPTk27dujFw4EDWr1/PkSNHOHDgAADff/89DRs2THf/Bg0aMHv2bGJjY4mKimLevHmAbc3B48eP07hxYz7++GMuXbpEdHR0jpyTiIiI3Jy6kMTh7dixg4EDB+Lk5ISrqyvjx4/n0qVLdOjQgcTEREJCQnjhhRfS3T8oKIhOnToREBBA8eLFCQkJASApKYlu3bpx6dIlTNNkwIABFCpUKIfOSkRERG7GuL4bPiPBwcHmtTUyJf94+OGHWbdu3Q1DNiKOxDCMLaZpBlsdR05Sm51/zJo1i06dOqmdljwjozZb0wZEREQkWzk7OxMYGGj/GTFiBGC7YdbPz4/AwEBiY2MZOHAgfn5+DBw4MNseGDNq1CiuXLmS5ceVnKNpAyIiIpKtPDw8iIyMvGH7tGnTeOutt+jWrRsAX331FRcuXMDZ2ZkhQ4ZkSyyjRo2iW7dueHp6ZsvxJfup51VERERy3Ndff82sWbN455136Nq1K2FhYURHR1OzZk1mzpyZpmzqR3ufO3eOsmXLArZ7EwYOHEhISAj+/v5MnDgRsK040KhRI9q3b0/lypXp2rUrpmkyevRoTp06RePGjWncuHGOnq9kHfW8Sp62ePFiXn75ZZKSkujduzdvvvlmmvdHjhzJ119/jYuLC8WKFeObb76hTJkyFkUrIpI3xcbGEhgYaH/91ltv0bt3b9auXUvr1q1p3749AN7e3vYe2sz0vE6ePBlfX182bdpEXFwc9evXp3nz5gBs3bqVXbt2UaJECerXr8+6desYMGAAI0eOZMWKFfY1ZMXxKHmVPCspKYl+/fqxdOlSSpUqRUhICGFhYVStWtVepkaNGmzevBlPT0/Gjx/PG2+8ccMVv4iI3J30pg3crSVLlrB9+3Z+/vlnwPaExf379+Pm5katWrUoVaoUAIGBgRw5coSHH344y2OQnKdpA5Jnbdy4kQoVKlC+fHnc3Nzo3Lkzc+bMSVOmcePG9nlPderU4cSJE1aEKiJyRy5evIiHlxedOnUCwMXVjRatw26xl+NxcXEhOTkZgKtXr9q3m6bJmDFjiIyMJDIyksOHD9t7Xt3d3e3lnJ2dSUxMzNmgJdsoeZU86+TJk5QuXdr+ulSpUpw8eTLd8pMnT+axxx7LidBERLKEr68vAwe+gVfhYvjWegLvgoV4561BVoeV5cqWLcuWLVsA7L2sAKGhoYwfP56EhAQA9u3bR0xMTIbH8vHxISoqKvuClWyn5FXStX37DlqGPcm6desAaBn2JNNmzLI4quzxww8/sHnzZgYOHGh1KCIimWYYBh8MeY/J48cS5BtL5JaN1K9f3+qwbnBtzuu1n+vvP7iV119/nfHjx1OjRg3OnTtn3967d2+qVq1KUFAQ1apVo0+fPrfsYX3++edp0aKFbthyYHpIgaRr+/btNGjSlKvuhXHyLET8kQimTp3KU126WB1apqxfv54hQ4bw+++/A/DRRx8BthsFUlu2bBn9+/dn1apVFC9ePMfjlJyhhxSIiDgOPaRA7oi/vz/bIzbzgK8brucPsuKPPxwmcQUICQlh//79HD58mPj4eGbMmEFYWNq5YFu3bqVPnz7MnTtXiauIiIgD0GoDkqEHHniA7RGbiI+Pp2DBglaHc1tcXFwYO3YsoaGhJCUl0atXL/z8/Hj33XcJDg4mLCyMgQMHEh0dTYcOHQDb+c6dO9fiyEVERCQ9mjYgIvmCpg2I5G63Wpd7woQJjBs3DmdnZ7y9vfnqq6+oWrUq06ZN49NPP7WX2759OxEREWnWlRXHk1GbreRVRPIFJa8iuVdSUhIVK1ZMsy739OnT06zLffnyZfsI4Ny5c/nyyy9ZvHhxmuPs2LGDtm3bcvDgwRyNX7Ke5ryKiIhIrpWZdblTT12LiYnBMIwbjjN9+nQ6d+6c7fGKtTTnVURERCx1s3W5//zzzxvKjRs3jpEjRxIfH88ff/xxw/szZ868IemVvEc9ryIiIuIQ+vXrx8GDB/n4448ZOnRomvf+/PNPPD09qVatmkXRSU5R8ioiIpLNGjdubF9z+ppRo0bx4osv3rR82bJl0yzGfzPDhw9P87pevXoAHDlyxJ7Abd68mQEDBtxp2HY9e/ZM82SrrFayZEmOHz9uf33ixAlKliyZbvnOnTsze/bsNNtmzJhBFwdazlHunJJXERGRbNalSxdmzJiRZtvdJlvXJ6/h4eE3lAkODmb06NF3XEdOycy63Pv377f//4IFC3jooYfsr5OTk5k1a5bmu+YTSl7zAF3Ri4jkbu3bt2fBggXEx8cDtrb01KlTnDx5kurVq1OtWjUGDRp0033btm1LzZo18fPz46uvvgLgzTfftD9ytWvXrgB4e3vfsO/KlStp3bo1AC1btrQ/ntXX15fvvvuOpKQkBg4cSEhICP7+/kycOBEA0zR56aWXqFSpEk2bNuXMmTNZ/pmklnpd7ipVqtCxY0f7utzX1t4eO3Ysfn5+BAYGMnLkSL777jv7/qtXr6Z06dKUL18+W+OUXMI0zUz/1KxZ05TcZ+LEiWbPnj3TbKtdu7a5atWqm5YvU6aMefbs2QyP6eXlddPthw8fNv38/O4s0HT06NHD/Omnn7L0mCLXAzabt9He5YUftdm5S6tWrczZs2ebpmmaH330kfnMM8+YpUuXNs+cOWMmJCSYjRs3Nn/77TfTNNO20+fPnzdN0zSvXLli+vn5mefOnTNN88Z2+trr1O30ihUrzFatWqUpt3nzZrN69ermxYsXzYkTJ5offvihaZqmefXqVbNmzZrmoUOHzF9++cVs2rSpmZiYaJ48edL09fVVOy05KqM2Wz2veYCu6EVEcr/UUwdmzJhBmTJlaNSoEcWKFcPFxYWuXbuyevXqG/YbPXo0AQEB1KlTh+PHj6cZPr9d586d4+mnn+bHH3/E19eXJUuWMHXqVAIDA6lduzbnz59n//79rF69mi5duuDs7EyJEiVo0qTJHdcpktW0VFYeUKRIEWrVqsWiRYt4/PHHmTFjBk2bNmXQoEFs2bKFwoUL07x5c2bPnk3btm3T7PvNN99QpEgRYmNjCQkJoV27dowYMYKxY8cSGRmZ6RgWLlwIwJYtW3jmmWdo27YtkydPxtfXl02bNhEXF0f9+vVp3rw5W7duZe/evezevZt//vmHqlWr0qtXryz8REREcp/HH3+cV155hYiICK5cuUJgYOAtF9NfuXIly5YtY/369Xh6etKoUSOuXr16R/UnJSXRuXNn3n33Xfv0L9M0GTNmDKGhoWnKXmvTRXIj9bzmEbqiFxHJ3by9vWncuDG9evWiS5cu1KpVi1WrVnHu3DmSkpKYPn06DRs2TLPPpUuXKFy4MJ6envz1119s2LDB/p6rqysJCQmZrv/NN9/E398/zU1NoaGhjB8/3n6cffv2ERMTQ4MGDZg5cyZJSUmcPn2aFStW3OXZi2Qd9bzmEbqiFxHJ/bp06cITTzzBjBkzuP/++xkxYgSNGzfGNE1atWrF448/nqZ8ixYtmDBhAlWqVKFSpUrUqVPH/t7zzz+Pv78/QUFBTJs27ZZ1f/bZZ/YbngA++OADevfuzZEjRwgKCsI0TYoVK8bs2bN54okn+OOPP6hatSoPPPAAdevWzdLPQeSupDcZ9mY/mvyfu3Xs2NEMCAgw3333XfPUqVPmAw88YJ49e9ZMTEw0H330UfuNAtduBJg9e7bZunVr0zRNc8+ePaa7u7u5YsUK0zRNs1ChQmZ8fLz92Le6EeD11183X3nllTTxTJw40Xz88cftx9m7d68ZHR1t/vLLL2bz5s3NxMRE89SpU2ahQoV0I4BkO3TDloiIw8iozVbPax6iK3oRERHJ6wxbcps5wcHB5ubNm7MxHBGR7GEYxhbTNIOtjiMnqc0WEUeVUZutG7ZEREQc3OLFi6lUqRIVKlRgxIgRVocjkq2UvIqIiDiwpKQk+vXrx6JFi9i9ezfTp09n9+7dVoclkm2UvIqIiDiwjRs3UqFCBcqXL4+bmxudO3dmzpw5Voclkm2UvMotaThKRCT3OnnyJKVLl7a/LlWqFCdPnrQwIpHspeRVMqThKBEREclNlLxKhjQcJSKSu5UsWZLjx4/bX584cYKSJUtaGJFI9lLyKhnScJSISO4WEhLC/v37OXz4MPHx8cyYMYOwsDCrwxLJNnpIgYiIiANzcXFh7NixhIaGkpSURK9evfDz87M6LJFso+RVMqThKBGR3K9ly5a0bNnS6jBEcoSmDUiGNBwlIiIiuYl6XiVDGo4SERGR3ETJq9yShqNEREQkt9C0ARERERFxGEpeRURERMRhKHkVEREREYeh5FVEREREHIaSVxERERFxGEpeRURERMRhKHkVEREREYdhmKaZ+cKGcRY4mn3hiIhkmzKmaRazOoicpDZbRBxYum32bSWvIiIiIiJW0rQBEREREXEYSl5FRERExGEoeRURERERh6HkVUREREQchpJXEREREXEYSl5FRERExGEoeRURERERh6HkVUREREQchpJXEREREXEY/wfKeBCrHpplfwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "WWT.plot_sankey_report(margin=.6, gap=.6)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "A=WWT.products()['WWTEffluent']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Effluent and solids mass flows"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Effluent volume (Liter) : 988.4007936272117\n",
      "Screen Rejects mass (kg): 0.2857142857142857\n",
      "RawWWTSoilds mass (kg): 0.0\n",
      "DewateredWWTSolids mass (kg): 0\n",
      "DeriedWWTSolids mass (kg): 3.6033275744680844\n",
      "\n"
     ]
    }
   ],
   "source": [
    "products = WWT.products()\n",
    "print( \"\"\"\n",
    "Effluent volume (Liter) : {}\n",
    "Screen Rejects mass (kg): {}\n",
    "RawWWTSoilds mass (kg): {}\n",
    "DewateredWWTSolids mass (kg): {}\n",
    "DeriedWWTSolids mass (kg): {}\n",
    "\"\"\".format(products['WWTEffluent'].vol,\n",
    "           products['WWTScreenRejects'].mass,\n",
    "           products['RawWWTSolids'].mass,\n",
    "           products['DewateredWWTSolids'].mass,\n",
    "           products['DriedWWTSolids'].mass))        "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Sensitivity to the Mixed liquor suspended solids (MLSS)  (Default: 1000 mg/L)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, 'Percent of Incoming PFAS that \\n remains in the Effluent (%)')"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZoAAAEGCAYAAABcolNbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAA5B0lEQVR4nO3dd3xUZfb48c9JI5RA6CSkUAWpIQQSQF1XrKuAIioaFHVtq+u665av/nR7dVfddRfLYluVIIqAYC/Y0QQJvUgRSAg10otAyvn9cW90iJnJhMzkppz36zWvzNw2ZyaTnHme+9zziKpijDHGhEuE1wEYY4xp3CzRGGOMCStLNMYYY8LKEo0xxpiwskRjjDEmrKK8DiDUOnTooN26dfM6DGOMaVDy8/O/UtWO4Th2o0s03bp1Y9GiRV6HYYwxDYqIFITr2NZ1ZowxJqws0RhjjAkrSzTGGGPCyhKNMcaYsLJEY4wxJqws0RhjjAkrSzTGGGPCqtElmm37v+bI8VKvwzDGGONqdIlm96HjTHj0M4r2HvE6FGOMMTTCRNOtfUu27D3CuCkLyNu42+twjDGmyWt0iSYuNoqXbxtFm+bRZD+Rx7TcsFVVMMYYE4RGl2gAenZsxZzbRnFa7w7c+/JK7pmzguOl5V6HZYwxTVKjTDQAbZpH8+TkYdz8vR7k5BUy6ck8dh865nVYxhjT5DTaRAMQGSHcfcGp/OuKNJZt2cfYKQtYtW2/12EZY0yT0qgTTYWLh3Rl5i0jKCtXJjz6Ga8t3+51SMYY02Q0iUQDMCgpnnm3j+LUhDhum76YB95eS3m5eh2WMcY0ek0m0QB0iovl+ZuyuDwjif+8t4Gbnsvn4NESr8MyxphGrUklGoBmUZHcd+kgfjemH++v3cU1Ty3kaEmZ12EZY0yj1eQSDYCIcO2o7ky5cghLCvdxz5yVqFo3mjHGhEOTTDQVLhiYwE/P7s2sxUU8+ckmr8MxxphGqUknGoCfnNWbCwZ04S+vr+HDdcVeh2OMMY1Ok080ERHC/ZcN5pTOcdw+fTEbiw95HZIxxjQqTT7RALRsFsXj12QQFRnBDc8u4oCNRDPGmJCxRONKbteCR7LTKdx9hDueX0KZXWNjjDEhYYnGR1aP9vxubH/eX1vMP95a63U4xhjTKER5HUB9MykrlTXbD/DYh1/St0scFw/p6nVIxhjToFmLpgq/HdOf4d3b8X+zlrO8aJ/X4RhjTINmiaYKMVERPJqdTodWzbjp2Xx2HTjqdUjGGNNgWaLxo32rZjx+TQb7vy7h5mn5VqbGGGNOkiWaAPoltubBywezpHAf975sZWqMMeZkWKKpxgUDE7hjdG9eyi/iqQWbvQ7HGGMaHBt1FoQ7Rvfmix0H+PNrq2kZE8nE4Sleh2SMMQ2GtWiCEBEhPHh5GqN6deCu2Sv47dyVlJSVex2WMcY0CPUm0YjIHSKyUkRWichP3WXtROQdEVnv/mzrVXwtm0Xx9LXDuPH07jzzWQHXPLmQvYePexWOMcY0GPUi0YjIAOBGYDgwGLhIRHoBdwHzVbU3MN997JmoyAjuubAfD1w2mPzCvYx9+BO+2HHAy5CMMabeqzbRiMh9wSyrpVOBPFU9oqqlwIfAeGAc8Iy7zTPAxSF+3pNy6dAkXrx5BMdKyhn/yKe8uXKH1yEZY0y9FUyL5pwqll0Q4jhWAqeLSHsRaQH8AEgGOqvqdnebHUDnqnYWkZtEZJGILCourps5ZdKS43nl9tM4pXMct0zL51/vrqPcCnEaY8x3+E00IvIjEVkB9BGR5T63TcDyUAahqmuA+4C3gTeBpUBZpW0UqPI/uapOVdUMVc3o2LFjKEMLqHPrWGbclMWl6Un869313JqzmMPHSuvs+Y0xpiEI1KKZDowB5rk/K25DVXVSqANR1SdVdaiqngHsBdYBO0UkAcD9uSvUz1tbsdGR3H/ZIH59UT/eXr2DSx/9lC17jngdljHG1Bt+E42q7lfVzap6paoWAF/jtChaiUjILyQRkU7uzxSc8zPTcZLcZHeTycDcUD9vKIgIPzytO89cP5zt+48ydsonfPrlV16HZYwx9UIwgwHGiMh6YBPOSfrNwBthiGWWiKwGXgFuU9V9wN+Ac9znP9t9XG+d3rsjc28bRftWzbj6yYU88+lmK1tjjGnygqkM8CcgC3hXVYeIyPeBKrvORGSEu+50IAGnFbQSeA2Ypqr7/T2Jqp5exbLdwOggYqw3unVoyZxbR/KzF5by23mrWLP9AH8YN4CYqHoxktwYY+pcMP/9Stx/+BEiEqGq7wMZlTcSkTeAG4C3gPNxEk0/4F4gFpgrImNDFnk9FhcbzdSrM/jx93sx4/MtXPV4LsUHj3kdljHGeEKq69oRkXdxrl/5K9AB54T8MFUdWWm7Dqoa8MREMNvUVkZGhi5atCicT1Ejry7fxi9mLqNtiximXp3BwKQ2XodkjDHfISL5qvqdRkQoBNOiGYfTBfYznKHHX+KMPjtBVQlEREa753ii/W3T2F00KJFZPxpJhAgTHvuUuUu3eh2SMcbUqWoTjaoeVtUyVS1V1WdU9d9uV1pAIvIAMAqnpEy9HC1WV/ontmHuj0cxOCmeO2Ys5b43v6DMLu40xjQRwYw6G+8WtdwvIgdE5KCIfKfAl4g8ICLxPotSgD8Cf3bvN2kdWjVj2g2ZXJWZwqMffMmNzy7iwNESr8MyxpiwC6br7O/AWFVto6qtVTVOVVtXsd1sYIaI/EREIoFngfeBz4DHQxdywxUTFcFfLhnIny4ewEfrirnk4QVsLD7kdVjGGBNWwSSanW6JmIBUdYGqng/swRl5Jqp6pqpmqepDtQ20MZmUlcq0GzLZe6SEcQ8v4IO19a7ggTHGhEygWmfjRWQ8sEhEXhCRKyuWucsrbx8lIhfijEq7GBgsIvNEZHDYom/Asnq0Z+5to0hq24Lr//c5Uz/60i7uNMY0SoEu2PQdWXYEONfnseJ0lfl6GaebrAWQraqTRSQR+IOIqKreGIJ4G5Xkdi2Y9aMR/HLmcv7y+hes2X6Qv44fSGx0pNehGWNMyPhNNKp6HYCIjFLVBb7rRGRUFbukqupFIhID5LrH2AbcICJpoQu5cWkRE8WUq4bQ9704HnhnHRuLD/HfqzPo0ibW69CMMSYkgjlH858gl00Vkc9w6qE96LtCVZfWPLSmQ0S4fXRvpl49lA27DjFmyicsLtzrdVjGGBMSfisDuHXLRgI/Bf7ps6o1cImq1stzL/WtMkBNrd1xkBufXcSO/Uf58yUDuCwj2euQjDFNgFeVAWKAVjjda3E+twPAhCqCvFdE2vo7mIicJSIX1S7cxq9Plzjm3jaKYd3b8suXlvP7V1ZRWlbudVjGGHPSAp2j+RD4UET+585HU50VwKsichRYDBTjFNPsDaQB7wJ/qXXETUDbljE8c91w/vz6Gp5esJn1Ow8x5aohxLeI8To0Y4ypsWqLatb4gCK9cUrPVEwTsAb4SFW/DukT+dHQu84qe3HRFu6ds5KE+FgevyaDUzrHeR2SMaYRCmfXWcgTjdcaW6IByC/Yy83P5fP18VL+eUUa5/bv4nVIxphGxuvqzcZjQ1Pb8srto+jZqRU3PZfPf+avt4s7jTENRqDKADe63WCI42m3qOZyEUmvuxANQEKb5rx48wguTkvkgXfWcdv0xRw5Xup1WMYYU61ALZo7gM3u/SuBQUB34E7Ab+2yqi7m9HOBp6mh2OhI/nlFGndf0Jc3V+7g0kc/o2jvEa/DMsaYgAIlmlJVrahjfxHwrKruVtV3gZYB9gv2Ak9zEkSEm7/Xk6euHUbR3iOMnbKAvI3VTg9kjDGeCZRoykUkQURigdE4w5MrNK+8sYiMEJGfAx1F5E6f2+8AK94VYmf26cTLt40ivkU02U/kMS03mBHoxhhT9wIlmt8Ai3C6z+ap6ioAEfkesLGK7Wt0gaepvZ4dW/HybaM4vXcH7n15Jf9vzgqOl9rFncaY+iXg8GYRiQLiVHWvz7IWQKSqHvSzT2qQF3iGRWMc3lydsnLlH2+t5bEPv2R4t3Y8MimdDq2aeR2WMaYB8Wx4s6qWViQZd+TZaJzzLRsC7NZMRKaKyNsi8l7FLYQxm0oiI4S7LujLQxPTWFa0j3FTFrBq236vwzLGGCCI62hEJEtE/g0UAHOBj4C+AXaZCSwB7gV+6XMzYTYurSsv3TKSclUuffRTXl2+zeuQjDEm4HU0fxGR9cCfgeXAEKBYVZ/x7UqrQqmqPqqqC1U1v+IW4riNHwOT2jD3x6Pon9iGH09fwv1vraW83C7uNMZ4J1CL5gZgJ/Ao8Jyq7saZWbM6r4jIre6ItXYVt1AEa4LTKS6W6TdmckVGMlPe38BNz+Vz8GhJ9TsaY0wYBJqPJhI4B+dizdHA+8DZQLKq+r0kXUQ2VbFYVbVH7cOtXlMcDOCPqvLsZwX84dXV9OjQkicmZ5DaPtAlUMaYpsqTwQCqWqaqb6rqZKAn8DKwANgqItMD7Ne9iludJBlzIhFh8shuPHf9cIoPHWPslAV8sv4rr8MyxjQxAc/R+Dw8Q1VnqeoEnPll3gywXwt3ErSp7uPeNuGZt0b26sC8206jS+tYrnkqjyc/2WRFOY0xdSbQOZrzfe7fV3FHVQ+o6rMB9nsaOI4zDTTAVuBPJx2hCYmU9i2YfetIzunXmT++uppfvrScY6VlXodljGkCwjFNQE9V/TtQAqCqRwAJw/OYGmrZLIpHs4dyx+jevJRfxMSpuew6cNTrsIwxjZzfqZyBTiJyJ06SqLj/DVV90M9+x0WkOe4INRHpCRwLRbCm9iIihJ+dcwp9u8Tx85nLGDPlE/57dQZpyfFeh2aMaaQCtWgex6lV1srnvu/Nn9/inMNJFpEcYD7wq5BEa0LmgoEJzPrRSKIjI7j8v58xe3GR1yEZYxqpaqdyFpEOqlqjoUoi0h7IwmkN5dZ0/9qw4c01s+fwcW7NySd34x5uPL07/3d+X6IibeJVY5oaT4Y3i8hFIlIMLBeRIhEZ6W/bSvudAfQHDuJUbu7nLjP1ULuWMTz3w0wmj0jl8Y83cd3/Pmf/Ebu40xgTOoEu2FwOXK6qX4hIJvB3Vf1etQcUecXnYSwwHMhX1bNCEXB1rEVz8p5fWMhv5q6kY6tmTB7ZjQlDk2hvVaCNaRLC2aIJNBigVFW/AFDVPBEJdF7mG6o6xvexiCQD/zrpCE2duXJ4Cqd0bsV9b6zlr298wQNvr+OCgV2YlJVKRmpbRGzwoDGm5oIZdVbl4wCjziorAk4NZkMR+RlOjTUFVgDXAQnADKA9kA9crarHg3xuU0NDU9vx4i0jWLfzINPzCpmVX8Tcpdvo0zmO7KwULh7Sldax0V6HaYxpQAJ1nf020I6q+ns/+/2Hb4tvRgBpwGZVnRQwEJGuwCdAP1X9WkReBF4HfgDMVtUZIvIYsExVH/V3HOs6C60jx0t5Zdk2puUWsmLrfppHRzIuLZFJWakM6NrG6/CMMSESzq6zaked1fiAIpN9HpbiJJkFQezXFcgFBuMMIngZZ5K1HKCLqpaKyAjgd6p6nr/jWKIJn+VF+8jJLWTusq0cLSlncFIbsjNTGZuWSGx0pNfhGWNqoUElmtoQkTtw5r/5GngbuANneHQvd30y8IaqDvB3DEs04bf/6xLmLC4iJ6+Q9bsOMaBra6ZenUFifHOvQzPGnCTPpnKuCRFZISLLq7itcEewVbd/W2Ac0B1IBFpyYr21QPveJCKLRGRRcXFxrV6HqV6b5tFcO6o7b//sDB6bNJSCr44wdsonLNq8x+vQjDH1UKDBADV1GU5L5GSdDWxS1WIAEZkNjALiRSTKnQMnCadI5wlUdSowFZwWTS1iMDUgIpw/oAu9OrXkxmfzufLxXP44bgATh6d4HZoxph6pNtFUrnHm2o9zbcxSn2XTVTVdRJ5T1atPIpZCIEtEWuAkrNHAIpwJ1ybgjDybDMw9iWObMOrVKY6Xbx3F7TOWcNfsFazefoBfX9SPaKswYIwhuBZNhnuruBDzImA5cIuIzHQrNQPEiMhVwEgRGV/5IKo6O9CTuNfqvAQsxhlEsASnlfIaMENE/uQuezKImE0da9MimqevHcZ9b37B1I82sn7nIR7OTqddyxivQzPGeCyYWmcfAT9Q1UPu41Y4//zPx2nV9HOXnwZkA5cD8yodRlX1+hDHXiUbDOC92YuLuGv2CjrFNeOJyRn07dLa65CMMdXwqjJAhU6cWOa/BOjsXuviuzxBVX8kIkvccyamiRqfnkSPjq246dlFjH/kUx68fDDnD0jwOixjjEeC6UTPAfJE5LfuRZwLgOki0hJY7bPd3e7PW0Ico2mA0pLjeeX20zilcxy3TFvMP99ZR3m5jdMwpikK6joaERnGt1MzL1DV7/RNicg7OBUBhgEfV16vqmNrF2pwrOusfjlaUsY9c1Yya3ER5/XvzIOXp9GyWSgHOxpjQsHrrjNwTtBvrdheRFJUtbDSNhcC6cBzwAMhi9A0aLHRkdx/2SD6Jbbmz6+tZvwjn/L4NRmktG/hdWjGmDoSzGCA23FmzdwJlOFMZqaqOsjP9h0rroXxWVZxHUzYWYum/vp4fTE/nr4EEXjkqnRG9urgdUjGGJfXlQHuAPqoan9VHaSqA6tKMiLyCYCqFovIc5VWLwxBrKaBO713R+beNooOrZpx9VMLeebTzdSnEkjGmPAIJtFswblAszotfe73r7TOJjIxAHTr0JI5t47k+3068tt5q7hr1gqOlZZ5HZYxJoyCOUezEfhARF7DZ5hzFfPRBPpqal9bzTfiYqOZenUGD76zjinvb2BD8SEemzSUjnE2m6cxjVEwiabQvcW4N3/iReQSnFZSvE91AAFs4hJzgogI4Rfn9aFvQhy/mLmMsVM+YerVGQxMso+KMY1NyKYJEJGnA61X1etC8kTVsMEADc+qbfu56dl8vjp0jL9PGMS4tK5eh2RMk+PJ8GYR+Zeq/lREXqGKrq/K18XUVSIxjU//xDbM/fEobp22mDtmLGXN9oP88rw+REbYqT1jGoNAXWcVI8fur4tATNPWoVUzpt2Qye9eWcVjH37J2h0HeOjKIbSOjfY6NGNMLdWrGTZDwbrOGr5puQX8bt4qUtq34N8ThzCgq523MSbcPL2ORkQuEpElIrJHRA6IyEERORCOYIwBmJSVyrQbMtl3pISL/vMJ4x9ZwKz8Io6W2DBoYxqiYCoDbADGAys0iOaPO3HZz4EUVb1RRHrjXPD5aigCro61aBqP/UdKeGlxETl5BWwsPkx8i2gmpCdxVWYKPTq28jo8YxqVcLZogkk07wOjVbU8qAOKvADkA9eo6gA38Xyqqmm1DTYYlmgaH1Uld+MepuUV8NbKHZSWKyN7tmdSVirn9OtsM3kaEwJeF9X8FfC6iHxI4As2K/RU1StE5Ep3uyMiYsOHzEkTEUb0bM+Inu3ZdfAoMxcVMT2vkFtzFtMxrhlXZCRzZWYKXeObex2qMaYKwSSaPwOHgFgCX7BZ4biINMcdEi0iPTlx4jRjTlqnuFhu+34vbvleTz5aV8y03AIe/mADj3ywgdGnduaP4wbQpU2s12EaY3wEk2gSVXVADY75W+BNIFlEcoBRwLUnEZsxfkVGCN/v24nv9+1E0d4jvPD5Fp76ZBNjpnzCf68eSnpKW69DNMa4guncfl1Ezg32gKr6Ds7ggWuB54EMVf3gpKIzJghJbVvw83P7MOe2UbSIiWTif3OZuWiL12EZY1zBJJofAW+KyFF3aHMww5tjgb3AAaCfiJxR20CNqc4pneOYe9sohndvxy9fWs7vX1lFaVlQY1iMMWFUbdeZqsbV5IAich9wBbAKqPgrV+CjGkdnTA3Ft4jhf9cN4y+vf8FTCzaxfuchplw1hPgWwZxeNMaEQ1BTOYvIWKCiVfJBNdfEXIxz3YwNADCeiIqM4Ddj+tE3IY5756xk3MMLePyaDE7pXKPvTMaYEAmmMsDfcGbZXO3e7hCRvwbYZSNgBaqM5y7PSOb5m7I4cryMSx5ewNurdngdkjFNUjAXbC4H0iou2BSRSGBJ5emcReQ/OF1kXYHBwHxOvO7mJ6ENvWp2waapbMf+o9z83CKWFe3n5+ecwo/P6oVd2mXMiby+YBMgHtjj3vdX4bDiv3s+MK/SusZVudM0KF3axPLCzSO4e/YKHnhnHWt2HOD+ywbTIibYj78xpjaC+Uv7K7DELUUjOOdq7qq8kao+AyAid6jqQ77rROSOEMRqzEmLjY7kwcsHc2pCHH974ws2Fh/msUlD6dahpdehGdPoBTVNgIgkAMPchwtV1W9nt4gsVtX0SsuWqOqQWkUaJOs6M9X5YO0ubn9+CYePlXJW305kZ6VyRu+ONtGaadK8Lqp5CfCequ53H8cDZ6rqy5W2uxK4CjgN+NhnVRxQrqqjQxe2f5ZoTDC27/+aabkFvPD5Fr46dJykts25cngKl2ck0zGumdfhGVPnvE40SytXXq6qhSIiqUB3nK423661g8ByVS0NScTVsERjauJ4aTlvr95BTm4hn23cTXSkcF7/LkzKSiWzezsbNGCaDK8HA1Q1BPo7+6lqAVAAjKhtUMbUlZioCC4alMhFgxLZsOsQ0/MKeSl/C68u306vTq3IzkxhfHoSbZrbiH1jTlYwLZqngH3Aw+6i24B2qnptWCM7SdaiMbX19fEyXl2+jZy8QpZu2UdsdARjByfyq/P70qGVdauZxsnrrrOWwK+Bs91F7wB/UtXD4QiotizRmFBauXU/OXmFzFpcRMdWzZh6zVD6J/ob4W9Mw+VpomloLNGYcFhRtJ+bnlvE3iPHuf+ywVw0KNHrkIwJqXAmmmBK0JwiIlNF5G0Rea/iFmD7USLyjoisE5GNIrJJRDaGNmxj6tbApDbM/fEo+ie24cfTl3D/W2spL29cX9KMCZdgBgPMBB4DngDKgtj+SeBnOBUCgtnemAahU1ws02/M5Dcvr2LK+xv4YscB/nlFGnGxNlDAmECCSTSlqvpoDY65X1XfONmAjKnPmkVF8rdLB9IvsTV/eHU14x/5lMevybAKA8YEEMzEZ6+IyK0ikiAi7SpuAbZ/X0T+ISIjRCS94haqgI3xmogweWQ3nrt+OMWHjjHu4QV8vL7Y67CMqbeCGXW2qYrFqqo9/Gz/vp/tz6rmefoAL/gs6gH8BnjWXd4N2Axcrqp7/R3HBgOYulS4+wg3PruI9bsO8v9+cCo/PK27XeRpGqQmN+rMnYpgK5CJc93OHlX9m4jcBbRV1f/zt68lGlPXDh0r5c4XlvL26p1cmp7Eny8ZQGx0pNdhGVMjnlQGEJHxgXZU1dmVtp+kqtNE5E4/2z9Yg7hGA1+qaoGIjAPOdJc/A3wA+E00xtS1Vs2ieGzSUP41fz3/nr+eL4sP8d+rh9K5dazXoRlTLwQaDDAmwDoFZldaVnE2NBTz5U4Ennfvd1bV7e79HUDnEBzfmJCKiBDuPOcUTu0Sx50vLuPMf3zA2MGJZGelMCgp3uvwjPFUves6E5EYYBvQX1V3isg+VY33Wb9XVdtW2ucm4CaAlJSUoQUFBXUZsjEn2LDrEE98vJG5S7fxdUkZg5LakJ2ZwpjBiTbZmqm3mtQ5Grer7DZVPdd9vBZnWoLt7rw4H6hqH3/72zkaU18cOFrCnMVbyckrYN3OQ8TFRnFpehLZmSn07hyKhr8xodPUEs0M4C1Vfdp9/A9gt89ggHaq+it/+1uiMfWNqrKoYC85uQW8vmIHx8vKGd69HdmZKZw/oAvNomzggPGeJ4lGRC5T1Zki0l1VqxriHPpgnAKehUAPn4nW2gMvAik40xBcrqp7/B3DEo2pz3YfOsZL+UVMX1hIwe4jtG8Zw2UZyVw1PIWU9i28Ds80YV4lmsWqml7V1MwBDyhyB/A0zoRnTwBDgLtU9e1QBFwdSzSmISgvVz7Z8BU5eQW8u2YX5aqc0bsj2ZkpnNW3E1GRwVxLbUzoeJVo3sEZXTaME6dmBkBVx/rZb5mqDhaR84CbcaYYeK4myao2LNGYhmbH/qPM+LyQGQu3sOPAURLaxDJxWApXDEumSxsbIm3qhleJJgZIB54Dbqi8XlU/9LPfclUdJCIP4Zy4n1PV1M/hYonGNFSlZeXM/2IXOXmFfLSumMgI4exTOzEpK5VRPTsQEWEVB0z4eD3xWUdVLRaRVgCqeqia7Z8GugLdgcFAJE7CGRqakAOzRGMag4Ldh5m+sJCZi4rYc/g43dq34KrMFCYMTaZdyxivwzONkNeJZgBOq6YdIEAxMFlVV/rZPgJIAzaq6j73ZH5XVV0eysD9sURjGpNjpWW8uXIHObmFLNy8h5jICCZkJPHrC/vRPMZGq5nQ8aQEjY+pwJ2q+r4bzJnuspFVbayq5SKyE+gnInZ1mjG10CwqknFpXRmX1pV1Ow/y7GebyckrZHnRPqZenUFifHOvQzSmWsG0aJap6uDqlvmsuw+4AljNtxOfqb/BA6FmLRrT2M1fs5M7ZiwlNjqCxyYNJaNboFk7jAmOp1M5AxtF5Nci0s293QsEmpr5YqCPqv5AVce4tzpJMsY0BaNP7czLt42kVbMornw8lxkLC70OyZiAgkk01wMdcYpozgI6uMv82QjY3LbGhFGvTnHMve00snq0567ZK/jN3JWUlJV7HZYxVar2HIo7ydhPanDMI8BSEZkPHPM5Tk2OYYypRpsW0Tx97TDue/MLHv94E+t3HuLh7HQblWbqnXCcrJ/n3owxYRYVGcE9F/ajb5fW3D1nBWOnfMITkzPo26W116EZ8416V1SztmwwgGmqlm7Zx03PLuLQsVIevHww5w9I8Dok04B4OhhAREYFuexF9+cKEVle+RaacI0x/qQlx/PK7adxSuc4bpm2mH++s47y8sb1RdI0TMEMb/5OUU0/yxLcOWNSqzqOqtbJbGTWojFN3dGSMu6Zs5JZi4s4r39n/jp+kJ23MdXy5IJNERmBc1FmRxG502dVa5yyMieomG65rhKKMaZqsdGR3H/ZIPoltubPr63m/bXzuXBgApOyUkhPaYuI1UwzdSvQYIAYoJW7je90gAeACf52EpEs4D/Aqe4xIoHDqmpnJ42pIyLCD0/rzum9OzAtt4DZi7cyZ8lW+naJIzsrlYvTEomLtasQTN0IpusstSatFBFZBEwEZgIZwDXAKap6d20CDZZ1nRnzXYePlTJv2Tam5RawatsBWsQ4pW0mZaXQP7GN1+GZesDropqnAL8AuuHTAlLVs/xsv0hVMyqmC3CX2TQBxtQDqsqyov3k5BYwb9k2jpWWk5YcT3ZmCmMGJxIbbYU6myqvE80y4DEgn29rl6Gq+X62/wg4G2d2zR3AduBaf7XRQs0SjTHB2X+khFmLi8jJK+DL4sO0jo1iwtBkrspMoVenVl6HZ+qY14kmvyZzybijznbhlKH5GdAGeERVN9Qm0GBZojGmZlSV3I17mJZXwNurdlBSpozo0Z7srBTO7deFmCibVrop8DrR/A4ncczhxJIye8IRUG1ZojHm5BUfPMaLi7YwPa+Qrfu+pkOrZlwxLImJw1JIbtfC6/BMGHmdaDZVsVhVtYef7S8C/gik4pzTEXf7Ohl1ZonGmNorK1c+Wl9MTm4B732xCwW+36cT2ZkpnNmnE5E2rXSj42miqfEBRTYA44EV6kF9G0s0xoTW1n1fM2NhITM+30LxwWN0jW/OlcOTuXxYMp3iYr0Oz4SI1y2aFsCdQIqq3iQivXHmm3nVz/bvA6NV1ZOa5ZZojAmPkrJy3lm9k5y8AhZs2E1UhHBu/85MykxlRM/2diFoA+f1VM5P44w4q5i6eSvONTJVJhrgV8DrIvIhJ57TebAWcRpjPBYdGcEPBibwg4EJbCw+xPS8QmbmF/H6ih306NCSqzJTmDA0ifgWVu7GnCiYFk3FdTHfXAtTzVTObwOHgBXAN60aVf196ML2z1o0xtSdoyVlvL5iO9NyC1hcuI9mURFcNCiR7KwUhiTHWyunAfG6RXNcRJoD6gbTE5+WShUSVXVAKIIzxtRvsdGRjE9PYnx6Equ3HSAnr4CXl2xl1uIiTk1ozaSsFMaldaVVs3BMfWUaimBaNOcA9wL9gLeBUTgXYH7gZ/u/A++q6tuhDTU41qIxxluHjpUyd+lWpuUWsmb7AVrGRHLxkK5kZ6bSL9FKHtZXno86E5H2QBbOUOVcVf0qwLYHgZbAcfdmw5uNaYJUlSVb9pGTW8iry51yN+kp8WRnpnLhoAQrd1PPeD3q7BLgPVXd7z6OB85U1ZfDEVBtWaIxpv7Zd+Q4L+UXMT2vkI1fHSa+RTQT0pO4KjOFHh2t3E194HWiWaqqaZWW+S2SKc7Zv2ygu6r+UUSSgQRVXRiimAOyRGNM/aWqfPblbnLyCnlr1Q5Ky5VRvdqTnZnKOf06Ex1p5W684vVggKp+84H2ewRntNlZOBUCDgEPA8NqHJ0xplEREUb26sDIXh3YdfAoMxc5rZxbcxbTMa4ZE4clM3F4Cl3jm3sdqgmhYFo0TwH7cJIFwG1AO1W91s/2i1U1Pdjh0KFmLRpjGpaycuXDdbvIyS3kvbW7EOCsvp3IzkzljFM6WrmbOuJ1i+Z24NfACzhDnN/BSTb+lIhIJN8Oh+6Iz/U0xhjjKzJCOKtvZ87q25mivUeYsXALMz7fwrtrPiepbXOuHJ7C5RnJdIxr5nWo5iQFbNG4CeNdVf1+0AcUyQauANKBZ3Cmfb5XVWfWMtagWIvGmIbveOm35W4+/XI30ZHCef27kJ2ZSlaPdnYhaBh41qJR1TIRKReRNhWjzgIRkQhgE04ZmtE4Q5svVtU1IYnWGNMkxERFcOGgBC4clMCXbrmbl/KLeHX5dnp2bEl2ZiqXpifRpkW016GaIARzjmYuMASny+xwxXJV/Ymf7ets2uaqWIvGmMbpaEkZry7fTk5eAUsK9xEbHcGYQYlkZ6UyOKmNtXJqyevhzZOrWq6qz/jZ/n7gM2C2TRNgjAmHVdv2k5NXyMtLtnLkeBn9E1uTnZnKuLREWlq5m5NSHyoDNMeZJmBtENtWVAYoBY5ilQGMMWFy8GgJc5duY1puAV/sOEirZlFcMqQr2Vkp9O1i5W5qwusWzRjgfiBGVbuLSBrwB1UdG46AassSjTFNj6qyuHAfObkFvLpiO8dLy8lIbUt2VgoXDLByN8HwOtHk41x8+YHPdTErw1Gh2S1v8wQwAGd49PXAWpyh1d2AzcDlqrrX3zEs0RjTtO09fJxZi4vIyStk01eHadsimssykrlqeArdOrT0Orx6K5yJJph6DyVVjDgL13UxDwFvqmpfYDCwBrgLmK+qvYH57mNjjKlS25Yx3HB6D+bf+T1ybshkRM/2PPXJJs68/wOufjKPN1dup6TMLu2rS8G0aJ7k23/wlwI/AaJV9ZaQBiLSBlgK9PAdRCAia3GKeG4XkQScllUff8exFo0xprJdB47ywudbeH5hIdv2H6Vz62ZcMSyFicOSSbRyN4D3XWctgHuAc91FbwF/UtWjIQ3EOfczFViN05rJB+4AtqpqvLuNAHsrHvvsexNwE0BKSsrQgoKCUIZmjGkkysqV97/YRU5eAR+sK0aA0ad2JjszhTN6dySiCZe78STRiEgscAvQC2da5idVtTQcQbjPlwHkAqNUNU9EHgIOALf7JhYR2auqbf0dx1o0xphgbNlzhOcXFvLioi18deg4ye2ac9XwVC7LSKJDq6ZX7sarRPMCUAJ8DFwAbFbVn4YjCPf5uuBMqtbNfXw6TnddL6zrzBgTJsdLy3lr1Q5y8grI3biH6Ejh/AEJZGemkNm96ZS78aoETT9VHegG8CQQ1vlkVHWHiGwRkT7u9TqjcbrRVgOTgb+5P+eGMw5jTNMSExXBmMGJjBmcyIZdB8nJK2RWfhGvLNtGr06tyM5MYXx6Em2aW7mbkxWoRbNYVdP9PQ5LMM55mieAGGAjcB3OyLgXgRSgAGd48x5/x7AWjTGmtr4+XsYry7eRk1fIsi3flruZlJXK4OR4r8MLC6+6zsr4traZAM2BI9Txlf41ZYnGGBNKK7fuJyevgJeXbOPrkjIGdm1DdmYKY9MSaRHTeMrdeF6CpiGxRGOMCYcDR0uYu2Qr03ILWbvzIHHNohif3pWrMlPp0yXO6/BqzRJNDViiMcaEk6qSX7CXabkFvL5iB8fLyhnWrS2TslI5f0AXmkU1zHI3lmhqwBKNMaau7Dl8nJfyt5CTV0jB7iO0axnDZUOTuCozhdT2DavcjSWaGrBEY4ypa+XlyoIvv2JabgHvrtlFWblyeu8OZGemcvapnYiKDKbal7cs0dSAJRpjjJd27P+23M2OA065m4nDUpg4PJmENvW33I0lmhqwRGOMqQ9Ky8p574tdTMsr5KN1xURGCKP7diI7K5XTe3Wod+VuvLpg0xhjzEmKiozg3P5dOLd/Fwp3H2G6W+7m7dU7SWnXgqsyU7hsaBLtm0C5G2vRGGNMHTlWWsabK3eQk1fIwk17iImM4IKBXcjOTGVYt7aelruxrrMasERjjGkI1u/8ttzNwWOlnNK5FdmZqVyS3pXWsXVf7sYSTQ1YojHGNCRHjpfyyjKn3M3yov00j45kXFoi2ZmpDExqU2dxWKKpAUs0xpiGannRPnJyC5m7bCtHS8oZnNSG7MxUxgxOpHlMeC8EtURTA5ZojDEN3f6vS5izuIicvELW7zpEXGwUl6YnkZ2ZQu/O4Sl3Y4mmBizRGGMaC1Xl881OuZs3Vm6npEwZ3r0dk7JSOa9/55CWu7FEUwOWaIwxjdFXh44xc1ER0xcWsGXP17RvGcNlGclkZ6aQ3K5FrY9viaYGLNEYYxqz8nLl4w1OuZv5a3aiwBm9O5KdmcJZfU++3I0lmhqwRGOMaSq27/+aGQu3MOPzQnYeOEZCm9hvyt10bh1bo2NZoqkBSzTGmKampKyc+Wt2kZNXwMfrvyIyQjjn1M5kZ6Uwqmdw5W6sBI0xxhi/oiMjOH9AF84f0IXNXx3mebfczZurdtCtfUW5m2TatozxJD5r0RhjTCN0tKSi3E0Bn2/eS0xUBBcOTCA7M4Whqd8td2NdZzVgicYYY060dsdBcvIKmL14K4eOldK3SxzZmSlcPKQrcW65G0s0NWCJxhhjqnb4mFPuZlpeASu3HqBFTCTj0rqSnZnCwKR4O0djjDGmdlo2i2Li8BQmDk9h2ZZ95OQVMGdJEc8vLAzr81qiMcaYJmhwcjyDk+O558J+zF5cxPX3he+56v9E1sYYY8KmTfNorhvVPazPYYnGGGNMWFmiMcYYE1aWaIwxxoSVJRpjjDFhZYnGGGNMWFmiMcYYE1aWaIwxxoSVJRpjjDFh1ehqnYlIMVDgdRwh1gH4yusg6il7b6pm74t/9t5UrY+qxoXjwI2uBI2qdvQ6hlATkUXhKnbX0Nl7UzV7X/yz96ZqIhK2asTWdWaMMSasLNEYY4wJK0s0DcNUrwOox+y9qZq9L/7Ze1O1sL0vjW4wgDHGmPrFWjTGGGPCyhKNMcaYsLJE4wERSRaR90VktYisEpE73OXtROQdEVnv/mzrLhcR+beIbBCR5SKS7nOsye7260VkslevKdREJFJElojIq+7j7iKS574HL4hIjLu8mft4g7u+m88x7naXrxWR8zx6KSEjIvEi8pKIfCEia0RkhH1mHCLyM/dvaaWIPC8isU31MyMiT4nILhFZ6bMsZJ8TERkqIivcff4tIlJtUKpqtzq+AQlAuns/DlgH9AP+DtzlLr8LuM+9/wPgDUCALCDPXd4O2Oj+bOveb+v16wvRe3QnMB141X38IjDRvf8Y8CP3/q3AY+79icAL7v1+wDKgGdAd+BKI9Pp11fI9eQa4wb0fA8TbZ0YBugKbgOY+n5Vrm+pnBjgDSAdW+iwL2ecEWOhuK+6+F1Qbk9dvit0UYC5wDrAWSHCXJQBr3fv/Ba702X6tu/5K4L8+y0/YrqHegCRgPnAW8Kr7gf4KiHLXjwDecu+/BYxw70e52wlwN3C3zzG/2a4h3oA27j9TqbS8yX9m3ESzxf2nGOV+Zs5ryp8ZoFulRBOSz4m77guf5Sds5+9mXWcec5vtQ4A8oLOqbndX7QA6u/cr/pAqFLnL/C1v6P4F/Aoodx+3B/apaqn72Pd1fvMeuOv3u9s3tvemO1AMPO12KT4hIi2xzwyquhW4HygEtuN8BvKxz4yvUH1Ourr3Ky8PyBKNh0SkFTAL+KmqHvBdp87XhSY39lxELgJ2qWq+17HUM1E43SGPquoQ4DBOF8g3mvBnpi0wDicZJwItgfM9Daoe8+JzYonGIyISjZNkclR1trt4p4gkuOsTgF3u8q1Ass/uSe4yf8sbslHAWBHZDMzA6T57CIgXkYrafL6v85v3wF3fBthN43tvioAiVc1zH7+Ek3jsMwNnA5tUtVhVS4DZOJ+jpv6Z8RWqz8lW937l5QFZovGAO0rjSWCNqj7os2oeUDG6YzLOuZuK5de4I0SygP1uM/gt4FwRaet+qzvXXdZgqerdqpqkqt1wTtS+p6rZwPvABHezyu9NxXs2wd1e3eUT3RFG3YHeOCcxGyRV3QFsEZE+7qLRwGrsMwNOl1mWiLRw/7Yq3psm/ZmpJCSfE3fdARHJct/ra3yO5Z/XJ62a4g04DafpuhxY6t5+gNNPPB9YD7wLtHO3F+BhnFEwK4AMn2NdD2xwb9d5/dpC/D6dybejznrg/NFvAGYCzdzlse7jDe76Hj773+O+Z2sJYmRMfb8BacAi93PzMs5oIPvMOK/p98AXwErgOZyRY03yMwM8j3OuqgSnJfzDUH5OgAz3ff4SmEKlASpV3awEjTHGmLCyrjNjjDFhZYnGGGNMWFmiMcYYE1aWaIwxxoSVJRpjjDFhZYmmARIRFZFpPo+jRKRYvq10PFZE7vJ/hKCf58yKY/pbHqrnaixE5FoRmVLDfTaLSIdaPu//RGSCe/8JEelX29hEJMHn93ym+7m7wWd9mrvsF5Vj8Nkmwq3wu9Kt+Pu5e40KInK9u2y5u36cu/x+ETkrQFz/EpEzgn0dQbzOiSJyj7/3R0Terah2bE6OJZqG6TAwQESau4/PwefqXFWdp6p/q4tA6uK5RCQynMdvbFT1BlVdHYJD3Qk87vN4JXC5z+MrcaodB3IFTlmYQao6ELgE2CciSTjXrJymqoNwqgEvd/f5D5XK61QQkfZAlqp+VMPXEsgFwJsB1j+HU/HZnCRLNA3X68CF7v0rcS7SAk785ioic0XkGvf+zSKS494/V0Q+E5HFIjLTrbuGiJwvznwni4Hx1QVR6bm6u8dcISJ/EpFD7vITWkYiMkVErnXvjxanSOQKcebRaOYu3ywi97lxXFbpOS9zvwEvE5GPKsfhPn7Vfd5I95t2xTfqn7nrPxCRh0RkqbtuuLu8pRvHQjeucT7Hny0ib4ozP8fffZ7rOhFZJyILcUqfVCzvKCKz3G/xn4vIKHd5exF5W5z5U57AuWiu8vvqL+40Ecl1WwFzqvqm7b62jGpi+857WIVLOfEfcAEQKyKdRURw6om94WffCgnAdlUtB1DVIlXdC3QCDgKH3OWHVHWTe78AaC8iXaqLyf2c/NX9PS4SkXQReUtEvhSRW9xtIkTkEfdz/Y6IvC7ftv4E50LYxQFewzycvzFzkizRNFwzcMplxAKDcKo/V+Um4Dcicjrwc+B2cbpp7gXOVtV0nKvN73SP9TgwBhgKVPWHHshDOEUfB+JcmRyQ+3z/A65w94kCfuSzyW5VTVfVGZV2/Q1wnqoOBsZW8zRpQFdVHeA+x9M+61qoahrOt9Wn3GX34JQkGQ58H/iHOFWSK451BTAQuEKcCewScK5KH4VT8cG3y+oh4J+qOgznH+QT7vLfAp+oan9gDpBSg7ifBf7PbQWscI9VpWpiC/geitO9tVdVj1Va9RJO4h+J88+58vrKXgTGuIngAREZ4i5fBuwENonI0yIyptJ+i/FJjD5G4VRm9lXo/h4/xvk8TcBpIf3eXT8ep2x+P+BqnCkDKgwBlmmAK9fdxNhMnNaUOQmWaBooVV2O88dzJU7rxt92O3H+qbwP/FxV9+D8EfYDFojIUpzaR6lAX5zihOvdP7xpVR/Vr1F827J6Lojt+7jPt859/AzOpE0VXvCz3wLgfyJyI1Bdt9pGoIeI/EdEzgd8q2Q/D+B2w7QWkXicmk53ue/LBzjlSioSwXxV3a+qR3FqaaUCmcAH6hR0PF4p5rOBKe6x5rnP0cp9jdPc534N2BtM3CLSBohX1Q/dbSq/X5UFiq269zABZ1qCyl7ESTQntKL9UdUinN/z3TjTPswXkdGqWobTIpqAM/HfP0Xkdz677sLpcgsmrnnuzxU4E3cdVNVi4Jj7Oz0NmKmq5erUjHvfZ99gWmWB4jFBsETTsM3DmYejuj/4gTjVaSv+UAR4R1XT3Fs/Vf1hiGKq6pthKSd+1mKDPNbhKp9A9RacFlkykO9+06zyOdxvo4NxksYtfNuqqCpWxXlvLvV5b1JUdY273vfbexlOCyyQCJzzCRXH6qqqh6rZp+I1Boq71vy8h76+porfk/uPugTnvOD8IJ/rmKq+oaq/BP4CXOwuV1VdqKp/xSmgeqnPbrFuDJVVFVfF76WcE39H5VT/OzoXeDuIl+EvHhMESzQN21PA71V1hb8N3HMPF+B0EfzC7RLJBUaJSC93m5YicgpOUcJuItLT3b2m/dILcP5hAGT7LC8A+olTFTcep7ouOIULu1XEgdOt8SHVEJGeqpqnqr/B+XabDGwG0tz++GSg4pxLByBCVWfh/GNN9znUFe42p+FUrd2PU7X2drfvHp+uHn/ygO+5512iOfF80tvA7T5xp7l3PwKucpddgFMcs/Jr/E7cbnx73W5QqP798hubn/fQ1zqcFnNVfoPTfVcW4LkrniddRBLd+xE43bwFIpIoPvPT43QVFvg8PgVn8EFla4BeVSwPZAFwqfvZ6IxTrBW3hRilqrureQ2C0428uYbPa1zVZXtTj7ndEv/2t16cE+uP41Re3SYiP8dJTmfhzKn+vLsNwL2quk5EbgJeE5EjOH3ecTUI6Q5guoj8Hz6lw1V1i4i8iPOPYxOwxF1+VESuA2aKMy/I5zhzu1fnHyLSG6f1MZ9vRz5twunSWsO3J3e74sxKWfGl6m6f4xwVkSVANE6lWoA/4szwudzdZxNwkb9AVHW72+XzGbAPpxJ3hZ8AD4vIcpy/tY9wWie/x3nvVwGf4pS5r8xf3JOBx0SkBU732nUnGZu/97Bi38PuCfVeqrqh0rpP/T0n8F8R+Zd7f4v7Wh/3+ZwtxKn42xm4301CR3GSXcXJ+2icZLKoiuO/BtxMzVp4s/h26oAtOJ+N/TitsncrbXutiFzs8zgLJ8nk6rezdZoasurNJmxE5JCqtvI6jqqIyAfAL1S1qn9mBhCRS4ChqnqvB8+brqq/9rP+E+AiVd1Xg2O2UtVDbhdhxQi8PwFPqGpuNfs+BMxT1aC6Cs13WYvGGFMlVZ3j0UirKOCBAOt/jjNAY18Njvmq220bA/zRPdd0Q+BdvrHSkkztWIvGGGNMWNlgAGOMMWFlicYYY0xYWaIxxhgTVpZojDHGhJUlGmOMMWH1/wFhcSduvPcSWQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "MLSS = np.linspace(500,10000,30)\n",
    "PFAS_inEffluent = []\n",
    "IncomingWaste.set_flow('LFLeachate', 1000)\n",
    "IncomingWaste.calc()\n",
    "for i in MLSS:\n",
    "    WWT.InputData.BioTrtmnt['sol_cont']['amount'] = i\n",
    "    WWT.calc(IncomingWaste.Inc_flow)\n",
    "    PFAS_inEffluent.append(round(WWT.report(normalized=True)['WWT Effluent']['PFOA']))\n",
    "plt.plot(MLSS,PFAS_inEffluent)\n",
    "plt.xlim(500,10000)\n",
    "plt.xlabel('Mixed liquor suspended solids (MLSS) (mg/L)')\n",
    "plt.ylabel('Percent of Incoming PFAS that \\n remains in the Effluent (%)')"
   ]
  }
 ],
 "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
}
