{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Thermal Reactivation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import PFAS_SAT_ProcessModels as pspd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "from IPython.display import Image\n",
    "import pandas as pd\n",
    "pd.set_option('display.max_colwidth', 0)"
   ]
  },
  {
   "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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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></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            \n",
       "1   0.97    fraction            \n",
       "2   0.97    fraction            \n",
       "3   0.97    fraction            \n",
       "4   0.90    fraction            \n",
       "5   0.97    fraction            \n",
       "6   0.97    fraction            \n",
       "7   0.97    fraction            \n",
       "8   0.98    fraction            \n",
       "9   0.99    fraction            \n",
       "10  0.33    fraction            \n",
       "11  0.33    fraction            \n",
       "12  0.33    fraction            \n",
       "13  0.33    fraction            \n",
       "14  0.33    fraction            \n",
       "15  0.33    fraction            \n",
       "16  0.33    fraction            \n",
       "17  0.33    fraction            \n",
       "18  0.33    fraction            \n",
       "19  0.33    fraction            \n",
       "20  0.15    fraction            \n",
       "21  0.90    fraction            "
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ThermalReactivation = pspd.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.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Moisture flow</td>\n",
       "      <td>kg</td>\n",
       "      <td>100.0</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.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>PFOA</td>\n",
       "      <td>μg</td>\n",
       "      <td>900000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>PFOS</td>\n",
       "      <td>μg</td>\n",
       "      <td>1500000.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>PFBA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>PFPeA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>PFHxA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>PFHpA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>PFNA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>PFDA</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>PFBS</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>PFHxS</td>\n",
       "      <td>μg</td>\n",
       "      <td>0.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.0    \n",
       "2   Moisture flow  kg   100.0    \n",
       "3   VS flow        kg   829.8    \n",
       "4   Carbon flow    kg   630.0    \n",
       "5   PFOA           μg   900000.0 \n",
       "6   PFOS           μg   1500000.0\n",
       "7   PFBA           μg   0.0      \n",
       "8   PFPeA          μg   0.0      \n",
       "9   PFHxA          μg   0.0      \n",
       "10  PFHpA          μg   0.0      \n",
       "11  PFNA           μg   0.0      \n",
       "12  PFDA           μg   0.0      \n",
       "13  PFBS           μg   0.0      \n",
       "14  PFHxS          μg   0.0      "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "IncominWaste = pspd.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": "iVBORw0KGgoAAAANSUhEUgAAAWAAAADrCAYAAABXYUzjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAodElEQVR4nO3deXhN5/7//+fOPIp5ShxECDJPSBEzMYW0huAYqspRU4ejqJ6qtoYe2tIa6ocqjibUFKoiNeugxphVDoIEUVNkTvbO/f3Dr/t8UrSiSVZ29vtxXbmuZu37vtd77UtflrXudS+dUgohhBClz0LrAoQQwlxJAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsBBCaMSqKI2rVq2q6tWrV0KlCCFE+XT06NHbSqlqv99epACuV68eR44cKb6qhBDCDOh0uiuP2y6XIIQQQiMSwEIIoREJYCGE0IgEsBBCaEQCWAghNCIBLIQQGpEAFkIIjUgACyGERiSAhRBCIxLAQgihEQlgIYTQiASwEEJoRAJYFCtLS0v8/f3x9vamZ8+e3L9/v1jHnzlzZqHfn3vuuWcaJykpia+++qrI/YYNG8b69eufaZ9C/J4EsChW9vb2JCQkcPr0aSpXrszChQuLdfzfB/CPP/74TOM8awALUZwkgEWJCQ0NJSUlBYCLFy8SHh5OUFAQrVu35vz58wBs3bqV5s2bExAQQMeOHUlNTQUgIyODF198ER8fH3x9fdmwYQOTJ08mOzsbf39/Bg0aBICTkxMAUVFRbNu2zbjv385Uk5KSaN26NYGBgQQGBhoDe/LkyRw4cAB/f38++eQTDAYDEydOJCQkBF9fX5YsWQKAUoqxY8fi6elJx44duXXrVul8ecI8KKWe+icoKEgJ8UccHR2VUkrp9XrVp08ftX37dqWUUu3bt1cXLlxQSil18OBB1a5dO6WUUnfv3lUFBQVKKaWWLl2qXn/9daWUUm+++aaaMGGCcdy7d+8WGv/3+9u4caMaMmSIUkqp3Nxc5ebmprKyslRmZqbKzs5WSil14cIF9duf4T179qju3bsbx1myZIl6//33lVJK5eTkqKCgIHXp0iW1YcMG1bFjR6XX61VKSopycXFRX3/9dXF8VcKMAEfUYzK1SAuyC/FnfjtDTUlJoUmTJnTq1ImMjAx+/PFH+vbta2yXm5sLQHJyMv379+fGjRvk5eVRv359AHbu3ElMTIyxfaVKlf5wv127dmXChAnk5uYSFxdHWFgY9vb2pKWlMXbsWBISErC0tOTChQuP7R8fH8/JkyeN13fT0tJITExk//79DBgwAEtLS2rXrk379u3/0vcjxP8llyBEsfrtGvCVK1dQSrFw4UIKCgqoWLEiCQkJxp9z584BMG7cOMaOHcupU6dYsmQJOTk5z7RfOzs72rZty44dO1i7di39+/cH4JNPPqFGjRqcOHGCI0eOkJeX99j+Sik+++wzY32XL1+mc+fOz/YlCPGUJIBFiXBwcODTTz/lo48+wsHBgfr16/P1118DD8PuxIkTwMMzTVdXVwBWrlxp7N+pU6dCN/Du3bsHgLW1Nfn5+Y/dZ//+/VmxYgUHDhwgPDzcOH6tWrWwsLBg9erVGAwGAJydnUlPTzf27dKlC4sXLzaOfeHCBTIzMwkLC2Pt2rUYDAZu3LjBnj17iuX7EQIkgEUJCggIwNfXl+joaNasWcPy5cvx8/PDy8uL2NhYAN5991369u1LUFAQVatWNfZ9++23uXfvHt7e3vj5+RmDb+TIkfj6+hpvwv1fnTt3Zt++fXTs2BEbGxsAXnnlFVauXImfnx/nz5/H0dERAF9fXywtLfHz8+OTTz5hxIgRNG3alMDAQLy9vRk1ahR6vZ7IyEgaNmxI06ZNGTJkCKGhoSX9tQkzont4ffjpBAcHK3kppxBCFI1OpzuqlAr+/XY5AxZCCI1IAAshhEYkgIUQQiMSwKLMSk1NZeDAgbi7uxMUFERoaCibNm0yfv7qq6/i6upKQUFBoX6rVq3C29sbHx8fAgICmDt3bmmXLsRTkQAWZZJSit69exMWFsalS5c4evQoMTExJCcnA1BQUMCmTZuoU6cO+/btM/bbvn078+bNIz4+nlOnTnHw4EFcXFy0Ogwh/pAEsCiTdu/ejY2NDf/4xz+M2+rWrcu4ceMA2Lt3L15eXowePZro6Ghjm1mzZjF37lxq164NgK2tLS+//HLpFi/EU5IAFmXSmTNnCAwMfOLn0dHRDBgwgMjISLZt22Z8gOL06dMEBQWVVplC/CUSwMIkjBkzBj8/P0JCQsjLy+Pbb7+ld+/eVKhQgebNm7Njxw6tSxSiyGQxHlEmeXl5sWHDBuPvCxcu5Pbt2wQHB7Njxw7u37+Pj48PAFlZWdjb29OjRw+8vLw4evSoLJojTIKcAYsyqX379uTk5LB48WLjtqysLODh5Ydly5aRlJREUlISly9f5rvvviMrK4spU6YwceJEbt68CUBeXh7Lli3T5BiE+DMSwKJM0ul0bN68mX379lG/fn2aNWvG0KFDmT59OnFxcXTv3t3Y1tHRkVatWrF161a6devG2LFj6dixI15eXgQGBvLgwQMNj0SIJ5O1IIQQooTJWhBCCFHGSAALIYRGJICF5uLi4vD09MTDw4PZs2c/8nlubi79+/fHw8OD5s2bk5SUZPxs1qxZeHh44OnpWWgq2vz58/H29sbLy4t58+YZt584cYLQ0FB8fHzo2bOnXB8W2nrci+Ke9CMv5RTFTa/XK3d3d3Xx4kWVm5urfH191ZkzZwq1WbhwoRo1apRSSqno6GjVr18/pZRSZ86cUb6+vionJ0ddunRJubu7K71er06dOqW8vLxUZmamys/PVx06dFCJiYlKKaWCg4PV3r17lVJKLV++XL399tuleLTCXPGEl3LKGbDQ1KFDh/Dw8MDd3R0bGxuioqKMb8v4TWxsLEOHDgWgT58+7Nq1C6UUsbGxREVFYWtrS/369fHw8ODQoUOcO3eO5s2b4+DggJWVFW3atGHjxo3Aw1cNhYWFAQ9fe/R/5xoLUdokgIWmUlJSqFOnjvF3Nzc3UlJSntjGysoKFxcX7ty588S+3t7eHDhwgDt37pCVlcW3337LtWvXAAq9Dunrr782bhdCCxLAotxp0qQJkyZNonPnzoSHh+Pv74+lpSUAX3zxBYsWLSIoKIj09HTju+OE0IIEsNCUq6trobPQ5ORk41uSH9dGr9eTlpZGlSpV/rDvSy+9xNGjR9m/fz+VKlWiUaNGADRu3Jj4+HiOHj3KgAEDaNCgQUkfohBPJAEsNBUSEkJiYiKXL18mLy+PmJgYIiIiCrWJiIgwvrJ+/fr1tG/fHp1OR0REBDExMeTm5nL58mUSExNp1qwZALdu3QLg6tWrbNy4kYEDBxbaXlBQwAcffFBouUshSpssxiM0ZWVlxYIFC+jSpQsGg4Hhw4fj5eXFO++8Q3BwMBEREbz00ksMHjwYDw8PKleuTExMDPDwem6/fv1o2rQpVlZWLFy40Hip4YUXXuDOnTtYW1uzcOFCKlasCDxcR2LhwoUAPP/887z44ouaHLcQII8iCyFEiZNHkYUQooyRABYm6/bt27w0chSHDh3SuhQhnokEsDBJZ8+exTcwmA0/JdK+UzjR0TFalyREkUkAC5OTn5/Pc61ac/v2HQryssDGgYEDB/DTTz9pXZoQRSIBLEyOtbU132yJZe7sGaRf+JkXunVg+/btBAQEaF2aEEUi09CESWrVqhV2dnYANGrUiPDwcI0rEqLo5AxYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoREJYCGE0IgEsDBZ2dnZANja2mpciRDPRgJYmKxz584BUKNGDY0rEeLZSAALk/WvadMB6Nmzp8aVCPFsJICFSdq/fz+3bl5nwOChVKxYUetyhHgmEsDC5Pzyyy+0adMGgLfe/KfG1Qjx7CSAhclIS0tj/qef0rhxYwCio2Pw9vYu8f0uXLKU+Z8tKPH9CPMjb0U2Ie++9wE3b6YW65hLFpdusIwaPbbIffL1+aTcuMmOb2KN22JjY4mIiCjO0p5o27btxH0TS/26fyu1fQrzoFNKPXXj4OBgdeTIkRIsRzyJXq+nUVMfrmcqHJqEFdu493YtLbax/ozO2o6KYYOL3O/3NTpVqIitvUNxlfWnHty/i2PzfuhPbOWH/Xvx9fUttX2L8kGn0x1VSgX/frucAZsIKysrdu34lsCQ5lhXdsPePahYxq0Q3KtYxilJDg1boPLz0NloM9/XTmeBlXNVMl1q0Klrd04eOyJT30SxkAA2IfXr12fr5o2E9+iFZZ8PsKn6N61LKhVWLmUj7BybtiHj/nU6d+vJzz/sx87OTuuShImTm3AmplWrViyY9zEZW2ZgyErTuhyz4xjan+R8BwYNGUZRLt8J8TgSwCZo2LChvDx0IBnfzEbp87Uux6zodBY4dhrPrp9P8u5772tdjjBxEsAmas6Hs2nRtD4ZuxfLmVgps7C2xannFD7+bBFr167TuhxhwiSATZSFhQXr135FtbybZB7eqHU5ZsfKqTLOPd/ipVH/QGYGiWclAWzCHB0d2Rm3DXVmO1kXftK6HLNjU6MBDu1fIbx7T5KTk7UuR5ggCWAT5+bmRtw3W8javYi81Etal2N2HBqFglc4nbr2IDMzU+tyhImRAC4HQkJCWPb5YtK3zsSQcU/rcsyOY8gL3LKqRt+oQRQUFGhdjjAhEsDlRFRUf14dM4qMb2ZSkJ+rdTlmRafT4dRhDD+ducSkKVO1LkeYEAngcuS9d6fRLsSHzJ2fycyIUqazssapx2SWrFjFypWrtC5HmAgJ4HJEp9OxZtWXuFpmkHkwRutyzI6lgwtOPacyZsJr/PDDD1qXI0yABHA5Y29vT/y3W7FI3Evmuf1al2N2bKrVxbHTeHr0iiQpKUnrckQZJwFcDtWsWZPvtm8je99Scq//onU5Zse+QTCWAZF06NKNBw8eaF2OKMMkgMspPz8/1qxcQfo3s9E/+FXrcsyOY2BP0irUJ7JPfwwGg9bliDJKArgc69WrF2+9+QYZW2dSkJetdTlmRafT4dh2JMeTbjH+tTe0LkeUURLA5dyUSW/SrW0omTvmoZTMUS1NOksrHLu9yep1m/h8yRKtyxFlkARwOafT6fhy+VLqOysyv1+tdTlmx9LOCeeIqfxz0lvs2bNH63JEGSMBbAZsbGyI+yYWm2uHyDi1U+tyzI51ZVccu75B7xf6kpiYqHU5ogyRADYTVatWZWfct+T9uIqc5DNal2N27Ov6YdN8AB06d+XePXlcXDwkAWxGmjZtyrroNWRsm0P+/Ztal2N2HP3CyarhS49ez5OfLwvpCwlgsxMeHs4H775DxpYZFOTK6l2lzTFsGOduZTFq9Bh5XFxIAJujCePH0qd7JzK2f4QqkDmqpUlnYYlT13+yYftOPpn3qdblCI1JABcjJycnTfb7zjvvsHPn099c0+l0fL5oAU1qOJK574tCnyUvHs715WO4/sVYUtf+y7i8pXH7inFcXzGOnORzAGRd+IkrH/Yg/8414xhKFXB35xKuL3+F68vHcGPla3LJ4/+wsHXgwY0kJv7zDerWrYufnx8fffTRMy9lOXPmzGKr7f79+yxatKjI/d59913mzp1bbHWYCwngcuC9996jY8eORepjbW3NttiNOP56hsyE7YU+qzFgJrWHL8CmpgdpB9cV3v7iZ9R+8TPs3JoAkHluP7ZuTck8+791J7LOHcCQfodawxdQ+6WFVIucioWdNn85lVU6K1uqDZjFvQcZLFiwgO3btzN9+vRnGutJAayUKnKoP2sAi2cjAVwC9u7dS9u2benTpw+NGzdm0KBBxut9hw8f5rnnnsPPz49mzZqRnp5OTk4OL774Ij4+PgQEBBjni3755Zf07t2bTp06Ua9ePRYsWMDHH39MQEAALVq04O7duwAMGzaM9evXA1CvXj2mTZtGYGAgPj4+nD9/HoBff/2VTp064eXlxYgRI6hbty56vZ5d8dvJP7yW7KSER47Dro43+feuP/E4C/KyyU0+Q5WuE8g8/78ANmTcxdKpMjrdwz9eVhWqYikB/Ag7t6bYthxK/0GDmTVrFgsWLEAphcFgYOLEiYSEhODr68uS//8hjhs3bhAWFoa/vz/e3t4cOHCAyZMnk52djb+/P4MGDSIpKQlPT0+GDBmCt7c3165dK/Qvs/Xr1zNs2DAAUlNTiYyMxM/PDz8/P3788UcmT57MxYsX8ff3Z+LEiQDMmTPHWMu0adOMY82YMYNGjRrRqlUrfvlF1hx5FlZaF1BeHT9+nDNnzlC7dm1atmzJDz/8QLNmzejfvz9r164lJCSEBw8eYG9vz/z589HpdJw6dYrz58/TuXNnLly4AMDp06c5fvw4OTk5eHh48OGHH3L8+HFee+01Vq1axauvvvrIvqtWrcqxY8dYtGgRc+fOZdmyZUyfPp327dszZcoU4uLiWL58OQAeHh5s3vA1PXq/gFXfwmdS2RcPYVOtnvH31Oi3wMICnaU1tYZ8TFbiQezqB2Fd2RVLO2dyb/4X25oeODRuzc01b5KTfAa7un44ebXDpkaDEvuuTZmjdwcy7qXwj7ETMBgM3Lp1i9jYWFxcXDh8+DC5ubm0bNmSzp07s3HjRrp06cLUqVMxGAxkZWXRunVrFixYQEJCAgBJSUkkJiaycuVKWrRo8Yf7Hj9+PG3atGHTpk0YDAYyMjKYPXs2p0+fNo4XHx9PYmIihw4dQilFREQE+/fvx9HRkZiYGBISEtDr9QQGBhIUFFTC31b5IwFcQpo1a4abmxsA/v7+JCUl4eLiQq1atQgJCQGgQoUKAHz//feMGzcOgMaNG1O3bl1jALdr1w5nZ2ecnZ1xcXGhZ8+eAPj4+HDy5MnH7vv5558HICgoiI0bNxr3sWnTJuDhTIhKlSoZ27dt25aP58zmn/96H1SBMWhtqtWjcuvBxnY1BszE0sHF+HvWuf04B0UA4NAkjKyz+7Ct6YFVhaq4vryEnCsnyLl6ktSYqVTtNRn7ev5/4Rstvxxb/Z2kbf8mKysLpRTx8fGcPHnS+K+atLQ0EhMTCQkJYfjw4eTn59O7d2/8/f0fO17dunX/NHwBdu/ezapVDxePt7S0xMXF5ZE5yvHx8cTHxxMQEABARkYGiYmJpKenExkZiYODAwARERHPevhmTQK4hNja2hr/29LSEr1e/5fHsbCwMP5uYWHxxDF/a1OU/Y58+WVOnj7LokWLqDn4Y6ycK/9he0N2OjlXTpL3axKgA1UA6KjYbjg6nQ6dlTX2DYKxbxCMhUNFshMPSgA/gU5ngW2zAdw79yPLvliBUorPPvuMLl26PNJ2//79bNu2jWHDhvH6668zZMiQR9o4Ojr+bnyd8b9zcnKKVJtSiilTpjBq1KhC2+fNm1ekccTjyTXgUuTp6cmNGzc4fPgwAOnp6ej1elq3bs2aNWsAuHDhAlevXsXT07NY992yZUvWrXt4Qy0+Pv6xT2PN/3gutlaWZP2w8k/nqGb98gOOXu1wG70Ct9Ff4PbKl1hVrEFu8hlyb/4Xffod4OGMiPxfL2NZoXqxHk95YshK496eZTgF9GT2R/OoXr06ixcvNj6sceHCBTIzM7ly5Qo1atTg5ZdfZsSIERw7dgx4eEP1jx7sqFGjBufOnaOgoMD4ryCADh06sHjx4oc1GAykpaXh7OxMenq6sU2XLl344osvyMjIACAlJYVbt24RFhbG5s2byc7OJj09na1btxb792IO5Ay4FNnY2LB27VrGjRtHdnY29vb27Ny5k1deeYXRo0fj4+ODlZUVX375ZaEz3+Iwbdo0BgwYwOrVqwkNDaVmzZo4OzsXamNpaUn16tWwyLhGxrEtOAb1euJ4mef24dK8T6FtDo2eI/PsPhwatuBu3Gcow8NQsKnViApBPYr1eEyd0udxfcU4MBjAwgJHr/ZUaNabvJttWPXVe0T1eZ7AwECUUlSrVo3Nmzezd+9e5syZg7W1NU5OTsbLByNHjsTX15fAwEBmzJjxyL5mz55Njx49qFatGsHBwcYwnT9/PiNHjmT58uVYWlqyePFiQkNDadmyJd7e3nTt2pU5c+Zw7tw5QkNDgYdTLf/zn/8QGBhI//798fPzo3r16sbLaqJodEV5Gic4OFgdOXKkBMsRJSU3NxdLS0usrKz46aefGD16tPFGy+9duXKFgOBmWLcdjX0D+R+rtGWeO4Du0GpOHDtCrVq1tC5HFAOdTndUKRX8++1yBmwmrl69Sr9+/SgoKMDGxoalS5c+sW3dunXZtmUzHcO7Y9nn/UIzIUTJc2zSmoz7KXTu1pOEIz9jaWmpdUmihEgAm4mGDRty/Pjxp24fGhrK87178m3CNmw6jSnBysTj6CrU4NeL+2S9iHJObsKJx9q0aRObvonDvkV/rUsxOznJ58j5fgW7dmzHyurJ50hxcXF4enri4eHB7NmzH/n8ypUrdOjQAV9fX9q2bUtycjIAe/bswd/f3/hjZ2fH5s2bS+pwxB+Qa8Dl0PDhw/nmm2+oXr06p0+fLnL/48eP07pdB5x7vYNtrYYlUKF4En3aLe6vncS6/3xJt27dntjOYDDQqFEjvvvuO9zc3AgJCSE6OpqmTZsa2/Tt25cePXowdOhQdu/ezYoVK1i9uvBbUe7evYuHhwfJycnGOb2i+D3pGrCcAZdDw4YNIy4u7pn63rhxg87demDfdpSEbykryM0ifcsHvDt18h+GL8ChQ4fw8PDA3d0dGxsboqKiiI2NLdTm7NmztG/fHnj4QM/vP4eHjyZ37dpVwlcjEsDlUFhYGJUr//GDFI+TnZ1Np649UJ4dcGzcqgQqE0+iCgxkxH1EZJd2vPH6a3/aPiUlhTp16hh/d3NzIyUlpVAbPz8/45OQmzZtIj09nTt37hRqExMTw4ABA4rhCMSzkAAWwMMnnvoPGswNVRHH5v20LsfsZB5YiWdVO5YuWVzoybW/Yu7cuezbt4+AgAD27duHq6troRkVN27c4NSpU4994k6UDpkFIQB4+51pHDh+ngovvF9sASCeTubJHdjfSGDbscNYW1s/VR9XV1euXfvfGszJycm4uroWalO7dm3jGXBGRgYbNmygYsWKxs/XrVtHZGTkU+9TFD85AxZ89VU0n36+DKcek9FZ2WhdjlnJuXqSvINfsSt+e5EuG4WEhJCYmMjly5fJy8sjJibmkQVxbt++bVwPeNasWQwfPrzQ59HR0XL5QWMSwGbu559/ZuQrY3Hu+RaWjpX+vIMoNvn3rpOx/SM2rIuhUaNGReprZWXFggUL6NKlC02aNKFfv354eXnxzjvvsGXLFuDhutSenp40atSI1NRUpk6dauyflJTEtWvXaNOmTbEekygamYZWDg0YMIC9e/dy+/ZtatSowfTp03nppZceaXf16lUCgpthGTYSB4/mGlRqvgw5GTxY+yaz3pnMK6NHa12OKGHyKLIZiY6O/tM2GRkZdAzvjoVPDwnfUqYMejK/ncPA53tJ+Jo5uQRhhgoKCnih3wDu2LviENxb63LMilKKzL3L8P1bFRZ8+onW5QiNyRmwGXpj4pscvnAN58h3ZcZDKcs6vo0K9xPZsuuQLLIjJIDNzfIvvmD5f9biEvVvdJYy/ag0ZV86iv7YBnYd/tn4Oiph3iSAzcj+/fsZ//pEKvadgaW9BEBpyrt9lcz4+cR9E0v9+vW1LkeUEXIN2ExcvHiRiMgXcOryKtZV6vx5B1FsDFlpZGyZwYJ5H9OqlTziLf5HAtgMpKWl0bFLN6yD+2JfP1DrcsyK0ueTse1DXh46kGHDhmpdjihjJIDLOb1eT8/IPjyo7IljQHetyzErSikyd39Oiyb1mPPho+v1CiHXgMu5MeMncDrlPs69xmlditnJOrKJKrnXWb/2Jyws5FxHPEoCuBxbsHAhMZu24dL/Q3QWMuWpNGUl/kzB6W/ZdfQwjo6OWpcjyigJ4HIqPj6eyW9Pw6XfLCzsnLQux6zkpV4ia9cC9ny3Azc3N63LEWWY/LuoHDp//jx9+g/Eqes/sa5UW+tyzIoh4x7pW2ey7PPFhISEaF2OKOMkgMuZO3fu0DG8G7ahg7D7m4/W5ZiVgvxcMr6ZyatjRhEVJS8zFX9OArgcycvLo1tEJDm1A3H07ax1OWZFKUXmzs9oG+zNe+9O07ocYSLkGnA5oZTipZH/IPFuPs49h2hdjtnJPLgWV8sMvlr9rayvIZ6aBHA5MWfuR2zdeYAKfWfKjIdSlnnuABaJe4g/dgR7e3utyxEmRAK4HNi6dSvvzfqQiv3/jYWtvF68NOVe/4Xsff8fP+zbQ82aNbUuR5gYuQZs4k6ePMnAwUNx7j4ZK5fqWpdjVvQPfiVj24esWbkCPz8/rcsRJkgC2ISlpqbSqWt37MJewta1sdblmJWCvBwyts5kysTX6dWrl9blCBMlAWyicnJy6NI9AoN7axybttW6HLOiVAGZOz6hW9tQpkx6U+tyhAmTADZBSin+PvRFruXa4ficvFa8tGV+/x/qOyu+XL5UZjyIv0Ruwpmg6e9/wM6DCVTo8wE6nfwdWpoyT+/C5trPxB07go2NjdblCBMnAWxi1q37mrnzF1Ix6t9YWNtpXY5ZyUk+Q+4PK9n3wwGqVq2qdTmiHJAANiFHjhxh+MhROPeehpVzFa3LeWrJi4ZhSL9dovuwsSvB+bc6Syp0eoXs/V+wPnoNTZs2Lbl9CbMiAWwi9Ho93Xr2wmBfmfxze8g/t6dYxr1zKLZYxnlaLl5hRe6jz85An5cLFhYYUi+iz83ivfffZ0BUFAC2trZUqlSpuEs1iuw3gN3bPmbu3I8IDw8vsf0I8yMBbCKsrKyY99Ecfv3112Id99VSDGBHR0emv/x8kfvp9Xpupt7iP9Frya1RB321Rvx73kIsLKyY+tbkEqi0sOCgQNzr1eXVCbKovSheOqXUUzcODg5WR44cKcFyhHiygoICduzYwYh/vEJW1SYYrp3ko5nv8vKIEVqXJsQf0ul0R5VSwb/fLrfQhcmwsLCga9euHD/8MxXuJWJRL4RXX3uDBw8eaF2aEM9EAliYnOrVq7Ntyyb0Fw5gV8eL1atXa12SEM9EAliYJG9vb7y9vcizcWHu/M8oyqU0IcoKCWBhssaNHomdPp3Um6ncvHlT63KEKDIJYGGy6tWrB9lp2Dq5cO/ePa3LEaLIJICFybK1tUUV6LGwsiE/P1/rcoQoMglgIYTQiASwEEJoRAJYCCE0IgEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwMIkrVmzhnXr1hl/j4mJITo6WsOKhCg6eS29MDn5+fnUdffg1u07VO0yBpSB2zs+p1bNGiRdvIClpaXWJQpRiLyWXpQb1tbWHD9yCK+mTcm/eJD8iz/j5+fLscMHJXyFSZEAFiapRo0a/PzDfroHNySiRVN+3L+HatWqaV2WEEUiASxMyubNm9HpdJw/fx47Oztmz3iP9Pt3sLW11bo0IYpMAliYlOjoaFq1amW84Va7dm3Wr1//SDu9Xl/apQlRZBLAwmRkZGTw/fffs3z5cmJiYgBISkrC29sbgC+//JKIiAjat29Phw4dtCxViKdipXUBQjyt2NhYwsPDadSoEVWqVOHo0aNUqVKlUJtjx45x8uRJKleurFGVQjw9OQMWJiM6OpqoqCgAoqKiHjvvt1OnThK+wmTIGbAwCXfv3mX37t2cOnUKnU6HwWBAp9MxZsyYQu0cHR01qlCIopMzYGES1q9fz+DBg7ly5QpJSUlcu3aN+vXrc+3aNa1LE+KZSQALkxAdHU1kZGShbS+88AKzZs3SqCIh/jp5FFkIIUqYPIoshBBljASwEEJoRAJYCCE0IgEsSpWlpSX+/v7Gn9mzZz913+vXr9OnT58i73PEiBGcPXu2yP2EKGlyE06UKicnJzIyMrQuQ4hSJTfhRJlWr149pkyZgr+/P8HBwRw7dowuXbrQoEEDPv/8c6Dwug9nzpyhWbNm+Pv74+vrS2JiIpmZmXTv3h0/Pz+8vb1Zu3YtAG3btuW3E4fo6Gh8fHzw9vZm0qRJxv07OTkxdepU/Pz8aNGiBampqaX8DQhzJAEsSlV2dnahSxC/hSTA3/72NxISEmjdujXDhg1j/fr1HDx4kGnTpj0yzueff86ECRNISEjgyJEjuLm5ERcXR+3atTlx4gSnT58mPDy8UJ/r168zadIkdu/eTUJCAocPH2bz5s0AZGZm0qJFC06cOEFYWBhLly4t0e9BCJBHkUUps7e3JyEh4bGfRUREAODj40NGRgbOzs44Oztja2vL/fv3C7UNDQ1lxowZJCcn8/zzz9OwYUN8fHx44403mDRpEj169KB169aF+hw+fJi2bdsaF24fNGgQ+/fvp3fv3tjY2NCjRw8AgoKC+O6774r3wIV4DDkDFmXGb4uqW1hYFFpg3cLC4pH1fQcOHMiWLVuwt7enW7du7N69m0aNGnHs2DF8fHx4++23ee+9955639bW1uh0OuDhjUJZT1iUBglgYZIuXbqEu7s748ePp1evXpw8eZLr16/j4ODA3//+dyZOnMixY8cK9WnWrBn79u3j9u3bGAwGoqOjadOmjUZHIIRcghCl7LdrwL8JDw8v0lS036xbt47Vq1djbW1NzZo1eeuttzh8+DATJ07EwsICa2trFi9eXKhPrVq1mD17Nu3atUMpRffu3enVq9dfPSQhnplMQxNCiBIm09CEEKKMkQAWQgiNSAALIYRGJICF5uLi4vD09MTDw+OxN+Ryc3Pp378/Hh4eNG/enKSkJAC+++47goKC8PHxISgoiN27dxv7HD16FB8fHzw8PBg/fjy/3ev4+uuv8fLywsLCArmfIbQmASw0ZTAYGDNmDNu3b+fs2bNER0c/snDO8uXLqVSpEv/973957bXXjI8QV61ala1bt3Lq1ClWrlzJ4MGDjX1Gjx7N0qVLSUxMJDExkbi4OAC8vb3ZuHEjYWFhpXeQQjyBBLDQ1KFDh/Dw8MDd3R0bGxuioqKIjY0t1CY2NpahQ4cC0KdPH3bt2oVSioCAAGrXrg2Al5cX2dnZ5ObmcuPGDR48eECLFi3Q6XQMGTLE+MhxkyZN8PT0LNVjFOJJJICFplJSUqhTp47xdzc3N1JSUp7YxsrKChcXF+7cuVOozYYNGwgMDMTW1paUlBTc3Nz+cEwhygJ5EEOYvDNnzjBp0iTi4+O1LkWIIpEzYKEpV1fXQq+WT05OxtXV9Ylt9Ho9aWlpVKlSxdg+MjKSVatW0aBBA2P75OTkPxxTiLJAAlhoKiQkhMTERC5fvkxeXh4xMTHGVdF+ExERwcqVKwFYv3497du3R6fTcf/+fbp3787s2bNp2bKlsX2tWrWoUKECBw8eRCnFqlWr5JFjUTYppZ76JygoSAlR3LZt26YaNmyo3N3d1QcffKCUUupf//qXio2NVUoplZ2drfr06aMaNGigQkJC1MWLF5VSSr3//vvKwcFB+fn5GX9SU1OVUkodPnxYeXl5KXd3dzVmzBhVUFCglFJq48aNytXVVdnY2Kjq1aurzp07a3DEwtwAR9RjMlXWghBCiBIma0EIIUQZIwEshBAakQAWQgiNSAALIYRGJICFEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiMSwEIIoZEiLUep0+l+Ba6UXDlCCFEu1VVKVfv9xiIFsBBCiOIjlyCEEEIjEsBCCKERCWAhhNCIBLAQQmhEAlgIITQiASyEEBqRABZCCI1IAAshhEYkgIUQQiP/D3KCxPtv1yVyAAAAAElFTkSuQmCC\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": "iVBORw0KGgoAAAANSUhEUgAAAq8AAAExCAYAAABS5UszAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABK+0lEQVR4nO3dd3gU1dvG8e9JIyEJvRcFpEgPHUQQkF4iCAiIUhRRFOwIWLD9KCp2BLGDhSIovUuxIk2kKoh0EKQnIT3n/WNDXoK0QJLJZO/PdeXS7M7O3LsJT56ZOTPHWGsREREREXEDH6cDiIiIiIhcKTWvIiIiIuIaal5FRERExDXUvIqIiIiIa6h5FRERERHXUPMqIiIiIq6h5jUDGGOsMaas0zmuhDHmOmNMpDHG1+ksIiJOUM0WcReval6T/8Gf/UoyxkSf833Pi7ymiTFmfzpmWGGM6Zde67tW1tq91toQa21ieq/bGPOZMSYu+fM9boxZYoy5Mfm5F4wx8ef9TJ4657XtjTGrjTFRxphjxpgvjTElLrCNF5L/8NRLY7bixphZybn2G2MeuMzydxpj9iTnmWmMyXfOc/mMMd8mP7fHGHPnOc/1McYknvc+I40xxdKY91LbaJL8+3zu+nunZf0iWZFq9n+5tWYbYwKMMa8n19tIY8xuY8xbacj2sjFmkzEmwRjzwmWWvVDWMuc8H2aMWWeMOZP837CLfAZnv36/0pxXso3k52saY75PXv9hY8wjad2GN/Oq5jX5H3yItTYE2At0OOexL53Ol029mvx5lwCOAJ+d89zUc38m1tpXAYwxXYCvgLeAAkBlIBb40RiT9+yLjTEG6AUcT/5vWnwB7AIKA+2AkcaYphda0BhTGZgA3J28/Blg3DmLvAfEJT/XExif/JqzfjnvfYZYaw9eYDsvXKIoX24bB89b/8TLfQAiWZ1qtiMyqmYPA2oDdYFQoAmwPg25/gKeAuZd4fLnZ/07OWsAMAvP34C8wERgVvLjZ7163murX2gDyTs2TS7w+CW3YYwpACzE83clP1AWWHyF70vwsub1YowxOYwxbxljDiZ/vZX8WDCwACh27hEzY0xdY8wvxpiTxphDxpix5/3iX+l2myTvhT5hjDmSvK6+5zwflLynuscYc8oY86MxJij5uXBjzJbkDCuMMRXPed1uY8xgY8zG5L3gj40xhY0xC4wxEcaYpWcLijGmlPEcufRL/n5F8h7uT8nLLk7+h3Z23b2S8xwzxjyXvK3ml3uv1tozeIpblct8JgZ4HfiftfYra220tfYfoB8QCTx2zuKNgKLAw0D3K/0ZGGNC8BTOEdbaeGvt78B04J6LvKQnMMda+721NhJ4DrjdGBOa/DvSGXjOWhtprf0RmI2n0U0X6bmN83/eyY+lHFkyxvgm/84dNcbsMsYMPH95EaepZruyZtcBvrXWHrQeu621k67og/fkmWitXQBEXOlrLqIJ4Ae8Za2Ntda+Axig2TWuNy3beBxYZK39Mvn5CGvttgutyFzgTMK5P8Pk37mJxpgTxphtxpinzl8+O1Lz6vEMUB8IA6rj2TN81lobBbQh9VGtg0Ainn+QBYAGwK3Ag1e57SJAbqA4cC/wnvn/PdUxQC3gJiAfnr3OJGNMeWAy8ChQEJgPzDmvGHcGWgDlgQ54CvrTycv74Gn4LuZOoC9QCAgAngQwxlTCc8SxJ56m8Wzuy0puGHsCv11m0QrAdcDX5z5orU0CZiS/p7N6A3OAacnfdzhne3caYzZeLM55/z37/xcr0pWBlNNG1tqdeI6Clk/+SrDWbj9n+d+TX5NermQbhYzn1NMuY8ybyX/Er8Z9eH7nw4CaQMerXI9IRlLNTs0NNXsV8Lgx5kFjTNXkpvfc7Y0zxpx7RutadTCeoQ9bjDEDznm8MrDRWmvPeWwj6VuzL7eN+sBxY8zPyTtBc4wx113ltp4HSgFl8HzWd13lelxFzatHT+Ala+0Ra+2/wItc4qiWtXadtXaVtTbBWrsbz6H/W65y2/HJ24631s7Hs6dawRjjg+dI4CPW2gPW2kRr7c/W2ligGzDPWrvEWhuPp2AG4SmYZ71rrT1srT0A/AD8aq39zVobA3wL1LhEpk+ttduttdF4GsOw5Me74DkC+aO1Ng4YDtiLrOOsJ40xJ/Gc8gkB+pzz3B3JRyHOfhXD88cF4NAF1nXo7PPGmJxAV+Cr5M9gOucMHUg+AlDtQoGstRHAT8BzxphAY0xNPH84cl7kPYQAp8577BSeU18hwOmLPHdW/fPe586LbOdiLreNP/D8jIri2bOvBbyRxm2cdQfwtrV2v7X2BDD6KtcjkpFUs1PL8jUbGAW8gudntxY4YM4Zm2+tfdBae7U7FOebBlTE0/jfBww3xvRIfu5S9fysJ897n2kdhnW5bZTAc/DlETyN/y48OzdX4w5gpLX2hLV2P/DOVa7HVdS8ehQD9pzz/Z7kxy7IGFPeGDPXGPOPMeY0MJL//weaVsestQnnfH8Gzy9+ASAQuFCjkypv8h7uPlLvUR8+5/+jL/B9yCUy/XOBPGe3u++c7Z4Bjl1iPQBjrLV5rLVFrLXhyUctz5qW/NzZr4PA0eTnil5gXUXPeb4TkIDnCAbAl0AbY0zBy+Q5qydQOvn9jMczNulip1oigVznPZYLz+mrSz131qrz3ucNZ59I/j06mfzHYigw9JyCOfcKto+19h9r7VZrbZK1dheeoz2dL/P+LybVz/i8/xfJKlSzU8vyNTu5mX/PWtsQyAOMAD45d/hEekmuhwfP7kAAb+Np5OHKavaY895nSpN9blML3AzMPeexoVe4jWg8QyjWJO+cvAjcZIzJfRVv1ytrtppXj4PA9ed8f13yY3DhvdTxeI52lbPW5sJzasdcYLlrcRSIAW64wHOp8iaffikJHEjnDOc7hGeP8ex2g/AMNk9Pf+JpIrue+2DyUY3OwHfJD/XGU6D3GmP+wXPKyh/P6bPLstbusda2t9YWtNbWw/OHZ/VFFt+C59Tk2SxlgBzA9uQvP2NMuXOWr578mivJ0f5sgcRzlHP0OQWzffJiad2G5fL/toPO+f9zC2aqnzGe3yuRrEY1+8pkpZqdwnrGxb4HnAAqpXOeC7H8/897C1DtvGEL1bjymp3nnJr9I9D+nMfOnqm63DY2kvr39HJHw4PO+97ra7aaV4/JwLPGmILJA92H4zkSB5693/zn7RGF4jmNG2k8txEZQDpL3jP/BHjDeC448DXGNDDG5MBzSqSdMeZWY4w/8ASeKzt/Tu8c55mOZxzRTcljtV4gnf8AJI8RehLPz+PO5NP6RYCP8Oy5vmmMKY5nzFp7PKfHwvA0c69whXcdMMZUNJ4LrgKMMXcBLbn4qfYv8bzvRsljSV8CvrGeQfZRwDfAS8aYYGNMQ+A24POref8XcrltGGOaGmOuNx4l8TTBsy6z2nuMMT7GmBZ4impo8u/SNOAR47mVWB5gSHq9D5F0pJp9ZbJEzQYwxjxqPBcfBRlj/JKHDIRy+TG1JL/e3xgTiKdv8UvezgXvdWuMuc0Ykze5JtbFM174bE1cgWcM9MPGc5HfwOTHl13N+7+Iy23jU6CT8dxOyx/PRcA/WmvPH2pwVj5jTMfkmn0PniPXuZPf/zRgWPL7LQ4MvMg6shU1rx7/wzMGZyOwCc/tO/4HYK39A0+h/Nv8/xifJ/Ec4YsAPgSmZlCuJ5PzrMFzO6hXAB9r7Z94BmW/i2dvvwOeW8jEZVAOAKy1W4BBwBQ8e3uReG6lEpvO25mKZ/zaY3hOcW3Fs+fZ0Fp7LPm5DdbaxcmnzP+xnqtb38Gzt1vFGNPTGHOpPelWwN949vwfAFonj50DUu4v2eic9/0Anib2CJ6Ce+7YrAeT8x3B87syIPk1ZzUw/73Pa500fiyX2kYNPH8Eo5L/u4lLX9wBnqND/+JpTl/E88e/PJ7f58V4/i38hmdYRgKeQiySVahmX4EsVLPBM5zhdTxDHI4CDwGd7f/fwup9Y8z7l9jMh3hOt/fAc8FedPI2ST6wEHnOst3xjNmNACYBr9jk2wcmf+Yd8RzoOIlnnHLH834WT51Xr4+SBpfbhrV2GZ6j//Pw/DzKcumzhv/gGdv6L56/Xe8A7+M5AvsSniPfu4CleHZY0vXnmxUZay93tFrkwoznatSTeE7F7XI4jlwBY0wpPEXO/7xxexdbvg3wvrX2+sstKyJZm2q2+xjPfWS/sNb+Z5Keiyw/AOhurb3aCxJdQUdeJU2MMR2MMTmTT5+PwXOUYbezqSS9JJ/Sa5t8Wq84ntuwfOt0LhG5OqrZ2ZsxpqgxpmHykIIKeIakZPuareZV0uo2PBcfHATK4dnD0+H77MPgGUZwAs+wgW14hhSIiDupZmdvAXhu/RaBZ0ztLFLPAJktadiAiIiIiLiGjryKiIiIiGukqXk1xizMqCAiIhnJG+uXN75nEckeLlW//NKyoly5crWqXbu2xhmIiBudP8VutqeaLSIudtGanabmtVy5cqxdu/ba44iIZDJjzA6nM2Q21WwRcatL1WyNeRURERER11DzKiIiIiKuoeZVRERERFxDzauIiIiIuIaaVxERERFxDTWvIiIiIuIaal5FRERExDXUvIqIiIiIa6h5FRERERHXUPMqIiIiIq6h5lVEREREXEPNq4iIiIi4hppXEREREXENNa8iIiIi4hpqXkVERETENdS8ioiIiIhrqHkVEREREddQ8yoiIiIirqHmVTKEr68vYWFhVKlShQ4dOnDy5Ml0Xf/IkSNTfX/TTTdd1Xp2797NV199lebX9enTh+nTp1/VNkVEshrVbHETNa+SIYKCgtiwYQObN28mX758vPfee+m6/vML4c8//3xV67naQigikp2oZoubqHmVDNegQQMOHDgAwM6dO2ndujW1atWiUaNG/PHHHwDMmTOHevXqUaNGDZo3b87hw4cBiIyMpG/fvlStWpVq1aoxY8YMhg4dSnR0NGFhYfTs2ROAkJAQALp37868efNStn12b3v37t00atSImjVrUrNmzZTCOXToUH744QfCwsJ48803SUxMZPDgwdSpU4dq1aoxYcIEAKy1DBw4kAoVKtC8eXOOHDmSOR+eiEgmU82WLM9ae8VftWrVsiJXIjg42FprbUJCgu3SpYtdsGCBtdbaZs2a2e3bt1trrV21apVt2rSptdba48eP26SkJGuttR9++KF9/PHHrbXWPvXUU/aRRx5JWe/x48dTrf/87X3zzTe2V69e1lprY2NjbYkSJeyZM2dsVFSUjY6OttZau337dnv2d3n58uW2Xbt2KeuZMGGCffnll6211sbExNhatWrZv//+286YMcM2b97cJiQk2AMHDtjcuXPbr7/+Oj0+KskkwFqbhnqXHb5Us+VKqWZLVnOpmu3ndPMs2dPZvewDBw5QsWJFWrRoQWRkJD///DNdu3ZNWS42NhaA/fv3061bNw4dOkRcXBylS5cGYOnSpUyZMiVl+bx5815yu23atOGRRx4hNjaWhQsX0rhxY4KCgjh16hQDBw5kw4YN+Pr6sn379gu+fvHixWzcuDFlbNSpU6fYsWMH33//PT169MDX15dixYrRrFmza/p8RESyEtVscRMNG5AMcXb81J49e7DW8t5775GUlESePHnYsGFDyte2bdsAGDRoEAMHDmTTpk1MmDCBmJiYq9puYGAgTZo0YdGiRUydOpVu3boB8Oabb1K4cGF+//131q5dS1xc3AVfb63l3XffTcm3a9cuWrZseXUfgoiIS6hmi5uoeZUMlTNnTt555x1ef/11cubMSenSpfn6668BT9H5/fffAc/ecvHixQGYOHFiyutbtGiR6sKBEydOAODv7098fPwFt9mtWzc+/fRTfvjhB1q3bp2y/qJFi+Lj48Pnn39OYmIiAKGhoURERKS8tlWrVowfPz5l3du3bycqKorGjRszdepUEhMTOXToEMuXL0+Xz0dEJCtRzRY3UPMqGa5GjRpUq1aNyZMn8+WXX/Lxxx9TvXp1KleuzKxZswB44YUX6Nq1K7Vq1aJAgQIpr3322Wc5ceIEVapUoXr16ikFqH///lSrVi1l8P+5WrZsycqVK2nevDkBAQEAPPjgg0ycOJHq1avzxx9/EBwcDEC1atXw9fWlevXqvPnmm/Tr149KlSpRs2ZNqlSpwv33309CQgKdOnWiXLlyVKpUiV69etGgQYOM/thERByhmi1ZnfGMib0ytWvXtmvXrs3AOCIiGcMYs85aW9vpHJlJNVtE3OpSNVtHXkVERETENdS8ioiIiIhrqHmVLO/w4cPceeedlClThlq1atGgQQO+/fbblOcfffRRihcvTlJSUqrXTZo0iSpVqlC1alVq1KjBmDFjMju6iIhXUt2WjKTmVbI0ay0dO3akcePG/P3336xbt44pU6awf/9+AJKSkvj2228pWbIkK1euTHndggULeOutt1i8eDGbNm1i1apV5M6d26m3ISLiNVS3JaOpeZUsbdmyZQQEBPDAAw+kPHb99dczaNAgAFasWEHlypUZMGAAkydPTllm1KhRjBkzhmLFigGQI0cO7rvvvswNLyLihVS3JaOpeZUsbcuWLdSsWfOiz0+ePJkePXrQqVMn5s2bl3Kvv82bN1OrVq3MiikiIslUtyWjqXkVV3nooYeoXr06derUIS4ujvnz59OxY0dy5cpFvXr1WLRokdMRRUTkHKrbkt78nA4gcimVK1dmxowZKd+/9957HD16lNq1a7No0SJOnjxJ1apVAThz5gxBQUG0b9+eypUrs27dOs1nLSKSyVS3JaPpyKtkac2aNSMmJobx48enPHbmzBnAc+rpo48+Yvfu3ezevZtdu3axZMkSzpw5w7Bhwxg8eDD//PMPAHFxcXz00UeOvAcREW+iui0ZTc2rZGnGGGbOnMnKlSspXbo0devWpXfv3rz44ossXLiQdu3apSwbHBzMzTffzJw5c2jbti0DBw6kefPmVK5cmZo1a3L69GkH34mIiHdQ3ZaMpulhRcQraHpYERH30PSwIiIiIpItqHkVEREREddQ8yoiIiIirqHmVbK8hQsXUqFCBcqWLcvo0aP/83xsbCzdunWjbNmy1KtXj927dwNw7NgxmjZtSkhICAMHDszk1CIi3kk1WzKamlfJ0hITE3nooYdYsGABW7duZfLkyWzdujXVMh9//DF58+blr7/+4rHHHmPIkCEABAYG8vLLLzNmzBgnoouIeB3VbMkMal4lS1u9ejVly5alTJkyBAQE0L17d2bNmpVqmVmzZtG7d28AunTpwnfffYe1NuUWLIGBgU5EFxHxOqrZkhnUvEqWduDAAUqWLJnyfYkSJThw4MBFl/Hz8yN37twcO3YsU3OKiIhqtmQONa8iIiIi4hpqXiVLK168OPv27Uv5fv/+/RQvXvyiyyQkJHDq1Cny58+fqTlFREQ1WzKHmlfJ0urUqcOOHTvYtWsXcXFxTJkyhfDw8FTLhIeHM3HiRACmT59Os2bNMMY4EVdExKupZktm8HM6gMil+Pn5MXbsWFq1akViYiL33HMPlStXZvjw4dSuXZvw8HDuvfde7r77bsqWLUu+fPmYMmVKyutLlSrF6dOniYuLY+bMmSxevJhKlSo5+I5ERLIv1WzJDMZae8ULa55sEXGrS82TnV2pZouIW12qZmvYgIiIiIi4hppXEREREXENNa8iIiIi4hpqXkVERETENdS8ioiIiIhrqHkV1ztw4ADht3dlwYIFTkcREZHLiImJ4f4HB/HKa2NIyx2PRM5S8yqutnbtWqrXrM0P++Poemcv3nzrHRVDEZEs6vDhw9S/+Ra+/v53Rr3zIXfe3ZvY2FinY4nLqHkV15o+fQY33dyI+FI3EVLvDoIa9eGpYcPo1/8BkpKSnI4nIiLn2LJlC9Vq1OKvE4kE3dyboCb9mbloOQ1vacbx48edjicuohm2xLX+3Pk35StVIeLoJvaO/Ybc+QpQoWIl9h04qKkGRUSymAMHD1GwSFHyJyaxbUI//HMEUb5SZeLj44mLi3M6nriIZtgS11u4cCFt2rTh7j73MunTj5yOI1mUZtgSyRoiIyMJDQ2lZJly7N253ek4kkVphi0RERERyRbUvIqIiIiIa6h5FRERERHXUPMqIiIiIq6h5lVEREREXEPNq4iIiIi4hppXEREREXENNa8iIiIi4hpqXkVERETENdS8ioiIiIhrqHkVEREREddQ8yoiIiIirqHmVURERERcQ82riIiIiLiGmlcRERERcQ01ryIiIiLiGmpeRURERMQ11LyKiIiIiGuoeRURERER11DzKiIiIiKuoeZVXC8+Ph4APz9fh5OIiMjlnK3Zvr5+DicRt1LzKq534sQJAAoVyO9wEhERuZyTJ08CkDtPbmeDiGupeRXXm7tgMQAVypdzOImIiFzOkiVLAahYobzDScSt1LyKqx09epSvp3wJQNeuXR1OIyIil2Kt5bU33wbgkYcGOJxG3ErNq7jaRx9/AkD/AQ8REhLicBoREbmUX3/9lb/+2ELZGytTr149p+OIS6l5FdfatGkTw4YOAeDRQQ85nEZERC7l1KlT9Li7NwCDH3sEY4zDicSt1LyKK/3+++80aNgIgKeGPk3FihUdTiQiIhdz7NgxmjRvxe6/ttPwlmb07dvH6UjiYmpexVX279/PM88+R1hYGFERp7ird19Gj/yf07EYPOxZ7rnvfqy1TkcREckyIiIiGDduPDdWrsqGtb9SqWoY82d/i7+/v6O5Zs+eTf2bb+HMmTOO5pCrk21vsrZv3z6WL1/udAxJBwkJCZw8eZLFy1ayaN7slMefeW44L7/4QpY49bR16zYWzp9DxQrlGfzkE07HEXGd2NhYpk6d6nQMSQdJSUmcPn2aDRs38+nHH6Y83qptB2ZMm0xwcLCD6Tz27NnD6lW/0LVHT+Z8OwMfHx3LcxOTliNFtWvXtmvXrs3AOOnnjTff4sknniB/tVswuhGyq/27fsl/Hrv33nu5/vrrHUhzYZ9/OZl/i9Yj/ve5TJn0Ke3bt3c6kpzHGLPOWlvb6RyZyU01e926ddSuXZs85WrhH5rP6ThyDS5Uszt06ECdOnUcSHNhq1atYtVhS9Lx/dx3R3teHT3S6UhynkvV7Gzb1T08aCCz585nc0Qgwc3udzqOXIOQRF+itq3EL1dBct5QG+Prz4xtkbBti9PRUtgCYeSq3pr4EpXpcVcvfv5hJVWrVnU6lohr1KpVi+eef4F3P5tKYJOn8fHP4XQkuUr5i1bl+MJ38QnOQ86ydfENDOGHY/DDwqxTsyGUoOqN8A3Oy/ufDKZKpYr06nW306HkCmXbI6/gubIxrHY9Iks3JbimjoRJ5ojcshy/36axcf1aChUq5HQcSaYjr1mftZbO3Xqw8o/DhLR5IksMCZLsL+7f3ZyeMZzF8+dw0003OR1Hkl2qZmfrQR65c+fmu0XziV83g+i/1zkdR7xESOWmJJS+mVbtwomJiXE6johrGGP4atJnFPeLJOqXyU7HES8RULAUwS0G0f62TuzZs8fpOHIFsnXzClCmTBnmzJxB1OK3iTu61+k44iWCb+rB3pgc3N3nHt2BQCQNAgMDWTx/Dj5/rSRq60qn44iXCLqhDr5ht9GsZRsiIiKcjiOXke2bV4BGjRrx7puvEzl7BIlnTjkdR7yAMT4Et3yYpb/8xsv/G+F0HBFXKVKkCEsXzifm+4+IPfin03HES+SsGc7J0NJ07HwHiYmJTseRS/CK5hWgb98+9OvVg6h5r2IT452OI17Axz+Q4PbDePWtsXz99XSn44i4SrVq1fhy0mdEzB1Nwul/nY4jXsAYQ0jT/vy26zCPPP6k03HkErymeQUY8+or1L3xOqK+e1+nciVT+IXmJ7TDMPre15/169c7HUfEVcLDw3lmyJNEzhlBUly003HECxhfP4LbPsWkqTOY8MEHTseRi/Cq5tXHx4fpU78if+xBzqz91uk44iVyFClLzqYDaNm2PQcPHnQ6joirDH1qMO2aNiRq0ZtYm+R0HPECvkGhhIY/yxNPDWPFihVOx5EL8KrmFSAkJISlC+eRtGkeZ3b86nQc8RI5K9wEFVvSvHU7TUcokgbGGD796APK5DJE/TDJ6TjiJfzzFSe49eN07NyVv/76y+k4ch6va14BSpYsyYK5sznz3VjijvztdBzxEsF1u3DYpwB33HkXSUk6giRypQICAlgwZyY59q8lctNSp+OIlwgqFYZ/3e7c2rINJ06ccDqOnMMrm1eAunXr8uH4cUTMHklipH4pJeMZYwhp/hA/bdzBsGeeczqOiKsUKFCApYvmE/fzJGL2bXY6jniJ4OqtiSpUhfYdOxMfr4u9swqvbV4BevToziMP9idy7iiS4mOdjiNewPj5E9JuKOM++owvvvjS6TgirlKxYkW+nvIVkfPHEH/yH6fjiJcIbtyXbYejeODBgbrYO4vw6uYV4OUXX6BJ7cpELR2rX0rJFL7BeQgNf5oHBj7ML7/84nQcEVdp1aoVI158nsjZ/yMpNsrpOOIFjI8vIW2eZPq8Jbz9zlin4whqXj3TEX4+kWI+p4laNdXpOOIlAgqWImeLQbQL76jpCEXS6JGHB9K1fSsi54/BJulm8pLxfHLkJCT8GZ594SUWLlzodByv5/XNK0BQUJBnOsIdy4na9oPTccRL5LyhDr5h4dzaqq2mIxRJo/fHjaVS0VAiV37sdBTxEv55ihDSbjB39OjJ1q1bnY7j1dS8JitatCiL588leuUHxB7a7nQc8RI5a97GyZDr6dSlm6YjFEkDPz8/5s6cQejRbURtmO90HPESgSUqk6Nhb5q3bsvRo0edjuO11LyeIywsjC8++yR5OkL9UkrGM8YQ3PR+1u08xGNPPuV0HBFXyZMnD0sXzSd+zTSid29wOo54ieAqtxJXsh6t299GbKwu9naCmtfzdOzYkWFPPkbk3JEkxcU4HUe8gPH1I6TdED77ahoffPih03FEXKVs2bLM+mY6kQvfIP7YfqfjiJcIvvkudkX60Ofe+3SxtwPUvF7A00OH0PaW+pqOUDLN2ekIHx88lJUrVzodR8RVbrnlFt58bTQRs/9HYvRpp+OIFzDGh+CWj7Jg5SpGvfKq03G8jprXCzg7HWHpUEvUj184HUe8hH/+EgS3fpzbbu+i6QhF0qj/fffRu0dXIue9ik3UzeQl4/kEBBLS4RlGvvYGs2bNcjqOV1HzehE5cuRg4dxZBOz7lajN3zkdR7xEUKkw/Ot049ZWbTl58qTTcURc5a3XX6PWDUWJWv6BTuVKpvDLVYDQ9kPp2bsvGzZscDqO11DzegkFChRg6cL5xP40kZj9W5yOI14iOKwNUQUq0aFTFxISEpyOI+Iavr6+fDt9Knmj9nJmvY6ESebIUbQ8Qbf0p2Xb9vzzj2Z+ywxqXi+jUqVKTP3qCyLnvabpCCXTBN9yD1sOnmbAQ4OcjiLiKqGhoSxdOI/EDbM5s3ON03HESwRXbERSuaa0aNOe6Ohop+Nke2per0CbNm14+fnniJwzgqTYM07HES9gfHwJafsk0+Ys0nSEIml0/fXXM2/2TKIWv0Pcv7udjiNeIrh+Nw4m5aLH3b01bCWDqXm9Qo8+MogubVsQuUDTEUrm8MkRTEj4Mzzz/IssWrTI6TgirtKgQQPeH/sOEbNHkhh10uk44gWMMQQ3H8jKtVsY/sKLTsfJ1tS8XiFjDO+PG0vFwsFEff+Z03HES/jnKUJI2yfp2v1Otm3b5nQcEVe5++67eLBfHyLnjcYmxDkdR7yAj38OQtoP4633JjBlylSn42Rbal7TwN/fn7kzZ5DzyEaifl/odBzxEoElqxBwUy+at2rLsWPHnI4j4iqjRrxMw2rliFz6nk7lSqbwDclLaPjT9HtgAKtXr3Y6Trak5jWN8ubNy3eLFhC/egrRe353Oo54iZCqzYktUZvW7W8jLk5HkESulI+PD9O++oLCSceIWj3d6TjiJQIKlSHnrQNp0+E29u3b53ScbEfN61UoV64cM2d8TdSC14k/fsDpOOIlghv14u/Tlr79+usIkkga5MyZkyUL5sK2xUT9+ZPTccRL5CxXD58qbWneuh2RkZFOx8lW1LxepSZNmvD6q6OSpyOMcDqOeAFjfAhu9Rjzlv3E6FdfczqOiKsUL16cRfPmEL38fWL/0Qx2kjly1u7EsRzF6NLtTpKSNN18elHzeg3u79+fXt06EzX/VWyibiYvGc8nIIiQ8GcY8coYZs+e7XQcEVepVasWEz/+kMg5o0iI0PhxyXjGGIJvfYBf/9jD4CFDnY6Tbah5vUZvvzGGGqULE7VC0xFK5vDLVTBlOsKNGzc6HUfEVTp37szgRwcSOWckSfExTscRL2B8/QlpN4QPJ37Fp59+5nScbEHN6zU6Ox1hntO7OLN+jtNxxEvkKFaBwEb30qJNOw4fPux0HBFXee7ZZ2hxU02iFr+DtTqVKxnPN2duQsKfYdBjT/Djjz86Hcf1TFqOFtauXduuXbs23UOEhIQ4Mph5+PDhNG7cmObNm1/zukqUKMGhf/7BhOTHP29xCrR7HN+QvOwffw8+AUHg49lPyNfiQQJLVCTu3z0cXzqBxMhjYJMIrtyM3Dd1xxiTss6Dnw7CP18JCt425JrzSfYU+fNXFI/awa8/fU9gYCC+vr5UrVqV+Ph4/Pz86NWrF4899hg+PmnfTx05ciRPP/10uuQ8efIkX331FQ8++GCaXvfCCy8QEhLCk08+ec0ZjDHrrLW1r3lFLpJRNRvcX7djYmLIlTs3SX5B+ATnxTc47zXV7cSoExxb8A4Jp/+FpER8cxeicFfdqF5Si/57HXHLxrJu9SrKlCkDoLp9EZeq2V595PWll15Kl8YVwM/PjzmzZxNIAr65CnJq1bSU5wr3GEmxvu9SrO+7BJaoSFJ8LP9+8zK563eh+H0TKNr3XWIP/EHkb/NSXhN/dB8kJRG7fwtJcTq1JRcW3KAH++ODuTN5OsKgoCA2bNjAli1bWLJkCQsWLODFF6/uD+jIkSMv+Li1Ns0XHpw8eZJx48ZdVQ6Rc6VX3Q4MDKRI4cLkyxVM7vqdCShS9prq9skfviSwVBjF7hlLsX7jyXtLn2vOKNlPUJla+NXszK2t2nL69GnPY6rbaZalmtcVK1bQpEkTunTpwo033kjPnj1TxpGuWbOGm266ierVq1O3bl0iIiKIiYmhb9++VK1alRo1arB8+XIAPvvsMzp27EiLFi0oVaoUY8eO5Y033qBGjRrUr1+f48ePA9CnTx+mT/fc969UqVI8//zz1KxZk6pVq/LHH38A8O+//9KiRQsqV65Mv379uP766zl69OgF89etW5fx775N0p61xB/de9H3GbV1JTmKVySodE0AfPwDydfiAU6t+v97EEZtW0lwlaYElqpJ9F+rrvGTlezKGENwi0EsX7OJF156OdVzhQoV4oMPPmDs2LFYa0lMTGTw4MHUqVOHatWqMWHCBAAOHTpE48aNCQsLo0qVKvzwww8MHTqU6OhowsLC6NmzJ7t376ZChQr06tWLKlWqsG/fPkJCQlK2NX36dPr06QPA4cOH6dSpE9WrV6d69er8/PPPDB06lJ07dxIWFsbgwYMBeO2111KyPP/88ynrGjFiBOXLl+fmm2/mzz//zOBPUK6Vm+u2j48P06dOJub7T/DNmYv4Ewcv+j4vV7cTo47jG1ogZfmAQqWv9aOVbCq4Znsi8pYn/PauJCSkvthbdfvK+GX4FtLot99+Y8uWLRQrVoyGDRvy008/UbduXbp168bUqVOpU6cOp0+fJigoiLfffhtjDJs2beKPP/6gZcuWbN++HYDNmzfz22+/ERMTQ9myZXnllVf47bffeOyxx5g0aRKPPvrof7ZdoEAB1q9fz7hx4xgzZgwfffQRL774Is2aNWPYsGEsXLiQjz/++JL57777Lt5860227tyLTYgH4PDkp8HHB+PrT9FebxB/dA8BRcqmep1/3qLYuGiSYs/gkyMnUX/8QOFuLxN/bD8R6+YQXKlJuny+kv34+OcguP0w3nh38H8KYZkyZUhMTOTIkSPMmjWL3Llzs2bNGmJjY2nYsCEtW7bkm2++oVWrVjzzzDMkJiZy5swZGjVqxNixY9mwYQMAu3fvZseOHUycOJH69etfMs/DDz/MLbfcwrfffktiYiKRkZGMHj2azZs3p6xv8eLF7Nixg9WrV2OtJTw8nO+//57g4GCmTJnChg0bSEhIoGbNmtSqVSsjPjZJR26u25UqVWLyF5Po1KUrgRWbpjye1rodWqMd/85+lYj1cwm6Pozgqs3xC82fPh+wZDvBt9zLxlkvM+iRx//znOr25WW55rVu3bqUKFECgLCwMHbv3k3u3LkpWrQoderUASBXrlwA/PjjjwwaNAiAG2+8keuvvz6lCDZt2pTQ0FBCQ0PJnTs3HTp0AKBq1aoXvUL79ttvBzy3U/nmm29StvHtt98C0Lp1a/LmzXvR7E2bNk0Zu1KgcBHWffce4Dn95Jsz9xV/BrGHduAblAu/XIXwDcnPsflvkxgdgW9Q6BWvQ7yLX0g+Qjs8zaHPHmHt2rXUrv3fYUKLFy9m48aNKUetTp06xY4dO6hTpw733HMP8fHxdOzYkbCwsAtu4/rrr79sAQRYtmwZkyZNAjxjuXLnzs2JEyf+k2Xx4sXUqFEDgMjISHbs2EFERASdOnUiZ86cAISHh1/xZyDOyQ51u1rlyvx1aBtJsWeAtNftoDK1KH7/R8TsWkf03+s49NkjFLv3vTStQ7yH8fUjpO1gvpr6FPEJF7/Vpur2hWW55jVHjhwp/+/r6/ufI0lXsx4fH5+U7318fC66zrPLXO12ly9fToECntNGUVFR1K7fkJNxUf9Zzr/AdcTu25zqsfiT/2ACgjxHXbetJP7YfvaPvweApLgznPnzJ0LDWqc5k3iPgMI3gK8/rdt1YMO6NZQoUYK///4bX19fChUqhLWWd999l1atWv3ntd9//z3z5s2jT58+PP744/Tq1es/ywQHB6f6/tyLC2Ni0jYu21rLsGHDuP/++1M9/tZbb6VpPZI1ZIe6ba2lV997mbfojQsud7m6DeAbFEpwpSYEV2rCkekvErNvM8EVGqY5k3gHn8AQQsKf5fQH/Vm6dGnKWG7V7cvLUmNeL6ZChQocOnSINWvWABAREUFCQgKNGjXiyy+/BGD79u3s3buXChUqpOu2GzZsyLRpnkH8ixcv/s+eyMUEBwezdOE8iIsheuea1M9VakLs/q1E794AQFJ8LCeWTiBXvc5Ym8SZP36k6D1jKTHgE0oM+IRCtz9H1Lbv0/V9SfZkfHyhchuat27H7t27eeCBBxg4cCDGGFq1asX48eOJj/cMZ9m+fTtRUVHs2bOHwoULc99999GvXz/Wr18PgL+/f8qyF1K4cGG2bdtGUlJSylEugFtvvZXx48cDkJiYyKlTpwgNDSUi4v9nomvVqhWffPJJytXqBw4c4MiRIzRu3JiZM2cSHR1NREQEc+bo9nNu5ba6bYzh4w/ep1xef2zMf++icKm6DRC95/eU+8YmxZ4h/sQh/HIVTKd3JNmVf95iGL8AOt/RnT///JN///1XdfsKZLkjrxcSEBDA1KlTGTRoENHR0QQFBbF06VIefPBBBgwYQNWqVfHz8+Ozzz5LteeeHp5//nl69OjB559/ToMGDShSpAihoVd2+r548eIULlyIUz98SkCh0p4jY3jGKBbs/BzHl7zP8SXjISmJ4MpNCa3Znth9m/ENyZ9qrFSOkpWJP7aXhMjj+IXkS9f3J9mLTYjj9JYVnDx9hGrVqvPcc8/yxBNPANCvXz92795NzZo1sdZSsGBBZs6cyYoVK3jttdfw9/cnJCQk5dRR//79qVatGjVr1mTEiBH/2dbo0aNp3749BQsWpHbt2ikF7e2336Z///58/PHH+Pr6Mn78eBo0aEDDhg2pUqUKbdq04bXXXmPbtm00aNAA8Nx26YsvvqBmzZp069aN6tWrU6hQoZRTzuI+bqzbAQEBLJgzk0KFCnNm6wpCa9+W8tyl6jZA3D9/cWLJ++DjCzaJ0OotyVG0fLq+L8mebGICUQmGKlWrUqF8+ZQjqaC6fTFZ4j6vWVlsbCy+vr74+fnxyy+/MGDAgJTBy1dq2rSvuffBQeTu/pqaT8kUNiGe098Mp1/Xdrz2yiin42QJus+r97jWuv3nn39St0FDglo/QeB11TIuqMg5Ild+Sjm/o3y/bAn+/v5Ox3Gc7vN6Dfbu3UudOnWoXr06Dz/8MB9++GGa13HHHV154uGHiJozCpsQlwEpRVIzfv6EtB/Ke+9/wA8//OB0HJFMda11u0KFCkyfOpnIBa8Tf+JQBqUUSS24US+2HTrFqFdedTpKlueKYQNOKleuHL/99ts1r+f2jrcxatRogqMjdPsUyRQxezYSEhxEuXLlnI4ikqnSo27XrVuX0JBg4o/uxT9v0XRKJnJxCScOknB0L01vaex0lCxPR14zwZEjR2jZtj25WjyoxlUyReyh7USv/IDF8+dSpEiRNL124cKFVKhQgbJlyzJ69Oj/PL9nzx5uvfVWqlWrRpMmTdi/fz/guWo7LCws5SswMJCZM2emx9sRyVQJCQncdntX4otUI6hsXafjiBdIjD5NxOwRvPX6qzRq1CjNr/e2uq3m9Tz33HMPhQoVokqVKumyvpiYGFq27UBC6Zs10YBkioTTR4mYO5ovPvvkovf+u5jExEQeeughFixYwNatW5k8eTJbt25NtcyTTz5Jr1692LhxI8OHD2fYsGGA536ZGzZsYMOGDSxbtoycOXPSsmXL9HpbIheU3jUbYODDj7Fx33GCm/RLdWshkYxgE+OJnPcKfXveQb97703z672xbqt5PU+fPn1YuHBhuqzLWstdvfuyLzaQ4Jt6pMs6RS4lKS6GyDkjGfbkY3Ts2DHNr1+9ejVly5alTJkyBAQE0L17d2bNmpVqma1bt9KsWTPAU/jOfx480w62adMm5abVIhklPWs2wHvjxjH529mEtB2M8dXIOslY1lqilk2gTrkSvPHa1Y119ca6reb1PI0bNyZfvvS5I8BL/xvB0lUbCG75MMboo5aMZW0SUYvfos0t9Xh66JCrWseBAwcoWbJkyvclSpTgwIEDqZapXr16ykxG3377LRERERw7dizVMlOmTKFHD+2wScZLz5q9ZMkShjwznJDwZ/EJDLn8C0SuUdS6WeSL3s83X0/B19f3qtbhjXVbHVUG+frr6Yx5eywh7Yfh4x/odBzxAlE/fUnpkCQ++/jDDD3VOWbMGFauXEmNGjVYuXIlxYsXT1V0Dx06xKZNmy44I4xIVvXnn3/SpVsPQto8gX/eYk7HES9w5q/VJG2cw9KF8wgJydidpexWt3VOJAOsW7eOvvf1J7Tj87pASzJF1ObvCNi7ioXr117TDd+LFy/Ovn37Ur7fv38/xYsXT7VMsWLFUvbgIyMjmTFjBnny5El5ftq0aXTq1En3KRTXOH78OM1btSWg/p26r6tkirgjuzizdCxLF87juuuuu6Z1eWPd1pHXdHbgwAFatetAzqYDyFGkrNNxxAvE7N9C7E8TWbpwPgUKFLimddWpU4cdO3awa9cu4uLimDJlCuHh4amWOXr0KElJSQCMGjWKe+65J9XzkydPds2pJ5H4+HjahXciumgYwdXccdRJ3C0x6gQRc0Yy4b13qV+//jWvzxvrtprXdHTmzBlatGkPFVuSs8JNTscRL5Bw6jCR815j6ldfUKlSpWten5+fH2PHjqVVq1ZUrFiRO+64g8qVKzN8+HBmz54NwIoVK6hQoQLly5fn8OHDPPPMMymv3717N/v27eOWW2655iwiGc1ay739H+DP43EEN+rtdBzxAjYhjsi5oxnY/1569rwzXdbpjXVb08Oep0ePHqxYsYKjR49SuHBhXnzxRe69gltXJCUl0aFTZ1btjSKk5cO6vYpkuKTYM5yaNoQXBj/KY48+7HScLE/Tw2ZPV1uzAca8/gb/e3McubqOwidH1r/CWtzNWkvkordoWCYPM6dPw8dHxw8v5VI1W2NezzN58uSret2wZ57j5007yXX7S2pcJcPZpEQiF77O7W2a8+gjg5yOI+KYq63Z8+bN44URo8nT7RU1rpIpon79mqKcYOqXc9S4XiM1r+ng88+/YNxHn5Gn+6sYP3cMdhZ3i/r+MyoWDOKD8e9pZ0kkjTZv3kz3nncT0v5p/HIXcjqOeIGoP37E/LmUJevXEhQU5HQc11Pzeo1++eUXBgx6hFydX8I3OI/TccQLRP2+kJyHNzJ3/RrXXBkqklUcOXKE5q3bEtj4HgJLVHQ6jniB2EM7iF4xgR+Wf0fRokWdjpMtqHm9Bnv27KFdeEeCWwwioGApp+OIF4je8ztxv07m519/IW/evE7HEXGVmJgYWrULJ6FUQ0IqNXU6jniBhIijRMwdxaRPPqJGjRpOx8k2NOjiKkVERHBryzb4ht1G0A11nI4jXiD++AGiFrzOzBlfU65cOafjiLiKtZZefe9lb0wOghumz1XeIpdydrruIY8/wu233+50nGxFzetVSExMpGPnOzgZWpqcNcMv/wKRa5QYE0nE7BGMeWUkTZvqiJFIWv1vxEgW/7RO03VLpvBM1/02rW6uw7NPD3M6Trajf8FX4dEnBvPbrsMEN+2vi2Ukw9nEBKLmvcLdd3TigfvvdzqOiOvMmDGDV998l5AOmq5bMkfUz19xfc54Jn36sfqEDKAxr2n0wYcfMnHKdHJ3exXjq49PMpa1lqgVH1CjdGHeefN1p+OIuM769evpfe99ydN1X9sMdCJXImrLcvx3/8yia5yuWy5O3VcarFixgscHDyX3HaPwDQp1Oo54gaj1c8l9ehfffrcKX19fp+OIuMrBgwdp1bY9QU3u13Tdkili9m8j5sdPWfXj9xQsWNDpONmWmtcr9Ndff9Gxc1eCWz+Of77iTscRLxC9cy2Jv33Dd2tXkytXLqfjiLjK2em67Y0tCLnxZqfjiBdIOHWYyPmvMu3Lz6lcubLTcbI1Na9X4MSJE9zasg3+dbsTVCrM6TjiBeL+3UPUkndYNG82pUqVcjqOiKskJSXRrefd/GPyEVKvq9NxxAskxZ4hYvYIXnx2GG3btnU6TranC7YuIz4+nvYdOxNVqArB1Vs7HUe8QOKZU0TMGcG4d96iYcOGTscRcZ2nnx3Ojxv+JKT5Q7pYRjJcynTdrZvx+GOPOh3HK6h5vYwBDw1i2+Eoghv3dTqKeAGbEE/k3NE80LcXvXrd7XQcEdf54osvee/DTwhpPxTjF+B0HPECUT9MpEKBQD54f5x2ljKJhg1cwltvv8vXcxeTu9tojI8ulpGMZa0l8rv3aFC5DK+MGuF0HBHX+eWXX3jgoUHk6vIyvsGagU4yXtTGRQQd2sA8TdedqdS8XsTChQt59oWXyN1tND45gp2OI14gas03FEr4l6+nzMTHRydFRNLi7HTdOTVdt2SSmD0biVv1FT+t+pl8+fI5HcerqHm9gK1bt3JHj56EtBuCf54iTscRL3Bm+y+wZQFL168lOFg7SyJpERERQfPW7fCpHk7OsnWdjiNeIP74ASIXjGHWjGmUL1/e6TheR83reY4ePUrz1m3J0bA3gSV0qwvJeHGHd3Jm2ThWLF1M8eK6DZtIWiQmJnJ71+4cz1mSkFq3OR1HvMDZ6bpfHfU/br31VqfjeCWdmzxHXFwcrdvfRlzJugRX0S+kZLyEyONEzBnJJx9MoHbt2k7HEXGdx558irV/HSCk2f26WEYynE1MIGr+a/TsfBsPDhjgdByvpeY1mbWWPvfex64IQ/DNuspbMl5SfCyRc0bx+KAHueMO3YtSJK0++PBDPvtqGiHthmB8dbGMZCzPdN0fUe26/Ix9502n43g1DRtINuqVV5m/4hdydR2JMerpJWNZa4la8i7N6lXjheHPOR1HxHVWrlzpma676wh8gzQDnWS8M7/NJdepv5j93a+artthal6BWbNmMeq1N8jd7VV8AoKcjiNeIOqXyZTwj+KrSZ/pVKdIGu3cuZPbbu/ima47f0mn44gXiP57HQnrv+G7Nb9quu4swOub1w0bNnBX776EhD+LX64CTscRLxC1dSW+O79n8fq1BAYGOh1HxFVOnjzpma67zh2arlsyRdzRvUQtfpuFc2dRunRpp+MIXt68/vPPP7Rs257AW+4jR7EKTscRLxB78E9ifviYn1Yup3Dhwk7HEXGVhIQEOnTqQmT+ioSEaf54yXiJZ04ROXsEY996g5tvvtnpOJLMawd3RkdH07JtB5LKNSG4YmOn44gXSDj9LxFzR/PlxE+pVq2a03FEXGfAwEFsOXia4Cb3Oh1FvIBNiCdy3ivc1/tO+vTp7XQcOYdXNq/WWnr26sOBxBCC63d3Oo54gaS4aCLnjOCZIU8SHh7udBwR13nn3bFMm72QkLZParpuyXDWWiKXjad+xVK89spop+PIebyyeR3+wossX7OZ4OaDdLGMZDhrk4ha9CbtmjZk6FODnY4j4jqLFy/m6eEvENLhGU3XLZkias03FIz7h+lTv9J03VmQ1415nTJlKm+9N4E83V/Dxz+H03HEC0T9MIkbcvvw6UcfaGdJJI3++OMPunTrQUjbp/DPW9TpOOIFzmz/BbtlAUvXrdF03VmUVzWvq1evpt8DA8jV6UV8Q/I6HUe8QOSmpeQ4sJYF69cSEBDgdBwRVzl27Bi3tmpDwE13E1iyitNxxAvEHf6bM8vGsXzJIkqUKOF0HLkIr2le9+/fT5v24eS89SECCpdxOo54gZh9m4n7eRI//vIT+fPndzqOiKucna47tnhtQqq2cDqOeIHEyBNEzBnJR++Pp06dOk7HkUvwioEcUVFRNG/dDlOlLTnL1Xc6jniB+JP/EDl/DF9P+Yobb7zR6TgirmKtpW+//uw8ZTVdt2SKpPhYIueO5NGH7qd7925Ox5HLyPbNa1JSEp279eBoQFGC69zudBzxAkmxUUTO/h8jXnyeVq1aOR1HxHVeeW0M85b/REjrx3RnAclw1lqilr5L0zpVeemF552OI1cg2zevg4cM49dtewi+9QFdLCMZziYlEjl/DHeEt+aRhwc6HUfEdWbPns3/Rr9GSIenNV23ZIqoVVMo7hvJl5qu2zWy9ZjXTz/9jA8nfknu7q9ifP2djiNXac8r7f/zmF9A1rpThDE+5O/6MnHbv6dysVyMH/uu05FEXGfjxo307NWH0PBn8ctVyOk4cpWOTH+R6J1rUj3m5x8AWawxzNukDzYwNz47VrB4/VqCgrSz5BbZtnn98ccfuf+BAeRqMYCEU0dIOHXE6UhylQrfORqbEE9S9CkSDv3JmS3LCAurwUcTxlG2bFmn4wHQsWt3fvjxcwr5nGHOutX4+WXbf1oiGeLw4cM0b90Wn/JNwPgQe2iH05HkKuWq35XQOh1Jio4g4dheErYtJ3dwDiaMG0uzZs2cjgfAuHHjGP72J5jTh/hp5XKKFCnidCRJg2z7F3bfgUOULlcedn8Hu51OI9ciISGRiIjTREWcJmelpuTvPpqde3+nwc2NeWPMq/S/7z6nI+Ln64fP8T0sXbeGPHnyOB1HxHUiIyPJl78A9sxOWLfT6ThyDZKSLFGREZw69i8hN9QmR4tBxMWdoXuve+nR9XbGjX3H8R38gIAAzuzfxsxvv6V69eqOZpG0M9baK164du3adu3atRkYR+Ti9u7dy3vj32fCBx/hV7gs/rU6E7XoTYYPeZwnn3jc0Wy/rl6Dj0G3V8nCjDHrrLW1nc6RmVSzxUmnT59m4qRJvP7WWE5GxxPY9AHi1nxNw0olmTl9Gr6+zl2Mt3//ftatW89tt2m67qzqUjVbzau4TlxcHI89+RRfTJ9FzuaDiJj3Ct9M/pyWLVs6HU2yMDWvIs6w1jJx4iQGPvo4IW2eJHbN1wzo3o6R/3vZ6WiShV2qZmf7uw1I9hMQEMB777zFoH69iF4xgZyN+3FXn3s4ffq009FEROQ8xhj69OnN15O/IHLh6+So25W3x47jt99+czqauJSaV3Gtl198gQbVbyTx9D8kFijHhA8+dDqSiIhcRJs2bXj+mWEkbJxHYN07GPqc7qkqV0fNq7iWMcZTCDcvxr9aG9569z3SMgxGREQy14MDHiBm3xZyXFeNH1auZN++fU5HEhdS8yquVq9ePYoUyIuNiyEiKpqdO3WVsohIVhUcHMxdPe8i9o+VhJStzdKlS52OJC6k5lVczRjDA/36krTrVwKLlkUXp4iIZG333duHpF2rSchbmp9/Xe10HHEhNa/ieiVLlsQnNpKEkELs3bvX6TgiInIJxYsXJy7yJH55i7Jj526n44gLqXkV1wsODsYmxGCNH/Hx8U7HERGRSwgODiYhLgbj46uaLVdFzau4njEGA5C1ps0WEZELMOZssVbRlquj5lVEREREXEPNq4iIiIi4hppXEREREXENNa8iIiIi4hpqXkVERETENdS8ioiIiIhrqHkVEREREddQ8yoiIiIirqHmVURERERcQ82riIiIiLiGmlcRERERcQ01ryIiIiLiGmpeRURERMQ11LyKiIiIiGuoeRURERER11DzKiIiIiKuoeZVRERERFxDzauIiIiIuIaaVxERERFxDT+nA4hcrc2bN7N582Y2bNiQ8tjGjRuZMmUKxYsXp1GjRs6FExGRVA4fPszy5cuJiYlJeezfI/8wZcoUAgICuP322x1MJ26i5lVc6+2x4/hownjy3ViPgBubEBhSiOUb5/D1XXdTo1Zt1v36i9MRRUQk2bLly7mzRw9CS1QgT73bCShSlpN/FaV3/wfxS4qnWbNm5MmTx+mY4gJqXsW1xr37NgkJCXyz4DuCilTAJzCU2LgYGtx0M3NnfeN0PBEROUeP7t05ffo0Tzw1DIpWwidnLpKC85MvdyhLF85X4ypXTM2ruJa/vz+ffDiBqm+9w/CXhuAbFMrtbZvzwfj38Pf3dzqeiIic5/7+/SlfrhwdO3clJndRbigYwqL1aylYsKDT0cRF1LyKqxljePyxR6hWtTIHDh6k1913Y4xxOpaIiFxE06ZNWfvrL3zz7UwefeRhcuTI4XQkcRk1r5ItNG/e3OkIIiJyhcqVK8eQpwY7HUNcSrfKEhERERHXUPMqIiIiIq6h5lVEREREXEPNq4iIiIi4hppXcaWZM2dijOGPP/4A4ODBg3Tp0sXhVCIiciGq2ZKe1LyKK02ePJmbb76ZyZMnA1CsWDGmT5/+n+USEhIyO5qIiJxHNVvSk5pXcZ3IyEh+/PFHPv74Y6ZMmQLA7t27qVKlCgCfffYZ4eHhNGvWjFtvvdXJqCIiXk81W9Kb7vMqrjNr1ixat25N+fLlyZ8/P+vWrSN//vypllm/fj0bN24kX758DqUUERFQzZb0pyOv4jqTJ0+me/fuAHTv3j3lNNS5WrRooSIoIpIFqGZLetORV3GV48ePs2zZMjZt2oQxhsTERIwxPPTQQ6mWCw4OdiihiIicpZotGUFHXsVVpk+fzt13382ePXvYvXs3+/bto3Tp0uzbt8/paCIich7VbMkIal7FVSZPnkynTp1SPda5c2dGjRrlUCIREbkY1WzJCMZae8UL165d265duzYD44iIZAxjzDprbW2nc2Qm1WwRcatL1WwdeRURERER11DzKiIiIiKuoeZVRERERFxDzas4wtfXl7CwsJSv0aNHX/Frr3ZO7H79+rF169Y0v05ExNupZktWogu2xBEhISFERkY6HUO8iC7YErl6qtmS2XTBlrhGqVKlGDZsGGFhYdSuXZv169fTqlUrbrjhBt5//30g9ZzYW7ZsoW7duoSFhVGtWjV27NhBVFQU7dq1o3r16lSpUoWpU6cC0KRJE87+IZ88eTJVq1alSpUqDBkyJGX7ISEhPPPMM1SvXp369etz+PDhTP4ERETcQzVbnKDmVRwRHR2d6hTU2WIFcN1117FhwwYaNWpEnz59mD59OqtWreL555//z3ref/99HnnkETZs2MDatWspUaIECxcupFixYvz+++9s3ryZ1q1bp3rNwYMHGTJkCMuWLWPDhg2sWbOGmTNnAhAVFUX9+vX5/fffady4MR9++GGGfg4iIm6gmi1ZiaaHFUcEBQWxYcOGCz4XHh4OQNWqVYmMjCQ0NJTQ0FBy5MjByZMnUy3boEEDRowYwf79+7n99tspV64cVatW5YknnmDIkCG0b9+eRo0apXrNmjVraNKkCQULFgSgZ8+efP/993Ts2JGAgADat28PQK1atViyZEn6vnERERdSzZasREdeJcvJkSMHAD4+Pin/f/b7hISEVMveeeedzJ49m6CgINq2bcuyZcsoX74869evp2rVqjz77LO89NJLV7xtf39/jDGA5wKF87cnIiKpqWZLZlPzKq72999/U6ZMGR5++GFuu+02Nm7cyMGDB8mZMyd33XUXgwcPZv369aleU7duXVauXMnRo0dJTExk8uTJ3HLLLQ69AxER76GaLelBwwbEEWfHT53VunXrNN165axp06bx+eef4+/vT5EiRXj66adZs2YNgwcPxsfHB39/f8aPH5/qNUWLFmX06NE0bdoUay3t2rXjtttuu9a3JCKSbalmS1aiW2WJiFfQrbJERNxDt8oSERERkWxBzauIiIiIuIaaV3GFhQsXUqFCBcqWLXvBcVaxsbF069aNsmXLUq9ePXbv3g3AkiVLqFWrFlWrVqVWrVosW7Ysk5OLiHgf1WzJSGpeJctLTEzkoYceYsGCBWzdupXJkyf/Z77rjz/+mLx58/LXX3/x2GOPpczAUqBAAebMmcOmTZuYOHEid999txNvQUTEa6hmS0ZT8ypZ3urVqylbtixlypQhICCA7t27M2vWrFTLzJo1i969ewPQpUsXvvvuO6y11KhRg2LFigFQuXJloqOjiY2NzfT3ICLiLVSzJaOpeZUs78CBA5QsWTLl+xIlSnDgwIGLLuPn50fu3Lk5duxYqmVmzJhBzZo1U91EW0RE0pdqtmQ03edVvMKWLVsYMmQIixcvdjqKiIhchmq2XIqOvEqWV7x4cfbt25fy/f79+ylevPhFl0lISODUqVPkz58/ZflOnToxadIkbrjhhswLLiLihVSzJaOpeZUsr06dOuzYsYNdu3YRFxfHlClTCA8PT7VMeHg4EydOBGD69Ok0a9YMYwwnT56kXbt2jB49moYNGzoRX0TEq6hmS0ZT8ypZnp+fH2PHjqVVq1ZUrFiRO+64g8qVKzN8+HBmz54NwL333suxY8coW7Ysb7zxRsqtWcaOHctff/3FSy+9RFhYGGFhYRw5csTJtyMikq2pZktG0/SwIuIVND2siIh7aHpYEREREckW1LyKiIiIiGuoeRURERER11DzKiIiIiKuoeZVRERERFxDzauIiIiIuIaaVxERERFxDTWvIiIiIuIaal5FRERExDXUvIqIiIiIa6h5FRERERHXUPMqIiIiIq6h5lVEREREXEPNq4iIiIi4hppXEREREXENNa8iIiIi4hpqXkVERETENdS8ioiIiIhrqHkVEREREddQ8yoiIiIirmGstVe+sDH/AnsyLo6ISIa53lpb0OkQmUk1W0Rc7KI1O03Nq4iIiIiIkzRsQERERERcQ82riIiIiLiGmlcRERERcQ01ryIiIiLiGmpeRURERMQ11LyKiIiIiGuoeRURERER11DzKiIiIiKuoeZVRERERFzj/wAEz78EpDN8PAAAAABJRU5ErkJggg==\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.7.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
