;
loadct,39
map=def_map(/npolar)
map.limit(0)=30.
labels=def_labels(/off)
labels.gridon=1  &  labels.dlon=90.
;sm=def_sm(/off)
sm=def_sm()
sm.thresh=0.1
;sm.method=0
;
print,'Enter year required:'
read,iyr
densfd=rd1yr(iyr,x=x,y=y,frac=frac)
;
; Remove missing data and data south of 30N
;
keeplist=where((densfd ne -9.99) and (y ge 30.),nkeep)
if nkeep gt 0 then begin
  frac=frac(keeplist)
  densfd=densfd(keeplist)
  x=x(keeplist)
  y=y(keeplist)
  inter_const,densfd,x,y,map=map,$
    maxdist=3000.,$
    gs=[2.,2.],$
    sm=sm,labels=labels,$
    shade=1,sh_thresh=0.,$
    levels=findgen(21)*0.5-5.,/follow,$
    title=string(iyr,format='(I4)'),$
    xtitle='Tree ring density anomaly'
endif
;
end