An Unconditionally Stable Scheme for the Shallow Water Equations

Moshe Israeli Computer Science Department, Technion, Haifa, Israel

Search for other papers by Moshe Israeli in
Current site
Google Scholar
PubMed
Close
,
Naomi H. Naik Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York

Search for other papers by Naomi H. Naik in
Current site
Google Scholar
PubMed
Close
, and
Mark A. Cane Lamont-Doherty Earth Observatory, Columbia University, Palisades, New York

Search for other papers by Mark A. Cane in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

A finite-difference scheme for solving the linear shallow water equations in a bounded domain is described. Its time step is not restricted by a Courant–Friedrichs–Levy (CFL) condition. The scheme, known as Israeli–Naik–Cane (INC), is the offspring of semi-Lagrangian (SL) schemes and the Cane–Patton (CP) algorithm. In common with the latter it treats the shallow water equations implicitly in y and with attention to wave propagation in x. Unlike CP, it uses an SL-like approach to the zonal variations, which allows the scheme to apply to the full primitive equations. The great advantage, even in problems where quasigeostrophic dynamics are appropriate in the interior, is that the INC scheme accommodates complete boundary conditions.

Corresponding author address: Dr. Naomi Naik, Lamont-Doherty Earth Observatory, Columbia University, 106 D Oceanography, Palisades, NY 10964.

Email: naomi@lamont.ldeo.columbia.edu

Abstract

A finite-difference scheme for solving the linear shallow water equations in a bounded domain is described. Its time step is not restricted by a Courant–Friedrichs–Levy (CFL) condition. The scheme, known as Israeli–Naik–Cane (INC), is the offspring of semi-Lagrangian (SL) schemes and the Cane–Patton (CP) algorithm. In common with the latter it treats the shallow water equations implicitly in y and with attention to wave propagation in x. Unlike CP, it uses an SL-like approach to the zonal variations, which allows the scheme to apply to the full primitive equations. The great advantage, even in problems where quasigeostrophic dynamics are appropriate in the interior, is that the INC scheme accommodates complete boundary conditions.

Corresponding author address: Dr. Naomi Naik, Lamont-Doherty Earth Observatory, Columbia University, 106 D Oceanography, Palisades, NY 10964.

Email: naomi@lamont.ldeo.columbia.edu

Save
  • Cane, M. A., and E. S. Sarachik, 1979: Forced baroclinic ocean motion. Part III: An enclosed ocean. J. Mar. Res.,37, 355–398.

  • ——, and R. J. Patton, 1984: A numerical model for low-frequency equatorial dynamics. J. Phys. Oceanogr.,14, 1853–1863.

  • Golub, G. H., and C. F. Van Loan, 1989: Matrix Computations. 2d ed. The Johns Hopkins University Press, 642 pp.

  • Roe, P., 1994: Linear bicharacteristic schemes without dissipation. ICASE Tech. Rep. 94-65, pp. [Available from Institute for Computer Applications in Science and Engineering, Mail Stop 32C, NASA LARC, Hampton, VA 23681-0001.].

  • Staniforth, A., and J. Côté, 1991: Semi-Lagrangian integration schemes for atmospheric models—A review. Mon. Wea. Rev.,119, 2206–2223.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 185 43 10
PDF Downloads 149 35 9