Dealing with the Bottom Intersection Problem in Isopycnic-Coordinate Models Using the Flooding Algorithm


Article Preview

In isopycnic-coordinate models, special algorithm is required to compute the horizontal pressure force on steep bottom slopes. The “extrapolation” algorithm which constructs an “underground” estimate of pressure force using extrapolation was used in previous studies to deal with the bottom intersection problem. However, extrapolation will definitely introduce some errors, and would not be the best way to deal with this problem. In this paper, the flooding algorithm which is traditionally used to simulate transient flooding of coastlines caused by tides or tsunamis is used in the isopycnic-coordinate model to deal with the bottom intersection problem. The two algorithms were compared in the test case with a two-dimensional isolated Gaussian topography. The results of comparison indicate that the flooding algorithm performs better in dealing with the bottom intersection problem especially with the situation where interfaces near the bottom topography are going upwards.



Advanced Materials Research (Volumes 588-589)

Edited by:

Lawrence Lim




C. B. Miao and X. Q. Lv, "Dealing with the Bottom Intersection Problem in Isopycnic-Coordinate Models Using the Flooding Algorithm", Advanced Materials Research, Vols. 588-589, pp. 2125-2130, 2012

Online since:

November 2012




[1] Bleck, R., Smith, L.T.: A Wind-Driven Isopycnic Coordinate Model of the North and Equatorial Atlantic Ocean 1. Model Development and Supporting Experiments. Journal of Geophysical Research 95: 3273-3285 (1990).


[2] Boris, J.P., Book, D.L.: Flux-Corrected Transport. Journal of Computational Physics 11: 38-69 (1973).

[3] Zalesak, S.T.: Fully multidimensional flux-corrected transport algorithms for fluids. Journal of Computational Physics 31(3): 335-362 (1979).


[4] Bleck, R., Boudra, D.: Wind-Driven Spin-Up in Eddy-Resolving Ocean Models Formulated in Isopycnic and Isobaric Coordinates. Journal of Geophysical Research 91: 7611-7621 (1986).


[5] Sun, S., Bleck, R., Chassignet, E.P.: Layer Outcropping in Numerical Models of Stratified Flows. Journal of Physical Oceanography 23(8): 1877-1884 (1993).


[6] Bogue, N.M., Huang, R.X., Bryan, K.: Verification Experiments with an Isopycnal Coordinate Ocean Model. Journal of Physical Oceanography 16(5): 985-990 (1986).

[7] Miao, C., Chen, H., Lv, X.: An isopycnic-coornidate internal tide model and its application to the South China Sea. Chinese Journal of Oceanology and Limnology: 1339-1356 (2011).


[8] Niwa, Y., Hibiya, T.: Numerical study of the spatial distribution of the M2 internal tide in the Pacific Ocean. Journal of Geophysical Research 106: 22441-22449 (2001).


[9] Hallberg, R.: Stable Split Time Stepping Schemes for Large-Scale Ocean Modeling. Journal of Computational Physics 135(1): 54-65 (1997).


[10] Bleck, R., et al.: Salinity-driven Thermocline Transients in a Wind- and Thermohaline-forced Isopycnic Coordinate Model of the North Atlantic. Journal of Physical Oceanography 22(12): 1486-1505 (1992).


[11] Lahaye, S., et al.: A numerical scheme for modeling tidal wetting and drying. Journal of Geophysical Research 116(C3): C03028 (2011).

[12] Kämpf, J. LongWaves in a Channel, in Ocean Modelling for Beginners: Using Open-Source Software, Springer: Berlin, 2009; 65-89.

[13] Fringer, O.B., Armfield, S.W., Street, R.L.: Reducing numerical diffusion in interfacial gravity wave simulations. International Journal for Numerical Methods in Fluids 49: 301-329 (2005).


[14] Thacker, W.C.: Some exact solutions to the nonlinear shallow-water wave equations. Journal of Fluid Mechanics 107: 499-508 (1981).


[15] Smith, R.B.: Linear Theory of Stratified Flow past an Isolated Mountain in Isosteric Coordinates. Journal of the Atmospheric Sciences 45(24): 3889-3896 (1988).


[16] Carter, G.S., Merrifield, M.A.: Open boundary conditions for regional tidal simulations. Ocean Modelling 18(3-4): 194-209 (2007).