Density varies only by about a few percent in the ocean.
Hence, at the top of the domain, the azimuthal velocity is . A no-slip condition is used for the tangential velocity while the normal velocity is zero on the top, side and bottom walls. The equations for the resolved field are obtained by filtering Eqns. The filtered equations are transformed into curvilinear coordinates while the Cartesian velocity components are retained. Hence, a non-staggered grid is an attractive alternative and is used here.