load rocky
load bath40
!rm rose.met
[m,n]=size(wlp);
nind=length(ind);
is=[1 3 6 12 14 16 19 21 24];   %surface 0 -10
im=[2 4 7 13 15 17 20 22 25];   % mid 10-30
inb=[5 8 9 23];                 % near bot 
ib=[10 11 ];
%
%  get coastline and 40 m isobath
load coast
load bath40
%  set up plot
lat_mid=42.3;
lat_max=42.6;
lon_mid=-70.4;
dlon=1.4;
lat_fac=cos(lat_mid*pi/180);
dlat=dlon*lat_fac;
% factor for roses
%
yfac=1./111;   % km => deg
xfac=yfac/lat_fac;

for k=4:6,    % fall 90- spring 91
hold off
axis('square')
axis([lon_mid-.5*dlon lon_mid+.5*dlon lat_mid-.5*dlat lat_mid+.5*dlat]);
plot(coast(:,1),coast(:,2));
hold on
mplot(bath40,'.')
hold on
for i=im,
   i1=ind(k);
   i2=ind(k+1);
   ii=i1:i2;
   str1=date_lbl(gregorian(jdlp(i1)),gregorian(jdlp(i2)));
   z=denan(wlp(ii,i));
  if(length(z)>60),
   [z]=shot(z,1,.5);
   cdf_groups(i,:)
   nl=int2str(length(z));
   [roze]=rose3(z,[0 100 1000]);
   xroze=real(roze)*xfac+lonlat(i,1);
   yroze=imag(roze)*yfac+lonlat(i,2); 
   plot(xroze,yroze,'g-')
   hold on
  end
end
title([str1 ' 10-30 m']);
if (exist('rose.met')~=2)
  meta rose
else
  meta
end
hold off
pause
end
