There is an issue when calculating the vapor fluxes after iterating surface temperature (lines 475-500). Glacier ice is modified (line 497) when snow is on surface.