;
restore,filename='reglists.idlsave'
;
; Ignore 1991 and 1992 (last two values) from all timeseries
;
bs=0.1
for i = 0 , nreg-1 do begin
  restore,filename='rwidts_'+regname(i)+'.idlsave'
  keeplist=where(x lt 1991,nkeep)
  y=wrwid(keeplist)
  x=x(keeplist)
  yhist=tim_hist(y,binsize=bs,binvals=binvals)
  yhist=[0.,yhist]  &  binvals=[binvals(0)-bs,binvals]
  pause
  plot,binvals,yhist,psym=10,$
    title='Distribution of normalised ring width anomalies for region: '+$
      regname(i),$
    xtitle='Ring width anomaly',$
    ytitle='Frequency'
endfor
;
end