function [Ir,Jr,nr] = region2_gfe_r
% Description:      Loads the numerical values of the coefficients and
%                   exponents of the residual part of the dimensionless 
%                   Gibbs free energy for Region 2 (Table 11, page 14)                 
% -------------------------------------------------------------------------
% Syntax:           [Ir,Jr,nr] = region2_gfe_r
% -------------------------------------------------------------------------
% Output:           Ir          Array for the exponent parameter Ir 
%                   Jr          Array for the exponent parameter Jr  
%                   nr          Array for the coefficient parameter nr

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