; ; Computes lists of grid-boxes in each region ; ;---------------------------------------------------- ; ; Get chronology locations ; ncid=ncdf_open('tree_dens_nh.nc') ncdf_varget,ncid,'longitude',xlon ncdf_varget,ncid,'latitude',ylat ncdf_close,ncid ; ; Get region info ; restore,'reglists.idlsave' ; ; Get IPCC instrumental dataset grid ; g=def_grid(/ipcc) ; ; For each region, locate the grid-boxes with data in ; boxlists=fltarr(g.nx,g.ny,nreg) boxnumbs=fltarr(g.nx,g.ny,nreg) for ireg = 0 , nreg-1 do begin statx=xlon(treelist(0:ntree(ireg)-1,ireg)) staty=ylat(treelist(0:ntree(ireg)-1,ireg)) statval=fltarr(ntree(ireg))+1. fdout=gridit(g.nx,g.ny,g.x,g.y,statx,staty,statval,nstat=nstat) boxlists(*,*,ireg)=fdout(*,*) boxnumbs(*,*,ireg)=nstat(*,*) endfor ; save,filename='regboxes.idlsave',nreg,boxlists,boxnumbs,regname ; end