On the Budyko–Sellers climate model with mushy region