ii=input('file number: ');
fid=fopen(char(files(ii)));
lidar=read_CSM(fid,0.);
fclose(fid);
dfig=plot(lidar.azimuth(:,1));
azbear=input(' Azimuth: ');
ind=find(abs(lidar.azimuth(:,1)-36)<=.1);
lidar.elevation(ind,1);
ind2=ind(2:end-4);
lidar.elevation(ind2,1);
hfig=pcolor(sqrt(lidar.xd(ind2,:).^2+lidar.yd(ind2,:).^2),3.2808*lidar.zd(ind2,:),rws(ind2,:));
set(hfig,'EdgeColor','none')
set(gca,'XLim',[0 1000])
set(gca,'YLim',[0 1000])