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	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
10	p = genpath('gcmfaces/'); addpath(p);%TBE
11	p = genpath('MITprof/'); addpath(p);%TBE
12	p = genpath('matlab/'); addpath(p);%TBE
13
14	fil='release1/nctiles_climatology/THETA';
15	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	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