Programme Matlab masse volumique
Programme Matlab masse volumique
Script Matlab de la fonction ies80 permettant de calculer la masse volumique de l'eau en fonction de sa salinité, de sa température et de la pression:
function rho=ies80(s,t,p)
% 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 vers la page cadre de l'étude et position du problème