We propose new crystal melting models whose statistical partition functions are A-model amplitudes of the resolved conifold, with or without D-branes. The crystal is bounded by two walls whose distance corresponds to the Kahler modulus of the geometry. These statistical models are derived, via geometric transition, from Chern-Simons theory on the three-sphere in a novel matrix model formulation.