Provost C., Bernier C. & Blayo L.E. (1994). A comparison of Two Numerical methods for Integrating a Quasi-geostrophic Multilayer Model of Ocean Circulations: Finite Difference Methods. Journal of computational physics, 110, 341–359.