plotsta

 PLOTSTA Plots time series output locations on map
 and numbers the stations 
function [x,y,d,xi,yi]=plotsta(ecomcdf,tscdf,iplot);

0001 function [x,y,d,xi,yi]=plotsta(ecomcdf,tscdf,iplot);
0002 % PLOTSTA Plots time series output locations on map
0003 % and numbers the stations
0004 %function [x,y,d,xi,yi]=plotsta(ecomcdf,tscdf,iplot);
0005 [d,x,y]=kslice(ecomcdf,'depth');
0006 loc=mcvgt(tscdf,'loc');
0007 loc=loc.';  %columnize
0008 [m,n]=size(loc);
0009 for i=1:m;
0010   xi(i)=x(loc(i,1),loc(i,2));
0011   yi(i)=y(loc(i,1),loc(i,2));
0012 end
0013 if nargin==3,
0014   pslice(x,y,d);
0015   line(xi,yi,'linestyle','none','color','black',...
0016    'marker','.','markersize',20);
0017   line(xi,yi,'linestyle','none','color','white',...
0018    'marker','.','markersize',14);
0019   for i=1:m;
0020     text(xi(i),yi(i),int2str(i),'color','white');
0021   end 
0022 end