{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Composting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import PFAS_SAT as ps\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from IPython.display import Image\n",
    "import pandas as pd\n",
    "pd.set_option('display.max_colwidth', 0)\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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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>NaN</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            NaN        \n",
       "14  10e-6g/kg            NaN        \n",
       "15  10e-6g/kg            NaN        \n",
       "16  10e-6g/kg            NaN        \n",
       "17  10e-6g/kg            NaN        \n",
       "18  10e-6g/kg            NaN        \n",
       "19  10e-6g/kg            NaN        \n",
       "20  10e-6g/kg            NaN        \n",
       "21  10e-6g/kg            NaN        \n",
       "22  10e-6g/kg            NaN        \n",
       "23  1:TRUE,0:FALSE       NaN        \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             NaN        \n",
       "30  kg TS/kg C           [64,24]    \n",
       "31  1:TRUE,0:FALSE       NaN        \n",
       "32  fraction             NaN        \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           NaN        \n",
       "40  days                 [48]       \n",
       "41  m/yr                 NaN        \n",
       "42  NaN                  NaN        \n",
       "43  fraction             NaN        "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Composting = ps.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</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>800</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>VS flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>120</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>76</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.93</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    \n",
       "2   Moisture flow  kg   800    \n",
       "3   VS flow        kg   120    \n",
       "4   Carbon flow    kg   76     \n",
       "5   PFOA           μg   934.565\n",
       "6   PFOS           μg   3678.93\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 = ps.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": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deViVdf7/8edhV1ASF0RxQxSVw2ERGFBAIZdyoZzc0DTX3MqydLLpW1q/ctRxcq0Ya9KxXFDT0TLNDddyABU3xAXBBdDU1FhkO3x+f/jtfCNAQdGb5f24Lq7Lc5/7/tzv+4QvT/fy/uiUUgghhHjyzLQuQAghaioJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjViUZ+UGDRqoli1bPqZShBCiejp8+PANpVTDPy4vVwC3bNmSuLi4iqtKCCFqAJ1Od7Gk5XIKQgghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAolLq2rUrP/zwQ5FlCxYsYOLEiSWun5KSgl6vv++YKSkprFq1yvQ6Li6OyZMnA7B8+XJeeeUVACIjI1mxYsWjlA/cu2voxo0bjzyOqL4kgEWlFBERwZo1a4osW7NmDREREQ895h8D2NfXl0WLFhVbb/z48QwfPvyh9yNEWUkAi0qpf//+fPfdd+Tm5gL3wjMtLY2goCCmTZuGXq/Hw8ODqKioYtumpKQQHByMj48PPj4+/PjjjwBMnz6d/fv34+Xlxfz589mzZw99+vQptv3MmTOZN28eaWlpeHl5mX7Mzc25ePEi169f54UXXsDPzw8/Pz8OHjwIwM2bN+nRowfe3t6MGzcO6bUtHkQCWFRK9evXx9/fn23btgH3vv0OGjSIDRs2EB8fz7Fjx9i5cyfTpk0jPT29yLaNGjVix44dHDlyhKioKNNphtmzZxMcHEx8fDxTpkx5YA1NmjQhPj6e+Ph4xo4dywsvvECLFi147bXXmDJlCrGxsXzzzTeMGTMGgPfff5+goCCOHj1KeHg4ly5dquBPRVQ35XoSTogn6bfTEM899xxr1qzhyy+/5KuvviIiIgJzc3McHR3p0qULsbGxGAwG03b5+fm88sorxMfHY25uztmzZx+pjoMHD/LFF1+wf/9+AHbu3ElCQoLp/V9//ZWMjAz27dvHhg0bAOjduzf16tV7pP2K6k8CWFRazz//PG+88QZHjhzh7t27+Pj4lOni2Pz583F0dOTYsWMUFhZiY2Pz0DWkp6czevRoNm/ejJ2dHQCFhYX89NNP1KpVq9j6Op3uofclah45BSEqLTs7O7p27cqoUaNMF99CQkKIiorCaDRy/fp19u3bh7+/f5Ht7ty5g5OTE2ZmZnz11VcYjUYA6tSpQ0ZGRpn3n5+fz8CBA5kzZw5t27Y1Le/RowdLliwxvY6PjzfVtnLlSgC2bt3KrVu3Hu7ARY0hASwqtYiICI4dO8bgwYMB6NevHwaDAU9PT8LCwpg7dy6NGzcuss3EiRP597//TUBAAGfPnsXW1hYAg8GAhYUFnp6ezJ8//4H7/vHHH4mNjWXGjBmmC3FpaWksWrSIuLg4DAYDHTp0IDIyEoAZM2awb98+fHx82L59O82bN6/gT0NUN7ryXKn19fVV0g1NCCHKR6fTHVZK+f5xuXwDFkIIjUgACyGERiSARY20bds23NzccHV1Zfbs2VqXI2ooCWBR4xiNRiZNmsTWrVtJSEhg9erVRe7rFeJJkQAWNU5MTAyurq64uLhgZWXF4MGD2bRpk9ZliRpIAljUOKmpqTRr1sz02tnZmdTUVA0rEjWVBLCocUq69VKeYBNakAAWNY6zszOXL182vb5y5QpNmjTRsCJRU0kAC5OrV68yePBgWrduTYcOHejVq9dDNbJZsGAB2dnZD1VDfHw833///UNtW1Z+fn6cO3eO5ORk8vLyWLNmDeHh4Y91n0KURAJYAPf+t7xfv3507dqVpKQkEhISmDVrFteuXSv3WJU9gC0sLFiyZAk9e/akffv2DBw4EHd398e6TyFKpJQq80/Hjh2VqJ527dqlgoODiy0vLCxUU6dOVe7u7kqv16s1a9YopZSKjo5WXbp0US+88IJyc3NTQ4YMUYWFhWrhwoXK0tJS6fV61bVrV6WUUuPHj1cdO3ZUHTp0UO+9955p7JiYGBUYGKgMBoPy8/NTt2/fVs2aNVMNGjRQnp6epn0JUdUBcaqETJUAFkoppRYuXKhef/31YsvXr1+vunXrpgoKCtTVq1dVs2bNVFpamoqOjlZ169ZVly9fVkajUQUEBKj9+/crpZRq0aKFun79ummMmzdvKqWUKigoUF26dFHHjh1Tubm5qlWrViomJkYppdSdO3dUfn6+WrZsmZo0adITOGIhnpzSAlhOQYj7OnDgQIkN0AH8/f1xdnbGzMwMLy8vUlJSShxj7dq1+Pj44O3tzalTp0hISODMmTM4OTnh5+cHQN26dbGwkPbUomaRABYAuLu7c/jw4WLL1X265VlbW5v+bG5uTkFBQbF1kpOTmTdvHrt27eL48eP07t2bnJwclFJy65eo8SSABQBhYWHk5uby+eefm5bFxsZSr169BzZA/6PfNz7/9ddfsbW1xd7enmvXrrF161YA2rVrR1pamunbdEZGBgUFBeVumi5EVSYBLIB7DyJs3LiRHTt20Lp1a9zd3Zk5cyZDhgx5YAP0P3r55Zd59tlnCQ0NxdPTE29vb9zd3Rk1ahSdO3cGwMrKiqioKF599VU8PT3p3r07OTk5hIaGkpCQgJeXV4kzHgtRnUhDdiGEeMykIbsQQlQyEsBCCKERCWAhhNCIBLCocI/yKLIQNYkEsKhwDxPARqPxMVVTuvT09Pve5yzE4yYBLEo1d+5cFi1aBMCUKVMICwsDYNeuXbz44otMmDABX19f3N3dmTFjBgCLFi0iLS2N0NBQQkNDAdi+fTuBgYH4+PgwYMAAMjMzAWjZsiUffPABQUFBrFu37okdl1KK92a+T1NnZwYPHUZubu4T27cQRZT0fHJpP9ILomb56aefVP/+/ZVSSgUFBSk/Pz+Vl5enZs6cqSIjI0vs8aBU0V4Q169fV8HBwSozM1MppdTs2bPV+++/b1pvzpw5T/qw1OHDhxVg+lm2fPkTr0HULJTSC0Ievhel6tixI4cPHyYjIwNra2t8fHyIi4tj//79LFq0iLVr17J06VIKCgpIT08nISEBg8FQZIxDhw6RkJBgegAjLy+PwMBA0/uDBg16oscE4O3tTXp6Ok5OTuzatYvg4OAnXoMQgASwKJ2lpSUtW7Zk2bJldOrUCYPBQHR0NElJSdSqVYt58+aZHlceMWIEOTk5xcZQStG9e3dWr15d4j5sbW0f92EUo9PpTE/zNWrUCEtLyydegxAg54DFA4SEhDBv3jxCQkIIDg4mMjISLy+vUns8QNFeEAEBARw8eJDz588DkJ2d/VCzbAhRHUkAi/sKDg4mPT2dwMBAHB0dsbGxITg4uNQeD1C0F0TDhg1Zvnw5ERERGAwGAgICSExM1PCIhKg8pBeEqLF0Oh0nTpxAr9drXYqo5qQXhBBCVDISwEIIoREJYCGE0IgEsChm27ZtuLm54erqyuzZs4u9r5Ri8uTJuLq6YjAYOHLkyAO3fffddzEYDHh5edGjRw/S0tIAyM/P56WXXsLDw4P27dvzt7/97fEfoBCVRUlPZ5T2I0/CVX8FBQXKxcVFJSUlqdzcXGUwGNSpU6eKrLNlyxb1zDPPqMLCQvXTTz8pf3//B257584d0/YLFy5U48aNU0optXLlSjVo0CCllFJZWVmqRYsWKjk5+QkcqVKAOnHixBPZl6jZkFmRRVnExMTg6uqKi4sLVlZWDB48mE2bNhVZZ9OmTQwfPhydTkdAQAC3b98mPT39vtvWrVvXtH1WVpZpQk6dTkdWVhYFBQXcvXsXKyurIusKUZ1JAIsiUlNTadasmem1s7MzqampZVrnQdu+8847NGvWjJUrV/LBBx8A0L9/f2xtbXFycqJ58+ZMnToVBweHx3V4QlQqEsCiCFXCfeF/nD6+tHUetO1HH33E5cuXGTp0KEuWLAHufeM2NzcnLS2N5ORk/vGPf3DhwoVHPQwhqgQJYFGEs7Mzly9fNr2+cuUKTZo0KdM6ZdkWYMiQIXzzzTcArFq1imeeeQZLS0saNWpE586dkYd9RE0hASyK8PPz49y5cyQnJ5OXl8eaNWsIDw8vsk54eDgrVqxAKcWhQ4ewt7fHycnpvtueO3fOtP3mzZtp164dAM2bN2f37t0opcjKyuLQoUOm94So7qQbmijCwsKCJUuW0LNnT4xGI6NGjcLd3Z3IyEgAxo8fT69evfj+++9xdXWldu3aLFu27L7bAkyfPp0zZ85gZmZGixYtTONNmjSJkSNHotfrUUoxcuTIYi0thaiupBeEqLGkF4R4UqQXhBBCVDISwEIIoREJYCGE0IgEsBBCaEQCWFSodevW8+2332ldhhBVggSwqBBKKV5/402GDHuJPw8YwIz3P9C6JCEqPbkPWFSY+OMnKcjNBiDu6DGNqxGi8pNvwKJC6HQ69uz8AQC7OnXY8p9vNK5IiMpPAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjUgAiwpnbm6udQlCVAkSwKLCtWrVSusSHig/Px8AR0dHjSsRNZkEsKgwFy9eBGDB/PkaV/Jgr776KgANGzbUuBJRk0kAiwqhlKJly5YAdOnSRdtiHiA3N5d//vOfhIaGal2KqOEkgMUjS01NxdLSEoAff/xR42ru7+rVq9jY2ACwZcuWCht3zPiJ/PDDDxU2nqgZdEqpMq/s6+ur4uLiHmM54vbt27z25rRiy1d8+YUG1ZTPyJEj6dSpk9ZllGjs2LFFXs+ZMwcHB4cKG//1N6eiChWHY/9Lu3btKmxcUT3odLrDSinfYsslgCuXpKQk9AZPzF0CsW7S1rT89t5/U5ibpWFlpbOs24Cn2vgCOq1LKdX1w1sBqOXYEjvn9hU+vrK2RdVtjHXCdxw/Ekf9+vUrfB+i6pIArkK2bNnCoGEjeWrQHCzsG2ldjiiHzH3LaaXSObh3N1ZWVlqXIyqJ0gJYzgFXQr1792bmO9PJ2Pwhhf87yaWoGmyDh3PhV8XIMS9Tni83omaSAK6k3nxjCs/36Ermtn+gCo1alyPKSKczw7bnFLbsPsjsuX/XuhxRyUkAV1I6nY4vlkbiVt+arP0rtC5HlIOZVS3swv/KR3PmsXnzZq3LEZWYBHAlZmlpyZbNG6mVfpSs49u1LkeUg0XdRtTpM52hw0dw/PhxrcsRlZQEcCXn4ODAru1byTu0kpxL8he5KrFu4oZNyBi6P9uba9euaV2OqIQkgKuAtm3b8s3aNWRu/Qf5t9K0LkeUg22HLhhbh9CjV19ycnK0LkdUMhLAVUS3bt2YO+v/kbn5Q4w5mVqXI8rBNnAwV/JrM2TYS3JnhChCArgKmThhAhH9wsn6/u8oY4HW5Ygy0unMsO0+md0xJ5jx/gdalyMqEQngKmbJovkYmjmQtfdfWpciysHM0hq7vm8zf8lnREWt1bocUUlIAFcxFhYWbNqwjjq3zpJ15DutyxHlYGHnQJ2+f2X0uPHIE6UCJICrJHt7e3Zv30rBkW+4e+Gw1uWIcrBybE3tsAk807svV65c0bocoTEJ4CqqVatWfPufDWRtX0jejUtalyPKoXbbTuD+DN2e6U1WVuVssCSeDAngKiwoKIjF8/9B5uaPMGbf0bocUQ62fi9w3bIRAwYPpbCwUOtyhEYkgKu4kSNHMGZ4BJlb5qAK8rUuR5SRTqfD7ulJ/HTqAn+Z/letyxEasdC6gMrCzs6OzMwnf3/te++9R0hICN26dXvoMebNnUPC6TPsWByBeV1H0Okwt61Hg95vYG5XjyufjcLMqhaY3fv31qH7RGz+tyfur7H/4dbef9Ps1a8xs7YFoDA/h5vbFpP/cwqgMLO2o9HA9++NISqMzsISuz7T+efyaeg7tGfEiJfuu/7Vq1d5/fXXiY2NxdrampYtW7JgwQLatm173+0et9u3b7Nq1SomTpyoaR1VkfQD/l9aBXBFyczMpJ6DA3X/NIA6nYdwa++/Ufk5OHQbx5XPRuH00nzMa9sX2y59xRR05pbYGXpg53HvH4E7P63FmH0Hh6fvzSKRf/MKFvaO6Cwsn+gx1RR51y9yZ/3/sP37b+ncuXOJ6yil6NSpEy+99BLjx48HID4+noyMDIKDg59kucWkpKTQp08fTp48qWkdlZn0Ay6jPXv20LVrV/r370+7du0YOnSo6eml2NhYOnXqhKenJ/7+/mRkZJCTk8PIkSPx8PDA29ub6OhoAJYvX87zzz9P3759adWqFUuWLOHjjz/G29ubgIAAfvnlFwBGjBjB+vXrAWjZsiUzZszAx8cHDw8PEhMTAbh+/Trdu3fHx8eHcePG0aJFC27cuFGkbjs7Oxo7OlKYuJPsc//Fppn+gY8t599KR+Xl8FTwMLIS9pqWG7NuYVHn/2Z0sKzvLOH7GFk1bIFdj9fo81w/UlJSSlwnOjoaS0tLU/gCeHl5ERQUxLRp09Dr9Xh4eBAVFQXc+z3u0qULAwcOpG3btkyfPp2VK1fi7++Ph4cHSUlJwL3fv/HjxxMcHEzbtm357rt7tzaW9nt96tQp/P398fLywmAwcO7cOaZPn05SUhJeXl5Mm1Z8Oi1ROjkFUYKjR49y6tQpmjRpQufOnTl48CD+/v4MGjSIqKgo/Pz8+PXXX6lVqxYLFy4E4MSJEyQmJtKjRw/Onj0LwMmTJzl69Cg5OTm4uroyZ84cjh49ypQpU1ixYgWvv/56sX03aNCAI0eO8OmnnzJv3jy++OIL3n//fcLCwnj77bfZtm0bS5cuLbFuc3NzolZ+xZ8HDkY56bFq2NL03rXVfwUzM3TmljgN/xiArNN7qd0+BOtm7uRvScWYdRtz26ew8+jOtbXvknXmIDYtvLDTh2Hp0LSCP2Xxe7Va+2K83Y+wHs8SH/df6tatW+T9kydP0rFjx2Lbbdiwgfj4eI4dO8aNGzfw8/MjJCQEgGPHjnH69GkcHBxwcXFhzJgxxMTEsHDhQhYvXsyCBQuAe99g9+7dS1JSEqGhoZw/f55PPvkEKP57HRkZyWuvvcbQoUPJy8vDaDQye/ZsTp48SXx8/GP+lKof+QZcAn9/f5ydnTEzM8PLy4uUlBTOnDmDk5MTfn5+ANStWxcLCwsOHDjAsGHDAGjXrh0tWrQwBXBoaCh16tShYcOG2Nvb07dvXwA8PDxK/abz5z//GYCOHTua1jlw4ACDBw8G4JlnnqFevXql1v7mm2/S4Km65F2Ixc6ju2m5Y8QsmoxcbApfgOzT+7BtH4JOZ0bttoFknzkAgJWjC03HfYG9/wsU5mSQvuIN8m9cLu/HKMrJ1qcvv9q78PwLAzEay9aE/8CBA0RERGBubo6joyNdunQhNjYWAD8/P5ycnLC2tqZ169b06NEDKP77N3DgQMzMzGjTpg0uLi4kJiaW+nsdGBjIrFmzmDNnDhcvXqRWLbku8CgkgEtgbW1t+rO5uTkFBQUopdDpik86eb9z6L8fx8zMzPTazMyMgoKSezn8ts5v+33QPv4oOjqalJQU3nrrL2TvWEhhfm6J6+X9nEz+rTSuRb3Llc9GkZ24n6yEff9Xr1Utart1on6Pidi5d+Xuhdgy1yAejk6nw7bry8Sn/Myrr79R5D13d3cOHy7+0E1F/P798fdap9OVOu6QIUPYvHkztWrVomfPnuzevfvBByZKJQFcRu3atSMtLc307SIjI4OCggJCQkJYuXIlAGfPnuXSpUu4ublV6L6DgoJYu/Ze/4Dt27dz69atB27z/96fSVdfPVk7Fpf4ftbpvTzVeQjOE7689zNpBQWZNym48zM5VxJMHdeUMZ+8G5cxryuTgz4JOnMLbHv9ha/XbeSzyEjT8rCwMHJzc/n8889Ny2JjY6lXrx5RUVEYjUauX7/Ovn378Pf3L9c+161bR2FhIUlJSVy4cAE3N7dSf68vXLiAi4sLkydPJjw8nOPHj1OnTh0yMjIq5gOoYSSAy8jKyoqoqCheffVVPD096d69Ozk5OUycOBGj0YiHhweDBg1i+fLlRb55VIQZM2awfft2fHx82Lp1K05OTtSpU+e+2+h0OlZ99W+aWmSiSpjOPuv0fmq3DSyyrHabALJO76PgdjrXVk0n7V+TSF/2GtaNXantVvLVeVHxzGvVoU74/zBt+jumb5g6nY6NGzeyY8cOWrdujbu7OzNnzmTIkCEYDAY8PT0JCwtj7ty5NG7cuFz7c3Nzo0uXLjz77LNERkZiY2NT6u91VFQUer0eLy8vEhMTGT58OPXr16dz587o9Xq5CFdOchtaFZCbm4u5uTkWFhb89NNPTJgwocwXPK5evYrBxxfl9yK27bW9XUmUz92Lx8j54WNiD/342O71HTFiBH369KF///6PZXxxT2m3ocldEFXApUuXGDhwIIWFhVhZWRX539AHady4Mdu//46grmFY2DfCuknFnh4Rj0+tFp4Y/QfTrWcvjh2Jve/FV1E1yTfgGuKt6dP5YvNe6oT/j9aliHIwZt3m+rKJ7Nn5AwEBAVqXIx6SPIhRg8XExPBJ5FKsA4ZoXYooB1WQT+aWObw6aWKJ4btt2zbc3NxwdXVl9uzZxd6/desW/fr1w2Aw4O/vX+RJtfnz5+Pu7o5eryciIqLYfHXz5s1Dp9OZHviJiYnBy8sLLy8vPD092bhxYwUfbc0kAVzJjBo1ikaNGqHX6ytkvMuXL/Ns3+eo/fQrWDVyqZAxxeOnlCJz5yd00rdmzt8+Kva+0Whk0qRJbN26lYSEBFavXk1CQkKRdWbNmoWXlxfHjx9nxYoVvPbaawCkpqayaNEi4uLiOHnyJEajkTVr1pi2u3z5Mjt27KB58+amZXq9nri4OOLj49m2bRvjxo0r9VZKUXYSwJXMiBEj2LZtW4WMlZmZSbdnemOm70XtNn+qkDHFk5EVs55GxuusXf01ZmbF/5rGxMTg6uqKi4sLVlZWDB48mE2bNhVZJyEhgaeffhq4dxtlSkoK165dA6CgoIC7d+9SUFBAdnY2TZo0MW03ZcoU5s6dW+T+4Nq1a2Nhce+SUU5OTon3xIvykwCuZEJCQnBwcHjkcQoLC+k/aAg3rZtQ27dfBVQmnpTsMz/C6e3s3LYFW1vbEtdJTU2lWbNmptfOzs6kpqYWWcfT05MNGzYA9wL74sWLXLlyhaZNmzJ16lSaN2+Ok5MT9vb2pqfkNm/eTNOmTfH09Cy2z//+97+4u7vj4eFBZGSkKZDFw5MArqamvTWdmMRL2D49Xr6tVCG5V8+THf0ZP2z5lqZNS++/UdLF8z/+d54+fTq3bt3Cy8uLxYsX4+3tjYWFBbdu3WLTpk0kJyeTlpZGVlYWX3/9NdnZ2Xz00Ud88EHJMzf/6U9/4tSpU8TGxvK3v/2t2HljUX7yT1g1tGzZcj7/9yrsB89FZy5dzKqKgoybZH73N75c+s8SG+/8nrOzM5cv/19/jitXrhQ5jQD3+pUsW7YMuBfYrVq1olWrVvzwww+0atWKhg0bAvf6j/z44494enqSnJxs+vZ75coVfHx8iImJKfJwR/v27bG1teXkyZP4+ha7sC/KQb4BVzMHDhzg1SlvYhf+Ton9f0XlVJifQ+Z3f+PNyZMYOHDAA9f38/Pj3LlzJCcnk5eXx5o1awgPDy+yzu3bt8nLywPgiy++ICQkhLp169K8eXMOHTpEdnY2Sil27dpF+/bt8fDw4OeffyYlJYWUlBScnZ05cuQIjRs3Jjk52XTR7eLFi5w5c4aWLVtW+OdQ08g34GrkwoUL9H3+z9j2eA2rBs0fvIGoFJQqJGv7IroFeDHj3bLdp21hYcGSJUvo2bMnRqORUaNG4e7uTuT/9o8YP348p0+fZvjw4Zibm9OhQwf+9a9/AfdOJfTv3x8fHx8sLCzw9vbm5Zdfvu/+Dhw4wOzZs7G0tMTMzIxPP/2UBg0aPNqBC3kQo7KJiIhgz5493LhxA0dHR95//31Gjx79wO1+/fVXPDv6k9kqFFufPk+gUlFRMg+uxPluEocO7MXGxkbrcsRjII8iVxGrV68u9zYFBQWE/3kAGfXaYifhW6VkJezB/MJ+fjgSJ+FbA8k54Grg1dfe4Pilm9h2HaN1KaIcclMTydn/JTu3fY+jo6PW5QgNSABXcZ9+9hmrNmzCrtc0dGbmWpcjyqjg15/J2DKbVSuW4+HhoXU5QiMSwFXYzp07+ctf38Uu/H8ws7HTuhxRRoW52WRunsV7b79lmqZK1EwSwFXUmTNneGHgYOyefRPLek0evIGoFFShkcwf5tPn6SCmTX1T63KExiSAq6BffvmFbj17YRUwFJvmBq3LEeWQdWAFrk+Z8+Xn/5QnFIXcBVHV5Ofn0/u5P3PXyRs7Qw+tyxHlkHV8OzapR9h6JBYrKyutyxGVgARwFaKUYsy4CZy5mUudPsO1LkeUQ86lE+QeWsn+Hw9Qv359rcsRlYScgqhC/vHxfDZt34NdzzfkjocqJP9WOplb5/FN1GratWundTmiEpEAriK2bNnCzI9mY9f3Hcysa2tdjiijwpxMMjd/yOwP36d79+5alyMqGQngKuDkyZMMfnE4dXr/BQv7RlqXI8pIFRrJ3DqPQc/34pVJk7QuR1RCEsCV3M8//0z3Z3tjEzwS66bttS5HlEPWnn+hb/oUny5epHUpopKSi3CVWG5uLs/0eY78lp2x6xCqdTmiHDKPbqHOL4l8uyNGZo4QpZLfjEpKKcWwEaO4eNcKu9AIrcsR5XA3+SgFcevYGXMIe3vpySxKJ6cgKqkPP5rF9oOHse0xGZ1O/jNVFfk3L5P1w3y+/c8GWrdurXU5opKTb8CV0DfffMPc+YuxHzwHM8ua3aIw/5dUVKHxkccx3rmGztIaKPr0mWX9ZlABT6TpzC1RhQVkbP6IhR//neDg4EceU1R/EsCVzLlz5xgydBiWrf9E9ul9puW3or/UsA6d+6EAABIeSURBVCqo7dYZ6yZuRZZlHt9J/s1L2Bl6YFnf+bHs93Eft01tWyytrB95nJzcXGo1aMbIoQMZPWpUBVQmagIJ4EqmcePGvDx+XLHlS6I1KOZ3WlneIdRQt8iyJdGXAMg8vh2Avs89T4uKnifM8HqFDJOdfZevV/ybvNx7M/lu2LCBfv36VcjYAI2cnPFyb83Hf59bYWOK6k+mJBKP7Pz587Rp0waAlJQUWrRooXFFJVNKsW/fPrp27QrA7t27CQ2tmLtLEhMTcXZ2xs5O2oKK4kqbkkgCWFQIpRRmZmamP1dmv/2D0dCxMT9fTde6HFEDlBbAcnldVAidTsdv/zinpqZqXM39ubq6MnbcBK5fu8qRI0e0LkfUYBLAosJ07NgRgPfee0/jSh5szKgRwL1ZRYTQigSwqHApKSlal/BAjRrd66lx+84djSsRNZkEsKhwlf0cMGA6X03lL1VUYxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmhEAljUOOs3bMAvMAiAxZ98Qt8/D9C4IlFTSQCLKuHq1asMHjyY1q1b06FDB3r16sXZs2dLXDclJQW9Xg/Anj176NOnT5H3vQwGatnYYOfqR25uHr16PF3iOAsWLCA7O7tcdZa0PyFKIwEsKj2lFP369aNr164kJSWRkJDArFmzuHbt2kON5+rqSvzhGAaFdWTb998xYfz4Etd7mAAWojwkgEWlFx0djaWlJeN/F5ReXl4EBQUxbdo09Ho9Hh4eREVF3XecrKwsRo0ahZ+fH6GhofTt9QxhYWEYjUamTp2Kh4cHBoOBxYsXs2jRItLS0ggNDTXNG7d9+3YCAwPx8fFhwIABZGZmArBt2zbatWtHUFAQGzZseHwfhKh2ZFZkUemdPHnSNNvG723YsIH4+HiOHTvGjRs38PPzIyQkpNRxPvroI8LCwvjyyy+5ffs2/v7+dOvWjRUrVpCcnMzRo0exsLDgl19+wcHBgY8//pjo6GgaNGjAjRs3+PDDD9m5cye2trbMmTOHjz/+mL/85S+MHTuW3bt34+rqyqBBgx7nRyGqGQlgUWUdOHCAiIgIzM3NcXR0pEuXLsTGxmIwGEpcf/v27WzevJl58+YBkJOTw6VLl9i5cyfjx4/HwuLeXwcHB4di2x46dIiEhAQ6d+4MQF5eHoGBgSQmJtKqVSvTrNAvvvgiS5cufRyHK6ohCWBR6bm7u7N+/fpiy8s784ZSim+++QY3N7diy3U63QO37d69O6tXry6yPD4+/oHbClEaOQcsKr2wsDByc3P5/PPPTctiY2OpV68eUVFRGI1Grl+/zr59+/D39y91nJ49e7J48WJTcB89ehSAHj16EBkZSUFBAQC//PILAHXq1CEjIwOAgIAADh48yPnz5wHIzs7m7NmztGvXjuTkZJKSkgCKBbQQ9yMBLCo9nU7Hxo0b2bFjB61bt8bd3Z2ZM2cyZMgQDAYDnp6ehIWFMXfuXBo3blzqOO+++y75+fkYDAb0ej3vvvsuAGPGjKF58+amsVatWgXAyy+/zLPPPktoaCgNGzZk+fLlREREYDAYCAgIIDExERsbG5YuXUrv3r0JCgqiRYsWT+QzEdWDrjz/G+fr66vi4uIeYzmiqtPpdISGhrJ7926tSxGi0tDpdIeVUr5/XC7fgIUQQiMSwEIIoREJYCGE0IgEsKiSytMb4n4e5XHj+Ph4vv/++4faVgiQABZVUEX2hpAAFlqSABZVTnl7Q+zZs4euXbvSv39/2rVrx9ChQ1FKldjvYcKECfj6+uLu7s6MGTNM48fGxtKpUyc8PT3x9/fnzp07vPfee0RFReHl5fXAPhRClEgpVeafjh07KiHuB1ChoaGPdR8LFy5Ur7/+erHl69evV926dVMFBQXq6tWrqlmzZiotLU1FR0erunXrqsuXLyuj0agCAgLU/v37lVJKtWjRQl2/ft00xs2bN5VSShUUFKguXbqoY8eOqdzcXNWqVSsVExOjlFLqzp07Kj8/Xy1btkxNmjTpsR6rqB6AOFVCpso3YFFtlNYbAsDf3x9nZ2fMzMzw8vIiJSWlxDHWrl2Lj48P3t7enDp1ioSEBM6cOYOTkxN+fn4A1K1b19Q3QohHIQEsqhx3d3cOHz5cbLm6z0NF1tbWpj+bm5ubHjv+veTkZObNm8euXbs4fvw4vXv3Jicnp0y9IoR4GBLAosqpqN4QULTfw6+//oqtrS329vZcu3aNrVu3AtCuXTvS0tJM36YzMjIoKCgosq0QD0MCWFQ5FdUbAor2e/D09MTb2xt3d3dGjRplaj1pZWVFVFQUr776Kp6ennTv3p2cnBxCQ0NJSEiQi3DioUkvCFGhpBeEEMVJLwghhKhkJICFEEIjEsCiQhQWFtK6bTvg3oMSL40eq3FFQlR+EsCiQpiZmTHltVexsqmNvUNDxo8ZpXVJZSJTzwstyUU4UaGOHz9O/fr1adq0qdallEnLli2Ji4ujQYMGZd7GaDRibm7+GKsS1Y1chBNPhMFg0CR8586dy6JFiwCYMmUKYWFhAOzatYsXX3yxxB4PJfWC2L59O4GBgfj4+DBgwAAyMzOBe0H9wQcfEBQUxLp165748YnqSQJYVAshISHs378fgLi4ODIzM8nPz+fAgQMEBwfz0UcfERcXx/Hjx9m7dy/Hjx9n8uTJNGnShOjoaKKjo7lx4wYffvghO3fu5MiRI/j6+vLxxx+b9mFjY8OBAwcYPHiwVocpqhl5oF1UCx07duTw4cNkZGRgbW2Nj48PcXFx7N+/n0WLFrF27VqWLl1KQUEB6enpJCQkYDAYioxx6NAhEhISTA9g5OXlERgYaHp/0KBBT/SYRPUnASyqBUtLS1q2bMmyZcvo1KkTBoOB6OhokpKSqFWrFvPmzTM9rjxixAhycnKKjaGUonv37qVOLW9ra/u4D0PUMHIKooarqJklKtrt27f59NNPy7VNSEgI8+bNIyQkhODgYCIjI/Hy8iq1xwMU7QUREBDAwYMHOX/+PADZ2dmV4rMQ1ZcEcA2mKnBmiYr2MAEcHBxMeno6gYGBODo6YmNjQ3BwcKk9HqBoL4iGDRuyfPlyIiIiMBgMBAQEkJiYWNGHJsT/KalJcGk/0pC9etm1a5cKDg4utrywsFBNnTpVubu7K71er9asWaOUUio6OlqFhISoAQMGqDZt2qi33npLff3118rPz0/p9Xp1/vx5pZRSL730kho3bpwKCgpSbdq0Ud9++61SSqm7d++qESNGKL1er7y8vNTu3buVUkqdPHlS+fn5KU9PT+Xh4aHOnj2rBg0apGxsbJSnp6eaOnXqE/pEhHg8KKUhu5wDrsFOnjxJx44diy3fsGED8fHxHDt2jBs3buDn50dISAgAx44d4/Tp0zg4OODi4sKYMWOIiYlh4cKFLF68mAULFgCQkpLC3r17SUpKIjQ0lPPnz/PJJ58AcOLECRITE+nRowdnz54lMjKS1157jaFDh5KXl4fRaGT27NmcPHmS+Pj4J/eBCPGEySkIUcz9Zpbw8/PDyckJa2trWrduTY8ePQDw8PAoMsvEwIEDMTMzo02bNri4uJCYmMiBAwcYNmwYcK/HbosWLTh79iyBgYHMmjWLOXPmcPHiRWrVqvXEj1kILUgA12CPOrOEmZmZ6bWZmVmRWSb+OIOETqcrddwhQ4awefNmatWqRc+ePcvVynLbtm24ubnh6urK7NmzSzyWyZMn4+rqisFg4MiRI6b3Ro0aRaNGjdDr9UW2WbduHe7u7piZmfH7Jz/z8vIYOXIkHh4eeHp6smfPnjLXKURJJIBrsIqcWeKP1q1bR2FhIUlJSVy4cAE3NzdCQkJYuXIlAGfPnuXSpUu4ublx4cIFXFxcmDx5MuHh4Rw/frxMs00YjUYmTZrE1q1bSUhIYPXq1SQkJBRZZ+vWrZw7d45z586xdOlSJkyYYHpvxIgRbNu2rdi4er2eDRs2mE67/Oa3z+nEiRPs2LGDN998k8LCwnJ9LkL8ngRwDVaRM0v8kZubG126dOHZZ58lMjISGxsbJk6ciNFoxMPDg0GDBrF8+XKsra2JiopCr9fj5eVFYmIiw4cPp379+nTu3Bm9Xs+0adNK3EdMTAyurq64uLhgZWXF4MGD2bRpU5F1Nm3axPDhw9HpdAQEBHD79m3S09OBe7etOTg4FBu3ffv2uLm5FVuekJDA008/DUCjRo146qmnkN4o4pGUdGWutB+5C0KUxUsvvaTWrVv32Pezbt06NXr0aNPrFStWFJsmvnfv3qYp6JVSKiwsTMXGxppeJycnK3d39xLH79KlS5F1//nPf6r+/fur/Px8deHCBWVvb6/Wr19fUYcjqjHkLghR3agSzin/8dxzWdYpq1GjRnH69Gl8fX1p0aIFnTp1kunpxSOR3x5R4ZYvX/5E9uPs7Mzly5dNr69cuUKTJk3KvU5ZWVhYMH/+fNPrTp060aZNm4caSwiQc8A1lrm5OV5eXuj1evr27cvt27cfapwxY8YUu/AF90L4lVdeeej67OzsHriOn58f586dIzk5mby8PNasWUN4eHiRdcLDw1mxYgVKKQ4dOoS9vT1OTk4PVVN2djZZWVkA7NixAwsLCzp06PBQYwkBEsA1Vq1atYiPj+fkyZM4ODiYHpIory+++EKzELKwsGDJkiX07NmT9u3bM3DgQNzd3YmMjCQyMhKAXr164eLigqurK2PHji3yeHNERASBgYGcOXMGZ2dn/vWvfwGwceNGnJ2d+emnn+jduzc9e/YE4Oeff8bHx4f27dszZ84cvvrqqyd/0KJ6KenEcGk/chGu+rC1tTX9+bPPPlMTJkwwvZ47d67y9fVVHh4e6r333lNKKZWZmal69eqlDAaDcnd3Nz2e/PsLVV9++aVq06aNCgkJUWPGjDFdEPvjRbnf9p2RkaHCwsKUt7e30uv16j//+U+xddLS0lRwcLDy9PRU7u7uat++fY/j4xDisUIuwomSGI1Gdu3axejRo4F7M0KcO3eOmJgYlFKEh4ezb98+rl+/TpMmTdiyZQsAd+7cKTJOeno6M2bM4PDhw9jb2xMaGoq3t/d9921jY8PGjRupW7cuN27cICAggPDw8CIXyVatWkXPnj155513MBqNMn+bqFbkFEQNdffuXby8vKhfvz6//PIL3bt3B+4F8Pbt2/H29sbHx4fExETOnTuHh4cHO3fu5K233mL//v3Y29sXGe+///0vXbt2pWHDhlhZWZWpeblSir/+9a8YDAa6detGampqsU5sfn5+LFu2jJkzZ3LixAnq1KlTcR+CEBqTAK6hfjsHfPHiRfLy8kzngJVSvP3228THxxMfH8/58+cZPXo0bdu25fDhw3h4ePD222/zwQcfFBuztNu7LCwsTE+MKaXIy8sDYOXKlVy/fp3Dhw8THx+Po6NjsUbpISEh7Nu3j6ZNmzJs2DBWrFhRkR+DEJqSAK7h7O3tWbRoEfPmzSM/P5+ePXvy5ZdfmiajTE1N5eeffyYtLY3atWvz4osvMnXq1CI9FQD+9Kc/sWfPHm7evEl+fn6RiStbtmxp6jmxadMm8vPzgXunMRo1aoSlpSXR0dFcvHixWH0XL16kUaNGjB07ltGjRxfbrxBVmZwDFnh7e+Pp6cmaNWsYNmwYp0+fNs2FZmdnx9dff8358+eZNm0aZmZmWFpa8tlnnxUZw8nJiZkzZxIYGIiTkxM+Pj4YjUYAxo4dy3PPPYe/vz9PP/20aWqfoUOH0rdvX3x9ffHy8qJdu3bFatuzZw9///vfsbS0xM7OTr4Bi2pFp+7T+eqPfH19lTz7LoQQ5aPT6Q4rpXz/uFxOQQghhEYkgIUQQiMSwEIIoREJYFFhHjQ7hRCiKAlgUSHKMjuFEKIoCWBRIcoyO4UQoigJYFEhUlNTadasmem1s7MzqampGlYkROUnASwqREn3kz/szBNC1BQSwKJCVOTME0LUFBLAokKUZXYKIURR0gtCVIjfz05hNBoZNWoU7u7uWpclRKUmvSCEEOIxk14QQghRyZTrG7BOp7sOFG/aKoQQ4n5aKKUa/nFhuQJYCCFExZFTEEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihkf8P/wcXSnuNiysAAAAASUVORK5CYII=\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": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd3gU1dvG8e9JQihJQBHB0KSDlBBaXlC6FKmKIgIqKjZAxYo/rFgBERuiYkVRgVhApEoRBQSkBoSIKEWBBKSTQvp5/9glBkiFLJNN7s915YKdnZnzzJZnnzlzZsZYaxERERER8QY+TgcgIiIiIpJbKl5FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl5FRERExGuoeM1nxpgSxhhrjKnsdCy5YYypY4w55nQcIiJOUd4W8S5Fong1xsRm+EszxpzM8PjmHJa9xhjzVz7GstoYc0t+re98WWu3W2sv8sS6jTHTjTGJ7tf5iDFmgTGmtvu5scaY5DPem+EZlu1jjFlnjIkzxhwyxkwxxgRn0sY17h+d4Wc+l4v4hhljdhpjYowxc4wxFbKZ9ytjzH5jzAljzDZjzKAMz9Vzx5BxWx7P4nU49bfmHOJtYYyJMMbEG2PWGGMaZnjubmPMRnd8e4wxLxtjisT3Wwon5e2seWvedu8kTDDG7HMvu9MY80ou4ypmjPnJGHPQGHPcne+657BMmDHmF3db0caYoRmea26MWele1x5jzMgMzw0xxqScsZ2xxphL8vhaXmqMme1+PXYZY/pmeC7UGLPBGHPU/Tr/YIypk5f1F2VF4sfNWht46g/4B+iVYdqXTsdXyL3oft2rAjHAhxme+yzje2OtnQDg/mGaDIwDLgFCAF9gmTGm9Bnrvw044v4314wxXYGngW5AOeAAMCWbRZ4HqlprSwM3AK9lLB6B1DO2ZdwZy794xvNhWcS12hjTMpPpJYFZwAfAxcDXwExjjJ97luLAfbheryuBXkCeC3qRgkJ521GeytujgCuApkAQ0AnYnMuYUoH7gWBrbRngAeCrrApKd9E8F3gLKAvUBZZmmOUrYL77uU7Ao8aYLhme/+mM7Qy01h7OpJ2xGQvfM3wAHAXKA3cCn5zaEcD1me7jbr88sBj4IqcXQVyKRPGaE2NMSWPMO+49s73GmFfde3mXADOBGhn3vIwxVxljfnXvsUUZY97IUETkpd1rjDF/GWOedO9N7svYo2CMCXDvpe5xt/XzqXaMMTcYYyKNMceMMYszfCEwrh7CR4wxW90xv2eMCTbGLDKunrkFp5KJcfUapmRYdrUxZpT73xPGmHnGmIszPH+XMeYfd7yPu9tqndO2WmtjgelAw+zmc2/feOBZa+1X1toEa20UcKqn8/4M85YBrgWGAI3PKCZz0guYZq39w1qbCLwMdDHGVMoi/q3W2qRTDwED1MhDe+erM5BgrX3XHe9ruJJ/a3d8E621K621ydbaPbhe66syW9GZ77l7WnrPkjHGz/25O2yM2WGMGX7m/CJOU972yrzdAvjWWnvAuuzM7Y6ItTbNWrvFWptijDFAGq6d9kxzNjAC+M4dT5K19oS1dps7Xl+gCvCltTbVWvsHsBpokJtYcsP9+vcCnrHWxllrfwR+AG52b88Ra+3f9r/bnKYCtbJZ32nvmXEVzR9leHxO77G3UvHq8jyuvcRGQDOgPfC4ey+rD7DzjD2vZFxfxrJAG1wf0LvOse3LcRVCFd3rnGSMCXQ/NwGoh+sLXxZXT6F1F2mfAsNw7bH9DHx/RiLuA7QD6gP9cfXaPQJUAAKBoWRtIK4vWDBwEfAguA5zAK8D/YDK7r9yudlId9IdAGzMYdaGwGW4ehbTWWtTgRm4irhT+gEHgW9w7VEPyriMMeYPY8z1WYXk/sv4+FT7WW3Dx8aYk8BWYAewKMPTvu4fsT3GmA+NMWWzWs85agBsOvXAWpsGbCHrZNvWHee5uB/XZ6chEAb0zX52EUcob5/OG/L2auB/xnVY/qzc5c6xr+cQ0yIgAfgFmAf8lsWsLYETxjXE6l9jzHenOifccb0N3Obe4WmA6zO0JIftzIt6QIy19u8M0zaRIWcbY4ob19jlROBVYMy5NHQ+77G3UvHqcjMwylp7yFp7AHgJuDWrma21a6y1a917bDuAj3AlnHMRD4xx95jNxNWrV8sYUwxXMfaAtXa/u63l7i/dAGCmtfYnd2/gaFwf1OYZ1vume3v+AVYCv1hrf7PWnsSVEJtkE9OH1tod1to4XIVhqHt6P1x7zavdvX9Pk/Nn6Cn3l/MPXIeQMv5Y3OrugTj1V9a9HRbXYfwzRXP6F/I2XL2nFpgK3OLeowbAWlvXWjsji7jmAwONMQ2MMaWAZ9ztlspqQ6y1d+L6AWmP6zU81RMbjeswWFXg/3D90EzO7HXI8Pd+Vu1kIRA4fsa047h6X09jXOO66gFv5rGNU/oBr1tro90/+mcOgRApCJS3T+cNeft5XHnpNmCDu8d8wKkZrbV3WmsfyS4wa21nXHmvN7AkQ8/lmSq727kXV24+c2jYd7jeq5O4OgImWmszDmFod8Z25rUzIMecba1NtK6xyxcBD5OhgyKPzuU99mqFeuNyw3344TIg497R32R9KAJjTH1jzHxjzAFjzAngWc59L+eguxftlHhcH/pgwA/YmckyFTPG606M+86IOWMSOZnJ40Cytj+TeE61uydDuyc4+8t5ppettRdZa4OttX3O2Av93P3cqb8jwCFcPRqZnTwV7H4eY0xNXIfFTx1y+hbXWNDOmSx3FmvtHOAV4Htcr/FvuIrRvTksl2qt/Rmog2sME9ba49baje7nonCNNe1hjCmRyetw6u9e93YUz5ggcf2QLcow7SH38rHAmeN9S+Maj5bOGNMP17iya6y153o28mnv8xn/F3Gc8namCnzedhf7b1lrW+HK168DU9z5PNfcwwBmA9eb08epZnQS+NqdmxOAF4AOxjXcpAKuXtsncA09qAbcYIwZnGH5n8/Yzow9posy5OyHgFEZcvY37tlylbPd2xMDTALCMw73yINzeY+9WpEvXt17bftxHQY6pSqupAKuvckzfQhsAGpa1wk8L3D6Iej8EA2kkPm4yigyxOvubazEfzF7SjSuvdlT7ZYGyuRzG1twJewbM050b2Mf/jusc+oErUXGmP3Adlw/GqcNHciOtfYNa21Na+1luIYAJAHbcrm4H5BVwj01JjbHz8SpPe9Tf8A6oHOGaad6T7cCjU8tZ1xXEmhIhqEBxpjewESgm3WP7cqGj7uX6JSMZy6f9j7jGhsmUmAob+c5poKSt9NZa+Otta/jOmRe7xzbzS4Pb+bsz8Gpx7VxHdKf7u50+BvXkIdsr16QvhJrO2fI2W8Cz2fI2aeGWW0DShtjqmZYtDFZD+fywdUre1k2TWfsEMkyZ3voPS5Qinzx6jYN157TJcaY8sBT/HfW3wGgfIbxTOD6gB231sa6x8rcnd8BWWuTcR3ieMsYU8EY42uMae1OBuFAH2NMW3cBMhI4jKvw8aSvcO2dtjDG+ONK/mk5LJMn1toU4H/Ai8aYG909kxWBz3AlqonuXpdbgSdxHRo79TcQuM6cfUWCsxjXSRVXGJfqwHvAePce8JnzVjTG9HUv42eM6YnrigM/up9vZYyp5V5XeVzJbKH7UF9+WQSUdI8VK47rEFMcsMIdQ1dcQxWutdbmND4NXD/ag40xPu7tqQsEGdf4u6+Ah40xlxnXyS+P5eN2iOQX5e3cKRB5G8AY86gxpo1xXTKrmDHmHlzDEnI8XG6MaWiM6eJe1t/dS9oCWJ7FIpOBfu7l/HF9Pn505+XfgQB3XjfGNRa2b27iyC1r7VFgDvCCMaaUMaY9cA3uo4XGmG7GmBD3Z6QMrt+NKODPbFY7yL3ttXCdrBzk/ix5/D0uaFS8ujwLROLaI4rANRD81Di/TbgOLf9t/hvf8zBwlzEmFngHV1LyhOG4TgzaiCvJvQgY97icO4H3cZ2wdDWuosWjZ4S7i6IRuM7k3Ydrb+84rj3n/GznM1zbNxLXZUZODchv7T4U3g7XCQ/vuceV7bfW7sc1zmsfrvE/GNeZ8jdk0UwpXF/4WFzv9xJcY+ZwL/u8MWZmhvkfxJVYDuO6MsFQa+0P7ufq4LrMSSyuz88xzu4Bfsacfr3AbIcnZPKanOS/Kyscw3Uyx3UZ3vPncO1pL8rQxsxMV+aSiutkl4O4LrH1knu7quP6oVmJ6zuxFlcCztf3WCQfKG/nQgHK27jbnIBr5+Jf4A5ceWwvgDHmU2NMVmP1fXHlqIPuZe8BbrDWbnEv28kYcyhDPPNwvfYLcfXSB+M+YmddY/lvxDVs4BiwHleuezVDe+3N2dd5bZTHl+VuXENTDuE6We9Oa+2p4rQsrt+sE7gK1oq4hntl93lIxDX0ZJZ72zoD/S7Ue1yQmKzHOotkz7jG5hwBKlpro52OR3LHGFMP2GKtzdVlgowxfYCx1tq6no1MRDxNeds7GdfwuL7W2hW5mLfQv8fqeZU8Mcb0Nq4B74G4Btv/Wli/HEWVMSbIfXjO17jGaz2Na49eRLyQ8nbhV9TeYxWvklc34joEsxfXyQbZ3qZRvJIPMBbXYae1uE5yeSnbJUSkIFPeLvyK1HusYQMiIiIi4jXU8yoiIiIiXiNPxasxZoGnAhER8STlLxGRwiFXZxufUrp06a7NmzfXOAMR8UYnnA7gQitXrpytVq2a02GIiOTZ+vXrD1lrL83suTwVr7Vr12bdOk9fT1lEJP8ZY7K7+HehVK1aNeVsEfFKxpi/s3pOY15FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl5FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl5FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl5FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl7F49q3b88PP/xw2rQ333yTYcOGZTr/7t27adiwYbbr3L17N1OnTk1/vG7dOoYPHw7Ap59+yv333w/ApEmTmDJlyvmED0C1atU4dOjQea9HRKSgU86Wgk7Fq3jcgAEDmD59+mnTpk+fzoABA855nWcmwubNmzNhwoSz5hsyZAiDBg0653ZERIoa5Wwp6FS8isf17duXOXPmkJiYCLiSWFRUFK1bt2bEiBE0bNiQRo0aER4eftayu3fvpk2bNjRt2pSmTZuycuVKAEaOHMny5csJDQ3ljTfe4KeffqJnz55nLf/cc88xfvx4oqKiCA0NTf/z9fXl77//5uDBg9xwww20aNGCFi1a8MsvvwBw+PBhunTpQpMmTbj33nux1nrwFRIRKTiUs6WgU/EqHnfJJZcQFhbGggULANce/E033cSMGTOIiIhg06ZNLF68mBEjRhAdHX3asuXLl2fRokVs2LCB8PDw9MNMY8eOpU2bNkRERPDwww/nGEPFihWJiIggIiKCu+++mxtuuIHLL7+cBx98kIcffpi1a9fy7bffctdddwHw/PPP07p1azZu3Ejv3r35559/8vlVEREpmJSzpaDzczoAKRpOHYa69tprmT59Op988gmff/45AwYMwNfXlwoVKtCuXTvWrl1LSEhI+nLJycncf//9RERE4Ovry/bt288rjl9++YWPPvqI5cuXA7B48WIiIyPTnz9x4gQxMTEsW7aMGTNmANCjRw8uvvji82pXRMSbKGdLQabiVS6I6667jkceeYQNGzZw8uRJmjZtmqtB+W+88QYVKlRg06ZNpKWlUaJEiXOOITo6mjvvvJPvv/+ewMBAANLS0li1ahUlS5Y8a35jzDm3JSLizZSzpSDTsAG5IAIDA2nfvj2DBw9OH/Tftm1bwsPDSU1N5eDBgyxbtoywsLDTljt+/DjBwcH4+Pjw+eefk5qaCkBQUBAxMTG5bj85OZl+/frxyiuvUKdOnfTpXbp0YeLEiemPIyIi0mP78ssvAZg/fz5Hjx49tw0XEfFCytlSkKl4lQtmwIABbNq0if79+wPQp08fQkJCaNy4MR07dmTcuHFcdtllpy0zbNgwPvvsM1q2bMn27dsJCAgAICQkBD8/Pxo3bswbb7yRY9srV65k7dq1jBo1Kv0EgKioKCZMmMC6desICQmhfv36TJo0CYBRo0axbNkymjZtysKFC6latWo+vxoiIgWbcrYUVCYvZ+Q1b97crlu3zoPhiIh4hjFmvbW2udNxXEjK2SLirbLL2ep5FRERERGvoeJVRERERLyGilcRERER8RoqXsXrLViwgLp161KrVi3Gjh3rdDgiIpIN5Ww5XypexaulpqZy3333MX/+fCIjI5k2bdppF7AWEZGCQzlb8oOKV/Fqa9asoVatWtSoUQN/f3/69+/PrFmznA5LREQyoZwt+UHFq3i1ffv2UaVKlfTHlStXZt++fQ5GJCIiWVHOlvyg4lW8WmbXKdYtAkVECiblbMkPKl7Fq1WuXJk9e/akP967dy8VK1Z0MCIREcmKcrbkBxWv4tVatGjBn3/+ya5du0hKSmL69On07t3b6bBERCQTytmSH/ycDkDkfPj5+TFx4kS6du1KamoqgwcPpkGDBk6HJSIimVDOlvxgMht/khXdJ1tEvFV298kurJSzRcRbZZezNWxARERERLyGilcvtn//fvr370/NmjWpX78+3bt3Z/v27Xlez5tvvkl8fPw5xRAREcG8efPOaVkRkaJEOVskf6h49VLWWvr06UP79u3ZsWMHkZGRjB49mgMHDuR5XUqEIiKepZwtkn9UvHqppUuXUqxYMYYMGZI+LTQ0lNatWzNixAgaNmxIo0aNCA8PB+Cnn36iffv29O3bl3r16nHzzTdjrWXChAlERUXRoUMHOnToAMDQoUNp3rw5DRo0YNSoUenrX7t2LVdeeSWNGzcmLCyM48eP8+yzzxIeHk5oaGh6WyIicjrlbJF8ZK3N9V+zZs2sFAxvvfWWfeihh86a/s0339hOnTrZlJQUu3//flulShUbFRVlly5dakuXLm337NljU1NTbcuWLe3y5cuttdZefvnl9uDBg+nrOHz4sLXW2pSUFNuuXTu7adMmm5iYaKtXr27XrFljrbX2+PHjNjk52U6ePNned999F2CLRc4PsM7mId8Vhj/l7IJDOVskb7LL2ep5LWRWrFjBgAED8PX1pUKFCrRr1461a9cCEBYWRuXKlfHx8SE0NJTdu3dnuo6vvvqKpk2b0qRJE7Zu3UpkZCR//PEHwcHBtGjRAoDSpUvj56crrYmInA/lbJG8U/HqpRo0aMD69evPmm6zufRZ8eLF0//v6+tLSkrKWfPs2rWL8ePHs2TJEjZv3kyPHj1ISEjAWqtb+ImInCPlbJH8o+LVS3Xs2JHExEQ+/PDD9Glr167l4osvJjw8nNTUVA4ePMiyZcsICwvLdl1BQUHExMQAcOLECQICAihTpgwHDhxg/vz5ANSrV4+oqKj0HoGYmBhSUlJOW1ZERDKnnC2Sf1S8eiljDDNnzmTRokXUrFmTBg0a8NxzzzFw4EBCQkJo3LgxHTt2ZNy4cVx22WXZruuee+6hW7dudOjQgcaNG9OkSRMaNGjA4MGDueqqqwDw9/cnPDycBx54gMaNG9O5c2cSEhLo0KEDkZGRGvwvIpIN5WyR/KM7bIlIkaA7bImIeA/dYUtERERECgUVryIiIiLiNVS8ioiIiIjXUPEq53WrwYLAWsvWrVuzveSMiEhhURhy9rfffsu2bducDkW8lIpXOadEmJqa6qFo8iYxMZGBt95G49AmDLhlEImJiU6HJCLiUd6cs5OTkxl81z3cPvRBwlpdxYIFC5wOSbyQitdCZNy4cUyYMAGAhx9+mI4dOwKwZMkSbrnllkzvf53ZfbIXLlxIq1ataNq0KTfeeCOxsbEAVKtWjRdeeIHWrVvz9ddfO7CFZ3vy6WeY/uXnpKYkEz71C/735FNOhyQikitFMWevXLmSTz/5iNiD+4g5doRh9w93OiTxQipeC5G2bduyfPlyANatW0dsbCzJycmsWLGCNm3a8PLLL7Nu3To2b97Mzz//zObNmxk+fDgVK1Zk6dKlLF26lEOHDvHSSy+xePFiNmzYQPPmzXn99dfT2yhRogQrVqygf//+Tm3maR556EF+/PFHABYvXsyIRx52OCIRkdwpijm7bdu26b3GgwYN4vetvzkckXgj3ei4EGnWrBnr168nJiaG4sWL07RpU9atW8fy5cuZMGECX331FR988AEpKSlER0cTGRlJSEjIaetYvXo1kZGR6Re6TkpKolWrVunP33TTTRd0m3JSqVIlKlWqBECHDh3w8dH+mIh4h6KYs40xlCxZEnDd/jbjLXBFckvFayFSrFgxqlWrxuTJk7nyyisJCQlh6dKl7Nixg5IlSzJ+/Pj02xHefvvtJCQknLUOay2dO3dm2rRpmbYREBDg6c0QESkSlLNFzo26qQqZtm3bMn78eNq2bUubNm2YNGkSoaGhWd7/Gk6/T3bLli355Zdf+OuvvwCIj49n+/btjmyLiEhhp5wtkncqXguZNm3aEB0dTatWrahQoQIlSpSgTZs2Wd7/Gk6/T/all17Kp59+yoABAwgJCaFly5a6nImIiIcoZ4vkncnLtTF1n2wpqIwxpKamasyrZCm7+2QXVsrZUlAZY7j77rv54IMPnA5FCqjscrZ+6UVERETEa6h4FRERERGvoeK1kFmwYAF169alVq1ajB079qznrbUMHz6cWrVqERISwoYNG3Jc9rnnnqNSpUqEhoYSGhrKvHnzLsi2iIgURp7I00eOHKFz587Url2bzp07c/ToUQC+/PLL9NwdGhqKj48PERERnt9IEU+y1ub6r1mzZlYKrpSUFFujRg27Y8cOm5iYaENCQuzWrVtPm2fu3Ln2mmuusWlpaXbVqlU2LCwsx2VHjRplX3311Qu+PXkB2NTUVKfDkAIMWGfzkO8Kw59ydsHjqTw9YsQIO2bMGGuttWPGjLGPP/74WW1v3rzZVq9e3cNbmDuAvfvuu50OQwqw7HK2el4LkTVr1lCrVi1q1KiBv78//fv3Z9asWafNM2vWLAYNGoQxhpYtW3Ls2DGio6NztayIiJwfT+XpWbNmcdtttwFw22238d13353V9rRp0xgwYIDnN1LEw1S8FiL79u2jSpUq6Y8rV67Mvn37cjVPTstOnDiRkJAQBg8enH44SkRE8sZTefrAgQMEBwcDEBwczL///ntW2+Hh4SpepVBQ8VqI2Ewue2aMydU82S07dOhQduzYQUREBMHBwTz66KP5FLGISNHiqTydk19//ZVSpUrRsGHDXEYqUnDp9rCFSOXKldmzZ0/6471791KxYsVczZOUlJTlshUqVEiffvfdd9OzZ09PbYKISKHmyTwdHR1NcHAw0dHRlC9f/rR1Tp8+Xb2uUmio57UQadGiBX/++Se7du0iKSmJ6dOn07t379Pm6d27N1OmTMFay+rVqylTpgzBwcHZLhsdHZ2+/MyZM7XnLiJyjjyVp3v37s1nn30GwGeffca1116bvr60tDS+/vpr+vfvf+E2VMSD1PNaiPj5+TFx4kS6du1KamoqgwcPpkGDBkyaNAmAIUOG0L17d+bNm0etWrUoVaoUkydPznZZgMcff5yIiAiMMVSrVo3333/fsW0UEfFmnsrTI0eOpF+/fnz88cdUrVqVr7/+Or3NZcuWUblyZWrUqHHhN1jEA3R7WCkUdHtYyYluDytScOj2sJIT3R5WRERERAoFFa8iIiIi4jVUvIqIiIiI11DxKiIiIiJeQ8WriIiIiHgNFa8iIiIi4jVUvIqIiIiI11DxKiIiIiJeQ8VrEXbgwAGuveFGpk8PdzoUERHJhblz5/Ll1KlOhyHiKBWvRdRvv/1GSJNm/LhlH7cNvpOnnx1FXu62JiIiF461lhdfepkbb76NIQ+OYOh9w0lJSXE6LBFHqHgtopb/sgofXz9it6/C2DTmLVjodEgiIpKF+Ph4xr/2GiePHyb2UBSTJ39EdHS002GJOELFaxE1bMg9RO/ZDcCY0S+zYc0qjDHOBiUiIpkKCAjg+NEjVK5cGYCE+HiqVKnicFQizlDxKiIiIiJeQ8WriIiIiHgNFa8iIiIi4jVUvIqIiIiI11DxKiIiIiJeQ8WriIiIiHgNFa8iIiIi4jVUvIqIiIiI11DxKiIiIiJeQ8WriIiIiHgNFa8iIiIi4jVUvIqIiIiI11DxKiIiIiJeQ8WriIiIl7DWOh2CiONUvBZhp5JguXLlHI7k/MTFxQFgjHE4EhERzzpy5IjTIeSb0qVLOx2CeCkVr0XY8uXLAQgNDXU4kvPz/fffAypeRaTwO3nyJK1bt3Y6jPMSHx8PwNVXX+1wJOKtVLwWYd179ACgUaNGDkdyfgYOHOh0CCIiHrdy5UoAvv76a4cjOT933HEHAN26dXM4EvFWKl6LqG+++Ya42Fjefvttp0M5L4sWLQJg1apVDkciIuI5qampXHXVVQBcdtllDkdz7qKiovjqq6/o37+/06GIF1PxWsQkJSUxduxYbrzxRurUrcd9993ndEjnbN68eXTp0oWQkBBatmzpkTY2btzIgFsGkZiY6JH1i4jkZP/+/fj5+QFw6NAhh6M5dxs2bKBSpUoATJ061SNt9LzuBj76+GOPrFsKDj+nAyjo9u/fz+bNm8+afuLECW688cYLEoOvry8+Pj7nNabziiuuIHr/fv49cACAuvXqsWD+fPbs2ZNfYV4Qa9eu5bvvvuOLL74AXD0Qs2fP5p9//vFIe8uXL2f6l58DMPXzzzSuVsQLbNy4kbS0tLOmf/jhh/z222+5WocxBmNMeu7N+JfZtOymn4vExES2bdvGrl270qd99dVXbN269ZzW55TVq1czceLE035rfvzxx/RzLvLbyl9WsPCH+dSpXZu2bdt6pA1xnsnLZTeaN29u161b58FwCp6nnx3F6NEvU652Ewz/JaGkuOMc27PdwcjOnY+vHwEXeecVBmIO7z/tcdAlnj98VqJhZ5J2/sqIoXfw1BMjPd6eeIYxZr21trnTcVxIRTFn//XXX9SuXZvSFbNG+XUAACAASURBVKpQrETAac8d/nubQ1GdH+PjS7ka3nluwsG/ItL/X+riCgRcEuzR9qx/KUzdDiT9/D7r16ymZs2aHm1PPCe7nK2e1xyMfHwEX30zg6PBIZRqcX369JJAGefCKrLKOtSuf/2rGTP+cerXq0ufPn0cikJEclKrVi2GPfAg0xYsJ+DaZzG+//3MBToYV1F1uUPtpsUdoVPX7kSsX0OZMvq1Lmw05jUHgYGBLPlhHmm/zSX+z1+dDkcc4le6HEE9R3LrHXcSERGR8wIi4pgJb7xGk2rlifvpQ13Uv4gKbNKDmLL16HndDaSkpDgdjuQzFa+5UKVKFebP+Z74xRNJOrDT6XDEIcWD61Cy3T106daD6Ohop8MRkSz4+voy85twypzYSdyGOU6HIw4JaH8nW6OOM+yB4U6HIvlMxWsuhYWF8cF77xAzezSpsUedDkccEnBFG9LqdKRL916cPHnS6XBEJAulS5dmyQ/zSN04g5M7ita4X3ExPr4EdnuM8O/m8fbEd5wOR/KRitc8GDhwAMOH3k3snDGkJevSSUVVQMubiEorzcBbb9MhSZECrFq1asyZNZO4RRNIOvi30+GIA3xKBBLY+2meeGYUCxcudDocyScqXvPopReep22z+sQtnqjCpYgyxhDQ6X5+WreVZ0Y953Q4IpKNq666incnvEns7NGkxh93OhxxQLGLgwnsPoK+/QeybZt3XnFCTqfiNY+MMUz/YgoVfU4Qtzrc6XDEIT7FihPY8wneevcDpk2b7nQ4IpKNQYNu5d47biV2zlhsSrLT4YgDSlRpSPFWt9Cpa3cOHz7sdDhynlS8noOSJUuycN5sfP5cStzvnrnQshR8voEXE9T7Se4eOow1a9Y4HY6IZOOVMS/TqkEN4n58V0fNiqiARp1JqNSMbr2uIykpyelw5DyoeD1HwcHBLJw3h5M/f0BitHferEDOn3/5GpS6+n669eztdXcrEylKfHx8+Hr6l1ya/C9xa2c4HY44JKD1rfx1LI3Bd9+rnRgvpuL1PISGhvL55I+JmTOWlBPee79pOT+lav8fPo160OmaHsTGxjodjohkISAggMUL5sLW+cRvX+V0OOIA4+NL4DUPM+fHFYwb/5rT4cg5UvF6nvr06cPIRx8idvZo0pISnA5HHFKqeR8OF69I35sGZnpPdREpGCpVqsSCubOJ//Fdkg7scDoccYCPf0kCez3Ji2PGMXv2bKfDkXOg4jUfPPXESK5pG0bsD29grQqXosgYQ8DVQ1iz7R8ee/x/TocjItlo3rw5n3zwPjGzx5ASe8TpcMQBfqXLE9TjfwwcdDubN292OhzJIxWv+cAYw2effESNIEvcii+cDkccYnyLEdDjcT6aMo1PPpnsdDgiko1+/W7kkQeGEqfrdhdZxSvVo0SbO+ncrQcHDhxwOhzJAxWv+aR48eLMn/0dxf5ZTdyWJU6HIw7xLVWGwN5PMfyRx1ixYoXT4YhINp579hk6hIUQt+htnbxTRAXUb0dqzbZ07dGbhAQN/fMWKl7z0aWXXsriBfNI/OUzEvZudToccYh/uaoEdHmQXtddz86dO50OR0SyYIxh6pRPqVwsjrhV05wORxwS0Ko/e5JKcvOg27UT4yVUvOazBg0aED71C2Lnvkrysf1OhyMOKVmjGX5Nb+Dqrt05flx39REpqEqUKMHCebPx3bGMuMifnQ5HHGCMDwGdh7Pk1808/+JLTocjuaDi1QO6devGC88+Tezsl0lLjHc6HHFIqSY9iLmoNtdefyMpKSlOhyMiWahQoQKL5s8lYfnHJEb94XQ44gCfYsUJ7PUEr014h6+++trpcCQHKl495OGHhnN9t07Ezh+PTUt1OhxxgDGGgPZ3sXnPEe4f/rDT4YhINkJCQvjys8nu63YfdDoccYBfYFmCej3J4HvuZd26dU6HI9lQ8eohxhg+eO8drihfirhlnzodjjjE+PoR2H0E02Z+zzvvvut0OCKSjd69e/PU/x5zHTVLOul0OOIA/wo1KdVxKNf06MW+ffucDkeyoOLVg4oVK8acWTModWAzcZsWOB2OOMSnRCCBvZ/mf089y+LFi50OR0SyMfLxEfTocBVxum53kVWqzpXQ4Bo6XdODuLg4p8ORTPg5HUBmAgMDHbnN5rPPPkvbtm3p1KnTea+rWrVqBAUF4ePjQ9VKl7Fp1Zf4XBTM4Xlv4eNfEnxc+w1lOw+jROUrSDr4N0cWv09qzCHAEtCgI2Wu7I8xJn2d/377Iqnxxwi+Vbe08zbFLq5IYLdHuaFff9as+oW6devmetn9+/fz0EMPsXbtWooXL061atV48803qVOnjgcjztmxY8eYOnUqw4YNczQOKRgKU942xpCybw8xSz6gdKch7H1vcJ7zdmrcUQ7Pn+AagpCWim+Z8lS48fnzjlEujIAWN/Dvwr3c2P9m5syagY9P7vv6lLM9r0AWr0554YUX8nV9S5cupVy5cjz55JOUL1+ehfNfgzRDhQGj8S1VJn2+tOREDs54kbJdhlGyelPSkhM4OHMMsRvnEtS0p2uehFiSDuzAFCtB8rH9FLvosnyNVTyvRNUQUlsOpFPX7mzasJayZcvmuIy1lj59+nDbbbcxffp0ACIiIjhw4ECBSITvvvtuoUiE4r08lbcffvhhPvzoY+Iq1ALIc94+tvxLSlQLpXTzawFI+ndXvsYpnmWMIfDq+1g141lGPvk048aOztVyytkXRoEeNvDTTz/Rvn17+vbtS7169bj55pvTr8G2du1arrzySho3bkxYWBgxMTEkJCRwxx130KhRI5o0acLSpUsB+PTTT7nuuuvo1asX1atXZ+LEibz++us0adKEli1bcuSI6/aAt99+O9988w3g2gMfNWoUTZs2pVGjRmzbtg2AgwcP0rlzZ5o2bcq9997L5ZdfzqFDh7LdjrZt25KQkMD4V0aTFn+M1ITTeyfiIn+meKUrKFm9KQA+xUpQtvMQjq/+5r95/lhJyVphBFzRlvjfl+XDqytOCAjpysngUHr07kNycnKO8y9dupRixYoxZMiQ9GmhoaG0bt2aESNG0LBhQxo1akR4eDjg+s60a9eOfv36UadOHUaOHMmXX35JWFgYjRo1YscO173cb7/9doYMGUKbNm2oU6cOc+bMAcjyO7R161bCwsIIDQ0lJCSEP//8k5EjR7Jjxw5CQ0MZMWJEfr9U4qUKS97u2rUrTZuEkrhyCjY16aznc8rbqXFH8A0qlz6/f/nq5/OyigOMXzECe47kvU8+Y8qUz3O1jHL2hVGgi1eAjRs38uabbxIZGcnOnTv55ZdfSEpK4qabbuKtt95i06ZNLF68mJIlS/LOO+8A8NtvvzFt2jRuu+229DtmbNmyhalTp7JmzRqeeuopSpUqxcaNG2nVqhVTpkzJtO1y5cqxYcMGhg4dyvjx4wF4/vnn6dixIxs2bKBPnz78888/OW7DnDlzaNSoEUPuvZfAgFIc+OxBoj65n+gpjwCQfOhv/C+rddoyxS4OxiadTL/UVvzvPxNwRVsC6rcjTsWrVwtocxt/HEniznuG5HhB7C1bttCsWbOzps+YMYOIiIj0z/+IESOIjo4GYNOmTbz11lv89ttvfP7552zfvp01a9Zw11138fbbb6evY/fu3fz888/MnTuXIUOGkJCQkOV3aNKkSTz44INERESwbt06KleuzNixY6lZsyYRERG8+uqr+fgKibcrLHm7VatWfD19Kjb+BPu/fJyoyQ/kOm8HNenB4fkT2D/tCY6vDCcl5vA5v57iHN9SZQjq9RTDhj/EypUrc5xfOfvCKPDFa1hYGJUrV8bHx4fQ0FB2797NH3/8QXBwMC1atACgdOnS+Pn5sWLFCm699VYA6tWrx+WXX8727dsB6NChA0FBQVx66aWUKVOGXr16AdCoUSN2796dadvXX389AM2aNUufZ8WKFfTv3x+Aa665hosvvjjL2Dt06EBoaCgnTpzgiSeeAKBs2bJc2epKSl9en8tOG7tqMl8JkBp3lOSj0RSv3IBiZSthfHxJOph5zFLwGR9fArs+wveLfua11984p3WsWLGCAQMG4OvrS4UKFWjXrh1r164FoEWLFgQHB1O8eHFq1qxJly5dgLM/6/369cPHx4fatWtTo0YNtm3bluV3qFWrVowePZpXXnmFv//+m5IlS57fiyCFWmHK2127duWiiy+iRDEfLhs4luBBr2eYO+u8XbJGMyrd+xFBjbuSfGQv0Z8+SGq8bljijfwvvZyAzsPpeW2fLD93OVHOzl8FvngtXrx4+v99fX1JSUnBWnvaiUynZNeLlXE9Pj4+6Y99fHyyvID8qXlOtZtTG2daunQpERERTJkyhYsuuih9+pTJH1HmxC7iNswGoFi5qiTt//O0ZZOP7cf4l8SneCnifl9OWkIs+ybdyd73BpNy/ABxvy/PdRxS8PgUL0Vgr6d47uWxzJ07N8v5GjRowPr168+anh+f9TO/Q8aYLNc7cOBAvv/+e0qWLEnXrl358ccfs2xfpLDl7dJBQVzXtSOx8/67bndOeRvAt2QQAfXbU67noxQPrk3Cni25jkMKlpI1m+MTeh1Xd+3OiRMnspxPOfvCKPDFa2bq1atHVFRU+l5LTEwMKSkptG3bli+//BKA7du3888//+TprO7caN26NV999RUACxcu5OjRo3leR1BQEEt+mEfqxpmc3LGOgPrtSdwbycndEYDrRICji9+n9P/dAEDc7z9Tod/zVB76CZWHfsJlt72pca+FgF+Z8gR2f5z+N9/Kli2Z/6h17NiRxMREPvzww/Rpa9eu5eKLLyY8PJzU1FQOHjzIsmXLCAsLy1P7X3/9NWlpaezYsYOdO3dSt27dLL9DO3fupEaNGgwfPpzevXuzefNmgoKCiImJOfcXQIoUb8/br74ylgYVSxP38ycAOebtk39vIi3ZNfwhLTGe5KPR+JW+ND82RxwS0LQXx0tXp0/fm0hNzfzmQ8rZF4ZXFq/+/v6Eh4fzwAMP0LhxYzp37kxCQgLDhg0jNTWVRo0acdNNN/Hpp5+etkeTH0aNGsXChQtp2rQp8+fPJzg4mKCgoDyvp1q1asyZNZO4RRNIORbNpTc8w/GV09n34b1Ef3I//pfVJqhpT1KOHyD1xEH8K9ZLX7bYRZfhU7yUbmNYCJSofAUl2g6m0zXd+ffff8963hjDzJkzWbRoETVr1qRBgwY899xzDBw4kJCQEBo3bkzHjh0ZN24cl12WtytQ1K1bl3bt2tGtWzcmTZpEiRIlsvwOhYeH07BhQ0JDQ9m2bRuDBg3ikksu4aqrrqJhw4ZeP/hfPM/b87afnx+zZ35D4KFI4iLm41OseJZ5GyBp/1/s/+xhoj65n/1fPEZQ4y4UD3b2bHM5P667Jt7Dxt3/8uAjj2U5j3K255m8HE5p3ry5Leq3TEtMTMTX1xc/Pz9WrVrF0KFDiYiIOOf1TZnyOfc9NpKL+r962mVYpGiJXfEFVRJ3sXrFz/n+w52Z22+/nZ49e9K3b1+Pt1VQGGPWW2ubOx3HhaSc7ZKfefuvv/6i+f+1onjnhyhZLTSfIxVvkJoQy/Hpj/PqC09x7z33XJA2lbNP55U9r076559/aNGiBY0bN2b48OGnHRo4F4MG3cqQOwYRO2csNiXnSydJ4RRw1UB2n7A885wuYi6S3/Izb9eqVYtZM74hbsHrJB/em49RirfwLRFIUO+nePTxJ9IvTSUXlnpeC4AtW7bQrHkYl971Pn6BOV+4Xgqf5GP7OR4+km+nf0HXrl2dDqdQUs+r5BdrLV2u6c6GuDIEtb7V6XDEIcd+mEi3huWZ9vlnTodSKKnntQA7ePAgXbr3pHTnoSpci6i0xHhiZ7/Mi6OeOa/CdcGCBdStW5datWoxduzYs54/evQoffr0ISQkhLCwsNNOEnvjjTdo0KABDRs2ZMCAAenX2RSRs7351tus2fIHAS1ucDoUcUhc5FJKHNzKW6959/VSvZWK1xwMHjyY8uXL07Bhw3xfd2JiItf0vJakqi0JaNAx39cvBZ9NSyV2/nj6du/MQw8+cM7rSU1N5b777mP+/PlERkYybdo0IiMjT5tn9OjRhIaGsnnzZqZMmcKDDz4IwL59+5gwYQLr1q1jy5YtpKampt/WUMTbeDJnA8yfP59nnn+RwF5PpV8SS4qWxH2/k7DsExbNn0v58uU91o6nP8veTMVrDm6//XYWLFiQ7+u11nLb4LvYHedLwFU35/v6xTvELZvMFRUCmPTuxEyvgZlba9asoVatWtSoUQN/f3/69+/PrFmzTpsnMjKSq6++GnBdtmj37t0cOHAAgJSUFE6ePElKSgrx8fFUrFjx3DdKxEGeytng+g7dNPAWAns8TrGL8namuBQOKcf/JWbuOKZ9MYVGjRp5tC1Pfpa9nYrXHLRt25ayZfP/cP7LY8ayYPkaAro8hDF6G4qiuE0LCPh3C3NnzaBYsWLnta59+/ZRpUqV9MeVK1dm3759p83TuHFjZsyYAbiK3b///pu9e/dSqVIlHnvsMapWrUpwcDBlypRJv8OLiLfxVM4+dOgQna7pTvGrbqNE5fr5vn4p+E4N8Rr15P/o2bOnx9vz1Ge5MFDV5ICZM2cy9rU3Cez1JD7+JZwORxxw8u9NJP06jSUL559297VzldmJl2f25I4cOZKjR48SGhrK22+/TZMmTfDz8+Po0aPMmjWLXbt2ERUVRVxcHF988cV5xyRSWKQP8aryfwQ0vNrpcMQBNi2V2B9ep3fndjz26CNOh1Pk+TkdQFGzceNGbr3jToKufRa/oHJOhyMOSD6yj7j5rzHnu2+pVatWvqyzcuXK7NmzJ/3x3r17zzr0X7p0aSZPngy4it3q1atTvXp1fvjhB6pXr86ll7ru/nP99dezcuVKbrnllnyJTcSbWWu5/c672RXrQ2APfSeKqrjlU6hT1p+PP5h0XkO8JH+o5/UCio6Opkv3npRsfy/Fg2s7HY44IPVkDDHfv8Rr48bQvn37fFtvixYt+PPPP9m1axdJSUlMnz6d3r17nzbPsWPHSEpKAuCjjz6ibdu2lC5dmqpVq7J69Wri4+Ox1rJkyRKuuOKKfItNxJuNHvsK85f9qiFeRVjc5oWUiN7AvO9nnvcQL8kf6nm9QE6ePEmX7r2wdToSWK+10+GIA2xqCnHzxjHophvy/a4sfn5+TJw4ka5du5KamsrgwYNp0KABkyZNAmDIkCH8/vvvDBo0CF9fX+rXr8/HH38MwP/93//Rt29fmjZtip+fH02aNOGeC3TXGJGC7LvvvmPM+De46KZxGuJVRCX8s5mk1V+yYtUvGn9agOgmBTkYMGAAP/30E4cOHaJChQo8//zz3HnnnXlah7WW6/r2Y8VfRwi85mEdciiCrLXELnmPJpdYfpg3G19fX6dDKnJ0k4KiIT9yNkBERASt23ck6NpnKB5cxwORSkGXfDSKY+Ej+e7r6XTu3PmCt59fn2VvlV3OVs9rDqZNm3be63j62VEs2/A7pW94UYVrERW/YTYXxexi5tJfVbiKeFB+5Oz9+/e7hni1u0eFaxGVlhBL7PcvMW70i44UrpA/n+XCSsWrh02dOo0Jkz7iov7jMH7+TocjDji5Yx2pEd+xZO2vBAUFOR2OiGTj5MmTdO7Wk7TaHQi8oo3T4YgDbFoqsfPGM6BPb+4bNszpcCQTKl49aM2aNdwz7H5KX/88vgEXOx2OOCDp4G7iFk1g4bzZVKtWzelwRCQb1loG3HobUWmlCWx5k9PhiEPifvqIkCoXM3HCG06HIllQ8eohe/bsoVvP3pTqdB/+5as7HY44IDXuGDHfj+a9t9/iyiuvdDocEcnBs889z8/rthKkIV5FVtyGOQQd3c6sRb/i56cSqaDSO+MBsbGxXN21Oz6NelKq1v85HY44wKYkEzv3FYbeeRu33qprQ4oUdNOnh/PmO+9zUf9X8SlW3OlwxAEnd64nef23LFm7mjJlyjgdjmRDF63LZ2lpadzQbwCHS1aiVPPrnA5HHOC6ssA7XNmwJmNHv+R0OCKSgzVr1nDXkGEE9X4S30AN8SqKkg/tIW7hW8z+7ltq1KjhdDiSAxWv+ezREY+z5o89BHYcosNORVTc2m8pn3KQr6Z9gY+PvmIiBdmePXvo1utaSl19H/7lVbQURanxx4mZ/TJvv/EabdroJD1voGED+ejjTz7h4y/CKXPTKxhf3YWjKIrfvhIif2Dx+rUEBAQ4HY6IZCM2NpZO1/TAp2F3StXWEK+iyKYmEzd3HHfd2p877rjd6XAkl1S85pPly5fz4CMjKNP3ZXxLaaxMUZR0YAfxP77Hz0sWUalSJafDEZFspKWl0femgRwuXpGA5n2cDkccYK0lbskkwupVZfy4V5wOR/JAxzTzwc6dO+nd5wYCuj5EsXJVnA5HHJASc5iY2aP55IP3adasmdPhiEgORvxvJGu2/UPA1RriVVTFr5vJJYlRfBM+VUO8vIx6Xs/T8ePHubpLN/ya9aVk9aZOhyMOSEtOIHbOGB4dfh/9+t3odDgikoPJkz/lwynTNMSrCIv/81fSfpvL4vVrCQwMdDocySMVr+chJSWFXn36cqJsXYKa9HA6HHGAtWnELZxAp5ahjHrmaafDEZEcrFixggceflRDvIqwpH93Er94Ij8uWkCVKjpa6o3UT34e7hv+IFv2HiOw/V1OhyIOiVs1nSrFE/jis8k69ChSwO3cuZNe111PQJcHNcSriEqNPUrM96P5cNK7hIWFOR2OnCP1vJ6jie+8w/SZc12HnXx8nQ5HHBAX+TO+O5bxw4Z1lChRwulwRCQbJ06c4Oqu3fFregMla2hcelGUlpxI7JwxPDjsHgYM6O90OHIeVLyeg0WLFjHy6VGU6TcGnxIaK1MUJe7bRsLyj1m57CcqVKjgdDgiko2UlBR6X38jMRfXIbBpT6fDEQdYa4lbPJH2zRvw4vPPOR2OnCcNG8ijP/74g743DSCw22MUu7ii0+GIA1JO/EvM3LFMnfIpjRo1cjocEcnBAw8+wuZ/jhCgIV5FVtzqcCr6nGDq559piFchoJ7XPDh8+DBXd+2Gf8ubKVFVRUtRlJZ0ktjvR/PsE/+jV69eTocjIjl49733mDpjFmVuGqchXkVU3O/L8flzKQvXr6VkyZJOhyP5QMVrLiUlJdG9dx8SgpsSGNLF6XDEATYtlbgFb9Dz6taMeOxRp8MRkRwsXryYx598RkO8irDE6O2c/PkDVvz0I8HBwU6HI/lEwwZywVrLXfcO5c8jyQS0GeR0OOKQuBWfU/MiHz758H0ddhIp4P744w9u6NefwG6PaohXEZVy4hAxc8byxaefEBoa6nQ4ko/U85oL4197ne8XLaP0jaN12CkP4revIm7bivxboU0jftty1/99/OCM+jGwQUdKeOAs4uIV65G4ewMl9q1n/oa1+Pv753sbIpJ/jhw5Qqeu3d1DvEKcDkcckJaUQOyc0Tzx2MNcd911Tocj+UzFaw7mzJnDE08+yUVdHyDp353/PWHTSD6yL+cVWAvWYm0a2DT34zSse/qpaRmftxmmu+Z1/z8Xghpfg0+p0mdNT40/jo9/KYzfhbubTMyGOST8vckzK09LOWtS7OaFtKpcPF+bid63lx2rpuGXFMPylSu45JJL8nX9IpK/kpOT6XHt9cRdWp/Aem1JS0rwQCs2Q05253j4L6ef+n8uGD9/fIqX8kCMeWPTUnP3m5ZHaQkx2OQkjO8Z5YaPL8XKVsr39oyvH8a/JLE/vEH3di15cuT/8r0NcZ6K1xwcOX6CBo2bwP4VsP+/6Qei9nEkaq9zgWUhZt0sgi4qe/b0Y0dOexxY5mKPH/ouBhTLJJbzkZSYgK9/SWzpy/AljYC0WH5cuID69evnazunTJ06lVsHDeK7efOoV6+eR9oQkfwTFRXFpg3rSE1LI27LktOeS0o46VBU2StW/OzrRCcnnl50G2Pw88/fnfPs2rsQMvutOl8JJ09S6oq21AqyTP7oAw3xKqSMzeXeIUDz5s3tunXrPBiOeNLevXvp2rUrkZGRAMTExHjlPZ0jIyN5a+I7fPHlVChfB/+ju1i/ZjU1atTI97bi4+P5/fffadZMFzX3dsaY9dba5k7HcSEpZ5/uxIkTJCUl5TifMSZXf56UnJzMtm3beOKJJ1i6dCkAM2bMoGvXrh5tN7+lpaURFRXFgAED2bBhPQB//fUXNWvW9Eh7l1WphsHy28b1lCtXziNtyIWRXc5W8VoEHTlyJP3wd1pamtfumUZERNChUxcSS1ehTlk/Itb96nRIUoCpeBVvFhoayqZNm1i1ahUtW7Z0OpxzkpiYmH43woMHD3qkuPzgo49pfWUrjx2Nkwsnu5ytqw0UQWXLluXw4cMAjBkzxuFozl1oaCjfhE/DNyaaHbv+ZsOGDU6HJCLiEREREQC0atXK4UjOXfHixdN7vht76Oz/e+66U4VrEaDitYgqW7YsNWvW5KmnnnI6lPPSsWNHKpYvR0pAOV5/622nwxER8Zjdu3cDsGmTh06EvQCKFSvGuHHjiNq3L1dDOEQyo+K1CFuyxHUyQ2pqqsORnDtjDEPvHoyfjw+Lf1zqdDgiIh5z+eWXA9C/f3+HIzk/jzzyCAAffPCBw5GIt1LxWoSdSoTR0dEOR3J+qlatip8PxJ447nQoIiIe5e/vz7Zt25wO47z4+rqul+7NPcjiLBWv4vWHbvz8/DDGYNPSnA5FRMSjAgICnA4h35w8WTAvXSYFn4pXEREREfEaKl5FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl5FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl5FRERExGuoeBURERERr6HiVURERES8hopXEREREfEaKl7Fqy1cuJAlS5YAkJKcxGeffcaOfdhpfgAAFfRJREFUHTscjkpERDKTkJDAunXrANi1axe///67wxGJN1LxKl7t4cceZ+J7k/AP6U6p2i25/fbb+W7WbKfDEhGRTCxZ8iMtWrQAYOXKlVzZug2pqakORyXeRsWreLVlS5fQvMX/kbh1Ecn7tjLmlXE88vCDToclIiKZ6N69Gy+PGUvJoIu4pEIwixf+gK+vr9NhiZfxczoAkfNxySWXsHzpYp5/6WVahbWgZ8+eTockIiJZMMbw5Mj/0aXT1VSqVIng4GCnQxIvpOJVvJ6/vz8vv/C802GIiEguNW/e3OkQxItp2ICIiIiIeA0Vr3LB7d+/n/79+1OzZk3q169P9+7d2b59e6bz7t69m4YNGwLw008/nfOwgDfffJP4+Pg8LXM+7YmIFBbK2VLQqHiVC8paS58+fWjfvj07duwgMjKS0aNHc+DAAY+2ey6JUESkqFPOloJIxatcUEuXLqVYsWIMGTIkfVpoaCitW7dmxIgRNGzYkEaNGhEeHp7teuLi4hg8eDAtWrSgSZMmzJo1C4DU1FQee+wxGjVqREhICG+//TYTJkwgKiqKDh060KFDB8B1fdhWrVrRtGlTbrzxRmJjYwFYsGAB9erVo3Xr1syYMcNDr4KIiHdQzpaCSCdsyQW1ZcsWmjVrdtb0GTNmEBERwaZNmzh06BAtWrSgbdu2Wa7n5ZdfpmPHjnzyySccO3aMsLAwOnXqxJQpU9i1axcbN27Ez8+PI0eOULZsWV5//XWWLl1KuXLlOHToEC+99BKLFy8mICCAV155hddff53HH3+cu+++mx9//JFatWpx0003efKlEBEp8JSzpSBS8SoFwooVKxgwYAC+vr5UqFCBdu3asXbtWkJCQjKdf+HChXz//feMHz8ecN215Z9//mHx4sUMGTIEPz/XR7ts2bJnLbt69WoiIyO56qqrAEhKSqJVq1Zs27aN6tWrU7t2bQBuueUWPvjgA09sroiIV1POFiepeJULqkGDBnzzzTdnTbfW5mk99v/bu/egqOr3D+DvxQsYKd9wvGCYtoCgu+wu6CKbsgimZBZNVy9kMmk5adpltKb6mY4/a7wwTWql0w1qMiUtf1qmqYiKTMXFVsSVkE3JkBRLTVJgWZ/fH07rF4FNBFzO7vs1szPu2f2cfc7+8fbhcDiPCL788kuEh4c32q5Sqf517ZgxY7Bu3boG2y0Wy7+uJSLyJsxs6oh4zSvdVImJiaitrcUHH3zg3Jafn4/bbrsNmZmZcDgcqKqqwr59+xATE9PsfpKSkrBq1SpngP70008AgLFjx2LNmjWor68HAPz5558AgO7du+PChQsAgNjYWOTm5qKsrAwAcPHiRZSWliIiIgLHjh2DzWYDgEZBSUTkbZjZ1BGxeaWbSqVSYdOmTdi5cydCQkKg0WiwcOFCTJ48GTqdDnq9HomJiVi2bBn69u3b7H7mz58Pu90OnU4HrVaL+fPnAwCmT5+OO+64w7mvzz//HADw9NNPY9y4cUhISECvXr2QkZGBSZMmQafTITY2FiUlJfDz88P777+P8ePHY+TIkRgwYMBN+U6IiDoqZjZ1RKqWnPofNmyYFBQUtGM5dLOpVCrYbDao1Wp3l0LUrlQqVaGIeNVYH2a25wkMDMTZs2db/Gt7IqVxldk880pEREREisHmlYiIiIgUg80rERERESkGm1dyu5bMzXalNeMELRYLvv322xtaS0TkTZjZ5G5sXsmt2nJuNoOQiKh9MbOpI2DzSm7V0rnZe/bswahRo/DII48gIiICKSkpEJEmZ2E/88wzGDZsGDQaDRYsWODcf35+Pu666y7o9XrExMTg/PnzeP3115GZmQmDwfCvM7qJiLwVM5s6BBG57sfQoUOFPAsAsdlsbvv8FStWyPPPP99o+8aNG+Xuu++W+vp6+f3336V///5y8uRJyc7Olh49esiJEyfE4XBIbGys5OTkiIjIgAEDpKqqyrmPP/74Q0RE6uvrJT4+Xg4ePCi1tbVy5513Sl5enoiInD9/Xux2u6Snp8usWbNuwhGTuwAokBbknSc8mNme57bbbpMr/3W7BzObbhZXmc0zr9QhNTc3GwBiYmIQHBwMHx8fGAwGHD9+vMl9fPHFF4iOjkZUVBQOHz4Mq9WKn3/+GUFBQTAajQCAHj16OGdqExHRjWFm083E5pXcSqPRoLCwsNF2cXEDbl9fX+e/O3Xq5Bwr+N+OHTuGtLQ0ZGVloaioCOPHj0dNTc11zdEmIqKmMbOpI2DzSm7VVnOzgYazsP/66y/4+/sjICAAp06dwrZt2wAAEREROHnypPOMwIULF1BfX99gLRERNY2ZTR0Bm1dyq7aamw00nIWt1+sRFRUFjUaDJ598EiNGjAAAdO3aFZmZmZg9ezb0ej3GjBmDmpoaJCQkwGq18uJ/IiIXmNnUEahcneq/Fudkex6VSgWbzQa1Wu3uUojalas52Z6Kme15AgMDcfbsWZe/pifyBK4ym2deiYiIiEgx2Lx6qeLiYqSlpQEAVq1ahdzcXDdXREREruTn5+Ps2bMAgO+//97N1RC5D5tXL7U3Jxfz5s0DcGXKycb/2+LmioiIqDl///034szxzucj4+Jw+vRpN1bUeq2ZsEXejc2rl5r1zAykp2egq68fFv3vYry1bIm7S2o1BiEReSp/f39k7dqJ//Tshe7/CcS2b79F79693V1Wq9xIZjscjnaqhpSEzasXS02ditqaS5j/P695xH30GIRE5MlGjBiBksOHcKS4CGPHjnV3OU7Lli3DypUrAQAvvPACEhMTAQBZWVl4/PHHmxz72tR42B07dsBkMiE6OhqPPvooqqurAQADBw7EokWLMHLkSGzYsMENR0gdDZtX6nAYhERETevTpw9uv/12d5fRgNlsRk5ODgCgoKAA1dXVsNvt2L9/P+Li4vDGG2+goKAARUVF2Lt3L4qKijBnzhz069cP2dnZyM7OxpkzZ7B48WLs2rULBw4cwLBhw/DWW285P8PPzw/79+/HxIkT3XWY1IGweaUOh0FIRKQcQ4cORWFhIS5cuABfX1+YTCYUFBQgJycHcXFxTY59vdYPP/wAq9WKESNGwGAw4JNPPkF5ebnz9QkTJtzMQ6IOjgOCqcO5Ngijo6OdQbhy5Up88cUXeP/991FfX4/KykpYrVbodLoG+/jvIASAuro6mEwm5+sMQiKittGlSxcMHDgQ6enpuOuuu6DT6ZCdnQ2bzYZu3bohLS3NOYUrNTUVNTU1jfYhIhgzZgzWrVvX5Gf4+/u392GQgvDMK3U41wZhXFxcoyC8dv71tf4JQovFAovFAqvVio8++sj5OoOQiKjtmM1mpKWlwWw2Iy4uDmvWrIHBYGh27CvQcDxsbGwscnNzUVZWBgC4ePEiSktL3XIs1PGxeaUOiUFIRKQccXFxqKyshMlkQp8+feDn54e4uLhmx74CDcfD9urVCxkZGZg0aRJ0Oh1iY2NRUlLixiOijoyXDVCH9M+1rSaTCf7+/k0GoVqtbjIIg4KCkJ2d7QzC2tpaAMDixYsxaNAgdx0SEZHHGj16NOx2u/P5f58syMjIaHLN7NmzMXv2bOfzxMRE5OfnN3rf8ePH26xO8gyqlsxH5pxsIlIqV3OyPRUzm4iUylVm87IBIiIiIlIMNq8K8/vvv2PixIkICQnBkCFDcO+993aIaznPnTuH9957z91lEBF1OMxtorbF5lVBRAQPPvggRo0aBZvNBqvVijfffBOnTp1yd2ntGoLbt29HeHg4QkNDsWRJ4zG2IoI5c+YgNDQUOp0OBw4cAACcOHECCQkJGDx4MDQaDVasWNEu9RERNcdbcvtGc9rV2nnz5iEiIgI6nQ4PPvggzp071ya1kgcQket+DB06VMh9srKyJC4urtH2y5cvy9y5c0Wj0YhWq5X169eLiEh2draYzWZ59NFHJSwsTF5++WX57LPPxGg0ilarlbKyMhERmTp1qsyYMUNGjhwpYWFh8vXXX4uIyKVLlyQ1NVW0Wq0YDAbZvXu3iIgUFxeL0WgUvV4vkZGRUlpaKhMmTBA/Pz/R6/Uyd+7cNjvm+vp6UavVYrPZpLa2VnQ6nRw+fLjBe7Zu3Sr33HOPXL58Wb7//nuJiYkREZGTJ09KYWGhiIj89ddfEhYW1mgteQ8ABdKCvPOEBzPb/bwht1uT067Wfvfdd2K320VE5KWXXpKXXnrphmsk5XGV2bzbgIIUFxdj6NChjbZ/9dVXsFgsOHjwIM6cOQOj0Qiz2QwAOHjwII4cOYLAwECo1WpMnz4deXl5WLFiBVatWoW3334bwJW/5ty7dy9sNhsSEhJQVlaGd999FwBw6NAhlJSUYOzYsSgtLcWaNWvw3HPPISUlBXV1dXA4HFiyZAmKi4thsVja9Jjz8vIQGhoKtVoNAJg4cSI2b96MIUOGON+zefNmPPHEE1CpVIiNjcW5c+dQWVmJoKAgBAUFAbhyG63BgwejoqKiwVoiovbkDbndmpw+fvx4s2vHjh3rXB8bG4uNGze2qk7yHLxswAPs378fkyZNQqdOndCnTx/Ex8c7bzdiNBoRFBQEX19fhISEOMMgMjKywe1HHnvsMfj4+CAsLAxqtRolJSXYv38/pkyZAgCIiIjAgAEDUFpaCpPJhDfffBNLly5FeXk5unXr1m7HVlFRgf79+zufBwcHo6KiosXvOX78OH766ScMHz683WolIrpenpTbrcnp61kLAB9//DHGjRvXZjWTsrF5VRCNRoPCwsJG28XF7c58fX2d//bx8XE+9/HxQX19vfM1lUrVYJ1KpWp2v5MnT8aWLVvQrVs3JCUlYffu3S06jpZoqoZra/2391RXV+Phhx/G22+/jR49erR9kUREzfCG3G5NTl/P2jfeeAOdO3dGSkpKKyslT8HmVUESExNRW1uLDz74wLntn3nRmZmZcDgcqKqqwr59+xATE9OifW/YsAGXL1+GzWbDL7/8gvDwcJjNZqxduxbAlRtO//rrrwgPD8cvv/wCtVqNOXPmIDk5GUVFRQ2mW7Wl4OBgnDhxwvn8t99+Q79+/a77PXa7HQ8//DBSUlLw0EMPtXl9RESueENutyan/23tJ598gm+++QZr165t1NSS92LzqiAqlQqbNm3Czp07ERISAo1Gg4ULF2Ly5MnQ6XTQ6/VITEzEsmXL0Ldv3xbtOzw8HPHx8Rg3bhzWrFkDPz8/zJw5Ew6HA5GRkZgwYQIyMjLg6+uLzMxMaLVaGAwGlJSU4IknnkDPnj0xYsQIaLVazJs3r82O2Wg04ujRozh27Bjq6uqwfv16JCcnN3hPcnIyPv30U4gIfvjhBwQEBCAoKAgigmnTpmHw4MF48cUX26wmIqLr5Q253ZqcdrV2+/btWLp0KbZs2YJbbrnlhusjD9TcX3I19eBfrnqmqVOnyoYNG9xdRrO2bt0qYWFholarZfHixSIisnr1alm9erWIXPmr3ZkzZ4parRatViv5+fkiIpKTkyMAJDIyUvR6vej1etm6davbjoPcC7zbAHmQjpbbN5rTza0VEQkJCZHg4GBnfs+YMePmHhS5lavM5nhYQmpqKu677z488sgj7i6FqN1wPCx5EuY2eTpXmc1bZREyMjLcXQIREbUAc5u8Ga95VYhOnTrBYDBAq9Xi/vvvv+FJI9OnT4fVam20PSMjA88+++wN13frrbfe8FoiIk/DzCZqP2xeFaJbt26wWCwoLi5GYGCg80bULfXhhx/yJv1ERO2MmU3Ufti8KpDJZGpwE+fly5fDaDRCp9NhwYIFAIC///4b48ePh16vh1arRWZmJgBg1KhR+OcauPT0dAwaNAjx8fHIzc117i81NbXBJJN/fkKvrq7G6NGjER0djcjISGzevLlRbZWVlTCbzc4zDjk5OW3/BRARKQgzm6ht8ZpXhXE4HMjKysK0adMAADt27MDRo0eRl5cHEUFycjL27duHqqoq9OvXD1u3bgUAnD9/vsF+KisrsWDBAhQWFiIgIAAJCQmIiopy+dl+fn7YtGkTevTogTNnziA2NhbJyckN7r33+eefIykpCa+99hocDgcuXrzYxt8AEZFyMLOJ2h7PvCrEpUuXYDAY0LNnT/z5558YM2YMgCtBuGPHDkRFRSE6OholJSU4evQoIiMjsWvXLrz88svIyclBQEBAg/39+OOPGDVqFHr16oWuXbtiwoQJ/1qDiODVV1+FTqfD3XffjYqKCpw6darBe4xGI9LT07Fw4UIcOnQI3bt3b7svgYhIIZjZRO2HzatC/HP9VHl5Oerq6pzXT4kIXnnlFVgsFlgsFpSVlWHatGkYNGgQCgsLERkZiVdeeQWLFi1qtM/mppV07twZly9fdu6/rq4OALB27VpUVVWhsLAQFosFffr0QU1NTYO1ZrMZ+/btw+23344pU6bg008/bcuvgYhIEZjZRO2HzavCBAQEYOXKlUhLS4PdbkdSUhI+/vhjVFdXAwAqKipw+vRpnDx5Erfccgsef/xxzJ07FwcOHGiwn+HDh2PPnj34448/YLfbsWHDBudrAwcOdM7i3rx5M+x2O4Arv8bq3bs3unTpguzsbJSXlzeqr7y8HL1798ZTTz2FadOmNfpcIiJvwswmanu85lWBoqKioNfrsX79ekyZMgVHjhyByWQCcOVC/c8++wxlZWWYN28efHx80KVLF6xevbrBPoKCgrBw4UKYTCYEBQUhOjoaDocDAPDUU0/hgQceQExMDEaPHg1/f38AQEpKCu6//34MGzYMBoMBERERjWrbs2cPli9fji5duuDWW2/lT/FE5PWY2URtixO2iMgrcMIWEZFyuMpsXjZARERERIrB5pWIiIiIFIPNKxEREREpBptXL7Z9+3aEh4cjNDQUS5YscXc5RETkAjOb6Ao2r17K4XBg1qxZ2LZtG6xWK9atWwer1erusoiIqAnMbKKr2Lx6qby8PISGhkKtVqNr166YOHFik3OviYjI/ZjZRFexefVSFRUV6N+/v/N5cHAwKioq3FgRERE1h5lNdBWbVy/V1P19mxs9SERE7sXMJrqKzauXCg4OxokTJ5zPf/vtN/Tr18+NFRERUXOY2URXsXn1UkajEUePHsWxY8dQV1eH9evXIzk52d1lERFRE5jZRFd1dncB5B6dO3fGO++8g6SkJDgcDjz55JPQaDTuLouIiJrAzCa6StXUdTTN4ZxsIlIqV3OyPRUzm4iUylVm87IBIiIiIlIMNq9EREREpBhsXomIiIhIMdi8EhEREZFisHklIiIiIsVg80pEREREisHmlYiIiIgUg80rERERESkGm1ciIiIiUgw2r0RERESkGGxeiYiIiEgx2LwSERERkWKoROT636xSVQEob79yiIjazQAR6eXuIm4mZjYRKVizmd2i5pWIiIiIyJ142QARERERKQabVyIiIiJSDDavRERERKQYbF6JiIiISDHYvBIRERGRYrB5JSIiIiLFYPNKRERERIrB5pWIiIiIFIPNKxEREREpxv8D6abkMC3nD2cAAAAASUVORK5CYII=\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": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEGCAYAAABCa2PoAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3gU5fbA8e9JQui9Q5BQQgkdAtJBBUUFUawoipWr0m1X/F3r9dr1qoAggl6xYAVBFFGRDkqH0AkdpErvJJzfHzO5d42bzSbZzW6S83meediZeWfmbODhZN4qqooxxhiTXRGhDsAYY0zeYAnFGGNMQFhCMcYYExCWUIwxxgSEJRRjjDEBERXqAEKlXLlyGhsbG+owjDEmV1m6dOlBVS3v7Vy+TSixsbEsWbIk1GEYY0yuIiLb0ztnVV7GGGMCwhKKMcaYgLCEYowxJiAsoRhjjAkISyjGGGMCImwSioi8LyL7RWR1OudFRN4WkSQRWSUizT3OdRORDe65x3MuamOMManCJqEA/wG6+Th/JRDnbv2AUQAiEgmMdM/HA71FJD6okRpjjPmLsEkoqjoHOOSjSE9gvDp+BUqJSGWgFZCkqltU9RzwmVs2KA4cP8s/p67lzPmUYD3CGGNypbBJKH6oCuz02N/lHkvv+F+ISD8RWSIiSw4cOJClIOZuOsC4eVu5Zcyv7D92Jkv3MMaYvCg3JRTxckx9HP/rQdUxqpqgqgnly3udOSBDvZrHMLpPCzbsPU7PkfNZvftolu5jjDF5TW5KKLuAah77McDvPo4HTbeGlfjqgTYIcMPoBXyfuCeYjzPGmFwhNyWUKcAdbm+v1sBRVd0DLAbiRKSGiEQDt7hlg6pBlZJMHtCe+MolePCTZbz18yZsOWVjTH4WNpNDisgEoDNQTkR2AU8DBQBUdTTwPXAVkAScAu5yzyWLyABgOhAJvK+qa3Ii5vLFCzKhX2uGTUzk3z9vZOP+47x2QxMKR0fmxOONMSasSH79rTohIUEDNduwqjJmzhZe+mE9DauU5L07EqhUslBA7m2MMeFERJaqaoK3c7mpyitsiQh/61SLsXcksOXACa4ZMY8VO4+EOixjjMlRllAC6LL6FZn4YDsKFojg5ncXMnnF7lCHZIwxOcYSSoDVrVScyf3b06RaKQZ/toJXp6/nwoX8Wa1ojMlfLKEEQZmi0Xx8z8X0blWNkTM3c9/4JRw/cz7UYRljTFBZQgmS6KgIXriuEf/s2YDZGw9w3TsL2HrwZKjDMsaYoLGEEkQiwu1tYvnonos5dPIcPUfMY/bGrE35Yowx4c4SSg5oU6ssk/u3o0qpwtz1wSLGzNlsgyCNMXmOJZQcUq1MESY+2JZuDSvxwvfrGfr5Cpux2BiTp1hCyUFFoqMYeWtzHrm8Dt+s+J0bRy9kz9HToQ7LGGMCwhJKDhMRBlwax3t3JLD14El6DJ/P0u2+loExxpjcwRJKiHSNr8ikB9tSrGAkt4z5lc8W7Qh1SMYYky2WUEIorqIzCLJ1zbI8PjGRf3yTyLnkC6EOyxhjssQSSoiVLFKAD+5syd861uTjX3dw63u/sv+4rQRpjMl9LKGEgajICIZdVZ/hvZux5vdj9Bg+j2U7Doc6LGOMyRRLKGGkR5MqTHywLQWjIrn53YVMsHYVY0wuYgklzNSvXIIpA9rRplY5hk1M5IlJiZxNtvEqxpjwZwklDJUqEs0Hd7bkwc61+PS3HfQe8yv7jlm7ijEmvFlCCVOREcJj3erxzm3NWb/3ON2Hz7PxKsaYsGYJJcxd1agykx5sR5FoZ7zKJ79tt3nAjDFhyRJKLlC3UnGm9G9Pu9rl+L9Jq/n716tsHjBjTNjJMKGISA1/jpngKlmkAOP6tmTgpbX5Yskubhy9kJ2HToU6LGOM+S9/3lC+9nLsq0AHYjIWGSE8fHld3rsjgW0HT9JjxDzm2PoqxpgwkW5CEZF6InI9UFJEenlsdwKFcixC8xdd4ysyZWB7KhYvRN8PFjHil022br0xJuR8vaHUBboDpYAeHltz4L7gh2Z8qVGuKJP6t+WaJlV47ceN9PtoCUdP27r1xpjQkYx6DIlIG1Vd6NfNRNoAfYAOQGXgNLAa+A74WFWPZi/cwElISNAlS5aEOoxsU1U+XLCN579bR0zpwoy+vQX1KpUIdVjGmDxKRJaqaoLXc34klELAPUADPKq6VPXuNOWmAb8Dk4ElwH63fB3gEpy3mzdUdUqWv0kA5ZWEkmrJtkM8+Mkyjp05z8vXN6Zn06qhDskYkwf5Sij+NMp/BFQCrgBmAzHAcS/lblfVe1R1iqr+rqrJqnpCVZep6uuq2hlYkMXvYDKQEFuGqYPa07hqKQZ/toJnpqyxqfCNMTnKn4RSW1WfBE6q6ofA1UCjtIVU9WDaYyJymYj0EJEC6ZXxKNtNRDaISJKIPO7lfGkRmSQiq0RkkYg09Dg3VETWiMhqEZngvlXlOxWKF+KT+y7mnvY1+M+CbfR+71f2HrUpW4wxOcOfhJLa0nvE/U+8JBCb0UUi8jrQBWiNUw3mq2wkMBK4EogHeotIfJpiTwArVLUxcAfwlnttVWAQkKCqDYFI4BY/vleeVCAygie7xzO8dzPW7TnG1W/PZX5SunncGGMCxp+EMkZESgNPAlOAtcAraQuJyGsiUtLj0EU4SeAf7mdfWgFJqrpFVc8BnwE905SJB2YAqOp6IFZEKrrnooDCIhIFFMFpy8nXejSpwpQB7ShdNJrbx/1mXYuNMUGXYUJR1bGqelhVZ6tqTVWtoKqjvRSdBHwuIgPdN47xwK/ACmBMBo+pCuz02N/lHvO0EugFICKtgOpAjKruBl4DdgB7gKOq+qO3h4hIPxFZIiJLDhzI+wMCa1cozuT+7ejhdi2++8PFHD55LtRhGWPyKH+mXikoIreKyBMi8lTqlracqs5X1W7AEeAH99jFqtpEVd/O6DFejqX9dfoloLSIrAAGAsuBZPftqSdQA6gCFBWRPt4eoqpjVDVBVRPKly+fQUh5Q9GCUbx5c1P+eW1DFiT9Qffh81ix80iowzLG5EH+VHlNxvkPOxk46bH9iYhEicjVwD7gOqCZiEwRkcZ+PGMXUM1jP4Y01VaqekxV71LVpjhtKOWBrTjtNFtV9YCqngcmAm39eGa+ISLc3ro6X97fBoAbRy/go4XbbNZiY0xARflRJsZ988jINzjVW0WA21S1r4hUAZ4TEVVVX6PrFwNx7qSTu3Ea1W/1LCAipYBTbhvLvcAcVT0mIjuA1iJSBGcg5WU442BMGk2qlWLqwPY89MUKnpy8hiXbD/PCdY0oWtCffwbGGOObP/+TLBCRRqqamEG56qraXUSicdpOUNXfgXtFpKmvC1U1WUQGANNxemm9r6prROR+9/xooD4wXkRScDoG3OOe+01EvgKW4bxFLSfjNpt8q3TRaMb1bck7s5J446eNrPn9GKP7NKd2heKhDs0Yk8ulO1JeRBJx2jGigDhgC3AWp71D3e67nuUH4vwnr8DrqvpxEOPOtrw2Uj4r5icdZNCE5Zw+n8KLvRrZ6HpjTIayNPWKiFT3dVNV3R6A2ELGEopj79EzDPh0GUu2H6Z3q4t4ukc8hQpEhjosY0yYytLUK6q63U0az6d+9jzm5SH/cHtcpRfEpSLSPStfwARPpZKFmNCvNfd3qsWERTu47p0FbD34lz4XxhiTIX/aUBp47rhjTFp4KZcITBWRMzjtGQdwJoeMA5oCPwMvZCtaExQFIiN4/Mp6tKpRmoe+WEn3t+fy0vWN6dGkSqhDM8bkIr4W2BomIseBxiJyzN2O48wi/JepVFR1sqq2A+4H1uA0rh8DPgZaqepQVc37owlzsUvrVeS7QR2oW6k4Aycs5x/fJNra9cYYv/kzff2Lqjosh+LJMdaGkr7zKRd4dfoGxszZQoMqJXjntuZUL1s01GEZY8JAtqavz4vJxPhWIDKCJ66qz9g7Eth1+DTd357HtMQ9oQ7LGBPm/Bkpb/KpLvEV+W5Qe2pWKMYDnyzjmSlrOJtsVWDGGO8CnlBEpJ0/x0zuEFO6CF/+rc1/11i5cfRCtv9hvcCMMX/lq1G+SOrCWO5+XXchq14Z3HO4n8dMLhEd5ayx8u7tLdh28CTd357Htyvz/QoBxpg0fHUb/gFn5PsmEakNLAQ+AbqLSMu0bSsi0gZnUsbyIvKQx6kSOD2+TC53RYNKNKhSgkETljNwwnIWbD7IU90bUDja/nqNMb6rvEqr6ib3c19ggqoOxFlV0dsAxWigGE6SKu6xHQNuCFjEJqRiShfh87+14YHOtZiwaCc9R85j077joQ7LGBMGfL2hePYnvhR4FUBVz4nIhb8UVp0NzBaR/6ROyyIiEUAxVT0WwJhNiBWIjODv3erRumZZHvp8BT1GzOO5axpyY0IMIt6WtjHG5Ae+3lBWucv6DgVqAz/Cf6eR9+VFESkhIkVxZgXeICKPBiZcE0461SnPtMEdaFG9NI99vYohn6/g+JnzoQ7LGBMivhLKfcBBIBa4XFVPucfjcZbcTU+8+0ZyLfA9znryt2c/VBOOKpQoxPi7L+aRy+vw7crf6TF8Hom7joY6LGNMCPiaHPK0qr6kqoNVdaXHqZ1AJR/3LOD2DrsWmOyuomhLA+ZhkRHCgEvj+KxfG84mX6DXqPm8P2+rrQhpTD7j1zgUESknIg+IyBxgFlDRR/F3gW1AUWCOOw2+taHkA61qlOH7QR3oVKc8z01dy33jl3Do5LlQh2WMySG+1kMpjrM2/K1AHWAScLOqxmT6ISJRqpqcnUADzebyCh5V5YP523hp2npKFy3Av29qStva5UIdljEmALI6l9d+nHEo/wJqqerDQIa/bopISRF5Q0SWuNvrOG8rJp8QEe5uX4OJD7alaMEobhv3G69OX8/5lL90DjTG5CG+EsoTOOuZjAKGiUgtP+/5PnAcuMndjgEfZCdIkzs1rFqSqQPbc1OLaoycuZkbRy9k56FTGV9ojMmVfDXK/1tVLwauwVlH/hugioj8XUTq+LhnLVV9WlW3uNuzQM3Ahm1yiyLRUbx8Q2NG3NqMzQdOcNVbc5m8YneowzLGBIE/09dvUdV/qWojoCVQEpjm45LTItI+dcedGPJ0tiM1uVr3xlX4flAH4ioWY/BnK3jky5WcPBtWzWrGmGzy1SjfWlV/zfQNRZoCH+IkHgEOAX1VdVV2Ag00a5QPjeSUC7w9YxMjZiZRvWxR3r6lGY1iSoY6LGOMn7LaKP+Oxw0W+vswVV2hqk2AxkAjVW0WbsnEhE5UZAQPXV6XT+9rzZnzKfQaNZ/35mzhwgUbs2JMbucroXhOylTI3xuKSFkReRtnvMpMEXlLRMpmMT6TR7WuWZZpgztwab0K/Ov7dfT9YBH7jp0JdVjGmGzwlVAiRKS0mwxSP5dJ3Xxc9xlwALgeZ5bhA8DngQvZ5BWlikQzuk8LXriuEYu3HaLbm3OYvmZvqMMyxmSRrzaUbcAF/vymkkpV1WvPLbd+rUWaY0vSq3MLFWtDCS9J+08w5PPlrN59jN6tLuLJ7vUpEu1rMmxjTChkqQ1FVWNVtaaq1vCy+eoGPFNEbhGRCHe7CfjOjyC7icgGEUkSkce9nC8tIpNEZJWILBKRhh7nSonIVyKyXkTWuYt9mVykdoViTHygHX/rVJPPFu+g+9s2yaQxuY2vN5QKOIMbawOrgJf8WddERI7jjIxPHRYdAaQuQq6qWsLLNZHARqArsAtYDPRW1bUeZV4FTqjqsyJSDxipqpe55z4E5qrqWBGJBoqo6hFfcdobSvhakHSQh75YyR8nz/Lw5XXp16EmERG2zoox4SCrvbzG4ySC4TgrL77tz8NUtbiqRqhqlLtFuMeKe0smrlZAkjvm5RxOO0zPNGXigRnuM9YDsSJSUURKAB2Bce65cxklExPe2tYuxw9DOtClfkVemrae28b+xp6jNpTJmHDnK6FUUtX/U9Xp7tK/jf29qYg0FpFrRKRX6pbBJVVxpsVPtcs95mkl0Mu9fyugOhCDMwr/APCBiCwXkbHu4l7e4uqXOsfYgQMH/P06JgRKFYnmndua88oNjVm56wjd3pzL94l7Qh2WMcYHn92GPXt2AZFp9tO76H2c+byuB3q4m7c16P90mZdjaeviXgJKi8gKYCCwHEjGWca4OTBKVZvhvFX9pQ0GQFXHqGqCqiaUL18+g5BMqIkINyVU47tBHYgtW4QHP1nGI1+utFUhjQlTvrrRlASW8uf/7Je5fyrpz8/VWlXjMxnHLqCax34M8LtnAbf95i5wMh2w1d2KALtU9Te36Fekk1BM7lSjXFG+eqAtb/28iXdmJfHrlj/4981NaRnrq/e6MSanBaOX10IRyWxCWQzEiUgNt1H9FmCKZwG3J1e0u3svMEdVj6nqXmCniNR1z12Gs5a9yUMKREbwyBV1+fL+NkSIcPO7C3nlh/WcS7Yp8Y0JF8Ho6P8hTlLZC5zFecNRVU23DUZVk0VkADAdiATeV9U1InK/e340UB8YLyIpOAnjHo9bDAQ+cRPOFtw3GZP3tKhehu8Hd+Cf367lnVmbmb3xAG/e3JS4isVDHZox+V663YazfEORJOAhIJH/dR1GVbcH9EHZZN2Gc7/pa/YybGIiJ88m8/iV9ejbJta6FxsTZL66DQfjDWWHqk7JuJgx2XNFg0o0u6gUj3+dyLPfruWX9ft59YYmVCrp99RzxpgAyvANJZ0eXcdV1WtXGxF5BygFfItT5QWAqk7MRpwBZ28oeYeq8umiHTw/dR3RURH867qGdG9cJdRhGZMnZXVgY6plOOM8NgKb3M9bRWSZiLTwUr4wTiK5HP+7DRuTZSLCbRdX57tB7YktV5QBny5n6OcrOHrauhcbk5P8eUMZDUxS1enu/uVAN+AL4C13meBcx95Q8qbzKRcY8UsSI2YmUaF4QV65oTEd4mzMkTGBkt03lITUZAKgqj8CHd3VHAt6eViMO4njfhHZJyJfi0hMlqM3JhMKREYwtGsdJj7QlsLRkdw+bhFPTV7NqXO23LAxweZPQjkkIn8Xkeru9hhw2J3Q0dsggA9wxpBUwZk+5Vv3mDE5pkm1Unw/qAN3tYtl/MLtXP32PJbtOBzqsIzJ0/xJKLfijFz/BpgMXOQeiwRu8lK+vKp+oKrJ7vYfwOocTI4rVCCSp3s04NP7LuZc8gVuGLWA16ZvsMGQxgRJhglFVQ+q6kB3bfimqjpAVQ+4s/omebnkoIj0EZFId+sD/BH40I3xT9ta5Zg2pAPXN49hxMwkrh05n/V7M1yJwRiTSRkmFBGpIyJjRORHEfkldfNxyd04by57gT04ywDfHZhwjcmaEoUK8OqNTRhzewv2Hz/DNcPnM3r2ZlIuBHZgrzH5mT+9vFYCo3EmikxJPa6qS4MbWnBZL6/8648TZ3liUiLT1+wjoXppXr+pCdXLel3xwBiTRnZ7eSWr6ihVXaSqS1M3Lw95JXXurTTHh4rIy1mI25igKFusIKP7tOCNm5qwYd9xur05l/ELt3HB3laMyRZ/Esq3IvKgiFROXQslndHz3YExXo6/BVydrSiNCTARoVfzGKYP6UjLGmV4avIa+oz7jZ2HToU6NGNyLX8SSl/gUWABTrXXUsBbXZGq6l+6z7jHbMY+E5aqlCrMh3e15MVejVi58wjd3pzDp7/tINCTphqTH/jTy8vf9VBOiUhc2oPuMVsQ3IQtEaF3q4uYPrQjTaqV4olJidzx/iJ+P2L/bI3JjHQb5UXkUlX9Jb314NNO9igiVwLDgedx3mIAEoBhwBBV/T5gUQeANcobby5cUD75bTsvfL+eqAjhqR7x3NAiBmeRUGNMVqev7wT8gjO5Y1oK/CmhqOo0EbkWp3psoHt4NXC9qiZmOmpjQiAiQri9TSyd6lTgka9W8uhXq/hh9V5e6NWIiiVsWnxjfAn4Alu5hb2hmIxcuKD8Z8E2Xpm+noJRkTx7TQN6Nq1ibysmX8vWAlsiUgq4A4j1LK+qgwIVoDHhKCJCuLt9DTrXLc+jX61iyOcrmLpqDy9c15AK9rZizF/408vre5xkksj/ennl6kGNxmRGzfLF+OJvbfjH1fWZu+kAXd6YzddLd1lPMGPS8Gek/DJVbZ5D8eQYq/IyWbHlwAn+/vUqFm87zCV1y/Nir8a25LDJV7I7Uv4jEbnPj4GNqQ+rIyIzRGS1u99YRP6RxdiNCSs1yxfj835teLpHPAu3/EHXf8/mi8U77W3FGPxLKOeAV4GF+B7YmOo9nK7C5wFUdRVwS/bCNCZ8REQId7WrwfQhHYmvXILHvl5F3w8Ws9vGrZh8zp+E8hBQW1VjMxjYmKqIqi5Kc8yWyzN5TvWyRZlwX2v+2bMBS7Yd4op/z2HCIhtlb/IvfxLKGiAzExwdFJFaOGNVEJEbcKaxNybPSR23Mn1IRxrHlGTYxERuH7fI5gQz+ZI/jfKTgAbATOBs6vH0ug2LSE2cSSLbAoeBrUAfVd0WmJADwxrlTaCpKp8u2sEL361DgceuqMsdbWKJiLBxKybv8NUo709C6evtuKp+mMF1RYEIVT3ub6A5yRKKCZbdR07zxMREZm88QIvqpXn5+kbUrlA81GEZExDZSijuDaKBOu7uBlU976NsQeB6/joQ8rlMxBx0llBMMKkqk5bv5rmpazl1NoXBXeLo17EmBSL9qWU2Jnxlq9uwiHQGNgEjgXeAjSLS0cclk4GeOA3xJz22jJ7TTUQ2iEiSiDzu5XxpEZkkIqtEZJGINExzPlJElovI1IyeZUywpa638tPQTnRtUJFXp2+g54j5rN59NNShGRM0/lR5LQVuVdUN7n4dYIKqtkin/GpVbejtnI9nRAIbga7ALmAx0FtV13qUeRU4oarPikg9YKSqXuZx/iGc2Y1LqGr3jJ5pbygmJ01fs5d/fLOaQyfP0a9jTQZfFkehApGhDsuYTMvuwMYCqckEQFU3AgV8lF8gIo0yGWMrIElVt6jqOeAznLccT/HADDeG9UCsiFQEEJEYnFUhx2byucbkiCsaVOLnoZ24vnlVRs3azFVvzWXxtkOhDsuYgPInoSwRkXEi0tnd3sPLXF4ikigiq4D2wDK3+mqVx3FfqgI7PfZ3ucc8rQR6uc9qBVQHYtxzbwKPAX9ZMTJNjP1EZImILDlw4EAGIRkTWCWLFOCVG5rw8T0Xcy7lAjeOXsiT36zm+Jl0mySNyVUynG0YeADoDwzCWcp3Dk5bSloZVjP54K1fZdq6uJeAt0RkBc5ElcuBZBHpDuxX1aVue0+6VHUM7rr3CQkJNvrMhET7uHJMH9KR137cwH8WbOPndfv4Z8+GdImvGOrQjMkWf9pQigJnVDXF3Y8ECqqq15FbIvKRqt6e0bE059sAz6jqFe7+MABVfTGd8oIzvqUxzjQvt+N0AigElAAmqmofX9/L2lBMOFi24zDDvk5kw77jXN24Mk/3iKdCcZts0oSv7LahzAAKe+wXBn72Ub5BmodHAl4b8D0sBuJEpIbbRfkWYEqa+5RyzwHcC8xR1WOqOkxVY1Q11r3ul4ySiTHhovlFpfl2YHse7lqHn9bso8vrNtmkyb38SSiFVPVE6o77uUjaQiIyTESOA41F5Ji7HQf243QlTpeqJgMDgOnAOuALVV0jIveLyP1usfrAGhFZD1wJDPYjdmPCXnRUBAMvi2PakA7UcyebvPW939h2MMPe9saEFX+qvOYDA1V1mbvfAhihqm3SKf+iqg4LeKQBZlVeJhxduKB8tngnL36/jnMpFxjcJY77OtiASBM+sjv1Skucbry/u4cqAzeraq5etdESigln+46d4enJa/hhzV7qVy7By9c3onFMqVCHZUxApl4pANTF6Y213tfUK7mFJRSTG/ywei9PTV7NwRNnuatdDR7qWoeiBf3pnGlMcPhKKP7+y2zJ/+bmaiYiqOr4AMVnjElHt4aVaFOrLC//sJ5x87YyLXEPz1kXYxOm/Kny+gioBawAUtzDmt709e41kUBF/jw55I5sRxtA9oZicpul2w8xbGIiG/ed4MqGlXjmmgZULGFdjE3Oym4byjogXv3sxygiA4GngX38b+S6qmpj/0MOPksoJjc6l3yB9+Zu4e0Zm4iOjOCxbnW59eLqRNqaKyaHZHccymqgUiaeNxioq6oNVLWRu4VVMjEmt4qOiqD/JbWZPqQjTaqV4snJa7h+1ALW7TkW6tCM8SuhlAPWish0EZmSuvkovxOwObqNCaLYckX56J5WvHlzU3YeOkWP4fN4adp6Tp9LyfhiY4LEn0b5ZzJ5zy3ALBH5jj8vGfxGJu9jjPFBRLi2WVU61SnPi9PWMXr2Zr5L/J3nr21EpzrlQx2eyYcyTCiqOjuT99zhbtHuZowJotJFo3nlhib0ah7DE5MS6fv+Iro3rsyT3eOt0d7kqHQb5d1pU7ydFJxG9hLBDCzYrFHe5EVnk1N4d/YWRsxMomBkBI9cUZc+ra3R3gROtgc2+vmQN1V1iIh8i5dEpKrXBORBAWIJxeRl2w6e5MnJq5m76SCNqpbkhesa0SimZKjDMnlATiWUFu6aJJ28nc9C1VlQWUIxeZ2qMnXVHp6bupY/TpzljjaxPHR5HUoU8rXgqjG+5UhCyW0soZj84tiZ87w+fQPjf91O+WIFeapHPFc3qoyzrJAxmZOlcSgiUjB4IRljckqJQgV4tmdDJvdvR8UShRjw6XL6frCY7X/Y9PgmsHyNQ1kI/516xRiTyzWOKcU3/dvx7DUNWLb9MF3/PYe3Z2zizHkbu2ICw1e34WgR6Qu0FZFeaU+q6sSMbi4iEUAxVbVhvMaEgcgIoW/bWLo1rMQ/p67ljZ82MnHZLp7t2dDGrphs8/WGcj/QGigF9EizdU/vIhH5VERKuGvRrwU2iMijgQvZGJNdFUsUYsStzfn4nouJEKHv+4t44OOl/H7kdKhDM7mYP5ND3qOq4/y+ocgKVW0qIrfhrCX/d2BpuM3nZY3yxjjOJqcwdu5Whv+yiQgRBl8Wx93ta9gqkcar7E4O+ZGIDBKRr9xtoLvgVnoKuOevBSa7i3Hlz65kxuQCBaMi6X9JbX4a2om2tcrx4rT1XPXWXAnTNe0AABghSURBVH7d8keoQzO5jD8J5R2cN4133K05MMpH+XeBbUBRYI6IVAesDcWYMFetTBHG9k1g7B0JnD6fwi1jfmXo5yvYf/xMqEMzuYQ/VV4rVbVJRscyuEeUqiZnMcagsCovY9J3+lwK78xK4t3ZWygYFcHDl9ehT+vqRFk1WL6X3SWAU0Sklqpudm9Wk/+t3OjtYQWB6/nfksGpnvM7YmNMSBWOjuThy+tyXbOqPD1lDc98u5bPl+ziuZ4NaBlbJtThmTDlz68bjwIzRWSWiMwGfgEe9lF+MtATSAZOemzGmFymZvlijL+7FaNua87RU+e4cfRCpxrsmFWDmb/ya+oV962jLs5Mw+tV9ayPsqtVtWHgQgwOq/IyJnNOnUtm5Mwk3puzleioCIZ0iaNv21jrDZbPZLeXF6p6VlVXqepKX8nEtUBEGmU6SmNMWCsSHcWjV9Rj+tCOJMSW5vnv1nH123NZuNl6gxlHMH61aA8sFZENIrJKRBJFZFUQnmOMCYEa5YrywZ0tGXN7C06dS6H3e78ycMJy9hy1QZH5nT+N8pl1ZVYuEpFuwFtAJDBWVV9Kc7408D5QCzgD3K2qq0WkGjAeqARcAMao6lvZiN8YkwER4fIGlehYpzyjZm1m1OzNzFi3j4GXxnFP+xpER1k1WH6U4d+6iMzw81jqCo7H09l8PSMSGImTjOKB3iISn6bYE8AKd8T9HTjJB5zG/4dVtT7OVDH9vVxrjAmCQgUiGdq1Dj+7gyJf/mE93d6cw6wN+0MdmgkBX9PXFxKRMkA5ESktImXcLRao4uWST90/lwJL3D+Xeuz70gpIUtUtqnoO+Aynp5ineGAGgKquB2JFpKKq7lHVZe7x48A6oGoGzzPGBNBFZZ1BkR/c2RIF7vxgMfd+uJhtB62DZ37iq8rrb8AQnOSxFKeHFzij3kemLayq3d0/a2QhjqrATo/9XcDFacqsBHoB80SkFVAdiAH2pRZwk10z4LcsxGCMyaZL6lWgbe2yfDB/G8NnbOLyf8/h3g416H9JbYoWDEYNuwkn6f4Nu+0Qb4nIQFUdnpmbuu0dcUAhj/vN8XWJtxDS7L/kxrMCSASW41R3pT6zGPA1MCS96fJFpB/QD+Ciiy7K+IsYYzKtYFQk93eqxXXNqvLytPW8M2szE5ftZthV9bimSRVbKTIP83ccSlvSjHxX1fHplL0XGIzz9rACp11joape6uP+bYBnVPUKd3+Y+4wX0ykvwFagsaoecyejnApMV9U3MvxC2DgUY3LK0u2HeHrKGlbvPkbL2NI8c00DGlQpGeqwTBZlaxyKu2LjazjdgVu6m9ebuQa7Zbar6iU4VVAHMnjMYiBORGqISDRwCzAlTRyl3HMA9wJz3GQiwDhgnb/JxBiTc1pUL8Pk/u15qVcjNh84SY/h8/i/SYkcPnku1KGZAPOnUjMBiFd/XmUcZ1T1jIggIgVVdb2I1PV1gaomi8gAYDpOt+H3VXWNiNzvnh8N1AfGi0gKzsJd97iXtwNuBxLd6jCAJ1T1ez/jNcYEWWSEcEuri7iyYWXenLGR8Qu3M3XVHoZ2ieO21tVttH0e4c9sw18Cg1R1j183FJkE3IXToH8pcBgooKpXZTPWgLIqL2NCZ8Pe4zw3dQ3zk/4grkIxnuoRT4c4W4I4N/BV5eVPQpkJNAUWAf+ddkVVr/HjwZ2AksAPbnfgsGEJxZjQUlV+WruP579bx45Dp+hSvyL/uLo+seWKhjo040N2E0onb8dVdbaPa0oD1fhzI/4yv6LNIZZQjAkPZ5NT/tvN+FzKBe5uV4MBl9ameCFfC8OaUMlWQnFvUB2IU9WfRaQIEOkOIvRW9p/AncAWnKlQANRXL69QsIRiTHjZf/wMr/6wgS+X7qJcsWgevaIuN7SoRmSEdTMOJ9l9Q7kPZ+xGGVWtJSJxwGhVvSyd8huARuFWxZWWJRRjwtOqXUd49tu1LN1+mIZVS/B0D1vUK5xkd/r6/jg9qY4BqOomoIKP8quBUpkN0hhjABrHlOKr+9vw1i1N+eOEs6jXgE+XsevwqVCHZjLgT7fhs6p6LnV0q4hE8ddR7J5eBJaLyGoy2YhvjDHgzGbcs2lVusZX5N3ZW3h3zmZ+XLuPe9vX4MFLalPMpnEJS/78rcwWkSeAwiLSFXgQ+NZH+Q+Bl3GmR7ngo5wxxvhUJDqKoV3rcHPLarw6fQPvzNrMF0t28egVdax9JQz504YSgTOI8HKcObem46xX4vVCEZmtql57hoUTa0MxJvdZsfMIz327hmU7jhBfuQT/6F6ftrXKhTqsfCW7jfJFcUa/p7j7kUBBVfVaoSkib+BUdU3hz1Ve1m3YGJNtqsrUVXt4adp6dh85zeXxFXniKhu/klOym1B+Bbqo6gl3vxjwo6q2Taf8TC+HrduwMSagzpxPYdy8rbwzM4lzKRe4s20sAy6No2RhG78STL4Sij9tKIVSkwmAqp5wx6J45U4IaYwxQVWoQCT9L6nNjS1ieO3HDYydt5Wvlu5iSJc63HrxRTY/WAj48xM/KSLNU3dEpAVwOr3CIlJRRMaJyDR3P15E7kmvvDHGZEeFEoV45YYmfDugPfUqleDpKWu44s05/LR2H/7PaWsCwZ+EMhj4UkTmishc4HNggI/y/8FpuE9dJngjzkSRxhgTNA2rluTT+y5m7B1Obcx945fQ+71fWb37aIgjyz98Vnm5PbyigXpAXZxeXutV9byPy8qp6hcei2Qlu1POG2NMUIkIXeIr0qlueSYs2sGbP2+ix4h5XNesKo9eUZfKJQuHOsQ8zecbiqpeAF5X1fOqulpVEzNIJuBUkZXFHfwoIq0B+xXBGJNjCkRGcEebWGY92pl+HWsyddUeLnltFq//uIETZ5MzvoHJEn+qvH4UkevF/4WgH8LpMlxLROYD44GBWQ3QGGOyqkShAgy7sj4zHupE1/hKDP8lic6vzmLCoh0kp9i460Dzp9vwcaAokILTGC843YBLeCkbgbOG/CL+V0W2wY+3mhxn3YaNyX+W7zjMv75bx5Lth6lTsRjDrqxP57rl8f/3ZZPt6esz+bCFqtomoDcNAksoxuRPqsoPq/fy8g/r2fbHKdrULMsTV9WnUUzJUIeWK2RrtmFx9BGRJ939aiLSysclma0iM8aYHCMiXNmoMj8O7cQzPeLZsO84PUbMY/Bny9l5yGY0zg5/qrxG4UzyeKmq1ndXY/xRVVumUz61iiwZOIOPKrJQsjcUYwzAsTPnGT1rM+PmbUUV7mwXS//OtSlZxEbce5Pd9VAuVtX+OMkBVT2M05XYK1UtrqoRqhqtqiXc/bBKJsYYk6pEoQI81q0esx7tzDVNq/De3C10fHUmY+du4WyyjXjIDH8Synl3QsjUbsDlsWnpjTF5TOWShXntxiZ8N7ADTaqV4vnv1tHljdlMWfk7Fy7YiHt/+JNQ3gYmARVE5F/APOCFoEZljDEhEl+lBOPvbsVH97SiWMECDJqwnJ4j5zM/6WCoQwt7fvXyEpF6wGU47SEzVHVdsAMLNmtDMcZkJOWCMnnFbl7/cSO7j5ymY53y/L1bXRpUyb89wrLUbVhECgH3A7VxVl8cp6p5ZoipJRRjjL/OnE/h41+3M2JmEkdPn+faplV5qGsdqpVJd+L1PCurCeVz4DwwF7gS2KaqeWaSR0soxpjMOnr6PKNnb+Z9t0fY7W2qM+CS2pQumm4/pTwnqwklUVUbuZ+jgEWq2txr4VzIEooxJqv2HD3Nmz9t4sulOykaHcX9nWtxd7saFI6ODHVoQZfVbsP/nS4lJ6q6RKSbiGwQkSQRedzL+dIiMklEVonIIhFp6O+1xhgTSJVLFublGxozfUhHLq5Zllenb6DzazPz/Rxhvt5QUoCTqbtAYeAUQRio6HZL3gh0BXYBi4HeqrrWo8yrwAlVfdbtJDBSVS/z51pv7A3FGBMoi7cd4sXv17FsxxFqlivKI1fU5cqGlfLkHGFZekNR1Uh3YGLq4MSoIA5UbAUkqeoWVT0HfAb0TFMmHpjhxrYeiBWRin5ea4wxQdMytgxfP9CW9+5IICpSePCTZfQcOZ95m/JXV+NwWXS5KrDTY3+Xe8zTSqAXgDuXWHUgxs9rca/rJyJLRGTJgQMHAhS6McY4c4R1ja/ItMEdef3GJvxx4hx9xv1Gn7G/sWrXkVCHlyPCJaF4ey9MWxf3ElBaRFbgrK+yHGe+MH+udQ6qjlHVBFVNKF++fHbiNcYYryIjhOtbxPDLI514qns8a/cc45oR83nwk6VsPnAi1OEFlc8lgHPQLqCax34M8LtnAVU9BtwFzgzIwFZ3K5LRtcYYk9MKRkVyd/sa3NSyGmPnbuG9OVuYvmYfNyXEMOiyuDy5HHG4vKEsBuJEpIaIRAO34Kz6+F8iUso9B3AvMMdNMhlea4wxoVKsYBRDutRh9mOXcEeb6ny1dBedX53FC9+v49DJc6EOL6DCIqG43ZIHANOBdcAXqrpGRO4XkfvdYvWBNSKyHmeg5WBf1+b0dzDGGF/KFSvI0z0a8MvDneneuApj526h4yszefPnjRw/E3aL2mZJwFdszC2s27AxJpQ27TvOGz9tZNrqvZQuUoD+l9SmT+vqFCoQ3oMjc3QJ4NzCEooxJhys2nWE137cyJyNB6hUohCDLovjxoQYCkSGRQXSX2R3gS1jjDFB0jimFOPvbsVn/VpTtXRhnpiUSJc3ZjN5xe5ctw6LJRRjjAkDrWuW5av72/D+nQkUiY5i8GcruOrtufy4Zi+5pSbJEooxxoQJEeHSehX5bmB7hvduxtnkC/T7aCnXjpzPnI0Hwj6xWEIxxpgwExEh9GhShZ+GduSVGxpz8MQ57nh/ETe/+yu/bfkj1OGlyxrljTEmzJ1NTuGLxTsZ/ksS+4+fpUNcOR7qWodmF5XO8Visl5cXllCMMbnNmfMpfLRwO6Nmb+bQyXN0qV+Bh7rWJb5KoOfrTZ8lFC8soRhjcqsTZ5P5z/ytvDtnC8fPJHN148oM7RJH7QrFg/5sSyheWEIxxuR2R0+dZ+y8Lbw/byunz6fQs2lVBl0WR41yRYP2TEsoXlhCMcbkFX+cOMu7c7YwfuE2zqcovZpVZeClcVxUtkjAn2UJxQtLKMaYvGb/8TOMnrWFj3/bzoULyo0JMfS/pDYxpQOXWCyheGEJxRiTV+09eoZRs5KYsGgninJzy2r0v6R2QKbMt4TihSUUY0xet/vIaUbOTOKLxTuJiBBubXURD3auRYUShbJ8T0soXlhCMcbkFzsPnWLEL0l8tWwXURHCsCvrcWe7Glm6l00OaYwx+Vi1MkV4+YbGzHioE90bV6FamcA31kP4LAFsjDEmyGLLFeX1m5oE7f72hmKMMSYgLKEYY4wJCEsoxhhjAsISijHGmICwhGKMMSYgLKEYY4wJCEsoxhhjAsISijHGmIDIt1OviMgBYHsWLy8HHAxgOLmBfef8wb5z/pCd71xdVct7O5FvE0p2iMiS9OayyavsO+cP9p3zh2B9Z6vyMsYYExCWUIwxxgSEJZSsGRPqAELAvnP+YN85fwjKd7Y2FGOMMQFhbyjGGGMCwhKKMcaYgLCEkkki0k1ENohIkog8Hup4gk1EqonITBFZJyJrRGRwqGPKCSISKSLLRWRqqGPJKSJSSkS+EpH17t93m1DHFEwiMtT9N71aRCaISNYXWg9jIvK+iOwXkdUex8qIyE8issn9s3QgnmUJJRNEJBIYCVwJxAO9RSQ+tFEFXTLwsKrWB1oD/fPBdwYYDKwLdRA57C3gB1WtBzQhD39/EakKDAISVLUhEAncEtqoguY/QLc0xx4HZqhqHDDD3c82SyiZ0wpIUtUtqnoO+AzoGeKYgkpV96jqMvfzcZz/ZKqGNqrgEpEY4GpgbKhjySkiUgLoCIwDUNVzqnoktFEFXRRQWESigCLA7yGOJyhUdQ5wKM3hnsCH7ucPgWsD8SxLKJlTFdjpsb+LPP6fqycRiQWaAb+FNpKgexN4DLgQ6kByUE3gAPCBW9U3VkSKhjqoYFHV3cBrwA5gD3BUVX8MbVQ5qqKq7gHnl0agQiBuagklc8TLsXzR71pEigFfA0NU9Vio4wkWEekO7FfVpaGOJYdFAc2BUaraDDhJgKpBwpHbZtATqAFUAYqKSJ/QRpX7WULJnF1ANY/9GPLoa7InESmAk0w+UdWJoY4nyNoB14jINpwqzUtF5OPQhpQjdgG7VDX17fMrnASTV3UBtqrqAVU9D0wE2oY4ppy0T0QqA7h/7g/ETS2hZM5iIE5EaohINE4j3pQQxxRUIiI49errVPWNUMcTbKo6TFVjVDUW5+/3F1XN87+5qupeYKeI1HUPXQasDWFIwbYDaC0iRdx/45eRhzsheDEF6Ot+7gtMDsRNowJxk/xCVZNFZAAwHadXyPuquibEYQVbO+B2IFFEVrjHnlDV70MYkwmOgcAn7i9LW4C7QhxP0KjqbyLyFbAMpyfjcvLoFCwiMgHoDJQTkV3A08BLwBcicg9Ocr0xIM+yqVeMMcYEglV5GWOMCQhLKMYYYwLCEooxxpiAsIRijDEmICyhGGOMCQhLKCZkROQ6EVERqReCZ28TkXI5+LzvRaRUBmXGpk68KSJP+HnfJ9LsL8h6lJknIp1FJNMDAkWkmYhkaq40EflMROIy+yyTc6zbsAkZEfkCqIwz6+kzOfzsbTgzzR7MwrVRqpoc+Kj+9IwTqlosUOWyGUu631dEngFOqOprmbznl8DzqroyE9d0Avqo6n2ZeZbJOfaGYkLCnRusHXAPHtOGu7/xzvJYl+MTdyRz6lvFsyKyTEQSU99sROQZEXnE4x6r3YksEZFvRGSpu+5FPz/iOiEir7vPmCEi5d3js0TkBRGZDQwWkfIi8rWILHa3dqnfS0Q+cONbJSLXe8ReTkRi3e/1oXv+KxEp4vGMBBF5CWcW3BUi8kl63yOdcifcP0VEXnV/FokicnNGP980P4e037eHiPwmzsSRP4tIRfdnfD8w1I2hQ3o/lzT3Lg40Tk0m7t/fhyLyo/tz6iUir7hx/yDO1D8Ac4Eu4swObMKRqtpmW45vQB9gnPt5AdDc/dwZOIozT1oEsBBo757bBgx0Pz8IjHU/PwM84nHv1UCs+7mM+2dh93hZj3uV8xKXAre5n58CRrifZwHveJT71COui3CmpgF4GXjTo1xpz+cBse4z2rnH30+N3X1Ggvv5RJq40vseacudcP+8HvgJZ0aHijijoSv7+vmmuU/a71ua/9Vo3Au8ns7P3uvPJc29LwG+9th/BpgHFMBZh+UUcKV7bhJwrUfZn4AWof73a5v3zTK9CZXeONPEgzMJY2+caTAAFqnqLgBxpnuJxfkPB5xJ/ACWAr38eM4gEbnO/VwNiAP+8FH+AvC5+/ljj+fhcRycyQXjPX65L+H+5t0FjzcuVT3s5Rk7VXW+xzMG4UylHsjv0R6YoKopOBMBzgZaAsfw/fP15Pl9Y4DPxZlIMBrYms5zvf5c1FlLJ1VlnKnyPU1T1fMikoiTBH9wjye68aXajzM7cH6bDTpXsIRicpyIlAUuBRqKiOL8B6Ii8phb5KxH8RT+/O/0rJfjyfy5+raQ+5zOOP/BtVHVUyIyK/VcJng2Mp70+Bzh3ve0Z2G3+iijhsm0532Wz+L38LbUQipfP19Pnt93OPCGqk5x43kmnWu8/lzSOM1f4z8LoKoXROS8uq8jOAneM75C7vUmDFkbigmFG4DxqlpdVWNVtRrOb7zts3i/bbhTrYtIc5w1LgBKAofd/4Tr4SxhnJEINz6AW/H+mzvAj8CA1B0RaZrOcW9rdV8k/1uvvXc6zzjv0Xbg63t4lvM0B7hZRCLddqCOwKJ0vos/SgK73c99PY4fB4p77Kf3c/G0DqidxTjqAHl9QtZcyxKKCYXeOHXjnr7G+Q88K74GyrjVNw8AG93jPwBRIrIK+Cfwqx/3Ogk0EJGlOG9Rz6VTbhCQ4Dasr8VpnAZ4HijtNoavxGkvSGsd0NeNqwwwykuZMcAqt7Hd1/fwLOdpErAKWAn8AjymzhT1WfUM8KWIzAU8e8Z9C1yX2ihP+j+X/1LV9UBJt4rQbyJSETit7kqDJvxYt2FjPEiQu+G6PaOmqmrDYD0jNxCRocBxVfV7LIp7zTFVHRe8yEx22BuKMSYURvHnthx/HAE+DEIsJkDsDcUYY0xA2BuKMcaYgLCEYowxJiAsoRhjjAkISyjGGGMCwhKKMcaYgPh/f5Na6QZw9/gAAAAASUVORK5CYII=\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.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
