gsw_enthalpy_first_derivatives

first derivatives of enthalpy (48-term equation)

Contents

USAGE:

[h_SA, h_CT, h_P] = gsw_enthalpy_first_derivatives(SA,CT,p)

DESCRIPTION:

Calculates the following three derivatives of specific enthalpy (h),
of seawater using the computationally-efficient 48-term expression for 
density in terms of SA, CT and p (IOC et al., 2010).  
 (1) h_SA, the derivative with respect to Absolute Salinity at
     constant CT and p, and
 (2) h_CT, derivative with respect to CT at constant SA and p.
 (3) h_P, derivative with respect to pressure (in Pa) at constant
     SA and CT.
Note that the 48-term equation has been fitted in a restricted range of 
parameter space, and is most accurate inside the "oceanographic funnel" 
described in IOC et al. (2010).  The GSW library function 
"gsw_infunnel(SA,CT,p)" is avaialble to be used if one wants to test if 
some of one's data lies outside this "funnel".
TEOS-10 
Click for a more detailed description of the first
derivatives of specific enthalpy.

INPUT:

SA  =   Absolute Salinity                                       [ g/kg ]
CT  =   Conservative Temperature                               [ deg C ]
p   =   sea pressure                                            [ dbar ]
        (i.e. absolute pressure - 10.1325 dbar)
SA & CT need to have the same dimensions.
p may have dimensions 1x1 or Mx1 or 1xN or MxN, where SA & CT are MxN.

OUTPUT:

h_SA   =  The first derivative of specific enthalpy with respect to
          Absolute Salinity at constant CT and p.
                                          [ J/(kg (g/kg))]  i.e. [ J/g ]
h_CT   =  The first derivative of specific enthalpy with respect to
          CT at constant SA and p.                          [ J/(kg K) ]
h_P    =  The first partial derivative of specific enthalpy with
          respect to pressure (in Pa) at fixed SA and CT.  Note that
          h_P is specific volume (1/rho).                     [ m^3/kg ]

EXAMPLE:

SA = [34.7118; 34.8915; 35.0256; 34.8472; 34.7366; 34.7324;]
CT =  [28.7856; 28.4329; 22.8103; 10.2600;  6.8863;  4.4036;]
p =  [     10;      50;     125;     250;     600;    1000;]
[h_SA, h_CT] = gsw_enthalpy_first_derivatives(SA,CT,p)
h_SA =
  -0.070229561740134
  -0.351186407715501
  -0.887088218758157
  -1.829825716345477
  -4.424057561846865
  -7.406091544182000
h_CT =
1.0e+003 *
   3.991899707135326
   3.992025647578570
   3.992210374218789
   3.992284037753561
   3.992685015634706
   3.993013551060776

AUTHOR:

Trevor McDougall.                      [ help@teos-10.org ]

VERSION NUMBER:

3.04 (10th December, 2013)

REFERENCES:

IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of
 seawater - 2010: Calculation and use of thermodynamic properties.
 Intergovernmental Oceanographic Commission, Manuals and Guides No. 56,
 UNESCO (English), 196 pp.  Available from the TEOS-10 web site.
  See Eqns. (A.11.18), (A.11.15) and (A.11.12) of this TEOS-10 Manual.
This software is available from http://www.TEOS-10.org