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

grid_load; gcmfaces_global; MITprof_global;

profIn=idma_float_plot('4900828');
[tmp1,tt]=sort(profIn.prof_date);

fil=[pwd '/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);