Import COPRA output and export ensemble data (time and proxy)
(N. Marwan, 2020-02-26)
Contents
Import COPRA output
The content of the COPRA output is stored in variable d. All variables can be accessed by d.VARNAME.
filename = 'BB-1_d13C_2015-03-18.mat';
load(filename)
The ensemble of time simulations
The ensemble of time simulations is in T_unc, starting with column 3. These times correspond to the depth values of each proxy measurement. The depth values are in d_prxy and also in the 1st column of T_unc.
depth = d.d_prxy; t_ensemble = d.T_unc(:,3:end);
The ensemble of proxy values on a fixed time axis
The final time axis is 1st column in T_cert. The proxy values are in v_prxy and also in 2nd column of T_unc. These values have to be interpolated using the ensemble of age models in T_unc. We have M realizations.
t = d.T_cert(:,1); proxy_ensemble = zeros(length(t), d.M); for i = 1:d.M proxy_ensemble(:,i) = interp1(d.T_unc(:,i+2),d.v_prxy,t); end
Remark
The ensemble of proxy values on a fixed time axis is a bit different from simply plotting the proxy values vs. the realizations of the time axis. In the first, the proxy values are interpolated to a new time axis, in the latter, the proxy values are untouched, but simply plotted at different time locations. The interpolated, therefore, look different.
subplot(1,2,1) plot(t_ensemble(:,1:100),d.v_prxy,'color',[0 0 0 .1]) xlabel('Time'), ylabel('Proxy') title('Proxy vs. realizations of time axis') subplot(1,2,2) plot(t, proxy_ensemble(:,1:100),'color',[0 0 0 .1]) xlabel('Time'), ylabel('Proxy') title('Proxy interpolated to same time axis')