;
; Outputs regional timeseries
;
; Get region info
;
restore,filename='reglists.idlsave'
regfn=['MEAN','n60','s60','wna','ena','cer','neu','seu']
missval=-99.9999
;
; For each region
;
;nreg=1
regname(0)='MEAN'
regfn(0)='mean'
for ireg = 0 , nreg-1 do begin
  ;
  ; Get adjusted density timeseries, adjusted width and actual temperatures
  ;
  print,'densadj_'+regname(ireg)+'.idlsave'
  restore,filename='densadj_'+regname(ireg)+'.idlsave'
  nt=n_elements(x)
  ;
  misslist=where(finite(densadj) eq 0,nmiss)
  if nmiss gt 0 then densadj(misslist)=missval
  ;
  openw,1,regfn(ireg)+'_mxd.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  MXD value'
  for iyr = 0 , nt-1 do begin
    printf,1,x(iyr),densadj(iyr),format='(I5,F9.4)'
  endfor
  close,1
  ;
endfor
;
end