nf=length(lidar_data.files);
for jfiles=184:nf
    ipnt=1;
    load(char(lidar_data.files(jfiles)));
    scan=scangen_2(lidar,'rhi');
    numscans=length(scan);
    for j=1:numscans
        hfig=figure('position',[100 100 900 700]);
pd=sind(lidar.azimuth).*lidar.xd+cosd(lidar.azimuth).*lidar.yd;
hf=pcolor(pd(scan(j).ind,:),lidar.zd(scan(j).ind,:),lidar.RWS(scan(j).ind,:))
caxis([-5 5])
set(hf,'edgecolor','none')
set(gca,'FontName','Arial')
set(gca,'FontSize',12)
xlabel('Cross Range (m)')
ylabel('Height AGL (m)')
title([datestr(scan(j).scantime),'Scan# ',num2str(j),' Lidar Azimuth: ',num2str(lidar.azimuth(1,1))])
aa=input('vortex present press y','s');
if aa=='y'
    lidar_data.scan(jfiles).ind(ipnt)=j;
    ipnt=ipnt+1;
end
close(hfig)
    end
end