; ; Now prepare for plotting ; loadct,39 multi_plot,nrow=3,layout='caption' if !d.name eq 'X' then begin window,ysize=800 !p.font=-1 endif else begin !p.font=0 device,/helvetica,/bold,font_size=18 endelse def_1color,20,color='red' def_1color,21,color='blue' def_1color,22,color='black' ; restore,'compbest_fixed1950.idlsave' ; plot,timey,comptemp(*,3),/nodata,$ /xstyle,xrange=[1881,1994],xtitle='Year',$ /ystyle,yrange=[-3,3],ytitle='Normalised anomalies',$ ; title='Northern Hemisphere temperatures, MXD and corrected MXD' title='Northern Hemisphere temperatures and MXD reconstruction' ; yyy=reform(comptemp(*,2)) ;mknormal,yyy,timey,refperiod=[1881,1940] filter_cru,5.,/nan,tsin=yyy,tslow=tslow oplot,timey,tslow,thick=5,color=22 yyy=reform(compmxd(*,2,1)) ;mknormal,yyy,timey,refperiod=[1881,1940] ; ; Apply a VERY ARTIFICAL correction for decline!! ; yrloc=[1400,findgen(19)*5.+1904] valadj=[0.,0.,0.,0.,0.,-0.1,-0.25,-0.3,0.,-0.1,0.3,0.8,1.2,1.7,2.5,2.6,2.6,$ 2.6,2.6,2.6]*0.75 ; fudge factor if n_elements(yrloc) ne n_elements(valadj) then message,'Oooops!' ; yearlyadj=interpol(valadj,yrloc,timey) ; ;filter_cru,5.,/nan,tsin=yyy+yearlyadj,tslow=tslow ;oplot,timey,tslow,thick=5,color=20 ; filter_cru,5.,/nan,tsin=yyy,tslow=tslow oplot,timey,tslow,thick=5,color=21 ; oplot,!x.crange,[0.,0.],linestyle=1 ; plot,[0,1],/nodata,xstyle=4,ystyle=4 ;legend,['Northern Hemisphere April-September instrumental temperature',$ ; 'Northern Hemisphere MXD',$ ; 'Northern Hemisphere MXD corrected for decline'],$ ; colors=[22,21,20],thick=[3,3,3],margin=0.6,spacing=1.5 legend,['Northern Hemisphere April-September instrumental temperature',$ 'Northern Hemisphere MXD'],$ colors=[22,21],thick=[3,3],margin=0.6,spacing=1.5 ; end