clear all
clf

filename='RTE87_Barotropic_Current_1987.nc';

Experiment=ncread(filename,'Experiment')
disp(' ')
Notes=ncread(filename,'Note')

time1=ncread(filename,'yd_m1m2');
U1=ncread(filename,'U_m1m2');

time2=ncread(filename,'yd_m1m3');
U2=ncread(filename,'U_m1m3');

time3=ncread(filename,'yd_m2m3');
U3=ncread(filename,'U_m2m3');

figure(1)
clf
h1a=axes('units','inches','Position',[1.25 5.5 6.0 2.0]);
plot([136 264],[0 0],'k-', time1,1000*U1,'LineWidth',1)
axis([136 264 -29.99 29.99])
set(gca,'FontSize',12,'XTickLabel',[])
text(245,23,'\bf 2 \rightarrow 1','FontSize',14)

h2a=axes('units','inches','Position',[1.25 3.25 6.0 2.0]);
plot([136 264],[0 0],'k-', time2,1000*U2)
axis([136 264 -29.99 29.99])
set(gca,'FontSize',12,'XTickLabel',[]);
ylabel('Current (mm/s)','FontSize',14)
text(245,23,'\bf 3 \rightarrow 1','FontSize',14)

h3a=axes('units','inches','Position',[1.25 1.0 6.0 2.0]);
plot([136 264],[0 0],'k-', time3,1000*U3)
axis([136 264 -29.99 29.99])
set(gca,'FontSize',12);
text(245,23,'\bf 2 \rightarrow 3','FontSize',14)
xlabel('Yearday 1987','FontSize',14)