;
; Plots maps of chronology location, tree type, country
;
ncid=ncdf_open('tree_dens_nh.nc')
ncdf_diminq,ncid,'station',dummy,nstat
ncdf_varget,ncid,'country',country
ncdf_varget,ncid,'tree',tree
ncdf_varget,ncid,'latitude',ylat
ncdf_varget,ncid,'longitude',xlon
ncdf_close,ncid
;
multiplot,1
if !d.name eq 'X' then begin
  window,xsize=700,ysize=700
endif else begin
  device,xsize=17,xoffset=2
  device,ysize=22,yoffset=4
endelse
;
map=def_map(/npolar)  &  map.limit(0)=30.
coast=def_coast(/off)
labels=def_labels(/off)
;
inter_boxfd,/nodata,map=map,coast=coast,labels=labels
map_continents,/coasts,color=!p.color,mlinethick=0.5
map_continents,/countries,color=!p.color,mlinestyle=1,mlinethick=0.5
;
cpl_usersym,/circle,/fill
plots,xlon,ylat,psym=8,color=!p.background
cpl_usersym,/circle
plots,xlon,ylat,psym=8,thick=2
;
plot,[0,1],/noerase,xstyle=4,ystyle=4,/nodata,$
  xmargin=[2,2],ymargin=[2,2],$
  title='Location of density chronologies'
;
pause,/force
;
inter_boxfd,/nodata,map=map,coast=coast,labels=labels
map_continents,/coasts,color=!p.color,mlinethick=0.5
map_continents,/countries,color=!p.color,mlinestyle=1,mlinethick=0.5
;
xyouts,xlon,ylat,string(tree(0:1,*)),alignment=0.5,charthick=0.5,$
  charsize=0.6*!p.charsize
;
plot,[0,1],/noerase,xstyle=4,ystyle=4,/nodata,$
  xmargin=[2,2],ymargin=[2,2],$
  title='Genus??? of density chronologies'
;
pause,/force
;
inter_boxfd,/nodata,map=map,coast=coast,labels=labels
map_continents,/coasts,color=!p.color,mlinethick=0.5
map_continents,/countries,color=!p.color,mlinestyle=1,mlinethick=0.5
;
xyouts,xlon,ylat,strcompress(string(country),/remove_all),$
  alignment=0.5,charthick=0.5,$
  charsize=0.6*!p.charsize
;
plot,[0,1],/noerase,xstyle=4,ystyle=4,/nodata,$
  xmargin=[2,2],ymargin=[2,2],$
  title='Country of density chronologies'
;
end