;
restore,filename='reglists.idlsave'
;
for i = 0 , nreg-1 do begin
  openw,1,'ts_'+regname(i)+'.out'
  ;
  ; Get data in
  ;
  restore,filename='instradj_'+regname(i)+'.idlsave'
  restore,filename='rwidadj_'+regname(i)+'.idlsave'
  restore,filename='densadj_'+regname(i)+'.idlsave'
  nyr=n_elements(densadj)
  ;
  ; Do ranking over available data
  ;
  keeplist=where(finite(densadj),nkeep)
  v1=densadj(keeplist)
  s1=sort(v1)
  rank1=intarr(nyr)-999
  k=0
  for j = 0 , nyr-1 do begin
    if finite(densadj(j)) then begin
      rank1(j)=where(s1 eq k)+1
      k=k+1
    endif
  endfor
  keeplist=where(finite(rwidadj),nkeep)
  v1=rwidadj(keeplist)
  s1=sort(v1)
  rank2=intarr(nyr)-999
  k=0
  for j = 0 , nyr-1 do begin
    if finite(rwidadj(j)) then begin
      rank2(j)=where(s1 eq k)+1
      k=k+1
    endif
  endfor
  keeplist=where(finite(instradj),nkeep)
  v1=instradj(keeplist)
  s1=sort(v1)
  rank3=intarr(nyr)-999
  k=0
  for j = 0 , nyr-1 do begin
    if finite(instradj(j)) then begin
      rank3(j)=where(s1 eq k)+1
      k=k+1
    endif
  endfor
  ;
  ; Set to missing flag
  ;
  missval=-9.999
  misslist=where(finite(densadj) eq 0,nmiss)
  if nmiss gt 0 then densadj(misslist)=missval
  misslist=where(finite(rwidadj) eq 0,nmiss)
  if nmiss gt 0 then rwidadj(misslist)=missval
  misslist=where(finite(instradj) eq 0,nmiss)
  if nmiss gt 0 then instradj(misslist)=missval
  ;
  ; Output it
  ;
  printf,1,'Timeseries (variance-adjusted & normalised over 1901-70) for region '+regname(i)
  printf,1,'      | Densities     | Widths        |        | Temperatures'
  printf,1,' Year | Rank  Anomaly | Rank  Anomaly | Nchron | Rank  Anomaly'
  for j = 0 , nyr-1 do begin
    printf,1,x(j),rank1(j),densadj(j),rank2(j),rwidadj(j),n(j),rank3(j),instradj(j), $
          format='(I5," |",I5,F9.3," |",I5,F9.3," |",I7," |",I5,F9.3)'
  endfor
  close,1
endfor
;
end