A robust and computationally efficient model of a two-dimensional coastal polynya