clear all; close all; clc
%This changes your Tpicks.mat file to one comprehensive pick file.

x_final=zeros(1,100);
amp1_final=zeros(1,100);
amp2_final=zeros(1,100);
start_timef_final=zeros(100,4);
picklength=0;
tick=0;

for i=161:174
    if i==164
        i
    elseif i==173
        i
    elseif i==169
        i
    else
        tt=strcat(['load Tpicks_',num2str(i),'.mat']);
        eval([tt])
    end
    amp1_final(1,tick+1:tick+length(amp1))=amp1;
    amp2_final(1,tick+1:tick+length(amp2))=amp2;
    x_final(1,tick+1:tick+length(x))=x;
    for j=1:length(x)
        start_timef_final(tick+j,:)=start_timef(j,:);
    end
        tick=tick+length(x);
        picklength=length(x)+picklength;
end

X_final(1,1:picklength)=x_final(1,1:picklength);
Amp1_final(1,1:picklength)=amp1_final(1,1:picklength);
Amp2_final(1,1:picklength)=amp2_final(1,1:picklength);
Start_time_final(1:picklength,:)=start_timef_final(1:picklength,:);

save start_times.mat 
        

pk2pk=(abs(Amp2_final-Amp1_final));
pk2pk_hist=hist(pk2pk);
bar(pk2pk_hist)
xlabel 'Peak to Peak Amplitude (Pa)';
ylabel 'Number';