%% biological saturation
figure(4)
subplot(121);
subplot('Position',[0.1 0.1 0.42 0.85])
%plot(Ra(:,1), biosat,'r.');
%plot(X, biosat, 'r.');
scatter(X, Y, 200, biosat, 'filled');

axis equal
title('Biological Saturation from O_2/Ar', 'fontsize', 18, 'fontweight', 'bold');
%xlabel('time')
%xlabel('Latitude (degree)', 'fontsize', 14, 'fontweight', 'bold');
xlabel('Longitude (degrees)', 'fontsize', 14, 'fontweight', 'bold');
%ylabel('Bio Sat (%)', 'fontsize', 14, 'fontweight', 'bold');
ylabel('Latitude (degrees)', 'fontsize', 14, 'fontweight', 'bold');
%labelxwithdate(gca, 15); 
grid on;
box on;
cb = colorbar;
caxis([-5 10]);
colormap('jet')
ylabel(cb, 'BioSat %', 'fontsize', 14, 'fontweight', 'bold');
set(gca, 'fontsize', 14);
% add station number
[stn_name, stn_lat, stn_lon] = getstationlabels (stnfname);
% add station number
hold on;
%scatter(stn_lon,stn_lat,50,'^','filled','MarkerFaceColor','k');
scatter(stn_lon,stn_lat,50,'k');
%text(stn_lon+0.1,stn_lat,stn_name);

xlim([-71.1 -70.6]); %use this pair for standard maps
ylim([39.6 40.5]);

%xlim([-71.6 -70.2]); %use this pair for full size maps
%ylim([39.2 40.5]);


text(stn_lon(1:2:end)+0.04,stn_lat(1:2:end),stn_name(1:2:end,:), 'fontsize', 14);

%% temperature
subplot(122)
subplot('Position',[0.55 0.1 0.42 0.85])
%plot(Ra(:,1), T,'b.');
%plot(X, T, 'b.');
scatter(X, Y, 200, Tthermistor,'filled');

axis equal
title('Temperature', 'fontsize', 18, 'fontweight', 'bold');
%xlabel('time')
xlabel('Longitude (degree)', 'fontsize', 14, 'fontweight', 'bold');
ylabel('Longitude (degree)', 'fontsize', 14, 'fontweight', 'bold')
%labelxwithdate(gca, 15); 
grid on;
box on;
cb = colorbar;
caxis([15 25]);
colormap('jet')
ylabel(cb, 'Temp (deg C)', 'fontsize', 14, 'fontweight', 'bold');
set(gca, 'fontsize', 14);
% add station number
[stn_name, stn_lat, stn_lon] = getstationlabels (stnfname);
% add station number
hold on;
%scatter(stn_lon,stn_lat,50,'^','filled','MarkerFaceColor','k');
scatter(stn_lon,stn_lat,50,'k');
%text(stn_lon+0.1,stn_lat,stn_name);

xlim([-71.1 -70.6]); %use this pair for standard maps
ylim([39.6 40.5]);

%xlim([-71.6 -70.2]); %use this pair for full size maps
%ylim([39.2 40.5]);

text(stn_lon(1:2:end)+0.04,stn_lat(1:2:end),stn_name(1:2:end,:), 'fontsize', 14);


