Scripts and data used in ‘The 405 kyr and 2.4 Myr eccentricity components in Cenozoic carbon isotope records’
climate of the past discussions
This repository contains the R-scripts used to generate the data composite
(scripts/read_data.R and scripts/make_datacomposite.R), as well as the data
sources used from the cited publications:
| file | reference |
|---|---|
data/320-U1334_d18O-d13C.tab |
Beddow et al., 2016 |
data/321-U1337_isotope_benth_foram.tab |
Holbourn et al., 2015 |
data/Pälike_clean.csv |
Pälike et al., 2006 |
data/Tian_U1337_all.csv |
Tian et al., 2013 |
data/Zachos2008.csv |
Zachos et al., 2001, 2008 |
The out/ directory contains the results of the R-scripts.
R-code used to generate the figures was omitted, because it would require substantial cleaning up and commenting before being publication-ready and possibly useful to others.
The LOSCAR source code in C can be obtained from the authors by sending a request to loscar.model@gmail.com.
The LOSCAR model version used here was 2.0.4. To see which changes were made
to the 2.0.4 LOSCAR model, see the patch file (Loscar-2.0.4.patch). In short,
it introduces astronomical forcing as a new input source, as well as noise. To
achieve this, we define some new input paramters (ORBFILE, ORBSTRENGTH, and
NOISELEVEL) and some new global variables for turning on/off the specific
components (NOISE, SCALE, VERBOSE, ADDNOISE, etc.). If you have access
to the Loscar source code, you can apply the patch by running the following comand
from your terminal:
# assuming you have Loscar-2.0.4 in your ~/Downloads
cd ~/Downloads/Loscar-2.0.4/
# download the patch file
wget https://raw.githubusercontent.com/japhir/cp-2018-42/master/Loscar-2.0.4.patch
# copy the emiss.c file to orb.c, our new file
cp emiss.c orb.c
# apply the patch to the whole directory
patch -p1 < ./Loscar-2.0.4.patch
# build the program
make
# make the binary executable
chmod +x loscar.x
# run the program
./loscar.x