Explicit iterative computation of diffusive vapour field in the 3-D snow matrix: preliminary results for low flux metamorphism