/[MITgcm]/MITgcm_contrib/gael/comm/course-idma2016/matlab/idma_float_plot.m
ViewVC logotype

Annotation of /MITgcm_contrib/gael/comm/course-idma2016/matlab/idma_float_plot.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.2 - (hide annotations) (download)
Fri Jan 8 04:35:49 2016 UTC (9 years, 6 months ago) by gforget
Branch: MAIN
Changes since 1.1: +37 -13 lines
- revise display
- add time-depth plots

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

  ViewVC Help
Powered by ViewVC 1.1.22