{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Landfill"
   ]
  },
  {
   "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",
    "While the model describes a generic landfill, there are in practice four types of landfills (MSW, hazardous waste, C&D waste, ash). Each of these types of landfill includes the same basic processes as illustrated in Figure below and each type may be represented by adjusting default parameters. They may be unlined (e.g., many C&D landfills), have a single liner system (e.g., most MSW landfills), or be double-lined (e.g., hazardous waste landfills). The landfills may or may not have active gas collection and control systems. If there is an active gas control system utilizing a flare or engine, then some fraction of the volatilized or aerosolized PFAS may be destroyed. The PFAS that is not released remains stored in the landfill.\n",
    "\n",
    "\n",
    "To calculate PFAS release, the disposed PFAS-containing waste material is assumed to be well mixed with the bulk MSW. The area occupied by the disposed waste is estimated based on a default parameter that represents the mass of total waste that can be disposed per unit area. The area occupied by the PFAS-containing waste mixed in with MSW is coupled with the leachate generation rate (volume/area-time) to estimate the volume of leachate produced. The partitioning of PFAS from the solid waste to the aqueous phase is modeled using a liquid-solid partition coefficient normalized to the amount of organic carbon. Model predictions are based on achievement of equilibrium. By default, it is assumed that no volatilization occurs, but a user can enter a fraction of PFAS that is volatilized.\n",
    "\n",
    "\n",
    "The partition coefficient is used to estimate the concentration of PFAS in the liquid and solids. The concentration in the liquid changes throughout the year as is leached to the leachate collection system (i.e., it is assumed that leaching is uniform throughout the year and continuously removes PFAS from the mixture). The leachate collection efficiency is used to calculate the fraction of leachate that is collected and subsequently treated. The landfill process models in the SAT is designed so that by changing default parameters, a variety of landfill processes can be represented. For example, a C&D landfill without a liner would be modeled with no collection efficiency, all the generated leachate would be released to groundwater. A double-lined hazardous waste landfill could be modeled by increasing the leachate collection efficiency. \n",
    "\n",
    "\n",
    "<img src=\"../Images/ProcessModels/Landfill_Diagram.jpg\" alt=\"Drawing\" style=\"width: 700px;\"/>\n",
    "\n",
    "\n",
    "### Assumptions and Limitations:\n",
    "\n",
    "1.\tThe organic carbon-normalized partition coefficient assumes that the organic carbon in different materials generally have the PFAS sorption capacity. However, material-specific coefficients were developed for activated carbon, \n",
    "2.\tVolatilization is assumed to be zero by default due to a lack of data. However, the user may assign a fraction of the PFAS to be volatilized.\n",
    "3.\tFuture work is also required to implement a dynamic (i.e., non-equilibrium) model to account for changes in the organic C content of over time as landfilled materials decompose, and the effectiveness of cover system improves. \n",
    "4.\tPFAS release is based on 1 year. A longer time horizon would result in additional PFAS release."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for Landfill model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Category</th>\n",
       "      <th>Dictonary_Name</th>\n",
       "      <th>Parameter Name</th>\n",
       "      <th>Parameter Description</th>\n",
       "      <th>amount</th>\n",
       "      <th>unit</th>\n",
       "      <th>Reference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>LF MSW Properties</td>\n",
       "      <td>LFMSW</td>\n",
       "      <td>frac_of_msw</td>\n",
       "      <td>Fraction of PFAS containing waste mixed with MSW</td>\n",
       "      <td>0.010</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>LF MSW Properties</td>\n",
       "      <td>LFMSW</td>\n",
       "      <td>bulk_dens</td>\n",
       "      <td>MSW bulk density - wet</td>\n",
       "      <td>890.000</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>LF MSW Properties</td>\n",
       "      <td>LFMSW</td>\n",
       "      <td>ts_cont</td>\n",
       "      <td>Total solids content of the bulk MSW in the landfill  - wet</td>\n",
       "      <td>0.746</td>\n",
       "      <td>fraction wet weight</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>LF MSW Properties</td>\n",
       "      <td>LFMSW</td>\n",
       "      <td>C_cont</td>\n",
       "      <td>MSW organic C content -dry</td>\n",
       "      <td>0.331</td>\n",
       "      <td>fraction TS</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>LF Water Balance</td>\n",
       "      <td>Water_Blnc</td>\n",
       "      <td>leach_gpad</td>\n",
       "      <td>Leachate generation rate per area</td>\n",
       "      <td>500.000</td>\n",
       "      <td>gal/acre-day</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>LF Water Balance</td>\n",
       "      <td>Water_Blnc</td>\n",
       "      <td>is_leach_col</td>\n",
       "      <td>Is leachate collected</td>\n",
       "      <td>1.000</td>\n",
       "      <td>1:TRUE/0:FALSE</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>LF Water Balance</td>\n",
       "      <td>Water_Blnc</td>\n",
       "      <td>frac_leach_col</td>\n",
       "      <td>Fraction of leachate collected</td>\n",
       "      <td>0.998</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Volatilization</td>\n",
       "      <td>Volatilization</td>\n",
       "      <td>frac_vol_loss</td>\n",
       "      <td>Fraction of PFAS lost to volatilization</td>\n",
       "      <td>0.000</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>LF Parameters</td>\n",
       "      <td>LF</td>\n",
       "      <td>lf_ton_area</td>\n",
       "      <td>Landfill capacity per acre</td>\n",
       "      <td>100000.000</td>\n",
       "      <td>Mg/acre</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log Koc (Landfill)</td>\n",
       "      <td>2.190</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log Koc (Landfill)</td>\n",
       "      <td>3.040</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log Koc (Landfill)</td>\n",
       "      <td>1.880</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log Koc (Landfill)</td>\n",
       "      <td>1.370</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log Koc (Landfill)</td>\n",
       "      <td>1.770</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log Koc (Landfill)</td>\n",
       "      <td>1.970</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log Koc (Landfill)</td>\n",
       "      <td>2.630</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log Koc (Landfill)</td>\n",
       "      <td>3.240</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log Koc (Landfill)</td>\n",
       "      <td>1.510</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log Koc (Landfill)</td>\n",
       "      <td>2.790</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log K (GAC)</td>\n",
       "      <td>4.870</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log K (GAC)</td>\n",
       "      <td>5.080</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log K (GAC)</td>\n",
       "      <td>4.490</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log K (GAC)</td>\n",
       "      <td>4.590</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log K (GAC)</td>\n",
       "      <td>4.680</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log K (GAC)</td>\n",
       "      <td>4.780</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log K (GAC)</td>\n",
       "      <td>4.970</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log K (GAC)</td>\n",
       "      <td>5.060</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log K (GAC)</td>\n",
       "      <td>4.680</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>Log partition coefficient for GAC</td>\n",
       "      <td>LogPartCoef_GAC</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log K (GAC)</td>\n",
       "      <td>4.880</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log K (SpentIER)</td>\n",
       "      <td>5.660</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log K (SpentIER)</td>\n",
       "      <td>7.050</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log K (SpentIER)</td>\n",
       "      <td>4.380</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log K (SpentIER)</td>\n",
       "      <td>4.780</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log K (SpentIER)</td>\n",
       "      <td>5.150</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log K (SpentIER)</td>\n",
       "      <td>5.440</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log K (SpentIER)</td>\n",
       "      <td>5.870</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log K (SpentIER)</td>\n",
       "      <td>6.250</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log K (SpentIER)</td>\n",
       "      <td>5.870</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>38</th>\n",
       "      <td>Log partition coefficient for Spent IER</td>\n",
       "      <td>LogPartCoef_SpentIER</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log K (SpentIER)</td>\n",
       "      <td>6.460</td>\n",
       "      <td>log L/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                   Category        Dictonary_Name  \\\n",
       "0   LF MSW Properties                        LFMSW                  \n",
       "1   LF MSW Properties                        LFMSW                  \n",
       "2   LF MSW Properties                        LFMSW                  \n",
       "3   LF MSW Properties                        LFMSW                  \n",
       "4   LF Water Balance                         Water_Blnc             \n",
       "5   LF Water Balance                         Water_Blnc             \n",
       "6   LF Water Balance                         Water_Blnc             \n",
       "7   Volatilization                           Volatilization         \n",
       "8   LF Parameters                            LF                     \n",
       "9   Log partition coefficient                LogPartCoef            \n",
       "10  Log partition coefficient                LogPartCoef            \n",
       "11  Log partition coefficient                LogPartCoef            \n",
       "12  Log partition coefficient                LogPartCoef            \n",
       "13  Log partition coefficient                LogPartCoef            \n",
       "14  Log partition coefficient                LogPartCoef            \n",
       "15  Log partition coefficient                LogPartCoef            \n",
       "16  Log partition coefficient                LogPartCoef            \n",
       "17  Log partition coefficient                LogPartCoef            \n",
       "18  Log partition coefficient                LogPartCoef            \n",
       "19  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "20  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "21  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "22  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "23  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "24  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "25  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "26  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "27  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "28  Log partition coefficient for GAC        LogPartCoef_GAC        \n",
       "29  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "30  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "31  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "32  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "33  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "34  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "35  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "36  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "37  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "38  Log partition coefficient for Spent IER  LogPartCoef_SpentIER   \n",
       "\n",
       "    Parameter Name  \\\n",
       "0   frac_of_msw      \n",
       "1   bulk_dens        \n",
       "2   ts_cont          \n",
       "3   C_cont           \n",
       "4   leach_gpad       \n",
       "5   is_leach_col     \n",
       "6   frac_leach_col   \n",
       "7   frac_vol_loss    \n",
       "8   lf_ton_area      \n",
       "9   PFOA             \n",
       "10  PFOS             \n",
       "11  PFBA             \n",
       "12  PFPeA            \n",
       "13  PFHxA            \n",
       "14  PFHpA            \n",
       "15  PFNA             \n",
       "16  PFDA             \n",
       "17  PFBS             \n",
       "18  PFHxS            \n",
       "19  PFOA             \n",
       "20  PFOS             \n",
       "21  PFBA             \n",
       "22  PFPeA            \n",
       "23  PFHxA            \n",
       "24  PFHpA            \n",
       "25  PFNA             \n",
       "26  PFDA             \n",
       "27  PFBS             \n",
       "28  PFHxS            \n",
       "29  PFOA             \n",
       "30  PFOS             \n",
       "31  PFBA             \n",
       "32  PFPeA            \n",
       "33  PFHxA            \n",
       "34  PFHpA            \n",
       "35  PFNA             \n",
       "36  PFDA             \n",
       "37  PFBS             \n",
       "38  PFHxS            \n",
       "\n",
       "                                          Parameter Description      amount  \\\n",
       "0   Fraction of PFAS containing waste mixed with MSW             0.010        \n",
       "1   MSW bulk density - wet                                       890.000      \n",
       "2   Total solids content of the bulk MSW in the landfill  - wet  0.746        \n",
       "3   MSW organic C content -dry                                   0.331        \n",
       "4   Leachate generation rate per area                            500.000      \n",
       "5   Is leachate collected                                        1.000        \n",
       "6   Fraction of leachate collected                               0.998        \n",
       "7   Fraction of PFAS lost to volatilization                      0.000        \n",
       "8   Landfill capacity per acre                                   100000.000   \n",
       "9   PFOA Log Koc (Landfill)                                      2.190        \n",
       "10  PFOS Log Koc (Landfill)                                      3.040        \n",
       "11  PFBA Log Koc (Landfill)                                      1.880        \n",
       "12  PFPeA Log Koc (Landfill)                                     1.370        \n",
       "13  PFHxA Log Koc (Landfill)                                     1.770        \n",
       "14  PFHpA Log Koc (Landfill)                                     1.970        \n",
       "15  PFNA Log Koc (Landfill)                                      2.630        \n",
       "16  PFDA Log Koc (Landfill)                                      3.240        \n",
       "17  PFBS Log Koc (Landfill)                                      1.510        \n",
       "18  PFHxS Log Koc (Landfill)                                     2.790        \n",
       "19  PFOA Log K (GAC)                                             4.870        \n",
       "20  PFOS Log K (GAC)                                             5.080        \n",
       "21  PFBA Log K (GAC)                                             4.490        \n",
       "22  PFPeA Log K (GAC)                                            4.590        \n",
       "23  PFHxA Log K (GAC)                                            4.680        \n",
       "24  PFHpA Log K (GAC)                                            4.780        \n",
       "25  PFNA Log K (GAC)                                             4.970        \n",
       "26  PFDA Log K (GAC)                                             5.060        \n",
       "27  PFBS Log K (GAC)                                             4.680        \n",
       "28  PFHxS Log K (GAC)                                            4.880        \n",
       "29  PFOA Log K (SpentIER)                                        5.660        \n",
       "30  PFOS Log K (SpentIER)                                        7.050        \n",
       "31  PFBA Log K (SpentIER)                                        4.380        \n",
       "32  PFPeA Log K (SpentIER)                                       4.780        \n",
       "33  PFHxA Log K (SpentIER)                                       5.150        \n",
       "34  PFHpA Log K (SpentIER)                                       5.440        \n",
       "35  PFNA Log K (SpentIER)                                        5.870        \n",
       "36  PFDA Log K (SpentIER)                                        6.250        \n",
       "37  PFBS Log K (SpentIER)                                        5.870        \n",
       "38  PFHxS Log K (SpentIER)                                       6.460        \n",
       "\n",
       "                   unit  Reference  \n",
       "0   fraction                        \n",
       "1   kg/m3                           \n",
       "2   fraction wet weight             \n",
       "3   fraction TS                     \n",
       "4   gal/acre-day                    \n",
       "5   1:TRUE/0:FALSE                  \n",
       "6   fraction                        \n",
       "7   fraction                        \n",
       "8   Mg/acre                         \n",
       "9   log L/kg OC          [2,3,4,5]  \n",
       "10  log L/kg OC          [2,3,4,5]  \n",
       "11  log L/kg OC          [2,3,4,5]  \n",
       "12  log L/kg OC          [2,3,4,5]  \n",
       "13  log L/kg OC          [2,3,4,5]  \n",
       "14  log L/kg OC          [2,3,4,5]  \n",
       "15  log L/kg OC          [2,3,4,5]  \n",
       "16  log L/kg OC          [2,3,4,5]  \n",
       "17  log L/kg OC          [2,3,4,5]  \n",
       "18  log L/kg OC          [2,3,4,5]  \n",
       "19  log L/kg                        \n",
       "20  log L/kg                        \n",
       "21  log L/kg                        \n",
       "22  log L/kg                        \n",
       "23  log L/kg                        \n",
       "24  log L/kg                        \n",
       "25  log L/kg                        \n",
       "26  log L/kg                        \n",
       "27  log L/kg                        \n",
       "28  log L/kg                        \n",
       "29  log L/kg                        \n",
       "30  log L/kg                        \n",
       "31  log L/kg                        \n",
       "32  log L/kg                        \n",
       "33  log L/kg                        \n",
       "34  log L/kg                        \n",
       "35  log L/kg                        \n",
       "36  log L/kg                        \n",
       "37  log L/kg                        \n",
       "38  log L/kg                        "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Landfill = pspd.Landfill()\n",
    "Landfill.InputData.Data[['Category','Dictonary_Name','Parameter Name', 'Parameter Description', 'amount', 'unit','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming Waste to landfill\n",
    "\n",
    "Here are the properties of the incoming waste to landfill."
   ]
  },
  {
   "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>254.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>746.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>VS flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>134.62</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>84.074</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1070.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>840.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.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   254.0 \n",
       "2   Moisture flow  kg   746.0 \n",
       "3   VS flow        kg   134.62\n",
       "4   Carbon flow    kg   84.074\n",
       "5   PFOA           μg   1070.0\n",
       "6   PFOS           μg   840.0 \n",
       "7   PFBA           μg   0.0   \n",
       "8   PFPeA          μg   0.0   \n",
       "9   PFHxA          μg   0.0   \n",
       "10  PFHpA          μg   0.0   \n",
       "11  PFNA           μg   0.0   \n",
       "12  PFDA           μg   0.0   \n",
       "13  PFBS           μg   0.0   \n",
       "14  PFHxS          μg   0.0   "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncomingWaste = pspd.IncomFlow()\n",
    "IncomingWaste.set_flow('MSW', 1000)\n",
    "IncomingWaste.calc()\n",
    "MSW = IncomingWaste.Inc_flow\n",
    "MSW.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in landfill\n",
    "* A substantial fraction of PFAS is stored in landfill for longterm."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Volatilized</th>\n",
       "      <th>Collected Leachate</th>\n",
       "      <th>Fugitive Leachate</th>\n",
       "      <th>Remaining in Landfill</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>PFOA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0180</td>\n",
       "      <td>0.0</td>\n",
       "      <td>99.9819</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0026</td>\n",
       "      <td>0.0</td>\n",
       "      <td>99.9974</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Volatilized  Collected Leachate  Fugitive Leachate  \\\n",
       "PFOA   0.0          0.0180              0.0                 \n",
       "PFOS   0.0          0.0026              0.0                 \n",
       "PFBA  NaN          NaN                 NaN                  \n",
       "PFPeA NaN          NaN                 NaN                  \n",
       "PFHxA NaN          NaN                 NaN                  \n",
       "PFHpA NaN          NaN                 NaN                  \n",
       "PFNA  NaN          NaN                 NaN                  \n",
       "PFDA  NaN          NaN                 NaN                  \n",
       "PFBS  NaN          NaN                 NaN                  \n",
       "PFHxS NaN          NaN                 NaN                  \n",
       "\n",
       "       Remaining in Landfill  \n",
       "PFOA   99.9819                \n",
       "PFOS   99.9974                \n",
       "PFBA  NaN                     \n",
       "PFPeA NaN                     \n",
       "PFHxA NaN                     \n",
       "PFHpA NaN                     \n",
       "PFNA  NaN                     \n",
       "PFDA  NaN                     \n",
       "PFBS  NaN                     \n",
       "PFHxS NaN                     "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Landfill.calc(Inc_flow=MSW)\n",
    "Landfill.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAArVElEQVR4nO3deXxM9/7H8ddkIpYkglorKog1yUwWidAKSYUQ0rpFkuq9lmovWtelVXopuqmWi1pa6mpxiyhqa1ENCWopQUTFVhUkRcWeRLaZ7++P/MxtCIJMTpbP8/GYx8Oc+Z7v+ZzJ5O3kzPl+j04phRBCiOJno3UBQghRXkkACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoxPZhGtesWVO5uLhYqRQhhCib9u/fn6qUqnXn8ocKYBcXF+Li4oquKiGEKAd0Ot2ZgpbLKQghhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsHggvV6Pp6en5ZGUlPRI/bRr1w6ApKQkli5dalkeFxfHP/7xj6IoVYhSRfcwtyRq3bq1kqHI5Y+DgwNpaWlF1l9sbCxTp07lu+++K7I+hSjJdDrdfqVU6zuXyxGweCQuLi6kpqYCeUewHTt2BODSpUsEBwfj5ubGoEGDaNiwoaWdg4MDAGPGjGHHjh14enoyffp0YmNj6d69O2azGRcXF65du2bZTtOmTbl48SKXLl3ihRdewNfXF19fX3bu3Fms+yuENUgAiwe6deuW5fRDz54979v23XffJSgoiCNHjtCrVy/Onj17V5vJkyfTvn174uPjGTFihGW5jY0Nzz33HKtXrwbg559/pmHDhtSpU4fhw4czYsQI9u3bx6pVqxg0aFDR7qQQGnio2dBE+VS5cmXi4+ML1fann36yBGhISAjVq1d/qG2Fh4fz3nvvMWDAAKKioggPDwcgOjqaxMRES7sbN26QlpZmOaoWojSSABaPxNbWFrPZDEBmZmaR9du2bVt+/fVXLl26xJo1axg3bhwAZrOZPXv2UKlSpSLblhBak1MQ4pG4uLiwf/9+AFatWmVZ/vTTT/PNN98AsHnzZq5evXrXuo6Ojty8ebPAfnU6HT179mTkyJG0bNmSJ554AoDOnTsza9YsS7vCHpELUZJJAItHMmHCBIYPH07r1q3R6/X5lm/evBl3d3dWrFhB3bp1cXR0zLeuwWBAr9djNBqZPn36XX2Hh4fz9ddfW04/AMycOZO4uDgMBgOtWrVi7ty51ts5IYqJXIYmilRWVhZ6vR5bW1t2797NkCFD5GhVlHv3ugxNzgGLInX27Fn69OmD2WzGzs6O+fPna12SECWWBLAoUk2bNuXgwYNalyFEqSDngEWBrHl518P2HRsby65du6xUjRDakQAWJZ4EsCirJIBFoZ06dYqQkBB8fHxo3749x44dA2D9+vW0adMGLy8vOnXqxMWLFwFIS0tjwIABeHh4YDAY8l2uNnbsWIxGI/7+/pb2BfWTlJTE3LlzmT59Op6enuzYsUOGJYuyQylV6IePj48S5YO9vf1dy4KCgtSJEyeUUkrt2bNHBQYGKqWUunLlijKbzUoppebPn69GjhyplFLqrbfeUsOHD7esf+XKFaWUUoBat26dUkqpUaNGqffff/++/UyYMEFNmTLF0k9kZKTasWOHUkqpM2fOqBYtWhTNTgthJUCcKiBT5Us4UShpaWns2rWL3r17W5ZlZWUBkJycTHh4OOfPnyc7O5tGjRoBecOHo6KiLO1vD0u2s7Oje/fuAPj4+PDjjz/et587ybBkUVZIAItCMZvNVKtWrcBreocNG8bIkSMJCwsjNjaWiRMn3revChUqoNPpgLy5hnNzcx+qHxmWLMoKOQcsCqVq1ao0atSIFStWAHmnrg4dOgTA9evXqV+/PgCLFi2yrBMcHMycOXMszwsalvxn9+rnzqHLMixZlBUSwKJAGRkZODs7Wx7Tpk1jyZIlLFiwAKPRiJubG2vXrgVg4sSJ9O7dGx8fH2rWrGnpY9y4cVy9ehV3d3eMRiMxMTH33ea9+unRowerV6+2fAknw5JFWSFDkYUQwsrkjhhCCFHCSAALIYRGJICFEEIjEsCiWGzatInmzZvj6urK5MmTtS5HiBJBAlhYnclk4rXXXmPjxo0kJiaybNmyfAMphCivJICF1e3duxdXV1caN26MnZ0dERERlkvYhCjPJICF1aWkpNCgQQPLc2dnZ1JSUjSsSIiSQQJYCCE0IgEsrK5+/fqcO3fO8jw5Odky5FiI8kwCWFidr68vJ0+e5PTp02RnZxMVFUVYWJjWZQmhOZkNTVidra0ts2fPpmPHjuj1el5++WXc3Ny0LksIzckRsLC63Nxc1n73Pb9f/IPaTzZgyJAhWpckRIkgASysbuGiRXzx+WfkZmXy887t/PONUVqXJESJIAEsrC4yIoLDhw8DsGDBAv79iYyEEwIkgEUxsLe3x93dHYAWLVpQq1YtjSsSomSQABZCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwKBZKKQDq1q2rcSVClBwSwKJYTJkyBQAXFxdtCxGiBJEAFsVi9OjRNGrUCBsb+cgJcZv8NgirCwoKAiAmJkbjSrSTlpbG8ePHtS5DlDASwMJqtmzZgk6nIyYmhpUrV9KwYUOtS9LMp59+ipdPa44ePap1KaIE0d3+cqQwWrdureLi4qxYTvlhMpnoFdGXa9evP3ZfsT9uKoKK7q9jcEih295ZT+069WhlMBZ1SaVK0m+n+P1mLk9UNJNwII6aNWtqXZIoRjqdbr9SqvWdy221KEZATk4Ox48f50y6DfaeoY/VV+Uml8k6f6KIKrtbFdc2HHH0LXT7qn5VyDz3C45e3bCp5IgCjlitulLC6EudpzzI+PkbuvZ4np3btmJnZ6d1VUJjcgSsoZSUFIw+vti0G0CVZu20LkcUA6XMpH03mS7erixZvBCdTqd1SaIY3OsIWM4Ba6h+/fr88P16MrZ+TvbFU1qXI4qBTmeDfZcRbIjZxeRPpmhdjtCYBLDGfHx8+PKLedxcP4nctCtalyOKgY1dZRzCxvLhx1NZu3at1uUIDUkAlwB9+vTmjX+8Rvr6jzDnZGldjigGtlVr4dh9DC/1H8ihQ4e0LkdoRAK4hJjwzjie9TeS/uNMlDJrXY4oBhWfbE6lgEEEdw3lwoULWpcjNCABXELodDq+XvQVDexukb47SutyRDGxbxmA2bUDnbv14NatW1qXI4qZBHAJUqlSJX74fh36U9tJT9ymdTmimNi3jSQl14G+f+vPw1yVJEo/CeASpk6dOkRv2kDmjgVkpRzTuhxRDHQ6HfbBw4jZd5gJ7753z3Z6vR5PT0/c3d3p0aMH165ds2pdgwYNIjEx8b5t5s6dy+LFi61aR1km1wGXUOvXryey38tUi/gY26q1tS5HFIPctCtcj3qL+bNnEBERftfrDg4OpKWlAdCvXz+aNWvG2LFji7tM8QjkOuBSpkePHox/+y3S1k/CnC3nBssDW4caOPb4F4MGD2Hfvn33bdu2bVtSUlIAOHXqFCEhIfj4+NC+fXuOHcv7y6l///4MGTIEf39/GjduTGxsLAMHDqRly5b079/f0teQIUNo3bo1bm5uTJgwwbK8Y8eO3D7gcnBwYOzYsRiNRvz9/bl48SIAEydOZOrUqZb2o0ePxs/Pj2bNmrFjxw4AMjIy6NOnD61ataJnz560adMGOZDLIwFcgo168w26Bz1D+qbpKLNJ63JEMbCr05gqQUMJ6R5GcnJygW1MJhNbtmwhLCwMgFdffZVZs2axf/9+pk6dytChQy1tr169yu7du5k+fTphYWGMGDGCI0eOcPjwYeLj4wH48MMPiYuLIyEhgW3btpGQkHDXNtPT0/H39+fQoUMEBAQwf/78AmvLzc1l7969zJgxg3fffReAzz77jOrVq5OYmMj777/P/v37H+ctKlMkgEswnU7Hl/Pn0aSaDek//VfrckQxqdKsLTq3rnQKCSU9Pd2y/NatW3h6elK3bl0uXrxIcHAwaWlp7Nq1i969e+Pp6cnf//53zp8/b1mnR48e6HQ6PDw8qFOnDh4eHtjY2ODm5kZSUhIA33zzDd7e3nh5eXHkyJECz/va2dnRvXt3IG/w0O117/SXv/zlrjY//fQTERERALi7u2MwGB73LSozJIBLODs7OzauX0OllP2kH/5R63JEMbH3/QupdnXoFfEiZnPedeGVK1cmPj6eM2fOoJRizpw5mM1mqlWrRnx8vOXx5ykvK1asCICNjY3l37ef5+bmcvr0aaZOncqWLVtISEggNDSUzMzMu+qpUKGCZd4KvV5Pbm5ugXXf3sb92oj/kQAuBZ544gmif9hA1q6vyTx7WOtyRDHQ6XTYBw1lT2ISb435V77XqlSpwsyZM/n3v/9NlSpVaNSoEStWrADy7r33MCPrbty4gb29PU5OTly8eJGNGzcW6X4APP3003zzzTcAJCYmcviwfIZvkwAuJVq0aMHK5UtJ2ziVnKvnH7yCKPV0thVwCB3NFwu/ZuHCRfle8/LywmAwsGzZMpYsWcKCBQswGo24ubk91PwSRqMRLy8vWrRowYsvvsjTTz9d1LvB0KFDuXTpEq1atWLcuHG4ubnh5ORU5NspjeQytFJm1uw5jJs0Fafwj7GpaK91OaIYZKee5cbKcWz6bi3PPPOM1uU8NJPJRE5ODpUqVeLUqVN06tSJ48ePl6v5kGVC9jJi2OuvcfjIEVZumILjc++gs9FrXZKwMruaT2HfeTg9nv8LB/b9TKNGjbQu6aFkZGQQGBhITk4OSik+++yzchW+9yNHwKVQbm4uQZ27kphWBYegV7UuRxSTtP3rqXpmG/FxP1O1alWtyxEPQQZilCG2trasX70SxyvHSDv4vdbliGJi792dG05NeO4vvTGZ5LrwskBOQdzDn4d9Fqfx48cTEBBAp06d7tvOycmJ6B824OPnz61qT1K5kdddbZI/H4iNXWXQ6dDbV6dm6Ej0DtX/t9wm7//fGsFDqeTckowTu7m0+kOeHPQ5FZ5oAOTdQufqlvlknjkE6NDZ2lHzudFUqFa3yPdd3J9Op8O+4yscWvsew4aP5LPZnxZJv+3atWPXrl1F0pd4OBLAJcx77917MpY7NWnShHWrV9E1rCe2vT+0hOaf1YmchL6KE1e3LeL6nm+o0env+Zb/WfrR7VR0bkV64naqte8LQMbRHZhuXqbewNnodDbk3khFZ1fpMfZQPA6d3hb7bm/xddRbeLi3ZMjgwY/dp4SvduQUxAPExsbSsWNHevXqRYsWLejbt69lysB9+/bRrl07jEYjfn5+3Lx5k8zMTAYMGICHhwdeXl7ExMQAsHDhQp5//nmCg4NxcXFh9uzZTJs2DS8vL/z9/blyJe92RP3792flypUAuLi4MGHCBLy9vfHw8LCM8b906RLBwcG4ubmxePFiKunh+pr3MN26cc/9qNTAnZyrv9/zdXP2LbKSj/BE1+GkH9tuWW5Ku4LeoQY6Xd5HxbZqTfSVHB7jHRWPS1/JAcewsYwaM5atW7c+dn8ODnk/z/t91oV1SAAXwsGDB5kxYwaJiYn89ttv7Ny5k+zsbMLDw/n00085dOgQ0dHRVK5cmTlz5qDT6Th8+DDLli2jX79+lpFFv/zyC99++y379u1j7NixVKlShYMHD9K2bdt7TulXs2ZNDhw4wJAhQyyTnrz77rsEBQVx5MgRevXqxZUrV3ix1/Okff8xypRTYD+3Tu3FrpaL5fnFZf/i96+GcX7xSAAyTu6hUiMfKtSoj76SI1kXfgWgSov2ZPy6l9+/GsaVrf+Rm4eWEBVq1Meh65v07NWHEydOFFm/BX3WhfVIABeCn58fzs7O2NjY4OnpSVJSEsePH6devXr4+voCULVqVWxtbfnpp5946aWXgLzBEw0bNrT8ggQGBuLo6EitWrVwcnKiR48eAHh4eDzy2PqQkBCqV6/O+xMn4NvUmfSt8/IdtdwOWnNWBk7+vS3L60RO4skBs6j3t2kAZBzdjn3LAACqtAwg4/8nhLetWpP6r8yjekA/dDobLkaN5VZS/OO+paIIVGpooEKbSJ7t0tXyF9TjKuizLqxHzgEXwp/H0D/OGPc7x+L/eZz+446t1+v1fLsiCi9ff67uX4ND655Awed672S6dZPMMwlkX0oCdKDMgI5qgQPR6XTobCtQuUlrKjdpjU2Vatw6uYfKLp6F3m9hPQ7GENKuphD63F/YvvVHKlSo8Fj9FdVnXRSOHAE/oubNm3P+/HnLvK03b94kNzeX9u3bs2TJEgBOnDjB2bNnad68eZFu+89j6zdv3szVq1eBvHN5W37YgPnQejJ+3Vvo/jKO78TeLRDnIV/hPORLnIcuxLZaHbKSj5B14Vdyb14G8q6IyLl0Gr1MEF+i2Af050RqFq8OHirnbEsZCeBHZGdnx/Llyxk2bBhGo5Hg4GAyMzMZOnQoZrMZDw8PwsPDWbhwYb6jiqIwYcIENm/ejLu7OytWrKBu3bo4OjoC8NRTT7Fh/VoyomeDuXBHL+lHt1GlWdt8y6o0a0d64jbMGde5tOo9fl8wlPNfvg42eqr6dC/S/RGPR2ejxz5kJN9u2sq06TO0Lkc8BBkJVwplZWWh1+uxtbVl9+7dDBkyxDK59m1Llixl8D/fpFrEJ+jtq2tTqChWudcvci1qNN8sWUS3bt20Lkf8icwFUYacPXuWPn36YDabsbOzK/DuBH37vsgviUeZFzWFqr0naVClKG62TnVwCB1NRN+/knBwPy4uLlqXJB5AArgUatq0KQcPHnxgO7PZDP8/ibYoJ1TeEGWd/NxLBTkHXIINHDiQ2rVr4+7u/tDr/ve/X/PZfxZi33WUFSoTJVHOtQtcXDqGFVFLadiwYZH1++cbbx47dgxPT0+8vLw4deoU7dq1AyApKcnyOY2NjbXcvkjcnwRwCda/f382bdr00Ovt3r2bIcOG4xj2L/T21Yq+MFHimLPSSVv3AXZ2dnTp0sVq21mzZg29evXi4MGDNGnSRIYxPyYJ4BIsICCAGjVqPNQ6Z86cITTseeyDh+Ub+SbKLmU2kbZhKn3CQgq8Dnj9+vW0adMGLy8vOnXqlO+W8gMHDqRjx440btyYmTNnWtb58MMPadasGc888wzHjx8HYMOGDcyYMYPPP/+cwMBA4H/DmMWjkQAuQ27evMmzXbqh93yOyk18tS5HFJP0bV/i9mRVPp89q8DXn3nmGfbs2cPBgweJiIjgk08+sbx27NgxfvjhB/bu3cu7775LTk4O+/fvJyoqivj4eDZs2GC51r1bt24MHjyYESNGWOY4EY9HvoQrI0wmEz17hXPNwQV77zCtyxHFJD1+Iw6piazfvxdb24J/nZOTkwkPD+f8+fNkZ2fnu6NGaGgoFStWpGLFitSuXZuLFy+yY8cOevbsSZUqVQAIC5PPk7XIEXAZ8c83RnHgt4vYB74q34CXE7eS4snZt5zoHzZQrVq1e7YbNmwYr7/+OocPH2bevHn5bjsvQ4+1JQFcBnwxfz6LolZi320UOr38UVMe5FxOJn3TNNZ+uxJXV9f7tr1+/Tr169cHYNGiRfdtC3nfPaxZs4Zbt25x8+ZN1q9fXyQ1i7tJAJdgkZGRtG3bluPHj+Ps7MyCBQvuahMbG8vIUWNwDBuHvrKjBlWK4ma6dYOb6z5g+tSP6dChQ77XMjIycHZ2tjymTZvGxIkT6d27Nz4+PtSsWfOB/Xt7exMeHo7RaKRr166WGf9E0ZOhyKXYr7/+Sus2bakY/E+ZnaycUKYcbq5+l5d6BDFz+r+1LkcUktyUs4y5du0az3buSgW/CAnfckIpRXrMF3g3qcf0qZ88eAVR4skJw1IoNzeX7s+/QHptdxyMIVqXI4pJxoG1VE8/y+qY3ej1eq3LEUVAArgUGjz0dRIvpOEYNlzrUkQxyTi1D1P8OqLj9lqmHhWln5yCKGVmfDqLFd9txqHrG+hs5CioPMi+lETGj7PYsH5toeZ42LRpE82bN8fV1ZXJkyff9fqZM2d49tlnMRgMdOzYkeTkZMtro0ePxt3dHXd3d5YvX25ZvnXrVry9vXF3d6dfv36Wy9WuXr1Kz549MRgM+Pn58csvvxTBHpcfEsClyKZNmxg38T0cwsZiU9Fe63JEMTClX+PmuknMnT0Tf3//B7c3mXjttdfYuHEjiYmJLFu2jMTExHxt3nzzTf72t7+RkJDA+PHjefvttwH4/vvvOXDgAPHx8fz8889MnTqVGzduYDab6devH1FRUfzyyy80bNjQcjnbpEmT8PT0JCEhgcWLFzN8uPxV9jAkgEuJo0eP0ieyLw6ho6hQra7W5YhioHKzSft+Mq+9MoCXXupbqHX27t2Lq6srjRs3xs7OjoiICNauXZuvTWJiIkFBQUDejWJvv56YmEhAQAC2trbY29tjMBjYtGkTly9fxs7OjmbNmgEQHBzMqlWr7uqrRYsWJCUlWeaaEA8mAVwKpKam0qlLNyo+3Y9Kzm5alyOKgVKKtOg5PG1oyqQP3iv0eikpKTRo0MDy3NnZmZSUlHxtjEYj3377LQCrV6/m5s2bXL58GaPRyKZNm8jIyCA1NZWYmBjOnTtHzZo1yc3N5fYlqCtXruTcuXN39bV3717OnDmT75SGuD8J4BIuOzubrj2eJ6uBH/buz2pdjigm6XtXUldd4ZulX2NjU7S/plOnTmXbtm14eXmxbds26tevj16vp3PnznTr1o127dpZBgHp9Xp0Oh1RUVGMGDECPz8/HB0dLVdhjBkzhmvXruHp6cmsWbPw8vKSKzQeglwFUYIppej/8iv8dgMcur+kdTmimKQf3wlHN/PjgTjLhDiFVb9+fcvRKeRNxHN7GPJtTz75pOWoNS0tjVWrVlnmkhg7dixjx44F4MUXX7Scdmjbti07duwA8u7EfeLECQCqVq3KV199BeR9Xhs1akTjxo0fco/LLzkCLsE++vgTNsTuxr7LP9Hp5EdVHmRd+JVbMXPZvOE7nnzyyYde39fXl5MnT3L69Gmys7OJioq6azaz1NTUvNtVAR999BEDBw4E8r7Au3z5MgAJCQkkJCTQuXNnAP7444+8+rKy+Pjjjxk8eDCQNyAoOzsbgP/85z8EBARQtWrVR9jz8kmOgEuotWvXMmnKNKqFf4KNXWWtyxHFIPfmZdLWf8SiBfPx9vZ+pD5sbW2ZPXs2Xbp0wWQyMXDgQNzc3Bg/fjytW7cmLCyM2NhY3n77bXQ6HQEBAcyZMweAnJwc2rdvD+Qd2X799deWKS6nTJnCd999h9lsZsiQIZYv3o4ePUq/fv3Q6XS4ubkVOF+JuDeZC6IEOnToEM90DMIhbBwV6zXTuhxRDMw5mdxYMZaRr7zE+HfGaV2OKGIyF0QpceHCBYK7hlIp4BUJ33JCKTPpm2fS+Wkf3hk3VutyRDGSAC5Bbt26ReduPTA3DcS+ZXutyxHFJH3XMp6qlMV/F34pk+mXM3IOuIRQStH3b/1JMTng4B+udTmimKQnxmB7+id+OBCX7+4UonyQI+ASYvzEd4nZ9wv2nYbJUVA5kZVylMwdXxG9aQO1a9fWuhyhAQngEiAqajkz5szDofvb2FSQo6DyIPf6H9z8/hOW/XcR7u7uWpcjNCIBrLF9+/YxaPBQHMP+hd6hutbliGJgzsogbf2HTPjXaLp37651OUJDEsAaSk5OJqR7GFWefQ272jJ6qDxQZhNpP0zjueAOvPnGSK3LERqTL+E0kpmZSaeQUDLtn8Remcg4seuR+7p5YANZ548XYXX5VWnqT5VmbQvdPvPsYTLPHsbRuzv6KtYdFZVzOZnru5c/uGEBVE6W5d862/+d+rGp5ECNTn9/7NruVMnFSMaeb2hew47/fDFXzvULCWCt2NnZ4eHhgfO163Ar4bH62nImvmiKuof0IzH4Nyj8/MNb9n8HwJUfZgNQu249PDx9rFLblu3fFUk/Kvd/YWxKy8KY9Xg/kzv99utJkn9aQm17Pd8f2EeFChWKtH9RSimlCv3w8fFRQhRWbGysAhSgli9frnU59/T7778rg8GgANW1a1erbOODDz5QVRyrqmPHjlmlf1GyAXGqgEyVc8DCajp06IBSii5duhAeHs7p06e1LqlA9erV49ChQ6xcuZKNGzcyadKkIt/GiBEjOHRgP82bNy/yvkXpJXNBiGKh0+l46qmnOHPmjNal3FeHDh3Yvn07D/N7IcSDyFwQQlNTpkzh7NmzlmkQS6o1a9YAEB0drW0holyQABbF4o033gAgKSlJ20IeoHr1vGux586dq3ElojyQABbF4vYlVxcuXNC4ksI5f/681iWIckACWAghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBC/ElCQt694M6dO8epU6c0rkaUdRLAQvy/c+fOYTQaLf9u36Gj3BlDWJUEsCjVLly4QEREBE2aNMHHx4du3bpx4sSJAtsmJSXh7u4OQGxsLN27d8/3eoMGDVi6dBmVHZyoVa8+G75bX+Ct42fMmEFGRsZD1VnQ9oSQABalllKKnj170rFjR06dOsX+/fv56KOPuHjx4iP3GRkZwZGEgyQmxOPp6Vlgm0cJYCEKIgEsSq2YmBgqVKjA4MGDLcuMRiPPPPMMo0aNwt3dHQ8PD5YvX37fftLT0xk4cCB+fn54eXmRkJBAzZo1MZlMvPnmm7i7u2MwGJg1axYzZ87k999/JzAwkMDAQAA2b95M27Zt8fb2pnfv3qSlpQGwadMmWrRogbe3N99++6313ghRatlqXYAQj+qXX37Bx8fnruXffvst8fHxHDp0iNTUVHx9fQkICLhnPx9++CFBQUF8+eWXXLt2DT8/Pzp16sTixYtJSkoiPj4eW1tbrly5Qo0aNZg2bRoxMTHUrFmT1NRUPvjgA6Kjo7G3t+fjjz9m2rRpvPXWW7zyyits3boVV1dXwsPDrflWiFJKAliUOT/99BORkZHo9Xrq1KlDhw4d2LdvHwaDocD2mzdvZt26dUydOhWAzMxMzp49S3R0NIMHD8bWNu/XpEaNGnetu2fPHhITE3n66acByM7Opm3bthw7doxGjRrRtGlTAF566SW++OILa+yuKMUkgEWp5ebmxsqVKx+7H6UUq1atonnz5o+0bnBwMMuWLcu3PD4+/rHrEmWfnAMWpVZQUBBZWVn5jiwTEhKoVq0ay5cvx2QycenSJbZv346fn989++nSpQuzZs2yXHJ28OBBAIKDg5k3bx65ubkAXLlyBQBHR0du3rwJgL+/Pzt37uTXX38F8s4nnzhxghYtWpCUlGS5lvjOgBYCJIBFKabT6Vi9ejXR0dE0adIENzc33n77bV588UUMBgNGo5GgoCA++eQT6tate89+3nnnHXJycjAYDLi5ufHOO+8AMGjQIJ566ilLX0uXLgXg1VdfJSQkhMDAQGrVqsXChQuJjIzEYDBYTj9UqlSJL774gtDQULy9valdu3axvCeidNE9zIXmrVu3VnFxcVYsR5RlOp2OnTt30q5dO61LEaJY6XS6/Uqp1nculyNgIYTQiASwEEJoRAJYlGoODg4lpu/Y2Fh27dplpWpEWSQBLEQRkQAWD0sCWJQ5p06dIiQkBB8fH9q3b8+xY8cAWL9+PW3atMHLy4tOnTpZ5oxIS0tjwIABeHh4YDAYWLVqlaWvsWPHYjQa8ff3t7QvqJ+kpCTmzp3L9OnT8fT0ZMeOHVy6dIkXXngBX19ffH192blzZ/G/GaJkU0oV+uHj46OEeFSA2rlzZ5H2aW9vf9eyoKAgdeLECaWUUnv27FGBgYFKKaWuXLmizGazUkqp+fPnq5EjRyqllHrrrbfU8OHDLetfuXLFUu+6deuUUkqNGjVKvf/++/ftZ8KECWrKlCmWfiIjI9WOHTuUUkqdOXNGtWjRomh2WpQ6QJwqIFNlJJwoU9LS0ti1axe9e/e2LMvKygIgOTmZ8PBwzp8/T3Z2No0aNQIgOjqaqKgoS/vq1asDYGdnZ5lC0sfHhx9//PG+/dwpOjqaxMREy/MbN26QlpZm1fPWonSRABZlitlsplq1agUOBR42bBgjR44kLCyM2NhYJk6ceN++KlSoYJkPWK/XW0bEFbYfs9nMnj17qFSp0uPskijD5BywKFOqVq1Ko0aNWLFiBZB3iu3QoUMAXL9+nfr16wOwaNEiyzrBwcHMmTPH8vzq1av33ca9+vnzEGWAzp07M2vWLMtzmR9C3EkCWJRqGRkZODs7Wx7Tpk1jyZIlLFiwAKPRiJubG2vXrgVg4sSJ9O7dGx8fH2rWrGnpY9y4cVy9ehV3d3eMRiMxMTH33ea9+unRowerV6+2fAk3c+ZM4uLiMBgMtGrVirlz51rnTRCllgxFFsVGhiKL8kqGIgshRAkjASyEEBqRABZCCI1IAIsyadOmTTRv3hxXV1cmT5581+tZWVmEh4fj6upKmzZtSEpKsrz20Ucf4erqSvPmzfnhhx8e2Ofs2bNxdXVFp9ORmppq1f0SZUxBozPu9ZCRcOJxYIWRcAXJzc1VjRs3VqdOnVJZWVnKYDCoI0eO5GszZ84c9fe//10ppdSyZctUnz59lFJKHTlyRBkMBpWZmal+++031bhxY5Wbm3vfPg8cOKBOnz6tGjZsqC5dumT1/ROlD/cYCSdHwMLqlFKYTCYATCYTZrPZqtvbu3cvrq6uNG7cGDs7OyIiIiyXot22du1a+vXrB0CvXr3YsmULSinWrl1LREQEFStWpFGjRri6urJ379779unl5YWLi4tV90mUTRLAwuq+/Oory52FAwICGPDyK1bdXkpKCg0aNLA8d3Z2JiUl5Z5tbG1tcXJy4vLly/dctzB9CvGwJICF1fV64QXaB3bCvpYzzi6NGfv2aK1LEqJEkLkghNU5OTmxdfNGli1bRmhoKDVq1LDq9urXr8+5c+csz5OTky1Dh+9s4+zsTG5uLtevX+eJJ56477oP6lOIhyVHwKJY2Nra8te//tXq4Qvg6+vLyZMnOX36NNnZ2URFRREWFpavTVhYmGUeh5UrVxIUFIROpyMsLIyoqCiysrI4ffo0J0+exM/Pr1B9CvGwJIBFmWNra8vs2bPp0qULLVu2pE+fPri5uTF+/HjWrVsHwMsvv8zly5dxdXVl2rRplsvK3Nzc6NOnD61atSIkJIQ5c+ag1+vv2SfAzJkzcXZ2Jjk5GYPBwKBBgzTbd1G6yFwQQghhZTIXhBBClDASwAKAwMDAfKO+AGbMmMGQIUMKbO/i4vLAUV+TJk3K9/z2LGhJSUm4u7sDEBcXxz/+8Y9HLduif//+rFy58rH7EaI4SQALACIjI/PdlgcgKiqKyMjIR+7zzgAu6I7BrVu3ZubMmY+8DSFKMwlgAeSNBvv+++/Jzs4G8o5Sf//9d1JSUvDw8MDd3Z3Rowu+fvf555/Hx8cHNzc3vvjiCwDGjBnDrVu38PT0pG/fvgAF3gstNjbWct+1bt264enpiaenJ05OTixatAiTycSoUaPw9fXFYDAwb948IG903euvv07z5s3p1KkTf/zxR5G/J0JYXUHjk+/1kLkgyrbQ0FC1Zs0apZRSH330kRowYIBq0KCB+uOPP1ROTo4KDAxUq1evVkqpfPMeXL58WSmlVEZGhnJzc1OpqalKqbvvWHz7+enTp5Wbm5tSSqmYmBgVGhqar11cXJzy8PBQ165dU/PmzbPcjTgzM1P5+Pio3377Ta1atUp16tRJ5ebmqpSUFOXk5KRWrFhhhXdFiMeHzAUhHuTPpyGioqJo2LAhHTt2pFatWtja2tK3b1+2b99+13ozZ87EaDTi7+/PuXPnOHny5CPXkJqayl//+leWLl2Kk5MTmzdvZvHixXh6etKmTRsuX77MyZMn2b59O5GRkej1ep588kmCgoIeeZtCaEUCWFg899xzbNmyhQMHDpCRkYGnp+cD14mNjSU6Oprdu3dz6NAhvLy8yMzMfKTtm0wmIiIiGD9+vOVLOqUUs2bNIj4+nvj4eE6fPk3nzp0fqX8hShoJYGHh4OBAYGAgAwcOJDIyEj8/P7Zt20Zqaiomk4lly5bRoUOHfOtcv36d6tWrU6VKFY4dO8aePXssr1WoUIGcnJxCb3/MmDEYDAYiIiIsy7p06cLnn39u6efEiROkp6cTEBDA8uXLMZlMnD9//oE30hSiJJK5IEQ+kZGR9OzZk6ioKOrVq8fkyZMJDAxEKUVoaCjPPfdcvvYhISHMnTuXli1b0rx5c/z9/S2vvfrqqxgMBry9vVmyZMkDtz116lTc3NwsR97vvfcegwYNIikpCW9vb5RS1KpVizVr1tCzZ0+2bt1Kq1ateOqpp2jbtm2Rvg9CFAcZCSeEEFYmI+GEEKKEkQAWQgiNSACLYvGgm2QKUR5JAAurM5lMvPbaa2zcuJHExESWLVtGYmKi1mUJoTkJYGF1hblJphDlkQSwsDq5oaUQBZMAFkIIjUgAC6srzE0yhSiPJICF1ckNLYUomAxFFlb35xtamkwmBg4caLmhpRDlmQxFFkIIK5OhyEIIUcJIAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNPNQdMXQ63SXgjPXKEUKIMqmhUqrWnQsfKoCFEEIUHTkFIYQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAa+T/fXx/DEkp+IgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "Landfill.plot_sankey()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABSy0lEQVR4nO3deZyN5f/H8dc1m2XGmiUhFJFhZhh7EbIvU8qaFkkLrRT1TUWrkihKSvrKN1sqWUNEEbKORCLZSXZmxpjt+v1xxvnNmDEGM3PPOfN+Ph7nwbnPfV/35z5nzud87vu+7vsy1lpERERERDyBj9MBiIiIiIhklopXEREREfEYKl5FRERExGOoeBURERERj6HiVUREREQ8hopXEREREfEYKl6zgTHGGmMqOx1HZhhjrjfGRBljfJ2ORUTECcrZIp4lTxWvyV/4848kY8zZFM97XmSZpsaY/VkYwzJjTJ+sau9qWWv3WmuDrLWJWd22MWaiMSYu+f09boz5wRhTLfm1ocaY+As+k0Eplu1gjFljjIk2xhwzxkw2xpRLZx1Dk3946l9mbF2NMSuNMTHGmGWZmP8eY8ye5Hi+M8YUT/Fa1AWPRGPMmOTXmib/rV04T8PLjDefMeZzY8xpY8w/xpgBKV6rboxZZ4w5kfxYbIypfjnti+RGytlpeWrONsYEGGPeM8bsT152tzHm/cuILcwYs9wYcyq5jZczudyS5N8IvxTTKhpjlibn/23GmBYpXktvO09mNs5MruMBY8z65Hy+3xgzPGV8cml5qnhN/sIHWWuDgL1AxxTTJjsdn5canvx+lwP+BSameG16ys/EWjscwBjTGZgCvA+UAIKBc8AKY0yx8wsbYwxwP3A8+d/LcTy5/bcvNaMxJhj4BLgPKA3EAGPPv37B39W1wFlgRoomDl6wnUHW2lXprKeXMWbiRcIYClQBKgDNgEHGmDbn2wc6A8VxvV+zgWmX2i6R3E452xHZlbP/A9QB6gGFgKbAhsuIawrwM648dxvQzxgTkdECyTs4/um8NBXYCFwDDAa+NsaUzGA7i16k/YnGmF4XWX1G6ygIPIPrvaoP3A48l9G2SGp5qni9mOSjWu8bYw4mP95PnhYIfA9cl2IP7DpjTD1jzCpjzEljzCFjzIfGmIArWG/T5L2uZ40x/ya39WCK1wsk76nuSd7bXGGMKZD8WoQxZktyDMuMMTenWG63MWagMea35L3gCcaY0saY740xZ5KPzBVLnrdiyr3S5LZeN8b8kjzvImNMiRRt358czzFjzMvJ62rBJVhrY3AlnxqXeE8M8B7whrV2irX2rLX2H6APEAX0TzF7Y6AM8BTQ/XI+A2vtYmvtV7gKv0vpCcyx1v5srY0CXgbuMsYUSmfeu3El/OWZjSWTHgBet9aesNb+AYwHegFYa09aa3db13B5BkgELnoK9MLPLPlIw5cpnl/RZyySU5SzPTJn1wVmWmsPWpfd1tpJmXrjXSoCk621idbancAKXEXyxeIqAgwBBl0w/SagNjAkOdZvgM24cneWuNQ6rLUfW2uXW2vjrLUHgMnALRdpK9XnnTzNfTbAGOOb/Dd31BizyxjzxIXzeyMVry6DgQZAGBCKa8/wJWttNNCW1EfODuIqDvrj2mtqiGuvqd8VrvtaoAhQFngI+Mj8/57qCCAcaIRrb3MQkJT8xZiKa8+tJDAfmHNBMr4baAncBHTEldBfTJ7fB1fBdzH3AA8CpYAAkvcIjetU9FhcxVyZFHFfkjEmKHm5jZeYtSpwPamPXGKtTQK+Sd6m8x4A5gBfJT/vmGJ99xhjfstMbJkQDGxKEctOIA7Xe3uhB4BJNgvHXU7+eyiTMobk/wdfMN9JIBYYA7x1heu64s9YJAcpZ6fmCTl7NTDAGNPPGFMzuehNub6xxpixXNz7wP3GGH9jTFVcn+PiDOZ/C/gY+OeC6cHA39baMymmpcmnV+ly19EE2HKF63oY1998GK6C+c4rbMejqHh16Qm8Zq3911p7BHgV1ynidFlr11trV1trE6y1u3GdUr7tCtcdn7zueGvtfFx7qlWNMT5Ab+Bpa+2B5L3Nldbac0A3YJ619gdrbTyuhFkAV8I8b4y19nDyXt1y4Fdr7UZrbSwwE6iVQUz/tdZut9aexVUYhiVP74zrCOQKa20c8ApwqSLtueSi6i8giOSjhcm6Jh+FOP+4DtePC8ChdNo6dP51Y0xBoAswJfk9+JoUXQeSjwCEXCK2zAoCTl0w7RSuU19uxpgKuP4Ovrhg3usu2M6TyUeILmf959d50fUnn9oqAjzBpX9wLuZKPmORnKacnVquz9nAMOAdXJ/dOuCAMeaB8zNaa/tZazPaoZibvD1ngW3ABGvt2vRmNMbUwXUkc0w6L2cmn1+4nUsziCs9mfrNSI61N67uFCMucx3ndQU+sNbut9aeIBNd4byBileX64A9KZ7vSZ6WLmPMTcaYucZ14cxpXHt4JS42/yUcs9YmpHgeg+sPvwSQH9h5qXiT93D3kXqP+nCK/59N53kQF5dyT/V8POfXuy/FemOAYxm0AzDCWlvUWnuttTYi+ajleV8lv3b+cRA4mvxamXTaKpPi9U5AAq4jGOA67dLWpO63lFWigMIXTCsMnLlg2n3ACmvtrgumH7xgO4smHyE6f7ThZPKPxVjgnhQJ8/yR46gU68xo/SS3Ow6YZIwpdZnbCVf2GYvkNOXs1HJ9zk4u5j+y1t4CFAXeBD5P2X3iYozrAtkFwGu43uPyQGtjTJpiN3knYiyunYiEC18nc/n8wu1slqL931Lk7HuAsSly9vkjx5n6zTDG3ImrqG9rrT3KlUn1GV/wf6+l4tXlIK4LYc67nv/vC5neXurHuPb8qlhrC+M6tWPSme9qHMV1CvjGdF5LFW/y6ZfywIEsjuFCh3B14j+/3gK4OqNnpT+B/biOqrolJ6S7gSXJkx7AlaD3GmP+wXXKyh9XMslqW3Cdmjwfyw1APmD7BfPdT9qjrhlKPtpQNPmoaT9cR5LPJ8yQ5HlO4HrvQ1MsGsrFTzP54LogIKPTgwVS/L9Iiv/nxGcscrWUszMnN+VsN+vqB/oRcALIzJ1RbgASrbWTko+e78d1UWq7dOYtjOtI5vTk34bzR2f3G2Ma48qbN5jU1yxklE8vjD0kRc6eAvRLkbPPF9OXXIdxXXA7HtdFiJszsepM5Wxcf1deT8Wry1TgJWNMyeSO7q8A5y9gOQxck9z5+7xCwGkgyrhuI9I3qwNK3jP/HBhpXBcc+BpjGhpj8uE6LdTeGHO7McYfeBbXlZ0rszqOC3wNdDTGNEruqzWULP4BsNZaXP21Xkrut5rfGHMt8BmupDTKGFMWV5+1DrhOj4XhSgzvkMm7DiS/n/kBP8AneT3pXZUKrqO6HY0xjZNP978GfJuyP5MxphGuYnHGRdq4WpNwvSfFkv/mHib5KmBjTEtjTK3kbSoMjMT1o/BHBu3da1wXuFTC9QNTKHn7s/0zFskCytmZkytyNoAx5hnjuuCtgDHGL7nLQCEy18Vpu6sJc48xxie5/W5Aetc1nMJ1NDIs+XG+wA3H1RVjOxAJDEmOtRMQgqt/bpa41DqMMc1x/a7cba1dk8lmeydve8vkts7n7K+Ap40xZY0xRYHns2o7cjMVry5v4OqD8xuuKwI3JE/DWrsNV6L8O0Ufn+dwHeE7g2vPaXo2xfVccjxrcd3a6R3Ax1r7J3Avrv48R3F17u+Y3Kcp21hrtwBP4trjPYTr1Mi/uJJwVq5nOq5T8P1xneLaimuv8xZr7bHk1yKttYustf+cfwCjgRBjTA1jTE9jTEZ70vfhOhX3Ma67FpzF9VkC7vtLNk6x3Y/hSjb/4kq4F56ueoALCtoUUl75fP5xuVe2DsF1OnIP8BPwrrV2QfJrRXH9jZ5KnudGoE1yX7mLiQJ24boo5C1c78F9OfUZi1wl5exMyEU5G1zdGd7D1cXhKPA4ruLtbwBjzDhjzLiLtH8auCu5/RO4CsPfSf7Mzf8P3HC9dUn5u3AkuZnDKd7v7riOzp7vI9rZuvpOn9ctnZx9ud2wMlrHy7iOns5P0f73l2jvxuRteR5XH+9XcF3cNx5YhOu7sBFXV7oEXBcpei1js+6iaMljjOtq1JO4TsVd2M9TciljzG6gj7U2oyt1z8+rz1jES+j77HmMMRVxHWjwv0gf3gvnbwuMs9ZWuNS8nkxHXuWyGGM6GmMKJp8+H4HrKMNuZ6OSrKTPWMR76Pvs3ZK7YbRL7opRFtdZuplOx5XdVLzK5boD18UHB3GN+NTd6vC9t9FnLOI99H32bgZXN4ITuLoN/IGrS4FXU7cBEREREfEYOvIqIiIiIh7jsopXY8yCS88lIpL7KH+JiHgHv8uZuXDhwq3r1KmjfgYi4olOOx1ATitRooStWLGi02GIiFy29evXH7XWpjtq5mUVr1WqVGHdunVZE5WISA4yxuxwOoacVrFiReVsEfFIxpg9F3tNfV5FRERExGOoeBURERERj6HiVUREREQ8hopXEREREfEYKl5FRERExGOoeBURERERj6HiVUREREQ8hopXEREREfEYKl69hK+vL2FhYe7H7t27r6idRo0aAbB7926mTJninr5u3TqeeuqprAhVRCTPU84WuXLG2syP9lqnTh2r0Vpyp6CgIKKiorKsvWXLljFixAjmzp2bZW2KOMkYs95aW8fpOHKScnbupZwtkrGMcraOvHqxihUrcvToUcC1F960aVMAjhw5QsuWLQkODqZPnz5UqFDBPV9QUBAAL7zwAsuXLycsLIxRo0axbNkyOnToQFJSEhUrVuTkyZPu9VSpUoXDhw9z5MgR7r77burWrUvdunX55ZdfcnR7RUQ8mXK2SOaoePUSZ8+edZ9+6tSpU4bzvvrqqzRv3pwtW7bQuXNn9u7dm2aet99+m8aNGxMZGUn//v3d0318fLjjjjuYOXMmAL/++isVKlSgdOnSPP300/Tv35+1a9fyzTff0KdPn6zdSBERL6GcLXLl/JwOQLJGgQIFiIyMzNS8K1ascCeyNm3aUKxYsctaV7du3Xjttdd48MEHmTZtGt26dQNg8eLFbN261T3f6dOniYqKch8ZEBERF+VskSun4tWL+fn5kZSUBEBsbGyWtduwYUP++usvjhw5wnfffcdLL70EQFJSEqtXryZ//vxZti4RkbxCOVskc9RtwItVrFiR9evXA/DNN9+4p99yyy189dVXACxatIgTJ06kWbZQoUKcOXMm3XaNMXTq1IkBAwZw8803c8011wDQqlUrxowZ454vs0cVREREOVsks1S8erEhQ4bw9NNPU6dOHXx9fVNNX7RoETVq1GDGjBlce+21FCpUKNWyISEh+Pr6EhoayqhRo9K03a1bN7788kv36SeA0aNHs27dOkJCQqhevTrjxo3Lvo0TEfEyytkimaNbZeVB586dw9fXFz8/P1atWkXfvn21xy1eT7fKEk+lnC15UUY5W31e86C9e/fStWtXkpKSCAgIYPz48U6HJCIiF6GcLZKaitc8qEqVKmzcuNHpMEREJBOUs0VSU59XD5adtzO53LaXLVvGypUrsykaERHPp5wtkjVUvEqWUCIUEfEcytniyVS8epmdO3fSpk0bwsPDady4Mdu2bQNgzpw51K9fn1q1atGiRQsOHz4MQFRUFA8++CA1a9YkJCQk1e1ZBg8eTGhoKA0aNHDPn147u3fvZty4cYwaNYqwsDCWL1+uYQdFRDJBOVvkClhrM/0IDw+3knsEBgammda8eXO7fft2a621q1evts2aNbPWWnv8+HGblJRkrbV2/PjxdsCAAdZaawcNGmSffvpp9/LHjx+31loL2NmzZ1trrR04cKB9/fXXM2xnyJAh9t1333W306NHD7t8+XJrrbV79uyx1apVy5qNFrlCwDp7GfnOGx7K2bmLcrZI5mWUs3XBlheJiopi5cqVdOnSxT3t3LlzAOzfv59u3bpx6NAh4uLiqFSpEuAaHnDatGnu+c8POxgQEECHDh0ACA8P54cffsiwnQtp2EERkYwpZ4tcGRWvXiQpKYmiRYume/+/J598kgEDBhAREcGyZcsYOnRohm35+/tjjAHA19eXhISEy2pHww6KiGRMOVvkyqjPqxcpXLgwlSpVYsaMGYCrS8imTZsAOHXqFGXLlgXgiy++cC/TsmVLPvroI/fz9IYdTOli7Vw4NKGGHRQRyZhytsiVUfHqwWJiYihXrpz7MXLkSCZPnsyECRMIDQ0lODiYWbNmATB06FC6dOlCeHg4JUqUcLfx0ksvceLECWrUqEFoaChLly7NcJ0Xa6djx47MnDnT3flfww6KiKSmnC2SNTQ8rIjkCRoeVkTEc2SUs3XkVUREREQ8hopXEREREfEYKl5FRERExGOoeM0DEhMTqVWrlvsegE2bNqVq1aqEhYURFhZG586d011u/fr11KxZk8qVK/PUU09xvn90t27d3MtWrFiRsLCwnNoUERGvVrFiRWrWrElYWBh16ri6+/Xq1YtKlSq5826jRo3SXXbXrl3Ur1+fypUr061bN+Li4gDYu3cvzZo1o1atWoSEhDB//vwc2x6R7KD7vOYBH3zwATfffDOnT592T5s8ebI7MV5M3759GT9+PPXr16ddu3YsWLCAtm3bMn36dPc8zz77LEWKFMm22EVE8pqlS5emujMAwLvvvnvRAw3nPf/88/Tv35/u3bvz2GOPMWHCBPr27csbb7xB165d6du3L1u3bqVdu3bs3r07G7dAJHvpyKuX279/P/PmzaNPnz6XtdyhQ4c4ffo0DRo0wBjD/fffz3fffZdqHmstX331FT169MjCiEVE5HJZa/nxxx/dBe4DDzzgztnGGPfBi1OnTnHdddc5FaZIllDx6uWeeeYZhg8fjo9P6o+6Z8+e7lNQAwcOTLPcgQMHKFeunPt5uXLlOHDgQKp5li9fTunSpalSpUr2BC8ikscYY2jVqhXh4eF8+umn7ukDBw505+yePXumWe7YsWMULVoUPz/XCdWUOXvo0KF8+eWXlCtXjnbt2qUajEDEE6nbgBebO3cupUqVIjw8nGXLlqV6LTPdBi5l6tSpOuoqIpKFVqxYQdmyZfn3339p2bIl1apVAzLXbeBipk6dSq9evXj22WdZtWoV9913H7///nuagxoinkLFqxf75ZdfmD17NvPnzyc2NpbTp09z7733pjtvYmIi4eHhAERERNC3b1/279/vfn3//v3uIQYBEhIS+Pbbb1m/fn32boSISB5yPs+WKlWKTp06sWbNmovO27p1aw4fPkydOnUYP348J0+eJCEhAT8/v1Q5e8KECSxYsACAhg0bEhsby9GjRylVqlT2b5BINtBulxcbNmwY+/fvZ/fu3UybNo3mzZvz5Zdfpjuvr68vkZGRREZG8tprr1GmTBkKFy7M6tWrsdYyadIk7rjjDvf8ixcvplq1aqm6FoiIyJWLjo7mzJkz7v8vWrSIGjVqXHT+hQsXEhkZyWeffYYxhmbNmvH1118D8MUXX7hz9vXXX8+SJUsA+OOPP4iNjaVkyZLZvDUi2UfFax6Vss9rixYt0p1n7Nix9OnTh8qVK3PjjTfStm1b92vTpk1TlwERkSx0+PBhbr31VkJDQ6lXrx7t27enTZs2QOo+r2FhYe7bYKX0zjvvMHLkSCpXrsyxY8d46KGHAHjvvfcYP348oaGh9OjRg4kTJ2KMydFtE8lK5vy9OzND42SLiKfKaJxsb6WcLSKeKqOcrSOvIiIiIuIxVLyKiIiIiMdQ8SoiIiIiHkPFq7jFx8czceJEjh075nQoIiKSCfPmzWPLli1OhyGSo1S8CgAnTpygaYvWPPHCEEJq12Hbtm1OhyQiIheRmJjIU/2fpdsDD9Pg1ibMnj3b6ZBEcoyKVyE6Opqw8Hr8um4DiQVLcvQs1K5bj99//93p0EREJB1detzLh6PfJ7FwGRKKlOeOO+9k4qT/OR2WSI7QCFtCwYIF6d3rfmJiYhj+ztt0jIjg5purU7p0aadDExGRdES0b8ONla5nxPDhADw3aBDVk4eSFfF2us+rpGKMYdWqVTRo0MDpUESylO7zKt7IGEP79u2ZO3eu06GIZCnd51VEREREvIKKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeJY2AgACnQxARkUxSzpa8RsWruCUmJgJQsWJFZwMREZFMU86WvEbFq7i9+OKLABQvXtzhSERE5FL+/vtvAB566CGHIxHJWSpexW348OGEh4c7HYaIiGRC06ZNAQgODnY2EJEcpuJVAJg3b16qf0VEJPc6ceIE+/bt491333U6FJEcp+JV2LJlCx06dKBWrVqULl3a6XC8TkxMDHv37nU6DBHxErGxse7uXc8++6zD0XifxMRE/vrrL6fDkAyoeM3jPvvsM2rUqIGvry8bNmxwOhyv9PIrQwitXZd9+/Y5HYqIeLiNGzdSoEABAE6dOoUxxuGIvM+8efOoVu1mlixZ4nQochF+TgeQleLj47HWZmv70dHRV9VGZhJNRvMUKFAAX1/fK17/6dOn+fPPP3n11VfdX8wGDRqwdOlSYmNjr7hdubiExASirR8t2rRn/a8rCQoKcjokkVwhISHBfZeT7GCt5cyZM9nWfmbky5ePfPnyXfHycXFx7Nixg6lTpzJixAj39BMnTpAvXz7OnTuXFWFKCvHx8QQUKcFdXbqxdvVKbrrpJqdDkgt4TfFqreX6ijdw+J+D+PhceXGXkcSE+Gxp93L5+V/5Pf0S4uPSTFu3fgOFChe5mpDkEoq2G8CxvRvp3O0e5s/5Dh8fnfQQ6dKjJ7O+mYGPb/b8FOWWnO3r53/Fy6a3Db5+/pQoWepqQpJLuKZeR5KKlKVF63Zs2rCWYsWKOR2SpOA1xasxhieeeJx3x35O4S7D8AnI73RIIqnYKvVZ8+1Qnhv0PCNH6CILkWee6MeiHxZTtMub+F9T3ulwRNKIOnGAdhGd+PnHH/D3v/KdEMlaXnX458UXnqdd04ZELxyFtUlOhyOSivH1J7D9ID6bNJX//nei0+GIOO62227j/RHvcGb2mySePe10OCJpBDZ+gO1Hz/HIY/2ytVuiXB6vKl6NMUycMJ5KhSzRK750OhyRNHwLFiEoYjBP9n+WFStWOB2OiOMe7tOHXvd0JWreO9jE3HGaX+Q84+NLYJsBfLvgR0aOet/pcCSZVxWv4BrjecHcWQTs+5Xo33WloOQ+ASWuJ7DV03S88y73CDkiedmoEcOpU7ks0T9+oqNbkuv45CtIoYjBDHljGPPnz3c6HMELi1eAEiVKsHjBfM798gWx+7c4HY5IGgVuCMev9t3c3rodp06dcjocEUf5+vry7YxpFIvZR/T6WU6HI5KGX5HSBLUbRPee97Fli+oKp3ll8QpQvXp1pk/5kqh57xJ/8h+nwxFJI7B2B84Uu4k77upCQkKC0+GIOKpQoUIsXjCPpE2ziflrjdPhiKSRv9zN5Gv8IC3atOPIkSNOh5OneW3xCtC2bVteH/IyUXPeJOlcjNPhiKQR2LQPv+07zpNPD3A6FBHHVahQgflzZhHzwxji/t3ldDgiaQRWb0Z8xUa0bh+he+w6yKuLV4Bnnn6Szu1aEvX9CGxS9t0MW+RKGB9fgtoNZMq3sxj78cdOhyPiuAYNGjDuozGcmfMWidEnnA5HJI3ARvewJ8af+x98SH20HeL1xasxhnFjP+Tm0oFE//xfp8MRScMnfxBBES8x6MWXWbx4sdPhiDju3nt78vjDvYma+zY2Ie3AKiJOMsaHwFZPs3DFWt54a5jT4eRJXl+8Avj7+zP3u28oeHgz0ZsWOB2OSBr+xa4jqO2z3N21O3/++afT4Yg47q03XuPWsKpELf5IR7ck1/EJyE9QxxcZPnI03377rdPh5Dl5ongFKFasGEsWfU/8mmmc3bPJ6XBE0sh/fQgBDXrSonU7jh8/7nQ4Io7y8fFh+uT/ca09TvSvM5wORyQNv0IlCOrwAvf37sPGjRudDidPyTPFK0CVKlX47psZRH//HvHHDzgdjkgagSGtOFumFu0jOhEfrxu2S95WsGBBfvh+LmbbD0Rv06AekvvkK1OFAk0fpVW7Dhw6dMjpcPKMPFW8AjRt2pT3hg/jzOw3SDx7xulwRNIIbHw/fx6Po8+jfXW6VPK86667joXz53J22SecO7TD6XBE0gisdiu26u20bNuBs2fPOh1OnpDnileARx95hPu73U30/OHYRN1fU3IX4+NLUOsBzFq0jPdGjnI6HBHH1a5dmy8mjCdq7tsknDnqdDgiaQTW78ohitGt530kJSU5HY7Xy5PFK8AHI0dQq1JpopZ+qqNbkuv45CtIUMfBDH3zbebNm+d0OCKOu/vuuxnU/0mi5gwjKT7W6XBEUjHGENTicZZv3Mbgl4c4HY7Xy7PFq6+vLzO/nk7RM7uI2TDH6XBE0vArUopC7QfR/d77+f33350OR8RxLw1+kVa3hBO98AOs1dEtyV2MXwBBHV7gw08nMHnyFKfD8Wp5tngFKFy4MEsWzicx8jvO7lzndDgiaeQrezP5k4cj/Pfff50OR8RRxhj+N/Fzri8QR/RKFQeS+/gGFqNQxxd5tN8TrF692ulwvFaeLl4BKlasyNxZM4n+YTRxR3Y7HY5IGoHVm5FQ8RbadLhDwxFKnpcvXz4WzpuN365fiN6y1OlwRNIIKFWJgi2foF3Enezdu9fpcLxSni9eARo1asTY0e+7hiOMOeV0OCJpBN5yD3vOBnBfr96X3Ufb19eXsLAwatSoQceOHTl58mT2BJmsT58+bN26NcN5xo0bx6RJk7I1DvFepUqVYvGC+cQu/5xzB/5wOhyRNApWro9PSAdatGlPVFTUZS2rnH1p5nJ+COvUqWPXrfPe0+sDn/8Pn82YR+G7XsP4+TsdjkgqSfGxnJ4xmGcfuY+XXxqc6eWCgoLcyfOBBx7gpptuYvDgzC/vLYwx6621dZyOIyd5e86eN28e3e57kKLd3sGvSCmnwxFJxVpL9OKPqF3Kl+/nzsLX1zdTyylnu2SUs3XkNYV3hr1Joxo3ErVEwxFK7uPjn5+gjv9h+KgxfPPNN1fURsOGDTlwwDVAx86dO2nTpg3h4eE0btyYbdu2AdCrVy/69u1LgwYNuOGGG1i2bBm9e/fm5ptvplevXu62+vbtS506dQgODmbIkP+/urZp06acL5iCgoIYPHgwoaGhNGjQgMOHDwMwdOhQRowY4Z7/+eefp169etx0000sX74cgJiYGLp27Ur16tXp1KkT9evXx5sLMbl87du3Z+jgFzgz+w2SzsU4HY5IKsYYAps/ytod+xkwcNAVtaGcnT4Vryn4+Pjw1dQvKZVwhOi1V1YciGQnv0IlCOr4Hx546GE2bNhwWcsmJiayZMkSIiIiAHjkkUcYM2YM69evZ8SIEfTr188974kTJ1i1ahWjRo0iIiKC/v37s2XLFjZv3kxkZCQAb775JuvWreO3337jp59+4rfffkuzzujoaBo0aMCmTZto0qQJ48ePTze2hIQE1qxZw/vvv8+rr74KwNixYylWrBhbt27l9ddfZ/369Ze1vZI3PDugP3e2akrUgvewSYlOhyOSivH1J6j98/x38ld8NmHCZS2rnH1xKl4vEBgYyOIF82DrQmK2r3Q6HJE08l1bmQJNH6V1+44cPHjwkvOfPXuWsLAwrr32Wg4fPkzLli2Jiopi5cqVdOnShbCwMB599NFUQxt27NgRYww1a9akdOnS1KxZEx8fH4KDg9m9ezcAX331FbVr16ZWrVps2bIl3T5TAQEBdOjQAYDw8HD3she666670syzYsUKunfvDkCNGjUICQnJ7FskeYgxhs8+HUfVa/IRvfwLp8MRScO3QGEKRQzm6WcH8dNPP11yfuXsS1Pxmo6yZcuyYO5sYn78mLjDO50ORyQN13CELWjZtgMxMRmfLi1QoACRkZHs2bMHay0fffQRSUlJFC1alMjISPfjjz/+/8KXfPnyAa6zEef/f/55QkICu3btYsSIESxZsoTffvuN9u3bExub9sbx/v7+GGMA10UICQnpj2h3fh0ZzSNyMf7+/sybPZMChyKJ/m2h0+GIpOF/TXmCWvfnjrs6s3NnxnWFcvalqXi9iDp16vD5p59wZs5bJEQddzockTQC63fhH1M808MRFixYkNGjR/Pee+9RsGBBKlWqxIwZMwDXhQWbNm3K9LpPnz5NYGAgRYoU4fDhw3z//fdXvB0Xc8stt/DVV18BsHXrVjZv3pzl6xDvUbx4cZYs+p641VOI3Zv2dKiI0wpUqoVfnS60aN2OU6cufWcj5eyLU/Gaga5duzDgyX5EzxlGUrzurym5y/nhCFds2s6LL72SqWVq1apFSEgIU6dOZfLkyUyYMIHQ0FCCg4OZNWtWptcdGhpKrVq1qFatGvfccw+33HLLlW7GRfXr148jR45QvXp1XnrpJYKDgylSpEiWr0e8x0033cQ3X00j6vv3iD9x6S41IjktqFZ7ThevSoc7787UEUvl7PTpVlmXYK2lS/d7WPrHIYLaPosxqvcld0mMPsnJaQP5eNS73HtvT6fDyTKJiYnEx8eTP39+du7cSYsWLfjzzz8JCAi4ovZ0q6y8Y+zHH/Of196mcLfh+OYPcjockVRsUiJnZr1O5+Z1+fTjsU6Hk2VyMmerErsEYwxffvFfygecJXrVNKfDEUnDN7AohSIG89gTT3nVcIQxMTHceuuthIaG0qlTJ8aOHXvFSVDyln59+9KjUwTR89/FJqoPteQuxseXoLbPMf27+Yz58COnw8kyOZmzdeQ1kw4fPkxI7TokhfcgsPptTocjkkbMX2tI+GkckevXcv311zsdTq6jI695S0JCAi1at+P3M/kJav6o0+GIpBF/4hCnvnqBb6dPoVWrVk6Hk+voyGsWKF26ND98P4/Y5RM4d/BPp8MRSaNg5Xr4hEZwe+t2nDlzxulwRBzl5+fHrG9nUOjEdqI3zHU6HJE0/IuVIajdQDp3u8c94IBkTo4feU057FlOeuWVV2jSpAktWrS4qnZmz57NnZ064Vu0DMbXH9/AYpRoPwDfoGLs/7g3PgEFwMe1T1C8ZT/yl7uZuCN7OL74ExKjjoFNIjC4OUUadXffjgLg4H+fxL94OUre8fxVxSd5m7WWqMUfEX6ZwxFeiUaNGrFypefcC1lHXq+cJ+ftXbt2UbtufU6fjce3YBEw5qrzdmL0CY59P5qE00cgKRHfIqUo3eXVrNpsyWOif1tEvi1z2LRhLddcc022rcebcnaeKV6zUrHixUnMX5Si3YdzatVX2PhYird4lP0f96bMA6NcCTJZUvw5Dn3+OMVb9aNApdokxcdyZOYwClauS6HarhsBxx/dx5FZb5MUe4brHv4Un4D8Tm2aeAGbGM+Zma9yX0RzPhj5ntPh5BoqXq+cp+ftFStW0KTJbZS6523yl6vOiZ++uKq8fWzBh/iXKE/hOncAEPfvLgJKVXJq88QLRP38X27gMCuWLVHf/mS5stvAsmXLaNq0KZ07d6ZatWr07NmT84X02rVradSoEaGhodSrV48zZ84QGxvLgw8+SM2aNalVqxZLly4FYOLEidx55520bNmSihUr8uGHHzJy5Ehq1apFgwYNOH7cdY/WXr168fXXXwNQsWJFhgwZQu3atalZs6b7cP2RI0do2bIlwcHB9OnThwoVKnD06NE0sRcpXJhWjesTvWAU+cpVz/CWLNFbfyJf2ZspUKk24BqfvnjLxzi1+uv/n+ePnwis0Yz8FWtz9i/vueBGnGF8/QlsP+iKhiO8HEFBrqu4M/oui3fx1Lx96623Urx4MaIWjCQx5hT5y9e4qrydGH0c30Il3POrcJWrFXjr/ew8mUTvhx/NtvzpTTnb0T6vGzdu5P3332fr1q38/fff/PLLL8TFxdGtWzc++OADNm3axOLFiylQoAAfffQRxhg2b97M1KlTeeCBB9yjQ/z+++98++23rF27lsGDB1OwYEE2btxIw4YNmTRpUrrrLlGiBBs2bKBv376MGDECgFdffZXmzZuzZcsWOnfuzN69ey8a++hR73FjUR9OL/+SgJIV3dMPT32Rg/99kkOTBgAQf3QPAddWTrWsf7Ey2LizJJ1zjYwUvW05gTc3IbB6E6K3XnroOJFLOT8c4TOZHI7waqX3XRbv5Kl5OygoiPu63kX0vOHE7Fh9VXm7UK32HPt+NP9M/Q+nVk4n4cyxq31bJY8zPr4Etn6GuUtWMDz5bzs7eXrOdrR4rVevHuXKlcPHx4ewsDB2797Nn3/+SZkyZahbty4AhQsXxs/PjxUrVnDvvfcCUK1aNSpUqMD27dsBaNasGYUKFaJkyZIUKVKEjh07AlCzZs0rHpe3TZs2FCtW7KKxt27dmrNnTuJz5h/8Cpd0Ty/d4y2ue3AMZe4fman34NyhHfgWKIxf4VLkrxBK3OG/STyri23k6vlfU57ANgMyNRzh1UrvuyzeyZPz9orlP5NwZBdxezZRuH5n9/TLzdsFbgin7KOfUSi0NfHH93No4tMkxlx6xCSRjPjkK0hQxxd5fdhwZs+ena3r8vSc7WjxmnL83asZH/fCcXxTjvGbXePyLl26lM2bN7Nh/Tri180gdt/v6c7nX+J64v75K9W0+JP/YAIK4JOvINF//ET8sf3s/7g3Bz55mKS4GGL+9Kw9IMm9ClQMw79uV25v1TZTwxFeqaz6Lkvu5+l5++DBA1Qscw2xv/9w0fkulbcBfAsUIrB6U0p0eJZ8Zapc9DdA5HL4FSlFofYv0PP+Xvz2W/YNc+zpOTvX3SqratWqHDp0iLVr1wJw5swZEhISaNy4MZMnTwZg+/bt7N27l6pVq2bpulOOy7to0SJOnDhxyWWqVavG19OmEDX/XUhKTPN6YPWmnNu/lbO7IwHXhQAnFn9C4fp3Y20SMdtWUKb3h5Tr+znl+n5OqbteJvqPn7N0uyRvCwxrR1SJmzM9HKHI5fKkvB0UFMSShfNJ+n0+MTt+TXeejPI2wNk9m0iKd3V/SDoXQ/yJQ6nOwIlcjXxlq5G/yUO0bNuew4cPOx1OrpTriteAgACmT5/Ok08+SWhoKC1btiQ2NpZ+/fqRlJREzZo16datGxMnTky155AVhgwZwqJFi6hRowYzZszg2muvpVChQpdcrlWrVgx7/VWSzp4kKS4m1Ws+/vkoeffLnFo5jQPjH+XQ508QcG0VCtXuwLl9W/ANuga/Qv9/a4x85YOJP7aXhKjjWbptkrcF3vYQWw6ept+TTzkdinghT8vb5cqVY8Hc2cQs+QiS0u7QZZS3AeL++Yt/vujPwc+f4J8vn6NQaCvylbkpS7dL8rbA6k1JvKExrdtHuPuJy//TCFspnDt3Dl9fX/z8/Fi1ahV9+/YlMjIy08s/0rcfX/+4jkJ3vITxyb77a4pciaRz0Zya/jxvvPgcTz7xuNPh5DjdKss7XU3enjZtOg8/8QxFu7+Lb9DF+8qKOMHaJKLmj6DpzWX4evrUVPeGzwty5a2ycqO9e/dSt25dQkNDeeqppxg/fvxlLT92zGiCrytM9PKJ2ROgyFXwyRdIUMfB/OflIe5bFol4uqvJ2927d+Ppfo8QNXcYNp1uXyJOMsaHwFZPseTXTbw57G2nw8lV/JwOIDepUqUKGzduvOLl/fz8CKkRzLalG7IwKpGs41uwCD75C7Fv/36nQxHJElebt2uHhRI3ZqzrmgWdMZNcxvjlw6dIGbb/9bfToeQqOvKahT4dP57/zZhJ/tvz3ilZyf1sUiJRC97j7vYtuS/59kUiednGjRu5v3cfCkcMxvhpVCPJfaJ/nUEZc5JxH37gdCi5Sp4uXnv37k2pUqWoUaPGVbe1bNkyBgx8gUIRL+Fb4NIXeYnktOifJ1KtZAE+/fijVH2nzo+6kpWGDh3qvon8tm3bCAsLo1atWuzcuZNGjRoBsHv3bvd3b9myZXTo0CHL4xDvkpU5+9ChQ7Rq14ECTR8lX5kqWRCdSNaK3rYCtv3AD9/PpWDBgu7pytl5vHjt1asXCxYsuOp2/vrrL+68uwuBbQbgX7xsFkQmkrWiNi2g4OHfmDfrW/z9/XN03d999x2dO3dm48aN3HjjjaxcuTJH1y/eI6ty9tmzZ2nZtgO26u0EVrs1CyITyVrn/vmLs8s+YdH8uVx33XU5um5PyNl5unht0qQJxYsXv6o2Tp48ye2t2+FfrzsFKoZlTWAiWSh2z2/E/zqVJYu+z3D0oZTmzJlD/fr1qVWrFi1atHDfa3Do0KH07t2bpk2bcsMNNzB69Gj3Mm+++SY33XQTt956K3/++ScA8+fP5/333+fjjz+mWbNmQPYcNZC8IStytrWW7vfezyFblMD6XbMoMpGsk3DmKFFzhvHFhPHUrl07U8vktZytC7auQkJCAh3uvJvoEtUJCm3jdDgiacQfP0DU9yOYM/NrqlTJ/KnRW2+9ldWrV2OM4bPPPmP48OG89957gOuU0tKlSzlz5gxVq1alb9++/Pbbb0ybNo3IyEgSEhKoXbs24eHhtGvXjscee4ygoCCee+657NpMkUx76ZUh/LzhDwrf/Xqeu/WQ5H5J8bFEzRnGwGee4O677870cnktZ6t4vQqP9XuCrf9EUSjiaadDEUkjMTaKM7PfZMQ7b7n3oDNr//79dOvWjUOHDhEXF0elSpXcr7Vv3558+fKRL18+SpUqxeHDh1m+fDmdOnVy98uKiIjI0m0RyQpTpkxl9LjPKNp9uC7QklzH2iSiF35Aq1vCefmlwZe1bF7L2Xm628DVeP+DMcyYu4igts9qQALJdWxiAtHzh3Nf10489uijl738k08+yRNPPMHmzZv55JNPUo3w4uljYkve9Ouvv/JIvyco1PFFfAM1IIHkPtG/TOH6AnH8b+Lnl31WIK/lbBWvV2DhwoW8NPQ1giIG45Mv0OlwRFKx1hK9bDxhFUoyetR7V9TGqVOnKFvWdfHhF198ccn5mzRpwnfffcfZs2c5c+YMc+bMuaL1imSHvXv30rbjHRRs8TgBpSpdegGRHBa9dSl+u39h4bzZVzSEcl7L2Xm6eO3RowcNGzbkzz//pFy5ckyYMOGSy/zxxx906X4PQe2ew7/otTkQpcjlid4wlyKn/+a7b77C1/fSZwViYmIoV66c+zFy5EiGDh1Kly5dCA8Pp0SJEpdso3bt2nTr1o3Q0FDatm1L3bp1s2JTRFK5kpwdFRVFizbt8QnpQMHK9XMgSpHLc+7AH8Qu/y+LF8ynVKlSl5xfORuMtTbTM+eFcbIzcvToUUJr1yUupBOBNW53OhyRNM7uXEfc0o/YuG4NFStWdDqcXCWjcbK9VV7P2UlJSbTtcAfrDscT1OIJXaAluU7CqX85Of15pv/vv7Rv397pcHKVjHJ2nj7yejni4uJo2/FOzpWvq8JVcqW4I3uI/mE0c2fNVOEqAjw7cBBrt+8jqPljKlwl10k6F8OZ2W8w5MXnVbheJhWvmWCtpddDD/P3aQi89T6nwxFJIzHmFGfmvMnY0e9zyy23OB2OiOMmfP45E76cTmD75zG+OTswh8ilnB+u+45WTXnu2QFOh+NxVLxmwtvD32X+slUEtn4GY/SWSe5iE+KJmvs2jz14P/ffn7M7V1k5XKdIVvn55595esBACkUMxrdAYafDEUkjevkkbromHxM+HZfjZwW8IW+rEruEWbNm8eY7Iwjq+CI+AQWcDkckFWst0T+OpWHwDbwz7M0cX39WDdcpklV27txJRKe7CWzdH/9ryjsdjkga0b8tosChjcyfPTPHh+sG78jbKl4zsGnTJu7t1ZtCHV7Ar3BJp8MRSSN67beUjP+XGdMm4+OT81/nrBiuUySrnDp1ihat2+FXpwsFKtVyOhyRNGL3/kbc6sksXjjfsdzpDXlbI2xdxD///EPLtu3J36QP+a6r6nQ4ImnEbF8FW75n8YZ1BAbqfsOStyUkJNCxU2dOF69KoVq6+EVyn/gTBzkzfwTfzZhG1aqqK66Gjrym4+zZs7Rq15GkKk0JvLmJ0+GIpBF3eCcxP45lwbw57htTi+Rljz/1NL8fOElQ0z5OhyKSRmJsFFGz32D4W6/TsmVLp8PxeCpeL2Ctpef9vTiQGERgg+5OhyOSRkLUcc7MeYvPP/2EOnXy1G1LRdL14UcfMW3mPILaPqfhuiXXcQ3X/S49OkXweL9+TofjFVS8XmDIq6+xdO1mAls8qfsCSq6TFH+OqDnDGPBkP7p27eJ0OCKOW7RoES+8NISgiJfwyR/kdDgiaUT/NIGQ8sX5cPQop0PxGipeU5g2bTqjPhxHYIf/4ON/+WMLi2Qnay3RP4yhef0Qhr7ystPhAFc2XKdIVtm2bRudu91DULuB+Bcr43Q4ImlEb5hLoRPbmfXtDPz8csdlRt6Qt3PHO5kLrF27lj6P9aVwp1fxC/Lsq/DEO0Wvmko5/2imTJqYa84KTJ061ekQJI86duwYLVq3I1/DnuQv77n3qxTvdfbv9cSv/4Yla1dTpEgRp8Nx84a8reIV2L9/P206RFCweT8CSt/gdDgiaURv/QnfnT+zaMM68ufP73Q4Io6Ki4ujXUQnYsvWJiikldPhiKQRd3Qv0Ys+4Ps533HDDaorslqeL16jo6Np0aY9JrgtBW9q6HQ4ImmcO/gnscsn8MtPSyldurTT4Yg4ylrLQ488xl8nEgnqcL/T4YikkRhziqjZbzJm1Hs0btzY6XC8Up7u85qUlESX7j054l+awLp3OR2OSBoJp49wZu7bTP7iv4SEhDgdjojj3h3xHnMWLyewTX/dWUByHZsYT/S84fS5vwcPPtjL6XC8Vp4uXge98CKrtu4i6PZ+uaYPoch5SXFniZrzJoOff46IiAinwxFx3Jw5c3ht2DsERWi4bsl9rLVELxlHvWrXM2L4O06H49XybLeBiRO/4NOJX1Kk+3CMX86PLSySEWuTiF44ivbNbuGFQQOdDkfEcZs3b+ae+x6gUMfB+BUu5XQ4ImnErJvJNecO8vX0lY4M152X5MnidcWKFTzxzAAKd34D34K55wpAkfOil0/ixiI+/PezT3VWQPK8f//9N3m47ofIV7aa0+GIpBGz41eSNs9j8fq1BAXpfsPZLc8Vr7t376bjnXcR2OppAkpc73Q4ImlEbV5MvgPr+H7DOgICApwOR8RRsbGxtGrXkYRKtxJUvanT4YikEffv38Qs+ZAfFy2gfPnyToeTJ+Sp4vX06dM0b9UW31p3UeCGcKfDEUkjdt/vxK2cxIpVv3DNNdc4HY6Io6y13PvAg+w7l5+g23s4HY5IGolRJzgz+y3GfzyWevXqOR1OnpFnOmUkJiZy591dOV3kRgJrd3A6HJE04k/+Q9T8d5kxbQrVqunUqMirr7/B4tWRBLZ6CmPyzM+VeIik+HNEzR3G0/0eoUeP7k6Hk6fkmWzw5NMDiNxzhMCmD6sPoeQ6SeeiiZr9Bm++OpTWrVs7HY6I4776agYjPviIoA7/wcdfA3NI7mKtJXrxhzStE8zrrw51Opw8J08Urx+PG8eXX39HYLtBGN881VNCPIBNSiRq/gi6RrTh6aeecDocEcetX7+e3o88SqGO/8GvkLrPSO4TvXo61/mcZsr/vtABMQd4fSW3ZMkSBr4wmCJdh+GbX1cASu4T/dMEgq8rzMcfjnE6FBHHHThwgNbtO1KwWV/yXVvZ6XBE0oj+Yzk+O5ayaP1aChTQ/Yad4NXF6/bt27mrSzeC2j6Hf/GyTocjkkZ05HyCjv7BnPVr8PPz6q+jyCXFxMTQsm0HuLkVBas2cjockTTOHdrO2Z8+ZcWyHylTpozT4eRZXvtreeLECVq0bod//R7kr6BhNSX3ObtrI/Frv2LJmtUULVrU6XBEHJWUlESXHj057FuSoHqdnQ5HJI2E00c5M/dtvpz4OWFhYU6Hk6d5ZZ/X+Ph42kV0IqZ0CEGhbZwORySN+GP7iFo4ilnffs2NN97odDgijnvhxZdYuXmnhuuWXCkpLpaouW/xn+f6c+eddzodTp7ndcWrtZZHHuvHn0djCWzSy+lwRNJIPHuaM7Pf5P0R73Dbbbc5HY6I4yZN+h8fT5hIUPvnNVy35Drnh+tu26Q+L77wvNPhCF5YvI56/wO+XfAjQW2exfj4Oh2OSCo2MZ7oecPpdU9XHu7Tx+lwRBy3atUq+j75NIUiBuMbWNTpcETSiF7xJZUKWSZOGK+zArmEV/V5nT9/Pq+89iZFu7+DT76CTocjkoq1lqgfP6FOlbKMGjHc6XBEHLdnzx7aR9xJUKunCChZ0elwRNKI/n0JAft+ZcGGdeTLl8/pcCSZ1xSv1lp69e6D9fElbsmHxGXDOo7s2JgNrV6+kjfVvuJlj2zf4P6/r38+ilcKzoqQrkrKmHJKycphObYui8W3fk8S/tlO8Zh9fPPVKnx9dVZA5KkBzxEVHU3RTbOJ2TQ7y9uPiznDqQN/ZXm7lyurcvbVtpVVHMnZVWrl6Pp8b7oVW7Qc5375gp9+WU6JEiVydP2SMa8pXo0x/LDwe06cOJFt69i+fTuTJk26ZBxX8/qlPPTQQ5QvX/6Kl583bx7/+9//OHLkCInx5ziyfQP33XcfvXr1uqq4rsZ7773H2bNns6y9pUuXpju9ZcuWANxyyy00a9Ysy9Z3KR+N/ZhZy7+kYMw/LF63hkKFCuXYukVys9EjR7Drycezrf24uDjeeOONS86XnaeCW7Rowa233nrFy//555+MGDGCnTt3Aq7CsVy5cpf8LcpOU6ZMYfv27VnW3po1a4iNjU0z/fw1ARUrVszR36hffvmF4R9/QWL0CWZM+ZLq1avn2Lolk6y1mX6Eh4db8R4HDx60VatWtYB99NFHnQ4nW+zcudMWLlzYAvadd95xJIYBzz5n/fwD7KpVqxxZv7gA6+xl5DtveChne5dz587ZRx55xAIWsElJSU6HlOXOnDljIyIiLGBr1qzpSAyzZs2ygB056gNH1i8uGeVsr7tgSzKvTJkybNu2jTFjxvDJJ5/w1ltvOR1Slrvhhhs4deoUzzzzDM8//zxTp07N8Rge7tOH7+fPp0GDBjm+bhHxHgEBAXzyySds3rwZwCtvkh8UFMSsWbNYvHgxmzdvpl27djkewy233MKXkyfzzNNP5vi6JXOMq7jNnDp16th169ZlYzjilP79+/P+++9z7tw5AgICnA4nW9x22238/PPPJCUl6YrRPMgYs95aW8fpOHKScrb3+vXXX2nQoAHz58+nbdu2ToeTLWbMmEHXrl3566+/dD/sPCijnK0jrwLAqFGjAOjXr5/DkWSfxYsXA/Dxxx87HImIyNWpX78+gCNHJnNKly5dAHQ/bElDxau4tWrVigkTJjgdRrbx9/enePHiPP549l0gIiKSU5YsWQLAsWPHHI4k+4wYMYIDBw6QmJjodCiSi6h4FbevvvoKgOjoaIcjyT5z5851OgQRkSzRvHlz4P/PnHmjAQMGAK77uIucp+JV3IoUKQLArl27HI4k+4SGhgKQlJTkcCQiIlnjr7+cv5dtdjl/fcKOHTscjkRyExWvkkZUVJTTIWSbggVdI68lJCQ4HImISNbw5px93pkzZ5wOQXIRFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIqIiIiIx/BzOgCRnGCtZf78+VhrAZg3bx7+/v60a9cOHx/tw4mI5DYrV67k2LFj7v/PmTOHRo0acc011zgcmThNxavkCfv27eOOO+7Ev1BxAO556DHiTx9n+/Y/ueGGGxyOTkREUrLWcs99D7B/3z4Afv51I4sWRfDl5Mn0vOceh6MTp+mQk+QJ119/PYsX/0CASaJwaCsCSGThwgUqXEVEciFjDMuX/ciNVW6iaHBjfI1l5KgPuKdHD6dDk1xAxavkGU2bNmX9mlV0qHkta1ev5Pbbb3c6JBERuYjy5cuz/teV9GhWixlTJtH/macwxjgdluQC6jYgeUrlypWZ/L8vnA5DREQyISgoiLFjPnA6DMlldORVRERERDyGilcRERER8RgqXkVERETEY6h4lVztn3/+oXv37tx4442Eh4fTrl07tm/fnu68u3fvpkaNGgAsW7aMDh06XNE633//fWJiYi5rmatZn4iIt1DOlpyg4lVyLWstnTp1omnTpuzcuZP169czbNgwDh8+nK3rvZJEKCKS1ylnS05R8Sq51tKlS/H39+exxx5zTwsNDeXWW29l4MCB1KhRg5o1azJ9+vQM24mOjqZ3797Uq1ePWrVqMWvWLAASExN57rnnqFGjBiEhIYwZM4bRo0dz8OBBmjVrRrNmzQBYtGgRDRs2pHbt2nTp0oWoqCgAFixYQLVq1ahduzbffvttNr0LIiKeQTlbcopulSW51u+//054eHia6d9++y2RkZFs2rSJo0ePUrduXZo0aXLRdt58802aN2/O559/zsmTJ6lXrx4tWrRg0qRJ7N69m8jISPz8/Dh+/DjFixdn5MiRLF26lBIlSnD06FHeeOMNFi9eTGBgIO+88w4jR45k0KBBPPzww/z4449UrlyZbt26ZedbISKS6ylnS05R8SoeZ8WKFfTo0QNfX19Kly7Nbbfdxtq1awkJCUl3/kWLFjF79mxGjBgBQGxsLHv37mXx4sU89thj+Pm5vgbFixdPs+zq1avZunUrt9xyCwBxcXE0bNiQbdu2UalSJapUqQLAvffey6effpodmysi4tGUsyWrqXiVXCs4OJivv/76qtux1vLNN99QtWrVK1q2ZcuWTJ06NdX0yMjIq45LRMSbKGdLTlGfV8m1mjdvzrlz51LtHf/2228ULVqU6dOnk5iYyJEjR/j555+pV6/eRdtp3bo1Y8aMwVoLwMaNGwFo2bIln3zyCQkJCQAcP34cgEKFCnHmzBkAGjRowC+//MJff/0FuPpibd++nWrVqrF792527twJkCZRiojkNcrZklNUvEquZYxh5syZLF68mBtvvJHg4GD+85//cM899xASEkJoaCjNmzdn+PDhXHvttRdt5+WXXyY+Pp6QkBCCg4N5+eWXAejTpw/XX3+9u60pU6YA8Mgjj9CmTRuaNWtGyZIlmThxIj169CAkJMR9+il//vx8+umntG/fntq1a1OqVKkceU9ERHIr5WzJKeb8nk1m1KlTx65bty4bwxGnGWNYtWoVDRo0cDoUkSxljFlvra3jdBw5STnb+xljaN++PXPnznU6FJEslVHO1pFXEREREfEYKl5FRERExGOoeJVcLSgoKNe0vWzZMlauXJlN0YiIeD7lbMkJKl5FMkmJUETEcyhney8Vr+Jxdu7cSZs2bQgPD6dx48Zs27YNgDlz5lC/fn1q1apFixYt3ONpR0VF8eCDD1KzZk1CQkL45ptv3G0NHjyY0NBQGjRo4J4/vXZ2797NuHHjGDVqFGFhYSxfvpwjR45w9913U7duXerWrcsvv/yS82+GiEgup5wtWc5am+lHeHi4Fe8G2FWrVjkdhltgYGCaac2bN7fbt2+31lq7evVq26xZM2uttcePH7dJSUnWWmvHjx9vBwwYYK21dtCgQfbpp592L3/8+HFrrWtbZ8+eba21duDAgfb111/PsJ0hQ4bYd999191Ojx497PLly6211u7Zs8dWq1YtazZasgWwzl5GvvOGh3K29wNs+/btnQ7DTTlbskpGOVsjbIlHiYqKYuXKlXTp0sU97dy5cwDs37+fbt26cejQIeLi4qhUqRIAixcvZtq0ae75ixUrBkBAQAAdOnQAIDw8nB9++CHDdi60ePFitm7d6n5++vRpoqKisrXPl4iIJ1HOluyg4lU8SlJSEkWLFk13qL8nn3ySAQMGEBERwbJlyxg6dGiGbfn7+2OMAcDX19c9aktm20lKSmL16tXkz5//ajZJRMRrKWdLdlCfV/EohQsXplKlSsyYMQNwdXvZtGkTAKdOnaJs2bIAfPHFF+5lWrZsyUcffeR+fuLEiQzXcbF2Ug5BCNCqVSvGjBnjfq6xs0VEUlPOluyg4lVytZiYGMqVK+d+jBw5ksmTJzNhwgRCQ0MJDg5m1qxZAAwdOpQuXboQHh5OiRIl3G289NJLnDhxgho1ahAaGsrSpUszXOfF2unYsSMzZ850d/4fPXo069atIyQkhOrVqzNu3LjseRNERDyEcrbkBA0PK6loeFjxVhoeVryRhocVb6XhYUVERETEK6h4FRERERGPoeJVRERERDyGilfxSAsWLKBq1apUrlyZt99+O83r586do1u3blSuXJn69euze/du92vDhg2jcuXKVK1alYULF16yzZ49e1K1alVq1KhB7969iY+Pz9ZtExHxVrGxsdSrV8998daQIUMAaNq0KVWrViUsLIywsDA6d+6c7vLr16+nZs2aVK5cmaeeeorz1+0MHTqUsmXLupefP39+jm2T5DwVr+JxEhMTefzxx/n+++/ZunUrU6dOTXXjaYAJEyZQrFgx/vrrL/r378/zzz8PwNatW5k2bRpbtmxhwYIF9OvXj8TExAzb7NmzJ9u2bWPz5s2cPXuWzz77LMe3WUTEG+TLl48ff/yRTZs2ERkZyYIFC1i9ejUAkydPJjIyksjISL7++ut0l+/bty/jx49nx44d7NixgwULFrhf69+/v3v5du3a5cj2iDNUvIrHWbNmDZUrV+aGG24gICCA7t27u2+9ct6sWbN44IEHAOjcuTNLlizBWsusWbPo3r07+fLlo1KlSlSuXJk1a9Zk2Ga7du0wxmCMoV69euzfvz/Ht1lExBsYY9wjWsXHxxMfH+8eeOBSDh06xOnTp2nQoAHGGO6//36+++67bIxWcisVr+JxDhw4QPny5d3Py5Urx4EDBy46j5+fH0WKFOHYsWMXXTYzbcbHx/O///2PNm3aZMdmiYjkCYmJiYSFhVGqVClatmxJ/fr1AddZrvOn/QcOHJhmuQMHDlCuXDn38wvz9IcffkhISAi9e/e+5MAG4tlUvIpkUr9+/WjSpAmNGzd2OhQREY/l6+tLZGQk+/fvZ82aNfz+++9A6m4D77777mW12bdvX3bu3ElkZCRlypTh2WefzY7QJZfwczoAkctVtmxZ9u3b536+f/9+99CAF85Trlw5EhISOHXqFNdcc02Gy2bU5quvvsqRI0f45JNPsmuzRETylKJFi9KsWbNU/VZTSkxMJDw8HICIiAj69u2bqttWyjxdunRp9/SHH36YDh06ZGPk4jQdeRUSEhJo1a4jDZs0B6Bbz/todNvt/P333w5Hlr66deuyY8cOdu3aRVxcHNOmTSMiIiLVPBEREe4xrr/++muaN2+OMYaIiAimTZvGuXPn2LVrFzt27KBevXoZtvnZZ5+xcOFCpk6dio+PvjIi4rynnh1EgybNAJg3bx4NmjRjxjczHY7q0o4cOcLJkycBOHv2LD/88APVqlVLd97zR2gjIyN57bXXKFOmDIULF2b16tVYa5k0aRJ33HEH4OoPe97MmTOpUaNGtm+LOEdHXoWEhARsUhLrN2ygWPM+/LvlR3wOHsDf39/p0NLl5+fHhx9+SOvWrUlMTKR3794EBwfzyiuvUKdOHSIiInjooYe47777qFy5MsWLF2fatGkABAcH07VrV6pXr46fnx8fffQRvr6+AOm2CfDYY49RoUIFGjZsCMBdd93FK6+84szGi4gAxQoHsWblcgo37EriyUP8tnE9JYoXdTqsSzp06BAPPPAAiYmJJCUl0bVrVzp06MCIESPo2bMnBQoUAKBEiRIsXrw4zfJjx46lV69enD17lrZt29K2bVsABg0aRGRkJMYYKlasqLNkXs6cv0daZmicbO+VmJjIM88OZNzHY2l4S2PmzPyaIkWKOB2WSJbJaJxsb6Wc7d2mTJlKn0cfo3jxa1i8cP5Fj2CKeKKMcraKV0nl999/p1q1avj56aC8eBcVr+KN9uzZQ5EiRShatKjToYhkqYxytioUSUX9hEREPEeFChWcDkEkx+nqExERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVw/SrFkzFi5cmGra+++/T9++fdOdv2LFihw9ejTDNt96661Uzxs1agTA7t273QMWrFu3jqeeeupKw3br1asXX3/99VW3IyLiCZSzRbKHilcP0qNHD6ZNm5Zq2rRp0+jRo8cVt3lhIly5cmWaeerUqcPo0aOveB0iInmRcrZI9lDx6kE6d+7MvHnziIuLA1x72gcPHuTAgQPUrFmTGjVq8Pzzz6e77J133kl4eDjBwcF8+umnALzwwgucPXuWsLAwevbsCUBQUFCaZZctW0aHDh0AaNeuHWFhYYSFhVGkSBG++OILEhMTGThwIHXr1iUkJIRPPvkEAGstTzzxBFWrVqVFixb8+++/Wf6eiIjkVsrZItnEWpvpR3h4uBVntW/f3n733XfWWmuHDRtmH3zwQVu+fHn777//2vj4eNusWTM7c+ZMa621FSpUsEeOHLHWWnvs2DFrrbUxMTE2ODjYHj161FprbWBgYKr2zz/ftWuXDQ4OttZau3TpUtu+fftU861bt87WrFnTnjx50n7yySf29ddft9ZaGxsba8PDw+3ff/9tv/nmG9uiRQubkJBgDxw4YIsUKWJnzJiRDe+KyKUB6+xl5DtveChnO085W+TKZJSzdeTVw6Q8DTVt2jQqVKhA06ZNKVmyJH5+fvTs2ZOff/45zXKjR48mNDSUBg0asG/fPnbs2HHFMRw9epT77ruPKVOmUKRIERYtWsSkSZMICwujfv36HDt2jB07dvDzzz/To0cPfH19ue6662jevPkVr1NExBMpZ4tkPRWvHuaOO+5gyZIlbNiwgZiYGMLCwi65zLJly1i8eDGrVq1i06ZN1KpVi9jY2Ctaf2JiIt27d+eVV15xXxxgrWXMmDFERkYSGRnJrl27aNWq1RW1LyLiTZSzRbKeilcPExQURLNmzejduzc9evSgXr16/PTTTxw9epTExESmTp3KbbfdlmqZU6dOUaxYMQoWLMi2bdtYvXq1+zV/f3/i4+Mzvf4XXniBkJAQunfv7p7WunVrPv74Y3c727dvJzo6miZNmjB9+nQSExM5dOgQS5cuvcqtFxHxLMrZIlnPz+kA5PL16NGDTp06MW3aNMqUKcPbb79Ns2bNsNbSvn177rjjjlTzt2nThnHjxnHzzTdTtWpVGjRo4H7tkUceISQkhNq1azN58uRLrnvEiBEEBwe7jx689tpr9OnTh927d1O7dm2stZQsWZLvvvuOTp068eOPP1K9enWuv/56GjZsmKXvg4iIJ1DOFslaxtUnNnPq1Klj161bl43hiIhkD2PMemttHafjyEnK2SLiqTLK2eo2ICIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIrbggULqFq1KpUrV+btt992OhwREcmAcrbkVSpeBXDdC/Dxxx/n+++/Z+vWrUydOpWtW7c6HZaIiKRDOVvyMhWvAsCaNWuoXLkyN9xwAwEBAXTv3p1Zs2Y5HZaIiKRDOVvyMhWvAsCBAwcoX768+3m5cuU4cOCAgxGJiMjFKGdLXqbiVUREREQ8hopXAaBs2bLs27fP/Xz//v2ULVvWwYhERORilLMlL1PxKgDUrVuXHTt2sGvXLuLi4pg2bRoRERFOhyUiIulQzpa8zM/pACR38PPz48MPP6R169YkJibSu3dvgoODnQ5LRETSoZwteZmx1mZ6Zo2TLSKeKqNxsr2VcraIeKqMcra6DYiIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMY63N/MzGHAH2ZF84IiLZpoK1tqTTQeQk5WwR8WAXzdmXVbyKiIiIiDhJ3QZERERExGOoeBURERERj6HiVUREREQ8hopXEREREfEYKl5FRERExGOoeBURERERj6HiVUREREQ8hopXEREREfEYKl5FRERExGP8H6RWkPVv9fAIAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "Landfill.plot_sankey_report(margin=.5, ft=6)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Sensitivity to Leachate generation rate per area (default: 250 gal/acre-day)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "ename": "NameError",
     "evalue": "name 'ps' is not defined",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mNameError\u001b[0m                                 Traceback (most recent call last)",
      "\u001b[1;32m<ipython-input-7-bb67ac9736d4>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mLandfill\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mps\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mLandfill\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m      2\u001b[0m \u001b[0mIncomingWaste\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mps\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mIncomFlow\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      3\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      4\u001b[0m \u001b[0mLF_leach\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mlinspace\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m100\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m50000\u001b[0m\u001b[1;33m,\u001b[0m\u001b[1;36m20\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      5\u001b[0m \u001b[0mPFAS_stored\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;33m[\u001b[0m\u001b[1;33m]\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mNameError\u001b[0m: name 'ps' is not defined"
     ]
    }
   ],
   "source": [
    "Landfill = ps.Landfill()\n",
    "IncomingWaste = ps.IncomFlow()\n",
    "\n",
    "LF_leach = np.linspace(100,50000,20)\n",
    "PFAS_stored = []\n",
    "IncomingWaste.set_flow('MSW', 1000)\n",
    "IncomingWaste.calc()\n",
    "for i in LF_leach:\n",
    "    Landfill.InputData.Water_Blnc['leach_gpad']['amount'] = i\n",
    "    Landfill.calc(IncomingWaste.Inc_flow)\n",
    "    PFAS_stored.append(Landfill.LF_storage.PFAS['PFOA']/Landfill.Inc_flow.PFAS['PFOA'] * 100)\n",
    "plt.plot(LF_leach,PFAS_stored)\n",
    "plt.ylim(98,101)\n",
    "plt.xlim(100,50000)\n",
    "plt.xlabel('Leachate generation rate per area (gal/acre-day)')\n",
    "plt.ylabel('Percent of Incoming PFAS that \\n remains in the landfill (%)')"
   ]
  }
 ],
 "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
}
