course-idma2016/matlab/idma_box_subset.m


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

%%demo: load a netcdf file
dir_nc='/net/harbor/raid6/ecco/ecco_for_las/version_4/release1/MITprof_inputs/';
ncload([dir_nc 'argo_feb2013_1992_to_2007.nc']);%TBE
clear;

%%set box location
box_l=[-180 -140]; box_L=[-10 10]-30; box_D=[0 300];%TBE
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: MITprof_load and MITprof_subset
dir_nc='/net/harbor/raid6/ecco/ecco_for_las/version_4/release1/MITprof_inputs/';
tmp1=MITprof_load([dir_nc 'argo_feb2013_1992_to_2007.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*;

%%load complete data set for select variables
dir_nc='/net/harbor/raid6/ecco/ecco_for_las/version_4/release1/MITprof_inputs/';
list_nc={'argo_feb2013_1992_to_2007.nc','argo_feb2013_2008_to_2010.nc','argo_feb2013_2011_to_2012.nc',...
         'climode_feb2013.nc','itp_feb2013.nc','seals_feb2013.nc'};
%note: ctd_feb2013.nc and xbt_feb2013.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	gforget	1.1
2			%%add MITprof routines to matlab path:
3			dir_MITprof=[pwd filesep 'c65r_MITprof' filesep];
4			p = genpath(dir_MITprof); addpath(p);%TBE
5
6			%%demo: load a netcdf file
7			dir_nc='/net/harbor/raid6/ecco/ecco_for_las/version_4/release1/MITprof_inputs/';
8			ncload([dir_nc 'argo_feb2013_1992_to_2007.nc']);%TBE
9			clear;
10
11			%%set box location
12			box_l=[-180 -140]; box_L=[-10 10]-30; box_D=[0 300];%TBE
13			fprintf('\n box definition: \n');
14			fprintf(' longitude range: %d %d \n',box_l);
15			fprintf(' latitude range: %d %d \n',box_L);
16			fprintf(' depth range: %d %d \n\n',box_D);
17
18			%%demo: MITprof_load and MITprof_subset
19			dir_nc='/net/harbor/raid6/ecco/ecco_for_las/version_4/release1/MITprof_inputs/';
20			tmp1=MITprof_load([dir_nc 'argo_feb2013_1992_to_2007.nc']);%TBE
21
22			tmp2=MITprof_subset(tmp1,'lon',box_l,'lat',box_L,'depth',box_D);%TBE
23			figure; set(gca,'FontSize',12);
24			plot(tmp1.prof_lon,tmp1.prof_lat,'.');
25			hold on; plot(tmp2.prof_lon,tmp2.prof_lat,'r.');
26			axis([-180 180 -90 90]); title('selected profile locations');
27			xlabel('longitude'); ylabel('latitude')
28			clear tmp*;
29
30			%%load complete data set for select variables
31			dir_nc='/net/harbor/raid6/ecco/ecco_for_las/version_4/release1/MITprof_inputs/';
32			list_nc={'argo_feb2013_1992_to_2007.nc','argo_feb2013_2008_to_2010.nc','argo_feb2013_2011_to_2012.nc',...
33			'climode_feb2013.nc','itp_feb2013.nc','seals_feb2013.nc'};
34			%note: ctd_feb2013.nc and xbt_feb2013.nc are left out here for simplicity
35			%note: it takes a couple minutes to load all of the data
36
37			list_var={'prof_T','prof_Testim','prof_Tweight',...
38			'prof_S','prof_Sestim','prof_Sweight'};%TBE
39
40			for vv=1:length(list_var);
41			tmp1=MITprof_stats_load(dir_nc,list_nc,'',list_var{vv});%TBE
42			tmp1=MITprof_subset(tmp1,'lon',box_l,'lat',box_L,'depth',box_D);%TBE
43			eval([list_var{vv} '=tmp1;']);%TBE
44			end;
45
46			clear tmp1 vv;
47