tm=[21.47111,21.57694,21.6925,21.85639,21.97139,22.11222]
ht=[5.4,6.7,7.127,7.121,6.45,5.59]

htind=intarr(6)
tmind=intarr(6)
mdopp=fltarr(6)
mdbz=fltarr(6)
mze=fltarr(6)
for i=0,5 do begin
	htind[i]=closest(height,ht[i],0)
	tmind[i]=closest(time,tm[i],0)
	wh=where(reflectivity[tmind[i]-10:tmind[i]+10,htind[i]] gt -900,nm)
	ind=tmind[i]-10+wh
	if nm gt 0 then begin
		mdbz[i]=10.0*alog10(mean(10.0^(reflectivity[ind,htind[i]]/10.0)))
		mze[i]=mean(10.0^(reflectivity[ind,htind[i]]/10.0))
		mdopp[i]=mean(velocity[ind,htind[i]])
	endif
endfor

;Calculate the atmospheric density correction factor for velocities
;----------------------------------------------
vcoeff=(exp(-height[htind]/8.5))^0.35
;vcoeff=1
dmed=0.84e-3*(mdopp*100.0*vcoeff)^3 - 6.60e-2*(mdopp*100.0*vcoeff)^2 + 6.10*(mdopp*100.0*vcoeff)-12.7
mdmnd=dmed/3.54
g3=((8.0e-5)*dmed^(-1.1))<(8.0e-5/(45.0^1.1))   ;Sets upper limit on G at 45 micron value
miwcd=mze/(dmed^(3.0)*g3)
x3=(0.00046*dmed^(-1.8))<(0.00046/(40.0^1.8))	;Sets upper limit on X at 40 micron value
mextd=mze/(x3*dmed^(4.0))

miwcs=0.0588*mze^0.63
mdmns=40.5*(0.0588)^(-0.5263)*mze^((1.0-(0.63))*0.5263)
mexts=mze/(x3*(mdmns*3.54)^(4.0))
end