Radiative forcing of emission profiles.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np\n",
    "\n",
    "def decay(t0,\n",
    "          t,\n",
    "          alpha,\n",
    "          tau):\n",
    "    if type(alpha) == float:\n",
    "        return np.exp(-(t-t0)/tau)\n",
    "    elif len(alpha) > 1: # probably CO2\n",
    "        assert len(alpha) == len(tau) + 1\n",
    "        return alpha[0] + sum(aa * np.exp(-(t-t0)/tt) for aa,tt in zip(alpha[1:],tau))\n",
    "\n",
    "def concentration(emission_profile, # list of coordinates (t,y), t is strictly increasing\n",
    "                  tau,\n",
    "                  alpha=1.0,\n",
    "                  step=1,\n",
    "                  c_0=0.0,\n",
    "                  t_max=2100):\n",
    "    \n",
    "    '''\n",
    "    \n",
    "    \n",
    "    \n",
    "    '''\n",
    "    \n",
    "    if type(alpha) == float:\n",
    "        CO2 = False\n",
    "    else:\n",
    "        CO2 = True\n",
    "    \n",
    "    n_points = len(emission_profile)\n",
    "    assert n_points >= 2\n",
    "    \n",
    "    # initialize\n",
    "    T0 = emission_profile[0][0]\n",
    "    c = np.zeros(shape=t_max-T0+1)\n",
    "    \n",
    "    for i in range(n_points-1):\n",
    "        \n",
    "        # extract the variables that we need\n",
    "        t0, y0 = emission_profile[i]\n",
    "        t1, y1 = emission_profile[i+1]\n",
    "        \n",
    "        # create the time series for the interval of interest\n",
    "        time = np.arange(t0,t1+1)\n",
    "        \n",
    "        # calculate the linear coefficients\n",
    "        a = (y1-y0)/(t1-t0)\n",
    "        b = (y0*t1-y1*t0)/(t1-t0)\n",
    "        \n",
    "        print(t0,t1,y0,y1,a,b)\n",
    "        \n",
    "        \n",
    "        if CO2:\n",
    "            # calculate the decay function\n",
    "            d = decay(t0,time,alpha,tau) # this should be fine (integrated until t_max)\n",
    "            \n",
    "            c_i = list(alpha[0]*(a*(time**2-t0**2)/2+b*(time-t0)) + \\\n",
    "                sum(aa*tt*(a*((tt-t0)*np.exp(-(time-t0)/tt)+time-tt) + b*(1-np.exp(-(time-t0)/tt))) for aa,tt in zip(alpha[1:],tau)))\n",
    "            \n",
    "        else:\n",
    "            # calculate the decay function\n",
    "            d = np.exp(-(time-t0)/tau)\n",
    "        \n",
    "            # calculate the dynamic concentrations\n",
    "            c_i = list(c_0*d + # any residue from a previous period \\\n",
    "                tau*(a*((tau-t0)*d+time-tau) + b*(1-d))) # plus the ongoing emissions\n",
    "        \n",
    "#    if t1 < t_max:\n",
    "        c_0 = c_i.pop()\n",
    "        time_to_end = np.arange(t1,t_max+1)\n",
    "        c_i.extend(c_0*decay(t1,time_to_end,alpha,tau))\n",
    "        # finally pop the last value, that's the start of the next period\n",
    "        \n",
    "        # store it\n",
    "        c += np.array([0.] * (t0-T0) + c_i)\n",
    "    \n",
    "    # if we want to calculate the concentration beyond the last t, then let's do that\n",
    "        print(c.shape)\n",
    "    \n",
    "    emissions = list(zip(*emission_profile))\n",
    "    \n",
    "    time = np.arange(T0,t_max+1)\n",
    "    concentration_series = pd.Series(c, index=time)\n",
    "    \n",
    "    emission_series = pd.Series(emissions[1],\n",
    "                                index=emissions[0]).reindex(time).interpolate(limit_area='inside')\n",
    "    \n",
    "    return emission_series, concentration_series\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "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>Chemical Formula</th>\n",
       "      <th>Lifetime (Years)</th>\n",
       "      <th>Radiative Efficiency (W m-2 ppb-1)</th>\n",
       "      <th>AGWP 20-year (W m-2 yr kg-1)</th>\n",
       "      <th>GWP 20-year</th>\n",
       "      <th>AGWP 100-year (W m-2 yr kg-1)</th>\n",
       "      <th>GWP 100-year</th>\n",
       "      <th>AGTP 20-year (K kg-1)</th>\n",
       "      <th>GTP 20-year</th>\n",
       "      <th>AGTP 50-year (K kg-1)</th>\n",
       "      <th>GTP 50-year</th>\n",
       "      <th>AGTP 100-year (K kg-1)</th>\n",
       "      <th>GTP 100-year</th>\n",
       "      <th>Molar mass (g mol-1)</th>\n",
       "      <th>Correction</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Acronym, Common Name or Chemical Name</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>(E)-1-Chloro-3,3,3-trifluoroprop-1-ene</th>\n",
       "      <td>trans-CF3CH=CHCl</td>\n",
       "      <td>0.07</td>\n",
       "      <td>0.04</td>\n",
       "      <td>1.370000e-13</td>\n",
       "      <td>5</td>\n",
       "      <td>1.370000e-13</td>\n",
       "      <td>1</td>\n",
       "      <td>1.090000e-15</td>\n",
       "      <td>2</td>\n",
       "      <td>1.540000e-16</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>1.120000e-16</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>130.49</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>(E)-HFC-1225ye</th>\n",
       "      <td>CF3CF=CHF(E)</td>\n",
       "      <td>0.01</td>\n",
       "      <td>0.01</td>\n",
       "      <td>7.250000e-15</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>7.250000e-15</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>5.770000e-17</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>8.140000e-18</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>5.950000e-18</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>132.03</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>(E)-HFC-1234ze</th>\n",
       "      <td>trans-CF3CH=CHF</td>\n",
       "      <td>0.04</td>\n",
       "      <td>0.04</td>\n",
       "      <td>8.740000e-14</td>\n",
       "      <td>4</td>\n",
       "      <td>8.740000e-14</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>6.980000e-16</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>9.820000e-17</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>7.180000e-17</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>114.04</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>(Trifluoromethyl) sulphur pentafluoride</th>\n",
       "      <td>SF5CF3</td>\n",
       "      <td>800.00</td>\n",
       "      <td>0.59</td>\n",
       "      <td>3.360000e-10</td>\n",
       "      <td>13500</td>\n",
       "      <td>1.600000e-09</td>\n",
       "      <td>17400</td>\n",
       "      <td>9.930000e-12</td>\n",
       "      <td>14500</td>\n",
       "      <td>1.100000e-11</td>\n",
       "      <td>17800</td>\n",
       "      <td>1.110000e-11</td>\n",
       "      <td>20200</td>\n",
       "      <td>196.06</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>(Z)-HFC-1225ye</th>\n",
       "      <td>CF3CF=CHF(Z)</td>\n",
       "      <td>0.02</td>\n",
       "      <td>0.02</td>\n",
       "      <td>2.140000e-14</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>2.140000e-14</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>1.700000e-16</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>2.400000e-17</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>1.760000e-17</td>\n",
       "      <td>&lt;1</td>\n",
       "      <td>132.03</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Sevoflurane (HFE-347mmz1)</th>\n",
       "      <td>(CF3)2CHOCH2F</td>\n",
       "      <td>2.20</td>\n",
       "      <td>0.32</td>\n",
       "      <td>1.980000e-11</td>\n",
       "      <td>795</td>\n",
       "      <td>1.980000e-11</td>\n",
       "      <td>216</td>\n",
       "      <td>2.060000e-13</td>\n",
       "      <td>302</td>\n",
       "      <td>8.240000e-13</td>\n",
       "      <td>38</td>\n",
       "      <td>1.640000e-14</td>\n",
       "      <td>30</td>\n",
       "      <td>200.05</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Sulphur hexafluoride</th>\n",
       "      <td>SF6</td>\n",
       "      <td>3200.00</td>\n",
       "      <td>0.57</td>\n",
       "      <td>4.370000e-10</td>\n",
       "      <td>17500</td>\n",
       "      <td>2.160000e-09</td>\n",
       "      <td>23500</td>\n",
       "      <td>1.290000e-11</td>\n",
       "      <td>18900</td>\n",
       "      <td>1.470000e-11</td>\n",
       "      <td>23800</td>\n",
       "      <td>1.540000e-11</td>\n",
       "      <td>28200</td>\n",
       "      <td>146.06</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Sulphuryl fluoride</th>\n",
       "      <td>SO2F2</td>\n",
       "      <td>36.00</td>\n",
       "      <td>0.20</td>\n",
       "      <td>1.710000e-10</td>\n",
       "      <td>6840</td>\n",
       "      <td>3.760000e-10</td>\n",
       "      <td>4090</td>\n",
       "      <td>4.580000e-12</td>\n",
       "      <td>6690</td>\n",
       "      <td>2.550000e-12</td>\n",
       "      <td>4140</td>\n",
       "      <td>9.010000e-13</td>\n",
       "      <td>1650</td>\n",
       "      <td>102.06</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Trifluoro(fluoromethoxy)methane</th>\n",
       "      <td>CH2FOCF3</td>\n",
       "      <td>4.40</td>\n",
       "      <td>0.33</td>\n",
       "      <td>6.820000e-11</td>\n",
       "      <td>2730</td>\n",
       "      <td>6.890000e-11</td>\n",
       "      <td>751</td>\n",
       "      <td>9.590000e-13</td>\n",
       "      <td>1400</td>\n",
       "      <td>9.270000e-14</td>\n",
       "      <td>150</td>\n",
       "      <td>5.720000e-14</td>\n",
       "      <td>105</td>\n",
       "      <td>118.03</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>Trifluoromethyl formate</th>\n",
       "      <td>HCOOCF3</td>\n",
       "      <td>3.50</td>\n",
       "      <td>0.31</td>\n",
       "      <td>5.370000e-11</td>\n",
       "      <td>2150</td>\n",
       "      <td>5.390000e-11</td>\n",
       "      <td>588</td>\n",
       "      <td>6.730000e-13</td>\n",
       "      <td>984</td>\n",
       "      <td>6.850000e-14</td>\n",
       "      <td>111</td>\n",
       "      <td>4.470000e-14</td>\n",
       "      <td>82</td>\n",
       "      <td>114.02</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>203 rows × 15 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                                         Chemical Formula  Lifetime (Years)  \\\n",
       "Acronym, Common Name or Chemical Name                                         \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene   trans-CF3CH=CHCl              0.07   \n",
       "(E)-HFC-1225ye                               CF3CF=CHF(E)              0.01   \n",
       "(E)-HFC-1234ze                            trans-CF3CH=CHF              0.04   \n",
       "(Trifluoromethyl) sulphur pentafluoride            SF5CF3            800.00   \n",
       "(Z)-HFC-1225ye                               CF3CF=CHF(Z)              0.02   \n",
       "...                                                   ...               ...   \n",
       "Sevoflurane (HFE-347mmz1)                   (CF3)2CHOCH2F              2.20   \n",
       "Sulphur hexafluoride                                  SF6           3200.00   \n",
       "Sulphuryl fluoride                                  SO2F2             36.00   \n",
       "Trifluoro(fluoromethoxy)methane                  CH2FOCF3              4.40   \n",
       "Trifluoromethyl formate                           HCOOCF3              3.50   \n",
       "\n",
       "                                         Radiative Efficiency (W m-2 ppb-1)  \\\n",
       "Acronym, Common Name or Chemical Name                                         \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene                                 0.04   \n",
       "(E)-HFC-1225ye                                                         0.01   \n",
       "(E)-HFC-1234ze                                                         0.04   \n",
       "(Trifluoromethyl) sulphur pentafluoride                                0.59   \n",
       "(Z)-HFC-1225ye                                                         0.02   \n",
       "...                                                                     ...   \n",
       "Sevoflurane (HFE-347mmz1)                                              0.32   \n",
       "Sulphur hexafluoride                                                   0.57   \n",
       "Sulphuryl fluoride                                                     0.20   \n",
       "Trifluoro(fluoromethoxy)methane                                        0.33   \n",
       "Trifluoromethyl formate                                                0.31   \n",
       "\n",
       "                                         AGWP 20-year (W m-2 yr kg-1)  \\\n",
       "Acronym, Common Name or Chemical Name                                   \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene                   1.370000e-13   \n",
       "(E)-HFC-1225ye                                           7.250000e-15   \n",
       "(E)-HFC-1234ze                                           8.740000e-14   \n",
       "(Trifluoromethyl) sulphur pentafluoride                  3.360000e-10   \n",
       "(Z)-HFC-1225ye                                           2.140000e-14   \n",
       "...                                                               ...   \n",
       "Sevoflurane (HFE-347mmz1)                                1.980000e-11   \n",
       "Sulphur hexafluoride                                     4.370000e-10   \n",
       "Sulphuryl fluoride                                       1.710000e-10   \n",
       "Trifluoro(fluoromethoxy)methane                          6.820000e-11   \n",
       "Trifluoromethyl formate                                  5.370000e-11   \n",
       "\n",
       "                                        GWP 20-year  \\\n",
       "Acronym, Common Name or Chemical Name                 \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene            5   \n",
       "(E)-HFC-1225ye                                   <1   \n",
       "(E)-HFC-1234ze                                    4   \n",
       "(Trifluoromethyl) sulphur pentafluoride       13500   \n",
       "(Z)-HFC-1225ye                                   <1   \n",
       "...                                             ...   \n",
       "Sevoflurane (HFE-347mmz1)                       795   \n",
       "Sulphur hexafluoride                          17500   \n",
       "Sulphuryl fluoride                             6840   \n",
       "Trifluoro(fluoromethoxy)methane                2730   \n",
       "Trifluoromethyl formate                        2150   \n",
       "\n",
       "                                         AGWP 100-year (W m-2 yr kg-1)  \\\n",
       "Acronym, Common Name or Chemical Name                                    \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene                    1.370000e-13   \n",
       "(E)-HFC-1225ye                                            7.250000e-15   \n",
       "(E)-HFC-1234ze                                            8.740000e-14   \n",
       "(Trifluoromethyl) sulphur pentafluoride                   1.600000e-09   \n",
       "(Z)-HFC-1225ye                                            2.140000e-14   \n",
       "...                                                                ...   \n",
       "Sevoflurane (HFE-347mmz1)                                 1.980000e-11   \n",
       "Sulphur hexafluoride                                      2.160000e-09   \n",
       "Sulphuryl fluoride                                        3.760000e-10   \n",
       "Trifluoro(fluoromethoxy)methane                           6.890000e-11   \n",
       "Trifluoromethyl formate                                   5.390000e-11   \n",
       "\n",
       "                                        GWP 100-year  AGTP 20-year (K kg-1)  \\\n",
       "Acronym, Common Name or Chemical Name                                         \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene             1           1.090000e-15   \n",
       "(E)-HFC-1225ye                                    <1           5.770000e-17   \n",
       "(E)-HFC-1234ze                                    <1           6.980000e-16   \n",
       "(Trifluoromethyl) sulphur pentafluoride        17400           9.930000e-12   \n",
       "(Z)-HFC-1225ye                                    <1           1.700000e-16   \n",
       "...                                              ...                    ...   \n",
       "Sevoflurane (HFE-347mmz1)                        216           2.060000e-13   \n",
       "Sulphur hexafluoride                           23500           1.290000e-11   \n",
       "Sulphuryl fluoride                              4090           4.580000e-12   \n",
       "Trifluoro(fluoromethoxy)methane                  751           9.590000e-13   \n",
       "Trifluoromethyl formate                          588           6.730000e-13   \n",
       "\n",
       "                                        GTP 20-year  AGTP 50-year (K kg-1)  \\\n",
       "Acronym, Common Name or Chemical Name                                        \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene            2           1.540000e-16   \n",
       "(E)-HFC-1225ye                                   <1           8.140000e-18   \n",
       "(E)-HFC-1234ze                                   <1           9.820000e-17   \n",
       "(Trifluoromethyl) sulphur pentafluoride       14500           1.100000e-11   \n",
       "(Z)-HFC-1225ye                                   <1           2.400000e-17   \n",
       "...                                             ...                    ...   \n",
       "Sevoflurane (HFE-347mmz1)                       302           8.240000e-13   \n",
       "Sulphur hexafluoride                          18900           1.470000e-11   \n",
       "Sulphuryl fluoride                             6690           2.550000e-12   \n",
       "Trifluoro(fluoromethoxy)methane                1400           9.270000e-14   \n",
       "Trifluoromethyl formate                         984           6.850000e-14   \n",
       "\n",
       "                                        GTP 50-year  AGTP 100-year (K kg-1)  \\\n",
       "Acronym, Common Name or Chemical Name                                         \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene           <1            1.120000e-16   \n",
       "(E)-HFC-1225ye                                   <1            5.950000e-18   \n",
       "(E)-HFC-1234ze                                   <1            7.180000e-17   \n",
       "(Trifluoromethyl) sulphur pentafluoride       17800            1.110000e-11   \n",
       "(Z)-HFC-1225ye                                   <1            1.760000e-17   \n",
       "...                                             ...                     ...   \n",
       "Sevoflurane (HFE-347mmz1)                        38            1.640000e-14   \n",
       "Sulphur hexafluoride                          23800            1.540000e-11   \n",
       "Sulphuryl fluoride                             4140            9.010000e-13   \n",
       "Trifluoro(fluoromethoxy)methane                 150            5.720000e-14   \n",
       "Trifluoromethyl formate                         111            4.470000e-14   \n",
       "\n",
       "                                        GTP 100-year  Molar mass (g mol-1)  \\\n",
       "Acronym, Common Name or Chemical Name                                        \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene            <1                130.49   \n",
       "(E)-HFC-1225ye                                    <1                132.03   \n",
       "(E)-HFC-1234ze                                    <1                114.04   \n",
       "(Trifluoromethyl) sulphur pentafluoride        20200                196.06   \n",
       "(Z)-HFC-1225ye                                    <1                132.03   \n",
       "...                                              ...                   ...   \n",
       "Sevoflurane (HFE-347mmz1)                         30                200.05   \n",
       "Sulphur hexafluoride                           28200                146.06   \n",
       "Sulphuryl fluoride                              1650                102.06   \n",
       "Trifluoro(fluoromethoxy)methane                  105                118.03   \n",
       "Trifluoromethyl formate                           82                114.02   \n",
       "\n",
       "                                         Correction  \n",
       "Acronym, Common Name or Chemical Name                \n",
       "(E)-1-Chloro-3,3,3-trifluoroprop-1-ene          1.0  \n",
       "(E)-HFC-1225ye                                  1.0  \n",
       "(E)-HFC-1234ze                                  1.0  \n",
       "(Trifluoromethyl) sulphur pentafluoride         1.0  \n",
       "(Z)-HFC-1225ye                                  1.0  \n",
       "...                                             ...  \n",
       "Sevoflurane (HFE-347mmz1)                       1.0  \n",
       "Sulphur hexafluoride                            1.0  \n",
       "Sulphuryl fluoride                              1.0  \n",
       "Trifluoro(fluoromethoxy)methane                 1.0  \n",
       "Trifluoromethyl formate                         1.0  \n",
       "\n",
       "[203 rows x 15 columns]"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a_x_i = pd.read_csv('GHG_property_table.csv',index_col=0)\n",
    "a_x_i['Lifetime (Years)'] = pd.to_numeric(a_x_i['Lifetime (Years)'],errors='coerce')\n",
    "a_x_i"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "12.4"
      ]
     },
     "execution_count": 38,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a_x_i['Lifetime (Years)']['Methane']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 107,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "2000 2200 40000000000000.0 40000000000000.0 0.0 40000000000000.0\n",
      "(201,)\n",
      "2000 2200 400000000000.0 400000000000.0 0.0 400000000000.0\n",
      "(201,)\n"
     ]
    }
   ],
   "source": [
    "atm_mass   = 5.1480e18 # kg\n",
    "atm_mol    = 29.80 # g/mol\n",
    "\n",
    "# CO2\n",
    "alpha_CO2 = [0.217278,0.224037,0.282381,0.276303]\n",
    "tau_CO2   = [394.409,36.5393,4.30365]\n",
    "rf_ppb_CO2 = a_x_i['Radiative Efficiency (W m-2 ppb-1)']['Carbon dioxide']\n",
    "MM_CO2 = a_x_i['Molar mass (g mol-1)']['Carbon dioxide']\n",
    "\n",
    "# CH4\n",
    "alpha_CH4 = 1.\n",
    "tau_CH4   = a_x_i['Lifetime (Years)']['Methane']\n",
    "rf_ppb_CH4 = a_x_i['Radiative Efficiency (W m-2 ppb-1)']['Methane']\n",
    "MM_CH4 = a_x_i['Molar mass (g mol-1)']['Methane']\n",
    "\n",
    "# profile in kg\n",
    "profile = ((2000, 20),\n",
    "           (2010, 30),\n",
    "           (2020, 35),\n",
    "           (2030, 38),\n",
    "           (2040, 40),\n",
    "           (2050, 40))\n",
    "\n",
    "profile = ((2000, 20),\n",
    "           (2010, 30),\n",
    "           (2020, 40))\n",
    "#\n",
    "#profile = ((2000, 20),\n",
    "#           (2020, 40))\n",
    "\n",
    "#profile = ((2000, 50),\n",
    "#           (2050, 30))\n",
    "\n",
    "profile_CO2 = ((2000, 40e12), # in Gt\n",
    "          (2200, 40e12))\n",
    "profile_CH4 = ((2000, .4e12), # in Gt\n",
    "          (2200, .4e12))\n",
    "\n",
    "a_CO2,b_CO2   = concentration(profile_CO2,\n",
    "                      alpha=alpha_CO2,\n",
    "                      tau=tau_CO2,\n",
    "                      c_0=275,\n",
    "                             t_max=2200)\n",
    "\n",
    "a_CH4,b_CH4   = concentration(profile_CH4,\n",
    "                      alpha=alpha_CH4,\n",
    "                      tau=tau_CH4,\n",
    "                      c_0=275,\n",
    "                             t_max=2200)\n",
    "\n",
    "results = pd.concat([a_CO2,\n",
    "                     a_CH4,\n",
    "          rf_ppb_CO2/MM_CO2 * atm_mol/atm_mass * 1e9 * b_CO2,\n",
    "          \n",
    "          rf_ppb_CH4/MM_CH4 * atm_mol/atm_mass * 1e9 * b_CH4],\n",
    "          keys = pd.MultiIndex.from_product([['Emissions','Radiative forcing'],['CO2','CH4']]),\n",
    "    axis=1)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 108,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <th colspan=\"2\" halign=\"left\">Emissions</th>\n",
       "      <th colspan=\"2\" halign=\"left\">Radiative forcing</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <th>CO2</th>\n",
       "      <th>CH4</th>\n",
       "      <th>CO2</th>\n",
       "      <th>CH4</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>3.602572e-11</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2001</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>0.069638</td>\n",
       "      <td>5.034378e-02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2002</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>0.135008</td>\n",
       "      <td>9.678697e-02</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2003</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>0.196889</td>\n",
       "      <td>1.396318e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2004</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>0.255902</td>\n",
       "      <td>1.791570e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2196</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>6.389718</td>\n",
       "      <td>6.497729e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2197</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>6.415285</td>\n",
       "      <td>6.497729e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2198</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>6.440825</td>\n",
       "      <td>6.497729e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2199</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>6.466337</td>\n",
       "      <td>6.497729e-01</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2200</th>\n",
       "      <td>4.000000e+13</td>\n",
       "      <td>4.000000e+11</td>\n",
       "      <td>6.491816</td>\n",
       "      <td>6.497729e-01</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>201 rows × 4 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         Emissions               Radiative forcing              \n",
       "               CO2           CH4               CO2           CH4\n",
       "2000  4.000000e+13  4.000000e+11          0.000000  3.602572e-11\n",
       "2001  4.000000e+13  4.000000e+11          0.069638  5.034378e-02\n",
       "2002  4.000000e+13  4.000000e+11          0.135008  9.678697e-02\n",
       "2003  4.000000e+13  4.000000e+11          0.196889  1.396318e-01\n",
       "2004  4.000000e+13  4.000000e+11          0.255902  1.791570e-01\n",
       "...            ...           ...               ...           ...\n",
       "2196  4.000000e+13  4.000000e+11          6.389718  6.497729e-01\n",
       "2197  4.000000e+13  4.000000e+11          6.415285  6.497729e-01\n",
       "2198  4.000000e+13  4.000000e+11          6.440825  6.497729e-01\n",
       "2199  4.000000e+13  4.000000e+11          6.466337  6.497729e-01\n",
       "2200  4.000000e+13  4.000000e+11          6.491816  6.497729e-01\n",
       "\n",
       "[201 rows x 4 columns]"
      ]
     },
     "execution_count": 108,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 109,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<AxesSubplot:title={'center':'Warming (W m-2)'}>"
      ]
     },
     "execution_count": 109,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEICAYAAABF82P+AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAAWJElEQVR4nO3df5xddX3n8dcHiERKCBj8RSaYkF9LQjDCgG6sW4q7SLAB1FVRa3FXSakPVmv1obEt2xTFByi73bKUYmgqlgVSH1iEVFiBLj+KwsIEkpSYTUxCMCPmJzWAgUDws3/cM3gZvjMJmTtzTyav5+Mxj7n3e8/5ns89+ea+7/ece89EZiJJUm8HtLsASVI9GRCSpCIDQpJUZEBIkooMCElSkQEhSSoyIFQ7EXFVRFw4gPX/OCL+ppU17eF2b4iIs6vbn4iI+4a6hoGKiH+IiNPbXYfqIfwehFohItYDbwRebGq+JjMvaE9FQysijgcWAdMzMyPiE8CnMvM397K/1wBPAOMz85nWVbrb7Z4M/HVmnjhU21R9HdTuAjSszMnMO9tdRJv8PnBdtu4d178Dlg5VOERE0HjD+GBEHBYRnZnZNRTbVn15iEmDrjrc8sOI+IuI+EVErIuIWVX7hojYHBHnNi1/TUR8tbp9ZET8Y7XekxHxzxFxQPXYlyLiZxHxdESsioh3V+3zI+J/NfV3ZkSsqPq4OyKObXpsfUR8ISKWR8T2iPj7iBi5u20XzAbu6WcffCMi7ouI0RFxYET8t4jYGhGPRcQFEZER0fyG7Qzg1oj4YEQs6dXX5yPie9XtgyPisoj4aURsqg7PvbZ67Iiq/i0R8a/V7Y6mfu6OiIsj4ofADuCY6qG7gff29Vy0/zAgNFTeDiwHxgDX0zgccxIwCfhd4IqIOLSw3ueBbuD1NA5h/TGQETEVuAA4KTNHAe8B1vdeOSKmADcAf1j1cSuwuDqE0+NDwOnABOB44BP9bbuwjd+o1l1VeOyAiLi66ve0zNwOnEcjUGYCJwBnF573GcD3gVuACc2hRmN/XVvdvhSYUvU1CRgL/NfqsQOAbwFvAY4GngWu6LWdjwNzgVHA41XbSuCthZq0nzEg1Erfq95t9/yc1/TYY5n5rcx8Efh7YBxwUWbuzMzbgedpvMD19gLwZuAtmflCZv5zdRjnReBgYFpEjMjM9Zm5trD+h4HvZ+YdmfkCcBnwWmBW0zKXZ+YTmfkksJjGi21/2+7t8Or3073aR9AIp9fROPy2o2r/EPCXmdmdmf8KXNK8UkQcA4zIzFWZubPaX79bPTYdGA/8Y3VY6Dzgc5n5ZGY+DXwNOAcgM7dl5nczc0f12MXAb/Wq8ZrMXJGZu6r90/M8Dkf7PQNCrXR2Zh7e9HN102Obmm4/C5CZvdtKM4hvAGuA26tDU/OqddfQmBXMBzZHxKKIOKqw/lH8+p0xmfkrYAONd9o9Njbd3tFUR3HbBb+ofo/q1T4JOAv488x8vldNG5rub+Dl3ktjptPj28BHq0D4OPCdKjheDxwCLOkJZeB/V+1ExCER8c2IeDwingLuBQ6PiAP72XbP8/hFoV37GQNCtZaZT2fm5zPzGGAO8Ec95xoy8/rqU0JvoXHo59JCF09UjwMvnYwdB/xsINvutdwvgbU0DvU0Wwn8J+C26pBYj58DHU33x/Var+fwUk//D9CYYb0L+Ci/Pry0lUawTm8K5dGZ2RNwnwemAm/PzMNonPgGiObyC0/9WGBZoV37GQNCtRYRvxMRk6oX9qdoHFp6MSKmRsSpEXEw8ByNF8oXC118B3hvRLw7IkbQeNHcCfxob7fdx+K38srDN2TmDTTOXdwZERObavpsRIyNiMOBLzVt87XAyTROFDf7OxrnD3Zl5n1V378Crgb+IiLeUK0/NiLeU60zisZ++UVEvA74s90958pvAbft4bIaxgwItdLiiHim6eemFvQ5GbgTeAa4H7gyM++mcf7hEhrvojcCb6DxQvwymbmKxvH7/1ktO4fG+YDney/7KrZdsgD4WBUmvWv4NnAR8H8iYjyNF/XbaZy0f4RGuOyiET7vBu7PzOd6dXMtcBy/nj30+BKNw2APVIeR7qQxawD4HzTOt2wFHqBx+KlfEXES8MvMfHB3y2r484tyUotExPU0zg9871WuNxu4KjPfEhFXAo9m5pW9lnktsBk4ITN/0qqaC7V8F1iYmbfudmENe35RTmqRzPzonixXvdj/No1ZxBtpHPrpmW0tpfFJqt7+AHhoMMMBIDM/MJj9a9/iDEIaYhFxCI0v1f0bGucIvg98NjOf6mP59TROLJ+dmY8MVZ2SASFJKvIktSSpqJbnII488sgcP358u8uQpH3KkiVLtmbm61vVX60CIiLmAHMmTZpEV5cXkpSkVyMiHt/9UnuuVoeYMnNxZs4dPXp0u0uRpP1erQJCklQfBoQkqciAkCQVGRCSpCIDQpJUZEBIkopq9T2IHuu2/JIPf/P+dpchSfu1Ws0gImJORCzYtWtXu0uRpP1eLS/W19nZmX6TWpJenYhYkpmdreqvVjMISVJ9GBCSpCIDQpJUZEBIkooMCElSkQEhSSoyICRJRQaEJKnIgJAkFRkQkqQiA0KSVFSrgOi5WN/27dvbXYok7fdqFRCZuTgz544ePbrdpUjSfq9WASFJqg8DQpJUZEBIkooMCElSkQEhSSoyICRJRQaEJKnIgJAkFRkQkqQiA0KSVGRASJKKDAhJUpEBIUkqMiAkSUUGhCSpqFYB4R8MkqT6qFVA+AeDJKk+ahUQkqT6MCAkSUUGhCSpyICQJBUZEJKkIgNCklRkQEiSigwISVKRASFJKjIgJElFBoQkqciAkCQVGRCSpCIDQpJUZEBIkooMCElSkQEhSSoyICRJRQaEJKnIgJAkFRkQkqQiA0KSVGRASJKKDAhJUlHLAyIijomIhRFxY1PbsRFxVUTcGBF/0OptSpJab48CIiL+NiI2R8SjvdpPj4hVEbEmIuYBZOa6zPxk83KZuTIzzwc+BHS2qnhJ0uDZ0xnENcDpzQ0RcSDwV8BsYBrwkYiY1lcHEXEmcB/wT3tVqSRpSO1RQGTmvcCTvZpPBtZUM4bngUXAWf30cUtmzgI+trfFSpKGzkDOQYwFNjTd7wbGRsSYiLgKeFtEfBkgIk6JiMsj4pvAraXOImJuRHRFRNeWLVsGUJYkqRUOGsC6UWjLzNwGnN+r8W7g7v46y8wFwAKAzs7OHEBdkqQWGMgMohsY13S/A3hiYOVIkupiIAHxEDA5IiZExGuAc4BbWlOWJKnd9vRjrjcA9wNTI6I7Ij6ZmbuAC4AfACuB72TmioEUExFzImLB9u3bB9KNJKkFIrN+h/s7Ozuzq6ur3WVI0j4lIpZkZsu+a+alNiRJRQaEJKnIgJAkFdUqIDxJLUn1UauAyMzFmTl39OjR7S5FkvZ7tQoISVJ9GBCSpCIDQpJUZEBIkopqFRB+ikmS6qNWAeGnmCSpPmoVEJKk+jAgJElFBoQkqciAkCQVGRCSpKJaBYQfc5Wk+jio3QU0y8zFwOLOzs7z2l2LpH3XCy+8QHd3N88991y7SxkUI0eOpKOjgxEjRgzqdmoVEJLUCt3d3YwaNYrx48cTEe0up6Uyk23bttHd3c2ECRMGdVu1OsQkSa3w3HPPMWbMmGEXDgARwZgxY4ZkdmRASBqWhmM49Biq52ZASNIg2bhxI+eccw4TJ05k2rRpnHHGGaxevZoVK1Zw6qmnMmXKFCZPnsxXvvIVMhOA6667juOPP57jjz+eWbNmsWzZsrbVb0BI0iDITN73vvdxyimnsHbtWn784x/zta99jU2bNnHmmWcyb948Vq9ezbJly/jRj37ElVdeCcCECRO45557WL58ORdeeCFz585t23MwICRpENx1112MGDGC888//6W2mTNnsnr1at75zndy2mmnAXDIIYdwxRVXcMkllwAwa9YsjjjiCADe8Y530N3dPfTFV2r1KaaImAPMmTRpUrtLkTRM/PniFfz4iada2ue0ow7jz+ZM73eZRx99lBNPPPEV7StWrHhF+8SJE3nmmWd46qmnOOyww15qX7hwIbNnz25N0XuhVjMIL/ctabjLzD5PMje333XXXSxcuJBLL710qEp7hVrNICSp1Xb3Tn+wTJ8+nRtvvLHYfu+9976sbd26dRx66KGMGjUKgOXLl/OpT32K2267jTFjxgxJvSW1mkFI0nBx6qmnsnPnTq6++uqX2h566CEmT57Mfffdx5133gnAs88+y2c+8xm++MUvAvDTn/6U97///Vx77bVMmTKlLbX3MCAkaRBEBDfddBN33HEHEydOZPr06cyfP5+jjjqKm2++ma9+9atMnTqVGTNmcNJJJ3HBBRcAcNFFF7Ft2zY+/elPM3PmTDo7O9v3HHo+e1snnZ2d2dXV1e4yJO2jVq5cybHHHtvuMgZV6TlGxJLMbFmiOIOQJBUZEJKkIgNCklRUq4DwDwZJUn3UKiD8opwk1UetAkKSVB8GhCQNkr4u933ccce9bLn58+dz2WWXvaztsssuIyLYunXrUJb8Ml5qQ5IGQc/lvs8991wWLVoEwNKlS9m0adNu192wYQN33HEHRx999GCX2S9nEJI0CPq63Pe4ceN2u+7nPvc5vv71r7f9r+I5g5A0vN02Dzb+S2v7fNMMmH1Jv4v0dblvgLVr1zJz5syX7m/cuJEvfOELANxyyy2MHTuWt771rS0rd28ZEJI0xCZOnMjSpUtfuj9//nwAduzYwcUXX8ztt9/ensJ6MSAkDW+7eac/WPq63Hd/1q5dy2OPPfbS7KG7u5sTTjiBBx98kDe96U2DUWa/PAchSYOgr8t9P/74432uM2PGDDZv3sz69etZv349HR0dPPzww20JBzAgJGlQ9He5732Fl/uWNOx4ue/WqNUMwmsxSVJ91CogvBaTJNVHrQJCklQfBoSkYamO51dbZaiemwEhadgZOXIk27ZtG5YhkZls27aNkSNHDvq2/KKcpGGno6OD7u5utmzZ0u5SBsXIkSPp6OgY9O0YEJKGnREjRjBhwoR2l7HP8xCTJKnIgJAkFRkQkqQiA0KSVGRASJKKDAhJUpEBIUkqMiAkSUUGhCSpyICQJBUZEJKkoloFhH9RTpLqo1YB4V+Uk6T6qFVASJLqw4CQJBUZEJKkIgNCklRkQEiSigwISVKRASFJKjIgJElFBoQkqciAkCQVGRCSpCIDQpJUZEBIkooMCElSkQEhSSoyICRJRQaEJKnIgJAkFRkQkqQiA0KSVGRASJKKDAhJUpEBIUkqMiAkSUUtD4iIOCYiFkbEjU1tZ0fE1RFxc0Sc1uptSpJab48CIiL+NiI2R8SjvdpPj4hVEbEmIuYBZOa6zPxk83KZ+b3MPA/4BPDhFtUuSRpEezqDuAY4vbkhIg4E/gqYDUwDPhIR03bTz59W60iSam6PAiIz7wWe7NV8MrCmmjE8DywCziqtHw2XArdl5sMDKViSNDQGcg5iLLCh6X43MDYixkTEVcDbIuLL1WP/Bfj3wH+MiPNLnUXE3IjoioiuLVu2DKAsSVIrHDSAdaPQlpm5DTi/V+PlwOX9dZaZC4AFAJ2dnTmAuiRJLTCQGUQ3MK7pfgfwxMDKkSTVxUAC4iFgckRMiIjXAOcAt7SmLElSu+3px1xvAO4HpkZEd0R8MjN3ARcAPwBWAt/JzBUDKSYi5kTEgu3btw+kG0lSC0Rm/Q73d3Z2ZldXV7vLkKR9SkQsyczOVvXnpTYkSUUGhCSpyICQJBXVKiA8SS1J9VGrgMjMxZk5d/To0e0uRZL2e7UKCElSfRgQkqQiA0KSVGRASJKKahUQfopJkuqjVgHhp5gkqT5qFRCSpPowICRJRQaEJKnIgJAkFRkQkqSiWgWEH3OVpPqoVUD4MVdJqo9aBYQkqT4MCElSkQEhSSoyICRJRQaEJKmoVgHhx1wlqT5qFRB+zFWS6qNWASFJqg8DQpJUZEBIkooMCElSkQEhSSoyICRJRQaEJKnooHYXULT1J/Ct97a7Cknar9VqBtHzTepdL+5qdymStN+LzGx3Da/Q2dmZXV1d7S5DkvYpEbEkMztb1V+tZhCSpPowICRJRQaEJKnIgJAkFRkQkqQiA0KSVGRASJKKDAhJUlEtvygXEU8Dq9pdxx44Etja7iL2gHW2lnW21r5Q575QI8DUzBzVqs7qeS0mWNXKbwMOlojoss7Wsc7Wss7W2RdqhEadrezPQ0ySpCIDQpJUVNeAWNDuAvaQdbaWdbaWdbbOvlAjtLjOWp6kliS1X11nEJKkNjMgJElFQxIQETEuIu6KiJURsSIiPlu1vy4i7oiIn1S/j2ha58sRsSYiVkXEe5raT4yIf6keuzwiol11RsR/iIglVT1LIuLUpr7urmpfWv28oY11jo+IZ5tquaqprzrtz4811bg0In4VETOrx9qxPz9Y3f9VRHT2WqdO47NYZw3HZ191Dvn43Isa6zY2vxER/y8ilkfETRFxeNM6rRubmTnoP8CbgROq26OA1cA04OvAvKp9HnBpdXsasAw4GJgArAUOrB57EPi3QAC3AbPbWOfbgKOq28cBP2vq626gsyb7czzwaB991WZ/9lp3BrCuzfvzWGBq723XcHz2VWfdxmdfdQ75+Hy1NdZwbJ4GHFS1X8ogvXYOyQwiM3+emQ9Xt58GVgJjgbOAb1eLfRs4u7p9FrAoM3dm5mPAGuDkiHgzcFhm3p+NZ/x3TesMeZ2Z+UhmPlG1rwBGRsTBraqnVXX2pW77s5ePADe0qpb+9FVnZq7MzNI3+ms1Pvuqs27js5/9WTSY+3OANdZhbN6embuqxR4AOqrbLR2bQ34OIiLG03hn83+BN2bmz6GxI4CeqdlYYEPTat1V29jqdu/2dtXZ7APAI5m5s6ntW9WU88JWHmrYyzonRMQjEXFPRLyraqvz/vwwr/xPONT7sy91G597og7jsz9tG597sS/rNjb/M40ZAbR4bA5pQETEocB3gT/MzKf6W7TQlv20t9SrqLNn+ek0pnm/39T8scycAbyr+vl4G+v8OXB0Zr4N+CPg+og4jPruz7cDOzLz0abmOu1Px2dr62zb+NzXx2ZE/AmwC7iup6mw+l6PzSELiIgYQeMJXpeZ/1A1b6qmPj3Tyc1Vezcwrmn1DuCJqr2j0N6uOomIDuAm4Pcyc21Pe2b+rPr9NHA9cHK76qymm9uq20toHJecQg33Z+Ucer1Da9P+7Evdxmd/y9dpfBa1a3y+2n1Zqc3YjIhzgd+hEVA9L/YtHZtD9SmmABYCKzPzvzc9dAtwbnX7XODmpvZzIuLgiJgATAYerA5HPB0R76j6/L2mdYa8zuqTA98HvpyZP2zq56CIOLK6PYLGP2LzO46hrvP1EXFgdfsYGvtzXd32Z7XOAcAHgUVNbe3an32p2/jsa/nDqdf47Gv5IR+fe/FvXquxGRGnA18CzszMHU2rtHZs5iCcee/9A/wmjenMcmBp9XMGMAb4J+An1e/XNa3zJzTeSayi6Ww70EnjH2AtcAXVt8HbUSfwp8Avm5ZdSuN4+m8AS6p+VgB/SfVJgjbV+YGqjmXAw8CcOu7Pap1TgAd69dOu/fk+Gu+8dgKbgB/UdHwW66zh+OyrziEfn3v5b16nsbmGxrmGnrarBmNseqkNSVKR36SWJBUZEJKkIgNCklRkQEiSigwISVKRASFJKjIgJElF/x850X8oS+jTMQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEICAYAAACktLTqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAAArdklEQVR4nO3dd3yV5f3/8deHEBJG2GETEgJhyCYOQIZUrLhnXV+lWkSr0lrrr1Vrv1pb/aq1rXugqDhwS0udqOwhS/bMYIURIKwMsq/fH+dEI5IFZ+Xk/Xw8ziPnXPfI5xxu3rnPfV/3dZtzDhERCS/1gl2AiIj4nsJdRCQMKdxFRMKQwl1EJAwp3EVEwpDCXUQkDCncpdYws8/NbJyf1h1lZuvNrJ0/1u9LZtbWzDaYWVSwa5HQpXCXCpnZvWb22TFtKRW0Xe3vepxzY51zU/y0+gnAXOfcHjMbYmZHzCyibKKZvVxB24v+KMbMnvB+rtlmttHMbiib5pzLBGZ5axY5LoW7VGYuMKws0Lx7tZHAoGPaunnnrTYzq+/jWk/WLcCb3ufLgAhgULnpw4Fdx7SNoIbvuwZygQuBZsA44CkzG1pu+tvemkWOS+EulVmKJ8wHeF+PwLPHuOmYtjTn3C4zu9F7uCDbzNLN7PvwMbNRZpZhZn80sz3Aa2b2oJl9YGZveZdZY2ZJ3m8Me81sh5mdU24ds81svPf5L81svncP96CZbTGzseXmTTCzud71fm1mz5nZW8d7k2YWByQCiwGcc0XAt973hpm1ARoA7x3TlkQF4W5mW83s/5nZajPLNbPJ3sMpn5erqUVFH7xz7gHn3EbnXKlzbjEwDxhSbpbFQFcz61LROqRuU7hLhZxzhXhCZIS3aQSekJl/TFtZwO0FLgCaAjcC/zKz8nu67YCWQBd+OKRwIZ495hbACuBLPNtlR+Ah4KVKSjwdzx+a1sDjwGQzM++0qcASoBXwIHB9JevpC6Q754rLtc095j3OP8773uKcy6hkvZcDY/D8EbgQ+By4z1tvPeA3lSz7PTNrCJwKrCtr89aaCvSvzjqk7lG4S1Xm8EOgDccT7vOOaZsD4Jz71DmX5jzmADO808uUAg845wqcc0e9bfOcc196w+oDIBZ41Lv3/C4Qb2bNK6htm3PuZedcCTAFaA+09e6Jnwr8r3Ou0Dk3H5heyXtsDmQf532f6f1jUfa+FwFnlGubU8k6AZ5xzmU653Z6l1/snFvhnCsApgEDq1i+zIvAKjx/+MrL9tYu8hMKd6nKXDwh1wKIdc6lAAuBod62Pt55MLOxZvatmR0ws0PAeXj2Usvsc87lH7P+zHLPjwL7vWFd9hqgSQW17Sl74pzLKzdvB+BAuTaAHZW8x4NAzDFt33rX1QfvNxbnXI53PWVtVR1vP/a9Hfu6CYCZvWhmOd7HfeVXYGZ/9/6+X7ifjvIXAxyqogapoxTuUpVFeE7qTQAWADjnjuA5uTgB2OWc2+LtlvcR8ATQ1jnXHPgMsHLrCtQQpLuBlmbWqFxb50rmX43n+PX3J3m9f4SW4jnM1N45t9E7aZ63rR8+OpnqnLvVOdfE+3ikrN3M/gKMBc7xfuaUm1Yfz4nsVb6oQcKPwl0q5T18sgy4C0+wlZnvbSsLuAZAFLAPKPae3DyHIHDObcNT84Nm1sDMhuA55l3R/BlACnDaMZPmAnfi+aZSZr63bY9zLs2HZf+Imd0LXAuMcc5lHWeW04Ct3vcq8hMKd6mOOUAbPMFWZp63bS6Acy4bzwnC9/Ec5riWyo9z+9t1eHqXZAF/w9PTpaCS+V/ipyddj/e+51PuffvRI0AckFLBIZvr8ByLFzku0806pC4ws/eAjc65ByqYHoWnt87PnHO7A1pcDXm7Yc4BBh7nHIYIoHCXMGVmpwIHgC14Dg/9GxjinFsRzLpEAiXUrhIU8ZV2wMd4+rlnAL9WsEtdoj13EZEwpBOqIiJhyC+HZVq3bu3i4+P9sWoRkbC0fPny/c65WF+tzy/hHh8fz7Jly/yxahGRsGRmPr1mQYdlRETCkMJdRCQMKdxFRMJQwPq5FxUVkZGRQX5+eF5QFx0dTadOnYiMjAx2KSIigQv3jIwMYmJiiI+P54f7KYQH5xxZWVlkZGSQkJAQ7HJERAJ3WCY/P59WrVqFXbADmBmtWrUK228lIlL7BPSYezgGe5lwfm8iUvvohKqISJCt33Wk6plqqE6F+549e7j66qtJTEykd+/enHfeeWzevJl169YxevRokpKS6N69O3/9618pG3Pn7bffpl+/fvTr14+hQ4eyapVufCMiJ885x6K0LMa9uoTznp5X9QI1VGfC3TnHpZdeyqhRo0hLS2P9+vU88sgjZGZmctFFF3HPPfewefNmVq1axcKFC3n++ecBSEhIYM6cOaxevZo///nPTJgwIcjvRERqs9JSxxdr93DJ8wu55uVvWbvzMHefk+Tz31NnhvydNWsWkZGR3Hrrrd+3DRgwgMmTJzNs2DDOOcdzR7hGjRrx7LPPMmrUKG6//XaGDh36/fxnnHEGGRkZAa9dRGq/guIS/r1iJy/NTSd9Xy5xLRvx10v6cOXgTkRHRjDRx78vKOH+l/+u8/kxpt4dmvLAhadUOH3t2rUMHjz4J+3r1q37SXtiYiI5OTkcOXKEpk2bft8+efJkxo4d67uiRSTsZecX8c6S7Uyev4XMIwWc0qEpz1wzkLF92lE/wn8HT+rMnntFnHMV9nQp3z5r1iwmT57M/PnzjzuviEh5+7ILeG3BFt78dhvZ+cUMTWzF36/oz/DurQPSuy4o4V7ZHra/nHLKKXz44YfHbZ8798f3Ok5PT6dJkybExMQAsHr1asaPH8/nn39Oq1atAlKviNRO27JymTQ3nQ+WZ1BUUsrYPu24ZUQi/Ts3D2gddeaE6ujRoykoKODll1/+vm3p0qV0796d+fPn8/XXXwNw9OhRfvOb3/CHP/wBgO3bt3PZZZfx5ptvkpTk+5MeIhIe1u48zB1Tv+OsJ2bzwbIMLh/UkW/uGsnz1w0OeLBDHTosY2ZMmzaNO++8k0cffZTo6Gji4+N58skn+c9//sPEiRO5/fbbKSkp4frrr+eOO+4A4KGHHiIrK4vbbrsNgPr162usehEBvN0Z07N4YXYa81L2ExNVnwkjErlpWDxtmkYHtTa/3EM1OTnZHRuAGzZsoFevXj7/XaGkLrxHEfF0Z/x6QybPz05j5Y5DxMZEcdOwBK47I46m0Sc2eKCZLXfOJfuqxjqz5y4icrKKSkqZvnIXL85JI2VvDp1bNuRvl/ThCm93xlCicBcRqUJ+UQnvLd3BpLnp7Dx0lJ7tYnjq6gGc37e9X7szngyFu4hIBQ4fLeKtb7fx6vwtZOUWMrhLCx66+BRG92wT8oMFVivczaw58ArQB3DATc65RX6sS0QkaPZlFzB5/hbe/nYb2QXFjEyK5fazunFaQstgl1Zt1d1zfwr4wjl3hZk1ABr5sSYRkaDYcSCPl+am8f6yDIpLShnbtz2/HplIn47Ngl1ajVUZ7mbWFBgB/BLAOVcIFPq3LBGRwNm0J5sXZqfy39W7iTDj8sEdmTAikYTWjYNd2gmrzpmArsA+4DUzW2Fmr5jZT96xmU0ws2Vmtmzfvn0+L9QXKhryt0+fPj+a78EHH+SJJ574UdsTTzyBmbF///5AliwifrR820HGT1nKz5+cy4z1mdw0LJ65fziL/7usX60OdqjeYZn6wCBgonNusZk9BdwD/Ln8TM65ScAk8PRz93WhJ6tsyN9x48bx7rvvArBy5UoyMzOrXHbHjh189dVXxMXF+btMEfEz5xxzU/bz/KxUFm85QPNGkdx5dnfGDYmnReMGwS7PZ6oT7hlAhnNusff1h3jCvVapaMjfrVu3Vrns7373Ox5//HEuvvhiP1YoIv5U4h1H/YU5qazdeYR2TaO5//xeXHNaHI2jwq/jYJXvyDm3x8x2mFkP59wm4GfA+pP6rZ/fA3vWnNQqfqJdXxj7aIWTKxryFyAtLY0BAwZ8/3rPnj3cfffdAEyfPp2OHTvSv39/n5YrIoFRWFzKtBUZvDQnnfT9uSS0bsxjl/flkoEdiaofWhce+VJ1/1xNBN729pRJB270X0mBl5iYyMqVK79//eCDDwKQl5fHww8/zIwZM4JTmIicsNyCYt5Zsp1X5m1hz5F8TunQlOeuHcS5fdoRUS+0+6j7QrXC3Tm3EvDZmAeV7WH7S0VD/lYmLS2NLVu2fL/XnpGRwaBBg1iyZAnt2rXzR5kicpIOHy1iysKtvLZgCwfzijg9oSWPX9EvYOOoh4rwO9BUgdGjR3Pffffx8ssvc/PNNwOeIX/z8vIqXKZv377s3bv3+9fx8fEsW7aM1q1b+71eEamZrBzPhUdvLvJcePSznm247axuDO7SItilBUWdCffKhvwVkdprz+F8Js1NZ+qSbRQUl3Je3/bcPqobvTs0rXrhMKYhf32oLrxHkVCxPSuPF+ak8dHyDEqc45IBHfn1qES6tWkS7NJOiIb8FZE6LXVvNs/PSuM/q3YRYcaVyZ24dWQinVtqVJTyFO4iUius3XmY52en8vnaPUTXj+DGofHcPKIrbYN8x6NQFdBwd86F7dlqfxzeEhFYvu0Az85MZdamfcRE1ef2Ud246cwEWobR1aT+ELBwj46OJisri1atWoVdwDvnyMrKIjpaexAivuCcY2FaFs/OTGVRehYtGkVy9zlJXD8knmYNT+w2dnVNwMK9U6dOZGRkEKqDip2s6OhoOnXqFOwyRGo15xwzN+7lmZmprNxxiDYxUdx/fi+uPT2ORg10FLkmAvZpRUZGkpCQEKhfJyK1SEmp4/O1u3luVhobdh+hU4vQvTdpbaE/hSISNEUlpfxn5S6en51K+r5cEmMb848r+3PRgA5Ehui9SWsLhbuIBFx+UQkfLs/gxTlpZBw8Sq/2dWvcl0BQuItIwOQVFjN18XYmzU1nb3YBA+Oa89DFp3BWj9C/4XRto3AXEb/LKSjmjUVbeWXeFg7kFjI0sRVPXjWAIYnh13suVCjcRcRvjuQXMWXBViYv2MKhvCJG9Yhl4uhuDO7SMtilhT2Fu4j43KG8Ql5d4Bl2Nzu/mLN7tWXi6G7079w82KXVGQp3EfGZA7mFTJ6fzpSF28gpKObcU9pxx+hu9OnYLNil1TkKdxE5aftzCnh5XjpvLtrG0aISzuvbnomju9GzXd0edjeYFO4icsL2HsnnpbnpvL14G4XFpVzYvwN3nNWN7m1jgl1anadwF5Ea2334KC/NSWfqku2UlHrGUr/9rES6xtbOsdTDkcJdRKot42AeL85J4/2lGZQ6x+WDOnHbWYl0adU42KXJMRTuIlKl7Vl5PD87lQ+XZ2AGv0jurBtkhDiFu4hUaMv+XJ6blcq0FTuJqGdcd3oct4xMpEPzhsEuTapQrXA3s61ANlACFPvyPn8iEnpS9+bw3KxU/rNyJ5ER9Rg3JJ5bRuquR7VJTfbcz3LO7fdbJSISdJv2ZPPMzBQ+XbOb6PoRjB/elfHDE2gTo1CvbXRYRkTYuOcIT3+Twmdr9tC4QQS/HpnIr85MoFWTqGCXJieouuHugBlm5oCXnHOTjp3BzCYAEwDi4uJ8V6GI+M2mPdk89c1mPluzh5io+kwc3Y1fnZlA80a6P2ltV91wH+ac22VmbYCvzGyjc25u+Rm8gT8JIDk5WXeLFglhmzOzeeprz+GXJlH1+c1oz02nFerho1rh7pzb5f2518ymAacBcytfSkRCzebMbJ76JoXP1uymcQPtqYezKsPdzBoD9Zxz2d7n5wAP+b0yEfGZFG+of7pmN40iI7h9lCfUWzRWqIer6uy5twWmeQfUrw9Mdc594deqRMQnUvdm89Q3qXyyeheNIiO4bVQi48/sqlCvA6oMd+dcOtA/ALWIiI+k7s3h6W9S+O/qXTSM9PR+GT+8Ky0V6nWGukKKhJG0fZ5Qn77KE+q3jEhkwgiFel2kcBcJA2n7cnjGG+pR9SOYMKIrE4Z3VT/1OkzhLlKLpe/L4ZmZnmECoupHcPPwrkwYoVAXhbtIrbRlfy7PfJPCv1fupEH9eoz3hnprhbp4KdxFapGt+3N5emYK/17hCfVfnZnAhBGJxMYo1OXHFO4itcCOA3k8/U0KH6/YSWSEcdOwBG4ZqVCXiincRULY7sNHeXZmKu8v24GZMW5IPLeO6qpRGqVKCneRELQ3O58XZqfx9uLtOOe46tTO3HFWd9o1U6hL9SjcRULIgdxCXpqbxpSFWykqcVwxqBN3jO6m29lJjSncRULA4aNFvDIvnVfnbyGvqIRLBnTktz/rTnxr3XhaTozCXSSIcgqKeW3+Fl6el86R/GLO79ueO8/uTve2McEuTWo5hbtIEBwtLOGNRVt5cU4aB/OKOLtXW+4ak0TvDk2DXZqECYW7SADlF5XwzpLtPDcrjf05BYxMiuWuMUn079w82KVJmFG4iwRAYXEpHyzfwbMzU9l9OJ8hXVvx4v8MIjm+ZbBLkzClcBfxo+KSUj5esZOnv0kh4+BRBndpwT+u7M/Qbq2DXZqEOYW7iB+UlDo+Wb2LJ79OYcv+XPp1asbfLunDyKRYvDe+EfErhbuIDznnmLE+k3/M2MTmzBx6toth0vWDGdO7rUJdAkrhLuIjC1L38/iXm1i14xBdYxvz7LUDOa9Pe+rVU6hL4CncRU7Siu0HeWLGJhakZtGhWTSPX96PywZ1pH5EvWCXJnWYwl3kBG3ak80/ZmxixvpMWjVuwP9e0Jvrzogjqn5EsEsTUbiL1NT2rDye/Hoz01bupEmD+vx+TBI3nZlA4yj9d5LQUe2t0cwigGXATufcBf4rSSQ07T2Sz9MzU3h3yQ4i6hkTRnTl1yMTad5IN5+W0FOTXY3fAhsAXR8tdcqhvEJenJPO6wu3UFziuPq0zkwc3Z22TTX8roSuaoW7mXUCzgceBu7ya0UiISK3oJjXFmzhpbnp5BQUc8mAjtx5dne6tNJIjRL6qrvn/iTwB6DCoerMbAIwASAuLu6kCxMJloLiEqYu3s5zs1LZn1PImN5tufucHvRop5EapfaoMtzN7AJgr3NuuZmNqmg+59wkYBJAcnKy81WBIoFSNlTAU1+nsPPQUYYmtuLlG3owMK5FsEsTqbHq7LkPAy4ys/OAaKCpmb3lnPsf/5YmEhjOOT5fu4d/zNhE2r5c+ndqxmOX9+PM7hr/RWqvKsPdOXcvcC+Ad8/9bgW7hIuFaft57PONrMo4TPc2TXjxfwbz81M0VIDUfuqYK3XS+l1HeOyLjczZvI8OzaL5+xX9uGxQJyI0VICEiRqFu3NuNjDbL5WIBMCOA3n886vN/HvlTppGR/Kn83px/ZAuREfqqlIJL9pzlzrhQG4hz85M5a1vt2EGt45M5NaRiTRrGBns0kT8QuEuYS2vsJjJ8zx91fMKi/lFcmfuPDuJds10AZKEN4W7hKWiklLeW7qDp75JYV92Aef0bssfzu1Btzbqqy51g8JdwkpZt8a/f7mJLftzOTW+BS/+zyAGd9G9SqVuUbhL2FiYtp/HvvDcLCOpbRMmj0tmdM826tYodZLCXWq98t0a26tbowigcJdaLONgHv+Y8UO3xvvO68kNQ+LVrVEEhbvUQoePFvH8rFReW7gVQ90aRY5H4S61RmFxKW99u42nZ6Zw+GgRlw3sxO/PSaJD84bBLk0k5CjcJeSV9YB57IuNbMvKY1i3Vtx3Xi9O6dAs2KWJhCyFu4S05dsO8vCn6/luu6cHzGs3nsqopFj1gBGpgsJdQtK2rFwe/2ITn67ZTWxMFI9e1pcrBneifkS9YJcmUiso3CWkHMwt5JmZqbz57VYiI+px59nduXl4VxpHaVMVqQn9j5GQkF9UwhuLtvLMzFRyC4q56tTO/O7sJNroJtQiJ0ThLkFVWur47+pdPP7FJnYeOsqoHrHcO7aX7lcqcpIU7hI0i9OzeOSzDazKOEzv9k11azsRH1K4S8Bty8rlkc828OW6TNo3i+YfV/bn0oEdqafhAkR8RuEuAXMkv4jnZqby2oKt1I8wfj8mifHDu9KwgYYLEPE1hbv4XUmp492l2/nnjM0cyCvkikGduPvnPWirk6UifqNwF79akLqfv36yno17sjktviWvX9Cbvp10ZamIvyncxS/S9+XwyGcb+HrDXjq3bMgL1w3i3D7tdGWpSIBUGe5mFg3MBaK883/onHvA34VJ7XQ4r4inZ6YwZeFWoiMj+OO5PblxmIbhFQm06uy5FwCjnXM5ZhYJzDezz51z3/q5NqlFiktKmbpkO//6ajOHjhZx9amduWtMD2JjooJdmkidVGW4O+cckON9Gel9OH8WJbXL7E17efjTDaTszWFI11b8+YLe9O7QNNhlidRp1TrmbmYRwHKgG/Ccc26xX6uSWiF1bw5/+3Q9szftI75VIyZdP5gxvdvquLpICKhWuDvnSoABZtYcmGZmfZxza8vPY2YTgAkAcXFxvq5TQkh2fhFPf5PCawu20rBBBPef34sbhsTToL5GbBQJFTXqLeOcO2Rms4FzgbXHTJsETAJITk7WYZswVFrq+Oi7DB77YhNZuQVcldyZu3/eg9ZNdFxdJNRUp7dMLFDkDfaGwNnAY36vTELKyh2HeGD6OlbtOMSguOa8+stk+nVqHuyyRKQC1dlzbw9M8R53rwe875z7xL9lSajYm53P419s4sPlGbSJieKfv+jPJQM0DoxIqKtOb5nVwMAA1CIhpLC4lCkLt/LUNykUFJdw68hE7hjdjSa6aYZIraD/qfITszft5aFP1pO+L5fRPdvw5wt6k9C6cbDLEpEaULjL97Zl5fLXTzbw9YZMElo35tVfJjO6Z9tglyUiJ0DhLhwtLOG5WalMmptOZIRxz1jPkAFR9TVkgEhtpXCvw5xzzFifyUP/Xc/OQ0e5dGBH7hnbU0PxioQBhXsdtS0rlwenr2PWpn30aBvDexPO4PSurYJdloj4iMK9jskvKuGF2Wm8MCeNyHrG/ef3YtzQeCIjdHWpSDhRuNchMzdm8uD09Ww/kMdF/Tvwp/N76RCMSJhSuNcBOw7k8dAn6/lqfSaJsY2ZOv50hnZrHeyyRMSPFO5hrKC4hJfnpvPsrFQMTy+Ym4YlaIAvkTpA4R6m5m7exwPT17Flfy7n9W3H/ef3pkPzhsEuS0QCROEeZvYczuehT9bx2Zo9JLRuzJSbTmNkUmywyxKRAFO4h4mSUsebi7byxIzNFJWU8vsxSUwY2VUXIonUUQr3MLB252Hum7aG1RmHGd69NX+7pA9dWmksGJG6TOFei+UWFPOvrzbz6oIttGwcxdPXDOTCfu11mzsRUbjXVl+vz+R//7OWXYfzufb0OP748540axQZ7LJEJEQo3GuZ3YeP8uD0dXy5LpMebWP46NqBDO7SMthliUiIUbjXEiWljjcWbeWJLzdR4hx/PLcn44cnaNgAETkuhXstUP6E6YikWP52cR/iWjUKdlkiEsIU7iEsv6iEf329mVfmbaFFowY8c81ALtAJUxGpBoV7iFqUlsW9H69ma1YeVyV35r7zeumEqYhUm8I9xBzJL+L/PtvIO0u2E9eyEW+PP51hGuRLRGpI4R5Cvlqfyf3/XsO+7AJuHp7AXWN60LCBrjAVkZqrMtzNrDPwBtAOKAUmOeee8ndhdcn+nAIenL6OT1bvpme7GCZdn0z/zs2DXZaI1GLV2XMvBn7vnPvOzGKA5Wb2lXNuvZ9rC3vOOaat2MlDn6wnr6CE349J4paRiRqSV0ROWpXh7pzbDez2Ps82sw1AR0DhfhIyDuZx37S1zN28j8FdWvDY5X3p1iYm2GWJSJio0TF3M4sHBgKLjzNtAjABIC4uzhe1haXSUsdbi7fx6OcbAfjLRadw/RldqFdP3RtFxHeqHe5m1gT4CLjTOXfk2OnOuUnAJIDk5GTnswrDyI4Defzhw9UsSs9iePfW/N9lfenUQhcjiYjvVSvczSwST7C/7Zz72L8lhR/nHG8v3s7/fbYBM+PRy/py1amddTGSiPhNdXrLGDAZ2OCc+6f/SwovGQfzuOejNcxP3c+Z3Vrz2BX96Kjb3YmIn1Vnz30YcD2wxsxWetvuc8595reqwoBzjneX7uDhTzfgnOORS/tyzWnaWxeRwKhOb5n5gBKpBnYdOsofP1rNvJT9DE1sxWOX96NzSx1bF5HA0RWqPuSc44NlGfz1k/WUOMdfLz6F605XTxgRCTyFu4/szc7nno/WMHPjXk5PaMnfr+ivYXlFJGgU7j7wxdo93PvxavIKS3jgwt6MGxKvvXURCSqF+0nIzi/iof+u54PlGfTp2JQnrxqgq0xFJCQo3E/Qki0HuOv9lew6dJSJo7sxcXR3jQkjIiFD4V5DBcUl/OurFF6am0Zcy0Z8cOtQBndpEeyyRER+ROFeA5v2ZHPneyvZsPsI15wWx/3n96JxlD5CEQk9SqZqKC11vLpgC49/sYmmDevzyg3JnN27bbDLEhGpkMK9CplH8rnr/ZUsSM3i7F5tefTyvrRuEhXsskREKqVwr8Q3GzK5+4NV5BeVarAvEalVFO7HkV9UwqOfb+T1hVvp1b4pz1wzkG5tmgS7LBGRalO4HyN1bw4T31nBht1HuHFYPH88tyfRkbpJtYjULgp3L+cc7y/bwYPT19OwQQSTxyXzs146aSoitZPCHTh8tIj7pq3h09W7GdatFf/8xQDaNo0OdlkiIieszof78m0H+M07K8k8ks8fz+3JLSO6alwYEan16my4l5Y6XpiTxj+/2kyH5tF8cOsQBsbpSlMRCQ91MtwP5Bbyu/dWMmfzPi7s34GHL+1D0+jIYJclIuIzdS7cl287wB1TV5CVU8jDl/bh2tPi1HddRMJOnQl35xyT52/h0c830qF5Qz6+bSh9OjYLdlkiIn5RJ8L98NEi7v5gFV+tz+Tnp7Tl71f212EYEQlrYR/uazIOc9vU5ew+lM+fL+jNTcPidRhGRMJelXeXMLNXzWyvma0NREG+4pzjzW+3cfkLCykpcbx/6xB+dWaCgl1E6oTq7Lm/DjwLvOHfUnwnv6iE+6at4ePvdnJWj1j++YsBtGjcINhliYgETJXh7pyba2bxAajFJzIO5nHrW8tZu/MIvzs7iYmju+miJBGpc3x2zN3MJgATAOLi4ny12hpZkLqfO6Z+R3Gp09gwIlKn+eyOzs65Sc65ZOdccmxsrK9WW93fzaS5aVw/eTGtm0Qx/Y4zFewiUqfV+t4yeYXF/OHD1Xyyejfn9W3H36/or/uaikidV6tTcOv+XG55czkpe7O5Z6xn0C/1hhERqUa4m9k7wCigtZllAA845yb7u7CqzNm8j4lTv6NePWPKTacxvHtgDwWJiISy6vSWuSYQhVSXc44pC7fy0CfrSWobw8s3JNO5ZaNglyUiElJq1WGZopJSHpi+jqmLtzOmd1uevGqAjq+LiBxHrUnGQ3mF3Pb2dyxMy+LXoxL5f+f0UP91EZEK1IpwT9uXw69eX8quQ/n848r+XD64U7BLEhEJaSEf7vNS9nHb29/RIKIeU28+neT4lsEuSUQk5IV0uL+3dDv3TVtLt9gmvDIuTE6cOgelJVBS6H0UHfO84Jj2InClnmVcyY+fl5Z6XruS40wv9/NH053ngfuhHk6mjRouy0+n+/Kz9e0KfbiqEK7ND6uT4AvJcHfO8a+vU3j6mxRGJMXy3LUDiQnW+OulpVBwGPIPw9FDnp/5hyH/0A9thblQlAuFeVB01PO86Kj3dbnnJQVQXEDo/k8y+P46gbLnJ9pWyfrgOMv48D34dHW+XF8o1wY+r0+CKuTCvaiklHs/XsOHyzO4cnAnHrmsL5ERPhsl4QfOQc5eOLQNDm2H7N2e1zl7IScTcvd5f+6n0jC2etAgBiIbQoNGEOl9NGgEjVp5Xzf0/oyGiAbeRyRERHl/lm8r97x+FNSrDxYB9ep5f0Z4fuePntfzPj92uvdn+elW76fBqgu/RILv9779fxhS4Z5TUMyv31rOvJT93Hl2d377s+4nd8Wpc57Q3rsB9m2ErLQfwvzQdijO//H89aOhSRto3Aaad4FOydA4Fhq2gOjmEN0MGnp/lr2OilE4ikjICZlw35ddwLhXl7ApM5vHL+/HL07tXLMVOOcJ753LIGMZ7FkD+zZ4Dp2UiW4OLbpAbE/ofg60iPeEePM4aNoeopoqqEUkLIREuGcczOP6yUvYczifyeOSGdWjTdULlZbA7pWQNhO2f+sJ9PxDnmmRjaF9P+hzOcT2gjY9PT+baIgCEakbgh7uaftyuP6VxeQUFPPW+NMY3KWSro7ZeyDlK0+gp8+Gowc87W16Q++LoGOy51BKbE/PMWYRkToqqOG+btdhbpi8BDN4d8IQendo+tOZjuyGDdNh3b9h+yLAQUx76DEWEkdDwkjtkYuIHCNo4b5s6wFufH0pMVH1eWv86XSNbfLDxJIi2PwlfDcFUr/29NVu0xvOug96nu95rmPjIiIVCkq4L0rL4qbXl9KuWTRvjT+djs0beibkHYClr3geOZmePfQz74J+V0FsUjBKFRGplQIe7ovTPcHesUVD3rn5DGJjouDQDlj0LHz3BhTleXqyJN8E3cZARNBPC4iI1DoBTc6yQzEdmkcz9ebTiY3IhS8fgiWTPIde+v4Chk6Etr0DWZaISNgJWLiv2H6Qca8uoV3TaN65aRBtVr0A8/4Fhdkw4FoYeQ80r2HfdhEROa6AhHvq3mxufH0prWOi+OD8CFpNHeO5YjTpXPjZA9pTFxHxMb+H+65DR7l+8hKa1Svgk66fEfPe69C0I1z7PiT93N+/XkSkTvJruB/MLeSGV5fQMT+Vt5u/QNSarXD6LTD6fs+YLCIi4hd+C/eC4hJunrKUIQen85fIN6hX2gp++SnED/PXrxQRES+/hftfPl7Odbsf5tKIBZDwM7hsEjRu7a9fJyIi5VRroHQzO9fMNplZqpndU9X8B7LzuHjtRE+wj74frvtQwS4iEkBV7rmbWQTwHDAGyACWmtl059z6ipZpnL2FQfUaU3rZZOr1u8J31YqISLVUZ8/9NCDVOZfunCsE3gUurmyBSEoouvZDBbuISJBUJ9w7AjvKvc7wtv2ImU0ws2VmtuxA/VgaJY3yUYkiIlJT1Qn34w2/+JObijrnJjnnkp1zya3btD/5ykRE5IRVJ9wzgPLjAnQCdvmnHBER8YXqhPtSoLuZJZhZA+BqYLp/yxIRkZNRZW8Z51yxmd0BfAlEAK8659b5vTIRETlh1bqIyTn3GfCZn2sREREfqdZFTCIiUrso3EVEwpDCXUQkDCncRUTCkDn3k+uRTn6lZtnAJp+v2LdaA/uDXUQ1qE7fUp2+pTp9p4dzzmc3uvDXkL+bnHPJflq3T5jZslCvEVSnr6lO31KdvmNmy3y5Ph2WEREJQwp3EZEw5K9wn+Sn9fpSbagRVKevqU7fUp2+49Ma/XJCVUREgkuHZUREwpDCXUQkDFX3BtmdzWyWmW0ws3Vm9ltve0sz+8rMUrw/W5Rb5l7vDbU3mdnPy7UPNrM13mlPm9nxbgZSYzWt0czGmNlyby3LzWx0uXXN9ta90vto44saT7DOeDM7Wq6WF8utyy+f5QnWeV25GleaWamZDfBOC8bneaX3damZJR+zTEC3zROpMwS3z4rqDLXts6I6A759VlLj381so5mtNrNpZta83DK+2zadc1U+gPbAIO/zGGAz0Bt4HLjH234P8Jj3eW9gFRAFJABpQIR32hJgCJ47PH0OjK1ODX6ocSDQwfu8D7Cz3LpmA8m+qMsHdcYDaytYl18+yxOp85hl+wLpQf48ewE9jv3dwdg2T7DOUNs+K6oz1LbP49YZjO2zkhrPAep72x/DT7lZrT1359xu59x33ufZwAY891G9GJjinW0KcIn3+cXAu865AufcFiAVOM3M2gNNnXOLnKfiN8otc1JqWqNzboVzruyOUuuAaDOL8kUtvqyzIv78LH1Q5zXAO76qpTIV1emc2+CcO95V0gHfNk+kzlDbPiv5PI8r1D7PYwRk+6ykxhnOuWLvbN/iubsd+HjbrPExdzOLx7NXsRho65zbXfZGgLKvMxXdVLuj9/mx7T5VzRrLuxxY4ZwrKNf2mvcr2p99+XXyBOtMMLMVZjbHzIZ72wLyWdawzjJX8dP/PIH+PCsS1G0Tql1neaGwfVYmlLbP6gj49llJjTfh2RMHH2+bNQp3M2sCfATc6Zw7Utmsx2lzlbT7TA1qLJv/FDxfjW4p13ydc64vMNz7uN6XNdawzt1AnHNuIHAXMNXMmhKAz7KGdZbNfzqQ55xbW645lD7PoG2boO0ziHWWzR/w7bOiGs3sT0Ax8HZZ03EWP+Fts9rhbmaR3gLfds597G3O9H5lKPsattfbXtFNtTP44StI+XafqGGNmFknYBpwg3MurazdObfT+zMbmAqc5qsaa1qn9ytalvf5cjzH4ZLw82dZ0zrLuZpj9oqC9HlWJCjb5gnUGWrb53GF4PZZlYBunxXVaGbjgAvw/GEpC2qfbpvV7S1jwGRgg3Pun+UmTQfGeZ+PA/5Trv1qM4syswSgO7DE+zU+28zO8K7zhnLLnJSa1ug9Q/0pcK9zbkG59dQ3s9be55F4/gHK/5UPdJ2xZhbhfd4Vz2eZ7s/P8kTq9C5TD7gSeLdcW7A+z4oEfNs8kTpDcPusaP5Q2z4rWyag22dFNZrZucAfgYucc3nlFvHttumqd9b3TDxfA1YDK72P84BWwDdAivdny3LL/AnPX/FNlDuzCyTj+fDSgGfxXiV7so+a1gjcD+SWm3clnuPHjYHl3vWsA57Ce8Y6SHVe7q1jFfAdcKG/P8uT+DcfBXx7zHqC9XleimePpwDIBL4M1rZ5InWG4PZZUZ2htn1W9u8e0O2zkhpT8RxbL2t70R/bpoYfEBEJQ7pCVUQkDCncRUTCkMJdRCQMKdxFRMKQwl1EJAwp3EVEwpDCXUQkDP1/nQxGHjJCx5cAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "results['Emissions'].plot(title='Emissions (kg/year)',logy=True)\n",
    "results['Radiative forcing'].plot(title='Warming (W m-2)')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}