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

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

PS=squeeze(mean(ps)); %make climatology

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

        