Annxe 1: Formule IES80

La formule IES80 est largement utilisée pour déterminer la masse volumique de l'eau de mer, en fonction de la température, la pression et la salinité.

Elle est présentée ci dessous sou forme de script matlab.

 

 function rho=ies80(s,t,p)

t=24.5
s=35;
p=6.5

% International Equation of State of sea water : rho=ies80(s,t,p)
%   rho  density     (kg/m3)
%    s   salinity
%    t   temperature (�C)
%    p   pressure    (bars)
%
%    rho,s,t,p can be vector or matrix
%
%    rho=ies80(s,t) use p=0


r0_coef=[999.842594,6.793952e-2,-9.09529e-3,1.001685e-4,-1.120083e-6,...
         6.536332e-9,8.24493e-1,-4.0899e-3,7.6438e-5,-8.2467e-7,...
         5.3875e-9,-5.72466e-3,1.0227e-4,-1.6546e-6,4.8314e-4];

r0=polyval(r0_coef(6:-1:1),t)+...
   polyval(r0_coef(11:-1:7),t).*s+...
   polyval(r0_coef(14:-1:12),t).*s.^1.5+...
   r0_coef(15)*s.^2;

if nargin==3

K_coef=[19652.21,148.4206,-2.327105,1.360447e-2,-5.155288e-5,3.239908,...
        1.43713e-3,1.16092e-4,-5.77905e-7,8.50935e-5,...
        -6.12293e-6,5.2787e-8,54.6746,-0.603459,1.09987e-2,...
        -6.1670e-5,7.944e-2,1.6483e-2,-5.3009e-4,2.2838e-3,...
        -1.0981e-5,-1.6078e-6,1.91075e-4,-9.9348e-7,...
        2.0816e-8,9.1697e-10];

K=polyval(K_coef(5:-1:1),t)+...
  polyval(K_coef(9:-1:6),t).*p+...
  polyval(K_coef(12:-1:10),t).*p.^2+...
  polyval(K_coef(16:-1:13),t).*s+...
  polyval(K_coef(19:-1:17),t).*s.^1.5+...
  polyval(K_coef(22:-1:20),t).*p.*s+...
  K_coef(23)*p.*s.^1.5+...
  polyval(K_coef(26:-1:24),t).*p.^2.*s;
 
rho=r0./(1-p./K);

else

rho=r0;

end


Retour à la page                                                                                  Sommaire