Adaptive finite element methods for Boussinesq equations