pro pmc_iceprops,clda,mod_vol,mod_scacoef,iwc,ice_colden

; calculate CIPS ice water content
; clda is cloud albedo in G (10^-6 str^-1)
; vol is volume per ice particle cm^3 (mod_vol is assumed to come in m^3)
; scacoef is scattering coeficient in cm^2/str (mod_scacoef is assumed m^2/str
; iwc is ice water contnet in micro-g m^-2  ; since it's a common unit

vol=mod_vol*1e6
scacoef=mod_scacoef*1e4

albedo=clda*1e-6 ; albedo in str^-1

; get momennts of the size distribution

;gauss_mom,cldr,gauss_width,3.,mom3
;gauss_mom,cldr,gauss_width,0.,mom0

iceden=0.92 *1000. ; kg m^-2
iceden=iceden *1000./1e6 ; g cm^-2

ice_colden=albedo/scacoef

iwc = ice_colden * vol * iceden  * 1e10

return
end