The Azimuth Project
Experiments in carbon cycle with Sage (Rev #4, changes)

Showing changes from revision #3 to #4: Added | Removed | Changed

Contents

Idea

Create some plots related to the Carbon cycle and do some simple analysis using Sage.

Details

Mauna loa

Diagram Diagrams

Code

The classic Keeling curve, plotted upto jan 2011:


Here is a plot of the whole series but showing season corrected values :

Draft Code

# 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]) #, aspect_ratio=1)