;
restore,'tchange.idlsave'
;
multi_plot,nrow=2
;
if !d.name eq 'X' then window,ysize=800
;
loadct,39
;
map=def_map(/npolar)  &  map.limit(0)=30.
labels=def_labels(/off)
sm=def_sm() & sm.thresh=0.4
;
levs=[-0.5,-0.4,-0.3,-0.2,-0.1]
levs=[levs,0.,-reverse(levs)]
labs=intarr(n_elements(levs))+1
thr=0.
labs(where(levs eq thr))=0
;
inter_confd,tempdiff,g.x,g.y,levels=levs,/follow,$
  c_labels=labs,$
  /shade,sh_thresh=thr,/hi_on,$
  map=map,labels=labels,sm=sm,fdsm=fdsm
;
restore,'pchange.idlsave'
;
levs=[0.5,0.6,0.7,0.8,0.9,0.95,1.,1.05,1.1,1.2,1.3,1.5,2.]
labs=intarr(n_elements(levs))+1
thr=1.
labs(where(levs eq thr))=0
;
inter_confd,precdiff,g.x,g.y,levels=levs,/follow,$
  c_labels=labs,$
  /shade,sh_thresh=thr,/hi_on,$
  map=map,labels=labels,sm=sm
;
end