{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Thermal Treatment"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import PFAS_SAT as ps\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from IPython.display import Image\n",
    "import pandas as pd\n",
    "pd.set_option('display.max_colwidth', 0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Model document\n",
    "\n",
    "Combustion systems include mass burn incinerators for MSW, hazardous solid waste, and medical waste, as well as cement kilns. The combustion process model in the SAT framework was designed so that by changing default parameters, a variety of conventional and emerging thermal treatment processes can be represented. In addition to MSW incinerators, other potential thermal processes include gasification/pyrolysis systems, and thermomechanochemical systems, which may include mechanical and chemical processing prior to combustion. \n",
    "To begin, the incoming materials are partially or fully combusted. During combustion, a fraction of the PFAS may be volatilized and/or destroyed, but data on destruction and release of volatilized PFAS and/or products of incomplete combustion are limited. The volatilized compounds may be transformed into chlorofluorocarbons, which deplete stratospheric ozone and are potent greenhouse gases. The incineration process produces combustion residuals that may include bottom ash, fly ash, and air pollution control (APC) residues. The mass of PFAS flows are modeled using a destruction and removal efficiency (DRE), and fraction of the DRE that remains in the residual. Based on these inputs, the model calculates the PFAS destroyed, released to the air and remaining in the ash. The ash is then sent to landfill.\n",
    "\n",
    "<img src=\"../Images/ProcessModels/ThemalTreatment.png\" alt=\"Drawing\" style=\"width: 500px;\"/>\n",
    "\n",
    "\n",
    "### Assumptions and Limitations\n",
    "1.\tThe model assumes a single DRE and fraction remaining on the residual for each material and type of PFAS. Data is relatively limited, but this can be readily updated as more data becomes available.\n",
    "2.\tThe model assumes there is a single ash stream with a common set of parameters.\n",
    "3.\tThe model assumes a single VS destruction efficiency, but this could vary by material.\n",
    "4.\tThe model does not consider transformations of PFAS."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for Thermal Treatment 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>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>Total destruction and removal efficiency</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>DRE of PFOA</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>DRE of PFOS</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>DRE of PFBA</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>DRE of PFPeA</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>DRE of PFHxA</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>DRE of PFHpA</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>DRE of PFNA</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>DRE of PFDA</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>DRE of PFBS</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>DRE of PFHxS</td>\n",
       "      <td>0.97</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>Fraction of PFOA remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>Fraction of PFOS remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>Fraction of PFBA remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>Fraction of PFPeA remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>Fraction of PFHxA remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>Fraction of PFHpA remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>Fraction of PFNA remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>Fraction of PFDA remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>Fraction of PFBS remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>Fraction of not destroyed and removed that remains in residuals</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>Fraction of PFHxS remained in residuals</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>Combustion parameters</td>\n",
       "      <td>frac_vs_to_res</td>\n",
       "      <td>Fraction of incoming Volatile Solids to Combustion Residuals</td>\n",
       "      <td>0.05</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>Combustion parameters</td>\n",
       "      <td>res_ts_cont</td>\n",
       "      <td>Solids content of residuals</td>\n",
       "      <td>0.85</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                           Category  \\\n",
       "0   Total destruction and removal efficiency                          \n",
       "1   Total destruction and removal efficiency                          \n",
       "2   Total destruction and removal efficiency                          \n",
       "3   Total destruction and removal efficiency                          \n",
       "4   Total destruction and removal efficiency                          \n",
       "5   Total destruction and removal efficiency                          \n",
       "6   Total destruction and removal efficiency                          \n",
       "7   Total destruction and removal efficiency                          \n",
       "8   Total destruction and removal efficiency                          \n",
       "9   Total destruction and removal efficiency                          \n",
       "10  Fraction of not destroyed and removed that remains in residuals   \n",
       "11  Fraction of not destroyed and removed that remains in residuals   \n",
       "12  Fraction of not destroyed and removed that remains in residuals   \n",
       "13  Fraction of not destroyed and removed that remains in residuals   \n",
       "14  Fraction of not destroyed and removed that remains in residuals   \n",
       "15  Fraction of not destroyed and removed that remains in residuals   \n",
       "16  Fraction of not destroyed and removed that remains in residuals   \n",
       "17  Fraction of not destroyed and removed that remains in residuals   \n",
       "18  Fraction of not destroyed and removed that remains in residuals   \n",
       "19  Fraction of not destroyed and removed that remains in residuals   \n",
       "20  Combustion parameters                                             \n",
       "21  Combustion parameters                                             \n",
       "\n",
       "    Parameter Name  \\\n",
       "0   PFOA             \n",
       "1   PFOS             \n",
       "2   PFBA             \n",
       "3   PFPeA            \n",
       "4   PFHxA            \n",
       "5   PFHpA            \n",
       "6   PFNA             \n",
       "7   PFDA             \n",
       "8   PFBS             \n",
       "9   PFHxS            \n",
       "10  PFOA             \n",
       "11  PFOS             \n",
       "12  PFBA             \n",
       "13  PFPeA            \n",
       "14  PFHxA            \n",
       "15  PFHpA            \n",
       "16  PFNA             \n",
       "17  PFDA             \n",
       "18  PFBS             \n",
       "19  PFHxS            \n",
       "20  frac_vs_to_res   \n",
       "21  res_ts_cont      \n",
       "\n",
       "                                           Parameter Description  amount  \\\n",
       "0   DRE of PFOA                                                   0.97     \n",
       "1   DRE of PFOS                                                   0.97     \n",
       "2   DRE of PFBA                                                   0.97     \n",
       "3   DRE of PFPeA                                                  0.97     \n",
       "4   DRE of PFHxA                                                  0.97     \n",
       "5   DRE of PFHpA                                                  0.97     \n",
       "6   DRE of PFNA                                                   0.97     \n",
       "7   DRE of PFDA                                                   0.97     \n",
       "8   DRE of PFBS                                                   0.97     \n",
       "9   DRE of PFHxS                                                  0.97     \n",
       "10  Fraction of PFOA remained in residuals                        0.33     \n",
       "11  Fraction of PFOS remained in residuals                        0.33     \n",
       "12  Fraction of PFBA remained in residuals                        0.33     \n",
       "13  Fraction of PFPeA remained in residuals                       0.33     \n",
       "14  Fraction of PFHxA remained in residuals                       0.33     \n",
       "15  Fraction of PFHpA remained in residuals                       0.33     \n",
       "16  Fraction of PFNA remained in residuals                        0.33     \n",
       "17  Fraction of PFDA remained in residuals                        0.33     \n",
       "18  Fraction of PFBS remained in residuals                        0.33     \n",
       "19  Fraction of PFHxS remained in residuals                       0.33     \n",
       "20  Fraction of incoming Volatile Solids to Combustion Residuals  0.05     \n",
       "21  Solids content of residuals                                   0.85     \n",
       "\n",
       "        unit  Reference  \n",
       "0   fraction NaN         \n",
       "1   fraction NaN         \n",
       "2   fraction NaN         \n",
       "3   fraction NaN         \n",
       "4   fraction NaN         \n",
       "5   fraction NaN         \n",
       "6   fraction NaN         \n",
       "7   fraction NaN         \n",
       "8   fraction NaN         \n",
       "9   fraction NaN         \n",
       "10  fraction NaN         \n",
       "11  fraction NaN         \n",
       "12  fraction NaN         \n",
       "13  fraction NaN         \n",
       "14  fraction NaN         \n",
       "15  fraction NaN         \n",
       "16  fraction NaN         \n",
       "17  fraction NaN         \n",
       "18  fraction NaN         \n",
       "19  fraction NaN         \n",
       "20  fraction NaN         \n",
       "21  fraction NaN         "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ThermalTreatment = ps.ThermalTreatment()\n",
    "ThermalTreatment.InputData.Data[['Category','Parameter Name', 'Parameter Description', 'amount', 'unit','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming MSW to Thermal Treatment "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Parameter</th>\n",
       "      <th>Unit</th>\n",
       "      <th>Amount</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Mass flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>1000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Solids flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>254</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>746</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>VS flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>134.62</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>84.074</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>1070</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>840</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "        Parameter Unit  Amount\n",
       "0   Mass flow      kg   1000  \n",
       "1   Solids flow    kg   254   \n",
       "2   Moisture flow  kg   746   \n",
       "3   VS flow        kg   134.62\n",
       "4   Carbon flow    kg   84.074\n",
       "5   PFOA           μg   1070  \n",
       "6   PFOS           μg   840   \n",
       "7   PFBA           μg   0     \n",
       "8   PFPeA          μg   0     \n",
       "9   PFHxA          μg   0     \n",
       "10  PFHpA          μg   0     \n",
       "11  PFNA           μg   0     \n",
       "12  PFDA           μg   0     \n",
       "13  PFBS           μg   0     \n",
       "14  PFHxS          μg   0     "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncominWaste = ps.IncomFlow()\n",
    "IncominWaste.set_flow('MSW', 1000)\n",
    "IncominWaste.calc()\n",
    "MSW = IncominWaste.Inc_flow\n",
    "MSW.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in Thermal Treatment"
   ]
  },
  {
   "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>Air Emission</th>\n",
       "      <th>Combustion Residuals</th>\n",
       "      <th>Destructed</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>PFOA</th>\n",
       "      <td>2.01</td>\n",
       "      <td>0.99</td>\n",
       "      <td>97.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFOS</th>\n",
       "      <td>2.01</td>\n",
       "      <td>0.99</td>\n",
       "      <td>97.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFPeA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHpA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFNA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFDA</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFBS</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>PFHxS</th>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Air Emission  Combustion Residuals  Destructed\n",
       "PFOA   2.01          0.99                  97.0      \n",
       "PFOS   2.01          0.99                  97.0      \n",
       "PFBA  NaN           NaN                   NaN        \n",
       "PFPeA NaN           NaN                   NaN        \n",
       "PFHxA NaN           NaN                   NaN        \n",
       "PFHpA NaN           NaN                   NaN        \n",
       "PFNA  NaN           NaN                   NaN        \n",
       "PFDA  NaN           NaN                   NaN        \n",
       "PFBS  NaN           NaN                   NaN        \n",
       "PFHxS NaN           NaN                   NaN        "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ThermalTreatment.calc(Inc_flow=MSW)\n",
    "ThermalTreatment.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deXxNd/7H8de92XcUEWJChFiyiSSSIoISW2JtUaXWUlo1MxQ1U6q1tLSl1hlrtRppTQlNpaklogZJEEqCjDZKErHUdiMkuTm/PzK984ugokmOm/t5Ph55PNxzv9/v+ZzL4+3knO/5Xo2iKAghhKh6WrULEEIIUyUBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEvPyNK5du7bSqFGjSipFCCGqpyNHjlxVFKXO/dvLFcCNGjUiJSWl4qoSQggToNFozj9ou1yCEEIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJICFEEIlEsBCCKESCWBRoS5dusTgwYNp0qQJLVu2pGfPnpw9e/YPjTlixAi2bNnyh2vbsGED2dnZhtdjxowhLS3tD48rxJOSABYVRlEU+vXrR1hYGOfOnSMtLY158+aRm5urdmlA2QBes2YNLVu2VLEiYeokgEWF2bt3LxYWFowfP96wzc/Pj/bt2zN16lS8vLzw9vYmOjoagISEBDp27MgLL7xAs2bNmD59Ops2bSIoKAhvb2/OnTtnGGfXrl106NCBZs2a8c033wAlgfraa68Z2vTu3ZuEhAT0ej0jRoww7O/jjz9my5YtpKSkMHToUPz8/MjPzycsLMywvGpUVBTe3t54eXkxbdo0w5j29vbMnDkTX19fgoODn5r/TET1IAEsKszJkydp06ZNme1ff/01qampHD9+nF27djF16lRycnIAOH78OEuWLOHHH3/ks88+4+zZsyQlJTFmzBiWLl1qGCMzM5N9+/YRGxvL+PHjuXv37kPrSE1NJSsri5MnT/Ljjz8ycuRIBg4cSEBAAJs2bSI1NRUbGxtD++zsbKZNm8aePXtITU0lOTmZbdu2AZCXl0dwcDDHjx8nNDSU1atXV9THJYQEsKh8P/zwA0OGDMHMzAxnZ2c6duxIcnIyAIGBgbi4uGBlZUWTJk3o1q0bAN7e3mRmZhrGeOGFF9BqtTRt2hR3d3dOnz790P25u7vz008/8frrrxMXF4ejo+Mj60tOTiYsLIw6depgbm7O0KFDSUxMBMDS0pLevXsD0KZNm1I1CfFHSQCLCtOqVSuOHDlSZruiKA/tY2VlZfizVqs1vNZqtRQVFRne02g0pfppNBrMzc0pLi42bPvtrLhmzZocP36csLAwli9fzpgxYx5Z96Pqs7CwMOzbzMysVE1C/FESwKLCdO7cmXv37pX6NT05OZmaNWsSHR2NXq/nypUrJCYmEhQUVK6xv/rqK4qLizl37hw//fQTnp6eNGrUiNTUVIqLi7lw4QJJSUkAXL16leLiYgYMGMC7777L0aNHAXBwcOD27dtlxm7bti379u3j6tWr6PV6oqKi6Nix4x/4JIR4POX6TjghHkWj0bB161YmT57MggULsLa2plGjRixevBidToevry8ajYYPPviAevXqPfIywv08PT3p2LEjubm5rFq1Cmtra9q1a0fjxo0NN8/8/f0ByMrKYuTIkYaz4/nz5wMl09nGjx+PjY0NBw8eNIzt4uLC/Pnz6dSpE4qi0LNnT/r06VOBn4wQD6Z51K9f9wsICFDkSzmFEKJ8NBrNEUVRAu7fLpcghBBCJRLAQgihEglgIYRQiQSwqBJmZmb4+fnh5eVFREQEN27ceKJxHrZ+w/1PxZWXvb39E/cV4klJAIsqYWNjQ2pqKidPnqRWrVosX778icaR9RtEdSIBLKpcSEgIWVlZhtcLFy4kMDAQHx8fZs2aBZQ8AtyrVy98fX3x8vIyrB/x/9dvWL9+Pc2aNaNjx44cOHDAMN79q6f9dnar0+no0qUL/v7+eHt7ExMTU6a2nJwcQkNDDWfr+/fvr/gPQIj/knnAokrp9Xp2797N6NGjAYiPjycjI4OkpCQURSEyMpLExESuXLlC/fr1iY2NBeDmzZulxsnJyWHWrFkcOXIEJycnOnXqROvWrR+5b2tra7Zu3YqjoyNXr14lODiYyMjIUk/ZffHFF4SHhzNz5kz0ej137typ4E9AiP+RM2BRJfLz8/Hz8+OZZ57h119/pWvXrkBJAMfHx9O6dWv8/f05ffo0GRkZeHt7s2vXLqZNm8b+/ftxcnIqNd7hw4cN6zdYWloyaNCg361BURTeeustfHx8eO6558jKyiqzullgYCDr169n9uzZ/Pjjjzg4OFTchyDEfSSARZX47Rrw+fPnKSgoMFwDVhSFGTNmkJqaSmpqKv/5z38YPXo0zZo148iRI3h7ezNjxgzmzJlTZsz714f4zf9fI0JRFAoKCgDYtGkTV65c4ciRI6SmpuLs7FxmVbXQ0FASExNp0KABw4YNY+PGjRX5MQhRigSwqFJOTk588sknLFq0iMLCQsLDw1m3bh06nQ4oeYz48uXLZGdnY2try0svvcSUKVMM6zn8pm3btiQkJHDt2jUKCwv56quvDO81atTIsChQTEwMhYWFQMlljLp162JhYcHevXs5f/58mfrOnz9P3bp1GTt2LKNHjy6zXyEqklwDFlWudevW+Pr6snnzZoYNG0Z6ejohISFAyQ2zzz//nP/85z9MnToVrVaLhYUFK1euLDWGi4sLs2fPJiQkBBcXF/z9/dHr9QCMHTuWPn36EBQURJcuXbCzswNg6NChREREEBAQgJ+fH82bNy9TW0JCAgsXLsTCwgJ7e3s5AxaVStaCEEKISiZrQQghxFNGAlgIIVQiASxUFxcXh6enJx4eHixYsKDM+4qiMGnSJDw8PPDx8Sl1Y+xhfY8fP05ISAje3t5ERERw69YtAAoKChg5ciTe3t74+vqSkJBQ6ccnxEMpivLYP23atFGEqEhFRUWKu7u7cu7cOeXevXuKj4+PcurUqVJtYmNjle7duyvFxcXKwYMHlaCgoN/tGxAQoCQkJCiKoihr165V/va3vymKoijLli1TRowYoSiKouTm5ir+/v6KXq+vqsMVJgpIUR6QqXIGLFSVlJSEh4cH7u7uWFpaMnjw4DKPCMfExDB8+HA0Gg3BwcHcuHGDnJycR/Y9c+YMoaGhAHTt2pV//etfAKSlpdGlSxcA6tatS40aNZAby0ItEsBCVVlZWTRs2NDw2tXVtdQ6EY9q86i+Xl5ebN++HSj5PrkLFy4A4OvrS0xMDEVFRfz8888cOXLE8J4QVU0CWKhKecA0yPufcHtYm0f1XbduHcuXL6dNmzbcvn0bS0tLAEaNGoWrqysBAQFMnjyZZ599FnNzmQ4v1CH/8oSqXF1dS52BXrx4kfr16z9Wm4KCgof2bd68OfHx8QCcPXvWsKiPubk5H3/8saHPs88+S9OmTSv+wIR4DHIGLFQVGBhIRkYGP//8MwUFBWzevJnIyMhSbSIjI9m4cSOKonDo0CGcnJxwcXF5ZN/Lly8DUFxczHvvvcf48eMBuHPnDnl5eQB8//33mJuby/rCQjVyBixUZW5uzrJlywgPD0ev1zNq1ChatWrFqlWrABg/fjw9e/bk22+/xcPDA1tbW9avX//IvgBRUVGGBX/69+/PyJEjgZJgDg8PR6vV0qBBAz777DMVjlqIEvIoshBCVDJ5FFkIIZ4yEsDCaF29epXRr4wjKSlJ7VKEeCISwMIopaWl4eMfwL8OZtC5a3eiojarXZIQ5SYBLIxOYWEhz7bvwNWr1yguuAOWtrz44hAOHjyodmlClIsEsDA6FhYWfLM9hkUL5nL77GEG9OzCzp07f/dLOYV42sg0NGGU2rdvj7W1NQDNmjWje/fuKlckRPnJGbAQQqhEAlgIIVQiASyEECqRABZCCJVIAAshhEokgIUQQiUSwEIIoRIJYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgYbTy8/MBsLKyUrkSIZ6MBLAwWunp6QA4OzurXIkQT0YCWBitv896B4CIiAiVKxHiyUgAC6OUmJjI5UvZDBn2MjVq1FC7HCGeiASwMDpnzpyhY8eOALz15hSVqxHiyUkAC6Nx8+ZNlnzyCc2bNwcgKmozXl5elb7f5f9YzZKlyyp9P8L0yLciG5HZc97j0qXcCh3zHyurNljGvfpaufsUFhWSlXOJ776JMWyLiYkhMjKyIkt7qNjYncR9E0Njtz9V2T6FadAoivLYjQMCApSUlJRKLEc8TFFREc1aepOdp2DbIrTCxr2+e3WFjfV7NBbW1AgdVu5+99do71gDKxvbiirrd9268St2bV+g6PgODiQm4OPjU2X7FtWDRqM5oihKwP3b5QzYSJibm7P7u2/xD2yLRS1XbNzbVMi4jgF9KmScymTbNBilsACNpTrzfa01WswdapPn5EzXHr04cTRFpr6JCiEBbEQaN27Mjm1f0713H8wGvodl7T+pXVKVMHd6OsLOrmVHdDey6dYzgsMHErG2tla7JGHk5CackWnfvj3LFn+Ebvtc9Hduql2OybELGcTFQluGDh9BeS7fCfEgEsBGaMSIlxn78ovovlmAUlSodjkmRaPRYtd1ErsPn2D2nHfVLkcYOQlgI7Xw/QUEt2yMbs9KOROrYloLK+wjZvDR0hVER3+pdjnCiEkAGymtVsuW6C+oU3CJvOSv1S7H5Jjb18Ih4i1GjxuPzAwST0oC2IjZ2dmxKy4W5dRO7pw9qHY5JsfSuQm2nSfQvVcEFy9eVLscYYQkgI2cq6srcd9s586eFRTk/qR2OSbHtlkItOpO1x69ycvLU7scYWQkgKuBwMBA1qxaye0d89DrrqtdjsmxCxzAZfM6PD94KMXFxWqXI4yIBHA1MXjwICZPHIfum3kUF95TuxyTotFosO8ykYOnfmLajJlqlyOMiARwNTJn9iw6BXqTt2upzIyoYhpzC+x7T+cf6zfy6acb1S5HGAkJ4GpEo9GwaeMGGpjpyDu0We1yTI6ZrRP2ETOZ+MafOXDggNrlCCMgAVzN2NjYEP/tDrQZCeSlJ6pdjsmxrOOGXddJ9O7Tj8zMTLXLEU85CeBqqF69eny/M5b8fau5l31G7XJMjk2TAMxa96NLeE9u3bqldjniKSYBXE35+vqy6dP13P5mAUW3rqhdjsmx84/gpmNj+g0chF6vV7sc8ZSSAK7G+vTpw1tv/hXdjnkUF+SrXY5J0Wg02IW9wrHMy0z681/VLkc8pSSAq7kZ096kZ1gIed8tRlFkjmpV0piZY9fzTT77ciur/vEPtcsRTyEJ4GpOo9GwYe1qGjso5P3wmdrlmBwza3scImcyZdpb7N27V+1yxFNGAtgEWFpaEvdNDJYXktD9uEvtckyORa0G2PX4K30HPE9GRoba5YiniASwiahduza74r6l4N8buXvxlNrlmBwbN18s2w6hS7ceXL8uj4uLEhLAJqRly5Z8GbUJXexCCm9cUrsck2Pn2507zj707tOfwkJZSF9IAJuc7t27897st9Ftn0vxPVm9q6rZhY4g/fIdxr06UR4XFxLApuiNSa8xsFdXdDs/RCmWOapVSaM1w77HFP61cxcfL/5E7XKEyiSAK5C9vb0q+3377bfZtevxb65pNBpWrVhGC2c78vatK/XexZWjyF47kex1r5Eb/XfD8paG7etfJ3v969y9mG7ocyt5G+cX9St1Rl1ceJcrOxaW9Fk7gUufvylzkf9La2XLrZxMpk75K25ubvj6+vLRRx898VKW8+bNq7Dabty4wYoVK8rdb/bs2SxatKjC6jAVEsDVwJw5c3juuefK1cfCwoLYmK+xu3KKvNSdpd5zHjKP+qOWYVnPg5uHviy9feRS6o9cirVrC8P2vPR9WLk0LfWtHLdTtmNmW4P6o5dTf/QKnukxCY3W/AmPsPrRmFtRZ8h8rt/SsWzZMr799lveeeedJxrrYQGsKEq5Q/1JA1g8GQngSpCQkEBYWBgDBw6kefPmDB061HC9Lzk5mWeffRZfX1+CgoK4ffs2d+/eZeTIkXh7e9O6dWvDfNENGzbQt29fIiIiaNy4McuWLeOjjz6idevWBAcH8+uvvwIwYsQItmzZAkCjRo2YNWsW/v7+eHt7c/r0aQCuXLlC165d8ff3Z9y4cbi5uVFUVMTu+J0UJkeTn5la5jisG3pReD37kcdaeD0HpeAuNToMIy9tn2G7Pu865g7PGF5bPOOKxtziD3yq1Y+1a0us2r3MoKHDmD9/PsuWLUNRFPR6PVOnTiUwMBAfHx/+8d+HOHJycggNDcXPzw8vLy/279/P9OnTyc/Px8/Pj6FDh5KZmUmLFi2YMGEC/v7+XLhwodRvZlu2bGHEiBEA5Obm0q9fP3x9ffH19eXf//4306dP59y5c/j5+TF16lQAFi5caKhl1qxZhrHmzp2Lp6cnzz33HGfOyJojT0JOSSrJsWPHOHXqFPXr16ddu3YcOHCAoKAgBg0aRHR0NIGBgdy6dQsbGxuWLFkCwI8//sjp06fp1q0bZ8+eBeDkyZMcO3aMu3fv4uHhwfvvv8+xY8f485//zMaNG5k8eXKZfdeuXZujR4+yYsUKFi1axJo1a3jnnXfo3LkzM2bMIC4ujn/+858AeHh4sO1fX9G77wDMny99JpV/LgnLOo0Mr3Oj3gKtFo2ZBS7DPwJKzn5tW4Ri1bAVhbFZ6PNuYGZXA3vvruR++XfyzhzA2s0Pe6/OWNRqUBkftVGz8+qC7noW4197g+LiYi5fvkxMTAxOTk4kJydz79492rVrR7du3fj6668JDw9n5syZ6PV67ty5Q4cOHVi2bBmpqSX/gWZmZnLmzBnWr1//u2eykyZNomPHjmzduhW9Xo9Op2PBggWcPHnSMF58fDwZGRkkJSWhKAqRkZEkJiZiZ2fH5s2bOXbsGEVFRfj7+9OmTZtK/7yqGwngShIUFISrqysAfn5+ZGZm4uTkhIuLC4GBgQA4OjoC8MMPP/D6668D0Lx5c9zc3AwB3KlTJxwcHHBwcMDJyYmIiAgAvL29OXHixAP33b9/fwDatGnD119/bdjH1q1bgZKZEDVr1jS0DwsL46OFC5jy93dBKTYErWWdRtTqMMzQznnIPMxsnUrt6056InX6zUSj0WLbLIQ7Z37Awb83ls7uNBi3hrs/HyP/fCo5G/+Cy0uLsKjd8A98qtWTXfuXyIz9gLy8PBRFIT4+nhMnThh+q7l58yYZGRkEBgYyatQoCgsL6du3L35+fg8cz83NjeDg4N/d7549e9i4sWTxeDMzM5ycnMrMUY6Pjyc+Pp7WrVsDoNPpyMjI4Pbt2/Tr1w9bW1sAIiMjn/j4TZkEcCWxsrIy/NnMzIyioiIURUGj0ZRp+6jpSP9/HK1Wa3it1WopKip6ZJ/f9vt7+wB4ZexYTpxMY8WKFdQb9hHmDrUe2R6g4PLPFF7PJjf67yUbioswd6qHg3/vkhotbbD1fBZbz2fRaDTk/5QsAfwAGo0Wq6AhXE//N2vWrUdRFJYuXUp4eHiZtomJicTGxjJs2DCmTp3K8OHDy7Sxs7O7b/z//Zu7e/duuWpTFIUZM2Ywbty4UtsXL178wH/LonzkGnAVat68OdnZ2SQnJwNw+/ZtioqKCA0NZdOmTQCcPXuWX375BU9Pzwrdd/v27fnyy5IbavHx8Q98GmvJR4uwMjfjzoFPH2uOal76Pmq0exHXV9eV/EzcSJHuGkU3L3P3Yhr6uzoAFH0hBVcvYOZYt0KPqbrQ37nJ9b1rsG8dwYIPF1O3bl1WrlxpeFjj7Nmz5OXlcf78eerWrcvYsWMZPXo0R48eBUpuqD7qwQ5nZ2fS09MpLi42/BYE0KVLF1auXFlSg17PrVu3cHBw4Pbt24Y24eHhrFu3Dp2u5O8yKyuLy5cvExoaytatW8nPz+f27dvs2LGjwj8XUyBnwFXI0tKS6OhoXn/9dfLz87GxsWHXrl1MmDCB8ePH4+3tjbm5ORs2bCh15lsRZs2axZAhQ4iOjqZjx464uLjg4OBQqo2ZmRl169ZBq7uA7uh27Nr0eeSYeen7cX5+dqlttk2DyUtPxMy+Jr/GrwBFAUXBpkkAtp7tKvSYjJlSVED2+tdBrwetFrtWnXEM6kvBpY5s/GIOgwf2x9/fH0VRqFOnDtu2bSMhIYGFCxdiYWGBvb294fLBK6+8go+PD/7+/sydO7fMvhYsWEDv3r1p2LAhXl5ehjBdsmQJr7zyCmvXrsXMzIyVK1cSEhJCu3bt8PLyokePHixcuJD09HRCQkKAkqmWn3/+Of7+/gwaNAg/Pz/c3Nzo0KFD1X141YimPE/jBAQEKCkpKZVYjqgs9+7dw8zMDHNzcw4ePMirr75quNFyv/Pnz9M6IAiLsFexaRJYxZWKvPT9aJI+4/jRFFxcXNQuR1QAjUZzRFGUgPu3yxmwifjll1944YUXKC4uxtLSktWrVz+0rZubG7Hbt/Fc916YDXy31EwIUfnsWnRAdyOLbj0jSE05jJmZmdoliUoiAWwimjZtyrFjxx67fUhICP37RvBtaiyWXSdWYmXiQTSOzlw5t0/Wi6jm5CaceKCtW7ey9Zs4bIIHqV2Kybl7MZ27P6xn93c7MTd/+DlSXFwcnp6eeHh4sGDBgjLvX79+nX79+uHj40NQUBAnT54E4MyZM/j5+Rl+HB0dWbx4caUdj3g4CeBqaNSoUdStWxcvL68n6n/s2DGGjRyNQ+8ZmDvUruDqxKMU3byM7tsPiN70Ga1atXpoO71ez8SJE9m5cydpaWlERUWRlpZWqs28efPw8/PjxIkTbNy4kTfeeAMAT09PUlNTSU1N5ciRI9ja2tKvX79KPS7xYBLA1dCIESOIi4t7or45OTl069kbm7BxWLk0reDKxKMU37vD7e3vMXvmdHr27PnItklJSXh4eODu7o6lpSWDBw8mJiamVJu0tDS6dOkClEyBzMzMJDc3t1Sb3bt306RJE9zc3Cr2YMRjkQCuhkJDQ6lV6/cfpLhffn4+XXv0RvHsgl3z9pVQmXgYpViPLu5D+oV34q9/+fPvts/KyqJhw/891OLq6kpWVlapNr6+voYnIZOSkjh//jwXL14s1Wbz5s0MGTKkAo5APAkJYAGUPPE0aOgwcpQa2LV9Qe1yTE7e/k/xrG3N6n+sfKwnzB50c+7+ftOnT+f69ev4+fmxdOlSWrduXeqackFBAdu3b+f555//4wcgnojMghAA/O3tWew/dhrHAe/KI6ZVLO/Ed9jkpBJ7NBkLi8dbMc7V1ZULFy4YXl+8eJH69euXauPo6Mj69euBksBu3LgxjRs3Nry/c+dO/P39cXZ2roCjEE9CzoAFX3wRxSer1mDfezoac0u1yzEpd385QcGhL9gdv7Ncl40CAwPJyMjg559/pqCggM2bN5dZEOfGjRsUFBQAsGbNGkJDQw0LQAFERUXJ5QeVyRmwiTt8+DCvTHgNx/7vYGZX8/c7iApTeD0b3c4P2fbVZpo1a1auvubm5ixbtozw8HD0ej2jRo2iVatWrFq1CoDx48eTnp7O8OHDMTMzo2XLlqxdu9bQ/86dO3z//feGtYaFOuRR5GpoyJAhJCQkcPXqVZydnXnnnXcYPXp0mXa//PILrQOCMAt9BVuPtipUarr0d3Xcin6T+W9PZ8Krr6pdjqhk8iiyCYmKivrdNjqdjue690Lr3VvCt4op+iLyvl3Ii/37SPiaOLkGbIKKi4sZ8MIQrtk0wDagr9rlmBRFUchLWIPPn55h2Scfq12OUJmcAZugv059k+SzF3DoN1tmPFSxO8dicbyRwfbdSbLIjpAANjVr161j7efROA3+AI2ZfElmVcr/6QhFR//F7uTDpWYjCNMlAWxCEhMTmfSXqdR4fi5mNhIAVang6i/kxS8h7puYUnNxhWmTa8Am4ty5c0T2G4B9+GQsnpHvZatK+js30W2fy7LFH9G+vTziLf5HAtgE3Lx5k+fCe2IR8Dw2jf3VLsekKEWF6GLfZ+zLLzJixMtqlyOeMhLA1VxRURER/QZyq5Yndq17qV2OSVEUhbw9qwhu0YiF75ddr1cIuQZczU2c9AYns27g0Od1tUsxOXdStvLMvWy2RB9Eq5VzHVGWBHA1tmz5cjZvjcVp0PtotDLlqSrdyThM8clv2X0kGTs7O7XLEU8pCeBqKj4+nul/m4XTC/PRWturXY5JKcj9iTu7l7H3++9wdXVVuxzxFJPfi6qh06dPM3DQi9j3mIJFzfq/30FUGL3uOrd3zGPNqpUEBgaqXY54ykkAVzPXrl3jue49sQoZivWfvNUux6QUF95D9808Jk8cx+DB8mWm4vdJAFcjBQUF9Izsx936/tj5dFO7HJOiKAp5u5YSFuDFnNmz1C5HGAm5BlxNKIrC6FfGk/FrIQ4Rw9Uux+TkHYqmgZmOLz77VtbXEI9NAriaWLjoQ3bs2o/j8/NkxkMVy0vfjzZjL/FHU7CxsVG7HGFEJICrgR07djBn/vvUGPQBWitbtcsxKfeyz5C/758c2LeXevXqqV2OMDJyDdjInThxgheHvYxDr+mYO9VVuxyTUnTrCrrY99n06Xp8fX3VLkcYIQlgI5abm0vXHr2wDh2NVYPmapdjUooL7qLbMY8ZU/9Cnz591C5HGCkJYCN19+5dwntFonfvgF3LMLXLMSmKUkzedx/TMyyEGdPeVLscYcQkgI2Qoii89PJILtyzxu5Z+Vrxqpb3w+c0dlDYsHa1zHgQf4jchDNC77z7HrsOpeI48D00Gvk/tCrlndyN5YXDxB1NwdLSUu1yhJGTADYyX375FYuWLKfG4A/QWlirXY5JuXvxFPcOfMq+A/upXbu22uWIakAC2IikpKQw6pVxOPSdhbnDM2qX89gurhiB/vbVSt2HpXUlzr/VmOHYdQL5ievYErWJli1bVt6+hEmRADYSRUVF9Izog96mFoXpeylM31sh415LiqmQcR6XU6vQcvcpypa5YksAAAkhSURBVNdRVHAPtFr0uecouneHOe++y5DBgwGwsrKiZs2aFV2qQb8XhrAn9iMWLfqQ7t27V9p+hOmRADYS5ubmLP5wIVeuXKnQcSdXYQDb2dnxztj+5e5XVFTEpdzLfB4VzT3nhhTVacYHi5ej1Zoz863plVBpaQFt/HFv5MbkN2RRe1GxNIqiPHbjgIAAJSUlpRLLEeLhiouL+e677xgzfgJ3ardAf+EEH86bzdgxY9QuTYhH0mg0RxRFCbh/u9xCF0ZDq9XSo0cPjiUfxvF6BtpGgUz+81+5deuW2qUJ8UQkgIXRqVu3LrHbt1J0dj/WDVvx2WefqV2SEE9EAlgYJS8vL7y8WlFg6cSiJUspz6U0IZ4WEsDCaL3+6itYF90m91Iuly5dUrscIcpNAlgYrUaNGkH+Tazsnbh+/bra5QhRbhLAwmhZWVmhFBehNbeksLBQ7XKEKDcJYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQiQSwEEKoRAJYCCFUIgEshBAqkQAWQgiVSAALIYRKJICFEEIlEsDCKG3atIkvv/zS8Hrz5s1ERUWpWJEQ5SdfSy+MTmFhIW7uHly+eo3a4RNB0XP1u1W41HMm89xZzMzM1C5RiFLka+lFtWFhYcGxlCRatWxJ4blDFJ47jK+vD0eTD0n4CqMiASyMkrOzM4cPJNIroCmRwS35d+Je6tSpo3ZZQpSLBLAwKlu3bkWj0XD69Gmsra1ZMHcOt29cw8rKSu3ShCg3CWBhVKKiomjfvj2bN28GoH79+mzZsqVMu6KioqouTYhykwAWRkOn03HgwAHWrl1rCODMzEy8vLwA2LBhA88//zwRERF069ZNzVKFeCzmahcgxOPatm0b3bt3p1mzZtSqVYujR49Sq1atUm0OHjzIiRMnymwX4mkkZ8DCaERFRTF48GAABg8e/MB5v127dpXwFUZDzoCFUbh27Rp79uzh5MmTaDQa9Ho9Go2GCRMmlGpnZ2enUoVClJ+cAQujsGXLFoYPH8758+fJzMzkwoULNG7cmIsXL6pdmhBPTAJYGIWoqCj69etXatuAAQOYN2+eShUJ8cfJo8hCCFHJ5FFkIYR4ykgACyGESiSAhRBCJRLAokqZmZnh5+dn+FmwYMFj983OzmbgwIHl3ueYMWNIS0srdz8hKpvchBNVyt7eHp1Op3YZQlQpuQknnmqNGjXirbfeIiQkhICAAI4ePUp4eDhNmjRh1apVQOl1H06dOkVQUBB+fn74+PiQkZFBXl4evXr1wtfXFy8vL6KjowEICwvjtxOHqKgovL298fLyYtq0aYb929vbM3PmTHx9fQkODiY3N7eKPwFhiiSARZXKz88vdQnit5AEaNiwIQcPHqRDhw6MGDGCLVu2cOjQId5+++0y46xatYo33niD1NRUUlJScHV1JS4ujvr163P8+HFOnjxJ9+7dS/XJzs5m2rRp7Nmzh9TUVJKTk9m2bRsAeXl5BAcHc/z4cUJDQ1m9enXlfhBCII8iiypmY2NDamrqA9+LjIwEwNvbG51Oh4ODAw4ODlhbW3Pjxo1SbUNCQpg7dy4XL16kf//+NG3aFG9vb6ZMmcK0adPo3bs3HTp0KNUnOTmZsLAww8LtQ4cOJTExkb59+2JpaUnv3r0BaNOmDd9//31FH7oQZcgZsHhq/LaoularLbXAularLbO+74svvsj27duxsbEhPDycPXv20KxZM44cOYK3tzczZsxgzpw5pfo86n6HhYUFGo0GKLlRKOsJi6ogASyM0k8//YS7uzuTJk0iMjKSEydOkJ2dja2tLS+99BJTpkzh6NGjpfq0bduWffv2cfXqVfR6PVFRUXTs2FGlIxBCLkGIKvbbNeDfdO/evVxT0X4THR3N559/joWFBfXq1ePtt98mOTmZqVOnotVqsbCwYOXKlaX6uLi4MH/+fDp16oSiKPTs2ZM+ffr84WMS4knJNDQhhKhkMg1NCCGeMhLAQgihEglgIYRQiQSwUF1cXByenp54eHg88IacoihMmjQJDw8PfHx8DLMbLly4QKdOnWjRogWtWrViyZIlhj6//vorXbt2pWnTpnTt2pXr168DJV9t1KlTJ+zt7Xnttdeq5gCFeAgJYKEqvV7PxIkT2blzJ2lpaURFRZVZOGfnzp1kZGSQkZHBP//5T1599VUAzM3N+fDDD0lPT+fQoUMsX77c0HfBggV06dKFjIwMunTpYgh2a2tr3n33XRYtWlS1ByrEA0gAC1UlJSXh4eGBu7s7lpaWDB48mJiYmFJtYmJiGD58OBqNhuDgYG7cuEFOTg4uLi74+/sD4ODgQIsWLcjKyjL0efnllwF4+eWXDY8c29nZ0b59e6ytravwKIV4MAlgoaqsrCwaNmxoeO3q6moI0fK0yczM5NixY7Rt2xaA3NxcXFxcgJL5v5cvX66sQxDiiUkAC1U9aB76b48EP24bnU7HgAEDWLx4MY6OjhVfpBCVRAJYqMrV1ZULFy4YXl+8eJH69es/dpvCwkIGDBjA0KFD6d+/v6GNs7MzOTk5AOTk5FC3bt3KPAwhnogEsFBVYGAgGRkZ/PzzzxQUFLB582bDqmi/iYyMZOPGjSiKwqFDh3BycsLFxQVFURg9ejQtWrTgL3/5S5k+n376KQCffvqpPHIsnk6Kojz2T5s2bRQhKlpsbKzStGlTxd3dXXnvvfcURVGUlStXKitXrlQURVGKi4uVCRMmKO7u7oqXl5eSnJysKIqi7N+/XwEUb29vxdfXV/H19VViY2MVRVGUq1evKp07d1Y8PDyUzp07K9euXTPsz83NTalZs6ZiZ2enNGjQQDl16lQVH7EwNUCK8oBMlbUghBCikslaEEII8ZSRABZCCJVIAAshhEokgIUQQiUSwEIIoRIJYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihEglgIYRQSbmWo9RoNFeA85VXjhBCVEtuiqLUuX9juQJYCCFExZFLEEIIoRIJYCGEUIkEsBBCqEQCWAghVCIBLIQQKpEAFkIIlUgACyGESiSAhRBCJRLAQgihkv8DTQBVbo9sa14AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ThermalTreatment.plot_sankey()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdeZxN9R/H8dd3NjNmi4RBWbKvQ4iUJfs2kYhky5IW7UIqbZZoUSRpIb+EKNmFsrXIOsoSUvYlss2+nt8f95oGY8wwM2fu3Pfz8fCouffcc9/33pnP/Zxzvud8jWVZiIiIiIi4Ag+7A4iIiIiIZJSaVxERERFxGWpeRURERMRlqHkVEREREZeh5lVEREREXIaaVxERERFxGWpes5gxxtcYYxljStidJSOMMeWNMWftziEiYhfVbRHX4hbNqzEmMtW/ZGNMTKqfu1/lsa2MMX9mYZb1xpgHs2p918uyrD2WZd2QHes2xswyxsQ53+fTxphlxphyzvvGGGMSLvlsnkj12I7GmE3GmChjzCljzHRjTEgaz9HK+aXzxKX3XSWbvzHma2PMAefj611l+ZuMMQudef42xtyX6r6+l7yOaOc6q6TxPlz4tyEzeZ3rqWOMCXeuf4Mxpmqq+3oaY/YYY84ZY04YYz41xvhn9jlEcgvV7Stz1brt3Eh43xhzxPnYv4wxb2YiW21jzM/OOnfIGDM0A4/xNcbsu/T3wRhzqzFmnbOe7jDGNEx130BjTOIlrzPSGHNjRrM615Pe98ZdxpgfnO/xSWPMTGNM4cys3525RfNqWVbAhX/AQaB9qttm2J0vj3vd+b7fAkQAH6e67/PUn41lWe8DOL+YpgJjgRuB6oAnsNYYE3TJ+nsBp53/zQwLWAM8AJzJwPJTnMsVBvoCn10o6JZlfXrJ79gzwC7Lsnakevzrl7zWumk9ifNL8rJG2hjjB8x35igAzAHmGWO8nIusAepblhUMlAWCgBEZeF0iuZLqtq2yq26PACoBtYBAoBnwWyZyfQUsBQo6H/usMabFVR4zHDiSxu1zgXXOdb0BfGuMSb1BsPqS1xlgWda/l67E2dBfqYm+4vcGcAMwESgJlAKSnctLBrhF83o1xhg/Y8wHxphjxpjDxphxxhhv51bWPKBM6i0vY0wDY8yvzq2/o8aYd1M1EZl53lbGmD+NMS84t7yOpN6jYBx7B993bmGeM8asufA8xphOxpidxpizxpiVqf4gMMYcN8Y849yajDTGfGiMCTHGrDDGnHduSQc5l61ojElM9dj1xpgRzv+eN8YsMcYUSHV/P2PMQWfe553PdefVXqtlWZHALKBqess5X99bwMuWZX1lWVasZVlHgZ7ORR5PtWwwcA8wEKhhUu2JzECeaMuy3rcs6yccRSO9TAWA9sBLlmVFWZb1A/AdcKW9P72AzzOaJYOaA7GWZU2yLCsOeBtH8b8TwLKsA6kKq8HxmsqmtaJLP3PnbSl7lowxXs7fu3+deyyeuHR5Ebupbrtk3a4DfG1Z1gnL4a+MbogYYzyBm4EZlmUlWZa1G1gPVEnnMeWBjjjqZerbqwPlcTTpsZZlzQT2AR0ykiWDedP93rAsa6FlWd9YlhVhWVYU8AHQIJ31XfSZGUfT/Emqn6/pM3ZVal4dXsWxlVgNuA1oDDzvbAY6An9dsuWVgOOPsSBwF45f0H7X+NwlcTQbxZzrnGyMCXDe9z5QEccffEHgRcByNmnTgEdxbNGtARZcUog7Ao2AykBXHHvtngGKAAHAI+lkegDHH1gIjq3DJwGMMaHAO0AXoITzX6GMvEhn0e0GbL3KolWBojj2LKawLCsJ+AZHE3dBF+Akji3oVfxXKC88525jzL0ZyXcVFYEIy7IOpLptG2kUTWexrAN8kQXPm1oV53MCYFlWMrA9dQZjTFNjzDngHNAGx+/PtXgcx+9OVaAucF/6i4vYQnX7Yq5Qt9cDQ4zjsHxa9fNTY8w7aT2Bc10TgF7OjZQqOD7379PJNQkYDMRdcnsVYI9lWTGpbkuzpl+HDH9vODUEdlzhvnRdz2fsqtS8OnQHRliWdcqyrBM4DiH0uNLClmVtsCxro3Prbx/wCY6Ccy2igdGWZSVYljUPx+HsssYYbxzN2CDLso47n2ud8w+4GzDPsqzVlmXFA6Nw/KLWTrXe8c7XcxD4GfjJsqzfnX+s84Ga6WT62LKsfc6twblAqPP2Lji2mtc79/69yNV/h4Ybx4kFu3EcQkr9ZdHDuQfiwr+CztdhASfSWNcxLv6D7AXMtCzLAr4EHnRunQNgWVYFy7K+uUq+jAjA0RCmdg7Hns9L9QJWWpZ16WGq4Ze81o+yOoNlWd87hw3cgqOQHczkc1zQBXjHsqxjzi/9sde4HpHspLp9MVeo268C43HUyS3OPebdLixoWVZfy7KeSSfXtzje3xgcG+8TLctKc9iBc70RlmUtTePujNT0Rpe8zsw2lhn+3jDG3AYMcf67FtfyGbu0PP3iMsIYY3BsMabeOjoAFE/nMZWNMUuN48SY88DLXPtWzknnXrQLonH80ocAXsBfaTymWOq8zsJ45JLMqYtITBo/B3Blx9PIc+F5D6V63vNc/sd5qZGWZd1gWVaIZVkdL9kK/Z/zvgv/TgOncOzRKJLGukKc92OMuRXHIZYLh5y+xjEWtHkaj7tekTjGkKYWhGMsWArn79KDpD1kYOQlr/Vh52PypS6QOL7IVqS67anMZACwLOsQjj3R17r396LP+ZL/F7Gd6naacn3ddjb771mWVR9HvX4HmO6s5+kyxhQBlgDDgHw4xol2MsY8lMayQTg2Zp6+wuoyUk/XXPI6Ux/lWpGqZj8FjEhVs+dm4jkwxlQCFgEPW5b16xXyXs21fMYuze2bV+deu+M4DgNdcAv/DfC20njYx8AW4FbLsoKA13D84WalY0AiUCaN+46SKq9zb2Nx0h6UntWZUi4l4ywQwVn8HNtxFOzOqW90vsaO/HeI6MIJWiuMMceBPTi+NC4aOpBF/gCCjDG3pLqtBpcf4rkbx/vxbUZXbFlWXOoCCWwCmqe6bbxz0R3O5wTAGOOB41DdlfYGeAHpfSF4OPcSXZD6RIWLPmcc48xEcg3V7Uxnyi11O4XlOO/gHRyH9Ctm4DnK4diTOsu5R/sAjmEKbdJYtjKOurXe+f3wJVDKOQ40BEfdLG+M8U31mLRqeposy2qeqmaPB15NVbMvDLO66veGs2lfAbxgWdZXGXjq1HmvWLOz6TPOVdy+eXWaiWPL6UbjuFTFcP7ba3UCKJxqPBM4dvufsywr0jnupn9WB7IsKwGYDrxnjClijPE0xtzpLAazgY7GmIbOBmQo8C+Oxic7fYVjS7eOMcYHR/FP92SnzLIsKxHHoZPXjTGdnXsmi+HYm+kFTHTudekBvIDj0NiFfw8AHczlVyRIk3PdF4qBzyWFLHWmMzi2jF8zxuQ3xjQGWvHfXt8LegFfXTKOKqusAPycY8Xy4dijEAX8CCmXyirh/P/SOD6b9MaCGeAhY4yHMaYdUAEIdI6/+wp42hhT1DhOfnkuG16PyPVS3c6YXFG3AYwxzxrHJaJ8jWPc6gAcwxK2XXHl/9kF+Btj7jMOxXGMx0/rsZtxbMxc+G54DMeeyVDghHOowV7gJWfWLjhOcJ1/zW/CJa72vWGMKQn8ALxpWdbUDK62pzHGxxhTFsfJyoHO36Vs/4xzGzWvDi8DO3FsEYUDP/HfOL9twALggPlvfM/TQD9jTCSOMwRnZ1OuJ3CcAbkVR5F7HTDOP7y+wEc4TlhqCtzjLCDZxrKsrTgGv8/DsbfgGI5DE5cOhr/e5/kcx+sbiuMyI78777rTsqyzOMapFQY+dI4rO25Z1nEc47yO4Bj/g3GcKd8pnac6gONQ3I04Tp6IMcYUdT72VWPMvFTL9sdxiPEUjpMu+lqWtffCnc4vyXu58lUGXjIXXy/wcAbfDgCcDfGFKyucxXEyR4dUn3l1YIMxJgrH5V+24Tgx5EqSnI85iaOwvwGMBErj+KL5GcffxEYcBThLP2ORLKC6nQG5qG7jfM73cWxc/AP0wVHHDgMYY6YZY8aTBuf4+844hg2cxdGgbgTGOR/bzBiTenhC6u+GM0CS8+cLTV1nHCdJnQVeATo6G84LGpvLr/NaLZNvS3rfGwNxNNijU63/1FXWF4fje2s+jt+r5kCXnPqMcxPjOPoiknnGcSmQ00Axy7KO2Z1HMsYYUxHYbllWhi4TZIzpCIyxLKtC9iYTkeymuu2anMMf7rMs68cMLJvnP2PteZVMMcaEGcf1FQNwDLb/Na/+cbgrY0ygMaaF85DnLTjOXJ13tceJSO6kup33udtnrOZVMqszjhMlDuM42SDdaRrFJXkAY3AcdtqI4ySXN2xNJCLXQ3U773Orz1jDBkRERETEZWjPq4iIiIi4jEw1r8aYZdkVREQkO7lj/XLH1ywieUN69StDZxtfEBQU1LJ27doaZyAirui83QFymmq2iLiwK9bsTDWv5cqVY9Om7L6esohI1jPG7L36UnmLaraIuKr0arbGvIqIiIiIy1DzKiIiIiIuQ82riIiIiLgMNa8iIiIi4jLUvIqIiIiIy1DzKiIiIiIuQ82riIiIiLgMNa8iIiIi4jLUvIqIiIiIy1DzKiIiIiIuQ82riIiIiLgMNa8iIiIi4jLUvIqIiIiIy1DzKiIiIiIuQ82riIiIiLgMNa8iIiIi4jLUvIqIiIiIy1DzKiIiIiIuQ82rZIvjx4/TtWtXbr31VipXrkybNm3Ys2fPda2zd+/ezJ0797qzTZs2jaNHj6b83K9fP3bu3Hnd6xURcWWq2+Iq1LxKlrMsi44dO9K4cWP27dvHzp07GTVqFCdOnLA7GnB5Efzkk0+oXLmyjYlEROylui2uRM2rZLlVq1bh7e3NwIEDU24LDQ3lzjvvZPDgwVStWpVq1aoxe/ZsAFavXk2jRo3o0qUL5cuXZ+jQocyYMYO6detSrVo19u3bl7KelStXctddd1G+fHkWLVoEOIra448/nrJMu3btWL16NUlJSfTu3Tvl+d59913mzp3Lpk2b6N69O6GhocTExNC4cWM2bdoEwMyZM6lWrRpVq1ZlyJAhKesMCAhg+PDh1KhRg3r16uWagi4ikhVUt8WVqHmVLLd9+3Zuu+22y27/5ptvCA8PZ9u2baxcuZLBgwdz7NgxALZt28Z7773H77//zv/+9z/27NnDhg0b6NevHxMmTEhZx/79+1mzZg2LFy9m4MCBxMbGXjFHeHg4R44cYfv27fz+++/06dOH++67j9q1azNjxgzCw8Px8/NLWf7o0aMMGTKEH374gfDwcDZu3Mi3334LQFRUFPXq1WPbtm00bNiQjz/+OKveLhER26luiytR8yo55scff6Rbt254enpSpEgRGjVqxMaNGwGoU6cOISEh5MuXj1tvvZUWLVoAUK1aNfbv35+yji5duuDh4UG5cuUoU6YMf/zxxxWfr0yZMvz1118MGjSIZcuWERQUlG6+jRs30rhxY2666Sa8vLzo3r07a9euBcDHx4d27doBcNttt12USUQkr1LdltxIzatkuSpVqrB58+bLbrcs64qPyZcvX8r/e3h4pPzs4eFBYmJiyn3GmIseZ4zBy8uL5OTklNsubNUXKFCAbdu20bhxYz744AP69euXbu708nl7e6c8t6en50WZRERcneq2uBI1r5Ll7r77buLi4i46RLNx40YKFCjA7NmzSUpK4uTJk6xdu5a6detmat1z5swhOTmZffv28ddff1GhQgVKlSpFeHg4ycnJHDp0iA0bNgBw6tQpkpOT6dSpE6+//jpbtmwBIDAwkIiIiMvWffvtt7NmzRpOnTpFUlISM2fOpFGjRtfxToiIuAbVbXElXnYHkLzHGMO8efN46qmnGDNmDL6+vpQqVYrx48cTGRlJjRo1MMYwduxYihYtmu4hpEtVqFCBRo0aceLECSZPnoyvry8NGjSgdOnSKQP2a9WqBcCRI0fo06dPytb96NGjAcelWwYOHIifnx+//PJLyrpDQkIYPXo0TZo0wbIs2rRpwz333JOF74yISO6kui2uxKS3y/1StWvXti6c3Sci4kqMMZsty6ptd46cpJotIq4qvZqtYQMiIiIi4jLUvIqIiIiIy1DzKjnK09OT0NBQqlatSvv27Tl79uw1redKUwNeeuHrzAoICLjmx4qI5DWq2ZIbqXmVHOXn50d4eDjbt2+nYMGCfPDBB9e0Hk0NKCKS/VSzJTdS8yq2qV+/PkeOHEn5edy4cdSpU4fq1aszYsQIwDFDStu2balRowZVq1ZNmZow9dSAU6dOpXz58jRq1IiffvopZX29e/dm7ty5KT9f2EKPjIykadOm1KpVi2rVqjF//vzLsh07doyGDRum7HFYt25d1r8BIiIuRDVbcgtdKktskZSUxPfff0/fvn0BWL58OXv37mXDhg1YlkVYWBhr167l5MmTFCtWjMWLFwNw7ty5i9Zz7NgxRowYwebNmwkODqZJkybUrFkz3ef29fVl3rx5BAUFcerUKerVq0dYWNhFF9L+8ssvadmyJcOHDycpKYno6OgsfgdERFyHarbkJtrzKjkqJiaG0NBQbrzxRk6fPk3z5s0BRyFcvnw5NWvWpFatWvzxxx/s3buXatWqsXLlSoYMGcK6desIDg6+aH2//vprytSAPj4+3H///VfNYFkWL7zwAtWrV6dZs2YcOXKEEydOXLRMnTp1mDp1Kq+88gq///47gYGBWfcmiIi4CNVsyY3UvEqOujB+6sCBA8THx6eMn7Isi2HDhhEeHk54eDh//vknffv2pXz58mzevJlq1aoxbNgwXnvttcvWeenUgxeknn7Qsizi4+MBmDFjBidPnmTz5s2Eh4dTpEiRlKkJL2jYsCFr166lePHi9OjRg+nTp2fl2yAi4hJUsyU3UvMqtggODub999/nrbfeIiEhgZYtW/LZZ58RGRkJOGZZ+eeffzh69Cj58+fnwQcf5LnnnkuZKvCC22+/ndWrV/Pvv/+SkJDAnDlzUu4rVapUylzd8+fPJyEhAXAcxipcuDDe3t6sWrWKAwcOXJbvwIEDFC5cmP79+9O3b9/LnldExJ2oZktuojGvYpuaNWtSo0YNZs2aRY8ePdi1axf169cHHAP1v/jiC/78808GDx6Mh4cH3t7efPjhhxetIyQkhFdeeYX69esTEhJCrVq1SEpKAqB///7cc8891K1bl6ZNm+Lv7w9A9+7dad++PbVr1yY0NJSKFStelm316tWMGzcOb29vAgICtBUvIm5PNVtyC00PKyJuQdPDioi4Dk0PKyIiIiJ5gppXEREREXEZal5FRERExGWoeZVcb9myZVSoUIGyZcsyZsyYy+63LIsnnniCsmXLUr169YvOMn3ooYcoXLgwVatWzcnIIiJuSzVbspuaV8nVkpKSeOyxx1i6dCk7d+5k5syZ7Ny586Jlli5dyt69e9m7dy9TpkzhkUceSbmvd+/eLFu2LKdji4i4JdVsyQlqXiVX27BhA2XLlqVMmTL4+PjQtWvXy+a1nj9/Pj179sQYQ7169Th79izHjh0DHBeuLliwoB3RRUTcjmq25AQ1r5KrHTlyhJtvvjnl5xIlSnDkyJFMLyMiItlPNVtygppXydXSug7xpVMLZmQZERHJfqrZkhPUvEquVqJECQ4dOpTy8+HDhylWrFimlxERkeynmi05Qc2r5Gp16tRh7969/P3338THxzNr1izCwsIuWiYsLIzp06djWRbr168nODiYkJAQmxKLiLgv1WzJCWpeJVfz8vJi4sSJtGzZkkqVKtGlSxeqVKnC5MmTmTx5MgBt2rShTJkylC1blv79+zNp0qSUx3fr1o369euze/duSpQowaeffmrXSxERyfNUsyUnmLTGnlyJ5skWEVeV3jzZeZVqtoi4qvRqtva8ioiIiIjLUPMqIiIiIi5DzauIiIiIuAw1ryIiIiLiMtS8ioiIiIjLUPMqLu/IkSOE3duZpUuX2h1FRESuIjY2locfHcSb495Kc7YtkatR8youbdOmTdSoVZt1h+Pp/EBP3h3/voqhiEgudeLECerd2Yg5a7cx+v2PeaBHL+Li4uyOJS5Gzau4rLlzv+aOO+8iodQdBNzeBb+7evP8sGH0GzCQ5ORku+OJiEgqO3bsoHrN2/jzTBJ+d/bCr/EAvv1uFQ0a3c3p06ftjicuxMvuACLXave+vyhfuSoRp37n4MRvCC5YiAqVKnPoyFGMMXbHExGRVI4cPcZNRUO4MSmZXR/1wzufH+UrVyEhIYH4+Hi744kL0Qxb4vKWLVtG69at6dG7L9OnfmJ3HMmlNMOWSO4QGRlJYGAgN5cpx8F9e+yOI7mUZtgSERERkTxBzauIiIiIuAw1ryIiIiLiMtS8ioiIiIjLUPMqIiIiIi5DzauIiIiIuAw1ryIiIiLiMtS8ioiIiIjLUPMqIiIiIi5DzauIiIiIuAw1ryIiIiLiMtS8ioiIiIjLUPMqIiIiIi5DzauIiIiIuAw1ryIiIiLiMtS8ioiIiIjLUPMqIiIiIi5DzauIiIiIuAw1ryIiIiLiMtS8ioiIiIjLUPMqLi8hIQEALy9Pm5OIiMjVXKjZnp5eNicRV6XmVVzemTNnAChc6Eabk4iIyNWcPXsWgOAbgm1OIq5Kzau4vEVLlwNQoXw5m5OIiMjVrFixEoBKFcrbnERclZpXcWmnTp1izqwZAHTu3NnmNCIikh7Lshj37nsAPPnYIzanEVel5lVc2ieffgbAgEceIyAgwOY0IiKSnl9//ZU//9hB2YpVuP322+2OIy5Kzau4rN9//51hQ4cA8NSgx2xOIyIi6Tl37hzdevQCYPDTT2KMsTmRuCo1r+KStm3bRv0GdwHw/NAXqFSpks2JRETkSv79918aN2vJ/j/30KDR3fTp09vuSOLC1LyKSzl8+DDDX3yJ0NBQoiLO8WCvPowZ9YbdsRg87EUe6v8wlmXZHUVEJNeIiIhg0qQPqVilGuGbfqVytVCWLJiHt7e3rbkWLFhAvTsbER0dbWsOuTZ59iJrhw4dYtWqVXbHkCyQmJjI2bNnWf7DGr5bvCDl9uEvvczrr76SKw497dy5i2VLFlKpQnkGP/es3XFEXE5cXByzZ8+2O4ZkgeTkZM6fP0/4b9uZ+unHKbe3bNOer7+aib+/v43pHA4cOMCG9b/QuVt3Fs77Gg8P7ctzJSYze4pq165tbdq0KRvjZJ133h3Pc88+y43VG2F0IWSXdnLListu69u3LyVLlrQhTdr+N2MmJ0NuJ2HbImZNn0q7du3sjiSXMMZstiyrtt05cpIr1ezNmzdTu3Ztbih3G96BBe2OI9chrZrdvn176tSpY0OatK1fv571JyySTx+mf5d2jB0zyu5Icon0anae7eqeGPQ4CxYtYXuEL/53P2x3HLkOAUmeRO1ag1fQTeS/tTbG05uvd0XCrh12R0thFQolqEYrEkpUoduDPfl53RqqVatmdywRl3Hbbbfx0ohXmDBtNr6NX8DDO5/dkeQa3RhSjdPLJuDhfwP5y9bF0zeAdf/CumW5p2ZDIH417sLTvwCTPxtM1cqV6Nmzh92hJIPy7J5XcJzZGFr7diJLN8G/lvaESc6I3LEKr61f8duWTRQuXNjuOOKkPa+5n2VZdLq/G2v+OEFA62dzxZAgyfviT+7n/Ncvs3zJQu644w6744hTejU7Tw/yCA4O5vvvlpCw+Wti/tpsdxxxEwFVmpBY+k5atg0jNjbW7jgiLsMYw5fTp1HcK5KoX2baHUfchM9NpfBvPoh293TkwIEDdseRDMjTzStAmTJlWPjt10Qtf4/4UwftjiNuwv+ObhyMzUeP3g/pCgQimeDr68vyJQvx+HMNUTvX2B1H3ITfrXXwDL2Hu1u0JiIiwu44chV5vnkFuOuuu5jw7ttELhhJUvQ5u+OIGzDGA/8WT7Dyl628/sZIu+OIuJSiRYuyctkSYtd+QtzR3XbHETeRv1YYZwNL06FTF5KSkuyOI+lwi+YVoE+f3vTr2Y2oxWOxkhLsjiNuwMPbF/92wxg7fiJz5sy1O46IS6levTozpk8jYtEYEs+ftDuOuAFjDAFNBrD17xM8+cxzdseRdLhN8wrw1tg3qVvxFqK+n6xDuZIjvAJvJLD9MPr0H8CWLVvsjiPiUsLCwhg+5DkiF44kOT7G7jjiBoynF/5tnmf67K/5aMoUu+PIFbhV8+rh4cHc2V9yY9xRojfNszuOuIl8RcuSv8kjtGjTjqNHj9odR8SlDH1+MG2bNCDqu3exrGS744gb8PQLJDDsRZ59fhirV6+2O46kwa2aV4CAgABWLltM8u+Lid77q91xxE3kr3AHVGpBs1ZtNR2hSCYYY5j6yRTKBBmi1k23O464Ce+CxfFv9QwdOnXmzz//tDuOXMLtmleAm2++maWLFhD9/UTi//nL7jjiJvzr3scJj0J0eeBBkpO1B0kko3x8fFi68FvyHd5E5O8r7Y4jbsKvVCjedbvStEVrzpw5Y3ccScUtm1eAunXr8vGHk4hYMIqkSP1SSvYzxhDQ7DF++m0vw4a/ZHccEZdSqFAhVn63hPifpxN7aLvdccRN+NdoRVThqrTr0ImEBJ3snVu4bfMK0K1bV558dACRi0aTnBBndxxxA8bLm4C2Q5n0yTS++GKG3XFEXEqlSpWYM+tLIpe8RcLZ43bHETfh37APu05EMfDRx3Wydy7h1s0rwOuvvkLj2lWIWjlRv5SSIzz9byAw7AUGPv4Ev/zyi91xRFxKy5YtGfnqCCIXvEFyXJTdccQNGA9PAlo/x9zFK3jv/Yl2xxHUvDqmI/zf5xTzOE/U+tl2xxE34XNTKfI3H0TbsA6ajlAkk5584nE6t2tJ5JK3sJJ1MXnJfh758hMQNpwXX3mNZcuW2R3H7bl98wrg5+fnmI5w7yqidq2zO464ify31sEzNIymLdtoOkKRTJo8aSKVQwKJXPOp3VHETXjfUJSAtoPp0q07O3futDuOW1Pz6hQSEsLyJYuIWTOFuGN77I4jbiJ/rXs4G1CSjvfdr+kIRTLBy8uLRd9+TeCpXUSFLwrfCsIAACAASURBVLE7jrgJ3xJVyNegF81ateHUqVN2x3Fbal5TCQ0N5YtpnzmnI9QvpWQ/Ywz+TR5m875jPP3c83bHEXEpN9xwAyu/W0LCxq+I2R9udxxxE/5VmxJ/8+20ancPcXE62dsOal4v0aFDB4Y99zSRi0aRHB9rdxxxA8bTi4C2Q5j25VdM+fhju+OIuJSyZcsy/5u5RC57h4R/D9sdR9yE/50P8nekB7379tfJ3jZQ85qGF4YOoU2jepqOUHLMhekInxk8lDVr1tgdR8SlNGrUiHfHjSFiwRskxZy3O464AWM88G/xFEvXrGf0m2PtjuN21Lym4cJ0hKUDLaJ+/MLuOOImvG8sgX+rZ7jn3vs0HaFIJg3o359e3ToTuXgsVpIuJi/Zz8PHl4D2wxk17h3mz59vdxy3oub1CvLly8eyRfPxOfQrUdu/tzuOuAm/UqF417mfpi3bcPbsWbvjiLiU8W+P47ZbQ4haNUWHciVHeAUVIrDdULr36kN4uMZd5xQ1r+koVKgQK5ctIe6nz4k9vMPuOOIm/ENbE1WoMu073kdiYqLdcURchqenJ/PmzqZA1EGit2hPmOSMfCHl8Ws0gBZt2nH8uGZ+ywlqXq+icuXKzP7yCyIXj9N0hJJj/Bs9xI6j53nksUF2RxFxKYGBgaxctpik8AVE79todxxxE/6V7iK5XBOat25HTEyM3XHyPDWvGdC6dWteH/ESkQtHkhwXbXcccQPGw5OANs/x1cLvNB2hSCaVLFmSxQu+JWr5+8Sf3G93HHET/vXu52hyEN169NKwlWym5jWDnnpyEPe1aU7kUk1HKDnDI58/AWHDGT7iVb777ju744i4lPr16zN54vtELBhFUpTGj0v2M8bg3+xx1mzawcuvvGp3nDxNzWsGGWOYPGkilYr4E7V2mt1xxE1431CUgDbP0bnrA+zatcvuOCIupUePB3m0X28iF4/BSoy3O464AQ/vfAS0G8b4Dz5i1qzZdsfJs9S8ZoK3tzeLvv2a/P/8RtS2ZXbHETfhe3NVfO7oSbOWbfj333/tjiPiUkaPfJ0G1csRufIDHcqVHOEZUIDAsBfoN/ARNmzYYHecPEnNayYVKFCA779bSsKGWcQc2GZ3HHETAdWaEVeiNq3a3UN8vPYgiWSUh4cHX335BUWS/yVqw1y744ib8ClchvxNH6d1+3s4dOiQ3XHyHDWv16BcuXJ8+/Ucopa+TcLpI3bHETfhf1dP/jpv0affAO1BEsmE/Pnzs2LpIti1nKjdP9kdR9xE/nK341G1Dc1atSUyMtLuOHmKmtdr1LhxY94eO9o5HWGE3XHEDRjjgX/Lp1n8w0+MGTvO7jgiLqV48eJ8t3ghMasmE3dcM9hJzshfuyP/5ivGffc/QHKyppvPKmper8PDAwbQ8/5ORC0Zi5Wki8lL9vPw8SMgbDgj33yLBQsW2B1HxKXcdtttfP7px0QuHE1ihMaPS/YzxuDfdCC//nGAwUOG2h0nz1Dzep3ee+ctapYuQtRqTUcoOcMr6KaU6Qh/++03u+OIuJROnTox+KnHiVw4iuSEWLvjiBswnt4EtB3Cx59/ydSp0+yOkyeoeb1OF6YjvOH830RvWWh3HHET+YpVwPeuvjRv3ZYTJ07YHUfEpbz04nCa31GLqOXvY1k6lCvZzzN/MAFhwxn09LP8+OOPdsdxeSYzewtr165tbdq0KctDBAQE2DKY+eWXX6Zhw4Y0a9bsutdVokQJjh0/jkdAIbwKFKNQ22fwDCjA4Q8fwsPHDzwc2wkFmz+Kb4lKxJ88wOmVH5EUcQqw8K9yN8F3dMUYk7LOf75+naTos4T0ePu680neFPnzlxSP2suvP63F19cXT09PqlWrRkJCAl5eXvTq1YunnnoKD4/Mb6eOGjWKF154IUtynj17li+//JJHH300U4975ZVXCAgI4LnnnrvuDMaYzZZl1b7uFbmQ7KrZ4Pp1OzY2lqDgYJK9/PDwL4Cnf4HrqttJUWf4d+n7JJ4/CclJeAYXpkhnXaheLhbz12bif5jI5g3rKVOmDIDq9hWkV7Pdes/ra6+9liWNK4CXlxcLFywgHwl4Bt3EufVfpdxXpNsoivWZQLE+E/AtUYnkhDhOfvM6wfXuo/iAKYT0mUDckT+I3Lo45THJsZHEn9hHcmwUCWePZ0lGyXv863fjcII/DzinI/Tz8yM8PJwdO3awYsUKlixZwquvXtsX6KhRo9K83bKsTJ94cPbsWSZNmnRNOURSy6q67evrS9EiRSgY5E9wvU74FC17XXX77LoZ+JYKpdhDEynW70MKNOp93Rkl7/ErcxtetTrRtGUbzp8/77hNdTvTclXzunr1aho3bsx9991HxYoV6d69e8o40o0bN3LHHXdQo0YN6tatS0REBLGxsfTp04dq1apRs2ZNVq1aBcC0adPo0KED7du3p3Tp0kycOJF33nmHmjVrUq9ePU6fPg1A7969mTvXcd2/UqVKMWLECGrVqkW1atX4448/ADh58iTNmzenVq1aPPzww5QsWZJTp06lmb9u3bp8OOE9kg9sIuHUwSu+zqida8hXvBJ+pWsB4OHtS8HmAzm3/r9rEEbt/hm/snXxr9SQ6F1rr/OdlbzKGIN/80Gs2vg7r7z2+kX3FS5cmClTpjBx4kQsyyIpKYnBgwdTp04dqlevzkcffQTAsWPHaNiwIaGhoVStWpV169YxdOhQYmJiCA0NpXv37uzfv59KlSrx6KOPUqtWLQ4dOkRAQEDKc82dO5fevXsDcOLECTp27EiNGjWoUaMGP//8M0OHDmXfvn2EhoYyePBgAMaNG5eSZcSIESnrGjlyJBUqVKBZs2bs3r07m99BuV6uXLc9PDyYO3smsWs/wzN/EAlnjl7xdV6tbidFncYzsFDK8j6FS1/vWyt5lH+tdkQUKE/YvZ1JTLz4ZG/V7YzxyvZnyKStW7eyY8cOihUrRoMGDfjpp5+oW7cu999/P7Nnz6ZOnTqcP38ePz8/3nvvPQB+//13/vjjD1q0aMGePXsA2L59O1u3biU2NpayZcvy5ptvsnXrVp5++mmmT5/OU089ddlzFypUiC1btjBp0iTeeustPvnkE1599VXuvvtuhg0bxrJly5gyZUq6+Xv0eJB3x7/Lzn0HsRITADgx8wXw8MB4ehPS8x0STh3Ap2jZix7nXSAEKz6G5LhoPPLlJ3rXGoIbdMPTvwAnvx1NcP0uWfH2Sh7k4Z0P/3bDeGfC4MsKYZkyZUhOTuaff/5h/vz5BAcHs3HjRuLi4mjQoAEtWrTgm2++oWXLlgwfPpykpCSio6O56667mDhxIuHh4QDs37+f3bt3M3Xq1KtuiT/xxBM0atSIefPmkZSURGRkJGPGjGH79u0p61u+fDl79+5lw4YNWJZFWFgYa9euxd/fn1mzZrF161YSExOpVasWt912W/a8cZJlXLluV65cmZlfTKfjfZ3xrdQk5fbM1u3Amm05uWAsEVsW4VcyFP9qzfAKvDEr3l7Jg/wb9eW3+a8z6MlnLrtPdfvqcl3zWrduXUqUKAFAaGgo+/fvJzg4mJCQEOrUqQNAUFAQAD/++CODBg0CoGLFipQsWTKlCDZp0oTAwEACAwMJDg6mffv2AFSrVu2KZ2jfe++9gONyKt98803Kc8ybNw+AVq1aUaBAgStmb9KkScrYlUJFirL5+w8Ax+Enz/zBlyxtLl+BU1LUGRLOHCNfiSoYYzAensSf3I/PTaWu+Bhxb14BBQls/wLHpj3Jpk2bqF37v2FCF/aCLV++nN9++y1lr9W5c+fYu3cvderU4aGHHiIhIYEOHToQGhqa5nOULFmSevXqXTXLDz/8wPTp0wHHWK7g4GDOnDlz0TLLly9n+fLl1KxZE4DIyEj27t1LREQEHTt2JH/+/ACEhYVl8p0QO+SFul29ShX+PLaL5LhoIPN126/MbRR/+BNi/95MzF+bOTbtSYr1/SCNdYiA8fQioM1gvpz9PAmJl19qU3U7fbmuec2XL1/K/3t6epKYmIhlWRedyHRBeiebpV6Ph4dHys8eHh6X7Z269DEXnvdqz3GpVatWUaiQ47BRVFQUtes14Gx81GXLeRe6hbhD2y+6LeHscYyPHx758nP+95Ukx0ZyZHJfR4b4aKJ2rVPzKunyKXIreHrTqm17wjdvpESJEvz11194enpSuHBhLMtiwoQJtGzZ8rLHrl27lsWLF9OjRw8GDx5Mz549L1vG39//op9T/03GxmbukkOWZTFs2DAefvjhi24fP358mn/rkrvlhbptWRY9+/Rl8XfvpLnc1eo2gKdfIP6VG+NfuTH/zH2V2EPb8a/QIMNZxL14+AYQEPYi56cMYOXKlSljuVW3ry5XjXm9kooVK3L06FE2btwIQEREBImJiTRs2JAZM2YAsGfPHg4ePEiFChWy9LnvvPNOvvrKMYh/+fLll22JXIm/vz8rly2G+Fhi9m28+L7KjYk7vJOY/Y5d8ckJcZxZ+RFBt3cCIGrXGop0eZUSj3xGiUc+o2iv8Rr3KhliPDyhSmuatWrL/v37GThwII8//jjGGFq2bMmHH35IQoJjOMuePXuIioriwIEDFC5cmP79+9O3b1+2bNkCgLe3d8qyaSlSpAi7du0iOTk5ZS8XQNOmTfnwww8BSEpK4vz58wQGBhIR8d9MdC1btuSzzz5LOVv9yJEj/PPPPzRs2JB58+YRExNDREQECxfq8nOuytXqtjGGT6dMplwBb6zYy6+icLW6HXNgW8p1Y5Pjokk4cwyvoJuy6iVJHuVdoBjGy4dOXbqye/duTp48qbqdAbluz2tafHx8mD17NoMGDSImJgY/Pz9WrlzJo48+ysCBA6lWrRpeXl5Mmzbtoi33rDBixAi6devG7NmzadSoESEhIQQGBmboscWLF6dIkcKcWzcVn8KlHXvGcIxRvKnTS5xeMZnTKz6E5GT8qzQhsFY7Es+dIOn8SXyKVUxZj/cNRfHIl5+4o7vJVyxri7zkLVZiPOd3rObs+X+oXr0GL730Is8++ywA/fr1Y//+/dSqVQvLsrjpppv49ttvWb16NePGjcPb25uAgICUQ0cDBgygevXq1KpVi5EjR172XGPGjKFdu3bcfPPNVK1aNaWgvffeewwYMIBPP/0UT09PPvzwQ+rXr0+DBg2oWrUqrVu3Zty4cezatYv69esDjssuffHFF9SqVYv777+f0NBQSpYsyV133ZVD75xkNVes2z4+Pixd+C2FCxcheudqAmvfk3JfenUbIP74n5xZMRk8PMFKJrBGC/KFlM/S1yV5k5WUSFSioWq1alQoX57evXvzzDOOsbCq22nLFdd5zc3i4uLw9PTEy8uLX375hUceeSRl8HJGffXVHPo+OojgruPwCiiYTUlF/mMlJnD+m5fp17kt494cbXecXEHXeXUf11u3d+/eTd36DfBr9Sy+t1TPxqQi/4lcM5VyXqdY+8MKvL297Y5jO13n9TocPHiQOnXqUKNGDZ544gk+/vjjTK+jS5fOPPvEY0QtHI2VGJ8NKUUuZry8CWg3lA8mT2HdunV2xxHJUddbtytUqMDc2TOJXPo2CWeOZVNKkYv539WTXcfOMfrNsXZHyfVcYtiAncqVK8fWrVuvez33driH0aPH4B8TocunSI6IPfAbAf5+lCtXzu4oIjkqK+p23bp1CQzwJ+HUQbwLhGRRMpErSzxzlMRTB2nSqKHdUXI97XnNAf/88w8t2rQjqPmjalwlR8Qd20PMmiksX7KIokWLZuqxy5Yto0KFCpQtW5YxY8Zcdv+ZM2fo2LEj1atXp27dumzf7jgDe/fu3YSGhqb8CwoKYvz48VnyekRyUmJiIvfc25mEotXxK1vX7jjiBpJizhOxYCTj3x57TWNG3a1uq3m9xEMPPUThwoWpWrVqlqwvNjaWFm3ak1j6TvwrN86SdYqkJ/H8KSIWjeGLaZ9d8dp/V5KUlMRjjz3G0qVL2blzJzNnzmTnzp0XLTNq1ChCQ0P57bffmD59Ok8++STgONQaHh5OeHg4mzdvJn/+/HTs2DHLXpdIWrK6ZgM8/sTT/HboNP6N++nSbZLtrKQEIhe/SZ/uXejXt2+mH++OdVvN6yV69+7NsmXLsmRdlmXxYK8+HIrzxf+OblmyTpH0JMfHErlwFMOee5oOHTpk+vEbNmygbNmylClTBh8fH7p27cr8+fMvWmbnzp00bdoUcFwOaf/+/Zw4ceKiZb7//ntuvfVWSpYsee0vRiQDsrJmA3wwaRIz5y0goM1gjKdG1kn2siyLqB8+ok65Erwz7trGurpj3VbzeomGDRtSsGDWXBHgtTdGsnJ9OP4tnsAYvdWSvSwrmajl42nd6HZeGDrkmtZx5MgRbr755pSfS5QowZEjRy5apkaNGikzGW3YsIEDBw5w+PDhi5aZNWsW3bppg02yX1bW7BUrVjBk+MsEhL2Ih2/A1R8gcp2iNs+nYMxhvpkzC09Pz2tahzvWbXVU2WTOnLm89d5EAtoNw8Pb1+444gaifppB6YBkpn368TUf6kzr0nmXrmvo0KGcOXOG0NBQJkyYQM2aNfHy+m8PVXx8PAsWLKBz587XlEHEDrt37+a++7sR0PpZvAsUszuOuIHoPzeQ/NtCVi5bTEDAtW8suWPd1jGRbLB582b69B9AYIcROkFLckTU9u/xObieZVs2XdcF30uUKMGhQ4dSfj58+DDFil38RR4UFMTUqVMBR9EsXbo0pUuXTrl/6dKl1KpViyJFilxzDpGcdPr0aZq1bINPvQd0XVfJEfH//E30yomsXLaYW2655brW5Y51W3tes9iRI0do2bY9+Zs8Qr6iZe2OI24g9vAO4n76nJXLllCoUKHrWledOnXYu3cvf//9N/Hx8cyaNYuwsLCLljl79izx8Y7rFX/yySc0bNiQoKCglPtnzpzpMoeeRBISEmgb1pGYkFD8q18+f7xIVkuKOkPEwlF89MEE6tWrd93rc8e6rT2vWSg6OprmrdtBpRbkr3CH3XHEDSSeO0Hk4nHM+fILKleufN3r8/LyYuLEibRs2ZKkpCQeeughqlSpwuTJkwEYOHAgu3btomfPnnh6elK5cmU+/fTTlMdHR0ezYsUKPvroo+vOIpLdLMui74CB7D4dT2C7XnbHETdgJcYTuWgMjw/oS/fuD2TJOt2xbmt62Et069aN1atXc+rUKYoUKcKrr75K3wxcuiI5OZn2HTux/mAUAS2e0OVVJNslx0Vz7qshvDL4KZ5+6gm74+R6mh42b7rWmg3w1tvv8Ma7kwjqPBqPfPmzOam4O8uyiPxuPA3K3MC3c7/Cw0MHv9OTXs3WntdLzJw585oeN2z4S/z8+z6C7n1NjatkOys5ichlb3Nv62Y89eQgu+OI2OZaa/bixYt5ZeQYbrj/TTWukiOifp1DCGeYPWOhGtfrpOY1C/zvf18w6ZNp3NB1LMbL2+444gai1k6j0k1+TPnwA20siWTS9u3b6dq9BwHtXsAruLDdccQNRP3xI2b3SlZs2YSfn5/dcVyemtfr9Msvv/DIoCcJ6vQanv432B1H3EDUtmXkP/Ebi7ZsxNtbG0simfHPP//QrFUbfBs+hG+JSnbHETcQd2wvMas/Yt2q7wkJCbE7Tp6g5vU6HDhwgLZhHfBvPgifm0rZHUfcQMyBbcT/OpOff/2FAgUK2B1HxKXExsbSsm0YiaUaEFC5id1xxA0kRpwiYtFopn/2CTVr1rQ7Tp6hQRfXKCIigqYtWuMZeg9+t9axO464gYTTR4ha+jbffj2HcuXK2R1HxKVYlkXPPn05GJsP/wZZc5a3SHouTNc95Jknuffee+2Ok6eoeb0GSUlJdOjUhbOBpclfK+zqDxC5TkmxkUQsGMlbb46iSRPtMRLJrDdGjmL5T5s1XbfkCMd03e/R8s46vPjCMLvj5Dn6C74GTz07mK1/n8C/yQCdLCPZzkpKJGrxm/To0pGBDz9sdxwRl/P1118z9t0JBLTXdN2SM6J+/pKS+ROYPvVT9QnZQGNeM2nKxx/z+ay5BN8/FuOpt0+yl2VZRK2eQs3SRXj/3bftjiPicrZs2UKvvv2d03Vf3wx0IhkRtWMV3vt/5rvrnK5brkzdVyasXr2aZwYPJbjLaDz9Au2OI24gassigs//zbzv1+Pp6Wl3HBGXcvToUVq2aYdf44c1XbfkiNjDu4j9cSrrf1zLTTfdZHecPEvNawb9+eefdOjUGf9Wz+BdsLjdccQNxOzbRNLWb/h+04aL5qAWkau7MF23VbE5ARXvtDuOuIHEcyeIXDKWr2b8jypVqtgdJ09T85oBZ86coWmL1njX7YpfqVC744gbiD95gKgV7/Pd4gWUKlXK7jgiLiU5OZn7u/fguClIwO2d7Y4jbiA5LpqIBSN59cVhtGnTxu44eZ5O2LqKhIQE2nXoRFThqvjXaGV3HHEDSdHniFg4kknvj6dBgwZ2xxFxOS+8+DI/hu8moNljOllGsl3KdN2t7uaZp5+yO45bUPN6FY88NohdJ6Lwb9jH7ijiBqzEBCIXjWFgn5707NnD7jgiLueLL2bwwcefEdBuKMbLx+444gai1n1OhUK+TJk8SRtLOUTDBtIx/r0JzFm0nOD7x2A8dLKMZC/Lsoj8/gPqVynDm6NH2h1HxOX88ssvDHxsEEH3vY6nv2agk+wX9dt3+B0LZ7Gm685Ral6vYNmyZbz4ymsE3z8Gj3z+dscRNxC18RsKJ55kzqxv8fDQQRGRzLgwXXd+TdctOST2wG/Er/+Sn9b/TMGCBe2O41bUvKZh586ddOnWnYC2Q/C+oajdccQNRO/5BXYsZeWWTfj7a2NJJDMiIiJo1qotHjXCyF+2rt1xxA0knD5C5NK3mP/1V5QvX97uOG5HzeslTp06RbNWbcjXoBe+JXSpC8l+8Sf2Ef3DJFavXE7x4roMm0hmJCUlcW/nrpzOfzMBt91jdxxxAxem6x47+g2aNm1qdxy3pGOTqcTHx9Oq3T3E31wX/6r6hZTslxh5moiFo/hsykfUrl3b7jgiLufp555n059HCLj7YZ0sI9nOSkokask4une6h0cfecTuOG5LzauTZVn07tufvyMM/nfqLG/JfskJcUQuHM0zgx6lSxddi1Iks6Z8/DHTvvyKgLZDMJ46WUayl2O67k+ofsuNTHz/XbvjuDUNG3Aa/eZYlqz+haDOozBGPb1kL8uyiFoxgbtvr84rL79kdxwRl7NmzRrHdN2dR+LppxnoJPtFb11E0Lk/WfD9r5qu22ZqXoH58+czetw7BN8/Fg8fP7vjiBuI+mUmJbyj+HL6NB3qFMmkffv2cc+99zmm677xZrvjiBuI+WsziVu+4fuNv2q67lzA7ZvX8PBwHuzVh4CwF/EKKmR3HHEDUTvX4LlvLcu3bMLX19fuOCIu5ezZs47puut00XTdkiPiTx0kavl7LFs0n9KlS9sdR3Dz5vX48eO0aNMO30b9yVesgt1xxA3EHd1N7LpP+WnNKooUKWJ3HBGXkpiYSPuO9xF5YyUCQjV/vGS/pOhzRC4YycTx73DnnXfaHUec3HZwZ0xMDC3atCe5XGP8KzW0O464gcTzJ4lYNIYZn0+levXqdscRcTmPPD6IHUfP49+4r91RxA1YiQlELn6T/r0eoHfvXnbHkVTcsnm1LIvuPXtzJCkA/3pd7Y4jbiA5PobIhSMZPuQ5wsLC7I4j4nLenzCRrxYsI6DNc5quW7KdZVlE/vAh9SqVYtybY+yOI5dwy+b15VdeZdXG7fg3G6STZSTbWVYyUd+9S9smDRj6/GC744i4nOXLl/PCy68Q0H64puuWHBG18Rtuij/O3NlfarruXMjtxrzOmjWb8R98xA1dx+Hhnc/uOOIGotZN59ZgD6Z+MkUbSyKZ9Mcff3Df/d0IaPM83gVC7I4jbiB6zy9YO5aycvNGTdedS7lV87phwwb6DXyEoI6v4hlQwO444gYif19JviObWLplEz4+PnbHEXEp//77L01btsbnjh743lzV7jjiBuJP/EX0D5NYteI7SpQoYXccuQK3aV4PHz5M63Zh5G/6GD5FytgdR9xA7KHtxP88nR9/+Ykbb7zR7jgiLuXCdN1xxWsTUK253XHEDSRFniFi4Sg+mfwhderUsTuOpMMtBnJERUXRrFVbTNU25C9Xz+444gYSzh4ncslbzJn1JRUrVrQ7johLsSyLPv0GsO+cpem6JUckJ8QRuWgUTz32MF273m93HLmKPN+8Jicn0+n+bpzyCcG/zr12xxE3kBwXReSCNxj56ghatmxpdxwRl/PmuLdYvOonAlo9rSsLSLazLIuolRNoUqcar70ywu44kgF5vnkdPGQYv+46gH/TgTpZRrKdlZxE5JK36BLWiiefeNzuOCIuZ8GCBbwxZhwB7V/QdN2SI6LWz6K4ZyQzNF23y8jTY16nTp3Gx5/PILjrWIynt91x5BodeLPdZbd5+eSuK0UY48GNnV8nfs9aqhQL4sOJE+yOJOJyfvvtN7r37E1g2It4BRW2O45co3/mvkrMvo0X3ebl7QO5rDEs0Lg3lm8wHntXs3zLJvz8tLHkKvJs8/rjjz/y8MBHCGr+CInn/iHx3D92R5JrVOSBMViJCSTHnCPx2G6id/xAaGhNPvloEmXLlrU7HgAdOndl3Y//o7BHNAs3b8DLK8/+aYlkixMnTtCsVRs8yjcG40Hcsb12R5JrFFSvM4F1OpAcE0HivwdJ3LWKYP98fDRpInfffbfd8QCYNGkSL7/3Geb8MX5as4qiRYvaHUkyIc9+wx46cozS5crD/u9hv91p5HokJiYREXGeqIjz5K/chBu7jmHfwW3Uv7Mh77w1lgH9+9sdES9PLzxOH2Dl5o3ccMMNdscRcTmRkZEUvLEQVvQ+2LzP7jhyHZKTLaIiIzj370kCbq1NvuaDLcfOPQAAC+RJREFUiI+PpmvPvnTrfC+TJr5v+wa+j48P0Yd38e28edSoUcPWLJJ5xrKsDC9cu3Zta9OmTdkYR+TKDh48yAcfTuajKZ/gVaQs3rd1Iuq7d3l5yDM89+wztmb7dcNGPAy6vEouZozZbFlWbbtz5CTVbLHT+fPn+Xz6dN4eP5GzMQn4NhlI/MY5NKh8M9/O/QpPT/tOxjt8+DCbN2/hnns0XXdulV7NVvMqLic+Pp6nn3ueL+bOJ3+zQUQsfpNvZv6PFi1a2B1NcjE1ryL2sCyLzz+fzuNPPUNA6+eI2ziHR7q2ZdQbr9sdTXKx9Gp2nr/agOQ9Pj4+fPD+eAb160nM6o/I37AfD/Z+iPPnz9sdTURELmGMoXfvXsyZ+QWRy94mX93OvDdxElu3brU7mrgoNa/isl5/9RXq16hI0vnjJBUqx0dTPrY7koiIXEHr1q0ZMXwYib8txrduF4a+pGuqyrVR8youyxjjKITbl+NdvTXjJ3xAZobBiIhIznr0kYHEHtpBvluqs27NGg4dOmR3JHFBal7Fpd1+++0ULVQAKz6WiKgY9u3TWcoiIrmVv78/D3Z/kLg/1hBQtjYrV660O5K4IDWv4tKMMQzs14fkv3/FN6QsOjlFRCR369+3N8l/byCxQGl+/nWD3XHEBal5FZd388034xEXSWJAYQ4ePGh3HBERSUfx4sWJjzyLV4EQ9u7bb3cccUFqXsXl+fv7YyXGYhkvEhIS7I4jIiLp8Pf3JzE+FuPhqZot10TNq7g8YwwGIHdNmy0iImkw5kKxVtGWa6PmVURERERchppXEREREXEZal5FRERExGWoeRURERERl6HmVURERERchppXEREREXEZal5FRERExGWoeRURERERl6HmVURERERchppXEREREXEZal5FRERExGWoeRURERERl6Hm9f/t3X1MVfcdx/HPQS4KctWwsAhi5sMFy+CiUyqwiIKxhaCiQk2oWtZo0/mQ4Oa0RE34w6S2/5joppVquxo1uZgRBbcISzeMJgZakCKhtglLi+HBuWpCDKZOoWd/bJBeHxrsuJz7k/crIfFefke+hz8+fnI89/wAAABgDMorAAAAjEF5BQAAgDEorwAAADAG5RUAAADGoLwCAADAGJRXAAAAGCPU6QGAH6utrU1tbW1qaWkZeq+1tVUVFRWaNm2aMjMzHZwOAPB9t27d0sWLF3X//v2h97751z9VUVGhsLAwFRQUODgdTEJ5hbEOHX5PH7x/VFEvpCnshSxNiPypLrb+WX/a8Jp+sSBVVz+pd3pEAMD/1F28qHWvvip33BxNSStQ2FSPev8Ro1+9uVWh3z3U0qVLNWXKFKfHhAEorzDWe384pP7+fp2t+bvCp85RyAS3/v3gvjJ+uUh/qT7r9HgAgO95tahId+/e1e/e2i3F/FwhEZP03cSfKGqyW3+rvUBxxbBRXmEsl8ulPx5/X96Dv1fZvlKNC3erIG+Zjh09IpfL5fR4AIBH/PrNN5UQH6/VhWt1f3KMZkdH6q/NTYqOjnZ6NBiE8gqjWZalHb/drhRvkrp7elT82muyLMvpsQAAT5Gdna2mT+p19lyVfrO9ROPHj3d6JBiG8ornwrJly5weAQAwTPHx8Sp9a5fTY8BQPCoLAAAAxqC8AgAAwBiUVwAAABiD8goAAABjUF5hpHPnzsmyLH355ZeSpJ6eHr3yyisOTwUAeBIyGyOJ8goj+Xw+LVq0SBUVFZKk2NhYVVZWPrauv79/tEcDADyCzMZIorzCOH19fbpy5Yo+/PDDoSDs6OhQcnKyJOnEiRNau3atVq5cqZdfftnJUQFgzCOzMdJ4ziuMU1VVpdzcXCUkJCgqKkrNzc2KioryW1NfX6/W1tbH3gcAjC4yGyONK68wjs/nU1FRkSSpqKhIPp/vsTUvvfQSIQgAQYDMxkjjyiuMcufOHdXV1amtrU2WZWlgYECWZWnr1q1+6yZOnOjQhACAQWQ2AoErrzBKZWWliouLdePGDXV0dKizs1MzZ85UV1eX06MBAB5BZiMQKK8wis/n05o1a/zeKyws1P79+x2aCADwNGQ2AsGybXvYi1NTU+2mpqYAjgMAgWFZ1lXbtlOdnmM0kdkATPVDmc2VVwAAABiD8goAAABjUF4BAABgDMorHDFu3DjNmzdv6Ovdd98d9rE/dk/sN954Q9evX3/m4wBgrCOzEUz4wBYcERkZqb6+PqfHwBjCB7aAH4/MxmjjA1swxowZM7Rnzx5lZGQoNTVVzc3NysnJ0ezZs1VeXi7Jf0/szz//XAsXLtS8efOUkpKi9vZ23bt3T8uXL9fcuXOVnJysM2fOSJKysrI0+A+5z+eT1+tVcnKySktLh35+ZGSk9u7dq7lz5yo9PV23bt0a5d8AAJiDzIYTKK9wxLfffuv3X1CDYSVJ06dPV319vTIzM/X666+rsrJSDQ0NKisre+zvKS8v1/bt29XS0qKmpibFxcWptrZWsbGxunbtmtra2pSbm+t3TE9Pj0pLS1VXV6eWlhY1NjaqqqpKknTv3j2lp6fr2rVrWrx4sY4fPx7YXwQAGIDMRjBhe1g4Ijw8XC0tLU/8Xn5+viTJ6/Wqr69PbrdbbrdbEyZMUG9vr9/ajIwMvf322+rq6lJBQYHi4+Pl9Xq1c+dOlZaWasWKFcrMzPQ7prGxUVlZWYqOjpYkrV+/XpcvX9bq1asVFhamFStWSJIWLFigjz/+eKRPHQCMQ2YjmHDlFUFn/PjxkqSQkJChPw++7u/v91u7bt06nT9/XuHh4crJyVFdXZ0SEhJ09epVeb1e7d69W/v27fM75ofu83a5XLIsS9J/P6Dw6M8DAPgjszHaKK8w2ldffaVZs2appKRE+fn5am1tVU9PjyIiIrRhwwbt3LlTzc3NfsekpaXp0qVLun37tgYGBuTz+bRkyRKHzgAAxg4yGyOB2wbgiMH7pwbl5uY+06NXBp05c0anT5+Wy+XS1KlTVVZWpsbGRu3atUshISFyuVw6evSo3zExMTF65513lJ2dLdu2lZeXp1WrVv3f5wQAzysyG8GER2UBGBN4VBYAmINHZQEAAOC5QHkFAACAMSivMEJtba3mzJkjj8fzxPusbNtWSUmJPB6PUlJShm747+zsVHZ2thITE5WUlKRDhw6N9ugAMOaQ2QgkyiuC3sDAgLZt26aamhpdv35dPp/vsf2ua2pq1N7ervb2dh07dkxbtmyRJIWGhurAgQP64osv1NDQoCNHjrBXNgAEEJmNQKO8Iuh9+umn8ng8mjVrlsLCwlRUVKTq6mq/NdXV1SouLpZlWUpPT1dvb69u3rypmJgYzZ8/X5LkdruVmJio7u5uJ04DAMYEMhuBRnlF0Ovu7tb06dOHXsfFxT0WZsNZ09HRoc8++0xpaWmBHRgAxjAyG4FGeUXQe9Lj3AZ3VBnumr6+PhUWFurgwYOaNGnSyA8JAJBEZiPwKK8IenFxcers7Bx63dXVpdjY2GGvefjwoQoLC7V+/XoVFBSMztAAMEaR2Qg0yiuC3osvvqj29nZ9/fXXevDggSoqKpSfn++3Jj8/XydPnpRt22poaNDkyZMVExMj27a1adMmJSYmaseOHQ6dAQCMHWQ2Ao3tYRH0QkNDdfjwYeXk5GhgYEAbN25UUlKSysvLJUmbN29WXl6eLly4II/Ho4iICH300UeSpCtXrujUqVPyer1DWxvu379feXl5jp0PADzPyGwEGtvDAhgT2B4WAMzB9rAAAAB4LlBeAQAAYAzKKwAAAIxBeQUAAIAxKK8AAAAwBuUVAAAAxqC8AgAAwBiUVwAAABiD8goAAABjUF4BAABgDMorAAAAjEF5BQAAgDEorwAAADAG5RUAAADGoLwCAADAGJRXAAAAGIPyCgAAAGNQXgEAAGAMyisAAACMQXkFAACAMSzbtoe/2LK+kXQjcOMAQMD8zLbtaKeHGE1kNgCDPTWzn6m8AgAAAE7itgEAAAAYg/IKAAAAY1BeAQAAYAzKKwAAAIxBeQUAAIAxKK8AAAAwBuUVAAAAxqC8AgAAwBiUVwAAABjjP98PTAovuVlWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ThermalTreatment.plot_sankey_report(margin=.5, gap=.3)"
   ]
  }
 ],
 "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
}
