format compact;clear all

var = 'PSL';

load([char(var) '.mat']);

X1 = X1/100;
X1_DJF = X1_DJF/100;
X1_JJA = X1_JJA/100;
X3 = X3/100;
X3_DJF = X3_DJF/100;
X3_JJA = X3_JJA/100;

figure;orient tall;orient landscape;colormap(ctb)
v=linspace(-5,5,41);
subplot('position',[.1 .5 .4 .3])
tmp=reshape(mean(X3-X1,1),im*jm,1);
tmp=max(tmp,-5);
tmp=min(tmp,5);
[c,h]=contourf(lon,lat,(reshape(tmp,im,jm))',v);
set(gca,'clim',[min(v) max(v)]);
set(h,'edgecolor','none')
hold on
plot(lon1,lat1,'k-')
[cc,hh]=contour(lon,lat,squeeze(mean(X1,1))',[900:10:1010],'b--','ShowText','on');
[cc,hh]=contour(lon,lat,squeeze(mean(X1,1))',[1020:10:1050],'b-','ShowText','on');
title(['(a) \Delta' char(var) ' (hPa) ANN'])
set(gca,'xtick',[0:60:360],'xminortick','on','xticklabel',[])
set(gca,'ytick',[-90:30:90],'yminortick','on','yticklabel',[])
box on

h=colorbar('horiz');
set(h,'position',[0.18 0.47 0.25 0.01])

var = 'TMQ';
load([char(var) '.mat']);


subplot('position',[.52 .5 .4 .3])

var = 'TMQ';
load([char(var) '.mat']);

v=[-10 -6 -4 -2 -1.5 -1 -.8 -.6 -.4 -.2 0 .2 .4 .6 .8 1 1.5 2 4 6 10]; 
tmp=reshape(mean(X3-X1,1),im*jm,1);
tmp=max(tmp,-10);
tmp=min(tmp,10);
[c,h]=contourf(lon,lat,(reshape(tmp,im,jm))',v);
set(gca,'clim',[min(v) max(v)]);
set(h,'edgecolor','none')
hold on
plot(lon1,lat1,'k-')
title(['(b) \Delta' char(var) ' (kg/m2) ANN'])
set(gca,'xtick',[0:60:360],'xminortick','on','xticklabel',[])
set(gca,'ytick',[-90:30:90],'yminortick','on','yticklabel',[])
box on

h=colorbar('horiz');
set(h,'position',[0.6 0.47 0.25 0.01])