README


Brovkin et al., 2016 QSR

This README explains links and scripts to repeat the analysis the paper is based on.

The paper is basiclaly showing model results based on three different models:

- CLIMBA (Climber + JSBACH)
- Genie
- Bern3D

and these data are compared to timeseries of reconstruced CO2 and 13C data.




co2_data:
--------------------------------------
The timeseries of modleld and reconstructed atmospheric CO2 and 13C


To repeat simulations of Genie, please contact Andy Ridgwell
To repeat simulations of Bern3D, please contact Fortunat Jost
To repeat simulations of JSBACH, please contact Tim Bruecher	


runing the models:
--------------------------------------

The data-output is saved within tar files at DKRZ.
The relevant data of Genie, Bern3D plus scripting and data of JSBACH

to download data, please use the archive system at DKRZ and navigate to the following system (via pftp)

pftp
cd mj0060/m300107/experiments/

There the data for the different mpdels and experiments can be found.

JSBACH war running on blizzard - a NEC SX6 supercomputer, that was shut down in spring 2016 and is replaced by MISTRAL, a Intel Cluster.
Both, JSBACH and Climber can be compiled and run on the new supercomputer, too.



Original model outut data 
--------------------------------------

model simulations of JSBACH and CLIMBER have been done by Tim Brücher using the account (m300107) at DKRZ

data are stored for each model and for each experiment seperately:




data anlysis:
--------------------------------------

cdo is used to merge the datasets of the models and to come up finaly with a timeseries of one variable for the full transient time series.

cdo is also used to calculate annual means . sums

The final processing of the data (like changes in the carbon) is calculated during plotting within ncl

The plotting routines are based on the default NCL libs (6.2) and a private lib to define colours, design and process data.



PLOTTING and Calculations:
--------------------------------------
-- final_bern.ncl
plotting data out of the BERN3d model

-- final_genie.ncl
plotting data out of GENIE model

-- paper.ncl:
Tihs paper is based on ncl scripting language to plot an process data (nc. format)
The script contains different subroutines for each figure and calculates data for the table2 Figure 1 is produced on data (Jplease contact ochen Schmidt)

there are some additional files (within libs and in the plotting folder) that are partly used to do the calculations