Skip to content

Instantly share code, notes, and snippets.

@ashao
Last active July 5, 2019 15:59
Show Gist options
  • Save ashao/dee3f445118eab60c9344dcbd1e94392 to your computer and use it in GitHub Desktop.
Save ashao/dee3f445118eab60c9344dcbd1e94392 to your computer and use it in GitHub Desktop.
Calculate the CO2 Annual Cycle at Barrow
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"from urllib.request import urlopen\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from scipy.signal import detrend"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"# Parse the Barrow CO2 record and only add entries after 1974 which doesn't have January and December\n",
"data_url = urlopen(\"https://cdiac.ess-dive.lbl.gov/ftp/trends/co2/barrsio.co2\")\n",
"data = []\n",
"for line in data_url:\n",
" first_entry = str(line).split('\\\\t')[0][-4:]\n",
" if first_entry.isdigit():\n",
" if(int(first_entry)>1974):\n",
" data.append([float(entry) for entry in str(line).replace('b','').replace('\\\\n','').replace('\\'','').split('\\\\t')])\n",
" \n",
"years = np.array(data)[:,0]\n",
"month_co2 = np.array(data)[:,1:-1]\n",
"annual_co2 = np.array(data)[:,-1]"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {},
"outputs": [],
"source": [
"month_anom = month_co2 - annual_co2[:,None]\n",
"month = range(1,13)"
]
},
{
"cell_type": "code",
"execution_count": 41,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[ 3.26090909 4.06545455 4.47424242 4.49 4.91666667\n",
" 2.89030303 -4.47575758 -10.78787879 -9.03848485 -3.74909091\n",
" 0.52545455 3.42333333]\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEGCAYAAABLgMOSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3deXxU9b3/8dcn+04IkAQIEJAlEAQkiCAuoAioVVRstbXWrUWr1l7rz6rd7b3db+1ttS7U3daCFakbgqJQULDsWwyyLyEJYc++zuf3xww2xiQkk5k5M5PP8/GYR2bOnJl5fx3MJ+ec7yKqijHGGNNREU4HMMYYE5qsgBhjjPGKFRBjjDFesQJijDHGK1ZAjDHGeCXK6QCB1LNnT83OznY6RrtUVlaSmJjodAy/sLaFrnBun7WtdevWrTuiqr2ab+9SBSQ7O5u1a9c6HaNdli1bxuTJk52O4RfWttAVzu2ztrVORPa1tN1OYRljjPGKFRBjjDFesQJijDHGK1ZAjDHGeMUKiDHGGK9YATHGGOMVKyDGGGO8YgXEGGOMV4K2gIjIXhHZIiIbReQLo//E7U8islNENovIWCdymuB23VOruO6pVU7HMCYsBftI9CmqeqSV5y4Fhnhu5wBPeH4aY4wJgKA9AmmHmcCL6vYxkCoivZ0OZYwxXYUE65K2IrIHOA4o8JSqzmn2/FvAr1X1Q8/j94EHVHVts/1mA7MBMjIy8ubOnRuI+J1WUVFBUlKS0zH8IhBtU1WO1Sj/t76G+Ej4wYQEv37eKeH8vUF4t8/a1ropU6asU9VxzbcH8ymsSapaJCLpwHsisk1Vlzd5Xlp4zReqoafwzAEYN26chspkaTaxW/vV1Dey/VA5BcVlFBS7f24rKedkdT0AKXFRnHXOJLrFR/vsM1sTzt8bhHf7rG0dF7QFRFWLPD9LRWQBMB5oWkAKgX5NHmcBRYFLaAJNVSk+WcO2kv8UioLiMvYcqcTl+dMhPjqSYZnJXHZmb4b3TubpFXs4cKyKa59YybM3n02/tMAciRjTFQRlARGRRCBCVcs996cBP2+22xvA3SIyF/fF85OqWhzgqCHvVA+lebdPDNjnnThRzen+GDp1VLGtuJxPiss+KxqnjioAsrrHM7x3Cpef2ZvhvVPI6Z3CgLQEIiL+c3D69uZiEmIiKTpRzdWPr+SZm8Yxul+qn1pnTNcSlAUEyAAWiAi4M76sqotE5A4AVX0SWAhcBuwEqoBbHMpqOkFVKSmr+dzpp+ZHFQkx/zmqGNE7mZzeKQzLTCYlrn2npLrFR/PY187i5ufWcN2cVfzx+rOYnpvpx1YZ0zUEZQFR1d3A6Ba2P9nkvgJ3BTJXOFJVGlzK4fJaFEUVXKq41P3cqcdf2I7nsYvPva7F/Vyexygnq+s5Was8/Gb+Z9cqTlR17KjCG4PTk1lw5yS++eJa7vjrOn50+QhunZSN548UY4wXgrKAGP84XF7rPi1UUs6nJWV8eqiCLYUncCmc/YslAc0yb80BhmUmc+lI744q2v05TU7N9UqOZe63JnDvvI3891ufsP9oJT+5IpfIThYnY7oqKyBBqL3XCVpTUdvA9kPlbC85VSzK2X6onKOVdZ/t0yMxhmGZyfRKjiUuKpJvnj8QESFChAgBEb7wOELEsw2EU9vls+ciPD9p9ljgs9f97M18aioreP/B6Z0+qvBGfEwkj98wll+9U8BfVuyh8Hg1f/rqWSTG2v8KxnSU/V8TwuoaXOw5Usm2kjK2H3IXim0l5RQer/5sn4SYSIZkJDN1eAZDM5PJyUxmWGYyPZNigf9cRL9xYnZAMqfEReOqEUeKxykREcIPLx9B/x6J/PT1rXzlqVU8e/PZZKTEOZbJmFBkBSQEuFzKwRPVfFpSzqeeU1DbS8rZfaSC+kb3leaoCGFQr0TG9Evl+rP7MTQjmZzMFLK6xzv6yzqY3ThhAFmp8dz18nqu/vNHPHvL2eRkpjgdy5iQYQUkyDQ0uiirrud4jfLg/M1sKylnx6FyKusaP9unb2o8OZnJXDQ8nZzMZIZmJDOoVyKxUZEd/rxAdd9t+nnLli0L6Ge2ZUpOOv+4YyK3Pr+Ga59YxeM3jOWCob2cjmVMSLACEgRcLmXN3mO8ubmId7aUfHatYnF+CcMyk7k2L4thme6LzEMzkkj28YXmri63Tzf+edckbnluDbc8v4b/uWokXx3f3+lYxgQ9KyAOUVU2HjjBm5uKeXtLEYfKaomLjmDq8Aw+LSnHVVvJkgcvsW6mAdK7Wzz/uGMid7+8gYde28L+Y1XcP22Ynf4zpg1WQAJIVfmkuIy3Nhfz5qYiCo9XExMZwYXDenHF6D5cnJNOYmyUuxdWfZUVjwBLjovmmZvG8ZM38nli2S72H6vi918eTVx0x08NGtMVWAEJgJ2lFby5qYg3Nxex+3AlkRHCeYN78l9Th3LJiIyATPJn2icqMoJfXDWSAWkJ/OqdbZScrOEv3xhHWmKM09GM8UpnhwW0xQpIO3gzX9T+o1W8ubmItzYXU1BchghMGNiDb543iBkjM+0XUhATEW6/8Az6pSVw77yNXP34Rzx389kM6hWeU30b4y0rID5UfLKatzcX8+bmYjYdOAHA2P6p/PSKEVx+Zm/S2znOINh6KnVVl53Zm4yUOGa/uJZrnljJnBvHMX5gmtOxjAkaVkA66UhFLe9sKebNTcWs3nsMgJF9U3jo0hwuH9WbrO42fXgoyxvQnQV3TuLm51fz9af/ze++PIqZY/o6HcuYoGAFxAsnq+pZlF/MW5uL+WjnEVwKQ9KTuO+SoXxpdB8G9kx0OqLxof49Enjt2+dy+0vr+O7cjew/WsXdFw22Tg4mJJRV11NY4eJkdb3Pr7daAWmnRpfyzw0HeXNTEct3HKa+URnQI4E7Jw/mitF9GJaZ7HRE40epCTG8eNt4Hpy/hd+/t539x6r45TVnEh0Z4XQ0Y1pVeLyKHaUVRAD+6JFuBeQ0VJWdpRUcq6pj7b7j9OkWxy2TBnLFqD6M7Jtif4V2IbFRkTzyldH0S0vgT+/voOhkNY/fkGe96ExQqq5rZPaL61Cgb5L4ZQCyFZDTODWLbHpyLH/+2ljG9u9ug8u6MBHhe5cMpX9aAg+9tplrn1jJc7ecbde6TFBRVe5/dRMFJWUMTU+Cuiq/fI4df7fDoF5JZPdIZFx2mhUPA8C1eVm8cOt4SspquOrPK9lceMLpSMZ85sl/7eatzcV8f3oOqQn+GzJgBcQYL517Rk8W3HkucdERfOWpVVzyyL/41b+rT/9CY/xo6bZSfrt4G1eM7sMdFw7y62dZATGmE04tlTssM4UdpRUcr1GnI5kubPfhCu6Zu4ERvVP47axRiAjzbp/IQ+fE++XzrIAY00mnlspNjY+mtFopOVnjdCTTBZXX1POtF9cSHRnBUzfmER/j/zncrIAY4wPxMZH0T3NfSH/3kxKH05iuxuVS7p23kX1Hq3j8hrEB69QRdAVERPqJyFIRKRCRfBH5bgv7TBaRkyKy0XP7iT8zzbt9YsAXXjKhJz4mkpgI9zouxgTSH5ZsZ0lBKT+5YgQTBvUI2OcGXQEBGoD7VHU4MAG4S0RGtLDfClUd47n9PLARjfmiebdPZFp2NB/vPsZxz6Jgxvjbwi3FPPrBTq4/ux83ThgQ0M8OugKiqsWqut5zvxwoAGzyIRMS8jIiaXQp728rdTqK6QIKisu475VNjO2fysMzcwM+sFlUg7fXiIhkA8uBkapa1mT7ZGA+UAgUAf9PVfNbeY/ZwGyAjIyMvLlz5/o3tI9UVFSQlBSe04eHc9vKyyv46boIBqRE8N2x7Zt9OZSE83cXam2rqFMeXlVNvQt+NjGO1LjWjwc627YpU6asU9VxX3hCVYPyBiQB64BrWnguBUjy3L8M2NGe98zLy9NQsXTpUqcj+E24t+2nr2/VIT9cqBU19U7H8blw/+5CRX1Do351ziod8sOFumH/8dPu39m2AWu1hd+pQXcKC0BEonEfYfxNVV9r/ryqlqlqhef+QiBaRHoGOKYxLZqem0ldg4t/bT/sdBQTpn6xsICVu47yq6vPZEy/VMdyBF0BEfdJvGeAAlV9pJV9Mj37ISLjcbfjaOBSGtO6s7O7k5YYw6Kt1hvL+N6r6wp57qO93DppILPyshzNEoyTKU4CbgS2iMhGz7YfAP0BVPVJ4Frg2yLSAFQD13sOs4xxXFRkBFOHp/POlhJqGxqJjfL/gC7TNWzYf5wfLNjCpME9+MFlOU7HCb4CoqofAm12JVDVx4DHApPImI6bMTKTV9YWsnLXUaYMS3c6jgkDpWU13PHXdWSkxPLYV8cSFQRr0TifwJgwdO4ZPUmMieRdG1RofKC2oZHb/7qO8poG/vKNcXRP9N8Mux1hBcQYP4iLjmRKTjrv5h+i0WVnV433VJWf/DOfDftP8PsvjyYnM8XpSJ+xAmKMn0zPzeRoZR3r9h13OooJYS99vI95aw/wnYsGc+mZvZ2O8zlWQIzxkyk56cRERlhvLOO1VbuO8vCbnzB1eDr3Th3qdJwvsAJijJ8kxUZx/pCeLM4vwToJmo4qPF7FXS+vZ2DPRP5w3ZigXA3VCogxfjQ9N5ODJ6rJLyo7/c7GeFTXNTL7xXXUN7qYc2MeyXHRTkdqkRUQY/xo6ogMIsSmeDftp6rc/+omCkrKePSrZzGoV/DOz2UFxBg/SkuMYfzANLsOYtrtiX/t4q3NxXx/eg6Tg3wMkRUQY/xsRm4mO0or2HW4wukoJsgt3VbK7xZ/yhWj+3DHhYOcjnNaVkCM8bNpuZmAncYybdt9uIJ75m5gRO8UfjtrVMDX9vCGFRBj/KxPajyjs7qx2E5jmVaU1dTzrRfXEh0ZwVM35hEfExrzp1kBMSYApuVmsqnwJEUnqp2OYoKMy6XcO3cj+45W8fgNY8nqnuB0pHazAmJMAMwY6T6NZXNjmeYeeW87728r5adXjGDCoB5Ox+kQKyDGBMAZvZIYnJ7E4vxDTkcxQWThlmIeW7qT68/ux9cnDHA6TocF3XTuxoSrGbmZPL5sJ8cq60gLktlUjTOue2oVVXUN7CytZGz/VB6emRsSF82bsyMQYwJkxshMXApLCuwopKurb3Sx/VAFKfFRPPn1vJBddMwKiDEBktsnhb6p8dYbq4trdCk7Syuoa3Tx1I3jSE+JczqS16yAGBMgIsL03ExW7DxCRW2D03GMQ15dd4CymgayeyQypl+q03E6xQqIMQE0PTeDugYXyz4tdTqKcUBVXQO/f3c7SbFR9EoK/etgVkCMCaBx2Wn0SIyx3lhd1DMr9lBaXkv/tPiQvGjeXNAWEBGZISKfishOEXmwhedjRWSe5/l/i0h24FMa0zGREcIlIzL4oOAQNfWNTscxAXS4vJYn/7WL6bkZQTs9e0cFZTdeEYkE/gxcAhQCa0TkDVX9pMlutwHHVXWwiFwP/Aa4LvBpjemY6SMzmbvmACt3HeGinAyn45gA+eP726lpcPHAjJygnqK9I4L1CGQ8sFNVd6tqHTAXmNlsn5nAC577rwIXSzgcE5qwd+4ZPUiKjWLxVjuN1VXsOlzB31cf4Gvj+4dN8YAgPQIB+gIHmjwuBM5pbR9VbRCRk0AP4EjTnURkNjAbICMjg2XLlvkpsm9VVFSETNaOsrbByDRl4aYDTEs7SmQQLlXaGvvuvPOn9TVEi3J2/GFH/vv5q23BWkBa+j+q+aLS7dkHVZ0DzAEYN26cTp48udPhAmHZsmWEStaOsrZBZVoxd728nsTsUSE1/5F9dx23Zu8x1i9axX2XDOXKi4f4/P3bw19tC9ZTWIVAvyaPs4Ci1vYRkSigG3AsIOmM6aTJw3oRExVhKxWGOVXllwsLyEiJ5ZvnB/8CUR0VrAVkDTBERAaKSAxwPfBGs33eAG7y3L8W+EBVv3AEYkwwSoyN4oIhvXg3vwT7Zxu+3tlawob9J/jeJUNDZo2PjgjKAqKqDcDdwGKgAHhFVfNF5OcicqVnt2eAHiKyE/ge8IWuvsYEs+m5GRSdrGHLwZNORzF+UNfg4jeLtjE0I4lr8/qd/gUh6LTXQERkPvAs8I6quvwfyU1VFwILm237SZP7NcCXA5XHGF+bOjyDyAhhcX4Jo7JCe0oL80Uv/3sf+45W8dzNZ4dUR4mOaM8RyBPA14AdIvJrEcnxcyZjuoTuiTGcMzDNroOEobKaev74/g7OPaMHk4f1cjqO35y2gKjqElW9ARgL7AXeE5GVInKLiITHcEpjHDJjZCa7Dleys7Tc6SjGh55ctovjVfU8dOnwsJiypDXtugYiIj2Am4FvAhuAP+IuKO/5LZkxXcC0Ee6lbm1urPBRdKKaZz7cw1Vj+nBmVjen4/jVaQuIiLwGrAASgCtU9UpVnaeq3wHCZ0ilMQ7I7BbHmH6pLLa10sPGI+9tRxXumzbM6Sh+154jkMdUdYSq/kpVi5s+oarj/JTLmC5jem4mmwtPcvBEtdNRTCd9UlTG/PWF3Dwpm35pCU7H8btWC4iIXCMi1wCpp+43vQUwozFhbXque0JFW6kw9P160TZS4qK5a/Jgp6MERFvdeK9o4zkFXvNxFmO6pEG9khiakcTi/BJuPW+g03GMl1bsOMzy7Yf50eXD6ZbQNfoXtVpAVPWWQAYxpiubkZvJY0t3crSilh5JsU7HMR3kcim/XLiNrO7x3DhxgNNxAqa9vbAuF5Hvi8hPTt38HcyYrmT6yExcCksKrDdWKFqw4SAFxWXcP30YsVHhN2VJa9rTC+tJ3As1fQf3DLhfBrpOiTUmAEb0TiGre7wNKgxBNfWN/P7dTzmzbzeuGNXH6TgB1Z4jkHNV9Ru4V/97GJjI52fKNcZ0kogwIzeTj3Yepbym3uk4pgOe+2gvRSdreOiyHCLCdMqS1rSngJzqW1glIn2AesCu9BnjY9NHZlLX6GLpp4edjmLa6VhlHY8v3clFOemce0ZPp+MEXHsKyFsikgr8DliPezqTuf4MZUxXNLZ/d3omxdqgwhDy6Ac7qKxr4MFLu+YUgaedjVdV/9tzd76IvAXEqarNP22Mj0VGCJeMyOCNjQepqW8kLrrrXIwNRfuOVvLXj/fxlXH9GJqR7HQcR7RnOvdI4HIg+9T+IoKqPuLfaMZ0PTNGZvL31fv5aOcRLh6e4XQc04bfLv6UqIgI7r1kqNNRHNOeU1hv4p5IsQeQ3ORmjPGxiYN6kBwXZb2xgtyG/cd5e3Mx3zp/IBkpcU7Hccxpj0CALFUd5fckxhhioiK4OCedJQWHaGh0ERUZlIuGdmmqyq8WbqNnUgyzLzzD6TiOas+/zndEZJrfkxhjAPfkiser6lm995jTUUwLlhSUsnrvMb47dShJse35Gzx8taeAfAwsEJFqESkTkXIRKfN3MGO6qguH9SI2KoJ3bY2QoNPQ6OLX7xQwqFci159tw+HaU0B+j3vwYIKqpqhqsqqm+DmXMV1WQkwUFw7txaKtJbhc6nQc08S8tQfYdbiSB2bkEG2nF9tVQHYAW1XV/iUbEyDTczMpKath80HrMR8sKmob+MN7Ozg7uzvTRlgPOWjfRfRiYJmIvAPUntroj268IvI73NPI1wG7gFtU9UQL++0FyoFGoMEWtjLh5uLh6URFCIvzSxjTL9XpOAb4y/LdHKmoZc438sJ6nfOOaM8RyB7gfSAG/3fjfQ8Y6en1tR14qI19p6jqGCseJhylJsQwYVAPFm8twQ7+nVdaVsOc5bu5/MzejO3f3ek4QaM9I9EfBhCRZPdDrfBXGFV9t8nDj4Fr/fVZxgS76SMz+fE/t7KztIIhXXSkc7D4w5IdNLhc3D89/Nc57wg53V83IjISeAlI82w6AnxDVfP9GkzkTWCeqv61hef2AMdxr4z4lKrOaeN9ZgOzATIyMvLmzg2NabwqKipISkpyOoZfWNva53iNi3uXVXPNkGiuPCPGJ+/ZWV3xuztY4eJHH1YzdUAUNwwPzcW+Ovu9TZkyZV2LZ3tUtc0bsBL36aJTjycDK0/3ujbebwmwtYXbzCb7/BBYgKfAtfAefTw/04FNwAXt+ey8vDwNFUuXLnU6gt9Y29rv6j9/qJf/ablP37MzuuJ3d9vzq3XkTxbp0YrawAbyoc5+b8BabeF3ansuoieq6tImBWeZiCR2vIZ99vqpbT0vIjcBXwIu9gRv6T2KPD9LRWQBMB5Y7m0mY4LV9NxMfvXONg4cq6JfWoLTcbqcj3cfZUlBKd+fMYy0xOA4Cgwm7bmIvltEfiwi2Z7bj3BfWPc5EZkBPABcqapVreyT6Lkeg6eQTcN9BGNM2JmemwnAu5/YoMJAc69zXkDvbnHcOsmWQGpJewrIrUAv4DXcp5V6Abf4Kc9juHt4vSciGz3L6SIifURkoWefDOBDEdkErAbeVtVFfspjjKOyeyaSk5nMYptcMeDe2lLM5sKT3DdtmE2t34r29MI6DtwTgCyo6uBWthcBl3nu7wZGByKPMcFgem4mf/pgB4fLa+mVHJoXcUNNbUMjv1u8jZzMZK4+q6/TcYLWaY9ARGSoiMwRkXdF5INTt0CEM8a41whRhSUFdhorUF5atY8Dx6r5wWXDiexi65x3RHsuov8DeBJ4GvfIb2NMAOVkJtM/LYHF+SV8dXx/p+OEvZNV9Tz6wU7OH9KTC4b2cjpOUGtPAWlQ1Sf8nsQY0yIRYcbITJ77aA9lNfWkxEU7HSmsPb5sJ2U19V12nfOOaNeKhCJyp4j0FpG0Uze/JzPGfGZ6bgb1jcrSbaVORwlrhcereG7lXq4+qy+5fbo5HSfotecI5CbPz/ubbFNgkO/jGGNacla/7vRKjmVxfgkzx9hFXX/5/bvbAfh/02zKkvZoTy8s6wBtjMMiIoRpIzJ4bf1BauobrVupj1331CpKjlSxr7ySOy48gz6p8U5HCgntWhFFREaKyFdE5Bunbv4OZoz5vBkjM6mub2TFjiNORwk7qsrhaqV7QjR3Tuna65x3RHu68f4UeNRzmwL8FrjSz7mMMc1MGNSDlLgoFtmgQp87WV1PVQN856Ih1kmhA9pzBHItcDFQoqq34B7EZ6OZjAmw6MgIpg7P4P1th6hvdDkdJ2yoKgeOVxMdAV+fMMDpOCGlPQWkWlVdQIOIpACl2AV0YxwxLTeTE1X1rN5zzOkoYWPVrqNU1TWSFifERNk65x3Rnv9aa0UkFfgLsA5Yj3sOKmNMgF04tBdx0REszrfTWL7yzId7iIoQUmyy3Q47bQFR1TtV9YSqPglcAtzkOZVljAmw+JhIJg9NZ3F+CS6XLXXbWbsOV/D+tlIyUmKJsHXOO6xDx2uquldVN/srjDHm9KaPzOBQWS2bCk84HSXkPffRHmIiI/jnXefx0DnWdbej7ISfMSHmopwMoiKERXYaq1OOV9bx6rpCrjqrj81y7CUrIMaEmG7x0Uw8owfPf7SXrzy50uk4Ievl1fupqXdx63k2VtpbVkCMCUEzRmZS2+Ciut4myPZGXYOLF1bu5fwhPcnJTHE6TshqtYCIyJki8rGIHPCsB9K9yXPWC8sYB10yIgOAY5X1DicJTW9vKaK0vNaOPjqprSOQJ4CfAWcC23EvI3tqjL8N1TTGQenJcSTHRnG0shZV643VEarK0yv2MDg9iQuH2HofndFWAUlS1UWeLrz/C9wNLBKRCbhn4zXGOKhnUgw19S42FZ50OkpI+Xj3MfKLyrh10kAibLXBTmmrgIiIfDYhvqouBWYBLwE23t8Yh6UlxSAC89cVOh0lpDzz4R66J0RzzVibFr+z2iogvwGGN93gGQNyMfCavwKJyM9E5KCIbPTcLmtlvxki8qmI7BSRB/2Vx5hgFRURQVpCDG9sKqK2wS6mt8eeI5W8v+0QN04YYFPi+0CrBURVX1bVjwFEJElEEj3b96vqt/yc6w+qOsZzW9j8SRGJBP4MXAqMAL4qIiP8nMmYoDLv9ok8ct0YTlbX80GBrVTYHs99tIfoiAi+PtFOovhCm914ReTbIrIf2AccEJF9InJnYKK1aTywU1V3q2odMBeY6XAmYwLuvME9yUiJZf56O411Oieq6vjH2kKuHNOH9OQ4p+OEhVZXJBSRHwHnApNVdbdn2yDgjyKSpqr/48dcd3sWrVoL3Keqx5s93xc40ORxIXBOS28kIrOB2QAZGRksW7bM92n9oKKiImSydpS1zbfyerhYvK2UNxYvJSXWvxeFQ/m7e3t3HdX1jYyKPdJiG0K5bafjt7apaos34FMgroXt8cD21l7XnhuwBNjawm0mkAFE4j46+gXwbAuv/zLwdJPHNwKPnu5z8/LyNFQsXbrU6Qh+Y23zre0lZTrggbf06RW7/f5Zofrd1TU06jm/WKJf+8uqVvcJ1ba1R2fbBqzVFn6ntrkmuqrWtLCtWkQ6tZqNqk5tz34i8hfgrRaeKgT6NXmcBRR1JpMxoWpIRjKjsrrx6rpCbrOBcS1auKWYkrIafnnNSKejhJW2roEUisjFzTeKyEVAsb8CiUjvJg+vxn1k0twaYIiIDBSRGOB64A1/ZTIm2M0am0VBcRmfFJU5HSXoqGfg4KBeiUwemu50nLDSVgG5B3hKRJ4Xke+IyN0i8gIwB/egQn/5rYhsEZHNuNdgvxdARPqIyEIAVW3wZFgMFACvqGq+HzMZE9SuHN2H6Eixi+ktWLP3OFsOnrSBg37Q6iksVc0XkZHA14BcQIDlwO0tndryFVW9sZXtRcBlTR4vBL7QxdeYrqh7YgwX52Tw+saDPHhpDtGRNk/qKU+v2E1qQjSzxmY5HSXstNULazCQoarPNtt+vogUqeouv6czxrTbrLwsFuWXsHz7YS4enuF0nKCw90gl7xUc4q7Jg4mPsYGDvtbWnyn/B5S3sL3a85wxJohMHtaLHokxdhqriedX7iUqQviGDRz0i7YKSLa2sHytqq4Fsv2WyBjjlejICK4c04cln5RyoqrO6TiOO1ldzytrD3DF6D6kp9jAQX9oq4C09V/cFg82JgjNGptFXaOLNzdZr/a5q/dTVddoXZv9qK0CskZEvjDnlYjcBqzzXyRjjLdy+6SQk5nMq+sPOh3FUfWNLp5fuZeJg3qQ26fb6V9gvMY9crMAABCDSURBVNLWQML/AhaIyA38p2CMA2Jwj88wxgQZEeHavCz+5+0CdpZWMDg9yelIjnhnawnFJ2v475k2cNCf2pqN95Cqngs8DOz13B5W1YmqWhKYeMaYjpo5pi+REV13TIiq8syK3QzsmchFOTZw0J9O21lcVZeq6qOe2weBCGWM8V6v5FguHNqLBesP0ujqeouHrtt3nE2FJ7l1UrYNHPQzG21kTBiaNTaLkrIaVu464nSUgHt6xR66xUczK88GDvqbFRBjwtDFw9NJiYvqcsvd7j9axbuflHDDOf1JiGlzrljjA1ZAjAlDcdGRXDG6D4vySyivqXc6TsA8t3IPESJ8Y2K201G6BCsgxoSpa/OyqKl3sXCL3ybPDiplNfW8ssY9cDCzmw0cDAQrIMaEqTH9UhnUK5H567rGmJB5qw9QaQMHA8oKiDFhSkSYNTaL1XuPsf9oldNx/KrBM3DwnIFpjOxrAwcDxQqIMWHsmrF9ESHsx4Qsyi/h4IlqO/oIMCsgxoSx3t3imXRGT17bUIgrjMeEPL1iD9k9Emwa+wCzAmJMmJuV15cDx6pZvfeY01H8Yt2+42w8cIJbJg0k0gYOBpQVEGPC3PTcTJJiw3dMyDMf7iYlLoprbeBgwFkBMSbMJcREcdmZmSzcUkxVXYPTcXzqwLEqFm0t4WvnDCAx1gYOBpoVEGO6gFljs6isa2RxfnjNg/r8yr1EiHDTubbioBOsgBjTBZydnUa/tPiwGhNSXlPPvDUHuHxUb3p3szXunBBUBURE5onIRs9tr4hsbGW/vSKyxbPf2kDnNCbUREQI15yVxUe7jlB0otrpOD4xb80BKmobrOuug4KqgKjqdao6RlXHAPOB19rYfYpn33EBimdMSJs1NgtVWLAh9I9CGhpdPPfRXsZnpzEqK9XpOF1WUBWQU0REgK8Af3c6izHhon+PBMZnpzF/XSGqoT0m5N1PDnHwRDW32tGHoyQY/yGJyAXAI60dXYjIHuA4oMBTqjqnjfeaDcwGyMjIyJs7d64fEvteRUUFSUnhuRyptc05ywvreXZrHT+aEMfg1MgOvz5Y2vc/H1dzslb5zQXxRIhvxn4ES9v8obNtmzJlyroWfx+rakBvwBJgawu3mU32eQK4r4336OP5mQ5sAi5oz2fn5eVpqFi6dKnTEfzG2uacsuo6HfajhfqD1zZ79fpgaN+6fcd0wANv6bMf7vbp+wZD2/yls20D1moLv1MD3nFaVae29byIRAHXAHltvEeR52epiCwAxgPLfZnTmHCUHBfNjNxM3txUxI+/NIK46I4fhTjtmQ/3kBwXxZfH9XM6SpcXjNdApgLbVLXFYbMikigiyafuA9NwH8EYY9phVl4WZTUNvF9Q6nSUDis87hk4OL4/STZw0HHBWECup9nFcxHpIyILPQ8zgA9FZBOwGnhbVRcFOKMxIevcM3qSmRIXkjP0vrByLwA3nZvtaA7jFnQlXFVvbmFbEXCZ5/5uYHSAYxkTNiIjhKvH9mXO8t2UlteQnhwaq/dV1DYwd/UBLjuzN31SbeBgMAjGIxBjjJ/NGptFo0t5fUOR01Ha7ZU1Byi3gYNBxQqIMV3Q4PQkxvRLZf760BgT0uhSnlu5h3EDujOmnw0cDBZWQIzpomblZbGtpJz8ojKno5zWe5+UcOCYrTgYbKyAGNNFXTGqNzGRESFxMf3pFXvolxbPtNxMp6OYJqyAGNNFpSbEMHVEOm9sLKK+0eV0nFZtPHCCtfuOc/O5tuJgsLECYkwXNmtsFkcr61j26WGno7TqmQ/3kBwbxVfG2YqDwcYKiDFd2AVDe9EzKZZX1x1wOsrnXPfUKq57ahVFJ6pZuKWY68f3Izku2ulYphkrIMZ0YdGREVw1pg8fbCvleGWd03G+4IWVe1FVGzgYpKyAGNPFzcrLor5ReWNTcI0JaXQpL6/ez6Vn9iare4LTcUwLrIAY08UN753CiN4pQdcb63B5LeU1NnAwmFkBMcYwKy+LzYUn2XGo3OkogHuZiZKyGsb2T2Vs/+5OxzGtsAJijGHmmD5ERQivBslRSGl5LbUNLm47b5DTUUwbrIAYY+iZFMvkYb3454aDNLqcm9qkpr6RH/9zK3uPVpEcG8X03AzHspjTswJijAHg2rwsDpXVsmKHM2NC9hyp5JrHV/LSx/vITIkjp3cyUZH2KyqY2bdjjAFgSk46qQnRzF9/MOCf/frGg3zpTysoOlnNMzeNY0CPBJ+tdW78xwqIMQaA2KhIrhzdh3fzSyirqQ/IZ9bUN/LQa5v57tyN5PROYeE953PxcDttFSqsgBhjPjNrbBa1DS7e3lzs98/aWVrBVX/+iL+vPsC3J5/B3NkTbKGoEGMFxBjzmVFZ3RicnsT8df7tjTV/XSFXPPohpeW1PH/L2TwwI4dou94RcoJuSVtjjHNEhFljs/jNom3sPVJJds9En75/VV0DP309n3+sK2T8wDT+dP1ZZHb74pK6826f6NPPNf5hJd8Y8zlXn9WXCIHXfDwmZPuhcmY+9hGvri/kOxcN5uVvntNi8TChw5ECIiJfFpF8EXGJyLhmzz0kIjtF5FMRmd7K6weKyL9FZIeIzBORmMAkNyb8ZXaL47whvZi//iAuH4wJUVVeWXOAKx/7kONVdbx06zncN22YddENA059g1uBa4DlTTeKyAjgeiAXmAE8LiKRLbz+N8AfVHUIcBy4zb9xjelaZo3ty8ET1Xy852in3qeytoHvvbKJ78/fzNj+3Vl4z/mcN6Snj1IapzlSQFS1QFU/beGpmcBcVa1V1T3ATmB80x1ERICLgFc9m14ArvJnXmO6mum5mSTHRjF/nfdjQgqKy7ji0Q95feNB7p06lJduO4f0FDtlFU6C7SJ6X+DjJo8LPdua6gGcUNWGNvb5jIjMBmYDZGRksGzZMp+F9aeKioqQydpR1rbQcFYveGtTIZekHSMuyj2orz3tU1X+VdjAXwvqSIoW7h8Xx/Cog6xYHvgBih0RTt9dc/5qm98KiIgsATJbeOqHqvp6ay9rYVvzk7Dt2ec/T6jOAeYAjBs3TidPntzarkFl2bJlhErWjrK2hYbE7GMsf3IVld2HMCMvi+ueWsWJE9UsfmByq68pr6nnBwu28mZ+EecP6ckfrhtDz6TYwIXuhHD67przV9v8VkBUdaoXLysE+jV5nAU0X+XmCJAqIlGeo5CW9jHGdNK4Ad0Z0COB+esLmZV3+vXItx48yd0vr2f/sSrunz6Mb194BhERNh1JOAu2bhBvANeLSKyIDASGAKub7qCqCiwFrvVsuglo7YjGGOOlU2NCVu46SuHxqlb3U1VeXLWXax5fSU29i7mzJ3LXlMFWPLoAp7rxXi0ihcBE4G0RWQygqvnAK8AnwCLgLlVt9LxmoYj08bzFA8D3RGQn7msizwS6DcZ0BVef5b68uKCVCRZPVtdz18vr+cnr+Uwa3IOF3z2f8QPTAhnROMiRi+iqugBY0MpzvwB+0cL2y5rc302z3lnGGN/rl5bAhEFpvLbhIL2SPj/catOBE9z99/UUn6jhoUtz+Nb5g+yoo4sJtlNYxpggM2tsFnuOVFJR6+74qKo88+Eern1yJS6Xe9qR2+16R5dkBcQY06ZLz+xNfHQkRyrqaHQps19ax3+/9QkXDk3n7XvOI2+ArVneVQXbOBBjTJBJio3i0pGZvL6xiGOi7Pu0lB9/aQS3TspGbNGnLs2OQIwxp3VtXhaNqgjw6h3nctt5A614GDsCMcac3sQzejAsI4nGmkpG90t1Oo4JEnYEYow5LREhNSGGSLtQbpqwAmKMMcYrVkCMMcZ4xQqIMcYYr1gBMcYY4xUrIMYYY7xiBcQYY4xXrIAYY9pl3u0TeeiceKdjmCBiBcQYY4xXrIAYY4zxihUQY4wxXrECYowxxitWQIwxxnjFCogxxhivWAExxhjjFSsgxhhjvGIFxBhjjFdEVZ3OEDAichjY53SOduoJHHE6hJ9Y20JXOLfP2ta6Aaraq/nGLlVAQomIrFXVcU7n8AdrW+gK5/ZZ2zrOTmEZY4zxihUQY4wxXrECErzmOB3Aj6xtoSuc22dt6yC7BmKMMcYrdgRijDHGK1ZAjDHGeMUKSBARkX4islRECkQkX0S+63QmXxORSBHZICJvOZ3F10QkVUReFZFtnu9wotOZfEVE7vX8m9wqIn8XkTinM3WGiDwrIqUisrXJtjQReU9Ednh+dncyo7daadvvPP8uN4vIAhFJ9cVnWQEJLg3Afao6HJgA3CUiIxzO5GvfBQqcDuEnfwQWqWoOMJowaaeI9AXuAcap6kggErje2VSd9jwwo9m2B4H3VXUI8L7ncSh6ni+27T1gpKqOArYDD/nig6yABBFVLVbV9Z775bh/AfV1NpXviEgWcDnwtNNZfE1EUoALgGcAVLVOVU84m8qnooB4EYkCEoAih/N0iqouB4412zwTeMFz/wXgqoCG8pGW2qaq76pqg+fhx0CWLz7LCkiQEpFs4Czg384m8an/A74PuJwO4geDgMPAc55TdE+LSKLToXxBVQ8C/wvsB4qBk6r6rrOp/CJDVYvB/ccckO5wHn+5FXjHF29kBSQIiUgSMB/4L1UtczqPL4jIl4BSVV3ndBY/iQLGAk+o6llAJaF7CuRzPNcCZgIDgT5Aooh83dlUxhsi8kPcp8r/5ov3swISZEQkGnfx+JuqvuZ0Hh+aBFwpInuBucBFIvJXZyP5VCFQqKqnjhhfxV1QwsFUYI+qHlbVeuA14FyHM/nDIRHpDeD5WepwHp8SkZuALwE3qI8GAFoBCSIiIrjPoReo6iNO5/ElVX1IVbNUNRv3BdgPVDVs/opV1RLggIgM82y6GPjEwUi+tB+YICIJnn+jFxMmHQSaeQO4yXP/JuB1B7P4lIjMAB4ArlTVKl+9rxWQ4DIJuBH3X+cbPbfLnA5l2u07wN9EZDMwBvilw3l8wnNU9SqwHtiC+/dGSE/7ISJ/B1YBw0SkUERuA34NXCIiO4BLPI9DTittewxIBt7z/F550iefZVOZGGOM8YYdgRhjjPGKFRBjjDFesQJijDHGK1ZAjDHGeMUKiDHGGK9YATHGR0REReSlJo+jROSwtzMPe2b3vbPJ48nhOIuxCV1WQIzxnUpgpIjEex5fAhzsxPulAneedi9jHGIFxBjfegf3jMMAXwX+fuoJz3oT//SsyfCxiIzybP+ZZw2HZSKyW0Tu8bzk18AZnoFfv/NsS2qy5sjfPCPDjXGEFRBjfGsucL1nwaVRfH425YeBDZ41GX4AvNjkuRxgOjAe+KlnTrQHgV2qOkZV7/fsdxbwX8AI3DMAT/JnY4xpixUQY3xIVTcD2biPPhY2e/o84CXPfh8APUSkm+e5t1W1VlWP4J7EL6OVj1itqoWq6gI2ej7LGEdEOR3AmDD0Bu71MyYDPZpsb+l006m5hGqbbGuk9f8327ufMX5nRyDG+N6zwM9VdUuz7cuBG8Ddowo4cpr1XspxT4BnTFCyv16M8TFVLcS9PnpzP8O9YuFmoIr/TB3e2vscFZGPRGQr7ovzb/s6qzGdYbPxGmOM8YqdwjLGGOMVKyDGGGO8YgXEGGOMV6yAGGOM8YoVEGOMMV6xAmKMMcYrVkCMMcZ45f8DY6ldGU+y4P8AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.errorbar(month,month_anom.mean(axis=0),yerr=month_anom.std(axis=0))\n",
"plt.xlabel('Month'); plt.ylabel('CO2 anomaly');plt.grid();\n",
"print(detrend_month_anom.mean(axis=0))"
]
},
{
"cell_type": "code",
"execution_count": 42,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEGCAYAAABLgMOSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nO3deXxU9bnH8c+TfQ9JIIGwhT0BhCgRVKyGTa1VqHVtUakb9ra2vdrb1t62V9vb9vZ207a3tUVFcWlxqRbqUspiXBCVfZMlYQ8QQoAEQsj+3D/moDEmYUhmcmZ53q/XvDLnzJk532NwnvzOOb/fT1QVY4wx5mxFuB3AGGNMcLICYowxplOsgBhjjOkUKyDGGGM6xQqIMcaYTolyO0B36tmzp+bk5LgdwysnT54kMTHR7Rh+YccWvEL5+OzY2rd69eoKVe3Ven1YFZCcnBxWrVrldgyvFBUVUVhY6HYMv7BjC16hfHx2bO0TkT1trbdTWMYYYzrFCogxxphOsQJijDGmU6yAGGOM6RQrIMYYYzrFCogxxphOsQJijDGmU6yAmJBVVlXLE8t3UVFd53YUY0JSWHUkNOFhW9kJ5ry1k4Xr99PQpDy9Yg/P3DmB7B7xbkczJqRYC8SEBFVlxY4j3PbEB1z+8Fu8tvEgMycM5I8zz+PwiTqu/9MKdlecdDumMSHFWiAmqDU1K//cVMact3awvrSKjMQY7ps2nFsuGEhaYgwAA9ITuHXuB1z/5xU8c8cERvROdjm1MaHBCogJSqfqm3hh9T4ee3sXe4/WkJORwE+vGc215/UjLjryE9uO7pvKc7Mv4ObH3+fGOSuYd9t4xvbv4VJyY0KHFRATVI5U1/HUij08tWI3x2oaOHdAD/7zylymjexNZIS0+75hWcm8cPdFzHz8PWY+9j6PzypgwuCM7gtuTAiyAmKCwp4jJ3ns7V08v2ofdY3NTM3L5O5Lh1AwMA2R9gtHSwMyEjxF5LH3uHXuB/z5lnEUjsj0c3JjQpcVEBPQ1u2rZM5bO/jnpjKiIiK45ty+3HXJIIZmdu46Ru/UOJ6/+0JunfsBdz21it/edC5XntPHx6mNCQ9WQEzAaW5WiraX8+c3d/L+rqMkx0Vx96VDuO2iHDJT4rr8+RlJsfzlrgu4/cmV3POXNfziurFcN66fD5IbE16sgJiAUdfYxIJ1B3j0rZ0Ul1eTnRrHDz6Xx03jB5AU69t/qqnx0Tx9x3hmP7Wa/3hhPSfrGpl1UY5P92FMqLMCYlx3vLaBv7y/lyeW7+LQ8Tpyeyfz0I1juWpMNtGR/uuqlBATxWOzCvj6X9fywMLNVNc18rVJQ/22P2NCjRUQ45qDVad4Yvlu/vL+XqrrGpk4NINfXDeWS4b19PrCeFfFRUfyx5nn8e0X1vPLRduormvkO5eP6Lb9GxPMrICYbrfvRDP3Pb+OhesOoMCV5/Th7ksGM7pvqit5oiMj+M0N+STERvFI0Q6qaxv50fRRRHRwW7AxwaKypp4Pyhop9MNnWwExXVbb0ERlTQOVp+o9P2saOH6qxfKpBqqc149U17O17BTx0WXcfMFA7rh4EP3TE9w+BCIihJ9+fjRJsVHMeWsnJ+sb+cW1Y4jy4yk0Y/ytsamZe/6ylvd21DHzilP09fF4cFZADOC586m6vtHzRV/T+su/nqpTDZ8qBpU1DVSdaqCusbndz42KEFLjo0lNiKZHfDR9UuMYmVzLD28q/GiokUAhInzvs7kkx0bx68Xbqalr4rdfzCc2KvLMbzYmAP3vP7fyTkkFt4+O8XnxACsgIae2oemjL3vPT8+X/+nHR+ud5eMttmnW9j83PjqSHgnRpMZH0yMhmkE9E+kRH+NZlxD98fP4j7fpkRBDYkzkp64nFBUVBVzxOE1E+PqUYSTERvHfr3zIXU+t5s83jyM+xoqICS5/X7ufR9/exawLB3JJaoVf9mEFJAA1NjVzvF7Zebi61Rf9J4tA1an6TxWF+g5aAxECKfGeloCnVRDDgPQEUuOjSEuIafHl7ykGPZyWQ2p8dNj9FX7HxYNIio3k/pc2MmvuBzz25QJS4qLdjmWMVzbtr+K7f9vA+EHp/OCqkSx/+y2/7McKiIuqTjWw43A1JeXV7Civ/uj53qM1ntbAsjfbfF9iTCQ9EmI+KgZDeiV99Fd/ivMzNd7TKmi5Pjk2yi4Mn4Ubzx9AQkwU9z63jpmPvs9Tt48P2JaTMadVVNcx+6lVZCTG8MeZ5/n1VngrIH6mqhysqv24UHz08ySHT3w8U15MZAQ5PRMYmZ3C58b0obJsHwVj8ugRH/OJopAaH+3XfxDmk64em01ibCRfeWYNN87xDAfvi97wxvhDQ1MzX312DUdO1vO3f7uInkmxft2fqwVERK4AfgtEAo+p6s9bvR4LPAWMA44AN6rqbhGZBvwciAHqgW+r6rJuDd9KfWMze46c/ESBOF0wauqbPtouJS6KoZlJFA7vxZDMJIb2SmJIZhL90+I/ccdPUVEZhefa8BqBYHJuFk/edj53zVv10ZwigXDnmDGt/eSVD/lg11EevjG/W26Ld62AiEgk8AdgGlAKrBSRhar6YYvN7gCOqepQEbkJ+F/gRqACuFpVD4jIaGAR0Lc7ch+vbXBON31cIHaUV7PnaA1NLa5CZ6fGMSQziRsK+jM0M4khvZIYmplEz6QY66QWhC4a0pNn7pzArLkfcMOfV/DMnRMY0ivJ7VjGfOT5lfuYt2IPd31mEJ8/t1u+Dl1tgYwHSlR1J4CIzAdmAC0LyAzgQef5i8D/iYio6toW22wG4kQkVlXr8IOfv76VDaWVlJRXU97itFN0pJCTkciI3slceU6fjwrF4F6JJPp47CbjvnMHpPHc3Rdyy+Pvc8OfVvD0HRMYmZ3idixjWLP3GD/4+yY+M6wn370it9v2K6od3Lvpzx2LXAdcoap3Osu3ABNU9Z4W22xytil1lnc421S0+pyvqOrUdvYzG5gNkJWVNW7+/PlnnfWXK09xqhGykyLokyjOzwh6xUuHkxh1RXV1NUlJofkXbrAfW9nJZn6xspbaRuW+cXEMTfv4DrVgP7YzCeXjC9Zjq6xt5sEVtURHwAMXxpMU8+nvpK4e26RJk1arasGnXlBVVx7A9Xiue5xevgX4fattNgP9WizvADJaLI9y1g3xZp/jxo3TYPHGG2+4HcFvQuHY9h09qZf+Ypnm/fB1faf48EfrQ+HYOhLKxxeMx1bb0KjX/OEdzf3B67rlYFW723X12IBV2sZ3qpu385QC/Vss9wMOtLeNiEQBqcBRZ7kf8DJwq6ru8HtaY1rol5bA83dfSP+0BG57ciVLPjzkdiQTZlSVBxZsZs3eSn59w1hye3f/6VQ3C8hKYJiIDBKRGOAmYGGrbRYCs5zn1wHLVFVFpAfwKvA9VV3ebYmNaSEzJY75sy8gr3cyX3lmNQvXt/77xxj/eeb9vcxfuY+vTRri2qyarhUQVW0E7sFzB9UW4HlV3SwiPxaR6c5mjwMZIlIC3Afc76y/BxgK/FBE1jkPm9zadLu0xBieuXMC5w1M45vz1/Lmvga3I5kw8MGuo/xo4WYm52Zy37QRruVw9VYhVX0NeK3Vuv9q8bwWz7WS1u/7CfATvwc0xgvJcdHMu208s59exbwPK7inuo4MP3fgMuHrQOUpvvrsagakJ/DQjfl+u5HHG9al2RgfiI+J5DuX59KsULTtsNtxTIiqbWji7qdXU9vQzJxbx5Ea7+74bFZAjPGR0X1T6BErLN1qF9SN76kq33tpIxv3V/HQjfkMzUx2O5IVEGN8RUQY2yuSt7ZXdDgqsjGdMXf5bl5eu5/7pg1n2sgst+MAVkCM8an8zEiq6xr5YNdRt6OYELK8pIKfvbaFy0dlcc+koW7H+YgVEGN8aGRGJLFRESzZYqexjG/sO1rD1/6yhsE9E/n1DfkBNSWDFRBjfCg2Upg4tCdLtx46PVqCMZ1WU9/IXU+torlZefTWApICbIw9KyDG+NiUvEz2HT1FSXm121FMEFNVvv3iBrYfOsHvv3QeOT0T3Y70KVZAjPGxybmePq1Lt5a7nMQEs0fe3MGrGw7ynStyuXR4L7fjtMkKiDE+1ic1nlHZKSy16yCmk97YVs4vF23j6rHZ3H3JYLfjtMsKiDF+MCU3k9V7jnHsZL3bUUyQ2VVxkm/8dS15vVP4xbVjAnoCOisgxvjBlLwsT6/07XYay3ivus5z0TwqQvjzLeOIj4k885tcZAXEGD84p28qvZJjWbrFCojxTnOzcu9z69hVcZI/fOk8+qcnuB3pjKyAGOMHERHC5BGZvLn9MA1N1ivdnNnvlhWz+MNDfP/KPC4a2tPtOF6xAmKMn0zOy+REbSMrd1uvdNOxf20u4+ElxVx7Xj9um5jjdhyvWQExxk8uHtqTmKgIO41lOlR86AT3PreOsf1S+ek1owP6onlrVkCM8ZPE2CguHJzBMusPYtpRdaqB2U+vJj4mkj/dMo646MC+aN6aFRBj/GhqXia7Kk6y47D1Sjef1NSsfHP+WkqP1fDIzePokxrvdqSzZgXEGD+a5PRKX2ansUwrv/7XNoq2HebB6aM4Pyfd7TidYgXEGD/ql5ZAbu9kG53XfMIrGw7wx6IdfHH8AGZOGOh2nE6zAmKMn03Jy2TVnmNU1TS4HcUEgNJjNXz7hQ2MG5jGg9NHuh2nS6yAGONnU/KyaGpW65VuAPjDGyU0NSu/++K5xEYF10Xz1qyAGONnY/v1ICMxxm7nNew7WsMLq0q5aXx/+vYIvovmrVkBMcbPIiOESbmZFG0rp9F6pYe1/1tWQkSE8NXCwJmWtivaLSAicvwMjxMisr07wxoTrKbmZXK8tpFVe465HcW4ZO+RGl5cU8qXxg+gd2qc23F8oqMWyA5VTengkQyc7K6gxgSzi4f1IiYywjoVhrHfLysmKkL4t8IhbkfxmY4KyLVevN+bbYwJe0mxUUwYnG6384ap3RUneWntfr40YQBZKaHR+oAOCoiq7hSRSBFZ0tE2/ollTOiZkpvJzsMn2VVhDfdw8/tlJZ7Wx6Wh0/qAM1xEV9UmoEZEUrspjzEha0peFoBNdRtmdlWc5OW1pdx8wUAyQ6j1ARDlxTa1wEYRWUyLax6q+g2/pTImBPVPT2B4VhLLtpZz52cCd55r41u/W1pMTFQEXwmx1gd4V0BedR7GmC6akpfFo2/t5HhtAylx0W7HMX6243A1C9bt587PDKZXcqzbcXzOmwJyBHhNVe0GdmO6aEpuJo8U7eCt7Ye5aky223GMn/1uaTGxUZHMviQ0W5zedCS8CSgWkV+ISJ4vdy4iV4jINhEpEZH723g9VkSec15/X0RyWrz2PWf9NhG53Je5jPGXcwekkZYQbb3Sw0BJ+QkWrj/ArRcNpGdS6LU+wIsCoqo3A+cCO4AnRGSFiMwWkeSu7FhEIoE/AJ8FRgJfFJHWI4vdARxT1aHAQ8D/Ou8diaewjQKuAP7ofJ4xAS0yQpg0IpM3tpXT1KxuxzF+9NulJcRHR3L3JaF37eM0r4YyUdXjwN+A+UAf4BpgjYh8vQv7Hg+UqOpOVa13PntGq21mAPOc5y8CU8Qz3+MMYL6q1qnqLqDE+TxjAt6UvCwqaxpYs9d6pYeq7YdO8MqGA8y6KIf0xBi34/jNGa+BiMjVwO3AEOBpYLyqlotIArAF+H0n990X2NdiuRSY0N42qtooIlVAhrP+vVbv7dtO/tnAbICsrCyKioo6Gbd7VVdXB03WsxXuxyYNSqTAE4tWcXJEcH25hPvvzlt/XFdLbASMijhIUVGZTz6zK/z1e/PmIvr1wEOq+lbLlapaIyK3d2Hfbc0c37pN39423rzXs1J1DjAHoKCgQAsLC88ionuKiooIlqxny44Nnt3zHsXH6ygsvNT/oXzIfndntq3sBCsXvcVXC4dw1WW5XQ/mA/76vXU0mOIiEbkX+Fnr4nGaqi7twr5Lgf4tlvsBB9rbRkSigFTgqJfvNSZgTc7Nori8mr1HatyOYnzst0u3kxgTxV1h0Neno2sgs4BjwIMiskZEHhGRGSKS5KN9rwSGicggEYnBc1F8YattFjo5AK4DlqmqOutvcu7SGgQMAz7wUS5j/G5qnmeu9KVbrVd6KPnwwHFe21jG7RNz6JEQXKcnO6OjsbDKVPVJVb0JKACeAsYBi0RkiYh8pys7VtVG4B5gEZ5rKc+r6mYR+bGITHc2exzIEJES4D7gfue9m4HngQ+BfwJfc4ZdMSYoDMxIZEivRLudN8T8dul2kuOiuOPi0G99gHfXQHA6Ea5wHv8lIj2BLve9UNXXgNdarfuvFs9r8VyDaeu9PwV+2tUMxrhlal4Wc5fv4kRtA8nWKz3obT5QxaLNh/jmlGGkJoTH7/OMt/E6p5geEpGXRGShiCwE5qrqs92Qz5iQNSUvi4Ym5e3iCrejGB94eEkxyXFR3H7xILejdBtvWiB/x3MqaSFgw5kY4yPnDehBarynV/qV5/RxO47pgk37q1j84SHunTqc1PjwaH2Al6Pxqurv/J7EmDATFRnBpBG9PuqVHhnR1t3pJhg8vGQ7KXFR3HZxjttRupU3PdF/KyIPiMiFInLe6YffkxkTBibnZXH0ZD3r9lW6HcV00obSSpZsKeeuzwwOuxGWvWmBnAPcAkzm41NY6iwbY7rg0uG9iIwQlm45xLiBaW7HMZ3w8JJieiRE8+WJOW5H6XbeFJBrgMHOeFXGGB9KjY/m/Jw0lm0t5ztXBEavZeO9dfsqWba1nG9fPiIs76Tz5hTWeqCHv4MYE66m5mWxtewEpcesV3qweWjxdtISopl1UY7bUVzhTQHJArY6Q5ssbHErrzHGBybnenqlL9tqnQqDyeo9x3hz+2FmXzKEpFivutSFHG+O+gG/pzAmjA3ulcTgnoks2VLOrRfmuB3HeOnhJdtJT4zh1gsHuh3FNd5MKPUmsA3PQIYpwDZnnTHGRybnZvLejiOcrGt0O4rxwuo9R3m7uIK7LxlMYpi2PsC7nuh34hmo8At4BjR8r4vDuBtjWpmSl0V9U7P1Sg8SDy0upmdSDLeEcesDvLsG8m3gXFX9sqrOwjOg4nf9G8uY8FKQk0ZyXBTLbHTegLdy91HeKangK5cOISEmfFsf4F0BKQVOtFg+wSdnEjTGdFF0ZASFIzJZtvUwzTZXekB7aPF2eibFMnNCeLc+wLsCsh94X0QeFJEH8EwlWyIi94nIff6NZ0z4mJKbSUV1HRv2V7kdxbTj/Z1HeHfHEf6tcAjxMZFux3GdNwVkB54BFU//WbQAOAgkOw9jjA8UjuhFhMDSLXYaK1A9tGQ7mcmxzJwwwO0oAeGMJ/BU9UcAIpLiWdQTZ3iLMaYTeiTEUDAwnSVbyvnWZSPcjmNaeXdHBe/tPMoDV48kLtpaH+DdXVgFIrIR2ABsFJH1IjLO/9GMCT9T8jLZcvA4BypPuR3FtKCqPLy4mKyUWL443lofp3lzCmsu8FVVzVHVHOBrwBN+TWVMmJqSlwXAUuuVHlDe3XGED3Yf5auFQ6310YI3BeSEqr59ekFV3+GTd2UZY3xkSK9EBmYksMyugwQMVeWhxdvpnRLHjef3dztOQPGmgHwgIn8WkUIRuVRE/ggU2bwgxvieiDAlN4vlO45QU2+90gPBOyUVrNpzjK9NttZHa970gsl3frYeE+sibF4QY3xuSl4mc5fvYnnJEaaNzHI7Tlg73frITo3jhoJ+bscJON7chTWpO4IYYzzOz0knOTaKpVsOWQFx2VvFFazZW8lPrxlNbJS1Plpr9xSWiFx1pjd7s40x5uzEREVwyfBeLNtabr3SXXS69dG3RzzXj7NrH23pqAXySxHZD0g7rwvwM+AVn6cyJsxNycvk1Y0H2XSgijH9bD43NxRtP8y6fZX8zxfOISbKm8vF4aejAnII+M0Z3l/swyzGGEfhiEynV3q5FRAXnG599EuL57pxdu2jPe0WEFUt7MYcxpgW0hNjOG9AGku3HuLeacPdjhN2lm0tZ0NpFf977TlER1rroz32X8aYADU5L5NN+49TVlXrdpSwoqo8vKSYAekJfOE8a310xAqIMQFqqtMr3eZK715LtpSzcX8VX5881FofZ2D/dYwJUMMyk+iXFm+TTHUjT+tjOzkZCVxzbl+34wS8dq+BiMgXOnqjqr7k+zjGmNNEhKl5WcxfuZfahibrBd0N/vXhITYfOM6vrx9LlLU+zqij/0JXO487gMeBmc7jMeBm/0czxkzOzaS2oZnlJTZXur81O9c+BvVMZEZ+tttxgkK7BURVb1PV2/AMVzJSVa9V1WuBUV3dqYiki8hiESl2fqa1s90sZ5tiEZnlrEsQkVdFZKuIbBaRn3c1jzGBasLgdBJjIm103m6w5lATWw4e5xtThlrrw0ve/FfKUdWDLZYPAV29r/B+YKmqDgOWOsufICLpeMbfmgCMBx5oUWh+paq5wLnARBH5bBfzGBOQYqMi+cywXizbUo6q9Ur3l+Zm5e8l9Qzulcj0sXbtw1veFJAiEVkkIl92WgGvAm90cb8zgHnO83nA59vY5nJgsaoeVdVjwGLgClWtUdU3AFS1HlgD2L12JmRNycuk7Hgtmw8cdztKyHp9Uxml1co3pwwjMqK9wTdMa+LNXzXOBfXPOItvqerLXdqpSKWq9mixfExV01pt8x9AnKr+xFn+IXBKVX/VYpseeArIVFXd2c6+ZgOzAbKyssbNnz+/K9G7TXV1NUlJSW7H8As7trNzvE755hs1fH5oNDOGxvj0s89WKP7uVJX/ereWhsYmfnZJIhESegWkq7+3SZMmrVbVgtbrvRnO/fQdV2d115WILAF6t/HS9739iLaitPj8KOCvwO/aKx4AqjoHmANQUFCghYWFXu7eXUVFRQRL1rNlx3b2ntixnJ21SmHhxT7/7LMRir+7d0sq2Hfife4YHcvkSaE5+Li/fm/ezIn+BecidpWIHBeREyJyxra0qk5V1dFtPBYAh0Skj/P5fYC2rhCWAi2HwOwHHGixPAcoVtWHz5TFmGA3NS+L9aVVlB+3Xum+Nnf5LnomxTChj1d/T5sWvLkG8gtguqqmqmqKqiarakoX97sQmOU8nwUsaGObRcBlIpLmXDy/zFmHiPwESAX+vYs5jAkKk3MzAXhjm92N5Uu7Kk6ydGs5MycMJCYy9E5d+Zs3BeSQqm7x8X5/DkwTkWJgmrOMiBSIyGMAqnoU+G9gpfP4saoeFZF+eE6DjQTWiMg6EbnTx/mMCSi5vZPp2yOeJVusgPjSvHd3Ex0RwcwLBrgdJSh502ZbJSLPAX8H6k6v7EpPdFU9AkxpY/0q4M4Wy3OBua22KaX9OUqMCUkiwuTcTF5cXWq90n2k6lQDz6/ax9Vjs8lMjuNDtwMFIW9aIClADZ5TSKd7p9tMhMZ0syl5mZxqaGLFziNuRwkJz6/cR019E7dNzHE7StDyZk7027ojiDGmYxcMziAhJpJlW8qZNCLT7ThBrbGpmSff3c2EQemM7pvqdpyg5c1dWMNFZKmIbHKWx4jID/wfzRjTUlx0JBcP7cnSLYesV3oXLf7wEPsrT3H7xYPcjhLUvDmF9SjwPaABQFU3ADf5M5Qxpm1T8jI5UFXL1rITbkcJanOX76J/evxHc66YzvGmgCSo6get1jX6I4wxpmOTnNt5l26xOUI6a0NpJSt3H+PLFw2yYUu6yJsCUiEiQ3B6gYvIdcDBjt9ijPGHzOQ4xvZLtdF5u+CJ5btJio3ihgIbQq+rvCkgXwP+DOSKyH48nfe+4tdUxph2TcnLYt2+Siqq6868sfmEQ8dreWXDAa4v6EdyXLTbcYJehwVERCKAAlWdCvQCclX1YlXd0y3pjDGfMjk3E1WbK70znnlvD43Nym0X2cVzX+iwgKhqM3CP8/ykqtqVO2NcNio7hd4pcSyzXulnpbahiWff38u0vCwGZCS4HSckeHMKa7GI/IeI9HdmEkx3JnsyxrhARJicl8nbxYepa2xyO07Q+Pva/Rw9WW+37vqQNwXkdjzXQd4CVjuPVf4MZYzp2NS8TE7WN/H+zqNuRwkKqsrc5bsY2SeFCYPs719f8aaA5KnqoJYPPAMZGmNcctGQnsRFR9jtvF5aXnKE7Yequf3iQUgIThjlFm8KyLterjPGdJO46EguHd6L1zaV0dRsvdLP5PScH1eP7eN2lJDSbgERkd4iMg6IF5FzReQ851EI2BUoY1w2I78vh0/U8Z4NrtihnYerWba1nJsvGEhslI1i7EsdDaZ4OfBlPDMB/pqPh1A/Dvynf2MZY85kcm4mSbFRLFx3gIlDe7odJ2A9+e5uYiIjmDlhoNtRQk67BURV5wHzRORaVf1bN2YyxnghLjqSy0Zl8dqmg/z486Psr+s2VNU08MKqUqbnZ9MrOdbtOCHnjNdArHgYE7hm5PflRG0jRdsOux0lID23ai+nGmzOD3/x5iK6MSZATRySQUZiDAvXHXA7SsBpbGpm3rt7uGBwOqOybc4Pf7ACYkwQi4qM4HNj+rBkyyGq62yQ7Jb+dXrOj4nWcdBfzjQWVoozEm/r9WP8F8kYczZm5GdT19jMvzaXuR0loDz+zi4GpCcwxeb88JuObuO9AdgK/E1ENovI+S1eftLfwYwx3jlvQBr90uJZYKexPrJuXyWr9xzjtok5NueHH3XUAvlPYJyq5gO3AU+LyBec1+w3YkyAEBGmj83mnZIKjtgQ7wA8sXwXybFRXF/Q3+0oIa2jAhKpqgcBnBkJJwHfF5Fv4EwuZYwJDNPzs2lqVl7baHO9lVXV8uqGg9xwfn+SYjvq6ma6qqMCcqLl9Q+nmEwCZgCj/B3MGOO93N4pjMhKttNYwNPv7aZZlS9flON2lJDXUQH5t9avq+px4Ao8I/QaYwLI9PxsVu05xr6jNW5Hcc2p+ib+8v5epo3Mon+6jbjkb+0WEFVdr6rFIjJIRK4Skc+JyGBVbVDVZ7szpDHmzKaPzQbgHxvCtxXy93X7OVbTYLfudpOO7sJKEZHngaV4Whx3AktE5AURSemugMYY7/RPT+C8AT3CtlOhqjL3nV2Myk5hvM350S06OoX1O+BDYKiqflZ21YgAABI5SURBVEFVrwGGABuB/+uOcMaYszMjvy9by06wrSz8Zp9+p6SC4vJqbp9oc350l44KyERVfdCZFx0A9fgxcKH/oxljztaV5/QhMkJYuH6/21G63dx3dtEzKZarbM6PbtNRAbESbkyQ6ZUcy0VDMli4/gCq4XO3fUl5NW9sO8wtNudHt+qogCwXkf+SVm1BEfkh8J5/YxljOmtGfl/2HT3F2n2VbkfpNk++u8sz58cFA9yOElY6KiBfB84BSkTkbyLyoojsAMY6rxljAtDlo7KIiYoIm4vplTX1/G31fmbkZ9Mzyeb86E4d3cZ7XFWvBy7DM/bVU8BlqnqdqnbpTxsRSReRxSJS7PxMa2e7Wc42xSIyq43XF4rIpq5kMSbUJMdFMzUvk1c2HKCxqfnMbwhy81fuc+b8sFt3u1tHt/FeLiLXqeoOVf2Hqi5U1R0iMlNEpnVxv/cDS1V1GJ7bhO9vY//pwAPABGA88EDLQuOMy1XdxRzGhKTpY7OpqK5nRYjPl97Q1My8d3dz0ZAMRmZb74Lu1tEprB8Bb7axfinw4y7udwYwz3k+D/h8G9tcDixW1aOqegxYjKcXPCKSBNwH/KSLOYwJSYUjMkmOjQr5oU0WbS7jYFWtdRx0SUcjjSWo6qfmyVTVMhFJ7OJ+s1oM1HhQRDLb2KYvsK/FcqmzDuC/gV8DZxyzQURmA7MBsrKyKCoq6kLs7lNdXR00Wc+WHVv3yO8Jr6wr5bL0o8RE+uamykA6PoCH3ztFZoIQcehDisq3dOmzAu3YfMlfx9ZRAYkTkShV/cQ0ZyISDcSf6YNFZAnQu42Xvu9ltrb+xauI5OPp3HiviOSc6UNUdQ4wB6CgoEALCwu93L27ioqKCJasZ8uOrXtE9j3M249/QFNmLoXn+KZvRCAd39q9xyj557s8ePVIJvugBRJIx+Zr/jq2jk5hvQQ82rK14Tz/k/Nah1R1qqqObuOxADgkIn2cz+wDlLfxEaVAy8H8+wEH8HRiHCciu4F3gOEiUnSmPMaEmwsHZ9AzKZaF60PzNNYTy3eTHBvFdTbnh2s6KiA/AA4Be0RktYisAXYDh53XumIhcPquqlnAgja2WQRcJiJpzsXzy4BFqvqIqmarag5wMbBdVQu7mMeYkBMVGcFVY/qwdGs5x2sb3I7jUwerTvHaxoPcaHN+uKqj23gbVfV+PK2AL+P5oh+gqveralf/Nf4cmCYixcA0ZxkRKRCRx5z9H8VzrWOl8/ixs84Y46UZ+dnUNzazaFNozZf+1Io9NKsyy+b8cJU3pVvxfMlfjOcaxDvAI6pa29mdquoRYEob61fhGfX39PJcYG4Hn7MbGN3ZHMaEuvz+PRiQnsDC9QdCZnrX03N+XDayt8354bKOTmGd9hSeGQh/j2cU3jzgaX+GMsb4xun50peXVHD4RGjMl/7S2lKqTjVw+8V2667bvCkgI1T1DlV9w3nMBob7O5gxxjdm5GfTrPBqCEw01dzsmfNjdN8Uzs9pcwAL0428KSBrReSC0wsiMgFY7r9IxhhfGpaVTF6fFBaEwN1Yb5dUsOPwSZvzI0B4U0AmAO+KyG7n1tkVwKUislFENvg1nTHGJ6aPzWbt3kr2Hgnu+dLnvrOLXsmxfG6MzfkRCLwpIFcAg4BLnccg4ErgKuBq/0UzxvjK1c4kS8E8X3pJ+Qne3G5zfgSSM96Fpap7uiOIMcZ/+qUlcH5OGn9fu5+vFg4JytM/TyzfTUxUBDMn2JwfgcKbFogxJgRMH5tNcXk1W4NwvvTKmnr+tqaUa/L7kmFzfgQMKyDGhImP50sPvtNYf/1gH7UNzdx2cY7bUUwLVkCMCRMZSbF8ZlhPFq47QHNz8MyXfnrOj4lDM8jtbXN+BBIrIMaEkRn52eyvPMWavcfcjuK11zeVUXbc5vwIRFZAjAkj00b2JjYqIqgmmpr7zi5yMhKYNKKtaYOMm6yAGBNGkmKjmDoyi9c2HqQhCOZLX7P3GOv2VXLbxEFERATfnWOhzgqIMWFmxthsjpysZ3lJhdtRzmjuO7tIjoviunH93I5i2mAFxJgwc+mIXqTERbEwwE9jHag8xeubyrjp/P4k2pwfAckKiDFhJjYqks+O7sOizWWcqm9yO067nlqxB1Xl1gtz3I5i2mEFxJgwNCM/m5P1TSzb2tZs0u6rqW/krx/s5fJRNudHILMCYkwYmjA4g8zkWBas2+92lDY9sXy3zfkRBKyAGBOGIiOEq8ZkU7TtMFU1gTVf+htby/n1v7ZxxajeFAy0OT8CmRUQY8LUjPxs6pua+efmg25H+cj2Qyf4+l/Xkts7hd/cODYoB30MJ1ZAjAlTY/qlkpOREDBjYx2pruOOeSuJj4nksVkFJMTYnVeBzgqIMWFKRJie35d3dxyh/Hitq1nqGpv4yjOrKT9ex6O3FpDdI97VPMY7VkCMCWPTx2ajCv/Y4N5pLFXlBy9vYuXuY/zy+rHk9+/hWhZzdqyAGBPGhmYmMSo7xdXTWI++vZMXVpfyjSnDmD4227Uc5uxZATEmzM3Iz2b9vkp2V5zs9n0v+fAQ//P6Vj53Th/+fcqwbt+/6RorIMaEuavHZiNCt7dCthw8zjfnr+Wcvqn86vqxNlhiELICYkyY65Maz/k56SxYtx/V7plo6vCJOu6ct4qkuCgevbWA+JjIbtmv8S0rIMYYZuRns+PwST48eNzv+6ptaOLup1dx5GQdj916PlkpcX7fp/EPKyDGGK4c3YeoCPH7CL2qyvde2siavZX85oZ8zumX6tf9Gf+yAmKMIS0xhkuG92Lhev/Ol/7Hoh28vHY/35o2nCvP6eO3/ZjuYQXEGAN4TmMdrKpl5e6jfvn8f246yC8XbWP62GzumTzUL/sw3csKiDEGgKl5WcRHR/rlbqxN+6u497n15PfvwS+uG2NjXIUIVwqIiKSLyGIRKXZ+tjnkpojMcrYpFpFZLdbHiMgcEdkuIltF5NruS29MaEqMjWLayCxe3XiQ+kbfzZdefryWu55aRVpCNHNuHUdctN1xFSrcaoHcDyxV1WHAUmf5E0QkHXgAmACMBx5oUWi+D5Sr6nBgJPBmt6Q2JsRNH5tNZU0D75Qc9snn1TY0cdfTq6msaeDRWQVkJtsdV6HErQIyA5jnPJ8HfL6NbS4HFqvqUVU9BiwGrnBeux34HwBVbVbVCj/nNSYsXDK8F6nx0Szwwd1Yqsq3X9zA+n2VPHxTPqOy7Y6rUCPd1XHoEzsVqVTVHi2Wj6lqWqtt/gOIU9WfOMs/BE4BjwEbgReAQmAHcI+qHmpnX7OB2QBZWVnj5s+f7/sD8oPq6mqSkpLcjuEXdmyB7clNdaw42MjvJiUQG/XJaxVnc3wLSup5uaSB64ZHc9XgGH9E9alQ+N21p6vHNmnSpNWqWtB6vd8G3BeRJUDvNl76vrcf0cY6xZO5H7BcVe8TkfuAXwG3tPUhqjoHmANQUFCghYWFXu7eXUVFRQRL1rNlxxbY4gYcoWjOe9T1GsHlrQY39Pb4Xt1wkJdL1vCFc/vyyxuCY2KoUPjdtcdfx+a3AqKqU9t7TUQOiUgfVT0oIn2A8jY2K8XTwjitH1AEHAFqgJed9S8Ad/giszEGxuek0zsljoXr9ndqdNwNpZV864V1jBuYxv9ce05QFA/TOW5dA1kInL6rahawoI1tFgGXiUiac/H8MmCRes65/YOPi8sU4EP/xjUmfERECFeP7cOb2w9TWVN/Vu8tq/LccZWRGMufbxlHbJTdcRXK3CogPwemiUgxMM1ZRkQKROQxAFU9Cvw3sNJ5/NhZB/Bd4EER2YDn1NW3ujm/MSFtRn5fGpqU1zeVef2eU/VN3PXUKqprG3lsVgE9k2L9mNAEAlcmHVbVI3haDq3XrwLubLE8F5jbxnZ7gEv8mdGYcDYqO4XBPRNZsG4/Xxw/4IzbNzcr33phHZsOVPHoLQXk9UnphpTGbdYT3RjzKZ750rN5f9dRyqrOPF/6w0uLeW1jGd/7bC5TR2Z1Q0ITCKyAGGPadHq+9Fc2dNwnZMG6/fxuaTHXj+vHXZ8Z3E3pTCCwAmKMadPgXkmM6ZfaYafCtXuP8e0XNzA+J52fXDPa7rgKM1ZAjDHtmj42m437q9hxuPpTrx2oPMXsp1eTlRLLIzefZ3dchSErIMaYdl01xpkvvVUrpKa+kTvnreJUfROPzzqfDLvjKixZATHGtKt3ahwXDMrgH+sPfDRfenOzcu9z69hadpzff+lchmclu5zSuMUKiDGmQzPys9lZcZJN+z3zpf968TYWbT7E9z83kkkjMl1OZ9xkBcQY06HPju5DdKSwYN1+3j3QyB/e2MEXx/fn9ok5bkczLnOlI6ExJnikJkRz6fBMXlxTSvWpBi4YnM6PptsdV8ZaIMYYL8zI90w0lR4vPDJzHDFR9tVhrAVijPHCZaOy+PrkofRt2E9aYuDP7WG6h/0ZYYw5o9ioSL512Qh6J9pXhvmY/WswxhjTKVZAjDHGdIoVEGOMMZ1iBcQYY0ynWAExxhjTKVZAjDHGdIoVEGOMMZ1iBcQYY0ynyOkhmsOBiBwG9ridw0s9gQq3Q/iJHVvwCuXjs2Nr30BV7dV6ZVgVkGAiIqtUtcDtHP5gxxa8Qvn47NjOnp3CMsYY0ylWQIwxxnSKFZDANcftAH5kxxa8Qvn47NjOkl0DMcYY0ynWAjHGGNMpVkCMMcZ0ihWQACIi/UXkDRHZIiKbReSbbmfyNRGJFJG1IvKK21l8TUR6iMiLIrLV+R1e6HYmXxGRe51/k5tE5K8iEud2pq4QkbkiUi4im1qsSxeRxSJS7PxMczNjZ7VzbL90/l1uEJGXRaSHL/ZlBSSwNALfUtU84ALgayIy0uVMvvZNYIvbIfzkt8A/VTUXGEuIHKeI9AW+ARSo6mggErjJ3VRd9iRwRat19wNLVXUYsNRZDkZP8uljWwyMVtUxwHbge77YkRWQAKKqB1V1jfP8BJ4voL7upvIdEekHfA54zO0sviYiKcAlwOMAqlqvqpXupvKpKCBeRKKABOCAy3m6RFXfAo62Wj0DmOc8nwd8vltD+Uhbx6aq/1LVRmfxPaCfL/ZlBSRAiUgOcC7wvrtJfOph4DtAs9tB/GAwcBh4wjlF95iIJLodyhdUdT/wK2AvcBCoUtV/uZvKL7JU9SB4/pgDMl3O4y+3A6/74oOsgAQgEUkC/gb8u6oedzuPL4jIVUC5qq52O4ufRAHnAY+o6rnASYL3FMgnONcCZgCDgGwgUURudjeV6QwR+T6eU+XP+uLzrIAEGBGJxlM8nlXVl9zO40MTgekishuYD0wWkWfcjeRTpUCpqp5uMb6Ip6CEgqnALlU9rKoNwEvARS5n8odDItIHwPlZ7nIenxKRWcBVwEz1UQdAKyABREQEzzn0Lar6G7fz+JKqfk9V+6lqDp4LsMtUNWT+ilXVMmCfiIxwVk0BPnQxki/tBS4QkQTn3+gUQuQGgVYWArOc57OABS5m8SkRuQL4LjBdVWt89blWQALLROAWPH+dr3MeV7odynjt68CzIrIByAd+5nIen3BaVS8Ca4CNeL43gnrYDxH5K7ACGCEipSJyB/BzYJqIFAPTnOWg086x/R+QDCx2vlf+5JN92VAmxhhjOsNaIMYYYzrFCogxxphOsQJijDGmU6yAGGOM6RQrIMYYYzrFCogxPiIiKiJPt1iOEpHDnR152Bnd96stlgtDcRRjE7ysgBjjOyeB0SIS7yxPA/Z34fN6AF8941bGuMQKiDG+9TqeEYcBvgj89fQLznwTf3fmZHhPRMY46x905nAoEpGdIvIN5y0/B4Y4Hb9+6axLajHnyLNOz3BjXGEFxBjfmg/c5Ey4NIZPjqb8I2CtMyfDfwJPtXgtF7gcGA884IyJdj+wQ1XzVfXbznbnAv8OjMQzAvBEfx6MMR2xAmKMD6nqBiAHT+vjtVYvXww87Wy3DMgQkVTntVdVtU5VK/AM4pfVzi4+UNVSVW0G1jn7MsYVUW4HMCYELcQzf0YhkNFifVunm06PJVTXYl0T7f+/6e12xvidtUCM8b25wI9VdWOr9W8BM8FzRxVQcYb5Xk7gGQDPmIBkf70Y42OqWopnfvTWHsQzY+EGoIaPhw5v73OOiMhyEdmE5+L8q77OakxX2Gi8xhhjOsVOYRljjOkUKyDGGGM6xQqIMcaYTrECYowxplOsgBhjjOkUKyDGGGM6xQqIMcaYTvl/Ip0a2djxHJkAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"mon_trends = np.polyfit(years,month_anom,1)[0,:]\n",
"plt.plot(month,mon_trends)\n",
"plt.xlabel('Month'); plt.ylabel('pCO2 trend [ppm/yr]'); plt.grid()"
]
},
{
"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.8"
}
},
"nbformat": 4,
"nbformat_minor": 2
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment