Agreed, very nice.
I've heard the speed issue mentioned as well, although slow compared to what, I wonder? I don't see offhand any fundamental reason why it would be dramatically slower than, say, a standard FLIP-based fluid sim. It amounts to advecting Lagrangian particles combined with solving the elastoplasticity equations on a Cartesian grid, with a snow-specific constitutive model.