function [p]=idma_float_plot(nameFloat); % % IDMA_FLOAT_PLOT extracts the time series of one float % (e.g. '4900828') from argo_feb2013_2008_to_2010_model.nc % and return the result in MITprof format (p) % % Example: [p]=idma_float_plot('4900828'); % % Other float examples: '4900829', '4900830', '4900831' 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); %compute display time and depth axis: tim=1992+(p.prof_date-datenum([1992 1 1]))/365; dep=-p.prof_depth; %plot location of profiles: figureL; subplot(2,2,1); set(gca,'FontSize',12); plot(p.prof_lon(tt),p.prof_lat(tt),'.-'); grid on; title(['float ' nameFloat ' path']); %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(:,kk); o(p.prof_Tweight(:,kk)==0)=NaN; m=p.prof_Testim(:,kk); plot(tim(tt),o(tt),'b.-'); hold on; plot(tim(tt),m(tt),'r.-'); grid on; title(['T at ' d 'm']); legend('Argo','ECCOv4'); subplot(2,2,4); set(gca,'FontSize',12); kk=2; d=num2str(round(p.prof_depth(kk))); o=p.prof_S(:,kk); o(p.prof_Sweight(:,kk)==0)=NaN; m=p.prof_Sestim(:,kk); plot(tim(tt),o(tt),'b.-'); hold on; plot(tim(tt),m(tt),'r.-'); grid on; title(['S at ' d 'm']); legend('Argo','ECCOv4'); %plot S-T diagram: subplot(2,2,3); set(gca,'FontSize',12); S=p.prof_S; S(p.prof_Sweight(:,kk)==0)=NaN; T=p.prof_T; T(p.prof_Tweight(:,kk)==0)=NaN; plot(S(:),T(:),'b.'); hold on; plot(p.prof_Sestim(:),p.prof_Testim(:),'r.'); hold on; xlabel('salinity'); ylabel('temperature'); grid on; title('T-S profiles'); legend('Argo','ECCOv4'); figureL; %display temperature subplot(2,1,1); set(gca,'FontSize',12); o=p.prof_T; o(p.prof_Tweight==0)=NaN; [y,x]=meshgrid(dep,tim); pcolor(x(tt,:),y(tt,:),o(tt,:)); colorbar; shading flat; title(['Argo T']); %display salinity subplot(2,1,2); set(gca,'FontSize',12); o=p.prof_S; o(p.prof_Sweight==0)=NaN; [y,x]=meshgrid(dep,tim); pcolor(x(tt,:),y(tt,:),o(tt,:)); colorbar; shading flat; title(['Argo S']);