{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Advanced 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",
    "\n",
    "\n",
    "<img src=\"../Images/ProcessModels/AdvWWT.png\" alt=\"Drawing\" style=\"width: 700px;\"/>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for Advanced 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>Dictonary_Name</th>\n",
       "      <th>Parameter Name</th>\n",
       "      <th>Parameter Description</th>\n",
       "      <th>amount</th>\n",
       "      <th>unit</th>\n",
       "      <th>Reference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Treatment option</td>\n",
       "      <td>Trtmnt_optn</td>\n",
       "      <td>frac_to_RO</td>\n",
       "      <td>Fraction of incoming PFAS-containing water to RO</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Treatment option</td>\n",
       "      <td>Trtmnt_optn</td>\n",
       "      <td>frac_to_GAC</td>\n",
       "      <td>Fraction of  incoming PFAS-containing water to GAC</td>\n",
       "      <td>1.000000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Treatment option</td>\n",
       "      <td>Trtmnt_optn</td>\n",
       "      <td>frac_to_IonEx</td>\n",
       "      <td>Fraction of  incoming PFAS-containing water to Ion Exchange</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>RO paramteres</td>\n",
       "      <td>RO</td>\n",
       "      <td>ts_rem_med</td>\n",
       "      <td>Solids content of RO concentrate wet</td>\n",
       "      <td>0.017663</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>RO paramteres</td>\n",
       "      <td>RO</td>\n",
       "      <td>VS_rem_med</td>\n",
       "      <td>VS content of RO concentrate - dry</td>\n",
       "      <td>0.700000</td>\n",
       "      <td>kg VS/kg TS</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>RO paramteres</td>\n",
       "      <td>RO</td>\n",
       "      <td>C_rem_med</td>\n",
       "      <td>C of RO concentrate - dry</td>\n",
       "      <td>0.001207</td>\n",
       "      <td>fraction TS</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>RO paramteres</td>\n",
       "      <td>RO</td>\n",
       "      <td>frac_effl_rem_med</td>\n",
       "      <td>RO concentrate fraction relative to total water treated</td>\n",
       "      <td>0.150000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>RO Removal Efficiency</td>\n",
       "      <td>RO_RemEff</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>GAC paramteres</td>\n",
       "      <td>GAC</td>\n",
       "      <td>ts_rem_med</td>\n",
       "      <td>Solids content of removal media - wet</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>GAC paramteres</td>\n",
       "      <td>GAC</td>\n",
       "      <td>VS_rem_med</td>\n",
       "      <td>VS content of removal media  - dry</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>GAC paramteres</td>\n",
       "      <td>GAC</td>\n",
       "      <td>C_rem_med</td>\n",
       "      <td>C content of removal media  - dry</td>\n",
       "      <td>0.700000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>GAC paramteres</td>\n",
       "      <td>GAC</td>\n",
       "      <td>GAC_dens</td>\n",
       "      <td>GAC bed density</td>\n",
       "      <td>435.000000</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>GAC paramteres</td>\n",
       "      <td>GAC</td>\n",
       "      <td>Bed_vol_ratio</td>\n",
       "      <td>Bed volumes: V(H2O)/V(GAC)</td>\n",
       "      <td>5000.000000</td>\n",
       "      <td>ratio</td>\n",
       "      <td>Woodard et al. (2017)</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>GAC Removal Efficiency</td>\n",
       "      <td>GAC_RemEff</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS removal efficiency</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>Ion Exchange paramteres</td>\n",
       "      <td>IonEx</td>\n",
       "      <td>ts_rem_med</td>\n",
       "      <td>Solids content of removal media - wet</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>Ion Exchange paramteres</td>\n",
       "      <td>IonEx</td>\n",
       "      <td>VS_rem_med</td>\n",
       "      <td>VS content of removal media  - dry</td>\n",
       "      <td>0.900000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>Ion Exchange paramteres</td>\n",
       "      <td>IonEx</td>\n",
       "      <td>C_rem_med</td>\n",
       "      <td>C content of removal media  - dry</td>\n",
       "      <td>0.500000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>Ion Exchange paramteres</td>\n",
       "      <td>IonEx</td>\n",
       "      <td>IonExRes_dens</td>\n",
       "      <td>Ion Exchange Resin bed density</td>\n",
       "      <td>500.000000</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td>Ion Exchange paramteres</td>\n",
       "      <td>IonEx</td>\n",
       "      <td>Bed_vol_ratio</td>\n",
       "      <td>Bed volumes: V(H2O)/V(Ion Exchange Resin)</td>\n",
       "      <td>20000.000000</td>\n",
       "      <td>ratio</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>38</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>40</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>44</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>45</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>46</th>\n",
       "      <td>Ion Exchange Removal Efficiency</td>\n",
       "      <td>IonEx_RemEff</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS removal efficiency</td>\n",
       "      <td>0.990000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                           Category Dictonary_Name     Parameter Name  \\\n",
       "0   Treatment option                 Trtmnt_optn    frac_to_RO          \n",
       "1   Treatment option                 Trtmnt_optn    frac_to_GAC         \n",
       "2   Treatment option                 Trtmnt_optn    frac_to_IonEx       \n",
       "3   RO paramteres                    RO             ts_rem_med          \n",
       "4   RO paramteres                    RO             VS_rem_med          \n",
       "5   RO paramteres                    RO             C_rem_med           \n",
       "6   RO paramteres                    RO             frac_effl_rem_med   \n",
       "7   RO Removal Efficiency            RO_RemEff      PFOA                \n",
       "8   RO Removal Efficiency            RO_RemEff      PFOS                \n",
       "9   RO Removal Efficiency            RO_RemEff      PFBA                \n",
       "10  RO Removal Efficiency            RO_RemEff      PFPeA               \n",
       "11  RO Removal Efficiency            RO_RemEff      PFHxA               \n",
       "12  RO Removal Efficiency            RO_RemEff      PFHpA               \n",
       "13  RO Removal Efficiency            RO_RemEff      PFNA                \n",
       "14  RO Removal Efficiency            RO_RemEff      PFDA                \n",
       "15  RO Removal Efficiency            RO_RemEff      PFBS                \n",
       "16  RO Removal Efficiency            RO_RemEff      PFHxS               \n",
       "17  GAC paramteres                   GAC            ts_rem_med          \n",
       "18  GAC paramteres                   GAC            VS_rem_med          \n",
       "19  GAC paramteres                   GAC            C_rem_med           \n",
       "20  GAC paramteres                   GAC            GAC_dens            \n",
       "21  GAC paramteres                   GAC            Bed_vol_ratio       \n",
       "22  GAC Removal Efficiency           GAC_RemEff     PFOA                \n",
       "23  GAC Removal Efficiency           GAC_RemEff     PFOS                \n",
       "24  GAC Removal Efficiency           GAC_RemEff     PFBA                \n",
       "25  GAC Removal Efficiency           GAC_RemEff     PFPeA               \n",
       "26  GAC Removal Efficiency           GAC_RemEff     PFHxA               \n",
       "27  GAC Removal Efficiency           GAC_RemEff     PFHpA               \n",
       "28  GAC Removal Efficiency           GAC_RemEff     PFNA                \n",
       "29  GAC Removal Efficiency           GAC_RemEff     PFDA                \n",
       "30  GAC Removal Efficiency           GAC_RemEff     PFBS                \n",
       "31  GAC Removal Efficiency           GAC_RemEff     PFHxS               \n",
       "32  Ion Exchange paramteres          IonEx          ts_rem_med          \n",
       "33  Ion Exchange paramteres          IonEx          VS_rem_med          \n",
       "34  Ion Exchange paramteres          IonEx          C_rem_med           \n",
       "35  Ion Exchange paramteres          IonEx          IonExRes_dens       \n",
       "36  Ion Exchange paramteres          IonEx          Bed_vol_ratio       \n",
       "37  Ion Exchange Removal Efficiency  IonEx_RemEff   PFOA                \n",
       "38  Ion Exchange Removal Efficiency  IonEx_RemEff   PFOS                \n",
       "39  Ion Exchange Removal Efficiency  IonEx_RemEff   PFBA                \n",
       "40  Ion Exchange Removal Efficiency  IonEx_RemEff   PFPeA               \n",
       "41  Ion Exchange Removal Efficiency  IonEx_RemEff   PFHxA               \n",
       "42  Ion Exchange Removal Efficiency  IonEx_RemEff   PFHpA               \n",
       "43  Ion Exchange Removal Efficiency  IonEx_RemEff   PFNA                \n",
       "44  Ion Exchange Removal Efficiency  IonEx_RemEff   PFDA                \n",
       "45  Ion Exchange Removal Efficiency  IonEx_RemEff   PFBS                \n",
       "46  Ion Exchange Removal Efficiency  IonEx_RemEff   PFHxS               \n",
       "\n",
       "                                          Parameter Description        amount  \\\n",
       "0   Fraction of incoming PFAS-containing water to RO             0.000000       \n",
       "1   Fraction of  incoming PFAS-containing water to GAC           1.000000       \n",
       "2   Fraction of  incoming PFAS-containing water to Ion Exchange  0.000000       \n",
       "3   Solids content of RO concentrate wet                         0.017663       \n",
       "4   VS content of RO concentrate - dry                           0.700000       \n",
       "5   C of RO concentrate - dry                                    0.001207       \n",
       "6   RO concentrate fraction relative to total water treated      0.150000       \n",
       "7   PFOA removal efficiency                                      0.990000       \n",
       "8   PFOS removal efficiency                                      0.990000       \n",
       "9   PFBA removal efficiency                                      0.990000       \n",
       "10  PFPeA removal efficiency                                     0.990000       \n",
       "11  PFHxA removal efficiency                                     0.990000       \n",
       "12  PFHpA removal efficiency                                     0.990000       \n",
       "13  PFNA removal efficiency                                      0.990000       \n",
       "14  PFDA removal efficiency                                      0.990000       \n",
       "15  PFBS removal efficiency                                      0.990000       \n",
       "16  PFHxS removal efficiency                                     0.990000       \n",
       "17  Solids content of removal media - wet                        0.900000       \n",
       "18  VS content of removal media  - dry                           0.900000       \n",
       "19  C content of removal media  - dry                            0.700000       \n",
       "20  GAC bed density                                              435.000000     \n",
       "21  Bed volumes: V(H2O)/V(GAC)                                   5000.000000    \n",
       "22  PFOA removal efficiency                                      0.900000       \n",
       "23  PFOS removal efficiency                                      0.900000       \n",
       "24  PFBA removal efficiency                                      0.900000       \n",
       "25  PFPeA removal efficiency                                     0.900000       \n",
       "26  PFHxA removal efficiency                                     0.900000       \n",
       "27  PFHpA removal efficiency                                     0.900000       \n",
       "28  PFNA removal efficiency                                      0.900000       \n",
       "29  PFDA removal efficiency                                      0.900000       \n",
       "30  PFBS removal efficiency                                      0.900000       \n",
       "31  PFHxS removal efficiency                                     0.900000       \n",
       "32  Solids content of removal media - wet                        0.900000       \n",
       "33  VS content of removal media  - dry                           0.900000       \n",
       "34  C content of removal media  - dry                            0.500000       \n",
       "35  Ion Exchange Resin bed density                               500.000000     \n",
       "36  Bed volumes: V(H2O)/V(Ion Exchange Resin)                    20000.000000   \n",
       "37  PFOA removal efficiency                                      0.990000       \n",
       "38  PFOS removal efficiency                                      0.990000       \n",
       "39  PFBA removal efficiency                                      0.990000       \n",
       "40  PFPeA removal efficiency                                     0.990000       \n",
       "41  PFHxA removal efficiency                                     0.990000       \n",
       "42  PFHpA removal efficiency                                     0.990000       \n",
       "43  PFNA removal efficiency                                      0.990000       \n",
       "44  PFDA removal efficiency                                      0.990000       \n",
       "45  PFBS removal efficiency                                      0.990000       \n",
       "46  PFHxS removal efficiency                                     0.990000       \n",
       "\n",
       "           unit              Reference  \n",
       "0   fraction                            \n",
       "1   fraction                            \n",
       "2   fraction                            \n",
       "3   fraction                            \n",
       "4   kg VS/kg TS                         \n",
       "5   fraction TS                         \n",
       "6   fraction                            \n",
       "7   fraction                            \n",
       "8   fraction                            \n",
       "9   fraction                            \n",
       "10  fraction                            \n",
       "11  fraction                            \n",
       "12  fraction                            \n",
       "13  fraction                            \n",
       "14  fraction                            \n",
       "15  fraction                            \n",
       "16  fraction                            \n",
       "17  fraction                            \n",
       "18  fraction                            \n",
       "19  fraction                            \n",
       "20  kg/m3                               \n",
       "21  ratio        Woodard et al. (2017)  \n",
       "22  fraction                            \n",
       "23  fraction                            \n",
       "24  fraction                            \n",
       "25  fraction                            \n",
       "26  fraction                            \n",
       "27  fraction                            \n",
       "28  fraction                            \n",
       "29  fraction                            \n",
       "30  fraction                            \n",
       "31  fraction                            \n",
       "32  fraction                            \n",
       "33  fraction                            \n",
       "34  fraction                            \n",
       "35  kg/m3                               \n",
       "36  ratio                               \n",
       "37  fraction                            \n",
       "38  fraction                            \n",
       "39  fraction                            \n",
       "40  fraction                            \n",
       "41  fraction                            \n",
       "42  fraction                            \n",
       "43  fraction                            \n",
       "44  fraction                            \n",
       "45  fraction                            \n",
       "46  fraction                            "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "AdvWWT = pspd.AdvWWT()\n",
    "AdvWWT.InputData.Data[['Category','Dictonary_Name','Parameter Name', 'Parameter Description', 'amount', 'unit','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming Waste Water 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>5.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>995.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>2.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>100000.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   5.0     \n",
       "2   Moisture flow  kg   995.0   \n",
       "3   Volume flow    L    1000.0  \n",
       "4   Carbon flow    kg   2.5     \n",
       "5   PFOA           μg   100000.0\n",
       "6   PFOS           μg   100000.0\n",
       "7   PFBA           μg   100000.0\n",
       "8   PFPeA          μg   100000.0\n",
       "9   PFHxA          μg   100000.0\n",
       "10  PFHpA          μg   100000.0\n",
       "11  PFNA           μg   100000.0\n",
       "12  PFDA           μg   100000.0\n",
       "13  PFBS           μg   100000.0\n",
       "14  PFHxS          μg   100000.0"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "InputFlow = pspd.IncomFlow()\n",
    "InputFlow.set_flow('ContaminatedWater',1000)\n",
    "InputFlow.calc()\n",
    "ContaminatedWater = InputFlow.Inc_flow\n",
    "ContaminatedWater.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in Advanced WWT (Using GAC)"
   ]
  },
  {
   "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>Effluent RO</th>\n",
       "      <th>RO Concentrate</th>\n",
       "      <th>Effluent GAC</th>\n",
       "      <th>Spent GAC</th>\n",
       "      <th>Effluent IonExchange</th>\n",
       "      <th>Spent IER</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>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>10.0</td>\n",
       "      <td>90.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Effluent RO  RO Concentrate  Effluent GAC  Spent GAC  \\\n",
       "PFOA   0.0          0.0             10.0          90.0        \n",
       "PFOS   0.0          0.0             10.0          90.0        \n",
       "PFBA   0.0          0.0             10.0          90.0        \n",
       "PFPeA  0.0          0.0             10.0          90.0        \n",
       "PFHxA  0.0          0.0             10.0          90.0        \n",
       "PFHpA  0.0          0.0             10.0          90.0        \n",
       "PFNA   0.0          0.0             10.0          90.0        \n",
       "PFDA   0.0          0.0             10.0          90.0        \n",
       "PFBS   0.0          0.0             10.0          90.0        \n",
       "PFHxS  0.0          0.0             10.0          90.0        \n",
       "\n",
       "       Effluent IonExchange  Spent IER  \n",
       "PFOA   0.0                   0.0        \n",
       "PFOS   0.0                   0.0        \n",
       "PFBA   0.0                   0.0        \n",
       "PFPeA  0.0                   0.0        \n",
       "PFHxA  0.0                   0.0        \n",
       "PFHpA  0.0                   0.0        \n",
       "PFNA   0.0                   0.0        \n",
       "PFDA   0.0                   0.0        \n",
       "PFBS   0.0                   0.0        \n",
       "PFHxS  0.0                   0.0        "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "AdvWWT.InputData.Trtmnt_optn['frac_to_RO']['amount'] = 0\n",
    "AdvWWT.InputData.Trtmnt_optn['frac_to_GAC']['amount'] = 1\n",
    "AdvWWT.InputData.Trtmnt_optn['frac_to_IonEx']['amount'] = 0\n",
    "AdvWWT.calc(Inc_flow=ContaminatedWater)\n",
    "AdvWWT.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADsCAYAAABKZHxbAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAvx0lEQVR4nO3deVxVdf7H8dflIsjughtL4Ipy2beJsARTcdzGylzCxjLTXLMZMZ1mFNPSn5mZS2jZL3NK1CxzKdRUUDEagTQXUikFd3Jn3y7n9wfD/UWKggIHuJ/n48FDzrnnnPs5N3h3+J7v+X41iqIghBCi7pmoXYAQQhgrCWAhhFCJBLAQQqhEAlgIIVQiASyEECqRABZCCJVIAIs6odFo9BqN5sjvvmb8d/3jGo3mxH/XWWg0mnf+u/yORqOJ0mg002qhlqkajcaypo8rRHVpqtMP2N7eXnF1da29akSjdfjwYXx9fe9Yn5GRgbW1NS1btjRs5+Pjg0aj4dKlS5iYmNC2bdsareXYsWN069YNU1PTGj2uEJVJSUm5pihKqz+ur9ZPoKurK8nJyTVXlTAa1tbWd/zsrF69munTpwPg5uZGdnY2hw8fRq/XM3PmTH7++Wesra2ZNm0aoaGhLFq0iICAAK5du0ZAQADp6eno9XpmzJhBfHw8hYWFTJw4kXHjxhEfH09UVBT29vYcP34cf39/PvvsM5YtW8a0adPQ6/U0b96cuLg4NT4OYWQ0Gk3G3dbLJYCoE/n5+fj4+BiWZ86cyZgxY0hISGDAgAEMGTIEKAvqI0eOABAVFXXf43788cfY2dmRlJREYWEhISEh9OnTByi7mj5x4gQODg6EhIRw8OBBpkyZwuLFi4mLi8Pe3r6mT1OIapEAFnXCwsLCEKw1adeuXRw9epRNmzYBcPv2bdLS0jAzMyMoKAgnJycAfHx8SE9Pp3v37jVegxAPSgJYNAimpqaUlpYCUFBQYFivKArLli0jPDy8wvbx8fGYm5sblrVaLSUlJXVTrBBVJL0gRIPg6upKSkoKgOFqFyA8PJzo6GiKi4sBOH36NLm5ufc8lo2NDdnZ2bVXrBBVJAEs6kR5G3D514wZM6q1/7Rp04iOjsbX15dr164Z1o8ZMwZ3d3f8/Pzw8PBg3Lhx973SHTt2LH379iUsLOyBzkWImlKtbmgBAQGK9IIQQojq0Wg0KYqiBPxxvVwBCyGESiSAhRBCJRLAQgihEumGJh6ItbU1OTk5D7RvaGgoly9fxsLCAoBOnTpV6NlQFWvWrCE5OZnly5c/UA1C1AcSwEIVn3/+OQEBd9yTEMKoSBOEeCiKohAZGYmHhweenp5s2LABKHsQIjQ0lCFDhtC1a1ciIiK4X4+bv/zlL6xduxaAVatWERERAUBSUhKPPfYY3t7eBAUFGfrwXrp0ib59+9K5c2fDmBIA48ePJyAgAJ1Ox+zZsw3rXV1dmT17Nn5+fnh6enLy5EkArl69Su/evdHpdIwZMwYXFxdDV7fPPvuMoKAgfHx8GDduHHq9voY+OSEo+wWq6pe/v78ihKIoipWVlaIoirJp0yalV69eSklJiXLlyhXF2dlZuXTpkhIXF6fY2toq58+fV/R6vfLoo48qBw4cUBRFUXr06KF06dJF8fb2Vry9vZVp06YpiqIoV65cUTp27Kjs379f6dy5s3L9+nWlsLBQad++vXLo0CFFURTl9u3bSnFxsfLJJ58o7du3V27duqXk5+crjzzyiHLu3DlFURTl+vXriqIoSklJidKjRw/lp59+UhRFUVxcXJSlS5cqiqIoK1asUF566SVFURRl4sSJyttvv60oiqLExsYqgHL16lUlNTVVGTBggFJUVKQoiqKMHz9e+fTTT2v9sxWND5Cs3CVTpQlCPJSEhARGjBiBVqulTZs29OjRg6SkJGxtbe85FsPdmiDatGnDm2++SVhYGJs3b6ZFixYcO3aMdu3aERgYCICtra1h+yeffBI7OzsA3N3dycjIwNnZmY0bN/Lhhx9SUlLC5cuXSU1NxcvLC4Cnn34aAH9/f7766ivDOWzevBmAvn370rx5cwD27NlDSkqK4b3z8/Np3bp1zX+IwmhJAIta8yBjMRw7doyWLVty6dKlBzr+2bNnWbRoEUlJSTRv3pwXXnihwtgR5ftUpR5FURg1ahTz58+/by1CPAhpAxYP5fHHH2fDhg3o9XquXr3K/v37CQoKeqBjHTp0iNjYWA4fPsyiRYs4e/Ysbm5uXL58maSkJACys7PvGZxZWVlYWVlhZ2dHZmYmsbGx933fkJAQNm7cCJSNrnbz5k2g7Ap706ZN/PbbbwDcuHGDjIy7DusqxAORK2DxUJ566ikSExPx9vZGo9GwcOFC2rZta7jBVZmIiAhDNzR7e3u++eYbXn75ZT755BMcHBx49913GT16NHv37mXDhg1MnjyZ/Px8LCws2L17d6XH9fb2xtfXl65du+Ls7ExISMh9z2H27NmMGDGCf//73wQHB9O2bVtsbGywt7dn3rx59OnTh9LSUpo0acKKFStwcXGp3ockRCVkLAhh9AoLC9FqtZiampKYmMj48eNrZexiYbwqGwtCroCF0Tt37hxDhw6ltLQUMzMzPvroI7VLEkZCAlgYvc6dO3P48GG1yxBGSG7CiQZjx44duLm50alTJxYsWKB2OUI8NAngBk6r1VYY6Lw8mA4cOIBOp8PHx4f8/HwiIyPR6XRERkYSFRXFokWLaryWJUuWkJeXV+PHBdDr9UycOJHY2FhSU1OJiYkhNTW1Vt5LiLoiTRANXGWTXX7++efMnDmTkSNHAvDhhx9y48YNtFptlWYbfhBLlixh5MiRWFpa1vixDx06RKdOnejQoQMAw4cPZ8uWLbi7u9f4ewlRV+QKuBFavXo1Gzdu5F//+hcREREMGjSInJwc/P39DWM1lAsNDaW8Z8u1a9dwdXUFyq44IyMjCQwMxMvLi1WrVgGVj/GwdOlSLl26RFhYWK1M9XPx4kWcnZ0Ny05OTly8eLHG30eIuiRXwA1c+Vxr5WbOnMmYMWNISEhgwIABDBkyBCgbPrL8SrkqV8Aff/wxdnZ2JCUlUVhYSEhICH369AHg8OHDnDhxAgcHB0JCQjh48CBTpkxh8eLFxMXFYW9vX9OnKUSjJAHcwFXWBPGwdu3axdGjRw3j9N6+fZu0tDTMzMzuOcZDbXF0dOT8+fOG5QsXLuDo6Fir7ylEbZMANnKmpqaUlpYCVBgzQVEUli1bRnh4eIXt4+PjH2iMh4cVGBhIWloaf582nb7hvVm/fj3r1q2r9fcVojZJG7CRc3V1JSUlBaDCrBTh4eFER0dTXFwMwOnTp8nNzb3nsWxsbAxj9dYGC2tbFr/7Dn369KFTly7odLpaey8h6oIEcANX3gZc/jVjxoxq7T9t2jSio6Px9fU1DEIOMGbMGNzd3fHz88PDw4Nx48bd90p37Nix9O3bt1Zuwmm1Wh4NDjYsDx06vMbfQ4i6JmNBiAZFo9Ewbtw4Vq5cqXYpQlRZZWNByBWwEEKoRAJYCCFUIgFshMofX/bw8GDgwIHcunXL8NqJEyfo2bMnbm5udO7cmblz5953Mk0hxIORADZC5X2Hjx8/TosWLVixYgVQdkNv0KBBzJgxg1OnTvHTTz/x/fff88EHH6hcsRCNkwSwkQsODjY80rtu3boKT7xZWlqyfPlyGXlMiFoiAWzE9Ho9e/bsYdCgQUBZ84O/v3+FbTp27EhOTg5ZWVlqlChEoyYBbITK+w63bduWzMxMevfurXZJQhglCWAjVN4GnJGRgaIohjZgd3d3w1Nx5c6cOYO1tTW2trZqlCpEoyYBbMQsLS1ZunQp7777LiUlJURERJCQkGCYdTg/P58pU6Ywffp0lSsVonGSADZyvr6+eHl5ERMTg4WFBVu2bGHevHm4ubnh6elJYGAgkyZNUrtMIRolGQ3NCOXk5FRY3rZtm+F7T09P4uPj67giIYyTXAELIYRKJICFEEIlEsBCCKESCWDBjh07cHNzo1OnTvLUmxB1SALYyOn1eiZOnEhsbCypqanExMSQmpqqdllCGAUJYCN36NAhOnXqRIcOHTAzM2P48OFs2bJF7bKEMAoSwEbu4sWLODs7G5adnJwMg/MIIWqXBLAQQqhEHsQwco6Ojpw/fx69Xo9Go+HChQs4OjqqXZYQRkEC2MgFBgaSmprKI64dsbaxQYueL774Qu2yhDAKEsBGztTUlCbmTTn76y8AtO/QEZ1Op3JVQhgHaQMW7Ivba/g+4cB+FSsRwrhIAAtDL4jmzZvj4OCgcjVCGA8JYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEsDNq1a6d2CVXSUOoU4n4kgAVZWVkALFy4UOVK7m3fvn0ATJw4UeVKhKgZEsACT09PAPr3769yJZVTFIXQ0FAA7O3t1S1GiBoiAWzEsrOzcXV15dy5c+zcuVPtcip1/vx5TEzKflQzMjJUrkZ9iYmJPPZEKLdv31a7FPGQZEaMGpSRkcEH0SvvWL/wfxaoUE3V+fn5s2dvHHv2xqldSgV//NyeHTqUFR9Eq1RN/fHzz6n8cDCBgU8NYe+uWExN5de4odIoilLljQMCApTk5ORaLKdh27NnD71698bapx+mNi0N62/tX6tiVZUzc3DDstOf1C6jUtk/7UR/OxPLLo9h1raT2uXUK01dfShMXMfQXo+y8oNlapcj7kOj0aQoihJwx3oJ4Jo19623eXfVWmyffQuTJk3VLkc0YqUFOdze8Dpv/2s6EydMULsccQ+VBbC0Adewf/5jJuHdA8jd+T6KUqp2OaIRM2lqjfWgN3j9jVl89913apcjHoAEcA3TaDSs/eRjXCyLyf0+Ru1yRCPXpLkD1n+expBhIzh16pTa5YhqkgCuBebm5uz8ZitN0g+Sm1q/bmyJxqfpI56YPRrBk+F/5vr162qXI6pBAriWtGrVit07vqXgwCcUXvxZ7XJEI2fl1YcCBz/6DXqKoqIitcsRVSQBXIt0Oh0bPv832d8spOT2b2qXIxo5q+5/Je1mMWPGjac6N9eFeiSAa1m/fv2IemMG2VvnUVqYp3Y5ohHTmGixDv8bW7/bz6J3F1d7/7feegudToeXlxc+Pj785z//qbHabt26xQcffFBjx2ssJIDrwN//9hpPhYeRs3MxSqle7XJEI2Zibon1wH8w5+0FbN++vcr7JSYmsn37dn788UeOHj3K7t27cXZ2rrG6JIDvTgK4Dmg0Gj5aFY1bS3NyD9TPhzJE42Fq1xqb/jMY8fwojh07VqV9Ll++jL29Pebm5kDZeBsODg64uroyffp0PD09CQoK4pdffgHg6tWrPPPMMwQGBhIYGMjBgwcBiIqKYvTo0YSGhtKhQweWLl0KwIwZM/j111/x8fEhMjKyFs66YZIAriNNmjThmy1fYXH5MLlHd6ldjmjkzB270vTx0fT+c39+++3+9x/69OnD+fPn6dKlCxMmTDCMPAdgZ2fHsWPHmDRpElOnTgXg1Vdf5bXXXiMpKYkvv/ySMWPGGLY/efIkO3fu5NChQ8yZM4fi4mIWLFhAx44dOXLkCO+8806Nn29DJQ+R16EWLVqwZ1csgY8+hrZZW5o+4qV2SaIRs3IPJefmRcL7DyLxQDxNm1b+ZKa1tTUpKSkcOHCAuLg4hg0bxoIFZWNxjBgxwvDva6+9BsDu3btJTU017J+VlUVOTg5QNqqeubk55ubmtG7dmszMzFo6w4ZPAriOdenShS83rmfws8PRDp1Pk+YOapckGjGrx0Zw7ttFPP/CaDbGfI5Go6l0W61WS2hoKKGhoXh6evLpp58CVNin/PvS0lJ++OGHu4Z6eTNG+TFLSkpq6nQaHWmCUEGvXr1Y+PZccrbOQ1+Qo3Y5ohHTaEyw6jOF3YmHmTvvrUq3O3XqFGlpaYblI0eO4OLiAsCGDRsM/wYHBwNlTRbLli2rsP292NjYkJ2d/aCn0WhJAKtkwvjxPPf0X8j99h0UvVwhiNpj0qQpVgNmsnDJcr74YtNdt8nJyWHUqFG4u7vj5eVFamoqUVFRANy8eRMvLy/ef/993nvvPQCWLl1KcnIyXl5euLu7s3LlncOw/l7Lli0JCQnBw8NDbsL9jlGMhmZtbW1on6pLs2bN4oknnqBXr153fV2v1/NkeD+OZ5lj3XPcPY91IXo0JmYWoNGgtWqOff+/obVu/v/r/ztgeYveE2jq1I2804lc3fwWDmOiadKyrDuRopRyc89HFGT8BGjQmJph/5fXadKsbY2et6ifCq/8QvbXc9i35zv8/f2rtI+rqyvJyckyC8lDqmw0NGkDrkVvvvnmPV/XarVs/eoLvP2DyPlxO1Z+A+65fZsRb6O1tOPmvk+5/cNGWvQaV2H97+X+vB9zJ3dyU/fT7PEIAPJ+PoA++zrtRi9HozGhJOsaGjMZMtNYmLfthD5sPOH9B3L0x2QcHOT+g9qMqgkiPj6e0NBQhgwZQteuXYmIiDA8spmUlMRjjz2Gt7c3QUFBZGdnU1BQwIsvvoinpye+vr7ExZUNrLNmzRoGDx5M7969cXV1Zfny5SxevBhfX18effRRbty4AcALL7zApk1lf/K5uroye/Zs/Pz88PT05OTJkwAUFhbi2MaeG3s/InPjLC5Ev4g+795TzTR19qD45qVKXy8tyqfwwgla/vlVck/uN6zX59xAa90CjabsP7uprT3aptYP+GmKhsjS7THo1odeffuTl3f/JzPT09Pl6rcWGVUAAxw+fJglS5aQmprKmTNnOHjwIEVFRQwbNoz333+fn376id27d2NhYcGKFSvQaDQcO3aMmJgYRo0aRUFBAQDHjx/nq6++IikpiTfeeANLS0sOHz5McHAwa9fe/WELe3t7fvzxR8aPH8+iRYsAmDNnDv3792dffDxK5mn0WVfvew75vx7CrJWrYTkz5h9c+mQyl9f+DYC8tB9o2t6fJi0c0Ta1ofBKWed5y66Pk/fLIS59Mpkbe1dTlPnrw3yUooGyChpCpok9Q58bSWmpjFmtJqML4KCgIJycnDAxMcHHx4f09HROnTpFu3btCAwMBMDW1hZTU1MSEhIYOXIkAF27dsXFxYXTp08DEBYWho2NDa1atcLOzo6BAwcCZTMMp6en3/W9n376aQD8/f0N2yQkJDB8+HC6d+/OyhXLQaNBn3/3u8XlQVtamIfdo88a1rcZ8TYOLy6j3V/Lnv/P+3k/Vt2eAMCy2xPkpZZ1qje1tcfx5VU0f2IUGo0JmevfID/9yAN+kqKh0mg0WPeayMGjacx841+VbpeZmclzzz1Hhw4d8Pf3Jzg4mM2bNxtenzp1Ko6OjneE+Nq1a/Hw8DD85Vh+sSHuZHRtwDXVR/H3xzExMTEsm5iYVHrM8m0qe98XXhjFK+NfIe+7ZTR5dh4abZMKr9+trfeP9PnZFGQcpehqOqABpRTQ0CxsNBqNBo1pEyw6BmDRMQATy2bkp/2AhatPlc9bNA4a0yZY95/BB6un4+HejeefH1nhdUVRGDx4MKNGjWLdunVA2aSzW7duBcr6AW/evBlnZ2f27dtHWFgYALGxsSxZsoRdu3bh4OBAYWFhpX8RCiO8Ar4bNzc3Ll++TFJSElA2XXtJSQmPP/44n3/+OQCnT5/m3LlzuLm51eh7h4SEsHHjRgB27dpFYUEBfp0dyd2z8oGGFMw7dRArXRhO4z/Bafz/4jRhDabN2lB44QSFV36hJLtswG5FKaX46lm0tq1r9HxEw6G1aobNoH8wfvKrJCYmVnht7969mJmZ8corrxjWubi4MHnyZKDsfopOp2P8+PHExPz/zC/z589n0aJFhht85ubmvPzyy3VwNg2TBDBgZmbGhg0bmDx5Mt7e3vTu3ZuCggImTJhAaWkpnp6eDBs2jDVr1lS48q0Js2fPZteuXXh4ePDFF1/Qtm1bvli/jpaFl8hL3nz/A/xB7s/7sOwSXGGdZZfHyE3dR2neba5++SaXPp7A5f+dBCZabP3v3fNCNG5mrVyx6j2Z/oMGk5GRYVh/4sQJ/Pz8Kt0vJiaGESNG8NRTT/HNN99QXFwMlN0bqWoXN2Ek/YDrs8LCQrRaLaampiQmJjJ+/HiOHDnChQsX8PYPRNv9ZSw719+p40XjcPObRQwL6caHK8uGjFy6dClnz541PHgxceJEEhISMDMz4+DBg7Rv356TJ09iY2PD008/zejRoxkwYAAtWrTg7Nmz2Nndu6nM2Eg/4Hrq3LlzDB06lNLSUszMzPjoo48AcHJyYsf2rYT17oupXSvMWndQuVLRWOWfScHk8glmvv6/hnU6nY4vv/zSsLxixQquXbtGQEAAO3fu5NatW3h6egKQl5eHhYUFAwYMQKfTkZKSQs+ePev8PBoiaYJQWefOnTl8+DA//fQTSUlJhp4YAIGBgaxe+QHZW9++b99gIR5E0bVz5O56n+1bNtO+fXvD+p49e1JQUEB0dLRhXXm/4ZiYGFavXk16ejrp6emcPXuW7777jry8PGbOnElkZCRXrlwpO35REatXr67bk2pAJIDrOXf3bpQW5VNalK92KaKR0efdJmfrW6x4/z1CQkIqvKbRaPj666/Zt28f7du3JygoiFGjRjFnzhx27NhB//79DdtaWVnRvXt3tm3bRr9+/Zg0aRK9evVCp9Ph5+dHVlZWXZ9agyFtwHVo9OjRbN++ndatW3P8+PH7bp+ZmYmXXwCl/iOwcu9RBxUKY6GUFJP11SxeGtKPRQsXqF1Oo1dZG7BcAdehF154gR07dlRp24KCAvr0G4i+4xMSvqJGKYpC7t5oHtW1Z+GCt9Uux6hJANehJ554ghYtWtx3O0VReO75UVwotsQqeHgdVCaMSW7SV9gXXWHT+nWYmEgEqEk+/Xoo6s257D10DKveUwwD5whRE/JOJ8KJWHbv+AYrK6sq77djxw7c3Nzo1KmTYaqi38vIyODJJ5/Ey8uL0NBQLly4UJNlN1ry213PbNiwkcXLPsB64ExMmtTsQx/CuBVlniFv7wfEbt+Kk5NTlffT6/VMnDiR2NhYUlNTiYmJqTAfHMC0adP461//ytGjR5k1axYzZ86s6fIbJQngeiQ5OZmXxr2CzcB/YGp9/6YKIaqqJOcG2dve4uNVKyt0dayKQ4cO0alTJzp06ICZmRnDhw9ny5YtFbZJTU019P0NCwu743VxdxLA9cSFCxfo238glj0nYNamo9rliEaktLiQ3O3zeW3SeIYNG1rt/S9evIizs7Nh2cnJiYsXL1bYxtvbm6+++gqAzZs3k52dzfXr1x+ucCMgAVyHRowYQXBwMKdOncLJyYmPP/4YgNzcXHr17Q+6P98xjoMQD0NRFHJ3LyMs0JM5s2fV2vssWrSIffv24evry759+3B0dESr1dba+zUW8ihyHfr9qFHlSktLeXZ4BFebtMY68GkVqhKNWe4P63HU5vD52m/vOSX9vTg6OnL+/HnD8oULF3B0dKywjYODg+EKOCcnhy+//JJmzZo9cN3GQq6AVfb6zDdIPHEG6ycnPvAviBB3k/vzfkzS4tn17TYsLCwe+DiBgYGkpaVx9uxZioqKWL9+PYMGDaqwzbVr1wwDs8+fP5/Ro0c/VO3GQgJYRZ9+upaVn6zFesAMNKZN7r+DEFVUeOkUBftXs3vHt7Rt+3CzXpuamrJ8+XLCw8Pp1q0bQ4cORafTMWvWLMMA7fHx8bi5udGlSxcyMzN54403auI0Gj15FFklBw8eJLz/IGyfmYtZKxe1yxGNSEnWVW6tn866NavvuFIV6pBHkeuR9PR0BvzlKaz6vCrhK2pUaVE+Odve4o3Xp0n4NgASwHUsKyuLJ8P7ofV9GosOMnOAqDmKUkruziX0C32MGdMj1S5HVIEEcB3S6/U8NWQYt207YOUnUwGJmpWb8G862MKajz+SG7oNhHRDq0NTXvs7h9N/w+Yvs+QXRNSonGO7MT+fROyPSZiZmaldjqgiCeA6snLVKv69cTN2wxei0crHLmpOwfnjFH2/lgPfJ2Bvb692OaIapAmiDsTFxTHt9X9gM+gNtE2t1S5HNCLFt66Q8+07fLF+Hd26dVO7HFFNEsC1LC0tjcHPPIv1n6fRpIXj/XcQoopKC3PJ2foW86JmEx4ernY54gHI38K16ObNmzzZ5880+dMImrp4qV2OqIJr3ywh/9dD1dqnNL/qc541C32hmhXdXVMXHwp/WMezA/sw9dXJNXJMUfckgGtJcXExA/7yNHltvLD27qt2OaKK9DnXqxWo1fWiX8uHPsbJkyfZ8dk0Huv+BCtXLK+BqoRaJIBrgaIojB0/kZ+v5mMz8DW1yxHV0GbY3IfaX5+fRd7x3RQd3sZzw59l8TsLqzXzRFV8//33XL91m+1ff4mpqfwKN2TyKHItWPze+8xZtBS7oQswMbdUuxyhAn1+NnnxH9FKf40Dcbtp1aqV2iUJFcmjyHUkNjaWWW/Ow2bQGxK+RkxrYYN139e4ZufGEz17U1BQoHZJoh6Sv19qUGpqKk8/O5SmAUMoyb5GSfY1AHIOx1JamKtKTeaO3TB31lVYp8+6RsG5Y1h27V6vR2EryjxDye1MLDs/CvX4wRWNxgQTq+Y0adb2D+s1WHX/K799u5CZb/yT995dpFKFor6SJoga9MMP/2HcpCn88SM99mP17qrXNE+/oArLf6zHpWMXbO2a1WFFVfPHOnU+AfVyGvUSvZ5L58+hbdYWrfegO2Y10efe4ua/p5CUmCB9dY1UZU0QEsBG7LPPPuP5558H4MqVK7Rp00bliu6kKApz585l9uzZQFnvkvp446mkpIRt27Yx6dW/UdTxCaz+VHHuteyEzxji1ZJVH0ivBWMkASzuqrS01DB3V3V+FurarVu3aN68Oe3bt+fMmTNql1OpzMxMvP0CUbqPxaK9r2F9SdZv3PrsNbJu3aiXV/GidslNOHFXJiYmnDx5EoDTp0+rXE3lmjVrxrJlyzh79iwlJSVql1OpNm3a8PbcKPTHvqmw3tS2NSbmFvz6668qVSbqIwlggZubGwAvvfSSypXc28SJEwGIjo5WuZJ7Gz58ONlnj1JamFdhfdM2HThx4oRKVYn6SAJYGJRfCddX5UN41vc6LS0tsbFrRmlBToX1miZNyc1VpzeMqJ8kgIWBXq9Xu4QqaQh1mmi1KEppxZUak3rdzi7qngSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJICFEEIlEsBCCKESCWAhhFCJBLAQQqhEAlgIIVQiASyEECqRABaiBl28eBHbZi24deM6JmYWFV5Tmtry0stjmTD5VZWqE/WNBLAQldBqtfj4+Bi+FixYAMCBAwfQ6XT4+PiQn59PZGQkOp2OyMhIVq5ciZtbF5ra2aMU//9EnIqiUGpmjbWNLaNHPV/tWpYsWUJeXt79NxQNSv2b20WIesLCwoIjR47csf7zzz9n5syZjBw5EoAPP/yQGzduoNVqiYqKYuiQZ9CamjHrzdex7j8d87adydmzgtYlV9l9OAVHR8dq17JkyRJGjhyJpaXMtN2YyBWwENWwevVqNm7cyL/+9S8iIiIYNGgQOTk5+Pv7s2HDBqBs3OK/vfYqHZzbkb1lHjdjpuHXugl5WTdxdHREr9cTGRlJYGAgXl5erFq1CoD4+HhCQ0MZMmQIXbt2JSIiAkVRWLp0KZcuXSIsLIywsDA1T1/UMLkCFqIS+fn5+Pj4GJZnzpzJmDFjSEhIYMCAAQwZMgQAa2trw5VyVFSUYfsWLVrwz3/+k6ysbAYNGkhQUNns1B9//DF2dnYkJSVRWFhISEgIffr0AeDw4cOcOHECBwcHQkJCOHjwIFOmTGHx4sXExcVhb29fJ+cu6oYEsBCVqKwJojo6dOhAQEAA165dM6zbtWsXR48eZdOmTQDcvn2btLQ0zMzMCAoKwsnJCQAfHx/S09Pp3r37Q9Ug6i8JYCFqiampKaWlZbNiFBRUvCG3bNkywsPDK2wfHx+Pubm5YVmr1dbrCUjFw5M2YCFqiaurKykpKQCGq12A8PBwoqOjKS4uBspmo77fXHE2NjZkZ2fXXrFCFRLAQlSivA24/GvGjBnV2n/atGlER0fj6+tboQlizJgxuLu74+fnh4eHB+PGjbvvle7YsWPp27ev3IRrZDTVmSQwICBASU5OrsVyhFo0Gg3Nmzfnxo0bapdyTxqNhnHjxrFy5Uq1SxGiyjQaTYqiKAF/XC9XwEIIoRIJYCGEUIkEsBAPITMzk+eee44OHTrg7+9PcHAwmzdvNrw+depUHB0dDb0hyq1duxYPDw88PT3x9fVl0aJFdV26qAckgIV4QIqiMHjwYJ544gnOnDlDSkoK69ev58KFCwCUlpayefNmnJ2d2bdvn2G/2NhYlixZwq5duzh27Bg//PADdnZ2ap2GUJEEsBAPaO/evZiZmfHKK68Y1rm4uDB58mSgrF+vTqdj/PjxxMTEGLaZP38+ixYtwsHBAQBzc3Nefvnlui1e1AsSwEI8oBMnTuDn51fp6zExMYwYMYKnnnqKb775xtDv9/jx4/j7+9dVmaIekwAWooZMnDgRb29vAgMDKSoq4ttvv2Xw4MHY2trypz/9iZ07d6pdoqhn5FFkIR6QTqfjyy+/NCyvWLGCa9euERAQwM6dO7l16xaenp4A5OXlYWFhwYABA9DpdKSkpNCzZ0+1Shf1hFwBC/GAevbsSUFBAdHR0YZ15YOmx8TEsHr1atLT00lPT+fs2bN899135OXlMXPmTCIjI7ly5QoARUVFrF69WpVzEOqSABbiAWk0Gr7++mv27dtH+/btCQoKYtSoUcyZM4cdO3bQv39/w7ZWVlZ0796dbdu20a9fPyZNmkSvXr3Q6XT4+fmRlZWl4pkItcijyAKQR5GFqE3yKLIQQtQzEsCiXiqfENPDw4OBAwdy69Ytw2s3btygZ8+euLm50blzZ+bOnUt1/pITor6QABb1UvlsFMePH6dFixasWLHCMOfa5s2bGTp0KKdOneKnn37i+++/54MPPlC5YiGqTwJY1HvBwcGcP3+eqX+PBKCkpISLVzIBsLS0ZPny5YYp49W0Y8cO3Nzc6NSp013r2b9/P35+fpiamlYYoF0YLwlgUa/p9Xr27NnD4MGDOZx8CMdHXPhzv368OXuWYZuOHTuSk5Ojak8CvV7PxIkTiY2NJTU1lZiYGFJTUyts88gjj7BmzRqee+45laoU9Y08iCHqpfLZKC5evEi3bt3o3bs3Wq2Woc88jYuLCxqNRu0SKzh06BCdOnWiQ4cOAAwfPpwtW7bg7u5u2MbV1RUAExO57hFl5CdB1EvlbcAZGRkoisKKFSsAcHd3N8yzVu7MmTNYW1tja2urRqkAXLx4EWdnZ8Oyk5MTFy9eVK0e0TBIAIt6zdLSkqVLl/Luu+9SUlJCREQECQkJ7N69Gyi7Up4yZQrTp09XuVIhqk8CWNR7vr6+eHl5ERMTg4WFBVu2bGHevHm4ubnh6elJYGAgkyZNUrVGR0dHzp8/b1i+cOECjo6OKlYkGgJpAxb1Uk5OToXlbdu2Gb739PQkPj6+jiu6t8DAQNLS0jh79iyOjo6sX7+edevWqV2WqOfkCliIGmBqasry5csJDw+nW7duDB06FJ1Ox6xZs9i6dSsASUlJODk58cUXXzBu3Dh0Op3KVQu1yVgQAmg4Y0EI0RDJWBBCCFHPSACLh3LlyhWGDx9Ox44d8ff3p1+/fpw+fVq1epYsWWIYk7c61qxZw6VLl2qhIiEqJwEsHpiiKDz11FOEhoby66+/kpKSwvz588nMzFStpnsFsF6vr3Q/CWChBglg8cDi4uJo0qRJhVmBvb296d69O5GRkXh4eODp6WkYRCc+Pp7Q0FCGDBlC165diYiIMIxilpSUxGOPPYa3tzdBQUFkZ2ej1+uJjIwkMDAQLy8vVq1adc/jLF26lEuXLhEWFkZYWBgA1tbW/P3vf8fb25vExETefPNNAgMD8fDwYOzYsSiKwqZNm0hOTiYiIgIfHx/y8/NJSUmhR48e+Pv7Ex4ezuXLl+v40xVGQVGUKn/5+/sronEClObNm1drn/fff1+ZOnXqHes3bdqk9OrVSykpKVGuXLmiODs7K5cuXVLi4uIUW1tb5fz584per1ceffRR5cCBA0phYaHSvn175dChQ4qiKMrt27eV4uJiZdWqVcrcuXMVRVGUgoICxd/fXzlz5kylx1EURXFxcVGuXr1a4bw2bNhgWL5+/brh+5EjRypbt25VFEVRevTooSQlJSmKoihFRUVKcHCw8ttvvymKoijr169XXnzxxWp9NkL8HpCs3CVTpR+wqHEJCQmMGDECrVZLmzZt6NGjB0lJSdja2hIUFISTkxMAPj4+pKenY2dnR7t27QgMDAQwPFK8a9cujh49ahg57Pbt26SlpWFmZnbX43Tv3v2OWrRaLc8884xhOS4ujoULF5KXl8eNGzfQ6XQMHDiwwj6nTp3i+PHj9O7dGyhrumjXrl0Nf0pCyIMY4iHodLpqD6tobm5u+F6r1VJSUlLptoqisGzZMsLDwyusj4+Pr/JxmjZtilarBaCgoIAJEyaQnJyMs7MzUVFRFBQU3PV9dTodiYmJ1To3IapL2oDFA+vZsyeFhYV8+OGHhnVHjx6lWbNmbNiwAb1ez9WrV9m/fz9BQUGVHsfNzY3Lly+TlJQEQHZ2NiUlJYSHhxMdHU1xcTEAp0+fJjc395412djYkJ2dfdfXysPW3t6enJycCv/z+P1+bm5uXL161RDAxcXFnDhx4n4fhxDVJlfA4oFpNBo2b97M1KlT+Z//+R+aNm2Kq6srS5YsIScnB29vbzQaDQsXLqRt27acPHnyrscxMzNjw4YNTJ48mfz8fCwsLNi9ezdjxowhPT0dPz8/FEWhVatWfP311/esaezYsfTt2xcHBwfi4uIqvNasWTNefvllPDw8aNu2raHJA+CFF17glVdewcLCgsTERDZt2sSUKVO4ffs2JSUlTJ06VZ5cEzVOnoQTgDwJJ0RtkifhjMxbb72FTqfDy8sLHx8f/vOf/9TYsW/duiVzsAlRAySAG6HExES2b9/Ojz/+yNGjR9m9e3eFwcL/aOPGjQDcvHnTMHDMvUgAC1EzJIAbocuXL2Nvb2/oKWBvb4+DgwOurq5Mnz4dT09PgoKC+OWXXwCI/vBjw74vvDiagwcPAhAVFcXo0aMJDQ2lQ4cOLF26FIAZM2bw66+/4uPjQ2RkZJ2d1/0mvRSioZEAboT69OnD+fPn6dKlCxMmTGDfvn2G1+zs7Dh27BiTJk1i6tSpALSxb87TQ4by+j/e4MeUZMaMGWPY/uTJk+zcuZNDhw4xZ84ciouLWbBgAR07duTIkSO88847dXJOVZn0UoiGRgK4EbK2tiYlJYUPP/yQVq1aMWzYMNasWQPAiBEjDP+Wd7Pau3cvv6adYsc32xk8eDBZWVmGAdH79++Pubk59vb2tG7dWrVxHn4/6aWZmZlh0kshGjLphtZIabVaQkNDCQ0NxdPTk08//RSgwmzC5d+Xlpbyww8/0LRp0zuOU50HJ2rT3Sa9rMkbi0KoQa6AG6FTp06RlpZmWD5y5AguLi4AhoFxNmzYQHBwMFDWZLFs2bIK29/LvR52EEJUnQRwI5STk8OoUaNwd3fHy8uL1NRUoqKigLKeDl5eXrz//vu89957ACxdupTk5GS8vLxwd3dn5cqV9zx+y5YtCQkJwcPDo85uwsmkl6IxkgcxjIirqyvJycnY29urXUq1lZSU0KVLF/bs2YOjoyOBgYGsW7dOnk4TDUJlD2JIG7BoEH4/6aVer2f06NESvqLBkytgIYSoZfIosgDKuqgBpKenY2FhgY+Pj+Fr7dq1QFlThaenJ15eXvTo0YOMjAw1Sxai0ZImCCNW/jDF3cTFxWFvb8/s2bOZN28eH330Ud0WJ4QRkCtgcU/BwcFcvHhR7TKEaJQkgI1Y+XgO5V8HDhy4Y5sdO3YwePDgui9OCCMgTRBG7F5NEGFhYdy4cQNra2vmzp1bt4UJYSTkCljcVVxcHBkZGfj4+DB79my1yxGiUZIAFpUyNTVlyZIlrF27VmbKEKIWSAAbsT+2AZeP9/t77dq1Y8SIEaxYsUKFCoVo3ORBDCGEqGXyIIYQQtQzEsBCCKESCWAhhFCJBLCQyS6FUIkEsJGTyS6FUI8EsJGTyS6FUI8EsJG722SXMviOEHVDAlgIIVQiAWzkZLJLIdQjAWzkAgMDSUtL4+zZsxQVFbF+/XoGDRqkdllCGAUZjtLIyWSXQqhHxoIQQohaVtlYENUKYI1GcxWQGRqFEKJ6XBRFafXHldUKYCGEEDVHbsIJIYRKJICFEEIlEsBCCKESCWAhhFCJBLAQQqhEAlgIIVQiASyEECqRABZCCJVIAAshhEr+DyHVu+0ZpumjAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "AdvWWT.plot_sankey(gap=0.8)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABVy0lEQVR4nO3de3yP9f/H8cd7m81so5jjyCFSxsxhSnw1OUxOKSriG0kHlOobfcmvqBT1VckhKZ2/hVJyyuGrEMXXiMi+RYWcj2EnO75/f3zmkznMxrZrn8+e99ttN/t8rutzXc/rs3ntdV3X+3NdxlqLiIiIiIgn8HE6gIiIiIhIbql5FRERERGPoeZVRERERDyGmlcRERER8RhqXkVERETEY6h5FRERERGPoea1ABhjrDGmttM5csMYc5UxJsEY4+t0FhERJ6hmi3iWYtW8Zv2HP/2VaYxJPuNx7wu8JtoYsycfM6wwxgzIr+VdLmvtH9baYGttRn4v2xjzvjEmNev9PWaM+Y8x5tqsaaONMWln/UyePOO1nY0x64wxicaYo8aYj40xVc+zjtFZf3iuz2O2O40x3xtjkowxK3Ix/93GmF1Zeb40xpQ9Y1pZY8ycrGm7jDF3nzGtnzEm46ztTDDGVMlj3pzWEZ31+3zm8vvmZfkiRZFq9rk8tWYbY/yNMa8YY/ZkvXanMWZCHrI9b4zZYoxJN8aMvsi8xhjzUlaOo1nfmzOmRxpjNmTV/w3GmMgLvAenv37Mbc5cruN872WtvK6jOCtWzWvWf/hga20w8AfQ5YznPnY6n5d6Oev9rgocAt4/Y9qsM38m1tqXAYwxPYBPgAlAKBAOpACrjTFXnn5xVjG6BziW9W9eHMta/riLzWiMCQemAX8HKgJJwBtnzDIFSM2a1huYmvWa09actZ3B1tp951nP6ByK8sXWse+s5X9wse0SKepUsx1RUDV7BNAUaAaEANHAD3nI9SvwJLAwF/M+AHQDGgIRQBfgways/sBc4N/AlcAHwNys5097+aztbHi+lWTt2ESf5/ncrOPs9/L3XGyXZClWzeuFGGMCjDETjDH7sr4mZD0XBCwCqpx5xMwY08wYs8YYc9wYs98YM/msX8rcrjc6ay/0CWPMoaxl3XvG9MCsPdVdxpgTxpjVxpjArGldjTFbszKsMMZcd8brdhpjhhljNmftBb9jjKlojFlkjIk3xiw7XVCMMTWM68ilX9bjFVl7uN9lzbvUGBN6xrLvycpz1BjzdNa62l5sW621SbiKW/2LvCcGeAUYY639xFqbbK09AAwAEoDHz5j9b0BlYAjQMy8/A2vtMmvtp8A5TeR59AbmW2u/tdYmAE8DtxtjQrJ+R7oDT1trE6y1q4F5uBrdfJGf6zj75531nPvIkjHGN+t37ogxZocx5uGz5xdxmmq2R9bsKGCOtXafddlprf0wV2+8K88H1tpFQHwuZu8LvGKt3WOt3ZuVr1/WtGjAD5hgrU2x1k4EDHBzbrPkQr6tw7jO3q0+6zn3MBdjTDljzHxjzEljTKwxZszZ83sjNa8uI4EbgEhce2rNgP+z1iYCt5D9qNY+IAPXf8hQoDnQBhh0ieuuBJQBwoD7gCnmrz3V8UAT4EagLK69zkxjzDXADOAxoDzwFTD/rGLcHWgHXINrr3MR8FTW/D64Gr4LuRu4F6gA+ANDAYwx9XAdceyNq2k8nfuijDHBWa/beJFZ6wJXAZ+d+aS1NhP4PGubTusLzAc+zXrc5Yz13W2M2ZybbLkQDrhPG1lrf8N1FPSarK90a+22M+b/Mes1+SU366hgjDmY1XC+lvVH/FLcj+t3PhJojOvohUhRo5qdnSfU7LXAP4wxg4wxDbKa3jPX94Yx5swzWpcjW80me70MBzZba+0Z0zeTvzU7N+voYlxDM7YaYwZexrqmAIm4fi/7Zn15PTWvLr2B56y1h6y1h4FnyeGolrV2g7V2rbU23Vq7E9cp5Zsucd1pWetOs9Z+hWtPta4xxgfoDzxqrd1rrc2w1n5vrU0B7gIWWmv/Y61Nw1UwA3EVzNMmWWsPZu11rgL+a63daK09BcwBGuWQ6T1r7TZrbTKuxjAy6/keuI5ArrbWpgLPAPYCyzhtqDHmOK5TPsH8tfcLcGfWUYjTX1Vw/XEB2H+eZe0/Pd0YUwq4A/gk6z2YzRlDB7KOAERcJFtuBQMnznruBK5TX8HAyQtMO+2Gs7bzt0tYf07r+BnXz6gyrj37JsCreVzHaXcCr2cdsfiTXAyrEHGAanZ2Rb5mA2OBl3D97NYDe80ZY/OttYOstZe6Q3G2s2v2CSA4q2HOqZ6fNvSs7czrMKyLreNT4DpcOyb3A88YY3rlcR0Y14f2ugOjrLVJ1to4XEMUvJ6aV5cqwK4zHu/Keu68jDHXGGMWGGMOGGNOAi/y13/QvDpqrU0/43ESrl/8UKAkcL5GJ1verD3c3WTfoz54xvfJ53kcnEOmA+fJc3q9u89YbxJwNIflAIy31l5hra1kre2addTytE+zpp3+2gccyZpW+TzLqnzG9NuAdFxHMAA+Bm4xxpS/SJ5LkQCUPuu50rhOX+U07bS1Z23n1acnZP0eHc/6YzEcGH5GwVyQi/VjrT1grY2z1mZaa3fgOtrT/RK3NdvP+KzvRYoK1ezsinzNzmrmp1hrWwBXAC8A7545fCIfnV0zSwMJWUdCc1Ozx5+1ne4m+8ymFmgJLDjjueEXWH+2dWTV632nd3CA13HtaORVeVzDE4pdzVbz6rIPqH7G46v4ayzk+fZSp+I62lXHWlsa16kdc575LscR4BRw9XmmZcubtTdZDdibzxnOth/XIP7T6w0EyuXzOn4B9uA6quqWdVSjO/B11lN9cRXoP4wxB3CdsiqB6/RZftuK69Tk6Sy1gABgW9aXnzGmzhnzN8x6zUVZazufLpC4jnKOO6Ngds6aLa/rsFz8/3bgGd+XOeP7bD9jXL9XIkWNanbuFKWa7WZd42KnAH8C9fI5D5xVs8leL7cCEWcNW4gg9zX7ijNq9mqg8xnPnT5Tldd1WHL+fXTXa2PMmfX6MK6DOMWuZqt5dZkB/J8xpnzWQPdncH1KEFx7v+XO+oUJwXUaN8G4LiNyOeNVzitrz/xd4FXj+sCBrzGmuTEmANcph07GmDbGmBLAE7g+2fl9fuc4y2xc43RuzBqrNZp8/gOQtWc8FNfP425jTEljTCVgOq4919eMMWG4xqx1xnV6LBJXcXqJXF51IOv9LIlrr9Unaz0lLjD7x7i2+29ZY0mfA76w1sZb1xi7L4DnjDFBxpgWwK3AR5ey/edzsXUYY1obY6obl2q4muC5F1lsf2OMjzGmHa6iGpK1/Z8CjxpjwowxVwD/zK/tEMlHqtm5UyRqNoAx5jHj+sBboDHGL2vIQAgXH1NL1utLZNVsH1w78yXNha91+yGu8bVhWUMbnuCvqyaswDUGeohxfcjv4aznv8njpuckx3UYY241xlyZVbOb4RrPnFPNbmBcHzr05a8PwJW1rsulfQGMNsaUyvrdzuuVdzySmleXMbjG4GwGtuC6fMcYAGvtz7gK5e/mrzE+Q3Ed4YsH3gZmFVCuoVl5YnFd2uklwMda+wvQB5iEa2+/C65LyKQWUA4ArLVbgUeAmbj26BNwXUolJZ/XMwvX+LXHcZ3iisO159nCWns0a9oma+3SrFPmB6zr060Tce3t1jfG9DbG5LQn/Xdcp+Km4rpqQTKunyXgvr7k387Y7odwNbGHcBXcM8dmDcrKdwjX78rArNec1tyce53XqDy+LTmtoxGuP4KJWf9uIecPd4Dr6NBhXM3ps7j++F+D6z1Yiuv/wkZcwzLScRVikaJCNTsXilDNBtdwhldwDXE4AgwGutusS0QZY940xryZw2rexlWne+H6wF5y1jrJOrCQcMa803B9mHcL8BOuy2tNy8qaiuuDqPcAx3GNU+521s/iybPq9RHyIBfr6IlrTHE8rkb7JZvz5Q3jcP1+H8J15Pwz/vodfhjX2bMDuA5ozCCff75FkbH2fGdYRC7OuD6NehzXqbgdDseRXDDG1AB2ACXOGrd3oflvAd601la/2LwiUrSpZnseY0w/YIC1tmUu538JqHTmOF1vpCOvkifGmC5ZpyeCcH1idguw09lUkl+yTul1zDqtFwaMwvVJZxHxQKrZ3s0Yc60xJuKMIQj3UQxqtppXyatbcX34YB9QB+hpdfjemxhcwwj+xDVs4H+4hhSIiGdSzfZuIbjGvSbiGkrwChf/zIPH07ABEREREfEYOvIqIiIiIh4jT82rMWZxQQURESlIxbF+FcdtFhHvkFP98svLgkqXLh3TtGlTjTMQEU909i12vZ5qtoh4sAvW7Dw1r3Xq1GH9+vWXH0dEpJAZY7Y7naGwqWaLiKfKqWZrzKuIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1r5Irvr6+REZGur/GjRsHwKpVqwgPDycyMpLk5GSGDRtGeHg4w4YNY/To0YwfPz7fs0yYMIGkpKR8X66IiLdQzRZvlqfbwzZt2tTqsivFU3BwMAkJCec8/9BDD9GyZUv69OkDQJkyZTh27Bi+vr6MHj2a4OBghg4dmq9ZatSowfr16wkNDc3X5Yp3M8ZssNY2dTpHYVLNLr5Us8XT5VSzdeRVLtn06dP59NNPefrpp+nduzddu3YlISGBJk2aMGvWrGzzRkdHu683eeTIEWrUqAFARkYGw4YNIyoqioiICKZNmwbAihUriI6OpkePHlx77bX07t0bay0TJ05k3759tG7dmtatWxfq9oqIeDLVbPEWebpJgRRfycnJREZGuh+PGDGCAQMGsHr1ajp37kyPHj0A197+pk2bABg9evRFl/vOO+9QpkwZYmNjSUlJoUWLFrRv3x6AjRs3snXrVqpUqUKLFi347rvvGDJkCK+++irLly/XXryIyAWoZos3U/MquRIYGOgucPlp6dKlbN68mdmzZwNw4sQJtm/fjr+/P82aNaNq1aoAREZGsnPnTlq2bJnvGUREvI1qtngzNa9SKPz8/MjMzATg1KlT7uettUyaNImYmJhs869YsYKAgAD3Y19fX9LT0wsnrIhIMaeaLUWZxrxKoahRowYbNmwAcO+xA8TExDB16lTS0tIA2LZtG4mJiTkuKyQkhPj4+IILKyJSzKlmS1Gm5lVy5fT4qdNfw4cPz9Prhw4dytSpU2nUqBFHjhxxPz9gwADq1atH48aNqV+/Pg8++OBF99YfeOABOnTooMH/IiIXoJot3kyXyhKRYkGXyhIR8Ry6VJaIiIiIeAU1ryIiIiLiMdS8ioiIiIjH0KWyiqEL3TYwN6Kjo9m/fz+BgYEA1K5dO9snUXPj/fffZ/369UyePPmSMoiIFDeq2yJ/UfMqefbxxx/TtGmx+tyLiIhHU90Wb6JhA8WYtZZhw4ZRv359GjRo4L639YXuUZ2TW2+9lQ8//BCAadOm0bt3bwBiY2O58cYbadiwIc2aNXNf62/fvn106NCBOnXq8OSTT7qXM3DgQJo2bUp4eDijRo1yP1+jRg1GjRpF48aNadCgAT///DMAhw8fpl27doSHhzNgwACqV6/uvqzLv//9b5o1a0ZkZCQPPvggGRkZ+fTOiYg4Q3VbBNd/hNx+NWnSxIrnCwoKstZaO3v2bNu2bVubnp5uDxw4YKtVq2b37dtnly9fbkuXLm13795tMzIy7A033GBXrVplrbX2pptustdcc41t2LChbdiwoR06dKi11toDBw7Yq6++2n777be2Tp069ujRozYlJcXWrFnTrlu3zlpr7YkTJ2xaWpp97733bM2aNe3x48dtcnKyveqqq+wff/xhrbX26NGj1lpr09PT7U033WR//PFHa6211atXtxMnTrTWWjtlyhR73333WWutHTx4sH3xxRettdYuWrTIAvbw4cM2Li7Odu7c2aamplprrR04cKD94IMPCvy9laILWG/zUO+84Us123uobktxk1PN1rCBYmz16tX06tULX19fKlasyE033URsbCylS5fO8R7V5zv9VLFiRZ577jlat27NnDlzKFu2LFu2bKFy5cpERUUBULp0aff8bdq0oUyZMgDUq1ePXbt2Ua1aNT799FPeeust0tPT2b9/P3FxcURERABw++23A9CkSRO++OIL9zbMmTMHgA4dOnDllVcC8PXXX7Nhwwb3upOTk6lQoUL+v4kiIoVIdVtEY17lAi7lHtVbtmyhXLly7Nu375KWv2PHDsaPH09sbCxXXnkl/fr1y3ZP7dOvyU0eay19+/Zl7NixF80iIuINVLeluNCY12Lsb3/7G7NmzSIjI4PDhw/z7bff0qxZs0ta1rp161i0aBEbN25k/Pjx7Nixg7p167J//35iY2MBiI+Pz7F4nTx5kqCgIMqUKcPBgwdZtGjRRdfbokULPv30UwCWLl3Kn3/+CbiOEMyePZtDhw4BcOzYMXbt2nVJ2yYiUlSobovoyGuxdtttt7FmzRoaNmyIMYaXX36ZSpUquQfVX0jv3r3dl1wJDQ1l4cKF3H///bz33ntUqVKFV155hf79+/PNN98wa9YsHnnkEZKTkwkMDGTZsmUXXG7Dhg1p1KgR1157LdWqVaNFixYX3YZRo0bRq1cvPvroI5o3b06lSpUICQkhNDSUMWPG0L59ezIzMylRogRTpkyhevXqeXuTRESKENVtETD2Ip9GPJPuky1FTUpKCr6+vvj5+bFmzRoGDhzIpk2bnI4lRVBO98n2VqrZUhSpbktu5FSzdeRVPNoff/zBnXfeSWZmJv7+/rz99ttORxIRkRyobsvlUvMqHq1OnTps3LjR6RgiIpJLqttyufSBLSk0ixcvpm7dutSuXZtx48Y5HUdERHKgmi1FlZpXKRQZGRkMHjyYRYsWERcXx4wZM4iLi3M6loiInIdqthRlal6lUKxbt47atWtTq1Yt/P396dmzJ3PnznU6loiInIdqthRlal6lUOzdu5dq1aq5H1etWpW9e/c6mEhERC5ENVuKMjWvIiIiIuIx1LxKoQgLC2P37t3ux3v27CEsLMzBRCIiciGq2VKUqXktQL6+vkRGRrq/Tn9ac9WqVYSHhxMZGUlycjLDhg0jPDycYcOGMXr0aMaPH5/vWSZMmEBSUlK+Lze3oqKi2L59Ozt27CA1NZWZM2fStWtXx/KIiJxNNfsvqtlSlOk6rwUoMDDwvHcN+fjjjxkxYgR9+vQB4K233uLYsWP4+voyevToAskyYcIE+vTpQ6lSpQpk+Rfj5+fH5MmTiYmJISMjg/79+xMeHu5IFhGR81HN/otqthRlOvJayKZPn86nn37K008/Te/evenatSsJCQk0adKEWbNmZZs3Ojqa07d2PHLkCDVq1ABclzAZNmwYUVFRREREMG3aNABWrFhBdHQ0PXr04Nprr6V3795Ya5k4cSL79u2jdevWtG7dulC390w1atTg+htb8vXXXzNy5EjHcoiI5FZxrtnNmzfnb9E3M336dNVsKVJ05LUAJScnExkZ6X48YsQIBgwYwOrVq+ncuTM9evQAIDg42L23n5u9+HfeeYcyZcoQGxtLSkoKLVq0oH379gBs3LiRrVu3UqVKFVq0aMF3333HkCFDePXVV1m+fDmhoaH5vZm5smjRIm6/405sSCXmz4/iq/lzufHGGx3JIiJyPqrZf9m+fTtt2t/C4eRMPvzgAyZNnMBDDz7oSBaRs6l5LUAXOgV1uZYuXcrmzZuZPXs2ACdOnGD79u34+/vTrFkzqlatCkBkZCQ7d+6kZcuW+Z4hr775djUhpctweP+vlK9Zm1Xfr1XzKiJFimr2X7Zs3QrGcOrgDgCWLf9WzasUGWpeizA/Pz8yMzMBOHXqlPt5ay2TJk0iJiYm2/wrVqwgICDA/djX15f09PTCCXsR/xr7AsMef5SKFSuy+/ftTscREcl33lSzb+/Wjdu7dcMYwzvvvEP//v2djiTipjGvRViNGjXYsGEDgHuPHSAmJoapU6eSlpYGwLZt20hMTMxxWSEhIcTHxxdcWBGRYk41W6RwqHktQKfHT53+Gj58eJ5eP3ToUKZOnUqjRo04cuSI+/kBAwZQr149GjduTP369XnwwQcvurf+wAMP0KFDB0cH/4uIFGWq2SKewVhrcz1z06ZN7elPUork1aFDh6hYsSJ5+Z0TyS/GmA3W2qZO5yhMqtlyuTRsQJySU83WkVcRERER8RhqXkVERETEY6h5LUJO35qwfv36dOnShePHj7unbd26lZtvvpm6detSp04dnn/+eZ1+FxFxmOq2SOFT81qEnL7G4E8//UTZsmWZMmUK4PoQQdeuXRk+fDi//PILP/74I99//z1vvPGGw4lFRIo31W2RwqfmtYhq3rw5e/fuBeCTTz7JdkeWUqVKMXnyZMaNG+dkRBEROYPqtkjhUPNaBGVkZPD111/TtWtXwHXqqUmTJtnmufrqq0lISODkyZNORBQRkTOobosUHjWvRcjpawxWqlSJgwcP0q5dO6cjiYhIDlS3RQqfmtci5PTYqV27dmGtdY+dqlevnvuuLaf9/vvvBAcHU7p0aSeiiogIqtsiTlDzWgSVKlWKiRMn8sorr5Cenk7v3r1ZvXo1y5YtA1x7+kOGDOHJJ590OKmIiIDqtkhhUvNaRDVq1IiIiAhmzJhBYGAgc+fOZcyYMdStW5cGDRoQFRXFww8/7HRMERHJorotUjj8nA4gf0lISMj2eP78+e7vGzRowIoVKwo5kYiI5ER1W6Tw6ciriIiIiHgMNa8iIiIi4jHUvIqIiIiIx1DzWsQtXryYunXrUrt2bd2ZRUSkiFPNFil4al6LsIyMDAYPHsyiRYuIi4tjxowZxMXFOR1LRETOQzVbpHCoeS3C1q1bR+3atalVqxb+/v707NmTuXPnOh1LRETOQzVbpHCoeS3C9u7dS7Vq1dyPq1atyt69ex1MJCIiF6KaLVI41LyKiIiIiMdQ81qEhYWFsXv3bvfjPXv2EBYW5mAiERG5ENVskcKh5rUIi4qKYvv27ezYsYPU1FRmzpxJ165dnY4lIiLnoZotUjh0e9gizM/Pj8mTJ9O2bVsA+vfvT3h4uMOpRETkfE7X7Hbt2pGZmcl9992nmi1SAHTktYg7fvwEe/btp1zFKjz00ENOxxERkRxUqVKFk4nJJKak0aVLF6fjiHglNa9F2IKFX9G7992knkomds1qbrq5rdORRETkAg4cOEDzG1tw+MA+Du3bww3Nb+TQoUNOxxLxOmpei7Abrm/GRx99BECzZs14/dXxDicSEZELqVixIu++M939+N13plO+fHkHE4l4JzWvRVhoaCh9+vQBoF27drRp08bhRCIiciHGGHr16uX+vmfPnhhjHE4l4n3UvIqIiIiIx1DzKiIiIiIeQ82riIiIiHgMNa8iIiIi4jHUvIqIiIiIx1DzKiIiIiIeQ82riIiIiHgMNa8iIiIi4jHUvIqIiIiIx1DzKiIiIiIeQ82riIiIiHgMNa8iIiIi4jHUvIqIiIiIx1Dz6iGMMU5HuGyZmZlORxARKRTeULNP86ZtEe+g5tVDVK5c2ekIl2337t1ORxARKRRhYWFOR8g33vD3R7yLmtcibu3atQB0797d4SSXb/DgwU5HEBEpUBkZGQAMGTLE4SSX74svvgCgVatWDicRyU7NaxHXvHlzACpWrOhwkstjrSU2NpYnnnjC6SgiIgVm+PDhAF5R604fNClVqpTDSUSyU/NahE2aNAmABQsWOJzk8pUpUwaAl156yeEknispKYmUlBSnY4jIBcTFxTF+/Hhuu+02jx8nOnLkSAD++9//OpzEc2VkZHDy5EmnY3glNa9FUEJCAm3btmXIkCEMGjSITp06OR3pkm3duhVjDPHx8axbtw5fX1+nI3msxs2ac8+992GtdTqKiJzBWsuzzz5LeHg48Nfpdk/0559/0qhRI1588UWeffZZmjVr5nQkjzXokUdp3Kw58fHxTkfxOn5OB7hcmzdvZvHixeed9s9//rOQ0+SvSpUqUb16dV5++WWno+TZ2e99v379WL58OcuXL3cokec7uH8fc3fvZcyLY3l65FNOxxG5JMeOHWP69OnnnebpNRtcZ5e8oWZfd911lCxZ0iO3pajYEBvLrr0HuK3HXSz5ar4O3uQjj29e5y34imeeeZormnQCX4/fHABKXFmF4GubkwKM+3Kd03EuS/C1N1LiysrM/SWeub949rY4zad+B8pdF81Lrz5F+HXXcvvttzsdSSTP4uLi+Oc//0lIeDR+pcs5HSffXHnD7WCMx9fsgIq1KFWrEQfw/L8/jitxFRX7PMSG/0zk8aFPMvG1V5xO5DVMXk5BNm3a1K5fv74A4+RdRkYGMR27sOkoBLUZ5PHjjEQuJmX/duLnPseq5V/TqFEjp+N4DGPMBmttU6dzFKaiWLMBnn/hRV6Z9gGl73gRnxIlnY4jUqAykuM5MXMY48c8wwP33+90HI+RU832+DGvvr6+zJk9iysT/yDph7lOxxEpcAGV6xAY/SDtO3Zm//79TscRybP/e2oEMS2jSFzyOtbq5iXi3XwDQwjp+n/8Y9hwVq5c6XQcr+DxzStASEgIyxYvJGPTPJJ+i3U6jkiBC7q2JbZuG9rd0pnk5GSn44jkiTGGD997h+ql0kj8/hOn44gUuBLlqhLU4R/censPfv31V6fjeDyvaF4BqlevzsJ5X5K4dCKph3c6HUekwAVdfyf77RXc1fvvuvWueJyAgACWLJxHiZ3fk7hVH+QU7xdYI5ISUXfRJqYjx48fdzqOR/Oa5hVcF/R/c/JE4ue9SEbicafjiBQoYwzB7R5m1Q8/8/So0U7HEcmz8uXLs2zxV5xa/R6n9vzP6TgiBS4o8hYSQq+jc7fupKenOx3HY3lV8wrw97/3YdCAfiQsHIdNT3U6jkiBMn7+BHcZzsQ3p/PJJzOcjiOSZ+Hh4cz6+CMSvnqZ9BOHnI4jUuCCb7qPuP3xDBz8iNNRPJbXNa8AY194nhYRdUhYNkUXdBev5xt0JSFdnuKBQQ/rbjjikTp27Miz/zeC+HljyExJcjqOSIEyPr4EdxzKp/OX8PrEyU7H8Uhe2bz6+Pjw6Sf/ppI9RuK62U7HESlw/hVqUqrtYG7pcit//PGH03FE8uwfjz/G7R1uJmHxK9jMDKfjiBQon4AggruOZOSo0SxZssTpOB7HK5tXgFKlSvGfRQvgf0tJ/OU7p+OIFLhSta/Hp0En2nboREJCgtNxRPLEGMNbb75B3dCSJK76wOk4IgWuxBWVCO44jDt63s3//qcx33nhtc0rQJUqVViycD7Jy98k5YAuTSHeL6jpbRwNDOP2O3rqCgTicUqUKMHCuV8QuH8TiZt1NEq8X8lq9fG/8R7adujI0aNHnY7jMby6eQVo0qQJ709/m/j5Y0mPP+J0HJECZYwh+OaHWL99D08Me9LpOCJ5VrZsWb5euojUtZ9watdmp+OIFLjgBm1JCWtKh863kpqqD5rnhtc3rwA9enTnycceJmH+WDLTTjkdR6RAGd8SBHX6J+/8exbvvPuu03FE8uyaa67hi89mkbBoPGnH9jodR6TABf3tHn4/abl3wAP6oHkuFIvmFeDp/xtJ+xZNSFw6UbcjFK/nG1iakK4jGfKPYXz77beXvJwXXniB8PBwIiIiiIyMzNerGRw/fpw33ngj35Yn3qVNmza8PHYMCfNfIOOUxnCLdzPGh6CYx1n4zXeMe/lfl7yc4lKzi03zaozho/ff5aqSKSR+r+thivcrUa4awTGP0fW27vz22295fv2aNWtYsGABP/zwA5s3b2bZsmVUq1Yt3/IVpUIoRdOggQO5+/ZbSfzqX9gMXdBdvJuPfyDBXZ/ihZfGM3fu3Dy/vjjV7GLTvMJftyP027GaBN2OUIqBwJqNKdH0DtrGdOTEiRN5eu3+/fsJDQ0lICAAgNDQUKpUqUKNGjV48sknadCgAc2aNXPfp/vw4cN0796dqKgooqKi+O4711U+Ro8eTf/+/YmOjqZWrVpMnDgRgOHDh/Pbb78RGRnJsGHD8nGrxZtMnvgaEVeVI3HFdJ1OFa/nV7oCIZ2H06fvvWzenLcx38WpZher5hWgQoUK/GfRQlJWvUvKXl2aQrxfUKNOnCxbly639cjT7Qjbt2/P7t27ueaaaxg0aBArV650TytTpgxbtmzh4Ycf5rHHHgPg0Ucf5fHHHyc2NpbPP/+cAQMGuOf/+eefWbJkCevWrePZZ58lLS2NcePGcfXVV7Np0yb+9a9LP00m3s3X15d5X3xG6RO/krRxgdNxRApcQJW6lGw1gHa3dOLgwYO5fl1xqtnFrnkFaNCgATP+/SHxC1/S7QilWAiOHsBPe48zeMijuX9NcDAbNmzgrbfeonz58tx11128//77APTq1cv975o1awBYtmwZDz/8MJGRkXTt2pWTJ0+6rzfbqVMnAgICCA0NpUKFCnkqyCKlS5fm6yVfkf7DFyT/vsHpOCIFLqjeTWRc3Yr2Hbtw6lTuPmhenGq2n9MBnNK5c2dGPTWcFya8QOk7xuITUMrpSCIFxvj4EnzLUGbO+icNwuvx8ODBuXqdr68v0dHRREdH06BBAz74wHXxeGPMX8vO+j4zM5O1a9dSsmTJc5Zz+jTW6WXm5QiwCEDNmjWZ/+UXdOh8K749xuAfepXTkUQKVFDzXuxZ9Ap3/70vn386M1vdvZDiUrOL5ZHX04Y+8Q+6truJhMWv6naE4vV8SgYT3HUkw/9vFEuXLr3o/L/88gvbt293P960aRPVq1cHYNasWe5/mzdvDrhOWU2aNCnb/DkJCQkhPj4+r5shxVjLli2ZPOFVEua9QEZS3sZwi3gaYwxB7R7hm3VbGPXscxedvzjV7GLdvBpjeOetN6lTtgSJqz50Oo5IgStxZRWCbxlKj7vu5ueff85x3oSEBPr27Uu9evWIiIggLi6O0aNHA/Dnn38SERHB66+/zmuvvQbAxIkTWb9+PREREdSrV48333wzx+WXK1eOFi1aUL9+fccH/4vn6NevL/f3vZuEhS9h09OcjiNSoHxKBBDcZQSvTZ7KzJmzcpy3ONVsk5dPbzZt2tSuX7++AOM449ixY0Q0bkpKvS4ERbR3Oo5IgUvcvJSAuPn8uCGWcuXK5em1NWrUYP369YSGhhZQuoJhjNlgrW3qdI7C5K01OzMzk45duhG7P4Xgdo/k6nSqiCdLPfgbJ+eMZvl/lhAVFZWn13pjzS7QMa/BwcHuwb+F6ZlnnqFVq1a0bds2V/OXLVuWr5csolnzFvheUYmSV0W4p+2Z2h8f/0AwBt+gKwnt9A98g6/863kf18Hrsu0GUbLqdaQe3sWxZdPISDgKNpOg8Jspc2PPbMV133uPUKJsVcrf+s/83XCRXAqKaE/C8b107Hobq5Yvw9/f3+lIUgR4Ss328fFh9qxPaHL9jRyN/YLgZt3d0/KzZmck/snRRRNJP3kYMjPwLVOBinc8WyDvgUhO/CteTambB3FL565s2hBL1apVnY7kKK/8wNZzz118bMjZ6taty+xZM+h2R0987xpHiSuruKdV7PUivqXK8OfKDzix9lPKtn0w2/OnZaalcPiL5ynbfhCBNRuTmXaKw3PGkrBxISGNOwOQdmQ3ZGaSsmcrmamn8PE/d6C0SGEIankP2xeM474HHuLD997J9dGrnTt3FmwwKXYupWYHBwfz9ZKvaNgkiqQrwyhV5wb3tPyq2cdXfUzJGpGUbnorAKmHdlzmlopculLXNCfh+D7adujEhv9+T1BQUK5e5401u1DGvK5YsYLo6Gh69OjBtddeS+/evd0Xm46NjeXGG2+kYcOGNGvWjPj4eE6dOsW9995LgwYNaNSoEcuXu24o8P7779OtWzfatWtHjRo1mDx5Mq+++iqNGjXihhtu4NixYwD069eP2bNnA67D5aNGjaJx48Y0aNDAPc7v8OHDtGvXjvDwcAYMGED16tVp1KgRL7/4PAnzXiDzPLcjLFmtPml/7rvgdibGrSQg7DoCazYGwKdEScq2e4gTa2f/Nc//VhJUvzUlazQm+de1+fDuilwa4+NLcMzjzF+2in+NfyVXrzl48CB33303tWrVokmTJjRv3pw5c+a4pz/22GOEhYWRmZn9Fswffvgh9evXd/+fHj9+fL5ui+QvT6nZJUuWZPGCeSR9PYXUg7+fsx2XW7MzEo/hG/LXqVb/CjUv960VuSxBUbdzuEQF7ujZ+5w6ez7eWrML7QNbGzduZMKECcTFxfH777/z3XffkZqayl133cXrr7/Ojz/+yLJlywgMDGTKlCkYY9iyZQszZsygb9++7uuc/fTTT3zxxRfExsYycuRISpUqxcaNG2nevDkffnj+D12Fhobyww8/MHDgQPcP4Nlnn+Xmm29m69at9OjRgz/++AOAwYMG0eu2LiR8Nf6cKxAk/7YO//I13I8PzniKfe89wv4P/wFA2pFd+Feqne01Ja6sjE1NJjMlCYDEn1cRdF0rguq1IjFuJSJO8gkoRXCXp3hu7EvMnz8/x3mttXTr1o1WrVrx+++/s2HDBmbOnMmePXsA1zjEOXPmUK1atWwXx160aBETJkxg6dKlbNmyhbVr11KmTJkLrUaKCE+p2VFRUUx/cyrx818kI+HPbMu53Jod0qgTRxdN5MCMEZz4fhbp8Ufz5b0VuVTGGILbDGbN1t95cvhTOc7rzTW70JrXZs2aUbVqVXx8fIiMjGTnzp388ssvVK5c2T34uHTp0vj5+bF69Wr69OkDwLXXXkv16tXZtm0bAK1btyYkJITy5ctTpkwZunTpArhuPHChQ+O33347AE2aNHHPs3r1anr27AlAhw4duPLKK93zT574GvWrXkHiiunAXwUvMyWJMjfc4Z6vYq8XqXLvJCrf82qu3oOU/dvxDSyNX+kKlKzekNSDv5ORXDQuOyHFl1+ZCoR0Gs7df+/Lli1bLjjfN998g7+/Pw899JD7uerVq/PII48ArqN14eHhDBw4kBkzZrjnGTt2LOPHj6dKFddQnICAAO6///4C2hrJL55Us3v2vIvHBj9IwoIXAZtvNTuwVhPCHpxOSMMY0o7tYf/7j+oSXeI441eC4M7Dmfb+R3zwwYWvlOTNNbvQmtf8uuDtmcvx8fFxP/bx8bngMk/Pk9v1+vn5MX/ObEL+3IZNTXIXvNDOT+BTMviCrysRehWpB37N9lza8QMY/0B8AkqR+L+VpB3dw56p/dk77X4yU5NI+uW7i+YRKWgBYddSstV9tLulE4cOnf+uc1u3bqVx48YXXMaMGTPo1asXt912GwsXLiQtzXUZo59++okmTZoUSG4pOJ5UswGeGz2K6Kb1scknqdDzhXyp2QC+gSEE1YsmtPMTBFSuw6ndP+Uqj0hB8i1VhpAuIxn86ON89935+whvrtmOXue1bt267N+/n9jYWADi4+NJT0/nb3/7Gx9//DEA27Zt448//qBu3br5uu4WLVrw6aefArB06VL+/DP76aYyZcrw9ZKvsKlJnPpjc66WGVQvmpQ9cSTv3AS4Pgzw57JplL6+O9ZmkvTzair3n0zVge9SdeC7VLj9aRL/922+bpfIpQqqF01y6eqMeXFsruYfPHgwDRs2JCoqitTUVL766iu6detG6dKluf7661myZEkBJ5bCVpRrtjGGTz76AF9jSdrwZa6WmVPNBkje9SOZaa7hD5kpSaT9uR+/0uXzZ4NELpN/+er4N+3OA4OH5Gp+b6rZjl5twN/fn1mzZvHII4+QnJxMYGAgy5YtY9CgQQwcOJAGDRrg5+fH+++/n23vPT+MGjWKXr168dFHH9G8eXMqVapESEhItnlq1apFxQoVOPnNNPxDa1AitFqOy/QpEUD57k9z7D9vcuw/UyEzk6Dw1oQ07kzK7p/wDS6HX8hf19QMqBZO2tE/SE84hl9w2XzdPpG8St6xEXPwZx595P3zTg8PD+fzzz93P54yZQpHjhyhadOmLFmyhOPHj9OgQQMAkpKSCAwMpHPnzoSHh7NhwwZuvvnmwtgMKUBFvWYHBgZSsUIFkn9fTWL5qwm67m85LjOnmg2QeuBX/vzPm+DjCzaTkIbtCah8Tb5ul8ilSj95mLSNc3nx3bfOO92ba3axvUlBSkoKvr6++Pn5sWbNGgYOHHjBW6O9++57PDb8acr0fDnbZVZEvEXa0d2c+GwkX82bQ6tWrc47j7WWG264gX79+jFw4EAA/vjjD1q1asWNN95Ily5d6NWrFwCJiYnUrFmTnTt3smLFCp5++mkWLlxIpUqVSE1N5cMPP2TAgAGFtn2gmxR4urzU7B9//JEWN7UmpOv/EVAlf48AixQFmanJnPzsKZ4cfB9PDT//NeO9uWZ75XVec+OPP/7gzjvvJDMzE39/f95+++0Lztu//72s/v575q9+n6D2jxZiSpGCl5F8kvh5LzDhlZcv2LiC67Tsl19+yeOPP87LL79M+fLlCQoK4tlnn+Xxxx/PdmvBoKAgWrZsyfz587nrrrs4ePAgbdu2xVqLMYb+/fsXxqaJF8lLzW7YsCFvTJzAw0OfIqD/tEJMKVLwrM0kcckEOkY3Z8Q/n7zgfN5cs4vtkde82LdvHxGNm+LT/F5K1b3R6Tgi+cZmpHFyzmju6dqG11/N3bVePZWOvBYfaWlptLq5HdvTyxF8U9H6oytyuRJWfUDNzH18t/Ibr747Yk4129EPbBWW/v37U6FCBerXr5/n1yYlJdG2Qye4rr0aV/Eq1loSv5lGVJ2qvPqvl52OI5LNpdZtay0DHhzIL0dTCPpb3wJKJ+KMxJ++JmB3LIvmf+nVjevFFIvmtV+/fixevDjPr8vMzOTOu/tw0CeUoGY9CiCZiHMSN8ylbPIevvhsJr6+vk7HEcnmUuv2K6++xtylKwiO+QfGR7/X4j1O7dlKyncfsGzJV4SGhl78BV6sWDSvrVq1omzZvH+af8TIp/lu83aC2w7O9X3fRTxB0q/ryNw8n2WLFxIcfOHrYIo45VLq9sKFCxn9wjiCu4x0X6dVxBukHT9AwsJ/8emMj7nuuuucjuO4YtG8Xop///tj3pj+PsGdhmP8SjgdRyTfpB7aQdKyyXw170uuuuqqS17O4sWLqVu3LrVr12bcuHHnTN+1axdt2rQhIiKC6Oho9y0JRQrCTz/9RM8+9xDS6Un8ylRwOo5IvslMSSRh3guMGf0MHTp0uOTleFPNVvN6HmvWrOGhh4cQ0vUpfIOucDqOSL7JSPyT+PkvMm3KJG644YZLX05GBoMHD2bRokXExcUxY8YM4uLiss0zdOhQ7rnnHjZv3swzzzzDiBEjLje+yHkdOnSIth06UvJv9xIQpqNS4j1sZgYJi17hjs7teezRRy55Od5Ws9W8nmXXrl106tqNUu0ewb98DafjiOQbm55KwoJxPPzAffTuffdlLWvdunXUrl2bWrVq4e/vT8+ePZk7d262eeLi4twXuW7duvU500XyQ0pKCh0630p6jRYE1WvtdByRfJW48l3qVQrmzTcmX9ZyvK1mq3k9Q3x8PG1iOuIb2ZVSV0c5HUck31hrSVg2hb81upYXnn/2spe3d+9eqlX7645zVatWZe/evdnmadiwIV988QUAc+bMIT4+nqNHj172ukVOs9by93792ZXsT1CLy9shEylqEjctIujwVhZ8+Tl+fpd3WX5vq9nFonnt1asXzZs355dffqFq1aq8884758yTkZHBbT3u4nhwdUo1vtWBlCIFJ/G/n1HJHmPmvz/Ex6dw/tuPHz+elStX0qhRI1auXElYWJiuaiC5lpu6PebFsSz9bgNB7YdgTLH4cybFRPLOTaTFzuLrpYu44oorCmWdnlSzi8UdtmbMmHHReR57Yhg//H6AkG6jdGUB8SqJP6/G/LKMZT+sp1Sp/PkEdlhYGLt373Y/3rNnD2FhYdnmqVKlinsvPiEhgc8//7zQirB4vovV7S+++IKXX51ImZ4v4VOiZCGlEil4aUf3kLj4VRZ8+Tm1a9fOl2V6W83Wrirw1ttv88GMzwjq+CTGt1j081JMpOzfTvKKaSz9agGVK1fOt+VGRUWxfft2duzYQWpqKjNnzqRr167Z5jly5AiZmZkAjB07tsjdXlA818aNG7mn/wCCu4zAL6R4X+9SvEtGcjzx81/g1X+NIzo6Ot+W6201u9g3rytXruQfw4YT0vX/8A0McTqOSL5Jjz9C/IKxfPjudBo1apSvy/bz82Py5MnExMRw3XXXceeddxIeHs4zzzzDvHnzAFixYgV169blmmuu4eDBg4wcOTJfM0jxtG/fPtp37Exg9IMEVMqfo1IiRYHNSCfxq5fp27MHD9x/f74u29tqtrHW5npmb7tP9q+//krT65sT0O4xAmtEOh1HJN9kpp3i5GdP8cSDfXl65FNOxykScrpPtrfytpqdlJREVPOWHLyyAcE33Ol0HJF8Y60l8es3iCwHS76aX2THmhamnGp2sT1Hfvz4cdrEdKRE1F1qXMWrWJtJ4pLXiWkZxf89VXSv0yeSF5mZmdzV++8cMGUJvv4Op+OI5KukH+ZxRcIu5ixfq8Y1F4pl85qenk7nbt1JDK1HcOQtTscRyVeJ339C9VJpfPjeO/rwoXiNkU+PYvWP2yh9+3P6vRavkvxbLBmb5vL1+nWEhGj4Ym4Uy+b1oUEPE3cggZCujzodRSRfJW5dTomd37Pkh/UEBAQ4HUckX3z88SdMfutdruj5MsbP3+k4Ivkm9fBOEv8zif8sWkD16tWdjuMxil3z+vrEyXy2YCll7hqH8dGhefEep/b8j1Or32Pt6m8pX76803FE8sXatWt5cPAjlO7+nG7XLV4lI/E48fNeZOqk12nevLnTcTxKsbrawJIlSxg56lmCu47EJyDI6Tgi+Sb9xEESFr7ErI8/Ijw83Ok4Ivnijz/+oGOXWynV9mHdrlu8yunbdQ8a0I+//72P03E8TrFpXv/3v/9xR8+7Ce44lBJXVHI6jki+yUxJIn7eC4z+vxF07NjR6Tgi+eL07bp9GnalVO1mTscRyTeu23W/QYuGdRj7wvNOx/FIxWLYwJEjR2gb05GAFn0pWa2+03FE8o3NzCBh8SvcFtOaJ/7xuNNxRPJFRkYG3e/sxbFS1Qhuott1i3dJXDebiplH+PSTuYV2u25v4/XNa2pqKrd06UZKtSiC67dxOo5Ivkpc9QF1Q0vy9rSp+gS2eI1/DHuS9dv3EnKbbtct3iXpl+/hf0v5z4bYfLtdd3Hk1c2rtZZ+993P7ychuPPfnY4jkq8SNy8hcP8mFv4QS4kSJZyOI5Ivpr/zDu99/Cller6M8dXvtXiPlAO/krR8Kiu//g9hYWFOx/FoXt28jn3pZb5a/j2l7xyLMTo0L97j1K7NpK79hO/Wfk/ZsmWdjiOSL1auXMljTzxJmTtewDewtNNxRPJNevxR4ueP5f3pb9OkSROn43g8r+3o5s6dy4svv+K6soB/oNNxRPJN2rG9xC8azxefzeKaa65xOo5Ivvjtt9+49fYeBHX4ByXKVXM6jki+yUw7RcKCsQx77GF69OjudByv4JVHXn/88Uf69OtPSJeR+JXW9S6LM5uRxr53Bl/y69P/3Jfn1/iVDSOkcedLXueF+PgHElj7ehLmv8C/xo6hTRuN4RbvcOLECdq0v4USUXfqdt3C4Tkvknp45yW//lLq9pVtH7zk9eUkOLw1iV+/QbvmjXjm/0YWyDqKI69rXg8cOEC7WzpRstUAAqrUdTqOOMxmZl5SIbsc6cf20uPq/D+p8cmM90iK/Zy+d3Rj0MCB+b58ESecvl13Quh1BEfqUm8Cp/bGkZl4vFDXWRA1e8WKFez6YQH1alXlo/ff1YcP85FXNa/Jycm079iFzDrRBF/Xyuk4UgT4lAig+j8X5MuyrLWkH9tDypalJMUtp9ddd/L6a68QHBycL8u/mGXfLOeqq65i8sTXCmV9IoVh0CND2LrvJCG3DnE6ihQR1R7+d74tKyPxOElbvyHtpyU0uO4aPv7gXWrUqJFvy8/JvQMe4MSiRSxZOI+SJUsWyjqLC68Z82qtpfc9/dibEUzQDT2djiNeyBhDiXLVCI6+j9D+bzJ3w07qXBtObGxsoax/1fJlzPviM3x9dVtj8Q6TJk9h1txFBHccqtt1S4HwDbqCkGa3c2W/N/jFtwYNIhvz/gcfFMq6J772Chv+u4YKFSoUyvqKE69pXp8Z/SzLY7cQ1PYRHZqXAudTMpjg9kPIaNaHm9t3YPXq1QW+zrCwMEqX1iewxTssXbqUEU+PIriLbtctBc/4+BLcrDsh3cfwyNARvDbh9QJfZ0hICFWrVi3w9RRHXtG8zpw5iwlTphHUeQQ+JQKcjiPFSKlrbqRU+8fp2KUb27dvdzqOiEf4+eef6dHzboI7DqPElZWdjiPFiH/56pTp8QLPjBnH7NmznY4jl8jjx7zGxsZy74D7Kd1xKMYYMhL/dE9LPbwLrHUw3fmVCK123lNkGUknMSUCvKIBT48/im/QFV5xKtCUKOn6usAR/cCajcho2p2effoSu2a1bvcnkoOjR4/SJuYW/CI6U6JsWLaanX78IJmpyQ6mOz+/MhXwCTj3bkiZqcmQmYFPycIZ916QMpJO4OMfiPHzdzrKZTO+JTABpS54fXe/MhUI7jiMAQ8OpFWrVjqt74E8vnn9ZsVK/P0DOLVsEqfOmnbyzyOOZMqN0leGnvPc2XmDQsrg6+eZd5g5e1vOt72ewFpLyqkkSpYuh1/99pSq3/a8F08PatyZHbPX8sbUqTw8+NIvzSXi7X7++WcSExKwmxdycvPCbNM8vWYbHx9CynjmTUPO3paQK8p67M190lJPYY0vQfXbENAghhJlz72bVUCVuqRd25oHBw9hzmczHUgpl8PYPByZbNq0qV2/fn0BxpFjx45x++23s3LlSgB2797tsWNm1q1bx/XXXw+Av78/p06d8sjxyNZa1q1bx6uvT2bBV18R0nEoJa+KOGe+lL3/w+fbN9i941cdfS2CjDEbrLVNnc5RmFSzC156ejrDhg1jwoQJAHz55Zfceuutzoa6RPv376dTp05s3LgRgCNHjlCuXDmHU12aHTt2MPmNN3nr7ekERN1BcJMu58yTmZrMkbfvY9v/tnrs31lvllPN1l/YIqZs2bKsWLGClJQUAKpVq0ZCQoLDqS5Ns2bNsNayceNGUlNTufLKK52OdEmMMVx//fXM+uQj5n3+KQlf/YuU/dvOmc+/yrUk2xIsW7bMgZQi4gQ/Pz9ee+01rLW0aNGCbt26sWbNGqdjXZLKlSvzww8/EB8fD0BoaChpaWkOp7o0NWvW5JV/vcSWTRsouW0piT8uPmceH/9AStWLZsrUNx1IKJdDzWsR5e/vT2pqKgDt2rVzOM3liYyMZN26dZw4cYJ169Y5HeeytGnThmlTJpHy7fRzphljsLVu5Mt5+XNdWRHxLKevOnLjjTc6nOTyBAcHuw+a9OvXz9kwl6lGjRosX7aEhNUfknHq3ANBfnVa8vmX8x1IJpdDzWsRVqJECR599FHWrl1LXoZ3FEVRUVEA7mEEnqxnz56USDlJyv5zry7gX6k23/+3cK77KiJFz88//wzAtm3nnp3xJEFBQdx666188sknTke5bHXq1CEmJoakn74+Z5p/xVrs/G2b+2CReAY1r0Xc+PHjAViwwPOP5r399ttOR8gXvr6+PHj/faRuW3XONP8KNdn281YHUolIUVC3ruu25J5+xBLg3/923elqy5YtDie5fI8MfBDz+/fnPO9ToiRBZSt6/M5GcaPmtYjz83NdEGLz5s0OJ7l8bdq0cTpCvql9dS18U06e87xPQBApSUkOJBKRosQbavbpW1//9NNPDie5fDVr1iQt4dh5p/kFBrvH+YpnUPPqIZKTi961D/MqKMh77qLj7++PsRnnTvDxwdrMwg8kIkVKkhftxHrD3x9/f38yM85TswGMD5mZqtueRM2riIiIiHgMNa8iIiIi4jHUvIqIiIiIx1DzKiIiIiIeQ82riIiIiHgMNa8iIiIi4jHUvIqIiIiIx1DzKiIiIiIeQ82riIiIiHgMNa8iIiIi4jHUvIqIiIiIx1DzKiIiIiIeQ82riIiIiHgMNa8iIiIi4jHUvIrkQWZmJtdcF07/AQ9gA0qfdx7/oNJcUa48T44YWcjpRETkbL3+3pd6DRriW6rM+WcoWZqYWzrR6ua2hRtMLpmaV5E8MMZwb7++ZAI+tVucbw5KRbTH19eXnnd0L+x4IiJyln597iYjIwPf69qcd7q55ibSMzK4/95+hRtMLpmf0wFEPIkxhhH/fJJ619alT997yWg1gKB6NwFg09NIWDaFSplHWLbpB6pUqeJwWhERiYmJYf1/19C2Q0cSEo4Q9Ld7MMYHay1JP8wjY9Ncli9bSvPmzZ2OKrmkI68il+DWW29l9crlmPWfkPD9J2Qk/snJL56meY0Q1q/9To2riEgRct1117H5h/XUtPtJWDCOjFMJJH79BmX+WM3G9evUuHoYNa8il6hhw4Zs/mE9YYnbOfj2AzxwV1fmfTGbUqVKOR1NRETOUq5cOb5b+Q0dmtThwJv9aFjWsnH9f6levbrT0SSP1LxKsebr60tkZKT7a9y4cQCsWrWK8PBwIiMjSU5OZtiwYYSHhzNs2DBGjx7N+PHjAahUqRL//e5bYv+7hpfGvoCPz6X/l5owYQJJSUn5sl0iIt7ocmu2v78///7gPdasXsXSRQsICQm55Cyq2c7RmFcp1gIDA9m0adM5z3/88ceMGDGCPn36APDWW29x7NgxfH19GT16dLZ5S5YsScOGDS87y4QJE+jTp4+O3IqIXEB+1GxjDE2bNr3sLKrZztGRV5GzTJ8+nU8//ZSnn36a3r1707VrVxISEmjSpAmzZs3KNm90dDTr168H4MiRI9SoUQOAjIwMhg0bRlRUFBEREUybNg2AFStWEB0dTY8ePbj22mvp3bs31lomTpzIvn37aN26Na1bty7U7RUR8WSq2cWPjrxKsZacnExkZKT78YgRIxgwYACrV6+mc+fO9OjRA4Dg4GD33v7Ze/Hn884771CmTBliY2NJSUmhRYsWtG/fHoCNGzeydetWqlSpQosWLfjuu+8YMmQIr776KsuXLyc0NDS/N1NExCuoZguoeZVi7kKnoC7X0qVL2bx5M7NnzwbgxIkTbN++HX9/f5o1a0bVqlUBiIyMZOfOnbRs2TLfM4iIeBvVbAE1ryKXxc/Pj8zMTABOnTrlft5ay6RJk4iJick2/4oVKwgICHA/9vX1JT09vXDCiogUc6rZ3kFjXkUuQ40aNdiwYQOAe48dXBfFnjp1KmlpaQBs27aNxMTEHJcVEhJCfHx8wYUVESnmVLO9g5pXKdZOj586/TV8+PA8vX7o0KFMnTqVRo0aceTIEffzAwYMoF69ejRu3Jj69evz4IMPXnRv/YEHHqBDhw4a/C8icgGq2QJgrLW5nrlp06b29Kf0pPAYYxg5ciRjxoxxOsplOXToEBUrViQvv3Mi+cUYs8Fae/nXx/EgqtnOMMZgjHGfnvZkxhjeeecd+vfv73QUKWZyqtk68ioiIiIiHkPNq4iIiIh4DDWvIhdx8OBB7r77bmrVqkWTJk1o3rw5c+bMcU9/7LHHCAsLO+cU4Ycffkj9+vVp0KABjRo1ct+eUERECo5qtvdT8yqSA2st3bp1o1WrVvz+++9s2LCBmTNnsmfPHgAyMzOZM2cO1apVY+XKle7XLVq0iAkTJrB06VK2bNnC2rVrKVOmjFObISJSLKhmFw9qXkVy8M033+Dv789DDz3kfq569eo88sgjgOsagOHh4QwcOJAZM2a45xk7dizjx4+nSpUqAAQEBHD//fcXbngRkWJGNbt4UPMqkoOtW7fSuHHjC06fMWMGvXr14rbbbmPhwoXuawT+9NNPNGnSpLBiiogIqtnFhZpXkTwYPHgwDRs2JCoqitTUVL766iu6detG6dKluf7661myZInTEUVEJItqtnfS7WFFchAeHs7nn3/ufjxlyhSOHDlC06ZNWbJkCcePH6dBgwYAJCUlERgYSOfOnQkPD2fDhg3cfPPNTkUXESl2VLOLBx15FcnBzTffzKlTp5g6dar7uaSkJMB1+mn69Ons3LmTnTt3smPHDv7zn/+QlJTEiBEjGDZsGAcOHAAgNTWV6dOnO7INIiLFhWp28aDmVSQHxhi+/PJLVq5cSc2aNWnWrBl9+/bl2WefZfHixXTq1Mk9b1BQEC1btmT+/Pl07NiRhx9+mLZt2xIeHk7jxo05efKkg1siIuL9VLOLB90e1gPo9rAil0+3h5XCotvDilw+3R5WRERERLyCmlcRERER8RhqXkVERETEY6h5lUKxfv163n//fQCmT5/O1q1bnQ10iRYvXkzdunWpXbs248aNO2f6t99+S+PGjfHz82P27NkOJBQRuXx79uxxf9p+zpw5fP311w4nujSq2d5JzasUimnvvMc///lPAO6//35mfuZ5RSIjI4PBgwezaNEi4uLimDFjBnFxcdnmueqqq3j//fe5++67HUopInL5vl292n171AULFvDsC+c2fkWdarb30k0KpFC8OXkiJfz8eHPKJJ4Z/Syjnv4/pyPl2bp166hduza1atUCoGfPnsydO5d69eq556lRowYAPj7aLxQRz3V3z56kpaQy4P77aRdzC5/N/NjpSHmmmu299NOSAuHr60tkZCT169enS5cuxMfH88ak1zly5Ah3dL+dNm3aULduXerUqcPzzz/vEZfP2rt3L9WqVXM/rlq1Knv37nUwkYhI/ji7Zh8/fpy+fe/h8KGDvPTi83Tp0kU1W4oMNa9SIAIDA9m0aRM//fQTZcuWZcqUKe7nu3btyvDhw/nll1/48ccf+f7773njjTccTiwiUnxdqGYHBATQrVs31WwpUtS8SoFr3ry5e2/3k08+oUWLFrRv3x6AUqVKMXny5PMOpC9qwsLC2L17t/vxnj17CAsLczCRiEj+U82Wok7NqxSojIwMvv76a7p27QrA1q1badKkSbZ5rr76ahISEor8rfiioqLYvn07O3bsIDU1lZkzZ7q3S0TEG6hmiydQ8yoFIjk5mcjISCpVqsTBgwdp166d05Eum5+fH5MnTyYmJobrrruOO++8k/DwcJ555hnmzZsHQGxsLFWrVuWzzz7jwQcfJDw83OHUIiIXp5qtmu1JdLUBKRCnx08lJSURExPDlClTGDJkCPXq1ePbb7/NNu/vv/9OcHAwpUuXdiht7nXs2JGOHTtme+65555zfx8VFcWePXsKO5aIyGVRzRZPoiOvUqBKlSrFxIkTeeWVV0hPT6d3796sXr2aZcuWAa69/SFDhvDkk086nFRERFSzxROoeZUC16hRIyIiIpgxYwaBgYHMnTuXMWPGULduXRo0aEBUVBQPP/yw0zFFRATVbCn6NGxACkRCQkK2x/Pnz3d/36BBA1asWFHIiURE5EJUs8WT6MiriIiIiHgMNa8iIiIi4jHUvBZjBw4coGfPnlx99dU0adKEjh07sm3bNsfyTJgwgaSkpDy/7v3332ffvn0FkEhEpGhR3RZR81psWWu57bbbiI6O5rfffmPDhg2MHTuWgwcPOpYppyKYkZFxwdepCIpIcaC6LeKi5rWYWr58OSVKlOChhx5yP9ewYUNatmzJsGHDqF+/Pg0aNGDWrFkArFixgujoaHr06MG1115L7969sdYCros833jjjTRs2JBmzZoRHx9PRkYGw4YNIyoqioiICKZNm5bjciZOnMi+ffto3bo1rVu3BiA4OJgnnniChg0bsmbNGp577jmioqKoX78+DzzwANZaZs+ezfr16+nduzeRkZEkJyezYcMGbrrpJpo0aUJMTAz79+8v5HdXRCT/qW6LZLHW5vqrSZMmVgofYEeOHJmvy3z99dftY489ds7zs2fPtm3btrXp6en2wIEDtlq1anbfvn12+fLltnTp0nb37t02IyPD3nDDDXbVqlU2JSXF1qxZ065bt85aa+2JEydsWlqanTZtmn3++eettdaeOnXKNmnSxP7+++8XXI611lavXt0ePnw423bPmjXL/fjo0aPu7/v06WPnzZtnrbX2pptusrGxsdZaa1NTU23z5s3toUOHrLXWzpw509577735+daJhwLW2zzUO2/4Us12BmCNMfm+XNVtKU5yqtm6VJZks3r1anr16oWvry8VK1bkpptuIjY2ltKlS9OsWTOqVq0KQGRkJDt37qRMmTJUrlyZqKgoAPcdV5YuXcrmzZuZPXs2ACdOnGD79u34+/ufdzktW7Y8J4uvry/du3d3P16+fDkvv/wySUlJHDt2jPDwcLp06ZLtNb/88gs//fST+9aGGRkZVK5cOZ/fJRGRokN1W4obNa/FVHh4uLtA5VZAQID7e19fX9LT0y84r7WWSZMmERMTk+35FStW5Ho5JUuWxNfXF4BTp04xaNAg1q9fT7Vq1Rg9ejSnTp0673rDw8NZs2ZNnrZNRKSoU90WcdGY12Lq5ptvJiUlhbfeesv93ObNm7niiiuYNWsWGRkZHD58mG+//ZZmzZpdcDl169Zl//79xMbGAhAfH096ejoxMTFMnTqVtLQ0ALZt20ZiYmKOmUJCQoiPjz/vtNMFLzQ0lISEhGwF/MzX1a1bl8OHD7uLYFpaGlu3br3Y2yEiUuSpbou46MhrMWWMYc6cOTz22GO89NJLlCxZkho1ajBhwgQSEhJo2LAhxhhefvllKlWqxM8//3ze5fj7+zNr1iweeeQRkpOTCQwMZNmyZQwYMICdO3fSuHFjrLWUL1+eL7/8MsdMDzzwAB06dKBKlSosX74827QrrriC+++/n/r161OpUiX36S6Afv368dBDDxEYGMiaNWuYPXs2Q4YM4cSJE6Snp/PYY48RHh5+2e+ZiIiTVLdFXIzN+uRhbjRt2tSuX7++AOPI+RhjGDlyJGPGjHE6iojHMsZssNY2dTpHYVLNdoYxBmMMmZmZTkcR8Vg51WwNGxARERERj6HmVUREREQ8hprXIiwzM9M9WD45Ofm8n9IUEZGi43TNttZe0m1TReTi1LwWYZ9//gXBwcEAvPrqq4RHRDob6DItXryYunXrUrt2bcaNG+d0HBGRfHXgwAFCQkLcj4ODgzl06JCDiS6ParYUVWpei7DOnTvRrfsdhFSqQcgV5XhzyiSnI12yjIwMBg8ezKJFi4iLi2PGjBnExcU5HUtEJN9UrFiRMS+OJejK8gRdGcpLL/+L8uXLOx3rkqhmS1Gm5tUBL7zwAuHh4URERBAZGcl///vf884XGBjIF5/NYsrLzxG79jv33Udycvz4cd544438jnzZ1q1bR+3atalVqxb+/v707NmTuXPnOh1LROSicluzjTE8NfyffPnpJ8ybPYthQ5/AGJPjslWzRfJO13ktZGvWrGHBggX88MMPBAQEcOTIEVJTUy84vzGGv//977le/ulCOGjQoPyIm2/27t1LtWrV3I+rVq16wT8AIiJFRV5rNkDbtm1zvXzVbJG805HXQrZ//35CQ0Pdt9oLDQ2lSpUq1KhRgyeffJIGDRrQrFkzfv31VwAOHz5M9+7diYqKIioqiu+++w6A0aNH079/f6Kjo6lVqxYTJ04EYPjw4fz2229ERkYybNgwZzZSRMRLqGaLFD1qXgtZ+/bt2b17N9dccw2DBg1i5cqV7mllypRhy5YtPPzwwzz22GMAPProozz++OPExsby+eefM2DAAPf8P//8M0uWLGHdunU8++yzpKWlMW7cOK6++mo2bdrEv/71r8LevAsKCwtj9+7d7sd79uwhLCzMwUQiIhenmu2imi1FiYYNFLLg4GA2bNjAqlWrWL58OXfddZf7U5y9evVy//v4448DsGzZsmyD5E+ePElCQgIAnTp1IiAggICAACpUqMDBgwcLeWtyLyoqiu3bt7Njxw7CwsKYOXMmn3zyidOxRERypJqtmi1Fj5pXB/j6+hIdHU10dDQNGjTggw8+AMg2sP/095mZmaxdu5aSJUues5zTp7FOLzM9Pb2Ak186Pz8/Jk+eTExMDBkZGfTv31/3rRYRj6CarZotRYuGDRSyX375he3bt7sfb9q0ierVqwMwa9Ys97/NmzcHXKesJk2alG3+nISEhBAfH5/PqfNHx44d2bZtG7/99hsjR450Oo6IyEWpZqtmS9Gj5rWQJSQk0LdvX+rVq0dERARxcXGMHj0agD///JOIiAhef/11XnvtNQAmTpzI+vXriYiIoF69erz55ps5Lr9cuXK0aNGC+vXra/C/iMhlUs0WKXqMtTbXMzdt2tSuX7++AOMUXzVq1GD9+vWEhoY6HUXEKxljNlhrmzqdozCpZhcc1WyRgpVTzdaRVxERERHxGPrAVhGxc+dOpyOIiEguqWaLOEdHXkVERETEY6h5LUKCg4MB1x59YGAgkZGR7q8PP/wQcI2zatCgAREREdx0003s2rXLycgiIsWWaraIMzRsoIg6fceV81m+fDmhoaGMGjWKMWPG8PbbbxduOBERyUY1W6Tw6MirB2vevDl79+51OoaIiOSCarZI/lDzWkT99ttv2U5BrVq16px5Fi9eTLdu3Qo/nIiIZKOaLVJ4NGygiMrpFFTr1q05duwYwcHBPP/884UbTEREzqGaLVJ4dOTVAy1fvpxdu3YRGRnJqFGjnI4jIiI5UM0WyV9qXj2Un58fEyZM4MMPP+TYsWNOxxERkRyoZovkHzWvRdTZ46cmTpx4zjyVK1emV69eTJkyxYGEIiJymmq2SOEx1tpcz6z7ZIuIp8rpPtneSjVbRDxVTjVbR15FRERExGOoeRURERERj6HmtYhbvHgxdevWpXbt2owbN87pOCIikgPVbJGCp+a1CMvIyGDw4MEsWrSIuLg4ZsyYQVxcnNOxRETkPFSzRQqHmtcibN26ddSuXZtatWrh7+9Pz549mTt3rtOxRETkPFSzRQqHmtcibO/evVSrVs39uGrVqrovtohIEaWaLVI41LyKiIiIiMdQ81qEhYWFsXv3bvfjPXv2EBYW5mAiERG5ENVskcKh5rUIi4qKYvv27ezYsYPU1FRmzpxJ165dnY4lIiLnoZotUjj8nA4gF+bn58fkyZOJiYkhIyOD/v37Ex4e7nQsERE5D9VskcKh28OKSLGg28OKiHgO3R5WRERERLyCmlcRERER8RhqXkVERETEY6h5FRERERGPoeZVRERERDyGmlcRERER8RhqXkVERETEY+TpOq/GmMPAroKLIyJSYKpba8s7HaIwqWaLiAe7YM3OU/MqIiIiIuIkDRsQEREREY+h5lVEREREPIaaVxERERHxGGpeRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpeRURERMRjqHkVEREREY/x/2sTQmgtr1cVAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "AdvWWT.plot_sankey_report(margin=.6, gap=.8)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Adv WWT: Spent GAC"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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>87.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Solids flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>78.30</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>8.70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>VS flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>70.47</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>54.81</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>90000.00</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Parameter Unit    Amount\n",
       "0   Mass flow      kg   87.00   \n",
       "1   Solids flow    kg   78.30   \n",
       "2   Moisture flow  kg   8.70    \n",
       "3   VS flow        kg   70.47   \n",
       "4   Carbon flow    kg   54.81   \n",
       "5   PFOA           μg   90000.00\n",
       "6   PFOS           μg   90000.00\n",
       "7   PFBA           μg   90000.00\n",
       "8   PFPeA          μg   90000.00\n",
       "9   PFHxA          μg   90000.00\n",
       "10  PFHpA          μg   90000.00\n",
       "11  PFNA           μg   90000.00\n",
       "12  PFDA           μg   90000.00\n",
       "13  PFBS           μg   90000.00\n",
       "14  PFHxS          μg   90000.00"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "AdvWWT.SpentGAC.report()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Flow Type: SpentGAC \n",
      "GAC mass: 87.0 kg\n"
     ]
    }
   ],
   "source": [
    "print(\"Flow Type: {} \\nGAC mass: {} kg\".format(AdvWWT.SpentGAC.FlowType,AdvWWT.SpentGAC.GAC_mass))"
   ]
  }
 ],
 "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
}
