course-idma2015/computing/idma_float_plot.m

function [p]=idma_float_plot(nameFloat);
%
%four Argo floats in subtropical Pacific:
% nameFloat='2900828';
% nameFloat='2900829';
% nameFloat='2900830';
% nameFloat='2900831';

addpath gcmfaces;
addpath MITprof;
MITprof_global

%load Argo data set:
dirIn='release1/MITprof/';
nameFile='argo_feb2013_2008_to_2010_model.nc';
prof=MITprof_load([dirIn nameFile]);

%isolate one instrument time series:
p=MITprof_subset(prof,'descr',nameFloat);

%find increasing time order, for plotting:
[tmp1,tt]=sort(p.prof_date);

%plot location of profiles:
figureL;
subplot(2,2,1); set(gca,'FontSize',12); 
plot(p.prof_lon(tt),p.prof_lat(tt),'.-'); grid on; 
t0=num2str(min(p.prof_YYYYMMDD));
t1=num2str(max(p.prof_YYYYMMDD));
title(['float ' nameFloat ' from ' t0 ' to ' t1]);
% axis([-180 -120 10 40]); 

%plot T time series at select depths:
subplot(2,2,2); set(gca,'FontSize',12); 
kk=2; d=num2str(round(p.prof_depth(kk)));
o=p.prof_T(tt,kk); m=p.prof_Testim(tt,kk);
plot(o,'b.-'); hold on; plot(m,'r.-'); grid on; 
title(['T at ' d 'm, argo (blue) and ecco v4 (red)']);

subplot(2,2,4); set(gca,'FontSize',12); 
kk=10; d=num2str(round(p.prof_depth(kk)));
o=p.prof_T(tt,kk); m=p.prof_Testim(tt,kk);
plot(o,'b.-'); hold on; plot(m,'r.-'); grid on; 
title(['T at ' d 'm, argo (blue) and ecco v4 (red)']);

%plot S-T diagram:
subplot(2,2,3); set(gca,'FontSize',12); 
plot(p.prof_S(:),p.prof_T(:),'b.'); hold on;
plot(p.prof_Sestim(:),p.prof_Testim(:),'r.'); hold on;
xlabel('salinity'); ylabel('temperature');
grid on; title('argo (blue) and ecco v4 (red)');

1	gforget	1.1	function [p]=idma_float_plot(nameFloat);
2			%
3			%four Argo floats in subtropical Pacific:
4			% nameFloat='2900828';
5			% nameFloat='2900829';
6			% nameFloat='2900830';
7			% nameFloat='2900831';
8
9			addpath gcmfaces;
10			addpath MITprof;
11			MITprof_global
12
13			%load Argo data set:
14			dirIn='release1/MITprof/';
15			nameFile='argo_feb2013_2008_to_2010_model.nc';
16			prof=MITprof_load([dirIn nameFile]);
17
18			%isolate one instrument time series:
19			p=MITprof_subset(prof,'descr',nameFloat);
20
21			%find increasing time order, for plotting:
22			[tmp1,tt]=sort(p.prof_date);
23
24			%plot location of profiles:
25			figureL;
26	gforget	1.2	subplot(2,2,1); set(gca,'FontSize',12);
27	gforget	1.1	plot(p.prof_lon(tt),p.prof_lat(tt),'.-'); grid on;
28			t0=num2str(min(p.prof_YYYYMMDD));
29			t1=num2str(max(p.prof_YYYYMMDD));
30			title(['float ' nameFloat ' from ' t0 ' to ' t1]);
31			% axis([-180 -120 10 40]);
32
33			%plot T time series at select depths:
34	gforget	1.2	subplot(2,2,2); set(gca,'FontSize',12);
35			kk=2; d=num2str(round(p.prof_depth(kk)));
36	gforget	1.1	o=p.prof_T(tt,kk); m=p.prof_Testim(tt,kk);
37	gforget	1.2	plot(o,'b.-'); hold on; plot(m,'r.-'); grid on;
38	gforget	1.1	title(['T at ' d 'm, argo (blue) and ecco v4 (red)']);
39	gforget	1.2
40			subplot(2,2,4); set(gca,'FontSize',12);
41			kk=10; d=num2str(round(p.prof_depth(kk)));
42	gforget	1.1	o=p.prof_T(tt,kk); m=p.prof_Testim(tt,kk);
43	gforget	1.2	plot(o,'b.-'); hold on; plot(m,'r.-'); grid on;
44	gforget	1.1	title(['T at ' d 'm, argo (blue) and ecco v4 (red)']);
45	gforget	1.2
46			%plot S-T diagram:
47			subplot(2,2,3); set(gca,'FontSize',12);
48			plot(p.prof_S(:),p.prof_T(:),'b.'); hold on;
49			plot(p.prof_Sestim(:),p.prof_Testim(:),'r.'); hold on;
50			xlabel('salinity'); ylabel('temperature');
51			grid on; title('argo (blue) and ecco v4 (red)');
52