-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathgsw_SP_from_SK.m
51 lines (44 loc) · 1.6 KB
/
gsw_SP_from_SK.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
function SP = gsw_SP_from_SK(SK)
% gsw_SP_from_SK Practical Salinity from Knudsen Salinity
%==========================================================================
%
% USAGE:
% SP = gsw_SP_from_SK(SK)
%
% DESCRIPTION:
% Calculates Practical Salinity from Knudsen Salinity.
%
% INPUT:
% SK = Knudsen Salinity [parts per thousand, ppt]
%
% OUTPUT:
% SP = Practical Salinity (PSS-78) [ unitless ]
%
% AUTHOR:
% Trevor McDougall and Paul Barker [ [email protected] ]
%
% VERSION NUMBER: 3.05 (27th January 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 http://www.TEOS-10.org
% See Appendix A.3 of this TEOS-10 Manual.
%
% The software is available from http://www.TEOS-10.org
%
%==========================================================================
%--------------------------------------------------------------------------
% Check variables
%--------------------------------------------------------------------------
if ~(nargin == 1)
error('gsw_SP_from_SK: Requires only one input')
end %if
%--------------------------------------------------------------------------
% Start of the calculation
%--------------------------------------------------------------------------
SP = (SK - 0.03).*(1.80655/1.805);
% This line ensures that SP is non-negative.
SP(SP < 0) = 0;
end