next up previous
Next: 4. Seamount Test Case Up: Improvement of Estuarine and Previous: 2.3 Sixth-Order Compact Difference

3. Pressure Gradient Error in the tex2html_wrap_inline528-Coordinate System

Coastal ocean models usually use a terrain-following tex2html_wrap_inline528-coordinate system to handle the effects of bottom topography. Here the water column is divided into the same number of grid cells independence of depth . Let (tex2html_wrap_inline562 denote Cartesian coordinates and (tex2html_wrap_inline564 sigma coordinates. In most sigma coordinate ocean models the relationship between the two coordinate systems are:
equation339
where z and tex2html_wrap_inline528 increase vertically upward such that tex2html_wrap_inline570 at the surface and z=-H and tex2html_wrap_inline574 at the bottom. H=H(x,y) is the bottom topography.

A problem has long been recognized in computing the horizontal pressure gradient in the tex2html_wrap_inline528-coordinate system (e.g., Gary, 1973; Haney, 1991; and Mellor et al., 1994, McCalpin, 1994): the horizontal pressure gradient becomes a difference between two terms, which leads to a large truncation error at a steep topography.



Peter Chu
Thu Aug 24 17:01:43 PDT 2000