Version 3.03

Added

isentropic compressibility (kappa) for the 48-term equation of state.
gsw_rho_first derivatves and gsw_rho_first derivatves_CT_exact.
gsw_alpha_on_beta and gsw_alpha_on_beta_CT_exact.
gsw_t_from_pt0
gsw_t_from_entropy
gsw_adiabatic_lapse_rate_from_CT.
gsw_deltaSA_from_rho_t_exact.
Remaned gsw_entropy_t_exact to gsw_entropy_from_t, and added ths inverse
function (gsw_t_from_entropy)

Renamed gsw_adiabatic_lapse_rate_t_exact to
gsw_adiabatic_lapse_rate_from_t

Toolbox/gsw_Abs_Pressure_from_p.m

@@ -19,7 +19,7 @@
 % AUTHOR: 
 %  Trevor McDougall and Paul Barker                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (13th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_C3515.m

@@ -17,7 +17,7 @@
 % AUTHOR: 
 %  Trevor McDougall and Paul Barker                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (13th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  Culkin and Smith, 1980:  Determination of the Concentration of Potassium  

Toolbox/gsw_CT_first_derivatives.m

@@ -34,7 +34,7 @@
 % AUTHOR: 
 %  Trevor McDougall and Paul Barker                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (11th March, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 
@@ -79,12 +79,9 @@
 %--------------------------------------------------------------------------
 
 cp0 = 3991.86795711963;             % from Eqn. 3.3.3 of IOC et al. (2010).
-n0 = 0; 
-n2 = 2;
-pr0 = zeros(size(SA));
 abs_pt = 273.15 + pt; 
 
-CT_pt = - (abs_pt.*gsw_gibbs(n0,n2,n0,SA,pt,pr0))./cp0;
+CT_pt = - (abs_pt.*gsw_gibbs_pt0_pt0(SA,pt))./cp0;
 
 %--------------------------------------------------------------------------
 

Toolbox/gsw_CT_freezing.m

@@ -31,7 +31,7 @@
 % AUTHOR: 
 %  Trevor McDougall, Paul Barker and Rainer Feistal    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (13th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_CT_from_entropy.m

@@ -22,7 +22,7 @@
 % AUTHOR:  
 %  Trevor McDougall and Paul Barker.                   [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (13th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of

Toolbox/gsw_CT_from_pt.m

@@ -22,7 +22,7 @@
 % AUTHOR: 
 %  David Jackett, Trevor McDougall and Paul Barker     [ [email protected] ]
 %  
-% VERSION NUMBER: 3.02 (13th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_CT_from_rho.m

@@ -11,11 +11,11 @@
 %  Calculates the Conservative Temperature of a seawater sample, for given
 %  values of its density, Absolute Salinity and sea pressure (in dbar), 
 %  using the computationally-efficient 48-term expression for density in 
-%  terms of SA, CT and p (McDougall et al., 2011)
+%  terms of SA, CT and p (IOC et al., 2010)
 %
 %  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 McDougall et al. (2011).  The GSW library function 
+%  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".  
 %
@@ -43,19 +43,14 @@
 % AUTHOR:
 %  Trevor McDougall & Paul Barker                      [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 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 http://www.TEOS-10.org
 %
-%  McDougall T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 
-%   computationally efficient 48-term expression for the density of 
-%   seawater in terms of Conservative Temperature, and related properties
-%   of seawater.  To be submitted to J. Atm. Ocean. Technol., xx, yyy-zzz.
-%
 %  The software is available from http://www.TEOS-10.org
 %
 %==========================================================================

Toolbox/gsw_CT_from_rho_exact.m

@@ -34,7 +34,7 @@
 % AUTHOR:
 %  Trevor McDougall & Paul Barker                      [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (13th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of

Toolbox/gsw_CT_from_t.m

@@ -25,7 +25,7 @@
 % AUTHOR: 
 %  David Jackett, Trevor McDougall and Paul Barker     [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (13th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_CT_maxdensity.m

@@ -12,12 +12,12 @@
 %  This function returns the Conservative temperature at which the density
 %  of seawater is a maximum, at given Absolute Salinity, SA, and sea 
 %  pressure, p (in dbar).  This function uses the computationally-efficient
-%  48-term expression for density in terms of SA, CT and p (McDougall et
-%  al., 2013).
+%  48-term expression for density in terms of SA, CT and p (IOC et
+%  al., 2010).
 %
 %  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 McDougall et al. (2013).  The GSW library function 
+%  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".  
 %
@@ -36,7 +36,7 @@
 % AUTHOR: 
 %  Trevor McDougall & Paul Barker                      [ [email protected] ]
 %
-% VERSION NUMBER: 3.01 (3rd April, 2011)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 
@@ -45,11 +45,6 @@
 %   UNESCO (English), 196 pp.  Available from http://www.TEOS-10.org
 %    See section 3.42 of this TEOS-10 Manual.  
 %
-%  McDougall T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 
-%   computationally efficient 48-term expression for the density of 
-%   seawater in terms of Conservative Temperature, and related properties
-%   of seawater.  To be submitted to J. Atm. Ocean. Technol., xx, yyy-zzz.
-%
 %  McDougall T.J. and S.J. Wotherspoon, 2012: A simple modification of 
 %   Newton’s method to achieve convergence of order "1 + sqrt(2)".
 %   Submitted to Applied Mathematics and Computation. 

Toolbox/gsw_CT_maxdensity_exact.m

@@ -28,7 +28,7 @@
 % AUTHOR: 
 %  Trevor McDougall & Paul Barker                      [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_CT_second_derivatives.m

@@ -44,7 +44,7 @@
 % AUTHOR: 
 %  Trevor McDougall and Paul Barker                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 
@@ -53,11 +53,6 @@
 %   UNESCO (English), 196 pp.  Available from http://www.TEOS-10.org.  
 %    See appendix A.12 of this TEOS-10 Manual.    
 %
-%  McDougall T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 
-%   computationally efficient 48-term expression for the density of 
-%   seawater in terms of Conservative Temperature, and related properties
-%   of seawater.  To be submitted to J. Atm. Ocean. Technol., xx, yyy-zzz.
-%
 %  This software is available from http://www.TEOS-10.org
 %
 %==========================================================================

Toolbox/gsw_C_from_SP.m

@@ -43,7 +43,7 @@
 % AUTHOR:
 %  Trevor McDougall, Paul Barker and Rich Pawlowicz    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (16th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  Hill, K.D., T.M. Dauphinee and D.J. Woods, 1986: The extension of the

Toolbox/gsw_IPV_vs_fNsquared_ratio.m

@@ -21,11 +21,11 @@
 %  IPV_vs_fNsquared_ratio is evaluated at the mid pressure between the 
 %  individual data points in the vertical.  This function uses the 
 %  computationally-efficient 48-term expression for density in terms of 
-%  SA, CT and p (McDougall et al., 2013). 
+%  SA, CT and p (IOC et al., 2010). 
 %
 %  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 McDougall et al. (2013).  The GSW library function 
+%  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".  
 %
@@ -53,7 +53,7 @@
 % AUTHOR:  
 %  Trevor McDougall and Paul Barker                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 
@@ -62,11 +62,6 @@
 %   UNESCO (English), 196 pp.  Available from http://www.TEOS-10.org
 %    See Eqn. (3.20.5) of this TEOS-10 Manual. 
 %
-%  McDougall T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 
-%   computationally efficient 48-term expression for the density of 
-%   seawater in terms of Conservative Temperature, and related properties
-%   of seawater.  To be submitted to J. Atm. Ocean. Technol., xx, yyy-zzz.
-%
 %   The software is available from http://www.TEOS-10.org
 %
 %==========================================================================

Toolbox/gsw_Nsquared.m

@@ -23,7 +23,7 @@
 %  efficient 48-term expression for density in terms of SA, CT and p.  The    
 %  48-term equation has been fitted in a restricted range of parameter 
 %  space, and is most accurate inside the "oceanographic funnel" described 
-%  in McDougall et al. (2013).  The GSW library function 
+%  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".  
 %
@@ -49,7 +49,7 @@
 % AUTHOR:  
 %  Trevor McDougall and Paul Barker                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  Griffies, S. M., 2004: Fundamentals of Ocean Climate Models. Princeton, 
@@ -61,11 +61,6 @@
 %   UNESCO (English), 196 pp.  Available from http://www.TEOS-10.org
 %    See section 3.10 and Eqn. (3.10.2) of this TEOS-10 Manual. 
 %
-%  McDougall T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 
-%   computationally efficient 48-term expression for the density of 
-%   seawater in terms of Conservative Temperature, and related properties
-%   of seawater.  To be submitted to J. Atm. Ocean. Technol., xx, yyy-zzz.
-%
 %   The software is available from http://www.TEOS-10.org
 %
 %==========================================================================

Toolbox/gsw_P0.m

@@ -15,7 +15,7 @@
 % AUTHOR: 
 %  Trevor McDougall and Paul Barker                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_R_from_SP.m

@@ -43,7 +43,7 @@
 % AUTHOR:
 %  Trevor McDougall, Paul Barker and Rich Pawlowicz    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  Culkin and Smith, 1980:  Determination of the Concentration of Potassium  

Toolbox/gsw_SA_CT_plot.m

@@ -14,11 +14,11 @@ function gsw_SA_CT_plot(SA,CT,p_ref,isopycs,title_string)
 %  point for p = 0 dbar assuming the seawater is completely saturated with
 %  dissolved air and user defined potential density contours.  This 
 %  function uses the computationally efficient 48-term expression for 
-%  density in terms of SA, CT and p (McDougall et al., 2011).  
+%  density in terms of SA, CT and p (McDougall et al., 2013).  
 %
 %  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 McDougall et al. (2011).  The GSW library function 
+%  described in McDougall et al. (2013).  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".  
 %
@@ -47,7 +47,7 @@ function gsw_SA_CT_plot(SA,CT,p_ref,isopycs,title_string)
 % MODIFIED:
 %  Paul Barker & Trevor McDougall
 %
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  McDougall T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 

Toolbox/gsw_SA_Sstar_from_SP.m

@@ -36,7 +36,7 @@
 % AUTHOR: 
 %  Paul Barker and Trevor McDougall                    [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (7th January, 2013)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_SA_from_SP.m

@@ -33,7 +33,7 @@
 % AUTHOR: 
 %  David Jackett, Trevor McDougall & Paul Barker       [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (7th January, 2013)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_SA_from_Sstar.m

@@ -31,7 +31,7 @@
 % AUTHOR: 
 %  David Jackett, Trevor McDougall and Paul Barker     [ [email protected] ]
 %
-% VERSION NUMBER: 3.02 (7th January, 2013)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 

Toolbox/gsw_SA_from_rho.m

@@ -10,11 +10,11 @@
 %  Calculates the Absolute Salinity of a seawater sample, for given values
 %  of its density, Conservative Temperature and sea pressure (in dbar). 
 %  This function uses the computationally-efficient 48-term expression for 
-%  density in terms of SA, CT and p (McDougall et al., 2013).
+%  density in terms of SA, CT and p (IOC et al., 2013).
 %
 %  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 McDougall et al. (2013).  The GSW library function 
+%  described in IOC et al. (2013).  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".  
 %
@@ -37,7 +37,7 @@
 % AUTHOR: 
 %  Trevor McDougall & Paul Barker                      [ [email protected] ]
 %      
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 
@@ -46,11 +46,6 @@
 %   UNESCO (English), 196 pp.  Available from http://www.TEOS-10.org
 %    See section 2.5 of this TEOS-10 Manual. 
 %
-%  McDougall, T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 
-%   computationally efficient 48-term expression for the density of 
-%   seawater in terms of Conservative Temperature, and related properties
-%   of seawater.  To be submitted to J. Atm. Ocean. Technol., xx, yyy-zzz.
-%
 %  Millero, F. J., R. Feistel, D. G. Wright, and T. J. McDougall, 2008: 
 %   The composition of Standard Seawater and the definition of the 
 %   Reference-Composition Salinity Scale, Deep-Sea Res. I, 55, 50-72. 

Toolbox/gsw_SA_from_rho_CT.m

@@ -16,11 +16,11 @@
 %  Calculates the Absolute Salinity of a seawater sample, for given values
 %  of its density, Conservative Temperature and sea pressure (in dbar). 
 %  This function uses the computationally-efficient 48-term expression for 
-%  density in terms of SA, CT and p (McDougall et al., 2013).
+%  density in terms of SA, CT and p (IOC et al., 2013).
 %
 %  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 McDougall et al. (2013).  The GSW library function 
+%  described in IOC et al. (2013).  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".  
 %
@@ -43,7 +43,7 @@
 % AUTHOR: 
 %  Trevor McDougall & Paul Barker                      [ [email protected] ]
 %      
-% VERSION NUMBER: 3.02 (15th November, 2012)
+% VERSION NUMBER: 3.03 (29th April, 2013)
 %
 % REFERENCES:
 %  IOC, SCOR and IAPSO, 2010: The international thermodynamic equation of 
@@ -52,11 +52,6 @@
 %   UNESCO (English), 196 pp.  Available from http://www.TEOS-10.org
 %    See section 2.5 of this TEOS-10 Manual. 
 %
-%  McDougall T.J., P.M. Barker, R. Feistel and D.R. Jackett, 2013:  A 
-%   computationally efficient 48-term expression for the density of 
-%   seawater in terms of Conservative Temperature, and related properties
-%   of seawater.  To be submitted to J. Atm. Ocean. Technol., xx, yyy-zzz.
-%
 %  Millero, F. J., R. Feistel, D. G. Wright, and T. J. McDougall, 2008: 
 %   The composition of Standard Seawater and the definition of the 
 %   Reference-Composition Salinity Scale, Deep-Sea Res. I, 55, 50-72.