course-idma2016/matlab/idma_resample.m

function [profIn]=idma_resample(profIn);
%
% IDMA_RESAMPLE subsamples ECCO temperature fields to observed locations and times
%
% Example: 
%    [p]=idma_float_plot('4900828');
%    grid_load;
%    p=idma_resample(p);

p = genpath('gcmfaces/'); addpath(p);%TBE
p = genpath('MITprof/'); addpath(p);%TBE
p = genpath('matlab/'); addpath(p);%TBE

fil='release1/nctiles_climatology/THETA';
fld=read_nctiles([fil '/THETA'],'THETA');
THETA=mean(fld,4);
THETAfil=fil;

fldIn.tim='const';
fldIn.fil=THETAfil;
fldIn.name='prof_THETAconst';
fldIn.fld=THETA;
profIn=MITprof_resample(profIn,fldIn);

fldIn.tim='monclim';
fldIn.fil=THETAfil;
fldIn.name='prof_THETAclim';
fldIn.fld=[];
profIn=MITprof_resample(profIn,fldIn);

if isdir('release1/nctiles/THETA/');
fldIn.tim='monser';
fldIn.fil='release1/nctiles/THETA';
fldIn.name='prof_THETA';
fldIn.fld=[];
profIn=MITprof_resample(profIn,fldIn);
end;

1	gforget	1.2	function [profIn]=idma_resample(profIn);
2			%
3			% IDMA_RESAMPLE subsamples ECCO temperature fields to observed locations and times
4			%
5			% Example:
6			% [p]=idma_float_plot('4900828');
7			% grid_load;
8			% p=idma_resample(p);
9	gforget	1.1
10			p = genpath('gcmfaces/'); addpath(p);%TBE
11			p = genpath('MITprof/'); addpath(p);%TBE
12			p = genpath('matlab/'); addpath(p);%TBE
13
14	gforget	1.2	fil='release1/nctiles_climatology/THETA';
15	gforget	1.1	fld=read_nctiles([fil '/THETA'],'THETA');
16			THETA=mean(fld,4);
17			THETAfil=fil;
18
19			fldIn.tim='const';
20			fldIn.fil=THETAfil;
21			fldIn.name='prof_THETAconst';
22			fldIn.fld=THETA;
23			profIn=MITprof_resample(profIn,fldIn);
24
25			fldIn.tim='monclim';
26			fldIn.fil=THETAfil;
27			fldIn.name='prof_THETAclim';
28			fldIn.fld=[];
29			profIn=MITprof_resample(profIn,fldIn);
30
31	gforget	1.2	if isdir('release1/nctiles/THETA/');
32			fldIn.tim='monser';
33			fldIn.fil='release1/nctiles/THETA';
34			fldIn.name='prof_THETA';
35			fldIn.fld=[];
36			profIn=MITprof_resample(profIn,fldIn);
37			end;
38	gforget	1.1