function [n1,edges,total]=disp_extrem_4_cross_speed(hurricane)

rep='wavelength analysis\same_psd_length\unfiltered\extrem';

 mu_total=[];
 sigma_total=[];
path=['C:\Documents and Settings\sylvie lorsolo\My Documents\research\',hurricane,'\data\',rep,'\unwrinkled\half\'];
files=dir([path,'\*.mat']);

%For each trans at each height
total=[];

for fil=1:length(files)
    load([path,files(fil).name]);
    for i=2:length(wave)
        [lignes,colonnes]=size(wave(i).p_to_t);
        total=[total;wave(i).p_to_t];
    end
end
    
edges=0:100:2000;
n=histc(total,edges);
x_axis=50:100:1200;
n=n(1:length(x_axis),:);
n1=n;
mu=[];
sigma=[];
for i=1:colonnes

    u_nzeros=find(n1(:,i) ~=0);
    somme=sum(n1(u_nzeros,i));
    n1(u_nzeros,i)=n1(u_nzeros,i)./somme;
end
