Monte Carlo simulation of ice models

We propose a number of Monte Carlo algorithms for the simulation of ice models and compare their efficiency. One of them, a cluster algorithm for the equivalent three colour model, appears to have a dynamic exponent close to zero, making it particularly useful for simulations of critical ice models. We have performed extensive simulations using our algorithms to determine a number of critical exponents for the square ice and F models.