{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Composting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import PFAS_SAT_ProcessModels as pspd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from IPython.display import Image\n",
    "import pandas as pd\n",
    "pd.set_option('display.max_colwidth', 0)\n",
    "import warnings\n",
    "warnings.filterwarnings('ignore')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Model document\n",
    "\n",
    "Composting can accept specific fractions of MSW including food, yard, and some paper wastes. It is also capable of treating dried and stabilized WWT solids. The primary purpose of composting is to produce a stabilized organic product that can be used as a soil amendment to enhance water retention, carbon and nutrient content, and erosion control among other potential benefits. In addition to finished compost, composting systems also produce residual solid materials from pre- and post-screening. Those materials can be combusted or landfilled. Compost also may produce leachate that needs to be managed through wastewater treatment. It is also possible that PFAS are volatilized during aerobic composting due to the elevated temperatures and stripping from active aeration, but there is currently no data on this potential pathway.\n",
    "\n",
    "In the composting process model, feedstocks (e.g., WWT solids or food waste) are mixed with amendments (wood chips by default) to increase porosity to allow air flow through the system. Active composting may be under cover or open to the atmosphere. In the latter case, there will be contact water that may contain PFAS that will need to be managed.\n",
    "\n",
    "Model predictions are based on achievement of equilibrium. 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 PFAS leaches through the compost piles (i.e., it is assumed that annual precipitation is uniform throughout the year and continuously removes PFAS from the mixture). The model then calculates a fraction of precipitation that becomes contact water. The contact water can either be collected and managed or be released to surface or ground- water depending on the composting process and applicable regulations. If active composting is under a cover, then no contact water is included. The PFAS that is not removed in the contact water remains in the compost. By default, the compost is cured after active composting. The curing model is the same as the active composting model, except it is not covered and contact water is not managed. Any PFAS that does not leave the curing piles through leaching or run-off remains in the finished compost that will be either land applied or landfilled. \n",
    "\n",
    "The behavior of the finished compost in a landfill and land application process are described in their respective sections. There are several different compost processes (e.g., windrows, static piles, in-vessel). The compost process model in the SAT framework is designed so that by changing default parameters, any compost process can be represented. \n",
    "\n",
    "<img src=\"../Images/ProcessModels/Compost_Diagram.png\" alt=\"Drawing\" style=\"width: 700px;\"/>\n",
    "\n",
    "\n",
    "### Assumptions and Limitations\n",
    "1.\tOrganic carbon partitioning coefficients determined from soils/sediments are used for compost partitioning. \n",
    "2.\tWood chips are currently the only amendment that is built into the material flow properties. \n",
    "3.\tWe assume that the mass of solid loss per carbon loss is similar both for active and curing composting stages. \n",
    "4.\tVolatilization is assumed to be zero by default. However, the user may assign a fraction of the PFAS to be volatilized.\n",
    "5.\tFuture work is required to implement a dynamic (i.e., non-equilibrium) model to account for changes in the organic C content over time as composted materials decompose, and to account for episodic precipitation events. \n",
    "6.\tWe do not consider the possible conversion of \"precursor\" compounds to commonly measured PFAAs."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for Composting model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Category</th>\n",
       "      <th>Dictonary_Name</th>\n",
       "      <th>Parameter Name</th>\n",
       "      <th>Parameter Description</th>\n",
       "      <th>amount</th>\n",
       "      <th>unit</th>\n",
       "      <th>Reference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA Log Koc (Composting)</td>\n",
       "      <td>2.19</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS Log Koc (Composting)</td>\n",
       "      <td>3.04</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA Log Koc (Composting)</td>\n",
       "      <td>1.88</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA Log Koc (Composting)</td>\n",
       "      <td>1.37</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA Log Koc (Composting)</td>\n",
       "      <td>1.77</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA Log Koc (Composting)</td>\n",
       "      <td>1.97</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA Log Koc (Composting)</td>\n",
       "      <td>2.63</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA Log Koc (Composting)</td>\n",
       "      <td>3.24</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS Log Koc (Composting)</td>\n",
       "      <td>1.51</td>\n",
       "      <td>log L/kg OC</td>\n",
       "      <td>[2,3,4,5]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Log partition coefficient</td>\n",
       "      <td>LogPartCoef</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS Log Koc (Composting)</td>\n",
       "      <td>2.79</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>Amendment Material Properties</td>\n",
       "      <td>AmndProp</td>\n",
       "      <td>mass_ratio</td>\n",
       "      <td>Mass of amendment to mass of feedstock</td>\n",
       "      <td>0.70</td>\n",
       "      <td>kg TS/kg TS</td>\n",
       "      <td>[62]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Amendment Material Properties</td>\n",
       "      <td>AmndProp</td>\n",
       "      <td>ts_cont</td>\n",
       "      <td>TS content of amendmente _wet</td>\n",
       "      <td>0.55</td>\n",
       "      <td>fraction wet weight</td>\n",
       "      <td>[63]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Amendment Material Properties</td>\n",
       "      <td>AmndProp</td>\n",
       "      <td>C_cont</td>\n",
       "      <td>Organic C content - dry</td>\n",
       "      <td>0.58</td>\n",
       "      <td>fraction TS</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>PFOA concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>PFOS concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>PFBA concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>PFPeA concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>PFHxA concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>PFHpA concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>PFNA concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>PFDA concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>PFBS concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>Amendment PFAS Concentration</td>\n",
       "      <td>AmndPFAS</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>PFHxS concentration - dry mass</td>\n",
       "      <td>0.00</td>\n",
       "      <td>10e-6g/kg</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>is_covered</td>\n",
       "      <td>Is active composting (AC) covered?</td>\n",
       "      <td>0.00</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>bulk_dens</td>\n",
       "      <td>Wet bulk density (AC pile)</td>\n",
       "      <td>530.00</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>wind_ht</td>\n",
       "      <td>Windrow height (AC pile)</td>\n",
       "      <td>2.50</td>\n",
       "      <td>m</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>wind_wid</td>\n",
       "      <td>Windrow width (AC pile)</td>\n",
       "      <td>4.50</td>\n",
       "      <td>m</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>frac_C_lost</td>\n",
       "      <td>Fraction of organic C lost during AC</td>\n",
       "      <td>0.50</td>\n",
       "      <td>kg C loss/kg C</td>\n",
       "      <td>[22]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>ts_end</td>\n",
       "      <td>Total solids content at the end of AC</td>\n",
       "      <td>0.55</td>\n",
       "      <td>kg TS/kg</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>frac_sol_to_cw</td>\n",
       "      <td>Fraction of solids in contact water (AC)</td>\n",
       "      <td>0.00</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>sol_loss_per_C_loss</td>\n",
       "      <td>Total solids loss per loss per C loss (AC)</td>\n",
       "      <td>2.00</td>\n",
       "      <td>kg TS/kg C</td>\n",
       "      <td>[64,24]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>is_cw_col</td>\n",
       "      <td>Is AC contact water collected?</td>\n",
       "      <td>1.00</td>\n",
       "      <td>1:TRUE,0:FALSE</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>32</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>frac_cw_col</td>\n",
       "      <td>Fraction contact water collection during AC</td>\n",
       "      <td>0.98</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>33</th>\n",
       "      <td>Active composting Input</td>\n",
       "      <td>AComp</td>\n",
       "      <td>ac_time</td>\n",
       "      <td>Activate composting time</td>\n",
       "      <td>70.00</td>\n",
       "      <td>days</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>34</th>\n",
       "      <td>Curing Input</td>\n",
       "      <td>Curing</td>\n",
       "      <td>frac_C_lost</td>\n",
       "      <td>Fraction of organic C lost during curing</td>\n",
       "      <td>0.10</td>\n",
       "      <td>kg C/kg C</td>\n",
       "      <td>[42]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>35</th>\n",
       "      <td>Curing Input</td>\n",
       "      <td>Curing</td>\n",
       "      <td>ts_end</td>\n",
       "      <td>Total solids content at the end of curing</td>\n",
       "      <td>0.65</td>\n",
       "      <td>kg TS/kg</td>\n",
       "      <td>[42]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>36</th>\n",
       "      <td>Curing Input</td>\n",
       "      <td>Curing</td>\n",
       "      <td>bulk_dens</td>\n",
       "      <td>Wet bulk density of curing pile</td>\n",
       "      <td>530.00</td>\n",
       "      <td>kg/m3</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>37</th>\n",
       "      <td>Curing Input</td>\n",
       "      <td>Curing</td>\n",
       "      <td>wind_ht</td>\n",
       "      <td>Windrow height of curing pile</td>\n",
       "      <td>2.00</td>\n",
       "      <td>m</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>38</th>\n",
       "      <td>Curing Input</td>\n",
       "      <td>Curing</td>\n",
       "      <td>wind_wid</td>\n",
       "      <td>Windrow width of curing pile</td>\n",
       "      <td>4.00</td>\n",
       "      <td>m</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>39</th>\n",
       "      <td>Curing Input</td>\n",
       "      <td>Curing</td>\n",
       "      <td>sol_loss_per_C_loss</td>\n",
       "      <td>Total solids loss per loss per C loss (Curing)</td>\n",
       "      <td>2.00</td>\n",
       "      <td>kg TS/kg C</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>40</th>\n",
       "      <td>Curing Input</td>\n",
       "      <td>Curing</td>\n",
       "      <td>curing_time</td>\n",
       "      <td>Curing composting time</td>\n",
       "      <td>30.00</td>\n",
       "      <td>days</td>\n",
       "      <td>[48]</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>41</th>\n",
       "      <td>Precipitation Data</td>\n",
       "      <td>Precip</td>\n",
       "      <td>ann_precip</td>\n",
       "      <td>Annual precipitation rate</td>\n",
       "      <td>1.13</td>\n",
       "      <td>m/yr</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>42</th>\n",
       "      <td>Volatilization</td>\n",
       "      <td>Volatilization</td>\n",
       "      <td>frac_vol_loss</td>\n",
       "      <td>Fraction of PFAS lost to volatilization</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>43</th>\n",
       "      <td>Screen</td>\n",
       "      <td>Screen</td>\n",
       "      <td>frac_rmvd</td>\n",
       "      <td>Fraction of incoming feedstock removed by screen</td>\n",
       "      <td>0.00</td>\n",
       "      <td>fraction</td>\n",
       "      <td></td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                         Category  Dictonary_Name       Parameter Name  \\\n",
       "0   Log partition coefficient      LogPartCoef     PFOA                  \n",
       "1   Log partition coefficient      LogPartCoef     PFOS                  \n",
       "2   Log partition coefficient      LogPartCoef     PFBA                  \n",
       "3   Log partition coefficient      LogPartCoef     PFPeA                 \n",
       "4   Log partition coefficient      LogPartCoef     PFHxA                 \n",
       "5   Log partition coefficient      LogPartCoef     PFHpA                 \n",
       "6   Log partition coefficient      LogPartCoef     PFNA                  \n",
       "7   Log partition coefficient      LogPartCoef     PFDA                  \n",
       "8   Log partition coefficient      LogPartCoef     PFBS                  \n",
       "9   Log partition coefficient      LogPartCoef     PFHxS                 \n",
       "10  Amendment Material Properties  AmndProp        mass_ratio            \n",
       "11  Amendment Material Properties  AmndProp        ts_cont               \n",
       "12  Amendment Material Properties  AmndProp        C_cont                \n",
       "13  Amendment PFAS Concentration   AmndPFAS        PFOA                  \n",
       "14  Amendment PFAS Concentration   AmndPFAS        PFOS                  \n",
       "15  Amendment PFAS Concentration   AmndPFAS        PFBA                  \n",
       "16  Amendment PFAS Concentration   AmndPFAS        PFPeA                 \n",
       "17  Amendment PFAS Concentration   AmndPFAS        PFHxA                 \n",
       "18  Amendment PFAS Concentration   AmndPFAS        PFHpA                 \n",
       "19  Amendment PFAS Concentration   AmndPFAS        PFNA                  \n",
       "20  Amendment PFAS Concentration   AmndPFAS        PFDA                  \n",
       "21  Amendment PFAS Concentration   AmndPFAS        PFBS                  \n",
       "22  Amendment PFAS Concentration   AmndPFAS        PFHxS                 \n",
       "23  Active composting Input        AComp           is_covered            \n",
       "24  Active composting Input        AComp           bulk_dens             \n",
       "25  Active composting Input        AComp           wind_ht               \n",
       "26  Active composting Input        AComp           wind_wid              \n",
       "27  Active composting Input        AComp           frac_C_lost           \n",
       "28  Active composting Input        AComp           ts_end                \n",
       "29  Active composting Input        AComp           frac_sol_to_cw        \n",
       "30  Active composting Input        AComp           sol_loss_per_C_loss   \n",
       "31  Active composting Input        AComp           is_cw_col             \n",
       "32  Active composting Input        AComp           frac_cw_col           \n",
       "33  Active composting Input        AComp           ac_time               \n",
       "34  Curing Input                   Curing          frac_C_lost           \n",
       "35  Curing Input                   Curing          ts_end                \n",
       "36  Curing Input                   Curing          bulk_dens             \n",
       "37  Curing Input                   Curing          wind_ht               \n",
       "38  Curing Input                   Curing          wind_wid              \n",
       "39  Curing Input                   Curing          sol_loss_per_C_loss   \n",
       "40  Curing Input                   Curing          curing_time           \n",
       "41  Precipitation Data             Precip          ann_precip            \n",
       "42  Volatilization                 Volatilization  frac_vol_loss         \n",
       "43  Screen                         Screen          frac_rmvd             \n",
       "\n",
       "                               Parameter Description  amount  \\\n",
       "0   PFOA Log Koc (Composting)                         2.19     \n",
       "1   PFOS Log Koc (Composting)                         3.04     \n",
       "2   PFBA Log Koc (Composting)                         1.88     \n",
       "3   PFPeA Log Koc (Composting)                        1.37     \n",
       "4   PFHxA Log Koc (Composting)                        1.77     \n",
       "5   PFHpA Log Koc (Composting)                        1.97     \n",
       "6   PFNA Log Koc (Composting)                         2.63     \n",
       "7   PFDA Log Koc (Composting)                         3.24     \n",
       "8   PFBS Log Koc (Composting)                         1.51     \n",
       "9   PFHxS Log Koc (Composting)                        2.79     \n",
       "10  Mass of amendment to mass of feedstock            0.70     \n",
       "11  TS content of amendmente _wet                     0.55     \n",
       "12  Organic C content - dry                           0.58     \n",
       "13  PFOA concentration - dry mass                     0.00     \n",
       "14  PFOS concentration - dry mass                     0.00     \n",
       "15  PFBA concentration - dry mass                     0.00     \n",
       "16  PFPeA concentration - dry mass                    0.00     \n",
       "17  PFHxA concentration - dry mass                    0.00     \n",
       "18  PFHpA concentration - dry mass                    0.00     \n",
       "19  PFNA concentration - dry mass                     0.00     \n",
       "20  PFDA concentration - dry mass                     0.00     \n",
       "21  PFBS concentration - dry mass                     0.00     \n",
       "22  PFHxS concentration - dry mass                    0.00     \n",
       "23  Is active composting (AC) covered?                0.00     \n",
       "24  Wet bulk density (AC pile)                        530.00   \n",
       "25  Windrow height (AC pile)                          2.50     \n",
       "26  Windrow width (AC pile)                           4.50     \n",
       "27  Fraction of organic C lost during AC              0.50     \n",
       "28  Total solids content at the end of AC             0.55     \n",
       "29  Fraction of solids in contact water (AC)          0.00     \n",
       "30  Total solids loss per loss per C loss (AC)        2.00     \n",
       "31  Is AC contact water collected?                    1.00     \n",
       "32  Fraction contact water collection during AC       0.98     \n",
       "33  Activate composting time                          70.00    \n",
       "34  Fraction of organic C lost during curing          0.10     \n",
       "35  Total solids content at the end of curing         0.65     \n",
       "36  Wet bulk density of curing pile                   530.00   \n",
       "37  Windrow height of curing pile                     2.00     \n",
       "38  Windrow width of curing pile                      4.00     \n",
       "39  Total solids loss per loss per C loss (Curing)    2.00     \n",
       "40  Curing composting time                            30.00    \n",
       "41  Annual precipitation rate                         1.13     \n",
       "42  Fraction of PFAS lost to volatilization           0.00     \n",
       "43  Fraction of incoming feedstock removed by screen  0.00     \n",
       "\n",
       "                   unit  Reference  \n",
       "0   log L/kg OC          [2,3,4,5]  \n",
       "1   log L/kg OC          [2,3,4,5]  \n",
       "2   log L/kg OC          [2,3,4,5]  \n",
       "3   log L/kg OC          [2,3,4,5]  \n",
       "4   log L/kg OC          [2,3,4,5]  \n",
       "5   log L/kg OC          [2,3,4,5]  \n",
       "6   log L/kg OC          [2,3,4,5]  \n",
       "7   log L/kg OC          [2,3,4,5]  \n",
       "8   log L/kg OC          [2,3,4,5]  \n",
       "9   log L/kg OC          [2,3,4,5]  \n",
       "10  kg TS/kg TS          [62]       \n",
       "11  fraction wet weight  [63]       \n",
       "12  fraction TS          [48]       \n",
       "13  10e-6g/kg                       \n",
       "14  10e-6g/kg                       \n",
       "15  10e-6g/kg                       \n",
       "16  10e-6g/kg                       \n",
       "17  10e-6g/kg                       \n",
       "18  10e-6g/kg                       \n",
       "19  10e-6g/kg                       \n",
       "20  10e-6g/kg                       \n",
       "21  10e-6g/kg                       \n",
       "22  10e-6g/kg                       \n",
       "23  1:TRUE,0:FALSE                  \n",
       "24  kg/m3                [48]       \n",
       "25  m                    [48]       \n",
       "26  m                    [48]       \n",
       "27  kg C loss/kg C       [22]       \n",
       "28  kg TS/kg             [48]       \n",
       "29  fraction                        \n",
       "30  kg TS/kg C           [64,24]    \n",
       "31  1:TRUE,0:FALSE                  \n",
       "32  fraction                        \n",
       "33  days                 [48]       \n",
       "34  kg C/kg C            [42]       \n",
       "35  kg TS/kg             [42]       \n",
       "36  kg/m3                [48]       \n",
       "37  m                    [48]       \n",
       "38  m                    [48]       \n",
       "39  kg TS/kg C                      \n",
       "40  days                 [48]       \n",
       "41  m/yr                            \n",
       "42  NaN                             \n",
       "43  fraction                        "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Composting = pspd.Comp()\n",
    "Composting.InputData.Data[['Category','Dictonary_Name','Parameter Name', 'Parameter Description', 'amount', 'unit','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming Dewatered WWT Solids to Composting"
   ]
  },
  {
   "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>200.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>800.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>VS flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>120.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>76.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>934.565</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>3678.926</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>146.226</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>365.103</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>511.518</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>88.188</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>373.149</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1552.18</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>125.018</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>107.265</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   200.0   \n",
       "2   Moisture flow  kg   800.0   \n",
       "3   VS flow        kg   120.0   \n",
       "4   Carbon flow    kg   76.0    \n",
       "5   PFOA           μg   934.565 \n",
       "6   PFOS           μg   3678.926\n",
       "7   PFBA           μg   146.226 \n",
       "8   PFPeA          μg   365.103 \n",
       "9   PFHxA          μg   511.518 \n",
       "10  PFHpA          μg   88.188  \n",
       "11  PFNA           μg   373.149 \n",
       "12  PFDA           μg   1552.18 \n",
       "13  PFBS           μg   125.018 \n",
       "14  PFHxS          μg   107.265 "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncominWaste = pspd.IncomFlow()\n",
    "IncominWaste.set_flow('Dewatered WWT Solids', 1000)\n",
    "IncominWaste.calc()\n",
    "DewateredWWTSolids = IncominWaste.Inc_flow\n",
    "DewateredWWTSolids.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in Composting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Volatilized</th>\n",
       "      <th>Compost</th>\n",
       "      <th>Contact water</th>\n",
       "      <th>Collected Contact water</th>\n",
       "      <th>Compost Residuals</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>PFOA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>97.92</td>\n",
       "      <td>0.50</td>\n",
       "      <td>1.58</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>99.70</td>\n",
       "      <td>0.07</td>\n",
       "      <td>0.23</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>95.94</td>\n",
       "      <td>0.98</td>\n",
       "      <td>3.08</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>88.83</td>\n",
       "      <td>2.81</td>\n",
       "      <td>8.36</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>94.88</td>\n",
       "      <td>1.24</td>\n",
       "      <td>3.88</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>96.65</td>\n",
       "      <td>0.81</td>\n",
       "      <td>2.55</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>99.23</td>\n",
       "      <td>0.18</td>\n",
       "      <td>0.59</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>0.0</td>\n",
       "      <td>99.81</td>\n",
       "      <td>0.05</td>\n",
       "      <td>0.15</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>91.39</td>\n",
       "      <td>2.13</td>\n",
       "      <td>6.48</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>0.0</td>\n",
       "      <td>99.46</td>\n",
       "      <td>0.13</td>\n",
       "      <td>0.41</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Volatilized  Compost  Contact water  Collected Contact water  \\\n",
       "PFOA   0.0          97.92    0.50           1.58                      \n",
       "PFOS   0.0          99.70    0.07           0.23                      \n",
       "PFBA   0.0          95.94    0.98           3.08                      \n",
       "PFPeA  0.0          88.83    2.81           8.36                      \n",
       "PFHxA  0.0          94.88    1.24           3.88                      \n",
       "PFHpA  0.0          96.65    0.81           2.55                      \n",
       "PFNA   0.0          99.23    0.18           0.59                      \n",
       "PFDA   0.0          99.81    0.05           0.15                      \n",
       "PFBS   0.0          91.39    2.13           6.48                      \n",
       "PFHxS  0.0          99.46    0.13           0.41                      \n",
       "\n",
       "       Compost Residuals  \n",
       "PFOA   0.0                \n",
       "PFOS   0.0                \n",
       "PFBA   0.0                \n",
       "PFPeA  0.0                \n",
       "PFHxA  0.0                \n",
       "PFHpA  0.0                \n",
       "PFNA   0.0                \n",
       "PFDA   0.0                \n",
       "PFBS   0.0                \n",
       "PFHxS  0.0                "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Composting.calc(Inc_flow=DewateredWWTSolids)\n",
    "Composting.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAx0ElEQVR4nO3deVxVdf7H8ddlR1AQNVQwEUlwLtzLrqiokFvuNS4g5ZY1LpOtmk6TVpOlxZi5TGbLWL+fC2o4OpnLqIhLOiwKhIgiiguiASUCynY5vz/8zX3EACUKHpbP8/Hg8fCe5Xs+53Z5czr3fL9fjaIoCCGEePhM1C5ACCFaKglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVmNVl4/bt2ysuLi4NVIoQQjRPiYmJeYqidPjv5XUKYBcXFxISEuqvKiHq4Pr167z00kvEx8djb2+Po6MjK1asoEePHqrWdfPmTTZu3Mjs2bNVrUM0XhqN5lJNy+UWhGgSFEXhySefZODAgWRmZpKYmMj777/PjRs31C6Nmzdv8re//U3tMkQTJAEsmoSYmBjMzc2ZOXOmcZler6dfv37MmzcPT09PvLy8iIqKAuDQoUMMGDCAMWPG4OrqyoIFC9iwYQOBgYF4eXmRmZkJwNSpU5k5cyb+/v706NGDb7/9FoCSkhKmTZuGl5cXPj4+xMTEAHD69GkCAwPx9vZGp9ORkZHBggULyMzMxNvbm3nz5j3kd0Y0ZXW6BSGEWlJTU/Hz86u2PDo6mqSkJJKTk8nLyyMgIID+/fsDkJyczJkzZ3BwcMDV1ZUZM2YQFxfHxx9/zKpVq1ixYgUAWVlZxMXFkZmZSUhICOfPn2fNmjVoNBp++OEH0tPTGTJkCOfOnWPt2rW8+OKLREREUFZWhsFgYOnSpaSmppKUlPQQ3xHRHMgVsGjSjh49Snh4OKampjg6OjJgwADi4+MBCAgIoFOnTlhaWtK9e3eGDBkCgJeXF1lZWcY2JkyYgImJCY899hiurq6kp6dz9OhRnn76aQA8PDzo2rUr586dIygoiPfee49ly5Zx6dIlrK2tH/o5i+ZDAlg0CVqtlsTExDrtY2lpafy3iYmJ8bWJiQkVFRXGdRqNpsp+//36lyZNmsTOnTuxtrZm+PDhHDx4sE41CfFLEsCiSQgNDaW0tJR169YZl6WkpGBvb09UVBQGg4Hc3FwOHz5MYGBgndreunUrlZWVZGZmcuHCBdzd3QkODmbDhg0AnDt3jsuXL+Pu7s6FCxdwdXVl7ty5jBkzhpSUFFq3bk1hYWG9nq9oGSSARZOg0WjYvn07+/fvp3v37mi1WhYuXMikSZPQ6XTo9XpCQ0P54IMP6NixY53afvTRRwkMDOSJJ55g7dq1WFlZMXv2bCorK/Hy8mLixImsX78eS0tLtmzZgqenJ97e3qSmpjJ58mTatWtH37598fT0lC/hRJ1o6jIesL+/vyLPAYvmZOrUqYwcOZJx48apXYpoxjQaTaKiKP7/vVyugIUQQiXyGJpo0davX692CaIFkytg0eiZmpri7e2Np6cno0aN4ubNm3VuIyEhgblz59a4zsXFhby8vPuq7a233iIyMvK+9hVCAlg0etbW1iQlJZGamoqDgwNr1qypcxv+/v6sXLmyAaoT4v5JAIsmJSgoiOzsbAAyMzMZNmwYfn5+BAcHk56eDtx9rMzT0xO9Xm/sFXfo0CFGjhwJQH5+PkOGDEGr1TJjxgz+80V0VlYWnp6exmNFRkby1ltvAfDZZ58REBCAXq/n97//Pbdv365W28qVK/nd736HTqcjLCyswd4D0XxIAIsmw2AwcODAAUaPHg3A888/z6pVq0hMTCQyMtI4Gtk777zD3r17SU5OZufOndXaefvtt+nXrx+nT5/mySef5PLly7957Keeeor4+HiSk5Pp2bMnX3zxRbVtli5dyqlTp0hJSWHt2rUPeLaiJZAv4USjd+fOHby9vcnOzqZnz54MHjyYoqIivv/+e8aPH2/crrS0FIC+ffsydepUJkyYwFNPPVWtvcOHDxMdHQ3AiBEjaNu27W/WkJqayp///Gdu3rxJUVERQ4cOrbaNTqcjIiKCsWPHMnbs2Ps8W9GSyBWwaPT+cw/40qVLKIrCmjVrqKysxN7enqSkJOPPmTNnAFi7di3vvvsuV65cwc/Pj/z8/Hs6jpmZGZWVlcbXJSUlxn9PnTqV1atX88MPP7B48eIq6/5j165dzJkzh5MnTxIQEFClu7MQNZEAFk1Gq1atWLlyJX/9619p1aoV3bp1Y+vWrcDd8YKTk5OBu/eGe/XqxTvvvEOHDh24cuVKlXb69+/Pxo0bAdi9ezc///wzAI6Ojvz444/k5+dTWlpqHJoSoLCwkE6dOlFeXm7sovxLlZWVXLlyhZCQEJYtW0ZBQQFFRUUN8j6I5kNuQYgmxcfHB51Ox6ZNm9iwYQOzZs3i3Xffpby8nLCwMPR6PfPmzSMjIwNFUXj88cfR6/XExsYa21i8eDHh4eFotVr69OnDo48+CoC5uTmLFi0iMDAQJycnPDw8jPv85S9/oVevXnTo0IFevXpVG/vBYDDw9NNPU1BQgKIozJ07F3t7+4fynoimS7oiCyFEA5OuyEII0chIAAshhEokgIUQQiUSwMLo+vXrhIWF0b17d/z8/Bg+fDjnzp2rczsrVqyosafYvUhKSuK77767r33rYs+ePbi7u+Pm5sbSpUsb/HhC1EQCWAD1O+17Yw9gg8HAnDlz2L17N2lpaWzatIm0tLQGPaYQNZEAFsD9Tfs+cOBAxo0bh4eHBxERESiKwsqVK7l27RohISGEhIQAMGvWLPz9/dFqtSxevNjYfnx8PH369EGv1xMYGEhBQQGLFi0iKioKb29v47HqW1xcHG5ubri6umJhYUFYWBg7duxokGMJ8WvkOWAB3N+076dOneL06dN07tyZvn37cuzYMebOncvy5cuJiYmhffv2ACxZsgQHBwcMBgOPP/44KSkpeHh4MHHiRKKioggICODWrVu0atWKd955h4SEBFavXt1g55qdnU2XLl2Mr52dnfn3v//dYMcTojZyBSx+1a9N+x4YGIizszMmJiZ4e3tXmer9l7Zs2YKvry8+Pj6cPn2atLQ0zp49S6dOnQgICACgTZs2mJnJ9YBoWSSABfDg076bmprWOPbBxYsXiYyM5MCBA6SkpDBixIgax1F4mJycnKp0T7569SpOTk4qViRaKglgAdTvtO+/nKb91q1b2NjYYGdnx40bN9i9ezcA7u7u5OTkGK+mCwsLqaioeChTvAcEBJCRkcHFixcpKytj8+bNxiEuhXiYJIAFUL/Tvj///PMMGzaMkJAQ9Ho9Pj4+eHh4MGnSJPr27QuAhYUFUVFRvPDCC+j1egYPHkxJSQkhISGkpaU16JdwZmZmrF69mqFDh9KzZ08mTJiAVqttkGMJ8WtkLAghhGhgMhaEEEI0MhLAQgihEglgIYRQiQSwqHcP0hVZiJZEAljUu/sJYIPB0EDV1C4nJ4e6fAktRH2TABa1+vDDD1m5ciUAL7/8MqGhoQAcPHiQiIiIGsd4qGksiH379hEUFISvry/jx483zpXm4uLC66+/jq+vr3Fut4dBURQWvfU2Ts7OhEU8Y5xNWYiHTlGUe/7x8/NTRMtx/PhxZdy4cYqiKEq/fv2UgIAApaysTHnrrbeUtWvXKvn5+YqiKEpFRYUyYMAAJTk5WVEURenatauSm5urKIqi5ObmKsHBwUpRUZGiKIqydOlS5e233zZut2zZsod9WkpiYqICGH/+vn79Q69BtCxAglJDpkrne1ErPz8/EhMTuXXrFpaWlvj6+pKQkMCRI0dYuXIlW7ZsYd26dVRUVJCTk0NaWho6na5KGydOnCAtLc3YAaOsrIygoCDj+okTJz7Uc4K7E3vm5OTQqVMnDhw4QHBw8EOvQQiQ0dDErzA3N6dbt26sX7+ePn36oNPpiImJ4fz581hbWxMZGUl8fDxt27Zl6tSpNY7xoCgKgwcPZtOmTTUew8bGpqFPoxqNRmPszffII49gbm7+0GsQAuQesPgNwcHBREZG0r9/f4KDg1m7di0+Pj61jvEAVceC6N27N8eOHeP8+fMAFBcX39csG0I0RxLA4lcFBweTk5NDUFAQjo6OWFlZERwcXOsYD1B1LIgOHTqwfv16wsPD0el0BAUFkZ6eruIZCdF4yFgQosXSaDT88MMPeHp6ql2KaOZkLAghhGhkJICFEEIlEsBCCKESCWBRzZ49e3B3d8fNzY2lS5dWW19aWsrEiRNxc3OjV69eVeaCe//993Fzc8Pd3Z29e/cCUFJSQmBgIHq9vtrMyAcOHMDX1xdvb2/69etnfFpCiBahpt4Ztf1IT7jmr6KiQnF1dVUyMzOV0tJSRafTKadPn66yzZo1a5Q//OEPiqIoyqZNm5QJEyYoiqIop0+fVnQ6nVJSUqJcuHBBcXV1VSoqKpTKykqlsLBQURRFKSsrUwIDA5Xjx48riqIojz32mJKWlmZsd8qUKQ/pTBUFUH744YeHdjzRclFLTzi5AhZVxMXF4ebmhqurKxYWFoSFhbFjx44q2+zYsYMpU6YAMG7cOA4cOICiKOzYsYOwsDAsLS3p1q0bbm5uxMXFodFosLW1BaC8vJzy8nI0Gg1w90mEW7duAVBQUEDnzp0f4tkKoS7pCSeqyM7OpkuXLsbXzs7O/Pvf/651GzMzM+zs7MjPzyc7O5vevXtX2Tc7Oxu4O9qZn58f58+fZ86cOfTq1QuAzz//nOHDh2NtbU2bNm04ceJEQ5+iEI2GXAGLh8LU1JSkpCSuXr1KXFwcqampAHz00Ud89913XL16lWnTpvHKK6+oXKkQD48EsKjCycmJK1euGF9fvXoVJyenWrepqKigoKCAdu3a3dO+9vb2hISEsGfPHnJzc0lOTjZeDU+cOJHvv/++oU5NiEZHAlhUERAQQEZGBhcvXqSsrIzNmzczevToKtuMHj2ar776CoBt27YRGhqKRqNh9OjRbN68mdLSUi5evEhGRgaBgYHk5uZy8+ZNAO7cucO//vUvPDw8aNu2LQUFBcaxIf71r3/Rs2fPh3q+QqhJ7gGLKszMzFi9ejVDhw7FYDAwffp0tFotixYtwt/fn9GjR/Pss8/yzDPP4ObmhoODA5s3bwZAq9UyYcIEfve732FmZsaaNWswNTUlJyeHKVOmYDAYqKysZMKECYwcORKAzz77jN///veYmJjQtm1bvvzySzVPX4iHSsaCEC2WjAUhHhYZC0IIIRoZCWAhhFCJBLAQQqhEAlgIIVQiASzq1dat2/jnP79VuwwhmgQJYFEvFEXhpVdeZdIzU3hq/HgWv/2O2iUJ0ejJc8Ci3iSlpFJRehuAhFPJKlcjROMnV8CiXmg0Gg7tvzv+r23r1uz6xzcqVyRE4ycBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJIBFvTM1NVW7BCGaBAlgUe+6deumdgm/qby8HABHR0eVKxEtmQSwqDeXLl0CYMVHH6lcyW974YUXAOjQoYPKlYiWTAJY1AtFUXBxcQFgwIAB6hbzG0pLS/n0008JCQlRuxTRwkkAiweWnZ2Nubk5AN9//73K1fy669evY2VlBcCuXbvqrd0ZM2ezd+/eemtPtAwaRVHueWN/f38lISGhAcsRN2/e5MVX51Vb/vWXn6tQTd1MmzaNPn36qF1GjZ577rkqr5ctW4aDg0O9tf/Sq6+hVCokxv8bDw+PemtXNA8ajSZRURT/asslgBuXzMxMPHV6TF2DsOzcw7j8ZuxXVJYWq1hZ7czbtMf+MX9Ao3YptcpN3A2AtaMLts496719xdIGpU1HLNO+JeVkAu3atav3Y4imSwK4Cdm1axcTn5mG/cRlmNk9onY5og6KDq+nm5LDsdiDWFhYqF2OaCRqC2C5B9wIjRgxgrfeWEDhznep/P9JLkXTYBM8mQu3FKbNeJ66XNyIlkkCuJF69ZWXGTtkIEV7/opSaVC7HHGPNBoTbIa+zK6Dx1j6wYdqlyMaOQngRkqj0fD5urW4t7Ok+MjXapcj6sDEwhrb0X9iybJIdu7cqXY5ohGTAG7EzM3N2bVzO9Y5pyhO2ad2OaIOzNo8QuuRC4iYPJWUlBS1yxGNlARwI+fg4MCBfbspO7GBksvyi9yUWHZ2x6r/DAY/MYIbN26oXY5ohCSAm4AePXrwzZbNFO3+K+U/X1O7HFEHNr8bgKF7f4YMH0VJSYna5YhGRgK4iRg0aBAfvPcXina+i6GkSO1yRB3YBIVxtbwVk56ZIk9GiCokgJuQ2bNmEf7kaIq/+xDFUKF2OeIeaTQm2Ayey8G4H1j89jtqlyMaEQngJmb1yo/QdXGgOPYLtUsRdWBibontqIV8tPoToqK2qF2OaCQkgJsYMzMzdkRvpfXP5yg++a3a5Yg6MLN1oPWoP/HsH2YiPUoFSAA3SXZ2dhzct5uKk99w50Ki2uWIOrBw7E6r0FkMGzGKq1evql2OUJkEcBPVrVs3/vmPaIr3fUxZ3mW1yxF10KpHH9AOY9CwERQXN84BlsTDIQHchPXr149VH/2Vop1LMNwuULscUQc2Ab8n1/wRxodFUFlZqXY5QiUSwE3ctGlTmTE5nKJdy1AqytUuR9wjjUaD7eNzOH76AvMX/EntcoRKZDjKX7C1taWo6OE/Y7to0SL69+/PoEGD7mv/yspKho8aS0JOGTczT2FiYQ0aDaY2bWk/4hVMbdty9ZPpd5eb3P2b6zB4NlbOPbl97ji525fQecYnmLfrAoCiVPLzgc8ouZQMaNCYWdB+zOuY23esr1MW/89wu4Cbm+ex6oP3mDp1yq9ue/36dV566SXi4+Oxt7fH0dGRFStW0KNHj1/dr6HdvHmTjRs3Mnv2bFXraMxkPOB7oFYA14eioiL8evXh/IWLdH7uU0xb2fFz7Fco5SU4DPoDVz+ZTqcpH2Hayq7Kfrk7lmEoysfqUT32wREAFKfFcvvsMdqPXYBGY0LFrTw0FlaYWtmqcWrNXlnuJQq2/Zl93/2Tvn371riNoij06dOHKVOmMHPmTACSk5O5desWwcHBD7PcarKyshg5ciSpqamq1tGYyXjAdXDo0CEGDhzIuHHj8PDwICIiwtiDKT4+nj59+qDX6wkMDKSwsJCSkhKmTZuGl5cXPj4+xMTEALB+/XrGjh3L4MGDcXFxYfXq1SxfvhwfHx969+7NTz/9BMDUqVPZtm0bAC4uLixevBhfX1+8vLxIT08HIDc3l8GDB6PVapkxYwZdu3YlLy/PWLOtrS0H9n4HZbe5c+EkAFZdPH+163Jl2R1Kr56m3RMvUpx+2LjcUPQTprYOaDR3Px5mbdpL+DYgiw5dsR3yIiPHPElWVlaN28TExGBubm4MXwC9Xk+/fv2YN28enp6eeHl5ERUVBdz9DA8YMIAxY8bg6urKggUL2LBhA4GBgXh5eZGZmQnc/ezNnDkTf39/evTowbff3n20sbbP9OnTpwkMDMTb2xudTkdGRgYLFiwgMzMTb29v5s2rPp2WqJ2Z2gU0VqdOneL06dN07tyZvn37cuzYMQIDA5k4cSJRUVEEBARw69YtrK2t+fjjj9FoNPzwww+kp6czZMgQzp07B0BqaiqnTp2ipKQENzc3li1bxqlTp3j55Zf5+uuveemll6odu3379pw8eZK//e1vREZG8vnnn/P2228TGhrKwoUL2bNnD198Ub0jhrOzM46OjhQc/hyLDl25kxmHRQcX4/obm/4EJiZoTM3pNHk5tzNOYNXND3MHJ0ytWlN6/TyWHd1o5RHM9Q3zKbl6Gquuemy1IVg4dm+ot1oA1t39Mdx8ktAhT5CU8G/atGlTZX1qaip+fn7V9ouOjiYpKYnk5GTy8vIICAigf//+wN0r5DNnzuDg4ICrqyszZswgLi6Ojz/+mFWrVrFixQrg7hVsXFwcmZmZhISEcP78edasWVPjZ3rt2rW8+OKLREREUFZWhsFgYOnSpaSmppKUlNTQb1OzI1fAtQgMDMTZ2RkTExO8vb3Jysri7NmzdOrUiYCAAADatGmDmZkZR48e5emnnwbAw8ODrl27GgM4JCSE1q1b06FDB+zs7Bg1ahQAXl5etV7tPPXUUwD4+fkZtzl69ChhYWEADBs2jLZt29a4r4WFBe3s23Djf17BUHwTu97jjescw9+j87RVdJq8HIDbZw5j0/PuL2urnv25nRYL3L3idXruU9r2n4JGY8KNzW9wJyvpft5GUQc2vqO4ZefK2N9PwGC4t0H4jx49Snh4OKampjg6OjJgwADi4+MBCAgIoFOnTlhaWtK9e3eGDBkCVP/sTZgwARMTEx577DFcXV1JT0+v9TMdFBTEe++9x7Jly7h06RLW1tb1+ya0MBLAtbC0tDT+29TUlIqK+xt74ZftmJiYGF+bmJjU2uZ/trnf4548eZKFf/oTFiU/gal5jdsY7hRScimF/D0rufrJdG7FRVOcftR4q0VjZo51d3/ahkynTe8J3Mk4Uec6RN1oNBpsBj5PUtaPvPDSK1XWabVaEhPr1unmXj97Gk3VyVT/+/UvTZo0iZ07d2Jtbc3w4cM5ePBgnWoSVUkA14G7uzs5OTnGK4zCwkIqKioIDg5mw4YNAJw7d47Lly/j7u5er8fu27cvW7bcHUNg3759/Pzzz7+6/V/efouB/p4U/2tVjSNw3T57DBttCM6z/o7zrC9xnr0eM3tHSq+epvT6eSoK84G7T0SU517EtI1MDvowaEzNsBk+n//dup1P1q41Lg8NDaW0tJR169YZl6WkpGBvb09UVBQGg4Hc3FwOHz5MYGBgnY65detWKisryczM5MKFC7i7u9f6mb5w4QKurq7MnTuXMWPGkJKSQuvWrSksLKyfN6CFkQCuAwsLC6KionjhhRfQ6/UMHjyYkpISZs+eTWVlJV5eXkycOJH169dXufqoD4sXL2bfvn14enqydetWOnbsSOvWrWvdXqPRsPF/vsLJrIjiE1HV1hefiaVVj6Aqy1r16ENxWiyVtwvI/eYdrn0xm5wv/wgmprTxG1mv5yNqZ2rdmtaj/8y8BW8YrzA1Gg3bt29n//79dO/eHa1Wy8KFC5k0aRI6nQ69Xk9oaCgffPABHTvW7XHBRx99lMDAQJ544gnWrl2LlZVVrZ/pLVu24Onpibe3N6mpqUyePJl27drRt29fPD095Uu4OpLH0JqI0tJSTE1NMTMz4/jx48yaNeuevvS4fv06Ol9/lICnsemp7uNKom7uXEqmZO9y4k9832DP+k6dOpWRI0cybty4Bmlf3FXbY2jyFEQTcfnyZSZMmEBlZSUWFhZ89tln97Rfx44d2ffdt/QbGIqZ3SNYdq7fWyOi4Vh31WMIDGPQ0OEkn4yv9YtX0XTJFXAL8fqCBXy+M5bWo/+sdimiDgzFN8n9+2wO7d9L79691S5H3CfpiNGCxcXFsWbtOix7T1K7FFEHSkU5RbuW8cKc2TWG7549e3B3d8fNzY2lS5dWW3/p0iUef/xxdDodAwcOrDL85fz589FqtfTs2ZO5c+dW+6J29OjReHp6Gl+/+eab6HQ6vL29GTJkCNeuydyE9UECuJGZPn06jzzySJUP/4O4cuUKT4waQ6vH/4jFI6710qZoeIqiULR/DX08u7Ps/SXV1hsMBubMmcPu3btJS0tj06ZNpKWlVdnmtddeY/LkyaSkpLBo0SIWLlwIwPfff8+xY8dISUkhNTWV+Ph4YmNjjftFR0dja1u15+O8efNISUkhKSmJkSNH8s47MrVSfZAAbmSmTp3Knj176qWtoqIiBg0bgYnncFo91qte2hQPR3HcNh4x5LJl0/9iYlL91zQuLg43NzdcXV2xsLAgLCyMHTt2VNkmLS2N0NBQ4G6HoP+s12g0lJSUUFZWRmlpKeXl5Tg6OgJ3PzPLly/nz3+ueqvqlz3ziouLf/VZYXHvJIAbmf79++Pg4PDA7VRWVjJu4iTyLTvTyv/JeqhMPCy3z34PZ/axf88ubGxsatwmOzubLl26GF87OzuTnZ1dZRu9Xk90dDQA27dvp7CwkPz8fIKCgggJCaFTp0506tSJoUOH0rNnT+DurYZXX32VVq1aVTvmG2+8QZcuXdiwYYNcAdcTCeBmat7rC4hLv4zN4zPlaqUJKb1+ntsxn7B31z9xcnJ6oLYiIyOJjY3Fx8eH2NhYnJycMDU15fz585w5c4arV6+SnZ3NwYMHOXLkCElJSWRmZvLkkzX/wV6yZAlXrlwhIiKC1atXP1Bt4i4J4Gbo739fz2dfbcRmxHw0tXRFFo1PRWE+Rd++z5frPq1x4J1fcnJy4sqVK8bXV69erRbYnTt3Jjo6mlOnTrFkyd37yPb29mzfvp3evXtja2uLra0tTzzxBMePH+f48eMkJCTg4uJCv379OHfuHAMHDqx27IiICL755psHP2EhAdzcHD16lBdefhXb0W9UG/tXNF6V5SUUffs+r86dw4QJ439z+4CAADIyMrh48SJlZWVs3ryZ0aNHV9kmLy/PON3R+++/z/Tp04G7Pd9iY2OpqKigvLyc2NhYevbsyaxZs7h27RpZWVkcPXqUHj16cOjQIQAyMjKM7e7YsQMPD496OvOWTQK4Gblw4QKjxj6FzZAXsWj/qNrliHukKJUU71vJoN7eLH7z3p7TNjMzY/Xq1cb7txMmTECr1bJo0SJ27twJ3B0T2N3dnR49enDjxg3eeOMNAMaNG0f37t3x8vJCr9ej1+uNo/TVZsGCBXh6eqLT6di3bx8ff/zxg520AKQjRqMTHh7OoUOHyMvLw9HRkbfffptnn332N/e7desWer9AirqFYOMr4zY0JUXHNuB8J5MTR2OxsrJSuxzRAKQrchOxadOmOu9TUVHB6KfGU9i2B7YSvk1KcdohTC8cYe/JBAnfFkhuQTQDL7z4CimX87EZOEPtUkQdlGanU3LkS/bv+c74HK5oWSSAm7i/ffIJG6N3YDt8HhoTU7XLEfeo4taPFO5aysav1+Pl5aV2OUIlEsBN2P79+5n/pzexHf1nTGTSzCajsvQ2RTvfY9HC13/zyy/RvEkAN1Fnz57l9xPCsH3iVczbdla7HHGPlEoDRXs/YuTj/Zj32qtqlyNUJgHcBP30008MGjoci94RWD2qU7scUQfFR7/Gzd6ULz/7VHooCnkKoqkpLy9nxJinuNPJB1vdELXLEXVQnLIPq+yT7D4Zj4WFhdrliEZAArgJURSFGX+Yxdn8UlqPnKx2OaIOSi7/QOmJDRz5/ijt2rVTuxzRSMgtiCbkr8s/Yse+Q9gOfUWeeGhCyn/OoWh3JN9EbZIuvKIKCeAmYteuXby1ZCm2o97AxLL6UIGicaosKaJo57ssffdtBg8erHY5opGRAG4CUlNTCXt6Mq1HzMfM7hG1yxH3SKk0ULQ7koljh/PHOXPULkc0QhLAjdyPP/7I4CdGYBU8DUunnmqXI+qg+NAXeDrZ87dVK9UuRTRS8iVcI1ZaWsqwkWMod+mL7e9C1C5H1EHRqV20/imdf/4rDjMz+TUTNZNPRiOlKArPTJ3OpTsW2IaEq12OqIM7F09RkbCV/XEnsLOTMZlF7eQWRCP17pL32HcsEZshc9Fo5D9TU1Gef4XivR/xz39E0717d7XLEY2cXAE3Qt988w0ffLQKu7BlmJi37CEKy3/KRqk0PHA7hoIbaMwtgaq9z8zbdYF66JGmMTVHqaygcOcSPl7+IcHBwQ/cpmj+JIAbmYyMDCZFPIN5917cPnPYuPznmC9VrApauffFsrN7lWVFKfspz7+MrW4I5u2cG+S4DX3eVq1sMLewfOB2SkpLsW7fhWkRE3j2/6f+EeK3SAA3Mh07duT5mX+otnx1jArF/EI38wJCdG2qLFsdcxmAopR9AIwaM5auLi71e2DdS/XSzO3bd/jfr7+irLQEgOjo6Fpn/70fj3RyxlvbneUfflBvbYrmT6YkEg/s/PnzPPbYYwBkZWXRtWtXlSuqmaIoHD582DjT78GDBwkJqZ+nS9LT03F2dsbWVoYFFdXVNiWRBLCoF4qiYGJiYvx3Y/afPxgdHDvy4/UctcsRLUBtASxfr4t6odFo+M8f5+zsbJWr+XVubm4894dZ5N64zsmTJ9UuR7RgEsCi3vj5+QGwaNEilSv5bTOmTwXuzioihFokgEW9y8rKUruE3/TII3fH1LhZUKByJaIlkwAW9a6x3wMGjPerafylimZMAlgIIVQiASyEECqRABZCCJVIAAshhEokgIUQQiUSwEIIoRIJYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgAixZnW3Q0AUH9AFi1Zg2jnhqvckWipZIAFk3C9evXCQsLo3v37vj5+TF8+HDOnTtX47ZZWVl4enoCcOjQIUaOHFllvbdOh7WVFbZuAZSWljF8yOM1trNixQpu375dpzprOp4QtZEAFo2eoig8+eSTDBw4kMzMTBITE3n//fe5cePGfbXn5uZGUmIcE0P92PPdt8yaObPG7e4ngIWoCwlg0ejFxMRgbm7OzF8EpV6vp1+/fsybNw9PT0+8vLyIior61XaKi4uZPn06gYGBhISEMGr4MEJDQzEYDLz22mt4enqi0+lYtWoVK1eu5Nq1a4SEhBjnjdu3bx9BQUH4+voyfvx4ioqKANizZw8eHh74+voSHR3dcG+EaHZkVmTR6KWmphpn2/il6OhokpKSSE5OJi8vj4CAAPr3719rO0uWLCE0NJQvv/ySmzdvEhgYyKBBg/j666/JysoiKSkJMzMzfvrpJxwcHFi+fDkxMTG0b9+evLw83n33Xfbv34+NjQ3Lli1j+fLlzJ8/n+eee46DBw/i5ubGxIkTG/KtEM2MBLBoso4ePUp4eDimpqY4OjoyYMAA4uPj0el0NW6/b98+du7cSWRkJAAlJSVcvnyZ/fv3M3PmTMzM7v46ODg4VNv3xIkTpKWl0bdvXwDKysoICgoiPT2dbt26GWeFfvrpp1m3bl1DnK5ohiSARaOn1WrZtm3bA7ejKArffPMN7u7u97Xv4MGD2bRpU5XlSUlJD1yXaLnkHrBo9EJDQyktLa1yZZmSkoK9vT1RUVEYDAZyc3M5fPgwgYGBtbYzdOhQVq1aZZwy6dSpUwAMHjyYTz/9lIqKCgB++uknAFq3bk1hYSEAvXv35tixY5w/fx64ez/53LlzeHh4kJWVRWZmJkC1gBbi10gAi0ZPo9Gwfft29u/fT/fu3dFqtSxcuJBJkyah0+nQ6/WEhobywQcf0LFjx1rbefPNNykvL0en06HVannzzTcBmDFjBo8++qixrY0bNwLw/PPPM2zYMEJCQujQoQPr168nPDwcnU5nvP1gZWXFunXrGDFiBL6+vsbJPoW4F5q6TKDo7++vJCQkNGA5oqnTaDSEhIRw8OBBtUsRotHQaDSJiqL4//dyuQIWQgiVSAALIYRKJICFEEIlEsCiSarL2BC/5kG6GyclJfHdd9/d175CgASwaILqc2wICWChJglg0eTUdWyIQ4cOMXDgQMaNG4eHhwcREREoilLjeA+zZs3C398frVbL4sWLje3Hx8fTp08f9Ho9gYGBFBQUsGjRIqKiovD29v7NcSiEqJGiKPf84+fnpwjxawAlJCSkQY/x8ccfKy+99FK15du2bVMGDRqkVFRUKNevX1e6dOmiXLt2TYmJiVHatGmjXLlyRTEYDErv3r2VI0eOKIqiKF27dlVyc3ONbeTn5yuKoigVFRXKgAEDlOTkZKW0tFTp1q2bEhcXpyiKohQUFCjl5eXK3//+d2XOnDkNeq6ieQASlBoyVa6ARbNR29gQAIGBgTg7O2NiYoK3tzdZWVk1trFlyxZ8fX3x8fHh9OnTpKWlcfbsWTp16kRAQAAAbdq0MY4bIcSDkAAWTY5WqyUxMbFO+1haWhr/bWpqaux2/EsXL14kMjKSAwcOkJKSwogRIygpKXngeoWojQSwaHLqa2wIqDrew61bt7CxscHOzo4bN26we/duANzd3cnJyTFeTRcWFlJRUVFlXyHuhwSwaHLqa2wIqDreg16vx8fHBw8PDyZNmmQcetLCwoKoqCheeOEF9Ho9gwcPpqSkhJCQENLS0uRLOHHfZCwIUa9kLAghqpOxIIQQopGRABZCCJVIAIt6UVlZSfceHsDdjhJTnn1O5YqEaPwkgEW9MDEx4eUXX8DCqhV2Dh2YOWO62iXdE5l6XqhJvoQT9SolJYV27drh5OSkdin3xMXFhYSEBNq3b3/P+xgMBkxNTRuwKtHcyJdw4qHQ6XSqhO+HH37IypUrAXj55ZcJDQ0F4ODBg0RERNQ4xkNNY0Hs27ePoKAgfH19GT9+PEVFRcDdoH799dfx9fVl69atD/38RPMkASyaheDgYI4cOQJAQkICRUVFlJeXc+TIEfr378+SJUtISEggJSWF2NhYUlJSmDt3Lp07dyYmJoaYmBjy8vJ499132b9/PydPnsTf35/ly5cbj9GuXTtOnjxJWFiYWqcpmhnp0C6aBT8/PxITE7l16xaWlpb4+vqSkJDAkSNHWLlyJVu2bGHdunVUVFSQk5NDWloaOp2uShsnTpwgLS3N2AGjrKyMoKAg4/qJEyc+1HMSzZ8EsGgWzM3N6datG+vXr6dPnz7odDpiYmI4f/481tbWREZGEh8fT9u2bZk6dWqNYzwoisLgwYNrnVrexsamoU9DtDByC0I0G8HBwURGRtK/f3+Cg4NZu3YtPj4+tY7xAFXHgujduzfHjh3j/PnzABQXF9/XLBtC3CsJYNFsBAcHk5OTQ1BQEI6OjlhZWREcHFzrGA9QdSyIDh06sH79esLDw9HpdAQFBZGenq7iGYnmTh5DE0KIBiaPoQkhRCMjASyEECqRABZN2p49e3B3d8fNzY2lS5dWW19aWsrEiRNxc3OjV69exqmI8vPzCQkJwdbWlj/+8Y9V9nnjjTfo0qULtra2VZZfunSJxx9/HJ1Ox8CBA7l69WqDnZdoGSSAW5CQkBD27t1bZdmKFSuYNWtWjdu7uLiQl5f3q22+9957VV736dMHgKysLDw9PYG7HSPmzp17v2UbTZ06lW3bthlfGwwG5syZw+7du0lLS2PTpk2kpaVV2eeLL76gbdu2nD9/npdffpnXX38dACsrK/7yl78QGRlZ7TijRo0iLi6u2vLXXnuNyZMnk5KSwqJFi1i4cOEDn5No2SSAW5Dw8HA2b95cZdnmzZsJDw+/7zb/O4C///77atv4+/sbuwnXp7i4ONzc3HB1dcXCwoKwsDB27NhRZZsdO3YwZcoUAMaNG8eBAwdQFAUbGxv69euHlZVVtXZ79+5Np06dqi1PS0szdnEOCQmpdiwh6koCuAUZN24cu3btoqysDLh7lXrt2jWys7Px8vLC09PTeIX438aOHYufnx9ardY4F9uCBQu4c+cO3t7eREREAFT733aAQ4cOMXLkSACGDx+Ot7c33t7e2NnZ8dVXX2EwGJg3bx4BAQHodDo+/fRT4G7HiD/+8Y+4u7szaNAgfvzxxyrtZmdn06VLF+NrZ2dnsrOza93GzMwMOzs78vPz6/zeAej1eqKjowHYvn07hYWF992WECAB3KI4ODgQGBho7IiwefNmBg0axOuvv87BgwdJSkoiPj6ef/zjH9X2/fLLL0lMTCQhIYGVK1eSn5/P0qVLsba2JikpiQ0bNtxTDd999x1JSUl88cUXdO3albFjx/LFF19gZ2dHfHw88fHxfPbZZ1y8eJHt27dz9uxZ0tLS+Prrr2u8un6YIiMjiY2NxcfHh9jYWJycnGRUNPFAJIBbmF/ehti8eTNdu3Zl4MCBdOjQATMzMyIiIjh8+HC1/VauXIler6d3795cuXKFjIyM+64hLy+PZ555ho0bN2JnZ8e+ffv4+uuv8fb2plevXuTn55ORkcHhw4cJDw/H1NSUzp07G//3/z+cnJy4cuWK8fXVq1erjcT2y20qKiooKCigXbt291V3586diY6O5tSpUyxZsgQAe3v7+2pLCJAAbnHGjBnDgQMHOHnyJLdv38bb2/s39zl06BD79+/n+PHjJCcn4+PjU+NYCvfCYDAQFhbGokWLjF/SKYrCqlWrSEpKIikpiYsXLzJkyJDfbCsgIICMjAwuXrxIWVkZmzdvZvTo0VW2GT16NF999RUA27ZtIzQ0FI1Gc1+15+XlUVlZCcD777/P9OlNY9B50XhJALcwtra2hISEMH36dMLDwwkMDCQ2Npa8vDwMBgObNm1iwIABVfYpKCigbdu2tGrVivT0dE6cOGFcZ25uTnl5+T0ff8GCBeh0uipDOg4dOpRPPvnE2M65c+coLi6mf//+REVFYTAYyMnJISYmpkpbZmZmrF69mqFDh9KzZ08mTJiAVqtl0aJF7Ny5E4Bnn32W/Px83NzcWL58eZVH1VxcXHjllVdYv349zs7Oxico5s+fj7OzM7dv38bZ2Zm33noLuPuHyN3dnR49enDjxg3eeOONez5vIWqkKMo9//j5+Smi6du+fbsCKGfOnFEURVE2btyoeHp6KlqtVpk/f75xu65duyq5ublKSUmJMmzYMMXDw0MZM2aMMmDAACUmJkZRFEWZP3++4uHhoUyaNElRFEWxsbFRFEVRLl68qGi1WkVRFCUmJkYZMWKEoiiKAiharVbR6/WKXq9XduzYoRgMBmXhwoXGGgYOHKjcvHlTqaysVObMmaP06NFDGTRokPLEE08oW7dufVhvkxD1BkhQashUGQtCCCEamIwFIYQQjYwEsBBCqEQCWNSb3xqXQQhRlQSwqBf3Mi6DEKIqCWBRL+5lXAYhRFUSwKJe3Mu4DEKIqiSAhRBCJRLAol7cy7gMQoiqJIBFvbiXcRmEEFWZqV2AaB5+OS6DwWBg+vTpaLVatcsSolGTrshCCNHApCuyEEI0MhLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiV1Ok5YI1GkwtcarhyhBCiWeqqKEqH/15YpwAWQghRf+QWhBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKo5P8A34mmysDSJW0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "Composting.plot_sankey(margin=1.7, offset=.4, gap=0.8)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Sankey Diagrams for PFOA and PFOS"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABWmUlEQVR4nO3dd3xUVfrH8c9JAgESQpEiBgRDaKZCQiQiJVSpiqIEsGBdhBXLimLDhoIu+lMEF3FVdKUpgnREqmBZCBgQItKVLkUgCaSf3x8TsiAQCCk3k3zfr9e8yNxy7nNnhmeee++Ze4y1FhERERERd+DhdAAiIiIiIpdKxauIiIiIuA0VryIiIiLiNlS8ioiIiIjbUPEqIiIiIm5DxauIiIiIuA0Vr4XAGGONMYFOx3EpjDFXG2OSjDGeTsciIuIE5WwR91Kqitfs//CnH1nGmFNnPO9/gXXaGmP2FGAMy40x9xdUe/llrf3dWutrrc0s6LaNMRONMWnZr+9RY8w3xpjG2fNeNMak/+U9efKMdbsbY1YbY5KNMUeMMZOMMbXPs40Xs794rstjbP7GmFnZce0xxgzMZdkYY8zPxphj2bHMNMb4X2A/Tz88s+e1zf6sJf3lEZ3HeL2NMR8ZY04YYw4YYx4/Y16L7Nf2qDHmkDHmC2NMrby0L1IcKWefy11ztjGmrDHmzex8m2SM2WWMeTsPsS3Lzm8njDHrjTE3XWT5ZsaYb7O3ddAY88gZ88KNMSuNMcez43n+jHlFkbOvNcbEGWP+zH4sNsZcm5f2S7tSVbxm/4f3tdb6Ar8DPc6YNsnp+EqoN7Jf79rAH8DEM+ZNO/M9sda+AWCM6Q1MBt4GqgFBQCqwyhhT5fTKxhgD3AUczf43Lz4DdgI1gW7Aa8aYmAssmwB0ttZWBq4CtgL/Ot9+nvE484tl31/m+Vprf/jrRowxA4wxEy8Qw4tAA6AuEAM8aYy5MXteFWACUC97fiLw8YV3XcQ9KGc7orBy9tNAJBAFVATaAuvyENcjQC1rrR/wIPDZhQ7SjTHVgIXA+8AVQCCw6IxFJgPfAlWBNsAgY0zPM+YXds7eB/TO3n41YDYw9SL7L2coVcXrhWQfIb1tjNmX/Xg7e5oPsAC46oyjr6uMMVHGmB+yz8TtN8aMNcaUvYztts0+6vuHMeaP7LbuOWN++ewj1d+yjxBXGWPKZ8/raYzZlB3DcmNMkzPW22WMGWqM2ZB9FPyhMaamMWaBMSYx+yivSvay9YzrzKVX9vPlxphXjDHfZS+7KDsRnG77rux4jhhjns/eVoeL7au19iSuhBF8kdfEAG8CI6y1k621p6y1B4D7gSTgsTMWbwXUAoYAsZf6HhhjfHElzlettenW2vXAdODeC8R+0Fq774xJmbiSYVG6G3jFWvuntfYX4ANgQHZ8C6y1X1hrT2S/zmOBlhdq6K/vmXGdUfnsjOeX9R6LFBXlbLfM2c2BmdbafdZll7X200t64V3xbLDWZpx+CpQB6lxg8ceBr621k6y1qdbaxOy8eVo9YJK1NtNaux1YhavgLki55exj2ftvAcNFvlOUs8+l4tXlWaAFEA6E4ToyfM5amwx04eyjsH24PmiP4TpiigbaA4Muc9tXApUAf+A+YJz535HqaCACuB7XEdqTQJYxpiEwBXgUqA7MB+b8JRnfCnQEGgI9cCX0Z7KX98BV8F1IP+AeoAZQFngCXJc6gPeA/riKxtNxX1R2wdgf+OkiizYCrga+OHOitTYL+DJ7n067G5gDfJ79vMcZ2+tnjNlwoXD+8u/pvy+YpI2rn9kx4BSu1+ONvywyyLgus601xtx6oXYuR/bnoRaw/ozJ67lwsm0NbLrMbV32eyxShJSzz+YOOftH4HFjzCBjTEh20Xvm9t4zxrx3kZjmGmNSgP8Cy4G4CyzaAjhqjPk++yBjjjHm6jPmvw3cZYwpY4xphOszsfgi+3nJLjVnZ3+npADvAq9d5rZKZc5W8erSH3jZWvuHtfYQ8BJw54UWttautdb+aK3NsNbuwnVpos1lbjs9e9vp1tr5uI5UGxljPHCdCXzEWrs3+wjxe2ttKtAHmGet/cZam44rYZbHlTBPezf7jOFeYCXwX2vtT9baFGAm0DSXmD621m6x1p7CVRiGZ0/vDcyx1q6y1qYBw3EdAefmiez/oNsAX7KPPLPdnn0W4vTjKlxfLgD7z9PW/tPzjTEVgNuAydmvwXTO6DqQfQYg9HwBWWsTge+A540x5YwxzXB9cVS40E5k9zOrnL3954DNZ8weg+vyUA3geWCiMebMM59X/WU/j2WfIbpUvtn/Hj9j2nFcl97OYowJxfW+DM1D+2e6nPdYpKgpZ5+t2OdsYCTwOq73Lg7Ya4y5+/SC1tpB1tpcDyistd1x5b2uwKLsAvl8auM6ufEIrsJ6J66Dh9Pm4nptTuHK5R9aa9ecMb9Icnb2d0ol4O9c/CDhQkplzvZyOoBi4irgtzOe/5Y97byyj6LfwtV/pwKu13HtZW77yBmXQgBO4vrgVwPKAdsvFq+1NssYs5uzj7YOnvH3qfM89+XCDpwnntPb3X3Gdk8aY47k0g7AaGvtcxeY97m19o4zJxhj/LL/rIUr4ZypFnA4++9eQAauMxgAk4DFxpjq2V9mF9MfGIdrf3bg6gN70ctG1tqjxphPgPXGGP/sL8Mz+23NN8ZMAm7BVSCD6yzQOT82A9fZBlxnTcB1xsTLGHNz9vPfswvwpOznfriO0k//nfiXtgJxna15xFq78mL7cgGX8x6LFDXl7LMV+5xtXb8DGIfrTHV5XIX+R8aY1X+5pJ+r7OJ/gTHmEWPMNmvt7PMsdgpXF4U12TG+BBw2xlQCPHH1h/07rm4RVwLTjTEHrbWnz/wWSc7O3p9kY8x44JAxpom19o9LeyVylMqcrTOvLvtwdao+7ersaXD+I5h/4Tpaa2Bdncef4exL0AXhMK4Pff3zzDsr3uzLL3WAvQUcw1/tx3VEe3q75XF1hi9IvwJ7cJ1VzZF9VuNWYEn2pLtxJejfjTEHcF2yKsP/kkqurLW/WWu7W2urW2uvw/XFs/oSY/TCdZbV7wLzT/djupQ4BllrK2cfgQ/CdSa5cvYjNHuZP3G99mFnrBrGGV0DjDF1cV32esVa+59L2HT5M/6udMbfRfEei+SXcvalKU45O4d19YsdB/wJXO6v7L04/2sNsIGzPwdn/h0AZFprP80++bAH14+lul7KRgsqZ/+FB66Dqtwu9ytnn0HFq8sU4DljTPXsju7DcZ2JA9fR7xXZR2ynVQROAEnGdRuRhwo6oOzLIR8BbxnXDw48jTHRxhhvXJeFuhlj2htjygD/wPXLzu8LOo6/mA70MMZcn91X60UK+AvAWmtx9dd6LrvfajljzJXAv3EVi/9nXLepag90x3V5LBxXYnidS7zrgDGmiTGmonHdvuUOoBOuMzPnW/YWY0wjY4yHMaZ69nI/WWuPZs/vbYzxzZ7fCbgD169HC9KnuF6TKtmfuQfI/hVw9uuxFBhrrR1/ie3dYVw/cLkG1xdMxezPUqG/xyIFQDn70hSLnA1gjHnUuH7wVt4Y45XdZaAil3C53BjT2BjTJXvdMtk5uzWw4gKrfAz0Mq5bYpXB1Z1rlbX2OLDF1aTpl52zr8TVreNCv5G4XLnl7I7GmKbZnxE/XN8pfwK5nYFWzj6DileXEbj64GwAfsZ1+44RANbazbgS5Q7zvz4+T+A6w5eI6xeE0wopriey41mD63ZQrwMe1tpfcRVI7+I62u+B6xYyaYUUBwDW2k3Aw7iOUvfjujTyB64kXJDbmYar/9pjwBFct6oqD7S01h7JnhdvrV1krT1w+oGr72moMSbYGNPfGJPbj5Y64+ou8CcwELjxzO4GxvUr5VbZT/1xXWZKxPV+ZOHqtnDaI7jOoBwD/gk8YK1dfsb8M3/5fPqR1x91vYDrcuRvuBL2P621C7Pn3Y/rbMKLZ27jIu0l4brEtwDXDwVaAXcW1Xsskk/K2ZegGOVscHVneBNXF4fDwGDgVmvtDgBjzPjsy+fnY3AVZX8Ah3Dl3D6nu2wZY1qdmfOstUtxnV2fl71OINlX5ay1J3B163oMV/6PBzaS/fnJVtg5uzKuz+jx7GXq4/oOSjlPO6cpZ5/BuA6aRPLOuH6NegzXpbi/9nWSYsoYswu431p70V/X6j0WKTn0/9k9KWefS2deJU+MMT2MMRWM65eXo3GdZdjlbFRSkPQei5Qc+v9c8pXG91jFq+TVTbh+fLAP1+2hYq1O35c0eo9FSg79fy75St17rG4DIiIiIuI2dOZVRERERNxGnopXY8zCiy8lIlL8KH+JiJQMeRphy8/Pr3NkZKT6GYiIOzrhdABFrVq1arZevXpOhyEikmdr1649bK2tfr55eSpeGzRoQFxcXMFEJSJShIwxW52OoajVq1dPOVtE3JIx5rcLzVOfVxERERFxGypeRURERMRtqHgVEREREbeh4lVERERE3IaKVxERERFxGypeRURERMRtqHgVEREREbeh4lVERERE3IaKVxERERFxGypeRURERMRtqHgVRxw4cIDY2Fjq169PREQEXbt2ZcuWLU6HxbFjx3jvvfecDkNEpFhRzpbiRMWrFDlrLb169aJt27Zs376dtWvXMnLkSA4ePOh0aEqEIiJ/oZwtxY2KVylyy5Yto0yZMgwcODBnWlhYGDfccANDhw4lODiYkJAQpk2bBsDy5ctp06YNN910EwEBAQwbNoxJkyYRFRVFSEgI27dvB2DAgAEMHDiQyMhIGjZsyNy5cwFISUnhnnvuISQkhKZNm7Js2TIANm3aRFRUFOHh4YSGhrJ161aGDRvG9u3bCQ8PZ+jQoUX8yoiIFD/K2VLceDkdgJQ+GzduJCIi4pzpM2bMID4+nvXr13P48GGaN29O69atAVi/fj2//PILVatWJSAggPvvv5/Vq1fzzjvv8O677/L2228DsGvXLlavXs327duJiYlh27ZtjBs3DmMMP//8M5s3b6ZTp05s2bKF8ePH88gjj9C/f3/S0tLIzMxk1KhRbNy4kfj4+CJ8RUREii/lbCludOZVio1Vq1bRt29fPD09qVmzJm3atGHNmjUANG/enFq1auHt7U39+vXp1KkTACEhIezatSunjdtvvx0PDw8aNGhAQEAAmzdvZtWqVdxxxx0ANG7cmLp167Jlyxaio6N57bXXeP311/ntt98oX758ke+ziIi7Us4Wp6h4lSIXFBTE2rVr87SOt7d3zt8eHh45zz08PMjIyMiZZ4w5a72/Pj9Tv379mD17NuXLl6dr164sXbo0TzGJiJQGytlS3Kh4lSLXrl07UlNTmTBhQs60DRs2ULlyZaZNm0ZmZiaHDh3i22+/JSoqKk9tf/HFF2RlZbF9+3Z27NhBo0aNaNWqFZMmTQJgy5Yt/P777zRq1IgdO3YQEBDAkCFDuOmmm9iwYQMVK1YkMTGxQPdXRMSdKWdLcaPiVYqcMYaZM2eyePFi6tevT1BQEE8//TT9+vUjNDSUsLAw2rVrxxtvvMGVV16Zp7avvvpqoqKi6NKlC+PHj6dcuXIMGjSIrKwsQkJC6NOnDxMnTsTb25vPP/+c4OBgwsPD2bhxI3fddRdXXHEFLVu2JDg4WJ3/RURQzpbix1hrL3nhyMhIGxcXV4jhiFy+AQMG0L17d3r37u10KFIMGWPWWmsjnY6jKClnS3GmnC25yS1n68yriIiIiLgN3SpLSoyJEyc6HYKIiFwi5Wy5XDrzKiIiIiJuQ8WrFClPT0/Cw8MJDg6mR48eHDt2LM9txMXFMWTIkPPOq1evHocPH76s2F588UVGjx59WeuKiJREytlSHKl4lSJVvnx54uPj2bhxI1WrVmXcuHF5biMyMpIxY8YUQnQiInIm5WwpjlS8imOio6PZu3cvANu3b+fGG28kIiKCVq1asXnzZsB1D8Dg4GDCwsJyhh1cvnw53bt3B+DIkSN06tSJoKAg7r//fk7fPWPXrl0EBwfnbGv06NG8+OKLAHzwwQc0b96csLAwbr31Vk6ePHlObGPGjOHaa68lNDSU2NjYQnsNRETchXK2FBcqXsURmZmZLFmyhJ49ewLw4IMP8u6777J27VpGjx7NoEGDAHj55Zf5+uuvWb9+PbNnzz6nnZdeeokbbriBTZs20atXL37//feLbvuWW25hzZo1rF+/niZNmvDhhx+es8yoUaP46aef2LBhA+PHj8/n3oqIuDflbClOdLcBKVKnTp0iPDycvXv30qRJEzp27EhSUhLff/89t912W85yqampALRs2ZIBAwZw++23c8stt5zT3rfffsuMGTMA6NatG1WqVLloDBs3buS5557j2LFjJCUl0blz53OWCQ0NpX///tx8883cfPPNl7m3IiLuTTlbiiOdeZUidbr/1G+//Ya1lnHjxpGVlUXlypWJj4/Pefzyyy8AjB8/nhEjRrB7924iIiI4cuTIJW3Hy8uLrKysnOcpKSk5fw8YMICxY8fy888/88ILL5w177R58+YxePBg1q1bR/Pmzc8ai1tEpLRQzpbiSMWrOKJChQqMGTOGN998kwoVKnDNNdfwxRdfAGCtZf369YCrX9V1113Hyy+/TPXq1dm9e/dZ7bRu3ZrJkycDsGDBAv78808AatasyR9//MGRI0dITU1l7ty5OeskJiZSq1Yt0tPTc8bPPlNWVha7d+8mJiaG119/nePHj5OUlFQor4OIiDtQzpbiRN0GxDFNmzYlNDSUKVOmMGnSJB566CFGjBhBeno6sbGxhIWFMXToULZu3Yq1lvbt2xMWFsaKFSty2njhhRfo27cvQUFBXH/99Vx99dUAlClThuHDhxMVFYW/vz+NGzfOWeeVV17huuuuo3r16lx33XUkJiaeFVdmZiZ33HEHx48fx1rLkCFDqFy5cpG8JiIixZVythQX5vQv/S6FxskWEXeV2zjZJZVytoi4q9xytroNiIiIiIjbUPEqIiIiIm5DxauIiIgUmYULF9KoUSMCAwMZNWqU0+GIG1LxKm5PiVBExD1kZmYyePBgFixYQEJCAlOmTCEhIcHpsMTNqHgVt6ZEKCLiPlavXk1gYCABAQGULVuW2NhYZs2a5XRY4mZUvIpbUyIUEXEfe/fupU6dOjnPa9euzd69ex2MSNyRildxa0qEIiIipYuKVzd24MABYmNjqV+/PhEREXTt2pUtW7bkuZ23336bkydPXlYM8fHxzJ8//7LWFREpTZSzwd/f/6xRt/bs2YO/v79j8Yh7UvHqpqy19OrVi7Zt27J9+3bWrl3LyJEjOXjwYJ7bUiIUESlcytkuzZs3Z+vWrezcuZO0tDSmTp1Kz549HYtH3JOKVze1bNkyypQpw8CBA3OmhYWFccMNNzB06FCCg4MJCQlh2rRpACxfvpy2bdvSu3dvGjduTP/+/bHWMmbMGPbt20dMTAwxMTEAPPTQQ0RGRhIUFMQLL7yQ0/6aNWu4/vrrCQsLIyoqiuPHjzN8+HCmTZtGeHh4zraKkhKhiLgD5WwXLy8vxo4dS+fOnWnSpAm33347QUFBRR6HuDlr7SU/IiIirBQP77zzjn300UfPmT59+nTboUMHm5GRYQ8cOGDr1Klj9+3bZ5ctW2b9/Pzs7t27bWZmpm3RooVduXKltdbaunXr2kOHDuW0ceTIEWuttRkZGbZNmzZ2/fr1NjU11V5zzTV29erV1lprjx8/btPT0+3HH39sBw8eXAR7fGHz5s2zDRo0sAEBAXbEiBGOxiLFFxBn85DvSsJDObv4UM4WyZvccraX08WzFKxVq1bRt29fPD09qVmzJm3atGHNmjX4+fkRFRVF7dq1AQgPD2fXrl3ccMMN57Tx+eefM2HCBDIyMti/fz8JCQkYY6hVqxbNmzcHwM/Pr0j3Kzddu3ala9euTochIpJnpTFni+SXug24qaCgINauXZundby9vXP+9vT0JCMj45xldu7cyejRo1myZAkbNmygW7dupKSk5DteEZHSTDlbpOCoeHVT7dq1IzU1lQkTJuRM27BhA5UrV2batGlkZmZy6NAhvv32W6KionJtq2LFiiQmJgJw4sQJfHx8qFSpEgcPHmTBggUANGrUiP3797NmzRoAEhMTycjIOGtdERE5P+VskYKj4tVNGWOYOXMmixcvpn79+gQFBfH000/Tr18/QkNDCQsLo127drzxxhtceeWVubb14IMPcuONNxITE0NYWBhNmzalcePG9OvXj5YtWwJQtmxZpk2bxsMPP0xYWBgdO3YkJSWFmJgYEhISHOv8LyLiDpSzRQqOcfWJvTSRkZE2Li6uEMMRESkcxpi11tpIp+MoSsrZIuKucsvZOvMqIiIiIm5DxauIiIiIuA0VryIiIiLiNlS8Sr6GGiwOrLVs2rSJvPTfFhFxVyUhZ3/55Zds3rzZ6VDETal4lctKhJmZmYUUTd6kpqbS7867CQtvSt877iI1NdXpkERECpU75+z09HTuvf9BBjz0CFHRLVm4cKHTIYkbUvFagvzzn/9kzJgxADz22GO0a9cOgKVLl9K/f//zjn99vnGyFy1aRHR0NM2aNeO2224jKSkJgHr16vHUU0/RrFkzvvjiCwf28FzPPPc8Uyf9h8yMdKZN/oynnnnW6ZBERC5JaczZ33//PRM/+jdJh/aSeOwog/4+xOmQxA2peC1BWrVqxcqVKwGIi4sjKSmJ9PR0Vq5cSevWrXn11VeJi4tjw4YNrFixgg0bNjBkyBCuuuoqli1bxrJlyzh8+DAjRoxg8eLFrFu3jsjISN56662cbVxxxRWsW7eO2NhYp3bzLI8/+ghLly4FYPHixQx9/DGHIxIRuTSlMWe3bt0656zxXXfdxS+bfnY4InFHXk4HIAUnIiKCtWvXcuLECby9vWnWrBlxcXGsXLmSMWPGnHf869DQ0LPa+PHHH0lISMi50XVaWhrR0dE58/v06VOk+3Qx/v7++Pv7AxATE4OHh47HRMQ9lMacbYyhfPnygGv42zOHwBW5VCpeS5AyZcpwzTXXMHHiRK6//npCQ0NZtmwZ27Zto3z58owePZo1a9ZQpUoVBgwYcN7xr621dOzYkSlTppx3Gz4+PoW9GyIipYJytsjl0WmqEqZVq1aMHj2a1q1b06pVK8aPH0/Tpk0vOP41nD1OdosWLfjuu+/Ytm0bAMnJyWzZssWRfRERKemUs0XyTsVrCdOqVSv2799PdHQ0NWvWpFy5crRq1eqC41/D2eNkV69enYkTJ9K3b19CQ0OJjo7W7UxERAqJcrZI3pm83BtT42RLcWWMITMzU31e5YJyGye7pFLOluLKGMMDDzzAhAkTnA5Fiqnccra+6UVERETEbah4FRERERG3oeK1hFm4cCGNGjUiMDCQUaNGnTM/NTWVPn36EBgYyHXXXceuXbty5o0cOZLAwEAaNWrE119/nTO9Xr16hISEEB4eTmRkqbrqKiJS4AojT1+ozVatWhEeHk54eDhXXXUVN998c2HumkjRsNZe8iMiIsJK8ZWRkWEDAgLs9u3bbWpqqg0NDbWbNm06a5lx48bZv/3tb9Zaa6dMmWJvv/12a621mzZtsqGhoTYlJcXu2LHDBgQE2IyMDGuttXXr1rWHDh0q2p3JI8BmZmY6HYYUY0CczUO+KwkP5ezipzDy9KW0aa21t9xyi/3kk08KfycvAWAfeOABp8OQYiy3nK0zryXI6tWrCQwMJCAggLJlyxIbG8usWbPOWmbWrFncfffdAPTu3ZslS5ZgrWXWrFnExsbi7e3NNddcQ2BgIKtXr3ZiN0RESqzCyNOX0uaJEydYunSpzrxKiaDitQTZu3cvderUyXleu3Zt9u7de8FlvLy8qFSpEkeOHMl1XWMMnTp1IiIiQr8MFRHJh8LI05fS5ldffUX79u3x8/MrjN0SKVIaYUsuatWqVfj7+/PHH3/QsWNHGjduTOvWrZ0OS0RELtGUKVO4//77nQ5DpEDozGsJ4u/vz+7du3Oe79mzB39//wsuk5GRwfHjx7niiityXff0vzVq1KBXr17qTiAicpkKI09frM3Dhw+zevVqunXrVli7JVKkVLyWIM2bN2fr1q3s3LmTtLQ0pk6dSs+ePc9apmfPnnzyyScATJ8+nXbt2mGMoWfPnkydOpXU1FR27tzJ1q1biYqKIjk5OWcYwuTkZBYtWkRwcHCR75uISElQGHn6Ym1Onz6d7t27U65cuSLdV5HCom4DJYiXlxdjx46lc+fOZGZmcu+99xIUFMTw4cOJjIykZ8+e3Hfffdx5550EBgZStWpVpk6dCkBQUBC333471157LV5eXowbNw5PT08OHjxIr169ANcZgH79+nHjjTc6uZsiIm6rMPI0cN42T5s6dSrDhg1zZH9FCoOGh5USQcPDysVoeFiR4kPDw8rFaHhYERERESkRVLyKiIiIiNtQ8SoiIiIibkPFq4iIiIi4DRWvIiIiIuI2VLyKiIiIiNtQ8SoiIiIibkPFq4iIiIi4DRWvpdjBgwe56dbbmDp1mtOhiIjIJZg3bx6TJk92OgwRR6l4LaV+/vlnQptGsHTjXu6+9z6eG/4CeRltTUREio61lldGvMpt/e9m4CNDeWjwEDIyMpwOS8QRKl5LqZXf/YCHpxdJW37A2CzmL1zkdEgiInIBJ0+eZPSbb3Lq+BGSDu/j44//zf79+50OS8QRKl5LqUEDH2T/7l0AjHztVdat/gFjjLNBiYjIefn4+HD8z6PUrl0bgJSTJ6lTp47DUYk4Q8WriIiIiLgNFa8iIiIi4jZUvIqIiIiI21DxKiIiIiJuQ8WriIiIiLgNFa8iIiIi4jZUvIqIiIiI21DxKiIiIiJuQ8WriIiIiLgNFa8iIiIi4jZUvIqIiIiI21DxKiIiIiJuQ8WriIiIiLgNFa8iIiJuwlrrdAgijlPxWoqdToLVqlVzOJL8SU5OBsAY43AkIiKF6+jRo06HUGD8/PycDkHclIrXUmzlypUAhIeHOxtIPs2ePRtQ8SoiJd+pU6e44YYbnA4jX06ePAlA+/btHY5E3JWK11Ksa7duAISEhDgcSf7069fP6RBERArd999/D8AXX3zhcCT5c8899wDQpUsXhyMRd6XitZSaPn06yUlJvPvuu06Hki/ffPMNAD/88IPDkYiIFJ7MzExatmwJwJVXXulwNJdv3759fP7558TGxjodirgxFa+lTFpaGqNGjeK2226jYaPGDB482OmQLtv8+fPp1KkToaGhtGjRolC28dNPP9H3jrtITU0tlPZFRC7mwIEDeHl5AXD48GGHo7l869atw9/fH4DJkycXyja633wr//7ww0JpW4oPL6cDKO4OHDjAhg0bzpl+4sQJbrvttiKJwdPTEw8Pj3z16WzSpAn7Dxzgj4MHAWjUuDELFyxg9+7dBRVmkVizZg1fffUVn332GeA6AzFnzhx+//33QtneypUrmTrpPwBM/s8n6lcr4gZ++uknsrKyzpn+wQcf8PPPP19SG8YYjDE5uffMx/mm5Tb9cqSmprJ582Z27tyZM+3zzz9n06ZNl9WeU3788UfGjh171nfN0qVLc35zUdC+/24Vi75eQMMGDWjdunWhbEOcZ/Jy243IyEgbFxdXiOEUP88Nf4HXXnuVag2aYvhfEkpLPs6x3VscjOzyeXh64VPZPe8wkHjkwFnPK15R+JfPygV3JG3Hfxn60D08+/SwQt+eFA5jzFprbaTTcRSl0pizt23bRoMGDfCrWYcy5XzOmnfkt80ORZU/xsOTagHu+duEQ9vic/6uUKUmPlfUKtTt2bIVMI1iSFvxPmtX/0j9+vULdXtSeHLL2TrzehHDnhzK59Nn8GetUCo0vyVnenmgknNhlVpVHdpu2WvbM3L0k1zbuBG9evVyKAoRuZjAwEAGPfwIUxauxOem4RjP/33N+ToYV2lV16HtZiUfpUPnrsSvXU2lSvq2LmnU5/UifH19WfL1fLJ+nsfJrf91OhxxiJdfNSp2H8ad99xHfHy80+GISC7G/N+bNK1Xg+TlH+im/qWUb9NuJFZtTPebbyUjI8PpcKSAqXi9BHXq1GHB3NmcXDyWtIM7nA5HHOJdqyHl2zxIpy7d2L9/v9PhiMgFeHp6MnP6NCqd2EHyurlOhyMO8Wl7H5v2HWfQw0OcDkUKmIrXSxQVFcWEf40jcc5rZCb96XQ44hCfJq3IatiOTl17cOrUKafDEZEL8PPzY8nX88n8aQantpeufr/iYjw88e3yBNO+ms+7Y8c5HY4UIBWvedCvX1+GPPQASXNHkpWuWyeVVj4t+rAvy49+d96tS5IixVi9evWYO2smyd+MIe3Qb06HIw7wKOeLb8/nePr5F1i0aJHT4UgBUfGaRyNefonWEdeSvHisCpdSyhiDT4e/szxuE8+/8KLT4YhILlq2bMl7Y94mac5rZJ487nQ44oAyVWrh23UovWP7sXmze95xQs6m4jWPjDFM/exTrvI4QfKP05wORxziUcYb3+5P8857E5gyZarT4YhILu66607+ds+dJM0dhc1IdzoccUC5OsF4R99Bh85dOXLkiNPhSD6peL0M5cuXZ9H8OXhsXUbyL4Vzo2Up/jx9q1Cx5zM88NAgVq9e7XQ4IpKL10e+SnRQAMlL39NVs1LKJ6QjKf4RdOlxM2lpaU6HI/mg4vUy1apVi0Xz53JqxQRS97vnYAWSf2VrBFCh/d/p0r2n241WJlKaeHh48MXUSVRP/4PkNTOcDkcc4nPDnWw7lsW9D/xNBzFuTMVrPoSHh/Ofjz8kce4oMk6473jTkj8VGlyHR0g3OtzYjaSkJKfDEZEL8PHxYfHCebBpASe3/OB0OOIA4+GJ742PMXfpKt4Y/abT4chlUvGaT7169WLYPx4lac5rZKWlOB2OOKRCZC+OeF9F7z79zjumuogUD/7+/iycN4eTS98j7eB2p8MRB3iULY9vj2d4ZeQbzJkzx+lw5DKoeC0Azz49jBtbR5H09f9hrQqX0sgYg0/7gaze/DtPPPmU0+GISC4iIyP5aML7JM4ZSUbSUafDEQd4+dWgYren6HfXADZs2OB0OJJHKl4LgDGGTz76NwEVLcmrPnM6HHGI8SyDT7cn+fenU/joo4+dDkdEcnH77bfx+MMPkaz7dpda3v6NKdfqPjp26cbBgwedDkfyQMVrAfH29mbBnK8o8/uPJG9c4nQ44hDPCpXw7fksQx5/glWrVjkdjojk4sXhzxMTFUryN+/qxzullM+1bcis35rO3XqSkqKuf+5CxWsBql69OosXzif1u09I2bPJ6XDEIWWrXY1Pp0focfMt7Nixw+lwROQCjDFM/nQitcskk/zDFKfDEYf4RMeyO608/e8aoIMYN6HitYAFBQUxbfJnJM37J+nHDjgdjjikfEAEXs1upX3nrhw/rlF9RIqrcuXKsWj+HDy3f0tywgqnwxEHGOOBT8chLPnvBl56ZYTT4cglUPFaCLp06cLLw58jac6rZKWedDoccUiFpt1IrNyAm265jYyMDKfDEZELqFmzJt8smEfKyg9J3fer0+GIAzzKeOPb42neHDOOzz//wulw5CJUvBaSxx4dwi1dOpC0YDQ2K9PpcMQBxhh82t7Pht1H+fuQx5wOR0RyERoayqRPPs6+b/chp8MRB3j5VqVij2e498G/ERcX53Q4kgsVr4XEGMOEf42jSY0KJH870elwxCHG0wvfrkOZMnM24957z+lwRCQXPXv25NmnnnBdNUs75XQ44oCyNetTod1D3NitB3v37nU6HLkAFa+FqEyZMsydNYMKBzeQvH6h0+GIQzzK+eLb8zmeenY4ixcvdjocEcnFsCeH0i2mJcm6b3epVaHh9RB0Ix1u7EZycrLT4ch5eDkdwPn4+vo6Mszm8OHDad26NR06dMh3W/Xq1aNixYoYY7ja/0rW/zAJj8q1ODL/HTzKlgcP13FD1Y6DKFe7CWmHfuPo4vfJTDoCNgufoHZUuj4WY0xOm/s+fpgyVWtT/SbdBN/dlKlyFb5d/sGtt8ey+ofvaNSo0SWve+DAAR599FHWrFlD5cqVqVmzJm+//TYNGzYsxIgv7tixY0yePJlBgwY5GocUDyUpbwNk7N1D4pIJ+HUYyJ5/3ZvnvJ2Z/CdHFoxxdUHIysSzUg1q3vZSvmOUouHT/Fb+WLSH22L7M3fWDDw8Lv1cn3J24SuWxatTXn755QJtb9myZVSrVo1nnnmGGjVqsGjBm5BlqNn3NTwrVMpZLis9lUMzXqFqp0GUv6YZWekpHJo5kqSf5lGxWXcA0g/vhqwsUvdsIistBY+y5Qo0Vil85a4OJbNFPzp07sr6dWuoWrXqRdex1tKrVy/uvvtupk6dCsD69es5ePBgsUiE7733XolIhOK+CitvP/bYY3zw7w9JrhkIkOe8fWzlJMrVC8cv8iYA0v7YWaBxSuEyxuDbfjA/zBjOsGee441Rr13SesrZRaNYdxtYvnw5bdu2pXfv3jRu3Jj+/fvn3INtzZo1XH/99YSFhREVFUViYiIpKSncc889hISE0LRpU5YtWwbAxIkTufnmm+nYsSP16tVj7NixvPXWWzRt2pQWLVpw9KhreMABAwYwffp0wHUE/sILL9CsWTNCQkLYvHkzAIcOHaJjx44EBQVx//33U7duXQ4fPpzrfrRu3ZqUlBRGv/4aWSePkZly9tmJ5IQVePs3ofw1zQDwKFOOqh0HcvzH6f9b5pcV+ATHUK5eM05t+7EAXl1xgk9oZ07VCqdbz16kp6dfdPlly5ZRpkwZBg4cmDMtLCyMG264gaFDhxIcHExISAjTpk0DXP9n2rRpw0033URAQADDhg1j0qRJREVFERISwvbtrrHcBwwYwMCBA4mMjKRhw4bMnTsX4IL/hzZt2kRUVBTh4eGEhoaydetWhg0bxvbt2wkPD2fo0KEF/VKJmyopebtz5840axpO6vefYjPTzpl/sbydmXwUz4rVcpYvW+Oa/Lys4gDjVQbf7sP410ef8Omn/7mkdZSzi0axLl4BfvrpJ95++20SEhLYsWMH3333HWlpafTp04d33nmH9evXs3jxYsqXL8+4ceMwxvDzzz8zZcoU7r777pwRMzZu3MiMGTNYs2YNzz77LBUqVOCnn34iOjqaTz/99LzbrlatGuvWreOhhx5i9OjRALz00ku0a9eOTZs20bt3b37//feL7sPcuXMJCQlh4N/+hq9PBQ5+8gj7Pvo7+z99HID0w79R9srAs9YpU6UWNu1Uzq22kjevxKdJa3yuba17Ebo5n1Z38+vRNO57cOBFb4i9ceNGIiIizpk+Y8YM4uPjcz7/Q4cOZf/+/YDrKH/8+PH88ssv/Oc//2HLli2sXr2a+++/n3fffTenjV27drF69WrmzZvHwIEDSUlJueD/ofHjx/PII48QHx9PXFwctWvXZtSoUdSvX5/4+Hj++c9/FuyLJG6tpOTt6Ohovpg6GXvyBAcmPcm+jx++5LxdsWk3jiwYw4EpT3P8+2lkJB7Jz0sqDvGsUImKPZ5l0JBH+f777y+6vHJ20Sj2xWtUVBS1a9fGw8OD8PBwdu3axa+//kqtWrVo3rw5AH5+fnh5ebFq1SruuOMOABo3bkzdunXZsmULADExMVSsWJHq1atTqVIlevToAUBISAi7du0677ZvueUWACIiInKWWbVqFbGxsQDceOONVKlS5YKxx8TEEB4ezokTJ3j66acBqFq1KtdHX49f3Wu58s43L+k1SN2/Fc/yfnj51aBc3TDSDu4g81TiJa0rxY/x8MS38+PM/mYFb771f5fVxqpVq+jbty+enp7UrFmTNm3asGbNGgCaN29OrVq18Pb2pn79+nTq1Ak497N+++234+HhQYMGDQgICGDz5s0X/D8UHR3Na6+9xuuvv85vv/1G+fLl8/ciSIlWkvJ2586dqVylMuXKeHBlv1HUuuutS3oNygdE4P+3f1MxrDPpR/ewf+IjZJ7UgCXuqGz1uvh0HEL3m3pd8HN3McrZBavYF6/e3t45f3t6el72zd7PbMfDwyPnuYeHxwXbPL3M5W532bJlxMfH8+mnn1K5cuWc6Z9+/G8qndhJ8ro5AJSpdjVpB7adtW76sQOYsuXx8K5A8i8rSD+yhz3/upe97z9AVtpJTv76XZ7jkeLDw7sCvj2e5cVXRzFv3rwLLhcUFMTatWvz1PalftbP/DHg+Z6fqV+/fsyePZvy5cvTtWtXli5dmqeYpHQpaXnbr2JFbu7cjqT5/7tv98XyNoBn+Yr4XNuWat3/gXetBqTs3pjneKR4KF8/Eo/wm2nfuSsnTpy44HLK2UWj2Bev59OoUSP279+fc9SSmJhIRkYGrVq1YtKkSQBs2bKF33//PU+/6r4ULVu25PPPPwdg0aJF/Pnnn3luo2LFiiz5ej6ZP83k1PY4fK5tS+qeBE7tigdcPwT4c/H7+F13K9ZmcXLzKmrdO5baD31E7Yc+osYtz5P8y7cFuVviAK9KNfDt+iSx/e9k48bzf6m1a9eO1NRUJkyYkDNtw4YNVK5cmWnTppGZmcmhQ4f49ttviYqKytP2v/jiC7Kysti+fTs7duygUaNGF/w/tGPHDgICAhgyZAg33XQTGzZsoGLFiiQm6gqAXBp3z9v/fH0UQVf5kbziI4Bc8zbAqd/Wk5Xu6v6QlXqS9D/34+VXvQD2Rpzi06wHx/2uoVfvPmRmnn/wIeXsouGWxWvZsmWZNm0aDz/8MGFhYXTs2JGUlBQGDRpEVlYWISEh9OnTh4kTJ551RFMQXnjhBRYtWkRwcDBffPEFV155Zc6tVfKiXr16zJ01k+RvxpBxbD/Vb32e499PZe8Hf2P/R3+n7JUNqNisO6m7N+HpewVeFa/IWde7ThDpR34nI+loQe6aOKBc7SaUa30vHW7syh9//HHOfGMMM2fOZPHixdSvX5+goCCefvpp+vXrR2hoKGFhYbRr14433niDK6+8Mk/bvvrqq4mKiqJLly6MHz+ecuXKXfD/0Oeff05wcDDh4eFs3LiRu+66iyuuuIKWLVsSHBzs9p3/pfC5e9728vJizszp+B5OIDl+AR5lvC+YtwHSDmzjwCePse+jv3PgsyeoGNYJ71rO/tpc8sc1auKD/LTrDx55/IkLLqOcXfjMxX4wcqbIyEhb2odMS01NxdPTEy8vL3744Qceeugh4uPjL7u9Tz/9D4OfGEbl2H+edRsWKV2SVn1GndSd/LhqRYF/cZ/PgAED6N69O7179y70bRUXxpi11tpIp+MoSsrZLgWZt7dt20bkddF4d3yU8vXCCzROcQ+ZKUkcn/ok/3z5Wf724INFsk3l7LO55ZlXJ/3+++80b96csLAwhgwZwgcffJCv9u66604G3nMXSXNHYTMufuskKZl8WvZj1wnL8y/qJuYiBa0g83ZgYCCzZkwneeFbpB/ZU4BRirvwLOdLxZ7P8o8nn865NZUULZ15LQY2btxIRGQU1e9/Hy/fi9+4Xkqe9GMHOD5tGF9O/YzOnTs7HU6JpDOvUlCstXS6sSvrkitR8YY7nQ5HHHLs67F0Ca7BlP984nQoJZLOvBZjhw4dolPX7vh1fEiFaymVlXqSpDmv8soLz+e7cF24cCGNGjUiMDCQUaNGnTP/t99+o3379oSGhtK2bVv27PnfmaMnn3ySoKAgmjRpwpAhQy56D1qR0urtd95l9cZf8Wl+q9OhiEOSE5ZR7tAm3nnTve+X6q5UvF7EvffeS40aNQgODi7wtlNTU7mx+02kXd0Cn6B2Bd6+FH82K5OkBaPp3bUjjz7ycL7ayszMZPDgwSxYsICEhASmTJlCQkLCWcs88cQT3HXXXWzYsIHhw4fn3H/4+++/57vvvmPDhg1s3LiRNWvWsGKFBsMQ91OYORtgwYIFPP/SK/j2eDbnllhSuqTu/YWUbz/imwXzqFGjRqFtp7A/y+5MxetFDBgwgIULFxZ4u9Za7r73fnYle+LTsn+Bty/uIfnbj2lS04fx743N9Z59l2L16tUEBgYSEBBA2bJliY2NZdasWWctk5CQQLt2rgOlmJiYnPnGGFJSUkhLSyM1NZX09HRq1qyZr3hEnFBYORtc/3/69LsD325PUqZy3n4pLiVDxvE/SJz3BlM++5SQkJBC3VZhfpbdnYrXi2jdujVVqxb85fxXR45i4crV+HR6FGP0NpRGyesX4vPHRubNmkGZMmXy3d7evXupU6dOzvPatWuzd+/es5YJCwtjxowZAMycOZPExESOHDlCdHQ0MTEx1KpVi1q1atG5c2eaNGmS75hEilph5ezDhw/T4caueLe8m3K1ry3w9qX4O93F64VnnqJ79+6Fvr3C+iyXBKqaHDBz5kxGvfk2vj2ewaNsOafDEQec+m09af+dwpJFC84afa2wjR49mhUrVtC0aVNWrFiBv78/np6ebNu2jV9++YU9e/awd+9eli5dysqVK4ssLpHiLKeLV53r8Alu73Q44gCblUnS12/Rs2MbnvjH406HU+p5OR1AafPTTz9x5z33UfGm4XhVrOZ0OOKA9KN7SV7wJnO/+pLAwMACa9ff35/du3fnPN+zZw/+/v5nLXPVVVflnHlNSkriyy+/pHLlynzwwQe0aNECX19fALp06cIPP/xAq1atCiw+EXdkrWXAfQ+wM8kD3253OB2OOCR55ac0rFqWDyeMz3cXL8k/nXktQvv376dT1+6Ub/s3vGs1cDoccUDmqUQSZ4/gzTdG0rZt2wJtu3nz5mzdupWdO3eSlpbG1KlT6dmz51nLHD58mKysLABGjhzJvffeC7hGblmxYgUZGRmkp6ezYsUKdRsQAV4b9ToLvv2vuniVYskbFlFu/zrmz55ZIF28JP/0P7GInDp1ik5de2AbtsOn8Q1OhyMOsJkZJM9/g7v63Fooo7J4eXkxduzYnP6qt99+O0FBQQwfPpzZs2cDsHz5cho1akTDhg05ePAgzz77LAC9e/emfv36hISEEBYWRlhYGD169CjwGEXcyVdffcXI0f+Hb3d18SqtUn7fQNqPk1jy9QL1Py1GNEjBRfTt25fly5dz+PBhatasyUsvvcR9992Xpzastdzc+3ZWbTuK742P6ZJDKWStJWnJv2h6heXr+XPw9PR0OqRSR4MUlA4FkbMB4uPjuaFtOyre9DzetRoWQqRS3KX/uY9j04bx1RdT6dixY5Fvv6A+y+4qt5ytPq8XMWXKlHy38dzwF/h23S/43fqKCtdS6uS6OVRO3MnMZf9V4SpSiAoiZx84cMDVxavNgypcS6mslCSSZo/gjddecaRwhYL5LJdUKl4L2eTJUxgz/t9Ujn0D41XW6XDEAae2x5EZ/xVL1vyXihUrOh2OiOTi1KlTdOzSnawGMfg20Q8WSyOblUnS/NH07dWTwYMGOR2OnIeK10K0evVqHhz0d/xueQlPnypOhyMOSDu0i+RvxrBo/hzq1avndDgikgtrLX3vvJt9WX74tujjdDjikOTl/ya0ThXGjvk/p0ORC1DxWkh2795Nl+49qdBhMGVrXON0OOKAzORjJM5+jX+9+w7XX3+90+GIyEUMf/ElVsRtoqK6eJVayevmUvHPLcz65r94ealEKq70zhSCpKQk2nfuikdIdyoEXud0OOIAm5FO0rzXeei+u7nzTt0bUqS4mzp1Gm+Pe5/Ksf/Eo4y30+GIA07tWEv62i9ZsuZHKlWq5HQ4kgvdKquAZWVlcevtfTlS3p8KkTc7HY44wHVngXFcH1yfUa+NcDocEbmI1atXc//AQVTs+QyevuriVRqlH95N8qJ3mPPVlwQEBDgdjlyEitcC9o+hT7L61934thuoy06lVPKaL6mRcYjPp3yGh4f+i4kUZ7t376ZLj5uo0H4wZWuoaCmNMk8eJ3HOq7z7f29qVEE3oW4DBejDjz7iw8+mUanP6xhPjcJRGp3c8j0kfM3itWvw8fFxOhwRyUVSUhIdbuyGR3BXKjRQF6/SyGamkzzvDe6/M5Z77hngdDhyiVS8FpCVK1fyyONDqdT7VTwrqK9MaZR2cDsnl/6LFUu+wd/f3+lwRCQXWVlZ9O7TjyPeV+ET2cvpcMQB1lqSl4wnqvHVjH7jdafDkTzQNc0CsGPHDnr2uhWfzo9Splodp8MRB2QkHiFxzmt8NOF9IiIinA5HRC5i6FPDWL35d3zaq4tXaXUybiZXpO5j+rTJ6uLlZnTmNZ+OHz9O+05d8IroTflrmjkdjjggKz2FpLkj+ceQwdx++21OhyMiF/HxxxP54NMp6uJVip3c+l+yfp7H4rVr8PX1dTocySMVr/mQkZFBj169OVG1ERWbdnM6HHGAtVkkLxpDhxbhvPD8c06HIyIXsWrVKh5+7B/q4lWKpf2xg5OLx7L0m4XUqaOrpe5I58nzYfCQR9i45xi+be93OhRxSPIPU6njncJnn3ysS48ixdyOHTvocfMt+HR6RF28SqnMpD9JnP0aH4x/j6ioKKfDkcukM6+Xaey4cUydOc912cnD0+lwxAHJCSvw3P4tX6+Lo1y5ck6HIyK5OHHiBO07d8Wr2a2UD1C/9NIoKz2VpLkjeWTQg/TtG+t0OJIPKl4vwzfffMOw516g0u0j8SinvjKlUerezaSs/JDvv11OzZo1nQ5HRHKRkZFBz1tuI7FKQ3ybdXc6HHGAtZbkxWNpGxnEKy+96HQ4kk/qNpBHv/76K7379MW3yxOUqXKV0+GIAzJO/EHivFFM/nQiISEhTocjIhfx8COPs+H3o/ioi1eplfzjNK7yOMHk/3yiLl4lgM685sGRI0do37kLZVv0p9zVKlpKo6y0UyTNfo3hTz9Fjx49nA5HRC7ivX/9i8kzZlGpzxvq4lVKJf+yEo+ty1i0dg3ly5d3OhwpACpeL1FaWhpde/YipVYzfEM7OR2OOMBmZZK88P/o3v4Ghj7xD6fDEZGLWLx4MU8+87y6eJViqfu3cGrFBFYtX0qtWrWcDkcKiLoNXAJrLff/7SG2Hk3Hp9VdTocjDkle9R/qV/bgow/e12UnkWLu119/5dbbY/Ht8g918SqlMk4cJnHuKD6b+BHh4eFOhyMFSGdeL8HoN99i9jff4nfba7rslAcnt/xA8uZVBdegzeLk5pWuvz284C/1o29QO8oVwq+Iva9qTOqudZTbu5YF69ZQtmzZAt+GiBSco0eP0qFz1+wuXqFOhyMOyEpLIWnuazz9xGPcfPPNTocjBUzF60XMnTuXp595hsqdHybtjx3/m2GzSD+69+INWAvWYm0W2Kzs51nY7Omnp505354x3bVs9t+XoGLYjXhU8DtneubJ43iUrYDxKrrRZBLXzSXlt/WF03hWxjmTkjYsIrq2d4FuZv/ePWz/YQpeaYms/H4VV1xxRYG2LyIFKz09nW433UJy9WvxbdyarLSUQtiKPSMnZ+d4+F9OP/33JTBeZfHwrlAIMeaNzcq8tO+0PMpKScSmp2E8/1JueHhSpqp/gW/PeHphypYn6ev/o2ubFjwz7KkC34Y4T8XrRRw9foKgsKZwYBUc+N/0g/v2cnTfHucCu4DEuFlUrFz13OnHjp713LdSlUK/9F0GKHOeWPIjLTUFz7LlsX5X4kkWPllJLF20kGuvvbZAt3Pa5MmTufOuu/hq/nwaN25cKNsQkYKzb98+1q+LIzMri+SNS86al5ZyyqGoclfG+9z7RKennl10G2PwKluwB+e5ba8onO+7Kr9STp2iQpPWBFa0fPzvCeriVUIZe4lHhwCRkZE2Li6uEMORwrRnzx46d+5MQkICAImJiW45pnNCQgLvjB3HZ5MmQ42GlP1zJ2tX/0hAQECBb+vkyZP88ssvRETopubuzhiz1lob6XQcRUk5+2wnTpwgLS3tossZYy7pUZjS09PZvHkzTz/9NMuWLQNgxowZdO7cuVC3W9CysrLYt28fffv2Y926tQBs27aN+vXrF8r2rqxTD4Pl55/WUq1atULZhhSN3HK2itdS6OjRozmXv7Oystz2yDQ+Pp6YDp1I9atDw6pexMf91+mQpBhT8SruLDw8nPXr1/PDDz/QokULp8O5LKmpqTmjER46dKhQissJ//6QG66PLrSrcVJ0csvZuttAKVS1alWOHDkCwMiRIx2O5vKFh4czfdoUPBP3s33nb6xbt87pkERECkV8fDwA0dHRzgaSD97e3jlnvsMK6df/D95/nwrXUkDFaylVtWpV6tevz7PPPut0KPnSrl07rqpRjQyfarz1zrtOhyMiUmh27doFwPr1hfRD2CJQpkwZ3njjDfbt3XtJXThEzkfFaym2ZInrxwyZmZkOR3L5jDE89MC9eHl4sHjpMqfDEREpNHXr1gUgNjbW4Ujy5/HHHwdgwoQJDkci7krFayl2OhHu37/f4Ujy5+qrr8bLA5JOHHc6FBGRQlW2bFk2b97sdBj54unpul+6O59BFmepeBW3v3Tj5eWFMQableV0KCIihcrHx8fpEArMqVPF89ZlUvypeBURERERt6HiVURERETchopXEREREXEbKl5FRERExG2oeBURERERt6HiVURERETchopXEREREXEbKl5FRERExG2oeBURERERt6HiVURERETchopXEREREXEbKl5FRERExG2oeBURERERt6HiVdzaokWLWLJkCQAZ6Wl88sknbN++3eGoRETkfFJSUoiLiwNg586d/PLLLw5HJO5Ixau4tceeeJKx/xpP2dCuVGjQggEDBvDVrDlOhyUiIuexZMlSmjdvDsD333/P9Te0IjMz0+GoxN2oeBW39u2yJUQ2v47UTd+QvncTI19/g8cfe8TpsERE5Dy6du3CqyNHUb5iZa6oWYvFi77G09PT6bDEzXg5HYBIflxxxRWsXLaYl0a8SnRUc7p37+50SCIicgHGGJ4Z9hSdOrTH39+fWrVqOR2SuCEVr+L2ypYty6svv+R0GCIicokiIyOdDkHcmLoNiIiIiIjbUPEqRe7AgQPExsZSv359IiIi6Nq1K1u2bDnvsrt27SI4OBiA5cuXX3a3gLfffpuTJ0/maZ38bE9EpKRQzpbiRsWrFClrLb169aJt27Zs376dtWvXMnLkSA4ePFio272cRCgiUtopZ0txpOJVitSyZcsoU6YMAwcOzJkWFhbGDTfcwNChQwkODiYkJIRp06bl2k5ycjL33nsvUVFRNG3alFmzZgGQmZnJE088QXBwMKGhobz77ruMGTOGffv2ERMTQ0xMDOC6P2x0dDTNmjXjtttuIykpCYCFCxfSuHFjmjVrxowZMwrpVRARcQ/K2VIc6QdbUqQ2btxIRETEOdNnzJhBfHw869ev5/DhwzRv3pzWrVtfsJ1XX32Vdu3a8dFHH3Hs2DGioqLo0KEDn376Kbt27SI+Ph4vLy+OHj1K1apVeeutt1i2bBnVqlXj8OHDjBgxgsWLF+Pj48Prr7/OW2+9xZNPPskDDzzA0qVLCQwMpE+fPoX5UoiIFHvK2VIcqXiVYmHVqlX07dsXT09PatasSZs2bVizZg2hoaHnXX7RokXMnj2b0aNHA65RW37//XcWL17MwIED8fJyfbSrVq16zro//vgjCQkJtGzZEoC0tDSio6PZvHkz11xzDQ0aNADgjjvuYMKECYWxuyIibk05W5yk4lWKVFBQENOnT893O9ZavvzySxo1anRZ63bs2JEpU6acNT0+Pj7fcYmIlCTK2VIcqc+rFKl27dqRmpp61tHxhg0bqFy5MtOmTSMzM5NDhw7x7bffEhUVdcF2OnfuzLvvvou1FoCffvoJgI4dO/L++++TkZEBwNGjRwGoWLEiiYmJALRo0YLvvvuObdu2Aa6+WFu2bKFx48bs2rWL7du3A5yTKEVEShvlbCmOVLxKkTLGMHPmTBYvXkz9+vUJCgri6aefpl+/foSGhhIWFka7du144403uPLKKy/YzvPPP096ejqhoaEEBQXx/PPPA3D//fdz9dVX57Q1efJkAB588EFuvPFGYmJiqF69OhMnTqRv376EhobmXH4qV64cEyZMoFu3bjRr1owaNWoUyWsiIlJcKWdLcWROHwVdisjISBsXF1eI4UhRM8awfft2AgICnA5FpFAZY9Zaa0vVsD7K2SVP1apV+fPPP8nLd7eIO8otZ+vMq4iIiIi4DRWvIiIiIuI2VLyKiIiIiNtQ8SqOy8u42bnJz3CC8fHxzJ8//7LWFREpTZSzxWkqXsVRBTluthKhiEjhUs6W4kDFqzgqr+NmL1++nLZt29K7d28aN25M//79sdaedyzshx56iMjISIKCgnjhhRdy2l+zZg3XX389YWFhREVFcfz4cYYPH860adMIDw+/6BjdIiKllXK2FAvW2kt+REREWClZALt9+3bHtv/OO+/YRx999Jzp06dPtx06dLAZGRn2wIEDtk6dOnbfvn122bJl1s/Pz+7evdtmZmbaFi1a2JUrV1prra1bt649dOhQThtHjhyx1lqbkZFh27RpY9evX29TU1PtNddcY1evXm2ttfb48eM2PT3dfvzxx3bw4MFFsMfiFCDO5iHflYSHcnbJU6VKFev66naGcrYUldxyts68SrF0oXGzAaKioqhduzYeHh6Eh4eza9eu87bx+eef06xZM5o2bcqmTZtISEjg119/pVatWjRv3hwAPz+/nDG1RUTk8ihnS1FS8SqOCgoKYu3atXlax9vbO+dvT0/PnGEFz7Rz505Gjx7NkiVL2LBhA926dSMlJSXf8YqIlGbK2VIcqHgVRxXUuNlw9ljYJ06cwMfHh0qVKnHw4EEWLFgAQKNGjdi/f3/OGYHExEQyMjLOWldERM5POVuKAxWv4qiCGjcbzh4LOywsjKZNm9K4cWP69etHy5YtAShbtizTpk3j4YcfJiwsjI4dO5KSkkJMTAwJCQnq/C8ikgvlbCkOjM3D+MgaJ7vkMcawfft2AgICnA5FpFDlNk52SaWcXfJUrVqVP//8k7x8d4u4o9xyts68ioiIiIjbUPFaSm3cuJHRo0cD8O677/Ldd985HJGIiORmzZo1/PnnnwD88MMPDkcj4hwVr6XUipXfMXToUMA1ysn0r2Y7HJGIiFxIcnIyrVq3yXl+Q6tW/PHHHw5GlH/5GWFLSjcVr6XU4If+xscfT6SsdzlefmUEb70xyumQ8k2JUERKKh8fH5Ys/obKV1SnYuWqLJg/nxo1ajgdVr5cTs7OzMwspGjEnah4LcUGDLib1JRTPP/csxhjnA4n35QIRaQka9myJZs3/cwvGzfQqVMnp8PJ8c9//pMxY8YA8Nhjj9GuXTsAli5dSv/+/c877Ov5hoddtGgR0dHRNGvWjNtuu42kpCQA6tWrx1NPPUWzZs344osvHNhDKW5UvEqxo0QoInJ+NWvWxN/f3+kwztKqVStWrlwJQFxcHElJSaSnp7Ny5Upat27Nq6++SlxcHBs2bGDFihVs2LCBIUOGcNVVV7Fs2TKWLVvG4cOHGTFiBIsXL2bdunVERkby1ltv5WzjiiuuYN26dcTGxjq1m1KMqHiVYkeJUETEfURERLB27VpOnDiBt7c30dHRxMXFsXLlSlq1anXeYV//6scffyQhIYGWLVsSHh7OJ598wm+//ZYzv0+fPkW5S1LMaYBgKXb+mgibNWuWkwjHjBnD559/zoQJE8jIyGD//v0kJCQQGhp6VhtnJkKAtLQ0oqOjc+YrEYqIFIwyZcpwzTXXMHHiRK6//npCQ0NZtmwZ27Zto3z58owePZo1a9ZQpUoVBgwYcN5hX621dOzYkSlTppx3Gz4+PoW9G+JGdOZVip2/JsJWrVqdkwgvNv716UQYHx9PfHw8CQkJfPjhhznzlQhFRApOq1atGD16NK1bt6ZVq1aMHz+epk2bXnDYVzh7eNgWLVrw3XffsW3bNsB1d4UtW7Y4si9S/Kl4lWJJiVBExH20atWK/fv3Ex0dTc2aNSlXrhytWrW64LCvcPbwsNWrV2fixIn07duX0NBQoqOj2bx5s4N7JMWZug1IsdSqVSteffVVoqOj8fHxOW8irFOnznkT4em+r6cTYWpqKgAjRoygYcOGTu2SiEiJ1b59e9LT03Oen3myYOLEiedd5+GHH+bhhx/Oed6uXTvWrFlzznK7du0qsDilZDB5GR9Z42SLiLvKbZzskko5W0TcVW45W90GRERERMRtqHgVEREREbeh4lWKvYULF9KoUSMCAwMZNercYWxTU1Pp06cPgYGBXHfddTn9o7755hsiIiIICQkhIiKCpUuXFnHkIiKlw+XmaYCRI0cSGBhIo0aN+PrrrwHYvXs3MTExXHvttQQFBfHOO+8U1a6IG1DxKsVaZmYmgwcPZsGCBSQkJDBlypRzbnD94YcfUqVKFbZt28Zjjz3GU089BUC1atWYM2cOP//8M5988gl33nmnE7sgIlKi5SdPJyQkMHXqVDZt2sTChQsZNGgQmZmZeHl58eabb5KQkMCPP/7IuHHjzju4gZROKl6lWFu9ejWBgYEEBARQtmxZYmNjmTVr1lnLzJo1i7vvvhuA3r17s2TJEqy1NG3alKuuugqAoKAgTp06lXPnARERKRj5ydOzZs0iNjYWb29vrrnmGgIDA1m9ejW1atWiWbNmgOs2iE2aNGHv3r1Fvm9SPKl4lWJt79691KlTJ+d57dq1z0lgZy7j5eVFpUqVOHLkyFnLfPnllzRr1gxvb+/CD1pEpBTJT56+lHV37drFTz/9xHXXXVeIeyHuRMVrMRYTE5PT/+e0t99+m4ceeui8y9erV4/Dhw/n2uZrr7121vPrr78ecCWH4OBgAOLi4hgyZMjlhp1jwIABTJ8+Pd/t5NemTZt46qmneP/9950ORURE8iApKYlbb72Vt99+Gz8/P6fDkWJCxWsx1rdvX6ZOnXrWtKlTp9K3b9/LbvOvxev3339/zjKRkZGMGTPmsrdRkPz9/dm9e3fO8z179uDv73/BZTIyMjh+/DhXXHFFzvK9evXi008/pX79+kUXuIiUOqX1hEN+8nRu66anp3PrrbfSv39/brnllsvZJSmhVLwWY71792bevHmkpaUBrmS1b98+9u7dS0hICMHBwTmd3v/q5ptvJiIigqCgICZMmADAsGHDOHXqFOHh4fTv3x8AX1/fc9Zdvnw53bt3B6Br166Eh4cTHh5OpUqV+OSTT8jMzGTo0KE0b96c0NDQnDOa1lr+/ve/06hRIzp06MAff/yR79egefPmbN26lZ07d5KWlsbUqVPp2bPnWcv07NmTTz75BIDp06fTrl07jDEcO3aMbt26MWrUqLNG4hIRKQyl9YRDfvJ0z549mTp1KqmpqezcuZOtW7cSFRWFtZb77ruPJk2a8PjjjzuxW1KcWWsv+REREWGlaHXr1s1+9dVX1lprR44cae+55x5bp04d+8cff9j09HQbExNjZ86caa21tm7duvbQoUPWWmuPHDlirbX25MmTNigoyB4+fNhaa62Pj89Z7Z9+vnPnThsUFGSttXbZsmW2W7duZy0XFxdnQ0JC7LFjx+z7779vX3nlFWuttSkpKTYiIsLu2LHDfvnll7ZDhw42IyPD7t2711aqVMl+8cUX+X4N5s2bZxs0aGADAgLsiBEjrLXWPv/883bWrFnWWmtPnTple/fubevXr2+bN29ut2/fbq219pVXXrEVKlSwYWFhOY+DBw/mOx5xT0CczUO+KwkP5eyideTIEVu9enWbmppqrXXl1Tp16tjJkyfb4OBgGxQUZJ988smc5c/M2TfddJNt1qyZvfbaa+37779vrbX2qaeesh4eHjYsLMz269fPWnvxnN2lS5ecfOfn52cnTpxoMzIy7BNPPGEjIyNtSEiIHT9+vLXW2qysLDt48GDbsGFD2759e9ulS5fLztmXm6ettXbEiBE2ICDANmzY0M6fP99aa+3KlSstYENCQnL2Z968eZcVm7in3HK2EmEx99lnn9nY2FhrrbVhYWH2xRdftHfeeWfO/H//+9/2scces9aenQhfeOEFGxoaakNDQ62fn5/94YcfrLWXV7weOnTINmnSxP7888/WWmtvvfVW26BBg5yEUq9ePfv111/bRx55xH744Yc56/Xq1atAileRgqDiVYqCTjiIFIzccra6DRRzN910E0uWLGHdunWcPHmS8PDwi66zfPlyFi9ezA8//MD69etp2rQpKSkpl7X9zMxMYmNjGT58eE7/Kmst7777LvHx8cTHx7Nz5046dep0We2LiJQkZ3YdmDp1KnXr1qVt27ZUr14dLy8v+vfvz7fffnvOemPGjCEsLIwWLVqwe/dutm7detkxHD58mDvvvJPJkydTqVIlFi1axKeffkp4eDjXXXcdR44cYevWrXz77bf07dsXT09PrrrqKtq1a3fZ2xQpSipeizlfX19iYmK499576du3L1FRUaxYsYLDhw+TmZnJlClTaNOmzVnrHD9+nCpVqlChQgU2b97Mjz/+mDOvTJkypKenX/L2hw0bRmhoKLGxsTnTOnfuzL/+9a+cdrZs2UJycjKtW7dm2rRpZGZmsn//fpYtW5bPvRcRcS864SBS+FS8uoG+ffuyfv16+vbtS61atRg1ahQxMTGEhYURERHBTTfddNbyN954IxkZGTRp0oRhw4bRokWLnHkPPvggoaGhOT/YupjRo0ezaNGinB9tzZ49m/vvv59rr72WZs2aERwczN/+9jcyMjLo1asXDRo04Nprr+Wuu+4iOjq6QF8HEZHiTiccRAqfcXUruDSRkZE2Li6uEMMRESkcxpi11tpIp+MoSsrZzvjqq6/o1asXv/zyC40bN2bKlCm89tprWGvp1q0br7/+OuC6VVZcXBwVK1bk5ptvZteuXTRq1Ihjx47x4osv0rZtW5566ilmz55Ns2bNmDRpEr6+viQlJbFr1y66d+/Oxo0bWb58OaNHj2bu3LkYYwgKCsLLywuAl19+me7du/Pcc88xZ84crLVUr16dr776Cj8/Px5++GG++eYbrr76asqUKcO9995L7969nXz5RIDcc7aKVxEpFVS8ioi4j9xytroNiIiIiIjbUPEqIiIiIm5DxWsptnDhQho1akRgYCCjRo1yOhwREbkI5W0RFa+lVmZmJoMHD2bBggUkJCQwZcoUEhISnA5LREQuQHlbxEXFaym1evVqAgMDCQgIoGzZssTGxjJr1iynwxIRkQtQ3hZxUfFaSu3du5c6derkPK9duzZ79+51MCIREcmN8raIi4pXEREREXEbKl5LKX9/f3bv3p3zfM+ePfj7+zsYkYiI5EZ5W8RFxWsp1bx5c7Zu3crOnTtJS0tj6tSp9OzZ0+mwRETkApS3RVy8nA5AnOHl5cXYsWPp3LkzmZmZ3HvvvQQFBTkdloiIXIDytoiLhocVkVJBw8OKiLgPDQ8rIiIiIiWCilcRERERcRsqXkVERETEbah4FRERERG3oeJVRERERNyGilcRERERcRsqXkVERETEbah4FRERERG3oeJVRERERNyGilcRERERcRsqXkVERETEbah4FRERERG3oeJVRERERNyGilcRERERcRsqXkVERETEbah4FRERERG3oeJVRERERNyGilcRERERcRvGWnvpCxtzCPit8MIRESk0da211Z0OoigpZ4uIG7tgzs5T8SoiIiIi4iR1GxARERERt6HiVURERETchopXEREREXEbKl5FRERExG2oeBURERERt6HiVURERETchopXEREREXEbKl5FRERExG2oeBURERERt/H/g9nS7xFGcvIAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "Composting.plot_sankey_report(margin=1, offset=.4, gap=0.8)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Sensitivity to precipitation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(0, 0.5, 'Percent of Incoming PFAS that \\n remains in the Compost (%)')"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEGCAYAAABCa2PoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAA4cUlEQVR4nO3dd3gVZfr/8fcnCb0KAaQXQRSpgnQExYIdbCuuir2i6Opafu7q2t21IrZV194LAmJBpAooEqSGGhFpIhGUKjX374+Z+D3GlENykpNyv67rXJkz88zMfQ5cufPM02RmOOeccwWVEO8AnHPOlQ6eUJxzzsWEJxTnnHMx4QnFOedcTHhCcc45FxNJ8Q4gXpKTk61Zs2bxDsM550qU2bNn/2xmdbI7VmYTSrNmzUhJSYl3GM45V6JI+iGnY/7IyznnXEx4QnHOORcTnlCcc87FhCcU55xzMeEJxTnnXEwUm4Qi6UVJGyQtzOG4JD0hKU3SfEmHRxwbIml5+BpSdFE755zLVGwSCvAyMCCX4ycArcLX5cAzAJJqAXcC3YCuwJ2SDijUSJ1zzv1JsUkoZjYV2JRLkdOAVy3wNVBTUn3geGC8mW0ys1+A8eSemAokbcM2Hvl8KfsyfNp/55yLVGwSShQaAqsj3q8J9+W0/08kXS4pRVJKenp6voL4dMGPjJiYxuWvprB15558XcM550qjkpRQCszMnjOzLmbWpU6dbGcOyNO1/Vtxz2mHMXlZOmc8M4NVG3fEOErnnCuZSlJCWQs0jnjfKNyX0/5Cc36PZrx6cVd+2rKL056axswVGwvzds45VyKUpIQyBrgg7O3VHdhsZj8C44DjJB0QNsYfF+4rVL1aJjPqml4cUKU8f31hJm9/s6qwb+mcc8VasUkokt4CvgJaS1oj6RJJV0q6MizyCbACSAOeB64GMLNNwD3ArPB1d7iv0DVPrsKHV/eiZ8tkbh25gLs+SmXvvoyiuLVzzhU7MiubvZW6dOlisZpteO++DO77ZDEvTV/JkQfXYcTgTtSoVC4m13bOueJE0mwz65LdsWJTQynJkhITuPOUw3jg9HbMSPuZQU9P5/uft8c7LOecK1KeUGJocNcmvH5pN37ZvpuBT01netrP8Q7JOeeKjCeUGOveojZjhvamXvUKXPDiN7z61UrK6mNF51zZ4gmlEDSuVZkPrurJUa3rcMfoVG4buYDde72x3jlXunlCKSTVKpbjufO7MPSolrw9azXnPv816Vt3xTss55wrNJ5QClFCgrjp+NY8eW4nUtdt4dQnpzF/za/xDss55wqFJ5QicHL7Brx/VQ8SJM569itGzSnUgfzOORcXnlCKyGENajBmaC86NK7J9e/M5YFPFvuMxc65UsUTShGqXbUCb1zajfO7N+W/U1dw8cuz2LzDZyx2zpUOnlCKWLnEBO4Z2Jb7B7Vjxnc/M/Dp6aRt2BbvsJxzrsA8ocTJud2a8OZl3dm6cw+DnprOhMU/xTsk55wrEE8ocXREs1qMHtqbpsmVufTVFEZMWE6Gt6s450ooTyhx1rBmJd67oiendWjAI+OXceXrs30lSOdcieQJpRioVD6Rx/7SkTtObsOEJRsY+NR0vkv3dhXnXMniCaWYkMTFvZvz+iXd+GXHHgY+OZ0vFnm7inOu5PCEUsz0OKg2H13bm2bJVbj01RQe/2KZt6s450oETyjFUMOalXjvyh6cfnhDHv9iOZe/Npst3q7inCvmPKEUUxXLJfLIWR2469TDmLw0aFfx8SrOueLME0oxJokhPZvxxqXd2PLbHgY+NZ1xqevjHZZzzmXLE0oJ0K1F0K5yUJ0qXPHabB75fKnPA+acK3Y8oZQQ9WtU4p0renB2l0aMmJjGhS99wy/bd8c7LOec+12eCUXSv6PZ5wpfxXKJ/OfMDjxwejtmrtjEySOmsWDN5niH5ZxzQHQ1lGOz2XdCrANx0RvctQnvXdkDM+OMZ2fw7qzV8Q7JOedyTiiSrpK0AGgtaX7E63tgftGF6LLToXFNxl7Xh67NanHzB/O5beR8du7ZF++wnHNlWFIux94EPgUeAG6N2L/VzDZld4KkisDJQB+gAfAbsBD42MxSYxKx+12tKuV55eKuPPL5Up6e/B2p67bwzHmdaVizUrxDc86VQTKLrreQpLpAxcz3ZrYqy/G7CJLJZGA2sCEsfzBwVLh9o5kVi9pNly5dLCUlJd5hxMy41PXc9O48khLFiMGH07tVcrxDcs6VQpJmm1mXbI/llVAknQI8SlDj2AA0BRab2WFZyp1kZh/ncp26QBMzKxa/xUtbQgFYkb6NK1+fTdqGbdx4XGuu6nsQCQmKd1jOuVIkt4QSTaP8vUB3YJmZNQf6A19nLZRdMpFUUVL18PiGvJKJpAGSlkpKk3RrNsebSpoQtuVMltQo4th/JKVKWizpCUll7jdpizpV+fDqXpzYrj4PjVvKFa/PZvNvPmWLc65oRJNQ9pjZRiBBUoKZTQKyzU6RJF0KjAI+kPRAFOUTgacIepC1AQZLapOl2MPAq2bWHriboH0HST2BXkB7oC1wBNA3is9W6lSpkMSIwZ244+Q2TFqygVOfnEbqOu9a7JwrfNEklF8lVQWmAm9IGg5sz1pI0qlZdh1jZgPM7FjgxCju0xVIM7MVZrYbeBs4LUuZNsDEcHtSxHEjaKMpD1QAygFldu73zKnw37miO7v2ZDDo6Rm8M2tV3ic651wBRJNQTiPorXUD8BnwHXBKNuXaSRotqWP4fr6kFyQ9D0TTw6shEDmgYk24L9I84PRwexBQTVJtM/uKIMH8GL7GmdnirDeQdLmkFEkp6enpUYRUsnVuWoux1/XmiGYHcMsHC/j7e/P4bbd3LXbOFY7cug0DYGaRtZFXcil3n6QDgbvD9ot/AtWASjHs2XUT8KSkCwlqTGuBfZJaAocCmW0q4yX1MbMvs8T4HPAcBI3yMYqpWEuuWoFXL+7G8C+W8cTENBas3cwz53WmeXKVeIfmnCtlopl65XRJyyVtlrRF0lZJW3Iovh24HniS4Bf3YGBZlLGsBRpHvG8U7vudma0zs9PNrBNwe7jvV4Laytdmts3MthGMn+kR5X1LvcQE8bfjWvPSRUewfstOThkxjc8W/hjvsJxzpUw0j7z+A5xqZjXMrLqZVTOz6lkLSboX+AAYCxxlZqcCc4FPJF0QxX1mAa0kNZdUHjgHGJPlHsmSMmO+DXgx3F4F9JWUJKkcQYP8nx55lXVHta7Lx9f14aC6Vbny9W+5d+wi9uzLiHdYzrlSIpqE8lN27RHZONnMjiPoVnwBgJmNAY4DDsjrZDPbCwwFxhEkg3fNLFXS3REN/v2ApZKWAfWA+8L97xO07SwgaGeZZ2YfRRFzmdOwZiXevaI7F/RoygvTvmfwc1+zfvPOeIflnCsFchzYKCmz8bsvcCBBF+BdmcfNbGSW8q8DO4BKwM9mdkMhxBszpXFg4/4aPXctt41cQKVyiQw/p5OPrnfO5SlfI+UlvZTLNc3MLs7mnHYE41aW5CvSIuQJJZC2YStXvf4taenbuPboVgzr34pEH13vnMtBbgklx15eZnZReHIvM5ue5YK9srlJbzOblksQ1QmmXlkYdeSu0LWsW43RQ3vxj1ELeWLCclJWbuLxczpSt1rFvE92zrkI0bShjIhy3xmSZki6Q9JJkrpKOlLSxZJeI2is92lwi6HK5ZN49OyO/OfM9ny76hdOHD6NGWk/xzss51wJk2MNRVIPoCdQR9LfIg5VBxKzljezGyTVAs4AzgLqEwyIXAz8N7faiysezu7SmA6NanL1G7M5738zGdb/YIYe3dIfgTnnopLbwMbyQNWwTLWI/VuAM7M7IVwn5fnw5Uqg1gdWY8zQ3tz+4QIe+2IZs8JHYMlVK8Q7NOdcMRfN9PVNzeyHIoqnyHijfO7MjHdmrebOManUqFSOJwZ3onuL2vEOyzkXZwWavr40JhOXN0mc07UJo67pRdUKSZz7/Nc8NSmNjIwyMWONcy4fommU3y+S/vRsJLt9rmQ4tH51xlzbm5PaN+ChcUu58OVZ/LxtV94nOufKnJgnFOCrKPe5EqJqhSSeOKcj9w1qy9crNnLi8C+9F5hz7k9yTCiSLpPUKtyWpJfCySHnSzo8m/IHSuoMVJLUSdLh4asfULmwPoArGpL4a7emjL6mF9UqJvHX/83k0c+XstfnAnPOhXLr5TUMeDncHkywGmJzoBMwHOiTpfzxwIUEswQ/AmT2Nd0K/L+YROvi7tD61fno2t7cMTqVJyam8fWKTQwf3JH6NXyIkXNlXW5Tr8w1s47h9pvATDMbHr7/1sz+VEsJj51hZh8UUrwx4728Cu7DOWu4/cOFVEhK4OGzOtD/0HrxDsk5V8jy28srQ1J9SRUJZhD+IuJYbn+ONpJUPXxM9oKkbyUdl4+4XTE3qFMjxl7bm/o1KnHJKyncM3YRu/f6IzDnyqrcEsodQAqwEhhjZqkAkvoCK3I572Iz20IwbX1t4HzgwZhE64qdFnWqMvLqngzp0ZT/TfueM5+dwQ8bt+d9onOu1MkxoZjZWKApcKiZXRZxaBbB4lc5yWw7ORF4NUxEPndHKVaxXCJ3ndaWZ8/rzMqft3PSE9P4aN66eIflnCtiuXYbNrO9ZvYL/N7Tqz/BxJBpuZw2W9LnBAllnKRqgD8HKQMGtD2QT4b14eB6Vbn2rTnc+sF8duzeG++wnHNFJJo15btLegL4ARgNTAUOyeWUS4BbgSPMbAfBnGAXxSBWVwI0OqAy71zRg6v6HcQ7Kas5ZcQ0Fq3bEu+wnHNFILdxKPdLWk6wzO58gu7C6Wb2SmatJTtmlkHQdfgfkh4GeprZ/BjH7YqxcokJ3DLgEF6/pBtbd+5l4FPTeWn69+Q1b5xzrmTLrYZyKfAT8AzwmpltBPL8jSDpQYIxLIvC13WS7o9BrK6E6dUymc+uP5I+rZK566NFXPJKCht92hbnSq3cxqEkAscSDGrsD0wCjgEam1mOD8YlzQc6hjWVzOvMMbP2MY69QHwcStExM1796gfu+2QxNSuV49GzO/r69c6VUPkah2Jm+8zsMzMbAhwEjAKmA2vDgY65qRmxXWP/wnWljSSG9GzGqKuDaVvOf3EmD366hD0+bYtzpUqubSgRb480sw/M7EygFfBZLtd8AJgj6WVJrwCzCdphXBnXpkF1xl7bh3OOaMKzU77jzGd8zIpzpUluj7x+n14lt6lWcji3PnAEQZvLLDNbH4tgY8kfecXXpwt+5JYP5pNhcM/AwxjUqVG8Q3LORaFAC2zlUw+gX/jqUUj3cCXYCe3q8+n1R3Jo/Wrc8M48hr09hy0798Q7LOdcAeQ223BdSX8jGOWeuf07M3s0u5MkPQ20BN4Kd10h6RgzuyYWAbvSo2HNSrx1WXeenvwdwycsJ2XlLzx+TkeOaFYr3qE55/Iht0ded+Z2opndlcN5Swima7HwfQKQamaHFjDWmPJHXsXLt6t+4fq357Lmlx0MPaol1/ZvRbnEwqpAO+fyK7dHXjnWUDIThqRkM9uf5fnSgCYEI+sBGpP7VC2ZQQ4gWGclEXjBzB7Mcrwp8CJQB9gEnGdma8JjTYAXwnsZcKKZrdyPmF2cHd7kAD4Z1oc7w3VWpi7/meHndKRp7SrxDs05F6XcenmdLCkdmC9pjaSeUV6zGrBY0mRJkwkGN1aXNEbSmBzulQg8BZwAtAEGS2qTpdjDBJNNtgfuJuhNlulV4KGwFtQV2BBlrK4YqVohiUfO7sCIwZ1Ykb6NE4d/yXspq32EvXMlRG5tKPcDfcxsiaRuwH+AvlFc8458xNEVSDOzFQCS3gZOI0hGmdoAme04kwjGxRAmniQzGw9gZtvycX9XjJzSoQGHNz2Av70zl7+/P5/Jy9K5f2A7alQuF+/QnHO5yO0h9V4zWwJgZjMJah55MrMpZjYFmAMsyHxF7M9OQ2B1xPs14b5I84DTw+1BQDVJtYGDgV8ljZQ0R9JDYY3nTyRdLilFUkp6eno0H8fFScOalXjzsu7cPKA14xauZ8DwqXz13cZ4h+Wcy0VuCaWupL9lvrJ5n63wl/Z6ggklUwgGNsai9fsmoK+kOQQ1pbXAPoJaVp/w+BFAC4K17f/EzJ4zsy5m1qVOnToxCMkVpsQEcXW/loy8uicVyyVy7gtf8+CnS9i1d1+8Q3POZSO3hPI8Qa0k85X1fU7+DrQ1s2Zm1sLMmptZizziWEvQoJ6pUbjvd2a2zsxON7NOwO3hvl8JajNzzWxFOMfYKCDqQZiu+GvfqCZjr+3NOUc05tkp3zHwqRks+2lrvMNyzmWRZy+vfPgO2LGf58wCWklqTpBIzgHOjSwgKRnYFE46eRtBj6/Mc2tKqmNm6cDRxKZG5IqRKhWSeOD09hx9SD1u/WA+J4+Yxi0DDuGins1ISPAFQZ0rDnJrlM+v24AZkmYCv89VbmbX5XSCme2VNBQYR9Bt+EUzS5V0N5BiZmMIRt0/IMkIFvm6Jjx3n6SbgAmSRPCI7flC+FyuGDi2TT06Nj6SWz+Yzz1jFzFxyU88fFYH6teoFO/QnCvzchzYmO8LSt8A0wga43+fTtbMXonpjQrIBzaWbGbG27NWc/dHiyiXKO4d1I5TOzSId1jOlXr5GthYAOXMLMdGe+diQRKDuzahR4vaXP/OXK57aw4TFv/E3ae1pUYl717sXDzkmVBy6NG1GZhtZnOzOfappMuBj/jjI69N+Q3SuZw0S67C+1f2+H0+sG++38QjZ3WgZ0tfwMu5opbnI69wMa0uBAkC4GSCLsHNgPfM7D9Zyn+fzWUsip5eRcofeZU+81b/yg3vzGXFz9u5tHdzbjq+NRXLZTskyTmXT7k98oomoUwlmBtrW/i+KvAxMICglpJ1ipQSwRNK6bRj917u/2Qxr3+9ilZ1q/Lo2R1p18gXDXUuVgq6HkpdIh5dAXuAemb2W5b9mTcrJ+k6Se+Hr6GS/KG2KxKVyydx78B2vHzREWzZuYdBT09n+BfLfblh54pANAnlDWCmpDvDKe2nA29KqsIf59rK9AzQGXg6fHUO9zlXZPq1rsu464/kxHb1eeyLZZz5zAzSNvg0b84Vpqi6DUs6AsicbXi6meX4rEjSPDPrkNe+ePNHXmXHx/N/5PZRC/ht9z5uPeEQhvTwwZDO5VcslgD+FngP+BDYEK4/kpN9kg6KuHkLgjm3nIuLk9rX5/Prj6RXy2Tu+mgRf31hJmt//S3eYTlX6kTTKH8tcCfwE0FiEEGvrfY5lO8PvASsCMs2BS4ys0kxjLvAvIZS9pgZ78xazT1jF5EgcccpbTizcyOCCRacc9EoaC+vNKCbmUU9d7ikCkDr8O1SM/tT4328eUIpu1Zv2sGN787jm5WbOLZNPe4f1I461SrEOyznSoSCPvJaTTCQMa+bnCfpfAAz22Vm881sPnC2pHPzON25ItO4VmXeurw7t594KFOWpXPcY1MYO39dvMNyrsSLpobyP4Laxsf8ceT7o1nKzQT6Z10xMewNNtXMOscq6FjwGooDWP7TVm58bx7z12zmpPb1uee0ttSqUj7eYTlXbBW0hrIKGA+UJ/f1UMplt/yumW0HfByKK5Za1avGyKt68vfjW/N56nqOe2wK41LXxzss50qkPOfy2o91USpJqhImkN9JqkaQjJwrlpISE7jmqJYcfUhdbnx3Hle8NpuBHRvwr1MPo2Zl/6/rXLRyrKFIejz8+ZGkMVlf2ZzyP+B9SU0jrtEMeDs85lyxdmj96owe2ovrj2nF2Pk/ctxjU5m45Kd4h+VciZFbDeW18OfD0VzIzB6WtA2YGs73BbANeNDMfKS8KxHKJSZw/TEHc8yh9bjpvXlc/HIKZ3VuxD9PaUP1iv7k1rncxHyBLfj9MRdmVmwX/vZGeZeXXXv3MWJCGs9M+Y661Srw4Bnt6XtwnXiH5VxcFahRXtLJkuZI2iRpi6Stkrbkdo6ZbS3OycS5aFRISuSm41sz8qqeVKmQxJAXv+GW9+ezZeeeeIfmXLEUTS+vx4EhQG0zq25m1cyseuGG5Vzx0aFxTcZe25ur+h3Ee7NXc/xjU5m0dEO8w3Ku2Il2YONCK4xnY86VEBXLJXLLgEP48OpeVKuYxEUvzeLGd+exeYfXVpzLFM2a8jcDn0iaQi4DGzNJqgzcCDQxs8sktQJam9nYWATsXDx1aFyTj67tzZMT03h68nd8uTyd+we145g29eIdmnNxF00N5T5gB1CR3Ac2ZnqJIPH0CN+vBe4tQIzOFSsVkhK58bjWjL6mF7WqlOfSV1O4/u05/LJ9d7xDcy6uoqmhNDCztvtxzYPM7C+SBgOY2Q75dK6uFGrbsAZjhvbmqUlpPDUpjWlpG7lvUFuOP+zAeIfmXFxEU0P5RNJx+3HN3ZIqAQYQro1S7GYbdi4WyiclcMOxBzN6aC/qVqvAFa/N5tq35rDJayuuDIomoVwFfCZpZ9hlOK9uw3cCnwGNJb0BTCBoh3Gu1DqsQQ1GD+3FjccezGcLf+SYR6cwZt46vC+LK0sKa2BjbaA7wQJbX5vZzzG/SQH5wEZXWJau38rNH8xn3upfOebQutw7sB0H1qgY77Cci4kCLwEs6VRJD4evk6M4pSLwC7AFaCPpyOjDda5ka31gMIPxP046lGlpP3Pso1N465tVXltxpV40I+UfBIYBi8LXMEkP5FL+38B04Hbg7+HrpmiCkTRA0lJJaZJuzeZ4U0kTJM2XNFlSoyzHq0taI+nJaO7nXGFJTBCX9mnBuOuPpG3DGtw2cgHnPj+THzZuz/tk50qoaBbYmg90NLOM8H0iMCeXNeWXAu33d9nf8LrLgGOBNcAsYLCZLYoo8x4w1sxekXQ0wVr150ccHw7UATaZ2dDc7uePvFxRMTPenrWa+z9ezJ6MDG46rjUX9WpOYoJ3fnQlT4EfeQE1I7Zr5FF2BflbUKsrkGZmK8xsN8G096dlKdMGmBhuT4o8LqkzUA/4PB/3dq7QSGJw1yaM/1tferdM5t6PF3P6MzNYut6nu3OlSzQJ5QFgjqSXJb0CzCYY7PgHkkZIeoJgEORcSf+V9ETmK4r7NCSY5iXTmnBfpHnA6eH2IKCapNqSEoBHyOPRmqTLJaVISklPT48iJOdi58AaFXn+gi48MbgTqzft4OQRX/LY+GXs2rsv3qE5FxPRrNj4lqTJwBHhrlvMLLs1UjOfH80Gsi7AFavWyJuAJyVdCEwlGIW/D7ga+MTM1uQ2htLMngOeg+CRV4xici5qkji1QwN6t0zmro9SGT5hOR8v+JF/n9GOzk1rxTs85wokz4QiaRAw0czGhO9rShpoZqMiy5nZK+HxYWY2PMs1hkURy1qgccT7RuG+yHusI6yhhIt4nWFmv0rqAfSRdDVQFSgvaZuZ/alh37nioFaV8gw/pxMDOzbkH6MWcuazX3Fet6bcPKA11XwhL1dCRdMoP9fMOmbZN8fMOuVQ/lszOzza8hFlkgga5fsTJJJZwLlmlhpRJpmgwT1D0n3APjO7I8t1LgS6eKO8Kym279rLw58v5eUZK6lXrSL3DGzLsT7ZpCumCtoon12ZP9VsJA2W9BHQPMv685OATXndxMz2AkOBccBi4F0zS5V0t6RTw2L9gKWSlhE0wP+pLce5kqZKhSTuPOUwRl7Vk5qVy3HZqylc88a3bNi6M96hObdfoqmhvAj8CjwV7roGqGVmF2Yp1xRoTtCIH/moaSswP0wYxYbXUFxxtGdfBs9NXcHwCcupmJTA7ScdytldGuPzq7riIrcaSjQJpQrwT+CYcNd44F4zK9EjtDyhuOLsu/Rt3DZyAd98v4nuLWrxwOntaZ5cJd5hOVewhFJaeUJxxV1GhvFOymru/2Qxu/ZmMKx/Ky7r04LySdEOH3Mu9gpaQzmYoLtuMyLaTszs6BjGWOQ8obiSYsOWndw5JpVPF67n4HpVuX9QO7o08y7GLj4KmlDmAc8SjC/5fQSWmc2OZZBFzROKK2m+WPQTd4xeyLrNOzm3WxNuGXAINSp5F2NXtHJLKNGs2LjXzJ7Zj5v1Av4FNA2vL8DMrEW013DO/dkxberR46DaPDp+GS9N/57xi37izlPacFK7+t5o74qFaGoo/wI2AB8SsfKimWXbFVjSEuAG/lyj2VjwcGPHayiuJFu4djO3jVzAgrWbOap1He4+rS2Na1WOd1iuDCjoI6/vs9mdY41D0kwz67b/YRYtTyiupNu7L4NXvvqBRz5fihnccGwrLu7VnKREb7R3hadIe3mF66ckAiP5Y43m25jeqIA8objSYt2vv3HH6FS+WPwTh9avzgOnt6Nj45rxDsuVUvlKKJJOz/ZAyMxG5nDepOyLF69eYZ5QXGliZoxLXc+dY1LZsHUX53Vryk3Ht/ZGexdz+W2UPyWXY0ZQA/nzAbOj9iM251wMSGJA2/r0apnMI58v49WvVvJZ6nr+eXIbTmnvjfauaMTskZek88zsdUl/y+64mT0akxvFiNdQXGm2YM1mbh+1gPlrNtOnVTJ3n9bWR9q7mIjFio3RyPzfWi2Hl3OuiLRrVIMPr+7F3acdxtxVv3L841MZ/sVyX8zLFSqfesW5Um7Dlp3c8/FiPpq3jhbJVbhnYFt6tUyOd1iuhMpXDUXSWeHP5oUVmHOu8NWtXpERgzvx2iVdyTDjry/MZNjbc3x6fBdzuT3yui38+UFRBOKcK1x9WtXhs+uPZFj/Vny6YD39H5nCa1+tZF9G2XxK4WIvt27D4wl6cx0BfJn1uJmd+qeTShB/5OXKshXp2/jn6IVMT9tIu4Y1uGdgWx+74qKS33Eo5YHDgdeAS7MeN7MpOZw3DHiJYGGtF4BOwK1m9nm+oi8knlBcWWdmjJ3/I/eMXUT6tl2cc0QTbhnQmpqVy8c7NFeMFXTqlTpmli6pKoCZbcuj/Dwz6yDpeOAKgsW5Xsu6zny8eUJxLrB15x4e/2I5L89YSY1K5bh1wCGc2bkRCQk+dsX9WUG7DdeTNAdIBRZJmi2pbW73C3+eSJBIUiP2OeeKmWoVy/HPk9sw9tretEiuws0fzOes/37FonVb4h2aK2GiSSjPAX8zs6Zm1gS4MdyXk9mSPidIKOMkVQMyCh6qc64wHVq/Ou9e0YOHzmzPyp+3c/KIL7nro1S27twT79BcCRHVAltm1iGvfRHHEoCOwAoz+1VSbaChmc2PUcwx4Y+8nMvZ5h17eOjzJbwxcxXJVSvwj5MO5dQODXwKF1fgR14rJP1TUrPw9Q9gRU6FzSwD+AloI+lI4DCgZj7ids7FSY3K5bh3YDtGXd2L+jUqMuztuZzz3NcsXb813qG5YiyaGsoBwF1Ab4JuxF8Cd5nZLzmU/zfwF2AR/7fAlhW3bsZeQ3EuOvsyjLdnreKhcUvZunMvF/Royg3HHkz1ij6TcVlU1OuhLAXam9muPAvHkScU5/bPL9t389DnS3nrm1XUrlKB2044hEGdGnpvsDKmqCaHzLQC8D9dnCtlDqhSnvsHtWP0Nb1odEAlbnxvHmf99ytS122Od2iumCiMGsoHQAdgAn9csfG6mN6ogLyG4lz+ZWQY789ew4OfLeHXHbs5r3tTbjy2NTUq+9+SpV2BaiiSekWzL8IY4B5gBjA74pXXfQZIWiopTdKt2RxvKmmCpPmSJktqFO7vKOkrSanhsb/kdS/nXMEkJIizj2jMpBv7cX73prz+9Q8c9chk3pm1igyfG6zMiqZR/tuso9yz21egIKREYBlwLLAGmAUMNrNFEWXeA8aa2SuSjgYuMrPzJR1M0Oi/XFIDguR1qJn9mts9vYbiXOykrtvMnaNTSfnhFzo0qsG/Tj2MTk0OiHdYrhDkawlgST2AnkCdLKswVgcSsyn/rpmdLWkBQW+wPzCz9rnE2BVIM7MV4bXeBk4j6CmWqQ2QGcckYFR43WUR91gnaQNQB/g1l/s552LosAY1eO/KHnw4Zy0PfrqEQU/P4MzOjbh5QGvqVqsY7/BcEcltTfnyQNWwTOSKi1uAM7MpPyz8eXI+4mgIrI54vwbolqXMPOB0YDgwCKgmqbaZbcwsIKlrGPd3+YjBOVcAkjj98EYcd9iBjJi4nBenfc9nC9czrH8rhvRsRvmkwugD5IqTaB55NTWzHwo1COlMYICZXRq+Px/oZmZDI8o0AJ4EmgNTgTOAtpmPtiTVByYDQ8zs6xzuczlwOUCTJk06//BDoX4s58q0FenbuGfsIiYtTadFnSrcecph9D24TrzDcgVU0G7DFSQ9J+lzSRMzX7ncrLukWZK2SdotaZ+kvGaZWws0jnjfKNz3OzNbZ2anm1kn4PZw36/hPasDHwO355RMwvLPmVkXM+tSp47/x3auMLWoU5WXLurKixd2ISPDGPLiN1z2agqrNu6Id2iukOT2yCvTe8CzBGub7MujLAS1iHPC87oAFwAH53HOLKBVuNzw2vD8cyMLSEoGNoVTu9wGvBjuLw98CLxqZu9HEZ9zrggdfUg9erVM5sVpKxkxcTnHPDaFy/u04OqjDqJy+Wh+BbmSIpoayl4ze8bMvjGz2Zmv3E4wszQg0cz2mdlLwIA8yu8FhgLjgMXAu2aWKuluSZlTtvQDlkpaBtQD7gv3nw0cCVwoaW746hjF53LOFZEKSYlc1e8gJt3Uj5Pa1efJSWkc/fAURs9dS6zHwrn4iaYN5V/ABoJaQORAxU05lJ8KHENQo1kP/AhcmNPsxPHi3Yadi5+UlZv410epLFy7hc5ND+COk9vQwZcgLhEKumLj99nsNjNrkUP5pgQJqBxwA1ADeDqstRQbnlCci6/M0fb/GbeUn7ftCroZH9+autW9m3FxVqSTQ5YUnlCcKx627tzDk5PSeHHa95RPTGDo0a24uHczKiT9abibKwYKOvVKZUn/kPRc+L6VpBzHmkg6WdIcSZskbZG0NYpeXs65MqpaxXLcdsKhjL+hLz1bJvPvz5Zw7KNTGZe63ttXSphoGuVfAnYTjJqHoBfWvbmUfxwYAtQ2s+pmVs3MqhcoSudcqdcsuQrPX9CF1y7pSoWkBK54bTbn/W+mL+pVgkSTUA4ys/8AewDMbAeQ2wIIq4GF5n9aOOfyoU+rOnw6rA93nXoYC9du4YThU/nHqAVs3Fasl1hyRDcOZbekSoTzc0k6iIjeXtm4GfhE0hT+2Cvs0YIE6pwrO5ISExjSsxmndmjA418s4/WZqxg9dx3XHt2SIT29faW4iqaGcifwGdBY0hsE65zcnEv5+4AdQEWCOcAyX845t18OqFKeu05ry7jr+9Cl6QHc/8kSjntsKp8t9PaV4iiqXl6SagPdCR51fW1mP+dSdqGZtY1diIXDe3k5V/JMWZbOvWMXsXzDNrq3qMU/TmpD24Y14h1WmVLQXl6DCEbLf2xmY4G9kgbmcsonko7LX6jOOZezvgcH7Sv3DGzL0vVbOeXJadz8/jw2bNkZ79Ac0Q1snGtmHbPsmxNO0phd+a1AFYKeYbsJajVW3Hp6eQ3FuZJt8297eHLicl6esZJyiQlcc1RLLundnIrlvH2lMBV0tuHsyuTYmB92E04ws4rebdg5V1hqVCrH7Se1YfwNfendMpmHxi2l/yPB/GC+DHF8RJNQUiQ9Kumg8PUouawRr8B5kv4Zvm8cLnzlnHMx1yy5Cs9d0IU3L+tGzcrlGPb2XAY9PZ1ZK7OdbtAVomgSyrUEj67eAd4GdgLX5FL+aaAH/zf9/DbgqQLE6Jxzeep5UDIfDe3Nw2d14Kctuzjr2a+48rXZrPx5e7xDKzNyHYciKREYa2ZH7cc1u5nZ4ZLmAJjZL+GaJc45V6gSEsSZnRtxUrv6PP/lCp6d8h0TlvzE+d2bcV3/ltSs7L+KClOuNRQz2wdkSNqffnl7wkSUORCyDpCR/xCdc27/VCqfyHX9WzH5pn6c2bkRL8/4nr4PTeaFL1ewe6//Oios0fTyGg10AsYDv9cdzey6HMr/FfgLcDjwCnAm8A8zey9GMceE9/JyruxYsn4L93+yhKnL0mlauzK3DDiEE9oeiJTbLFIuOwVdD2VIdvvN7JVsyiYQDIDcBPQn6DI8wcwW72/Qhc0TinNlz5Rl6dz/8WKW/rSVzk0P4P+deAidm9aKd1glSoHXQwnn8mpiZkujKJvjGJXixBOKc2XT3n0ZvDd7DY+OX0b61l2c0PZAbh5wCM2Tq8Q7tBKhoCPlTwHmEsznhaSOksbkcsoESWfI65LOuWIoKTGBwV2bMPmmftxwzMFMWZbOsY9O4c7RC31G4wKK5pHXbOBoYHJmzSO3+boiRsrvJehi7CPlnXPF1oatOxn+xXLenrWaSuUSuarfQVzcqzmVyvuI++wUdKT8HjPbnGVfjt0kIkbKl/eR8s654q5utYrcN6gd467vQ/cWtXlo3FKOengy76WsZp+PuN8v0SSUVEnnAonh8r8jgBmFHJdzzhWplnWr8cKQLrxzeXfq1ajI39+fz0lPfMmUZek+VX6Uoh0pfxjBYllvApuB6wsxJueci5tuLWoz6uqejBjcie279zLkxW84738zWbAm64Mal1WObSiSKgJXAi2BBcD/zGxvEcZWqLwNxTmXl1179/HmzFWMmJjGpu27OaVDA/5+XGua1K4c79DiJl/dhiW9Q7CO/JfACcBKM7u+sIIsap5QnHPR2rpzD89NXcELX37P3owM/tqtKdce3ZLaVSvEO7Qil9+EssDM2oXbScA3ZnZ44YVZtDyhOOf214YtO3l8wnLeCXuEXX5kCy7t05zK5XOdFrFUyW8vrz2ZG6XpUZdzzuVX3eoVuX9QO8ZdfyS9Wybz6Phl9H1oMq9//QN79vkcYbkllA6StoSvrUD7zG1JWwojGEkDJC2VlCbp1myON5U0QdJ8SZMlNYo4NkTS8vCV7XQxzjkXCy3rVuXZ8zvzwVU9aVa7Mv8YtZDjH5vKJwt+LNM9wqKaeqUohDMULwOOBdYAs4DBZrYoosx7BNPpvyLpaOAiMztfUi0gBehCMMvxbKCzmf2S0/38kZdzLhbMjIlLNvDvz5aw7KdttG9Ug1sGHEKvlsnxDq1QFHRgY1HpCqSZ2Qoz202wmNdpWcq0ASaG25Mijh8PjDezTWESGQ8MKIKYnXNlnCT6H1qPT4cdycNndWDjtt389YWZnPfCTOav+TXe4RWp4pRQGgKrI96vCfdFmgecHm4PAqpJqh3luUi6XFKKpJT09PSYBe6cc4nh4l4Tb+rLHSe3YdGPWzj1yelc88a3fJe+Ld7hFYnilFCicRPQN1wNsi+wFtgX7clm9pyZdTGzLnXq1CmsGJ1zZViFpEQu7t2cKX/vx7D+rZi8dAPHPTaV20bOZ/3mnfEOr1AVp4SyFmgc8b5RuO93ZrbOzE4PJ6m8Pdz3azTnOudcUapWsRw3HHswU24+igt6NOX92Wvo+9AkHvh0Mb/u2B3v8ApFcWqUTyJolO9PkAxmAeeaWWpEmWRgk5llSLoP2Gdmd4SN8rMJVokE+JagUX5TTvfzRnnnXFFavWkHj32xjA/nrKVqhSSu7HsQF/ZsRpUKJWsMS4lolA/HugwFxgGLgXfNLFXS3ZJODYv1A5ZKWgbUA+4Lz90E3EOQhGYBd+eWTJxzrqg1rlWZR8/uyKfD/m9W474PTeKl6d+za2/UT+6LtWJTQylqXkNxzsXTt6t+4eFxS5nx3UYa1qzEsP6tOP3whiQlFpu/87NVImoozjlXlhze5ADevKw7b1zajeRqFbj5g/kc9/hUPp7/IxkldB0WTyjOORdHvVomM+rqnjx3fmeSEsQ1b37LKU9OY9LSDSVu1L0nFOecizNJHHfYgXw67Ege/0tHtu7cy0UvzeLs/37FzBUb4x1e1LwNxTnnipndezN4N2U1T0xYzoatu+jTKpkbj2tNx8Y14x1a/qavL+08oTjnirude/bx+tc/8PTk79i0fTfHHFqXvx3bmjYNqsctJk8o2fCE4pwrKbbt2ssrM1by3ynfsWXnXk5qX58bjmlFy7rVijwWTyjZ8ITinCtpNv+2hxe+XMGL077ntz37GNipIdf3P7hIlyT2hJINTyjOuZJq47Zd/HfqCl6ZsZJ9GcZZXRpz7dEtaVCzUqHf2xNKNjyhOOdKup+27OSpSWm89c0qhBjctTFXH9WSetUrFto9PaFkwxOKc660WPPLDp6cmMZ7s9eQlCDO696UK/seRJ1qFWJ+L08o2fCE4pwrbX7YuJ0RE9MY+e0ayiclMKRHMy4/sgW1q8YusXhCyYYnFOdcabUifRtPTFjO6HnrqFQukYt6NeOyPi2oWbl8ga/tCSUbnlCcc6Vd2oatPP7FcsbO/5GqFZK4uHdzLundnBqVyuX7mp5QsuEJxTlXVixZv4XHxy/ns9T1VKuYxL/PaM+J7ern61q5JZSStbKLc865/XbIgdV59vzOLFy7mce/WE6z2lUK5T6eUJxzroxo27AGLwzJtnIREz7bsHPOuZjwhOKccy4mPKE455yLCU8ozjnnYsITinPOuZjwhOKccy4mPKE455yLCU8ozjnnYqLMTr0iKR34IZ+nJwM/xzCcksA/c9ngn7lsKMhnbmpmdbI7UGYTSkFISslpLpvSyj9z2eCfuWworM/sj7ycc87FhCcU55xzMeEJJX+ei3cAceCfuWzwz1w2FMpn9jYU55xzMeE1FOecczHhCcU551xMeELZT5IGSFoqKU3SrfGOp7BJaixpkqRFklIlDYt3TEVBUqKkOZLGxjuWoiCppqT3JS2RtFhSj3jHVNgk3RD+n14o6S1JFeMdU6xJelHSBkkLI/bVkjRe0vLw5wGxup8nlP0gKRF4CjgBaAMMltQmvlEVur3AjWbWBugOXFMGPjPAMGBxvIMoQsOBz8zsEKADpfyzS2oIXAd0MbO2QCJwTnyjKhQvAwOy7LsVmGBmrYAJ4fuY8ISyf7oCaWa2wsx2A28Dp8U5pkJlZj+a2bfh9laCXzQN4xtV4ZLUCDgJeCHesRQFSTWAI4H/AZjZbjP7Na5BFY0koJKkJKAysC7O8cScmU0FNmXZfRrwSrj9CjAwVvfzhLJ/GgKrI96voZT/co0kqRnQCZgZ51AK2+PAzUBGnOMoKs2BdOCl8DHfC5KqxDuowmRma4GHgVXAj8BmM/s8vlEVmXpm9mO4vR6oF6sLe0JxUZFUFfgAuN7MtsQ7nsIi6WRgg5nNjncsRSgJOBx4xsw6AduJ4WOQ4ihsNziNIJk2AKpIOi++URU9C8aNxGzsiCeU/bMWaBzxvlG4r1STVI4gmbxhZiPjHU8h6wWcKmklwSPNoyW9Ht+QCt0aYI2ZZdY83ydIMKXZMcD3ZpZuZnuAkUDPOMdUVH6SVB8g/LkhVhf2hLJ/ZgGtJDWXVJ6gEW9MnGMqVJJE8Gx9sZk9Gu94CpuZ3WZmjcysGcG/70QzK9V/uZrZemC1pNbhrv7AojiGVBRWAd0lVQ7/j/enlHdEiDAGGBJuDwFGx+rCSbG6UFlgZnslDQXGEfQKedHMUuMcVmHrBZwPLJA0N9z3/8zsk/iF5ArBtcAb4R9KK4CL4hxPoTKzmZLeB74l6Mk4h1I4BYukt4B+QLKkNcCdwIPAu5IuIVjC4+yY3c+nXnHOORcL/sjLOedcTHhCcc45FxOeUJxzzsWEJxTnnHMx4QnFOedcTHhCcXEjaaAkk3RIHO69UlJyEd2rQdhFNa9yM8KfzSSdG0X5P5ST1EXSEwWLdv9IulBSg3ycd72kC/ajfHlJU8N5t1wx5QnFxdNgYFr4s8TY319qZrbOzM6MolzmSO1mQJ4JJWs5M0sxs+v2J7ZohLNs5+RCgqlL9ud6ScDFwJvRnhNOxjoB+Mv+3MsVLU8oLi7CucF6A5cQMW24pH6SJkeszfFGOJI5s1Zxl6RvJS3IrNlI+pekmyKusTCcyBJJoyTNDte9uDyKuLZJeiwsP0FSnXD/ZEmPS0oBhknqLGlKeO1xEVNZtJT0haR5YZwHhTWJheHxCyWNDq+3XNKdkfcONx8E+kiaG67Z0UzSl+H1vpXUM4dy/RSu36JgzYtRkuZL+lpS+4jv6sXw/iskZZuAwu/hEUnzgB6S7pA0K/xun1PgTKALwYDIuZIq5fS9ZHE08K2Z7Y34bh+TlKJgLZYjJI0Mv597I84bBfw1r39DF0dm5i9/FfmL4BfD/8LtGUDncLsfsJlgnrQE4Cugd3hsJXBtuH018EK4/S/gpohrLwSahdu1wp+Vwv21I66VnE1cBvw13L4DeDLcngw8HW6XC2OuE77/C8GsCRDMxDwo3K5IMC16M2BhuO9Cgtlta0fE1CU8ti3iOxgbEVNloGK43QpIyaHc7++BEcCd4fbRwNyI72oGUAFIBjYC5XL4Hs6OeF8rYvs14JSI76VLXt9LlmvflfnvGHGNf4fbwwimka8fxrgm4t8sEUiP9/9df+X88ueRLl4GEyzqBMEkjIOBzBl+vzGzNQAKpntpRvBoDIJJ/AjLnh7Ffa6TNCjcbkzwC3ljLuUzgHfC7dcj7kfE/tZAW2B8WHlKBH6UVA1oaGYfApjZzvAzZL3HeDPbGB4bSVBTS8klpnLAk5I6AvuAg3Mpm6k3cEYYx0RJtSVVD499bGa7gF2SNhBMX74my/n7CCYEzXSUpJsJklstIBX4KMs52X4v2cRWnz/Pm5U5J94CINXC6dUlrSD4d9toZvsk7ZZUzYK1eVwx4wnFFTlJtQj+am4nyQh+8Zikv4dFdkUU38cf/5/uymb/Xv74+LZieJ9+BLPK9jCzHZImZx7bD5FzE23P/AgEv/T+sExumFD295rZvc/qBuAngpUUE4CdUd4nJ7l9v5l2mtk+AAVL4z5NUBNZLelfZP89Zvu9ZOO3bM7PjCkjS3wZWeKrQME/vysk3obi4uFM4DUza2pmzcysMfA90Cef11tJON26pMMJ1rgAqAH8EiaTQwiWMM5LQhgfBA3e07IpsxSoo3DddUnlJB0W/tW8RtLAcH8FSZWzOf/YsI2jEsFqedOzHN8KRCanGsCPZpZBMFFnYg7lIn1J2N4QJtafLf/r2GT+8v9ZQdtXZAeDyBiy/V6yud5ioOX+BiGpNsHn2LO/57qi4QnFxcNg4MMs+z4g/729PgBqSUoFhgLLwv2fAUmSFhM0YH8dxbW2A13DRvSjgbuzFrCgx9GZwL/DRuu5/N9aGucTPGabT9CecGA29/gmjHk+8IGZZX3cNR/YFzbs30BQOxgS3usQ/q+mlLVcpH8BncM4HuT/pivfbxYsB/w8QXvPOIJlHDK9DDwbPppMJOfvJdKnBEsO76+jgI/zcZ4rIj7bsHMRJG0zs6qFeP0LCR4dDS2se5QEkj4Ebjaz5ftxzkjgVjNblmdhFxdeQ3HOxcOtBI3zUVGwTssoTybFm9dQnHPOxYTXUJxzzsWEJxTnnHMx4QnFOedcTHhCcc45FxOeUJxzzsXE/wepGPqqXTjsqgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "precip = np.linspace(0,10,30)\n",
    "frac_remain = []\n",
    "for i in precip:\n",
    "    Composting.InputData.Precip['ann_precip']['amount'] = i\n",
    "    Composting.calc(Inc_flow=DewateredWWTSolids)\n",
    "    frac_remain.append(sum(Composting.report()['Compost'])/sum(DewateredWWTSolids.PFAS))\n",
    "\n",
    "plt.plot(precip,frac_remain)\n",
    "plt.xlabel('Annual precipitation rate (m)')\n",
    "plt.ylabel('Percent of Incoming PFAS that \\n remains in the Compost (%)')"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
