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: ratio for Sage List_plot:

I think it can be better with setting it default aspect ratio to one, so we’ll try that as well.

Here is a plot of the whole series but showing season corrected values : (also default aspect ratio)

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:

Trend - season corrected aslo with SAGE aspect ratio one:

Draft Code

Trying to see if TimeSeries, can be used better instead of this

# 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])