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];

colormap(jet);
axes('position',P1);
[cc h]=contourf(idist,RZ,gT,[7:.5:12.5 12.63]);  hold on;
set(h(:),'edgecolor','w');
caxis([7.5 12.7])
plot([7 7],[67 105],'k--');
axis([min(Rdist) max(Rdist) min(RZ) 250]);
t=text(1,240,['Temperature (' 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);
[cc,h]=contourf(idist,RZ,gS,[33.32:.02:33.38 33.4:.2:35 33.]); hold on
caxis([33.32 34]); 
set(h(:),'edgecolor','w');
plot([7 7],[67 105],'k--');
axis([min(Rdist) max(Rdist) min(RZ) 250]);
t=text(1,240,['Salinity (psu)']); 
set(t,'fontweight','bold','fontsize',12); 
set(gca,'fontsize',12,'tickdir','out');
set(gca,'ydir','reverse'); % ylabel('Depth (m)'); 
set(gca,'XTicklabel',[],'YaxisLocation','right');

return

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)');

