64 |
|
|
65 |
for vv=1:length(listVar); |
for vv=1:length(listVar); |
66 |
varCur=listVar{vv}; |
varCur=listVar{vv}; |
67 |
|
|
68 |
%get the cost square root: |
%get the cost square root: |
69 |
[MITprof]=MITprof_stats_load([dirData],listData,varCur); |
[MITprof]=MITprof_stats_load([dirData],listData,varCur); |
70 |
|
if ~isfield(MITprof,'prof_point'); |
71 |
|
loc_tile=gcmfaces_loc_tile(90,90,MITprof.prof_lon,MITprof.prof_lat); |
72 |
|
MITprof.prof_point=loc_tile.point; |
73 |
|
end; |
74 |
%mask out values that are not in year range: |
%mask out values that are not in year range: |
75 |
ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1); |
ii=find(MITprof.prof_date<date0|MITprof.prof_date>date1); |
76 |
MITprof.prof(ii,:)=NaN; |
MITprof.prof(ii,:)=NaN; |