NCCOOS Trac Projects: Top | Web | Platforms | Processing | Viz | Sprints | Sandbox | (Wind)

root/gliderproc/trunk/MATLAB/seawater/sw_salrt.m

Revision 495 (checked in by cbc, 11 years ago)

Initial import of Stark code.

Line 
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
Note: See TracBrowser for help on using the browser.