; ; Plots changes in correlations between regional MXD and regional monthly & ; seasonal precip. and temperature. Partial correlations for precip., full ; for temp. ; for trv = 0 , 1 do begin case trv of 0: begin fnadd='mxd' titadd='(a)' end 1: begin fnadd='trw' titadd='(b)' end endcase ; ; Restore all correlations ; restore,filename='reg'+fnadd+'_moncorr_early.idlsave' ; allr,ncorr,nvar,nreg,regname,varname,corrname,allp,lowr,lowp,thalf allr1=lowr allp1=lowp restore,filename='reg'+fnadd+'_moncorr_late.idlsave' ; allr,ncorr,nvar,nreg,regname,varname,corrname,allp,lowr,lowp,thalf allr2=lowr allp2=lowp restore,filename='reg'+fnadd+'_moncorr.idlsave' ; allr,ncorr,nvar,nreg,regname,varname,corrname,allp,lowr,lowp,thalf ; ; Combine temp r with prec partial r ; regr=allr regr(*,*,0)=lowr(*,*,1) regr(*,*,1)=lowp(*,*,0) varname=reverse(varname) conname=['',' constant Temp'] earr=allr1 earr(*,*,0)=allr1(*,*,1) earr(*,*,1)=allp1(*,*,0) latr=allr2 latr(*,*,0)=allr2(*,*,1) latr(*,*,1)=allp2(*,*,0) ; ; Now generate three extra composite regions ; newr=fltarr(nreg+3,ncorr,nvar) newe=fltarr(nreg+3,ncorr,nvar) newl=fltarr(nreg+3,ncorr,nvar) newr(0:nreg-1,*,*)=regr(*,*,*) newe(0:nreg-1,*,*)=earr(*,*,*) newl(0:nreg-1,*,*)=latr(*,*,*) regname=[regname,'HILAT','LOLAT','NH'] ; rlist=[0,2,3,7,8] newr(nreg,*,*)=total(regr(rlist,*,*),1)/float(n_elements(rlist)) newe(nreg,*,*)=total(earr(rlist,*,*),1)/float(n_elements(rlist)) newl(nreg,*,*)=total(latr(rlist,*,*),1)/float(n_elements(rlist)) ; rlist=[1,4,5,6] newr(nreg+1,*,*)=total(regr(rlist,*,*),1)/float(n_elements(rlist)) newe(nreg+1,*,*)=total(earr(rlist,*,*),1)/float(n_elements(rlist)) newl(nreg+1,*,*)=total(latr(rlist,*,*),1)/float(n_elements(rlist)) ; newr(nreg+2,*,*)=total(regr,1)/float(nreg) newe(nreg+2,*,*)=total(earr,1)/float(nreg) newl(nreg+2,*,*)=total(latr,1)/float(nreg) ; regr=newr earr=newe latr=newl nreg=nreg+3 ; ; Now prepare for plotting ; if trv eq 0 then begin loadct,39 multi_plot,nrow=3 if !d.name eq 'X' then begin window,ysize=800 !p.font=-1 endif else begin !p.font=0 device,/helvetica,/bold,font_size=12,xsize=8. endelse def_1color,20,color='mgrey' def_1color,21,color='palegrey' def_1color,22,color='orange' !x.margin=[8,0] !y.margin=[7,3] endif ; for i = 9 , 9 do begin for ivar = 0 , nvar-2 do begin if ivar eq 1 then yta='Partial ' else yta=' ' if ivar eq 0 then $ ymt='Correlations between '+titadd+' and monthly/seasonal TEMPERATURES!C'+$ 'High-latitude regions!C1901-1950 (BARS) 1935-1984 (DOTS)' $ else $ ymt='Partial correlations between '+titadd+' and monthly/seasonal PRECIPITATION!C'+$ 'High-latitude regions!C1901-1950 (BARS) 1935-1984 (DOTS)' ymt='(All time series decadally-filtered)!C'+ymt x=findgen(28) yy=reform(earr(i,*,ivar)) y=fltarr(28)*!values.f_nan yloc=[indgen(16)+1,indgen(5)*2+18] y(yloc)=yy(*) plot,x,y,/nodata,$ /ystyle,yrange=[-0.5,1.],ytitle=yta+'Correlation',$ xstyle=5,xrange=[7.,27.],title=titadd ;ymt ; if !d.name eq 'PS' then device,font_size=14 for j = -10,10,2 do oplot,!x.crange,replicate(float(j)*0.1,2),linestyle=1 oplot,!x.crange,[0.,0.] ; oplot,[7.5,7.5],!y.crange,thick=4 oplot,[17.,17.],!y.crange,thick=4 cpl_barts,x(8:27),y(8:27),outline=-1,/overplot,bar_color=[21,20],thick=2 if !d.name eq 'PS' then device,font_size=9 for j = 7 , ncorr-1 do begin x1=x(yloc(j)) y1ear=earr(i,j,ivar) y1lat=latr(i,j,ivar) ; plots,x1,y1ear,psym=def_sym(6) plots,x1,y1lat,psym=def_sym(10),symsize=1.3 x1=x(yloc(j))+0.28 y1=!y.crange(0)-0.05 xyouts,x1,y1,corrname(j),align=1,orient=90. endfor ; xyouts,7.,-0.85,align=1.,'Previous!Cyear' ; xyouts,8.,-0.85,align=0.,'Same!Cyear' xyouts,17.5,-0.35,align=0.,'Seasons' if !d.name eq 'PS' then device,font_size=12 ;!y.margin=[3,7] endfor endfor ; endfor ; end