  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.
