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 |
gforget |
1.2 |
%compute display time and depth axis: |
23 |
|
|
tim=1992+(p.prof_date-datenum([1992 1 1]))/365; |
24 |
|
|
dep=-p.prof_depth; |
25 |
|
|
|
26 |
gforget |
1.1 |
%plot location of profiles: |
27 |
|
|
figureL; |
28 |
|
|
subplot(2,2,1); set(gca,'FontSize',12); |
29 |
|
|
plot(p.prof_lon(tt),p.prof_lat(tt),'.-'); grid on; |
30 |
gforget |
1.2 |
title(['float ' nameFloat ' path']); |
31 |
gforget |
1.1 |
|
32 |
|
|
%plot T time series at select depths: |
33 |
|
|
subplot(2,2,2); set(gca,'FontSize',12); |
34 |
|
|
kk=2; d=num2str(round(p.prof_depth(kk))); |
35 |
gforget |
1.2 |
o=p.prof_T(:,kk); |
36 |
|
|
o(p.prof_Tweight(:,kk)==0)=NaN; |
37 |
|
|
m=p.prof_Testim(:,kk); |
38 |
|
|
plot(tim(tt),o(tt),'b.-'); hold on; |
39 |
|
|
plot(tim(tt),m(tt),'r.-'); grid on; |
40 |
|
|
title(['T at ' d 'm']); legend('Argo','ECCOv4'); |
41 |
gforget |
1.1 |
|
42 |
|
|
subplot(2,2,4); set(gca,'FontSize',12); |
43 |
gforget |
1.2 |
kk=2; d=num2str(round(p.prof_depth(kk))); |
44 |
|
|
o=p.prof_S(:,kk); |
45 |
|
|
o(p.prof_Sweight(:,kk)==0)=NaN; |
46 |
|
|
m=p.prof_Sestim(:,kk); |
47 |
|
|
plot(tim(tt),o(tt),'b.-'); hold on; |
48 |
|
|
plot(tim(tt),m(tt),'r.-'); grid on; |
49 |
|
|
title(['S at ' d 'm']); legend('Argo','ECCOv4'); |
50 |
gforget |
1.1 |
|
51 |
|
|
%plot S-T diagram: |
52 |
|
|
subplot(2,2,3); set(gca,'FontSize',12); |
53 |
gforget |
1.2 |
S=p.prof_S; S(p.prof_Sweight(:,kk)==0)=NaN; |
54 |
|
|
T=p.prof_T; T(p.prof_Tweight(:,kk)==0)=NaN; |
55 |
|
|
plot(S(:),T(:),'b.'); hold on; |
56 |
gforget |
1.1 |
plot(p.prof_Sestim(:),p.prof_Testim(:),'r.'); hold on; |
57 |
|
|
xlabel('salinity'); ylabel('temperature'); |
58 |
gforget |
1.2 |
grid on; title('T-S profiles'); legend('Argo','ECCOv4'); |
59 |
gforget |
1.1 |
|
60 |
gforget |
1.2 |
figureL; |
61 |
|
|
%display temperature |
62 |
|
|
subplot(2,1,1); set(gca,'FontSize',12); |
63 |
|
|
o=p.prof_T; |
64 |
|
|
o(p.prof_Tweight==0)=NaN; |
65 |
|
|
[y,x]=meshgrid(dep,tim); |
66 |
|
|
pcolor(x(tt,:),y(tt,:),o(tt,:)); |
67 |
|
|
colorbar; shading flat; title(['Argo T']); |
68 |
|
|
%display salinity |
69 |
|
|
subplot(2,1,2); set(gca,'FontSize',12); |
70 |
|
|
o=p.prof_S; |
71 |
|
|
o(p.prof_Sweight==0)=NaN; |
72 |
|
|
[y,x]=meshgrid(dep,tim); |
73 |
|
|
pcolor(x(tt,:),y(tt,:),o(tt,:)); |
74 |
|
|
colorbar; shading flat; title(['Argo S']); |
75 |
gforget |
1.1 |
|