Initializes all the constants and parameters in the GERG-2008 model.
This function sets up the GERG-2008 model by initializing various constants, parameters, and coefficients used in the model. It includes the initialization of molar masses, critical densities, critical temperatures, exponents, coefficients for pure fluid equations, mixture parameters, and ideal gas parameters.
The function does not take any parameters and does not return any value.
The function also includes code to produce nearly exact values for n0(1) and n0(2), which is not called in the current code but is included for reference. This routine must be called once before any other routine.
776{
777
778
779
781 double T0, d0;
782
784 Rs = 8.31451;
786 o13 = 1.0 / 3.0;
787
788 for (
int i = 1; i <=
MaxFlds; ++i){
790 }
792
793
815
816
817 for(
int i = 1; i <=
MaxFlds; ++i){
820 }
836
837
838 Dc[1] = 10.139342719;
840 Dc[3] = 10.624978698;
847 Dc[10] = 2.705877875;
848 Dc[11] = 2.315324434;
849 Dc[12] = 2.056404127;
855 Dc[18] = 17.87371609;
858 Dc[21] = 13.407429659;
859
860
882
883
884 for(
int i = 1; i <=
MaxFlds; ++i){
885 Vc3[i] = 1 / pow(
Dc[i], o13) / 2;
886 Tc2[i] = sqrt(
Tc[i]);
899 }
900 for (int i = 1; i <= 4; ++i){
901 if (i != 3){
926 }
927 }
928
929
930
931 noik[1][1] = 0.57335704239162;
932 noik[1][2] = -1.676068752373;
933 noik[1][3] = 0.23405291834916;
934 noik[1][4] = -0.21947376343441;
935 noik[1][5] = 0.016369201404128;
936 noik[1][6] = 0.01500440638928;
937 noik[1][7] = 0.098990489492918;
938 noik[1][8] = 0.58382770929055;
939 noik[1][9] = -0.7478686756039;
940 noik[1][10] = 0.30033302857974;
941 noik[1][11] = 0.20985543806568;
942 noik[1][12] = -0.018590151133061;
943 noik[1][13] = -0.15782558339049;
944 noik[1][14] = 0.12716735220791;
945 noik[1][15] = -0.032019743894346;
946 noik[1][16] = -0.068049729364536;
947 noik[1][17] = 0.024291412853736;
948 noik[1][18] = 5.1440451639444E-03;
949 noik[1][19] = -0.019084949733532;
950 noik[1][20] = 5.5229677241291E-03;
951 noik[1][21] = -4.4197392976085E-03;
952 noik[1][22] = 0.040061416708429;
953 noik[1][23] = -0.033752085907575;
954 noik[1][24] = -2.5127658213357E-03;
955
956 noik[2][1] = 0.59889711801201;
957 noik[2][2] = -1.6941557480731;
958 noik[2][3] = 0.24579736191718;
959 noik[2][4] = -0.23722456755175;
960 noik[2][5] = 0.017954918715141;
961 noik[2][6] = 0.014592875720215;
962 noik[2][7] = 0.10008065936206;
963 noik[2][8] = 0.73157115385532;
964 noik[2][9] = -0.88372272336366;
965 noik[2][10] = 0.31887660246708;
966 noik[2][11] = 0.20766491728799;
967 noik[2][12] = -0.019379315454158;
968 noik[2][13] = -0.16936641554983;
969 noik[2][14] = 0.13546846041701;
970 noik[2][15] = -0.033066712095307;
971 noik[2][16] = -0.060690817018557;
972 noik[2][17] = 0.012797548292871;
973 noik[2][18] = 5.8743664107299E-03;
974 noik[2][19] = -0.018451951971969;
975 noik[2][20] = 4.7226622042472E-03;
976 noik[2][21] = -5.2024079680599E-03;
977 noik[2][22] = 0.043563505956635;
978 noik[2][23] = -0.036251690750939;
979 noik[2][24] = -2.8974026866543E-03;
980
981 noik[4][1] = 0.63596780450714;
982 noik[4][2] = -1.7377981785459;
983 noik[4][3] = 0.28914060926272;
984 noik[4][4] = -0.33714276845694;
985 noik[4][5] = 0.022405964699561;
986 noik[4][6] = 0.015715424886913;
987 noik[4][7] = 0.11450634253745;
988 noik[4][8] = 1.0612049379745;
989 noik[4][9] = -1.2855224439423;
990 noik[4][10] = 0.39414630777652;
991 noik[4][11] = 0.31390924682041;
992 noik[4][12] = -0.021592277117247;
993 noik[4][13] = -0.21723666564905;
994 noik[4][14] = -0.28999574439489;
995 noik[4][15] = 0.42321173025732;
996 noik[4][16] = 0.04643410025926;
997 noik[4][17] = -0.13138398329741;
998 noik[4][18] = 0.011492850364368;
999 noik[4][19] = -0.033387688429909;
1000 noik[4][20] = 0.015183171583644;
1001 noik[4][21] = -4.7610805647657E-03;
1002 noik[4][22] = 0.046917166277885;
1003 noik[4][23] = -0.039401755804649;
1004 noik[4][24] = -3.2569956247611E-03;
1005
1006 noik[5][1] = 1.0403973107358;
1007 noik[5][2] = -2.8318404081403;
1008 noik[5][3] = 0.84393809606294;
1009 noik[5][4] = -0.076559591850023;
1010 noik[5][5] = 0.09469737305728;
1011 noik[5][6] = 2.4796475497006E-04;
1012 noik[5][7] = 0.2774376042287;
1013 noik[5][8] = -0.043846000648377;
1014 noik[5][9] = -0.2699106478435;
1015 noik[5][10] = -0.06931341308986;
1016 noik[5][11] = -0.029632145981653;
1017 noik[5][12] = 0.01404012675138;
1018
1019 noik[6][1] = 1.04293315891;
1020 noik[6][2] = -2.8184272548892;
1021 noik[6][3] = 0.8617623239785;
1022 noik[6][4] = -0.10613619452487;
1023 noik[6][5] = 0.098615749302134;
1024 noik[6][6] = 2.3948208682322E-04;
1025 noik[6][7] = 0.3033000485695;
1026 noik[6][8] = -0.041598156135099;
1027 noik[6][9] = -0.29991937470058;
1028 noik[6][10] = -0.080369342764109;
1029 noik[6][11] = -0.029761373251151;
1030 noik[6][12] = 0.01305963030314;
1031
1032 noik[7][1] = 1.0626277411455;
1033 noik[7][2] = -2.862095182835;
1034 noik[7][3] = 0.88738233403777;
1035 noik[7][4] = -0.12570581155345;
1036 noik[7][5] = 0.10286308708106;
1037 noik[7][6] = 2.5358040602654E-04;
1038 noik[7][7] = 0.32325200233982;
1039 noik[7][8] = -0.037950761057432;
1040 noik[7][9] = -0.32534802014452;
1041 noik[7][10] = -0.079050969051011;
1042 noik[7][11] = -0.020636720547775;
1043 noik[7][12] = 0.005705380933475;
1044
1045 noik[8][1] = 1.0963;
1046 noik[8][2] = -3.0402;
1047 noik[8][3] = 1.0317;
1048 noik[8][4] = -0.1541;
1049 noik[8][5] = 0.11535;
1050 noik[8][6] = 0.00029809;
1051 noik[8][7] = 0.39571;
1052 noik[8][8] = -0.045881;
1053 noik[8][9] = -0.35804;
1054 noik[8][10] = -0.10107;
1055 noik[8][11] = -0.035484;
1056 noik[8][12] = 0.018156;
1057
1058 noik[9][1] = 1.0968643098001;
1059 noik[9][2] = -2.9988888298061;
1060 noik[9][3] = 0.99516886799212;
1061 noik[9][4] = -0.16170708558539;
1062 noik[9][5] = 0.11334460072775;
1063 noik[9][6] = 2.6760595150748E-04;
1064 noik[9][7] = 0.40979881986931;
1065 noik[9][8] = -0.040876423083075;
1066 noik[9][9] = -0.38169482469447;
1067 noik[9][10] = -0.10931956843993;
1068 noik[9][11] = -0.03207322332799;
1069 noik[9][12] = 0.016877016216975;
1070
1071 noik[10][1] = 1.0553238013661;
1072 noik[10][2] = -2.6120615890629;
1073 noik[10][3] = 0.7661388296726;
1074 noik[10][4] = -0.29770320622459;
1075 noik[10][5] = 0.11879907733358;
1076 noik[10][6] = 2.7922861062617E-04;
1077 noik[10][7] = 0.46347589844105;
1078 noik[10][8] = 0.011433196980297;
1079 noik[10][9] = -0.48256968738131;
1080 noik[10][10] = -0.093750558924659;
1081 noik[10][11] = -6.7273247155994E-03;
1082 noik[10][12] = -5.1141583585428E-03;
1083
1084 noik[11][1] = 1.0543747645262;
1085 noik[11][2] = -2.6500681506144;
1086 noik[11][3] = 0.81730047827543;
1087 noik[11][4] = -0.30451391253428;
1088 noik[11][5] = 0.122538687108;
1089 noik[11][6] = 2.7266472743928E-04;
1090 noik[11][7] = 0.4986582568167;
1091 noik[11][8] = -7.1432815084176E-04;
1092 noik[11][9] = -0.5423689552545;
1093 noik[11][10] = -0.13801821610756;
1094 noik[11][11] = -6.1595287380011E-03;
1095 noik[11][12] = 4.8602510393022E-04;
1096
1097 noik[12][1] = 1.0722544875633;
1098 noik[12][2] = -2.4632951172003;
1099 noik[12][3] = 0.65386674054928;
1100 noik[12][4] = -0.36324974085628;
1101 noik[12][5] = 0.12713269626764;
1102 noik[12][6] = 3.071357277793E-04;
1103 noik[12][7] = 0.5265685698754;
1104 noik[12][8] = 0.019362862857653;
1105 noik[12][9] = -0.58939426849155;
1106 noik[12][10] = -0.14069963991934;
1107 noik[12][11] = -7.8966330500036E-03;
1108 noik[12][12] = 3.3036597968109E-03;
1109
1110 noik[13][1] = 1.1151;
1111 noik[13][2] = -2.702;
1112 noik[13][3] = 0.83416;
1113 noik[13][4] = -0.38828;
1114 noik[13][5] = 0.1376;
1115 noik[13][6] = 0.00028185;
1116 noik[13][7] = 0.62037;
1117 noik[13][8] = 0.015847;
1118 noik[13][9] = -0.61726;
1119 noik[13][10] = -0.15043;
1120 noik[13][11] = -0.012982;
1121 noik[13][12] = 0.0044325;
1122
1123 noik[14][1] = 1.0461;
1124 noik[14][2] = -2.4807;
1125 noik[14][3] = 0.74372;
1126 noik[14][4] = -0.52579;
1127 noik[14][5] = 0.15315;
1128 noik[14][6] = 0.00032865;
1129 noik[14][7] = 0.84178;
1130 noik[14][8] = 0.055424;
1131 noik[14][9] = -0.73555;
1132 noik[14][10] = -0.18507;
1133 noik[14][11] = -0.020775;
1134 noik[14][12] = 0.012335;
1135
1136 noik[16][1] = 0.88878286369701;
1137 noik[16][2] = -2.4879433312148;
1138 noik[16][3] = 0.59750190775886;
1139 noik[16][4] = 9.6501817061881E-03;
1140 noik[16][5] = 0.07197042871277;
1141 noik[16][6] = 2.2337443000195E-04;
1142 noik[16][7] = 0.18558686391474;
1143 noik[16][8] = -0.03812936803576;
1144 noik[16][9] = -0.15352245383006;
1145 noik[16][10] = -0.026726814910919;
1146 noik[16][11] = -0.025675298677127;
1147 noik[16][12] = 9.5714302123668E-03;
1148
1149 noik[17][1] = 0.90554;
1150 noik[17][2] = -2.4515;
1151 noik[17][3] = 0.53149;
1152 noik[17][4] = 0.024173;
1153 noik[17][5] = 0.072156;
1154 noik[17][6] = 0.00018818;
1155 noik[17][7] = 0.19405;
1156 noik[17][8] = -0.043268;
1157 noik[17][9] = -0.12778;
1158 noik[17][10] = -0.027896;
1159 noik[17][11] = -0.034154;
1160 noik[17][12] = 0.016329;
1161
1162 noik[19][1] = 0.87641;
1163 noik[19][2] = -2.0367;
1164 noik[19][3] = 0.21634;
1165 noik[19][4] = -0.050199;
1166 noik[19][5] = 0.066994;
1167 noik[19][6] = 0.00019076;
1168 noik[19][7] = 0.20227;
1169 noik[19][8] = -0.0045348;
1170 noik[19][9] = -0.2223;
1171 noik[19][10] = -0.034714;
1172 noik[19][11] = -0.014885;
1173 noik[19][12] = 0.0074154;
1174
1175 noik[21][1] = 0.85095714803969;
1176 noik[21][2] = -2.400322294348;
1177 noik[21][3] = 0.54127841476466;
1178 noik[21][4] = 0.016919770692538;
1179 noik[21][5] = 0.068825965019035;
1180 noik[21][6] = 2.1428032815338E-04;
1181 noik[21][7] = 0.17429895321992;
1182 noik[21][8] = -0.033654495604194;
1183 noik[21][9] = -0.13526799857691;
1184 noik[21][10] = -0.016387350791552;
1185 noik[21][11] = -0.024987666851475;
1186 noik[21][12] = 8.8769204815709E-03;
1187
1188 coik[3][1] = 0;
doik[3][1] = 1;
toik[3][1] = 0;
noik[3][1] = 0.52646564804653;
1189 coik[3][2] = 0;
doik[3][2] = 1;
toik[3][2] = 1.25;
noik[3][2] = -1.4995725042592;
1190 coik[3][3] = 0;
doik[3][3] = 2;
toik[3][3] = 1.625;
noik[3][3] = 0.27329786733782;
1191 coik[3][4] = 0;
doik[3][4] = 3;
toik[3][4] = 0.375;
noik[3][4] = 0.12949500022786;
1192 coik[3][5] = 1;
doik[3][5] = 3;
toik[3][5] = 0.375;
noik[3][5] = 0.15404088341841;
1193 coik[3][6] = 1;
doik[3][6] = 3;
toik[3][6] = 1.375;
noik[3][6] = -0.58186950946814;
1194 coik[3][7] = 1;
doik[3][7] = 4;
toik[3][7] = 1.125;
noik[3][7] = -0.18022494838296;
1195 coik[3][8] = 1;
doik[3][8] = 5;
toik[3][8] = 1.375;
noik[3][8] = -0.095389904072812;
1196 coik[3][9] = 1;
doik[3][9] = 6;
toik[3][9] = 0.125;
noik[3][9] = -8.0486819317679E-03;
1197 coik[3][10] = 1;
doik[3][10] = 6;
toik[3][10] = 1.625;
noik[3][10] = -0.03554775127309;
1198 coik[3][11] = 2;
doik[3][11] = 1;
toik[3][11] = 3.75;
noik[3][11] = -0.28079014882405;
1199 coik[3][12] = 2;
doik[3][12] = 4;
toik[3][12] = 3.5;
noik[3][12] = -0.082435890081677;
1200 coik[3][13] = 3;
doik[3][13] = 1;
toik[3][13] = 7.5;
noik[3][13] = 0.010832427979006;
1201 coik[3][14] = 3;
doik[3][14] = 1;
toik[3][14] = 8;
noik[3][14] = -6.7073993161097E-03;
1202 coik[3][15] = 3;
doik[3][15] = 3;
toik[3][15] = 6;
noik[3][15] = -4.6827907600524E-03;
1203 coik[3][16] = 3;
doik[3][16] = 3;
toik[3][16] = 16;
noik[3][16] = -0.028359911832177;
1204 coik[3][17] = 3;
doik[3][17] = 4;
toik[3][17] = 11;
noik[3][17] = 0.019500174744098;
1205 coik[3][18] = 5;
doik[3][18] = 5;
toik[3][18] = 24;
noik[3][18] = -0.21609137507166;
1206 coik[3][19] = 5;
doik[3][19] = 5;
toik[3][19] = 26;
noik[3][19] = 0.43772794926972;
1207 coik[3][20] = 5;
doik[3][20] = 5;
toik[3][20] = 28;
noik[3][20] = -0.22130790113593;
1208 coik[3][21] = 6;
doik[3][21] = 5;
toik[3][21] = 24;
noik[3][21] = 0.015190189957331;
1209 coik[3][22] = 6;
doik[3][22] = 5;
toik[3][22] = 26;
noik[3][22] = -0.0153809489533;
1210
1211 coik[15][1] = 0;
doik[15][1] = 1;
toik[15][1] = 0.5;
noik[15][1] = 5.3579928451252;
1212 coik[15][2] = 0;
doik[15][2] = 1;
toik[15][2] = 0.625;
noik[15][2] = -6.2050252530595;
1213 coik[15][3] = 0;
doik[15][3] = 2;
toik[15][3] = 0.375;
noik[15][3] = 0.13830241327086;
1214 coik[15][4] = 0;
doik[15][4] = 2;
toik[15][4] = 0.625;
noik[15][4] = -0.071397954896129;
1215 coik[15][5] = 0;
doik[15][5] = 4;
toik[15][5] = 1.125;
noik[15][5] = 0.015474053959733;
1216 coik[15][6] = 1;
doik[15][6] = 1;
toik[15][6] = 2.625;
noik[15][6] = -0.14976806405771;
1217 coik[15][7] = 1;
doik[15][7] = 5;
toik[15][7] = 0;
noik[15][7] = -0.026368723988451;
1218 coik[15][8] = 1;
doik[15][8] = 5;
toik[15][8] = 0.25;
noik[15][8] = 0.056681303156066;
1219 coik[15][9] = 1;
doik[15][9] = 5;
toik[15][9] = 1.375;
noik[15][9] = -0.060063958030436;
1220 coik[15][10] = 2;
doik[15][10] = 1;
toik[15][10] = 4;
noik[15][10] = -0.45043942027132;
1221 coik[15][11] = 2;
doik[15][11] = 1;
toik[15][11] = 4.25;
noik[15][11] = 0.424788402445;
1222 coik[15][12] = 3;
doik[15][12] = 2;
toik[15][12] = 5;
noik[15][12] = -0.021997640827139;
1223 coik[15][13] = 3;
doik[15][13] = 5;
toik[15][13] = 8;
noik[15][13] = -0.01049952137453;
1224 coik[15][14] = 5;
doik[15][14] = 1;
toik[15][14] = 8;
noik[15][14] = -2.8955902866816E-03;
1225
1226 coik[18][1] = 0;
doik[18][1] = 1;
toik[18][1] = 0.5;
noik[18][1] = 0.82728408749586;
1227 coik[18][2] = 0;
doik[18][2] = 1;
toik[18][2] = 1.25;
noik[18][2] = -1.8602220416584;
1228 coik[18][3] = 0;
doik[18][3] = 1;
toik[18][3] = 1.875;
noik[18][3] = -1.1199009613744;
1229 coik[18][4] = 0;
doik[18][4] = 2;
toik[18][4] = 0.125;
noik[18][4] = 0.15635753976056;
1230 coik[18][5] = 0;
doik[18][5] = 2;
toik[18][5] = 1.5;
noik[18][5] = 0.87375844859025;
1231 coik[18][6] = 0;
doik[18][6] = 3;
toik[18][6] = 1;
noik[18][6] = -0.36674403715731;
1232 coik[18][7] = 0;
doik[18][7] = 4;
toik[18][7] = 0.75;
noik[18][7] = 0.053987893432436;
1233 coik[18][8] = 1;
doik[18][8] = 1;
toik[18][8] = 1.5;
noik[18][8] = 1.0957690214499;
1234 coik[18][9] = 1;
doik[18][9] = 5;
toik[18][9] = 0.625;
noik[18][9] = 0.053213037828563;
1235 coik[18][10] = 1;
doik[18][10] = 5;
toik[18][10] = 2.625;
noik[18][10] = 0.013050533930825;
1236 coik[18][11] = 2;
doik[18][11] = 1;
toik[18][11] = 5;
noik[18][11] = -0.41079520434476;
1237 coik[18][12] = 2;
doik[18][12] = 2;
toik[18][12] = 4;
noik[18][12] = 0.1463744334412;
1238 coik[18][13] = 2;
doik[18][13] = 4;
toik[18][13] = 4.5;
noik[18][13] = -0.055726838623719;
1239 coik[18][14] = 3;
doik[18][14] = 4;
toik[18][14] = 3;
noik[18][14] = -0.0112017741438;
1240 coik[18][15] = 5;
doik[18][15] = 1;
toik[18][15] = 4;
noik[18][15] = -6.6062758068099E-03;
1241 coik[18][16] = 5;
doik[18][16] = 1;
toik[18][16] = 6;
noik[18][16] = 4.6918522004538E-03;
1242
1243 coik[20][1] = 0;
doik[20][1] = 1;
toik[20][1] = 0;
noik[20][1] = -0.45579024006737;
1244 coik[20][2] = 0;
doik[20][2] = 1;
toik[20][2] = 0.125;
noik[20][2] = 1.2516390754925;
1245 coik[20][3] = 0;
doik[20][3] = 1;
toik[20][3] = 0.75;
noik[20][3] = -1.5438231650621;
1246 coik[20][4] = 0;
doik[20][4] = 4;
toik[20][4] = 1;
noik[20][4] = 0.020467489707221;
1247 coik[20][5] = 1;
doik[20][5] = 1;
toik[20][5] = 0.75;
noik[20][5] = -0.34476212380781;
1248 coik[20][6] = 1;
doik[20][6] = 3;
toik[20][6] = 2.625;
noik[20][6] = -0.020858459512787;
1249 coik[20][7] = 1;
doik[20][7] = 5;
toik[20][7] = 0.125;
noik[20][7] = 0.016227414711778;
1250 coik[20][8] = 1;
doik[20][8] = 5;
toik[20][8] = 1.25;
noik[20][8] = -0.057471818200892;
1251 coik[20][9] = 1;
doik[20][9] = 5;
toik[20][9] = 2;
noik[20][9] = 0.019462416430715;
1252 coik[20][10] = 2;
doik[20][10] = 2;
toik[20][10] = 1;
noik[20][10] = -0.03329568012302;
1253 coik[20][11] = 3;
doik[20][11] = 1;
toik[20][11] = 4.5;
noik[20][11] = -0.010863577372367;
1254 coik[20][12] = 3;
doik[20][12] = 2;
toik[20][12] = 5;
noik[20][12] = -0.022173365245954;
1255
1256
1257
1258 dijk[3][1] = 1;
tijk[3][1] = 0;
cijk[3][1] = 0;
eijk[3][1] = 0; bijk[3][1] = 0;
gijk[3][1] = 0;
nijk[3][1] = -9.8038985517335E-03;
1259 dijk[3][2] = 4;
tijk[3][2] = 1.85;
cijk[3][2] = 0;
eijk[3][2] = 0; bijk[3][2] = 0;
gijk[3][2] = 0;
nijk[3][2] = 4.2487270143005E-04;
1260 dijk[3][3] = 1;
tijk[3][3] = 7.85;
cijk[3][3] = 1;
eijk[3][3] = 0.5; bijk[3][3] = 1;
gijk[3][3] = 0.5;
nijk[3][3] = -0.034800214576142;
1261 dijk[3][4] = 2;
tijk[3][4] = 5.4;
cijk[3][4] = 1;
eijk[3][4] = 0.5; bijk[3][4] = 1;
gijk[3][4] = 0.5;
nijk[3][4] = -0.13333813013896;
1262 dijk[3][5] = 2;
tijk[3][5] = 0;
cijk[3][5] = 0.25;
eijk[3][5] = 0.5; bijk[3][5] = 2.5;
gijk[3][5] = 0.5;
nijk[3][5] = -0.011993694974627;
1263 dijk[3][6] = 2;
tijk[3][6] = 0.75;
cijk[3][6] = 0;
eijk[3][6] = 0.5; bijk[3][6] = 3;
gijk[3][6] = 0.5;
nijk[3][6] = 0.069243379775168;
1264 dijk[3][7] = 2;
tijk[3][7] = 2.8;
cijk[3][7] = 0;
eijk[3][7] = 0.5; bijk[3][7] = 3;
gijk[3][7] = 0.5;
nijk[3][7] = -0.31022508148249;
1265 dijk[3][8] = 2;
tijk[3][8] = 4.45;
cijk[3][8] = 0;
eijk[3][8] = 0.5; bijk[3][8] = 3;
gijk[3][8] = 0.5;
nijk[3][8] = 0.24495491753226;
1266 dijk[3][9] = 3;
tijk[3][9] = 4.25;
cijk[3][9] = 0;
eijk[3][9] = 0.5; bijk[3][9] = 3;
gijk[3][9] = 0.5;
nijk[3][9] = 0.22369816716981;
1267
1268 dijk[4][1] = 1;
tijk[4][1] = 2.6;
cijk[4][1] = 0;
eijk[4][1] = 0; bijk[4][1] = 0;
gijk[4][1] = 0;
nijk[4][1] = -0.10859387354942;
1269 dijk[4][2] = 2;
tijk[4][2] = 1.95;
cijk[4][2] = 0;
eijk[4][2] = 0; bijk[4][2] = 0;
gijk[4][2] = 0;
nijk[4][2] = 0.080228576727389;
1270 dijk[4][3] = 3;
tijk[4][3] = 0;
cijk[4][3] = 0;
eijk[4][3] = 0; bijk[4][3] = 0;
gijk[4][3] = 0;
nijk[4][3] = -9.3303985115717E-03;
1271 dijk[4][4] = 1;
tijk[4][4] = 3.95;
cijk[4][4] = 1;
eijk[4][4] = 0.5; bijk[4][4] = 1;
gijk[4][4] = 0.5;
nijk[4][4] = 0.040989274005848;
1272 dijk[4][5] = 2;
tijk[4][5] = 7.95;
cijk[4][5] = 0.5;
eijk[4][5] = 0.5; bijk[4][5] = 2;
gijk[4][5] = 0.5;
nijk[4][5] = -0.24338019772494;
1273 dijk[4][6] = 3;
tijk[4][6] = 8;
cijk[4][6] = 0;
eijk[4][6] = 0.5; bijk[4][6] = 3;
gijk[4][6] = 0.5;
nijk[4][6] = 0.23855347281124;
1274
1275 dijk[1][1] = 3;
tijk[1][1] = 0.65;
cijk[1][1] = 0;
eijk[1][1] = 0; bijk[1][1] = 0;
gijk[1][1] = 0;
nijk[1][1] = -8.0926050298746E-04;
1276 dijk[1][2] = 4;
tijk[1][2] = 1.55;
cijk[1][2] = 0;
eijk[1][2] = 0; bijk[1][2] = 0;
gijk[1][2] = 0;
nijk[1][2] = -7.5381925080059E-04;
1277 dijk[1][3] = 1;
tijk[1][3] = 3.1;
cijk[1][3] = 1;
eijk[1][3] = 0.5; bijk[1][3] = 1;
gijk[1][3] = 0.5;
nijk[1][3] = -0.041618768891219;
1278 dijk[1][4] = 2;
tijk[1][4] = 5.9;
cijk[1][4] = 1;
eijk[1][4] = 0.5; bijk[1][4] = 1;
gijk[1][4] = 0.5;
nijk[1][4] = -0.23452173681569;
1279 dijk[1][5] = 2;
tijk[1][5] = 7.05;
cijk[1][5] = 1;
eijk[1][5] = 0.5; bijk[1][5] = 1;
gijk[1][5] = 0.5;
nijk[1][5] = 0.14003840584586;
1280 dijk[1][6] = 2;
tijk[1][6] = 3.35;
cijk[1][6] = 0.875;
eijk[1][6] = 0.5; bijk[1][6] = 1.25;
gijk[1][6] = 0.5;
nijk[1][6] = 0.063281744807738;
1281 dijk[1][7] = 2;
tijk[1][7] = 1.2;
cijk[1][7] = 0.75;
eijk[1][7] = 0.5; bijk[1][7] = 1.5;
gijk[1][7] = 0.5;
nijk[1][7] = -0.034660425848809;
1282 dijk[1][8] = 2;
tijk[1][8] = 5.8;
cijk[1][8] = 0.5;
eijk[1][8] = 0.5; bijk[1][8] = 2;
gijk[1][8] = 0.5;
nijk[1][8] = -0.23918747334251;
1283 dijk[1][9] = 2;
tijk[1][9] = 2.7;
cijk[1][9] = 0;
eijk[1][9] = 0.5; bijk[1][9] = 3;
gijk[1][9] = 0.5;
nijk[1][9] = 1.9855255066891E-03;
1284 dijk[1][10] = 3;
tijk[1][10] = 0.45;
cijk[1][10] = 0;
eijk[1][10] = 0.5; bijk[1][10] = 3;
gijk[1][10] = 0.5;
nijk[1][10] = 6.1777746171555;
1285 dijk[1][11] = 3;
tijk[1][11] = 0.55;
cijk[1][11] = 0;
eijk[1][11] = 0.5; bijk[1][11] = 3;
gijk[1][11] = 0.5;
nijk[1][11] = -6.9575358271105;
1286 dijk[1][12] = 3;
tijk[1][12] = 1.95;
cijk[1][12] = 0;
eijk[1][12] = 0.5; bijk[1][12] = 3;
gijk[1][12] = 0.5;
nijk[1][12] = 1.0630185306388;
1287
1288 dijk[2][1] = 3;
tijk[2][1] = 1.85;
cijk[2][1] = 0;
eijk[2][1] = 0; bijk[2][1] = 0;
gijk[2][1] = 0;
nijk[2][1] = 0.013746429958576;
1289 dijk[2][2] = 3;
tijk[2][2] = 3.95;
cijk[2][2] = 0;
eijk[2][2] = 0; bijk[2][2] = 0;
gijk[2][2] = 0;
nijk[2][2] = -7.4425012129552E-03;
1290 dijk[2][3] = 4;
tijk[2][3] = 0;
cijk[2][3] = 0;
eijk[2][3] = 0; bijk[2][3] = 0;
gijk[2][3] = 0;
nijk[2][3] = -4.5516600213685E-03;
1291 dijk[2][4] = 4;
tijk[2][4] = 1.85;
cijk[2][4] = 0;
eijk[2][4] = 0; bijk[2][4] = 0;
gijk[2][4] = 0;
nijk[2][4] = -5.4546603350237E-03;
1292 dijk[2][5] = 4;
tijk[2][5] = 3.85;
cijk[2][5] = 0;
eijk[2][5] = 0; bijk[2][5] = 0;
gijk[2][5] = 0;
nijk[2][5] = 2.3682016824471E-03;
1293 dijk[2][6] = 1;
tijk[2][6] = 5.25;
cijk[2][6] = 0.25;
eijk[2][6] = 0.5; bijk[2][6] = 0.75;
gijk[2][6] = 0.5;
nijk[2][6] = 0.18007763721438;
1294 dijk[2][7] = 1;
tijk[2][7] = 3.85;
cijk[2][7] = 0.25;
eijk[2][7] = 0.5; bijk[2][7] = 1;
gijk[2][7] = 0.5;
nijk[2][7] = -0.44773942932486;
1295 dijk[2][8] = 1;
tijk[2][8] = 0.2;
cijk[2][8] = 0;
eijk[2][8] = 0.5; bijk[2][8] = 2;
gijk[2][8] = 0.5;
nijk[2][8] = 0.0193273748882;
1296 dijk[2][9] = 2;
tijk[2][9] = 6.5;
cijk[2][9] = 0;
eijk[2][9] = 0.5; bijk[2][9] = 3;
gijk[2][9] = 0.5;
nijk[2][9] = -0.30632197804624;
1297
1298 dijk[5][1] = 2;
tijk[5][1] = 1.85;
cijk[5][1] = 0;
eijk[5][1] = 0; bijk[5][1] = 0;
gijk[5][1] = 0;
nijk[5][1] = 0.28661625028399;
1299 dijk[5][2] = 3;
tijk[5][2] = 1.4;
cijk[5][2] = 0;
eijk[5][2] = 0; bijk[5][2] = 0;
gijk[5][2] = 0;
nijk[5][2] = -0.10919833861247;
1300 dijk[5][3] = 1;
tijk[5][3] = 3.2;
cijk[5][3] = 0.25;
eijk[5][3] = 0.5; bijk[5][3] = 0.75;
gijk[5][3] = 0.5;
nijk[5][3] = -1.137403208227;
1301 dijk[5][4] = 1;
tijk[5][4] = 2.5;
cijk[5][4] = 0.25;
eijk[5][4] = 0.5; bijk[5][4] = 1;
gijk[5][4] = 0.5;
nijk[5][4] = 0.76580544237358;
1302 dijk[5][5] = 1;
tijk[5][5] = 8;
cijk[5][5] = 0;
eijk[5][5] = 0.5; bijk[5][5] = 2;
gijk[5][5] = 0.5;
nijk[5][5] = 4.2638000926819E-03;
1303 dijk[5][6] = 2;
tijk[5][6] = 3.75;
cijk[5][6] = 0;
eijk[5][6] = 0.5; bijk[5][6] = 3;
gijk[5][6] = 0.5;
nijk[5][6] = 0.17673538204534;
1304
1305 dijk[6][1] = 2;
tijk[6][1] = 0;
cijk[6][1] = 0;
eijk[6][1] = 0; bijk[6][1] = 0;
gijk[6][1] = 0;
nijk[6][1] = -0.47376518126608;
1306 dijk[6][2] = 2;
tijk[6][2] = 0.05;
cijk[6][2] = 0;
eijk[6][2] = 0; bijk[6][2] = 0;
gijk[6][2] = 0;
nijk[6][2] = 0.48961193461001;
1307 dijk[6][3] = 3;
tijk[6][3] = 0;
cijk[6][3] = 0;
eijk[6][3] = 0; bijk[6][3] = 0;
gijk[6][3] = 0;
nijk[6][3] = -5.7011062090535E-03;
1308 dijk[6][4] = 1;
tijk[6][4] = 3.65;
cijk[6][4] = 1;
eijk[6][4] = 0.5; bijk[6][4] = 1;
gijk[6][4] = 0.5;
nijk[6][4] = -0.1996682004132;
1309 dijk[6][5] = 2;
tijk[6][5] = 4.9;
cijk[6][5] = 1;
eijk[6][5] = 0.5; bijk[6][5] = 1;
gijk[6][5] = 0.5;
nijk[6][5] = -0.69411103101723;
1310 dijk[6][6] = 2;
tijk[6][6] = 4.45;
cijk[6][6] = 0.875;
eijk[6][6] = 0.5; bijk[6][6] = 1.25;
gijk[6][6] = 0.5;
nijk[6][6] = 0.69226192739021;
1311
1312 dijk[7][1] = 1;
tijk[7][1] = 2;
cijk[7][1] = 0;
eijk[7][1] = 0; bijk[7][1] = 0;
gijk[7][1] = 0;
nijk[7][1] = -0.25157134971934;
1313 dijk[7][2] = 3;
tijk[7][2] = -1;
cijk[7][2] = 0;
eijk[7][2] = 0; bijk[7][2] = 0;
gijk[7][2] = 0;
nijk[7][2] = -6.2203841111983E-03;
1314 dijk[7][3] = 3;
tijk[7][3] = 1.75;
cijk[7][3] = 0;
eijk[7][3] = 0; bijk[7][3] = 0;
gijk[7][3] = 0;
nijk[7][3] = 0.088850315184396;
1315 dijk[7][4] = 4;
tijk[7][4] = 1.4;
cijk[7][4] = 0;
eijk[7][4] = 0; bijk[7][4] = 0;
gijk[7][4] = 0;
nijk[7][4] = -0.035592212573239;
1316
1317
1318 dijk[10][1] = 1;
tijk[10][1] = 1;
cijk[10][1] = 0;
eijk[10][1] = 0; bijk[10][1] = 0;
gijk[10][1] = 0;
nijk[10][1] = 2.5574776844118;
1319 dijk[10][2] = 1;
tijk[10][2] = 1.55;
cijk[10][2] = 0;
eijk[10][2] = 0; bijk[10][2] = 0;
gijk[10][2] = 0;
nijk[10][2] = -7.9846357136353;
1320 dijk[10][3] = 1;
tijk[10][3] = 1.7;
cijk[10][3] = 0;
eijk[10][3] = 0; bijk[10][3] = 0;
gijk[10][3] = 0;
nijk[10][3] = 4.7859131465806;
1321 dijk[10][4] = 2;
tijk[10][4] = 0.25;
cijk[10][4] = 0;
eijk[10][4] = 0; bijk[10][4] = 0;
gijk[10][4] = 0;
nijk[10][4] = -0.73265392369587;
1322 dijk[10][5] = 2;
tijk[10][5] = 1.35;
cijk[10][5] = 0;
eijk[10][5] = 0; bijk[10][5] = 0;
gijk[10][5] = 0;
nijk[10][5] = 1.3805471345312;
1323 dijk[10][6] = 3;
tijk[10][6] = 0;
cijk[10][6] = 0;
eijk[10][6] = 0; bijk[10][6] = 0;
gijk[10][6] = 0;
nijk[10][6] = 0.28349603476365;
1324 dijk[10][7] = 3;
tijk[10][7] = 1.25;
cijk[10][7] = 0;
eijk[10][7] = 0; bijk[10][7] = 0;
gijk[10][7] = 0;
nijk[10][7] = -0.49087385940425;
1325 dijk[10][8] = 4;
tijk[10][8] = 0;
cijk[10][8] = 0;
eijk[10][8] = 0; bijk[10][8] = 0;
gijk[10][8] = 0;
nijk[10][8] = -0.10291888921447;
1326 dijk[10][9] = 4;
tijk[10][9] = 0.7;
cijk[10][9] = 0;
eijk[10][9] = 0; bijk[10][9] = 0;
gijk[10][9] = 0;
nijk[10][9] = 0.11836314681968;
1327 dijk[10][10] = 4;
tijk[10][10] = 5.4;
cijk[10][10] = 0;
eijk[10][10] = 0; bijk[10][10] = 0;
gijk[10][10] = 0;
nijk[10][10] = 5.5527385721943E-05;
1328
1329
1337 fij[1][6] = 0.771035405688;
1339 fij[4][5] = 0.13042476515;
1340 fij[4][6] = 0.260632376098;
1341 fij[4][7] = 0.281570073085;
1342 fij[5][6] = -0.0551609771024;
1343 fij[5][7] = 0.0312572600489;
1344 fij[6][7] = -0.0551240293009;
1345
1346
1347 for(
int i = 1; i <=
MaxFlds; ++i){
1349 for (
int j = i + 1; j <=
MaxFlds; ++j){
1353 }
1354 }
1355
1356
1372
1373
1374 n0i[1][3] = 4.00088;
n0i[1][4] = 0.76315;
n0i[1][5] = 0.0046;
n0i[1][6] = 8.74432;
n0i[1][7] = -4.46921;
n0i[1][1] = 29.83843397;
n0i[1][2] = -15999.69151;
1375 n0i[2][3] = 3.50031;
n0i[2][4] = 0.13732;
n0i[2][5] = -0.1466;
n0i[2][6] = 0.90066;
n0i[2][7] = 0;
n0i[2][1] = 17.56770785;
n0i[2][2] = -2801.729072;
1376 n0i[3][3] = 3.50002;
n0i[3][4] = 2.04452;
n0i[3][5] = -1.06044;
n0i[3][6] = 2.03366;
n0i[3][7] = 0.01393;
n0i[3][1] = 20.65844696;
n0i[3][2] = -4902.171516;
1377 n0i[4][3] = 4.00263;
n0i[4][4] = 4.33939;
n0i[4][5] = 1.23722;
n0i[4][6] = 13.1974;
n0i[4][7] = -6.01989;
n0i[4][1] = 36.73005938;
n0i[4][2] = -23639.65301;
1378 n0i[5][3] = 4.02939;
n0i[5][4] = 6.60569;
n0i[5][5] = 3.197;
n0i[5][6] = 19.1921;
n0i[5][7] = -8.37267;
n0i[5][1] = 44.70909619;
n0i[5][2] = -31236.63551;
1379 n0i[6][3] = 4.06714;
n0i[6][4] = 8.97575;
n0i[6][5] = 5.25156;
n0i[6][6] = 25.1423;
n0i[6][7] = 16.1388;
n0i[6][1] = 34.30180349;
n0i[6][2] = -38525.50276;
1380 n0i[7][3] = 4.33944;
n0i[7][4] = 9.44893;
n0i[7][5] = 6.89406;
n0i[7][6] = 24.4618;
n0i[7][7] = 14.7824;
n0i[7][1] = 36.53237783;
n0i[7][2] = -38957.80933;
1381 n0i[8][3] = 4;
n0i[8][4] = 11.7618;
n0i[8][5] = 20.1101;
n0i[8][6] = 33.1688;
n0i[8][7] = 0;
n0i[8][1] = 43.17218626;
n0i[8][2] = -51198.30946;
1382 n0i[9][3] = 4;
n0i[9][4] = 8.95043;
n0i[9][5] = 21.836;
n0i[9][6] = 33.4032;
n0i[9][7] = 0;
n0i[9][1] = 42.67837089;
n0i[9][2] = -45215.83;
1383 n0i[10][3] = 4;
n0i[10][4] = 11.6977;
n0i[10][5] = 26.8142;
n0i[10][6] = 38.6164;
n0i[10][7] = 0;
n0i[10][1] = 46.99717188;
n0i[10][2] = -52746.83318;
1384 n0i[11][3] = 4;
n0i[11][4] = 13.7266;
n0i[11][5] = 30.4707;
n0i[11][6] = 43.5561;
n0i[11][7] = 0;
n0i[11][1] = 52.07631631;
n0i[11][2] = -57104.81056;
1385 n0i[12][3] = 4;
n0i[12][4] = 15.6865;
n0i[12][5] = 33.8029;
n0i[12][6] = 48.1731;
n0i[12][7] = 0;
n0i[12][1] = 57.25830934;
n0i[12][2] = -60546.76385;
1386 n0i[13][3] = 4;
n0i[13][4] = 18.0241;
n0i[13][5] = 38.1235;
n0i[13][6] = 53.3415;
n0i[13][7] = 0;
n0i[13][1] = 62.09646901;
n0i[13][2] = -66600.12837;
1387 n0i[14][3] = 4;
n0i[14][4] = 21.0069;
n0i[14][5] = 43.4931;
n0i[14][6] = 58.3657;
n0i[14][7] = 0;
n0i[14][1] = 65.93909154;
n0i[14][2] = -74131.45483;
1388 n0i[15][3] = 2.47906;
n0i[15][4] = 0.95806;
n0i[15][5] = 0.45444;
n0i[15][6] = 1.56039;
n0i[15][7] = -1.3756;
n0i[15][1] = 13.07520288;
n0i[15][2] = -5836.943696;
1389 n0i[16][3] = 3.50146;
n0i[16][4] = 1.07558;
n0i[16][5] = 1.01334;
n0i[16][6] = 0;
n0i[16][7] = 0;
n0i[16][1] = 16.8017173;
n0i[16][2] = -2318.32269;
1390 n0i[17][3] = 3.50055;
n0i[17][4] = 1.02865;
n0i[17][5] = 0.00493;
n0i[17][6] = 0;
n0i[17][7] = 0;
n0i[17][1] = 17.45786899;
n0i[17][2] = -2635.244116;
1391 n0i[18][3] = 4.00392;
n0i[18][4] = 0.01059;
n0i[18][5] = 0.98763;
n0i[18][6] = 3.06904;
n0i[18][7] = 0;
n0i[18][1] = 21.57882705;
n0i[18][2] = -7766.733078;
1392 n0i[19][3] = 4;
n0i[19][4] = 3.11942;
n0i[19][5] = 1.00243;
n0i[19][6] = 0;
n0i[19][7] = 0;
n0i[19][1] = 21.5830944;
n0i[19][2] = -6069.035869;
1393 n0i[20][3] = 2.5;
n0i[20][4] = 0;
n0i[20][5] = 0;
n0i[20][6] = 0;
n0i[20][7] = 0;
n0i[20][1] = 10.04639507;
n0i[20][2] = -745.375;
1394 n0i[21][3] = 2.5;
n0i[21][4] = 0;
n0i[21][5] = 0;
n0i[21][6] = 0;
n0i[21][7] = 0;
n0i[21][1] = 10.04639507;
n0i[21][2] = -745.375;
1395 th0i[1][4] = 820.659;
th0i[1][5] = 178.41;
th0i[1][6] = 1062.82;
th0i[1][7] = 1090.53;
1396 th0i[2][4] = 662.738;
th0i[2][5] = 680.562;
th0i[2][6] = 1740.06;
th0i[2][7] = 0;
1397 th0i[3][4] = 919.306;
th0i[3][5] = 865.07;
th0i[3][6] = 483.553;
th0i[3][7] = 341.109;
1398 th0i[4][4] = 559.314;
th0i[4][5] = 223.284;
th0i[4][6] = 1031.38;
th0i[4][7] = 1071.29;
1399 th0i[5][4] = 479.856;
th0i[5][5] = 200.893;
th0i[5][6] = 955.312;
th0i[5][7] = 1027.29;
1400 th0i[6][4] = 438.27;
th0i[6][5] = 198.018;
th0i[6][6] = 1905.02;
th0i[6][7] = 893.765;
1401 th0i[7][4] = 468.27;
th0i[7][5] = 183.636;
th0i[7][6] = 1914.1;
th0i[7][7] = 903.185;
1402 th0i[8][4] = 292.503;
th0i[8][5] = 910.237;
th0i[8][6] = 1919.37;
th0i[8][7] = 0;
1403 th0i[9][4] = 178.67;
th0i[9][5] = 840.538;
th0i[9][6] = 1774.25;
th0i[9][7] = 0;
1404 th0i[10][4] = 182.326;
th0i[10][5] = 859.207;
th0i[10][6] = 1826.59;
th0i[10][7] = 0;
1405 th0i[11][4] = 169.789;
th0i[11][5] = 836.195;
th0i[11][6] = 1760.46;
th0i[11][7] = 0;
1406 th0i[12][4] = 158.922;
th0i[12][5] = 815.064;
th0i[12][6] = 1693.07;
th0i[12][7] = 0;
1407 th0i[13][4] = 156.854;
th0i[13][5] = 814.882;
th0i[13][6] = 1693.79;
th0i[13][7] = 0;
1408 th0i[14][4] = 164.947;
th0i[14][5] = 836.264;
th0i[14][6] = 1750.24;
th0i[14][7] = 0;
1409 th0i[15][4] = 228.734;
th0i[15][5] = 326.843;
th0i[15][6] = 1651.71;
th0i[15][7] = 1671.69;
1410 th0i[16][4] = 2235.71;
th0i[16][5] = 1116.69;
th0i[16][6] = 0;
th0i[16][7] = 0;
1411 th0i[17][4] = 1550.45;
th0i[17][5] = 704.525;
th0i[17][6] = 0;
th0i[17][7] = 0;
1412 th0i[18][4] = 268.795;
th0i[18][5] = 1141.41;
th0i[18][6] = 2507.37;
th0i[18][7] = 0;
1413 th0i[19][4] = 1833.63;
th0i[19][5] = 847.181;
th0i[19][6] = 0;
th0i[19][7] = 0;
1416
1417
1418 bvij[1][2] = 0.998721377;
gvij[1][2] = 1.013950311;
btij[1][2] = 0.99809883;
gtij[1][2] = 0.979273013;
1419 bvij[1][3] = 0.999518072;
gvij[1][3] = 1.002806594;
btij[1][3] = 1.02262449;
gtij[1][3] = 0.975665369;
1420 bvij[1][4] = 0.997547866;
gvij[1][4] = 1.006617867;
btij[1][4] = 0.996336508;
gtij[1][4] = 1.049707697;
1421 bvij[1][5] = 1.00482707;
gvij[1][5] = 1.038470657;
btij[1][5] = 0.989680305;
gtij[1][5] = 1.098655531;
1422 bvij[1][6] = 1.011240388;
gvij[1][6] = 1.054319053;
btij[1][6] = 0.980315756;
gtij[1][6] = 1.161117729;
1423 bvij[1][7] = 0.979105972;
gvij[1][7] = 1.045375122;
btij[1][7] = 0.99417491;
gtij[1][7] = 1.171607691;
1424 bvij[1][8] = 1;
gvij[1][8] = 1.343685343;
btij[1][8] = 1;
gtij[1][8] = 1.188899743;
1425 bvij[1][9] = 0.94833012;
gvij[1][9] = 1.124508039;
btij[1][9] = 0.992127525;
gtij[1][9] = 1.249173968;
1426 bvij[1][10] = 0.958015294;
gvij[1][10] = 1.052643846;
btij[1][10] = 0.981844797;
gtij[1][10] = 1.330570181;
1427 bvij[1][11] = 0.962050831;
gvij[1][11] = 1.156655935;
btij[1][11] = 0.977431529;
gtij[1][11] = 1.379850328;
1428 bvij[1][12] = 0.994740603;
gvij[1][12] = 1.116549372;
btij[1][12] = 0.957473785;
gtij[1][12] = 1.449245409;
1429 bvij[1][13] = 1.002852287;
gvij[1][13] = 1.141895355;
btij[1][13] = 0.947716769;
gtij[1][13] = 1.528532478;
1430 bvij[1][14] = 1.033086292;
gvij[1][14] = 1.146089637;
btij[1][14] = 0.937777823;
gtij[1][14] = 1.568231489;
1431 bvij[1][15] = 1;
gvij[1][15] = 1.018702573;
btij[1][15] = 1;
gtij[1][15] = 1.352643115;
1433 bvij[1][17] = 0.997340772;
gvij[1][17] = 1.006102927;
btij[1][17] = 0.987411732;
gtij[1][17] = 0.987473033;
1434 bvij[1][18] = 1.012783169;
gvij[1][18] = 1.585018334;
btij[1][18] = 1.063333913;
gtij[1][18] = 0.775810513;
1435 bvij[1][19] = 1.012599087;
gvij[1][19] = 1.040161207;
btij[1][19] = 1.011090031;
gtij[1][19] = 0.961155729;
1436 bvij[1][20] = 1;
gvij[1][20] = 0.881405683;
btij[1][20] = 1;
gtij[1][20] = 3.159776855;
1437 bvij[1][21] = 1.034630259;
gvij[1][21] = 1.014678542;
btij[1][21] = 0.990954281;
gtij[1][21] = 0.989843388;
1438 bvij[2][3] = 0.977794634;
gvij[2][3] = 1.047578256;
btij[2][3] = 1.005894529;
gtij[2][3] = 1.107654104;
1439 bvij[2][4] = 0.978880168;
gvij[2][4] = 1.042352891;
btij[2][4] = 1.007671428;
gtij[2][4] = 1.098650964;
1440 bvij[2][5] = 0.974424681;
gvij[2][5] = 1.081025408;
btij[2][5] = 1.002677329;
gtij[2][5] = 1.201264026;
1441 bvij[2][6] = 0.98641583;
gvij[2][6] = 1.100576129;
btij[2][6] = 0.99286813;
gtij[2][6] = 1.284462634;
1442 bvij[2][7] = 0.99608261;
gvij[2][7] = 1.146949309;
btij[2][7] = 0.994515234;
gtij[2][7] = 1.304886838;
1443 bvij[2][8] = 1;
gvij[2][8] = 1.154135439;
btij[2][8] = 1;
gtij[2][8] = 1.38177077;
1444 bvij[2][9] = 1;
gvij[2][9] = 1.078877166;
btij[2][9] = 1;
gtij[2][9] = 1.419029041;
1445 bvij[2][10] = 1;
gvij[2][10] = 1.195952177;
btij[2][10] = 1;
gtij[2][10] = 1.472607971;
1446 bvij[2][11] = 1;
gvij[2][11] = 1.40455409;
btij[2][11] = 1;
gtij[2][11] = 1.520975334;
1447 bvij[2][12] = 1;
gvij[2][12] = 1.186067025;
btij[2][12] = 1;
gtij[2][12] = 1.733280051;
1448 bvij[2][13] = 1;
gvij[2][13] = 1.100405929;
btij[2][13] = 0.95637945;
gtij[2][13] = 1.749119996;
1449 bvij[2][14] = 1;
gvij[2][14] = 1;
btij[2][14] = 0.957934447;
gtij[2][14] = 1.822157123;
1450 bvij[2][15] = 0.972532065;
gvij[2][15] = 0.970115357;
btij[2][15] = 0.946134337;
gtij[2][15] = 1.175696583;
1451 bvij[2][16] = 0.99952177;
gvij[2][16] = 0.997082328;
btij[2][16] = 0.997190589;
gtij[2][16] = 0.995157044;
1452 bvij[2][17] = 1;
gvij[2][17] = 1.008690943;
btij[2][17] = 1;
gtij[2][17] = 0.993425388;
1453 bvij[2][18] = 1;
gvij[2][18] = 1.094749685;
btij[2][18] = 1;
gtij[2][18] = 0.968808467;
1454 bvij[2][19] = 0.910394249;
gvij[2][19] = 1.256844157;
btij[2][19] = 1.004692366;
gtij[2][19] = 0.9601742;
1455 bvij[2][20] = 0.969501055;
gvij[2][20] = 0.932629867;
btij[2][20] = 0.692868765;
gtij[2][20] = 1.47183158;
1456 bvij[2][21] = 1.004166412;
gvij[2][21] = 1.002212182;
btij[2][21] = 0.999069843;
gtij[2][21] = 0.990034831;
1457 bvij[3][4] = 1.002525718;
gvij[3][4] = 1.032876701;
btij[3][4] = 1.013871147;
gtij[3][4] = 0.90094953;
1458 bvij[3][5] = 0.996898004;
gvij[3][5] = 1.047596298;
btij[3][5] = 1.033620538;
gtij[3][5] = 0.908772477;
1459 bvij[3][6] = 1.076551882;
gvij[3][6] = 1.081909003;
btij[3][6] = 1.023339824;
gtij[3][6] = 0.929982936;
1460 bvij[3][7] = 1.174760923;
gvij[3][7] = 1.222437324;
btij[3][7] = 1.018171004;
gtij[3][7] = 0.911498231;
1461 bvij[3][8] = 1.060793104;
gvij[3][8] = 1.116793198;
btij[3][8] = 1.019180957;
gtij[3][8] = 0.961218039;
1462 bvij[3][9] = 1.024311498;
gvij[3][9] = 1.068406078;
btij[3][9] = 1.027000795;
gtij[3][9] = 0.979217302;
1463 bvij[3][10] = 1;
gvij[3][10] = 0.851343711;
btij[3][10] = 1;
gtij[3][10] = 1.038675574;
1464 bvij[3][11] = 1.205469976;
gvij[3][11] = 1.164585914;
btij[3][11] = 1.011806317;
gtij[3][11] = 1.046169823;
1465 bvij[3][12] = 1.026169373;
gvij[3][12] = 1.104043935;
btij[3][12] = 1.02969078;
gtij[3][12] = 1.074455386;
1466 bvij[3][13] = 1;
gvij[3][13] = 0.973386152;
btij[3][13] = 1.00768862;
gtij[3][13] = 1.140671202;
1467 bvij[3][14] = 1.000151132;
gvij[3][14] = 1.183394668;
btij[3][14] = 1.02002879;
gtij[3][14] = 1.145512213;
1468 bvij[3][15] = 0.904142159;
gvij[3][15] = 1.15279255;
btij[3][15] = 0.942320195;
gtij[3][15] = 1.782924792;
1471 bvij[3][18] = 0.949055959;
gvij[3][18] = 1.542328793;
btij[3][18] = 0.997372205;
gtij[3][18] = 0.775453996;
1472 bvij[3][19] = 0.906630564;
gvij[3][19] = 1.024085837;
btij[3][19] = 1.016034583;
gtij[3][19] = 0.92601888;
1473 bvij[3][20] = 0.846647561;
gvij[3][20] = 0.864141549;
btij[3][20] = 0.76837763;
gtij[3][20] = 3.207456948;
1474 bvij[3][21] = 1.008392428;
gvij[3][21] = 1.029205465;
btij[3][21] = 0.996512863;
gtij[3][21] = 1.050971635;
1475 bvij[4][5] = 0.997607277;
gvij[4][5] = 1.00303472;
btij[4][5] = 0.996199694;
gtij[4][5] = 1.01473019;
1476 bvij[4][6] = 1;
gvij[4][6] = 1.006616886;
btij[4][6] = 1;
gtij[4][6] = 1.033283811;
1477 bvij[4][7] = 0.999157205;
gvij[4][7] = 1.006179146;
btij[4][7] = 0.999130554;
gtij[4][7] = 1.034832749;
1478 bvij[4][8] = 1;
gvij[4][8] = 1.045439935;
btij[4][8] = 1;
gtij[4][8] = 1.021150247;
1479 bvij[4][9] = 0.993851009;
gvij[4][9] = 1.026085655;
btij[4][9] = 0.998688946;
gtij[4][9] = 1.066665676;
1480 bvij[4][10] = 1;
gvij[4][10] = 1.169701102;
btij[4][10] = 1;
gtij[4][10] = 1.092177796;
1481 bvij[4][11] = 1;
gvij[4][11] = 1.057666085;
btij[4][11] = 1;
gtij[4][11] = 1.134532014;
1482 bvij[4][12] = 1.007469726;
gvij[4][12] = 1.071917985;
btij[4][12] = 0.984068272;
gtij[4][12] = 1.168636194;
1483 bvij[4][13] = 1;
gvij[4][13] = 1.14353473;
btij[4][13] = 1;
gtij[4][13] = 1.05603303;
1484 bvij[4][14] = 0.995676258;
gvij[4][14] = 1.098361281;
btij[4][14] = 0.970918061;
gtij[4][14] = 1.237191558;
1485 bvij[4][15] = 0.925367171;
gvij[4][15] = 1.10607204;
btij[4][15] = 0.932969831;
gtij[4][15] = 1.902008495;
1487 bvij[4][17] = 1;
gvij[4][17] = 1.201417898;
btij[4][17] = 1;
gtij[4][17] = 1.069224728;
1489 bvij[4][19] = 1.010817909;
gvij[4][19] = 1.030988277;
btij[4][19] = 0.990197354;
gtij[4][19] = 0.90273666;
1492 bvij[5][6] = 0.999243146;
gvij[5][6] = 1.001156119;
btij[5][6] = 0.998012298;
gtij[5][6] = 1.005250774;
1493 bvij[5][7] = 0.999795868;
gvij[5][7] = 1.003264179;
btij[5][7] = 1.000310289;
gtij[5][7] = 1.007392782;
1494 bvij[5][8] = 1.040459289;
gvij[5][8] = 0.999432118;
btij[5][8] = 0.994364425;
gtij[5][8] = 1.0032695;
1495 bvij[5][9] = 1.044919431;
gvij[5][9] = 1.019921513;
btij[5][9] = 0.996484021;
gtij[5][9] = 1.008344412;
1496 bvij[5][10] = 1;
gvij[5][10] = 1.057872566;
btij[5][10] = 1;
gtij[5][10] = 1.025657518;
1497 bvij[5][11] = 1;
gvij[5][11] = 1.079648053;
btij[5][11] = 1;
gtij[5][11] = 1.050044169;
1498 bvij[5][12] = 1;
gvij[5][12] = 1.102764612;
btij[5][12] = 1;
gtij[5][12] = 1.063694129;
1499 bvij[5][13] = 1;
gvij[5][13] = 1.199769134;
btij[5][13] = 1;
gtij[5][13] = 1.109973833;
1500 bvij[5][14] = 0.984104227;
gvij[5][14] = 1.053040574;
btij[5][14] = 0.985331233;
gtij[5][14] = 1.140905252;
1501 bvij[5][15] = 1;
gvij[5][15] = 1.07400611;
btij[5][15] = 1;
gtij[5][15] = 2.308215191;
1503 bvij[5][17] = 1;
gvij[5][17] = 1.108143673;
btij[5][17] = 1;
gtij[5][17] = 1.197564208;
1504 bvij[5][18] = 1;
gvij[5][18] = 1.011759763;
btij[5][18] = 1;
gtij[5][18] = 0.600340961;
1505 bvij[5][19] = 0.936811219;
gvij[5][19] = 1.010593999;
btij[5][19] = 0.992573556;
gtij[5][19] = 0.905829247;
1508
1509
1510 bvij[6][7] = 0.999120311;
gvij[6][7] = 1.00041444;
btij[6][7] = 0.999922459;
gtij[6][7] = 1.001432824;
1511
1512 bvij[6][8] = 1;
gvij[6][8] = 1.002284353;
btij[6][8] = 1;
gtij[6][8] = 1.001835788;
1513 bvij[6][9] = 1;
gvij[6][9] = 1.002779804;
btij[6][9] = 1;
gtij[6][9] = 1.002495889;
1514 bvij[6][10] = 1;
gvij[6][10] = 1.010493989;
btij[6][10] = 1;
gtij[6][10] = 1.006018054;
1515 bvij[6][11] = 1;
gvij[6][11] = 1.021668316;
btij[6][11] = 1;
gtij[6][11] = 1.00988576;
1516 bvij[6][12] = 1;
gvij[6][12] = 1.032807063;
btij[6][12] = 1;
gtij[6][12] = 1.013945424;
1517 bvij[6][13] = 1;
gvij[6][13] = 1.047298475;
btij[6][13] = 1;
gtij[6][13] = 1.017817492;
1518 bvij[6][14] = 1;
gvij[6][14] = 1.060243344;
btij[6][14] = 1;
gtij[6][14] = 1.021624748;
1519 bvij[6][15] = 1;
gvij[6][15] = 1.147595688;
btij[6][15] = 1;
gtij[6][15] = 1.895305393;
1521 bvij[6][17] = 1;
gvij[6][17] = 1.087272232;
btij[6][17] = 1;
gtij[6][17] = 1.161390082;
1523 bvij[6][19] = 1.012994431;
gvij[6][19] = 0.988591117;
btij[6][19] = 0.974550548;
gtij[6][19] = 0.937130844;
1526 bvij[7][8] = 1;
gvij[7][8] = 1.002728434;
btij[7][8] = 1;
gtij[7][8] = 1.000792201;
1527 bvij[7][9] = 1;
gvij[7][9] = 1.01815965;
btij[7][9] = 1;
gtij[7][9] = 1.00214364;
1528 bvij[7][10] = 1;
gvij[7][10] = 1.034995284;
btij[7][10] = 1;
gtij[7][10] = 1.00915706;
1529 bvij[7][11] = 1;
gvij[7][11] = 1.019174227;
btij[7][11] = 1;
gtij[7][11] = 1.021283378;
1530 bvij[7][12] = 1;
gvij[7][12] = 1.046905515;
btij[7][12] = 1;
gtij[7][12] = 1.033180106;
1531 bvij[7][13] = 1;
gvij[7][13] = 1.049219137;
btij[7][13] = 1;
gtij[7][13] = 1.014096448;
1532 bvij[7][14] = 0.976951968;
gvij[7][14] = 1.027845529;
btij[7][14] = 0.993688386;
gtij[7][14] = 1.076466918;
1533 bvij[7][15] = 1;
gvij[7][15] = 1.232939523;
btij[7][15] = 1;
gtij[7][15] = 2.509259945;
1535 bvij[7][17] = 1;
gvij[7][17] = 1.084740904;
btij[7][17] = 1;
gtij[7][17] = 1.173916162;
1536 bvij[7][18] = 1;
gvij[7][18] = 1.223638763;
btij[7][18] = 1;
gtij[7][18] = 0.615512682;
1537 bvij[7][19] = 0.908113163;
gvij[7][19] = 1.033366041;
btij[7][19] = 0.985962886;
gtij[7][19] = 0.926156602;
1539 bvij[7][21] = 1;
gvij[7][21] = 1.214638734;
btij[7][21] = 1;
gtij[7][21] = 1.245039498;
1540 bvij[8][9] = 1;
gvij[8][9] = 1.000024335;
btij[8][9] = 1;
gtij[8][9] = 1.000050537;
1541 bvij[8][10] = 1;
gvij[8][10] = 1.002995876;
btij[8][10] = 1;
gtij[8][10] = 1.001204174;
1542 bvij[8][11] = 1;
gvij[8][11] = 1.009928206;
btij[8][11] = 1;
gtij[8][11] = 1.003194615;
1543 bvij[8][12] = 1;
gvij[8][12] = 1.017880545;
btij[8][12] = 1;
gtij[8][12] = 1.00564748;
1544 bvij[8][13] = 1;
gvij[8][13] = 1.028994325;
btij[8][13] = 1;
gtij[8][13] = 1.008191499;
1545 bvij[8][14] = 1;
gvij[8][14] = 1.039372957;
btij[8][14] = 1;
gtij[8][14] = 1.010825138;
1546 bvij[8][15] = 1;
gvij[8][15] = 1.184340443;
btij[8][15] = 1;
gtij[8][15] = 1.996386669;
1548 bvij[8][17] = 1;
gvij[8][17] = 1.116694577;
btij[8][17] = 1;
gtij[8][17] = 1.199326059;
1550 bvij[8][19] = 1;
gvij[8][19] = 0.835763343;
btij[8][19] = 1;
gtij[8][19] = 0.982651529;
1553 bvij[9][10] = 1;
gvij[9][10] = 1.002480637;
btij[9][10] = 1;
gtij[9][10] = 1.000761237;
1554 bvij[9][11] = 1;
gvij[9][11] = 1.008972412;
btij[9][11] = 1;
gtij[9][11] = 1.002441051;
1555 bvij[9][12] = 1;
gvij[9][12] = 1.069223964;
btij[9][12] = 1;
gtij[9][12] = 1.016422347;
1556 bvij[9][13] = 1;
gvij[9][13] = 1.034910633;
btij[9][13] = 1;
gtij[9][13] = 1.103421755;
1557 bvij[9][14] = 1;
gvij[9][14] = 1.016370338;
btij[9][14] = 1;
gtij[9][14] = 1.049035838;
1558 bvij[9][15] = 1;
gvij[9][15] = 1.188334783;
btij[9][15] = 1;
gtij[9][15] = 2.013859174;
1560 bvij[9][17] = 1;
gvij[9][17] = 1.119954454;
btij[9][17] = 1;
gtij[9][17] = 1.206043295;
1561 bvij[9][18] = 1;
gvij[9][18] = 0.95667731;
btij[9][18] = 1;
gtij[9][18] = 0.447666011;
1562 bvij[9][19] = 0.984613203;
gvij[9][19] = 1.076539234;
btij[9][19] = 0.962006651;
gtij[9][19] = 0.959065662;
1565 bvij[10][11] = 1;
gvij[10][11] = 1.001508227;
btij[10][11] = 1;
gtij[10][11] = 0.999762786;
1566 bvij[10][12] = 1;
gvij[10][12] = 1.006268954;
btij[10][12] = 1;
gtij[10][12] = 1.001633952;
1567 bvij[10][13] = 1;
gvij[10][13] = 1.02076168;
btij[10][13] = 1;
gtij[10][13] = 1.055369591;
1568 bvij[10][14] = 1.001516371;
gvij[10][14] = 1.013511439;
btij[10][14] = 0.99764101;
gtij[10][14] = 1.028939539;
1569 bvij[10][15] = 1;
gvij[10][15] = 1.243461678;
btij[10][15] = 1;
gtij[10][15] = 3.021197546;
1571 bvij[10][17] = 1;
gvij[10][17] = 1.155145836;
btij[10][17] = 1;
gtij[10][17] = 1.233272781;
1572 bvij[10][18] = 1;
gvij[10][18] = 1.170217596;
btij[10][18] = 1;
gtij[10][18] = 0.569681333;
1573 bvij[10][19] = 0.754473958;
gvij[10][19] = 1.339283552;
btij[10][19] = 0.985891113;
gtij[10][19] = 0.956075596;
1576 bvij[11][12] = 1;
gvij[11][12] = 1.006767176;
btij[11][12] = 1;
gtij[11][12] = 0.998793111;
1577 bvij[11][13] = 1;
gvij[11][13] = 1.001370076;
btij[11][13] = 1;
gtij[11][13] = 1.001150096;
1578 bvij[11][14] = 1;
gvij[11][14] = 1.002972346;
btij[11][14] = 1;
gtij[11][14] = 1.002229938;
1579 bvij[11][15] = 1;
gvij[11][15] = 1.159131722;
btij[11][15] = 1;
gtij[11][15] = 3.169143057;
1581 bvij[11][17] = 1;
gvij[11][17] = 1.190354273;
btij[11][17] = 1;
gtij[11][17] = 1.256123503;
1583 bvij[11][19] = 0.828967164;
gvij[11][19] = 1.087956749;
btij[11][19] = 0.988937417;
gtij[11][19] = 1.013453092;
1586 bvij[12][13] = 1;
gvij[12][13] = 1.001357085;
btij[12][13] = 1;
gtij[12][13] = 1.000235044;
1587 bvij[12][14] = 1;
gvij[12][14] = 1.002553544;
btij[12][14] = 1;
gtij[12][14] = 1.007186267;
1588 bvij[12][15] = 1;
gvij[12][15] = 1.305249405;
btij[12][15] = 1;
gtij[12][15] = 2.191555216;
1590 bvij[12][17] = 1;
gvij[12][17] = 1.219206702;
btij[12][17] = 1;
gtij[12][17] = 1.276565536;
1591 bvij[12][18] = 1;
gvij[12][18] = 0.599484191;
btij[12][18] = 1;
gtij[12][18] = 0.662072469;
1595 bvij[13][14] = 1;
gvij[13][14] = 1.00081052;
btij[13][14] = 1;
gtij[13][14] = 1.000182392;
1596 bvij[13][15] = 1;
gvij[13][15] = 1.342647661;
btij[13][15] = 1;
gtij[13][15] = 2.23435404;
1598 bvij[13][17] = 1;
gvij[13][17] = 1.252151449;
btij[13][17] = 1;
gtij[13][17] = 1.294070556;
1600 bvij[13][19] = 1;
gvij[13][19] = 1.082905109;
btij[13][19] = 1;
gtij[13][19] = 1.086557826;
1603 bvij[14][15] = 1.695358382;
gvij[14][15] = 1.120233729;
btij[14][15] = 1.064818089;
gtij[14][15] = 3.786003724;
1605 bvij[14][17] = 1;
gvij[14][17] = 0.87018496;
btij[14][17] = 1.049594632;
gtij[14][17] = 1.803567587;
1606 bvij[14][18] = 1;
gvij[14][18] = 0.551405318;
btij[14][18] = 0.897162268;
gtij[14][18] = 0.740416402;
1607 bvij[14][19] = 0.975187766;
gvij[14][19] = 1.171714677;
btij[14][19] = 0.973091413;
gtij[14][19] = 1.103693489;
1611 bvij[15][17] = 1;
gvij[15][17] = 1.121416201;
btij[15][17] = 1;
gtij[15][17] = 1.377504607;
1617 bvij[16][18] = 1;
gvij[16][18] = 1.143174289;
btij[16][18] = 1;
gtij[16][18] = 0.964767932;
1620 bvij[16][21] = 0.999746847;
gvij[16][21] = 0.993907223;
btij[16][21] = 1.000023103;
gtij[16][21] = 0.990430423;
1622 bvij[17][19] = 0.795660392;
gvij[17][19] = 1.101731308;
btij[17][19] = 1.025536736;
gtij[17][19] = 1.022749748;
1624 bvij[17][21] = 1;
gvij[17][21] = 1.159720623;
btij[17][21] = 1;
gtij[17][21] = 0.954215746;
1625 bvij[18][19] = 1;
gvij[18][19] = 1.014832832;
btij[18][19] = 1;
gtij[18][19] = 0.940587083;
1627 bvij[18][21] = 1;
gvij[18][21] = 1.038993495;
btij[18][21] = 1;
gtij[18][21] = 1.070941866;
1631
1632 for(
int i = 1; i <=
MaxFlds; ++i){
1637 for (
int j = i + 1; j <=
MaxFlds; ++j){
1642 }
1643 }
1644
1645 for (
int i = 1; i <=
MaxMdl; ++i){
1646 for (
int j = 1; j <=
MaxTrmM; ++j){
1650 }
1651 }
1652
1653
1654 T0 = 298.15;
1655 d0 = 101.325 /
RGERG / T0;
1656 for (
int i = 1; i <=
MaxFlds; ++i){
1657 n0i[i][3] =
n0i[i][3] - 1;
1658 n0i[i][2] =
n0i[i][2] + T0;
1659 for (int j = 1; j <= 7; ++j){
1660 n0i[i][j] = Rsr *
n0i[i][j];
1661 }
1662 n0i[i][2] =
n0i[i][2] - T0;
1663 n0i[i][1] =
n0i[i][1] - log(d0);
1664 }
1665 return;
1666
1667
1668
1669
1670
1671
1672
1673
1674
1675
1676
1677
1678
1679
1680
1681
1682
1683
1684
1685
1686}
static double th0i[MaxFlds+1][7+1]
static double xold[MaxFlds+1]
static double n0i[MaxFlds+1][7+1]
static double cijk[MaxMdl+1][MaxTrmM+1]
static double gtij[MaxFlds+1][MaxFlds+1]
static double noik[MaxFlds+1][MaxTrmP+1]
static int kpol[MaxFlds+1]
static double gvij[MaxFlds+1][MaxFlds+1]
static int coik[MaxFlds+1][MaxTrmP+1]
static double Dc[MaxFlds+1]
static int kpolij[MaxMdl+1]
static double btij[MaxFlds+1][MaxFlds+1]
static int kexpij[MaxMdl+1]
static int doik[MaxFlds+1][MaxTrmP+1]
static int mNumb[MaxFlds+1][MaxFlds+1]
static double nijk[MaxMdl+1][MaxTrmM+1]
static double Tc2[MaxFlds+1]
static double toik[MaxFlds+1][MaxTrmP+1]
static double bvij[MaxFlds+1][MaxFlds+1]
static double eijk[MaxMdl+1][MaxTrmM+1]
static double tijk[MaxMdl+1][MaxTrmM+1]
static int dijk[MaxMdl+1][MaxTrmM+1]
static double fij[MaxFlds+1][MaxFlds+1]
static int kexp[MaxFlds+1]
static double Vc3[MaxFlds+1]
static double Tc[MaxFlds+1]
static double gijk[MaxMdl+1][MaxTrmM+1]