% dosalt   calculate salinity from Seabird Cond and Temp output
% function S = dosalt(C, T)
% Uses the CSIRO Sewater toolbox
% S = salinity
% C = microcat conductivity S/m (Siemens / meter)
% T = microcat temperature, deg. C
% P is assumed to be 14.7 psi = 10.14 db

% written by Marinna Martini 7/18/06

function S = dosalt(C, T)

% make an assumption for pressure, nominal atmospheric
P = 10.14; %db

% convert C from S /m to mS/cm
C = C.*10;

% get the ratio using in-situ conductivity
%   SW_C3515   Conductivity at (35,15,0)
%   USAGE:  c3515 = sw_c3515
%     Returns conductivity at S=35 psu , T=15 C [ITPS 68] and P=0 db).
%     c3515  = Conductivity   [mmho/cm == mS/cm] 
cndr = C / sw_c3515;    
 
%   USAGE: S = sw_salt(cndr,T,P)
%     Calculates Salinity from conductivity ratio. UNESCO 1983 polynomial.
%     cndr = Conductivity ratio     R =  C(S,T,P)/C(35,15,0) [no units]
%     T    = temperature [degree C (IPTS-68)]
%     P    = pressure    [db]
%     S    = salinity    [psu      (PSS-78)]

S = sw_salt(cndr, T, P);

return