; ; 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