Notation for equations

Symbols for variables

\(z\) refers to geopotential elevation (or height), increasing upward and with \(z=0\) defining the resting ocean surface. Much of the ocean has \(z < 0\).

\(x\) and \(y\) are the Cartesian horizontal coordinates. MOM6 uses generalized orthogonal curvilinear horizontal coordinates. However, the equations are simpler to write using Cartesian coordinates, and it is very straightforward to generalize the horizontal coordinates using the methods in Chapters 20 and 21 of [23].

\(\lambda\) and \(\phi\) are the geographic coordinates on a sphere (longitude and latitude, respectively).

Horizontal components of velocity are indicated by \(u\) and \(v\) and vertical component by \(w\).

\(p\) is the hydrostatic pressure.

\(\Phi\) is the geopotential. In the absence of tides, the geopotential is given by \(\Phi = g z,\) whereas more general expressions hold when including astronomical tide forcing.

The thermodynamic state variables can be salinity, \(S\), and potential temperature, \(\theta\). Alternatively, one can choose the Conservative Temperature if using the TEOS10 equation of state from [73].

\(\rho\) is the in-situ density computed as a function \(\rho(S,\theta,p)\) for non-Boussinesq ocean or \(\rho(S,\theta,p=-g \, \rho_o \, z)\) for Boussinesq ocean. See Young (2010) [71] or Section 2.4 of Vallis (2017) [65] for reasoning behind the simplified pressure used in the Boussinesq equation of state.

Vector notation

The three-dimensional velocity vector is denoted \(\mathbf{v}\) and it is decomposed into its horizontal and vertical components according to

\[\begin{align} \mathbf{v} = \mathbf{u} + \hat{\mathbf{z}} \, w = \hat{\mathbf{x}} \, u + \hat{\mathbf{y}} \, v + \hat{\mathbf{z}} \, w, \end{align}\]

where \(\hat{\mathbf{z}}\) is the unit vector pointed in the upward vertical direction and \(\mathbf{u} = (u, v, 0)\) is the horizontal component of velocity normal to the vertical.

The three-dimensional gradient operator is denoted \(\nabla\), and it is decomposed into its horizontal and vertical components according to

\[\begin{align} \nabla = \nabla_z + \hat{\mathbf{z}} \, \partial_z = \hat{\mathbf{x}} \, \partial_x + \hat{\mathbf{y}} \, \partial_y + \hat{\mathbf{z}} \, \partial_z. \end{align}\]