3_plotting.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plotting\n",
    "\n",
    "This notebook explains how to plot LCIA results. For this, we first choose activities and methods to perform an LCIA. Then, we plot the result using pandas."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Imports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "import brightway2 as bw"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Choosing activities\n",
    "\n",
    "In this example, we want to compare three ecoinvent activities. Don't forget to activate the brightway project into which you have imported ecoinvent previously."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# activate project\n",
    "bw.projects.set_current('ab')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Next, we look through the ecoinvent database and choose three activities to compare. In this case, electricity production from wind, PV and natural gas in Germany."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Excluding 217 filtered results\n",
      "Excluding 116 filtered results\n",
      "Excluding 562 filtered results\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "('electricity production, wind, >3MW turbine, onshore' (kilowatt hour, DE, None),\n",
       " 'electricity production, photovoltaic, 3kWp slanted-roof installation, multi-Si, panel, mounted' (kilowatt hour, DE, None),\n",
       " 'electricity production, natural gas, combined cycle power plant' (kilowatt hour, DE, None))"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "db_name = 'apos37'\n",
    "act1 = bw.Database(db_name).search(\n",
    "    'electricity production, wind, >3MW', \n",
    "    filter={'location':'DE'}\n",
    ")[0]\n",
    "act2 = bw.Database(db_name).search(\n",
    "    'electricity production, photovoltaic, 3kWp slanted-roof installation, multi-Si', \n",
    "    filter={'location':'DE'}\n",
    ")[0]\n",
    "act3 = bw.Database(db_name).search(\n",
    "    'electricity production, natural gas, combined cycle power plant', \n",
    "    filter={'location':'DE'}\n",
    ")[0]\n",
    "\n",
    "act1, act2, act3"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Choosing methods\n",
    "\n",
    "Next, we choose methods for the impact assessment. In this case ILCD 2.0 including long-term impacts."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[('ILCD 2.0 2018 midpoint', 'climate change', 'climate change biogenic'),\n",
       " ('ILCD 2.0 2018 midpoint', 'climate change', 'climate change fossil'),\n",
       " ('ILCD 2.0 2018 midpoint',\n",
       "  'climate change',\n",
       "  'climate change land use and land use change'),\n",
       " ('ILCD 2.0 2018 midpoint', 'climate change', 'climate change total'),\n",
       " ('ILCD 2.0 2018 midpoint',\n",
       "  'ecosystem quality',\n",
       "  'freshwater and terrestrial acidification'),\n",
       " ('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater ecotoxicity'),\n",
       " ('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater eutrophication'),\n",
       " ('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'marine eutrophication'),\n",
       " ('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'terrestrial eutrophication'),\n",
       " ('ILCD 2.0 2018 midpoint', 'human health', 'carcinogenic effects'),\n",
       " ('ILCD 2.0 2018 midpoint', 'human health', 'ionising radiation'),\n",
       " ('ILCD 2.0 2018 midpoint', 'human health', 'non-carcinogenic effects'),\n",
       " ('ILCD 2.0 2018 midpoint', 'human health', 'ozone layer depletion'),\n",
       " ('ILCD 2.0 2018 midpoint', 'human health', 'photochemical ozone creation'),\n",
       " ('ILCD 2.0 2018 midpoint', 'human health', 'respiratory effects, inorganics'),\n",
       " ('ILCD 2.0 2018 midpoint', 'resources', 'dissipated water'),\n",
       " ('ILCD 2.0 2018 midpoint', 'resources', 'fossils'),\n",
       " ('ILCD 2.0 2018 midpoint', 'resources', 'land use'),\n",
       " ('ILCD 2.0 2018 midpoint', 'resources', 'minerals and metals')]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "methods = sorted(\n",
    "    [m for m in bw.methods if \"ILCD 2.0 2018 midpoint\" in str(m) and \n",
    "     'no LT' not in str(m) and 'incl. biogenic' not in str(m)]\n",
    ")\n",
    "methods"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculating the LCIA scores\n",
    "\n",
    "Next, we define a function that does the LCIA calculations for us. We saw how to do this in a previous notebook, so I won't go into details here."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "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>'electricity production, wind, &gt;3MW turbine, onshore' (kilowatt hour, DE, None)</th>\n",
       "      <th>'electricity production, photovoltaic, 3kWp slanted-roof installation, multi-Si, panel, mounted' (kilowatt hour, DE, None)</th>\n",
       "      <th>'electricity production, natural gas, combined cycle power plant' (kilowatt hour, DE, None)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'climate change', 'climate change biogenic')</th>\n",
       "      <td>8.200771e-05</td>\n",
       "      <td>4.388205e-04</td>\n",
       "      <td>6.532204e-05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'climate change', 'climate change fossil')</th>\n",
       "      <td>3.148162e-02</td>\n",
       "      <td>9.537874e-02</td>\n",
       "      <td>4.300353e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'climate change', 'climate change land use and land use change')</th>\n",
       "      <td>1.100778e-04</td>\n",
       "      <td>2.096568e-04</td>\n",
       "      <td>1.278100e-05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'climate change', 'climate change total')</th>\n",
       "      <td>3.167371e-02</td>\n",
       "      <td>9.602722e-02</td>\n",
       "      <td>4.301134e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater and terrestrial acidification')</th>\n",
       "      <td>2.068838e-04</td>\n",
       "      <td>5.807309e-04</td>\n",
       "      <td>4.974314e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater ecotoxicity')</th>\n",
       "      <td>1.550215e-01</td>\n",
       "      <td>1.618358e-01</td>\n",
       "      <td>2.567960e-02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater eutrophication')</th>\n",
       "      <td>2.981222e-05</td>\n",
       "      <td>6.173620e-05</td>\n",
       "      <td>3.639684e-06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'marine eutrophication')</th>\n",
       "      <td>4.760436e-05</td>\n",
       "      <td>1.247935e-04</td>\n",
       "      <td>1.278583e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'terrestrial eutrophication')</th>\n",
       "      <td>5.081589e-04</td>\n",
       "      <td>1.184353e-03</td>\n",
       "      <td>1.383222e-03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'human health', 'carcinogenic effects')</th>\n",
       "      <td>9.542035e-09</td>\n",
       "      <td>6.680907e-09</td>\n",
       "      <td>1.735982e-09</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'human health', 'ionising radiation')</th>\n",
       "      <td>2.678976e-03</td>\n",
       "      <td>9.685159e-03</td>\n",
       "      <td>6.938518e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'human health', 'non-carcinogenic effects')</th>\n",
       "      <td>2.094892e-08</td>\n",
       "      <td>3.554039e-08</td>\n",
       "      <td>2.651108e-09</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'human health', 'ozone layer depletion')</th>\n",
       "      <td>2.548350e-09</td>\n",
       "      <td>1.204897e-08</td>\n",
       "      <td>7.988278e-08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'human health', 'photochemical ozone creation')</th>\n",
       "      <td>1.474750e-04</td>\n",
       "      <td>3.937023e-04</td>\n",
       "      <td>5.045739e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'human health', 'respiratory effects, inorganics')</th>\n",
       "      <td>2.661299e-09</td>\n",
       "      <td>5.979152e-09</td>\n",
       "      <td>9.775673e-10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'resources', 'dissipated water')</th>\n",
       "      <td>1.417678e-02</td>\n",
       "      <td>1.322061e-01</td>\n",
       "      <td>2.711740e-02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'resources', 'fossils')</th>\n",
       "      <td>4.191031e-01</td>\n",
       "      <td>1.454789e+00</td>\n",
       "      <td>7.332433e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'resources', 'land use')</th>\n",
       "      <td>1.094371e+00</td>\n",
       "      <td>2.264004e+00</td>\n",
       "      <td>1.477621e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>('ILCD 2.0 2018 midpoint', 'resources', 'minerals and metals')</th>\n",
       "      <td>1.138813e-05</td>\n",
       "      <td>2.880110e-05</td>\n",
       "      <td>3.271545e-07</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                                    'electricity production, wind, >3MW turbine, onshore' (kilowatt hour, DE, None)  \\\n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       8.200771e-05                                 \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       3.148162e-02                                 \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       1.100778e-04                                 \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       3.167371e-02                                 \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       2.068838e-04                                 \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       1.550215e-01                                 \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       2.981222e-05                                 \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       4.760436e-05                                 \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       5.081589e-04                                 \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'car...                                       9.542035e-09                                 \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'ion...                                       2.678976e-03                                 \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'non...                                       2.094892e-08                                 \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'ozo...                                       2.548350e-09                                 \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'pho...                                       1.474750e-04                                 \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'res...                                       2.661299e-09                                 \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'dissip...                                       1.417678e-02                                 \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'fossils')                                       4.191031e-01                                 \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'land u...                                       1.094371e+00                                 \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'minera...                                       1.138813e-05                                 \n",
       "\n",
       "                                                    'electricity production, photovoltaic, 3kWp slanted-roof installation, multi-Si, panel, mounted' (kilowatt hour, DE, None)  \\\n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       4.388205e-04                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       9.537874e-02                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       2.096568e-04                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       9.602722e-02                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       5.807309e-04                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       1.618358e-01                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       6.173620e-05                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       1.247935e-04                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       1.184353e-03                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'car...                                       6.680907e-09                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'ion...                                       9.685159e-03                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'non...                                       3.554039e-08                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'ozo...                                       1.204897e-08                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'pho...                                       3.937023e-04                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'res...                                       5.979152e-09                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'dissip...                                       1.322061e-01                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'fossils')                                       1.454789e+00                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'land u...                                       2.264004e+00                                                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'minera...                                       2.880110e-05                                                                            \n",
       "\n",
       "                                                    'electricity production, natural gas, combined cycle power plant' (kilowatt hour, DE, None)  \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       6.532204e-05                                            \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       4.300353e-01                                            \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       1.278100e-05                                            \n",
       "('ILCD 2.0 2018 midpoint', 'climate change', 'c...                                       4.301134e-01                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       4.974314e-04                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       2.567960e-02                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       3.639684e-06                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       1.278583e-04                                            \n",
       "('ILCD 2.0 2018 midpoint', 'ecosystem quality',...                                       1.383222e-03                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'car...                                       1.735982e-09                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'ion...                                       6.938518e-04                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'non...                                       2.651108e-09                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'ozo...                                       7.988278e-08                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'pho...                                       5.045739e-04                                            \n",
       "('ILCD 2.0 2018 midpoint', 'human health', 'res...                                       9.775673e-10                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'dissip...                                       2.711740e-02                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'fossils')                                       7.332433e+00                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'land u...                                       1.477621e-01                                            \n",
       "('ILCD 2.0 2018 midpoint', 'resources', 'minera...                                       3.271545e-07                                            "
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "def calculate_lcia_results(activities, methods):\n",
    "    \"\"\"\n",
    "    Calculates LCIA results for all activities in label_act and all methods in label_method\n",
    "    :param label_act: dict containing label:activity pairs\n",
    "    :param label_m: dict containing label:method pairs\n",
    "    :return: dataframe containing lcia results\n",
    "    \"\"\"\n",
    "\n",
    "    # prepare result store\n",
    "    results = {}\n",
    "\n",
    "    # start calculation\n",
    "    for act in activities:\n",
    "        # define functional unit (final demand)\n",
    "        fu = {bw.get_activity(act): 1}\n",
    "        # instantiate lca object\n",
    "        lca = bw.LCA(fu)\n",
    "        # do LCI (matrix inversion)\n",
    "        lca.lci()\n",
    "        # extend results object\n",
    "        results[str(act)] = {}\n",
    "        for  method in methods:\n",
    "            # activate current lcia method\n",
    "            lca.switch_method(method)\n",
    "            # calculate lcia result\n",
    "            lca.lcia()\n",
    "            # store result\n",
    "            results[str(act)][str(method)] = lca.score\n",
    "\n",
    "    # convert dict of results to dataframe\n",
    "    return pd.DataFrame(results)\n",
    "\n",
    "df = calculate_lcia_results([act1, act2, act3], methods)\n",
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's define some short-hands for the methods and activities to make the table a little more overseeable."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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>wind</th>\n",
       "      <th>PV</th>\n",
       "      <th>nat gas</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>CC\\nbio</th>\n",
       "      <td>8.200771e-05</td>\n",
       "      <td>4.388205e-04</td>\n",
       "      <td>6.532204e-05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CC\\nfossil</th>\n",
       "      <td>3.148162e-02</td>\n",
       "      <td>9.537874e-02</td>\n",
       "      <td>4.300353e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CC\\nLUC</th>\n",
       "      <td>1.100778e-04</td>\n",
       "      <td>2.096568e-04</td>\n",
       "      <td>1.278100e-05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>CC\\ntotal</th>\n",
       "      <td>3.167371e-02</td>\n",
       "      <td>9.602722e-02</td>\n",
       "      <td>4.301134e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EQ\\nacid</th>\n",
       "      <td>2.068838e-04</td>\n",
       "      <td>5.807309e-04</td>\n",
       "      <td>4.974314e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EQ\\ntox</th>\n",
       "      <td>1.550215e-01</td>\n",
       "      <td>1.618358e-01</td>\n",
       "      <td>2.567960e-02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EQ\\nFW</th>\n",
       "      <td>2.981222e-05</td>\n",
       "      <td>6.173620e-05</td>\n",
       "      <td>3.639684e-06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EQ\\nmar</th>\n",
       "      <td>4.760436e-05</td>\n",
       "      <td>1.247935e-04</td>\n",
       "      <td>1.278583e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>EQ\\nterr</th>\n",
       "      <td>5.081589e-04</td>\n",
       "      <td>1.184353e-03</td>\n",
       "      <td>1.383222e-03</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>HH\\nCE</th>\n",
       "      <td>9.542035e-09</td>\n",
       "      <td>6.680907e-09</td>\n",
       "      <td>1.735982e-09</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>HH\\nrad</th>\n",
       "      <td>2.678976e-03</td>\n",
       "      <td>9.685159e-03</td>\n",
       "      <td>6.938518e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>HH\\nNCE</th>\n",
       "      <td>2.094892e-08</td>\n",
       "      <td>3.554039e-08</td>\n",
       "      <td>2.651108e-09</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>HH\\nODP</th>\n",
       "      <td>2.548350e-09</td>\n",
       "      <td>1.204897e-08</td>\n",
       "      <td>7.988278e-08</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>HH\\nPOC</th>\n",
       "      <td>1.474750e-04</td>\n",
       "      <td>3.937023e-04</td>\n",
       "      <td>5.045739e-04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>HH\\nresp</th>\n",
       "      <td>2.661299e-09</td>\n",
       "      <td>5.979152e-09</td>\n",
       "      <td>9.775673e-10</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>res\\nwater</th>\n",
       "      <td>1.417678e-02</td>\n",
       "      <td>1.322061e-01</td>\n",
       "      <td>2.711740e-02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>res\\nfoss</th>\n",
       "      <td>4.191031e-01</td>\n",
       "      <td>1.454789e+00</td>\n",
       "      <td>7.332433e+00</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>res\\nland</th>\n",
       "      <td>1.094371e+00</td>\n",
       "      <td>2.264004e+00</td>\n",
       "      <td>1.477621e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>res\\nMM</th>\n",
       "      <td>1.138813e-05</td>\n",
       "      <td>2.880110e-05</td>\n",
       "      <td>3.271545e-07</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                    wind            PV       nat gas\n",
       "CC\\nbio     8.200771e-05  4.388205e-04  6.532204e-05\n",
       "CC\\nfossil  3.148162e-02  9.537874e-02  4.300353e-01\n",
       "CC\\nLUC     1.100778e-04  2.096568e-04  1.278100e-05\n",
       "CC\\ntotal   3.167371e-02  9.602722e-02  4.301134e-01\n",
       "EQ\\nacid    2.068838e-04  5.807309e-04  4.974314e-04\n",
       "EQ\\ntox     1.550215e-01  1.618358e-01  2.567960e-02\n",
       "EQ\\nFW      2.981222e-05  6.173620e-05  3.639684e-06\n",
       "EQ\\nmar     4.760436e-05  1.247935e-04  1.278583e-04\n",
       "EQ\\nterr    5.081589e-04  1.184353e-03  1.383222e-03\n",
       "HH\\nCE      9.542035e-09  6.680907e-09  1.735982e-09\n",
       "HH\\nrad     2.678976e-03  9.685159e-03  6.938518e-04\n",
       "HH\\nNCE     2.094892e-08  3.554039e-08  2.651108e-09\n",
       "HH\\nODP     2.548350e-09  1.204897e-08  7.988278e-08\n",
       "HH\\nPOC     1.474750e-04  3.937023e-04  5.045739e-04\n",
       "HH\\nresp    2.661299e-09  5.979152e-09  9.775673e-10\n",
       "res\\nwater  1.417678e-02  1.322061e-01  2.711740e-02\n",
       "res\\nfoss   4.191031e-01  1.454789e+00  7.332433e+00\n",
       "res\\nland   1.094371e+00  2.264004e+00  1.477621e-01\n",
       "res\\nMM     1.138813e-05  2.880110e-05  3.271545e-07"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "label_method = {\n",
    "    \"('ILCD 2.0 2018 midpoint', 'climate change', 'climate change biogenic')\": 'CC\\nbio', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'climate change', 'climate change fossil')\": 'CC\\nfossil', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'climate change', 'climate change land use and land use change')\": 'CC\\nLUC', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'climate change', 'climate change total')\": 'CC\\ntotal', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater and terrestrial acidification')\": 'EQ\\nacid', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater ecotoxicity')\": 'EQ\\ntox', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'freshwater eutrophication')\": 'EQ\\nFW', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'marine eutrophication')\": 'EQ\\nmar', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'ecosystem quality', 'terrestrial eutrophication')\": 'EQ\\nterr', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'human health', 'carcinogenic effects')\": 'HH\\nCE', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'human health', 'ionising radiation')\": 'HH\\nrad', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'human health', 'non-carcinogenic effects')\": 'HH\\nNCE', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'human health', 'ozone layer depletion')\": 'HH\\nODP', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'human health', 'photochemical ozone creation')\": 'HH\\nPOC', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'human health', 'respiratory effects, inorganics')\": 'HH\\nresp', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'resources', 'dissipated water')\": 'res\\nwater', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'resources', 'fossils')\": 'res\\nfoss', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'resources', 'land use')\": 'res\\nland', \n",
    "    \"('ILCD 2.0 2018 midpoint', 'resources', 'minerals and metals')\": 'res\\nMM'\n",
    "}\n",
    "\n",
    "label_act = {\n",
    "    \"'electricity production, wind, >3MW turbine, onshore' (kilowatt hour, DE, None)\": \"wind\",\n",
    "    \"'electricity production, photovoltaic, 3kWp slanted-roof installation, multi-Si, panel, mounted' (kilowatt hour, DE, None)\": \"PV\",\n",
    "    \"'electricity production, natural gas, combined cycle power plant' (kilowatt hour, DE, None)\": \"nat gas\",\n",
    "}\n",
    "\n",
    "df = df.rename(columns=label_act, index=label_method)\n",
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Much better! \n",
    "\n",
    "## Plotting the result\n",
    "\n",
    "Now, let's plot the result. Luckily, pandas offers this out of the box."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16,\n",
       "        17, 18]),\n",
       " <a list of 19 Text major ticklabel objects>)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAzYAAAHrCAYAAAAdcAiZAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3debgedX03/vfHEAmbpGq0VYqhuKFGokQU2+BaqmKt60+pUFELT59qi9al6vNUafu0WpfWUheKVaQVoSriitbWoqJiSwLBBHFDY4tajbQiW5Tl+/tj5sAhnCQnyZlzMjmv13Xd17ln7rlnPmfubd4z3/lOtdYCAAAwZreb6wIAAAB2lGADAACMnmADAACMnmADAACMnmADAACM3m5zXcBkd77zndvSpUvnugwAAGAntXr16h+11pZsOn6nCjZLly7NqlWr5roMAABgJ1VV35lqvKZoAADA6Ak2AADA6Ak2AADA6O1U59gAAMB8cf311+fyyy/Pxo0b57qUndKiRYuy3377ZeHChdOaXrABAIA5cPnll2efffbJ0qVLU1VzXc5OpbWWK664IpdffnkOOOCAaT1HUzQAAJgDGzduzJ3udCehZgpVlTvd6U7bdDRLsAEAgDki1Gzetq4bwQYAABg959gAAMBOYOkrPj6j81v/uiN36PlPeMIT8t73vjeLFy+e3vLWr88Tn/jErFu3boeWu70EGwAA4DbOOeecuS5hm2iKBgAA89DrX//6nHTSSUmSF7/4xXn0ox+dJPn0pz+do48+OkuXLs2PfvSjrF+/PgcddFCOO+643P/+988RRxyR6667LkmyevXqHHzwwTnssMPy1re+dc7+l0SwAQCAeenwww/PeeedlyRZtWpVrr766lx//fX5/Oc/n5UrV95q2m984xt5wQtekEsuuSSLFy/OWWedlSR57nOfm5NOOinnn3/+rNe/KcEGAADmoUMOOSSrV6/OVVddld133z2HHXZYVq1alfPOO+82weaAAw7I8uXLb37e+vXrc+WVV+bHP/5xHvGIRyRJjjnmmFn/HyZzjg0AAMxDCxcuzNKlS3Pqqafm4Q9/eB74wAfm3HPPzWWXXZaDDjroVtPuvvvuN99fsGBBrrvuurTWdqruqh2xAQCAeerwww/PG9/4xhx++OFZuXJlTj755CxfvnxagWXx4sXZd9998/nPfz5Jcvrppw9d7hY5YgMAADuBHe2eeXusXLkyf/Znf5bDDjsse+21VxYtWnSbZmhbcuqpp+Z5z3te9txzz/zar/3agJVuXbXW5rSAyVasWNFWrVo112UAAMDgLr300ts0+eLWplpHVbW6tbZi02k1RQMAAEZPsAEAAEbPOTYAALAFy05bttnH1j5n7SxWwpY4YgMAAIyeYAMAAIyeYAMAAIyec2wAAGBncOK+Mzy/K7c6yYIFC7Js2bLccMMNOeigg3LaaaflCU94Ql75ylfe6ro0b37zm/P1r389b3vb22a2xhnkiA0AAMxTe+yxR9asWZN169bl9re/fU4++eQcddRROfPMM2813ZlnnpmjjjpqjqqcHsEGAADIypUr881vfjNPf/rT87GPfSw//elPkyTr16/P9773vfzKr/zKHFe4ZYINAADMczfccEM+8YlPZNmyZbnTne6UQw89NJ/85CeTdEdrnvnMZ6aq5rjKLRNsAABgnrruuuuyfPnyrFixIvvvv3+e//znJ8mtmqONoRlaovMAAACYtybOsdnUk5/85PzBH/xBLrzwwlx33XV58IMfPAfVbRtHbAAAgFvZe++988hHPjLPe97zRnG0JnHEBgAAdg7T6J55Nh111FF56lOfepse0nZWgg0AAMxTV1999WYfe8pTnpLW2ixWs2M0RQMAAEZPsAEAAEZPsAEAAEZvsGBTVfepqjWTbj+pqhcNtTwAAGD+GqzzgNba15IsT5KqWpDku0nOHmp5AADA/DVbTdEek+Sy1tp3Zml5AADAPDJb3T0/K8kZUz1QVccnOT5J9t9//1kqBwAAdi7LTls2o/Nb+5y1Mzq/JPnQhz6Ue9/73rnf/e434/PeUYMfsamq2yd5UpL3T/V4a+2U1tqK1tqKJUuWDF0OAACwnT70oQ/lK1/5ylyXMaXZaIr2+CQXttZ+MAvLAgAApmH9+vU56KCDctxxx+X+979/jjjiiFx33XVJkne84x15yEMekoMPPjhPe9rTcu211+aLX/xiPvKRj+RlL3tZli9fnssuu+xW87vsssvysIc9LA95yEPy6le/OnvvvXeS7iKgj3nMY/LgBz84y5Yty4c//OEkyTXXXJMjjzwyBx98cB7wgAfkH//xH3fo/5mNYHNUNtMMDQAAmDvf+MY38oIXvCCXXHJJFi9enLPOOitJ8tSnPjUXXHBBLr744hx00EF55zvfmYc//OF50pOelDe84Q1Zs2ZNDjzwwFvN64QTTsgJJ5yQCy64IHe7291uHr9o0aKcffbZufDCC3PuuefmJS95SVpr+eQnP5m73e1uufjii7Nu3bo87nGP26H/ZdBgU1V7JvnVJB8ccjkAAMC2O+CAA7J8+fIkySGHHJL169cnSdatW5eVK1dm2bJlOf3003PJJZdsdV7nn39+nvGMZyRJfvM3f/Pm8a21vOpVr8oDH/jAPPaxj813v/vd/OAHP8iyZcvyL//yL/nDP/zDnHfeedl333136H8ZNNi01q5trd2ptXblkMsBAAC23e67737z/QULFuSGG25Ikhx77LF5y1vekrVr1+Y1r3lNNm7cuN3LOP3007Nhw4asXr06a9asyV3vetds3Lgx9773vbN69eosW7Ysr3zlK/Mnf/InO/S/zFZ3zwAAwEhcddVV+YVf+IVcf/31Of30028ev88+++Sqq66a8jkPe9jDbm7KduaZZ948/sorr8xd7nKXLFy4MOeee26+853uCjDf+973sueee+boo4/OS1/60lx44YU7VPNsdfcMAABswRDdM2+vP/3TP81DH/rQ3OMe98iyZctuDjPPetazctxxx+Wkk07KBz7wgVudZ/PmN785Rx99dN70pjflyCOPvLlp2bOf/ez8+q//elasWJHly5fnvve9b5Jk7dq1ednLXpbb3e52WbhwYd7+9rfvUM3VWtuhGcykFStWtFWrVs11GQAAcLMtXV9mR8LIpZdemoMOOmi7n7+zufbaa7PHHnukqnLmmWfmjDPOuLkHtO011TqqqtWttRWbTuuIDQAAsMNWr16dF77whWmtZfHixXnXu941q8sXbAAAgB22cuXKXHzxxXO2fJ0HAADAHNmZTgvZ2WzruhFsAABgDixatChXXHGFcDOF1lquuOKKLFq0aNrP0RQNAADmwH777ZfLL788GzZsmOtSdkqLFi3KfvvtN+3pBRsAAJgDCxcuzAEHHDDXZewyNEUDAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGT7ABAABGb9BgU1WLq+oDVfXVqrq0qg4bcnkAAMD8tNvA8//rJJ9srT29qm6fZM+BlwcAAMxDgwWbqrpDksOTHJskrbWfJfnZUMsDAADmryGbov1Skg1JTq2qi6rq76pqr00nqqrjq2pVVa3asGHDgOUAAAC7qiGDzW5JHpzk7a21ByW5JskrNp2otXZKa21Fa23FkiVLBiwHAADYVQ0ZbC5Pcnlr7d/64Q+kCzoAAAAzarBg01r7ryT/WVX36Uc9JslXhloeAAAwfw3dK9rvJTm97xHtW0meO/DyAACAeWjQYNNaW5NkxZDLAAAAGPQCnQAAALNBsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZPsAEAAEZvtyFnXlXrk1yV5MYkN7TWVgy5PAAAYH4aNNj0HtVa+9EsLAcAAJinNEUDAABGb+hg05J8qqpWV9XxU01QVcdX1aqqWrVhw4aBywEAAHZFQwebX26tPTjJ45O8oKoO33SC1toprbUVrbUVS5YsGbgcAABgVzRosGmtfa//+8MkZyc5dMjlAQAA89Ngwaaq9qqqfSbuJzkiybqhlgcAAMxfQ/aKdtckZ1fVxHLe21r75IDLAwAA5qnBgk1r7VtJDh5q/gAAABN09wwAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIzeVoNNdY6uqlf3w/tX1aHDlwYAADA90zli87YkhyU5qh++KslbB6sIAABgG+02jWke2lp7cFVdlCSttf+pqtsPXBcAAMC0TeeIzfVVtSBJS5KqWpLkpkGrAgAA2AbTCTYnJTk7yV2q6s+SfD7Jnw9aFQAAwDbYYlO0qrpdkm8neXmSxySpJE9urV06C7UBAABMyxaDTWvtpqp6U2vtsCRfnaWaAAAAtsl0mqJ9qqqeVlU1eDUAAADbYTq9ov1Bkr2S3FhVG/txrbV2h+HKAgAAmL6tBpvW2j6zUQgAAMD2ms4Rm1TVk5Ic3g9+prX2seFKAgAA2DZbPcemql6X5IQkX+lvJ/TjAAAAdgrTOWLzhCTLW2s3JUlVnZbkoiSvGLIwAACA6ZpOr2hJsnjS/X2HKAQAAGB7TeeIzWuTXFRV56a7QOfhSV45aFUAAADbYDq9op1RVZ9J8pB0weYPW2v/NXRhAAAA0zWdzgOekuTa1tpHWmsfTrKxqp48fGkAAADTM51zbF7TWrtyYqC19uMkrxmuJAAAgG0znWAz1TTTuv5NklTVgqq6qKpc+wYAABjEdILNqqr6y6o6sKp+qar+KsnqbVjGCUku3b7yAAAAtm46web3kvwsyT8meX+SjUleMJ2ZV9V+SY5M8nfbWyAAAMDWTKdXtGvSX4yzqhYk2asfNx1vTvLyJPtsboKqOj7J8Umy//77T3O2AAAAt5hOr2jvrao7VNVeSS5J8rWqetk0nvfEJD9srW2x2Vpr7ZTW2orW2oolS5ZMu3AAAIAJ02mKdr/W2k+SPDnJOUn2T3LMNJ73y0meVFXrk5yZ5NFV9Z7tLRQAAGBzphNsFlbVwnTB5sOtteuTtK09qbX2ytbafq21pUmeleRfW2tH71C1AAAAU5hOsPnbJOuT7JXkc1V1jyQ/GbIoAACAbTGdzgNOSnLSxHBV/UeSR23LQlprn0nymW2sDQAAYFqmfaHNCa21luSGAWoBAADYLtNpigYAALBTm053z7tPZxwAAMBcmc4Rm/OnOQ4AAGBObPYcm6r6+SR3T7JHVT0oSfUP3SHJnrNQGwAAwLRsqfOAX0tybJL9krwptwSbnyR51bBlAQAATN9mg01r7bQkp1XV01prZ81iTQAAANtkOufYHFJViycGqurnqur/DVgTAADANplOsHl8a+3HEwOttf9J8oThSgIAANg20wk2CyZ371xVeyTR3TMAALDT2FLnARPek+TTVXVqkpbkeUlOG7QqAACAbbDVYNNae31VrU3ymHQ9o/1pa+2fBq8MAABgmqZzxCattU8k+cTAtQAAAGyXrZ5jU1UPq6oLqurqqvpZVd1YVT+ZjeIAAACmYzqdB7wlyVFJvpFkjyS/neRvhiwKAABgW0y3Kdo3q2pBa+3GJKdW1RcHrgsAAGDaphNsrq2q2ydZU1WvT/L9JHsNWxYAAMD0Tacp2jH9dC9Mck2SX0zytCGLAgAA2BbT6e75O/0Rm6VJPpjka621nw1dGAAAwHRtNdhU1ZFJTk5yWbrr2BxQVf+r7wIaAABgzk3nHJs3JXlUa+2bSVJVByb5eFzXBgAA2ElM5xybH06Emt63kvxwoHoAAAC22XSO2FxSVeckeV+SluQZSS6oqqcmSWvtgwPWBwAAsFXTCTaLkvwgySP64Q1J7pjk19MFHcEGAACYU9PpFe25s1EIAADA9ppOr2gHJPm9dN093zx9a+1Jw5UFAAAwfdNpivahJO9M8tEkNw1bDgAAwLabTrDZ2Fo7afBKAAAAttN0gs1fV9VrknwqyU8nRrbWLhysKgAAgG0wnWCzLMkxSR6dW5qitX4YAABgzk0n2DwlyS+11n42dDEAAADb43bTmObiJIuHLgQAAGB7TeeIzV2TfLWqLsitz7HR3TMAALBTmE6wec3gVQAAAOyArQab1tpnZ6MQAACA7bXZYFNVV6Xr/ew2DyVprbU7DFYVAADANthssGmt7TObhQAAAGyv6fSKBgAAsFMTbAAAgNETbAAAgNETbAAAgNETbAAAgNETbAAAgNETbAAAgNETbAAAgNETbAAAgNETbAAAgNETbAAAgNEbLNhU1aKq+vequriqLqmqPx5qWQAAwPy224Dz/mmSR7fWrq6qhUk+X1WfaK19acBlAgAA89Bgwaa11pJc3Q8u7G9tqOUBAADz16Dn2FTVgqpak+SHSf65tfZvU0xzfFWtqqpVGzZsGLIcAABgFzVosGmt3dhaW55kvySHVtUDppjmlNbaitbaiiVLlgxZDgAAsIualV7RWms/TvKZJI+bjeUBAADzy5C9oi2pqsX9/T2SPDbJV4daHgAAMH8N2SvaLyQ5raoWpAtQ72utfWzA5QEAAPPUkL2ifTnJg4aaPwAAwIRZOccGAABgSIINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAweoINAAAwervNdQEAAMA0nbjvFh67cvbq2Ak5YgMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIyeYAMAAIzeYMGmqn6xqs6tqkur6pKqOmGoZQEAAPPbbgPO+4YkL2mtXVhV+yRZXVX/3Fr7yoDLBAAA5qHBjti01r7fWruwv39VkkuT3H2o5QEAAPPXrJxjU1VLkzwoyb9N8djxVbWqqlZt2LBhNsoBAAB2MYMHm6raO8lZSV7UWvvJpo+31k5pra1ora1YsmTJ0OUAAAC7oEGDTVUtTBdqTm+tfXDIZQEAAPPXkL2iVZJ3Jrm0tfaXQy0HAABgyCM2v5zkmCSPrqo1/e0JAy4PAACYpwbr7rm19vkkNdT8AQAAJsxKr2gAAABDEmwAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDRE2wAAIDR222uCwAAgJ3CiftOPf6A/We3DraLIzYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDoCTYAAMDo7TbXBQAAALdY+oqPb/ax9YtmsZCRccQGAAAYPcEGAAAYvcGCTVW9q6p+WFXrhloGAABAMuwRm3cnedyA8wcAAEgyYLBprX0uyX8PNX8AAIAJc36OTVUdX1WrqmrVhg0b5rocAABghOY82LTWTmmtrWitrViyZMlclwMAAIzQnAcbAACAHSXYAAAAozdkd89nJDk/yX2q6vKqev5QywIAAOa33YaacWvtqKHmDQAAMJmmaAAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOgJNgAAwOjtNtcFcGvLTlu22cfWPmftLFYybtYjAMD8ItgAADBvLH3Fxzf72PpFs1gIM05TNAAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPQEGwAAYPR2m+sCgOEtfcXHN/vY+kW/ufknnnjlANUAAMw8wQYAAHYBy05bNuX4tc9ZO8uVzA3BBgCAYZ247xYe0zqAmeEcGwAAYPQcsZkLW9prccD+s1fH2FmPAAD0HLEBAABGb9BgU1WPq6qvVdU3q+oVQy4LAACYvwZrilZVC5K8NcmvJrk8yQVV9ZHW2leGWubOZMvd685iISNnPQIAMB1DnmNzaJJvtta+lSRVdWaS30gyL4IN7Ao2121kMn+6jgRg+ja3Q3JLOyP91jBTqrU2zIyrnp7kca213+6Hj0ny0NbaCzeZ7vgkx/eD90nytRks485JfjSD8xuCGmfGGGpMxlGnGmeGGmeGGmfGGGpMxlGnGmeGGmfGfK3xHq21JZuOHPKITU0x7jYpqrV2SpJTBimgalVrbcUQ854papwZY6gxGUedapwZapwZapwZY6gxGUedapwZapwZary1ITsPuDzJL04a3i/J9wZcHgAAME8NGWwuSHKvqjqgqm6f5FlJPjLg8gAAgHlqsKZorbUbquqFSf4pyYIk72qtXTLU8jZjkCZuM0yNM2MMNSbjqFONM0ONM0ONM2MMNSbjqFONM0ONM0ONkwzWeQAAAMBsGfQCnQAAALNBsAEAAEZvlwk2VfXzVXVmVV1WVV+pqnOq6t797Zyq+mZVXVpV76uqu6pRjQPXeGNVrZl0e0U//vZV9ea+9m9W1ceqan81qlGNdfUmw8dW1Vv6+ydW1Us3eXx9Vd1ZjWpU485bI8y2Ia9jM2uqqpKcneS01tqz+nHLk9w1ybuS/EFr7aP9+EclWZLkB2pU44Cua60tn2L8nyfZJ8m9W2s3VtVzk3y4qg5prd00uyWqcYaoEWAr+t/v2pm/W9Q4M+ayxl3liM2jklzfWjt5YkRrbU2SeyU5f2JDtx9/bmttnRrVONuqas8kz03y4tbajUnSWjs1ydVJHjuXtU1Q48xQI0BSVUura0HxtiQXJvmjqrqgqr5cVX/cT7NXVX28qi6uqnVV9Uw1qnF77RJHbJI8IMnqbRg/F9Q4M8ZQY5LsUVVrJg2/NsmlSf6jtfaTTaZdleR+ST41W8X11Dgz1DgzNq3xjrn1tc9eXFVHTxq+2+yUdStqnBlqnBljqDFJ7pNuJ8qHkjw9yaFJKslHqurwdC0rvtdaOzJJqmpfNapxe+0qwQZ2Nrdp+lNVByeZqn/1mp2SbkONM0ONM+NWNVbVsUlWTHr8r1prb5z0+PrZK+1mapwZapwZY6gxSb7TWvtSVb0xyRFJLurH752utcV5Sd5YVX+R5GOttfPUqMbttas0RbskySHbMH4uqHFmjKHGzflmkntU1T6bjH9wur3kOwM1zgw1AnSu6f9Wkte21pb3t3u21t7ZWvt6ut/vtUleW1WvVqMat9euEmz+NcnuVXXcxIiqeki6H+6HV9WRk8Y/rqqWqVGNs621dk2S05L8ZVUtSJKq+q0kG5N8YS5rm6DGmaFGgNv4pyTPq6q9k6Sq7l5Vd6mquyW5trX2niRvTLdzRY1q3C67RFO01lqrqqckeXN13ZhuTLI+yYuSPLEf/+Yk1yf5cpIT1KjGgW3a9vmTrbVXJHllkjck+VpV7ZFkQ5LDWmtTNQlSoxrnU43ALqy19qmqOijJ+VWVdB2VHJ3knkneUFU3pfvt/t9qVOP2Kr9fMDeq6ueTfDLJ21prp8x1PVNR48xQIwAMT7ABAABGb1c5xwYAAJjHBBsAAGD0Rh9sqrvS6W2uLl9Vf1dV95ujmn6/uquvnj4D8zqnqhb396/e8eqmXMbVmwzfZp1W1YlV9dJJwy+tqq9Wd+XYi/velGa6rsVV9btbmWZpVf3mNOY15ftkLlTVk/qOD6Z6bJDXeJNlbHW97gyq6saqWjPptrSqLqqq5f3ju1XVNTXpAnRVtbqq5rInmJ3GWF7nyarq56vqzKq6rKq+0n//3buqrtvkvTDj3zdDqKpHVtXHZmE5rareNGn4pVV14qTh3+q/qy/p1+tL+/HvrqpvT1qvXxywxv2q6sNV9Y3+9f3rqrp9v46u7D/bX6uqz1XVEyc978Sq+m5f37qqetKANd44aTnvr6o9t1T7pOcd2tf9tf538e8mnsvmVXa8VBMAABFwSURBVNWL5nI9zeS22hBmantgtrZ/+u+hf5g0vFtVbZj4DqyqY/tpHjNpmqf0454+U3WMPthsTmvtt1trX5mjxf9ukie01p69ozNqrT2htfbjGahpxlTV7yT51SSHttYekOTwDHNRv8Xp1uWWLE2y1WCzM2mtfaS19ro5LGE663VncN2kPvCXt9bWJ/likof3jx+c5GsTw1W1V5JfSnLxXBS7PapqyJ4pt+l1rs7tNjc8tKqqJGcn+Uxr7cDW2v2SvCrJXZNctsl74e9nq67N1Tqb62YafprkqVV1500fqKrHp+s18ojW2v3Tda965aRJXjZpvT580+fPhP61/WCSD7XW7pXk3uku2Pdn/STntdYe1Fq7T5LfT/KWyRs/6S40uTzJM5K8a8B1P/Gd84AkP0vyO1urvarumuT9Sf6wr/+gdJ1wbHp9qBm3E74Pt9WLkmxTsKm+a/oZMmPbaiTprmHzgOp62Ey67cTvbjLN2iRHTRp+Vmb4N3vMH4jJdquq06rqy1X1garas6o+U1UrkqSqjqqqtf1emL8YspCqOjndxtVHquolVfWhvq4vVdUD+2keMWkP2UVVtU9V/UK/x2dib9HKftr1U/1YzbFXJfnd1tpPkqS1dmVr7bQBlvO6JAf26+QN/W1d/1o+c9I0K/tpXtzvmTivqi7sb4P8UE+lf61XV7dX9Ph+3OP6Oi6uqk/3446tqrf09w+oqvOr6oKq+tNZKnWr67WqTqr+wllV9Wv9e3Nn+L74Qm4JNg9PcnKSiStvH5rkwtbajUMX0b/PJvbMrquq06vqsVX1her26h7a377Yf8a/WFX36Z97bHV7gz+a5FMDlrnp6/yy/n325ar640n/x6VV9bYkF6b7LE0e/sUB69vUo5Jc31o7eWJEa21Nkv+cxRo2a4p19c6qWtV/3v940nSP698bn0/y1Fkq74YkpyR58RSPvTLJS1tr30uS1trG1to7ZqmuCY9OsrG1dmpfw43pan1eNtmw7V/zP0nywk1n0lq7NN3/Ohu/ieel65p2s7VXd7ThBUlOa62d3z/eWmsfaK39YIiipngf/tEUn+u9qurj/e/Ouknf6+ur6i+q6t/72z1nqKaXV9Xv9/f/qqr+tb//mKp6T1W9fdPPSj/93ZKcW1Xn9uOO6H8PL+y/IyeugbK+ql7df6aeMUM1z+i22pCqau+q+nS/XtZW1W/04yfeC+/o1+2nqg8VVXVI//qfn+49Ols+kWTieoJHJTljk8fPS3JoVS3sX997JlmTmdRaG/Ut3R77luSX++F3JXlpks8kWZHug/MfSZaku27PvyZ58sA1rU/3xfs3SV7Tj3t0kjX9/Y9Oqnfvvq6XJPk//bgFSfaZPK/+/tUD1Xv1JsNLk6zbZNyJ/XrdJ8n/zOJru66//7Qk/9yvm7v2r+kvJHlkko9Nes6eSRb19++VZNXm/qcB6r1j/3ePJOv6Ov8zyQGbPH5skrf09z+S5Lf6+y8Y6jXejvW6Z5JL0m1sfi3JgbPxmm9S543pvvDWJDl7Uu3f6u+fkeS+Sc7t35f/J8mfzOJ784Yky9LtIFqd7runkvxGkg8luUOS3frpH5vkrEmv/+UT74dZep2PSLfhW329H0t3pHVpkpuSPGzSc24enuXX+/fT7Zmf6v+4btJ7YU2SlXNQ36brauLzvCDd780DkyzqP/P36tf1+zLp+2nA2q7u32/rk+yb7rv6xP6x/06y72ae9+4k3560Xk+f5df2ov6xj20yfnmSS/v7J6YLZkny0CTfS9+j6xDrsf+7W5IPp7vGxpZqf2C6ozm/MRfvwy18rp+W5B2TnrNv/3d9btnO+K2Zem/2tby/v39ekn9PsjDJa5L8r6k+K5Pqmdi+uXOSzyXZqx/+wySvnjTdywdYl+szQ9tqs/B+vMOk9fTN/jVfmu53aHn/2PuSHN3f/3KSR/T335CBt38m6u0/Ex9I9124JpO20dJv+yT5y3TXHXx2/x55d5Knz1QdO8Me2Jnwn621iStlvyfJr0x67CHpmjZsaK3dkOT0dB/82fArSf4hSVpr/5rkTlW1b7q9zn/Z77FY3Nd1QZLnVtcuellr7apZqnEqm+sDvKX7MM1FH+G/kuSM1tqNrdsT9tl0r+2mFiZ5R1WtTdc8YDbPs/r9qro4yZfS7ek+PsnnWmvfTpLW2n9P8Zxfzi17NP5hiseHNuV6ba1dm+S4dKHnLa21y+agtslN0Z6SJK1rjnb76q65ct90oeuCdBs8D0/XVG22fLu1tra1dlO6EPjp1n17r033g7NvkvdX17b5r5Lcf9Jz/3kz74ehHNHfLkq3l/e+6Ta+k+Q7rbUvTZp20+GdwaZN0c6bozomr5v/r6ouTLdO75/uu+a+6d4X3+jfC++ZrcJadwT979NtiG+LyU3RhmqSs7nfjS2Nn+zF1V1g9o1Jntmv2yFMXMh2VbqdPO/cSo1zdb2Miffh5j7Xa5M8tj86s7K1Nrnp4RmT/h42Q/WsTnJIVe2Trlnk+el2LK9MF3Sm+qxs6mH9+C/0r8Fzktxj0uP/OEO1TmVn31arJH9eVV9O8i9J7p5uR2TSfd9MHPFYnWRpX/vi1tpn+/Gztm3RWvtyut+/o5Kcs5nJzkzXBO1Zue0RnR22qwSbTb9cJg8Pce7HdE217Na68yt+O92e/S9V1X1ba59LF7i+m+Qfam5Pjr0iyc9tMu6OSX7U/3heU1W/NMs1Tfd1fHGSH6Q7/2JFkttvefKZUVWPTLdX/rDW2sHpvsAvzvR++ObyYlJbWq/L0r0X7jZLtUzX+UmenuT7/QbOl9IFxEP7+7Plp5Pu3zRp+KZ0e9j+NMm5rWuv/+vp9mBNuGZWKrxFJXntpA3Ye7bW3rmZWma7tgmXJDlkjpY9XdckXRPSdEdFHtNae2CSj+eW13cuP89vTvL8JHtNGrczrNdL0n0f36yq7pBuB9BUO00elOTSScN/1b9vVw4caifvTPm91trPsvXa52L9TnxGp/xct9a+3te0Nslrq29W3Gubub/dWmvXpzv68dx0O5fOS3e0/8B0R1s391mZrNLt8Jn4X+7XWnv+pMeH/F7a2bfVnp2u1dEhrTvX7Ae5ZR1O/h26Md1vz1yG7qRrifLGbCa0tNb+PckD0h2t+/pML3xXCTb7V9XEnoejknx+0mP/luQRVXXn6k46OyrdXunZ8Ll0b8iJDd8ftdZ+UlUH9nt6/yLdnqH7VtU9kvywdW2f35nuBM850Vq7Osn3qz95s6rumORxuWW9vjbJW/sv91TVHao/p2SGXZVbTsD8XJJnVtWCqlqS7ovl3zeZJun2kn+/34t+TLpDxbNh33RN9K6tqvum2/u0e7r33gHJzetxU19It9ci6d8rs2Cr67V/P74k3QbG46vqobNU23R8IV2APb8fPj9ds4r/ajtXRxv75pYTJ4+dg+VPfp3/Kd05ARNt1u9eVXeZg5q25F+T7F5Vx02MqKqH5NZ7bXcWd0i3oXVldSePP74f/9UkB1TVgf3wUVM9eSj9UcD3pQs3E16b5PX9Uc5U1e4T50PMok8n2XNiI7D/LX5TuiYo106esD+/4Y+SvHWWa9yczdbeH9l+S5LnTP6OrKqjJ9b3wKb8XFfV3ZJc21p7T7oNzMnbE8+c9Pf8zJzPpQswn0sXbH4nXVOkzX1Wklt/R30pyS9Xf95PdedK33sG69ta7Tvzttq+/TKvr6pHZSvfif3v4JVVNdF6abY7R3hXumbha7cwzSvTna8944bskWc2XZrui+Vvk3wjydvT7SFNa+37VfXKdO3wK8k5rbUPz1JdJyY5tT98eG26Q6tJ8qL+zXljkq+kO9nqWUleVlXXp2unOJtHbPasqssnDf9lv/y31i1diP7xpOZIb0/X3vSCvt7r033Rz6jW2hXVnYy9Lt06+nJuOQry8tbaf1XVFUlu6JuAvTvJ25KcVVXPSPeaz9be50+m60Hny+maR30pyYZ0zdE+WN2J9z9M10vIZCckeW9VnZDkrNkodGvrNd3eoH9Of8JxVT0/ybur6iGttY2zUeNWfCFd066Jk3W/329szGYztOl4fZLTquoP0m20z6opXuf3Jjm/qpLuO+bodN9BO4XWWquqpyR5c3XdoW9Mtxf4Rek7QZg0+btaayfNQZlJktbaxVV1Ubq99d9K955Ma21jv5Pn41X1o3Q7gx4wy+W9KZNOvG+tndNvUP5LdS9+S7fhMeENVfV/Jw0f2h+pmDGTXtu3VdUfpdupek66DZvD0nVacVG6c/t+mOT3W2ufnskattdWak9r7QdV9awkb+x3FtyUbkP5g7NQ26eq6qDc9nN9z3Sv603pfp//96Sn7V5V/9b/HzMZvM9Ld57j+a21a6pqY7re7qb8rPROSfKJqvp+a+1RVXVskjOqavf+8f+bZMb36E/hxOzc22qnJ/loVa1KFxa/Oo3nPDddD4LXpgvAs6a1dnmSv97KNJ8Yavk1XFNVAAB2BlW1PsmK1tqP5roWGMqu0hQNAACYxxyxAQAARs8RGwAAYPQEGwAAYPQEGwAAYPQEG4B5rKqunqPlHttfb2NH57O4qn53JmoCYNwEGwDmwrFJdjjYJFmcZPBgU1W7ynXfAHZZgg0AqapHVtVnq+p9VfX1qnpdVT27qv69qtZW1YH9dO+uqpOr6rx+uif245f24y7sbw+fNO+X9/O4uJ/v05OsSHJ6Va2pqj02qeWeVfUv/fQXVtWBVbV3VX26H15bVb/RT/669BfwrKo39M9/WVVdUFVfrqo/njTfP6qqr1bVP1fVGVX10n788qr6Uj/92VX1c/34z1TVn1fVZ5P8n6r6dlUt7B+7Q1WtnxgGYO7ZAwXAhIOTHJTkv9NdJfzvWmuHVtUJSX4vyYv66ZYmeUSSA5OcW1X3THfF+F9trW2sqnslOSPJiqp6fJInJ3loa+3aqrpja+2/q+qFSV7aWls1RR2nJ3lda+3sqlqUbifcz5I8pbX2k6q6c5IvVdVHkrwiyQNaa8uTpKqOSHKvJIcmqSQfqarD011R/GlJHpTut+/CJKv75f19kt9rrX22qv4kyWsm/a+LW2uP6Oe9NMmRST6U7grkZ7XWrt+O9QzAAAQbACZc0Fr7fpJU1WVJPtWPX5vkUZOme19r7aYk36iqbyW5b5JvJ3lLVS1PcmOSe/fTPjbJqa21a5OktfbfWyqgqvZJcvfW2tn99Bv78QuT/HkfUm5Kcvckd51iFkf0t4v64b3TBZ19kny4tXZdP7+P9n/3TRdePttPf1qS90+a3z9Ouv93SV6eLtg8N8lxW/pfAJhdgg0AE3466f5Nk4Zvyq1/Lza9snNL8uIkP0h31Od2STb2j9UU029JbWb8s5MsSXJIa+36qlqfZNFmnv/a1trf3mpk1Yu3oYbJrpm401r7Qt/k7hFJFrTW1m3nPAEYgHNsANhWz6iq2/Xn3fxSkq8l2TfJ9/sjOcckWdBP+6kkz6uqPZOkqu7Yj78q3VGUW2mt/STJ5VX15H763fvn7pvkh32oeVSSe2xmPv/UL2/v/vl3r6q7JPl8kl+vqkX9Y0f2y7syyf9U1cr++cck+Ww27+/TNbM7datrCYBZ5YgNANvqa+k2/u+a5Hf682reluSsqnpGknPTH+lorX2yb562qqp+luScJK9K8u4kJ1fVdUkOm2gi1jsmyd/257tcn+QZ6c67+WhVrUqyJslX+/lfUVVfqKp1ST7RWntZVR2U5PyqSpKrkxzdWrugPyfn4iTfSbIqyZX98p7T17JnunOLnruF//30JP8vXbgBYCdSrW1LCwEA5rOqeneSj7XWPjDXtWyrqtq7tXZ1H2A+l+T41tqF2ziPpyf5jdbaMYMUCcB2c8QGgPnilKq6X7pzc07bjlDzN0ken+QJQxQHwI5xxAYAABg9nQcAAACjJ9gAAACjJ9gAAACjJ9gAAACjJ9gAAACj9/8DUoZJFTnlZRgAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1008x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plot data\n",
    "df.plot.bar(xlabel='Impact category', ylabel='Impact score', figsize=(14,8))\n",
    "\n",
    "# rotate xticks for horizontal impact category names\n",
    "plt.xticks(rotation=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Some of the bars are not visible because the impact categories have wildly different units (and hence magnitudes). Let's normalize each category by the largest score. While we're at it, let's also increase the font size and add a grid."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15, 16,\n",
       "        17, 18]),\n",
       " <a list of 19 Text major ticklabel objects>)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA04AAAH7CAYAAAD7OM9uAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzde5xVdb3/8deHi4IICIJoIqA//alcUgvLuyOZ2tHwePJWeQTMSzdTT3mio0cR7fZLzTI7almoddS8a5yUCndq5knLlFDzEkiaKSIgKBeR7++PvcGZYe9Ze2b2bZjX8/GYx8xe37XWfs/ew8CbtdZ3RUoJSZIkSVJpPeodQJIkSZIancVJkiRJkjJYnCRJkiQpg8VJkiRJkjJYnCRJkiQpg8VJkiRJkjL0qneAWhoyZEgaNWpUxfb35ptv0q9fv4rtrxrMWBldISN0jZxmrAwzVoYZK6MrZISukdOMlWHGyuiuGf/whz+8llIausFASqnbfLz//e9PlXTfffdVdH/VYMbK6AoZU+oaOc1YGWasDDNWRlfImFLXyGnGyjBjZXTXjMCjqUiX8FQ9SZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDN3qPk6SJElSe6xcuZKFCxeycuVK1qxZU7H9Dhw4kKeeeqpi+6uGjS1j79692WqrrRgwYECHnsviJEmSJBWxdOlSXnnlFYYOHcrWW29Nr169iIiK7HvZsmX079+/Ivuqlo0pY0qJFStW8NJLLwF0qDx5qp4kSZJUxGuvvcbw4cMZNGgQvXv3rlhpUu1FBJttthnbbrstr776aof2YXGSJEmSili9ejV9+/atdwxVUN++fXn77bc7tK3FSZIkSSrBo0wbl868nxYnSZIkScpgcZIkSZKkDDWfVS8ihgNfBsYDuwF9ge1TSvPL2LYPcCFwArAF8Cfgyyml+6sWWJIkSWpl1NSZdX3++d84vCbPM2PGDKZMmcK8efMYNWpURfY5atQompqamDFjRkX2Vyv1OOK0I3AssBh4oJ3bXgOcApwHHAG8DNwbEbtXNKEkSZIkDj/8cH73u9+xzTbb1DtK3dXjPk73p5SGAUTEycAh5WwUEbsBnwBOSin9uLDsN8BcYDowsTpxJUmSpO5p6NChDB06tN4xGkLNjzillNZ2cNOJwNvATc32tQa4ETg0IjatQDxJkiRpo/Too48SETz44IPrl11++eVEBOeee+76Zc8++ywRwb333suMGTOICObPn79+fNSoUZxwwgnceOON7LrrrvTr14/x48e32O863/nOdxg1ahR9+vRh/PjxPPBAe084axxdaXKIMcC8lNJbrZbPBTYhfwqgJEmSpCLe9773scUWWzB79uz1y2bPnk3fvn03WNazZ0/22Wefkvt64IEHuOSSS7jwwgu56aabeOeddzjiiCNYsmTJ+nWuueYazjzzTA466CDuuOMOJk+ezMc//nEWL15cnW+wyiKlVL8nz5+q9wPKmBwiImYBA1JKe7VafjDwS+CAlNIGFTYiTgVOBRg2bNj7b7zxxk5lbsodWXIs13Rnp/ZdDcuXL2fzzTdf//j0F04vue7lIy+vRSTA17FSusLraMbqaP4z2Sg/j63VM+Pke94sOTbjsH7rv279Z7tRdIWfya6QsTX/3HRMV3ivq5Vx4MCB7Lhj8f+bH/fV+s5NNuecAzq03fHHH8+yZcuYOXMma9euXX/06Morr2TBggVsk/7B8Z+Zyvy//Z2Hf34dM266iyn/No15D/+cLUfnn3Ps2LG88cYbPP744wwaNAiAP/7xjzQ1NfHDH/6QY489lrVr1zJmzBh22WUXbr/99vXPf+uttzJlyhQ+8YlPcOWVV3boe+i/7LmSY8v6Zx9Lee6551i6dGnJ8YMOOugPKaXxrZfX4xqnjgqgWMtr8y5WKaWrgasBxo8fn5qamjqXIld6qNP7roJcLtcy17Wl161p/lzpIV/HdsiVHmqY1zFXesiMHdfiZ7JRfh5bqWvGe0rPdtX8+Tb4s90ocqWHGiZvrvRQw2RsxT83HZQrPdQw73Wu9FBnMj711FP079+/w9tXU0dzffjDH2bq1Kn07t2bJ598kqVLl3Luuefy4x//mD/96U9ss9vW5H73B046bsPpA9Y9Z0Sw9957M2LEiPVje+2VP7axcOFC+vfvz4IFC3jppZeYPn16i6wnnHACp5xyCr179+74a7us9FA5++zTpw977LFHu5+2KxWn14ERRZYPajYuSZIkqYQJEyawatUqHnroIR577DF22203hg0bxn777cd9993HiD4f4JWFizho3z3b3M/gwYNbPN500/x0AytXrgTg5ZdfBmDYsGEt1uvVqxdbbrllpb6dmupKxWkucFREbNbqOqfRwGqg9DE7SZIkSYwbN44hQ4Ywe/ZsHnvsMSZMmADkC9XPfvYzthsQbLJJb/bdc7dOPc+66ctfeeWVFsvXrFnDokWLOrXveulKk0PcBfQGjlm3ICJ6AccBs1JKq+oVTJIkSeoKIoIDDzyQX/7ylzzwwAMtitNjjz3G7b+YzQf3GMtmfft26nmGDx/Odtttx89+9rMWy2+99VbWrFnTqX3XS12KU0QcHRFHA+8vLPpIYdmBhfGREbEmIs5bt01K6U/kpyK/LCJOjogPkZ+KfHvg/Bp/C5IkSVKXNGHCBH7/+9/z1ltvsf/++wP5GfcGDBjAfQ89ykH7bDAvQrv16NGD888/n3vvvZcpU6Zw7733csUVV/DFL36RAQMGdHr/9VCvU/VubvX4+4XPvwGayE/40JMNi90U4KvARcAWwOPAYSmlP1YtqSRJktTK/G8c3qntly1btn4ig7mvzS253pghYzr1PMUcdNBBAIwfP359ienRowcHHHAAd911Fwft0/b1TeX61Kc+xfLly7n00ku54YYbGDt2LDfeeCMnnHBCRfZfa3UpTimlrJnw5lNktryU0grg3wofkiRJktpp1113pdgtie688074+2Mtlk0+biKTW82w1/xmuM0V2+cZZ5zBGWecUdb2ja4rXeMkSZIkSXVhcZIkSZKkDF1pOnJJ6t6mDVz/ZRO8e8PH7Yvd4k6SJFWSR5wkSZIkKYPFSZIkSZIyWJwkSZIkKYPFSZIkSZIyWJwkSZIkKYPFSZIkSZIyWJwkSZIkKYP3cZLUaaOmziw5Nr9PDYNIklQrze6t1xH9m309pkPPv7RDzztjxgymTJmy/vHmm2/ODjvswCmnnMJxxx3He0Z+gFM+cRTf//pXim5/zTXXcPLJJ3PffffR1NTUoQxdlcVJkiRJ6mZuvvlmhg8fzhtvvMHNN9/M6aefzquvvspHDtqHm+6axWUXfIlNNum9wXbXXXcdI0eO5MADD6xD6vryVD1JkiSpm9l9993Za6+9OOSQQ/jBD35AU1MTl112GZOO+SivL1nKzF8/sME28+fP54EHHuDEE08kIuqQur4sTpIkSVI3t+eee7Js2TL2/+AeDN5iINffuuFp+Ndffz0pJU488cQ6JKw/i5MkSZLUzc2bN4+ePXuyeb++HH/kIcz89YO8vrjldVQ/+clP2Geffdhxxx3rlLK+LE6SJElSN/POO++wZs0aFi9ezFVXXcVtt93GEUccwWZ9+zLpmI+yevXb3HTXrPXrP/yHJ3jmmWeYNGlSHVPXl5NDSJIkSd3MLrvssv7rHj168MlPfpLLLrsMVr7AB/YYy647bc/1t87kM5OOAeC6W2ay6aabcuyxx9Yrct1ZnCRJkqRu5vbbb2f48OH079+fkSNH0qdP4f4hf38BgBOPPoKvfP1ynpu3gBHbbsNNd83iyCOPZIsttqhj6vryVD1JkiSpmxk7dizjx49n5513frc0NfOvHzucHj16cP2tM5n56wd4fcnSbn2aHnjESZIkSVIr226zFR/a7wP85LZf8MRTzzFs6JYccsgh9Y5VVx5xkiRJkrSBScccwV9feJE7783xyaM+Qq9e3fuYS/f+7iVJkqSOmLY0e502LFu2jP79+wMw97W5JdcbM2RMp56nM/7lIxMY0H9z3li2nBOPPqJuORqFxUmSJEnqJiZPnszkyZPLWrdv3z4sffr+6gbqQjxVT5IkSZIyWJwkSZIkKYPFSZIkSZIyWJwkSZIkKYPFSZIkSZIyWJwkSZIkKYPFSZIkSZIyWJwkSZIkKYPFSZIkSZIyWJwkSZIkKUOvegeQJEmSuppx146r6/PPmTSnrs+/Ti6XI5fLcd5559Gjx8Z9TGbj/u4kSZIkVU0ul+OCCy5g7dq19Y5SdRYnSZIkScpgcZIkSZK6iWnTphERPPvssxx++OFsvvnmjBw5kunTp7c4arRy5SrOOv9ixk44hs132petd/8wH/3oR3n66adb7OuCCy4AoHfv3kQEEdHm8y9cuJCPf/zjDBgwgEGDBjFlyhTuuusuIoJcLrd+vVmzZvFP//RPbLPNNmy22WaMHTuWSy65hHfeeafF/v779l+wxyEfZ/Od9mXgLgcwbtw4rrrqqgq8UhvyGidJkiSpmznqqKOYMmUKZ511FnfffTfnn38+2223HVMO3R2AVatXs+zNtzj3jJPZZqshvL5kKd+/8R722msvnn76abbeemtOPvlkXnzxRa655hoefPBBevbsmfm8//Iv/8KcOXP4+te/zo477sitt97K6aefvsF6f/3rX/nQhz7E6aefTp8+fXj00UeZNm0aCxcu5BtfOA6AB3//GCecfi5fOOl4vnXumaxNa3l64RqWLFlS2RerwOIkSZIkdTNf/OIXmTJlCgAHH3wws2fP5oYbblhfnAYO6M8PLz5v/frvvPMOh378MwwbNowbbriBs846i+HDhzN8+HAAPvjBD9KrV9vVYtasWTz44IPcdNNNHHvssQAceuihTJw4kQULFrRY99Of/vT6r1NK7L///qxevZqLL76Yr33+GHr06MHDf5zDFgP7c9n0s9eve8h79ujEq9I2T9WTJEmSupnDDz+8xeOxY8duUF5+dtcsPnjEiWyx6wH0GrEn/fr1Y/ny5fzlL3/p0HM+/PDD9OzZk6OOOqrF8qOPPnqDdV9++WVOO+00Ro4cySabbELv3r0599xzWbJkCa++9joAe+42hsVL3uCE08/h57+8nyVLl3UoV7ksTpIkSVI3M3jw4BaPN910U1auXLn+8d2zfsNxn5nKrjtuz39/72v878+v45FHHmHo0KEt1muPl19+mUGDBtG7d+8Wy4cNG9bi8dq1a5k4cSI///nPOffcc5k9ezaPPPII55xzDgArV60G4MC938/NV/0//vb3Vzjq5C8x9L0f4uCDD+aJJ57oUL4snqonSZIkqYUb75rFjqO2Y8ZlF6xf9vbQsbz++usd3uc222zD4sWLefvtt1uUp1deeaXFes8//zyPPvoo119/PSeccML65XffffcG+zz6iIM5+oiDWf7mW+QeepQv/7+rOeyww3jxxRcrfl8pjzhJkiRJauGtFSvp1avlZA/XX3/9BrPabbrppgCsWLEic5977bUX77zzDrfffnuL5TfffHPL537rLYAW5ertt9/mpz/9acl9b95vM4748AGcdtppvPzyyyxatCgzT3t5xEmSJElSC4c17c0d99zHWedfzBEH788f5jzFd2fcyhZbbNFivdGjRwNwySWX8JGPfISePXsyfvz4ovs85JBD2G+//Tj11FN57bXX2HHHHbnlllt4/PHHAdYfIdp1110ZOXIk55xzDj179qR37958+9vf3mB/533rv3hl4SIO2ndP3jNsCC++/Crf/e6P2H333Rk6dGglXw7A4iRJkiS125xJczq1/bJly+jfvz8Ac1+bW3K9MUPGdOp5OuqUT/4Lf/v7K/zopju56ie3sefuo7n77rs3mNjhiCOO4LOf/Szf//73mT59OiklUkol93vbbbdx+umn8+Uvf5mePXsyceJELrzwQiZPnszAgQMB2GSTTbjjjjv4/Oc/z4knnsjgwYM56aSTGDFiBKeccsr6fX1wj7F890c3cNa0S3h9yVK22nIwh3zkcC688MKqvCYWJ0mSJKmbmDZtGtOmTdtg+YwZM/Jf/P0xIH/056Ivf46Lvvy5d1d6zx7Mnz+/xXY9e/bkiiuu4Iorrijr+YcOHcqNN97YYtnnPvc5NttsM3bZZZf1y3bffXcefPDBDbY/+eST12c8/OD9Ofzg/VuuUMXpyC1OkiRJkmpixowZLF26lDFjxrB69WruuecerrzySs4+++z110s1KouTJEmSpJro168fl112Gc8//zyrVq1i++2352tf+xpnn3129sZ1ZnGSJEmSVBPHHHMMxxxzTL1jdIjTkUuSJElSBouTJEmSVEJbM8Sp6+nM+2lxkiRJkorYZJNNyrqxq7qOFStWtLixbnt4jZO0kRp37biSY52994RU0rSB679sAsgVHmw/ovZZJKmThgwZwosvvsiQIUPo378/vXr1IiLqHUsdkFJixYoVvPTSSwwbNqxD+7A4SZIkSUUMHDiQTTfdlIULF7Jo0SLWrFlTsX2vXLmSPn36APCP5f8ouV6PhTU+QWzJq6XHlj5Vuxxt6WDG3r17M2zYMAYMGNChp7U4SZIkSSX06dOH7bbbruL7zeVy7LFH/matx157bMn1an6WyLS92hhbWrscbalTRq9xkiRJkqQMFidJkiRJymBxkiRJkqQMFidJkiRJymBxkiRJkqQMFidJkiRJymBxkiRJkqQMFidJkiRJymBxkiRJkqQMFidJkiRJymBxkiRJkqQMFidJkiRJymBxkiRJkqQMFidJkiRJymBxkiRJkqQMNS9OEbFdRNwSEUsj4o2IuC0iRpS57YiIuDYiFkTEWxHxTERcFBH9qp1bkiRJUvfVq5ZPFhGbAbOBVcAkIAEXAfdFxHtTSm+2sW0/4FdAb+A/gQXAnsAFwE7AcdVNL0mSJKm7qmlxAk4BdgB2Tik9BxARTwDPAqcBl7ax7b7kC9KhKaVZhWX3RcRg4EsRsVlK6a3qRZckSZLUXdX6VL2JwMPrShNASmke8FvgyIxtNyl8fqPV8iXkv4+oVEhJkiRJaq7WxWkM8Ociy+cCozO2/RX5I1PfjIjREbF5REwAzgCubOs0P0mSJEnqjFqfqjcYWFxk+evAoLY2TCmtjIj9gFvJF611fgh8vtR2EXEqcCrAsGHDyOVy7YzcUlMbY53ddzUsX7687Fy1zN/UxpivY/maOrhdo7zGjZKjqY2xRskIHXu/a52/qQPb1Dpj8+drz5/tWmpqY6xR8ja1MdYoGVsr9/2uZ/5GzNjUxlijvNdNbYw1SsbWGvG9hq7xWja1MVbNjLUuTpCfEKK1zNPsIqIPcBOwFfCv5CeH+ABwHrAG+EzRJ0vpauBqgPHjx6empqYOhV4vV3qo0/uuglwu1zLXtaXXrWn+XOkhX8d2yHVss4pnvGdmY+ToqFzpoYbJCB16v2ueP9f+TaqSsY2fyebPt8Gf7UaRKz3UMHlzpYcaJmMrLd7vRvk93kpDZsyVHmqY9zpXeqhhMrbSkO81dI3XMld6qJoZa12cFpM/6tTaIIofiWruU+QL5o4ppecLy+6PiKXA1RFxZUrp8YollSRJkqSCWl/jNJf8dU6tjQaezNh2HLC4WWla5/eFz7t2MpskSZIkFVXr4nQXsFdE7LBuQUSMIj/V+F0Z2/4DGBQRO7Za/sHC55cqlFGSJEmSWqh1cfoBMB+4MyKOjIiJwJ3A34Cr1q0UESMjYk1EnNds2xnAMuB/ImJSRBwUEWcDFwN/ID+luSRJkiRVXE2LU2HK8AnAM8D1wE+BecCElNLyZqsG0LN5vpTSfGAv4E/ARcD/kL+h7tXAh1NKa2vwLUiSJEnqhmo+q15KaQHwsYx15lNkpr2U0pPAsdVJJkmSJEnF1WM6ckmSJAFMG9jiYRO8O9Xy9iNqm0VSm2p9jZMkSZIkdTkWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnKYHGSJEmSpAwWJ0mSJEnK0KveASRJqqlpA9d/2QSQaz62tLZZJEldhkecJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMpRdnCKiX0R8ISJuiYj7ImKnwvLjI2KX6kWUJEmSpPrqVc5KEbEdkAOGA08DY4H+heGDgIOBk6uQT5IkSZLqrtwjTpcAq4CdgPcD0WzsN8ABFc4lSZIkSQ2jrCNOwIeBU1NKCyKiZ6uxl4BtKxtLkiRJkhpHuUecNgGWlRgbCLxdmTiSJEmS1HjKLU5PAB8rMfYR4A+ViSNJkiRJjafcU/W+BdwSEQD/XVg2OiKOBD4FTKxCNkmSJElqCGUVp5TSbRHxWeAbwEmFxdeRP33v8ymle6qUT5IkSZLqrtzpyAcCPwauB/YGtgIWAQ+llEpd+yRJkiRJG4XM4hQRvciXpKNSSncDv6p6KkmSJElqIJmTQ6SU1gCvAO9UP44kSZIkNZ5yZ9X7CXByNYNIkiRJUqMqd1a9+cAnIuIR4E7gZSA1XyGl9KPKRpMkSZKkxlBucbqi8Hlb4P1FxhNgcZIkSZK0USr3VL3tMz52KPcJI2K7iLglIpZGxBsRcVtEjGjH9rtGxM0R8VpErIiIv0TEGeVuL0mSJEntVe59nF6oxJNFxGbAbGAVMIn8kaqLgPsi4r0ppTczth9f2D5H/pqrpcBOwOaVyCdJkiRJxZR7qh4AETEWOBAYTH6K8vtTSn9uxy5OIX90aueU0nOFfT4BPAucBlzaxnP3AK4Ffp1SOqrZ0H3t+R4kSZIkqb3KvQFuL2AG8HEgmg2liPhvYHJKqZzpyicCD68rTQAppXkR8VvgSNooTkATMBr4dDmZJUmSJKlSyr3G6XzgWOA88tc09S18Pg84rvC5HGOAYkeo5pIvRW3Zr/C5T0Q8HBFvR8SrEfHdiOhb5vNLkiRJUruVe6reCcCFKaWvNlv2AvDViOgJTCFfrrIMBhYXWf46MChj2/cUPt8EfA+YCowHpgPbAUcV2ygiTgVOBRg2bBi5XK6MmKU1tTHW2X1Xw/Lly8vOVcv8TW2M+TqWr6mD2zXKa9woOZraGGuUjNCx97vW+Zs6sE0jvcaNkqWpjTEzdly5v8v9Pd5SUxtjjfJeN7Ux1igZW2vEn0foGq9lUxtj1cxYbnF6D/C7EmMPAee04zlTkWVRZFlr646O/SSltO4IV65Q3L4REaNTSk9u8GQpXQ1cDTB+/PjU1NTUjqhF5EoPdXrfVZDL5Vrmurb0ujXNnys95OvYDrmObVbxjPfMbIwcHZUrPdQwGaFD73fN8+fav0lVMvozWX250kMNk7GVFr/L/T1evlzpoYZ5r3OlhxomYysN+fMIXeO1zJUeqmbGck/V+zuwb4mxfQrj5VhM/qhTa4MofiSquUWFz79stXxW4fPuZWaQJEmSpHYp94jTT4FzImJt4euXga2B48kfbfpmmfuZS/46p9ZGAxscLSqyLWx4xGrd0aq1ZWaQJEmSpHYptzhNIz+N+AWFr9cJ4IbC8nLcBVwcETuklP4KEBGjyB/Nmpqx7S/I3//pMODnzZYfWvj8aJkZJElSGUZNLX3K4/w+NQwiSQ2g3BvgrgE+ERFfBQ4gf7rd68Bvil1X1IYfAJ8H7oyIc8kfPboQ+Btw1bqVImIk8DwwPaU0vZBhUUR8HfjPiHiD/I1wx5Of0e/a5lOcS5IkSVIltesGuCmlubx7yly7pZTejIgJwLeB68kfsfo1cGZKaXmzVQPoyYbXYE0HlgGfBb5E/pTBb5EvX5IkSZJUFeXeAHcKMDKlNK3I2DRgXkqpjflA3pVSWgB8LGOd+RSZaS+llMjfJLetG+VKkiRJUkWVO6veGbw7q11rrwJnViaOJEmSJDWecovTjpQ+Re8p4P9UJo4kSZIkNZ5yi9MaYEiJsaEVyiJJkiRJDanc4vR74NMlxj4NPFKZOJIkSZLUeMqdVe+rwK8i4n+BHwIvAdsCJwPvAz5cnXiSJEmSVH/l3sfpNxFxNHAZze63BMwHPpZSylU+miRJkiQ1hrLv45RSupP8jWt3BrYEXkspPVO1ZJIkSZLUINp1A1yAlNJf1n0dEVumlEpNUy5JkiRJG4WyJoeIiFMi4uxmj8dFxIvAqxHxaERsXbWEkiRJklRn5c6qdzqwotnjS4El5G98OxCYXuFckiRJktQwyj1VbwTwNEBEDAQOBP45pfQ/EbEI+HqV8kmSJElS3ZV7xKknsLbw9X5AAnKFx38DtqpsLEmSJElqHOUWp2eBwwtfHw88lFJ6q/D4PcDrlQ4mSZIkSY2i3FP1Lgauj4hJwCDgmGZjBwFPVDqYJEmSJDWKcm+A+98RsQD4IPBISun+ZsOvAHdVI5wkSZIkNYL23AD3QeDBIsvPr2giSZIkSWow5V7jJEmSJEndlsVJkiRJkjJYnCRJkiQpg8VJkiRJkjKUVZwiYkRE9C4x1isiRlQ2liRJkiQ1jnKPOM0D9igxtlthXJIkSZI2SuUWp2hjrDewtgJZJEmSJKkhlbyPU0RsAQxutmjbiNih1Wp9gUnAP6qQTV3cqKkzS47N71PDIJIkSVIntXUD3DOA84FU+LilxHpRWE+SJEmSNkptFac7gPnki9GPgIuA51utswp4MqX0RFXSSZIkSVIDKFmcUkqPA48DRH2rYWAAACAASURBVEQCfp5SWlSrYJIkSZLUKMqdHOJ3wNhiAxFxQETsVLlIkiRJktRYyi1OlwEfLTF2BPDtysSRJEmSpMZTbnEaD9xfYux+YM/KxJEkSZKkxlNuceoPrCwx9jYwsDJxJEmSJKnxlFuc/gp8qMTYBPKz70mSJEnSRqnc4nQdcFZEfC4iNgWIiE0j4nPAmcC11QooSZIkSfXW1n2cmruY/HVMlwPfiYjXgcHki9etwDerE0+SJEmS6q+s4pRSegc4OiImAB8GtgReA2allHLViydJkiRJ9VfuEScAUkqzgdlVyiJJkiRJDanca5wkSZIkqdsquzhFxKkR8VhEvBUR77T+qGZISZIkSaqnsopTRJxIfmKIR4A+wI+BnwBvAM8D06sVUJIkSZLqrdwjTmcCXwc+U3j8/ZTSJGAHYAWwqArZJEmSJKkhlDs5xE7A/cDawscmACmlxRHxVeCrwPeqklCSKmTU1JlFl8/vU+MgkiTVUam/D8G/E9tS7hGnFUCPlFIC/kH+SNM6y4H3VDqYJEmSJDWKco84zQF2BH4FPAD8R0TMA9YA04Cnq5JOkiRJkhpAucXpat49yvSf5AvUg4XHy4B/rnAuSZIkSWoYZRWnlNJNzb5+LiLGAPsAfYGHUkqvVSmfJEmSJNVduUecWkgpvQn8ssJZJEmSJKkhlV2cIqIncCKwN7At8BLwEHB9Sskb4EqSJEnaaJV7A9yRwFzgGuAwYKvC5x8Bfy6MS5IkSdJGqdzpyL8HDAD2SymNSCntmVIaAewPDAQur1ZASZIkSaq3covTBOArKaWHmi9MKf0W+I/CuCRJkiRtlMotTsuBV0uMvQq8VZk4kiRJktR4yi1OPwE+XWLsNOC6ysSRJEmSpMZT7qx6zwHHRMQc4FbgFWAYcDTQH/hFRJy0buWU0o8qHVSSJEmS6qXc4nRF4fNwYEyR8e83+zqRn21PkiRJkjYK5Ran7auaQpIkSZIaWFnFKaX0QrWDSJIkSVKjKveIUwsRscGkEimltZ2PI0mSpFobNXVmybH5fWoYRGpgZc2qFxF9I+IbEfF8RKwC3m71sbqKGSVJkiSprso94vR94JPA3cCNWJQkSZIkdSPlFqeJwJdSSt+tZhhJkiRJakTl3gB3FfBUNYNIkiRJUqMqtzjNAI6vYg5JkiRJaljlnqr3n8B/RcQs4F5gcesVUkre9FaSJEnSRqnc4vR+8tc5bQUcXGQ8ARYnSZIkSRulcovTlcAi4BTgaZxVT5IkSVI3Um5x2gU4OqX0P9UMI0mSJEmNqNzJIf4C9KtmEEmSJElqVOUWp6nAuRExspphJEmSJKkRlXuq3rnkJ4Z4JiKeYcNZ9VJK6cCKJpMkSZKkBlHuEad3yE8K8RDwWuFx84+15T5hRGwXEbdExNKIeCMibouIEe3MTUR8JSJSRDzY3m0lSZIkqT3KOuKUUmqqxJNFxGbAbGAVMIn8NOYXAfdFxHtTSm+WuZ8dgHOAVyuRS5IkSZLaUu6pepVyCrADsHNK6TmAiHgCeBY4Dbi0zP38F/BTYGdq/z1IkiRJ6mZKlo7CUZ2ypZT+WsZqE4GH15WmwnbzIuK3wJGUUZwi4hPA+4CPA7e1J6MkSZIkdURbR2ueI38qXbl6lrHOGODOIsvnAsdkbRwRg4BvA/+eUno9ItoRT5IkSZI6pq3iNKUKzzeYDWfkA3gdGFTG9t8CngFmVDCTJEmSJLWpZHFKKV1bpecsdhQr89BRROwPnAi8L6VU9pGwiDgVOBVg2LBh5HK5cjctqqmNsc7uuxqWL19edq5Gyd8oOZpr1NexqYPbNcpr3Cg52tJIGZs6sE2t8zd1YJtGeo0bJUtTG2ONkrEtjZqx3N/l/h4vX6PkaGpjrFEyttaIP49ZGiVLUxtj1cxY64kVFpM/6tTaIIofiWruKuAa4MWI2KKwrBfQs/B4RUppVeuNUkpXA1cDjB8/PjU1NXUwekGu9FCn910FuVyuZa426nDF898zs0Ob+Tq2Q65jm3Xb97oDORvq5zHX/k1qnj/X/k2qkrGr/EyWkis9VNOMXf11bKXF73J/j7fU1d/rXOmhhsnYSl1/Hn2/O6TWxWku+eucWhsNPJmx7a6Fj08XGVsMnAVc1ql0kiRJkupu1NTS5W5+nxoGaabWxeku4OKI2GHdLHwRMQrYF5iase1BRZZdRn5SitPJT2YhSZIkSRVX6+L0A+DzwJ0RcS75650uBP5G/lQ8ACJiJPA8MD2lNB0gpZRrvbOIWAL0KjYmSZIkSZXSo5ZPllJ6E5hAfma868nfxHYeMCGltLzZqkH+SFJN80mSJElSMbU+4kRKaQHwsYx15lPGTHsppabKpJIkSZKk0jyiI0mSJEkZLE6SJEmSlMHiJEmSJEkZLE6SJEmSlKHmk0NIkkprxBv+SZIkjzhJkiRJUiaPOEmS2sWjYpKk7sgjTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRmcHEJqcF6IL0mSVH8ecZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDBYnSZIkScpgcZIkSZKkDL3qHUCSJGljNmrqzJJj8/vUMIikTvGIkyRJkiRlsDhJkiRJUgaLkyRJkiRlsDhJkiRJUgaLkyRJkiRlsDhJkiRJUganI5ckqQ6colqSuhaPOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElSBouTJEmSJGWwOEmSJElShl71DiA1rGkD13/ZBJBrNrb9iNpmkSRJUl1ZnDZG/oNfkiRJqihP1ZMkSZKkDBYnSZIkScpQ8+IUEdtFxC0RsTQi3oiI2yIi8/yxiBgfEVdHxNMR8VZELIiIn0bE9rXILUmSJKn7qmlxiojNgNnALsAk4F+BnYD7IqJfxubHA2OA7wIfAaYC7wMejYjtqhZakiRJUrdX68khTgF2AHZOKT0HEBFPAM8CpwGXtrHtN1NKC5sviIjfAvMK+z2vKoklSZIkdXu1Lk4TgYfXlSaAlNK8QgE6kjaKU+vSVFj2QkQsBLatRlhJkiR1beOuHVdybM6kOTVMoq6u1tc4jQH+XGT5XGB0e3cWEbsCWwFPdTKXJEmSJJVU6+I0GFhcZPnrwKD27CgiegFXAguBazofTZIkSZKKq8cNcFORZdGB/XwP2Ac4PKVUrIzldxxxKnAqwLBhw8jlch14qnc1tTHW2X1XSlMHt2uU/I2So6mD2zVK/rY0SsZGydGWrpCxLV0hfyNlbKQspZix45YvX15WtkbN31yjZGyUHE0d3K7W+ZtyR777NcC6p9++9ATTjfIaQ2NlKaWaGWtdnBaTP+rU2iCKH4kqKiK+Tr4MTUopzWpr3ZTS1cDVAOPHj09NTU1lhy0qV3qo0/uulFzHNqt4/ntmNkaOjsp1bLNGeR3b0igZa/5edyBnV8jYlqrk34gz1vT9NmPN5XK5d7NdW3q9Rvkd2ZZGydgw73WuY5vVPH+u/Zv4O7J9qpmx1sVpLvnrnFobDTxZzg4i4hzyU5F/IaV0fQWzSZIkSVJRtb7G6S5gr4jYYd2CiBgF7FsYa1NEfAG4CDgnpXR5lTJKkiRJUgu1Lk4/AOYDd0bEkRExEbgT+Btw1bqVImJkRKyJiPOaLTseuAy4B5gdEXs1+2j3jHySJEmSVK6anqqXUnozIiYA3wauJz8pxK+BM1NKy5utGkBPWha7wwrLDyt8NPcbOn5d4AZGTS19TuX8PpV6FkmSJEldRc1n1UspLQA+lrHOfFrNtJdSmgxMrlYuSZIkSSqlHtORS5LUkMZdO67o8jmT5tQ4iSSp0dT6GidJkiRJ6nIsTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRm8j5MkSZI6pNS9z8D7n2nj4xEnSZIkScrgESdJdeP/VEqSpK7CI06SJEmSlMHiJEmSJEkZPFVPkiRtnKYNbPGwCSBXeLD9iNpmkdTlWZwqqNT1Gl6rIUmSJHVtnqonSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUoVe9A0iSJEnlGDV1ZtHl8/vUOIi6JY84SZIkSVIGi5MkSZIkZbA4SZIkSVIGi5MkSZIkZbA4SZIkSVIGi5MkSZIkZbA4SZIkSVIGi5MkSZIkZbA4SZIkSVIGi5MkSZIkZbA4SZIkSVIGi5MkSZIkZbA4SZIkSVKGXvUOoI4ZNXVmybH5fWoYpIvzdZSkyht37biSY3MmzalhEkmqHIuTJEldiKVEkurDU/UkSZIkKYNHnCRJkiRlKnXEu7sc7faIkyRJkiRl8IiTGo7n70uSJKnRWJwkqQ0WeUmSBBYnSZIkqWK81cnGy2ucJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJCmDxUmSJEmSMlicJEmSJClDr3oHkCRJ6qhRU2eWHJvfp4ZBJG30POIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRksTpIkSZKUweIkSZIkSRlqXpwiYruIuCUilkbEGxFxW0SMKHPbPhHxrYh4OSJWRMTvIuKAameWJEmS1L3V9Aa4EbEZMBtYBUwCEnARcF9EvDel9GbGLq4BDgfOBv4KfA64NyL2Tin9qXrJJUmSJDW6cdeOK7p8zqQ5nd53TYsTcAqwA7BzSuk5gIh4AngWOA24tNSGEbEb8AngpJTSjwvLfgPMBaYDE6sbXZIkSVJ3VetT9SYCD68rTQAppXnAb4Ejy9j2beCmZtuuAW4EDo2ITSsfV5IkSZJqX5zGAH8usnwuMLqMbeellN4qsu0mwI6djydJkiRJG4qUUu2eLGI1cGlKaWqr5RcBU1NKJU8djIhZwICU0l6tlh8M/BI4IKX0QJHtTgVOLTzcGfhL576LFoYAr1Vwf9VgxsroChmha+Q0Y2WYsTLMWBldISN0jZxmrAwzVkZ3zTgypTS09cJaX+ME+QkhWosytouObJtSuhq4uoz9t1tEPJpSGl+NfVeKGSujK2SErpHTjJVhxsowY2V0hYzQNXKasTLMWBlmbKnWp+otBgYXWT6oMNaW19vYdt24JEmSJFVcrYvTXPLXKrU2GniyjG23L0xp3nrb1cBzG24iSZIkSZ1X6+J0F7BXROywbkFEjAL2LYxlbdsbOKbZtr2A44BZKaVVlQ5bhqqcAlhhZqyMrpARukZOM1aGGSvDjJXRFTJC18hpxsowY2WYsZlaTw7RD3gcWAGcS/6apQuB/sB7U0rLC+uNBJ4HpqeUpjfb/kbgUPI3wJ0HfAY4AtgnpfTHmn0jkiRJkrqVmh5xSim9CUwAngGuB35KvgBNWFeaCgLoWSTfFODHwEXATGA74DBLkyRJkqRqqukRJ0mSJEnqimp9jVOXExF7R8TPIuLvEbE6IhZFxC8jYlJE9Gy23uiI+HFEvBARqyJiaUQ8EBFfiIg+ZjRjJ7NNjohU4mNJq3X3jIhbI+KVQr75EXFFRLynGtnMaMYGz7jBzdEjoldhbFrhcVPh8cEl9vViRMwwoxnNWN+MUr3V4z5OXUZEnAlcCswGvgy8QH7680OA/wKWAHdGxDHkTz38M/lrtp4F+gEHAheQP/XwO2Y0YwUcA7zYatmaZt/Hv5I/nfVB4Azg78CuwL8DR0fEh1JKf65iPjOasdEySpJUERanEiLiAPL/kP5eSukLrYbvjIhLgX4RsRNwHfAL4JiU0ppm6/1PRFwM/F8zmrFC/pRSKjr1fkTsDPwAuAM4NqW0tjB0f0TcAvwv8LOIeG+r7GY048acUdJGLCI2rdPMymUzY2U0QkZP1SttKvmb6v57scGU0vMppSeAs8gX0M8W+8s/pbQwpfRbM5qxBs4kP6nK6c3+kQpASmkR8B/k/7d/Yh2yrWPGyjCjpG4nIqYVThEcGxH3RsRy8v8Bs1lEfDMi5kX+VPt5EXFORPRotu3mEXF5RCyI/GnDr0TEryJiFzOasVwecSoi8tezNAF3pJRWZqx+MPBISunlqgdrxoyV0RUyttIz8vcva25t4R+mHwIebSPfTGAt+e/jNjOasRtn7Fl0TehRZN1aMGNlmLEyukLGO4FrgG+SPwhwLzCa/Cn0c4C9gP8EBgNfLGzzbfL/UfMf5E+z35L8fUS3MKMZy2VxKm4I0Jf8dS5ZtgP+UN04RZmxMrpCxuaeLrJsJvn7mbWZL6X0ZkQsBEZWKds6ZqwMM1ZGsYyl3Fu1FG0zY2WYsTK6QsbvppS+A+uvpdwPODCldH9h/NcRAXB+RHwzpfQqsDfw05TSNc32c7sZzdgeFiepazmKDS/GX1JsxRKC/P/yV5MZK8OMlVEsY0/g4SLrfg74fZHlMysdqhUzVoYZK6MrZGz+j+DDyP/n50Otjn7NIn/fz72Au4BHgMkR8Vph7LGU0jtmNGN7WJyKWwSsoLz/Cf1bmetVmhkroytkbO7PpS7GJ/8X3ahSG0ZEP/JH2F6qQq7mzFgZZqyMDTK2cWrRMymlR1svjIi3q5LsXWasDDNWRlfI2PwU4K3I/91c6jm3LHw+HfgHcBLwVeD1iLgOOCel9JYZzVgOi1MRKaU1EZEDPhzZM3j8Cjg5IrZOKf2jNgnNWCldIWM7/Br4VERsU+K6ksPJnx/8m9rGasGMlWFGSd1Zavb1ImAecGyJdecDpJSWA18BvhIRI4GjgW8Aq8nfhsSMZszkrHqlfYN8c/1WscGI2D4i3kv+ArR3gO9Hs5ukNltvSETsa0Yz1sB3yJ/2dHnz2WUAImIw8DXy/wNTzXORs5ixMswoSXn3kL+mcnlK6dEiH6+13iCl9EJK6RLykwuMNaMZy+URpxJSSvdHxL8Bl0bErsAMYAH5m6J+CDgZ+ERK6c6IOBH4CfBwRFzJuzdF3R84DZgOVHyaajN2n4zN7B4RQ4osfzSl9FREnAb8kPzFkleSP8S9C/mp1ocDh5Qxe6AZzbgxZZS0cfspMIX875lLgMeBTYD/Q35mtX9OKb0VEb8jf/3LHGA5+RvX7wZca0Yzlsvi1IaU0mUR8Xvy9++5mPz5+MuAR8n/I/nuwno3R8STwNnA+cDW5K+beYL8FInXbLh3M5qxQ24usXwo8FpKaUZEPE3+UPT3yB9JC/KHtvdIKT1V5XxmNGOjZZS0EUspvR0Rh5K/J+OpwPbAm8Dz5CepWF1Y9X7yp3hNJf/v378CZ6WUvmtGM5YrUkrZa0nqsiLia+T/h//olNId9c5TjBkrw4ySJFWPxUnayEVEkD+8fRRwWEqp4S7GN2NlmFGSpOqxOEmSJElSBmfVkyRJkqQMFidJkiRJymBxaiUipkVEitJ3ySYimgrrNFUxx0kR8WxErI6IJdV6niLPOyMi5jd7PKrwvU4uc/vJhfV3LDJW8rVt9jwnt1oeEfHJiPh1RCyKiLcj4sWIuDEiDmrH9/XPhSnH260z2xa2X/eajOroPjr4vPMjYkYZ67V4zzv4XJ16jWqh2ftQ7ONfC5/PK7LdU4Wx0a2Wjyks/3ztvov66QrvcSkRsXdE/Cwi/l74nbooIn4ZEZMiomez3+mlPrao9/eQpb2/q0vsY92fkSURMajVWK/C2LRWy0dHxI8j4oWIWBURSyPigYj4QkT0aZWt1MfuHcj6/9s783i7prOPf38ZiKHGUjFUqpSUt9UWbxSRKm0+lFJT3taQGt8qamhLhEoJQREdDK20DaVoUXNJSSTSCPGaqzEmxhgiMaQZ8bx/PGvn7uy7z7nn3HvOvcHz/Xz255y9pr322ms/a3rWs78h6e/pWc6X9JSks0vyPT13nfdS+PsknVUmk3PtVHa8Jel+Sd+tM39FefOupEckHZlvAyUtJ2lI8pubym9CtevVUuaB01lyS13UZ6uFrE538jU7JI8K78/nSvzzMnvH5Ja9u/MkrdxGmq36qLUS5sjbx4PA1sATzUhc0trA72ixV9+Z3zk5Hf9wZZcj/8js1fgm8suAXwOz8I+f7Y3b8V/VzN6uIbndgR2B89uRlY7E7Ur2AN7ppGt9mMpob+ClgtsTwDNA/7yjpDXw7w7NTX75dz4LO6E52Vzq+DA948VIOgbP81jcLPrz+DfavgFcDLwFZDLkaGBKSTLvNj+nSxUr42V1YrVAkvYG/gQ8jrcd2Xfttgd+jpuez7cnI/BvqxR5qp7MSToJOAO4Af/O3izgKynPe0r6mpm9mItyBzAs5WcV4MvAocCRkvY3s7KPMG+Lf/B8tRT2Skm9zOwP9eSVFnmzUvr/a2BN4Gepc/cPoC9wHi5LegHfAa6QtL2ZHV6493rL/ONO0+VWF/fZPuq8C+yPf+4lzwHJ7xMlcRYBe9H68zDV4tRMDJzagZm9A0xu4iU2AroDl5nZxCZepxVm9mxnXq8NhuCVfy8zu67gd6Wkb+AvSFCCmT3U1XlYSnnYzJ4pOkqaAAyS1NPMsnrVH+9U/y39vyQXpT/eYXusyfltKGlCQmb2XhfnY1kzW1DiLqCnmS0siVbvNfrjHabfmNnRBe8bJZ2PdzpXS27/NrNmyvaKVCqPLmIMcJSkC8zs1bIAkjYCLgf+DuxdqE+3SToXKM4UP9fR8pVrGgwHLjCzY3Ne4yX9Dfi/lK+8RsLMwnX/LumX+IDqSkmfM7PiZMp92T1JGgP8GzgGqHfglJc3Y9JM9zHAz/ABzheBbc0sP2C/TdJjwAWSJpnZZSkf7SnzmlnK6uBSS0k5dVmf7WPA9cB+kn5myZqdpOWAPYHrgMEV4uxPbuAkaT18cuHyCnFqJlT1KtNX0ri0dD5D0mmSukG5qp6cYyU9mZZqZ0j6jaSV6rmoXLXq7nR6V7rOaEk9JQ2Xqx0sTL/DJfXMxe0h6XRJz8pVF2ZKmihp21yY70p6SNKctLz/mKTD89dXB9W2GoGkZYDjgVtLBk0AmNkYM5tbQ1qjgQOBdXLLtNOT38aS/iZXx5gnabKkgTXG7SVppKTHU3m+KulmSZt04L43lPQnSdNSfp6TdLEK6icp7PZydaO3Jf1HrupxcM6/laqepK9LejDVj2fzz74DeR5N+8v3kBR+95xbd7mqyrOSOjQzVCfjgeXxmeuM/sCk5Ne/EH47YGImzNuDWlQLNpF0R3qOL0j6fvLfX9LUVL/GSfpsLu4gSWMlvZH8H5J0YMk1TNIZkk6UNA3/QOB/1ZnP0VR+xp9MdfRludrQVEmHFeJnKhL9Jf1VrspyX/KbLukKuarL1JS/XerJXxVOxAe3Py3zNLNnzezRBl2rZnLPfbP03OcAf5Grn90mbz/mJtlyvHywm4+/vKSL5GpncyTdBKzbwCwOT79Dq4Q5Fp98PaJsEG5mb5jZPxuYp4yf4s90SMk1pwFnAQMk/Xe1RMxsDnAEsBz+kfNqYd8DHgLard6TYwrwCUmfBvYDRhUGTRm/wle486t+DSvzKnVwebnK4zR5X2OapKFK/Z8Ud0VJv06yaoGk1yTdmW/7cnJnqFy9fl6S663UMiVtkcLn+ypHJbfhObeNktvOktaQ9Fu5iuZcSS9K+rOkdXLhR1NBbiX/Dsmu3DXuTqed3mdrD3J10XslzZK3zZMl7VIIk6naHS7v/85IYW+WtG4hbDPl0Z+A9fEV4Iw98IFqad8QHxz1l7R+zm1/4AUaoCESK06VuQGfWRoBfBNfJvwAX+4v4wxckF8I3Ax8Hl9G/6J8uf2DGq97Oj5j9ivgh7ha4Bu4qto+wJnARFxV8GRgAyDThT4BF6xDgYdx1YAtSLOp6WW8IqX9E3zgvAmuurC0sQWerzK1jno5HVgD2BLYLbktkC+vT8SXbo/EVxZ+CNwq6Vtm9vdKcdPvsviS73BgBl7ORwCTJW1Saaa2DdbG1TqOAWbjz/ck4Db8mQMg6du40Pgn3ujPBDbFBUwpkvqmdB4ABqX8DwNWxFVS2ku7y9fMRslXDkdJmmJmL+Pv2tb4LGwzVKS6a8l9dmZm79MiUPvTsqLcH/gLcA/e+G5gZs/JBzDrACMblKe/ApcC5+J16A/y2eUBeMepJz47/Wcg6xBuAFyLdxQ/SHkdJWk5M7uEJRmMfzn9x/hX1l+pM3+VnvFKeB1cDq9L03B5ebF8VvbXhXSuBK7CV5Lzz+BrwOa4qtHrwPQ689eKNNgYANxgZrWqznRT6z2YWf1oBjfis6Jn489wE+AuXJ1rPi4Hh+Fln+9A/xbYFy+vKcBOeN1oFDOA3wDHSDrXzJ4vCbMjMMXMZtSRbofKN8XdHrixyjO9CS/PHch1cMsws0ckvQJsU8PlP4OrdXaUz+Dydlu881faxpmZSboZOEFS71TO7SnztsjXwW74KlzWf3kM6IfL5NXwyUxwubcb3jY9DayOl2GxL3EA3lk9Em9vTsMHFxuZ2axcuAfxst0BbzNI/+elX3Ju7+PyeG38HRmC95HWTvn7Z2p/51Ol/W6g7Pow9tn6AKNwOdsD2BW4RdLOqd+TZwg+eXgQrmJ6Hl4Wji+I0AAAFTtJREFU2+fCNFMePY+3zfvjzx28Xv0NmFMhzj34vX0PL39S/CuAju/1MrM4cgf+AhlwYsH9UrwDuAreGBswIPmthr/Aowtx9kvhdqszDzsW0t8snQ8rhDs5uX8hnd8CXF8l3R8Ds9q49mhgeu68T7rG4BrzPjiF37BK2fYo8cuuc0g63zedf7NBz3U08FLB7VzgvXxe8YbsSeDBanErXKM7vmLxLnBsSZn0aUe+e+ANrAFfSm7ChcIDQLcqcafn6yQu7GYCK+Tc1sNn+KfXm7cGlu8quHAchwvj94AhjXjuFepm8ZiYC/M8cEv6vzKpg5POX8neA1yP3YAtOpin7J04IOe2aiqDN4GVcu5Hp7Drl6TTLdWVS4FHCn6W8r5cE57xKbjs26jgfmmqaz0KZT+yQj2dC6zV4Of9qXTNETWEHVChbhjweBPqYvbcf1QljNIzHYpPonRL7hunellsoy6mDlld4ZrZc9oQb9feAv6Q/HqQa4fwTu1VNabbp0r5zmnkM8X3CBlwUa5+XVEl/L24imbx2Syb7nlN4NTkdkE7ynLjlM6q+CTX+/jE7AmZf5U0/jeF2areMm9PHcQ7lwb0L4QdircTa6bzx4Hz20jfaN3e9MHV608vCX8jMC7974avKp6Xwq+Y3K8GJle4Xne8PTNgj5z7aErabxogu3JxuqzPVs+zruCXtR1j8AmJ4js7viRPBqydzpsij1hSFh2Ey8BeQG+8fdyJFrm9Y6FO98AH6f9O7lsl942o0ket9QhVvcr8pXB+NT4zv1lJ2H64kL2iJM57LDkybw+ZilAx/ew8S38KsHNaHt9Wru6WZwqwqlwt5lv6EFiKajL9cSG8eL+L+cznVcDmqkHNUtI+cgtNb+HP+j94Pdm4PRmStIykk5LKwDy80chmWTbO/a6Pq3jUupIJPuN1m5n9J3Mw30DdDHUaqLF8zewtfAZuO3y28x589rNZ7IHPQGbHwTm/CcA2crWUbfHyz9RoJtLyLvbHB8iN2ke2eJbPzGbjqy6TzfdTZkxNv+vBYrWVqyS9nPK5CN8oX1b3bjezeQ3Ka56B+Kz+tKR20iOtCtyBz0J/vhC+bBM++L22Z4W20fyQJevGlvgkTrNYojwk9U7qR8/jHdVF+Ir2KngHHnzFsRvlbVTDMF8ROA84QFK75FkJw2ldvtvVEV8NykcxTStxn4+X/2v4ysoFtGEsowJTUzqzgIvwCayDqO1emnG/RfJ1cCA+eTSp8D6PwVe9+6VwU4DBqa3aQgVV0hzF9mY6vpq/dUnYccDWcquAm+N1/hx8hSirIwNwIy8ASPqBXEV9Dt7+vpC8aqmvjZJdZSzVfTZJX5F0i6TX8HJbhA9Eysrt1sJ5tqf30+m3M+TRX/E+9q74KtKr+Mp8NS4HNpG0Jb5CNdnMnm5EZkJVrzKvVThfp8Qv21i8xPK5mb0n6c2cf3spTR+vPHn/M3Fhvx8u6OdIuhb4iZnNNLPxcos8R5GEgKTxwHHWOXr+mU5299x/cm75MJlFpPWbmJ/VKO/4voo3WKtSxSqdpF2Ba/Al+Z/js1Qf4Opw7TUJOwJ/Pqfhy+Pv4rrC1+fSXD39Fjczt0VvWtddkttn6s5p29RTvpPxlajPA7+sc0BYL49biXGIxAT8/fkC3vjdby2bgO8BfpT+9wcmWeNUuGYXzhdWcAPoJWlF3BrXXLwz92zy/wHeKSvSSNWePGviM4KVjLSsXjivlI9m5O9NfIa+HhnylJk90IS8VGLxfafB+k24ytEwvMM9D7cKNpSW9793+q3URjWSkbTIo+8V/F6kfvn8fAfLdyZeJn2qhMn8XqwSJs96tExK5OmHz6TPBl6wFoMx9bIHLqvfxe9/PoCkLH99cNlXRla+maxvT5m3Rf7dWzOl39b7fBQuxw/CtynMknQ5MNSW3Hdcqb3ZtMR9LN45/irwJXzl/DVJE4GvSXoBX3EcB74HCldhOx9XY5uNd+AnU1v72yjZVcZS22eTG0m4C98/dxQ+2HwPVznsWxJlVuE8aw87TR6Z2buSbsBXRPsAV5rZB1LleQUze0bSvfjE6F60tsrXbmLgVJlP4XsC8ucAL9O63LKKtRbwr8wxzV6sjjfgHSGfft7q3Vrp902AJNjPBs6WtBbwLVyoLE+aNTWza4FrU8drQAp/u6R1m9xZBZ9FB+8YTCv4rZ1+s5ftAVxNZFfczGczmEVLGeZZC5+BLAqMIoOAZ8xscOYg3/jZkYHyIOByM8tviF2xEGZm+l2H+phBSz3OU+bWCOop31PxZfRHgZGSxlltZuYbTX6fU39ys5v4wOlXaQZrA1xHvKvYGu/gbGc5K04l+0cyymbUG8Gb+Hv9owr+xQ5hpXw0PH9p4upuYCctvdbC8vf9WXx/w/5mtnimOk3Q5Mk6Y5XaqMZlzmyOpBH4ytMvCt53AodIWquzVgvTM52AP9NeVr7PKdvLMrbEbwnkhgrWpvxd/j9rjOXJShM1d+MTbbvhqxzFvAlv/6aaWbYnsRllnq+Db+Jt8z4Vwk6HxYY1hgBD0gb8vfC9lgtxFcSMSu3NyyXuj+Ft2w74wCl7fmNTfl5M6WcaEoOAu8ws23eFpHomABslu8pYmvtsA3E19H0sZ0lS0vLtSAs6Tx5djq9+dQP+p444F+IDw2salZFQ1atMUXAMwjeiPV4SdjI+Ch9UcN8XH2SN72BesvjF9LMZwFZWQszsVTMbhQvaVuqFZjbHzG7BN/X1pvXsSjO4O/3uWeK3Jz7zMjnlbyHeWH9LUll4JO1Ux8u+AN8Emmc80E+5jyAmlYN9gYesxTBBWVxw4VZsWPenZfWsPSxP6xmw7xfOn8IbsENUbcqlNffiagErZA5p9qmWjdFt0e7ylbQdPts2FO8orILrR3c6ZvYkPis4ELeud0/O+1F8hezkdN6V32/K6v3iuiK3vPjtJl6z7Bnfjm9WfsHMHig5uvr7R2fhsq3Y6Qe8oyXpC52bpYqUPdOetF7puQ/vcJe1Uc3gIryjO7zgPhJfkbmoTFVLbrGsEbKlyC/wZ3pm0SN1nE8AJphZVcMQqSN6Ib5q+9sm5LMq5oZw/ozL8S1LghyNr8Cfk3Nrdpnfjq/AzanwPs8sRjCz583sPHzgU+xrFNubPvhK3r0l6RjeZuyEq+blB05fwlfu7sutaNXSVkLl9ruZsmtp7rOVyZnP0f5+QGfJo3/g6oCXmNm/2gqcuAZfxT/LljRG0iFixakyhybViSm4pZVD8I1+bxX7qmY2S/49kCGS/oOravXFG5qJtNYRrQsz+5ekq4BhaUZ5Ej7jfAq+UfRRAEk3Ao/gVl1m48JmIKlRkHQaLUvdr+AqYEfj35l4oyN5LGGgpOKM2Nv4zN4ISWumfPQAdsY3zf487e/IGIF/4+IaucnPm/GZnHXxgdZ3cHWvWngCWE3SD/DVrPl4IzQY+IekU/FO8RH4dzB2qRbXzB7DBe/ukkbimzy/gpdnRywv3Q4cKP+GxzPpHr+aD2BmJv+o5/XAWEmX4FZ8+uKbd0+tkPZw/AOMYyT9AlgGVzFsxJJ6u8o3dfavxOvCueneDsPN4t5h6fslncw9+Ayq4e8aAEk1YBI+uJtH+YdSO4tJeHlemMp2BXxANxOfTWwGlZ7xvsA96T14MuVlE3w1rJkDuTYxswmSjgPOl1uVHI2rpqwKfB2X69+l5QO4fdN+iSKP5fdqNIl/4/tLzpD0Pt6xObYYyMyelPRnIPtERmbFaudmZMrMFqS243cF96clHYDv25ic5FD2MdbtcJl+GkvuodxAUj9a81StHRszu0vSz/D774PPKs/GP2p7Iv4s9y9E+2S6rvD3I/sA7hrA/+RWdDqbo/DB0Vj5N5iyD+DuiavB/d7M/pgFbmeZ10P2Ade7JJ2H9yeWwVdDdwN2N7O5SQXqJnywNAffs/NFXG09zzxa2ptl8fbmHSpbIx2LD2Yzy3ng/Zl3cMubp+XC3o5bHDwJuB9fqdqrJM1K7XfTZNdS3me7E5/wvTw94974c3mBdiymdJY8Smrxta40ZXFm4wPuxmIdtNbxUTtoscqxGV5Z5+Ez0KfTYtVoADkLKslNeCP3JL6cPAMXACu1Iw87lqTfE+/4Po83qM+n8565MMfjKzaZbv+T6X56Jv9dcJWAGfgszIu4GdK1c2mMpjFW9cqOx/EX83hc4M7HZ/seAA6ukJ5w/d9xuGBZhOt7X4ULt1rLdIUUZ3bKy/TkvjFu4ehtWla8BtYYt1t6Bq+k+xiPC77pLGnNLiuTPjXk85P4psrZ6bgS30Dd6hngDcU4vOGagwvg7+f8l8hHrm49lJ7/c3hDu8Qzb+d7067yxTd9vgH0LqQ3Ct8T0G7LN1XqZtU0cQMBRs7yX87vpOQ3rkF5GkaJpUlKLIHR2oLQDulZzsPVQY6mxHpSijO8AXmt9IxXxTsh03DZ9zre6TmmlrIvu9dGH/jkw19x2Zdt1B+Dy5ZuVLeqZ3TQemIdz31zfLJtLi7nTsMHd0vID3zW+OJ0H3PwTuw2NNCqXsG9B77SbbS2FLYpLQPShfi7fg8+SbJsCtOnjfLdqx15HYi3Z7NxefY0vhq1Wkn9yq6T7Vmagq9Grl/rs2lUWZaEWx6XK4/h7/K7qQ7sVyVOm2XegTrYi5Y9dgtSHZuS3DJLc2fjsudt3CDSY8DRhXQM3/90UqrL81MeN6+Sp74p3uSC+4207hMtl96BN1KZ3YLv1V2ijlJBbiW/DsmuXJgu67PV86xz5/uk5zsf314yiMp9v0MKaQ0oudeGy6Mayz3LSyureh19L6sdSgkFQRAEQRAEQYeRZMAZZnZym4GD4ENE7HEKgiAIgiAIgiBogxg4BUEQBEEQBEEQtEGo6gVBEARBEARBELRBrDgFQRAEQRAEQRC0QQycgiAIgiAIgiAI2iAGTkEQBEEQBEEQBG0QA6cgCIKPKZIGSzJJG3Z1XtqLpM0lDZO02sfhukEQBEHXEQOnIAiC4MPM5sCpQGcPYLrqukEQBEEXEQOnIAiCIPgIIqmnJHV1PoIgCD4qxMApCIIgWIykuyVNlDRQ0sOS5kl6SNJ/S+oh6UxJMyTNkjRa0gq5uH2S6t8Rks6X9LqkuZJukdSncJ1BksZKekPSnHSNA0vy00PSCZKekDQ/hb9d0iaSBgN/TEGfTte24rVK0jxU0oPp3mZLGi/pqzn/nyf/tyXNTPnsl/Ovet2U5yGSpkpaIOkVSedJ6lXIxwaSbktl9HoKc1jxHtIAaLik6ZIWpt/hknpWKPtzJL0CLAC+nNy/XVIOoyW9JKl7tfIKgiAInB5dnYEgCIJgqWND4BfAGcAc4BzgpnT0AAYDfVOY14GfFuIPAR4Gvg+sCZwJjJG0qZktSmE2AK4FzgI+APoDoyQtZ2aX5NK6GtgduAC4E+iVwvYGbgWGAycDewMvpTgzKt2YpHOB44Hf46p2HwD9gE8Dk1KwdYCRKb0VgP2ACZK2MLNHa7juFcCuwNkpzb7A6UAfYM+Uj2WAf6T7OSKV4yHAXiXZvgzYJ5XjRGDrdO0NgO8Wwg4FpgCHAd2BJ9L54cCNuXJYJaV5jpm9X6m8giAIghxmFkccccQRx8fwwAdABmyYc7sbWARskHPbLYW7sxD/emBa7rxPCvcE0C3nvk1yP7hCPrrhA7JLgUdy7jukeEfXcw9Vwm4IvA+cX0cZdU95exL4ZVvXBbZL7gcU3L+X3DdP54el861yYQQ8ktz7JLfN0vmwQnonJ/cvFMr+QdLH7Qt5fR9YP+d2NPAesG5X18M44ogjjg/LEap6QRAEQZGnzOy53PnU9HtHIdxUYN2SfTTXmtkH2YmZ/RNfldk6c5O0kaSrJL2MD9QW4SsuG+fS+QY+GLi0IzeTY0d8kPa7aoEk7ShpnKQ38cHFIuBzhbxVYiCwELguqez1kNQDGJP8+6fffsALZnZ/FtHMDLiukF4W/oqCe3a+fcH9hpROnquBt4BDc26HA7ea2UsEQRAENREDpyAIgqDI7ML5wiruPfBVmTyvlaT5Gq4Ch6QVcTW1LwIn4qs0WwJ/AJbNxVkdmGVm8+rMfyVWT78VBwuSvgzchqsoHowPcLbEV4J6VYqXY01gmRR/Ue54vZCH3jm3PMWyy6z2FdUPXy34UyEcZjYf35N1cBrIbQd8HrikGDYIgiCoTOxxCoIgCBrNpyq4PZz+bw2sD2xnZhOzAGllJs9MYLW076kRg6eZ6XcdXPWujD3xVabvWMt+LCStiq/atMWbwHx8MFjGK+l3Bj54KVIsu1npdy3g2Zz7Wrnr5SmuNmVcDBwHfBvYA5hO6xXEIAiCoAqx4hQEQRA0mr0kLW5fJG0DrAvcm5yWT7/FgUnR8tsYfN/PIVWutSD9LldDvu7EjUEcViXM8vh+oMUDEEk74MYjarnu7fjK1Mpm9kDJkQ2cJgOflrRV7joiGY/IMT79Diq4fy/9TqhyL4sxs2fx8vwJboDi0rw6ZRAEQdA2seIUBEEQNJpPADdI+i2wBjACeBq4PPlPAt4BLpR0Km657mR8RWjlLBEzGyfpOuB8SesBY4Ge+L6fW83sbtwQBcAPJV2GD8YeNbNMvZBces9KGgkcJ+kTuJXA94GtgKlmdg0+8DkGGC3pj/jeplOAlwvJVbru3ZKuAq6VdD5wPz5Y6wPsDJxgZk8Bo4ETgOslDQXewAeIq6Z0P0h5/ldKb1hakZuEr9idAlxlbuWvVi7CLestwtUigyAIgjqIFacgCIKg0YwAnsEHBxfhlt6+mam+mdkbuLpYd9wk+QhgFK0NIICvtAzDTZLfhHf4NyXt5TGzR5L/rrip7inA2pUyZmY/xs1/98MNMVwJfA14IfnfgVuc2wa4BTgIOCDdTz6datfdL/nthQ9UrgWOxAePr6X4C3HjF4/ie40uA14ELkxpvJ273IG4afOD8P1XB6fzVt+9aoNbgbnAjWb2aluBgyAIgiVRa+M7QRAEQVA/6aOt04BDzWxU1+bmw4mkW4C+ZvbZJqS9E66ut6OZ3dXo9IMgCD7qhKpeEARBEHQBko7Dre89jas37g3sAvygwdf5LP6x3JHAgzFoCoIgaB8xcAqCIAiCrmEBcCxueKI7bunvEDP7fYOvcwquPvgIrnYYBEEQtINQ1QuCIAiCIAiCIGiDMA4RBEEQBEEQBEHQBjFwCoIgCIIgCIIgaIMYOAVBEARBEARBELRBDJyCIAiCIAiCIAjaIAZOQRAEQRAEQRAEbRADpyAIgiAIgiAIgjb4f7fQTHMuGv3LAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 1008x576 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# normalize\n",
    "df = (df.T / df.abs().max(axis=1)).T\n",
    "\n",
    "# activate grid and set font size\n",
    "plt.rcParams.update({\n",
    "    'font.size': 16,\n",
    "    'axes.grid': True,\n",
    "    'axes.grid.which': 'major',\n",
    "    'axes.axisbelow': True,\n",
    "})\n",
    "\n",
    "# plot data\n",
    "df.plot.bar(xlabel='Impact category', ylabel='Impact score', figsize=(14,8))\n",
    "\n",
    "# rotate xticks for horizontal impact category names\n",
    "plt.xticks(rotation=0)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "There you go! A nice plot to compare the three technologies."
   ]
  }
 ],
 "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.8.6"
  },
  "latex_envs": {
   "LaTeX_envs_menu_present": true,
   "autoclose": false,
   "autocomplete": true,
   "bibliofile": "biblio.bib",
   "cite_by": "apalike",
   "current_citInitial": 1,
   "eqLabelWithNumbers": true,
   "eqNumInitial": 1,
   "hotkeys": {
    "equation": "Ctrl-E",
    "itemize": "Ctrl-I"
   },
   "labels_anchors": false,
   "latex_user_defs": false,
   "report_style_numbering": false,
   "user_envs_cfg": false
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}