% PCO2_UWAY_PLOT.M

start_time = datenum('08-Feb-2002 00:00');
end_time = datenum('09-Feb-2002 00:00');
pco2_num = [];
% 
% for i = [39]     %These are Julian days - Jan 1 = 1 (Duh!). Change this for the days you want to process.
%     in_file = sprintf('pco2_2002%03d_co2exe.txt', i);     %open the pCO2 input file
%     data = load(in_file);
%     pco2_num = [pco2_num; data];
%     clear data;
% end

x = find(~isnan(pco2(:,13)) & (pco2(:,1) >= start_time & pco2(:,1) <= end_time));
subplot(2,2,1)
scatter(pco2(x,4),pco2(x,3),35,pco2(x,13),'filled');
colorbar;
title('Surface pCO_{2} [\mumol mol^{-1}]: Northern patch, Survey 3');
xlabel('Longitude [deg]');
ylabel('Latitude [deg]');
subplot(2,2,2)
scatter(pco2(x,4),pco2(x,3),35,pco2(x,2),'filled');
colorbar;
title('Fluorescence [volts]: Northern patch, Survey 3');
xlabel('Longitude [deg]');
ylabel('Latitude [deg]');
subplot(2,2,3)
scatter(pco2(x,4),pco2(x,3),35,pco2(x,10),'filled');
colorbar;
title('Temperature [C^o]: Northern patch, Survey 3');
xlabel('Longitude [deg]');
ylabel('Latitude [deg]');
subplot(2,2,4)
scatter(pco2(x,4),pco2(x,3),35,pco2(x,5),'filled');
colorbar;
title('Salinity [PSU]: Northern patch, Survey 3');
xlabel('Longitude [deg]');
ylabel('Latitude [deg]');
