; ; Plots the final long volcano record (3 series mean) ; multi_plot,nrow=5 if !d.name eq 'X' then begin window,ysize=900 initx endif else begin device,xoffset=2,xsize=17 endelse ; restore,'longvolcadj.idlsave' ; keeplist=where(x ge 1000) x=x(keeplist) densadj=densadj(keeplist) ; ;filter_cru,20,tsin=densadj,tslow=tslow,/nan cpl_barts,x,densadj,$ xtitle='Year',/xstyle,$ ytitle='Normalised adjusted density record',$ zeroline=0.,yrange=[-5,5],/ystyle ;oplot,x,tslow,thick=2 oplot,!x.crange,[0.,0.] oplot,!x.crange,[-2.,-2.],linestyle=1 oplot,!x.crange,[2.,2.],linestyle=1 ; end