files=uipickfiles;

nfiles=length(files)

for j=1:nfiles
    load(char(files(j)))
    vd = interp1(prof.heights,prof.v,lidar.zd);
wd = interp1(prof.heights,prof.w,lidar.zd);
ud = interp1(prof.heights,prof.u,lidar.zd);
Vh_mod=(sind(lidar.azimuth).*ud + cosd(lidar.azimuth).*vd).*cosd(lidar.elevation) + wd.*sind(lidar.elevation);
h2=figure('Position',[500 300 600 600]);
hd=pcolor(lidar.xd,lidar.yd,lidar.rws_filt)
axis equal
set(hd,'edgecolor','none')
caxis([-25 25])
xlim([-2000 2000])
ylim([-2000 2000])
F(j)=getframe(h2);
close(gcf)
end