; ; Plots the final long volcano record (3 series mean) ; loadct,39 def_1color,r,g,b,10,color='red' def_1color,r,g,b,11,color='dblue' multi_plot,nrow=1 if !d.name eq 'X' then begin window,ysize=300,xsize=700 initx endif else begin device,yoffset=8,ysize=8 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,bar_color=[10,11] ;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