pro minnis_cloud_type, Cloud_Amount, column_type

; send in the 4 minnis types (low, middle, high, all) for a given point

column_type=' '


high_frac=cloud_amount[2]/100.
mid_frac=cloud_amount[1]/100.
low_frac=cloud_amount[0]/100.

if low_frac lt 0.1 and mid_frac lt 0.1 and high_frac ge 0.8 then column_type='high_clds'
if low_frac lt 0.1 and high_frac lt 0.1 and mid_frac ge 0.8 then column_type='mid_clds'
if mid_frac lt 0.1 and high_frac lt 0.1 and low_frac ge 0.8  then column_type='low_clds'
if mid_frac lt 0.1 and high_frac gt 0.8 and low_frac ge 0.8 then column_type='high_low'

if low_frac lt 0.1 and mid_frac lt 0.1 and high_frac gt 0.5 and high_frac lt 0.8 then column_type='broken_high_clds'
if low_frac lt 0.1 and high_frac lt 0.1 and mid_frac gt 0.5 and mid_frac lt 0.8 then column_type='broken_mid_clds'
if mid_frac lt 0.1 and high_frac lt 0.1 and low_frac gt 0.5 and low_frac lt 0.8 then column_type='broken_low_clds'

if low_frac eq 0.0 and mid_frac eq 0.0 and high_frac gt 0.1 and high_frac lt 0.5 then column_type='scattered_high_clds'
if low_frac eq 0.0 and high_frac eq 0.0 and mid_frac gt 0.1 and mid_frac lt 0.5 then column_type='scattered_mid_clds'
if mid_frac lt 0.0 and high_frac lt 0.0 and low_frac gt 0.1 and low_frac lt 0.5 then column_type='scattered_low_clds'

frac=max([high_frac,mid_frac])
frac=max([frac,low_frac])

if column_type eq ' ' and frac gt 0.8 then column_type='overcast_clds'
if column_type eq ' ' and frac gt 0.4 and frac lt 0.8 then column_type='broken_clds'
if column_type eq ' ' and frac gt 0.0 and frac lt 0.4 then column_type='scattered_clds'
if column_type eq ' ' and high_frac+mid_frac+low_frac eq 0.0 then column_type='clear'

;print, 'column cloud type is ',column_type

return
end