We extend our study for the warm cloud model in [13] to the analysis of a more general cloud model including the ice microphysics in [28]. The moisture variables comprise water vapor, cloud condensates (cloud water, cloud ice), and cloud precipitations (rain, snow), with respective mass ratios qv, qc and qp. A typical assumption in [13] for the calculation of condensation rate is that the warm clouds are exactly at water saturation with no supersaturation in general. When the ice microphysics are included, the situation becomes more complicated. We have to consider both the saturation mixing ratio with respect to water (qvw) and the saturation with respect to ice (q vi) when the temperature T is below the freezing point Tw but above the threshold T i for homogeneous ice nucleation. A remedy, acceptable from the physical and mathematical viewpoints, is to define the overall saturation mixing ratio qvs as a convex combination of qvw and q vi. Under this setting, supersaturation can still be avoided and we have the constraint qv ≤ qvs with qvs depending itself on the state. Mathematically, we are led to a system of equations and inequations involving some quasi-variational inequalities for which we prove the global existence and regularity of solutions.