This document contains a history of changes made to get the CSIRO Mk3L-COAL model running.
Source code and scripts to run the model can be found in various folders.  For more information
contact Andrew Lenton at CSIRO (Andrew.Lenton@csiro.au).

The original code I got was r1006

I have updated it and got more stuff but I did not get the latest
version of CABLE to work.  Therefore, I have gone back to r1006 version
in the CABLE repository. 

Combined CABLE(r1006) with the latest version of Mk3L with OBGC.
The code is stored in coal_1. 

To do the conversion I used both the CABLE version in repo and Pak version.

The simple test case setup for CABLE is now working with Mk3L 

Scripts used to build and run code 

1. compile 
2. setup_shine_cable
3. run_shine_cable -> in directory where setup put the required files


I will redo the scripts to make the code work with ocean, atm, coupled models

Now testing the new code

1. atmos works setup_shine_cable
2. coupled fails!
3. ocean only works

June 20, 2011

Okay, I have now modified the code and makefile to build the coupled model
with some optimization. Further I have added the scripts to compute the 
change in carbon uptake with time from the ocean and land carbon uptake. 
Finally, I have modified the scripts used to compile, setup and run COAL 
which are:

1. compile_coal
2. setup_coal
3. run_coal

From a coupling perspective I may break the atmCO2conc.f90 routine into two,
one to calculate ocean and land carbon uptake and one to calculate the atmos
concentration.  Now thinking about it I could just modify the code so that 
if the land or ocean output is missing it just sets the value to zero. 
I like this idea and I will do this.

July 10, 2011
New version of cable added to the repository. The code compiles and runs 
but additional files are needed to make it work with the new version of 
cable (cable_jul_qian  directory).  I have tested the code and it ran for 
a year.  I have copied the cspin1 restart files for the climate model into
this directory for the restart file and I will run the coupled system with
the new land-surface scheme to see how it affects the climate system.

July 12, 2011 
New pre-industrial run started from my cspin1 simulation.  I would like 
to run for 100 years to assess the potential drift. For this run there is
no change in the atmospheric co2, it is held constant at 280.

Made a fix to the CABLE to use atmospheric co2 rather than a fixed value.

Have no run the coupled model for 100 years with the new land scheme but
no carbon interaction.

July 18, 2011
I notice the ocean air-sea gas fluxes were not getting the wind and ice
information.  The changes to radin.f and icedriver.f that I made in 
mk3l v1.2 had not been incorporated into COAL.  I forgot I had
changed these files.  I have now add the information back and I will
continue the couple spinup.

I have started another 100 year run with the corrected air-sea fluxes.  
I start the OBGC from the end of the last run.  Run has finished and
the OBGC code is now working properly.  The coupled system
could do with a longer spinup.