input_data.hurricane='\TKE\katrina';
input_data.jour='28\10km\';
input_data.path1=[];
[u10,cd, U10_full, Cd_full]=jobfile_cd_curve(5);
input_data.u10=u10;
input_data.cd=cd;
close all
input_data.hurricane='\TKE\frances';
input_data.jour='01\10km\';
for i=1:4:5
    figure
    input_data.file_number=i;
    [u_star,u_10_mean,Cd_interp,x_km]=u_star_calculation(input_data,'radar', 'hghg' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star,x_km,u_10_mean);
    
    %change program to use powell Cd data
    [u_star2,u_10_mean2,Cd_interp2,x_km2]=u_star_calculation(input_data,'radar','powell' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star2,x_km2,u_10_mean2);
    plot(x_km,u_10_interp,'.-c')
end
input_data.hurricane='\TKE\rita';
input_data.jour='20\10km\';
for i=1
    figure
    input_data.file_number=i;
    [u_star,u_10_mean,Cd_interp,x_km]=u_star_calculation(input_data,'radar', 'hghg' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star,x_km,u_10_mean);
    
    %change program to use powell Cd data
    [u_star2,u_10_mean2,Cd_interp2,x_km2]=u_star_calculation(input_data,'radar','powell' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star2,x_km2,u_10_mean2);
    plot(x_km,u_10_interp,'.-c')
end

input_data.hurricane='\TKE\rita';
input_data.jour='21\10km\';
for i=1:4:5
    figure
    input_data.file_number=i;
    [u_star,u_10_mean,Cd_interp,x_km]=u_star_calculation(input_data,'radar', 'hghg' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star,x_km,u_10_mean);
    
    %change program to use powell Cd data
    [u_star2,u_10_mean2,Cd_interp2,x_km2]=u_star_calculation(input_data,'radar','powell' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star2,x_km2,u_10_mean2);
    plot(x_km,u_10_interp,'.-c')
end

input_data.hurricane='\TKE\rita';
input_data.jour='22\10km\';
for i=1
    figure
    input_data.file_number=i;
    [u_star,u_10_mean,Cd_interp,x_km]=u_star_calculation(input_data,'radar', 'hghg' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star,x_km,u_10_mean);
    
    %change program to use powell Cd data
    [u_star2,u_10_mean2,Cd_interp2,x_km2]=u_star_calculation(input_data,'radar','powell' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star2,x_km2,u_10_mean2);
    plot(x_km,u_10_interp,'.-c')
end

input_data.hurricane='\TKE\rita';
input_data.jour='23\10km\';
for i=1:4:5
    figure
    input_data.file_number=i;
    [u_star,u_10_mean,Cd_interp,x_km]=u_star_calculation(input_data,'radar', 'hghg' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star,x_km,u_10_mean);
    
    %change program to use powell Cd data
    [u_star2,u_10_mean2,Cd_interp2,x_km2]=u_star_calculation(input_data,'radar','powell' );
    [u_sl_top, u_10_interp,x_km,z_axis,u_star] = U10_SL_calculation(input_data,u_star2,x_km2,u_10_mean2);
    plot(x_km,u_10_interp,'.-k')
end