function liu_illingworth_cirrus_regression, dbz, temp_in, locale

; From Liu and Illingworth, 2000, JTech, vol39 p1130-1146.
; dbz is the water equivalent radar reflectivity factor with Z in mm^6/m^3
; temp is air temperature in Kelvins

temp=[200.,222.,228.,234.,240.,246.,252.,258.,264.,273.]
a_mid=[0.2093,0.3451,0.2136,0.1574,0.1619,0.1204,0.1044,0.09247,0.2001]
b_mid=[0.677,0.802,0.768,0.76,0.835,0.827,0.895,0.839,0.937]
a_tropics=[0.1854,0.1827,0.1716,0.1648,0.1440,0.1192,0.1215,0.1254,0.1235]
b_tropics=[0.658,0.677,0.705,0.723,0.757,0.774,0.819,0.767,0.797]

index=-1
for j=0,n_elements(temp)-2 do begin
	if temp_in ge temp[j] and temp_in lt temp[j+1] then index=j
endfor

;index=0

if index gt -1 and index lt n_elements(a_mid) then begin
if locale ne 'tropics' then begin
	a=a_mid[index] & b=b_mid[index]
endif else begin
	a=a_tropics[index] & b=b_tropics[index]
endelse



iwc=a*((10.^(dbz/10.))^b)

endif else begin

iwc=-9999.

endelse

return, iwc
end