The code works in 2D if the initial domain is not completely dry. With some wet cells, the model works without problems