clear all; load avgUD_GF 
figure(2); clf; paper;
P1 = [.17 .52 .33 .37];
P2 = [.52 .52 .33 .37];
P3 = [.17 .14 .33 .37];
P4 = [.52 .14 .33 .37];

axes('position',P1);
% pcolor(idist,RZ,gT); shading interp; colorbar;  hold on
[cc h]=contour(idist,RZ,gT,[7:1:13],'k'); clabel(cc,h,'fontsize',12); hold on
contour(idist,RZ,gT,[7.5:1:13],'k'); 
[cc h]=contour(idist,RZ,gT,[12.63 12.63],'k--');
plot([7 7],[67 105],'k:');
% clabel(cc,h,'fontsize',12);
% axis([min(Rdist) max(Rdist) min(RZ) max(RZ)]);
axis([min(Rdist) max(Rdist) min(RZ) 250]);
p=plot(Rdist,repmat(min(RZ),length(Rdist),1),'ko');
set(p,'markerfacecolor','k')
t=text(1,240,['Temperature (' setstr(176) 'C)']); 
set(t,'fontweight','bold','fontsize',12); % caxis([8.9 12.7]); 
t=text(1,30,['>12.63' setstr(176) 'C']); 
set(t,'fontweight','bold','fontsize',12); % caxis([8.9 12.7]); 
set(gca,'fontsize',12,'tickdir','out');
set(gca,'ydir','reverse'); ylabel('Depth (m)'); 
set(gca,'XTicklabel',[]);

axes('position',P2);
% pcolor(idist,RZ,gS); shading interp; colorbar; 
[cc,h]=contour(idist,RZ,gS,[33:.2:35],'k'); 
clabel(cc,h,'fontsize',12); hold on;
contour(idist,RZ,gS,[33.1:.2:35],'k'); 
[cc,h]=contour(idist,RZ,gS,[33.32:.02:33.38],'k--'); 
plot([7 7],[67 105],'k:');
% clabel(cc,h,'fontsize',12);
% axis([min(Rdist) max(Rdist) min(RZ) max(RZ)]);
axis([min(Rdist) max(Rdist) min(RZ) 250]);
p=plot(Rdist,repmat(min(RZ),length(Rdist),1),'ko');
set(p,'markerfacecolor','k')
t=text(1,240,['Salinity (psu)']); 
set(t,'fontweight','bold','fontsize',12); % caxis([8.9 12.7]); 
set(gca,'fontsize',12,'tickdir','out');
set(gca,'ydir','reverse'); % ylabel('Depth (m)'); 
set(gca,'XTicklabel',[],'YaxisLocation','right');

axes('position',P3);
% pcolor(idist,RZ,gO2); shading interp; colorbar; 
[cc,h]=contour(idist,RZ,gO2,[1:1:6],'k'); 
clabel(cc,h,'fontsize',12); hold on
contour(idist,RZ,gO2,[1.5:1:6],'k'); 
[cc,h]=contour(idist,RZ,gO2,[5.6:.1:6],'k--'); 
clabel(cc,h,'fontsize',12); hold on
plot([7 7],[67 105],'k:');
% axis([min(Rdist) max(Rdist) min(RZ) max(RZ)]);
axis([min(Rdist) max(Rdist) min(RZ) 250]);
p=plot(Rdist,repmat(min(RZ),length(Rdist),1),'ko');
set(p,'markerfacecolor','k')
t=text(1,240,['Oxygen (ml/l)']); 
set(t,'fontweight','bold','fontsize',12); % caxis([8.9 12.7]); 
set(gca,'fontsize',12,'tickdir','out');
set(gca,'ydir','reverse'); ylabel('Depth (m)'); 
xlabel('Distance (km)');

axes('position',P4);
% pcolor(idist,RZ,gpH); shading interp; colorbar 
[cc,h]=contour(idist,RZ,gpH,[7:.1:8.2],'k'); 
clabel(cc,h,'fontsize',12); hold on;
% [cc,h]=contour(idist,RZ,gpH,[7.05:.1:8.5],'k--'); 
plot([7 7],[67 105],'k:');
% axis([min(Rdist) max(Rdist) min(RZ) max(RZ)]);
axis([min(Rdist) max(Rdist) min(RZ) 250]);
p=plot(Rdist,repmat(min(RZ),length(Rdist),1),'ko');
set(p,'markerfacecolor','k')
t=text(1,240,['pH']); 
set(t,'fontweight','bold','fontsize',12); % caxis([8.9 12.7]); 
set(gca,'fontsize',12,'tickdir','out','YaxisLocation','right');
set(gca,'ydir','reverse'); % ylabel('Depth (m)'); 
xlabel('Distance (km)');

