course-idma2016/matlab/idma_box_subset.m

function [prof_T,prof_Testim,prof_Tweight,...
          prof_S,prof_Sestim,prof_Sweight]=idma_box_subset(box_l,box_L,box_D);
%
% IDMA_BOX_SUBSET extracts profiles located in a region defined
%    by ranges of longitude (box_l), latitude (box_L), and 
%    depth (box_D) accross data collections. Reurns 
%    observed profiles (prof_T and prof_S) and ECCO v4
%    profiles (prof_Testim and prof_Sestim) in MITprof format;
%
% Example: box_l=[-180 -140]; box_L=[-10 10]-30; box_D=[0 300];
%          [prof_T,prof_Testim,prof_S,prof_Sestim]=idma_box_subset(box_l,box_L,box_D);

demoIO=1;%set this to one to activate additional demonstration of IO

%%add MITprof routines to matlab path:
dir_MITprof=[pwd filesep 'MITprof' filesep];
p = genpath(dir_MITprof); addpath(p);%TBE

%%set box location
fprintf('\n box definition: \n');
fprintf(' longitude range: %d %d \n',box_l);
fprintf(' latitude  range: %d %d \n',box_L);
fprintf(' depth     range: %d %d \n\n',box_D);

%%demo: load a netcdf file
if demoIO;
dir_nc='release1/MITprof_release1/';
ncload([dir_nc 'argo_feb2013_1992_to_2007_model.nc']);%TBE
clear prof_*;
end;

%%demo: MITprof_load and MITprof_subset
if demoIO;
dir_nc='release1/MITprof_release1/';
tmp1=MITprof_load([dir_nc 'argo_feb2013_1992_to_2007_model.nc']);%TBE

tmp2=MITprof_subset(tmp1,'lon',box_l,'lat',box_L,'depth',box_D);%TBE
figure; set(gca,'FontSize',12);
plot(tmp1.prof_lon,tmp1.prof_lat,'.');
hold on; plot(tmp2.prof_lon,tmp2.prof_lat,'r.');
axis([-180 180 -90 90]); title('selected profile locations');
xlabel('longitude'); ylabel('latitude') 
clear tmp*;
end;

%%load complete data set for select variables
dir_nc='release1/MITprof_release1/';
list_nc={'argo_feb2013_1992_to_2007_model.nc','argo_feb2013_2008_to_2010_model.nc','argo_feb2013_2011_to_2012_model.nc',...
         'climode_feb2013_model.nc','itp_feb2013_model.nc','seals_feb2013_model.nc'};
%note: ctd_feb2013_model.nc and xbt_feb2013_model.nc are left out here for simplicity
%note: it takes a couple minutes to load all of the data

list_var={'prof_T','prof_Testim','prof_Tweight',...
'prof_S','prof_Sestim','prof_Sweight'};%TBE

for vv=1:length(list_var);
  tmp1=MITprof_stats_load(dir_nc,list_nc,'',list_var{vv});%TBE
  tmp1=MITprof_subset(tmp1,'lon',box_l,'lat',box_L,'depth',box_D);%TBE
  eval([list_var{vv} '=tmp1;']);%TBE
end;

clear tmp1 vv;

1	function [prof_T,prof_Testim,prof_Tweight,...
2	prof_S,prof_Sestim,prof_Sweight]=idma_box_subset(box_l,box_L,box_D);
3	%
4	% IDMA_BOX_SUBSET extracts profiles located in a region defined
5	% by ranges of longitude (box_l), latitude (box_L), and
6	% depth (box_D) accross data collections. Reurns
7	% observed profiles (prof_T and prof_S) and ECCO v4
8	% profiles (prof_Testim and prof_Sestim) in MITprof format;
9	%
10	% Example: box_l=[-180 -140]; box_L=[-10 10]-30; box_D=[0 300];
11	% [prof_T,prof_Testim,prof_S,prof_Sestim]=idma_box_subset(box_l,box_L,box_D);
12
13	demoIO=1;%set this to one to activate additional demonstration of IO
14
15	%%add MITprof routines to matlab path:
16	dir_MITprof=[pwd filesep 'MITprof' filesep];
17	p = genpath(dir_MITprof); addpath(p);%TBE
18
19	%%set box location
20	fprintf('\n box definition: \n');
21	fprintf(' longitude range: %d %d \n',box_l);
22	fprintf(' latitude range: %d %d \n',box_L);
23	fprintf(' depth range: %d %d \n\n',box_D);
24
25	%%demo: load a netcdf file
26	if demoIO;
27	dir_nc='release1/MITprof_release1/';
28	ncload([dir_nc 'argo_feb2013_1992_to_2007_model.nc']);%TBE
29	clear prof_*;
30	end;
31
32	%%demo: MITprof_load and MITprof_subset
33	if demoIO;
34	dir_nc='release1/MITprof_release1/';
35	tmp1=MITprof_load([dir_nc 'argo_feb2013_1992_to_2007_model.nc']);%TBE
36
37	tmp2=MITprof_subset(tmp1,'lon',box_l,'lat',box_L,'depth',box_D);%TBE
38	figure; set(gca,'FontSize',12);
39	plot(tmp1.prof_lon,tmp1.prof_lat,'.');
40	hold on; plot(tmp2.prof_lon,tmp2.prof_lat,'r.');
41	axis([-180 180 -90 90]); title('selected profile locations');
42	xlabel('longitude'); ylabel('latitude')
43	clear tmp*;
44	end;
45
46	%%load complete data set for select variables
47	dir_nc='release1/MITprof_release1/';
48	list_nc={'argo_feb2013_1992_to_2007_model.nc','argo_feb2013_2008_to_2010_model.nc','argo_feb2013_2011_to_2012_model.nc',...
49	'climode_feb2013_model.nc','itp_feb2013_model.nc','seals_feb2013_model.nc'};
50	%note: ctd_feb2013_model.nc and xbt_feb2013_model.nc are left out here for simplicity
51	%note: it takes a couple minutes to load all of the data
52
53	list_var={'prof_T','prof_Testim','prof_Tweight',...
54	'prof_S','prof_Sestim','prof_Sweight'};%TBE
55
56	for vv=1:length(list_var);
57	tmp1=MITprof_stats_load(dir_nc,list_nc,'',list_var{vv});%TBE
58	tmp1=MITprof_subset(tmp1,'lon',box_l,'lat',box_L,'depth',box_D);%TBE
59	eval([list_var{vv} '=tmp1;']);%TBE
60	end;
61
62	clear tmp1 vv;
63