%build the surface U from the individual matlab saved decadal means
off=0;
for n=1:36;
n
fil=['U' num2str(n) '.mat'];
load(['F:\trace\U\mat_dec\' fil], 'lat','lon', 'U1');
[NT, NN,N]=size(U1);

U(off+1:off+NT, 1:NN, 1:N)=U1;
NTNT(n)=NT;
clear U1
off=off+NT;
end

save('F:\trace\U\U_dec.mat','U', 'lat', 'lon')