This companion paper deals with numerical particularities of partial differential equations underlying 1D snow models. In this first part we neglect mechanical settling and demonstrate that the nonlinear coupling between diffusive transport (heat and vapor), phase changes and ice mass conservation contains a wave instability that may be relevant for weak layer formation. Numerical requirements are discussed in view of the underlying homogenization scheme.