In order to solve the coupled system of energy balances and stochastic ice nucleation equations described above, we have used an adapted version of the algorithm presented earlier. Importantly, we note that repetitive simulations of the process are required to capture the variability resulting from the stochasticity of ice nucleation. The individual simulations are independent of each other, thus enabling parallelization of the computations.