;$Id$
;-----------------------------------------------------------------------------------------------------------
;
;-----------------------------------------------------------------------------------------------------------
pro get_prior_2d, x, isfc, lut, attrs, prior_prob_clear, prior_prob_ice, prior_prob_water

    ix = min([attrs.nbins_x -1,max([0,fix((x - attrs.x_min) / attrs.x_bin)])])

    probc = lut[isfc].clear[ix,*]
    probi = lut[isfc].ice[ix,*]
    probw = lut[isfc].water[ix,*]
    obs = lut[isfc].obs_prob[ix,*]
    idx = where(obs gt 0,cc)
    prior_prob_clear = total(probc[idx] * obs[idx])/total(obs[idx])
    prior_prob_water = total(probw[idx] * obs[idx])/total(obs[idx])
    prior_prob_ice = total(probi[idx] * obs[idx])/total(obs[idx])
end