cd ('C:\Documents and Settings\sylvie lorsolo\My Documents\research\data')
% [sweepfile,pathname] = uigetfile('*.','Select the sweepfile');
% data=wk1read(sweepfile);
slant_range=wk1read('slant_range.wk1');
az=wk1read('az.wk1');
az=az';
radar_long=-77.03749847412109*ones(length(slant_range),1);
radar_lat=35.07916641235352*ones(length(slant_range),1);


 % Be careful this work only for elevation angles smaller than 10
% el=sweepfile(30:32);
% el=str2num(el);
el=0.5;
% compute horizontal range
horizontal_range=(slant_range*cos(deg2rad(el)))/1000;
horizontal_range=horizontal_range*ones(1,length(az));
    ray=ones(length(slant_range),1)*az;
    for i=1:length(az)
    [lat,lon] = reckon(radar_lat,radar_long,km2deg(horizontal_range(:,i)),ray(:,i));
      lati(:,i)=lat;
      longi(:,i)=lon;
 end
    
    


