; ; Plots correlations between regional temperatures and regional MXD series ; based on varying numbers of stations. ; restore,filename='regbest_fixed1950.idlsave' ; nreg,regname,nsites,regr,locr,bestx,besty,bestmxd,regtemp,timey ; fixedyr,corrper,doreg,maxntree ; ; Prepere for plotting ; loadct,39 multi_plot,nrow=6,ncol=2,layout='large' !p.multi(4)=1 if !d.name eq 'X' then begin window,ysize=850 !p.font=-1 endif else begin !p.font=0 device,/helvetica,/bold,font_size=9 endelse ; ; Now plot each region separately ; for ireg = 0 , nreg-1 do begin ; pause plot,locr(*,ireg),thick=3,$ title=regname(ireg),$ /xstyle,xtickformat='nolabels',$ /ystyle,yrange=[-0.45,0.95],ymargin=[0,2],ytitle='Local correlation' for j = -4 , 8 , 2 do oplot,!x.crange,replicate(j*0.1,2),linestyle=(j ne 0) plot,regr(*,ireg,0),thick=3,$ /xstyle,xtitle='Chronology',$ /ystyle,yrange=[-0.45,0.95],ymargin=[4,0],ytitle='Regional correlation' for j = -4 , 8 , 2 do oplot,!x.crange,replicate(j*0.1,2),linestyle=(j ne 0) oplot,regr(*,ireg,1),thick=6,linestyle=2 oplot,regr(*,ireg,2),thick=6 ; endfor ; end