gsw_melting_ice_into_seawater

the resulting SA and CT when ice is melted
into seawater

Contents

USAGE:

[SA_final, CT_final] = ...
       gsw_melting_ice_into_seawater(SA,CT,p,saturation_fraction,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)). 

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 )
saturation_fraction = the saturation fraction of dissolved air in 
             seawater. saturation_fraction must be between 0 and 1.          
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 and saturation_fraction 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 mixture of the melted ice 
             and the orignal seawater                           [ g/kg ]
CT_final  =  Conservative Temperature of the mixture of the melted 
             ice and the orignal seawater                      [ deg C ]

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;]
saturation_fraction = [1; 0.8; 0.6; 0.5; 0.4; 0;]
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] = ...
  gsw_melting_ice_into_seawater(SA,CT,p,saturation_fraction,w_Ih,t_Ih)
SA_final =
   32.767939199999994
  34.014676604999998
  34.269397295999994
  34.425548880000001
  34.409033862000001
  34.471559675999998
CT_final =
  -0.298448911022612
   0.215263001418311
  -0.074341719211558
   0.207796293045473
  -0.123785388299875
  -0.202531182809225

AUTHOR:

Trevor McDougall & Paul Barker       [ help@teos-10.org ]

VERSION NUMBER:

3.04 (3rd January, 2014)

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.
The software is available from http://www.TEOS-10.org