pro condensed_water_parameterization

common cond_input, lwp,ze,temp,vdop,vap,z

common cond_output, iwc, lwc, frac
common cond_output2, re_liq, re_ice

ccm3_liquid_profile=fltarr(n_elements(z)) & liquid_wt=fltarr(n_elements(z))
lwc=fltarr(n_elements(z)) & iwc=fltarr(n_elements(z))
re_liq=fltarr(n_elements(z)) & re_ice=fltarr(n_elements(z))

for j=0,n_elements(z)-1 do begin
    if ze[j] gt -9999. then begin
ccm3_ice_frac=ccm3_cloud_iceliquid_fraction_param_function(temp[j])
;print, z[j], ' is height going into liquid param'
ccm3_liquid_profile[j]=(1.-ccm3_ice_frac)*ccm3_liquid_water_param_function(z[j], vap)
    endif
endfor

; now convert the liquid profile to a weighting function
    if lwp gt 0.0 and lwp lt 1.e4 then begin
for j=0,n_elements(z)-1 do begin
    if ze[j] gt -9999. then begin
if total(ccm3_liquid_profile) gt 0.0 then begin
     liquid_wt[j]=ccm3_liquid_profile[j]/total(ccm3_liquid_profile)
endif else begin
    liquid_wt[j]=0.0
endelse
lwc[j]=(liquid_wt[j]*(lwp))/90. ; assume 90m range gates

if 10.*alog10(ze[j]) lt 0.0 then begin
    re_liq[j]=19.5*exp(0.0384*(10.*alog10(ze[j])))
endif else begin
    re_liq[j]=8.4
endelse

;re_liq[j]=5.0

;print, liquid_wt[j],lwc[j],re_liq[j],((3./2.)*lwc[j]/re_liq[j])*90.,' is the liquid wt, lwc, and re'


    endif
endfor
    endif

; now determine the ice water content from the regression data.  Assume all water warmer than -15c and the regression below applies
; to all ice at colder temps

bz=-2.0 & ba=1.7
;IDL> print, r1   this is based on most data from 2000 iop with temp not in cc
;     0.509637     0.640465    0.0502367
;IDL> print, rmul1
 ;    0.660524
;IDL> print, mlr_coeff1
; 1.95209e-008
; 1.64473e-007
;-4.37153e-010
;IDL> print, mlr_const1
; 1.07064e-007

; 9/17/2004
;IDL> print, mlr_coeff1
; 8.76098e-008
; 2.67518e-010
;IDL> print, mlr_const1
;-1.52032e-008

coef=-1.52032e-008 & coef=[coef,8.76098e-008] & coef=[coef,2.67518e-010]; & coef=[coef,-4.37153e-010]
bz=-2.4 & ba=1.0 & bm=1.8 & bi=0.6

for j=0,n_elements(z)-1 do begin
    if ze[j] gt -9999. and temp[j] lt 273. then begin

	if temp[j] gt 268. and 10.*alog10(ze[j]) gt 10. then begin
		zzz=10.^(10./10.)
	endif else begin
		zzz=ze[j]
	endelse

    if vdop[j] lt 0. then vdop[j]=0.5
    ;z_temp=((ze[j]))^0.6
    ;v_temp=vdop[j]*((zzz))^(ba/(6.+bz))
    v_temp=vdop[j]*(((zzz)))^((ba/((6.+bz))))
    ;t_temp=611.*exp(((2.5e6)/461.5)*((1./273.)-(1./temp[j])))
    ;t_temp=temp[j]
    t_temp=611.*exp(((2.5e6)/461.5)*((1./273.)-(1./temp[j])))

    iwc[j]=coef[0]+(coef[1]*v_temp)+(coef[2]*t_temp)  ; iwc in g/cm3, z in mm6/m3, vdop in cm/s, temp in K
    iwc[j]=iwc[j]*1.e6 ; converts to g/m3
    ;re_ice[j]=15.
    re_ice[j]=ccm3_cloud_ice_re_param_function(exp(-z[j]/5500.))
    ;print, 're ice is ',re_ice[j],exp(-z[j]/5500.)

    endif else begin
        iwc[j]=0.0
    re_ice[j]=0.

    endelse
endfor



return
end