[num,txt,raw] = xlsread('A.xls') ;
date = datevec(num2str(num(:,1)),'yyyymmdd') ;
years = date(:,1) ;
months = date(:,2) ;
%%Get average
[Y,ia,ib] = unique(years) ;
[M,ia,ib] = unique(months) ;
iwant1 = zeros(length(Y),length(M)) ; % Average of max T
iwant2 = zeros(length(Y),length(M)) ; % Average of min T
iwant3 = zeros(length(Y),length(M)) ; % Average of PRCP
for i = 1:length(Y)
idx_Y = years==Y(i) ;
data = num(idx_Y,2:end) ;
YM = months(idx_Y) ;
for j = 1:length(M)
idx_M = YM==M(j) ;
iwant1(i,j) = mean(data(idx_M,1)) ;
iwant2(i,j) = mean(data(idx_M,2)) ;
iwant3(i,j) = mean(data(idx_M,3)) ;
end
end
