Isobaric heat capacity
gsw_cp_t_exact(SA, t, p)
Absolute Salinity [ g/kg ]
in-situ temperature (ITS-90) [ degC ]
sea pressure [dbar], i.e. absolute pressure [dbar] minus 10.1325 dbar
heat capacity [ J/(kg*K) ]
The present R function works with a wrapper to a C function contained within the GSW-C system (Version 3.05-4 dated 2017-08-07, available at https://github.com/TEOS-10/GSW-C, as git commit '5b4d959e54031f9e972f3e863f63e67fa4f5bfec'), which stems from the GSW-Fortran system (https://github.com/TEOS-10/GSW-Fortran) which in turn stems from the GSW-Matlab system (https://github.com/TEOS-10/GSW-Matlab). Consult http://www.teos-10.org to learn more about these software systems, their authorships, and the science behind it all.
# NOT RUN {
SA <- c(34.7118, 34.8915, 35.0256, 34.8472, 34.7366, 34.7324)
t <- c(28.7856, 28.4329, 22.8103, 10.2600, 6.8863, 4.4036)
p <- c( 10, 50, 125, 250, 600, 1000)
cp_t_exact <- gsw_cp_t_exact(SA, t, p)
expect_equal(cp_t_exact/1e3, c(4.002888003958537, 4.000980283927373, 3.995546468894633,
3.985076769021370, 3.973593843482723, 3.960184084786622))
# }
Run the code above in your browser using DataLab