function [perct,data, x_axis, z_axis]=per_radius_disp(data,percent, x_axis,z_axis,c)


for i=1:length(data)
    figure
    
    turb_radius = data(i).turb;
    data(i).profile=nanmean(turb_radius);
    data(i).std=nanstd(turb_radius);
    perct=percent(i).perct1;
    contourf(x_axis,z_axis,perct',25);
    shading flat
    %caxis([0 70])
    colorbar
    title(['radius ',num2str(i)])
    xlabel('TKE (m^2s^{-2})')
    ylabel('z (m)')
end
%figure
% errorbar(z_axis,data(1).profile,  data(i).std,'.-r')
% hold on
% errorbar(z_axis,data(2).profile, data(i).std,'.-b')
% errorbar(z_axis,data(3).profile,data(i).std, '.-g')
% errorbar(z_axis,data(4).profile, data(i).std,'.-k')
% figure(1)
% plot(data(1).profile,z_axis,c)
% xlabel('TKE (m^2s^{-2})')
% ylabel('z (m)')
% legend('Rita 21','Rita 22','Katrina 28','Isabel 12', 'Isabel 14','Frances 01','Fabian 03')
% hold on
% figure(2)
% plot(data(2).profile,z_axis, c)
% xlabel('TKE (m^2s^{-2})')
% ylabel('z (m)')
% legend('Rita 21','Rita 22','Katrina 28','Isabel 12', 'Isabel 14','Frances 01','Fabian 03')
% hold on
% figure(3)
% plot(data(3).profile,z_axis, c)
% xlabel('TKE (m^2s^{-2})')
% ylabel('z (m)')
% legend('Rita 21','Rita 22','Katrina 28','Isabel 12', 'Isabel 14','Frances 01','Fabian 03')
% hold on
% figure(4)
% plot(data(4).profile,z_axis, c)
% xlabel('TKE (m^2s^{-2})')
% ylabel('z (m)')
% legend('Rita 21','Rita 22','Katrina 28','Isabel 12', 'Isabel 14','Frances 01','Fabian 03')
% hold on
% figure(5)
% plot(data(5).profile,z_axis, c)
% xlabel('TKE (m^2s^{-2})')
% ylabel('z (m)')
% legend('Rita 21','Katrina 28','Fabian 03')
% hold on