Abstract
Two different types of numerical ocean circulation models are used in a classical idealized problem, the thermally induced circulation in an ocean basin bounded by two meridians to the east and west and by the equator and a line of constant latitude. A simple scaling theory exists for predicting poleward heat transport and the strength of meridional overturning as a function of vertical diffusivity and other external factors. However, previous studies have indicated conflicting results, and other scaling laws have been proposed. Experiments with two widely used types of numerical models, one based on depth coordinates and the other based on isopycnal layers, provide insight into the discrepancies of previous studies. In the numerical experiments vertical diffusivity is varied over a range of 200. The source of the difficulty in previous studies is in part traced to applying a fixed restoring coefficient at the upper boundary and considering the buoyancy forcing at the surface fixed irrespective o... Abstract Two different types of numerical ocean circulation models are used in a classical idealized problem, the thermally induced circulation in an ocean basin bounded by two meridians to the east and west and by the equator and a line of constant latitude. A simple scaling theory exists for predicting poleward heat transport and the strength of meridional overturning as a function of vertical diffusivity and other external factors. However, previous studies have indicated conflicting results, and other scaling laws have been proposed. Experiments with two widely used types of numerical models, one based on depth coordinates and the other based on isopycnal layers, provide insight into the discrepancies of previous studies. In the numerical experiments vertical diffusivity is varied over a range of 200. The source of the difficulty in previous studies is in part traced to applying a fixed restoring coefficient at the upper boundary and considering the buoyancy forcing at the surface fixed irrespective o...