for k=1:20
input_data.file_number=k;
[tke_profile,x,z,vertical_tke]=tke_smartr(input_data);
reslt(1:length(vertical_tke),k)=vertical_tke';
end