; loadct,39 multi_plot,nrow=2,ncol=1,layout='large' if !d.name eq 'X' then begin window,ysize=850 fac=1. endif else begin fac=2. endelse def_1color,10,color='vvlgreen' ; map=def_map(/npolar) & map.limit(0)=30. labels=def_labels(/off) & labels.gridon=0 & labels.dlon=30. labels.dlat=10. coast=def_coast(/get_device) coast.double=0 coast.fill=1 coast.fillcolor=10 coast.thick=[0.5,0.5] ; inter_boxfd,/nodata,map=map,$ labels=labels,coast=coast ; xcir=findgen(361)-180. ycir=replicate(30.,361) map_plots,xcir,ycir ; x1=[180.,131.,75.,38.,-40.,-99.] y1=[30.,30.,30.,30.,30.,30.] x2=[180.,131.,75.,38.,-40.,-99.] y2=[90.,90.,90.,90.,90.,90.] n=n_elements(x1) for i = 0 , n-1 do begin map_plots,[x1(i),x2(i)],[y1(i),y2(i)],thick=5 endfor ; xs=findgen(79)-40. ys=replicate(50.,79) map_plots,xs,ys,thick=5 ; xs=findgen(82)-180. ys=replicate(56.,82) map_plots,xs,ys,thick=5 ; ; Plots maps of chronology location ; ncid=ncdf_open('tree_dens_nh.nc') ncdf_diminq,ncid,'station',dummy,nstat ncdf_varget,ncid,'country',country ncdf_varget,ncid,'tree',tree & tree=string(tree) ncdf_varget,ncid,'statid',statid & statid=string(statid) ncdf_varget,ncid,'statabbr',statabbr & statabbr=string(statabbr) ncdf_varget,ncid,'statloc',statloc & statloc=string((statloc <126B) >32B) ; get rid of control codes ncdf_varget,ncid,'yrstart',yrst ncdf_varget,ncid,'yrend',yren ncdf_varget,ncid,'tree',tree & tree=string(tree) ncdf_varget,ncid,'latitude',ylat ncdf_varget,ncid,'longitude',xlon ncdf_close,ncid ; restore,filename='reglists.idlsave' iall=where(regname eq 'ALL') i=iall(0) ; ; Output table of sites used ; ns=ntree[i] sl=treelist[0:ns-1,i] openw,1,'b98sens_site.txt' printf,1,'Table of tree-ring sites used in' printf,1 printf,1,'Briffa KR et al. (1998) Reduced sensitivity of recent tree growth to temperature at high northern latitudes. Nature 391, 678-682.' printf,1 printf,1,'Number of sites: '+txt(ns) printf,1 printf,1,' ID Name Abbr Genus Longitude Latitude Start End' for is = 0 , ns-1 do printf,1,strtrim(statid[sl[is]],2),strtrim(statloc[sl[is]],2),$ strtrim(statabbr[sl[is]],2),strtrim(tree[sl[is]],2),xlon[sl[is]],ylat[sl[is]],$ yrst[sl[is]],yren[sl[is]],$ format='(A6,A26,A9,A6,2F10.2,2I6)' close,1 ; x=xlon(treelist(0:ntree(i)-1,i)) y=ylat(treelist(0:ntree(i)-1,i)) cpl_usersym,/circle,/fill plots,x,y,psym=8,color=!p.background,symsize=1.05,noclip=0 cpl_usersym,/circle plots,x,y,psym=8,thick=1,symsize=1.05,noclip=0 print,'Number of sites:',n_elements(x) ; regname(0:2)=['ESIB','CSIB','WSIB'] rnx=[170,95,50,-35,-33,-48,-128,-160] rny=[40,45,42,53,33,33,75,36] for ir = 0 , 7 do $ xyouts,rnx(ir),rny(ir),regname(ir) ; end