; ; Computes lists of tree numbers in each region ; ncid=ncdf_open('tree_dens_nh.nc') ncdf_varget,ncid,'longitude',xlon ncdf_varget,ncid,'latitude',ylat ncdf_close,ncid ; ; Remove Nepal from them all too, individually ********************** ; ;keeplist=where(ylat ge 30.,maxtrees) ;xlon=xlon(keeplist) ;ylat=ylat(keeplist) ; ; Define arrays ; maxtrees=n_elements(xlon) regname=['ALL','NORTH_OF_60','SOUTH_OF_60','WNA','ENA','CER','NEUR','SEUR',\$ 'ARCTIC','EXTRA-ARCTIC'] nreg=n_elements(regname) ntree=intarr(nreg) treelist=fltarr(maxtrees,nreg) & treelist(*,*)=!values.f_nan ; ; 1) All (except Nepal) ; i=0 keeplist=where(ylat ge 30.,nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 2) North of 60 ; i=1 keeplist=where(ylat ge 60.,nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 3) South of 60 ; i=2 keeplist=where( (ylat lt 60.) and (ylat ge 30.),nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 4) WNA ; i=3 keeplist=where( ( (xlon lt -105.) and (ylat lt 55.) ) \$ or ( (xlon lt -107.) and (ylat ge 55.) ) , nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 5) ENA ; i=4 keeplist=where( ( (xlon lt -40.) and (xlon ge -105.) and (ylat lt 55.) ) \$ or ( (xlon lt -40.) and (xlon ge -107.) and (ylat ge 55.) ) \$ , nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 6) CER ; i=5 keeplist=where( (xlon ge 60.) and (ylat ge 30.),nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 7) NEUR ; i=6 keeplist=where( (xlon ge -15.) and (xlon lt 60.) and (ylat ge 54.8) , nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 8) SEUR ; i=7 keeplist=where( (xlon ge -15.) and (xlon lt 35.) and (ylat lt 54.8) , nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 9) ARCTIC ; i=8 keeplist=where(ylat ge 66.6,nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; ; 10) OUTSIDE THE ARCTIC ; i=9 keeplist=where(ylat lt 66.6,nkeep) ntree(i)=nkeep treelist(0:nkeep-1,i)=keeplist(*) ; save,filename='reglists.idlsave',ntree,treelist,nreg,regname ; end