First Derivatives of Enthalpy wrt CT
First Derivatives of Enthalpy wrt CT
gsw_enthalpy_first_derivatives_CT_exact(SA, CT, p)
SA |
Absolute Salinity [ g/kg ] |
CT |
Conservative Temperature [ degC ] |
p |
sea pressure [dbar], i.e. absolute pressure [dbar] minus 10.1325 dbar |
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.
a list containing h_SA [ (J/kg)/(g/kg) ], the derivative
of enthalpy wrt Absolute Salinity, and h_CT [ (J/kg)/degC ],
the derivative of enthalpy wrt Conservative Temperature.
The HTML documentation suggests that this function returns 3 values, but there are only 2 returned values in the C code used here (and the matlab code on which that is based). Also, the d/dSA check values given the HTML are not reproduced by the present function. This was reported on Mar 18, 2017 as https://github.com/TEOS-10/GSW-Matlab/issues/7. See https://github.com/TEOS-10/GSW-R/issues/34
Other things related to enthalpy: gsw_CT_from_enthalpy,
gsw_dynamic_enthalpy,
gsw_enthalpy_CT_exact,
gsw_enthalpy_diff,
gsw_enthalpy_first_derivatives,
gsw_enthalpy_ice,
gsw_enthalpy_t_exact,
gsw_enthalpy,
gsw_frazil_properties_potential_poly,
gsw_frazil_properties_potential,
gsw_pot_enthalpy_from_pt_ice_poly,
gsw_pot_enthalpy_from_pt_ice,
gsw_pot_enthalpy_ice_freezing_poly,
gsw_pot_enthalpy_ice_freezing,
gsw_pt_from_pot_enthalpy_ice_poly,
gsw_pt_from_pot_enthalpy_ice,
gsw_specvol_first_derivatives_wrt_enthalpy,
gsw_specvol_first_derivatives
SA <- c(34.7118, 34.8915, 35.0256, 34.8472, 34.7366, 34.7324)
CT <- c(28.7856, 28.4329, 22.8103, 10.2600, 6.8863, 4.4036)
p <- c( 10, 50, 125, 250, 600, 1000)
d <- gsw_enthalpy_first_derivatives_CT_exact(SA, CT, p)
expect_equal(d$h_SA, c(-0.070224183838619, -0.351159869043798, -0.887036550157504,
-1.829626251448858, -4.423522691827955, -7.405211691293971))
expect_equal(d$h_CT/1e3, c(3.991899712269790, 3.992025674159605, 3.992210402650973,
3.992283991748418, 3.992685275917238, 3.993014370250710))Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.