{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Concentration de CO2 dans l'atmosphère depuis 1958"
]
},
{
"cell_type": "code",
"execution_count": 69,
"metadata": {},
"outputs": [],
"source": [
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"import pandas as pd\n",
"import isoweek"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Les données de la concentration de CO2 dans l'atmosphere sont disponibles sur le [site Web de l'institut Scripps](https://scrippsco2.ucsd.edu/data/atmospheric_co2/primary_mlo_co2_record.html). Nous les récupérons sous forme d'un fichier en format CSV dont chaque ligne correspond aux observations sur une semaine specifique. Nous téléchargeons ici le jeu de données complet, qui commence en 1958 et se termine avec une semaine récente."
]
},
{
"cell_type": "code",
"execution_count": 70,
"metadata": {},
"outputs": [],
"source": [
"data_url = \"https://scrippsco2.ucsd.edu/assets/data/atmospheric/stations/in_situ_co2/weekly/weekly_in_situ_co2_mlo.csv\""
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Nous avons recupere les observations hebdomadaires"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Date du telechargement : **16/04/2024**"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Nous allons charger les donnees a partir d'un fichier CSV locale"
]
},
{
"cell_type": "code",
"execution_count": 71,
"metadata": {},
"outputs": [],
"source": [
"data_file = \"weekly_co2_measures.csv\"\n",
"\n",
"import os\n",
"import urllib.request\n",
"if not os.path.exists(data_file):\n",
" urllib.request.urlretrieve(data_url, data_file)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Les 44 premières lignes du fichier CSV sont des commentaires, que nous ignorons en précisant `skiprows=44`"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"De plus, on renomme les colonnes du dataset pour avoir quelque chose de plus comprehensible (**Date** et **Concentration**)"
]
},
{
"cell_type": "code",
"execution_count": 72,
"metadata": {
"scrolled": false
},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Date
\n",
"
Concentration
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
1958-03-29
\n",
"
316.19
\n",
"
\n",
"
\n",
"
1
\n",
"
1958-04-05
\n",
"
317.31
\n",
"
\n",
"
\n",
"
2
\n",
"
1958-04-12
\n",
"
317.69
\n",
"
\n",
"
\n",
"
3
\n",
"
1958-04-19
\n",
"
317.58
\n",
"
\n",
"
\n",
"
4
\n",
"
1958-04-26
\n",
"
316.48
\n",
"
\n",
"
\n",
"
5
\n",
"
1958-05-03
\n",
"
316.95
\n",
"
\n",
"
\n",
"
6
\n",
"
1958-05-17
\n",
"
317.56
\n",
"
\n",
"
\n",
"
7
\n",
"
1958-05-24
\n",
"
317.99
\n",
"
\n",
"
\n",
"
8
\n",
"
1958-07-05
\n",
"
315.85
\n",
"
\n",
"
\n",
"
9
\n",
"
1958-07-12
\n",
"
315.85
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Date Concentration\n",
"0 1958-03-29 316.19\n",
"1 1958-04-05 317.31\n",
"2 1958-04-12 317.69\n",
"3 1958-04-19 317.58\n",
"4 1958-04-26 316.48\n",
"5 1958-05-03 316.95\n",
"6 1958-05-17 317.56\n",
"7 1958-05-24 317.99\n",
"8 1958-07-05 315.85\n",
"9 1958-07-12 315.85"
]
},
"execution_count": 72,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# nous travaillons avec la version locale des donnees\n",
"raw_data = pd.read_csv(data_file, skiprows=44, sep=r'\\s+', engine='python', parse_dates=[0], names = ['Date', 'Concentration']) \n",
"\n",
"# Appercu des donnees (10 premieres lignes)\n",
"raw_data.head(10)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" Y a-t-il des points manquants dans ce jeux de données ? Non"
]
},
{
"cell_type": "code",
"execution_count": 73,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(15, 6))\n",
"plt.plot(raw_data['Date'], raw_data['Concentration'])\n",
"plt.title('Concentration de CO2 au fil du Temps')\n",
"plt.xlabel('Annee')\n",
"plt.ylabel('Concentration de CO2')\n",
"plt.grid(True)\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Réalisation du graphique pour montrer l'oscillation lente (200 dernières lignes du tableau)"
]
},
{
"cell_type": "code",
"execution_count": 75,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA4wAAAGDCAYAAABgA55eAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xd8nGeV9//PUe+SJdmSVWy5xL0lcXpInAohjUBgqRtq4FlYYHcpmwfYH7CUpbNs4GEpW0iAJKRASHeKU5w4sR333q3eLI16nev3x8zII2nUbI1mJH3fr5df0dz3PfecmXsUzZnrus4x5xwiIiIiIiIiA8VEOgARERERERGJTkoYRUREREREJCQljCIiIiIiIhKSEkYREREREREJSQmjiIiIiIiIhKSEUUREREREREJSwigiIhFlZnPMrMXMYifgsTaY2cfD/TjRyMwWm9k2M2s2s8+a2S/N7Gv+fevMrGwM55q2r6OIyHSjhFFEJILM7P1mtsWfMFWa2ZNmdnmk4wrFzErMzJlZ3Fme57iZXRu47Zw76ZxLc871nn2U4WNmCWb2dTM7ZGat/ufxX2ZWEnTMTWb2hn9/vZn93syKgvbfaGavmFmjmVWZ2a/NLH2CnsKXgA3OuXTn3M+cc59yzv3rBD02AP73eeCf18zag25/YCJjERGR0VHCKCISIWb2j8BPge8AecAc4BfArZGM62ycbTIZ5R4EbgHeD2QCq4GtwDUAZnY78Afg34FcYDnQCbxiZjP858gEvgUUAEuBIuAHExT/XGDPBD1WSP4vBtKcc2nASeDmoG2/j2RsIiISmhJGEZEIMLNM4JvAp51zDzvnWp1z3c65vzrnvug/JtHMfmpmFf5/PzWzRP++dWZWZmb/ZGY1/tHJjwSdP9nMfmRmJ8zM4x/VSvbvu9jMXvWPcu0ws3VB99tgZv9qZhv9UxefMbNc/+6X/P9t9I8IXWJmH/Yf+xMzOwV83cwWmNnz/hG2Ov8oW5b//PfgS4z/6j/HlwaOXJpZgZk9amanzOywmX0iKL6vm9kDZvY7f3x7zGztMK/zdWa23/8a3A3YgP0fNbN9ZtZgZk+b2dwhznMtcB1wq3Nus3Ouxznncc793Dn3WzMz4EfAt5xzv3fOtTvnqoCPAy3APwA45/7gnHvKOdfmnGsAfg1cNkz8/2xmR/zPda+Z3Tbgtbg36PaQI8Bm9jxwFXC3/3VfZGb/Y2bfGuqxR/s6jiWOUTxOrJl9zcyOhnjvLDGzHjP7mJmV+99fH/W/D3f7388/DjrXp/zvw/80syb/63dF0P5PmG+UuNn/eO8ea7wiItOBEkYRkci4BEgCHhnmmK8AFwNr8I1mXQh8NWh/Pr4Rq0LgY8DPg0ayfgicD1wKZOObjug1s0LgcXyjXNnAF4CHzGxm0HnfD3wEmAUk+I8BCHzYzvKPCL3mv30RcNR//LfxJRPf5fQoWjHwdQDn3IfoP7L0/RDP+49Amf/+twPfMbNrgvbfAtwHZAGPAneHevH8ie5D/tcsFzhCUHJmZu8A/i/wTmAm8LL/sUO5FnjDOVc6xP7F+BLhPwVvdM55/TFcN8T9rmD4Ub8jwFvwXedvAPea2exhjg/JOXc1vuf3Gf/rfnC09x3pdRxnXwSuBy7HN/raDfwkaH8ssAqYj+89+h/43p9X+rd/xMwuCjr+CmAHkAP8G/BnM8vw/578ALjGOZeO7zXeHabnJCIyqSlhFBGJjBygzjnXM8wxHwC+6Zyrcc7V4ksYPhS0v9u/v9s59wS+kazFZhYDfBT4nHOu3DnX65x71TnXCXwQeMI594RzzuucWw9sAd4edN7/ds4ddM61Aw/gS1iHU+Gc+w//qFu7c+6wc269c67TH/eP8X2gH5GZFeNLFr7snOtwzm0HfjPgeb/ij78XuAdfMh3K24G9zrkHnXPd+Kb/VgXt/yTwXefcPv91+A6wZohRxhygcpjQA6OwoY6pDNrfx8yuA+4A/mWokzrn/uScq/Bfq/uBQ/i+OJhII72O4+mTwD/7n3MHvvf83/hHcAO+6X9vPeq//TvnXL1z7iTwKnBu0LGlzrlf+H9Hfofvi4i3Bu1fYWZJ/t+TfWF6TiIik5oSRhGRyKgHckeYtlcAnAi6fcK/re8cAxLONiANX3KShG8kaKC5wLv90/cazawRX4IWPGoVnAwEzjmcfqNuZjbLzO7zTxtsAu4lRMI0hALglHOuOWjbCXyjqEPFlzTE61gQHJtzzg2IdS7w70Gvwyl8o6OFDFZP/9dooDr/f0MdMztoP+CbFoxvvePtw432mdnfmtn2oBhXMPrXcryM9DqOC39SWAw8EfR8t+H7rJLjP6zXOVcfdLd2oHrA7eD368DKryeAAv904A8AnwWq/FOgF47fsxERmTqUMIqIRMZrQAfwjmGOqcCX1ATM8W8bSZ3/3AtC7CsF7nHOZQX9S3XO/dsozutGuf27/m2rnHMZ+EY1bZjjg1UA2da/cugcoHwU8Q1UiS8BAfolJAGlwCcHvBbJzrlXQ5zrWeBCC6p4OsABfMlJv3Vw/tHedwHPBW07F99U2o86555jCP6Rzl8DnwFynHNZ+KZNBl7LViAl6C75Q53rLI30Oo5LHP5EtBy4esA1SXLO1Y10/yEMvF59v0POucedc9fgS4hPAv/vDB9DRGRKU8IoIhIBzjkPvqmIPzezd5hZipnFm9kNZhZY1/dH4KtmNtO/juxf8I3WjXRuL/BfwI/NV0Am1l8YJNF//5vN7K3+7UnmK6AzVCIUrBbw4ls/Npx0fNNjG/1rJr84YH/1UOfwrxF8FfiuP7ZV+NZnnkkFzceB5Wb2Tv8I5Gfpn8z8ErjLzJaDrxDRUIVPnHPPAuuBR8zsfDOLM7N0f2GVj/qTnS/gu17vN1/RoXx802kz8K/DM7MVwFPA3zvn/jpC/Kn4kuta/30/gm+EMWA7cIX5+lhmAneN+pUZm5Fex/GM45fAv/mnJgdGq28+i/MV+69RnJl9EF/C+IyZFZqvxUkKvkq2LUBUt3UREYkUJYwiIhHinPsx8I/4ionU4hvx+gzwZ/8h38K3vnAnsAt4079tNL7gv89mfFMtvwfE+BOyW/EVewk85hcZxd8D51wbvqI2G/1TBi8e4tBvAOcBHnzJxsMD9n8XX2LVaGZfGHhn4H1ACb6RoEeA/8+/1nJM/KNS78ZX7KQeOAfYGLT/EXyvy33+qbO7gRuGOeXtwBPA/f7nthtYi2/0Ef8aww/hq4haB+wFkoHLgqZR/hO+Aju/tdP9B0MWvXHO7cVXefU1fEn2ygHxr/fHshNfe4/HRvO6jNUoXsfxjOP7+F7P582sGd+XB+edxflewrem8RS+IlK3+b+sicWX2Fbhe04XAH9/Fo8jIjJlme9LUREREZGpw8w+hW+N6LWRjkVEZDLTCKOIiIiIiIiEpIRRREREREREQtKUVBEREREREQlJI4wiIiIiIiISkhJGERERERERCSku0gGcjdzcXFdSUhLpMEaltbWV1NTUSIch6FpEE12L6KFrER10HaKHrkV00HWIHroW0WPgtdi6dWudc25muB5vUieMJSUlbNmyJdJhjMqGDRtYt25dpMMQdC2iia5F9NC1iA66DtFD1yI66DpED12L6DHwWpjZiXA+nqakioiIiIiISEhKGEVERERERCQkJYwiIiIiIiISkhJGERERERERCUkJo4iIiIiIiISkhFFERERERERCUsIoIiIiIiIiISlhFBERERERkZDCnjCaWayZbTOzx/y3f2Bm+81sp5k9YmZZQceuMrPXzGyPme0ys6RwxyciIiIiIiKhTcQI4+eAfUG31wMrnHOrgIPAXQBmFgfcC3zKObccWAd0T0B8IiIiIiIiEkJYE0YzKwJuBH4T2Oace8Y51+O/uQko8v98PbDTObfDf1y9c643nPGJiIiIiIjI0Mw5F76Tmz0IfBdIB77gnLtpwP6/Avc75+41s88D5wOzgJnAfc6574c4553AnQB5eXnn33fffWGLfzy1tLSQlpYW6TAEXYtoomsRPXQtooOuQ/TQtYgOug7RY7pdC+ccJ5u9zEmPwcwiHU4/A6/FVVddtdU5tzZcjxcXrhOb2U1AjXNuq5mtC7H/K0AP8PugWC4HLgDagOfMbKtz7rng+znnfgX8CmDt2rVu3bpBp45KGzZsYLLEOtXpWkQPXYvooWsRHXQdooeuRXTQdYge0+1a/PDpA9z96mHuu/NiLp6fE+lw+pnoaxHOKamXAbeY2XHgPuBqM7sXwMzuAG4CPuBOD3GWAS865+qcc23AE8B5YYxPRERERESkn79sL+fuFw4DsPVEQ4SjibywJYzOubucc0XOuRLgvcDzzrkPmtnbgC8Dt/gTw4CngVVmluIvgHMlsDdc8YmIiIiIiATbdrKBLz64kwvnZTMnO4WdZY2RDiniItGH8W58axrXm9l2M/slgHOuAfgxsBnYDrzpnHs8AvGJiIiIiMg0U9HYzp33bCUvI5FffvB8zpuTxY5ST6TDiriwrWEM5pzbAGzw/7xwmOPuxddaQ0REREREZEL0eh2fvGcr7V29/P7jF5GdmsCqoiz+vL2CmqYOZmVM3/bwkRhhFBERERERiRrbTjawq9zDv9y8jEV56QCsLs4EYEfZ9B5lVMIoIiIiIiLT2vq91cTHGjesyO/btmx2JrExxo7S6b2OUQmjiIiIiIhMa+v3VnPx/BzSk+L7tiUnxLI4L50d07zwjRJGERERERGZtg7XtHC0rpXrl+UN2re6OJOdZR5OdwKcfpQwioiIiIjItLV+bzUA14ZIGFcVZeFp7+ZEfdugfdOFEkYREREREZm2ntlbxcrCTGZnJg/at7ooC2BaT0tVwigiIiIiItNSTXMH20sbuS7E6CLAorw0kuJj2DmNK6UqYRQRERERkWnpuX01OMeQCWNcbAzLCzKndaVUJYwiIiIiIjItrd9bTdGMZJbkpw95zKqiTHZXeOjp9U5gZNFDCaOIiIiIiEw7rZ09vHK4juuW5WFmQx63pjiLjm4vh2paJjC66KGEUUREREREpp2XD9XS1ePl+mX5wx63KlD4ZppOS1XCKCIiIiIiU07pqTa++8Q+uoeYSvrM3mqyUuK5oGTGsOcpyUkhIymOHdO08I0SRhERERERmXIe21nJf750lNeO1A/a193r5fn9NVy9eBZxscOnRGbGqqIsdk7T1hpKGEVEREREZMopa2gD4Kk9VYP2bTpaT2NbN29bMfx01IDVxZnsr2qmo7t3XGOcDJQwioiIiIjIlFPW0A7AM3uq6fW6fvue3F1FSkIsVyyaOapz3biygO+9axXOjXzsVKOEUUREREREppyyhjZSEmKpa+lk28mGvu29Xscze6q4askskuJjR3WuZQUZ3H5+EckJozt+KlHCKCIiIiIiU4pzjrKGdm5dU0BCbAxP7T49LXXL8VPUtXRxwyino053ShhFRERERGRKqW3ppLPHy5L8DC4/J5en9lTh/PNJn9xdRWJcDFctnhXhKCcHJYwiIiIiIjKlBNYvFs1I5q3L8yhraGdPRRNer+PpPVVcsWgmqYlxEY5yclDCKCIiIiIiU8rphDGFa5fmEWPwzJ4qdpQ1Uunp0HTUMVBaLSIiIiIik1JjWxepiXHED+ilGGipUTQjmdTEOC6cl81Te6ro7PESH2tcszQvEuFOShphFBERERGRSaeju5erfriB/7fhyKB9pafayU5N6Jt2+tbl+RysbuH+LaVcuiCXzOT4iQ530lLCKCIiIiIik86mo/U0tHWz5UTDoH1lDW0UzUjuu/3W5b4pqI1t3ZqOOkZKGEVEREREZNJ5YX8NAPsqmwbtK29o75cwFmQls6ookxiD65ZpOupYaA2jiIiIiEQF5xx3P38YB6wqyqSly0U6JIlSzjme21+DGdQ2d1LX0kluWiIAXq+jrLF9UGL4xbcu5mB1Czn+42R0lDCKiIiISFSob+3iR+sP9tv2sz0v8ts7LqA4OyVCUUk0OlzTQllDOzetms1jOyvZV9nEW86ZCfh6MHb1ePuNMAK85ZyZfcfI6GlKqoiIiIhEhYpGXyuEH79nNX/4+EW8e1E8J0+18dNnD0U4Mok2z/mno/7duoVA/2mppyuk6kuG8aCEUURERESiQkVjBwCL8tK5dGEuN85P4P0XzuXP28s5Wd8W4egkmjy/v4alszNYVpBBfkYS+yqb+/YFejAWZycPdXcZAyWMIiIiIhIVAiOMBVmnP+h/8sr5xMYYv9hwOFJhSZTxtHWz9UQDVy/xTS9dOjt9wAij731UmKURxvGghFFEREREokKlp53EuBhmpJzukZeXkcR7Lyjmwa1lfVMNZXp78VAtvV7H1Ut8RW2Wzs7gcE0LnT29AJSeaiM3LYHkhNhIhjllhD1hNLNYM9tmZo/5b//AzPab2U4ze8TMsgYcP8fMWszsC+GOTURERESiR0VjB4VZyZhZv+2funIBZvDLFwc3aJfp5/l91WSnJrCm2JdGLJ2dQY/Xcai6BfCNMBZq/eK4mYgRxs8B+4JurwdWOOdWAQeBuwYc/xPgyQmIS0RERESiSIWnndlZSYO2F2Ql8+61xTywuYxKT3sEIpNo0et1bDhYy7pFM4mN8X2xsHR2BnC68E1ZQxvFM7R+cbyENWE0syLgRuA3gW3OuWeccz3+m5uAoqDj3wEcBfaEMy4RERERiT6VjR3Mzgz9Qf//XLkAr3P854tHJzgqiSbbTjbQ2NbNVUtm9W2bl5tKUnwM+yqb8Xod5Y3tqpA6jsI9wvhT4EuAd4j9H8U/mmhmqcCXgW+EOSYRERERiTLdvV6qmzv6FbwJVpydwm3nFvKHN07S3tU7wdFJtHh+fw2xMcYVi073U4yNMRbnZ7Cvsonq5g66e92gHoxy5uLCdWIzuwmocc5tNbN1IfZ/BegBfu/f9A3gJ865loHz1gfc707gToC8vDw2bNgwzpGHR0tLy6SJdarTtYgeuhbRQ9ciOug6RA9dizPX1On48dYO3jYvnotnj/6jZl27F+egueoEGzZUAIOvw6zeHrp6vNz35AbmZaqgyUSKlt+J9dvbmZdhbHt9Y7/tWa6TLaU9PPrcqwCcKj3Eho5jkQgx7Cb6WoQtYQQuA24xs7cDSUCGmd3rnPugmd0B3ARc45xz/uMvAm43s+8DWYDXzDqcc3cHn9Q59yvgVwBr165169atC+NTGD8bNmxgssQ61elaRA9di+ihaxEddB2ih67FmXtsZwXHm7bxq52dLF26lFvXFI7qfpuPn4IXX2PdRWu40j96NPA6zKlt4efbXyS1cBHr1haHI3wZQrT8Ttz12nNcsiCHdevW9Nt+MvE4L/5lD73Zc4ED3LjuYhbMTItIjOE20dcibAmjc+4u/AVt/COMX/Ani2/DN/X0SudcW9Dxbwn8bGZfB1oGJosiIiIiEt12lXtIiI3hvLlZ/MP92/E6x23nFo14v0APxsIQRW8C5ub41qrtD2rSLtNHT6+X6qYOCkKscw0Uvlm/txqAwiGmNsvYhXOEcSh3A4nAev/U003OuU9FIA4RERERGWe7yz0szk/nvz98IR/738384wM7aGjtpiArieqmTqqaOlg4M413nd8/iaxo7AAYsugN+NaqLcpL50B105DHyNRV09yJ1xGyku6S/HQAtp1sZGZ6IknxmrI8XiYkYXTObQA2+H9eOIrjvx7eiERERERkvDnn2F3exI2rZpOcEMtv77iAj/9uM998bG+/45LiY3jHuYV9bRHAN8KYmRxPauLwH08X56XzwoGasMQv0a3S4/tSIdQIY3pSPMXZyZSealdLjXEWiRFGEREREZmCSk+142nvZmVhJgDJCbH814cvYOuJBjKS4snLSOKF/TV86aGdHKtrZeGs02vMKj3tzM4cejpqwJLZGfxpaxm1zZ3MTE8M23OR6BPowRlqhBFgaX4GpafUUmO8hbuthoiIiIhME7vKPQB9CSNAYlwsly7IZUVhJjPTE1nh37enwtPvvhWNQ7fUCBaYenigSusYp5vKEaYtB9YxqqXG+FLCKCIiIiLjIlDwZlFe+pDHLJyVRkJsDHsr+69DrPC0UzBMwZuAxf6EcX+V1jFONxWedlITYslICj1JcllBIGHUCON4UsIoIiIiIuNiV3kji/PTSYgb+iNmQlwM5+SlsbfidMLX1tVDY1v3sAVvAnLTEslNS2S/RhinncrGDvIzkxiqZ/tF87K5eH42ly7ImeDIpjYljCIiIiJy1gIFb1YETUcdyvKCDPZWNBFoxx2okDqaEUbwTUvVlNTpp9LTPuy05ayUBO678xJKclMnMKqpTwmjiIiIiJy1gQVvhrNsdgb1rV1UN3UCp4uZhKp+GcqS/HQOVjfT63VnHrBMOpWejlEVRpLxpYRRRERERM5aqII3Q1nuP2Zvpe8+lX0jjKNLGBfnp9PZ4+V4feuZhCqTUFePl9qWzlFNW5bxpYRRRERERM7arnIP8bHGovy0EY8NVLPcU+5bx1je2I4Z5GWMdkqq7/6aljp9VDd14Nzopy3L+FHCKCIiIiJnbXe5h8X56STGxY54bFpiHCU5KX2VUis97cxMSxy2WE6wc/LSiDHYX6lKqdNFpWf4lhoSPkoYRUREROSsOOfYVe4Z1XTUgOUFmezxV0qtaOxg9iinowIkxcdSkpuqSqnTSN86V40wTjgljCIiIiJyVsoafAVvRlMhNWBZQQYnT7XR1NHt68E4xmImS/MzOFCthHG6CFTSzdcI44RTwigiIiIiZ2UsBW8CAk3W91U0UdnYMeqCNwGL89M5Ud9Ga2fPmO4n0aWls4ev/nkXR2tbhj2uytNOelIcaYlxExSZBChhFBEREZGzEih4szg/fdT3We4vfLPxSD3t3b1jbpcQeKyDGmWc1L712F7u3XSSnz13aNjjKjwdo267IuNLCaOIiIiIjMlLB2v56bMH+euOCvZWNLH9ZCOL8kZX8CZgVkYSuWmJPLevGhh9S42ApaqUOuk9v7+a+zaXkpOawBO7q2ho7Rry2EpPO7O1fjEiNKYrIiIiImPy9Uf3cLSufw/E915QPObzLC/I4MWDtcDYE8aiGcmkJMSq8M0k1dDaxZcf2sWS/HS+965V3PrzjTy8rZyPXT4v5PGVjR2sLMya4CgFlDCKiIiIyBjUt3RytK6Vf7h2Edcvz+NIbQsn6tu4ceXsMZ9rWXDCOMYpqTExvimw+9RaY1L66l9209jWxf9+5EKWFWSwpjiLP75xko9eVoKZ9Tu2o7uX+tauMb9HZHxoSqqIiIiIjNqWEw0AXH5ODktnZ3DTqgI+fdVCSnJTx3yu5f7CN/GxRm5a4pjvv7ooix1ljXT29I75vhI5j+6o4PGdlXz+2kV9xY/ef+EcDte09L2/glUFejCOcRRaxocSRhEREREZtS3HT5EQFzOmFhpDWeYvfJOfmURMjI1w9GCXL8ylo9vLmycaQ+6vaGzHOXdWMcr4+7cn9rG6OItPXjG/b9tNq2eTlhjHH14/Oej4ykDCqBHGiFDCKCIiIiKjtvl4A6uLMsdU4GYoJTmppCbEMvsMq19eND+b2Bhj4+G6QftO1Ldyxfdf4N5NJ842TBlH9S2dVHg6uHnVbOJiT6ciKQlxvOPcAh7fVUljW//iN5WedkAJY6QoYRQRERGRUWnv6mV3uYe1Jdnjcr6YGOP9F805o/WPAOlJ8awpzuKVEAnjE7uq6PE67tl0QqOMUSRQ1TZUC5b3XziXrh4vD79Z3m/76RFGTUmNBCWMIiIiIjIqO8oa6fE6LiiZMW7n/MqNy7jj0pIzvv9lC3PZWdaIp7273/andlcSH2scrG5ha4h1cRIZ+4dJGJcVZLDaX/wmOMmvaGxnRko8yQlnP6otY6eEUURERERGZcvxUwCcP2d8RhjHw+ULc/E6eO1Ifd+2soY2dpR5+NSVC3zr4t4YvC5OIuNAVTPZqQnMHKLI0fsvLOZQTUu/UeNKT4dGFyNICaOIiIiIjMrm4w0szksnMyU+0qH0WVOcRUpCbL91jE/trgLgnecV+dbF7azE09Y91CnkLDV3dPPNv+4dNMobyv7qZhbnpQ9qnRFw65pC8jOS+Omzh/pGGSsa2ynI0vrFSFHCKCIiIiIj6vU63jzRwPnjOB11PCTExXDRvOxBCeOS/HTm5abyvgvn0Nnj5ZFtZRGMcmp75VAd/7XxGPdvHn4k1+t1HKpuDjkdNSApPpZPX7WArScaePmQ75pqhDGylDCKiIiIyIgOVDXT3NkzrusXx8tlC3M5WtdKeWM7NU0dbD3ZwA0rfIV0lhdksrookz++UariN2FyrL4VgAe3lg37Gpc2tNHW1cuSYRJGgPdcUExBZhI/efYgbV09eNq7yVeF1IhRwigiIiIiI9pywrd+ce3c6Fm/GHD5ObkAbDxcx9N7qnAObliZ37f/fRfO4UB1M2+eDN2vUc7Oibo2AA5Wt7CzzDPkccMVvAmWGBfL3121kG0nG7l/cymApqRGkBJGERERERnRluMN5GckUTQj+qYGLs5LJzctkY2H63hydxXzZ6Zyzqy0vv03ry4gLTGOP6r4TVgcq29lSX46iXExPLh16Km/gZYai/KGTxgB3rO2mMKsZH70zEFALTUiSQmjiIiIiIxoy/FTnF8yY8hiJZFkZly+MIcNB2p5/dgpbliR3y/O1MQ4bl1TwGM7K2jqUPGb8XaivpUVhZm8dXk+j+6ooKO7N+RxB6qamZOdQmpi3IjnTIiL4TNXL6SlsweAAiWMEaOEUURERESGVd7YToWngwvmRt/6xYDLFubiae+m1+v61i8Gu3l1AR3dXjYfOxWB6Kautq4eqps6mZebyu3nF+Fp7+a5fTUhj91f1TTidNRgt59f1DeinZcZug2HhJ8SRhEREREZVqD/4tqS6Fu/GHDZQt86xuLsZJYXZAzav6ook9gYY3up1jGOp+P+9Ytzc1K4bGEu+RlJPLi1dNBxHd29HK9vG7HgTbD42Bi+fdtKPnb5PBLjYsctZhmbkceDRURERGRa23K8gdSE2DF92J9oBVnJXLcsj4vn54ScNpuSEMeivHQljOPshL9CaklOKrExxjtBlsAlAAAgAElEQVTPK+SXLx6hpqmDWRmnC9Ucrmmh1+vGNMIIcOWimVy5aOa4xixjE/YRRjOLNbNtZvaY//YPzGy/me00s0fMLMu//Toz22pmu/z/vTrcsYmIiIjIyLaVNrC6OIu42OienPbrv13Lxy6fN+T+NcVZbC9txOtVe43xEmipUZKbCvimkXodPLKtvN9xgYI30fylg4Q2Eb/1nwP2Bd1eD6xwzq0CDgJ3+bfXATc751YCdwD3TEBsIiIiIjKM9q5e9lU2c96c6F2/OFrnzsmiuaOHo3WtkQ5lyjhe10puWiJp/kI282emcf7cGfxpQE/GA9XNJMTFUJKTGqlQ5QyFNWE0syLgRuA3gW3OuWeccz3+m5uAIv/2bc65Cv/2PUCSmWl1q4iIiEgE7SxrpNfrOHdOVqRDOWvnFvueg6aljp/j9W3My03pt+0DF83hcE0LD2w5vZZxf1UzC2emRf0otQxmwZn/uJ/c7EHgu0A68AXn3E0D9v8VuN85d++A7bcDn3LOXRvinHcCdwLk5eWdf99994Ur/HHV0tJCWlrayAdK2OlaRA9di+ihaxEddB2ih67FaU8c7eKBg9387OoUMhImtqXGeF8Hr3N8+rk2Lp4dxx3LNS4xFkNdi8+/0MaK3Fg+vvL06+l1ju+90UFps5fvvCWZrMQYPv9CG8tyYrlzlV73szXwWlx11VVbnXNrw/V4YSt6Y2Y3ATXOua1mti7E/q8APcDvB2xfDnwPuD7UeZ1zvwJ+BbB27Vq3bt2gU0elDRs2MFlinep0LaKHrkX00LWIDroO0UPX4rQ/lm6hJKeZW66/asIfOxzX4bwjm6hp62bdureM63mnulDXoq2rh8annuaSFQtYt25hv31zV7Rww7+/zPq6LL592woan1rPlasXsu7KBRMY9dQ00f9/CueY8GXALWZ2HLgPuNrM7gUwszuAm4APuKAhTv8U1keAv3XOHQljbCIiIiIyAuccb55s5NwpsH4xYE1xFvurmmnvCt1cXkYvuKXGQAtmpvHZqxfy+K5Kfv7CYYAxV0iV6BC2hNE5d5dzrsg5VwK8F3jeOfdBM3sb8GXgFudcW+B4f7XUx4G7nHMbwxWXiIiIiIxOeWM7tc2dU2L9YsCa4hn0eh27KzyRDmXSOx7UUiOUO69YwOK8dH798jEAluQP7o8p0S8Sq07vxremcb2ZbTezX/q3fwZYCHzNv327mc2KQHwiIiIiAmw76SsOMxUqpAasCRS+OanCN2fr+ICWGgMlxMXwb+9aiRlkJseTl6H1i5NR2NYwBnPObQA2+H9eOMQx3wK+NRHxiIiIiMjI3jzZQFJ8zJSaSjgzPZHCrGRVSh0HA1tqhHLunBl88a2Lae7owWxiiybJ+JiQhFFEREREJp9tJxtZVZhF/BRrhbBmTpZGGMdBqJYaofzdupDjRTJJTK3ffhEREREZF509veytaOLcuVNn/WLAucVZlDe2U9PcEelQJrXjda3MHWL9okwdShhFREREZJDd5U109Xo5t3jqrF8M0DrGs9fa2UNNcyfzhli/KFOHEkYRERERGWTbyQYAzptCFVIDVhRmEhdjWsd4Fk7UD91SQ6YWJYwiIiIiMsi20kYKs5KZlZEU6VDGXVJ8LEtnZyhhPAsjtdSQqUMJo4iIiIgMsu1Ew5TqvzjQmuIsdpQ20t3rjXQok9JILTVk6lDCKCIiEkYv7K/htl9spKsn9IfStq4eevSBVaJMlaeDCk/HlOq/ONClC3Jo7erVKOMZGk1LDZkalDCKiIiE0atH6th2spFd5Z5B+5xz3Pwfr3D9T1/qWy8mEg1eO1oHMKVHGC9dmEtsjPHSwdpIhzIpHa8bXUsNmfyUMIqIiIRRhcdXtn/z8VOD9p2ob+NIbStlp9p51/97lQcPdtHZ0zvRIYoM8vtNJynJSWF10dRNGDOT41lTnKWE8Qwdr1dLjelCCaOIiEgYVTa2A7D52OCEcdPRegAe+NQlvPv8Yh472s2td2+k9FTbhMYoEmxvRRNbTjTwwYvnEhNjkQ4nrK44ZyY7yz2cau2KdCgR55xj8/FT/Pqlo3z6D29y+fee50O/fT3ksc0d3WqpMY0oYRQREQmjKv8I45YTDXi9rt++14+dIjctkdVFmXzv9lV8/rxEyhva+eKDO3DOhTqdSNjds+kESfExvPv84kiHEnZXLMrFOXj5kEYZ/7y9nHf/8jW+/cQ+tp9sJDk+lpcP1eFp7x507J6KJgCWFWRMdJgSAcMmjGaWYWYLQmxfFb6QREREpoZer6O6uZPZmUl42rs5VNPSt885x6aj9Vw0Pxsz3yjOmllxfOXGpWw6eooHtpRGKmyZIu7ffJKrfrhhTNOcmzq6+fO2cm5ZXUBmSnwYo4sOq4qyyEqJ56WDdZEOJeKO1bYSY/DGV65h4z9fzVdvWgbAnhDrr3eV+batLMyc0BglMoZMGM3sPcB+4CEz22NmFwTt/p9wByYiIjLZ1TR30Ot13LK6AOi/jvHkqTYqPR1cPD+n333+5oJiLpqXzbcf30dNc8eExitTR3NHN9976gDH6lrZX9k86vs9tLWM9u5e/vaSkvAFF0ViY4zLFuby8qHaaT+qX+npYGZ6IrPSfX03V/mTwZ2hEsZyDwWZSeSmJU5ojBIZw40w/l/gfOfcGuAjwD1m9k7/vqk9oV1ERGQcVPqno144L5tZ6Yn9EsbXj/p+vmR+dr/7mBnffedKOnq8fOPRvRMXrEwpv375WN+6vFAf+ENxznHPphOcOyeLFdNo5OjKc2ZS09zJ/qrRJ9ZTUVVTB/mZyX23Z6QmUJyd3DeaGGxXuYeVRdPnPTLdDZcwxjrnKgGcc28AVwFfMbPPAtP7KxgREZFRqGz0JYyzM5O5oCSbLcdPt87YdLSe3LQEFsxMG3S/+TPT+Nw15/D4rkrW762esHhlaqhr6eQ3Lx/l7SvzyU5NYOco+wy+eqSeo7WtfOjiuWGOMLq8ZVEuwLSvllrl6SA/o/+I4arCLHaW93//NHV0c6yuVdNRp5HhEsbm4PWL/uTxKuBWYHm4AxMREZnsKj2+CqkFWUlcUDKD8sZ2yhvbT69fnJfTt35xoDuvmM+S/HS+9ufdtHb2TGTYMsnd/fxhOnu8/NP1i1lVlMnOECNEofzuteNkpybw9pWzwxtglJmdmcyivDReGlD4pvRUGx3d06fNTZWng9lBI4wAK4syKT3VTkNQFdnd/hHrlVO45Yr0N1zC+H8G7nfONQFvAz4azqBERESmgkpPB8nxsWQmx7O2xDf1dMvxU5Q1tFPh6eDiAdNRg8XHxvCv71hBVVMHD79ZNlEhyyRXeqqN379+gvesLWLBzDRWFWZyqKaZtq7hv3Sob+lk/d5q3rO2mKT42AmKNnpccc5MNh9roK2rh55eLz9+5gBX/OAFvvPEvkiHNiFaOnto7uwhPzOp3/ZQ6xhV8Gb6GTJhdM7tcM4dMrN5ZnaTmd1oZvOdc93Oud9PZJAiIiKTUaWnndmZSZgZS2dnkJYYxxvHTvGav//iRQMK3gy0du4MVhZmcs+mE9O+IIeMzk+ePUiMGZ+95hzAVwXU6063QRjKtpONeB1cs3TWRIQZda5YNJOuXi+PbCvnfb/exM+eP0xOaiKPvFk+YrI9FQTa/+Rn9E8Yl/uTwl1lp6el7ir3UJiVTHZqwsQFKBE1XJXUDDN7AHgO34jix4FnzexPZqamKyIiIiOoaOxgdpbvA1hsjHHe3BlsOd7ApqP1ZKcmcM6swesXg5kZH7pkLgerW3j92KlhjxU5UtvCI9vK+fClJX1TC1cV+z7w7xhhHeOOskZiY4zl07Sv3oXzskmMi+Erj+xmX2Uz//7eNfziA+fR3NnDE7uqIh1e2PUljANGGDOT45mXm9pvWvOucg+rVPBmWhluSurPgL3AQufcO51ztwELgF3A3RMRnIiIyGQ2cE3QhSUzOFDdzEsHa7k4qP/icG5eVUBmcjz3vHYinKHKFPCXbeUY8LG3zOvbNis9idmZSewaoVLqjjIPi/LSSUmIC3OU0SkpPpbbzi3kgpIZPPb3l3PrGt/P83NTue+Nk5EOL+yqmgIFupIG7VtZmNn3/vG0dXOivm1aVdGV4RPGy5xzX3fOeQMbnM83gUvCH5qIiMjk1dPrpaa5o98HsMA6xrqWLi6aN/x01IDkhFjes7aIp/dUUd2kvozTXempNmqbOwdtd87x2K5KLp6f09dHL2Bl4fCFb5xz7ChtZE3x9E4C/u1dq/jTpy6lJDcV8I3w/80FxWw50cDhmqndcqPKX6ArL2NwwriqKJNKTweNnV52V3j6tsn0MVzCqF6LIiIiZ6i6uROvo98I45riLOJjfX9eLx5h/WKwD1w0lx6v44/TYKRDhveR/9nMJ+/ZMmj7/qpmjta2cuOqwRVOVxdncayuFU97d8hznqhvw9PezWpVvRzknecVERdj3L+5NNKhhFWlp4MZKfEhCx6t8r8vjnu8fSONKwqUME4nwyWMG83sX2zAfBkz+xqwKbxhiYiIRI97XjvOx/5n85gKzwS+sQ+sYQTftLeVhZnMSIkfcf1isJLcVK5cNJM/vnGS7l7vyHeQKcnT3s3hmhbePNnItpMN/fY9vrOSGIO3Lc8fdL/AaFCoBuzgW7/oO04J40Az0xO5dmkeD71ZTlfP1P3dq27qIH9AS42A5QUZmMExj5ddZR6Ks5OZoYI308pwCePfAyuBw2b2kJk9aGZHgNX+fSIiItPC+n01PLe/hleP1I/6PhWNodcEffWmZfzoPauJiRnbRJ4PXTyX6iZf6wOZnvYErUP8743H+352zvH4rkouXZBLTlrioPutKvQlggMbsAdsL20kKT6GRXmj/xJjOvmbC4s51drFs/um7u9epacj5PpFgNTEOBbOTON4k5ed5Y197yeZPoZrq9HknHs3cD3wP8DvgOudc7c754YvtSUiIjKFHKlpAeC/Xjk26vsEqg4ObIR93pwZXL0kb8wxXLVkFoVZyfzvq8fHfF+ZGgK98N55XiFP7Krse4/trWziWF3o6agAmSnxzM1JYWfpECOMpY2sLMwkLna4cYTp64pzZlKQmcR9U3haapWnI+T6xYCVRZkcbOil9FS7Ct5MQ8O11Xirmd3unDvinPurc+5R59wRM/uAmV03kUGKiIhESntXL+WN7cxIief5AzUcq2sd1f0qPO2kJsSSkTQ+VSdjY4yPXFbC68dO8eSuynE5p0wuu8o9FM1I5vPXLKLXOe7ZdBzwTUeNjTHeGmI6asCqoix2lg3+vr+718ueiiatXxxGbIzx7rXFvHyolk1HRz/LYLLo7OmlvrVryBFGgFWFmbT721Gq4M30M9xXSd8AXgyx/Tngm+EJR0REJLocrfONLn7umnOIi7FRj/BVNnaQn5k0qtYZo3XHpSWsKMzga3/ZzanWrnE7r0wOu8p8/e/m5KRw3dI8/vD6STq6e/3TUXOGbaS+uiiTCk/HoAqrB6qa6ezxsrpYCeNw3nfhHHJSE3nvrzbxN//5GhsO1IxpTXM0q2nyvScG9mAMtjLoCwUVvJl+hksYU5xztQM3OueqgNTwhSQiIhI9jtT6RhQvWZDLzasK+NOWUpo6QlebDFbZ1EFBVugiEmcqPjaGH9y+Gk97N9/4655xPbdEt8a2Lk6eOt3/7iOXzaOhrZtvPb6XE/Vt3DTEdNSAQEGbXQPWMQYK3qxRwjis/MwkXvziOr5641JO1Lfx4f/ezHv+87UpUQin0j+1OX+YKanLCzKIMZibk0JmSvxEhSZRYriEMcnMBs2jMbN4YHz/AoqIiESpIzUtfR+UPnLZPFq7enlgFGuZKhvbh53idaaWzs7g01ct5C/bK1QAZxrZXd4EnC5gc/H8bJbOzuDeTSeJizGuXzb0dFQ4/YF/x4B1jDtKG8lOTaBohj7ajSQ1MY6Pv2U+L33pKj5z1UI2H2/gYPXk789Y1RS6QFewpPhYFs+I4S3n5E5UWBJFhksYHwZ+bWZ9o4n+n3/p3yciIjLlHaltoTg7xdcSoyiTC0pm8L+vHafXO/R0tK4eL7UtnUOWqT9bf7duIUvy0/nKI7vwtI082imTX6DC6YrCDMDXVP4jl5UAcNnC3BHbHKQmxrFwVhrbSgeMMJb6prmO59TpqS4hLoYbVvoS9BP1bRGO5uwFWgANNyUV4IsXJPHNW1ZMREgSZYZLGL8KVAMnzGyrmW0FjgO1/n2jYmaxZrbNzB7z3/6Bme03s51m9oiZZQUde5eZHTazA2b21jN6RiIiIuPoSG0r83NPr8T4yGXzKD3VznPDlNivae7AOSgIwwgj+D6w/vDdq6lv7eLrmpo6Lewu9zAnO4WslNOJ4S2rC3jLObl9ieNIrluWx0sHa3lkWxkALZ09HKxpVsGbM1CS4/t/wvH60RXBimaVng5SE2JJTxp+qmmM2ZjbAcnUMFxbjR7n3D8DxcCH/f/mOOf+2Tk3lq8zPwfsC7q9HljhnFsFHATuAjCzZcB7geXA24BfmFnsGB5HRERkXHm9jmN1LSyYebo/3fXL8ijMSuY/nj885ChjYE3Q7HFewxhsRWEmn7lqIY9sK+fhN8vC9jgSHXaWeVg5oDplUnws93zsItYtnjWqc3z+2kVcMj+HLz+0izdPNrC73INzWr94JlIT45iZnsjxUVZNjmbVTR0jji7K9DZiwx3nXLtzbpf/X/tYTm5mRcCNwG+CzveMc85fmJdNQJH/51uB+5xznc65Y8Bh4MKxPJ6IiMh4qvC009HtZcGs0wljXGwMX3rbYnaVe7h/iLWMFY2+P5fhWMMY7O+vXsiF87L56p93c7S2JayPJZHT0NpFWUM7K8+y/118bAy/+MB5zM5M4s7fbeWp3VWA2iScqXk5qVNmhHFgv1iRYOHu0PpT4EvAUCWkPgo86f+5EAj+y1vm3yYiIhIRgQqpwSOM4JsKeNG8bL7/9H4aQrS3CDRUD3fCGBcbw7+/dw2JcTF85g/b6OzpDevjSWTsKvcVqlk1Dg3TZ6Qm8Ns71tLZ3cv/vHqc4uxkctISz/q801FJbgrHp8QaRo0wyvDGp5twCGZ2E1DjnNtqZutC7P8K0AP8PrApxGkGzfUxszuBOwHy8vLYsGHDeIUcVi0tLZMm1qlO1yJ66FpED10L8DpHzIDCH88c963AqDq0gw0n+u+7ucDL5uPd/MN/v8CHV/T/wL15byfJcbB108YxxXCm1+GOJTH89M0mPv2rZ/ngMn34Hw/R9Dvx1yO+LyUaju1iQ9n4rCH7xIpYfrK1h4KErqh5nqFE03UYyOvpora5myeffYHkuMm5ts/rHNVNHXQ2VI/4OkfztZhuJvpajJgwmq9s1geA+c65b5rZHCDfOffGCHe9DLjFzN4OJAEZZnavc+6DZnYHcBNwjTvd9bQM33rJgCKgYuBJnXO/An4FsHbtWrdu3bqRnkJU2LBhA5Ml1qlO1yJ66FpEj+l+LUpPtXHdT17kN397AZcHlY1f/8guslIqufm6dSGrSB5lL//96jH+8R0X9fW5A/jDyS0UZbeybt2VY4rjTK/DOqA5ZS+/feUYH7j6XK5aMro1bTK0aPqduK90KyU5Tdx43VXjds51wEXn11M4I5miGSnjdt7xFk3XYaC2nEoePPQmxUvP6+uPOdlUN3Xgffo5Llq1mHUXzx322Gi+FtPNRF+L0UxJ/QVwCfA+/+1m4Ocj3ck5d5dzrsg5V4KvmM3z/mTxbcCXgVucc8Hj+I8C7zWzRDObB5wDjJSUioiInLXXj52io9vLg1v7r0k8UtvC/NzUIVsOfP66c8hJTeRrf9mDN6gATlVTR1gL3oTy5bctYWZ6In/ZXj6hjyvht6vcE5aE5KL5OVGdLEa7QKXUydxao3KCps/L5DaahPEi59yngQ4A51wDMHyzn+HdDaQD681su5n90n/ePcADwF7gKeDTzjktxhARkbDb4e9N9+y+Gjq6T//pOVLbOmj9YrCMpHj+79uXsKO0kW8/sY/uXt+S/YrGjrC11BhKQlwMa4qz2FnuGflgmTTqWzopb2xXYZooNDfHl2xP5sI3gfXWeRlKGGVoo1nD2O1vb+EAzGwmQxexCck5twHY4P954TDHfRv49ljOLSIicrZ2lDWSlhhHS2cPrxyq49pleXjau6lt7uxXITWU284tZMuJBn77yjHePNnAj969mrqWzogUkVhZmMmz+6pp7ugesaeaTA6BgjcrC9X6ItqkJsYxa5K31qjyTExFZ5ncRjPC+DPgEWCWmX0beAX4TlijEhERmSCdPb3sq2ziby4oJiMpjid2VQL0takYboQRwMz4zm0r+Y/3ncvh6hbe/rOXASiIQJn6lUWZOAd7Kpom/LFlfDjnOFzTwksHa3lgSyn3bjoBwPLCjAhHJqGUTPLWGpVNHSTExpCdejaTB2WqG3GE0Tn3ezPbClyDr5LpO5xz+8IemYhErbufP8TS2RlcszQv0qGInLV9lc109zouKJmBp72bp/dU0dnTy9G+lhqpozrPzasLWF2UxWfv28b20sa+6WoTKdCnb1eZh4vn50z448vZu39zKf/88K5+295yTi4ZGjGOSiW5KbxwoDbSYZyxak8HeZmJQ67TFoFhEkYzyw66WQP8MXifc+5UOAMTkejU3evlZ88d5srFM5UwypQQWL+4ujiLxLhYHtxaxsbDdRypbSE+1ijOHn3iNycnhT996hJ2lDZy/twZ4Qp5SLlpiRRkJvVNY5TJxTnH/752gsV56fzrO1aQn5HErIxEkuJjIx2aDKEkN5XaLWW0dPaQlhi2bnVhU+npYHbGxM+GkMlluCmpW4Et/v/WAgeBQ/6ft4Y/NBGJRkdqW+jq9VJ6avJWhZOpobOnl3+4fzuHqpvP6jw7ShuZlZ5IfkYSly3MJT0pjid2VXGktoU52SnEx45m9cZp8bExrC3Jjtg39iuLMpUwTlI7yjzsq2ziQ5fM5cJ52czJSVGyGOUClVIn6zrGqqaOiKy3lsllyL+Czrl5zrn5wNPAzc65XOdcDr7+iQ9PVIAiEl32VfrWRpU1tHO6jeroff3RPVz0nWf5+P9u4dEjXbx8qJZe79jPI/LGsVM8sq2cJ3ZVndV5tpc1sqooCzMjIS6G65bl8cyeKg5UNY+4fjEarSzM5FhdK00d3ZEORcboj6+fJDk+llvXFEQ6FBmlydxawzlHpUcJo4xsNF+bXuCceyJwwzn3JDC2TsQiMmXsq/SN5rR09tDYNvYPpC8fqsUwjta18PChbj702zf44TMHxjtMmQY2Hq4H4LC/OM2Z8LR3c7S2lTXFp1sWvH3FbJo6ejhe3zZihdRotLLIV01zt0YZJ5Xmjm4e3VHBLasLVOF2EpnMrTUa27rp6vGSr5YaMoLRJIx1ZvZVMysxs7lm9hWgPtyBiUh02htUfbG0YWzfqDrnKGto5+bVs3n+n9bx82tSWLd4Jg9tLdMoo4zZq0fqAM5qSmogqVpdfLplwVsW5fatRZqsI4zgK3wjk8dftlfQ3t3L+y6aE+lQZAwmc2uNbaUNABTO0BpGGd5oEsb3ATPxtdZ4xP/z+8IZlIhEJ+cc+yqb+hpIl55qH9P9a1s66ezx9hURSY03bj+/iJrmTt44pjpaMnqNbV3sKveQGBfD0brWM/7CYbu/4M2qoB53iXGxXLt0FjD6CqnRJDs1gcKsZK1jjGKvHKrrtw7cOccfXj/J0tkZrC7KHOaeEo1Kcidfa41er+N7Tx5gTrbvi1uR4YyYMDrnTjnnPuecO9c5d55z7vOqkCoyPdU2d1Lf2sV1/uqoJ8dY+KaswZdgFgV9m3n1klkkx8fy2M6K8QtUprxNR+txDm5dU0BXz5kXYdpR2sj83FQyU/pPAbzj0hIunp/NkvzJ2ftulQrfRK1Hd1Twwd++zg3//jJ/3lYOwK5yD3srm3j/hcVqbzAJleSkcHySrWF8YEspB6qbueuGJSTGqbCSDG9spd9EZFrb6y94c8G8bLJS4sc8JfV0wni6TUFKQhzXLsvjyd1V9PR6xy9YmdI2Hq4nJSGWd51XBMDhmjNbx7ijrLHfdNSAc+fM4L47LyE5YXJ+kFpZlMmJ+jY8Z7DOWMJne2kjX/zTDs6fO4Ml+el8/v7t/NMDO/ivV46RFB/DrecWRjpEOQMluanUNnfS0tkT6VBGpaWzhx89c4ALSmbwthX5kQ5HJgEljCIyaoGCN0vzMyiekTLmUZ3A8YVZ/ddL3LRqNqdau3j1iJZHy+hsPFLHhfOyWVrgGwE8dAYJY5Wng+qmzr4p1lNJYB3j7gqNMk6EI7UtvPMXG/nh0wc4MkQRpipPB3f+bgsz0xP51YfO5747L+bvr17Iw9vK+PP2Cm5eVUCGit1MSpOttcYvXjhMXUsXX71xmUa0ZVQmX4dREYmYfZVNFGYlk5kST3F2cl8COVplDe3kpCaQOqC58ZWLZpKeGMdfd1RwxSKtpZDhVXk6OFrbyvsvnENGUjx5GYlnNMIYWL8YaoRxsgskjDvLPFy2MDfC0Ux9T+2u4s2TjWwvbeTuFw6zujiLt6/IZ3F+OgtmppGdmsAnfreF1s4e7vnYZeSkJQLwT9cv5pIFOfz8hcN88sr5EX4WcqaCW2usKIzuL6DKGtr4zSvHuO3cwin5/z4JjxETRjObCXwCKAk+3jn30fCFJSLRaF9lE0tn+0Z0imek8OzeGrxeR0zM6L6hLGto67d+MSApPpbrlufx9J4qvnXbCq2nmEa6e7088mY5t6wpGHWD8o2HfdVRL13gS4QWzkrjcM3YK6XuLGskLsZYNntyrlMcTlZKAsXZyWqtMUG2nfSthb3vzov5y/YKHnqzjO8+ub9vf2yM4XWO3/ztWhbnp/e776ULcvveyzI5leROntYa33/qAAZ88a2LIx2KTCKjGWH8C/Ay8CzQG95wRCQSenq9PLytnEgbkeUAACAASURBVHeeW0hcbOiZ6h3dvRypbeEG/3qH4uwUunq9VDd3MDvzdBLoq6TazLKCwR/Cyxva+xLOgW5eXcDDb5bz8sE6rl2WNw7PSiaDZ/ZU86WHdlLhaefz1y4a1X02Hq4jOzWBJf4P3ufMSudPW0pxzo1petWOskaWzs4YdaI62awqzGJneWOkw5jynHNsL23gikUzmZWRxCeumM8nrphPfUsnR2pbOVLbwpGaFtaWZHPNUv2/bSpKSZgcrTWO1bXy6I4K/m7dAgqy1EpDRm80axhTnHNfds494Jx7KPAv7JGJyIR5+XAdX3pwJy8dqh3ymIPVzXgdp0cY/a0xBrbWeHpPFW//2cv9+jUCeL2+HoxF2aH/SF2+MJeslHhVS51mNvp7Kf7ni0ep8nSMeLxzjo1H6rhkQU7fyPaCWWm0dvVSOYr7B59nZ5lnSq5fDFhZlEnpqXYa27oiHcqUVtbQTl1LF+fOmdFve05aIhfOy+Z9F87hqzctU3GRKW4ytNb43WvHiYsxPnxpSaRDkUlmNAnjY2b29rBHIiIRc7TW90fuWN3QRWz2+Suknp6S6kv8Bha+2XTU13VnR1n/kY3alk66er39KqQGi4+N4YYV+azfW01HtyYzTBevHalneUEGvV7HD54+MGh/d6+33/vhSG0r1U2dXBY0he+cWWnA2CqlljW009zRw/L/n737Dm+zvBo//r0lW957j3hm2ImzB9kkISlhj0ILlFLogI630L7t274tnXS3tKWl/bUFWqDQt5RdGiBAErITMknixM7wSrz33tL9+0OyE8eOrSSWH0s6n+vSRSw9kk54ZEVH97nPSfTghNGxl+rgaVlldKWDjr2ws2U/mFfLSQxjb3EDD792hKaO8deduK2rl5f3lXLt9ARiQ/2NDke4GWcSxoewJ42dSqkWx6V5xHsJIdxGXxlNyTDfjuZVtBBkMZPiWFlMighAqcGzGA+cbgAYtMJY6hjBMdQexj43zEikrdvKO0crL/4vIdxOeWMHRbVt3DI7ifuWpvHKgVKOlJ7dc1fV3MmNf9jB/B9v4Odv51Pd0slOx4rkkolR/cdNdCSMF9Mpte8LkKyEkBGOdF9zUyOICPTlhb2njQ7Fox083YC/r2nQ3kThXb5+9WQ+szSdf+45zVW/3sJ/DpWjtR7zONq7e4f80vXVA6W0dPXyKVldFJdgxD2MWmt5BxTCwxU5EsbhBg8fq2hmSnxIfxmgn4+ZuBD/AbMY27t7OepIFPs+kPfpm8E4YZiEcWFGFGlRgTy7s5ibZsk8Mk+3yzFGZXFmNMmRAby0r5Qfv3mMF+5fyKnqVu59ei+N7d0smRjNE1sL+NuOIiICfUkKD+j/4gIgKshCRKDvRa0w5le2oBRMifPcf+L8fc3cuSCFP28p4Ex9e38ZuRhdH55pZEZSOL4X2P8tvEOgxYfvXj+VW2Yn8a1Xj/Dlfx7kwOkGvn/DNJc+b1N7D394/yTHq+x7ZcsaO4gI9OWlzy/u/zJNa82zu0qYnhTGnBRZCRcXz6l3N6XUjUqpRx2X610dlBBibBWNsMJob2TTPKhhzYTIAErP2cN4uLQJq02TEhlIXkUzNtvZb1f7SlcvVJIKYHLsrehrTy88286Cuv7mNaH+vnx1zWQ+KKrn0XeP89E/7aTbauNfDyziiXvmselrK7htbjINbT2smRo3oLmNUopJsSEX1Sk1v7KZ1MjAQSNePM0nF6WilOLZncVGh+KRunqtHC1rZrZ8CBcOOUlhvP6lJdy5YALP7CweVG0z2l45UMqT24qob+tifloEX109GbNJ8dln9/bvX95xqo5T1a18anGazF0Ul2TEhFEp9XPsZanHHJeHHNcJITxAZ4+VssYOLD4mShs66LHaBh3Tt99rUMIYEThghbGvHPWuK1Jo67YOuK20oYPoYL8RO1LeNm8CIX4+PL2j6HL+WmKc01qzs6CWRRlnm9fcOX8CE2OD+eP7BUSH+PHqFxb3zzRLiw7ip7dM58D31vDwddmDHi8zNpiT1a1Ol4DlVbSQFe954zTOlxAWwLXTE/jXvjO0dvUaHY7HOVbeTLfVxizZvyjOYTYp/veabMIDfPnxm8dcWpq641QtaVGBrPvyMh67YzYPrZ7EXz45l/LGTr74jwP0WG08s7OYyCAL189IcFkcwrM5s8J4LbBGa/03rfXfgLWO64QQHqDEUYa6MCMKq6OT6fnOb3jTZ0JkIJXNnXT12vdLHChpICMmiMWZ9v1l536zWtrQMez+xT7Bfj58bP4E3jxc4VTXTOGeiuvaqWjqZFHm2b2IPmYTj94+kzvmT+CVzy8esoQy2M9nyNK/SbHBNLb3UNc2ckfQ9u5eiuvaPHr/4rnuW5JGS2cvr+wvNToUj9NXCXF+h1QhwgJ8+crqyewsqGNTfrVLnqPHamN3YR1LJg6c4zk3NZKf3TqdnQV1PPjPg2zMr+LOBRM8doSQcD1nC+7P/erMc1vKCeGFimrt+75WTokBhh48nFdh3++VdV5ThwmRgWgN5Y2daK05cLqRuSkRTI4LwaTs+x77lDY4v4fqU4vSsGrN87tLLvWvJca5vuY1i89JGAFmTQjn5x+dQUSQ5aIer7/xTdXI+xhPVLWi9eAvQDzVnJQIZk0I55mdxQPKxMXlO3i6kYQwf+LDpOukGOyuK1LIiAniJ2/lDVm9c7kOnWmkrdvK0vMSRoCPzk3mgSszeDu3EpNS3L0wddSfX3gPZxLGnwEHlVLPKKWeBfYDP3VtWEKIsVLo2L+4YkosACVDDB4+VtFEWlTQoP1e547WKK5rp76tm7mpEfj7msmMCe5fmbTaNGWNzq0wAqREBbI6O45/fFAiIzY81M6COuJD/UmPDhqVx5sU5xitUTNywti/Yu4FJal97luSRlFtG5tPuGalw1sdPNMg+xfFBfmaTTx8bTaFNW38wwVfgG4/VYtSDKjUONc3rs7izgUT+OyydBLCnPv3V4ihjJgwaq3/CSwEXnVcFmmtX3B1YEKIsVFc20ZMiB9pUYEEWcxDdkrNLWvu30t2rr4VwzMN7ewvse9fnJtqL82amhjaX5Ja3dJJj1U7nTACfHpJOg3tPfz7w7KL/juJ8c1m0+wuqGPxxKhRa8AQH+pPsJ8Pp6rsjW+sNs13Xj/CjX/YTu953+znVzQTZDFf1OvR3V07PYG4UD+e3lFsdCgeo7a1izP1HbJ/UQxrVVYsSyZG8djGkzS1Dz2f8XhlC3/aXHDRex13nKplelIY4YFDV2SYTYqf3TqDb10zeN+3EBfjggmjUirL8d85QAJQCpwBEh3XCSE8QFFtG+lRQSilSI0KGtQptba1i7LGDmYMkTDGhfrja1acrrcnjKH+PmTG2Fd6shNCKW/qpLG9u39f5HAdUs+3MCOSrPgQnt5RbMgsK+E6x6taqGvrZnHm4DKqS6WUIjM2mFM1rXT32vjyPw/w/O7THC5tYp/jy4w+eZUtA0bEeANfs4l7FqWx7WSty7s2eosPT8v+RTEypRQPXzuV5o4evvnK4UFl4bWtXdz79B5+sT6/v2O5M1q7ejl4unHIclQhRttwK4z/7fjvr4e4POriuIQQY6Sotq2/LDAtOrC/CU6fvkHq05MHJ4xmkyIp3D5a40BJA7NTIvo/hE917A87VtFMqaNb6nAzGM+nlH3ERn5lC0fKmka+g3AbO/vnLw5dRnWpJsYEc7yyhc/9fR9vHanka2sm4+djYn1uZf8xWmvyhxgR4w3uviKVIIuZP20pMDoUj3DwTAM+JkVOorR2EMObmhjKt67JZv3RSn72dl7/9T1WG1/6xwFqW7sA2FVY5/Rj7imqo9emJWEUY+KCCaPW+n7HH6/RWq8894J0SRXCIzR39lDb2k16jD1hTI0K4kxD+4ASvsOlTSjFkCWpYC9Lzato5kR1S385KpxtKHKsvJkzjlmNieEXVwL4kWnxmBRsyJN9V55kV0Et6dFBF/16GMmkuGBqW7vZerKGn986nS9fNYllk2J492hl/yp1eVMnzZ29ZHlhwhgW6Mvdi1J583D5Ra1kiKF9eKaR7IRQAizSeVKM7LPL0rlnUSpPbivi77uKAfjpW3l8UFTPL2+bQXyoP7sKnE8Yt5+sw8/HxJxUWeEWrudM05udTl4nhHAzxY4PjWlRjhXGqEB6rJqKc8ZZHClrJDMmmOALDDhPjgiksLYNre3dGPvEhPgRE+JHXkULpQ3txIaMPIPxfJFBFuamRrDhWNXF/tXEONXe3csHhfUXbNJwORZnRhEZZOH3d8zmjgUpAKzNiae8qbN/lTq/v+GNd4zUON9nl2bgazbx582yyng5rDbNoTNNsn9ROE0pxfdvmMbq7Fh+8MZRvv3aEZ7eUcynl6Rzy+xkFmVGsbuwzuktGDtO1bIgPVJGZYgxMdwexnil1FwgQCk1Wyk1x3FZATi/EUkIMW71rTJknLPCCANHaxwubRpy/2KfFEfjG5OCmRMGHpedEOooSe1weqTG+VZnx3GsopmyxsHzIYX7+d3Gk7R09fLROcmj/tgzksPZ/53V3DAzsf+61dmxmE2qvyw1v9LeFGeKlyaMMSF+fHz+BF49WEq5/E5dEptN8/BrR2jt6mXZJCkHFM4zmxS/v3M2OUlh/N8Hp1mYEcm3rs0C7J1Oa1u7OVk9cqfn6pZOjle1DJq/KISrDLfCeDX2vYrJwG84u3/xv4Fvuz40IYSrFda0odTZpC+tP2G07zmsau6kuqVryP2LfSZE2ssKp8SHEuLvO+C2qQmhnKpuoai27ZI7Uq6eGgfAxjxZZXR3+ZXN/HVbEXfMnzCgfHk0nd91NTzQwsKMSN45ak8Yj1U0MyEyYNBr1ZvcvzwDreGJrYVGh+J2tNZ8/42jvLD3DA+umshHpsUbHZJwM4EWH/76qfn818qJ/PGuOfia7R/FF2XYqy6cKUvdcco+x1b2L4qxMtwexmcd+xXvPW8P441a61fHMEYhhIsU17WRFB7QX9JiLxs19c9iPOxoeDNjuITR0fl0zhCzyKYmhvaXuF5qwpgZE0xGdJDsY3Rz9lWZXEIDfPnm2qwxfe610+IpqGnjVHUL+RXNZHnR/MWhJEcEcvPsJF7Ye7q/2YYYmdaaH63L47ndJTxwZQZfXTPZ6JCEm4oJ8ePrV08hKtiv/7oJkYEkRwSws6B2xPtvP1lHeKBvf3M5IVzNmTmMryilrlNKfUMp9b2+i7NPoJQyK6UOKqXWOX6+XSl1VCllU0rNO+c4X6XUs0qpI0qpPKXUty7trySEcNa5HVIBTCZFamRQ/wrjkdJGTAqmJlw4YcyMDSYlMpC1OYO/aZ+acLbs72JGapxv9dQ4dhXU0tI59AwrMf69uO8M+0sa+Pa12UQEDT0zzFX6VoFeP2hv9uKt+xfP9YUVmXT12vjb9iKjQ3ELWmt+sf44f9tRxH1L0vjftVmjNkNUiD6LM6PYXVg/aPTGubTW7DhVy5LMaK8aDSSMNWLCqJT6M/Bx4MuAAm4HUi/iOR4C8s75ORe4Fdh63nG3A35a6+nAXOABpVTaRTyPEOIiaK0pqhmYMAKkRgX2z2I8VNrE5LiQYbsABvv5sPUbK1k2KWbQbenRwfj72t9mJlxOwpgdR49Vs+3kyN+8ivGnrrWLn72dzxXpkXx0TtKYP39cqD+zU8J5dlcxNo1XjtQ4X2ZMMNdOT+C5XSU0yxcxI/rthpP8eUsBdy9M4XvXT5VkUbjEoswomjp6OFZx4Vmpf9tRTGVzJ1dOGfxvrhCu4kyX1MVa63uABq31D4FFwARnHlwplQxcBzzVd53WOk9rfXyIwzUQpJTyAQKAbkCmCwvhIrWt3bR09Q5KGNOigyipb8dm0xwpa2L6MA1vRmI2KabE2VdzLrUkFezlruGBvtIt1U39Yn0+7d29/OSWHMM+aK+dFk9LZy+AV47UGMrnl2fS0tXLC3tOGx3KuPaHTSf5/caTfHzeBB650bjXsPB8izLsexJ3X2Ae4+sHy/jRumOsnRbvksZhQlzI0H3yB+rrr9+ulEoE6oB0Jx//MeAbgDP1Py8DNwEV2LuwflVrXX/+QUqp+4H7AeLi4ti8ebOToRirtbXVbWL1dN5yLlq6NS3dmsTgob8XOtFgtR9XXsDmzSX913fX9dDda+PJ1zdR39ZNQEf1Zf3/CqcLBZw8vIfi88pnLuZcTA3XvJNbxsaYBsxShjPqXPV70W3VvH6gnSVJPpQe20/psVF/CqeEt9lni1rMUHRkDyXj9EP/WL8/ZUea+NPG42T0nsZHfq8GaG1t5ZtPv8e/jnezONGHq6Pq2Lp1i9FheR1v+Te7T1yg4j97TjDROvCLnMM1vfzuQBdZkSZuTWxmmwGvRW87F+PZWJ8LZxLG/yilwoFfAQewrwQ+OdKdlFLXA9Va6/2OURwjWQBYgUQgAtimlNqgtR7Qxk1r/QTwBMC8efP0ihXOPLTxNm/ejLvE6um85Vx8+7UjbMqrZve3Vwx5e/XeM/DBYW5atah/nAaA5VQtzxz9gFJTHFDCR1fOZ+ZlzBpLzWnj0JlG1sweXIp4MeeiPaqCnf84QEj6TBakR15yPGJorvq92JRfRbdtH59eM4flk40toXr21DYCLGZWrVxsaBzDGev3Jx1fzX3P7KU5fBK3yorFAA8/Y08Wr5+RwGMfn4WP2ZmiLDHavOXf7D5XNRzhjQ/LWbpsef9r7uDpBv608QOmxIfyrwcWGtbl2dvOxXg21udi2Hc/pZQJ2Ki1btRav4J972KW1tqZpjdLgBuVUsXAC8AqpdTzwxx/F7Bea92jta4GdgDzhjleCDGMU9WtVDZ30tbVO+TthbVt+JoVSeEDS0VTHSWqbx2pwNesyEq4vAYh6dFB3DxEsnixlk+OwWI2sUHGa7iVDXnVBFnMXJFhfJL/l0/O5bGPzzI6jHFlxZQYJscF88TWQqcHhnuD53eX8I/8bq6eFsdvJVkUY2hRRhStXb3kljfTa7Xx5NZC7nryA2JD/Xj20wu8eiSQMM6w74Baaxv22Yt9P3dprZuceWCt9be01sla6zTgDmCT1vruYe5yGntSqZRSQcBCIN+Z5xJCDHba0en0QgPvi2pbSYkMHPRBKCHUH4uPibq2bqbEh+Dnc+GGN2Mp2M+HhZlRvLK/lHeOVsqHWzegtWZTXjXLJsWMi9fRhMhAJkReevMlT6SU4nPLMsivbGGrNJUC4MW9Z/jO67nMjDHz+J1n5+QJMRYWOuYx/n1XMTf+YQc/eSuPxZlRvHD/QmJC/Ia/sxAu4sy74LtKqY+qUdrlrZS6RSlVir15zptKqXccN/0RCMbeRXUv8LTW+vBoPKcQ3qazx0pls337cWlD+5DHFNe2D2p4A/bRGimOD9XTky69FNUV/ndtFpFBFh54bj+feOoD8iulL9Z4drS8mcrmTq7KjjU6FDGMm2YlERfqxxNbC4wOxXCvHSzlm68eZvnkGP5rth8WH0kWxdiKCfFjclwwrx4oo7a1iz99Yg5PfWoeCWGX3jhOiMvlzDvhfwMvAV1KqWalVItS6qI+pWmtN2utr3f8+TXHyqOf1jpOa3214/pWrfXtWutpWuupWutfXfTfRggBDEwSSxsGrzDabJqiusEjNfqkRdkTxhnJl94h1RWmJoby9kPLeOSmaRyraOba323j2Z3FRoclLmBjXjVKwcosSRjHM4uPifuWpLPjVB25ZU4VEXmkzcer+dqLh1iUEcUTn5yLrzQBEgb56urJfHFFJhu+diXXTE+QzrzCcCMmjFrrEK21SWtt0VqHOn6WnuRCjGOn64dPGMubOujutZEeHTzk/fua4FzOSA1X8TGbuGdRGpu/voI5KRH8ZUuBlKeOU5vyq5g1IZzoYCmjGu/uXJBCkMXMb9474ZW/T129Vr7/xlEyY4J56lPz8Pc1voRaeK9rpifwjbVZhMp+RTFOjJgwKqU2OnOdEGL8KHHsXwwP9B2yJLW41n5dWvTQ+7muyopl+eQYpsRfXsMbVwoPtHD7vGTKmzrJq2gxOhxxnurmTg6VNnGVrC66hbAAX766ZjKb8qv5+66Ske/gYf6+s4SSuna+e/1UAi3ONJAXQgjvccF3RaWUP/Z5iNFKqQigbz08FPvoCyHEOHW6vp1Ai5npSWFDrjAW1bYCkHGBFcbFE6NZPDHapTGOhr5Sx035VUxNlMKH8eT949UAXJUdZ3AkwlmfWZrOzoI6fvJmHvPSIpiWOP4qDFyhvq2b3286yYopMYaPfhFCiPFouBXGB4D9QJbjv32Xf2NvUCOEGKdO17WT4ugIOVTCWFjbRoCvmbhQ9y4VjA3xZ+aEcDbkVRsdijjPhrxqksIDyBrHq9RiIKUUv7ptBhFBvnz5/w5ecCSPp/ndhhO0d1t5+Npso0MRQohx6YIJo9b6d1rrdODrWusMrXW64zJTa/2HMYxRCDGEyqZOZv7wXfYV1w+67XS9I2GMCKS+rXvQB7+iWnvDG0/YSH9VViyHShupaekyOhTh0NljZfvJWlZlxXrEa8ybRAX78djHZ1Nc18Z3/51rdDgud6q6lec/OM1dC1KYFCdfbgghxFCcaXrzuFJqsVLqLqXUPX2XsQjOHe0rrueev+2ho9tqdCjCw+0vaaCpo4etJ2oGXG+zaU7Xt5MaFUhyhL0N9/mzGItrL9wh1d1clR2L1mdLIIXxdhXW0dFjlXEabmpRZhRfXjWJVw+UsT63wuhwXOqnb+UR6GvmK6snGR2KEEKMW840vXkOeBRYCsx3XOa5OC631d1rY+uJGn6xPt/oUISH65tBeOS8NvjVLV109dpIiQrqTxjPbXzT3WvjTEOHxySMUxNCSQjzZ2NeldGheI3OHis/fzufgprWQbdZbZpndhQTaDH3D6AW7ufBqyaRFB7Ai/tKjQ7FZQ6cbmBTfjVfWjWRKOnkK4QQF+RMK7B5wFTtjX22L8HiidHcuziNZ3YW85GpcW7ROES4p7yKswmj1rq/9K9vpEZKZCDJEfYuqOfuYzzT0I7Vpj0mYVRKsSorltcOltHVa8XPR9rhu9r2k7X8eUsBrx0s5cUHFvWPYdFa86N1x9hyooZHbpomowncmNmkWJsTz3O7Smjp7CHEA9v7v/FhOX4+Ju5emGp0KEIIMa6NuMII5ALxrg7Ek3xzbRYZ0UH8z8uHae7sMToc4aHyKlrwMSlqW7upaOrsv76krg2A1MhAooMt+PmYBiSMxbX229M8JGEEe1lqe7eV3YWD93OK0benuB6L2UR3r427nvyAckfJ81+3F/HMzmI+uzSdexalGRukuGxrc+Lpttp4/3jNyAe7GatN89aRClZOiSXYT8ZoCCHEcJxJGKOBY0qpd5RSb/RdXB2YOwuwmHn0YzOpaOrgR/85ZnQ4wgM1dfRQ1tjBKsdYiXPLUs/Ut2NSkBgegFKK5IiAASWpRY6EMcODEsbFmdH4+5rYJGWpY+KDonpmTQjn75++guaOHj7x1Ac8v7uEn7yVx9pp8Xxbuk16hLkpEcSE+HnkPsZ9xfVUt3Rx3YwEo0MRQohxz5mE8QfAzcBPgV+fcxHDmJMSwRdWZPLS/lI2HJMPsWJ05TvKUW+dk4TZpDhSejZhLKlvJzE8AIuP/dc7OWLgaI3C2jbCA32JCLKMbdAu5O9rZunEGDbkVSPV867V1tVLblkTC9IjmZ4cxjOfnk9VcyffeT2XWRPCeeyOWZhM0hnVE5hMio9MjeP9/Bo6ezyrkdubRyrw9zX1f+kmhBDiwpzpkroFKAZ8HX/eCxxwcVwe4aGrJpOdEMr33zgqH2LFqOrbvzg7JYJJscEDVhj7Rmr0sa8wnk0Yi2o8p0Pqua7KjqWssYMTVYMbsYjRs7+kAatNsyA9EoC5qZH87d753DAzkafumSf7Fj3MNTkJdPRY2XLCc8pS7eWolazKiiVIylGFEGJEznRJ/RzwMvAXx1VJwOuuDMpTWHxM3Lc4jbLGDspaJWEUoyevooXIIAuxIX7MSA7rb3wDcLrOPlKjT/J5sxiL69pIj/LAhDErFqXgZ2/n0d7tHQPHjbCnqB6zSTEnNaL/uoUZUTx+52zpNOmBrsiIJCzAl3dyK40OZdTsKaqntrWL66YnGh2KEEK4BWdKUr8ELAGaAbTWJwGp4XDS0kn2Lqm5tZ5VziOMlV/ZTHZCCEoppieFUd/WTXlTJ61dvdS1dTPhvBVGsM9ibO/upaKp0yNXGGND/fnRTTlsPVHDXU9+QF1rl9EheaQ9RfXkJIZKoxAv4Ws2sTo7jg15VXT32owOZ1S8eaScAF8zK7NijA5FCCHcgjMJY5fWurvvB6WUDyDLZU5KDA9gYmwwuXWSMIrRYbVpjle1kBUfCsD05HAAjpQ2crrO3twmNfJsQnjuLMbiWvvt6TGelzAC3L0wlT/dPZe8imY++qed/R1jxejo7LHy4ZnG/nJU4R2uyYmnubOXXYV1Rody2XqtNtbnVrIqO5ZAi3zpIYQQznAmYdyilPo2EKCUWgO8BPzHtWF5lmWTojleb/W4pgHi8pU2tPP55/bT0NY98sEORbVtdPbYyE6wJ4xZ8SH4mBSHS5sGzGDsc+4sxmJHApXmgSWpfa6eFs//fe4KGjt6uPX/7aSquXPkOwmnHDrTSLfVxoL0KKNDEWNo6aRogixm1ntAWaq9HLWb66dLd1QhhHCWMwnj/wI1wBHgAeAt4DuuDMrTLJ8UQ48N9hU3GB2KGGfeP17D+qOVPLe7xOn75FfaG95kJ4QA9g6hk+NCOFLWxOl6e0KYcs4exnNnMfaN1PDEktRzzU2N5Jn7FlDX1s1WD2rWYbQ9RfY5l/PTIkY4UngSf18zK7Niee9YJVabexcYrTtSQaDFzIopsrNGCCGc5UzCGAD8TWt9u9b6NuBvjuuEk67IiMRHwbaT8sFVDFRYY+/o+dzuEqf3B+VVNONjUkyMDe6/rq/xHjqsuAAAIABJREFUTUldO2EBvoQF+Pbfdu4sxsKaNuJC/byiM+CMpDCCLGaOljcbHYrH2FNcT1Z8COGBnjOSRThnbU48ta3d7CuuNzqUS9bjKEe9KjuOAIt08xVCCGc5kzBuZGCCGABscE04ninQ4sOkCBNbT9YaHYoYZwpq2vD3NVHT0sWbR8qduk9eRQuZMcH4+Zz9wDM9OYzG9h52FdQN6JDap28WY1Ftq8evLvYxmRRTE0PJPWfkiLh0PVYb+0sauEL2L3qlFVNisfiYWH/UfctStxyvob6tmxtnSndUIYS4GM4kjP5a6/7BZo4/D/5EKoaVE20mr6KZ6hbZTyXOKqhuZc3UeDJjgnh6R7FT8zrzK5r7y1H7TE8KA6Cwtm3A/sU+fbMYi+vavSZhBJiWGMbR8ma3L6MbD46WN9PebZX9i14q2M+H5ZNieCe30m3nCr+8v5SoIAsrpkh3VCGEuBjOJIxtSqk5fT8opeYCHcMcL4aQE21fDdpxSlYZhV1Ht5Wyxg4mxQZz75J0Dpc2sb9k+H2uje328Rl9DW/6TIkPwdesAC6QMNpnMda3dXtVwjg9KYyOHitFta0jHyyGtafI3iFzfrrsX/RWa3PiKW/q5HCp+63a17d1szG/iptnJ+FrduajjxBCiD7OvGt+BXhJKbVNKbUN+BfwX64Ny/NMCDERFWRh2wlJGIVdoSOJyYwJ5qNzkgj19+HpHcXD3ievogWArPMSRj8fc/+YjaFLUs9WladHBw+63VPlOFZec8tkH+Pl2lNUT0Z0ELEh/kaHIgyyOjsWH5PibTfslvrGh2X0WDW3zU02OhQhhHA7IyaMWuu9QBbwBeCLQLbWer+rA/M0JqVYOimarSdr3bacR4yughp7x9KMmCACLT7cuSCF9UcrKWu88AJ+XsXADqnn6kuOJgyxwnjudenR3lNRnhkThL+viSOyj/Gy9Fpt7CmqZ36a7F/0ZuGBFhZlRrE+t8Lt/h17+UApOUmhg6ozhBBCjMzZuoz5wAxgNnCnUuoe14XkuZZNiqG2tYv8yhajQxHjQEF1K0qdHXHxyUWpaK35+67iC94nv7KZ6GDLkKs8V6RH4mse2D21T98Ko0kNnVB6Kh+ziewEaXxzuQ6eaaS5s5flk2Xvl7dbmxNPcV07J6rcp8w7r6KZ3LJmbpsjq4tCCHEpRkwYlVLPAY8CS7EnjvOBeS6OyyMtmxQNwNM7inhx7xn+ur2IP75/atgVJeG5CmvbSI4IwN/Xvr81OSKQtTnx/POD0zS0dQ95n7yKlv7S0/PdNCuRLf+zcshkMirIgr+vieSIwAHdVb1BTmIYx8qbsUnjm0u2Kb8aH5Ni2eRoo0MRBlszNQ6l4O3cCqNDcdor+0vxNStunJVkdChCCOGWnBnGNg+Yqt2t/mQcigv1JycplBf3lfLivtL+63cX1vH3Ty9AKWVgdGKsFVS3khkzcDXwK6sn8+7RKn75znF+duv0AbflVTSTV9HM55ZnDPl4SikSw4cekaqUIi0q6IK3e7KcpFCe211CSb13dYgdTZvyqpmfFkmov+/IBwuPFhviz7zUCNbnVvKV1ZONDmdEPVYbr39YxlVZcUQGyfxQIYS4FM4kjLlAPOA+XyeOY89/5goqmzsJ9vMh2M+Hl/eX8uM389h8ooaVU2KNDk+MEZtNU1jbyqLMgSMKJseFcO/iNP66o4g75k9g5oRwALp7bfz3i4cID/Tls0vTL+k5H79zdv9qpjeZltjX+KZJEsZLUNbYwfGqFh6+NtvoUMQ4sTYngR+tO0ZxbRtp4/x3asvxGmpbu6XZjRBCXAZn9jBGA8eUUu8opd7ou7g6ME8VHmghKz6U5IhAwgMt3LMojfToIH7yZh49VpvR4YkxUt7UQWePbdAKI8BDqycRHezH9/6d219G+buNJ8iraOZnt84gKtjvkp5zUlyIV+1f7DM5LgSL2URuuexjvBSb8qsBWJUtX2gJu6unxQGw/uj47paqteYfH5QQHWzhSpm9KIQQl8yZhPEHwM3AT4Ffn3MRo8DiY+Jb12RxqrqVF/acNjocMUb6OqRmxgz+dj7E35eHr83mUGkT/9p3hgOnG/jT5gJun5vMmqlxYx2q27P4mJgSHyKNby7RprwqUqMCyRjnK0li7CRHBDIjOYy3j4zvbql/3V7E+8dr+PTSdJm9KIQQl8GZsRpbgHwgxHHJc1wnRsmaqXEszIjktxtO0tTRY3Q4YgwUVNs7DGYMscII9gY2C9Ij+eX6fL724iESwgL43g1TxzJEj5KTFEpuWfO4/nA7HnV0W9lZUMfKKbGyx1oMcOPMRA6VNvH55/dT29pldDiDbMqv4qdv5XHt9Hg+vzzT6HCEEMKtOdMl9WPAHuB24GPAB0qp21wdmDdRSvGd66bS0N7NH98/ZXQ4YgwU1LQS6u9DdPDQTRiUUjxy0zSaO3spqm3jV7fPIEQajlyynKQwmjp6KG2QjsQXY2dBLV29Nq6SclRxnvuWpPOta7J4P7+Gj/x2K+vHUdfUE1UtPPjPD5maGMqvb5+FySRfdgghxOVwpkbjYWC+1vpTWut7gAXAd519AqWUWSl1UCm1zvHz7Uqpo0opm1Jq3nnHzlBK7XLcfkQpNXg+gIfKSQrjtjnJPL2jiJNVMqfR0xXWtJEZGzzsqk1WfCg/vjmHR26axuJMGWdwOXIcjW+Oyj7Gi7Ipv5pAi5kF6ZFGhyLGGbNJ8cCVmax7cClJ4QF8/vkD/OCNo0aHRX1bN595di+BFjNP3jOPAIv3NfoSQojR5kzCaNJaV5/zc52T9+vzEJB3zs+5wK3A1nMPUkr5AM8Dn9daTwNWAF5Vn/k/V08hLMCXz/19H43tQ8/hE56hoGbwSI2h3LkghXsWpbk+IA83JT4Es0lxRPYxOk1rzfv51SydGO11szuF8ybHhfDqFxfziStSeGZnMXuL6w2N58frjlHV3MUT98wjIcz7xggJIYQrOJP4rXd0SL1XKXUv8CbwtjMPrpRKBq4Dnuq7Tmudp7U+PsThHwEOa60POY6r01pbnXkeTxEb6s9fPjmX8sZOvviPA9I11UM1d/ZQ3dLlVMIoRoe/r5lJscHkljUbHYrbyK9sobypU8pRxYh8zSYevi6b+FB/frzuWH9357FW39bNusMV3LUghVmOkURCCCEun3KmCYRS6lZgKaCArVrr15x6cKVeBn6GvVnO17XW159z22bHdfscP38FmAvEAjHAC1rrXw7xmPcD9wPExcXNfeGFF5wJxXCtra0EBzuXIOwo6+HJI92snODDPVMt0mxilF3MuXCFwkYrj+zu5MHZfsyJc2YUqucay3Px1yNdfFjTy+9XBsrv1BDOPxfrCrp5+WQPj60IINxfOkyOFaPfny5H379d98/wY3Hi2L+3vVXUzYvHe/jJkgCSQi7/NevO58KTyHkYP+RcjB/nn4uVK1fu11rPG+Yul+WC7+hKqYlAnNZ6h9b6VeBVx/XLlVKZWuuC4R5YKXU9UK213q+UWuFkLEuB+UA7sFEptV9rvfHcg7TWTwBPAMybN0+vWOHMQxtv8+bNOBvrCsD8dj5/3lLAitlT+NTiNBdG5n0u5ly4Qt3+Uth9iBtXLmRirHe/8Y7luSixFLPtjaNkzVkopWpDOPdcHC1vYvP2vcxMDuPmtUuNDczLGP3+dDmW2zS763fwn5IuvnrbsjHdP2izab6/dzML0kL4xA2LRuUx3flceBI5D+OHnIvxY6zPxXBfwT0GDNV9pd1x20iWADcqpYqBF4BVSqnnhzm+FNiita7VWrcDbwFznHgej/SNq6ewOjuWH607RkOb7Gf0JAU1rfiYFKlRgUaH4lVykkIBpCx1BNtO1vDxv+zGYlY8evtMo8MRbsRkUnznumwqmjr56/ZCABrbu/n52/nMeuRd7v/7Po5Xuqap246CWkrq2vnEwhSXPL4QQniz4RLGNK314fOvdJSQpo30wFrrb2mtk7XWacAdwCat9d3D3OUdYIZSKtDRAOdK4NhIz+OpTCbFF1Zk0mvT7C6sMzocMYoKa9pIiQqUQdJjLDshFJOCXGl8c0GvHijlvqf3khwRwKtfXMKkuBCjQxJu5oqMKNZOi+f/bS7g0XeOs+wX7/OXrQXMnhDOroI61v5uK1/914eU1LWN6vM+v7uEyCALa3PiR/VxhRBCDJ8wDjfS4pLruZRStyilSoFFwJtKqXcAtNYNwG+AvcCHwAGt9ZuX+jyeYEZyOEEWMzsLJGH0JM52SBWjK9DiQ2ZMsCSMQ9Bas66gm/9+8RBXZETy4ucXER/mNVONxCj732uy6LHa+MP7p1iUGcX6h5bz9H0L2PqNldy/PIO3cyu4+rGtlDeOzlzUyqZONuRVc/u8ZOnoK4QQLjDcrvS9SqnPaa2fPPdKpdRngP0X8yRa683AZsefXwOGbJqjtX4e+2gNgb3z3IL0SHYU1BodihglvVYbxXVtXJUdZ3QoXiknKYyd8vs0gNWm+f4bubx8soebZyXyy9tmYvGR1W9x6dKig/j7p68gyM/MjOSz3Uojgix865psrs1J4KY/7mBPUT03z0667Of7194zWG2auxZIOaoQQrjCcAnjV4DXlFKf4GyCOA+wALe4OjBhtzgzmveP51HZ1Cnf+HuA1z8sp8eqmZkcZnQoXmlaYiivHSyjuqWT2BD5ferssfLgPw/y7rEqrk335Tcfm4XJJB1kxeVblBl1wdumJYbi72viUGnjZSeMvVYb/9xzmmWTokmNCrqsxxJCCDG0C36NrLWu0lovBn4IFDsuP9RaL9JaV45NeKLvH91dhbIq4u4a2rr56Vt5zE2N4Oppss/GCNOT7In60XJpfNPa1ctdT+7mvbwqfnjjND42xSLJohgTPmYTOYlhHCm9/PLwt3MrqWzu5O6FqaMQmRBCiKGMWHektX5fa/2447JpLIISZ01NCCU80Jcdp2Qfo7v7xfp8mjp6+MktOfLB3CBTEx2dUkfhg6q7e/1gGQdON/L4nbNldI8Yc9OTwzha3kyv1XbJj9HY3s0j646RFR/CVVmxoxidEEKIc8lGlXHOZFIsyohiV0EdWmujwxGXaF9xPS/sPcNnl6aTFR9qdDheK8Tfl/ToIHLLJWHcfrKWpPAArpueYHQowgvNSA6jo8dKQc2ld0t9ZN0x6tu6efT2mfhI12khhHAZeYd1A4szoyhr7OB0fbvRoYhL0GO18fBruSSFB/DQ6klGh+P1cpLCvH4Wo9Wm2VlQy5KJUSglq91i7E1PsjfDOVzaeEn335RfxasHyvjiikxykmRPuBBCuJIkjG5gUWY0gIzXcFN/217E8aoWfnDjNAItw/WZEmMhJzGUssYOGtq6jQ7FMLllTTR39rJkYrTRoQgvlREdRJDFzJFLGHPT3NnDt1/NZXJcMP+1aqILohNCCHEuSRjdQGZMEHGhfuw4JY1v3E1nj5U/vH+Kq7JiWTNVRmmMB32rEd5clrrd8V4iCaMwismkyEkK49Al7Cf+ybo8qls6+dVtM2XuohBCjAFJGN2AUorFmdGyj9ENbciroqWzl88sTTc6FOEwra/xjReXpe44VUt2QijRwX5GhyK82MwJ4eRVNNPd63zjm8Oljfxr3xnuX57JzAnhI99BCCHEZZOE0U0syoyirq2bE1WtRociLsIr+0tJDPNnYcaFZ5KJsRUeaGFCZAC5l1AK5wk6uq3sK25g6UR5TQpjTU8Ko7vXxomqFqfv88yOYoIsZr60MtOFkQkhhDiXJIxuYrFjHqOUpRrvVHUri3+2ccRmDdUtnWw9WcvNs5NkjMY4k5MY5rUlqftK6um22qQcVRhuRrK9PNzZfYy1rV2sO1zBbXOTCfH3dWVoQgghziEJo5tIjggkNSpQEsZx4JUDpZQ3dfKL9fnDHvfGh+VYbZpb5ySNUWTCWTlJYZTUtdPU0WN0KGNu+8lafM2KBemRRocivFxKZCCh/j4cPm8f46b8KjYcqxp0/D8/OE231cY9MjdUCCHGlCSMbmRtTjwb86slaTSQ1po3D1fg72tix6k6dhZc+Fy8cqCMmclhTIwNGcMIhTP6Gt8c9cJVxu2napmTEiEde4XhlFLMSA4fUK1RUtfGF54/wBf+sX9A2XiP1cbzH5SwbFI0mTHBRoQrhBBeSxJGN/KVqyaTERPE1148RGO7944EMFJuWTOn69t5+Nps4kL9+PW7J4ZsRHSsvJm8imY+OjfZgCjFSPoa3xz1ssY39W3dHC1vZqmUo4pxYnpyGMcrW+jssaK15nv/Poqv2UREoIUHXzhIR7cVgHeOVlLV3MW9sroohBBjThJGNxJgMfP7O2ZT19bFt149Ih1TDbDucDm+ZsWNM5P48qpJ7C9pYPOJmkHHvXawFF+z4voZiQZEKUYSHexHQpi/1+1j7FsRXzJJEkYxPsxMDqPXpsmvbOGtI5VsOVHD1z4ymd9+fBaFNW385K1jADy7s5iUyEBWTIk1OGIhhPA+kjC6mZykML72kSm8nVvJS/tLjQ7Hq2itWXe4gqUTowkL9OVj8yaQHBHAr989PiB577XaeP3DclZOiSUyyGJgxGI4M5LD2HGqjvbuXqNDGTM7TtUS4ufDDEdJrhBGm55sH42xs6CWH/7nKDlJoXxyYSpLJkbzuWXpPL/7NI9vPMne4gbuWZSKWRqICSHEmJOE0Q3dvyyDRRlR/PCNoxTXthkdjtc4VNpEWWMH1zlWDS0+Jr6yejK5Zc28c7Sy/7htp2qpaeni1jlSjjqefW5ZBrWtXTyxtdDoUMaE1pptJ2tZmBmFj1ne+sX4kBjmT1SQhcfeO0lNaxc/uXl6/+vz61dPITshlF+/d4IAXzO3z51gcLRCCOGd5FODGzKZFL/+2ExMJsVv3jthdDhe401HOeqaqXH9190yO4nMmCC+83ou1z++jXk/3sCnn9lLeKAvK7NiDIxWjGReWiTXTU/gL1sKqWzqNDocl7HaNOsOl3P949spbehgVZaU9InxQynF9OQwuq02PrkwlZkTwvtv8/Mx8/s7ZuHva+Jj85IJC5RRGkIIYQRpk+emEsMDuDYngbdyK+jutWHxkdzflfq6oy6fFENYwNkPLWaT4vs3TOPX7x4nMsjC9KQw4kL9WTIxGj8fs4ERC2d8c20W7x2r4tF3j/Po7TONDmfUvX2kgl+sz6e4rp2M6CB+8dHpskojxp1VWbGU1LXz9aunDLptUlwI27+5ivAASRaFEMIokjC6sdVT4/jXvjPsLa6XIdwudvBMI+VNnUN+oFk+OYblk2U10R2lRAVy39I0nthayL2L0/rHbXiCquZOHnzhIJkxwfzpE3P4yLR42f8lxqV7FqXxyYWpKDX06zM62G+MIxJCCHEuWZZyY0snRuPnY+K9IQYci9H15uEKLGYTq88pRxWe4UsrJxIRaOFH646htaautYsX957hKy8c5Eip+3ZRfWZnMVab5i+fnMs10xMkWRTj2oWSRSGEEMaTFUY3FmAxs2xSNBvyqvj+DVPlH1wX6bHaeOtIBcsnxxDqL2VRnibU35evrpnMd1/P5brfbye/shmbtpcbbztZy2tfXEJKVKDRYV6U1q5ent9dwtqceFKjgowORwghhBBuTFYY3dzq7DhKGzo4XtVidCge60+bC6ho6uSuK2Tvl6e6c/4E5qdFYNOaL6+axJsPLuWdryyn16a595k9NLZ3Gx3iRXlhz2laOnv53LIMo0MRQgghhJuTFUY3tyrb3vFww7EqsuJDDY7G8+RXNvP4ppNcPyOBVVlSjuqpfMwmXvr84kHXP3nPPO5+6gPuf24/z31mgVs0Muqx2nh6RzEL0iKZnRJhdDhCCCGEcHOywujmYkP8mTkhnPfyqo0OxeP0Wm38z0uHCfX35Yc3TjM6HGGABemR/Or2Gewpqud/XjqM1trokEb01pEKyho7uH+5rC4KIYQQ4vJJwugB1mTHcuhMI9UtnjtLzghPbCvkSFkTj9yUQ5R06fNaN81K4sGrJvHGoXLyKsZ36bfWmr9sKSQzJkjmLQohhBBiVEjC6AH6OnduklXGUXOquoXH3jvJNTnxXDcjwehwhMHWZNt/x0ob2g2OZHg7C+o4VtHM55ZlYJKuqEIIIYQYBZIweoApcSEkRwSwIU/Ga4yWH/7nGEF+Zh65KcfoUMQ4EB/mD0BF0/hexX9uVwnRwRZunp1kdChCCCGE8BCSMHoApRSrs+PYdrKWjm6r0eGMS71WG398/xRNHT0jHtvW1cuugjruWJBCTIiUogqICrLga1bjOmHs7rWx7WQNH5kWj7/v+G/OI4QQQgj3IAmjh1gzNY6uXhvbT9UaHcq4dKi0kV+9c5y3jlSMeOze4np6bZrFmVFjEJlwByaTIj7Mn8qmDqNDuaD9JQ20dVu5cnKM0aEIIYQQwoNIwugh5qdFEuznw/vHZR/jUIpr7XvPimrbRjx2V2EdvmbFvNRIV4cl3EhCaADl43iFccuJGnxMSr7oEEIIIcSocnnCqJQyK6UOKqXWOX6+XSl1VCllU0rNG+L4FKVUq1Lq666OzZNYfEwsmRjF5vxqt2j9P9ZK6u0JY2FN64jH7iqoY/aECAIsUtYnzrKvMI7vhHFuagQh/r5GhyKEEEIIDzIWK4wPAXnn/JwL3ApsvcDxvwXednVQnmjFlFjKmzo5WT1yUuRtSursK4uFI6wwNnf2kFvWxEJZpRHnSQi3J4w22/j7QqaquZO8imaunCLlqEIIIYQYXS5NGJVSycB1wFN912mt87TWxy9w/M1AIXDUlXF5qhWOD4ubpSx1kOI6+wrj6bp2eof5wL+nsB6bRsr6xCAJof50W23Ut3cbHcogW0/UAMj+RSGEEEKMOh8XP/5jwDeAkJEOVEoFAd8E1gAXLEdVSt0P3A8QFxfH5s2bRyVQV2ttbR2TWJODFa/tPsFk2xmXP5c7OVXZhr8ZOq2akrq2C56LF/O68DVBS/FhNp+WOXauNla/F6OhrqoXgHUbt5MWNnblygeqekkJNREdcOHv9176sJMwP0X18QPUnLi01607nQtPJudh/JBzMT7IeRg/5FyMH2N9LlyWMCqlrgeqtdb7lVIrnLjLD4Hfaq1blbrwBx6t9RPAEwDz5s3TK1Y489DG27x5M2MR63XtefxtRxHzFi0l2M/V3we4h6b2HtrWv8s1OfG8nVtJi/a/4Ln4xaFtzE/3Zc2qhWMbpJcaq9+L0RBZ2sjjB3eQOHEaK6bFj8lz7i2u5/frdxEZZOHPd89iQfrgRkxWm+ahLe+xeloCK1fOvOTncqdz4cnkPIwfci7GBzkP44eci/FjrM+FK0tSlwA3KqWKgReAVUqp54c5/grgl47jvwJ8Wyn1Xy6MzyNdOSWGHqtmp4zX6FdSb9+3uDIrFoCKtqFLUhvausmraJZyVDGk+DB/ACqbx67xzW/fO0F0sIXwQF8+8dRuXtw3uHLgUGkjTR09/SXpQgghhBCjyWUJo9b6W1rrZK11GnAHsElrffcwxy/TWqc5jn8M+KnW+g+uis9TzUuNJMhiZrNjT5M4u39xZnI4EYG+VLbZhjxud2EdAIskYRRDiA7yw9esqBijTqkfFNaxs6COz1+ZyWtfXMLCjCi+8fJhfvpWHr3Ws6/hLcdrMClYOjF6TOISQgghhHcZ8zmMSqlblFKlwCLgTaXUO2Mdgyezj9eIZsvxGhmv4VDi6IyaEhlIRkzwBRPGXYV1BFrMzEgOH8vwhJswmRRxof5UNHaMyfP9dsMJYkL8uHthKmEBvjx973zuWZTKE1sLue3PuyhwjIjZcqKGmRPCiQiyjElcQgghhPAuY5Iwaq03a62vd/z5NcfKo5/WOk5rffUQx/9Aa/3oWMTmiVZmxVLW2CHjNRyK69qJD/UnwGImPTqIyvahE+ldBXXMT4vE1zzm36MIN5EQ5j8mK4y7CurYXVjPF67MxN/X3mDHx2zikZty+P2dsymqbePa323j8Y0nOVTaKN1RhRBCCOEy8snYA8l4jYFK6tpIjQoEICMmiKYuTUtnz4Bjqlvs8yulHFUMJyEswOV7GLXW/HbDCWJD/LjripRBt984M5H3vrqcpROj+fV7J9BaxmkIIYQQwnUkYfRACWEBTIkLYfNx2ccIUFLfTlpUEAAZ0cEAFDnKVPvsLqwHYFGGJIziwvpWGF1Z7r2roI49RfV8ccXZ1cXzxYb689Sn5vGr22bwiStSpIxaCCGEEC4jCaOHWjElhr3F9TR19Ix8sAdr6+qlpqWLlHNWGAEKawYmjNtP1hDi78O0xNAxj1G4j/gwf7p7bdS3dbvsOR7fdIr4UH/uWDB4dfFcSilunzeBn9wyHbNJZoYKIYQQwjUkYfRQ189IxKbhgef20dljNTocw5Q4OqT2rTCmRgWigMJzVhh7rTbeO1bFyimx+Mj+RTGMhLAAAJftYyypa2NXYR2fXJR6wdVFIYQQQoixJJ+OPdT05DB+87GZfFBUzwPP7aer1zuTxpI6e2LYt4fRz8dMdICisOZsQ6A9RfU0tPdwTc7YDGMX7ivBMYvRVQnjawfLUApumZ3kkscXQgghhLhYkjB6sJtmJfHzW6ez5UQND/7zID3WocdJeLKSevsKY1/CCBAfZBqwh/Ht3Er8fU1cKYPPxQj6EsbKptEfraG15tUDZSzOjCIxPGDUH18IIYQQ4lJIwujhPj4/hR/cMJV3jlbxzVcOGx3OmCupayMqyEKIv2//dfFBiqLaNrTW2Gyad45WsmJyLIEWHwMjFe4gOtgPH5Oi3AUrjPtKGjhd386ts5NH/bGFEEIIIS6VJIxe4N4l6TxwZQavHiij+LzuoJ6uuLZ9wOoi2FcY27utVDZ3cuB0A9UtXVwzXcpRxchMJkVcqD+VLkgYXz1QSqDFzFopjRZCCCHEOCIJo5e4Z1EaAG8eqTA2kDFWUtfW3/CmT0KQ/WVfVNPG27msBqsuAAAUmklEQVSVWMwmVmXFGhGecEOJ4f5UjHJJamePlXWHK1ibE0+Qn6x0CyGEEGL8kITRSySFBzAnJZw3D3tPwtjZY6WiuZPU8xLGuED7CIKCmlbW51aybFL0gJJVIYYTHxYwYtObp7YV8vL+Unqd3Df83rEqWjp7+egcKUcVQgghxPgiCaMXuW5GIscqmgd0CPVkpQ3taA1p0QNLUiP8FQG+Zt44VE5ZY4eUAIqLkhDmT0VTJ1rrIW8vbWjnx2/m8fWXDnH1Y1t583AFNtvQx/Z59UApCWH+LMyIckXIQgghhBCXTBJGL3KtY5/eW15Sllpca++QmhI5MGE0KUV6dBB7ixvwMSnWTI0zIjzhphLC/OnutdHQ3jPk7ZvyqwH43vVTMSnFl/7vANc/vp3HNpzg/ePV1Ld1Dzi+uqWTrSdruWV2EmaTcnn8QgghhBAXQzbLeJGEsADmpUaw7nAF/7VqktHhuFyxYwbj+XsYAdJjgjhW0cyizCjCAy1jHZpwY32jNcobO4gMGvzaee9YFRkxQXx6aTqfWpzGvz8s48ltRfxu40n6FiXjQ/0JC/AlyM9MR48Nq01zq5SjCiGEEGIckoTRy1w3I4Ef/ucYp6pbmRgbbHQ4LlVS106ovw/hgYP3J2ZG25NIKUcVFys+zD4jsbKpk5yksAG3tXT2sLuwjvuWpANgNilunZPMrXOSae3qJbesiUNnGjlR1UprVw+tXb1YNdy5YILH/z4KIYQQwj1JwuhlrslJ4JF1x3jrSAUPXuXZq4zFdW2kRQeh1OAyv4WZUbx6sIy10yRhFBcn0bHCOFSn1K0naumxalZnDy5zDvbzYWFGlOxTFEIIIYRbkT2MXiY+zJ/5qZGsO1xudCguV1LXPmj/Yp/FmdFs/+YqooL9xjgq4e6igv3wMakhO6VuzKsiItCXOSnhBkQmhBBCCDH6JGH0QtfNSOBEVSsnqlqMDsVlNhyr4nR9O7MmyAd3MbrMJkVcqD+V5yWMvVYbm45XszIrFh+zvLUKIYQQwjPIpxovdE1OPErBukOeucpY39bN/756hOyEUO5ZlGZ0OMIDJYT5U35eSer+kgYa23tYM0Q5qhBCCCGEu5KE0QvFhvqzbFIMj79/im++fJiali6jQxo1Wmu+8/oRmjq6+c3HZmLxkZe4GH3xYYNXGDfkVWExm1g2OcagqIQQQgghRp98mvZSf7hrNp9dms4rB0pZ+ehm/rylgO5em9FhXbY3DpXz1pFKvrpmMtkJoUaHIzxUSmQgZxo6ePVAaf91G/OqWZQZRbCf9BITQgghhOeQhNFLhfr78vB1U3n3q8tZmBHJz9/O51fv5Bsd1mWpbOrku6/nMiclnAeWZxodjvBgn12WwbzUCP77xUN89/Vc8iubKaxtY3V2rNGhCSGEEEKMKkkYvVxGTDBPfWo+101P4OX9pW67yqi15tuvHaHbauPXH5uF2TR4lIYQoyUyyMI/PnsFn1uWznO7S7j9z7sAuEr2LwohhBDCw0jCKAC4fV4yDe09bMqvNjqUS/LesSo25VfztTVTSI8OMjoc4QV8zCYevm4qf7hrNlabZmZyGInhAUaHJYQQQggxqmSzjQBg2aQY4kL9eHl/KWtz3GuYfWePlUfWHWNSbDD3LkkzOhzhZa6fkcj8tEiULGoLIYQQwgPJCqMA7LPlbpmdzPvHq92ua+r/21xAaUMHj9yUg6/MvxMGiAv1JzbE3+gwhBBCCCFGnXy6Fv1um5uE1ab594dlRofitJK6Nv68pYAbZiayKDPK6HCEEEIIIYTwKJIwin4TY0OYNSGcl/aVorU2OhynPPKfY/iaFA9fm210KEIIIYQQQngcSRjFALfNTeZ4VQtHy5uNDmVE7+dXszG/modWTyI+TMoBhRBCCCGEGG2SMIoBbpiRiMXHxMv7S0c+2GDP7iomMcyfexenGx2KEEIIIYQQHkkSRjFAWKAvH5kax+sfltHVazU6nAtqaOtm+8labphlT3CFEEIIIYQQo08+aYtBbpubTGN7D+8dqzI6lAtaf7SSXpvmhhmJRocihBBCCCGEx3J5wqiUMiulDiql1jl+vl0pdVQpZfv/7d17kJX1fcfx95fl5oKKgDAoCMjFG6IgaeKFChgvUTTxEi/V2naMiU61MTaZ0TE6NRpbtSZGOzaxqW01GGw1VlGj1VCiMlEjdxAvQZAghov1AqiA8Osf56Fz3D677MK5PMu+XzNn9uzv/J7n+Z3z4Qt89znP2YgYVzbv+IiYFRELsq+Tqr025Rs/Ym8G9d6Nf525rN5LadZj81cytG8PDtlnj3ovRZIkSdpl1eIM4zeBxWXfLwTOAJ5tMm8tcGpK6VDgz4D7arA25WjoFPz5UUN5+a33mL/i/Xov5/9Zs24jv1nyLqeOHkD429IlSZKkqqlqwxgRA4FTgJ9uG0spLU4pvdZ0bkppTkppZfbtIqB7RHSr5vrUvK+OG0iPrg38SwHPMv5y4TtsTTD5MN+OKkmSJFVTVPP37UXEg8DfArsD304pTS57bEY29nLOdmcBl6SUvpjz2NeBrwP079//iKlTp1Zp9ZW1fv16evbsWe9ltMmUxRuZvvxTbjt2N3p1L87lrje9+DEfbU7ceEzjDm3fHrPYVZlFcZhFMZhDcZhFMZhDcZhFcTTNYuLEibNSSuNa2GSndK7WjiNiMrA6pTQrIia0YbtDgJuBE/IeTyndDdwNMG7cuDRhQqt3XVczZsygvax1myGjNvDMbTN4s2EgV04YWe/lALDy/Y95/cnp/PXxI5kwYcQO7aM9ZrGrMoviMItiMIfiMItiMIfiMIviqHUW1TxtdDRwWkQsA6YCkyLiZy1tkL2F9WHgwpTSkiquTa0wpG8PjjuwH1NeeItPNhfjV2w8seAdwLejSpIkSbVQtYYxpXR1SmlgSmkIcC4wPaV0QXPzI6IX8DhwdUppZrXWpbb5i6OH8u6GTUybt3L7k2tg2ryVHLrvngzt26PeS5EkSZJ2eTW/MC0iTo+IFcCRwOMR8VT20GXAcODaiJib3frVen36rKOG9eGA/rtzz8xlVPN619ZYtnYD81Z8wOTRA+q6DkmSJKmjqEnDmFKase0Db1JKD2dnHrullPqnlE7Mxm9MKfVIKR1edltdi/WpeRHBReOHsvidD7nvhbfqto5Nn27lyn+fy25dGjjtcN+OKkmSJNVCcT76UoV11tiBHHdgP2547BXmLH+vLmu4ftoiZi9/n1u/OpoBe+5WlzVIkiRJHY0No7arU6fgtrMPo/8e3fnLKbN5b8Ommh5/6kvLmfLici45dhiTR3t2UZIkSaoVG0a1Sq/Grtx1/ljWrt/EFQ/MZevW2lzPOHv5e1z3yCLGj+jLd048oCbHlCRJklRiw6hWGz2wF9edejC/fn0NP3zm9ap/CM66TzZz6c9m0X/Pbtx53hgaOkVVjydJkiTps2wY1Sbnf34/zhi7L3dO/x0X3/syq9d9UrVj3f/iclZ9uJHbzxlDr8auVTuOJEmSpHw2jGqTiODvzzqM755yEM++sZYTfvgsj82v/O9o3PjpFv75+aUcNawPRwzeq+L7lyRJkrR9Noxqs06dgq+N358n/mo8g/v04LL75/C9aa9U9BgPz36b1es2cumEYRXdryRJkqTWs2HUDhverycPXXIkFx45mHtmLq3YmcYtWxN3P/smo/bdg2OG963IPiVJkiS1nQ2jdkrnhk5cO/lgxu7Xi6sfWsDydz/a6X3+16I/8ObaDVxy7DAi/KAbSZIkqV5sGLXTujR04kfnjiECLv/5bDZ9unWH95VS4se/XsLgPo18adSACq5SkiRJUlvZMKoiBvVu5JazRjNvxQfc+tSrO7yf37z5LvNWfMA3/niYv0ZDkiRJqrPO9V6Adh0njRrAn35hMP/03FIG9W7kgs8PptN2mr5PNm9h0coPWbJmPUvWrOfpRavo27MbZ4zdt0arliRJktQcG0ZV1DWnHMTStRu47pFFTJu3kptOP5QR/XfPnTv91VVc+5+LePv9jwHo2tCJIX0bufErh9C9S0Mtly1JkiQphw2jKqp7lwbuu+iP+I9ZK7jpicWcfMdzXDx+f447qB/D996dPRu7sOrDT7h+2iKeWPAHRvTryV3nj+WQffZg4F6Nvg1VkiRJKhAbRlVcRHD2uEEcd2A/vv/4Yu6asYS7ZiwBoG/Prny8aQubtya+c+IBXDx+f7p29lJaSZIkqYhsGFU1fXp24wfnHM63jh/JG6vXsWT1BpasWc/mLYnLJw1nSN8e9V6iJEmSpBbYMKrqBvVuZFDvRiYdWO+VSJIkSWoL3wsoSZIkScplwyhJkiRJymXDKEmSJEnKZcMoSZIkScplwyhJkiRJymXDKEmSJEnKZcMoSZIkScplwyhJkiRJymXDKEmSJEnKZcMoSZIkScplwyhJkiRJymXDKEmSJEnKZcMoSZIkScoVKaV6r2GHRcQa4K16r6OV+gJr670IAWZRJGZRHGZRDOZQHGZRDOZQHGZRHE2zGJxS2rtaB2vXDWN7EhEvp5TG1XsdMosiMYviMItiMIfiMItiMIfiMIviqHUWviVVkiRJkpTLhlGSJEmSlMuGsXburvcC9H/MojjMojjMohjMoTjMohjMoTjMojhqmoXXMEqSJEmScnmGUZIkSZKUy4axGRExKCL+OyIWR8SiiPhmNt47Ip6OiDeyr3tl48dHxKyIWJB9nVS2ryOy8d9FxB0REc0cs9l5EXF2RLySreX+aj//IqlTFt+PiN9HxPom41dmOcyPiF9FxOBqPvciqXAOua9vzjGtiRx1ysKayFGpLCKiMSIej4hXs/38XQvHtC5y1CkL66KJCv/99GREzMv28+OIaGjmmNZEjjplYU3kqGQWZft8NCIWtnDMytVFSslbzg0YAIzN7u8OvA4cDNwCXJWNXwXcnN0fA+yT3R8FvF22r5eAI4EAfgl8qZlj5s4DRgBzgL2y7/vV+/XpAFl8ITvu+ibjE4HG7P6lwAP1fn3aaQ65r2/OMa2J4mRhTVQxC6ARmJjd7wo818LfT9ZFcbKwLqqUQ/b9HtnXAB4Czm3mmNZEcbKwJqqcRTZ2BnA/sLCFY1asLur+AraXG/AIcDzwGjCgLPzXcuYG8C7QLZvzatlj5wE/aeYPUu687A/T1+r9GhTlVu0smmzf7H+is2KeWe/Xo73l0IbX15ooSBZtyKxD10Slssge+xFwcc64dVGQLJrMsS6qmAPQBZgGnJOzjTVRkCyazLMmqpQF0BN4nlLDmdswVroufEtqK0TEEEp/uF8E+qeU3gHIvvbL2eRMYE5KaSOwL7Ci7LEV2VhTLc0bCYyMiJkR8UJEnLTjz6Z9q1EWrXURpZ/YdDg7mUNrWROtUKMsWqvD1gRULouI6AWcCvwqZxvrohVqlEVrddi6qEQOEfEUsBpYBzyYs4010Qo1yqK1OmxNQEWyuAG4DfiohcNUtC46b29CRxcRPSmder8ipfRh5F/yVj7/EOBm4IRtQznTUt6mLczrTOn08QRgIPBcRIxKKb2/3SewC6lhFq1ZywXAOODYHdm+PatADq0+VM6YNVGmhlm0Zi0dtiagcllERGfg58AdKaU38zbNGbMuytQwi9aspcPWRaVySCmdGBHdgSnAJODpppvm7M6aKFPDLFqzlg5bE7DzWUTE4cDwlNK3ssaz2U1zxna4LjzD2IKI6EIp1CkppV9kw6siYkD2+ABKP2nZNn8g8DBwYUppSTa8glIY2wwEVkZEQ0TMzW7fa25e2T4eSSltTiktpXT6ekQln2vR1TiL7a3li8A1wGlVOktTWBXKobl9WxNtUOMstreWDlsTUPEs7gbeSCndns21Ltqgxllsby0dti4q/fdTSukT4FHgy9ZE29Q4i+2tpcPWBFQsiyOBIyJiGaW3pY6MiBlVr4tqvC93V7hR6szvBW5vMn4rn7049Zbsfi9gHnBmzr5+S+ki4G0XnZ7czDFz5wEnAf+W3e8L/B7oU+/XaFfOomx+04u2xwBLgBH1fl3acw7Nvb6tzcuaqH0Wzc3ryDVR6SyAGyn9Z6LTdo5pXRQki7L51kWFc6B0nda2a7s6Aw8AlzVzTGuiIFmUbWNNVCGLJtsOoeUPvalYXdT9BSzqDTiG0qnb+cDc7HYy0IfStQxvZF97Z/O/C2womzuX7FOHKJ16X5gVyj8A0cwxc+dlQf8AeAVYQDOfTLWr3uqUxS2UfgKzNfv6N9n4M8Cqsv0+Wu/Xp53mkPv65hzTmihOFtZEFbOg9NPfBCwuG8/9UALrolBZWBfVy6E/pf/wzgcWAXcCnZs5pjVRnCysiSpm0WSfQ2i5YaxYXWzbUJIkSZKkz/AaRkmSJElSLhtGSZIkSVIuG0ZJkiRJUi4bRkmSJElSLhtGSZIkSVIuG0ZJkoCIOD0iUkQcWO+1SJJUFDaMkiSVnAc8D5xb74VIklQUNoySpA4vInoCRwMXkTWMETEhImZExIMR8WpETImIyB5bFhHXR8TsiFiw7axkRPSIiHsi4rcRMScivpyNN0TErdn4/Ij4Rp2eqiRJbWLDKEkSfAV4MqX0OvA/ETE2Gx8DXAEcDOxPqancZm1KaSzwj8C3s7FrgOkppc8BE4FbI6IHpUb0g2z8c8DFETG02k9KkqSdZcMoSVLp7ahTs/tTs+8BXkoprUgpbQXmAkPKtvlF9nVW2fgJwFURMReYAXQH9svGL8zGXwT6ACOq8UQkSaqkzvVegCRJ9RQRfYBJwKiISEADkIAngI1lU7fw2X83N+aMB3BmSum1JscI4PKU0lOVfwaSJFWPZxglSR3dWcC9KaXBKaUhKaVBwFLgmB3Y11PA5WXXOo4pG780Irpk4yOzt6pKklRoNoySpI7uPODhJmMPAX+yA/u6AegCzI+Ihdn3AD8FXgFmZ+M/wXf5SJLagUgp1XsNkiRJkqQC8gyjJEmSJCmXDaMkSZIkKZcNoyRJkiQplw2jJEmSJCmXDaMkSZIkKZcNoyRJkiQplw2jJEmSJCmXDaMkSZIkKdf/As9a9RBC4RVFAAAAAElFTkSuQmCC\n",
"text/plain": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(15, 6))\n",
"plt.plot(raw_data['Date'][-200:], raw_data['Concentration'][-200:])\n",
"plt.title('Concentration de CO2 au fil du Temps')\n",
"plt.xlabel('Annee')\n",
"plt.ylabel('Concentration de CO2')\n",
"plt.grid(True)\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Pour séparer les deux phénomènes, nous procederons de la facon suivante :\n",
"- On calcule les valeurs moyennes des concentraiton en CO2 par an et on obtent un graphique qui represente la contribution lente\n",
"- puis la différence entre les données initiales et les valeurs moyennes de chaque année est calculée pour étudier l'oscillation périodique"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Ajoutons une colonne a notre dataset pour les annees et affichons un appercu des donnees"
]
},
{
"cell_type": "code",
"execution_count": 83,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
"
],
"text/plain": [
" Date Concentration Annee\n",
"3361 2024-01-27 422.27 2024\n",
"3362 2024-02-03 423.97 2024\n",
"3363 2024-02-10 424.62 2024\n",
"3364 2024-02-17 423.58 2024\n",
"3365 2024-02-24 425.10 2024\n",
"3366 2024-03-02 424.50 2024\n",
"3367 2024-03-09 425.09 2024\n",
"3368 2024-03-16 425.71 2024\n",
"3369 2024-03-23 425.94 2024\n",
"3370 2024-03-30 425.02 2024"
]
},
"execution_count": 85,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"# Affichons les 10 dernieres lignes\n",
"raw_data.tail(10)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Contribution lente : on calcule la valeur moyenne par an et on stocke le resultat dans une nouvelle colonne appellee Moyenne"
]
},
{
"cell_type": "code",
"execution_count": 105,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"\n",
"
\n",
" \n",
"
\n",
"
\n",
"
Date
\n",
"
Concentration
\n",
"
Annee
\n",
"
Moyenne
\n",
"
Oscilation
\n",
"
\n",
" \n",
" \n",
"
\n",
"
0
\n",
"
1958-03-29
\n",
"
316.19
\n",
"
1958
\n",
"
315.474
\n",
"
0.716
\n",
"
\n",
"
\n",
"
1
\n",
"
1958-04-05
\n",
"
317.31
\n",
"
1958
\n",
"
315.474
\n",
"
1.836
\n",
"
\n",
"
\n",
"
2
\n",
"
1958-04-12
\n",
"
317.69
\n",
"
1958
\n",
"
315.474
\n",
"
2.216
\n",
"
\n",
"
\n",
"
3
\n",
"
1958-04-19
\n",
"
317.58
\n",
"
1958
\n",
"
315.474
\n",
"
2.106
\n",
"
\n",
"
\n",
"
4
\n",
"
1958-04-26
\n",
"
316.48
\n",
"
1958
\n",
"
315.474
\n",
"
1.006
\n",
"
\n",
"
\n",
"
5
\n",
"
1958-05-03
\n",
"
316.95
\n",
"
1958
\n",
"
315.474
\n",
"
1.476
\n",
"
\n",
"
\n",
"
6
\n",
"
1958-05-17
\n",
"
317.56
\n",
"
1958
\n",
"
315.474
\n",
"
2.086
\n",
"
\n",
"
\n",
"
7
\n",
"
1958-05-24
\n",
"
317.99
\n",
"
1958
\n",
"
315.474
\n",
"
2.516
\n",
"
\n",
"
\n",
"
8
\n",
"
1958-07-05
\n",
"
315.85
\n",
"
1958
\n",
"
315.474
\n",
"
0.376
\n",
"
\n",
"
\n",
"
9
\n",
"
1958-07-12
\n",
"
315.85
\n",
"
1958
\n",
"
315.474
\n",
"
0.376
\n",
"
\n",
" \n",
"
\n",
"
"
],
"text/plain": [
" Date Concentration Annee Moyenne Oscilation\n",
"0 1958-03-29 316.19 1958 315.474 0.716\n",
"1 1958-04-05 317.31 1958 315.474 1.836\n",
"2 1958-04-12 317.69 1958 315.474 2.216\n",
"3 1958-04-19 317.58 1958 315.474 2.106\n",
"4 1958-04-26 316.48 1958 315.474 1.006\n",
"5 1958-05-03 316.95 1958 315.474 1.476\n",
"6 1958-05-17 317.56 1958 315.474 2.086\n",
"7 1958-05-24 317.99 1958 315.474 2.516\n",
"8 1958-07-05 315.85 1958 315.474 0.376\n",
"9 1958-07-12 315.85 1958 315.474 0.376"
]
},
"execution_count": 105,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"raw_data['Moyenne'] = raw_data.groupby('Annee')['Concentration'].transform('mean')\n",
"# Affichons les 10 premières lignes\n",
"raw_data.head(10)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Realisation du graphique qui montre l'evolution de la concentration moyenne de CO2 par an en fonction du temps"
]
},
{
"cell_type": "code",
"execution_count": 97,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAp8AAAFNCAYAAABG9gSxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3XecVNX5x/HPs43ey9J7B0EQUKxgiS2WxF5iNComMZZEjRrNT40xmlgSY4z+bNFYgtiVgAVxVVSKiPTO0neXuizLsmx7fn/cSzLuj10GmNnZ8n2/XvPavefee+aZM3dnnj33nnPN3RERERERqQpJiQ5AREREROoOJZ8iIiIiUmWUfIqIiIhIlVHyKSIiIiJVRsmniIiIiFQZJZ8iIiIiUmWUfIrUUWbWxczyzSw50bHUFmY22szWJToO2Tcz+72ZbTaz7PJ/C2aWYWZXRVmP3nOR/aTkU+okM7vYzL4Ov3CyzGySmR2d6Lj2xsy6mZmbWcpB1rPKzE7cs+zua9y9sbuXHnyU/++53MxyImM2sxQz22hmmly4AmY20swmmlmumW01sxlmdkXE+uZm9kSYMBWY2bxy6+uZ2bNmttrMdpjZbDM7NTGvpvoys87ATcAAd28Xz7+FSmL4Tfj5k29mhWZWGrG8oKriEEkEJZ9S55jZr4C/AH8A0oEuwN+BsxIZ18E42MQ0TnKByMTnNGBbgmKp9sxsFDAF+BToBbQCfkbYhmaWBkwGugKjgGbALcAD4TENkAKsBY4L1/8WGG9m3aog/up4DFakK7DF3TcmKgB3/0OY8DYGfgp8tWfZ3QcmKi6RKuHueuhRZx4EX8j5wHmVbFOPIDndED7+AtQL140G1hH0mmwEsoArIvZtADwMrAa2A1OBBuG6I4AvCZKyOcDoiP0ygHuBL4AdwIdA63DdGsDDuPMJEo/Lw23/DGwFfg/0JEhetgCbgZeB5mEdLwJlwK6wjl8D3cJ6U8JtOgDvhvUtB66OiO9uYDzwzzC+BcDwStrQgTuB1yLKXgfuCD52/lO21+cE2gEFQKuIbQ8DNgGp4fJPgEUECe0HQNdyz/9TYFm4/nHAwnWXh+/LQ+G6TODUcsfIs+F7uz5s2+QKXmcD4PmwnoUEyeC6cq/vjTDuTOD6StpsKvB4JeuvJDjmGpUrvyB8T5tWsN9c4JwK1u05jh4jOF4XAydErL8ibOMdwErgmoh1own+Fm4FsoEX91J/hcdkuH4VcHMY43bgVaB+BbHuqy4HekUsPw/8fi/1nEjwd1AWttvz/P+/hQzgqgN8z6OKYy/vw9S9lA8KX/O28H04O2LdOOBR4CNgZxhzW4J/pHMJ/kYPidg+m+DvfjHB39tT/PdzrR3wfrjfFmBKNJ+neuhxoI+EB6CHHlX5AE4BSvZ8yVSwze+AaeEHeRuChPHecN3ocP/fAakEvXkFQItw/ePhl0BHIBk4kiCZ7Rh+qJ9GcMbhpHC5TbhfBrAC6BN+uWUAD4TrvvPFGJZdHsZxHUFvVwOC3rKTwudrA3wG/CVin1XAiRHL5b9wPw2/uOoDhxIkTCeE6+4GCsP4k4H7gWmVtKGHX5w5QPPwkROWecR2lT3nROBnEdv+GXgs/P1sgmS1f/j67wS+LPf8E8Ln7RLWe0pE2xUDV4ev5WcE/2TsSU7fBv4XaBQeAzOISLrKvc4HgM+BlkBnYD5hIhK+z7OA/wHSgB4ECdzJe6mnIVAKjKmkTccBL+ylPCU8FvZWb3r4vvWroM49x9EvCY7nCwiSwJbh+tMJkj4j6E0tAIaV+1v4I8Ex12Av9UdzTM4gSNJbEiRYP60g1n3VFXXSF8YemTB2I/rks8L3fH/jKPc+TC1X1pTgH6BLwuN0BEHS2CvieMgGhhD8/U8Nj68Lwu0fBCZF1JcNzA7bug0wE7gz4m/r0fBYSgOOPdjPWj30qOyR8AD00KMqH+EHefY+tlkBnBaxfDKwKvx9NEGvSWQiuJGgVzMpXDdkL3XeSrmeIYLeuh+Hv2fs+SIIl38OvB/+/p0vxrDscmDNPl7H2cDsiOVVVJB8hl+ipUCTiPX3A8+Hv98NTI5YNwDYVclzO0Gy8AxwDUEv5NNhmYfb7Os5LwC+CH9PDr88R4bLk4ArI/ZLIkiMukY8/9ER68cDt0W03fKIdQ3D7dsRJGu7iUikgIuATyp4nSsJk9pweSz/TT4PL/8eAbcD/9hLPR3DGPaaJIbbTCb8h2Qv67KBS8qVpYb7/G8ldV5OROIdls0AflTB9m8DN0T8LRRRQU/lfhyTl0Ys/wl48gDrqqrks8L3fH/jKPc+lE8+fwx8VK7sBeDW8PdxhP+Mhcu3lGuPEUR81oXHyOURyz8EFkS0+2tAj2jfSz30OJhHTbpGRyQWtgCtzSzF3Usq2KYDwWnzPVaHZf+po9y+BUBjoDVBD96KvdTZFTjPzM6IKEsFPolYzt5LnZVZG7lgZm2BvwLHAE0IErJor7HsAGx19x0RZauB4ZXEV38f7QjBafr7CXrObt3P53wHeNLMehD0CG939xnhuq7Ao2b2cMS+RpDE7XnvKmvP/6xz9wIzI1zfkuB9yQrLIGjH77R1udcQuS7yuOkKdDCz3IiyZIJes/K2EZwGbk9wWnRvNofrvyO81rJ1uH5PWRLBpRZFwC8qqG+P9e7u5V5Dh7CeU4G7CNo/iSBRnxex7SZ3L6yo4iiPyfLvUwf24iCP71iq7D2Ppa7AseWOnxS++5pzIn7ftZfl8p8h5ePe09b3EZzN+cTMioG/u/sjBxG7SKU04Ejqmq8ITkOeXck2Gwg++PfoEpbty+aw7p57WbeWoOezecSjkbs/EEW9HmX5/WHZYHdvClxKkJDtqx4IXl9LM2sSUdaF4JrHg/E5QcKUTnBaMOrnDJOa8QS91T8iSKb2WEtwKjyyPRu4+5cHGe9agp7P1hH1NvWKB4BkEfTgRsYfWVdmuRibuPtp5Stx9wKCY/OcSmKbDJxqZo3KlZ8TxjwNwIKs+VmCNj/H3YsrqROgo0Vk2uFr2GBm9QiuV30ISHf35gSXQkR7TMG+j8n9sa+6CgiS4z3aHeDz7Etl73ks41gLfFju+Gns7jceYH3w/+PeAODu2939BnfvSnA83WlmRx3E84hUSsmn1Cnuvp3gGrzHzexsM2toZqlmdqqZ/Snc7F8EH75tzKx1uP1LUdRdBjwHPGJmHcws2cxGhV/iLwFnmNnJYXn9cH7ATlGEvYmgV6zHPrZrQjCAItfMOhKchouUU1Ed7r6W4NrW+8PYBhMMcHk5ivgqFPaonQGcWa53Ldrn/CfBKckz+e578CRwu5kNBDCzZmZ23sHEGsaURTDY62Eza2pmSWbW08yOq2CX8WEcLcL38rqIdTOAPDO71cwahO/7IDMbUUFdvwYuN7NbzKxV+LqGmNm4cP2LBAN8Xgun30o1s5MJegPvDo9tgCcIroU9w913RfGy2wLXh/WdF+47keDav3oEx19J2Av6vSjqi7SvYzKWdX0LXBy28ykE16jGQ2XveSzjeBsYamYXhO9NmpkdYWZ9DiL2682sffi5dhvBAC/M7Ewz6x7+E7Kd4HKYKpt2SuoeJZ9S54Snk35FMEhlE0EPwy8IPuwhGN38NcEI3HnAN2FZNG4O95lJMDjgj0BSmGidBfwm4jlvIYq/wbBX7D7gi3D+xyMq2PQeYBjBl8e/gTfLrb+fIKnONbOb97L/RQTXvm0A3gLucveP9hVfFPEvcPeK5i2s9Dnd/QuCxPsbd18VUf4WQduOM7M8gkEfsZrP8jKCxGshwSnO19nL6e7QPQSnLzMJktb/9M56MGfkGQQDqTIJesafIRhN//+EvbbHh4+VZrZnRPLEcP1ugpHaa4HpQB7wCHCHuz8IYGZdCa6xPRTIjpg38pJKXu90oHcY333Aue6+Jbwc4nqCZGsbcDHBzAT7Y1/HZCzruoGgvXMJesvfJj4qfM9jGYe7byO43vwKgt7WDQSfQ6kHFHVgHMGlPssIPqf2/MPdn+A61x0EA7kecvdpB/E8IpWycp0RIiLViplNAV5x92cSHUttY2aXEwysqZY3WJDYMbNsgn8syl/+IlLlNOBIRKqt8BT1MGrwDQBEROS7dNpdRKolM3uBYJDNjeVGxIuISA2m0+4iIiIiUmXU8ykiIiIiVUbJp4iIiIhUmRo94Kh169berVu3RIexX3bu3EmjRuXniJaDoTaNLbVn7KlNY0vtGXtq09iqq+05a9asze7eZl/b1ejks1u3bnz99deJDmO/ZGRkMHr06ESHUauoTWNL7Rl7atPYUnvGnto0tupqe5pZVLeb1Wl3EREREakySj5FREREpMoo+RQRERGRKqPkU0RERESqjJJPEREREakySj5FREREpMoo+RQRERGRKqPkU0RERESqjJJPEREREakySj5FREREapkNubt4deaaRIexVzX69poiIiIi8l1Lc3Zw2bMz2FlUwon902nVuF6iQ/oO9XyKiIiI1BIzV23l3Ce+pMyd8deMqnaJJ6jnU0RERKRW+GBBNtf/azYdWzTgnz8ZSacWDRMd0l4p+RQRERGp4V6ZvoY7357H4E7Nee7yEbRslJbokCqk5FNERESkhnJ3Hv14GX+ZvIzRfdvw90uG0TCteqd31Ts6EREREdmr0jLnt+/M55XpazhnWCceOOcQUpOr/3CeuEdoZslmNtvMJoTLD5rZYjOba2ZvmVnziG1vN7PlZrbEzE6Od2wiIiIiNVFhcSk/f3kWr0xfw89G9+Sh8wbXiMQTqma0+w3Aoojlj4BB7j4YWArcDmBmA4ALgYHAKcDfzSy5CuITERERqTE25O7i0mem8+HCHO46YwC3ntIPM0t0WFGLa/JpZp2A04Fn9pS5+4fuXhIuTgM6hb+fBYxz993ungksB0bGMz4RERGRmuTdORs45S+fsTArj8cuGsoVR3VPdEj7Ld7XfP4F+DXQpIL1PwFeDX/vSJCM7rEuLBMRERGp0/IKi7nrnQW8NXs9Q7s05y8XHErXVo0SHdYBiVvyaWbfBza6+ywzG72X9XcAJcDLe4r2Uo3vZb+xwFiA9PR0MjIyYhVylcjPz69xMVd3atPYUnvGnto0ttSesac2ja1Yt+eSraU8NXc323Y7Z/dK5YweRWTOm0lmzJ6hasWz5/Mo4EwzOw2oDzQ1s5fc/VIz+zHwfeAEd9+TYK4DOkfs3wnYUL5Sd38KeApg+PDhPnr06Di+hNjLyMigpsVc3alNY0vtGXtq09hSe8ae2jS2YtWeRSVl/GXyUp6YuYIuLRvy9E8OZViXFgcfYILFLfl099v572Ci0cDNYeJ5CnArcJy7F0Ts8i7wipk9AnQAegMz4hWfiIiISHW1fGM+N746m/nr87hgeGd+e8YAGterHTNkJuJV/A2oB3wUjsya5u4/dfcFZjYeWEhwOv5ady9NQHwiIiIiCVFYXMo/vljFox8vpUFqMk9eehinDGqX6LBiqkqST3fPADLC33tVst19wH1VEZOIiIhIdeHuvD8/mz9MWsTarbs4aUA6vz97EOlN6yc6tJirHf23IiIiIjXU/PXb+d2EhczI3Erf9Ca8dOXhHN27daLDihslnyIiIiIJsHFHIQ99sITXZq2jRcM07vvBIC4Y3pmUGnKnogOl5FNERESkChUWl/LcF5k8PmU5RaVlXH1MD64d04tmDVITHVqVUPIpIiIiUgXcnUnzs/nDxEWs2xZc13nHaf3p1rpmThZ/oJR8ioiIiMRZ5HWd/do14eWrDueoXrX3us7KKPkUERERiZONeYU89GFwXWfLhmn84QeHcMGIziQn7e3GjnWDkk8RERGRGCssLuXZqZn8/ZPgus6xx/Tg2uN70bR+3biuszJKPkVERERixN2ZkV3CnY98yrptu/jegHR+Uwev66yMkk8RERGRGNiYV8htb85jyuLd9G/flFeuHsyRPevmdZ2VUfIpIiIicpD+PTeLO96ex66iUi7ul8a9lx1dp6/rrIySTxEREZEDlFtQxP+8s4B352xgSOfmPHzeENYt/FqJZyWUfIqIiIgcgIwlG7n1jblsyS/ippP68LPRPUlJTmLdwkRHVr0p+RQRERHZDzt3l3DfxEW8Mn0NfdIb8+yPRzCoY7NEh1VjKPkUERERidKs1Vv55atzWLutgGuO7cEvT+pD/dTkRIdVoyj5FBEREdkHd+fZqZncP2kxHZrX59WxoxjZvWWiw6qRlHyKiIiIVCJ/dwm/fn0OE+dlc8rAdjx43mCaaLL4A6bkU0RERKQCy3J2cM1Ls1i9pYA7TuvPVcd0x0wj2Q+Gkk8RERGRvXjn2/Xc/uY8Gqal8PJVh3NEj1aJDqlWUPIpIiIiEqGopIw/TFzE81+uYnjXFjx+yTDSm9ZPdFi1hpJPERERkVDW9l1c+/I3fLMml58c1Z3bT+tHanJSosOqVZR8ioiISJ1XWua88c06/jhpMbuKS/nbxUP5/uAOiQ6rVlLyKSIiInWWu5OxZBMPTFrMkpwdDOnUjIfPH0Kvtk0SHVqtpeRTRERE6qQ5a3O5f9Iipq3cStdWDXn84mGcdkg7jWaPMyWfIiIiUqes3rKTBz9YwoS5WbRqlMY9Zw7kopFdSEvRtZ1VQcmniIiI1Alb8nfz2JTlvDx9NSlJSVx/fC+uPraHJoyvYko+RUREpFYrK3Ne/XotD0xaTP7uEs4f3plfntibtpo+KSGUfIqIiEittTg7jzvems+s1dsY2b0lvz97EH3SNZgokZR8ioiISK1TUFTCo5OX8czUTJo1SOWh84ZwzrCOGkxUDSj5FBERkVpl8sIc7np3Aetzd3HhiM7ceko/WjRKS3RYElLyKSIiIrXChtxd3PPeAj5YkEOf9Ma89tNRjOjWMtFhSTlKPkVERKTG+2TJRn7x8jeUunPbqf248ujuui1mNaXkU0RERGq012et49Y35tKvXROevPQwOrdsmOiQpBIV/ktgZp3NbJyZfW5mvzGz1Ih1b1dNeCIiIiJ75+48kbGCm1+bw6gerXj1mlFKPGuAyvqjnwMygOuA9sCnZtYqXNc1znGJiIiIVKiszPndhIX88f3FnDmkA89dPoLG9XRCtyao7F1q4+5Phr9fZ2aXAp+Z2ZmAxz80ERERkf9vd0kpvxo/h3/PzeLKo7tzx2n9SUrSFEo1RWXJZ6qZ1Xf3QgB3f8nMsoEPgEZVEp2IiIhIhB2FxVzz4iy+XLGFO07rz9XH9kh0SLKfKjvt/gxweGSBu08GzgPmxzMoERERkfI25hVywf9OY0bmVv58wRAlnjVUhcmnu//Z3T/dS/lsdz8p2icws2Qzm21mE8Lllmb2kZktC3+2iNj2djNbbmZLzOzk/X0xIiIiUjsty9nBD5/4klVbdvLs5SP4wdBOiQ5JDlClE2CZ2alm9pmZbTazTWb2qZmdtp/PcQOwKGL5NuBjd+8NfBwuY2YDgAuBgcApwN/NLHk/n0tERERqkbzCYv4wcRGn/fVzdhWV8q+rj+C4Pm0SHZYchAqv+TSzq4FrgF8DX4fFw4EHzKyTuz+1r8rNrBNwOnAf8Kuw+CxgdPj7CwQj6m8Ny8e5+24g08yWAyOBr/bvJYmIiEhNV1rmvD5rLQ9+sIQtO4s4/7DO3HxyX9o0qZfo0OQgmfveB66b2ULgaHffWq68FTDV3fvvs3Kz14H7gSbAze7+fTPLdffmEdtsc/cWZvY3YJq7vxSWPwtMcvfXy9U5FhgLkJ6efti4ceP24+UmXn5+Po0bN050GLWK2jS21J6xpzaNLbVn7FW3Nl26rZSXFxWxOq+M3s2TuKR/Gt2a1ZyTodWtPavKmDFjZrn78H1tV9lodyufeAK4+xazfU9nYGbfBza6+ywzG73PHWBvlf6/zDjscX0KYPjw4T56dDRVVx8ZGRnUtJirO7VpbKk9Y09tGltqz9irLm26IXcXD0xazLtzNtC+WX0evbAfZw7pQDR5R3VSXdqzuqos+cwzsyHuPiey0MyGADuiqPso4MzwGtH6QFMzewnIMbP27p5lZu2BjeH264DOEft3AjZE+0JERESkZiooKuHpzzJ54tPluMP1J/Tmp8f1oGGaJo2vjSp7V28C3jWzfwCzCHohRwA/Bi7dV8XufjtwO0DY83mzu19qZg+GdTwQ/nwn3OVd4BUzewToAPQGZhzAaxIREZEaoLTMeWPWOh7+aAk5ebs57ZB23H5qf90is5arMPl096lmNhK4Fric4LT4AuAId88+iOd8ABhvZlcCawjmDcXdF5jZeGAhUAJc6+6lB/E8IiIiUk19tnQTf5i4iMXZOzi0c3Mev3gYw7u1THRYUgUqG+3eBmjl7v9TrnygmZW6+6Zon8TdMwhGtePuW4ATKtjuPoKR8SIiIlILLcrK4/5Ji/ls6SY6t2zA4xcP47RD2tW46zrlwFV22v0x4Im9lHcC7gAujktEIiIiUuvk5BXy8IdLeG3WOprWT+XO0/vzo1FdqZdSc0axS2xUlnweUsEdjj4ws4fjGJOIiIjUEqVlzjOfr+Qvk5dRWuZcdXR3fjGmN80apiY6NEmQypLPyo4KHTEiIiJSqTVbCrjptW+ZuWob3xuQzp2nD6BLKw0mqusqSz6Xmdlp7j4xstDMTgVWxjcsERERqancnXEz13LvhIUkJxl/vmAIZx/aUdd1ClB58vlLYIKZnU8w1RIEt9ccBXw/3oGJiIhIzbNxRyG3vTGPKYs3cmTPVjx03hA6NG+Q6LCkGqlsqqWlZnYIwcCiQWHxp8A17l5YFcGJiIhIzTFxXhZ3vDWPgqJS7jpjAD8e1Y2kJPV2yndVeusAd98N/KOKYhEREZEaaPuuYu5+dwFvzV7P4E7NeOT8IfRq2yTRYUk1pftWiYiIyAH7eFEOd749n407dnPjib25dkwvUpOTEh2WVGNKPkVERGS/bdqxm3veW8CEuVn0TW/Ck5cexpDOzRMdltQAUSWfZtYA6OLuS+Icj4iIiFRj7s7rs9bx+38vYldRKTed1IdrjutJWop6OyU6+0w+zewM4CEgDehuZocCv3P3M+MdnIiIiFQfa7YU8Ju35jF1+WZGdGvB/T8cTK+2jRMdltQw0fR83g2M5L/3Zv/WzLrFLSIRERGpVkpKy/jHF6t4+KMlpCQlce/Zg7hkZBeNZJcDEk3yWeLu2zUxrIiISN0zc9VW7p2wkLnrtnNi/7bce/Yg2jfTvJ1y4KJJPueb2cVAspn1Bq4HvoxvWCIiIpIo7s4Xy7fw2JRlTM/cSuvG9fjbxUM5/ZD2ukuRHLRoks/rgDuA3cC/gA+Ae+MZlIiIiFQ9d+eTJRt5bMpyZq/JJb1pPf7n+wO4aGQXGqQlJzo8qSX2mXy6ewFB8nlH/MMRERGRqlbmzqR5WTw2ZTkLs/Lo1KIB9/1gEOce1ol6KUo6JbYqTD7N7D3AK1qv0e4iIiI1m7szcV42932xiw3539C9dSMePHcwZw/tqIniJW4q6/l8KPz5Q6Ad8FK4fBGwKo4xiYiISJxtyd/NHW/N5/0F2XRsbPz1ouCazmSNYJc4qzD5dPdPAczsXnc/NmLVe2b2WdwjExERkbj4cEE2v3lrHnm7Srjt1H70KVvD8UM6JDosqSOi6VNvY2Y99iyYWXegTfxCEhERkXjIKyzmpvFzGPviLNo2qc+71x3FT4/rSZJGsEsVima0+y+BDDNbGS53A8bGLSIRERGJuS+Wb+aW1+aQnVfIdcf34rrje+uWmJIQ0Yx2fz+c37NfWLTY3XfHNywRERGJhV1FpTwwaREvfLWaHm0a8cbPjmRolxaJDkvqsGh6PgmTzTlxjkVERERixN3JWLKJ301YSObmnVxxVDd+fXI/zdcpCRdV8ikiIiI1x7SVW3jwgyXMWr2NLi0b8spVh3Nkr9aJDksEUPIpIiJSa8xdl8uDHyzh82WbSW9aj/t+MIjzh3fWnJ1Srewz+bTgJq6XAD3c/Xdm1gVo5+4z4h6diIiI7NOynB08/OFS3l+QTYuGqdxxWn9+NKor9VN1il2qn2h6Pv8OlAHHA78DdgBvACPiGJeIiIjsw5otBfzl46W8NXs9jdJSuPHE3lx5dHea1E9NdGgiFYom+Tzc3YeZ2WwAd99mZmlxjktEREQqsHJTPo9/soK3v11PSpIx9pge/PS4nrRopK9nqf6iST6LzSyZ8D7vZtaGoCdUREREqtDSnB38bcpyJszdQGpyEpeN6spPj+tJetP6iQ5NJGrRJJ9/Bd4C2prZfcC5wJ1xjUpERET+Y+GGPP72yTImzc+mQWoyVx/Tg6uO6UGbJvUSHZrIfotmkvmXzWwWcAJgwNnuvijukYmIiNRxc9fl8tePlzN5UQ5N6qVw7ehe/OTo7rTU6XWpwaKdamkZkLdnezPr4u5r4haViIhIHTZr9TYem7KMjCWbaFo/GEh0xZHdadZQA4mk5otmqqXrgLuAHKCUoPfTgcHxDU1ERKRumbZyC49NWcYXy7fQslEat5zcl8tGddXodalVoun5vAHo6+5b4h2MiIhIXePuTF2+mcc+Xs6MVVtp3bged5zWn0uO6ELDNN0LRmqfaI7qtcD2eAciIiJSl+y59/qjHy/j27W5tGtan7vPGMCFI7tocnip1aJJPlcCGWb2b2D3nkJ3fyRuUYmIiNRi89dv53/emc83a3Lp2LwBvz97EOcN70S9FCWdUvtFk3yuCR9p4SMqZlYf+AyoFz7P6+5+l5kdCjwJ1AdKgJ/vuVWnmd0OXElwben17v7BfrwWERGRam17QTEPfbiEl6evpkXDNO7/4SGce1gn3Xtd6pRoplq65wDr3g0c7+75ZpYKTDWzSQS36LzH3SeZ2WnAn4DRZjYAuBAYCHQAJptZH3cvPcDnFxERqRbKypzXZ63jgfcXk1tQxGWjuvHLk/rQrIEGEkndE81o9z7AzUC3yO3d/fjK9nN3B/LDxdTw4eGjaVjeDNgQ/n4WMM7ddwOZZrYcGAl8FeVrERERqXbmr9/Ob9+Zz+w1uQzv2oJ7zhrJwA7NEh2WSMJEc9r9NYLT5M8QnA6PWnhbzllAL+Bxd59uZjcCH5jZQ0AScGS4eUdgWsTu68IyERGRGie3oChnK0XrAAAgAElEQVQ8xb6GVo3SePi8IfxwWEfMLNGhiSSUBR2UlWxgNsvdDzuoJzFrTnCLzuuAscCn7v6GmZ0PjHX3E83sceArd38p3OdZYKK7v1GurrFhHaSnpx82bty4gwmtyuXn59O4ceNEh1GrqE1jS+0Ze2rT2Kru7VnmzufrS3h9SRH5xXBi1xTO7pVGo9Tqm3RW9zataepqe44ZM2aWuw/f13bR9Hy+Z2Y/J0geI0e7b402GHfPNbMM4BTgxwRzh0LQq/pM+Ps6oHPEbp347yn5yLqeAp4CGD58uI8ePTraMKqFjIwMalrM1Z3aNLbUnrGnNo2t6tye89YFp9i/XVvA8K4t+N1ZgxjQoem+d0yw6tymNZHas3LRJJ8/Dn/eElHmQI/KdjKzNkBxmHg2AE4E/kiQUB4HZADHE9y6E+Bd4BUze4RgwFFvYEZ0L0NERCRxcguKePCDJbwyYw2tGtXTKXaRSkQz2r37AdbdHnghvO4zCRjv7hPMLBd41MxSgELCU+juvsDMxgMLCaZgulYj3UVEpDorK3PGf72WP76/mO27irn8yGAUe1PdDlOkQtGMdm8I/Aro4u5jzaw3we02J1S2n7vPBYbupXwqsNdrSN39PuC+aAIXERFJpP+eYs9lRLfgFHv/9tX/FLtIokVz2v0fBCPW94xKX0dwrWalyaeIiEhttCgrj6c/W8lb366nVaN6PHL+EH4wVKfYRaIVTfLZ090vMLOLANx9l+kvTERE6hB3Z+ryzTz12Uo+X7aZhmnJXHV0d647obdOsYvsp2iSz6JwwJADmFlPIka9i4iI1FbFpWVMmLuBpz7LZFFWHm2a1OOWk/ty6eFdadZQSafIgYgm+bwbeB/obGYvA0cBl8cxJhERkYTaUVjMuBlree6LTLK2F9K7bWP+dO5gzjq0A/VSkhMdnkiNFs1o9w/NbBZwBGDADe6+Oe6RiYiIJMC7czbw27fns31XMaN6tOIPPziE4/q0ISlJV5yJxEI0o91fB54DJrl7WfxDEhERqXrbdxVz1zvzefvbDQzt0px7zhzI4E7NEx2WSK0TzWn3J4ErgMfM7DXgeXdfHN+wREREqs5XK7Zw0/hvydmxm1+d1Iefj+5JSnJSosMSqZWiOe0+GZhsZs2Ai4CPzGwt8DTwkrsXxzlGERGRuNhdUsojHy7lqc9X0q1VI9742ZEc2lm9nSLxFE3PJ2bWCrgU+BEwG3gZOJrg1puj4xWciIhIvCzJ3sGNr37Loqw8Lj68C3ee3p+GaVF9LYrIQYjmms83gX7Ai8AZ7p4VrnrVzL6OZ3AiIiKxVlbm/OPLVfzx/cU0rZ/CM5cN58QB6YkOS6TOiOZfvL+5+5S9rXD34TGOR0REJG6+XrWV+yYuYvaaXE7s35YHzhlM68b1Eh2WSJ0STfL5uZldDxwbLn8KPKlrPUVEpKbI3LyTP05azPsLsklvWo+HzhvCOcN0S0yRRIgm+XwCSAX+Hi7/KCy7Kl5BiYiIxMLWnUX89eNlvDRtNWkpSdx0Uh+uPKa7ru0USaBo/vpGuPuQiOUpZjYnXgGJiIgcrMLiUp7/chWPT1nOzqISLhzZhRtP7E3bJvUTHZpInRdN8llqZj3dfQWAmfUASuMbloiIyP4rKinjvTkbeOSjpazP3cUJ/dpy26n96J3eJNGhiUgomuTzFuATM1tJcHvNrgSTzouIiFQL2dsLeWX6al6ZsZbN+bsZ2KEpD547mCN7tU50aCJSTjSTzH9sZr2BvgTJ52J33x33yERERCrh7kxbuZV/frWKDxfmUObOmL5tuWxUV47trXuxi1RXFSafZvbDClb1NDPc/c04xSQiIlKhXSXOi9NW8+JXq1iak0/zhqlcdXR3Ljm8K11aNUx0eCKyD5X1fL4OfBs+IOj13MMBJZ8iIlIlCotL+XzZZt6fn82/5xRQWDqfQR2b8qdzB3PmkA7UT01OdIgiEqXKks9zgAuAwcA7wL/cfXmVRCUiInXejsJiPlmyiQ/mZ/PJko0UFJXStH4Kh6WncNPZIxnaubnm6RSpgSpMPt39LeAtM2sEnAU8HN7j/Q53/7SqAhQRkbpj284iPlqYw/sLspm6bDNFpWW0blyPHwztyCmD2nFEj1Z88flnDOvSItGhisgBima0eyGwHcgDugCaJE1ERGJqxaZ8nvk8kze+WUdRSRmdWjTgslFdOWVQO4Z2aUGyBg+J1BqVDTgaA1wEjAQmA4+6+9dVFZiIiNRu7s70zK088/lKJi/aSFpKEucM68Qlh3dhYIemOqUuUktV1vP5MTAXmArUAy4zs8v2rHT36+Mcm4iI1EIlpWVMmp/N05+vZO667bRslMYNJ/TmR6O60rpxvUSHJyJxVlnyqYnkRUQkZnaXlPLStDU8NzWT9bm76N66Eff9YBDnDOuk0eoidUhlA45eqMpARESk9lq7tYCfv/wN89ZvZ2S3ltx95kBO6NdWE8GL1EHRDDgSERE5YFMW5/DLV+dQ5s7//ugwTh7YLtEhiUgCKfkUEZG4KC1zHvloCY9/soIB7ZvyxKXD6NqqUaLDEpEEU/IpIiIxt2nHbm4YN5svV2zhguGdueesgbquU0SAKJJPM2sDXA10i9ze3X8Sv7BERKSmmrlqK7945RtyC4r507mDOX9450SHJCLVSDQ9n+8AnxPM9Vka33BERKSmcneenZrJ/ZMW06lFA978+QgGdmiW6LBEpJqJJvls6O63xj0SERGpkTbk7uLjRTm8NzeLGZlb+d6AdB46fwhN66cmOjQRqYaiST4nmNlp7j4x7tGIiEi15+4s2JDHRwtzmLwohwUb8gDo3roRd50xgMuP7Ka7E4lIhaJJPm8AfmNmRUBxWObu3jR+YYmISHVSVuZ8sWIzHy4IEs6s7YWYwWFdWnDbqf04aUA6Pds0TnSYIlID7DP5dPcmVRGIiIhUP8WlZbw3ZwN/z1jB8o35NEhN5pjerfnlSX04vl9b3Q5TRPZbVFMtmdmZwLHhYoa7T4hfSCIikmiFxaW8Nmsd//vpCtZt20W/dk149MJDOXlgO02ZJCIHJZqplh4ARgAvh0U3mNnR7n7bPvarD3wG1Auf53V3vytcdx3wC6AE+Le7/zosvx24kmBU/fXu/sEBvSoRETkg+btLeGX6ap7+PJNNO3YztEtz7jlzIMf3a6vrOEUkJqLp+TwNONTdywDM7AVgNlBp8gnsBo5393wzSwWmmtkkoAFwFjDY3XebWduw3gHAhcBAoAMw2cz6uLumdxIRibPcgiL+8cUqnv9yFdt3FXN0r9Y8euGhjOrRSkmniMRUtHc4ag5sDX+PatI2d3cgP1xMDR8O/Ax4wN13h9ttDLc5CxgXlmea2XJgJPBVlDGKiMh+Kiwu5bkvMnkiYwU7Cks4aUA6147pxaGdmyc6NBGppaJJPu8HZpvZJ4ARXPt5ezSVm1kyMAvoBTzu7tPNrA9wjJndBxQCN7v7TKAjMC1i93VhmYiIxFhpmfPGrHU88tFSsvMKOaFfW24+uS/922siExGJLws6KPexkVl7gus+DZju7tn79SRmzYG3gOuAccAUgimcRgCvAj2AvwFfuftL4T7PAhPd/Y1ydY0FxgKkp6cfNm7cuP0JJeHy8/Np3FjTkcSS2jS21J6xV53a1N2Zs6mU15YWsT7f6dEsifP7ptGvZc0ZRFSd2rO2UJvGVl1tzzFjxsxy9+H72q7Cnk8z6+fui81sWFi0LvzZwcw6uPs30Qbj7rlmlgGcEtbzZnhafoaZlQGtw/LIGwB3Ajbspa6ngKcAhg8f7qNHj442jGohIyODmhZzdac2jS21Z+xVlzadvWYb909azIzMrXRv3Yi/n9WXUwe1q3HXdFaX9qxN1KaxpfasXGWn3X9F0MP48F7WOXB8ZRWbWRugOEw8GwAnAn8kuA70eCAjPAWfBmwG3gVeMbNHCAYc9QZm7N/LERGR8pbl7ODPk5cycV42rRunce9ZA7lwZBdSk5MSHZqI1EEVJp/uPjb89VR3L4xcF06jtC/tgRfC6z6TgPHuPsHM0oDnzGw+UAT8OOwFXWBm44GFBFMwXauR7iIiB275xh08+vFyJszdQMPUZG44oTdXH9uDxvWiHWsqIhJ70XwCfQkMi6LsO9x9LjB0L+VFwKUV7HMfcF8UMYmISAWWb8znsSnLeHfOBhqkJvPT43py9TE9aNkoLdGhiYhUes1nO4LR5g3MbCjBYCOApkDDKohNRET2w8pN+fz14yDprJeSzNhjezD2mB600i0wRaQaqazn82TgcoKBP49ElO8AfhPHmEREJErbdhaxKCuP12et4+1v15OWksRVx/Rg7LE9dN91EamWKrvm8wWCazbPKT/dkYiIVK2S0jIyN+9kUfYOFmXlsSgrj8VZO8jOCy7Jr5+axJVHd2fssT1p00RJp4hUX/u85tPd3zCz0wlue1k/ovx38QxMRKSuKiktY9nGfOauy2Xuuu3MW7+dJdk72F1SBkBKktGrbWNG9WxF//ZN6NeuKYM7NaN5Q13TKSLV3z6TTzN7kuAazzHAM8C5aAokEZGYcHcyN+9k7rrtzAmTzQUbtlNYHCSaTeqlcEinZlw2qiv92jWlf/um9GrbmLQUTZMkIjVTNKPdj3T3wWY2193vMbOHgTfjHZiISG2WtX0Xb36znte+XsuqLQVAcOp8UIdmXDyyK4M7NWNwp2Z0a9WIpKSaNQm8iEhlokk+98zxWWBmHYAtQPf4hSQiUjvtLill8sKNjP96LZ8v20SZw8juLbn62B4M69KC3m0bk6KJ30Wklosm+XwvvDf7g8A3BHc3ejquUYmI1CKr80q5+90FvP3tenILimnfrD7XjunFuYd1omurRokOT0SkSlWafJpZEvCxu+cCb5jZBKC+u2+vkuhERGqotVsLeG/uBt79dgOLswtJS1nD9wakc/7wzhzVqzXJOpUuInVUpcmnu5eF13iOCpd3A7urIjARkZome3shE+ZuYMLcLL5dmwvA0C7NubR/Gjefd5xGo4uIEN1p9w/N7BzgzfAe7CIiEtqcv5tJ87J4b24WM1dtxR0GdmjKbaf24/RD2tO5ZUMyMjKUeIqIhKJJPn8FNAJKzKyQ4Dab7u5N4xqZiEg1VVJaRsaSTYybuYZPlmyitMzp1bYxvzyxD98f3J4ebRonOkQRkWormknmm1RFICIi1d363F28OnMt42euJTuvkDZN6jH22B6cdWgH+qY3wUzXcYqI7Es0k8x/7O4n7KtMRKQ2KiktY8rijfxrxhoylm4C4Njebbj7zIGc0L8tqZoaSURkv1SYfJpZfYI7G7U2sxYEp9sBmgIdqiA2EZGE2V1SytOfreTFaavJydtN2yb1+MWYXpw/vDOdWzZMdHgiIjVWZT2f1wA3EiSas/hv8pkHPB7nuEREEmbmqq3c9sZcVmzayXF92nDvWV04vl9bTQAvIhIDFSaf7v4o8KiZXefuj1VhTCIiCbGjsJg/vb+EF6etpmPzBrzwk5Ec16dNosMSEalVohlw9JiZHQl0i9ze3f8Zx7hERKrURwtz+O3b88nZUchPjurOTd/rQ6N60UwIIiIi+yOaAUcvAj2Bb4HSsNgBJZ8iUuNt2rGbu99bwL/nZtE3vQlPXDqMoV1aJDosEZFaK5p/64cDAzTBvIjUJmVlzuvfrOO+fy9iV1EpN53Uh2uO60laiq7rFBGJp2iSz/lAOyArzrGIiMRdcWkZ783ZwJOfrmBpTj4jurXg/h8OpldbTQwvIlIVokk+WwMLzWwGEfd1d/cz4xaViEiMFRSV8OrMtTzzeSbrc3fRN70Jf75gCGcN6UhSkiaHFxGpKtEkn3fHOwgRkXjZtrOIf361mue/zGRbQTEjurXg3rMHMqZvW92RSEQkAaIZ7f6pmXUFerv7ZDNrCCTHPzQRkQO3dmsBz3+5in/NWENBUSkn9GvLT0f3ZES3lokOTUSkTotmtPvVwFigJcGo947Ak4Buryki1cqqzTuZND+b9+dnMWfddpKTjLOGdOCa43rSt12TRIcnIiJEd9r9WmAkMB3A3ZeZWdu4RiUiEgV3Z9nGfCbNy2bS/CwWZ+8AYEinZtx6Sj/OPLQDHZs3SHCUIiISKZrkc7e7F+25NsrMUgjm+RQRSYi1Wwt4deZaJs7PYuWmnZjB8K4t+O33B3DKoHZKOEVEqrFoks9Pzew3QAMzOwn4OfBefMMSEfkud+erFVv4x5ermLwoBwOO6NGKK47sxskD29G2af1EhygiIlGIJvm8DbgSmAdcA0wEnolnUCIiexQUlfDW7PW88OUqlubk07JRGj8f3ZNLj+hK+2bq4RQRqWmiST4bAM+5+9MAZpYclhXEMzARqdvWbi3gxWmrGTdjDXmFJQzs0JQHzx3MGUM6UD9VE26IiNRU0SSfHwMnAvnhcgPgQ+DIeAUlInVTWZnz2bJNvDRtDR8vziHJjFMGteOKI7txWNcWmpdTRKQWiCb5rO/uexJP3D0/nOtTRCQmtu4sYvzXa3ll+hrWbC2gdeM0rh3di0uO6KJT6yIitUw0yedOMxvm7t8AmNlhwK74hiUitZ27882abbz41WomzsumqLSMkd1bcsvJfTl5YDvSUpISHaKIiMRBNMnnjcBrZrYhXG4PXBC/kESkNtszgOjFr1azOHsHTeqlcNHIzlxyRFf6pGsieBGR2i6a22vONLN+QF/AgMXuXhz3yESkVsneXsg/v1rFy9PXsH1XMQPaN+X+Hx7CmUM60KheNP8Hi4hIbRDtJ/4IoFu4/VAzw93/GbeoRKTWmL9+O89OzeS9ORsoc+d7A9px5THdGa4BRCIidVI093Z/keCe7t8CpWGxA0o+RWSvysqcjxdv5NmpK5m2ciuN0pL50aiuXHFkd7q00nhFEZG6LJqez+HAAHffr1tqmll94DOgXvg8r7v7XRHrbwYeBNq4++aw7HaCCe1Lgevd/YP9eU4RSaw1WwqYOD+LV2euJXPzTjo0q88dp/XngpGdaVo/NdHhiYhINRBN8jkfaAdk7Wfdu4Hjw6mZUoGpZjbJ3aeZWWfgJGDNno3NbABwITAQ6ABMNrM+7l66t8pFpHpYvjGf9+dnMXFeNguz8gAY2qU5j100lFMHtSMlWaPWRUTkv6JJPlsDC81sBkFCCYC7n1nZTmFP6Z75QVPDx57e0z8DvwbeidjlLGCcu+8GMs1sOTAS+CqKGEWkirg7S3J2MHFeNu/Pz2JpTvBnPqxLc+48vT8nD2xH55Y6tS4iInsXTfJ594FWHt6KcxbQC3jc3aeb2ZnAenefU26wQUdgWsTyurBMRBLM3Zm/Po9J87OYND+bzM07MYOR3Vpy9xkDOGVQe9o1q5/oMEVEpAawaC7lNLN0ghHvADPcfeN+PYlZc+At4AbgaeB77r7dzFYBw919s5k9Dnzl7i+F+zwLTHT3N8rVNRYYC5Cenn7YuHHj9ieUhMvPz6dx48aJDqNWUZvG1p72LHNnZW4ZX+eU8HVOKZt3OUkG/VsmMTw9hWHpKTSrp9Hq0dAxGltqz9hTm8ZWXW3PMWPGzHL34fvaLprR7ucTDAzKIJjn8zEzu8XdX482GHfPNbMMglPr3YE9vZ6dgG/MbCRBT2fniN06ARvKVYW7PwU8BTB8+HAfPXp0tGFUCxkZGdS0mKs7tWnslJU5T789hay8Nrw/P5vsvEJSk42je7Xm1EHtOWlAOi0apSU6zBpHx2hsqT1jT20aW2rPykVz2v0OYMSe3k4zawNMBipNPsPtisPEswFwIvBHd28bsc0q/tvz+S7wipk9QjDgqDcw4wBek4jsp9IyZ8LcDTw2ZTnLNxZSL2UNx/Vpw62H9OX4fuk0a6CR6iIiEhvRJJ9J5U6zbwGiGb7aHnghvO4zCRjv7hMq2tjdF5jZeGAhUAJcq5HuIvFVUlrGO99u4PFPlrNy8076pjdh7OB63HDOaN11SERE4iKab5f3zewD4F/h8gXApH3t5O5zgaH72KZbueX7gPuiiElEDkJxaRlvfbOev32ynDVbC+jXrglPXDKMkwe247PPPlXiKSIicRPNvd1vMbMfAkcTXPP5lLu/FffIRCTmikrKeH3WOv6esZx123YxqGNTnvrRYZzYP52kJA0eEhGR+Ksw+TSzXkC6u3/h7m8Cb4blx5pZT3dfUVVBisiBc3cWZe3g7W/X8/bs9WzcsZshnZvzu7MGMqZvW91fXUREqlRlPZ9/AX6zl/KCcN0ZcYlIRGJiQ+4u3vl2A2/PXs+SnB2kJBmj+7bh0iO6clyfNko6RUQkISpLPruF121+h7t/bWbd4haRiByw7buKmTQvi7dmr2d65lYADuvagnvPHsTph7SnpaZJEhGRBKss+azsdiUNYh2IiBy4+eu389zUTCbMy6KopIwerRtx00l9OOvQjnRppVtdiohI9VFZ8jnTzK5296cjC83sSoJbZopIApWWOR8tzOa5qauYsWorDdOSuXBEZ849rBOHdGym0+oiIlItVZZ83gi8ZWaX8N9kcziQBvwg3oGJyN5t31XM+Jlref7LVazP3UWnFg248/T+nDe8syaDFxGRaq/C5NPdc4AjzWwMMCgs/re7T6mSyETkOzI37+T5LzJ5bdY6CopKGdm9Jb/9/gBOGpBOsqZJEhGRGiKaeT4/AT6pglhEZC8WZeXx+CfLmTgvi5SkJM4Y0oErjurGoI7NEh2aiIjIftNtTESqqTlrc/nbJ8v5aGEOjdKSGXtsT35ydDfaNqlsLKCIiEj1puRTpJqZvnILf/tkOZ8v20yzBqnceGJvLj+yG80bapokERGp+ZR8ilQDpWXO1OWbeXzKcmas2krrxmncdmo/Lj2iK411n3UREalF9K0mkgDuzopNO/lqxWa+WL6Fr1ZuYfuuYto3q8/dZwzgghFdaJCWnOgwRUREYk7Jp0gV2ZC7iy9XbOHL5Zv5YsVmcvJ2A9CxeQO+NyCdY/q04eSB6dRLUdIpIiK1l5JPkThwd9ZsLWBG5lZmZG5l5qqtrNpSAECrRmmM6tmKI3u25qherejSsqEmhBcRkTpDyadIDJSVOcs25jMjcwvTw2RzT89mi4apDO/WkkuP6MpRvVrTN70JSZqXU0RE6iglnyL7oaS0jDVbC1i5aScrNuWzYlM+KzftZGnODvIKSwBo17Q+h3dvxcjuLRnZvSW92jRWsikiIhJS8ilSidIy57Wv1zJl8UZWbMpnzdYCikv9P+tbN65HjzaNOH1wBw7r2oLDu7ekU4sGOo0uIiJSASWfIhX4Zs027npnAfPWb6d760b0TW/CyQPb0aNNY3q2aUSPNo11L3UREZH9pORTpJwt+bv54/uLGf/1Oto1rc9jFw3l+4PbqzdTREQkBpR8ioRKSst4efoaHv5wCbuKS7nmuB5cf3xvGmmSdxERkZjRt6oI8PWqrfz2nQUsysrj6F6tufvMgfRq2zjRYYmIiNQ6Sj6lziouLWNm5lbGf72Wt7/dQIdm9XnikmGcMqidTrGLiIjEiZJPqVN2FBbz2dLNfLQwmymLN5JXWEK9lCSuHdOTa8f0omGa/iRERETiSd+0Uutlby/ko0U5fLQwh2krtlBUWkaLhql8b2A7ThqQzjG9WyvpFBERqSL6xpVao7i0jFWbd7IkZwdLs3cEP3Pyydy8E4BurRry4yO7ctKAdhzWtQXJmvhdRESkyin5lBrJ3VmwIY+pyzeTMaeQB779jBWb8v8zAXySQbfWjejXrgnnDe/ESf3T6dW2sa7lFBERSTAln1Jj5BYU8fmyzWQs2cRnyzaxaUdw7/RW9Y3BXetzXN829GvXhD7pTejZpjH1U5MTHLGIiIiUp+RTqq2yMmfe+u18unQTGUs28u3aXMocmjVI5ZjerRndty3H9m7Nwm+mMXr0yESHKyIiIlFQ8inVzvZdxbwxax0vTltN5uadmMHgjs34xfG9Oa5PGw79v/buPraq+77j+PvrB4yxHeMHcMF2MMQhDRASwDBCUIG0aZN0K8mSdklaES1Rn4aydtpDwzJpmqJWS1NV05SqW9SmImo31oSmZWtaRrPRLIEAMc+QQDAmYEMAY2NsMMYP3/1xfg4XYhts7gPGn5d05XO/95zD735k469/55x7ykdfcL7m7hSOVURERAZGzadcNd794BQvrn+fVzbX09bRxawJBSxdVMmim8ZQlJuV6uGJiIhIHKj5lJTq6Opmze6jLF93gA21jWRlpLH4tvEsub2CaaX5qR6eiIiIxJmaT0m6tnNdbDnYxLqaE7xUfYijp9opK8hm2T0f5wtV5RTkjEj1EEVERCRB1HxKwjWdPsfb7zex6UAjG2sb2VnfTGe3YwbzK4v5zv23sPCmsfrcTRERkWFAzafElbtT23CaLQdPsvlg1HDuPdoKwIj0NKaX5fPlT0xiTkUhMycUkJ+dmeIRi4iISDKp+ZQrcvLMObYeOsmWgyfZeih6NLd1AJCblcHMCQV87tbxzK4o5Nby0frsTRERkWEuYc2nmY0EXgeywr/zsrv/vZk9C/wRcA6oAf7U3U+GbZYBjwNdwJ+7++pEjU8Gp7vb2VDbyKpt9WzY38j+cOtKM7ipJI97pn2MGdeP5rbyAirH5upQuoiIiFwgkTOf7cCd7t5qZpnAG2b2G2ANsMzdO83sGWAZ8C0zmwI8BEwFxgO/M7PJ7t6VwDHKZXrnyCl+uaWeVdsOc6T5LDkj0plXWcyDVWXcVj6a6WWjyc3SRLqIiIj0L2Hdgrs70BqeZoaHu/t/x6z2FvBgWF4MrHD3dqDWzPYBc4D1iRqj9O/wyTZ+tfUwv9xSz56jLWSkGQsmj+Fv772ZT91cQvYIHUIXERGRgUnoVJWZpQPVQCXwA3ffcNEqjwH/EZZLiZrRHnWhJkl0ur2TX+84wsrqOjbUNgIwa0IBTy+eymenj6dQH4MkIiIiV8CiCemuvooAAAs6SURBVMoE/yNmo4FXgCfcfWeoPQVUAX/s7m5mPwDWu/tPw+s/Bl5195UX7esrwFcASkpKZq1YsSLh44+n1tZWcnNzUz2MC7g7e5u6eaO+k40fdNLeBSWjjDtKM5g7LoOxo9JSPcR+XY2ZDmXKM/6UaXwpz/hTpvE1XPNctGhRtbtXXWq9pJyk5+4nzWwtcDew08weBf4Q+KSf737rgPKYzcqAw73s63ngeYCqqipfuHBhAkcef2vXriVRYz7b0cXxlnaa2zrIz86kODer30PjHzSfZeXmOl6urqO2ITqP874Z5Xxhdhkzry/AbGhcLJTITIcj5Rl/yjS+lGf8KdP4Up79S+TV7mOAjtB4ZgOfAp4xs7uBbwEL3P1MzCargH8zs+8TXXB0I7AxUeMbilrbO9lU28jeoy0cb2nnWEt7+HqW4y3tnDrb+ZFtckakU5yXRXFuFmNysyjOG0FRThbb6k7y+t7jdDvMmVjI0kWV3HvLxxg1QhcNiYiISOIkstMYBywP532mAT939/8KFxJlAWvCzNpb7v41d99lZj8HdgOdwNLhfqX72Y4uNh9sYn3NCd7c18C2uma6uqOJ4pGZaYzNG8nYvCwml+Qxv7KYMXlZjM0bSf6oTJrbOmhobaeh5RzHW9tpaGmn5ngrG2rbaTrTwbj8kfzZwkoenFVGRXFOit+piIiIDBeJvNp9OzCjl3plP9t8G/h2osZ0tTvb0cWuw6d4a3/UbL79fhPnOrtJTzOml+Xz9QU3MO+GIqaV5ZOXlTHow+IdXd1kpNmQOawuIiIi1w4dY02RtnNd7D7SzM76U+yob2ZnfTPvHWv9cGbz5nHXsWTuBOZVFjG7opC8kfG7DWVm+tV9AZGIiIhcu9R8JkHP/c7X1Zzg1e3tfGfL79l3rJXQZ1KcO4JppfncNaWEqePzmTOxUB9pJCIiItckNZ8JcqjxDOtrTrB+/wnW1TRw9FQ7APlZRtXEUdwzbRzTSvO5pTSfkuuydAhcREREhgU1n1eou9s52nKWuqY2DjScZtOBRtbvP8GhxjYgmtW8/YZibp9UxLwbiqjdsZFFi2aneNQiIiIiqaHm8zLtP97Kjvpm6praqGs6Q11TG4caz3D45FnOdXV/uF5+diZzJxXy+B0TmVdZzI1jcy+Y1TygGU4REREZxtR8XqaXquv44doaAIpyRlBWOIqppfncPW0cZQXZlBeOoqwgm4qiHNLT1GCKiIiI9EbN52X60twJ3D+jlLKCbH0Qu4iIiMggqYu6TKWjs1M9BBEREZEhTx/4KCIiIiJJo+ZTRERERJJGzaeIiIiIJI2aTxERERFJGjWfIiIiIpI0aj5FREREJGnUfIqIiIhI0qj5FBEREZGkUfMpIiIiIkmj5lNEREREksbcPdVjGDQzOw68n+pxDFAx0JDqQVxjlGl8Kc/4U6bxpTzjT5nG13DNc4K7j7nUSkO6+RyKzOxtd69K9TiuJco0vpRn/CnT+FKe8adM40t59k+H3UVEREQkadR8ioiIiEjSqPlMvudTPYBrkDKNL+UZf8o0vpRn/CnT+FKe/dA5nyIiIiKSNJr5FBEREZGkUfMZB2b2gpkdM7OdMbVbzWy9me0ws/80s+tiXpseXtsVXh8Z6rPC831m9s9mZql4P6k2kDzN7ItmtjXm0W1mt4XXlGcwwEwzzWx5qL9jZstitlGmDDjPEWb2k1DfZmYLY7ZRnoCZlZvZ/4bvt11m9o1QLzSzNWb2XvhaELPNspDbHjP7TExdmTLwTM2sKKzfambPXbSvYZ/pIPK8y8yqQ27VZnZnzL6GfZ64ux5X+AA+AcwEdsbUNgELwvJjwNNhOQPYDtwanhcB6WF5I3A7YMBvgHtS/d6u9jwv2u4WYH/Mc+U5iEyBR4AVYXkUcACoUKaDznMp8JOwPBaoBtKU5wV5jgNmhuU8YC8wBfgu8GSoPwk8E5anANuALGAiUKP/R6840xxgPvA14LmL9jXsMx1EnjOA8WF5GlCvPM8/NPMZB+7+OtB4Ufkm4PWwvAZ4ICx/Gtju7tvCtifcvcvMxgHXuft6j747XwTuS/zorz4DzDPWw8C/AyjPCw0wUwdyzCwDyAbOAaeU6XkDzHMK8FrY7hhwEqhSnue5+xF33xyWW4B3gFJgMbA8rLac8/ksJvoDqd3da4F9wBxlet5AM3X30+7+BnA2dj/KNDKIPLe4++FQ3wWMNLMs5RlR85k4O4HPheXPA+VheTLgZrbazDab2d+EeilQF7N9XahJpK88Y/0JoflEeV6OvjJ9GTgNHAEOAt9z90aU6aX0lec2YLGZZZjZRGBWeE159sLMKohmjTYAJe5+BKJf/kQzxxDldChms57slGkvLjPTvijTiwwizweALe7ejvIE1Hwm0mPAUjOrJpqiPxfqGUSHNr4Yvt5vZp8kmn6/mD6K4Ly+8gTAzP4AOOPuPefgKc9L6yvTOUAXMJ7okOZfmtkklOml9JXnC0S/YN4G/glYB3SiPD/CzHKBlcA33f1Uf6v2UvN+6sPWADLtcxe91IZtpgPN08ymAs8AX+0p9bLasMszI9UDuFa5+7tEh9gxs8nAZ8NLdcDv3b0hvPYq0bljPwXKYnZRBhxGgH7z7PEQ52c9IcpZefajn0wfAX7r7h3AMTN7E6gC/g9l2qe+8nT3TuAvetYzs3XAe0ATyvNDZpZJ9Ev9Z+7+i1A+ambj3P1IOFx5LNTruPDoR092+rmPMcBM+6JMg4HmaWZlwCvAEnevCWXliWY+E8bMxoavacDfAf8SXloNTDezUeGcugXA7jBd32Jmc8OVb0uAX6Vg6FelfvLsqX0eWNFTU56X1k+mB4E7LZIDzAXeVab96yvP8LOeE5bvAjrdXT/zMcL7/zHwjrt/P+alVcCjYflRzuezCngonEM3EbgR2KhMzxtEpr1SppGB5mlmo4FfA8vc/c2elZVnkOornq6FB9GM2xGgg+ivmseBbxBdDbcX+EfCB/qH9b9EdALyTuC7MfWqUKsBnovdZjg9BpHnQuCtXvajPAeRKZALvBS+R3cDf61MryjPCmAP0QUKvwMmKM+P5Dmf6NDjdmBreNxL9GkgrxHNFL8GFMZs81TIbQ8xVwsr0yvK9ADRhXSt4ft6ijIdXJ5Ef4Cejll3KzBWeUYP3eFIRERERJJGh91FREREJGnUfIqIiIhI0qj5FBEREZGkUfMpIiIiIkmj5lNEREREkkbNp4hInJjZ/WbmZvbxVI9FRORqpeZTRCR+HgbeILrjloiI9ELNp4hIHIR7Pt9B9IHzD4XaQjNba2Yvm9m7ZvazcFcTzOyAmf2DmW02sx09s6VmlmNmL5jZJjPbYmaLQz3dzJ4N9e1m9tU+hiIiclVT8ykiEh/3Ab91971Ao5nNDPUZwDeBKcAkoga1R4O7zwR+CPxVqD0F/I+7zwYWAc+G23M+DjSH+mzgy+HWkiIiQ4qaTxGR+HgYWBGWV4TnEN1zvM7du4lusVcRs80vwtfqmPqngSfNbCuwFhgJXB/qS0J9A9Ft/W5MxBsREUmkjFQPQERkqDOzIuBOYJqZOZBOdB/oV4H2mFW7uPD/3fZe6gY84O57Lvo3DHjC3VfH/x2IiCSPZj5FRK7cg8CL7j7B3SvcvRyoBeYPYl+rgSdizg2dEVP/upllhvrkcDheRGRIUfMpInLlHgZeuai2EnhkEPt6GsgEtpvZzvAc4EfAbmBzqP8rOnolIkOQuXuqxyAiIiIiw4RmPkVEREQkadR8ioiIiEjSqPkUERERkaRR8ykiIiIiSaPmU0RERESSRs2niIiIiCSNmk8RERERSRo1nyIiIiKSNP8PS2lv/r2W1SwAAAAASUVORK5CYII=\n",
"text/plain": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.figure(figsize=(11, 5))\n",
"plt.plot(raw_data['Annee'], raw_data['Moyenne'])\n",
"plt.title('Concentration Moyenne de CO2 par an au fil du Temps')\n",
"plt.xlabel('Annee')\n",
"plt.ylabel('Concentration Moyenne de CO2')\n",
"plt.grid(True)\n",
"plt.show()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Oscilation périodique : on calcule cette valeur en faisant la difference entre la valeur moyenne annuelle et la valeur actuelle.\n",
"On stocke le resultat dans une nouvelle colonne que l'on appelle \"Oscilation\""
]
},
{
"cell_type": "code",
"execution_count": 99,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"