% This script runs CORELOGINTEGRATE.m and DATADJUST.m with the data used
% in Gilbert and Burke (2008) as an example.
%
%
% -------------------------------------------------------------------------
%
% Gilbert, L.A. and A. Burke, 2008, Depth-shifting cores
% incompletely recovered from the upper oceanic crust, IODP Hole 1256D,
% G-cubed, 2008gc002010.
%
% Revised 5/8/08, AB and LAG
% for further information contact aburke@mit.edu or lgilbert@williams.edu

% Load the data
load GRA.dat %GRA density data
load log.dat %Logging density data
load pieces.dat %Piece depth information

% This runs corelogintegrate.m for the above data.
[newdepths,GRAbefore,GRAafter] = corelogintegrate(GRA,log,pieces,0.2);

% Load other physical data
load MAD.dat %MAD density data

% This runs datadjust.m for MAD density data 
[MADbefore,MADafter] = datadjust(MAD, newdepths, log);