-
Notifications
You must be signed in to change notification settings - Fork 1
/
annualmean.m
41 lines (32 loc) · 1.05 KB
/
annualmean.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
function indicator = annualmean(data, var, output)
% annual mean for variable "var"
dates = data.dates;
min_year = min(dates(:, 1));
max_year = max(dates(:, 1));
fields = fieldnames(data);
time_dim = (size(data.(fields{1})) == size(dates,1));
for ifield = 1:length(fields)
if ~strcmp(fields{ifield}, 'dates')
idx = ifield;
end
end
indicator.data = NaN(size(min_year:max_year, 2), size(data.(fields{idx}), find(~time_dim)));
indicator.dates = NaN(size(min_year:max_year, 1));
if isfield(data, var)
F = getfield(data, var);
else
error(strcat('this indicator requires a data struct with .', var))
end
for y = min_year:max_year
indicator.data(y-min_year +1, :) = squeeze(mean(F(dates(:, 1) == y, :), 1));
indicator.dates(y - min_year +1, :) = y;
end
indicator.units = 'deg C';
indicator.var = strcat(['Annual mean of ', var]);
indicator.title = strcat(['Annual mean of ', var]);
if exist('output','var')
assert(strcmp(output, 'matrix'))
% output parameter exist, so output only the data matrix
indicator = indicator.data;
end
end