%Compute mask
y(1) = min(bZ);
y(2) = max(bZ);

y(3) = y(2);
y(4) = y(1);

x1(1) = datenum('15-Apr-2005 14:59:23')
x1(3) = datenum('20-Apr-2005 20:59:23')

x1(2) = x(1);
x2(4) = x(3);

  
f = figure('Visible','on');
hold on;
[C,h] = contourf(time,bZ,BTemp,20);
set(h,'LineStyle','none')
foo = axis;
foo(3) = 0;
axis(foo);
set(gca,'xtick',[time(1):7:time(end)])
datetick('x','keeplimits','keepticks');
pl

set(gca,'YDir','reverse');
colorbar;
ylabel('depth (m)');
xlabel('GMT');
title('Temperature (deg C)');

set(f,'Units','pixels');
rect = get(f,'Position');
rect(3) = 550;
rect(4) = 200;
set(f,'Position',rect);
set(f,'PaperPositionMode','auto');

print -dpng Temp
print -depsc2 Temp

f = figure('Visible','on');
[C,h] = contourf(time,bZ,BSal,20);
set(h,'LineStyle','none')
foo = axis;
foo(3) = 0;
axis(foo);
set(gca,'xtick',[time(1):7:time(end)])
datetick('x','keeplimits','keepticks');
set(gca,'YDir','reverse') ;
colorbar;
ylabel('depth (m)');
xlabel('GMT');
title('Salinity (PSU)');

set(f,'Units','pixels');
rect = get(f,'Position');
rect(3) = 550;
rect(4) = 200;
set(f,'Position',rect);
set(f,'PaperPositionMode','auto');

print -dpng Sal
print -depsc2 Sal

f = figure('Visible','on');
[C,h] = contourf(time,bZ,BFluor,20);
set(h,'LineStyle','none')
foo = axis;foo(3) = 0;axis(foo);
caxis([0 5]);
set(gca,'xtick',[time(1):7:time(end)])
datetick('x','keeplimits','keepticks');
set(gca,'YDir','reverse') ;
colorbar;

ylabel('depth (m)');
xlabel('GMT');
title('Fluorometer (0-5V)');

set(f,'Units','pixels');
rect = get(f,'Position');
rect(3) = 550;
rect(4) = 200;
set(f,'Position',rect);
set(f,'PaperPositionMode','auto');

print -dpng Fluor
print -depsc2 Fluor


f = figure('Visible','on');
[C,h] = contourf(time,bZ,BOxygen,20);
set(h,'LineStyle','none')
foo = axis; foo(3) = 0; axis(foo);
caxis([0 5]);
set(gca,'xtick',[time(1):7:time(end)])
datetick('x','keeplimits','keepticks');
set(gca,'YDir','reverse') ;
colorbar;

ylabel('depth (m)');
xlabel('GMT');
title('Oxygen (0-5V)');

set(f,'Units','pixels');
rect = get(f,'Position');
rect(3) = 550;
rect(4) = 200;
set(f,'Position',rect);
set(f,'PaperPositionMode','auto');

print -dpng Oxygen
print -depsc2 Oxygen

f = figure('Visible','on');
[C,h] = contourf(time,bZ,BTrans,20);
set(h,'LineStyle','none')
foo = axis; foo(3) = 0; axis(foo);
set(gca,'xtick',[time(1):7:time(end)])
datetick('x','keeplimits','keepticks');
set(gca,'YDir','reverse') ;
colorbar;

ylabel('depth (m)');
xlabel('GMT');
title('Transmissometer (0-5V)');

set(f,'Units','pixels');
rect = get(f,'Position');
rect(3) = 550;
rect(4) = 200;
set(f,'Position',rect);
set(f,'PaperPositionMode','auto');

print -dpng Transmissometer
print -depsc2 Transmissometer


!cp Temp.png W:/MOOSVP/.
!cp Sal.png W:/MOOSVP/.
!cp Fluor.png W:/MOOSVP/.
!cp Oxygen.png W:/MOOSVP/.
!cp Transmissometer.png W:/MOOSVP/.

disp(['Done at ' datestr(clock)]);