%load rocky
!rm rose.met
[m,n]=size(wlp);
nind=length(ind);
for i=n:n,
%for i=1:1,
%  for k=1:nind-1,
    for k=1:1,
   i1=ind(k);
%  i2=ind(k+1);
   i2=ind(nind);
   ii=i1:i2;
   str1=date_lbl(gregorian(jdlp(i1)),gregorian(jdlp(i2)));
   axis([-15 15 -15 15]);
  z=denan(wlp(ii,i));
  if(length(z)>60),
   nl=int2str(length(z));
   rose2(z,[0 100 1000]);
   title([cdf_groups(i,:) str1]);
   ylabel('percent')
   xlabel([nl ' samples']);
%    title(['bb wind stress x10 ' str1]);
  if (exist('rose.met')~=2)
   meta rose
  else
   meta
  end
  pause
  end
end
end
