; ; Plot locations of instrumental data boxes in each region and ; those with sufficient data and with MXD sites in them ; ; Prepare for plotting ; loadct,39 multi_plot,nrow=2,layout='large' if !d.name eq 'X' then begin window,xsize=650,ysize=950 endif else begin device,xsize=15,xoffset=3 device,ysize=24,yoffset=2 endelse def_1color,40,color='lgreen' def_1color,41,color='deepblue' def_1color,42,color='red' ; ; Define map projection ; map=def_map(/npolar) & map.limit(0)=25. coast=def_coast(/get_device) & coast.double=0 & coast.coasts=1 labels=def_labels(/off) ; ; Get regional data of instrumental boxes ; restore,filename='reg_mxdboxes.idlsave' ; nreg,regname,boxlists,boxnumbs,g,boxregs,boxprec ; ; Get regional data of MXD sites ; restore,filename='allmxd.idlsave' ; nchron,idno,idname,location,country,tree,yrstart,yrend,statlat,statlon,$ ; mxd,fraction,timey,nyr restore,filename='reg_mxdlists.idlsave' ; ntree,treelist,nreg,regname cpl_usersym,/circle ; for i = 0 , nreg-1 do begin ; pause ; fd1=boxlists(*,*,i) ; boxes with MXD sites ; fd2=boxprec(*,*,i) ; boxes within region fd3=fd2+finite(fd1) mlist=where( (finite(fd1) ne 0) and (finite(fd2) eq 0),nmonly) if nmonly gt 0 then fd3(mlist)=3. ; inter_boxfd,fd3,g.x,g.y,map=map,coast=coast,labels=labels,$ title=regname(i)+': PRECIPITATION & MXD DATA',$ levels=[0.5,1.5,2.5,3.5],c_colors=[40,41,42] n=ntree(i) x=statlon(treelist(0:n-1,i)) y=statlat(treelist(0:n-1,i)) plots,x,y,psym=8,thick=2,symsize=0.5,color=!p.background ; fd2=boxtemp(*,*,i) ; boxes within region fd3=fd2+finite(fd1) mlist=where( (finite(fd1) ne 0) and (finite(fd2) eq 0),nmonly) if nmonly gt 0 then fd3(mlist)=3. ; inter_boxfd,fd3,g.x,g.y,map=map,coast=coast,labels=labels,$ title=regname(i)+': TEMPERATURE & MXD DATA',$ levels=[0.5,1.5,2.5,3.5],c_colors=[40,41,42] plots,x,y,psym=8,thick=2,symsize=0.5,color=!p.background ; endfor ; end