i=1;
turb=[];
refle=[];
radi=[];
verti=[];
 hurricane='TKE\rita';
  jour='23';
 filename=dir(['C:\Documents and Settings\lorsolo\My Documents\research\HRD research\Gamache_soft\wfiles\',hurricane,'\',jour,'\*.w']) ;
 fichier=filename(i).name
 M=length(filename);

while i<(M-1)
    i
    fichier=filename(i).name
    [ws,radwind,w,dbz, x_fill, y_fill, z_fill, type_data]=read_w_files(hurricane,jour,i);
    
    j=1;
    
    if type_data == 0
        
        refle=[refle;dbz];
        verti=[verti;w];
        radi=[radi;radwind];
%         figure(i)
%         h=contourf(x_fill,z_fill,dbz',25);
%         velcolor=load(['C:\Documents and Settings\lorsolo\My Documents\research\programs\carbon43.dat']);
%         colormap(velcolor)
%         axis([-93000 93000 0 13000])
%          hold on
%         
%         figure(i+1)
%         h=contourf(x_fill,z_fill,radwind',25);
%         velcolor=load(['C:\Documents and Settings\lorsolo\My Documents\research\programs\carbon43.dat']);
%         colormap(velcolor)
%         caxis([-22 22])
%         colorbar
%         axis([-93000 93000 0 13000])
%         hold on
%         
%         
%         
%         figure(i+2)
%         h=contourf(x_fill,z_fill,w',25);
%         velcolor=load(['C:\Documents and Settings\lorsolo\My Documents\research\programs\carbon43.dat']);
%         colormap(velcolor)
%         caxis([-8 8])
%         colorbar
%         axis([-93000 93000 0 13000])
%         hold on
        
    else 
        
        turb=[turb;dbz.^2];
   
        figure(i+2)
        h=contourf(x_fill,z_fill,dbz.^2',25);
        velcolor=load(['C:\Documents and Settings\lorsolo\My Documents\research\programs\carbon43.dat']);
        colormap(velcolor)
        caxis([0 16])
        colorbar
        axis([-93000 93000 0 13000])
        hold on
        
        
    end
    
    [ws,radwind,w,dbz, x_fill, y_fill, z_fill]=read_w_files(hurricane,jour,i+2);
    
    
    
    if type_data == 0
%         refle=[refle;dbz];
%         radi=[radi;radwind];
%         verti=[verti;w];
%         figure(i)
%         contourf(-x_fill,z_fill,dbz',25);
%         caxis([-15 42])
%         colorbar
% %         saveas(gcf,['C:\Documents and Settings\lorsolo\My Documents\research\HRD research\Gamache_soft\wfiles\',hurricane,'\',jour,'\', fichier(1:13),'_dbz.fig'])
% %         
%         figure(i+1)
%         contourf(-x_fill,z_fill,radwind',25);
%         caxis([-22 22])
%         colorbar
% %         saveas(gcf,['C:\Documents and Settings\lorsolo\My Documents\research\HRD research\Gamache_soft\wfiles\',hurricane,'\',jour,'\',fichier(1:13),'_rad.fig'])
% %         
%         figure(i+2)
%         contourf(-x_fill,z_fill,w',25);
%         caxis([-8 8])
%         colorbar
%         saveas(gcf,['C:\Documents and Settings\lorsolo\My Documents\research\HRD research\Gamache_soft\wfiles\',hurricane,'\',jour,'\',fichier(1:13),'_vert.fig'])
else 
        turb=[turb;dbz.^2];
        figure(i+2)
        contourf(-x_fill,z_fill,dbz.^2',25);
        caxis([0 16])
        colorbar
        %saveas(gcf,['C:\Documents and Settings\lorsolo\My Documents\research\HRD research\Gamache_soft\wfiles\',hurricane,'\',jour,'\',fichier(1:13),'_tke.fig'])
    end
    
    
    
    if rem(i,2)==0
        i=i+3
    else
        i=i+1
    end
    
end