; 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