files=dir(['C:\Documents and Settings\sylvie lorsolo\Desktop\ess\frances\*dat']); 

el=[0.5;0.8;1.1;1.4;1.7;2.0;2.3;2.6;2.9;3.2;3.5;3.8;4.1;4.5;5.0;5.5;6.0;6.5;7.0;7.5;8.0;8.5];%[0.5;0.8;1.1;1.4;1.7;2.0;2.3;2.6;2.9;3.2;3.5;3.8;4.1;4.5;5.0;5.5;6.0;6.5;7.0;7.5;8.0;8.5];%[0.75 2.25 3.75 5.25 6.75 8.25];
eli=0.5:0.2:8.5;
cut=324;
rhi2=[];
rhii=[];
for fil=1:22;
        sweepfile=files(fil).name;
        mod_data=load(['C:\Documents and Settings\sylvie lorsolo\Desktop\ess\frances\',sweepfile]);
        mod_az=mod_data(length(mod_data),:);
        mod_az=mod_az';
        vel=mod_data(1:length(mod_data)-1,:);
        dif=abs(cut-mod_az);
        azimuth=find(dif==min(dif));
        mod_az(azimuth);
        if length(azimuth) >1
            break
        end
        rhi2=[rhi2 vel(:,azimuth)];
        
    end
     for rg=1:200
         temp=interp1(el,rhi2(rg,:),eli,'linear');
         rhii=[rhii;temp];
     end
rg=33.35:66.7:149*66.7+33.35;
 [th,r]=meshgrid(deg2rad(eli'),rg);
 [x,y]=pol2cart(th,r);

 pcolor(x,y,rhii(1:150,:))
 shading flat
 velcolor=load(['C:\Documents and Settings\sylvie lorsolo\My Documents\research\programs\carbon43.dat']);
  colormap(velcolor);

 v=[-4 4];
 caxis(v);
 colorbar