Re: [ferret_users] hi

["Paulo B. Oliveira" <pbo@xxxxxxxxx>]

Attached is a script to compute salinity. Please check it and report if
you find any error.

Regards,

Paulo
---------------------------------------------
On Wed, 2011-04-27 at 16:30 +0630, simi wrote:
> hi all,
>  Does anybody has a script for converting conductivity data into salinity data.
> thanks
> simi
> 
> --
> Center for Atmospheric and Oceanic Science,
> Indian Institute of Science,
> Bangalore - 560 012, INDIA
> http://caos.iisc.ernet.in/
> 
> 

\canc mode verify

! Salinity of sea water as a function of Conductivity, Temperature and
! pressure - UNESCO 1983 polynomial.
!
! Input: CO = Conducticity [S/m]
!        TE = temperature [degree C]
!        PR = Pressure [dbar]
!
! REF.
!    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.

let R  = CO*10/42.914
let rt = 0.6766097+((2.00564e-2)+((1.104259e-4)+((-6.9698e-7)+1.0031e-9*TE)*TE)*TE)*TE
let Rp = 1.0+(PR*((2.070e-5)-(6.370e-10*PR)+3.989e-15*PR*PR))/(1.0+3.426e-2*TE+4.464e-4*TE*TE+((4.215e-1)-3.107e-3*TE)*R)
let rrt = R/(Rp*rt)
let del_T = TE-15.0
let Rtx = rrt^0.5
let del_S = (del_T/(1+0.0162*del_T))*(0.0005+(-0.0056+(-0.0066+(-0.0375+(0.0636+(-0.0144*Rtx))*Rtx)*Rtx)*Rtx)*Rtx)
let Sa = 0.0080+(-0.1692+(25.3851+(14.0941+(-7.0261+2.7081*Rtx)*Rtx)*Rtx)*Rtx)*Rtx+del_S

set mode verify