;
;input the 9 grids at various binning
; wedge data should be in dar
histMinCnt=5
stpAr=[.3,.6,.9,1.2]
nstp=n_elements(stpar)
nscans=9
sclp=1000. ; to mm
minR=0.
maxR=4.
minA=-4.
maxA=4.
bin=.2
start=1
for iscan=0,nscans-1 do begin
	print,"start scan:",iscan
	for istp=0,nstp-1 do begin
		stp=stpAr[istp]
		wedge_getgrid,wi,dar[iscan],stp,gdI
	    ii=where(gdi.cntar ge histMinCnt)
		herr=histogram(gdI.avgar[ii]*sclp,bin=bin,min=minA,max=maxA,loc=hlocE)
		hrms=histogram(gdI.rmsar[ii]*sclp,bin=bin,min=minR,max=maxR,loc=hlocR)
		if start then begin
			a={scan : 0l,$
               stp  : 0.,$
			   nx   : 0,$  ; dim covering grid
			   ny   : 0,$
			   avgAll:0.,$
               rmsAll:0.,$
               herr : herr,$
               hrms : hrms $
			}
			hAr=replicate(a,nstp,nscans)
		    start=0
		endif
		har[istp,iscan].scan=dar[iscan].scannum
		har[istp,iscan].stp = stp
		har[istp,iscan].herr= herr
		har[istp,iscan].hrms= hrms
		har[istp,iscan].nx=  gdI.nx
		har[istp,iscan].ny=  gdI.ny
		har[istp,iscan].avgall=gdi.avgall
		har[istp,iscan].rmsall=gdi.rmsall
	endfor
endfor
end