
for i=1:5
    figure(i+1)
    plot(nanmean((turb_out(i).turb)),z_axis,coul)
    hold on
    legend([input_data.hurricane(5:length(input_data.hurricane))])
end


% for i=1:5
%    
%      max_bin=12;
% bins=0:max_bin;
% bins=[bins inf];
%      [li_turb,col_turb]=size(turb_out(i).turb);
%     temp=reshape(turb_out(i).turb, li_turb*col_turb,1);
%         u_nan=find(isnan(temp)~=1);
%         p=((histc(temp, bins))/length(u_nan))*100;
%     
% figure
% bar(x_axis,p)
%     contourf(x_axis,z_axis,p',50);
%      shading flat
%  end