1 |
|
---|
2 |
function rt = sw_salrt(T) |
---|
3 |
|
---|
4 |
% SW_SALRT Conductivity ratio rt(T) = C(35,T,0)/C(35,15,0) |
---|
5 |
%========================================================================= |
---|
6 |
% SW_SALRT $Id: sw_salrt.m,v 1.1 2003/12/12 04:23:22 pen078 Exp $ |
---|
7 |
% Copyright (C) CSIRO, Phil Morgan 1993. |
---|
8 |
% |
---|
9 |
% USAGE: rt = sw_salrt(T) |
---|
10 |
% |
---|
11 |
% DESCRIPTION: |
---|
12 |
% Equation rt(T) = C(35,T,0)/C(35,15(IPTS-68),0) used in calculating |
---|
13 |
% salinity. |
---|
14 |
% UNESCO 1983 polynomial. |
---|
15 |
% |
---|
16 |
% INPUT: |
---|
17 |
% T = temperature [degree C (ITS-90)] |
---|
18 |
% |
---|
19 |
% OUTPUT: |
---|
20 |
% rt = conductivity ratio [no units] |
---|
21 |
% |
---|
22 |
% AUTHOR: Phil Morgan 93-04-17, Lindsay Pender (Lindsay.Pender@csiro.au) |
---|
23 |
% |
---|
24 |
% DISCLAIMER: |
---|
25 |
% This software is provided "as is" without warranty of any kind. |
---|
26 |
% See the file sw_copy.m for conditions of use and licence. |
---|
27 |
% |
---|
28 |
% REFERENCES: |
---|
29 |
% Fofonoff, P. and Millard, R.C. Jr |
---|
30 |
% Unesco 1983. Algorithms for computation of fundamental properties of |
---|
31 |
% seawater, 1983. _Unesco Tech. Pap. in Mar. Sci._, No. 44, 53 pp. |
---|
32 |
%========================================================================= |
---|
33 |
|
---|
34 |
% Modifications |
---|
35 |
% 03-12-12. Lindsay Pender, Converted to ITS-90. |
---|
36 |
|
---|
37 |
% CALLER: sw_salt |
---|
38 |
% CALLEE: none |
---|
39 |
|
---|
40 |
% rt = rt(T) = C(35,T,0)/C(35,15,0) |
---|
41 |
% Eqn (3) p.7 Unesco. |
---|
42 |
|
---|
43 |
T68 = T * 1.00024; |
---|
44 |
|
---|
45 |
c0 = 0.6766097; |
---|
46 |
c1 = 2.00564e-2; |
---|
47 |
c2 = 1.104259e-4; |
---|
48 |
c3 = -6.9698e-7; |
---|
49 |
c4 = 1.0031e-9; |
---|
50 |
|
---|
51 |
rt = c0 + (c1 + (c2 + (c3 + c4.*T68).*T68).*T68).*T68; |
---|
52 |
|
---|