; restore,filename='reglists.idlsave' ; ; Ignore 1991 and 1992 (last two values) from all timeseries ; for i = 0 , nreg-1 do begin restore,filename='densts_'+regname(i)+'.idlsave' keeplist=where(finite(wdens) and (x lt 1991),nkeep) y=wdens(keeplist) x=x(keeplist) print,'Moments of normalised density anomalies for region: '+$ regname(i) print,' Mean Variance Skewness Kurtosis' print,moment(y),format='(4F10.4)' print,'----------------------------------------' endfor ; end