Experiments in carbon cycle with Sage in The Azimuth Project

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