c=readfromfile('co2.recent.txt') y=c(0,*) n=n_elements(y) nym=n*12 b=reform(c(1+indgen(12),*),nym) ym=fltarr(nym) i=0 for y=y(0),y(n-1) do for m=1,12 do begin ym(i)=y+(m-1)/12. I=I+1 endfor plot,/xs,/ynoz,ym,b sm=6 j=indgen(nym-2*sm)+sm oplot,ym(j),(smooth(b,2*sm+1))(j),th=3 gettwogifs,out='~/np/wmc/a' end