function [time,d_data] = GetData(dataname);
%  Extract and view data from .mat files (for example data1.mat)
%  For example   [time,d_data] = GetData

load(dataname);

%  Extract time and displacement data
%  I'm not quite sure which rows we want...  try these

time = trace_x(1,:)';     % observation times
d_data = trace_y(4,:)';   % displacement
N = length(time);

figure(1)
plot(time,d_data)
xlabel('Time (s)')
ylabel('Displacement (m?)')


%  truncate data to start at max displacement
[maxd,index] = max(d_data)
time = time(index:N,1) - time(index,1);
d_data = d_data(index:N,1);


%  center the truncated displacement data at the average
d_avg = mean(d_data)
d_data = d_data - d_avg;


%  plot new data
figure(2)
plot(time,d_data)
xlabel('Time (s)')
ylabel('Displacement (m?)')