%%------------------------
%% Plot 
%%------------------------
idx=find(S>29 & S<35); sax = [min(S(idx)) max(S(idx))];
fax = [min(F) max(F)]; 
XL = 'mission time';

figure(1); clf; giffig

subplot(411); plot(time,depth,'k'); ylabel('depth'); set(gca,'Ydir','reverse');
xl=get(gca,'xlim'); set(gca,'box','off','YaxisLocation','left','tickdir','out'); 
hold on
P = get(gca,'position'); axes('position',P);
plot(time,pitch,'c'); hold on; plot(xl,[0 0],'k--');
set(gca,'YaxisLocation','right','Ycolor','c','box','off','tickdir','out');
ylabel('Pitch'); set(gca,'color','none','Ylim',[min(pitch) max(pitch)]);
set(gca,'xlim',xl,'XTickLabel',[]);

subplot(412); 
plot(time,T,'r'); set(gca,'ylim',[min(T) max(T)]);
ylabel('temperature'); hold on;
plot(time,Biotemp,'m');
set(gca,'Ycolor','r','box','off','YaxisLocation','left','tickdir','out');
set(gca,'xlim',xl,'XTickLabel',[]);
P = get(gca,'position'); axes('position',P);
plot(time,S,'b'); set(gca,'YaxisLocation','right','Ycolor','b','box','off');
set(gca,'Ylim',sax); 
ylabel('salinity'); set(gca,'color','none','tickdir','out');
set(gca,'xlim',xl,'XTickLabel',[]);

subplot(413); 
plot(time,Biolum,'m'); set(gca,'ylim',[min(biolum) max(biolum)]);
ylabel('biolum');
set(gca,'Ycolor','r','box','off','YaxisLocation','left','tickdir','out');
set(gca,'xlim',xl,'XTickLabel',[]);
P = get(gca,'position'); axes('position',P);
plot(time,Bioflow,'k'); set(gca,'YaxisLocation','right','box','off');
ylabel('biolum flow'); set(gca,'color','none','tickdir','out');
set(gca,'xlim',xl,'XTickLabel',[]);

subplot(414), plot(time,F,'g'); 
set(gca,'YaxisLocation','left','Ycolor','g','box','off','tickdir','out');
set(gca,'Ylim',fax); ylabel('fluorometer (volts)'); 
set(gca,'xlim',xl); hold on;
P = get(gca,'position'); axes('position',P)
plot(time,OBS,'r'); ylabel('OBS');
xlabel(XL); 
set(gca,'box','off','xlim',xl,'color','none','Yaxislocation','right','Ycolor','r','tickdir','out');

