 %vers='v1'; 
  hurricane='frances';
  rep='meanless\rhi\';
files=dir(['C:\Documents and Settings\sylvie lorsolo\My Documents\research\',hurricane,'\data\modified\',rep,'*.dat']); 
    %[sweepfile,pathname] = uigetfile('*.','Select the sweepfile');
  for fil=73%1:length(files)
   sweepfile=files(fil).name;
  
 
data=load(['C:\Documents and Settings\sylvie lorsolo\My Documents\research\',hurricane,'\data\modified\',rep,sweepfile]);
el=data(length(data),:);
el=el';
eli=0.5:0.1:25;

rhi=data(1:length(data)-1,:);
rg=33.35:66.7:399*66.7+33.35;
rhii=[];
for rg=1:200
         temp=interp1(el,rhi(rg,:),eli,'linear');
        rhii=[rhii;temp];
    end
     rg=33.35:66.7:199*66.7+33.35;
  [th,r]=meshgrid(deg2rad(eli'),rg);
  [x,y]=pol2cart(th,r);
%  [th,r]=meshgrid(deg2rad(el'),rg);
%   [x,y]=pol2cart(th,r);
  figure
  h=pcolor(x,y,rhii(1:200,:));
  shading flat
  
   a=min(min(data(1:200,:)));
   b=max(max(data(1:200,:)));
%    v=[-8 5];
%    caxis(v)
 velcolor=load(['C:\Documents and Settings\sylvie lorsolo\My Documents\research\programs\carbon43.dat']);
 colormap(velcolor);
colorbar
caxis([-10 6]);
       %saveas(h,['C:\Documents and Settings\sylvie lorsolo\My Documents\research\frances\images\original\rhi\',sweepfile,'.jpg'])
  % saveas(h,['C:\Documents and Settings\sylvie lorsolo\Desktop\',sweepfile,'.jpg'])
    end