%
% Routine to plot tidal ellipse
%
hold off
load /home2/ocn/jmanning/gbr/posn_surf.dat
p=posn_surf;
load /home2/ocn/jmanning/futil/tides/file2_surf.dat
e=file2_surf;
%
georbat;
axis([-69.4 -68.2 40.4 41.2]);
set(gca,'xtick',[-69.4 -69.2 -69 -68.8 -68.6 -68.4]);
set(gca,'xticklabels',['     ';' 69.2';'     ';' 68.8';'     ';' 68.4']);
set(gca,'ytick',[40.4 40.6 40.8 41.0 41.2]);
set(gca,'yticklabels',[' 40.4';'     ';' 40.8';'     ';'41.2 ']);
hold on;
%
[nr,nc]=size(e);
for i=1:nr;
%
% Note: tide program gives angle anticlockwise from east
%
ang=-1*e(i,4);
%
ellipse(p(i,3),p(i,2),e(i,2),-1*(e(i,3)),ang,20,7);
end;
%
for i=[1:2:7 8];
  h=text(p(i,3),p(i,2),num2str(p(i,1)));
  set(h,'fontsize',18);
end;
%h=text(p(8,3),p(8,2),'6');
%set(h,'color','c')
%hx=text(p(3,3),p(3,2),'3');
%set(hx,'color','c')
set(gcf,'paperposition',[.5 2 6 5.26]);



