as=readfromfile('$WMCDATA/a-s.txt',ign=1,col=13)	; South Pole data, from READER
as1=as(indgen(12)+1,*)					; Extract monthly data, removing year
as2=deseasonalise(as1)					; Deseasonalise

t=as1							; Make time info
for i=0,n_elements(t(0,*))-1 do t(*,i)=as(0,i)+indgen(12)/12.

; Now we'll plot various looks at the dataset
!p.multi=[0,2,6]

; First, the raw data. With seasonal cycle, its not useful...
lund_plot,as1,tis=t,/nop

; Second, deseasonalised
lund_plot,as2,tis=t,/nop

; Thirs, annual mean
lund_plot,avg(as2,0),tis=as(0,*),/nop

; Fourth, Jan
lund_plot,as1(0,*),tis=as(0,*),/nop

; Fifth, JJA
lund_plot,avg(as1([6,7,8]-1,*),0),tis=as(0,*),/nop

; Sixth, MAM
lund_plot,avg(as1([3,4,5]-1,*),0),tis=as(0,*),/nop

gettwogifs,out='ex1'