The adaptive EVP method for solving the sea ice momentum equation

Stability and convergence of the modified EVP implementation of the visco-plastic sea ice rheology by Bouillon et al., Ocean Modell., 2013, is analyzed on B- and C-grids. It is shown that the implementation on a B-grid is less restrictive with respect to stability requirements than on a C-grid. On C...

Full description

Bibliographic Details
Published in:Ocean Modelling
Main Authors: Kimmritz, Madlen, Danilov, Sergey, Losch, Martin
Format: Article in Journal/Newspaper
Language:unknown
Published: ELSEVIER SCI LTD 2016
Subjects:
Online Access:https://epic.awi.de/id/eprint/40430/
https://epic.awi.de/id/eprint/40430/1/adaptiveEVP_r2.pdf
http://www.sciencedirect.com/science/article/pii/S1463500316300038
https://hdl.handle.net/10013/epic.47507
https://hdl.handle.net/10013/epic.47507.d001
Description
Summary:Stability and convergence of the modified EVP implementation of the visco-plastic sea ice rheology by Bouillon et al., Ocean Modell., 2013, is analyzed on B- and C-grids. It is shown that the implementation on a B-grid is less restrictive with respect to stability requirements than on a C-grid. On C-grids convergence is sensitive to the discretization of the viscosities. We suggest to adaptively vary the parameters of pseudotime subcycling of the modified EVP scheme in time and space to satisfy local stability constraints. This new approach generally improves the convergence of the modified EVP scheme and hence its numerical efficiency. The performance of the new “adaptive EVP” approach is illustrated in a series of experiments with the sea ice component of the MIT general circulation model (MITgcm) that is formulated on a C-grid