1 |
gforget |
1.1 |
|
2 |
|
|
p = genpath('gcmfaces/'); addpath(p);%TBE |
3 |
|
|
p = genpath('MITprof/'); addpath(p);%TBE |
4 |
|
|
p = genpath('matlab/'); addpath(p);%TBE |
5 |
|
|
|
6 |
|
|
grid_load; gcmfaces_global; MITprof_global; |
7 |
|
|
|
8 |
|
|
profIn=idma_float_plot('4900828'); |
9 |
|
|
[tmp1,tt]=sort(profIn.prof_date); |
10 |
|
|
|
11 |
|
|
fil=[pwd '/release1/nctiles_climatology/THETA']; |
12 |
|
|
fld=read_nctiles([fil '/THETA'],'THETA'); |
13 |
|
|
THETA=mean(fld,4); |
14 |
|
|
THETAfil=fil; |
15 |
|
|
|
16 |
|
|
fldIn.tim='const'; |
17 |
|
|
fldIn.fil=THETAfil; |
18 |
|
|
fldIn.name='prof_THETAconst'; |
19 |
|
|
fldIn.fld=THETA; |
20 |
|
|
profIn=MITprof_resample(profIn,fldIn); |
21 |
|
|
|
22 |
|
|
fldIn.tim='monclim'; |
23 |
|
|
fldIn.fil=THETAfil; |
24 |
|
|
fldIn.name='prof_THETAclim'; |
25 |
|
|
fldIn.fld=[]; |
26 |
|
|
profIn=MITprof_resample(profIn,fldIn); |
27 |
|
|
|
28 |
|
|
|