Search Results
You are looking at 1 - 2 of 2 items for
- Author or Editor: T. Ezer x
- Refine by Access: All Content x
Abstract
In a recent paper by Mellor et al., it was found that, in two-dimensional (x, z) applications with finite horizontal viscosity and zero diffusivity, the velocity error, associated with the evaluation of horizontal density or pressure gradients on a sigma coordinate grid, prognostically disappeared, leaving behind a small and physically insignificant distortion in the density field. The initial error is numerically consistent in that it decreases as the square of the grid increment size. In this paper, we label this error as a sigma error of the first kind.
In three-dimensional applications, the authors have encountered an error that did not disappear and that has not been understood by us or, apparently, others. This is a vorticity error that is labeled a sigma error of the second kind and is a subject of this paper. Although it does not prognostically disappear, it seems to be tolerably small. To evaluate these numerical errors, the authors have adopted the seamount problem initiated by Beckman and Haidvogel. It represents a stringent test case, as evidenced by their paper, wherein the model is initialized with horizontal isopycnals, zero velocity, and no forcing; then, any velocities that develop must be considered errors.
Two appendices are important adjuncts to the paper, the first providing theoretical confirmation and understanding of the numerical results, and the second delving into additional errors related to horizontal or isosigma diffusion. It is, however, shown that satisfactory numerical solutions are obtained with zero diffusivity.
Abstract
In a recent paper by Mellor et al., it was found that, in two-dimensional (x, z) applications with finite horizontal viscosity and zero diffusivity, the velocity error, associated with the evaluation of horizontal density or pressure gradients on a sigma coordinate grid, prognostically disappeared, leaving behind a small and physically insignificant distortion in the density field. The initial error is numerically consistent in that it decreases as the square of the grid increment size. In this paper, we label this error as a sigma error of the first kind.
In three-dimensional applications, the authors have encountered an error that did not disappear and that has not been understood by us or, apparently, others. This is a vorticity error that is labeled a sigma error of the second kind and is a subject of this paper. Although it does not prognostically disappear, it seems to be tolerably small. To evaluate these numerical errors, the authors have adopted the seamount problem initiated by Beckman and Haidvogel. It represents a stringent test case, as evidenced by their paper, wherein the model is initialized with horizontal isopycnals, zero velocity, and no forcing; then, any velocities that develop must be considered errors.
Two appendices are important adjuncts to the paper, the first providing theoretical confirmation and understanding of the numerical results, and the second delving into additional errors related to horizontal or isosigma diffusion. It is, however, shown that satisfactory numerical solutions are obtained with zero diffusivity.
Abstract
Much has been written of the error in computing the horizontal pressure gradient associated with sigma coordinates in ocean or atmospheric numerical models. There also exists the concept of “hydrostatic inconsistency” whereby, for a given horizontal resolution, increasing the vertical resolution may not be numerically convergent.
In this paper, it is shown that the differencing scheme cited here, though conventional, is not hydrostatically inconsistent; the sigma coordinate, pressure gradient error decreases with the square of the vertical and horizontal grid size. Furthermore, it is shown that the pressure gradient error is advectively eliminated after a long time integration. At the other extreme, it is shown that diagnostic calculations of the North Atlantic Ocean using rather coarse resolution, and where the temperature and salinity and the pressure gradient error are held constant, do not exhibit significant differences when compared to a calculation where horizontal pressure gradients are computed on z-level coordinates. Finally, a way of canceling the error ab initio is suggested.
Abstract
Much has been written of the error in computing the horizontal pressure gradient associated with sigma coordinates in ocean or atmospheric numerical models. There also exists the concept of “hydrostatic inconsistency” whereby, for a given horizontal resolution, increasing the vertical resolution may not be numerically convergent.
In this paper, it is shown that the differencing scheme cited here, though conventional, is not hydrostatically inconsistent; the sigma coordinate, pressure gradient error decreases with the square of the vertical and horizontal grid size. Furthermore, it is shown that the pressure gradient error is advectively eliminated after a long time integration. At the other extreme, it is shown that diagnostic calculations of the North Atlantic Ocean using rather coarse resolution, and where the temperature and salinity and the pressure gradient error are held constant, do not exhibit significant differences when compared to a calculation where horizontal pressure gradients are computed on z-level coordinates. Finally, a way of canceling the error ab initio is suggested.
