; ; Plots regional and large-scale averages of correlations between ; monthly & seasonal precip. and temperature. ; ; Restore all correlations ; restore,filename='tp_moncorr.idlsave' ; allr,ncorr,nchron,corrname,statlat,statlon ; ; Restore regional breakdown ; restore,filename='reg_mxdlists.idlsave' ; ntree,treelist,nreg,regname ; ; For each region, average the site-by-site correlations together ; regr=fltarr(nreg,ncorr) for i = 0 , nreg-1 do begin n=ntree(i) oner=allr(treelist(0:n-1,i),*) totr=total(oner,1)/float(n) regr(i,*)=totr(*) endfor ; ; Now generate three extra composite regions ; newr=fltarr(nreg+3,ncorr) newr(0:nreg-1,*)=regr(*,*) regname=[regname,'HILAT','LOLAT','NH'] ; rlist=[0,2,3,7,8] oner=total(regr(rlist,*),1)/float(n_elements(rlist)) newr(nreg,*)=oner(*) ; rlist=[1,4,5,6] oner=total(regr(rlist,*),1)/float(n_elements(rlist)) newr(nreg+1,*)=oner(*) ; oner=total(regr,1)/float(nreg) newr(nreg+2,*)=oner(*) ; regr=newr nreg=nreg+3 ; ; Now prepare for plotting ; loadct,39 multi_plot,nrow=2 if !d.name eq 'X' then begin window,ysize=800 !p.font=-1 endif else begin !p.font=0 device,/helvetica,/bold,font_size=14 endelse def_1color,20,color='deepblue' def_1color,21,color='red' !x.margin=[10,10] ; for i = nreg-3 , nreg-2 do begin x=findgen(12) y=reform(regr(i,0:11)) plot,x,y,/nodata,$ /ystyle,yrange=[-0.6,0.6],ytitle='Correlation',$ xstyle=5,xrange=[-1,12],title=regname(i)+' r(temp,prec)' for j = -6 , 6 , 2 do oplot,!x.crange,replicate(float(j)*0.1,2),linestyle=1 oplot,!x.crange,[0.,0.] cpl_barts,x,y,outline=-1,/overplot,bar_color=[21,20],thick=2 ; if !d.name eq 'PS' then device,font_size=8 for j = 0 , 11 do begin x1=x(j)+0.25 y1=y(j) if y1 lt 0 then y1=+0.05 else y1=y1+0.05 xyouts,x1,y1,corrname(j),align=0.,orient=90. endfor ; if !d.name eq 'PS' then device,font_size=14 endfor ; end