gibbstable2.m

WTP主要是calculate basic thermodynamic properties for water

function [J0,n0,I,J,n] = gibbstable2
% -----------------------------------------------------------------
% Author: Philippe Daigle
% Date: 6 juin 2002
% Numerical values of the coefficients 
% and exponents of the dimensionless 
% Gibbs fre energy for region 2
% -----------------------------------------------------------------
% i = 1
J0(1) = 0;
n0(1) = -0.96927686500217*10^1;%-0.96927686500217*10^1;%
% --------------------------------
% i = 2
J0(2) = 1;
n0(2) = 0.10086655968018*10^2;%0.10087275970006*10^2%
% --------------------------------
% i = 3
J0(3) = -5;
n0(3) = -0.56087911283020*10^-2;
% --------------------------------
% i = 4
J0(4) = -4;
n0(4) = 0.71452738081455*10^-1;
% --------------------------------
% i = 5
J0(5) = -3;
n0(5) = -0.40710498223928;
% --------------------------------
% i = 6
J0(6) = -2;
n0(6) = 0.14240819171444*10^1;
% --------------------------------
% i = 7
J0(7) = -1;
n0(7) = -0.43839511319450*10^1;
% --------------------------------
% i = 8
J0(8) = 2;
n0(8) = -0.28408632460772;
% --------------------------------
% i = 9
J0(9) = 3;
n0(9) = 0.21268463753307*10^-1;

% ====================================================

% i = 1
I(1) = 1;
J(1) = 0;
n(1) = -0.17731742473213*10^-2;
% --------------------------------
% i = 2
I(2) = 1;
J(2) = 1;
n(2) = -0.17834862292358*10^-1;
% --------------------------------
% i = 3
I(3) = 1;
J(3) = 2;
n(3) = -0.45996013696365*10^-1;
% --------------------------------
% i = 4
I(4) = 1;
J(4) = 3;
n(4) = -0.57581259083432*10^-1;
% --------------------------------
% i = 5
I(5) = 1;
J(5) = 6;
n(5) = -0.50325278727930*10^-1;
% --------------------------------
% i = 6
I(6) = 2;
J(6) = 1;
n(6) = -0.33032641670203*10^-4;
% --------------------------------
% i = 7
I(7) = 2;
J(7) = 2;
n(7) = -0.18948987516315*10^-3;
% --------------------------------
% i = 8
I(8) = 2;
J(8) = 4;
n(8) = -0.39392777243355*10^-2;
% --------------------------------
% i = 9
I(9) = 2;
J(9) = 7;
n(9) = -0.43797295650573*10^-1;
% --------------------------------
% i = 10
I(10) = 2;
J(10) = 36;
n(10) = -0.26674547914087*10^-4;
% --------------------------------
% i = 11
I(11) = 3;
J(11) = 0;
n(11) = 0.20481737692309*10^-7;
% --------------------------------
% i = 12
I(12) = 3;
J(12) = 1;
n(12) = 0.43870667284435*10^-6;
% --------------------------------
% i = 13
I(13) = 3;
J(13) = 3;
n(13) = -0.32277677238570*10^-4;
% --------------------------------
% i = 14
I(14) = 3;
J(14) = 6;
n(14) = -0.15033924542148*10^-2;
% --------------------------------
% i = 15
I(15) = 3;
J(15) = 35;
n(15) = -0.40668253562649*10^-1;
% --------------------------------
% i = 16
I(16) = 4;
J(16) = 1;
n(16) = -0.78847309559367*10^-9;
% --------------------------------
% i = 17
I(17) = 4;
J(17) = 2;
n(17) = 0.12790717852285*10^-7;
% --------------------------------
% i = 18
I(18) = 4;
J(18) = 3;
n(18) = 0.48225372718507*10^-6;
% --------------------------------
% i = 19
I(19) = 5;
J(19) = 7;
n(19) = 0.22922076337661*10^-5;
% --------------------------------
% i = 20
I(20) = 6;
J(20) = 3;
n(20) = -0.16714766451061*10^-10;
% --------------------------------
% i = 21
I(21) = 6;
J(21) = 16;
n(21) = -0.21171472321355*10^-2;
% --------------------------------
% i = 22
I(22) = 6;
J(22) = 35;
n(22) = -0.23895741934104*10^2;
% --------------------------------
% i = 23
I(23) = 7;
J(23) = 0;
n(23) = -0.59059564324270*10^-17;
% --------------------------------
% i = 24
I(24) = 7;
J(24) = 11;
n(24) = -0.12621808899101*10^-5;
% --------------------------------
% i = 25
I(25) = 7;
J(25) = 25;
n(25) = -0.38946842435739*10^-1;
% --------------------------------
% i = 26
I(26) = 8;
J(26) = 8;
n(26) = 0.11256211360459*10^-10;
% --------------------------------
% i = 27
I(27) = 8;
J(27) = 36;
n(27) = -0.82311340897998*10^1;
% --------------------------------
% i = 28
I(28) = 9;
J(28) = 13;
n(28) = 0.19809712802088*10^-7;
% --------------------------------
% i = 29
I(29) = 10;
J(29) = 4;
n(29) = 0.10406965210174*10^-18;
% --------------------------------
% i = 30
I(30) = 10;
J(30) = 10;
n(30) = -0.10234747095929*10^-12;
% --------------------------------
% i = 31
I(31) = 10;
J(31) = 14;
n(31) = -0.10018179379511*10^-8;
% --------------------------------
% i = 32
I(32) = 16;
J(32) = 29;
n(32) = -0.80882908646985*10^-10;
% --------------------------------
% i = 33
I(33) = 16;
J(33) = 50;
n(33) = 0.10693031879409;
% --------------------------------
% i = 34
I(34) = 18;
J(34) = 57;
n(34) = -0.33662250574171;
% --------------------------------
% i = 35
I(35) = 20;
J(35) = 20;
n(35) = 0.89185845355421*10^-24;
% --------------------------------
% i = 36
I(36) = 20;
J(36) = 35;
n(36) = 0.30629316876232*10^-12;
% --------------------------------
% i = 37
I(37) = 20;
J(37) = 48;
n(37) = -0.42002467698208*10^-5;
% --------------------------------
% i = 38
I(38) = 21;
J(38) = 21;
n(38) = -0.59056029685639*10^-25;
% --------------------------------
% i = 39
I(39) = 22;
J(39) = 53;
n(39) = 0.37826947613457*10^-5;
% --------------------------------
% i = 40
I(40) = 23;
J(40) = 39;
n(40) = -0.12768608934681*10^-14;
% --------------------------------
% i = 41
I(41) = 24;
J(41) = 26;
n(41) = 0.73087610595061*10^-28;
% --------------------------------
% i = 42
I(42) = 24;
J(42) = 40;
n(42) = 0.55414715350778*10^-16;
% --------------------------------
% i = 43
I(43) = 24;
J(43) = 58;
n(43) = -0.94369707241210*10^-6;
% --------------------------------