% South Sound Oxygen(ssoxy5.m)
% Survey 5 (9909B664)

clf
% orient landscape;
orient tall;
load d:\ssound\ss.mat % Permanent data for all cruises
load d:\ssound\sstrans5.mat % Temporary data for this cruise

subplot(211);
adspuget;
hold on
s=size(oxy);
% Carr Inlet
surf(x(1:s(1),67:75),y(1:s(1),67:75),dep(1:s(1),67:75)*-1,oxy(1:s(1),67:75)),shading('interp');
% Main & Case Inlet
surf(x(1:s(1),[45:59 61:67 76]),y(1:s(1),[45:59 61:67 76]),dep(1:s(1),[45:59 61:67 76])*-1,oxy(1:s(1),[45:59 61:67 76])),shading('interp');
% Budd Inlet
surf(x(1:s(1),[4:9 11 13 14 3 56]),y(1:s(1),[4:9 11 13 14 3 56]),dep(1:s(1),[4:9 11 13 14 3 56])*-1,oxy(1:s(1),[4:9 11 13 14 3 56])),shading('interp');
% Eld Inlet
surf(x(1:s(1),[14:17]),y(1:s(1),[14:17]),dep(1:s(1),[14:17])*-1,oxy(1:s(1),[14:17])),shading('interp');
% Totten Inlet
surf(x(1:s(1),[24 25 23 22 21 20 27 28]),y(1:s(1),[24 25 23 22 21 20 27 28]),dep(1:s(1),[24 25 23 22 21 20 27 28])*-1,oxy(1:s(1),[24 25 23 22 21 20 27 28])),shading('interp');
% Oakland Bay
surf(x(1:s(1),[28:36]),y(1:s(1),[28:36]),dep(1:s(1),[28:36])*-1,oxy(1:s(1),[28:36])),shading('interp');
% Pickering Passage
surf(x(1:s(1),[28 37 40:44 49]),y(1:s(1),[28 37 40:44 49]),dep(1:s(1),[28 37 40:44 49])*-1,oxy(1:s(1),[28 37 40:44 49])),shading('interp');

caxis([0 20]), colorbar;
%colorbar;

view(10,55)
box on
axis([-123.10 -122.55 47.00 47.45 -175 0]);
set(gca,'xlim',[-123.10 -122.55]);
set(gca,'xtick',[-123.1:0.1:-122.5]);
set(gca,'xticklabel',[123.1 123.0 122.9 122.8 122.7 122.6 122.5]);
set(gca,'zlim',[-175 0]);
set(gca,'ztick',[-175:25:0]);
set(gca,'zticklabel',[175:-25:0]);
zlabel('Depth (m)');
ylabel('North Latitude (deg)');
WL=xlabel('West Longitude (deg)');
T1=title({'South Sound Oxygen','Sept. 20-23, 1999'});
set(T1,'FontName','Times','FontWeight','Bold','FontSize',18);

subplot(212);
adspuget;
hold on
s=size(oxy);
% Carr Inlet
surf(x(1:s(1),67:75),y(1:s(1),67:75),dep(1:s(1),67:75)*-1,oxy(1:s(1),67:75)),shading('interp');
% Main & Case Inlet
surf(x(1:s(1),[45:59 61:67 76]),y(1:s(1),[45:59 61:67 76]),dep(1:s(1),[45:59 61:67 76])*-1,oxy(1:s(1),[45:59 61:67 76])),shading('interp');
% Budd Inlet
surf(x(1:s(1),[4:9 11 13 14 3 56]),y(1:s(1),[4:9 11 13 14 3 56]),dep(1:s(1),[4:9 11 13 14 3 56])*-1,oxy(1:s(1),[4:9 11 13 14 3 56])),shading('interp');
% Eld Inlet
surf(x(1:s(1),[14:17]),y(1:s(1),[14:17]),dep(1:s(1),[14:17])*-1,oxy(1:s(1),[14:17])),shading('interp');
% Totten Inlet
surf(x(1:s(1),[24 25 23 22 21 20 27 28]),y(1:s(1),[24 25 23 22 21 20 27 28]),dep(1:s(1),[24 25 23 22 21 20 27 28])*-1,oxy(1:s(1),[24 25 23 22 21 20 27 28])),shading('interp');
% Oakland Bay
surf(x(1:s(1),[28:36]),y(1:s(1),[28:36]),dep(1:s(1),[28:36])*-1,oxy(1:s(1),[28:36])),shading('interp');
% Pickering Passage
surf(x(1:s(1),[28 37 40:44 49]),y(1:s(1),[28 37 40:44 49]),dep(1:s(1),[28 37 40:44 49])*-1,oxy(1:s(1),[28 37 40:44 49])),shading('interp');

%caxis([28 30]), colorbar;
colorbar;

view(-65,55)
box on
axis([-123.10 -122.55 47.00 47.45 -175 0]);
set(gca,'xlim',[-123.10 -122.55]);
set(gca,'xtick',[-123.1:0.1:-122.5]);
set(gca,'xticklabel',[123.1 123.0 122.9 122.8 122.7 122.6 122.5]);
set(gca,'zlim',[-175 0]);
set(gca,'ztick',[-175:25:0]);
set(gca,'zticklabel',[175:-25:0]);
zlabel('Depth (m)');
ylabel('North Latitude (deg)');
WL=xlabel('West Longitude (deg)');

print -djpeg90 d:\ssound\ssoxy5.jpg