{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Landfill"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import PFAS_SAT as ps\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from IPython.display import Image\n",
    "import pandas as pd\n",
    "pd.set_option('display.max_colwidth', 0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Model document\n",
    "\n",
    "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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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             NaN        \n",
       "1   kg/m3                NaN        \n",
       "2   fraction wet weight  NaN        \n",
       "3   fraction TS          NaN        \n",
       "4   gal/acre-day         NaN        \n",
       "5   1:TRUE/0:FALSE       NaN        \n",
       "6   fraction             NaN        \n",
       "7   fraction             NaN        \n",
       "8   Mg/acre              NaN        \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             NaN        \n",
       "20  log L/kg             NaN        \n",
       "21  log L/kg             NaN        \n",
       "22  log L/kg             NaN        \n",
       "23  log L/kg             NaN        \n",
       "24  log L/kg             NaN        \n",
       "25  log L/kg             NaN        \n",
       "26  log L/kg             NaN        \n",
       "27  log L/kg             NaN        \n",
       "28  log L/kg             NaN        \n",
       "29  log L/kg             NaN        \n",
       "30  log L/kg             NaN        \n",
       "31  log L/kg             NaN        \n",
       "32  log L/kg             NaN        \n",
       "33  log L/kg             NaN        \n",
       "34  log L/kg             NaN        \n",
       "35  log L/kg             NaN        \n",
       "36  log L/kg             NaN        \n",
       "37  log L/kg             NaN        \n",
       "38  log L/kg             NaN        "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Landfill = ps.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</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>746</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</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>840</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>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   \n",
       "2   Moisture flow  kg   746   \n",
       "3   VS flow        kg   134.62\n",
       "4   Carbon flow    kg   84.074\n",
       "5   PFOA           μg   1070  \n",
       "6   PFOS           μg   840   \n",
       "7   PFBA           μg   0     \n",
       "8   PFPeA          μg   0     \n",
       "9   PFHxA          μg   0     \n",
       "10  PFHpA          μg   0     \n",
       "11  PFNA           μg   0     \n",
       "12  PFDA           μg   0     \n",
       "13  PFBS           μg   0     \n",
       "14  PFHxS          μg   0     "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncomingWaste = ps.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": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deVxUZf//8dfAiLKoaS6hqKiQpOwKKRKCiruUK5KV5XKLmZWVlpVp3eZt6lfNpdu0XFOh3NBCRUpDTW9AGdQbUSpxQXNFBASB4fz+4OfcoaCozByWz/Px4PFozlznms+Z8O3xnHNdl0ZRFIQQQpiemdoFCCFEdSUBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEu3DNG7QoIFib29vpFKEEKJqOnz48FVFURrevf2hAtje3p74+Pjyq0oIIaoBjUZzpqTtcglCCCFUIgEshBAqkQAWQgiVSAALIYRKJICFEEIlEsBCCKESCWAhhFCJBLAQQqhEAlgIIVQiASweyNzcHHd3d8NPamrqI/WzdOlS1qxZA8CqVau4cOGC4b3Ro0eTlJRUHuUKUWloHmZJog4dOigyFLn6sbGxISsrq1z79Pf3Z+7cuXTo0KFc+xWiItJoNIcVRbnnl13OgMUjWbVqFW+88Ybhdb9+/di7dy8A3377LU8//TT+/v6MGTPG0G769OnMnTuXjRs3Eh8fz/Dhw3F3dycnJwd/f3/i4+P597//zeTJk4t9zoQJEwD47rvv8Pb2xt3dnbFjx6LX6013wEIYgQSweKCcnBzD5YcBAwbct+2FCxf45z//yaFDh9i9ezfJycn3tBk8eDAdOnRg3bp16HQ6LC0ti723efNmw+vw8HCCg4M5ceIE4eHhHDhwAJ1Oh7m5OevWrSu/gxRCBQ81G5qoniwtLdHpdGVqGxsbS5cuXahfvz4AQ4YM4dSpU2X+rIYNG9KqVSsOHTqEo6MjJ0+epHPnzixZsoTDhw/j5eUFFP2l0KhRo4c/GCEqEAlg8Ui0Wi2FhYWG17m5uQA8zD2F0gQHB/P999/j5OTEgAED0Gg0KIrCiBEj+Ne//vXY/QtRUcglCPFI7O3t0el0FBYWcu7cOWJjYwHw9vbm119/JT09nYKCAjZt2lTi/rVr1yYzM7PE9wYOHMjWrVvZsGEDwcHBAHTr1o2NGzdy+fJlAK5fv86ZMyVOsSpEpSFnwOKRdO7cmZYtW+Li4oKzszOenp4ANG3alA8//JBnn32WJk2a0LZtW+rWrXvP/q+++iqhoaFYWlpy8ODBYu/Vq1ePtm3bkpSUhLe3NwBt27ZlxowZ9OjRg8LCQmrUqMGSJUto0aKF8Q9WCCORx9BEucvKysLGxoaCggIGDBjAyJEjH3jzToiqTB5DEyYzffp03N3dcXZ2pmXLlrzwwgtqlyREhSSXIES5mzt3rtolCFEpyBmwKJGNjY3R+ra3t+fq1atlbq/T6YiMjDRaPUKoRQJYVHgSwKKqkgAWZXblyhUGDRqEl5cXXl5eHDhwACgafOHj44OHhwc+Pj6cPHkSAL1ez3vvvYeLiwuurq4sWrTI0NeiRYvw9PTExcXFMFqupH7y8vL45JNPCA8Px93dnfDwcLKzsxk5ciReXl54eHgQERFh+i9DiPKgKEqZf9q3b6+I6sHa2vqebSEhIcq+ffsURVGUM2fOKE5OToqiKEpGRoaSn5+vKIqi7N69Wxk4cKCiKIry1VdfKQMHDjS8d+3aNUVRFKVFixbKwoULFUVRlCVLliijRo26bz8rV65Uxo8fb6hjypQpytq1axVFUZT09HTF0dFRycrKKsejF6J8AfFKCZkqN+FEmUVHRxebMvLmzZtkZmaSkZHBiBEjSElJQaPRkJ+fb2gfGhqKVlv0a3ZneDIUDbYAaN++vWHuh9L6uVtUVBTbtm0z3OzLzc3l7NmzPPPMM+V/0EIYkQSwKLPCwkIOHjxYbPIcgAkTJhAQEMCWLVtITU3F398fKPrXlUajKbGvmjVrAkVzDRcUFAAwderUEvu5m6IobNq0iTZt2pTPgQmhErkGLMqsR48eLF682PD6zgQ9GRkZNG3aFCiaPvLv7ZcuXWoI2OvXr9+3/9L6uXvYcs+ePVm0aJFh3omEhIRHPyghVCQBLEp069Yt7OzsDD/z5s1j4cKFxMfH4+rqStu2bVm6dCkAkydPZsqUKXTu3LnYHL2jR4+mefPmuLq64ubmxvr16+/7maX1ExAQQFJSkuEm3NSpU8nPz8fV1RVnZ2emTp1qnC9BCCOTochCCGFkMhRZCCEqGAlgIYRQiQSwEEKoRAJYmMTOnTtp06YNDg4OzJo1S+1yhKgQJICF0en1esaPH8+OHTtISkpiw4YNxQZ0CFFdSQALo4uNjcXBwYFWrVphYWHBsGHDZP4GIZAAFiaQlpZGs2bNDK/t7OxIS0tTsSIhKgYJYGF0JT1rXtoQZSGqEwlgYXR2dnacO3fO8Pr8+fM0adJExYqEqBgkgIXReXl5kZKSwunTp8nLyyMsLIygoCC1yxJCdTIbmjA6rVbL4sWL8ff3x9zcnFGjRtGuXTu1yxJCdXIGLIyuoKCAiB9/4sKlyzRq0oxx48apXZIQFYIEsDC6VatXs+zfX1FwO5f/HIjh7XcnqV2SEBWCBLAwupBhwzh27BgA3377Lf83W0bCCQESwMIErK2tcXZ2BsDJyYmGDRuqXJEQFYMEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJICFSdxZmPOpp55SuRIhKg4JYGESc+bMAcDe3l7dQoSoQCSAhUm8//77tGzZEjMz+ZUT4g750yCMrmvXrgDs2bNH5UrUk5WVxcmTJ9UuQ1QwEsDCaH7++Wc0Gg179uxh48aNtGjRQu2SVPPll1/i0b4DJ06cULsUUYFo7twcKYsOHToo8fHxRiyn+tDr9QweNpwbGRmP3dfe3TvLoaL78w/sVea2d9fTqLEtbV3dyrukSiX1zz+4kFnAkzULOXokngYNGqhdkjAhjUZzWFGUDndv16pRjID8/HxOnjzJmWwzrN37PlZflq2vcfviqXKq7F5WDs/y39peZW5fx9uK3HPHqe3RB7NatVGA/xqtukrCzYvGzV249Z/v6d3/BQ78+gsWFhZqVyVUJmfAKkpLS8OtvRdmPq9h9bSP2uUIE1CUQrJ+nEVPTwfWrVmFRqNRuyRhAqWdAcs1YBU1bdqUXT9t59Yv/ybv0h9qlyNMQKMxw7rnRCL3/Mas2XPULkeoTAJYZe3bt2fFsq/J3D6TgqzrapcjTMDMwhKboI/4/Iu5REREqF2OUJEEcAUwdOgQ3n1zPNnb/0Vh/m21yxEmoK3TkNr9PuClV0eSmJiodjlCJRLAFcS0qR/TraMb2bsXoiiFapcjTKBmkzbU8htNYO++/PXXX2qXI1QgAVxBaDQavlu9kmYWOWQfDFO7HGEi1s/4UejQhR59+pOTk6N2OcLEJIArkFq1arHrp22Y/xFDdtKvapcjTMS6UwhpBTYMf+VVHuapJFH5SQBXMI0bNyZ6ZyS5+77ldlqy2uUIE9BoNFgHTmBP3DGmffpZqe3Mzc1xd3fH2dmZ/v37c+PGDaPWFR8fz5tvvvnAdj4+8gjlo5IAroBcXFxYv2YVmT99QcHNy2qXI0zArEZNrPtNYf7ipYSFhZfYxtLSEp1Ox/Hjx6lfvz5Lliwxak0dOnRg4cKFD2z322+/GbWOqkwCuILq378/n0yZTNb2mRTmybXB6kBrU5/a/T9kdOg44uLi7tu2U6dOpKWlGV7PmTMHLy8vXF1dmTZtGgCpqak4OTkxevRonJ2dGT58ONHR0XTu3BlHR0diY2MBiI2NxcfHBw8PD3x8fAyTBu3du5d+/foBMH36dEaOHIm/vz+tWrUqFsw2NjaG9v7+/gwePBgnJyeGDx9uuKQSGRmJk5MTvr6+vPnmm4Z+qzsJ4Aps0nvv0q+rL9k756MU6tUuR5iAReNWWHV9nV79gjh//nyJbfR6PT///DNBQUEAREVFkZKSQmxsLDqdjsOHDxMTEwPA77//zltvvcXRo0dJTk5m/fr17N+/n7lz5zJz5kwAnJyciImJISEhgc8++4wPP/ywxM9NTk5m165dxMbG8umnn5Kfn39Pm4SEBBYsWEBSUhJ//vknBw4cIDc3l7Fjx7Jjxw7279/PlStXyuOrqhIkgCswjUbDiuVf0/oJM7L3r1W7HGEiVk93QtOuN9179SU7O9uwPScnB3d3d5588kmuX79OYGAgUBTAUVFReHh44OnpSXJyMikpKQC0bNkSFxcXzMzMaNeuHd26dUOj0eDi4kJqaioAGRkZDBkyBGdnZyZOnMh//1vyzB19+/alZs2aNGjQgEaNGnHp0qV72nh7e2NnZ4eZmRnu7u6kpqaSnJxMq1ataNmyJQAhISHl+XVVahLAFZyFhQU7tm+lVtphso/tVrscYSLWXgO5atGYwcNepLCw6LnwO9eAz5w5Q15enuEasKIoTJkyBZ1Oh06n4/fff2fUqFEA1KxZ09CnmZmZ4bWZmRkFBQUATJ06lYCAAI4fP8727dvJzc0tsaa/92Vubm7Y/0Ft5MmO0kkAVwJPPvkk0bsiuf3bd+SePaZ2OcIENBoN1l1f51BSKpM/KH5JoG7duixcuJC5c+eSn59Pz549WbFiBVlZWUDRJE+XL5f95m1GRgZNmzYFYNWqVeV2DHc4OTnx559/Gs64w8NLvslYHUkAVxJOTk5sDF9P1o655KdfVLscYQIabQ1s+r7PslXfsWrV6mLveXh44ObmRlhYGD169ODFF1+kU6dOuLi4MHjwYDIzM8v8OZMnT2bKlCl07twZvb787zVYWlry1Vdf0atXL3x9fWncuDF169Yt98+pjGQ6ykpm0eIlfDxzLnWDv8CsprXa5QgTyLt6lpsbP2bnjxH4+vqqXc4jycrKwsbGBkVRGD9+PI6OjkycOFHtskxGpqOsIia8MZ7g53uTFTlHnoyoJiwaNMe6x1v0f2Egp0+fVrucR7J8+XLc3d1p164dGRkZjB07Vu2SKgQ5A66ECgoK6NqjN0lZVth0/Yfa5QgTyTq8nTpnfkUX/x/q1KmjdjniIcgZcBWi1WrZvmUjta8nk5Xwk9rlCBOx9uzHzbqteX7gEKNcqxWmJ2vClcDGxsZwR9mUPvnkE/z8/OjevfsD29atW5foXZG09+5IzhNNsGzpUez98/8eiZmFJWg0mFvXo0HfdzC3qfe/7WZFf/fWD3ydWnbPAHAzbivpv66m2YTvDNeXC/NzubZzEfmXUwEFs5o2NBr6aVEfwqQ0Gg3W/mNIjPiMCW+9w1eLvyyXfn18fGQ4sUokgCuQzz4rfSKWkrRu3ZptWzbRO2gA2iGfU+PJZsXebxwyE3OruqT/upqMQ99Tv/vYYtvvln3iV2raOnLr1EFsXIr+EsiM34a51RM0HFX0zGn+tfNozOTXRi0acy3WfSbzXdhkXJyfYVxo6GP3KeGrHrkEcR/3G9seFxeHj48Pbm5ueHt7k5mZSW5uLq+99houLi54eHiwZ88eoOjZyhdeeIH+/fvTsmVLFi9ezLx58/Dw8KBjx45cv160FNGrr77Kxo0bAbC3t2fatGl4enri4uJCcnLRzGhXrlwhMDAQT09Pxo4dy8svv8yM6VPJ3PY5+pybJR5HrWbO5KdfuO+x5qdfRMnL5YnnXi42FaY+Ox1t7ScNr2s8aYdGW+MRv1FRHsxr2VA76CMmffARv/zyy2P3V5a5HIRxSAA/QElj2/Py8ggODubLL78kMTGR6OhoLC0tDSOTjh07xoYNGxgxYoRhVNHx48dZv349sbGxfPTRR1hZWZGQkECnTp1Ys2ZNiZ/doEEDjhw5wrhx45g7dy4An376KV27duXIkSMMGDCAs2fP8vJLL/Ha8KFk/fQFiv7e8fk5f8Ri0dDe8PrShg+5sHICF9e8Y9iWfeJXrJ7xo2azduRfT0OfXTTVoY1LIBn/2cTFte+SHrOW/Otpd3cvVFCjflNser/HgMFDOXXqVLn1W9LvuzAeCeAHKGls+8mTJ7G1tcXLywuAOnXqoNVq2b9/Py+//DJQNHCiRYsWhj8cAQEB1K5dm4YNG1K3bl369+8PUGxM/t0GDhwIFC3ceafN/v37GTZsGAC9evWiXr16AMybMxsvRzuyf/nacNZyJ2gLb9+ibschhn4bh8ykyWuLsH1lnmHbrRMxWD/jh0ZjhtXTnbh1cj9QNDlM07HfUNd7EIW5mVxc8w75V8899vcqHl+tFq7UeDaEbj17G/4V9bhK+n0XxiMX8x6gtLHtGo3mnrb3++daWcbkl7bP38fdl/YZ5ubmbP4hDA+vjqQf3gqUfq33bnmXT5OffoFL4VOLNhQWoK37FLU9i6YMNLOwxKqND1ZtfNBoNOT8GUeNBs3u06MwFRu3XmSlp9H3+YHE/LKbGjUe7/JQWeZ7EOVHzoAfgZOTExcuXDDM2ZqZmUlBQQF+fn6sW7cOgFOnTnH27FnatGlTrp/t6+vL999/DxTNgpWenm54z8bGhp93RVKYuB2loOyrK2ef+JUnOr+I3bgVRT/j11CQdY2CjMvknk9Cn1v0RIiizyfv6jnM6zQq12MSj8fa71VOXb3NP0Jfl2u2lYwE8COwsLAgPDycCRMm4ObmRmBgILm5ubz++uvo9XpcXFwIDg5m1apVxc4oysO0adOIiorC09OTHTt2YGtrS+3atQ3vN2/enMjtESg5meRfPVumPrNP7MPq6U7Ftlk5diT7RAwFNy5yaf0HXPh2PBdXvkXNpxywatO5XI9JPB6NmTnWvd5h885fmDd/gdrliIcgI+Eqmdu3b2Nubo5Wq+XgwYOMGzcOnU53T7t169YT+vZ7PDFsNubW9VSoVJhaQcYlboS9z/frVtOnTx+1yxF/U9pIOLkGXMmcPXuWoUOHUlhYiIWFBcuXLy+x3fDhL3I86QRfh82hzpCZJq5SqEFbtzE2fd9n2PCXOZpwGHt7e7VLEg8gAVzJODo6kpCQUKa2hYWFUMLNQlGFKUVDlEu6SSwqHrkGXIGNHDmSRo0a4ezs/ND7rl37HV99swrr3pOMUJmoiPJv/MWl9R/wQ9h6WrRoUW79Tp8+3fAcenJyMu7u7nh4ePDHH38YlqRPTU01/J7+fTFPcX8SwBXYq6++ys6dOx96v4MHDzJuwlvUDvoQc+snjFCZqGgKb2eTtW0GFhYW9OzZ02ifs3XrVp5//nkSEhJo3bq1DGN+TBLAFZifnx/169d/qH3OnDlD36AXsA6cUGz0m6i6lEI9WZFzGRrUq8TngLdv386zzz6Lh4cH3bt3Nyymeb+l5j///HPatGlD9+7dDcvUR0ZGsmDBAr755hsCAgKA/w1jFo9GArgKyczMpFvPPpi7P49lay+1yxEmkv3rCto1qcO/Fy8q8X1fX18OHTpEQkICw4YNY/bs2Yb3Slpq/vDhw4SFhZGQkMDmzZsNz7v36dOH0NBQJk6caJjnRDweuQlXRej1egYMDuaGjT3WnkFqlyNMJFu3A5urSWw/HItWW/If5/PnzxMcHMzFixfJy8szLA8P/1tqvmbNmoal5vft28eAAQOwsrICIChIfp+MRc6Aq4i3353EkT8vYR3wD7kDXk3kpOrIjwsnelckTzxR+rX+CRMm8MYbb3Ds2DG+/vrrYsvOlzb0WH6HTEMCuApYtnw5q8M2Yt1nEhpz+UdNdZB/7TzZO+cRsXkjDg4O923792XnV69efd+2UHTvYcuWLeTk5JCZmcn27dvLpWZxLwngCiwkJIROnTpx8uRJ7Ozs+Pbbb+9ps3fvXt6Z9AG1gz7G3LJ2Cb2Iqkafc5PMbTOYP/cLunTpUuy9W7duYWdnZ/iZN28e06dPZ8iQITz33HM0aNDggf17enoSHByMu7s7gwYN4rnnnjPWoVR7MhS5Evv999/p8Gwnaga+jaW9u9rlCBNQ9PlkbvmUl/p3ZeH8/1O7HFFGsihnFXPjxg269ehNDe9hEr7VhKIoZO9ZhmdrW+bPnf3gHUSFJxcMK6GCggL6vTCI7EbO2Lj1UrscYSK3jkRQL/ssW/YcxNzcXO1yRDmQAK6EQl9/g6S/sqgd9JbapQgTufVHHHrdNqLjY4tNPyoqN7kEUcks+HIRP/wYhU3vd9GYyVlQdZB3JZVbuxcRuT2iTHM87Ny5kzZt2uDg4MCsWbPueT89PZ0BAwbg6uqKt7c3x48fN7z35Zdf4uzsTLt27Viw4H9zCycmJtKpUydcXFzo378/N28WLQCbl5dnWIjWzc2NvXv3Pv4BVyMSwJXIzp07+Xj6Z9gEfYRZTWu1yxEmoM++Qea2mSxdvJCOHTs+uL1ez/jx49mxYwdJSUls2LCBpKSkYm1mzpyJu7s7R48eZc2aNbz1VtG/pI4fP87y5cuJjY0lMTGRH3/8kZSUFABGjx7NrFmzOHbsGAMGDGDOnDkAhulQjx07xu7du3n33XeLZuETZSIBXEmcOHGCoSHDsek7iRpPPKV2OcIElII8sn6axfgxr/HSS8PLtE9sbCwODg60atUKCwsLhg0bRkRERLE2SUlJdOvWDShaXis1NZVLly5x4sQJOnbsiJWVFVqtli5durBlyxYATp48iZ+fHwCBgYFs2rTpnr4aNWrEE088gTwpVXYSwJXA1atX6d6zDzU7j6CWXTu1yxEmoCgKWdFL6OzqyMwZn5V5v7S0NJo1+9+CqXZ2dqSlpRVr4+bmxubNm4GiwD5z5gznz5/H2dmZmJgYrl27xq1bt4iMjOTcuaIVsJ2dndm2bRsAP/zwg2G7m5sbERERFBQUcPr0aQ4fPmx4TzyYBHAFl5eXR+/+L3C7mTfWzt3ULkeYSHbsRp5SrvP9+u8wMyv7H9OSnuu/e1jxBx98QHp6Ou7u7ixatAgPDw+0Wi3PPPMM77//PoGBgfTq1Qs3NzfD/BIrVqxgyZIltG/fnszMTCwsLICiOavt7Ozo0KEDb7/9Nj4+PqXOSSHuJd9UBaYoCq+OGsOfN8Gm30tqlyNMJPvkATgRxe4j8YYJccrKzs6u2Bno+fPnadKkSbE2derUYeXKlUDR71jLli0NE/SMGjWKUaNGAfDhhx9iZ2cHFF2qiIqKAopW/P7pp58A0Gq1zJ8/39C3j48Pjo6OD1VzdSZnwBXYv76YTeTeg1j3fBuNRv5XVQe3//qdnD1LiYr88Z7gLAsvLy9SUlI4ffo0eXl5hIWF3TOb2Y0bN8jLywPgm2++wc/Pjzp16gBw+fJloGjtwc2bNxMSElJse2FhITNmzCA0NBQoGvqcnZ0NwO7du9FqtbRt2/YRjrx6kjPgCioiIoKZc+bxRPBszCws1S5HmEBB5jWytv+L1d8ux9PT85H60Gq1LF68mJ49e6LX6xk5ciTt2rVj6dKlAISGhnLixAleeeUVzM3Nadu2bbE5RgYNGsS1a9eoUaMGS5YsoV69ohW1N2zYwJIlSwAYOHAgr732GlAUzD179sTMzIymTZuydu3ax/kKqh2ZC6ICSkxMxNe/KzZBH1PT9mm1yxEmUJify80fPuKdMS/xydSP1S5HlDOZC6KS+Ouvvwjs3ZdafmMkfKsJRSkkO2ohPTq3Z+rHH6ldjjAhCeAKJCcnhx59+lPoGID1MzIFYHWR/dsGmte6zdpVK2Qi9GpGrgFXEIqiMPyVV0nT22DTMVjtcoSJZCftQXt6P7uOxBdbnUJUD3IGXEF8Mv1T9sQdx7r7BDkLqiZup50gd99KondG0qhRI7XLESqQAK4AwsLCWbDka2z6TcGshpwFVQcFGZfJ/Gk2G9auxtnZWe1yhEokgFUWFxfH6NDXqR30IeY29dQuR5hA4e1bZG3/nGkfvk+/fv3ULkeoSAJYRefPn6dXvyCsuo3HolErtcsRJqAU6snaNY/nA7vw3rvvqF2OUJnchFNJbm4u3Xv1Jde6CdaKnlunfnvkvjKPRHL74slyrK44K8eOWD3dqcztc88eI/fsMWp79sPcqo7R6oKi1YEzDoY/0r5K/m3Df2u0/7v0Y1bLhvrdxz52bXerZe/GrUPf06a+Bd8sWyrX+oUEsFosLCxwcXHB7kYG5Bx9rL5+PqMrp6pKlv3fPXRsVvb5h38+/CMA13ctBqDRU7a4uLc3Sm0/x/xYLv0oBf8LY33WbdxuP97/k7v9+XsK5/evo5G1OT8diaNGjRrl2r+opBRFKfNP+/btFSHKau/evQqgAEp4eLja5ZTqwoULiqurqwIovXv3NspnzJgxQ7GqXUdJTk42Sv+iYgPilRIyVa4BC6Pp0qULiqLQs2dPgoODOX36tNollcjW1pbExEQ2btzIjh07mDlzZrl/xsSJE0k8cpg2bdqUe9+i8pK5IIRJaDQamjdvzpkzZ9Qu5b66dOlCTExMifPqCvGoZC4Ioao5c+Zw9uzZCr9e2NatWwGIjo5WuRJRHUgAC5N49913AUhNTVW3kAe4M/3inekbhTAmCWBhEnceufrrr79UrqRsLl68qHYJohqQABZCCJVIAAshhEokgIUQQiUSwEIIoRIJYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEH9z9GjRWnDnzp3jjz/+ULkaUdVJAAvx/507dw43NzfDfz/XxV9WxhBGJQEsKrW//vqLYcOG0bp1a9q2bUufPn04depUiW1TU1NxdnYGYO/evfTr16/Y+82aNWP9+g1Y2tSloW1TIn/cXuLS8QsWLODWrVsPVWdJnyeEBLCotBRFYcCAAfj7+/PHH3+QlJTEzJkzuXTp0iP3GRIyjP8eTSDpqA53d/cS2zxKAAtREglgUWnt2bOHGjVqEBoaatjm7u6Or68vkyZNwtnZGRcXF8LDw+/bT3Z2NiNHjsTLywsPDw+OHj1KgwYN0Ov1vPfee7i4uODq6sqiRYtYuHAhFy5cICAggICAAACioqLo1KkTnp6eDBkyhKysLAB27tyJk5MTvr6+bN682fVW4dkAAAsKSURBVHhfhKi0tGoXIMSjOn78OO3bt79n++bNm9HpdCQmJnL16lW8vLzw8/MrtZ/PP/+crl27smLFCm7cuIG3tzfdu3dnzZo1nD59moSEBLRaLdevX6d+/frMmzePPXv20KBBA65evcqMGTOIjo7G2tqaL774gnnz5jF58mTGjBnDL7/8goODA8HBwcb8KkQlJQEsqpz9+/cTEhKCubk5jRs3pkuXLsTFxeHq6lpi+6ioKLZt28bcuXMByM3N5ezZs0RHRxMaGopWW/THpH79+vfse+jQIZKSkujcuTMAeXl5dOrUieTkZFq2bImjoyMAL730EsuWLTPG4YpKTAJYVFrt2rVj48aN92x/2CcXFEVh06ZNtGnT5p7tJd2Eu7tNYGAgGzZsKLZdp9M9cF8h5BqwqLS6du3K7du3Wb58uWFbXFwc9erVIzw8HL1ez5UrV4iJicHb27vUfnr27MmiRYsMwZ2QkABAjx49WLp0KQUFBQBcv34dgNq1a5OZmQlAx44dOXDgAL///jsAt27d4tSpUzg5OXH69GnDs8R3B7QQIAEsKjGNRsOWLVvYvXs3rVu3pl27dkyfPp0XX3wRV1dX3Nzc6Nq1K7Nnz+app54qtZ+pU6eSn5+Pq6srzs7OTJ06FYDRo0fTvHlzQ1/r168H4B//+Ae9e/cmICCAhg0bsmrVKkJCQnB1daVjx44kJydTq1Ytli1bRt++ffH19aVFixYm+U5E5aJ5mH+udejQQYmPjzdiOaIq02g0HDhwAB8fH7VLEcKkNBrNYUVROty9Xc6AhRBCJRLAQgihEglgUanZ2NgYrW97e3uuXr1a5vY6nY7IyEij1SOqHglgIcqJBLB4WBLAosq5cuUKgwYNwsvLCy8vLw4cOABAbGwsPj4+eHh44OPjw8mTJwFKHHJ8x6JFi/D09MTFxYXk5ORS+8nLy+OTTz4hPDwcd3d3wsPD7xniHBERYfovQ1RsiqKU+ad9+/aKEI8KUA4cOFCufVpbW9+zLSQkRNm3b5+iKIpy5swZxcnJSVEURcnIyFDy8/MVRVGU3bt3KwMHDlQURVG++uorZeDAgYb3rl27piiKorRo0UJZuHChoiiKsmTJEmXUqFH37WflypXK+PHjDXVMmTJFWbt2raIoipKenq44OjoqWVlZ5Xj0orIA4pUSMlVGwokqJzo6mqSkJMPrmzdvkpmZSUZGBiNGjCAlJQWNRkN+fr6hfWlDjgcOHAhA+/btDRPqlNbP3Uob4vzMM8+U/0GLSkkCWFQ5hYWFHDx4EEtLy2LbJ0yYQEBAAFu2bCE1NRV/f3/g/kOOa9asCYC5ublhRNzUqVNL7OduSilDnIW4Q64BiyqnR48eLF682PBap9MBRWeuTZs2BWDVqlXF2pc05Lg0pfXz9yHKUPoQZyHukAAWldqtW7ews7Mz/MybN4+FCxcSHx+Pq6srbdu2ZenSpQBMnjyZKVOm0LlzZ/R6vaGP0oYcl6a0fgICAkhKSjLchCttiLMQd8hQZGEyMhRZVFcyFFkIISoYCWAhhFCJBLAQQqhEAlhUSTt37qRNmzY4ODgwa9ase95XFIU333wTBwcHXF1dOXLkyAP3vX79OoGBgTg6OhIYGEh6ejoA165dIyAgABsbG9544w3jH5yoMiSARZWj1+sZP348O3bsICkpiQ0bNhQbmAGwY8cOUlJSSElJYdmyZYwbN+6B+86aNYtu3bqRkpJCt27dDOFcq1Yt/vnPfxoGXAhRVhLAwugURTE8rqXX6yksLDTq58XGxuLg4ECrVq2wsLBg2LBh98zDEBERwSuvvIJGo6Fjx47cuHGDixcv3nffiIgIRowYAcCIESPYunUrANbW1vj6+lKrVi2jHpeoeiSAhdGtWLnSMMzXz8+P10aNMernpaWl0axZM8NrOzs70tLSytTmfvteunQJW1tbAGxtbbl8+bIxD0NUAxLAwugGDxrEcwHdsW5oh519Kz6a8r5RP6+kZ9vvHmpcWpuy7CtEeZG5IITR1a1bl1+idrBhwwb69u1bbLIbY7Czs+PcuXOG1+fPn6dJkyZlapOXl1fqvo0bN+bixYvY2tpy8eJFGjVqZNTjEFWfnAELk9Bqtbz88stGD18ALy8vUlJSOH36NHl5eYSFhREUFFSsTVBQEGvWrEFRFA4dOkTdunWxtbW9775BQUGsXr0agNWrV/P8888b/VhE1SZnwKLK0Wq1LF68mJ49e6LX6xk5ciTt2rUzzAkRGhpKnz59iIyMxMHBASsrK1auXHnffQE++OADhg4dyrfffkvz5s354YcfDJ9pb2/PzZs3ycvLY+vWrURFRdG2bVvTH7yoVGQuCCGEMDKZC0IIISoYCWABgL+/P7t27Sq2bcGCBbz++usltk9NTcXZ2fm+faamphab2jE+Pp4333wTKJpH986osaVLl7JmzZrHKR94+FWMhVCbBLAAICQkhLCwsGLbwsLCCAkJeeQ+7w7gDh06sHDhwnvahYaG8sorrzzy5whRWUkACwAGDx7Mjz/+yO3bt4Gi8Lxw4QK+vr5MmjQJZ2dnXFxcCA8Pv2ff1NRUnnvuOTw9PfH09OS3334Dim5a7du3D3d3d+bPn8/evXvp16/fPftPnz6duXPncuHCBdzd3Q0/5ubmnDlzptRVjq9du0aPHj3w8PBg7NixJT7DK0RFJgEsAHjyySfx9vZm586dQNHZb3BwMJs3b0an05GYmEh0dDSTJk3i4sWLxfZt1KgRu3fv5siRI4SHhxsuM8yaNYvnnnsOnU7HxIkTH1hDkyZN0Ol06HQ6xowZw6BBg2jRogVvvfUWEydOJC4ujk2bNjF69GgAPv30U3x9fUlISCAoKIizZ8+W87cihHHJY2jC4M5liOeff56wsDBWrFjB2rVrCQkJwdzcnMaNG9OlSxfi4uJwdXU17Jefn88bb7yBTqfD3NycU6dOPVYdBw4c4JtvvmHfvn1A6ascx8TEGFYq7tu3L/Xq1XuszxXC1CSAhcELL7zAO++8w5EjR8jJycHT07NMN8fmz59P48aNSUxMpLCw8LEmpbl48SKjRo1i27Zt2NjYAKWvcgwyTFhUbnIJQhjY2Njg7+/PyJEjDTff/Pz8CA8PR6/Xc+XKFWJiYvD29i62X0ZGBra2tpiZmbF27VrDzGd3rxL8IPn5+QwdOpQvvviCp59+2rC9tFWO/fz8WLduHVA0veSd+XmFqCwkgEUxISEhJCYmMmzYMAAGDBhgWC24a9euzJ49m6eeeqrYPq+//jqrV6+mY8eOnDp1CmtrawBcXV3RarW4ubkxf/78B372b7/9RlxcHNOmTTPciLtw4UKpqxxPmzaNmJgYPD09iYqKonnz5uX8bQhhXDISTgghjExGwgkhRAUjASyEECqRABYm8aBFMoWojiSAhdGVZZFMIaojCWBhdGVZJFOI6kgCWBhdWRbJFKI6kgAWRicLXQpRMglgYXRlWSRTiOpIAlgYXVkWyRSiOpLJeITR3W+hSyGqMxmKLIQQRiZDkYUQooKRABZCCJVIAAshhEokgIUQQiUSwEIIoRIJYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqeagVMTQazRXgjPHKEUKIKqmFoigN7974UAEshBCi/MglCCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQyf8DR/1rCtR6dCgAAAAASUVORK5CYII=\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": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd3gU1dvG8e9JIYQkNBGkqIAiSA1VBcSgNGlSREBU+AEKIthFxYavBQVEBbGj2CiKIr0XpUoNikhVepEaUghp5/1jlxhISIEkk93cn+vKBdmdOfPMbvbZZ86cmWOstYiIiIiIeAIfpwMQEREREcksFa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8WriIiIiHgMFa8iIiIi4jFUvGYzY0xBY4w1xpRzOpbMMMbcYIw55XQcIiJOUd4W8Sz5ong1xkSl+EkyxpxJ8XuPDNZtZYzZmY2xrDbG3Jdd7V0ua+12a23RnGjbGDPJGHPW/TqfMMbMNcZUcj/3ljEm/oL35tEU63Y0xqwzxkQbY44ZY742xpROYxut3F86j174XAaxBRljfjTG7HGvf3MGy19pjJnhjucfY8zdKZ7rc8F+xLjbrJbG63DuZ01W4nW3U98YE+5uf40xpnqK5x4wxmw3xkQYY44YY8YZY4Kyug2RvEJ5++I8NW+7DxJGG2MOuNf92xjzdhZiq2eMWenOc/uMMc9lYp2CxphdF/49GGOuM8Ysc+fTP40xTVI8198Yk3DBfkYZY67IbKzudtL73rjVGLPY/RofNcZMNMaUzEr7+Vm+KF6ttcHnfoC9QLsUj33ndHxe7jX3634NEAl8luK5r1K+N9ba0QDuL6YvgeHAFUBNwBf41RhT+IL2ewIn3P9mhQV+Ae4FTmZi+U/dy5UE+gBfnEvo1tpxF/yNPQn8Za39M8X6r12wrw3S2oj7SzJVIW2MCQSmueMoBvwATDXG+LkX+QW4xVpbBLgeKAy8kon9EsmTlLcdlVN5+xXgRqAOEAI0A37PQlzfA3OA4u51nzLGtMhgnReAA2k8PgVY5m7rdeBnY0zKA4KlF+xnsLX2+IWNuAv6ixXRF/3eAIoCHwDXAuWBJPfykgn5onjNiDEm0Bgz1hhzyBiz3xgzwhjj7z7KmgpUTHnkZYxpZIz5zX30d9AY826KIiIr221ljNlpjBniPvI6kLJHwbh6B0e7jzAjjDG/nNuOMaazMWaLMeaUMWZhig8ExpjDxpgn3UeTUcaYj4wxpY0xC4wxp91H0oXdy1YxxiSkWHe1MeYV97+njTGzjTHFUjzf1xiz1x3vYPe2Gme0r9baKGASUD295dz7NxJ42Vr7vbU21lp7EHjAvcjAFMsWAe4C+gO1TIqeyEzEE2OtHW2tXYEraaQXUzGgHfCStTbaWrsYmAdcrPenJ/BVZmPJpOZArLX2Q2vtWeAdXMm/MYC1dk+KxGpw7dP1aTV04Xvufiy5Z8kY4+f+uzvu7rF49MLlRZymvO2Rebs+8KO19oh1+TuzByLGGF/gauA7a22itXYbsBqols46NwAdceXLlI/XBG7AVaTHWmsnAruADpmJJZPxpvu9Ya2dYa39yVobaa2NBsYCjdJp77z3zLiK5s9T/H5J77GnUvHq8iquo8QaQF0gDBjsLgY6An9fcOQVj+vDWBy4FdcfaN9L3Pa1uIqNMu42PzbGBLufGw1UwfWBLw68CFh3kTYeGIDriO4XYPoFibgjcBtQFeiGq9fuSaAUEAw8nE5M9+L6gJXGdXT4GIAxJhQYBdwDlHP/lMjMTrqTbndgYwaLVgeuwtWzmMxamwj8hKuIO+ce4CiuI+gl/Jcoz21zmzGmU2biy0AVINJauyfFY5tII2m6k2V94Nts2G5K1dzbBMBamwRsThmDMeYOY0wEEAG0xvX3cykG4vrbqQ40AO5Of3ERRyhvn88T8vZq4FnjOi2fVv4cZ4wZldYG3G2NAXq6D1Kq4XrfF6UT14fAM8DZCx6vBmy31p5J8ViaOf0yZPp7w60J8OdFnkvX5bzHnkrFq0sP4BVr7TFr7RFcpxDuv9jC1to11tq17qO/XcDnuBLOpYgBhllr4621U3Gdzr7eGOOPqxgbZK097N7WMvcHuDsw1Vq71FobB7yJ6w+1Xop233Pvz15gJbDCWvuH+8M6DaidTkyfWWt3uY8GpwCh7sfvwXXUvNrd+/ciGf8NvWBcFxZsw3UKKeWXxf3uHohzP8Xd+2GBI2m0dYjzP5A9gYnWWgtMAO5zH50DYK2tbK39KYP4MiMYV0GYUgSuns8L9QQWWmsvPE31wgX7+kl2x2CtXeQeNnANrkS2N4vbOOceYJS19pD7S3/4JbYjkpOUt8/nCXn7VeA9XHlyg7vHvPu5Ba21fay1T6YT18+4Xt8zuA7eP7DWpjnswN1upLV2ThpPZyan33bBfma1sMz094Yxpi7wrPvnUlzKe+zRvHrnMsMYY3AdMaY8OtoDlE1nnarGmDnGdWHMaeBlLv0o56i7F+2cGFx/9KUBP+DvNNYpkzJed2I8cEHMKZPImTR+D+biDqcRz7nt7kux3dOk/nBe6A1rbVFrbWlrbccLjkK/cT937ucEcAxXj0apNNoq7X4eY8x1uE6xnDvl9COusaDN01jvckXhGkOaUmFcY8GSuf+W7iPtIQNvXLCv/dzrBKRMkLi+yBakeOzxrMQAYK3dh6sn+lJ7f897ny/4v4jjlLfTlOfztrvYf99aewuufD0K+Nqdz9NljCkFzAaeBwJwjRPtbIzpncayhXEdzDxxkeYyk09/uWA/U57lWpAiZz8OvJIiZ0/JwjYwxtwIzAT6WWt/u0i8GbmU99ij5fvi1d1rdxjXaaBzruG/Ad42jdU+AzYA11lrCwP/h+uDm50OAQlAxTSeO0iKeN29jWVJe1B6dseUfCsZd4Ioks3b2IwrYXdJ+aB7Hzvy3ymicxdoLTDGHAa24/rSOG/oQDbZChQ2xlyT4rFapD7Fczuu1+PnzDZsrT2bMkEC64DmKR57z73on+5tAmCM8cF1qu5ivQF+QHpfCD7uXqJzUl6ocN77jGucmUieobyd5ZjySt5OZl3XHYzCdUq/Sia2UQlXT+okd4/2HlzDFFqnsWxVXHlrtfv7YQJQ3j0OtDSuvHmDMaZginXSyulpstY2T5Gz3wNeTZGzzw2zyvB7w120LwCGWGu/z8SmU8Z70ZydQ+9xnpLvi1e3ibiOnK4wrltVvMB/vVZHgJIpxjOBq9s/wlob5R5382B2B2StjQe+Bt43xpQyxvgaYxq7k8FkoKMxpom7AHkOOI6r8MlJ3+M60q1vjCmAK/mne7FTVllrE3CdOnnNGNPF3TNZBldvph/wgbvX5X5gCK5TY+d+7gU6mNR3JEiTu+1zyaDABYksZUwncR0Z/58xppAxJgxoxX+9vuf0BL6/YBxVdlkABLrHigXg6lGIBpZD8q2yyrn/XwHXe5PeWDAD9DbG+Bhj2gKVgRD3+LvvgSeMMVcZ18UvT+fA/ohcLuXtzMkTeRvAGPOUcd0iqqBxjVt9CNewhE0Xbfw/fwFBxpi7jUtZXOPx01p3Pa6DmXPfDY/g6pkMBY64hxrsAF5yx3oPrgtcp13yi3CBjL43jDHXAouBt621X2ay2QeMMQWMMdfjulg5xP23lOPvcV6j4tXlZWALriOicGAF/43z2wRMB/aY/8b3PAH0NcZE4bpCcHIOxfUorisgN+JKcq8Bxv3B6wN8guuCpTuAu9wJJMdYazfiGvw+FVdvwSFcpyYuHAx/udv5Ctf+PYfrNiN/uJ9qbK09hWucWkngI/e4ssPW2sO4xnkdwDX+B+O6Ur5zOpvag+tU3BW4Lp44Y4y5yr3uq8aYqSmWfRDXKcZjuC666GOt3XHuSfeXZCcufpeBl8z59wvcn8mXAwB3QXzuzgqncF3M0SHFe14TWGOMicZ1+5dNuC4MuZhE9zpHcSX214E3gAq4vmhW4vpMrMWVgLP1PRbJBsrbmZCH8jbubY7GdXDxL/A/XHlsP4AxZrwx5j3S4B5/3wXXsIFTuArUtcAI97rNjDEphyek/G44CSS6fz9X1HXBdZHUKWAo0NFdcJ4TZlLf57VGFl+W9L43+uMqsIelaP9YBu2dxfW9NQ3X31Vz4J7ceo/zEuM6+yKSdcZ1K5ATQBlr7SGn45HMMcZUATZbazN1myBjTEfgLWtt5ZyNTERymvK2Z3IPf7jbWrs8E8t6/XusnlfJEmNMe+O6v2IwrsH2v3nrhyO/MsaEGGNauE95XoPrytWpGa0nInmT8rb3y2/vsYpXyaouuC6U2I/rYoN0p2kUj+QDvIXrtNNaXBe5vO5oRCJyOZS3vV++eo81bEBEREREPIZ6XkVERETEY2SpeDXGzM2pQEREcpLyl4iId8jU1cbnFC5cuGW9evU0zkBEPNFppwPIbSVKlLDly5d3OgwRkSxbv379MWvtlWk9l6XitVKlSqxbl9P3UxYRyX7GmB0ZL+Vdypcvr5wtIh7JGLPnYs9pzKuIiIiIeAwVryIiIiLiMVS8ioiIiIjHUPEqIiIiIh5DxauIiIiIeAwVryIiIiLiMVS8ioiIiIjHUPEqIiIiIh5DxauX8PX1JTQ0NPln9+7dl9TOxx9/zNdffw3A+PHjOXjwYPJzffv2ZcuWLdkRrohIvqacLXLpjLWZn+21Xr16VrO15E3BwcFERUVla5thYWGMHDmSevXqZWu7Ik4wxqy31uarP2bl7LxLOVskfenlbPW8erHx48czcODA5N/btm3L0qVLARg3bhw33HADYWFhPPjgg8nLDR06lJEjRzJlyhTWrVtHjx49CA0N5cyZM4SFhbFu3To++ugjBg8efN52Bg0aBMC3335LgwYNCA0NpV+/fiQmJubeDouIeDDlbJHMUfHqJc6cOZN8+qljx47pLnvw4EFee+01Vq9ezYIFC9i6dWuqZe6++27q1avHd999R3h4OIGBgec999NPPyX/PnnyZLp27cpff/3F5MmTWbFiBeHh4fj6+vLdd99l306KiHgJ5WyRS+fndACSPQIDAwkPD8/UsmvWrOG2226jePHiAHTp0oXt27dneltXXnklFStWZPXq1VSqVIlt27bRqFEjxo4dy/r166lfvz7gSs4lS5bM+s6IiHg55WyRS6fi1Yv5+fmRlJSU/HtsbCwAWRnnfDFdu3bl+++/p0qVKnTs2BFjDNZaevbsybBhwy67fRGR/EY5WyRzNGzAi5UvX57w8HCSkpLYt28fa9asAaBBgwb88ssvnDx5koSEBH788cc01w8JCSEyMjLN5zp16sTPP//MxIkT6dq1KwB33HEHU6ZM4d9//wXgxIkT7NmzJwf2TETE+yhni2SOel69WKNGjahQoQI1atSgevXq1KlTB4CyZcsyZMgQbrrpJsqUKUPVqlUpUqRIqvV79epF//79CQwMZNWqVec9V6xYMapWrcqWLVto0KABAFWrVuX111+nRYsWJCUl4e/vz9ixY7n22mtzfmdFRDyccrZI5uhWWflUVFQUwcHBJCQk0LFjR3r37p3hRQMinky3yhJPppwt+Y1ulSWpDB06lNDQUKpXr06FChXo0KGD0yGJiMhFKGeL/EfDBvKpkSNHOh2CiIhkknK2yH/U8+rBgoODc6zt8uXLc+zYsUwvHx4ezuzZs3MsHhERT6ecLZI9VLxKtlAiFBHxHMrZ4slUvHqZo0eP0rlzZ+rXr0/9+vVZsWIF4LrJdcOGDalduzYNGzZk27ZtACQmJvL0009To0YNatasyZgxY5LbGjNmDHXq1KFGjRrJM7qk1U5cXBwvv/wykydPJjQ0lMmTJxMdHU3v3r2pX78+tWvXZtq0abn/YoiI5HHK2SKXwFqb6Z+6detayTuCgoJSPda9e3e7bNkya621e/bssVWqVLHWWhsREWHj4+OttdYuWLDAdurUyVpr7Ycffmg7deqU/Nzx48ettdZee+21dvTo0dZaa8eOHWv79OmTbjtffvmlfeSRR5LjeP755+0333xjrbX25MmTtlKlSjYqKiob914ka4B1Ngv5zht+lLPzFuVskcxLL2frgi0vs3DhQrZs2ZL8++nTp4mMjCQiIoKePXuyY8cOjDHEx8cnL9+/f3/8/Fx/CuemHwTXTa0B6tatmzwv9sXaudD8+fOZPn168kUGsbGx7N27lxtvvDH7d1pExEMpZ4tknYpXL5OUlMSqVasIDAw87/FBgwbRtGlTpk6dyu7duwkLCwNcPe/GmDTbCggIAMDX15eEhAQAXnrppTTbuZC1lh9//JHKlStnz46JiHgh5WyRrNOYVy/TokULPvjgg+Tfw8PDAdfRd9myZQEYP378ect//PHHyYnuxIkT6bZ/sXYunJawZcuWjBkzJnlO7o0bN176TomIeCnlbJGsU/HqwWJiYihXrlzyz6hRoxg9ejTr1q2jZs2aVK1alY8//hiAwYMH8/zzz9OoUSMSExOT2+jbty/XXHMNNWvWpFatWkyYMCHdbV6snaZNm7Jly5bkwf8vvfQS8fHx1KxZk+rVq/PSSy/lzIsgIuIhlLNFsoemhxWRfEHTw4qIeA5NDysiIiIiXkHFq4iIiIh4DBWvIiIiIuIxVLzmA4mJidSuXZu2bdsC0KtXLypUqEBoaCihoaE0bNgwzfX++ecfbrrpJipVqkTXrl2Ji4sDYMSIEcnrVq9eHV9f3wyveBURkYyVL1+eGjVqEBoaSr16ruF+l5uzIyIiaNeuHbVq1aJatWp8+eWXubY/IjlBxWs+8P7776e60fSIESMIDw8nPDyclStXprnes88+yxNPPMGOHTsoVqwY48aNA+CZZ55JXnfYsGHcdttt590oW0RELt2SJUsIDw8n5cV2l5Ozx44dS9WqVdm0aRNLly7lqaeeSi5sRTyRilcvt3//fmbNmkXfvn2ztJ61lsWLF3P33XcD0LNnT37++edUy02cOJHu3btnS6wiInJp0svZxhgiIyOx1hIVFUXx4sWTZ+gS8UQqXr3c448/zvDhw/HxOf+tfuaZZ5JPQfXo0SPVesePH6do0aLJCa5cuXIcOHDgvGViYmKYO3cunTt3zrkdEBHJR4wxtGjRgrp16/Lpp58mP345OXvgwIH89ddflClThho1avD++++n+k4Q8SQ69PJiM2fOpGTJktStW5elS5ee99yIESOSj9DTktb9fy+cknDGjBk0atRIQwZERLLJihUrKFOmDP/++y/NmzenSpUqwOXl7Hnz5hEaGsrixYvZtWsXzZs359Zbb6Vw4cI5sxMiOUyHXl5sxYoVTJ8+nfLly9OtWzcWL17Mfffdd9HlW7ZsSWhoKH379qVEiRKcOnUqeQrC/fv3U6ZMmfOWnzRpkoYMiIhko3N5tmTJknTs2JE1a9ZcdNnM5uwvv/ySTp06YYzh+uuvp0KFCmzdujXnd0Ykh6h49WLDhg1j//797N69m0mTJnH77bfz7bffXnT5efPmER4ezueff44xhqZNmzJlyhQAvvrqK+66667kZSMiIvjll1/Oe0xERC5ddHQ0kZGRyf+fP38+1atXv+jymc3Z11xzDYsWLQLgyJEjbNu2jYoVK+bw3ojkHBWv+VTK8VOhoaFpXnn69ttvM2rUKK6//nqOHz9Onz59kp+bOnUqLVq0ICgoKDfDFhHxWkeOHKFx48bUqlWLBg0a0KZNG1q1agVcXs5+6aWXWLlyJTVq1OCOO+7g7bffpkSJErm6byLZyaQ1TuZiNE+2iHiq9ObJ9lbK2SLiqdLL2ep5FRERERGPoeJVRERERDyGilcRERER8RgqXiVZfHw848eP5/jx406HIiIimTBr1iz+/PNPp8MQyVUqXgWAkydPEtasJQOfe4WaderpHoAiInlYYmIijz7xFF17PsjNjZswffp0p0MSyTUqXoXo6GhC6zbgt3UbSCx0JcfOQJ36Ddi8ebPToYmISBq6dL+PD0a/R2Lh0iQUuZq7OnRg/NffOB2WSK7Q9LBCoUKF6N3rAWJiYhj+9lu0a9+eG2+sSqlSpZwOTURE0tC+TSuuq3ANI4cPB+DpwYOp6p5KVsTb6T6vch5jDKtWreLmm292OhSRbKX7vIo3MsbQpk0bZs6c6XQoItlK93kVEREREa+g4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKV0mlQIECTocgIiKZpJwt+Y2KV0mWmJgIQPny5Z0NREREMk05W/IbFa+SbMiQIQAUL17c4UhERCQjf//9NwB9+vRxOBKR3KXiVZINHz6cunXrOh2GiIhkQlhYGADVqlVzNhCRXKbiVQCYNWvWef+KiEjedfLkSfbt28eIESOcDkUk16l4Ff7880/atm1L7dq1KVWqlNPheJ2YmBj27t3rdBgi4iViY2OTh3c99dRTDkfjfRITE9m5c6fTYUg6VLzmc59//jnVq1fH19eXDRs2OB2OV3rp5VeoVac++/btczoUEfFwGzduJDAwEICIiAiMMQ5H5H1mzZpFlSo3smjRIqdDkYvwczqA7BQfH4+1Nkfbj46Ovqw2MpNo0lsmMDAQX1/fS97+6dOn2bZtG6+++mryB/Pmm29myZIlxMbGXnK7cnEJiQlEWz+atWrD+t9WEhwc7HRIInlCQkJC8l1OcoK1lsjIyBxrPzMCAgIICAi45PXj4uLYsWMHEydOZOTIkcmPnzx5koCAAM6ePZsdYUoK8fHxFChSgk5durJ29UpuuOEGp0OSC3hN8Wqt5ZryFTly+CA+Ppde3KUnMSE+R9rNKj//S7+nX0J8XKrH1q3fQEjhIpcTkmSgaOsnOb53I3d3vZfZM37Gx0cnPUS6dO/BtB9/wMc3Z76K8krO9vXzv+R109oHXz9/SlxZ8nJCkgxc0aAdSUXK0qxlazZtWEuxYsWcDklS8Jri1RjDwIGPMOLDLyjcZRg+BQo6HZLIeWylm1jz01CeHvwso0bqIguRxwcOYP6ChRTt8gb+V1ztdDgiqUSdPEDr9h35dfEC/P0v/SBEspdXdf8Mee5ZWofdQvS8d7E2yelwRM5jfP0JajOYz7+eyJdfjnc6HBHH3Xbbbbw38m0ip79B4pnTTocjkkrQrT3ZfuwsD/UfkKPDEiVrvKp4NcYwftxnVAixRC//1ulwRFLxLVSE4PYvMOiJp1i+fLnT4Yg47sG+fel17z1EzXobm5g3TvOLnGN8fAlq9SQ/zV3MqHffczoccfOq4hVcczzPnTmNAvt+I3qzrhSUvKdAiWsIavEY7Tp0Sp4hRyQ/e3fkcOpdX5boxZ+od0vyHJ+AQoS0f4FXXh/G7NmznQ5H8MLiFaBEiRIsnDubsyu+Inb/n06HI5JKYMW6+NXpzB0tWxMREeF0OCKO8vX15acfJlEsZh/R66c5HY5IKn5FShHcejDdetzPn3+qrnCaVxavAFWrVmXyhG+JmjWC+FOHnQ5HJJWgOm2JLHYDd3XqQkJCgtPhiDgqJCSEhXNnkbRpOjE71zgdjkgqBcvdSMCt/6NZq9YcPXrU6XDyNa8tXgHuvPNOXnvlJaJmvEHS2RinwxFJJSisL7/vO8Ggx550OhQRx1177bXMnjGNmAVjiPv3H6fDEUklqGpT4ss3pGWb9rrHroO8ungFePyxQdzdujlRc0Zik3LuZtgil8L4+BLc+hkm/DSNDz/6yOlwRBx388038/HYMUTOeJPE6JNOhyOSSlDDe9kT488D/+ujMdoO8fri1RjDxx9+wI2lgoj+9UunwxFJxadgMMHtX2TwkJdYuHCh0+GIOO6++3rwyIO9iZr5FjYh9cQqIk4yxoegFo8xb/laXn9zmNPh5EteX7wC+Pv7M/PnHyl05A+iN811OhyRVPyLlSH4zqfofE83tm3b5nQ4Io578/X/o3FoZaIWjlXvluQ5PgUKEtxuCMNHjeann35yOpx8J18UrwDFihVj0fw5xK+ZxJk9m5wORySVgtfUpMDNPWjWsjUnTpxwOhwRR/n4+DD5u2+4yp4g+rcfnA5HJBW/kBIEt32OB3r3ZePGjU6Hk6/km+IVoFKlSvz84w9Ez3mH+BMHnA5HJJWgmi04U7o2bdp3JD5eN2yX/K1QoUIsmDMTs3UB0Vs1qYfkPQGlKxEY1o8Wrdty6NAhp8PJN/JV8QoQFhbGO8OHETn9dRLPRDodjkgqQbc+wLYTcfTt97BOl0q+V6ZMGebNnsmZpZ9w9tAOp8MRSSWoSmNs5Ttofmdbzpw543Q4+UK+K14B+j30EA907Uz07OHYRN1fU/IW4+NLcMsnmTZ/Ke+MetfpcEQcV6dOHb4a9xlRM98iIfKY0+GIpBJ00z0cohhde9xPUlKS0+F4vXxZvAK8P2oktSuUImrJp+rdkjzHJ6AQwe1eYOgbbzFr1iynwxFxXOfOnRn8xCCiZgwjKT7W6XBEzmOMIbjZIyzbuJUXXnrF6XC8Xr4tXn19fZk6ZTJFI/8hZsMMp8MRScWvSElC2gym230PsHnzZqfDEXHciy8MoUWjukTPex9r1bsleYvxK0Bw2+f44NNxfPfdBKfD8Wr5tngFKFy4MIvmzSYx/GfO7FrndDgiqQSUvZGC7ukI//33X6fDEXGUMYZvxn/BNYFxRK9UcSB5j29QMULaDaHfgIGsXr3a6XC8Vr4uXgHKly/PzGlTiV4wmriju50ORySVoKpNSSjfiFZt79J0hJLvBQQEMG/WdPz+WUH0n0ucDkcklQIlK1Co+UBat+/A3r17nQ7HK+X74hWgYcOGfDj6Pdd0hDERTocjkkpQo3vZc6YA9/fqneUx2r6+voSGhlK9enXatWvHqVOncihKl3Xr1vHoo49muFzDhg1zNA7xXiVLlmTh3NnELvuCswf+cjockVQKXX8TPjXb0qxVG6KiorK0rnJ2xlS8uj3wwP30/98D7ukIdX9NyVtc0xE+yvwV63n9jTeztG5gYCDh4f0Hg6gAACAASURBVOFs3ryZ4sWLM3bs2ByK0qVevXqMHj06w+VWrlyZo3GId6tevTqTvvuGyFnDSYjQkBrJe4LqduBEYDk6delGYmJiptdTzs6YitcU3h72Bg2rX0fUIk1HKHmPj39Bgts9z/B3x/Djjz9eUhu33HILBw78N0HHiBEjqF+/PjVr1uSVV1xXyO7evZsqVarQt29fqlevTo8ePVi4cCGNGjWiUqVKrFmzBoA1a9bQsGFDateuTcOGDZOntV26dClt27YFYOjQofTu3ZuwsDAqVqx4XoIMDg5OXj4sLIy7776bKlWq0KNHj+TP3+zZs6lSpQqNGzfm0UcfTW5XBKBNmzYMfeE5Iqe/TtLZGKfDETmPMYag2/uxdsd+nnxm8CW1oZydNhWvKfj4+PD9xG8pmXCU6LWXVhyI5CS/kBIEt3uenn0eZMOGDVlaNzExkUWLFtG+fXsA5s+fz44dO1izZg3h4eGsX7+eX3/9FYCdO3fy2GOP8fvvv7N161YmTJjA8uXLGTlyJG++6er5rVKlCr/++isbN27k//7v/xgyZEia2926dSvz5s1jzZo1vPrqq2nOHLZx40bee+89tmzZwt9//82KFSuIjY2lX79+zJkzh+XLl3P06NEs7a/kD089+QQdWoQRNfcdbFLme7dEcoPx9Se4zbN8+d33fD5uXJbWVc6+OBWvFwgKCmLh3FmwZR4x2/NOF7nIOQFXXU9gWD9atmnHwYMHM1z+zJkzhIaGcsUVV3DixAmaN28OuBLh/PnzqV27NnXq1GHr1q3s2OGawahChQrUqFEDHx8fqlWrxh133IExhho1arB7924AIiIi6NKlC9WrV+eJJ57gzz//THP7bdq0ISAggBIlSlCyZEmOHDmSapkGDRpQrlw5fHx8CA0NZffu3WzdupWKFStSoUIFALp3734pL5d4OWMMn3/6MZWvCCB62VdOhyOSim9gYULav8BjTw3ml19+yXB55eyMqXhNQ9myZZk7czoxiz8i7sgup8MRScU1HWEzmt/ZlpiY9E+Xnhs/tWfPHuLi4pLHT1lref755wkPDyc8PJydO3fSp08fwHVF9zk+Pj7Jv/v4+JCQ4JqV7qWXXqJp06Zs3ryZGTNmEBub9o3jU7bl6+ubvH5Gy2jojmSWv78/s6ZPJfBQONG/z3M6HJFU/K+4muCWT3BXp7vZtSv9ukI5O2MqXi+iXr16fPHpJ0TOeJOEqBNOhyOSStBNXThsimd6OsIiRYowevRoRo4cSXx8PC1btuSLL75IvhL2wIEDWbqXbEREBGXLlgVg/Pjxl7QP6alSpQp///13cq/B5MmTs30b4j2KFy/OovlziFs9gdi9vzsdjkgqgRVq41evC81atiYiIuM7GylnX5yK13Tcc08Xnhw0gOgZw0iK1/01JW85Nx3h8k3bGfLiy5lap3bt2tSqVYtJkybRokUL7r33Xm655RZq1KjB3XffTWRkZKa3P3jwYJ5//nkaNWqUpStpMyswMJAPP/yQVq1a0bhxY0qVKkWRIkWyfTviPW644QZ+/H4SUXPeIf5kxkNqRHJbcO02nC5embYdOqfZo3kh5ey0max089arV8+uW5e/ZqKy1tKl270s+esQwXc+hTGq9yVvSYw+xalJz/DRuyO4774eToeTraKioggODsZayyOPPEKlSpV44oknLqktY8x6a229bA4xT8uPORvgw48+4vn/e4vCXYfjWzDY6XBEzmOTEomc9hp3316fTz/60OlwslVu5WxVYhkwxvDtV19ydYEzRK+a5HQ4Iqn4BhUlpP0L9B/4qNdNR/jZZ58RGhpKtWrViIiIoF+/fk6HJB5gwMMP071je6Jnj8AmZty7JZKbjI8vwXc+zeSfZzPmg5y9h2tuy62crZ7XTDpy5Ag169QjqW53gqre5nQ4IqnE7FxDwi8fE75+Lddcc43T4eQ56nnNXxISEmjWsjWbIwsSfLsOeiTviT95iIjvn+OnyRNo0aKF0+HkOep5zQalSpViwZxZxC4bx9mD25wORySVQtc3wKdWe+5o2TpL46BEvJGfnx/TfvqBkJPbid4w0+lwRFLxL1aa4NbPcHfXe9m6davT4XiUXO95DQ4OzvI8v9nh5ZdfpkmTJjRr1uyy2pk+fTodOnbEt2hpjK8/vkHFKNHmSXyDi7H/o974FAgEH9cxQfHmAyhY7kbiju7hxMJPSIw8BliCqt1OkYbdMMYkt/vvj6+RGHOK0ve/c1nxSf5mrSVq4VjqlvRlzsxp+Pr65ti2GjZsmKemC8yIel4vnSfn7X/++Yc69W/i9Jl4fAsVAWMuO28nRp/k+JzRJJw+CkmJ+BYpSakur2bTXkt+E/37fAL+nMGmDWu54oorcmw73pSz803xmp2KFS9OYsGiFO02nIhV32PjYynerB/7P+pN6Z7vuhKkW1L8WQ598QjFWwwgsEIdkuJjOTp1GIWur09IHde0aUmxURz8YiDGvyAluwzFv+hVTu2aeAGbGE/k1Fe5v/3tvD9KB0PnqHi9dJ6et5cvX06TJrdR8t63KFiuKid/+eqy8vbxuR/gX+JqCte7C4C4f/+hQMkKTu2eeIGoX7+kIkdYvnQRBQoUcDqcPCFPDhtIb27ctWvX0rBhQ2rVqkWDBg2IjIwkNjaW//3vf9SoUYPatWuzZMkSwHWvsg4dOtCuXTsqVKjABx98wKhRo6hduzY333wzJ0647tHaq1cvpkyZAkD58uV55ZVXqFOnDjVq1Ejurj969CjNmzenTp069OvXj2uvvZZjx46lir1I4cK0uPUmoue+S0C5qunekiV6yy8ElL2RwAp1ANf89MWb9ydi9ZT/ltm2ksDrGxB0YxNi/vo1G15dyc+Mrz9BbQZf0nSEWZGZea7Fu3hq3m7cuDHFixcjau4oEmMiKHh19cvK24nRJ/ANKZG8vApXuVxBjR9g16kkej/YL8fypzflbEfHvKY1N25cXBxdu3bl/fffZ9OmTSxcuJDAwMDkGSb++OMPJk6cSM+ePZNnh9i8eTMTJkxgzZo1vPDCCxQqVIiNGzdyyy238PXXX6e57RIlSrBhwwYefvhhRo4cCcCrr77K7bffzoYNG+jYsSN79+69aOyj332H64r6cHrZtxS4snzy40cmDuHgl4M49PWTAMQf20OBq64/b13/YqWxcWdIOuuaGSnmr18IurEJQVVvI1rFq2SDc9MRPp7J6QgvV1qfZfFOnpq3g4ODuf+eTkTPGk7MjtWXlbdDarfh+JzRHJ74PBErJ5MQefyyXlMR4+NLUMvHmbloOcPdf9s5ydNztqPFa1pz427bto3SpUtTv359AAoXLoyfnx/Lly/n/vvvB1yzOFx77bVs374dgKZNmxISEsKVV15JkSJFaNeuHcB5c/peqFOnTgDUrVs3eZnly5fTrVs3AFq1akWxYsUuGnvLli05E3kKn8jD+BW+MvnxUt3fpMz/xlD6gVEpljapG3BLjD5J/MlDBJSrhn/xshgfX+KOph2zSFb4X3E1Qa2ezNR0hJcrrc+yeCdPztvLl/1KwtF/iNuzicI33Z38eFbzdmDFupTt9zkhtVoSf2I/h8Y/RmJMxjMmiaTHJ6AQwe2G8Nqw4UyfPj1Ht+XpOdvR4vVic+OmvJDpnPS6tDMzp+/F1kk5b29Wus2XLFnCH3/8wYb164hf9wOx+zanuZx/iWuIO7zjvMfiTx3GFAjEJ6AQ0X8tIyk2igMf92H/R71JiDhC9F/LMh2HSHoCy4fiX/8e7mhxZ6amI7xUmZkLW7yDp+ftgwcPUL70FcRuXnDR5TLK2wC+gSEEVQ2jRNunCChd6aLfASJZ4VekJCFtnqPHA734/fecm+bY03N2nrtVVpUqVTh48CBr164FIDIykoSEBJo0acJ3330HwPbt29m7dy+VK1fO1m03btyY77//HoD58+dz8uTJTMU7ZdIEomaPgKTU060FVQ3j7P4tnNkdDrguBDi58BMK39QZgOi/fqHUPa9S7uEvKPfwF1zV8z2Ne5VsFRTamqgSN2Z6OkKRrPKkvB0cHMyiebNJ2jybmB2/pblMRnn7zJ5NJMW7hj8knY0h/uSh887AiVyOgLJVKNikD83vbMORI0ecDidPynPFa4ECBZg8eTKDBg2iVq1aNG/enNjYWAYMGEBiYiI1atSga9eujB8//rwjh+zwyiuvMH/+fOrUqcOcOXMoXbo0ISEhGa7XokULhr32KklnTpEUF3Pecz7+AVzZ+SUiVk7iwGf9OPTFQApcVYmQOm1JiDhC4umjFChTJXl5/6JX4RNQSPeSlWwVdFsf/jx4mgGDHnU6FPFCnpa3y5Urx9yZ04lZNBaSUh/QpZe3AeIO7+TwV09w8IuBHP72aUJqtSCg9A3Zul+SvwVVDSOx4q20bNM+eZy4/EczbKVw9uxZfH198fPzY9WqVTz88MOEh4dnev2HHh7AlMXrCLnrRYxPzt1fU+RSJJ2NJmLys7w+5GkGDXzE6XBynW6V5Z0uJ29PmjSZBwc+TtFuI/ANvvhYWREnWJtE1OyRhN1YmimTJ6Y5NMeb5clbZeVFe/fupX79+tSqVYtHH32Uzz77LEvrfzhmNNXKFCZ62ficCVDkMvgEBBHc7gWef+mV5FsWiXi6y8nb3bp15bEBDxE1cxg2jWFfIk4yxoegFo+y6LdNvDHsLafDyVP8nA4gL6lUqRIbN2685PX9/PyoWb0aW5dsyMaoRLKPb6Ei+BQMYd/+/U6HIpItLjdv1wmtRdyYD13XLOiMmeQxxi8AnyKl2b7zb6dDyVPU85qNPv3sM775YSoF78h/p2Ql77NJiUTNfYfObZpz/333OR2OiOM2btzIA737Urj9Cxg/zWokeU/0bz9Q2pzi4w/edzqUPCVfF6+9e/emZMmSVK9e/bLbWrp0KU8+8xwh7V/ENzDji7xEclv0r+OpcmUgn3409ryxU+dmXclOQ4cOTb6J/NatWwkNDaV27drs2rWLhg0bArB79+7kz97SpUtp27Zttsch3iU7c/ahQ4do0botgWH9CChdKRuiE8le0VuXw9YFLJgzk0KFCiU/rpydz4vXXr16MXfu3MtuZ+fOnXTo3IWgVk/iX7xsNkQmkr2iNs2l0JHfmTXtJ/z9/XN12z///DN33XUXGzdu5LrrrmPlypW5un3xHtmVs8+cOUPzO9tiK99BUJXG2RCZSPY6e3gnZ5Z+wvzZMylTpkyubtsTcna+Ll6bNGlC8eLFL6uNU6dOcUfL1vg36EZg+dBsikwk+8Tu+Z343yayaP6cdGcfSmnGjBncdNNN1K5dm2bNmiXfa3Do0KH07t2bsLAwKlasyOjRo5PXeeONN6hcuTLNmjVj2zbXrd5mz57Ne++9x+eff07Tpk2BnOk1kPwhO3K2tZZu9z3AIVuUoJvuyabIRLJPQuQxomYM46txn1GnTp1MrZPfcrYu2LoMCQkJtO3QmegSVQmu1crpcERSiT9xgKg5I5kxdQqVKmX+1Gjjxo1ZvXo1xhg+//xzhg8fzjvvvAO4TiktWbKEyMhIKleuzMMPP8zvv//OpEmT2LhxIwkJCdSpU4e6devSunVr+vfvT3BwME8//XRO7aZIpr348iv8uuEvCnd+Ld/dekjyvqT4WKJmDOOZxwfSuXPnTK+X33K2itfL0H/AQLYcjiKk/WNOhyKSSmJsFJHT32Dk228mH0Fn1v79++natSuHDh0iLi6OChUqJD/Xpk0bAgICCAgIoGTJkhw5coRly5bRsWPH5HFZ7du3z9Z9EckOEyZMZPTHn1O023BdoCV5jrVJRM97nxaN6vLSiy9kad38lrPz9bCBy/He+2P4YeZ8gu98ShMSSJ5jExOInj2c++/pSP9+/bK8/qBBgxg4cCB//PEHn3zyyXkzvFxsTmz1Ykle9ttvv/HQgIGEtBuCb5AmJJC8J3rFBK4JjOOb8V9kOZ/mt5yt4vUSzJs3jxeH/h/B7V/AJyDI6XBEzmOtJXrpZ4ReeyWj333nktqIiIigbFnXxYdfffVVhss3adKEqVOncubMGSIjI5kxY8YlbVckJ+zdu5c7291FoWaPUKBkhYxXEMll0VuW4Ld7BfNmTb+kKZTzW87O18Vr9+7dueWWW9i2bRvlypVj3LhxGa7z119/0aXbvQS3fhr/olflQpQiWRO9YSZFTv/Nzz9+j69vxmcFYmJiKFeuXPLPqFGjGDp0KF26dOHWW2+lRIkSGbZRp04dunbtSmhoKJ07d+bWW2/Njl0ROc+l5OyoqCiatWqDT822FLr+plyIUiRrzh74i9hlX7Jw7mxKliyZ4fLK2WCstZleOD/Mk52eY8eOUatOfeJqdiSo+h1OhyOSypld64hbMpaN69ZQvnx5p8PJU9KbJ9tb5fecnZSUxJ1t72LdkXiCmw306NOk4p0SIv7l1ORnmfzNl7Rp08bpcPKU9HJ2vu55zYq4uDjubNeBs1fXV+EqeVLc0T1ELxjNzGlTVbiKAE89M5i12/cRfHt/Fa6S5ySdjSFy+uu8MuRZFa5ZpOI1E6y19OrzIH+fhqDG9zsdjkgqiTERRM54gw9Hv0ejRo2cDkfEceO++IJx304mqM2zGN/cnZhDJCPnpuu+q0UYTz/1pNPheBwVr5nw1vARzF66iqCWj2OMXjLJW2xCPFEz36L//x7ggQdy9+AqO6frFMkuv/76K489+Qwh7V/AN7Cw0+GIpBK97GtuuCKAcZ9+nOtnBbwhb6sSy8C0adN44+2RBLcbgk+BQKfDETmPtZboxR9yS7WKvD3sjVzffnZN1ymSXXbt2kX7jp0JavkE/ldc7XQ4IqlE/z6fwEMbmT19aq5P1w3ekbdVvKZj06ZN3NerNyFtn8Ov8JVOhyOSSvTan7gy/l9+mPQdPj65/3HOjuk6RbJLREQEzVq2xq9eFwIr1HY6HJFUYvf+Ttzq71g4b7ZjudMb8rZm2LqIw4cP0/zONhRs0peAMpWdDkcklZjtq+DPOSzcsI6gIN1vWPK3hIQE2nW8m9PFKxNSWxe/SN4Tf/IgkbNH8vMPk6hcWXXF5VDPaxrOnDlDi9btSKoURtCNTZwORySVuCO7iFn8IXNnzUi+MbVIfvbIo4+x+cApgsP6Oh2KSCqJsVFETX+d4W++RvPmzZ0Ox+OpeL2AtZYeD/TiQGIwQTd3czockVQSok4QOeNNvvj0E+rVy1e3LRVJ0wdjxzJp6iyC73xa03VLnuOarnsE3Tu255EBA5wOxyuoeL3AK6/+H0vW/kFQs0G6L6DkOUnxZ4maMYwnBw3gnnu6OB2OiOPmz5/Pcy++QnD7F/EpGOx0OCKpRP8yjppXF+eD0e86HYrXUPGawqRJk3n3g48Javs8Pv5Zn1tYJCdZa4leMIbbb6rJ0Jdfcjoc4NKm6xTJLlu3buXurvcS3PoZ/IuVdjockVSiN8wk5OR2pv30A35+eeMyI2/I23njlcwD1q5dS9/+D1O446v4BXv2VXjinaJXTaScfzQTvh6fZ84KTJw40ekQJJ86fvw4zVq2JuCWHhS82nPvVyne68zf64lf/yOL1q6mSJEiToeTzBvytopXYP/+/bRq255Ctw+gQKmKTocjkkr0ll/w3fUr8zeso2DBgk6HI+KouLg4WrfvSGzZOgTXbOF0OCKpxB3bS/T895kz42cqVlRdkd3yffEaHR1Ns1ZtMNXupNANtzgdjkgqZw9uI3bZOFb8soRSpUo5HY6Io6y19HmoPztPJhLc9gGnwxFJJTEmgqjpbzDm3Xe49dZbnQ7HK+XrMa9JSUl06daDo/6lCKrfyelwRFJJOH2UyJlv8d1XX1KzZk2nwxFx3IiR7zBj4TKCWj2hOwtInmMT44meNZy+D3Tnf//r5XQ4XitfF6+DnxvCqi3/EHzHgDwzhlDknKS4M0TNeIMXnn2a9u3bOx2OiONmzJjB/w17m+D2mq5b8h5rLdGLPqZBlWsYOfxtp8Pxavl22MD48V/x6fhvKdJtOMYv9+cWFkmPtUlEz3uXNk0b8dzgZ5wOR8Rxf/zxB/fe35OQdi/gV7ik0+GIpBKzbipXnD3IlMkrHZmuOz/Jl8Xr8uXLGfj4kxS++3V8C+WdKwBFzole9jXXFfHhy88/1VkByff+/fdf93TdfQgoW8XpcERSidnxG0l/zGLh+rUEB+t+wzkt3xWvu3fvpl2HTgS1eIwCJa5xOhyRVKL+WEjAgXXM2bCOAgUKOB2OiKNiY2Np0bodCRUaE1w1zOlwRFKJ+/dvYhZ9wOL5c7n66qudDidfyFfF6+nTp7m9xZ341u5EYMW6Tocjkkrsvs3Erfya5atWcMUVVzgdjoijrLXc1/N/7DtbkOA7ujsdjkgqiVEniZz+Jp999CENGjRwOpx8I98MykhMTKRD53s4XeQ6guq0dTockVTiTx0mavYIfpg0gSpVdGpU5NXXXmfh6nCCWjyKMfnm60o8RFL8WaJmDuOxAQ/RvXs3p8PJV/JNNhj02JOE7zlKUNiDGkMoeU7S2Wiipr/OG68OpWXLlk6HI+K477//gZHvjyW47fP4+GtiDslbrLVEL/yAsHrVeO3VoU6Hk+/ki+L1o48/5tspPxPUejDGN1+NlBAPYJMSiZo9knvat+KxRwc6HY6I49avX0/vh/oR0u55/EI0fEbynujVkynjc5oJ33ylDjEHeH0lt2jRIp557gWK3DMM34K6AlDynuhfxlGtTGE++mCM06GIOO7AgQO0bNOOQk0fJuCq650ORySV6L+W4bNjCfPXryUwUPcbdoJXF6/bt2+nU5euBN/5NP7Fyzodjkgq0eGzCT72FzPWr8HPz6s/jiIZiomJofmdbeHGFhSq3NDpcERSOXtoO2d++ZTlSxdTunRpp8PJt7z22/LkyZM0a9ka/5u6U/BaTaspec+ZfzYSv/Z7Fq1ZTdGiRZ0OR8RRSUlJdOnegyO+VxLc4G6nwxFJJeH0MSJnvsW3478gNDTU6XDyNa8c8xofH0/r9h2JKVWT4FqtnA5HJJX44/uImvcu036awnXXXed0OCKOe27Ii6z8Y5em65Y8KSkulqiZb/L800/QoUMHp8PJ97yueLXW8lD/AWw7FktQk15OhyOSSuKZ00ROf4P3Rr7Nbbfd5nQ4Io77+utv+GjceILbPKvpuiXPOTdd951NbmLIc886HY7ghcXru++9z09zFxPc6imMj6/T4YicxybGEz1rOL3uvYcH+/Z1OhwRx61atYqHBz1GSPsX8A3S8BnJe6KXf0uFEMv4cZ/prEAe4VVjXmfPns3L//cGRbu9jU9AIafDETmPtZaoxZ9Qr1JZ3h053OlwRBy3Z88e2rTvQHCLRylwZXmnwxFJJXrzIgrs+425G9YREBDgdDji5jXFq7WWXr37Yn18iVv0AXE5sI2jOzbmQKtZd+UNdS553aPbNyT/39c/gOIVqmVHSJclZUy55crrc2+wvcXie1MPEg5vp3jMPn78fhW+vjorIPLok08TFR1N0U3Tidk0Pdvbj4uJJOLAzmxvN6uyK2dfblvZxZGcXal2rm7P94bG2KLlOLviK35ZsYwSJUrk6vYlfV5TvBpjWDBvDidPnsyxbWzfvp2vv/46wzgu5/mM9OnTh6uvvvqS1581axbffPMNR48eJTH+LEe3b+D++++nV69elxXX5XjnnXc4c+ZMtrW3ZMmSNB9v3rw5AI0aNaJp06bZtr2MjP3wI6Yt+5ZCMYdZuG4NISEhubZtkbxs9KiR/DPokRxrPy4ujtdffz3D5XLyVHCzZs1o3LjxJa+/bds2Ro4cya5duwBX4ViuXLkMv4ty0oQJE9i+fXu2tbdmzRpiY2NTPX7umoDy5cvn6nfUihUrGP7RVyRGn+SHCd9StWrVXNu2ZJK1NtM/devWteI9Dh48aCtXrmwB269fP6fDyRG7du2yhQsXtoB9++23HYnhyaeetn7+BeyqVasc2b64AOtsFvKdN/woZ3uXs2fP2oceesgCFrBJSUlOh5TtIiMjbfv27S1ga9So4UgM06ZNs4Ad9e77jmxfXNLL2V53wZZkXunSpdm6dStjxozhk08+4c0333Q6pGxXsWJFIiIiePzxx3n22WeZOHFirsfwYN++zJk9m5tvvjnXty0i3qNAgQJ88skn/PHHHwBeeZP84OBgpk2bxsKFC/njjz9o3bp1rsfQqFEjvv3uOx5/bFCub1syx7iK28ypV6+eXbduXQ6GI0554okneO+99zh79iwFChRwOpwccdttt/Hrr7+SlJSkK0bzIWPMemttPafjyE3K2d7rt99+4+abb2b27NnceeedToeTI3744Qfuuecedu7cqfth50Pp5Wz1vAoA7777LgADBgxwOJKcs3DhQgA++ugjhyMREbk8N910E4AjPZO5pUuXLgC6H7akouJVkrVo0YJx48Y5HUaO8ff3p3jx4jzySM5dICIiklsWLVoEwPHjxx2OJOeMHDmSAwcOkJiY6HQokoeoeJVk33//PQDR0dEOR5JzZs6c6XQIIiLZ4vbbbwf+O3PmjZ588knAdR93kXNUvEqyIkWKAPDPP/84HEnOqVWrFgBJSUkORyIikj127nT+XrY55dz1CTt27HA4EslLVLxKKlFRUU6HkGMKFXLNvJaQkOBwJCIi2cObc/Y5kZGRTocgeYiKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRjqHgVEREREY+h4lVEREREPIaKVxERERHxGCpeRURERMRj+DkdgEhusNYye/ZsrLUAzJo1C39/f1q3bo2Pj47hRETympUrV3L8+PHk/8+YMYOGDRtyxRVXOByZOE3FkBFe5AAAEmFJREFUq+QL+/bt4667OuAfUhyAe/v0J/70CbZv30bFihUdjk5ERFKy1nLv/T3Zv28fAL/+tpH589vz7Xff0ePeex2OTpymLifJF6655hoWLlxAAZNE4VotKEAi8+bNVeEqIpIHGWNYtnQx11W6gaLVbsXXWEa9+z73du/udGiSB6h4lXwjLCyM9WtW0bbGVaxdvZI77rjD6ZBEROQirr76atb/tpLuTWvzw4SveeLxRzHGOB2W5AEaNiD5yvXXX89333zldBgiIpIJwcHBfDjmfafDkDxGPa8iIiIi4jFUvIqIiIiIx1DxKiIiIiIeQ8Wr5GmHDx+mW7duXHfddVStWpXWrVuzffv2NJfdvXs31atXB2Dp0qW0bdv2krb53nvvERMTk6V1Lmd7IiLeQjlbcoOKV8mzrLV07NiRsLAwdu3axZYtW3jzzTc5cuRIjm73UhKhiEh+p5wtuUXFq+RZS5Yswd/fn/79+yc/FhoaSuPGjXnmmWeoXr06NWrUYPLkyem2Ex0dTe/evalfvz61a9dm2rRpACQmJvL0009To0YNatasyZgxYxg9ejQHDx6kadOmNG3aFID58+dzyy23UKdOHbp06UJUVBQAc+fOpUqVKjRu3Jiffvoph14FERHPoJwtuUW3ypI8a/PmzdStWzfV4z/99BPh4eFs2rSJY8eOUb9+fZo0aXLRdt544w1uv/12vvjiC06dOkWDBg1o1qwZX3/9Nf/88w8bN27Ez8+PEydOULx4cUaNGsWSJUsoUaIEx44d4/XXX2fhwoUEBQXx9ttvM2rUKAYPHsyDDz7I4sWLuf766+natWtOvhQiInmecrbkFhWv4nGWL19O9+7d8fX1pVSpUtx2222sXbuWmjVrprn8/PnzmT59OiNHjgQgNjaWvXv3snDhQvr374+fn+tjULx48VTrrl69mi1bttCoUSMA4uLiuOWWW9i6dSsVKlSgUqVKANx33318+umnObG7IiIeTTlbspuKV8mzqlWrxpQpU1I9bq3NUjvWWn788UcqV66c6vGMZmux1tK8eXMmTpx43uPh4eGa6UVEJAXlbMktGvMqedbtt9/O2bNn+eyzz5IfW7t2LcWKFWPy5MkkJiZy9OhRfv31Vxo0aHDRdlq2bMmYMWOSE+jGjRsBaNGiBR9//DEJCQkAnDhxAoCQkBAiIyMBuPnmm1mxYgU7d+4EICYmhu3bt1OlShX++ef/27v3oKjq/4/jr/1iXkYKTVMxFFPIC7gsNyVvgQ46ij+8jkpN2njpomY/+46pf3ip/git0YbIL+NYiU4K8/OSZYLCOI6O4iDR8ksZcyshRSokUfulsi7n90fTFsqiFrAceD5m9o+zez5n32f/eM+Ls4d9n9d3330nSXc1SgBobejZaCqEVzRbFotFe/fuVU5Ojvr166eQkBCtXbtWzzzzjKxWq8LCwjR69GitX79ePXr08HicVatWyel0ymq1KjQ0VKtWrZIkzZ8/X71793Yfa8eOHZKkF154QePHj1dcXJwee+wxbd26VUlJSbJarYqJidHZs2fVvn17bd68WQkJCRoxYoQCAwOb5DMBgOaKno2mYnmQy/lRUVFGQUFBI5YDb7NYLMrLy1NMTIy3SwEalMVi+dIwjChv19GU6Nktn8ViUUJCgvbv3+/tUoAGVV/P5sorAAAATIPwCgAAANMgvKJZ8/X1bbRj9+nTR5cvX77v/e12uw4cONBo9QCA2dGz0RQIr8B9ohECgHnQs1suwitMp6KiQtOmTVN0dLSio6N1/PhxSVJ+fr6GDRum8PBwDRs2TN98842kukcK/uH9999XRESEBg8erLNnz3o8TnV1tVavXq3MzEzZbDZlZmZ6HGEIAPgTPRsNzjCM+35ERkYaaNkkGXl5ed4uw61jx453PZeUlGQcO3bMMAzDKC0tNQYMGGAYhmFcvXrVcDqdhmEYRk5OjjF16lTDMAxj06ZNxtSpU92vVVZWGoZhGIGBgUZKSophGIbxwQcfGPPmzav3OB9//LGxaNEidx0rV640tm/fbhiGYVy5csUIDg42fv311wY8ezQkSQXGA/S7lvCgZ7d8koyEhARvl+FGz0ZDqa9nM2ELppObm6vi4mL39rVr13T9+nVdvXpVc+bMkcPhkMVikdPpdO/vaaTg1KlTJUmRkZHas2ePJHk8zp08jTAcOHBgw580AJgUPRsNjfAK06mpqVFeXp46dOhQ6/lXXnlFcXFx2rt3r0pKShQbGyup/pGC7dq1kyT5+Pi4p7asWrWqzuPcyfAwwhAA8Cd6Nhoa97zCdMaOHavU1FT3tt1ul/T7X9+PP/64JGnr1q219q9rpKAnno7z1xGEkucRhgCAP9Gz0dAIr2jWfvvtNwUEBLgfGzZsUEpKigoKCmS1WjVo0CClpaVJkl5//XWtXLlSw4cPl8vlch/D00hBTzwdJy4uTsXFxe6b/z2NMASA1oqejabAeFjUwnhYtFSMh0VLxHhYtFSMhwUAAECLQHgFAACAaRBeAQAAYBqEV5hSdna2+vfvr6CgICUnJ9/1umEYWrJkiYKCgmS1WlVYWHjPtb/88ovi4+MVHBys+Ph4XblyRZKUk5OjyMhIDR48WJGRkTp8+HDjnyAAtEA3b97UkCFDFBYWppCQEK1Zs0aS9Pzzz+uJJ56QzWaTzWbTsGHD6lx//vx5DR06VMHBwZo5c6aqq6slSUeOHJGfn597/Ztvvtlk54SmR3iF6bhcLi1atEhZWVkqLi7Wzp07a/0AtiRlZWXJ4XDI4XBo8+bNevnll++5Njk5WWPGjJHD4dCYMWPcwbZr1676/PPP9fXXXys9PV3PPfdc054wALQQ7dq10+HDh1VUVCS73a7s7GydPHlSkvTOO+/IbrfLbrfrxIkTda5fvny5li5dKofDoc6dO+vDDz90vzZy5Ej3+tWrVzfJ+cA7CK8wnfz8fAUFBalv375q27atZs2addeM6n379mn27NmyWCyKiYlRVVWVysvL6127b98+zZkzR5I0Z84cffrpp5Kk8PBw9ezZU5IUEhKimzdv6tatW014xgDQMlgsFvn6+kqSnE6nnE6nx4EEdzIMQ4cPH9b06dMl1e7TaF0IrzCdsrIy9erVy70dEBCgsrKy+9qnvrU//fST/P39JUn+/v76+eef73rv3bt3Kzw83D3lBQDwYFwul2w2m7p166b4+HgNHTpUkrRs2TL31/7PPvvsXesqKyvVqVMn99jYO3t/Xl6ewsLCNH78eJ05c6ZpTgZewXhYmE5dv01851/unva5n7WenDlzRsuXL9ehQ4fus1IAwJ18fHxkt9tVVVWlKVOm6PTp05J+v23gj6uqdamvf0dERKi0tFS+vr46cOCAJk+eLIfD0TgnAK/jyitMJyAgQBcuXHBvX7x40f21/r32qW9t9+7dVV5eLkkqLy9Xt27dau03ZcoUbdu2Tf369WuU8wKA1qRTp06KjY1Vdna2x33GjRsnm82m+fPnq2vXrqqqqnKPjf1r/37kkUfctyNMmDBBTqdTly9fbvyTgFcQXqHbt29r7IT/0lOjRkuSZj77nIY9PUbff/+9lyurW3R0tBwOh86fP6/q6mplZGQoMTGx1j6JiYnatm2bDMPQyZMn5efnJ39//3rXJiYmKj09XZKUnp6uSZMmSZKqqqqUkJCgt99+W8OHD2/akwWAOiz59+uKGRUnSfriiy8UMypO/7N7r5erureKigpVVVVJkm7cuKHc3FwNGDDA4/4HDx6U3W7Xli1bZLFYFBcXp127dkmq3ad//PFH95XZ/Px81dTUqEuXLo18NvAWwit0+/ZtGTU1+rKwUJ1Hz9fP/1ej8ktleuihh7xdWp3atGmj1NRUjRs3TgMHDtSMGTMUEhKitLQ098zsCRMmqG/fvgoKCtKCBQu0adOmetdK0ooVK5STk6Pg4GDl5ORoxYoVkqTU1FR9++23euutt9z3Y9V1PywANJXOj/gq/8QxPfLUDHUcOFL/+9WX6vpoJ2+XdU/l5eWKi4uT1WpVdHS04uPjNXHiREm173m12Wzun8H6q3Xr1mnDhg0KCgpSZWWl5s2bJ0natWuXQkNDFRYWpiVLligjI+O+bwmD+VjquofEE+Zkt1wul0v//e9lSvvPJj01fKQ+37tLfn5+3i4LaDD1zcluqejZLduOHTs1/8WX9OijXZR78EC9VzABs6mvZxNeUcvp06c1YMAA939zAi0F4RUtUWlpqfz8/NSpU/O/6go8iPp6NgkFtYSGhnq7BADAfQoMDPR2CUCT455XAAAAmAbhFQAAAKZBeAUAAIBpEF4BAABgGoRXAAAAmAbhFQAAAKZBeAUAAIBpEF4BAABgGoRXE4mNjdXBgwdrPffee+9p4cKFde5fUlJyz6EDJSUl2rFjh3u7oKBAS5YskSRt3bpVixcvliSlpaVp27Zt/6R8SVKfPn10+fLlf3wcAGju6NlA4yC8mkhSUpIyMjJqPZeRkaGkpKS/fcw7G2FUVJRSUlLu2u+ll17S7Nmz//b7AEBrQ88GGgfh1USmT5+u/fv369atW5J+b2KXLl3SiBEjtGzZMoWGhmrw4MHKzMy8a21JSYlGjhypiIgIRURE6MSJE5KkFStW6NixY7LZbNq4caOOHDmiiRMn3rV+7dq1evfdd3Xp0iXZbDb3w8fHR6WlpaqoqNC0adMUHR2t6OhoHT9+XJJUWVmpsWPHKjw8XC+++KIMw2jETwgAmg96NtA4CK8m0qVLFw0ZMkTZ2dmSfv8LfubMmdqzZ4/sdruKioqUm5urZcuWqby8vNbabt26KScnR4WFhcrMzHR/zZScnKyRI0fKbrdr6dKl96yhZ8+estvtstvtWrBggaZNm6bAwEC9+uqrWrp0qU6dOqXdu3dr/vz5kqQ33nhDI0aM0FdffaXExET98MMPDfypAEDzRM8GGkcbbxeAB/PH11CTJk1SRkaGPvroI23fvl1JSUny8fFR9+7d9fTTT+vUqVOyWq3udU6nU4sXL5bdbpePj4/OnTv3j+o4fvy4tmzZomPHjkmScnNzVVxc7H792rVrun79uo4ePao9e/ZIkhISEtS5c+d/9L4AYCb0bKDhEV5NZvLkyXrttddUWFioGzduKCIi4r5uyt+4caO6d++uoqIi1dTUqH379n+7hvLycs2bN0+fffaZfH19JUk1NTXKy8tThw4d7trfYrH87fcCADOjZwMNj9sGTMbX11exsbGaO3eu+6b/UaNGKTMzUy6XSxUVFTp69KiGDBlSa93Vq1fl7++vf/3rX9q+fbtcLpck6eGHH9b169fv+/2dTqdmzJihdevW6cknn3Q/P3bsWKWmprq37Xa7u7ZPPvlEkpSVlaUrV678vRMHABOiZwMNj/BqQklJSSoqKtKsWbMkSVOmTJHValVYWJhGjx6t9evXq0ePHrXWLFy4UOnp6YqJidG5c+fUsWNHSZLValWbNm0UFhamjRs33vO9T5w4oVOnTmnNmjXufwC4dOmSUlJSVFBQIKvVqkGDBiktLU2StGbNGh09elQRERE6dOiQevfu3cCfBgA0b/RsoGFZHuQ/CaOiooyCgoJGLAcAGofFYvnSMIwob9fRlOjZAMyqvp7NlVcAAACYBuEVAAAApkF4BQAAgGkQXuGWnZ2t/v37KygoSMnJyd4uBwBQD3o2WivCKyRJLpdLixYtUlZWloqLi7Vz585aP2ANAGg+6NlozQivkCTl5+crKChIffv2Vdu2bTVr1izt27fP22UBAOpAz0ZrRniFJKmsrEy9evVybwcEBKisrMyLFQEAPKFnozUjvEKSVNfv/TIiEACaJ3o2WjPCKyT9/lf7hQsX3NsXL15Uz549vVgRAMATejZaM8IrJEnR0dFyOBw6f/68qqurlZGRocTERG+XBQCoAz0brVkbbxeA5qFNmzZKTU3VuHHj5HK5NHfuXIWEhHi7LABAHejZaM0sdd034wlzsgGYVX1zslsqejYAs6qvZ3PbAAAAAEyD8AoAAADTILwCAADANAivAAAAMA3CKwAAAEyD8AoAAADTILwCAADANAivAAAAMA3CKwAAAEyD8AoAAADTILwCAADANAivAAAAMA3CKwAAAEyD8AoAAADTILwCAADANCyGYdz/zhZLhaTSxisHABpNoGEYj3m7iKZEzwZgYh579gOFVwAAAMCbuG0AAAAApkF4BQAAgGkQXgEAAGAahFcAAACYBuEVAAAApkF4BQAAgGkQXgEAAGAahFcAAACYBuEVAAAApvH/323Cc5pIAZsAAAAASUVORK5CYII=\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": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, 'Percent of Incoming PFAS that \\n remains in the landfill (%)')"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaoAAAEKCAYAAABDkxEYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxcdb3/8den+1660zYTaKEUytJtUhZBWRUrIFuTKioKws8drz+vVy5XAVdcr179cbEiV66KTUrZQbZSiihI0r2UFiqUTtrSFOi+J/n8/vieodPpzCRNk8wkeT8fjzwy853vOec7Z+aRT845n/P5mrsjIiJSqDrlewAiIiK5KFCJiEhBU6ASEZGCpkAlIiIFTYFKREQKmgKViIgUtIIJVGZ2l5nVmNmylLaBZvaUmb0W/R4QtQ8ys7lmtt3Mfp1jnRmXFxGRtqNgAhXwe+DCtLZvAnPcfQwwJ3oOsBv4FvD1BtaZbXkREWkjCiZQuftzwLtpzR8F7o4e3w1cGvXd4e7PEwJWLhmXFxGRtqNLvgfQgGHuvh7A3deb2dCWWt7MrgeuB+jdu/fk448/vqljFhHpcObPn/+2uw9piXUXeqBqNe4+A5gBEI/HvaqqKs8jEhFpO8zszZZad8Gc+stig5kNB4h+17Ty8iIikmeFHqgeAq6OHl8NPNjKy4uISJ4VTKAysz8DLwBjzazazK4FbgMuMLPXgAui58n+q4GfA5+O+o+L2u80s3jULevyIiLSNhTMNSp3/1iWl87L0v/oLO2fTXn8TrblRUSkbSiYIyoREZFMFKhERKSgKVCJiEhBU6ASEZGCpkAlIiIFTYFKREQKmgKViIgUNAUqEREpaApUIiJS0BSoRESkoClQiYhIQVOgEhGRgqZAJSIiBU2BSkRECpoClYiIFDQFKhERKWgKVCIiUtAUqEREpKApUImISEFToBIRkYKmQCUiIgVNgUpERAqaApWIiBQ0BSoRESloClQiIlLQFKhERKSgKVCJiEhBK5hAZWZ3mVmNmS1LaRtoZk+Z2WvR7wEpr91oZqvMbKWZfSjLOm8xs7Vmtij6mdoa70VERJpPwQQq4PfAhWlt3wTmuPsYYE70HDMbB0wHToyWud3MOmdZ73+6+4To57EWGbmIiLSYgglU7v4c8G5a80eBu6PHdwOXprTPdPc97v4GsAqY0ioDFRGRVtVgoDKzUY1payHD3H09QPR7aNQ+Ekik9KuO2jL5kpktiU4tDsjSR0REClRjjqhmZ2i7t7kHcogsQ5tnaPtv4BhgArAe+FnWFZpdb2ZVZla1cePG5hmliIgcti7ZXjCz4wnXgPqb2eUpL/UDerT0wCIbzGy4u683s+FATdReDcRS+hUB69IXdvcNycdm9lvgkWwbcvcZwAyAeDyeKeiJiEge5DqiGgtcBBwBXJzyMwm4ruWHBsBDwNXR46uBB1Pap5tZ9+g05BjgpfSFo+CWdBmwLL2PiIgUtqxHVO7+IPCgmZ3u7i+09EDM7M/A2cBgM6sGbgZuAyrM7FpgDTAtGtvLZlYBLAdqgS+6e120njuBO9y9CvixmU0gnBZcDfyfln4fIiLSvMw991kuM+sBXEs4DfjeKT93v6Zlh5Y/8Xjcq6qq8j0MEZE2w8zmu3u8JdbdmGSKPwBHAh8C5hGuB21ricGIiIika0ygOtbdvwXscPe7gY8AJ7fssERERILGBKp90e/NZnYS0B84uqGFzKx3jmoRIiIijZI1mSLFjOhG2W8Rsu36AN9O72RmnQhlja4CSoA9QHcz2wg8Bsxw99eaa+AiItIxNBio3P3O6OE8YHSOrnOBp4EbgWXuXg+hsCxwDnCbmd3v7n88vCGLiEhH0mCgMrPuwBWE033v9Xf376R1Pd/d96W14e7vEqpbzDazroc1WhER6XAac+rvQWALMJ9wOi+j9CAVpbV/AugJ3OPu72QKZCIiIrk0JlAVuXv69BuN8UtgAbAbeAA4qwnrEBGRDq4xWX9/N7MG09HN7B4zOyalaSDwJ+DPgKqWi4hIk+QqSruUUHqoC/AZM3udcOrPAHf3U9IW+Q/ge2a2Dvgu8FNClmAP4JbmH7qIiHQEuU79XXQoK3L314GPm9mZQDnwKHBBsgafiIhIU2Q99efub7r7m8D3ko9T29L7m9kAM/siMA4oJSRgPGFmhxTwREREUjXmGtWJqU+iahOTM/R7gHBqsAfwB3f/X8K0IJPN7KHDHaiIiHRMua5R3Qj8O9DTzLYmm4G9RBMMphkE3ENIR/8UgLvvAm5NmxdKRESk0XLNR/VD4Idm9kN3v7ER67oZeAqoA76Ztq71hzVKERHpsBpTQqkxQQp3n02oQCEiItJsGnONqlHMbEZUXT3Ta73N7Bozu6q5ticiIh1DYypTNNbtwLejm4OXARsJiRVjgH7AXYQbgEVERBotVzJFL2Bfsj6fmY0FpgJvuvt96f3dfRFQamZ9gDgwHNgFvOLuK1ti8CIi0v7lOqJ6HLgWeM3MjgVeIBwRXWRmJdmuXbn7duDZ5h6oiIh0TLmuUQ1ImejwauDP7v5l4MMcYtUKERGRpsoVqDzl8bmE1HPcfS9Q35KDEhERScp16m+Jmf0UWAscCzwJYGZHtMbAREREIHegug64gTCz7wfdfWfUPo5QGf0AZvYwBx6FHcDdL2n6MEVEpKPKVZliF3BbhpcSwPsytB8UvERERA5Xo+6jMrPBwDTgY8BI4P70Pu4+r3mHJiIikvs+qr7AZcDHgeMIwWm0uxdl6Z+caDGjDBMtioiINCjXEVUN8BJh5t7n3d3N7LIc/ZWyLiIizS5Xevq/E0og/Tdwo5kdk2tF6ZMrZphsMSczu8vMasxsWUrbQDN7ysxei34PSHntRjNbZWYrzexDWdaZdXkREWkbcs3w+5/ufipwCWEeqgeAEWb2b2Z2XHp/M3s++r3NzLam/GxLmc8ql98DF6a1fROY4+5jgDnRc8xsHDCdMKnjhcDt0YSO6TIuLyIibUeD1dPd/XV3/767nwyUAP2Bv2Tompwssa+790v56evu/RqxneeAd9OaPwrcHT2+G7g0pX2mu+9x9zeAVcCUDKvNtnxOW3btY09tXWO6iohIC8saqMzstPQ2d1/q7v/u7plOA86KlpvTjOMblpx0Mfo9NGofSUiTT6qO2hq7/EHM7HozqzKzqjXv7uT0Hz7D9x5ZzmsbtjXLGxERkabJlUxxOzAJwMxecPfTG1hXJzO7GTjOzL6W/qK7/7zpwzyIZWjLmnHYGO4+A5gBMPakCX7a6IHc/cJq7nz+DSYVH8H0kmI+cspwendvzplRRESkIblO/aUGgx6NWNd0YDch+PXN8NMUG8xsOED0uyZqrwZiKf2KgHWHsHxOfXt04farJvPCjedx09QT2LJrH9+YvYQp33+aG+9bwqLEZtwPKy6KiEgj5To86BRlyXVKefxe8HL3A64nRXNO/cjMlrh7pmtYTfEQoXL7bdHvB1Pa7zGznwMjCJMzvnQIyzfK4D7due79o/nsWaNYsGYTM19K8MDCdfz5pQTHH9mXspIYl00cyRG9ujXt3YmISIMs25GBma0mVEnPeJrN3Uen9T/odF/aAjlP/ZnZn4GzgcHABuBmQqZhBVAMrAGmJQOkmd0EXAPUAl9NBkczuxO4w92rzGxQtuVzicfjXlVVlfG1bbv38fDi9ZRXrmFx9Ra6denEhSceyfSSGKeNHkSnTpl2l4hI+2Zm89093iLrbq5TWNH1KYCxhOzAh6LnFwPPuftnm2VDrSBXoEq1fN1WKqoS3Legmq27a4kN7ElZPMaVk2Mc2b8xZ0tFRNqHvAQqMxtKuOn3WGAJcJu7N3g/lJk9CVzh7tui532BWe6efo9UwWpsoErava+OJ15+i5kvJXjh9XfoZHDO2KGUlcQ45/ihdO3c4F0AIiJtWr4C1ePAfOA5Qnmkvu7+6QZXaLYCGO/ue6Ln3YHF7n58cw26pR1qoEq1+u0dlFcluHd+NRu37WFI3+5cObmI0niMUYN7N/NIRUQKQ74C1SJ3n5DyfIG7T2pwheHaUSmhiK0TCttWuPsPmmfILe9wAlVSbV09c1dupLwywdyVNdTVO6eOGsj0KTE+fNJwenTNVEhDRKRtylegWkxIbkhmB8xNfZ4rKcHMJgNnRk+fc/eFzTPc1tEcgSrVhq27uXd+NRVVCd58Zyd9e3ThsokjKSuJceKI/s22HRGRfMlXoFrNIWT9pS3bGRhGSvq7u685rJG2ouYOVEn19c4/3niX8so1PLbsLfbW1nPyyP6UlsT46IQR9OvRtdm3KSLSGtpE1t97KzT7MiG1fANQRwh03pbmo2qpQJVqy859PLBoLTMrE7yyfis9unZi6snDmV5STMnRAzBTmruItB1tLVCtAk5193eadcWtqDUCVZK7s2ztVmZWruGhRevYtqeW0YN7U1oS44pJRQzp271VxiEicjjaWqCaC1zg7rXNuuJW1JqBKtXOvbU8tvQtKioTvLT6Xbp0Ms47YSjTS4p5/3FD6KybiUWkQLVkoGqJCquvA8+a2aPAnmRjMxelbZd6devClZOLuHJyEatqtr93M/ETL2/gyH49mBYPae6xgb3yPVQRkVbT4BGVmQ3M0LzN3fdl6X9zpnZ3v/XQh5cf+TqiymRfXT1zXqmhvHIN817dSL3DmccOprQkxgfHDVOau4gUhLye+ouy/2LAJkJixBHAekIl8uvcfX5LDCyfCilQpVq3edd7ae7Vm3ZxRK+u76W5H39kg3NTioi0mHwHqjuA+939iej5BwnTv1cAv4ymq0/tPwT4BmGa+PcK3rn7uc079JZTqIEqqb7e+fs/32Fm5RqefHkDe+vqGR87guklMS4eP4I+mjNLRFpZvgNVVfrGk23p1Sui154EyoGvA58jTK+x0d3/rXmH3nIKPVCl2rRjL/cvXMvMyjW8umE7vbp15iMnD2f6lBiTipXmLiKtI9/JFO+a2b8BM6PnZcCm6Kbe+gz9B7n778zsBnefB8wzs3nNNF5JM6B3N645cxSfed/RLEpsprwywcOL1zFrfjXHDu3D9GjOrEF9lOYuIm1TY46oBhNu4D2TcI3qeeBWYAtQ7O6r0vq/6O6nmdkTwH8RZt69192PaYHxt4i2dESVyY49tTy6ZD0zK9ewYM1munY2Lhg3jNJ4jLPGKM1dRJpfW7uP6iLgr4QEjF8B/YBb3f2hnAsWkLYeqFK9tmEbMytDmvumnfsY0b8H0+IxpsWLKBqgNHcRaR75vkZ1HOF609EcWLuvzSRHHKr2FKiS9tTW8fTyGmZWruH5VW8DIc29rCTGBeOG0b2L0txFpOnyHagWA3cQ5qaqS7anp6Wb2a8I03pk5O5fOayRtqL2GKhSVW/ayayqau6dX83azbsY0Ksrl08qoqwkxnHD+uZ7eCLSBuU7UM1398kNrsjs6lyvu/vdhzi2vGnvgSqprt55ftXbVFQmeHL5W+yrcyYWH0FZPMZFSnMXkUOQ70B1C+Hm3vs5sCRS1vmo2rqOEqhSvbN9D/cvXEt5ZYLXakKa+8WnjKC0JMak4iOU5i4iOeU7UL2RoTnnfFRtXUcMVEnuzoI1m6moTPDwknXs3FvHmKF9KFOau4jk0Kay/tqDjhyoUm3fU8sji9dRXpVgYUqae1lJMWceO1hp7iLynnzN8Huuuz9jZpdnet3d72uJARUCBaqDrXxr23vV3JXmLiLp8hWobnX3m83sfzK87O5+TZbljgP+Gxjm7ieZ2SnAJe7+vWYbdQtToMpuT20dTy3fQHll4oA09+klxZw/bqjS3EU6qDZ16i8ql/SvwG/cfWLUtszdT2rWDbUgBarGSaa5z6pKsG7LbqW5i3Rg+U6mOAL4FAff8Jvxvigzq3T3EjNbmBKoDipeW8gUqA5NMs29vHINTy3foDR3kQ4o30VpHwNeBJaSuQhturfN7Biim3/N7ErC/FXSTnXuZHzguCF84LghB6S5f/O+pXznkeVcdMpwykqKleYuIk3SmCOqBe4+qdErNBsNzADOIEy2+AbwCXdf3eRBmt0AXEcoivtbd/+FmY0nVMzoA6wGrnL3rRmWXQ1sI1TVqG1MxNcR1eHLlOauau4i7Ve+T/39C7AdeIRDuOHXzHoDndx922EN0OwkwhQjU4C9wOPA54F7gK+7+zwzuwYY5e7fyrD8aiDu7m83dpsKVM1Lae4i7V++T/3tBX4C3MT+Wn4OZLzh18y6A1cQXdNKnupx9+80cYwnAC+6+85o/fOAy4CxwHNRn6eAJ4CDApXkX5/uXZg+pZjpU4p5dcM2yqNq7o8tfYsR/XtwZTzGtMlFxAYqzV1EDtaYI6p/Aqc29ojEzB4nzFWVXsT2Z00aoNkJwIPA6cAuYA5QBUwGfuTuD5rZ1whTiRyUahZV1thECK6/cfcZWbZzPXA9QHFx8eQ333yzKcOVRkpWcy+vSvDX1zYCIc29NB7jgyeqmrtIW5PvU38PAdOTRzQNrrAFUtHN7Frgi4RTkMsJAes3hIkZBwEPAV9x90EZlh3h7uvMbCjhyOvL7v5cer9UOvXXuqo37eTe+dXMqtpfzf3SiSMpK4lx/JH98j08EWmEfAeq+4ETgbkceI0qW3r6DOBX7r60GceZuv4fANXufntK23HAH919SgPL3gJsd/ef5uqnQJUfdfXO31a9TXlKNffxsSOYXhLjolOG07dH13wPUUSyyHegyjh9R/q0HWa2lHB6rQswBnidENgsdPdTmjxIs6HuXmNmxcCThNOAXaO2TsDvgWfd/a605d5L6IgePwV8x90fz7U9Bar8e3fH3ijNfQ2vbthOz66dozT3GJOPGqA0d5ECk/fKFGbWDTguerrS3fdl6HNUrnW4e5Mv+pjZXwmn+PYBX3P3OVHK+hejLvcBN7q7m9kI4E53nxqlyt8f9ekC3OPu329oewpUhcPdWZTYTHllgocXr2PH3jqOGdKbspIYl08qYrDS3EUKQr6PqM4G7ibcq2RADLg623UeM/uDu3+yobZCpkBVmHbsqeXRpespr0ww/81NdOlknH/CMMpKYrz/uCFKcxfJo3ynp/8M+KC7r4wGcxzwZ0LWXSYnpj4xs845+oo0Wu/uXSiNxyiNx1hVk0xzX8vjL7/Fkf16MC1eRGk8pjR3kXamMUdUS9KvL2VpuxH4d6AnkMwQNMJ9WDPc/cZmG3UL0xFV27G3tp5nVmxgZmWC517dSL3DGccMoqwkxodOPJIeXZXmLtIa8n3q7y5CksQfoqargC7u/pks/X/YloJSJgpUbdP6Lbu4t6qa8qoE1Zt20b9nVy6dMIKykmLGjVCau0hLyneg6k5IWjiTcIT0HHC7u+/JuWAbpkDVttXXOy+8/g4zKxM8sewt9tbVc/LI/pSVxLhkwgj6Kc1dpNnlO1D1Bna7e130vDPQvbE3ALdFClTtx+ade3lg4VpmViZY8dY2enTtxNSThlNaEuPUUQOV5i7STPIdqF4Eznf37dHzPsCT7n5GSwyoEChQtT/uztK1W5hZmeDhRevYtqeWUYN7My1exJWTihjar0e+hyjSpuU7UB006WFDEyFGR13DOHCixTWHOdZWo0DVvu3aW8djUZr7S6vfpXMn45yxQykriXHO2CF06dwp30MUaXPynZ6+w8wmufuCaDCTCbX2MjKzLwM3AxvYP9GiA02uTCHSnHp268wVk4u4YnIRr2/cTkVVNbMXVPP0KxsY0rc7V04Oae6jBvfO91BFhMYdUZUQ5oNaFzUNB8rcfX6W/qsI1dbfac6BtiYdUXU8++rqmbuihoqqBHNXbqSu3pkyaiBl8RhTTx5Oz25KcxfJpRBKKHUlzP9kwIpMJZRS+s4FLnD32mYbZStToOrYNmzdzewF1VRUJlj9zk76du/CJRNGUFYS4+SR/ZWAIZJBIQSqM4gmQky2ufv/Zun7O0JQe5QDq63//DDH2moUqARCAsY/3niXisoEjy1bz+599Rx/ZF+ml8S4dOJIjujVLd9DFCkY+U6m+ANwDLCI/RMheo5pPm7O1O7utx7GOFuVApWk27JrHw8tXkdFZYKla7fQrUsnPnTikZTFY5xxzCA6qc6gdHD5DlSvAOO8MYde7YQCleSyfN1WKqoS3L9wLVt27aNoQE9K4zGunFzEiCN65nt4InmR70A1izB77voG+v3C3b9qZg8TsvwO4O6XHNZIW5EClTTG7n11PPHyW1RUJfjbqncwg/ePGUJZSYzzTxhGty5Kc5eOI9+Bai4wAXiJA685XZLWb7K7zzezD2Raj7vPO/zhtg4FKjlUiXd3Mqsqwaz51azfspuBvbtx2cSRlJXEOG5Y33wPT6TF5TtQtfnAc6gUqKSp6uqd517bSEVlgqdf2cC+Omdi8RGUxWNcNH4Efbo35tZFkbYn71l/HY0ClTSHd7bv4f6FaymvTPBazXZ6devMR04eTllJjMlHDVCau7QreQlUZraNDNeaCPdSubu323kTFKikObk7CxObKX8pwSNL1rFjbx3HDOlNWUmMyycVMbhP93wPUeSwtdkjKjPrBPRx960ttpEWoEAlLWXHnloeXbKe8qoE89/cRJdOxnknhDqD7x+jOoPSdrWpQGVm9wCfI9xzNR/oD/zc3X/SrBtqQQpU0hpW1WwLdQbnV/POjr0M67e/zuBRg1RnUNqWthaoFrn7BDO7CpgM/BswP33q+kKmQCWtaW9tPc+s2EB5ZYJ5r26k3uG00QMpK4nx4ZOG06Or6gxK4cvXNaruTZnF18xeJqSz3wP82t3nmdlidx9/eENtPQpUki/rt+xi9vxqKqqqWfPuTvr26MJHJ4xgekkxJ43sn+/hiWSVr2k+XgAmmdkf3P2Th7DO3wCrgcXAc2Z2FNCmrlGJ5Mvw/j350rlj+MLZx/LiG+9QUZlgVlU1f3xxDeOG96M0XqQ6g9Lh5DqiWgb8BPg28K/pr7v7fY3eiFmXtlRNXUdUUki27NzHQ4vXMrMywcvrtqrOoBSkfJ36OxO4CigFHkp72d39mizLdQeu4OBq699phvG2CgUqKVTL1m6hoirBAwvXsnV3LUUDejJtcoxpcdUZlPzKd2WKa939d41eodnjwBZCxl+y2jru/rOmDrK1KVBJoUvWGSyvTPD3f4Y6g2eNGUJZPMYF41RnUFpfvgNVN0K6+fujpnnAHdkmTzSzZe5+UrOOspUpUElbsuadncyaH65lvbU11Bm8dEKoMzj2SNUZlNaR70B1J9AVuDtq+iRQ5+6fzdJ/BvArd1/abIM0uwG4jlAV47fu/gszGw/cAfQhJG9clenGYjO7EPgl0Bm4091va2h7ClTSFmWqMzg+FuoMXjx+OH17dM33EKUdy3egOii1PFe6uZktB44F3iBUW0+WXGrSfVRmdhIwE5gC7AUeBz5PSH//epT+fg0wyt2/lbZsZ+BV4AKgGqgEPubuy3NtU4FK2rr0OoM9u3ZmalRnsORo1RmU5pev9PSkOjM7xt3/GQ1mNCnXnjL4cLOMbL8TgBfdfWe0/XnAZYTp7p+L+jwFPAF8K23ZKcAqd389WnYm8FEgZ6ASaesG9enOZ88azbVnjmJhYjMVlQkeXryO2QuqGT24N9PiMa6YPJKhfXvke6giDWpMoPpXYK6ZvU44OjoK+Ex6JzPrF51629a8Q2QZ8H0zGwTsAqYCVVH7JcCDwDQglmHZkUAi5Xk1cGqmjZjZ9cD1AMXFxc01dpG8MjMmFQ9gUvEAvnXROB5dup6KygQ/enwFP31yJeeMDXUGzxmrOoNSuBoMVO4+x8zGEI5gDFiRpWLFPcBFhGw/j/q+txpgdFMG6O6vmNmPCEdN2wk3EtcC1wD/ZWbfJqTP782weKbzGxnPdbr7DGAGhFN/TRmrSCHr3b0LpfEYpfEYq2q2M2t+gtnz1/L0KxsY0rc7V0wqojRexOghffI9VJEDtLn5qMzsB0C1u9+e0nYc8Ed3n5LW93TgFnf/UPT8RgB3/2GubegalXQU++rqmbuihoqqBHNXbqSu3ply9EBKS2JMPflIenXTRI/SOG2qKC2AmQ0AxgDvnQB39+eyL9Hg+oa6e42ZFQNPAqcDXaO2TsDvgWfd/a605boQkinOA9YSkik+7u4v59qeApV0RBu27mb2gmpmVVXzxts76NO9CxePH0FZSYzxRf2VgCE5talAZWafBW4AioBFwGnAC+5+7mGs86/AIGAf8LXodOQNwBejLvcBN7q7m9kIQhr61GjZqcAvCOnpd7n79xvangKVdGTuzktvvEt5VYLHlq5n9756xg7ry7R4EZdPKmJgb9UZlIPlOz19jruf11BbymtLgRJCpt4EMzseuNXdy5pr0C1NgUok2LZ7Hw8vDhM9Lk5spmtn44JxwyiNxzhrzBA6q86gRPKSnm5mPYBewODoVF7yG9kPGJFjnbvdfbeZJacKWWFmY5tvyCLSWvr26MrHTy3m46cWs+KtrVRUVnP/wmoeW/oWI/r34MrJRUyLx4gN7JXvoUo7lqso7Q3AVwlBaS37A9VWQnWIX2dZ7n5C+vpXgXOBTYTrSVObd+gtR0dUItntqa3j6eU1lFcl+OtrG3GHM44ZRFlJjA+deKQmeuyg8n3q78vu/qsmrdzsA4Sp6B9390zp4wVJgUqkcdZu3sW9VdXMmp+getMu+vXowqUTR1Iaj2mixw4m78kUZnYGB0/b8b85+g8g3ICb2n/B4Qy0NSlQiRya+nrnhdffobwyweMvv8Xe2npOHNGP0niMSyeMpH8v1Rls7/J9RPUH4BhCBl+ydJK7+1ey9P8u8GngdaA+pX+Ts/5amwKVSNNt3rmXBxeto7wywfL1YaLHC088krKSGKeP1kSP7VW+A9UrwDhvZB67ma0ETm5Lp/rSKVCJNI/0iR5jA8NEj1dO1kSP7U2+A9Us4Cvuvr5RKzSbDXze3WuaYXx5oUAl0rxyTfR4/rihdO+iBIy2Lt+Bai4wAXiJMG0HAO5+SZb+cUKh2GWN6V+IFKhEWk7i3Z3Mqkowa34167fsZkCvrlw2sYjSkiKOP7JfvocnTZTvQPWBTO3uPi9L/5eB3wBL2X+NKmv/QqRAJdLy6uqdv762kVlV1QMd3LkAABcuSURBVDy5/K0w0WNRf0pLYlw8fgT9NNFjm1IIWX9HAWPc/Wkz6wV0dveM03mY2Tx3zxjc2goFKpHW9e6Ovdy/cC0VlQlWbthGj66dmHrycErjMU4dNVB1BtuAfB9RXUeYp2mgux8TTflxR44SSj8nnPJ7iANP/Sk9XURycneWVG+hvCrBw4vWsW1PLUcP6sW0eEjAGNZPEz0WqnwHqkWEmXL/4e4To7al7n5ylv5zMzQrPV1EDsmuvXX8Zdl6yisT/OONd+lkcM7YoUyLxzjvhKF01USPBSXfU9Hvcfe9yUPvaOqMrNHN3c9pprGJSAfWs1tnLp8UKravfnsHFVUJ7p1fzZwVNQzu043Lo4kejx3aN99DlRbWmCOqHwObgU8BXwa+ACx395uy9B8G/AAY4e4fNrNxwOnu/rtmHXkL0hGVSGGqravnudc2Ul6ZYM4rNdTWO5OKj6A0HuOi8SPo010TPeZLvk/9dQKuBT5IKEz7BGG+p4wLmtlfgP8BbnL38dER2MJspwoLkQKVSOF7e/se7l+wlvKqBKtqttOrW2c+cvJwykpiTD5qgBIwWlm+A1VvwtQdddHzzkB3d9+ZpX+lu5eY2cKUa1qL3H1CM4+9xShQibQd7s6CNZuZVZXg4cXr2LG3jtFDelMaj3H5pJEM7asEjNbQkoGqMVcj5wCptU56Ak/n6L/DzAYRXccys9OALU0eoYhIDmbG5KMGcNsVp/DSTefz4ytPYVDvbtz2lxWc/sNn+OzdVTy1fAO1dfUNr0wKUmNO6PZw9+3JJ+6+PbqXKpuvEVLTjzGzvwFDgCsPb5giIg3r3b0LpfEYpfEY/9y4nYqqBLPnr+XpVzYwpG93Lp8UpiA5ZkiffA9VDkFjTv39Dfhy8j4oM5sM/NrdT8/QtxNwGqHc0ljCNa2V7r6vuQfeknTqT6T92FdXz7MrN1JRleCZFTXU1TvxowZQWhLjIycPp7cSMJpFvq9RxYFyYF3UNBwoc/f5Wfq/kCmItSUKVCLtU8223dy3IFTAeP3tHfTu1pmLThlBaUmMScVHKAHjMOQtUKUcIVWy/whpRa4jJDO7FVgC3NfYqUEKjQKVSPvm7sx/cxPllQkeXbqenXvrOHZoH0rjRVw2sYghfbvne4htTr6PqA7pCMnMtgG9gVpgNyG4ubu3mbLIClQiHcf2PbU8uiRM9LhgzWa6dDLOPX4oZSUxPnDcELqoAkaj5DtQtfkjpEOlQCXSMa2q2UZFVTX3Lajm7e17Gdq3O1dMLmLa5CJGKwEjp3wHquQRUh2wizZ4hHSoFKhEOrZ9dfU8s6KGisoEc1fWUO8w5eiBlJbEmHrykfTqpgSMdHmf5qOjUaASkaQNW3cze0E1s6qqeePtHfTp3oWLxw9nWjzGxJgSMJLyfURlwFXAKHf/rpnFgOHu/lJLDKgQKFCJSDp356U33qWiqprHlq5n1746xgztQ2k8xmWTRjK4T8dOwMh3oPpvwky957r7CWY2AHjS3UtaYkCFQIFKRHLZtnsfjyxZT0VVgoVRAsZ5J4QEjPeP6ZgJGPme5uNUd59kZgsB3H2TmXVricFkY2Y3ANcRro/91t1/YWYTgDuAHoQMwy9kOsozs9XANsI1ttqW2pEi0nH07dGVj00p5mNTinl1wzZmVSW4b8Fannh5A8P6deeKSUVMi8cYNbh3vofaLjTmiOofwBlAZRSwhhCOqCa2ygDNTgJmEiZv3As8Dnwe+H/Af7r7X8xsKvANdz87w/Krgbi7v93YbeqISkQO1d7aKAGjKsGzyQSMUQMpjXeMBIx8H1H9F3A/MNTMvk+o2/cfLTGYLE4AXkxWazezecBlhKK3yczD/uyvnCEi0uq6denEhScdyYUnHXlAAsbXZy3mlode5uLxwymNx5igBIxD1qisPzM7HjiPcOptjru/0tIDS9n2CcCDwOmE9Pg5QBVwO2FuLCNUgT/D3d/MsPwbwCZCYPuNu89oaJs6ohKR5uDuVK4OFTCSCRjHDQsJGJdObF8JGHlJpjCzHsDngGOBpcDv3L22JQbREDO7FvgisB1YTghYnYF57j7bzEqB6939/AzLjnD3dWY2FHiKUGD3uQz9rgeuByguLp785psHxTwRkSbLlIBx/gnDKC0pahcJGPkKVOXAPuCvwIeB1e7+1ZYYxKEwsx8A1cAPgSPc3aMU+i0N3YRsZrcA2939p7n66YhKRFrSaxu2URElYLyzY+97CRil8RhHt9EEjHwFqqXJ6eOj6eRfcvdJLTGIhpjZUHevMbNi4EnCacC/A59392fN7Dzgx+4+OW253kAnd98WPX4K+I67P55rewpUItIakgkYs6pSKmCMGkhZPMaH21gCRr6SKd6rkO7utXm++Dc7mjV4H/DFKEX+OuCXURDdTXTazsxGAHe6+1RgGHB/NPYuwD0NBSkRkdaSLQHj/85azM1KwHhPriOqOmBH8ilhCvqdqNafiEiLSSZgVFQleHTJ/goYZSWFnYChWn+tTIFKRArBtt37eDRKwFhQ4BUwFKhamQKViBSa1zZsY9b8/VOQFFoFDAWqVqZAJSKFKnUKkmdf3UhdvRfEFCQKVK1MgUpE2oINW3dz34K1zKpK8PrbO+jdrTMXjx9BaUnrT0GiQNXKFKhEpC1xd6re3ERFZYJHl65n5946jh3ah9J4EZdNLGJI35ZPwFCgamUKVCLSVm3fU8ujS9ZRUVXN/Dc30aWTce7xQymNxzh7bMslYChQtTIFKhFpD1bVbGNWVTWzF6zl7e17GNI3mYBRxDFD+jTrthSoWpkClYi0J/vq6pm7ooaKqmrmrqyhrt6JHzWA0niMj5wynN7dDz8BQ4GqlSlQiUh7VbMtJGBUVCV4feMOenXrzEWnhAoYk48a0OQEDAWqVqZAJSLtnbuzYE2YguSRJSEBY/SQ3pTGY1w+cSRD+/U4pPUpULUyBSoR6Uh27Knl0aXrmVWVoHL1Jjp3Ms4ZO4Rp8RjnHj+Uro1IwFCgamUKVCLSUb2+cTuz5lcze341Ndv2MLhPNy6fVERpvIhjh/bNupwCVStToBKRjq62rp55r26koirBnFdqqK13JhYfQWk8xkWnDKdvj64H9FegamUKVCIi+729fQ8PLFxLeWWC12q207NrZ6aePJzSeBFTRg3EzPI2H5WIiAiD+3Tns2eN5tozR7EosZmKqmoeXryO2QuqOXpQL6bFYy26fQUqERFpFDNjYvEAJhYP4NsXjeMvy9ZTXpngJ0+sbNnt6tTfwXTqT0Sk8Va/vYNRQ/q02Km/wpl1S0RE2qSjW3g+LAUqEREpaApUIiJS0BSoRESkoClQiYhIQVOgEhGRgqZAJSIiBU2BSkRECpoClYiIFDQFKhERKWgKVCIiUtDaRKAysxvMbJmZvWxmX43aJpjZi2a2yMyqzGxKlmUvNLOVZrbKzL7ZuiMXEZHDVfCBysxOAq4DpgDjgYvMbAzwY+BWd58AfDt6nr5sZ+D/AR8GxgEfM7NxrTV2ERE5fAUfqIATgBfdfae71wLzgMsAB/pFffoD6zIsOwVY5e6vu/teYCbw0VYYs4iINJO2MB/VMuD7ZjYI2AVMBaqArwJPmNlPCQH3jAzLjgQSKc+rgVMzbcTMrgeuj55uN7OWnWClbRgMvJ3vQRQI7Yv9tC/2077Yb2xLrbjgA5W7v2JmPwKeArYDi4Fa4PPAv7j7bDMrBX4HnJ+2uGVaZZbtzABmNNvA2wEzq2qp+WXaGu2L/bQv9tO+2M/MWmwSv7Zw6g93/527T3L39wPvAq8BVwP3RV1mEU7zpasGUudILiLzKUIRESlQbSJQmdnQ6HcxcDnwZ0LA+UDU5VxC8EpXCYwxs1Fm1g2YDjzU8iMWEZHmUvCn/iKzo2tU+4AvuvsmM7sO+KWZdQF2E11fMrMRwJ3uPtXda83sS8ATQGfgLnd/OU/voS3SqdD9tC/2077YT/tivxbbF+ae8ZKNiIhIQWgTp/5ERKTjUqASEZGCpkDVzpnZXWZWY2bLUtoGmtlTZvZa9HtAyms3RuWmVprZh1LaJ5vZ0ui1/zIzi9q7m1l51P4PMzu6Nd/foTCzmJnNNbNXonJcN0TtHW5/mFkPM3vJzBZH++LWqL3D7YskM+tsZgvN7JHoeYfcF2a2OnoPi5Ip53nfF+6un3b8A7wfmAQsS2n7MfDN6PE3gR9Fj8cR7lPrDowC/gl0jl57CTidcG/aX4APR+1fAO6IHk8HyvP9nnPsi+HApOhxX+DV6D13uP0RjbtP9Lgr8A/gtI64L1L2ydeAe4BHoucdcl8Aq4HBaW153Rd53yn6afkf4GgODFQrgeHR4+HAyujxjcCNKf2eiL5ow4EVKe0fA36T2id63IVwl77l+z03cr88CFzQ0fcH0AtYQKja0iH3BeEeyzmEW12Sgaqj7ovVHByo8rovdOqvYxrm7usBot9Do/ZMJadGRj/VGdoPWMZDLcYtwKAWG3kziU43TCQcSXTI/RGd6loE1ABPuXuH3RfAL4BvAPUpbR11XzjwpJnNt1BaDvK8L9rKfVTSOrKVnMpViqrRZaoKhZn1AWYDX3X3rdGp84xdM7S1m/3h7nXABDM7ArjfwkwF2bTbfWFmFwE17j7fzM5uzCIZ2trFvoi8z93XWSi08JSZrcjRt1X2hY6oOqYNZjYcIPpdE7VnKzlVHT1Obz9gGQs3X/cnlLkqSGbWlRCk/uTuyRJcHXZ/ALj7ZuBZ4EI65r54H3CJma0mzLBwrpn9kY65L3D3ddHvGuB+Qnm6vO4LBaqO6SFCrUSi3w+mtE+PsnJGAWOAl6JD/W1mdlqUufOptGWS67oSeMajk8+FJhr774BX3P3nKS91uP1hZkOiIynMrCehoPMKOuC+cPcb3b3I3Y8mXNx/xt0/QQfcF2bW28z6Jh8DHyTMYJHffZHvC3f6adkfQl3E9YTyU9XAtYTzwXMI9RHnAANT+t9EyNxZSZSlE7XHoy/sP4Ffs7+qSQ9CUeBVhCyf0fl+zzn2xZmEUwxLgEXRz9SOuD+AU4CF0b5YBnw7au9w+yJtv5zN/mSKDrcvgNGELL7FwMvATYWwL1RCSURECppO/YmISEFToBIRkYKmQCUiIgVNgUpERAqaApWIiBQ0Bap2wsy2t+C6V5vZ4EPoP8HMprbUeFqDmf172vO/t/L2zzazM1pzm/liZj3NbJ6ZdW7i8r83sytTnn/MzG5qvhEe0lg+bWa/buKyXzKzzzT3mNoDBSppCRMI9ycVrEb8UTwgULl7sweN6K78bM4GWjVQNSVQNPAeGusa4D4PJZ2aw4XA44ezgqYGzcN0F/CVPGy34ClQtWNR9YHZZlYZ/bwvap9iZn+3MPfO381sbNTe2cx+Gs0hs8TMvpyyui+b2YLoteOzrcfMugHfAcoszGdTFt3tflc0hoVm9tEMY+1kZrdbmBvpETN7LPlfsoV5beZFRTKfSCnl8qyZ/cjCvEqvmtlZKe/jJ9H2lpjZ/4naz7YwH9U9wNKo7YFovS9bVIDTzG4Dekbj/1PUtj36bdG6l0X7oixl3c+a2b1mtsLM/hTdkZ/+Pp81sx+Y2TzgBjO72MKcPAvN7GkzG2ahYO7ngH+JxnBWts8ybd2fNrMHzexxC3MD3Zzy2iei/bTIzH6T/ENsZtvN7Dtm9g9C1evU9V0XbWtxtO1eUfvvzeznZjYX+FG2z9fMjjazv0bfmwWW/QjxKqKqBQ18D74dbWOZmc3Isn+N8I/Sgkzfz5Tvx0HfcwtnDr5tZs8D08zsmGhfzo/ex/GZBm9mn4m+f/MI5ZiS7Zk+204W5nQakvJ+V5nZYHffCaw2sylZ9lPHle87ofXTbHeUb8/Qdg9wZvS4mFA6CKAf0CV6fD4wO3r8eUIdvORrA6Pfq4EvR4+/ANzZwHo+Dfw6ZRw/AD4RPT6CMA9U77SxXgk8Rvjn6UhgU9TWFfg7MCTqVwbcFT1+FvhZ9Hgq8HT0+HrgP6LH3YEqwlw5ZwM7gFEp202+x56Eu+gHZdqfyefAFcBTQGdgGLCGMKXB2YQq0EXRe3ghue/T1vMscHvK8wHsv2P/synv5xbg6w19lmnr/jShCsmglPcTB04AHga6Rv1uBz4VPXagNMt3alDK4++lfAd+DzzC/nmHMn6+hOlDekTtY4CqDNvoBrzV0Pcg9bOKHv8BuDhlPMk+k4D/PYzv+TdStjEHGBM9PpVQ6id9/MOj78CQ6L38jei7n+OzvZlQEBlCiaLZKeu7Cfi/+f57Umg/qp7evp0PjEv5x7OfhTpe/YG7zWwM4Q9V15T+d3govY+7pxaKTBZwnQ9cHj3Otp50HyQU/fx69LwH0R/blD5nArPcvR54K/pvHWAscBKhijOEALE+y7iOTtneKbb/ukV/wh/KvYQ6ZG+kLP8VM7ssehyL+r2T5X0kx/lnD6epNkT/RZcAW6N1VwNYmD7jaOD5DOsoT3lcBJRbOErsBryRoT9k+SzdfVtav6fc/Z1oDPdF460FJgOV0fI92V9UtI7wRzuTk8zse4Tg04cwj1DSLN9/qi7b57sO+LWZTYi2c1yGbQwGNqc8z/Y9ADjHzL5BCIADCSV+Hk5b34WESfqgad/zcnivwv4ZwKyUfd49w/hPBZ51943RcuUp7zPbZ3sX4QjyF4TTnv+Tsr4aIOORW0emQNW+dSJMULYrtdHMfgXMdffLotNMzyZfInu5/T3R7zr2f2++m2U96Qy4wt1X5hhrtrk2DHjZ3U/P8nqmcRnhv//UP6xYmMJhR9rz8wn7aKeZPUv4I5tL1jlBUsaSPp50O1Ie/wr4ubs/FI3nlizLZPwsM0j//JJTLtzt7jdm6L/bs18b+j1wqbsvNrNPE44ak1LfQ8bP18xuATYA46Px786wjV0cuM8z7l8z60E4Eoy7eyJad6bP6oOEo17I/v3M9T1Pvq9OwGZ3n5A2js6Ef4ogFFddkGNdGT/baPwbzOxcQqC7KmWZHoR9Iil0jap9exL4UvJJ9J8thP8010aPP53W/3MWXSA3s4ENrD/berYRpnpPeoJwjcui9U7MsK7ngSuic/bD2P9HcSUwxMxOj5btamYnNjCuJ4DPW5jSAzM7zkIl6Ezj3xQFqeMJU7En7Usun+Y5wvW3ztF1hvcTCms2Veo+vDqlPX0fZvss011gZgMtVES/lHAqag5wpYX5hYheP6oRY+sLrI/2w1U5+mX7fPsD66Ojo08SjoYP4O6bgM5RIILs34Pk629HRztXksbM+hNO5yWPiJv8PXf3rcAbZjYt6mNmNt7d69x9QvTzbcLEm2eb2aBoP01LWU22zxbgTuCPQEXaPwrHEU7ZSgoFqvajl5lVp/x8jZBBFI8uGC8nXKAH+DHwQzP7Gwf+8biTcL59iZktBj7ewDazrWcu4TTVIgvJBt8lnHZZYmbLoufpZhOquy8DfkP4A7DF3fcS/ij9KBrTIhrOhrsTWE64oJ5cX6ajm8eBLma2JBrTiymvzYjG+6e0Ze4nVBxfDDxDuKbxVgPjyeUWwumlvxKm5E56GLgs2odnkf2zTPc84frNIsK1jyp3Xw78B2HW1iWEa2zDGzG2bxE+h6cIU4Bkk+3zvR242sxeJPwB3pFl+ScJp/wg+/dgM/BbQhLMA0BlhvVcADyd8vxwv+dXAddGfV4GDkoC8jCdxS2Ea5JPE46wkm4h82cL4WisDwee9oOQjPE0cgBVT5eCYWZ93H27mQ0iHKW87zCDQIcSnZ6Lu/uXGupbSKIjsK+5+yej5036HpjZnYREnxcb6ptvZhYH/tPdz0ppO2A/yH66RiWF5BELk/l1A76rINUxuPtCC7cNdI5OgzXpe+Dun23RgTYTM/smIfMw/XTqYMJRrKTREZWIiBQ0XaMSEZGCpkAlIiIFTYFKREQKmgKViIgUNAUqEREpaP8fvx9d5apBT1MAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "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.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
