;
; Outputs regional timeseries
;
; Get region info
;
restore,filename='reglists.idlsave'
regfn=['all','n60','s60','wna','ena','cer','neu','seu']
missval=-99.9999
;
; For each region
;
for ireg = 0 , nreg-1 do begin
  ;
  ; Get adjusted density timeseries, adjusted width and actual temperatures
  ;
  restore,filename='densadj_'+regname(ireg)+'.idlsave'
  restore,filename='rwidadj_'+regname(ireg)+'.idlsave'
  restore,filename='nc.instrts_'+regname(ireg)+'.idlsave'
  nt=n_elements(x)
  ;
  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(winstr) eq 0,nmiss)
  if nmiss gt 0 then winstr(misslist)=missval
  ;
  openw,1,regfn(ireg)+'_tdwi.dat'
  printf,1,regname(ireg),format='(A15,5X,"Region name")'
  printf,1,nt,format='(I5,15X,"Number of years")'
  printf,1,missval,format='(F9.4,11X,"Missing value")'
  printf,1,' Year  Density    Width    April      May     June     July   August     Sept'
  for iyr = 0 , nt-1 do begin
    printf,1,x(iyr),densadj(iyr),rwidadj(iyr),reform(winstr(iyr,3:8)),$
      format='(I5,8F9.4)'
  endfor
  close,1
  ;
  openw,1,regfn(ireg)+'_temp.dat'
  printf,1,regname(ireg),format='(A15,5X,"Region name")'
  printf,1,nt-256,format='(I5,15X,"Number of years")'
  printf,1,missval,format='(F9.4,11X,"Missing value")'
  printf,1,' Year  Jan-Dec'
  for iyr = 256 , nt-1 do begin
    printf,1,x(iyr),reform(winstr(iyr,0:11)),$
      format='(I5,12F9.4)'
  endfor
  close,1
  ;
endfor
;
end