%make the climatological surface pressure
clear all; close all;
start_year=1970;
end_year=2008;

for yr=start_year:1:end_year;
    tic
    yr
    for mn=1:12;
        load([ 'E:\data\data\means\' num2str(yr) '_' num2str(mn) '_means.mat'], 'lat','lon','Plev' ,'T');
        NN=length(lat); N=length(lon);
        NL=length(Plev);
        t(yr-start_year+1, mn,1:NL, 1:NN,1:N)= T;
    end
    toc
end

T=squeeze(mean(t)); %make climatology

save('E:\data\data\means\T.mat', 'T','Plev', 'lat', 'lon')

        