Newer
Older
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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
function dS = sw_salds(Rtx,delT)
% SW_SALDS Differiential dS/d(sqrt(Rt)) at constant T.
%=========================================================================
% SW_SALDS $Revision: 1.3 $ $Date: 1994/10/10 05:46:08 $
% Copyright (C) CSIRO, Phil Morgan 1993.
%
% USAGE: dS = sw_salds(Rtx,delT)
%
% DESCRIPTION:
% Calculates Salinity differential dS/d(sqrt(Rt)) at constant T.
% UNESCO 1983 polynomial.
%
% INPUT: (all must have same dimensions)
% Rtx = sqrt(Rt) where Rt defined in sw_salt.m
% delT = T-15 [degree C (IPTS-68)]
%
% OUTPUT:
% dS = S differential dS/d(sqrt(Rt)) at constant T.
%
% AUTHOR: Phil Morgan 93-04-21 (morgan@ml.csiro.au)
%
% DISCLAIMER:
% This software is provided "as is" without warranty of any kind.
% See the file sw_copy.m for conditions of use and licence.
%
% REFERENCES:
% Fofonoff, P. and Millard, R.C. Jr
% Unesco 1983. Algorithms for computation of fundamental properties of
% seawater, 1983. _Unesco Tech. Pap. in Mar. Sci._, No. 44, 53 pp.
%=========================================================================
% CALLER: sw_cndr.m
% CALLEE: none
%-------------
% CHECK INPUTS
%-------------
if nargin~=2
error('sw_salds.m: must have 2 input arguments')
end %if
[m1,n1] = size(Rtx);
[m2,n2] = size(delT);
if ~(m1==m2 | n1==n2)
error('sw_salds.m: Rtx and delT must have the same shape')
end %if
%-------
% BEGIN
%-------
a0 = 0.0080;
a1 = -0.1692;
a2 = 25.3851;
a3 = 14.0941;
a4 = -7.0261;
a5 = 2.7081;
b0 = 0.0005;
b1 = -0.0056;
b2 = -0.0066;
b3 = -0.0375;
b4 = 0.0636;
b5 = -0.0144;
k = 0.0162;
dS = a1 + (2*a2 + (3*a3 + (4*a4 + 5*a5.*Rtx).*Rtx).*Rtx).*Rtx + ...
(delT./(1+k*delT))* ...
(b1 + (2*b2 + (3*b3 + (4*b4 + 5*b5.*Rtx).*Rtx).*Rtx).*Rtx);
return
%-----------------------------------------------------------------------