;
; Plots various versions of the NH reconstruction using age-banded.
;
nver=9
doreg='ALL'
docol=0          ; 0=b&w 1=col
;
; First, allsites
;
restore,filename='bandalltemp_calibrated.idlsave'
; Gets: nyr,nhtit,timey,nhtemp
if nhtit ne doreg then message,'Re-calibrate the age-banded!'
;
nyrnh=nyr
timeynh=timey
allts=fltarr(nyrnh,nver) & allts(*,*)=!values.f_nan
allts(*,0)=nhtemp(*)
;
; Next, various versions of the regional combination
;
for iver = 2 , 7 do begin
  ;
  fn='ver'+string(iver,format='(I1)')+'_nhrecon.idlsave'
  print,fn
  restore,filename=fn
  ; Gets: calregts,nreg,nyr,regname,tempnyr,temptimey,tempregts,timey
;  if regname(nreg-1) ne doreg then message,'Re-calibrate the age-banded!'
  if regname(nreg-1) ne doreg then print,regname(nreg-1)
  allts(*,iver-1)=calregts(*,nreg-1)
  ;
endfor
;
; Now for the PCR reconstruction
;
;restore,filename='bandtempNH_calmultipcr.idlsave'
restore,filename='bandtempALL_calmultipcr.idlsave'
; Gets: nyr,nhtit,yrmxd,prednh,fullnh
if nhtit ne doreg then message,'Re-calibrate the age-banded!'
;allts(*,7)=fullnh(*)
allts(*,7)=prednh(*)
;
thalf=25.       ; filter for high or low-pass
;
if docol ne 0 then begin
  lcol=[10,11,12,13,14,15,16,17]
  lthi=[ 2, 4, 4, 4, 4, 4, 4, 4]
endif else begin
  lcol=[10,10,10,10,10,10,10,10]
  lthi=[ 1, 1, 1, 1, 1, 1, 1, 4]
endelse
;
loadct,39
multi_plot,nrow=3,layout='large'
if !d.name eq 'X' then begin
  window,ysize=900
  !p.font=-1
endif else begin
  !p.font=0
  device,/helvetica,bold=0,font_size=14
endelse
def_1color,10,color='black'
def_1color,11,color='mdyellow'
def_1color,12,color='lblue'
def_1color,13,color='orange'
def_1color,14,color='deepblue'
def_1color,15,color='red'
def_1color,16,color='vlpurple'
def_1color,17,color='lgreen'
;
plot,[0,1],/nodata,$
  /xstyle,xrange=[1400,1960],xtitle='Year',$
;  /ystyle,yrange=[-0.85,0.3],ytitle='Temperature anomaly  (!Uo!NC wrt 1961-90)'
  /ystyle,yrange=[-1.1,0.4],ytitle='Temperature anomaly  (!Uo!NC wrt 1961-90)'
oplot,!x.crange,[0.,0.],linestyle=1
;
for iver = 0 , nver-1 do begin
  y=allts(*,iver)
  filter_cru,thalf,/nan,tsin=y,tslow=ylow
  oplot,timeynh,ylow,thick=lthi(iver),color=lcol(iver)
endfor
;
y=tempregts(*,nreg-1)
filter_cru,thalf,/nan,tsin=y,tslow=ylow
oplot,temptimey,ylow,thick=7,linestyle=1
;
if docol ne 0 then begin
  legend,['1','2','3','4','5','6','7','8','Observed'],$
    thick=[lthi,5],$
    color=[lcol,!p.color],/horiz,pspacing=2.5
endif
;
end