Frost boundaries are calculated as model outputs in a separate subroutine as isotherms of 0 oC. The somewhat less simplistic assumption of a linear heat change between adjacent layers, give these isotherms a strong dependence on the choice of layer thickness. Not more than two frost layers are allowed to occur simultaneously for output purposes.