Within M\+O\+M6, there is a wrapper for the equation of state, so that all calls look the same from the rest of the model. The equation of state code has to calculate not just in situ density, but also the compressibility and various derivatives of the density. There is also code for computing specific volume and the freezing temperature.\hypertarget{Equation_of_State_Linear_EOS}{}\section{Linear Equation of State}\label{Equation_of_State_Linear_EOS}
Compute the required quantities with uniform values for $\alpha = \frac{\partial \rho}{\partial T}$ and $\beta = \frac{\partial \rho}{\partial S}$, (D\+R\+H\+O\+\_\+\+DT, D\+R\+H\+O\+\_\+\+DS in M\+O\+M\+\_\+input, also uses R\+H\+O\+\_\+\+T0\+\_\+\+S0).\hypertarget{Equation_of_State_Wright_EOS}{}\section{Wright Equation of State}\label{Equation_of_State_Wright_EOS}
Compute the required quantities using the equation of state from \cite{wright1997}. This equation of state is in the form\+: \[ \alpha(s, \theta, p) = A(s, \theta) + \frac{\lambda(s, \theta)}{P(s, \theta) + p} \] where $A, \lambda$ and $P$ are functions only of $s$ and $\theta$ and $\alpha = 1/ \rho$ is the specific volume. This form is useful for the pressure gradient computation as discussed in \hyperlink{Discrete_PG_section_PG}{Pressure Gradient Term}.\hypertarget{Equation_of_State_NEMO_EOS}{}\section{N\+E\+M\+O Equation of State}\label{Equation_of_State_NEMO_EOS}
Compute the required quantities using the equation of state from \cite{roquet2015}.\hypertarget{Equation_of_State_UNESCO_EOS}{}\section{U\+N\+E\+S\+C\+O Equation of State}\label{Equation_of_State_UNESCO_EOS}
Compute the required quantities using the equation of state from \cite{jackett1995}.\hypertarget{Equation_of_State_TEOS-10_EOS}{}\section{T\+E\+O\+S-\/10 Equation of State}\label{Equation_of_State_TEOS-10_EOS}
Compute the required quantities using the equation of state from \href{http://www.teos-10.org/}{\tt T\+E\+O\+S-\/10}.\hypertarget{Equation_of_State_TFREEZE}{}\section{Freezing Temperature of Sea Water}\label{Equation_of_State_TFREEZE}
There are three choices for computing the freezing point of sea water\+:

\begin{DoxyItemize}
\item Linear The freezing temperature is a linear function of the salinity and pressure\+: \[ T_{Fr} = (T_{Fr0} + a\,S) + b\,P \] where $T_{Fr0},a,b$ are contants which can be set in M\+O\+M\+\_\+input (T\+F\+R\+E\+E\+Z\+E\+\_\+\+S0\+\_\+\+P0, D\+T\+F\+R\+E\+E\+Z\+E\+\_\+\+DS, D\+T\+F\+R\+E\+E\+Z\+E\+\_\+\+DP).\end{DoxyItemize}
\begin{DoxyItemize}
\item Millero The \cite{millero1978} equation is used, but modified so that it is a function of potential temperature rather than {\itshape in situ} temperature\+: \[ T_{Fr} = S(a + (b \sqrt{\max(S,0.0)} + c\, S)) + d\,P \] where $a,b, c, d$ are fixed contants.\end{DoxyItemize}
\begin{DoxyItemize}
\item T\+E\+O\+S-\/10 The T\+E\+O\+S-\/10 package is used to compute the freezing conservative temperature \mbox{[}degC\mbox{]} from absolute salinity \mbox{[}g/kg\mbox{]}, and pressure \mbox{[}Pa\mbox{]}. This one must be used if you are using the N\+E\+MO or T\+E\+O\+S-\/10 equation of state. \end{DoxyItemize}
