Contents
USAGE:
[SA_final, CT_final, w_Ih_final] = ...
gsw_melting_ice_into_seawater(SA,CT,p,w_Ih,t_Ih)
DESCRIPTION:
Calculates the Absolute Salinity and Conservative Temperature that
results when a given mass of ice melts and is mixed into a known mass
of seawater (whose properties are (SA,CT,p)).
When the mass fraction w_Ih_final is calculated as being a positive
value, the seawater-ice mixture is at thermodynamic equlibrium.
This code returns w_Ih_final = 0 when the input bulk enthalpy, h_bulk,
is sufficiently large (i.e. sufficiently "warm") so that there is no ice
present in the final state. In this case the final state consists of
only seawater rather than being an equlibrium mixture of seawater and
ice which occurs when w_Ih_final is positive. Note that when
w_Ih_final = 0, the final seawater is not at the freezing temperature.
INPUT:
SA = Absolute Salinity of seawater [ g/kg ]
CT = Conservative Temperature of seawater (ITS-90) [ deg C ]
p = sea pressure [ dbar ]
( i.e. absolute pressure - 10.1325 dbar )
w_Ih = mass fraction of ice, that is the mass of ice divided by the
sum of the masses of ice and seawater. That is, the mass of
ice divided by the mass of the final mixed fluid.
w_ice must be between 0 and 1. [ unitless ]
t_Ih = the in-situ temperature of the ice (ITS-90) [ deg C ]
SA, CT, w_Ih, and t_Ih must have the same dimensions.
p may have dimensions 1x1 or Mx1 or 1xN or MxN, where SA, CT, w_Ih and
t_Ih are MxN.
OUTPUT:
SA_final = Absolute Salinity of the seawater in the final state,
whether or not any ice is present. [ g/kg ]
CT_final = Conservative Temperature of the seawater in the the final
state, whether or not any ice is present. [ deg C ]
w_Ih_final = mass fraction of ice in the final seawater-ice mixture.
If this ice mass fraction is positive, the system is at
thermodynamic equilibrium. If this ice mass fraction is
zero there is no ice in the final state which consists
only of seawater which is warmer than the freezing
temperature. [unitless]
EXAMPLE:
SA = [34.7118; 34.8915; 35.0256; 34.8472; 34.7366; 34.7324;]
CT = [ 4.7856; 2.4329; 1.8103; 1.2600; 0.6886; 0.4403;]
p = [ 10; 50; 125; 250; 600; 1000;]
w_Ih = [0.0560; 0.02513; 0.02159; 0.01210; 0.00943; 0.00751;]
t_Ih = [-4.7856; -4.4329; -3.8103; -4.2600; -3.8863; -3.4036;]
[SA_final, CT_final, w_Ih_final] = ...
gsw_melting_ice_into_seawater(SA,CT,p,w_Ih,t_Ih)
SA_final =
32.767939199999994
34.014676604999998
34.269397295999994
34.425548880000001
34.409033862000001
34.471559675999998
CT_final =
-0.298448911022612
0.215263001418312
-0.074341719211557
0.207796293045473
-0.123785388299875
-0.202531182809225
w_Ih_final =
0
0
0
0
0
0
AUTHOR:
Trevor McDougall & Paul Barker [ help@teos-10.org ]
VERSION NUMBER:
3.05 (16th February, 2015)
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.
McDougall, T.J., P.M. Barker, R. Feistel and B.K. Galton-Fenzi, 2014:
Melting of Ice and Sea Ice into Seawater and Frazil Ice Formation.
Journal of Physical Oceanography, 44, 1751-1775.
The software is available from http://www.TEOS-10.org
