; ; 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