-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathgsw_specvol_anom_standard_CT_exact.m
104 lines (93 loc) · 3.71 KB
/
gsw_specvol_anom_standard_CT_exact.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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
function specvol_anom_CT_exact = gsw_specvol_anom_standard_CT_exact(SA,CT,p)
% gsw_specvol_anom_standard_CT_exact specific volume anomaly
%==========================================================================
%
% USAGE:
% specvol_anom_CT_exact = gsw_specvol_anom_standard_CT_exact(SA,CT,p)
%
% DESCRIPTION:
% Calculates specific volume anomaly from Absolute Salinity, Conservative
% Temperature and pressure. The reference value of Absolute Salinity is
% SSO and the reference value of Conservative Temperature is equal to
% 0 degress C.
%
% Note that this function uses the full Gibbs function. There is an
% alternative to calling this function, namely
% gsw_specvol_anom_standard(SA,CT,p), which uses the computationally
% efficient 75-term expression for specific volume in terms of SA, CT and
% p (Roquet et al., 2015).
%
% INPUT:
% SA = Absolute Salinity [ g/kg ]
% CT = Conservative Temperature (ITS-90) [ deg C ]
% p = sea pressure [ dbar ]
% ( i.e. absolute pressure - 10.1325 dbar )
%
% SA & CT need to have the same dimensions.
% p may have dimensions 1x1 or Mx1 or 1xN or MxN, where SA & CT are MxN.
%
% OUTPUT:
% specvol_anom_CT_exact = specific volume anomaly [ m^3/kg ]
%
% AUTHOR:
% Trevor McDougall and Paul Barker [ help@teos-10.org ]
%
% VERSION NUMBER: 3.05 (29th 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 Eqn. (3.7.3) of this TEOS-10 Manual.
%
% Roquet, F., G. Madec, T.J. McDougall, P.M. Barker, 2015: Accurate
% polynomial expressions for the density and specifc volume of seawater
% using the TEOS-10 standard. Ocean Modelling.
%
% The software is available from http://www.TEOS-10.org
%
%==========================================================================
%--------------------------------------------------------------------------
% Check variables and resize if necessary
%--------------------------------------------------------------------------
if ~(nargin == 3)
error('gsw_specvol_anom_standard_CT_exact: Requires three inputs')
end %if
[ms,ns] = size(SA);
[mt,nt] = size(CT);
[mp,np] = size(p);
if (mt ~= ms | nt ~= ns)
error('gsw_specvol_anom_standard_CT_exact: SA and CT must have same dimensions')
end
if (mp == 1) & (np == 1) % p scalar - fill to size of SA
p = p*ones(size(SA));
elseif (ns == np) & (mp == 1) % p is row vector,
p = p(ones(1,ms), :); % copy down each column.
elseif (ms == mp) & (np == 1) % p is column vector,
p = p(:,ones(1,ns)); % copy across each row.
elseif (ns == mp) & (np == 1) % p is a transposed row vector,
p = p.'; % transposed then
p = p(ones(1,ms), :); % copy down each column.
elseif (ms == mp) & (ns == np)
% ok
else
error('gsw_specvol_anom_standard_CT_exact: Inputs array dimensions arguments do not agree')
end %if
if ms == 1
SA = SA.';
CT = CT.';
p = p.';
transposed = 1;
else
transposed = 0;
end
%--------------------------------------------------------------------------
% Start of the calculation
%--------------------------------------------------------------------------
t = gsw_t_from_CT(SA,CT,p);
specvol_anom_CT_exact = gsw_specvol_anom_standard_t_exact(SA,t,p);
if transposed
specvol_anom_CT_exact = specvol_anom_CT_exact.';
end
end