Recently, I’m looking into the Meridional Heat Transport across the 35$^{\circ}$S in our high-resolution CCSM simulation following Dong et al. 2011. To calculate the heat transport, the following formula is used:

\begin{equation} \int \rho C_{p}v\theta dx dz \end{equation}

The tricky part is whether one should use temporal and spatial-varying values of $C_{p}$ (seawater heat capacity) and $\rho$ (seawater density). Since density and heat capacity are not among the standard outputs, I went on looking for some standard libraries to calculate them using the fields I have: temperature, pressure, and salinity. Soon, I found the Python version of **Gibbs SeaWater Oceanographic Package**.