{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Thermal Reactivation"
   ]
  },
  {
   "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",
    "The thermal reactivation of GAC uses the same thermal treatment model, except it calculates the PFAS remaining in the reactivated GAC instead of an ash residual stream.\n",
    "\n",
    "<img src=\"../Images/ProcessModels/ThemalTreatment.png\" alt=\"Drawing\" style=\"width: 500px;\"/>\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 does not consider transformations of PFAS. \n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Input Parameters for Thermal Reactivation model"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Category</th>\n",
       "      <th>Dictonary_Name</th>\n",
       "      <th>Parameter Name</th>\n",
       "      <th>Parameter Description</th>\n",
       "      <th>amount</th>\n",
       "      <th>unit</th>\n",
       "      <th>Reference</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Total destruction and removal efficiency</td>\n",
       "      <td>DRE</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>DRE</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>DRE</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>DRE</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>DRE</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>DRE of PFHxA</td>\n",
       "      <td>0.90</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>DRE</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>DRE</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>DRE</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>DRE</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>DRE of PFBS</td>\n",
       "      <td>0.98</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>DRE</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>DRE of PFHxS</td>\n",
       "      <td>0.99</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFOA</td>\n",
       "      <td>Fraction of PFOA remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFOS</td>\n",
       "      <td>Fraction of PFOS remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFBA</td>\n",
       "      <td>Fraction of PFBA remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFPeA</td>\n",
       "      <td>Fraction of PFPeA remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFHxA</td>\n",
       "      <td>Fraction of PFHxA remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFHpA</td>\n",
       "      <td>Fraction of PFHpA remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFNA</td>\n",
       "      <td>Fraction of PFNA remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFDA</td>\n",
       "      <td>Fraction of PFDA remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFBS</td>\n",
       "      <td>Fraction of PFBS remained in reactivated GAC</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 GAC</td>\n",
       "      <td>Frac_to_GAC</td>\n",
       "      <td>PFHxS</td>\n",
       "      <td>Fraction of PFHxS remained in reactivated GAC</td>\n",
       "      <td>0.33</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>Reactivation parameters</td>\n",
       "      <td>Reac_param</td>\n",
       "      <td>frac_vs_to_GAC</td>\n",
       "      <td>Fraction of incoming  Volatile Solids to reactivated GAC</td>\n",
       "      <td>0.15</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>Reactivation parameters</td>\n",
       "      <td>Reac_param</td>\n",
       "      <td>GAC_ts_cont</td>\n",
       "      <td>Solids content of reactivated GAC</td>\n",
       "      <td>0.90</td>\n",
       "      <td>fraction</td>\n",
       "      <td>NaN</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                     Category Dictonary_Name  \\\n",
       "0   Total destruction and removal efficiency                   DRE             \n",
       "1   Total destruction and removal efficiency                   DRE             \n",
       "2   Total destruction and removal efficiency                   DRE             \n",
       "3   Total destruction and removal efficiency                   DRE             \n",
       "4   Total destruction and removal efficiency                   DRE             \n",
       "5   Total destruction and removal efficiency                   DRE             \n",
       "6   Total destruction and removal efficiency                   DRE             \n",
       "7   Total destruction and removal efficiency                   DRE             \n",
       "8   Total destruction and removal efficiency                   DRE             \n",
       "9   Total destruction and removal efficiency                   DRE             \n",
       "10  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "11  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "12  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "13  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "14  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "15  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "16  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "17  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "18  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "19  Fraction of not destroyed and removed that remains in GAC  Frac_to_GAC     \n",
       "20  Reactivation parameters                                    Reac_param      \n",
       "21  Reactivation parameters                                    Reac_param      \n",
       "\n",
       "    Parameter Name                                     Parameter Description  \\\n",
       "0   PFOA            DRE of PFOA                                                \n",
       "1   PFOS            DRE of PFOS                                                \n",
       "2   PFBA            DRE of PFBA                                                \n",
       "3   PFPeA           DRE of PFPeA                                               \n",
       "4   PFHxA           DRE of PFHxA                                               \n",
       "5   PFHpA           DRE of PFHpA                                               \n",
       "6   PFNA            DRE of PFNA                                                \n",
       "7   PFDA            DRE of PFDA                                                \n",
       "8   PFBS            DRE of PFBS                                                \n",
       "9   PFHxS           DRE of PFHxS                                               \n",
       "10  PFOA            Fraction of PFOA remained in reactivated GAC               \n",
       "11  PFOS            Fraction of PFOS remained in reactivated GAC               \n",
       "12  PFBA            Fraction of PFBA remained in reactivated GAC               \n",
       "13  PFPeA           Fraction of PFPeA remained in reactivated GAC              \n",
       "14  PFHxA           Fraction of PFHxA remained in reactivated GAC              \n",
       "15  PFHpA           Fraction of PFHpA remained in reactivated GAC              \n",
       "16  PFNA            Fraction of PFNA remained in reactivated GAC               \n",
       "17  PFDA            Fraction of PFDA remained in reactivated GAC               \n",
       "18  PFBS            Fraction of PFBS remained in reactivated GAC               \n",
       "19  PFHxS           Fraction of PFHxS remained in reactivated GAC              \n",
       "20  frac_vs_to_GAC  Fraction of incoming  Volatile Solids to reactivated GAC   \n",
       "21  GAC_ts_cont     Solids content of reactivated GAC                          \n",
       "\n",
       "    amount      unit  Reference  \n",
       "0   0.97    fraction NaN         \n",
       "1   0.97    fraction NaN         \n",
       "2   0.97    fraction NaN         \n",
       "3   0.97    fraction NaN         \n",
       "4   0.90    fraction NaN         \n",
       "5   0.97    fraction NaN         \n",
       "6   0.97    fraction NaN         \n",
       "7   0.97    fraction NaN         \n",
       "8   0.98    fraction NaN         \n",
       "9   0.99    fraction NaN         \n",
       "10  0.33    fraction NaN         \n",
       "11  0.33    fraction NaN         \n",
       "12  0.33    fraction NaN         \n",
       "13  0.33    fraction NaN         \n",
       "14  0.33    fraction NaN         \n",
       "15  0.33    fraction NaN         \n",
       "16  0.33    fraction NaN         \n",
       "17  0.33    fraction NaN         \n",
       "18  0.33    fraction NaN         \n",
       "19  0.33    fraction NaN         \n",
       "20  0.15    fraction NaN         \n",
       "21  0.90    fraction NaN         "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ThermalReactivation = ps.ThermalReactivation()\n",
    "ThermalReactivation.InputData.Data[['Category','Dictonary_Name','Parameter Name', 'Parameter Description', 'amount', 'unit','Reference']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Incoming Spent GAC to Thermal Reactivation "
   ]
  },
  {
   "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>900</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>100</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>VS flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>829.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Carbon flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>630</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>900000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1.5e+06</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   900    \n",
       "2   Moisture flow  kg   100    \n",
       "3   VS flow        kg   829.8  \n",
       "4   Carbon flow    kg   630    \n",
       "5   PFOA           μg   900000 \n",
       "6   PFOS           μg   1.5e+06\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('Spent GAC', 1000)\n",
    "IncominWaste.calc()\n",
    "SpentGAC = IncominWaste.Inc_flow\n",
    "SpentGAC.report()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## PFAS balance in Thermal Reactivation"
   ]
  },
  {
   "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>Reactivated GAC</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  Reactivated GAC  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": [
    "ThermalReactivation.calc(Inc_flow=SpentGAC)\n",
    "ThermalReactivation.report(normalized=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3deVhV1f7H8fdhklHUHFLxhxKODAcZVErRNAIHnE3NnK9zmfdeLc17czZLM83xllOmF/R6UzQEycwh0wSHvCgqj4kJGo6gTMKB9fuDX+f+yBEDNofzfT0Pz9PZZ629v/vk82Gx99rr6JRSCCGEKH8WWhcghBDmSgJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjViVpXLNmTdWwYcMyKkUIISqn48eP31RK1fr99hIFcMOGDYmPjy+9qoQQwgzodLrLD9sulyCEEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaEQCWJQqS0tLfHx88PT0JCwsjPT09FLd//79+/nhhx+Mr1evXs3GjRufaV8bNmzg6tWrJeqTnJyMp6fnMx1PiN+TABalys7OjlOnTpGQkECNGjVYsWJFqe7/9wE8duxYhgwZ8kz7epYAFqI0SQCLMhMYGEhqaqrx9cKFCwkICMDb25sZM2YYt/fs2RM/Pz88PDz47LPPjNtjYmLw9fVFr9fTqVMnkpOTWb16NZ988gk+Pj4cOnSImTNnsmjRIhITE2nVqpWxb3JyMt7e3gDMnj2bgIAAPD09GT16NEoptm3bRnx8PIMGDcLHx4ecnByOHz9O+/bt8fPzIyQkhGvXrgFw/Phx9Ho9gYGBpf4LRZg5pdRT//j5+SkhHsfBwUEppZTBYFB9+/ZV0dHRSiml9uzZo0aNGqUKCwtVQUGB6tq1qzpw4IBSSqlbt24ppZTKzs5WHh4e6ubNm+r69evKxcVF/fzzz8XazJgxQy1cuNB4vP//Wq/Xq4sXLyqllFqwYIGaM2dOsb5KKfXGG2+onTt3KqWUat++vYqLi1NKKZWXl6cCAwPV9evXlVJKRUREqOHDhyullPLy8lL79+9XSik1efJk5eHhUXofmDALQLx6SKbKCFiUqpycHHx8fHjuuee4ffs2wcHBAMTGxhIbG0vLli3x9fXl3LlzJCUlAfDpp5+i1+tp06YNV65cISkpiaNHjxIUFESjRo0AqFGjxhOP/dprr7F161YAtmzZQv/+/QH47rvvaN26NV5eXuzbt48zZ8480Pf8+fMkJCQQHByMj48Pc+fOJSUlhYyMDNLT02nfvj0AgwcP/uMfkhD/p0TfiCHEk/x2DTgjI4Nu3bqxYsUKJk6ciFKKadOmMWbMmGLt9+/fz969ezly5Aj29vZ06NCB3NxclFLodLoSHbt///7069eP3r17o9PpaNy4Mbm5uYwfP574+HgaNGjAzJkzyc3NfaCvUgoPDw+OHDlSbHt6enqJ6xDiackIWJQJZ2dnPv30UxYtWkR+fj4hISGsW7eOzMxMAFJTU7l+/ToZGRlUr14de3t7zp07x9GjR4Gi68cHDhzg0qVLANy+fRsAJycn7t2799BjvvDCC1haWjJnzhzj6Pe3sK1ZsyaZmZls27bN2P7/76tp06bcuHHDGMD5+fmcOXOGatWq4ezszPfffw/A5s2bS/VzEuZNRsCizLRs2RK9Xk9ERASDBw8mMTGRwMBAABwdHdm0aROhoaGsXr0ab29vmjZtSps2bQCoVasWn332Gb1796awsJDatWvzzTffEBYWRt++fYmMjGTZsmUPHLN///5MmTLFGNzVqlVj1KhReHl50bBhQwICAoxthw0bxtixY7Gzs+PIkSNs27aNiRMnkpGRgcFgYNKkSXh4eLB+/XpGjBiBvb09ISEh5fDJCXOhK7o+/HT8/f2VfCmnEEKUjE6nO66U8v/9drkEIYQQGpEAFkIIjUgACyGERiSARYWVlpbG66+/jpubG35+fgQGBrJ9+3bj+2+//Tb169ensLCwWL/o6Gj8/f1p3rw5zZo1Y/LkyeVduhBPRQJYVEhKKXr27ElQUBA///wzx48fJyIigpSUFAAKCwvZvn07DRo04ODBg8Z+CQkJvPnmm2zatInExEQSEhJwc3PT6jSEeCwJYFEh7du3DxsbG8aOHWvc5urqyltvvQUUPd3m6enJuHHjCA8PN7b56KOPmD59Os2aNQPAysqK8ePHl2/xQjwlCWBRIZ05cwZfX99Hvh8eHs7AgQPp1asXX3/9Nfn5+UDRCNjPz6+8yhTiD5EAFiZhwoQJ6PV6AgICyMvLY/fu3fTs2ZOqVavSunVrYmNjtS5RiBKTJ+FEheTh4cG///1v4+sVK1Zw8+ZN/P39iYmJISMjAy8vLwCys7Oxt7ena9eueHh4GJePFKKikxGwqJA6duxIbm4uq1atMm7Lzs4Gii4/rFmzhuTkZJKTk7l06RKxsbFkZ2czZcoU5s+fz4ULF4Cim3WLFy/W5ByEeBIJYFEh6XQ6duzYwYEDB2jUqBGtWrVi6NChzJo1iz179tC1a1djWwcHB9q2bcuuXbvw9vZmyZIlDBw4kObNm+Pp6WlcWF2IikbWghBCiDIma0EIIUQFIwEshBAakQAWmouJiaFp06a4u7uzYMGCB95XSjFx4kTc3d3x9vbmxIkTT+z7008/ERgYiJeXF2FhYdy9exeAvLw8hg8fjpeXF3q9nv3795f5+QnxSA/7orhH/ciXcorSZjAYlJubm7p48aK6f/++8vb2VmfOnCnWJioqSoWGhqrCwkJ15MgR1apVqyf29ff3N36R5tq1a9Xf/vY3pZRSy5cvV8OGDVNKKZWWlqZ8fX1VQUFBeZ2uMFPIl3KKiujYsWO4u7vj5uaGjY0NAwYMIDIyslibyMhIhgwZgk6no02bNqSnp3Pt2rXH9j1//jxBQUEABAcHG+cUnz17lk6dOgFQu3ZtqlWrhtxYFlqRABaaSk1NpUGDBsbXLi4upKamPlWbx/X19PRk586dAPzrX//iypUrAOj1eiIjIzEYDFy6dInjx48b3xOivEkAC02ph0yD/P23ED+qzeP6rlu3jhUrVuDn58e9e/ewsbEBYMSIEbi4uODv78+kSZN48cUXsbKSB0KFNuRfntCUi4tLsRFoSkoK9erVe6o2eXl5j+zbrFkz4/oQFy5cICoqCihaHe2TTz4x9nnxxRdp3Lhx6Z+YEE9BRsBCUwEBASQlJXHp0iXy8vKIiIige/fuxdp0796djRs3opTi6NGjODs7U7du3cf2vX79OlD0KPLcuXONy1pmZ2eTlZUFwDfffIOVlRUtWrQoxzMW4r9kBCw0ZWVlxfLlywkJCaGgoIARI0bg4eHB6tWrARg7dixdunRh9+7duLu7Y29vz/r16x/bF4rWi1ixYgUAvXv3Zvjw4UBRMIeEhGBhYUH9+vX58ssvNThrIYrIo8hCCFHG5FFkIYSoYCSAhcm6efMmI0eP4dixY1qXIsQzkQAWJuns2bN4+/rz7yNJdAwOJTw8QuuShCgxCWBhcvLz83mxbTtu3rxFYV422Njz+usDOXLkiNalCVEiEsDC5FhbW/P1zkgWLZjHvQs/0qdLJ6Kjo2nZsqXWpQlRIjINTZiktm3bYmtrC0CTJk0IDQ3VuCIhSk5GwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFycrJyQGgSpUqGlcixLORABYmKzExEYA6depoXIkQz0YCWJisv8+YBUBYWJjGlQjxbCSAhUk6ePAg13+9ysDBQ6lWrZrW5QjxTCSAhck5f/487du3B+C9dyZrXI0Qz04CWJiMjIwMln76Kc2aNQMgPDwCT0/PMj/uin98ztJly8v8OML8yLcim5CZs+fy669ppbrPf6wq32AZM+7NEvfJN+STeu1X9nwdadwWGRlJ9+7dS7O0R4qKiibm60gauf5PuR1TmAedUuqpG/v7+6v4+PgyLEc8isFgoEkLL65mKeybB5Xafu98+3mp7etJdNa2VAsaXOJ+v6/RsWo1qtjZl1ZZT3Q3/TYOrV/D8NMuDh/cj7e3d7kdW1QOOp3uuFLK//fbZQRsIqysrPh2z258A1pjXcMFOze/UtlvVf8epbKfsmTfuA0qPw+djTbzfW11Flg51STLuQ7Bnbty+kS8TH0TpUIC2IQ0atSIXTu+IrRbDyz7zsWm5v9oXVK5sHKuGGHn0KI9melXebVLGD8ePoitra3WJQkTJzfhTEzbtm1ZvmQxmTvnUZCdoXU5ZschsD8p+fYMGjKMkly+E+JhJIBN0LBhQxk19HUyv16AMuRrXY5Z0ekscAieyLc/nmbm7DlalyNMnASwiVr44QLatGhE5r5VMhIrZxbWVXAMm8biZSvZsmWr1uUIEyYBbKIsLCzYtuWf1Mr7lay4r7Qux+xYOdbAKew9Ro4Zi8wMEs9KAtiEOTg4sDcmCnUmmuwLR7Qux+zY1HkB+47jCe0aRkpKitblCBMkAWziXFxciPl6J9n7VpKX9rPW5Zgd+yaB4BFKcOduZGVlaV2OMDESwJVAQEAAa1av4t6u+RRk3tG6HLPjENCH61a16DdgEIWFhVqXI0yIBHAlMWBAfyZNGEPm1/MpzL+vdTlmRafT4dhpAkfO/My706ZrXY4wIRLAlcjsmTN4OcCLrL3LZGZEOdNZWePYbSr/WL+RL77YqHU5wkRIAFciOp2OzRs3UN8yk6yjEVqXY3Ys7Z1xDJvOhLf/zOHDh7UuR5gACeBKxs7Ojtjdu7BI2k9W4kGtyzE7NrVccQieSLcevUhOTta6HFHBSQBXQs8//zzfREeRc+Bz7l89r3U5ZsfuBX8sW/aiU0gX7t69q3U5ogKTAK6k9Ho9m79Yz72vF2C4e0PrcsyOg28YGVUb0atvfwoKCrQuR1RQEsCVWI8ePXjvnb+SuWs+hXk5WpdjVnQ6HQ4dRnMy+ToT//xXrcsRFZQEcCU37d136NIhkKw9S1BK5qiWJ52lFQ5d3uHLrdtZ/Y9/aF2OqIAkgCs5nU7HhrWf08hJkfX9l1qXY3YsbR1x6j6dye++x3fffad1OaKCkQA2AzY2NsR8HYnNlWNk/mev1uWYHesa9XHo/Fd69ulHUlKS1uWICkQC2EzUrFmTvTG7yfthI7kpZ7Qux+zYueqxaT2QTq925s4deVxcFJEANiMtWrRga/hmMqMWkp/+q9blmB0HfSjZdbzp1qM3+fmykL6QADY7oaGhzJ35Ppk751F4X1bvKm8OQcNIvJ7NmHET5HFxIQFsjt6e+CZ9uwaTGf0xqlDmqJYnnYUljp0n8+/ovXyy5FOtyxEakwAuRY6Ojpoc9/3332fv3qe/uabT6Vi9cjnN6ziQdWBdsfdSVo3g6toJXF33Jmlb/m5c3tK4ff1bXF3/FrkpicY+d+N2cHlRr2Ij6sL8XG7sWljUZ+14ft30jsxF/j8WVey5ey2ZKZP/iqurK3q9nsWLFz/zUpbz588vtdrS09NZuXJlifvNnDmTRYsWlVod5kICuBKYPXs2r7zySon6WFtbExX5FQ43zpB1KrrYe3UGzqfeiOXYPO9OxtGtxbcPX0a94cuwdWlu3J6VeIAqdRsX+1aOe/E7sbSvRr2RK6g3ciXPdZ6IzsLqGc+w8tFZVaHWwA+4czeT5cuXs3v3bmbNmvVM+3pUACulShzqzxrA4tlIAJeB/fv306FDB/r27UuzZs0YNGiQ8XpfXFwcL774Inq9nlatWnHv3j1yc3MZPnw4Xl5etGzZ0jhfdMOGDfTs2ZOwsDAaNWrE8uXLWbx4MS1btqRNmzbcvn0bgGHDhrFt2zYAGjZsyIwZM/D19cXLy4tz584BcOPGDYKDg/H19WXMmDG4urpiMBj4Njaa/Lgt5CSfeuA8bBt4kn/n6mPPNf/ONVReLtXaDSbr7AHj9oKsO1g5PWd8bf2cCzor6z/wqVY+ti4tqPLSUPoPGswHH3zA8uXLUUpRUFDAlClTCAgIwNvbm3/830Mc165dIygoCB8fHzw9PTl06BBTp04lJycHHx8fBg0aRHJyMs2bN2f8+PH4+vpy5cqVYn+Zbdu2jWHDhgGQlpZGr1690Ov16PV6fvjhB6ZOncrFixfx8fFhypQpACxcuNBYy4wZM4z7mjdvHk2bNuWVV17h/HlZc+RZyJCkjJw8eZIzZ85Qr149XnrpJQ4fPkyrVq3o378/W7ZsISAggLt372JnZ8fSpUsB+M9//sO5c+d49dVXuXDhAgAJCQmcPHmS3Nxc3N3d+fDDDzl58iR//vOf2bhxI5MmTXrg2DVr1uTEiROsXLmSRYsWsWbNGmbNmkXHjh2ZNm0aMTExfPbZZwC4u7uz49//olvPPlj1Kz6Syrl4DJtaDY2v08LfAwsLdJbW1B2yGCga/do3D6JKAw/yo1IpyErH0qEajl7BpG39O1nnD2Pr6oOjZ0esa9Qvi4/apDl4diLzTipj33ybwsJCrl+/TmRkJM7OzsTFxXH//n1eeuklXn31Vb766itCQkKYPn06BQUFZGdn065dO5YvX86pU0W/QJOTkzl//jzr169/4kh24sSJtG/fnu3bt1NQUEBmZiYLFiwgISHBuL/Y2FiSkpI4duwYSim6d+/OwYMHcXBwICIigpMnT2IwGPD19cXPz6/MP6/KRgK4jLRq1QoXFxcAfHx8SE5OxtnZmbp16xIQEABA1apVAfj+++956623AGjWrBmurq7GAH755ZdxcnLCyckJZ2dnwsLCAPDy8uL06dMPPXbv3r0B8PPz46uvvjIeY/v27UDRTIjq1asb23fo0IHFCxcw+e9zQBUag9amVkNqtBtsbFdn4Hws7Z2LHSs78SC1ek1Hp7PAvkkg2ee/x8m3GzZ13Kg/Zg25l06Sc/kU1zb+hbpvLMK6ZoM/8KlWTg5t3yA56iOysrJQShEbG8vp06eNf9VkZGSQlJREQEAAI0aMID8/n549e+Lj4/PQ/bm6utKmTZsnHnffvn1s3Fi0eLylpSXOzs4PzFGOjY0lNjaWli1bApCZmUlSUhL37t2jV69e2NvbA9C9e/dnPn9zJgFcRqpUqWL8b0tLSwwGA0opdDrdA20fNx3p/+/HwsLC+NrCwgKDwfDYPr8d90nHABg9ahSnE86ycuVKnh+8GCunGo9tD5B3/RL5d66StuXvRRsKDVg5P4+Tb7eiGm3ssG/6IvZNX0Sn05Hzc5wE8EPodBZUaTWQO4k/sGbdepRSLFu2jJCQkAfaHjx4kKioKAYPHsyUKVMYMmTIA20cHBx+t////pvLzc0tUW1KKaZNm8aYMWOKbV+yZMlD/y2LkpFrwOWoWbNmXL16lbi4OADu3buHwWAgKCiIzZs3A3DhwgV++eUXmjZtWqrHbtu2LVu3Ft1Qi42NfejTWEsXL6KKlSXZh794qjmqWYkHqPbS67iMW1f0M2EjhsxbGDKuk5tyloLcTABUQT55N69gWbV2qZ5TZVGQncGd79bg2DKMBR8voXbt2qxatcr4sMaFCxfIysri8uXL1K5dm1GjRjFy5EhOnDgBFN1QfdyDHXXq1CExMZHCwkLjX0EAnTp1YtWqVUU1FBRw9+5dnJycuHfvnrFNSEgI69atIzOz6P9lamoq169fJygoiO3bt5OTk8O9e/fYtWtXqX8u5kBGwOXIxsaGLVu28NZbb5GTk4OdnR179+5l/PjxjB07Fi8vL6ysrNiwYUOxkW9pmDFjBgMHDmTLli20b9+eunXr4uTkVKyNpaUltWvXwiLzCpknduLg1+Ox+8xKPESdfjOLbbNv3IasxINYOlbnduxKUAqUwu4Ff+ybvlSq52TKlCGPq+vfgoICsLDAwaMjVVv1JO/X9mz852wG9O2Nr68vSilq1arFjh072L9/PwsXLsTa2hpHR0fj5YPRo0fj7e2Nr68v8+bNe+BYCxYsoFu3bjRo0ABPT09jmC5dupTRo0ezdu1aLC0tWbVqFYGBgbz00kt4enrSuXNnFi5cSGJiIoGBgUDRVMtNmzbh6+tL//798fHxwdXVlXbt2pXfh1eJ6EryNI6/v7+Kj48vw3JEWbl//z6WlpZYWVlx5MgRxo0bZ7zR8nuXL1+mpX8rrDuMw+6FgHKuVGQlHkJ37Et+OhFP3bp1tS5HlAKdTndcKeX/++0yAjYTv/zyC6+99hqFhYXY2Njw+eefP7Ktq6srUTt38EpoVyz7zik2E0KUPYfm7chMT+XVLmGciv8RS0tLrUsSZUQC2Ew0btyYkydPPnX7wMBAevcMY/epKGyCJ5RhZeJhdFXrcOPiAVkvopKTm3DiobZv3872r2Owa9Nf61LMTm5KIrnfr+fbPdFYWT16jBQTE0PTpk1xd3dnwYIFD7x/584devXqhbe3N61atSIhIQGA8+fP4+PjY/ypWrUqS5YsKbPzEY8mAVwJjRgxgtq1a+Pp6flM/U+ePMng4SNx6jYNK6eapVydeBxDxnUyd3/Els1f4uHh8ch2BQUFTJgwgejoaM6ePUt4eDhnz54t1mb+/Pn4+Phw+vRpNm7cyNtvvw1A06ZNOXXqFKdOneL48ePY29vTq1evMj0v8XASwJXQsGHDiImJeaa+165d49Uu3bDrMIYqdRuXcmXicQrvZ3Nv51xmTp9Kly5dHtv22LFjuLu74+bmho2NDQMGDCAyMrJYm7Nnz9KpUyegaApkcnIyaWlpxdp8++23vPDCC7i6upbuyYinIgFcCQUFBVGjxpMfpPi9nJwcgjt3QzXthEOztmVQmXgUVVhAZszH9Ap5mb/+5c9PbJ+amkqDBv99qMXFxYXU1NRibfR6vfFJyGPHjnH58mVSUlKKtYmIiGDgwIGlcAbiWUgAC6Doiaf+gwZzTVXDofVrWpdjdrIOfUHTmrZ8/o9VT/WE2cNuzv2+39SpU7lz5w4+Pj4sW7aMli1bFrumnJeXx86dO+nXr98fPwHxTGQWhADgb+/P4NDJc1TtM0ceMS1nWaf3YHftFFEn4rC2froV41xcXLhy5YrxdUpKCvXq1SvWpmrVqqxfvx4oCuxGjRrRqFEj4/vR0dH4+vpSp06dUjgL8SxkBCz45z/D+XT1Ghy7TUVnZaN1OWYl95fT5B39J9/GRpfoslFAQABJSUlcunSJvLw8IiIiHlgQJz09nby8PADWrFlDUFCQcQEogPDwcLn8oDEZAZu5H3/8kdHj36Rq71lYOlR/cgdRavLvXCUz+mN2/CuCJk2alKivlZUVy5cvJyQkhIKCAkaMGIGHhwerV68GYOzYsSQmJjJkyBAsLS1p0aIFa9euNfbPzs7mm2++Ma41LLQhjyJXQgMHDmT//v3cvHmTOnXqMGvWLEaOHPlAu19++YWW/q2wDBqNvXtrDSo1XwW5mdzd8g4fvD+V8ePGaV2OKGPyKLIZCQ8Pf2KbzMxMXgntioVXNwnfcqYKDGTtXsjrvXtI+Jo5uQZshgoLC+nz2kBu2dXH3r+n1uWYFaUUWfvX4P0/z7H800+0LkdoTEbAZuivU94h7sIVnHrNlBkP5Sz7ZBRV05PY+e0xWWRHSACbm7Xr1rF20xacB3yEzlK+JLM85fx8HMOJf/Nt3I/FZiMI8yUBbEYOHjzIxL9MoVq/eVjaSQCUp7ybv5AVu5SYryOLzcUV5k2uAZuJixcv0r1XHxxDJmH9nHwvW3kqyM4gc+c8li9ZTNu28oi3+C8JYDOQkZHBKyFdsPbvh10jX63LMSvKkE9m1IeMGvo6w4YN1bocUcFIAFdyBoOBsF59uVujKQ4tu2pdjllRSpG1bzVtmjdk4YcPrtcrhFwDruQmTHybhNR0nHq8pXUpZic7fjvP3b/Kti1HsLCQsY54kARwJbZ8xQoitkfh3P9DdBYy5ak8ZSf9SGHCbr49HoeDg4PW5YgKSgK4koqNjWXq32bg/NoHWNg6al2OWclL+5nsb5fz3Td7cHFx0bocUYHJ30WV0Llz5+jb/3UcO0/Gunq9J3cQpaYg8w73ds1nzepVBAQEaF2OqOAkgCuZW7du8UpoF6oEDsL2f7y0LsesFObfJ/Pr+UyaMIYBA+TLTMWTSQBXInl5eXTp3ovcer44eL+qdTlmRSlF1t5ldPD3ZPbMGVqXI0yEXAOuJJRSjBw9lqTb+TiFDdG6HLOTdXQL9S0z+eeXu2V9DfHUJIAriYWLPmbX3kNU7TdfZjyUs6zEQ1gkfUfsiXjs7Oy0LkeYEAngSmDXrl3M/uBDqvX/CIsq9lqXY1buXz1PzoHPOHzgO55//nmtyxEmRq4Bm7jTp0/z+uChOHWdipVzba3LMSuGuzfIjPqQzV+sR6/Xa12OMEESwCYsLS2N4M5dsQ0aSZX6zbQux6wU5uWSuWs+06b8hR49emhdjjBREsAmKjc3l5Cu3Slwa4dDiw5al2NWlCoka88ndOkQyLR339G6HGHCJIBNkFKKN4YO58p9WxxelK8VL29Z32+ikZNiw9rPZcaD+EPkJpwJmjVnLnuPnqJq37nodPI7tDxlJXyLzZUfiTkRj42NjdblCBMnAWxitm79F4uWrqDagI+wsLbVuhyzkptyhvuHv+DA4UPUrFlT63JEJSABbELi4+MZMXoMTj1nYOX0nNblPLWUlcMouHezTI9hY1uG8291llQNHk/OwXVsC99MixYtyu5YwqxIAJsIg8FAl7AeFNjVID/xO/ITvyuV/d46Flkq+3lazh5BJe5jyMnEkHcfLCwoSLuI4X42s+fMYeCAAQBUqVKF6tWrl3apRr1eG8i+qMUsWvQxoaGhZXYcYX4kgE2ElZUVSz5eyI0bN0p1v5PKMYAdHByYNap3ifsZDAZ+TbvOpvAt3K/TAEOtJny0ZAUWFlZMf29qGVRanL+fL24NXZn0tixqL0qXTin11I39/f1VfHx8GZYjxKMVFhayZ88e/jR2PNk1m1Nw5TQfz5/JqD/9SevShHgsnU53XCnl//vtcgtdmIWyfosAAAhkSURBVAwLCws6d+7MybgfqXonCYuGAUz681+5e/eu1qUJ8UwkgIXJqV27NlE7t2O4cAjbBh58+eWXWpckxDORABYmydPTE09PD/JsnFm0dBkluZQmREUhASxM1lvjRmNruEfar2n8+uuvWpcjRIlJAAuT1bBhQ8jJoIqjM3fu3NG6HCFKTAJYmKwqVaqgCg1YWNmQn5+vdTlClJgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWBhkjZv3szWrVuNryMiIggPD9ewIiFKTr6WXpic/Px8XN3cuX7zFjVDJoAq4Oae1dR9vg7JFy9gaWmpdYlCFCNfSy8qDWtra07GH8OjRQvyLx4l/+KP6PXenIg7KuErTIoEsDBJderU4cfDB+nq35jubVrww8HvqFWrltZlCVEiEsDCpGzfvh2dTse5c+ewtbVlwbzZ3Eu/RZUqVbQuTYgSkwAWJiU8PJy2bdsSEREBQL169di2bdsD7QwGQ3mXJkSJSQALk5GZmcnhw4dZu3atMYCTk5Px9PQEYMOGDfTr14+wsDBeffVVLUsV4qlYaV2AEE9rx44dhIaG0qRJE2rUqMGJEyeoUaNGsTZHjhzh9OnTD2wXoiKSEbAwGeHh4QwYMACAAQMGPHTeb3BwsISvMBkyAhYm4datW+zbt4+EhAR0Oh0FBQXodDrGjx9frJ2Dg4NGFQpRcjICFiZh27ZtDBkyhMuXL5OcnMyVK1do1KgRKSkpWpcmxDOTABYmITw8nF69ehXb1qdPH+bPn69RRUL8cfIoshBClDF5FFkIISoYCWAhhNCIBLAQQmhEAliUK0tLS3x8fIw/CxYseOq+V69epW/fviU+5p/+9CfOnj1b4n5ClDW5CSfKlaOjI5mZmVqXIUS5kptwokJr2LAh7733HoGBgfj7+3PixAlCQkJ44YUXWL16NVB83YczZ87QqlUrfHx88Pb2JikpiaysLLp27Yper8fT05MtW7YA0KFDB34bOISHh+Pl5YWnpyfvvvuu8fiOjo5Mnz4dvV5PmzZtSEtLK+dPQJgjCWBRrnJycopdgvgtJAEaNGjAkSNHaNeuHcOGDWPbtm0cPXqU999//4H9rF69mrfffptTp04RHx+Pi4sLMTEx1KtXj59++omEhARCQ0OL9bl69Srvvvsu+/bt49SpU8TFxbFjxw4AsrKyaNOmDT/99BNBQUF8/vnnZftBCIE8iizKmZ2dHadOnXroe927dwfAy8uLzMxMnJyccHJywtbWlvT09GJtAwMDmTdvHikpKfTu3ZvGjRvj5eXF5MmTeffdd+nWrRvt2rUr1icuLo4OHToYF24fNGgQBw8epGfPntjY2NCtWzcA/Pz8+Oabb0r71IV4gIyARYXx26LqFhYWxRZYt7CweGB939dff52dO3diZ2dHSEgI+/bto0mTJhw/fhwvLy+mTZvG7Nmzi/V53P0Oa2trdDodUHSjUNYTFuVBAliYpJ9//hk3NzcmTpxI9+7dOX36NFevXsXe3p433niDyZMnc+LEiWJ9WrduzYEDB7h58yYFBQWEh4fTvn17jc5ACLkEIcrZb9eAfxMaGlqiqWi/2bJlC5s2bcLa2prnn3+e999/n7i4OKZMmYKFhQXW1tasWrWqWJ+6devywQcf8PLLL6OUokuXLvTo0eMPn5MQz0qmoQkhRBmTaWhCCFHBSAALIYRGJICFEEIjEsBCczExMTRt2hR3d/eH3pBTSjFx4kTc3d3x9vY2zm64cuUKL7/8Ms2bN8fDw4OlS5ca+9y+fZvg4GAaN25McHAwd+7cAYq+2ujll1/G0dGRN998s3xOUIhHkAAWmiooKGDChAlER0dz9uxZwsPDH1g4Jzo6mqSkJJKSkvjss88YN24cAFZWVnz88cckJiZy9OhRVqxYYey7YMECOnXqRFJSEp06dTIGu62tLXPmzGHRokXle6JCPIQEsNDUsWPHcHd3x83NDRsbGwYMGEBkZGSxNpGRkQwZMgSdTkebNm1IT0/n2rVr1K1bF19fXwCcnJxo3rw5qampxj5Dhw4FYOjQocZHjh0cHGjbti22trbleJZCPJwEsNBUamoqDRo0ML52cXExhmhJ2iQnJ3Py5Elat24NQFpaGnXr1gWK5v9ev369rE5BiGcmASw09bB56L89Evy0bTIzM+nTpw9LliyhatWqpV+kEGVEAlhoysXFhStXrhhfp6SkUK9evaduk5+fT58+fRg0aBC9e/c2tqlTpw7Xrl0D4Nq1a9SuXbssT0OIZyIBLDQVEBBAUlISly5dIi8vj4iICOOqaL/p3r07GzduRCnF0aNHcXZ2pm7duiilGDlyJM2bN+cvf/nLA32++OILAL744gt55FhUTEqpp/7x8/NTQpS2qKgo1bhxY+Xm5qbmzp2rlFJq1apVatWqVUoppQoLC9X48eOVm5ub8vT0VHFxcUoppQ4dOqQA5eXlpfR6vdLr9SoqKkoppdTNmzdVx44dlbu7u+rYsaO6deuW8Xiurq6qevXqysHBQdWvX1+dOXOmnM9YmBsgXj0kU2UtCCGEKGOyFoQQQlQwEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGinRcpQ6ne4GcLnsyhFCiErJVSlV6/cbSxTAQgghSo9cghBCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjfwvAjTqz+piJ+4AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ThermalReactivation.plot_sankey()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdeZyN5f/H8dc1G7MwSHZfEdkZa6Qs2dcoa0UkQtFG0SYViRZFSJv8KpSSfS2kInuyJCm7LNlmX6/fH+eYxj7DzNxz5ryfj4dHzTn38j7nzHzO514vY61FRERERMQT+DgdQEREREQktdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPOazowxOY0x1hhTzOksqWGMucUYc9rpHCIiTlHdFvEsXtG8GmMiUvxLMsZEp/j5vqvM28IY82c6ZllrjLk/vZZ3vay1f1hr82TEso0xM4wxse73+aQxZrExpoz7udHGmPgLPptBKebtYIzZYIyJNMacMMZMM8YUvsQ6Wri/dAZd+Fwq8g0wxvxljAk3xsw3xhS8wrQ3GmPmufP8bYzpeMHzLY0xf7ifX57yS/CC9+Hcv3XXkLeWMWaLMSbKGLPOGFMpxXP9jDEJF6yjblrXIZJVqG5fnqfWbfdGwrvGmEPuef8yxryehmyvG2O2GWMSjTFDrzLtpbIWSfH8zcaY1e56ut0YUz/Fc5eqpxHGmBtSm9W9nKt9b4QYY6YYY/41xpwxxixLy/K9mVc0r9bakHP/gP1A2xSPfe50vmzuFff7/j8gHPggxXOfpvxsrLXvAri/mD4BxgA3AFUAX+AHY0zuC5b/AHDS/d9UM8Y0B54HWgL5gaPAtCvMMgU4BRQAegMfpyjohYGZwBD3snYAn10w/ysXvNbal8m11hhT5xKPBwJz3DnyAl8Bs40xfikmW3nBOtZc+V0QybpUtx2VUXV7OFAeqA7kApoAW9OQaxfwFJDaJu/CrIdTPDcLWA3kA14FvjXGpNwguLCehlhr/71wBe4m+XKN9GW/N9ymAjmAW9w5rtiQy3+8onm9GmNMoDHmPWPMEWPMQWPMWGOMv3srazZQKuWWlzGmnjHmF/eW0mFjzNsXNBGpXW8LY8yfxphnjTHH3Vuj96V4Pti9lXrAva5V59ZjjLnHGLPDGHPauPb0lUkx3z/GmCfdW5MRxphJxpjCxphlxpiz7i3p3O5pyxljElLMu9YYM9z937PGmIXGmLwpnn/IGLPfnfdp97puv9prtdZGADOASleazv363gBetNZ+aa2NcRecHu5JHk0xbShwF9APqGpS7IlMhbbAdGvtLmttLDASaGaMKXqJTHnd079grY201n4PLAHOfVadgA3W2jnW2mjgReA2Y8xNachzNU2BGGvtRHfeN3EV/6u+9xe68DN3P5a8Z8kY4+f+vfvXGLPHGDPowulFnKa67ZF1uxbwtbX2qHX5Ky0bItbaj621S4CI1M5zmbxVcDWMr7izTgf2AO2vZ7kXrOOK3xvuDI2BAdbaf621idbajVdY3nmfmXE1zR+m+PmaPmNPpebVZQSurcTKQA2gIfC0eyurA/DXBVte8bj+GPMBd+D6BX3oGtddAjBAEfcyJxtjQtzPvQuUw/UHnw/XnkLrbtKmAgNwbdGtAuZeUIg7AA2ACkBXXHvtngQKAiFA/ytkuhfXH1hhIA/wGIAxJgx4C+gMFHP/y5+aF+kuut2AzVeZtBJQCNeexWTW2kTgG1xN3DmdgeO4tqBX8F+hPLfOXcaYuy8Xyf0v5c/n1n+hckC4tXZfisd+BSq6/7+i++dzWU/j2lNUkfRz4TqSgG0XrKOOcR2q22WMGWaMuda/70dx/e5UAmoDHa88uYgjVLfP5wl1ey3wjHEdlr+oPhpjPjLGvJWabKnU0bhOffjNGJPys64I/OHe2XBOypqeHq72vVEHV8P8urtu/2qMaXctK7qez9hTqXl1uQ8Ybq09Ya09iusQQvfLTWytXWetXe/eUtoDfIir4FyLKOA1a228tXY2YIHSxhh/XM3YQGvtP+51rXYXg27AbGvtSmttHDAK1y9qzRTLHed+PfuBn4GfrLW/uf9Y5wDVrpDpA2vtHmttJK7GMMz9eGdcW81r3Xv/nufqv0PPGdeFBbtwHUJKWUC6u/dAnPuXz/06LK7D+Bc6wvl/kA/g2ntqgS+A+40xvueetNaWtdZ+c5lci4B7jTEVjTFBwAvu9QZdYtoQ4MwFj53BteczNc+D+31I8e/9y+S6nKutYzmuL5ACQBegF+4vr2vQGXjLWnvE/aU/5hqXI5KRVLfP5wl1ewQwDlft3uTeY97t3ITW2t7W2ievki21PsPVQN4IPAK8Zozp4H4uNTW7wQWvc3sa13+1dRTDtdF1CNdG0GBgujGmdBrXA9f2GXu0bP3iUsMYY3BtMabcOtoHXHT4OMU8FYwxi4wxR40xZ3EdJr7WrZzj7r1o50Th+qUvDPgBf11iniIp87oL46ELMqcsItGX+DmEy/vnEnnOrfdAivWe5eI/zguNtNbmsdYWttZ2uGAr9P/cz537dxI4gWuPxqUunirsfh5jzM1APeDcIaevcZ0L2vQS813EWjsfeB2Yi+s9/g2IAw5eYvII4MJzbXPjOhcsNc/Df+/DuX8Pu19HjpQFEtcX2bIUjz2emnVYa/+01u611iZZa7fg+mK81j2m533OF/y/iONUty8py9dtd7P/jrW2Lq56/RYwzV3P05W1dluKDYgfgPf4ryampmavuuB1Ju+VdZ/Kca5mPw4MT1GzZ6VyHdG4PqfXrbVx1tpluDZYGl/Dy72Wz9ijeX3z6t5r9w+uw0Dn/A9XUQHX1uSFPgA2ATdba3MDL3P+Iej0cARIAEpd4rnDpMjr3ttYlP8yZ5QjuLYWz603NxCazuvYhqtgd0r5oPs1dgC+cz907gKtZcaYf4A/cH1pnHfqwJVYa9+21t5srS2E6wKAOOD3S0z6O5DbGPO/FI9VBc5tiW93/3wuayiuz+eqW+rW2tiUBRLYADRN8di4y6zDB9ee1sutw3Ll30kf916ic1JeqHDe5wwUv9rrEMlMqttpzpRV6nYya22UtfYtIBbXHtKMlrImbgduMcbkTPF8ypp+5QVZ2zRFzR4HjEhRs881yFf73tjqznSp39XLSZn3sjU7gz7jLMXrm1e36bi2nG4wxhQAnuO/q8WPAgVSnM8Ert3+Z6y1Ee7zdvqkdyBrbTyuq9/fMcYUNMb4GmNudxeDmUAHY0x9dwMyFPgXV+OTkb4E7jGuWzYF4Cr+SVeZJ02stQnAM8ArxphO7j2TRYBPcTWnE9x7XboDz+I6NHbu371Ae3PxHQkuYlwXVZQ3LiWBScAb1trwC6e11p4C5gMvG2OCjDENgRb8t9d3FlDLGNPWXQxHAD9ba/de+ztxkWVAoPtcsRzAE0Ak8KP79bQyxtzo/v9KwDBchxkvxwAPGmN8jDFtgLJALvf5d18CTxhjChnXxS+D0/F1iKQX1e3UyRJ1G8AY85Qx5g7jumWWvzGmL67TEn697MJTcM+TE1fv4udeziX7GOO6bVeou8bXxXXqwBx33q3AbuAFd9bOQGmuXDPTJBXfG8tx3SlnsPv3pBGu82CXX2GxPYwxAe5TC+7CVbP9yYTPOKtR8+ryIq7bG20HtgA/8d95fr/iOrS8z/x3fs8TwEPGmAhchyJmZlCuQbhO6N6Mq8i9Ahj3H15v4H1cFyw1Bu5yF5AMY63djOt2ULNx7S04guvQRGw6r+dTXK9vKK7bjPzmfup267oYqgGuczsnuQ8L/WOt/QdXE3kI1/k/GNeV8vdcZjVBuP7gI3B93t/hOmcO97wjjDGzU0zfB9chxhO4Lrroba3d7c57GNd5pm/hKkaVgAvvCfmCOf9+gZc6PeFK70k0/91Z4TSuiznap/jMWwLbjTGRuArw57iu/r2cRFwXuxzHVdRfxXXHhZK4vmh+xvU3sR5XAU7Xz1gkHahup0IWqtu41/kuro2LY7jOzW9vrT0IYIyZaowZx+X9H67D7R1wva/R/FfvmxhjTqSYtjuwF9dh+o9xnR+d8jPvBNTHVU9fAjq4G85zGpqL7/NaObXvh9uVvjdicV002Ak4C4wHurrPx76cWFynnsxxv/6mQOfM+oyzEuM6+iKSdsZ1K5CTQBFr7RGn80jqGGPKAdustam6TZBxXeQw2lpbNmOTiUhGU932TO7T4zpaa39MxbTZ/jPWnldJE2NMO+O6v2IIrj2Nv2TXPw5vZYzJZYxp5j6U9T9cV67Ovtp8IpI1qW5nf972Gat5lbTqhOtCiYO4Lja44jCN4pF8gNG4Djutx3WRy6tXnENEsjLV7ezPqz5jnTYgIiIiIh5De15FRERExGOkqXk1xizOqCAiIhnJG+uXN75mEckerlS/UnW18Tm5c+duXrNmTZ1nICKe6KzTATKbaraIeLDL1uw0Na9lypRhw4aMvp+yiEj6M8bsdjpDZlPNFhFPdaWarXNeRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpeRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpeRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpeRURERMRjqHkVEREREY+h5lVEREREPIaaVxERERHxGGpeRURERMRjqHkVEREREY+h5lUyhK+vL2FhYVSqVIm2bdty+vTpdF3+ypUr+fnnn5N/njx5MtOmTbumZU2dOpXDhw+naZ69e/dSqVKla1qfiEhWo5otnkTNq2SIwMBAtmzZwrZt28iXLx/vvfdeui7/wkLYr18/evTocU3LupZCKCKSnahmiydR8yoZrm7duhw6dCj557Fjx1KrVi2qVKnC8OHDkx9v3749NWrUoGLFikyZMiX58cWLF1O9enWqVq1K48aN2bt3L5MnT+btt98mLCyM1atX89JLL/HGG2+wc+dOateunTzv3r17qVKlCgAvv/wytWrVolKlSvTt2xdrLbNmzWLDhg3cd999hIWFER0dzcaNG2nQoAE1atSgefPmHDlyBICNGzdStWpV6tatm+6FXUQkq1DNlizPWpvqfzVq1LAiqREcHGyttTYhIcF27NjRLlq0yFpr7ZIlS2yfPn1sUlKSTUxMtK1bt7arVq2y1lr777//WmutjYqKshUrVrQnTpywx44ds8WKFbN//fXXedMMHz7cjh07Nnl9KX+uWrWq3bNnj7XW2tGjR9tXXnnlvHmttfb++++3c+fOtdZa26BBA7t+/XprrbVxcXG2bt269tixY9Zaa2fMmGF79eplrbW2cuXKduXKldZaawcPHmwrVqyYfm+YZDhgg01DvcsO/1SzJbVUsyWruVLN1p5XyRDR0dGEhYVxww03cPLkSZo2bQrA0qVLWbp0KdWqVaN69er8/vvv7N69G4B3332XqlWrUqdOHQ4cOMDu3btZu3Yt9evXp2TJkgDky5fvquvu3LkzX375JQAzZ86kS5cuAKxYsYJbb72VypUr8/3337N9+/aL5t21axfbtm2jadOmhIWF8eqrr3Lw4EHOnDnD6dOnadCgAQDdu3e//jdJRCSLUM0WT+LndADJns6dP3XmzBnatGnDe++9x6BBg7DWMmzYMB5++OHzpl+5ciXLly9nzZo1BAUF0bBhQ2JiYrDWYoxJ07q7dOlCp06duPvuuzHGUKZMGWJiYhgwYAAbNmygePHivPTSS8TExFw0r7WWihUrsmbNmvMeP336dJpziIh4CtVs8STa8yoZKjQ0lHfffZc33niD+Ph4mjdvzscff0xERAQAhw4d4tixY5w5c4a8efMSFBTE77//ztq1awHXuVerVq3i77//BuDkyZMA5MqVi/Dw8Euu8+abb8bX15dXXnkleQv+XNHLnz8/ERERzJo1K3n6lMsqW7Ysx48fTy6E8fHxbN++nTx58hAaGsqPP/4IwOeff56u75OISFagmi2eQHteJcNVq1aNqlWrMmPGDLp3787OnTupW7cuACEhIXz22We0aNGCyZMnU6VKFcqWLUudOnUAuPHGG5kyZQp33303SUlJFChQgGXLltG2bVs6duzInDlzGD9+/EXr7NKlC0OGDEkuoHny5KFPnz5UrlyZm266iVq1aiVP27NnT/r160dgYCBr1qxh1qxZDBo0iDNnzpCQkMDjjz9OxYoV+eSTT3jwwQcJCgqiefPmmfDOiYhkPtVsyeqM65zY1KlZs6bdsGFDBsYREckYxpiN1tqaTufITKrZIuKprlSzddqAiIiIiHgMNa8iIiIi4jHUvEqWd/ToUe69915KlSpFjRo1qFu3LrNnz05+/rHHHqNo0aIkJSWdN9+iRYuoWbMm5cuXp1y5cgwePDizo4uIeCXVbclIal4lS7PW0r59e+rXr89ff/3Fxo0bmTFjBgcPHgQgKSmJ2bNnU7x4cX744Yfk+bZt28ajjz7KZ599xs6dO9m2bRulSpVy6mWIiHgN1W3JaGpeJUv7/vvvCQgIoF+/fsmPlShRgoEDBwKum1hXqlSJ/v37M3369ORpxowZw3PPPUe5cuUA8PPzY8CAAZkbXkTEC6luS0ZT8ypZ2vbt26levfpln58+fTrdunWjQ4cOzJ8/n/j4eMC1BV+jRo3MiikiIm6q25LR1LyKR3nkkUeoWrUqtWrVIi4ujoULF9K+fXty587NrbfeytKlS52OKCIiKahuS3rTIAWSpVWsWJGvv/46+ef33nuPEydOULNmTRYvXsyZM2eoXLkyAFFRUQQFBdG6dWsqVqzIxo0bqVq1qlPRRUS8kuq2ZDTteZUs7c477yQmJoZJkyYlPxYVFQW4Dj19+OGH7N27l7179/L333+zdOlSoqKiGDJkCKNGjeKPP/4AXBcIvPXWW468BhERb6K6LRlNzatkacYYvv32W1atWkXJkiWpXbs2DzzwACNGjGDJkiW0bt06edrg4GBuv/125s2bR5UqVRg3bhzdunWjfPnyVKpUiSNHjjj4SkREvIPqtmQ0DQ8rIl5Bw8OKiHgODQ8rIiIiItmCmlcRERER8RhqXkVERETEY6h5lSxv8eLFlC1bltKlSzN69OiLnrfWMmjQIEqXLk2VKlXYtGlT8nMPPvggBQoUoFKlSpkZWUTEa6lmS0ZT8ypZWmJiIo888giLFi1ix44dTJ8+nR07dpw3zaJFi9i9eze7d+9mypQp9O/fP/m5nj17snjx4syOLSLilVSzJTOoeZUsbd26dZQuXZpSpUoREBBA165dmTNnznnTzJkzhx49emCMoU6dOpw+fTr59ir169cnX758TkQXEfE6qtmSGdS8SpZ26NAhihcvnvxzsWLFOHToUJqnERGRjKeaLZlBzatkaZe6D7ExJs3TiIhIxlPNlsyg5lWytGLFinHgwIHknw8ePEiRIkXSPI2IiGQ81WzJDGpeJUurVasWu3fv5u+//yYuLo4ZM2bQrl2786Zp164d06ZNw1rL2rVrCQ0NpXDhwg4lFhHxXqrZkhnUvEqW5ufnx4QJE2jevDnly5enc+fOVKxYkcmTJzN58mQAWrVqRalSpShdujR9+vRh4sSJyfN369aNunXrsmvXLooVK8ZHH33k1EsREcn2VLMlM5hLnXtyORonW0Q81ZXGyc6uVLNFxFNdqWZrz6uIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1r+LxDh06RLu7O7Fo0SKno4iIyFXExMTw8ICBvD72jUuOtiVyNWpexaNt2LCBqtVrsvpgHJ3u7cHb495VMRQRyaKOHj1Kndsb8NUPv/Laux9wb/cHiI2NdTqWeBg1r+KxZs36mttuv4P4m24j5NbOBN7Rk6eHDeOhvv1ISkpyOp6IiKSwfft2qlSrwZ+nEgm8/QECG/bl2yUrqNfgTk6ePOl0PPEgfk4HELlWu/b8xS0VKhF+4jf2T/iG0Hz5KVu+AgcOHcYY43Q8ERFJ4dDhI9xYqDA3JCax8/2H8M8RyC0VKhIfH09cXJzT8cSDaIQt8XiLFy+mZcuWdO/Zm2mffOh0HMmiNMKWSNYQERFBrly5KF6qDPv3/OF0HMmiNMKWiIiIiGQLal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppXEREREfEYal5FRERExGOoeRURERERj6HmVUREREQ8hppX8Xjx8fEA+Pn5OpxERESu5lzN9vX1cziJeCo1r+LxTp06BUCB/Dc4nERERK7m9OnTAITmCXU4iXgqNa/i8eYvWgpA2VvKOJxERESuZtmy5QCUL3uLw0nEU6l5FY924sQJvprxOQCdOnVyOI2IiFyJtZaxb78DwGOP9Hc4jXgqNa/i0T786GMA+vZ/hJCQEIfTiIjIlfzyyy/8+ft2SperyK233up0HPFQal7FY/32228MG/oMAI8PfMThNCIiciVnzpyhW/cHABjyxGMYYxxOJJ5Kzat4pF9//ZW69e4A4Omhz1K+fHmHE4mIyOX8+++/NGzSnL1//kG9BnfSq1dPpyOJB1PzKh7l4MGDPPf8C4SFhREZfob7H+jF6FGvOh2LIcOe58E+D2OtdTqKiEiWER4ezsSJkyhXsTJbNvxChcphLJw7G39/f0dzzZ07lzq3NyAqKsrRHHJtsu1N1g4cOMCKFSucjiHpICEhgdOnT7P0+1UsWTA3+fHnXniRV0a8lCUOPe3YsZPFC+dRvuwtDBn8lNNxRDxObGwsM2fOdDqGpIOkpCTOnj3Llq3b+OSjD5Ifb96qLV9/OZ3g4GAH07ns27ePdWvX0Knbfcyb/TU+PtqX50lMWvYU1axZ027YsCED46Sft94ex+CnnuKGKg0wuhGyRzu+adlFj/Xu3ZsSJUo4kObS/u/z6RwvfCvxv85nxrRPaNOmjdOR5ALGmI3W2ppO58hMnlSzN27cSM2aNclTpgb+ufI5HUeuw6Vqdtu2balVq5YDaS5t7dq1rD1qSTp5kD6d2zBm9CinI8kFrlSzs21XN2jgo8ydv5Bt4TkJvvNhp+PIdQhJ9CVy5yr8ct9I0M01Mb7+fL0zAnZudzpaMps/jNxVWxBfrCLd7u/Bz6tXUblyZadjiXiMGjVq8MLwlxg/dSY5Gz6Lj38OpyPJNbqhcGVOLh6PT3AegkrXxjdnCKv/hdWLs07NhlwEVr0D3+C8TP54CJUqlKdHj+5Oh5JUyrZ7XsF1ZWNYzVuJKNmI4OraEyaZI2L7Cvw2f8nWTRsoUKCA03HETXtesz5rLfd06caq348S0vKpLHFKkGR/ccf3cvbrF1m6cB633Xab03HE7Uo1O1uf5BEaGsp3SxYSv/Frov/a6HQc8RIhFRuRUPJ2mrduR0xMjNNxRDyGMYYvpk2lqF8EkWumOx1HvETAjTcR3HQgbe7qwL59+5yOI6mQrZtXgFKlSjHv26+JXPoOcSf2Ox1HvETwbd3YH5OD7j0f1B0IRNIgZ86cLF04D58/VxG5Y5XTccRLBN5cC9+wu7izWUvCw8OdjiNXke2bV4A77riD8W+/ScTckSRGnXE6jngBY3wIbjaI5Ws288qrI52OI+JRChUqxPLFC4n54UNiD+9yOo54iaDq7TidqyTt7+lMYmKi03HkCryieQXo1asnD/XoRuSCMdjEeKfjiBfw8c9JcJthjBk3ga++muV0HBGPUqVKFT6fNpXw+aNJOHvc6TjiBYwxhDTqy+a/j/LYk4OdjiNX4DXNK8AbY16ndrn/EfndZB3KlUzhl+sGcrUdRq8+fdm0aZPTcUQ8Srt27XjumcFEzBtJUly003HECxhfP4JbPc20mV/z/pQpTseRy/Cq5tXHx4dZM7/ghtjDRG2Y7XQc8RI5CpUmqFF/mrVqw+HDh52OI+JRhj49hNaN6hG55G2sTXI6jngB38Bc5Gr3PE89PYyVK1c6HUcuwauaV4CQkBCWL15A0m8LiNr9i9NxxEsElb0NyjejSYvWGo5QJA2MMXzy4RRK5TZErp7mdBzxEv75ihLc4kna39OJP//80+k4cgGva14BihcvzqL5c4n6bgJxx/5yOo54ieDaHTnqk5/O995PUpL2IImkVkBAAIvmfUuOgxuI+G2503HESwTeFIZ/7a40btaSU6dOOR1HUvDK5hWgdu3afDBpIuFzR5EYoV9KyXjGGEKaPMJPW3cz7LkXnI4j4lHy58/P8iULift5GjEHtjkdR7xEcNUWRBaoRJv29xAfr4u9swqvbV4BunXrymMD+hIx/zWS4mOdjiNewPj5E9J6KBM/nMpnn33udBwRj1K+fHm+mvEFEQvfIP70P07HES8RXL8XO49G0m/Ao7rYO4vw6uYV4JURL9GwZkUil0/QL6VkCt/gPORq9yz9Hh3EmjVrnI4j4lGaN2/OyBHDiZj7KkmxkU7HES9gfHwJaTmYWQuW8c67E5yOI6h5dQ1H+H+fUsTnLJFrZzodR7xEwI03EdR0IK3btddwhCJp9NigR+nUpjkRC9/AJulm8pLxfHIEEdLuOZ5/6WUWL17sdByv5/XNK0BgYKBrOMLdK4jcudrpOOIlgm6uhW9YOxo3b6XhCEXSaPLECVQonIuIVR85HUW8hH+eQoS0HkLnbvexY8cOp+N4NTWvboULF2bpwvlEr5pC7JE/nI4jXiKo+l2cDilBh45dNByhSBr4+fkx/9uvyXViJ5FbFjodR7xEzmIVyVHvAZq0aMWJEyecjuO11LymEBYWxmdTP3YPR6hfSsl4xhiCGz3Mxj1HeGLw007HEfEoefLkYfmShcSv/5LovVucjiNeIrhSY+KK30qLNncRG6uLvZ2g5vUC7du3Z9jgJ4iYP4qkuBin44gXML5+hLR+hqlffMmUDz5wOo6IRyldujRzvplFxOK3iP/3oNNxxEsE334/f0f40LN3H13s7QA1r5fw7NBnaNWgjoYjlExzbtTK1LEAACAASURBVDjCJ4cMZdWqVU7HEfEoDRo04O2xowmf+yqJ0WedjiNewBgfgps9zqJVa3nt9TFOx/E6al4v4dxwhCVzWSJ//MzpOOIl/G8oRnCLJ7nr7o4ajlAkjfr26cMD3ToRsWAMNlE3k5eM5xOQk5C2zzFq7FvMmTPH6TheRc3rZeTIkYPF8+cQcOAXIrd953Qc8RKBN4XhX6sLjZu34vTp007HEfEo494cS42bCxO5YooO5Uqm8Mudn1xthnLfA73YskXnXWcWNa9XkD9/fpYvXkjsT58Sc3C703HESwSHtSQyfwXaduhIQkKC03FEPIavry+zZ80kb+R+ojZpT5hkjhyFbyGwQV+atWrDP/9o5LfMoOb1KipUqMDMLz4jYsFYDUcomSa4wYNsP3yW/o8MdDqKiEfJlSsXyxcvIHHLXKL2rHc6jniJ4PJ3kFSmEU1btiE6OtrpONmemtdUaNmyJa8Mf4GIeSNJio1yOo54AePjS0irwXw5b4mGIxRJoxIlSrBg7rdELn2XuON7nY4jXiK4ThcOJ+WmW/cHdNpKBlPzmkqPPzaQjq2aErFIwxFK5vDJEUxIu+d4bvgIlixZ4nQcEY9St25dJk94l/C5o0iM1PnjkvGMMQQ3eZRVG7bz4ksjnI6Tral5TSVjDJMnTqB8wWAif5jqdBzxEv55ChHSajCdut7Lzp07nY4j4lG6d7+fAQ/1JGLBaGxCnNNxxAv4+OcgpM0wxr33PjNmzHQ6Tral5jUN/P39mf/t1wQd20rkr4udjiNeImfxSgTc1oMmzVvx77//Oh1HxKO8NvIV6lUpQ8Ty93QoVzKFb0hecrV7lof69WfdunVOx8mW1LymUd68efluySLi180get+vTscRLxFSuQmxxWrSos1dxMVpD5JIavn4+PDlF59RMOlfItfNcjqOeImAAqUIavwoLdvexYEDB5yOk+2oeb0GZcqU4duvvyJy0ZvEnzzkdBzxEsF39OCvs5ZeD/XVHiSRNAgKCmLZovmwcymRu35yOo54iaAyt+JTqRVNWrQmIiLC6TjZiprXa9SwYUPeHPOaezjCcKfjiBcwxofg5k+w4PufGD1mrNNxRDxK0aJFWbJgHtErJhP7j0awk8wRVLMD/+YoQscu95KUpOHm04ua1+vwcN++9OhyD5ELx2ATdTN5yXg+AYGEtHuOka+/wdy5c52OI+JRatSowacffUDEvNdICNf545LxjDEEN+7HL7/vY8gzQ52Ok22oeb1O77z1BtVKFiRypYYjlMzhl/vG5OEIt27d6nQcEY9yzz33MOTxR4mYN4qk+Bin44gXML7+hLR+hg8+/YJPPpnqdJxsQc3rdTo3HGGes38TtWme03HES+QoUpacd/SmacvWHD161Ok4Ih7lheefo+lt1Ylc+i7W6lCuZDzfoFBC2j3HwCee4scff3Q6jsczadlbWLNmTbthw4Z0DxESEuLIycwvvvgi9evXp0mTJte9rGLFinHkn3/wCcmPX94i5G/9JL4heTk46UF8AgLBx7WdkK/pAHIWK0/c8X2cXP4+ieEnAEtwxTsJva0rxpjkZR77+hUSo05TuPub151PsqeIn7+gaORufvnpB3LmzImvry+VK1cmPj4ePz8/HnjgAR5//HF8fNK+nTpq1CieffbZdMl5+vRpvvjiCwYMGJCm+V566SVCQkIYPHjwdWcwxmy01ta87gV5kIyq2eD5dTsmJobcoaEk+QXiE5wX3+C811W3EyNP8e+id0k4exySEvENLUDBTrpRvZwv+q+NxH0/gY3r1lKqVCkA1e3LuFLN9uo9ry+//HK6NK4Afn5+zJs7lxzE45v7Rs6s/TL5uYLdRlGk13iK9BpPzmLlSYqP5fg3rxBapyNF+06hcK/xxB76nYjNC5LnSYqJIO7oHpJiIok//U+6ZJTsJ7huNw7GB3OvezjCwMBAtmzZwvbt21m2bBkLFy5kxIhr+wIdNWrUJR+31qb5woPTp08zceLEa8ohklJ61e2cOXNSqGBB8uUOJrTOPQQUKn1ddfv06s/JeVMYRR6cQJGHJpG3Qc/rzijZT2CpGvhVv4fGzVtx9uxZ12Oq22mWpZrXlStX0rBhQzp27Ei5cuW47777ks8jXb9+PbfddhtVq1aldu3ahIeHExMTQ69evahcuTLVqlVjxYoVAEydOpX27dvTtm1bSpYsyYQJE3jrrbeoVq0aderU4eTJkwD07NmTWbNc9/276aabGD58ONWrV6dy5cr8/vvvABw/fpymTZtSvXp1Hn74YUqUKMGJEycumb927dpMGv8OSfs2EH9i/2VfZ+SOVeQoWp7AktUB8PHPSb6m/Tiz9r97EEbu+pnA0rUJLl+fqJ0/XOc7K9mVMYbgpgNZsf43Xnr5lfOeK1CgAFOmTGHChAlYa0lMTGTIkCHUqlWLKlWq8P777wNw5MgR6tevT1hYGJUqVWL16tUMHTqU6OhowsLCuO+++9i7dy/ly5dnwIABVK9enQMHDhASEpK8rlmzZtGzZ08Ajh49SocOHahatSpVq1bl559/ZujQoezZs4ewsDCGDBkCwNixY5OzDB8+PHlZI0eOpGzZsjRp0oRdu3Zl8Dso18uT67aPjw+zZk4n5oeP8Q3KTfypw5d9nVer24mRJ/HNlT95+oACJa/3rZVsKrh6G8Lz3kK7uzuRkHD+xd6q26njl+FrSKPNmzezfft2ihQpQr169fjpp5+oXbs2Xbp0YebMmdSqVYuzZ88SGBjIO++8A8Bvv/3G77//TrNmzfjjjz8A2LZtG5s3byYmJobSpUvz+uuvs3nzZp544gmmTZvG448/ftG68+fPz6ZNm5g4cSJvvPEGH374ISNGjODOO+9k2LBhLF68mClTplwxf/fu9/P2uLfZsWc/NiEegKPTnwUfH4yvP4V7vEX8iX0EFCp93nz+eQtj46JJio3CJ0cQUTtXEVqvG77BeTn+7WuE1u2cHm+vZEM+/jkIbjOMt8YPuagQlipViqSkJI4dO8acOXMIDQ1l/fr1xMbGUq9ePZo1a8Y333xD8+bNee6550hMTCQqKoo77riDCRMmsGXLFgD27t3Lrl27+OSTT666JT5o0CAaNGjA7NmzSUxMJCIigtGjR7Nt27bk5S1dupTdu3ezbt06rLW0a9eOH374geDgYGbMmMHmzZtJSEigevXq1KhRI2PeOEk3nly3K1SowPTPptGhYydylm+U/Hha63auaq05PncM4ZvmE1gijODKTfDLdUN6vL2SDQU36M3WOa8w8LEnL3pOdfvqslzzWrt2bYoVKwZAWFgYe/fuJTQ0lMKFC1OrVi0AcufODcCPP/7IwIEDAShXrhwlSpRILoKNGjUiV65c5MqVi9DQUNq2bQtA5cqVL3uF9t133w24bqfyzTffJK9j9uzZALRo0YK8efNeNnujRo2Sz13JX7AQG797D3AdfvINCr1ganPxAtwSI08Rf+oIOYpVxBiD8fEl7vheAm686bLziHfzC8lHrrbPcmTqY2zYsIGaNf87TejcXrClS5eydevW5L1WZ86cYffu3dSqVYsHH3yQ+Ph42rdvT1hY2CXXUaJECerUqXPVLN9//z3Tpk0DXOdyhYaGcurUqfOmWbp0KUuXLqVatWoAREREsHv3bsLDw+nQoQNBQUEAtGvXLo3vhDghO9TtKhUr8ueRnSTFRgFpr9uBpWpQ9OEPifl7I9F/beTI1Mco0vu9SyxDBIyvHyGthvDFzKeJT7j4Vpuq21eW5ZrXHDlyJP+/r68vCQkJWGvPu5DpnCtdbJZyOT4+Psk/+/j4XLR36sJ5zq33auu40IoVK8if33XYKDIykpp16nE6LvKi6fzz/4/YA9vOeyz+9D+YgEB8cgRx9rflJMVEcGhyb1eGuCgid65W8ypXFFDwZvD1p0XrtmzZuJ5ixYrx119/4evrS4ECBbDWMn78eJo3b37RvD/88AMLFiyge/fuDBkyhB49elw0TXBw8Hk/p/ybjIlJ2y2HrLUMGzaMhx9++LzHx40bd8m/dcnaskPdttbSo1dvFix565LTXa1uA/gG5iK4QkOCKzTk2KwRxBzYRnDZeqnOIt7FJ2cIIe2e5+yUvixfvjz5XG7V7avLUue8Xk65cuU4fPgw69evByA8PJyEhATq16/P559/DsAff/zB/v37KVu2bLqu+/bbb+fLL10n8S9duvSiLZHLCQ4OZvniBRAXQ/Se9ec/V6EhsQd3EL3XtSs+KT6WU8vfJ/et9wAQuXMVBTuPoFj/jynW/2MKPTBO571KqhgfX6jYkiYtWrN371769evHo48+ijGG5s2bM2nSJOLjXaez/PHHH0RGRrJv3z4KFChAnz596N27N5s2bQLA398/edpLKViwIDt37iQpKSl5LxdA48aNmTRpEgCJiYmcPXuWXLlyER7+30h0zZs35+OPP06+Wv3QoUMcO3aM+vXrM3v2bKKjowkPD2fePN1+zlN5Wt02xvDRlMmUyeuPjbn4LgpXq9vR+35Nvm9sUmwU8aeO4Jf7xvR6SZJN+ectgvEL4J7OXdm1axfHjx9X3U6FLLfn9VICAgKYOXMmAwcOJDo6msDAQJYvX86AAQPo168flStXxs/Pj6lTp5635Z4ehg8fTrdu3Zg5cyYNGjSgcOHC5MqVK1XzFi1alIIFC3Bm9ScEFCjp2jOG6xzFG+95gZPLJnNy2SRISiK4YiNyVW9DwpmjJJ49TkCRcsnL8c9TCJ8cQcQe3kWOIulb5CV7sQlxnN2+ktNnj1GlSlVeeOF5nnrqKQAeeugh9u7dS/Xq1bHWcuONN/Ltt9+ycuVKxo4di7+/PyEhIcmHjvr27UuVKlWoXr06I0eOvGhdo0ePpk2bNhQvXpxKlSolF7R33nmHvn378tFHH+Hr68ukSZOoW7cu9erVo1KlSrRs2ZKxY8eyc+dO6tatC7huu/TZZ59RvXp1unTpQlhYGCVKlOCOO+7IpHdO0psn1u2AgAAWzfuWAgUKErVjJblq3pX83JXqNkDcP39yatlk8PEFm0Suqs3IUfiWdH1dkj3ZxAQiEwyVKlem7C230LNnT5580nUurOr2pWWJ+7xmZbGxsfj6+uLn58eaNWvo379/8snLqfXll1/Re8BAQruOxS8kXwYlFfmPTYjn7Dcv8lCn1ox9/TWn42QJus+r97jeur1r1y5q161HYIunyPm/KhmYVOQ/Eas+oYzfCX74fhn+/v5Ox3Gc7vN6Hfbv30+tWrWoWrUqgwYN4oMPPkjzMjp37sRTgx4hct5r2IS4DEgpcj7j509Im6G8N3kKq1evdjqOSKa63rpdtmxZZs2cTsSiN4k/dSSDUoqcL/iOHuw8cobXXh/jdJQszyNOG3BSmTJl2Lx583Uv5+72d/Haa6MJjg7X7VMkU8Ts20pIcCBlypRxOopIpkqPul27dm1yhQQTf2I//nkLp1MykctLOHWYhBP7adSgvtNRsjztec0Ex44do1mrNuRuOkCNq2SK2CN/EL1qCksXzqdQoUJpmnfx4sWULVuW0qVLM3r06IueP3XqFB06dKBKlSrUrl2bbdtcV2Dv2rWLsLCw5H+5c+dm3Lhx6fJ6RDJTQkICd93difhCVQgsXdvpOOIFEqPPEj53JOPeHHNN54x6W91W83qBBx98kAIFClCpUqV0WV5MTAzNWrUloeTtBFdomC7LFLmShLMnCJ8/ms+mfnzZe/9dTmJiIo888giLFi1ix44dTJ8+nR07dpw3zahRowgLC2Pr1q1MmzaNxx57DHAdat2yZQtbtmxh48aNBAUF0aFDh3R7XSKXkt41G+DRQU+w9cBJghs+pFu3SYazifFELHidXvd15qHevdM8vzfWbTWvF+jZsyeLFy9Ol2VZa7n/gV4ciM1J8G3d0mWZIleSFBdDxLxRDBv8BO3bt0/z/OvWraN06dKUKlWKgIAAunbtypw5c86bZseOHTRu3Bhw3Q5p7969HD169LxpvvvuO26++WZKlChx7S9GJBXSs2YDvDdxItNnzyWk1RCMr86sk4xlrSXy+/epVaYYb429tnNdvbFuq3m9QP369cmXL33uCPDyqyNZvnYLwc0GYYzeaslY1iYRuXQcLRvcyrNDn7mmZRw6dIjixYsn/1ysWDEOHTp03jRVq1ZNHslo3bp17Nu3j4MHD543zYwZM+jWTRtskvHSs2YvW7aMZ557kZB2z+OTM+TqM4hcp8iNc8gXfZBvvpqBr6/vNS3DG+u2OqoM8tVXs3jjnQmEtBmGj39Op+OIF4j86XNKhiQx9aMPrvlQ56VunXfhsoYOHcqpU6cICwtj/PjxVKtWDT+///ZQxcXFMXfuXDp16nRNGUScsGvXLjp26UZIy6fwz1vE6TjiBaL+XEfS1nksX7yAkJBr31jyxrqtYyIZYOPGjfTq05dc7YfrAi3JFJHbviNg/1oWb9pwXTd8L1asGAcOHEj++eDBgxQpcv4Xee7cufnkk08AV9EsWbIkJUuWTH5+0aJFVK9enYIFC15zDpHMdPLkSZo0b0VAnXt1X1fJFHHH/iZq+QSWL17A//73v+taljfWbe15TWeHDh2ieeu2BDXqT45CpZ2OI14g5uB2Yn/6lOWLF5I/f/7rWlatWrXYvXs3f//9N3FxccyYMYN27dqdN83p06eJi3Pdr/jDDz+kfv365M6dO/n56dOne8yhJ5H4+Hhat+tAdOEwgqtcPH68SHpLjDxF+LxRvP/eeOrUqXPdy/PGuq09r+koKiqKpi3bQPlmBJW9zek44gUSzhwlYsFYvvriMypUqHDdy/Pz82PChAk0b96cxMREHnzwQSpWrMjkyZMB6NevHzt37qRHjx74+vpSoUIFPvroo+T5o6KiWLZsGe+///51ZxHJaNZaevftx66TceRq84DTccQL2IQ4IuaP5tG+vbnvvnvTZZneWLc1POwFunXrxsqVKzlx4gQFCxZkxIgR9E7FrSuSkpJo2+Ee1u6PJKTZIN1eRTJcUmwUZ758hpeGPM4Tjw9yOk6Wp+Fhs6drrdkAb7z5Fq++PZHcnV7DJ0dQBicVb2etJWLJOOqVysO3s77Ex0cHv6/kSjVbe14vMH369Guab9hzL/Dzb3vIfffLalwlw9mkRCIWv8ndLZvw+GMDnY4j4phrrdkLFizgpZGjydPldTWukikif/mKwpxi5ufz1LheJzWv6eD//u8zJn44lTxdx2D8/J2OI14g8oeplL8xkCmT3tPGkkgabdu2ja73dSekzbP4hRZwOo54gcjff8TsWs6yTRsIDAx0Oo7HU/N6ndasWUP/gY+R+56X8Q3O43Qc8QKRvy4m6OhW5m9aj7+/NpZE0uLYsWM0adGKnPUfJGex8k7HES8Qe2Q30SvfZ/WK7yhcuLDTcbIFNa/XYd++fbRu157gpgMJuPEmp+OIF4je9ytxv0zn51/WkDdvXqfjiHiUmJgYmrduR8JN9Qip0MjpOOIFEsJPED7/NaZ9/CHVqlVzOk62oZMurlF4eDiNm7XEN+wuAm+u5XQc8QLxJw8RuehNvv36K8qUKeN0HBGPYq2lR6/e7I/JQXC99LnKW+RKzg3X/cyTj3H33Xc7HSdbUfN6DRITE2l/T2dO5ypJUPV2V59B5DolxkQQPnckb7w+ikaNtMdIJK1eHTmKpT9t1HDdkilcw3W/Q/Pba/H8s8OcjpPt6C/4Gjz+1BA2/32U4EZ9dbGMZDibmEDkgtfp3rkD/R5+2Ok4Ih7n66+/Zszb4wlpq+G6JXNE/vwFJYLimfbJR+oTMoDOeU2jKR98wKczZhHaZQzGV2+fZCxrLZErp1CtZEHefftNp+OIeJxNmzbxQO8+7uG6r28EOpHUiNy+Av+9P7PkOofrlstT95UGK1eu5MkhQwnt/Bq+gbmcjiNeIHLTfELP/s3s79bi6+vrdBwRj3L48GGat2pDYMOHNVy3ZIqYgzuJ+fET1v74AzfeeKPTcbItNa+p9Oeff9L+nk4Et3gS/3xFnY4jXiB6zwYSN3/DdxvWnTcGtYhc3bnhum25poSUu93pOOIFEs4cJWLhGL78/P+oWLGi03GyNTWvqXDq1CkaN2uJf+2uBN4U5nQc8QJxx/cRuexdliyYy0033eR0HBGPkpSURJf7uvOPyUfIrZ2cjiNeICk2ivC5Ixnx/DBatWrldJxsTxdsXUV8fDxt2t9DZIFKBFdt4XQc8QKJUWcInzeSie+Oo169ek7HEfE4zz7/Ij9u2UVIk0d0sYxkuOThulvcyZNPPO50HK+g5vUq+j8ykJ1HIwmu38vpKOIFbEI8EfNH069XD3r06O50HBGP89lnn/PeBx8T0mYoxi/A6TjiBSJXf0rZ/DmZMnmiNpYyiU4buIJx74znq/lLCe0yGuOji2UkY1lrifjuPepWLMXrr410Oo6Ix1mzZg39HhlI7o6v4BusEegk40VuXULgkS0s0HDdmUrN62UsXryY5196mdAuo/HJEex0HPECkeu/oUDCcb6a8S0+PjooIpIW54brDtJw3ZJJYvZtJW7tF/y09mfy5cvndByvoub1Enbs2EHnbvcR0voZ/PMUcjqOeIGoP9bA9kUs37SB4GBtLImkRXh4OE1atManajuCStd2Oo54gfiTh4hY9AZzvv6SW265xek4XkfN6wVOnDhBkxatyFHvAXIW060uJOPFHd1D1PcTWbl8KUWL6jZsImmRmJjI3Z26cjKoOCE17nI6jniBc8N1j3ntVRo3bux0HK+kY5MpxMXF0aLNXcQVr01wJf1CSsZLiDhJ+LxRfDzlfWrWrOl0HBGP88Tgp9nw5yFC7nxYF8tIhrOJCUQuHMt999zFgP79nY7jtdS8ullr6dm7D3+HG4Jv11XekvGS4mOJmPcaTw4cQOfOuhelSFpN+eADpn7xJSGtn8H46mIZyViu4bo/pMr/bmDCu287Hcer6bQBt9deH8PClWvI3WkUxqinl4xlrSVy2XjuvLUKL734gtNxRDzOqlWrXMN1dxqJb6BGoJOMF7V5PrnP/Mnc737RcN0OU/MKzJkzh9fGvkVolzH4BAQ6HUe8QOSa6RTzj+SLaVN1qFMkjfbs2cNdd3d0Ddd9Q3Gn44gXiP5rIwmbvuG79b9ouO4swOub1y1btnD/A70Iafc8frnzOx1HvEDkjlX47vmBpZs2kDNnTqfjiHiU06dPu4brrtVZw3VLpog7sZ/Ipe+weP4cSpYs6XQcwcub13/++YdmrdqQs0EfchQp63Qc8QKxh3cRs/ojflq1goIFCzodR8SjJCQk0LZDRyJuKE9ImMaPl4yXGHWGiLkjmTDuLW6//Xan44ib157cGR0dTbNWbUkq05Dg8vWdjiNeIOHsccLnj+bzTz+hSpUqTscR8Tj9Hx3I9sNnCW7Y2+ko4gVsQjwRC16nzwP30rPnA07HkRS8snm11nJfj54cSgwhuE5Xp+OIF0iKiyZi3kiee2Yw7dq1czqOiMd5d/wEvpy7mJBWgzVct2Q4ay0R30+iTvmbGPv6aKfjyAW8snl98aURrFi/jeAmA3WxjGQ4a5OIXPI2rRvVY+jTQ5yOI+Jxli5dyrMvvkRI2+c0XLdkisj133Bj3D/MmvmFhuvOgrzunNcZM2Yy7r33ydN1LD7+OZyOI14gcvU0bg714ZMPp2hjSSSNfv/9dzp26UZIq6fxz1vY6TjiBaL+WIPdvojlG9druO4syqua13Xr1vFQv/7k7jAC35C8TscRLxDx23JyHNrAok0bCAgIcDqOiEf5999/ady8JQG3dSdn8UpOxxEvEHf0L6K+n8iKZUsoVqyY03HkMrymeT148CAt27QjqPEjBBQs5XQc8QIxB7YR9/M0flzzEzfccIPTcUQ8yrnhumOL1iSkclOn44gXSIw4Rfi8UXw4eRK1atVyOo5cgVecyBEZGUmTFq0xlVoRVKaO03HEC8Sf/oeIhW/w1YwvKFeunNNxRDyKtZZeD/Vlzxmr4bolUyTFxxIxfxSPP/IwXbt2cTqOXEW2b16TkpK4p0s3TgQUJrjW3U7HES+QFBtJxNxXGTliOM2bN3c6jojHeX3sGyxY8RMhLZ7QnQUkw1lriVw+nka1KvPyS8OdjiOpkO2b1yHPDOOXnfsIbtxPF8tIhrNJiUQsfIPO7Vrw2KBHnY4j4nHmzp3Lq6PHEtL2WQ3XLZkicu0MivpG8LmG6/YY2fqc108+mcoHn35OaNcxGF9/p+PINdr3epuLHvMLyFp3ijDGhxs6vULcHz9QsUhuJk0Y73QkEY+zdetW7uvRk1ztnscvdwGn48g1OjZrBNF71p/3mJ9/AGSxxjBvw57YnKH47F7J0k0bCAzUxpKnyLbN648//sjD/fqTu2l/Es4cI+HMMacjyTUqeO9obEI8SdFnSDiyi6jt3xMWVo0P359I6dKlnY4HQPtOXVn94/9RwCeKeRvX4eeXbf+0RDLE0aNHadKiFT63NATjQ+yR3U5HkmuUu04nctVqT1J0OAn/7idh5wpCg3Pw/sQJ3HnnnU7HA2DixIm8+M7HmLNH+GnVCgoVKuR0JEmDbPsNe+DQEUqWuQX2fgd7nU4j1yMhIZHw8LNEhp8lqEIjbug6mj37f6Xu7fV5640x9O3Tx+mI+Pn64XNyH8s3ridPnjxOxxHxOBEREeS7IT82ag9s3ON0HLkOSUmWyIhwzvx7nJCba5Kj6UDi4qLo2qM33TrdzcQJ7zq+gR8QEEDUwZ18O3s2VatWdTSLpJ2x1qZ64po1a9oNGzZkYByRy9u/fz/vTZrM+1M+xK9gafxr3EPkkrd58ZknGfzUk45m+2XdenwMur1KFmaM2Witrel0jsykmi1OOnv2LJ9Om8ab4yZwOjqenI36Ebf+K+pVKM63s77E19e5i/EOHjzIxo2buOsuDdedVV2pZqt5FY8TFxfHE4Of5rNZcwhqWGSiMgAACz1JREFUMpDwBa/zzfT/o1mzZk5HkyxMzauIM6y1fPrpNB59/ElCWg4mdv1X9O/amlGvvuJ0NMnCrlSzs/3dBiT7CQgI4L13xzHwoR5Er3yfoPoPcX/PBzl79qzT0URE5ALGGHr2fICvpn9GxOI3yVG7E+9MmMjmzZudjiYeSs2reKxXRrxE3arlSDz7D4n5y/D+lA+cjiQiIpfRsmVLhj83jIStC8hZuzNDX9A9VeXaqHkVj2WMcRXCbUvxr9KScePfIy2nwYiISOYa0L8fMQe2k+N/VVi9ahUHDhxwOpJ4IDWv4tFuvfVWCuXPi42LITwymj17dJWyiEhWFRwczP333U/s76sIKV2T5cuXOx1JPJCaV/Foxhj6PdSLpL9/IWfh0ujiFBGRrK1P754k/b2OhLwl+fmXdU7HEQ+k5lU8XvHixfGJjSAhpAD79+93Oo6IiFxB0aJFiYs4jV/ewuzes9fpOOKB1LyKxwsODsYmxGCNH/Hx8U7HERGRKwgODiYhLgbj46uaLddEzat4PGMMBiBrDZstIiKXYMy5Yq2iLddGzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjHUPMqIiIiIh5DzauIiIiIeAw1ryIiIiLiMdS8ioiIiIjH8HM6gMi12rZtG9u2bWPLli3Jj23dupUZM2ZQtGjR/2/vfmOquu84jn8OclGQq4aFRQhmrV6wDC66SgUWaaGxQlBRoSb4jxk1rmqCm8MSNeGBSWmfNNFNK9U6jZpczIiCW4SlG6ZNDFQQkahtwqIY/jhXm5AGU6fQswcbpNc/De24nPuT9ysh8V5+R76HBx8/OZ57fsrMzHRwOgDAd929e1cXLlzQgwcPht/78l//VFVVlcLCwlRQUODgdDAJ5RXG2n/gA3304SFFvZSmsJeyNCnyp7rQ/mf9ae06/WJeqi5/1uj0iACA/2m4cEGrV62SO262pqUVKGy6R33/iNGvNm9V6LeP9Prrr2vatGlOjwkDUF5hrA/+sF8DAwM6U/d3hU+frZBJbv374QNl/HKB/lJ7xunxAADfsaqoSF9//bV+9/YuKebnComYom8n/0RRU936W/15iitGjPIKY7lcLv3xyIfy7vu9yveWaUK4WwV5C3X40EG5XC6nxwMAPObXmzcrIT5eywtX6sHUGM2KjtRfW1sUHR3t9GgwCOUVRrMsSzt+u10p3iT19PaqeN06WZbl9FgAgGfIzs5Wy2eNOnO2Rr/ZXqKJEyc6PRIMQ3nFc2HhwoVOjwAAGKH4+HiVvb3T6TFgKB6VBQAAAGNQXgEAAGAMyisAAACMQXkFAACAMSivMNLZs2dlWZa++OILSVJvb6/efPNNh6cCADwNmY3RRHmFkXw+nxYsWKCqqipJUmxsrKqrq59YNzAwMNajAQAeQ2ZjNFFeYZz+/n5dvHhRR48eHQ7Czs5OJScnS5KOHz+ulStXaunSpVq0aJGTowLAuEdmY7TxnFcYp6amRrm5uUpISFBUVJRaW1sVFRXlt6axsVHt7e1PvA8AGFtkNkYbV15hHJ/Pp6KiIklSUVGRfD7fE2veeOMNQhAAggCZjdHGlVcY5auvvlJDQ4OuXbsmy7I0ODgoy7K0detWv3WTJ092aEIAwBAyG4HAlVcYpbq6WsXFxbp9+7Y6OzvV1dWlF198Ud3d3U6PBgB4DJmNQKC8wig+n08rVqzwe6+wsFAVFRUOTQQAeBYyG4Fg2bY94sWpqal2S0tLAMcBgMCwLOuybdupTs8xlshsAKb6vszmyisAAACMQXkFAACAMSivAAAAMAblFY6YMGGC5s6dO/z13nvvjfjYH7sn9qZNm3Tjxo0ffBwAjHdkNoIJH9iCIyIjI9Xf3+/0GBhH+MAW8OOR2RhrfGALxnjhhRe0e/duZWRkKDU1Va2trcrJydGsWbNUWVkpyX9P7OvXr2v+/PmaO3euUlJS1NHRofv372vx4sWaM2eOkpOTdfr0aUlSVlaWhv4h9/l88nq9Sk5OVllZ2fDPj4yM1J49ezRnzhylp6fr7t27Y/wbAABzkNlwAuUVjvjmm2/8/gtqKKwkacaMGWpsbFRmZqbWr1+v6upqNTU1qby8/Im/p7KyUtu3b1dbW5taWloUFxen+vp6xcbG6urVq7p27Zpyc3P9junt7VVZWZkaGhrU1tam5uZm1dTUSJLu37+v9PR0Xb16Va+++qqOHDkS2F8EABiAzEYwYXtYOCI8PFxtbW1P/V5+fr4kyev1qr+/X263W263W5MmTVJfX5/f2oyMDL3zzjvq7u5WQUGB4uPj5fV6VVpaqrKyMi1ZskSZmZl+xzQ3NysrK0vR0dGSpDVr1ujTTz/V8uXLFRYWpiVLlkiS5s2bp48//ni0Tx0AjENmI5hw5RVBZ+LEiZKkkJCQ4T8PvR4YGPBbu3r1ap07d07h4eHKyclRQ0ODEhISdPnyZXm9Xu3atUt79+71O+b77vN2uVyyLEvSfz+g8PjPAwD4I7Mx1iivMNrNmzc1c+ZMlZSUKD8/X+3t7ert7VVERITWrl2r0tJStba2+h2TlpamTz75RPfu3dPg4KB8Pp9ee+01h84AAMYPMhujgdsG4Iih+6eG5Obm/qBHrww5ffq0Tp06JZfLpenTp6u8vFzNzc3auXOnQkJC5HK5dOjQIb9jYmJi9O677yo7O1u2bSsvL0/Lli37v88JAJ5XZDaCCY/KAjAu8KgsADAHj8oCAADAc4HyCgAAAGNQXmGE+vp6zZ49Wx6P56n3Wdm2rZKSEnk8HqWkpAzf8N/V1aXs7GwlJiYqKSlJ+/fvH+vRAWDcIbMRSJRXBL3BwUFt27ZNdXV1unHjhnw+3xP7XdfV1amjo0MdHR06fPiwtmzZIkkKDQ3V+++/r88//1xNTU06ePAge2UDQACR2Qg0yiuC3qVLl+TxeDRz5kyFhYWpqKhItbW1fmtqa2tVXFwsy7KUnp6uvr4+3blzRzExMXr55ZclSW63W4mJierp6XHiNABgXCCzEWiUVwS9np4ezZgxY/h1XFzcE2E2kjWdnZ26cuWK0tLSAjswAIxjZDYCjfKKoPe0x7kN7agy0jX9/f0qLCzUvn37NGXKlNEfEgAgicxG4FFeEfTi4uLU1dU1/Lq7u1uxsbEjXvPo0SMVFhZqzZo1KigoGJuhAWCcIrMRaJRXBL1XXnlFHR0dunXrlh4+fKiqqirl5+f7rcnPz9eJEydk27aampo0depUxcTEyLZtbdy4UYmJidqxY4dDZwAA4weZjUBje1gEvdDQUB04cEA5OTkaHBzUhg0blJSUpMrKSknSW2+9pby8PJ0/f14ej0cRERE6duyYJOnixYs6efKkvF7v8NaGFRUVysvLc+x8AOB5RmYj0NgeFsC4wPawAGAOtocFAADAc4HyCgAAAGNQXgEAAGAMyisAAACMQXkFAACAMSivAAAAMAblFQAAAMagvAIAAMAYlFcAAAAYg/IKAAAAY1BeAQAAYAzKKwAAAIxBeQUAAIAxKK8AAAAwBuUVAAAAxqC8AgAAwBiUVwAAABiD8goAAABjUF4BAABgDMorAAAAjGHZtj3yxZb1paTbgRsHAALmZ7ZtRzs9xFgiswEY7JmZ/YPKKwAAAOAkbhsAAACAMSivAAAAMAblFQAAAMagvAIAAMAYlFcAAAAYg/IKAAAAY1BeAQAAYAzKKwAAAIxBeQUAAIAx/gMQdQ4VWejt8QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 864x360 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "ThermalReactivation.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
}
