# Contents * cc {:toc} ## Idea Create some plots related to the [[Carbon cycle]] and do some simple analysis using [[Sage]]. ## Details ### Mauna loa ### Diagrams The classic Keeling curve, plotted upto jan 2011 in Sage with default aspect ratio for Sage List_plot: <div align = "center"> <img src = "http://www.azimuthproject.org/azimuth/files/mloaco2.png"> </div> Here is a plot of the whole series but showing season corrected values : (also default aspect ratio) <div align = "center"> <img src = "http://www.azimuthproject.org/azimuth/files/mloaco2trend.png"> </div> I think it can be better with setting it default aspect ratio to one, so we'll try that as well and add some major grid lines: <div align = "center"> <img src = "http://www.azimuthproject.org/azimuth/files/mloaco2aspect1.png"> </div> Trend - season corrected aslo with SAGE aspect ratio one: <div align = "center"> <img src = "http://www.azimuthproject.org/azimuth/files/mloaco2trendaspect1.png"> </div> ### Draft Code Trying to see if TimeSeries, can be used better instead of this <pre> # based on http://wiki.sagemath.org/interact/web made by Marshall Hampton CC 3.0 import urllib2 as U import time current_year = time.localtime().tm_year co2data = U.urlopen('ftp://ftp.cmdl.noaa.gov/ccg/co2/trends/co2_mm_mlo.txt').readlines() datalines = [] for a_line in co2data: if a_line.find('Creation:') != -1: cdate = a_line if a_line[0] != '#': temp = a_line.replace('\n','').split(' ') temp = [float(q) for q in temp if q != ''] datalines.append(temp) @interact def azimviewer(month = [1,2,6,12,24],trend= True): if trend: sel_data = [[q[2],q[5]] for q in datalines] else: sel_data = [[q[2],q[4]] for q in datalines] lp = list_plot(sel_data, plotjoined=True, rgbcolor=(1,0,0)) lp.show(axes = True, figsize = [7,5]) </pre>