;Calculate the N-dbz relationship as per Frisch, 2000.


;Make bins 1 dbz wide.
bins=findgen(51)-60
binz=10^((bins-180.)/10.)
ns=fltarr(n_elements(bins),10)
rs=ns
rd=ns
binn=fltarr(n_elements(bins))
bind=binn
binr=binn
jwh=where(n gt 0)

for i=0,n_elements(bins)-1 do begin
	iwh=where(dbz[jwh] gt bins[i]-.5 and dbz[jwh] le bins[i]+.5)
	if iwh[0] eq -1 then begin
		ns[i,*]=0
		rs[i,*]=0
		rd[i,*]=0
		binn[i]=0
	endif else begin
		ns[i,*]=findgen(10)/9*(max(n[jwh[iwh]])-min(n[jwh[iwh]]))+min(n[jwh[iwh]])
		rs[i,*]=(binz[i]/(ns[i,*]*1e+6))^(1/6.)*0.47*(1e+6)
		for j=0,9 do begin
			rd[i,j]=sqrt(total((re[jwh[iwh]]-rs[i,j])^2)/(n_elements(iwh)-1))
		endfor
		wh=where(rd[i,*] eq min(rd[i,*]),num)
		if num gt 1 then begin
			binn[i]=mean(ns[i,wh])
			bind[i]=mean(rd[i,wh])
			binr[i]=mean(rs[i,wh])
		endif else begin
			binn[i]=ns[i,wh[0]]
			bind[i]=rd[i,wh[0]]
			binr[i]=rs[i,wh[0]]
		endelse
	endelse
endfor
end