;
; Computes correlations on full, high and low pass timeseries of density
; anomalies, ring width anomalies and instrumental anomalies.
; NOW USES DEC-MAR TEMPERATURES
;
restore,filename='reglists.idlsave'
;
; Read in timeseries and compute correlations
;
perst=1881.
peren=1960.
;
openw,1,'corr_djfm2all.out'
thalf=10.
printf,1,'Correlations between timeseries'
printf,1,'TEMPERATURES ARE DEC-JAN-FEB-MAR'
printf,1
printf,1,' TRW-MXD vs. instr density vs. instru width vs. instru'
printf,1,' Region Full High Low Full High Low Full High Low'
;
for i = 0 , nreg-1 do begin
print,i
;
restore,filename='densadj_'+regname(i)+'.idlsave'
y1=densadj
filter_cru,thalf,tsin=y1,tslow=ylow1,tshigh=yhi1,/nan
;
restore,filename='rwidadj_'+regname(i)+'.idlsave'
y2=rwidadj
filter_cru,thalf,tsin=y2,tslow=ylow2,tshigh=yhi2,/nan
;
restore,filename='instrdjfm_'+regname(i)+'.idlsave'
y3=instrdjfm
filter_cru,thalf,tsin=y3,tslow=ylow3,tshigh=yhi3,/nan
;
keeplist=where(finite(y1+y2+y3) and (x ge perst) and (x le peren),nkeep)
ts1=fltarr(2,nkeep)
ts1(0,*)=y1(keeplist)
ts1(1,*)=y2(keeplist)
ts2=y3(keeplist)
r1=m_correlate(ts1,ts2)
;
ts1=fltarr(2,nkeep)
ts1(0,*)=yhi1(keeplist)
ts1(1,*)=yhi2(keeplist)
ts2=yhi3(keeplist)
r2=m_correlate(ts1,ts2)
;
ts1=fltarr(2,nkeep)
ts1(0,*)=ylow1(keeplist)
ts1(1,*)=ylow2(keeplist)
ts2=ylow3(keeplist)
r3=m_correlate(ts1,ts2)
;
keeplist=where(finite(y1+y3) and (x ge perst) and (x le peren),nkeep)
ts1=y1(keeplist)
ts2=y3(keeplist)
r4=correlate(ts1,ts2)
;
ts1=yhi1(keeplist)
ts2=yhi3(keeplist)
r5=correlate(ts1,ts2)
;
ts1=ylow1(keeplist)
ts2=ylow3(keeplist)
r6=correlate(ts1,ts2)
;
keeplist=where(finite(y2+y3) and (x ge perst) and (x le peren),nkeep)
ts1=y3(keeplist)
ts2=y2(keeplist)
r7=correlate(ts1,ts2)
;
ts1=yhi3(keeplist)
ts2=yhi2(keeplist)
r8=correlate(ts1,ts2)
;
ts1=ylow3(keeplist)
ts2=ylow2(keeplist)
r9=correlate(ts1,ts2)
;
printf,1,regname(i),r1,r2,r3,r4,r5,r6,r7,r8,r9,$
format='(A11,3(2X,3F6.2))'
;
endfor
;
printf,1,' '
printf,1,'Correlations carried out over the period ',perst,peren
printf,1,' '
printf,1,'To separate low and high frequency components, a gaussian weighted'
printf,1,'filter was used with a half-width (years) of ',thalf
;
close,1
;
end