load bbvel
% get ustarw, wind only during times when 1 mab currents are available
ind=find(wlp~=NaN);
ustarwx=ustarw(ind);
windx=bb_str(ind);
wlpx=wlp(ind);
%
ind=find(ustarwx<2.5);
ind2=find(ustarwx>2.5);
%
axis([-15 15 -15 15]);
rose2(windx(ind));
if exist('wind1.met')==2,
!rm wind1.met
end
meta wind1
axis([-15 15 -15 15]);
rose2(windx(ind2));
if exist('wind2.met')==2,
!rm wind2.met
end
meta wind2
%
axis([-15 15 -15 15]);
rose2(wlpx(ind));
if exist('w1.met')==2,
!rm w1.met
end
meta w1
axis([-15 15 -15 15]);
rose2(wlpx(ind2));
if exist('w2.met')==2,
!rm w2.met
end
meta w2
%
compass2([mean(windx(ind)) mean(windx(ind2))])
if exist('windmean.met')==2
!rm windmean.met
end
meta windmean
compass2([mean(wlpx(ind)) mean(wlpx(ind2))])
if exist('wmean.met')==2
!rm wmean.met
end
meta wmean

