course-idma2016/matlab/idma_float_plot.m

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

p = genpath('MITprof/'); addpath(p);%TBE

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