; ; Plots location of MXD grid boxes with complete 1697-1976 data, ; differentiating those that were infilled. ; doinfill=1 ; 0=don't mark infilled boxes, 1=do, 2=do, with same colour trv=0 ; selects tree-ring-variable: 0=MXD 1=TRW case trv of 0: begin fnadd='mxd' end 1: begin fnadd='trw' end endcase titadd=strupcase(fnadd) ; restore,fnadd+'_infill.idlsave' ; shyr=[1950,1400,1583,1697] nshow=n_elements(shyr) ; titpt='('+['a','c','e','g']+') ' ; loadct,39 multi_plot,nrow=4,ncol=1,layout='large' if !d.name eq 'X' then begin window,ysize=750 !p.font=-1 endif else begin !p.font=0 device,/helvetica,/bold,font_size=14 endelse def_1color,cr,cg,cb,10,color='black' def_1color,cr,cg,cb,11,color='mgrey' def_1color,cr,cg,cb,12,color='mdgrey' ; map=def_map(/npolar) & map.limit(0)=25. map.xmargin=[0,0] & map.ymargin=[0,0] labels=def_labels(/off) coast=def_coast(/get_device) & coast.thick=[0.5,2] ;coast.color=12 ; for i = 0 , nshow-1 do begin ; iyr=where(mxdyear eq shyr(i)) fd1=reform(mxdfd(*,*,iyr)) fd2=reform(mxdfd2(*,*,iyr)) fd=fd1 fd(*,*)=!values.f_nan ; fd(where(finite(fd2)))=1. fd(where(finite(fd2)))=doinfill fd(where(finite(fd1)))=2. ; inter_boxfd,fd,g.x,g.y,levels=[0.5,1.5,2.5],c_colors=[11,10],$ labels=labels,map=map,coast=coast ; xyouts,-55,25,titpt[i]+string(shyr(i),format='(I4)') ; endfor ; end